我最近做的也是接触问题,就是两个管接触 但是半径不同 属于间隙接触 一直没有作通啊 不知道怎样设置接触对关键字啊命令如下:
finish
/clear
/units,si !定义国际单位制
! ************* 初始化参数********************************
/prep7
et,1,solid45 !进行结构分析的体单元类型
mp,ex,1,2020 !定义材料的弹性模量,1号材料对应 换热管钢管
mp,nuxy,1,0.3 !定义材料的泊松比
mp,alpx,1,1.618E-5 !定义热膨胀系数
mp,ex,2,2e4 !定义材料的弹性模量,2号材料对应 管板锻件
mp,nuxy,2,0.3 !定义材料的泊松比
mp,alpx,2,1.655E-5 !定义热膨胀系数
mp,ex,3,2e5 !定义材料的弹性模量,3号材料对应 接触单元材料号
mp,nuxy,3,0.3 !定义材料的泊松比
mp,alpx,3,1.655E-5 !定义热膨胀系数
! **********建立模型**********
CSYS,1
CYLIND,50,60,-5,240,-20,20, !柔性体
CYLIND,60.2,70.2,0,100,-20,20, !刚性体
LSEL,S,LOC,X,50
LSEL,A,LOC,X,60
LSEL,A,LOC,X,60.2
LSEL,A,LOC,X,70.2
LSEL,U,LOC,Y,-20
LSEL,U,LOC,Y,20
LESIZE,ALL,,5 !圆弧线每份5度
ALLSEL,ALL
LSEL,S,LOC,Y,-20
LSEL,A,LOC,Y,20
LSEL,U,LOC,z,0
LSEL,U,LOC,z,100
LSEL,U,LOC,z,-5
LSEL,U,LOC,z,240
LESIZE,ALL,5 !轴向剖单元长度
ALLSEL,ALL
LSEL,S,LENGTH,,10
LESIZE,ALL,,,4 !径向剖分数
VSEL,S,LOC,X,50,60
TYPE,1
MAT,1 !柔性体
VMESH,ALL
ALLSEL,ALL
VSEL,S,LOC,X,60.2,70.2
TYPE,1
MAT,2 !刚性体
VMESH,ALL
ALLSEL,ALL
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
!----------------------
CSYS,1
allsel
NSEL,S,LOC,X,60
CM,N_cont,NODE !创建接触节点组
allsel
NSEL,S,LOC,X,60.2
CM,N_targ,NODE !创建目标节点组
allsel
NSEL,S,LOC,Y,-20
NSEL,A,LOC,Y,20
CM,N_A,NODE !创建约束面A节点组
allsel
NSEL,S,LOC,Z,0
NSEL,A,LOC,Z,100
NSEL,U,LOC,x,50,60
NSEL,A,LOC,Z,140
NSEL,A,LOC,Z,-5
CM,N_D,NODE !创建约束面D节点组
allsel
NSEL,S,LOC,X,70.2
CM,N_E,NODE !创建约束面E节点组
!=========================================建立接触对
ET,2,170 !目标单元
ET,3,173
MP,MU,3,0.3
MAT,3
R,3,!,,,,-0.15 !接触刚度 最大渗透范围 真实缝隙
REAL,3
KEYOPT,3,9,0
KEYOPT,3,10,2
R,3,!,,,,-0.15
RMORE,0.22,0.2
RMORE,
RMORE,,0
RMORE,0
!-------------! 以下生成目标面 实常数3单元类型2
NSEL,S,,,N_TARG
TYPE,2
ESLN,S,0
ESURF
!-------------! 以下生成接触面 实常数3单元类型3
NSEL,S,,,N_CONT
TYPE,3
ESLN,S,0
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,2
ESEL,A,TYPE,,3
ESEL,R,REAL,,3
/PSYMB,ESYS,1
EPLOT
!=========================================以下进行加载
CSYS,1
allsel
CMSEL,S,N_A
NROTAT,ALL !转换到节点坐标系 边界A
D,ALL,UY !边界A
allsel
CMSEL,S,N_D
D,ALL,UZ !边界D
CMSEL,S,N_E
NROTAT,ALL
D,ALL,UX !边界E
allsel
NSEL,S,LOC,X,50
NSEL,R,LOC,Z,0,100
SF,ALL,PRES,14 !施加胀接压力 胀管内壁边界条件B
/SOL
allsel !全选
solve !求解
finish
ALLSEL,ALL
!SAVE,'Example-solve','db'
/POST1
ALLSEL,ALL
RSYS,1
AVPRIN,0,0
AVRES,2,
/EFACET,1
LAYER,0
FORCE,TOTAL
/EFACET,1
PLNSOL, U,X, 0,1.0 |