张海涛2013424064第三次作业修改张海涛2013424064第三次作业修改.doc

《张海涛2013424064第三次作业修改张海涛2013424064第三次作业修改.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、四连杆机构动力学分析上次第三次作业的修改：由于上次作业，matlab中各杆的转动惯量使用计算，但是ADAMS中的杆件不是严格长方体杆件，导致ADAMS计算的各杆的转动惯量与matlab不一样，导致matlab的仿真结果与ADAMS相差过大。由算出的各杆的转动惯量分别为2116、69026、40016，而ADAMS算出的转动惯量为3246.2937888297、78327.27911、58297.131592，明显可以看出matlab中算出的转动惯量小，所以matlab算出的各杆转动速度比ADAMS快，曲线提前一定相位，如下图所示。图0转动惯量使用时的仿真曲线

2、对比目录一、动力学分析31、外力矢量32、外力所做的虚功33、系统的约束方程34、雅克比矩阵35、46、系统的动力学方程4二、说明5三、Matlab与ADAMS仿真曲线对比8附录一：求雅克比矩阵的Matlab程序16附录二：求的Matlab程序16附录三：用欧拉迭代法进行动力学分析的Matlab程序17使用ADAMS建立如图1所示的四连杆机构，二杆长150mm，三杆长500mm，四杆长450mm，二杆所受的转矩为40，三个杆都受重力，二杆初始角度为90度。运用欧拉法，用Matlab对该系统进行动力学分析，并与ADAMS仿真结果进行对比。二杆力矩三杆四杆图1四

3、杆机构一、动力学分析1、外力矢量2、外力所做的虚功3、系统的约束方程4、雅克比矩阵由系统的约束方程，求雅克比矩阵得。具体程序见附录一，可得雅克比矩阵：5、具体程序见附录二，可得：6、系统的动力学方程其中：使用欧拉迭代法就可以计算出各个曲线，循环运算1666667次，具体Matlab程序见附录三。二、说明1、由于本程序采用欧拉迭代法求解，所以计算过程有误差。当仿真时间越长时，与ADAMS曲线相差越大，但曲线形状一致，只有相位提前。2、Matlab的仿真时间间隔越小，仿真结果越真实。但当间隔非常小时，再降低时间间隔仿真结果相差不大。3、当ADMAS的仿真步数很大

4、时，再增加仿真步数对结果影响很小。4、当改变ADAMS的误差时，仿真结果也有差别，如图2-6所示。当误差小于10-6时，再减小误差，曲线差别很小。5、由于在ADAMS中只能由三个点测角度，三杆的角度无法测量。所以再建立两个辅助杆，并设这两个辅助杆的质量和转动惯量都为零，用于测量三杆的角度，如图2中辅助杆1和辅助杆2。再就可以测量由E、F、G点组成的角度作为三杆的角度。点F点G点E辅助杆2辅助杆1图2辅助测量杆图3误差为1*10-3二杆质心X方向位移变化曲线图4误差为5*10-3二杆质心X方向位移变化曲线图5误差为1*10-4二杆质心X方向位移变化曲线图6误差

5、为1*10-6二杆质心X方向位移变化曲线图7误差为1*10-9二杆质心X方向位移变化曲线三、Matlab与ADAMS仿真曲线对比局部放大图8铰链1处作用力变化曲线图9上图局部放大的曲线对比局部放大图10铰链2处作用力变化曲线图11上图局部放大的曲线对比图12铰链3处作用力变化曲线图13铰链4处作用力变化曲线图14二杆位置变化曲线图15三杆位置变化曲线图16四杆位置变化曲线图17二杆速度变化曲线图18三杆速度变化曲线图19四杆速度变化曲线图20二杆加速度变化曲线图21三杆加速度变化曲线图22四杆加速度变化曲线附录一：求雅克比矩阵的Matlab程序clcclea

6、rsymsl2l3l4l5;symsm2m3m4;symsx4x5x6x7x8x9x10x11x12g=[x4-l2/2*cos(x6);x5-l2/2*sin(x6);x4+l2/2*cos(x6)-x7+l3/2*cos(x9);x5+l2/2*sin(x6)-x8+l3/2*sin(x9);x7+l3/2*cos(x9)-x10+l4/2*cos(x12);x8+l3/2*sin(x9)-x11+l4/2*sin(x12);x10+l4/2*cos(x12)-l5;x11+l4/2*sin(x12);]cq=jacobian(g,[x4x5x6x7x8

7、x9x10x11x12])m=[m200000000;0m20000000;00m2/12*(l2)^2000000;000m300000;0000m30000;00000m3/12*(l3)^2000;000000m400;0000000m40;00000000m4/12*(l4)^2]o=zeros(8,8);cq1=cq.';A=[mcq1;cqo]附录二：求的Matlab程序clcclearsymsl2l3l4l5;symsm2m3m4;symsx4x5x6x7x8x9x10x11x12symsy4y5y6y7y8y9y10y11y12h=[y4+y

8、6*l2/2*sin(x6);y5-y6*l2/2*