热分析

我用一块平板，在上面堆焊，模拟焊接温度分布，求解过程出现这种错误提示，弄了半天，也不知是怎么回事，请问各位这个是什么原因？
以下是file里的命令：不好意思，我的金币数不多，跪谢了！
BATCH  
/COM,ANSYS RELEASE 14.0    UP20111024       16:42:48    05/05/2014
RESUME,'ciaochicun25','db','H:/ANSYS'   
/COM,ANSYS RELEASE 14.0    UP20111024       16:42:48    05/05/2014
/SOLU   
!*  
TUNIF,20,   
FLST,2,18,5,ORDE,6  
FITEM,2,1   
FITEM,2,-4  
FITEM,2,7   
FITEM,2,-9  
FITEM,2,11  
FITEM,2,-21
/GO
!*  
SFA,P51X,1,CONV,30,20   
FLST,2,2,5,ORDE,2   
FITEM,2,5   
FITEM,2,10  
/GO
!*  
!*  
SFA,P51X,1,HFLUX, %GAOSI%   
!*  
TIME,L/V
AUTOTS,-1   
NSUBST,50,50,50,1   
KBC,0   
!*  
TSRES,ERASE
LSWRITE,1,  
!*  
TIME,50
AUTOTS,-1   
NSUBST,30,30,30,1   
KBC,0   
!*  
TSRES,ERASE
LSWRITE,2,  
!*  
TIME,1100   
AUTOTS,-1   
NSUBST,105,105,105,1
KBC,0   
!*  
TSRES,ERASE
LSWRITE,3,  
! solve
/BATCH  
/COM,ANSYS RELEASE 14.0    UP20111024       21:30:09    05/05/2014
/input,menust,tmp,'',,,,,,,,,,,,,,,,1   
/COM,ANSYS RELEASE 14.0    UP20111024       21:30:10    05/05/2014
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1   
/REPLOT,RESIZE  
WPSTYLE,,,,,,,,0
/REPLOT,RESIZE  
EPLOT   
APLOT   
LPLOT   
APLOT   
!*  
/NOPR   
KEYW,PR_SET,1   
KEYW,PR_STRUC,1
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,1
KEYW,PR_CFD,0   
/GO
!*  
/COM,   
/COM,Preferences for GUI filtering have been set to display:
/COM,  Structural   
/COM,  Thermal  
!*  
/PREP7  
!*  
LPLOT   
/REPLOT
LPLOT   
/REPLOT
!*  
FINISH  
! /EXIT,NOSAV   
/BATCH  
/COM,ANSYS RELEASE 14.0    UP20111024       13:57:33    05/06/2014
/input,menust,tmp,'',,,,,,,,,,,,,,,,1   
/COM,ANSYS RELEASE 14.0    UP20111024       13:57:33    05/06/2014
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1   
/REPLOT,RESIZE  
WPSTYLE,,,,,,,,0
APLOT   
APLOT   
/REPLOT
/PREP7  
!*  
!*  
K,1,0,0,0,  
*SET,I        , 160
*SET,H        , 5000000000E-03  
*SET,W        ,0.1  
*SET,Q        , U*I*YITA
*SET,H        , 0.006   
*SET,QM       , 27284509.2  
/PNUM,KP,1  
/PNUM,LINE,0
/PNUM,AREA,1
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0   
!*  
/PNUM,ELEM,0
/REPLOT
!*  
K,2,0,0.1,0,
K,3,-0.05*0.15,0.1,0,   
K,4,-0.05*0.3,0.1,0,
K,5,-0.05*5,0.1,0,  
KDELE,       5  
K,5,-0.05*0.5,0.1,0,
K,6,-0.05,0.1,0,
K,7,-0.05,0,0,  
K,8,-0.05*0.5,0,0,  
K,8,-0.05*0.3,0,0,  
KDELE,       8  
K,8,-0.05*0.5,0,0,  
K,9,-0.05*0.3,0,0,  
K,10,-0.05*0.15,0,0,
K,10,-0.05*0.15,0,0,
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   
FLST,2,4,3  
FITEM,2,2   
FITEM,2,1   
FITEM,2,10  
FITEM,2,3   
A,P51X  
FLST,2,4,3  
FITEM,2,3   
FITEM,2,4   
FITEM,2,9   
FITEM,2,10  
A,P51X  
FLST,2,4,3  
FITEM,2,5   
FITEM,2,4   
FITEM,2,9   
FITEM,2,8   
A,P51X  
FLST,2,4,3  
FITEM,2,5   
FITEM,2,6   
FITEM,2,7   
FITEM,2,8   
A,P51X  
ESIZE,0.0012,0,
MSHAPE,0,2D
MSHKEY,0
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
ESIZE,0.0025,0,
CM,_Y,AREA  
ASEL, , , ,       2
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
ESIZE,0.005,0,  
CM,_Y,AREA  
ASEL, , , ,       3
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
ESIZE,0.0065,0,
CM,_Y,AREA  
ASEL, , , ,       4
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
TYPE,   2   
EXTOPT,ESIZE,4,0,   
EXTOPT,ACLEAR,0
!*  
EXTOPT,ATTR,0,0,0   
MAT,1   
REAL,_Z4
ESYS,0  
!*  
FLST,2,4,5,ORDE,2   
FITEM,2,1   
FITEM,2,-4  
VEXT,P51X, , ,0,0,0.006,,,,
NUMMRG,ALL, , , ,LOW
FINISH  
/SOL
!*  
ANTYPE,4
!*  
TRNOPT,FULL
LUMPM,0
!*  
TUNIF,20,   
FLST,2,18,5,ORDE,6  
FITEM,2,1   
FITEM,2,-4  
FITEM,2,7   
FITEM,2,-9  
FITEM,2,11  
FITEM,2,-21
/GO
!*  
SFA,P51X,1,CONV,30,20   
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNCCSYS   
*SET,_FNCNAME,'GAOSI'   
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*SET,_FNC_C1(1),Qm  
*SET,_FNC_C2(1),V   
*SET,_FNC_C3(1),R   
*SET,_FNCCSYS,0
! /INPUT,gaosi.func,,,1

获得机床从开始运行到6个小时后的温升和主轴受热变形的位移情况，进行热误差补偿建模。

和最后一个时间点上的位移值基本相同。并不是从一开始的几微米慢慢增加到几十微米。为什么会是这样

我也遇到了相同的问题，bfv,1,hgen,20000     ！定义电机发热功率
bfv,14,hgen,642857   ！定义上轴承发热功率
bfv,20,hgen,353294   ！定义下轴承发热功率
CMSEL,S,CM_AXLE         ！定义主轴表面换热系数
sfa,all,,conv,45,20   
cmsel,s,cm_base          ！定义基座换热系数
sfa,all,,conv,3,20      
cmsel,s,cm_column        ！定义立柱换热系数
sfa,all,,conv,4,20
cmsel,s,cm_lowercover     ！定义轴承端盖换热系数
sfa,all,,conv,13,20
cmsel,s,cm_uppercover
sfa,all,,conv,10,20
cmsel,s,cm_lower_bear    ！定义轴承换热系数
sfa,all,,conv,40,20
cmsel,s,cm_upper_bear
sfa,all,,conv,35,20
cmsel,s,cm_motor           ！定义电机换热系数
sfa,all,,conv,150,20
cmsel,s,cm_spindle_box      ！定义主轴箱换热系数
sfa,all,,conv,20,20
cmsel,s,cm_worktable        ！定义工作台换热系数
sfa,all,,conv,4.5,20
allsel
antype,trans         
timint,on
time,14400   
deltim,300
kbc,1
tunif,20
outres,all,all
solve
finish
/PREP7  
etchg,tts   !转为结构分析
/SOL
antype,static
da,30,all,0
*do,jj,300,14400,300
     Time,jj
    ldread,temp,,,jj,,'file','rth'
solve

我个人认为需要注意两点：
1 保证各种单位的统一。因为你的位移是微米级别，与机床尺寸不是一个数量级；
2 在保证1的前提下，问题可能出现在最后一段：
*do,jj,300,14400,300       !可能出现问题的地方
     Time,jj
    ldread,temp,,,jj,,'file','rth'
solve
*enddo
我的建议：
将前面的 “time,14400 $ deltim,300”等有关时间、时间步的参数统统设为默认。
这样TIME变为一个简单的计数器，默认从1开始，后面依次+1。不容易出错。
那么最后一段需要相应的改为：
*do,jj,1,20,1       !假设laststep为20
     Time,jj
    ldread,temp,,,jj,,'file','rth'
solve
*enddo
我前面也做过类似的分析，一般都是按上述方式处理的，比较顺利。

意思就是前面的time,14400和deltim,300都删掉还是怎么改？

果与实际情况是不一样的。
    楼主应该首先检查一下温度的计算结果，看不同时间的结构温度差别有多大，如果温度差别本来就不大，那么位移结果相差也不会很大。
    另外，楼主的位移计算部分，除了温度载荷外，没有设置任何约束条件，也可能导致位移结果不正确。需要添加合适的约束条件，然后重新计算。

对一个面进行了全约束，和实际约束基本是相符的，我表达的位移结果相差不大是指从计算一开始到计算结束，所有的位移都在一个数量上，比如从开始到结束所有的位移都是0.312e-4到0.313e-4之间，只是后面的小数多了而已。

你检查你的温度结果了吗？不同时间的温度值相差很大吗？

我检查了结果的，不同时间，三个方向的位移相差很小，并且不同时刻的位移结果和最后一个时刻的位移结果相差很小，差值在5微米之内，这不符合实际的结果:(