10-Linear Strain Triangle and Other Types of 2d Elements

Linear Strain Triangle
and other types of 2D

elements
By
S. Ziaei Rad
Linear Strain Triangle (LST or T6)

This element is also called quadratic triangular element.
Quadratic Triangular Element

There are six nodes on this element: three corner nodes and
three midside nodes. Each node has two degrees of freedom
(DOF) as before. The displacements (u, v) are assumed to be
quadratic functions of (x, y),
where bi (i = 1, 2, ..., 12) are constants. From these, the

strains are found to be,

which are linear functions. Thus, we have the linear strain
triangle(LST), which provides better results than the CST.
In the natural coordinate system we defined earlier, the six
shape functions for the LST element are,
in which
. Each of these six shape functions
represents a quadratic form on the element as shown in the
figure.
Shape Function N1 for LST

Displacements can be written as,
The element stiffness matrix is still given by

but here
is quadratic in x and y. In general, the integral
has to be computed numerically.
Linear Quadrilateral Element (Q4)
Linear Quadrilateral Element
Linear Quadrilateral Element (Q4)

There are four nodes at the corners of the quadrilateral
shape. In the natural coordinate system
, the
four shape functions are,
Note that
at any point inside the element, as
expected. The displacement field is given by
which are bilinear functions over the element.
Isoparametric Element
If we use the same parameters (shape functions) to express
Geometry, we are using an isoparametric formulation.
u N1
=
v 0
0
N1
N2
0
x = N1 x1 + N 2 x2 + N 3 x3 + N 4 x4
y = N1 y1 + N 2 y2 + N 3 y3 + N 4 y4
0
N2
N3
0
0
N3
x N1
=
y 0
N4
0
0
N1
u1
v
1
u 2

0 v2

N 4 u3
v3

u 4
v
4
N2
0
0
N2
N3
0
0
N3
N4
0
x1
y
1
x2

0 y2

N 4 x3
y3

x4
y
4
f x y f

x or
= x y f

f
y
x y
J11 J12
[J ] =
=
x
y
J 21 J 22

1
{ } =
=
det J
u v
+
y x
J 22
0
J 21
f
f

x

= [J ] f
f

y

A1
J12
0
0
J 21
J11
J 22
f
f

x
1
f = [J ]

f
y

u

u
0

J11
v
J12

v

A2
u1
u
v

1

0
(1 )
0
(1 + )
0
0 u 2
(1 + )
(1 )
u
v
1 (1 )
0
(
1
)
0
(
1
)
0
(
1
)
0
+
+

2
=
0
(1 )
0
(1 + )
0
(1 )
(1 + ) u3
v 4 0

0
(
1
)
0
(
1
)
0
(
1
)
0
(
1
+
+
v3
u 4

v

4
{ } = [ A1][ A2]{d } = [ B]{d }

1 1
[ K ]e = t [ B]T [ E ][ B] det J dd
1 1
Quadratic Quadrilateral Element

(Q8)
This is the most widely used element for 2-D problems due
to its high accuracy in analysis and flexibility in modeling.

(Q8)
N i = Fi ( , )Gi ( , )
Fi ( , ) Give a value of zero along the sides of the element
That the given node does not contact
Gi ( , ) Select such that when multiply by Fi, it will produce

A value of unity at node i and a value of zero at
other neighboring nodes.
Example: Consider N3
F3 ( , ) = (1 + )(1 + )
G3 ( , ) = c1 + c2 + c3
G3 (1,0) = 0; G3 (0,1) = 0; N 3 (1,1) = F3 (1,1)G3 (1,1) = 1

(Q8)
G3 (1,0) = 0 c1 + c2 = 0
G3 (0,1) = 0 c1 + c3 = 0
N 3 (1,1) = 1 4(c1 + c2 + c3 ) = 1
c1 = 1 / 4
c2 = 1 / 4
c3 = 1 / 4
1
N 3 = (1 + )(1 + )( + 1)
4

(Q8)
There are eight nodes for this element, four corners nodes
and four midside nodes. In the natural coordinate system
the eight shape functions are,
Quadratic Quadrilateral Element (Q8)

Again, we have
at any point inside the element.
The displacement field is given by
which are quadratic functions over the element. Strains

and stresses over a quadratic quadrilateral element are
linear functions, which are better representations.
Notes:
Q4 and T3 are usually used together in a mesh with linear
elements.
Q8 and T6 are usually applied in a mesh composed of
quadratic elements.
Quadratic elements are preferred for stress analysis,
because of their high accuracy and the flexibility in
modeling complex geometry, such as curved boundaries.
Example 3.2
A square plate with a hole at the center and under pressure
in one direction.
The dimension of the plate is 10 in. x 10 in., thickness is

0.1 in. and radius of the hole is 1 in. Assume E = 10x10E6 psi,
v= 0.3 and p = 100 psi. Find the maximum stress in the plate.
Example 3.2
FE Analysis:
From the knowledge of stress concentrations, we should
expect the maximum stresses occur at points A and B on the
edge of the hole. Value of this stress should be around 3p (=
300 psi) which is the exact solution for an infinitely large plate
with a hole.
We use the ANSYS FEA software to do the modeling
(meshing) and analysis, using quadratic triangular (T6 or LST)
linear quadrilateral (Q4) and quadratic quadrilateral (Q8)
elements. Linear triangles (CST or T3) is NOT available in
ANSYS.
The stress calculations are listed in the following table,
along with the number of elements and DOF used, for
comparison.
Example 3.2
Discussions:
Check the deformed shape of the plate
Check convergence (use a finer mesh, if possible)
Less elements (~ 100) should be enough to achieve
the same accuracy with a better or smarter mesh
Well redo this example in next chapter employing the
symmetry conditions.
Example 3.2
FEA Mesh (Q8, 493 elements)

FEA Stress Plot (Q8, 493 elements)
Transformation of Loads
Concentrated load (point forces), surface traction (pressure
loads) and body force (weight) are the main types of loads
applied to a structure. Both traction and body forces need to be
converted to nodal forces in the FEA, since they cannot be
applied to the FE model directly. The conversions of these
loads are based on the same idea (the equivalent-work concept)
which we have used for the cases of bar and beam elements.
Traction on a Q4 element
Suppose, for example, we have a linearly varying traction q
on a Q4 element edge, as shown in the figure. The traction is
normal to the boundary. Using the local (tangential) coordinate
s, we can write the work done by the traction q as,
where t is the thickness, L the side length and

of displacement normal to the edge AB.
the component
For the Q4 element (linear displacement field), we have
The traction q(s), which is also linear, is given in a similar way,
Thus, we have,
and the equivalent nodal force vector is,
Note, for constant q, we have,
For quadratic elements (either triangular or quadrilateral),

the traction is converted to forces at three nodes along the edge,
instead of two nodes.
Traction tangent to the boundary, as well as body forces,
are converted to nodal forces in a similar way.
s = 0 @ = 1
L
ds = d
2
L
s = (1 + )
2
s = L @ =1
q1
q2
L
s
1
1
v = (1 )v1 + (1 + )v2
2
2
1
1
q = (1 )q1 + (1 + )q2
2
2
1
(
1
)
1
2
1
Wq = t [v1 v2 ]
(1 )
1
(1 + ) 2
1
2
q1
q1
L
s
(
1
)
1
2
q1 L
Wq = t [v1 v2 ]
d
1
(1 + ) q1 2
1
2
f1 tL 2 1 q1
=
q
f
1
2
6
2
2
1
q1 L
(1 + ) d
2
q2 2
f1 tLq1 1
=
1
f
2
1- Point Load considered in a usual manner by having a
Structural node at the point.
2- Traction Force As it was seen in the previous example
First the force and the deflection along the side express by
Use of shape functions and then numerical integration will
Be used to calculate the equivalent nodal forces.
3- Body Force A Body force which is a distributed force
Per unit volume, contribute to the global force vector F.
Assume {F } = [ f x f y ]T as constant within each element.
T
T e
u
FdV
=
d
f Where the 8*1 element body force is given by
1 1
f x
e
T
f = te N det Jdd
f y
1 1
Stress Calculation
The stress in an element is determined by the following
relation,
where B is the strain-nodal displacement matrix and d is the

nodal displacement vector which is known for each element
once the global FE equation has been solved.
Stresses can be evaluated at any point inside the element
(such as the center) or at the nodes. Contour plots are usually
used in FEA software packages (during post-process) for users
to visually inspect the stress results.
The von Mises Stress

The von Mises stress is the effective or equivalent stress for
2-D and 3-D stress analysis. For a ductile material, the stress
level is considered to be safe, if
where
is the von Mises stress and
the yield stress of the
material. This is a generalization of the 1-D (experimental)
result to 2-D and 3-D situations.
The von Mises stress is defined by
in which
and
are the three principle stresses at the
considered point in a structure.
The von Mises Stress

For 2-D problems, the two principle stresses in the plane
are determined by
Thus, we can also express the von Mises stress in terms of

the stress components in the xy coordinate system. For plane
stress conditions, we have,
Averaged Stresses:
Stresses are usually averaged at nodes in FEA software
packages to provide more accurate stress values. This option
should be turned off at nodes between two materials or other
geometry discontinuity locations where stress discontinuity does
exist.
Discussions
1) Know the behaviors of each type of elements:
T3 and Q4: linear displacement, constant strain and stress;
T6 and Q8: quadratic displacement, linear strain and stress.
2) Choose the right type of elements for a given problem:
When in doubt, use higher order elements or a finer mesh.
3) Avoid elements with large aspect ratios and corner angles:
where Lmax and Lmin are the largest and smallest characteristic
lengths of an element, respectively.
Elements with Bad Shapes Elements with Nice Shapes
Discussions
4) Connect the elements properly:
Dont leave unintended gaps or free elements in FE models.
Improper connections (gaps along AB and CD)

10-Linear Strain Triangle and Other Types of 2d Elements

Uploaded by

Copyright:

Available Formats

10-Linear Strain Triangle and Other Types of 2d Elements

Uploaded by

Document Information

Original Description:

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

10-Linear Strain Triangle and Other Types of 2d Elements

Uploaded by

Copyright:

Available Formats

What are the different types of elements discussed in the document?

What are the different types of elements discussed in the document?

What are some guidelines for choosing appropriate elements for analysis?

What are some guidelines for choosing appropriate elements for analysis?

Linear Strain Triangle

and other types of 2D

Linear Strain Triangle (LST or T6)

Quadratic Triangular Element

Linear Strain Triangle

where bi (i = 1, 2, ..., 12) are constants. From these, the

Linear Strain Triangle

Linear Strain Triangle

Shape Function N1 for LST

The element stiffness matrix is still given by

Linear Quadrilateral Element (Q4)

Linear Quadrilateral Element

Linear Quadrilateral Element (Q4)

which are bilinear functions over the element.

{ } = [ A1][ A2]{d } = [ B]{d }

Quadratic Quadrilateral Element

Quadratic Quadrilateral Element

Quadratic Quadrilateral Element

Gi ( , ) Select such that when multiply by Fi, it will produce

Quadratic Quadrilateral Element

Quadratic Quadrilateral Element

Quadratic Quadrilateral Element (Q8)

at any point inside the element.

The displacement field is given by

which are quadratic functions over the element. Strains

The dimension of the plate is 10 in. x 10 in., thickness is

FEA Mesh (Q8, 493 elements)

where t is the thickness, L the side length and

For the Q4 element (linear displacement field), we have

The traction q(s), which is also linear, is given in a similar way,

and the equivalent nodal force vector is,

Note, for constant q, we have,

For quadratic elements (either triangular or quadrilateral),

where B is the strain-nodal displacement matrix and d is the

The von Mises Stress

The von Mises Stress

Thus, we can also express the von Mises stress in terms of

Elements with Bad Shapes Elements with Nice Shapes

Improper connections (gaps along AB and CD)

You might also like