哈工大机械原理大作业2凸轮机构27题
机械原理大作业2凸轮(18题)
1、运动分析题目如图1所示直动从动件盘形凸轮机构，其原始参数见表1。
图 1表 1序号 升程 （mm ） 升程运动角（） 升程运动规律 升程许用压力角（） 回程运动角（）回程运动规律 回程许用压力角（） 远休止角（） 近休止角（）18 100150正弦加速度 30 100等减等加速 6040702、凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图：2.1从动件运动方程：（1）从动件升程运动方程升程段采用正弦加速度运动规律，运动方程为：()1212112100sin 5/6251001251cos 05/656210012sin 55/6s v a ϕϕππωπϕϕππωϕπ⎫⎡⎤⎛⎫=-⎪⎪⎢⎥⎝⎭⎣⎦⎪⎪⎡⎤⎪⎛⎫⎛⎫=-≤≤⎬ ⎪⎪⎢⎥⎝⎭⎝⎭⎣⎦⎪⎪⨯⎛⎫⎪= ⎪⎝⎭⎪⎭（2）从动件远休止运动方程在远休止s Φ段，即5/619/18πϕπ≤≤时，100s h mm ==,0v =,0a =。
（3）从动件回程运动方程升程段采用等减等加运动规律，运动方程为：()221221220019100518()94001919/184/3518()94005()9s v a πϕπωπϕπϕππωπ⎫⎪⎡⎤=--⎪⎢⎥⎣⎦⎪⎪⎪⎪⎡⎤=--≤≤⎬⎢⎥⎣⎦⎪⎪⎪⎪=-⎪⎪⎭()221221220029()518()94002()4/329/1853()94005()9s v a πϕπωπϕπϕππωπ⎫⎪=-⎪⎪⎪⎪⎪=--≤≤⎬⎪⎪⎪⎪=⎪⎪⎭（4）从动件近休止运动方程在近休止s 'Φ段，即29/182πϕπ≤≤时，0s =,0v =,0a =。
2.2推杆位移、速度、加速度线图：（1）推杆位移线图图 2 推杆位移线图（2）推杆速度线图图 3 推杆速度线图（3）推杆加速度线图图 4 推杆加速度线图3、凸轮机构的ds s d ϕ-线图，并由此确定凸轮的基圆半径和偏距：图 5 凸轮机构的dss d ϕ-线图 4滚子半径的确定及凸轮理论廓线和实际廓线的绘制 4.1凸轮的理论轮廓方程为：00()cos sin (02)()sin cos x s s e y s s e ϕϕϕπϕϕ=+-⎫≤≤⎬=++⎭式中，220031.45s r e mm =-=（1）推程凸轮轮廓方程：11231100sin cos 18sin 5/625(05/6)112(31100sin )sin 18cos 5/625x y ϕϕϕϕππϕπϕϕϕϕππ⎫⎧⎫⎡⎤⎛⎫=+--⎨⎬⎪⎪⎢⎥⎝⎭⎣⎦⎪⎩⎭≤≤⎬⎡⎤⎛⎫⎪=+-+ ⎪⎢⎥⎪⎝⎭⎣⎦⎭（2）远休止凸轮轮廓方程：131cos 18sin (5/619/18)131sin 18cos x y ϕϕπϕπϕϕ=-⎫≤≤⎬=+⎭（3）回程凸轮轮廓方程：222220019131cos 18sin 518()9(19/184/3)20019131sin 18cos 518()9x y πϕϕϕππϕππϕϕϕπ⎫⎧⎫⎪⎪⎪⎡⎤=---⎪⎨⎬⎢⎥⎣⎦⎪⎪⎪⎩⎭⎪≤≤⎬⎧⎫⎪⎪⎪⎡⎤⎪=--+⎨⎬⎢⎥⎪⎣⎦⎪⎪⎪⎩⎭⎭22222002931()cos 18sin 518()9(19/184/3)2002931()sin 18cos 518()9x y πϕϕϕππϕππϕϕϕπ⎫⎧⎫⎪⎪⎪=+--⎪⎨⎬⎪⎪⎪⎩⎭⎪≤≤⎬⎧⎫⎪⎪⎪⎪=+-+⎨⎬⎪⎪⎪⎪⎩⎭⎭（4）近休止凸轮轮廓方程：31cos 18sin (29/182)31sin 18cos x y ϕϕπϕπϕϕ=-⎫≤≤⎬=+⎭4.2凸轮理论轮廓曲线为：图 6 凸轮理论轮廓由上图可编程可求其最小曲率半径为min 10.309110mm ρ=≈，所以滚子半径min 1037r r mm ρ=-∆=-=。
哈尔滨工业大学机械原理大作业_连杆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字形”。
机械原理大作业2-1120810417-凸轮
机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电工程学院班级：1208104完成者：学号：1120810417指导教师：林琳刘福利设计时间：2014年6月2日哈尔滨工业大学一、设计题目如下图所示为直动从动件盘形凸轮机构，据此设计该凸轮机构：二、原始参数 序号升程升程运动角 升程运动规律 升程许用压力角 回程运动角 回程运动规律 回程许用压力角 远休止角 近休止角 15 90mm150°正弦加速度30°100°余弦加速度60°55°55°三、推杆升程方程和推杆回程方程： 在这里取ω=1rad/s. （1）推杆升程方程：650,)512sin(215690)(πφφππφφ≤≤⎥⎦⎤⎢⎣⎡-=s 650),512cos(108)(πφφφπφν≤≤-=650,512sin 2.259)(πφφπφ≤≤=a（2）推杆回程方程：36613641,)05.059cos(145)(πφππφφ≤≤⎥⎦⎤⎢⎣⎡-+=s ω36613641,)05.059sin(181)(πφππφφν≤≤⎥⎦⎤⎢⎣⎡---= 36613641),05.059cos(8.145)(≤≤--=φππφφa四、matlab 程序及曲线图像注：每一段都为完整程序，可直接运行。
1.推杆位移曲线clear allp1=0:pi/360:(5*pi/6-pi/360); w=1;s1=90*(6*p1/(5*pi)-1/(2*pi)*sin(12*p1/5)); p2=5*pi/6:pi/360:(41*pi/36-pi/360); s2=90*ones(1,length(p2));p3=41*pi/36:pi/360:(61*pi/36-pi/360); s3=45*(1+cos(9*p3/5-1*pi/20)); p4=61*pi/36:pi/360:2*pi; s4=0*p4;p=[p1,p2,p3,p4]; s=[s1,s2,s3,s4];plot(p,s)xlabel('Φ(角度)');ylabel('S(位移)'); title('推杆位移曲线');2.推杆速度曲线clear allp1=0:pi/360:(5*pi/6-pi/360);w=1;v1=108*w/pi*(1-cos(12*p1/5));p2=5*pi/6:pi/360:(41*pi/36-pi/360);v2=0*p2;p3=41*pi/36:pi/360:(61*pi/36-pi/360);v3=-81*w*sin(9*p3/5-1*pi/20);p4=61*pi/36:pi/360:2*pi;v4=0*p4;p=[p1,p2,p3,p4];v=[v1,v2,v3,v4];plot(p,v)xlabel('Φ(角度)');ylabel('V(速度)'); title('推杆速度曲线');3.推杆加速度曲线clear allp1=0:pi/360:(5*pi/6-pi/360);w=1;a1=36*36*w^2/5/pi*sin(12*p1/5);p2=5*pi/6:pi/360:(41*pi/36-pi/360);a2=0*p2p3=41*pi/36:pi/360:(61*pi/36-pi/360);a3=-18*81*w^2/10*cos(9*p3/5-1*pi/20);p4=61*pi/36:pi/360:2*pi;a4=0*p4;p=[p1,p2,p3,p4];a=[a1,a2,a3,a4];plot(p,a)xlabel('Φ(角度)');ylabel('a(加速度)'); title('推杆加速度曲线');4.凸轮机构的ds/dφ-s线图clear allp1=0:pi/360:(5*pi/6-pi/360);w=1;s1=90*(6*p1/(5*pi)-1/(2*pi)*sin(12*p1/5)); p2=5*pi/6:pi/360:(41*pi/36-pi/360);s2=90*ones(1,length(p2));p3=41*pi/36:pi/360:(61*pi/36-pi/360);s3=45*(1+cos(9*p3/5-1*pi/20));p4=61*pi/36:pi/360:2*pi;s4=0*p4;p=[p1,p2,p3,p4];s=[s1,s2,s3,s4];p1=0:pi/360:(5*pi/6-pi/360);w=1;v1=108*w/pi*(1-cos(12*p1/5));p2=5*pi/6:pi/360:(41*pi/36-pi/360);v2=0*p2;p3=41*pi/36:pi/360:(61*pi/36-pi/360);v3=-81*w*sin(9*p3/5-1*pi/20);p4=61*pi/36:pi/360:2*pi;v4=0*p4;p=[p1,p2,p3,p4]; v=[v1,v2,v3,v4]; vx=-v; hold on plot(vx,s)%直线Dtdty=-100:0.01:100; x=-69; hold onplot(x,y,'-r'); % 直线Dt’dt’ x=-100:0.01:100; y=-0; hold onplot(x,y,'-r'); grid on hold offtitle('ds/d φ-s 曲线');曲线为升程阶段的类速度-位移图，根据升程压力角与回城压力角做直线与其相切，, 其直线斜率分别为：K 1=)30150tan(+=0 K 2=)60150tan(-为∞；两直线方程为： }{0,69=-=y x进而确定凸轮偏距和基圆半径：在轴心公共许用区内取轴心位置，能够满足压力角要求，由图可得:取s0=200mm ，e=30；r0=（2002 ＋502）1/2=206.2mmclear allp1=0:pi/360:(5*pi/6-pi/360);w=1;s1=90*(6*p1/(5*pi)-1/(2*pi)*sin(12*p1/5));p2=5*pi/6:pi/360:(41*pi/36-pi/360);s2=90*ones(1,length(p2));p3=41*pi/36:pi/360:(61*pi/36-pi/360);s3=45*(1+cos(9*p3/5-1*pi/20));p4=61*pi/36:pi/360:2*pi;s4=0*p4;p=[p1,p2,p3,p4];s=[s1,s2,s3,s4];s0=200;e=30;x=(s0+s).*cos(p)-e*sin(p);y=(s0+s).*sin(p)+e*cos(p);plot(x,y)title('凸轮理论轮廓');6.凸轮实际轮廓工作轮廓曲率半径ρ、理论轮廓曲率半径ρ与滚子半径r三者存在如下关系aρa=ρ+r，不妨最终设定滚子半径为30mm，这时滚子与凸轮间接触应力最小，可提高凸轮寿命。
机械原理大作业二-凸轮机构设计..
机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机械设计制造及其自动化班级：1208104完成者：郑鹏伟学号：**********指导教师：林琳刘福利设计时间：2014年6月4日哈尔滨工业大学一、 设计题目：凸轮的机构运动简图如下图所示:序 号 升程 （mm ） 升程运动角（°） 升程运 动规律 升程 许用压力角（°）回程运动角 （°） 回程运动规律回程许用压力角（°） 远休止角（°） 近休止角（°） 14 90120余弦 加速 度3590等减等加速657575二、 凸轮推杆升程、回程运动方程及推杆位移，速度加速度线图：（1）凸轮推杆升程运动方程： 根据题意知：00120759075s s Φ=Φ='Φ='Φ=（1）从动件升程运动方程（设为1rad sω=）122212s [1cos()]2sin()2cos()2h h h a πϕπωπυϕπωπϕ=-Φ=ΦΦ=ΦΦ（2）从动件远休止运动方程在远休止s Φ段，即213312πϕπ≤≤时，90,0,0s h mm a υ====。
（3）从动件回程运动方程升程段采用等减等加运动规律，运动方程为：①当回程0002s s ϕ'ΦΦ+Φ≤≤Φ+Φ+134()123πϕπ≤≤时： 20s 201022122[-+]4[()]4s hs h h h a ϕωυϕω=-ΦΦ'Φ=--Φ+Φ'Φ=-'Φ（）②当回程0002s s ϕ'Φ'Φ+Φ+≤≤Φ+Φ+Φ419()312πϕπ≤≤时：20020100202122[)]4[)]4s s h s h h a ϕωυϕω'=Φ+Φ+Φ-'Φ'=-Φ+Φ+Φ-'Φ='Φ（（（4）从动件近休止运动方程在近休止s 'Φ段，即19212πϕπ≤≤时，s 0,0,0a υ===。
哈工大、机械原理大作业、凸轮机构设计20题
Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：能源科学与工程学院班级：1102301设计者：刘平成学号：1110200724指导教师：唐德威设计时间：2013年6月7日凸轮机构设计1．设计题目(1) 凸轮机构运动简图：(2)凸轮机构的原始参数表2-1.凸轮机构原始参数 序号 升程（mm ）升程运动角 升程运动规律升程许用压力角20 110 120° 正弦加速度35°回程运动角回程运动规律 回程许用压力角 远休止角近休止角 90°正弦加速度 65°90°60°（二）凸轮运动方程及相关图像、程序凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图： ○1 凸轮推杆升程、回程方程 πϕπϕϕs)650(πϕ≤≤140)(2=ϕs 511()69πφπ≤≤pi))*5708)/(23.2289)/1.-(sin(2+57083.2289)/1.-(-140(1)(3ϕπϕϕ=s1116()99πφπ≤≤)2914(πϕπ≤≤ 0)(4=ϕs ○2速度方程/2.0944;/2.09440))cos(2-140(1)(1πϕϕ=v 16(2)9πφπ≤≤ 0)(2=ϕv 511()69πφπ≤≤ 708;5708))/1.53.2289)/1.-(cos(2-140(1)(3ϕπϕ=v 1116()99πφπ≤≤ 0)(4=ϕv 16(2)9πφπ≤≤○3加速度方程 .0944^2;/2.0944)/2sin(2280)(1πϕπϕ=a )650(πϕ≤≤0)(2=ϕa 511()69πφπ≤≤08^25708)/1.573.2289)/1.-(sin(2280)(3ϕππϕ=a 1116()99πφπ≤≤ 0)(4=ϕa 16(2)9πφπ≤≤推杆位移、速度、加速度线图matlab编程clear,clcpu=0*pi/180:0.0001:120*pi/180; %升程运动角范围pf=120*pi/180:0.0001:210*pi/180; %远休止角范围pd=210*pi/180:0.0001:300*pi/180; %回程运动角范围pn=300*pi/180:0.0001:2*pi; %近休止角范围h=110e-3; %升程w=10; %凸轮角速度p0=120*pi/180; %升程运动角p01=90*pi/180; %回程运动角ps=90*pi/180; %远休止角%----------推程-----------------------------------------su=h.*(pu./p0-sin(2.*pi.*pu./p0)/(2*pi)); %推杆位移vu=h*w/p0*(1-cos(2*pi*pu./p0)); %推程速度au=2*pi*h*w^2/p0^2*sin(2*pi*pu./p0); %推程加速度%------------远休止----------------------------nf=size(pf);sf=h*ones(nf); %推杆位移vf=zeros(nf); %推程速度af=zeros(nf); %推程加速度%---------------回程------------------------------T=pd-(p0+ps);sd=h/2*(1+cos(pi/p01*T)); %回程位移vd=-pi*h*w/(2*p01)*sin(pi/p01*T); %回程速度ad=-pi^2*h*w^2/(2*p01^2)*cos(pi/p01*T); %回程加速度%--------------------近休止---------------------------------nn=size(pn);sn=zeros(nn); %推杆位移vn=zeros(nn); %推程速度an=zeros(nn); % 推程加速度%------画出推杆位移、速度、加速度线图---------------p=[pu,pf,pd,pn];s=[su,sf,sd,sn];subplot(2,3,1),hold onplot(p,s*1e3,'linewidth',2),xlabel('\phi/rad'),ylabel('s/mm'),grid on,title('推杆位移'),axis([0,2*pi,1.1*min(s)*1e3,1.1*max(s)*1e3]) subplot(2,3,2) v=[vu,vf,vd,vn];plot(p,v,'linewidth',2),xlabel('\phi/rad'),ylabel('v/m/s'),grid on,title('推杆速度'),axis([0,2*pi,1.1*min(v),1.1*max(v)]) subplot(2,3,3) a=[au,af,ad,an];plot(p,a,'linewidth',2),xlabel('\phi/rad'),ylabel('a/m/s^2'),grid on,title('推杆加速度'),axis([0,2*pi,1.1*min(a),1.1*max(a)]) hold off（三）凸轮机构s d ds-ϕ图像及程序代码 %--------------------求ds/d_phi-------------------subplot(2,3,4),plot(v/w*1e3,s*1e3,'linewidth',2),xlabel('ds/d\phi/mm'),ylabel('s/mm'),axis equal,grid on,title('ds/d\phi —s')%---------------------凸轮轴心许用区域--------------------------- alpha_up=35*pi/180; %升程许用压力角 alpha_down=65*pi/180; %回程许用压力角 p1=pi/2-alpha_up; %推程斜率角 p2=alpha_down-pi/2; %回程斜率角 ku=tan(p1); %推程切线斜率 kd=tan(p2); %回程切线斜率 R2=[cos(-p2),-sin(-p2);sin(-p2),cos(-p2)];%推程旋转矩阵 R1=[cos(-p1),-sin(-p1);sin(-p1),cos(-p1)];%推程旋转矩阵 nu=size(pu); for i=1:nu(2)Temp=R1*[vu(i)/w;su(i)];vut(i)=Temp(1); %旋转推程ds/dp-s 曲线 sut(i)=Temp(2); endnd=size(pd); for i=1:nd(2)Temp=R2*[vd(i)/w;sd(i)];vdt(i)=Temp(1); %旋转回程ds/dp-s 曲线 sdt(i)=Temp(2); endfor j=1:nu(2)if sut(j)==min(sut)temu=j; %旋转推程ds/dp-s 曲线后求最低点 end endfor j=1:nd(2)if sdt(j)==min(sdt)temd=j; %旋转回程ds/dp-s曲线后求最低点endendt1=1.2*min(vd/w):0.01:1.2*max(vu/w); %切线定义域t2=min(vd/w)/6:0.01:1.2*max(vu/w);t3=0:0.01:1.2*max(vu/w);s1=ku*(t2-vu(temu)/w)+su(temu); %推程切线s2=kd*(t1-vd(temd)/w)+sd(temd); %回程切线s3=tan(-p1)*t3; %推程起点压力角限制线subplot(2,3,5) %画图hold on,axis equal,grid onplot(v/w*1e3,s*1e3,'linewidth',2)plot(t2*1e3,s1*1e3,'linewidth',1,'color','r')plot(t1*1e3,s2*1e3,'linewidth',1,'color','r')plot(t3*1e3,s3*1e3,'linewidth',1,'color','r')xlabel('ds/d\phi/mm'),ylabel('s/mm'),hold off,title('ds/d\phi—s，轴向许用范围')（四）确定凸轮的基圆半径和偏距、绘制凸轮机圆、偏距圆、理论轮廓曲线---------------画理论廓线图-------------------------e=36e-3;s0=52e-3;r0=sqrt(s0^2+e^2);x=(s0+s).*cos(p)-e.*sin(p);y=(s0+s).*sin(p)+e.*cos(p);x1=r0*cos(p);y1=r0*sin(p);subplot(2,3,6)plot(x*1e3,y*1e3,'linewidth',1),axis equal,grid on,hold on,title('廓线图')plot(x1*1e3,y1*1e3,'linewidth',1,'color','r')%-------------求最小曲率半径-----------------------nx=size(x);nx1=nx(2)-2;dydp=diff(y)./diff(p);%求微分dxdp=diff(x)./diff(p);d2ydp2=diff(dydp)./diff(p(1:nx1+1));d2xdp2=diff(dxdp)./diff(p(1:nx1+1));rho=(dxdp(1:nx1).^2+dydp(1:nx1).^2).^1.5./abs((dxdp(1:nx1).*d2ydp2(1:nx1)-dydp(1:nx1).*d2x dp2(1:nx1)));%理论廓线曲率半径rhomin=min(rho);%最小曲率半径rr=rhomin-3e-3;%----------------实际廓线图----------------X=x(1:nx(2)-1)-rr*dydp./(dxdp.^2+dydp.^2).^0.5;%求实际廓线坐标Y=y(1:nx(2)-1)+rr*dxdp./(dxdp.^2+dydp.^2).^0.5;plot(X*1e3,Y*1e3,'linewidth',2,'color','k')%画实际廓线图Legend('理论廓线','基圆','实际廓线'),axis([1.1*min(x)*1e3,1.1*max(x)*1e3,1.1*min(y)*1e3,1.1*max(y)*1e3])得到基圆半径311mm、偏距36mm。
哈工大机械原理大作业凸轮机构设计第题
哈工大机械原理大作业-凸轮机构设计(第题)————————————————————————————————作者：————————————————————————————————日期：机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：机电学院班级：1208103完成者：xxxxxxx学号：11208103xx指导教师：林琳设计时间：2014.5.2哈尔滨工业大学凸轮机构设计一、设计题目如图所示直动从动件盘形凸轮机构，其原始参数见表，据此设计该凸轮机构。
序号 升程（mm ） 升程运动角（°） 升程运动规律 升程许用压力角（°） 回程运动角（°） 回程运动规律 回程许用压力角（°）远休止角（°） 近休止角（°）3 50 150 正弦加速度 30 100 余弦加速度60 30 80二、凸轮推杆升程、回程运动方程及其线图1 、凸轮推杆升程运动方程(650πϕ≤≤) 升程采用正弦加速度运动规律，故将已知条件mm h 50=，650π=Φ带入正弦加速度运动规律的升程段方程式中得：⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512sin 215650ϕππϕS ； ⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-=512cos 1601ππωv ； ω⎪⎭⎫ ⎝⎛=512sin 14421ϕπωa ； 2、凸轮推杆推程远休止角运动方程（πϕπ≤≤65） mm h s 50==；0==a v ；3、凸轮推杆回程运动方程（914πϕπ≤≤） 回程采用余弦加速度运动规律，故将已知条件mm h 50=，95'0π=Φ，6s π=Φ带入余弦加速度运动规律的回程段方程式中得：⎥⎦⎤⎢⎣⎡-+=)(59cos 125πϕs ； ()πϕω--=59sin451v ； ()πϕω-=59cos 81-a 21；4、凸轮推杆回程近休止角运动方程（πϕπ2914≤≤） 0===a v s ；5、凸轮推杆位移、速度、加速度线图根据以上所列的运动方程，利用matlab 绘制出位移、速度、加速度线图。
哈工大机械原理连杆和凸轮大作业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级组。
哈工大机械原理大作业凸轮结构设计3
仅供个人参考Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮结构设计院系：机电工程学院班级：设计者：学号：指导教师：设计时间：哈尔滨工业大学1、设计题目2、凸轮机构推杆升程、回程运动方程，推杆位移、速度、加速度线图。
（1）推杆各行程运动方程（设定角速度为s s rad /2/1πω==）①从动件推程运动方程(650πϕ≤≤) ⎥⎦⎤⎢⎣⎡⎪⎪⎭⎫ ⎝⎛-=ϕφπφω0012cos 1h v ; 代入数据，可得：②从动件远休程运动方程(πϕπ≤≤65) ③从动件回程运动方程(914πϕπ≤≤) 代入数据，可得:④从动件进休程运动方程(πϕπ2914≤≤) （2）、推杆位移、速度、加速度线图①推杆的位移线图如下②推杆的速度线图如下③推杆的加速度线图如下 3、凸轮机构的dss d ϕ-线图，并依次确定凸轮的基圆半径和偏距凸轮机构的s d ds -ϕ线图如下图所示 由图中范围选定点（-10，-50）为凸轮转轴O 点，则mm r 99.501050220=+=取基圆半径为r 0 =51mm ，偏距e = 10mm 。
4、滚子半径的确定及凸轮理论廓线和实际廓线的绘制由程序计算得凸轮理论轮廓线最小曲率半径mm r 51min = .由滚子半径选择范围∆-<min ρr r ,mm 5~3=∆得到滚子半径mm r r 46≤.又因为凸轮整体尺寸较小，此范围明显过大，故适当减小滚子半径,这里取半径为mm r r 12= .得到图线为：附录1.求位移、速度、加速度的程序（matlab ）function f = tulunh=50;x1=150;t1=30;x2=100;t2=80;w=2*pi;x1=x1*pi/180;x2=x2*pi/180;t1=t1*pi/180;t2=t2*pi/180;%升程x=0:0.001:x1;s = h*(x/x1-sin(2*pi*x/x1)/(2*pi));v = h*w*(1-cos(2*pi*x/x1))/x1;a = 2*pi*h*w*w*sin(2*pi*x/x1)/(x1*x1);subplot(3,1,1),plot(x,s),hold onsubplot(3,1,2),plot(x,v),hold onsubplot(3,1,3),plot(x,a),hold on%远休x = x1:0.001:x1+t1;s = h;v=0;a=0;subplot(3,1,1),plot(x,s),hold onsubplot(3,1,2),plot(x,v),hold onsubplot(3,1,3),plot(x,a),hold on%回程x= x1+t1:0.001:x1+t1+x2;s = h*(1+cos(pi*(x-(x1+t1))/x2))/2;v = -pi*h*w*sin(pi*(x-(x1+t1))/x2)/(2*x2);a = -pi*pi*h*w*w*cos(pi*(x-(x1+t1))/x2)/(2*x2*x2);subplot(3,1,1),plot(x,s),hold onsubplot(3,1,2),plot(x,v),hold onsubplot(3,1,3),plot(x,a),hold on%近休x=x1+t1+x2:0.001:x1+x2+t1+t2;s = 0;v = 0;a = 0;subplot(3,1,1),plot(x,s),xlabel('φ/rad'),ylabel('S/mm'),title('位移-转角图线'),hold onsubplot(3,1,2),plot(x,v),xlabel('φ/rad'),ylabel('v/(mm/s)'),title('速度-转角图线'),hold onsubplot(3,1,3),plot(x,a),xlabel(φ/rad'),ylabel('a/(mm/s^2)'),title('加速度-转角图线'),hold on2.绘制凸轮机构s d ds -ϕ线图 function f= jiyuan;x1=150;t1=30;x2=100;t2=80;h=50;x1=x1*pi/180;x2=x2*pi/180;t1=t1*pi/180;t2=t2*pi/180;x= 0:0.001:150*pi/180;%升程 v/ws = h*(x/x1-sin(2*pi*x/x1)/(2*pi));k =-h*(1-cos(2*pi*x/x1))/x1;plot(k,s,'r'),hold on ;x=180*pi/180:0.001:280*pi/180;%回程 v/ws = h*(1+cos(pi*(x-(x1+t1))/x2))/2;k = pi*h*sin(pi*(x-(x1+t1))/x2)/(2*x2);plot(k,s,'g'),hold on ;%回程切线for i=-11:1:-11;f=@(k)k*tan(pi/6)+i;k =-40:0.1:50;s=f(k);plot(k,s),hold on ;end%升程切线for i=-45:0.2:-45;f=@(k)-k*tan(60*pi/180)+i;k =-40:0.1:50;s=f(k);plot(k,s),hold on ;endgrid onf=@(k)k*tan(50*pi/180);k=-50:0.1:0;s=f(k);plot(k,s),hold onxlabel('ds/d φ');ylabel('s(φ)');title('类速度-位移图线 ');plot(-10,-50,’o ’);3.绘制凸轮轮廓曲线function f= lunkuo;h=50;x1=150;t1=30;x2=100;t2=80;x1=x1.*pi./180;x2=x2.*pi./180;t1=t1.*pi./180;t2=t2.*pi./180;s0=51;e=10;rr=12;%升程x=0:pi/200:150.*pi/180;s = h.*(x./x1-sin(2.*pi.*x./x1)./(2.*pi));X1=(s0+s).*cos(x)-e.*sin(x);Y1=(s0+s).*sin(x)+e.*cos(x);%实际轮廓X11=X1-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y11=Y1-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X1,Y1,'r',X11,Y11,'r'),hold on;%远休x=150.*pi/180:pi/180:180.*pi/180;s=50;X2=(s0+s).*cos(x)-e.*sin(x);Y2=(s0+s).*sin(x)+e.*cos(x);X22=X2-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y22=Y2-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X2,Y2,'g',X22,Y22,'g'),hold on;%回程x=180.*pi/180:pi/180:280.*pi/180;s = h.*(1+cos(pi.*(x-(x1+t1))./x2))./2;X3=(s0+s).*cos(x)-e.*sin(x);Y3=(s0+s).*sin(x)+e.*cos(x);X33=X3-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y33=Y3-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X3,Y3,'k',X33,Y33,'k'),hold on;%近休x=280*pi/180:pi/180:2*pi;s=0;X4=(s0+s).*cos(x)-e.*sin(x);Y4=(s0+s).*sin(x)+e.*cos(x);X44=X4-(rr.*(cos(x).*(s + s0) - e.*sin(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);Y44=Y4-(rr.*(sin(x).*(s + s0) + e.*cos(x)))./((sin(x).*(s + s0) + e.*cos(x)).^2 + (cos(x).*(s + s0) - e.*sin(x)).^2).^(1./2);plot(X4,Y4,'b',X44,Y44,'b'),hold on;x=0:pi/200:2*pi;X4=(s0+s).*cos(x)-e.*sin(x); Y4=(s0+s).*sin(x)+e.*cos(x); plot(X4,Y4,'b');text(-40,90,'理论轮廓线');text(-40,72,'实际轮廓线');text(-5,55,'基圆');grid on;axis equal仅供个人用于学习、研究；不得用于商业用途。
哈工大机械原理大作业——凸轮——2号
哈工大机械原理大作业——凸轮——2号————————————————————————————————作者：————————————————————————————————日期：Harbin Institute of Technology机械原理大作业课程名称：机械原理设计题目：凸轮机构设计一、设计题目(1)凸轮机构运动简图：(2)凸轮机构的原始参数序号升程升程运动角升程运动规律升程许用压力角回程运动角回程运动规律回程许用压力角远休止角近休止角14 90°120°余弦加速度35°90°3-4-5多项式65°80°70°(1) 推杆升程、回程运动方程如下：A.推杆升程方程：设为1rad sω=升程位移为：()()1cos451cos1.52hsπψψψ⎡⎤⎛⎫=-=-⎢⎥⎪Φ⎝⎭⎣⎦23ψπ≤≤升程速度为：()()1100sin67.5sin1.52hvπωπψψωψ⎛⎫==⎪ΦΦ⎝⎭23ψπ≤≤升程加速度为：()()2221100cos101.25cos1.52haπωπψψωψ⎛⎫==⎪ΦΦ⎝⎭23ψπ≤≤B.推杆回程方程：回程位移为：()()345111110156s h T T T ψ⎡⎤=--+⎣⎦1029918ψπ≤≤ 回程速度为：()()22111103012h v T T T ωψ=--+'Φ 1029918ψπ≤≤ 回程加速度为：()()221111260132h a T T T ωψ=--+'Φ 1029918ψπ≤≤其中：()010s T ψ-Φ+Φ='Φ1029918ψπ≤≤ (2) 利用Matlab 绘制推杆位移、速度、加速度线图 A. 推杆位移线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300); x3=linspace(10*pi/9,29*pi/18,300); x4=linspace(29*pi/18,2*pi,300); T1=(x3-10*pi/9)/(pi/2); s1=45*(1-cos(1.5*x1)) s2=90;s3=90*(1-(10*T1.^3-15*T1.^4+6*T1.^5)); s4=0;plot(x1,s1,'r',x2,s2,'r',x3,s3,'r',x4,s4,'r') xlabel('角度ψ/rad'); ylabel('位移s/mm') title('推杆位移线图') gridaxis([0,7,-10,100]) 得到推杆位移线图：B.推杆速度线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);v1=67.5*1*sin(1.5*x1);v2=0;v3=-30*90*1*T1.^2/(pi/2).*(1-2*T1+T1.^2);v4=0;plot(x1,v1,'r',x2,v2,'r',x3,v3,'r',x4,v4,'r') xlabel('角度ψ/rad');ylabel('速度v/(mm/s)')title('推杆速度线图')Grid得到推杆速度线图：C.推杆加速度线图clcclearx1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);a1=101.25*1.^2.*cos(1.5*x1);a2=0;a3=-60.*90.*T1./((pi/2).^2).*(1-3*T1+2*T1.^2); a4=0;plot(x1,a1,'r',x2,a2,'r',x3,a3,'r',x4,a4,'r') xlabel('角度ψ/rad');ylabel('加速度a/')title('推杆加速度线图')Grid得到推杆加速度线图：三、凸轮机构的ds/dψ-s线图，并依次确定凸轮的基圆半径和偏距.1、凸轮机构的ds/dψ--s线图：x1=linspace(0,2*pi/3,300);x2=linspace(2*pi/3,10*pi/9,300);x3=linspace(10*pi/9,29*pi/18,300);x4=linspace(29*pi/18,2*pi,300);T1=(x3-10*pi/9)/(pi/2);s1=45*(1-cos(1.5*x1))s2=90;s3=90*(1-(10*T1.^3-15*T1.^4+6*T1.^5));s4=0;v1=67.5*1*sin(1.5*x1);v2=0;v3=-30*90*1*T1.^2/(pi/2).*(1-2*T1+T1.^2);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曲线')gridaxis([-120,80,-10,100])得到ds/dψ—s曲线：2、确定凸轮的基圆半径和偏距：在dssdϕ-线图中，右侧曲线为升程阶段的类速度-位移图，作直线D t d t与其相切，且与位移轴正方向呈夹角[α1]=350, 故该直线斜率：32sin2=tan5533cos2okϕϕ⨯=⨯通过编程求其角度。
哈工大机械原理大作业2号凸轮
一、如图2-1所示直动从动件盘形凸轮机构，其原始参数见表2-1。
从表2-1中选择一组凸轮机构的原始参数，据此设计该凸轮机构。
表2-1 凸轮机构原始参数升程（mm）升程运动角（）升程运动规律升程许用压力角（）回程运动角（）回程运动规律回程许用压力角（）远休止角（）近休止角（）50 120 余弦加速度35 90 等减等加速65 60 90二、确定凸轮推杆升程、回程运动方程，并绘制推杆位移、速度、加速度线图升程：（0<φ<2π/3）假设w=10;s1=25-25*cos(1.5*φ);v1=37.5*w*sin(1.5*φ);a1=56.25*w^2* cos(1.5φ);远休止：（2π/3<φ<π）s2=50;v2=0;a2=0;回程・1（π<φ<1.25π）s31=50-(100/(π^2/4))*( φ-π).^2;v31=-200*w/(π^2/4)*( φ-π);a31=-200*w^2/(π^2/4);回程・2图2-1（1.25π<φ<1.5π）s32=100/(π^2/4)*(3*π/2-φ).^2; v32=-200*w/(π^2/4)*( 3*π/2-φ); a32=200*w^2/(π^2/4);近休止（1.5π<φ<2π）s4=0;v4=0;a4=0;三、求ds/dφ与s的关系：x=0:(pi/1000):(2*pi/3);s1=25-25*cos(1.5*x);ns1=37.5*sin(1.5*x);y=(2*pi/3):(pi/1000):(pi);s2=50;ns2=0;z1=(pi):(pi/1000):(1.25*pi);s31=50-(100/(pi^2/4))*(z1-(pi)).^2;ns31=-(100/(pi^2/4))*(2*z1-2*pi)z2=(1.25*pi):(pi/1000):(1.5*pi);s32=100/(pi^2/4)*(3*pi/2-z2).^2;ns32=100/(pi^2/4)*(-3*pi+2*z2)m=(1.5*pi):(pi/1000):(2*pi);s4=0;ns4=0;plot(ns1,s1,'b',ns2,s2,'b',ns31,s31,'b',ns32,s32,'b',ns4,s4,'b')凸轮机构的s s dd 线图-80-60-40-200204005101520253035404550凸轮中心应该取在公共区下方。
哈工大—机械原理凸轮大作业
一、题目要求及机构运动简图如图1所示直动从动件盘形凸轮机构。
其原始参数见表1。
图一凸轮运动简图表一凸轮原始参数二、计算流程框图凸轮机构分析建立数学模型位移方程速度方程加速度方程速度线图位移线图加速线图ds/dΨ-s曲线升程压力角回程压力角确定轴向及基圆半径压力角图确定滚子半径实际轮廓理论轮廓轮廓图结束三、建立数学模型1.从动件运动规律方程首先，由于设计凸轮轮廓与凸轮角速度无关，所以不妨设凸轮运动角速度为w = 1rad/s。
（1）推程运动规律(0 < φ < 90°)s=φ）v=φφ）a=φφ）式中：h=65mm，Φ0=π/2（2）远休程运动规律(90°< φ < 190°)s = 65mmv = 0a = 0（3）回程运动规律(190°< φ < 240°)(190°< φ < 196.25°)(196.25°< φ < 233.75°)(233.75°< φ < 240°) 回程运动中的速度和加速度为位移对时间t的倒数：（4）近休程运动规律(240°< φ < 360°)s = 0v = 0a = 02.从动件位移、速度、加速度线图（1）位移线图（2）速度线图（3）加速度线图（4）位移、速度、加速度线图MATLAB源程序%% 已知条件h = 65; %mmphi_0 = 90./180*pi; %radalpha_up_al = 35./180*pi; %升程许用压力角phi_00 = 50./180*pi;alpha_down_al = 70./180*pi; %回程许用压力角phi_s = 100./180*pi;phi_ss = 120./180*pi;w = 1;%% 绘制从动件位移、速度、加速度线图% 推程阶段t_up = 0 : 0.5 : 90;t_up1 = t_up./180*pi;syms t_up1phi_up s_up v_up a_upphi_up = w.*t_up1;s_up = h./2.*(1 - cos(pi.*phi_up./phi_0));v_up = diff(s_up,t_up1);a_up = diff(v_up,t_up1);s_up1 = double(subs(s_up,t_up./180*pi));v_up1 = double(subs(v_up,t_up./180*pi));a_up1 = double(subs(a_up,t_up./180*pi));% 远休程t_s = 90 : 0.5 : (90+100);t_s1 = t_up./180*pi;s_s(1:201) = h;v_s(1:201) = 0;a_s(1:201) = 0;% 回程阶段1t_down1 = (90+100) : 0.5 : (90+100+50/8);t_down11 = t_down1./180*pi;syms t_down11phi_down1s_down1v_down1a_down1phi_down1 = w.*t_down11;s_down1 = h - h./(4+pi).*(pi.*(phi_down1 - phi_0 - phi_s)./phi_00 - ...sin(4.*pi.*(phi_down1 - phi_0 - phi_s)./phi_00)./4);v_down1 = diff(s_down1,t_down11);a_down1 = diff(v_down1,t_down11);s_down11 = double(subs(s_down1,t_down1./180*pi));v_down11 = double(subs(v_down1,t_down1./180*pi));a_down11 = double(subs(a_down1,t_down1./180*pi));% 回程阶段2t_down2 = (90+100+50/8) : 0.5 : (90+100+7*50/8);t_down22 = t_down2./180*pi;syms t_down22phi_down2s_down2v_down2a_down2phi_down2 = w.*t_down22;s_down2 = h - h./(4+pi).*(2+pi.*(phi_down2 - phi_0 - phi_s)./phi_00 - 9.*sin(pi./3 + 4.*pi.*(phi_down2 - phi_0 - phi_s)./(3.*phi_00))./4);v_down2 = diff(s_down2,t_down22);a_down2 = diff(v_down2,t_down22);s_down22 = double(subs(s_down2,t_down2./180*pi));v_down22 = double(subs(v_down2,t_down2./180*pi));a_down22 = double(subs(a_down2,t_down2./180*pi));% 回程阶段3t_down3 = (90+100+7*50/8) : 0.5 : (90+100+50);t_down33 = t_down3./180*pi;syms t_down33phi_down3s_down3v_down3a_down3phi_down3 = w.*t_down33;s_down3 = h - h./(4+pi).*(4+pi.*(phi_down3 - phi_0 - phi_s)./phi_00 - …sin(4.*pi.*(phi_down3 - phi_0 - phi_s)./phi_00)./4);v_down3 = diff(s_down3,t_down33);a_down3 = diff(v_down3,t_down33);s_down33 = double(subs(s_down3,t_down3./180*pi));v_down33 = double(subs(v_down3,t_down3./180*pi));a_down33 = double(subs(a_down3,t_down3./180*pi));% 近休程t_ss = (90+100+50) : 0.5 : 360;s_ss(1:241) = 0;v_ss(1:241) = 0;a_ss(1:241) = 0;% 绘图位移t = [t_up t_s t_down1 t_down2 t_down3 t_ss];phi = w .* t ./ 180 .*pi;s = [s_up1 s_s s_down11 s_down22 s_down33 s_ss];v = [v_up1 v_s v_down11 v_down22 v_down33 v_ss];a = [a_up1 a_s a_down11 a_down22 a_down33 a_ss];figure('Name','从动件位移-时间线图');plot(t,s,'k','linewidth',1.0);grid on;title('从动件位移-时间线图');xlabel('转角\phi / 度');ylabel('位移h/mm');% 绘图速度figure('Name','从动件速度-时间线图');plot(t,v,'k','linewidth',1.0);grid on;title('从动件速度-时间线图');xlabel('转角\phi / 度');ylabel('速度v/mm*s^{-1}');% 绘图加速度figure('Name','从动件加速度-时间线图');plot(t,a,'k','linewidth',1.0);grid on;title('从动件加速度-时间线图');xlabel('转角\phi / 度');ylabel('加速度a/mm*s^{-2}');3.绘制ds/dΦ线图并确定基圆半径和偏距（1）绘制ds/dΦ线图及源程序①MATLAB源程序：%% 绘制ds/dphi-s线图，确定基圆半径和偏距ds_dphi = v ./ w;figure('Name','凸轮ds/dphi - s线图');plot(ds_dphi,s,'k','linewidth',1.5);hold on;axis([-150 150 -70 70]);grid on;title('凸轮ds/dphi - s线图');xlabel('ds/dphi / (mm*s^{-2})');ylabel('s/mm');% 三条临界线x = linspace(-150,150,301);k_up = tan(pi/2 - alpha_up_al);y_up = k_up.*x - 66;plot(x,y_up,'linewidth',1.5);k_down = - tan(pi/2 - alpha_down_al);y_down = k_down.*x - 24.7;plot(x,y_down,'linewidth',1.5);x0 = linspace(0,150,151);k0 = - tan(alpha_up_al);y0 = k0.*x0;plot(x0,y0,'--');% 由图像选取凸轮基圆半径为r0 = sqrt(23^2 + 34^2) = 41 mm,偏距e = 23mm plot(23,-34,'or');r0 = 41;e = 23;plot(linspace(0,23,10),linspace(0,-34,10),'r',linspace(0,23,10),linspace(-34,-34,10),'r',li nspace(23,23,10),linspace(0,-34,10),'r','linewidth',1.0);（2）确定基圆半径和偏距在凸轮机构的ds/d φ-s 线图里再作斜直线Dt-dt 与升程的[ds/d φ-s]曲线相切并使与纵坐标夹角为升程许用压力角[α]，则Dt-dt 线的右下方为选择凸轮轴心的许用区。
哈工大机械原理大作业2凸轮机构设计
机械原理大作业（二）作业名称：凸轮机构设计设计题目：23题院系：班级：设计者：学号：指导教师：设计时刻：哈尔滨工业大学机械设计1.运动分析题目：设计直动从动件盘形凸轮机构，其原始参数见下表2.确信凸轮机构推杆升程、回程运动方程（设定角速度为ω=10 rad/s）升程：0°< Φ < 120°由公式可得：s=60-60*cos(3*Φ/2)；v=90*ω*sin(3*Φ/2)；a=135*ω2 *cos(3*Φ/2)；远停止：120°< Φ < 200°由公式可得：s=120；v=0；a=0；回程：200°< Φ < 290°由公式可得：s=h[1-(10T23-15T24+6T25)]v=(-30hω1/Φ0')T22(1–2T2+T22)a=(-60hω12/Φ0'2)T2(1–3T2+2T22)式中：T2=(Φ-Φ0-Φs)/ Φ0'近停止： 290°< Φ < 360°由公式可得：s=0；v=0；a=0；3.绘制推杆位移、速度、加速度线图(设ω=10rad/s）1) 推拉位移曲线代码：%推杆位移曲线;x=0:(pi/1000):(2*pi/3);s1=60-60*cos(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);s2=120;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;s3=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));m=(29*pi/18):(pi/1000):(2*pi);s4=0;plot(x,s1,'b',y,s2,'b',z,s3,'b',m,s4,'b'); xlabel('角度(rad)');ylabel('行程(mm)');title('推杆位移曲线');grid;2）推杆速度曲线代码：%推杆速度曲线;w=10;x=0:(pi/1000):(2*pi/3);v1=90*w*sin(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);v2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;v3=(-30*120*w/(pi/2))*T2.^2.*(1-2*T2.^2+T2.^2); % v3=-120*w*sin(2*z-20*pi/9);m=(29*pi/18):(pi/1000):(2*pi);v4=0;plot(x,v1,'r',y,v2,'r',z,v3,'r',m,v4,'r'); xlabel('角度(rad)');ylabel('速度(mm/s)');title('推杆速度曲线(w=10rad/s)');grid;3）凸轮推杆加速度曲线代码：%凸轮推杆加速度曲线;w=10;x=0:(pi/1000):(2*pi/3);a1=135*w^2*cos(3*x/2);y=(2*pi/3):(pi/1000):(10*pi/9);a2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;a3=(-60*120*w^2/(pi/2)^2)*T2.*(1-3*T2.^2+2*T2.^2); m=(29*pi/18):(pi/1000):(2*pi);a4=0;plot(x,a1,'m',y,a2,'m',z,a3,'m',m,a4,'m');xlabel('角度(rad)');ylabel('加速度(mm/s^2)');title('凸轮推杆加速度曲线(w=10rad/s)');grid;4）绘制凸轮机构的dd/dd−d线图，并依次确信凸轮的基圆半径和偏距代码：%dd/dd−d线图，确信e，s0;x=0:(pi/1000):(2*pi/3);s1=60-60*cos(1.5*x);ns1=90*sin(1.5*x);y=(2*pi/3):(pi/1000):(10*pi/9);s2=120;ns2=0;z=(10*pi/9):(pi/1000):(29*pi/18);T2=(z-10*pi/9)*2/pi;s3n=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));ns3=-120*10*3*T2.^2+120*15*4*T2.^3-120*6*5*T2.^4 ;m=(29*pi/18):(pi/1000):(2*pi);s4=0;ns4=0;x1=0:pi/36000:pi/2;s1n=60-60*cos(1.5*x1);v1=90*sin(1.5*x1);m1=diff(s1n);%求切线1n1=diff(v1);z=m1./n1;for i=1:length(z);if abs(z(i)+tan(-55*pi/180))<0.001;breakendendb11=s1n(i)-z(i)*v1(i);x1=-300:200;y01=z(i)*x1+b11;%切线1k1=z(i);plot(x1,y01)x3=10*pi/9:pi/36000:14*pi/9;%求切线2s3n=120*(1-(10*T2.^3-15*T2.^4+6*T2.^5));v3=-120*10*3*T2.^2+120*15*4*T2.^3-120*6*5*T2.^4 ;m3=diff(s3n);n3=diff(v3);p=m3./n3;for o=1:length(p);if abs(p(o)-tan(-25*pi/180))<0.01;breakendendo;b33=s3n(o)-p(o)*v3(o);x3=-300:700;y03=p(o)*x3+b33;%切线2plot(x3,y03);sym uv[u,v]=solve('u= 1.4281*v-81.7665','u=-0.4663*v-59.6715');%v=11.66332347972972972972972972973 x%u=-65.110107738597972972972972972973 yplot(ns1,s1,'m',ns2,s2,'b',ns3,s3n,'b',ns4,s4,'b',x1,y01,'g',x3,y03,'g',v,u,'*'); xlabel('ds/d¦µ');ylabel('S');axis([-300,200,-300,300]);title('s0,e 的确信');grid;确信凸轮基圆半径与偏距：偏距e=90mm，d020mm；基圆半径为d0=150mm。
哈尔滨工业大学机械原理大作业凸轮设计
大作业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结构相同，只运动分析过程与其相反。
机械原理大作业-凸轮作业24题(作业试题)
升程
（mm）
升程运
动角
（ ）
升程运
动规律
升程
许用
压力角
（ ）
回程运
动角
（ ）
回程运
动规律
回程
许用
压力角
（ ）
远休
止角
（ ）
近休
止角
（ ）
120
150
正弦加速度
40
100
余弦加速度
60
50
60
2.凸轮推杆升程，回程运动方程及推杆位移，速度，加速度，运动线图
（1）推杆升程，回程方程运动方程如下：
x2=5*pi/6:0.001:pi;
y2=0;
x3=pi:0.001:14*pi/9;
y3=-117*sin(1.8*x3-1.8*pi);
x4=14*pi/9:0.001:2*pi;
y4=0;
plot(x1,y1,x2,y2,x3,y3,x4,y4);
C．推杆加速度线图
Matlab程序：
x1=0:0.001:5*pi/6;
y=60+60*cos(9*(t-10*pi/9)/5);
holdon
plot(x,y,'-r');
t=15*pi/9:0.01:2*pi;
x=0;
y=0;
holdon
plot(x,y,'-r');
gridon
holdoff
如上图所示，在这三条直线所围成的公共许用区域，只要在公共许用区域内选定凸轮轴心O的位置，凸轮基圆半径r0和偏距e就可以确定了。
y3=60+60*cos(9*(x3-10*pi/9)/5);
哈工大机械原理大作业二凸轮机构
void main()
{
float fi=0,fi0=90,h=100,S=0,v=0,a=0,fis=110,fi01=80,fis1=80,T2=0,w1=10,A=0,B=0,v0=0,v01=0,v02=0,f11=30,f31=70,A11=0,A12=0,B11=0,B12=0,A21=0,A22=0,B21=0,C11=0,C21=0,W=0;
{
W=fi*2*PI/360.0;
if(fi<=90)
{
S=h*(W/fi0-(1/(2*PI))*sin((2*PI/fi0)*W));
v=(h*w1/fi0)*(1-cos((2*PI/fi0)*W));
a=(2*PI*h*w1*w1/(fi0*fi0))*sin(2*PI*W/fi0);
}
elseif(fi>=90&&fi<200)
{
A=90;
A=A*2*PI/360.0;
S=h*(A/fi0-(1/(2*PI))*sin((2*PI/fi0)*A));
v=(h*w1/fi0)*(1-cos((2*PI/fi0)*A));
a=(2*PI*h*w1*w1/(fi0*fi0))*sin(2*PI*A/fi0);
}
elseif(fi>=200&&fi<=280)
}
printf("%f\n",a);
}
}
此程序计算推杆位移，速度，加速度线图。
程序二：
#include "stdio.h"
#include "math.h"
哈工大机械原理大作业凸轮 - 黄建青
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机械原理大作业二课程名称：机械原理设计题目：凸轮机构设计院系：能源学院班级： 1302402 设计者：黄建青学号： 1130240222 指导教师：焦映厚陈照波设计时间： 2015年06月23日凸轮机构设计说明书1. 设计题目设计直动从动件盘形凸轮机构，机构运动简图如图1，机构的原始参数如表1所示。
图1 机构运动简图表1 凸轮机构原始参数序号升程升程运动角升程运动规律升程许用压力角回程运动角回程运动规律回程许用压力角远休止角近休止角13 45mm 50°余弦加速度35°90°抛-直-抛70°100°120°计算流程框图：2. 凸轮推杆升程，回程运动方程及推杆位移、速度、加速度线图2.1 确定凸轮机构推杆升程、回程运动方程 设定角速度为ω=1 rad/s(1) 升程：0°<φ<50°由公式可得 )]cos(1[20ϕπΦh s -=计算推程、回程的推杆s 、v 、a运动线图及凸轮dss d ϕ-线图 确定凸轮机构基圆半径和偏距计算曲率半径和压力角，确定滚子半径确定凸轮的理论廓线和实际廓线)s i n (2001ϕπωπΦΦh v =)cos(2020212ϕπωπΦΦh a =(2) 远休止：50°<φ<150° 由公式可得s = 45 v = 0a = 0(3) 回程：150°<φ<240° 由公式得：()()2200002000000022000000,2(1)(1)1,12(1)(1),2(1)s s s s s s s s s Φhn s h ΦΦΦΦΦΦn Φn ΦΦn h n s h ΦΦΦΦΦΦn Φn n ΦΦΦn hn s ΦΦΦΦΦn Φn ϕϕϕϕϕϕ'⎧=---+<≤++⎪'-⎪⎪⎡⎤''-⎪=----++<≤++⎨⎢⎥'-⎣⎦⎪⎪'---⎪'=-++<≤++'-⎪⎩2010000001000200100000n (),(1)(1)n ,(1)(1)n (1),(1)s s s s s s s s Φh v ΦΦΦΦΦΦn Φn ΦΦn h v ΦΦΦΦn Φn n ΦΦΦn h v ΦΦΦΦΦn ΦΦn ωϕϕωϕϕωϕ⎧'=---+<≤++⎪'-⎪⎪''-⎪=-++<≤++⎨'-⎪⎪'---'⎪=--++<≤++''-⎪⎩220100200000220100020n ,(1)(1)0,(1)n ,(1)s s s s s s Φh a ΦΦΦΦn Φn ΦΦn a ΦΦΦΦn n Φn h a ΦΦΦΦΦn Φn ωϕϕωϕ'⎧=-+<≤++⎪'-⎪⎪''-=++<≤++⎨⎪⎪'-'=-++<≤++⎪'-⎩式中 Φ0——推程运动角； Φs ——远休止角；Φ0‘——回程运动角。
(完整word版)哈工大机械原理大作业凸轮DOC
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机械原理大作业二课程名称：机械原理设计题目: 凸轮机构设计院系：班级：设计者：学号：指导教师：哈尔滨工业大学一、设计题目如右图所示直动从动件盘形凸轮机构,选择一组凸轮机构的原始参数,据此设计该凸轮机构。
凸轮机构原始参数序号升程（mm）升程运动角升程运动规律升程许用压力角27130150正弦加速度30°回程运动角回程运动规律回程许用压力角远休止角近休止角100°余弦加速度60°30°80°二. 凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图凸轮推杆升程运动方程:)]512sin(2156[130s ϕππϕ-= )512sin(4.374)]512cos(1[156v 211ϕπϕπωω=-=a％ t 表示转角,s 表示位移t=0：0.01:5*pi/6;%升程阶段s= [(6＊t)/(5*pi ）- 1/（2＊pi ）*sin(12＊t/5）］＊130; hold on plot(t ，s ）; t= 5＊pi/6:0。
01:pi; ％远休止阶段s=130； hold on plot(t,s ）;t=pi ：0.01:14*pi/9；%回程阶段s=65*［1+cos(9*（t-pi ）/5)]； hold on plot(t ，s ）;t=14*pi/9：0.01:2＊pi ；s=0；hold onplot（t,s）；grid onhold off％t表示转角,令ω1=1t=0：0。
01:5*pi/6;%升程阶段v=156*1*［1-cos（12＊t/5）]/pi hold onplot(t，v)；t= 5＊pi/6：0。
01:pi;v=0hold onplot(t，v)；t=pi:0.01:14＊pi/9;％回程阶段v=—117*1*sin（9＊(t—pi)/5) hold onplot（t，v)；t=14＊pi/9:0。
机械原理大作业凸轮..
Harbin Institute of Technology机械原理大作业二课程名称：机械原理设计题目：凸轮结构设计院系：机电工程学院班级：1308108设计者：仲星光学号：1130810816指导教师：林琳设计时间：2015年6月7日一、设计题目如图所示直动从动件盘形凸轮机构，其原始参数见表，据此设计该凸轮机构。
行程（mm)升程运动角（°）升程运动规律升程许用压力角（°）回程运动角（°）回程运动规律回程许用压力角（°）远休止角（°）近休止角（°）6080余弦加速度30 60 摆抛摆60 100 120二．数学计算方法（设计书中所给出的回程运动方程有误）1.理论轮廓和工作轮廓廓线方程正偏置平面凸轮，反转法，小滚子中心在反转运动中的轨迹即为凸轮的理论轮廓，推杆所在直线一直与偏心圆相切，由几何关系可得盘形凸轮理论轮廓线方程为：滚子从动件盘形凸轮的实际轮廓线是以理论轮廓上各点为圆心，以小滚子半径为半径的圆族的包络线，理论廓线B点相对应的实际廓线B'点的方程为：2.压力角3.凸轮轮廓曲率半径由数学分析知，凸轮理论轮廓曲线上任一点的曲率半径的计算公式为：（以上方程来自《基于MatLab语言的机构设计与分析》上海科学技术出版社）三．MatLab程序程序设计流程1.主程序：推杆位clc;clear;h=60;%行程phi01=80;%推程运动角phis1=100;%远休止角phi02=60;%回程运动角phis2=120;%近休止角alpha1=pi/6;%升程许用压力角alpha2=pi/3;%回程许用压力角omiga=1;%凸轮角速度设为1%计算凸轮转角、从动件位移，速度和加速度%绘制从动件位移、速度、加速度线图[psi,s,v,a]=CanShuJiSuan(h,phi01,phis1,phi02,phis2,... omiga,alpha1,alpha2);e=22;r0=57;%由ds/dpsi-s图像确定基圆半径和偏距%计算凸轮轮廓线曲率半径及压力角%绘制凸轮理论轮廓线上的压力角线图和曲率半径图[ang,rou,DxDpsi,DyDpsi]=YaLiJiao_QuLvBanJin( r0,psi,s, v,a,e,omiga );Rr=9;%由曲率半径最小值确定小滚子半径%计算凸轮轮廓曲线%绘制理论和实际轮廓曲线LunKuoXian( r0,psi,s,e,Rr,DxDpsi,DyDpsi );2.余弦加速运动规律子函数function [ s1,v1,a1,psi1 ] = Yuxian( phi01,h,omiga )%计算余弦加速度运动规律psi1=linspace(0,phi01,round(phi01));s1=(h/2).*[1 - cos(pi.*psi1./phi01)];v1=(pi*h*omiga./(2.*(phi01*pi/180))) .*sin(pi.*psi1./phi01);a1=pi^2*h*omiga^2/(2*(phi01*pi/180)^2)*cos(pi.*psi1./p hi01)end3.远休止程子函数function [ s2,v2,a2,psi2 ] = YuanXiu( phi01,phis1,h ) psi2=linspace(phi01+1,phi01+phis1,round(phis1));s2=h*psi2./psi2;v2=0*psi2;a2=0*psi2;end4.摆抛摆子函数function [ s3,v3,a3,psi3 ] =BaiPaoBai( phi01,phis1,phi02,h,omiga )%计算回程摆线-抛物线-摆线运动规律%计算当phi01+phis1<psi<=phi01+phis1+phi02/8时，从动件的位移、速度、加速度psi3_1=linspace(phi01+phis1+phi02/8/100,phi01+phis1+ph i02/8,100);s3_1=h-h*(2*(psi3_1-phi01-phis1)/phi02-sin(4*pi*(psi3_ 1-phi01-phis1)/...phi02)/(2*pi))/(2+pi);v3_1=-2*h*omiga*(1-cos(4*pi*(psi3_1-phi01-phis1)/phi02 ))/((2+pi)*...phi02*pi/180);a3_1=-8*h*pi*omiga^2*sin(4*pi*(psi3_1-phi01-phis1)/phi 02)/((2+pi)*...(phi02*pi/180)^2);%计算当phi01+phis1<psi<=phi01+phis1+3*phi02/8时，从动件的位移、速度、加速度psi3_2=linspace(phi01+phis1+phi02/8+phi02/4/100,phi01+ phis1+3*phi02/8,100);s3_2=h-h*(4*pi*((psi3_2-phi01-phis1).^2)/(phi02^2)-(pi -2)*(psi3_2-phi01...-phis1)/phi02+pi/16-1/(2*pi))/(2+pi);v3_2=-h*omiga*(8*pi*(psi3_2-phi01-phis1)/phi02-pi+2)/( (2+pi)*phi02*pi/180);a3_2=-8*h*pi*omiga^2/((2+pi)*(phi02*pi/180)^2)*psi3_2. /psi3_2;%计算当phi01+phis1<psi<=phi01+phis1+5*phi02/8时，从动件的位移、速度、加速度psi3_3=linspace(phi01+phis1+3*phi02/8+phi02/4/100,phi0 1+phis1+5*phi02/8,100);s3_3=h-h*(2*(pi+1)*(psi3_3-phi01-phis1)/phi02-pi/2-sin (4*pi*(psi3_3-phi01...-phis1)/phi02-pi)/(2*pi))/(2+pi);v3_3=-2*h*omiga*(pi+1-cos(4*pi*(psi3_3-phi01-phis1)/ph i02-pi))/((2+pi)...*phi02*pi/180);a3_3=-8*h*pi*omiga^2*sin(4*pi*(psi3_3-phi01-phis1)/phi 02-pi)/((2+pi)*...(phi02*pi/180)^2);%计算当phi01+phis1<psi<=phi01+phis1+7*phi02/8时，从动件的位移、速度、加速度psi3_4=linspace(phi01+phis1+5*phi02/8+phi02/4/100,phi0 1+phis1+7*phi02/8,100);s3_4=h-h*(-4*pi*((psi3_4-phi01-phis1).^2)/(phi02^2)+(7 *pi+2)*(psi3_4-...phi01-phis1)/phi02-33*pi/16+1/(2*pi))/(2+pi);v3_4=-h*omiga*(-8*pi*(psi3_4-phi01-phis1)/phi02+7*pi+2 )/((2+pi)*phi02*pi/180);a3_4=8*h*pi*omiga^2/((2+pi)*(phi02*pi/180)^2)*psi3_4./ psi3_4;%计算当phi01+phis1<psi<=phi01+phis1+phi02时，从动件的位移、速度、加速度psi3_5=linspace(phi01+phis1+7*phi02/8+phi02/8/100,phi0 1+phis1+phi02,100);s3_5=h-h*(2*(psi3_5-phi01-phis1)/phi02+pi-sin(4*pi*(ps i3_5-phi01-phis1).../phi02-2*pi)/(2*pi))/(2+pi);v3_5=-2*h*omiga*(1-cos(4*pi*(psi3_5-phi01-phis1)/phi02 -2*pi))/((2+pi)...*phi02*pi/180);a3_5=-8*h*pi*omiga^2*sin(4*pi*(psi3_5-phi01-phis1)/phi 02-2*pi)/((2+pi)...*(phi02*pi/180)^2);%输出回程阶段凸轮的转角、从动件的位移速度加速度的数组s3=[s3_1,s3_2,s3_3,s3_4,s3_5];v3=[v3_1,v3_2,v3_3,v3_4,v3_5];a3=[a3_1,a3_2,a3_3,a3_4,a3_5];psi3=[psi3_1,psi3_2,psi3_3,psi3_4,psi3_5];end5.近休止程子函数function [ s4,v4,a4,psi4 ] =JinXiu( phi01,phis1,phi02,phis2 )psi4=linspace(phi01+phis1+phi02+1,phi01+phis1+phis2+ph i02,round(phis2));s4=0*psi4;v4=0*psi4;a4=0*psi4;End6.参数计算及图像绘制[s1,v1,a1,psi1]=YuXian(phi01,h,omiga);%推程[s2,v2,a2,psi2]=YuanXiu(phi01,phis1,h);%远休程[s3,v3,a3,psi3]=BaiPaoBai(phi01,phis1,phi02,h,omiga);%回程[s4,v4,a4,psi4]=JinXiu(phi01,phis1,phi02,phis2);%近休程psi=[psi1,psi2,psi3,psi4];%凸轮转角s=[s1,s2,s3,s4];%从动件位移v=[v1,v2,v3,v4];%从动件速度a=[a1,a2,a3,a4];%从动件加速度%绘制从动件位移、速度、加速度线图figure(1)subplot(3,1,1)plot(psi,s);grid onxlabel('凸轮转角(度)');ylabel('位移（mm）');subplot(3,1,2);plot(psi,v);grid onxlabel('凸轮转角(度)');ylabel('速度（mm/s）');subplot(3,1,3)plot(psi,a);grid onxlabel('凸轮转角(度)');ylabel('加速度(mm/s^2)');%绘制ds/dphi-s线图，并确定凸轮基圆半径和偏距DsDphi=v/omiga;figure(2)plot(DsDphi,s);hold onx1=-150:100;y1=tan(pi/2-alpha1)*(x1-65.26)+22.33; plot(x1,y1);hold ony2=tan(pi/2+alpha2)*(x1+107.6)+22.36; plot(x1,y2);7.轮廓线绘制if Rr==0x=x0;y=y0;elseA=sqrt(DxDpsi.^2+DyDpsi.^2);x=x0+Rr*DyDpsi./A;y=y0-Rr*DxDpsi./A;endfigure(4)plot(r0.*cos(psi.*pi/180),r0.*sin(psi.*pi/180),'-.',.. .x0,y0,'--',x,y,e*cos(psi.*pi/180),e*sin(psi.*pi/180)); grid on;legend('基圆','凸轮理论轮廓','凸轮实际轮廓','偏距圆');axis equalend8.压力角曲线及曲率半径function[ang,rou,DxDpsi,DyDpsi]=YaLiJiao_QuLvBanJin( r0,psi,s, v,a,e,omiga )s0=sqrt(r0.^2-e.^2);rs1=s0+s;ang=abs(atan((v/omiga-e)./rs1))*180/pi;DxDpsi=(v./omiga-e).*sin(psi.*pi/180)+...(s0+s).*cos(psi.*pi/180);DyDpsi=(v./omiga-e).*cos(psi.*pi/180)-...(s0+s).*sin(psi.*pi/180);DDxDpsi=(a./(omiga^2)-(s0+s)).*sin(psi.*pi/180)... +(2*v./omiga-e).*cos(psi.*pi/180);DDyDpsi=(a./(omiga^2)-(s0+s)).*cos(psi.*pi/180)... +(2*v./omiga-e).*sin(psi.*pi/180);A=(DxDpsi.^2+DyDpsi.^2).^1.5;B=abs(DxDpsi.*(DDyDpsi)-DyDpsi.*(DDxDpsi));rou=A./B;%凸轮理论轮廓曲率半径%绘图figure(3)axis onplotyy(psi,rou,psi,ang);%绘制曲率半径图[AX]=plotyy(psi,rou,psi,ang);set(AX(1),'yTick',[0:20:200]);set(AX(2),'yTick',[0:10:100]);end四．程序运行结果与生成图像。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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机械原理大作业说明书课程名称：机械原理设计题目：凸轮机构设计院系：能源学院班级：1002101班设计者：学号：指导教师：赵永强设计时间：6月10日-6月24日哈尔滨工业大学一、设计题目 如右图所示直动从动件盘形凸轮机构，选择一组凸轮机构的原始参数，据此设计该凸轮机构。
序号 升程（mm ）升程运动角 升程运动规律 升程许用压力角27 130150 正弦加速度 30°回程运动角 回程运动规律 回程许用压力角 远休止角近休止角100° 余弦加速度 60° 30° 80°二. 凸轮推杆升程、回程运动方程及推杆位移、速度、加速度线图凸轮推杆升程运动方程：)]512sin(2156[130s ϕππϕ-=)512sin(4.374)]512cos(1[156v 211ϕπϕπωω=-=a% t 表示转角，s 表示位移 t=0:0.01:5*pi/6; %升程阶段s=130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]; hold on plot(t,s); t= 5*pi/6:0.01:pi; %远休止阶段 s=130; hold on plot(t,s);t=pi:0.01:14*pi/9; %回程阶段s=65*[1+cos(9*(t-pi)/5)]; hold onplot(t,s);t=14*pi/9:0.01:2*pi;%近休止阶段s=0;hold onplot(t,s);grid onhold off% t表示转角，令ω1=1 t=0:0.01:5*pi/6;%升程阶段v=156*1*[1-cos(12*t/5)]/pi hold onplot(t,v);t= 5*pi/6:0.01:pi;%远休止阶段v=0t=pi:0.01:14*pi/9;%回程阶段v=-117*1*sin(9*(t-pi)/5) hold onplot(t,v);t=14*pi/9:0.01:2*pi;%近休止阶段v=0hold ont=0:0.001:5*pi/6;a=374.4*sin(12*t/5)/pi; hold onplot(t,a);t=5*pi/6:0.01:pi;a=0;hold onplot(t,a);t=pi:0.001:14*pi/9;a=-210.6*cos(9*(t-pi)/5);t=14*pi/9:0.001:2*pi; a=0;hold on三. 绘制凸轮机构的sdds-ϕ线图% t表示转角，x(横坐标)表示速度ds/dφ，y（纵坐标）表示位移s t=0:0.001:5*pi/6;% 升程阶段x= 156*1*(1-cos(12*t/5))/pi;y= 130*((6*t)/(5*pi)-1/(2*pi)*sin(12*t/5));hold onplot(x,y,'-r');t= 5*pi/6:0.01:pi;%远休止阶段x=0;y=130;hold onplot(x,y,'-r');t=pi:0.001:14*pi/9;x=-117*1*sin(9*(t-pi)/5);y=65*(1+cos(9*(t-pi)/5));hold onplot(x,y,'-r');t=14*pi/9:0.01:2*pi;%近休止阶段x=0;y=0;hold onplot(x,y,'-r');grid onhold off四.按许用压力角确定凸轮基圆半径和偏距1. 求切点转角（1）在图-4中，右侧曲线为升程阶段的类速度-位移图，作直线D t d t与其相切，且位移轴正方向呈夹角[ 1]=300,则切点处的斜率与直线D t d t的斜率相等，因为k Dtdt=tan300,右侧曲线斜率可以表示为,所以，，通过编程求其角度。
编码：求得转角t =1.1123进而求的切点坐标（x,y）=(93.8817, 45.8243)（2）在图-4中，左侧曲线为回程阶段的类速度-位移图，作直线D’t d’t与其相切，且位移轴正方向呈夹角[ 1]=600,则切点处的斜率与直线D’t d’t的斜率相等，因为k Dtdt=tan300同理求得切点坐标（x,y）=( -110.0654, 42.3144)2. 确定直线方程直线D t d t：y =tan(pi/3)(x-93.8817)=45.8243;直线D t’d t’: y =-tan(pi/3)(x+84.3144)+110.0654;3. 绘图确定基圆半径和偏距% 直线Dtdtx=-125:1:150;y= tan(pi/3)*(x-93.8798)+45.8243;hold onplot(x,y);% 直线Dt’dt’x=-125:1:150;y=-tan(pi/6)*(x+110.0654)+34.3144;hold onplot(x,y);%直线Ddx=0:1:150;y=tan(2*pi/3)*x;hold onplot(x,y);t=0:0.001:5*pi/6;x= 156*1*[1-cos(12*t/5)]/pi;y= 130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]; hold onplot(x,y,'-r');t= 5*pi/6:0.01:pi;%远休止阶段x=0;y=130;hold onplot(x,y,'-r');t=pi:0.001:14*pi/9;% 回程阶段x=-117*1*sin(9*(t-pi)/5);y=65*[1+cos(9*(t-pi)/5)];hold onplot(x,y,'-r');t=14*pi/9:0.01:2*pi;%近休止阶段x=0;y=0;hold onplot(x,y,'-r');grid onhold off如图，在这三条直线所围成的公共许用区域，只要在公共许用区域内选定凸轮轴心O的位置，凸轮基圆半径r0和偏距e就可以确定了。
现取轴心位置为x=20，y=-125，则可得偏距e=20，基圆半径=127五.绘制凸轮理论轮廓线编码：%凸轮的理论轮廓，t表示转角，x表示横坐标，y表示纵坐标t=0:0.0001:5*pi/6;x=(125+130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]).*cos(t)-20*sin(t);y=(125+130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]).*sin(t)+20*cos(t);hold onplot(x,y);t= 5*pi/6:0.0001:pi;x=(125+130).*cos(t)- 20*sin(t);y=(125+130).*sin(t)+ 20*cos(t);hold onplot(x,y);t=pi:0.0001:14*pi/9;x=(125+65*[1+cos(9*(t-pi)/5)]).*cos(t)- 20*sin(t);y=(125+65*[1+cos(9*(t-pi)/5)]).*sin(t)+ 20*cos(t);hold onplot(x,y);x=(125).*cos(t)- 20*sin(t);y=(125).*sin(t)+ 20*cos(t);hold onplot(x,y);%基圆t=0:0.001:2*pi;x=20.1074*cos(t);y=20.1074*sin(t);hold onplot(x,y);% 偏心圆t=0:0.001:2*pi;x=13.3509*cos(t);y=13.3509*sin(t);hold onplot(x,y);grid onhold off六、在理论廓线上分别绘出基圆与偏距圆：编码：%凸轮的理论轮廓，t表示转角，x表示横坐标，y表示纵坐标x=(125+130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]).*cos(t)-20*sin(t); y=(125+130*[(6*t)/(5*pi)-1/(2*pi)*sin(12*t/5)]).*sin(t)+20*cos(t); hold onplot(x,y);t= 5*pi/6:0.0001:pi;x=(125+130).*cos(t)- 20*sin(t);y=(125+130).*sin(t)+ 20*cos(t);hold onplot(x,y);t=pi:0.0001:14*pi/9;x=(125+65*[1+cos(9*(t-pi)/5)]).*cos(t)- 20*sin(t);y=(125+65*[1+cos(9*(t-pi)/5)]).*sin(t)+ 20*cos(t);hold onplot(x,y);t= 14*pi/9:0.0001:2*pi;x=(125).*cos(t)- 20*sin(t);y=(125).*sin(t)+ 20*cos(t);hold onplot(x,y);%基圆t=0:0.001:2*pi;x=127*cos(t);y=127*sin(t);hold onplot(x,y);% 偏心圆t=0:0.001:2*pi;x=20*cos(t);y=20*sin(t);hold onplot(x,y);grid onhold off七.确定滚子半径1. 绘制曲率半径图% 凸轮理论轮廓半径，t表示转角，p表示曲率半径,%dxi表示dx/dφ, dyi表示dy/dφ,i=1,2,3,4h=130; %升程t0=pi*5/6; % 升程角t01=pi*5/9; % 回程角ts=pi/6; %远休止角ts1=pi*4/9; %近休止角e=20; %偏距s0=125;% 升程阶段t=linspace(0,pi*5/6,1000);s=h*(t/t0-sin(2*pi*t/t0)/(2*pi));dx1 =(h/t0-h*cos(2*pi*t/t0)).*cos(t)-(s0+s).*sin(t)- e*cos(t); dy1=(h/t0-h*cos(2*pi*t/t0)).*sin(t)+(s0+s).*cos(t)- e*sin(t); p=sqrt(dx1.^2+dy1.^2);hold onplot(t,p);% 远休止阶段t=linspace(pi*5/6,pi,1000);s=h;dx2 =- sin(t).*(s + s0) - e*cos(t);dy2 =cos(t).*(s + s0) - e*sin(t);p=sqrt(dx2.^2+dy2.^2);hold onplot(t,p);% 回程阶段t=linspace(pi,pi*14/9,1000);s=0.5*h*(1+cos(pi*(t-(t0+ts))/t01));dx3 =-0.5*h*pi/(2*t01)*sin((pi/t01)*(t-(t0+ts))).*cos(t)- sin(t).*(s + s0) - e*cos(t); dy3 =-0.5*h*pi/(2*t01)*sin((pi/t01)*(t-(t0+ts))).*sin(t)+ cos(t).*(s + s0) - e*sin(t); p=sqrt(dx3.^2+dy3.^2);hold onplot(t,p);%近休止阶段t=linspace(pi*14/9,pi*2,1000);s=0;dx4 =- sin(t).*(s + s0) - e*cos(t);dy4 =cos(t).*(s + s0) - e*sin(t);p=sqrt(dx4.^2+dy4.^2);hold onplot(t,p);hold offtitle('曲率半径ρ','FontSize',20);grid on八. 绘制实际轮廓线% 凸轮理论轮廓半径，t表示转角，p表示曲率半径,%dxi表示dx/dφ, dyi表示dy/dφ,i=1,2,3,4h=130; %升程t0=pi*5/6; % 升程角t01=pi*5/9; % 回程角ts=pi/6; %远休止角ts1=pi*4/9; %近休止角e=20; %偏距s0=125;rr=10; %滚子半径% 升程阶段t=linspace(0,pi*5/6,1000);s=h*(t/t0-sin(2*pi*t/t0)/(2*pi));x1=(s0+s).*cos(t)-e*sin(t);y1=(s0+s).*sin(t)+e*cos(t);dx1 =(h/t0-h*cos(2*pi*t/t0)).*cos(t)-(s0+s).*sin(t)- e*cos(t); dy1=(h/t0-h*cos(2*pi*t/t0)).*sin(t)+(s0+s).*cos(t)- e*sin(t);X1=x1-rr*dy1./(sqrt(dx1.^2+dy1.^2));Y1=y1+rr*dx1./(sqrt(dx1.^2+dy1.^2));hold onplot(x1,y1);plot(X1,Y1);% 远休止阶段t=linspace(pi*5/6,pi,1000);s=h;x2=(s+s0).*cos(t)-e*sin(t);y2=(s+s0).*sin(t)+e*cos(t);dx2 =- sin(t).*(s + s0) - e*cos(t);dy2 =cos(t).*(s + s0) - e*sin(t);X2=x2-rr*dy2./(sqrt(dx2.^2+dy2.^2));Y2=y2+rr*dx2./(sqrt(dx2.^2+dy2.^2));hold onplot(x2,y2);plot(X2,Y2);% 回程阶段t=linspace(pi,pi*14/9,1000);s=0.5*h*(1+cos(pi*(t-(t0+ts))/t01));x3=(s+s0).*cos(t)-e*sin(t);y3=(s+s0).*sin(t)+e*cos(t);dx3 =-0.5*h*pi/(2*t01)*sin((pi/t01)*(t-(t0+ts))).*cos(t)- sin(t).*(s + s0) - e*cos(t); dy3 =-0.5*h*pi/(2*t01)*sin((pi/t01)*(t-(t0+ts))).*sin(t)+ cos(t).*(s + s0) - e*sin(t);X3=x3-rr*dy3./(sqrt(dx3.^2+dy3.^2));Y3=y3+rr*dx3./(sqrt(dx3.^2+dy3.^2));hold onplot(x3,y3);plot(X3,Y3);%近休止阶段t=linspace(pi*14/9,pi*2,1000);s=0;x4=(s+s0).*cos(t)-e*sin(t);y4=(s+s0).*sin(t)+e*cos(t);dx4 =- sin(t).*(s + s0) - e*cos(t);dy4 =cos(t).*(s + s0) - e*sin(t);X4=x4-rr*dy4./(sqrt(dx4.^2+dy4.^2));Y4=y4+rr*dx4./(sqrt(dx4.^2+dy4.^2));hold onplot(x4,y4);plot(X4,Y4);hold offgrid ontitle('凸轮实际轮廓线','FontSize',20);。