“ANÁLISIS ESTRUCTURAL” ANÁLISIS MATRICIAL PARA EL CÁLCULO DE CERCHAS
Presentado por: Juan Carlos Nontoa Andrés Ricardo Rocha Carlos EduardoTorres Presentado a: Ricardo Correa Uribe
UNIVERSIDAD SANTO TOMAS BOGOTA D.C. 21 de Abril de 2010
INTRODUCCION Las armaduras de acero o de distintos tipos de material, constituyen un elemento de gran utilidad dentro del campo de la ingeniería estructural. Su diseño permite distribuir las fuerzas producidas por diferentes cargas a lo largo de su estructura interna y así poder llevarlas a sus respectivos apoyos unas ves definidas. Las diferentes clases de armaduras tienen varios tipos de análisis dependiendo de su diseño y de su función a futuro, en este trabajo de investigación se busca analizar las armaduras por medio de el método matricial de rigidez y así comparar los resultados con otros métodos utilizados dentro del campo de la ingeniería.
OBJETIVOS •
Analizar el método matricial para solucionar armaduras.
•
Desarrollar y explicar paso a paso el proceso de solución de armaduras con el método matricial de rigidez.
•
Solucionar y mostrar diferentes problemas en donde se utiliza el método matricial para armaduras.
METODO MATRICIAL PARA CÁLCULO DE CERCHAS Armaduras Estáticamente Determinadas: Definición: Una armadura es una estructura consistente en un número finito de barras conectadas en uniones por pasadores sin fricción en los cuales pueden aplicarse las fuerzas externas. En la Figura 3.1 se muestra una unión de armadura típica en la situación ideal, donde un pasador se inserta en los extremos de barras de ojo. Las barras tienen libertad de girar sobre el pasador.
Como las fuerzas externas sólo se aplican en los nudos, cada barra en una armadura puede sujetarse solamente a un par de fuerzas axiales iguales y opuestas; de aquí el nombre de miembros de dos fuerzas. Este hecho puede confirmarse considerando las peculiaridades de las barras fuera de las
uniones, como se muestra en la Fig. 3.Ib. Así, las fuerzas internas desconocidas son las fuerzas axiales F¡, para ; = 1 a j = NF, donde NF es el número total de miembros. Para orientar una armadura bidimensional en su propio plano, debe haber al menos un soporte articulado fijo y un soporte de rodillos; o como una alternativa poco usual, tres soportes de rodillos que no tengan más de dos reacciones paralelas. En un soporte articulado fijo, hay dos componentes desconocidas de la reacción de la articulación, y para un soporte de rodillos, hay solamente una reacción desconocida. De aquí que el número total de componentes desconocidas de las reacciones Nfí debe ser igual o mayor que tres. Cualquier sistema de fuerza externa actuando sobre una armadura bidimensional puede describirse por las magnitudes numéricas en las dos direcciones de referencia en los nudos, excepto aquellas a lo largo de las cuales ya hay componentes desconocidas de las reacciones. Si NP es el número total de direcciones posibles en las cuales las fuerzas externas puedan actuar, y N/ el número total de uniones, entonces (3.1.1)
NP + NR = 2(NJ)
Como el número total disponible de ecuaciones de equilibrio estático para una armadura bidimensional es 2 ( NJ ) y el número total de fuerzas des conocidas es (NF + Nñ), el requerimiento necesario para determina ción estática es (3.1.2) NF+NR = 2(NJ) Comparando las Ees. (3.1.1) y (3.1.2), se obtiene (3.1.3) NP = NF La Ec. (3.1.3) prueba que, para armaduras estáticamente determinadas, el número de direcciones posibles en el cual pueden actuar las fuerzas externas es igual al número de barras. Como ilustraciones, las cuatro armaduras de la Fig. 3.2 satisfacen la condición NP = NF, aunque solamente las primeras tres son estática mente estables.
Método Matricial de Nudos
Se ha demostrado que, para una armadura estáticamente determinada, el número de direcciones posibles en las cuales pueden actuar las fuerzas externas, JVP, debe ser igual al número de fuerzas axiales, NF. Sean P ( l ) a P(NP) las fuerzas externas conocidas en los nudos y F ( l ) a F(NF) las fuerzas axiales internas desconocidas. Como conexión de estos dos grupos de variables tenemos a las NP ecuaciones de estática, una en cada dirección de P, como sigue:
La matriz [A] como se define en la Ec. (3.3.2) se denomina la matriz estática debido a que sus renglones son las ecuaciones de estática para los nudos. Como la matriz [A] es cuadrada, tal como debe ser para una arma dura estáticamente determinada, la Ec. (3.3.2) puede escribirse como Esta última ecuación puede programarse en una computadora digital donde las matrices { P } y [ A ] , de fuerza externa j estática, respectivamente, fuesen los datos de entrada y la matriz de fuerzas internas {F} fuese la salida. El significado físico de la matriz [A- 1] consiste en que cada una de sus columnas dé las fuerzas axiales debidas a cada fuerza externa uni taria. Como ejercicio académico, la matriz [A-1] puede determinarse di rectamente del método convencional de nudos y secciones, y ser verificada confirmando que [A] [A- 1] =[I]. Si existe la inversa de la matriz estática [A], la armadura es estática mente estable, pero si la matriz estática [A] es singular, la armadura es estáticamente inestable. Armaduras Estáticamente Indeterminadas: Grados de Libertad;
Típicamente, una estructura reticular consiste de miembros conectados en nudos o uniones. Una armadura es una estructura reticular en la cual las cargas se aplican solamente en uniones perfectamente lisas, con pasadores, de barras rectas de sección transversal uniforme. En el capí tulo previo el símbolo NP se ha utilizado para indicar el número de di recciones posibles en las cuales las fuerzas pueden ser aplicadas en los nudos. Se hace notar aquí que NP es también el número de desplazamien tos nodales desconocidos; de hecho, hay una correspondencia biunívoca entre fuerzas nodales y desplazamientos. Como los desplazamientos nodales representan los modos en los cua les una estructura puede libremente responder a cualquier solicitación externa, el número de desplazamientos nodales desconocidos (o de di recciones posibles de fuerzas nodales) será llamado de aquí en adelante el grado de libertad NP. El Método de la Matriz de Desplazamientos En contraste con el método de la fuerza redundante, sólo relaciona do con estructuras estáticamente indeterminadas, el método de despla zamientos se aplica por igual a estructuras determinadas o indetermi nadas. La razón es que en este método, los desplazamientos nodales se toman como las incógnitas primarias, y hay siempre tantas ecuaciones de equilibrio en las direcciones de los grados de libertad como desplaza mientos nodales desconocidos, independientemente del grado de indeter minación de la estructura.
Es conveniente usar matrices para expresar las ecuaciones lineales involucradas en las aplicaciones del método de desplazamientos; de aquí el nombre de método de la matriz de desplazamientos. Se definirán dos matrices importantes: la matriz de fuerzas internas vs. la matriz de desplazamientos nodales [C] y la matriz de fuerza externa vs. la matriz de desplazamientos nodales [ K ], Considérese la armadura de la Fig. 4.9. Su grado de libertad ÍVP es cuatro, y las direcciones de fuerzas nodales posibles P1 a P t y los des plazamientos nodales correspondientes desconocidos X, a X 4 se numeran en sucesión en la Fig. 4.9a. El número de fuerzas en barras NF es cinco• a estas fuerzas en barras y las elongaciones correspondientes se les lla ma FJ a F 3 y e-¡, a e;, en la Fig. 4.9b.
La matriz de fuerzas internas vs. la matriz de desplazamientos nodales [C] expresa las fuerzas en las barras F( l ) a F(NF) en términos de los desplazamientos nodales X(l) a X(NP). Para la armadura de la
Fig. 4.9, en la cual Si = EA1/Ll a S5 = EA5/L-, son las fuerzas que alargan las barras respectivas en una unidad y entonces representan la rigidez de la barra. La matriz [C] puede establecerse por columnas mediante las fuerzas en las barras de las Figs. 4.9c, d, e y f. Nótese que solamente la componente del desplazamiento del nudo en la dirección de la barra causa un cambio en la longitud de ésta. La matriz de fuerza externa vs. la matriz de desplazamientos nodales [K] expresa las fuerzas externas P ( l - ) a P(NP) en términos de los desplazamientos nodales X ( l ) a X(JVP). Para la armadura de la Fig. 4.9,
EJERCICIOS 1. Resuelva la estructura mostrada. El material es de acero estructural con E = 2040 T/cm2 3
5T 3m
1
2
Áreas: A1-3 = 100 cm2 A3-2 = 150 cm2 A3-4 = 30 cm2 A1-4 = 40 cm2 A2-4 = 40 cm2
4 4m
4m
Luego de numerar los nudos, el siguiente paso es darle sentido alas barras.
3
1
2 4
ELEMENT O 1 ------ 3 1 ------ 4 3 ------ 2 4 ------ 2 4 ------ 3
Ø 36,87 0 -36,87 0 90
COSØ 0,8 1 0,8 1 0
SENØ 0,6 0 -0,6 0 1
(COSØ)^ 2 0,64 1 0,64 1 0
(SENØ)^ 2 0,36 0 0,36 0 1
COSØSEN Ø 0,48 0 -0,48 0 0
Utilizando la ecuación matricial:
c2 cs − c 2 − cs s 2 − cs − s 2 AE cs [K] = 2 L − c − cs c 2 cs 2 cs s 2 − cs − s u1
v1
u3
v3
u1 v1
u4 v4
K1−3
0.48 − 0.64 − 0.48 0.64 0.48 0.36 − 0.48 − 0.36 E = 0.48 5 − 0.64 − 0.48 0.64 0.36 − 0.48 − 0.36 0.48 U3
K 3− 2
u2
v2
0.64 − 0.48 − 0.64 0.48 0.48 − 0.36 3E − 0.48 0.36 = 0.64 − 0.48 10 − 0.64 0.48 0.48 − 0.36 − 0.48 0.36 u4 v4
K 4 −3
v3
0 0 E 0 1 = 10 0 0 0 − 1
K1− 4
u3 v3 0 0 0 − 1 0 0 0 1
Por superposición, la nueva matriz será:
1 E 0 = 10 − 1 0
0 −1 0 0 0 1 0 0
u4 v4
K 3− 2
1 E 0 = 10 − 1 0
0 0 0 0
u2 v2
0 −1 0 0 0 1 0 0
0 0 0 0
U1
K=
E 10
La ecuación matricial quedara así:
V1
U2
V2
U3 V3 U4 1,28 -0,96 -1
1 0,96 1,28
V4 U1
0,96 -0,72 1,44 1,92 1,44 -1
0,96 0,72
V1
1 1,92 -1,44 1,08 1,44 -1,08 1,28 -0,96 -1,92 1,44 1,28 1,92 0,96 -0,72 1,44 1,08 0,96 1,44 -1
U2 V2
0,96 -1,44
U3
1
-1
V3
0,72 1,08
-1
1 1 -1
U4 1
X1 0 0 − 1.28 − 0.96 − 1 0 2.28 0.96 Y 0.96 0.72 0 0 − 0.96 − 0.72 0 0 1 X2 = 0 0 0 2.92 − 1.44 − 1.92 1.44 − 1 0 0 − 1.44 1.08 1.44 − 1.08 0 0 Y2 = 10 0 X 3 = 4 E − 1.28 − 0.96 − 1.92 1.44 3.20 − 0.48 0 Y3 = 3 − 0.96 − 0.72 1.44 − 1.08 − 0.48 2.80 0 X = 0 −1 0 −1 0 0 0 2 4 0 0 0 0 −1 0 0 Y4 = −20
V4 u1 = 0 v = 0 1 u2 v2 = 0 0 u3 − 1 v3 0 u4 1 v4
Ahora se reordena y se llega a la nueva matriz con las incógnitas dadas: T
u2 2.92 − 1.92 1.44 − 1 0 0 u 3 10 − 1.92 3.20 − 0.48 0 0 4 v3 = 1.44 − 0.48 2.80 0 − 1 3 E 0 0 2 0 0 u4 −1 v 0 0 −1 0 1 − 20 4
cm.
1 − 1.333 1 − 1.333 0 0.1307 2 1 0.826 − 0.580 0.500 − 0.580 4 0.0645 10 = − 1.333 − 0.580 1.468 − 0.667 1.468 3 = − 0.1337 E 0.500 − 0.667 1 − 0.667 0 0.0654 1 − 1.333 − 0.580 1.468 − 0.667 2.468 − 20 − 0.2317 Ahora se procede a calcular las reacciones. 0.1307 − 1.28 − 0.96 − 1 0 0.0645 − 4 X1 0 Y = 204 0 − 0.96 − 0.72 0 0 − 0.1337 = 7.01 1 Y2 − 1.44 1.44 − 1.08 0 0 0.0654 10.01 T − 0.2317 Para finalizar se calculan las fuerzas internas. 0.0645 − 0 s1−3 = 408[ 0.8 0.6] = −11.68T (C ) − 0.1337 − 0 0.0654 − 0 s1− 4 = 204[1 0] = 13.34T (T − 0.2317 − 0 0.1307 − 0.0645 s3− 2 = 612[ 0.8 0.6] = −16.68T (T ) + 0.1337 0 0.1307 − 0.0645 s4 − 2 = 204[1 0] = 13.32T (T ) + 0.2317 0 0.0654 − 0.0654 s4 −3 = 204[ 0 1] = 20.00T (T ) − 0.1337 + 0.2317 2. Resuelva completamente la estructura mostrada
3 2m
Áreas: A1-2 = 20√2 cm2 A1-3 = 20 cm2 A3-2 = 20 cm2
2
1 4T 2m 5T
Y
3
2
1
x
Formula matricial: X1 u1 Y v 1 1 X 2 u 2 = [ K ] Y2 v2 X3 u3 Y3 v3 Cuadro para determinar las matrices de rigidez. ELEMENT O 1 -------- 2 1 -------- 3 2 -------- 3
Ø 180 135 90
COS -1 √2/2 0
SEN 0 √2/2 1
(COS)^2 (SEN)^2 COSSEN 1 0 0 0.5 0.5 -0.5 0 1 0
De la ecuación 2.0 se tiene: u 1 v1 u 2 v2
u1
v1
u3
v3
[ K ]1−2
1 E 0 = 10 − 1 0
0 −1 0 0 0 1 0 0
0 0 0 0
[ K ]1−3
0.5 − 0.5 − 0.5 0.5 0.5 − 0.5 E − 0.5 0.5 = 0.5 − 0.5 10 − 0.5 0.5 0.5 − 0.5 − 0.5 0.5
u 2 v2 u 3 v3
[ K ] 2−3
0 0 E 0 1 = 10 0 0 0 − 1
0 0 0 − 1 0 0 0 1
Ahora se ensambla la matriz por superposición. u1 1+ 0,5 K=
E 10
v1
u2 -0,5
-0,5
v2
u3
-1
0,5
-1
v3 -0,5
0,5 u1
0,5
-0,5 v1
1
u2 1
-0,5
0,5
0,5
-0,5
-1 v2 0,5
-1
-0,5 u3 1+ -0,5 0,5 v3
Ahora efectuando las operaciones y reemplazando, la expresión resultante será:
X1 1.5 − 0.5 − 1 0 − 0.5 0.5 u1 Y − 0 .5 0 .5 0 0 0.5 − 0.5 v1 1 X 2 E −1 0 1 0 0 0 u 2 = 0 0 1 0 − 1 v2 Y2 10 0 X3 − 0 .5 0 .5 0 0 0.5 − 0.5 u3 Y3 0.5 − 0.5 0 − 1 − 0.5 1.5 v3
X1 1.5 − 0.5 0 u1 − 1 − 0 .5 0 .5 u 2 = 0 Y = E − 0.5 0.5 0 v + E 0 0.5 − 0.5 u3 = 0 1 10 1 10 Y2 0 0 0 1 v2 0 1 v3 = 0 −1
u1 1.5 − 0.5 0 X 1 1 1 0 4 −10 E 10 10 v = − 0.5 0.5 0 Y = 1 3 0 − 5 = −110 1 E 1 E E v2 0 0 0 1 0 0 0 1 Y2
El siguiente paso es obtener las reacciones. 0 0 u1 0 0 u 2 = 0 X 2 −1 1 E X = E − 0 .5 0 .5 0 v1 + 0 0.5 − 0.5 v 2 = 0 3 10 10 Y3 0.5 − 0.5 − 1 v 2 0 − 0.5 1.5 v 3 = 0 0 0 −10 E 1T −1 E = − 0.5 0.5 0 −110 E = − 5T 10 0.5 − 0.5 − 1 0 5T
Las reacciones para la armadura son: ;
X 2 = 1T
X 3 = −5T
Y3 = 5T
Y como paso final se calculan las fuerzas internas en las barras. 0 + 10 E E [ − 1 0] = −1T (C ) 10 0 + 110 E 0 + 10 E E 2 ] = [− 2 2 2 = 5 2T (T ) 10 0 + 110 E
s1−2 = s1−3
s2 −3 =
0 − 0 E [ 0 1] =0 10 0 − 0
3. Calcular los desplazamientos en los nudos y las fuerzas en las barras de esta cercha. E y A son constantes para todas. Ecuación completa:
P1 P 2 = P3 P 4
K11
K12 K 22
K 13 K 23 K 33
0 ∆ 1 K 24 ∆ 2 K 34 ∆ 3 K 44 ∆ 4
Condiciones de contorno: ∆1 = 0 i = 1, 4
Ecuación final:
P2 K 22 P = K 3 32 Rigidez del sistema:
–
Elemento 2,1
K 23 ∆ 2 K33 ∆ 3
0−0 =0 15 15 − 0 m= =1 15 EA 0 0 K 21 = − 15 0 1
l=
–
Elemento 2,3 l = −1 m=0 K 23 = −
–
EA 1 0 20 0 0
Elemento 2,4 4 5 3 m= 5
l=−
K 24 = −
−12 16 EA 25 25 9 25 −12 25 25
0.0753 −0.0192 − [ K 21 + K 23 + K 24 ] = K 22 = EA −0.0192 0.081
–
Elemento 3,2 l =1 m=0 K 32 = −
–
Elemento 3,1
EA 1 0 20 0 0
4 5 3 m= 5
l=
12 16 EA 25 25 K 31 = − 9 25 12 25 25
–
Elemento 3,4 l=0 m =1 K 34 = −
EA 0 0 15 0 1
0.0753 0.0192 − [ K 32 + K 31 + K 34 ] = K 33 = EA 0.0192 0.081 –
Ecuación final:
0 ∆ 2 X 8 0.0753 − 0.0192 − 0.05 0 − 0.0192 0.081 0 0 ∆ 2Y = EA 0 − 0.05 0 0.0753 0.0192 ∆ 3 X 0 0.0192 0.081 ∆ 3Y −5 0 –
Desplazamientos en los nudos:
∆ 2 X 246.3 ∆ 2Y = 1 58.4 ∆ 3 X EA 190.0 −106.8 ∆ 3Y –
Fuerzas en los elementos: Pij =
a. Tensión
EA RIJ (∆i − ∆ j ) L
P12 =
−246.3 1 EA = 3,89 [ 0 − 1] 15 −58.4 EA
b. Compresión P23 = –
246.3 190.0 1 EA = −2.82 [ −1 0] 20 58.4 106.8 EA
Análogamente P13 = 3.52 P24 = −6.48 P34 = −7.12
–
Cálculo de reacciones a partir de las fuerzas en los elementos:
−4 0 ¨− 2.84k 5 P1 = P12 = P13 = R13P12 + R13P13 = 0.389 + 3.52 = −3 −1 −6.02k 5 4 0 − 5.18K 5 P4 = − 7.12 + − 6.48 = −3 −1 11.00K 5 –
Equilibrio de la estructura: ∑ PX = 8 − 2.84 − 5.18 ≈ 0 ∑ PY = −5 − 6.02 + 11.0 ≈ 0
4. Calcular las fuerzas en todos los elementos de esta cercha debidas a la carga mostrada. E, A, son constantes en todas partes. Las condiciones de contorno:
La ecuación final:
P1 = K11 ∆1 Donde
K11 = ∑ 1
EA B L
l 2 l m B= 2 l m m Xi − X j l= Lij m=
Rigideces de los elementos: Elemento 1,2:
Yi − Y j Lij
3−0 =1 3 0−0 m= =0 3 EA 41.66 0 K12 = − 0 125 0
l=
Elemento 1,3: l=
3 5
m=−
4 5
K13 = −
EA 9 −12 125 −12 16
3 5 4 m=− 5
l =−
Elemento 1,4:
K14 = −
Elemento 1,5:
Elemento 1,6 3 l =− 5 4 m= 5
:
Elemento 1,7:
EA 9 12 125 12 16
−12 EA 9 K15 = − 125 −12 16
l =0 m =1 K15 = −
EA 0 0 125 0 12.5
3 5 4 m= 5
l=
K15 = −
EA 9 12 125 12 16
K11 = − [ K12 + K13 + K14 + K 15 + K 16 + K 17] =
0 EA 77.66 76.5 125 0
La ecuación matricial de rigidez final es:
0 ∆1x 100 EA 77.66 200 = 125 0 76.5 ∆1 y Los desplazamientos de los nudos libres son: 1 ∆1x 125 77.66 ∆ = 1 y EA 0
0 100 1 160.9 = 1 200 EA 326.8 76.5
Las fuerzas en los elementos son: EA Rij ( ∆ i − ∆ j ) L 160.9 1 EA k P12 = [ 1 0] EA = 53.6 ( tensión ) 326.8 3 Pij =
P13 =
Análogamente
EA 3 5 5
4 160.9 1 − = −33k ( compresión ) 5 326.8 EA
P14 = −71.6k ( compresión ) P15 = 33k ( tensión ) P16 = 32.7 k ( tensión ) P17 = 71.6 k ( tensión ) 5. Una armadura tridimensional se define por los siguientes valores: índices de los elementos, coordenadas de los nudos y condiciones de contorno. La cercha está sometida a tres condiciones de carga diferentes. Determinar las fuerzas en los elementos para cada condición de carga y verificar que los resultados de la tercera condición de carga son iguales a la suma de los dos primeros. E e 1 son constantes para todos los elementos.
Después de que introducen las condiciones de contorno, la ecuación final toma la siguiente forma:
Donde
Elemento 1,2:
Elemento 1,3:
Elemento 1,4:
Elemento 1,5:
Elemento 1,6:
Los desplazamientos del nudo 1 (el único nudo libre en esta estructura) Son:
Como puede verse de los anteriores resultados,
Las fuerzas en los elementos son:
Elemento 1,2:
Todas las otras fuerzas en las barras pueden calcularse análogamente.
Se comprueba el equilibrio del nudo 1
Que es igual a las cargas aplicadas en el nudo 1.
6. Se pide resolver completamente la estructura mostrada en la figura para cualquier hipótesis de carga, y en particular para las cargas mostradas. Todos los elementos están hechos del mismo material.
Solución Se comienza por numerar los nudos y asignar a los elementos un sentido positivo arbitrario, pero sobre teniendo que el sentido positivo va del nudo i al nudo j, esto se ha hecho en la figura:
.
Obsérvese que el nudo libre se ha dejado de primero, seguido del apoyo inferior que ofrece una restricción y de último se ha dejado el apoyo superior que esta restringido en ambos sentidos
La formulación matricial del problema será así x1
=
k
*
u1
y1 x2 y2 x3 y3
v1 u2 v2 u3 v3
Para determinar fácilmente las matrices de rigidez de los elementos, conviene elaborar el siguiente cuadro elemento 1-2 1-3 2-3
ө 180º 135º 90º
cos ө -1 -√2/2 0
sen ө 0 √2/2 1
cos² ө 1 ½ 0
sen² ө 0 ½ 1
cos ө sen ө 0 -½ 0
Y de ahí:
K(1-2)
K(1-3)
K(2-3)
=
=
=
E/10
U1 1 0 -1 0
V1 0 0 0 0
U2 -1 0 1 0
V2 0 0 0 0
U1 V1 U2 V2
E/10
U1 1/2 -1/2 -1/2 1/2
V1 -1/2 1/2 1/2 -1/2
U3 -1/2 1/2 1/2 -1/2
V3 1/2 -1/2 -1/2 1/2
U1 V1 U3 V3
E/10
U2 0 0 0 0
V2 0 1 0 -1
U3 0 0 0 0
V3 0 -1 0 1
U1 V1 U3 V3
Ensamblando ahora la matriz por superposición K = K(1-2) + K(1-3) + K(2-3) Se expande cada una de las matrices y luego se suman o simplemente se van sumando los términos respectivos en las casillas correspondientes Nótese que en este caso particular, AE/L resulta constante para todos los elementos lo cual facilita la superposición. Efectuando las operaciones indiciadas y reemplazando la expresión resultante queda:
= X1 Y1 X2 Y2
E/10
U1 3/2 -1/2 -1 0
V1 - 1/2 ½ 0 0
U2 -1 0 1 0
V2 0 0 0 1
U3 - 1/2 ½ 0 0
V3 ½ -1/2 0 -1
U1 V1 U2 V2
X3 Y3
-1/2 1/2
½ - 1/2
0 0
0 -1
½ -1/2
-1/2 3/2
U3 V3
Considerando ahora las condiciones de los apoyos, se procede a reordenar la expresión anterior intercambiando la tercera fila con las cuarta y por consiguiente las columnas respectivas, resulta entonces
X1 Y1 Y2 X2 X3 Y3
=
E/10
U1 3/2 -1/2 0 -1 -1/2 1/2
V1 - 1/2 1/2 0 0 ½ - 1/2
U2 -1 0 1 0 0 0
V2 0 0 0 1 0 -1
U3 - 1/2 ½ 0 0 ½ -1/2
V3 1/2 -1/2 0 0 -1/2 3/2
U1 V1 V2 U2 U3 V3
La ecuación anterior sirve de base para resolver el caso general. Para el caso particular en consideración procediendo directamente:
X1 Y1 Y2
=
U1 V1 V2
E/10
=
3/2 -1/2 0
-1/2 1/2 0
0 0 1
U1 V1 V2
-1 1 0
10/E
+ 1 3 0
E/10
-1 0 0
0 0 1
-1/2 1/2 0 4 -5 0
1/2 -1/2 1 =
U2=0 U3=0 V3=0 -10/E -110/E 0
Para obtener las reacciones a partir de la primera ecuación X1 X3 Y3 =
=
E/10
10/E
-1 -1/2 1/2 -1 -1/2 1/2
0 1/2 -1/2
0 0 -1
10 1/2 -1/2
U1 V1 V2
+ 0 0 -1
E/10 -10/E -110/E 0
-1 0 0
0 1/2 -1/2 =
0 -1/2 3/2
U2=0 V2=0 V3=0
1T→ -5T← 5T↑
Colocando estos valores en un diagrama de cuerpo libre de la estructura total se comprueba que la estructura esta en equilibrio.
7. Resuelva completamente la estructura mostrada. El material es acero estructural con E= 2040 T/cm². Las áreas están dadas entre paréntesis en cm².
Solución Se empieza por numerar los nudos y asignarle un sentido a las barras (figura superior derecha). Obsérvese que se han numerado de últimos los nudos libres para minimizar el reordenamiento de la matriz de rigidez. A continuación se elabora el cuadro de funciones trigonométricas elemento 1-3 1-4 3-2 4-2 4-3
ө 69.87º 0º -36.87º 0º 90 º
cos ө 0.8 1 0.8 1 0
sen ө 0.6 0 -0.6 0 1
cos² ө 0.64 1 0.64 1 0
sen² ө 0.36 0 0.36 0 1
cos ө sen ө 0.48 0 -0.48 0 0
Utilizando este cuadro y las ecuaciones:
K(1-3)
K(1-4)
=
=
E/5
U1 0.64 0.48 -0.64 -0.48
V1 0.48 0.36 -0.48 -0.36
U3 -0.64 -0.48 0.64 0.48
V3 -0.48 -0.36 0.48 0.36
U1 V1 U3 V3
E/10
U1 1 0 -1 0
V1 0 0 0 0
U4 -1 0 1 0
V4 0 0 0 0
U1 V1 U4 V4
V3 -0.48 0.36 0.48 -0.36
U2 -0.64 0.48 0.34 -0.48
V2 0.48 -0.36 -0.48 0.36
U3 V3 U2 V2
V4
U2
V2
K(3-2)
=
3E/10
U3 0.64 -0.48 -0.64 0.48
K(4-2)
=
E/100
U4
K(4-3)
=
E/10
1 0 -4 0
0 0 0 0
-1 0 1 0
0 0 0 0
U4 V4 U2 V2
U4 0 0 0 0
V4 0 1 0 -1
U3 0 0 0 0
V3 0 -1 0 1
U4 V4 U3 V3
Ensamblando por superposición, la ecuación matricial queda de esta manera X1 Y1 Y2 X2=0 X3=4 Y3=3 X4=0 Y4=-20
2.28 0.96 0 0 =E/10 -1.28 -0.96 -1 0
0.96 0.72 0 0 -0.96 -0.72 0 0
0 0 1.08 -1.44 1.44 -1.08 0 0
0 0 -1.44 2.92 -1.92 1.44 -1 0
-1.28 -0.96 1.44 -1.92 3.20 -0.48 0 0
-0.96 -0.72 -1.08 1.44 -0.48 2.80 0 -1
-1 0 0 -1 0 0 2 0
0 0 0 0 0 -1 0 1
U1=0 V1=0 V2=0 U2 U3 V3 U4 V4
Efectuando la partición en la forma establecida e invirtiendo
=10/E
2 1 -1.333 1 -1.333
1 0.826 -0.580 0.500 -0.580
-1.333 -0.580 1.468 -0.667 1.468
1 0.500 -0.667 1 -0.667
-1.333 -0.580 1.468 -0.667 2.468
-0.96 -0.72 -1.08
-1 0 0
0 4 3 0 -20
0.1307 0.0645 -0.1337 0.654 -0.2317
=
Ahora se pueden calcular las reacciones X1 Y1 Y2
=
204
0 0 -1.44
-1.28 -0.96 1.44
0.1307 0.0645 -0.1337 0.0654 -0.2317
0 0 0
=
-4 7.01 10.01
Finalmente se calculan las fuerzas internas 0.0645 S(1-3) S(1-4)
= =
408
0.8
0.6
204 1
-0
0
-0.1337
-0
0.0654
-0
=
11.68T
=
13.34T
S(3-2)
S(4-2)
S(4-3)
0
=
=
612
204
204
0.8
1
0
-0.2317
-0
0.1307
-0.0645
-0.6 0
0.1337
0.1307
-0.0654
0 0
0.2317
0.0645
-0.0654
1 -0.1337
0.2317
Los siguientes diagramas prueban la bondad de estas respuestas:
8. Resuelva la estructura que se muestra a continuación
Cuadro de áreas
=
16.68T
=
13.32T
=
20.00T
Barra 1-2 1-3 1-4 1-6 2-3 2-4 2-5 3-5 3-6 4-5 4-6 5-6
Área cm² 20 20 40 50 20 50 40 50 40 10 10 10
Cargas aplicadas en toneladas nudo 1 2 3
Fx 10.00 5.00 -4.00
Fy 15.00 -3.00 -2.00
Fz -12.00 -10.00 -6.00
y 6000 6000 6000 0 0 0
Z 4800 2400 4800 6000 0 6000
Coordenadas de los nudos nudo 1 2 3 4 5 6
x 2250 3750 5250 0 3750 7500
Se evalúan las matrices de rigidez de todos los miembros, referidas a coordenadas generales, cuidando de calcular únicamente los términos de las columnas indispensables, es decir, las que corresponden a los nudos libres
K(1-2)=
U1 1.169 0 -6.671 -4.169 0 6.671
V1 0 0 0 0 0 0
W1 -6.671 0 10.674 6.671 0 -10.674
U2 -4.169 0 6.671 4.169 0 -6.671
V2 0 0 0 0 0 0
W2 6.671 0 -10.674 -6.671 0 10.674
U1 V1 W1 U2 V2 W2
K(1-3)=
U1 14.000 0 0 -14.000 0 0
V1 0 0 0 0 0 0
W1 0 0 0 0 0 0
U3 -14.000 0 0 14.000 0 0
V3 0 0 0 0 0 0
W3 0 0 0 0 0 0
U1 V1 W1 U3 V3 W3
K(2-3)=
U2 4.169 0 6.671 -4.169 0 -6.671
V2 0 0 0 0 0 0
W12 6.671 0 10.674 -6.671 0 -10.674
U3 -4.169 0 -6.674 4.169 0 6.671
V3 0 0 0 0 0 0
W3 -6.671 0 -10.674 6.671 0 10.674
U2 V2 W2 U3 V3 W3
K(1-6)
U1 5.523 -6.312 1.262 -5.523 6.312 -1.262
V1 -6.312 7.214 -1.443 6.312 -1.214 1.443
W1 1.262 -1.443 0.288 -1.262 1.443 -0.288
U1 V1 W1 U6 V6 W6
K(2-4)
U2 2.951 4.721 -2.833 -2.951 -4.721 2.833
V2 4.721 7.554 -4.533 -4.721 -7.554 4.533
W2 -2.833 -4.533 2.720 2.830 4.533 -2.720
U2 V2 W2 U4 V4 W4
K(3-5)
U3 0.492 1.969 1.575 -0.492 -1.969 -1.575
V3 1.969 7.877 -6.302 -1.969 -7.877 -6.302
W3 1.575 6.302 -1.575 -1.262 -6.302 -5.041
U3 V3 W3 U5 V5 W5
U1 1.535 4.093
V1 4.093 10.915
W1 -0.819 -2.183
U1 V1
K(1-4)
-0.819 -1.535 -4.093 0.819
-2.183 -4.093 -10.915 2.183
0.437 0.819 2.183 -0.437
W1 U4 V4 W4
K(3-6)
U3 1.535 -4.093 0.819 -1.535 4.093 -0.819
V3 -4.093 10.915 -2.183 4.093 -10.915 2.183
W3 0.819 -2.183 0.437 -0.819 2.183 -0.437
U3 V3 W3 U6 V6 W6
K(2-5)
U2 0 0 0 0 0 0
V2 0 11.206 4.482 0 -11.206 -4.482
W2 0 4.482 1.793 0 -4.482 -1.793
U2 V2 W2 U5 V5 W5
El siguiente paso es el ensamble de la matriz de rigidez de la estructura Como los tres apoyos 4,5 y 6 son de segundo genero, no tienen ningún desplazamiento y basta ensamblar las porciones [Knn] y [Kan] de la matriz de rigidez total. Para luego aplicar las ecuaciones simplificadas
Y resolviendo para los desplazamientos desconocidos mediante la inversión de [Knn]
Conociendo los desplazamientos, de la ecuación anterior se pueden averiguar las reacciones
Conviene ahora verificar el equilibrio de la estructura, conociendo el cuerpo libre total, se obtiene
Cuya magnitud es Mo =√ (18²+0+33²)=37.6 T-mm < 0.04 T-mm Comprobándose así que helero de cierre en la estructura total es insignificante. Luego se calculan las fuerzas internas que debe soportar cada barra.
S12=
14.841
0.5300
0
-0.8481
2.2263 -0.7276 -2.7318
-0.8043 -0.0331 4.4634
=-10.609T
S13=
S23=
S45=
14.000
14.841
2.968
1.000
0
0.5300
0
0.5300
0
0.7508 0.3669 -1.7722
-0.8043 -0.0331 4.4634
=-0.750T
0.8481
0.7508 0.3669 -1.7722
-2.2263 0.7276 2.7318
=0.471T
0 0 0
=0T
0
-0.8481
S46= S56=
S16=
S24=
S35=
2.800
1.000
0
0
2.968
0.5300
0
0.8481
13.204
13.226
13.412
0.6512
-0.4724
-0.1916
0.744 2
0.755 8
0.766 4
0.1488
0.4534
-0.6131
S14= 12.885
-0.3451
0.920 4
0.1841
0 0 00 0 0
0
-0.8043
0
-0.0331
0
4.4634
0
-2.2263
0
0.7276
0
2.7318
0
-0.7508
0
-0.3669
0
1.7722
0 0
-0.8043 -0.0331
=0T =0T
=2.152T
=23.019T
=-8.872T
=14.557T
S36=
S25=
12.885
12.999
0.3451
0
0.920 4
0.928 5
0.1841
-0.3714
0
4.4634
0
-0.7508
0
-0.3669
0
1.7722
0
-2.2263
0
0.7276
0
2.7318
=5.216T
=-21.970T
Finalmente se verifica el equilibrio de cada nudo. miembro 1-2 1-3 1-4 1-6 2-3 2-4 2-5 3-5 3-6 4-5 4-6 5-6
Sx Sy Sz -5.623* 0 8.997 -0.750 0 0 -5.024 -13.398 2.680 1.401 -1.602 0.320 0.250 0 0.400 -10.874 -17.398 10.437 0 20.399 8.160 1.700 6.800 5.439 1.800 -4.801 0.960 0 0 0 0 0 0 0 0 0 *nudo inicial , **C= compresión , T= tensión
esf C** C T T T T C CV T
Con estos valores y los de las cargas aplicadas se comprobó el equilibrio de los nudos superiores, obteniéndose: nudo ∑Fx ∑Fy ∑Fz 1 0.004 0.000 -0.003 2 -0.001 0.001 0.000 3 0.000 -0.001 -0.001 9.
Como los desplazamientos en los apoyos son nulos, basta con considerar la porción siguiente X1 Y2 X2 Y2 X3=4 Y3=3 X4=0 Y4=-20
-1.28 -0.96 -1.92 1.44 3.20 -0.48 0 0
=E/10
-0.96 -0.72 1.44 -1.08 -0.48 2.80 0 -1
-1 0 -1 0 0 0 2 0
0 0 0 0 0 -1 0 1
U3 V3 U4 V4
Tomando la parte inferior y reordenado para invertir por partes X3 Y3 Y4 X4
= = = =
4 3 -20 0
[0] = E/10 [ 0 0 0 ] +
=E/10
3.20 -0.48 0 0
-0.48 2.80 -1 0
0 -1 1 0
0 0 0 2
U3 V3 V4 U4
U3 V3 2E/10 [U4] ; U4 = 0 V4
Por consiguiente: 4 3 -20 Y despejando: U3 V3 V4 =
=E/10
3.20 -0.48 0
-0.48 2.80 -1
0 -1 1
U3 V3 V4
=10/E
0.3266 0.0868 0.0868
0.0868 0.5788 0.5788
0.0868 0.5788 1.5788
4 3 -20
-8.49 x 10 ^-4 -0.0465 -0.1446
Reemplazando estos valores X1 Y1 X2 Y2
=204
-1.28 -0.96 -1.92 1.44
-0.96 -0.72 1.44 -1.08
-1 0 -1 0
0 0 0 0
-0.0008 -0.0465 0 -0.1446
=
9.33 7.00 -13.34 10.00
Haciendo el diagrama de cuerpo libre de la estructura total:
En cuanto a fuerzas interna -0.0008 S(1-3)
S(1-4)
S(3-2)
S(4-2)
-3)
=
=
=
=
=
408
204
612
204
204
0.8
1
0.8
1
0
-0
0.6 -0.0465
-0
0
-0
0 -0.1446
-0
0
0.0008
-0.6 0
0.0465
0
-0
0
0.1446
-0.0008
-0
-0.0465
0.1446
0
1
=
-11.64T
=
0
=
-16.68T
=
0
=
20.01T
Verificándose el equilibrio de todos los nudos con un error de cierre máximo de 0.03T
10. Una armadura de dos elementos
Parámetros: Número de Nodes N = 3, Número de Elementos M = 2, propiedades de los elementos: EA = L =1 para ambos elementos
Elemento numero e1:
Cartografía de Matrix: nodo local de un nodo = global (2), el nodo local b = nodo global (1). Por lo tanto,
Rigidez Ampliado elemento de la matriz:
Elemento numero e2:
Cartografía de Matrix: nodo local de un nodo = global (3), nodo local b = nodo global (1). Por lo tanto,
Rigidez Ampliado elemento de la matriz:
Rigidez total matriz:
Tenga en cuenta las condiciones de carga: PX1 = 0 y AP1 = -30. Mundial de nodos (2) y (3) son fijos. La ecuación de equilibrio por lo tanto:
Desde UX1 y Uy1 son los desplazamientos única incógnita, la ecuación matricial anterior se puede simplificar, descartando los elementos de matriz que se multiplican por el desplazamiento cero en los nudos fijos, es decir,
Resolver la ecuación ahora la parte superior de la matriz particionada de la siguiente manera:
a continuación, calcular la reacción de las fuerzas de la partición utilizando como fondo
Una nota interesante que resulte de este sencillo ejemplo numérico es que la matriz que se utilizará para la determinación de los desplazamientos desconocidos se puede obtener mediante la eliminación de todas las columnas y filas asignado a la fija grados nodales de libertad, desde el mundial de la matriz. Inestabilidad de las condiciones de Fronteras.
Consideremos un caso similar ahora con la idéntica condición de carga, pero nodo global (3) está en un rodillo en la dirección x, por lo tanto, UX2 = Uy2 Uy3 = = 0, pero Ux3 se deja como un desconocida. Sabemos que esto es físicamente un sistema inestable, pero permite comprobar numéricamente. Al eliminar las filas y columnas asociadas con grados de libertad UX2, Uy2 y Uy3, la ecuación de equilibrio rendimientos
Se puede demostrar fácilmente que la ecuación matriz anterior es singular, calculando el valor de ∆ como factor determinante.
11. 1) La cercha simple tiene m= 9 barras, se desea calcular los desplazamientos verticales de los nudos del cordón superior a y b, bajo la acción de dos condiciones de carga dependientes, la numeración de las barras y las dimensiones están a continuación:
Cada una de las barras tiene una sección transversal Ai= 1 cm2 y un modulo de elasticidad E= 30000 Kg/cm2. Se comienza con un cálculo de los esfuerzos axiles Si en las nueve barras bajo las dos condiciones de carga. Este conjunto de esfuerzos axiales constituirá la matriz S para el problema. Igualmente se calculan los esfuerzos axiales en todas las barras para cada una de las cargas unidad que actúan en los nudos a y b. Este conjunto de esfuerzos axiales constituirá la matriz s para el problema en cuestión. Finalmente, se calculan los coeficientes de flexibilidad para cada elemento y se escribe la matriz de flexibilidad diagonal. Estas tres matrices básicas para el caso que nos ocupa son las siguientes: Sij
Pi =10E
Si
Ahora se halla los desplazamientos requeridos Δj por multiplicación matricial del siguiente modo:
Efectuando en dos pasos las multiplicaciones indicadas, se obtiene:
Así con E= 30000 Kg/cm2, los desplazamientos debidos a la carga P= 9Kg que actúa según se indico: Δ´a= 0,0318 cm Δ´b= 0,0525 cm Mientras que las debidas a la carga Q=10 kg son: Δ´´a= 0,0237 cm Δ´´b= 0,0296 cm 12. Considerando la cercha estáticamente indeterminada. Esta estructura plana tiene siete barras, de las que tres son superabundantes. Se elige como tales las barras 5, 6 y 7, de modo que el sistema principal sea:
Como primer paso se hace un análisis estático de este sistema principal bajo la acción de las cargas P dadas y, asimismo, bajo la acción de un par de fuerzas unidad en lugar de cada una de las barras superabundantes. Después, con las secciones transversales A1=A4=A5= 5 cm2, A2=A3=8,66 cm2 y A6 = A7=2,89cm2, la matriz geométrica, la matriz de flexibilidad y la matriz de carga serán:
Donde, por conveniencia se ha tomado E= 1 Kg/cm2 Transponiendo ahora s y s´ y postmultiplicándola por f y después por s, se obtiene:
Como este resultado es una matriz diagonal, se convierte fácilmente y se tiene:
Premultiplicando ésta por la matriz s, se obtiene:
La postmultiplicación de ésta por s’f da:
Finalmente, restando ésta de la matriz unidad I y completando las operaciones indicadas:
Así, los valores finales, en unidades kilogramos, de las fuerzas axiales en las barras son:
13.
Las fuerzas de as reacciones se pueden calcular desde esta nueva solución:
14) Determine el desplazamiento horizontal y vertical en la junta (3) de el ensamble.
Use la rigidez hallada en la matriz del ensamble anterior y deje Q=KD
Solucionando
15) Determine la fuerza guiándose por el problema 1