finish/clear!!!!!!!!!!!!!!!!环境设置!!!!!!!!!!!!!!!!!!!!!!!!!/filn,flange/title,flang_bolt!设定结构基本参数,单位N,mm,MPa.密度为tons/mm3*SET,PN,4!工作压力!结构尺寸参数遵循图形*SET,DN,900*SET,D,1115*SET,D1,1055*SET,D2,1010*SET,D3,1000*SET,D4,997*SET,dta,86*SET,H,170*SET,xh,42*SET,a,3*SET,a1,3*SET,dta1,24*SET,dta2,36*SET,R,15*SET,xd,33*SET,h1,225*SET,h2,50*SET,L,2.5*((DN/2+0.5*dta1)*dta1)**0.5!筒体长*SET,N,360/96!法兰分度数*SET,nut_rad,50.85/2!螺母半径*SET,nut_h,25.6!螺母高度*SET,bolt_h,250!螺栓高度*SET,forc_bolt,35000!螺栓预紧力/prep7!进入前处理器et,1,solid45!实体单元et,2,TARGE170!3D接触目标面单元170et,3,CONTA174!3D接触接触面单元174keyopt,3,9,0!包含初始化穿透*AFUN,DEG!设定角度单位为度数!定义法兰材料1mp,ex,1,2.16e5!设定弹性模量mp,nuxy,1,0.3!定义泊松比mp,dens,1,7.85e-9!定义密度!定义螺栓材料2mp,ex,2,2.11e5!设定弹性模量mp,nuxy,2,0.3!定义泊松比mp,dens,2,7.85e-9!定义密度!定义摩擦系数为材料3mp,mu,3,0.3!!!!!!!!!!!!!!创建上法兰!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!创建关键点k,1,DN/2,a+Hk,2,DN/2+dta1,a+Hk,3,DN/2+dta1,a+xh+dtak,4,DN/2+dta2,a+dtak,5,D1/2,a+dtak,6,D/2,a+dtak,7,D/2,ak,8,D1/2,ak,9,D4/2,ak,10,D4/2,0k,11,DN/2,0k,12,DN/2,a+dtak,13,DN/2,a+dta+xh!通过关键点连线*do,ii,1,12l,ii,ii+1*enddol,1,13lfillt,2,3,R!创建倒角lfillt,3,4,R!创建倒角l,13,14l,12,17l,17,9l,5,8!连线构面al,1,2,16,13al,16,14,3,15,17,12al,17,18,9,10,11al,4,19,8,18al,5,6,7,19cm,flang_up,area!创建上法兰面积组件!!!!!!!!!!!!!上法兰创建完毕!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!创建下法兰!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!创建关键点k,18,DN/2,a-a1k,19,D3/2,a-a1k,20,D3/2,ak,21,D1/2,ak,22,D/2,ak,23,D/2,a-a1-dtak,24,D1/2,a-a1-dtak,25,D3/2,a-a1-dtak,26,DN/2+dta2,a-a1-dtak,27,DN/2+dta1,a-a1-dta-xhk,28,DN/2+dta1,a-a1-Hk,29,DN/2,a-a1-Hk,30,DN/2,a-a1-dta-xhk,31,DN/2,a-a1-dta!通过关键点连线*do,ii,18,30l,ii,ii+1*enddol,18,31lfillt,27,28,Rlfillt,28,29,Rl,30,35l,31,32l,21,24l,19,25!连线构面al,29,30,31,36al,35,36,32,37,34,28al,20,39,27,37,33al,21,22,38,26,39al,23,24,25,38asel,u,area,,flang_up!选择除flang_up之外的所有面cm,flang_down,area!创建下法兰面积组件!!!!!!!!!!!!!下法兰创建完毕!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!创建椭圆封头!!!!!!!!!!!!!!!!!!!!!!!!!!!!!k,,0,a+H+h2!创建封头所处位置关键点circle,36,1!创建单位长圆KWPAVE,36!将工作平面移置于新设定点CSYS,4!激活随工作平面的坐标系lsscale,40,43,1,DN/2+dta1,h1+dta1,,,1!对单位圆按照椭圆比例进行缩放(外表面)lsscale,40,43,1,DN/2,h1,,,1!对单位圆按照椭圆比例进行缩放(内表面)ldele,49,51,1,1!删除多余的线段ldele,45,47,1,1ldele,40,43,1,1wpstyl,defa!恢复工作平面至默认状态l,45,41!连接封头端、尾l,42,46l,1,45!连接封头与法兰l,2,41!由线构成封头面al,40,42,1,43al,44,40,48,41asel,u,area,,flang_down!选择构成封头的所有面cm,coping,area!创建封头面积组件!!!!!!!!!!!!!封头创建完毕!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!创建筒体!!!!!!!!!!!!!!!!!!!!!!!!!!!!!allselk,,DN/2,a-a1-H-L!创建筒体底端位置点k,,DN/2+dta1,a-a1-H-La,28,29,37,38!连接筒体与下法兰,构成面asel,u,area,,flang_upasel,u,area,,flang_downasel,u,area,,coping!排除筒体之外的所有面cm,pot,area!创建筒体面积组件!!!!!!!!!!!!!筒体创建完毕!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!布尔操作,将面与面之间粘接!!!!!!!!!!!!!asel,s,area,,flang_upasel,a,area,,copingaglue,all!将封头与上法兰面粘接asel,s,area,,flang_downasel,a,area,,potaglue,all!将筒体与下法兰面粘接kdele,36!删除建模过程中辅助点!!!!!!!!!!!!!布尔操作完毕!!!!!!!!!!!!!!!!!!!!!!!!!!allselvrotat,all,,,,,,42,46,N!将所有的面绕y轴旋转N度,创建体!!!!!!!!!!!!!创建螺栓!!!!!!!!!!!!!KL,76,0.5,,!设定关键点平分线76,该点为螺栓中心所在KWPAVE,77!将工作平面移动至该点wprota,0,90!将工作平面绕x轴旋转90度cyl4,,,xd/2,,,,bolt_h!创建螺栓!将与螺栓相交的体组成一个组件,以便选取vsel,s,volu,,4,5,1vsel,a,volu,,9,10,1cm,buer,volu!创建与螺栓相交的体积组件buerallselbopt,numb,offvsbv,buer,14!将组件buer与螺栓进行布尔操作,挖出螺栓孔洞cyl4,,,xd/2,,,,bolt_h!创建螺栓,体积编号为4VGEN,,4,,,,,-(bolt_h/2-dta-a),,,1!将螺栓平移至合适位置rpr4,6,,,nut_rad,,-nut_h!创建螺母vgen,2,5,,,,,2*dta+a1+nut_h,,,0!复制螺母至合适位置!!!!!!!!!!!!!螺栓创建完毕!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!创建各部分体积组件方便选择!!!!!!!!!!!!!vsel,s,volu,,4,5,1vsel,a,volu,,9vptn,allcm,bvolu,volu!创建螺栓组件allselwpstyl,defavsel,s,loc,y,-H-L,a1vsel,u,volu,,bvolucm,flang_dvolu,volu!创建下法兰组件allselvsel,s,loc,y,a1,Hvsel,u,volu,,bvolucm,flang_upvolu,volu!创建上法兰组件allselvsel,s,loc,y,H+a,H+h1+dta1cm,cop_volu,volu!创建封头组件!!!!!!!!!!!!!!!划分网格!!!!!!!!!!!!!!!!!!!!!!!!!!!!mat,1!指定材料esize,5!指定单元大小cmsel,s,cop_volu!选中封头组件vsweep,all!对封头分网esize,5!指定单元大小cmsel,s,flang_dvolu!选中下法兰组件vsweep,all!对下法兰分网cmsel,s,flang_upvolu!选中上法兰组件vsweep,all!对上法兰分网esize,4!指定单元大小mat,2!指定材料cmsel,s,bvolu!选中螺栓组件vsweep,all!对螺栓分网!!!!!!!!!!!!!网格划分完毕!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!对模型节点进行压缩合并!!!!!!!!!!!!!!!!!!!!!!合并上法兰节点allselvsel,s,loc,y,0,2*Hvsel,u,volu,,bvolunslv,s,1nummrg,node!合并下法兰节点allselvsel,s,loc,y,-(H+L),a1vsel,u,volu,,bvolunslv,s,1nummrg,node!!!!!!!!!!!!!节点合并完毕!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建立螺栓预紧单元!!!!!!!!!!!!!!!!!!!!!allsel!在0号坐标系下,对螺栓1(14号体积)建立预紧单元,secid为1,name为bolt1!预紧单元方向为y轴,位置为y向0位置处,同时将预紧单元存储于bolt_1单元组件中,预紧单元pretention179无需预先定义psmesh,1,bolt1,,volu,14,0,y,0,,,,bolt_1!!!!!!!!!!!!!螺栓预紧单元建立完毕!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!创建接触!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!创建凹面接触单元cmsel,s,flang_upvolu!选择上法兰体积组件aslv!选择依附于体上的面asel,r,loc,y,0!通过坐标选择法兰上的接触面cm,as_up,area!建立上法兰接触面组件as_upallselcmsel,s,flang_dvolu!选择下法兰体积组件aslv!选择依附于体上的面asel,r,loc,y,a-a1!通过坐标选择法兰上的接触面asel,r,loc,x,0,D3/2cm,as_down,area!建立下法兰接触面组件as_downr,1!定义接触实常数1real,1!指定实常数1mat,3!指定材料3cmsel,s,as_up!选择上法兰接触面组件as_uptype,2!指定单元接触类型2nsla,s,1!选择接触面上的所有节点,包括边界点esln,s,0!选择依附于这些节点的单元esurf,all!构建接触单元allselcmsel,s,as_down!选择下法兰接触面组件as_downtype,3!指定单元接触类型3nsla,s,1!选择接触面上的所有节点,包括边界点esln,s,0!选择依附于这些节点的单元esurf,all!构建接触单元!创建凸面接触单元allselcmsel,s,flang_upvolu!选择上法兰体积