24小时热门版块排行榜    

查看: 635  |  回复: 2

wuhaitao9029

金虫 (小有名气)

[求助] 使用生死单元的喷涂模拟中死单元的数值为什么会变化呢?

做的等离子喷涂模拟,三种材料,各自赋有初始温度条件,我时间步按照1s设置,感觉不收敛,我就给了一个初始的0.001s,但是我用生死单元的时候,就会出现,还没激活的死单元的温度条件发生了变化。求大家帮忙看下,是哪里出现的问题。
/PREP7  
!*  
ET,1,PLANE13
!*  
KEYOPT,1,1,4
KEYOPT,1,2,0
KEYOPT,1,3,1
KEYOPT,1,4,0
KEYOPT,1,5,0
!*  
MP,DENS,1,2715  
MP,KXX,1,221.5  
MP,C,1,921  
MP,EX,1,0.689E11
MP,PRXY,1,0.33  
MP,ALPX,1,23.6E-6   
MP,DENS,2,7320  
MP,KXX,2,4.3
MP,C,2,501  
MP,EX,2,2.25E11
MP,PRXY,2,0.30  
MP,ALPX,2,14E-6
MP,DENS,3,6037  
MP,KXX,3,1.5
MP,C,3,500  
MP,EX,3,0.53E11
MP,PRXY,3,0.25  
MP,ALPX,3,7.2E-6
*SET,Sub_length,0.025   
*SET,Sub_height,0.002   
*SET,Coat_height,0.0005
*SET,Metallic_height,0.0001
*SET,YSZ_height, Coat_height- Metallic_height   
*SET,Coat_layer,10  
*SET,Sub_layer,5
*SET,Single_height, Coat_height /Coat_layer
*SET,Ta,300
*SET,Ha,8   
*SET,Tg_metallic,723
*SET,Tg_YSZ,748
*SET,Hg,320
*SET,dt,1e-7
*SET,Inittemp_sub,300   
*SET,Inittemp_metallic,1673
*SET,Inittemp_YSZ,2478  
*SET,Single_time,10
K,1,0,0
K,2,25e-3,0
K,3,25e-3,2e-3  
K,4,0,2e-3  
K,5,0,2.1e-3
K,6,25e-3,2.1e-3
K,7,25e-3,2.5e-3
K,8,0,2.5e-3
A,1,2,3,4   
A,3,4,5,6   
A,5,6,7,8   
Aglue,all   
Numcmp,all  
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
APLOT   
/PNUM,KP,0  
/PNUM,LINE,0
/PNUM,AREA,1
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/PNUM,DOMA,0
/NUMBER,0   
!*  
/PNUM,ELEM,0
/REPLOT
!*  
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
/DIST,1,0.924021086472,1
/REP,FAST   
FINISH  
/SOL
FINISH  
/PREP7  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       1, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       2
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       3
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       3, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
LESIZE,1,,,250  
LESIZE,3,,,250  
LESIZE,6,,,250  
LESIZE,9,,,250  
LESIZE,2,,,5
LESIZE,4,,,5
LESIZE,7,,,2
LESIZE,5,,,2
LESIZE,8,,,8
LESIZE,10,,,8   
MSHAPE,0,2D
MSHKEY,1
!*  
FLST,5,3,5,ORDE,2   
FITEM,5,1   
FITEM,5,-3  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
/FOC,1,-0.3,,,1
/REP,FAST   
NSEL,S,LOC,Y,0  
SF,ALL,CONV,Ha,Ta   
ALLSEL,ALL  
NSEL,S,LOC,Y,0,Sub_height   
IC,ALL,TEMP, Inittemp_sub   
ALLSEL,ALL  
NSEL,S,LOC,Y,Sub_height+dt,Sub_height+Metallic_height   
IC,ALL,TEMP, Inittemp_metallic  
ALLSEL,ALL  
NSEL,S,LOC,Y,Sub_height+ Metallic_height+dt,Sub_height+Coat_height  
IC,ALL,TEMP, Inittemp_YSZ   
ALLSEL,ALL  
FINISH  
/SOLU   
!*  
ANTYPE,4
!*  
TRNOPT,FULL
LUMPM,0
!*  
NSUBST,500,1000,200
OUTRES,ERASE
OUTRES,ALL,50   
NSEL,S,LOC,Y, Sub_height, Sub_height +Coat_height   
ESLN,S,1
EKILL,ALL   
ESEL,S,LIVE
EPLOT   
*DO,i, 1,Coat_layer,1   
TIME,Single_time*(i-1)+0.001
NSEL,S,LOC,Y, Sub_height+Single_height*(i-1),Sub_height +Single_height*i
ESLN,S,1
EALIVE,ALL  
ALLSEL,ALL  
NSEL,S,LOC,Y, Sub_height +Single_height*i   
SF,ALL,CONV,Hg,Tg_metallic  
ALLSEL,ALL  
NSEL,S,LOC,Y, Sub_height +Single_height*(i-1)   
SFDELE,ALL,CONV
ALLSEL,ALL  
SOLVE   
TIME,Single_time*i  
SOLVE   
ESEL,S,LIVE
EPLOT   
*ENDDO  
ESEL,S,LIVE
EPLOT   
FINISH
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

htbbzzg

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
在 ANSYS 中,死单元实际上也参加计算,只是 ANSYS 把死单元的主要材料性能设置的很小,以基本不影响周围的单元 (反之则不然),因此,分析后,死单元的自由度 (热 - 温度;结构 - 位移,等) 也是会改变的。
2楼2013-11-20 18:48:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuhaitao9029

金虫 (小有名气)

引用回帖:
2楼: Originally posted by htbbzzg at 2013-11-20 18:48:30
在 ANSYS 中,死单元实际上也参加计算,只是 ANSYS 把死单元的主要材料性能设置的很小,以基本不影响周围的单元 (反之则不然),因此,分析后,死单元的自由度 (热 - 温度;结构 - 位移,等) 也是会改变的。

那么怎么能将死单元的改变降到最低呢?
3楼2013-11-20 19:16:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuhaitao9029 的主题更新
信息提示
请填处理意见