【正文】 (a1/2)/cosf(a1/2)))*$ (2*m(i)*k(i)+b1*m1*m1*j(i))/absf((b1*cosf(i)/$ cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))))elsept(i)=point/c1*i/360,0,(c1/(m1*m1*sinf(a1/2)/cosf(a1/2)))*$ (2*m(i)*k(i)+b1*m1*m1*j(i))/(b1*cosf(i)/$ cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i)))))endif輸出曲線圖形如下： CD桿角加速度曲線其輸出函數值同前利用print函數輸出，這里不再詳述。3滑塊D的速度曲線= l3ω3，而ω3=（l1ω1cos(?1?3)）/s，即滑塊D的速度實際上是CD桿的角速度乘以c1，所以其相關分析內容和前面CD桿角速度分析一樣。這里不再詳述。對其速度公式進行分析，得最大值為：c1*b1*m1/(m2b1)。即得其系數T為c1/(c1*b1*m1/(m2b1))。當i=90時，VD=(c1/(c1*b1*m1/(m2b1)))*c1 *b1*m1/(b1+m2)。當i=270時，VD= (c1/(c1*b1*m1/(m2b1)))*c1*b1*m1/(m2b1)。據上述，程序代碼如下：amp。wcs=obj5ifthen/i==90pt(i)=point/c1*i/360,0,(c1/(c1*b1*m1/(m2b1)))*c1*b1*m1/(b1+m2)elseif/i==270pt(i)=point/c1*i/360,0,(c1/(c1*b1*m1/(m2b1)))*c1*b1*m1/(m2b1)elsept(i)=point/c1*i/360,0,(c1/(c1*b1*m1/(m2b1)))*c1*m1*cosf(iatanf((m2+b1*sinf(i))/(b1*cosf(i))))/(cosf(i)/cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i)))))endif圖形輸出如下： 滑塊D速度曲線其輸出函數值同前利用print函數輸出，這里不再詳述。4滑塊D的加速度曲線，可以知道，滑塊D的加速度實際上就是CD桿的角加速度乘以一個CD桿的長度，而CD桿的長度是個定量，所以其滑塊D的加速度曲線實際上和CD桿的曲線走勢一樣。分析過程完全一樣，這里亦不再詳述。對于其最大值也是CD桿的最大角加速度乘以CD桿長。程序代碼如下:ifthen/i==270pt(i)=point/c1*i/360,0,0elseif/i==90pt(i)=point/c1*i/360,0,0elseif/i270pt(i)=point/c1*i/360,0,(c1/(c1*m1*m1*sinf(a1/2)/cosf(a1/2)))*c1*absf((2*m(i)*k(i)+b1*m1*m1*j(i))/(b1*cosf(i)/cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))))elseif/i90,if/i270pt(i)=point/c1*i/360,0,(c1/(c1*m1*m1*sinf(a1/2)/cosf(a1/2)))*c1*(2*m(i)*k(i)+b1*m1*m1*j(i))/absf((b1*cosf(i)/cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))))elsept(i)=point/c1*i/360,0,(c1/(c1*m1*m1*sinf(a1/2)/cosf(a1/2)))*c1*(2*m(i)*k(i)+b1*m1*m1*j(i))/(b1*cosf(i)/cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i)))))，其輸出函數值同前利用print函數輸出，這里不再詳述。5E點的X坐標經過數學模型的分析可知，E點的X坐標值即是D點的X坐標值，所以E點的X坐標值的變化即是D點的變化。即：c1*cosf(atanf((m2+b1*sinf(i)) /(b1 *cosf(i))))。對上式分析，最大值在CD桿的極限位置，得：c1/(c1*sinf(a1/2))在區(qū)間[90,270]上為負值。所以要分區(qū)間分析。程序代碼如下： 滑塊D加速度amp。wcs=obj5ifthen/i=270point/c1*i/360,0,(c1/(c1*sinf(a1/2)))*c1*$ cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))elseif/i90,if/i270pt(i)=point/c1*i/360,0,(c1/(c1*sinf(a1/2)))*$ c1*cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))elsept(i)=point/c1*i/360,0,(c1/(c1*sinf(a1/2)))*$ c1*cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))endif程序輸出圖形如下： E點X坐標變化曲線其輸出函數值同前利用print函數輸出，這里不再詳述。6E點的速度曲線：VE=l3ω3sin?3即E的速度就是D速度的X軸分量。即此式多了一個變量?3，首先對E（即輸出桿件）的運動過程的物理模型進行分析，在區(qū)間[(180+a1/2),(360a1/2)]中為負向轉動，在[0,(180+a1/2)]和[(360a1/2),360]中為正向轉動。首先分析其速度最大值點，由于刨床機構的急回特性，輸出刨刀機構需要有此特性，再對其數學公式進行分析。得到其最大值點為i=270，此時最大速度值為：c1*b1*m1/(m2b1)，可得系數T為c1/(c1*b1*m1/(m2b1))，當i=90時，(c1/(c1*b1*m1/(m2b1)))* (c1*b1*m1/(m2b1))，當i=270時，(c1/(c1*b1*m1/(m2b1)))*(c1*b1*m1/(b1+m2))。程序代碼如下：amp。wcs=obj5ifthen/i==270pt(i)=point/c1*i/360,0,(c1/(c1*b1*m1/(m2b1)))*(c1*b1*m1/(m2b1))elseif/i==90pt(i)=point/c1*i/360,0,(c1/(c1*b1*m1/(m2b1)))*(c1*b1*m1/(b1+m2))elseif/i(360a1/2)pt(i)=point/c1*i/360,0,(c1/(c1*b1*m1/(m2b1)))*absf(sinf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))*c1*m1*cosf(iatanf((m2+b1*sinf(i))/(b1*cosf(i))))/(cosf(i)/cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))))elseif/i(180+a1/2),if/i=(360a1/2)pt(i)=point/c1*i/360,0,(c1/(c1*b1*m1/(m2b1)))*absf(sinf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))*c1*m1*cosf(iatanf((m2+b1*sinf(i))/(b1*cosf(i))))/(cosf(i)/cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))))elsept(i)=point/c1*i/360,0,(c1/(c1*b1*m1/(m2b1)))*absf(sinf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))*c1*m1*cosf(iatanf((m2+b1*sinf(i))/(b1*cosf(i))))/(cosf(i)/cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))))endif圖形如下： E點速度曲線其輸出函數值同前利用print函數輸出，這里不再詳述。7E點的加速度曲線= H/2/（sinβ）α3 sin ? 3即輸出桿件的加速度實際是D的加速度的水平分量。對前面的E點的速度曲線分析可以得其加速度曲線，可知加速度曲線可分區(qū)間[0,90],[90,270],[270,360]分析。在90和270兩個點加速度皆為0。加速度最大值在360a1/2，和180+a1/2兩點。此時最大值為：c1*m1*m1*sinf(a1/2)即得系數：c1/(c1*m1*m1*sinf(a1/2))。此時X軸的密度需有所調整：c1*i/250帶入相關參數并輸出源代碼如下：ifthen/i==270pt(i)=point/c1*i/250,0,0elseif/i==90pt(i)=point/c1*i/250,0,0elseif/i270pt(i)=point/c1*i/250,0,(c1/(c1*m1*m1*sinf(a1/2)))*(c1*sinf(atanf((m2+b1*sinf(i))/(b1*cosf(i)))))*absf((2*m(i)*k(i)+b1*m1*m1*j(i))/(b1*cosf(i)/cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))))elseif/i90,if/i270pt(i)=point/c1*i/250,0,(c1/(c1*m1*m1*sinf(a1/2)))*(c1*sinf(atanf((m2+b1*sinf(i))/(b1*cosf(i)))))*(2*m(i)*k(i)+b1*m1*m1*j(i))/absf((b1*cosf(i)/cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i))))))elsept(i)=point/c1*i/250,0,(c1/(c1*m1*m1*sinf(a1/2)))*(c1*sinf(atanf((m2+b1*sinf(i))/(b1*cosf(i)))))*(2*m(i)*k(i)+b1*m1*m1*j(i))/(b1*cosf(i)/cosf(atanf((m2+b1*sinf(i))/(b1*cosf(i)))))endif輸出圖形如下： E點加速度曲線其輸出函數值同前利用print函數輸出，這里不再詳述。 本章小結本章為整個設計中最重要的部分，也是整個設計耗時最多的階段。本章詳細介紹了刨床機構運動過程的實現，運動曲線的輸出過程，詳細介紹了各部分代碼的作用。本章首先介紹本設計的設計思想，接著開始運動初始準本工作：建立坐標系及其分析過程。然后給出了本次設計思想的流程圖。從流程圖中可以看出設計的主要步驟。接著，就是確定了運動方式。最后開始運動程序的建立，運動曲線的輸出。在建立刨床機構的運動過程中，詳細介紹了每個運動部件的運動的實現方法，并且詳細的介紹了各代碼的作用。在運動曲線輸出階段，首先說明了輸出曲線的對后續(xù)的動力學分析的意義作用，接著詳細闡述了每個公式的得到過程，有詳細的計算過程。最后對每個參數的輸出都分別作了詳細的講解，并且輸出了各運動參數的輸出曲線圖。第五章 優(yōu)化并完善設計軟件 查詢功能的實現 建立查詢窗口f10:choose/39。分析結果查詢39。,39。查詢39。,39。不查詢39。,respjump/f10:,trm:,f20:,f30:,respf20:param/39。輸入查詢角度39。,39。角度39。,m7,resp 設計思想在此設計過程中，需要為產品使用者提供查詢功能，以便更加完善產品的功能，也便于使用者查詢相關數據。經過分析，查詢功能實現主要有以下兩點：第一就是當使用者選擇“查詢”時，在函數輸出曲線上，對相關點給以標示，第二就是在print輸出函數值的TXT窗口中有詳細描述，便于查詢。對于第一點，本設計采用在相關點作一垂線至X軸，對于第二點，就是在函數輸出值的窗口中有相應角度于其值對應。當使用者輸出查詢角度后，需要使刨床桿件轉到相應角度，當點擊“繼續(xù)查詢”后，需使其回歸到原點，便于下一次查詢，直至點擊不查詢?yōu)橹埂．斒褂谜哌x擇“重新運動分析”時，需要清除屏幕上的所有曲線即實體。并跳出原始參數輸入窗口，便于使用者選擇其他桿件分析。 實現查詢功能在此查詢功能中的桿件轉動過程和前一章的運動分析過程一樣，其設計思路一樣，這里不再詳述，只輸出代碼。1曲柄的移動amp。wcs=obj3mat3(1..12)=matrix/zxrot,m7obj(23+m7)=transf/mat3,obj(23),move2CD桿的轉動ifthen/m7270n(m7)=atanf((m2+b1*sinf(m7))/(b1*cosf(m7)))elseif/m790,if/m7270n(m7)=180absf(atanf((m2+b1*sinf(m7))/(b1*cosf(m7))))elseif/m7==90n(m7)=90elseif/m7==270n(m7)=