| 查看: 789 | 回复: 2 | |||
[求助]
螺栓法兰平垫片连接
|
|
请问我用六个螺栓模拟管道法兰密封，在加载预紧力情况下螺栓的应力各不相同，况且弯曲应力较大与薄膜应力差不多 正常弯曲应力不会这么大！这是咋回事！请各位大神不吝赐教！！！命令流如下： 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 |
回复此楼» 猜你喜欢
极限拉扯，如何抉择
已经有14人回复
-
申请科研项目纯靠实力能成功吗？
已经有6人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有12人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
fatewu
主管区长 (文坛精英)
小木虫浪漫体验师~\(^o^)/~
-
专家经验: +790 - CE-EPI: 7
- 应助: 521 (博士)
- 贵宾: 19.629
- 金币: 19729
- 散金: 169099
- 红花: 648
- 沙发: 2560
- 帖子: 43961
- 在线: 6049.8小时
- 虫号: 3063713
- 注册: 2014-03-17
- 专业: 未来学
- 管辖: 化学化工区
2楼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
200