Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 648 vistas

    Metodo Jacobi

    Cargado por

    Omar Caqui Zegarra
    El documento describe el método de Jacobi para calcular los modos de vibración de un sistema. El método diagonaliza la matriz dinámica mediante rotaciones sucesivas de la matriz, lo que permite obtener las frecuencias y períodos de vibración del sistema de manera independiente. Se presenta un ejemplo numérico completo del procedimiento aplicado a una matriz de 3x3.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Descargue como PDF, TXT o lea en línea desde Scribd

    Metodo Jacobi

    Cargado por

    Omar Caqui Zegarra
    648 vistas14 páginas
    El documento describe el método de Jacobi para calcular los modos de vibración de un sistema. El método diagonaliza la matriz dinámica mediante rotaciones sucesivas de la matriz, lo que permite obtener las frecuencias y períodos de vibración del sistema de manera independiente. Se presenta un ejemplo numérico completo del procedimiento aplicado a una matriz de 3x3.

    Derechos de autor

    © Attribution Non-Commercial (BY-NC)

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    El documento describe el método de Jacobi para calcular los modos de vibración de un sistema. El método diagonaliza la matriz dinámica mediante rotaciones sucesivas de la matriz, lo que permite obtener las frecuencias y períodos de vibración del sistema de manera independiente. Se presenta un ejemplo numérico completo del procedimiento aplicado a una matriz de 3x3.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    648 vistas14 páginas

    Metodo Jacobi

    Cargado por

    Omar Caqui Zegarra
    El documento describe el método de Jacobi para calcular los modos de vibración de un sistema. El método diagonaliza la matriz dinámica mediante rotaciones sucesivas de la matriz, lo que permite obtener las frecuencias y períodos de vibración del sistema de manera independiente. Se presenta un ejemplo numérico completo del procedimiento aplicado a una matriz de 3x3.

    Copyright:

    Attribution Non-Commercial (BY-NC)
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 14

    CURSO

    INGENIERIA SISMO
    RESISTENTE I
    Mtodos para calculo de modos de vibracin.Mtodos de iteracin matricial.- Mtodo de
    Jacobi

    Ing. Omart Tello Malpartida

    Mtodos de calculo de Modos de


    Vibracin
    METODOS DE CALCULO DE MODOS DE VIBRACION

    Metodos de Iteracin Matricial :


    Metodo de Jacobi
    Conocido tambien como metodo de diagonalizacion por rotaciones sucesivas, este
    procedimiento permite determinar, simultaneamente, todas las frecuencias y
    modos de sistemas complejos de hasta 200 grados de libertad.
    Consiste esencialmente, en diagonalizar la matriz dinamica o su inversa con el
    objeto de obtener, en la ecuacion matricial caracteristica, una serie de
    expresiones independientes, faciles de resolver.

    Ingeniera Sismo Resistente I

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi
    Datos:

    5 2.5 0
    K := 80 2.5 3.5 1
    0 1 1

    0
    0
    0.4077

    M := 1 0
    0.4077
    0
    0

    0
    0.2039

    Calculo de Matriz Diagonal Inferior [ L], tal que M=L.L

    M1 , 1

    L := 0

    0
    M2 , 2
    0

    M3 , 3
    0

    0
    0.639 0

    L = 0
    0.639
    0
    0

    0
    0.452

    Ingeniera Sismo Resistente I

    0
    0.639 0

    L= 0
    0.639
    0
    0

    0
    0.452

    0
    1.566 0

    1.566
    0
    = 0
    0

    0
    2.215

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi
    Calculo de Matriz Dinamica [ G*]
    1

    G := L

    ( )
    1

    K L

    981.114 490.557
    G = 490.557 686.779 277.467

    0
    277.467 392.349

    Ingeniera Sismo Resistente I

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi ( 1ra Rotacin)


    Ciclos de rotacin:
    Primera Rotacin
    Los ciclos de rotacin se inician eliminando , por ejemplo el termino g1,2 (r=1, s=2) ; as tenemos:
    Tag 2 = 2 g(r,s) /[g

    (r,r) -

    (s,s)]

    =0.5 arctag( 2 g(r,s) /[g (r,r) - g (s,s)])

    G 1 , 2 = 490.557
    G 1 , 1 = 981.114
    G 2 , 2 = 686.779

    2 G 1 , 2

    G1 , 1 G2 , 2

    1 := 0.5 atan
    1 = 0.64
    s := sin ( 1 )

    s = 0.597

    c := cos ( 1 )

    c = 0.802

    Ingeniera Sismo Resistente I

    c s 0
    s c 0

    0 0 1
    Ing. Omart Tello Malpartida

    Mtodo de Jacobi ( 1ra Rotacin)


    c s 0
    R1 := s c 0
    0 0 1

    0.802 0.597 0
    R1 = 0.597 0.802 0
    0

    0
    1

    Luego el triple producto escalar resulta ser : G1=R1.G.R1

    G1 := R1 GR1

    Ingeniera Sismo Resistente I

    165.628
    0
    1346.103

    321.79 222.61
    G1 = 0
    165.628 222.61 392.349

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi( 2da Rotacin)


    En forma similar para el termino g1,3 (r=1, s=3) ; se obtiene:
    G1 1 , 3 = 165.628
    G1 1 , 1 = 1.346 10

    G1 3 , 3 = 392.349
    2 G1 1 , 3

    G1 1 , 1 G1 3 , 3

    2 := 0.5 atan
    2 = 0.167

    cos ( 2 ) 0 sin ( 2 )
    R2 :=
    1
    0
    0
    sin ( 2 ) 0 cos ( 2 )

    0.986
    R2 =
    0
    0.166

    0.166

    0.986

    Luego el triple producto escalar resulta ser : G2 =R2.G1.R2


    T

    G2 := R2 G1 R2

    1374.047
    G2 = 37.034
    0

    Ingeniera Sismo Resistente I

    219.507
    321.79
    219.507 364.405
    37.034

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi( 3ra Rotacin)


    En forma similar para el termino g3,2 (r=3, s=2) ; se obtiene:
    G2 3 , 2 = 219.507
    G2 3 , 3 = 364.405
    G2 2 , 2 = 321.79
    2 G2 3 , 2

    G2

    G2
    2,2
    3,3

    3 := 0.5 atan
    3 = 0.737

    0
    0
    1
    R3 := 0 cos ( 3 ) sin ( 3 )
    0 sin ( 3 ) cos ( 3 )

    0
    0
    1
    R3 = 0 0.74 0.672
    0 0.672 0.74

    Luego el triple producto escalar resulta ser : G3 =R3.G2.R3


    T

    G3 := R3 G2R3

    1374.047 27.423 24.89


    G3 = 27.423 122.558
    0
    24.89

    0
    563.637

    Ingeniera Sismo Resistente I

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi( 4ta Rotacin)


    Terminado un ciclo se continua, reiniciando con la posicin de la primera rotacin:
    Para el termino g1,2 (r=1, s=2) ; se obtiene:
    G3 1 , 2 = 27.423
    G3 1 , 1 = 1.374 10

    G3 2 , 2 = 122.558
    2 G3 1 , 2

    G3

    G3
    2,2
    1,1

    4 := 0.5 atan
    4 = 0.022
    s := sin ( 4 )

    s = 0.022

    c := cos ( 4 )

    c = 1

    c s 0
    R4 := s c 0
    0 0 1

    1
    R4 = 0.022
    0

    0.022
    1
    0

    0
    1

    Luego el triple producto escalar resulta ser : G4 =R4.G3.R4


    T

    G4 := R4 G3 R4

    0
    24.884
    1374.648

    G4 =
    0
    121.958 0.545
    24.884

    0.545 563.637

    Ingeniera Sismo Resistente I

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi( 5ta Rotacin)


    En forma similar para el termino g1,3 (r=1, s=3) ; se obtiene:
    G4 1 , 3 = 24.884
    G4 1 , 1 = 1.375

    10

    G4 3 , 3 = 563.637
    2 G4 1 , 3

    G4 1 , 1 G4 3 , 3

    5 := 0.5 atan
    5 = 0.031

    cos ( 5 ) 0 sin ( 5 )
    1
    0
    0
    R5 :=
    sin ( 5 ) 0 cos ( 5 )

    1
    R5 =
    0
    0.031

    0.031

    Luego el triple producto escalar resulta ser : G5 =R5.G4.R5

    G5 := R5 G4 R5

    1375.41
    G5 = 0.017

    Ingeniera Sismo Resistente I

    0.017
    121.958
    0.545

    0.545

    562.874
    0

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi( 6ta Rotacin)


    En forma similar para el termino g3,2 (r=3, s=2) ; se obtiene:
    G5 3 , 2 = 0.545
    G5 3 , 3 = 562.874
    G5 2 , 2 = 121.958
    2 G5 3 , 2

    G5
    G5
    3,3
    2,2

    6 := 0.5 atan

    6 = 1.236 10

    0
    0
    1
    R6 := 0 cos ( 6 ) sin ( 6 )
    0 sin ( 6 ) cos ( 6 )

    0
    1

    R6 = 0
    1

    0 1.236 10

    1.236 10
    3

    Luego el triple producto escalar resulta ser : G5 =R5.G4.R5

    G6 := R6 G5R6

    0
    1375.41 0.017

    G6 = 0.017 121.957
    0
    0

    0
    562.875

    Ingeniera Sismo Resistente I

    Se puede observar que se


    obtenido la matriz diagonalizada,
    con la precisin esperada, por lo
    tanto se culmina el proceso de
    rotaciones

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi( Frecuencia y


    Periodos de Vibracin)

    FRECUENCIAS

    PERIODOS

    W1 :=

    G61 , 1

    W1 = 37.087

    T1 := 2

    W1

    T1 = 0.169

    3ro

    W2 :=

    G62 , 2

    W2 = 11.043

    T2 := 2

    W2

    T2 = 0.569

    1ro

    W3 :=

    G63 , 3

    W3 = 23.725

    T3 := 2

    W3

    T3 = 0.265

    2do

    Ingeniera Sismo Resistente I

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi
    ( Modos de Vibracin)

    U := R1R2R3R4R5R6

    0.77 0.37 0.52


    U = 0.62 0.65 0.45
    0.17 0.67 0.73

    Modos de Vibracion:

    ( 1)

    D := L

    Ingeniera Sismo Resistente I

    1.2 0.58 0.82


    D = 0.97 1.02 0.7
    0.39 1.47 1.61

    Ing. Omart Tello Malpartida

    Mtodo de Jacobi
    ( Modos de Vibracin)
    Modos de Vibracion Normalizados:
    i := 1 , 2 .. 3

    d1i := Di , 1

    1.202
    d1 = 0.966
    0.386

    d2i := Di , 2

    0.58
    d2 = 1.015
    1.473

    d3i := Di , 3

    0.82
    d3 = 0.699
    1.608

    Ingeniera Sismo Resistente I

    1 :=

    d1
    D1 , 1

    2 :=

    d2
    D1 , 2

    3 :=

    d3
    D1 , 3

    1.000
    1 = 0.804
    0.321

    3ro

    1.000
    2 = 1.751
    2.541

    1ro

    1.000
    3 = 0.853
    1.962

    2do

    Ing. Omart Tello Malpartida

    También podría gustarte