기계진동학 텀프로젝트
- 최초 등록일
- 2020.09.14
- 최종 저작일
- 2020.09
- 8페이지/ 한컴오피스
- 가격 5,000원
* 본 문서(hwp)가 작성된 한글 프로그램 버전보다 낮은 한글 프로그램에서 열람할 경우 문서가 올바르게 표시되지 않을 수 있습니다.
이 경우에는 최신패치가 되어 있는 2010 이상 버전이나 한글뷰어에서 확인해 주시기 바랍니다.
목차
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
참고 자료
없음