壓力容器接管處施加彎矩的處理和內(nèi)部加內(nèi)壓.doc

! * 設(shè)置環(huán)境變量 *FINISH/CLEAR/FILN, E65/TITLE, /UNITS,SI !采用國際單位制!*參數(shù)定義*DI=2000 ! 設(shè)備內(nèi)徑TC=16 ! 筒體及封頭壁厚HC=5*NINT(SQRT(DI/2*TC)/10)*10 ! 模型中筒體長(zhǎng)度（含封頭直邊高度)DN=630 ! 接管外徑TN=8 ! 接管壁厚LN=350 ! 接管外伸高度 NT=2 ! 厚度方向剖分?jǐn)?shù)NME=30 ! 橢圓封頭經(jīng)向剖分?jǐn)?shù)NH=20 ! 接管環(huán)向剖分?jǐn)?shù)M=1E6 ! 接管端部彎矩PC=0.2 ! 最高壓力PAX1=-PC*DN*2/(DN+2*TN)*2-DN*2) ! 接管端部軸向平衡面載荷! * 前處理 */PREP7 ET,1,95 ! 定義單元類型ET,2,21 ! 定義單元類型MP,EX,1,2E5 ! 定義材料的彈性模量MP,NUXY,1,0.3 ! 定義材料的泊松比R,1,1 ! 定義質(zhì)量單元的實(shí)常數(shù)! * 建立模型 *CSWPLA,11,1,0.5,1, ! 定義橢圓坐標(biāo)系K,1,-DI/2, ! 定義橢圓的一個(gè)端點(diǎn)K,2,DI/2, ! 定義橢圓的另一個(gè)端點(diǎn)L,1,2 ! 生成橢圓線CSYS,0 ! 激活總體直角坐標(biāo)系K,3,KX(2)+TC, ! 定義封頭壁厚關(guān)鍵點(diǎn)L,2,3 ! 生成封頭壁厚線ADRAG,2, , , , , ,1 ! 由封頭壁厚沿橢圓線生成橢圓面BLC4,DN/2-TN, ,TN,DI/4+TC+LN, ! 生成中心接管截面 APTN,ALL ! 面域互分ADELE,4,5,1 ! 刪除多余面*GET,KPMAX,KP,NUM,MAX ! 提取最大節(jié)點(diǎn)號(hào)K,KPMAX+1,0,0,0 ! 生成中心線一個(gè)關(guān)鍵點(diǎn)K,KPMAX+2,0,HC,0 ! 生成中心線另一個(gè)關(guān)鍵點(diǎn)VROTAT,ALL,KPMAX+1,KPMAX+2,180 ! 旋轉(zhuǎn)生成橢圓封頭及中心接管ASEL,S,LOC,Y,0 ! 選擇橢圓封頭端面*GET,KPMAX,KP,NUM,MAX ! 提取最大節(jié)點(diǎn)號(hào)K,KPMAX+1,0,0,0 ! 生成筒體中心線一個(gè)關(guān)鍵點(diǎn)K,KPMAX+2,0,-HC,0 ! 生成筒體中心線另一個(gè)關(guān)鍵點(diǎn)L,KPMAX+1,KPMAX+2 ! 生成筒體中心線 VDRAG,ALL, , , , , ,48 ! 生成筒體 LDELE,1,1 ! 刪除多余線LDELE,48,1 ! 刪除多余線!*劃分網(wǎng)格*LSEL,S,LENGTH,0,TC ! 選擇厚度方向的線段LESIZE,ALL, , ,NT, , , ,1 ! 設(shè)定厚度方向剖分?jǐn)?shù) VSEL,S,LOC,Y,-HC,0 ! 選擇筒體VSWEEP,ALL ! 剖分筒體VSEL,S,LOC,X,0,DI ! 選擇X軸正向的部分封頭及部分中心接管體素VSEL,U,LOC,Y,-HC,0 ! 去掉筒體 LOCAL,11,1,90 ! 定義局部柱坐標(biāo)ASLV ! 按體選面LSLA ! 按面選線LSEL,U,LOC,Z,0 ! 去掉筒體線段LESIZE,ALL, , ,NME, , , , ,1 ! 設(shè)定橢圓封頭經(jīng)線剖分?jǐn)?shù) ALLS!vSEL,S,1!vSEL,a,4!VSWEEP,ALL!ALLS!vSEL,S,3!vSEL,a,6!VSWEEP,ALL!ALLS!vSEL,S,2!vSEL,a,5!VMESH,ALL ! ASEL,S,30 ! 選擇結(jié)構(gòu)的內(nèi)表面ASEL,A,35ASEL,A,5ASEL,A,20ASEL,A,14ASEL,A,15ASEL,A,27ASEL,A,28ASEL,A,24ASEL,A,11CM,AREAINNER,AREA ! 定義結(jié)構(gòu)的內(nèi)表面為組件 CSYS,0 ! 激活總體直角坐標(biāo)系A(chǔ)SEL,S,LOC,Y,KY(8) ! 選擇接管端面NSLA,1 ! 選擇接管端面上的節(jié)點(diǎn)CM,NMOMENT,NODE ! 定義接管端面上的節(jié)點(diǎn)組件NSEL,R,LOC,X,KX(8) ! 選擇接管端面上的一個(gè)節(jié)點(diǎn)NSEL,R,LOC,Z,KZ(8)*GET,NMAST,NODE,NUM,MIN ! 提取該節(jié)點(diǎn)號(hào)TYPE,2 ! 設(shè)定質(zhì)量單元E,NMAST ! 定義質(zhì)量單元ALLS FINI ! 退出前處理!* 求解 */SOLU ! 進(jìn)入求解器ASEL,S,LOC,Y,-HC ! 選擇筒體端面DA,ALL,UY ! 約束軸向位移NSLA,1 ! 按面選節(jié)點(diǎn)NSEL,R,LOC,X,0 ! 再選擇X為0的節(jié)點(diǎn)D,ALL,UX ! 約束X方向位移NSLA,1 ! 按面選節(jié)點(diǎn)NSEL,R,LOC,Z,0 ! 再選擇Z為0的節(jié)點(diǎn)D,ALL,UZ ! 約束Z方向位移ASEL,S,LOC,Z,0 ! 選擇對(duì)稱面DA,ALL,SYMM ! 施加對(duì)稱約束! * 以下定義載荷、加載并求解 *CMSEL,S,AREAINNER ! 選擇內(nèi)表面組件AREAINNER SFA,ALL,1,PRES,PC ! 施加內(nèi)壓ASEL,S,LOC,Y,KY(8) ! 選擇接管端面SFA,ALL,1,PRES,PAX1 ! 對(duì)接管施加端部平衡面載荷CMSEL,S,NMOMENT ! 選擇接管端面上的節(jié)點(diǎn)組件CERIG,NMAST,ALL ! 定義剛性平面F,NMAST,MZ,M ! 施加彎矩ALLS ! 全選 SOLVE ! 求解SAVE ! 保存FINI ! 退出求解器! * 后處理 */POST1 ! 進(jìn)入后處理/GRAPHICS,FULLPLNSOL,S,INT,1 !顯示應(yīng)力云圖!*應(yīng)力評(píng)定*NODE1=NODE(KX(13),KY(13),KZ(13) ! 定義內(nèi)外壁節(jié)點(diǎn)NOD