Pauta Control 2

CI42G Dinámica de Estructuras

Prof: Rubén Boroschek Krauskopf.
Aux: Francisco Hernández Prado.

Viernes 2 de Noviembre de 2007

P1. Para la estructura que se muestra en la figura. Determine la Ec. de movimiento del
sistema, para ello asuma que β1=1% y que β2=3%. Determine adicionalmente los períodos
naturales de la estructura.

A = 2m
EI = 1000 kgf ⋅ m
m=0
EA = GA = ∞
γ = 100 kgf m 2
r = 0.5 m
kgf ⋅ m
kθ = 3000
rad

Solución:

• Definición de GDL estáticos:

Los GDL dinámicos son los que poseen masa, en este caso los GDL 1 y 2. Se concluye que
hay que condensar la matriz de rigidez a los GDL dinámicos.
 • Matriz de Masa:

⎡1 0⎤
⎢ 2 ⎥⋅
γ ⋅ π ⋅ r 2 ⎡8.014 0 ⎤
[M ] = ⎢ =⎢
0.501⎥⎦
r ⎥
0 g ⎣ 0
⎣⎢ 4 ⎦⎥

• Matriz de Rigidez:

⎡12 EI 6 EI ⎤
⎢ 3 0 ⎥
⎢ Lc Lc 2 ⎥ ⎡1500 1500 0 ⎤
⎢ 6 EI 4 EI 4 EI 2 EI ⎥ ⎢ ⎥
[K ] = ⎢ 2 + ⎥ = ⎢1500 3000 500 ⎥
⎢ Lc Lc Lv Lv ⎥ ⎢ 0 500 4000 ⎥⎦
⎢ 2 EI 4 EI ⎥ ⎣
⎢ 0 + kθ ⎥
⎣ Lv Lv ⎦

• Matriz de Rigidez Condensada:

A partir de la matriz anterior es posible identificar las sub-matrices asociadas a la

condensación como:

⎡1500 1500 ⎤ ⎡ 0 ⎤
[ K aa ] = ⎢ ⎡⎣ K pp ⎤⎦ = [ 4000]
T
⎥ ; ⎡⎣ K ap ⎤⎦ = ⎡⎣ K pa ⎤⎦ = ⎢ ⎥ ;
⎣1500 3000 ⎦ ⎣500 ⎦

Luego la matriz de rigidez condensada a los GDL dinámicos es igual a:

−1 ⎡1500 1500 ⎤
[ K cond ] = [ K aa ] − ⎡⎣ K ap ⎤⎦ ⋅ ⎡⎣ K pp ⎤⎦ ⎡⎣ K pa ⎤⎦ = ⎢ ⎥
⎣1500 2937.5⎦

• Períodos de la Estructura:

Partir del problema de valores propios es posible determinar las frecuencias angulares al
cuadrado. De las cuales es posible determinar
 {ω } = eigenvals ([ M ] ⋅ [ K ]) = ⎛⎜ 5961.6 ⎞⎟ ⎢ rad
2 −1 ⎡ 90.1 2
⎤
n 2 ⎥
⎠⎣ ⎝ s ⎦
⎛ 9.492 ⎞ ⎡ rad ⎤ 2 ⋅ π ⎛ 0.662 ⎞
⇒ {ωn } = ⎜ ⎟⎢ ⇒ Tn = = seg
⎝ 77.211⎠ ⎣ s ⎦
⎥ ωn ⎜⎝ 0.081 ⎟⎠

• Matriz de Amortiguamiento vía Rayleigth:

Para obtener una matriz de amortiguamiento clásica vía Rayleigth se debe cumplir:

⎡ 1 ⎤
⎢ 9.492 ⎥
⎛ 0.01 ⎞ 1 9.492 ⎛a⎞
⎜ ⎟ = ⋅⎢ ⎥ ⋅⎜ ⎟
⎝ 0.03 ⎠ 2 ⎢ 1 77.211⎥ ⎝ ⎠
b
⎣⎢ 77.211 ⎥⎦
⎛ a ⎞ ⎛ 0.12167 ⎞
⇒⎜ ⎟=⎜ ⎟
⎝ b ⎠ ⎝ 0.00076 ⎠

Luego la matriz de amortiguamiento es igual a:

⎡ 2.110 1.135 ⎤
[C ] = a ⋅ [ M ] + b ⋅ [ K ] = ⎢ ⎥
⎣1.135 2.284 ⎦

Una forma alternativa sería emplear Penzien-Wilson. (todas las unidades están en kgf, m,
s).

Finalmente la ecuación de movimiento es:

⎡8.014 r ( t ) ⎫⎪ ⎡ 2.110 1.135 ⎤ ⎧⎪ r1 ( t ) ⎫⎪ ⎡1500 1500 ⎤ ⎧⎪ r1 ( t ) ⎫⎪ ⎧0⎫

0 ⎤ ⎧⎪ 
⎢ 0 ⎥ ⋅⎨ 1 ⎬+ ⎢ ⋅⎨ ⎬+ ⋅⎨ ⎬ = ⋅⎨ ⎬
⎣ r2 ( t ) ⎭⎪ ⎣1.135 2.284 ⎥⎦ ⎪⎩ r2 ( t ) ⎭⎪ ⎢⎣1500 2937.5⎥⎦ ⎩⎪r2 ( t ) ⎭⎪ ⎩0⎭
0.501⎦ ⎩⎪

(todas las unidades están en kgf, m, s).

P2. Para la estructura que se muestra en la figura:

i) Complete la matriz de formas modales. ( [ Φ ] ).

ii) Determine el corte basal máximo asumiendo un movimiento horizontal en la base
descrito por el espectro de desplazamientos definido en la figura. (Combinación
SSRS).
iii) Determine la máxima aceleración del piso inferior. ¿Relativa ó absoluta?.
iv) Determine el máximo esfuerzo de la biela del piso superior.

⎡ 0.826 0.414 ⎤ ⎧ 0.720 ⎫

[Φ ] = ⎢ Tn = ⎨ ⎬ seg k1 ≠ k 2
⎣ 0.563 ? ⎥⎦ ⎩0.227 ⎭

Solución:

i) Complete la matriz de formas modales. ( [ Φ ] ).

Primero determinamos la matriz de masa:

⎡300 0 ⎤ 1 kgf ⋅ s 2 ⎡30.612 0 ⎤ kgf ⋅ s 2

[M ] = ⎢ ⎥ ⋅ ⋅ = ⎢ 0 ⋅
⎣ 0 200 ⎦ 9.8 m ⎣ 20.408⎥⎦ m

Empleando la propiedad de ortogonalidad de modos se tiene que cumplir que:

 {φi } ⋅ [ M ] ⋅ {φ j } = 0
T

⎡30.612 0 ⎤ ⎧0.414 ⎫
⇒ {0.826 0.563} ⋅ ⎢ ⋅⎨ ⎬=0
⎣ 0 20.408⎥⎦ ⎩ x ⎭
⎧ 12.673 ⎫
⇒ {0.826 0.563} ⋅ ⎨ ⎬=0
⎩20.408 ⋅ x ⎭
⇒ 10.468 + 11.490 ⋅ x = 0 ⇒ x = −0.911

De esta forma la matriz de formas modales queda:

⎡0.826 0.414 ⎤
[Φ] = ⎢ ⎥
⎣ 0.563 −0.911⎦

ii) Determine el corte basal máximo asumiendo un movimiento horizontal en la base

descrito por el espectro de desplazamientos definido en la figura. (Combinación
SSRS).

Determinamos las masas modales:

⎡ 27.4 0 ⎤ kgf ⋅ s 2
[ Mm] = [Φ ] ⋅ [ M ] ⋅ [ Φ ] = ⎢
T

⎣ 0 22.2 ⎥⎦ m

Luego los factores de participación modal; como es un edificio de corte asumimos un vector
{r} = {1 1} .
T

{Lm } = [Φ ] ⋅ [ M ] ⋅ {r}
T

T
⎡0.826 0.414 ⎤ ⎡30.612 0 ⎤ ⎧1⎫ ⎧ 36.776 ⎫ kgf ⋅ s 2
=⎢ ⋅
⎥ ⎢ ⋅⎨ ⎬ = ⎨ ⎬
⎣ 0.563 −0.911⎦ ⎣ 0 20.408⎥⎦ ⎩1⎭ ⎩−5.918⎭ m

De igual forma se determina las masas modales efectivas como:

⎧ Lm 2 ⎫ ⎧49.411⎫ kgf ⋅ s 2
{Meff } = ⎨ ⎬=⎨ ⎬
⎩ Mm ⎭ ⎩ 1.579 ⎭ m

Se evalúan el espectro de desplazamiento para los períodos de la estructura:

 ⎧0.108⎫
Sd = ⎨ ⎬m
⎩0.034 ⎭

Determinamos las frecuencias angulares a partir de los períodos:

⎧ 2 ⋅ π ⎫ ⎧ 8.727 ⎫ rad
{ωn } = ⎨ ⎬=⎨ ⎬
⎩ Tn ⎭ ⎩27.679 ⎭ s

El Pseudo espectro de aceleración se puede estimar a partir del espectro de desplazamiento

como:

⎧ 8.225 ⎫ m
{PSa} = {Sd ⋅ ωn2 } = ⎨ ⎬ 2
⎩26.087 ⎭ s

Luego el corte basal máximo por modos será:

⎧ 406.633⎫
{Qm } = {Meff ⋅ PSa} = {Meff ⋅ Sd ⋅ ωn2 } = ⎨ ⎬ kgf
⎩ 41.190 ⎭

El corte basal máximo se obtiene combinando los máximos cortes por modos, (se emplea
SSRS):
Qmax = 406.6332 + 41.1902 = 408.714 kgf

iii) Determine la máxima aceleración del piso inferior. ¿Relativa ó absoluta?.

La máxima aceleración de los grados de libertad por modos se obtiene como:

Lmi
{vm } = φi ⋅ ⋅ PSai
Mmi
Lm1 ⎧0.826 ⎫ 36.776 ⎧9.152 ⎫ m
{v1 } = φ1 ⋅ ⋅ PSa1 = ⎨ ⎬⋅ ⋅ 8.255 = ⎨ ⎬ 2
Mm1 ⎩0.563⎭ 27.4 ⎩6.238⎭ s
Lm ⎧ 0.414 ⎫ −5.918 ⎧−2.879 ⎫ m
{v2 } = φ2 ⋅ 2 ⋅ PSa2 = ⎨ ⎬⋅ ⋅ 26.087 = ⎨ ⎬ 2
Mm2 ⎩−0.911⎭ 22.2 ⎩ 6.335 ⎭ s

Haciendo combinación según SSRS, se obtiene la máxima aceleración del piso inferior:

m
v2− max = 6.2382 + 6.3352 = 8.891
s2
 iv) Determine el máximo esfuerzo de la biela del piso superior.

El máximo desplazamiento de los grados de libertad por modos se obtiene como:

Lmi
{vm } = φi ⋅ ⋅ Sdi
Mmi
Lm1 ⎧0.826 ⎫ 36.776 ⎧0.120 ⎫
{v1 } = φ1 ⋅ ⋅ Sd1 = ⎨ ⎬⋅ ⋅ 0.108 = ⎨ ⎬m
Mm1 ⎩0.563⎭ 27.4 ⎩0.082 ⎭
Lm ⎧ 0.414 ⎫ −5.918 ⎧ −0.00375⎫
{v2 } = φ2 ⋅ 2 ⋅ Sd 2 = ⎨ ⎬⋅ ⋅ 0.034 = ⎨ ⎬m
Mm2 ⎩−0.911⎭ 22.2 ⎩ 0.00826 ⎭

Luego el desplazamiento relativo máximo del piso superior por modos es:

{∆v1 } = 0.120 − 0.082 = 0.038m

{∆v2 } = −0.00375 − 0.00826 = −0.012m
Combinando ambos desplazamientos modales se obtiene:

{∆v ( ) } =
max
1
0.0382 + 0.0122 = 0.0398m

El ángulo con respecto a la horizontal de la biela y su largo son:

⎛ 2.5 ⎞
α = a tan ⎜ ⎟ ⇒ α = 39.8º Lbiela = 32 + 2.52 = 3.905m
⎝ 3 ⎠

Luego el máximo esfuerzo de la biela será su rigidez.

EA
Fbiela = ⋅ cos(39.8º ) ⋅ 0.0398m = 6.271 kgf
Lbiela
P3. Para la estructura que se muestra en la figura, determine:

i) Matriz de masa ( [M ] ) y de Rigidez ( [K ] ).

ii) Formas Modales ( [Φ] ) normalizadas por la masa modal y períodos {Tn } .
iii) Vector de influencia sísmico ( {r } ) ante un sismo horizontal en la base.

kgf
K H = 2000
m
kgf
Kv = 4000
m
kgf ⋅ m
Kθ = 1000
rad
kgf
γ acero = 7850 3
m

Solución:

i) Matriz de masa ( [M ] ) y de Rigidez ( [K ] ).

Se definen los GDL, en este caso se colocan en el centro de masas, para tener una matriz de
masa desacoplada. (La solución puede ser hecha imponiendo 3 GDL en cualquier posición).
Luego la matriz de masas es igual a:

⎡ ⎤
⎢1 0 0 ⎥ ⎡6.729 0 0 ⎤
⎢ ⎥ γ acero ⋅ a ⋅ b ⋅ e ⎢
[ M ] = ⎢0 1 0 ⎥⋅ =⎢ 0 6.729 0 ⎥⎥
⎢ 2 ⎥
g
⎢⎣ 0 0 1.082 ⎥⎦
⎢0 0 a + b ⎥
2

⎢⎣ 12 ⎦ ⎥

Donde a = 0.7 m , b = 1.2 m y e = 0.01m . (Unidades en kgf, m y seg).

La matriz de Rigidez será:

⎡ KH ⋅ b ⎤
⎢ KH 0 ⎥
2
⎢ ⎥ ⎡ 2000 0 1200 ⎤
KV ⋅ a ⎥=⎢ 0
[K ] = ⎢⎢ 0 Kv ⎥ ⎢ 4000 1400 ⎥⎥ (Unidades en kgf, m y
2
⎢ ⎥ ⎢1200 1400 2210 ⎥⎦
⎢ KH ⋅ b Kv ⋅ a Kv ⋅ a 2 K H ⋅ b 2 ⎥ ⎣
Kθ + +
⎢⎣ 2 2 4 4 ⎥⎦
seg).

ii) Formas Modales ( [Φ] ) normalizadas por la masa modal y períodos {Tn } .

Períodos de la estructura:

⎛ 2298.9 ⎞
( ) ⎡ ⎤
2

{ωn } = eigenvals [ M ] ⋅ [ K ] = ⎜ 114.9 ⎟⎟ ⎢ rad

2 −1 ⎜
2 ⎥
⎜ 490.1 ⎟ ⎣ s ⎦
⎝ ⎠
⎛ 47.947 ⎞ ⎛ 0.131 ⎞
⎜ ⎟ ⎡ rad ⎤ 2 ⋅π ⎜ ⎟
⇒ {ωn } = ⎜ 12.039 ⎟ ⎢ ⎥ ⇒ Tn = = ⎜ 0.522 ⎟ seg
⎜ 22.137 ⎟ ⎣ s ⎦ ωn ⎜ ⎟
⎝ ⎠ ⎝ 0.284 ⎠

Formas modales:

⎡-0.088 -0.728 0.383 ⎤

[Φ] = ⎢⎢-0.121 -0.288 -0.826 ⎥⎥ (los que me dio mi calculadora).
⎢⎣-0.989 0.622 0.414 ⎥⎦
Masas Modales:

⎡1.208 0 0 ⎤
⎢ 0 ⎥⎥ (los que me dio mi calculadora).
[Mm ] = [Φ] ⋅ [ M ] ⋅ [ Φ ] = ⎢0
T
4.545
⎢⎣0 0 5.759 ⎥⎦

Luego las formas modales normalizadas son:

{φi }
{φi} =
i
Mmi
⎡-0.080 -0.342 0.160 ⎤
⇒ ⎣Φ ⎦ = ⎢⎢-0.110 -0.135 -0.344 ⎥⎥
⎡ i ⎤
⎢⎣-0.900 0.292 0.173 ⎥⎦

iii) Vector de influencia sísmico ( {r } ) ante un sismo horizontal en la base.

⎛ 1⎞
{r } = ⎜⎜ 0 ⎟⎟
⎜0⎟
⎝ ⎠