Scholarly article on topic 'Abolladura de estructuras de barras de sección tipo doble-T'

Abolladura de estructuras de barras de sección tipo doble-T Academic research paper on "Educational sciences"

CC BY-NC-ND
0
0
Share paper
Keywords
{"Abolladura pandeo local" / "Carga crítica" / Inestabilidad / "Placa delgada" / "Local buckling" / "Critical load" / Instability / "Thin plate"}

Abstract of research paper on Educational sciences, author of scientific article — M. Cacho-Pérez, A. Lorenzana Ibán

Resumen En este trabajo se expone una metodología de varias etapas para determinar, con mínimo esfuerzo computacional el factor de carga crítico que produce el fenómeno de inestabilidad local de abolladura de las paredes de los perfiles tipo doble-T. Para ello, se estudia a nivel local como placa delgada cada una de las paredes que forman el perfil. Se constata que es suficiente con analizar las zonas más solicitadas (alas y alma) de la estructura. Abstract A stepped methodology for the evaluation, with limited computational effort critical buckling load factor which produces the instability of the walls of double-T sections is proposed in this work. For this purpose, each wall of the cross-section studied in a local level as a thin plate. It is ratified that it is enough to check the most compressed flange and the web of the most loaded beams of the structure.

Academic research paper on topic "Abolladura de estructuras de barras de sección tipo doble-T"

G ll'PW^M AIIIIII .E IN PRESS

Rev. int. métodos numér. cálc. diseño ing. 2016;xxx(xx):xxx-xxx

ELSEVIER

Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería

www.elsevier.es/rimni

Abolladura de estructuras de barras de sección tipo doble-T

M. Cacho-Pérez * y A. Lorenzana Ibán

Escuela de Ingenierías Industriales, Universidad de Valladolid, Paseo del Cauce 59,47011, Valladolid, Spain INFORMACIÓN DEL ARTÍCULO RESUMEN

Historia del artículo: Recibido el 10 de octubre de 2014 Aceptado el 4 de enero de 2016 On-line el xxx

Palabras clave: Abolladura pandeo local Carga crítica Inestabilidad Placa delgada

En este trabajo se expone una metodología de varias etapas para determinar, con mínimo esfuerzo computacional el factor de carga crítico que produce el fenómeno de inestabilidad local de abolladura de las paredes de los perfiles tipo doble-T. Para ello, se estudia a nivel local como placa delgada cada una de las paredes que forman el perfil. Se constata que es suficiente con analizar las zonas más solicitadas (alas y alma) de la estructura.

© 2016 CIMNE (Universitat Politécnica de Catalunya). Publicado por Elsevier España, S.L.U. Este es un artículo Open Access bajo la licencia CC BY-NC-ND (http://creativecommons.org/licenses/by-nc-nd/4.0/).

Keywords: Local buckling Critical load Instability Thin plate

Local buckling of beam structures with double-T type sections

ABSTRACT

A stepped methodology for the evaluation, with limited computational effort critical buckling load factor which produces the instability of the walls of double-T sections is proposed in this work. For this purpose, each wall of the cross-section studied in a local level as a thin plate. It is ratified that it is enough to check the most compressed flange and the web of the most loaded beams of the structure.

© 2016 CIMNE (Universitat Politécnica de Catalunya). Published by Elsevier España, S.L.U.

This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

1. Introducción

Para los elementos estructurales en los que el esfuerzo predominante es el momento flector el mejor aprovechamiento del material se consigue cuando el fallo en el mismo se produce por agotamiento de su capacidad resistente. Sin embargo, también debe considerarse la posibilidad de fallo por pandeo local de los elementos que conforman la viga, es lo que se conoce como fenómeno de abolladura[19,22,23,26].

Las secciones transversales de los elementos a flexión buscan maximizar el valor de su módulo resistente para reducir el valor de las tensiones y desplazamientos máximos. Por estos motivos son muy adecuadas y se emplean con frecuencia las

* Autor para correspondencia. Correos electrónicos: cacho@eii.uva.es (M. Cacho-Pérez), ali@eii.uva.es (A. Lorenzana Ibán).

secciones tipo doble-T u otras de pared delgada, bien sean abiertas o cerradas[20,21].

Por lo tanto, resulta de gran interés el estudio de placas planas cargadas en su plano, pues es el caso del alma y de las alas de vigas y pilares de este tipo de perfil y ya que la mayoría de los elementos estructurales se pueden considerar formados por un conjunto de placas planas. De esta manera, cuando una placa se encuentra sometida a esfuerzos de compresión, flexión, cortante o combinación de ellos, puede que se produzca el fallo local por pandeo antes de cualquier otro tipo de fallo.

Timoshenko (1921) fue el primero en presentar una solución para el pandeo local de placas rectangulares simplemente apoyadas en sus cuatro bordes sometidas a cortante. Lo hizo aplicando el método energético, el cual resultó ser una excelente herramienta para resolver un problema que no podía ser resuelto analíticamente como problema de autovalores. En anos posteriores diversos investigadores como Southwell y Skan, Seydel, Stein y Neff, Galam-bos y otros, continuaron con esas investigaciones para obtener

http://dx.doi.org/10.1016/j.rimni.2016.01.001

0213-1315/© 2016 CIMNE (Universitat Politécnica de Catalunya). Publicado por Elsevier España, S.L.U. Este es un artículo Open Access bajo la licencia CC BY-NC-ND (http://creativecommons.org/licenses/by-nc-nd/4.0/).

ll'llHWi^M lili IIIII.E IN PRESS

M. Cacho-Pérez, A. Lorenzana lbán / Rev. int. métodos numér. cálc. diseño ing. 2016;xxx(xx):xxx-xxx

tensiones críticas de abolladura más precisas y teniendo en cuenta las diferentes condiciones de contorno de los paneles.

Todas las expresiones deducidas por estos investigadores se fundamentan en que existe una perfecta proporcionalidad entre tensiones y deformaciones de acuerdo a la ley de Hooke generalizada. Sin embargo, estas condiciones se sobrepasan en situaciones de diseno usuales de las estructuras metálicas.

Al contrario del caso del estudio del pandeo de la pieza prismática aislada sometida a compresión axil, no existe una teoría general que estudie el comportamiento de las placas en estas condiciones proporcionando unas relaciones de fácil manejo para el diseno, por ello en lugar de usar teorías más rigurosas como la de Stowell (1948), que utiliza un diagrama tensión-deformación curvo con resultados que se ajustan bien a los de los ensayos, se prefiere en general emplear teorías simplificadas como la de Bleich (1952).

Mientras que en la columna con carga axial la diferencia entre la carga de pandeo y su resistencia posterior es tan pequena que la carga de pandeo se puede considerar su límite resistente, en el caso de las placas la resistencia post-pandeo puede ser mucho mayor que su carga precrítica. La placa presenta, a diferencia del caso de la columna aislada, un comportamiento totalmente distinto al alcanzar su carga de pandeo, ya que con el aumento de sus deformaciones y desplazamientos modifica su comportamiento mecánico al aparecer en ella el efecto membrana. Por este efecto forma se produce una redistribución de los esfuerzos con la aparición de las fuerzas estabilizadoras de membrana, permitiendo a la placa recuperar su estabilidad en una nueva configuración deformada. Si se siguen aumentando las cargas, la placa en su nueva situación plastifica en las zonas donde su estado tensional es máximo hasta alcanzar su resistencia última[24,25j.

En cuanto a la influencia de los rigidizadores longitudinales, diversos han sido los autores que han propuesto soluciones para la variación del coeficiente de abolladura en función del tipo y número de rigidizadores empleados (Crate y Lo (1948), Hoglund (1977) o Beg (2003)). Finalmente, autores como Lee et al. (1996) o Estrada (2005), han estudiado también factores como la influencia del grado de empotramiento que proporcionan las alas en función de los espesores de alas y almas, proponiendo expresiones que toman como base a los resultados clásicos citados anteriormente[18].

Por lo tanto, queda claro que los problemas asociados al fenómeno de pandeo en elementos estructurales de acero son a veces la causa más importante de fallo mecánico[1,3,4,7,9], sobre todo hoy día, que el material cada vez es más resistente y las estructuras proyectadas cada vez más esbeltas.

La formulación matricial del Método Directo de Rigidez (MDR) de estructuras de barras para el caso de no linealidad geométrica en hipótesis de linealidad material facilita la estimación de la carga crítica de pandeo global teórica de forma sistemática con la ayuda del ordenador[10-13]. Sin embargo, lógicamente no permite ningún tipo de comprobación sobre la posible inestabilidad local de los elementos que forman los perfiles de tipo doble-T, las alas y el alma, de las estructuras analizadas. Entre otros motivos, porque la tipología estructural necesaria para realizar dicho análisis es otra, tipo placa en lugar de tipo barra[6,15-17].

Por este motivo, es muy interesante y es el objetivo de este trabajo llevar a cabo dicha comprobación a partir de los resultados disponibles tras un análisis matricial del sistema estructural de barras, desplazamientos y esfuerzos a nivel de sección. Y mediante la teoría de Resistencia de Materiales pasar dichos resultados como condiciones de contorno al problema de placas, con el fin de poder estimar el valor del factor crítico de carga que provoca el fenómeno de abolladura[2,5,8,14].

Este trabajo se ha organizado de la manera siguiente: en primer lugar, tras esta introducción, se presenta la metodología empleada. A continuación, se muestran varios ejemplos de aplicación y en el

Figura 1. Perfil laminado serie IPE.

último apartado se presentan las principales conclusiones obtenidas.

2. Metodología

Se pretende realizar la comprobación a pandeo local de los elementos más comprometidos de cualquier sistema espacial de barras de perfil tipo doble-T. Se asume comportamiento elástico-lineal del material y se adopta un estado proporcional de carga, siendo \ dicho factor de proporcionalidad. Para alcanzar dicho objetivo se establece la siguiente metodología:

1) Análisis lineal de la estructura. La estructura se modela con elementos tipo barra esbelta y se resuelve mediante MDR. Así, se obtienen los grados de libertad de la estructura, desplazamientos y giros desconocidos en los nudos para el valor nominal de las cargas (k = 1).

2) Análisis de estabilidad. Se realiza un análisis no lineal del modelo de barras para calcular el factor de carga crítico de pandeo global (kcr1) y se determina la correspondiente carga crítica asociada a los posibles modos de pandeo del elemento tipo barra: pandeo por flexión, pandeo por flexión y/o torsión, pandeo lateral o vuelco, etc.

3) Se pasa del modelo de barras al modelo de placas para analizar el elemento más solicitado de la estructura frente a posible abolladura. Para lo cual se trasladan ciertos resultados del modelo de barras obtenidos en el primer análisis, como condiciones de contorno en desplazamientos del modelo de placas.

4) Análisis para obtener los esfuerzos en el modelo de placas.

5) Análisis no lineal del modelo de placas, para estimar las cargas críticas asociadas al pandeo local (kcr2) del elemento más solicitado.

3. Resultados y discusión

A continuación se muestran varios ejemplos a los que se aplica la metodología de este trabajo. Los dos primeros ejemplos a modo de validación del método y un tercer ejemplo de aplicación.

Para todos los cálculos los datos son: longitud (L = 4 m), sección tipo doble-T

(h = 300 mm, b = 150 mm, e1 = 10.7 mm, e = 7.1 mm), corresponde al perfil comercial 1PE300 (véase fig. 1), de acero (E = 2.1 -1011 Pa, v = 0.3, oy = 275MPa).

''ИгЯИ^И RIIILE IN PRESS

M. Cacho-Pérez, A. Lorenzana Ibán | Rev. int. métodos numér. cálc. diseño ing. 2016;xxx(xx):xxx-xxx

Figura 2. Abolladura del ala.

3.1. Ejemplo 1. Abolladura del ala comprimida

El primer ejemplo objeto de estudio es un pilar empotrado-libre sometido a un momento concentrado de valor M0 = 105 N • m en su extremo libre. Para el valor nominal de las cargas (X = 1) se lleva a cabo un análisis lineal de la estructura mediante un modelo de barras 3D (se emplean 16 elementos) y se obtiene entre otros resultados un desplazamiento máximo de valor SMAX = 0.0456 m. A continuación se realiza un análisis de estabilidad global del modelo de barras, que considera flexión según ambos planos y torsión no uniforme, y tras el análisis del equilibrio en la configuración deformada se obtiene un valor crítico del factor de carga de Xcr1 =0.63, que corresponde al nivel de carga que provoca el pandeo global o vuelco de la estructura. Del análisis lineal se conocen los grados de libertad de toda la estructura, por lo tanto se malla el pilar con elementos tipo placa (4x4,16 elementos placa por pared del perfil, dos alas y el alma, en total 48 elementos placa) y se pasan los resultados de desplazamientos y giros del modelo de barra 3D al modelo placa, con objeto de determinar los esfuerzos de tipo placa. Por último, se hace un análisis no lineal del modelo placa para estimar a qué nivel de carga se pueden producir inestabilidades locales y se comprueba que para un factor de carga de Xcr2 = 11.45 se produce la abolladura el ala comprimida tal y como puede verse en la fig. 2.

A modo de comprobación, se ha resuelto el problema con el programa comercial de elementos finitos Simulation Mechanical 2015 de Autodesk (abreviado ASM), para ello se ha mallado el pilar completo con elementos placa (plate del software indicado), se ha empleado un total de 1218 elementos y se ha obtenido un valor

crítico del factor de carga de X^^M? = 11.45 para el fenómeno de inestabilidad por abolladura del ala.

3.2. Ejemplo 2. Abolladura del alma por cortadura

El segundo ejemplo corresponde a un pilar empotrado-libre sometido a una carga transversal puntual de valor P0 = 105 N en su extremo libre. Del análisis lineal se obtiene un desplazamiento máximo de valor SMAX = 0.122 m (se emplean 16 elementos) y del análisis de estabilidad global se obtiene un valor crítico del factor de carga de Xcr1 = 1.40, que corresponde al nivel de carga que provoca el pandeo con deformaciones de flexión y/o torsión. Ahora, se malla el pilar con elementos tipo placa (16x4, 16 elementos placa en dirección longitudinal por 4 en dirección transversal y por pared del perfil, en total 192 elementos placa) y se lleva a cabo un análisis lineal del modelo placa para determinar los esfuerzos. A continuación, tras el análisis de estabilidad, se estima que el nivel de carga que puede producir inestabilidad local es de Xcr2 = 3.41 y corresponde a la abolladura del alma (véase la fig. 3).

De nuevo, y a modo de comprobación, se ha resuelto el problema con el programa ASM (se ha empleado un total de 1218 elementos) y se ha obtenido un valor crítico del factor de carga de X^?™ = 3.42 para el fenómeno de inestabilidad por abolladura del alma.

3.3. Ejemplo 3. Pórtico espacial

Como ejemplo de aplicación se resuelve una estructura espacial de barras tipo pórtico 3D (fig. 4). Las cargas aplicadas se suponen concentradas de valor (P). Se asume que todas las barras tienen la

'iMlUMi^M AI IIIII.E IN PRESS

M. Cacho-Pérez, A. Lorenzana lbán / Rev. int. métodos numér. cálc. diseño ing. 2016;xxx(xx):xxx-xxx

Figura 4. Pórtico 3D.

misma longitud (L), las bases de los pilares están perfectamente empotradas y todos los nudos son rígidos. Todas las barras son de perfil tipo IPE300 y su disposición espacial es la que se indica.

Para resolver este problema lo primero que se lleva a cabo es un análisis lineal de la estructura, modelo de barras (se emplean 16 elementos/barra), mediante MDR y se obtienen los grados de libertad de la estructura, desplazamientos y giros desconocidos en los nudos. Para el valor nominal de carga (P0 = 105 N; k = 1) se obtiene un desplazamiento máximo de valor SMAX = 0.0219 m.

En segundo lugar, se realiza un análisis de inestabilidad del modelo de barras para calcular el factor de carga crítico de pandeo global kcr1 = 2.08 y determinar la correspondiente carga crítica.

A continuación, se pasa del modelo de barras al modelo de placas para analizar el pilar más solicitado. Para ello, se trasladan los resultados del modelo de barras obtenidos en el primer análisis como condiciones de contorno en desplazamientos del modelo de placa. Con el modelo placa se hace un análisis lineal para obtener los esfuerzos de placa y después un análisis de inestabilidad para

estimar las cargas críticas asociadas al pandeo local. Se emplea una malla de 16x4 elementos placa por pared del pilar analizado (las dos alas y el alma), por lo tanto se emplean 16 elementos placa en dirección longitudinal por 4 en dirección transversal y por pared del perfil, en total 192 elementos placa. Se estima un valor crítico del factor de carga de kcr2 = 10.78 para el fenómeno de abolladura del alma del pilar analizado.

También para este ejemplo de aplicación se ha querido contrastar los resultados con los obtenidos con software ASM de elementos finitos. Para ello se ha mallado toda la estructura con elementos placa (se emplean un total de 9476 elementos) y el programa estima que la abolladura del alma del pilar más solicitado (véase la fig. 5) se produce para un valor del factor de carga de k^M = 10.83, valor muy próximo al estimado con la metodología de este trabajo.

Por lo tanto, para el pórtico 3D analizado se comprueba que la metodología de este trabajo permite contar con una herramienta útil para la comprobación del fenómeno de inestabilidad local que corresponde a la posibilidad de abolladura del alma del pilar más solicitado.

4. Conclusiones

La principal conclusión consiste en que este trabajo permite analizar un fenómeno de inestabilidad tan complejo como es la abolladura de este tipo de perfiles mediante una metodología sencilla, con mínimos requerimientos computacionales disponibles en cualquier oficina técnica.

Por lo tanto, es suficiente discretizar la estructura mediante elementos barra 1 D y emplear un sencillo método matricial (MDR) para resolver la estructura y posteriormente, en la siguiente etapa del análisis, calcular y comprobar la inestabilidad local de los elementos estructurales aisaldos mediante la metodología expuesta en este trabajo.

Anexo. Formulación de placa delgada

Para el análisis del fenómeno de abolladura, se emplea la teoría de placas delgadas (de Kirchhoff) con posibilidad de grandes desplazamientos [17]. La hipótesis fundamental de dicha teoría consiste en despreciar los efectos de las deformaciones transversales en base al pequeno espesor de la placa.

En base a esta hipótesis, se pueden expresar los desplazamientos de cualquier punto material de la placa como sigue:

■Л /0Л

v - z 0y (1)

De tal forma que las deformaciones de Green-Lagrange no nulas del modelo de placas son:

1 í \ ' + 2 (wy)

y + 2 V 'У

\u,y + v,x + w,xw,y )

í w,;

= Ep-z-Kb

Figura 5. Pórtico 3D. Abolladura del alma.

(asociado al valor crítico del factor de carga kcr,2 = 10.83).

donde EP son las deformaciones de membrana en el plano y Kb representa los cambios de curvatura por flexión.

'MllPPI^M RlIiLE IN PRESS

M. Cacho-Pérez, A. Lorenzana Ibán / Rev. int. métodos numér. cálc. diseño ing. 2016;xxx(xx):xxx-xxx

Figura 6. Elemento placa normalizado (4 nodos, 5 grados de libertad/nodo).

(Tx Ty

\Txy, (Mx My

a ■ dz =

yy= \Mxy,

donde t es el espesor de la placa. Referencias

ayy dz

J- /2 ayy axy

Matriz de rigidez tangente

La expresión de interés de la matriz de rigidez tangente se indica a continuación (fig. 6):

; <x,p = 1.. .4 (3)

donde cada submatriz tiene por expresión:

¡■1 !• 1

:<M)a/3 ■

KM )ap (KM)

(Kt)afi = 1 (KM?ap (KM)

(KM)afi = abJ J (BPa)TDpB^dV; = 1...4 (KM)afi = abJ j (Bba)TDbBfádm = 1...4

(KM)afi = ab j J (BPa)TD%dm a,P = 1...4

(KLc)afi =ab

Tx Txy

-1 -1 y se define como:

00 (1-v

Gpd^dv; a,fi = 1...4

Db = — Dp 12

\Na¿/

Ga —

NW N -

Na,x Na,x Na,x

nw N°x N9y ' Na,y Na,y Na,y /

Nwt 1 Ne\ 1 №yt a a a¿ a

1 Nw 1N0X 1 №

\ bNa,n b' b' /

Es importante indicar que el término BL supone un acoplamiento entre el comportamiento de membrana y de flexión de la placa, debido a la no linealidad geométrica del problema planteado.

Por último, se definen los esfuerzos de placa, fuerzas de membrana (Tp) y momentos flectores (Mb):

[1 [2 [3 [4

[5 [6 [7 [8

[10 [11

[12 [13

[14 [15 [16

[20 [21

[23 [24

M. Cacho-Pérez, A. Lorenzana, Cálculo directo de la carga crítica de pandeo de pórticos. Parte II, Rev. Int. Mét. Num. Cálc. Dise. Ing. 26 (2010) 31-38. W.F. Chen, T. Atsuta, Theory of beam-columns. Volume 1: in-plane behaviour and design, J. Ross Pub Inc (2007).

A. Lorenzana, M. Cacho-Pérez, Cálculo directo de la carga crítica de pandeo de pórticos, Parte I. Rev. Int. Mét. Num. Cálc. Dise. Ing. 25 (3) (2009) 247-258. M.A. Ortega, J.L. Romero, Rosa.F E., A historical study of the problem of straight prismatic elements subjected to compression, Informes de la Construcción 59 (507) (1997) 69-81.

M. Cacho-Pérez, A. Lorenzana, Optimización global de pórticos 2 D con barras de sección variable, DYNA 85 (8) (2010) 667-675.

W.F. Chen, T. Atsuta, Theory of beam-columns. Volume 2: space behaviour and design, J. Ross Pub Inc (2007).

Z. Lei,T.G. Shu, Lateral buckling ofweb-tapered I-beams: A new theory, J. Constr. Steel Res. 64 (12) (2008) 1379-1393.

M. Kovac, Lateral-torsional buckling of web-tapered I-beams. 1D and 3D FEM approach, 23rd Czech and Slovak Conference on Steel Structures and Bridges 40(2012)217-222.

S.E. Kim, C.M. Uang, S.H. Choi, K.Y., An. Practical advanced analysis of steel frames considerin lateral-torsional buckling, Thin-Walled Structures 44 (7) (2006) 709-720.

N. Hoveidae, B. Rafezy, Global buckling prevention condition of all-steel buckling restrained braces, J. Theor. App. Mech. 51 (4) (2013) 891-892. M.A. Serna, A. Lopez, I. Puente, Yong.F D.J., Equivalent uniform moment factors for lateral-torsional buckling ofsteel members, J. Constr. Steel Res. 62 (6) (2006) 566-580.

J.D. Aristizabal-Ochoa, Column stability and minimum lateral bracing: Effects of shear deformations, J. Eng. Mech. 130 (10) (2004) 1223-1232.

A. Andrade, D. Camotim, P.B. Dinis, Lateral-torsional buckling of singly symmetric web-tapered thin-walled I-beams: 1 D model vs. shell FEA, Comput. Struct. 85 (17-18) (2007) 1343-1359.

W.F. Chen, E.M. Lui, Stability design ofsteel frames, CRC Press, Inc. (2000). W. McGuire, R.H. Gallagher, Matrix Structural Analysis, Wiley, New York, 1979. H. Kardestuncer, Introducción al análisis estructural con matrices, McGraw-Hill, 1975.

O.C. Zienkiewicz, R.L. Taylor, El método de los elementos finitos. Volumen 2: mecánica de sólidos, 5a ed., CIMNE, Barcelona, 2004.

N. Guerrero, M.E. Marante, R. Picon, J.F. Lopez, Analysis of steel hollow structural beams subjected to biaxial bending by lumped damage mechanics, Rev. Int. Mét. Num. Cálc. Dise. Ing. 25 (1) (2009)3-27.

C. Basaglia, D. Camotim, Buckling Analysis of Thin-Walled Steel Structural Systems Using Generalized Beam Theory (GBT), International Journal of Structural Stability and Dynamics 15 (1) (2015) 1540004.

D. Camotim, C. Basaglia, On the behaviour, failure and direct strength design of thin-walled steel structural systems, Thin-Walled Structures 81 (2014) 50-66. C.K. Iu, M.A. Bradford, Second-orderelastic finite element analysis ofsteel structures using a single element per member, Engineering Structures 32(9) (2010) 2606-2616.

J. Becque, K.J.R. Rasmussen, Experimental Investigation of the Interaction of Local and Overall Buckling of Stainless Steel I-Columns, Journal of Structural Engineering-ASCE 135 (11) (2009) 1340-1348.

M. Lecce, K.J.R. Rasmussen, Design of wide-flange stainless steel sections, Advanced Steel Construction 5 (2) (2009) 164-174.

P.B. Dinis, D. Camotim, N. Silvestre, FEM-based analysis of the local-plate/distortional mode interaction in cold-formed steel lipped channel columns, Computers & Structures 85 (19-20) (2007) 1461-1474. K. Kheradmandia, M. Azhari, M.A. Bradford, Local instability and free vibration of supported plates using the basic plate triangle, Computational Mechanics 40 (1) (2007) 73-84.

B. Young, J.T. Tan Channel columns undergoing local, distortional, and overall buckling, Journal ofStructural Engineering-ASCE 128 (6) (2002) 728-736.

Txy Ty