T4 Adema PDF

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

PONTIFICIA UNIVERSIDAD CATÓLICA DE CHILE

Escuela de Ingenierı́a
Departamento de Ingenierı́a Estructural y Geotécnica
ICE 3743 Análisis Sı́smico

Tarea 4

Andrés Adema Yusta


8 de Octubre de 2018
Problema 1. Dada una función rampa-periódica (o serrucho):

t
f (t) = p0
τ
, que se repite como muestra la Figura 1

Figura 1: Señal de serrucho

1.1) Serie de Fourier Real

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.

Se sabe que la Serie de Fourier real viene dada por


X ∞
X
p(t) = a0 + an cos(nΩt) + an sin(nΩt)
n=1 n=1

, donde


Ω =
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 τ

Luego al realizar la integración, se obtiene

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

Luego se grafica para n=1,2,3,4 y 10 términos, y se muestra a continuación

1
Figura 2: Serie de Fourier considerando n=1

Figura 3: Serie de Fourier considerando n=2

2
Figura 4: Serie de Fourier considerando n=3

Figura 5: Serie de Fourier considerando n=4

3
Figura 6: Serie de Fourier considerando n=10

1.2) Espectro de Fourier


Dibuje o bosqueje el espectro de Fourier para dicha función.

Se sabe que el espectro de Fourier viene dado po

Z τ Z τ
p̄ = p(t)cos(ωt)dt + p(t)sin(ωt)dt
0 0

cuya solución corresponde a


p0
p̄ =
ω

Figura 7: Espectro de Fourier

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.

2.1) Matrices de masa M y de rigidez K global

La estructura de la Figura 8 tiene las siguientes propiedades y supuestos.

(1) Los marcos de muro A y B tienen rigideces en su plano dadas por kx


(2) Los marcos de muro C y D tienen rigideces en su plano dadas por ky
(3) La rigidez de los muros en su eje débil es despreciable.

(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 .

(5) Las losas tienen dimensiones lx xly =26mx20m

 
416 −192 0
−192 kN
kx = 384 −192
mm
0 −192 192
 
256 −112 0
−112 kN
ky = 208 −96
mm
0 −96 96

Se pide determinar las matrices M y K global.

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 )

(1) Obtener las coordenadas del centro de gravedad de los muros

M arcoA = [0, 6]; M arcoB = [0, −6]; M arcoC = [−3, 0]; M arcoD = [0, 0];

(2) Obtener las matrices de transformación de los muros

En este caso, se tiene que los DOFs del marco y losas del piso n

una = unx − θnz ya

, lo que puede escribirse como

 
u1x
 u2x 
 
 u3x 
 
  u1y 
 

ua = I 0 −ya I  u2y 


 u3y 
 
 θ1z 
 
theta2z 
theta3z

Luego la matriz de transformación para los muros horizontales corresponde a

 
axa = I 0 −ya I

Para los muros verticales, se tiene

unc = unc + θnz xc

, lo que puede escribirse como

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

donde para un marco orientado en ’x’

ki = aTxi kxi axi

y para un marco orientado en ’y’

ki = aTyi kxi ayi

Con esto, se obtuvo la matriz K, la que corresponde a

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

En cuanto a la matriz de masa, esta corresponde a

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.

Se determinaron las frecuencias y los periodos modales

 
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

Y los modos corresponden a las columnas de la siguiente matriz

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.

Figura 10: Modo 1

9
Figura 11: Modo 2

Figura 12: Modo 3

Problema 3. Respuestas sı́smicas máximas

(3.1) Análisis tiempo-historia real en 2 direcciones

La estructura de la Figura 8 está sometida al sismos de constitución del 2010:

(a) El canal L actúa en la dirección ± x


(b) El canal T actúa en la dirección ± y

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

 1414,740 −393,296 −64,936 0 0 0 0 0 0 


−393,296 1240,964 −425,764 0 0 0 0 0 0
 −64,936 −425,764 698,387 0 0 0 0 0 0 
0 0 0 1109,264 −295,732 −44,607 −913,912 236,330 34,2836 

0 0 0 −295,732 919,000 −291,804 242,061 −738,443 223,8092  N
C=


 0 0 0 −44,607 −291,804 495,159 34,851 228,195 −390,0797   m
 0 0 0 −913,912 242,061 34,851 8,30·104 −2,30·104 −3782,6572 
0 0 0 236,330 −738,443 228,195 −2,30·104 7,26·104 −2,48·104
0 0 0 34,284 223,809 −390,080 −3782,657 −2,48·104 4,07·104

Con esto, y luego de cargar el registro en Matlab, se obtuvieron las siguientes respuestas para cada piso.

Figura 13: Desplazamiento en x, piso 1

11
Figura 14: Desplazamiento en x, piso 3

Figura 15: Desplazamiento en x, piso 3

12
Figura 16: Desplazamiento en y, piso 1

Figura 17: Desplazamiento en y, piso 3

13
Figura 18: Desplazamiento en y, piso 3

Figura 19: Giro en Z, piso 1

14
Figura 20: Giro en Z, piso 2

Figura 21: Giro en Z, piso 3

(3.2) Análisis tiempo-historia real en 2 direcciones

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 creo el vector de fuerzas

p(t) = Ln [ügL ; ügT ]

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.

Figura 22: Comparación Desplazamiento en x, piso 1

16
Figura 23: Comparación Desplazamiento en x, piso 3

Figura 24: Comparación Desplazamiento en x, piso 3

17
Figura 25: Comparación Desplazamiento en y, piso 1

Figura 26: Comparación Desplazamiento en y, piso 3

18
Figura 27: Comparación Desplazamiento en y, piso 3

Figura 28: Comparación Giro en Z, piso 1

19
Figura 29: Comparación Giro en Z, piso 2

Figura 30: Comparación Giro en Z, piso 3

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

(3.3) Desplazamientos máximos de las losas

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

Figura 32: Desplazamiento máximo esquina superior derecha, eje x

Figura 33: Desplazamiento máximo esquina superior derecha, eje y

22
Esquina inferior derecha:

dxmax = 0,0909m
dymax = −0,1385m

Figura 34: Desplazamiento máximo esquina inferior derecha, eje x

Figura 35: Desplazamiento máximo esquina inferior derecha, eje y

23
Esquina superior derecha:

dxmax = 0,0976m
dymax = 0,1343m

Figura 36: Desplazamiento máximo esquina superior derecha, eje x

Figura 37: Desplazamiento máximo esquina superior derecha, eje y

24
Esquina superior izquierda:

dxmax = −0,0814m
dymax = 0,1302m

Figura 38: Desplazamiento máximo esquina superior izquierda, eje x

Figura 39: Desplazamiento máximo esquina superior izquierda, eje y

25

También podría gustarte