有限元-齿轮.docxVIP

有限元作业一直齿齿轮，齿数为20，模数为3，齿厚20mm，压力角为20。，试分析其齿的受力状况。Ansys命令流：finish/clear,all/prep7m=3 !模数z=20 !齿数pi=3.14 a=20*pi/180 !压力角r=0.5*m*z !分度圆半径rb=r*cos(20*pi/180) !基圆半径ha=1*m !齿顶高hf=(1+0.25)*m !齿根高ra=r+ha !齿顶圆半径rf=r-hf !齿根圆半径ad=acos(rb/ra) !齿顶圆压力角csys,1 $n=20 !设置柱坐标系，并设20个点*do,i,1,nai=ad/(n-1)*(i-1) !求得i点的压力角ri=rb/cos(ai) !求得i点的rcita=tan(ai)-ai !求得i点的θk,i,ri,cita*180/pi*enddo*do,i,1,n-1l,i,i+1*enddocita=(tan(20*pi/180)-(20*pi/180))*180/picita=cita+360/(4*z)lgen,,all,,,,-cita,,,,1csys,0 !设置直角坐标系,准备对称生成线k,n+1,kx(1)-(rb-rf),ky(1) !创建关键点，采用渐开线始点的切线l,1,n+1lcomb,allnumcmp,all !合并所有线，并压缩图素编号lsymm,y,all !关于x-z平面对称csys,1 !设置柱坐标系k,5,rf,360/(2*z) !创建关键点5k,6,rf,-360/(2*z) !创建关键点6kp1=knear(6) !得到距离关键点 6 最近的点l,6,kp1 !创建线kp1=knear(5) !得到距离关键点 5 最近的点l,5,kp1 !创建线lfillt,1,3,0.38*mlfillt,2,4,0.38*m !齿根角倒角ksel,s,loc,x,ra !选择齿顶的两个关键点*get,kp1,kp,0,num,min !得到选择集中的最小关键点号kp2=kpnext(kp1)l,kp1,kp2 !创建线allselnumcmp,alllgen,z,all,,,,360/z !复制单齿廓线Z次，形成整个齿廓线nummrg,kp !合并关键点al,all !创建齿廓面cyl4,,,r/5 !创建轴孔圆面asba,1,2 !齿廓面减轴孔面numcmp,allet,1,plane42 !定义单元类型及材料特性mp,ex,1,2.1e5mp,prxy,1,0.3mp,dens,1,7.85e-6 amesh,allcsys,1nsel,s,loc,x,r/5nrotat,alld,all,allallselcsys,-1f,734,fx,-1500*9.8*cos(20*pi/180)f,734,fy,-1500*9.8*sin(20*pi/180)/solusolve单元应力云图：X方向应力云图：Y方向应力云图：