找回密码
 注册
Simdroid-非首页
查看: 44|回复: 1

对称生成的单元,在哪一步要合并

[复制链接]
发表于 2011-3-29 17:26:00 | 显示全部楼层 |阅读模式 来自 湖北武汉
我用对称的方法生成的几个齿的单元,最后进行接触分析,得到如图的结果,一块和另外一块是完全分开的

请问在哪一步用什么命令可以解决这个问题


小弟先谢谢大家了,帮我看看吧

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
 楼主| 发表于 2011-3-29 17:26:18 | 显示全部楼层 来自 湖北武汉
Simdroid开发平台
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
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 16:46 , Processed in 0.027882 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表