多层多道堆焊问题

xufujia

再次求助各位大侠，程序能够运行 就是每一步的结果都是室温，感觉热流密度没有加上，换成平均热流密度值也不行，麻烦各位高手给小弟 指点一下 不胜感激！！！
/CONFIG,NRES,12000
/FILNAM,U50
/TITLE,PARALLEL LONG RASTER  
/PREP7
/units,si  !采用SI国际单位
RF=0.005
P=3000
L=160E-3
W=110E-3
TS=11E-3
LD=80E-3
WD=50E-3
TD=6E-3
NNY=1.5E-3
NNX=1E-3
NNZ=2E-3
WW=(W-WD)/2
LL=(L-LD)/2
BEAD=10E-3
beadnum=5
layer=4                                !层数
Temp=20                                        !环境温度
InitTemp=20
CoffConv=30                                !对流换热系数
                             !初始温度
Velocity=0.003                                !焊接速度
stepdis=0.004                                !每个载荷步位移                                    
StepTime=stepdis/Velocity                !每个载荷步时间
TotalTime=(LD+stepdis)/Velocity        !载荷持续时间(焊接一层)
stepnum=LD/stepdis        !载荷步数
ET,1,SOLID70
toffst,273    !以下定义材料参数
MPTEMP,,,,,,,,  
MPTEMP,1,0                              
MPDATA,DENS,1,,7850
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,300
MPTEMP,3,900
MPTEMP,4,1400   
MPTEMP,5,2000   
MPDATA,KXX,1,,18.6  
MPDATA,KXX,1,,21.4  
MPDATA,KXX,1,,28.4  
MPDATA,KXX,1,,33.9  
MPDATA,KXX,1,,48
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,600
MPTEMP,3,800
MPTEMP,4,1400   
MPTEMP,5,2000   
MPDATA,C,1,,502
MPDATA,C,1,,612
MPDATA,C,1,,635
MPDATA,C,1,,659
MPDATA,C,1,,670
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,1300   
MPTEMP,3,1410   
MPTEMP,4,1440   
MPTEMP,5,1550   
MPTEMP,6,2000   
MPDATA,ENTH,1,,7.88e7   
MPDATA,ENTH,1,,6.131e9  
MPDATA,ENTH,1,,7.347e9  
MPDATA,ENTH,1,,9.145e9  
MPDATA,ENTH,1,,1.03e10  
MPDATA,ENTH,1,,1.272e10
NROPT,FULL,,ON
K,1,0,0,0
K,2,WD,0,0
K,3,WD,TS,0
K,4,0,TS,0
K,5,0,TS+TD,0
K,6,WD,TS+TD,0
K,7,WD+WW,0,0
K,8,WD+WW,TS,0
K,9,-WW,0,0
K,10,-WW,TS,0
A,1,2,3,4
A,2,7,8,3
A,9,1,4,10
A,6,5,4,3
/VIEW,1,,77,,77,,77
AGEN,2,1,3,1,,,LL
AGEN,2,1,3,1,,,-LD
*DO,I,1,4,1
VOFFST,I,-LD
*ENDDO
*DO,I,5,7,1
VOFFST,I,-LL
*ENDDO
*DO,I,8,10,1
VOFFST,I,-LL
*ENDDO
VGLUE,all
NUMMRG,ALL
LSEL,S,LENGTH,,LD
LESIZE,ALL,NNZ,,,,1
ALLSEL,ALL
LSEL,S,LENGTH,,WD
LESIZE,ALL,NNX,,,,1
ALLSEL,ALL
LSEL,S,LENGTH,,TD
LESIZE,ALL,NNY,,,,1
ALLSEL,ALL
LSEL,S,LINE,,7,8,1
LSEL,A,LINE,,20,21,1
LESIZE,ALL,,,10,0.4,,1
ALLSEL,ALL
LSEL,S,LINE,,5,9,4
LSEL,A,LINE,,18,22,4
LSEL,A,LINE,,134,139,5
LSEL,A,LINE,,120,125,5
LSEL,A,LINE,,121,124,3
LSEL,A,LINE,,135,138,3
LESIZE,ALL,,,10,2.5,,1
ALLSEL,ALL
LSEL,S,LINE,,115,116,1
LSEL,A,LINE,,118,119,1
LSEL,A,LINE,,114,117,3
LSEL,A,LINE,,122,123,1
LESIZE,ALL,,,10,1.33,,1
ALLSEL,ALL
LSEL,S,LINE,,132,133,1
LSEL,A,LINE,,129,130,1
LSEL,A,LINE,,128,131,3
LSEL,A,LINE,,136,137,1
LESIZE,ALL,,,10,1.33,,1
LSEL,S,LENGTH,,TS
LSEL,R,LOC,X,1E-3,-100E-3
LESIZE,ALL,,,5,0.2,1
ALLSEL,ALL
LSEL,S,LENGTH,,TS
LSEL,R,LOC,X,0,-100E-3
LESIZE,ALL,,,5,5,1
ALLSEL,ALL
LSEL,S,LENGTH,,TS
LSEL,R,LOC,X,1E-3,100E-3
LESIZE,ALL,,,5,0.2,1
ALLSEL,ALL
VMESH,1,16,15
VMESH,13,15,2
VMESH,18,19,1
VMESH,11,12,1
VMESH,14,17,3
ACLEAR,ALL
NUMCMP,ALL
ALLSEL
/solu  !进入运算器
antype,trans  !定义求解选项
trnopt,full
nropt,auto, ,
solcontrol,on
autots,on  !自动步长
timint,on   
wpoffs,45E-3,12.5E-3,0
wprota,0,0,90
wprota,0,-90,0
CSYS,4  !激活4号坐标系
CSWPLA,11,0,1,1,  !在4号坐标系内建立局部坐标系
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'abj'
*SET,_FNCCSYS,0
! /INPUT,lengque.func   
*DIM,%_FNCNAME%,TABLE,6,28,1,,,,%_FNCCSYS%  
!   
! Begin of equation: 10+0.6*5.67*10^(-8)*(({TEMP}+273)^2+293^2)*({TEMP}+273+   
! 293)  
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 0.6, 0, 0, 0   
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 5.67, 0, 0, -1  
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -1, 2, -2
*SET,%_FNCNAME%(0,7,1), 0.0, -1, 0, 8, 0, 0, -4
*SET,%_FNCNAME%(0,8,1), 0.0, -2, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 10, 0, 0, -2
*SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, -1, 17, -2  
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, -3, 3, -4   
*SET,%_FNCNAME%(0,12,1), 0.0, -2, 0, 273, 0, 0, 5   
*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 1, 5, 1, -2
*SET,%_FNCNAME%(0,14,1), 0.0, -2, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,15,1), 0.0, -4, 0, 1, -3, 17, -2  
*SET,%_FNCNAME%(0,16,1), 0.0, -2, 0, 293, 0, 0, 0   
*SET,%_FNCNAME%(0,17,1), 0.0, -3, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,18,1), 0.0, -5, 0, 1, -2, 17, -3  
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 1, -4, 1, -5   
*SET,%_FNCNAME%(0,20,1), 0.0, -3, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,21,1), 0.0, -1, 0, 273, 0, 0, 5   
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, 5, 1, -1
*SET,%_FNCNAME%(0,23,1), 0.0, -1, 0, 293, 0, 0, -2  
*SET,%_FNCNAME%(0,24,1), 0.0, -4, 0, 1, -2, 1, -1   
*SET,%_FNCNAME%(0,25,1), 0.0, -1, 0, 1, -3, 3, -4   
*SET,%_FNCNAME%(0,26,1), 0.0, -2, 0, 10, 0, 0, -1   
*SET,%_FNCNAME%(0,27,1), 0.0, -3, 0, 1, -2, 1, -1   
*SET,%_FNCNAME%(0,28,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 10+0.6*5.67*10^(-8)*(({TEMP}+273)^2+293^2)*({TEMP}+273+293)
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNCCSYS   
*SET,_FNCNAME,'XU'  
*DIM,_FNC_C1,,3
*DIM,_FNC_C2,,3
*SET,_FNC_C1(2),P
*SET,_FNC_C2(2),RF  
*SET,_FNCCSYS,11
! /INPUT,xu.func,,,1
*DIM,%_FNCNAME%,TABLE,6,21,3,,,,%_FNCCSYS%  
!   
! Begin of equation: sqrt({X}^2+{Y}^2)  
*SET,%_FNCNAME%(0,0,1), 0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*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, 2  
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, 2, 17, -1
*SET,%_FNCNAME%(0,3,1),   0, -1, 0, 2, 0, 0, 3  
*SET,%_FNCNAME%(0,4,1), 0.0, -3, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,5,1), 0.0, -1, 0, 1, -2, 1, -3
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 16, 1, -1, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, 99, 0, 1, -1, 0, 0
*SET,%_FNCNAME%(0,8,1),   0
*SET,%_FNCNAME%(0,9,1),   0
*SET,%_FNCNAME%(0,10,1),   0
*SET,%_FNCNAME%(0,11,1),   0
*SET,%_FNCNAME%(0,12,1),   0
*SET,%_FNCNAME%(0,13,1),   0
*SET,%_FNCNAME%(0,14,1),   0
*SET,%_FNCNAME%(0,15,1),   0
*SET,%_FNCNAME%(0,16,1),   0
*SET,%_FNCNAME%(0,17,1),   0
*SET,%_FNCNAME%(0,18,1),   0
*SET,%_FNCNAME%(0,19,1),   0
*SET,%_FNCNAME%(0,20,1),   0
*SET,%_FNCNAME%(0,21,1),   0
! End of equation: sqrt({X}^2+{Y}^2)
!   
! Begin of equation: P/{PI}/RF/RF*exp(-3*({X}^2+{Y}^2)/(RF^2))  
*SET,%_FNCNAME%(0,0,2), 0.0055, -999
*SET,%_FNCNAME%(2,0,2), 0.0
*SET,%_FNCNAME%(3,0,2), %_FNC_C1(2)%
*SET,%_FNCNAME%(4,0,2), %_FNC_C2(2)%
*SET,%_FNCNAME%(5,0,2), 0.0
*SET,%_FNCNAME%(6,0,2), 0.0
*SET,%_FNCNAME%(0,1,2), 1.0, -1, 0, 3.14159265358979310, 0, 0, 17   
*SET,%_FNCNAME%(0,2,2), 0.0, -2, 0, 1, 17, 4, -1
*SET,%_FNCNAME%(0,3,2),   0, -1, 0, 1, -2, 4, 18
*SET,%_FNCNAME%(0,4,2), 0.0, -2, 0, 1, -1, 4, 18
*SET,%_FNCNAME%(0,5,2), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,6,2), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,7,2), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,8,2), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,9,2), 0.0, -3, 0, 1, -4, 3, -1
*SET,%_FNCNAME%(0,10,2), 0.0, -1, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,11,2), 0.0, -4, 0, 1, 2, 17, -1   
*SET,%_FNCNAME%(0,12,2), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,13,2), 0.0, -5, 0, 1, 3, 17, -1   
*SET,%_FNCNAME%(0,14,2), 0.0, -1, 0, 1, -4, 1, -5   
*SET,%_FNCNAME%(0,15,2), 0.0, -4, 0, 1, -3, 3, -1   
*SET,%_FNCNAME%(0,16,2), 0.0, -1, 0, 2, 0, 0, 18
*SET,%_FNCNAME%(0,17,2), 0.0, -3, 0, 1, 18, 17, -1  
*SET,%_FNCNAME%(0,18,2), 0.0, -1, 0, 1, -4, 4, -3   
*SET,%_FNCNAME%(0,19,2), 0.0, -1, 7, 1, -1, 0, 0
*SET,%_FNCNAME%(0,20,2), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,21,2), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: P/{PI}/RF/RF*exp(-3*({X}^2+{Y}^2)/(RF^2))
csys,0  
NSEL, S, LOC, Y, 0, TS                      !基底初始温度
IC,ALL,TEMP,InitTemp
ALLSEL,ALL
NSEL, S, LOC, Z, 40E-3                      !基体四个侧面
NSEL, a, LOC, Z, -120E-3
NSEL, a, LOC, X, -40E-3
NSEL, a, LOC, X, 80E-3
NSEL, a, LOC, Y, 0
NSEL, a, LOC, Y, TS                         !基底底面，换热边界条件
SF, ALL, CONV, CoffConv, Temp
ALLSEL,ALL

NSEL, S, LOC, Y, TS, TS+TD        !杀死层单元
ESLN, S, 1, ALL
EKILL,ALL
ALLSEL,ALL
ESEL,S,LIVE
EPLOT
wpoffs,-stepdis,0,0   !偏移一个步长
*DO, i, 1, layer, 1
  *DO, j, 1, beadnum, 1
     *DO, t, StepTime, TotalTime, StepTime
         TIME, TotalTime*(j-1)+t+(i-1)*5*TotalTime                            !载荷步结束时间
                                          
        nsubst,4
        CSYS,4                               !激活4号坐标系
        CSDELE,ALL                           !删除局部坐标系
         wpoffs,stepdis,0,0                   !偏移一个步长
         CSWPLA,11,0,1,1,   
         CSWPLA,12,1,1,1,  
        csys,11
        NSEL, S, LOC, Y, -0.005, 0.005
        NSEL, R, LOC, X, 0, 0.004            !沿着焊接方向每次激活两个格
        NSEL, R, LOC, Z, -NNY,0
        ESLN, S, 1
        EALIVE,ALL
        ALLSEL, ALL
        csys,12                              !激活12号局部坐标系，首先删除要加载热源位置的对流边界，然后加载面热源--热流密度函数
        nsel,S,loc,x,0,0.0051
        NSEL,R, LOC, Z, 0
        sfdele,all,conv
        nsel,R,loc,x,0,0.005
        nplot
        csys,11
        sf,all,hflux,%XU%
        allsel,all
        csys,0
        ESEL,S,LIVE                          !激活单元的表面，如果包含左端面，指定为对流换热边条
        NSLE,S,1
        NSEL, R, LOC, Y, TS, TS+4*NNY
        NSEL, R, LOC, Z, 0
        SF, ALL, CONV, CoffConv, Temp
        ALLSEL,ALL
      
        ESEL,S,LIVE                           !激活单元的表面，如果包含右端面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,Y,TS, TS+4*NNY
        NSEL, R, LOC, Z, -LD
        SF, ALL, CONV, CoffConv, Temp
        ALLSEL,ALL        
        ESEL,S,LIVE                             !激活单元的表面，如果包含右端面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,Y,TS, TS+4*NNY
        NSEL, R, LOC, X, 0
        SF, ALL, CONV, CoffConv, Temp
        ALLSEL,ALL
        
        ESEL,S,LIVE                             !激活单元的表面，如果包含右端面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,Y,TS, TS+4*NNY
        NSEL, R, LOC, X, WD
        SF, ALL, CONV, CoffConv, Temp
        ALLSEL,ALL
        outres,all,3
        OUTPR,NSOL,LAST,
        SOLVE
        SAVE
   
        sfdele,all,hflux
        ESEL,S,LIVE                            !激活单元的上表面，指定为对流换热边条
        NSLE,S,1
        csys11
        NSEL, R, LOC, Z, 0
        SF, ALL, CONV, CoffConv, Temp
        ALLSEL,ALL
        ESEL,S,LIVE
        EPLOT
           *ENDDO
            wpoffs,-(LD+stepdis),10*NNX,0
               
     *ENDDO
      wpoffs,0,-WD,NNY
  *ENDDO  

ESEL,S,LIVE
EPLOT
!===============================================================================================!冷却
!==============================================================================!~100s
        *DO, k, 1, 2, 1
        TIME,TotalTime*beadnum*layer+50*k                                !载荷步结束时间
        NSUBST, 5
        KBC, 1
        SOLVE
        SAVE
        *ENDDO
!==============================================================================!~1000s
        *DO, k, 1, 9, 1
        TIME,TotalTime*beadnum*layer+100+100*k                                !载荷步结束时间
        NSUBST, 5
        KBC, 1
        SOLVE
        SAVE
        *ENDDO
!==============================================================================!~3000s
        *DO, k, 1, 10, 1
        TIME,TotalTime*beadnum*layer+1000+200*k                                !载荷步结束时间
        NSUBST, 5
        KBC, 1
        SOLVE
        SAVE
        *ENDDO
!==============================================================================!~10000s
        *DO, k, 1, 14, 1
        TIME,TotalTime*beadnum*layer+3000+500*k                                !载荷步结束时间
        NSUBST, 5
        KBC, 1
        SOLVE
        SAVE
        *ENDDO
save
FINISH

jx04319

问一下 热源形式除了高斯热源 在焊接中还有其他形式的吗