ppt

Report
MANE 4240 & CIVL 4240
Introduction to Finite Elements
Prof. Suvranu De
Four-noded
rectangular element
Reading assignment:
Logan 10.2 + Lecture notes
Summary:
• Computation of shape functions for 4-noded quad
• Special case: rectangular element
• Properties of shape functions
• Computation of strain-displacement matrix
• Example problem
•Hint at how to generate shape functions of higher order
(Lagrange) elements
Finite element formulation for 2D:
Step 1: Divide the body into finite elements connected to each
other through special points (“nodes”)
py
v3
3
px
4
3
u3
u 1 
v4
2
v 
v2
Element ‘e’
v
1 

1
4
u 2 
u
u4
ST
 
v1
2 u2
v 2 
y
d 
x
y
u3 

Su
u1
v 3 
1
 
v
x
x
u 4 
u
v 
 4
Summary: For each element
Displacement approximation in terms of shape functions
uNd
Strain approximation in terms of strain-displacement matrix
ε Bd
Stress approximation
  DB d
Element stiffness matrix
k   e B D B dV
T
V
Element nodal load vector
f   e N X dV   e N T S dS
V
ST

 

 

T
f
T
b
f
S
Constant Strain Triangle (CST) : Simplest 2D finite element
v1
1
(x1,y1)
v2
y
u1
v3
(x3,y3)
v
(x,y)
2 (x2,y2)
u
u2
3
u3
u (x,y)  N1(x,y) u1  N2(x,y) u 2  N3(x,y) u 3
v (x,y)  N1(x,y) v1  N2(x,y) v2  N3(x,y) v3
x
• 3 nodes per element
• 2 dofs per node (each node can move in x- and y- directions)
• Hence 6 dofs per element
Formula for the shape functions are
v1
v3
1
u1
(x3,y3)
(x1,y1)
u3
v2 v
u
3
y
(x,y)
where
2 (x2,y2)
u2
a1  b1 x  c1 y
2A
a2  b2 x  c2 y
N2 
2A
a3  b3 x  c3 y
N3 
2A
N1 
x
1 x1
1
A  area of triangle  det 1 x 2
2
1 x 3
y1 
y2 
y3 
a1  x2 y3  x3 y2
b1  y2  y3
c1  x3  x2
a2  x3 y1  x1 y3
a3  x1 y2  x2 y1
b2  y3  y1
b3  y1  y2
c2  x1  x3
c3  x2  x1
Approximation of displacements
uNd
u (x, y)  N1
u

v (x, y)  0
0
N1
N2
0
0
N2
Approximation of the strains
 x

   y

 xy
u


x



 v


 y



 u
v

 y
x

 N1(x, y)
0


x

N1(x, y)
B 0

y
 N (x, y) N (x, y)
1
 1
x
 y
N 2(x, y)
x
0
N 2(x, y)
y




  Bd




0
N 2(x, y)
y
N 2(x, y)
x
N3
0
 u1 
v 
 1
0  u 2 
 

N 3  v 2 
u 3 
 
 v 3 

0 
b1 0 b2 0 b3 0 

N 3(x, y) 1 
0
  0 c1 0 c2 0 c3 
y  2 A
c1 b1 c2 b2 c3 b3 
N 3(x, y) N 3(x, y)

y
x 
N 3(x, y)
x
Element stiffness matrix
t
k   e B D B dV
T
V
Since B is constant
A
k  B D B e dV  B D B At
T
T
V
Element nodal load vector
f   e N X dV   e N T S dS
V
ST

 

 

T
f
T
b
f
S
t=thickness of the element
A=surface area of the element
Class exercise
For the CST shown below, compute the vector of nodal loads due to surface traction
f S   e N T S dS
T
ST
1
f S  t
l13
y
fS2y
fS3y
fS2x
2
(0,0)
fS3x
3
py=-1
(1,0)
x
T
N
e
along 23
T S dS
Class exercise
f S  t
l13
1
y
fS2y
T
2
fS3x
3
py=-1
along 23
T S dS
0
TS   
 1
fS3y
fS2x
N
e
x
The only nonzero nodal loads are
1
f S2 y  t 
x10
f S3 y  t 
x 0
N 2 along 23 p y dx
N3 along 23 p y dx
a2  b2 x  x3 y1  x1 y3    y3  y1 x
 a  b x  c2 y 
N 2 along 23   2 2



2A
2A
2A

 y 0
y1  y1 x
y1 (1  x)
y1 (1  x)



1 x1 y1 
1 x1 y1  y1 ( x3  x2 )
det 1 x 2 y2  det 1 x 2 0 
1 x 3 y3 
1 x 3 0 
 1 x
(can you derive this simpler?)
1
 f S2 y  t 
x 0
N 2 along 23 p y dx
1
 t  (1  x)(1) dx
x 0
t

2
Now compute
1
f S3 y  t 
x 0
N3 along 23 p y dx
4-noded rectangular element with edges parallel to the
coordinate axes:
(x4,y4)
4
3 (x3,y3)
4
v
u (x, y)   Ni (x, y)ui
u
2b
(x,y)
i 1
4
v (x, y)   N i (x, y)vi
i 1
y
1
(x1,y1)
2a
2 (x2,y2)
x
• 4 nodes per element
• 2 dofs per node (each node can move in x- and y- directions)
•8 dofs per element
Generation of N1:
At node 1
y
l1 ( x) 
x  x2
x1  x2
has the property
4
l1 ( x2 )  0
2b
N1
1
2
2a
l1 ( x1 )  1
l1(y)
3
Similarly
1
y  y4
y1  y 4
has the property
l1 ( y1 )  1
x
l1 ( y 4 )  0
l1(x)
1
l1 ( y) 
Hence choose the shape function at node 1 as
 x  x2
N1  l1 ( x)l1 ( y )  
 x1  x 2
 y  y 4

 y1  y 4

1
 
x  x2  y  y 4 
 4ab
Using similar arguments, choose
1
x  x2  y  y 4 
4ab
1
x  x1  y  y3 
N2  
4ab
1
x  x4  y  y 2 
N3 
4ab
1
x  x3  y  y1 
N4  
4ab
N1 
Properties of the shape functions:
1. The shape functions N1, N2 , N3 and N4 are bilinear functions
of x and y
2. Kronecker delta property
 1 at node' i'
Ni ( x, y)  
0 at other nodes
3. Completeness
4
N
i 1
i
1
4
N x
i 1
i
4
N
i 1
i
i
x
yi  y
3. Along lines parallel to the x- or y-axes, the shape functions
are linear. But along any other line they are nonlinear.
4. An element shape function related to a specific nodal point is
zero along element boundaries not containing the nodal point.
5. The displacement field is continuous across elements
6. The strains and stresses are not constant within an element
nor are they continuous across element boundaries.
The strain-displacement relationship
x 
 
   y 
 xy 
 
u 1 
v 
 N 1(x, y)
 1 
N 3 (x, y)
N 2 (x, y)
N 4 (x, y)
0
0
0
0

 u 2 

x

x

x

x

 

N
(x,
y)

N
(x,
y)

N
(x,
y)

N
(x,
y)
3
1
2
4
 v 2 

0
0
0
0
 

y
y
y
y  u 3 
 N (x, y) N (x, y) N (x, y) N (x, y) N (x, y) N (x, y) N (x, y) N (x, y)
3
3
1
2
2
4
4
 1
 v 3 
x
y
x
y
x
y
x  u 
 y
  4 
B
v 
 4
0
y3  y
0
y  y2
0
y1  y
0 
 y  y4
1 

B
0
x

x
0
x

x
0
x

x
0
x

x
2
1
4
3

4ab 
 x  x 2 y  y 4 x1  x y3  y x  x 4 y  y 2 x3  x y1  y 
Notice that the strains (and hence the stresses) are NOT constant within an element
Computation of the terms in the stiffness matrix of 2D elements (recap)
v4
v3
4
The B-matrix (strain-displacement) corresponding to this element is
3
u4
u3
u1
v
y
v1 (x,y)
v2
u
2
1 u1
u2
 N1 (x,y)

x


0


 N1 (x,y)

y

v1
u2
v2
0
N 2 (x,y)
x
0
N1 (x,y)
y
N1 (x,y)
x
0
N 2 (x,y)
y
N 2 (x,y)
y
N 2 (x,y)
x
u3
N 3 (x,y)
x
0
N 3 (x,y)
y
x
We will denote the columns of the B-matrix as


 N1 (x,y) 


0




x


 N1 (x,y)  ; and so on...
B u1  
0
 ; B v1  

y
 N (x,y) 


1



N
(x,y)


1
y




x
v3
0
N 3 (x,y)
y
N 3 (x,y)
x
u4
v4
N 4 (x,y)
x



N 4 (x,y) 

y

N 4 (x,y) 

x

0
N 4 (x,y)
y
0
The stiffness matrix corresponding to this element is
k   e B D B dV
T
which has the following form
V
u1
v1
 k11
k
 21
 k31

k
k   41
 k51

 k61
k
 71

 k81
k12
k 22
k32
k 42
k52
k62
k72
k82
u2
k13
k 23
k33
k 43
k53
k63
k73
k83
v2
u3
k14
k 24
k34
k 44
k54
k64
k74
k84
k15
k 25
k35
k 45
k55
k65
k75
k85
v3
k16
k 26
k36
k 46
k56
k66
k76
k86
u4
v4
k17
k 27
k 37
k 47
k57
k 67
k 77
k87
k18 
k 28 

k 38 

k 48 
k58 

k 68 
k 78 

k88 

u1
v1
u2
v2
u3
v3
u4
v4
The individual entries of the stiffness matrix may be computed as follows
k11   e Bu1 D Bu1 dV; k12   e Bu1 D Bv1 dV; k13   e Bu1 D Bu2 dV,...
T
T
V
V
V
k21   e Bv1 D Bu1 dV; k21   e Bv1 D Bv1 dV;.....
T
V
T
V
T
Notice that these formulae are quite general (apply to all kinds
of finite elements, CST, quadrilateral, etc) since we have not
used any specific shape functions for their derivation.
Example
1000 lb
300 psi
y
3
4
2 in
1
2
Thickness (t) = 0.5 in
E= 30×106 psi
n=0.25
x
3 in
(a) Compute the unknown nodal displacements.
(b) Compute the stresses in the two elements.
This is exactly the same problem that we solved in last class, except
now we have to use a single 4-noded element
Realize that this is a plane stress problem and therefore we need to use


1
n
0
3.2 0.8 0 


E
0.8 3.2 0  107 psi


D
n
1
0


1 n 2 
1 n  
0 1.2
0 0
  0
2 

Write down the shape functions
1
x  x2  y  y4   ( x  3)( y  2)
N1 
4ab
6
1
x  x1  y  y3    x( y  2)
N2  
4ab
6
1
x  x4  y  y2   xy
N3 
4ab
6
1
x  x3  y  y1    ( x  3) y
N4  
4ab
6
x
y
0
0
3
0
3
2
0
2
We have 4 nodes with 2 dofs per node=8dofs. However, 5 of these are fixed.
The nonzero displacements are
u2
u3
v3
Hence we need to solve
u2
u3
v3
 k11
k
 21
k31
k12
k13  u2   0 
   
k 23  u3    0 
  
k33  
 v3   f 3 y 
k 22
k32
Need to compute only the relevant terms in the stiffness matrix
k11   e Bu2 D Bu2 dV; k12   e Bu2 D Bu3 dV; k13   e Bu2 D Bv3 dV
T
V
T
V
T
V
k21   e Bu3 D Bu2 dV; k22   e Bu3 D Bu3 dV; k13   e Bu3 D Bv3 dV
T
V
T
V
T
V
k31   e Bv3 D Bu2 dV; k22   e Bv3 D Bu3 dV; k13   e Bv3 D Bv3 dV
T
V
T
V
T
V
Compute only the relevant columns of the B matrix
B u2
B u3
 N 2   (2  y ) 
 x   6 

 

 0  0 
 N 2    x 
 y  
6 

 

 N 3   y 
 x   6 

  
  0   0
 N 3   x 
 y   6 

  

  
 0  0
 N   x 

  
B v3   3    
 y   6 
 N 3   y 

6

 x 
 
k11   e Bu2 D Bu2 dV
T
V
2 y

8
7
5 2
(0.1067

10

0.533

10
)(
)

3.33

10
x  dxdy


6

x 0 y 0 
3
 0.5 
2
 0.656 107
Similarly compute the other terms
How do we compute f3y
f3 y  1000 f S3 y
f S3 y  t 
3
x 0
N3 edge  N3
4 3
N 3 along (300) dx
y 2
x
 xy 
  
 6  y 2 3
edge 3 4
x
 (0.5)(300) 
dx
x 0 3
3
 150
2
 225 lb
3
 f3 y  1000 f S3 y  1225lb
4
3
How about a 9-noded rectangle?
Corner nodes
y 5
a
a
2
b
6
b
3
9
7
 x(a  x)   y(b  y ) 
 x(a  x)   y(b  y) 
N1  
N


2
2
2




2a 2   2b 2 
 2a   2b 

 x(a  x)   y(b  y) 
 x(a  x)   y(b  y) 
N 3  

N

4
2
2


 2a 2   2b 2 
2
a
2
b






1
8
4
x
Midside nodes
2
2
 a 2  x 2   y (b  y ) 
 x(a  x)   b  y 
N5  
N 6  



2
2

a
2
b
2a 2   b 2 




2
2
 a 2  x 2   y (b  y ) 
 x(a  x)   b  y 
N7  
N8  
 


2
2
2
2

a
2
b

 2a   b



Center node
 a 2  x 2   b2  y 2 
N9  


2
2
 a
  b

Question: Can you generate the shape functions of a 16-noded rectangle?
Note: These elements, whose shape functions are generated by multiplying the
shape functions of 1D elements, are said to belong to the “Lagrange” family

similar documents