Forma Modales
Forma Modales
Forma Modales
son solamente las que generan las masas mx y m2 al moverse lateralmente, entonces en dinmica se habla de un sistema de dos grados de
libertad, que son precisamente los desplazamientos laterales 1 y 2.
Esto no implica que en los restantes grados de libertad los giros y
desplazamientos correspondientes se anulen; adems, la matriz de rigideces de la estructura, que sera de 10 x 10, se pueden transformar
a una de 2 x 2 (expresada en funcin de los grados de libertad 1 y 2),
denominada matriz de rigideces lateral, mediante el proceso de condensacin esttica descrito tambin en la seccin 1.2.1, vase la
expresin 1.19.
En lo que sigue de este captulo al hablar de cierto nmero de grados de libertad, se alude slo a aquellos en que existen fuerzas generalizadas de inercia.
Como el ttulo del captulo lo seala, aqu se tratan slo las ideas
fundamentales. Existen presentaciones mucho ms completas de la
dinmica estructural, por ejemplo, las hechas en las referencias 35
a 39.
5.
Conceptos fundamentales
de dinmica estructural
Desde el punto de vista dinmico, los grados de libertad que interesan son aquellos en los que se consideran fuerzas generalizadas de
inercia; es decir, fuerzas iguales a masa por aceleracin y momentos
iguales al momento de inercia de masa por aceleracin angular. Por
ejemplo, en la figura 5.1 se muestra un marco que, de acuerdo con la
definicin de grados de libertad dada en la seccin 1.2.1 de este manual, tiene 10 grados de libertad si se ignoran las deformaciones
ai
Elemento elstico
Base
Amortiguador
Visco so.
fuerzas de rigidez y de amortiguamiento son, respectivamente, proporcionales al desplazamiento u y a la velocidad u de la masa con
respecto a su base. Sean k y c las correspondientes constantes de proporcionalidad, que se supone que no cambian con el tiempo. Este
conjunto constituye un sistema lineal de un grado de libertad, con
amortiguamiento viscoso o lineal.
En estas circunstancias, usando el principio de D'Alambert la
ecuacin de equilibrio dinmico es
mx+c
Considerando que *
escnbir
+ ku = 0
= s0 +
m+c
M,
+ ku=-ms0
(5.1)
s0
(5.2)
y
5,3 Vibraciones libres de un s i s t e m a
1 de la figura 5 . 2 .
! 0 + w W0
T)
'i
(5.3)
(5>4)
Cuando la base del sistema est en movimiento, es necesario resolver la ecuacin 5.2, considerando no nulo el segundo trmino.
La solucin est dada por la siguiente expresin, que proporciona
el valor de u en un instante t (referencia 37).
1
t
u(t) = jjra J 0 S o 0") exp { - a) (/ - r)} sen w c (/ - r) dr
"
1PerOd
'
SC
btene
61 e S P C C t r 0 d e l a r e s p u e s l a
en
1
3E
HS.
, '
U,.
I "JNJ
pfp
tC
tir.j""
t z i
En la figura 5.5 se presentan los espectros de la aceleracin absoluta (Sq + u) y de la velocidad relativa del temblor mencionado
Es frecuente obtener el espectro de desplazamiento D y en lugar de
las velocidades y aceleraciones dibujar las cantidades V = D y
A D, que se denominan espectros de seudovelocidades y de
seudoaceleraciones, respectivamente. De acuerdo con su definicin
es posible trazar estos tres ltimos espectros (D, V y A) en una sola
grfica, con rayado logartmico en cuatro direcciones, segn se hace
en la figura 5.5. Este tipo de grfica se debe a Newmann y ha sido
empleado en muchas publicaciones (referencias 40, 41 y 42) Es interesante notar que la fuerza mxima en el elemento elstico k es
igual a
k D = m D = m (J2 D = m A
m
Conceptos fundamentales
de dinmica estructural3
< <oN
O
o
+ s 0 ), m 2 (u2 + s 0 ) y r
("3 + ^o-
PERIODO
(SEG)
PERIODO
(SEG)
* , figura 5.4 ( ^ o s
0, 2, 5,
Conceptos fundamentales
de dinmica
estructural
i (')"]
(5.11)
Ea = Cu
u"3
2
(( 0l )
JI
Zi
z2
, Z
+ Cu + Ku = Ms0
(5.12)
y sustituyendo 5 . 1 4 y 5 . 1 5 e n 5 . 1 3 s e llega a:
m1
m2
o l
m3J
M Z6
+ KZO
= 0
rmzi'
+ &kijZ,)d
= 0
de donde
- [i]
u =
L =
9
r* <i
0
So
o =
(5.16)
(5.15)
^ kU z
mz
(5.17)
So
0 + w2 6 =
cuya solucin es
6 = o s e n w ( / - r)
M + &
m 0
De acuerdo con lo anterior existen modos de vibracin que satisfacen las condiciones de la expresin 5.15. Estos son tales que el movimiento de cada masa es armnico simple con periodo natural
T = 2ir/u; w se llama frecuencia natural circular.
Derivando dos veces la ecuacin 5.18 se tiene
-ttf 2 a sen
= Z O (i)
(5.18)
(5.13)
(5.14)
- r) = -w2
(5.19)
(5.20)
5.3.5 Ejemplo
'20
re
Considrese la estructura mostrada en la figura 5.7, tomada de la referencia 74. Las matrices de masas y de rigideces de esta estructura
son:
Presenta
a , v l e t e ( 2 treS ^ d C a S f 1 3 f l 8 U r a 5 - 6 ) C U y a s o i u c i * i n d u c e
a /i valores de , es decir a * frecuencias de vibracin
que corresponden a otros tantos periodos naturales 2tt/ w
Los valores de a,* son reales y positivos, y sus races cuadradas son
las frecuencias naturales. Se acostumbra numerar a las w en orden
creciente, es decir la primera frecuencia co (llamada frecuencia fundamental) es el menor valor, y la ltima
el mayor.
Si cada valor de la frecuencia co, se reemplaza en 5.19 es posible obtener valores Z, diferentes de cero (cada uno de estos vectores se llama
modo de vibracin). Para cada modo no se obtienen soluciones nies decir que no estn
cas sino solamente valores relativos entre las
definidas las amplitudes de las vibraciones de las masas, sino las relaciones entre todas ellas.
Se demuestra que los modos de vibracin tienen las siguientes propiedades:
K =
= Osi jr
= o si j*r
1
0
m3J
kx 0+ k2 0
k2
- k 2 k2 + k
-
0
-
k3
k3
W 3 = 200
*3=
80
w 2 = 400
Arj = 200
(5.21)
W\ = 400
Z]KZr
m,
kj=
(5.22)
200
W-
K = 80
(400 -
2.5
3.5
- 1.0
.o
i.o
5.0 -
0.407750
o'
80
122 x 0.40775)
0.0
2.5
200
(280 -
2.5
3.5 -
0.0
- 0.0
haciendo y =
w2
0.40775
1.0
80
1.0
1.0 - 0.203875
122 x
(80 -
80
200\.
- 200
Zn
Zn
+ 230.2545
- 80
y2 =
1.525
7.030
80
Z21 =
Z31
122 x 0.203875)
Z21
Z21
Z21
80
+ 55.1273
Tx = 0.5686 seg
T2 = 0.2650seg
T3 = 0.1694 seg
a?M)Z = 0
= 0
Z31 = 0
Z31 = 0
Z, =
1.000
1.751
2.541
0.853
- 1.969
respectiva-
1.000
=
- 0.804
0.321
wf = 1375.2,
1.000
2
y, = 17.190
562.4,
<4 =
Zn
184.386 = 0
122.0,
Zn
Z21
f =
80
80
ecuacin siguiente:
y3 - 25.751 .y2 + 157.885 y -
80
200
r,
Mtodos
frecuencias
de vibrar
Conceptos fundamentales
de dinmica
estructural
200
200
fton/cm)
{ 0.40775
0.71397
0
0.40775
0.203875
0.51805 }
Z \ K = { 1.000
1.751
2.541
400
200
[-
- 200
280
-
80
80 1
il
iiS
u,
<
llll4n*w
a'.
f'"
.
V
,ton-seg2
cm
'
1
2
3
4
5
6
X
F/co2
V/co2
AY/u,2
Y/co2
OJ2
1
2
3
4
5
6
X
F co2
V/u2
A Y/oj2
Y/co2
w2
1
2
3
4
5
6
X
F/oj2
V/u2
A Y/co2
Y/co2
co2
^ y v
v -
V V
KJ
0.408
0.408
0.204
1.000
0.408
2.000
0.816
3.000
0.612
1.836
0.00918
1.428
0.00714
0.612
0.00765
0.00918
109
0.01632
123
0.02397
125
1.000
0.408
1.780
0.726
2.610
0.532
1.664
0.00837
80j
xi"
ZlKZx
... A
Rengln
0.40775
0
1.258
0.00629
0.532
0.00665
0.00837
119
0.01466
121
0.2131
1.000
0.408
1.750
0.714
2.550
0.520
1.642
0.00821
122
1.234
0.00617
0.520
0.0065
0.00821
121.8
0.01438
1.000
1.752
S FX
0.02088
122.1
121.7
0.024475
0.000201
2.543
ft
121
_2
9868
= 2 f f / w = 0 . 5 6 8 6 seg
Los pasos en que consiste el mtodo (tabla 5.1) son los siguientes:
) Supngase una forma para el modo. Esta es la que aparece en el
rengln 1 de la tabla. Para comenzar, es usualmente apropiado
suponer valores iguales al nmero de orden del piso (de abajo
hacia arriba).
b) Obtngase la fuerza de inercia en cada masa correspondiente a
la configuracin supuesta. Esta fuerza sera
como se
desconoce u2 se calculan los productos MX = F/<J, que forman el segundo rengln de la tabla.
c) Con las fuerzas de inercia calclense las fuerzas cortantes en loi
entrepisos, tambin divididas entre w2, es decir se calcula V/c,
como se anota en el tercer rengln de la tabla.
d) A partir de las fuerzas cortantes y de las rigideces de entrepiso,
obtnganse las deformaciones de entrepiso tambin divididas
I3)
de vibrar
LM(Y/o)2
Xx.
Fx.
f2/k2.
x2.
g) Reptanse los pasos (d) a (/) con el tercer resorte y la tercera masa.
h) Continese el proceso hasta llegar a la ltima masa. Si se satisface
el equilibrio entre la fuerza cortante del ltimo resorte y la fuerza
de inercia de la ltima masa, la frecuencia escogida y las amplitudes calculadas corresponden a un modo natural de vibracin. Por
lo general se obtendr un residuo.
Representando en una grfica los residuos obtenidos contra los
distintos valores de w2 supuestos, se obtendr una curva cuyos ceros corresponden a las frecuencias naturales.
Ntese que un cambio de signo en los residuos correspondientes a
dos valores de w2 indica que hay una frecuencia comprendida entre
dichos valores y se puede interpolar, por ejemplo linealmente, para obtener una mejor aproximacin de la frecuencia buscada.
Cuando se est probando un valor de w2 suficientemente prximo al
correspondiente a un modo de vibrar (cuando el residuo es pequeo), se
encuentra que una aproximacin ms precisa de dicha frecuencia es (referencia 44)
^
. t
V A X
(5.24)
L FX
= J-MZ
(5.25)
= MX.
c) Calclese el vector Y resolviendo el sistema de ecuaciones siguientes (que proviene de la expresin 5.25).
X'
d) Si el vector Y es igual al vectorX. multiplicado por una constante, entonces se tiene una forma modal y la constante es igual a
1/w2. En la prctica se busca que Y sea aproximadamente igual a
una constante por Xyse calcula o?2 con la relacin siguiente (que
es una manera de escribrir el cociente de Rayleigh)
<J =
YrM
(5.27)
Si se considera que Y no es lo suficientemente parecida a X, se empieza otra vez en el paso a) con un vector X que sea proporcional a ZSe demuestra (referencia 45) que as el proceso converge rpidamente
al primer modo.
El mtodo sirve tambin para determinar modos superiores de vibracin si es que los pasos anteriores se aplican empleando en vez de K la
matriz K' con un corrimiento de origen, es decir
K' = K - n M
(5.28)
= +
por
r
X'
Y
T
Y
MY
(5.29)
0.40775
0
0
0
0.40775
0
4 0 0 - 2 0 0
200
280
0
0
0.20388
0
- 80
80
200
280
80
0 y\
80 y2
80
X
x'z
(en ton/cm)
200
= 2 ^
80
x{
200
y* = y i +
Grado de
libertad
Tabla 5.4 Mtodo de iteracin inversa con corrimiento (n = 490.5, converge al segundo modo).
0.40775
0.81550
0.009174
0.01631
0.40775
0.72498
0.008326
0.01461
Grado de
libertad
0.02396
0.02127
0.40775
0.71560
0.52091
0.008221
0.01440
0.02092
1.000
1.000
1.000
X'
0.40775
0.40775
0.20388
0.01198
0.009939
1.000
0.830
-2.149
'
Y
0.40775
0.33843
-0.43813
0.01453
0.01249
-0.02805
1.000
0.860
X'
0.40775
0.35067
0.01370
0.01167
0.02574
-1.931
(
-0.39368
-0.02698
-1.969
0.851
1.000
A
Nota: Los valores de X, salvo para la primera iteracin, son proporcionales a los de
de la iteracin anterior
X'=
Y
MX
- I -
w 2 = p + m (ecuacin 5.59)
Y my
=
J
YTMY
0.01370 X 0.40775+0.01167
(ecuacin 5.27)
P =
X 0.35067 + 0.02698
X 0.39368
p = 72.4,
!~485/2)'
em
nCeS ^ 0nvergencia
s e r el v a l o r
S U P e r
K' = K -
r 200 --
pMM ==
I -
200
200
80
80
0"
80
20
x{
*3
La solucin es:
2*i +
yi
-V2 = / i
5
4xj
200
-
200
' $
C
u,
I
c
. .i
1
a,
II
= ~ 4 J>2 -
20
Puede notarse que el mtodo de iteracin inversa da, para el primer modo, los mismos resultados que el mtodo de Newmark. De
hecho ambos mtodos son equivalentes, en este caso, al procedimiento conocido como de Stodola-Vianello, que consiste bsicamente en
iterar empleando la expresin siguiente (vase por ejemplo la referenca Z),
Z
K-^MZ
(5.30)
que es otra forma de escribir la ecuacin 5.19, y que equivale a resolver el sistema de ecuaciones 5.26.
Sin embargo, el mtodo de iteracin inversa se puede aplicar
cualesquiera que sean las matrices de masas y rigideces, y no slo a
sistemas sencillamente acoplados; adems, como se ha visto, empleado con corrimientos, sirve para calcular cualquier modo de vibrar
Por tales motivos el mtodo de iteracin inversa constituye la base
de varios afcontmos (como el de iteracin de subespacios y el de bsqueda del determinante) apropiados para computadoras. En la referencia 45, se trata con ms amplitud este mtodo y sus variantes, y se
presentan inclusive los correspondientes programas para compu-
(1.5)
4/2,
i
hx + h2
LK, C1
ZKn
A,C1
Mtodos
aproximados
Rx
24 E
2/ij +
h2
(1.6)
P a r a el segundo entrepiso:
S u p o n i e n d o las c o l u m n a s e m p o t r a d a s en la cimentacin
48 E
R2 =
4/2,
(177)
+ h2
h2 + h3
LKr
Cl
12
48 E
R2 =
h7
+ h,
(1.8)
2 h + h7
EA:C 2
EA:ti
P a r a entrepisos intermedios:
48E
Rn =
(1.9)
K
Ak
"
LKr
+>
LK
rm
+, h n
hD
LK
tn
En estas ecuaciones
rigidez del entrepiso en cuestin.,
rigidez (I/L) de las vigas desnivel sobre el enirepiso n.
7
rigidez (I/L) de las columnas del entrepiso n.
m, n, o ndices que identifican tres niveles consecutivos de abajo
hacia arriba.
hn
altura del entrepiso n.
Rn
X tn
Para el entrepiso superior, si se acepta que la cortante del penltimo piso es el doble que la del ltimo, se encuentra que es aplicable la
frmula para entrepisos intermedios, poniendo 2h m en vez de h m y
haciendo h0 = 0.
7
t