向左走向右走 发表于 2008-6-13 10:48

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

采用热焓法处理相变,轴对称模型,选用均匀网格,发生固液相变的温度区间是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 发表于 2008-6-13 12:51

回复 楼主 的帖子

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

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

向左走向右走 发表于 2008-6-13 19:52

回复 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.

yagn 发表于 2008-6-14 09:40

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

向左走向右走 发表于 2008-6-14 09:43

回复 4楼 的帖子

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

[ 本帖最后由 向左走向右走 于 2008-6-14 09:49 编辑 ]
页: [1]
查看完整版本: 求助处理固液相变问题时提示的错误,一直无法解决,附有命令流