ru_u1.m

来自「钢、混凝土结构时频反应谱程序」· M 代码 · 共 62 行

%该程序为vidic主程序
clear;%清除工作空间变量；
clc;format compact;format long g;
%per=[0:0.2:1.0 1.5:0.5:2 3:1:6];%设定自振周期范围,分别采用不同的密度；
per=0:0.1:3;
sp=length(per);
uc=2;
mag=0.4*9.81;
mvg=mag/0.6/9.81;
Tc=2*pi*2.0/2.5*mvg/mag;
To=0.65*Tc*uc^0.3;
for i=1:sp
    T=per(i);
    if(T<=To)
        Ru(i)=(uc-1)*T/To+1;
    else
        Ru(i)=uc;
    end 
end
figure(1)
plot(per,Ru,'kd-.') %绘制出对应曲线图；
hold on

uc=4;
To=0.65*Tc*uc^0.3;
for i=1:sp
    T=per(i);
    if(T<=To)
        Ru(i)=(uc-1)*T/To+1;
    else
        Ru(i)=uc;
    end 
end
plot(per,Ru,'kd-.') %绘制出对应曲线图；
hold on

uc=6;
To=0.65*Tc*uc^0.3;
for i=1:sp
    T=per(i);
    if(T<=To)
        Ru(i)=(uc-1)*T/To+1;
    else
        Ru(i)=uc;
    end 
end
plot(per,Ru,'kd-.') %绘制出对应曲线图；
hold on

uc=8;
To=0.65*Tc*uc^0.3;
for i=1:sp
    T=per(i);
    if(T<=To)
        Ru(i)=(uc-1)*T/To+1;
    else
        Ru(i)=uc;
    end 
end
plot(per,Ru,'kd-.') %绘制出对应曲线图；
hold on