[구조동역학] (matlab) 모드중첩법을 이용한 다자유도의 해석

소개글

[구조동역학] (matlab) 모드중첩법을 이용한 다자유도의 해석

매틀랩을 이용하여 3층 전단빌딩에 대해 엘센트로 지진을 해석한 자료

다자유도 해석 방법중 모드 중첩법을 이용하여 해석하였음.

컴파일 실행환경

MATLAB 7.1

본문내용

% close all
clear all
figure
clc


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 모드중첩법을 이용한 MDOF 해석
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


m=[1 0 0 ; 0 1 0 ; 0 0 1];
% k=[1 -1 0 ; -1 3 -2 ; 0 -2 5]*100;
k=[5 -2 0 ; -2 3 -1 ; 0 -1 1]*100;
xi=0.1;
[vec,value]=eig(k,m)
w=sqrt(value);
% c=2*xi*w;
c=inv((vec)`) * (vec`*m*vec) * 2 * xi * w * inv(vec);


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%% 단자유도로 해석하기 위한 m,c,k,p 변환 --- 1 ----

modalM=vec`*m*vec;
modalC=2*modalM(1)*0.1*w;
modalK=vec`*k*vec;
M1=modalM(1);
C1=modalC(1);
K1=modalK(1);

q0=0;
qdot0=0;

load elcen_ns.dat;
t=elcen_ns(:,1);
ga=elcen_ns(:,2);
pt=-M1*elcen_ns(:,2);
P=[pt pt pt] ;
Pt=(vec`*P`)` ;
R0 = Pt(:,1) ;

참고 자료

없음