准双曲面齿轮参数化计算程序
- 格式:pdf
- 大小:810.22 KB
- 文档页数:4
基于Excel、Matlab与UG的准双曲面齿轮精确建模研究准双曲面齿轮是一种特殊的齿轮,在机械传动领域中广泛使用。
准双曲面齿轮的特点是具有精确的啮合性能、传动能力强等优点。
因此,如何精确地建模准双曲面齿轮对于机械设计非常重要。
本文将介绍使用Excel、Matlab和UG软件进行准双曲面齿轮精确建模的研究。
首先,我们在Excel中建立了准双曲面齿轮的基本几何参数模型,包括齿轮的模数、齿数、压力角、齿宽等参数。
然后,通过Matlab编写程序,利用几何算法生成齿面曲线。
在程序中,我们使用了双曲线函数来描述准双曲面齿轮的齿面曲线。
通过调整参数,可以得到精确的齿面曲线。
接下来,我们将生成的齿轮齿面曲线导入到UG软件中进行建模。
首先,我们生成一个旋转体,即齿轮的基本体形。
然后,在基本体形上使用齿轮齿面曲线进行网格划分,从而得到具有精确齿形的齿轮模型。
通过使用Excel、Matlab和UG软件进行准双曲面齿轮精确建模,我们成功地解决了传统建模方法中难以精确建模的问题。
这种方法具有如下优点:1. 精度高。
利用双曲线函数生成齿面曲线,可以得到精确的齿形。
通过将齿面曲线导入到UG软件中,可以获得具有高精度的齿轮模型。
2. 灵活性好。
齿轮的基本几何参数模型可以根据需求灵活调整。
齿面曲线也可以通过调整参数进行优化,得到更加精确的齿形。
3. 使用方便。
本方案基于Excel、Matlab和UG软件,这些软件都是非常常用的工具。
因此,使用本方案进行准双曲面齿轮建模相对于其他方法更加容易上手。
综上所述,通过Excel、Matlab和UG软件进行准双曲面齿轮精确建模,可以得到具有高精度、灵活性好和使用方便等优点的齿轮模型。
这种方法可以被广泛应用于机械设计和制造领域。
数据分析是对已有数据进行收集、处理、分析、解释、推断以及展示等一系列过程的总称,是在统计学和数据挖掘等领域中常用的方法。
在现代社会,各行各业都需要利用数据分析来支持决策和管理。
准双曲面齿轮和螺旋锥齿轮设计的统一算法准双曲面齿轮传动是锥齿轮传动中的普遍形式,螺旋锥齿轮是它的一种特殊情况.当准双曲面齿轮的偏置距E12=0时,就成为螺旋锥齿轮传动.在外形和加工方法上,准双曲面齿轮与螺旋锥齿轮无本质区别,切齿计算方法差别也不大[1,2].在实际设计中,它们的几何计算方法却不相同.当偏置距E12趋近于零时,现行的准双曲面齿轮的几何计算公式误差增大,甚至失效.因此螺旋锥齿轮设计的几何计算不能采用准双曲面齿轮几何计算公式和计算方法.在CAD软件开发中必须对这两种锥齿轮分别进行处理. 作者提出一种适合于准双曲面齿轮和螺旋锥齿轮设计的统一几何计算方法,其特点是当偏置距E12较大时,它与准双曲面齿轮现行计算结果一致;当偏置距E12为零时,得到正确的螺旋锥齿轮几何参数;当E12较小时,计算误差很小.因此在锥齿轮CAD软件开发中,可将这两种锥齿轮甚至包括直齿锥齿轮统一处理.1 分度锥参数基本公式准双曲面齿轮与螺旋锥齿轮几何计算中最大的区别在于分度锥参数的确定方法.分析现行准双曲面齿轮几何计算公式可知,当偏置距E12趋近于零时,齿轮的偏置角η,ε,ε′也趋近于零,因而导致公式计算误差增大甚至失效.作者在分析过程中发现,虽然E12趋近于零时,齿轮的偏置角η,ε,ε′也趋近于零,但它们属于同阶无穷小.即极限和存在.令式中e1和e2为偏置角系数.根据偏置角系数,可给出分度锥参数基本公式为式中k为放大系数;上面这组基本公式不仅适合于准双曲面齿轮,也适合于螺旋锥齿轮,不会因E12=0而失效.2 分度锥参数的求解上面给出的基本公式是一组非线性方程组,其中有5个参数是在几何计算前确定的.根据传动和强度等要求先确定齿轮的偏置距E12,轴交角ζ=90°-Σ,齿轮齿数z1和z2,大齿轮中点端面模数mt2,小齿轮中点螺旋角β1.则上面基本公式中的已知参数为i12=z2/z1, r2=mt2z2/2,及E12,ζ,β1. 由于基本公式是非线性方程组,在此采用迭代法求解.即初选k和e1值,按下面步骤进行迭代:若|k*-k|≤ξ(由计算精度确定的某一小量),则可进行下面的迭代;否则改变k初值重新迭代.式中rc为刀盘半径. 若|k0-kc|>ξ,则改变e1初值重新迭代,直到|k0-kc|≤ξ为止.迭代完毕,便得到了所有的分度锥参数.然后根据齿宽、齿高系数、变位系数和齿根倾斜类型,按准双曲面齿轮的方法进行其它所有几何尺寸参数的计算.3 算例作者采用上面的统一公式和算法分别对准双曲面齿轮和螺旋锥齿轮两种情况进行了大量的计算分析.表1是偏置距E12=0的螺旋锥齿轮算例结果;表2是偏置距E12=30mm的准双曲面齿轮算例结果.大量的计算分析结果表明:当E12=0时,上面方法所确定的分度锥参数与现行螺旋锥齿轮几何计算结果一致;当E12≠0时,上面方法与现行准双曲面齿轮几何计算结果一致;特别是当E12非常小时,本方法所得结果比较精确.因此,可用上面方法将这两种锥齿轮的几何计算方法统一起来.这对CAD软件开发特别有利.表1 螺旋锥齿轮参数表2 准双曲面齿轮参数4 统一设计中的问题现行准双曲面齿轮和螺旋锥齿轮的标准参数,如模数、齿高系数、变位系数等都定义在大端.这对准双曲面齿轮会导致理论啮合节点偏离齿宽中点而与螺旋锥齿轮不同.因此建议将标准参数定义在齿宽中点,这样也可以与强度计算方法一致[3].此外,现行准双曲面齿轮标准参数中的螺旋角是小齿轮螺旋角,而标准参数中的模数是大齿轮端面模数.建议标准参数取大齿轮螺旋角和法向模数,这样更合理.5 结论大量算例和实际应用表明,作者提出的几何计算方法是可行的.作者已经根据此原理开发了CAD应用软件,并用于实际设计中.这样就使准双曲面齿轮和螺旋锥齿轮甚至直齿锥齿轮设计中的几何计算方法的统一有了依据.结果也在一定程度上揭示了准双曲面齿轮和螺旋锥齿轮理论上的本质联系.对锥齿轮的标准化、系列化和CAD技术也有一定的意义.。
/*圆弧齿双曲面齿轮的几何尺寸计算*/#include<math.h>#include<stdio.h>#define pi 3.141592654#define PF printfmain(){double s[150];double z1=7,z2=46,f2=55,e=38,d2=322,rd=114.3,beta1=44;doubles67l,s67r,s68l,s68r,s123l,s123r,s124l,s124r,s125l,s125r,s126l,s126r;int i;double O2,Q2,Qs,TR,QTR,O2T,Q2T;PF("\n圆弧齿双曲面齿轮的几何尺寸计算\n\n");s[1]=z1;s[2]=z2; PF("主动锥齿轮齿数 z1=7 从动锥齿轮齿 z2=46\n");s[3]=s[1]/s[2];s[4]=43.4; PF("大齿轮齿面宽 F2=55\n");s[5]=28; PF("小齿轮轴线偏移距 E=38\n");s[6]=280; PF("大齿轮分度圆直径 d2=322\n");s[7]=152.4; PF("刀盘名义半径 rd=114.3\n");s[8]=beta1*pi/180; PF("小齿轮螺旋角预选值 beta1=44\n");s[9]=tan(s[8]);s[10]=1.2*s[3];s[11]=sin(atan(1/s[10]));s[12]=(s[6]-s[4]*s[11])/2.0;s[13]=s[5]*s[11]/s[12];s[14]=cos(asin(s[13]));s[15]=s[14]+s[9]*s[13];s[16]=s[3]*s[12];s[17]=s[15]*s[16];s[18]=1.300;s[19]=s[12]/s[10]+s[17];s[20]=s[5]/s[19];s[21]=sqrt(1+s[20]*s[20]);s[22]=s[20]/s[21];s[23]=asin(s[22]);s[24]=(s[5]-s[17]*s[22])/s[12];s[25]=tan(asin(s[24]));s[27]=cos(atan(s[26]));s[28]=s[24]/s[27];s[29]=cos(asin(s[28]));s[30]=(s[15]-s[29])/s[28];s[31]=s[28]*(s[9]-s[30]);s[32]=s[3]*s[31];s[33]=s[24]-s[22]*s[32];s[34]=tan(asin(s[33]));s[35]=s[22]/s[34];s[36]=atan(s[35]); PF("小齿轮节锥角 r1=%f\n",s[36]*180/pi);s[37]=cos(s[36]);s[38]=s[33]/s[37];s[39]=asin(s[38]);s[40]=cos(s[39]);s[41]=(s[15]+s[31]-s[40])/s[38];s[42]=atan(s[41]); PF("小齿轮中点螺旋角 beta1=%f\n",s[42]*180/pi);s[43]=cos(s[42]);s[44]=s[42]-s[39]; PF("大齿轮中点螺旋角 beta2=%f\n",s[44]*180/pi);s[45]=cos(s[44]);s[46]=tan(s[44]);s[47]=s[22]/s[33];s[48]=atan(1/s[47]);PF("大齿轮节锥角 r2=%f\n",s[48]*180/pi);s[49]=sin(s[48]);s[50]=cos(s[48]);s[51]=(s[17]+s[12]*s[32])/s[37];s[52]=s[12]/s[50];s[53]=s[51]+s[52];s[54]=s[12]*s[45]/s[49];s[55]=s[43]*s[51]/s[35];s[56]=(s[41]*s[55]-s[46]*s[54])/s[53];s[57]=atan(s[56]);s[58]=cos(s[57]);s[59]=s[41]*s[56]/s[51];s[60]=s[46]*s[56]/s[52];s[61]=s[54]*s[55];s[62]=(s[54]-s[55])/s[61];s[63]=s[59]+s[60]+s[62];s[64]=(s[41]-s[46])/s[63];s[65]=s[64]/s[58];s[66]=s[7]/s[65];s67r=1-s[3];s68l=s[5]/s[34]-s[17]*s[35];s68r=s[35]*s[37];s[69]=s[37]+s[40]*s67l;s[70]=s[49]*s[51];s[71]=s[12]*s[47]-s[70]; PF("大齿轮节锥顶点到小齿轮轴线的距离Z=%f\n",s[71]);s[72]=s[12]/s[49];s[73]=0.5*s[6]/s[49]; PF("大锥齿轮节距 A0=%f\n",s[73]);s[74]=s[73]-s[72];s[75]=4*s[12]*s[45]/s[2];s[76]=s[12]*s[46]/s[7];s[77]=s[49]/s[45]-s[76];s[78]=pi*42.5/180;s[79]=sin(s[78]);s[80]=s[78]/2;s[81]=cos(s[80]);s[82]=tan(s[80]);s[83]=s[77]/s[82];s[84]=10560*s[83]/s[2];s[85]=0.170;s[86]=1.150-s[85];s[87]=s[75]*s[85];s[88]=s[75]*s[86]+0.05;O2=3438*s[87]/s[72];Q2=3438*s[88]/s[72];Qs=O2+Q2;TR=s[84]/Qs-s[18];if(TR>0){QTR=s[18]*Qs/60*pi/180; PF("TR=%f 为倾根锥母线收缩齿\n",TR);} else{QTR=s[84]/60*pi/180; PF("TR=%f 为双重收缩齿\n",TR);}O2T=s[85]*QTR;Q2T=QTR-O2T;s[89]=O2T; PF("大齿轮齿顶角 O2=%f\n",s[89]*180/pi*60);s[90]=sin(s[89]);s[91]=Q2T; PF("大齿轮齿根角 Q2=%f\n",s[91]*180/pi*60);s[92]=sin(s[91]);s[93]=s[87]+s[74]*s[90]; PF("大齿轮齿顶高 hd2=%f\n",s[93]);s[94]=s[88]+s[74]*s[92]; PF("大齿轮齿根高 hg2=%f\n",s[94]);s[95]=0.15*s[75]+0.05; PF("径向间隙 C=%f\n",s[95]);s[96]=s[93]+s[94]; PF("大齿轮齿全高 h=%f\n",s[96]);s[97]=s[96]-s[95]; PF("大齿轮齿工作高 hg=%f\n",s[97]);s[98]=s[48]+s[89]; PF("大齿轮面锥角 r02=%f\n",s[98]*180/pi);s[99]=sin(s[98]);s[100]=cos(s[98]);s[101]=s[48]-s[91]; PF("大齿轮根锥角 rR2=%f\n",s[101]*180/pi); s[102]=sin(s[101]);s[103]=cos(s[101]);s[104]=1/tan(s[101]);s[105]=s[93]*s[50]/0.5+s[6]; PF("大齿轮外圆直径 do2=%f\n",s[105]);s[106]=s[70]+s[74]*s[50];s[107]=s[106]-s[93]*s[49]; PF("大齿轮外缘至小齿轮轴线的距离Xo2=%f\n",s[107]);s[108]=(s[72]*s[90]-s[87])/s[99];s[109]=(s[72]*s[92]-s[88])/s[102];s[110]=s[71]-s[108]; PF("大齿轮面锥顶点到小齿轮轴线的距离Zo=%f\n",s[110]);s[111]=s[71]+s[109]; PF("大齿轮根锥顶点到小齿轮轴线的距离ZR=%f\n",s[111]);s[112]=s[12]+s[70]*s[104];s[113]=s[5]/s[112];s[114]=sqrt(1-s[113]*s[113]);s[115]=s[113]/s[114];s[116]=s[103]*s[114];s[117]=asin(s[116]); PF("小齿轮面锥角 ro1=%f\n",s[117]*180/pi);s[118]=cos(s[117]);s[119]=tan(s[117]);s[120]=(s[102]*s[111]+s[95])/s[103];s[121]=(s[5]*s[113]-s[120])/s[114]; PF("小齿轮面锥顶点到大齿轮轴线的距离 Go=%f\n",s[121]);s[122]=s[38]*s67l/s[69];s123l=atan(s[122]);s123r=cos(s123l);s124l=s[39]-s123l;s124r=cos(s124l);s125l=s[117]-s[36];s125r=cos(s125l);s126l=s[113]*s67r-s68r;s126r=-s[113]*s67r-s68r;s[127]=s123r/s124r;s[128]=s68l+s[87]*s68r;s[129]=s[118]/s125r;s[130]=s[74]*s[127];s[131]=s[128]+s[130]*s[129]+s[75]*s126l; PF("小齿轮外缘到大齿轮轴线的距离 BR=%f\n",s[131]);s[132]=s[4]*s[127]-s[130];s[133]=s[128]-s[132]*s[129]+s[75]*s126r; PF("小齿轮前缘至大齿轮轴线的距离 B1=%f\n",s[133]);s[134]=s[121]+s[131];s[135]=s[119]*s[134]/0.5; PF("小齿轮的外圆直径 do1=%f\n",s[135]);s[136]=s[70]*s[100]/s[99]+s[12];s[137]=s[5]/s[136];s[138]=asin(s[137]);s[139]=cos(s[138]);s[140]=(s[99]*s[110]+s[95])/s[100];s[141]=(s[5]*s[137]-s[140])/s[139]; PF("小齿轮根锥顶点到大齿轮轴线的距离 GR=%f\n",s[141]);s[142]=s[100]*s[139];s[143]=asin(s[142]); PF("小齿轮根锥角 rR1=%f\n",s[143]*180/pi);s[144]=cos(s[143]);s[145]=tan(s[143]);s[146]=0.239894; PF("最小齿侧间隙允许值 Bmin=0.239894\n");s[147]=0.324936; PF("最大齿侧间隙允许值 Bmax=0.324936\n");s[148]=s[90]+s[92];s[149]=s[96]-s[4]*s[148];s[150]=s[73]-s[4];PF("计算完毕!!!\n\n");/*用来检查每个数的数值*/for(i=1;i<151;i++){PF("\ns[%d]=%f\n",i,s[i]);}PF("%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f,%f\n",s67l,s67r,s68l,s68r,s123l,s12 3r,s124l,s124r,s125l,s125r,s126l,s126r);}。