24小时热门版块排行榜

返回列表

傻帽1004

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 25.1
帖子: 38
在线: 15.8小时
虫号: 3667366
注册: 2015-01-28
性别: GG
专业: 固体力学

[求助] 螺栓法兰平垫片连接

请问我用六个螺栓模拟管道法兰密封，在加载预紧力情况下螺栓的应力各不相同，况且弯曲应力较大与薄膜应力差不多  正常弯曲应力不会这么大！这是咋回事！请各位大神不吝赐教！！！命令流如下：
finish
/clear
/filn,flange
/title,flang_bolt
*afun,deg
n=360/2
d1=76
d2=107.95
d3=146.05
d4=177.5
d5=111
d6=84.15
d7=87.4
d8=104.8
h1=7
h2=54.7
h3=25
h4=117
h5=4.5
tong_h=63
bd1=30
bd2=27
bolt_h=190
nut_rad=20.5
nut_h=27
r=10
forc_bolt=26243
pc=4.5
pa=20
/prep7
et,1,solid45
et,2,targe170
et,3,conta174
keyopt,3,9,0
et,4,inter195
!定义法兰筒体材料1!
MP,EX,1,2.13E5
MP,NUXY,1,0.286
!定义螺栓材料2!
MP,EX,2,2.14E5
MP,NUXY,2,0.28
!定义摩擦系数为材料3!
MP,MU,3,0.3
!定义垫片材料4!
DALTA0=0.00e-3
STIFF0=0.0e7
SCAP=1.0e-5
TB,GASKET,4,,,PARA
TBDATA,1,DELTA0,STIFF0,SCAP
!定义垫片加载曲线1!
tb,gasket,4,1,6,comp
TBTEMP,20.000
tbpt,,0,0
tbpt,,0.2,10
tbpt,,0.38,20
tbpt,,0.5,30
tbpt,,0.82,60
tbpt,,1.04,80
!定义垫片实际卸载曲线1!
tb,gasket,4,1,8,nunl
TBTEMP,20
tbpt,,1.04,80
tbpt,,1.035,70
tbpt,,1.03,60
tbpt,,1.01,40
tbpt,,1.00,30
tbpt,,0.98,20
tbpt,,0.92,10
tbpt,,0.80,0
!定义螺母材料10!
MP,EX,10,2.15E5
MP,NUXY,10,0.286
!上法兰!
k,1,d1,0
k,2,d2,0
k,3,d2,h1
k,4,d4,h1
k,5,d4,h2
k,6,d5,h2
k,7,d6,h4-h3
k,8,d6,h4
k,9,d1,h4
*do,i,1,8
l,i,i+1
*enddo
l,9,1
lfillt,6,7,r
lfillt,6,5,r
l,3,13
al,1,2,12,11,6,10,7,8,9
al,12,5,4,3
kwpave,12
wprota,,90
asbw,all
kwpave,10
asbw,all
kwpave,13
asbw,3
wpstyl,defa
k,18,d1,tong_h+ky(9)
k,19,d6,tong_h+ky(9)
a,18,19,8,9
NUMCMP,AREA
agen,2,1,6,1,,h5/2,,,,1
k,20,d7,h5/2
k,21,d8,h5/2
kwpave,20
csys,4
wprota,,,90
asbw,6
kwpave,21
asbw,8
wpstyl,defa
allsel
arsym,y,all
!创建垫片!
k,45,d7,-h5/2
k,46,d8,-h5/2
a,20,21,46,45
aglue,all
!创建回转体!
k,1001,0,0
k,1002,0,100
vrotat,all,,,,,,1001,1002,n
kdele,1001,1002,1
!!!创建螺栓孔1!!!!
wprota,,,15
k,1000,d3,ky(5)
kwpave,1000
wprota,0,90
csys,4
cyl4,,,bd1/2,,,,bolt_h
!复制螺栓孔！
wpstyl,defa
wprota,,,15
csys,5
vgen,6,35,,,,30,,,,0
vsel,s,volu,,2
vsel,a,volu,,10
cm,buer1,volu
allsel
BOPT,NUMB,OFF
vsbv,buer1,35
vsel,s,volu,,41,42,1
cm,buer2,volu
allsel
vsbv,buer2,36
vsel,s,volu,,2
vsel,a,volu,,10
cm,buer3,volu
allsel
vsbv,buer3,37
vsel,s,volu,,19
vsel,a,volu,,27
cm,buer4,volu
allsel
vsbv,buer4,38
vsel,s,volu,,2
vsel,a,volu,,10
cm,buer5,volu
allsel
vsbv,buer5,39
vsel,s,volu,,27
vsel,a,volu,,19
cm,buer6,volu
allsel
vsbv,buer6,40
wpstyl,defa
!创建螺栓!!
kwpave,1000
wprota,0,90
csys,4
cyl4,,,bd2/2,,,,bolt_h
wpstyl,defa
vgen,,19,,,,-(ky(187)-ky(31))/2,,,,1
!螺母!
kwpave,1000
wprota,0,90
csys,4
rpr4,6,,,nut_rad,,-nut_h
wpstyl,defa
vgen,2,27,,,,-ky(5)+ky(31)-nut_h,,,,0
!复制螺栓！
allsel
wpstyl,defa
wprota,,,15
csys,5
vgen,6,19,27,8,,30,,,,0
vgen,6,37,,,,30,,,,0
!分割螺栓!
vsel,s,volu,,19
vsel,a,volu,,27
vsel,a,volu,,37
vptn,all
cm,bvolu1,volu
vsel,s,volu,,38,39,1
vsel,a,volu,,48
vptn,all
cm,bvolu2,volu
allsel
vsel,s,volu,,40,41,1
vsel,a,volu,,49
vptn,all
cm,bvolu3,volu
allsel
vsel,s,volu,,42,43,1
vsel,a,volu,,50
vptn,all
cm,bvolu4,volu
allsel
vsel,s,volu,,44,45,1
vsel,a,volu,,51
vptn,all
cm,bvolu5,volu
allsel
vsel,s,volu,,46,47,1
vsel,a,volu,,52
vptn,all
cm,bvolu6,volu
allsel
wpstyl,defa
!划分网格！
!下法兰与垫片接触处!
mat,1
type,1
lesize,144,,,35
lesize,49,,,10
lesize,43,,,4
vmesh,14
mat,4
type,4
lesize,59,,,1
vmesh,17
mat,1
type,1
lesiz,24,,,10
vmesh,6
lesize,240,,,35
vmesh,31
mat,4
type,4
vmesh,34
mat,1
type,1
vmesh,23
!上法兰!
!倒角处!
vsel,s,volu,,7
vsel,a,volu,,9
vsel,a,volu,,4
vsel,a,volu,,1
esize,7
vsweep,all
vsel,s,volu,,24
vsel,a,volu,,26
vsel,a,volu,,21
vsel,a,volu,,18
esize,7
vsweep,all
vsel,s,volu,,12
vsel,a,volu,,15
vsel,a,volu,,8
vsel,a,volu,,16
esize,7
vsweep,all
vsel,s,volu,,25
vsel,a,volu,,29
vsel,a,volu,,32
vsel,a,volu,,33
esize,7
vsweep,all
vsel,s,volu,,10
vsel,a,volu,,2
vsel,a,volu,,35
vsel,a,volu,,36
lesize,5,,,8
lesize,35,,,8
lesize,162,,,8
lesize,208,,,8
vsweep,all
vsel,s,volu,,5
vsel,a,volu,,22
esize,10
vsweep,all
vsel,s,volu,,3
vsel,a,volu,,20
esize,20
vsweep,all
vsel,s,volu,,13
vsel,a,volu,,30
esize,10
vsweep,all
vsel,s,volu,,28
vsel,a,volu,,11
esize,20
vsweep,all
!螺栓!
cmsel,s,bvolu1,volu
vsel,u,,,58
vsel,u,,,59
mat,2
type,1
esize,3
vsweep,all
cmsel,s,bvolu2,volu
vsel,u,,,62
vsel,u,,,63
mat,2
type,1
esize,3
vsweep,all
cmsel,s,bvolu3,volu
vsel,u,,,66
vsel,u,,,67
mat,2
type,1
esize,3
vsweep,all
cmsel,s,bvolu4,volu
vsel,u,,,70
vsel,u,,,71
mat,2
type,1
esize,3
vsweep,all
cmsel,s,bvolu5,volu
vsel,u,,,74
vsel,u,,,75
mat,2
type,1
esize,3
vsweep,all
cmsel,s,bvolu6,volu
vsel,u,,,78
vsel,u,,,79
mat,2
type,1
esize,3
vsweep,all
vsel,S,,,58,78,4
vsel,A,,,59,79,4
mat,10
type,1
esize,3
vsweep,all
VSEL,S,VOLU,,17,34,17
eslv
emodif,all,type,4
!创建上法兰接触!
!上法兰接触!
allsel
asel,s,area,,382
asel,a,area,,370
asel,a,area,,358
asel,a,area,,346
asel,a,area,,334
asel,a,area,,322
cm,b_up,area
allsel
asel,s,area,,183
asel,a,area,,195
cm,fnag_up,area
allsel
R,2
REAL,2
CMSEL,S,b_up
TYPE,2
NSLA,S,1
ESLN,S,0
ESURF,ALL
allsel
cmsel,s,fnag_up
TYPE,3
NSLA,S,1
ESLN,S,0
ESURF,ALL
ALLSEL
!下法兰接触!
asel,s,area,,325
asel,a,area,,337
asel,a,area,,349
asel,a,area,,361
asel,a,area,,373
asel,a,area,,385
cm,b_down,area
allsel
asel,s,area,,175
asel,a,area,,170
cm,flang_down,area
allsel
ET,7,TARGE170
ET,8,CONTA174
KEYOPT,8,9,0
R,3
REAL,3
MAT,3
CMSEL,S,b_down
TYPE,7
NSLA,S,1
ESLN,S,0
ESURF,ALL
ALLSEL
CMSEL,S,flang_down
TYPE,8
NSLA,S,1
ESLN,S,0
ESURF,ALL
allsel
NUMMRG,NODE
!预紧单元！
allsel
psmesh,1,bolt1,,volu,27,0,y,0,,,,bolt_1
psmesh,2,bolt2,,volu,54,0,y,0,,,,bolt_2
psmesh,3,bolt3,,volu,39,0,y,0,,,,bolt_3
psmesh,4,bolt4,,volu,41,0,y,0,,,,bolt_4
psmesh,5,bolt5,,volu,43,0,y,0,,,,bolt_5
psmesh,6,bolt6,,volu,45,0,y,0,,,,bolt_6
allsel
!求解!
finish
/sol
allsel
LSCLEAR,ALL
NROPT,full
PRED,ON
LNSRCH,ON
AUTOTS,1
NSUBST,10,100,1
OUTRES,ALL,ALL
EQSLV,PCG,1E-8
NROP,UNSYM
WPSTYL,DEFA
csys,0
ASEL,S,LOC,Y,ky(34)
DA,ALL,UY                                  ! 约束轴向位移
NSLA,,1                                  ! 按面选节点
NSEL,R,LOC,X,0                            ! 再选择X为0的节点
D,ALL,UX                                  ! 约束X方向位移
NSLA,,1                                  ! 按面选节点
NSEL,R,LOC,Z,0                            ! 再选择Z为0的节点
D,ALL,UZ                                  ! 约束Z方向位移
ASEL,S,LOC,Z,0                            ! 选择对称面
DA,ALL,SYMM                               ! 施加对称约束
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
time,10
allsel
lswrite,1
time,20
lswrite,2
!上法兰内表面！
allsel
asel,s,area,,104
asel,a,area,,31
asel,a,area,,41
asel,a,area,,114
asel,a,area,,20
asel,a,area,,93
asel,a,area,,107
asel,a,area,,34
asel,a,area,,120,121,1
asel,a,area,,47,48,1
cm,fup_area,area
!垫片!
allsel
asel,s,area,,162
asel,a,area,,89
cm,gas_area,area
!下法兰内表面!
allsel
asel,s,area,,155,156,1
asel,a,area,,82,83,1
asel,a,area,,71
asel,a,area,,144
asel,a,area,,51
asel,a,area,,124
asel,a,area,,147
asel,a,area,,74
asel,a,area,,136
asel,a,area,,63
cm,fdo_area,area
!内压!
allsel
cmsel,s,fup_area
cmsel,a,gas_area
cmsel,a,fdo_area
sfa,all,1,pres,pc
allsel
!等效载荷!
asel,s,loc,y,ky(18)
sfa,all,1,pres,-pa
allsel
time,30
lswrite,3
lssolve,1,3,1

回复此楼