下面这个例子是利用ANSYS进行齿轮接触分析的一段代码，本站没有验证，仅供参考。

FINISH
/CLEAR
/PREP7
/TITILE, WWW.MEKENICAL.COM
!单位是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, , ,D, , , ,0
!*********************************************************
!指定各条线划分的精度
NUMCMP,LINE
*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*2
*ENDDO

ET,1,PLANE42
MP,EX,1,30000
MP,NUXY,1,0.3
MSHAPE,0,2D
MSHKEY,0
AMESH,1   !划分网格
AMESH,3

!下面这些是手工操作局部细化网格时产生的命令，就是为了局部细化
!*****************************************************
FLST,5,12,1,ORDE,12
FITEM,5,134
FITEM,5,-135
FITEM,5,140
FITEM,5,-141
FITEM,5,146
FITEM,5,-147
FITEM,5,1016
FITEM,5,-1017 
FITEM,5,1022
FITEM,5,-1023 
FITEM,5,1280
FITEM,5,-1281 
CM,_Y,NODE
NSEL, , , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
CMDELE,_Y 
NREFINE,_Y1, , ,2,1,1,1
CMDELE,_Y1
!********************************************************

ET,2,CONTA175       ! 2-D CONTACT ELEMENTS
ET,3,TARGE169       ! 2-D TARGET ELEMENTS
MP,EX,1,30000       ! SMALLER CYLINDER PROPERTIES
MP,NUXY,1,0.25
MP,EX,2,29120       ! LARGER CYLINDER PROPERTIES
MP,NUXY,2,0.30

!ET,2,CONTAC48       !2D接触单元
!R,1,1E5, ,0.01, , , ,   !刚度1E5,TOLS=0.01
!RMORE, ,

!先选择6组节点
!FLST,5,4,1,ORDE,4   !这些也都是手工选取接触节点时产生的命令
!FITEM,5,1023 
!FITEM,5,1721
!FITEM,5,1731
!FITEM,5,1825
!NSEL,S, , ,P51X
!CM,TARGET1,NODE
!FLST,5,4,1,ORDE,4 
!FITEM,5,134
!FITEM,5,394
!FITEM,5,935
!FITEM,5,938
!NSEL,S, , ,P51X
!CM,CONTACT1,NODE
!
!FLST,5,6,1,ORDE,6 
!FITEM,5,1016
!FITEM,5,-1017 
!FITEM,5,1719
!FITEM,5,1757
!FITEM,5,1875
!FITEM,5,1882
!NSEL,S, , ,P51X
!CM,TARGET2,NODE
!FLST,5,5,1,ORDE,5 
!FITEM,5,141
!FITEM,5,390
!FITEM,5,393
!FITEM,5,925
!FITEM,5,934
!NSEL,S, , ,P51X
!CM,CONTACT2,NODE
!
!FLST,5,3,1,ORDE,3 
!FITEM,5,1280
!FITEM,5,1745
!FITEM,5,-1746 
!NSEL,S, , ,P51X
!CM,TARGET3,NODE
!FLST,5,4,1,ORDE,4 
!FITEM,5,147
!FITEM,5,395
!FITEM,5,409
!FITEM,5,915
!NSEL,S, , ,P51X
!CM,CONTACT3,NODE

LSEL,S,,,273 
NSLL,S,1   
CM,CONTACT1,NODE

LSEL,S,,,4 
NSLL,S,1   
CM,TARGET1,NODE

LSEL,S,,,275
NSLL,S,1
CM,CONTACT2,NODE
LSEL,S,,,2 
NSLL,S,1   
CM,TARGET2,NODE

LSEL,S,,,277       
NSLL,S,1           
CM,CONTACT3,NODE

LSEL,S,,,1 
NSLL,S,1   
CM,TARGET3,NODE


!三个对称的接触面

CMSEL,S,CONTACT1,NODE
TYPE,2
ESURF
CMSEL,S,TARGET1,NODE
TYPE,3
ESURF

CMSEL,S,CONTACT2,NODE
TYPE,2
ESURF

CMSEL,S,TARGET2,NODE
TYPE,3
ESURF

CMSEL,S,CONTACT3,NODE
TYPE,2
ESURF

CMSEL,S,TARGET3,NODE
TYPE,3
ESURF


!GCGEN,CONTACT1,TARGET1, , ,TOP,
!GCGEN,TARGET1,CONTACT1, , ,TOP,
!GCGEN,CONTACT2,TARGET2, , ,TOP,
!GCGEN,TARGET2,CONTACT2, , ,TOP,
!GCGEN,CONTACT3,TARGET3, , ,TOP,
!GCGEN,TARGET3,CONTACT3, , ,TOP,

!给两个内圆加约束，其中左边轮是主动的，故只约束其径向（允许有旋转），右边的全部约
!束
CSYS,1
FLST,5,4,4,ORDE,2 
FITEM,5,181
FITEM,5,-184
LSEL,S, , ,P51X
NSLL,S,1
FLST,2,52,1,ORDE,2
FITEM,2,1285
FITEM,2,-1336 
NROTAT,ALL
D,ALL, , , , , ,UX, , , , ,

ALLSEL,ALL
WPOFF,D

FLST,5,4,4,ORDE,2 
FITEM,5,365
FITEM,5,-368
LSEL,S, , ,P51X
NSLL,S,1
FLST,2,52,1,ORDE,2
FITEM,2,271
FITEM,2,-322
NROTAT,ALL
D,ALL, , , , , ,ALL, , , , , 
ALLSEL,ALL
!**************************************************************************
!求解
WPOFF,-D
CSYS,1
FINISH
/SOLU
!因为有可能初始位置不是刚好接触的，所以先旋转一个角度求解，然后去掉旋转自由度加上力矩再求解
FLST,5,4,4,ORDE,2 
FITEM,5,181
FITEM,5,-184
LSEL,S, , ,P51X
NSLL,S,1
FLST,2,52,1,ORDE,2
FITEM,2,1285
FITEM,2,-1336 
D,P51X, ,0.3, , , ,UY, , , , , !旋转一个角度，使接触面处于初始接触的位置
ALLSEL,ALL
SOLVE

FLST,5,4,4,ORDE,2 
FITEM,5,181
FITEM,5,-184
LSEL,S, , ,P51X
NSLL,S,1
FLST,2,52,1,ORDE,2
FITEM,2,1285
FITEM,2,-1336
DDELE,ALL,UY       !去掉第一次求解时加的旋转自由度
F,ALL,FY,100000/50   !给节点加上力矩，大小为：100000/50 N*MM
SOLVE
FINISH