کد تحلیل دینامیکی سازه (عمران / سازه)

رشته : عمران گرایش سازه من کد تحلیل دینامیکی یک سازه چند درجه آزاد برشی رو تحت بار زلزله در تراز پایه میخوامکه البته کد رو دارم فقط ماتریس جرم و میرایی و سختی اش تعریف نشده آ نا رو میخوام تعریف بشه و کد نوشته شده توضیح داده بشه. function Result=MDOF\_simulation(M,C,K,f,fs)n=size(f,1);dt=1/fs; %sampling rate\[Vectors, Values\]=eig(K,M);Freq=sqrt(diag(Values))/(2\*pi); % undamped natural frequencysteps=size(f,2);Mn=diag(Vectors'\*M\*Vectors); % uncoupled massCn=diag(Vectors'\*C\*Vectors); % uncoupled dampingKn=diag(Vectors'\*K\*Vectors); % uncoupled stifnesswn=sqrt(diag(Values));zeta=Cn./(2\*sqrt(Mn.\*Kn)); % damping ratiowd=wn.\*sqrt(1-zeta.^2);fn=Vectors'\*f; % generalized input force matrixt=\[0:dt:dt\*steps-dt\];%forced vibrationfor i=1:1:n h(i,:)=(1/(Mn(i)\*wd(i))).\*exp(-zeta(i)\*wn(i)\*t).\*sin(wd(i)\*t); %transfer function of displacement hd(i,:)=(1/(Mn(i)\*wd(i))).\*(-zeta(i).\*wn(i).\*exp(-zeta(i)\*wn(i)\*t).\*sin(wd(i)\*t)+wd(i).\*exp(-zeta(i)\*wn(i)\*t).\*cos(wd(i)\*t)); %transfer function of velocity hdd(i,:)=(1/(Mn(i)\*wd(i))).\*((zeta(i).\*wn(i))^2.\*exp(-zeta(i)\*wn(i)\*t).\*sin(wd(i)\*t)-zeta(i).\*wn(i).\*wd(i).\*exp(-zeta(i)\*wn(i)\*t).\*cos(wd(i)\*t)-wd(i).\*((zeta(i).\*wn(i)).\*exp(-zeta(i)\*wn(i)\*t).\*cos(wd(i)\*t))-wd(i)^2.\*exp(-zeta(i)\*wn(i)\*t).\*sin(wd(i)\*t)); %transfer function of acceleration qq=conv(fn(i,:),h(i,:))\*dt; qqd=conv(fn(i,:),hd(i,:))\*dt; qqdd=conv(fn(i,:),hdd(i,:))\*dt; q(i,:)=qq(1:steps); % modal displacement qd(i,:)=qqd(1:steps); % modal velocity qdd(i,:)=qqdd(1:steps); % modal acceleration endx=Vectors\*q; %displacementv=Vectors\*qd; %vecloitya=Vectors\*qdd; %vecloity% Free vibrationxi=zeros(n,1); % displacement initial conditionvi=zeros(n,1); % velocity initial conditionxno=Vectors'\*M\*xi./Mn;vno=Vectors'\*M\*vi./Mn;for i=1:1:nAA=(vno(i)+xno(i).\*zeta(i).\*wn(i))./wd(i);BB=xno(i);qf(i,:)=exp(-zeta(i)\*wn(i)\*t).\*(AA.\*sin(wd(i)\*t)+BB.\*cos(wd(i)\*t));qdf(i,:)=wd(i)\*exp(-zeta(i)\*wn(i)\*t).\*(AA.\*cos(wd(i)\*t)-BB.\*sin(wd(i)\*t))-zeta(i).\*wn(i).\*exp(-zeta(i)\*wn(i)\*t).\*(AA.\*sin(wd(i)\*t)+BB.\*cos(wd(i)\*t));qddf(i,:)=wd(i)^2\*exp(-zeta(i)\*wn(i)\*t).\*(-AA.\*sin(wd(i)\*t)-BB.\*cos(wd(i)\*t))-2\*zeta(i).\*wn(i).\*wd(i).\*exp(-zeta(i)\*wn(i)\*t).\*(AA.\*cos(wd(i)\*t)-BB.\*sin(wd(i)\*t))+zeta(i)^2.\*wn(i)^2\*exp(-zeta(i)\*wn(i)\*t).\*(-AA.\*sin(wd(i)\*t)-BB.\*cos(wd(i)\*t));endx=x+Vectors\*qf;v=v+Vectors\*qdf;a=a+Vectors\*qddf;Result.Displacement=x;Result.Velocity=v;Result.Acceleration=a;Result.Parameters.Freq=Freq;Result.Parameters.DampRatio=zeta\*100;Result.Parameters.ModeShape=Vectors;end که من این کد رو دارم ولی این کد کامل نیستاولش باید ماتریس های جرم و سختی و میرایی و رکورد زلزله تعریف بشنو برام کلا توضیح داده بشه صورت پروژه:یک سازه چند درجه ازاد برشی را تحلیل دینامیکی کنید تحت بار زلزله در تراز پایه.