/TITLE,Steady-state thermal analysis of pipe junction
/UNITS,BIN ! 英制单位；Use U. S. Customary system of units (inches)
! /SHOW, ! Specify graphics driver for interactive run
ET,1,90 ! Define 20-node, 3-D thermal solid element MP,DENS,1,.285 ! Density = .285 lbf/in^3 MPTEMP,,70,200,300,400,500 ! Create temperature table MPDATA,KXX,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12
! 指定与温度相对应的数据材料属性；导热系数；Define conductivity values
MPDATA,C,1,,.113,.117,.119,.122,.125
! Define specific heat values（比热）
MPDATA,HF,2,,426/144,405/144,352/144,275/144,221/144
! Define film coefficient;除144是单位问题，上面的除12也是单元问题
! Define parameters for model generation RI1=1.3 ! Inside radius of cylindrical tank RO1=1.5 ! Outside radius Z1=2 ! Length
RI2=.4 ! Inside radius of pipe RO2=.5 ! Outside pipe radius Z2=2 ! Pipe length
CYLIND,RI1,RO1,,Z1,,90 ! 90 degree cylindrical volume for tank
WPROTA,0,-90 ! 旋转当前工作的平面；从Y到Z旋转-90度；；Rotate working plane to pipe axis
CYLIND,RI2,RO2,,Z2,-90 ! 角度选择在了第四象限；90 degree cylindrical volume for pipe
WPSTYL,DEFA ! 重新安排工作平面的设置；另外
WPSTYL,STAT to list the status of the working plane；；Return working plane to default setting
BOPT,NUMB,OFF ! 关掉布尔操作的数字警告信息；Turn off Boolean numbering warning
VOVLAP,1,2 ! 交迭体；Overlap the two cylinders
/PNUM,VOLU,1 ! 体编号打开；Turn volume numbers on /VIEW,,-3,-1,1
/TYPE,,4 ! 精确面的显示；Precise hidden display
/TITLE,Volumes used in building pipe/tank junction VPLOT
VDELE,3,4,,1 ! 修剪一些体与体相关的体的因素都删掉；Trim off excess volumes
! Meshing 网格划分
ASEL,,LOC,Z,Z1 ! Select area at remote Z edge of tank ASEL,A,LOC,Y,0 ! Select area at remote Y edge of tank CM,AREMOTE,AREA ! 为面建立数组；Create area component called AREMOTE
/PNUM,AREA,1 /PNUM,LINE,1
/TITLE,Lines showing the portion being modeled APLOT
/NOERASE ! 预防抹去
LPLOT ! Overlay line plot on area plot /ERASE
ACCAT,ALL ! 连接面和线的准备映射；Concatenate areas and lines at remote tank edges LCCAT,12,7 LCCAT,10,5
LESIZE,20,,,4 ! 4 divisions through pipe thickness LESIZE,40,,,6 ! 6 divisions along pipe length
LESIZE,6,,,4 ! 4 divisions through tank thickness ALLSEL ! Restore full set of entities
ESIZE,.4 ! Set default element size线的默认划分数
MSHAPE,0,3D ! Choose mapped brick mesh MSHKEY,1 ! 映射网格
SAVE ! Save database before meshing
VMESH,ALL ! Generate nodes and elements within volumes
/PNUM,DEFA ! 重新安排数字规格 /TITLE,Elements in portion being modeled EPLOT
FINISH
/COM, *** Obtain solution *** /SOLU
ANTYPE,STATIC ! Steady-state analysis type
NROPT,AUTO ! 自动选择牛顿-拉普森 Program-chosen Newton-Raphson option
TUNIF,450 ! 给结点统一的温度；Uniform starting temperature at all nodes
CSYS,1 ! 1 — Cylindrical with Z as the axis of rotation
NSEL,S,LOC,X,RI1 ! Nodes on inner tank surface
SF,ALL,CONV,250/144,450 ! 为结点指定表面载荷；对流；Convection（对流）； load at all nodes
CMSEL,,AREMOTE ! 选择子集组合；Select AREMOTE component
NSLA,,1 ! Nodes belonging to AREMOTE
D,ALL,TEMP,450 ! 设定边界温度条件 Temperature constraints at those nodes
WPROTA,0,-90 ! Rotate working plane to pipe axis CSWPLA,11,1 ! 在工作区声明本地的圆柱体系； Define local cylindrical c.s at working plane
NSEL,S,LOC,X,RI2 ! Nodes on inner pipe surface
SF,ALL,CONV,-2,100 ! 这里的-2表示材料2；；Temperature-dep. convection load at those nodes
ALLSEL
/PBC,TEMP,,1 ! 边界符号的显示Temperature b.c. symbols on
/PSF,CONV,,2 ! Convection symbols on 箭头显示 /TITLE,Boundary conditions NPLOT
WPSTYL,DEFA CSYS,0
AUTOTS,ON ! Automatic time stepping NSUBST,50 ! Number of substeps KBC,0 ! Ramped loading (default)
OUTPR,NSOL,LAST ! 显示最后一次的结点约束；Optional command for solution printout SOLVE
FINISH
/COM, *** Review results *** /POST1
/EDGE,,1 ! Displays only the \"edges（刀口, 利刃, 锋, 优势, 边缘, 优势, 尖锐）\" of an object；Edge display /PLOPTS,INFO,ON ! Legend column on
/PLOPTS,LEG1,OFF ! Legend header off 圆柱数列的头部
/WINDOW,1,SQUARE ! SQUA, form largest square window within the current graphics area；Redefine window size /TITLE,Temperature contours at pipe/tank junction
PLNSOL,TEMP ! Plot temperature contours
CSYS,11
NSEL,,LOC,X,RO2 ! Nodes and elements at outer radius of pipe ESLN ! 选择单元 NSLE ! 选择结点
/SHOW,,,1 ! 向量显示；Vector mode /TITLE,Thermal flux vectors at pipe/tank junction
PLVECT,TF ! Plot thermal flux（热通量） vectors
FINISH
因篇幅问题不能全部显示，请点此查看更多更全内容