기계진동학 텀프로젝트

목차

1. 자유물체도
2. 수식풀이
3. 결과값 분석
4. MATLAB CODE

본문내용

1. 자유물체도
물체에 작용하는 x 방향의 힘이 상쇄되므로
,

2. 수식풀이
운동방정식
해를 로 가정하면
교재 P.272, 식 (3.82) 및 (3.83)으로부터
,

가진력의 출력 ,
속도
힘은 이 계에 가해지는 힘이 y 방향의 힘 이므로

주기 에서의 출력은 다음의 식으로 주어진다.

여기서 다음 조건을 만족하는 w, e, M, m, k, c값을 구한다.
조건 1)
조건 2)
조건 3)
조건 4)

3. 결과값 분석
m=20kg, M=1000kg, k=100N/m, c=1.0kg/s, e=0.1m로 가정할 때
조건에 맞는 결과 값 w=188.49rad/sec, X=0.002, P=1.07KW 을 도출할 수 있었다.

w=188.49rad/sec
→
조건 1) 에 만족
X=0.002m
→
조건 2) 에 만족
P= 1.070W
→
조건 3) 에 만족

→
조건 4) 에 만족

4. MATLAB CODE

function mv(e,k,c,M,m)
m=20; M=1000; k=100; c=1.0; e=0.1;
fid=fopen('result.txt','w');

for w=40*pi:pi:60*pi
Wn=sqrt(k/M);
Cc=2*sqrt(M*k);
zt=c/Cc;
r=w/Wn;
st=atan((2*zt*r)/(1-r^2));
X=m*e*w^2/(sqrt((k-M*w^2)^2+(c*w)^2));
x_3=0;

for t=0:0.0001:2*pi/w;
x=X*sin(w*t-st);
x_1=X*w*cos(w*t-st);
x_2=(m*e*w^2)*sin(w*t);
x_3=x_3+abs(M*x_1*x_2);
P=(abs(-X*e*m*w^3*(cos(2*w*t-st)/(4*w)-(t*sin(st))/2)))/(2*pi/w);
fprintf(fid,'%.7d %.2d %.2d\n',w,X,P);
end
end

참고 자료

없음