地大 Matlab机械原理第二次作业2好学长不留名.docx
- 文档编号:25556784
- 上传时间:2023-06-09
- 格式:DOCX
- 页数:20
- 大小:193.49KB
地大 Matlab机械原理第二次作业2好学长不留名.docx
《地大 Matlab机械原理第二次作业2好学长不留名.docx》由会员分享,可在线阅读,更多相关《地大 Matlab机械原理第二次作业2好学长不留名.docx(20页珍藏版)》请在冰豆网上搜索。
地大Matlab机械原理第二次作业2好学长不留名
(关注文库shlimin2012上传文档,包含几乎所有课程报告)
在图10示的正弦机构中,已知lAB=100mm,h1=120mm,h2=80mm,W1=10rad/s(常数),滑块2和构件3的重量分别为G2=40N和G3=100N,质心S2和S3的位置如图所示,加于构件3上的生产阻力Fr=400N,构件1的重力和惯性力略去不计。
试用解析法求机构在Φ1=60°、150°、220°位置时各运动副反力和需加于构件1上的平衡力偶Mb。
图10机构原理示意图
解:
分别对三个构件进行受力分析如图11:
(A)构件3受力图
(B)构件2受力图
(C)构件1受力图
图11构件受力图
(1)滑块2:
(1)
构件3:
(2)
(2)确定惯性力:
(3)
(3)各构件的平衡方程:
构件3:
(4)
构件2:
(5)
构件1:
(6)
总共有八个方程,八个未知数。
归纳出一元八次方程矩阵:
=
(7)
由AX=B进而可得:
X=A\B
接下来运用软件工具进行分析。
1、进行MATLAB编程分析
首先编写函数F用于实现上述运算功能:
functiony=F(x)
%
%inputparameters
%
%x
(1)=lAB
%x
(2)=h1
%x(3)=h2
%x(4)=W1
%x(5)=G2
%x(6)=G3
%x(7)=Fr
%x(8)=theta1
%
%outputparameters
%
%y
(1)=FR23
%y
(2)=FR4'
%y(3)=FR4
%y(4)=FR12x
%y(5)=FR12y
%y(6)=FR41x
%y(7)=FR41y
%y(8)=Mb
%
A=[10000000;
01-100000;
-x
(1)*cos(x(8))/x(3)0100000;
00010000;
-10001000;
000-10100;
0000-1010;
-x
(1)*cos(x(8))0000001];
B=[x(7)-(x(6)/10)*x
(1)*x(4)^2*sin(x(8));0;0;(x(5)/10)*x
(1)*x(4)^2*cos(x(8));-(x(5)/10)*x
(1)*x(4)^2*sin(x(8));0;0;0];
y=A\B;
接下来运行程序计算Φ1=60°的各未知量值:
lAB=0.1;
h1=0.120;
h2=0.08;
W1=10;
G2=40;
G3=100;
Fr=400;
th1=60*pi/180;
x=[lABh1h2W1G2G3Frth1];
y=F(x)
y=
313.3975
195.8734
195.8734
20.0000
278.7564
20.0000
278.7564
15.6699
得到:
Φ1=60°时FR23=FR32=313.3975N;FR4=FR4’=195.8734N;FR12x=20.0000N;FR12y=278.7564N;FR41x=20.0000N;FR41y=278.7564N;Mb=15.6699N*m。
运行程序计算Φ1=150°的各未知量值:
>>th1=150*pi/180;
>>x=[lABh1h2W1G2G3Frth1];
>>y=F(x)
y=
350.0000
-378.8861
-378.8861
-34.6410
330.0000
-34.6410
330.0000
-30.3109
得到:
Φ1=150°时FR23=FR32=350.0000N;FR4=FR4’=-378.8861N;FR12x=-34.6410N;FR12y=330.0000N;FR41x=-34.6410N;FR41y=330.0000N;Mb=-30.3109N*m。
运行程序计算Φ1=220°的各未知量值:
>>th1=220*pi/180;
>>x=[lABh1h2W1G2G3Frth1];
>>y=F(x)
y=
464.2788
-444.5727
-444.5727
-30.6418
489.9903
-30.6418
489.9903
-35.5658
得到:
Φ1=220°时FR23=FR32=464.2788N;FR4=FR4’=-444.5727N;FR12x=-30.6418N;FR12y=489.9903N;FR41x=-30.6418N;FR41y=489.9903N;Mb=-35.5658N*m。
接下来取Φ1=0~360°范围分析其受力:
h1=0.120;
h2=0.08;
W1=10;
G2=40;
G3=100;
Fr=400;
th1=linspace(0,2*pi,36);
x=zeros(length(th1),8);
forn=1:
36
x(n,:
)=[lABh1h2W1G2G3Frth1(n)];
end
p=zeros(8,length(th1));
fork=1:
36
p(:
k)=F(x(k,:
));
end
>>p
p=
Columns1through8
400.0000382.1443364.8625348.7101334.2061321.8169311.9404304.8943
500.0000470.0039426.9963374.1872314.6014250.8119184.7735117.7719
500.0000470.0039426.9963374.1872314.6014250.8119184.7735117.7719
40.000039.357237.449434.338030.122924.939618.954712.3607
400.0000375.0020350.8075328.1941307.8886290.5436276.7166266.8521
40.000039.357237.449434.338030.122924.939618.954712.3607
400.0000375.0020350.8075328.1941307.8886290.5436276.7166266.8521
40.000037.600334.159729.935025.168120.065014.78199.4218
Columns9through16
300.9050300.1007302.5072308.0472316.5427327.7205341.2215356.6116
50.4893-16.8300-84.1427-151.3378-217.9780-283.0943-345.0675-401.6200
50.4893-16.8300-84.1427-151.3378-217.9780-283.0943-345.0675-401.6200
5.3693-1.7946-8.9008-15.7210-22.0359-27.6425-32.3607-36.0388
261.2670260.1410263.5101271.2661283.1597298.8087317.7101339.2563
5.3693-1.7946-8.9008-15.7210-22.0359-27.6425-32.3607-36.0388
261.2670260.1410263.5101271.2661283.1597298.8087317.7101339.2563
4.0391-1.3464-6.7314-12.1070-17.4382-22.6475-27.6054-32.1296
Columns17through24
373.3963391.0361408.9639426.6037443.3884458.7785472.2795483.4573
-449.9252-486.8273-509.1470-514.0376-499.3489-463.9495-407.9684-332.9190
-449.9252-486.8273-509.1470-514.0376-499.3489-463.9495-407.9684-332.9190
-38.5585-39.8390-39.8390-38.5585-36.0388-32.3607-27.6425-22.0359
362.7548387.4505412.5495437.2452460.7437482.2899501.1913516.8403
-38.5585-39.8390-39.8390-38.5585-36.0388-32.3607-27.6425-22.0359
362.7548387.4505412.5495437.2452460.7437482.2899501.1913516.8403
-35.9940-38.9462-40.7318-41.1230-39.9479-37.1160-32.6375-26.6335
Columns25through32
491.9528497.4928499.8993499.0950495.1057488.0596478.1831465.7939
-241.6872-138.3782-28.034983.7439191.2451289.0952372.6779438.4701
-241.6872-138.3782-28.034983.7439191.2451289.0952372.6779438.4701
-15.7210-8.9008-1.79465.369312.360718.954724.939630.1229
528.7339536.4899539.8590538.7330533.1479523.2834509.4564492.1114
-15.7210-8.9008-1.79465.369312.360718.954724.939630.1229
528.7339536.4899539.8590538.7330533.1479523.2834509.4564492.1114
-19.3350-11.0703-2.24286.699515.299623.127629.814235.0776
Columns33through36
451.2899435.1375417.8557400.0000
484.2616509.2386513.9257500.0000
484.2616509.2386513.9257500.0000
34.338037.449439.357240.0000
471.8059449.1925424.9980400.0000
34.338037.449439.357240.0000
471.8059449.1925424.9980400.0000
38.740940.739141.114140.0000
P矩阵的每一行分别是八个未知量在Φ1等于36个分量下的值。
整理表格如下:
表一
Φ1
FR23/FR32
FR4
FR4’
FR12X
FR12Y
FR41X
FR41Y
Mb
rad
N
N
N
N
N
N
N
N*m
0
0.1795
0.3590
0.5386
0.7181
0.8976
1.0771
1.2566
1.4362
1.6157
1.7952
1.9747
2.1542
2.3338
2.5133
2.6928
2.8723
3.0518
3.2314
3.4109
3.5904
3.7699
3.9494
4.1290
4.3085
4.4880
4.6675
4.8470
5.0265
5.2061
5.3856
5.5651
5.7446
5.9241
6.1037
6.2832
400.0000
382.1443
364.8625
348.7101
334.2061
321.8169
311.9404
304.8943
300.9050
300.1007
302.5072
308.0472
316.5427
327.7205
341.2215
356.6116
373.3963
391.0361
408.9639
426.6037
443.3884
458.7785
472.2795
483.4573
491.9528
497.4928
499.8993
499.0950
495.1057
488.0596
478.1831
465.7939
451.2899
435.1375
417.8557
400.0000
500.0000
470.0039
426.9963
374.1872
314.6014
250.8119
184.7735
117.7719
50.4893
-16.8300
-84.1427
-151.3378
-217.9780
-283.0943
-345.0675
-401.6200
-449.9252
-486.8273
-509.1470
-514.0376
-499.3489
-463.9495
-407.9684
-332.9190
-241.6872
-138.3782
-28.0349
83.7439
191.2451
289.0952
372.6779
438.4701
484.2616
509.2386
513.9257
500.0000
500.0000
470.0039
426.9963
374.1872
314.6014
250.8119
184.7735
117.7719
50.4893
-16.8300
-84.1427
-151.3378
-217.9780
-283.0943
-345.0675
-401.6200
-449.9252
-486.8273
-509.1470
-514.0376
-499.3489
-463.9495
-407.9684
-332.9190
-241.6872
-138.3782
-28.0349
83.7439
191.2451
289.0952
372.6779
438.4701
484.2616
509.2386
513.9257
500.0000
40.0000
39.3572
37.4494
34.3380
30.1229
24.9396
18.9547
12.3607
5.3693
-1.7946
-8.9008
-15.7210
-22.0359
-27.6425
-32.3607
-36.0388
-38.5585
-39.8390
-39.8390
-38.5585
-36.0388
-32.3607
-27.6425
-22.0359
-15.7210
-8.9008
-1.7946
5.3693
12.3607
18.9547
24.9396
30.1229
34.3380
37.4494
39.3572
40.0000
400.0000
375.0020
350.8075
328.1941
307.8886
290.5436
276.7166
266.8521
261.2670
260.1410
263.5101
271.2661
283.1597
298.8087
317.7101
339.2563
362.7548
387.4505
412.5495
437.2452
460.7437
482.2899
501.1913
516.8403
528.7339
536.4899
539.8590
538.7330
533.1479
523.2834
509.4564
492.1114
471.8059
449.1925
424.9980
400.0000
40.0000
39.3572
37.4494
34.3380
30.1229
24.9396
18.9547
12.3607
5.3693
-1.7946
-8.9008
-15.7210
-22.0359
-27.6425
-32.3607
-36.0388
-38.5585
-39.8390
-39.8390
-38.5585
-36.0388
-32.3607
-27.6425
-22.0359
-15.7210
-8.9008
-1.7946
5.3693
12.3607
18.9547
24.9396
30.1229
34.3380
37.4494
39.3572
40.0000
400.0000
375.0020
350.8075
328.1941
307.8886
290.5436
276.7166
266.8521
261.2670
260.1410
263.5101
271.2661
283.1597
298.8087
317.7101
339.2563
362.7548
387.4505
412.5495
437.2452
460.7437
482.2899
501.1913
516.8403
528.7339
536.4899
539.8590
538.7330
533.1479
523.2834
509.4564
492.1114
471.8059
449.1925
424.9980
400.0000
40.0000
37.6003
34.1597
29.9350
25.1681
20.0650
14.7819
9.4218
4.0391
-1.3464
-6.7314
-12.1070
-17.4382
-22.6475
-27.6054
-32.1296
-35.9940
-38.9462
-40.7318
-41.1230
-39.9479
-37.1160
-32.6375
-26.6335
-19.3350
-11.0703
-2.2428
6.6995
15.2996
23.1276
29.8142
35.0776
38.7409
40.7391
41.1141
40.0000
按照上述分析结果输出图像:
输出FR23和FR32的图像,程序如下:
plot(th1,p(1,:
),'--')
ylabel('FR23,FR32/N');
xlabel('\theta1/rad');
图12FR23和FR32图像
输出FR4和FR4’的图像,程序如下:
plot(th1,p(3,:
),'--')
xlabel('\theta1/rad');
ylabel('FR4,FR4''/N');
图13FR4和FR4’图像
输出FR12x的图像,程序如下:
plot(th1,p(4,:
),'--')
ylabel('FR12x/N');
xlabel('\theta1/rad');
图14FR12x的图像
输出FR12y的图像,程序如下:
plot(th1,p(5,:
),'--')
ylabel('FR12y/N');
xlabel('\theta1/rad');
图15FR12y的图像
输出FR41X的图像,程序如下:
plot(th1,p(6,:
),'--')
xlabel('\theta1/rad');
ylabel('FR41x/N');
图16FR41X变化图像
输出FR41Y的图像,程序如下:
plot(th1,p(7,:
),'--')
xlabel('\theta1/rad');
ylabel('FR41y/N');
图17FR41Y变化图像
输出Mb的图像,程序如下:
plot(th1,p(8,:
),'--')
xlabel('\theta1/rad');
ylabel('Mb/N*m');
图18Mb变化图像
参考文献
【1】孙恒,陈作模.机械原理【M】.7版.北京:
高等教育出版社,2006
【2】曲秀全.基于MATLAB/Simulink平面连杆机构的动态仿真.黑龙江:
哈尔滨工业大学出版社.2007
【3】贺超英.MATLAB应用与实验教程.北京:
电子工业出版社,2010
【4】王宏.MATLAB6.5及其在信号处理中的应用.北京:
清华大学出版社,2004
【5】同济大学数学系.线性代数.北京:
高等教育出版社,2009
【6】同济大学数学系.高等数学.北京:
高等教育出版社,2008
【7】雷培.刘云霞.基于MATLAB的四连杆机构运动分析.机械工程与自动化.2009年4月,第2期
【8】张建民,白景岭.理论力学.武汉:
中国地质大学出版社,2006
【9】叶修梓,陈超强.SolidWorks高级装配教程.北京:
机械工业出版社
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地大 Matlab机械原理第二次作业2好学长不留名 Matlab 机械 原理 第二次 作业 好学 留名