ansys齿轮接触分析命令流

/prep7
!单位是n,和mm
!先建立两个齿轮模型（没又齿根过渡线，只适合基圆半径小于齿根圆的情况）

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 !齿根高
d=m*z1 !分度圆直径
db=d*cos(ang*pi/180) !基圆直径
da=d+2*ha !齿顶圆直径
df=d-2*hf !齿根圆直径
x=0
s=pi*m/2+2*x*m*tan(ang*pi/180) !分度圆齿厚
theta_s=tan(ang*pi/180)-ang*pi/180
fai_s=s/d
alfa_a=acos(db/da)*180/pi !齿顶圆压力角（角度表示）
alfa_f=acos(db/df)*180/pi !齿根圆压力角（角度表示）
dfr=0.38*m !齿根圆角半径
!b=0.012 !齿宽
deata_ang=360/z1 !齿轮两齿间的夹角

*dim,alfa,array,10
*dim,rk,array,10
*dim,theta,array,10
*dim,fai,array,10

*do,i,1,10
alfa(i)=(alfa_f+((alfa_a-alfa_f)/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

*do,i,1,10
fai(i)=(theta_s+fai_s-theta(i))*180/pi !角度
*enddo

CSYS,1

*do,i,1,10
k,i,rk(i),fai(i)
*enddo

ksel,all
BSPLIN,ALL !绘制齿廓线

csys,0

LSYMM,Y,1
K,30000

LARC,10,12,30000,da/2
csys,1
LSEL,ALL
LGEN,Z1,ALL,,,0,-deata_ang


kpoint=10+1
*do,i,1,z1-1,1
larc,kpoint,kpoint+2,30000,df/2
kpoint=kpoint+4
*enddo

larc,187,1,30000,df/2
LSEL,ALL
AL,ALL
CYL4, , ,50
ASBA, 1, 2
csys,0
AGEN,2,3, , ,1.00*d, , , ,0
!*********************************************************
ET,1,PLANE82
MP,EX,1,500000
MP,NUXY,1,0.3

MP,EX,2,400000
MP,NUXY,2,0.3

TYPE, 1
MAT, 1
REAL,
ESYS, 0
SECNUM,
!*
SMRT,3
MSHAPE,0,2D
MSHKEY,0
AMESH,3
!*

MAT, 2
AMESH,1

FLST,5,90,2,ORDE,73
FITEM,5,5
FITEM,5,-6
FITEM,5,61
FITEM,5,-62
FITEM,5,83
FITEM,5,-84
FITEM,5,89
FITEM,5,-90
FITEM,5,145
FITEM,5,-147
FITEM,5,180
FITEM,5,-181
FITEM,5,187
FITEM,5,211
FITEM,5,-212
FITEM,5,216
FITEM,5,-217
FITEM,5,242
FITEM,5,-243
FITEM,5,248
FITEM,5,-249
FITEM,5,393
FITEM,5,-394
FITEM,5,444
FITEM,5,-446
FITEM,5,461
FITEM,5,-464
FITEM,5,483
FITEM,5,-484
FITEM,5,525
FITEM,5,584
FITEM,5,587
FITEM,5,624
FITEM,5,635
FITEM,5,641
FITEM,5,643
FITEM,5,681
FITEM,5,703
FITEM,5,-704
FITEM,5,732
FITEM,5,737
FITEM,5,774
FITEM,5,-775
FITEM,5,864
FITEM,5,-865
FITEM,5,882
FITEM,5,-883
FITEM,5,895
FITEM,5,903
FITEM,5,-904
FITEM,5,922
FITEM,5,941
FITEM,5,961
FITEM,5,967
FITEM,5,-968
FITEM,5,982
FITEM,5,1000
FITEM,5,-1003
FITEM,5,1114
FITEM,5,-1117
FITEM,5,1127
FITEM,5,1163
FITEM,5,-1166
FITEM,5,1171
FITEM,5,-1176
FITEM,5,1266
FITEM,5,1270
FITEM,5,1310
FITEM,5,-1312
FITEM,5,1343
FITEM,5,-1346
FITEM,5,1350
FITEM,5,-1351
CM,_Y,ELEM
ESEL,

, , ,P51X
CM,_Y1,ELEM
CMSEL,S,_Y
CMDELE,_Y
!*
!*
EREF,_Y1, , ,1,0,1,1
CMDELE,_Y1


CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.2
MAT,1
R,3
REAL,3
ET,2,169
ET,3,172
KEYOPT,3,9,0
KEYOPT,3,10,2
R,3,
RMORE,
RMORE,,0
RMORE,0
! Generate the target surface
LSEL,S,,,1
LSEL,A,,,2
LSEL,A,,,3
LSEL,A,,,4
LSEL,A,,,5
LSEL,A,,,6
LSEL,A,,,7
LSEL,A,,,8
LSEL,A,,,9
LSEL,A,,,10
LSEL,A,,,128
LSEL,A,,,130
LSEL,A,,,131
LSEL,A,,,132
LSEL,A,,,133
LSEL,A,,,134
LSEL,A,,,135
LSEL,A,,,136
LSEL,A,,,137
LSEL,A,,,138
LSEL,A,,,178
LSEL,A,,,179
LSEL,A,,,180
CM,_TARGET,LINE
TYPE,2
NSLL,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
LSEL,S,,,264
LSEL,A,,,265
LSEL,A,,,266
LSEL,A,,,267
LSEL,A,,,268
LSEL,A,,,269
LSEL,A,,,270
LSEL,A,,,271
LSEL,A,,,272
LSEL,A,,,273
LSEL,A,,,274
LSEL,A,,,275
LSEL,A,,,276
LSEL,A,,,277
LSEL,A,,,278
LSEL,A,,,279
LSEL,A,,,280
LSEL,A,,,281
LSEL,A,,,282
LSEL,A,,,283
LSEL,A,,,284
LSEL,A,,,285
LSEL,A,,,286
LSEL,A,,,287
LSEL,A,,,288
CM,_CONTACT,LINE
TYPE,3
NSLL,S,1
ESLN,S,0
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
NROPT,UNSYM
LOCAL,11,1,180,0,0, , , ,1,1,

EPLOT
CSYS,1
LSEL,ALL
FLST,5,4,4,ORDE,2
FITEM,5,181
FITEM,5,-184
LSEL,S, , ,P51X
NSEL,ALL
NSLL,S,1
NPLOT
FLST,2,88,1,ORDE,4
FITEM,2,181
FITEM,2,-184
FITEM,2,635
FITEM,2,-718
NROTAT,P51X
FINISH
/SOL
FLST,2,88,1,ORDE,4
FITEM,2,181
FITEM,2,-184
FITEM,2,635
FITEM,2,-718
!*
/GO
D,P51X, , , , , ,UX, , , , ,
FLST,2,88,1,ORDE,4
FITEM,2,181
FITEM,2,-184
FITEM,2,635
FITEM,2,-718
!*
/GO
D,P51X, ,0.6, , , ,UY, , , , ,

ALLSEL,ALL
CSYS,11,
EPLOT
LSEL,ALL
FLST,5,4,4,ORDE,2
FITEM,5,365
FITEM,5,-368
LSEL,S, , ,P51X
NSEL,ALL
NSLL,S,1
NPLOT
FINISH
/PREP7
FLST,2,88,1,ORDE,4
FITEM,2,2815
FITEM,2,-2818
FITEM,2,3269
FITEM,2,-3352
NROTAT,P51X
FINISH
/SOL
FLST,2,88,1,ORDE,4
FITEM,2,2815
FITEM,2,-2818
FITEM,2,3269
FITEM,2,-3352
!*
/GO
D,P51X, ,0, , , ,UX, , , , ,
ALLSEL,ALL
EPLOT
!*
!*
OUTRES,ALL,ALL,
!*
TIME,1
AUTOTS,-1
DELTIM,0.05,0.05,0.05,1
KBC,0
!*
TSRES,ERASE
solve
/POST1
!*
/EFACET,1
/DSCALE,ALL,1.0
PLNSOL, S,EQV, 0,1.0