求助！杀死单元后进行温度场分析为什么结果显示单元上还有温度分布？

ytyzwwj0420

求助！杀死单元后进行温度场分析为什么结果显示单元上还有温度分布？
图4中的那个体2是一个竖向焊缝，进行温度场分析前 将体2所附的单元都杀死了 然后在体1上施加移动热源进行温度场分析 但是为什么结果显示体2所附的单元上还是有温度分布？   命令流在2楼

ytyzwwj0420

finish
/clear
/filename,cy
/title,cy
/CONFIG,NRES,50000
/prep7
/units,si
wpstyl,,,,,,,,1
csys,4
!!!!!!!定义几何参数
*set,d0,0.09            !!!主管直径
*set,t0,0.008           !!!主管厚度
*set,L0,12*d0           !!!主管长度
*set,d1,0.025           !!!支管直径
*set,t1,0.008           !!!支管厚度
*set,L1,0.3             !!!支管长度
*set,lbc,0.05           !!!肋板长
*set,lbg,0.05           !!!肋板高
*set,lbh,0.008          !!!肋板厚
*set,hf,0.008           !!!焊缝
!!!!!!!!!!!!!!!!!!!!!!热源
pi=3.1415926
I=160
U=23
efficient=0.78
lsize=0.002
v=0.004
tmg=lbg/v
tmc=lbc/v
tinc=lsize/v
t=0
!!!!!!!!!!!!!!!!!!!!!!!定义材料参数
et,1,solid70
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,7850
MPTEMP,,,,,,,,  
MPTEMP,1,18
MPTEMP,2,300
MPTEMP,3,600
MPTEMP,4,900   
MPTEMP,5,1200   
MPDATA,KXX,1,,53.3
MPDATA,KXX,1,,44
MPDATA,KXX,1,,34
MPDATA,KXX,1,,27.3
MPDATA,KXX,1,,27.3
MPTEMP,,,,,,,,  
MPTEMP,1,18
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,18
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
MPTEMP,,,,,,,,  
MPTEMP,1,0
MPDATA,ALPX,1,,1.75E-005
MPTEMP,,,,,,,,  
MPTEMP,1,18
MPTEMP,2,300
MPTEMP,3,700
MPTEMP,4,1100
MPTEMP,5,1500
MPTEMP,6,2000
MPDATA,EX,1,,2.06E+11   
MPDATA,EX,1,,1.7E+11   
MPDATA,EX,1,,2.68E+10      
MPDATA,EX,1,,5.0E+9
MPDATA,EX,1,,3.0E+9
MPDATA,EX,1,,1.4E+9
MPDATA,PRXY,1,,0.25
MPDATA,PRXY,1,,0.25
MPDATA,PRXY,1,,0.25
MPDATA,PRXY,1,,0.25
MPDATA,PRXY,1,,0.25
MPDATA,PRXY,1,,0.25
tb,bkin,1,5              
tbtemp,18,1
tbdata,1,330e6,
tbtemp,500,2
tbdata,1,213e6,
tbtemp,800,3
tbdata,1,153e6,
tbtemp,1200,4
tbdata,1,73e6,
tbtemp,1600,5
tbdata,1,13e6,
mp,ex,2,2.06e16           !!!加载板
mp,prxy,2,0.3
mp,dense,2,7850
MP,KXX,2,21.4  
MP,C,2,502
MP,ALPX,2,1.75E-005
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建模zhuzhiguan
cyli,d0-t0,d0,-L0/2,L0/2
wprota,,90
cyli,d1-t1,d1,0,L1+d0
wprota,,,90
wpoffs,,,lbh/2
blc4,,,-(d0+hf+lbg),(d1+hf+lbc),-lbh
vsel,all
vovl,all,
vsel,s,,,11,12
vsel,a,,,4,16,6
vdele,all,,,1
!!!!!!!!!!!!!!生成加载板
wpcsys,-1,0
wprota,,90
wpoffs,,,L1+d0-t0
cyli,0,d1,0,t0
vsel,s,,,1
vsel,a,,,14
vovl,all
alls
!!!!!!!!!!!!1/2结构
wpcsys,-1,0
vsbw,all
vsel,s,loc,z,0,-L0/2
vdel,all,,,1
allse
vplot
!!!切割
wprota,,90
vsbw,all
lwplan,-1,68
vsbw,all
wpoffs,,,-d0
vsbw,all
wpoffs,,,lbc+d0
vsbw,all
lwplan,-1,65
vsbw,all
wpoffs,,,-d1
vsbw,all
wpoffs,,,lbg+d1
vsbw,all
wpcsys,-1,0
wprota,,,90
vsel,s,loc,y,0,-d0
vsel,r,loc,x,-(d1+hf+lbc),-L0/2
vsel,a,loc,y,-(d0+hf+lbg),-(d0+L1)
vsel,a,loc,y,0,d0
vsbw,all
alls
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!注意角度的变化
wprota,,30
vsbw,all
wprota,,-60
vsbw,all
alls
vsel,s,,,37
vdele,all,,,1
vsel,all
vglue,all
vsel,s,,,36
cm,shuhanfeng,volu
vsel,s,,,35
cm,henghanfeng,volu
vsel,s,,,34
cm,leiban,volu
vsel,s,,,35,36
cm,hanfeng,volu
wpcsys,-1,0
vsel,s,loc,y,-(L1+d0-t0),-(L1+d0)
cm,jzb,volu
alls
!cmsel,s,jzb
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!mesh
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!mesh on weld
wpcsys,-1,0
cmsel,s,shuhanfeng,volu
type,1
mat,1
esize,0.003
mshkey,1
vsweep,all
cmsel,s,henghanfeng,volu
type,1
mat,1
esize,0.003
mshkey,1
vsweep,all
cmsel,s,leiban,volu
type,1
mat,1
esize,0.004
mshkey,1
vsweep,all
alls
vplot
!!!!!!!!!!!!!!!!!!!!!!!!!zhuguan
vsel,s,,,65
type,1
mat,1
esize,0.006
mshkey,1
vsweep,all,104,271
vsel,s,,,57
type,1
mat,1
esize,0.006
mshkey,1
vsweep,all,103,237
vsel,s,,,14,16,2
vsel,a,,,23
vsel,a,,,6
type,1
mat,1
esize,0.006
mshkey,1
vsweep,all
vsel,s,,,11
type,1
mat,1
esize,0.006
mshkey,1
vsweep,all,49,42
vsel,s,,,12
type,1
mat,1
esize,0.006
mshkey,1
vsweep,all,43,52
vsel,s,,,8
type,1
mat,1
esize,0.006
mshkey,1
vsweep,all
vsel,s,,,21,25,4
vsel,a,,,58,60,2
vsel,a,,,66
type,1
mat,1
esize,0.006
mshkey,1
vsweep,all
vsel,s,,,2,18,8
vsel,a,,,9
vsel,a,,,27
vsel,a,,,41
type,1
mat,1
esize,0.008
mshkey,1
vsweep,all
vsel,s,loc,y,0,d0
vsel,r,loc,z,0,d1+hf+lbc+d0
type,1
mat,1
esize,0.01
mshkey,1
vsweep,all
vsel,s,loc,z,d1+hf+lbc+d0,L0/2
type,1
mat,1
mshkey,1
esize,0.015
vsweep,all
!!!!!!!!!!!!!!!!!!!!!!!!!zhiguan
vsel,s,loc,y,-d0,-(d0+hf+lbg)
vsel,r,loc,z,0,d1
type,1
mat,1
mshkey,1
esize,0.006
vsweep,all
vsel,s,loc,y,-(d0+hf+lbg),-(d0+hf+lbg+d1)
type,1
mat,1
mshkey,1
esize,0.009
vsweep,all
vsel,s,loc,y,-(d0+hf+lbg+d1),-(d0+L1-t0)
type,1
mat,1
mshkey,1
esize,0.015
vsweep,all
!!!!!!!!!!!!!!!!加载板
vsel,s,loc,y,-(d0+L1-t0),-(d0+L1)
vsel,u,,,3,20,17
type,1
mat,2
mshkey,1
esize,0.006
vsweep,all
vsel,s,,,3
type,1
mat,2
esize,0.006
mshkey,1
vsweep,all,24,25
vsel,s,,,20
type,1
mat,2
esize,0.006
mshkey,1
vsweep,all,25,30
alls
!!!!!!!!!!!!!!!!!!!!!!!!进入求解
OUTPR,nsol
OUTRES,nsol
/solu
antype,4                           ! 指定分析类型为瞬态分析
trnopt,full                        ! 指定瞬态分析的求解方法为完全法
nropt,full
timint,on
tintp,0.005,,,1,0.5,0.2            ! 设置瞬态积分参数
autots,on                          ! 打开自动时间步
kbc,1                              ! 阶跃载荷变化
solcontrol,on
lnsrch,on  
tunif,30
vsel,all
asel,s,ext
sfa,all,,conv,110,30                !!!施加对流边界条件
alls
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!与主管焊接
allsel
lwplan,-1,205,0.5       !!!!!!!!!!!!!!!!!!!加热中心(heng向)
!!!!!!!!!!!!!!!!!!!!!热源形状系数
a=0.008           !!!对于X形状系数
b=0.0045           !!!对于y形状系数
c=0.0045           !!!对于z形状系数
heat=6*sqrt(3.0)*I*U*efficient/(a*b*c*pi*sqrt(pi))
tmg=2.5
tmc=2.5
allsel
!  *do,t,1e-4,15,tinc*0.25
  *do,t,1e-4,15,tinc
  *if,t,le,tmc+1e-4,then              !!!注意 有一个1e-5的差
    cmsel,s,shuhanfeng
    eslv,s
    ekill,all
    time,t
    hengfeng.mac
    alls
    deltim,0.005,0.005,0.2            
    solve
    nsel,all
    bfdele,all,hgen
  *else  
    cmsel,s,shuhanfeng
    eslv,s
    ealive,all
    allsel
    lwplan,-1,148,0.5       !!!!!!!!!!!!!!!!!!!加热中心
!      *do,i,1e-4,tmg+10,tinc*0.25
       *do,i,1e-4,tmg+2,tinc
        *if,i,le,tmg+1e-4,then  
          time,i+tmc+1e-4
          shufeng.mac
          alls
          deltim,0.005,0.005,0.2         
          solve
          nsel,all
          bfdele,all,hgen
        *else
          time,i+tmc+1e-4
          alls
          deltim,0.005,0.005,0.2         
          solve
        *endif
    *enddo
  *exit
  *endif
*enddo

fly_tangwei

哥们，你问的问题解决了么？ 知道是什么原因么 可以分享解释一下么，有同样的迷惑？
  发现版主大哥 好长时间都不上论坛解决问题啊

lwpnmamtf

生死单元技术中杀死单元相当于给单元乘以一个很小的因子，使其表现为零，但之前的温度分布没有置零。
所以在进行后处理如etable命令时，要记得将死单元排除后再进行。

fly_tangwei

lwpnmamtf 发表于 2012-6-25 22:37
生死单元技术中杀死单元相当于给单元乘以一个很小的因子，使其表现为零，但之前的温度分布没有置零。
所以 ...

怎么排除死单元呢？用什么命令？
但是我发现有些例子的死单元被杀后温度还会上升，这是什么原因呢？

lwpnmamtf

排除死单元可以用esel命令，将生存单元选择后在反选。
死单元温度上升，可能是后续杀死的单元温度较高，它与之前杀死的单元间温度进行了热量交换导致。
但死亡单元不再属于我们关注的范围，不需要去理会