斜齿轮参数化建模

本文介绍了ANSYS APDL斜齿轮参数化建模一则实例。 finish

/clear

PI=ACOS(-1)

mn=3

an=20

bita=14

aa=180

han=1

cn=0.25

T=286470

KAVH=2.36

z1=19

z2=96

xn1=0.5067

xn2=0.2646

b1=50

b2=50

kr1=10

kr2=100

mt=mn/COS(bita*PI/180)

at=ATAN(TAN(an*PI/180)/COS(bita*PI/180))

hat=han*COS(bita*PI/180)

ct=cn*COS(bita*PI/180)

xt1=xn1*COS(bita*PI/180)

xt2=xn2*COS(bita*PI/180)

a=mt/2*(z1+z2)

*IF,xn1+xn2,EQ,0,THEN

aa=a

*ELSE

aa=aa

*ENDIF

att=ACOS(a/aa*COS(at))

cigmat=xt1+xt2-(aa-a)/mt

r1=mt*z1/2

r2=mt*z2/2

rb1=r1*COS(at)

rb2=r2*COS(at)

ra1=r1+(hat+xt1-cigmat)*mt

ra2=r2+(hat+xt2-cigmat)*mt

rf1=r1-(hat+ct-xt1)*mt

rf2=r2-(hat+ct-xt2)*mt

rj1=r1*COS(at)/COS(att)

rj2=r2*COS(at)/COS(att)

s1=mn*(PI/2+2*xn1*TAN(an*PI/180))

s2=mn*(PI/2+2*xn2*TAN(an*PI/180))

inva=TAN(an*PI/180)-an*PI/180

aj1=ACOS(rb1/rj1)

aj2=ACOS(rb2/rj2)

invaj1=TAN(aj1)-aj1

invaj2=TAN(aj2)-aj2

sb1=(s1/r1+2*inva)*180/PI

sb2=(s2/r2+2*inva)*180/PI

sj1=(s1/r1-2*(invaj1-inva))*180/PI

sj2=(s2/r2-2*(invaj2-inva))*180/PI

aat1=ACOS(rb1/ra1)

aat2=ACOS(rb2/ra2)

ea=(z1*(TAN(aat1)-TAN(att))+z2*(TAN(aat2)-TAN(att)))/2/PI b=MAX(b1,b2)

eb=b*SIN(bita*PI/180)/PI/mn

!主动齿轮建模程序

/PREP7

ET,1,185

ESIZE,r1/10

MP,EX,1,ex_mat1

MP,PRXY,1,prxy_mat1

MP,MU,1,mu_mat

CSYS,1

K,1,0,0,0

K,2,kr1,90,0

CSYS,0

*DO,i,0.03,SQRT((ra1**2-rb1**2)/rb1**2),0.01

K,i*100,rb1*SIN(i)-rb1*i*COS(i),rb1*COS(i)+rb1*i*SIN(i) !BSPLINE,i*100,i*100-1

k1=i*100

*ENDDO

ksel,s,,,2,k1,1

bspline,all

CSYS,1

ss=k1+1

K,ss,kr1,90-sb1,0

CSYS,0

WPROT,-sb1,0,0

CSYS,4

*DO,i,0.01,SQRT((ra1**2-rb1**2)/rb1**2),0.01

K,i*100+ss,-rb1*SIN(i)+rb1*i*COS(i),rb1*COS(i)+rb1*i*SIN(i) k2=i*100+ss

*ENDDO

ksel,s,,,ss,k2,1

bspline,all

ALLSEL

LARC,k1,k2,1,ra1

LARC,2,ss,1,kr1

AL,ALL

CSYS,1

K,k2+1,rf1,90-sb1/2+180/z1,0

K,k2+2,rf1,90-sb1/2-180/z1,0

LARC,k2+1,k2+2,1,rf1

K,k2+3,kr1,90-sb1/2+180/z1,0

K,k2+4,kr1,90-sb1/2-180/z1,0

LARC,k2+3,k2+4,1,kr1

LSTR,k2+1,k2+3

LSTR,k2+2,k2+4

AL,5,6,7,8

*IF,90-sb1/2+180/z1,LT,90,THEN

APTN,1,2

ADELE,3,4,1,1

AADD,ALL

LFILLT,16,21,0.38*mn

LFILLT,19,22,0.38*mn

ADELE,ALL,,,0

LDELE,2,,,0

LDELE,4,,,0

LDELE,7,8,1,0

LCOMB,11,15

LCOMB,14,18

AL,ALL

LESIZE,21,(ra1-rb1)/10

LESIZE,22,(ra1-rb1)/10

*ELSE

AADD,1,2

LFILLT,10,13,0.38*mn

LFILLT,12,14,0.38*mn

ADELE,ALL,,,0

LDELE,2,,,0

LDELE,5,,,0

LDELE,15,16,1,0

LSEL,S,,,4

LSEL,A,,,9,11,2

LCOMB,ALL

ALLSEL

AL,ALL

LESIZE,10,(ra1-rb1)/10

LESIZE,12,(ra1-rb1)/10

*ENDIF

K,k2+5,rb1,90-sb1/2,0

K,k2+6,rb1,90-sb1/2-b1/rb1*TAN(bita*PI/180)*COS(at)/PI*180,-b1 L,k2+5,k2+6

agen,z1,1,,,0,360/z1,0

aadd,all

vdrag,z1+1,,,,,,2

vsweep,all

asel,s,,,z1*7+2,z1*7+1+z1,1

NSLA,S,1

*GET,nN,NODE,,COUNT

NROTAT,ALL

D,ALL,,0,,,,UX,UZ

F,ALL,FY,T*KAVH/kr1/nN

ALLSEL

SAVE

FINISH