response_n.m

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

function [u,rd,ad,rv,av,ra,aa,rep]=response_n(para,sa,per,xg,dt,at,eb)
t=0:dt:at;
st=length(t);
sp=length(per);
disp=zeros(st,sp);
rd=zeros(st,sp);
ad=zeros(st,sp);
rv=zeros(st,sp);
av=zeros(st,sp);
ra=zeros(st,sp);
aa=zeros(st,sp);
u=zeros(st,sp);
for np=1:sp
	if (np==1)
		rd(:,np)=0;
		ad(:,np)=0;
		rv(:,np)=0;
		av(:,np)=0;
		ra(:,np)=-xg(1);
		aa(:,np)=0;
        u(:,np)=0;
        perin=0;
        rep(1)=0;%此处处理有问题，因为基频越高则耗能越大；
	else
        perin=per(np);
        ssa=sa(np);
		[yu,d,v,a,ep]=nonlinear(para,ssa,np,perin,xg,t,dt,eb); 
		dis=abs(d);
		vel=abs(v);
		acc=abs(a+xg);
        rep(np)=ep;
		for nt=1:st
			rd(nt,np)=d(nt);
			ad(nt,np)=dis(nt);
			rv(nt,np)=v(nt);
			av(nt,np)=vel(nt);
			ra(nt,np)=a(nt);
			aa(nt,np)=acc(nt);
            u(nt,np)=yu(nt);
		end
	end
end