求助啊求助啊，齿轮分析啊

jevons1212

小弟做了一个齿轮分析，APDL做的，我把.igs和APDL的文本文件都上传到附件里，大家帮我看看吧，反正最后的分析是没法进行下去了

求助啊，大家一定要帮帮我啊

我写了很久才搞出来这个，但是还是有问题。建模的时候看到的问题并不明显，但是显然在某些问题上，设置是有问题的。

jevons1212

我把命令流贴上来，大家看看

FINISH
/CLEAR,START

/FILENAME,gear_meshing
/TITLE,two_gears_contact analysis
WPSTYLE,,,,,,,,0
/AUX15  
!*  
IOPTN,IGES,NODEFEAT
IOPTN,MERGE,YES
IOPTN,SOLID,YES
IOPTN,SMALL,YES
IOPTN,GTOLER, DEFA  
IGESIN,'as','igs',' '   
LPLOT   
!*  

/PREP7  

ET,1,PLANE42
ET,2,SOLID95

TYPE,1

KL,-8,0.2,,
LANG,12,4,90                 !在线8上找一点做对面线的垂线
LSBL,8,14,,DELETE,KEEP        !切开线8
KL,1,0.2,,
LANG,5,9,90
LSBL,1,17,,DELETE,KEEP

WPCSYS,-1,0                !移动实体到合适位置所做的工作
!CSKP,11,CART,5,11,2,        !在小圆的圆心处定义坐标系11                移动元件的方法：定义一个新坐标系，回到原坐标系，在使用LTRAN等命令
!CSYS,4                        !回到工作平面坐标系
!LTRAN,11,ALL,,,,0,1        !转换，移动

LESIZE,7,,,5,,,,,1        !将各线段划分为几份，便于划分网格
LESIZE,14,,,5,,,,,1
LESIZE,13,,,8,,,,,1
LESIZE,15,,,8,,,,,1
LESIZE,16,,,5,,,,,1
LESIZE,12,,,10,,,,,1
LESIZE,9,,,2,,,,,1
LESIZE,10,,,10,,,,,1
LESIZE,11,,,3,,,,,1

LESIZE,6,,,5,,,,,1
LESIZE,17,,,5,,,,,1
LESIZE,5,,,5,,,,,1
LESIZE,19,,,5,,,,,1
LESIZE,8,,,10,,,,,1
LESIZE,18,,,5,,,,,1
LESIZE,4,,,3,,,,,1
LESIZE,3,,,10,,,,,1
LESIZE,2,,,2,,,,,1

AL,7,13,14,15
AL,10,11,12,14,16,9
AL,5,6,19,17
AL,3,4,8,17,18,2

LCCAT,2,3
LCCAT,17,8
LCCAT,9,10
LCCAT,12,14

MSHAPE,0,2D
MSHKEY,1
AMESH,1,2
AMESH,3,4

LSEL,R,LCCA
LDELE,ALL                !删除连接线（拉伸之前需要做）

*AFUN,RAD
WPROT,10

/TRIAD,LTOP
WPCSYS,-1,0
CSYS,4
WPROT,10
ASEL,S,,,3,4,,1                !选择面和一下的单元
ARSYM,Y,ALL,,,,0,0        !镜像
WPROT,10
ASEL,S,,,3,6,,1
ARSYM,Y,ALL,,,,0,0
WPCSYS,-1,0
ASEL,S,,,3,10,,1
ARSYM,Y,ALL,,,,0,0

WPOFFS,-216
WPROT,-2
ASEL,S,,,1,2,,1
ARSYM,Y,ALL,,,,0,0
ASEL,S,,,1,2,,1
ASEL,A,,,19,20,,
WPROT,-2
ARSYM,Y,ALL,,,,0,0
ASEL,S,,,1,2,,1
ASEL,A,,,19,24,,
WPROT,-4
ARSYM,Y,ALL,,,,0,0
WPROT,8
ASEL,S,,,1,2,,1
ASEL,A,,,19,32,,
ARSYM,Y,ALL,,,,0,0                !到这里，平面网格划分已经结束

NUMMRG,ALL,,,,LOW        !对实体进行合并
NUMCMP,ALL                !压缩编号
/VIEW,1,1,1,1
CSYS,0
WPCSYS,-1,0

K,,,,
K,,,,9,
K,,,,-9,
K,,,,8.5,
K,,,,-8.5,
LSTR,171,172
LSTR,171,174
LSTR,171,173
LSTR,171,175

ASEL,S,,,1,2,1
ASEL,A,,,19,48
CM,a1,AREA
ASEL,S,,,3,18
CM,a2,AREA

TYPE,2
ESIZE,4,0

CMSEL,S,a1
VDRAG,ALL,,,,,,204
CMSEL,S,a1
VDRAG,ALL,,,,,,202
CMSEL,S,a2
VDRAG,ALL,,,,,,203
CMSEL,S,a2
VDRAG,ALL,,,,,,201                !到这里，三维网格已经划分完毕

ESEL,S,TYPE,,1                        !选择类型1的单元
ACLEAR,ALL                        !删除2D单元

!-------------------选择需要覆盖接触单元的节点

ASEL,S,,,429
ASEL,A,,,517
ASEL,A,,,422
ASEL,A,,,510
ASEL,A,,,374
ASEL,A,,,462
ASEL,A,,,389
ASEL,A,,,477
NSLA,S,1
CM,TARGET_node,NODE
ALLSEL,ALL
ASEL,S,,,151
ASEL,A,,,311
ASEL,A,,,143
ASEL,A,,,303
ASEL,A,,,215
ASEL,A,,,55
ASEL,A,,,228
ASEL,A,,,68
NSLA,S,1
CM,CONTACT_node,NODE
ALLSEL,ALL
EPLOT


MP,MU,2,0.2        !MU是设置摩擦因数
MAT,1                !定义单元材料属性指示器
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,1.2e11
MPDATA,PRXY,1,,0.3

R,4                !定义单元实常数
REAL,4                !定义单元实常数属性指示器
ET,3,170        !TARGE170
ET,4,174        !CONTA174或者CONTA173


R,4,,,0.1,0.1,0,,                !这里定义为3或者4，最后的效果都是对接触对的单元进行设置实常数
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0        !这里的实常数设置在多处看到这样设置，可取


KEYOPT,4,4,0
KEYOPT,4,5,0
NROPT,UNSYM
KEYOPT,4,7,0
KEYOPT,4,8,0
KEYOPT,4,9,0                        !在接触分析中，单元选项的设置应该对CONTACT单元进行设置
KEYOPT,4,10,1
KEYOPT,4,11,0
KEYOPT,4,12,0
KEYOPT,4,2,0
KEYOPT,3,5,0

TYPE,3
NSEL,S,,,TARGET_node                !选择面上的节点，在节点上覆盖接触单元
ESLN,S,0,
ESURF,,,,

TYPE,4
NSEL,S,,,CONTACT_node                !选择面上的节点，在节点上覆盖接触单元
ESLN,S,0,
ESURF

ALLSEL
CSYS,0

!-------------------------------------------------------

FINISH
APLOT
/SOLU

ANTYPE,STATIC

!----------------施加约束
CSYS,0
WPCSYS,-1,0
DA,173,ALL
DA,186,ALL
DA,202,ALL
DA,195,ALL
DA,155,ALL
DA,167,ALL
DA,146,ALL
DA,139,ALL
DA,49,ALL
DA,63,ALL
DA,82,ALL
DA,74,ALL
DA,115,ALL
DA,129,ALL
DA,104,ALL
DA,96,ALL
DA,333,ALL
DA,346,ALL
DA,362,ALL
DA,355,ALL
DA,315,ALL
DA,327,ALL
DA,306,ALL
DA,299,ALL
DA,209,ALL
DA,223,ALL
DA,242,ALL
DA,234,ALL
DA,275,ALL
DA,289,ALL
DA,264,ALL
DA,256,ALL

CSWPLA,11,1
WPCSYS,-1,11

DA,526,UX
DA,534,UX
DA,515,UX
DA,501,UX
DA,460,UX
DA,468,UX
DA,493,UX
DA,479,UX
DA,438,UX
DA,446,UX
DA,427,UX
DA,413,UX
DA,372,UX
DA,380,UX
DA,405,UX
DA,391,UX
DA,526,UZ
DA,534,UZ
DA,515,UZ
DA,501,UZ
DA,460,UZ
DA,468,UZ
DA,493,UZ
DA,479,UZ
DA,438,UZ
DA,446,UZ
DA,427,UZ
DA,413,UZ
DA,372,UZ
DA,380,UZ
DA,405,UZ
DA,391,UZ
!--------------------------约束就这些

!--------------------------------施加扭矩

NSEL,S,LOC,X,20
/PREP7
NROTAT,ALL                !----旋转节点坐标系，使之与激活坐标系平行
FINISH                        !----如果不旋转，力的方向就会出错

/SOLU
F,ALL,FY,-200

!-----------------------------时间载荷步设置
NSUBST,10,100,2
AUTOTS,ON
TIME,4

ALLSEL,ALL
APLOT

NLGEOM,ON                !打开大位移选项

SOLVE
FINISH
/POST1