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

Matriz de Amortiguamiento

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 2

DETERMINACIÓN DE LA MATRIZ DE AMORTIGUAMIENTO POR RAYLEIGH

1) Ensamblamos la matriz de rigidez "K" (Matriz de rigidez lateral condensada)

10235.21 -5661.22 686.16 -16.56


-5661.22 9934.87 -5605.13 645.19
KL = 686.16 -5605.13 9617.69 -4671.76 Tn/m
-16.56 645.19 -4671.76 4031.25

La matriz de masas ha quedado asi:

3.270 0.000 0.000 0.000


0.000 3.234 0.000 0.000
M= 0.000 0.000 3.234 0.000 Tn.seg²/m
0.000 0.000 0.000 3.050

La matriz modal ha quedado así: Las frecuencia naturales para cada modo son:

1.0000 1.0000 1.0000 1.0000 ω1 = 11.42 rad/seg.


2.0736 1.1219 -0.2328 -1.4100 ω2 = 34.75 rad/seg.
Φ= 2.8962 0.0684 -0.8985 1.1704 ω3 = 57.80 rad/seg.
3.3602 -1.1152 0.6545 -0.4662 ω4 = 76.28 rad/seg.

2) Ahora determinamos a0 y a1 a partir de la ecuación…………………… 1 1/ω𝑖 ω𝑖 𝑎0 ξ𝑖


{ } = { } ........ (1)
Para ello especificaremos para el primer y tercer modo ξ = 5% 2 1/ω𝑗 ω𝑗 𝑎1 ξ𝑗
Entonces reemplazando las matrices en la ecuación (1), tenemos: Ver ecuación 11.4.9 de Chopra

1
11.42 a0 0.05
11.42
= 2.00 *
1
57.80 a1 0.05
57.80

0.0876 𝑎0 + 11.42 𝑎1 = 0.10


0.0173𝑎0 + 57.80 𝑎1 = 0.10

Solucionando tenemos:
0.10 11.42 0.0876 0.10 4.8649 a0 = 0.95342
0.10 57.80 0.0173 0.10 4.6383 a1 = 0.00144
0.0070
3) Ahora calculamos la matriz de amortiguamiento:

17.9051 -8.1788 0.9913 -0.0239


-8.1788 17.4362 -8.0978 0.9321
C = a0m + a1K = 0.9913 -8.0978 16.9780 -6.7494 Tn.seg./m
-0.0239 0.9321 -6.7494 8.7322

𝑎0 1 𝑎1
4) Ahora se procederá a calcular ξ2 y ξ4, a partir de la ecuación: ξ𝑛 = + ω , resumiendo:
2 ω𝑛 2 𝑛
ξ1 = 0.050 Ver ecuación 11.4.8 de Chopra
ξ2 = 0.039
ξ3 = 0.050
ξ4 = 0.061

UNIVERSIDAD NACIONAL DEL SANTA//INGENIERÍA ANTISÍSMICA 2020-1 Ing. IVÁN LEÓN MALO CIP 93535
DETERMINACIÓN DE LA MATRIZ DE AMORTIGUAMIENTO POR CAUGHEY

1) La serie de Caughey para un sistema de 4 grados de libertad es:

C = 𝑎0 𝒎 + 𝑎1 𝒌 + 𝑎2 𝒌𝒎−1 𝒌 + 𝑎3 𝒎(𝒎−1 𝒌)3 Ver ecuación 11.4.11 de Chopra

2) Determine a0, a1, a2, a3, a partir de la ecuación siguiente:

𝐽−1
1
Ver ecuación 11.4.14 de Chopra ξ𝑛 = ෍ 𝑎𝑙 ω𝑛2𝑙−1 Dónde: J, es número grados de libertad
2
𝑙=0

𝑎0 1 𝑎1 𝑎2 3 𝑎3 5
ξ𝑛 = + ω + ω 𝑛+ ω 𝑛
2 ω𝑛 2 𝑛 2 2

Ó matricialmente:

0.0876 11.42 1.49E+03 1.94E+05 a0 0.05


0.0288 34.75 4.20E+04 5.07E+07 a1 0.05
0.0173 57.80 1.93E+05 6.45E+08 x a2 = 2* 0.05
0.0131 76.28 4.44E+05 2.58E+09 a3 0.05

Solucionando tenemos:
Det = 274070994714354.00

0.10 11.42 1.49E+03 1.94E+05 0.0876 0.10 1.49E+03 1.94E+05


0.10 34.75 4.20E+04 5.07E+07 0.0288 0.10 4.20E+04 5.07E+07
0.10 57.80 1.93E+05 6.45E+08 0.0173 0.10 1.93E+05 6.45E+08
0.10 76.28 4.44E+05 2.58E+09 0.0131 0.10 4.44E+05 2.58E+09

0.0876 11.42 0.10 1.94E+05 0.0876 11.42 1.49E+03 0.10


0.0288 34.75 0.10 5.07E+07 0.0288 34.75 4.20E+04 0.10
0.0173 57.80 0.10 6.45E+08 0.0173 57.80 1.93E+05 0.10
0.0131 76.28 0.10 2.58E+09 0.0131 76.28 4.44E+05 0.10

a0 = 0.78162
a1 = 0.00283
a2 = -5.56514E-07
a3 = 4.66801E-11

3) Calculamos la matriz de amortiguamiento reemplazando los valores de a0, a1, a2, a3, en la serie:

C = 𝑎0 𝒎 + 𝑎1 𝒌 + 𝑎2 𝒌𝒎−1 𝒌 + 𝑎3 𝒎(𝒎−1 𝒌)3

17.3986 -5.6067 -0.4072 -0.2550


-5.6067 15.9195 -6.0215 -0.6870
C = -0.4072 -6.0215 15.3040 -6.1716 Tn.seg./m
-0.2550 -0.6870 -6.1716 9.3022

UNIVERSIDAD NACIONAL DEL SANTA//INGENIERÍA ANTISÍSMICA 2020-1 Ing. IVÁN LEÓN MALO CIP 93535

También podría gustarte