请问在ansys中如何在同一表面上同时施加热流密度和对流载荷
如题!谢谢! 使用表面效应单元施加对流回复 #2 kitty827 的帖子
我是初学者,能给个例子吗? 谢谢!! 谁有例子共享一下啊,我也不知道怎么加,郁闷啊 若要同时施加对流荷载和热流密度,需要用到表面效应单元,surf151和surf152,前者用于二维热分析,后者为三维热分析。一般情况下,将热对流边界施加到表面效应单元,而将热流密度施加到实体单元上
对表面效应单元:设有孤立节点,则对流系数施加到表面效应单元,流体温度施加到孤立节点上,如果没有孤立节点,则对流系数和流体温度都施加在表面效应单元
生成表面效应单元的一般顺序:
lsel,s,,,
type,1!表面效应单元
mat,1
esurf,N!生成带孤立节点N的表面单元 按照楼上的说法,我现在碰到的问题是这样的,假如命令流是这样的:
/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,!施加对流面荷载
楼上可看那样如何!
回复 #7 不倒臂膀 的帖子
楼上的,谢谢你在建模中提醒,我没用多长时间,所以不太熟悉,我的命令流是这样的:/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,!施加对流面荷载
是直接加完热流密度,就加对流了,但我的必须是,先求解"结晶器"(边界是热流密度),然后在求解"二冷"(边界是对流),不知道怎么写?楼上能否在指点一下.小弟先谢了!
回复 #7 不倒臂膀 的帖子
并且,我试了一个照你说的例子,好像也不行,同样也是提示: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
回复 #10 不倒臂膀 的帖子
恩,楼上说的很对,我跟据你的,完整了一下(基本没变你的),但结果很出人意料命令流是这样的:
/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
我发不上来附件,结果只显示一条线,希望楼上的能帮解决一下,谢谢您! 并且,好像有效单元没有作用 楼上难道没发现求解之前,出现一个warning吗?这是因为你没全选,改正如下:
esel,s,mat,,2
sfe,all,1,conv,1,800
allsel!全选
solve
结果如图
当然不是你所说的一条线,你看再看看! xiexiefenxiang!!!:) :)
回复 #13 不倒臂膀 的帖子
不倒臂膀大哥,别笑话我啊,我是接学不久,真是不太懂,老师上去就让我学命令流,我就看别人的例子,然后自己写,所以才会这样的,谢谢啊!我还有一个问题,假如我这个模型是动态的,也就是说,对流和热流密度不是同时加,在第一个载荷步(time 1)加热流密度,然后在第二个载荷步(time 2)加对流,这样也同样可以类似的这样做吗? 我没试过,从技术上讲应该是可以的!你可以在模型上试试
我接触Ansys也不是很久,你老师让你直接学习命令流是对的,若用gui方式,对复杂模型不实用,最好的是两者互为使用。当然,我们都需要再努力!
回复 #16 不倒臂膀 的帖子
我试了一下,提示覆盖(如果我没写错),但跟着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 调好使了,不倒臂膀大哥
回复 #18 crazy_leo 的帖子
不倒臂膀大哥, 可以留个联系方式吗?我也有类似的问题想请教您。我的QQ号码是:231456023 谢谢,不倒臂膀大哥
页:
[1]
2