Estudio de la carga crítica de una columna

Trabajo realizado por estudiantes
Título	Estudio de la carga crítica de una columna. Grupo 14-B
Asignatura	Ecuaciones Diferenciales
Curso	Curso 2013-14
Autores	Manuel Campayo Fernández, Alejandro Arturo Serrano Leo, David Suárez Cuesta, Arturo Rodríguez Gárate, Jorge Martín Sebastián, Joaquín Sánchez Molina
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

Se nos plantea analizar cual es la carga crítica de una columna, para lo cual realizamos el siguiente estudio. Cogemos una columna homogénea y de sección transversal circular. La colocamos de forma vertical y aplicamos sobre ella una carga P. El eje de simetría de la columna, que se extiende desde x=0 (extremo superior de la viga) hasta x=L (extremo inferior), sufrirá con la carga P desplazamientos horizontales (eje Y) respecto a la posición de equilibrio. Para medir estos desplazamientos utilizaremos una función y(x) denominada deflexión, la cual (si los desplazamientos son pequeños) es posible aproximar mediante la resolución de la ecuación de la curva elástica que proviene del equilibrio de momentos y que mostramos a continuación:: [math] y''(x) = \frac{M(x)}{EI(x)} \left\{\begin{matrix} E=Módulo De Elasticidad Lineal. Lo Supondremos Constante. \\ I(x)=Momento De Inercia De La Sección Trasversal Restecto Al Centro \\ M(x)=Momento Flector \end{matrix}\right. [/math]

En este caso el momento depende de la deflexión, por lo que [math] M(x)=-PI(x) [/math]

1 Problema de contorno para y(x)

El problema físico planteado no se resuelve mediante una ecuación diferencial en derivadas parciales, sino con un problema de contorno, que en nuestro caso, aplicándole todas las condiciones antes explicadas queda de la siguiente manera:: [math] (PA) = \left\{\begin{matrix} E*I(x)*y''(x)=M(x) \\ y(0)=y(L)=0 \end{matrix}\right\}. [/math]

La condición y(0)=y(L)=0 viene determinada por el hecho de que la viga no sufre ningún tipo de desplazamiento horizontal en sus extremos.

Sustituyendo en lo anterior la condición de [math] M(x)=-PI(x) [/math] , el problema de contorno que buscamos queda:: [math] (PA) = \left\{\begin{matrix} E*I(x)*y''(x)+P*y(x)=0 \\ y(0)=y(L)=0 \end{matrix}\right\}. [/math]

2 Valores de P que hacen estable a la columna

Una vez tenemos el problema de contorno planteado, queremos obtener para qué valores de P la columna es estable. Para ello, suponemos que E=1, la sección de la columna circular es de radio constante R=1 m, su densidad es ρ=1[math]Kg/m^3[/math] y L=10 m.

Antes de nada, se debe aclarar que una columna es estable si la única solución posible de la elástica es la trivial [math] y(x)=0 [/math].

El problema (PA) es homogéneo, por tanto siempre verificará la solución trivial. Ahora bien, este problema es estable sólo si verifica esa y ninguna otra, por tanto, cuando haya infinitas soluciones, el problema será inestable. Esto se traduce en un pandeo de la columna y en el posterior derrumbamiento de la estructura. Por ende, sólo tendremos que determinar los valores de P que hacen que el problema de contorno (PA) tenga infinitas soluciones.

Sustituyendo los valores de L, E, R y densidad que nos dan en el enunciado, el problema de autovalores a determinar queda:: [math] (PA) = \left\{\begin{matrix} \frac{1}{4} π*y''(x)+P*y(x)=0 \\ y(0)=y(10)=0 \end{matrix}\right\}. [/math]

Resolviendo el problema de autovalores obtenemos los valores de la carga P que hacen que (PA) tenga infinitas soluciones, provocando la inestabilidad de la barra y su posterior pandeo. La solución del problema de autovalores es la siguiente:

Los autovalores son : [math] P_k = \frac{k^2*π^3}{400} [/math] y las autofunciones asociadas son : [math] Φ_k = A*sin(\frac{k*π*x}{10}) [/math]

con k=1,2,...,N y A número real.

Por tanto, los valores de P que hacen la columna estable son todos los P menores o iguales a 0 y aquellos que no verifican:: [math] P_k = \frac{k^2*π^3}{400} [/math] con k=1,2,...,N.

Así la columna sería estable.

3 Función de las cargas críticas mínimas para las cuales la columna deja de ser estable

En el apartado anterior, estudiábamos el problema para R=1. Sin embargo, en este caso existen varios valores de R a estudiar, concretamente los pertenecientes al intervalo [1,5]. Si resolvemos el problema de autovalores (PA) en función del radio (aportado por la expresión del momento de inercia del círculo), obtenemos la expresión de los autovalores en función del radio de la sección del pilar.

Los autovalores son : [math] P_k = \frac{k^2*π^3*R^4}{400} [/math] y las autofunciones asociadas son : [math] Φ_k = A*sin(\frac{k*π*x}{10}) [/math]

con k=1,2,...,N y A número real.

Introduciendo esa expresión en Matlab (código 1) y evaluándola en el intervalo [1,5], obtenemos los valores de P_crit. Para ello, en vez de tomar únicamente cinco valores de R, hemos tomado un salto más pequeño, en concreto h=0.1. Por tanto, el código utilizado es:

%Representamos los valores de la carga critica en funcion del radio de la sección.
clear; clf; clc;
r1=1;
r2=5;
h=0.01;
N=(r2-r1)/h;
R=r1:h:r2;
P=zeros(1,N+1)
for k=1:(N+1)
P(1,k)=(pi^3*R(1,k).^4)/400
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')

Introduciendo la función en Matlab, la gráfica que se obtiene es la siguiente:

Al haber tomado un valor tan pequeño del módulo de elasticidad, E=1, y resolviendo, los primeros valores de la gráfica son muy pequeños, pero todos distintos de cero.

Nota: la carga crítica (P_crit) se obtiene particularizando la expresión de los autovalores para k=1 y físicamente se corresponde con la carga que provoca el pandeo del pilar. El resto de los autovalores (k=2,3, …) tienen valores más elevados y su significado físico es más complejo, el cual tampoco es interesante determinar puesto que la sobrecarga que hay que introducir para que se produzca es mucho más elevada que en el caso del primer autovalor. Estudiamos el problema para el caso k=1.