求助处理固液相变问题时提示的错误，一直无法解决，附有命令流

向左走向右走

采用热焓法处理相变，轴对称模型，选用均匀网格，发生固液相变的温度区间是2摄氏度。
曾经改变过模型的几何大小，网格大小，还有热源的大小，但是还是一直出现下面的错误提示，我真是不知道改该如何解决了，哪请高手们帮我看看是哪里出现了问题，非常感谢！
出错文件是
*** ERROR ***                           CP =       8.891   TIME= 10:26:48
The value of TEMP at node 4 is 2.786279066E+18.  It is greater than the
current limit of 1000000.  This generally indicates rigid body motion   
as a result of an unconstrained model.  Verify that your model is      
properly constrained.

我用GUI做的，生成的命令流为
/BATCH  
/COM,ANSYS RELEASE 10.0    UP20050718       10:17:34    06/12/2008
/TITLE,laser heating
!*  
/NOPR   
/PMETH,OFF,0
KEYW,PR_SET,1   
KEYW,PR_STRUC,0
KEYW,PR_THERM,1
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0   
KEYW,MAGEDG,0   
KEYW,MAGHFE,0   
KEYW,MAGELC,0   
KEYW,PR_MULTI,0
KEYW,PR_CFD,0   
/GO
!*  
/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Thermal  
!*  
/PREP7  
!*  
ET,1,PLANE55
!*  
KEYOPT,1,1,0
KEYOPT,1,3,1
KEYOPT,1,4,0
KEYOPT,1,8,0
KEYOPT,1,9,0
!*  
TOFFST,273  
!*  
!*  
!*  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,KXX,1,,206   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2707
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPTEMP,2,695
MPTEMP,3,697
MPTEMP,4,1000   
MPDATA,ENTH,1,,0
MPDATA,ENTH,1,,1.6857e9
MPDATA,ENTH,1,,2.7614e9
MPDATA,ENTH,1,,3.6226e9
RECTNG,0,4e-3,0,2e-3,   
FLST,5,4,4,ORDE,2   
FITEM,5,1   
FITEM,5,-4  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1,20e-6, , , , , , ,1  
!*  
MSHKEY,0
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
ALLSEL,ALL  
FINISH  
/SOL
!*  
ANTYPE,4
!*  
TRNOPT,FULL
LUMPM,0
!*  
NSUBST,100,1000000,1
OUTRES,ERASE
OUTRES,ALL,1
AUTOTS,1
LNSRCH,1
RESCONTRL,DEFINE,ALL,1,1
TIME,10e-9  
TUNIF,0,
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNCCSYS   
*SET,_FNCNAME,'q'   
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*SET,_FNC_C1(1),150e-3  
*SET,_FNC_C2(1),500e-6  
*SET,_FNCCSYS,0
! /INPUT,heat source function-ns.func,,,1   
*DIM,%_FNCNAME%,TABLE,6,27,1,,,,%_FNCCSYS%  
!   
! Begin of equation: 0.0588*(E/(PI*a^2*10^(-8)))*exp(-{X}^2/a^2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 18
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 18, 17, -1   
*SET,%_FNCNAME%(0,3,1),   0, -1, 0, 3.14159265358979310, 0, 0, -2   
*SET,%_FNCNAME%(0,4,1), 0.0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,6,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,7,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 8, 0, 0, -4
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 10, 0, 0, -2   
*SET,%_FNCNAME%(0,11,1), 0.0, -4, 0, 1, -1, 17, -2  
*SET,%_FNCNAME%(0,12,1), 0.0, -1, 0, 1, -3, 3, -4   
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 1, 17, 4, -1   
*SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 0.0588, 0, 0, -2   
*SET,%_FNCNAME%(0,15,1), 0.0, -3, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, -1, 2, -2   
*SET,%_FNCNAME%(0,19,1), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,20,1), 0.0, -2, 0, 1, 2, 17, -1   
*SET,%_FNCNAME%(0,21,1), 0.0, -1, 0, 1, -4, 3, -2   
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 2, 0, 0, 18
*SET,%_FNCNAME%(0,23,1), 0.0, -4, 0, 1, 18, 17, -2  
*SET,%_FNCNAME%(0,24,1), 0.0, -2, 0, 1, -1, 4, -4   
*SET,%_FNCNAME%(0,25,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,26,1), 0.0, -2, 0, 1, -3, 3, -1   
*SET,%_FNCNAME%(0,27,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 0.0588*(E/(PI*a^2*10^(-8)))*exp(-{X}^2/a^2)  
!-->
LSEL,S, , ,       1
NSLL,S,1
FLST,2,201,1,ORDE,2
FITEM,2,1   
FITEM,2,-201
/GO
!*  
!*  
SF,P51X,HFLUX, %Q%  
ALLSEL,ALL  
/GST,1  
!*  
OUTRES,ALL,1,   
TIMINT,1
TINTP,0.005, , ,1, , ,  
!*  
/STATUS,SOLU
SOLVE   
FINISH

[ 本帖最后由 向左走向右走 于 2008-6-13 10:50 编辑 ]

16443

显热容法处理相变，相变区间取10K。10ns单脉冲激光作用铝靶材，轴对称几何模型4mm*2mm，光斑半径100微米，大网格为100微米，采用局部区域细化网格。发生熔化相变时的平台应当出现在熔点附近吧，但是我做出的结果这个相变平台却高出熔点100多K，下面是我的命令流，有哪位高手帮忙看看有什么办法解决吗？谢谢！

属于同一个问题的内容不要分开发帖

向左走向右走

主任，这个帖子主要是求助这个错误提示如何解决啊！已经很长时间了，一直没找到解决办法，希望大家帮帮忙！
出错文件是
*** ERROR ***                           CP =       8.891   TIME= 10:26:48
The value of TEMP at node 4 is 2.786279066E+18.  It is greater than the
current limit of 1000000.  This generally indicates rigid body motion   
as a result of an unconstrained model.  Verify that your model is      
properly constrained.

yagn

对这个问题比较感兴趣，关注中。。

向左走向右走

希望大家多提建议啊:@L :@P
我是在ANSYS中输入焓随温度变化的值来处理相变的，不知道哪里不对，总出现那个错误提示

[ 本帖最后由 向左走向右走 于 2008-6-14 09:49 编辑 ]