基于matlab的平面四连杆机构
设计
领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计
以及该机构的运动分析
基于matlab的平面四连杆机构设计以及该机构的运动仿真分析
摘要
四连杆机构因其结构方便灵活,能够传递动力并实现多种运动形式而被广泛应用于各个领域,因此对其进行运动分析具有重要的意义。传统的分析方法主要应用几何综合法和解析综合法,几何综合法简单直观,但是精确度较低;解析法精确度较高,但是计算工作量大。随着计算机辅助数值解法的发展,特别是MATLAB软件的引入,解析法已经得到了广泛的应用。对于四连杆的运动分析,若应用MATLAB则需要大量的编程,因此我们引入proe软件,我们不仅可以在此软件中建立实物图,而且还可以对其进行运动仿真并对其运动分析。
在设计四连杆时,我们利用解析综合法建立数学模型,再根据数学模型在MATLAB中编程可以求得其他杆件的长度。
针对范例中所求得的各连杆的长度,我们在proe软件中画出其三维图(如图4)并在proe软件中进行仿真分析得出的角加速度B,C的变化,从而得到两接触处所受到的力是成周期性变化的,可以B,C
看出两点处极易疲劳断裂,针对两点处的疲劳断裂,我们提B,CB,C
出了在设计四连杆中的一些建议。
关键字:解析法 MATLAB软件 proe软件 运动仿真
建立用解析法设计平面四杆机构模型
CDAB对于问题中所给出的连架杆的三个位置与连架杆的三个位
,,,,,,,,,,,112233置相对应,即三组对应位置为:,其中他们对应的值
,,,,,,135,112,90,82,45,52分别为:,为了便于写代数式,可作出AB与CD对应的关系,其图如下:
图—2 AB与CD三个位置对应的关系
通过上图我们可以通过建立平面直角坐标系并利用解析法来求解,其直角坐标系图如下:
图—3 平面机构直角坐标系
,,OXY00通过建立直角坐标系,如上图所示,其中与为AB杆与
θi
φi
iα
φα
a,b,c,dCD杆的初始角,各杆件的长度分别用矢量,表示,将各矢量分别在X轴与Y轴上投影的方程为
,,,*cos(),*cos(),,*cos()abdc,
,*sin(),*sin(),*sin()a,b,c,,
,在上述的方程中我们可以消除,从而可以得到与之间的关系,,如下:
2222 (1) (a,c,d,b),2cdcos(,),2accos(,,,),2abcos(,)
为便于化简以及matlab编程我们可以令:
2222,a,c,d,bH,1,ac2,
d,H, (2) ,2a,
d,H,3,c,
通过将(2)式代入(1)式中则可以化简得到如下等式:
(3) H,Hcos(,),cos(,,,),Hcos(,)123
我们可以通过(3)式将两连架杆对应的位置带入(3)式中,我们可以得到如下方程:
,,,,H,Hcos(),cos(,),Hcos(),1211131,,,,,H,Hcos(),cos(,),Hcos() (4) ,1222232
,,cos(),cos(,),cos()HH,,,H,1233333,
H,H,H联立(4)方程组我们可以求得,再根据(2)中的条件123
d以及所给定的机架的长度,我们可以求出其它杆件的长度为:
d,,a,H2,
d,,c (5) ,H3,
,222b,a,c,d,2acH1,,
四连杆设计范例:
在日常生活中,我们经常看到消防门总能自动关上,其实它是利用四连杆机构与弹簧组成的。已知门与另一连架杆三组对应的位置值
,,,,,,分别为:。1试设计一四连杆机构用于消防门上,135,112,90,82,45,52
2并对你所设计的四连杆机构提出一些建议。
范例的求解:
H,H,H,在matlab中,我们首先构建用于求的矩阵,并调用求123
H,H,H,可逆函数inv来求解线性方程组,即可得到的值,然后可以123
d根据机架的长度,可以求出a,b,c的长度(程序见附录一),求解的结果为:
H,0.7384,1,H,1.2162求解的系数为: ,2
,H,1.80973,
a,27.6293,
,b,57.2363,各杆的长度为: ,c,41.1104,
,d,50.0000,
范例的分析:
根据四连杆的判别方法可知:对上述设计所得到四连杆长度可知,是最短的并且最短杆与最长杆要小于其他两杆之和即;a
a,b,c,dd,故上述所设计得到的四连杆机构中存在曲柄,由于杆为机架且与曲柄相连,故可以进一步判定该四连杆机构为曲柄摇杆机a
构。
范例中四连杆机构的运动仿真分析:
针对问题二我们利用proe软件对已经设计好的四连杆进行仿真运动分析,则其在proe中的实物图如下:
图—4 四连杆实物图
b根据问题一中我们已经设计好的四连杆可知杆曲柄,杆为连a
20,/s杆,为了便于仿真分析,我们假设曲柄的角速度为,则可以a对图—3中的点进行角加速度分析,则可以利用proe软件的运动B,C
仿真分析可知:
B图—5 点处角加速度的变化
图—6 c点处角加速度的变化
通过对图5,6出角加速度的变化可知,,c两点的角加速度变B
化较大且呈现周期性变化,故所受到的力也是呈现周期性变化,因此在两点处极易疲劳断裂。 B,C
对于问题二,由于在两点处极易疲劳断裂,故我们可以从连B,C
杆本身以及工作环境入手,因此我们可以提出以下几条建议:
一 在连杆结构设计方面应尽量避免尖角,缺口和截面突变,以避免应力集中及由此引起的疲劳断裂。
二 降低零件表面的粗糙度,尤其是两点接触处的粗糙度,B,C
提高表面的质量,以及尽量减少能成为疲劳源的表面缺陷与表面B,C
伤损。
三 可以采用各种表面强化处理,尤其是两点接触处的表面B,C
硬度,对于两处应采用表面淬火以增强其表面硬度同时也可以保B,C
持其内部韧性,还可获得有益的表层残余应力,以抵消或降低产生疲劳裂纹的拉应力。
四 对于四连杆所受到的力是呈周期性变化,故不宜用于高速,重载的场合。
参考文献
[1] 杨柯桢.程光蕴.李仲生 机械设计基础第5版[M].北京:高等教育出版社.2009
[2] 郭仁生.机械
工程
路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理
设计分析和matlab应用第3版[M].北京:机械工业出版社.2012
[3] 王运炎.朱莉.机械工程材料第3版[M].北京:机械工业出版社.2010
[4] 李南南.吴清.曹辉林MATLAB简明教程[M].北京:清华大学出版社.2006
[5] 祝凌云.李斌.pro/ENGINEER运动仿真和有限元素分析[M].北京:人民邮电出版社.2004
附录一:
[f]=[45 90 135]*pi/180;
[p]=[52 82 112]*pi/180;
a1=[1 -cos(f(1)) cos(p(1))];
a2=[1 -cos(f(2)) cos(p(2))];
a3=[1 -cos(f(3)) cos(p(3))];
a=[a1;a2;a3];
b1=[cos(f(1)-p(1))];
b2=[cos(f(2)-p(2))];
b3=[cos(f(3)-p(3))];
b=[b1 b2 b3]';
r=inv(a)*b
d=50;
a=d/r(3)
c=d/r(2)
b=sqrt(a^2+c^2+d^2-2*a*c*r(1))