压力容器接管处施加弯矩的处理和内部加内压.doc

! * 设置环境变量 *FINISH/CLEAR/FILN, E65/TITLE, /UNITS,SI !采用国际单位制!*参数定义*DI=2000 ! 设备内径TC=16 ! 筒体及封头壁厚HC=5*NINT(SQRT(DI/2*TC)/10)*10 ! 模型中筒体长度（含封头直边高度)DN=630 ! 接管外径TN=8 ! 接管壁厚LN=350 ! 接管外伸高度 NT=2 ! 厚度方向剖分数NME=30 ! 椭圆封头经向剖分数NH=20 ! 接管环向剖分数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 ! 定义质量单元的实常数! * 建立模型 *CSWPLA,11,1,0.5,1, ! 定义椭圆坐标系K,1,-DI/2, ! 定义椭圆的一个端点K,2,DI/2, ! 定义椭圆的另一个端点L,1,2 ! 生成椭圆线CSYS,0 ! 激活总体直角坐标系K,3,KX(2)+TC, ! 定义封头壁厚关键点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 ! 提取最大节点号K,KPMAX+1,0,0,0 ! 生成中心线一个关键点K,KPMAX+2,0,HC,0 ! 生成中心线另一个关键点VROTAT,ALL,KPMAX+1,KPMAX+2,180 ! 旋转生成椭圆封头及中心接管ASEL,S,LOC,Y,0 ! 选择椭圆封头端面*GET,KPMAX,KP,NUM,MAX ! 提取最大节点号K,KPMAX+1,0,0,0 ! 生成筒体中心线一个关键点K,KPMAX+2,0,-HC,0 ! 生成筒体中心线另一个关键点L,KPMAX+1,KPMAX+2 ! 生成筒体中心线 VDRAG,ALL, , , , , ,48 ! 生成筒体 LDELE,1,1 ! 删除多余线LDELE,48,1 ! 删除多余线!*划分网格*LSEL,S,LENGTH,0,TC ! 选择厚度方向的线段LESIZE,ALL, , ,NT, , , ,1 ! 设定厚度方向剖分数 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 ! 定义局部柱坐标ASLV ! 按体选面LSLA ! 按面选线LSEL,U,LOC,Z,0 ! 去掉筒体线段LESIZE,ALL, , ,NME, , , , ,1 ! 设定椭圆封头经线剖分数 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 ! 选择结构的内表面ASEL,A,35ASEL,A,5ASEL,A,20ASEL,A,14ASEL,A,15ASEL,A,27ASEL,A,28ASEL,A,24ASEL,A,11CM,AREAINNER,AREA ! 定义结构的内表面为组件 CSYS,0 ! 激活总体直角坐标系ASEL,S,LOC,Y,KY(8) ! 选择接管端面NSLA,1 ! 选择接管端面上的节点CM,NMOMENT,NODE ! 定义接管端面上的节点组件NSEL,R,LOC,X,KX(8) ! 选择接管端面上的一个节点NSEL,R,LOC,Z,KZ(8)*GET,NMAST,NODE,NUM,MIN ! 提取该节点号TYPE,2 ! 设定质量单元E,NMAST ! 定义质量单元ALLS FINI ! 退出前处理!* 求解 */SOLU ! 进入求解器ASEL,S,LOC,Y,-HC ! 选择筒体端面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 ! 施加对称约束! * 以下定义载荷、加载并求解 *CMSEL,S,AREAINNER ! 选择内表面组件AREAINNER SFA,ALL,1,PRES,PC ! 施加内压ASEL,S,LOC,Y,KY(8) ! 选择接管端面SFA,ALL,1,PRES,PAX1 ! 对接管施加端部平衡面载荷CMSEL,S,NMOMENT ! 选择接管端面上的节点组件CERIG,NMAST,ALL ! 定义刚性平面F,NMAST,MZ,M ! 施加弯矩ALLS ! 全选 SOLVE ! 求解SAVE ! 保存FINI ! 退出求解器! * 后处理 */POST1 ! 进入后处理/GRAPHICS,FULLPLNSOL,S,INT,1 !显示应力云图!*应力评定*NODE1=NODE(KX(13),KY(13),KZ(13) ! 定义内外壁节点NOD