近场动力学理论在脆性材料破坏研究中的应用现状

王玲玲 曹俊鑫 赵银霜 程想 孔德文

摘 要：近场动力学假定一定范围内的物质点之间存在非局部相互作用，通过空间积分重构物质点的运动方程，克服了传统有限元方法位移场连续性条件的局限，在分析强非线性不连续问题时具有无网格属性的数值优势，已成为研究脆性材料破坏的新兴理论。本文简要介绍了近场动力学的基本内容及其理论框架，总结了近场动力学理论在脆性材料准静态裂纹扩展、动态裂纹扩展及冲击失效研究方面的应用现状。

关键词：近场动力学;脆性材料;裂纹扩展;数值模拟;冲击失效

中图分类号：O346.1

文献标志码：A

由于本征脆性，混凝土、陶瓷、玻璃等脆性材料的破坏模式与破坏机理研究尤为重要。在达到极限承载力前，脆性材料经历从微裂纹产生到扩展的损伤过程，较低的抗拉强度使得裂纹扩展成为脆性材料的主要破坏模式。因此，许多学者力争能够准确预测脆性材料或结构的承载力以及相应的裂纹扩展过程与路径。目前，数值模拟是研究材料与结构内部裂纹产生与扩展问题的主要方法，如有限元方法（finite element method，FEM）[1]、扩展有限元方法（extended finite element method，XFEM）[2]和粒子方法[3]等。上述方法可以有效预测材料出现的大部分裂纹问题，但在复杂的裂纹问题（如裂纹合并、裂纹分支和任意三维裂纹问题）研究方面存在一定的局限性，而近场动力学在很大程度上克服了连续介质力学（computational continuum mechanics，CCM）的局限性，能够有效解决复杂的裂纹问题。

近场动力学（peridynamics，PD）的基本思想是由SILLING[4]提出的，它通过空间积分方程的求解来描述物质点的运动，可以看作是经典连续力学的一种非局部形式，因此PD中不再需要CCM中连续位移场的假设。即使在材料中出现不连续或裂纹，PD的控制方程也可以保持有效性。PD的另一个优点是无需使用额外的裂纹扩展准则，重新网格化方法和裂纹跟踪方法，就可以自发预测损伤过程及裂纹扩展情况。这些自然特征使近场动力学在研究和预测材料中裂纹或不连续性问题时拥有强大的适用性，甚至可以使用近场动力学很好地模拟复杂的裂纹问题。

SILLING[5]首先提出的理论是键型近场动力学 （bond-based peridynamics，BBPD），其中，两个物质点间的相互作用仅取决于该物质点间键的变形。这种假设使泊松比受到限制，平面应力情况下的泊松比固定为1/3，平面应变和三维情况下的泊松比固定为1/4。为了解决这个问题，SILLING[6]在此基础上提出了更通用的态型近场动力学（state-based peridynamics，SBPD）。态型近场动力学模型中两个物质点之间的相互作用不仅取决于物质点间键的变形，还取决于连接到这两个物质点的其他键的变形。态型近场动力学模型包括常规态型近场动力学和非常规态型近场动力学两种类型。常规态型近场动力学模型中力矢量的方向与键的方向平行，而非常规态型近场动力学模型中力矢量的方向与键的方向不要求平行[7]。相对于键型近场动力学，态型近场动力学的另一个优势是其PD表面效应明显小于BBPD[8]。由于近场动力学具有模拟各种复杂裂纹问题的强大能力，引起了研究人员越来越多的关注，并应用于脆性材料破坏行为的研究工作中。

1 PD理论基础

近场动力学理论通过求解空间积分方程而非有限元方法中的偏微分方程来描述物质点之间的相互作用。在空间域R上，假设材料中的两个物质点在给定近场范围δ内时存在非局部的相互作用。物质点之间的相互作用会引起键的伸长或缩短（见图1），物质点x和x′之间的相互作用力f可以由下式表示[4]：

2 近场动力学理论在材料破坏行为方面的应用

近场动力学理论已成为解决不连续问题的一种新兴理论，用于研究材料破坏行为，并取得了一系列较好的研究成果。本文总结了近场动力学理论在脆性材料准静态裂纹扩展，动态裂纹扩展及冲击破坏等方面的研究成果，阐明了近场动力学理论在研究脆性材料破坏行为方面的优势。

2.1 准静态裂纹扩展研究

脆性材料的准静态破坏过程一般包括裂纹的萌生、扩展直至材料破坏，近场动力学理论可以较好地预测准静态裂纹扩展路径与模式。KILIC等[10]学者基于近场动力学理论，引入自适应动态松弛法[11]，模拟了单轴拉伸荷载下带孔平板的裂纹扩展过程，数值计算收敛性好，验证了自适应动态松弛法预测材料损伤行为的适用性，裂纹产生和扩展的预测结果如图2所示。

黄丹等[12]学者将人工阻尼和系统失衡判断准则引入到模型中，实现了PD框架下静力学的定量计算，结合外力分级加载方法，采用近场动力学方法定量地分析了准静态情况下材料的变形计算和破坏问题。并通过典型算例验证了该算法的准确性，采用近场动力学方法实现了由准静态变形、裂纹萌生和扩展直至结构破坏的全过程模拟，准确预测了破坏荷载、裂纹萌生部位与扩展路径。此外，RABCZUK等[13]學者基于近场动力学理论，提出了适用于颗粒状材料与岩石断裂行为研究的dual-horizon peridynamics （DH-PD）公式。与传统有限元法相比，DH-PD可以模拟裂纹的自然路径，无需在裂纹表面添加任何准则来处理复杂的断裂形式（如裂纹分支和合并）。同时引入了人工阻尼，开展了岩石剪切试验和巴西圆盘劈裂试验过程的模拟，DH-PD的计算结果与试验结果吻合较好。裂纹扩展问题的静态解是分析结构破坏荷载和破坏机理的有效方法。LI等[14]学者建立了适用于桁架和张拉整体结构非线性静力分析的近场动力学模型，通过应变能密度的等价关系求解近场动力学参数，利用自适应动态松弛方法来获得稳态解，并采用典型算例验证了模型的有效性和准确性。

近场动力学的运动方程是基于时间和位移的空间积分，模型中的物质点是动态变化的，不适合直接进行静力学计算与模拟准静态破坏过程。但可以通过引入人工阻尼的动态松弛方法，使物质点在外力作用下达到平衡状态，充分利用了PD方法在处理不连续问题上的优势来成功解决准静态破坏问题。此外，BREITENFELD等[15]学者用隐式算法实现了在非常规态型近场动力学框架内的线弹性静力学计算。GALVANETTO等[16]学者耦合了有限元与近场动力学网格以获得静态解。NI等[17]学者基于近场动力学数值程序，提出了两种全新的隐式算法来研究裂纹扩展问题。以上研究成果表明基于近场动力学方法可以较好地实现准静态裂纹扩展形态及过程的数值模拟。

2.2 动态裂纹扩展研究

脆性材料的动态裂纹扩展过程短暂，比准静态情况复杂，很难通过试验观测捕捉其裂纹萌生部位或裂纹扩展路径，而近场动力学方法可以很好地重现动态裂纹的扩展过程。LEE等[18]学者建立近场动力学数值模型，研究了不同加载速率下裂纹的萌生和扩展现象，模拟了双轴拉伸荷载下含有缺口的脆性聚合物的裂纹扩展和分叉情况，破坏模式与试验观察结果一致，并从裂纹扩展角度和弹性应变能等方面揭示了脆性聚合物的破坏机理。ZHANG等[19]学者基于近场动力学理论，研究了动态载荷作用下环形试件的裂纹扩展问题，采用试验结果验证了近场动力学模拟结果的正确性，并开展了动态荷载作用下带孔洞环形试件裂纹扩展的模拟，讨论了孔洞对试件破坏路径的影响。通过峰值应力和破坏路径模拟与试验结果对比，证明了近场动力学理论分析动态荷载作用下岩石类材料破坏过程的有效性。Kalthoff-Winkler冲击试验（见图3）是经典的动态断裂问题，成为了验证模型适用性与算法准确性的典型算例[20]。

ZHOU等[21]学者基于非常规态基近场动力学理论，引入应力破坏准则，即当相互作用的物质点之间的平均应力满足应力破坏准则时，相互作用的物质点之间的键就会断裂。将该方法用于动载荷作用下脆性材料的裂纹扩展和分叉现象研究中。利用Kalthoff-Winkler试验验证该数值方法的有效性和准确性，并模拟了动态双轴载荷下脆性材料（Homalite-100）的单个裂纹扩展与分叉，讨论了几何特征与载荷条件的影响情况。

近场动力学的非局部特性也存在一些不足，如存在表面效应、边界条件不够灵活且计算效率低等缺点，但可以通过近场动力学与有限元耦合来弥补这些不足。MACEK等[22]学者利用桁架元素在FEM框架中实现了PD建模。KILIC等[23]学者提出引入重叠区域的耦合方法，在重叠区域中同时使用PD和FEM方程。SHOJAEI等[24]学者将近场动力学与有限点法（finite point-peridynamic method，FPM）耦合，开发了一种无网格方法来解决动态裂纹扩展问题，该方法将求解域分为三部分：PD离散化部分，FPM离散化部分及两种方法之间切换的过渡部分，以完全无网格的形式实现耦合。用PD描述裂纹存在或可能扩展的区域，用FPM描述需要较少计算量的区域。通过数值方法成功模拟了动态断裂动力学问题。PANCHADHARA等[25]学者基于近场动力学理论，利用近场动力学数据来估算应力强度因子，研究加载速率对动态裂纹扩展的影响情况，并通过准静态拉伸试验和Kalthoff-Winkler试验验证了该方法的可行性，结果表明随着加载速率的增大，裂纹扩展速度加快，裂纹分叉数量增多。YANG等[26]学者将常规态基近场动力学和有限元法耦合，研究了动态载荷下脆性材料的断裂现象。模型中的PD网格与FEM网格之间的作用力始终是相互平衡的，避免了耦合模型中过渡区域的复杂性，大幅降低了PD的表面效应。通过对比FEM结果和常规态型近场动力学结果，验证了所提出的耦合模型的正确性。此外，IMACHI等[27]学者基于常规态型近场动力学建立新的键失效模型，研究动态荷载下脆性材料的裂纹扩展现象。该模型将连接的键变为过渡键，用阻尼检验了过渡键的有效性，并减弱了应力强度因子的数值振荡，裂纹路径的模拟结果与Kalthoff-Winkle试验结果吻合良好。LAI等[28]学者建立了非常规态型近场动力学模型，引入改进的JH-2本构模型并构建了态型近场动力学公式，模拟了边缘撞击和玻璃平板落球冲击试验，模拟结果与试验结果基本一致。

经典的连续理论假设一个物质点仅与其直接相邻的物质点产生相互作用。当材料内部发生裂纹萌生和扩展时，原本连续的位移场与应力场将不再连续，裂纹尖端会产生奇异性。采用基于局部理论的有限元方法解决这些问题时，需要加入断裂准则来模拟裂纹的扩展与分叉，且在裂纹扩展后需要重新划分网格，对网格的依赖性较大[29]。为了克服上述不足，WAGNER等[30]学者提出了扩展有限元方法，但当出现位移不连续且涉及多个裂纹相互作用和分叉时，仍需引入裂纹扩展准则。使用XFEM模拟三维裂纹扩展时，计算过程会变得复杂。光滑粒子流体动力学（smoothed particle hydrodynamics，SPH）方法能有效解決上述问题[31]，但是这类分子动力学方法存在如计算时间较长、计算量较低的明显不足等。而近场动力学方法求解裂纹扩展等不连续问题时展现出了明显的优势，无需预设裂纹路径就可以模拟裂纹的自然萌生和扩展，并允许多条裂纹相互作用，且表现出很高的计算精度。此外，将有限元法与近场动力学耦合[32-33]，利用近场动力学离散裂纹扩展的区域，利用有限元法离散其他不需要大量计算的区域，大幅提高了模型的计算效率，同时降低了近场动力学模型的表面效应。

2.3 冲击破坏行为研究

脆性材料的冲击破坏过程快速且复杂，准确揭示冲击裂纹的萌生和扩展机理难度较大，而近场动力学方法在脆性材料冲击断裂过程研究方面适用性较强。WU等[34]学者基于近场动力学建立了动态损伤模型来研究混凝土的冲击破坏行为。该模型引入键的动态破坏及伸长率的损伤演化，建立了混凝土材料键的拉压破坏与伸长率之间的关系，并预报了不同情况下试样的裂纹扩展路径和裂纹扩展速度，与试验结果吻合较好。GUO等[35]学者建立了常规态型近场动力学模型，该模型引入线性的本构来表征材料的力学行为，采用该模型模拟了Kalthoff-Winkler低速冲击试验，与试验结果一致。王涵等[20]学者建立非常规态型近场动力学热黏塑性模型，基于Kalthoff-Winkler冲击试验，考虑靶板温度的变化，研究了不同冲击速度对裂纹扩展角度与扩展速度的影响规律。

另外，近场动力学在材料热冲击行为研究中也表现出较强的适用性，在一定程度上克服了理论和传统数值方法难以揭示热冲击裂纹萌生和扩展机理的弊端[36-39]。KILIC等[40-41]学者将键的热拉伸引入到本构力函数中，建立了非耦合键基近场动力学模型，成功地将该模型应用到断裂问题研究中。OTERKUS等[42]学者建立了键基与态基近场动力学的热力耦合模型。KILIC等[43]学者基于近场动力学，研究了不同淬火温度条件下含初始裂纹玻璃平板的裂纹扩展问题。WANG等[44]学者采用近场动力学热力耦合模型，研究了热冲击下脆性材料的裂纹扩展与分叉，并验证了模型的准确性，模拟了红外辐射加热下陶瓷圆盘的裂纹扩展试验，破坏路径与试验现象基本吻合。此外，还分析了热流密度、预制缺口长度及近场范围等因素对裂纹破坏形态的影响，系统地研究了裂纹尖端应变能、裂纹分叉角和裂纹扩展速度对裂纹扩展的影响情况。XU等[45]学者建立了近场动力学模型，以预测热冲击下玻璃薄板的裂纹扩展。WANG等[46]学者在此基础上考虑了热力荷载下物质点相对位移对热弹性刚度与热传导的影响，基于键基近场动力学提出了改进的热力耦合模型，用以研究热冲击作用下脆性材料的损伤机理，模拟了淬火条件下陶瓷材料的裂纹扩展试验，数值结果（如图4与图5）与试验现象吻合良好。

D′ANTUONO等[47]學者建立常规状态型近场动力学热力耦合模型，研究了热冲击下陶瓷薄板与厚板的断裂行为， 得到二维有序平行裂纹和三维柱状蜂窝裂纹的模拟结果。GIANNAKEAS等[48]学者基于键基近场动力学，模拟了热冲击下陶瓷材料的裂纹扩展，预测热冲击裂纹的损伤演化，评价极端温度变化下陶瓷材料的热冲击响应，模拟裂纹成核的动态过程。但有限元与近场动力学的耦合模型不适合模拟热冲击裂纹扩展过程中温度场不连续变化现象。WANG等[49]学者基于常规态型近场动力学，考虑应变软化特性，引入微观导热系数，构建微观与宏观几何条件之间的关系，建立弱耦合热弹性模型，模拟了淬火条件下陶瓷板的裂纹扩展问题，裂纹路径与试验现象吻合较好，结果表明试件尺寸和淬火温度对应变能和裂纹扩展速度影响显著。

基于含有部分微分运动控制方程的数值方法可以成功预报冲击裂纹的分叉和扩展[50-53]，但如果不引入破坏准则就不能通过偏微分方程处理裂纹引起的位移不连续问题。近场动力学以积分形式重构了连续体的运动方程，允许位移场不连续，在处理不连续性问题方面具有强大数值优势。利用近场动力学方法能更加直观地获得冲击荷载作用下脆性材料裂纹扩展的整个过程，为不同工况下脆性材料裂纹扩展及分叉问题研究提供新途径。

3 结语

作为一种非局部理论，近场动力学避免了传统方法在解决不连续问题时存在的奇异性，成为研究脆性材料破坏行为的一种新兴理论。本文在简单介绍了近场动力学理论的基础上，综述了近场动力学理论在脆性材料准静态裂纹扩展、动态裂纹扩展、冲击损伤等方面的研究近展，总结如下：

1）与经典连续介质力学不同，近场动力学采用积分方程而不是位移分量的导数来表示。近场动力学允许材料内部自发的裂纹萌生与扩展路径自由，无需引入额外的裂纹扩展准则。在近场动力学理论中，内力是由连续体中任意两个物质点之间的本构关系来表示的，损伤是本构模型的一部分。

2）通过引入人工阻尼的动态松弛法，使材料点在外力作用下达到平衡状态，也可通过隐式算法得到静态解，证明了近场动力学方法模拟准静态裂纹扩展形态及过程的适用性。

3）近场动力学以积分形式重构连续介质的运动方程，允许位移场不连续，在处理不连续问题时有数值优势。采用近场动力学模型可以更好地展现脆性材料裂纹扩展的全过程，弥补试验中难以观测的不足。

4）近场动力学与有限元的耦合降低了近场动力学模型的表面效应，降低了计算成本，提高了计算效率。

参考文献：

[1]KEUI O， KIYOTSUGU O. FEM analysis of crack closure and delay effect in fatigue crack growth under variable amplitude loading[J]. Engineering Fracture Mechanics， 1977， 9（2）： 471-480.

[2] GREGOIRE D， MAIGRE H， RETHORE J， et al. Dynamic crack propagation under mixed-mode loading comparison between experiments and X-FEM simulations[J]. International Journal of Solids and Structures， 2007， 44（20）： 6517-6534.

[3] RABCZUK T， SONG J H， BELYTSCHKO T. Simulations of instability in dynamic fracture by the cracking particles method[J]. Engineering Fracture Mechanics， 2009， 76（6）： 730-741.

[4] SILLING S A. Reformulation of elasticity theory for discontinuities and long-range forces[J]. Journal of the Mechanics Physics of Solids， 2000， 48（1）： 175-209.

[5] SILLING S A， ASKARI E. A meshfree method based on the peridynamic model of solid mechanics[J]. Computers & Structures， 2005， 83（17）： 1526-1535.

[6] SILLING S A， EPTON M， WECKNER O， et al. Peridynamic states and constitutive modeling[J]. Journal of Elasticity， 2007， 88（2）： 151-184.

[7] YANG D， HE X Q， YI S H， et al. An improved ordinary state-based peridynamic model for cohesive crack growth in quasi-brittle materials[J]. International Journal of Mechanical Sciences， 2019， 153： 402-415.

[8] LE Q V， BOBARU F. Surface corrections for peridynamic models in elasticity and fracture[J]. Computational Mechanics， 2008，61（4）： 499-518.

[9] HA Y D， BOBARU F. Studies of dynamic crack propagation and crack branching with peridynamics[J]. International Journal of Fracture， 2010， 162（1）： 229-244.

[10]KILIC B， MADENCI E. An adaptive dynamic relaxation method for quasi-static simulations using the peridynamic theory[J]. Theoretical Applied Fracture Mechanics， 2010， 53（3）： 194-204.

[11]UNDERWOOD P. Dynamic relaxation[J]. Computational Methods for Transient Analysis， 1983， 1： 245-265.

[12]黃丹， 卢广达， 章青. 准静态变形破坏的近场动力学分析[J]. 计算力学学报， 2016， 33（5）： 657-662.

[13]RABCZUK T， REN H. A peridynamics formulation for quasi-static fracture and contact in rock[J]. Engineering Geology， 2017， 225： 42-48.

[14]LI H， ZHANG H， ZHENG Y G， et al. A peridynamic model for the nonlinear static analysis of truss and tensegrity structures[J]. Computational Mechanics， 2016， 57（5）： 843-858.

[15]BREITENFELD M S， GEUBELLE P H， WECKNER O， et al. Non-ordinary state-based peridynamic analysis of stationary crack problems[J]. Computer Methods in Applied Mechanics and Engineering， 2014， 272： 233-250.

[16]GALVANETTO U， TEO M， SHOJAEI A， et al. An effective way to couple FEM meshes and peridynamics grids for the solution of static equilibrium problems[J]. Mechanics Research Communications， 2016， 76： 41-47.

[17]NI T， ZACCARIOTTO M， ZHU Q Z， et al. Static solution of crack propagation problems in Peridynamics[J]. Computer Methods in Applied Mechanics Engineering， 2018， 346： 126-151.

[18]LEE J， HONG J. Dynamic crack branching and curving in brittle polymers[J]. International Journal of Solids Structures， 2016， 100： 332-340.

[19]ZHANG Y N， DENG H W， DENG J R， et al. Peridynamics simulation of crack propagation of ring-shaped specimen like rock under dynamic loading[J]. International Journal of Rock Mechanics and Mining Sciences， 2019， 123： 104093.

[20]王涵， 黄丹， 徐业鹏， 等. 非常规态型近场动力学热黏塑性模型及其应用[J]. 力学学报， 2018， 50（4）： 810-819.

[21]ZHOU X P， WANG Y T， QIAN Q H. Numerical simulation of crack curving and branching in brittle materials under dynamic loads using the extended non-ordinary state-based peridynamics[J]. European Journal of Mechanics/A Solids， 2016， 60： 277-299.

[22]MACEK R W， SILLING S A. Peridynamics via finite element analysis[J]. Finite Elements in Analysis and Design， 2007， 43（15）： 1169-1178.

[23]KILIC B， MADENCI E. Coupling of peridynamic theory and the finite element method[J]. Journal of Mechanics of Materials and Structures， 2010， 5（5）： 707-733.

[24]SHOJAEI A， MUDRIC T， ZACCARIOTTO M， et al. A coupled meshless finite point/Peridynamic method for 2D dynamic fracture analysis[J]. International Journal of Mechanical Sciences， 2016， 119： 419-431.

[25]PANCHADHARA R， GORDON P. Application of peridynamic stress intensity factors to dynamic fracture initiation and propagation[J]. International Journal of Fracture， 2016， 201（1）： 81-96.

[26]YANG D， HE X Q， YI S H， et al. Coupling of peridynamics with finite elements for brittle crack propagation problems[J]. Theoretical and Applied Fracture Mechanics， 2020， 107： 102505.

[27]IMACHI M， TANAKA S， BUI T Q， et al. A computational approach based on ordinary state-based peridynamics with new transition bond for dynamic fracture analysis[J]. Engineering Fracture Mechanics， 2019， 206： 359-374.

[28]LAI X， LIU L S， LI S F， et al. A non-ordinary state-based peridynamics modeling of fractures in quasi-brittle materials[J]. International Journal of Impact Engineering， 2018， 111： 130-146.

[29]MOES N， DOLBOW J， BELYTSCHKO T. A finite element method for crack growth without remeshing[J]. International Journal for Numerical Methods in Engineering， 1999， 46（1）： 131-150.

[30]WAGNER G J， MOES N， LIU W K， et al. The extended finite element method for rigid particles in Stokes flow[J]. International Journal for Numerical Methods in Engineering， 2001， 51（3）： 293-313.

[31]GINGOLD R A， MONAGHAN J J. Smoothed particle hydrodynamics： theory and application to non-spherical stars[J]. Monthly Notices of the Royal Astronomical Society， 1977， 181（3）： 375-389.

[32]BIE Y H， CUI X Y， LI Z C. A coupling approach of state-based peridynamics with node-based mothed finite element method[J].Computer Methods in Applied Mechanics and Engineering， 2018， 331： 675-700.

[33]ZACCARIOTTO M， MUDRIC T， TOMASI D， et al. Coupling of fem meshes with peridynamic grids[J].Computer Methods in Applied Mechanics and Engineering， 2018， 330： 471-497.

[34]WU L W， HUANG D， XU Y P， et al. A Rate-dependent dynamic damage model in peridynamics for concrete under impact loading[J]. International Journal of Damage Mechanics， 2020， 29（7）： 1035-1058.

[35]GUO J S， GAO W C. Study of the Kalthoff-winkler experiment using an ordinary state-based peridynamic model under low velocity impact[J]. Advances in Mechanical Engineering， 2019， 11（5）： 1-11.

[36]BOURDIN B， MARIGO J J， MAURINI C， et al. Morphogenesis and propagation of complex cracks induced by thermal shocks[J]. Physical Review Letters， 2014， 112（1）： 014301.

[37]LI J， SONG F， JIANG C. A non-local approach to crack process modeling in ceramic materials subjected to thermal shock[J]. Engineering Fracture Mechanics， 2015， 133： 85-98.

[38]HUANG X， TANG S B， TANG C A， et al. Numerical simulation of cracking behavior in artificially designed rock models subjected to heating from a central borehole[J]. International Journal of Rock Mechanics and Mining Sciences， 2017， 98： 191-202.

[39]BADNAVA H， MSEKH M A， ETEMADI E， et al. An h-adaptive thermo-mechanical phase field model for fracture[J]. Finite Elements in Analysis and Design， 2018， 138： 31-47.

[40]KILIC B， MADENCI E. Prediction of crack paths in a quenched glass plate by using peridynamic theory[J]. International Journal of Fracture， 2009， 156（2）： 165-177.

[41]KILIC B， MADENCI E. Peridynamic theory for thermomechanical analysis[J]. IEEE Transactions on Advanced Packaging， 2010， 33（1）： 97-105.

[42]OTERKUS S， MADENCI E， AGWAI A. Fully coupled peridynamic thermomechanics[J]. Journal of the Mechanics and Physics of Solids， 2014， 64（1）： 1-23.

[43]KILIC B， MADENCI E. Structural stability and failure analysis using peridynamic theory[J]. International Journal of Non-linear Mechanics， 2019， 44（8）： 845-854.

[44]WANG Y T， ZHOU X P， KOU M M. Numerical studies on thermal shock crack branching instability in brittle solids[J]. Engineering Fracture Mechanics， 2018， 204： 157-184.

[45]XU Z P， ZHANG G F， CHEN Z G， et al. Elastic vortices and thermally-driven cracks in brittle materials with peridynamics[J]. International Journal of Fracture， 2018， 209（1）： 203-222.

[46]WANG Y T， ZHOU X P， KOU M M. An improved coupled thermo-mechanic bond-based peridynamic model for cracking behaviors in brittle solids subjected to thermal shocks[J]. European Journal of Mechanics/A Solids， 2019， 73： 282-305.

[47]D′ANTUONO P， MORANDINI M. Thermal shock response via weakly coupled peridynamic thermo-mechanics[J]. International Journal of Solids and Structures， 2017， 129： 74-89.

[48]GIANNAKEAS L N， PAPATHANASIOU T K， BAHAI H. Simulation of thermal shock cracking in ceramics using bond-based peridynamics and FEM[J]. Journal of the European Ceramic Society， 2018， 38（8）： 3037-3048.

[49]WANG Y T， ZHOU X P， ZHANG T. Size effect of thermal shock crack patterns in ceramics： insights from a nonlocal numerical approach[J]. Mechanics of Materials， 2019， 137： 103-113.

[50]HUGHES T， TOYLOR R， SACKMAN J， et al. Finite element method for a class of contact-impact problems[J]. Computer Methods in Applied Mechanics and Engineering， 1976， 8（3）： 249-276.

[51]CHAUDHARY A B， BATHE K J. A solution method for static and dynamic analysis of three-dimensional contact problems with friction[J]. Computers & Structures， 1986， 24（6）： 855-873.

[52]MENOUILLARD T， BELYTSCHKO T. Analysis and computations of oscillating crack propagation in a heated strip[J]. International Journal of Fracture， 2010， 167（1）： 57-70.

[53]TANG S B， ZHANG H， TANG C A， et al. Numerical model for the cracking behavior of heterogeneous brittle solids subjected to thermal shock[J]. International Journal of Solids and Structures， 2016，80： 520-531.

（责任编辑：于慧梅）

Application Situation of Peridynamics Theory in the Study

of the Failure of Brittle Materials

WANG Lingling， CAO Junxin*， ZHAO Yinshuang， CHENG Xiang， KONG Dewen

（College of Civil Engineering， Guizhou University， Guiyang 550025， China）

Abstract：

Peridynamics assumes that there are non-local interactions between material points within a certain range， and the equations of motion of the material points are reconstructed through integral operator. The limitation of the continuity condition of displacement field of traditional finite element was eliminated. The numerical advantage of the meshless property of peridynamic in analysing strong nonlinear discontinuity problems makes it become an emerging theory for studying the failure of brittle materials. The basic content and theoretical framework of peridynamics were briefly introduced. The application situation of peridynamics theory in material failure research was summarized in terms of quasi-static crack propagation， dynamic crack propagation and impact failure.

Key words：

peridynamic; brittle material; crack propagation; numerical simulation; impact failure

收稿日期：2021-03-04

基金項目：国家自然科学基金资助项目（51968009）;贵州省科技计划资助项目（[2020]1Y244）;贵州大学资助项目（贵大培育[2019]28号）

作者简介：王玲玲（1985—），女，副教授，博士，研究方向：复合材料力学性能研究，E-mail：llwang1@gzu.edu.cn.

通讯作者：曹俊鑫，E-mail：1115006391@qq.com.