isoparametric elements structural mechanics displacement-based formulations
Post on 19-Dec-2015
237 views
TRANSCRIPT
![Page 1: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/1.jpg)
Isoparametric Elements
Structural MechanicsDisplacement-based Formulations
![Page 2: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/2.jpg)
Fundamental Dilemma
• A primary reason engineers go to FEA is complex geometry
• But elements give the most accurate results when they have regular shapes (isosceles triangles, squares)
• You should always minimize element distortion when you create a mesh (more on this later …)
• It is also important to understand how element shape is managed …interpolation (shape) functions
![Page 3: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/3.jpg)
Reduced Accuracy
• These elements work, but not well …
![Page 4: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/4.jpg)
Isoparametric Elements
• There are two roles of interpolation in FEA:– Defining the location of interior points within an element in terms of
nodal values (geometry interpolation)– Defining the displacement of interior points within an element in
terms of nodal values (result interpolation)
• There is no fundamental reason why both types of interpolation must be conducted in the same way
• But a common class of highly versatile elements does just that – Iso = same; the same basis for geometry and result interpolation
![Page 5: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/5.jpg)
Bilinear Quadrilateral (Q4)
• Interpolation involves the summation of nodal values multiplied by corresponding shapes functions
€
x
y
⎧ ⎨ ⎩
⎫ ⎬ ⎭=
N ix i∑N iy i∑
⎧ ⎨ ⎪
⎩ ⎪
⎫ ⎬ ⎪
⎭ ⎪= N[ ] c{ }
€
u
v
⎧ ⎨ ⎩
⎫ ⎬ ⎭=
N iui∑N iv i∑
⎧ ⎨ ⎪
⎩ ⎪
⎫ ⎬ ⎪
⎭ ⎪= N[ ] d{ }
€
c{ } = x1 y1 x2 y2 x3 y3 x4 y4[ ]T
d{ } = u1 v1 u2 v2 u3 v3 u4 v4[ ]T
N[ ] =N1 0 N2 0 N3 0 N4 0
0 N1 0 N2 0 N3 0 N4
⎡
⎣ ⎢
⎤
⎦ ⎥
- where -geometry interpolation field variable interpolation
nodal coordinates
nodal displacements
shape functions
![Page 6: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/6.jpg)
Shape Functions• Shape functions have a value of 1.0 at the “corresponding” node and a
value of 0.0 at all others (the function “belongs” to a node)• They span a normalized domain, typically [-1,1] over each spatial
dimension
€
N1 =1
41− a( ) 1−b( )
€
N2 =1
41+ a( ) 1−b( )€
N3 =1
41+ a( ) 1+ b( )
€
N4 =1
41− a( ) 1+ b( )
![Page 7: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/7.jpg)
Element Geometry Interpolation• Edges of adjacent elements match (no overlaps, gaps) as long as common
nodes are shared• There are consistent interior point locations defined by the interpolation
functions (e.g. you can define the “center” of an element)
![Page 8: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/8.jpg)
Example 2D Element
N1 = (3,2)N2 = (11,3)N3 = (10,10)N4 = (4,9)
€
xmid = 3( )1
41− 0( ) 1− 0( ) + 11( )
1
41+ 0( ) 1− 0( ) + 10( )
1
41+ 0( ) 1+ 0( ) + 4( )
1
41− 0( ) 1+ 0( ) = 7
€
ymid = 2( )1
41− 0( ) 1− 0( ) + 3( )
1
41+ 0( ) 1− 0( ) + 10( )
1
41+ 0( ) 1+ 0( ) + 9( )
1
41− 0( ) 1+ 0( ) = 6
N1N2
N3N4
X
Y
The shape functions establish a geometric equivalence between elements with different node locations…
![Page 9: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/9.jpg)
XY
Field Quantity Interpolation• We assume the field quantity (e.g. x-component of displacement) varies within the element
as a sum of node-weighted shape functions (visualized here as height above the plane)• If the field quantity actually does vary in this way (or close to it) then the element choice
(size, order) is justified
u1 = 2u2 = 3u3 = 4u4 = 5
€
umid = 2( )1
41− 0( ) 1− 0( ) + 3( )
1
41+ 0( ) 1− 0( ) + 4( )
1
41+ 0( ) 1+ 0( ) + 5( )
1
41− 0( ) 1+ 0( ) = 3.5
nodal values
At the element “center” …
![Page 10: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/10.jpg)
XY
Multiple Elements
N1 = (3,2)N2 = (11,3)N3 = (10,10)N4 = (4,9)
N5 = (0,-5)N6 = (9,-3)N3 = (11,3)N4 = (3,2)
element 1
element 2
……
u1 = 2u2 = 3u3 = 4u4 = 5u5 = 1u6 = 2.5
![Page 11: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/11.jpg)
C0 Continuous
• A key feature of these elements is their continuity• The value (C0) of an interpolated quantity is continuous across
element boundaries– No geometric “gaps” or “overlaps” along the edges of shared elements– Field quantities transition without any “steps” in value
• But the slope (C1) (derivative) of a field quantity is not continuous along element edges– As element density increases the transitions become less abrupt– It is possible to construct C1 shape functions, but they require “slope
nodes” and are not as computationally efficient as mesh refinement
![Page 12: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/12.jpg)
Higher-Order Elements• The overall scheme stays exactly the same• But there are more nodes and different shape function definitions
€
N1 =1
41− a( ) 1−b( ) −
1
2N8 + N5( )
€
N2 =1
41+ a( ) 1−b( ) −
1
2N5 + N6( )
€
N3 =1
41+ a( ) 1+ b( ) −
1
2N6 + N7( )
€
N4 =1
41− a( ) 1+ b( ) −
1
2N7 + N8( )
€
N5 =1
21− a2( ) 1−b( )
€
N6 =1
21+ a( ) 1−b2
( )
€
N7 =1
21− a2( ) 1+ b( )
€
N8 =1
21− a( ) 1−b2
( )
![Page 13: Isoparametric Elements Structural Mechanics Displacement-based Formulations](https://reader033.vdocument.in/reader033/viewer/2022061614/56649d2d5503460f94a03496/html5/thumbnails/13.jpg)
Multiple Elements