En esta investigación se obtiene el campo de desplazamientos y deformaciones para una viga de sección abierta y pared delgada formada por laminados compuestos incluyendo en su cinemática las deformaciones tangenciales por cortante por flexión y torsión. La ecuación de equilibrio del sistema se define a través de una formulación variacional, mostrando que incluir las deformaciones de corte por torsión hace posible formular elementos finitos de clase Co. Las relaciones tensión-deformación se obtienen en la sección transversal de una viga de sección abierta formada por laminados compuestos extendiendo la teoría de laminados de primer orden - FSDT: first-order shear deformation) y usando una «condición de esfuerzos nulos en el contorno». Se formulan 3 elementos finitos unidimensionales BSW (beam with shear and warping) con grado de continuidad Co para el estudio de vigas de sección abierta formadas por laminados compuestos. Mediante el uso de los elementos BSW se evalúa la influencia de la estrategia de integración ante el fenómeno de bloqueo por tensiones tangenciales y la velocidad de convergencia para los desplazamientos y giros. Se compara la integración de la matriz de rigidez cuando se aplican técnicas de integración exacta y reducida. Finalmente, se realizan 2 ejemplos donde se estudia el comportamiento a torsión pura y flexotorsión de una viga con sección en C empotrada en un extremo y se comparan los resultados obtenidos con otros de la literatura.
In this paper we derive the field of displacements and strains for thin-walled open composite beams with composite laminated material including in their kinematics flexural and torsional shear deformations effects. The equilibrium equations are defined through the variational formulation and show that is possible to formulate Co finite elements taking into account the torsional shear deformation. Stress-strain relationships for the cross-section of thin-walled composite beams are obtained by extending first-order laminate (FSDT: first-order shear deformation) theory and using a «free stress resultant condition at the boundary». Three different one-dimensional finite elements with Co continuity are formulated for the study of thin-walled open composite beams and they are labelled as BSW (beam with shear and warping). The influence of the integration strategy in the BSW elements is evaluated via the shear-locking phenomenon and the rate of convergence for displacements and rotations. The stiffness matrix integration is compared using exact and reduced integration methods. Examples of pure torsion and flexo-torsion in a cantilever composite beam are performed. Numerical results are compared to those reported by other authors.
El uso de vigas de sección abierta formadas por laminados compuestos experimentó un aumento importante en ingeniería civil, naval y mecánica desde comienzos de la década de 1990 gracias a las propiedades específicas que exhiben: alta relación rigidez/peso y resistencia mecánica/peso, alta resistencia a la degradación ambiental y baja conductividad térmica. Además, por su ligereza, estas vigas son muy adecuadas para estructuras donde el peso constituye una variable fundamental en el proceso de diseño [1–3].
Desde la teoría de Vlasov [4], aplicable a vigas de pared delgada en materiales isótropos, se han realizado numerosas investigaciones, entre las que destacan los estudios realizados por Bauld y Tzeng [5], quienes extendieron la teoría de Vlasov a laminados compuestos, aunque no consideraron las deformaciones transversales en sus análisis; Bank y Bednarczyk [6] y Barbero et al. [7] desarrollaron expresiones simples para el análisis a flexión, torsión y alabeo de vigas en compuestos; Wu y Sun [8] incluyeron en su cinemática las deformaciones tangenciales por flexión y torsión mostrando la utilidad de tenerlas en cuenta durante los análisis; Massa y Barbero [9] consideraron las deformaciones transversales por flexión en sus análisis; Kollar y Pluzsik [10] y Pluzsik y Kollar [11] mostraron la importancia de incluir las deformaciones tangenciales por flexión en los análisis de vigas con alabeo restringido; Lee y Lee [12] desarrollaron un modelo analítico para el estudio de vigas a flexión y torsión que incluye alabeo, pero no consideraron las deformaciones tangenciales, lo que los llevó a utilizar una combinación de elementos lagrangianos y hermíticos en sus análisis. Posteriormente, Lee [13] presentó un modelo general para el estudio de vigas sometidas a flexión incluyendo en su formulación deformaciones tangenciales por flexión y «deformaciones de corte por alabeo», lo que le permitió desarrollar elementos finitos lagrangianos. Sin embargo, Lee utilizó una fórmula para el cálculo del alabeo que solo es aplicable a secciones transversales doblemente simétricas y no contempló la torsión en sus estudios. Kim et al. [14] presentaron soluciones exactas para el estudio a torsión de vigas de sección abierta sin considerar las deformaciones tangenciales dentro de su formulación, lo que les llevó a usar elementos hermíticos de tercer orden en sus análisis. Back y Will [15] extendieron el modelo general de Lee [13] mejorando la fórmula para el cálculo del alabeo haciéndola válida también para secciones con una sola simetría, pero tampoco estudiaron la torsión. Finalmente, Feo y Mancusi [16] sugirieron una aproximación al estudio de vigas en materiales compuestos mediante el uso de elementos finitos hermíticos, contemplaron en su análisis las deformaciones tangenciales por flexión y resolvieron la función de alabeo a través de funciones de aproximación dependientes de la coordenada curvilínea de la línea media de la sección transversal.
En este trabajo se compatibilizan de forma adecuada los campos de deformaciones de los laminados con la sección transversal de viga partiendo de las investigaciones de Lee [13] y Back y Will [15], y extendiendo la hipótesis de « condición de esfuerzos nulos en el contorno » [7,14] a laminados de primer orden - FSDT [17]. Se formulan elementos finitos tipo BSW que incluyen las deformaciones tangenciales debidas al cortante por flexión y torsión para el análisis a flexión y torsión de vigas de sección abierta y laminados compuestos. A diferencia de Lee [13] y Back y Will [15], el valor de la función de alabeo se obtiene a través de un algoritmo de búsqueda por grafos [18], generalizando la formulación a cualquier configuración de sección abierta. Los elementos se han formulado para funciones de aproximación lineal (L-BSW), cuadrática (Q-BSW) y cúbica (C-BSW), comprobándose la validez de los resultados frente a los obtenidos usando elementos cuadriláteros S9R5 [19].
En este trabajo también se estudia la influencia del método de integración numérica (exacta y reducida) en el bloqueo producido por tensiones tangenciales y sus consecuencias en la convergencia de la solución. Se modeliza una viga en voladizo con sección C y se comparan los resultados obtenidos a flexión y torsión con otros de la literatura [14].
2Teoría general de vigas2.1CinemáticaEl desarrollo del modelo de viga requiere 2 sistemas ortogonales de coordenadas (fig. 1), el sistema cartesiano (x, y, z) donde el eje x corresponde al eje longitudinal de la viga, y el sistema local (x¯,n,s) para un punto P sobre la línea media de la sección transversal.
En el sistema local el eje x¯ es paralelo al eje x, por lo que en adelante se usarán indistintamente. El eje n es normal a la línea media de la sección en el punto P. Los 2 sistemas están relacionados por el ángulo α medido desde el eje y hasta el eje s en el sentido contrario a las agujas del reloj. La distancia entre el punto P y el centro de giro c de la sección transversal se evalúa como la suma de los vectores ρn y ρs, paralelos a los ejes s y n, respectivamente.
El campo de desplazamientos se obtiene suponiendo que el ángulo de giro por torsión es pequeño, el contorno de la sección transversal no se deforma en su plano, y las líneas rectas normales al eje longitudinal de la viga permanecen rectas antes de la deformación pero no necesariamente normales al eje después. Así, para un punto P con coordenadas (y, z, n) se tiene:
donde uo, vc y wc son los desplazamientos de sólido rígido en los ejes x, y y z, respectivamente; θx, θy y θz son las rotaciones alrededor de los ejes x, y y z; ϕω es la intensidad de alabeo por torsión a la que está sometida la sección; ωs es la función de alabeo primario definida desde el centro de esfuerzos cortantes c y evaluada sobre la línea media de la sección transversal; y ωn es la componente de alabeo secundario o alabeo en el espesor.
De acuerdo con Wu y Sun [8], las rotaciones θy y θz se definen en función de las deformaciones tangenciales debidas a flexión γxzo y γxyo, la intensidad de alabeo ϕω, en función de la variación en el ángulo de giro por torsión ∂xθx, y las deformaciones tangenciales generadas por la torsión γt:
De la expresión (2) se deduce que cuando no hay deformaciones tangenciales por torsión, o estas son despreciables (γt≈0), la intensidad del alabeo equivale a la variación del ángulo de torsión por unidad de longitud (ϕω=∂xθx) y se recupera la teoría de Vlasov [4].
La función normalizada de alabeo primario por torsión ωs y la función de alabeo secundario ωn se definen como:
La función de alabeo ωs se evalúa de forma automática para cada punto de la línea media de la sección transversal utilizando algoritmos de búsqueda por grafos [18,20].
2.2DeformacionesSiguiendo la teoría de pequeños desplazamientos, se define el campo de deformaciones como:
donde ¿xo es la deformación axial; γxyo y γxzo son las deformaciones tangenciales por flexión; γt es la deformación angular por torsión; κ son las curvaturas debido a flexión (κy, κz), torsión (κxs) y alabeo (κω). De las ecuaciones (1) y (4) se deducen las siguientes deformaciones y curvaturas [20]:
Para el caso en que se supone despreciable la contribución de las deformaciones angulares por torsión (γt≈0 y ϕω=∂xθx), el valor de la curvatura por torsión será κxs=2∂xθx, recuperándose las expresiones empleadas por Lee y Lee [12], Back y Will [15], Sheikh y Thomsen [21].
Agrupando los términos del campo de deformaciones 4 y aplicando un cambio de variables, resulta:
donde ¿¯xo es una deformación axial, γ¯xso y γ¯xno son deformaciones tangenciales, y κ¯x y κ¯xs son curvaturas de flexión y torsión en la línea media de la sección.2.3Equilibrio seccional mediante una formulación variacional
La energía potencial de deformación se define como:
Al reemplazar las deformaciones de la ecuación (4) se obtiene:
Así, la variación de la energía potencial de deformación del sistema se evalúa a partir de la ecuación (8), resultando la expresión en función de los esfuerzos seccionales [20]:
donde Nx es el esfuerzo axial, Mz y My son los momentos flectores en los ejes z e y, Mt es el momento torsor en el espesor de la pared de la sección, Bω es el momento de alabeo (bimomento), Qy y Qz son las fuerzas de corte en las direcciones y y z, y T es el torsor asociado a las deformaciones angulares, definidas en el área de la sección transversal. Las expresiones de estos esfuerzos son [20]: 2.4Ecuaciones constitutivas para el laminado compuesto
La ecuación constitutiva para un laminado de material compuesto en su sistema de coordenadas local (xˆ,yˆ,nˆ), siguiendo la teoría FSDT (first-order shear deformation)[17] se escribe como [20]:
donde Aij, Bij, Dij y Hij son los coeficientes de las matrices de rigidez membranal, de acoplamiento, de flexión y de tensiones tangenciales, definidas por [17,20]:
La función de ponderación f(nˆ) permite considerar la distribución parabólica de las tensiones en el espesor del laminado.
Los ejes xˆ y nˆ del sistema de coordenadas local del laminado (xˆ,yˆ,nˆ) son paralelos y tienen el mismo sentido que los ejes x¯ y n del sistema de coordenadas local en la línea media de la sección transversal (x¯,n,s). De esta forma:
Las relaciones tensión-deformación (11) deben simplificarse con el fin de hacerlas compatibles con las componentes del vector de deformaciones a nivel de la sección transversal de las vigas de pared delgada en laminados compuestos (6).
Lee y Lee [12], Sheikh y Thomsen [21] y Back y Will [15] sugieren una «condición de tensiones nulas en el contorno», donde simplifican las componentes de tensión a nivel de las capas del laminado suponiendo que las tensiones laterales (normal y de corte) son nulas. Esta condición, sin embargo, tiene la desventaja de no contemplar la interrelación entre las diferentes capas que conforman el laminado. Por otra parte, Barbero et al. [7] y Kim et al. [14] proponen una «condición de esfuerzos nulos en el contorno» para vigas laminadas esbeltas CLT en la que simplifican las componentes de los esfuerzos a nivel de laminado considerando que Nyˆ=Myˆ=0 en el contorno, lo que permite considerar la interrelación de las capas del laminado. Sin embargo, por tratarse de laminados esbeltos, no se captura de forma precisa el comportamiento a torsión. En este trabajo se usa la propuesta de Barbero et al. [7] y Kim et al. [14] y de forma adicional se impone que Qyˆ=0 en el borde (fig. 2), con el fin de incluir la interrelación entre las capas del laminado.
Al imponer esta nueva condición de esfuerzo cortante nulo en el contorno, es posible evaluar las deformaciones ¿yˆo, κyˆ y γyˆnˆo como una combinación lineal de las componentes de deformación conocidas (¿xˆo, γxˆyˆo, κxˆ, κxˆyˆ y γxˆnˆo).
Imponiendo que Nyˆ=Myˆ=Qyˆ=0 en el borde del laminado en las relaciones tensión-deformación para un laminado, se obtiene la siguiente relación entre las fuerzas internas y las deformaciones generalizadas de la sección transversal de la viga de laminados compuestos:
En forma compacta:
Los coeficientes de esta matriz de rigidez dependen de la geometría y de las propiedades de la sección transversal del elemento, tal como se muestra en el Apéndice 1.1.
3Formulación de elementos finitos BSWLos elementos finitos formulados (denominados BSW por beam with shear and warping) para el estudio de vigas de sección abierta en laminados compuestos son unidimensionales con grado de continuidad Co. Tienen 7 grados de libertad por nodo y consideran dentro de su cinemática los desplazamientos debidos a las deformaciones por corte inducidas por flexión y por torsión, así como los desplazamientos longitudinales generados por el alabeo restringido.
En esta investigación se han desarrollado 3 elementos: L-BSW (lineal de 2 nodos), Q-BSW (cuadrático de 3 nodos) y C-BSW (cúbico de 4 nodos). La formulación que se presenta sigue la metodología empleada por Oñate [22] para elementos de viga en 3 dimensiones, por Lee [13] para análisis a flexión de vigas de materiales compuestos, y por Back y Will [15] para vigas de materiales compuestos con sección T.
3.1DesplazamientosEl vector de movimientos de un punto u:
dondeyson, respectivamente, la matriz de funciones de forma (N), la matriz identidad de 7×7 (I7), el número de nodos del elemento (n), el vector de movimientos de un punto del elemento (a(e)) y de un nodo (ai(e)).3.2Matriz de deformación del elemento
La matriz de deformación B relaciona las deformaciones generalizadas de la sección transversal de la viga ¿ˆ con el vector de movimientos nodales a(e) o con los siguientes movimientos:
Es importante destacar que las deformaciones de la sección transversal de la viga (ecuaciones 4 y 5) se definen en términos de las primeras derivadas del campo de desplazamientos, permitiendo formular elementos finitos con un grado de continuidad Co.
La expresión del operador D se deduce de las ecuaciones (4) y (5). La matriz de deformación B se obtiene por:
Para un nodo i:
3.3Matriz de rigidez del elemento
La matriz de rigidez se obtiene del equilibrio del elemento. Su expresión tiene la forma usual [22,20]:
donde E˜ es la matriz constitutiva que relaciona las fuerzas internas y las deformaciones generalizadas de la sección transversal de la viga de laminados compuestos (13a). Los componentes de E˜ se incluyen en el Apéndice 1.1. El tamaño de la matriz K(e) es (7×n)×(7×n), siendo n el número de nodos del elemento.4Integración numérica de la matriz de rigidez
En este apartado se comparan los resultados obtenidos al aplicar 2 métodos diferentes de integración numérica (exacta y reducida) en las componentes a flexión y torsión de la matriz de rigidez de los elementos formulados (L-BSW, Q-BSW y C-BSW).
El estudio que se presenta incluye 2 tipos de análisis de una viga en voladizo con sección en doble T. El primer análisis evalúa el bloqueo por tensiones tangenciales de los elementos mediante la variación de la relación de esbeltez de viga. El segundo análisis estudia la convergencia de la solución numérica al variar el número de elementos usados en la discretización.
En este trabajo se adoptan la geometría y las propiedades usadas por Lee [13], Kim et al. [14], Back y Will [15] y Sheikh y Thomsen [21] en sus análisis.
La figura 3 muestra la sección y las propiedades de la viga de laminados compuestos usada en los análisis. Los laminados de la sección transversal son unidireccionales con 16 capas ([φ]16) de láminas en fibra de vidrio embebida en matriz epóxica (Gl.-Ep, Reddy [17]) de 0,13mm de espesor cada una [20].
Los resultados obtenidos con los elementos desarrollados se comparan con los resultados de modelos en Abaqus utilizando elementos cuadriláteros tipo S9R5 [19] con mallas de entre 648 y 960 elementos, dependiendo de la esbeltez de la viga analizada. En la figura 4 se muestra el modelo de 960 elementos S9R5 utilizado para verificar los resultados del análisis de torsión de una viga de 1m de longitud (esbeltez λ=L/20).
![Modelo de 960 elementos cuadriláteros S9R5 [19], usado para simular una viga I de 1.000mm de longitud sometida a torsión.](https://static.elsevier.es/multimedia/02131315/0000003000000004/v1_201410050003/S0213131514000030/v1_201410050003/es/main.assets/gr4.jpeg?xkr=ue/ImdikoIMrsJoerZ+w997EogCnBdOOD93cPFbanNdXKeWo41MsZkw/vID64cuEiPuJOyFxneHZSkizm96flfSXRuq9PKMsoqi+EW0rIT4RmdB0rXJVLjYZ01/kXBRmtd2Jz2dwhLgFDFVD1i8yuF7qPcGwRLohekQXheHWaFVDknP4zZ5Y6HSI94Fl/K9TfeRz03wVwnWSzsM0AFl1Y88uFYufS3UTvlRbebpR7sf+2R+wJjTDqOnIf2aUfOWdawIG4Ruw8uRaqIPzXjTJApGa+f4ZKhhPPRk7sws6eDd2zVF9GQQrfL6uzZphHE7pjgOLWXe8pwWmipna88lNkQ== "Modelo de 960 elementos cuadriláteros S9R5 [19], usado para simular una viga I de 1.000mm de longitud sometida a torsión.")
Modelo de 960 elementos cuadriláteros S9R5 [19], usado para simular una viga I de 1.000mm de longitud sometida a torsión.
El análisis a flexión se realiza calculando el desplazamiento vertical w del extremo libre de la viga en voladizo donde se aplica una carga vertical de 2kN. En el análisis a torsión (fig. 4) se aplica un momento torsor de 1kNmm en el extremo libre y se obtiene el giro por torsión θ de la sección transversal en el punto de aplicación del momento.
4.1Análisis de bloqueo por tensiones tangencialesSe entiende por bloqueo la incapacidad intrínseca que tiene un elemento de satisfacer la condición de deformación transversal nula cuando su longitud tiende al infinito [22]. Esto conduce a una sobrerrigidez numérica a medida que la longitud del elemento aumenta. Una alternativa para sortear este problema es subintegrar los términos de la matriz de rigidez que incluyen las deformaciones transversales (debidas a la flexión y a la torsión) utilizando un número de puntos de integración inferior al necesario para su cálculo exacto [22].
Concretamente, el número de puntos de integración para las integraciones exacta (E) y reducida (R) para los diferentes elementos estudiados es el siguiente: L-BSW: 2 (E), 1 (R), Q-BSW: 3 (E), 2 (R), C-BSW: 4 (E), 3 (R), donde el numero indica la cantidad de puntos de integración para el cálculo de los términos de la matriz de rigidez que contienen las deformaciones transversales. Es destacable que el resto de términos de la matriz de rigidez se calcula de manera exacta, utilizando un número de puntos de integración igual al de la integración reducida [22].
Con el ánimo de estudiar la capacidad para satisfacer la condición de deformaciones transversales nulas de los elementos implementados, se analizan modelos de un solo elemento finito BSW y se varía su longitud (manteniendo su altura constante) alcanzando relaciones de esbeltez λ entre 1 y 20.
Los resultados obtenidos aplicando integración numérica exacta y reducida en las componentes de flexión y torsión de la matriz de rigidez se normalizan respecto a los obtenidos usando los modelos con elementos cuadriláteros S9R5 [19], que se comparan en las tablas 1 y 2. La primera tabla corresponde al análisis a flexión y la segunda al análisis a torsión. Los datos incluidos corresponden a las vigas modelizadas con secuencia de apilado unidireccional [0]16.
Viga I de laminados [0]16 sometida a flexión. Comparación de resultados en análisis de bloqueo entre integración exacta e integración reducida.
| λ=L/h | wBSW/wS9R5 | |||
| L-BSW | Q-BSW | Q-BSW | C-BSW | |
| (Int. red) | (Int. exac) | (Int. red) | (Int. exac) | |
| 1 | 0,5520 | 0,5711 | 0,5713 | 0,5713 |
| 2 | 0,5706 | 0,6294 | 0,6312 | 0,6312 |
| 3 | 0,6180 | 0,7168 | 0,7237 | 0,7237 |
| 4 | 0,6536 | 0,7799 | 0,7957 | 0,7957 |
| 5 | 0,6783 | 0,8193 | 0,8468 | 0,8468 |
| 6 | 0,6952 | 0,8417 | 0,8828 | 0,8828 |
| 7 | 0,7073 | 0,8528 | 0,9083 | 0,9083 |
| 10 | 0,7272 | 0,8527 | 0,9510 | 0,9510 |
| 15 | 0,7397 | 0,8258 | 0,9771 | 0,9771 |
| 20 | 0,7442 | 0,8031 | 0,9870 | 0,9870 |
Viga I de laminados [0]16 sometida a torsión. Comparación de resultados en análisis de bloqueo entre integración exacta e integración reducida
| λ=L/h | θBSW/θS9R5 | |||
| L-BSW | Q-BSW | Q-BSW | C-BSW | |
| (Int. red) | (Int. exac) | (Int. red) | (Int. exac) | |
| 1 | 3,0683 | 3,4272 | 3,4411 | 3,4410 |
| 2 | 1,2825 | 1,5362 | 1,5756 | 1,5755 |
| 3 | 1,0575 | 1,2677 | 1,3420 | 1,3417 |
| 4 | 0,9747 | 1,1418 | 1,2491 | 1,2486 |
| 5 | 0,9314 | 1,0599 | 1,1937 | 1,1930 |
| 6 | 0,9078 | 1,0039 | 1,1578 | 1,1567 |
| 7 | 0,8937 | 0,9629 | 1,1309 | 1,1295 |
| 10 | 0,8849 | 0,8953 | 1,0830 | 1,0806 |
| 15 | 0,9118 | 0,8585 | 1,0457 | 1,0412 |
| 20 | 0,9535 | 0,8528 | 1,0302 | 1,0237 |
Las tablas 1 y 2 muestran que los resultados que más se aproximan a los obtenidos con las modelizaciones con elementos cuadriláteros S9R5[19] para flexión y torsión son los de los elementos C-BSW con integración exacta y Q-BSW con integración reducida, y cabe señalar que entre ellos no hay diferencias apreciables. Después destacan los resultados del elemento Q-BSW con integración exacta y, por último, los resultados de L-BSW con integración reducida. Los resultados para el elemento C-BSW con integración reducida fueron coincidentes con los de integración exacta.
Vale la pena remarcar que en ninguno de los casos estudiados se obtuvo el bloqueo de la solución.
4.2Análisis de convergenciaSe estudia la convergencia de los elementos formulados (L-BSW, Q-BSW y C-BSW) para los casos en que se utiliza integración numérica exacta y reducida de las componentes a flexión y torsión de la matriz de rigidez del elemento. El análisis se realiza en vigas de longitud y altura constante (esbeltez λ=20) variando el número de elementos BSW usados en cada modelización y se comparan con los resultados obtenidos usando un modelo de 960 elementos cuadriláteros S9R5[19].
Aunque el análisis de convergencia se ha realizado en vigas con secuencias de apilado [φ]16 (φ∈[015, 30, 45, 75, 90]), en este apartado solamente se muestran los resultados obtenidos para la secuencia [45]16, por ser la que presenta mayores diferencias porcentuales con respecto a los resultados obtenidos en modelos con elementos cuadriláteros S9R5[19]. En [20] se puede encontrar más información sobre el resto de resultados.
Los resultados de la integración numérica exacta y reducida se comparan en las tablas 3 y 4. La primera tabla corresponde al análisis a flexión y la segunda al análisis a torsión.
Viga de laminados [45]16 sometida a flexión. Convergencia de resultados utilizando integración exacta y reducida
| Núm. elementos | wBSW/wS9R5 | |||
| L-BSW | Q-BSW | Q-BSW | C-BSW | |
| (Int. red) | (Int. exac) | (Int. red) | (Int. exac) | |
| 1 | 0,7224 | 0,7426 | 0,9613 | 0,9613 |
| 2 | 0,9016 | 0,9178 | 0,9613 | 0,9613 |
| 3 | 0,9350 | 0,9471 | 0,9613 | 0,9613 |
| 4 | 0,9466 | 0,9552 | 0,9613 | 0,9613 |
| 5 | 0,9517 | 0,9583 | 0,9613 | 0,9613 |
| 6 | 0,9547 | 0,9598 | 0,9613 | 0,9613 |
| 7 | 0,9563 | 0,9603 | 0,9613 | 0,9613 |
| 8 | 0,9578 | 0,9608 | 0,9613 | 0,9613 |
Viga de laminados [45]16 sometida a torsión. Convergencia de resultados utilizando integración exacta y reducida
| Núm. elementos | θBSW/θS9R5 | |||
| L-BSW | Q-BSW | Q-BSW | C-BSW | |
| (Int. red) | (Int. exac) | (Int. red) | (Int. exac) | |
| 1 | 1,0224 | 0,8370 | 0,9924 | 0,9766 |
| 2 | 0,9862 | 0,9475 | 0,9889 | 0,9880 |
| 3 | 0,9880 | 0,9713 | 0,9889 | 0,9889 |
| 4 | 0,9880 | 0,9801 | 0,9889 | 0,9889 |
| 5 | 0,9880 | 0,9836 | 0,9889 | 0,9889 |
| 6 | 0,9889 | 0,9862 | 0,9889 | 0,9889 |
| 7 | 0,9889 | 0,9871 | 0,9889 | 0,9889 |
| 8 | 0,9889 | 0,9871 | 0,9889 | 0,9889 |
En líneas generales, todos los elementos estudiados convergen a la solución obtenida con la malla de elementos cuadriláteros S9R5.
En los análisis a flexión, los resultados del elemento Q-BSW con integración exacta se aproximan más a los de elementos C-BSW que los resultados de L-BSW con integración reducida, mientras que en los análisis a torsión sucede lo contrario.
En las figuras 5 y 6 se comparan los resultados obtenidos con los elementos L-BSW, Q-BSW y C-BSW; el parámetro de comparación es el número de nodos. No se incluyen los resultados de los modelos con elementos lineales L-BSW e integración exacta por presentar bloqueo de la solución.
Las figuras 5 y 6 muestran que los modelos con elementos L-BSW e integración reducida tienen un mejor comportamiento que los elementos Q-BSW con integración exacta para el mismo número de nodos.
También se observa que no hay diferencias importantes entre los resultados obtenidos con los elementos Q-BSW con integración reducida y con los elementos C-BSW.
5Análisis de una viga en voladizo con sección en CSe estudia el comportamiento a flexión y torsión de una viga con sección C formando laminados compuestos (fig. 7), empotrada en un extremo y sometida a carga vertical excéntrica que genera un momento torsor negativo de 1kNmm.
Las propiedades de la sección se muestran en la figura 8.
El análisis se realiza con 10 elementos de 2 nodos (lineales) e integración reducida L-BSW.
Al igual que en el apartado 4, los laminados que conforman la viga tienen 16 capas de láminas en fibra de vidrio embebida en matriz epóxica (Gl.-Ep, Reddy [17]) de 0,13mm de espesor cada una (Vargas [20]).
Los resultados obtenidos usando 10 elementos L-BSW se comparan con los resultados obtenidos con 400 elementos cuadriláteros S9R5 (tipo DK, 9 nodos, 5 grados de libertad e integración reducida [19]).
Se analiza la variación del giro por torsión para laminados unidireccionales [φ]16 y bidireccionales [±φ]4s cuando cambia el ángulo de orientación de las fibras que componen la sección transversal de la viga.
- •
Laminados unidireccionales [φ]16. La figura 9 muestra los resultados cuando se aplica la carga puntual excéntrica. Se observan diferencias entre los modelos comparados de un 4,3% en el giro por torsión y de un 8,5% en el desplazamiento vertical.
- •
Laminados [±φ]4s. La figura 10 muestra los resultados cuando se aplica la carga puntual excéntrica. Se observa que las modelizaciones con elementos unidimensionales (L-BSW) presentan diferencias de un 2,7% en el giro por torsión y de un 2,3% en el desplazamiento vertical con respecto a las modelizaciones con elementos cuadriláteros S9R5 [19].
![Variación de los resultados en el extremo libre de la viga C en voladizo para diferentes orientaciones de las fibras en laminados [±φ]4s (carga vertical excéntrica).](https://static.elsevier.es/multimedia/02131315/0000003000000004/v1_201410050003/S0213131514000030/v1_201410050003/es/main.assets/gr10.jpeg?xkr=ue/ImdikoIMrsJoerZ+w997EogCnBdOOD93cPFbanNdXKeWo41MsZkw/vID64cuEiPuJOyFxneHZSkizm96flfSXRuq9PKMsoqi+EW0rIT4RmdB0rXJVLjYZ01/kXBRmtd2Jz2dwhLgFDFVD1i8yuF7qPcGwRLohekQXheHWaFVDknP4zZ5Y6HSI94Fl/K9TfeRz03wVwnWSzsM0AFl1Y88uFYufS3UTvlRbebpR7sf+2R+wJjTDqOnIf2aUfOWdawIG4Ruw8uRaqIPzXjTJApGa+f4ZKhhPPRk7sws6eDd2zVF9GQQrfL6uzZphHE7pjgOLWXe8pwWmipna88lNkQ== "Variación de los resultados en el extremo libre de la viga C en voladizo para diferentes orientaciones de las fibras en laminados [±φ]4s (carga vertical excéntrica).")
En las referencias [22,20] se muestran otros ejemplos de aplicación del elemento L-BSW al cálculo de vigas de sección abierta con materiales compuestos.
6ConclusionesSe han formulado, implementado y validado elementos finitos tipo BSW para el análisis a flexión y torsión en vigas de sección abierta y constituidas por laminados compuestos. La cinemática de la viga incluye las deformaciones de corte por flexión y torsión. Los elementos usan algoritmos de búsqueda por grafos [18] para el cálculo de la función de alabeo primario de tal manera que son independientes de la sección de la viga. Los elementos desarrollados en esta investigación son de 3 tipos, según las funciones de interpolación de los movimientos: lineal L-BSW, cuadrático Q-BSW y cúbico C-BSW. Los resultados se han comparado con los obtenidos con elementos cuadriláteros S9R5 [19] y con resultados de la literatura.
El estudio de la influencia entre aplicar integración numérica exacta y reducida de las componentes a flexión y torsión de la matriz de rigidez del elemento ha mostrado que los elementos L-BSW presentan bloqueo con integración exacta, mientras que los de aproximación Q-BSW y C-BSW no. También se observa que se requieren al menos 5 elementos de aproximación cuadrática Q-BSW con integración exacta para obtener los mismos resultados que utilizando elementos de aproximación cúbica C-BSW.
El fenómeno de bloqueo se eliminó en los elementos de aproximación lineal L-BSW utilizando integración reducida. En este caso, esto equivale a calcular todos los términos de la matriz de rigidez con un solo punto de integración.
Los elementos de aproximación cuadrática mejoran sus resultados al aplicar integración reducida y se aproximan a los resultados obtenidos con elementos de aproximación cúbica C-BSW.
Los elementos C-BSW no presentan variaciones apreciables en sus resultados cuando cambia el tipo de integración numérica (reducida o exacta).
Los elementos de aproximación cuadrática Q-BSW convergen con 2 elementos a los resultados obtenidos con modelos de elementos de aproximación cúbica. Estos últimos proporcionan buenos resultados utilizando un solo elemento.
El elemento que mejor convergencia presenta a un menor número de nodos es el de aproximación cuadrática Q-BSW con integración reducida. Sus resultados no presentan diferencias apreciables con respecto a los obtenidos con elementos con aproximación cúbica C-BSW.
En el estudio de una viga con sección en C, la diferencia en los resultados utilizando los elementos de aproximación lineal L-BSW e integración reducida (10 elementos en cada modelo de viga) frente a los modelos de 400 elementos cuadriláteros S9R5 fue inferior al 10% en los casos de torsión pura y flexotorsión.
Se agradece al Ministerio de Fomento de España por subvencionar parte de este trabajo a través de las ayudas del Programa del Plan Nacional I+D en la Construcción: 2004/32 y C9/2006. Este trabajo fue también parcialmente financiado por el proyecto SAFECON del European Research Council.
