Cálculo de las Constantes de Torsión para Vigas de Secciones Transversales Arbitrarias Utilizando Elementos Finitos



==Cálculo de las Constantes de Torsión para Vigas de Secciones Transversales Arbitrarias Utilizando Elementos Finitos==

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Jesús Gerardo Valdés–Vázquez. Departamento de Ingeniería Civil y Ambiental, Universidad de Guanajuato. Av. Juárez 77, Zona Centro. Guanajuato, Gto., México. [mailto:valdes@ugto.mx valdes@ugto.mx]</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
Jesús Fernando Valdés–Vázquez. Departamento de Ingeniería Civil y Ambiental, Universidad de Guanajuato. Av. Juárez 77, Zona Centro. Guanajuato, Gto., México. valdesjfv[mailto:@ugto.mx @ugto.mx]</div>

===Resumen===

En el análisis estructural, es esencial que un ingeniero cuente con un profundo conocimiento de las propiedades de una sección transversal, tales como el área, los momentos de inercia y las constantes de torsión. En particular, el cálculo de las constantes de torsión presenta un desafío significativo, ya que, para secciones conocidas, existen fórmulas exactas o aproximadas que facilitan el proceso. Sin embargo, con las actuales tendencias en construcción que incorporan secciones compuestas o de geometría arbitraria, el cálculo de estas constantes se convierte en un reto considerable. Este trabajo se centra, principalmente, en explorar la obtención de las constantes de torsión mediante métodos numéricos, mejorando su precisión en comparación con las propuestas de otros investigadores.

===Palabras clave===

Constante de torsión. Torsión de Saint-Venant. Torsión por alabeo.

===Introducción===

El análisis estructural tridimensional de edificaciones, compuesta por vigas y columnas, demanda un profundo entendimiento de las propiedades geométricas inherentes a las secciones transversales de estos elementos. Entre estas propiedades se encuentran el área, los momentos de inercia y las constantes de torsión. Desde la magnificencia de antiguos templos hasta la imponente altura de los contemporáneos rascacielos, el cálculo de la torsión ha representado un desafío constante para los ingenieros a lo largo de la historia.

En los últimos siglos, el estudio de la torsión ha evolucionado, siendo iniciado por distinguidos pensadores. Johann Bernoulli (1667-1748) destacó como precursor, aportando a la comprensión de las fuerzas internas en estructuras y estableciendo las primeras bases para la teoría de vigas. Posteriormente, Daniel Bernoulli (1700-1782), junto con su hermano Johann, realizó contribuciones fundamentales a la teoría de la elasticidad y la mecánica, cimentando así los pilares para el estudio detallado de la torsión. Casi en paralelo, Leonhard Euler (1707-1783) desplegó una contribución esencial al establecer los cimientos matemáticos para la comprensión de fenómenos torsionales. En el siglo XIX, Adhémar Jean Claude Barre de Saint-Venant (1797-1886) formuló las ecuaciones fundamentales para la torsión, enriqueciendo la teoría de la elasticidad y la mecánica de sólidos deformables. En épocas más recientes, Stephen Timoshenko (1878-1972) emergió como una figura clave, realizando contribuciones significativas, especialmente en la teoría de vigas y torsión. En el ámbito práctico de la ingeniería estructural, (Euler, 1744), (Saint-Venant, 1855) y (Timoshenko, 1934) han dejado un legado influyente, proporcionando herramientas y teorías esenciales que han sido fundamentales para el diseño y análisis de estructuras en la ingeniería moderna.

En el análisis espacial de marcos conformados por vigas y columnas, se emplean predominantemente secciones circulares, rectangulares, tubos circulares, tubos cuadrados y secciones en forma de I. Todas estas tipologías cuentan con fórmulas exactas o aproximadas para el cálculo de la constante de torsión de Saint-Venant ( <math display="inline">J</math>), cuyos detalles se pueden encontrar disponibles en (Boresi, Schmidt, & Sidebottom, 1993). Sin embargo, en lo que respecta a la constante de torsión por alabeo o ''warping'' ( <math display="inline">{C}_{w}</math>), calcularla no resulta tan sencillo, llegando al punto en el cual su influencia en el análisis estructural contemporáneo tiende a ser pasada por alto en la mayoría de los softwares comerciales.

Excluyendo la torsión por alabeo, el análisis de estructuras con secciones arbitrarias o compuestas, como en el caso del análisis de tableros de puentes, plantea un desafío adicional. La determinación de la constante de torsión de Saint-Venant se vuelve todo un reto y, en ocasiones, incluso resulta imposible debido a la complejidad geométrica involucrada. En situaciones de esta índole, se recurre a métodos numéricos para abordar el cálculo de la torsión. Estos métodos pueden incluir la aplicación de técnicas como el método de elementos finitos, tal como se detalla en (Pilkey, 2002).

Propuestas más contemporáneas que emplean el método de los elementos finitos para la determinación de la constante de torsión de Saint-Venant ( <math display="inline">J</math>) y la constante de torsión por alabeo ( <math display="inline">{C}_{w}</math>) se describen detalladamente en la investigación llevada a cabo por (Haukaas, 2020). Su estudio se centra específicamente en el cálculo mediante elementos finitos de ambas constantes de torsión utilizando triángulos lineales. Nuestro artículo se basa en este último y presenta una mejora en la determinación de las inercias, las cuales son necesarias como paso previo al cálculo de dichas constantes.

===Conceptos teóricos y formulación===

Comencemos considerando una viga con sección transversal arbitraria que está empotrada en un extremo y sometida a una fuerza de torsión <math display="inline">T</math> en el otro extremo. Según el principio de Saint-Venant, la distribución de esfuerzos en una sección transversal a una distancia <math display="inline">x</math> considerable de los extremos depende principalmente de la fuerza <math display="inline">T</math> y no de la distribución de esfuerzos en dichos extremos. Saint-Venant observó que la deformación por torsión resulta de los desplazamientos de la sección. Siguiendo la formulación de (Boresi, Schmidt, & Sidebottom, 1993), si la fuerza de torsión <math display="inline">T</math> induce una rotación <math display="inline">\theta</math>  a una distancia <math display="inline">x</math> de uno de los extremos, las componentes del vector de desplazamiento pueden expresarse mediante:

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="text-align: center;vertical-align: top;"|<math display="inline">u=</math><math>\theta \cdot \Omega (y,z)</math>,          <math display="inline">v=-\theta xz</math>,          <math display="inline">w=</math><math>\theta xy</math>
|  style="text-align: right;vertical-align: top;"|(1)
|}


Aquí, <math display="inline">\Omega (y,z)</math> representa la función de alabeo. La determinación de esta función se lleva a cabo de manera que satisfaga las ecuaciones de elasticidad. Dado que estamos considerando un tensor de deformaciones infinitesimales, las ecuaciones de compatibilidad se cumplen de manera automática. Las deformaciones en un punto <math display="inline">x</math> de la viga se expresarán entonces como:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\epsilon }_{x}={\epsilon }_{y}={\epsilon }_{z}={\gamma }_{yz}=0</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\gamma }_{xy}=\theta \left( \frac{\partial \Omega }{\partial y}-z\right)</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(2)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\gamma }_{xz}=\theta \left( \frac{\partial \Omega }{\partial z}+y\right)</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|
|}


Al emplear la Ley de Hooke Generalizada, la cual establece la relación entre esfuerzos y deformaciones, derivamos de la ecuación previa las deformaciones distintas de cero de la siguiente manera:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{xy}=G\cdot {\gamma }_{xy}</math>
|}
|  rowspan='2' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(3)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{xz}=G\cdot {\gamma }_{xz}</math>
|}
|}


Si sustituimos la ecuación (3) en las ecuaciones de equilibrio interno encontradas en la mecánica de medios continuos y despreciamos las fuerzas másicas encontramos:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\frac{\partial {\tau }_{xy}}{\partial x}=0</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\frac{\partial {\tau }_{xz}}{\partial x}=0</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(4)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\frac{\partial {\tau }_{yx}}{\partial y}+\frac{\partial {\tau }_{zx}}{\partial z}=</math><math>0</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|
|}


En virtud de la conservación del momento angular, que conduce a la simetría del tensor de esfuerzos, se demuestra que los esfuerzos cortantes presentes en la ecuación (4) son independientes de la posición <math display="inline">x</math>. Como consecuencia, es necesario que exista una función de esfuerzos <math display="inline">\phi (y,z)</math>, reconocida como la función de esfuerzos de Prandtl, que satisface la relación:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{xy}=\frac{\partial \phi }{\partial z}</math>
|}
|  rowspan='2' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(5)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\tau }_{xz}=-\frac{\partial \phi }{\partial y}</math>
|}
|}


Mediante condiciones de contorno adecuadas y siguiendo los principios de la teoría descrita por (Boresi, Schmidt, & Sidebottom, 1993), es posible establecer la relación entre la fuerza de torsión <math display="inline">T</math>, los esfuerzos cortantes y la función de Prandtl <math display="inline">\phi</math> , lo que se traduce en:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
|-
| <math>T=\int_{}^{}\int_{}^{}y\cdot {\tau }_{xz}-z\cdot {\tau }_{xy}\, dy\, dz=2\, \int_{}^{}\int_{}^{}\phi dy\, dz</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(6)
|}


En términos sencillos, esta ecuación establece que el torque es equivalente a dos veces el volumen bajo la superficie formada por la función de esfuerzos y el plano de la sección transversal <math display="inline">yz</math>. Es crucial destacar que esta expresión fue derivada considerando deformaciones infinitesimales y aplicando relaciones de elasticidad lineal.

Aunque la solución específica para la función de Prandtl nos brinda la capacidad de determinar la constante de torsión en diversas secciones transversales, nuestro propósito principal radica en aplicar la teoría previamente expuesta para alcanzar soluciones numéricas en contextos donde las secciones transversales son de naturaleza arbitraria. En lugar de confinarnos a soluciones analíticas particulares, buscamos la versatilidad y precisión que ofrecen los métodos numéricos al abordar la complejidad inherente a secciones transversales no convencionales. Este enfoque nos permite extender la aplicabilidad de la teoría a una gama más amplia de situaciones prácticas en ingeniería estructural.

Con la finalidad de encontrar las constantes de torsión utilizando el método de los elementos finitos, comencemos con escribir las ecuaciones anteriores en forma matricial. También consideremos que la función de alabeo es debida a una rotación unitaria <math display="inline">\theta =</math><math>1</math>, de manera que el desplazamiento longitudinal debido a esta rotación es:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
|-
| <math>u=1\cdot \Omega (y,z)</math>
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(7)
|}


La expresión matricial de la ecuación (2) es, por tanto:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\mathit{\boldsymbol{\epsilon }}=\left[ \begin{matrix}{\gamma }_{xy}\\{\gamma }_{xz}\end{matrix}\right] =</math><math>\left[ \begin{matrix}\frac{\partial \Omega }{\partial y}\\\frac{\partial \Omega }{\partial z}\end{matrix}\right] +</math><math>\left[ \begin{matrix}-z\\y\end{matrix}\right] ={\mathit{\boldsymbol{\epsilon }}}_{w}+</math><math>{\mathit{\boldsymbol{\epsilon }}}_{0}</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(8)
|}


En esta formulación, el vector de deformaciones <math display="inline">\mathit{\boldsymbol{\epsilon }}</math> se divide en dos componentes distintas. La primera, <math display="inline">{\mathit{\boldsymbol{\epsilon }}}_{w}</math>, corresponde a las deformaciones generadas por el alabeo, mientras que la segunda, <math display="inline">{\mathit{\boldsymbol{\epsilon }}}_{0}</math>, representa las deformaciones inducidas por este fenómeno, pudiendo ser interpretada como una "carga externa". Por otro lado, la expresión (3) adquiere una forma matricial cuando se expresa de la siguiente manera:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\mathit{\boldsymbol{\sigma }}=\, \left[ \begin{matrix}{\tau }_{xy}\\{\tau }_{xz}\end{matrix}\right] =</math><math>\left[ \begin{matrix}G&0\\0&G\end{matrix}\right] \cdot \left[ \begin{matrix}{\gamma }_{xy}\\{\gamma }_{xz}\end{matrix}\right] =</math><math>\boldsymbol{C\cdot }\mathit{\boldsymbol{\epsilon }}</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(9)
|}


en la que <math display="inline">\boldsymbol{C}</math> representa la ecuación constitutiva. La ecuación anterior se puede representar por

{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="margin:auto;width: 100%;"
|-
| <math>\mathit{\boldsymbol{\sigma }}=\boldsymbol{C\cdot }\mathit{\boldsymbol{\epsilon =\, }}\boldsymbol{C\cdot }{\mathit{\boldsymbol{\epsilon }}}_{w}+</math><math>\, \boldsymbol{C\cdot }{\mathit{\boldsymbol{\epsilon }}}_{0}</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(10)
|}


En la que hemos sustituido la ecuación (8) en la ecuación (9). En esta última ecuación se aprecia de manera más evidente que la segunda expresión del lado derecho de la igualdad es equivalente a una “carga externa”.

===Discretización con elementos finitos===

La forma débil de las ecuaciones anteriores se puede expresar mediante el principio del trabajo virtual, en la que:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\delta {W}_{int}=\delta {W}_{ext}</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(11)
|}


Podemos escribir que el trabajo virtual interno es igual a:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\delta {W}_{int}=\int_{A}^{}\delta \mathit{\boldsymbol{\epsilon }}\cdot \mathit{\boldsymbol{\sigma }}\, dA</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(12)
|}


donde la variación de la deformación es igual a:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\delta \mathit{\boldsymbol{\epsilon }}=\delta \, \left( {\mathit{\boldsymbol{\epsilon }}}_{w}+\right. </math><math>\left. {\mathit{\boldsymbol{\epsilon }}}_{0}\right) =\delta {\mathit{\boldsymbol{\epsilon }}}_{w}</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(13)
|}


ya que <math display="inline">{\mathit{\boldsymbol{\epsilon }}}_{0}</math> es como una constante. Entonces, el trabajo virtual interno puede expresarse de la siguiente manera:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\delta {W}_{int}=\int_{A}^{}\delta {\mathit{\boldsymbol{\epsilon }}}_{w}\cdot \mathit{\boldsymbol{\sigma }}dA=</math><math>\int_{A}^{}\delta {\mathit{\boldsymbol{\epsilon }}}_{w}\cdot \left( \boldsymbol{C\cdot }\mathit{\boldsymbol{\epsilon }}\right) dA=</math><math>\int_{A}^{}\delta {\mathit{\boldsymbol{\epsilon }}}_{w}\cdot \left( \boldsymbol{C\cdot }{\mathit{\boldsymbol{\epsilon }}}_{w}\right) dA+</math><math>\, \int_{A}^{}\delta {\mathit{\boldsymbol{\epsilon }}}_{w}\cdot \left( \boldsymbol{C\cdot }{\mathit{\boldsymbol{\epsilon }}}_{0}\right) dA</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(14)
|}


En la ecuación (14) se puede ver que el último término que involucra a <math display="inline">{\mathit{\boldsymbol{\epsilon }}}_{0}</math> es constante, por lo que no representa un trabajo virtual interno y lo podemos mandar al lado derecho de la ecuación (11) convirtiéndose en trabajo virtual externo. De esta forma, el principio del trabajo virtual se expresa de manera siguiente:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\left[ \delta {W}_{int}=\int_{A}^{}\delta {\mathit{\boldsymbol{\epsilon }}}_{w}\cdot \boldsymbol{C}\boldsymbol{\cdot }{\mathit{\boldsymbol{\epsilon }}}_{w}dA\right] =</math><math>\left[ \delta {W}_{ext}=-\int_{A}^{}\delta {\mathit{\boldsymbol{\epsilon }}}_{w}\cdot \boldsymbol{C}\boldsymbol{\cdot }{\mathit{\boldsymbol{\epsilon }}}_{0}dA\right]</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(15)
|}


La ecuación (8) depende de las variables independientes <math display="inline">\Omega</math> , <math display="inline">y</math> y <math display="inline">z</math>, las cuales procederemos a discretizar de la siguiente manera:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>y\approx {y}^{h}=\sum _{J=1}^{{n}_{nodos}}{N}_{J}\left( y,z\right) \cdot {y}_{J}</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>z\approx {z}^{h}=\sum _{J=1}^{{n}_{nodos}}{N}_{J}\left( y,z\right) \cdot {z}_{J}</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(16)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\Omega \approx {\Omega }^{h}=\sum _{J=1}^{{n}_{nodos}}{N}_{J}\left( y,z\right) \cdot {u}_{J}</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|
|}


siendo <math display="inline">{N}_{J}(y,z)</math> las funciones de forma del nodo <math display="inline">J</math> del elemento finito bidimensional utilizado en la discretización de la sección transversal, <math display="inline">{y}_{J}</math> la coordenada en dirección del eje <math display="inline">y</math> de cada nodo <math display="inline">J</math> del elemento finito bidimensional y <math display="inline">{z}_{J}</math> su correspondiente coordenada en dirección del eje <math display="inline">z</math>. Por otro lado, ya que se esta asumiendo una rotación unitaria <math display="inline">\theta =</math><math>1</math>, entonces, de la ecuación (7), los desplazamientos son iguales a la función de alabeo, es decir <math display="inline">u=</math><math>\Omega (y,z)</math>, siendo <math display="inline">u</math> el desplazamiento axial. Como se describe en la ecuación (16), <math display="inline">{u}_{J}</math> es el desplazamiento axial del nodo <math display="inline">J</math> en el elemento finito. A continuación, expresaremos matricialmente las deformaciones de la ecuación (8) utilizando los valores ya discretizados en la ecuación (16), obteniendo las siguientes expresiones:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\mathit{\boldsymbol{\epsilon }}}_{0}\approx {\mathit{\boldsymbol{\epsilon }}}_{0}^{h}=</math><math>\left[ \begin{matrix}-{z}^{h}\\{y}^{h}\end{matrix}\right] =\sum _{J=1}^{{n}_{nodos}}\left[ \begin{matrix}-{N}_{J}\left( y,z\right) \cdot {z}_{J}\\{N}_{J}\left( y,z\right) \cdot {y}_{J}\end{matrix}\right] =</math><math>\sum _{J=1}^{{n}_{nodos}}\left[ \begin{matrix}0&-{N}_{J}\left( y,z\right) \\{N}_{J}\left( y,z\right) &0\end{matrix}\right] \cdot \left[ \begin{matrix}{y}_{J}\\{z}_{J}\end{matrix}\right]</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(17)
|}


mientras que

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\mathit{\boldsymbol{\epsilon }}}_{w}\approx {\mathit{\boldsymbol{\epsilon }}}_{w}^{h}=</math><math>\left[ \begin{matrix}\frac{\partial {\Omega }^{h}}{\partial y}\\\frac{\partial {\Omega }^{h}}{\partial z}\end{matrix}\right] =</math><math>\sum _{J=1}^{{n}_{nodos}}\left[ \begin{matrix}\frac{\partial {N}_{J}\left( y,z\right) }{\partial y}\cdot {u}_{J}\\\frac{\partial {N}_{J}\left( y,z\right) }{\partial z}\cdot {u}_{J}\end{matrix}\right] =</math><math>\sum _{J=1}^{{n}_{nodos}}\left[ \begin{matrix}\frac{\partial {N}_{J}\left( y,z\right) }{\partial y}\\\frac{\partial {N}_{J}\left( y,z\right) }{\partial z}\end{matrix}\right] \cdot \left[ {u}_{J}\right]</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(18)
|}


Sin embargo, la ecuación (18) se puede apreciar que se trata de las derivadas de las funciones de forma, cuya notación habitual en elementos finitos es mediante la matriz <math display="inline">\mathit{\boldsymbol{B}}</math>. Entonces, la ecuación (18) se puede simplificar y escribirla de la forma:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{\mathit{\boldsymbol{\epsilon }}}_{w}\approx {\mathit{\boldsymbol{\epsilon }}}_{w}^{h}=</math><math>\left[ \begin{matrix}\frac{\partial {\Omega }^{h}}{\partial y}\\\frac{\partial {\Omega }^{h}}{\partial z}\end{matrix}\right] =</math><math>\sum _{J=1}^{{n}_{nodos}}\left[ {B}_{J}\right] \cdot \left[ {u}_{J}\right]</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(19)
|}


La discretización de la ecuación (13) da como resultado la parte virtual, en la que se cambia de notación para el nodo <math display="inline">I</math> en lugar del nodo <math display="inline">J</math> y así evitar confusiones, por lo que:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\delta {\mathit{\boldsymbol{\epsilon }}}_{w}\approx {\delta \mathit{\boldsymbol{\epsilon }}}_{w}^{h}=</math><math>\sum _{I=1}^{{n}_{nodos}}\left[ {B}_{I}\right] \cdot \left[ \delta {u}_{I}\right]</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(20)
|}


Sustituyendo las expresiones anteriores en el trabajo virtual interno, encontramos:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\delta {W}_{int}={\left[ \delta {u}_{I}\right] }^{T}\cdot \int_{A}^{}\sum _{I=1}^{{n}_{nodos}}\sum _{J=1}^{{n}_{nodos}}{\left[ {B}_{I}\right] }^{T}\cdot \boldsymbol{C\cdot }\left[ {B}_{J}\right] \cdot \left[ {u}_{J}\right] \, dA</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(21)
|}


mientras que el trabajo virtual externo es igual a:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\delta {W}_{ext}={\left[ \delta {u}_{I}\right] }^{T}\cdot \int_{A}^{}\sum _{I=1}^{{n}_{nodos}}\sum _{J=1}^{{n}_{nodos}}{\left[ {B}_{I}\right] }^{T}\cdot \boldsymbol{C\cdot }\left[ \begin{matrix}0&-{N}_{J}\left( y,z\right) \\{N}_{J}\left( y,z\right) &0\end{matrix}\right] \cdot \left[ \begin{matrix}{y}_{J}\\{z}_{J}\end{matrix}\right] \, dA</math>
|}
|  rowspan='2' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(22)
|-
| 
{| style="text-align: center;vertical-align: top;margin:auto;width: 100%;"
|-
| <math>\delta {W}_{ext}={\left[ \delta {u}_{I}\right] }^{T}\cdot \int_{A}^{}\sum _{I=1}^{{n}_{nodos}}\sum _{J=1}^{{n}_{nodos}}{\left[ {B}_{I}\right] }^{T}\cdot \boldsymbol{C\cdot }{\mathit{\boldsymbol{\epsilon }}}_{0}^{h}\, dA</math>
|}
|}


Ya que los desplazamientos virtuales son arbitrarios, encontramos que las fuerzas internas son iguales a las fuerzas externas, cuya solución nos permite encontrar los desplazamientos <math display="inline">\mathit{\boldsymbol{u}}</math> que representa la incógnita del problema denominada función de alabeo. La expresión para resolver resulta ser:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\int_{A}^{}{\mathit{\boldsymbol{B}}}^{T}\cdot \mathit{\boldsymbol{C}}\cdot \mathit{\boldsymbol{B}}\, dA\cdot \mathit{\boldsymbol{u}}\mathit{\boldsymbol{=-}}\int_{A}^{}{\mathit{\boldsymbol{B}}}^{T}\cdot \mathit{\boldsymbol{C}}\cdot {\mathit{\boldsymbol{\epsilon }}}_{o}^{h}\, dA</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(23)
|}


Hay que destacar que, en planteamiento anterior, <math display="inline">\mathit{\boldsymbol{B}}=</math><math>\sum _{}^{}\left[ {\mathit{\boldsymbol{B}}}_{I}\right]</math> . La ecuación (23) resulta muy común para aquellos familiarizados con el análisis matricial de estructuras donde, la ecuación anterior resulta ser:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\mathit{\boldsymbol{K}}\cdot \mathit{\boldsymbol{u}}={\mathit{\boldsymbol{f}}}^{ext}</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(24)
|}


En la ecuación anterior resulta que la matriz <math display="inline">\mathit{\boldsymbol{K}}</math> es singular, ya que el sistema no cuenta con condiciones de contorno. La solución propuesta por (Haukaas, 2020) consiste en igualar a uno algún grado de libertad de los desplazamientos de la ecuación y después resolver para el resto de los desplazamientos. Las discretizaciones anteriores nos permiten encontrar el momento de torsión <math display="inline">T</math> definido en la ecuación (6), que escrita en forma matricial resulta:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>T=\, \int_{A}^{}\left[ \begin{matrix}-z&y\end{matrix}\right] \cdot \left[ \begin{matrix}{\tau }_{xy}\\{\tau }_{xz}\end{matrix}\right] \, dA=\, \int_{A}^{}\left[ \begin{matrix}-z&y\end{matrix}\right] \cdot \left[ \begin{matrix}G&0\\0&G\end{matrix}\right] \cdot \left[ \begin{matrix}{\gamma }_{xy}\\{\gamma }_{xz}\end{matrix}\right] \, dA</math>
|}
|  rowspan='2' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(25)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>T=\, \int_{A}^{}\left[ \begin{matrix}-{z}^{h}&{y}^{h}\end{matrix}\right] \cdot \left[ \begin{matrix}G&0\\0&G\end{matrix}\right] \cdot \left\{ \mathit{\boldsymbol{B}}\cdot \mathit{\boldsymbol{u}}+\left[ \begin{matrix}-{z}^{h}\\{y}^{h}\end{matrix}\right] \right\} \, dA</math>
|}
|}


===Aspectos computacionales===

Una vez encontrado el vector de los desplazamientos <math display="inline">\mathit{\boldsymbol{u}}</math>, se tiene que normalizar el sistema mediante la ecuación:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>\int_{A}^{}{\Omega }^{h}\left( y,z\right) \, dA=0</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(26)
|}


De esta manera, se calcula la constante de corrección por alabeo para corregir el campo de los desplazamientos. En nuestro caso por haber considerado una rotación unitaria, los desplazamientos coinciden con el resultado de nuestra función de alabeo <math display="inline">\Omega \left( y,z\right)</math> . Otra corrección al campo de los desplazamientos viene dada por el centro de cortante. De acuerdo con la torsión de St.-Venant y lo referido considerando la formulación por elementos finitos de (Haukaas, 2020) las coordenadas del centro de cortante se pueden calcular mediante:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{y}_{sc}={z}_{sc}\cdot \frac{{I}_{yz}}{{I}_{y}}-\, \frac{{\left( {z}_{K}\right) }^{T}\int_{A}^{}{N}_{I}^{T}\cdot {N}_{J}\, dA\, {u}_{K}}{{I}_{y}}</math>
|}
|  rowspan='2' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(27)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{z}_{sc}={y}_{sc}\cdot \frac{{I}_{yz}}{{I}_{z}}+\, \frac{{\left( {y}_{K}\right) }^{T}\int_{A}^{}{N}_{I}^{T}\cdot {N}_{J}\, dA\, {u}_{K}}{{I}_{z}}</math>
|}
|}


Recordando que el producto de inercia vale cero sobre los ejes principales, se pueden calcular los centros de cortante desacoplando las ecuaciones anteriores de la siguiente manera. Si la sección presenta algún eje de simetría paralelo al eje <math display="inline">y</math> o al eje <math display="inline">z</math>, entonces <math display="inline">{I}_{yz}=</math><math>0</math> y los centros de cortante valen:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{y}_{sc}=-\, \frac{{\left( {z}_{K}\right) }^{T}\int_{A}^{}{N}_{I}^{T}\cdot {N}_{J}\, dA\, {u}_{K}}{{I}_{y}}</math>
|}
|  rowspan='2' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(28)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{z}_{sc}=\, \frac{{\left( {y}_{K}\right) }^{T}\int_{A}^{}{N}_{I}^{T}\cdot {N}_{J}\, dA\, {u}_{K}}{{I}_{z}}</math>
|}
|}


En cualquier otro caso, los centros de cortante se pueden calcular si trabajamos con las inercias principales <math display="inline">{I}_{y}^{P}</math>, <math display="inline">{I}_{Z}^{P}</math> donde el producto principal de inercia <math display="inline">{I}_{yz}^{P}=</math><math>0</math>. De esta manera, los centros de cortante se calculan de acuerdo con:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{y}_{sc}=-\, \frac{{\left( {z}_{K}\right) }^{T}\int_{A}^{}{N}_{I}^{T}\cdot {N}_{J}\, dA\, {u}_{K}}{{I}_{y}^{P}}</math>
|}
|  rowspan='2' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(29)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{z}_{sc}=\, \frac{{\left( {y}_{K}\right) }^{T}\int_{A}^{}{N}_{I}^{T}\cdot {N}_{J}\, dA\, {u}_{K}}{{I}_{z}^{P}}</math>
|}
|}


Las ecuaciones anteriores revelan la necesidad de calcular los momentos de inercia para determinar los centros de cortante. Una de las primeras propuestas para este cálculo se atribuye a (Owen & Hinton, 1980), quienes sugirieron que la inercia podría obtenerse multiplicando el área de cada elemento finito por su distancia al cuadrado con respecto al eje de referencia. Aunque esta aproximación proporciona resultados cercanos al valor analítico, una metodología más precisa fue desarrollada por (Haukaas, 2020). Este último investigador consideró que las coordenadas de la sección transversal no estaban referidas al centroide de la sección, logrando así resultados prácticamente idénticos al valor analítico. No obstante, una desventaja de este enfoque surge cuando la malla no es simétrica, generando valores distintos de cero del producto de inercia para secciones transversales con un eje de simetría paralelo a uno de los ejes coordenados. Es importante señalar que la magnitud de esta discrepancia puede variar según el sistema de unidades elegido, lo que puede conducir a resultados precisos en algunas ocasiones (dentro de una tolerancia aceptable) y, en otras, introducir errores significativos no deseados.

Con el objetivo de superar esta limitación y lograr que el cálculo de las inercias no esté sujeto ni al sistema de unidades ni a mallas estructuradas, se sugiere calcularlas siguiendo el enfoque presentado por (Hally, 1987). Este autor propuso expresiones universales aplicables a cualquier polígono, proporcionando así una solución más generalizada y libre de dependencias respecto a unidades o estructuras de malla específicas, obteniéndose:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{I}_{z}=\frac{1}{12}\sum _{}^{}\left( {y}_{i}{z}_{i+1}-{y}_{i+1}{z}_{i}\right) \left( {y}_{i}^{2}+\right. </math><math>\left. {y}_{i}{y}_{i+1}+{y}_{i+1}^{2}\right)</math> 
|}
|  rowspan='3' style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(30)
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{I}_{y}=\frac{1}{12}\sum _{}^{}\left( {y}_{i}{z}_{i+1}-{y}_{i+1}{z}_{i}\right) \left( {z}_{i}^{2}+\right. </math><math>\left. {z}_{i}{z}_{i+1}+{z}_{i+1}^{2}\right)</math> 
|}
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{I}_{yz}=\frac{1}{24}\sum _{}^{}\left( {y}_{i}{z}_{i+1}-{y}_{i+1}{z}_{i}\right) \left( {y}_{i}{z}_{i+1}+\right. </math><math>\left. 2{y}_{i}{z}_{i}+2{y}_{i+1}{z}_{i+1}+{y}_{i+1}{z}_{i}\right)</math> 
|}
|}


Dado que el polígono que describe el contorno de la sección transversal en cuestión puede presentar huecos o intersecciones entre líneas, la aplicación de las ecuaciones (30) en estos escenarios queda excluida. En nuestro caso, el uso apropiado de las ecuaciones (30) implica aplicarlas individualmente a cada elemento finito de la sección transversal y luego sumar las contribuciones de cada uno para obtener el momento de inercia correcto. Es crucial destacar que las coordenadas utilizadas deben hacer referencia al centroide geométrico de la sección para garantizar la precisión buscada.

Con los cálculos anteriores se encuentra la función final de alabeo en cada nodo sumando a los desplazamientos la constante de corrección por alabeo más el centro de cortante <math display="inline">{y}_{sc}</math> multiplicado por las coordenadas en <math display="inline">z</math> del nodo referidas al centroide menos el centro de cortante <math display="inline">{z}_{sc}</math> multiplicado por las coordenadas en <math display="inline">y</math> del nodo también referidas al centroide.

===Cálculo de las constantes de torsión ===

Una vez que se ha obtenido la función final por alabeo <math display="inline">\Omega (y,z)</math>, el cálculo de la constante de torsión por alabeo <math display="inline">{C}_{w}</math> se obtiene fácilmente mediante:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{C}_{w}=\int_{A}^{}{\Omega }^{2}dA={u}_{K}\, \int_{A}^{}{N}_{I}^{T}\cdot {N}_{J}dA\, \, {u}_{K}</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(31)
|}


Por otro lado, para encontrar la constante de torsión de Saint-Venant <math display="inline">J</math>, recordemos que la relación existente entre el momento de torsión <math display="inline">T</math> y el giro de la sección <math display="inline">\theta</math>  es mediante:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>T=G\cdot J\cdot \theta</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(32)
|}


Ya que hemos supuesto que el giro es unitario, entonces podemos igualar la ecuación (25) con la ecuación (32), de donde despejamos <math display="inline">J</math> obteniedo:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>J=\, \int_{A}^{}\left[ \begin{matrix}-{z}^{h}&{y}^{h}\end{matrix}\right] \cdot \left\{ \mathit{\boldsymbol{B}}\cdot \mathit{\boldsymbol{u}}+\left[ \begin{matrix}-{z}^{h}\\{y}^{h}\end{matrix}\right] \right\} \, dA</math>
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(33)
|}


La determinación precisa de las constantes de torsión, representadas por <math display="inline">{C}_{w}</math> y <math display="inline">J</math>, para cualquier sección transversal abre las puertas a la innovación en el diseño estructural. Estas constantes proporcionan un entendimiento profundo de cómo una estructura resiste y transmite fuerzas torsionales. Este conocimiento no solo es esencial para garantizar la integridad estructural, sino que también brinda la oportunidad de explorar y crear formas geométricas únicas y estéticas que pueden embellecer significativamente las edificaciones modernas. Al comprender las características torsionales de diferentes secciones transversales, los ingenieros civiles y estructurales pueden experimentar con nuevas configuraciones, contribuyendo así a la creación de espacios habitacionales innovadores y visualmente atractivos que cumplen con los estándares contemporáneos de funcionalidad y estética. La investigación y aplicación cuidadosa de estas constantes no solo mejoran la eficiencia estructural, sino que también permiten la incorporación de elementos distintivos que enriquecen la identidad visual de las edificaciones, contribuyendo a un paisaje urbano más diverso y emocionante.

Estás constantes también permiten mejorar el comportamiento mecánico de los elementos estructurales, ya que particularmente hablando de <math display="inline">{C}_{w}</math>, se agregan nuevos miembros a la matriz de rigidez convencional de un elemento, que de acuerdo con (Haukaas, 2020) es:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>K=\left[ \begin{matrix}{k}_{aa}&{k}_{ab}\\{k}_{ba}&{k}_{bb}\end{matrix}\right]</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(34)
|}


En cada una de estas cuatro submatrices <math display="inline">({k}_{aa},\, {k}_{ab},\, {k}_{ba},\, {k}_{bb})</math> aparece un nuevo término que involucra la constante por alabeo <math display="inline">{C}_{w}</math> de la sección transversal y que da lugar a un nuevo grado de libertad conocido como bimomento en la sección. Cada una de estas cuatro matrices es igual a:

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{k}_{aa}=\left[ \begin{matrix}\frac{EA}{L}&0&0&0&0&0&0\\0&\frac{12E{I}_{z}}{{L}^{3}}&0&0&0&\frac{6E{I}_{z}}{{L}^{2}}&0\\0&0&\frac{12E{I}_{y}}{{L}^{3}}&0&-\frac{6E{I}_{y}}{{L}^{2}}&0&0\\0&0&0&\frac{12E{C}_{w}}{{L}^{3}}+\frac{6GJ}{5L}&0&0&\frac{6E{C}_{w}}{{L}^{2}}+\frac{GJ}{10}\\0&0&-\frac{6E{I}_{y}}{{L}^{2}}&0&\frac{4E{I}_{y}}{L}&0&0\\0&\frac{6E{I}_{z}}{{L}^{2}}&0&0&0&\frac{4E{I}_{z}}{L}&0\\0&0&0&\frac{6E{C}_{w}}{{L}^{2}}+\frac{GJ}{10}&0&0&\frac{4E{C}_{w}}{{L}^{3}}+\frac{2GJL}{15}\end{matrix}\right]</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(35)
|}


{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{k}_{ab}=\left[ \begin{matrix}-\frac{EA}{L}&0&0&0&0&0&0\\0&-\frac{12E{I}_{z}}{{L}^{3}}&0&0&0&-\frac{6E{I}_{z}}{{L}^{2}}&0\\0&0&-\frac{12E{I}_{y}}{{L}^{3}}&0&-\frac{6E{I}_{y}}{{L}^{2}}&0&0\\0&0&0&-\frac{12E{C}_{w}}{{L}^{3}}-\frac{6GJ}{5L}&0&0&\frac{6E{C}_{w}}{{L}^{2}}+\frac{GJ}{10}\\0&0&\frac{6E{I}_{y}}{{L}^{2}}&0&\frac{2E{I}_{y}}{L}&0&0\\0&-\frac{6E{I}_{z}}{{L}^{2}}&0&0&0&\frac{2E{I}_{z}}{L}&0\\0&0&0&-\frac{6E{C}_{w}}{{L}^{2}}-\frac{GJ}{10}&0&0&\frac{2E{C}_{w}}{L}-\frac{GJL}{30}\end{matrix}\right]</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(36)
|}


{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{k}_{ba}=\left[ \begin{matrix}-\frac{EA}{L}&0&0&0&0&0&0\\0&-\frac{12E{I}_{z}}{{L}^{3}}&0&0&0&\frac{6E{I}_{z}}{{L}^{2}}&0\\0&0&-\frac{12E{I}_{y}}{{L}^{3}}&0&\frac{6E{I}_{y}}{{L}^{2}}&0&0\\0&0&0&-\frac{12E{C}_{w}}{{L}^{3}}-\frac{6GJ}{5L}&0&0&-\frac{6E{C}_{w}}{{L}^{2}}-\frac{GJ}{10}\\0&0&-\frac{6E{I}_{y}}{{L}^{2}}&0&\frac{2E{I}_{y}}{L}&0&0\\0&\frac{6E{I}_{z}}{{L}^{2}}&0&0&0&\frac{2E{I}_{z}}{L}&0\\0&0&0&\frac{6E{C}_{w}}{{L}^{2}}+\frac{GJ}{10}&0&0&\frac{2E{C}_{w}}{L}-\frac{GJL}{30}\end{matrix}\right]</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(37)
|}


{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center;margin:auto;width: 100%;"
|-
| <math>{k}_{bb}=\left[ \begin{matrix}\frac{EA}{L}&0&0&0&0&0&0\\0&\frac{12E{I}_{z}}{{L}^{3}}&0&0&0&-\frac{6E{I}_{z}}{{L}^{2}}&0\\0&0&\frac{12E{I}_{y}}{{L}^{3}}&0&\frac{6E{I}_{y}}{{L}^{2}}&0&0\\0&0&0&\frac{12E{C}_{w}}{{L}^{3}}+\frac{6GJ}{5L}&0&0&-\frac{6E{C}_{w}}{{L}^{2}}-\frac{GJ}{10}\\0&0&\frac{6E{I}_{y}}{{L}^{2}}&0&\frac{4E{I}_{y}}{L}&0&0\\0&-\frac{6E{I}_{z}}{{L}^{2}}&0&0&0&\frac{4E{I}_{z}}{L}&0\\0&0&0&-\frac{6E{C}_{w}}{{L}^{2}}-\frac{GJ}{10}&0&0&\frac{4E{C}_{w}}{{L}^{3}}+\frac{2GJL}{15}\end{matrix}\right]</math> 
|}
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(38)
|}


El lector podrá advertir que todo lo concerniente a <math display="inline">{C}_{w}</math> representa un concepto novedoso, y se observará la presencia tanto de un renglón adicional como de una columna extra en cada submatriz, ampliando así la matriz local de la viga de 12x12 a una matriz local de 14x14. Este tipo específico de matriz rara vez se encuentra incluido en los softwares comerciales convencionales, lo cual genera un espacio propicio para el avance del conocimiento en el campo de la ingeniería estructural. La inclusión de estas dimensiones adicionales en la matriz local no solo expande las posibilidades de modelado, sino que también plantea desafíos que, al ser abordados y comprendidos, contribuirán al desarrollo y perfeccionamiento de herramientas y enfoques analíticos en el ámbito estructural.

===Ejemplos de validación===

====Sección circular====

En este ejemplo se determinarán las constantes de torsión de una sección circular que tiene un diámetro de <math display="inline">1\, m</math>. La respuesta analítica para la constante de torsión de Saint-Venant es <math display="inline">J=</math><math>\pi {D}^{4}/32</math>, que en nuestro caso valdrá <math display="inline">J=</math><math>0.09817\, {m}^{4}</math> mientras que la constante de torsión por alabeo es <math display="inline">{C}_{w}=</math><math>0</math>. Se añade la comparación del producto de inercia <math display="inline">{I}_{yz}</math> ya que es el valor que tiene más variación e importancia porque cuando no es cero, genera una rotación de los ejes. El resto de los valores obtenidos (ej.: área <math display="inline">A</math>, momentos de inercia <math display="inline">{I}_{y}</math>, <math display="inline">{I}_{z}</math>, centros de cortante <math display="inline">{y}_{sc}</math>, <math display="inline">{z}_{sc}</math>) aunque se calculan no se reportan por ser bastante parecidos. Los resultados numéricos con tres diferentes mallados (Fig. 1) y comparando con los resultados de (Haukaas, 2020) se muestran a continuación.

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  rowspan='2' style="border: 1pt solid black;text-align: center;vertical-align: top;"|
|  colspan='3'  style="border: 1pt solid black;text-align: center;"|Malla Gruesa
|  colspan='3'  style="border: 1pt solid black;text-align: center;"|Malla Media
|  colspan='3'  style="border: 1pt solid black;text-align: center;"|Malla Fina
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{I}_{yz}</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>J</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{C}_{w}</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{I}_{yz}</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>J</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{C}_{w}</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{I}_{yz}</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>J</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{C}_{w}</math>
|-
|  style="border: 1pt solid black;text-align: center;"|Haukaas
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.09e-05</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.09595</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.14e-14</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.45e-05</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.09737</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.81e-14</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.01e-06</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.09790</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.54e-14</span>
|-
|  style="border: 1pt solid black;text-align: center;"|Valdés
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.58e-11</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.09595</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.14e-14</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.83e-18</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.09737</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.81e-14</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.42e-19</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.09790</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.54e-14</span>
|}


{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;width: 33%;"|[[Image:Draft_Valdes-Vazquez_142602049-image1-c.png|96px]] 
|  style="border: 1pt solid black;text-align: center;width: 33%;"|[[Image:Draft_Valdes-Vazquez_142602049-image2-c.png|96px]] 
|  style="border: 1pt solid black;text-align: center;width: 33%;"|[[Image:Draft_Valdes-Vazquez_142602049-image3-c.png|96px]] 
|-
|  colspan='3'  style="border: 1pt solid black;text-align: center;"|Fig. 1 Mallas utilizadas
|}


Si el mismo ejemplo se repite, pero cambiando las unidades a centímetros, el diámetro será de  <math display="inline">100\, cm</math>. La constante de torsión de Saint-Venant valdrá <math display="inline">J=</math><math>9,817,477\, c{m}^{4}</math>. Se hará la comparación utilizando solamente una malla fina entre nuestra propuesta y la de Haukaas.

{| style="width: 41%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  rowspan='2' style="border: 1pt solid black;text-align: center;vertical-align: top;"|
|  colspan='3'  style="border: 1pt solid black;text-align: center;"|Malla Fina
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{I}_{yz}</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>J</math>
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<math>{C}_{w}</math>
|-
|  style="border: 1pt solid black;text-align: center;"|Haukaas
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">101.02</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">9,789,572</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.41e-06</span>
|-
|  style="border: 1pt solid black;text-align: center;"|Valdés
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.09e-11</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">9,789,572</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.41e-06</span>
|}


Como se puede ver en este último ejemplo, la propuesta de Haukaas genera productos de inercia muy diferentes de cero mientras que nuestra propuesta mejora este cálculo. Ya que se trata de una sección circular, la rotación a que da lugar el cálculo de las constantes de torsión no se ve afectada en este ejemplo debido a que, aunque los ejes se roten, las inercias no cambian. Por último, el error entre la solución analítica y la solución numérica con la malla más fina de este ejemplo es de solo el 0.28%.

====Sección rectangular====

En el segundo ejemplo se determinarán las constantes de torsión de una sección rectangular con unidades en centímetros y que tiene una base <math display="inline">b=</math><math>30\, cm.</math> y una altura <math display="inline">h=60\, cm.</math> La fórmula para calcular la torsión de Saint-Venant para este ejemplo es: <math display="inline">J=</math><math>0.229\cdot h\cdot {b}^{3}=370,980</math>. Esta fórmula no es exacta y considera que puede tener errores del 4%. Los resultados que obtenidos son:

{| style="width: 42%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  rowspan='2' style="border: 1pt solid black;text-align: center;"|
|  colspan='3'  style="border: 1pt solid black;text-align: center;"|Malla Fina
|-
|  style="border: 1pt solid black;text-align: center;"|<math>{I}_{yz}</math>
|  style="border: 1pt solid black;text-align: center;"|<math>J</math>
|  style="border: 1pt solid black;text-align: center;"|<math>{C}_{w}</math>
|-
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Haukaas</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.4860</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">375,969</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">14,863,370</span>
|-
|  style="border: 1pt solid black;text-align: center;"|Valdés
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.81e-11</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">375,969</span>
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">14,863,370</span>
|}


La malla utilizada en este ejemplo aparece en la siguiente figura.

{| style="margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;width: 100%;"|[[Image:Draft_Valdes-Vazquez_142602049-image4-c.png|102px]] 
|-
|  style="border: 1pt solid black;text-align: center;"|Fig. 2 Malla utilizada
|}


En este ejemplo, el error en el cálculo de la constante de torsión de Saint-Venant con la solución analítica (aproximada) y la solución numérica es del 1.34%. Se aprecia además que hay una importante diferencia en el cálculo del producto de inercia <math display="inline">{I}_{yz}</math>, resultando más exacto con nuestra propuesta.

====Sección arbitraria de puente====

Con la finalidad de mostrar las capacidades de nuestra propuesta, se combinado con el software de procesamiento de elementos finitos (GiD – The personal pre and post processor, 2023). Los resultados de una sección transversal que pudiera servir para un puente se muestran en las figuras 3-7. Por cuestión de programación, los ejes <math display="inline">{I}_{y},\, {I}_{z}</math> se cambian de nombre por <math display="inline">{I}_{x},\, {I}_{y}</math>.

{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;width: 100%;"|[[Image:Draft_Valdes-Vazquez_142602049-image5-c.png|600px]] 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Fig. 3 Área de la sección transversal <math display="inline">A</math>
|}


{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;width: 100%;"|[[Image:Draft_Valdes-Vazquez_142602049-image6-c.png|600px]] 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Fig. 4 Momento de inercia <math display="inline">{I}_{x}</math>
|}


{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;width: 100%;"|[[Image:Draft_Valdes-Vazquez_142602049-image7-c.png|600px]] 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Fig. 5 Momento de inercia <math display="inline">{I}_{y}</math>
|}


{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;width: 100%;"|[[Image:Draft_Valdes-Vazquez_142602049-image8-c.png|600px]] 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Fig. 6 Constante de torsión <math display="inline">J</math>
|}


{| style="width: 100%;border-collapse: collapse;" 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;width: 100%;"|[[Image:Draft_Valdes-Vazquez_142602049-image9-c.png|600px]] 
|-
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|Fig. 7 Constante de torsión <math display="inline">{C}_{w}</math>
|}


===Conclusiones===

El análisis de marcos espaciales implica la necesidad de comprender a fondo las propiedades de las secciones transversales bajo estudio, destacando especialmente el área, los momentos de inercia y las constantes de torsión. A medida que la construcción moderna incorpora cada vez más secciones compuestas o con geometrías caprichosas, el cálculo preciso de estas propiedades se convierte en un desafío fundamental para los ingenieros que se embarcan en el análisis de dichas estructuras. Para abordar esta complejidad, diversos investigadores han desarrollado técnicas numéricas avanzadas capaces de calcular con precisión las propiedades de cualquier tipo de sección transversal.

En este contexto, el presente trabajo no solo se suma a esta línea de investigación, sino que también destaca al presentar una mejora significativa en la precisión del cálculo de dichas propiedades. Esta mejora se respalda con ejemplos concretos que demuestran la eficacia del método propuesto. De este modo, se contribuye no solo a la comprensión teórica de las propiedades de las secciones transversales en marcos espaciales, sino también a la aplicación práctica de herramientas más precisas y avanzadas en el análisis estructural.

==Referencias==

Boresi, A. P., Schmidt, R. J., & Sidebottom, O. M. (1993). ''Advanced Mechanis of Materials.'' John Wiley & Sons, Inc.

Euler, L. (1744). ''Methodus Inveniendi Lineas Curvas Maximi Minimive Proprietate Gaudentes.''

GiD – The personal pre and post processor. (2023). ''Barcelona (Spain); CIMNE.'' accessed 2023 Dec 13.

Hally, D. (1987). ''Calculation of the Moments of Polygons. Technical Memorandum 87/209.'' Canadian National Defense.

Haukaas, T. (2020). ''Computational Cross-section Analysis.''

Owen, D. R., & Hinton, E. (1980). ''Finite Elements in Plasticity. Theory and Practice.'' Swansea, U.K.: Pineridge Press Limited.

Pilkey, W. D. (2002). ''Analysis and Design of Elastic Beams. Computational Methods.'' John Wiley & Sons.

Saint-Venant, A. C. (1855). ''Memoire sur la Torsion des Prismes.''

Timoshenko, S. (1934). ''Theory of Elasticity.''