Ejercicio 3.

Consideremos la barra mostrada en la figura para cada elemento i, A, L son el área transversal y
longitud, respectivamente. Cada elemento “i” está sometido a una fuerza de tracción “Ti” por
unidad de longitud y una fuerza de cuerpo “f” por unidad de volumen, las unidades Tp, f, Ap
etc, se suponen que son consistentes, el módulo de Young del material es E, Se aplica una carga
concentrada P2 en el nodo 2, ahora se ensamblara la matriz de rigidez estructural y el vector de
carga nodal.

Ensamble de la matriz de rigidez global y del vector de carga

La matriz de rigidez del elemento para cada elemento i se obtiene de la siguiente ecuación:

EA i 1 −1
K ( i )= [
Li −1 1 ]
¿
Tabla de conectividad con el elemento

Elemento 1 2

1 1 2

2 2 3

3 3 4

4 4 5

1
Las matrices de rigidez del elemento pueden desarrollarse usando la tabla de conectividad y
luego sumarse (o ensamblarse) para obtener la matriz de rigidez estructural como sigue:

[ ] [ ]
1 −1 0 0 0 0 0 0 0 0
EA 1 −1 1 0 0 0 EA 2 0 1 −1 0 0
K= 0 0 0 0 0 + 0 −1 1 0 0
L1 L2
0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0

[ ] [ ]
0 0 0 0 0 0 0 0 0 0
+ EA 3 0 0 0 0 0 EA 4 0 0 0 0 0
0 0 1 −1 0 + 0 0 0 0 0
L3 L4
0 0 −1 1 0 0 0 0 1 −1
0 0 0 0 0 0 0 0 −1 1

Resultado:

[ ]
A1 −A1
0 0 0
L1 L1
−A 1 A A −A2
( 1+ 2) 0 0
L1 L 1 L2 L2
−A2 A A −A 3
K=E∗ 0 ( 2+ 3) 0
L2 L 2 L3 L3
− A3 A A −A 4
0 0 ( 3+ 4)
L3 L3 L 4 L4
−A 4 A4
0 0 0
L4 L4

Vector Global:
 A 1 L1 f L1 T 1
+
2 2
¿
A L f LT + A 1 L1 f L2 T 2 A L f L T
¿( 1 1 + 1 1 ) + ( 1 1 + 2 2) A L f L T
2 2 2 2 2 2 4 4
+ 4 4
+ A 3 L3 f L3 T 3 A 3 L 3 f L3 T 3 2 2
+ ¿ ( + )
2 2 2 2
¿

[]
0
P2
F= [
A 4 L 4 f L4 T 4
2
+
2 ]
¿¿ + 0
0
0

Ejercicio 3.3

Considere la placa delgada (acero) en la figura, la placa tiene un espesor uniforme t=1 in,
6 lb
módulo de Young E=30 X 10 psi y densidad por peso p= 0.2836 . Además de su
¿3
propio peso, la placa está sometida a una carga concentrada p= 100 lb en su punto medio.

a) Modele la placa con los elementos finitos

b) Escriba expresiones para las matrices de rigidez del elemento y vectores fuerza de
cuerpo del elemento

c) Ensamble la matriz de rigidez estructural K y el vector de carga global F

d) Usando el método de eliminación, obtenga el vector de desplazamiento global Q

e) Evalué los elementos en cada elemento.

f) Determine la fuerza de reacción en el soporte

Solución:

a) Usando dos elementos, cada uno de 12 in de longitud

 Área promedio del elemento 1:

6+ 4.5
A 1= =5.25 ¿2
2
 Área promedio del elemento 2:

4.5+ 3
A 2= =3.75 ¿2
2

b) Expresiones de matrices de rigidez del elemento para los elementos:

30∗10 6∗5.25 1 −1
K 1=
12 [
−1 1 ]
6
K 2=
12 [
30∗10 ∗3.75 1 −1
−1 1 ]
f 1=
5.25∗12∗0.2836 1
2 1 []
f 2=
3.75∗12∗0.2836 1
2 1 []

c) Matriz de rigidez global K

[ ]
6 5.25 −5.25 0
30∗10
K= −5.25 9.00 −3.75
12
0 −3.75 3.75
 El vector F aplicada extremamente se ensambla a partir f1 y f2 y la carga
puntual p=100 lb

8.9334
¿
¿ 15.3144 +100 6.381
¿ ¿ ¿
F=¿

d) Utilizando el método de eliminación obtenemos la vector de desplazamiento

global Q

30∗106 9
[ ][ ] [
−3.75 ∗ Q2 = 115.3144
12 −3.75 3.75 Q3 6.3810 ]
Q2=0.9272 x 10−5 in

Q3=0.9953 x 10−5 in

Así que Q = [0,0.272∗10−5 ,0.9953∗10−5 ]t ∈¿

e) Obtenemos el esfuerzo en cada elemento:

1
o1 =30∗10 6 x
12 (
[ −1 1 ] 0
0.9272 x 10−5 )
¿ 23.18 psi

( )
−5
1
o2 =30∗10 6 x [ −1 1 ] 0.9272∗10−5
12 0.9953∗10

¿ 1.70 psi

f) Determinando la fuerza de reacción del soporte

−8.9334

( )
6 0
30∗10
R 1= ]
[ 5.25 −5.25 0 0.9272∗10−5
12 −5
0.9953∗10

=-130.6 lb
Nota: Evidentemente la reacción es igual y opuesta a la carga total hacia abajo
que actúa sobre la placa

Ejercicio 1.4
Una placa de acero se somete a una carga axial, la aproximación de las
desviaciones y las tensiones medias a lo largo de la placa, la placa tiene un
espesor de 1/16” con un módulo de elasticidad E = 29 x lo6 Ib/in

A 1 E 5∗0.0625∗29∗10 6
K 1= = =9,062.500lb/¿
L1 1
A 2 E 2∗0.0625∗29∗106
K 2=K 3 = = =906,250 lb / ¿
L2 4
A E 5∗0.0625∗29∗106
K4= 4 = =4,531.250 lb/ ¿
L4 2

Calculando la matriz:

k 1=
[ K 1 −K 1
−K 1 K 1 ]

Matriz de posición global:

[ ]
K 1 −K 1 0 0
k 1 G = −K 1 K 1 0 0
0 0 0 0
0 0 0 0
Matriz de rigidez por elemento:

[ ]
K 1 −K 1 0 0
2G

k 2=
[ K2
−K 2
−K 2
K2 ]
k = −K 1 K 1
0
0
0
0
0
0
0
0
0
0
 [ ]
0 0 0 0
3G

k 3=
[ K3
−K 3
−K 3
K3 ]
k =
0 K 3 −K 3
0 −K 3 K 3
0
0
0 0 0 0

[ ]
0 0 0 0
3G

k 4=
[ K4
−K 4
−K 4
K4 ]
k =
0
0
0
0 0 0
0 K 4 −K 4
0 −K 4 K 4

K G= K 1G + K 2 G + K 3 G + K 4 G

[ ]
K1 −K 1 0 0
−K 1 K 1+ K 2+ K 3 −K 2−K 3 0
K G=
0 −K 2−K 3 K 2 + K 3 + K 4 −K 4
0 0 −K 4 K4

[ ]
9062500 −9062500 0 0
G −9062500 10875000 −18125000 0
K =
0 −18125000 6343750 −4531250
0 0 −4531250 4531250

[][ ]
u1 0
u2 8827∗10−5
= −4
u3 5296∗10
u4 7062∗10−4

Calculando el esfuerzo en cada elemento:

u2−u1 (29∗106)∗( 8.827∗10−5 −0)

1
o =E
L
= ( ) 1
Lb
=2560 2
¿

u3−u 2 (29∗106 )∗( 5.296∗10−5 −8827∗10−5)

2
o =E
L
= ( ) 1
Lb
=3200 2
¿
 u 4−u 3 (29∗106 )∗(7062∗10−5−5296∗10−5 )
o 4=E ( L ) =
1
Lb
=2560 2
¿