【分享】借花献佛，金属浇铸的例子和帮人改进的例子

cswbww

版子里的一个朋友给我看的一个金属浇铸的例子，如下：
/clear
/prep7
et,1,55
et,2,55
kxx,1,1.99/3600
dens,1,0.00073
c,1,0.100
kxx,2,7.38e-6
dens,2,9.35e-5
c,2,1
k,1,-0.5
k,4,6
k,12
k,13,5.5
k,21,-0.5,6
k,22,,6
k,23,5.5,6
k,24,6,6
a,12,13,23,22
a,1,12,22,21
a,4,24,23,13
/pnum,line,1
lplot
lsel,s,,,5,6
lsel,a,,,9,10
lesize,all,,,8
lsel,all
lesize,all,,,20
/pnum,area,1
aplot
asel,s,,,1
aatt,2,2,2
asel,all
amesh,all
fini
/solu
antype,tran
timint,off
/pnum,line,1
lplot
lsel,s,,,7,8
nsll,s
sf,all,conv,100,70
lsel,s,,,5,10,5
nsll,s
sf,all,conv,1,70
lsel,s,,,1
nsll,s
sf,all,hflux,0
alls
esel,s,type,,2
nropt,full
ekill,all
alls
tunif,100
time,1
deltim,0.1
outres,all,all
solve
*SET,dtop,0
*SET,tim,1
*SET,delt,0.05
timint,on
tintp,,,,1
*do,ii,1,20
ddele,all,temp
esel,s,type,,2
nsle,s
*SET,tim,tim+delt
time,tim
*SET,dtop,dtop+0.3
nsel,r,loc,y,0,dtop
esln,s,1
ealive,all
nsel,r,loc,y,dtop
esln,s
nsle,s,1
d,all,temp,600
esel,s,live
nall
eall
solve
*enddo
fini
/post1
*SET,dtop,0
*do,ii,2,21
set,ii
*SET,dtop,dtop+0.3
nsel,r,loc,y,0,dtop
esln,s,1
esel,a,type,,1
nall
plns,temp
*enddo
  
然后，根据这个例子又改进成了金属入口冷却的例子，如下：
/clear
/prep7
et,1,55
et,2,55
kxx,1,1.99/3600
dens,1,0.00073
c,1,0.100
kxx,2,7.38e-6
dens,2,9.35e-5
c,2,1
k,1,-0.5
k,4,6
k,12
k,13,5.5
k,21,-0.5,6
k,22,,6
k,23,5.5,6
k,24,6,6
a,12,13,23,22
a,1,12,22,21
a,4,24,23,13
/pnum,line,1
lplot
lsel,s,,,5,6
lsel,a,,,9,10
lesize,all,,,8
lsel,all
lesize,all,,,20
/pnum,area,1
aplot
asel,s,,,1
aatt,2,2,2
asel,all
amesh,all
fini
/solu
antype,tran
timint,off
/pnum,line,1
lplot
lsel,s,,,7,8
nsll,s
sf,all,conv,100,70
lsel,s,,,5,10,5
nsll,s
sf,all,conv,1,70
lsel,s,,,1
nsll,s
sf,all,hflux,0
alls
esel,s,type,,2
nropt,full
ekill,all
alls
tunif,100
time,1
deltim,0.1       !modify
outres,all,all
solve
*SET,dtop,6
*SET,tim,1
*SET,dt,0.05
timint,on
tintp,,,,1
*do,ii,1,20
ddele,all,temp
esel,s,type,,2
nsle,s
*SET,tim,tim+dt
time,tim
*SET,dtop,dtop-0.3
nsel,r,loc,y,dtop,6
esln,s,1
ealive,all
nsel,r,loc,y,6
esln,s
nsle,s,1
d,all,temp,600
esel,s,live
nall
eall
solve
*enddo
fini
/post1
*SET,dtop,6
*do,ii,2,21
set,ii
*SET,dtop,dtop-0.3
nsel,r,loc,y,dtop,6
esln,s,1
esel,a,type,,1
nall
plns,temp
*enddo
  
希望大家看看，给点意见。还有，但愿那个朋友不会责怪我帖出来：（大家进步嘛。如果大家都有一颗博爱的心，那么请都来帖自己成功的例子吧~~~！
  
另，我改的后者，不知道有没有问题，大家看看。

sophiaboy

斑竹仁兄
  
你改的例子好像得不到结果
*do,ii,1,20  
...........
...........
*SET,dtop,dtop-0.3  
nsel,r,loc,y,dtop,6  
esln,s,1  
ealive,all  
他只是激活单元，瞬态分析的前一时刻温度不知道
  
nsel,r,loc,y,6  
esln,s  
nsle,s,1  
d,all,temp,600  
总是给入口赋定值是对的！
......
*enddo
循环的执行并没有使金属向下流动！
  
这其实和金属浇注到一个腔里有很大区别的（上贴的第一个代码）！
不知道斑竹注意到没有！
从口流入的金属：
在时间t(m)时各节点温度为T(i,j),在时间t(m+1)时，把上时刻各节点
的温度T(i,j)付给这些节点(这些节点的y坐标都小一个单位)，而空出
来的顶层节点赋初始定值！

sophiaboy

我写了一下代码，不知道怎么回事，运行得不好
前些日子好像给你邮寄过去了，我想你很忙

cswbww

sophiaboy，说的对，我没注意，我再写写。你的email我没收到，不知道怎么回事情。

sophiaboy

我再给你发一份吧，思路如上！

cswbww

ok,不过我可能会回的比较慢。

sophiaboy

谢谢你！

yjzhang_xaut

提个意见，考虑一下：金属浇注存在相变，那就必须定义材料的焓（ENTH），否则只是固体的导热，对吧？

boxsteel

焓的定义，ANSYS很简单，可以通过比热和密度的定义来实现，也可以自己按温度比热的积分计算出焓之间输入焓随温度的变化值

sophiaboy

这个问题比较好解决吧！

sophiaboy

cswbww  
我把二维改为三维的，出现错误如下！
  *** ERROR ***                           CP=     24.281   TIME= 07:37:04
  DOF (e.g.  Displacement) limit exceeded at time 5.E-02                  
   (load step 1 substep 1 equilibrium iteration 1)                        
   Max.absolute value= 3.343472078E+11 (limit= 1000000) at TEMP of node   
  11684 .                                                                 
   May be due to an unrestrained or unstable model.   
  
我的代码如下，很简单的，你们就帮忙看一下！
  
!***********
/prep7   
et,1,70
  
kxx,1,40               !materials NO1
dens,1,7500   
c,1,420
  
mp,kxx,2,413e3          !materials NO2
mp,dens,2,8900
mp,c,2,390
  
!***********
r1=0.09
r2=0.105   
h=0.70   
msh_arc=8  
msh_line=8   ! MUST BE THE SAME WITH msh_arc
msh_h=35   
space=2
  
!***********
cyl4,,,r1,0,,60,h
cyl4,,,r1,60,,120,h   
cyl4,,,r1,120,,180,h  
cyl4,,,r1,180,,240,h  
cyl4,,,r1,240,,300,h  
cyl4,,,r1,300,,360,h  
nummrg,all   
numcmp,all  
  
cyl4,,,r1,0,r2,60,h  
cyl4,,,r1,60,r2,120,h   
cyl4,,,r1,120,r2,180,h   
cyl4,,,r1,180,r2,240,h   
cyl4,,,r1,240,r2,300,h   
cyl4,,,r1,300,r2,360,h   
nummrg,all   
numcmp,all   
  
/pnum,line,1
lplot   
/ANG, 1 ,120,XS,1  
/REP,FAST  
  
!***************
vsel,,,,1,6
vatt,1,1,1
vsel,,,,7,12
vatt,2,1,1
allsel
   
!***************  
lsel,,,,10,25,5  
lsel,a,,,13,28,5
lsel,a,,,2,5,3   
lesize,all,,,msh_line,space  
  
lsel,,,,1,6,5   
lesize,all,,,msh_line,1/space  
  
lsel,,,,7,9  
lsel,a,,,14,29,5
lesize,all,,,msh_h!,space
  
lsel,,,,3,4  
lsel,a,,,30,31   
lsel,a,,,11,26,5
lsel,a,,,12,27,5
lesize,all,,,msh_arc
  
!************************
lsel,,,,32,33
lsel,a,,,36,37
lsel,a,,,40,55,5
lsel,a,,,43,58,5
lesize,all,,,msh_line
  
lsel,,,,34,35
lsel,a,,,60,61
lsel,a,,,41,56,5
lsel,a,,,42,57,5
lesize,all,,,msh_arc
  
lsel,,,,38,39
lsel,a,,,44,59,5
lesize,all,,,msh_h!,space
  
allsel   
vmesh,all   
  
!*******************solve the question*************  
/solu
antype,trans
timint,off
  
!*********add the side face hflux
allsel
csys,1
nsel,s,loc,x,0.105
sf,all,hflux,-1.5e6
  
alls
esel,s,mat,,1
nropt,full
ekill,all
alls
tunif,50
time,1
deltim,0.05
outres,basic,last
solve

sophiaboy

上面只是第一步！

limingtiger

why not plaster it, so we can share it,host

ythit

请问：
1、焓随温度变化应怎样求。难道还需知道比热随温度变化的公式才可以吗？
2、如知道一公式，在ansys中，可否直接给出公式，让ansys自行计算？
3、采用气体对一熔体冷却，气体的流量是可以改变的，这样的问题应怎样考虑计算。

terryzhou

好例子
谢谢

imlini

" 焓的定义，ANSYS很简单，可以通过比热和密度的定义来实现，也可以自己按温度比热的积分计算出焓之间输入焓随温度的变化值"
那是不是只要定义比热和密度，ANSYS可以自动计算出焓，不用再添加呢

imlini

" 焓的定义，ANSYS很简单，可以通过比热和密度的定义来实现，也可以自己按温度比热的积分计算出焓之间输入焓随温度的变化值"
那是不是只要定义比热和密度，ANSYS可以自动计算出焓，不用再添加呢