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,169 !TARGE169
ET,3,172
MP,MU,2,0.2 !MU是设置摩擦因数
MAT,1 !定义单元材料属性指示器
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,1.2e11
MPDATA,PRXY,1,,0.3
R,31 !定义单元实常数
REAL,31 !定义单元实常数属性指示器
R,31,,,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 !这里的实常数设置在多处看到这样设置,可取
R,32
REAL,32
R,32,,,0.1,0.1,0,,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
R,33
REAL,33
R,33,,,0.1,0.1,0,,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
R,34
REAL,34
R,34,,,0.1,0.1,0,,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,3,4,0
KEYOPT,3,5,1 !All selected contact pairs are initially open. Rigid body motion can
!occur. You may use auto CNOF/ICONT by setting KEYOPT(5) to close
!small gaps.
NROPT,UNSYM
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0 !在接触分析中,单元选项的设置应该对CONTACT单元进行设置
KEYOPT,3,10,1
KEYOPT,3,11,0
KEYOPT,3,12,0
KEYOPT,3,2,0
KEYOPT,2,5,0
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
!-------------------选择需要覆盖接触单元的节点
!---------------FIRST CONTACT PAIR
*SET,_REALID,31
LSEL,S,,,59
CM,TARGET1,LINE
TYPE,2
NSLL,S,1
ESLN,S,0,ALL
ESURE,ALL
LSEL,S,,,156
CM,_CONTACT1,LINE
TYPE,3
NSLL,S,1
ESLN,S,0,ALL
ESURE,ALL
!---------------SECOND CONTACT PAIR
*SET,_REALID,32
LSEL,S,,,50
CM,TARGET2,LINE
TYPE,2
NSLL,S,1
ESLN,S,0,ALL
ESURE,ALL
LSEL,S,,,149
CM,_CONTACT2,LINE
TYPE,3
NSLL,S,1
ESLN,S,0,ALL
ESURE,ALL
!---------------THIRD CONTACT PAIR
*SET,_REALID,33
LSEL,S,,,3
CM,TARGET3,LINE
TYPE,2
NSLL,S,1
ESLN,S,0,ALL
ESURE,ALL
LSEL,S,,,10
CM,_CONTACT3,LINE
TYPE,3
NSLL,S,1
ESLN,S,0,ALL
ESURE,ALL
!---------------FORTH CONTACT PAIR
*SET,_REALID,34
LSEL,S,,,23
CM,TARGET4,LINE
TYPE,2
NSLL,S,1
ESLN,S,0,ALL
ESURE,ALL
LSEL,S,,,87
CM,_CONTACT4,LINE
TYPE,3
NSLL,S,1
ESLN,S,0,ALL
ESURE,ALL
FINISH
/SOLU
ANTYPE,STATIC
!----------------施加约束
CSWPLA,11,1
WPCSYS,-1,11
NSEL,S,LOC,X,20
D,ALL,UX
D,ALL,UY,-0.2 !一定要有这个约束,否则完全算不出来,我草,真鸡巴搞。主动齿轮还需要这个约束,他妈的接触约束不了主动齿轮的运动范围,真鸡巴搞
/PREP7
NROTAT,ALL !----旋转节点坐标系,使之与激活坐标系平行
FINISH !----如果不旋转,力的方向就会出错
/SOLU
F,ALL,FY,-200
!--------------------------------施加扭矩
CSYS,0
WPOFF,-216
CSWPLA,12,1
WPCSYS,-1,12
NSEL,S,LOC,X,150
D,ALL,ALL
!-----------------------------时间载荷步设置
NLGEOM,ON !打开大位移选项
NSUBST,20,200,1
AUTOTS,ON
KBC,0
TIME,1
ALLSEL,ALL
APLOT
CSYS,0
SOLVE
FINISH
/POST1
PLNSOL,S,EQV,0 |