关于表面效应单元及对流换热系数的探讨

半个苹果

小弟最近要用到表面效应单元施加对流收到一个例子和大家分享，但有几个问题不明白？ 希望有经验的高手帮忙解释下， 命令流如下：
FINISH  
/CLEAR
!!!gauss加载
/title, the simulate of laser cladding process
*SET,P,2000，       !the power of the laser
*SET,V,0.008       !scanning speed   
*SET, RADIUS,0.003       !radius of the laser beam
*SET, R,0.003
*set,k,0.6          !the coefficient of the absorb
*SET,L,0.05          !the length of the block
/prep7
!熔覆层高5mm=0.005m
!基体高10mm=0.01m,长x=50mm=0.05,宽=20mm=0.02m
BLOCK,0,0.05,0, -0.0105,-0.01,0,
Vsel,all
vplot
!
/USER, 1
/VIEW, 1, -0.198026943995E-01,  0.382677277552    ,  0.923669829831  
/ANG,  1, -0.107453019382   
/REPLO
!
ET,1,SOLID70
ET,2,Surf152  
DOF,TEMP
KEYOPT,2,3,0
KEYOPT,2,4,1   
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,7,0
KEYOPT,2,8,4      !根据固体表面温度计算对流系数
KEYOPT,2,9,0
!用温度表方式定义材料性能
Mp,kxx,1,82.9
Mp,dens,1,8900   
Mp,c,1,471               
mptemp,1,20,300,600,750,900,1200,1500,1800,2000         
mpdata,hf,1,1,6,50,120,180,200,250,378,700,850
!
!
Mshkey,0        !自由网格化
Mshape,1,3d
type,1
esize,0.002,         
Vmesh,all
eplot
!  
!
asel,s,loc,y,0          !对流边界条件要施加在5个面的周围
asel,a,loc,z,-0.01
asel,a,loc,x,0
asel,a,loc,x,0.05
type,2         !定义属性
amesh,all
esel,s,type,,2         !选择表面效应单元
sf,all,conv,-1,25         !确定对流边界条件，表格加载
!sf,all,conv,200,25         !确定对流边界条件
!!底面也施加轻微的对流
allsel,all
asel,s,loc,y, -0.0105
type,2         !定义属性
amesh,all
esel,s,type,,2
sf,all,conv,3,25
allsel,all
!
nsel,all
/psf,hflux,,2,0        !以箭头方式显示载荷
!
finish
/SOLU  
!!!!!
antype,transient,new
TRNOPT,FULL  
Nropt,,,off            !关闭自适应下降
LUMPM,0
Solcontrol,on
Autots,on
Kbc,1
Outres,all,last,
Lnsrch,on
neqit,50          !每一子步中方程的迭代次数限值  
!load step1,initial conditions 25
time,0.001
deltim,0.001,0.001,0.001            
allsel,all
timint,off
tunif,25   !均匀的节点初始温度25，只对分析的第一个子步有效  
solve
EPLOT
!
! load step2………,apply moving heat flux
allsel
nsel,all
t=0
dt=0.5
local,11,1, 0.001+0.005,0,0,0,-90,0        
CSYS,11
esel,s,type,,1
nsel,r,loc,x, 0,0.003
nsel,r,loc,y,0,180
nsel,r,loc,z,0
sf,,hflux,7000000
!sf,,hflux,%flux2%     
csys,0  
t=t+dt
time,t   !每一载荷步结束时的时间
deltim,0.1,0.05,0.1,on   
kbc,1        !设定为阶越载荷
!neqit,25       !每个子步的默认次数
autots,on        !打开自动时间步长
timint,on,therm    !时间积分效果
TINTP,,,,1     
Lnsrch,on
!outres,all,all
outres,all,last   
!outres,all,10
allsel
solve
save




能否帮忙解释下 这个段命令流
希望具体些  ！

asel,s,loc,y,0          !对流边界条件要施加在5个面的周围
asel,a,loc,z,-0.01
asel,a,loc,x,0
asel,a,loc,x,0.05
type,2         !定义属性
amesh,all
esel,s,type,,2         !选择表面效应单元
sf,all,conv,-1,25         !确定对流边界条件，表格加载
!sf,all,conv,200,25         !确定对流边界条件
!!底面也施加轻微的对流
allsel,all
asel,s,loc,y, -0.0105
type,2         !定义属性
amesh,all
esel,s,type,,2
sf,all,conv,3,25
allsel,all

谢谢！！

Forthuth

asel,s,loc,y,0          !选择y=0这个面
asel,a,loc,z,-0.01     ！再选择z=-0.01这个面
asel,a,loc,x,0           ! 再选择x=0这个面
asel,a,loc,x,0.05       ! 再选择x=0.05这个面
type,2         !定义属性
amesh,all
esel,s,type,,2         !选择表面效应单元
sf,all,conv,-1,25         !确定对流边界条件，表格加载 -1是对流系数(不是太明白),25是外界温度
!sf,all,conv,200,25         !确定对流边界条件  200是对流系数，25是外界温度
!!底面也施加轻微的对流
allsel,all
asel,s,loc,y, -0.0105    !重新选择y=-0.0105这个面
type,2         !定义属性
amesh,all
esel,s,type,,2        ! 选择单元type2
sf,all,conv,3,25      ! 加载对流边界条件到刚刚选的这些面
allsel,all                ! 全部选择

大概就这些了，错了不管！

半个苹果

sf,all,conv,-1,25         !确定对流边界条件，表格加载 -1是对流系数(不是太明白),25是外界温度
!sf,all,conv,200,25         !确定对流边界条件  200是对流系数，25是外界温度
对楼上这条命令sf,all,conv,-1,25   翻译 我觉得 sf,all,conv,-1,25
-1是中的1 指的是材料编号！  但是后面 确定对流边界条件 为什么还要加个“！”

Forthuth

不知道，加个“！”应该就不管用了，你不妨直接删掉这条命令，看看结果有没有什么不一样，或者把“！”去掉再看看结果

aiyouguai

对流系数是正数才是正确的吧 楼主改一下

ansis

sf,all,conv,-1,25         !确定对流边界条件，表格加载 -1是对流系数(不是太明白),25是外界温度
!sf,all,conv,200,25         !确定对流边界条件  200是对流系数，25是外界温度
对楼上这条命令sf,all,conv,-1, ...
半个苹果 发表于 2010-5-20 09:55

前面加负号说明将采用前面所定义材料所对应的表格进行加载