新人求助，关于气割过程中温度场的影响

cumt410

小弟学ansys不久，最近需要研究气割对结构下部冻土的温度场影响，参照前人写了段命令流，动态加载不太会，选择了一部分节点施加一个热生成荷载，可是总是运行不正确，请各位前辈不惜赐教，在下感激不尽FINI
FINI
/clear
/config,nproc,4  !4核cpu计算


KXX1=1.96!(秒)
KXX2=1.96
KXX3=1.72
KXX4=1.72
KXX5=1.72
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
w=30/100                !地层含水量
wu=0/100                !土层未冻水含量
rou0=1.32*1000          !土干密度


ENTH1=0
ENTH2=81957
ENTH3=197764
ENTH4=292092


!!!!!!尝试修改比热!!!!!!!!
C1=1690
C2=1690
C3=1690
C4=1690
TEMP0=20   !初始土体温度
TIMENUM=60  !冻结天数
TDD=-0.4         !土层的结冰温度,℃
DDT=4             !相变区间的宽度,℃,相变区间越小,计算收敛的越慢,建议不小于1℃






!!!!!!!!!!!!数值计算区
!土层的结冰潜热计算
Q=rou0*(W-WU)*334.56*1000    !土层的结冰潜热,334.56（单位kj/kg）为水的结冰潜热,  J/立方米


!土的湿密度计算
ROUD=ROU0*(1+W)          !土层的湿密度,kg/立方米


CDU=1690               !未冻土的比热,
CDF=1160               !冻土的比热,J/(kg·℃),


!焓值计算
ENTH1=0                                 !-65℃的焓值,为相对值,取为0kJ/立方米
ENTH2=ROUD*CDF*(TDD-DDT/2+180)                !相变区冷端温度时的焓值,J/立方米
ENTH3=ENTH2+Q                               !相变区热端温度时的焓值,J/立方米
ENTH4=ENTH3+ROUD*CDU*(20-(TDD+DDT/2))      !20℃时的焓值,kJ/立方米
ENTH5=ENTH3+ROUD*CDU*(100-(TDD+DDT/2))      !100℃时的焓值,kJ/立方米


*DIM,STEPTIME,ARRAY,TIMENUM        !定义时间数组
*VFILL,STEPTIME(1),DATA,1,1,1,1,1,1,1,1,1,1      
!*VFILL,STEPTIME(11),DATA,1,1,1,1,1,1,1,1,1,1
!*VFILL,STEPTIME(21),DATA,1,1,1,1,1,1,1,1,1,1
!*VFILL,STEPTIME(31),DATA,1,1,1,1,1,1,1,1,1,1


*DIM,DJ_TIME,ARRAY,TIMENUM
!定义冻结时间数据数组（按时间步提取计算结果时用）
*VFILL,DJ_TIME(1),DATA,1,2,3,4,5,6,7,8,9,10
!*VFILL,DJ_TIME(11),DATA,11,12,13,14,15,16,17,18,19,20
!*VFILL,DJ_TIME(21),DATA,21,22,23,24,25,26,27,28,29,30
!*VFILL,DJ_TIME(31),DATA,31,32,33,34,35,36,37,38,39,40




!!!!!!!!!!!!!!!!!!!!!!!!!!!建模
/prep7


R1=10       !!!!!!!!!!!!盾壳内壁径
DD=1        !!!!!!!!!!!!盾壳厚度
R2=20        !!!!!!!!!!!!冻土体外边界
L=10         !!!!!!!!!!!!研究深度
DX=0.8       !!!!!!!!!!!!胸板厚度
HK=0.10      !!!!!!!!!!!!气焊孔直径
BL=0.5      !!!!!!!!!!!!胸板保留长度
HJ=1000      !!!!!!!!!!!!焊接热生成率
xx=0.5     !!!!!!!!!!!!!!割缝直径
TIMENUM=10    !!!!!!!!!!!定义加载时间




ET,1,solid70
MPTEMP,1,-20,TDD-DDT/2,TDD+DDT/2,20,100, ,   !为材料属性定义一个温度表  MPTEMP,STLOC,T1,T2,T3,T4,T5,T6.! 结冰温度±半个相变区间 !置，缺省值为最后填充位置+1
MPDATA,KXX,1,1,KXX1,KXX2,KXX3,KXX4,KXX5 , ,         !kJ/(m?d?℃)
MPDATA,ENTH,1,1,ENTH1,ENTH2,ENTH3,ENTH4,ENTH5 , ,     !kJ/立方米
MPDATA,C,1,1,C1,C2,C3,C4, , ,     
!kJ/kg?℃
mpdata,dens,1,1.85


ET,2,solid70  !!!!!!!!!!!!!!钢板的材料属性
MPDATA,DENS,2,,7.8
MPTEMP,,,,,,,,  
MP,KXX,2,1.296
MP,C,2,0.465






CYL4,0,0,R1,-20,R1+DD,20,L
CYL4,0,0,R1+DD,-20,R2,20,L
wpoff,0,0,L/2-DX
CYL4,0,0,R1/2,-20,R1,20,DX


csys,1




!csys,0
!wpoff,R1-BL,0,0
!CYL4,0,0,HK,,,,DX


!allsel
!Vglue,all   


!VOVLAP,all
NUMCMP,all




allsel
vglue,all
CSYS,0  
WPAVE,0,0,0
WPSTYLE,,,,,,,,0
allsel
lsel,s,loc,x,0,R1+DD
lesize,all,0.1


!allsel
!lsel,s,loc,x,R1,R2
!lesize,all,0.1






csys,1
vsel,s,loc,x,0,R1
cm,gb,volu
TYPE,2                  !给生成的单元激活一个单元号，即把ET中定义的单元确定
VSWEEP,ALL               !对所有的体进行扫掠式网格划分




allsel
cmsel,u,gb
TYPE,1                  !给生成的单元激活一个单元号，即把ET中定义的单元确定
!Vsweep,ALL
vmesh,all            
allsel
NUMCMP,all


!!!!!!!!!!!!!!!!!!!!!!!!!!加载荷载


/SOLU


ANTYPE,4       !分析类型为瞬态分析
LNSRCH,1       !打开线形搜索
TRNOPT,FULL    !瞬态分析采用完全法（默认）




TUNIF,20                 !设定初始温度
ALLSEL




csys,1
allsel
vsel,s,loc,x,0,R1+DD
NSLv,S,1
D,all,temp,20


allsel
vsel,s,loc,x,R1+DD,R2
NSLv,S,1
D,all,temp,-15




*DO,I,1,TIMENUM,1  


    TIME,DJ_TIME(I)
    AUTOTS,1                !使用自动时间步长跟踪
    DELTIM,STEPTIME(I),,,OFF  !在本荷载步中指定时间步长大小
    KBC,1                     !0为递增荷载 1为阶跃荷载
    EQSLV,jCG,0.1,2           !指定一个方程求解器
     NEQIT,150,  






!csys,0
!wpoff,0,0,L/2-DX
!wpoff,R1-BL,0,0
!nsel,s,loc,x,0,0.5   !选择要加载热源的节点
!nsel,r,loc,y,0,0.5
!nsel,r,loc,z,0,DX
!D,all,temp,300    !温度 荷载
!bf,all,HGEN,HJ    !或者热生成荷载


!选择要加载对流载荷的表面
ASEL,s,,,15,16,1
ASEL,a,,,5,6,1
NSLA,S,1
sf,all,conv,-1,20 !在节点上施加对流载荷(与空气的对流换热系数）








   ALLSEL


SOLVE
*ENDDO


FINISH  
SAVE,coolman,db,,all   


/POST1  
set,,,1,,15
plnsol,temp,,0