ELEMENTO FINITO TAREA 2

Ejercicio 3
Datos

9 kg No. de elementos
Ea   
20010 Pa b   7800
3
m Ne  6
es  0.5in
Grados de libertad
de   12in
Gn  3
PV   350kgf
Número de nodos
Lb   3m
Nn  7
di   11in  4  4
1  de   di  
  2 2 Dim   Gn Nn Ib         
At    de  di 4  2   2  
4  


Ee   Ea Ea Ea Ea Ea Ea T 
Ie   Ib Ib Ib Ib Ib Ib T
T
Ae   ( At At At At At At )

Le    1.2Lb 1.5Lb Lb 2 Lb Lb Lb T

Cómputo de Matrices
colocarlos de menor a mayor
Matrices posicion

 0   0 
1 1 2  
     1.2Lb   0 
2 3 2    
     1.2Lb   1.5Lb 
3 3 4    
y    1.5Lb 
ele   P ele  
4 4 5 x    2.2Lb 
   
 b 
   1.5L 
5 5 6  4.2Lb 
     4.2L   0.5L 
6 6 7
 b
 b
 5.2L   0.5Lb 
 b  

 
MP P ele  Ne   for i  1  Ne
   1   1   1   2   2   2 
pp( i)  3  P ele   2 3  P ele   1 3  P ele  3  P ele   2 3  P ele   1 3  P ele  
 i i i i i i
for j  1  Ne
 j
A  pp( j )
T  
T
A


GDL   MP P ele  Ne 
Longitudes

2 2
He( i)   x  2    x  1    y   y
 P  1  
  Pele     Pele 2  
 ele  i  ele  i
P
 i  i

Angu los

( i)   ang lex
 P  1  
  x  y y
  P ele 2    P  1    P  2  
 i  ele  i  ele  i  ele  i

 cos( ( i) ) sin( ( i) ) 0 0 0 0
 
 sin( ( i) ) cos( ( i) ) 0 0 0 0
 
R r ( i)  
 0 0 1 0 0 0
 0 0 0 cos( ( i) ) sin( ( i) ) 0
 
 0 0 0 sin( ( i) ) cos( ( i) ) 0
 
 0 0 0 0 0 1

Ie Ee Ae Le
 
KK Ie  Le  Ee  Ae  e   Ie 
4
 Ee 
Pa
 Ae 
2
 Le 
m
m m

for i  1  Ne

 2 2 
  Le  Ae  Lei Aei 
  i  i
0 0 
  0 0 
 Ie
i
Ie
i

 
 0 12 6.Le 0 12 6 Le 
 i i

 2 2
 0 6.Le 4  Le  0 6 Le 2  Le 
i  i i  i 
o( i)  
  L  2 A 2
 Lei Aei 
  ei ei 0 0   0 0 
 Ie Ie 
 i i 
 0 12  6Le 0 12 6 Le 
 i i 
 2 2
 0 6 Le 2  Le  0 6 Le 4  Le  
 i  i i  i 
Ee Ie
i i T
K ( i)  R r ( i) o( i) R r ( i)
3
 Lei
 
K ( e)


Ke( i)   KK Ie  Le  Ee  Ae  i 
 e Ae Le
MM  eLe Aee   e 
kg
 Ae 
2
 Le 
m
m
3
m
for i  1  Ne

140 0 0 70 0 0 
 
 0 156 22Le 0 54 13Le 
 i i

 2 
e Ae Le  0 22Lei 4  Lei 0 13Le 3 Le 
i i i   i i
p( i)   
420  70 0 0 140 0 0 
 0 54 13Le 0 156 22Le 
 i i 
 
 0 13Le 3 Le 2
0 22Le 4  Le  
 i i i  i 
T
M( i)  Rr ( i) p( i) Rr ( i)

M( e)

Me( i)   MM  e Le  Ae  i
Definiciones de Fuerza axiales Definiciones de Fuerza flectoras verticales

 0   0 
 
 e g Ae Le   
 2 2 2
  0 
   
0 0 
P a     P fv   
 0   P V 
  g A L   
 e5 e5 e5   0 
   
 0   0 

Definiciones de Fuerza distribuidas adimensional binario

1
 
0
 
P dis  
1
1
 
0
 
1
Fuerza Distribuida (peso)

Ae Le e
 
FFd P dis  Ae  Le  e  e   Ae 
2
 Le 
m e
  
kg
m
3
m
for i  1  Ne

 1 
 
 Lei 
 
F ( i) 
1 9.8
e  Ae Le  
 6  P
2  i i i  1  disi
 
 Lei 
 
 6 
F ( e)


F d( i)   FFd P dis  Ae  Le  e  i 

Fuerza axial

Pa
 
FFa P a  e   Pa 
N
ojo con la tranformación (giro)
for i  1  Ne

  0.5
F ( i)  P a  
i   0.5
F ( e)


F a ( i)   FFa P a  i 

T
 1 3  T
 2  3     1 2 
 4   3  3 
 L   4 
 ei   L 
1  
ei
 2 1 2 3 
   1          Lei   2 1  2 
   1  2  3  
   2   m 4
m 4
He(  i)   
 1 3  Hpe(  i)     
  2  3      m   1 2 
 4    3  3  
 Le   4 
 i   Le 
 2 1  i 
2 3   2 1
 m  4  1         2 
   m  4  1  2  3  
 

Fuerza flectoras verticales

F ( i)  He  0 iT  Pfv
i
F ( e)


F c( i)   FFc P fv i 
Unificacion del Vector Fueza
F t F a  F c  F d  e   for i  1  Ne

 F a ( i)
1

 F ( i)  F ( i) 
 c 1 d 1
 
 Fc( i) 2  Fd( i) 2 
F p( i)  
 F a ( i)
2 
 
 Fc( i) 3  Fd( i) 3 
 
 Fc( i) 4  Fd( i) 4 
 
F p( e)


F to( i)   F t F a  F c  F d  i 
Ensamblaje de Matices


ENSK GDL K.e  Dim    for c1  1  Dim

for f1  1  Dim

( Kc) 0
f1
 c1
KKm  Kc

for i  1  Ne

for f  1  6

for c  1  6

   GDL c  
 i 
Aux  KKm    GDL f  
i




 
 c 
 GDL i    GDL f  
 c
KKm  Aux  
 K.e( i)  f

i
KKm


KK   ENSK GDL  Ke  Dim 
 
ENSM GDL Ke  Dim   for c1  1  Dim

for f1  1  Dim

( Kc) 0
f1
 c1
KKm  Kc

for i  1  Ne

for f  1  6

for c  1  6

   GDL c  
 i 
Aux  KKm    GDL f  
i

   GDL c  
 i 
  GDL f  
 c
KKm   Aux  
 Ke( i)  f

i
KKm


MM   ENSK GDL  Me  Dim 


ENSF GDL F to  Dim   for c1  1  Dim

for f1  1  Dim

( Fc )  0
f1
for i  1  Ne

for f  1  6
Aux  Fc
 GDL f   i
Fc
 GDL f   i   Fto( i)  f  Aux
Fc


F F   ENSF GDL  F to  Dim 

Leyenda
Vector selector de
Penalizacion 0= no hay desplazamientos
desplazamientos
1 = tiene desplazamiento

q  0 0 0 1 1 1 1 1 1 1 0 1 1 1 1 1 1 1 0 0 0
T

5
Kc   10  
ma xKK
 
P EN KK  Dim q    for c1  1  Dim

for f1  1  Dim

( P c)  0
f1
 c1
P EN  Pc
for i  1  Dim

 P EN i  i  105 m a xK
 K   qi  1
 
P EN  KK  P EN

P EN

1 2 3 4 5 6
 
7 8 9 4 5 6
 

KF   PEN KK  Dim  q  7 8 9 10 11 12 
GDL 
 10 11 12 13 14 15 
 
1  13 14 15 16 17 18 
X   KF F F  
 16 17 18 19 20 21 

1
1 -3.412·10 -11
2 3.051·10-11
3 2.313·10-11  X1   X2   X3 
X  X  X 
4 6.141·10-6  4  5  6
5 1.586·10-4      
 X7   X8   X9 
6 2.574·10-4      
7 1.725·10-3 u    X10  m     X11  m d    X12  rad
     
X 8 1.579·10-4  X13   X14   X15 
     
9 2.395·10-4  X16   X17   X18 
10 1.725·10-3      
 X19   X20   X21 
11 7.385·10-11      
12 -4.043·10 -4
13 1.725·10-3
14 -1.312·10 -3
15 4.233·10-4
16 ...

 8   8 
 3.412  10  3.051  10   9 
 1.325  10 
   
 
 6.141  10 3   0.159 
 0.015 
   
 0.158   
 1.725  
0.014

u 1.725  m m    7.385  10 8  mm d   0.023  de g
   
 1.312   
 1.725  0.024
     
 3.488  10 3   1.308   0.037 
     
 8   4.434  10 8   9 
   5.772  10 
 2.325  10 
 N
Ftoo o    1      Kc 
 4 7 m

4
Ftoo o  1.733  10 N

Fprueba  e g Ae 7.7Lb  PV

4
Fp ru eba  2.402  10 N