| 查看: 2743 | 回复: 0 | ||
[求助]
ansys热分析 高斯热源模型 热流密度施加不上去 温度场不发生变化
|
|
我用ansys做的热分析，是对一个钢板进行的热流密度施加，但是热流密度好像施加不上去，命令流能运行，但是用后处理器查看温度场没有变化，用时间历程处理器查看加热区某一节点的温度场也没有发生变化，不知道为什么，希望哪位高手能指点一下，具体命令流如下： 【我是刚注册的号，金币不多，望请见谅！】 我也附加传了一份txt文本附件 finish /clear /filname,Fname /title,FEM model of heating /unit,SI !**********************modeling****************************** /prep7 et,1,solid70 !3D热实体 et,2,solid87 !三维10节点4面体热单元 et,3,surf152 !接触单元 keyopt,3,4,1 keyopt,3,8,1 !******************材料常数*********************************** MPTEMP,1,0,20,75,100,175,200 MPTEMP,7,225,300,325,375,400,450 MPTEMP,13,475,500,550,575,600,675 MPTEMP,19,700,720,725,775,800,900 MPTEMP,25,1000,1100,1200,1500 MPDATA,KXX,1,1,51.9,51.74,51.3,51.1,49.525,49 MPDATA,KXX,1,7,48.275,46.1,45.25,43.55,42.7,41.05 MPDATA,KXX,1,13,40.225,39.4,37.5,36.55,35.6,32.75 MPDATA,KXX,1,19,31.8,30.64,30.35,27.45,26,26.6 MPDATA,KXX,1,25,27.2,27.7,28.2,29.7 MPDATA,C,1,1,450,459.6,486,494.25,519,525.5 MPDATA,C,1,7,532,563.5,574,599,614.75,646.25 MPDATA,C,1,13,662,683.75,727.25,749,773.25,846 MPDATA,C,1,19,1139,1373.4,1432,950,931.03,855.17 MPDATA,C,1,25,779.31,703.45,627.59,400 MPDATA,DENS,1,1,7880,7880,7880,7880,7880,7800 MPDATA,DENS,1,7,7800,7800,7790,7790,7780,7750 MPDATA,DENS,1,13,7720,7710,7710,7710,7700,7700 MPDATA,DENS,1,19,7650,7650,7600,7600,7588,7500 MPDATA,DENS,1,25,7400,7320,7200,7258 MPDATA,EX,1,1,2.08e11,2.06e11,2.02e11,2.0e11,1.99e11,1.98e11 MPDATA,EX,1,7,1.92e11,1.9e11,1.87e11,1.86e11,1.85e11,1.8e11 MPDATA,EX,1,13,1.77e11,1.75e11,1.7e11,1.63e11,1.65e11,1.3e11 MPDATA,EX,1,19,1.2e11,1.38e11,1.35e11,0.9e11,0.8e11,0.6e11 MPDATA,EX,1,25,0.4e11,0.25e11,0.25e11,0.25e11 MPDATA,PRXY,1,1,0.29,0.29,0.3,0.31,0.325,0.33 MPDATA,PRXY,1,7,0.335,0.35,0.355,0.36,0.367,0.37 MPDATA,PRXY,1,13,0.38,0.39,0.395,0.4,0.405,0.42 MPDATA,PRXY,1,19,0.425,0.43,0.43,0.435,0.44,0.445 MPDATA,PRXY,1,25,0.45,0.475,0.475,0.475 MPDATA,ALPX,1,1,10e-6,10e-6,10e-6,11e-6,11e-6,11e-6 MPDATA,ALPX,1,7,11e-6,12e-6,12e-6,12e-6,12e-6,13e-6 MPDATA,ALPX,1,13e-6,13e-6,13e-6,14e-6,14e-6,14e-6,14e-6 MPDATA,ALPX,1,19,14e-6,14e-6,14e-6,14e-6,14e-6,14e-6 MPDATA,ALPX,1,25,14e-6,14e-6,15e-6,15e-6 block,0,500,0,400,0,-14 !plate block,0,40,0,250,0,-14 !flame region 1 block,0,90,0,280,0,-14 !transition region 5 block,0,90,0,400,0,-14 !up left region 2 block,90,500,0,280,0,-14 !upper right rigion 3 boptn,keep,no !布尔设置 vovl,all !体分开 numcmp,kp !压缩代号 numcmp,line numcmp,area numcmp,volu !*********************meshing************************ !meshing the flame rigion type,1 lsel,s,line,,21,22,1 lsel,a,line,,5,8,3 lesize,all,8,,,1 lsel,s,line,,6,7,1 lsel,a,line,,23,24,2 lesize,all,8,,,1 lsel,s,line,,1,11,10 lesize,all,,,3,2 lsel,a,line,,9,10,1 lesize,all,,,3,0.5 lsel,all vsel,s,volu,,1 vmesh,all vsel,all !meshing the linear rigion type,1 lsel,s,line,,25,26,1 lsel,a,line,,39,40,1 lesize,all,30,,,1 lsel,a,line,,29,30,1 lesize,all,30,,,1 lsel,s,line,,13,14,1 lsel,a,line,,33,34,1 lesize,all,30,,,1 lsel,s,line,,12,15,3 lsel,a,line,,27,28,1 lesize,all,40,,,1 lsel,s,line,,31,32,1 lsel,a,line,,37,38,1 lesize,all,40,,,1 lsel,s,line,,35,36,1 lesize,all,40,,,1 lsel,all vsel,s,volu,,2,4,1 vmesh,all vsel,all !meshing the transition region !过渡网格划分 type,2 lsel,s,line,,43,44,1 lesize,all,,,2,1 lsel,s,line,,41,42,1 lesize,all,,,3,1 mopt,expnd,1.5 mopt,tetexpnd,2 mopt,trans,1.5 lsel,all vsel,s,volu,,5 vmesh,all vsel,all !表面有效单元的划分 type,3 asel,s,area,,8 amesh,all asel,all numcmp,node numcmp,elem save finish !****************施加热流 heating flux ********************* finish /clear resume,Fname,db /filname,Fname1 /title,termperature filed of heating *set,vatl,817.3 !flow rate of acetylene(l/h) *set,ata,0.29 !heat transfer efficiency *set,v,4 !flame moving speed *set,tmptint,23 !initial temperature *set,k,3/40/40 !for heating source *set,q0,k*ata*11470*4.184*vatl/3600/3.1416 *set,timestep,8/v !time step length(time stayed on each element) *set,tend,210/v !the end time of heating !************************** define components ****************** !define the node component of flame region. asel,s,area,,8 !upper flame area. nsla,s,1 !all nodes. cm,flmrg,node !flame region component !define the upper surface asel,s,area,,12,17,5 asel,a,area,,22,26,4 !upper surface except flame area. nsla,s,1 !all nodes cm,upsfc,node !define the whole upper surface cmsel,s,flmrg cmsel,a,upsfc cm,whupsfc,node !define the below surface asel,s,area,,7,11,14 asel,a,area,,16,21,5 asel,a,area,,25 !below surface. nsla,s,1 cm,blsfc,node !define the around surface asel,s,area,,5,10,5 asel,a,area,,14,15,1 asel,a,area,,19,20,1 asel,a,area,,23 !around surface. nsla,s,1 cm,rdsfc,node !define the 152 surface element component esel,s,type,,3 cm,htflrg,elem !heat flux elements !**************************** heating ****************************** /solu antype,TRANS !solution setting tunif,tmptint !uniform temperature for the fist load step TREF,TMPTINT nsubst,2,5,1 !two substeps. kbc,0 !ramped load. outres,nsol,last !write node sesults at the last substep outpr,all,none !none print out. eqslv,sparse,le-5 !pcg solver 精度设置 !**************************** **************************** !load step loop during heating *set,Tstart,0.01 *set,Tend,tend *set,Tinc,timestep *do,TM,Tstart,Tend,Tinc *set,xorg,0 *set,yorg,v*(TM-Tstart) *set,zorg,0 clocal,11,1,xorg,yorg,zorg *if,TM-Tstart,le,Tend,then !flame heating csys,11 !flame round *set,R0,40 !the outside flame round radius *set,Rend,8 !the inner flame round rsdius *set,Rinc,-8 !flame round step. *do,r,R0,Rend,Rinc cmsel,s,flmrg !nodes in flame region surface. nsel,r,loc,x,0,r-0.05*8 esln esel,r,type,,3 *set,Rmid,r+Rinc/2 sfe,all,,hflux,,q0*exp(-1*k*Rmid*Rmid) *enddo cmsel,s,flmrg nsel,r,loc,x,0,R0 esln,s,1 cm,flame,elem !current flame round cmsel,s,htflrg !surface 152 elelements cmsel,u,flame !non current flame region sfedel,all,,hflux !delete last flame load *endif !delete heat flux after the last flame heating step. *if,TM-Tstrat,eq,Tend+timestep,then cmsel,s,flame sfedel,all,,hflux !delete heat flux *endif allsel time,TM csys,0 solve *enddo |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 热流施加 1.txt
2012-05-11 16:06:23, 7.07 K
» 猜你喜欢
本科211，293分请求调剂
已经有4人回复
-
求调剂
已经有3人回复
-
085400电子信息319求调剂（接受跨专业调剂）
已经有5人回复
-
268求调剂
已经有9人回复
-
一志愿双非085502，267分，过四级求调剂
已经有3人回复
-
一志愿085404，总分291，四级已过，求调剂
已经有4人回复
-
297求调剂
已经有12人回复
-
26调剂 086003
已经有3人回复
-
265求调剂
已经有20人回复
-
271分求调剂学校
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
找到一些相关的精华帖子，希望有用哦~
高斯计算的问题？？
已经有5人回复
- 高斯 gaussian 在 windows7 中充分利用 多核
已经有19人回复
- 急求用高斯09怎么计算电子激发态
已经有7人回复
- 关于高斯稳定结构的问题
已经有3人回复
- 高斯计算两个有机分子
已经有11人回复
- 高斯计算中分子轨道问题
已经有6人回复
- 请问牛顿法和高斯牛顿法有什么不同？非常感谢
已经有6人回复
- fluent中的体热源加载问题（实际上是热流密度，不是热生成率）
已经有19人回复
- 在进行高斯计算时如何消去虚频（转载）
已经有18人回复
- 高斯计算中zORA方法是否可以使用
已经有9人回复
- 高斯中键伸缩力常数与键角力常数的计算
已经有6人回复
- dmol 中的GGA-RPBE相当于高斯中的什么计算方法？
已经有5人回复
- 【求助】高斯计算过渡态不收敛，请问如何解决？
已经有8人回复
- 【讨论】ansys如何加载高斯分布的热流密度？
已经有13人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
