用ansys绘制螺旋锥齿轮.doc

用ansys绘制螺旋锥齿轮请教： 用ansys画螺旋锥齿轮，其中齿向线都是方程表示的，请教各位高手，可不可行？chenhuxing (陈胡兴)当然可行，我刚建完一个圆弧齿轮。 直接在ansys中建模的 你可以参照画螺栓的程序我建了一个标准直齿轮的自动生成程序， 当然可行，我刚建完一个圆弧齿轮。直接在ansys中建模的 我想对比一下，能给发一份吗？ xiaopj2002_xiaopj!仅适用于标准直齿圆齿轮 !一般情况下仅需更改 m 、z 两个参数 !特殊情况下angle1、ha、c 等参数可以更改 !对齿根过渡圆角的处理分两种情况：整圆弧过渡及0.38m倒角 !因为 circle 命令不受坐标系影响，所以必须有初始偏转。 !Pj2004.2 /title,gear m=3.0 !定义模数 z=21.0!齿数 !齿轮相关数据计算 pi=3.14159265358979 angle1=20*pi/180.0 !压力角 ha=1.0 !正常齿制的齿顶高系数 *if,m,LT,1,then !正常齿制的顶隙系数选择 c=0.35 *else c=0.25 *endif r=0.5*m*z !分度圆半径 rb=r*cos(angle1) !基圆半径 ra=0.5*(z+2*ha)*m !齿顶圆半径 rf=0.5*(z-2*ha-2*c)*m!齿根圆半径zz=2*(ha+c)/(1-cos(angle1)!以次判断用何种方法(基于对rf、rb大小的判断) *if,z,GE,zz,then rff=0.38*m !齿根圆角半径 *else rff=(2.0*pi*rb/z-m*cos(angle1)*(pi/2.0+z*(tan(angle1)-angle1)/2.0 !同上 *endif angle2=180.0/z !镜像旋转角 angle3=(pi/2.0/z+tan(angle1)-angle1)*180.0/pi !初始的偏转角(其值不能改变）angle4=(pi/2.0/z+tan(angle1)-angle1)*180.0/pi !对称偏转角/prep7 csys,4 !激活工作坐标系 wprot,-angle3,0,0 !初始偏转角度 K,1,0,0 !绘制原点（编号为1） *do,t,0,1,0.01 !描点 *SET,x,rb*(cos(t)+t*sin(t)!渐开线方程 *SET,y,rb*(sin(t)-t*cos(t) k,x,y,0 !开始描点*enddo flst,3,101,3 !连点成线 *do,t,2,102 fitem,3,t *enddo bsplin,p51x !用多意线连点成线 KDELE,3,101,1 !删多余点 NUMCMP,KP !重排点号 wprot,angle4,0,0 !X轴旋转angle4度 对称偏转角 LSYMM,Y,1, , , ,0,0 !镜像渐开线 circle,1,ra !画出齿顶圆线 LSBL,1,6 !对渐开线进行分割 LSBL,2,3 ldele,4 !删除多余的线与点 ldele,5 ldele,6 ldele,8 kdele,6,9,1 kdele,3 kdele,5 NUMCMP,line !重排线号、点号 NUMCMP,KP *if,z,GE,zz,then !判断是否进行第二次分割 circle,1,rf LSBL,1,3 !对渐开线进行分割 LSBL,2,6ldele,4 !删除多余的线与点 ldele,5 ldele,1 ldele,7 kdele,7,9,1 kdele,2 kdele,3 NUMCMP,line NUMCMP,KP larc,2,3,1,ra !联接齿顶 *else larc,4,5,1,ra !联接齿顶 FLST,3,3,4,ORDE,2 FITEM,3,1 FITEM,3,-3 LGEN,2,P51X, , , , , , ,0 !复制，重排线号 ldele,1,3,1,1 NUMCMP,line NUMCMP,KP *endif *do,t,1,z-1,1 !开始旋转镜像以上所形成的齿廓曲线 wprot,angle2,0,0 FLST,3,3,4,ORDE,2 FITEM,3,1 FITEM,3,-3 LSYMM,Y,P51X, , , ,0,0 *enddo WPCSYS,-1,0 !重新对齐工作坐标并开始做齿根圆圆弧 *if,z,GE,zz,then !方法选择 larc,4,9,1,rf !当rb<rf时，直接连接rf，并倒0.38m圆角 *do,t,7,4*z-5,4 larc,t,t+6,1,rf*enddo larc,4*z-1,5,1,rf lfillt,1,4*z,rff *do,t,0,z-2,1 lfillt,5+3*t,3*z+1+t,rff lfillt,4+3*t,3*z+2+t,rff *enddo lfillt,2,3*z+1,rff *else csys,2 !当rb>rf时，用计算圆角倒齿根圆角 *do,t,1,2*z,2 k,rf*1.1,angle2*t,0 *enddo csys,0 larc,2,6,4*z+2,rff *do,t,0,z-3,1 larc,8+4*t,10+4*t,4*z+3+t,rff *enddo larc,4,4*z,5*z+1,rff *endif finish save !okpj 2004.2 !对齿根的过渡还有点问题 !有些地方还相当累赘，帮忙减减肥。 !若有人人乐意帮忙把变位给加上，先谢！ !要能任意旋转更好xyni看看我写的，可以随意改变M,Z /prep7 m=4 !齿轮模数 z1=45 !齿轮齿数 pi=3.1415926 ang=20 !分度圆上的压力角 ha_c=1 ! 齿顶高系数 c_c=0.25 !顶隙系数 ha=ha_c*m !齿顶高 hf=(ha_c+c_c)*m !齿根高 d1=m*z1 !分度圆直径 db=d1*cos(ang*pi/180) !基圆直径 da=d1+2*ha !齿顶圆直径 df=d1-2*hf !齿根圆直径 alfa1=acos(db/da)*180/pi !齿顶圆压力角（角度表示） alfa2=acos(db/df)*180/pi !齿根圆压力角（角度表示） dfr=0.38*m !齿根圆角半径 b=12 !齿宽 *dim,alfa,array,10 *dim,rk,array,10 *dim,theta,array,10 *do,i,1,10 alfa(i)=(alfa2+(alfa1-alfa2)/9)*(i-1)*pi/180 !每条渐开线上生成10各关键点所对应的压力角大小（弧表示） *enddo *do,i,1,10 rk(i)=(db/2)/cos(alfa(i) theta(i)=tan(alfa(i)-alfa(i) !弧度表示 *enddo csys,1 *do,i,1,10 k,i,rk(i),theta(i)*180/pi *enddo *do,i,1,9 l,i,i+1 *enddo alls lcomb,all numcmp,all csys,1 n,1,m*z1,(tan(ang*pi/180)-ang*pi/180)*180/pi,0 csys,0 n,2,0,0,0 n,3,0,0,1 cs,11,0,1,2,3 csys,11 lsymm,z,all,0,0 csys,1 lgen,2,360/(z1*2),1 alls lgen,z1,all,360/z1, lplot csys,1 *do,i,2,z1*4-2,4 l,i,i+2 *enddo *do,i,3,z1*4-3-2,4 l,i,i+2 *enddo l,1,z1*4-1 *do,i,2,z1*2,2 lfillt,i,z1*3+i/2,dfr *enddo *do,i,3,z1*2-1,2 lfillt,i,z1*3+(i-1)/2,dfr *enddo lfillt,1,z1*4,dfr alls al,all, alls voffst,1,b/2 *get,line_count,line,count *dim,line_length,array,line_count *do,i,1,line_count *get,line_length(i),line,i,leng *enddo line_min=1000 *do,i,1,line_count *if,line_min,gt,line_length(i),then line_min=line_length(i) *endif *enddo *do,i,1,line_count lesize,i,line_min *enddoxiaopj xyni我粗略的看了看你的程序，很好，比我的苗条多了。而且齿根过渡上的处理也比我好。但是仍