【讨论】请高手帮忙！焊接温度场的分析！

aiboderen

平板堆焊
建模采用了完整模型的一半，焊缝按半球形，建模后是四分之一球形，热流服从高斯分布。采用plane55单元。
有很多问题，请高手帮忙！
FINI
/CLE
/PREP7   
!*   
ET,1,PLANE55
!*   
!*   
MPTEMP,,,,,,,,   
MPTEMP,1,0   
MPDATA,KXX,1,,20
MPTEMP,,,,,,,,   
MPTEMP,1,0   
MPDATA,C,1,,0.5  
MPTEMP,,,,,,,,   
MPTEMP,1,0   
MPDATA,KXX,2,,30
MPTEMP,,,,,,,,   
MPTEMP,1,0   
MPDATA,C,2,,0.6  
*AFUN,DEG   
CYL4,0,0,0,0,6,90
RECTNG,0,6,0,-20,   
RECTNG,6,100,0,-20,  
FLST,2,3,5,ORDE,2   
FITEM,2,1   
FITEM,2,-3   
AGLUE,P51X   
ALLSEL,ALL
ASEL,S , , ,       1
AATT,       1, ,   1,       0,  
LESIZE,1, , ,50, , , , ,0  
LESIZE,2, , ,30, , , , ,0
LESIZE,3, , ,30, , , , ,0
MSHAPE,1,2D  
MSHKEY,0
AMESH,1
CM,_Y,AREA   
ASEL, , , ,       4  
CM,_Y1,AREA  
CMSEL,S,_Y   
!*   
CMSEL,S,_Y1  
AATT,       2, ,   1,       0,   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
!*   
/AUTO,1  
/REP,FAST   
CM,_Y,AREA   
ASEL, , , ,       4  
CM,_Y1,AREA  
CMSEL,S,_Y   
!*   
CMSEL,S,_Y1  
AATT,       2, ,   1,       0,   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
!*   
FLST,5,2,4,ORDE,2   
FITEM,5,3   
FITEM,5,-4   
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
!*   
LESIZE,_Y1, , ,30, , , , ,1  
!*   
FLST,5,2,4,ORDE,2   
FITEM,5,12   
FITEM,5,-13  
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
!*   
FLST,5,2,4,ORDE,2   
FITEM,5,12   
FITEM,5,-13  
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
!*   
LESIZE,_Y1, , ,30, , , , ,1  
!*   
CM,_Y,AREA   
ASEL, , , ,       4  
CM,_Y1,AREA  
CHKMSH,'AREA'   
CMSEL,S,_Y   
!*   
AMESH,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2   
!*   
GPLOT   
CM,_Y,AREA   
ASEL, , , ,       5  
CM,_Y1,AREA  
CMSEL,S,_Y   
!*   
CMSEL,S,_Y1  
AATT,       2, ,   1,       0,   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
!*   
FLST,5,2,4,ORDE,2   
FITEM,5,14   
FITEM,5,-15  
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
!*   
LESIZE,_Y1, , ,30, , , , ,1  
!*   
FLST,5,2,4,ORDE,2   
FITEM,5,9   
FITEM,5,12   
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
!*   
LESIZE,_Y1, , ,30, , , , ,1  
!*   
CM,_Y,AREA   
ASEL, , , ,       5  
CM,_Y1,AREA  
CHKMSH,'AREA'   
CMSEL,S,_Y   
!*   
AMESH,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2   
!*   
FLST,2,2,5,ORDE,2   
FITEM,2,4   
FITEM,2,-5   
ACLEAR,P51X  
GPLOT   
CM,_Y,AREA   
ASEL, , , ,       4  
CM,_Y1,AREA  
CMSEL,S,_Y   
!*   
CMSEL,S,_Y1  
AATT,       2, ,   1,       0,   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
!*   
FLST,5,2,4,ORDE,2   
FITEM,5,3   
FITEM,5,-4   
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
!*   
LESIZE,_Y1, , ,30, , , , ,1  
!*   
FLST,5,2,4,ORDE,2   
FITEM,5,12   
FITEM,5,-13  
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
!*   
LESIZE,_Y1, , ,30, , , , ,1  
!*   
CM,_Y,AREA   
ASEL, , , ,       4  
CM,_Y1,AREA  
CHKMSH,'AREA'   
CMSEL,S,_Y   
!*   
AMESH,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2   
!*   
GPLOT   
CM,_Y,AREA   
ASEL, , , ,       5  
CM,_Y1,AREA  
CMSEL,S,_Y   
!*   
CMSEL,S,_Y1  
AATT,       2, ,   1,       0,   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
!*   
FLST,5,2,4,ORDE,2   
FITEM,5,9   
FITEM,5,12   
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
!*   
LESIZE,_Y1, , ,30, , , , ,1  
!*   
FLST,5,2,4,ORDE,2   
FITEM,5,14   
FITEM,5,-15  
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
!*   
LESIZE,_Y1, , ,30, , , , ,1  
!*   
CM,_Y,AREA   
ASEL, , , ,       5  
CM,_Y1,AREA  
CHKMSH,'AREA'   
CMSEL,S,_Y   
!*   
AMESH,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2   
!*   
  
ALLSEL,ALL  
FINISH   
/SOLU   
!*   
ANTYPE,4
!*   
TRNOPT,FULL  
LUMPM,0  
CSYS,1   
wpstyle,0.05,0.1,-1,1,0.003,1,2,,5  
ALLSEL,ALL
*DO,I,1,30,1
NSEL,S,LOC,X,0.2*(I-1) ,0.2*I  
NSEL,R,LOC,Y,0,90   
ESLN,R  
SFE,ALL,,HFLUX, ,900*exp(-1.3*(0.2*(I-1)*(I-1))), , ,  
ALLSEL,ALL
*ENDDO
ALLSEL,ALL
wpstyle,0.05,0.1,-1,1,0.003,0,2,,5   
OUTPR,BASIC,ALL,
TIME,60
SOLVE

ainite

也许别人很忙，没时间把命令流跑一遍，你干脆把出现的问题描述一下，如果是经常遇到的典型问题，就直接回复了。

aiboderen

问题是焊缝有彩云显示，而平板没有任何变化！
这是怎么回事？已经粘在一起了呀！

jnjl

１、你看看单位是否统一，比如导热系数的单位是否和后面其他物理量的单位一致？可能是你的导热系数太小了，所以热传导的作用就不明显了；
２、你把加载的热源加大十倍或者几十倍，看看结果是否还是平板没有变化。
如果还是不行，大家再帮你想想其他的原因。

ainite

后处理中分开选取焊缝区域或平板区域，再PLOT看看。

aiboderen

多谢二位了
我试一下

aiboderen

可能我的一些参数，比如热导率、密度等没有设好，我算的时候都是自己估计的数字！
那么，请问，关于这些参数，DENS、KXX、C都是在哪里查呀？
如果关于温度变化，变化的值又怎么算？
深表感谢！

aiboderen

再问：
     关于焊接后热过程中氢的扩散，它加的边界条件同温度场一样！也就是说可以按照设定温度场分析那样设置边界条件。
    因为氢的扩散区域非常小，所以，按恒温扩散考虑。假设那个区域温度是1000℃,区域大小已知。那么，怎么才能把氢的扩散彩云图算出来？

jnjl

用自己估计的数字当然不可以了，去图书馆找一些专业的书籍，金属材料手册之类的，应该都有介绍的。
氢的扩散彩云图？你是说后处理吧？如果和温度场差不多的话，就用plot ctrls->nodal solution 显示，想看具体节点的数值，用ｌｉｓｔ命令，你试试吧！

200802172

xuexi le     dd