[转帖]ANSYS热分析例子(Solidification)

VibInfo · 发表于 2005-9-26 09:41

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

/BATCH ! ! /GRA,POWER /GST,ON /filnam,flywheel ! specify jobname rad1 = 0.010 ! define inside radius (1 cm) rad2 = 0.050 ! define hub outside radius (5 cm) rad3 = 0.175 ! define inside radius of rim (17.5 cm) rad4 = 0.200 ! define outside radius of rim (20 cm) radf = 0.010 ! define radius of fillet at rim ( 1 cm) hgt1 = 0.015 ! define 1/2 thickness of hub (1.5 cm) hgt2 = 0.050 ! define 1/2 thickness of rim (5 cm) hgt3 = 0.003 ! define 1/2 thickness of disk at rim (0.3 cm) radm = 0.250 ! define radius of the sand mold (25 cm) hgtm = 0.100 ! define 1/2 height of the sand mold (10 cm) hside = 7.50 ! define outside wall film coefficient htop = 5.75 ! define outside top film coefficient /prep7 ! enter the preprocessor /title, Chapter 9 - Phase Change - Sand Casting of a Flywheel et,1,plane55,,,1 ! define element type 1, PLANE55, axisymmetric mp,kxx,1,0.346 ! define thermal conductivity of sand mp,dens,1,1520 ! define specific heat of sand mp,c,1,816 ! define specifi heat of sand mptemp ! clear the temperature table mptemp,1,100,200,300,400,530,800 ! define temperature table mpdata,kxx,2,1,206,215,228,249,268,290 ! define thermal conductivity of Al mptemp ! clear the temperature table mptemp,1,0,695,697,1000 ! define a temperature table aldens = 2707 ! density of AL alcs = 896 ! specific heat of solid AL allat = 395441 ! latent heat of AL alcl = 1050 ! specific heat of liquid AL cavg = (alcs+alcl)/2 ! average specific heat cstar = cavg+allat/(697-695) ! cstar hs = aldens*alcs*(695-0) ! enthaply of solid AL at melting temp hl = hs+aldens*cstar*(697-695) ! enthalpy of liquid AL at melting temp h1 = 0 ! enthalpy at 0 degrees C h2 = hs ! enthalpy at 695 degrees C h3 = hl ! enthalpy at 697 degrees C h4 = hl+aldens*alcl*(1000-697) ! enthalpy at 1000 degrees C ! mpdata,enth,2,1,h1,h2,h3,h4 ! define material property, enthalpy mpplot,kxx,2 ! plot KXX versus temperature mpplot,enth,2 ! plot enthalpy versus temperature ! k,1,rad1 ! create flywheel geometry k,2,rad4 k,3,rad4,hgt2 k,4,rad3,hgt2 k,5,rad3,hgt3 k,6,rad2,hgt1 k,7,rad1,hgt1 l,2,3 *repeat,5,1,1 l,1,7 k,8,rad2 k,9,rad3-radf l,1,8 l,8,9 l,9,2

VibInfo · 发表于 2005-9-26 09:41

LFILLT,3,4,radf /pnum,kp,1 lplot /pnum,kp,0 a,1,8,6,7 ! define areas a,8,9,11,6 a,9,2,10,11 a,10,2,3,4 ARSYM,Y,ALL ! reflect geometry across x axis nummrg,kp ! merge coincident keypoints RECTNG,0,radm,0,hgtm ! create sand mold areas rectng,0,radm,0,-hgtm RECTNG,0,rad1,0,hgt1 RECTNG,0,rad1,0,-hgt1 RECTNG,rad4,radm,0,hgt2 RECTNG,rad4,radm,0,-hgt2 RECTNG,0,radm,hgt2,hgtm RECTNG,0,radm,-hgtm,-hgt2 AOVLAP,ALL nummrg,kp ! merge coincident keypoints numcmp,all ! compress numbers /pnum,line,1 lplot /pnum,line,0 /pnum,area,1 aplot /pnum,area,0 eshape,2 ! specify quad element shapes esize,rad1/2 ! set element size mat,2 ! set attributes for AL real,2 type,1 amesh,1,8 ! mesh flywheel areas mat,1 ! set attributes for sand real,1 type,1 amesh,10,11 ! mesh the sand areas in the hole of the flywheel esize,rad1 ! set new element size amesh,9,12,3 ! mesh areas outside rim lsel,s,,,2 ! select and concatenate lines for quad meshing lsel,a,,,27 lsel,a,,,36 lccat,all lsel,s,,,23 lsel,a,,,35 lsel,a,,,39 lccat,all allsel ! amesh,13,14 ! mesh remaining areas eshape,0 ! set element shape for mixed quad/tri amesh,15,16 ! mesh remaining areas finish ! exit the preprocessor !

VibInfo · 发表于 2005-9-26 09:42

/solu ! Enter Solution antype,trans ! specify transient analysis timint,off ! turn off time integration time,1e-3 ! use a small value of time for this load step deltim,1e-3 ! run one substep to initialize temperature esel,s,mat,,1 ! select sand elements and nodes nsle,s,1 d,all,temp,25 ! set initial temperature of sand esel,s,mat,,2 ! select all aluminum elements and nodes nsle,s,1 d,all,temp,750 ! set initial temperature of aluminum lsel,s,loc,x,radm ! select lines on OD for convection loading sfl,all,conv,hside,,30 ! specify outside wall convection lsel,s,loc,y,hgtm ! select lines for top convection loading sfl,all,conv,htop,,30 ! specify top convection load 
 allsel ! select everything ! solve ! solve ddele,all,temp ! delete initial temperature specs timint,on ! turn on time integration tintp,,,,1 ! specify use EULER bakward diff. integration lnsrch,on ! turn on linesearch autots,on ! auto time stepping on neqit,99 ! limit number of equilibrium iterations to 99 time,2400 ! set final time to 40 minutes deltim,0.01,0.0001,100 ! use small initial time step outres,all,1 ! store all solutions on the results file solve ! solve the transient load step /psf,conv,,2 ! turn on plotting of convection eplot ! plot elements /psf,default ! reset plotting finish save ! save the database ! /post26 ! enter the time-history postprocessor solu,2,dtime,,deltime ! DELTIM for each step ! nsol,3,node((rad1+rad2)/2,0,0,),temp,,T1 ! define checking points T1-T4 nsol,4,node(rad2,0,0),temp,,T2 nsol,5,node(rad3,0,0),temp,,T3 nsol,6,node((rad3+rad4)/2,0,0),temp,,T4 nsol,7,node(0,hgtm,0),temp,,SandTOP ! define checking points for mold nsol,8,node(radm,hgtm,0),temp,,SandCRNR nsol,9,node(0,-hgtm,0),temp,,SandBOT ! plvar,2 ! plot results variables defined above plvar,3,4,5,6 plvar,7,8,9 *GET,nsets,VARI,,NSET ! get the number of datasets on the results file ! finish ! exit the time-history postprocessor ! /post1 ! enter the postprocessor set,first ! read in first set of results ! esel,s,mat,,2 ! select AL elements and nodes nsle,s ! /cval,1,695,900 ! set two uneven contours, at 695 and 900 ! *do,i,2,nsets ! use a do loop to view results set,next plnsol,temp *enddo ! finish ! exit the postprocessor ! /exit,nosa

xianglove · 发表于 2005-9-26 20:12

你知道"八趣通宝"吗? 八趣致力于与全球网民和各网站主一起创造互惠互利的双赢局面,给大家提供一个创造财富的机会. 赶快
申请加入吧. 用八趣通宝的理由: 一.无需任何点击,跑时间积分.就像网易POPO一样,对工作无任何影响.不像NEW**ar(新闻时报)或SABar等要及时保
存积分. 二.真正的无干扰.隐藏时,只在任务栏中有一个小"铜钱",继续赚钱,使用过程中始终没有弹出任何页面. 三.绿色小软件无需安装. 四.无需一分钱的投入,免费注册使用,因而无任何风险. 五.下线数与积分兑现金.通宝支持五级下线:10%, 5%, 3%, 3%,3% 自己挂机上网,(一天限制4小时),一个月也就是7元,但就是7元,也比NEW**ar(新闻时报)或SABar等容易赚
多了.如有下线收入相当可观了 . 公司的收入当然靠广告商了.他们有信心,用户会阅读并点击广告来获得额外的积分以增加收入.从而也达到
了广告的效果. 六. .如果不是"网虫",一定要想办法发展十来个下线. .自己上网时,就挂上 ,无意中也就能积点分,久而久之,也
能挣些上网费,不正白不挣,是吧? 公司现在是发展初期,上升空间非常巨大,你发展下线空间同样也很大, 赶快行动吧?!心动不如行动,早加入,早一
天的回报,多一天的收入. 如你还未注册,那就不要犹豫了,赶紧注册吧!! 注册后,可详阅读该网站说明. 注册网址:http://www.8qu.net/register.asp?net=ylwl2005

水天凝香 · 发表于 2006-5-13 12:41

好

漓江烟雨 · 发表于 2006-5-14 10:41

NICE!!!

hello2001 · 发表于 2006-5-30 21:48

研究一下！

sc108 · 发表于 2006-6-27 09:17

好东西，顶一下。 辛苦了。

214216709 · 发表于 2006-6-28 10:45

好

账号		自动登录	找回密码
密码			我要加入

[转帖]ANSYS热分析例子(Solidification)

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（VibInfo）[转帖]ANSYS热分析例子(Solidificat...

回复：（VibInfo）[转帖]ANSYS热分析例子(Solidificat...

kl