Module 5

Medium Answer Questions

2. List and describe steps involved in FEA.
Steps Involved in FEA
1. Problem Definition and Idealization
The physical problem is identified and simplified into a mathematical model.
Activities
• Define geometry
• Determine loading conditions
• Specify material properties
• Define assumptions
Example
A bridge beam may be idealized as a 2D or 3D structural model.
2. Discretization (Meshing)
The structure is divided into smaller finite elements connected at nodes.
Types of Elements
• 1D elements (bars, beams)
• 2D elements (triangles, quadrilaterals)
• 3D elements (tetrahedral, hexahedral)
Purpose
• Converts complex geometry into manageable parts.
Example
A plate divided into triangular mesh elements.
3. Selection of Element Type and Shape Functions
Appropriate finite elements and interpolation functions are selected.
Shape Functions
Used to approximate variation of:
• Displacement
• Temperature
• Pressure
within an element.
Example
• Linear triangular element
• Quadratic quadrilateral element
4. Formulation of Element Equations
Mathematical equations are developed for each element.
Includes
• Stiffness matrix
• Load vector
• Material behavior equations
General form:
[𝐾]𝑒 {𝑞}𝑒 = {𝐹}𝑒
where:
• [𝐾]𝑒 = element stiffness matrix
• {𝑞}𝑒 = nodal displacement vector
• {𝐹}𝑒 = force vector
5. Assembly of Global Equations
Description
All element equations are combined to form the global system equation.
Global Equation
 [𝐾]{𝑄} = {𝐹}
where:
• [𝐾]= global stiffness matrix
• {𝑄}= global nodal displacement vector
• {𝐹}= global force vector
6. Application of Boundary Conditions
Known constraints and loads are applied.
Types
• Displacement boundary conditions
• Force boundary conditions
Example
• Fixed support
• Roller support
• Applied pressure load
7. Solution of Equations
The global equations are solved to obtain unknown nodal values.
Methods
• Direct methods
• Iterative methods
Outputs
• Displacements
• Temperatures
• Pressures
8. Post-Processing
Results are interpreted and displayed graphically.
Calculated Results
• Stress
• Strain
• Deformation
• Heat flux
• Vibration modes
Tools
• Contour plots
• Vector diagrams
• Animations
9. Verification and Validation
The obtained results are checked for accuracy.
Methods
• Compare with theoretical solutions
• Experimental validation
• Mesh refinement study

3. State different types of finite elements and their applications.

1. One-Dimensional (1D) Elements
These elements have one significant dimension (length).
(a) Bar Element
Used for axial load carrying members.
Characteristics
• Two-node element
• Only axial deformation considered
Applications
• Trusses
• Tie rods
 • Structural members
Example
Bridge truss analysis.
(b) Beam Element
Used for bending analysis.
Characteristics
• Resists bending and shear
• Includes rotational degrees of freedom
Applications
• Building beams
• Machine frames
• Shafts
Example
Cantilever beam analysis.
(c) Frame Element
Combination of bar and beam behavior.
Characteristics
• Carries axial, bending, and torsional loads
Applications
• Building frames
• Aircraft structures
• Machine supports
2. Two-Dimensional (2D) Elements
Used when thickness is small compared to other dimensions.
(a) Triangular Element
Three-node or six-node element.
Characteristics
• Suitable for irregular geometries
Applications
• Stress analysis of plates
• Complex-shaped domains
Example
Constant Strain Triangle (CST).
(b) Quadrilateral Element
Four-node or higher-order element.
Characteristics
• Better accuracy than triangular elements
Applications
• Plate and shell analysis
• Structural mechanics
Example
Plane stress analysis of plates.
(c) Plane Stress Element
Used for thin structures loaded in-plane.
Applications
• Thin plates
• Sheet metal components
Example
Automobile body panels.
(d) Plane Strain Element
Used when strain in one direction is negligible.
Applications
 • Dams
• Long tunnels
• Retaining walls
(e) Axisymmetric Element
Used for rotationally symmetric structures.
Applications
• Pressure vessels
• Cylinders
• Pipes
Example
Rocket nozzles.
3. Three-Dimensional (3D) Elements
Used for fully three-dimensional stress analysis.
(a) Tetrahedral Element
Four-node or ten-node solid element.
Characteristics
• Suitable for complex geometries
Applications
• Biomedical models
• Automotive parts
Example
Engine block analysis.
(b) Hexahedral (Brick) Element
Eight-node solid element.
Characteristics
• High accuracy
• Efficient numerical performance
Applications
• Structural analysis
• Solid mechanics
Example
Machine component stress analysis.
(c) Wedge and Pyramid Elements
Transition elements used between meshes.
Applications
• Complex mesh generation
• Hybrid meshing
4. Shell Elements
Used for thin-walled curved structures.
Characteristics
• Combines membrane and bending behavior
Applications
• Aircraft fuselage
• Storage tanks
• Car bodies
5. Plate Elements
Used for flat structures subjected to transverse loads.
Applications
• Floor slabs
• Ship decks
• Machine covers
6. Spring Elements
Represents elastic connections.
Applications
• Suspension systems
• Vibration isolation
7. Mass Elements
Represents concentrated masses.
Applications
• Dynamic and vibration analysis
8. Thermal Elements
Used for heat transfer analysis.
Applications
• Heat exchangers
• Cooling systems
• Furnaces
9. Fluid Elements
Used in fluid flow and CFD analysis.
Applications
• Aerodynamics
• Pipe flow
• Turbomachinery
10. Contact Elements
Used to model interaction between surfaces.
Applications
• Gear contact
• Bearing analysis
• Metal forming
5. State two properties of stiffness matrix.
Two Properties of Stiffness Matrix
1. Symmetric Property
The stiffness matrix is symmetric, which means:
[K]=[K]T
or
𝐾𝑖𝑗 = 𝐾𝑗𝑖
Reason
This property arises from:
• Maxwell’s reciprocal theorem
• Conservation of energy principles
Example
𝑘 𝑘12
[𝐾] = [ 11 ]
𝑘12 𝑘22
Here,
𝑘12 = 𝑘21
2. Positive Definite Property
The stiffness matrix is positive definite for a stable structure.
For any nonzero displacement vector {𝑞}:
{𝑞}𝑇 [𝐾]{𝑞} > 0
Meaning
• Strain energy stored in the structure is always positive.
• Ensures structural stability.
Importance
• Guarantees unique solutions in FEA.
• Prevents rigid body instability.
6. What is the role of transformation matrix while analysing truss elements?
Role of Transformation Matrix
1. Converts Local Displacements to Global Displacements
The transformation matrix relates local and global displacement vectors.
{𝑞}𝑙𝑜𝑐𝑎𝑙 = [𝑇]{𝑞}𝑔𝑙𝑜𝑏𝑎𝑙
where:
• [𝑇]= transformation matrix
2. Converts Local Stiffness Matrix into Global Stiffness Matrix
The element stiffness matrix is first derived in local coordinates and then transformed into global coordinates using:
[𝐾]𝑔 = [𝑇]𝑇 [𝐾]𝑙 [𝑇]
where:
• [𝐾]𝑔 = global stiffness matrix
• [𝐾]𝑙 = local stiffness matrix
• [𝑇]= transformation matrix
3. Enables Assembly of Structure Stiffness Matrix
All elements in a truss may have different orientations.
Transformation allows all element matrices to be expressed in a common global coordinate system for proper
assembly.
4. Simplifies Element Formulation
Element equations are easier to derive in local coordinates because:
• Axial deformation occurs only along the member axis.
• The local stiffness matrix becomes simple.
Transformation transfers these simple equations into the global system.
5. Helps in Calculating Global Forces and Reactions
Nodal forces and support reactions are ultimately required in global coordinates.
Transformation ensures accurate force conversion.
7. Define plane stress element and give example.
A plane stress element is a two-dimensional finite element used for analyzing thin structures in which the stress
perpendicular to the plane is assumed to be negligible.
In plane stress condition:
𝜎𝑧 = 0
Example: Thin Steel Plate with a Hole
A thin steel plate subjected to tensile force in its plane can be analyzed using a plane stress element because:
• Thickness is very small
• Loads act only in the plane
• Out-of-plane stresses are negligible
This type of analysis is commonly performed in:
• Aircraft fuselage panels
• Structural plates
• Machine covers

8. Differentiate displacement and structural boundary conditions.

Aspect Displacement Boundary Condition Structural Boundary Condition
Definition Specifies displacement constraints Specifies applied loads or forces
Also Called Essential boundary condition Natural boundary condition
Known Quantity Displacement or rotation Force or stress
Purpose Restricts movement Produces deformation
Typical Examples Fixed support, roller support Point load, pressure load
Effect on Structure Controls motion Creates stress and strain
9. Differentiate shape function and natural coordinate of a finite element.
Difference Between Shape Function and Natural Coordinate in Finite Element Analysis
 Aspect Shape Function Natural Coordinate
A mathematical function used to interpolate field A local coordinate system used to
Definition variables (such as displacement, temperature, etc.) describe the position of points inside a
within an element. finite element.
Simplifies element geometry and
Purpose Approximates variation of variables inside the element.
numerical integration.
Defines location inside the element
Function Relates nodal values to values at any interior point.
independent of global coordinates.
Representat
Usually denoted by (N_1, N_2, N_3,\dots) Usually denoted by (\xi, \eta, \zeta)
ion
Nature Mathematical interpolation equation. Coordinate transformation parameter.
Dependenc Used as independent variables in element
Depends on coordinates and element type.
y formulation.
Role in Used to derive displacement field, strain matrix, and Used for coordinate transformation and
FEA stiffness matrix. Gaussian integration.
Example
for 1D [N_1=\frac{1-\xi}{2}, \quad N_2=\frac{1+\xi}{2}] For a 1D element: (-1 \le \xi \le 1)
Element
Physical Indicates contribution of each node to displacement at a Indicates relative position of a point
Meaning point. inside the element.
Mapping complex elements into standard
Application Interpolation of displacement, temperature, stress, etc.
parent elements.
10. How polynomials suits to formulate shape functions.
1. Easy Mathematical Representation
Polynomials are simple mathematical expressions that can easily represent variation within an element.
Example
For a 1D linear element:
𝑢(𝑥) = 𝑎1 + 𝑎2 𝑥
This polynomial smoothly approximates displacement variation.
2. Continuity and Smoothness
Polynomials provide continuous functions, ensuring smooth variation of displacement and strain between nodes.
This is important for:
• Structural compatibility
• Accurate stress calculation
3. Ability to Approximate Complex Functions
Higher-order polynomials can approximate curved or complex distributions accurately.
Examples
• Linear polynomial → straight-line variation
• Quadratic polynomial → curved variation
• Cubic polynomial → more complex variation
4. Satisfy Nodal Conditions
Polynomial-based shape functions can be constructed to satisfy:
Nodal Property
At node 𝑖:
𝑁𝑖 = 1
and at all other nodes:
𝑁𝑗 = 0
This ensures exact interpolation at nodes.
5. Easy Differentiation and Integration
FEA requires:
• Differentiation for strain calculation
• Integration for stiffness matrix formation
Polynomials are mathematically convenient for both operations.
6. Supports Different Element Orders
Polynomials can generate:
• Linear elements
• Quadratic elements
• Cubic elements
• Higher-order elements
depending on accuracy requirements.
12. Distinguish with example, plane stress and plain strain problems.
Difference Between Plane Stress and Plane Strain
Aspect Plane Stress Plane Strain
Definition Out-of-plane stress is negligible Out-of-plane strain is negligible
Zero Quantity (\sigma_z = 0) (\epsilon_z = 0)
Applicable Structures Thin bodies Very long bodies
Thickness/Length Condition Thickness is very small Length is very large
Out-of-Plane Stress Negligible May exist
Out-of-Plane Strain May exist Zero
Common Applications Thin plates, sheet metal Dams, tunnels
Example Thin steel plate under tension Long concrete dam
14. Differentiate between a Rod element and Bar element.
Difference Between Rod Element and Bar Element
Aspect Rod Element Bar Element
A finite element that resists only axial A structural element that primarily carries axial loads
Definition
tensile or compressive loads. and may be oriented in space.
Degrees of Usually, one degree of freedom per May have one or more DOF depending on
Freedom node (axial displacement only). formulation and application.
Type of Load Axial load only Axial load mainly
Stress Considered Only normal stress due to axial force Mainly axial stress
Bending Effect Neglected Neglected in basic bar element
Geometry Straight slender member Straight structural member
Orientation Usually aligned along one axis Can be inclined in 2D or 3D space
Element
Simpler mathematical model More general structural formulation
Formulation
Applications Springs, tie rods, truss members Truss structures, framework members
Example A straight steel tie rod under tension A truss member in a bridge structure
4. Explain the following terms briefly (i) Discretization process (ii) Global Stiffness Matrix (iii) Quad4
element
(i) Discretization Process
Definition
Discretization is the process of dividing a continuous body or structure into a finite number of smaller parts called
elements, which are connected at points called nodes.
This process converts a complex physical problem into a set of simpler finite element equations.
Purpose
• Simplifies complex geometry
• Enables numerical analysis using FEA
• Improves solution accuracy
Steps in Discretization
1. Divide the structure into finite elements.
2. Select element type.
3. Define nodes and connectivity.
4. Generate mesh.
Example
A rectangular plate may be divided into:
 • Triangular elements
• Quadrilateral elements
for stress analysis.
Importance
• Accuracy depends on mesh quality.
• Finer mesh generally gives more accurate results.
(ii) Global Stiffness Matrix
Definition
The global stiffness matrix is the combined stiffness matrix of the entire structure obtained by assembling the stiffness
matrices of all individual elements.
It represents the overall force-displacement relationship of the complete system.
General Equation
[𝐾]{𝑄} = {𝐹}
where:
• [𝐾]= Global stiffness matrix
• {𝑄}= Global displacement vector
• {𝐹}= Global force vector
Characteristics
• Square matrix
• Symmetric matrix
• Sparse matrix
• Size depends on total degrees of freedom
Formation
1. Derive element stiffness matrices.
2. Transform into global coordinates if necessary.
3. Assemble according to node connectivity.
Example
In truss analysis, stiffness matrices of all truss members are assembled to form the global stiffness matrix of the entire
truss.
Importance
• Used to determine unknown nodal displacements.
• Fundamental equation in FEA.
(iii) Quad4 Element
Definition
A Quad4 element is a four-node quadrilateral finite element used for two-dimensional analysis.
Each node usually has:
• Two displacement degrees of freedom:
o 𝑢in x-direction
o 𝑣in y-direction
Geometry
The element has:
• 4 corner nodes
• Quadrilateral shape
Characteristics
• Two-dimensional element
• Bilinear shape functions
• Suitable for plane stress and plane strain problems
Shape Functions
Shape functions are defined using natural coordinates (𝜉 , 𝜂):
1
𝑁1 = (1 − 𝜉)(1 − 𝜂)
4
 1
𝑁2 = (1 + 𝜉)(1 − 𝜂)
4
1
𝑁3 = (1 + 𝜉)(1 + 𝜂)
4
1
𝑁4 = (1 − 𝜉)(1 + 𝜂)
4
Applications
• Plate analysis
• Stress analysis
• Heat transfer problems
• Structural mechanics
Advantages
• Better accuracy than triangular elements
• Suitable for regular geometries
• Efficient numerical performance
Example
Stress analysis of a rectangular metal plate using quadrilateral mesh elements.