哈工大机械原理大作业一12题
哈尔滨工业大学机械原理大作业_连杆26题

1.运动分析题目如图所示机构,已知机构各构件的尺寸为AB=108mm,EF=320mm,BC=CE=CD=200mm,FG=162mm,AD=258mm,AG=514mm,DG=384mm,β=80º,构件1的角速度为ω1=10rad/s,试求构件2上点E的轨迹及构件5的角位移、角速度和角加速度,并对计算结果进行分析。
2.机构的结构分析,组成机构的基本杆组划分3.各基本杆组的运动分析数学模型(1)RR基本杆组:delt=0xB = xA + AB * Cos(f + delt)yB = yA + AB* Sin(f + delt)vxB = vxA - w * AB * Sin(f + delt)vyB = vyA + w * AB* Cos(f + delt)axB = axA - w ^ 2 * AB * Cos(f + delt):ayB = ayA - w ^ 2 * AB* Sin(f + delt)(2)RRR基本杆组Ci = lbc * Cos(fbc)Si = lbc * Sin(fbc)Cj = lcd * Cos(fcd)Sj = lcd * Sin(fcd)G1 = Ci * Sj - Cj * Siwbc = (Cj * (vxD - vxB) + Sj * (vyD - vyB)) / G1wcd = (Ci * (vxD - vxB) + Si * (vyD - vyB)) / G1vxC = vxB - wbc * lbc * Sin(fbc)vyC = vyB + wbc * lbc * Cos(fbc)G2 = axD - axB + wbc ^ 2 * Ci - wcd ^ 2 * CjG3 = ayD - ayB + wbc ^ 2 * Si - wcd ^ 2 * Sjebc = (G2 * Cj + G3 * Sj) / G1ecd = (G2 * Ci + G3 * Si) / G1axC = axB - ebc * lbc * Sin(fbc) - wbc ^ 2 * lbc * Cos(fbc)ayC = ayB + ebc * lbc * Cos(fbc) - wbc ^ 2 * lbc * Sin(fbc)EB = 2 * BC * Cos(febc)xE = xB + EB* Cos(fbc + febc)yE = yB + EB * Sin(fbc + febc)vxE = vxB – wbc * EB* Sin(fbc + febc)vyE = vyB + wbc * EB* Cos(fbc + febc)axE = axB - wbc ^ 2 * EB * Cos(fbc + delt) - ebc * EB * Sin(fbc + febc)ayE = ayB - wbc ^ 2 * leb * Sin(fbc + delt) + ebc * leb * Sin(fbc + febc) Ci = lef * Cos(fef)Si = lef * Sin(fef)Cj = lfg * Cos(ffg)Sj = lfg * Sin(ffg)G1 = Ci * Sj - Cj * Siwef = (Cj * (vxG - vxE) + Sj * (vyG - vyE)) / G1wfg = (Ci * (vxG - vxE) + Si * (vyG - vyE)) / G1vxF = vxE - wef * lef * Sin(fef)vyF = vyE + wef * lef * Cos(fef)G2 = axG - axE + wef ^ 2 * Ci - wfg ^ 2 * CjG3 = ayG - ayE + wef ^ 2 * Si - wfg ^ 2 * Sjeef = (G2 * Cj + G3 * Sj) / G1efg = (G2 * Ci + G3 * Si) / G1axF = axE - eef * lef * Sin(fef) - wef ^ 2 * lef * Cos(fef)ayF = ayE + eef * lef * Cos(fef) - wef ^ 2 * lef * Sin(fef)4.计算编程Dim xA As DoubleDim yA As DoubleDim vxA As DoubleDim vyA As DoubleDim axA As Double 'A '点加速度x轴分量Dim ayA As Double 'A '点加速度y轴分量Dim xB As Double 'B'点'x轴坐标Dim yB As Double 'B点y轴坐标Dim vxB As Double 'B点速度x轴分量Dim vyB As Double 'B点速度y轴分量Dim axB As Double 'B点加速度x轴分量Dim ayB As Double 'B点加速度y轴分量Dim xC As Double 'C点x轴坐标Dim yC As Double C'点y轴坐标Dim vxC As Double 'C点速度x轴分量Dim vyC As Double 'C点速度y轴分量Dim axC As Double 'C点加速度x轴分量Dim ayC As Double 'C点加速度y轴分量Dim xD As Double 'D点x轴坐标Dim yD As Double 'D点y轴坐标Dim vxD As Double 'D点速度x轴分量Dim vyD As Double 'D点速度y轴分量Dim axD As Double 'D点加速度x轴分量Dim ayD As Double 'D点加速度y轴分量Dim xE As Double 'E点x轴坐标Dim yE As Double 'E点y轴坐标Dim vxE As Double 'E点速度x轴分量Dim vyE As Double 'E点速度y轴分量Dim axE As Double 'E点加速度x轴分量Dim ayE As Double 'E点加速度y轴分量Dim xF As Double 'F点x轴坐标Dim yF As Double 'F点y轴坐标Dim vxF As Double 'F点速度x轴分量Dim vyF As Double 'F点速度y轴分量Dim axF As Double 'F点加速度x轴分量Dim ayF As Double 'F点加速度y轴分量Dim xG As Double 'G点x轴坐标Dim yG As Double 'G点y轴坐标Dim vxG As Double 'G点速度x轴分量Dim vyG As Double 'G点速度y轴分量Dim axG As Double 'G点加速度x轴分量Dim ayG As Double 'G点加速度y轴分量Dim delt As Double ' AB杆初始转角Dim lab As Double 'AB杆长Dim lbc As Double 'BC杆长Dim lcd As Double ' CD杆长Dim lce As Double 'CE杆长Dim lef As Double 'EF杆长Dim lfg As Double 'FG杆长Dim leb As Double 'ED杆长Dim f As Double 'AB杆转角Dim fbc As Double 'BC杆转角Dim fcd As Double 'CD杆转角Dim fce As Double 'CE杆转角Dim fef As Double 'EF杆转角Dim ffg As Double 'FG杆转角Dim fge As Double 'ge杆转角Dim w As Double 'AB杆角速度Dim wbc As Double ' BC角速度Dim wcd As Double 'CD角速度Dim wce As Double 'CE角速度Dim wef As Double 'EF角速度Dim wfg As Double 'FG角速度Dim e As Double 'AB杆角加速度Dim ebc As Double ' BC杆角加速度Dim ecd As Double 'CD杆角加速度Dim ece As Double 'CE杆角加速度Dim eef As Double 'EF杆角加速度Dim efg As Double 'FG杆角加速度Dim LBD As Double 'BD距离Dim leg As Double 'EG距离Dim JCBD As Double '角CBDDim jfeg As Double '角FEGDim fBD As Double 'BD转角Dim feg As Double 'EG转角Dim Ci As DoubleDim Cj As DoubleDim Si As DoubleDim Sj As DoubleDim G1 As DoubleDim G2 As DoubleDim G3 As DoubleDim val As DoubleDim pi As DoubleDim pa As DoubleDim febc As Double '角EBCDim i As DoubleDim fj1 As DoublePrivate Sub Command1_Click() '求点E的轨迹Picture1.Scale (-300, 400)-(10, -15)Picture1.Line (-300, 0)-(10, 0) 'XPicture1.Line (0, 400)-(0, -15) 'YFor i = -300 To 0 Step 50 'X轴坐标Picture1.DrawStyle = 2Picture1.Line (i, 400)-(i, 0)Picture1.CurrentX = i - 10: Picture1.CurrentY = 0 Picture1.Print iNext iFor i = 0 To 350 Step 50 'Y轴坐标Picture1.DrawStyle = 2Picture1.Line (0, i)-(-400, i)Picture1.CurrentX = -20: Picture1.CurrentY = i + 7 Picture1.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall RR1Call RRR1Call RR2Picture1.PSet (xE, yE)Next fj1End SubPrivate Sub Command2_Click() '求构件5的角位移Picture2.Scale (-20, 5)-(380, -0.5)Picture2.Line (-20, 0)-(380, 0) 'XPicture2.Line (0, 3)-(0, -0.5) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture2.DrawStyle = 2Picture2.Line (i, 3)-(i, 0)Picture2.CurrentX = i - 10: Picture2.CurrentY = 0 Picture2.Print iNext iFor i = -0.5 To 3 Step 0.5 'Y轴坐标Picture2.Line (0, i)-(380, i)Picture2.CurrentX = -25: Picture2.CurrentY = i Picture2.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Picture2.PSet (fj1, ffg)Next fj1End SubPrivate Sub Command3_Click() '求构件5的角速度Picture3.Scale (-20, 10)-(380, -10)Picture3.Line (-20, 0)-(380, 0) 'XPicture3.Line (0, 10)-(0, -10) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture3.DrawStyle = 2Picture3.Line (i, 10)-(i, -10)Picture3.CurrentX = i - 10: Picture3.CurrentY = 0 Picture3.Print iNext iFor i = -8 To 8 Step 2 'Y轴坐标Picture3.Line (0, i)-(380, i)Picture3.CurrentX = -20: Picture3.CurrentY = i Picture3.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Picture3.PSet (fj1, wfg)Next fj1End SubPrivate Sub Command4_Click() '求构件5的角加速度Picture4.Scale (-20, 300)-(380, -200)Picture4.Line (-20, 0)-(380, 0) 'XPicture4.Line (0, 300)-(0, -200) 'YFor i = 0 To 360 Step 30 'X轴坐标Picture4.DrawStyle = 2Picture4.Line (i, 300)-(i, -200)Picture4.CurrentX = i - 10: Picture4.CurrentY = 0 Picture4.Print iNext iFor i = -200 To 300 Step 50 'Y轴坐标Picture4.Line (0, i)-(380, i)Picture4.CurrentX = -25: Picture4.CurrentY = i + 5 Picture4.Print iNext iFor fj1 = 0 To 360 Step 0.01f = fj1 * paCall RR1Call RRR1Call RR2Call RRR2Picture4.PSet (fj1, efgNext fj1End SubPrivate Sub Form_Load() '赋初值lab = 108lce = 200lbc = 200lcd = 200lef = 320lfg = 162w = 10e = 0delt = 0xA = 0yA = 0vyA = 0axA = 0ayA = 0xD = -178.311284yD = 186.464704vxD = 0vyD = 0axD = 0ayD = 0xG = -514yG = 0vxG = 0vyG = 0axG = 0ayG = 0pi = 3.1415926pa = pi / 180febc = pa * 50End SubPrivate Sub RR1() 'RR基本杆组xB = xA + lab * Cos(f + delt)yB = yA + lab * Sin(f + delt)vxB = vxA - w * lab * Sin(f + delt)vyB = vyA + w * lab * Cos(f + delt)axB = axA - w ^ 2 * lab * Cos(f + delt) - e * lab * Sin(f + delt)ayB = ayA - w ^ 2 * lab * Sin(f + delt) + e * lab * Sin(f + delt)End SubPrivate Sub RR2() 'RR基本杆组leb = 2 * lbc * Cos(febc)xE = xB + leb * Cos(fbc + febc)yE = yB + leb * Sin(fbc + febc)vxE = vxB - wbc * leb * Sin(fbc + febc)vyE = vyB + wbc * leb * Cos(fbc + febc)axE = axB - wbc ^ 2 * leb * Cos(fbc + delt) - ebc * leb * Sin(fbc + febc) ayE = ayB - wbc ^ 2 * leb * Sin(fbc + delt) + ebc * leb * Sin(fbc + febc) End SubPrivate Sub RRR1() 'RRR基本杆组LBD = Sqr((xD - xB) ^ 2 + (yD - yB) ^ 2)If LBD > lbc + lcd And LBD < Abs(lbc - lcd) ThenIf MsgBox("RRR杆组杆长不符合要求", vbOKOnly, "提示") = 1 Then EndEnd IfElseEnd IfIf LBD < lbc + lcd And LBD > Abs(lbc - lcd) Then val = (lbc ^ 2 + LBD ^ 2 - lcd ^ 2) / (2 * lbc * LBD) JCBD = Atn(-val / Sqr(-val * val + 1)) + 2 * Atn(1) ElseEnd IfIf LBD = lbc + lcd ThenJCBD = 0ElseEnd IfIf LBD = Abs(lbc - lcd) ThenIf lbc > lcd ThenJCBD = 0ElseEnd IfIf lbc < lcd ThenJCBD = piElseEnd IfElseEnd IfIf xD > xB And yD >= yB Then '第一象限fBD = Atn((yD - yB) / (xD - xB))ElseEnd IfIf xD = xB And yD > yB ThenfBD = pi / 2ElseEnd IfIf xD < xB And yD >= yB Then '第二象限fBD = pi + Atn((yD - yB) / (xD - xB))ElseEnd IfIf xD < xB And yD < yB Then '第三象限fBD = pi + Atn((yD - yB) / (xD - xB))ElseEnd IfIf xD = xB And yD < yB ThenfBD = 3 * pi / 2ElseEnd IfIf xD > xB And yD <= yB Then '第四象限fBD = 2 * pi + Atn((yD - yB) / (xD - xB))ElseEnd Iffbc = fBD - JCBDxC = xB + lbc * Cos(fbc)yC = yB + lbc * Sin(fbc)If xC > xD And yC >= yD Then '第一象限fcd = Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC = xD And yC >= yD Thenfcd = pi / 2ElseEnd IfIf xC < xD And yC >= yD Then '第二象限fcd = pi + Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC < xD And yC < yD Then '第三象限fcd = pi + Atn((yC - yD) / (xC - xD))ElseEnd IfIf xC = xD And yC < yD Thenfcd = 3 * pi / 2ElseEnd IfIf xC > xD And yC <= yD Then '第四象限fcd = 2 * pi + Atn((yC - yD) / (xC - xD))ElseEnd IfCi = lbc * Cos(fbc)Si = lbc * Sin(fbc)Cj = lcd * Cos(fcd)Sj = lcd * Sin(fcd)G1 = Ci * Sj - Cj * Siwbc = (Cj * (vxD - vxB) + Sj * (vyD - vyB)) / G1 wcd = (Ci * (vxD - vxB) + Si * (vyD - vyB)) / G1 vxC = vxB - wbc * lbc * Sin(fbc)vyC = vyB + wbc * lbc * Cos(fbc)G2 = axD - axB + wbc ^ 2 * Ci - wcd ^ 2 * CjG3 = ayD - ayB + wbc ^ 2 * Si - wcd ^ 2 * Sj ebc = (G2 * Cj + G3 * Sj) / G1ecd = (G2 * Ci + G3 * Si) / G1axC = axB - ebc * lbc * Sin(fbc) - wbc ^ 2 * lbc * Cos(fbc)ayC = ayB + ebc * lbc * Cos(fbc) - wbc ^ 2 * lbc * Sin(fbc)End SubPrivate Sub RRR2() 'RRR基本杆组leg = Sqr((xG - xE) ^ 2 + (yG - yE) ^ 2)If leg > lef + lfg And leg < Abs(lef - lfg) ThenIf MsgBox("RRR杆组杆长不符合要求", vbOKOnly, "提示") = 1 Then EndElseEnd IfElseEnd IfIf leg < lef + lfg And leg > Abs(lef - lfg) Thenval = (lef ^ 2 + leg ^ 2 - lfg ^ 2) / (2 * lef * leg)jfeg = Atn(-val / Sqr(-val * val + 1)) + 2 * Atn(1)ElseEnd IfIf leg = lef + lfg Thenjfeg = 0ElseEnd IfIf leg = Abs(lef - lfg) ThenIf lef > lfg Thenjfeg = 0ElseEnd IfIf lef < lfg Thenjfeg = piElseEnd IfElseEnd IfIf xG > xE And yG >= yE Then '第一象限feg = Atn((yG - yE) / (xG - xE))ElseEnd IfIf xG = xE And yG > yE Thenfeg = pi / 2ElseEnd IfIf xG < xE And yG >= yE Then '第二象限feg = pi + Atn((yG - yE) / (xG - xE))ElseEnd IfIf xG < xE And yG < yE Then '第三象限feg = pi + Atn((yG - yE) / (xG - xE)) ElseEnd IfIf xG = xE And yG < yE Thenfeg = 3 * pi / 2ElseEnd IfIf xG > xE And yG <= yE Then '第四象限feg = 2 * pi + Atn((yG - yE) / (xG - xE)) ElseEnd Iffef = feg - jfegxF = xE + lef * Cos(fef)yF = yE + lef * Sin(fef)If xF > xG And yF >= yG Then '第一象限ffg = Atn((yF - yG) / (xF - xG))ElseEnd IfIf xF = xG And yF >= yG Thenffg = pi / 2ElseEnd IfIf xF < xG And yF >= yG Then '第二象限ffg = pi + Atn((yF - yG) / (xF - xG)) ElseEnd IfIf xF < xG And yF < yG Then '第三象限ffg = pi + Atn((yF - yG) / (xF - xG)) ElseEnd IfIf xF = xG And yF < yG Thenffg = 3 * pi / 2ElseEnd IfIf xF > xG And yF <= yG Then '第四象限ffg = 2 * pi + Atn((yF - yG) / (xF - xG)) ElseEnd IfCi = lef * Cos(fef)Si = lef * Sin(fef)Cj = lfg * Cos(ffg)Sj = lfg * Sin(ffg)G1 = Ci * Sj - Cj * Siwef = (Cj * (vxG - vxE) + Sj * (vyG - vyE)) / G1wfg = (Ci * (vxG - vxE) + Si * (vyG - vyE)) / G1vxF = vxE - wef * lef * Sin(fef)vyF = vyE + wef * lef * Cos(fef)G2 = axG - axE + wef ^ 2 * Ci - wfg ^ 2 * CjG3 = ayG - ayE + wef ^ 2 * Si - wfg ^ 2 * Sjeef = (G2 * Cj + G3 * Sj) / G1efg = (G2 * Ci + G3 * Si) / G1axF = axE - eef * lef * Sin(fef) - wef ^ 2 * lef * Cos(fef)ayF = ayE + eef * lef * Cos(fef) - wef ^ 2 * lef * Sin(fef)End Sub5.结果及分析图1 E点的运动轨迹(1)由图1所示,E点的运动轨迹呈稍倾斜“8字形”。
哈工大机械原理大作业

H a r b i n I n s t i t u t e o f T e c h n o l o g y机械原理大作业一课程名称:机械原理设计题目:连杆机构运动分析院系:机电学院班级: 1208105分析者:殷琪学号:指导教师:丁刚设计时间:哈尔滨工业大学设计说明书1 、题目如图所示机构,一只机构各构件的尺寸为AB=100mm,BC=4.28AB,CE=4.86AB,BE=8.4AB,CD=2.14AB,AD=4.55AB,AF=7AB,DF=3.32AB,∠BCE=139?。
构件1的角速度为ω1=10rad/s,试求构件2上点E的轨迹及构件5的角位移、角速度和角加速度,并对计算结果进行分析。
2、机构结构分析该机构由6个构件组成,4和5之间通过移动副连接,其他各构件之间通过转动副连接,主动件为杆1,杆2、3、4、5为从动件,2和3组成Ⅱ级RRR基本杆组,4和5组成Ⅱ级RPR 基本杆组。
如图建立坐标系3、各基本杆组的运动分析数学模型1) 位置分析2) 速度和加速度分析 将上式对时间t 求导,可得速度方程:将上式对时间t 求导,可得加速度方程:RRR Ⅱ级杆组的运动分析如下图所示 当已知RRR 杆组中两杆长L BC 、L CD 和两外副B 、D 的位置和运动时,求内副C的位置、两杆的角位置、角运动以及E 点的运动。
1) 位置方程由移项消去j ϕ后可求得i ϕ:式中,可求得j ϕ:E 点坐标方程:其中2) 速度方程两杆角速度方程为式中,点E 速度方程为3) 加速度方程两杆角加速度为式中,点E 加速度方程为RPR Ⅱ级杆组的运动分析(1) 位移方程(2)速度方程其中(3)加速度方程4、 计算编程利用MATLAB 软件进行编程,程序如下:% 点B 和AB 杆运动状态分析>>r=pi/180;w 1=10;e 1=0;l 1=100;Xa=0;Ya=0;Vax=0;Vay=0;aax=0;aay=0;f1=0:1: 360;% B 点位置Xb=Xa+l1*cos(r*f1);Yb=Ya+l1*sin(r*f1);% B点速度Vbx=Vax-w1*l1*sin(r*f1);Vby=Vay+w1*l1*cos(r*f1);% B点加速度abx=aax-l1*w1.^2.*cos(r*f1);aby=aay-l1*w1.^2.*sin(r*f1);% RRR2级杆组运动分析% 输入D点参数l2=428;l3=214;Xd=455;Yd=0;Vdx=0;Vdy=0;adx=0;ady=0;% 计算E点、2杆、3杆运动参数lbe=840;lce=486;a0=2*l2*(Xd-Xb);b0=2*l2*(Yd-Yb);c0=l2^2+(Xb-Xd).^2+(Yb-Yd).^2-l3^2;f2=2*atan((b0+sqrt(a0.^2+b0.^2-c0.^2))./(a0+c0)); % C点位置Xc=Xb+l2*cos(f2);Yc=Yb+l2*sin(f2);% 2杆、3杆运动参数计算dX=Xc-Xd;dY=Yc-Yd;for n=1:length(dX)if dX(n)>0&dY(n)>=0f3(n)=atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)>0f3(n)=pi/2;elseif dX(n)<0&dY(n)>=0f3(n)=pi+atan(dY(n)/dX(n));elseif dX(n)<0&dY(n)<0f3(n)=pi+atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)<0f3(n)=1.5*pi;elseif dX(n)>0&dY(n)<0f3(n)=2*pi+atan(dY(n)/dX(n));endendC2=l2*cos(f2);C3=l3*cos(f3);S2=l2*sin(f2);S3=l3*sin(f3);G1=C2.*S3-C3.*S2;w2=(C3.*(Vdx-Vbx)+S3.*(Vdy-Vby))./G1;w3=(C2.*(Vdx-Vbx)+S2.*(Vdy-Vby))./G1;G2=adx-abx+(w2.^2).*C2-(w3.^2).*C3;G3=ady-aby+(w2.^2).*S2-(w3.^2).*S3;e2=(G2.*C3+G3.*S3)./G1;% E点位置w=acos((l2^2+lbe^2-lce^2)/(2*l2*lbe));Xe=Xb+lbe*cos(f2-w);Ye=Yb+lbe*sin(f2-w);Vex=Vbx-lbe*w2.*sin(f2-w);Vey=Vby+lbe*w2.*cos(f2-w);aex=abx-lbe*(e2.*sin(f2-w)+w2.^2.*cos(f2-w));aey=aby+lbe*(e2.*cos(f2-w)-w2.^2.*sin(f2-w));% 计算杆5运动参数Xf=646.2912088;Yf=-268.9008617;l5=sqrt((Xe-Xf).^2+(Ye-Yf).^2);dX=Xe-Xf;dY=Ye-Yf;for n=1:length(dX)if dX(n)>0&dY(n)>=0f5(n)=atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)>0f5(n)=pi/2;elseif dX(n)<0&dY(n)>=0f5(n)=pi+atan(dY(n)/dX(n));elseif dX(n)<0&dY(n)<0f5(n)=pi+atan(dY(n)/dX(n));elseif dX(n)==0&dY(n)<0f5(n)=1.5*pi;elseif dX(n)>0&dY(n)<0f5(n)=2*pi+atan(dY(n)/dX(n));endendw5=(-Vex.*sin(f5)+Vey.*cos(f5))./l5;a5=(-aex.*sin(f5)+aey.*cos(f5))./l5;% 画出各参数曲线figure(1);plot(Xe,Ye,'k');xlabel('Xe/\mm');ylabel('Ye/mm');grid on;title('E点位置');figure(2);plot(f1,f5,'k');xlabel('f/\circ');ylabel('f5/\circ');grid on;title('5杆角位移');figure(3);plot(f1,w5,'k');xlabel('f/\circ');ylabel('w5/rad/s');grid on;title('5杆角速度');figure(4);plot(f1,a5,'k');xlabel('f/\circ');ylabel('a5/rad/s2');gridon;title('5杆角加速度');Warning: Unable to interpret TeX string "Xe/\mm"5、计算结果图一:E点的运动轨迹图二:5杆角位移图三:5杆角速度图四:5杆角加速度6、计算结果分析由E点位置图像可看出,构件4做周期往复运动,由图二、三、四可看出,构件5的角位移、角速度、角加速度均成周期性变化。
哈工大机械原理大作业直动从动件盘形凸轮机构满分完美版哈尔滨工业大学

end
%曲率半径
dx=diff(x);
dx(36001)=0;
dy=diff(y);
dy(36001)=0;
dydx=dy./dx;
ddy=diff(dydx);
ddy(36001)=0;
ddy=ddy./dx;
~7~
xlabel('凸轮转角φ/°');
ylabel('从动件加速度 a/(mm/s^2)');
title('加速度');
subplot(2,3,4);
plot(dsdPhi,s);
axis equal
axis([-30 50 -50 30]);
hold on
plot(x1,f1,'r');
hold on
axis equal
hold on
plot(x,y,'k')
legend('理论廓线',-1);
hold on
%轮廓图
%实际廓线
%理论廓线
~8~
哈尔滨工业大学
直动从动件盘型凸轮机构设计说明书
plot(x0,y0,':')
hold on
%基圆
plot(xe,ye,'k:')
%偏距圆
legend('实际廓线','理论廓线','基圆','偏距圆',-1);
subplot(2,3,1);
plot(Phi,s);
grid on
axis([0 360 0 27]);
哈工大-机械原理大作业-连杆机构运动分析

机械原理大作业(一)作业名称:连杆机构运动分析设计题目: 20院系:英才学院班级: XXXXXXX设计者:邵广斌学号: XXXXXXXXXX指导教师:林琳设计时间: 2013年05月19日哈尔滨工业大学机械设计1.运动分析题目如图所示机构,已知机构各构件的尺寸为150AB mm =,97β=︒,400BC mm =,300CD mm =,320AD mm =,100BE mm =,230EF mm =,400FG mm =,构件1的角速度为110/rad s ω=,试求构件2上点F 的轨迹及构件5上点G 的位移、速度和加速度,并对计算结果进行分析。
2. 机构分析该机构由原动件AB (Ⅰ级杆组)、BCD (RRR Ⅱ级杆组)和FG (RRP Ⅱ级杆组)组成。
3. 建立坐标系如图3,建立以定点A 为原点的平面直角坐标系A-xy 。
图1 运动机构结构图4. 运动分析数学模型4.1 原动件AB原动件AB 的转角: 10~2ψπ= 原动件AB 的角速度:110/rad s ω=原动件AB 的角加速度: 10α= 运动副A 的位置坐标: 0A x = 0A y =运动副A 的速度: 0xA v = 0yA v = 运动副A 的加速度: 0xA a = 0yA a =原动件AB 长度:150AB l mm =运动副B 的位置坐标: 1B A AB x x l cos ψ=+1B A AB y x l sin ψ=+运动副B 的速度: 11 xB xA AB v v l sin ωψ=-11 yB yA AB v v l cos ωψ=+运动副B 的加速度: 2 1111 xBxA AB AB a a l cos l sin ωψαψ=--21111yB yA AB AB a a l sin l cos ωψαψ=-+4.2 RRR Ⅱ级杆组BCD运动副D 的位置坐标: 320D x mm = 0D y = 运动副D 的速度: 0xD v = 0yD v = 运动副D 的加速度: 0xD a = 0yD a = 杆BC 长度: 400BC l mm = 杆CD 长度:300CD l mm =BC 相对于x 轴转角:200ψ=其中02BC D B A l x x =-() 0 2 BC D B B l y y =-()2220B B C C l C l D l D =+- 222())(BDD B D B l x x y y =-+- CD 相对于x 轴转角: 3C DC Dy y arctanx x ψ-=-求导可得BC 角速度2ω、角加速度2α以及CD 角速度3ω、角加速度3α。
哈工大机械原理大作业_凸轮机构运动分析12题

Harbin Institute of Technology机械原理大作业二课程名称:机械原理设计题目:凸轮结构设计设计时间:2013年6月25日哈尔滨工业大学1、设计题目如图所示直动从动件盘形凸轮机构,原始参数见下表,据此设计该凸轮机构。
2、推杆升程、回程运动方程对于不同运动规律的凸轮结构,其上升与下降的方式不一,但遵循同样的运动顺序:上升、远休止点恒定,下降、近休止点恒定。
因此,在设计时,仅需确定这四个阶段的角度与位置即可。
(1) 推杆升程运动方程式中(2) 推杆远休程运动方程在远休程段,即时,。
(3) 推杆回程运动方程式中 (4) 推杆近休程运动方程在远休程段,即时,。
3、推杆位移、速度、加速度线图(为方便作图和坐标的度量,取,用Matlab 作图)(1) 推杆位移线图升程 (mm )升程 运动角(°) 升程运动 规律 升程许用压力角(°) 回程 运动角(°) 回程运动 规律 回程许用压力角(°) 远休 止角(°) 近休 止角(°) 801503-4-5 多项式401003-4-5 多项式604070程序:x1=linspace(0,5*pi/6,1000);x2=linspace(5*pi/6,19*pi/18,1000);x3=linspace(19*pi/18,29*pi/18,1000);x4=linspace(29*pi/18,2*pi,1000);T1=x1/(5*pi/6);T2=(x3-19*pi/18)/(5*pi/9);s1=80*(10*T1.^3-15*T1.^4+6*T1.^5);s2=80;s3=80*(1-(10*T2.^3-15*T2.^4+6*T2.^5));s4=0;plot(x1,s1,'k',x2,s2,'k',x3,s3,'k',x4,s4,'k') xlabel('角度/ψrad');ylabel('位移s/mm');title('推杆位移线图');grid(2)推杆速度线图程序:x1=linspace(0,5*pi/6,1000);x2=linspace(5*pi/6,19*pi/18,1000);x3=linspace(19*pi/18,29*pi/18,1000);x4=linspace(29*pi/18,2*pi,1000);T1=x1/(5*pi/6);T2=(x3-19*pi/18)/(5*pi/9);v1=30*80*1*T1.^2.*(1-2*T1+T1.^2)/(5*pi/6);v2=0;v3=-30*80*1*T2.^2.*(1-2*T2+T2.^2)/(5*pi/9);v4=0;plot(x1,v1,'k',x2,v2,'k',x3,v3,'k',x4,v4,'k') xlabel('角度ψ/rad');ylabel('速度v/(mm/s)')title('推杆速度线图')grid(3)推杆加速度线图程序:x1=linspace(0,5*pi/6,1000);x2=linspace(5*pi/6,19*pi/18,1000);x3=linspace(19*pi/18,29*pi/18,1000);x4=linspace(29*pi/18,2*pi,1000);T1=x1/(5*pi/6);T2=(x3-19*pi/18)/(5*pi/9);a1=60*80*1*1*T1.*(1-3*T1+2*T1.^2)/((5*pi/6)^2);a2=0;a3=-60*80*1*1*T2.*(1-3*T2+2*T2.^2)/((5*pi/9)^2);a4=0;plot(x1,a1,'k',x2,a2,'k',x3,a3,'k',x4,a4,'k')xlabel('角度ψ/rad');ylabel('加速度a/ ')title('推杆加速度线图')grid4、凸轮机构的线图程序:x1=linspace(0,5*pi/6,1000);x2=linspace(5*pi/6,19*pi/18,1000);x3=linspace(19*pi/18,29*pi/18,1000);x4=linspace(29*pi/18,2*pi,1000);T1=x1/(5*pi/6);T2=(x3-19*pi/18)/(5*pi/9);s1=80*(10*T1.^3-15*T1.^4+6*T1.^5);s2=80;s3=80*(1-(10*T2.^3-15*T2.^4+6*T2.^5));s4=0;v1=30*80*1*T1.^2.*(1-2*T1+T1.^2)/(5*pi/6);v2=0;v3=-30*80*1*T2.^2.*(1-2*T2+T2.^2)/(5*pi/9);v4=0;plot(v1,s1,'r',v2,s2,'r',v3,s3,'r',v4,s4,'r')xlabel('ds/dψ');ylabel('(位移s/mm)')title('ds/dψ—s曲线')grid5、凸轮的基圆半径和偏距以凸轮机构的ds/dφ-s线图为基础,可分别作出三条限制线(推程许用压力角的切界线,回程许用压力角的限制线,起始点压力角许用线),以这三条线可确定最小基圆半径及所对应的偏距e,在其下方选择一合适点,即可满足压力角的限制条件。
哈工大机械原理大作业一连杆运动分析(02)

哈⼯⼤机械原理⼤作业⼀连杆运动分析(02)⼀.设计题⽬⼆. 结构分析与基本杆组划分1.机构的结构分析机构各构件都在同⼀平⾯内运动,活动构件数n=3 P L=4 P H=0则机构的⾃由度为: F = 3n -2P L –P H = 3×3-2×4 = 12.基本杆组划分(1)去除虚约束和局部⾃由度本机构中⽆虚约束或局部⾃由度,此步骤跳过。
(2)拆杆组。
从远离原动件(即杆1)进⾏拆分,就可以得到由杆2,3组成的RRRⅡ级杆组和Ⅰ级机构杆1。
如下图:(3)确定机构的级别由(2)知,机构为Ⅱ级机构。
三. 运动分析数学模型以A为原点建⽴坐标系,如图:原动件AB的转⾓:φ1=0--2π运动副A的位置坐标:x A=0 y A=0 运动副D的位置坐标:x D=d y D=0 则运动副B的位置坐标:x B = acosφ1 y B = asinφ1其中:t=0:0.001:2*pi;a=60;b=90;c=120;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd)曲柄a=50,55,60,65红蓝绿黄b=90,c=120,d=100 t=0:0.001:2*pi; a=50;b=90;c=120;d=100;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'r'); hold on;grid on;t=0:0.001:2*pi;d=100;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'b'); hold on;t=0:0.001:2*pi;a=60;b=90;c=120;d=100;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'g'); hold on;t=0:0.001:2*pi;a=65;b=90;c=120;d=100;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'y');a=50,55,60,65红蓝绿黄2.摇杆c=105,115,125,135红蓝绿黄a=50,b=90,d=100 t=0:0.001:2*pi; a=50;b=90;d=100;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'r'); hold on;grid on;t=0:0.001:2*pi;a=50;b=90;c=115;d=100;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'b');t=0:0.001:2*pi;a=50;b=90;c=125;d=100;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'g'); hold on;t=0:0.001:2*pi;a=50;b=90;d=100;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'y');c=105,115,125,135红蓝绿黄3.连杆b=80,90,100,110红蓝绿黄a=50,c=120,d=100 t=0:0.001:2*pi; a=50;b=80;c=120;d=100;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'r'); hold on;grid on;t=0:0.001:2*pi;a=50;b=90;c=120;d=100;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'b'); hold on;t=0:0.001:2*pi;a=50;b=100;c=120;d=100;xa=0;ya=0;xd=d;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'g'); hold on;t=0:0.001:2*pi;a=50;b=110;c=120;d=100;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'y');b=80,90,100,110红蓝绿黄4.机架d=85,95,105,115 红蓝绿黄a=50,b=90,c=120 t=0:0.001:2*pi; a=50;b=90;c=120;d=85;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'r'); hold on;grid on;t=0:0.001:2*pi;a=50;b=90;c=120;d=95;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'b'); hold on;t=0:0.001:2*pi;d=105;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'g'); hold on;t=0:0.001:2*pi;a=50;b=90;c=120;d=115;xa=0;ya=0;xd=d;yd=0;xb=a.*cos(t);yb=a.*sin(t);m=xd-xb;n=yd-yb;lbd=(m.^2+n.^2).^(1/2);a0=2*b.*(xd-xb);b0=2*b.*(yd-yb);c0=b.^2+lbd.^2-c.^2;dd=2*atan((b0+(a0.^2+b0.^2-c0.^2).^(1/2))./(a0+c0)); plot(dd,'y');d=85,95,105,115 红蓝绿黄。
哈工大机械原理大作业一连杆-15题

Harbin Institute of Technology机械原理大作业一课程名称:机械原理设计题目:连杆传动设计院系:船舶与海洋工程学院班级:1513102班分析者:刘康哲学号:*********指导教师:***设计时间:2017.06.051.题目:如图所示机构,已知机构各构件的的尺寸为l AB=100mm,l BD=400mm,l AE=400mm,l BC=200mm,曲柄AB的角速度为w1=10rad/s,试求构件2上点D的轨迹,构件5的角位移、角速度和角加速度,并对结果进行分析。
(题中构件尺寸满足l BD-l AB<l AE<l BD+l AB)。
2.机构的结构分析及基本杆组划分机构各构件都在同一平面内运动,可拆分成一个原动件(构件1),一个RRPⅡ级基本杆组(构件2、3)和一个RPR Ⅱ级基本杆组(构件4、5)。
活动构件数n=5,P L,P H,因此,机构的自由度为F=3*n-2* P L- P H =3*5-2*7=1原动件:RRPⅡ级基本杆组:RPR Ⅱ级基本杆组:3.坐标系的建立以A点为坐标原点,AE方向为x轴正方向,建立平面直角坐标系。
4.各基本杆组的运动分析数学模型(1)原动件:位置分析: X B = X A + L AB * cos(wt) Y B = Y A + L AB * sin(wt) 其中,X A = Y A = 0 ,w = 10rad/s , L AB = 100 mm 速度分析:将上式分别对时间求导,可得速度方程,V XB = V XA – w*L AB *sin(wt)V YB = V YA + w*L AB *cos(wt)其中,V XA = V YA = 0 加速度分析:将上式分别对时间求导,可得速度方程,a XA = a XA – α* L AB * sin(wt) – w*w*L AB *cos(wt)a YA = a YA + α* L AB * cos(wt) – w*w*L AB *sin(wt)其中,a XA = a YA = 0 α= dw/dt(2)RRP Ⅱ级基本杆组:位置分析: X C = LAB*cos(wt)+wt Lab Lbc 222sin *Y C = 0 可求得,β= arcsin (Y B / L BC ) 速度分析:将上式分别对时间求导,可得速度方程,V XC =-w* LAB*cos(wt) – (wt Lab Lbc 222sin * )-1/2 * sin(wt)*wV YC = 0加速度分析:将上式分别对时间求导,可得加速度方程, (3)RPR Ⅱ级基本杆组位置分析: X D = X C + L CD * cos βY D = - L CD * sin β可求得, θ = arctan (-YD / (LAE -XD) ) 速度分析:将上式分别对时间求导,可得速度方程,5.计算编程在MATLAB 中编写如下程序:function [ XD,YD,Q,Q1,Q2] = calculate(t) %UNTITLED4 Summary of this function goes here % Detailed explanation goes here Lab=100; Lbc=200; Lcd=200; Lae=400; w=10;XB=Lab*cos(w.*t); YB=Lab*sin(w.*t); P=asin(YB./Lbc); XC=sqrt(Lbc.^2-YB.^2)+XB; XD=Lcd.*cos(P)+XC; YD=-Lcd.*sin(P); Q=atan(YD./(Lae-XD)); Q1=diff(Q); Q2=diff(Q1); end5.计算结果t 时间 XD D 点横坐标YD D 点纵坐标Q构件5角度Q1 构件5角速度 Q2构件5角加速度5000.10050.003008D点运动轨迹:构件5的角位移线图:构件5的角速度图:构件5的角加速度线图:6.计算结果分析原动件1可做整周回转,同时构件5也可做整周运动。
哈工大机械原理大作业连杆

哈工大机械原理大作业连杆Modified by JACK on the afternoon of December 26, 2020Harbin Institute of Technology机械原理大作业一课程名称:机械原理设计题目:连杆机构运动分析院系:机电工程学院班级:设计者:学号:指导教师:设计时间:1.运动分析题目(11)在图所示的六杆机构中,已知:AB l =150mm, AC l =550mm, BD l =80mm, DE l =500mm,曲柄以等角速度1w =10rad/s 沿逆时针方向回转,求构件3的角速度、角加速度和构件5的位移、速度、加速度。
2.机构的结构分析建立以点A 为原点的固定平面直角坐标系A-x, y,如下图:机构结构分析该机构由Ⅰ级杆组RR (原动件1)、Ⅱ级杆组RPR (杆2及滑块3)和Ⅱ级杆组RRP (杆4及滑块5)组成。
3.建立组成机构的各基本杆组的运动分析数学模型原动件1(Ⅰ级杆组RR )由图所示,原动件杆1的转角a=0-360°,角速度1w =10rad/s ,角加速度1a =0,运动副A 的位置坐标A x =A y =0,速度(A ,A),加速度(A ,A ),原动件1的长度AB l =150mm 。
求出运动副B 的位置坐标(B x , B y )、速度(B ,B )和加速度(B ,B )。
杆2、滑块3杆组(RPR Ⅱ级杆组)已出运动副B 的位置(B x , B y )、速度(B ,B )和加速度(B ,B ),已知运动副C 的位置坐标C x =0, C y =550mm,速度,加速度,杆长AC l =550mm 。
求出构件2的转角b,角速度2w 和角加速度2a . 构件二上点D 的运动已知运动副B 的位置(B x , B y )、速度(B ,B )、加速度(B ,B ),已经求出构件2的转角b ,角速度2w 和角加速度2a ,杆BD 的长度BD l =80mm 。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
机械原理大作业(一)
作业名称:机械原理
设计题目:连杆机构运动分析
院系:机电工程学院
班级:
设计者:
学号:
指导教师:
设计时间: 2014年6月3日
哈尔滨工业大学机械设计
连杆机构运动分析
(12)题:图1-12所示的六连杆机构中,各构件尺寸分别为:AB l =200mm ,BC l =500mm ,CD l =800mm ,F x =400mm ,D x =350mm ,
D y =350mm ,1 =100rad/s ,求构件5上的F 点的位移、速度和加速度。
1.建立直角坐标系
以F 点为直角坐标系的原点建立直角坐标系X-Y ,如下图所示。
2.机构结构分析
该机构由I级杆组RR(原动件AB)、II级杆组RRR(杆2、3)、II级杆组PRP(杆5、滑块4)组成。
3.各基本杆组运动分析
1.I级杆组RR(原动件AB)
已知原动件AB的转角
ϕ2
π
=
~
原动件AB的角速度
ω
=
10
rad/
s
原动件AB 的角加速度
=α
运动副A 的位置
0,400=-=A A y x
运动副A 的速度
0,0==A A v v
运动副A 的加速度
0,0==A A a a
可得:
)cos(ϕAB A B l x x += )sin(ϕAB A B l y y +=
速度和加速度分析:
)sin(ϕ⨯⨯-=AB xA xB l w v v )
sin(ϕ⨯⨯+=AB yA yB l w v v
)sin()cos(2ϕϕAB AB xA xB el l w a a --= )()s i n (2ϕϕc o a el l w a a AB AB yA yB +-=
2.II 级杆组RRR (杆2、3)
杆2的角位置、角速度、角加速度
BC l =500mm ,CD l =800mm ,D x =350mm ,D y =350mm ,
)
B +A
C B +A B arctan(=ψ20
02
2
02
00-+
)arctan(=ψ3D
C D
C x x y y --
)(L *2A BC 0B D x x -= )
(L *2B BC 0B D y y -=
2B D 2B D 2
BD )y -y ()x -x (L +=
2
CD 2BD 2BC 0L L L C -+=
)
2cos(*L x x BC B C ψ+=
)2sin(*L x y BC B C ψ+=
求导可得C 点的角速度和角加速度。
3.II 级杆组PRP (杆5、滑块4)
)
-arctan(=D C D
C
x x y y --β
=F x
EF
D L x -+)tan(*y =y D F β
求导可得F 点的速度和加速度。
4.计算编程 t=0:pi/50000:pi/50; leng=length(t); a=0;
x=-400+200*cos(a+100*t); y=200*sin(a+100*t); Vx=diff(x); Vy=diff(y); Ax=diff(x,2); Ay=diff(y,2); plot(x,y)
%由B 点和D 点确定杆的角位置、角速度和角加速度: xb=200*cos(a + 100*t) - 400; yb=200*sin(a + 100*t); xd=350; yd=350;
A0=2*500*(xd-xb); B0=2*500*(yd-yb);
C0=500*500+(xd-xb).^2+(yd-yb).^2-800*800;
Q1=2.*atan(B0+(A0.*A0+B0.*B0-C0.*C0).^0.5/(A0+C0)); xc=xb+500*cos(Q1);
yc=yb+500*sin(Q1);
Q2=atan((yc-yd)./(xc-xd));
w1=diff(Q1);
w2=diff(Q2);
Vcx=diff(xc);
Vcy=diff(yc);
O1=diff(Q1,2);
O2=diff(Q2,2);
Acx=diff(xc,2);
Acy=diff(yc,2);
%已知杆CD的运动和移动副导路参数推出F点的运动:EF=250;
xf=0;
yf=350+350*tan(pi-Q2)-EF;
Vfx=0;
Vfy=diff(yf);
Afx=0;
Afy=diff(Vfy);
for(m=1:(leng-1))
t1(m)=t(m);
end
for(m=1:(leng-2)) t2(m)=t(m);
end
5、计算结果图令a=0,EF=250;F点位移:
F点y方向位移图
F点速度图:
F点加速度图:。