四连杆机运动学分析

栏杆机四杆机构运动学分析1四杆机构运动学分析1.1机构运动分析的任务、目的和方法曲柄摇杆机构是平面连杆机构中最基本的由转动副组成的四杆机构，它可以用来实现转动和摆动之间运动形式的转换或传递动力。对四杆机构进行运动分析的意义是：在机构尺寸参数已知的情况下，假定主动件（曲柄）做匀速转动，撇开力的作用，仅从运动几何关系上分析从动件（连杆、摇杆）的角位移、角速度、角加速度等运动参数的变化情况。还可以根据机构闭环矢量方程计算从动件的位移偏差。上述这些内容，无论是设计新的机械，还是为了了解现有机械的运动性能，都是十分必要的，而且它还是研究机械运动性能和动力性能提供必要的依据。机构运动分析的方法很多，主要有图解法和解析法。当需要简捷直观地了解机构的某个或某几个位置的运动特性时，采用图解法比较方便，而且精度也能满足实际问题的要求。而当需要精确地知道或要了解机构在整个运动循环过程中的运动特性时，采用解析法并借助计算机，不仅可获得很高的计算精度及一系列位置的分析结果，并能绘制机构相应的运动线图，同时还可以把机构分析和机构综合问题联系起来，以便于机构的优化设计。1.2机构的工作原理在平面四杆机构中，其具有曲柄的条件为：a.各杆的长度应满足杆长条件，即：最短杆长度+最长杆长度≤其余两杆长度之和。b.组成该周转副的两杆中必有一杆为最短杆，且其最短杆为连架杆或机架（当最短杆为连架杆时，四杆机构为曲柄摇杆机构；当最短杆为机架时，则为双曲柄机构）。三台设备测绘数据分别如下：第一组(2代一套)四杆机构L1=125.36mm，L2=73.4mm，L3=103.4mm，L4=103.52mm最短杆长度+最长杆长度(125.36+73.4)其余两杆长度之和(103.4+103.52)最短杆为连架杆，四杆机构为曲柄摇杆机构图1-1II-1型栏杆机机构测绘及其运动位置图第二组(2代二套)四杆机构L1=125.36mm，L2=50.1mm，L3=109.8mm，L4=72.85mm最短杆长度+最长杆长度(125.36+50.1)其余两杆长度之和(109.8+72.85)最短杆为连架杆，四杆机构为曲柄摇杆机构图1-2II-2型栏杆机机构测绘及其运动位置图第三组(3代)四杆机构L1=163.2mm，L2=64.25mm，L3=150mm，L4=90.1mm最短杆长度+最长杆长度(163.2+64.25)其余两杆长度之和(150+90.1)最短杆为连架杆，四杆机构为曲柄摇杆机构图1-3III型栏杆机机构测绘及其运动位置图在如下图1所示的曲柄摇杆机构中，构件AB为曲柄，则B点应能通过曲柄与连杆两次共线的位置。曲柄摇杆机构死点情况分析:在曲柄摇杆机构中,一般两连架杆一为主动件,一为从动件,我们知道,当从动件连架杆与连杆处于共线(拉直共线或重叠共线)位置时,机构的传动角为0,即机构处于死点位置,机构在死点位置上无法启动且具有运动不确定性,因而我们有必要对其进行详细研究.摇杆主动时曲柄摇杆机构有两个死点位置,而对曲柄主动时,有否死点位置的问题,基本没有涉及.有的资料上则直接说,曲柄主动时无死点位置.本文对此问题进行了分析研究,发现:曲柄主动时,最短杆长度+最长杆长度其余两杆长度之和,此时无死点位置.图1-4曲柄摇杆机构表1曲柄摇杆机构的死点个数及位置情况表1.3机构的数学模型的建立图1-5曲柄摇杆机构数学模型简图1.3.1建立机构的闭环矢量位置方程在用矢量法建立机构的位置方程时，需将构件用矢量来表示，并作出机构的封闭矢量多边形。如图1所示，先建立一直角坐标系。设各构件的长度分别为L1、L2、L3、L4，其方位角为、、、。以各杆矢量组成一个封闭矢量多边形，即ABCDA。其个矢量之和必等于零。即：式1式1为图1所示四杆机构的封闭矢量位置方程式。对于一个特定的四杆机构，其各构件的长度和原动件2的运动规律，即为已知，而=0，故由此矢量方程可求得未知方位角、。角位移方程的分量形式为：式2闭环矢量方程分量形式对时间求一阶导数（角速度方程）为：式3其矩阵形式为：2222224344334433cossincoscossinsinLLLLLL式4联立式3两公式可求得：式5式6闭环矢量方程分量形式对时间求二阶导数（角加速度方程）矩阵形式为：424432332222222232334244sinsinsincoscoscos434cos43cos34sin43sin3LLLLLLLLLL式7由式7可求得加速度：式8式9注：式1~式9中，Li(i=1,2,3,4）分别表示机架1、曲柄2、连杆3、摇杆4的长度；（i=1,2,3,4）是各杆与x轴的正向夹角，逆时针为正，顺时针为负，单位为rad;是各杆的角速度，，单位为rad/s;为各杆的角加速度，单位为。1.3.2求解方法（1）求导中应用了下列公式：式10（2）在角位移方程分量形式（式2）中，由于假定机架为参考系，矢量1与x轴重合，=0，则有非线性超越方程组：式11可以借助牛顿-辛普森数值解法或Matlab自带的fsolve函数求出连杆3的角位移和摇杆4的角位移。（3）求解具有n个未知量（i=1,2,…,n）的线性方程组：式12式中，系列矩阵是一个阶方阵：式13的逆矩阵为;常数项b是一个n维矢量：式14因此，线性方程组解的矢量为：式15式11是求解连杆3和摇杆4角速度和角加速度的依据。1.4基于MATLAB程序设计四连杆机构的解析法同样可以用MATLAB的计算工具来求值，并结合MATLAB的可视化手段，把各点的计算值拟合成曲线，得到四连杆机构的运动仿真轨迹。1.4.1程序流程图图1-6Matlab运动分析程序流程1.4.2M文件编写首先创建函数FoutBarPosition，函数fsolve通过他确定。functiont=fourbarposition(th)%求解其他两杆的θ_3,θ_4L1=163.2mm;L2=64.25mm;L3=150mm;L4=90.1mm;%给定已知量，各杆长L1,L2,L3,L4globalth21%给定初始θ_2t(1)=L2*cos(th21)+L3*cos(th(1))-L4*cos(th(2))-L1;t(2)=L2*sin(th21)+L3*sin(th(1))-L4*sin(th(2));主程序如下：disp'******平面四杆机构的运动分析******'L1=0.1632;L2=0.06425;L3=0.150;L4=0.091;%各杆长度globalth21%θ_2h=30;%给出转角步长30度th2=[0:h:360]*pi/180;%曲柄输入角度从0至360度，步长为pi/6th34=zeros(length(th2),2);%建立一个N行2列的零矩阵，第一列存放options=optimset('display','off');%θ_3,第二列存放θ_4form=1:length(th2)%建立for循环，求解θ_3，θ_4th21=th2(m);y3=fsolve('fourbarposition',[11]);%的非线性超越方程，结果保存在th34中th34(m,:)=y3;endy=L2*sin(th2)+L3*sin(th34(:,1)');%连杆3的C端点Y坐标值x=L2*cos(th2)+L3*cos(th34(:,1)');%连杆3的C端点X坐标值xx=[L2*cos(th2)];%连杆3的B端点X坐标值yy=[L2*sin(th2)];%连杆3的B端点Y坐标值figure(1)plot([x;xx],[y;yy],'k',[0L1],[00],'k--^',x,y,'ko',xx,yy,'ks')%绘制连杆3的几个位置点title('连杆3的几个位置点')xlabel('水平方向(m)')ylabel('垂直方向(m)')axisequal%XY坐标均衡h=5;%重新细分曲柄输入角度θ_2，步长为5度th2=[0:h:360]*pi/180;th34=zeros(length(th2),2);options=optimset('display','off');form=1:length(th2)%建立for循环，求解θ_3，θ_4th21=th2(m);y3=fsolve('fourbarposition',[11]);th34(m,:)=y3;endfigure(2)%plot(th2*180/pi,th34(:,1),th2*180/pi,th34(:,2))%绘制连杆3的角位移关于曲柄2的角位移图plot(th2*180/pi,th34(:,1)*180/pi,th2*180/pi,th34(:,2)*180/pi)%绘制摇杆4的角位移关于曲柄2的角位移图axis([03600170])%确定XY边界值grid%图形加网格xlabel('主动件转角\theta_2(度)')%横坐标标题ylabel('从动件角位移(度)')%纵坐标标题title('角位移线图')text(120,120,'摇杆4角位移')%设定显示范围text(150,40,'连杆3角位移')w2=5.23599;%设定曲柄角速度fori=1:length(th2)A=[-L3*sin(th34(i,1)),L4*sin(th34(i,2));L3*cos(th34(i,1)),-L4*cos(th34(i,2))];B=[w2*L2*sin(th2(i));-w2*L2*cos(th2(i))];w=inv(A)*B;w3(i)=w(1);%求解杆3角速度w4(i)=w(2);%求解杆4角速度%w3(i)=(w2*L2*sin(th34(i,2)-th2(i)))/(L3*sin(th34(i,1)-th34(i,2)));%w4(i)=(w2*L2*sin(th2(i)-th34(i,1)))/(L4*sin(th34(i,2)-th34(i,1)));endfigure(3)plot(th2*180/pi,w3,th2*180/pi,w4);%绘制角速度线图axis([0360-55])text(50,4.2,'摇杆4角速度(\omega_4)')text(220,3,'连杆3角速度(\omega_3)')gridxlabel('主动件转角\theta_2(度)')ylabel('从动件角速度(rad\cdots^{-1})')title('角速度线图')fori=1:length(th2)C=[L3*sin(th34(i,1)),-L4*sin(th34(i,2));L3*cos(th34(i,1)),-L4*cos(th34(i,2))];D=[w4(i)^2*L4*cos(th34(i,2))-w3(i)^2*L3*cos(th34(i,1))-w2^2*L2*cos(th2(i));w2^2*L2*sin(th2(i))+w3(i)^2*L3*sin(th34(i,1))-w4(i)^2*L4*sin(th34(i,2))];a=inv(C)*D;%s43=th34(:,2)-th34(:,1);%s23=th2'-th34(:,1);%a3=(-L2*w2^2.*cos(s23)-w3.^2*L3.*cos(-s43)+L4.*w4.^2)./(L4.*sin(s43));%a4=(L2*w2^2.*cos(s23)-w4.^2*L4.*cos(s43)+L3.*w3.^2)./(L4.*sin(s43));a3(i)=a(1);%求解杆3角加速度a4(i)=a(2);%求解杆4角加速度endfigure(4)plot(th2*180/pi,a3,th2*180/pi,a4);%绘制角加速度线图axis([0360-3060])text(30,18,'摇杆4角加速度(\alpha_4)')text(180,7,