宏观的热能的传递有以下几种基本形式:
传热问题的三大类边界条件
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。
设置输出热通量 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℃。
中截面温度场变化过程如下