许可优化
产品
解决方案
服务支持
关于
软件库
当前位置:服务支持 >  软件文章 >  ABAQUS传热分析与实践入门

ABAQUS传热分析与实践入门

阅读数 6
点赞 0
article_banner

传热问题基本概念

宏观的热能的传递有以下几种基本形式:

  • 热传导 :物体各部分之间不发生相对位移时,依靠分子、原子及自由电子等微观粒子的热运动而产生的热能传递。热传导遵循傅立叶定率, \Phi = -\lambda A\Delta T ;其中 \Phi 为热流量(单位: W), \lambda为导热系数(单位:W/(m·K)) , A 为面积, T 为温度。
  • 热对流 :流体流动导致的热量交换。对流传热遵循牛顿冷却公式, \Phi = h A\Delta T;其中 h 为对流换热系数(单位:W/(㎡·K))。
  • 热辐射 :通过电磁波辐射产生热量交换,理论上只要物体的温度高于绝对零度,物体就会不断的把热能变为辐射能,向外发出热辐射。热辐射遵循 Stefan-Boltzmann 定律 (四次方定律), \Phi = \varepsilon\sigma AT^4;其中 \varepsilon 为物体的发射率(黑度,其值总小于1), \sigma 为 Stefan-Boltzmann 常量,其值为 5.67\times 10^{-8}W/(m^2·K^4)
  • 相变:物质改变状态时会释放或吸收热量。主要形式有:冷凝-蒸发、沸腾-升华、熔化-凝固等,需要通过多相模型来模拟。


传热问题的三大类边界条件

  • 第一类边界条件,规定了边界上的温度值,也称为 Dirichlet 条件
  • 第二类边界条件, 规定了边界上的热流密度值,也称为 Neumann 条件
  • 第三类边界条件,规定了边界上物体与周围流体间的表面传热系数及周围流体的温度,也称为 Robin 条件

ABAQUS 传热分析

ABAQUS主要是有限元法 用来进行应力分析的软件,其热传导求解能力是从求解热应力问题中发展出来的。ABAQUS可以完成稳态热分析 瞬态热分析 ,允许在接触面有热流动(节点导热),可以考虑潜热项(相变)、强制对流(给定换热系数值)、热应力耦合分析功能、热传导壳单元(fluent中应避免壳单元)、空腔辐射 等功能。但 ABAQUS 不是专业热传导分析软件,不能做流体分析、自由对流换热、逆传热分析,对热冲击问题不能做自适应网格划分等。

ABAQUS 传热分析种类(此处不含电学分析)

类型 描述
非耦合传热分析 Uncoupled heat transfer analysis 包括传导、强制对流和边界辐射的传热问题可以在Abaqus/Standard中进行分析。这些分析中,温度场是在不知道所研究物体的应力/变形状态或电场的情况下计算的。Heat Transfer问题可以是瞬态的或稳态的,也可以是线性的或非线性的。
顺序耦合热应力分析 Sequentially coupled thermal-stress analysis 如果应力/位移解依赖于温度场,但没有反向依赖性,则可以在Abaqus/Standard中进行顺序耦合热-应力分析。通过首先求解热传导问题,然后读取温度场到预定义场进行应力分析。在应力分析中,温度随时间和位置变化,但不受应力分析解的影响。Abaqus允许在传热分析模型和热应力分析模型之间存在不同的网格。
完全耦合热应力分析 Fully coupled thermal-stress analysis 采用温度-位移耦合方法同时求解应力/位移场和温度场。当热解和力学解相互影响较大时,采用耦合分析的方法。例如,在快速金属加工问题中,材料的非弹性变形会引起放热。
在Abaqus/ standard中,传热方程采用后向差分格式积分,耦合系统采用牛顿法求解。这些问题可以是瞬态的或稳态的,也可以是线性的或非线性的。
在Abaqus/ explicit中,采用显式前向差分积分规则对传热方程进行积分,并采用显式中心差分积分规则获得应力响应。explicit 分析总是瞬态的。
腔辐射效应不能包括在全耦合热应力分析中。
空腔辐射 Cavity radiation 在Abaqus/ standard中,非耦合传热问题中可以包括空腔辐射效应。空腔可以是封闭的或者开放的,腔内反射和吸收可以模拟。反射因子是自动计算的,在分析过程中可以规定物体包围空腔的运动。空腔辐射问题是非线性的,可以是瞬态的也可以是稳态的。
绝热分析 Adiabatic analysis 绝热力学分析可以用于这样的情况:机械变形引起加热,但是时间很短,热传导不明显。所有的温度增加是在材料局部发生的,并且只影响局部材料点的机械属性 。绝热分析可以用Abaqus/Standard 和 Abaqus/Explicit进行,可以是静态的或动态的,也可以是线性的或非线性的。

稳态热传导案例

稳态传热分析即热交换达到平衡状态,节点温度不随时间变化,采用 Heat transfer:Steady-state 分析步。稳态分析中时间历程无意义,不代表真实的时间大小,选择 Steady-state 振幅曲线会自动调整为(Ramp)。

稳态传热分析模型
稳态分析设置

材料设置:热导率 Conductivity = 52 mW/mm/ ℃

边界条件:顶部和底部恒定热流 Heat flux = 50 mW/mm² 两端加热模式;圆柱外侧强制对流换热,表面环境温度 = 20℃;对流换热系数 Film coefficient = 1 mW/mm²/ ℃;部件初始温度10℃

对流换热边界
恒定热流边界

接下来选择单元,ABAQUS提供了热传导单元DC3D8。

HeatTransfer单元 DC3D8

设置输出热通量 HFL 和节点温度 NT ,求解后得到圆柱中心截面的温度场如下,圆柱底部中心点最高温度56.08℃。

中心截面温度场

瞬态热传导案例

瞬态热传导问题求解温度场是某一时间状态下的,采用 Heat transfer:Transient 分析步,振幅曲线自动调整为(Instantaneous),可以求解温度变化的实际过程,这里时间历程代表真实时间。

瞬态分析

瞬态分析材料设置必须加上比热容和密度:热导率 Conductivity = 52 mW/mm/ ℃,比热容 Specific heat = 434000000 mJ/ton/℃,密度 Density = 7.8e-9 ton/mm³

边界条件相同,注意这里默认热流加载是瞬态的,想要模拟恒定热流必须自定义幅值曲线,这里设置热流在1s后达到目标值并保持恒定。

定义热流载荷幅值曲线

设置求解400s, 大约260s后接近稳态结果,节点温度变化不再显著,中心节点温度56.03℃。

稳态解与瞬态解

中截面温度场变化过程如下

动图封面


免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删
相关文章
QR Code
微信扫一扫,欢迎咨询~

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 155-2731-8020
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

手机不正确

公司不为空