结构动力分析频域分析法模态叠加法.docVIP

结构动力分析频域分析法——模态叠加法 说明： %a:频域内的模态分析法基本原理可以参考相关文献 %b:程序输入的参数，在前几行也基本上数目，质量矩阵和刚度矩阵全部来自于ANSYS有限元模型 %c:本程序目前只是用beam44,beam18x基本单元构成的有限元模型动力分析计算 %d:本程序目前得到脉动响应结果偏小，原因未找到 %e:希望大家找出其中的原因，一起来讨论 clc;clear;flag=input((SSRS:输入0；CQC：输入1):);fs=input(样本采样频率：);dir_in=uigetfile({*.txt;*.dat},输入风荷载文件名);dof_in=uigetfile({*.txt;*.dat},节点约束矩阵(Size:7*..或..*7));mass_in=uigetfile({*.txt;*.dat},质量矩阵);stiff_in=uigetfile({*.txt;*.dat},刚度矩阵);nodeL_in=uigetfile({*.txt;*.dat},加载节点编号);node_degree=load(dof_in);[row,col]=size(node_degree);if row==7????node_degree=node_degree;%转置成node_sum*7的矩阵????node_sum=col;else????node_sum=row;endmass_array=load(mass_in);stifness_array=load(stiff_in);DOF_sum=sum(sum(node_degree(:,2:7)0));%计算总的自由度mass_array=reshape(mass_array,DOF_sum,DOF_sum);stifness_array=reshape(stifness_array,DOF_sum,DOF_sum);%-------------------------------------------------------------------------%----------------------导入荷载时程----------------------------------------node_load=load(nodeL_in);load_L=length(node_load);nstep=1000;%与节点位移方向对应load_dof=zeros(6,node_sum,nstep);fid=fopen(dir_in,r);for ii=1:nstep????data_in=fscanf(fid,%g,[3 load_L]);????load_dof(1:3,node_load,ii)=data_in;endfclose(fid);%-------------------------------------------------------------------------count=0;load_F=zeros(DOF_sum,nstep);mean_F=zeros(DOF_sum,1);std_F=zeros(DOF_sum,1);for ii=1:node_sum????for jj=2:7????????if node_degree(ii,jj)0.5?????????count=count+1;?????????if jj5??%只读XYZ方向的力?????????load_F(count,:)=load_dof(jj-1,ii,:);?????????mean_F(count)=mean(load_F(count,:));?????????std_F(count)=std(load_F(count,:));?????????load_F(count,:)=load_F(count,:)-mean_F(count);?????????end????????end?????????endend%-------------------------------------------------------------------------%静风响应和背景响应y_mean=stifness_array\mean_F;y_std=stifness_array\std_F;out=[y_mean y_std];save 静风背景响应.txt out -ascii%-----------------------生成广义荷载谱-------------------------------------%步A:自振频率