Analisis Estructural - Cálculo matricial de las deformaciones - Vigas continuas

download Analisis Estructural - Cálculo matricial de las deformaciones - Vigas continuas

of 12

Transcript of Analisis Estructural - Cálculo matricial de las deformaciones - Vigas continuas

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    1/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    ANALISIS ESTRUCTURAL

    CALCULO MATRICIAL DE LASDEFORMACIONES

    Rigideces de una barraVigas continuas

    Resistencia de Materials Prof: Mateu Moy Borrs

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    2/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    ANALISIS ESTRUCTURALCALCULO MATRICIAL DE LAS DEFORMACIONES

    Metodo de Rigidez

    Relaciona las acciones que inciden sobre una estructura con lasdeformaciones resultantes en la misma (giros y desplazamientos de los nudos)a travs de la Rigidez.

    Accin = Rigidez x Deformacin

    En el caso del metodo matricial restringido se considera que todas lasbarras son de longitud inalterable, es decir que no se acortan por el efecto delos esfuerzos normales. Este mtodo permite calcular la mayora de estructurasen edificacin sin que la simplificacin afecte de manera significativa a los

    resultados.

    Los nudos de la barra podrn desplazarse en el sentido de las flechas,es decir, perpendicularmente a la directriz inicial de la barra:

    Puede considerarse un desplazamiento perpendicular de los extremosya que el arco de circunferencia que describira es sumamente pequeo debidoa la pequeez de las deformaciones.

    El movivimento de la barra ir acompaado de unos giros en susextremos (dependiendo del grado de empotramiento de los mismos) y/o de undesplazamiento relativo de sus extremos.

    = desplazamiento relativo entre los nudos= giro que se produce en el nudo.

    Rigidez: Es la accin necesaria aplicar para ques e produzca una deformacinunitaria.

    Resistencia de Materials Prof: Mateu Moy Borrs

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    3/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    Si Def. = 1 Accin = Deformacin

    Criterio de signos en analisis extructural:

    Giros: antihorario (+)Horario (-)

    Movimientos: (+) (-)

    Clculo de las rigideces

    Supongamos una barra a la que se le produce un giro unidad (1 radin)

    en uno de sus extremos; cero en el otro y despalzamiento entre los nudostambin cero. Es decir, la barra puede girar en uno de sus extremos y tiene unempotramiento perfecto en el otro.

    1rad I = 1; D =0; = 0

    Para que dicha situacin de movimientos sea posible debern aparecer unosmomentos M y M.

    M

    M M+M

    Pero si solo apareciesen los momentos, la barra giraria debido al parM+M, por lo que para que la barra permanezca en posicin deben aparecerunas fuerzas que anulen dicho par.

    T

    MM

    T

    Las fuerzas que deben aparecer para que se produzca esta deformacinunitaria (giro 1 rad.) corresponden a las rigideces de la barra.

    A continuacin deduciremos el valor de M, M y T a partir de losteoremas de Mohr.

    Resistencia de Materials Prof: Mateu Moy Borrs

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    4/12

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    5/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    EI

    UV

    AB

    A

    BA == 0 (*3)0=AB

    AU

    Ahora podemos tomar dos de las tres condiciones halladas paraestablecer un sistema de dos ecuaciones con dos incgnitas y deducir el valorde M y M:

    EIAB = 0=AB

    AU

    EILMML

    =+

    2

    '

    2

    0

    3

    2*

    2

    '

    3

    *

    2

    =+ LLMLML

    L

    EIM

    2'= ,,

    L

    EIM

    4=

    Ahora ya conocemos los valores de M y M que producen la deformacin(giro) de la barra, pero nos falta determinar los valores de T para que losextremos se mantengan en su posicin inical:

    T LTL

    EI*

    6= ,,

    2

    6

    L

    EIT =

    L

    EI6

    T

    M, M y T seran por tanto las fuerzas necesarias aplicar a una barrapara producir una deformacin de un giro unitario (1 rad) en uno de susextremos y las denominaremos rigideces.

    L

    EI6

    L

    EI4

    L

    EI2

    L

    EI6

    Cada una de las rigideces relaciona la accin necesaria aplicar en esenudo (M o T) con la deformacin que se produce (giro () o desplazamiento()). Utilizaremos los subindices G y D segn el siguiente cuadro:

    Resistencia de Materials Prof: Mateu Moy Borrs

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    6/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    Giro (G) Desplazamiento (D)Accin M T

    Deformacin

    KGD

    KGG KGG

    KGD

    Como en la realidad el ngulo nunca ser 1 radian, el valor de lasfuerzas que apareceran en los nudos seran el giro que se produjesemultiplicado por la rigidez correspondiente:

    = i K* i KGG* iKGD* i(...)

    A continuacin analizaremos una barra en la que se produce undesplazamiento relativo unitario entre sus nudos extremos y suponiendo quesus extremos permanecen como empotramientos perfectos. Para ello sernecesario aplicar unos momentos (M , M) y unas fuerzas (T) compatibles conla deformacin propuesta:

    T1 M

    T M

    M (-)

    (+) M

    0==

    AB

    AB

    EI 0=AB M = M

    1=BA

    V 11

    *)3

    2*

    23

    1*

    2( =+

    EI

    LMLML 2

    6

    L

    EIM =

    Resistencia de Materials Prof: Mateu Moy Borrs

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    7/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    M+M T

    T2

    12

    L

    EI LT

    L

    EI*

    122= 3

    12

    L

    EIT=

    KGD KDD

    KDD KGD

    En el caso de que uno de los extremos de la barra fuese unaarticulacin, en ese nudo no aparecen momentos que se opongan a ladeformacin y las rigideces quedaran como sigue:

    KDDL

    EIKGG

    3=

    0

    1rad 0' =GGK

    KGG

    KDD 2

    3

    L

    EIKDD =

    1 KGD KDD 2

    3

    L

    EIKGD =

    KDD 0

    0)( =BGDK

    3

    3

    L

    EIKDD =

    En el caso general que nos encontremos con giro en los dos extremosde la barra y, adems, despalazmiento relativo entre los nudos

    T

    Mi Md T

    Resistencia de Materials Prof: Mateu Moy Borrs

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    8/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    Lo descompondremos en tres ecuaciones:

    - Primero el giro izquierdo = 1, el giro derecho = 0 y el desplazamiento = 0

    - Segundo: giro izquierdo = 0, giro derecho = 1 y desplazamiento = 0- Tercero: giro izquierdo =0, giro derecho = 0 y desplazamienmto = 1

    Con lo que se puede establecer la siguiente matriz:

    i=1 d=0 = 0

    i=0 d=1 = 0

    i=0 d=0 = 1

    Mi KGG KGG KGD iMd = KGG KGG KGD * dT KGD KGD KDD

    Accin = Rigidez * Deformacin

    De esta manera quedan relacionadas las acciones en los nudos de unabarra con las deformaciones en los mismos a travs de las rigideces de labarra.

    Ejemplo:

    Calcular la siguiente viga continua:

    2T 4T 5T 2T/m

    1m 2m 4m 2m 2m 6m

    A B C D

    El giro que se producir en los apoyos A, B y C estar condicionado porlas cargas a izquierda y derecha de los mismos y por las rigideces respectivasde las vigas que los acometen.

    El nudo A, al tener un voladizo en uno de sus lados, podemos considerarque gira libremente, lo cual simplifica el clculo sin que practicamente seaprecie en los resultados.

    En primer lugar calcularemos cada tramo de viga como si sus nudosfueran empotramientos, es decir, calcularemos sus momentos hiperestticos.

    Resistencia de Materials Prof: Mateu Moy Borrs

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    9/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    El voladizo del apoyo A lo sustituimos por el momento que producesobre el apoyo.

    El criterio de signos ser el expuesto al principio del tema.

    Tramo AB:

    4T 2T/m

    MBA mTQl

    l

    alPaM

    BA 06,115.182

    )(2

    2

    22

    =+

    =

    3mT

    Tramo BC:

    2T/m 5T

    MBC MCB mTPlQl

    MBC

    17.5812

    2

    =+= mTMCB 17.5=

    Tramo CD:

    2T/m mTQl

    MCD 00,6

    12

    2

    == mTMDC 00,6=

    M

    CDM

    DC

    Juntando los tres tramos nos quedara:

    -3 3 -11,06 5,17 -5,17 6,00 -6,00

    Estos resultados solo seran correctos en caso de que los apoyos B y Cfueran empotramientos perfectos.

    En nuestro caso, los momentos a izquierda y derecha de cada apoyodebern equilibrarse en funcin de las rigideces de las vigas que los acometen.Ello tambin afectar al empotramiento D.

    A continuacin vamos a calcular las rigideces de las barras en funcindel giro de cada nudo (solamente B y C ya que D es un empotramiento perfectoy no girar).

    Resistencia de Materials Prof: Mateu Moy Borrs

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    10/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    Si consideramos que la viga es de seccin y material constante,podemos sustituir el mdulo de elasticidad y el momento de inercia por 1 al sercomunes a todas las frmulas de las rigideces.

    Por otro lado, si imponemos la condicin de que en los apoyos de la viga

    no puede producirse ningun descenso, podremos prescindir de los valores KGDy KDD.

    En primer lugar impondremos un giro unitario al nudo B manteniendo elnudo C como empotramiento perfecto:

    0.50 0.50

    A B 1.00 C

    Tramo AB: M G

    G 50.06

    33===

    L

    EIKGG ,, 0' =GGK

    Tramo BC: M G

    G 00.14

    44===

    L

    EIK

    GG,, 50.0

    4

    22'

    ===L

    EIK

    GG

    A continuacin impondremos un giro unitario al nudo C manteniendo elnudo B como empotramiento perfecto:

    1.00 0.33

    A B 0.50 C 0.67

    Tramo BC: 00.14

    44===

    L

    EIK

    GG,, 50.0

    2' ==

    L

    EIK

    GG

    Tramo CD: 67.06

    44===

    L

    EIKGG ,, 33.0

    6

    22'

    ===L

    EIKGG

    Como ya se ha explicado anteriormente, los valores de las rigidecesseran los momentos necesarios a aplicar en los nudos para que estos girasenla unidad (1 radian).

    Resistencia de Materials Prof: Mateu Moy Borrs

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    11/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    Pero como en la realidad los nudos solo giran una pequea fraccin, losmomentos que aparezcan en los nudos por efecto de estos giros sern el valorde las rigideces por el giro real de cada nudo.

    El momento final en cada nudo ser la suma de estos momentos mas

    los momentos hiperestticos hallados inicialmente.

    Para ello plantearemos el equilibrio esttico de cada nudo, es decir, quela suma de todos los momentos en un nudo sea igual al momento puro quehaya en l.

    En nuestro caso: Momento puro en B=0Momento puro en C=0

    Por tanto las ecuaciones de equilibrio en B y C sern:

    0 = -11,06+5,17+(1+0,5)B+0,5c =-5,93+1,5B+0,5c B= 4.5460 = -5,17+6,00+0,5B+(1+0,67)c =0,83+0,5B+1,67c C= -1.858

    En consecuencia, los momentos reales que se producirn en cadaextremo de barra debido al giro de los nudos ser:

    mTKM BAB

    GG

    AB

    BA 273,250,0*546,4* ===

    mTKKM CBC

    GGB

    BC

    GG

    BC

    BC 617,350.0*858,100,1*546,4** ' ==+=

    mTKKM BBC

    GGC

    BC

    GG

    BC

    CB 415.050,0*546,400,1*858,1** ' =+=+=

    mTKM CCD

    GG

    CD

    CD 239,167,0*858,1* ===

    mTKM CCD

    GG

    CD

    DC 619,033,0*858,1*' ===

    Sumando estos momentos a los momentos hiperestticos de cada barra,tentremos el equilibrio de momentos en los nudos (apoyos) de la viga:

    -3,00 3,00 -11,06 5,15 -5,17 6,00 -6,002,273 3,617 0,415 -1,239 -0,619

    --------- --------- --------- --------- -----------8,79 8,79 -4,76 4,76 -6,619

    Una vez hallado el diagrama de momentos hiperestticos de la viga,representado por (---) en el dibujo siguiente ser necesario colgar del mismo losdiagramas isostticos de cada una de las vigas, teniendo en cuenta que el calculoisosttico del tramo AB se hara incluyendo el voladizo:

    Resistencia de Materials Prof: Mateu Moy Borrs

  • 8/6/2019 Analisis Estructural - Clculo matricial de las deformaciones - Vigas continuas

    12/12

    Universitat de les Illes Balears Arquitectura Tcnica Departament de Fsica

    A B C D

    A continuacin plantearemos el mismo problema en forma matricial:

    Las ecuaciones de equilibrio en los nudos de la estructura quedara de lasiguiente forma:

    0 -5,89 1,50 0,50 B= + *

    0 0,83 0,50 1,67 C

    Para obtener los momentos finales en los extremos de barrasplanteariamos:

    MAB -11,06 0,50 0MBC 5,17 1,00 0,50 BMCB = -5,17 + 0,50 1,00 *MCD 6,00 0 0,67 CMDC -6,00 0 0,33

    Resistencia de Materials Prof: Mateu Moy Borrs