请问在ansys中如何在同一表面上同时施加热流密度和对流载荷

smartboy

如题！谢谢!

kitty827

使用表面效应单元施加对流

smartboy

我是初学者,  能给个例子吗? 谢谢!!

crazy_leo

谁有例子共享一下啊,我也不知道怎么加,郁闷啊

不倒臂膀

若要同时施加对流荷载和热流密度，需要用到表面效应单元，surf151和surf152，前者用于二维热分析，后者为三维热分析。
一般情况下，将热对流边界施加到表面效应单元，而将热流密度施加到实体单元上
对表面效应单元：设有孤立节点，则对流系数施加到表面效应单元，流体温度施加到孤立节点上，如果没有孤立节点，则对流系数和流体温度都施加在表面效应单元
生成表面效应单元的一般顺序：
lsel,s,,,
type,1！表面效应单元
mat,1
esurf,N!生成带孤立节点Ｎ的表面单元

crazy_leo

按照楼上的说法,我现在碰到的问题是这样的，假如命令流是这样的:
/PREP7
ET,1,PLANE55
et,2,surf151
………………
生成了一个矩形
………………
TYPE,2
REAL,2
MAT,2
*GET,NN,NODE,,num,max
LSEL,S,,,2
NSLL,S,1
LSEL,a,,,3
NSLL,a,1
N,NN+1,0,0.2,0
ESURF,NN+1

！以上，我定义了一个表面效应单元
！以下，我想在这个表面效应单元上加对流载荷，我分

/SOLU
time 1
load hflux in 2 side and 3 side
solve

time 2
load conv in 2 side and 3 side
solve

!我在time 2的时候不会使用表面效应单元，所以总提示，conv override hfulx
我是这样的写的
ESEL,s,TYPE,,2
LSEL,s,,,2
NSLL,s,1
LSEL,a,,,3
NSLL,a,1
SF,all,conv,800,27
……………………

不知道哪错了

不倒臂膀

首先楼上的建模是否啰嗦了一下，
LSEL,S,,,2
NSLL,S,1
LSEL,a,,,3
NSLL,a,1
上面这几句完全可以用
LSEL,S,,,2,3
NSLL,S,1来代替，没必要一个个的选！
错误提示意思好像是conv把hflux给覆盖了，原因我估计是没用到表面单元，因为你是使用了时间步，在此完全没有必要，你可以
antype,0
nsel,s,,,
sf,all,hflux,!施加热流密度
esel,s,mat,,2
sfe,all,1,conv,1,!施加对流面荷载
楼上可看那样如何！

crazy_leo

楼上的,谢谢你在建模中提醒,我没用多长时间,所以不太熟悉,我的命令流是这样的:
/FILNAME,CONTINUOUS CASTING
/TITLE,CONTINUOUS CASTING


/PREP7
ET,1,PLANE55
et,2,surf151

KEYOPT,2,4,1
KEYOPT,2,5,1
KEYOPT,2,6,0
KEYOPT,2,7,1
KEYOPT,2,8,2
KEYOPT,2,9,1

r,2,1,5.67e-8

mptemp,1,211,800,911,1392
mpdata,c,1,,482.166587,622.093856,648.463851,762.733828
mpdata,kxx,1,,44.605510,31.300000,28.792510,17.926720
mpdata,enth,1,,0,7.5e9,9.6e9,1.1e10
mpdata,dens,1,,7779.765800,7637.110000,7610.225800,7493.727600
MP,EMIS,2,0.8

RECTNG,0,0.75,0,0.125
LSEL,S,,,1
LESIZE,ALL,,,50,0.2
LSEL,S,,,3
LESIZE,ALL,,,50,5
LSEL,S,,,2
LESIZE,ALL,,,20,0.2
LSEL,S,,,4
LESIZE,ALL,,,20,5
AMESH,1

TYPE,2
REAL,2
MAT,2
*GET,NN,NODE,,num,max
LSEL,S,,,2,3
NSLL,S,1
N,NN+1,0,0.2,0
ESURF,NN+1
ALLSEL                           
SAVE,124v925,db              
SAVE            


/SOLU

ANTYPE,TRANS
TIMINT,ON
bfunif,temp,1550

!结晶器-----------------------------------------------------------
*set,result1,0
*set,result11,0

*set,result2,0
*set,result21,0

*set,a,-6614.7
*set,b,10234
*set,c,-5811.9
*set,d,2300.9

*set,e,-6814
*set,f,10592
*set,g,-60503
*set,h,2023.1


*DO,i,1,42
!宽面
*set,a1,a*i*i*i/3600/60
*set,a2,b*i*i/3600+c*i/60
*set,result1,a1+a2+d
*set,result11,(result1+result11)
LSEL,S,,,3
NSLL,S,1
SF,ALL,HFLUX,result11

!窄面
*set,m,i
*set,a3,e*m*m*m/3600/60
*set,a4,f*m*m/3600+g*m/60
*set,result2,a3+a4+h
*set,result21,(result2+result21)
LSEL,a,,,2
NSLL,a,1
SF,ALL,HFLUX,result21

TIME,i
autos,1
KBC,1
NSUBST,2
ALLSEL,ALL
LSWRITE,i
*ENDDO

OUTRES,ALL,ALL
LSSOLVE,1,42,1



!二冷区---------------------------------------------------------------------------------
!zhc198111@163.com
antype,4,rest
ESEL,s,TYPE,,2
LSEL,s,,,2
NSLL,s,1
LSEL,a,,,3
NSLL,a,1
SF,all,conv,800,27
TIMINT,ON
TIME,300
deltim,15,15,30
autots,on
KBC,1
NSUBST,2
OUTRES,ALL,ALL
SOLVE

!辐射---------------------------------------------------------------------------------------
TIME,1026
deltim,15,15,30
TOFFST,273
AUTOTS,ON
ALLSEL
d,nn+1,temp,27
KBC,1
NSUBST,2
OUTRES,all,ALL
outpr,nsol,last
SOLVE
finish

/POST1
/PLOPTS,INFO,ON
/TITLE,CONTINUOUS CASTING
PLNSOL,TEMP
FINISH

你说的
"sf,all,hflux,!施加热流密度
esel,s,mat,,2
sfe,all,1,conv,1,!施加对流面荷载
是直接加完热流密度,就加对流了,但我的必须是,先求解"结晶器"(边界是热流密度),然后在求解"二冷"(边界是对流),不知道怎么写?楼上能否在指点一下.小弟先谢了!

crazy_leo

并且,我试了一个照你说的例子,好像也不行,同样也是提示:conv覆盖hflux

/FILNAME,124v925
/TITLE,CONTINUOUS CASTING


/PREP7
ET,1,PLANE55
et,2,surf151

KEYOPT,2,4,1
KEYOPT,2,5,1
KEYOPT,2,6,0
KEYOPT,2,7,1
KEYOPT,2,8,2
KEYOPT,2,9,1

r,2,1,5.67e-8

mptemp,1,211,800,911,1392
mpdata,c,1,,482.166587,622.093856,648.463851,762.733828
mpdata,kxx,1,,44.605510,31.300000,28.792510,17.926720
mpdata,enth,1,,0,7.5e9,9.6e9,1.1e10
mpdata,dens,1,,7779.765800,7637.110000,7610.225800,7493.727600
mpdata,hf,2,,800,800,800,800
MP,EMIS,2,0.8

RECTNG,0,0.75,0,0.125
LSEL,S,,,1
LESIZE,ALL,,,50,0.2
LSEL,S,,,3
LESIZE,ALL,,,50,5
LSEL,S,,,2
LESIZE,ALL,,,20,0.2
LSEL,S,,,4
LESIZE,ALL,,,20,5
AMESH,1

TYPE,2
REAL,2
MAT,2
*GET,NN,NODE,,num,max
LSEL,S,,,2
NSLL,S,1
LSEL,a,,,3
NSLL,a,1
N,NN+1,0,0.2,0
ESURF,NN+1
ALLSEL                           
SAVE,124v925,db              
SAVE            


/SOLU
!第一个载荷
ANTYPE,TRANS
bfunif,temp,1550
LSEL,S,,,2,3,1
NSLL,S,1
sf,all,hflux,-3000   !施加热流密度
esel,s,mat,,2
sfe,all,1,conv,800  !施加对流
TIME,550
KBC,1
ALLSEL,ALL
OUTRES,ALL,ALL
solve

不倒臂膀

楼上说的没错，原因在于我们一直在mat=2的单元上加对流和热流密度，表面效应单元没用上，现改正如下：
fini
/clear
/PREP7
ET,1,PLANE55
et,2,surf151

KEYOPT,2,4,1
KEYOPT,2,5,1
KEYOPT,2,6,0
KEYOPT,2,7,1
KEYOPT,2,8,2
KEYOPT,2,9,1

r,2,1,5.67e-8
mptemp,1,211,800,911,1392
mpdata,c,1,,482.166587,622.093856,648.463851,762.733828
mpdata,kxx,1,,44.605510,31.300000,28.792510,17.926720
mpdata,enth,1,,0,7.5e9,9.6e9,1.1e10
mpdata,dens,1,,7779.765800,7637.110000,7610.225800,7493.727600
mpdata,hf,2,,800,800,800,800
MP,EMIS,2,0.8

RECTNG,0,0.75,0,0.125
LSEL,S,,,1
LESIZE,ALL,,,50,0.2
LSEL,S,,,3
LESIZE,ALL,,,50,5
LSEL,S,,,2
LESIZE,ALL,,,20,0.2
LSEL,S,,,4
LESIZE,ALL,,,20,5
AMESH,1


*GET,NN,NODE,,num,max
LSEL,S,,,2,3
NSLL,s,1
TYPE,2
MAT,2
N,NN+1,0,0.2,0
ESURF,NN+1
ALLSEL   

esel,s,mat,,1
nsle,s,1
nsel,r,loc,y,0.125
nsel,a,loc,x,0.75
sf,all,hflux,-3000 !施加热流密度

esel,s,mat,,2
sfe,all,1,conv,1,800  !施加对流
另外楼上看好sfe后的选项，应该是sfe,all,1,conv,1,800，不是sfe,all,1,conv,800

crazy_leo

恩,楼上说的很对,我跟据你的,完整了一下(基本没变你的),但结果很出人意料
命令流是这样的:
/PREP7
ET,1,PLANE55
et,2,surf151

KEYOPT,2,4,1
KEYOPT,2,5,1
KEYOPT,2,6,0
KEYOPT,2,7,1
KEYOPT,2,8,2
KEYOPT,2,9,1

r,2,1,5.67e-8
mptemp,1,211,800,911,1392
mpdata,c,1,,482.166587,622.093856,648.463851,762.733828
mpdata,kxx,1,,44.605510,31.300000,28.792510,17.926720
mpdata,enth,1,,0,7.5e9,9.6e9,1.1e10
mpdata,dens,1,,7779.765800,7637.110000,7610.225800,7493.727600
mpdata,hf,2,,800,800,800,800
MP,EMIS,2,0.8

RECTNG,0,0.75,0,0.125
LSEL,S,,,1
LESIZE,ALL,,,50,0.2
LSEL,S,,,3
LESIZE,ALL,,,50,5
LSEL,S,,,2
LESIZE,ALL,,,20,0.2
LSEL,S,,,4
LESIZE,ALL,,,20,5
AMESH,1


*GET,NN,NODE,,num,max
LSEL,S,,,2,3
NSLL,s,1
TYPE,2
MAT,2
N,NN+1,0,0.2,0
ESURF,NN+1
ALLSEL   


/SOLU
TIMINT,ON
time,500
bfunif,temp,1550
esel,s,mat,,1
nsle,s,1
nsel,r,loc,y,0.125
nsel,a,loc,x,0.75
sf,all,hflux,-30000

esel,s,mat,,2
sfe,all,1,conv,1,800  
solve


/POST1
/PLOPTS,INFO,ON
/TITLE,CONTINUOUS CASTING
PLNSOL,TEMP
FINISH
我发不上来附件,结果只显示一条线,希望楼上的能帮解决一下,谢谢您!

crazy_leo

并且,好像有效单元没有作用

不倒臂膀

楼上难道没发现求解之前，出现一个warning吗？这是因为你没全选，改正如下：
esel,s,mat,,2
sfe,all,1,conv,1,800  
allsel！全选
solve
结果如图
当然不是你所说的一条线，你看再看看！

新元素

xiexiefenxiang!!! :)

crazy_leo

不倒臂膀大哥,别笑话我啊,我是接学不久,真是不太懂,老师上去就让我学命令流,我就看别人的例子,然后自己写,所以才会这样的,谢谢啊!

我还有一个问题,假如我这个模型是动态的,也就是说,对流和热流密度不是同时加,在第一个载荷步(time 1)加热流密度,然后在第二个载荷步(time 2)加对流,这样也同样可以类似的这样做吗?

不倒臂膀

我没试过，从技术上讲应该是可以的！你可以在模型上试试
我接触Ansys也不是很久，你老师让你直接学习命令流是对的，若用gui方式，对复杂模型不实用，最好的是两者互为使用。当然，我们都需要再努力！

crazy_leo

我试了一下,提示覆盖(如果我没写错),但跟着ansys又继续算下去了,最好还出来了一个结果,但我不知道,这个结果是不是覆盖后的结果,还是这两个载荷步,叠加后的结果.
/FILNAME,Test
/TITLE,TEST

/PREP7
ET,1,PLANE55                        ! define element type 1, PLANE55
KEYOPT,1,1,0                        ! set keyoption- film coefficient
KEYOPT,1,3,1                        ! set keyoption-for axisymmetric
KEYOPT,1,4,0                        ! parallel to global
KEYOPT,1,8,0                        ! exclude mass transport
KEYOPT,1,9,0                        ! std heat transfer
ET,2,SURF151                        ! define element type 2, SURF151
KEYOPT,2,3,1                        ! set keyoption-axisymmetric
KEYOPT,2,4,1                        ! exclude midside nodes
KEYOPT,2,5,1                        ! include extra node
KEYOPT,2,6,0                        ! use of bulk temp
KEYOPT,2,8,2                        ! HF at average T
KEYOPT,2,9,1                        ! exclude radiation

r,2,1,5.67e-8

mptemp,1,211,800,911,1392
mpdata,c,1,,482.166587,622.093856,648.463851,762.733828
mpdata,kxx,1,,44.605510,31.300000,28.792510,17.926720
mpdata,enth,1,,0,7.5e9,9.6e9,1.1e10
mpdata,dens,1,,7779.765800,7637.110000,7610.225800,7493.727600

MP,EMIS,2,0.8

RECTNG,0,0.75,0,0.125
LSEL,S,,,1
LESIZE,ALL,,,50,0.2
LSEL,S,,,3
LESIZE,ALL,,,50,5
LSEL,S,,,2
LESIZE,ALL,,,20,0.2
LSEL,S,,,4
LESIZE,ALL,,,20,5
AMESH,1


*GET,NN,NODE,,num,max
TYPE,2
MAT,2
N,NN+1,0,0.2,0
NSEL,S,LOC,X,0.75
NSEL,A,LOC,Y,0.125                       
NSEL,A,NODE,,NN+1
ESURF,NN+1
ALLSET
FINISH       

/SOLU
!第一阶段------------------------------------------------------------------------------
ANTYPE,TRANS
TOFFST,273,
time,100
bfunif,temp,1550
nsel,r,loc,y,0.125
nsel,a,loc,x,0.75
sf,all,hflux,-3000
solve

!第二阶段---------------------------------------------------------------------------------
antype,4,rest
TIMINT,ON
time,300
ESEL,S,TYPE,,2                        ! select all type 2 elements  
SFE,ALL,1,CONV,,800                                        ! apply film coeff to elements
ALLSEL,ALL                        ! select everything  
D,NN+1,TEMP,27                              ! apply temp to extra node                       
SOLVE                                           ! solve

crazy_leo

调好使了,不倒臂膀大哥

mlf

不倒臂膀大哥, 可以留个联系方式吗？
我也有类似的问题想请教您。我的QQ号码是：231456023

yidexian

谢谢,不倒臂膀大哥