哈工大机械原理大作业——连杆——24号

合集下载

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

哈尔滨工业大学机械原理大作业_连杆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字形”。

哈工大机械原理大作业-连杆

哈工大机械原理大作业-连杆

哈工大机械原理大作业-连杆本页仅作为文档页封面,使用时可以删除This document is for reference only-rar21year.MarchHarbin Institute of Technology机械原理大作业一课程名称:机械原理设计题目:连杆运动分析院系:机电工程学院班级:设计者:学号:指导教师:陈明设计时间:2013年6月25日1、运动分析题目在图1-10中所示的干草压缩机中,已知LAB=150mm,LBC=600mm,LCE=120mm,LCD=500mm,LEF=600mm,XD=400mm,YD=500mm,YF=600mm,曲柄1作等速转动,其转速n1=50r/min。

求在一个运动循环中活塞5的位移、速度和加速度的变化曲线。

图1-102、机构的结构分析(1)基本杆组的划分①AB即杆件1为原动件②DECB即杆件2、3为RRR型II级杆组,其中CE为同一构件上点。

③ EF 和滑块即4、5为RRP 型II 级杆组(2)、建立以点A 为原点的固定平面直角系3、确定已知参数和求解流程(1)原动件1(I 级杆组RR )如图所示,已知原动件1的转角 πϕ2~01=原动件杆1的角速度 s rad /236.51=ω原动件1的角加速度01=α运动副A 的位置坐标0,0==A A y x运动副A 的速度0,0==yA xA v v运动副A 的加速度0,0==yA xA a a原动件杆I 的长度mm l 1501=可求出B 的位置B 的速度B 的加速度(2)构件2、3(II 级杆组RRR )D 的位置500400==D D y xD 的速度 0==yD xD v vD 的加速度00==yD xD a a杆长mm l l CD j 500==,mm l l BC i 600== 由关系j j D i i B C l x l x x ϕϕcos cos +=+= j j D i i B C l y l y y ϕϕsin sin +=+= 其中)(20sin cos B D i i i x x l A C B A -==-+ ϕϕ)(2B D i y y l B -=222j BD i l l l C -+=222)()(B D B D BD y y X x l -+-= 可解得C A C B A B i +-++=222arctan 2ϕDC DC j x x y y --=arctan ϕ由上面两个式子可以得到两杆的角速度1)()(G y y s x x c w B D j B D j i i -+-==•ϕ132G s G c Gij i i +==ϕα其中i j j i s c s c G -=1,i i i l c ϕcos =,i i i l s ϕsin =,j j j l c ϕcos =,j j j l s ϕsin =可得E 的位置iB E iB E y y x x ϕϕsin 480cos 480+=+=E 的速度 i i yB E yE i i xB E xE v yv v xv ωϕωϕcos 480sin 480+==-==E 的加速度i i i i xB E xE a xa αϕωϕsin 480cos 4802--== i i i i yB E yE a ya ϕαϕωcos 480sin 4802+-== (3)、构件4、5杆组(II 级杆组RRP )在建立的坐标系中取一参考点K600,0==K K y x 则速度0,0==yK xK v v加速度0,0==yK xK a a杆长mm l EF 600=,设F 位移为s 由s x l x x K i EF E F +=+='cos ϕK i EF E F y l y y =+='sin ϕ由上面两个式子可以得到600600arcsin arcsin'E EF E k i y l y y -=-=ϕ 所以:F 点位移 )600600arcsin(cos 600E E F y x x s -+== 速度F F xv = 加速度F F xa =四、编程计算并输出结果(VB编程)主程序:Private Sub Command1_Click()Dim s5(3600) As DoubleDim v5(3600) As DoubleDim a5(3600) As Double Dim pi As DoubleDim pa As Doublepi = 3.1415926pa = pi / 180Dim i As LongDim f1(3600) As Double Dim RR1 As RRDim RR2 As RRDim RRR1 As RRRDim RRP1 As RRPSet RR1 = New RRSet RR2 = New RRSet RRR1 = New RRRSet RRP1 = New RRPFor i = 0 To 3600 Step 1 f1(i) = i * pa / 10RR1.delt = 0RR1.f = f1(i)RR1.w = 5.24RR1.e = 0RR1.L = 150RR1.xA = 0RR1.yA = 0RR1.vxA = 0RR1.vyA = 0RR1.axA = 0RR1.ayA = 0RR1.calRRR1.Li = 600RRR1.Lj = 500RRR1.xB = RR1.xBRRR1.yB = RR1.yBRRR1.vxB = RR1.vxBRRR1.vyB = RR1.vyBRRR1.axB = RR1.axBRRR1.ayB = RR1.ayBRRR1.xD = 400RRR1.yD = 500RRR1.vxD = 0RRR1.vyD = 0RRR1.axD = 0RRR1.ayD = 0RRR1.M = 1RRR1.calRRRRR2.delt = 0RR2.f = RRR1.fi RR2.w = RRR1.wi RR2.e = RRR1.ei RR2.L = 480RR2.xA = RR1.xB RR2.yA = RR1.yB RR2.vxA = RR1.vxB RR2.vyA = RR1.vyB RR2.axA = RR1.axB RR2.ayA = RR1.ayB RR2.calRRP1.Li = 600RRP1.Lj = 0RRP1.fj = piRRP1.wj = 0RRP1.ej = 0RRP1.xB = RR2.xB RRP1.yB = RR2.yB RRP1.vxB = RR2.vxB RRP1.vyB = RR2.vyB RRP1.axB = RR2.axB RRP1.ayB = RR2.ayB RRP1.xK = 0RRP1.yK = 600RRP1.vxK = 0RRP1.vyK = 0RRP1.axK = 0RRP1.ayK = 0RRP1.M = 1RRP1.cals5(i) = RRP1.ssv5(i) = RRP1.vssa5(i) = RRP1.assNext iPicture1.Scale (-30, 700)-(360, 580)Picture1.Line (0, 0)-(360, 0) 'XPicture1.Line (0, 580)-(0, 700) 'YFor i = 0 To 360 Step 10 'X轴坐标Picture1.DrawStyle = 2Picture1.Line (i, 700)-(i, 580)Picture1.CurrentX = i - 10: Picture1.CurrentY = 0 Picture1.Print iNext iFor i = 580 To 700 Step 10 'Y轴坐标Picture1.DrawStyle = 2Picture1.Line (0, i)-(360, i)Picture1.CurrentX = -10: Picture1.CurrentY = iPicture1.Print iNext iFor i = 0 To 3600 Step 1Picture1.PSet (i / 10, s5(i))Next iEnd SubRR:Public L As DoublePublic f As DoublePublic delt As DoublePublic w As DoublePublic e As DoublePublic xA As DoublePublic yA As DoublePublic vxA As DoublePublic vyA As DoublePublic axA As DoublePublic vxB As DoublePublic vyB As DoublePublic axB As DoublePublic ayB As DoublePublic Sub cal()xB = xA + L * Cos(f + delt)yB = yA + L * Sin(f + delt)vxB = vxA - w * L * Sin(f + delt)vyB = vyA + w * L * Cos(f + delt)axB = axA - w ^ 2 * L * Cos(f + delt) - e * L * Sin(f + delt) ayB = ayA - w ^ 2 * L * Sin(f + delt) + e * L * Cos(f + delt) End SubRRR:Public Li As DoublePublic Lj As DoublePublic fi As DoublePublic fj As DoublePublic wi As DoublePublic wj As DoublePublic ei As DoublePublic ej As DoublePublic xB As DoublePublic yB As DoublePublic vxB As DoublePublic vyB As DoublePublic axB As DoublePublic ayB As DoublePublic xC As DoublePublic yC As DoublePublic vxC As DoublePublic vyC As DoublePublic axC As DoublePublic ayC As DoublePublic xD As DoublePublic yD As DoublePublic vxD As DoublePublic vyD As DoublePublic M As DoublePublic Sub calRRR()Dim fDB As DoubleDim Ci As DoubleDim Cj As DoubleDim Si As DoubleDim Sj As DoubleDim G1 As DoubleDim G2 As DoubleDim G3 As DoubleDim LBD As DoubleDim JCBD As DoubleDim val As Doublepi = 3.1415926LBD = Sqr((xB - xD) ^ 2 + (yD - yB) ^ 2)If LBD < Li + Lj And LBD > Abs(Li - Lj) Thenval = (Li ^ 2 + LBD ^ 2 - Lj ^ 2) / (2 * Li * LBD) JCBD = Atn(-val / Sqr(-val * val + 1)) + 2 * Atn(1) End IfRRP:Public Li As DoublePublic Lj As DoublePublic fi As DoublePublic fj As DoublePublic wi As DoublePublic wj As DoublePublic ei As DoublePublic ej As DoublePublic xB As DoublePublic yB As DoublePublic vxB As DoublePublic vyB As DoublePublic axB As DoublePublic ayB As DoublePublic xK As DoublePublic vyK As DoublePublic axK As DoublePublic ayK As DoublePublic xC As DoublePublic yC As DoublePublic vxC As DoublePublic vyC As DoublePublic axC As DoublePublic ayC As DoublePublic xD As DoublePublic yD As DoublePublic vxD As DoublePublic vyD As DoublePublic axD As DoublePublic ayD As DoublePublic M As SinglePublic ss As DoublePublic vss As DoublePublic ass As DoublePublic Sub cal()Dim A0 As DoubleDim Q1 As DoubleDim Q2 As DoubleDim Q3 As DoubleDim Q4 As DoubleDim Q5 As DoubleDim val As DoubleDim pi As Doublepi = 3.14159216A0 = Lj + ((yK - yB) * Cos(fj) - (xK - xB) * Sin(fj)) val = A0 / Lifi = M * Atn(val / Sqr(-val * val + 1)) + fjxC = xB + Li * Cos(fi)yC = yB + Li * Sin(fi)ss = (xC - xK) * Cos(fj) + (yC - yK) * Sin(fj)xD = xK + ss * Cos(fj)yD = yK + ss * Sin(fj)Q1 = vxK - vxB - wj * (ss * Sin(fj) + Lj * Cos(fj))Q2 = vyK - vyB + wj * (ss * Cos(fj) - Lj * Sin(fj))Q3 = Li * Sin(fi) * Sin(fj) + Li * Cos(fi) * Cos(fj)wi = (-Q1 * Sin(fj) + Q2 * Cos(fj)) / Q3vss = -(Q1 * Li * Cos(fi) + Q2 * Li * Sin(fi)) / Q3vxC = vxB - wi * Li * Sin(fi)vyC = vyB + wi * Li * Cos(fi)vxD = vxK + vss * Cos(fj) - ss * wj * Sin(fj)vyD = vyK + vss * Sin(fj) + ss * wj * Cos(fj)Q4 = axK - axB + wi ^ 2 * Li * Cos(fi) - ej * (ss * Sin(fj) + Lj * Cos(fj)) - wj ^ 2 * (ss * Cos(fj) - Lj * Sin(fj)) - 2 * vss * wj * Sin(fj)Q5 = ayK - ayB + wi ^ 2 * Li * Sin(fi) + ej * (ss * Cos(fj) - Lj * Sin(fj)) - wj ^ 2 * (ss * Sin(fj) + Lj * Cos(fj)) + 2 * vss * wj * Cos(fj)ei = (-Q4 * Sin(fj) + Q5 * Cos(fj)) / Q3ass = (-Q4 * Li * Cos(fi) - Q5 * Li * Sin(fi)) / Q3axC = axB - ei * Li * Sin(fi) - wi ^ 2 * Li * Cos(fi)ayC = ayB + ei * Li * Cos(fi) - wi ^ 2 * Li * Sin(fi)axD = axK + ass * Cos(fj) - ss * ej * Sin(fj) - ss * wj ^ 2 * Cos(fj) - 2 * vss * wj * Sin(fj)ayD = ayK + ass * Sin(fj) + ss * ej * Cos(fj) - ss * wj ^ 2 * Sin(fj) + 2 * vss * wj * Cos(fj)End Sub五、计算结果数据如图:位移曲线:速度曲线:加速度曲线:六、计算结果分析主动件转角为0时,滑块的位移为628mm,随着转角的匀速增加,滑块位移先上升,速度为负快速下降,加速度为负且开始值较小并逐渐下降,到达最低点-2200左右,此时速度为0, 然后位移开始继续下降,速度继续下降,然后到达最小值-100左右,此时加速度为0,又开始上升(向左运动),速度正向增大,在后面一段时期速度继续增大,加速度也正向增大,然后速度下降,加速度下降.滑块分别在4.2s,11.4s速度达到正向最大,1.7s,7.7s速度达到负向最大,0.8s,2.7s,5.8s,9.6s加速度达到极值,滑块就是这样周期性的左右运动。

哈尔滨工业大学机械原理大作业-连杆机构运动分析

哈尔滨工业大学机械原理大作业-连杆机构运动分析
end;
%打印图像
figure(1);
plot(fi1,sF);
title('位移变化曲线');
figure(2);
plot(fi1,vF);
title('速度变化曲线');
figure(3);
plot(fi1,aF);
title('加速度变化曲线');
六、计算结果
图8:推杆位移变化曲线
图9:推杆速度变化曲线
xD=400;yD=500;vDx=0;vDy=0;aDx=0;aDy=0;
xK=0;yK=600;vKx=0;vKy=0;aKx=0;aKy=0;
l1=150;l2=600;l3=500;lBE=480;l4=600;
fi5=pi;
n1=50;
w1=2*pi*n1/60;
fi1=linspace(0,2*pi,1000);
aEy(i)=aBy(i)-w2(i)^2*lBE*sin(fi2(i))+a2(i)*lBE*cos(fi2(i));
%求F点的运动参数
A1(i)=(yE(i)-yK)*cos(fi5)-(xE(i)-xK)*sin(fi5);
fi4(i)=fi5-asin(A1(i)/l4);
xF(i)=xE(i)+l4*cos(fi4(i));
(2)速度和加速度分析
B点的速度
B点的加速度
2、由“RRRⅡ级杆组”,已知B点和D点的运动参数,可求得构件2、构件3的运动参数
图5
D点的坐标方程
D点的速度
D点的加速度
(1)位置方程
构件2、构件3的长度
先求出 和
其中
则可求得

哈工大机械设计大作业平面连杆机构设计说明书

哈工大机械设计大作业平面连杆机构设计说明书

机械设计基础大作业计算说明书题目:平面连杆机构设计学院:材料学院班号:学号:姓名:日期: 2014年9月30日哈尔滨工业大学机械设计基础大作业任务书题目:平面连杆机构设计设计原始数据及要求:()目录1设计题目2设计原始数据3设计计算说明书3.1计算极位夹角3.2设计制图3.3验算最小传动角4参考文献1 设计题目平面连杆机构的图解法设计2 设计原始数据设计一曲柄摇杆机构。

已知摇杆长度,摆角,摇杆的行程速比系数,要求摇杆靠近曲柄回转中心一侧的极限位置与机架间的夹角为,试用图解法设其余三杆的长度,并检验(测量或计算)机构的最小传动角。

()3 设计计算说明书3.1 计算极位夹角极位夹角代入数值3.2 设计制图3.2.1 在图纸上取一点作为点,从点垂直向上引出一条长为的线段,终点为;3.2.2 从点在左侧引出一条与夹角为的射线;3.2.3 以点为圆心,以为半径画圆,与射线交于点;3.2.4 分别从、两点向下引两条射线,射线与夹角为,两射线交于点,点即为曲柄的回转中心;3.2.5 以点为圆心以为半径画圆;3.2.6 过点向左侧引出一条射线,射线与夹角,与圆交于点;3.2.7 连接,并量取其长度,以为半径画圆,直线,与圆的交点分别为,;3.2.8 在图中量取,,3.3 验算最小传动角3.3.1 在处根据余弦定理3.3.2 在处根据余弦定理所以最小传动角4 参考文献[1]宋宝玉,王瑜,张锋主编.机械设计基础.哈尔滨:哈尔滨工业大学出版,2010.[2]王瑜主编.机械设计基础大作业指导书.哈尔滨:哈尔滨工业大学出版社,2014.。

哈工大机械原理大作业连杆

哈工大机械原理大作业连杆

哈工大机械原理大作业-连杆连杆是机械原理中常见的机构之一,也是机械工程中非常重要的部件。

它由两个旋转接头和一个连接两个旋转接头的杆件组成。

连杆广泛应用于各种机械设备中,如汽车发动机、泵、机床等。

本文将介绍连杆的工作原理、应用以及设计要点。

连杆的工作原理是将旋转运动转化为直线运动或将直线运动转化为旋转运动。

它通过两个旋转接头的运动将杆件上的一个点的运动转化为另一个点的运动。

连杆的运动有两种基本形式:一是曲柄连杆机构,二是摇杆连杆机构。

曲柄连杆机构中,一个旋转接头为曲柄,另一个旋转接头为连杆;摇杆连杆机构中,一个旋转接头为摇杆,另一个旋转接头为连杆。

连杆广泛应用于各种机械设备中。

在汽车发动机中,连杆将曲轴的旋转运动转化为活塞的直线运动,从而驱动汽缸的工作;在泵中,连杆将电机的旋转运动转化为柱塞的直线运动,从而产生压力;在机床中,连杆将电机的旋转运动转化为工作台的直线运动,从而实现加工。

设计连杆时需要考虑一些要点。

首先是连杆的材料选择和尺寸设计。

连杆需要承受较大的力和扭矩,因此需要选择具有较高强度和刚度的材料。

同时,根据应用需求和力学原理,设计连杆的尺寸,以确保其能够承受正常工作条件下的负荷。

其次是连杆的润滑和密封。

连杆在工作过程中需要润滑剂来减少摩擦和磨损,同时需要密封装置来防止润滑剂泄漏。

因此,设计连杆时需要考虑润滑剂的供给和密封装置的设计。

最后是连杆的制造和装配。

连杆的制造需要保证其精度和质量,以确保其运转平稳和可靠。

在装配过程中,需要按照设计要求进行装配,同时进行必要的调试和检测,以确保连杆的工作性能符合要求。

总之,连杆是机械工程中非常重要的部件,广泛应用于各种机械设备中。

设计和制造连杆需要考虑材料选择、尺寸设计、润滑和密封以及制造和装配等方面的要点。

通过合理的设计和制造,可以确保连杆的工作性能和可靠性,从而提高机械设备的工作效率和寿命。

哈工大机械原理连杆和凸轮大作业24题

哈工大机械原理连杆和凸轮大作业24题

班级学号机械原理大作业说明书题目 1、连杆机构运动分析2、凸轮机构设计学生姓名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 =-BxB 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 dt dx C B -= 1s i n 131c o s 13.....2ϕϕϕϕl l x dt x d C B --= 1cos 31..ϕϕl y dt dy C B += 1c o s 131s i n 13.....2ϕϕϕϕl l y dty d C B +-= 根据关系:1111d 122..11.αϕϕωϕϕ====dtd dt , 故可得出: D x =)1cos(4βϕ++l x CD y =)1sin(4βϕ++l y Cb. 如右图所示是由2个回转副和1个移动副组成的II级组。

哈工大-机械原理大作业-连杆机构运动分析

哈工大-机械原理大作业-连杆机构运动分析

机械原理大作业(一)作业名称:连杆机构运动分析设计题目: 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α。

哈工大机械原理大作业一连杆运动分析(02)

哈工大机械原理大作业一连杆运动分析(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 红蓝绿黄。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

Harbin Institute of Technology
机械原理设计说明书
课程名称:机械原理
设计题目:连杆机构运动分析
指导老师:陈明丁刚
班级:
设计者:
学号:
指导教师:
设计时间:
一、运动分析题目
如图1-6所示是曲柄滑块机构,各构件长度分别为a、b,偏距为e,连杆BC上
一点到铰链B的距离为l
m ,试研究各构件长度a、b、l
m
及偏距e的变化对点m
的轨迹的影响规律。

二、机构结构分析及基本杆组划分
1.除去虚约束力和局部自由度,计算机构的自由度并确定原动件
本机构中无虚约束或局部自由度。

机构各杆件都在同一平面运动,活动构件数
n=3,P
L =4,P
H
=0,则机构的自由度为:F=3n-2P
L
-P
H
=33-24-10=1。

原动
件为曲柄AB。

2.拆分杆组
从远离原动件(即杆AB,如图1)进行拆分,可得到由杆BC和C点处滑块组成的RRPⅡ级杆组(如图2),剩下的就是Ⅰ级机构杆AB。

3.确定机构的级别
由上可知,机构为Ⅱ级机构
⨯⨯⨯⨯⨯
三、各基本杆组的运动分析数学模型
(1)原动件AB (Ⅰ级杆组)
原动件AB 的转角为:i ϕ=0~2π;角速度为:s rad /101=ω 角加速度为:01=ε
假定运动副A 的位置坐标为:x A =0,y A =0 A 点与机架相连,
即该点的速度和加速度都为0。

原动件AB 长度为l i ,从而可求得运动副B 点的位置坐标:i
ϕcos l x x i A B +=,
i i A B l y y ϕsin +=
(2)杆BC 和C 点的滑块(RRP Ⅱ级杆组)
RRP Ⅱ级杆组是由两个构件两个转动副及一个外移动副组成的。

已知两杆长为l i 和l j (l j 杆垂直于滑块导路),外回转副B 的参数,滑块导路方向角和计算位移时的参考点K 的位置和导路的运动参数,求内运动副C 的运动参数。

位置方程:内回转副C 的位置方程为:
j j j K i i B C j j j K i i B C l s y l y y l s x l x x ϕϕϕϕϕϕcos sin sin sin cos cos ++=+=-+=+=
消去s
可得:j i
j
i l l A ϕϕ++=0arcsin
其中
j K B j K B y y x x A ϕϕcos )(sin )(0---=
为保证机构的存在,应满足装配条件i j l l A ≤+0,求得
i ϕ后,可求得x C 和y C ,而后求得滑块的位移s :
j j j K C j j j K C l y y l x x s ϕϕϕϕsin /)cos (cos /)sin (--=+-=
滑块D 点的位移方程为:
j
K D j K D s y y s x x ϕϕsin cos +=+=
(3)求M 点位置坐标
要画出点m 的轨迹图,需要求出点m 的位置坐标。

假定以A 点为原点,则点B 的位置坐标为:
1
1sin cos ααa y a x B B == 点C
的位置坐标为:
e
y b a x C C =+=21cos cos αα 而点m 是杆BC 上的一点,且到B 点的距
离为l m ,则B,C,m 三点满足关系:C
B m
B B
C B m m y y y y x x x x b l --=--= 从而求得m 点的位置坐
标为:m
m m m m l b
e a b l y l a x +-=+=12
1sin )1(cos cos ααα
四、建立坐标系
五、计算编程
利用matlab 编程,代码如下
f1=[0:pi/180:2*pi];%原动件AB 杆转角范围0~2*pi w1=1;%AB 杆的角速度 e1=0;%AB 杆的角加速度
xA=0;yA=0;%运动副A 的位置
e=200;%偏距e ,改变e 的大小可以观察偏距e 对m 点轨迹的影响 l1=50;%杆AB 的长,可通过修改l1来观察长度a 对m 点轨迹的影响 l2=100;%杆BC 的长,可通过修改l2来观察长度b 对m 点轨迹的影响
lm=30;%BC 杆上的m 点到B 点的距离,可通过修改lm 来观察长度l m 对m 点轨迹的影响
r1=w1*f1;%AB 杆与x 轴正方向的夹角
r2=asin((e -l1*sin(r1))/l2);%360。

-BC 杆与x 轴正方向所夹锐角 xm=lm*cos(r2)+l1*cos(r1);%m 点的横坐标
ym=(1-lm/l2)*l1*sin(r1)+e*lm/l2;%m 点的纵坐标 plot(xm,ym); hold off grid on
title('点m 的轨迹图','FontSize',20);
六、计算结果
(1)构件a长度对m运动轨迹的影响(b=100,e=20,l
m
=30)
(2)构件B长度对m运动轨迹的影响(a=20,e=20,l
m
=100)
(3)l
m
长度对m点运动轨迹的影响(a=50,b=100,e=20)
(3)e的大小对m点运动轨迹的影响(a=50,b=100,l
m =30
a=20 a=50 a=80 a=100 b=20 b=50 b=80 b=100 l=20 l=50 l=80 l=90
e=20 e=50 e=40 e=50
七、计算结果分析
a的长度对m运动轨迹的影响
当构件长度a增大时,m点轨迹图的横向的左右两端会变长,纵向的上下两端也会变宽,并且左端的凸起会略微向下移动,右端凸起会略微向右移,而上端和下端凸起都会向左移。

图形整体上下左右都变大。

b的长度对m运动轨迹的影响
当构件长度b增大时,m点轨迹图的横向的左右两端长度基本保持不变,纵向的上下两端会变宽,并且左端的凸起会略微向右和向下移动,右端凸起会略微向右和向下移,而上端和下端凸起都向右移。

图形整体上下变宽,向右移。

L m的长度对m运动轨迹的影响
长度l
增大时,m点轨迹图的横向的左右两端会略微增长,纵向的上下两端会变m
窄,并且左端和右端的凸起都会向上和向右移动,而上端的凸起会向右移和向下移,下端的凸起会向右和向上移。

图形整体上下变窄,向右移动。

e的长度对m运动轨迹的影响
e增大时,m点轨迹图的横向的左右两端会略微变长,纵向的上下两端不变,并且左端和右端的凸起都会向左和向上移动,而上端的凸起会向右移和向上移,下端的凸起会向左和向上移。

图形整体向上移动。

相关文档
最新文档