FEA Practice Problems
FEA Practice Problems
Topics: Physical space and isoparametric elastic finite element formulation in 1D and 2D, trian-
gle and quadrilaterial finite elements, numerical integration, convergence, recommended practices.
Covers Lectures 10-19, Practicums 3-6
1. The figure below shows one 4-node quadrilateral element in a finite element mesh, with nodal
coordinates given in mm. Assume 4 Gauss quadrature integration points (2 × 2).
(a) In isoparametric space (s, t), what is the location of the Gauss quadrature point closest
to node 2?
For the remaining parts, evaluate the listed quantities at the location s = − √13 , t = √1 .
3
(b) Find the location of the point in physical space (x, y).
(c) Compute [J], |J|, and [J]−1 .
(d) Compute the first row of [B].
(e) Compute the strain component εxx given the nodal displacements: u1 = 0, v1 = 0,
u2 = 0, v2 = −6.4 × 10−3 mm, u3 = 10.8 × 10−3 mm, v3 = −9.2 × 10−3 mm, u4 = 0,
and v4 = −0.8 × 10−3 mm.
(f) If the analytic solution to the displacement field is as follows:
Compute the strain components εxx , εyy , εxy at the location of interest.
1
Space for problem 1 solution
2
2. A three-node triangle element is shown with constant body force acting in the −y direction
of by = −8 × 10−3 N/cm3 distributed uniformly over the whole element as shown (dimensions
in cm). The element is of thickness h = 10 cm.
3
Space for problem 2 solution
4
psi,2 Poisson’s ratio of ν =
3. A thick cantilever plate with Young’s modulus of E = 10 × 106N/cm
0.25, and thickness of h = 10cm (plane strain assumption valid) has a distributed load T~
applied on the top edge, see figure (a). The traction load on the top surface acts only in the
−y direction with a constant load Ty = −50 N/cm
lb/in. Consider the finite element discretization
of the plate using four finite elements as shown in figure (b).
(a) Using 1-point Gauss quadrature evaluate the component of the global stiffness matrix
associated with u5 , the x-displacement of node 5 shown. The [B] matrix for any element
in the mesh evaluated at its center s = 0, t = 0 is
1 1 1
0 − 14
−4 0 4 0 4 0
−1
[B] = 0 − 14 0 − 41 0 1
0 1
cm
4 4
− 14 − 14 − 14 1
4
1
4
1
4
1
4 − 14
(b) How many Gauss quadrature points would be required to exactly compute the corre-
sponding force vector from the traction vector T~ on the top surface?
(c) Compute the component of the global force vector due to the traction at node 2.
(d) You would like to improve the solution convergence by using more elements to represent
the cantilever plate, what type of refinement would that be?
5
Space for problem 3 solution