24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 680  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +13 努力奋斗112 2026-04-07 13/650 2026-04-07 15:04 by -迷了路啊路
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 304求调剂 +7 素年祭语 2026-04-06 10/500 2026-04-06 23:31 by 素年祭语
[考研] 工科 22408 267求推荐 +4 wanwan00 2026-04-05 5/250 2026-04-06 22:47 by chenzhimin
[考研] 一志愿河北工业大学材料工程,初试344求专硕调剂 +6 15933906766 2026-04-05 6/300 2026-04-06 13:21 by 无际的草原
[考研] 277求调剂 +5 考研调剂lxh 2026-04-05 5/250 2026-04-05 19:03 by chy09050039
[考研] 材料0856 英一数二 323 求调剂 +14 袁sy 2026-04-01 14/700 2026-04-05 18:18 by cql1109
[考研] 0832食品科学与工程学硕282调剂 +6 鱼在水中游a 2026-04-02 9/450 2026-04-05 11:45 by flysky1234
[考研] 男生,一志愿沪9生物学071000,初试308求调剂 +3 刘墨墨 2026-04-04 3/150 2026-04-05 08:26 by barlinike
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 调剂 +11 JLLLLLLLLLL 2026-04-03 11/550 2026-04-04 22:21 by hemengdong
[考研] 309求调剂 +6 刘刘刘1231 2026-04-02 7/350 2026-04-04 13:41 by liucky
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 求调剂不挑专业 +3 xrh030412 2026-04-01 3/150 2026-04-03 14:40 by 氮气气气
[考研] 315求调剂 +6 顺理成张 2026-04-03 8/400 2026-04-03 14:04 by 百灵童888
[考研] 一志愿陕西师范大学生物学317分 +5 1563日。 2026-04-02 5/250 2026-04-03 06:58 by ilovexiaobin
[考研] 282求调剂 +13 呼吸都是减肥 2026-04-01 13/650 2026-04-02 14:10 by baoball
[考研] 材料求调剂 +10 呢呢妮妮 2026-04-01 13/650 2026-04-02 09:17 by olim
[考研] 279求调剂 +6 学而思兮知 2026-04-01 6/300 2026-04-02 09:16 by vgtyfty
信息提示
请填处理意见