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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

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!