Rayleigh Ritz
Rayleigh Ritz
Rayleigh Ritz
INTRODUCTION
We learned Direct Stiffness Method in Chapter 2
Limited to simple elements such as 1D bars
Strain energy
Beam
Interpolation Nodal
deflection
function
DOF
BEAM THEORY
Euler-Bernoulli Beam Theory
y
Neutral axis
Plane of loading
x
z
A
dv
u ( x, y ) u0 ( x) y
dx
u du0
d 2v
xx
y 2
x dx
dx
dv
dx
y(dv/dx)
Neutral axis
x
L
= dv/dx
y
v(x)
u du0
d 2v
xx
y 2
x dx
dx
d 2v
xx E
xx E
0 Ey 2
dx
P EA
0
d 2v
M EI 2
dx
d v 2
y dA
dx 2 A
d 2v
M EI 2
dx
Sign convention
+M
+Vy
y
x
+P
+M
+P
+Vy
C1
F1
C2
F2
C3
F3
GOVERNING EQUATIONS
Beam equilibrium equations
dVy
0
(
)
f
p
x
dx
V
y
y dx dx Vy 0
dVy
dx
dM
dx
dx pdx Vy dx 0
M M
dx
2
p( x)
Vy
dM
dx
d 4v
p
Vy
dVy
dx
Vy
dx
dM
dx
dx
dx
xx ( x, y )
d 2v
M EI 2
dx
M ( x) y
I
Bending stress
u v
v v
0
y x
x x
u ( x, y ) u0 ( x) y
dv
dx
POTENTIAL ENERGY
Potential energy U V
Strain energy
Strain energy density
2
1
1
1
d 2v
1 2 d 2v
2
U 0 xx
xx E (
xx ) E y 2 Ey 2
2
2
2
dx
2
dx
1 2 d 2v
1 d 2v
U L ( x) U 0 ( x, y, z ) dA Ey 2 dA E 2
2
2 dx
dx
A
A
1 d v
U L ( x) EI 2
2 dx
2
y dA
2
Moment of
inertia
Strain energy
2
1 L d 2v
U U L ( x) dx EI 2 dx
0
2 0
dx
L
NC
i 1
i 1
V p( x)v( x) dx Fv
i ( xi ) Ci
L
dv( xi )
dx
Potential energy
2
NC
NF
L
dv( xi )
1 L d 2v
U V EI 2 dx p( x)v( x) dx Fv
(
x
)
C
i
i
i
0
2 0 dx
dx
i 1
i 1
0
v
v*
10
RAYLEIGH-RITZ METHOD
1. Assume a deflection shape
v( x) c1 f1 ( x) c2 f 2 ( x)..... cn f n ( x)
(c1 , c2 ,...cn ) U V
3. Apply the principle of minimum potential energy to determine
the coefficients
11
x
L
E,I,L
Strain energy
2
1 L d 2v
C 2 EI 4
U EI 2 dx
2 0
4 L3
dx
p0 L4
76.5EI
Cexact
p0 L4
76.8EI
p0 L
p
x 0 x2
2
2
pL
Vy ( x) 0 p0 x
2
M ( x)
13
Deflection
0.8
v(x)/v_max
0.6
0.4
v-exact
0.2
v-approx.
0.0
0
0.2
0.4
Bending
moment
0.6
0.8
0.2
0.4
0.6
0.8
-0.04
-0.06
-0.08
-0.10
M_exact
-0.12
M_approx
Error increases
0.00
-0.02
x
x
-0.14
0.6
V_exact
Shear force
0.4
V_approx
0.2
0.0
-0.2
-0.4
-0.6
0
0.2
0.4
0.6
0.8
14
Assumed deflection
v( x) a bx c1x 2 c2 x3
E,I,L
Need to satisfy BC
v(0) 0, dv(0) / dx 0
v( x) c1 x 2 c2 x3
L
Strain energy U EI 2c1 6c2 x 2 dx
2 0
Potential of loads
L
V c1 , c2 p0 v( x)dx Fv( L) C
0
dv
( L)
dx
p0 L3
p0 L4
2
c1
FL 2CL c2
FL3 3CL2
3
4
15
c1
U
6 EI 2c1 6c2 x xdx EI 6 L2c1 12 L3c2
c2
0
L
PMPE:
0
c1
0
c2
p0 L3
EI 4 Lc1 6 L c2
FL2 2CL
3
p0 L4
2
3
EI 6 L c1 12 L c2
FL3 3CL2
4
2
1
5400 x 2 800 x3 300 x 4
24 EI
16
0.0
0.0
0.0
v-exact
0.0
v-approx.
0.0
0
0.2
0.4
0.6
0.8
Error increases
Bending
moment
500.00
M_exact
400.00
M_approx
300.00
200.00
100.00
0.00
-100.00
0
0.2
0.4
0.6
0.8
Shear force
600.0
500.0
400.0
300.0
V_exact
V_approx
200.0
0
0.2
0.4
0.6
0.8
17
Beam element
C2
p(x)
18
2}T
Scaling parameter s
Length L of the beam is scaled to 1 using scaling parameter s
x x1
1
, ds dx,
L
L
ds 1
dx Lds,
dx L
s
v2
v1
x
L
x1
s=0
x2
s=1
v(0) v1
dx
dv(0)
1
dx
ds dx
v(1) v2
v1 v(0) a0
1
dv
1 (0) a1
dx
L
v2 v(1) a0 a1 a2 a3
2
1
dv
(1) (a1 2a2 3a3 )
dx
L
a0 v1
a1 L
1
a2 3v1 2 L
1 3v2 L
2
a3 2v1 L
1 2v2 L
2
20
Shape functions
1.0
N1 ( s ) 1 3s 2 2 s 3
0.8
N 2 ( s ) L( s 2 s 2 s 3 )
0.6
N 3 ( s ) 3s 2 2 s 3
0.4
N 4 ( s ) L( s 2 s 3 )
0.2
Hermite polynomials
Interpolation property
N1
N3
N2/L
0.0
0.0
0.2
0.4
0.6
N4/L
0.8
1.0
-0.2
21
Approximation of curvature
Curvature is second derivative and related to strain and stress v1
!
!
d 2v 1 d 2v 1
! 1!
[
6
12
s
,
L
(
4
6
s
),
6
12
s
,
L
(
2
6
s
)]
" #
dx 2 L2 ds 2 L2
! v2 !
!
d 2v 1
$
2 !
%
B
q
{
}
2
2
B:
strain-displacement
vector
dx
L 14 41
dx 2
q {B }
L2 14 41
22
Linear
dM
d 3v EI
Vy
EI 3 3 [12 6 L 12 6 L]{q}
dx
dx
L
Constant
23
EXAMPLE INTERPOLATION
v1
v2
Cantilevered beam
Given nodal DOFs
1
{q} {0, 0, 0.1, 0.2}T
L
Deflection and slope at x = 0.5L
Parameter s = 0.5 at x = 0.5L
1
L
1
L
Shape functions: N1 ( 12 ) , N 2 ( 12 ) , N3 ( 12 ) , N 4 ( 12 )
2
8
2
8
Deflection at s = 0.5:
v( 12 ) N1 ( 12 )v1 N 2 ( 12 )
1 N 3 ( 12 )v2 N 4 ( 12 )
2
L
L
v L
1
1
0 0 v2
2 2 2 0.025
2
Slope at s = 0.5:
dN
dv 1 dv 1 dN1
dN
dN 4
v1
1 2 v2 3
2
dx L ds L ds
ds
ds
ds
1
1
v1 (6s 6s 2 )
1 1 4 s 3s 2 v2 (6 s 6 s 2 )
2 2 s 3s 2 0.1
24
L
L
EXAMPLE
A beam finite element with length L
L3
L2
v1 0,
1 0, v2
,
2
3EI
2 EI
Calculate v(s)
v( s) N1 ( s)v1 N 2 (s)
1 N3 (s)v2 N 4 (s)
2
Bending moment
d 2v EI d 2v EI
M ( s ) EI 2 2 2 2 & (6 12s )v2 L(2 6s )
2 '
dx
L ds
L
EI (
L3
L2 )
2 (6 12s )
L(2 6s)
3EI
2 EI
L
L(1 s ) ( L x) Bending moment cause by unit force at the tip
25
Strain energy
Sum of each elements strain energy
LT
NEL
e1
e
x2
NEL
e1
e
e
x2
EI e
x1
1 d 2v
EI
dx 3
2
L
2 dx
C1
p(x)
x
1
F1
x 11
1 d 2v
0 2 ds 2 ds
1
x 21
F2
C2
C3
C4
x12
x 22
F3
x1 3
x 23
F4
x1 4
C5
5
F5
x 24
26
d 2v d 2v d 2v
T
e
e T
q
{
} B B {q }
ds 2 ds 2 ds 2
14
41
41
14
Approximate element strain energy in terms of nodal DOFs
(e)
1 e ( EI
{q }T 3
2
L
e
L(4 6 s) 6 12 s
L(2 6 s) ' ds
27
U U
(e)
e1
1 NEL e T e e
{q } [k ]{q }
2 e1
Assembly
1
U {Q s }T [K s ]{Q s }
2
28
EXAMPLE ASSEMBLY
2EI
1
2L
v1
Two elements
Global DOFs
EI
{Qs }T {v1
1 v2
2 v3
3}
v2
F2
F3
v2
v3
6 L 12 6 L ) v2
( 12
2
2
EI 6 L 4 L 6 L 2 L
2
2
[k ] 3
L 12 6 L 12 6 L v3
2
2
6 L 2 L 6 L 4 L
3
3 3L ) v1
( 3 3L
2
2
EI 3L 4 L 3L 2 L
1
1
[k ] 3
3
3L v2
L 3 3L
2
2
3L 2 L 3L 4 L
2
0
0 )
3 3L
( 3 3L
3L 4 L2 3L 2 L2
0
0
3
3
L
15
3
L
12
6
L
EI
[K s ] 3
L 3L 2 L2 3L 8 L2 6 L 2 L2
0
0
12 6 L 12 6 L
0
6 L 2 L2 6 L 4 L2
0
29
F1
!C !
Concentrated forces and couples
ND
!! 1 !!
V Fi vi Ci
i
V v1
1 v2 ......
ND " F2 # {Q s }T {Fs }
i 1
!
!
!
!
$CND !%
Distributed load (Work-equivalent nodal forces) !
NEL
e
x2
NEL
V e p( x)v( x) dx V
e1
(e)
L
(e)
x1
(e)
(e)
e1
p(s) v N
1
e
x2
e p( x)v( x) dx L
x1
(e)
p(s)v(s) ds
0
1 N 2 v2 N 3
2 N 4 ds
1
1
1
1
v1 L( e ) p ( s ) N1ds
1 L( e ) p ( s ) N 2 ds v2 L( e ) p (s ) N 3ds
2 L( e ) p (s ) N 4ds
0
0
0
0
(e)
(e)
(e)
(e)
v1F1
1C1 v2 F2
2C2
30
F1 pL N1 ( s) ds pL (1 3s 2 2s 3 ) ds
pL
{F} "
$ 2
T
pL2
12
pL
2
pL2
#
12 %
p
pL/2
Equivalent
pL/2
pL2/12
pL2/12
31
6 L 12 6 L ) v1
F1
( 12
pL / 2
! ! 2
2
2 !
! ! !
EI 6 L 4 L 6 L 2 L !
1 ! ! pL /12 ! ! C1 !
" #"
# " #
3
L 12 6 L 12 6 L ! v2 ! ! pL / 2 ! ! F2 !
2
2
2
!
!
!
L
L
L
L
pL
6
2
6
4
/12 %! !$C2 !%
$ 2% $
One beam element has four variables
When there is no distributed load, p = 0
Applying boundary conditions is identical to truss element
PMPE
Potential energy has its minimum when
[K s ]{Q s } {Fs }
Applying BC
The same procedure with truss elements (striking-the-rows and
striking-he-columns)
[K ]{Q} {F}
v1
!
!
dM
d 3v
EI d 3v EI
! !
EI 3 3 3 3 [12 6 L 12 6 L] " 1 #
Vy ( s )
dx
dx
L ds
L
! v2 !
!
Constant
$
2 !
%
Shear force
When true moment is not linear and true shear is not constant, many
elements should be used to approximate it
My
Bending stress x
I
Shear stress for rectangular section
1.5Vy 4 y 2
xy ( y )
1 2
bh
h
34
y
EXAMPLE CLAMPED-CLAMPED
BEAM
v1
v2
x
1
2
1m
6 12 6 ) v1
( 12
6
4 6 2
1
(1)
[k ] 1000
12 6 12 6 v2
2 6 4
2
6
1m
F2 = 240 N
v2
v3
6 12 6 ) v2
( 12
6
4 6 2
2
(2)
[k ] 1000
12 6 12 6 v3
2 6 4
3
6
6 12 6
0
0 ) v1
F1
( 12
! !
6
4 6 2
0
0 !
1 ! ! C1 !
!
!
!
!
!240 !
12 6 24 0 12 6 v2
1000
" # " 0 #
2
0
8 6 2 !
2 ! !
!
6
!
!
!
0
F3 !
0 12 6 12 6 v3
!
! ! !
0
6
2 6 4 $
3 % $ C3 %
0
35
v2 0.01
2 0.0
At x = 0.5
( 12 )
dN3
1
v
2
L(1)
ds
0.015rad
s 12
At x = 1.0
(1)
dN 3
1
v
(1) 2
L
ds
0.0 rad
s 1
(0)
1
dN1
v
(2) 2
L
ds
0.0 rad
s 0
36
p0 = 120 N/m
EI = 1000 N-m2
L = 1m
6 12 6 ) v1
( 12
6
4 6 2
1
[K s ] 1000
12 6 12 6 v2
2 6 4
2
6
C = 50 N-m
Applying BC
(12 6 ) v2
60
1000
#
" # "
6 4 $
2 % $60 %
v2 0.01m
2 0.03 rad
38
F1 120 N
1000 6v2 2
2 C1 10
C1 10 N m
Bending moment
EI
B {q}
L2
EI
2 & (6 12 s)v1 L(4 6 s)
1 (6 12 s)v2 L(2 6 s)
2 '
L
1000[0.01(6 12s ) 0.03(2 6 s)]
60s N m
M ( s)
Shear force
EI
&12v1 6 L
1 12v2 6 L
2 '
L3
1000[12 (0.01) 6(0.03)]
60 N
Vy
39
0.000
FEM
Exact
-0.002
FEM
Exact
-0.005
-0.010
-0.004
-0.015
-0.006
-0.020
-0.008
-0.025
Deflection
0
0.2
0.4
0.6
0.8
10
0.2
0.4
0.6
0.8
0
FEM
Exact
FEM
Exact
-20
-10
-40
Vy
-20
Slope
-0.030
-0.010
-30
-60
-80
-40
-100
Bending moment
-50
-60
Shear force
-120
0
0.2
0.4
0.6
0.8
0.2
0.4
0.6
0.8
40
Frame structure
Can carry axial force, transverse shear force, and bending moment
(Beam + Truss)
Assumption
v1
u2
1
2
p
F
v2
u2
v2
u1
u1
v1
v2
u1
v1
41
sin * 0
0
cos * 0
0
0
1
0
0
0 cos *
0
0 sin *
0
0
0
0
0 ) f x1
0
0 ! f y1 !
! !
0
0 ! c1 !
" #
sin * 0 ! f x 2 !
cos * 0 ! f y 2 !
! !
0
1 $ c2 %
Local coordinates
{f } [T]{f }
{q} [T]{q}
v1
y
u1
1
v2
u2
1
Global coordinates
42
Beam bending
6 L 12 6 L ) v1
f y1
( 12
! ! !
2
2 !
EI 6 L 4 L 6 L 2 L !
1 ! ! c1 !
" #" #
L3 12 6 L 12 6 L ! v2 ! ! f y 2 !
2
2
6 L 2 L 6 L 4 L $!
2 !% !$ c2 !%
43
0
6 La2
4 L2 a2
0
6 La2
2 L2 a2
a1
0
0 12a2
0 6 La2
a1
0
0
12a2
0 6 La2
0 ) u1
f x1
! ! ! !
6 La2 ! v1 ! ! f y1 !
2 L2 a2 !
1 ! ! c1 !
" # " #
0 ! u2 ! ! f x 2 !
6 La2 ! v2 ! ! f y 2 !
! ! ! !
4 L2 a2 $
2 % $ c2 %
EA
L
EI
a2 3
L
a1
[k ]{q} {f }
[T]T [k ][T]{q} {f }
[k ]{q} {f }
[k ] [T]T [k ][T]
AE
u2 u1
L
EI
Bending moment M ( s) 2 B {q}
L
EI
Shear force Vy ( s ) 3 [12 6 L 12
L
Other method: V
6L
( 12
y1
!
!
2
! M 1 ! EI 6 L 4 L
"
# 3
V
! y 2 ! L 12 6 L
2
!$ M 2 !%
6L 2L
Axial force P
6 L]+q,
12 6 L ) v1
6 L 2 L2 !!
1 !!
" #
12 6 L ! v2 !
6 L 4 L2 !$
2 !%
45