TRABAJO PRÁCTICO DE INGENIERIA SISMORESISTENTE

𝑬𝒋𝒆𝒓𝒄𝒊𝒄𝒊𝒐 𝑵𝒓𝒐 𝟎𝟏

Se tiene una caja cuya masas P= 2 últimos dígitos de su dni en kg que esta soldado desde
una atura de 1 metro al centro de la viga de masa despreciable. La sección de la viga es
de 30 x 50 cm. El módulo de elasticidad E= es 27000 Mpa, L = # dígitos de su primer
apellido en metros. Calcular las propiedades dinámicas, la definición del movimiento
oscilatorio que genera la máxima deflexión vertical y las fuerzas máximas.

Hallamos los valores de la masa de la caja y la longitud de la viga.

Masa (m)
Julio = 02
Berly = 59
Kelvin = 03
Erickson = 37
Masa = 25.25 Kg
Longitud (L)
Julio Larico =6
Berly Condori =7
Kelvin Cuayla =6
Erickson Chicche =7
L = 6.5 m

pág. 1
Datos:
M = 25.25 Kg
L = 6.5 m
E = 27 000 Mpa
Sección viga = 0.3 x 0.5 m
H=1m
SOLUCION:
Suponiendo que cuando se suelta la caja y en el instante que hace contacto con la viga, ambas
quedan adheridas y empieza a vibrar la viga.

Hallamos la Inercia de la viga

𝑏ℎ3 0.3(0.5)3
𝐼= =
12 12
𝑰 = 𝟎. 𝟎𝟎𝟑𝟏𝟐𝟓 𝒎𝟒

Hallamos la rigidez
𝑃𝐿3
𝛿=
48𝐸𝐼
𝑃 48𝐸𝐼 48(27000)(0.003125)
𝐾= = 3 =
𝛿 𝐿 6.53
𝑲 = 𝟏𝟒. 𝟕𝟓 × 𝟏𝟎𝟔 𝑵/𝒎
Hallamos la frecuencia natural

𝐾 14.75 × 106 𝑁/𝑚

𝜔𝑛 = √ = √
𝑚 25.25 Kg

𝝎𝒏 = 𝟕𝟔𝟒. 𝟑 𝒓𝒂𝒅/𝒔

pág. 2
Hallamos el periodo
2𝜋 2𝜋
𝑇= =
𝜔𝑛 764.3

𝑻 = 𝟖. 𝟐𝟐 × 𝟏𝟎−𝟑 𝒔𝒆𝒈

Hallamos la frecuencia

𝑓 = 1/𝑇 = 1/8.22 × 10−3

𝒇 = 𝟏𝟐𝟏. 𝟔𝟒 𝒉𝒛

Sin amortiguamiento tenemos:

𝒖̇
𝒖(𝒕) = 𝜔𝟎 × 𝐬𝐢𝐧(𝜔𝑛 × 𝑡) + 𝒖𝟎 × 𝐜𝐨𝐬(𝜔𝑛 × 𝑡) ………………… (1)
𝑛

En el instante en t=0, entonces 𝑢0 = 0

Por Conservación de Energía hallaremos la velocidad de la caja junto con la viga después de
haber caído desde 1m.
𝐸2 = 𝐸1
𝐸𝑐𝑖𝑛é𝑡𝑖𝑐𝑎 = 𝐸𝑃𝑜𝑡𝑒𝑛𝑐𝑖𝑎𝑙
1
× 𝑚𝑣 2 = 𝑚𝑔ℎ
2
𝑣 2 = 2𝑔ℎ

𝑣 = √2𝑔ℎ

𝑣 = √2(9.81𝑚/𝑠 2 )(1𝑚)

𝒗 = 𝒖𝟎̇ = 𝟒. 𝟒𝟑 𝒎/𝒔
Reemplazamos en la ecuación (1) tenemos la ecuación de posición:
𝒖𝟎̇
𝒖(𝒕) = × 𝐬𝐢𝐧(𝜔𝑛 × 𝑡) + 𝒖𝟎 × 𝐜𝐨𝐬(𝜔𝑛 × 𝑡)
𝜔𝑛
4.43
𝒖(𝒕) = × sin(764.3 × 𝑡) + (0) × cos(764.3 × 𝑡)
764.3
𝒖(𝒕) = 𝟎. 𝟎𝟎𝟓𝟕𝟗𝟔 × 𝐬𝐢𝐧(𝟕𝟔𝟒. 𝟑 × 𝒕)

pág. 3
Derivamos la ecuación de posición para obtener la ecuación de la velocidad:
𝑑(𝑢(𝑡))
= 𝑢̇ (𝑡) = 0.005796 × 764.3 × cos(764.3 × 𝑡)
𝑑𝑡
𝒖̇ (𝒕) = 𝟒. 𝟒𝟑 × 𝐜𝐨𝐬(𝟕𝟔𝟒. 𝟑 × 𝒕)
Derivamos la ecuación de velocidad para obtener la ecuación de la aceleración:
𝑑(𝑢̇ (𝑡))
= 𝑢̈ (𝑡) = −4.43 × 764.3 × sin(764.3 × 𝑡)
𝑑𝑡
𝒖̈ (𝒕) = −𝟑𝟑𝟖𝟓. 𝟖𝟒𝟗 × 𝐬𝐢𝐧(𝟕𝟔𝟒. 𝟑 × 𝒕)

pág. 4
𝑬𝒋𝒆𝒓𝒄𝒊𝒄𝒊𝒐 𝑵𝒓𝒐 𝟎𝟐

Un aparato de aire acondicionado cuyo peso P= 2 ultimos dígitos de su codigo en lb se

atornilla en medio de dos vigas paralelas de acero simplemente apoyadas (figura 2). La
luz libre de las vigas es L = # dígitos de su primer nombre en pies. El momento de inercia
de la sección transversal de cada viga es de 10 pulg4. El motor de la unidad funciona a
200 rpm y, a esta velocidad, produce una fuerza vertical desbalanceada de 60 lb.
Desprecie el peso de las vigas y suponga 1% de amortiguamiento viscoso en el sistema.
Para el acero E = 30,000 ksi. Considere la fuerza desbalanceada y determine las
amplitudes de la respuesta de régimen de desplazamiento y aceleración en medio de la
viga

Figura 2. Aparato de aire acondicionado

Hallamos los valores del peso y la longitud de las vigas.
Peso (W)
Julio = 36
Berly = 45
Kelvin = 26
Erickson = 01
El promedio es entonces, W = 27 lb
Longitud (L)
Julio =5
Berly =5
Kelvin =6
Erickson =8
El promedio es entonces, L = 6 pies

pág. 5
Datos:
W = 27 lb
L = 6 pies
E = 30 000 ksi
I = 10 pulg4
ξ = 0.01
𝑃0 = 60 𝑙𝑏
200
𝜔= (2𝜋) = 20.94 𝑟𝑎𝑑/𝑠
60

SOLUCION:
Hallamos la rigidez para las dos vigas
48𝐸𝐼 48(30000)(10)
𝐾 = 2( 3
) = 2( )
𝐿 (6 × 12)3
𝑲 = 𝟕𝟕. 𝟏𝟔 𝒍𝒃/𝒑𝒖𝒍𝒈
Hallamos la Frecuencia natural

𝐾 𝐾 77.16
𝜔𝑛 = √ = √ =√
𝑚 𝑊/𝑔 27/387

𝝎𝒏 = 𝟑𝟑. 𝟐𝟔 𝒓𝒂𝒅/𝒔
Hallamos el periodo
2𝜋 2𝜋
𝑇= =
𝜔𝑛 33.26
𝑻 = 𝟎. 𝟏𝟗 𝒔𝒆𝒈
Hallar el Factor de modificación dinámica
𝜔 20.94 𝑟𝑎𝑑/𝑠
𝐵= =
𝜔𝑛 33.26 𝑟𝑎𝑑/𝑠
𝑩 = 𝟎. 𝟔𝟑

pág. 6
Reemplazamos en:
1
𝑅𝑑 =
√(1 − 𝐵2 )2 + (2𝜉𝐵)2
1
𝑅𝑑 =
√(1 − 0.632 )2 + (2 × 0.01 × 0.63)2

𝑹𝒅 = 𝟏. 𝟔𝟓𝟖
Entonces

𝑢(𝑡) = 𝑢𝑒𝑠𝑡 . 𝑅𝑑. 𝑒 𝑖(𝜔𝑡−𝜃)

𝑃0
𝑢𝑒𝑠𝑡 =
𝐾
𝑃0
𝑢 (𝑡 ) = . 𝑅𝑑. 𝑒 𝑖(𝜔𝑡−𝜃)
𝐾

La amplitud máxima será cuando

𝑖 (𝜔𝑡 − 𝜃) = 0

Finalmente tenemos:
𝑃0
𝑢 (𝑡 ) = . 𝑅𝑑
𝐾
60 𝑙𝑏
𝑢 (𝑡 ) = × 1.658
77.16 𝑙𝑏/𝑝𝑢𝑙𝑔

𝒖(𝒕) = 𝟏. 𝟐𝟗 𝒑𝒖𝒍𝒈
Para la Aceleración tenemos
𝑃0
𝑢̈ (𝑡 ) = . 𝑅𝑑. 𝜔2 . 𝑒 𝑖(𝜔𝑡−𝜃)
𝐾

pág. 7
Aplicando la misma condición tenemos
𝑃0
𝑢̈ (𝑡 ) = . 𝑅𝑑. 𝜔2
𝐾
60
𝑢̈ (𝑡 ) = × 1.658 × (20.94)2
77.16
𝒖̈ (𝒕) = 𝟓𝟔𝟓. 𝟑𝟐 𝒑𝒖𝒍𝒈/𝒔𝟐

𝐄𝐣𝐞𝐫𝐜𝐢𝐜𝐢𝐨 𝐍𝐫𝐨 𝟎𝟑

Determinar la respuesta del tanque elevado, mostrado en la figura 3, sometida a la excitación

con un periodo de 0,64 segundos como se muestra en la figura 4. Considerar una razón de
amortiguamiento del 5%. P=tres últimos dígitos de su código multiplicado por 10.

Figura 3 y 4. Deposito elevado propuesto y función de excitación

Definición de datos:

P = 10(X) Kp

pág. 8
 𝐂ó𝐝𝐢𝐠𝐨𝐬 𝐝𝐞 𝐥𝐚𝐮𝐦𝐧𝐨𝐬:

• Berly Antony Condori : 145 (se considerará el mayor valor)

• Julio larico Ramos : 016
• Kelvin Cuayla Flores : 026
• Erickson Chicche incahcuaman : 001

𝐏 = 10(145) = 𝟏𝟒𝟓𝟎 𝐊𝐩
Solución:

W 1450kp 𝐬𝟐
𝐚) 𝐌𝐚𝐬𝐚 𝐝𝐞𝐥 𝐬𝐢𝐬𝐭𝐞𝐦𝐚: m = = cm = 𝟏. 𝟒𝟖 𝐤𝐩 ∗ 𝐜𝐦
g 980 s

kp 𝐤𝐩
𝐛) 𝐂𝐨𝐞𝐟𝐢𝐜𝐢𝐞𝐧𝐭𝐞 𝐝𝐞 𝐫𝐢𝐠𝐢𝐝𝐞𝐳 = 1000 = 𝟏𝟎
m 𝐜𝐦

𝐬
𝐜) 𝐂𝐨𝐞𝐟𝐢𝐜𝐢𝐞𝐧𝐭𝐞 𝐝𝐞 𝐚𝐦𝐨𝐫𝐭𝐢𝐠𝐮𝐚𝐜𝐢𝐨𝐧: C = 2ξ√km = 2 ∗ 0.1 ∗ √10 ∗ 1.48 = 𝟎. 𝟕𝟕 𝐤𝐩 ∗
𝐜𝐦

Ecuacion de la recta para A,B,C:

(0.16,1450)
1450

(0.64,0)
(0,0)
(0.32,0)

-1450
(0.48,-1450)

pág. 9
 PUNTOS: (t1, P1) ; (t2, P2)

P2 − P1
P(t) − P1 = (t − t1)
t2 − t1

▪ RECTA A:

𝐏𝐔𝐍𝐓𝐎𝐒: (𝟎, 𝟎) ; (𝟎. 𝟏𝟔, 𝟏𝟒𝟓𝟎)

1450 − 0
P(t) − 0 = (t − 0)
0.16 − 0

𝐏(𝐭) = 𝟗𝟎𝟔𝟐. 𝟓𝐭
▪ RECTA B:

𝐏𝐔𝐍𝐓𝐎𝐒: (𝟎. 𝟑𝟐, 𝟎) ; (𝟎. 𝟏𝟔, 𝟏𝟒𝟓𝟎)

1450 − 0
P(t) − 0 = (t − 0.32)
0.16 − 0.32

𝐏(𝐭) = −𝟗𝟎𝟔𝟐. 𝟓𝐭 + 𝟐𝟗𝟎𝟎

▪ RECTA C:

𝐏𝐔𝐍𝐓𝐎𝐒: (𝟎. 𝟔𝟒, 𝟎) ; (𝟎. 𝟒𝟖, −𝟏𝟒𝟓𝟎)

−1450 − 0
P(t) − 0 = (t − 0.64)
0.48 − 0.64

𝐏(𝐭) = 𝟗𝟎𝟔𝟐. 𝟓𝐭 − 𝟓𝟖𝟎𝟎

pág. 10
Para la carga periódica por series de Fourier:

Se representa por:

𝛂 𝛂

𝐏(𝐭) = 𝐚𝐨 + ∑ 𝐚𝐣 𝐜𝐨𝐬𝐰𝐣𝐭 + ∑ 𝐛𝐣 𝐬𝐞𝐧𝐰𝐣𝐭

𝐣=𝟏 𝐣=𝟏

Calculo de :

𝟏 𝐓𝐨
𝐚𝐨 = ∫ 𝐏(𝐭)𝐝𝐭
𝐓𝐨 𝟎
𝟐 𝐓𝐨
𝐚𝐣 = ∫ 𝐏(𝐭)𝐜𝐨𝐬𝐰𝐣𝐭𝐝𝐭
𝐓𝐨 𝟎

𝟐 𝐓𝐨
𝐛𝐣 = ∫ 𝐏(𝐭)𝐬𝐞𝐧𝐰𝐣𝐭𝐝𝐭
𝐓𝐨 𝟎

Entonces:
𝟏 𝐓𝐨
𝐚𝐨 = ∫ 𝐏(𝐭)𝐝𝐭
𝐓𝐨 𝟎
1 To
ao = ∫ P(t)dt =
To 0

+ =0

𝟐 𝐓𝐨
𝐚𝐣 = ∫ 𝐏(𝐭)𝐜𝐨𝐬𝐰𝐣𝐭𝐝𝐭
𝐓𝐨 𝟎

0.16 0.48
2
𝐚𝐣 = [∫ 9062.5tcoswjtdt + ∫ (−9062.5t + 2900)coswjtdt
To 0 0.16
0.64
+ ∫ (9062.5t − 5800)coswjtdt]
0.48

pág. 11
 Tramo recta A:

0.16 0.16
∫ 9062.5tcoswjtdt = 9062.5 ∫ tcoswjtdt
0 0

Integrando por partes:

uv − ∫ vdu

u=t
du = dt

dv = cos wj t dt
1
v = sen wj t
wj

0.16
t 1
9062.5( sen wj t − ∫ sen wj t dt)
wj 0 wj

t 1 0.16
9062.5 ( sen wj t + 2 cos wj t ) |
wj wj 0

𝟎. 𝟏𝟔 𝟏 𝟏
𝟗𝟎𝟔𝟐. 𝟓 ( 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟏𝟔 + 𝟐 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟏𝟔 − )
𝐰𝐣 𝐰𝐣 𝐰𝐣𝟐

Tramo recta B:

0.48 0.48 0.48

∫ (−9062.5t + 2900)coswjtdt = −9062.5 ∫ tcoswjtdt + 2900 ∫ coswjtdt
0.16 0.16 0.16

A B

pág. 12
 0.48
𝐀 = −9062.5 ∫ tcoswjtdt
0.16

Integrando por partes: ya calculado del tramo de la recta A; tenemos:

t 1 0.48
−9062.5 ( sen wj t + 2 cos wj t ) |
wj wj 0.16

𝟎. 𝟒𝟖 𝟎. 𝟏𝟔 𝟏 𝟏
−𝟗𝟎𝟔𝟐. 𝟓 ( 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟒𝟖 − 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟏𝟔 + 𝟐 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟒𝟖 − 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟏𝟔 )
𝐰𝐣 𝐰𝐣 𝐰𝐣 𝐰𝐣𝟐

0.48
𝐁 = 2900 ∫ coswjtdt
0.16

1 0.48
2900 sen wj t |
wj 0.16

𝟐𝟗𝟎𝟎
(𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟒𝟖 − 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟏𝟔)
𝐰𝐣

Tramo recta C:

0.64 0.64 0.64

∫ (9062.5t − 5800)coswjtdt = 9062.5 ∫ tcoswjtdt − 5800 ∫ coswjtdt
0.48 0.48 0.48

A B

0.64
𝐀 = 9062.5 ∫ tcoswjtdt
0.48

Integrando por partes: ya calculado del tramo de la recta A; tenemos:

pág. 13
 t 1 0.64
9062.5 ( sen wj t + 2 cos wj t ) |
wj wj 0.48

𝟎. 𝟔𝟒 𝟎. 𝟒𝟖 𝟏 𝟏
𝟗𝟎𝟔𝟐. 𝟓 ( 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟔𝟒 − 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟒𝟖 + 𝟐 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟔𝟒 − 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟒𝟖 )
𝐰𝐣 𝐰𝐣 𝐰𝐣 𝐰𝐣𝟐

0.64
𝐁 = −5800 ∫ coswjtdt
0.48

1 0.64
−5800 sen wj t |
wj 0.48

𝟓𝟖𝟎𝟎
− (𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟔𝟒 − 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟒𝟖)
𝐰𝐣

Considerando:

✓ 0.64=To

2 To To 1 To 1
[9062.5 ( sen wj + 2 cos wj − 2)
To 4wj 4 wj 4 wj
3To 3To To To 1 3To
− 9062.5 ( sen wj − sen wj + 2 cos wj
4wj 4 4wj 4 wj 4
1 To 2900 3To To
− 2
cos wj )+ (sen wj − sen wj )
wj 4 wj 4 4
To 3To 3To 1
+ 9062.5 ( senwjTo − sen wj + 2 cos wj To
wj 4wj 4 wj
1 3To 5800To 3To
− 2 cos wj )− (sen wj To − sen wj )]
wj 4 j2π 4

Quedando:

pág. 14
 2 9062.5 To 9062.5 9062.5 3To 9062.5 To
[+ 2
cos wj − 2
− 2
cos wj + 2
cos wj
To wj 4 wj wj 4 wj 4
9062.5 9062.5 3To
+ + 2
cos wj To − 2
cos wj ]
wj wj 4

2 18125 To 9062.5 18125 3To 9062.5

[+ 2
cos wj − 2
− 2
cos wj + + cos wj To]
To wj 4 wj wj 4 wj2

j2π
Considerando : wj =
To

Factorizando: y simplificando

2 18125To2 j2π To 9062.5To2 18125To2 j2π 3To

[+ 2
cos − 2
− 2
cos
To (j2π) To 4 (j2π) (j2π) To 4
2
9062.5To j2π
+ cos To]
(j2π)2 To

18125To jπ 9062.5To 18125To j3π 9062.5To j2π

cos − − cos + cos To
2(jπ)2 2 2(jπ)2 2(jπ)2 2 2(jπ)2 To

18125To jπ 18125To j3π

2
cos − 2
cos
2(jπ) 2 2(jπ) 2

Si:

j es par:

• Para 2,6,10,14,…etc=0
• Para 4,8,12,16,…etc=0

j es impar = 0

𝐚𝐣 = 𝟎

pág. 15
 𝟐 𝐓𝐨
𝐛𝐣 = ∫ 𝐏(𝐭)𝐬𝐞𝐧𝐰𝐣𝐭𝐝𝐭
𝐓𝐨 𝟎

0.16 0.48
2
𝐛𝐣 = [∫ 9062.5tsenwjtdt + ∫ (−9062.5t + 2900)senwjtdt
To 0 0.16
0.64
+ ∫ (9062.5t − 5800)senwjtdt]
0.48

Tramo recta A:

0.16 0.16
∫ 9062.5tsenwjtdt = 9062.5 ∫ tsenwjtdt
0 0

Integrando por partes:

uv − ∫ vdu

u=t
du = dt

dv = sen wj t dt
1
v = − cos wj t
wj

0.16
t 1
9062.5(− cos wj t − ∫ − cos wj t dt)
wj 0 wj

t 1 0.16
9062.5 (− cos wj t + 2 sen wj t ) |
wj wj 0

𝟎. 𝟏𝟔 𝟏
𝟗𝟎𝟔𝟐. 𝟓 (− 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟏𝟔 + 𝟐 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟏𝟔 )
𝐰𝐣 𝐰𝐣

pág. 16
 Tramo recta B:

0.48 0.48 0.48

∫ (−9062.5t + 2900)senwjtdt = −9062.5 ∫ tsenwjtdt + 2900 ∫ senwjtdt
0.16 0.16 0.16

A B

0.48
𝐀 = −9062.5 ∫ tsenwjtdt
0.16

Integrando por partes: ya calculado del tramo de la recta A; tenemos:

t 1 0.48
−9062.5 (− cos wj t + 2 sen wj t ) |
wj wj 0.16

𝟎. 𝟒𝟖 𝟎. 𝟏𝟔 𝟏 𝟏
−𝟗𝟎𝟔𝟐. 𝟓 (− 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟒𝟖 + 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟏𝟔 + 𝟐 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟒𝟖 − 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟏𝟔 )
𝐰𝐣 𝐰𝐣 𝐰𝐣 𝐰𝐣𝟐

0.48
𝐁 = 2900 ∫ senwjtdt
0.16

1 0.48
−2900 cos wj t |
wj 0.16

−𝟐𝟗𝟎𝟎
(𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟒𝟖 − 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟏𝟔)
𝐰𝐣

pág. 17
 Tramo recta C:

0.64 0.64 0.64

∫ (9062.5t − 5800)senwjtdt = 9062.5 ∫ tsenwjtdt − 5800 ∫ senwjtdt
0.48 0.48 0.48

A B

0.64
𝐀 = 9062.5 ∫ tsenwjtdt
0.48

Integrando por partes: ya calculado del tramo de la recta A; tenemos:

t 1 0.64
−9062.5 (− cos wj t + 2 sen wj t ) |
wj wj 0.48

𝟎. 𝟔𝟒 𝟎. 𝟒𝟖 𝟏 𝟏
𝟗𝟎𝟔𝟐. 𝟓 (− 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟔𝟒 + 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟒𝟖 + 𝟐 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟔𝟒 − 𝐬𝐞𝐧 𝐰𝐣 𝟎. 𝟒𝟖 )
𝐰𝐣 𝐰𝐣 𝐰𝐣 𝐰𝐣𝟐

0.64
𝐁 = −5800 ∫ senwjtdt
0.48

1 0.64
5800 cos wj t |
wj 0.48

𝟓𝟖𝟎𝟎
(𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟔𝟒 − 𝐜𝐨𝐬 𝐰𝐣 𝟎. 𝟒𝟖)
𝐰𝐣

Considerando:

✓ 0.16=To/4
✓ 0.32=To/2
✓ 0.48=3To/4

pág. 18
 ✓ 0.64=To

En resumen tenemos:

2 To To 1 To
[9062.5 (− cos wj + 2 sen wj )
To 4wj 4 wj 4
3To 3To To To 1 3To
− 9062.5 (− cos wj + cos wj + 2 sen wj
4wj 4 4wj 4 wj 4
1 To 2900 3To To
− 2 sen wj )− (cos wj − cos wj )
wj 4 wj 4 4
To 3To 3To 1
+ 9062.5 (− cos wj To + cos wj + 2 sen wj To
wj 4wj 4 wj
1 3To 5800 3To
− 2 sen wj )+ (cos wj To − cos wj )]
wj 4 wj 4

Considerando:

2π
wj = j ∗
To

2 1 To 1 3To 1 To
[9062.5 ( 2 sen wj ) − 9062.5 ( 2 sen wj − 2 sen wj )
To wj 4 wj 4 wj 4
1 1 3To
+ 9062.5 ( 2 sen wj To − 2 sen wj )]
wj wj 4

2 9062.5To2 j2π To 9062.5To2 j2π 3To 9062.5To2 j2π To

[ sen . − sen + sen
To (j2π)2 To 4 (j2π)2 To 4 (j2π)2 To 4
2 2
9062.5To j2π 9062.5To j2π 3To
+ 2
sen To − 2
sen ]
(j2π) To (j2π) To 4

9062.5𝑇𝑜 jπ 9062.5𝑇𝑜 j3π 9062.5𝑇𝑜 jπ 9062.5𝑇𝑜

2
sen − 2
sen + 2
sen + sen j2π
2(jπ) 2 2(jπ) 2 2(jπ) 2 2(jπ)2
9062.5𝑇𝑜 j3π
− 2
sen
2(jπ) 2

pág. 19
 Quedando:

9062.5𝑇𝑜 jπ 9062.5𝑇𝑜 j3π 9062.5𝑇𝑜

2
sen − 2
sen + sen j2π
(jπ) 2 (jπ) 2 2(jπ)2

9062.5𝑇𝑜 jπ 9062.5𝑇𝑜 j3π

2
sen − 2
sen
(jπ) 2 (jπ) 2

Si

j es par: 0

j es impar:
18125𝑇𝑜 𝟏𝟏𝟔𝟎𝟎
J = 1,5,9,13,17 = ⋯ . . etc = 4n − 3 = =
(jπ)2 (𝐣𝛑)𝟐

Siempre ocurrirá que:

jπ j3π
sen = 1; sen = −1
2 2
18125𝑇𝑜 𝟏𝟏𝟔𝟎𝟎
J = 3,7,11,15,19 = ⋯ . . etc = 4n − 1 = − 2
=−
(jπ) (𝐣𝛑)𝟐

Siempre ocurrirá que:

jπ j3π
sen = −1; sen =1
2 2

Finalmente:

∞
𝟏𝟏𝟔𝟎𝟎
𝑷(𝒕) = ∑ 𝐬𝐞𝐧 𝐰𝐣𝐭
(𝐣𝛑)𝟐
𝒋=𝟒𝒏−𝟑

∞
𝟏𝟏𝟔𝟎𝟎
𝑷(𝒕) = ∑ − 𝐬𝐞𝐧 𝐰𝐣𝐭
(𝐣𝛑)𝟐
𝒋=𝟒𝒏−𝟏…

pág. 20
 La respuesta desplazamiento será para un To=0.64s:

∞
𝟏𝟏𝟔𝟎𝟎 𝟏 (𝟏 − 𝜷𝒋𝟐 )𝒔𝒆𝒏𝒘𝒋𝒕 − 𝟐𝝃𝜷𝒋𝒄𝒐𝒔𝒘𝒋𝒕
𝒖(𝒕) = ∑ (
𝛑𝟐 𝟏𝟎(𝐣)𝟐 (𝟏 − 𝜷𝒋𝟐 )𝟐 + (𝟐𝝃𝜷)𝟐
𝒋=𝟒𝒏−𝟑
∞
−𝟏𝟏𝟔𝟎𝟎 𝟏 (𝟏 − 𝜷𝒋𝟐 )𝒔𝒆𝒏𝒘𝒋𝒕 − 𝟐𝝃𝜷𝒋𝒄𝒐𝒔𝒘𝒋𝒕
+ ∑ (
𝛑𝟐 𝟏𝟎(𝐣)𝟐 (𝟏 − 𝜷𝒋𝟐 )𝟐 + (𝟐𝝃𝜷)𝟐
𝒋=𝟒𝒏−𝟏

𝑬𝒋𝒆𝒓𝒄𝒊𝒄𝒊𝒐 𝑵𝒓𝒐 𝟎𝟒
Calcular la respuesta dinámica del tanque de agua de la Figura 5, el cual está sujeto a una
carga explosiva cuyo histograma de fuerza se muestra en la misma figura. w = # dígitos de su
primer nombre por 10.

𝑤 = 𝑥 𝑘𝑙𝑏 − 𝑓 𝑥 10

𝑥 = 𝐵𝐸𝑅𝐿𝑌 + 𝐸𝑅𝐼𝐶𝐾𝑆𝑂𝑁 + 𝐽𝑈𝐿𝐼𝑂 + 𝐾𝐸𝐿𝑉𝐼𝑁

𝒘 = 𝟐𝟒 𝒌𝒍𝒃 − 𝒇 𝒙 𝟏𝟎 = 𝟐𝟒𝟎 𝒌𝒍𝒃 − 𝒇 = 𝟐𝟒𝟎𝟎𝟎𝟎 𝒍𝒃-f

𝑘𝑙𝑏 𝑙𝑏 − 𝑓
𝑘 = 2700 = 2700000
𝑓𝑡 𝑓𝑡

240000 𝑙𝑏 − 𝑓
𝑚= = 7462.6865 𝑙𝑏 − 𝑓
𝑙𝑏
32.16 2
𝑠

pág. 21
 𝜺=𝟎

𝑘 0.5 𝑟𝑎𝑑
𝑤𝑛 = ( ) = 19.021
𝑚 𝑠

𝑟𝑎𝑑
𝑤𝑑 = 𝑤𝑛 (1 − 𝜀 2 )0.5 = 19.021
𝑠

𝑘 0.5 𝑟𝑎𝑑
𝑇𝑛 = ( ) = 19.021
𝑚 𝑠

𝑟𝑎𝑑
𝑇𝑑 = 2𝑥0.025 = 0.050
𝑠

𝑇𝑑
= 0.1513 < 0.25 𝑖𝑚𝑝𝑢𝑙𝑠𝑜 𝑑𝑒 𝑐𝑜𝑟𝑡𝑎 𝑑𝑢𝑟𝑎𝑐𝑖𝑜𝑛
𝑇𝑛

𝐼𝑚𝑝𝑢𝑙𝑠𝑜 = 𝐴𝑟𝑒𝑎 𝑑𝑒𝑏𝑎𝑗𝑜 𝑑𝑒 𝑙𝑎 𝑐𝑢𝑟𝑣𝑎

𝑃 𝑙𝑏 − 𝑓
𝐼𝑚𝑝𝑢𝑙𝑠𝑜 = 𝑇𝑑 ∗ = 6000
2 𝑠
𝐼𝑚𝑝 𝑥 2𝜋
𝑈𝑚𝑎𝑥 = = 0.04226 𝑓𝑡
𝑘 𝑥 𝑇𝑛
𝐹𝑠 = 𝑈𝑚𝑎𝑥 𝑥 𝑘 = 114126.24 𝑙𝑏 − 𝑓

pág. 22
Para un pulso triangular simétrico calculamos las 3 ecuaciones para su rango:

Donde el tiempo varia para los limites de td

t RD1 RD2 RD3 R dinamica
0.000 0 0
0.0125 0.0046984 0.0046984
0.025 0.03727 0.03727 0.03727
0.0375 0.1146222 0.1146222
Rangos de evaluacion
0.050 0.2136616 0.2136616 0.2136616
0.0625 0.3053756 0.3053756 0≤t≤0.125
0.075 0.3799047 0.3799047
0.0875 0.4330546 0.4330546 0.25≤t≤0.05
0.100 0.4618343 0.4618343
0.1125 0.4646243 0.4646243
0.125 0.4412676 0.4412676 0.05≤t
0.1375 0.3930786 0.3930786
0.150 0.3227691 0.3227691
0.1625 0.2342958 0.2342958
0.175 0.1326374 0.1326374
0.1875 0.0235149 0.0235149
0.200 -0.086931 -0.086931
0.213 -0.192485 -0.192485
0.2250 -0.287206 -0.287206
0.238 -0.365766 -0.365766
0.250 -0.423741 -0.423741
0.2625 -0.457871 -0.457871
0.275 -0.466234 -0.466234
0.288 -0.44836 -0.44836
0.3000 -0.405254 -0.405254
0.313 -0.339343 -0.339343
0.325 -0.254335 -0.254335
0.3375 -0.155014 -0.155014
0.350 -0.04697 -0.04697
0.3625 0.0637171 0.0637171
0.375 0.1708186 0.1708186
0.388 0.2683073 0.2683073
0.400 0.3506971 0.3506971
0.4125 0.4133513 0.4133513
0.425 0.4527442 0.4527442
0.438 0.466659 0.466659
0.4500 0.4543125 0.4543125
0.463 0.4163996 0.4163996
0.475 0.3550538 0.3550538
0.4875 0.2737274 0.2737274
0.500 0.1769969 0.1769969

pág. 23
𝑬𝒋𝒆𝒓𝒄𝒊𝒄𝒊𝒐 𝑵𝒓𝒐 𝟎𝟓
En la figura, se muestra un pórtico, cuyas columnas son C1, C2y C3, E= 2.25x107 KN/m2. asumir
que las vigas tienen rigidez a la flexión infinita. Considere el peso estimado de cada nivel de
P=X ton/m2, donde X=número de letras de su segundo apellido. Determinar:
Las matrices de rigidez y de masa, Los periodos, frecuencias y modos. Normalizar y graficar los
modos respecto al desplazamiento del ultimo nivel.

Definición de datos

𝑃 = 𝑋 𝑡𝑜𝑛/𝑚2

𝑄𝑈𝐼𝑆𝑃𝐸 + 𝐹𝐿𝑂𝑅𝐸𝑆 + 𝐼𝑁𝐶𝐴𝐻𝑈𝐴𝑀𝐴𝑁 + 𝑅𝐴𝑀𝑂𝑆

𝑋= ≅7
4

𝑷 = 𝟕 𝒕𝒐𝒏/𝒎𝟐
𝐾𝑁
𝐸𝑐 = 2.25 𝑥 107
𝑚2

pág. 24
Solución
Idealización

Hallamos las masas de cada piso, ya que nos da el peso de cada nivel para lo cual
hacemos un pequeño metrado de cargas.

𝑡𝑜𝑛
(7 𝑚2 𝑥 ((8 + 6) ∗ (6.20 + 6.50 + 6.20)))
𝑚1 = 𝑚2 = 𝑚3 =
9.81

𝒎𝟏 = 𝒎𝟐 = 𝒎𝟑 = 𝟏𝟖𝟖. 𝟖𝟏 𝑻𝒏 − 𝒇 𝒔𝟐 /𝒎

𝑡𝑜𝑛
7 𝑚2 𝑥 ((6) ∗ (6.20 + 6.50 + 6.20))
𝑚4 =
9.81

𝒎𝟒 = 𝟖𝟎. 𝟗𝟐 𝑻𝒏 − 𝒇 𝒔𝟐 /𝒎

Compatibilizamos unidades de la Elasticidad del concreto:

𝐾𝑁 1 𝑇𝑛
𝐸𝑐 = 2.25𝑥107 ∗( )
𝑚2 9.807 𝐾𝑁

𝑇𝑛
𝐸𝑐 = 2.294 𝑥 106
𝑚2

pág. 25
 a) Armamos la matriz de masa.

b) Inercia en columnas: Debemos de tener en cuenta la dirección del Sismo, para este
caso es en dirección X.

Para columna C1

𝑰𝑿𝑿𝑪 = 𝟑. 𝟒 𝒙𝟏𝟎−𝟑 𝒎𝟒

Para columna C2

0.60 𝑥 0.303
𝐼𝑋𝑋𝐶 =
12

𝑰𝑿𝑿𝑪 = 𝟏. 𝟑𝟓 𝒙𝟏𝟎−𝟑 𝒎𝟒

pág. 26
Para columna C3

0.50 𝑥 0.503
𝐼𝑋𝑋𝐶 =
12

𝑰𝑿𝑿𝑪 = 𝟓. 𝟐𝟎𝟖 𝒙𝟏𝟎−𝟑 𝒎𝟒

c) Hallando la rigidez por fórmula (Debemos de tener en cuenta que tenemos 4

columnas en cada dirección y-y)

-Para primer piso: K1er Piso

12𝐸𝐼
𝐾1𝑒𝑟 𝑝𝑖𝑠𝑜 = 4𝑥
ℎ13

12𝑥2.294 𝑥 106 𝑥(3.40 + 1.35 + 5.208)𝑥10−3

𝐾1𝑒𝑟 𝑝𝑖𝑠𝑜 = 4𝑥
4.203

𝑲𝟏𝒆𝒓 𝒑𝒊𝒔𝒐 = 𝟏. 𝟒𝟖𝟎𝒙𝟏𝟎𝟒 𝑻𝒏 − 𝒇/𝒎

-Para segundo y tercer piso: K2do y K3er Piso, ya que están a una misma altura.

12𝐸𝐼
𝐾2𝑑𝑜 𝑝𝑖𝑠𝑜 = 4𝑥
ℎ3

12𝑥2.294 𝑥 106 𝑥(3.40 + 1.35 + 5.208)𝑥10−3

𝐾2𝑑𝑜 𝑝𝑖𝑠𝑜 = 4𝑥
3.003

𝑲𝟐𝒅𝒐 𝒑𝒊𝒔𝒐 = 𝑲𝟑𝒆𝒓 𝒑𝒊𝒔𝒐 = 𝟒. 𝟎𝟔𝟏𝒙𝟏𝟎𝟒 𝑻𝒏 − 𝒇/𝒎

pág. 27
-Para cuarto piso: K4to piso.

12𝐸𝐼
𝐾4𝑡𝑜 𝑝𝑖𝑠𝑜 = 4𝑥
ℎ43

12𝑥2.294 𝑥 106 𝑥(1.35 + 5.208)𝑥10−3

𝐾4𝑡𝑜 𝑝𝑖𝑠𝑜 = 4𝑥
3.003

𝑲𝟒𝒕𝒐 𝒑𝒊𝒔𝒐 = 𝟐. 𝟔𝟕𝟒𝒙𝟏𝟎𝟒 𝑻𝒏 − 𝒇/𝒎

d) Aplicando rigidez gráfica

Cuando u1=1 y u2=u3=u4=0

pág. 28
Cuando u2=1 y u1=u3=u4=0

Cuando u3=1 y u1=u2=u4=0

pág. 29
Cuando u4=1 y u1=u2=u3=0

e) Ensamblando la matriz de rigidez de las columnas

Reemplazando valores y operando:

 1.480 104 + 4.061 104 −4.061 10

4
0 0 
 4 4 4 4 
K =
 −4.061 10 4.061 10 + 4.061 10 −4.061 10 0 
 4 4 4 4
 0 −4.061 10 4.061 10 + 2.674 10 −2.674 10 
 4 4 
 0 0 −2.674 10 2.674 10 

pág. 30
  5.541  104 −4.061  10
4
0 0 
 4 4 4 
K=
 −4.061  10 8.122  10 −4.061  10 0 
 4 4 4
 0 −4.061  10 6.735  10 −2.674  10 
 4 4 
 0 0 −2.674  10 2.674  10 

f) Frecuencias y modos naturales de vibración

-Frecuencias Naturales:

|𝐾 − 𝑊𝑛2 𝑚| = 0

Donde:

𝑊𝑛2 = 𝜆

 5.541  104 −4.061  10

4
0 0 
   188.81 0 0 0 
4 4 4
 −4.061  10 8.122  10 −4.061  10 0  −   0 188.81 0 0 
=0
 4 4 4  0 0 188.81 0 
 0 −4.061  10 6.735  10 −2.674  10   
 4 4   0 0 0 80.92 
 0 0 −2.674  10 2.674  10 

 5.541  104 − 188.81  −4.061  10

4
0 0
 
8.122  10 − 188.81 
4 4 4
 −4.061  10 −4.061  10 0  =0
 
6.735  10 − 188.81 
4 4 4
 0 −4.061  10 −2.674  10 
 4
2.674  10 − 80.92  
4 
 0 0 −2.674  10

Hallamos la determinante.

𝐷𝑒𝑡 = 5.447𝑥108 𝜆4 − 7.684𝑥1011 𝜆3 + 3.279𝑥1014 𝜆2 − 4.222𝑥1016 𝜆 + 6.527𝑥1017 = 0

pág. 31
 Mediante la HP50g u otro programa que resuelva ecuación de 4to grado, obtenemos las
soluciones:

𝝀𝟏 = 17.825 → 𝑾𝒏𝟏 = √17.825 = 𝟒. 𝟐𝟐𝟐 𝒓𝒂𝒅/𝒔

𝝀𝟐 = 196.539 → 𝑾𝒏𝟐 = √196.539 = 𝟏𝟒. 𝟎𝟏𝟗 𝒓𝒂𝒅/𝒔

𝝀𝟑 = 472.600 → 𝑾𝒏𝟑 = √472.600 = 𝟐𝟏. 𝟕𝟑𝟗 𝒓𝒂𝒅/𝒔

𝝀𝟒 = 723.719 → 𝑾𝒏𝟒 = √723.719 = 𝟐𝟔. 𝟗𝟎𝟐 𝒓𝒂𝒅/𝒔

Hallamos los periodos

2𝜋
𝑻𝟏 = → 𝟏. 𝟒𝟖𝟖 𝒔
4.222

2𝜋
𝑻𝟐 = → 𝟎. 𝟒𝟒𝟖 𝒔
14.019

2𝜋
𝑻𝟑 = → 𝟎. 𝟐𝟖𝟗 𝒔
21.739

2𝜋
𝑻𝟒 = → 𝟎. 𝟐𝟑𝟒 𝒔
26.902

Modos de Vibración

-Modo 1 𝜆1 = 17.825 , 𝑊𝑛1 = 4.222 𝑟𝑎𝑑/𝑠

 5.541 4
 10 − 188.81  17.825 −4.061  10
4 
0 0
11
    
0

  21  =  0 
4 4 4
 −4.061  10 8.122  10 − 188.81  17.825 −4.061  10 0
 4 4 4   31   0 
 0 −4.061  10 6.735  10 − 188.81  17.825 −2.674  10    
 4 4   41   0 
 0 0 −2.67  10 2.674  10 − 80.92  17.825 

pág. 32
 𝐶𝑜𝑛𝑠𝑖𝑑𝑒𝑟𝑎𝑛𝑑𝑜 ∅41 = 1 𝑦 𝑟𝑒𝑠𝑜𝑙𝑣𝑖𝑒𝑛𝑑𝑜 𝑡𝑒𝑛𝑒𝑚𝑜𝑠:

 5.204  104 4 
−4.061  10 0 0
 11   0 
 
 21
4 4 4
P=
 −4.061  10 7.785  10 −4.061  10 0   =
0
 4 4 4  31  0
 0 −4.061  10 6.398  10 −2.674  10 
 4 
   
 0 0 −2.674  10
4
2.53  10   1  0

−2.674𝑥104 ∅31 + 2.530𝑥104 (1) = 0

∅𝟑𝟏 = 𝟎. 𝟗𝟒𝟔

−4.061𝑥104 ∅21 + 6.398𝑥104 (0.946)

= 2.674𝑥104
∅𝟐𝟏 = 𝟎. 𝟖𝟑𝟐

5.204𝑥104 ∅11 − 4.061𝑥104 (0.832) = 0

∅𝟏𝟏 = 𝟎. 𝟔𝟒𝟗

0.649
Modo 1 = [0.832]
0.946
1

-Modo 2 𝜆2 = 196.539 , 𝑊𝑛2 = 14.019 𝑟𝑎𝑑/𝑠

 5.541 4
 10 − 188.81  196.539 −4.061  10
4 
0 0
12
    
0

  22  =  0 
4 4 4
 −4.061  10 8.122  10 − 188.81  196.539 −4.061  10 0
 4 4 4   32   0 
 0 −4.061  10 6.735  10 − 188.81  196.539 −2.674  10    
 4 4   42   0 
 0 0 −2.674  10 2.674  10 − 80.92  196.539 

pág. 33
 𝐶𝑜𝑛𝑠𝑖𝑑𝑒𝑟𝑎𝑛𝑑𝑜 ∅42 = 1 𝑦 𝑟𝑒𝑠𝑜𝑙𝑣𝑖𝑒𝑛𝑑𝑜 𝑡𝑒𝑛𝑒𝑚𝑜𝑠:

 1.83  104 −4.061  10

4
0 0   12   0 
 
P=
 −4.061  10
4 4
4.411  10 −4.061  10
4
0   22  =
0
 4
−4.061  10 3.024  10
4 4
−2.674  10   32  0
 0
 4 
   
 0 0 −2.674  10
4
1.084  10   1  0

−2.674𝑥104 ∅32 + 1.084𝑥104 (1) = 0

∅𝟑𝟐 = 𝟎. 𝟒𝟎𝟔

−4.061𝑥104 ∅22 + 3.024𝑥104 (0.406)

= 2.674𝑥104
∅𝟐𝟐 = −𝟎. 𝟑𝟓𝟕

1.83𝑥104 ∅12 − 4.061𝑥104 (−0.357) = 0

∅𝟏𝟐 = 𝟎. 𝟕𝟗𝟑

0.793
Modo 2 = [−0.357]
0.406
1

-Modo 3 𝜆3 = 472.600 , 𝑊𝑛3 = 21.739 𝑟𝑎𝑑/𝑠

 5.541 4
 10 − 188.81  472.600 −4.061  10
4 
0 0
13
    
0

  23  =  0 
4 4 4
 −4.061  10 8.122  10 − 188.81  472.600 −4.061  10 0
 4 4 4   33   0 
 0 −4.061  10 6.735  10 − 188.81  472.600 −2.674  10    
 4 4   43   0 
 0 0 −2.674  10 2.674  10 − 80.92  472.600 

pág. 34
 𝐶𝑜𝑛𝑠𝑖𝑑𝑒𝑟𝑎𝑛𝑑𝑜 ∅43 = 1 𝑦 𝑟𝑒𝑠𝑜𝑙𝑣𝑖𝑒𝑛𝑑𝑜 𝑡𝑒𝑛𝑒𝑚𝑜𝑠:

 −3.382  104 −4.061  10

4
0 0   13   0 
 
P=
 −4.061  10
4 3
−8.012  10 −4.061  10
4
0   23  =
0
 4
−4.061  10 −2.188  10
4 4
−2.674  10 
 33  0
 0
 4 
   
 0 0 −2.674  10
4
−1.15  10   1  0

−2.674𝑥104 ∅33 − 1.150𝑥104 (1) = 0

∅𝟑𝟑 = −𝟎. 𝟒𝟐𝟗

−4.061𝑥104 ∅23 − 2.188𝑥104 (−0.429)

= 2.674𝑥104
∅𝟐𝟑 = −𝟎. 𝟒𝟐𝟖

−3.382𝑥104 ∅13 − 4.061𝑥104 (−0.428) = 0

∅𝟏𝟑 = 𝟎. 𝟓𝟏𝟒

0.514
Modo 3 = [−0.428]
−0.429
1

-Modo 4 𝜆4 = 723.719 , 𝑊𝑛4 = 26.902 𝑟𝑎𝑑/𝑠

 5.541 4
 10 − 188.81  723.719 −4.061  10
4 
0 0
14
    
0

  24  =  0 
4 4 4
 −4.061  10 8.122  10 − 188.81  723.719 −4.061  10 0
 4 4 4   34   0 
 0 −4.061  10 6.735  10 − 188.81  723.719 −2.674  10    
 4 4   44   0 
 0 0 −2.674  10 2.674  10 − 80.92  723.719 

pág. 35
𝐶𝑜𝑛𝑠𝑖𝑑𝑒𝑟𝑎𝑛𝑑𝑜 ∅44 = 1 𝑦 𝑟𝑒𝑠𝑜𝑙𝑣𝑖𝑒𝑛𝑑𝑜 𝑡𝑒𝑛𝑒𝑚𝑜𝑠:

 −8.124  104 
 14
4
−4.061  10 0 0   0 
 
 24
4 4 4
P=
 −4.061  10 −5.543  10 −4.061  10 0   0
 4 =
 0
4
−4.061  10 −6.93  10
4
−2.674  10   34  0
 4    
 0 0 −2.674  10
4
−3.182  10   1  0

−2.674𝑥104 ∅34 − 3.182𝑥104 (1) = 0

∅𝟑𝟒 = −𝟏. 𝟏𝟖𝟗

−4.061𝑥104 ∅24 − 6.93𝑥104 (−1.189)

= 2.674𝑥104
∅𝟐𝟒 = 𝟏. 𝟑𝟕𝟐

−8.124𝑥104 ∅14 − 4.061𝑥104 (1.372) = 0

∅𝟏𝟒 = −𝟎. 𝟔𝟖𝟓

−0.685
Modo 4 = [ 1.372 ]
−1.189
1

pág. 36