Pandeo

El pandeo es un fenómeno llamado inestabilidad elástica que puede darse en elementos comprimidos esbeltos, y, que se manifiesta por la aparición de desplazamientos importantes transversales a la dirección principal de compresión.

En ingeniería estructural el fenómeno aparece principalmente en pilares y columnas, y se traduce en la aparición de una flexión adicional en el pilar cuando se halla sometido a la acción de esfuerzos axiales de cierta importancia.

La aparición de deflexión por pandeo limita severamente la resistencia en compresión de un pilar o cualquier tipo de pieza esbelta. Eventualmente, a partir de cierto valor de la carga axial de compresión, denominada carga crítica de pandeo, puede producirse una situación de inestabilidad elástica y entonces fácilmente la deformación aumentará produciendo tensiones adicionales que superarán la tensión de rotura, provocando la ruina del elemento estructural. Además del pandeo flexional ordinario existe el pandeo torsional o inestabilidad elástica provocado por un momento torsor excesivo.

Existen diferentes maneras o modos de fallo por pandeo. Para un elemento estructural frecuentemente hay que verificar varios de ellos y garantizar que las cargas están lejos de las cargas críticas asociadas a cada modo o manera de pandear. Los modos típicos son:

Los pilares y barras comprimidas de celosías pueden presentar diversos modos de fallo en función de su esbeltez mecánica:

El pandeo local es el que aparece en piezas o elementos aislados o que estructuralmente pueden considerarse aislados. En este caso la magnitud de la carga crítica viene dada según el caso por la fórmula de Leonhard Euler o la de Engesser. La carga crítica de Euler depende de la longitud de la pieza, del material, de su sección transversal y de las condiciones de unión, vinculación o sujeción en los extremos. Para una pieza que puede considerarse biarticulada en sus extremos la carga crítica de Euler viene dada por:

(1)${displaystyle F_{ m {crit}}=pi ^{2}{frac {EI_{min }}{L^{2}}}=pi ^{2}{frac {EA}{lambda ^{2}}}}$

Siendo: F_crit, la carga crítica; E, Módulo de Young del material de que está hecha la barra; I_min, momento de inercia mínimo de la sección transversal de la barra; L, longitud de la barra y λ la esbeltez mecánica de la pieza. Cuando las condiciones de sujeción de los extremos son diferentes la carga crítica de Euler viene dada por una ecuación del tipo:

(2)${displaystyle F_{ m {crit}}=pi ^{2}{frac {EI_{min }}{(alpha L)^{2}}}qquad {egin{cases}alpha =0,5&{mbox{empotrado-empotrado}}\alpha approx 0,70&{mbox{articulado-empotrado}}\alpha =1&{mbox{articulado-articulado}}\alpha =2&{mbox{libre-empotrado}}end{cases}}}$

Al producto ${displaystyle alpha L;}$ se le llama 'longitud de pandeo'.

En una estructura compleja formada por barras y otros elementos enlazados pueden aparecer modos de deformación en los que los desplazamientos no sean proporcionales a las cargas y la estructura puede pandear globalmente sin que ninguna de las barras o elementos estructurales alcance su propia carga de pandeo. Debido a este factor, la carga crítica global de cierto tipo de estructuras (por ejemplo en entramados de cúpulas monocapa) es mucho menor que la carga crítica (local) de cada uno de sus elementos.

El tipo de estructura más simple que presenta 'pandeo global' para carga crítica diferente de la de sus elementos está formado por dos barras articuladas entre sí^[1]​ y a la cimentación, que se muestra en la figura.

Las ecuaciones que gobiernan el comportamiento de la estructura son:

${displaystyle P=2Nsin( heta -Delta heta ),}$

${displaystyle N=EA{frac {Delta L}{L}}}$

${displaystyle Lcos heta =(L-Delta L)cos( heta -Delta heta ),}$

Donde: N, esfuerzo axial de cada una de las barras; ΔL, acortamiento sufrido por las barras para adoptar la configuración deformada; Δθ = θ-θ', es la diferencia de ángulos mostrada en la figura; E, módulo de Young del material de las barras; A, área transversal de cada una de las barras; L, longitud inicial de cada una de las dos barras.

Substituyendo la segunda de las ecuaciones en la primera, despejando ΔL de la tercera y substituyendo su valor también su valor en la primera se llega a:

${displaystyle P=2EAleft(1-{frac {cos( heta )}{cos( heta -Delta heta )}} ight)sin( heta -Delta heta )}$

El valor de Δθ para el que se alcanza el máximo es precisamente la carga crítica global. Las cargas de pandeo global y local vienen dadas por:

${displaystyle P_{ m {{crit},G}}=2EAleft(1-cos ^{frac {2}{3}} heta ight)qquad qquad P_{ m {{crit},L}}={frac {pi ^{2}EI}{L}}}$

Cada una estas cargas presenta modos de fallo diferentes en la estructura. De entre los dos posibles modos de fallo por pandeo ocurrirá el que presente un ángulo de aparición mayor donde estos ángulos vienen dados por:

${displaystyle cos( heta _{G})=cos ^{frac {1}{3}}( heta )qquad qquad cos( heta _{L})={frac {pi ^{2}EI}{L}}}$

El 'plano de pandeo' se refiere al plano que contendrá el inicio de la deformada de una pieza sometida a compresión dominante. El plano de pandeo contiene el eje baricéntrico de la viga y sobre él la deflexión por pandeo es máxima. Para una pieza sometida solo a compresión sin momentos apreciables adicionales, el plano de pandeo coincidirá con el plano perpendicular sea paralela al eje de menor inercia de la sección.

Matemáticamente el pandeo local está asociado a una bifurcación tridente, es decir, cuando se plantean las ecuaciones exactas no lineales que describen la forma de una pieza prismática, incluyendo la carga axial y los parámetros relacionados con las imperfecciones, los posibles comportamientos cualitativos están separados unos de otros por una bifurcación tridente. Eso lleva que en estos casos la carga real que puede soportar una barra venga dada por la ley 2/3 de Koiter:^[2]​

${displaystyle P_{c}leq C(d)varepsilon ^{2/3}}$

Donde:

Si las imperfecciones tienen naturaleza estadística y vienen dadas por una distribución normal multivariante entonces la carga crítica tiene una distribución dada por:^[3]​

${displaystyle f(P_{c})={frac {3|P_{c}-P_{c}^{0}|^{1/2}}{{sqrt {2pi }}{hat {C}}^{3/2}}}exp left(-{frac {1}{2}}{frac {|P_{c}-P_{c}^{0}|^{3}}{hat {C}}} ight)}$

Esta distribución de probabilidad permite ajustar las curvas reales de pandeo observadas, ya que en condiciones normales una barra recta de sección constante tiene una resistencia inferior debido a paredicha por la teoría de Euler por el efecto de las imperfecciones.

En vigas de alas anchas o de escasa rigidez torsional, el pandeo flexional convencional puede ir acompañado de la aparición de una torsión de la sección, resultando un modo de fallo mixto conocido como pandeo torsional o pandeo lateral. El momento torsor crítico para el cual aparecería ese tipo de fallo viene dado por:^[4]​

${displaystyle M_{ m {cr}}={frac {pi }{alpha L}}{sqrt {left({frac {pi EI_{m}}{alpha L}} ight)^{2}{frac {I_{omega }}{I_{m}}}+EI_{m}GJ}}}$

Donde las nuevas magnitudes son:

Y el resto de magnitudes tienen el mismo significado que para el pandeo flexional puro. En piezas donde el momento de alabeo es despreciable puede usarse la expresión aproximada:

${displaystyle M_{ m {cr}}approx {frac {pi }{alpha L}}{sqrt {EI_{m}GJ}}}$

Una manera de encontrar la carga crítica de una estructura consiste en presuponer la forma cualitativa en que esta pandeará, parametrizando esa forma cualitativa mediante varios parámetros incógnita. Introduciendo esa forma cualitativa en la ecuación de la curva elástica y buscando que la solución parametrizada satisfaga las condiciones de contorno cualitativas, que normalmente se refieren a desplazamientos y giros de los nudos de las barras de la estructura, se obtienen relaciones entre los parámetros incógnita introducidos. El valor de la carga crítica es precisamente el que hace que dichas relaciones se cumplan.

El método de Euler para barras aisladas es un ejemplo de uso de este método. Por ejemplo para determinar la carga crítica de un pilar empotrado en su base y libre en el extremo tratamos de resolver la ecuación de la curva elástica bajo las siguientes condiciones:

${displaystyle {frac {d^{2}w(h)}{dh^{2}}}+{frac {Pw(h)}{EI}}={frac {Pdelta }{EI}}qquad {egin{cases}w'_{ m {base}}=0\w_{sup }=delta end{cases}}}$

La solución de esa ecuación, en función del parámetro de desplazamiento horizontal del pilar, resulta ser:

${displaystyle w(h)=delta left[1-cos left(h{sqrt {frac {P}{EI}}} ight) ight]}$

La condición de contorno en el extremo superior (donde h = H y w_sup = δ) solo se cumple para ciertos valores de P, que cumplen:

${displaystyle H{sqrt {frac {P_{n}}{EI}}}=(2n+1){frac {pi }{2}}qquad nin mathbb {Z} }$

El menor de estos valores es precisamente el valor aceptado para la carga crítica de Euler de un pilar empotrado en su base y libre en su parte superior:

${displaystyle P_{ m {crit}}=P_{0}={frac {pi ^{2}EI}{(2H)^{2}}}}$

Para estructuras de una cierta complejidad el método anterior resulta de muy difícil aplicación, ya que requiere integrar un número elevado de ecuaciones diferenciales para cada elemento lineal de la estructura.

Un método aproximado consiste en presuponer aproximadamente las deformaciones asociadas al pandeo, que satisfaga las condiciones de contorno en los extremos de las piezas, y en igualar la energía de deformación W_int con el trabajo exterior realizado por la fuerza que produce el fenómeno de pandeo W_ext durante la deformación, W_int = W_ext. Esas dos ecuaciones pueden escribirse en términos el campo de desplazamientos de los momentos flectores asociados. Para cada elemento lineal la energía de deformación y el trabajo exterior vienen dados por:

${displaystyle W_{ m {int}}={frac {1}{2}}int _{0}^{L}{frac {M_{f}^{2}(x)}{EI}}dx={frac {1}{2}}int _{0}^{L}EIleft[{frac {d^{2}w(x)}{dx^{2}}} ight]^{2}dxqquad qquad W_{ext}={frac {P_{ m {crit}}}{2}}int _{0}^{L}left[{frac {dw(x)}{dx}} ight]^{2}dx}$

Donde:

Cuanto más ajustado sea el campo de desplazamientos supuesto w(x) mejores resultados da el método anterior.

En ingeniería estructural existe una necesidad práctica de dimensionar los elementos lineales sometidos a compresión con la suficiente sección transversal como para que no fallen por pandeo. La sección transversal necesaria para que eso no ocurra es muchas veces mayor que la que sería necesaria para soportar un esfuerzo de tracción de la misma magnitud (entre 1,5 y 6 veces en la mayoría de casos). La mayoría de normas usan un coeficiente de reducción de la resistencia cuando el esfuerzo sobre el elemento lineal es de compresión y no de tracción. El Eurocódigo por ejemplo da para la resistencia de un pilar sometido a compresión y tracción simples las siguientes resistencias:

${displaystyle N_{t,u}=Af_{y},qquad N_{t,c}=chi A_{ef}f_{y}}$

Donde:

El mismo coeficiente se puede usar para estimar por exceso la tensión y determinar si un elemento es seguro. Así cuando un elemento está sometido a flexión o compresión compuestas la tensión de referencia para calcular si el elemento es seguro o no se toma aproximadamente como:

${displaystyle sigma _{r}={frac {|N_{x}|}{chi A}}+{frac {|M_{y}|}{W_{y}}}+{frac {|M_{z}|}{W_{z}}}leq sigma _{u}}$

Donde:

En situaciones donde las tensiones tangentes y el alabeo seccional de la sección sean importantes debe substituirse el miembro antes del signo de menor que, por la tensión de Von Mises y en la expresión de Navier debe contabilizarse el efecto del bimomento.

La carga crítica de un elemento estructural unidimensional esbelto corresponde a un esfuerzo axial por encima de la cual cualquier pequeña imperfección impide que exista un equilibrio estable. Para una pieza prismática recta muy esbelta, de material elástico y con extremos articulados, la carga crítica se aproxima mucho a la llamada carga crítica de Euler:

${displaystyle N_{ m {crit}}=N_{ m {Eul}}:={frac {pi ^{2}EI_{f}}{L^{2}}}}$

Donde:

En otros casos más complejos con otras condiciones en los extremos, con sección variable, etc, la carga crítica anterior debe ser corregida por un factor constante. En piezas de sección constante puede definirse además la longitud de pandeo ${displaystyle scriptstyle L_{k}}$ o ${displaystyle scriptstyle L_{pan}}$como:

${displaystyle L_{k}=L_{ m {pan}}:=pi {sqrt {frac {EI_{f}}{N_{ m {cr}}}}}=pi i_{c}{ar {lambda }}{sqrt {frac {E}{f_{y}}}}}$

Donde:

Si la pieza no es de sección constante no existe una manera de definir la longitud de pandeo, aunque el concepto de carga crítica sigue estando perfectamente definido.

En el enfoque moderno de la teoría de bifurcación corresponde a un punto del espacio de configuración tal que cualquier entorno de ese punto se interseca con más de una solución de las ecuaciones de comportamiento estructural. Los elementos bidimensionales comprimidos como los muros de carga, entre otros, también pueden sufrir pandeo, aunque en ese caso la carga crítica se define en términos de la carga compresiva sobre el borde de la misma, para la que aparecen fenómenos de pandeo.

Usualmente las diferentes normas tecnológicas preven una reducción de la resistencia de pilares y otras piezas en términos de su esbeltez mecánica. Cuanto más esbelto sea el elemento tanto mayor será la reducción de su resistencia debida al probable efecto de pandeo sobre el mismo. Existen varias maneras, todas ellas esencialmente equivalentes, de tratar esta reducción de la resistencia por efecto del pandeo, por ejemplo el eurocódigo y el CTE definen la esbeltez mecánica reducida ${displaystyle scriptstyle {ar {lambda }}}$ o razón entre la resistencia plástica de la sección de cálculo y la compresión crítica por pandeo, como:

${displaystyle {ar {lambda }}={sqrt {frac {Af_{y}}{N_{ m {cr}}}}}}$

Donde:

El factor de reducción de la resistencia por pandeo o ${displaystyle scriptstyle chi }$ (coeficiente ji), se de acuerdo con el CTE simplemente como:

${displaystyle {egin{cases}phi =0,5[1+alpha ({ar {lambda }}-0,2)+{ar {lambda }}^{2}]\chi ={cfrac {1}{phi +{sqrt {phi ^{2}-{ar {lambda }}^{2}}}}}end{cases}}}$

Donde en la fórmula anterior:

---