Deformación Axial en Barras Empotradas: Análisis bajo Cargas Puntuales y Distribuidas.

En el estudio riguroso de la RESISTENCIA DE MATERIALES II, comprender cómo se deforman los elementos estructurales bajo solicitaciones complejas es el primer paso para un diseño seguro. Desde los laboratorios de la Universidad de Cuenca hasta las aplicaciones de ingeniería en campo, el análisis de columnas, pilares y tirantes requiere ir más allá de las fórmulas estáticas básicas.

En Engineering Vault Pro, abordamos hoy un clásico del análisis estructural: el cálculo del desplazamiento axial total en el extremo libre de una barra empotrada que está sometida a dos tipos de cargas simultáneas. Este modelo es fundamental para evaluar elementos que soportan tanto fuerzas aplicadas externamente como su propio peso.

Antes de adentrarnos en la deducción matemática, es vital establecer el vocabulario técnico y las variables que gobernarán nuestro sistema:

• $u(L)$: Deformación o desplazamiento total en el extremo libre de la barra (ubicado en $x = L$).
• $F$: Carga axial puntual aplicada directamente en el extremo libre.
• $w$: Carga axial uniformemente distribuida a lo largo de toda la longitud (fuerza por unidad de longitud, representando frecuentemente el peso propio).
• $L$: Longitud total del elemento estructural.
• $E$: Módulo de Elasticidad o Módulo de Young, que define la rigidez intrínseca del material.
• $A_0$: Área de la sección transversal de la barra, la cual asumimos constante.

Para llegar a la ecuación final de desplazamiento, los ingenieros estructurales suelen utilizar el Principio de Superposición (calculando los efectos por separado y sumándolos) o la Integración Directa.

Optaremos por la integración directa. Este método es matemáticamente más puro y nos permite demostrar de dónde surge cada término del comportamiento elástico del material, analizando la fuerza interna sección por sección.

Imagina nuestra barra empotrada en el origen $x = 0$ y completamente libre en su extremo $x = L$. Para entender qué ocurre en su interior, realizamos un "corte" imaginario a una distancia genérica $x$ desde el empotramiento.

Nuestro objetivo es determinar la fuerza interna $N(x)$ que actúa exactamente en esa sección transversal. Para que este fragmento de la estructura se mantenga en equilibrio estático, la fuerza interna $N(x)$ debe resistir y ser igual a la suma de todas las fuerzas externas aplicadas a la derecha de nuestro corte (es decir, en el tramo desde $x$ hasta $L$).

En este tramo actúan dos solicitaciones:

1. La carga puntual $F$ en el extremo libre.
2. La porción de la carga distribuida $w$ contenida en ese segmento restante, cuya magnitud total equivale a multiplicar la carga por la longitud del segmento: $w(L - x)$.

Ensamblando estos componentes, la ecuación analítica de la fuerza interna en cualquier coordenada $x$ es:

$$N(x) = F + w(L - x)$$

Una vez conocida la distribución de la fuerza interna, recurrimos a la Ley de Hooke. En el régimen elástico lineal, sabemos que un diferencial de deformación $du$ que ocurre en un diferencial de longitud $dx$ está dictado por la fuerza, la geometría y el material:

$$du = \frac{N(x)}{EA_0} dx$$

Para hallar el desplazamiento total acumulado en el extremo libre $u(L)$, debemos integrar los diferenciales de deformación a lo largo de toda la barra, planteando los límites desde el empotramiento ($x = 0$) hasta la punta ($x = L$):

$$u(L) = \int_{0}^{L} \frac{N(x)}{EA_0} \,dx$$

El siguiente paso es sustituir nuestra función de fuerza interna $N(x)$ dentro de la integral definida:

$$u(L) = \int_{0}^{L} \frac{F + w(L - x)}{EA_0} \,dx$$

Dado que estamos analizando una sección prismática de un material homogéneo, tanto el módulo $E$ como el área $A_0$ son constantes. Podemos extraerlos de la integral por linealidad y separar el polinomio del numerador para facilitar el cálculo:

$$u(L) = \frac{1}{EA_0} \left( \int_{0}^{L} F \,dx + \int_{0}^{L} (wL - wx) \,dx \right)$$

Procedemos a integrar matemáticamente respecto a la variable $x$:

$$u(L) = \frac{1}{EA_0} \left[ Fx \Big|_0^L + \left( wLx - \frac{wx^2}{2} \right) \Bigg|_0^L \right]$$

Evaluamos la expresión aplicando el Teorema Fundamental del Cálculo en los límites superior ($L$) e inferior ($0$):

$$u(L) = \frac{1}{EA_0} \left[ FL + \left( wL(L) - \frac{wL^2}{2} \right) \right]$$

Simplificando los términos algebraicos, obtenemos:

$$u(L) = \frac{1}{EA_0} \left[ FL + wL^2 - \frac{wL^2}{2} \right]$$

Al restar los términos semejantes de la carga distribuida ($1 - 1/2 = 1/2$), el corchete se reduce a su forma más elegante:

$$u(L) = \frac{1}{EA_0} \left[ FL + \frac{wL^2}{2} \right]$$

Finalmente, al distribuir el denominador común, llegamos a la ecuación magistral que rige este fenómeno estructural:

$$u(L) = \frac{FL}{EA_0} + \frac{wL^2}{2EA_0}$$

La belleza de esta demostración radica en la claridad con la que la matemática refleja la realidad física. La ecuación final de desplazamiento no es más que la superposición perfecta de dos comportamientos aislados.

El primer término, $\frac{FL}{EA_0}$, representa exclusivamente la deformación axial aportada por la carga puntual aplicada en el extremo. Por su parte, el segundo término, $\frac{wL^2}{2EA_0}$, cuantifica la deformación generada por la carga distribuida a lo largo de toda la geometría.

Comprender esta separación de variables permite a los ingenieros programar rutinas de análisis más eficientes y auditar el comportamiento de estructuras sujetas a grandes cargas muertas. La automatización de estos principios es lo que separa a un cálculo tradicional de un diseño optimizado de alto nivel.