12.07.2015 Views

Placas Planas Rectangulares - Universidad Nacional de La Plata

Placas Planas Rectangulares - Universidad Nacional de La Plata

Placas Planas Rectangulares - Universidad Nacional de La Plata

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Facultad <strong>de</strong> Ingeniería<strong>Universidad</strong> <strong>Nacional</strong> <strong>de</strong> <strong>La</strong> <strong>Plata</strong>ESTRUCTURAS IIIPara alumnos <strong>de</strong> la carrera <strong>de</strong> Ingeniería Aeronáutica y Mecánica <strong>de</strong> la UNLPPLACAS PLANAS RECTANGULARESDE ESPESOR DELGADOAutores:Ing. Alejandro J. PatanellaIng. Marcos D. Actis-2008-


Estructuras IIIPLACAS PLANAS RECTANGULARES DE ESPESORDELGADOIntroducción - Hipótesis.<strong>La</strong>s placas <strong>de</strong> espesor pequeño representan un elemento estructural muy común eimportante en estructuras <strong>de</strong> uso aeroespacial ya que las unida<strong>de</strong>s más gran<strong>de</strong>s seencuentran cubiertas con este tipo <strong>de</strong> paneles.Para el <strong>de</strong>sarrollo <strong>de</strong> la ecuación general <strong>de</strong> una placa sometida a esfuerzos <strong>de</strong>flexión puros se consi<strong>de</strong>ran las siguientes hipótesis simplificativas:• Se limita al caso <strong>de</strong> que el material que las compone sea homogéneo, isotrópico ycompletamente elástico.• <strong>La</strong> placa no experimenta variaciones <strong>de</strong> espesor <strong>de</strong>bido a la <strong>de</strong>formación (σ z = 0).• Se consi<strong>de</strong>ra valida la hipótesis <strong>de</strong> <strong>de</strong>formaciones planas <strong>de</strong> Bernoulli (las normales a lasuperficie media se conservan normales a la superficie <strong>de</strong>formada).• <strong>La</strong> flecha para cualquier punto <strong>de</strong> la placa es muy pequeña con respecto a su espesor (w


Estructuras IIIPara pequeños <strong>de</strong>splazamientos la curvatura <strong>de</strong> la superficie media pue<strong>de</strong> hallarseen forma aproximada omitiendo las potencias <strong>de</strong> primer or<strong>de</strong>n ( ∂ w∂x , ∂ w). De esta forma, la∂y curvatura <strong>de</strong>l la superficie media en planos paralelos a xz e yz respectivamente es:21 w=− ∂ 21 w2y =− ∂ 2R ∂xR ∂yxy<strong>de</strong>;También se pue<strong>de</strong> cuantificar la torsión <strong>de</strong> la superficie media <strong>de</strong> la placa a través1Rxy2w=− ∂ ∂∂ xyA través <strong>de</strong> la curvatura y la torsión se pue<strong>de</strong>n expresar las <strong>de</strong>formaciones queexperimenta la placa. Dichas <strong>de</strong>formaciones en su expresión general queda <strong>de</strong>finida como;εx∂u= εy∂x∂v=∂y(<strong>de</strong>formaciones lineales en función <strong>de</strong> los corrimientos uen al dirección x y v en la dirección y)εxy∂u∂v= γxy= +∂y∂x(<strong>de</strong>formación angular, distorsión, en el plano xy)En el caso <strong>de</strong> flexión pura <strong>de</strong> barras prismáticas una solución rigurosa se obtiene apartir <strong>de</strong> suponer que las secciones transversales se mantienen planas luego <strong>de</strong> flexionarsey rotar <strong>de</strong> forma tal <strong>de</strong> seguir siendo perpendiculares a su eje neutro. Combinando estasuposición en las dos direcciones perpendiculares se pue<strong>de</strong> obtener la ecuación generalpara la flexión pura <strong>de</strong> placas.Figura 2a.Figura 2b.Página 2 <strong>de</strong> 17


Estructuras III<strong>La</strong> figura 2a representa una fina placa rectangular cargada con un momento flexoruniformemente distribuido (Mx, My) por unidad <strong>de</strong> longitud en sus bor<strong>de</strong>s. Se consi<strong>de</strong>ranpositivos los momentos que generen un esfuerzo <strong>de</strong> compresión en la cara superior <strong>de</strong> laplaca y uno <strong>de</strong> tracción en la superficie inferior. En la figura 2b se representa un elementorectangular tomado <strong>de</strong> la placa con unas dimensiones dx, dy y t, siendo esta última muypequeña comparada con las <strong>de</strong>más dimensiones. Debido a que este elemento esta extraído<strong>de</strong> la placa <strong>de</strong> la figura 2a, la misma solicitación que existía sobre la placa pue<strong>de</strong> serextendida a dicho elemento. Los lados laterales <strong>de</strong>l elemento se mantienen planos durantela flexión y la rotación <strong>de</strong>bido a la solicitación existente, <strong>de</strong> forma tal <strong>de</strong> que permanezcannormales a la superficie media <strong>de</strong>flectada; y <strong>de</strong>bido a la simetría existente la superficiemedia no sufre ningún cambio en sus dimensiones y es por eso que la misma es lasuperficie neutra o media <strong>de</strong> la placa.En función <strong>de</strong>l mecanismo <strong>de</strong> <strong>de</strong>formación anteriormente presentado los<strong>de</strong>splazamientos en las direcciones x, y y z se pue<strong>de</strong>n hallar con el siguienteprocedimiento:Un punto B <strong>de</strong> la superficie media (no <strong>de</strong>formada) se <strong>de</strong>splazó a la posición B 1(<strong>de</strong>formada) en la dirección <strong>de</strong>l eje z en una magnitud W. El elemento <strong>de</strong> superficie dzdyasociado a este punto ha rotado un ángulo igual a la pendiente <strong>de</strong> la superficie mediaflexionada en una dirección tal que este se mantiene normal esta. (ver figura 3)Figura 3.Este ángulo <strong>de</strong> rotación para pequeños <strong>de</strong>splazamientos es igual a ∂ w. De esta∂x forma los distintos <strong>de</strong>splazamientos en cada una <strong>de</strong> las direcciones en un punto a unadistancia z <strong>de</strong> la superficie media se pue<strong>de</strong>n escribir como;u =− z ∂w x∂xu z w y=− ∂ ∂y(los signos - indican <strong>de</strong>splazamientos negativosen x e y para <strong>de</strong>splazamientos positivos en z)uz = w( x, y)Página 3 <strong>de</strong> 17


<strong>La</strong>s <strong>de</strong>formaciones correspondientes se pue<strong>de</strong>n hallar a partir <strong>de</strong>;Estructuras IIIεxy2∂ux∂ w zεx= =− z2=∂x∂xRx∂u2y ∂ w zεy= =− z2=∂y∂yRy∂u∂u2x y ∂ w= γxy= + = 2z=∂x∂y∂∂ xy2zRxyYa que para el estudio <strong>de</strong> las placas se consi<strong>de</strong>ró un estado plano <strong>de</strong> tensiones, através <strong>de</strong> la ley <strong>de</strong> Hook se obtiene, para un material cuyo modulo <strong>de</strong> elasticidad sea E ycoeficiente <strong>de</strong> Poisson µ, las siguientes expresiones para las tensiones en el plano <strong>de</strong> laplaca;Ez ⎛ µ ⎞2∂σµ ∂ 21Ez ⎛ w w⎞x= ⎜ + ⎟ =− ⎜ + ⎟( 1 − µ2 ) ⎝ RxRy⎠ ( 1 − µ2 22) ⎝ ∂x∂y⎠Ez ⎛ µ ⎞2∂σµ ∂ 21Ez ⎛ w w⎞y= ⎜ + ⎟ =− ⎜ + ⎟( 1−µ2 ) ⎝ RyRx⎠ ( 1−µ2 22) ⎝ ∂y∂x⎠Estas tensiones normales están linealmente distribuidas a lo largo <strong>de</strong>l espesor <strong>de</strong> laplaca y su resultante, <strong>de</strong>bido a las condiciones <strong>de</strong> equilibrio existente, <strong>de</strong>be ser igual a Mxy My respectivamente, es <strong>de</strong>cir;∫h2σ xh2h2∫σ yh−2−zdydz=M dyzdxdz=M dxSustituyendo la expresión para las tensiones nos queda;yxMxh22222 3E ⎛∂w w⎞E ⎛ w w⎞2h=−2⋅ ⎜ 2+ µ ∂ ∂2 ⎟ ∫ z dz=−2⋅ ⎜ 2+ µ ∂ 2 ⎟ ⋅( 1−µ ) ⎝ ∂x∂ y ⎠ h ( 1−µ ) ⎝ ∂x∂y⎠ 12−2De la misma forma se opera para M y y se obtiene2⎛∂w w⎞Mx =−D⋅ ⎜ + µ ∂ 222 ⎟⎝ ∂x∂y⎠2wM =− wyD ⋅ ⎛∂⎜⎝ y+⎞µ ∂ 222 ⎟∂ ∂x⎠(Momentos por unidad <strong>de</strong> longitud para laflexión pura sobre la placa)Página 4 <strong>de</strong> 17


3EhSiendo D la rigi<strong>de</strong>z a flexión <strong>de</strong> la placa dada por D =2 .12( 1 − )µEstructuras IIIA<strong>de</strong>más <strong>de</strong> los momentos flectores existe sobre la placa un momentouniformemente distribuido <strong>de</strong> torsión a lo largo <strong>de</strong> las lados <strong>de</strong> la misma dados por M xy yM yx , cada uno <strong>de</strong> ellos <strong>de</strong>berá ser igual a la resultante <strong>de</strong> las fuerzas <strong>de</strong> corte existentes a lolargo <strong>de</strong> las lados <strong>de</strong>l elemento, es <strong>de</strong>cir,τxyM dx = ∫τ z dxdzxyh−2h2xy22 Gz ∂ w= τyx= = 2 GzR ∂∂ xyxyM dy = ∫τ z dydzyxh−2h2yxM = M = D( 1 − )xyyx2wµ ∂ ∂∂ xy(Momento por unidad <strong>de</strong> longitud para torsiónpura sobre la placa <strong>de</strong>bida a la flexión)Ecuación diferencial <strong>de</strong> la superficie <strong>de</strong>formada.En el <strong>de</strong>sarrollo anterior <strong>de</strong> la teoría <strong>de</strong> las pequeñas <strong>de</strong>formaciones para placas<strong>de</strong>lgadas se utilizó como hipótesis que en los límites <strong>de</strong> la placa sus bor<strong>de</strong>s se pue<strong>de</strong>nmover libremente en el plano <strong>de</strong> la misma. De esta forma, las fuerzas reactivas en losbor<strong>de</strong>s <strong>de</strong>bido a los vínculos son normales a la placa. Con el mismo criterio se pue<strong>de</strong>n<strong>de</strong>spreciar las <strong>de</strong>formaciones <strong>de</strong> la superficie media durante la flexión <strong>de</strong> la placa.qFigura 4.Página 5 <strong>de</strong> 17


Estructuras IIIViendo la figura anterior (fig. 4) se estudia un elemento dxdy <strong>de</strong> la superficiemedia, el cual posee distribuido en sus bor<strong>de</strong>s los momentos Mx, My y Mxy. Estosmomentos son los resultantes <strong>de</strong> las distribuciones lineales <strong>de</strong> tensiones <strong>de</strong>bido a flexión ytorsión a lo largo <strong>de</strong>l espesor <strong>de</strong> la placa. Si la placa es solicitada a través <strong>de</strong> cargasexternas normales a la superficie media a <strong>de</strong>más <strong>de</strong> los momentos anteriormenteexpresados existirá también un esfuerzo <strong>de</strong> corte vertical Qx y Qy, que actúan sobre loslados <strong>de</strong>l elemento.Qx=∫σ dzh−2h2xzQy=∫σ dzh−2h2yzLlamando como q a la carga transversal por unidad <strong>de</strong> área que actúa normal a lacara superior <strong>de</strong> la placa y consi<strong>de</strong>rando el equilibrio <strong>de</strong> fuerzas en la dirección <strong>de</strong>l eje z<strong>de</strong>l elemento se tiene,∂Q∂∂x dxdyQxy+∂y dydx + q dxdy = 0 o bien∂Qx∂x∂Qy+ + q = 0∂yTomando el equilibrio <strong>de</strong> momentos existentes en la dirección <strong>de</strong>l eje x,∂Mxy∂∂x dxdy My−∂y dydx + Q dxdy y= 0 o bien∂M∂xxy∂My− + Qy=0∂yRealizando el mismo procedimiento para la dirección y,∂M∂yxy∂Mx+ − Qx= 0∂xEliminando Qx y Qy <strong>de</strong> las ecuaciones anteriores se obtienen las relaciones <strong>de</strong>equilibrio entre momentos,2∂ M2∂x22∂ M ∂ M+2− 2 =− q∂y∂∂ xyx y xyPara representar esta ecuación en términos <strong>de</strong> la <strong>de</strong>flexión w <strong>de</strong> la placa, seconsi<strong>de</strong>ra que la ecuación <strong>de</strong>sarrollada para el caso <strong>de</strong> la placa sometida a una solicitación<strong>de</strong> flexión pura, es aproximadamente igual a la que correspon<strong>de</strong>ría al caso <strong>de</strong> placascargadas lateralmente, es <strong>de</strong>cir2⎛∂wMx =−D⋅ ⎜ 2+⎝ ∂xw⎞µ ∂ 22 ⎟∂y⎠Página 6 <strong>de</strong> 17


2wM =− yD ⋅ ⎛∂⎜⎝ ∂y+ 2w⎞µ ∂ 22 ⎟∂x⎠Estructuras IIIEsta consi<strong>de</strong>ración trae como consecuencia <strong>de</strong>spreciar los efectos sobre la flexión<strong>de</strong> las fuerzas <strong>de</strong> corte y <strong>de</strong> la tensión <strong>de</strong> compresión σ z . Reemplazando se tiene,∂4∂xw444∂ w ∂ w+4+ 22 2=∂y∂x∂yqDEcuación <strong>de</strong>Germain-<strong>La</strong>grange<strong>La</strong> dificultad existente en el análisis <strong>de</strong> la flexión en placas se concentra en laintegración <strong>de</strong> la ecuación anterior en términos <strong>de</strong> w.<strong>La</strong>s fuerzas <strong>de</strong> corte expresadas en términos <strong>de</strong> los <strong>de</strong>splazamientos se pue<strong>de</strong>nexpresar como,QQxy∂M2 2xy ∂Mx ∂ ⎛∂w ∂ w⎞= + =−D ⋅ ⎜ + ⎟∂y∂x∂ x2 2⎝ ∂x∂y⎠∂My∂M2 2xy ∂ ⎛∂w ∂ w⎞= − =−D ⋅ ⎜ + ⎟∂y∂x∂ y2 2⎝ ∂x∂y⎠El análisis realizado anteriormente es suficiente para plantear soluciones acualquier problema específico. El método general <strong>de</strong> trabajo consistirá, entonces, en hallarsoluciones aproximadas para la ecuación diferencial <strong>de</strong> cuarto or<strong>de</strong>n que satisfaga lascondiciones <strong>de</strong> bor<strong>de</strong> y las condiciones <strong>de</strong> cargas, <strong>de</strong>ntro <strong>de</strong> errores apropiadamenteacotados.Página 7 <strong>de</strong> 17


Estructuras IIIMÉTODOS PARA HALLAR UNA SOLUCIÓN SATISFACTORIA.Flexión <strong>de</strong> placas rectangulares.<strong>La</strong> ecuación diferencial general para la flexión dada anteriormente se resolveráutilizando métodos exactos como el <strong>de</strong> Serie <strong>de</strong> Fourier y el método alternativo <strong>de</strong> Levi,también se hará mención mas a<strong>de</strong>lante <strong>de</strong> métodos aproximados a través <strong>de</strong> diferenciasfinitas como ser el método <strong>de</strong> Marcus.Para cada uno <strong>de</strong> los métodos se hará el tratamiento para placas rectangulares condistintas condiciones <strong>de</strong> contorno, es <strong>de</strong>cir, <strong>de</strong> vinculación <strong>de</strong> la placa al mediocircundante. <strong>La</strong>s condiciones <strong>de</strong> contorno se pue<strong>de</strong>n clasificar <strong>de</strong> la siguiente forma:• Bor<strong>de</strong> empotrado:<strong>La</strong> <strong>de</strong>flexión a la largo <strong>de</strong>l bor<strong>de</strong> empotrado es cero y la tangente al plano <strong>de</strong> lasuperficie media <strong>de</strong>flectada es horizontal. Entonces si el eje x coinci<strong>de</strong> con el bor<strong>de</strong>empotrado tendremos;( ) w y==00⎛∂w⎞⎜ ⎟ = 0⎝ ∂ ⎠y y =0• Bor<strong>de</strong> simplemente apoyado:<strong>La</strong> <strong>de</strong>flexión a lo largo <strong>de</strong>l bor<strong>de</strong> simplemente apoyado es cero y el momento flexorparalelo a este lado también será nulo. Entonces si el eje x coinci<strong>de</strong> con el bor<strong>de</strong>simplemente apoyada tendremos;( ) w y==00 ( My )y=02 2⎛∂w ∂ w⎞= ⎜ 2+2 ⎟ = 0⎝ ∂x∂y⎠y=0• Bor<strong>de</strong> libre:El momento flexor, el momento torsor y la fuerza <strong>de</strong> corte a lo largo <strong>de</strong> estelado es nula. Si el lado libre coinci<strong>de</strong> con la línea recta correspondiente a un x=aL setiene,( M x )x=a= 0 ( M xy )x=a= 0 ( Q x )x=aComo ha sido probado por Kirchoff, dos condiciones <strong>de</strong> bor<strong>de</strong> son solo necesariaspara encontrar una única solución al problema <strong>de</strong> flexión. Aunque el también <strong>de</strong>mostróque las dos últimas condiciones <strong>de</strong> M xy y Q x pue<strong>de</strong>n ser reemplazadas por una sola, es<strong>de</strong>cir,= 03⎛∂w⎜ 3⎝ ∂x( µ )∂ 3w ⎞+ 2−2 ⎟ = 0∂∂ xy ⎠x=aPágina 8 <strong>de</strong> 17


Estructuras IIIY finalmente expresando la condición <strong>de</strong> M x en términos <strong>de</strong> w las condicionesfinales para un bor<strong>de</strong> libre queda como,2⎛∂w⎜ 2⎝ ∂x2∂ w⎞+ ⎟ =∂y⎠20x=a3⎛∂w⎜ 3⎝ ∂x( µ )∂ 3w ⎞+ 2−2 ⎟ = 0∂∂ xy ⎠x=aPágina 9 <strong>de</strong> 17


Estructuras IIIDistintos Ejemplos <strong>de</strong> <strong>Placas</strong> con distintas condiciones <strong>de</strong>Bor<strong>de</strong>.<strong>Placas</strong> rectangulares simplemente apoyadas.Se estudiará una placa rectangular <strong>de</strong> dimensiones a y b con los ejes x e y colocados<strong>de</strong> forma tal que su origen sea este en un vértice <strong>de</strong> la placa como se ve en la siguientefigura. <strong>La</strong> placa se encuentra simplemente apoyada a lo largo <strong>de</strong> sus cuatro lados y estasolicitada por una carga distribuida uniformemente q=f(x,y).xabqyqFigura 6.• Solución a través <strong>de</strong> series dobles <strong>de</strong> Fourier.Siempre se pue<strong>de</strong> expresar una función f(x,y) como una serie trigonométrica doble(Fourier), es <strong>de</strong>cir,don<strong>de</strong>:∞ ∞⎛q f x y a sin m π x ⎞ ⎛= =sin n π y ⎞( , ) ∑∑ mn⎜ ⎟ ⎜ ⎟⎝ a ⎠ ⎝ b ⎠m=1 n=1amn4=ab∫0a∫0b⎛f x y sin m π( , ) ⎜⎝ ax ⎞ ⎛sin n π y ⎞⎟ ⎜ ⎟dxdy⎠ ⎝ b ⎠<strong>La</strong>s condiciones <strong>de</strong> bor<strong>de</strong> son:o bienW=0 y M x =0 en x=0 y x=aW=0 y M y =0 en y=0 y y=bPágina 10 <strong>de</strong> 17


Estructuras IIIW=0 yW=0 y2∂ w2= 0∂xen x=0 y x=a2∂ w2= 0∂yen y=0 y y=bEstas condiciones <strong>de</strong> bor<strong>de</strong> son satisfechas si se toma la siguiente expresión para el<strong>de</strong>splazamiento total,∞ ∞⎛W C sin m π x ⎞ ⎛=sin n π y ⎞∑∑ mn⎜ ⎟ ⎜ ⎟⎝ a ⎠ ⎝ b ⎠m=1 n=1De la ecuación diferencial general <strong>de</strong> la <strong>de</strong>formada <strong>de</strong> la placa4 44∂ w ∂ w ∂ w q4+4+ 22 2= , hallo cada término reemplazando la serie propuesta∂x∂y∂x∂yDanteriormente y <strong>de</strong>rivando.Esta operación se hace para po<strong>de</strong>r hallar una relación en los coeficientes a mn y C mn .El calculo <strong>de</strong> las <strong>de</strong>rivadas parciales queda como,Análogamente,∞ ∞∂w⎛mπ ⎞ ⎛mπ x⎞⎛ πCsin n y ⎞= ∑∑mn⎜⎟ cos⎜⎟ ⎜ ⎟∂xm=1 n=1 ⎝ a ⎠ ⎝ a ⎠ ⎝ b ⎠2 ∞ ∞2 2∂ wm π mπ x π2= C⎛ ⎞ ⎛ ⎞ ⎛2sin n y ⎞∑∑mn⎜−⎟ sin ⎜ ⎟ ⎜ ⎟∂xm=1 n=1 ⎝ a ⎠ ⎝ a ⎠ ⎝ b ⎠3 ∞ ∞3 3∂ wm π mπ x π3= C⎛ ⎞ ⎛ ⎞ ⎛3sin n y ⎞∑∑mn⎜−⎟ cos ⎜ ⎟ ⎜ ⎟∂xm=1 n=1 ⎝ a ⎠ ⎝ a ⎠ ⎝ b ⎠4 ∞ ∞4 4∂ w⎛mπ ⎞ ⎛mπ x⎞⎛ π4C4sin n y ⎞= ∑∑mn⎜⎟ sin ⎜ ⎟ ⎜ ⎟∂xm=1 n=1 ⎝ a ⎠ ⎝ a ⎠ ⎝ b ⎠4 ∞ ∞4 4∂ w⎛nπ ⎞ ⎛mπ x⎞⎛ π4C4sin n y ⎞= ∑∑mn⎜⎟ sin ⎜ ⎟ ⎜ ⎟∂ym=1 n=1 ⎝ b ⎠ ⎝ a ⎠ ⎝ b ⎠Y la <strong>de</strong>rivada cruzada queda como,4∂ w∂x∂y2 2∞ ∞2 2 2 2⎛mπ ⎞ n π mπ x π= C22sin n ymn⎜⎟ ⎛m 1 n 1 ⎝ a ⎠ ⎝ ⎜⎞ ⎛ ⎞ ⎛ ⎞∑∑⎟ sin ⎜ ⎟ ⎜ ⎟= =a ⎠ ⎝ a ⎠ ⎝ b ⎠Página 11 <strong>de</strong> 17


Quedando le ecuación general <strong>de</strong> la siguiente forma:Estructuras III∞ ∞4 4⎛ m π ⎞ m xCsin n y ∞ ∞4 4⎛ π ⎞ ⎛ π ⎞ ⎛ n π ⎞ ⎛mπ x⎞⎛∑∑mn⎜4 ⎟ sin ⎜ ⎟ ⎜ ⎟ + Csin n π y ⎞∑∑mn⎜4 ⎟ sin ⎜ ⎟ ⎜ ⎟ +m=1 n=1 ⎝ a ⎠ ⎝ a ⎠ ⎝ b ⎠ m=1 n=1 ⎝ b ⎠ ⎝ a ⎠ ⎝ b ⎠∞ ∞2 2 2 2⎛ m π ⎞ n π mπ x+ 2 Csin n π ymn⎜2 ⎟ ⎛ 2m 1 n 1 ⎝ a ⎠ ⎝ ⎜⎞ ⎛ ⎞ ⎛ ⎞∑∑⎟ sin ⎜ ⎟ ⎜ ⎟ = q= =a ⎠ ⎝ a ⎠ ⎝ b ⎠ DReemplazando la aproximación <strong>de</strong> la carga distribuida q(x,y) propuesta en serie <strong>de</strong>Fourier tenemos que operando matemáticamente se pue<strong>de</strong> <strong>de</strong>spejar C mn siendo,Cmn=a⎛D ⎜m ⎝ aπ 4 2mn2n+b22⎞⎟⎠2El <strong>de</strong>splazamiento total queda,1W =π D∞∞∑∑4 2m=1 n=1⎛m⎜⎝ a2amnn+b22⎞⎟⎠2⎛sin m π x ⎞ ⎛sin n π⎜y ⎞⎟ ⎜ ⎟⎝ a ⎠ ⎝ b ⎠Como para nuestro caso habíamos consi<strong>de</strong>rado solo una carga distribuidauniformemente, los coeficientes a mn se pue<strong>de</strong>n hallar integrando <strong>de</strong> la siguiente manera,f(x,y) = q o = cte.amnabqsin m x sin n yoqo=4 ⎛ π ⎞ ⎛ π ⎞∫ ∫ ⎜ ⎟ ⎜ ⎟ dxdy =162ab ⎝ a ⎠ ⎝ b ⎠ π mn00Don<strong>de</strong> m y n son números naturales impares, lo cual al sustituir en la expresión <strong>de</strong>l<strong>de</strong>splazamiento total tenemos que;16qoW =π D∞∑∞∑⎛sin m π x ⎞ ⎛sin n π y ⎞⎜ ⎟ ⎜ ⎟⎝ a ⎠ ⎝ b ⎠mn m 2 2⎛ n ⎞⎜ 2+2 ⎟⎝ a b ⎠6 2m=135 , , n=135 , ,Siendo la máxima <strong>de</strong>flexión en el centro <strong>de</strong> la placa, dado su valor por,Wmaxqo=16 ∞ ∞( −∑ ∑1 )6π D m=n=⎛mn⎜m 2135 , , 135 , ,2⎝ am+ n − 12n+b22⎞⎟⎠2Página 12 <strong>de</strong> 17


Estructuras IIIEsta es una serie cuya convergencia es inmediata y se pue<strong>de</strong> tomar una buenaaproximación consi<strong>de</strong>rando solamente el primer termino, lo cual para una placa cuadradase convierte enWmax444 qoa qoa≅6= 0. 04546(siendo µ = 0.3)π D hEEsta última expresión difiere <strong>de</strong> la exacta en un 2 ½ % .<strong>La</strong> expresión para el momento flector queda como,Mx16q=4πo∞∑∞∑m=135 , , n=135 , ,⎛a n m 2⎜⎝ am22n+b22⎞⎟⎠2⎛sin m π x ⎞ ⎛sin n π⎜y ⎞⎟ ⎜ ⎟⎝ a ⎠ ⎝ b ⎠<strong>La</strong>s expresiones para los momentos flectores y torsores finalizan en series cuyaconvergencia es más difícil <strong>de</strong> hallar, es por eso que si no se posee una herramienta <strong>de</strong>cálculo potente no es posible hallar estas expresiones a través <strong>de</strong>l método <strong>de</strong> Series Dobles<strong>de</strong> Fourier.• Solución a través <strong>de</strong> series simples, Método <strong>de</strong> Levy:Este método se <strong>de</strong>sarrolló solo para el caso en que la carga q o sea constante. Levysugirió una solución <strong>de</strong> la forma,∞⎛W Y y sin m π=x ⎞∑ m( ) ⎜ ⎟⎝ a ⎠m=1Don<strong>de</strong> Y m es una función <strong>de</strong> y solamente. Cada termino <strong>de</strong> la serie satisface las2∂ wcondiciones <strong>de</strong> contorno W=0 y2= 0, en x=0 y x=a. Queda por <strong>de</strong>terminar la∂x2∂ wfunción Y m tal que satisfaga las condiciones <strong>de</strong> bor<strong>de</strong> restantes, es <strong>de</strong>cir W=0 y2= 0∂yen y=0 y y=b.Una aproximación más se pue<strong>de</strong> realizar si se consi<strong>de</strong>ra una solución <strong>de</strong> la forma,Don<strong>de</strong>w= w1 + w2qowD x 4 ax 3= − 2 +24a 3( x )Página 13 <strong>de</strong> 17


Estructuras IIIEsta expresión representa la <strong>de</strong>flexión <strong>de</strong> una larga “tira” con los lados mas largosen dirección <strong>de</strong>l eje y, que se encuentra cargada uniformemente por una carga q o ysimplemente apoyada en los bor<strong>de</strong>s cortos en x=0 y x=a, y libre en los otros dos.Aunque esta ultima expresión satisface las condiciones <strong>de</strong> bor<strong>de</strong> en x=0 y x=a, elproblema estará resuelto si se encuentra la solución <strong>de</strong>,4∂ w∂x2444+ ∂ w2∂2422 20∂y+ w∂x∂y=tomando como w 2 a la solución propuesta por Levy, para que satisfaga con w 1 lascondiciones <strong>de</strong> bor<strong>de</strong> W=0 y2∂ w2= 0 en∂yy=±b/2 en la expresión supuesta para w.(ver figura 7)b/2axb/2yFigura 7.Sustituyendo se obtiene,∞∑m=1⎛⎜Y⎝IVm2 24 42 m π m π ⎞ ⎛YY sin m π x ⎞II−2 m+4 m ⎟ ⎜ ⎟ = 0aa ⎠ ⎝ a ⎠don<strong>de</strong> por simetría m solo toma los valores impares (m=1,3,5,....).Esta ecuación pue<strong>de</strong> ser satisfecha para cualquier valor <strong>de</strong> x siYIVm2 24 42 m π m πII−2Ym+4Ym= 0aaSiendo su solución general,Ymqoa(y) =D4⎡ ⎛ mπy⎞ mπy⎛ mπy⎞⎢Amch⎜⎟ + Bmsh⎜⎟ + C⎢ ⎝ a ⎠ a ⎝ a ⎠⎢ mπy⎛ mπy⎞⎢+Dmch⎜⎟⎣ a ⎝ a ⎠m⎛ mπy⎞ ⎤sh⎜⎟ + ⎥⎝ a ⎠ ⎥⎥⎥⎦Página 14 <strong>de</strong> 17


Estructuras IIIDebido a que la <strong>de</strong>flexión es simétrica con respecto al eje x se tiene que C m =D m =0,entonces;ó,qoW =24Dqoa+D4 3 3( x − 2ax + a x)⎡⎢A⎣+⎛ mπy⎞ch⎜⎟ + B⎝ a ⎠mπy⎛ mπy⎞⎤⎛ mπx⎞sh⎜⎟⎥sin⎜⎟a ⎝ a ⎠⎦⎝ ⎠4∑ ∞mmm= 1,3,5aWq a 4 ∞ ⎡A ch m y B m y sh m y sin m xo4 ⎛ π ⎞ π ⎛ π ⎞⎤⎛ π ⎞= ∑ ⎢ 5 5+m⎜ ⎟ +m⎜ ⎟⎥⎜ ⎟D m=135 , , ⎣πm ⎝ a ⎠ a ⎝ a ⎠⎦⎝ a ⎠Sustituyendo las condiciones <strong>de</strong> bor<strong>de</strong>, W=0 y2∂ w2= 0 en y=±b/2, se obtiene,∂y4Am cham + Bm am sham+5 5= 0π mA + 2B cha + B a sha = 0( )m m m m m msiendo am=mπb.2 aDe estas ecuaciones se obtiene,Am( m)( )2( amtanh a + 2)=5 5π mcham2Bm=5 5π mcha( )mDe esta forma,4 ∞ ⎡4 qao 1 2W = ⎢5 ∑ 5 1 −π D m=135 , , m⎣⎢m( amtanh( am)+ ) ⎛ 2 amy⎞amch⎜⎟ +2 ch( a ) ⎝ b ⎠ 2 ch( a )m⎤2yb sh ⎛ 2a ysin m xm⎞⎜ ⎟⎥⎝ b ⎠⎦⎥ ⋅ ⎛ π ⎞⎜ ⎟⎝ a ⎠<strong>La</strong> <strong>de</strong>flexión máxima ocurre en el centro <strong>de</strong> la placa (x=a/2, y=0) y valeWmax4 ∞4 qao 1=5 ∑Dπ m=135 , ,( + 2)m−1( )⎡− 2 amtanh( am)⎢51−m ⎢ 2 ch( a )⎣m⎤⎥⎦⎥Página 15 <strong>de</strong> 17


Estructuras III<strong>La</strong> suma <strong>de</strong> los primeros términos <strong>de</strong> esta serie correspon<strong>de</strong> a la solución <strong>de</strong> lasección central <strong>de</strong> una tira cargada uniformemente, es <strong>de</strong>cir,Wmax( )m−1( ) 2 amtanh( am)5m 2 ch( a )4 4 ∞5 qao4 qao − 1 + 2= −3 ∑384 D Dπ m=135 , ,Esta serie converge rápidamente. Consi<strong>de</strong>rando a una placa cuadrada (a/b=1) setiene quea 1= π 3, a 23= π , etc..24 45 qao4 qaoWmax= −6( 0. 68562 − 0. 00025+......)384 D π D4qao≅ 0.00406DSe pue<strong>de</strong> observar que solo los dos términos <strong>de</strong> la serie se pue<strong>de</strong>n tomar en cuentapara hallar un resultado satisfactorio.Los momentos pue<strong>de</strong>n ser encontrados sustituyendo la expresión <strong>de</strong> W. Los valoresmáximos <strong>de</strong> estos momentos se hallan en x=a/2 y y=0.2qam−1⎛22µ= + − qa − 2 m⎜Am− B8m=135 , , ⎝ 1 − µ2∞qam−1⎛o2 222µ= µ + ( 1− µ ) qaoπ ∑ − 1 2 m⎜Am+ B8⎝ 1 − µ∞o2 2( M ) ( )x( 1 µ )oπ ∑ 1max( M ) ( )ymaxm=135 , ,mmm⎞⎟⎠⎞⎟⎠• Placa rectangular con dos lados simplemente apoyados, uno libre y elúltimo empotrada o simplemente apoyado.xbayFigura 8.x= 0 es un bor<strong>de</strong> simplemente apoyadox= a es un bor<strong>de</strong> simplemente apoyadoy= b es un bor<strong>de</strong> librey= 0 es un bor<strong>de</strong> empotrado.Página 16 <strong>de</strong> 17


Para este caso las condiciones <strong>de</strong> bor<strong>de</strong> son,Estructuras III2∂ wW=0 y2= 0 en x=0 y x=a∂x∂wW=0 y = 0 en y=0∂y 3∂µ ∂ 32 2⎛ ww ⎞⎛∂w ∂ w⎞⎜ 3+ 21 ( − )2 ⎟ = 0 ⎜ 2+2 ⎟ = 0 en y=b⎝ ∂y∂x∂y⎠⎝ ∂x∂y⎠Consi<strong>de</strong>rando una carga repartida q o , a través <strong>de</strong>l método <strong>de</strong> series simples <strong>de</strong> Levi,se tiene que,4 ∞4 qawD m sin m xo 1 π1=5 ∑ 5π m=135 , , ay∞⎛w Y y sin m π x ⎞2= ∑ m( ) ⎜ ⎟⎝ a ⎠m=135 , ,Es obvio que las dos primeras condiciones <strong>de</strong> bor<strong>de</strong> para x=0 y x=a son satisfechaspor w=w 1 +w 2 . Los coeficientes A m , B m , C m y D m <strong>de</strong> Y m <strong>de</strong>ben satisfacer las <strong>de</strong>máscondiciones <strong>de</strong> bor<strong>de</strong>.Tomando las condiciones para y=0 se tiene,A4=−5 5, C =−Dπ mm m mTomando las condiciones para y=b se tiene,Bm( + µ )( − µ ) ch² Bm + µ chBm − µ ( − µ ) BmshBm− ( − µ ² )( 3+ µ )( 1− µ ) ch²B + ( 1− µ ) B + ( 1+µ )4 3 1 2 1 1=5 5⋅2m2 2πmmCm( + µ )( − µ ) + µ ( + µ ) − µ ( − µ ) − ( − µ )2 2 2( 3+ µ )( 1− µ ) ch²B + ( 1− µ ) B + ( 1+µ )4 3 1 sh B ch B 1 shB 1 B chB 1 B=5 5⋅π mm m m m m mmm2Sustituyendo estas expresiones, se pue<strong>de</strong> hallar la <strong>de</strong>flexión <strong>de</strong> la placa, su máximovalor estará en el medio <strong>de</strong>l lado que se encuentra libre.Página 17 <strong>de</strong> 17

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!