哈工大机械原理大作业 第25题 连杆机构设计
机械原理大作业一连杆机构参考例子

机械原理大作业课程名称:机械原理设计题目:机械原理大作业院系:汽车工程学院车辆工程班级:1101201姓名:。
学号:。
指导教师:游斌弟大作业1 连杆机构运动分析1、运动分析题目如图所示机构,已知机构各构件的尺寸为280mm AB =,350mm BC =,320mm CD =,160mm AD =,175mm BE = 220mm EF =,25mm G x =,80mm G y =,构件1的角速度为110rad/s ω=,试求构件2上点F 的轨迹及构件5的角位移、角速度和角加速度,并对计算结果进行分析。
图 12、对机构进行结构分析该机构由I 级杆组RR (原动件1)、II 级杆组RRR (杆2、杆3)和II 级杆组RPR (滑块4及杆5)组成。
I 级杆组RR ,如图2所示;II 级杆组RRR ,如图3所示;II 级杆组RPR ,如图4所示。
图2 图 2图4 3、建立坐标系建立以点A为原点的固定平面直角坐标系4、各基本杆组运动分析的数学模型(1)同一构件上点的运动分析:如图5所示的构件AB,,已知杆AB 的角速度=10/rad s ω,AB 杆长i l =280mm,可求得B 点的位置B x 、B y ,速度xB v 、yB v ,加速度xB a 、yB a 。
=cos =280cos B i x l ϕϕ; =sin =280sin B i y l ϕϕ;图 3==-sin =-BxB i B dx v l y dt ωϕω; ==cos =;B yB i B dyv l x dt ωϕω222B 2==-cos =-BxB i d x a l x dt ωϕω;2222==-sin =-ByB i Bd y a l y dtωϕω。
图 4(2)RRRII 级杆组的运动分析如图6所示是由三个回转副和两个构件组成的II 级组。
已知两杆的杆长2l 、3l 和两个外运动副B 、D 的位置(B x 、B y 、D x 、D y )、速度( xB yB xD yD v v v v 、、、 ) 图6和加速度(xB yB xD yD a a a a 、、、)。
哈工大机械原理大作业1 连杆机构运动分析 牛头刨床 学号5号

大作业1 连杆机构运动分析运动分析题目图1-13 牛头刨床机构简图一、按比例画出机构运动简图见A3图纸二、机构的结构分析及基本杆组的划分(1)机构的结构分析该机构为平面连杆机构,活动构件数n=5;有A、B、C、E、F五个转动副和B、F两个移动副,没有高副,故低副数PL=7,高副数PH=0。
机构的自由度F=3n-2PL-PH=3*5-2*7=1。
机构有确定运动轨迹。
(2)基本杆组的划分图1-13中,原动件构件1为Ⅰ级杆组,移除之后按照拆杆组法对机构进行拆分,可得由滑块2和杆3组成的RPRⅡ级杆组,以及由杆4和滑块5组成的RRPⅡ级杆组。
拆分图如下:图1-13-1 Ⅰ级杆组图1-13-2 RPRⅡ级杆组图1-13-3 RRPⅡ级杆组三、各杆组的运动分析数学模型图1-13-1 Ⅰ级杆组 x B =x A +l AB cos ϕy B =y A +l AB sin ϕ图1-13-2 RPR Ⅱ级杆组 x E =x B +(l CE -s)cos ϕj y E =y B +(l CE -s)sin ϕjϕj =arctan(B 0/A 0)A 0=xB -x CB 0=y B -y CS =2020B A + 图1-13-3 RRP Ⅱ级杆组F 的位置 x F =x E -l EF cos ϕiy F =y E +l EF sin ϕi ϕi =arcsin(A 1/l EF )A 1=H-H 1-y EF 的速度F 的加速度四、建立坐标系以A 为原点建立坐标系,则A (0,0)五、编程使用MATLAB 编写程序如下lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);a=0:0.001:0.5;xF=subs(xF,t,a);plot(a,xF)title('位移x随时间t变化图像') xlabel('t(s)'),ylabel(X')lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);vF=diff(xF,t);a=0:0.001:0.5;vF=subs(vF,t,a);plot(a,vF)title('速度v随时间t变化图像') xlabel('t(s)'),ylabel(V')lAB=108;lCD=620;lEF=300;H1=350;H=635;xC=0;yC=-350;syms t;f=(255*pi/30)*t;xB=lAB*cos(f);yB=lAB*sin(f);A0=xB-xC;B0=yB-yC;S=sqrt(A0^2+B0^2);f1=atan(B0/A0);xE=xB+(lAB-S)*cos(f1);yE=yB+(lAB-S)*sin(f1);A1=H-H1-yE;f2=asin(A1/lEF);xF=xE-lEF*cos(f2);yF=yE+lEF*sin(f2);aF=diff(xF,t,2);a=0:0.001:0.5;aF=subs(aF,t,a);plot(a,aF)title('加速度a随时间t变化图像') xlabel('t(s)'),ylabel(A')六、位移、速度、加速度图像位移xE随时间t变化图像速度vE随时间t变化图像加速度aE随时间t变化图像七、运动分析由MATLAB编程后绘制所得图线分析可知,牛头刨床的刨头F行程约为500mm,在0.06s~0.17s之间速度和加速度波动较大,在0.17s~0.29s 之间速度和加速度波动较小,这就保证了刨头在空行程时有急回运动,而在工作行程是由较均匀的切削速度,运动连续,且运动性能良好稳定。
哈工大机械原理大作业——连杆机构运动分析16___2014

Harbin Institute of Technology机械原理大作业——连杆机构运动分析课程名称:机械原理院系:能源科学与工程学院班级:完成者:学号:题号: 16任课教师:丁刚完成内容:在完成题目计算要求的同时,扩展了内容,程序为该结构的通用程序,可解决机构在不同条件下的运动情况,文本最末为几种情况的分析哈尔滨工业大学16、如图所示机构,已知机构各构件的尺寸为,试求构件5的角位移、角速度和角加速度,并对计算结构进行分析。
(1)、结构分析从侧面看原机构为此机构分为级杆组(原动件1),级杆组RRP(2号套筒、3号杆),级杆组RRP(4号套筒、5号杆)(2)、建立坐标系(3)、各个杆组的运动分析采用逆推法,从RRP杆组(4号套筒、5号杆)开始分析已知,,,,现在假定已知,,其中,,,即消去,可得可求得,也可以通过书上3-23式求得通过正弦定理可求得再来看看角速度关系对于加速度,有如下关系其中到此4、5杆就分析完毕了,别忘记之前的假设,我假设了已知,,为求,,,现在来分析RRP杆组(2号套筒、3号杆)已知,,,已知,,,,其中,,,即消去,可得反解,即可求得,也可以通过书上3-23式求得通过正弦定理可求得继续,我们来看看角速度关系对于加速度,有如下关系其中现在,只需将所求得的,,和,,关联起来这是同一根杆,,,现在来看,,,由题目得,,和是未知的,但不影响整体,不然给一个初值,,当然,这是可以随意更改的。
基于以上的基本原理,matlab R2012b程序如下syms theta theta1 theta2 lamuda lamuda1 lamuda2 sigma sigma1 sigma2 beta beta1 beta2 l1 l11 l2 l21 t output itheta1=10;theta2=0;i=0;for theta3=60:420theta=theta3/180*pi;beta=asin((100/200)*sin(theta))+theta;l1=0.2*sin(beta)/sin(theta);beta1=(-theta1*(l1*sin(theta))*sin(theta)+theta1*(l1*cos(theta))*cos(theta))/(0.2*(sin(theta)*sin(b eta)+cos(theta)*cos(beta)));l11=-(theta1*(l1*sin(theta))*l1*cos(beta)+theta1*(l1*cos(theta))*l1*sin(beta))/(0.2*(sin(theta)*si n(beta)+cos(theta)*cos(beta)));C=(theta1^2)*0.2*cos(beta)-theta2*l1*sin(theta)-(theta1^2)*l1*cos(theta)-2*l11*theta1*sin(theta) ;D=(theta1^2)*0.2*cos(beta)+theta2*l1*sin(theta)-(theta1^2)*l1*cos(theta)+2*l11*theta1*sin(thet a);beta2=(-C*sin(theta)+D*cos(theta))/(0.2*(sin(theta)*sin(beta)+cos(theta)*cos(beta)));lamuda=beta-pi/2;lamuda1=beta1;lamuda2=beta2;sigma=asin((100/200)*sin(lamuda))+lamuda;l2=0.2*sin(sigma)/sin(lamuda);sigma1=(-lamuda1*(l2*sin(lamuda))*sin(lamuda)+lamuda1*(l2*cos(lamuda))*cos(lamuda))/(0.2 *(sin(lamuda)*sin(sigma)+cos(lamuda)*cos(sigma)));l21=-(lamuda1*(l2*sin(lamuda))*l2*cos(sigma)+lamuda1*(l2*cos(lamuda))*l2*sin(sigma))/(0.2* (sin(lamuda)*sin(sigma)+cos(lamuda)*cos(sigma)));A=(lamuda1^2)*0.2*cos(sigma)-lamuda2*l2*sin(lamuda)-(lamuda1^2)*l2*cos(lamuda)-2*l21*la muda1*sin(lamuda);B=(lamuda1^2)*0.2*cos(sigma)+lamuda2*l2*sin(lamuda)-(lamuda1^2)*l2*cos(lamuda)+2*l21*l amuda1*sin(lamuda);sigma2=(-A*sin(lamuda)+B*cos(lamuda))/(0.2*(sin(lamuda)*sin(sigma)+cos(lamuda)*cos(sigma )));i=i+1;output(i,1)=fix(theta/pi*180);output(i,2)=fix(sigma/pi*180);output(i,3)=fix(sigma1);output(i,4)=fix(sigma2);endoutputa=output(:,1);b=output(:,2);c=output(:,3);d=output(:,4);h1=plot(a,b);hold on;h2=plot(a,c);hold on;h3=plot(a,d);hold on;set(h1,'color',[1 0 0],'linewidth',2);set(h2,'color',[0 1 1],'linewidth',1);set(h3,'color',[0 0 1],'linewidth',2);m=legend('角位移','角速度','角加速度');x label('θ');title('平面连杆机构运动分析');figure;h1=plot(a,b);hold on;x label('θ');ylabel('角位移');title('平面连杆机构运动角度——角位移图');figure;h2=plot(a,c);hold on;x label('θ');ylabel('角速度');title('平面连杆机构运动角度——角速度图'); figure;h3=plot(a,d);hold on;x label('θ');ylabel('角加速度');title('平面连杆机构运动角度——角加速度图');汇总图各自的图像结果分析,上面的图形只是在一个初值,的条件下得出的,为了能解决所有问题,修改程序如下syms theta theta1 theta2 lamuda lamuda1 lamuda2 sigma sigma1 sigma2 beta beta1 beta2 l1 l11 l2 l21 t output iprompt={'输入:', '输入' ,'输入' };%设置提示字符串name='输入初值';%设置标题 numlines=1;%指定输入数据的行数 defAns={'60','10','0'};%设定默认值 Resize='on';%设定对话框尺寸可调节answer=inputdlg(prompt,name,numlines,defAns,'on');%创建输入对话框 h= str2num(answer{1}); theta1= str2num(answer{2}); theta2= str2num(answer{3}); i=0;for theta3=h:(360+h) theta=theta3/180*pi;beta=asin((100/200)*sin(theta))+theta; l1=0.2*sin(beta)/sin(theta);beta1=(-theta1*(l1*sin(theta))*sin(theta)+theta1*(l1*cos(theta))*cos(theta))/(0.2*(sin(theta)*sin(b eta)+cos(theta)*cos(beta)));l11=-(theta1*(l1*sin(theta))*l1*cos(beta)+theta1*(l1*cos(theta))*l1*sin(beta))/(0.2*(sin(theta)*si n(beta)+cos(theta)*cos(beta)));C=(theta1^2)*0.2*cos(beta)-theta2*l1*sin(theta)-(theta1^2)*l1*cos(theta)-2*l11*theta1*sin(theta) ;D=(theta1^2)*0.2*cos(beta)+theta2*l1*sin(theta)-(theta1^2)*l1*cos(theta)+2*l11*theta1*sin(thet a);beta2=(-C*sin(theta)+D*cos(theta))/(0.2*(sin(theta)*sin(beta)+cos(theta)*cos(beta)));lamuda=beta-pi/2;lamuda1=beta1;lamuda2=beta2;sigma=asin((100/200)*sin(lamuda))+lamuda;l2=0.2*sin(sigma)/sin(lamuda);sigma1=(-lamuda1*(l2*sin(lamuda))*sin(lamuda)+lamuda1*(l2*cos(lamuda))*cos(lamuda))/(0.2 *(sin(lamuda)*sin(sigma)+cos(lamuda)*cos(sigma)));l21=-(lamuda1*(l2*sin(lamuda))*l2*cos(sigma)+lamuda1*(l2*cos(lamuda))*l2*sin(sigma))/(0.2* (sin(lamuda)*sin(sigma)+cos(lamuda)*cos(sigma)));A=(lamuda1^2)*0.2*cos(sigma)-lamuda2*l2*sin(lamuda)-(lamuda1^2)*l2*cos(lamuda)-2*l21*la muda1*sin(lamuda);B=(lamuda1^2)*0.2*cos(sigma)+lamuda2*l2*sin(lamuda)-(lamuda1^2)*l2*cos(lamuda)+2*l21*l amuda1*sin(lamuda);sigma2=(-A*sin(lamuda)+B*cos(lamuda))/(0.2*(sin(lamuda)*sin(sigma)+cos(lamuda)*cos(sigma )));i=i+1;output(i,1)=fix(theta/pi*180);output(i,2)=fix(sigma/pi*180);output(i,3)=fix(sigma1);output(i,4)=fix(sigma2);endoutputa=output(:,1);b=output(:,2);c=output(:,3);d=output(:,4);h1=plot(a,b);hold on;h2=plot(a,c);hold on;h3=plot(a,d);hold on;set(h1,'color',[1 0 0],'linewidth',2);set(h2,'color',[0 1 1],'linewidth',1);set(h3,'color',[0 0 1],'linewidth',2);m=legend('角位移','角速度','角加速度');x label('θ');title('平面连杆机构运动分析');figure;h1=plot(a,b);hold on;xlabel('θ');y label('角位移');title('平面连杆机构运动角度——角位移图');figure;h2=plot(a,c);hold on;xlabel('θ');y label('角速度');title('平面连杆机构运动角度——角速度图');figure;h3=plot(a,d);hol d on;xlabel('θ');y label('角加速度');title('平面连杆机构运动角度——角加速度图');这样,在运行程序时就会弹出一个如下图所示的对话框,可以任意给定初值,解决不同问题。
机械设计基础.连杆机构

不同类型的连杆机构
四杆机构
由四个杆件和四个关节组成, 常用于平移和旋转运动的转 换。
曲柄滑块机构
由曲柄、连杆和滑块组成, 用于将旋转运动转化为直线 运动。
双曲柄机构
由两个曲柄、连杆和滑块组 成,常用于往复运动的转换。
连杆机构的运动分析
连杆机构的运动分析包括确定杆件的位置、速度和加速度,以及分析关联元 件在不同条件下的运动特性。
连杆机构的设计要点
1 强度和刚度
连杆材料的选择和尺寸设 计要满足强度和刚度的需 求。
2 运动关节的设计
3 减少运动噪音
关节的设计要考虑到磨损、 摩擦和润滑等因素。
采用合适的材料和润滑方 式可以有效减少连杆机构 的运动噪音。
常见的连杆机构应用举例
内燃机
连杆机构将活塞运动转化为曲轴的旋转运动,实现 发动机的工作。
起重机
连上下运动。
机器人臂
连杆机构用于机器人臂的关节,实现灵活的运动和 抓取功能。
自行车
脚踏板和齿轮系统通过连杆机构将人的脚力转化为 车轮的旋转运动。
总结
连杆机构是机械设计中不可或缺的一部分。通过了解其定义、原理、分析和设计要点,我们能够更好地应用它 们于实践中。
机械设计基础.连杆机构
在机械设计中,连杆机构是一种关键的元件。了解连杆机构的定义、原理以 及运动分析,让我们能够更好地设计和应用机械系统。
机械设计基础
机械设计是将工程原理和技术应用于机械系统的过程。它涉及材料选择、构 建模型、工程计算以及系统优化。
连杆机构的定义与原理
连杆机构是由连接杆件和关节组成的机械系统。它基于欧拉运动方程和级数法等原理,用于转换和传递运动。
哈工大机械原理大作业24题

班级1013102学号 6机械原理大作业说明书题目 1、连杆机构运动分析2、凸轮机构设计3、齿轮传动设计学生姓名1连杆机构运动分析1.设计题目:一、先建立如下坐标系:二、划分杆组如下,进行结构分析:该机构由I级杆组RR(如图1)、II级杆组RPR(如图2、3)和II级杆组RRP(如图4)组成。
(1)(2)(3)(4)三、运动分析数学模型:(1)同一构件上点的运动分析:如右图所示的原动件1,已知杆1的角速度=10/rad s ω,杆长1l =170mm,A y =0,A x =110mm 。
可求得下图中B 点的位置B x 、B y ,速度xB v 、yB v ,加速度xB a 、yB a 。
θcos 1l xB =,θsin 1l yB =θωυsin 1l xB -=,θωυcos 1l yB =,222B2==-cos =-B xB i d x a l x dt ωϕω2222==-sin =-B yB i B d y a l y dtωϕω。
(2)RPRII 级杆组的运动分析:a. 如右图所示是由2个回转副和1个移动副组成的II 级组。
已知两个外运动副C 、B 的位置(B x 、B y 、c x =110mm 、C y =0)、速度(xB υ,yB υ,xC υ=0,yC υ=0)和加速度(0,0,,==yC xC yB xB a a a a )。
可确定下图中D 点的位置、速度和加速度。
确定构件3的角位移1ϕ、角速度1ω、角加速度1α。
1sin 31..ϕϕl x dtdx C B-= 1sin 131cos 13.....2ϕϕϕϕl l x dt x d C B --= 1cos 31..ϕϕl y dtdy C B+=1cos 131sin 13.....2ϕϕϕϕl l y dt y d C B +-= 根据关系:1111d 122..11.αϕϕωϕϕ====dtd dt , 故可得出:D x =)1cos(4βϕ++l x C D y =)1sin(4βϕ++l y Cb. 如右图所示是由2个回转副和1个移动副组成的II级组。
机械原理大作业一连杆机构参考例子

机械原理大作业课程名称:机械原理设计题目:机械原理大作业院系:汽车工程学院车辆工程班级:1101201姓名:。
学号:。
指导教师:游斌弟大作业1 连杆机构运动分析1、运动分析题目如图所示机构,已知机构各构件的尺寸为280mm AB =,350mm BC =,320mm CD =,160mm AD =,175mm BE = 220mm EF =,25mm G x =,80mm G y =,构件1的角速度为110rad/s ω=,试求构件2上点F 的轨迹及构件5的角位移、角速度和角加速度,并对计算结果进行分析。
图 12、对机构进行结构分析该机构由I 级杆组RR (原动件1)、II 级杆组RRR (杆2、杆3)和II 级杆组RPR (滑块4及杆5)组成。
I 级杆组RR ,如图2所示;II 级杆组RRR ,如图3所示;II 级杆组RPR ,如图4所示。
图2 图 2图4 3、建立坐标系建立以点A为原点的固定平面直角坐标系4、各基本杆组运动分析的数学模型(1)同一构件上点的运动分析:如图5所示的构件AB,,已知杆AB 的角速度=10/rad s ω,AB 杆长i l =280mm,可求得B 点的位置B x 、B y ,速度xB v 、yB v ,加速度xB a 、yB a 。
=cos =280cos B i x l ϕϕ; =sin =280sin B i y l ϕϕ;图 3==-sin =-BxB i B dx v l y dt ωϕω; ==cos =;B yB i B dyv l x dt ωϕω222B 2==-cos =-BxB i d x a l x dt ωϕω;2222==-sin =-ByB i Bd y a l y dtωϕω。
图 4(2)RRRII 级杆组的运动分析如图6所示是由三个回转副和两个构件组成的II 级组。
已知两杆的杆长2l 、3l 和两个外运动副B 、D 的位置(B x 、B y 、D x 、D y )、速度( xB yB xD yD v v v v 、、、 ) 图6和加速度(xB yB xD yD a a a a 、、、)。
哈工大机械原理大作业连杆机构9
机械原理大作业大作业一:连杆机构运动分析学生姓名:学号:指导教师:丁刚完成时间:机电工程学院机械设计系制二〇一八年四月连杆机构运动分析1题目(9)图1 设计题目在图1所示的机构中,已知l AB=60mm,l BC=180mm,l DE=200mm,l CD=120mm,l EF=300mm,h=80mm,h1=85mm,h2=225mm,构件1以等角速度ω1=100rad/s 转动。
求在一个运动循环中,滑块5的位移、速度和加速度曲线。
2分析结构1、杆1为RR主动件,绕A以ω1 转动,自由度1.2、4杆和滑块5为RRP II级杆组.,自由度0.3、2,3杆组成II级杆组RRR,自由度0.总共自由度为F=5*3-2*7=1 .由上述的杆组类型,确认出所需运动分析数学模型:同一构件上的点、RRP、RRR。
3.杆组法对平面连杆机构进行运动分析3.1对主动件杆1 RR I级构件的分析主动杆1转角:φ= [0°,360°) δ=0°,则φ’=ω1=100 rad/s角加速度φ’’=0.已知h2=225mm, h=80mm, l AB=60mm 所以A(225mm,80mm)A点速度(0,0),加速度(0,0)B点位置(x A+l AB*cos(φ), Y A+l AB*sin(φ))B点速度(-l AB*sin(φ), l AB*cos(φ)),加速度(-l AB*cos(φ), -l AB*sin(φ))3.2RRRII 级杆组分析(模型参考教材P37-38)图3 如图所示两个构件组成II 级杆组。
已知了B 的位置(x B ,y B )= (x A +l AB *cos(φ), Y A +l AB *sin(φ)),速度(x ’B ,y ’B ) 和加速度(x ’’B ,y ’’B ), 已知运动副D (0,0), 还可知,x ’D =y ’D =0, x ’’D =y ’’D =0. l BC =180 mm, l CD = 120mm所以,x c =x D +l CD *cos(φi)= x B +l BC *cos(φj) y c =x D +l CD *sin(φi)= x B +l BC *sin(φi) 对于φ的求解: A 0=2*l CD (x B -x D ) B 0=2*l BC (y B -y D ) C 0=l CD 2+ l BD 2- l BC 2为了保证机构的装配正常:l BD ≤l CD + l BC AND l BD ≥Abs (l CD - l BC )可求3杆的转角φi=2*arctan((B 0±sqrt (A 02 + B 02- C 02))/(A 02+ C 02)),角速度w3=φi ’和角加速度α3= φi ’’3.3 同一构件上的点(模型参考书P35-36)Φiφjφi已知D(0,0),速度(0,0),加速度(0,0),3杆转角φi 角速度φi’角加速度φi’’,Φi和它的导数在3.2都有体现LDE= 200mm可求出E的坐标,速度,加速度.x E =x c+lCE*cos(φi)y E =x C+lCE*sin(φi)同样地,速度、加速度通过求导即可得出算式,可以编出程序。
哈工大机械原理大作业连杆14-2015版
Harbin Institute of Technology机械原理大作业连杆课程名称:机械原理设计题目:连杆机构运动分析院系:xx学院班级:xxxxxxxx姓名:xxxxxxxxx学号:xxxxxxxxxxx指导教师:焦映厚陈照波设计时间:2015年xxxxxxxxxx哈尔滨工业大学题目14如图所示机构,已知各机构尺寸为,,,,构件l的角速度为,试求构件2上点D的轨迹及构建5的角位移、角速度和角加速度,并对计算结果进行分析。
(题中构建尺寸满足)对机构进行结构分析该机构由原动件AB(Ⅰ级组),CD(RPRⅡ级杆组)和DE(RPRⅡ级杆组)组成,拆分如下:RR杆组RPR杆组RPR杆组数学模型计算对B点:对D点(∠BAC=,∠BCA=)流程框图开始编程及代码根据以上计算采用vb编程编程代码如下:Option ExplicitPrivate lab, lac, lbd, lce, lbc, lae, pi, sita As DoubleFunction arcsin(x As Double) As DoubleIf x >= -1 And x < -0.5 Then arcsin = -Atn(Sqr(1 - x * x) / x) - 2 * Atn(1)If x >= -0.5 And x <= 0.5 Then arcsin = Atn(x / Sqr(1 - x * x))If x > 0.5 And x <= 1 Then arcsin = -Atn(Sqr(1 - x * x) / x) + 2 * Atn(1) End FunctionFunction arccos(x As Double) As DoubleIf x >= -1 And x < -0.5 Then arccos = Atn(Sqr(1 - x * x) / x) + 4 * Atn(1)If x >= -0.5 And x <= 0.5 Then arccos = -Atn(x / Sqr(1 - x * x)) + 2 * Atn(1) If x > 0.5 And x <= 1 Then arccos = Atn(Sqr(1 - x * x) / x)End FunctionPrivate Sub Command1_Click()Line1.Visible = FalseLine2.Visible = FalseLine3.Visible = False'画坐标系lab = 200lac = 400lbd = 700lae = 800lce = 400pi = 3.1416Dim i, j As IntegerDim beta, xb, yb, xd, yd As DoublePicture1.BackColor = vbWhitePicture1.FillStyle = 0Picture1.Scale (-400, 440)-(1500, -440)Picture1.Line (-400, 0)-(1500, 0)Picture1.Line (0, -440)-(0, 440)For i = 0 To 22Picture1.CurrentX = -68Picture1.CurrentY = 440 - 40 * iPicture1.Print " " & 440 - 40 * i & " "Picture1.DrawStyle = 2Picture1.Line (-400, 440 - 40 * i)-(1500, 440 - 40 * i)NextFor j = 0 To 70Picture1.CurrentX = 1520 - 80 * jPicture1.CurrentY = -5Picture1.Print 1520 - 80 * jPicture1.DrawStyle = 2Picture1.Line (1520 - 80 * j, 440)-(1520 - 80 * j, -440) Next'画D点位移图线For sita = 0 To 2 * pi Step 0.001xb = lab * Cos(sita)yb = lab * Sin(sita)lbc = Sqr(lab ^ 2 + lac ^ 2 - 2 * lab * lac * Cos(sita)) beta = arcsin(lab * Sin(sita) / lbc)xd = xb + lbd * Cos(beta)yd = yb - lbd * Sin(beta)Picture1.PSet (xd, yd)NextEnd SubPrivate Sub Command2_Click()‘D点角位移图Line1.Visible = FalseLine2.Visible = FalseLine3.Visible = Falselab = 200lac = 400lbd = 700lae = 800lce = 400pi = 3.1416Dim i, j As IntegerDim beta, xb, yb, xd, yd, gama, gama1 As Double Picture1.BackColor = vbWhitePicture1.FillStyle = 0Picture1.Scale (-20, 360)-(360, -80)Picture1.Line (-0, 0)-(360, 0)Picture1.Line (0, 360)-(0, -80)For i = 0 To 22Picture1.CurrentX = -15Picture1.CurrentY = 360 - 20 * iPicture1.Print " " & 360 - 20 * i & " "Picture1.DrawStyle = 2Picture1.Line (0, 360 - 20 * i)-(360, 360 - 20 * i)NextFor j = 0 To 18Picture1.CurrentX = 360 - 20 * jPicture1.CurrentY = -5Picture1.Print 360 - 20 * jPicture1.DrawStyle = 2Picture1.Line (360 - 20 * j, 360)-(360 - 20 * j, -80)NextFor sita = 0 To 360 Step 0.01xb = lab * Cos(sita * pi / 180)yb = lab * Sin(sita * pi / 180)lbc = Sqr(lab ^ 2 + lac ^ 2 - 2 * lab * lac * Cos(sita * pi / 180)) beta = arcsin(lab * Sin(sita * pi / 180) / lbc)xd = xb + lbd * Cos(beta)yd = yb - lbd * Sin(beta)gama = arccos((xd - 800) / Sqr((xd - 800) ^ 2 + yd ^ 2))If yd <= 0 Thengama1 = 2 * pi - gamaElsegama1 = gamaEnd IfPicture1.PSet (sita, gama1 * 180 / pi)NextEnd SubPrivate Sub Command3_Click()‘画角角速度Line1.Visible = FalseLine2.Visible = FalseLine3.Visible = Falselab = 200lac = 400lbd = 700lae = 800lce = 400pi = 3.1416Dim beta, xb, yb, xd, yd, gama, gama1, beta0, xb0, yb0, xd0, yd0, gama0, gama10, lbc0, w, sita0 As DoubleDim i, j As IntegerPicture1.BackColor = vbWhitePicture1.FillStyle = 0Picture1.Scale (-20, 90)-(360, -540)Picture1.Line (-0, 0)-(360, 0)Picture1.Line (0, 360)-(0, -540)For i = 0 To 22Picture1.CurrentX = -15Picture1.CurrentY = 90 - 30 * iPicture1.Print " " & 90 - 30 * i & " "Picture1.DrawStyle = 2Picture1.Line (0, 90 - 30 * i)-(360, 90 - 30 * i)NextFor j = 0 To 18Picture1.CurrentX = 360 - 20 * jPicture1.CurrentY = -5Picture1.Print 360 - 20 * jPicture1.DrawStyle = 2Picture1.Line (360 - 20 * j, 360)-(360 - 20 * j, -540)NextFor sita = 0 To 360 Step 0.01xb = lab * Cos(sita * pi / 180)yb = lab * Sin(sita * pi / 180)lbc = Sqr(lab ^ 2 + lac ^ 2 - 2 * lab * lac * Cos(sita * pi / 180))beta = arcsin(lab * Sin(sita * pi / 180) / lbc)xd = xb + lbd * Cos(beta)yd = yb - lbd * Sin(beta)gama = arccos((xd - 800) / Sqr((xd - 800) ^ 2 + yd ^ 2))If yd <= 0 Thengama1 = 2 * pi - gamaElsegama1 = gamaEnd Ifsita0 = sita + 0.001xb0 = lab * Cos(sita0 * pi / 180)yb0 = lab * Sin(sita0 * pi / 180)lbc0 = Sqr(lab ^ 2 + lac ^ 2 - 2 * lab * lac * Cos(sita0 * pi / 180)) beta0 = arcsin(lab * Sin(sita0 * pi / 180) / lbc0)xd0 = xb + lbd * Cos(beta0)yd0 = yb - lbd * Sin(beta0)gama0 = arccos((xd0 - 800) / Sqr((xd0 - 800) ^ 2 + yd0 ^ 2))If yd0 <= 0 Thengama10 = 2 * pi - gama0Elsegama10 = gama0End Ifw = (gama10 - gama1) / (pi / 180000)Picture1.PSet (sita, w * 180 / pi)NextEnd SubPrivate Sub Command4_Click()‘画运动图Timer1.Enabled = TrueEnd SubPrivate Sub Command5_Click()‘画角加速度Line1.Visible = FalseLine2.Visible = FalseLine3.Visible = Falselab = 200lac = 400lbd = 700lae = 800lce = 400pi = 3.1416Dim w2, a, sita0, sita02, gama102, gama12, sita2, beta, xb, yb, xd, yd, gama, gama1, beta0, xb0, yb0, xd0, yd0, gama0, gama10, lbc0, w1, beta2, xb2, yb2, xd2, yd2, lbc2, lbc02, gama2, beta02, xb02, yb02, xd02, yd02, gama02 As DoubleDim i, j As IntegerPicture1.BackColor = vbWhitePicture1.FillStyle = 0Picture1.Scale (-20, 140000)-(360, -140000)Picture1.Line (-0, 0)-(360, 0)Picture1.Line (0, 140000)-(0, -140000)For i = 0 To 22Picture1.CurrentX = -15Picture1.CurrentY = 140000 - 10000 * iPicture1.Print " " & 140000 - 10000 * i & " "Picture1.DrawStyle = 2Picture1.Line (0, 140000 - 10000 * i)-(360, 140000 - 10000 * i)NextFor j = 0 To 18Picture1.CurrentX = 360 - 20 * jPicture1.CurrentY = -5Picture1.Print 360 - 20 * jPicture1.DrawStyle = 2Picture1.Line (360 - 20 * j, 360)-(360 - 20 * j, -540)NextFor sita = 0 To 360 Step 0.01xb = lab * Cos(sita * pi / 180)yb = lab * Sin(sita * pi / 180)lbc = Sqr(lab ^ 2 + lac ^ 2 - 2 * lab * lac * Cos(sita * pi / 180))beta = arcsin(lab * Sin(sita * pi / 180) / lbc)xd = xb + lbd * Cos(beta)yd = yb - lbd * Sin(beta)gama = arccos((xd - 800) / Sqr((xd - 800) ^ 2 + yd ^ 2))If yd <= 0 Thengama1 = 2 * pi - gamaElsegama1 = gamaEnd Ifsita0 = sita + 0.001xb0 = lab * Cos(sita0 * pi / 180)yb0 = lab * Sin(sita0 * pi / 180)lbc0 = Sqr(lab ^ 2 + lac ^ 2 - 2 * lab * lac * Cos(sita0 * pi / 180)) beta0 = arcsin(lab * Sin(sita0 * pi / 180) / lbc0)xd0 = xb + lbd * Cos(beta0)yd0 = yb - lbd * Sin(beta0)gama0 = arccos((xd0 - 800) / Sqr((xd0 - 800) ^ 2 + yd0 ^ 2))If yd0 <= 0 Thengama10 = 2 * pi - gama0Elsegama10 = gama0End Ifw1 = (gama10 - gama1) / (pi / 180000)sita2 = sita + 0.01xb2 = lab * Cos(sita2 * pi / 180)yb2 = lab * Sin(sita2 * pi / 180)lbc2 = Sqr(lab ^ 2 + lac ^ 2 - 2 * lab * lac * Cos(sita2 * pi / 180)) beta2 = arcsin(lab * Sin(sita2 * pi / 180) / lbc2)xd2 = xb + lbd * Cos(beta2)yd2 = yb - lbd * Sin(beta2)gama2 = arccos((xd2 - 800) / Sqr((xd2 - 800) ^ 2 + yd2 ^ 2))If yd <= 0 Thengama12 = 2 * pi - gama2Elsegama12 = gamaEnd Ifsita02 = sita2 + 0.001xb02 = lab * Cos(sita02 * pi / 180)yb02 = lab * Sin(sita02 * pi / 180)lbc02 = Sqr(lab ^ 2 + lac ^ 2 - 2 * lab * lac * Cos(sita02 * pi / 180)) beta02 = arcsin(lab * Sin(sita02 * pi / 180) / lbc02)xd02 = xb02 + lbd * Cos(beta02)yd02 = yb02 - lbd * Sin(beta02)gama02 = arccos((xd02 - 800) / Sqr((xd02 - 800) ^ 2 + yd02 ^ 2)) If yd02 <= 0 Thengama102 = 2 * pi - gama02Elsegama102 = gama02End Ifw2 = (gama102 - gama12) / (pi / 180000)a = (w2 - w1) / 0.01Picture1.PSet (sita, a * 180 / pi)NextEnd SubPrivate Sub Form_Load()Timer1.Enabled = FalseEnd SubPrivate Sub Timer1_Timer()‘动画Dim k As IntegerDim a As IntegerDim beta As Doublelab = 200lac = 400lbd = 700lae = 800lce = 400pi = 3.1416Dim i, j As IntegerDim xb, yb, xd, yd As Double'画坐标系Picture1.BackColor = vbWhitePicture1.FillStyle = 0Picture1.Scale (-400, 440)-(1500, -440) Picture1.Line (-400, 0)-(1500, 0) Picture1.Line (0, -440)-(0, 440)For i = 0 To 22Picture1.CurrentX = -68Picture1.CurrentY = 440 - 40 * i Picture1.Print " " & 440 - 40 * i & " " Picture1.DrawStyle = 2Picture1.Line (-400, 440 - 40 * i)-(1500, 440 - 40 * i) NextFor j = 0 To 70Picture1.CurrentX = 1520 - 80 * jPicture1.CurrentY = -5Picture1.Print 1520 - 80 * jPicture1.DrawStyle = 2Picture1.Line (1520 - 80 * j, 440)-(1520 - 80 * j, -440) NextText1.Text = Val(Text1.Text) + 1a = Text1.Textlbc = Sqr(lab ^ 2 + lac ^ 2 - 2 * lab * lac * Cos(sita)) beta = arcsin(lab * Sin(sita) / lbc)Line1.X1 = 0Line1.Y1 = 0Line1.X2 = lab * Cos(a)Line1.Y2 = lab * Sin(a)Line2.X1 = lab * Cos(a)Line2.Y1 = lab * Sin(a)Line2.X2 = Line2.X1 + lbd * Cos(beta)Line2.Y2 = Line2.Y1 - lbd * Sin(beta)Line3.Y1 = Line2.Y2Line3.X1 = Line2.X2Line3.X2 = 800Line3.Y2 = 0Line1.Visible = TrueLine2.Visible = TrueLine3.Visible = TrueFor sita = 0 To 2 * pi Step 0.001xb = lab * Cos(sita)yb = lab * Sin(sita)lbc = Sqr(lab ^ 2 + lac ^ 2 - 2 * lab * lac * Cos(sita)) beta = arcsin(lab * Sin(sita) / lbc)xd = xb + lbd * Cos(beta)yd = yb - lbd * Sin(beta)Picture1.PSet (xd, yd)NextEnd Sub运行上述程序点击“运动”按钮得运动动画,截图如下:D点的轨迹如下:角位移如下:角速度如下:角加速度如下:运动结果分析D点做整周的回转运动,若以向右为角度的起始处,在D转角位180度时角速度最小,转角位0处角速度最大。
机械原理课程设计 连杆机构的设计及运动分析
目录一设计任务---------------------------------------------------2 二设计过程---------------------------------------------------22.1设计思想-----------------------------2 2.2参数的定义---------------------------2 2.3数学模型-----------------------------3 2.4程序流程图---------------------------5 2.5源程序设计--------------------------------------------6 三设计结果----------------------------------------------------123.1 连杆运动示意图----------------------123.2 连杆参数的计算结果------------------133.3 位移、角速度、加速度曲线绘制--------------14 四课程设计总结---------------------------------------15 五参考文献-----------------------------------------------15 六中期检查报告---------------------------------------17一设计任务任务:连杆机构的设计及运动分析已知:中心距X1=50mm,X2=160mm,Y=280mm,构件3的上、下极限ϕ/3=60、ϕ//3=120,滑块的冲程H=160mm,比值CE/CD=1/2,EF/DE=1/4,各构件S重心的位置,曲柄每分钟转速N1=90r/min。
要求:1)建立数学模型;2)用C语言编写计算程序并运行;3)绘制从动件运动规律线图,并进行连杆机构的动态显示;4)用计算机打印出计算说明。
哈尔滨工业大学机械原理大作业凸轮设计
大作业1连杆机构运动分析1.题目(8)如图所示机构,已知机构各构件的尺寸为==100AC CE l l mm ,==200BC CD l l mm ,90BCD ∠=︒,构件1的角速度为1=10/w rad s ,试求构件5的角位移、角速度和角加速度,并对计算结果进行分析。
1.1机构的运动分析AB 为原动件,AB 转动通过转动导杆机构带动杆BCD 转动,BCD 转动通过转动导杆机构带动杆DE 摆动。
1.2 机构的结构分析杆组可以划分为一个RR I级杆组(杆1)、RRPII级杆组(滑块2,杆3)、RPRII 级杆组(滑块4,杆5)(1)RRI级杆组1:(2)RRPII级杆组2,3:(3)RPRII级杆组4,5:2.分析过程 2.1 建立坐标系建立以点E 为原点的固定平面直角坐标系x-E-y ,如图所示。
2.2 建立数学模型(1)构件1、2、3的分析原动件杆1的转角:1θ=0--360。
原动件杆1的角速度:1ω=.1θ=10/rad s 原动件杆1的角加速度:..1αθ==0运动副A 的坐标:0200A A x y mm =⎫⎬=⎭运动副A 的速度及加速度都为零。
构件1为BC (RRP Ⅱ级杆组)上滑块B 的导路 滑块B 的位置为:132cos cos B A C x x s x l θθ=+=+132sin sin B A C y y s x l θθ=+=+消去s,得:0212arcsinA l θθ=+式中:011()sin ()cos C A C A A x x y y θθ=---构件3的角速度i ω和滑块B 沿导路的移动速度D υ:.211213(Q sin Q cos )/Q ωϕθθ==-+ 1322323(Q cos Q sin )/Q D s l l υθθ⋅==-+式中:..11111211321212Q sin ;Q cos ;Q sin sin cos sin l l l θθθθθθθθ=-==+构件3的角加速度和滑块B 沿导路移动的加速度:..241513(Q sin Q cos )/Q αθθθ==-+..4325323(Q cos Q sin )/Q B s l l υθθ==-+式中:122......21142211111Q cos sin cos 2sin l l l s θθθθθθθθ=---- 122......21152211111Q sin cos sin 2cos l l l s θθθθθθθθ=+-+(2)构件3,4,5的分析构件3,4,5,由1个Ⅰ级基本杆组和一个RRP Ⅱ级杆组组成,与构件1,2,3结构相同,只运动分析过程与其相反。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
机械原理大作业一 课程名称: 机械原理 设计题目: 连杆机构运动分析 院 系: 机电工程学院
指导教师: 林琳 刘福利 设计时间: 2014.06
哈尔滨工业大学 题目:如图所示机构,已知机构各构件的尺寸为齿轮1,齿轮2,齿轮5的齿数分为z1=17,z2=z5=40,BC=70mm,CD=270mm,EF=80mm,DE=240mm,DG=500mm,h=160mm,齿轮1的角速度为w1=10rad/s,试求点D的轨迹及构件7上点G的位移,速度和加速度,并对计算结果进行分析。
1.建立直角坐标系 上述机构可简化为下图中所示机构,其中BC和EF杆的角速度相等(大小相同,方向相同),所以该机构只有一个自由度。可建立如图所示的坐标系 2机构结构分析 该机构可认为由两个I级杆组RR(杆BC、杆EF)、II级杆组RRR(杆3、4)、II级杆组RRP(杆6及滑块G)组成。如下图所示:
I级杆组RR(杆BC、杆EF)
I
II级杆组RRR(杆3、4) II级杆组RRP(杆6及滑块G) 3.各基本杆组运动分析数学模型 3.1 I级杆组BC
如上图中所示,已知运动副C的位置BC=70mm, 𝒙𝑩=0 , 𝒚𝑩=0 , 速度𝒗𝒙𝑩=0 , 𝒗𝒚𝑩=𝟎 ,加速度𝒂𝒙𝑩=𝟎,𝒂𝒚𝑩
=𝟎.角速度𝛚=−𝟏𝟎×
𝟏𝟕÷𝟒𝟎 𝐫𝐚𝐝/𝐬,角加速度𝛂=𝟎. 初始转角𝝋𝟎=0. 转角符号为𝛗.
1)位置分析 cosBCBClxx sinBCBClyy
2)速度和加速度分析 sinBCxAxClvv sinBCyAyClvv sincos2BCBCxBxCllaa
sincos2BCBCyByCllaa
将已知参数带入便可求得 3.2 I级杆组EF
如上图中所示,已知运动副E的位置EF=80mm,可令BF=350mm从而,角加速度/40/1710srad𝛂=𝟎, 初始转角𝛗𝟎=0.转角符号为𝛗. 1)位置分析 cosEFFElxx
sinEFFElyy
yCxCyCxCCaavvyx,,,,,C.0,0加速度,0,0速度,0,350yFxFyFxFFFaavvymmx2)速度和加速度分析 sinEFxExElvv sinEFyFyElvv sincos2EFEFxFxEllaa
sincos2EFFEyFyEllaa
将已知参数带入便可求得 3.3 II级杆组RRR(杆3、4) 如下图中的示意图所示,已知参数 DC=270mm,DE=240.可求出运动副D点的一系列的 参数 以及两杆的角位置φ3,φ4,角速度ω3,𝜔4 , 和角加速度α3,𝛼4
.
1)位置方程 4433coscoslxlxxECD (1)
4433sinsinlylyyECD (2) 连立(1)、(2)求得: 𝜑3002020200arctan2CACBAB 其中: )(230CExxlA
yExEyExEEaavvyx,,,,,EyCxCyCxCCaavvyx,,,,,CyExEyExEEaavvyx,,,,,EyDxDyDxDDaavvyx,,,,,D )(230CEyylB 22242230)()(,CECEBDBDyyxxllllC
求得分𝜑3后,带入(1)、(2),便可求得𝜑4 , 2)速度方程
1443/)]()([GvvSvvCyCyExCxE
1434/)]()([GvvSvvCyCyExCxE 其中: 333333sin,coslSlC
444444sin,coslSlC 444333sinsinlvlvvxDxCxD
444333coscoslvlvvyDyCyD 3)加速度方程
143423/)(GSGCG
133324/)(GSGCG 其中: jjiixCxECCaaG222
jjiiyCyESSaaG223
323cossiniiiixCxDllaa
323sincosiiiiyCyDlflaa 3.4 II级杆组RRP(杆6、滑块G) 如下图中的示意图所示已知运动副D点的 杆长DG=0500jilmml,,滑块导路方向角f(j)=0和计算位移s的参考点K的参数
Dyx,DyDxDyDxDDaavvyx,,,,,D 可求G点的
1)位置方程 )(sin)(cos)(cosGjfljfsxiflxxjKiD )(cos)(sin)(sinGjfljfsyiflyyjKiD
可求得)(arcsin)(jfllAifijk,其中)(cos)()(sin)(0jfyyjfxxAKDKD带入上式即可解的xG,yG,s 2)速度方程 321/)(cos)(sinQjfQjfQwi(
321/)(sin)(cosQiflQiflQviiD( 其中: ))(cos)(sin(1jfljfsvvQjjxDxK ))(sin)(cos(2jfljfsvvQjjyDyK )(cos)(cos)(sin)(sin3jfifljfiflQii
yGxGyGxGGaavvyx,,,,,
G
jjyKxKykxkKaavvmmyx,,0,0,0,0,160,0K )(siniflvviixDxG
)(cosiflvviiyDyG 3)加速度方程
354/))(cos)(sin(QjfQjfQi
354/))(sin)(cos(QjflQiflQsiii 其中: )(sin2))(sin)(cos())(cos)(sin()(cos224jfvjfljfsjfljfsjflaaQjsjjjjiixDxK)(cos2))(sin)(sin())(sin)(cos()(sin225jfvjfljfsjfljfsjflaaQjsjjjjiiyDyK)(sin2)(cos)(sin)(cos2jfvjfsjfsjfaaajsjjsxKxG )(cos2)(sin)(cos)(sin2jfvjfsjfsjfaaajsjjsyKyG
四、计算编程 采用VB编程,其窗口代码如下: Private s7(3600) As Double Private v7(3600) As Double Private a7(3600) As Double Private x3(3600) As Double Private y3(3600) As Double Private Sub Form_Load() Dim i As Single Dim pi As Double Dim pa As Double Dim RR1 As RR Dim RR2 As RR Dim RRR1 As RRR Dim RRP1 As RRP Set RR1 = New RR Set RR2 = New RR Set RRR1 = New RRR Set RRP1 = New RRP pi = 3.1415926 pa = pi / 180
For i = 0 To 3600 RR1.xA = 0: RR1.yA = 0 RR1.vxA = 0: RR1.vyA = 0 RR1.axA = 0: RR1.ayA = 0 RR1.f = i * pa / 10: RR1.delt = 0 RR1.w = -10 * 17 / 40: RR1.e = 0 RR1.L = 70 RR1.cal
RR2.xA = 350: RR2.yA = 0 RR2.vxA = 0: RR2.vyA = 0 RR2.axA = 0: RR2.ayA = 0 RR2.f = i * pa / 10: RR2.delt = pi RR2.w = -10 * 17 / 40: RR2.e = 0 RR2.L = 80 RR2.cal
RRR1.xB = RR1.xB: RRR1.yB = RR1.yB RRR1.vxB = RR1.vxB: RRR1.vyB = RR1.vyB RRR1.axB = RR1.axB: RRR1.ayB = RR1.ayB RRR1.xD = RR2.xB: RRR1.yD = RR2.yB RRR1.vxD = RR2.vxB: RRR1.vyD = RR2.vyB RRR1.axD = RR2.axB: RRR1.ayD = RR2.ayB RRR1.Li = 270: RRR1.Lj = 240 RRR1.M = 1 RRR1.cal
x3(i) = RRR1.xC: y3(i) = RRR1.yC RRP1.xB = RRR1.xC: RRP1.yB = RRR1.yC RRP1.vxB = RRR1.vxC: RRP1.vyB = RRR1.vyC RRP1.axB = RRR1.axC: RRP1.ayB = RRR1.ayC RRP1.Li = 500: RRP1.Lj = 0