【求助】模拟焊接温度分布,求解过程出现这种错误提示

我用一块平板,在上面堆焊,模拟焊接温度分布,求解过程出现这种错误提示,弄了半天,也不知是怎么回事,请问各位这个是什么原因?
以下是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
*DIM,%_FNCNAME%,TABLE,6,19,1,,,,%_FNCCSYS%
我也来说两句 查看全部回复

最新回复

  • zouyou (2015-10-24 10:50:36)

    获得机床从开始运行到6个小时后的温升和主轴受热变形的位移情况,进行热误差补偿建模。
  • yuanyuan (2015-10-24 10:51:05)

    和最后一个时间点上的位移值基本相同。并不是从一开始的几微米慢慢增加到几十微米。为什么会是这样
  • 女儿情 (2015-10-24 10:51:38)

    我也遇到了相同的问题,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
    *enddo
  • 跳跳哈里 (2015-10-24 10:54:02)

    我个人认为需要注意两点:
    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
  • 大花猫bb (2015-10-24 10:54:51)

    意思就是前面的time,14400和deltim,300都删掉还是怎么改?
  • qinqinai (2015-10-24 10:55:30)

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

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

    你检查你的温度结果了吗?不同时间的温度值相差很大吗?
  • 886爱 (2015-10-24 10:57:18)

    我检查了结果的,不同时间,三个方向的位移相差很小,并且不同时刻的位移结果和最后一个时刻的位移结果相差很小,差值在5微米之内,这不符合实际的结果:(
  • zouyou (2015-10-24 10:58:16)

    不同时间的温度结果是符合实际情况,从20度慢慢升到27°左右,最后趋于水平。我觉得不同时刻的位移结果应该是从最开始的几乎接近与零开始慢慢增加才是正确的,不知道我这样考虑是不是正确的
  • efp (2015-10-24 10:58:49)

    有两个物体你考虑一下:
    1  材料属性中,定义参考温度是 30 度,是否应该改为 20 度;
    2  好像需要对结构分析也设置一个参考温度 (默认是零度),也应该是 20 度。
        没有设置参考温度,那么默认是 零 度,后面分析时,温度从 0 -> 20,和从 0 ->27 ,其温度变化差不多吧?