Servicios Personalizados
Revista
Articulo
Indicadores
-
Citado por SciELO
-
Accesos
Links relacionados
-
Similares en SciELO
Compartir
Boletín Técnico
versión impresa ISSN 0376-723X
IMME v.45 n.3 Caracas nov. 2007
Método de holzer - van den dugen en el análisis de pórticos.
Jorge Cruz
Universidad Central de Venezuela, Facultad de Ingeniería, Escuela de Ingeniería Mecánica. Departamento de Diseño. e-mail: jorge_cruz@cantv.net.
Resumen
Este trabajo muestra la aplicación del método de Holzer Van den Dungen en su forma matricial en el cálculo de pórticos planos. Se presenta el desarrollo del método, las matrices de transferencia para cargas longitudinales y transversales en elementos de ejes rectos; las matrices de rotación, las ecuaciones de compatibilidad y equilibrio de nodos a los cuales concurren elementos con orientaciones diferentes.
Palabras Claves: Matriz de transferencia, cargas longitudinales y transversales, matriz de rotación.
Holzer-van den dugen method for plane frames.
ABSTRACT
This paper shows the application of the method of Holzer - Van den Dungen in matrix mode to calculate plane frames. The development of the method, the longitudinal and transversal loads transfer matrix in elements of straight axis appears, with the rotation matrix, the equations of compatibility and equilibrium of nodes to which elements with different directions concur.
Keywords: Longitudinal and transversal loads transfer matrix loads, rotation matrix.
Recibido: 26/05/06 Revisado: 27/03/07 Aceptado: 20/08/07
1. IntroducCión
Según Dimarogonas (1996), Holzer, Güembel y Tolle desarrollaron entre 1907 y 1921 el método conocido como método de Holzer para el cálculo de de frecuencias naturales torsionales de ejes y en 1928 Van den Dungen lo amplió para vibraciones laterales de barras. Este método toma en cuenta el hecho de que una gran cantidad de problemas de ingeniería, y en especial aquéllos que se modelan por líneas, el comportamiento de cada sección transversal está influenciado solamente por el comportamiento de las secciones vecinas. Esta forma de modelar el problema permite obtener las respuestas estáticas y dinámicas sin necesidad de establecer las ecuaciones que lo gobiernan. Según Dimarogonas (1996), Thomson desarrolló el método en forma matricial y lo denominó método de la matriz de transferencia.
En este paper se muestra cómo el método de la matriz de transferencia puede utilizarse en la resolución de pórticos estructurales.
2. LIMITACIONES
Las limitaciones del desarrollo que se presenta a continuación, son las habituales en los modelos lineales de la mecánica de sólidos:
-
Las secciones transversales son planas y permanecen planas después de la deformación.
-
El material es lineal y elástico, y se considera que se mantienen estas propiedades durante toda la historia de carga del pórtico.
-
Uno de los ejes centroidales principales de inercia de cada una de las secciones transversales está contenido en el plano del pórtico.
-
Las deformaciones son pequeñas.
-
Se desprecian las deformaciones por corte.
-
El pórtico se considera libre de esfuerzos residuales y de imperfecciones.
-
Las juntas son rígidas.
3. MATRIZ DE TRANSFERENCIA EN PÓRTICOS
Los pórticos estructurales se modelan mediante líneas, donde los valores de las variables de estado (fuerza axial, fuerza cortante, momento flector, desplazamientos y rotación) en una sección transversal solo dependen de lo que ocurre en las secciones transversales vecinas. Este tipo de problemas puede estudiarse mediante matrices de transferencia, las cuales dependerán de la orientación asociada a cada línea y de las cargas en ellas contenidas.
Supóngase el pórtico bidimensional mostrado en la figura 1. En éste existen cambios de orientación y cargas aplicadas. El pórtico puede dividirse en elementos, los cuales tienen la característica de ser prismáticos, de eje recto y la función de carga que los solicita, si la hay, debe ser continua. Los extremos de los elementos están conectados a los nodos, que representan los planos donde comienzan y terminan los elementos. Lo anterior indica que deben colocarse nodos donde existan cambios en la geometría de la sección transversal, cambios de material, presencia de cargas concentradas (axial, transversal o momento flector), cambios en la función de carga externa y/o cambios en la orientación de los elementos.
Figura 1. Pórtico.
Para describir completamente la situación de cada sección, se requiere conocer la magnitud de seis cantidades, denominadas variables de estado que son: el desplazamiento d en la dirección axial de la sección, el desplazamiento y en la dirección normal, la rotación q que gira la sección transversal respecto a su orientación original, la fuerza axial A, el momento flector M, y la fuerza cortante V.
La figura 2 muestra un elemento i del sistema entre los nodos i e i+1 en su condición deformada con la convención de signos habitual. Los desplazamientos axiales di y di+1 medidos desde las posiciones iniciales de las respectivas secciones se han omitido por claridad del dibujo. Las variables de estado forman el vector zi = [d y q A M V ] que describe el estado del sistema en el nodo i, y se denomina vector de estado. Nótese que estas variables de estado están asociadas con la orientación indeformada del elemento. Esto conduce a la necesidad de establecer dos sistemas de referencia, el primero, que en la literatura se llama global y respecto al cual se establecen las coordenadas de los nodos y la orientación de los elementos, y el segundo denominado local, que es solidario con cada elemento y en el que uno de los ejes, generalmente el eje x, se hace coincidir con el eje axial del elemento.
Figura 2. Elemento de la estructura en condición deformada.
Supóngase que en el nodo i el vector de estado es:
Zi = [di yi qi Ai Mi Vi ] (1)
de valores desconocidos. Si entre los nodos i e i+1 hay fuerzas distribuidas longitudinales representadas por una función fi(x), un cambio de temperatura uniforme DTi y fuerzas distribuidas transversales representadas por una función continua y creciente de grado n de la forma pxi = pi xn / li n (siendo pi la intensidad máxima de la carga sobre el elemento) que soliciten al elemento entre los nodos i e i+1, por la teoría simple de viga-columna, las variables de estado del nodo i+1 en función de las del nodo i serán: (las cargas longitudinales y transversales no se muestran en la figura 2)
(2)
Las tres últimas ecuaciones expresan el equilibrio de fuerzas y momentos, y las tres primeras, el desplazamiento longitudinal, el desplazamiento transversal y la rotación de la curva elástica debido a estos momentos y fuerzas.
Este conjunto de ecuaciones puede escribirse en forma matricial así:
(3)
con
(4)
donde el subíndice i indica las propiedades l, E, A, I del elemento i. La fila 7 se ha agregado por conveniencia computacional. Es importante notar que la matriz de transferencia Li fue obtenida conforme a un sistema de referencia local, es decir, uno de los ejes es el eje axial del elemento y el otro es perpendicular a dicho eje.
La ecuación (3) indica que el vector de estado del nodo i+1 es igual al vector de estado del nodo i premultiplicado por la matriz cuadrada Li de 7*7, la cual depende solo de las propiedades del elemento que generalmente son conocidas y de las cargas que lo solicitan. Esta matriz Li se denomina matriz de transferencia. Por cada elemento viga-columna existe una matriz de transferencia.
Si lo elementos i+1 e i+2 tienen la misma orientación, la repetición del proceso previo conduce a las ecuaciones siguientes:
(5)
(6)
(7)
En el caso de los nodos donde concurren dos elementos con orientaciones diferentes, como el mostrado en la figura 3, se requiere una matriz de transferencia mediante la cual se relacione el vector de estado correspondiente al nodo, respecto a los sistemas de referencia locales asociados a las direcciones de dichos elementos.
Figura 3. Elementos no colineales.
Para este fin se presenta la figura 4. En ella se tienen dos sistemas de referencia denominados XY y XY ambos con un origen de coordenadas común en O, que se considera como la posición de un punto de un sólido deformable antes de la deformación. El punto luego de la deformación pasa a la posición b. Las componentes del desplazamiento del punto respecto a los dos sistemas de referencia se relacionan mediante las ecuaciones.
Figura 4. Componentes del desplazamiento de un punto respecto a dos sistemas de coordenadas.
(8)
(9)
Que expresadas en forma matricial conducen a la conocida matriz de rotación expresada por la ecuación (10).
(10)
Figura 5. Componentes del desplazamiento de un nodo respecto a dos sistemas de referencia orientados en las direcciones longitudinales de los elementos que concurren a dicho nodo.
Como el desplazamiento del nodo es único, la transformación anterior permite relacionar las componentes de dicho desplazamiento del nodo respecto a los dos sistemas de referencia locales, cada uno de ellos asociado a la orientación de cada elemento, ambos con origen en el nodo, y que se diferencian exclusivamente en la rotación a, como se muestra en la Figura 5. En este caso, estas componentes del desplazamiento son los desplazamientos axiales y transversales de las secciones de los elementos que concurren al nodo. Pudiendo expresarse en forma matricial así:
(11)
Por la consideración de junta rígida se tiene
(12)
ampliándose la transformación anterior a
(13)
que en forma compacta se expresa así
(14)
con
(15)
Donde zdR y zdL son las componentes del desplazamiento y la rotación del nodo según la orientación de los sistemas de referencia locales asociados con cada uno de los elementos que concurren al nodo. Rd es la matriz de transferencia de desplazamiento.
Similarmente se realiza la transformación de fuerzas axiales y cortantes en los elementos que concurren al nodo, la cual se presenta en forma gráfica en la figura 6.
Figura 6. Fuerzas sobre el nodo de los elementos que concurren en él.
Por equilibrio se encuentra
(16)
Y en forma matricial
(17)
En forma compacta
(18)
con
(19)
La ecuación anterior expresa las componentes de la fuerza que actúan en un lado del nodo orientadas según el sistema de referencia local asociado al elemento situado en ese mismo lado del nodo, en función de las componentes de la fuerza que actúa en el otro lado del nodo orientadas de acuerdo al sistema de referencia asociado con el elemento de este otro lado del nodo.
Reuniendo todas las variables de estado y las matrices de rotación anteriormente halladas, se obtiene la siguiente expresión, que relaciona las variables de estado de ambos lados de un nodo que une elementos con orientaciones diferentes.
(20)
que expandida es
(21)
y de manera compacta se expresa
(22)
Nótese que el ángulo es positivo en sentido anti-horario, midiéndose a partir del eje x correspondiente al elemento que llega al nodo. Para este cálculo los dos ejes se colocan en un mismo origen (sobre el nodo). Ver figura 7.
Figura 7. Detalle del cálculo del ángulo entre elementos.
Para el caso de un nodo con una carga axial PA, una transversal PT y un momento flector M0, con los sentidos indicados en la figura 8, la matriz de transferencia que se obtiene por consideraciones de equilibrio es la indicada por la ecuación 23.
Figura 8. Nodo sometido a cargas concentradas.
La matriz de cargas puntuales es
(23)
De lo anterior se concluye que en un pórtico como el de la figura Nº 1, el sistema de ecuaciones que lo modela es de la forma
(24)
o
(25)
con
(26)
que consiste en la multiplicación de matrices cuadradas de 7 * 7 en las que sus términos son propiedades geométricas, cargas y orientaciones de los elementos del pórtico. Esto conduce a un sistema de 6 ecuaciones y 12 variables de estado que corresponden a los desplazamientos, rotaciones, fuerzas y momentos en las secciones inicial y final del pórtico. De estas doce variables de estado, cualesquiera que sean las condiciones de contorno, siempre 6 de ellas son conocidas, por lo que el sistema de ecuaciones obtenido es determinado.
4. NODOS CON MÁS DE DOS ELEMENTOS
Otros casos que pueden presentarse son los pórticos en los que a un nodo concurren más de dos elementos, (nodos 2 y 5 del pórtico de la figura 9). Una manera de resolver estos casos es descomponer la estructura es subestructuras. A cada una de estas subestructuras de les aplica el método de la matriz de transferencia, y a las ecuaciones obtenidas se les adicionan las de compatibilidad y equilibrio de los nodos comunes, obteniéndose nuevamente un sistema determinado de ecuaciones. Para ilustrar el procedimiento se presenta el pórtico de la figura 9. En la figura 10 se muestran las subestructuras.
Figura 9. Pórtico de dos pisos. Las flechas indican la orientación de los elementos.
Figura 10. Despiece del pórtico en subsistemas.
Las ecuaciones para cada subestructura son
(27)
a las cuales se le agregan las ecuaciones de equilibrio y compatibilidad de las juntas 2 y 5, que relacionan las variables de estado z21 z22 y z23, y los vectores de estado z55 z53 y z56. En los vectores de estado escritos de la forma zij se ha utilizado la convención siguiente: el primer subíndice (i) expresa el número del nodo y el segundo (j) indica a cual elemento que concurre el nodo corresponden el vector de estado.
Para el establecimiento de las ecuaciones de equilibrio y compatibilidad de los desplazamientos del nodo se elabora la figura 11 que corresponde a un nodo genérico al cual llegan dos elementos y salen dos elementos. Lo que se obtendrá a continuación puede extrapolarse al caso de que sea un número mayor o menor de elementos que concurran al nodo.
Figura 11. Nodo con cuatro barras concurrentes.
Los elementos que llegan al nodo son los que su sección final se une al nodo, mientras que los que salen del nodo son aquellos que su sección inicial está unida al nodo. Para el caso de la figura se considera que los elemento 1 y 2 llegan al nodo y los elementos 3 y 4 salen.
Para el equilibrio se pueden realizar dos sumatorias de fuerzas en dos direcciones ortogonales, que por costumbre se han indicado en la figura como las direcciones horizontal x y vertical y.
Las fuerzas del elemento 1 que actúan sobre el nodo y referidas a la orientación del elemento son zfL1. Dichas fuerzas se pueden expresar en la direcciones xy mediante la rotación de ejes anteriormente indicada, esto es
(28)
(29)
Recuérdese que el ángulo es entre el eje del elemento y la dirección del eje sobre el cual se desean proyectar las fuerzas y/o desplazamientos, ambos situados sobre el mismo origen. Ver figura N° 7.
En forma análoga para el elemento 2
(30)
(31)
Las fuerzas del nodo en las direcciones xy sobre los elementos 3 y 4 que salen del nodo son respectivamente
(32)
(33)
(34)
(35)
El equilibrio del nodo es
(36)
Ecuación que puede ser expandida a una sumatoria de fuerzas de los elementos que concurren al nodo.
La compatibilidad se expresa en forma análoga basándose en la misma figura 11. Descomponiendo los desplazamientos de las secciones de los elementos que concurren al nodo
(37)
(38)
(39)
(40)
(41)
(42)
(43)
(44)
Por lo que las ecuaciones de compatibilidad son
(45)
Estas ecuaciones de equilibrio y de compatibilidad conjuntamente con las ecuaciones de las subestructuras forman un sistema de ecuaciones determinado.
5. EJEMPLOS
Considérese el pórtico de la figura 12 con la carga mostrada. Este ejercicio, resuelto por el Método de la Pendiente-Deflexión, está tomado de Wang (1953). En la Figura 13 se puede observar los sistemas de referencia locales del pórtico de la Figura 12.
Figura 12. Pórtico a dos aguas.
Figura 13. Orientación de los sistemas de referencia locales del pórtico de la figura 12.
Las ecuaciones del sistema están representadas por
(46)
donde
(47)
con l = 200cm, I = 100000 cm4; A = 10000 cm2.
(48)
con q = -12 cos a2 * sen a2 ; p = -12 cos2 a2; l = 223.6 cm; I= 200000 cm4, A=1000 cm2.
Siendo a2 = arctan (1/2).
(49)
con q = 12 cos a3 * sen a3 ; p = -12 cos2 a3; l = 223.6 cm; I = 200000 cm4, A = 1000 cm2.
Siendo a3 = arctan (1/2).
Para las matrices de rotación, los ángulos que se indican en la figura 14, que son los que se sustituirán en las matrices de rotación
Figura 14. Sentidos de los ángulos para las matrices de rotación.
(50)
(51)
(52)
con a12 = - arctan (1/2); a23 = - 2 * arctan (2) y a34 = - arctan (1/2)
Las condiciones de contorno son
d1 = 0; y1 = 0; q1 = 0; d5 = 0; y5 = 0; q5 = 0. (53)
La solución del sistema de ecuaciones conduce a los resultados indicados en la figura 15.
Figura 15. Sentido de las reacciones.
Para ilustrar el caso de más de dos elementos concurriendo a un nodo considérese el pórtico de la figura 16. Este ejercicio, resuelto por el Método de la Distribución de Momentos, aparece en Hsieh (1973)
Figura 16. Pórtico de tres tramos.
La estructura se discretiza como se muestra en la figura 17.
Figura 17. Discretización del pórtico de la figura 16. Las flechas paralelas a los elementos indican el sentido de avance de la transferencia.
Las subestructuras a considerar se muestran en la figura 18, y los sistemas de referencia locales en la figura 19.
Figura 18. Subsistemas del pórtico de la figura 19.
Figura 19. Sistemas de referencia de cada subsistema.
Las ecuaciones de las subestructuras son
(54)
Siendo Li la matriz indicada por la ecuación (4) con las correspondientes inercias, áreas y longitudes y Qi la matriz de carga puntual indicada por (23) con la carga axial y el momento concentrado nulos.
Las ecuaciones de equilibrio y compatibilidad de desplazamientos en los puntos de unión de los subsistemas corresponden a los nodos 2 y 3. La indicación de los elementos que llegan y salen de los nodos se indican en la figura 20. Estas ecuaciones son
(55)
Figura 20. Indicación de cuales elementos llegan y salen del nodo.
Con a = 0° para los elementos 1, 2 y 3; y a= +90° para los elementos 4 y 5.
Los valores de las variables de estado, resuelto el sistema de ecuaciones, se resume en la Tabla 1, cuyos sentidos verdaderos se muestran en la Figura 21.
Tabla 1. Valores de las variables de estado del segundo ejemplo.
A1 | M1 | V1 | A21 | M21 | V21 |
431 kg | 338 kg-m | -356kg | 431 kg | -730 kg-m | -356 kg |
A22 | M22 | V22 | A24 | M24 | V24 |
0 | -1308 kg-m | 3654 kg | -4010 kg | 578 kg-m | 3654 kg |
A32 | M32 | V32 | A33 | M33 | V33 |
0 | 1308 kg-m | 3654 kg | -431 kg | 730 kg-m | -356 kg |
A35 | M35 | V35 | A4 | M4 | V4 |
4010 kg | 578 kg-m | -431 kg | -431 kg | -338 kg-m | -356 kg |
A5 | M5 | V5 | A6 | M6 | V6 |
-4010 kg | -284 kg-m | -431 kg | 4010 kg | -284 kg-m | -431 kg |
Figura 21. Sentidos reales de las fuerzas según resultados obtenidos.
6. CONCLUSIONES
Se obtuvieron las matrices de rotación para poder aplicar el método de la matriz de transferencia en elementos rectos no alineados.
Se obtuvieron las ecuaciones de compatibilidad y equilibrio basadas en matrices de rotación para los nodos donde concurren más de dos elementos para que el método de la matriz de transferencia pueda ser aplicado a cualquier tipo de pórtico plano de elementos rectos.
Los ejemplos mostrados demuestran la factibilidad de la aplicación del método en los pórticos planos.
7. REFERENCIAS
1. Dimarogonas, A. 1996. Vibrations for Engineers. Second Edition. Prentice Hall. New Jersey. 825 p. [ Links ]
2. Hsieh, Y. 1973. Teoría Elemental de Estructuras. Pág. 285-287. Prentice Hall Interamericana. [ Links ]
3. Wang, C. 1953. Statically Indeterminate Structures. Pág. 190-195. McGraw-Hill. [ Links ]