|
|
| (No se muestran 72 ediciones intermedias del mismo usuario) |
| Línea 1: |
Línea 1: |
| − | {{ TrabajoED | Carga crítica de una columna. Grupo 7B | [[:Categoría:Ecuaciones Diferenciales|Ecuaciones Diferenciales]]|[[:Categoría:ED13/14|Curso 2013-14]] | María Ramírez
| |
| − | Ignacio Posada
| |
| | | | |
| − | Antonio López-Mateos
| |
| − |
| |
| − | Pablo Bueno }}
| |
| − |
| |
| − | = Introducción =
| |
| − | En este proyecto numérico hemos estudiado la estabilidad y deflexión de una columna vertical, sometida a una carga constante P.
| |
| − |
| |
| − | La base de la teoría de las columnas es la formula de Euler, que fue publicada en 1757 por Leonardo Euler, un matemático suizo. Esta formula es válida solamente para columnas largas. También tenemos que suponer que esta columna es recta, homogénea, y de sección transversal constante en toda su longitud.
| |
| − | Se aplica la ley de Hooke y los esfuerzos son inferiores al limite de proporcionalidad del material, para evitar fracturas en la estructura de la columna.
| |
| − |
| |
| − | La carga critica de pandeo es la carga axial máxima que una columna puede soportar cuando esta a punto de pandearse.
| |
| − |
| |
| − | Un esquema representativo sería:
| |
| − |
| |
| − | [[Archivo:Columna1.JPG|marco|centro|Esquema de la deflexión de la columna]]
| |
| − |
| |
| − | = Planteamiento matemático (LISTO)=
| |
| − |
| |
| − | Examinando los datos iniciales del problema, nos encontramos ante una columna vertical larga y esbelta de sección transversal uniforme, y de longitud total igual a L. Su eje de simetría, por tanto, ocupa el intervalo espacial '''<math>x \epsilon [0, L]</math>''', donde '''<math>x=0</math>''' corresponde al
| |
| − | extremo superior de la columna, donde se aplica la carga, y '''<math>x=L</math>''' al extremo inferior de la columna. Sea '''<math>y(x)</math>''' la '''función elástica''' o función que expresa la curvatura del eje de la columna al aplicarle una fuerza vertical de compresión o carga '''<math>P</math>''', en su extremo superior (ver figura Apartado 1). Al comparar los momentos flectores en cualquier punto de la columna se obtiene:
| |
| − |
| |
| − | \begin{matrix}
| |
| − | y''(x)=\frac{M(x)}{E I(x)}
| |
| − | \end{matrix}
| |
| − |
| |
| − | Siendo:
| |
| − |
| |
| − | · '''<math>E</math>''': módulo de elasticidad de Young.
| |
| − |
| |
| − | · '''<math>I(x)</math>''': momento de inercia de una sección trasversal respecto a una recta vertical por el centro.
| |
| − |
| |
| − | · '''<math>M(x)</math>''': momento flector.
| |
| − |
| |
| − | En el caso de la columna, el momento depende de la deflexión, de manera que: ''<math>M(x)=-Py(x)</math>''.
| |
| − |
| |
| − | Por tanto, sustituyendo la fórmula del momento flector en la ecuación inicial:
| |
| − | \[\begin{array}{crl}
| |
| − | \\
| |
| − | y''(x)=\frac{M(x)}{E I(x)} \Longrightarrow y''(x)=\frac{-Py(x)}{E I(x)} \Longrightarrow y''(x)+\frac{Py(x)}{E I(x)}=0 \\
| |
| − | \\
| |
| − | \end{array}\]
| |
| − | Para poder resolver este problema de contorno faltarían las condiciones de contorno que serían:
| |
| − | <math>
| |
| − | \left\{\begin{matrix} y(0)=0 \\ y(L)=0\\ \end{matrix}\right.
| |
| − | </math>
| |
| − | ya que en los extremo de la viga: ''<math>x=0</math>'' y ''<math>x=L</math>'' tiene deflexión cero.
| |
| − |
| |
| − |
| |
| − | Por tanto, el problema de contorno planteado para resolver el ejercicio sería:
| |
| − |
| |
| − | <math>
| |
| − | \left\{\begin{matrix} y''(x)+\frac{Py(x)}{E I(x)}=0 \\ y(0)=0 \\ y(L)=0\\ \end{matrix}\right.
| |
| − | </math>
| |
| − |
| |
| − | = Estabilidad de la columna =
| |
| − | Vamos a estudiar ahora la estabilidad de la columna en función de la carga P aplicada. El problema que tenemos que resolver para determinar los valores de P que proporcionan estabilidad a la misma, o lo que es lo mismo, los puntos donde la columna no flecta,es el problema de contorno ya planteado en el Apartado 2:
| |
| − |
| |
| − | <math>
| |
| − | \left\{\begin{matrix} y''(x)+\frac{Py(x)}{E I(x)}=0 \\ y(0)=0 \\ y(L)=0\\ \end{matrix}\right.
| |
| − | </math>
| |
| − |
| |
| − |
| |
| − | Suponemos, según los datos proporcionados, que el módulo de elasticidad de Young es '''<math>E=1</math>''', que la sección de la columna circular es de radio constante '''<math>R=1\ m</math>''', que su densidad es '''<math>\rho=1\ Kg/m^3</math>''' y que la longitud total de la columna es '''<math>L=10\ m.</math>'''.
| |
| − |
| |
| − | Con el problema de contorno ya planteado, nos disponemos ahora a resolverlo por los procedimientos conocidos.
| |
| − | (Nota: para facilitar los cálculos, sustituimos las variables que dependen de '''x''' por variables independientes, como una herramienta de cálculo, esto es '''<math>y(x)=y</math>''', '''<math>I(x)=I</math>''').
| |
| − |
| |
| − | Por tanto, el problema quedaría
| |
| − |
| |
| − | \[\begin{array}{crl}
| |
| − | \\
| |
| − | E Iy''+Py=0 \\
| |
| − | \\
| |
| − | \end{array}\]
| |
| − |
| |
| − |
| |
| − | Resolución:
| |
| − |
| |
| − | \[\begin{array}{crl}
| |
| − | \\
| |
| − | E Iy''+Py=0 \Longrightarrow E Ir^2+P=0 \Longrightarrow r^2=\frac{-P}{E I}=0 \Longrightarrow r=\; \pm\;\sqrt{\frac{P}{E I}} i \\
| |
| − | \\
| |
| − | \end{array}\]
| |
| − | La solución general de la ecuación diferencial homogénea será
| |
| − | \[\begin{array}{crl}
| |
| − | \\
| |
| − | y\;= y(x)\; =\; A\; cos\; (\sqrt{\frac{P}{E I}} x) + \; B\; sin\; (\sqrt{\frac{P}{E I}} x) = 0\\
| |
| − | \\
| |
| − | \end{array}\]
| |
| − |
| |
| − | Teniendo en cuenta las condiciones de contorno iniciales '''<math>y(0)=0</math>''' y '''<math>y(L)=0</math>''' y
| |
| − | que, por definición, se dice que una columna es estable si la única solución posible de la función elástica '''<math>y(x)</math>''' es la trivial, '''<math>y(x)=0</math>''', sustituimos dichos valores en la solución general obtenida y llegamos a la siguiente ecuación
| |
| − | \[\begin{array}{crl}
| |
| − | \\
| |
| − | y(x)\; =\; A\; cos\; (0) + \; B\; sin\; (0) = 0 \\
| |
| − | \\
| |
| − | \end{array}\]
| |
| − |
| |
| − | = Dibujar la función que a cada radio R 2 [1; 5] le asocia el valor más pequeño de Pcr para el cual la columna deja de ser estable. El valor Pcr se conoce como carga crítica (obtenido por primera vez por L. Euler)(3) =
| |
| − |
| |
| − | == Proceso de obtencion de la ecuacion ==
| |
| − |
| |
| − | == Resolucion numerica ==
| |
| − |
| |
| − | = Punto 5 de despejar en funcion del radio =
| |
| − |
| |
| − | == Proceso de obtencion de la ecuacion ==
| |
| − |
| |
| − | == Comparacion de valores de Pcr segun el radio ==
| |
| − |
| |
| − | = Resolucion por Fourier =
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − |
| |
| − | [[Categoría:Ecuaciones Diferenciales]]
| |
| − | [[Categoría:ED13/14]]
| |
| − | [[Categoría:Trabajos 2013-14]]
| |
