ansys多物理耦合场有限元分析

《ansys多物理耦合场有限元分析》由会员分享，可在线阅读，更多相关《ansys多物理耦合场有限元分析（162页珍藏版）》请在金锄头文库上搜索。

1、ANSYS多物理耦合场有限元分析,ANSYS多物理耦合场有限元分析,结构-热耦合分析 流体-固体耦合分析,ANSYS中的典型物理量( 国际单位制 ),温度 热流量 热传导率 密度 比热 对流换热系数 热流 温度梯度 内部热生成,Degrees C ( or K ) Watts Watts/ ( meter.degree C ) kilogram/ ( meter3 ) ( Watt.sec ) / ( kilogram .degree C) Watt/ ( meter2.degree C ) Watt/ ( meter2 ) degree C / meter Watt/ ( meter3 ),

2、ANSYS热分析,热传递的类型,热传递有三种基本类型: 传导 - 两个良好接触的物体之间或一个物体内部不同部分之间由于温度梯度引起的能量交换。 对流 - 在物体和周围流体之间发生的热交换。 辐射 - 一个物体或两个物体之间通过电磁波进行的能量交换。 在绝大多数情况下，分析的热传导问题都带有对流和/或辐射边界条件。,ANSYS热分析,传导引起的热通量流由传导的傅立叶定律决定： 负号表示热量沿梯度的反向流动 (例如, 热量从热的部分流向冷的部分).,传导,ANSYS热分析,对流,对流引起的热通量由冷却牛顿定律得出: 对流一般作为面边界条件施加,ANSYS热分析,热力学第一定律,能量守恒要求系统的能

3、量改变与系统边界处传递的热和功数值相等。 能量守恒在一个微小的时间增量下可以表示为方程形式 将其应用到一个微元体上，就可以得到热传导的控制微分方程。,ANSYS热分析,单元类型 下表显示通常使用的热单元类型。 节点自由度是：TEMP。,常用的热单元类型,材料特性 至少需要 Kxx 稳态分析热传导系数。 如果是瞬态分析，则需要比热 (C) 。 优先设置为 “thermal” （热分析），在 GUI 方式中只显示热材料特性。 实常数 主要用于壳和线单元。,热分析有限元模型,ANSYS热分析,稳态热传递,如果热量流动不随时间变化的话，热传递就称为是稳态的。 由于热量流动不随时间变化, 系统的温度和热

4、载荷也都不随时间变化。 由热力学第一定律，稳态热平衡可以表示为: 输入能量 输出能量 = 0,ANSYS热分析,稳态热传递控制方程,对于稳态热传递，表示热平衡的微分方程为:,相应的节点处的有限元平衡方程为:,ANSYS热分析,热载荷和边界条件的类型,温度 自由度约束，将确定的温度施加到模型的特定区域。 均匀温度 可以施加到没有温度约束的所有节点上。可以在稳态或瞬态分析的第一个子步对所有节点施加初始温度而非约束。它也可以在非线性分析中用于估计随温度变化材料特性的初值。 热流率 是集中节点载荷。正的热流率表示热量流入模型。热流率同样可以施加在关键点上。此载荷通常用于不能施加对流和热通量的情况下。施

5、加该载荷到热传导率相差很大的区域上时应注意。,ANSYS热分析,热载荷和边界条件的类型,对流 施加在模型外表面上的面载荷，模拟模型表面与周围流体之间的热量交换。 热通量(热流密度) 同样是面载荷。当通过面的热流率已知的情况下使用。正的热流密度值表示热量流入模型。 热生成率 作为体载荷施加，代表体内生成的热，单位是单位体积内的热流率。,ANSYS热分析,热载荷和边界条件的类型,ANSYS 热载荷分为四大类: 1. DOF 约束 - 指定的 DOF (温度) 数值 2. 集中载荷 - 施加在点上的集中载荷(热流) 3. 面载荷 - 在面上的分布载荷(对流、热流密度) 4. 体载荷 - 体积或场载荷

6、（热生成）,ANSYS热分析,热载荷和边界条件注意事项 在 ANSYS中, 没有施加载荷的边界作为完全绝热处理。 通过施加绝热边界条件（缺省条件）得到对称边界条件。 如果模型某一区域的温度已知，就可以固定为该数值。 反作用热流率只在固定了温度自由度时才具有。,热载荷和边界条件的类型,ANSYS热分析,何为瞬态分析?,由于受随时间变化的载荷和边界条件,如果需要知道系统随时间的响应，就需要进行瞬态分析 。 热能存储效应在稳态分析中忽略，在此要考虑进去。时间，在稳态分析中只用于计数，现在有了确定的物理含义。 涉及到相变的分析总是瞬态分析。,时变载荷,时变响应,ANSYS热分析,除了导热系数 (k),

7、 还要定义密度 (r) 和 比热 (c ) 。 稳态分析和瞬态分析对明显的区别在于加载和求解 过程。,* MASS71热质量单元比较特殊，它能够存贮热能单不能传导热能。因此，本单元不需要热传导系数。,瞬态分析前处理考虑因素,ANSYS热分析,控制方程,回忆线性系统热分析的控制方程矩阵形式。热存储项的计入将静态系统转变为瞬态系统: 在瞬态分析中，载荷随时间变化. . . . . . 或,对于非线性瞬态分析, 时间 和 温度:,热存储项 = (比热矩阵) x (时间对温度的微分),ANSYS热分析,选择合理的时间步很重要，它影响求解的精度和收敛性。,如果时间步长 太小, 对于有中间节点的单元会形成

8、不切实际的振荡，造成温度结果不真实。,时间步大小建议,如果时间步长 太大, 就不能得到足够的温度梯度。,一种方法是先指定一个相对较保守的初始时间步长，然后使用自动时间步长按需要增加时间步。下面说明使用自动时间步长大致估计初始时间步长的方法。,ANSYS热分析,在瞬态热分析中大致估计初始时间步长，可以使用Biot和Fourier数。 Biot 数 是无量纲的对流和传导热阻的比率: 其中 D x是名义单元宽度, h是平均对流换热系数，K 是平均导热系数。Fourier 数 是无量纲的时间(Dt/t ) , 对于宽度为D x 的单元它量化了热传导与热存储的相对比率: 其中 r 和 c 是平均的密度和

9、比热。,时间步大小说明 (续),ANSYS热分析,如果Bi 1: 时间步长可以用Fourier 和 Biot数的乘积预测: 求解 D t 得到: (Again, where 0.1 b 0.5) 时间步长的预测精度随单元宽度的取值，材料特性的平均方法和比例因子b 而变化。,时间步大小说明 (续),ANSYS热分析,进行瞬态分析,ANSYS缺省情况下是稳态分析。使用下列求解菜单指定要进行瞬态分析:,“FULL” 是瞬态热分析唯一可以使用的选项。,7. 用户要输入求解选项，并不是只对热分析有效 (如求解器，N-R 选项等),ANSYS热分析,初始条件,初始条件 必须对模型的每个温度自由度定义，使得

10、时间积分过程得以开始。 施加在有温度约束的节点上的初始条件被忽略。 根据初始温度域的性质，初始条件可以用以下方法之一指定:,注: 如果没有指定初始温度，初始DOF数值为0。,ANSYS热分析,均匀初始温度,如果整个模型的初始温度为均匀且非0，使用下列菜单指定:,1,2,ANSYS热分析,非均匀的初始温度,如果模型的初始温度分布已知但不均匀，使用这些菜单将初始条件施加在特定节点上:,4. 用图形选取或输入点号的方法确定要建立初始温度的节点。 5. 单击 OK. 注: 当手动或借助于输入文件输入IC命令时，可以使用节点组元名来区分节点。,1,2,3,ANSYS热分析,非均匀初始温度 (续),注:

11、没有定义DOF初始温度的节点其初始温度缺省为TUNIF命令指定的均匀数值。,6.选择 DOF 标记 “TEMP”。 7. 指定初始温度数值。 8. 完成后单击OK。单击APPLY重复操作，将初始温度指定到其它节点上。,ANSYS热分析,由稳态分析得到的初始温度 (续),当模型中的初始温度分布是不均匀且未知的，单载荷步的稳态热分析可以用来确定瞬态分析前的初始温度。要这样做，按照下列步骤: 1. 稳态第一载荷步: 进入求解器，使用稳态分析类型。 施加稳态初始载荷和边界条件。 为了方便，指定一个很小的结束时间 (如1E-3 秒)。避免使用非常小的时间数值 ( 1E-10) 因为可能形成数值错误。 指

12、定其它所需的控制或设置 (如非线性控制)。 求解当前载荷步。,ANSYS热分析,施加瞬态分析控制和设置。 求解之前, 打开时间积分: 求解当前瞬态载荷步。 求解后续载荷步。时间积分效果保持打开直到在后面的载荷步中关闭为止。,由稳态分析得到的初始温度 (续),2. 后续载荷步为瞬态: 在第二个载荷步中，根据第一个载荷步施加载荷和边界条件。记住删除第一个载荷步中多余的载荷。,1,2,3,4,ANSYS热分析,打开/关闭时间积分效果,象刚刚说明的那样, 稳态分析可以迅速的变为瞬态分析，只要简单的在后续载荷步中将时间积分效果打开。 同样，瞬态分析可以变成稳态分析，只要简单的在后续载荷步中将时间积分效果

13、关闭。 结论: 从求解方法来说，瞬态分析和稳态分析的差别就在于时间积分。,ANTYPE,TRANS + TIMINT,OFF ANTYPE,STATIC ANTYPE,STATIC + TIMINT,ON ANTYPE,TRANS,ANSYS热分析,另外的时间积分例子,在本例中，不是在分析的开始关闭时间积分效果来建立初始条件，而是在分析的结束关闭时间积分来“加速”瞬态。,通常，分析的目标将将瞬态热现象中最严重的温度梯度定量。这些梯度通常在瞬态的初始阶段发生,并在系统进入稳态时随时间衰减。 当系统响应稳定后，后面的结果就没有意义了，分析可以简单的结束或如果稳态温度场也需要得到，就在最后载荷步关闭

14、时间积分效果。,注意改变到稳态边界时的突变。最后一个载荷步的终止时间可以是任意的,但必须比前面的瞬态载荷步时间数值要大。,ANSYS热分析,打开控制,打开控制 用于在当瞬态热分析接近稳态时让自动时间步 “打开”(增加)时间步长。在缺省情况下，如果连续3个子步间的最大温度变化都小于 0.1个温度单位，那么时间步长将迅速增加以提高效率。这个控制只能在求解控制中实现。用这些菜单改变设置:,3. 指定温度。 4. 指定门槛值。 5. 指定子步数。 6. 单击OK。,1,2,ANSYS热分析,阶跃还是渐变?,要准确模拟系统的瞬态响应，载荷必须以正确的幅值，在正确的时间和正确的速率施加。 回忆一下载荷在载

15、荷步中相对时间可以是阶跃的或渐变的: ANSYS 缺省是渐变加载的。渐变加载可以提高瞬态求解的适应性，如果有非线性时可以提高收敛性。参考第4章学习ANSYS如何处理渐变载荷。,ANSYS热分析,阶跃还是渐变? (续),要模拟阶跃载荷，将载荷在很短的时间内渐变施加到全值，然后在后续载荷步中保持不变。,问题: 对茶壶进行瞬态热分析。在底上施加热流模拟炉子的加热。热流载荷应该是阶跃的还是渐变的如果 . . . 1. 茶壶在一个刚燃着的炉子上 2. 茶壶载一个已经很热的炉子上,ANSYS热分析,什么是耦合场分析?,耦合场 分析考虑两个或两个以上的物理场之间的相互作用。这种分析包括直接和间接耦合分析。,

16、当进行直接耦合时, 多个物理场（如热电）的自由度同时进行计算。这称为直接方法，适用于多个物理场各自的响应互相依赖的情况。由于平衡状态要满足多个准则才能取得，直接耦合分析往往是非线性的。每个结点上的自由度越多，矩阵方程就越庞大，耗费的机时也越多。,下表列出了ANSYS中可以用作直接耦合分析的单元类型。不是所有单元都有温度自由度。,结构-热耦合分析,什么是耦合场分析? (续),间接耦合分析是以特定的顺序求解单个物理场的模型。前一个分析的结果作为后续分析的边界条件施加。有时也称之为序贯耦合分析。 本分析方法主要用于物理场之间单向的耦合关系。例如，一个场的响应（如热）将显著影响到另一个物理场（如结构）的响应，反之不成立。本方法一般来说比直接耦合方法效率高，而且不需要特殊的单元类型。 本章中我们只讨论涉及热的耦合现象。请注意并非所有ANSYS产品都支持所有耦合单元类型和分析选项。例如，ANSYS/Thermal产品只提供热电直接耦合。详