齿轮接触分析有问题

jevons1212

我在做齿轮的接触分析中出问题了，但是不知道是哪里有问题，很是苦恼

其实在另外一个板块也发了一个类似的帖子了，但是实在是急得不行，希望大家能帮我一下

我把APDL语言的贴出来，大家帮我看一下

模型最开始是用IGS倒进来的，我把IGS的文件放在附件里

版主，如果这样发帖子不合适，麻烦您告诉我一下，也是太急了，所以才这样发帖子，的确是有重复的嫌疑，不过还是请版主手下留情啊

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