热传导问题的有限元方法

,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,热传导问题的有限元方法,焊接过程的,ANSYS,仿真,目录,瞬态热传导方程的数值解法,1,焊接过程的仿真分析,2,环焊缝的,ANSYS,仿真实例,3,难点和工作安排,4,瞬态热传导方程的数值解法,瞬态温度场中,n,个节点温度,的有限元方程为,时间积分,模态叠加,求解一阶偏微分方程,瞬态热传导方程的数值解法,将求解的时间域 划分成若干个时间步长,在一定数目的 时间区域内，假设 和 的函数形式来近似方程的精确解,仅在相隔 的离散时间点上满足微分方程来代替时间域内任何时刻,t,都满足微分方程,进一步假设 时刻的解都已经求得，下一步要计算的是 时刻的温度场,瞬态热传导方程的数值解法,用加权余量法建立两点循环公式,在两个时间点 和 之间的 时间区域内，采取如下线性插值形式,其中，,用加权余量法建立两点循环公式,由于采用近似插值，在时间域 内，方程将产生余量，对于这一时间区域，典型的加权余量格式可以表示为如下形式,当求解初值问题时，如果已知一组参数 ，则可以利用上式近似确定另一组参数 。将插值函数及其导数代入加权余量表达式，经过整理，得到,其中,用加权余量法建立两点循环公式,其中,假定,P,采用与未知场函数,相同的插值表达式，得到,当,和,P,都已知时，就可以求得下一时刻的,这就是,两点循环公式，可以记成,其中,参数,的选择,前差分公式,中心差分公式,后差分公式,为常数,伽辽金型权函数,=0,=1/2,=1,=1/2,=2/3,=1/3,n,n+1,解的稳定性问题,解的稳定性一般利用不耦合的齐次方程来讨论,解析式为,其中,是任意常数，,用两点循环公式求解,令,则,避免发散,避免振荡,解的稳定性问题,避免发散,解的稳定性问题,（后差分、中心差分）无条件稳定,（前差分）,时，稳定,避免振荡,解的稳定性问题,热应力的计算,物体由于热膨胀只产生线应变，而剪切应变为零。这种由于热变形产生的应变可以看作是物体的初应变 ，对于三维问题，,物体存在初应变的情况下，应力应变关系可表示成,将上式代入虚位移原理的表达式，并进行有限元离散，得到,这里，,焊接过程的仿真分析,二十世纪七十年代以来，国内外很多学者都对数值模拟技术在焊接中的应用进行了研究，取得了不少成果。特别是“计算焊接力学（,Computational Weld Mechanics,）”的发展使焊接模拟有了更为坚实的理论基础。,例如，欧洲空中客车,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,Ext 1,-,Fname2,Ext 2,-,KPOS,Clab,KSHS,TOLOUT,TOLHGT,环焊缝的,ANSYS,仿真实例,力学,模拟选择SOLID,4,5,单元，整体模型共含有,5,125,个单元，,6,839,个节点,。,环焊缝的,ANSYS,仿真实例,模型中焊缝截面和,与焊缝接触的管道、四通部分网格最细，边长为,5mm,。焊缝,36,个体中，每个体沿焊接方向包含,1,个单元。,环焊缝的,ANSYS,仿真实例,动画演示,难点和工作安排,计算时间长、需要硬盘空间大,需要详细的焊接方案,热源模型的建立,材料属性,4,3,2,1,难点,难点和工作安排,从局部到整体,，先对转向架的一根侧梁进行仿真分析，计算温度场和位移场，再计算整个转向架。预计在,10,月份之前，完成侧梁的焊接模拟。,7.158.31,几何清理、建立热源模型、计算温度场,9.19.30,计算位移场、应力场,