关键词:
近场动力学
扩散方程
波动方程
δ收敛
影响函数
非局部理论
非局部性
摘要:
近场动力学(Peridynamics,PD)是连续介质力学基于非局部理论的重构,由于其控制方程中不包含对空间的微分项,因此特别适合用来处理腐蚀、裂纹等包含不连续性的问题。近些年来,近场动力学得到了越来越多的关注,并发展出了许多数值方法来求解PD模型。得到PD方程的解析解有利于更好地理解PD模型。目前,对PD解析解的研究主要集中在无限域,对有限域问题解析解的研究鲜有报道。此外,PD模型非局部性的定量衡量方法还有所欠缺,非局部行为长时间演化的特点尚不明确,PD的δ收敛特性研究尚不完善,对影响函数的研究也较为片面。因此,本文推导了有限域上的近场动力学解析解,并对非局部行为的演化进行分析,探究影响函数对近场动力学模拟结果的影响规律。本文的主要工作和结论如下:
(1)提出了近场动力学有限域问题的分离变量求解法,基于该方法推导获得了一维和二维有限域近场动力学瞬态扩散方程的级数形式解析解。定义了衡量近场动力学核函数非局部性的函数——非局部因子。与经典偏微分方程的解析解比较后发现,PD解析解可以通过在相应经典偏微分方程的分离变量级数解的时间项中添加非局部因子直接获得。非局部因子越接近1,PD模型的非局部性越弱,PD解越接近经典解。研究了近场域大小、核函数以及边界条件对非局部传热行为演化的影响,发现在扩散问题中,非局部效应随着时间的推移逐渐减弱。进一步分析了非局部边界条件和局部边界条件之间的关系,指导了数值模型中虚节点的应用。还对扩散模型的δ收敛特性进行了研究。结果表明,在扩散早期δ收敛可能不是单调的,某一时刻后δ收敛的单调性不会发生改变。
(2)基于上述方法推导获得了一维和二维有限域近场动力学波动方程的级数形式解析解。利用非局部因子得到了近场动力学的频散和群速度公式。研究了近场域大小、核函数和初始条件对非局部波动行为演化的影响,发现在波动问题中,非局部效应在时间上持续存在。对于初始条件含有不连续性的情况,PD解中不连续性的位置不会随时间发生改变。此外,初始条件的波长会影响非局部波动行为对非局部性的敏感性。利用PD波动方程解析解对弹性模型的δ收敛特性进行了研究,结果表明,弹性模型的δ收敛行为与扩散模型不同,表现出了关于δ的准周期特征。
(3)定义了非局部常数来定量衡量由影响函数导致的PD模型非局部性强弱程度,非局部常数越接近0,影响函数导致的非局部性越弱。通过解析分析证实了非局部常数与非局部行为之间的定量对应关系。将影响函数对PD数值模拟结果的影响拆分为非局部性和数值离散两方面,应用不同的数值算例分析了非局部常数与两者之间的关系。通过数值算例探究了影响函数对PD数值模拟结果的影响规律。结果表明,使用非局部常数小的影响函数倾向于获得接近经典解的结果,而使用非局部常数大的影响函数则更容易捕捉非均质因素的影响。
本文提出的PD解析解考虑了边界的影响,指明了局部边界条件和非局部边界条件之间的关系。此外,PD解析解还可以作为PD数值方法验证的参考依据。本文的研究成果有助于对PD非局部行为演化的理解,指导PD模型的参数选取。