[구조동역학] (matlab) 주파수영역방법 (frequency domain method)을 이용한 다자유도 지진해석

소개글

[구조동역학] (matlab) 주파수영역방법 (frequency domain method)을 이용한 지진해석

매틀랩을 이용하여 주파수 영역 방법으로 엘센트로 지진을 해석한 매트랩 코딩 자료입니다.

3층 전단빌딩의 대해 각 층별 변위, 속도, 가속도를 나타내었습니다.

컴파일 실행환경

MATLAB 7.1

본문내용

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Frquency Domain Method를 이용한 MDOF해석
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
close all
clc


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Frquency Domain Method를 이용한 응답
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% 물성치 입력
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);

% % dt = 0.02;
% % T = 10;
% % t = 0:dt:T ;
% % df = 1/T ;
% % fmax = 1/dt;
% % f = 0.1:df:fmax;
% % w = 2*pi*f;
% % c = 2*xi*w


%%% F(w)를 구하는 루틴
dt = 0.02;
T = 10;

참고 자료

없음