ANSYS辐射仿真模拟
2016-10-17 by:CAE仿真在线 来源:互联网
辐射传热过程是是借助于电磁波的能量传播过程,是由物体内部微观粒子在运动状态改变时所激发出来的。由于辐射传热引起的热流与物体表面热力学温度的4次方成正比,因此辐射传热分析是高度非线性的。借助于温度场数值模拟仿真技术,可以了解研究热辐射规律,对于炉内传热的合理设计十分重要,对于高温炉操作工的劳动保护也有积极意义。
本文基于大型有限元软件ANSYS对辐射传热过程温度场模拟仿真,随着ANSYS版本不断更新,核心技术不断完善,其稳态瞬态热分析、辐射热分析、相变分析、热应力分析和流体热分析功能不断强大,更能显示其计算精度与计算速度的良好兼顾性。
1 辐射传热过程温度场模拟仿真1.1研究对象本文研究的同轴圆柱体尺寸如图所示:
图1 研究模型
1.2基本假设在复杂的辐射传热过程实际条件下,抓住主要方面模拟实验情况,做一些合理化的假设,但同时又能保证其结果的准确性。本文做如下假设:
1)由于两个圆柱体足够长,将问题简化为平面问题;
2)考虑到整个辐射传热过程为封闭系统,不需设置空间节点。
1.3初始条件假设圆柱体是瞬时传热的。圆柱体为已知初始均匀温度场,即:
T(x,y,z,t=0)=T
T为圆柱体温度,即100°C.
1.4 边界条件传热是在圆柱体内径行的的,所以把外圆柱体当做边界条件。
外圆柱体的初始温度:100°C
辐射率:1
两圆柱体的辐射传热用Newton冷却定律描述:
式中:α为对流换热系数,α=65 W/m2·℃;Tf为液态金属的特征温度;Tw为砂型边界温度。
辐射传热后,两圆柱体之间的导热主要以不稳定导热方式进行。三维不稳定热传导方程为:
式中:ρ为密度,kg/m3;c为定压比热容,J/(kg·℃);t为温度,℃;T为时间,s;λ为热导率,W/(m·℃);Q为内热源密度(此处为金属液凝固时释放的潜热),W/m3。
因为整个辐射传热过程为封闭系统,所以不必考虑两圆柱体与外界的传热。
材料性能参数
密度ρ1 (kg/m3) |
密度ρ2 (kg/m3) |
导热系数K1 W/(m·℃) |
导热系数K2 W/(m·℃) |
比热容C1 J/(kg·℃) |
比热容C2 J/(kg·℃) |
7800 |
2700 |
70 |
30 |
460 |
320 |
根据圆柱体结构特征,定义其ANSYS单元类型为Thermal Solid实体单元;而后采用ANSYS中的实体建模工具构建其模型,建模完毕对铸件局部采用映射(Mapped)方式网格剖分。其他部分则利用智能网格划分器自由剖分,以达到所需部位网格细化的目的,从而兼顾计算精度和运算速度。铸件划分网格后的有限元模型如图2所示。
图2.划分网格
1.7加载求解指定分析类型为Steady-state;先作稳态分析,确定本文第1. 3节初始条件及本文第1. 4节边界条件,设置稳态分析时间步长为O.Ols,施加温度载荷、对流载荷,得到初始温度场分布,将其作为整个瞬态分析过程的初始温度场;打开时间积分选项,设置计算终止时间为16200s,进行瞬态分析,着重研究该时间段内温度场的变化规律。
1.8模拟结果及分析
图3
图4
图5
从上面分析可以看出,如果不采取措施,实验结果可能不准确。使结果不准确的主要因素是多方面的,例如它们的吸收率(或反射率)、物体的形状和大小及其相互间的位置与距离的影响。
因此可以采用隔热性能良好的材料,避免将由壁面 以对流和辐射两种方式散失于周围环境中的热量。
从本例建立有限元模型、设置材料性能参数等可知,可以选择适当壁厚、绝热性能良好的材料,来控制零件辐射传热过程温度场分布。
2 结语
本文在结合辐射传热过程的基础上,给出一种对其温度场应用ANSYS软件模拟仿真的简单方法。该方法充分利用零件的结构特征,选取最恰当的单元类型,不但计算简单,省时省力,而且误差较小,精度较高,从而达到了兼顾计算精度和计算时间的模拟要求;该模拟结果对于炉内传热的合理设计十分重要,对于高温炉操作工的劳动保护也有积极义。
依据ANSYS模拟温度场,对于炉内传热的合理设计十分重要,对于高温炉操作工的劳动保护也有积极意义。当某系统需要保温时,即使此系统的温度不高,辐射传热的影响也不能忽视。如保温瓶胆镀银,就是为了减少由辐射传热造成的热损失。
相关标签搜索:ANSYS辐射仿真模拟 Ansys有限元培训 Ansys workbench培训 ansys视频教程 ansys workbench教程 ansys APDL经典教程 ansys资料下载 ansys技术咨询 ansys基础知识 ansys代做 Fluent、CFX流体分析 HFSS电磁分析 Abaqus培训