mainprogram.asv

matlab编写结构计算振型分解法

%                  振型分解法
%%%%%%%%%%%%%%%%%%%  MainProgram  %%%%%%%%%%%%%%%%%%%%%%%%%
%结构原始参数输入
h0=[3.0,3.0,3.0];                %单位：m
m0=[2.561 2.545 0.559]*1e+6;    %单位：kg
k0=[5.43 9.03 8.23]*1e+8;        %单位：N/m
cn=length(m0);                     %层数cn

Tg=0.4;amax=0.16;                  %求解地震影响系数所需参数,查规范

r=0.9;n1=0.02;n2=1.0;              %一般结构的阻尼参数
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

[w,zhx]=plzhx(m0,k0,cn)            %求解自振频率w,振型zhx
[q,a]=zhxcyxsh(zhx,m0,w,cn);       %求解振型参与系数q，地震影响系数a
[F]=dzhli(a,q,zhx,m0,cn);          %求解水平地震作用F
[V]=cjdzhjli(F,cn)                 %求解层间地震剪力V
[u,Thelta]=cjwy(V,k0,h0,cn)