我这儿有个三维O型圈接触密封分析，有几个错误，望高手指点一下，万分感谢
浏览：28112
*** ERROR *** SUPPRESSED MESSAGE CP = 2639.000 TIME= 18:38:26
The value of UY at node 73232 is 130647280. It is greater than the
current limit of 1000000. This generally indicates rigid body motion
as a result of an unconstrained model. Verify that your model is
properly constrained.
*** ERROR *** SUPPRESSED MESSAGE CP = 2639.000 TIME= 18:38:26
*** MESSAGE CONTINUATION ---- DIAGNOSTIC INFORMATION ***
If one or more parts of the model are held together only by contact
verify that the contact surfaces are closed. You can check contact
status in the SOLUTION module for the converged solutions using
CNCHECK.
*** ERROR *** SUPPRESSED MESSAGE CP = 2639.016 TIME= 18:38:26
*** MESSAGE CONTINUATION ---- DIAGNOSTIC INFORMATION ***
Rigid body motion can also occur when net section yielding has
occurred resulting in large displacements for small increments of load
or when buckling has occurred. You can plot the time history curve
for node 73232 in the UY direction to check for stiffness (slope of
the curve) approaching zero.
*** WARNING *** SUPPRESSED MESSAGE CP = 2745.266 TIME= 18:40:59
Coefficient ratio exceeds 1.0e8 - Check results.
*** WARNING *** SUPPRESSED MESSAGE CP = 3178.703 TIME= 18:53:37
Contact element 199152 (real ID 8) status changes abruptly from
no-contact -> contact (with target element 192982).
*** WARNING *** SUPPRESSED MESSAGE CP = 4931.203 TIME= 19:33:40
There are 6321 small equation solver pivot terms.
*** WARNING *** SUPPRESSED MESSAGE CP = 5565.453 TIME= 19:50:23
There are 26379 small equation solver pivot terms.
*** WARNING *** SUPPRESSED MESSAGE CP = 6280.562 TIME= 20:09:14
There are 148640 small equation solver pivot terms.
*** WARNING *** SUPPRESSED MESSAGE CP = 6956.484 TIME= 20:26:37
There are 146334 small equation solver pivot terms.
*** WARNING *** SUPPRESSED MESSAGE CP = 7083.500 TIME= 20:30:51
Coefficient ratio exceeds 1.0e8 - Check results.
命令流如下
finish
/clear
/filn,tank
/title,analysis of non-circular tank
*set,b1,1850 ! 平盖长度
*set,b2,1440 ! 平盖宽度
*set,s1,50 ! 平盖厚度
*set,b3,1850 ! 法兰长度
*set,b4,1440 ! 法兰宽度
*set,s2,80 ! 法兰厚度
*set,h1,130 ! 法兰高度
*set,b5,1850 ! 筒体长度
*set,b6,1440 ! 筒体宽度
*set,h2,895 ! 筒体高度
*set,s3,40 ! 筒体厚度
*set,l1,1107 ! 筒体长度方向两个焊缝之间的距离
*set,l2,1321 ! 筒体宽度方向两个焊缝之间的距离
*set,rr1,160 ! 筒体内侧过渡圆角半径
*set,rr2,190
*set,r1,885 ! 筒体内圆半径
*set,r2,925 ! 筒体外圆半径
*set,rr3,200 ! 平盖过渡圆角半径
*set,d1,600 ! 平盖管口一的内径
*set,d2,710 ! 平盖管口一的外径
*set,t1,2 ! 管口一的外圆深度
*set,d3,650 ! 平盖接管螺栓中心圆直径
*set,c1,5.7 ! 密封槽的深度
*set,c2,10 ! 密封槽的宽度
*set,h3,650 ! 密封槽内侧距中心距的高度
*set,d4,1630 ! 平盖底部（与法兰连接处）内伸直径
*set,t2,2 ! 平盖底部（与法兰连接处）内伸高度
*set,n,72 ! 双头螺柱个数
*set,l,120 ! 双头螺柱长度
*set,d5,17 ! 双头螺柱的螺栓孔直径
*set,d6,16 ! 螺栓直径
*set,nut_rad,27.7/2 ! 螺母半径
*set,nut_h,13 ! 螺母高度
*set,d7,1740 ! 密封槽的内侧直径
*set,e1,400 ! 接口二距一口的方向坐标位移
*set,e2,420 ! 接口二距一口的方向坐标位移
*set,e3,1039 ! 密封槽内径过渡圆角坐标方向距离
*set,e4,920 ! 密封槽内径过渡圆角坐标方向距离
*set,b7,1110 ! 平盖长度方向直边长度
*set,a1,1000 ! 平盖底部折边长度
*set,a2,106.4 ! 平盖底部折边长度
*set,h4,608 ! 平盖底部折边高度
*set,h5,579.49 ! 平盖底部折边高度
*set,h6,1220/2 ! 平盖宽度方向最高点
*set,a3,107.54 ! 上法兰底部折边长度
*set,h7,609 ! 上法兰底部折边高度
*set,h8,580.18 ! 上法兰底部折边高度
*set,h9,1240/2 ! 上法兰宽度方向最高点
*set,rr5,50 ! 上法兰内径过渡圆角
*set,d8,1650 ! 上法兰内圆直径
*set,a4,1770 ! 上法兰长度方向内径长度（与筒体连接处）
*set,a5,1670 ! 上法兰长度方向倒角中心距距离
*set,a6,1360 ! 上法兰宽度方向内径长度（与筒体连接处）
*set,a7,1260 ! 上法兰宽度方向倒角中心距距离
*set,d9,1390 ! 上法兰宽度方向螺栓中心距距离
*set,d10,1800 ! 上法兰长度方向螺栓中心距距离
*set,a8,1360 ! 上法兰宽度方向内径长度（与平盖连接处）
*set,e5,500 ! 上法兰外径过渡圆角坐标的方向位移
*set,e6,510
*set,d11,7 ! 垫片直径
*set,rr6,110
*afun,deg ! 角度单位设定
*set,al1,30
*set,al2,2
*set,al3,75
*set,pi,0.3 ! 设计压力（内压）
*set,FORC_BOLT,86260.4
*set,py,2.60 ! 端面载荷
/prep7
et,1,solid45 ! 实体单元
et,2,HYPER58
!定义螺母材料1
mp,ex,1,1.92e5 ! 设定弹性模量
mp,nuxy,1,0.27 ! 定义泊松比
mp,dens,1,7.85e-9 ! 材料的密度
!定义双头螺柱材料2
mp,ex,2,2.09e5 ! 设定弹性模量
mp,nuxy,2,0.269 ! 定义泊松比
mp,dens,2,7.85e-9 ! 材料的密度
!定义平盖、筒体材料3
mp,ex,3,2.0E5 ! 设定弹性模量
mp,nuxy,3,0.25 ! 定义泊松比
mp,dens,3,7.85e-9 ! 材料的密度
!定义法兰材料4
mp,ex,4,2.06E5 ! 设定弹性模量
mp,nuxy,4,0.3 ! 定义泊松比
mp,dens,4,7.85e-9 ! 材料的密度
!定义O形密封圈材料5
mp,ex,5, ! 设定弹性模量
mp,nuxy,5,0.4997 ! 定义泊松比
TB,MOONEY,5,1
TBDATA,1,0.5516
TBDATA,2,0.1739
!定义摩擦系数为材料6
mp,mu,6,0.3
!定义摩擦系数为材料7
mp,mu,7,0.25
! 创建模型
!创建平盖模型
k,1,0,0,0
k,2,0,0,-h6
k,3,e5,0,-h4
k,4,e5+a2,0,-h5
k,5,e5+a2+h5*tan(al1),0,0
k,6,d4/2,0,0
k,7,0,0,-d4/2
l,1,2
l,2,3
l,3,4
l,4,5
larc,6,7,1,d4/2
lptn,4,5
ldele,7,,,1
ldele,9,,,1
ldele,11,,,1
l,1,6
lcomb,6,10
al,1,2,3,6,8,4
k,5,0,0,-h3
k,7,e3/2,0,-h3
l,2,5
l,5,7
k,8,e3/2,0,-e4/2
l,1,9
cskp,11,0,9,1,4
lgen,2,6,,,-(d7/2-d4/2),,,,,0
lgen,2,6,,,-(d7/2-d4/2+c2),,,,,0
lgen,2,6,,,-(b1/2-d4/2),,,,,0
csdele,all
k,16,d7/2,0,0
larc,11,16,1,d7/2
l,6,16
k,17,d7/2+c2,0,0
larc,13,17,1,d7/2+c2
l,16,17
k,18,b1/2,0,0
larc,15,18,1,b1/2
l,17,18
ldele,9
k,19,e3/2+rr2,0,-e4/2
larc,7,19,8,rr2
lovlap,9,10
ldele,20,,,1
ldele,22,,,1
ldele,23,,,1
lcomb,19,24
k,10,0,0,-h3-c2
k,19,e3/2,0,-h3-c2
l,5,10
l,10,19
k,20,e3/2+rr3,0,-e4/2
larc,19,20,8,rr3
lovlap,11,20
ldele,22,,,1
ldele,25,,,1
ldele,27,,,1
lcomb,24,26
k,12,0,0,-b2/2
k,20,b1/2,0,-b2/2+b1/2*tan(al2)
l,10,12
l,12,20
lextnd,12,14,100
lovlap,12,20
ldele,25,,,1
ldele,27,,,1
lfillt,22,26,rr3
al,2,5,7,19,21,13,14,8,6,3
al,9,10,24,23,15,16,13,21,19,7
al,11,22,12,26,17,18,15,23,24,10
asel,s,,,1,4,1
vext,all,,,,s1
wpoff,0,t2
wprota,0,90
vsel,s,,,2,4,1
vsbw,all
vdele,5,,,1
vdele,6,,,1
vdele,7,,,1
wpoff,0,0,-c1
vsbw,9
vdele,3,,,1
wpstyl,defa
wprota,0,-90
allsel
cyl4,0,0,d1/2,,,,s1
vsbv,1,3
wpoff,0,0,s1
cyl4,0,0,d2/2,,,,-s1
vovlap,1,4
vdele,5,,,1
wpoff,0,0,-t1
vsbw,3
vdele,4,,,1
wpstyl,defa
!创建上法兰
lgen,2,108,111,3,,-t2,,,,0
k,71,0,0,0
cskp,11,0,23,71,70
lgen,2,68,,,(b3/2-d8/2),,,,,0
csdele,all
k,74,d8/2,0,0
larc,73,74,71,d8/2
l,74,65
k,75,0,0,-h9
k,76,0,0,-b4/2
l,75,76
k,77,b3/2,0,-b4/2+b3/2*tan(al2)
l,76,77
lextnd,68,70,200
lovlap,68,94
ldele,104,,,1
ldele,107,,,1
lfillt,106,95,rr3
k,78,e6,0,-h7
k,79,e6+a3,0,-h8
l,75,78
l,78,79
l,72,79
lcomb,69,107
al,94,104,69,79,83,66,95,68,106,91
vext,22,,,,-s2
NUMMRG,ALL
NUMCMP,ALL
kdele,71
wpoff,0,-s2
k,71,0,-s2,0
cskp,11,0,85,71,86
lgen,2,122,,,(b3/2-a5/2),,,,,0
csdele,all
lgen,2,124,,,,,b4/2-a7/2,,,0
lextnd,141,90,200
lextnd,142,92,200
lovlap,141,142
ldele,144,,,1
ldele,146,,,1
lfillt,143,145,rr6
k,94,a5/2,-s2,0
larc,89,94,71,a5/2
csys,4
k,95,0,0,-b4/2
k,96,554.9057,,-700.6223
k,97,721.1309,,-600.7441
k,98,808.1206,,-450.0735
k,99,b3/2,0,0
l,91,95
l,95,96
l,97,98
larc,98,99,71,b3/2
lextnd,146,96,200
lextnd,147,97,200
lovlap,146,147
ldele,150,,,1
ldele,151,,,1
lfillt,149,152,rr3
l,94,99
al,145,141,143,142,147,148,152,146,149,144
vext,64,,,,-(h1-s2)
kwpave,105
cyl4,0,0,rr5
vdrag,76,,,,,,154,153,151,150
kwpave,101
vext,96,,,-100
vovlap,7,8,9,10,11,12
vdele,13,,,1
vdele,14,,,1
vdele,15,,,1
vdele,16,,,1
vdele,17,,,1
vdele,21,,,1
vdele,20,,,1
vdele,19,,,1
vdele,18,,,1
vdele,22,,,1
wpstyl,defa
vptn,3,23
!创建筒体
vext,135,,,,-h2
!创建密封圈
wpoff,d7/2+c2/2,0,0
wpoff,0,d11/2,0
cyl4,0,0,d11/2
vdrag,83,,,,,,9,10,12,14
wpstyl,defa
vext,103,,,-50
vadd,12,13
wprota,0,0,90
vsbw,14
vdele,12,,,1
wpstyl,defa
!创建螺栓孔
wprota,0,90
cyl4,38,-675,d5/2,,,,-86
vgen,7,12,,,76,,,,,0
cyl4,570,-670,d5/2,,,,-86
cyl4,641,-645,d5/2,,,,-86
cyl4,696,-594,d5/2,,,,-86
cyl4,734,-529,d5/2,,,,-86
cyl4,770,-464,d5/2,,,,-86
cyl4,806,-399,d5/2,,,,-86
cyl4,837,-330,d5/2,,,,-86
cyl4,861,-259,d5/2,,,,-86
cyl4,880,-187,d5/2,,,,-86
cyl4,892,-113,d5/2,,,,-86
cyl4,899,-38,d5/2,,,,-86
vsel,s,,,12
vsel,a,,,14,30,1
cm,bolt,volu
allsel
vsbv,6,bolt
!创建螺栓
cyl4,38,-675,d6/2,,,,-86
vgen,7,6,,,76,,,,,0
cyl4,570,-670,d6/2,,,,-86
cyl4,641,-645,d6/2,,,,-86
cyl4,696,-594,d6/2,,,,-86
cyl4,734,-529,d6/2,,,,-86
cyl4,770,-464,d6/2,,,,-86
cyl4,806,-399,d6/2,,,,-86
cyl4,837,-330,d6/2,,,,-86
cyl4,861,-259,d6/2,,,,-86
cyl4,880,-187,d6/2,,,,-86
cyl4,892,-113,d6/2,,,,-86
cyl4,899,-38,d6/2,,,,-86
!创建螺母
wpoff,0,0,-s1
rpr4,6,38,-675,nut_rad,,-nut_h
vgen,7,30,,,76,,,,,0
rpr4,6,570,-670,nut_rad,,-nut_h
rpr4,6,641,-645,nut_rad,,-nut_h
rpr4,6,696,-594,nut_rad,,-nut_h
rpr4,6,734,-529,nut_rad,,-nut_h
rpr4,6,770,-464,nut_rad,,-nut_h
rpr4,6,806,-399,nut_rad,,-nut_h
rpr4,6,837,-330,nut_rad,,-nut_h
rpr4,6,861,-259,nut_rad,,-nut_h
rpr4,6,880,-187,nut_rad,,-nut_h
rpr4,6,892,-113,nut_rad,,-nut_h
rpr4,6,899,-38,nut_rad,,-nut_h
!视镜孔处法兰联接螺栓的当量面（垫片）
wpoff,0,0,t2
cyl4,,,d3/2-d6/2,,d3/2+d6/2,-90,s1-t2
vovlap,1,49
wpstyl,defa
allsel
NUMMRG,ALL
NUMCMP,ALL
!划分网格
type,1
vsel,s,,,5
vsel,a,,,11
vsel,a,,,13,28,1
vsel,a,,,29
vsel,a,,,31,47,1
vptn,all
mat,1
vsel,s,,,105,122,1
cm,nut,volu
esize,6
vsweep,all
! 螺母划分
mat,2
vsel,s,,,51,68,1
vsweep,all
esize,12
vsel,s,,,69,103,2
vsweep,all
esize,18
vsel,s,,,70,104,2
cm,nut_bot,volu
vsweep,all
allsel
mat,3
!kwpave,26
!wprota,0,0,90
!vsbw,3
!kwpave,1
!vsbw,5
!wpstyl,defa
wpoff,0,c1+2
wprota,0,90
vsbw,4
vsbw,30
wpstyl,defa
kwpave,655
wprota,0,90
vsel,s,,,3
!vsel,a,,,12,13,1
vsel,a,,,48,50,1
cm,flat_head,volu
vsbw,flat_head
wpstyl,defa
lesize,111,,,6
lesize,34,,,6
lesize,55,,,6
lesize,33,,,6
vsel,s,,,4,5,1
vsel,a,,,11,13,2
vsel,a,,,1
vsweep,all
vsel,s,,,20,21,1
vsweep,all
vsel,s,,,14,19,1
vsweep,all
allsel
mat,4
vext,14,,,,-100
vovlap,7,3
vdele,23,,,1
vdele,22,,,1
vext,75,,,,100
vovlap,26,3
vdele,7,,,1
lesize,131,,,8
lesize,132,,,8
lesize,1181,,,6
vsel,s,,,22,25,1
vsweep,all
VSEL,S,,,6
VPLOT
ASLV
LSLA
VDRAG,65,,,,,,158
VDRAG,41,,,,,,118
VDRAG,205,,,,,,119
VDRAG,210,,,,,,120
BTOL,0.1
VSBA,27,64
VDELE,28,,,1
VDELE,6,,,1
VGLUE,ALL
lesize,232,,,6
lesize,222,,,6
lesize,157,,,6
vsel,s,,,6
vsel,a,,,27,28,1
vsel,a,,,30
vsweep,all
mat,3
vsel,s,,,2
vsweep,all
type,2
mat,5
esize,1
vsel,s,,,8,10,1
vsel,a,,,12
cm,O_seal,volu
vsweep,all
type,1
! 创建螺栓预紧单元
ALLSEL
MAT,2
PSMESH,1,BOLT1,,VOLU,70,0,Y,s1/2,,,,BOLT_1
PSMESH,2,BOLT1,,VOLU,72,0,Y,s1/2,,,,BOLT_2
PSMESH,3,BOLT1,,VOLU,74,0,Y,s1/2,,,,BOLT_3
PSMESH,4,BOLT1,,VOLU,76,0,Y,s1/2,,,,BOLT_4
PSMESH,5,BOLT1,,VOLU,78,0,Y,s1/2,,,,BOLT_5
PSMESH,6,BOLT1,,VOLU,80,0,Y,s1/2,,,,BOLT_6
PSMESH,7,BOLT1,,VOLU,82,0,Y,s1/2,,,,BOLT_7
PSMESH,8,BOLT1,,VOLU,84,0,Y,s1/2,,,,BOLT_8
PSMESH,9,BOLT1,,VOLU,86,0,Y,s1/2,,,,BOLT_9
PSMESH,10,BOLT1,,VOLU,88,0,Y,s1/2,,,,BOLT_10
PSMESH,11,BOLT1,,VOLU,90,0,Y,s1/2,,,,BOLT_11
PSMESH,12,BOLT1,,VOLU,92,0,Y,s1/2,,,,BOLT_12
PSMESH,13,BOLT1,,VOLU,94,0,Y,s1/2,,,,BOLT_13
PSMESH,14,BOLT1,,VOLU,96,0,Y,s1/2,,,,BOLT_14
PSMESH,15,BOLT1,,VOLU,98,0,Y,s1/2,,,,BOLT_15
PSMESH,16,BOLT1,,VOLU,100,0,Y,s1/2,,,,BOLT_16
PSMESH,17,BOLT1,,VOLU,102,0,Y,s1/2,,,,BOLT_17
PSMESH,18,BOLT1,,VOLU,104,0,Y,s1/2,,,,BOLT_18
! 建立约束方程
! 螺栓和法兰的约束
ALLSEL
CMSEL,S,nut_bot
ASLV
ASEL,R,LOC,Y,0
cm,acont,area
vsel,s,,,25
aslv
Asel,r,loc,y,0
nsla
esln
cm,econ,elem
cmsel,s,acont
nsla
ceintf
! 上法兰两端面的约束
ALLSEL
vsel,s,,,22,25,1
ASLV
ASEL,R,LOC,Y,-s2
cm,acont1,area
vsel,s,,,6
vsel,a,,,27,28,1
vsel,a,,,30
aslv
Asel,r,loc,y,-s2
nsla
esln
cm,econ1,elem
cmsel,s,acont1
nsla
ceintf
! 法兰下端和筒体的约束
ALLSEL
vsel,s,,,6
vsel,a,,,27,28,1
vsel,a,,,30
ASLV
ASEL,R,LOC,Y,-h1
cm,acont2,area
vsel,s,,,2
aslv
Asel,r,loc,y,-h1
nsla
esln
cm,econ2,elem
cmsel,s,acont2
nsla
ceintf
! 建立接触
! 建立螺母和平盖的接触
allsel
CMSEL,S,nut ! 选择螺母体积组件
ASLV ! 选择依附于体上的面
ASEL,R,LOC,Y,s1 ! 通过坐标选择螺栓上的接触面
CM,AS_NUT,AREA ! 创建螺母与平盖的接触面组件
ALLSEL
VSEL,S,,,4 ! 选择与螺母接触的平盖体积
ASLV ! 选择依附于体上的面
ASEL,R,LOC,Y,s1 ! 通过坐标选择平盖上的目标面
CM,AS_FLAT1,AREA ! 创建与螺母接触的平盖目标面
et,4,targe170 ! 3D接触目标面单元170
et,5,conta174 ! 3D接触接触面单元174
KEYOPT,5,5,1
KEYOPT,5,9,0
KEYOPT,5,12,2
R,4
REAL,4
MAT,6
CMSEL,S,AS_FLAT1
TYPE,4 ! 激活目标单元类型
NSLA,S,1 ! 选择目标面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建目标单元
ALLSEL
CMSEL,S,AS_NUT
TYPE,5 ! 激活接触单元类型
NSLA,S,1 ! 选择接触面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建接触单元
allsel
! 建立o形圈和沟槽左侧的接触
ALLSEL
VSEL,S,,,12 ! 选择与密封圈接触的沟槽体积
ASLV ! 选择依附于体上的面
ASEL,R,,,15,16,1
ASEL,a,,,18,19,1 ! 选择沟槽的目标面
CM,AS_groove1,AREA ! 创建与密封圈接触的沟槽目标面
allsel
CMSEL,S,O_seal ! 选择密封圈体积组件
ASLV ! 选择依附于体上的面
ASEL,R,,,81,82,1
ASEL,a,,,86,87,1
ASEL,a,,,91,92,1
ASEL,a,,,96,97,1 ! 选择密封圈上的接触面
CM,AS_seal1,AREA ! 创建密封圈与沟槽的接触面组件
et,6,targe170 ! 3D接触目标面单元170
et,7,conta174 ! 3D接触接触面单元174
KEYOPT,7,5,1
KEYOPT,7,9,0
R,5
REAL,5
MAT,7
CMSEL,S,AS_groove1
TYPE,6 ! 激活目标单元类型
NSLA,S,1 ! 选择目标面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建目标单元
ALLSEL
CMSEL,S,AS_seal1
TYPE,7 ! 激活接触单元类型
NSLA,S,1 ! 选择接触面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建接触单元
allsel
! 建立o形圈和沟槽右侧的接触
CMSEL,S,O_seal ! 选择密封圈体积组件
ASLV ! 选择依附于体上的面
ASEL,R,,,80,83,3
ASEL,a,,,85,88,3
ASEL,a,,,90,93,3
ASEL,a,,,95,98,3 ! 选择密封圈上的接触面
CM,AS_seal2,AREA ! 创建密封圈与沟槽的接触面组件
ALLSEL
VSEL,S,,,13 ! 选择与密封圈接触的沟槽体积
ASLV ! 选择依附于体上的面
ASEL,R,,,21,23,1
ASEL,a,,,25 ! 选择沟槽的目标面
CM,AS_groove2,AREA ! 创建与密封圈接触的沟槽目标面
et,8,targe170 ! 3D接触目标面单元170
et,9,conta174 ! 3D接触接触面单元174
KEYOPT,9,5,1
KEYOPT,9,9,0
R,6
REAL,6
MAT,7
CMSEL,S,AS_groove2
TYPE,8 ! 激活目标单元类型
NSLA,S,1 ! 选择目标面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建目标单元
ALLSEL
CMSEL,S,AS_seal2
TYPE,9 ! 激活接触单元类型
NSLA,S,1 ! 选择接触面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建接触单元
allsel
! 建立o形圈和平盖的接触
CMSEL,S,O_seal ! 选择密封圈体积组件
ASLV ! 选择依附于体上的面
ASEL,R,,,80,81,1
ASEL,a,,,85,86,1
ASEL,a,,,90,91,1
ASEL,a,,,95,96,1 ! 选择密封圈上的接触面
CM,AS_seal3,AREA ! 创建密封圈与平盖的接触面组件
ALLSEL
VSEL,S,,,1 ! 选择与密封圈接触的平盖体积
ASLV ! 选择依附于体上的面
ASEL,R,loc,y,c1+t2 ! 选择平盖的目标面
CM,AS_groove3,AREA ! 创建与密封圈接触的平盖目标面
et,10,targe170 ! 3D接触目标面单元170
et,11,conta174 ! 3D接触接触面单元174
KEYOPT,11,5,1
KEYOPT,11,9,0
R,7
REAL,7
MAT,7
CMSEL,S,AS_groove3
TYPE,10 ! 激活目标单元类型
NSLA,S,1 ! 选择目标面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建目标单元
ALLSEL
CMSEL,S,AS_seal3
TYPE,11 ! 激活接触单元类型
NSLA,S,1 ! 选择接触面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建接触单元
allsel
! 建立o形圈和沟槽内表面的接触
CMSEL,S,O_seal ! 选择密封圈体积组件
ASLV ! 选择依附于体上的面
ASEL,R,,,82,83,1
ASEL,a,,,87,88,1
ASEL,a,,,92,93,1
ASEL,a,,,97,98,1 ! 选择密封圈上的接触面
CM,AS_seal4,AREA ! 创建密封圈与平盖的接触面组件
ALLSEL
VSEL,S,,,24 ! 选择与密封圈接触的平盖体积
ASLV ! 选择依附于体上的面
ASEL,R,loc,y,0 ! 选择平盖的目标面
CM,AS_groove4,AREA ! 创建与密封圈接触的平盖目标面
et,12,targe170 ! 3D接触目标面单元170
et,13,conta174 ! 3D接触接触面单元174
KEYOPT,13,5,1
KEYOPT,13,9,0
R,8
REAL,8
MAT,7
CMSEL,S,AS_groove4
TYPE,12 ! 激活目标单元类型
NSLA,S,1 ! 选择目标面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建目标单元
ALLSEL
CMSEL,S,AS_seal4
TYPE,13 ! 激活接触单元类型
NSLA,S,1 ! 选择接触面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建接触单元
allsel
! 求解设置
!NSEL,S,LOC,Y,-h1-h2
!NROTAT,ALL
!CSYS,1
!D,ALL,UX
/SOLU
! 施加预紧力
CUTCONTROL,PLSLIMIT,10,
CUTCONTROL,CRPLIMITexp,0.2,0
CUTCONTROL,CRPLIMITimp,0.25,1
LSCLEAR,ALL ! 清空所有载荷
ALLSEL
CSYS,0 ! 激活笛卡尔坐标系
ASEL,S,LOC,x,0 ! 通过坐标选择为0的面
DA,ALL,SYMM ! 施加对称边界条件
ASEL,S,LOC,z,0 ! 通过坐标选择为0的面
DA,ALL,SYMM ! 施加对称边界条件
nSEL,S,LOC,y,-h1-h2 ! 通过坐标选择面
D,ALL,uy ! 施加Y向位移约束
nsel,s,,,7028
D,ALL,uy
NROPT,FULL ! 激活完全NEWTON-RAPHSON
PRED,ON ! 打开变形预测
LNSRCH,on ! 激活线性搜索
nlgeom,on
ANTYPE,0 ! 静态分析
AUTOTS,1 ! 自动时间步长
NSUBST,10,20,1 ! 设置载荷步
OUTRES,ALL,ALL ! 输出所有结果
ALLSEL
SLOAD,1,9,LOCK,FORC,FORC_BOLT,1,2 ! 在螺栓中施加螺栓预紧
SLOAD,2,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,3,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,4,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,5,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,6,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,7,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,8,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,9,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,10,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,11,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,12,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,13,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,14,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,15,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,16,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,17,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,18,9,LOCK,FORC,FORC_BOLT,1,2
ALLSEL
TIME,10 ! 设置时间步至10
LSWRITE,1 ! 写出第一个载荷步文件
!对容器壁施加内压
lsel,s,radius,,d1/2
asll
asel,u,loc,y,0
asel,u,loc,y,c1+t2
asel,u,loc,y,s1-t1
asel,a,,,362,372,10
asel,a,,,15,16,1
asel,a,,,364,368,4
asel,a,,,43
ASEL,A,,,30,34,2
ASEL,A,,,29
asel,a,,,39
CMSEL,a,AS_seal1
CMSEL,a,AS_groove1
ASEL,A,,,62,63,1
ASEL,A,,,66
ASEL,A,,,100
ASEL,A,,,27,28,1
ASEL,A,,,33
ASEL,A,,,36,37,1
ASEL,A,,,46
ASEL,A,,,60
ASEL,A,,,75
! 选中施加内压的面
SFA,ALL,1,PRES,Pc ! 对内壁施加压强,大小为Pc
ASEL,S,,,369
SFA,ALL,1,PRES,-PY ! 施加端面当量面载荷
ALLSEL
TIME,20 ! 设置时间步至20
LSWRITE,2 ! 写出第二个载荷步文件
LSSOLVE,1,2,1 ! 计算载荷步1-2
The value of UY at node 73232 is 130647280. It is greater than the
current limit of 1000000. This generally indicates rigid body motion
as a result of an unconstrained model. Verify that your model is
properly constrained.
*** ERROR *** SUPPRESSED MESSAGE CP = 2639.000 TIME= 18:38:26
*** MESSAGE CONTINUATION ---- DIAGNOSTIC INFORMATION ***
If one or more parts of the model are held together only by contact
verify that the contact surfaces are closed. You can check contact
status in the SOLUTION module for the converged solutions using
CNCHECK.
*** ERROR *** SUPPRESSED MESSAGE CP = 2639.016 TIME= 18:38:26
*** MESSAGE CONTINUATION ---- DIAGNOSTIC INFORMATION ***
Rigid body motion can also occur when net section yielding has
occurred resulting in large displacements for small increments of load
or when buckling has occurred. You can plot the time history curve
for node 73232 in the UY direction to check for stiffness (slope of
the curve) approaching zero.
*** WARNING *** SUPPRESSED MESSAGE CP = 2745.266 TIME= 18:40:59
Coefficient ratio exceeds 1.0e8 - Check results.
*** WARNING *** SUPPRESSED MESSAGE CP = 3178.703 TIME= 18:53:37
Contact element 199152 (real ID 8) status changes abruptly from
no-contact -> contact (with target element 192982).
*** WARNING *** SUPPRESSED MESSAGE CP = 4931.203 TIME= 19:33:40
There are 6321 small equation solver pivot terms.
*** WARNING *** SUPPRESSED MESSAGE CP = 5565.453 TIME= 19:50:23
There are 26379 small equation solver pivot terms.
*** WARNING *** SUPPRESSED MESSAGE CP = 6280.562 TIME= 20:09:14
There are 148640 small equation solver pivot terms.
*** WARNING *** SUPPRESSED MESSAGE CP = 6956.484 TIME= 20:26:37
There are 146334 small equation solver pivot terms.
*** WARNING *** SUPPRESSED MESSAGE CP = 7083.500 TIME= 20:30:51
Coefficient ratio exceeds 1.0e8 - Check results.
命令流如下
finish
/clear
/filn,tank
/title,analysis of non-circular tank
*set,b1,1850 ! 平盖长度
*set,b2,1440 ! 平盖宽度
*set,s1,50 ! 平盖厚度
*set,b3,1850 ! 法兰长度
*set,b4,1440 ! 法兰宽度
*set,s2,80 ! 法兰厚度
*set,h1,130 ! 法兰高度
*set,b5,1850 ! 筒体长度
*set,b6,1440 ! 筒体宽度
*set,h2,895 ! 筒体高度
*set,s3,40 ! 筒体厚度
*set,l1,1107 ! 筒体长度方向两个焊缝之间的距离
*set,l2,1321 ! 筒体宽度方向两个焊缝之间的距离
*set,rr1,160 ! 筒体内侧过渡圆角半径
*set,rr2,190
*set,r1,885 ! 筒体内圆半径
*set,r2,925 ! 筒体外圆半径
*set,rr3,200 ! 平盖过渡圆角半径
*set,d1,600 ! 平盖管口一的内径
*set,d2,710 ! 平盖管口一的外径
*set,t1,2 ! 管口一的外圆深度
*set,d3,650 ! 平盖接管螺栓中心圆直径
*set,c1,5.7 ! 密封槽的深度
*set,c2,10 ! 密封槽的宽度
*set,h3,650 ! 密封槽内侧距中心距的高度
*set,d4,1630 ! 平盖底部（与法兰连接处）内伸直径
*set,t2,2 ! 平盖底部（与法兰连接处）内伸高度
*set,n,72 ! 双头螺柱个数
*set,l,120 ! 双头螺柱长度
*set,d5,17 ! 双头螺柱的螺栓孔直径
*set,d6,16 ! 螺栓直径
*set,nut_rad,27.7/2 ! 螺母半径
*set,nut_h,13 ! 螺母高度
*set,d7,1740 ! 密封槽的内侧直径
*set,e1,400 ! 接口二距一口的方向坐标位移
*set,e2,420 ! 接口二距一口的方向坐标位移
*set,e3,1039 ! 密封槽内径过渡圆角坐标方向距离
*set,e4,920 ! 密封槽内径过渡圆角坐标方向距离
*set,b7,1110 ! 平盖长度方向直边长度
*set,a1,1000 ! 平盖底部折边长度
*set,a2,106.4 ! 平盖底部折边长度
*set,h4,608 ! 平盖底部折边高度
*set,h5,579.49 ! 平盖底部折边高度
*set,h6,1220/2 ! 平盖宽度方向最高点
*set,a3,107.54 ! 上法兰底部折边长度
*set,h7,609 ! 上法兰底部折边高度
*set,h8,580.18 ! 上法兰底部折边高度
*set,h9,1240/2 ! 上法兰宽度方向最高点
*set,rr5,50 ! 上法兰内径过渡圆角
*set,d8,1650 ! 上法兰内圆直径
*set,a4,1770 ! 上法兰长度方向内径长度（与筒体连接处）
*set,a5,1670 ! 上法兰长度方向倒角中心距距离
*set,a6,1360 ! 上法兰宽度方向内径长度（与筒体连接处）
*set,a7,1260 ! 上法兰宽度方向倒角中心距距离
*set,d9,1390 ! 上法兰宽度方向螺栓中心距距离
*set,d10,1800 ! 上法兰长度方向螺栓中心距距离
*set,a8,1360 ! 上法兰宽度方向内径长度（与平盖连接处）
*set,e5,500 ! 上法兰外径过渡圆角坐标的方向位移
*set,e6,510
*set,d11,7 ! 垫片直径
*set,rr6,110
*afun,deg ! 角度单位设定
*set,al1,30
*set,al2,2
*set,al3,75
*set,pi,0.3 ! 设计压力（内压）
*set,FORC_BOLT,86260.4
*set,py,2.60 ! 端面载荷
/prep7
et,1,solid45 ! 实体单元
et,2,HYPER58
!定义螺母材料1
mp,ex,1,1.92e5 ! 设定弹性模量
mp,nuxy,1,0.27 ! 定义泊松比
mp,dens,1,7.85e-9 ! 材料的密度
!定义双头螺柱材料2
mp,ex,2,2.09e5 ! 设定弹性模量
mp,nuxy,2,0.269 ! 定义泊松比
mp,dens,2,7.85e-9 ! 材料的密度
!定义平盖、筒体材料3
mp,ex,3,2.0E5 ! 设定弹性模量
mp,nuxy,3,0.25 ! 定义泊松比
mp,dens,3,7.85e-9 ! 材料的密度
!定义法兰材料4
mp,ex,4,2.06E5 ! 设定弹性模量
mp,nuxy,4,0.3 ! 定义泊松比
mp,dens,4,7.85e-9 ! 材料的密度
!定义O形密封圈材料5
mp,ex,5, ! 设定弹性模量
mp,nuxy,5,0.4997 ! 定义泊松比
TB,MOONEY,5,1
TBDATA,1,0.5516
TBDATA,2,0.1739
!定义摩擦系数为材料6
mp,mu,6,0.3
!定义摩擦系数为材料7
mp,mu,7,0.25
! 创建模型
!创建平盖模型
k,1,0,0,0
k,2,0,0,-h6
k,3,e5,0,-h4
k,4,e5+a2,0,-h5
k,5,e5+a2+h5*tan(al1),0,0
k,6,d4/2,0,0
k,7,0,0,-d4/2
l,1,2
l,2,3
l,3,4
l,4,5
larc,6,7,1,d4/2
lptn,4,5
ldele,7,,,1
ldele,9,,,1
ldele,11,,,1
l,1,6
lcomb,6,10
al,1,2,3,6,8,4
k,5,0,0,-h3
k,7,e3/2,0,-h3
l,2,5
l,5,7
k,8,e3/2,0,-e4/2
l,1,9
cskp,11,0,9,1,4
lgen,2,6,,,-(d7/2-d4/2),,,,,0
lgen,2,6,,,-(d7/2-d4/2+c2),,,,,0
lgen,2,6,,,-(b1/2-d4/2),,,,,0
csdele,all
k,16,d7/2,0,0
larc,11,16,1,d7/2
l,6,16
k,17,d7/2+c2,0,0
larc,13,17,1,d7/2+c2
l,16,17
k,18,b1/2,0,0
larc,15,18,1,b1/2
l,17,18
ldele,9
k,19,e3/2+rr2,0,-e4/2
larc,7,19,8,rr2
lovlap,9,10
ldele,20,,,1
ldele,22,,,1
ldele,23,,,1
lcomb,19,24
k,10,0,0,-h3-c2
k,19,e3/2,0,-h3-c2
l,5,10
l,10,19
k,20,e3/2+rr3,0,-e4/2
larc,19,20,8,rr3
lovlap,11,20
ldele,22,,,1
ldele,25,,,1
ldele,27,,,1
lcomb,24,26
k,12,0,0,-b2/2
k,20,b1/2,0,-b2/2+b1/2*tan(al2)
l,10,12
l,12,20
lextnd,12,14,100
lovlap,12,20
ldele,25,,,1
ldele,27,,,1
lfillt,22,26,rr3
al,2,5,7,19,21,13,14,8,6,3
al,9,10,24,23,15,16,13,21,19,7
al,11,22,12,26,17,18,15,23,24,10
asel,s,,,1,4,1
vext,all,,,,s1
wpoff,0,t2
wprota,0,90
vsel,s,,,2,4,1
vsbw,all
vdele,5,,,1
vdele,6,,,1
vdele,7,,,1
wpoff,0,0,-c1
vsbw,9
vdele,3,,,1
wpstyl,defa
wprota,0,-90
allsel
cyl4,0,0,d1/2,,,,s1
vsbv,1,3
wpoff,0,0,s1
cyl4,0,0,d2/2,,,,-s1
vovlap,1,4
vdele,5,,,1
wpoff,0,0,-t1
vsbw,3
vdele,4,,,1
wpstyl,defa
!创建上法兰
lgen,2,108,111,3,,-t2,,,,0
k,71,0,0,0
cskp,11,0,23,71,70
lgen,2,68,,,(b3/2-d8/2),,,,,0
csdele,all
k,74,d8/2,0,0
larc,73,74,71,d8/2
l,74,65
k,75,0,0,-h9
k,76,0,0,-b4/2
l,75,76
k,77,b3/2,0,-b4/2+b3/2*tan(al2)
l,76,77
lextnd,68,70,200
lovlap,68,94
ldele,104,,,1
ldele,107,,,1
lfillt,106,95,rr3
k,78,e6,0,-h7
k,79,e6+a3,0,-h8
l,75,78
l,78,79
l,72,79
lcomb,69,107
al,94,104,69,79,83,66,95,68,106,91
vext,22,,,,-s2
NUMMRG,ALL
NUMCMP,ALL
kdele,71
wpoff,0,-s2
k,71,0,-s2,0
cskp,11,0,85,71,86
lgen,2,122,,,(b3/2-a5/2),,,,,0
csdele,all
lgen,2,124,,,,,b4/2-a7/2,,,0
lextnd,141,90,200
lextnd,142,92,200
lovlap,141,142
ldele,144,,,1
ldele,146,,,1
lfillt,143,145,rr6
k,94,a5/2,-s2,0
larc,89,94,71,a5/2
csys,4
k,95,0,0,-b4/2
k,96,554.9057,,-700.6223
k,97,721.1309,,-600.7441
k,98,808.1206,,-450.0735
k,99,b3/2,0,0
l,91,95
l,95,96
l,97,98
larc,98,99,71,b3/2
lextnd,146,96,200
lextnd,147,97,200
lovlap,146,147
ldele,150,,,1
ldele,151,,,1
lfillt,149,152,rr3
l,94,99
al,145,141,143,142,147,148,152,146,149,144
vext,64,,,,-(h1-s2)
kwpave,105
cyl4,0,0,rr5
vdrag,76,,,,,,154,153,151,150
kwpave,101
vext,96,,,-100
vovlap,7,8,9,10,11,12
vdele,13,,,1
vdele,14,,,1
vdele,15,,,1
vdele,16,,,1
vdele,17,,,1
vdele,21,,,1
vdele,20,,,1
vdele,19,,,1
vdele,18,,,1
vdele,22,,,1
wpstyl,defa
vptn,3,23
!创建筒体
vext,135,,,,-h2
!创建密封圈
wpoff,d7/2+c2/2,0,0
wpoff,0,d11/2,0
cyl4,0,0,d11/2
vdrag,83,,,,,,9,10,12,14
wpstyl,defa
vext,103,,,-50
vadd,12,13
wprota,0,0,90
vsbw,14
vdele,12,,,1
wpstyl,defa
!创建螺栓孔
wprota,0,90
cyl4,38,-675,d5/2,,,,-86
vgen,7,12,,,76,,,,,0
cyl4,570,-670,d5/2,,,,-86
cyl4,641,-645,d5/2,,,,-86
cyl4,696,-594,d5/2,,,,-86
cyl4,734,-529,d5/2,,,,-86
cyl4,770,-464,d5/2,,,,-86
cyl4,806,-399,d5/2,,,,-86
cyl4,837,-330,d5/2,,,,-86
cyl4,861,-259,d5/2,,,,-86
cyl4,880,-187,d5/2,,,,-86
cyl4,892,-113,d5/2,,,,-86
cyl4,899,-38,d5/2,,,,-86
vsel,s,,,12
vsel,a,,,14,30,1
cm,bolt,volu
allsel
vsbv,6,bolt
!创建螺栓
cyl4,38,-675,d6/2,,,,-86
vgen,7,6,,,76,,,,,0
cyl4,570,-670,d6/2,,,,-86
cyl4,641,-645,d6/2,,,,-86
cyl4,696,-594,d6/2,,,,-86
cyl4,734,-529,d6/2,,,,-86
cyl4,770,-464,d6/2,,,,-86
cyl4,806,-399,d6/2,,,,-86
cyl4,837,-330,d6/2,,,,-86
cyl4,861,-259,d6/2,,,,-86
cyl4,880,-187,d6/2,,,,-86
cyl4,892,-113,d6/2,,,,-86
cyl4,899,-38,d6/2,,,,-86
!创建螺母
wpoff,0,0,-s1
rpr4,6,38,-675,nut_rad,,-nut_h
vgen,7,30,,,76,,,,,0
rpr4,6,570,-670,nut_rad,,-nut_h
rpr4,6,641,-645,nut_rad,,-nut_h
rpr4,6,696,-594,nut_rad,,-nut_h
rpr4,6,734,-529,nut_rad,,-nut_h
rpr4,6,770,-464,nut_rad,,-nut_h
rpr4,6,806,-399,nut_rad,,-nut_h
rpr4,6,837,-330,nut_rad,,-nut_h
rpr4,6,861,-259,nut_rad,,-nut_h
rpr4,6,880,-187,nut_rad,,-nut_h
rpr4,6,892,-113,nut_rad,,-nut_h
rpr4,6,899,-38,nut_rad,,-nut_h
!视镜孔处法兰联接螺栓的当量面（垫片）
wpoff,0,0,t2
cyl4,,,d3/2-d6/2,,d3/2+d6/2,-90,s1-t2
vovlap,1,49
wpstyl,defa
allsel
NUMMRG,ALL
NUMCMP,ALL
!划分网格
type,1
vsel,s,,,5
vsel,a,,,11
vsel,a,,,13,28,1
vsel,a,,,29
vsel,a,,,31,47,1
vptn,all
mat,1
vsel,s,,,105,122,1
cm,nut,volu
esize,6
vsweep,all
! 螺母划分
mat,2
vsel,s,,,51,68,1
vsweep,all
esize,12
vsel,s,,,69,103,2
vsweep,all
esize,18
vsel,s,,,70,104,2
cm,nut_bot,volu
vsweep,all
allsel
mat,3
!kwpave,26
!wprota,0,0,90
!vsbw,3
!kwpave,1
!vsbw,5
!wpstyl,defa
wpoff,0,c1+2
wprota,0,90
vsbw,4
vsbw,30
wpstyl,defa
kwpave,655
wprota,0,90
vsel,s,,,3
!vsel,a,,,12,13,1
vsel,a,,,48,50,1
cm,flat_head,volu
vsbw,flat_head
wpstyl,defa
lesize,111,,,6
lesize,34,,,6
lesize,55,,,6
lesize,33,,,6
vsel,s,,,4,5,1
vsel,a,,,11,13,2
vsel,a,,,1
vsweep,all
vsel,s,,,20,21,1
vsweep,all
vsel,s,,,14,19,1
vsweep,all
allsel
mat,4
vext,14,,,,-100
vovlap,7,3
vdele,23,,,1
vdele,22,,,1
vext,75,,,,100
vovlap,26,3
vdele,7,,,1
lesize,131,,,8
lesize,132,,,8
lesize,1181,,,6
vsel,s,,,22,25,1
vsweep,all
VSEL,S,,,6
VPLOT
ASLV
LSLA
VDRAG,65,,,,,,158
VDRAG,41,,,,,,118
VDRAG,205,,,,,,119
VDRAG,210,,,,,,120
BTOL,0.1
VSBA,27,64
VDELE,28,,,1
VDELE,6,,,1
VGLUE,ALL
lesize,232,,,6
lesize,222,,,6
lesize,157,,,6
vsel,s,,,6
vsel,a,,,27,28,1
vsel,a,,,30
vsweep,all
mat,3
vsel,s,,,2
vsweep,all
type,2
mat,5
esize,1
vsel,s,,,8,10,1
vsel,a,,,12
cm,O_seal,volu
vsweep,all
type,1
! 创建螺栓预紧单元
ALLSEL
MAT,2
PSMESH,1,BOLT1,,VOLU,70,0,Y,s1/2,,,,BOLT_1
PSMESH,2,BOLT1,,VOLU,72,0,Y,s1/2,,,,BOLT_2
PSMESH,3,BOLT1,,VOLU,74,0,Y,s1/2,,,,BOLT_3
PSMESH,4,BOLT1,,VOLU,76,0,Y,s1/2,,,,BOLT_4
PSMESH,5,BOLT1,,VOLU,78,0,Y,s1/2,,,,BOLT_5
PSMESH,6,BOLT1,,VOLU,80,0,Y,s1/2,,,,BOLT_6
PSMESH,7,BOLT1,,VOLU,82,0,Y,s1/2,,,,BOLT_7
PSMESH,8,BOLT1,,VOLU,84,0,Y,s1/2,,,,BOLT_8
PSMESH,9,BOLT1,,VOLU,86,0,Y,s1/2,,,,BOLT_9
PSMESH,10,BOLT1,,VOLU,88,0,Y,s1/2,,,,BOLT_10
PSMESH,11,BOLT1,,VOLU,90,0,Y,s1/2,,,,BOLT_11
PSMESH,12,BOLT1,,VOLU,92,0,Y,s1/2,,,,BOLT_12
PSMESH,13,BOLT1,,VOLU,94,0,Y,s1/2,,,,BOLT_13
PSMESH,14,BOLT1,,VOLU,96,0,Y,s1/2,,,,BOLT_14
PSMESH,15,BOLT1,,VOLU,98,0,Y,s1/2,,,,BOLT_15
PSMESH,16,BOLT1,,VOLU,100,0,Y,s1/2,,,,BOLT_16
PSMESH,17,BOLT1,,VOLU,102,0,Y,s1/2,,,,BOLT_17
PSMESH,18,BOLT1,,VOLU,104,0,Y,s1/2,,,,BOLT_18
! 建立约束方程
! 螺栓和法兰的约束
ALLSEL
CMSEL,S,nut_bot
ASLV
ASEL,R,LOC,Y,0
cm,acont,area
vsel,s,,,25
aslv
Asel,r,loc,y,0
nsla
esln
cm,econ,elem
cmsel,s,acont
nsla
ceintf
! 上法兰两端面的约束
ALLSEL
vsel,s,,,22,25,1
ASLV
ASEL,R,LOC,Y,-s2
cm,acont1,area
vsel,s,,,6
vsel,a,,,27,28,1
vsel,a,,,30
aslv
Asel,r,loc,y,-s2
nsla
esln
cm,econ1,elem
cmsel,s,acont1
nsla
ceintf
! 法兰下端和筒体的约束
ALLSEL
vsel,s,,,6
vsel,a,,,27,28,1
vsel,a,,,30
ASLV
ASEL,R,LOC,Y,-h1
cm,acont2,area
vsel,s,,,2
aslv
Asel,r,loc,y,-h1
nsla
esln
cm,econ2,elem
cmsel,s,acont2
nsla
ceintf
! 建立接触
! 建立螺母和平盖的接触
allsel
CMSEL,S,nut ! 选择螺母体积组件
ASLV ! 选择依附于体上的面
ASEL,R,LOC,Y,s1 ! 通过坐标选择螺栓上的接触面
CM,AS_NUT,AREA ! 创建螺母与平盖的接触面组件
ALLSEL
VSEL,S,,,4 ! 选择与螺母接触的平盖体积
ASLV ! 选择依附于体上的面
ASEL,R,LOC,Y,s1 ! 通过坐标选择平盖上的目标面
CM,AS_FLAT1,AREA ! 创建与螺母接触的平盖目标面
et,4,targe170 ! 3D接触目标面单元170
et,5,conta174 ! 3D接触接触面单元174
KEYOPT,5,5,1
KEYOPT,5,9,0
KEYOPT,5,12,2
R,4
REAL,4
MAT,6
CMSEL,S,AS_FLAT1
TYPE,4 ! 激活目标单元类型
NSLA,S,1 ! 选择目标面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建目标单元
ALLSEL
CMSEL,S,AS_NUT
TYPE,5 ! 激活接触单元类型
NSLA,S,1 ! 选择接触面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建接触单元
allsel
! 建立o形圈和沟槽左侧的接触
ALLSEL
VSEL,S,,,12 ! 选择与密封圈接触的沟槽体积
ASLV ! 选择依附于体上的面
ASEL,R,,,15,16,1
ASEL,a,,,18,19,1 ! 选择沟槽的目标面
CM,AS_groove1,AREA ! 创建与密封圈接触的沟槽目标面
allsel
CMSEL,S,O_seal ! 选择密封圈体积组件
ASLV ! 选择依附于体上的面
ASEL,R,,,81,82,1
ASEL,a,,,86,87,1
ASEL,a,,,91,92,1
ASEL,a,,,96,97,1 ! 选择密封圈上的接触面
CM,AS_seal1,AREA ! 创建密封圈与沟槽的接触面组件
et,6,targe170 ! 3D接触目标面单元170
et,7,conta174 ! 3D接触接触面单元174
KEYOPT,7,5,1
KEYOPT,7,9,0
R,5
REAL,5
MAT,7
CMSEL,S,AS_groove1
TYPE,6 ! 激活目标单元类型
NSLA,S,1 ! 选择目标面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建目标单元
ALLSEL
CMSEL,S,AS_seal1
TYPE,7 ! 激活接触单元类型
NSLA,S,1 ! 选择接触面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建接触单元
allsel
! 建立o形圈和沟槽右侧的接触
CMSEL,S,O_seal ! 选择密封圈体积组件
ASLV ! 选择依附于体上的面
ASEL,R,,,80,83,3
ASEL,a,,,85,88,3
ASEL,a,,,90,93,3
ASEL,a,,,95,98,3 ! 选择密封圈上的接触面
CM,AS_seal2,AREA ! 创建密封圈与沟槽的接触面组件
ALLSEL
VSEL,S,,,13 ! 选择与密封圈接触的沟槽体积
ASLV ! 选择依附于体上的面
ASEL,R,,,21,23,1
ASEL,a,,,25 ! 选择沟槽的目标面
CM,AS_groove2,AREA ! 创建与密封圈接触的沟槽目标面
et,8,targe170 ! 3D接触目标面单元170
et,9,conta174 ! 3D接触接触面单元174
KEYOPT,9,5,1
KEYOPT,9,9,0
R,6
REAL,6
MAT,7
CMSEL,S,AS_groove2
TYPE,8 ! 激活目标单元类型
NSLA,S,1 ! 选择目标面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建目标单元
ALLSEL
CMSEL,S,AS_seal2
TYPE,9 ! 激活接触单元类型
NSLA,S,1 ! 选择接触面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建接触单元
allsel
! 建立o形圈和平盖的接触
CMSEL,S,O_seal ! 选择密封圈体积组件
ASLV ! 选择依附于体上的面
ASEL,R,,,80,81,1
ASEL,a,,,85,86,1
ASEL,a,,,90,91,1
ASEL,a,,,95,96,1 ! 选择密封圈上的接触面
CM,AS_seal3,AREA ! 创建密封圈与平盖的接触面组件
ALLSEL
VSEL,S,,,1 ! 选择与密封圈接触的平盖体积
ASLV ! 选择依附于体上的面
ASEL,R,loc,y,c1+t2 ! 选择平盖的目标面
CM,AS_groove3,AREA ! 创建与密封圈接触的平盖目标面
et,10,targe170 ! 3D接触目标面单元170
et,11,conta174 ! 3D接触接触面单元174
KEYOPT,11,5,1
KEYOPT,11,9,0
R,7
REAL,7
MAT,7
CMSEL,S,AS_groove3
TYPE,10 ! 激活目标单元类型
NSLA,S,1 ! 选择目标面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建目标单元
ALLSEL
CMSEL,S,AS_seal3
TYPE,11 ! 激活接触单元类型
NSLA,S,1 ! 选择接触面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建接触单元
allsel
! 建立o形圈和沟槽内表面的接触
CMSEL,S,O_seal ! 选择密封圈体积组件
ASLV ! 选择依附于体上的面
ASEL,R,,,82,83,1
ASEL,a,,,87,88,1
ASEL,a,,,92,93,1
ASEL,a,,,97,98,1 ! 选择密封圈上的接触面
CM,AS_seal4,AREA ! 创建密封圈与平盖的接触面组件
ALLSEL
VSEL,S,,,24 ! 选择与密封圈接触的平盖体积
ASLV ! 选择依附于体上的面
ASEL,R,loc,y,0 ! 选择平盖的目标面
CM,AS_groove4,AREA ! 创建与密封圈接触的平盖目标面
et,12,targe170 ! 3D接触目标面单元170
et,13,conta174 ! 3D接触接触面单元174
KEYOPT,13,5,1
KEYOPT,13,9,0
R,8
REAL,8
MAT,7
CMSEL,S,AS_groove4
TYPE,12 ! 激活目标单元类型
NSLA,S,1 ! 选择目标面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建目标单元
ALLSEL
CMSEL,S,AS_seal4
TYPE,13 ! 激活接触单元类型
NSLA,S,1 ! 选择接触面上所有的节点，包括边界节点
ESLN,S,0 ! 选择依附于这些节点的单元
ESURF,ALL ! 构建接触单元
allsel
! 求解设置
!NSEL,S,LOC,Y,-h1-h2
!NROTAT,ALL
!CSYS,1
!D,ALL,UX
/SOLU
! 施加预紧力
CUTCONTROL,PLSLIMIT,10,
CUTCONTROL,CRPLIMITexp,0.2,0
CUTCONTROL,CRPLIMITimp,0.25,1
LSCLEAR,ALL ! 清空所有载荷
ALLSEL
CSYS,0 ! 激活笛卡尔坐标系
ASEL,S,LOC,x,0 ! 通过坐标选择为0的面
DA,ALL,SYMM ! 施加对称边界条件
ASEL,S,LOC,z,0 ! 通过坐标选择为0的面
DA,ALL,SYMM ! 施加对称边界条件
nSEL,S,LOC,y,-h1-h2 ! 通过坐标选择面
D,ALL,uy ! 施加Y向位移约束
nsel,s,,,7028
D,ALL,uy
NROPT,FULL ! 激活完全NEWTON-RAPHSON
PRED,ON ! 打开变形预测
LNSRCH,on ! 激活线性搜索
nlgeom,on
ANTYPE,0 ! 静态分析
AUTOTS,1 ! 自动时间步长
NSUBST,10,20,1 ! 设置载荷步
OUTRES,ALL,ALL ! 输出所有结果
ALLSEL
SLOAD,1,9,LOCK,FORC,FORC_BOLT,1,2 ! 在螺栓中施加螺栓预紧
SLOAD,2,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,3,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,4,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,5,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,6,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,7,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,8,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,9,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,10,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,11,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,12,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,13,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,14,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,15,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,16,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,17,9,LOCK,FORC,FORC_BOLT,1,2
SLOAD,18,9,LOCK,FORC,FORC_BOLT,1,2
ALLSEL
TIME,10 ! 设置时间步至10
LSWRITE,1 ! 写出第一个载荷步文件
!对容器壁施加内压
lsel,s,radius,,d1/2
asll
asel,u,loc,y,0
asel,u,loc,y,c1+t2
asel,u,loc,y,s1-t1
asel,a,,,362,372,10
asel,a,,,15,16,1
asel,a,,,364,368,4
asel,a,,,43
ASEL,A,,,30,34,2
ASEL,A,,,29
asel,a,,,39
CMSEL,a,AS_seal1
CMSEL,a,AS_groove1
ASEL,A,,,62,63,1
ASEL,A,,,66
ASEL,A,,,100
ASEL,A,,,27,28,1
ASEL,A,,,33
ASEL,A,,,36,37,1
ASEL,A,,,46
ASEL,A,,,60
ASEL,A,,,75
! 选中施加内压的面
SFA,ALL,1,PRES,Pc ! 对内壁施加压强,大小为Pc
ASEL,S,,,369
SFA,ALL,1,PRES,-PY ! 施加端面当量面载荷
ALLSEL
TIME,20 ! 设置时间步至20
LSWRITE,2 ! 写出第二个载荷步文件
LSSOLVE,1,2,1 ! 计算载荷步1-2
