FEM C4 Dam English PDF

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

Chapter 4

BEAM

BEAM ELEMENT

ConsiderChuyeå n vò theobeam
the following phöông vuoâwhich
element ng goùchastruïtwo
c y
degrees of freedom
thanh (ñoä voõng): (deflection and rotation) at each node v1 q1 v2 q3
as shown:
v  x  a1  a2 x  a3 x 2  a4 x 3 e
x
 1 x x2 x 3  a   P  x  a 1 q2 2 q4
  N1 N2 N3 N 4 qe   N qe L,EJ

 x2 x3 x2 x3
 1N  1  3  2 N 3  3  2
 L2 L3 L2 L3
in which:  2 2
 N  x 1  2 x  x  N  x   x  x 
 2  2  4  2 
 L L   L L 

1
BEAM ELEMENT

From Mechanics of
Materials and Structures
A' (Figure 5), we have the
y v
x relationship between the
B' axial displacement and
deflection:
dv
u  y
v dx
A du d 2v d2
y Thus:  x    y 2   y 2  N qe
dx dx  dx

B
 x   B qe

u=-y  v
x
Figure 5: Bernuli-Euler beam in flexure

BEAM ELEMENT

d2  6 12 x 4 6x 6 12 x 2 6x 
We have: 
B   y  N    y   2  3     2   2  3     2 
dx 2  L L   L L   L L   L L 
14 

Beam in flexure:     x  ,     x  ,  D  E


12 6L 12 6L
– The stiffness matrix of an element:
T
L
T EJ z 4 L2 6 L 2 L2
  K e    B   D   B  dV  E    B   B  dFdx 
 44  Ve  41 11 14  0 F
L3 12 6 L
– Nodal load vector of an element: ñx 4 L2
L nQ
T T
Pe    N  q  x  dx    N  xQ  
 i
Qi y
0 i 1
 41 M1
T P1 Q1 P3
nM
 dN  p(x)
   
xM i  M i
i 1  dx  x
P2 e P4
Qi , xQi are point load and its coordinate. xQ
1
xM1
Mi , xMi are moment and its coordinate.

nQ , nM are the number of point load and moment on the element.

2
BEAM ELEMENT

 Distributed load: p  x   const  p0

 x2 x3   po L 
1  3 2  2 3   2 
 L L 
2 3
 2

 x x   p0 L 
L x  2 
L L2   12 
Pe   2 3
 p0 dx   
0 
3
x
 2
x   po L 
 L2 L3   2 
 2 3
  2
 x x    p0 L 
 L L2   12 

BEAM ELEMENT

 Point load Q as shown in the following Figure (xQ = a)


 a2 a3  y
 1  3  2 
 L2 L3  P1 Q P3
 a 2
a   P1 
3

 a2  2   
T  L L   P2  x
Pe   N  a  Q  Q  2 3 
  P2 e P4
 3 a  2 a   P3 
 L2 x Q =a
L3   P4 
 2 3  L
 a a 
2
 L L 
 If the point load P is located at the middle and downward:

3
BEAM ELEMENT

d 2v d2
Determine the internal moment: M  EJ  EJ  N qe  EJ  N qe  Baäc I
dx 2 dx 2
 M  at node n1  
Set: M e    is the bending moment vector at two nodes of an element
 M  at node n2  

  N   x  0   
Then: M e  EJ    qe   S e qe
 21   N   x  L  
  24   41

EJ  6 L 4 L2 6 L 2 L2 
in which:  S e    is the matrix for calculating moment.
L3  6 L 2 L2 6 L 4 L2 

Note: M  EJ  N qeis moment due to nodal displacement. We should add the


q
moment (M0) due to load on element with fixed ends.

qL2 qL2 qL2


12 8 12

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy