平台印刷机-机械原理课程设计报告

----.可修编.机械原理课程设计设计题目平台印刷机主传动机构运动简图设计学院专业机械设计制造及其自动化指导老师班级学号设计者CONTENTS一、设计任务3二、机械运动方案的选择及其论证31．平台印刷机主要机构及功能3（1）主要机构：3(2)各机构功能：42．各机构形态学矩阵43．设计思路概述4三、原始数据51．平台印刷机设计数据：52．其它数据：5四、设计方法说明及计算公式61．曲柄滑块机构综合分析：6(1)机构的运动几何关系：6(2)参数选择：72．双曲柄机构的运动分析：8(1)曲柄滑块位移计算Ψ：8(2)由Ψ1求Ψ3：83．曲柄滑块机构的位置分析：94．函数平方逼近法设计双曲柄机构：105．凸轮机构的设计:12(1)凸轮机构从动件运动规律的确定12(2)绘制补偿凸轮轮廓13五、子程序名称：13六、编程框图及主程序14七、结果分析16八、心得感悟16九、参考资料17十．附录程序清单及运行结果171.程序部分17(1)C语言程序17(2)程序清单二（全局变量应用——参考样板——2010.7.8）20（3）VC屏幕截图222、结果记录及处理27（1）写入CC.txt的数据27（2）EXCEL表格及图表30一、设计任务工作原理：平台印刷机的工作过程由输纸，着墨，压印和收纸四部分组成，主运动是压印，由卷有空白纸X的滚筒与镶着铅字的版台之间纯滚动来完成。滚筒与版台表面之间的滑动会造成字迹模糊，是不允许的。因此，对运动的主要要求是：其一，版台的移动速度严格等于滚筒表面的圆周速度；其二，为了提高生产率，要求版台的运动有急回特性。有一台电动机驱动。需设计满足上述两个要求的传动机构。执行件的运动为滚筒连续转动和版台往返移动。（图1）图1二、机械运动方案的选择及其论证1．平台印刷机主要机构及功能（1）主要机构：1)传动机构I——从电动机到版台的运动链；2)传动机构II——从电动机到印刷滚筒的运动链；3)位移补偿机构。(2)各机构功能：1) 传动机构I——把电动机的旋转运动转化为版台的直线移动；2) 传动机构II——把电动机的旋转运动转化为印刷滚筒的旋转运动；3) 位移补偿机构——把凸轮的旋转运动转化为版台的直线移动。2．各机构形态学矩阵功能元实现功能的解（各个方案）传动机构I摩擦传动螺旋机构齿轮齿条凸轮机构曲柄滑块机构组合机构传动机构带传动机构双曲柄机构齿轮机构补偿机构凸轮机构3．设计思路概述由电动机到印刷机滚筒这一条运动链的两端皆为定轴转动，适宜于传递旋转运动的机构有四杆机构，齿轮机构，涡轮机构，根据平稳性和简易度，应选齿轮机构。由电动机到版台间的运动链须将转动变为移动，转动变移动的传动方式很多，常见的有摩檫传动，齿轮齿条传动，螺旋机构传动，凸轮机构，曲柄滑块机构以及组合机构等等。但前三种机构不具有急回特性。由于曲柄滑块比凸轮运动精确性高所以选择曲柄滑块机构两条传动链的终端——滚筒和版台的瞬时线速度须相等。当传动机构1和2本身不能满足时，可以另设补偿机构实现它。用速度补偿来设计其运动曲线可以选择凸轮机构.图2所示的方案可作为机构1的参考方案，双曲柄机构AoABBo与曲柄滑块机构BoCD串联，将曲柄AoA的转动改变为D点的往复移动。当齿条6固定不动时，中心为D的行星齿轮5将带动齿条7移动。齿条7固定在印刷版台的下面。齿条6的位置由补偿凸轮3通过磙子从动件控制。该凸轮和从动曲柄BoB为同一构件。而主动主动曲柄AoA与滚筒的转动同步。当齿条7（即版台）的工作移动速度偏离要求时，可通过凸轮廓线控制齿条6补充移动来补偿。另一传动机构2是一个双曲柄机构，它用于带动滚筒的转动，其各杆的尺寸为h1=145.0mm,h2=178.0mm,h3=175.0mm,h4=65.5mm.传动机构1，2之间通过一对齿轮传动建立其运动关系，齿轮传动的参数为I=1,m=4mm,z=105,见图4所示。三、原始数据1．平台印刷机设计数据：印刷纸X最大幅面（）415590生产率（X/小时）4400~4500滚筒直径()360版台往复行程长度（）795电动机功率（kw）3转速（r.p.m）14502．其它数据：传动机构I： 机架长55.0mm传动机构II：h1=145.0mm，h2=178.0mm，h3=175.0mm，h4=65.5mm齿轮传动的参数：i=1，m=4mm，z=105四、设计方法说明及计算公式分析其运动特点,准备运动的分析计算公式和C语言程序.计算数据,绘制运动关系图.进行机构设计1．曲柄滑块机构综合分析：(1)机构的运动几何关系：如图3所示普通的偏置曲柄滑块机构，标记曲柄长为r，连杆长l，偏心距e：图3引入参量λ=l/r,δ=e/l;则有cosα=(1-sin2α)1/2=(1-[(l/λ)sinΦ-δ]2)1/2,当（l/λ）和δ均较小时，将该示展开，得其近似值。cosα=1-sin2Φ/(2λ2)+δ/λ*sinΦ-δ2/2D点的位置方程为：X=r*cosΦ+l*cosα=r*cosΦ+l*(1-sinΦ2/(2λ2)+δ/λ*sinΦ-δ2/2)D点的位移（由D1点算起）为：S=((r+l)2-e2)1/2-X速度和加速度方程分别为：V=dS/dt=rω*(sinΦ+1/(2λ)*sin2Φ-δ*cosΦ)A=dV/dt=d2S/dt2=rω2(cosΦ+1/λ*cos2Φ+δ*sinΦ)+rε*(sinΦ+1/(2λ)*sin2Φ-δ*cosΦ)滑块的行程长度为：H=D1D2=((块机构的H大于对心的滑块行程长度（H=2r）.滑块的行程速比系数；K=(180+θ)/(180-θ)=(180+sin((δ*λ)/(λ-1))-sin((δ*λ)/(λ+1)))/(180-sin((δ*λ)/(λ-1))+sin((δ*λ)/(λ+1)))机构传动角分为工作行程和回程：工作行程的最小传动角γa=cos((r-e)/l)=cos(1/λ-δ)回程的最小传动角γa=cos((r+e)/l)=cos(r/λ+δ)(2)参数选择：a.Ha根据设计原始参数计算b.λa由以上式子可知，λ增大对γ有利。但是尺寸将跟随λ增大。在此取为3.6.c.δa由上式可知，在H和λ即定情况下，δ影响曲柄长r,印刷机中，δ取0.32.d.在确定λ和δ时要照顾eo.e=0则无急回效果，K随e增。同时，由(2—5)式知，当H即定时，e增大则必使r减少。l+r)2-e2)（1/2）-((l-r)2-e2)1/2=r(2+(δ2λ2)/(λ2-1)2．双曲柄机构的运动分析：图4图4所示，为了改善版台的速度特性，机构1设计成串联组合机构。第一级为双曲柄机构ABCD。其主动曲柄DC用一级速比I=1的圆柱齿轮与机构2（亦为双曲柄机构）的原动曲柄h相连接，即Φ=Ψ机构2的从动曲柄h则与滚筒同步转动。(1)曲柄滑块位移计算Ψ：Ψ=2*S/Rcir式中S—滑块位移；Rcir—滚筒半径。(2)由Ψ1求Ψ3：如图5所示，几何关系为：h4^2+h3^2+h1^2-h2^2-2h1*h4cosΨ—1-2h1*h3cosΨ1cosΨ3-2h1*h3sinΨ1sinΨ3=0A*sinΨ3+B*cosΨ3-C=0式中A=sinΨ1B=cosΨ1-K1C=K2-K3cosΨ1K1=h4/h1K2=(h1^2-h2^2+h3^2+h4^2)/(2h1*h3)K3=h4/h3可解得Φ3=2arctg((A+(A^2+B^2-C^2)^2)/(B+C))Φ3=2arctg((A-(A^2+B^2-C^2)图5 图63．曲柄滑块机构的位置分析：图5 图6由滑块位移求对应的曲柄转角Φ1杆长符号如图6所示，则有So=((r+l)^2-e^2)^(1/2)X=So-SX=r*cosΦ+lcosαe=r*sinΦ-lsinα移项，作平方和，消去α得X^2+r^2+e^2-2X*rcosΦ-2ersinΦ=l^22esinΦ+2X*cosΦ-(X^2+r^2+e^2-l^2)/r=0改写成A*sinΦ+B*cosΦ-C=0（1）其中A=2e,B=2X,C=(X^2+r^2+e^2-l^2)/r（2）式的解为：Φ=2*arctg((A+/-(A^2+B^2-C^2)^(1/2))/(B+C))考虑（2）式，并令Z=C/2=(X^2+r^2+e^2-l^2)/(2r)代入后得：Φ=2*arctg((e+/-(e^2+X^2-Z^2)^(1/2))/(X+Z))对照图4，可知：Φ1=π/2-Φ-β4．函数平方逼近法设计双曲柄机构：图7图7坐标设置如图7所示，当αo=Φo=00,l4=l0机构的位置方程l2+l32+l12-l22+2l3*cosΦ-2l1*cosα-2l1*l3*cos(Φ-α)=0除以2l1,得l3/l1*cosΦ+P1*cos(Φ-α)+(1+l32+l12-l22)/(2*l1)-cosα=0其中P0=l3/l1P1=-l3P2=(1+l32+l12-l22)/(2*l1)上述方程中只有3个待定参数，只能满足3个精确点（Φi，αi，I=1,2,3）.精确点外的机构函数与目标函数存在结构误差。Δy=P0*cosΦ+P1*cos(Φ-α)+P2-cosα!=0所谓平方逼近法设计，就是以结构误差的均方根值最小为目的，做逼近函数机构的设计。对于目标函数为连续函数时，我们建立I=∫[P*cosΦ+P*cos(Φ-α)+P-cosα]^2*dα而对于[α0，αm]中的若干离散点，则构造I=∑[P*cosΦi+P1*cos(Φi-αi)+P2-cosαi]^2可改写成I=∑[Po*fo(αi)+P1*f1(αi)+P2*f2(αi)-F(αi)]^2其中，fo(αi)=cosΦif1(αi)=cos(Φi-αi)f2(αi)=1F(αi)=cosαi为求I的极小值。令I/Pi=0,I=0,1,2,得∑[Po*fo(αi)+P1*f1(αi)-F(αi)]*fo(αi)=0∑[Po*fo(αi)+P1(αi)+P2*f2(αi)-F(αi)*f1(αi)=0∑[Po*fo(αi)+P1*f1(αi)+P2*f2(αi)-F(αi)]*f2(αi)=0令C00=∑f0(αi)*f0(αi)C01=C10=∑f0(αi)*f1(αi)C02=C20=∑f0(αi)*f2(αi)C11=∑f1(αi)*f(αi)C12=C21=∑f(αi)*f(αi)C22=∑f(αi)*f(αi)或缩写成Ckl=∑fk()*f()K=0,1,2;L=0,1,2;γ0=∑F(αi)*f0(αi)γ1=∑F(αi)*f1(αi)γ2=∑F(αi)*f2(αi)则有C00P0+C01P1+C02P2=γ0C10P0+C11P1+C12P2=γ1C20P0+C21P1+C22P2=γ2从中可以解出p0,p1,p2,然后由上式计算相对杆长。5．凸轮机构的设计:由于印刷机的版台和滚筒各由一套独立的运动传动系统驱动.为了保证印刷质量.在压印阶段.滚筒表面点的线速度必须和版台移动速度保持相等.在设计时.应尽可能满足这一要求.如果设计的结果不能完全满足这一要求.即在压印区滚筒表面点的线速度与版台移动速度尚有一定的差别.则采用凸轨机构进行运动补偿.如图4所示.若版台由双曲柄六杆机构和齿轮齿条串联而成的传动系统所驱动.则可将下齿条做成活动齿条.在曲柄BE上固连一补偿凸轮.此补偿凸轮通过滚子直动推杆推动活动的下齿条.经过齿轮与上齿条的传动.给版台附加一个运动.使版台移动速度与滚筒表面点的线速度完全一致.这样凸轮机构起到了运动补偿作用.由于凸轮轮廓形状和凸轮机构的工作性能取决于从动件运动规律.显然在印刷机的设计中.此凸轮机构从动件运动规律完全取决于所需补偿量的变化规律.(1) 凸轮机构从动件运动规律的确定1) 设版台传动系统未安装补偿凸轮机构.且下齿条为固定齿条.编主程序对此传动系统进行运动分析.可求出版台的位移曲线Sp-φ(φ为曲柄AB的转角,Sp为版台位移)及速度曲线Vp-φ:2) 编主程序对滚筒的传动系统进行运动分析.求出滚筒表面点的位移曲线Sc-φ(Sc为滚筒表面点转过的弧长)及速度曲线Vc-φ:3) 按上述速度曲线与位移曲线.选定同步区(即压印区).同步区的始点一般应为同速点(即版台运动速度与滚筒表面点速度相同的点).同步区应根据压印行程要求(即印刷纸X的长度)选在速度变化较小的区域.即∣△V∣=∣Vp-Vc∣较小的区域.4) 设同步区始点为O点.从O点开始.将同步区分成n段.得到(n＋1)个分点.则可依次求出:△Sp(i)=Sp(i)－Sp(o)i=(1,2,3,……，n)△Sc(i)=Sc(i)－Sc(o)i=(1,2,3,……，n)△S(i)=Sp(i)－Sc(i)i=(1,2,3,……，n)然后.在△S-φ坐标中描出△S(i)(i=1,2,3,……，n)各点.并将它们光滑联接.即得到曲线.此曲线就是补偿凸轮从动件在压印区的位移曲线.也就是从动件的推程运动规律.另外再选择合适的过渡曲线作为从动件的回程运动规律.(2) 绘制补偿凸轮轮廓选择合适的从动件结构型式并考虑其它因素来选择凸轮机构的基本参数.用图解法设计凸轮轮廓曲线.五、子程序名称：《1》源程序及改写部分1)voidCFI2)voidCPSI3)voidCRLE4)voidSLNPD5)voidSQU6)voidEQU7)voidPOS18)voidPOS29)voidPOS310)voidVEL111)voidACCE12)voidWRT1《2》版台位移voidPS六、编程框图及主程序开始开始定义定义λ,δ及最大行程调用CRLE确定曲柄滑块R,，E，调用CRLE确定曲柄滑块R,，E，L判断最小传动角判断最小传动角定义定义ψ10，φ30，H[5],S[11]及循环变量i=1i<=10i<=10S[i]=S[i]=130+((370-130)*(i-1))/9;真假调用CFI确定调用CFI确定φ3i调用CPSI确定调用CPSI确定φ1i定义DL=55.0(即l4定义DL=55.0(即l4)，调用SQU确定双曲柄机构定义所需变量及已知数据定义所需变量及已知数据起始FAI=0起始FAI=0FAI<=360假FAI<=360真调用POS3确定曲柄滑块位置角FAI[3]调用POS3确定曲柄滑块位置角FAI[3]SET[1]=(IFA+14.84722+90)*3.14159SET[1]=(IFA+14.84722+90)*3.14159调用POS2确定机构I位置角SET[4]调用POS2确定机构I位置角SET[4]确定SSET[3]确定SSET[3]调用POS1确定机构II位置角SSET调用POS1确定机构II位置角SSET确定ALP的值，确定SSET[1]确定ALP的值，确定SSET[1]调用VEL调用VEL1求出版台、滚筒的速度Vp、Vc调用SCI求出滚筒转过的弧长，调用PS求版台位移，并求出位移差调用WRT1输出结果调用WRT1输出结果结束结束七、结果分析一、此次设计，充分利用了各类机构的运动特性：平面四杆机构的急回特性，能够节省回程时间；利用齿条的将圆周运动变换为直线运动的性质，推动版台作直线运动；利用凸轮机构响应快速，适用于自动机和自动控制装置的特性，实现运动补偿。二、设计的机构基本能满足运动要求。八、心得感悟这次机械原理课程设计，题目在多个指导书上都见过。在这两周的时间里，我过得是相当紧X忙碌。一方面要进行课程设计，另一方面还要准备两门课程考试。人在这些天里感到十分忙碌，经常忙到凌晨2点，第二天早上6点30又要起床。紧X忙碌的两星期里，我付出了许多，有所失去，也有所得。有所失去，应该算是课程考试成绩的不XX想吧。由于首次接触课程设计，感觉颇为陌生，课程设计占用了大部分时间，因此课程考试成绩不太理想。但我更有所得。其一，通过这次真刀实枪的演练，我学到的知识终于成功投入了实践，特别是投入了具一定工程意义的实践。以前我们总抱怨知识学了未能用，这次有了用场。比如C语言、EXCEL等。其二，我巩固了相关知识，也学到了新的知识。比如，凸轮机构的反转法作图，四杆机构的函数平方逼近法，等等。其三，我对设计有了初步的了解，这对以后我们的发展将是很有益的。最后，这次课程设计让我重新认识了自己。比如，某些C程序段我坚持认为没问题，但运行总出错，这时候上网搜索答案，我才明白；我想用MATLAB进行曲线拟合，结果却忘了函数，又因为数组输入过于复杂，因此放弃。这些都说明了我的基本功还不够，我高看了自己。此后必须脚踏实地，不能再眼高手低了。九、参考资料十．附录程序清单及运行结果1.程序部分(1)C语言程序#include<stdio.h>#include<math.h>/*函数声明部分*/voidCFI(intN,floatGAMA10,floatFI0,float*H,float*S,float*FI,FILE*fp1);voidCPSI(intN,floatBATA,float*TL,float*S,float*PSI,FILE*fp1);voidCRLE(floatALAMT,floatER,floatHD0,floatHD,float*TL,FILE*fp1);voidSQU(intN,float*FI,float*PSI,floatDL,float*TL1,FILE*fp1);voidSLNPD(floatA[4][4],float*B,FILE*fp1);voidPOS1(float*TL,float*SET,intI1,intI2);voidPOS2(float*TL,float*SET,intI3,intI4);voidPOS3(float*TL,float*FAI);voidVEL(floatALP,float*VF,float*VCIR,float*TL1,float*TL2,float*A,float*TLL,float*SET,float*SSET,float*FAI);voidEQU(float*SE,floatZ,floatE,floatF,floatG,intI);voidSCI(intIFA,floatSSETR,float*SCIR);voidACCE(float*L1,float*SET,float*DSET,float*DDSET1);voidCheck(float*TL,FILE*fp1);voidVEL1(floatFAI,float*TL1,float*TL2,float*A,float*TLL,float*VF,float*VCIR,float*SCIR,float*SP);voidWRT1(floatFAI,floatSP,floatSCIR,floatVP,floatVCIR,FILE*fp);/*主函数部分*/voidmain(){ FILE*fp=fopen("F:\\CC.txt","a+");constfloatPI=3.1415926;/*宏定义*/ floatALAMT=3.55,ER=0.305,HDO=795/2.0;/*曲柄滑块机构的设计*/ floatHD,TL[4]; CRLE(ALAMT,ER,HDO,HD,TL,fp); floatTHITA,cd1,cd2,K; THITA=asin(TL[3]/(TL[2]-TL[1]))-asin(TL[3]/(TL[1]+TL[2]));/*最小传动角校核*/cd1=acos((TL[1]-TL[3])/TL[2])*180/PI;cd2=acos((TL[1]+TL[3])/TL[2])*180/PI;K=(PI+THITA)/(PI-THITA);printf("行程速比系数K=%6.4f\n\n",K);fprintf(fp,"行程速比系数K=%7.4f\n\n",K);printf("曲柄滑块机构工作行程最小传动角：γ1=%7.4f\n\n",cd1);fprintf(fp,"曲柄滑块机构工作行程最小传动角：γ1=%10.4f\n\n",cd1);printf("曲柄滑块机构回程最小传动角：γ2=%7.4f\n\n",cd2);fprintf(fp,"曲柄滑块机构回程最小传动角：γ2=%7.4f\n\n",cd2);/*校核后输出最小传动角*/constintN=11;/*宏定义*/ /*由滑块位移求对应的曲柄转角*/ floatBATA=14.84722*PI/180,S[N],PSI[N]; S[1]=130; S[10]=370; for(inti=2;i<N;i++)S[i]=S[1]+(S[10]-S[1])*(i-1)/9; CPSI(N,BATA,TL,S,PSI,fp);floatGAMA10=30.79111*PI/180,FI0=111.42795*PI/180,FI[N];/*由滑块的行程，通过传动机构II求得对应曲柄主动件的转角*/ floatH[]={0,145.0,178.0,175.0,65.5}; CFI(N,GAMA10,FI0,H,S,FI,fp); floatDL=55.0,TL1[5];/*通过函数平方逼近法设计传动机构I中的四杆机构*/ SQU(N,FI,PSI,DL,TL1,fp); floatFAI=0,VF,VCIR,SP,SCIR,fai[37]; intm=0; /*通过版台与滚筒的速度关系求得FAI角与两者的位移、速度和加速度的对应关系*/ floatTLL[]={0,0,127.6568f*PI/180,14.84722f*PI/180}; while(FAI<=360) { VEL1(FAI*PI/180,TL1,H,TL,TLL,&VF,&VCIR,&SCIR,&SP); WRT1(FAI,SP,SCIR,VF,VCIR,fp); FAI+=10;/*********************fai[]是补偿位移***********************/fai[m++]=SP-SCIR; printf("fai[%d]=%f\n",(m-1),fai[m-1]); fprintf(fp,"%f\n",fai[m-1]); }}voidVEL1(floatFAI,float*TL1,float*TL2,float*A,float*TLL,float*VF,float*VCIR,float*SCIR,float*SP){ floatFA[3]; FA[1]=FAI; POS3(A,FA); floatSET[4]; SET[1]=3.1415926f/2+FAI+TLL[3]; POS2(TL1,SET,1,2); floatALP=SET[3]-3.1415926f; floatSSET[4]; VEL(ALP,VF,VCIR,TL1,TL2,A,TLL,SET,SSET,FA); floatRCIR=180.0,BSSET1;staticconstfloatSSETO=30.79111f*3.1415926f/180;BSSET1=SSET[1]-SSETO;*SCIR=fabs(BSSET1)*RCIR;PS(A,FAI); }voidWRT1(floatFAI,floatSP,floatSCIR,floatVP,floatVCIR,FILE*fp){ printf("FAI=%10.4f°,SP=%10.4f,SCIR=%10.4f,VP=%10.4f,VCIR=%10.4f\n",FAI,SP,SCIR,VP,VCIR); fprintf(fp,"FAI=%10.4f°,SP=%10.4f,SCIR=%10.4f,VP=%10.4f,VCIR=%10.4f\n",FAI,SP,SCIR,VP,VCIR); return;}/*********************下面是版台位移子程序***********************/voidPS(float*A,float*FAI){*SP=2*(sqrt(pow(A[1]+A[2],2)-pow(A[3],2))-A[1]*cos(FA[1])-A[2]*cos(FA[2]));}(2)程序清单二（全局变量应用——参考样板——2010.7.8）注：程序一可以完成全部工作，程序二尝试使用新方法。但是由于结果未能全部得出，因时间关系未能修改，后头附录部分不含此程序#include<stdio.h>#include<math.h>#defineN11#definePI3.1415926FILE*fp1,fp2;inti,j,k,IFA,A[4][4],B[4],I1,I2,I3,I4,I;doubleALAMT=3.55;ER=0.305;HD0=397.5;doubleHD,K,gamag,gamah,GAMA10,FI0,DL,ALP,DALP=7.85,RCIR=180.0;doubleTL[5],M[5],TLL[4],TL1[5],L1[5],TL2[4],S[N],FI[N],FAI[3],SET[4],SSET[4];doubleH[5]={0,145.0,178.0,175.0,65.5};doubleBATA,PSI[N],*SE,Z,E,F,G;double*VF,DSET[4],*DDSET1,SP,AP,X,*SCIR,ACIRN,ACIRT,ACIR,*VCIR,SSETR;doublem,n,p;/*以上是预处理部分和全局变量定义部分*/voidmain()/*主函数部分*/{ fp1=fopen("f:\\cc.txt","a+");fp2=fopen("f:\\ccq.txt","a+");/*以追加方式打开，创建新文件*/ CRLE();/*求取K和最小传动角*/K=(PI+sin((ER*ALAMT)/(ALAMT-1))-sin((ER*ALAMT)/(ALAMT+1)))/(PI-sin((ER*ALAMT)/(ALAMT-1))+sin((ER*ALAMT)/(ALAMT+1)));gamag=acos(1/ALAMT-ER);gamah=acos(1/ALAMT+ER); fprintf(fp1,"K=%10.4f,gamag=%10.4f,gamah=%10.4f\n",K,gamag,gamah);printf("K=%10.4f,gamag=%10.4f,gamah=%10.4f\n",K,gamag,gamah);/*调用CFI函数*/FI0=111.42795/180.0*PI;GAMA10=30.79111/180.0*PI;CFI();for(j=1;j<N;j++){S[1]=130;S[N-1]=370;S[j]=S[1]+(S[N-1]-S[1])*(j-1)/9;printf("S[j]=%10.4f\n",S[j]);/*S数组从S[0]到S[10]*/fprintf(fp1,"S[j]=%10.4f\n",S[j]);}/*CPdI确定φ1i*/BATA=14.84722*PI/180;CPSI();/*调用SQU确定双曲柄机构*/DL=TL1[3]=55;SQU();for(k=0;k<=4;k++){printf("TL1[k]=%10.4f\n",i,TL1[k]);fprintf(fp1,"TL1[%d]=%10.4f\n",k,TL1[k]);}/*调用POS3确定曲柄滑块位置角FAI[3]*/DALP=7.85; for(IFA=0;IFA<=360;IFA=IFA+5) {DALP=7.85;SET[1]=(IFA+90.0)/180.0*PI+BATA;/*P17*/if(SET[1]>=(360*3.1415926/180)) SET[1]=SET[1]-360*PI/180;POS2();/*求出theta1.2.3*/for(i=1;i<=N-1;i++){printf("曲柄滑块位置角FAI[3]=%10.4f\n",FAI[3]);fprintf(fp1,"曲柄滑块位置角FAI[3]=%10.4f\n",FAI[3]);}ALP=SET[3]-PI;/*alpha*/FAI[1]=IFA/180.0*PI;POS3();SSET[3]=PI+TLL[2]-ALP;/*求伽马3*/for(i=1;i<5;i++)M[i]=TL[i];VEL();DSET[3]=DALP;ACCE();X=TL[1]*cos(FAI[1])+TL[2]*cos(FAI[2]);SP=2*(sqrt(pow((TL[1]+TL[2]),2)-pow(TL[3],2))-X);AP=TL[1]*pow(DSET[1],2)*(cos(FAI[1])+1/ALAMT*cos(2*FAI[1])+ER/ALAMT*sin(FAI[1]))+TL[1]*(*DDSET1)*(sin(FAI[1])+0.5*1/ALAMT*sin(2*FAI[1])-ER/ALAMT*cos(FAI[1]));SCI();ACCE();ACIRT=RCIR*(*DDSET1);ACIRN=RCIR*pow(DSET[1],2);ACIR=sqrt(pow(ACIRT,2)+pow(ACIRN,2));WRT();}fcloseall();}/*相对于前者，程序二具有更大的行程速比系数（1.18多），工作行程传动角高达2.2弧度多，传力性能好,这两者是优势；但是回程却不足1弧度，这是相对于前者的劣势。*/（3）VC屏幕截图2、结果记录及处理（1）写入CC.txt的数据HD0=397.5000,R0=197.9511,L0=702.7263,E0=60.3751HD=397.6061,R=198.0000,L=702.5000,E=60.50002HD=795.2122行程速比系数K=1.0343曲柄滑块机构工作行程最小传动角：γ1=78.7127曲柄滑块机构回程最小传动角：γ2=68.4094PSI=(DEGREE)133.3740140.7340147.9701155.2380162.6878170.4900178.8701188.1747199.0442213.0665ALFA0'=2.2280316(RADIAN),127.6568115(DEGREE)FI=(DEGREE)176.5468186.5587196.7699207.4588218.8301231.0278244.1497258.2710273.4776289.9133{C}{R}8.073707-4.371522-8.9307034.730534-4.3715222.7935495.000989-3.495667-8.9307035.00098910.000000-5.482457P0=1.164670,P1=-2.953947,P2=1.969152N=1.164670,U=0.394276,M=0.979352,A/A=1.000000L10=162.4671,L20=136.6159,L30=139.4962,L40=55.0000L1=162.5000,L2=136.5000,L3=139.5000,L4=55.0000FAI=0.0000°,SP=1.1507,SCIR=990.2909,VP=252.1273,VCIR=781.6365-989.140137FAI=10.0000°,SP=2.9182,SCIR=971.8409,VP=-415.8387,VCIR=808.7986-968.922729FAI=20.0000°,SP=19.8866,SCIR=953.3945,VP=-1112.0031,VCIR=846.2148-933.507874FAI=30.0000°,SP=51.0962,SCIR=934.8199,VP=-1807.5989,VCIR=891.2706-883.723755FAI=40.0000°,SP=94.9381,SCIR=916.0762,VP=-2475.0527,VCIR=941.4937-821.138123FAI=50.0000°,SP=149.2708,SCIR=897.1959,VP=-3088.2798,VCIR=994.5869-747.925049FAI=60.0000°,SP=211.5717,SCIR=878.2615,VP=-3622.3022,VCIR=1048.5588-666.689819FAI=70.0000°,SP=279.1143,SCIR=859.3773,VP=-4052.6907,VCIR=1101.8580-580.263062FAI=80.0000°,SP=349.1545,SCIR=840.6374,VP=-4355.6792,VCIR=1153.4521-491.482849FAI=90.0000°,SP=419.1063,SCIR=822.0983,VP=-4509.8535,VCIR=1202.8230-402.992004FAI=100.0000°,SP=486.6839,SCIR=803.7589,VP=-4499.6470,VCIR=1249.8972-317.074982FAI=110.0000°,SP=549.9942,SCIR=785.5508,VP=-4319.7173,VCIR=1294.9478-235.556641FAI=120.0000°,SP=607.5717,SCIR=767.3376,VP=-3978.3701,VCIR=1338.4949-159.765991FAI=130.0000°,SP=658.3586,SCIR=748.9167,VP=-3498.3093,VCIR=1381.2395-90.558167FAI=140.0000°,SP=701.6453,SCIR=730.0211,VP=-2914.1807,VCIR=1424.0325-28.375732FAI=150.0000°,SP=736.9883,SCIR=710.3148,VP=-2267.7480,VCIR=1467.884426.673523FAI=160.0000°,SP=764.1232,SCIR=689.3835,VP=-1602.2343,VCIR=1514.019774.739685FAI=170.0000°,SP=782.8859,SCIR=666.7106,VP=-957.2432,VCIR=1563.9883116.175354FAI=180.0000°,SP=793.1508,SCIR=641.6490,VP=-365.0146,VCIR=1619.8324151.501709FAI=190.0000°,SP=794.7888,SCIR=613.3765,VP=151.8728,VCIR=1684.3499181.412292FAI=200.0000°,SP=787.6499,SCIR=580.8452,VP=581.6070,VCIR=1761.4268206.804688FAI=210.0000°,SP=771.5698,SCIR=542.7314,VP=923.3455,VCIR=1856.3159228.838379FAI=220.0000°,SP=746.4041,SCIR=497.4214,VP=1186.5912,VCIR=1975.1353248.982758FAI=230.0000°,SP=712.0909,SCIR=443.1316,VP=1389.4678,VCIR=2121.0999268.959259FAI=240.0000°,SP=668.7411,SCIR=378.4326,VP=1555.3607,VCIR=2280.7900290.308594FAI=250.0000°,SP=616.7496,SCIR=303.7158,VP=1707.5105,VCIR=2393.7710313.033783FAI=260.0000°,SP=556.9097,SCIR=223.7006,VP=1862.0537,VCIR=2341.9021333.209106FAI=270.0000°,SP=490.5096,SCIR=147.8094,VP=2021.5856,VCIR=2064.3462342.700226FAI=280.0000°,SP=419.3802,SCIR=84.2898,VP=2172.3689,VCIR=1679.5896335.090332FAI=290.0000°,SP=345.8695,SCIR=34.7991,VP=2287.2510,VCIR=1341.2136311.070374FAI=300.0000°,SP=272.7411,SCIR=3.5051,VP=2333.2852,VCIR=1100.5013269.235992FAI=310.0000°,SP=203.0031,SCIR=34.1266