MateWiki - Contribuciones del usuario [es]

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T21:24:50Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L], estudiaremos la gráfica par un intervalo de b [-10,10] y para valores positivos de b[0,10]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

[[Archivo:Bpositiva1.png|centro]]

Como podemos comprobar en la figura 6 y 7 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, hemos tomado valores
{{matlab|codigo=
clear all clc
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1)
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(r,p)

}}

A partir de este programa, hemos creado los distintos valores máximos y mínimos de las cargas críticas.

Se puede deducir que para radios pequeños la carga crítica será menor, y para radios mayores la carga crítica será mayor. Esto es porque tanto la carga crítica como la inercia dependen del radio, y al ir multiplicados, radios mayores darán cargas críticas mayores.

También se puede comprobar que en el apartado 4 de este trabajo, se llega la misma conclusión, comparando los momentos de inercia del cuadrado y del rectángulo, solo que en este caso hay que añadir que el radio es otra variable y no es constante, es variable.

'''6.2.1 Valores de la carga crítica en función del parámetro b ''

Vamos a representar de las cargas críticas dependiendo del valor b , así como la gráfica de cargas críticas para n=1 para diferentes valores de b

{{matlab|codigo=
clear all clc

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(b,A)
figure
plot(b,p)
figure
plot(r,p)
}}
[[Archivo:Bypcr.jpg|centro]]

[[Archivo:Bypcr1.jpg|centro]]

Observamos que en ambas gráficas se cumple que a mayor valor de b mayor carga crítica.

Al igual que se ha hecho en apartados anteriores hemos elaborado un gráfico que represente los radios ( que en este caso dependerán de x , R(x)) y las cargas críticas para n=1

[[Archivo:Rfrenteap1.jpg|centro]]

Estudiando la figura 11 llegamos a la misma conclusión que anteriormente, un mayor radio implica una mayor carga crítica.

'''6.2.2 Valor de R para Prc máximo y mínimo'''

Hemos hallado los valores mínimos y máximos de '''r''', R(x) y a Pcr son directamente proporcionales, por ello hemos utilizados los comando max y min para obtener estos valores.

{{matlab|codigo=

xa=0;
xb=9;
ya=0;
yb=0;

%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((8.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end

for i=1:length(r)
H(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);

end

p= H(:,1);
for i=1:length(r)
for i=1:length(r)
Rmin(i)=min(r(i,:));
end
for i=1:length(r)
Rmax(i)=max(r(i,:));
end
A(i,:)=eig(1.2.*0.25*pi.*Rmin(i)^4.*K);
P(i,:)=eig(1.2.*0.25*pi.*Rmax(i)^4.*K);
end
t= A(:,1);
m= P(:,1);

figure(1)
plot(min(Rmin),min(t),'*')
xlabel(' Rmin')
ylabel(' Prc minimo')
figure(2)
plot(max(Rmax),max(m),'*')
xlabel('Rmax')
ylabel(' Pcr maximo')

rmin=(min(Rmin))
prcmin=(min(t))
rmax=(max(Rmax))
pcrmax=(max(m))
G=[t,m]
}}

[[Archivo:Maximosss.jpg|centro]]

La carga crítica máxima se alcanza en 2.0967e+17 con un radio máximo de 3.6760e+04

[[Archivo:Minimoos.jpg|centro]]

La carga crítica mínima se alcanza en 0 con un radio minimo de 1

Hemos representados las cargas mínimas y máximas en dos columnas ( su longitud es de 91x1) pero al repetirse los valores sólo hemos hecho captura a lo más relevante

[[Archivo:Columnaspcr.png|centro]]

Para comparar la ganancia ( pcr con radio variable entre pcr con radio constante) se ha creado un programa matlab
{{matlab|codigo=
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
R=linspace(1,5,N+1);
for i=1:length(R)
B(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=B(:,1);

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
t=A(:,1);
G=100*(t./p)

}}

El programa nos saca por pantalla el valor ( %) de la ganancia
[[Archivo:Tablaganancias.png|centro]]

[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Archivo:Tablaganancias.png

2016-05-03T21:22:49Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T20:01:34Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L], estudiaremos la gráfica par un intervalo de b [-10,10] y para valores positivos de b[0,10]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

[[Archivo:Bpositiva1.png|centro]]

Como podemos comprobar en la figura 6 y 7 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, hemos tomado valores
{{matlab|codigo=
clear all clc
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1)
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(r,p)

}}

A partir de este programa, hemos creado los distintos valores máximos y mínimos de las cargas críticas.

Se puede deducir que para radios pequeños la carga crítica será menor, y para radios mayores la carga crítica será mayor. Esto es porque tanto la carga crítica como la inercia dependen del radio, y al ir multiplicados, radios mayores darán cargas críticas mayores.

También se puede comprobar que en el apartado 4 de este trabajo, se llega la misma conclusión, comparando los momentos de inercia del cuadrado y del rectángulo, solo que en este caso hay que añadir que el radio es otra variable y no es constante, es variable.

'''6.2.1 Valores de la carga crítica en función del parámetro b ''

Vamos a representar de las cargas críticas dependiendo del valor b , así como la gráfica de cargas críticas para n=1 para diferentes valores de b

{{matlab|codigo=
clear all clc

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(b,A)
figure
plot(b,p)
figure
plot(r,p)
}}
[[Archivo:Bypcr.jpg|centro]]

[[Archivo:Bypcr1.jpg|centro]]

Observamos que en ambas gráficas se cumple que a mayor valor de b mayor carga crítica.

Al igual que se ha hecho en apartados anteriores hemos elaborado un gráfico que represente los radios ( que en este caso dependerán de x , R(x)) y las cargas críticas para n=1

[[Archivo:Rfrenteap1.jpg|centro]]

Estudiando la figura 11 llegamos a la misma conclusión que anteriormente, un mayor radio implica una mayor carga crítica.

'''6.2.2 Valor de R para Prc máximo y mínimo'''

Hemos hallado los valores mínimos y máximos de '''r''', R(x) y a Pcr son directamente proporcionales, por ello hemos utilizados los comando max y min para obtener estos valores.

{{matlab|codigo=

xa=0;
xb=9;
ya=0;
yb=0;

%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((8.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end

for i=1:length(r)
H(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);

end

p= H(:,1);
for i=1:length(r)
for i=1:length(r)
Rmin(i)=min(r(i,:));
end
for i=1:length(r)
Rmax(i)=max(r(i,:));
end
A(i,:)=eig(1.2.*0.25*pi.*Rmin(i)^4.*K);
P(i,:)=eig(1.2.*0.25*pi.*Rmax(i)^4.*K);
end
t= A(:,1);
m= P(:,1);

figure(1)
plot(min(Rmin),min(t),'*')
xlabel(' Rmin')
ylabel(' Prc minimo')
figure(2)
plot(max(Rmax),max(m),'*')
xlabel('Rmax')
ylabel(' Pcr maximo')

rmin=(min(Rmin))
prcmin=(min(t))
rmax=(max(Rmax))
pcrmax=(max(m))
G=[t,m]
}}

[[Archivo:Maximosss.jpg|centro]]

La carga crítica máxima se alcanza en 2.0967e+17 con un radio máximo de 3.6760e+04

[[Archivo:Minimoos.jpg|centro]]

La carga crítica mínima se alcanza en 0 con un radio minimo de 1

Hemos representados las cargas mínimas y máximas en dos columnas ( su longitud es de 91x1) pero al repetirse los valores sólo hemos hecho captura a lo más relevante

[[Archivo:Columnaspcr.png|centro]]

{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Archivo:Columnaspcr.png

2016-05-03T20:01:11Z

NoemiP:

Archivo:Minimoos.jpg

2016-05-03T19:58:44Z

NoemiP:

Archivo:Maximosss.jpg

2016-05-03T19:56:09Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T19:54:56Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L], estudiaremos la gráfica par un intervalo de b [-10,10] y para valores positivos de b[0,10]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

[[Archivo:Bpositiva1.png|centro]]

Como podemos comprobar en la figura 6 y 7 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, hemos tomado valores
{{matlab|codigo=
clear all clc
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1)
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(r,p)

}}

A partir de este programa, hemos creado los distintos valores máximos y mínimos de las cargas críticas.

Se puede deducir que para radios pequeños la carga crítica será menor, y para radios mayores la carga crítica será mayor. Esto es porque tanto la carga crítica como la inercia dependen del radio, y al ir multiplicados, radios mayores darán cargas críticas mayores.

También se puede comprobar que en el apartado 4 de este trabajo, se llega la misma conclusión, comparando los momentos de inercia del cuadrado y del rectángulo, solo que en este caso hay que añadir que el radio es otra variable y no es constante, es variable.

'''6.2.1 Valores de la carga crítica en función del parámetro b ''

Vamos a representar de las cargas críticas dependiendo del valor b , así como la gráfica de cargas críticas para n=1 para diferentes valores de b

{{matlab|codigo=
clear all clc

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(b,A)
figure
plot(b,p)
figure
plot(r,p)
}}
[[Archivo:Bypcr.jpg|centro]]

[[Archivo:Bypcr1.jpg|centro]]

Observamos que en ambas gráficas se cumple que a mayor valor de b mayor carga crítica.

Al igual que se ha hecho en apartados anteriores hemos elaborado un gráfico que represente los radios ( que en este caso dependerán de x , R(x)) y las cargas críticas para n=1

[[Archivo:Rfrenteap1.jpg|centro]]

Estudiando la figura 11 llegamos a la misma conclusión que anteriormente, un mayor radio implica una mayor carga crítica.

'''6.2.2 Valor de R para Prc máximo y mínimo'''

Hemos hallado los valores mínimos y máximos de '''r''', R(x) y a Pcr son directamente proporcionales, por ello hemos utilizados los comando max y min para obtener estos valores.

{{matlab|codigo=

xa=0;
xb=9;
ya=0;
yb=0;

%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((8.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end

for i=1:length(r)
H(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);

end

p= H(:,1);
for i=1:length(r)
for i=1:length(r)
Rmin(i)=min(r(i,:));
end
for i=1:length(r)
Rmax(i)=max(r(i,:));
end
A(i,:)=eig(1.2.*0.25*pi.*Rmin(i)^4.*K);
P(i,:)=eig(1.2.*0.25*pi.*Rmax(i)^4.*K);
end
t= A(:,1);
m= P(:,1);

figure(1)
plot(min(Rmin),min(t),'*')
xlabel(' Rmin')
ylabel(' Prc minimo')
figure(2)
plot(max(Rmax),max(m),'*')
xlabel('Rmax')
ylabel(' Pcr maximo')

rmin=(min(Rmin))
prcmin=(min(t))
rmax=(max(Rmax))
pcrmax=(max(m))
G=[t,m]
}}
{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T19:30:57Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L], estudiaremos la gráfica par un intervalo de b [-10,10] y para valores positivos de b[0,10]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

[[Archivo:Bpositiva1.png|centro]]

Como podemos comprobar en la figura 6 y 7 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, hemos tomado valores
{{matlab|codigo=
clear all clc
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1)
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(r,p)

}}

A partir de este programa, hemos creado los distintos valores máximos y mínimos de las cargas críticas.

Se puede deducir que para radios pequeños la carga crítica será menor, y para radios mayores la carga crítica será mayor. Esto es porque tanto la carga crítica como la inercia dependen del radio, y al ir multiplicados, radios mayores darán cargas críticas mayores.

También se puede comprobar que en el apartado 4 de este trabajo, se llega la misma conclusión, comparando los momentos de inercia del cuadrado y del rectángulo, solo que en este caso hay que añadir que el radio es otra variable y no es constante, es variable.

'''6.2.1 Valores de la carga crítica en función del parámetro b ''

Vamos a representar de las cargas críticas dependiendo del valor b , así como la gráfica de cargas críticas para n=1 para diferentes valores de b

{{matlab|codigo=
clear all clc

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(b,A)
figure
plot(b,p)
figure
plot(r,p)
}}
[[Archivo:Bypcr.jpg|centro]]

[[Archivo:Bypcr1.jpg|centro]]

Observamos que en ambas gráficas se cumple que a mayor valor de b mayor carga crítica.

Al igual que se ha hecho en apartados anteriores hemos elaborado un gráfico que represente los radios ( que en este caso dependerán de x , R(x)) y las cargas críticas para n=1

[[Archivo:Rfrenteap1.jpg|centro]]

Estudiando la figura 11 llegamos a la misma conclusión que anteriormente, un mayor radio implica una mayor carga crítica.

'''6.2.2 Valor de R para Prc máximo y mínimo'''

Hemos hallado los valores mínimos y máximos de '''r''', R(x) y a Pcr son directamente proporcionales, por ello hemos utilizados los comando max y min para obtener estos valores.

{{matlab|codigo=
xa=0;
xb=9;
ya=0;
yb=0;

%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((8.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end

for i=1:length(r)
H(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);

end

p= H(:,1);
for i=1:length(r)
for i=1:length(r)
Rmin(i)=min(r(i,:));
end
for i=1:length(r)
Rmax(i)=max(r(i,:))
end
A(i,:)=eig(1.2.*0.25*pi.*Rmin(i)^4.*K);
P(i,:)=eig(1.2.*0.25*pi.*Rmax(i)^4.*K);
end
t= A(:,1);
m= P(:,1);

for i=1:length(t)
pmax(i)=max(t(i,:));
end
for i=1:length(m)
pmin(i)=min(m(i,:));
end

figure(1)
plot(b,r)
figure(2)
plot(r,pmax)
figure(3)
plot(r,pmin)
figure(4)
plot(p,r)
}}

{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T19:30:09Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L], estudiaremos la gráfica par un intervalo de b [-10,10] y para valores positivos de b[0,10]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

[[Archivo:Bpositiva1.png|centro]]

Como podemos comprobar en la figura 6 y 7 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, hemos tomado valores
{{matlab|codigo=
clear all clc
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1)
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(r,p)

}}

A partir de este programa, hemos creado los distintos valores máximos y mínimos de las cargas críticas.

Se puede deducir que para radios pequeños la carga crítica será menor, y para radios mayores la carga crítica será mayor. Esto es porque tanto la carga crítica como la inercia dependen del radio, y al ir multiplicados, radios mayores darán cargas críticas mayores.

También se puede comprobar que en el apartado 4 de este trabajo, se llega la misma conclusión, comparando los momentos de inercia del cuadrado y del rectángulo, solo que en este caso hay que añadir que el radio es otra variable y no es constante, es variable.

'''6.2.1 Valores de la carga crítica en función del parámetro b ''

Vamos a representar de las cargas críticas dependiendo del valor b , así como la gráfica de cargas críticas para n=1 para diferentes valores de b

{{matlab|codigo=
clear all clc

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(b,A)
figure
plot(b,p)
figure
plot(r,p)

[[Archivo:Bypcr.jpg|centro]]

[[Archivo:Bypcr1.jpg|centro]]

Observamos que en ambas gráficas se cumple que a mayor valor de b mayor carga crítica.

Al igual que se ha hecho en apartados anteriores hemos elaborado un gráfico que represente los radios ( que en este caso dependerán de x , R(x)) y las cargas críticas para n=1

[[Archivo:Rfrenteap1.jpg|centro]]

Estudiando la figura 11 llegamos a la misma conclusión que anteriormente, un mayor radio implica una mayor carga crítica.

'''6.2.2 Valor de R para Prc máximo y mínimo'''

Hemos hallado los valores mínimos y máximos de '''r''', R(x) y a Pcr son directamente proporcionales, por ello hemos utilizados los comando max y min para obtener estos valores.

{{matlab|codigo=
xa=0;
xb=9;
ya=0;
yb=0;

%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((8.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end

for i=1:length(r)
H(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);

end

p= H(:,1);
for i=1:length(r)
for i=1:length(r)
Rmin(i)=min(r(i,:));
end
for i=1:length(r)
Rmax(i)=max(r(i,:))
end
A(i,:)=eig(1.2.*0.25*pi.*Rmin(i)^4.*K);
P(i,:)=eig(1.2.*0.25*pi.*Rmax(i)^4.*K);
end
t= A(:,1);
m= P(:,1);

for i=1:length(t)
pmax(i)=max(t(i,:));
end
for i=1:length(m)
pmin(i)=min(m(i,:));
end

figure(1)
plot(b,r)
figure(2)
plot(r,pmax)
figure(3)
plot(r,pmin)
figure(4)
plot(p,r)
}}

{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T19:24:07Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L], estudiaremos la gráfica par un intervalo de b [-10,10] y para valores positivos de b[0,10]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

[[Archivo:Bpositiva1.png|centro]]

Como podemos comprobar en la figura 6 y 7 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, hemos tomado valores
{{matlab|codigo=
clear all clc
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1)
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(r,p)

}}

A partir de este programa, hemos creado los distintos valores máximos y mínimos de las cargas críticas.

Se puede deducir que para radios pequeños la carga crítica será menor, y para radios mayores la carga crítica será mayor. Esto es porque tanto la carga crítica como la inercia dependen del radio, y al ir multiplicados, radios mayores darán cargas críticas mayores.

También se puede comprobar que en el apartado 4 de este trabajo, se llega la misma conclusión, comparando los momentos de inercia del cuadrado y del rectángulo, solo que en este caso hay que añadir que el radio es otra variable y no es constante, es variable.

'''6.2.1 Valores de la carga crítica en función del parámetro b ''

Vamos a representar de las cargas críticas dependiendo del valor b , así como la gráfica de cargas críticas para n=1 para diferentes valores de b

{{matlab|codigo=
clear all clc

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(b,A)
figure
plot(b,p)
figure
plot(r,p)

[[Archivo:Bypcr.jpg|centro]]

[[Archivo:Bypcr1.jpg|centro]]

Observamos que en ambas gráficas se cumple que a mayor valor de b mayor carga crítica.

Al igual que se ha hecho en apartados anteriores hemos elaborado un gráfico que represente los radios ( que en este caso dependerán de x , R(x)) y las cargas críticas para n=1

[[Archivo:Rfrenteap1.jpg|centro]]

Estudiando la figura 11 llegamos a la misma conclusión que anteriormente, un mayor radio implica una mayor carga crítica.

{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T19:22:11Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L], estudiaremos la gráfica par un intervalo de b [-10,10] y para valores positivos de b[0,10]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

[[Archivo:Bpositiva1.png|centro]]

Como podemos comprobar en la figura 6 y 7 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, hemos tomado valores
{{matlab|codigo=
clear all clc
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1)
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(r,p)

}}

A partir de este programa, hemos creado los distintos valores máximos y mínimos de las cargas críticas.

Se puede deducir que para radios pequeños la carga crítica será menor, y para radios mayores la carga crítica será mayor. Esto es porque tanto la carga crítica como la inercia dependen del radio, y al ir multiplicados, radios mayores darán cargas críticas mayores.

También se puede comprobar que en el apartado 4 de este trabajo, se llega la misma conclusión, comparando los momentos de inercia del cuadrado y del rectángulo, solo que en este caso hay que añadir que el radio es otra variable y no es constante, es variable.

'''6.2.1 Valores de la carga crítica en función del parámetro b ''

Vamos a representar de las cargas críticas dependiendo del valor b , así como la gráfica de cargas críticas para n=1 para diferentes valores de b

{{matlab|codigo=
clear all clc

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(b,A)
figure
plot(b,p)
figure
plot(r,p)

[[Archivo:Bypcr.jpg|centro]]

[[Archivo:Bypcr1.jpg|centro]]

Observamos que en ambas gráficas se cumple que a mayor valor de b mayor carga crítica.

Al igual que se ha hecho en apartados anteriores hemos elaborado un gráfico que represente los radios ( que en este caso dependerán de x , R(x)) y las cargas críticas para n=1

[[Archivo:Rfrenteap1.jpg|centro]]

Estudiando la figura 11 llegamos a la misma conclusión que anteriormente, un mayor radio implica una mayor carga crítica.

'''6.2.2 Valor mínimo de Pcr en función del radio'''
[[Archivo:Figura8pcr.png|centro]]

El valor mínimo que alcanza la carga crítica es '''0.1148'''

'''6.2.3 Valor máximo de Pcr en función del radio'''

[[Archivo:Figura9radiosmax.png|centro]]

El valor máximo de la carga crítica es '''11483'''

{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Archivo:Rfrenteap1.jpg

2016-05-03T19:20:45Z

NoemiP:

Archivo:Rfrenteap.jpg

2016-05-03T19:15:15Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T19:13:28Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L], estudiaremos la gráfica par un intervalo de b [-10,10] y para valores positivos de b[0,10]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

[[Archivo:Bpositiva1.png|centro]]

Como podemos comprobar en la figura 6 y 7 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, hemos tomado valores
{{matlab|codigo=
clear all clc
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1)
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(r,p)

}}

A partir de este programa, hemos creado los distintos valores máximos y mínimos de las cargas críticas.

Se puede deducir que para radios pequeños la carga crítica será menor, y para radios mayores la carga crítica será mayor. Esto es porque tanto la carga crítica como la inercia dependen del radio, y al ir multiplicados, radios mayores darán cargas críticas mayores.

También se puede comprobar que en el apartado 4 de este trabajo, se llega la misma conclusión, comparando los momentos de inercia del cuadrado y del rectángulo, solo que en este caso hay que añadir que el radio es otra variable y no es constante, es variable.

'''6.2.1 Valores de la carga crítica en función del parámetro b ''

Vamos a representar de las cargas críticas dependiendo del valor b , así como la gráfica de cargas críticas para n=1 para diferentes valores de b

{{matlab|codigo=
clear all clc

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(b,A)
figure
plot(b,p)
figure
plot(r,p)

[[Archivo:Bypcr.jpg|centro]]

[[Archivo:Bypcr1.jpg|centro]]

Observamos que en ambas gráficas se cumple que a mayor valor de b mayor carga crítica.

'''6.2.2 Valor mínimo de Pcr en función del radio'''
[[Archivo:Figura8pcr.png|centro]]

El valor mínimo que alcanza la carga crítica es '''0.1148'''

'''6.2.3 Valor máximo de Pcr en función del radio'''

[[Archivo:Figura9radiosmax.png|centro]]

El valor máximo de la carga crítica es '''11483'''

{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Archivo:Bypcr1.jpg

2016-05-03T19:11:05Z

NoemiP:

Archivo:Bypcr.jpg

2016-05-03T19:10:02Z

NoemiP:

Archivo:Bpositiva1.png

2016-05-03T18:52:02Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T11:38:00Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

Como podemos comprobar en la figura 6 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, hemos tomado valores de b positivos:
{{matlab|codigo=
clear all clc
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
;
%Parte 2. Paso y tamaño de paso

h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1)
L=9;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end

for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*r(i)^4.*K);
%La matriz A es EI/M*K
end
end
p=A(:,1);
figure
plot(b,r)
figure
plot(r,p)

}}

A partir de este programa, hemos creado los distintos valores máximos y mínimos de las cargas críticas.

Se puede deducir que para radios pequeños la carga crítica será menor, y para radios mayores la carga crítica será mayor. Esto es porque tanto la carga crítica como la inercia dependen del radio, y al ir multiplicados, radios mayores darán cargas críticas mayores.

También se puede comprobar que en el apartado 4 de este trabajo, se llega la misma conclusión, comparando los momentos de inercia del cuadrado y del rectángulo, solo que en este caso hay que añadir que el radio es otra variable y no es constante, es variable.

'''6.2.1 Valores de la carga crítica en función del radio'''

[[Archivo:Figura7pyr2.png|centro]]

Como podemos comprobar en la figura 7, para cada valor de R(x), variará nuestra carga máxima. Podemos extrapolarlo a dos casos importantes: Cuando el radio es mínimo (entre 0 y 10) y para el radio máximo (0.5 a 3.5x〖10〗^4)

'''6.2.2 Valor mínimo de Pcr en función del radio'''
[[Archivo:Figura8pcr.png|centro]]

El valor mínimo que alcanza la carga crítica es '''0.1148'''

'''6.2.3 Valor máximo de Pcr en función del radio'''

[[Archivo:Figura9radiosmax.png|centro]]

El valor máximo de la carga crítica es '''11483'''

{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Archivo:Figura9radiosmax.png

2016-05-03T11:04:07Z

NoemiP:

Archivo:Figura8pcr.png

2016-05-03T10:59:10Z

NoemiP:

Archivo:Figura7pyr2.png

2016-05-03T10:56:45Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T10:55:23Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

Como podemos comprobar en la figura 6 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, así como las gráficas de los valores máximos y mínimos:
{{matlab|codigo=
clear all clc

%Parte 1. Condiciones de Contorno
xa=0;xb=9;ya=0;yb=0;
%Parte 2. Paso y tamaño de paso
h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((8.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
for i=1:length(r)
Rsm(i)=min(r(i,:));
end
for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*Rsm(i)^4.*K);
end
p= A(:,1);

for i=1:length(p)
pmax(i)=max(p(i,:));
end
for i=1:length(p)
pmin(i)=min(p(i,:));
end
for i=1:length(r)
rmax(i)=max(r(i,:));
end
for i=1:length(r)
rmin(i)=min(r(i,:));
end
}}

A partir de este programa, hemos creado los distintos valores máximos y mínimos de las cargas críticas.

Se puede deducir que para radios pequeños la carga crítica será menor, y para radios mayores la carga crítica será mayor. Esto es porque tanto la carga crítica como la inercia dependen del radio, y al ir multiplicados, radios mayores darán cargas críticas mayores.

También se puede comprobar que en el apartado 4 de este trabajo, se llega la misma conclusión, comparando los momentos de inercia del cuadrado y del rectángulo, solo que en este caso hay que añadir que el radio es otra variable y no es constante, es variable.

'''6.2.1 Valores de la carga crítica en función del radio'''

[[Archivo:Figura7pyr.png|centro]]

{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Archivo:Figura7pyr.png

2016-05-03T10:54:34Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T10:51:08Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

Como podemos comprobar en la figura 6 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, así como las gráficas de los valores máximos y mínimos:
{{matlab|codigo=
clear all clc

%Parte 1. Condiciones de Contorno
xa=0;xb=9;ya=0;yb=0;
%Parte 2. Paso y tamaño de paso
h=0.1;
N=(xb-xa)/h;
x=xa:h:xb;
b=linspace(1,10,N+1);
L=
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((8.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5);
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
for i=1:length(r)
Rsm(i)=min(r(i,:));
end
for i=1:length(r)
A(i,:)=eig(1.2.*0.25*pi.*Rsm(i)^4.*K);
end
p= A(:,1);

for i=1:length(p)
pmax(i)=max(p(i,:));
end
for i=1:length(p)
pmin(i)=min(p(i,:));
end
for i=1:length(r)
rmax(i)=max(r(i,:));
end
for i=1:length(r)
rmin(i)=min(r(i,:));
end
}}
{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T10:49:58Z

NoemiP:

{{ TrabajoED | Carga crítica de una columna. Grupo 5-A | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED15/16|Curso 2015-16]] | Noemi Palomino

Maria Lucia Cueva }}

'''1. Introducción'''

En este tema se estudia el fenómeno de inestabilidad de piezas simples solicitadas a compresión centrada. Se trata por lo tanto de la primera aproximación a los fenómenos de pandeo de pilares, en la que se restringe el enfoque el caso de cortante puro.
Así mismo, en este tema se tratará de forma general, la carga crítica a la que se somete una columna se sección circular constante de una pieza simple e ideal. Aunque la información puede ser extrapolada a los pilares con sección variable, su caso se analizará más adelante.
Para la exposición del pandeo por composición pura de piezas simples se comenzará con los fundamentos teóricos del problema.

'''2. La pieza ideal'''

Se denomina pieza ideal a una abstracción matemática de los pilares reales, basada en las siguientes hipótesis:

• El material tiene una rigidez constante caracterizada por un módulo de deformación '''E'''.

• El material es infinitamente elástico, sin que se alcance nunca un límite de elasticidad, o tensión de agotamiento.

• La pieza no tiene imperfecciones de ninguna clase, ajustándose de manera perfecta a su directriz teórica.

• No existen tensiones residuales dentro de la sección

Con estas simplificaciones se pueden plantear de manera sencilla las ecuaciones que gobiernan el problema de inestabilidad, y su resolución matemática

'''3. Carga crítica'''

Vamos a analizar la posible inestabilidad de la columna biapoyada de la figura siguiente. Es evidente que la solución trivial de la pieza consiste en permanecer recta, sin deformación. Sin embargo, la posible aparición de una flecha y(x), hace que el axil (P) provoque un momento flector que a su vez puede incrementar la flecha. Para que la deformación sea posible, debe encontrarse una solución de equilibrio sobre la pieza deformada, entre los esfuerzos solicitantes y la respuesta de la sección.
Esta solución de equilibrio es un análisis de segundo orden. El planteamiento del equilibrio en cualquier sección de la pieza deformada es inmediato:

Momento solicitante:
[[Archivo:Formulamomento.png]]

'''EI''': Módulo de elasticidad lineal (o módulo de Young) que depende de las propiedades elásticas del material, y que supondremos constante

'''I(x)''': Momento de inercia de la sección transversal respecto del centro

'''3.1 Estudio de la pieza de equilibrio'''

Para que la pieza esté en equilibrio, se debe resolver la ecuación diferencial de segundo orden. Para resolverla, suponemos que los desplazamientos son pequeños y '''aproximamos y(x)''', resolviendo la ecuación de la curva elástica que proviene del equilibrio de momentos:

[[Archivo:3.1foto.png|centro]]

Se trata de una ecuación diferencial homogénea que admite solución general:
[[Archivo:Solgeneral.png|sinmarco|centro]]

[[Archivo:Simplificado.png|600x600px|centro]]

'''3.2 Problema de contorno de la pieza de equilibrio. Concepto de carga crítica'''

Para obtener la solución se deben particularizar las condiciones de contorno:
[[Archivo:Condicionescontorno.png|sinmarco|centro]]

Lo que da lugar a:

[[Archivo:Ayb.png|sinmarco|centro]]
Este sistema tiene una solución trivial que es '''A=0''', con lo que la pieza está en equilibrio en la situación sin deformar, basado en la inexistencia de momentos flectores.
Sin embargo, también admite una solución cuando la función trigonométrica se anula:
[[Archivo:EcuacionSeno.png|sinmarco|centro]]

En este caso, el parámetro A queda indeterminado, pudiendo adoptar cualquier valor, este es un problema de contorno homogéneo, por ello tendremos que determinar los valores de P que hacen que el problema de contorno tenga infinitas soluciones. Esto supone que la pieza admite infinitas posiciones de equilibrio en situación deformada, cuando se alcanzan unos determinados valores del axil aplicado (carga crítica). Hemos conseguido resolver el problema de autovalores asociado, dados por la ecuación siguiente:

[[Archivo:Pcrr.png|centro]]

A partir de esta ecuación podemos observar que la deformada de la pieza vendrá dada por los distintos valores de cargas críticas, y que estos dependen a su vez de '''n (n=1,2,3…)'''

[[Archivo:Fotosnodos1.png|centro]]

En la siguiente figura se ha hecho una representación de los diferentes axiles críticos. Se muestran los casos de los tres axiles críticos con los valores de n=1,2,3. La deformada que representa en la figura se le llama '''modo de pandeo'''.

Con ello se explica que para efectos prácticos, la carga crítica relevante es la del primer modo, ya que una vez alcanzada la flexión, la pieza seguirá deformándose ante cualquier perturbación transversal. Se denomina '''carga crítica de Euler''':

'''4. Estabilidad de P'''

Se dice que una columna es estable cuando la único solución posible es la trivial '''y(x)=0.'''

'''4.1 Ejemplo de estabilidad de una columna con sección circular'''

Sea

[[Archivo:Datosproblemas.png|sinmarco|centro]]

Para estudiar la estabilidad de la columna, necesitamos saber cuáles son los valores de P, para ello sustituimos en (1), los datos enunciados:

[[Archivo:Pcrdatos.png|500x500px|centro]]

Podemos deducir que habrá estabilidad cuando la carga '''no tenga un valor superior a 0,0470 KN'''

'''4.1.1 Gráfica de estabilidad en función del radio'''

Para comprobar que la carga crítica que se ha calculado es correcta, vamos a realizar un programa que asocia a cada radio el valor más pequeño de la carga crítica para el cual la columna deja de ser estable.

{{matlab|codigo=

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)=1.2*(pi^3*R(1,k).^4)/(4*(9^2));
end
plot(R,P,'b-')
xlabel('Radio (m)');
ylabel('Carga critica (kN)')
}}

[[Archivo:Radios2.jpg|800x800px|sinmarco|centro]]

La carga crítica, se obtiene particularizando con n=1, como hemos comentado anteriormente, ya que es la carga que provoca el pandeo del pilar. En el apartado anterior, solo hay un radio, sin embargo para este ejemplo se ha utilizado un intervalo de radios (de 1 a 5)

En la gráfica se observa que un radio mayor permite un mayor valor para la carga crítica.

'''4.2 Ejemplo de estabilidad de una columna con sección cuadrada'''

Otro ejemplo particular, es el estudio de una columna con una sección cuadrada. Este ejemplo se resuelve de la misma forma que el apartado 4.1 solo que el momento de inercia cambia:

[[Archivo:Momentoinerciarectanfula.png|centro]]

Vamos a suponer que la base (B), tiene la misma longitud que el radio de la circunferencia, por lo que nuestro nuevo valor para la carga crítica será:

[[Archivo:Pcrrectangulo.png|centro]]

Podemos comprobar que en una sección cuadrada la carga crítica disminuye, esto se debe al momento de inercia

'''5.Estabilidad por diferencias finitas'''

Sea nuestra ecuación diferencial

[[Archivo:Ecuacione2.png|centro]]

Escribimos una aproximación por diferencias finitas del apartado 4.1:

El método de Diferencias Finitas es un método de carácter general que permite la resolución aproximada de ecuaciones diferencias en derivadas parciales definidas en recintos finitos.

En este caso vamos a calcular una aproximación de la carga crítica. Vamos a calcular el problema matricialmente de forma que podamos representarlo numéricamente
[[Archivo:Numericamente.png|centro]]

Aproximando:

[[Archivo:Aprox2.png|centro]]

De esta forma, la ecuación diferencial nos quedará de la forma:

[[Archivo:Ecuacionaprox.png|centro]]

Y pasándolo a forma matricial:

[[Archivo:Taylo2.png|centro]]

[[Archivo:Matrices2.png|centro]]

Por lo tanto, la matriz de coeficientes nos queda de la forma:

[[Archivo:Matrizcoef.png|centro]]

Se trata por lo tanto, un problema en el que tendremos que aproximar nuestro autovalor (que es la carga crítica) y encontrar una solución de éste:

[[Archivo:Problemaautovalor.png|centro]]

Numéricamente, se emplea el comando “eig”, que nos permite calcular los valores propios de la matriz A:

{{matlab|codigo=
%Calcular numéricamente por diferencias finitas
%los valores de P(cr). Para ellos se realizará
%un problema de autovalores

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo del problema de contorno
plot(R,p,'r')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura3).png|centro]]
{{matlab|codigo=
clear all, clc

%Caso 1. Calculado a partir del autovalor Pr

L=input('introduce una longitud: ');
r=0:0.1:5;
%función del pcr
Pcr=((pi^3).*(r.^4).*1.2)./(4.*(L^2))

%Caso calculado por aproximación de diferencias finitas
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*0.25*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

%Dibujo de la comparación de los 2 programas
plot(R,p,'r',r,Pcr,'*')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Figura4).png|centro]]

En la figura 4 vemos que las gráficas coinciden. Esto es porque la aproximación numérica difieren a partir del tercer decimal, por lo que podemos considerar ese error despreciable.

Es decir, calculando la carga crítica con n=1 por autovalores o aproximándolo por diferencias finitas, la carga crítica difiere en '''〖10〗^(-3) m''', considerando así los valores calculados correctos.

'''5.1 Comparación de la Carga crítica en sección circular y rectangular'''

Uno de los estudios más interesantes es la comparación de la carga crítica en función de su sección

{{matlab|codigo=

%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);

%Parte 4. Bucle
for i=1:length(R)
A(i,:)=eig(1.2.*(0.25)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
p=A(:,1);

for i=1:length(R)
A(i,:)=eig(1.2.*(1/12)*pi*R(i)^4.*K);
%La matriz A es EI/M*K
end
m=A(:,1);
%Dibujo del problema de contorno
plot(R,p,'r',R,m,'b')
xlabel('Radio (m)')
ylabel('Carga crítica (kN)')
}}

[[Archivo:Comparacionciryrec1.png|centro]]

En la figura 5 Observamos que la sección rectangular admite menos carga que la circular

'''6. Pandeo en sección variable'''

En ocasiones, a la hora de dimensionar diferentes columnas nos encontramos con el caso de columnas de sección variable, y por tanto de inercia variable, en los que no se puede aplicar de la misma forma el cálculo del pandeo.

Un ejemplo de secciones variables pueden ser soportes diseñados con un solo perfil o con varios perfiles enlazados entre sí de manera continua, y con la sección con un área ligeramente variable.

El cálculo de estos soportes se realiza obteniendo un valor del radio, y determinado como se indica en el siguiente apartado

'''6.1 Pandeo con R(x) variable'''

A continuación, se va a analizar el pandeo de un cambio de sección en la columna planteando el caso de radio variable más sencillo: un polinomio

[[Archivo:Seccionvariable.png|centro]]

Vamos a comprobar que R(x) se puede escribir en función de un único parámetro '''b'''

[[Archivo:Integral1.png|centro]]

'''6.1.1 Gráfica de R(x) para diferentes valores de b'''

Como R(x) es variable, vamos a tener distintos valores del radio en función de b como hemos indicado en el apartado 1.1. Ahora, se va a dibujar numéricamente los radios de las curvas de [0,L]

{{matlab|codigo=
%Dibujar R(x) para diferentes valores de b
%Parte 1. Condiciones de Contorno
xa=0;
xb=9;
ya=0;
yb=0;
%Parte 2. Paso y tamaño de paso
N=100;
h=(xb-xa)/N;
x=xa:h:xb;
%Parte 3.
%%Matriz K
K=(1/h^(2))*(2*diag(ones(N-1,1))-diag(ones(N-2,1),1)-diag(ones(N-2,1),-1));
%%Vector de longitudes (Radio de 1 a 5)
R=linspace(1,5,N+1);
b=linspace(-10,10,N+1);
L=9;
%Parte 4. Bucle

for k=1:length(b)
a=(((-4.*b(k)*((L/2)).^3)/3)-sqrt((4.*b(k)*(L/2).^3)/3).^2+((4.*(L/2)^5)/5).*((25.83/pi)+b(k).^2*L))/((4*(L/2)^5)/5)
for j=1:length(x)
r(k,j)=a.*((x(j)-L/2)^2)+b(k);
end
end
%Dibujo del problema de contorno
plot(b,r)
}}

[[Archivo:Grafica6ryb.png|centro]]

Como podemos comprobar en la figura 6 tendremos un radio distinto para distintos valores de b. Sin embargo, para cálculos posteriores solo se tendrá en cuenta los valores positivos. También se puede comprobar que la gráfica de radios siempre será positiva para cualquier valor de b.

'''6.2 Valores de la carga crítica para los diferentes valores de b'''

Otro de los estudios interesantes es comprobar para que valores de R(x) la carga Crítica será máxima o mínima. En el cálculo de piezas de hormigón resulta de importancia, ya que podremos saber cuál es la carga máxima que necesitaremos aplicar a una columna para evitar que produzca pandeo.

A continuación, se detallará el programa que se ha utilizado para su cálculo, así como las gráficas de los valores máximos y mínimos:

{{ Beta }}
[[Categoría:Ecuaciones Diferentiales]]
[[Categoría:ED15/16]]
[[Categoría:Trabajos 2015-16]]

Archivo:Grafica6ryb.png

2016-05-03T10:48:09Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T10:46:28Z

NoemiP:

Archivo:Gradicaryb.png

2016-05-03T10:37:07Z

NoemiP:

Archivo:Radioyb.jpg

2016-05-03T10:35:52Z

NoemiP:

Archivo:Integral1.png

2016-05-03T10:30:53Z

NoemiP:

Archivo:Seccionvariable.png

2016-05-03T10:29:01Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T09:53:24Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T09:52:42Z

NoemiP:

Archivo:Comparacionciryrec1.png

2016-05-03T09:49:39Z

NoemiP:

Archivo:Comparacionciryrec.png

2016-05-03T09:29:47Z

NoemiP:

Archivo:Problemaautovalor.png

2016-05-03T09:11:21Z

NoemiP:

Archivo:Matrizcoef.png

2016-05-03T09:07:31Z

NoemiP:

Archivo:Matrices2.png

2016-05-03T09:05:37Z

NoemiP:

Archivo:Taylo2.png

2016-05-03T09:04:07Z

NoemiP:

Archivo:Numericamente.png

2016-05-03T08:58:23Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T08:55:26Z

NoemiP:

Archivo:Ecuacione2.png

2016-05-03T08:54:33Z

NoemiP:

Archivo:Ecuacione.png

2016-05-03T08:51:25Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T08:48:21Z

NoemiP:

Archivo:Fotosnodos1.png

2016-05-03T08:47:06Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-03T08:42:59Z

NoemiP:

Archivo:3.1foto.png

2016-05-03T08:40:22Z

NoemiP:

Carga crítica de una columna, Trabajo 9, (Grupo 5)

2016-05-02T22:13:17Z

NoemiP:

Archivo:Figura4).png

2016-05-02T22:09:12Z

NoemiP:

Archivo:Figura3).png

2016-05-02T22:06:53Z

NoemiP: