移动热源

hai00011

*SET,%_FNCNAME%(0,21,1), 0.0, 99, 0, 1, -2, 0, 0
变量应该是以字母开头吧,其他的看不懂啊

jun722

这样的程序都是读入函数生成的，没法看懂啊

xgmSS

什么是读入函数生成的？

jun722

先在函数编辑器中定义一个函数，然后保存，最后在读入这个函数，就会生成上面的很多程序

zhanglinghnsm

你是在所有单元上加载了随时间变化的热源函数，任何时刻所有的单元上都有热源，只是不同的时间有不同值，当然不会有移动的动画。

如果你是在每一个时间步都更新载荷（删除上一步的载荷，加这一步的载荷），那么就会看到动画。

另外，建议写程序时多谢注释，不然过段时间自己都看不懂。

zhanglinghnsm

我是做激光熔覆的，我们的方向很类似，以后多交流交流。
贴上我的程序，共同学习。

!上表面上没有对流换热边界条件
!单位制:米、秒、摄氏度
!
/CLEAR,START
/FILNAME,temp,0
/COM,ANSYS RELEASE 10.0    UP20050718       20:15:52    07/14/2007
/CONFIG, NRES, 5000
/PREP7  
!*
!===============================================================================================!指定单元
ET,1,SOLID70
!*  
!*  
!===============================================================================================!材料属性
MPTEMP,,,,,,,,  
MPTEMP,1,0
MPTEMP,2,1350
MPTEMP,3,1400
MPTEMP,4,1420   
MPTEMP,5,1440
MPTEMP,6,1490
MPTEMP,7,2000   
MPDATA,ENTH,1,,0
MPDATA,ENTH,1,,6.5e9
MPDATA,ENTH,1,,7.2e9
MPDATA,ENTH,1,,8e9
MPDATA,ENTH,1,,9e9  
MPDATA,ENTH,1,,9.9e9
MPDATA,ENTH,1,,12.84e9

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,800
MPTEMP,5,1440
MPTEMP,6,2000   
MPDATA,KXX,1,,14.7  
MPDATA,KXX,1,,18
MPDATA,KXX,1,,20.8  
MPDATA,KXX,1,,26.3  
MPDATA,KXX,1,,34.7
MPDATA,KXX,1,,51.5

MPTEMP,,
MPTEMP,1,0  
MPTEMP,2,200
MPTEMP,3,400
MPTEMP,4,800
MPTEMP,5,1200
MPTEMP,6,2000   
MPDATA,C,1,,494
MPDATA,C,1,,536
MPDATA,C,1,,569
MPDATA,C,1,,644
MPDATA,C,1,,669
MPDATA,C,1,,729

MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,7700  

!===============================================================================================!定义常量
Width_Base=0.025                                !宽度
Height_Base=0.02                                !基底高度
Length=0.05                                            !长度

Width_Clad=0.0015                                !宽度
Height_Deposition=0.0025                        !覆层高度

Layer=10                                        !层数
Height_Clad=Height_Deposition/layer
dt=0.00001                                        !小量

Temp=50                                                !环境温度
Init_Temp=200                                        !初始温度
Coff_Conv=30                                        !对流换热系数

!===============================================================================================!定义常量
Velocity=0.003                                !扫描速度
StepDis=0.001                                !每个载荷步位移

LaserPower=700                                !激光功率
Radius=0.0015                                !激光光斑半径
Area=3.14159265*(Radius**2)                !激光光斑面积
Factor=0.25                                !吸收因子


StepTime=StepDis/Velocity                !每个载荷步时间
TotalTime=(Length+Radius*2)/Velocity        !载荷持续时间(扫描一层)
StepNum=(Length+Radius*2)/StepDis        !载荷步数
!===============================================================================================!建模

BLOCK,0,Length,0,Height_Deposition,0,Width_Clad,                !覆层
BLOCK,0,Length,0,-Height_Base,0,Width_Clad,                           !基底
BLOCK,0,Length,0,-Height_Base,Width_Clad,Width_Base  
VGLUE,ALL
NUMCMP,ALL

!===============================================================================================!划分网格

LESIZE,18,,,15,0.125
LESIZE,22,,,15,8
LESIZE,24,,,15,8
     
LESIZE,26,,,15,8
LESIZE,28,,,15,8


LESIZE,20,,,15,8
LESIZE,21,,,15,8
LESIZE,23,,,15,8
     
LESIZE,25,,,15,8
LESIZE,27,,,15,8

LSEL, S, LOC, Y, dt, Height_Deposition-dt,                        !覆层高度方向的单元数目
LESIZE, ALL, , , Layer,

LSEL, S, LOC, X, dt, Length-dt,                                !长度方向的单元数目
LESIZE, ALL, , , Length/StepDis,

LSEL, S, LOC, Z, dt, Width_Clad-dt,                        !宽度方向的单元数目
LESIZE, ALL, , , 5,

VSEL,ALL                                                !网格划分
TYPE,1
MAT,1
MSHAPE,0,3D
MSHKEY,1
VMESH,ALL


!===============================================================================================!基底边界条件、初始条件
NSEL, S, LOC, Y, -Height_Base, 0                        !基底初始温度
IC,ALL,TEMP,Init_Temp
ALLSEL,ALL

NSEL, S, LOC, Y, -Height_Base, 0                        !基底侧面，对流换热边界条件
NSEL, R, LOC, Z, Width_Base
SF, ALL, CONV, Coff_Conv, Temp
ALLSEL,ALL

NSEL, S, LOC, Y, -Height_Base, 0                        !基底左端面，对流换热边界条件
NSEL, R, LOC, X, 0
SF, ALL, CONV, Coff_Conv, Temp
ALLSEL,ALL

NSEL, S, LOC, Y, -Height_Base, 0                        !基底右端面，对流换热边界条件
NSEL, R, LOC, X, Length
SF, ALL, CONV, Coff_Conv, Temp
ALLSEL,ALL

NSEL, S, LOC, Y, 0                                               !基底上表面，对流换热边界条件
NSEL, R, LOC, Z, Width_Clad, Width_Base
SF, ALL, CONV, Coff_Conv, Temp
ALLSEL,ALL

FINISH

/SOLU

!===============================================================================================!瞬态分析参数设置
ANTYPE,4                                !分析类型：瞬态
!*  
TRNOPT,FULL                                 !求解方法：完全的N-R方法
LUMPM,0                                 !不使用集中质量矩阵
!*  
!NSUBST, 1                                !指定载荷子步数目
DELTIM,0.01,0.001,0.05                         !载荷子步(默认子步时间步长、最小、最大)——载荷步为0.333
CNVTOL,HEAT, ,0.01,1,0.000001,           !收敛准则：控制热流
OUTRES,NSOL,LAST                        !结果输出：所有

!===============================================================================================!杀死单元
NSEL, S, LOC, Y, 0, Height_Deposition        !杀死熔覆层单元
ESLN, S, 1, ALL
EKILL,ALL
ALLSEL,ALL

ESEL,S,LIVE
EPLOT

ESEL,S,LIVE                                        !激活单元的上表面，指定为对流换热边条
NSLE,S,1
NSEL,R,LOC,Y,HeightBase
NSEL,R,LOC,Z,0,Radius
SF,ALL,CONV,Coff_Conv,Temp
ALLSEL,ALL
!===============================================================================================!预热
*DO, i, 1, 2

m=mod(i,2)

*IF,m,EQ,1,THEN                                        !如果为奇数层，向右扫描

        *DO, k, 1, StepNum, 1

        TIME,TotalTime*(i-1)+StepTime*k                !载荷步结束时间
        KBC, 1                                        !载荷步内载荷随时间分布：常数

        LeftX=StepDis*(k-1)
        RightX=StepDis*k


        NSEL, S, LOC, Y, 0
        ESEL, S, LIVE
        ESLN, R, 0
        NSLE, S, 1
        NSEL, R, LOC, x, RightX-2*Radius, RightX
        NSEL, R, LOC, Z, 0, Radius                        !激活单元的上表面，加热流密度
        ESLN, S, 1
        SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,
        ALLSEL, ALL

        SOLVE

        SFEDELE,ALL,4,HFLUX                                !删除热流密度载荷
        ALLSEL, ALL

        ESEL,S,LIVE                                        !激活单元的上表面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,Y,0
        NSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDis
        SF,ALL,CONV,Coff_Conv,Temp
        ALLSEL,ALL
         
  ESEL,S,LIVE
        EPLOT

        *ENDDO

*ELSE                                                !如果为偶数层，向左扫描

        *DO, k, 1, StepNum, 1

        TIME,TotalTime*(i-1)+StepTime*k                !载荷步结束时间
        KBC, 1                                        !载荷步内载荷随时间分布：常数

        LeftX=Length-StepDis*k
        RightX=Length-StepDis*(k-1)


        NSEL, S, LOC, Y, 0
        ESEL, S, LIVE
        ESLN, R, 0
        NSLE, S, 1
        NSEL, R, LOC, x, RightX-2*Radius, RightX
        NSEL, R, LOC, Z, 0, Radius                        !激活单元的上表面，加热流密度
        ESLN, S, 1
        SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,
        ALLSEL, ALL

        SOLVE

        SFEDELE,ALL,4,HFLUX                                !激活单元的上表面，删除载荷
        ALLSEL, ALL

        ESEL,S,LIVE                                        !激活单元的上表面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,Y,0
        NSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDis
        SF,ALL,CONV,Coff_Conv,Temp
        ALLSEL,ALL
       
        ESEL,S,LIVE
        EPLOT

        *ENDDO

*ENDIF

*ENDDO

ESEL,S,LIVE
EPLOT
!===============================================================================================!熔覆

*DO, i, 1, Layer, 1

m=mod(i,2)

*IF,m,EQ,1,THEN                                        !如果为奇数层，向右扫描

        *DO, k, 1, StepNum, 1

        TIME,TotalTime*(i+1)+StepTime*k                !载荷步结束时间
        KBC, 1                                        !载荷步内载荷随时间分布：常数

        LeftX=StepDis*(k-1)
        RightX=StepDis*k

        NSEL, S, LOC, Y, Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i
        NSEL, R, LOC, x, LeftX, RightX
        ESLN, S, 1
        EALIVE,ALL
        ALLSEL, ALL

        NSEL, S, LOC, Y, Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i
        NSEL, R, LOC, x, RightX-2*Radius, RightX
        NSEL, R, LOC, Z, 0, Radius                        !激活单元的上表面，加热流密度
        ESLN, S, 1
        SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,
        ALLSEL, ALL

        ESEL,S,LIVE                                        !激活单元的表面，如果包含左端面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,X,0
        NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i
        SF,ALL,CONV,Coff_Conv,Temp
        ALLSEL,ALL

        ESEL,S,LIVE                                        !激活单元的表面，如果包含右端面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,X,Length
        NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i
        SF,ALL,CONV,Coff_Conv,Temp
        ALLSEL,ALL

        ESEL,S,LIVE                                        !激活单元的侧面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,Z,Width_Clad
        NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i
        SF,ALL,CONV,Coff_Conv,Temp
        ALLSEL,ALL

        SOLVE

        SFEDELE,ALL,4,HFLUX                                !删除热流密度载荷
        ALLSEL, ALL
       
        ESEL,S,LIVE                                        !激活单元的上表面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,Y,Height_Deposition/Layer*i
        NSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDis
        SF,ALL,CONV,Coff_Conv,Temp
        ALLSEL,ALL

        NSEL,S,LOC,Y,Height_Deposition/Layer*(i-1)        !激活单元的下表面，删除对流换热边条
        ESEL, S, LIVE
        ESLN, R, 0
        NSLE, S, 1
        NSEL, R, LOC, x, LeftX-2*StepDis,RightX-2*StepDis
        NSEL, R, LOC, Z, 0, Radius
        NSEL, U, LOC, Y, Height_Deposition/Layer*i
  ESLN, S, 1
        SFDELE, ALL, CONV
        ALLSEL,ALL
         
  ESEL,S,LIVE
        EPLOT

        *ENDDO

*ELSE                                                !如果为偶数层，向左扫描

        *DO, k, 1, StepNum, 1

        TIME,TotalTime*(i+1)+StepTime*k                !载荷步结束时间
        KBC, 1                                        !载荷步内载荷随时间分布：常数

        LeftX=Length-StepDis*k
        RightX=Length-StepDis*(k-1)

        NSEL, S, LOC, Y, Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i
        NSEL, R, LOC, x, LeftX, RightX
        ESLN, S, 1
        EALIVE,ALL
        ALLSEL, ALL

        NSEL, S, LOC, Y, Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i
        NSEL, R, LOC, x, LeftX, LeftX+2*Radius
        NSEL, R, LOC, Z, 0, Radius                        !激活单元的上表面，加热流密度
        ESLN, S, 1
        SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,
        ALLSEL, ALL


        ESEL,S,LIVE                                        !激活单元的表面，如果包含左端面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,X,0
        NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i
        SF,ALL,CONV,Coff_Conv,Temp
        ALLSEL,ALL

        ESEL,S,LIVE                                        !激活单元的表面，如果包含右端面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,X,Length
        NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i
        SF,ALL,CONV,Coff_Conv,Temp
        ALLSEL,ALL

        ESEL,S,LIVE                                        !激活单元的侧面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,Z,Width_Clad
        NSEL,R,LOC,Y,Height_Deposition/Layer*(i-1), Height_Deposition/Layer*i
        SF,ALL,CONV,Coff_Conv,Temp
        ALLSEL,ALL

        SOLVE

        SFEDELE,ALL,4,HFLUX                                !激活单元的上表面，删除载荷
        ALLSEL, ALL

        ESEL,S,LIVE                                        !激活单元的上表面，指定为对流换热边条
        NSLE,S,1
        NSEL,R,LOC,Y,Height_Deposition/Layer*i
        NSEL, R, LOC, x, LeftX+2*StepDis, RightX+2*StepDis
        SF,ALL,CONV,Coff_Conv,Temp
        ALLSEL,ALL

        NSEL,S,LOC,Y,Height_Deposition/Layer*(i-1)        !激活单元的下表面，删除对流换热边条
        ESEL, S, LIVE
        ESLN, R, 0
        NSLE, S, 1
        NSEL, R, LOC, x, LeftX+2*StepDis,RightX+2*StepDis
        NSEL, R, LOC, Z, 0, Radius
        NSEL, U, LOC, Y, Height_Deposition/Layer*i
  ESLN, S, 1
        SFDELE, ALL, CONV
        ALLSEL,ALL
       
        ESEL,S,LIVE
        EPLOT

        *ENDDO

*ENDIF

*ENDDO

ESEL,S,LIVE
EPLOT

acailaogu

靠，，楼上的程序不错啊！！看来是个高手！！
以后多多请教了！

wxpand2004

我也是做熔覆的，你的程序先学习了
希望以后可以多多指点！
7# zhanglinghnsm

laokaner

学习学习，不错的很啊

zhjg216

很好，学习一下

huake1

老天，不活了，我也是这个，还没弄你那么多，楼主加油。。。