Frame Element PDF

Download as pdf or txt
Download as pdf or txt
You are on page 1of 7

Frame Element Stiffness Matrices

CEE 421L. Matrix Structural Analysis


Department of Civil and Environmental Engineering
Duke University

Henri P. Gavin
Fall, 2012

Truss elements carry axial forces only. Beam elements carry shear forces
and bending moments. Frame elements carry shear forces, bending moments,
and axial forces. This document picks up with the previously-derived truss
and beam element stiffness matrices in local element coordinates and proceeds
through frame element stiffness matrices in global coordinates.

1 Frame Element Stiffness Matrix in Local Coordinates, k

A frame element is a combination of a truss element and a beam element.


The forces and displacements in the local axial direction are independent of the
shear forces and bending moments.

EA
      





N1 









q1 



  L 


 u1 




 
 
 
   


 
 
 
   


 
 
 
   

      

12EI
     
V1 q2 0 u2

 
 
 
   
      

L3

 
 
 
   


 
 
 
   


 
 
 
   






















 sym 







 
    
6EI 4EI  
M1 q3 0 u3

 
 
 
   

L2
     
L

 
 
 
   

 
 
 
   

 
= =

 
     
− EA EA
      
N2 q4 u4
   
0 0 
      


 
 
 
   


















 L L 








 
 
 
   


 
 
 
   


 
 
 
   

− 12EI 6EI 12EI
     
V2 q5 0 − 0 u5

 
 
 
   
      












  L3 L2 L3  


 
 
 
   


 
 
 
   


 
 
 
   


 
 
 
   

 
6EI 2EI
0 − 6EI 4EI
      
M2 q6 0 u6

 
 
 
 
 

L2 L2
     
L L
2 CEE 421L. Matrix Structural Analysis – Duke University – Fall 2012 – H.P. Gavin

2 Relationships between Local Coordinates and Global Coordinates: T

The geometric relationship between local displacements, u, and global dis-


placements, v, is

u1 = v1 cos θ + v2 sin θ u2 = −v1 sin θ + v2 cos θ u3 = v3

or, u = T v.
The equilibrium relationship between local forces, q, and global forces, f ,
is

q1 = f1 cos θ + f2 sin θ q2 = −f1 sin θ + f2 cos θ q3 = f3

or, q = T f , where, in both cases,


 
c s 0 x2 − x1
c = cos θ =
 



−s c 0 0 

 L

0 0 1 
T=
 
 

 c s 0 

y2 − y 1
 



0 −s c 0 

 s = sin θ =
0 0 1 L

The coordinate transformation matrix, T, is orthogonal, T−1 = TT .

CC BY-NC-ND H.P. Gavin


Frame Element Stiffness Matrices 3

3 Frame Element Stiffness Matrix in Global Coordinates: K

Combining the coordinate transformation relationships,

q = ku
Tf = kTv
f = TT k T v
f = Kv

which provides the force-deflection relationships in global coordinates. The


stiffness matrix in global coordinates is K = TT k T


EA 2 EA
− EA − EA

2
 L c L cs L c L cs 
12EI 2 12EI



+ L3 s − L3 cs − L2 s − L3 s + L3 cs − 6EI
6EI 12EI 2 12EI
L2 s



 
 
 
EA 2


 L s − EA
L cs − EA
L s
2 


12EI 2 6EI



+ L3 c L2 c + 12EI 12EI 2
L3 cs − L3 c
6EI
L2 c



 
 
 
4EI 6EI


 L L2 s − 6EI
L2 c
2EI
L



K= 





EA 2 EA
L c L cs
 
 
 
12EI 2 12EI 6EI



+ L3 s − L3 cs L2 s






sym 


EA 2
L s
 
 
 
12EI 2



+ L3 c − 6EI
L2 c



 
 
 
4EI
L

CC BY-NC-ND H.P. Gavin


4 CEE 421L. Matrix Structural Analysis – Duke University – Fall 2012 – H.P. Gavin

4 Frame Element Stiffness Matrices for Elements with End-Releases

Some elements in a frame may not be fixed at both ends. For example,
an element may be fixed at one end and pinned at the other. Or, the element
may be guided on one end so that the element shear forces at that end are zero.
Or, the frame element may be pinned at both ends, so that it acts like a truss
element. Such modifications to the frame element naturally affect the elements
stiffness matrix.
Consider a frame element in which a set of end-coordinates r are released,
and the goal is to find a stiffness matrix relation for the primary p retained
coordinates. The element end forces at the released coordinates, qr are all zero.
One can partition the element stiffness matrix equation as follows
   

qp 
kpp kpr 
up 
qr krp krr ur

The element displacement coordinates at the released coordinates do not equal


the structural displacements at the collocated structural coordinates, since the
coordinates r are released. Since the element end forces at the released coor-
dinates are all zero (qr = 0), the element end displacements at the released
coordinates must be related to the displacements at the primary (retained)
coordinates as:
−1
ur = −krr krp up

The element stiffness matrix equation relating qp and up is

−1
h i
qp = kpp − kpr krr krp up

The rows and columns of the released element stiffness matrix corresponding
to the released coordinates, r, are set to zero. The rows and columns of the
released element stiffness matrix corresponding to the retained coordinates, p,
−1
are [kpp − kpr krr krp ]. This is the element stiffness matrix that should assemble
into the structural coordinates collocated with the primary (retained) coordi-
nates p. The following sections give examples for pinned-fixed and fixed-pinned
frame elements. Element stiffness matrices for many other end-release cases
can be easily computed.

CC BY-NC-ND H.P. Gavin


Frame Element Stiffness Matrices 5

4.1 Pinned-Fixed Frame Element in Local Coordinates, k . . . (r = 3)


EA
0 − EA

 L
0 L 0 0 
 
 
 
3EI


 L3 0 0 − 3EI
L3
3EI 
L2


 
 
 



0 0 0 0 
k= 





EA


 L 0 0 



sym 


3EI


 L3 − 3EI 
L2 

 
 
 
3EI
L

4.2 Pinned-Fixed Frame Element in Global Coordinates, K = TT kT


EA 2 EA
0 − EA − EA

2 3EI
 L c L cs L c L cs − L2 s 
3EI 2 3EI



+ L3 s − L3 cs − 3EI
L3 s
2
+ 3EI
L3 cs



 
 
 
EA 2


 L s 0 − EA EA 2
L cs − L s
3EI
L2 c



3EI 2



+ L3 c + 3EI 3EI 2
L3 cs − L3 c



 
 
 



0 0 0 0 


K= 





EA 2 EA EI
L c L cs L2 s
 
 
 
3EI 2 3EI



+ L3 s − L3 cs 





sym 


EA 2


 L s − 3EI
L2 c



3EI 2
L3 c
 
 
 
 
 
 
3EI
L

CC BY-NC-ND H.P. Gavin


6 CEE 421L. Matrix Structural Analysis – Duke University – Fall 2012 – H.P. Gavin

4.3 Fixed-Pinned Frame Element in Local Coordinates, k . . . (r = 6)


EA
− EA

 L
0 0 L 0 0 
 
 
 
3EI 3EI


 L3 L2 0 − 3EI
L3 0 


 
 
 
3EI


 L 0 − 3EI
L2 0 


k= 





EA


 L 0 0 





sym 


3EI


 L3 0 


 
 
 
0

4.4 Fixed-Pinned Frame Element in Global Coordinates, K = TT kT


EA 2 EA EA 2
− EA

 L c L cs − 3EI
L2 s − L c L cs 0 
3EI 2 3EI 3EI 2 3EI



+ L3 s − L3 cs − L3 s + L3 cs 


 
 
 
EA 2


 L s
3EI
L2 c − EA EA 2
L cs − L s 0 


3EI 2 3EI 3EI 2



+ L3 c + L3 cs − L3 c 


 
 
 
3EI 3EI


 L L2 s − 3EI
L2 c 0



K= 





EA 2 EA


 L c L cs 0 


3EI 2 3EI



+ L3 s − L3 cs 





sym 


EA 2


 L s 0 


3EI 2



+ L3 c 


 
 
 
0

CC BY-NC-ND H.P. Gavin


Frame Element Stiffness Matrices 7

5 Notation
u = Element deflection vector in the Local coordinate system
q = Element force vector in the Local coordinate system
k = Element stiffness matrix in the Local coordinate system
... q = k u
T = Coordinate Transformation Matrix
... T−1 = TT
v = Element deflection vector in the Global coordinate system
... u = T v
f = Element force vector in the Global coordinate system
... q = T f
K = Element stiffness matrix in the Global coordinate system
... K = TT k T
d = Structural deflection vector in the Global coordinate system
p = Structural load vector in the Global coordinate system
Ks = Structural stiffness matrix in the Global coordinate system
... p = Ks d

Local Global
Element Deflection u v
Element Force q f
Element Stiffness k K
Structural Deflection - d
Structural Loads - p
Structural Stiffness - Ks

For frame element stiffness matrices including shear deformations, see:


J.S. Przemieniecki, Theory of Matrix Structural Analysis, Dover Press, 1985.
(... a steal at $12.95)

CC BY-NC-ND H.P. Gavin

You might also like