地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab編程

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程一、地震波生成反應(yīng)譜

圖1

***********讀入地震記錄***********

ElCentro(:,1)*9.8067;%單位統(tǒng)一為

%單位

%初始化各儲(chǔ)存向量

%相對(duì)位移

%相對(duì)速度

%絕對(duì)加速度

%阻尼比

%結(jié)構(gòu)周期

%地震記錄的步長(zhǎng)

%結(jié)構(gòu)自振頻率

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

%根據(jù)地震記錄,計(jì)算不同的反應(yīng)

Displace=zeros(1,N);%初始化各儲(chǔ)存向量，避免下次不同周期計(jì)算時(shí)引用到前一個(gè)周期的結(jié)果

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

%繪制地震記錄圖

%繪制位移反應(yīng)譜legend('ζ=0.02')

%繪制速度反應(yīng)譜legend('ζ=0.02')

%繪制絕對(duì)加速度反應(yīng)譜

legend('ζ=0.02')

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

%%規(guī)范反應(yīng)譜取值程序

%%%%烈度選擇

%%%%場(chǎng)地類別，設(shè)計(jì)地震分組選擇

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

%%%%%阻尼比

%%%%%分段位置

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

T1～T2

T3～6.0

6.0～

%%%%%%反應(yīng)譜值

%%%主程序%%%%%%%%確定需要控制的反應(yīng)譜

Sa（T）(T=T1,...,TM)的坐標(biāo)點(diǎn)數(shù)

M，

ceita=0.035;%%%阻尼比：0.035

變換的頻率差：2*pi*0.005(可以保證長(zhǎng)周期項(xiàng)

項(xiàng)三角級(jí)數(shù))；%%%%頻率變化范圍

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

%%%保證率%%%%%%人造地震動(dòng)持續(xù)時(shí)間

40s，時(shí)間間隔：0.02s

%%%%上升段

%%%%%平穩(wěn)段;

%%%%衰減段參數(shù)

5%的要求,需要時(shí)程動(dòng)力分析

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

%%%%%%%修改

%%%%%循環(huán)體函數(shù)

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

%%%%%%%修改

地震工程學(xué)-反應(yīng)譜和地震時(shí)程波的相互轉(zhuǎn)化matlab

編程

%%%%%%%修改

地震工程學(xué)-反應(yīng)