Trabajo de Investigacin Matemtica Aplicaciones de las integrales mltiples en el anlisis de elementos estructurales
Autores: Vctor Araya Christian Eltit Blas Larran Mattias Manent Pablo Ramos Profesor: Eduardo Saez Fecha: 19-11-10
ndice
Introduccin y marco terico........................................................................................................................ 3 Integrales mltiples ...................................................................................................................................... 5 Aplicaciones de las integrales dobles............................................................................................................ 6 Mecnica Clsica ........................................................................................................................................... 8 Objetivos de la investigacin ........................................................................................................................ 9 Procedimiento............................................................................................................................................... 9 Casos de perfiles a estudiar ........................................................................................................................ 10 Distribuciones de carga a estudiar .............................................................................................................. 15 Calculo de Esfuerzos ................................................................................................................................... 19 Conclusiones ............................................................................................................................................... 23 Bibliografa .................................................................................................................................................. 23
Introduccin y marco terico
Optimizar y abaratar costos es uno de los muchos objetivos de las empresas del mundo actual, en particular en el rea de la ingeniera. Una de las herramientas que ms se utiliza para lograr estos objetivos es el clculo, que tras un anlisis detallado de un problema, entrega una solucin adecuada a ste. En la construccin cobra gran importancia lo anterior, ya que los costos son abismantes y hay claramente muchas reas donde estos pueden ser reducidos, sin dejar de lado los mltiples factores de seguridad con los que son regulados. Las vigas son elementos estructurales fundamentales en la construccin, donde cumplen mltiples funciones y ser su fin lo que la regular en trminos de sus dimensiones, material del que este hecha y forma, porque esas propiedades influirn en la resistencia de distintas cargas a las cuales son sometidas.
Diagrama 1: Viga simplemente apoyada, sometida a carga puntual.
La viga puede estar sometida a varios estados de esfuerzos de flexin y traccin, que le darn la forma para resistirlos, ya que no queremos que nuestra viga colapse y genere una catstrofe. Sin embargo estos estados hacen trabajar de distinta forma a la viga, habiendo partes de ellas donde no est sometida a esfuerzos, o son muy pequeos, con lo que se puede ahorrar material eliminando esas partes ya que no son tiles y aumentan los costos de la estructura sin afectar los rendimientos de esta. Una viga puede fallar por muchas razones, y cada una de ellas est regulada por las normas chilenas NCh. Estas fallas pueden ser por traccin (estiramiento de la viga), compresin (cuando se comprime la viga), flexin, deflexin (curvatura), pandeo, entre otros.
3
Nos centraremos en la falla por flexin, cuya norma considera esfuerzos normales y de corte admisibles. Estos esfuerzos vienen dados por las siguientes formulas:
Frmula 1
Frmula 2
La frmula 1 (Navier) representa el esfuerzo normal, que es el de traccin y compresin de la viga, el de la ecuacin representa el momento de inercia del perfil de la viga en torno al eje z, M es el momento interno de la viga e y es la distancia del eje neutro al lugar de referencia que se est calculando el esfuerzo. Veremos esto con ms detalle posteriormente. En la Frmula 2 (Jourasky) lo que se calcula es el esfuerzo de corte, y sus variables son: V es el corte interno de la viga. Q es el primer momento, b es el ancho del perfil donde se est trabajando y finalmente el I es el momento de inercia de la viga. Para el clculo de las variables como lo son el primer momento o el momento de inercia es til usar las integrales mltiples. La explicacin de cmo aplicarlas se explica a continuacin.
Integrales mltiples
Consideremos que P={ } es una particin, por una familia de curvas, de una regin , cerrada y acotada. Diremos que: d = supremo {d . } es el dimetro de entre sus puntos mas alejados, donde
Sea S(P,f)= una Doble suma de Riemann de una funcin f en la particin P , donde es un punto arbitrario de cada . Luego se define (P) como la norma de la particin y analizando el caso particular cuando (P)0 se tiene que: d =0 , 0 Tiende a un punto. La regin quedara dividida en regiones infinitesimales, y adems se puede plantear que: Si el existe, entonces =
Integrales Iteradas: Sea f : una funcin real y continua , definida en una regin =*a,b+x*c,d+, la integral iterada de la funcin f sobre puede escribirse como dos integrales en sentido unidimensional de Newton. =
Teorema de Fubini: Sea cerrada y acotada f: , Entonces: 1)f es integrable en el sentido Riemann 2)La integral de Riemann coincide con las integrales iteradas =
5
Aplicaciones de las integrales dobles
Las integrales mltiples pueden tener distintas interpretaciones, tanto fsicas como geomtricas, dependiendo del sentido que se le d a la suma. rea de una figura plana: En el caso que se considere f(x,y)=1 , la integral quedara representada como: = Masa: Considerando ahora que la funcin f representa densidad variable , f(x,y)=(x,y), y que = , la integral para calcular la masa de una lmina plana que ocupa una regin D, M(D), queda representada como: = = , = lo cual represente el area de una regin D.
Para el caso de la masa de un slido en el espacio, se integra sobre una regin R considerando como densidad volumtrica variable, entonces la expresin: = representa la masa de un slido en el espacio.
Centro de masa: Consideramos el 1 momento de inercia como: =M , donde y es la distancia de una particula de masa dm , despejando , se obtiene la segunda coordenada del centro de masa.
= = ( ) son las coordenadas del centro de masa. En el caso que la densidad no sea uniforme la expresin para , queda expresada como: = , anlogamente =
De la misma manera si quiere analizar el centro de masa de un slido tridimensional de densidad no uniforme, se tendr la siguiente igualdad: =M las coordenadas del centro de masa en el espacio: = = =
Momento de inercia: La integral que representa al momento de inercia, de una placa D con eje de rotacin = dm O bien, en el caso que la placa no tenga densidad uniforme la expresin para calcular el momento de inercia estar dada por: =
Mecnica Clsica
Adems del clculo, requerimos en nuestro anlisis posterior un conocimiento de la mecnica clsica: Para analizar la resistencia del material, hay que introducirnos en la viga y estudiar su comportamiento interno Cmo se comportan las fuerzas internas a lo largo de la viga? Por equilibrio esttico, la suma de las fuerzas debe ser cero para que no se desplace, asimismo la sumatoria de los momentos debe ser cero (tendencia a rotar del objeto) ya que no debe rotar. Diagrama 2: Comportamiento interno
Para determinar las fuerzas y momentos internos, primero se corta la viga. Por accin y reaccin, y estabilidad de la viga se agrega al extremo una fuerza y un momento que son los que simulan el comportamiento del lado cortado de la viga (diagrama 2). Se repite el procedimiento a lo largo de toda la viga. Con lo anterior se crean diagramas de fuerzas internas, que son el comportamiento de cada magnitud fsica a lo largo de la viga representados en una funcin (N(x) fuerza axial, V(x) fuerza de corte y M(x) momento interno). Estas aplicaciones sin duda nos facilitan el anlisis que haremos a continuacin de los distintos perfiles de vigas y la resistencia de stos.
Objetivos de la investigacin
La siguiente investigacin tiene como objetivo estudiar el comportamiento de una viga bajo distintos estados de esfuerzo. Estudiaremos para esto distintas cargas y perfiles que son recurrentes en el rea de la ingeniera civil. Como se trata de un rea profesional en que los recursos deben ser usados con de manera eficiente, buscaremos un perfil de viga ptimo, econmico y resistente utilizando las integrales mltiples, y aplicando conocimientos de la mecnica de slidos. Considerando lo anterior, la investigacin consistir de una aplicacin matemtica concreta al rea del clculo de elementos estructurales, como una aproximacin de lo que se realiza en las oficinas de clculo de los ingenieros civiles.
Procedimiento
El procedimiento consistir en calcular las distintas variables por separado, es decir, los momentos y fuerzas internas y los momentos de inercia de las vigas, y luego aplicar las formulas mencionadas anteriormente para obtener el valor de los esfuerzos en los elementos.
Casos de perfiles a estudiar
Uno de los puntos ms importantes en el anlisis de elementos estructurales es la determinacin del momento de inercia de las secciones transversales de estos. Aquello se debe a que este valor influir en gran parte de los clculos de los esfuerzos en el material.
Como vimos anteriormente, el momento de inercia de una figura plana viene dado por la integral donde r es la distancia desde el elemento diferencial de rea al eje respectivo. Cuando se trata de flexin simtrica, debemos determinar este momento con respecto al eje z arbitrario, el cual ubicamos en el centro geomtrico del elemento. A continuacin analizaremos tres casos de figuras simtricas que son recurrentes en elementos estructurales, y utilizaremos integrales mltiples para determinar sus momentos de inercia. Para este efecto, seleccionamos las siguientes figuras: (1) una viga rectangular, (2) una viga doble T, (3) un tubo utilizado como viga estructural. Diremos que los tres elementos son del mismo material, y sus dimensiones estn dadas por las mismas constantes.
Figura (1)
Figura (2)
Figura (3)
10
Luego de estas consideraciones, podemos comenzar a calcular los momentos de inercia. Para la figura (1), vemos que el centro geomtrico, y por consiguiente el eje z, se ubica en la mitad de la altura, por lo que la distancia al eje viene dada por la coordenada y, por lo que la distancia al cuadrado al eje ser . Entonces ubicamos la figura en el plano zy de la siguiente forma:
Por consiguiente, la integral queda: Falta fijar los lmites de integracin. Integraremos utilizando coordenadas cartesianas, por lo que el elemento diferencial de rea ser , y los limites deben cubrir toda el rea de la seccin transversal. Entonces, para la variable z los lmites van desde 0 a b, y para la variable y, desde -h/2 a h/2. Finalmente la integral queda (por el teorema de Fubini): ( ) ( ( )
11
Para la figura (2), nuevamente por simple inspeccin ubicamos el eje z en la mitad de la altura de la viga. Nuevamente se cumple que la distancia al cuadrado viene dada por la coordenada y, por lo que el momento de inercia est determinado por: No obstante, el caso de los lmites de integracin es algo ms complejo que el anterior. Debemos considerar la viga como tres elementos separados y proceder a plantear tres integrales para cada uno de estos. Es decir, calcular el momento de inercia del alma (parte central) y luego el de las alas (parte superior e inferior de una viga).
.-
El alma tiene dimensiones (h-2b) de alto y b de ancho. Entonces para el alma, tenemos que el momento de inercia con respecto al eje z es
( )
( )
( (
))
12
Para las alas, el nico problema consiste en determinar los lmites de integracin. Sin embargo esto es ms simple de lo que parece, ya que las referencias no cambian, basta con utilizar la altura donde comienzan y la donde terminan los elementos rectangulares.
Entonces para el ala superior, el momento de inercia es: ( ) ( ( )
( )
Y para el ala inferior:
( )
( (
Debido a que usamos el mismo eje de referencia para el clculo de las integrales, y como las integrales en realidad son bsicamente sumas de cosas, podemos decir que el momento de inercia total es simplemente la suma de los momentos de inercia de los tres elementos que conforman la viga. TOTAL:
13
El caso de la figura (3), como se trata reas circulares en el plano, conviene claramente analizarlo por medio de coordenadas polares.
Para saber cul ser la distancia r hasta el eje z, el cual coincida con la coordenada y en cartesianas, debemos considerar la transformacin que lleva coordenadas polares a rectangulares: { Entonces la distancia al cuadrado desde un punto al eje z en coordenadas polares es
Los lmites de integracin vendrn dados por los radios de las circunferencias, que determinan el rea cubierta. En este caso van desde Entonces por Fubini, la integral queda: hasta . La variable debe cubrir desde 0 a 2 .
Ntese que para pasar a coordenadas polares en la integral, se debi multiplicar por el jacobiano de la transformacin.
14
Distribuciones de carga a estudiar
Estudiaremos las 3 distribuciones de cargas ms comunes en las vigas simplemente apoyadas. Estas son, viga sometida a cargas puntuales, las vigas con carga distribuida por toda su extensin, y las vigas sometidas a cargas distribuidas pero en segmentos de la viga. El largo de la viga es de 3 [m] La primera muestra una viga con dos carga puntuales a 0,5[m] del centro, esto suele presentare en las vigas transversales de las pasarelas. Siendo el corte mximo 1500[kgf] y el momento mximo 1500[kgf x m]
1500 [kgf]
1500 [kgf]
Diagrama de Corte:
Corte 1500
1000
500
x
0.5 500 1.0 1.5 2.0 2.5 3.0
1000
1500
15
Diagrama de Momento:
Momento
1500
1000
500
0 0.0
x
0.5 1.0 1.5 2.0 2.5 3.0
La segunda muestra una viga con una carga distribuida de 3000[kgf/m] en segmentos a 0,5 [m] y 0,5 [m] a del centro, esto suele presentarse, cuando las vigas sostiene otras viga, por ejemple en puentes. Siendo el corte mximo 1500[kgf] y el momento mximo 1125 [kgf x m].
3000[kgf/m]
3000[kgf/m]
16
Diagrama de Corte:
Corte 1500
1000
500
x
0.5 500 1.0 1.5 2.0 2.5 3.0
1000
1500
Diagrama de Momento:
Momento 1400 1200 1000 800 600 400 200 0 0.0
x
0.5 1.0 1.5 2.0 2.5 3.0
17
La tercera muestra una viga que soporta una carga uniforme distribuida por toda su extensin, por ejemplo cuando una viga sostiene una losa. Siendo el corte mximo 1500[kgf/m] y el momento mximo 1125 [kgf x m].
1000[kgf/m]
Diagrama de Corte:
Corte 1500
1000
500
x
0.5 500 1.0 1.5 2.0 2.5 3.0
1000
1500
Diagrama de Momento:
Momento 1400 1200 1000 800 600 400 200 0 0.0
x
0.5 1.0 1.5 2.0 2.5 3.0
18
Calculo de Esfuerzos
Para calcular los esfuerzos en las vigas, se consideraran nueve casos. Dentro de cada tipo de distribucin de carga ms comn, se analizaran los tres perfiles mencionados anteriormente. Para hacer el anlisis de resistencia de materiales se trabaja con los esfuerzos mximos que actan sobre el elemento en estudio, as se podr concluir si la viga en este caso, puede o no ser sometida a una cierta distribucin de carga. Sea: V M
Corte mximo que acta sobre la viga. Momento mximo que acta sobre la viga. Momento de inercia de la seccion respecto al eje z. Primer momento de area respecto al eje z. Altura maxima de la seccion transversal de la viga. Ancho de la seccion transversal de la viga donde actua el esfuerzo de corte maximo. Esfuerzo axial maximo en la viga. Esfuerzo de corte maximo en la viga.
Las dimensiones b y h utilizadas anteriormente para calcular los momentos de inercia seran utilizados con valores de vigas comunes. [ ] [ ] Se utilizaran las formulas de Navier y de Jouraski para obtener tanto del Esfuerzo Axial maximo, como el Esfuerzo de corte maximo.
Caso 1.1 Viga con cargas puntuales / Seccion transversal Rectangular [ [ [ ] ] ] [ [ [ ] ] ]
19
Caso 1.2 Viga con cargas puntuales / Seccion transversal Doble T [ [ [ ] ] ] [ [ ] ] [ ]
-Caso 1.3 Viga con cargas puntuales / Seccion transversal anular [ [ [ ] ] ] [ [ ] [ ] ]
Caso 2.1 Viga con carga distribuida por segmentos / Seccion transversal Rectangular [ [ [ ] ] ] [ [ ] ] [ ]
20
Caso 2.2 Viga con carga distribuida por segmentos / Seccion transversal Doble T [ [ [ ] ] ] [ [ ] ] [ ]
Caso 2.3 Viga con carga distribuida por segmentos / Seccion transversal anular [ [ ] ] [ ] [ [ ] [ ] ]
Caso 3.1 Carga distribuida uniformemente a lo largo de la viga/ Seccion transversal Rectangular [ [ [ ] ] ] [ [ ] ] [ ]
21
Caso 3.2 Carga distribuida uniformemente a lo largo de la viga / Seccion transversal Doble T [ [ [ ] ] ] [ [ ] ] [ ]
Caso 3.3 Carga distribuida uniformemente a lo largo de la viga / Seccion transversal anular [ [ ] [ ] ] [ [ ] [ ] ]
[ .
22
Conclusiones
En conclusin, podemos afirmar que se pudo comparar de una forma efectiva los distintos tipos de perfil seleccionados, sometidos a cargas recurrentes en el rea de la ingeniera civil. En la seccin anterior, los resultados concluyeron que el perfil que mejor resista las cargas, produciendo un menor esfuerzo en el material, fue el perfil doble T, el cual precisamente es uno de los ms usados actualmente en la construccin. Le sigue el perfil anular, y luego el perfil rectangular. Si considersemos vigas del mismo material y sometidas a la misma carga, sta sera la preferencia. No obstante, estos resultados no son del todo concluyentes. Por un lado, se debe analizar tambin, en estudios futuros, qu perfil de viga es el que requiere menor material en su construccin, ya que el factor econmico es muy importante en el sector ingenieril. La viga doble T resiste mejor a los esfuerzos, pero usa considerablemente ms material que sus contrapartes rectangular y anular. Tambin se deben considerar otros factores relativos a la construccin. Por ejemplo, un buen trmino medio sera el perfil anular, que resiste bien los esfuerzos y no representa un gasto excesivo de material. Sin embargo, la forma del perfil trae dificultades a la hora de hacer las uniones en los extremos, por lo que estas vigas se utilizan principalmente cuando no es necesario efectuar estas uniones. Considerando lo anterior, el estudio realizado consisti en un anlisis aproximado de lo que se realiza en las oficinas de clculo, ya que en estos lugares se deben considerar muchos otros factores que por distintas razones no fue posible incluir en ste compilado. No obstante lo anterior, la realizacin de este artculo requiri por parte de los autores un manejo ptimo de los temas tratados, tanto en el rea de las matemticas cmo de la mecnica de slidos, por lo que cumpli su cometido, el cual era acercar a los futuros profesionales al rea en el cual eventualmente se desarrollarn, y mostrar que la matemtica y el clculo estn presentes en prcticamente todos los aspectos de la profesin ingenieril.
23
Bibliografa
GERE J.M, TIMOSHENKO S.P, Mecnica de Materiales, Mxico, Unin Tipogrfica Editorial, 1974 STEWART J, Clculo. Trascendentes Tempranas, Mxico, International Thomson Editores, 2007 HALLIDAY D, RESNICK R, Fsica, Mxico, Compaa Editorial Continental, 1980
Saez, Eduardo. MAT024 Matemticas IV, Segundo Semestre 2010, Apuntes de Clases. Tobar, Ren. CIV132 Resistencia de Materiales, Segundo Semestre 2010, Apuntes de Clases.
24
