热传导问题的有限元方法.ppt

热传导问题的有限元方法 焊接过程的ANSYS仿真 目录 瞬态热传导方程的数值解法 瞬态温度场中n个节点温度 的有限元方程为 求解一阶偏微分方程 瞬态热传导方程的数值解法 将求解的时间域划分成若干个时间步长在一定数目的时间区域内 假设和的函数形式来近似方程的精确解仅在相隔的离散时间点上满足微分方程来代替时间域内任何时刻t都满足微分方程进一步假设时刻的解都已经求得 下一步要计算的是时刻的温度场 瞬态热传导方程的数值解法 用加权余量法建立两点循环公式在两个时间点和之间的时间区域内 采取如下线性插值形式其中 用加权余量法建立两点循环公式 由于采用近似插值 在时间域内 方程将产生余量 对于这一时间区域 典型的加权余量格式可以表示为如下形式当求解初值问题时 如果已知一组参数 则可以利用上式近似确定另一组参数 将插值函数及其导数代入加权余量表达式 经过整理 得到 其中 用加权余量法建立两点循环公式 其中 假定P采用与未知场函数 相同的插值表达式 得到 当 和P都已知时 就可以求得下一时刻的 这就是 两点循环公式 可以记成 其中 参数 的选择 前差分公式 中心差分公式 后差分公式 为常数 伽辽金型权函数 0 1 2 1 1 2 2 3 1 3 n n 1 解的稳定性问题 解的稳定性一般利用不耦合的齐次方程来讨论 解析式为 其中 是任意常数 用两点循环公式求解 令 则 避免发散 避免振荡 解的稳定性问题 避免发散 解的稳定性问题 后差分 中心差分 无条件稳定 前差分 时 稳定 避免振荡 解的稳定性问题 热应力的计算 物体由于热膨胀只产生线应变 而剪切应变为零 这种由于热变形产生的应变可以看作是物体的初应变 对于三维问题 物体存在初应变的情况下 应力应变关系可表示成 将上式代入虚位移原理的表达式 并进行有限元离散 得到 这里 焊接过程的仿真分析 二十世纪七十年代以来 国内外很多学者都对数值模拟技术在焊接中的应用进行了研究 取得了不少成果 特别是 计算焊接力学 ComputationalWeldMechanics 的发展使焊接模拟有了更为坚实的理论基础 例如 欧洲空中客车340飞机开发中 飞机机身的铝合金蒙皮壁板的纵向加强筋采用激光束焊接 主要问题是保持低变形和减少残余应力 要考虑接头类型的变化 焊接顺序 冷却条件 装夹模式及纵向预载荷等措施 决定这些措施及组合 就需要采用焊接热力数值模拟技术 焊接过程的仿真分析 焊接热力耦合分析耦合分析是指在有限元分析的过程中考虑多种物理场的交叉作用和相互影响 耦合分析最终可归结为两种不同的方法 直接耦合和顺序耦合 直接耦合 顺序耦合 包含所有必须自由度的耦合单元类型仅仅通过一次求解就能得出耦合场分析结果 按照顺序进行两次相关场分析把第一次场分析的结果作为第二次场分析的载荷 焊接过程的仿真分析 当两种物理场相互作用不明显 或者一种物理场对另一种物理场有决定性影响 而后一种物理场对前一种物理场影响较小时 进行两种物理场的完全耦合分析会使分析的问题过于复杂化 这时就可以考虑使用顺序耦合分析 顺序耦合分析具有很高的效率和灵活性 焊接过程的塑性变形热和相变潜热与焊接热输入相比 可以忽略不计 焊接热分析的温度场决定了焊接结构分析的应力场和变形场 而焊接力学场对温度场的影响较小 因此 一般进行顺序耦合热力分析 将焊接热分析各载荷步的温度场结果作为力学分析的热载荷 进行求解 焊接过程的仿真分析 焊接过程的仿真分析 焊接过程中的非线性现象 当位移显著的改变结构的刚度时 则被视为几何非线性 影响焊接热力的材料热物理参数和力学参数均与焊接热循环过程有关 是温度的非线性函数 焊接过程中会出现一系列相变 由于材料的状态不同 其本构关系也要随之变化 这种现象称为状态非线性问题 例如在低温区域使用弹塑性材料模型 而高温区域使用弹粘塑性本构模型 几何非线性 材料非线性 状态非线性 焊接过程仿真分析中存在的问题 焊接过程中的很多复杂现象之间的关系难以用准确的数学模型统一描述 移动的热源伴随着金属的熔融从而带来结构约束的不断变化 这种变化的约束关系大大增加了分析的难度 焊接结构三维分析模型的自由度数目庞大 分析规模受计算速度 内存和硬盘的限制 材料在高温阶段的热物理参数和力学参数严重缺乏 而且高温力学参数降低到很小的值 这种材料非线性影响了求解的效率 造成收敛困难 材料在较高温度区域和较低温度区域呈现不同的本构关系 多载荷步与多子步分析使得求解误差不断积累 焊接过程仿真分析的简明求解 将三维模型简化为二维甚至一维 简化构件几何和加载 将非线性热弹性 粘塑性模型简化为线性热弹性 将瞬态过程简化为准稳态过程 使热过程和力学过程分离 忽略缺陷和裂纹的形成 忽略高温发生的熔化 凝固相 以及随后在低屈服应力的相变过程 对屈服规律进行简化 简化坡口形状和焊层结构 用给定温度范围内与温度无关的平均值取代与温度相关的材料特征值 环焊缝的ANSYS仿真实例 有限元分析的步骤 结论 环焊缝的ANSYS仿真实例 结论 环焊缝长约为478mm 分为36个体 每一个体的长度约为13mm 近似熔池的长度 环焊缝的ANSYS仿真实例 生死单元让单元 死掉 并不是删除这些单元 而是将这些单元的刚度矩阵乘以一个非常小的因子 一般默认值是1 0E 6 即让其具有非常小的刚度 在热学分析过程中 所以 杀死 的单元节点同样被约束以温度载荷 直到模拟焊接过程填充经过该节点时 使其 复活 即加载相应的高温载荷 在随后的应力分析过程中 生死单元 在达到凝固温度时才被 激活 通过焊料和母材各自的热膨胀系数设定的熔化温度和环境温度作为判断 生死 的参考温度 环焊缝的ANSYS仿真实例 热学模拟选择SOLID70单元 整体模型共含有67276个单元 81190个节点 环焊缝的ANSYS仿真实例 模型中焊缝截面和与焊缝接触的管道 四通部分网格最细 边长为3mm 焊缝36个体中 每个体沿焊接方向包含5个单元 环焊缝的ANSYS仿真实例 焊接热源模型当模拟焊接热输入过程时 应用生死单元 每次激活一个体 沿焊缝方向包含5个单元 施加1600 加载时间持续10s 四通与管道由5道焊缝相连 每道焊缝之间的冷却时间为1200s 环焊缝的ANSYS仿真实例 力学模型的网格划分热学部分的网格划分较为密集 这是由于如果其网格密度过于粗大 就会导致低温现象出现 即这个模型最低温度远远低于常温20 与实际情况差别较大 从而造成误差 然而 对于结构模拟分析 倘若网格过于密集会导致计算时间过长 在结构分析中将模型重新划分网格 环焊缝的ANSYS仿真实例 力学模型的网格划分对于重新划分的网格 若想在节点施加热学部分的温度载荷 就需要用到ANSYS中的BFINT命令 对体载荷进行插值运算 该方法在很多热力耦合研究中被采纳 BFINT Fname1 Ext1 Fname2 Ext2 KPOS Clab KSHS TOLOUT TOLHGT 环焊缝的ANSYS仿真实例 力学模拟选择SOLID45单元 整体模型共含有5125个单元 6839个节点 环焊缝的ANSYS仿真实例 模型中焊缝截面和与焊缝接触的管道 四通部分网格最细 边长为5mm 焊缝36个体中 每个体沿焊接方向包含1个单元 环焊缝的ANSYS仿真实例 动画演示 难点和工作安排 计算时间