论文部分内容阅读
传热和热力耦合是科学研究和工程实际中普遍存在的自然现象,在材料结构设计、航空装备防护、自然灾害预测等领域具有广泛的应用。材料与结构一体化热力耦合模型与算法研究是力学、数学、物理学与材料科学等多学科交叉的研究领域,也是工业生产和装备制造中需要着重解决的基础问题。因此,建立多物理场耦合模型,发展高效的数值算法,研究材料和结构在热力荷载下的耦合热弹塑性动力响应和损伤破坏行为的机理,探索热力耦合行为的内在规律,对于材料设计和优化以及结构性能安全评估具有重要意义。广义插值物质点法作为一种粒子类无网格方法,在多场耦合和损伤破坏等问题的数值模拟方面具有一定的优势,在材料和结构的动力问题中得到了广泛的应用,已成为科学与工程计算中常用的数值方法之一。然而在热力耦合问题中,由于尚未发展热传导问题的物质点法模型,一般采用物质点法进行动力分析,结合其他方法进行传热分析。为了解决耦合场算法的不一致的问题,本文建立了热力耦合分析的广义插值物质点法统一求解框架。首先提出热传导问题的广义插值物质点法模型;然后发展耦合热弹塑性动力问题的广义插值物质点法及其轴对称形式;最后针对热软化和应变局部化问题建立基于Cosserat理论的热弹塑性广义插值物质点法模型。本文工作的具体内容如下:第一、针对瞬态热传导问题,提出了显式广义插值物质点法。该方法基于热传导方程的等效积分弱形式,推导了广义插值物质点法离散方程,阐述了不同温度边界的施加方法,并给出了算法流程。数值算例验证了算法的有效性和准确性,收敛性分析表明该算法随着物质点背景网格尺寸的加密和网格内物质点数量的增加而逐渐收敛。针对物质点上温度结果在空间分布的震荡问题,提出了“再映射”的方法对物质点温度进行光滑处理。此外,发展了指定温度边界的多级网格技术,有效减小了温度震荡和计算误差,使边界处的计算结果更加真实。第二、针对稳态热传导问题,建立了隐式广义插值物质点法计算格式。该方法基于稳态热传导方程的弱形式,建立了热传导矩阵的广义插值物质点法数学列式,并采用Newton-Raphson方法对热平衡方程进行迭代求解。数值算例验证了算法的在不同边界下的有效性和准确性。计算结果还表明算法随着背景网格尺寸的减小和网格内物质点数量的增加而逐渐收敛。此外,该方法可以通过增加网格内物质点的数量取得比有限元更加精确的结果。第三、考虑温度变化对结构变形的影响,发展了弱耦合热弹塑性动力问题的广义插值物质点法。该方法考虑温度变化对变形的单向耦合作用,在每一个时间步内将温度场和位移场单独求解,设计了耦合场的交错求解策略。通过代表性的算例验证了耦合算法的准确性。进一步结合基于分岔分析的位移不连续本构模型,考虑随温度变化的材料属性,模拟了预应力铝薄膜结构局部受激光照射而发生的热弹塑性破坏过程,显示了方法在处理这类问题中的适用性。第四、考虑温度和变形之间的相互耦合作用,发展了强耦合热弹塑性动力问题的广义插值物质点法。基于质量守恒、动量守恒以及能量守恒方程,该方法考虑了应变率做功对温度的影响,在热传导方程中引入弹性耦合项和塑性耦合项,采用显式时间积分设计了耦合场分析的交错求解策略。通过热弹性冲击和热弹性振动问题验证了方法的准确性。进一步考虑热塑性耦合效应,对金属杆往复弯曲生热问题进行了参数分析。此外,给出了一种基于温度变化的弹塑性本构模型,结合位移不连续破坏本构,模拟了重力和热对流共同作用的热破坏过程。第五、针对轴对称结构,建立了轴对称强耦合热弹塑性动力问题的广义插值物质点法模型。该方法推导了轴对称框架下的插值形函数,对径向坐标小于零的情况进行了坐标截断处理,修正后的形函数消除了原始形函数在靠近对称轴位置的奇异性。通过厚壁圆筒在热力荷载下的响应验证了算法在热弹性问题中的准确性与有效性。结合Johnson-Cook本构模型,将所发展的方法用于Taylor杆撞击实验过程中变形和升温现象的模拟,得到的结果与实验数据和文献报道吻合良好。通过与三维全尺寸模型计算时间的对比,体现了轴对称算法在时间上的效率优势。第六、针对材料热软化和应变局部化问题,发展了基于Cosserat连续体理论热弹塑性问题的隐式广义插值物质点法。该方法以Cosserat理论作为正则化机制,首先建立了弹塑性分析的隐式广义插值物质点法模型。结合稳态热传导分析,发展了基于Cosserat理论的热弹塑性问题广义插值物质点分析方法,采用Newton-Raphson迭代算法构造了增量格式的耦合场交错求解策略。通过分片试验验证了方法的准确性。通过矩形板的应变局部化分析,探讨了 Cosserat特征长度对剪切带形状的影响。将耦合算法应用于二维平板热力荷载作用下的热软化和应变局部化现象的模拟,结果显示本文方法克服了传统理论广义插值物质点法在应变局部化问题中的网格依赖性。最后,附录介绍了基于上述算法开发的热弹塑性问题广义插值物质点法分析程序CTGIMP。该程序支持OpenMP并行计算,实现了广义插值物质点法统一框架下材料与结构稳态和瞬态传热分析、动力分析以及耦合热弹塑性分析等功能。