24小时热门版块排行榜    

查看: 750  |  回复: 2

傻帽1004

铁虫 (初入文坛)

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

请问我用六个螺栓模拟管道法兰密封,在加载预紧力情况下螺栓的应力各不相同,况且弯曲应力较大与薄膜应力差不多  正常弯曲应力不会这么大!这是咋回事!请各位大神不吝赐教!!!命令流如下:
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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatewu

主管区长 (文坛精英)

小木虫浪漫体验师~\(^o^)/~

优秀区长优秀区长优秀区长优秀版主优秀版主优秀版主优秀版主优秀版主优秀版主

★ ★ ★ ★ ★
傻帽1004(Lancashire代发): 金币+5, 感谢交流。 2015-10-21 16:27:45
建议最好能够上传附件,这样太长,太浪费阅读时间了
【化学化工区欢迎您,小晨欢迎您~\(^o^)/~】何以飘零去,何以少团栾,何以别离久,何以不得安~\(^o^)/~走自己的路,相信一切皆有可能!
2楼2015-08-31 20:22:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

傻帽1004

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by fatewu at 2015-08-31 20:22:24
建议最好能够上传附件,这样太长,太浪费阅读时间了

已经上传,谢谢!!!

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : YU+NEI.txt
  • 2015-08-31 20:30:17, 8.08 K
3楼2015-08-31 20:30:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 傻帽1004 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
请填处理意见