Alabeo seccional

El alabeo unitario o alabeo seccional es una función ω(y,z) que predice la forma deformada de la sección transversal de un prisma mecánico y que define varias características geométricas importantes relacionadas con el cálculo de tensiones en caso de flexión, torsión y cortante combinados. [Este alabeo unitario tiene dimensiones de longitud al cuadrado (L²)].

Para un prisma mecánico de sección constante A, el alabeo unitario es una función ${displaystyle omega (y,z);}$ definida sobre dicha sección transversal, que es solución del siguiente problema de Von Neumann:

(1)${displaystyle {egin{cases}{cfrac {partial ^{2}omega }{partial y^{2}}}+{cfrac {partial ^{2}omega }{partial z^{2}}}=0&forall (y,z)in A\mathbf {n} cdot {oldsymbol { abla }}omega ={cfrac {1}{2}}{cfrac {d}{d{ar {s}}}}left[(y-y_{C})^{2}+(z-z_{C})^{2} ight]&forall (y,z)in partial Aend{cases}}}$

Donde:

En el problema de torsión pura de Saint-Venant para una pieza prismática la hipótesis cinemática lleva a que los desplazamientos están relacionados con los giros del eje baricéntrico alrededor de sí mismo por la siguiente condición:

(2)${displaystyle {egin{Bmatrix}u_{x}(x,y,z)\u_{y}(x,y,z)\u_{z}(x,y,z)end{Bmatrix}}={egin{bmatrix}omega (y,z)&0\0&-z+z_{C}\0&y-y_{C}end{bmatrix}}{egin{Bmatrix}{cfrac {d heta _{x}}{ds}}\ heta _{x}end{Bmatrix}}}$

Calculando a partir de ellos las deformaciones y aplicando después las ecuaciones de Lamé-Hooke se llega a que la relación entre tensiones y giros sobre el eje son:

(3)${displaystyle {egin{cases}sigma _{xx}={frac {2G}{(1-2 u )}}left[(1- u )varepsilon _{xx}+ u (varepsilon _{yy}+varepsilon _{zz}) ight]=0\sigma _{xy}=Gleft[{partial u_{x} over partial y}+{partial u_{y} over partial x} ight]=Gleft[{frac {partial omega }{partial y}}-left(z-z_{C} ight) ight]{frac {d heta _{x}}{ds}}\sigma _{xz}=Gleft[{partial u_{x} over partial z}+{partial u_{z} over partial x} ight]=Gleft[{frac {partial omega }{partial z}}+left(y-y_{C} ight) ight]{frac {d heta _{x}}{ds}}end{cases}}}$

El equilibrio de fuerzas sobre el eje longitudinal de la pieza prismática o viga requiere que:

(4)${displaystyle {frac {partial sigma _{xx}}{partial x}}+{frac {partial sigma _{xy}}{partial y}}+{frac {partial sigma _{xz}}{partial z}}=0}$

Donde substiuyendo las ecuaciones (3) en (4) se llega precisamente a la ecuación del alabeo unitario (1).

Puede demostrarse que la solución de la anterior ecuación puede encontrarse fácilmente introduciendo una función de alabeo auxiliar relacionada con la anterior y con las coordenadas (y_C, z_C) del centro de cortante. La función auxiliar ${displaystyle omega _{0}(y,z);}$ satisface la ecuación:^[1]​

${displaystyle {egin{cases}{cfrac {partial ^{2}omega _{0}}{partial y^{2}}}+{cfrac {partial ^{2}omega _{0}}{partial z^{2}}}=0&forall (y,z)in A\mathbf {n} cdot {oldsymbol { abla }}omega _{0}={cfrac {1}{2}}{cfrac {d}{d{ar {s}}}}left(y^{2}+z^{2} ight)&forall (y,z)in partial Aend{cases}}}$

En términos de esta función auxiliar se pueden encontrar tanto la función de alabeo como las coordenadas del centro de cortante:^[1]​

${displaystyle omega (y,z)=omega _{0}(y,z)-z_{C}y+y_{C}zqquad y_{C}=-{frac {I_{z}I_{y{ar {omega }}}-I_{yz}I_{z{ar {omega }}}}{I_{z}I_{y}-I_{yz}^{2}}}qquad z_{C}={frac {I_{y}I_{z{ar {omega }}}-I_{yz}I_{y{ar {omega }}}}{I_{z}I_{y}-I_{yz}^{2}}}}$

Donde ${displaystyle I_{y},I_{z},I_{yz},}$ son los momentos de área y el producto de inercia. Y donde ${displaystyle I_{y{ar {omega }}},I_{z{ar {omega }}},}$ son los productos de inercia sectoriales definidos como:

${displaystyle I_{y{ar {omega }}}=int _{A}zomega _{0}(y,z) dydzqquad I_{z{ar {omega }}}=int _{A}yomega _{0}(y,z) dydz}$

En general si una sección no es circular o circular hueca presentará alabeo seccional diferente de cero. Esto puede probarse rigurosamente calculando el alabeo seccional de una sección elíptica, que depende de la diferencia de cuadrados de las longitudes de los semiejes, si estos son iguales como sucede en un círculo la función de alabeo se anula.

En el caso general la sección de alabeo es complicada y requiere resolver un problema de Von Neumann. Para algunos casos sencillos cuando la sección es maciza y el contorno viene expresado por una función de tipo f(y, z) = 0 siendo el laplaciano de f constante el problema de buscar la función de alabeo puede simplificarse notablemente mediante la función de Prandtl, ya que esta función basta encontrar una función de Prandtl que se anule sobre el contorno. Esto es precisamente lo que sucede con la secciones elíptica y triangular, sin embargo con secciones más complicadas como una sección rectangular el cálculo es más complicado.

En una sección triangular equilátera cualquiera de las tres alturas del triángulo constituye un eje de simetría, por lo que para una sección triangular equilátera el centro de cortante coincide con el centro geométrico o baricentro del triángulo. La función de alabeo considerando coordenadas (y, z) con el origen de coordenadas sobre el centro geométrico viene dada por:^[2]​

${displaystyle omega (y,z)=-{frac {3z^{2}y-y^{3}}{2h}}}$

Donde hemos considerado que uno de los lados es paralelo al eje Y, y h es la altura del triángulo.

En una sección elíptica existen dos ejes de simetría, el semieje mayor y el semieje menor, lo cual implica que el centro de cortante coincida con el centro geométrico de la sección. Tomando coordenadas de la sección (y, z) con origen en el centro geométrico de la sección la función de alabeo unitario viene dada por:^[3]​

${displaystyle omega (y,z)=-{frac {a^{2}-b^{2}}{a^{2}+b^{2}}}yz}$

Donde a y b son, respectivamente, las longitudes del semieje mayor y el semieje menor de la elipse. Puede verse que en el caso particular de un círculo de radio r (donde a = b = r) el alabeo seccional unitario es nulo, en consonancia con la teoría de la torsión de Saint-Venant para secciones circulares.

En una sección rectangular, donde el centro de cortante coincide con centro geométrico, la función de alabeo puede calcularse en términos de la función de Prandtl^[4]​ que a su vez puede obtenerse por integración de Laplace mediante separación de variables:

${displaystyle {egin{cases}{cfrac {partial omega }{partial y}}=z+{cfrac {1}{G heta }}{cfrac {partial Phi }{partial z}}=z+{cfrac {8b}{pi ^{2}}}sum _{k=0}^{infty }left[{cfrac {(-1)^{k+1}}{(2k+1)^{2}}}left({cfrac {sinh {frac {(2k+1)pi z}{b}}}{cosh {frac {(2k+1)pi a}{b}}}} ight)cos {frac {(2k+1)pi y}{b}} ight]\{cfrac {partial omega }{partial z}}=-y-{cfrac {1}{G heta }}{cfrac {partial Phi }{partial y}}=-y+{cfrac {8b}{pi ^{2}}}sum _{k=0}^{infty }left[{cfrac {(-1)^{k}}{(2k+1)^{2}}}left(1-{cfrac {cosh {frac {(2k+1)pi z}{b}}}{cosh {frac {(2k+1)pi a}{b}}}} ight)sin {frac {(2k+1)pi y}{b}} ight]end{cases}}}$

El momento de alabeo es la magnitud definida por la siguiente integral:^[5]​

${displaystyle I_{omega }=int _{A}omega ^{2}(y,z) dydz}$

Para una sección I o H el módulo de alabeo viene dado por:^[6]​

${displaystyle I_{omega }={frac {h^{2}I_{min}}{4}}}$

Donde h denota la altura total del perfil e I_min el momento de inercia mínimo.

---