..:: ANSYS专区-激光,耦合,热应力,烧结-激光烧结热应力耦合分析，{Simwe仿真论坛(forum.simwe.com)，CAE/CAD/CAM/，FEA/FEM/有限元分析论坛-(手机验证注册)}

feilong_065191 发表于 2010-5-18 09:58:51

激光烧结热应力耦合分析

不知道怎么回事，我在调用温度场的时候老是出错
*** WARNING ***
Element convections are not valid for element 995.
The SFE command is ignored.
这样的WARNING一大堆，但是如果我单个调用温度场的结果时却没有这样的WARNING，当然我没有一个个调用，只不过用LDREAD,TEMP,,,,,SLS_temp,RTH调用默认的。
温度场我粉层用的是solid70，基板solid90
以下是我的apdl，跪求高手解答，谢谢：
/CONFIG,NRES,1e6
/units,si!m, kg, s, K
/filname,SLS_temp,0
/title,SLS thermal analysis
!选择分析问题的类型：结构、热
!--------------Parameters---------------!
*SET,Block_x,88e-4 !几何模型尺寸
*SET,Block_y,88e-4
*SET,Block_z,8e-4
*SET,numMesh,4e-4 !划分网格的大小
*SET,Tuniform,300
*SET,hc,6
*SET,Tamb,300
!---------------------------------------!
/PREP7
!1定义单元类型
ET,1,SOLID70
ET,2,SURF152
DOF,TEMP
KEYOPT,2,4,1
KEYOPT,2,8,3
ET,3,SOLID90
!2实常数
!3材料热物性数据
!定义材料1属性(304L块材)
!!!!!!!!!!!!!!!!!!!!材料一(粉末304L不锈钢)的热物性参数!!!!!!!!!!!!!!!!
MPTEMP,1,300,473,573,673,773,873
MPTEMP,7,973,1073,1170,1370,1500
MPDATA,DENS,1,1,7694,7890,7804,7716,7631,7539
MPDATA,DENS,1,7,7492,7444,7396,7301,7273!(Kg/m3)
MPDATA,KXX,1,1,12.6,13.1,14.0,15.6,17.2,19.2
MPDATA,KXX,1,7,20.7,22.8,24.0,26.3,31.7
MPDATA,C,1,1,446,487,522,564,580,613
MPDATA,C,1,7,622,635,646,666,682
MPDATA,ENTH,1,1,0.96E9,1.60E9,2.37E9,3.18E9,4.01E9,4.87E9
MPDATA,ENTH,1,7,5.33E9,5.77E9,6.24E9,6.50E9,6.65E9
!!!!!!!!!!!!!!!!!!!!材料二(基板45钢)的热物性参数!!!!!!!!!!!!!!!!
MP,ALPX,2,1.06E-5
MP,DENS,2,7870
MPTEMP,1,200,300,400,500,600,800
MPTEMP,7,900,1000,1100,1200,1300
MPDATA,KXX,2,1,62.1,55.5,49.2,43,38.2,29.4
MPDATA,KXX,2,7,28.2,29.3,32.2,32.2,32.2
MPDATA,C,2,1,518.8,552.3,602.5,661.1,753.1,832.6
MPDATA,C,2,7,656.9,569.0,606.7,640.2,736.4
MPDATA,ENTH,2,1,6.4E8,9.9E8,1.4E9,1.8E9,2.3E9,3.25E9
MPDATA,ENTH,2,7,3.8E9,4.3E9,4.9E9,5.5E9,6.1E9
CSYS,0
BLOCK,0,Block_x,0,Block_y,0,Block_z
BLOCK,-12e-4,100e-4,-12e-4,100e-4,0,-8e-4 !建立模型
VGLUE,1,2!合并
NUMCMP,ALL
LESIZE,1,numMesh
LESIZE,2,numMesh
LESIZE,3,numMesh
LESIZE,4,numMesh
LESIZE,5,numMesh
LESIZE,6,numMesh
LESIZE,7,numMesh
LESIZE,8,numMesh
LESIZE,9,numMesh
LESIZE,10,numMesh
LESIZE,11,numMesh
LESIZE,12,numMesh
MAT,1
VMESH,1
ASEL,S,LOC,Z,0.0008
TYPE,2
AMESH,ALL
ESEL,S,TYPE,,2
SF,ALL,CONV,hc,Tamb
/PSF,CONV,2 !划分粉层网格
MAT,2
type,3
ESIZE,0.0004
MSHAPE,1,3d
VMESH,2
ALLSEL,ALL                !划分基板网格
FINISH
!6其他
!*************************************************************************************************!
!开始，进入处理器，定义下列属性－
!1对流换热
!2模拟激光加热（指定某一单元温度值），移动策略（根据扫描方式确定）
!*************************************************************************************************!
/SOLU
ANTYPE,4
TRNOPT,FULL
LUMPM,0
NROPT,FULL,,ON
TUNIF,Tuniform
TREF,Tuniform
*SET,InitX,3*numMesh
*SET,InitY,3*numMesh
*SET,numZ,Block_z
*SET,NumTime,0
*SET,NumLoadStep,0
V=0.01
TIMC=numMesh/V
hf_cent=336000000
!1对流换热
ALLSEL,ALL
ASEL,S,,,3,6,1
NSLA,S,1
SF,ALL,CONV,hc,Tamb
/PSF,CONV,2
SFTRAN
NROPT,FULL
TIMINT,ON
!2模拟激光加热，移动热源"S"
*SET,Flag,1
*DO,i,1,16,1
*DO,j,1,16,1
!*定义载荷步选项
*SET,NumTime,NumTime+TIMC
Time,NumTime
AUTOTS,ON
            NSUBST,1,,,ON
            KBC,1
            !*定义单元热流密度
NSEL,S,LOC,Z,numZ
NSEL,R,LOC,X,InitX+numMesh
NSEL,R,LOC,Y,InitY-numMesh
ESEL,S,TYPE,,1
            ESLN,R,0
            ESEL,R,TYPE,,1
            *GET,meax,ELEM,ALL,NUM,MAX
            *SET,SEL_ELEM,meax
            SFE,SEL_ELEM-1,6,HFLUX,,hf_cent
            SFE,SEL_ELEM-22,6,HFLUX,,hf_cent*0.6128
            SFE,SEL_ELEM-23,6,HFLUX,,hf_cent*0.7828
            SFE,SEL_ELEM,6,HFLUX,,hf_cent*0.7828
            NSEL,S,LOC,Z,numZ
            NSEL,R,LOC,X,InitX+numMesh
            NSEL,R,LOC,Y,InitY+numMesh
            ESEL,S,TYPE,,1
            ESLN,R,0
            ESEL,R,TYPE,,1
            *GET,emax,ELEM,ALL,NUM,MAX
            *SET,SEL_ELE,emax
            SFE,SEL_ELE,6,HFLUX,,hf_cent*0.6128
            SFE,SEL_ELE-1,6,HFLUX,,hf_cent*0.7828
            SFE,SEL_ELE-22,6,HFLUX,,hf_cent*0.7828
            SFE,SEL_ELE-23,6,HFLUX,,hf_cent
            ALLSEL,ALL
            NSEL,S,LOC,Z,numZ
            NSEL,R,LOC,X,InitX-numMesh
            NSEL,R,LOC,Y,InitY+numMesh
            ESLN,R,0
            ESEL,R,TYPE,,1
            *GET,eamx,ELEM,ALL,NUM,MAX
            *SET,SEL_EL,eamx
            SFE,SEL_EL,6,HFLUX,,hf_cent*0.7828
            SFE,SEL_EL-1,6,HFLUX,,hf_cent*0.6128
            SFE,SEL_EL-22,6,HFLUX,,hf_cent
            SFE,SEL_EL-23,6,HFLUX,,hf_cent*0.7828
            ALLSEL,ALL
            NSEL,S,LOC,Z,numZ
            NSEL,R,LOC,X,InitX-numMesh
            NSEL,R,LOC,Y,InitY-numMesh
            ESLN,R,0
            ESEL,R,TYPE,,1
            *GET,eaxm,ELEM,ALL,NUM,MAX
            *SET,SEL_E,eaxm
            ESEL,S,TYPE,,1
            SFE,SEL_E,6,HFLUX,,hf_cent
            SFE,SEL_E-22,6,HFLUX,,hf_cent*0.7828
            SFE,SEL_E-1,6,HFLUX,,hf_cent*0.7828
            SFE,SEL_E-23,6,HFLUX,,hf_cent*0.6128
            /PSF,HFLUX,,2
            !*写载荷步
ALLSEL,ALL
*SET,NumLoadStep,NumLoadStep+1
LSWRITE,NumLoadStep
!*清理上一载荷步的内容
            SFEDELE,ALL,6,HFLUX
      *SET,InitX,InitX+Flag*numMesh
*ENDDO
*SET,InitX,InitX-Flag*numMesh
*SET,Flag,Flag*(-1)
*SET,InitY,InitY+numMesh
*ENDDO
ALLSEL,ALL
NEQIT,80
LSSOLVE,1,NumLoadStep,1
FINISH
/POST1
PLNSOL,TEMP,0
/UI,COPY,SAVE,JPEG,GRAPHY,COLOR,NORM,LANDSCAPE,NO,100 !自动保存图片
/POST26
nsol,2,992,temp,,temp992 !取出节点992的温度值并赋给变量号2
/gropt,axnsc,1.5,!缩放图标字体的大小
plvar,2!图形显示变量2
/IMAGE,SAVE,sdirection,BMP
SAVE
FINISH

/COM,STRUCTURAL!指定为结构分析范畴
ii=0
/PREP7
ETCHG,TTS
MPTEMP !清除原来的温度场
MPTEMP,1,300,473,573,673,773,873
MPTEMP,7,973,1073,1170,1370,1500
MPDATA,EX,1,1,2.00E11,1.95E11,1.90E11,1.85E11,1.73E11,1.65E11
MPDATA,EX,1,7,1.55E11,1.44E11,1.03E11,0.71E11,0.11E11
MPDATA,PRXY,1,1,0.3,0.3,0.3,0.3,0.3,0.3
MPDATA,PRXY,1,7,0.3,0.3,0.3,0.3,0.3
MPDATA,ALPX,1,1,14.0E-6,14.7E-6,15.9E-6,16.8E-6,17.7E-6,18.3E-6
MPDATA,ALPX,1,7,18.7E-6,19.0E-6,19.2E-6,19.5E-6,20.0E-6!(1/K)
TB,BKIN,1,6
TBTEMP,300
TBDATA,1,278E6,1.96E9!Yield (Pa); Tangent (Pa)
TBTEMP,473
TBDATA,1,254E6,1.92E9
TBTEMP,573
TBDATA,1,208E6,1.90E9
TBTEMP,673
TBDATA,1,182E6,1.82E9
TBTEMP,773
TBDATA,1,159E6,1.69E9
TBTEMP,873
TBDATA,1,140E6,1.61E9
MPTEMP
MPTEMP,1,200,300,400,500,600,800
MPTEMP,7,900,1000,1100,1200,1300
MPDATA,EX,2,1,1.94E11,1.87E11,1.79E11,1.7E11,1.57E11,1.29E11
MPDATA,EX,2,7,1.11E11,0.9E11,0.69E11,0.477E11,0.33E11
MPDATA,PRXY,2,1,0.3,0.3,0.3,0.3,0.3,0.3
MPDATA,PRXY,2,7,0.3,0.3,0.3,0.3,0.3
FINISH
/SOLU
NLGEOM,ON
DA,7,UZ,0
DA,8,UY,0
DA,10,UX,0
TREF,300
*DO,ii,1,256
LDREAD,TEMP,ii,,,,SLS_temp,RTH
TIME,ii*0.04
NSUBST,1,,,ON
*ENDDO
LSSOLVE,1,256,1
FINISH
/POST1
PLNSOL,S,EQV
/UI,COPY,SAVE,JPEG,GRAPHY,COLOR,NORM,LANDSCAPE,NO,100
PLDISP,1
/UI,COPY,SAVE,JPEG,GRAPHY,COLOR,NORM,LANDSCAPE,NO,100
SAVE
FINISH

feilong_065191 发表于 2010-5-18 14:56:27

我怪小弟我是新手，本来应该在标题上写清楚这是个求助帖，哎···
期待高手回复中···
不知道怎么回事，调用温度场结果的时候会出错，若用*DO循环调用全部*.rth结果时，就会不停的弹出警告窗口，*** WARNING ***
Element convections are not valid for element 995.
The SFE command is ignored.
当然不止这一个单元，好多单元···
跪求高手解答

feilong_065191 发表于 2010-5-23 18:17:49

好歹有个人回个撒，不要让我的希望变成失望啊

lixiaoyu0323 发表于 2010-9-25 16:36:04

LZ,练到你这程度得要修炼多久啊

htcpro 发表于 2010-11-2 15:24:26

没看懂。。。。。。。。。

373684251 发表于 2010-11-2 16:08:14

不是我做的方向，不懂。

htcpro 发表于 2011-1-4 15:05:17

如何处理材料熔化后凝固的问题

jiangyiyoung 发表于 2011-4-20 13:33:51

2# feilong_065191

lz，你的温度场怎么和我一样，有4000-5000度？这样的温度场和不合理啊？能否交流一下？

feilong_065191 发表于 2011-12-5 20:44:40

大学毕业一年半，如今看到自己当年的程序，感慨万千，怎么当初我就那么有才呢，现在叫我看都难看懂，真是悲剧

页: [1]

Simwe仿真论坛(forum.simwe.com)，CAE/CAD/CAM/，FEA/FEM/有限元分析论坛---(手机验证注册)'s Archiver

激光烧结热应力耦合分析