/PREP7
/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
