四连杆机构分析代码动力学--精简

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)最短杆为连架杆，四杆机构为曲柄摇杆机构第二组(2代二套)四杆机构L1=125.36mm,L2=50.1mm,L3=109.8mm,L4=72.85mm最短杆长度+最长杆长度(125.36+50.1)≤其余两杆长度之和(109.8+72.85)最短杆为连架杆，四杆机构为曲柄摇杆机构第三组(3代)四杆机构L1=163.2mm,L2=61.6mm,L3=150mm,L4=90mm最短杆长度+最长杆长度(163.2+61.6)≤其余两杆长度之和(150+90)最短杆为连架杆，四杆机构为曲柄摇杆机构在如下图1所示的曲柄摇杆机构中，构件AB为曲柄，则B点应能通过曲柄与连杆两次共线的位置。1.3机构的数学模型的建立2图1机构结构简图1.3.1建立机构的闭环矢量位置方程在用矢量法建立机构的位置方程时，需将构件用矢量来表示，并作出机构的封闭矢量多边形。如图1所示，先建立一直角坐标系。设各构件的长度分别为L1、L2、L3、L4，其方位角为、、、。以各杆矢量组成一个封闭矢量多边形，即ABCDA。其个矢量之和必等于零。即：式1式1为图1所示四杆机构的封闭矢量位置方程式。对于一个特定的四杆机构，其各构件的长度和原动件2的运动规律，即为已知，而=0，故由此矢量方程可求得未知方位角、。角位移方程的分量形式为：式2闭环矢量方程分量形式对时间求一阶导数（角速度方程）为：式3其矩阵形式为：2222224344334433cossincoscossinsinLLLLLL式4联立式3两公式可求得：式5式6闭环矢量方程分量形式对时间求二阶导数（角加速度方程）矩阵形式为：3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）的线性方程组：式124式中，系列矩阵是一个阶方阵：式13的逆矩阵为;常数项b是一个n维矢量：式14因此，线性方程组解的矢量为：式15式11是求解连杆3和摇杆4角速度和角加速度的依据。2平面连杆机构的动力学分析机构的动力分析，主要是在运动学分析的基础上，由已知工作阻力，求出运动副的约束反力和驱动力(或力矩)，为选择和设计轴承、零部件强度的计算及选择原动机提供理论依据。本章以机构的组成原理为出发点，主要以应用最为广泛的平面连杆II级机构为分析对象，用复数向量推导出曲柄原动件、RRR杆组、RRP杆组、RPR杆组、PRP杆组和RPP杆组的动力学矩阵数学模型，并编制相应仿真M函数。在Matlab／simulink仿真平台，可以搭建所有平面连杆II级机构的动力学仿真模型并进行动力学仿真。2.1曲柄的动力学矩阵表达式曲柄AB复向量的模ri为常数,幅角θi为变量。质心到转动副A的距离为rci，质量为mi，绕质心的转动惯量为Ji,作用于质心上的外力为Fxi和Fyi,外力矩为Mi,曲柄与机架联接，转动副A的约束反力为Rxa和Rya,驱动力矩为M1.由理论力学可得：由运动学知识可推得：5合并整理得2.2RRRII级杆组动力学矩阵表达式如图所示，RRRII级杆组，分别以2个构件BC(长度为)和CD(长度为）为受力分析对象进行受力分析，其受力情况同曲柄，只是不受驱动力矩，则转动副B，C，D的约束反力推导如下。构件BC受力分析得对构件CD受力分析得由运动学可推得代入、化简、合并写成矩阵为63基于MATLAB程序设计四连杆机构的解析法同样可以用MATLAB的计算工具来求值，并结合MATLAB的可视化手段，把各点的计算值拟合成曲线，得到四连杆机构的运动仿真轨迹。3.1程序流程图74栏杆机各机型的分析结果4.12代1机构尺寸参数各构件的尺寸为r1=73.4mm,r2=103.4mm,r3=103.52mm,r4=125.36mm；质心为rc1=mm，rc2＝mm．rc3＝mm质量为m1＝kg，m2＝kg．m3＝kg；转动惯量为J1＝kg•m2，J2＝kg•m2，J3＝kg•m2，构件3的工作阻力矩M3＝N•m．顺时针方向，其他构件所受外力和外力矩(弹簧拉力大小及位置)构件1以等角速度5.326rad／s逆时针方向回转2代1型运动时间0.6s曲柄角速度t*180/12θ12为曲柄两极限点的转角范围1.352.21812=183.1θ31为摇杆两极限点转角范围θ31=169.6412-79.2815=90.36-0.06-0.04-0.0200.020.040.060.080.10.120.14-0.06-0.04-0.0200.020.040.060.080.1连杆3的几个位置点水平方向(m)垂直方向(m)图4-1连杆3的空间位置点050100150200250300350020406080100120140160主动件转角2(度)从动件角位移(度)角位移线图摇杆4角位移连杆3角位移406080100120140160180200220020406080100120140160180摇杆4角位移连杆3角位移X:216Y:169.6主动件转角2(度)从动件角位移(度)角位移线图X:35Y:79.28图4-2连杆3和摇杆4的角位移曲线图4-3工作区间内连杆3和摇杆4的角位移曲线曲柄两极限点的转角范围1.352.21812=183.1摇杆两极限点转角范围θ31=169.6412-79.2815=90.368050100150200250300350-10-505摇杆4角速度(4)连杆3角速度(3)X:110Y:3.898主动件转角2(度)从动件角速度(rads-1)角速度线图406080100120140160180200220-5-4-3-2-1012345摇杆4角速度(4)连杆3角速度(3)X:222Y:3.835X:121Y:3.934X:31Y:-0.721主动件转角2(度)从动件角速度(rads-1)角速度线图X:31Y:-4.293图4-4连杆3和摇杆4角速度曲线图4-5工作区间内连杆3和摇杆4角速度曲线工作区间内摇杆角速度最大值:曲柄转角121度，摇杆转角126.6188，摇杆角速度3.9335工作区间内摇杆角速度最小值:曲柄转角35度，摇杆转角79.2817，摇杆角速度-0.0196角速度变化量：3.9335+0.0196=3.9531050100150200250300350-100-80-60-40-20020406080100摇杆4角加速度(4)连杆3角加速度(3)主动件转角2(度)从动件角加速度(rads-2)角加速度线图X:15Y:81.09406080100120140160180200220-40-20020406080摇杆4角加速度(4)连杆3角加速度(3)X:31Y:57.51从动件角加速度从动件角加速度(rads-2)角加速度线图X:180Y:-22.81图3-6连杆3和摇杆4角加速度曲线图3-7工作区间内连杆3和摇杆4角加速度曲线工作区间内摇杆角加速度最大值:曲柄转角35度，摇杆转角79.2817，摇杆角加速度49.5924工作区间内摇杆角加速度最小值:曲柄转角181度，摇杆转角164.2275，摇杆角速度-22.8097角速度变化量：49.5924+22.8097=72.40214.22代2机构尺寸参数9各构件的尺寸为r1=50.1mm,r2=109.8mm,r3=72.85mm,r4=125.36mm；质心为rc1=mm，rc2＝mm．rc3＝mm质量为m1＝kg，m2＝kg．m3＝kg；转动惯量为J1＝kg•m2，J2＝kg•m2，J3＝kg•m2，构件3的工作阻力矩M3＝N•m．顺时针方向，其他构件所受外力和外力矩(弹簧拉力大小及位置)构件1以等角速度3.38594rad／s逆时针方向回转2代2型运动时间0.9s曲柄角速度t*180/12θ12为曲柄两极限点的转角范围4.2620112=174.6θ31为摇杆两极限点转角范围θ31=162.9068-75.5916=87.3152下表为曲柄转动一周,各参数变化量,角度间隔5度-0.04-0.0200.020.040.060.080.10.120.14-0.06-0.04-0.0200.020.040.060.08连杆3的几个位置点水平方向(m)垂直方向(m)图4-21连杆3的空间位置点050100150200250300350020406080100120140160摇杆4角位移连杆3角位移X:200Y:162.9X:25Y:75.61X:310Y:62.59主动件转角2(度)从动件角位移(度)角位移线图X:140Y:8.126406080100120140160180200220020406080100120140160180摇杆4角位移连杆3角位移X:199Y:162.9主动件转角2(度)从动件角位移(度)角位移线图X:35Y:76.43图4-22连杆3和摇杆4角位移曲线图4-23工作区间内连杆3和摇杆4角位移曲线曲柄两极限点的转角范围4.2620112=174.610摇杆两极限点转角范围θ31=162.9068-75.5916=87.3152050100150200250300350-5-4-3-2-1012345摇杆4角速度(4)连杆3角速度(3)X:335Y:-3.057X:115Y:2.423X:230Y:1.755主动件转角2(度)从动件角速度