T4 Adema PDF
T4 Adema PDF
T4 Adema PDF
Escuela de Ingenierı́a
Departamento de Ingenierı́a Estructural y Geotécnica
ICE 3743 Análisis Sı́smico
Tarea 4
t
f (t) = p0
τ
, que se repite como muestra la Figura 1
Determine la expansión de la serie de Fourier de la función. Grafique la serie usando: 1 (valor constante),
2, 3, 4, y 10 términos.
∞
X ∞
X
p(t) = a0 + an cos(nΩt) + an sin(nΩt)
n=1 n=1
, donde
2π
Ω =
T1
Z τ +T1
1
a0 = p(t)dt
T1 τ
Z τ +T1
2
an = p(t)cos(nΩt)dt
T1 τ
Z τ +T1
2
bn = p(t)sin(nΩt)dt
T1 τ
p0
a0 =
2
−p0 ∗ (2sin(πn)2 − 2nπsin(2πn))
an =
2n2 π 2
2
2p0 ∗ ( sin(2πn)
4 − nπcos(2πn)
2 )
bn =
n π2
2
1
Figura 2: Serie de Fourier considerando n=1
2
Figura 4: Serie de Fourier considerando n=3
3
Figura 6: Serie de Fourier considerando n=10
Z τ Z τ
p̄ = p(t)cos(ωt)dt + p(t)sin(ωt)dt
0 0
4
Problema 2. Descomposición modal de un edificio en 3D
El siguiente problema se basa en un edificio 3D de tres pisos como muestra la Figura 8. Los muros de
la estructura pueden asumirse como rı́gidos en la dirección axial (vertical). Las losas son lo suficientemente
rı́gidas para asumir que se comportan como diafragmas rı́gidos. Es decir, cada losa tiene 3 grados de libertad:
dos traslaciones ux y uy , y un giro θz .
Figura 8: Edificio muro 3D de 3 pisos: (a) Vista isométrica 1; (b) Modelo de planta.
(4) Las losas tienen una masa sı́smica distribuida deq = 0,5 ton ton
m2 , excepto por el techo que tiene qroof = 0,3 m2 .
416 −192 0
−192 kN
kx = 384 −192
mm
0 −192 192
256 −112 0
−112 kN
ky = 208 −96
mm
0 −96 96
Para comenzar el problema, se leyó el anexo B del enunciado, el cual indica que para obtener la matriz
de rigidez de la estructura, se debe seguir el siguiente procedimiento.
5
Figura 9: Análisis simplificado de un edificio de múltiples pisos: (a) Losa número n con los marcos y su
distancia perpendicular al centro de masa; (b) Marco A (marco tipo en la dirección X); (c) Marco C (marco
tipo en la dirección Y )
M arcoA = [0, 6]; M arcoB = [0, −6]; M arcoC = [−3, 0]; M arcoD = [0, 0];
En este caso, se tiene que los DOFs del marco y losas del piso n
u1x
u2x
u3x
u1y
ua = I 0 −ya I u2y
u3y
θ1z
theta2z
theta3z
axa = I 0 −ya I
6
u1x
u2x
u3x
u1y
ua = 0 I xc I u2y
u3y
θ1z
theta2z
theta3z
ayc = 0 I xc I
Luego se ensambla la matriz K, sumando las matrices de rigidez transformadas de los marcos
X
K= ki
i
8,32·105 −3,84·105
0 0 0 0 0 0 0
−3,84·105 7,68·105 −3,84·105 0 0 0 0 0 0
−3,84·105 3,84·105
0 0 0 0 0 0 0
0 0 0 5,12·105 −2,24·105 0 −7,68·105 3,36·105 0 kN
K= −2,24·105 4,16·105 −1,92·105 3,36·105 −6,24·105 2,88·105
0 0 0
m
0 0 0 0 −1,92·105 1,92·105 0 2,88·105 −2,88·105
−7,68·105 3,36·105 3,23·107 −1,48·107
0 0 0 0 0
0 0 0 3,36·105 −6,24·105 2,88·105 −1,48·107 2,95·107 −1,47·107
0 0 0 0 2,88·10 −2,88·105
5
0 −1,47·107 1,47·107
7
donde
b2 + h2
I0n = mn
12
Por lo tanto,
260,000 0 0 0 0 0 0 0 0
0 260,000 0 0 0 0 0 0 0
0 0 156,000 0 0 0 0 0 0
0 0 0 260,000 0 0 0 0 0
0 0 0 0 260,000 0 0 0 0
M =
0
ton
0 0 0 0 156,000 0 0 0
0
0 0 0 0 0 2,33 · 104 0 0
0
0 0 0 0 0 0 2,33 · 104 0
4
0 0 0 0 0 0 0 0 1,40 · 10
(2.2) Determine las frecuencias naturales wn y sus modos asociados φn . Haga gráficas para los primeros 3
modos.
12,5654
16,2284
20,1320
33,5529
rad
wn 41,9186 s
=
45,1486
54,0716
56,1894
72,7073
0,5000
0,3872
0,3121
0,1873
Tn 0,1499 s
=
0,1392
0,1162
0,1118
0,0864
8
0 0 −0,022 0 0 0 0,049 0 −0,0314
0
0 −0,042 0 0 0 0,009 0 0,0443
0
0 −0,051 0 0 0 −0,048 0 −0,0386
−0,009 0,019 0 −0,021 0,036 −0,023 0 0,034 0
−0,017
V = 0,038 0 −0,006 0,018 0,017 0 −0,038 0
−0,021 0,048 0 0,022 −0,048 −0,007 0 0,028 0
−9,60 · 10−4
−0,002 0 −0,005 −0,003 −0,003 0 −0,002 0
−8,95 · 10−4 1,56 · 10−4
−0,004 −0,002 0 0,004 0 0,002 0
−0,005 −0,002 0 0,005 0,002 −0,004 0 −0,002 0
Ahora se procede a entregar los gráficos de los primeros 3 modos, amplificados 80 veces para que se note
la diferencia.
9
Figura 11: Modo 2
Considere amortiguamiento modal igual a 5 % en todos los modos. Use su función IntFiniteDiffMDOF() para
obtener la respuesta u(t).
10
Lo primero que se hizo fue calcular la matriz de amortiguamiento, la cual se obtiene usando amortigua-
miento modal, donde
X 2ξn wn
c = m( φn φTn )m
Mn
Luego al sumar la contribución de cada modo, se obtiene la matriz de amortiguamiento de la estructura
la que corresponde a
Con esto, y luego de cargar el registro en Matlab, se obtuvieron las siguientes respuestas para cada piso.
11
Figura 14: Desplazamiento en x, piso 3
12
Figura 16: Desplazamiento en y, piso 1
13
Figura 18: Desplazamiento en y, piso 3
14
Figura 20: Giro en Z, piso 2
Use la función IntExcitInterp() para obtener la respuesta modal qn (t) para todos los modos n = 1 ..9.
Compare sus desplazamientos (y rotaciones) con los resultados de la pregunta anterior.
Para este caso, se calculó el factor de participación modal de cada modo, ası́ como el amortiguamiento de
cada modo utilizando
Ln = φT M φ
c = 2ξwn
15
Además, se calculó la masa modal de cada modo, la cual corresponde a 1 para todos los modos. Por lo
tanto, el factor de participación modal corresponde a
Ln
Γn =
mn
Γn = Ln
Luego se utilizó la función IntExcitInterp() y se obtuvo la respuesta modal qn . Para obtener la respuesta
en coordenadas globales de la estructura, se hizo
X
u(t) = φn qn
Luego se graficó la respuesta encontrada con este método, con la respuesta encontrada con el tiempo
historia modal. Los gráficos se muestran a continuación.
16
Figura 23: Comparación Desplazamiento en x, piso 3
17
Figura 25: Comparación Desplazamiento en y, piso 1
18
Figura 27: Comparación Desplazamiento en y, piso 3
19
Figura 29: Comparación Giro en Z, piso 2
A priori, los gráficos parecen ser exactamente iguales. Sin embargo, al hacer un zoom al gráfico, se puede
apreciar una diferencia en el cuarto decimal, como se observa en la siguiente imagen.
20
Figura 31: Zoom Desplazamiento en x, piso 1
Las losas actúan como diafragmas rı́gidos: el desplazamiento de cualquier punto de un determinado nivel se
puede calcular a partir de los grados-de-libertad (DOF) de la losa. Determine los desplazamientos absolutos
máximos en x e y de las 4 esquinas de la losa superior (techo), e indı́quelos en un dibujo/bosquejo de la losa.
Explique las diferencias en los desplazamientos de las 4 esquinas (si las hay).
En la siguiente página, se entregan gráficos con los desplazamientos máximos obtenidos, para cada una de
las esquinas, amplificando los desplazamientos 20 veces para que se note en el gráfico. Además, se entregan
los valores de los desplazamientos máximos por esquina. Cabe mencionar, que estos desplazamientos no
necesariamente se generan en el mismo periodo de tiempo.
21
Esquina inferior izquierda:
dxmax = −0,0877m
dymax = −0,1335m
22
Esquina inferior derecha:
dxmax = 0,0909m
dymax = −0,1385m
23
Esquina superior derecha:
dxmax = 0,0976m
dymax = 0,1343m
24
Esquina superior izquierda:
dxmax = −0,0814m
dymax = 0,1302m
25