[Escribir el título del documento]
[Año]
[Escriba aquí una descripción breve del documento. Una descripción breve es un resumen corto del contenido del documento. Escriba aquí una descripción breve del documento. Una descripción breve es un resumen corto del contenido del documento.]
UNIVERSIDAD AUTÓNOMA DE CHIHUAHUA
Facultad De Ingeniería
Materia: Ecuaciones Diferenciales
Fecha de entrega 25/11/2014
Deflexión De Una Viga
Se resolverá la deflexión de una viga con valores en la frontera por medio de la transformada de Laplace.
Catedrática: Ing. Karla Bojorquez Gutiérrez
Alumnos Matricula
Irving Martin Salcido Medrano 282103
Alejandro Romo González 281972
Marcos Gallegos Hernández 281834
Introducción
La transformada de Laplace es un operador lineal sumamente útil a la hora de resolver de manera más eficaz las ecuaciones diferenciales de orden superior las cuales tienen una gran aplicación en diversas ramas de las matemáticas y de la física.
En esta ocasión abordaremos una de sus aplicaciones en ingeniería civil, la cual se centra en las vigas. Estas últimas son un elemento fundamental en la construcción, no solo soportan presión y peso, sino también flexión y tensión, además han ayudado a construir muchas estructuras incluso en el mundo antiguo.
El problema que resolveremos se trata sobre la deflexión de una viga con valores en la frontera, por lo tanto para iniciar necesitamos saber más sobre estos elementos estructurales lineales o unidimensionales de la construcción.
Después nos adentraremos en el modelo matemático del fenómeno que analizaremos (deflexión de una viga), veremos los métodos con los cuales vamos a solventar el problema que se nos plantea.
Por último nos centraremos en la resolución del problema en particular, realizando los cálculos necesarios y explicándolos de manera detallada, una vez hecho esto mostraremos la gráfica que se genera con el ejercicio caso particular que se nos ha planteado, para así pasar a la conclusión del problema.
Vigas
La viga es una estructura horizontal que puede sostener carga entre 2 apoyos sin crear empuje lateral en estos. El uso más importante de estas, es quizás el que se aplica a la estructura de puentes. Son un elemento estructural, es decir, forma parte del diseño de una estructura rigiéndose por los principios de la resistencia de materiales y de la ingeniería.
Son elementos lineales, los cuales también son llamados prismas mecánicos o unidimensionales. Estos son alargados y son sometidos a un estado de tensión plana.
Por lo anterior dicho se sabe que las vigas estas sometidas a una tensión por lo tanto con los distintos materiales se comportara de una forma diferente por ejemplo el acero hace a las vigas mas rígidas, las de aluminio son más flexibles y las de madera tienen mayor elasticidad, no obstante cualquier viga se romperá cuando se aplica una cantidad de presión excesiva.
El uso de vigas está sumamente extendido por lo cual se utilizan en la construcción desde rascacielos a estadios, aunque también se pueden utilizar en la construcción residencial por lo tanto existe una gran variedad de tipos y clasificación de las vigas, que van desde el tipo de material, forma en que se colocan hasta el tipo de soporte necesario para cada estructura.
Deflexión de una viga
Se entiende por deflexión a la deformación que sufre un elemento por el efecto de las flexiones internas, tornándose en una curvatura o desviación de un curso o línea horizontal.
La problemática de la deflexión de vigas y ejes en una determinada estructura, es un tema muy importante para los ingenieros y en especial para los ingenieros civiles que al construir una edificación pueden presentar inconvenientes debido a la deflexión de las vigas de la obra que se esté realizando.
Para prevenir daños como la deflexión de las vigas es necesario identificar adecuadamente cada uno de los factores que pueden llegar a tener un gran impacto en una edificación en el futuro, lo que ocasionaría un colapso de la edificación.
La deflexión, matemáticamente es una función y(x) que está gobernada por una ecuación diferencial lineal de cuarto orden.
Considérese una viga de longitud L de un material homogéneo y que también es uniforme en su sección transversal tal como se muestra en la figura a. El eje de simetría se indica por la línea punteada.
Si se aplica una carga a la viga en un plano vertical que contiene al eje de simetría, la viga experimenta una distorsión y la curva que conecta los centroides de las secciones transversales se llama curva de deflexión o curva elástica, esto se muestra en la figura b.
La curvatura de deflexión se aproxima a la forma de una viga. El eje x coincide con el eje de simetría de la viga y que la deflexión y(x), medida desde este eje es positiva si es hacia abajo.
En la teoría de elasticidad se muestra que el momento de flexión M(x) en un punto x lo largo de la viga se relaciona con la carga por unidad de longitud w(x) mediante la ecuación:
d2Mdx2=w(x)
Además, el momento de flexión M(x) es proporcional a la curvatura de κ de la curva elástica.
Mx=EIκ
Donde E e I son constantes; E es el modulo de Young de elasticidad del material de la viga e I es el momento de inercia de una sección transversal de la viga. El producto se llama rigidez flexional de la viga.
La curvatura está dada por κ=y''/[1+(y')2)32. Cuando la deflexión y(x) es pequeña, la pendiente se acerca a cero, y por tanto y''/[1+(y')2)32 se acerca a uno. Si se permite que κ y'', la ecuación Mx=EIκ se convierte en M=EIy''. La segunda derivada de esta última expresión es:
d2Mdx2=EId2dx2y''=EId4ydx4
Por lo tanto la deflexión y(x) satisface la ecuación diferencial de cuarto orden.
EId4ydx4=w(x)
Las condiciones de frontera (contorno o región donde está definida la ecuación diferencial) asociadas con esta ultima ecuación depende de los apoyos extremos de la viga.
Empotrada en ambos extremos
Las condiciones en la frontera en x=0 y x=L son y=0, y'=0
Libres
Las condiciones en la frontera del lado en voladizo son x=0 y x=L son y''=0, y'''=0
Apoyados
Las condiciones en la frontera en x=0 y x=L son y=0, y''=0
Planteamiento del problema
El problema que se nos plantea en el proyecto es el siguiente:
Una viga esta empotrada en su extremo izquierdo y simplemente apoyado en el derecho.
Encuentre la deflexión y(x) cuando la carga está dada por:
wx=w0, 0
Graficar la solución cuando w0=48EI y la viga mide un metro de largo. Considere que la viga es de granito y el momento de inercia es de I=1kgm2
Existen varias maneras de resolver este tipo de situaciones, a nosotros se nos indico que debe de ser resuelta por el método de la transformada de Laplace y una función llamada función escalón unitario o función de Heaviside, la cual se explicara mejor.
Función escalón unitario
En ingeniería es común encontrar funciones que están ya sea "activadas" o "desactivadas". Es conveniente entonces definir una función especial la cual se llama función escalón unitario o función de Heaviside en honor al matemático inglés Oliver Heaviside. Es una función matemática que tiene como característica, el tener un valor de 0 para todos los números negativos (desactivado) y de 1 para todos los positivos (activado), esta función se aplica en la fuerza externa que actúa en un sistema mecánico, o un voltaje aplicado a un circuito.
La función escalón unitario U(t-a) se define como.
Ut-a=0, 0 t
Para nuestro problema solo tomaremos el eje t positivo tal como muestra la figura.
Cando una función f definida para t 0 se multiplica por u=(t-a), la función escalo unitario "desactiva" una parte de la grafica de esa función.
La función escalón unitario también se puede usar para escribir funciones definidas por tramos en forma compacta, tenemos 2 tipos distintos los cuales son:
ft=gt, 0 t
Estas se puede escribir como: ft=gt-gtut-a+htu(t-a)
Y análogamente tenemos.
ft=0, 0 t
ft=gt[ut-a-u(t-b)]
Transformada de Laplace
La transformada de Laplace puede ser usada para resolver ecuaciones diferenciales lineales y ecuaciones integrales, comúnmente a problemas con coeficientes constantes.
Cuando se resuelven estas ecuaciones usando la técnica de la transformada, se cambia de una ecuación diferencial a una ecuación algebraica.
Sea f una función definida para t 0, la transformada de Laplace de f(t) se define como
Lf(t)=0+ e-stf(t)dt
Transformada de una derivada
Si necesitamos encontrar la transformada de Laplace de una derivada de orden n de una funcion f(t) se aplica lo siguiente.
L{dndtnf(t)}=SnFs-Sn-1f0-Sn-2f'0-Sn-3f''0-…-fn-1(0)
Resolución del problema
La función es a tramos por lo tanto convertimos w(x)=w0, 0
a su forma compacta a partir de la función escalón, por lo cual para dejar un cero en la primera parte tal como en la función de Heaviside, de la función a tramos extraemos w0 y realizamos la operación necesaria para hacer que la segunda parte de la función siga siendo cero.
wx=wo0, 0
La segunda parte de la función a tramos es la que se multiplica por la función de escalón unitario dando como resultado lo siguiente.
w0-w0u(x-L2)
Ya que nuestra viga esta empotrada en x=0 y apoyado simplemente en x=L tenemos las siguientes condiciones:
y0=0, y'0=0, y''L=0 Y y'''L=0
Igualamos w(x) a la ecuación diferencial de cuarto orden que satisface la deflexión de una viga.
EIy''''=w0[1-u(x-L2)]
Aplicamos la transformada de una derivada a la ecuación anterior.
EILy''''=w0L[1-u(x-L2)]
EIS4YS-S3y0-S2y'0-Sy''0-y'''(0)
Aplicando linealidad en la segunda parte de la ecuación.
w0L1-Lux-L2
Resolviendo.
L1=1S
Con el segundo teorema de traslación.
Lux-L2 ft=1, a=L2 e-LS21S
La expresión total seria.
EIS4YS-S3y0-S2y'0-Sy''0-y'''(0)=w01S-e-LS2S
Dadas nuestras condiciones en la frontera y0 y y'0 serán cero por lo tanto la ecuación se simplifica a lo siguiente.
EIS4YS-Sy''0-y'''=w01S-e-LS2S
De aquí hacemos y''=C1 y y'''=C2 y despejamos Y(S), a continuación se muestran el resultado de las operaciones antes mencionadas.
EIS4YS-SC1-C2=w01S-e-LS2S
S4YS-SC1-C2=w0EI1S-e-LS2S
S4YS=SC1-C2+w0EI1S-e-LS2S
YS=SC1S4-C2S4+w0EI1S-e-LS2SS4
YS=C1S3+C2S4+w0EI1S5-e-LS2S5
Aplicamos la transformada inversa para encontrar y(x).
L-1Y(S)=L-1C1S3+C2S4+w0EI1S5-e-LS2S5
L-1Y(S)=y(x)
Aplicando linealidad.
L-1C11S3 n=2, 2!=2 C12x2
L-1C21S4 n=3, 3!=6 C26x3
L-11S5 n=4, 4!=24 124x4
L-1-eLS21S5 n=4, 4!=24, a=L2 -124x-L24
La expresión general es.
yx=C12x2+c26x3+w0EI124x4-124(x-L2)4u(x-L2)
Tomamos nuestras segundas condiciones yL=0 y y''L=0 por lo tanto igualamos a cero nuestra expresión general y sustituimos nuestras variables por la longitud L.
Aplicamos la primera condición.
yL=C12L2+C26L3+w0EI124L4-124(L-L2)4u(L-L2)=0
Ya que L>L2 no "desactiva" la función por lo cual no se multiplica por la función escalón unitario, entonces:
yL=C12L2+C26L3+w0EI124L4-124(L-L2)4=0
yL=C12L2+C26L3+w0EI124L4-1384L4=0
yL=C12L2+C26L3+5w0L4128EI=0
Aplicamos la segunda condición por lo tanto sacamos la segunda derivada.
y(x)'=C1x+C22x2+w0EI16x3-16x-L23
y(x)''=C1+C2x+w0EI12x2-12x-L22
y(L)''=C1+C2L+w0EI12L2-12L-L22=0
y(L)''=C1+C2L+w0EI12L2-18L2=0
y(L)''=C1+C2L+3w0L28EI=0
Entonces tenemos nuestro sistema de ecuaciones para averiguar el valor de las constantes.
C12L2+C26L3+5w0L4128EI=0
C1+C2L+3w0L28EI=0
Resolviendo el sistema de ecuaciones, despejamos "C1" en la segunda ecuación.
C1=-C2L-3w0L28EI
Sustituimos el valor obtenido en la primera ecuación.
(-C2L-3w0L28EI)2L2+C26L3+5w0L4128EI=0
(-C2L2-3w0L216EI)L2+C26L3+5w0L4128EI=0
-C2L32-3w0L416EI+C26L3+5w0L4128EI=0
Dejamos los términos con "C2" de un lado y del otro, los términos independientes.
-C22L3+C26L3=3w0L416EI-5w0L4128EI
Sacamos "C2" y "w0L4EI" como factor común.
C2(-L32+L36)=3w0L416EI-5w0L4128EI
C2(-L33)=w0L4EI(316-5128)
C2(-L33)=w0L4EI(19128)
Ahora despejamos "C2".
C2=w0L4EI(19128)(-L33)
Realizamos la división y obtenemos el valor de "C2".
C2=-(57128)w0LEI
Lo sustituimos en la segunda ecuación.
C1=-(-(57128)w0LEI)L-3w0L28EI
Resolvemos.
C1=57128w0L2EI-(38)w0L2EI
C1=w0L2EI(57128-38)
C1=(9128)w0L2EI
Sustituimos los valores de C1 y C2 en la expresión general.
yx=(9128)w0L2EI2x2+-(57128)w0LEI6x3+w0EI124x4-124(x-L2)4u(x-L2)
Simplificando
yx=9w0L2256EIx2-19w0L256EIx3+w0EI124x4-124(x-L2)4u(x-L2)
Factorizando
yx=9w0L2256EIx2-19w0L256EIx3+w024EIx4-x-L24ux-L2
A continuación se observa la representación gráfica de la deflexión de la viga, con el caso particular de w0=48EI, I=1kgm2 con la viga de granito y un metro de largo.
Bibliografía:
Dennis G. Zill y Michael R. Cullen. Ecuaciones diferenciales con problemas con valores en la frontera. 7th ed.; 2009.
F.P. Beer & F- Russell. . Mecánica de materiales. 2nd ed.; 1982