RRR—PRP六杆机构的Simulink运动学仿真.docVIP

RRR—PRP六杆机构的Simulink运动学仿真 RRR—PRP六杆机构的Simulink运动学仿真 题目：由原动件(曲柄1)和1个RRR杆组、PRP杆组所组成的RRR—PRP六杆机构，各构件的尺寸为r1=140mm，r2=420mm，r3=420mm，ReD=420mm，LmD=0mm，ReA=0mm，LmA=420mm.，OK=180 mm复数向量坐标如下图1所示，构件l以等角速度10 rad／s逆时针方向回转，试求构件5和构件6的加速度、速度和位移。  图1RRR—PRP六杆机构 MATLAB版本为7.11.0.584 (R2010b) 将机构进行拆杆组分解 曲柄AB+第一个RRR杆组（BCD）+第二个PRP杆组（EFG） 2、确定初始条件 首先对机构ABCD部分进行初始位移、速度分析： 在MATLAB中建立位移函数 function -1.834352636620214 2.886984215567561 即杆2的初始角速度=-1.834352636620214； 杆3的初始角速度=2.886984215567561； 对机构EFG部分进行初始位移、速度分析: 构件5的位移： s5=(0.18-0.14*cos(45*pi/180))/cos(-y1(1)) s5 = 0.083806056003369 即构件5的位移s5= 0.083806056003369 构件6的位移： s6=0.42+0.14*sin(45*pi/180)-(0.18-0.14*cos(45*pi/180))*tan(-y1(1)) s6 = 0.497509214688385 即构件6的位移s6= 0.497509214688385 通过SolidWorks 2009 软件建立平面机构来分析速度 构件5的速度： v5=-y2(1)*0.72475802426 v5 = 1.329461792712988 v55=v5*sin(50pi/180) v55 = 1.029361907794855 v555=v55/cos(y1(1)*pi/180) v555 = 1.064955345590595 即构件5的速度ds5= 1.064955345590595 构件6的速度： v6=v5*sin(24pi/180) v6 = 0.549335282591239 v66=0.549335282591239/cos(y1(1)*pi/180) v66 = 0.568330284312065 即构件6的速度ds6= 0.568330284312065 3、建立Simulink仿真模型 图2 RRR—PRP六杆机构Simulink仿真模型 RRR—PRP六杆机构Simulink仿真模型如上图2所示，共有3个MATLAB函数模块。 crank.m函数模块的输入参数为曲柄AB的角位移、角速度和角加速度，输出参数曲柄端部(转动副B)的加速度的水平分量和垂直分量。 RRRki．m函数模块的输入参数为构件2（BC）和构件3（CD）的角位移、角速度和转动副B的加速度。 PRPki．m函数模块的输入参数为构件5和构件6的长度、位移和速度，构件2的角位移、角速度和角加速度和转动副B的加速度。 长度单位设置为m，角度单位为rad。设置仿真时间为1.4s,求解器选用ode45，步长选用变步长。 主要初值设置为： 构件1长度 r1=0.14 构件2长度 r2=0.42 构件3长度 r3=0.42 构件i长度 ri=0 构件j长度 rj=0 构件1角位移 theta-1= 0.785398163397448 构件1角速度 dtheta-1= 10 构件2角位移 theta-2= -0.259269450610624 构件2角速度 dtheta-2= -1.834352636620214 构件3角位移 theta-3=1.367111049465176 构件3角速度 dtheta-3= 2.886984215567561 构件6角位移 theta-s6= 1.570796326794897 构件5角位移 s5= 0.083806056003369 构件5角速度 ds5= 1.064955345590595 构件6位移 s6= 0.497509214688385 构件6速度 d6= 0.568330284312065 4、RRR-PRP六杆机构MATLAB仿真结果 在MAYLAB中执行 subplot(3,2,1); plot(tout,simout(:,4))