Research Article  Open Access
Peridynamic Model for the Numerical Simulation of Anchor Bolt Pullout in Concrete
Abstract
Predictive simulation of anchor pullout from concrete structures is not only a serious problem in structural mechanics but also very important in structural design safety. In the finite element method (FEM), the crack paths or the points of crack initiation usually need to be assumed in advance. Otherwise, some special crack growth treatment or adaptive remeshing algorithm is normally used. In this paper, an extended peridynamic method was introduced to avoid the difficulties found in FEM, and its application on anchor bolt pullout in plain concrete is studied. In the analysis, the interaction between the anchor bolt and concrete is represented by a modified shortrange force and an extended bondlevel model for concrete is developed. Numerical analysis results indicate that the peak pullout load obtained and the crack branching of the anchoring system agreed well with the experimental investigations.
1. Introduction
Anchor bolts are very important components of load transfer in a wide range of civil engineering structures such as dams, nuclear power plants, highways, and bridges. A better understanding of the pullout behavior of anchor bolts can contribute not only to the optimization of the design of the anchor system, but also to the improvement of the durability and stability of a structure. Therefore, the pullout behavior of anchor bolts in concrete structures has become a major concern in the past three decades and a lot of experimental studies have been performed [1–6]. Likewise, various numerical methods have been adopted to analyze the failure mechanism and the progressive damage of anchor bolts in concrete structures.
Among previous works, most of the researchers focused on the finite element method (FEM) [7–11]. One of the early methods [7] includes modeling of concrete as the fournode cracked element. In this approach, the tensile fracture behavior of the concrete for a solid body containing an internal discontinuous surface is formulated by deriving an incremental formulation. As a result, the crack is distributed in each grid and the localization of microcracking cannot be obtained. In order to better simulate the modeI fracture, Alfaiate et al. [8] utilized the interface elements which are inserted along the interelement boundaries among the concrete elements to model the cracking path. Although multiple cracks can be obtained and no special remeshing technique is required with such approach, the crack direction is still restricted and laddershaped crack paths are formed.
Etse [9] predicted the distribution of equivalent fracture strain at peak load of the anchor system by adopting a fracture energybased plasticity formulation. The propagation of crack is described in terms of an equivalent plastic softening process, but still it is difficult to obtain the final cracking pattern by this approach. Xu et al. [10] simulated the crack patterns and mechanical behavior of the anchor bolt pullout in concrete by using the Mohr–Coulomb criterion with tension cutoff. In this approach, the heterogeneity of the concrete material is modeled by randomly assigning strength and elastic modulus to the elements according to Weibull’s distribution and the ongoing cracking process is represented by groups and alignments of failed elements. Feenstra [11] used the smeared crack approach to study twodimensional pullout problems. In this approach, it is assumed that a crack can be distributed over a special band in the model and its influence on the mechanical behavior of a material can be represented by adjusting the constitutive matrix irrespective of the real displacement discontinuities within the band. It is obvious that these aforementioned approaches had overcome some deficiencies of crack propagation problems, but still these approaches require remeshing or redefining of the geometry to model the progressive crack growth. Furthermore, the accuracy of results heavily relies on the complex adaptive meshing algorithms [12]. Therefore, the low efficiency and accuracy of simulation are still a major problem [12–14].
Recent developments of meshfree (or meshless) methods such as diffuse element method (DEM), material point method (MPM), and element free Galerkin method (EFGM) are invented to circumvent the meshdependence problem and relieve the volumetric locking for suitable choice of support size of shape function [15, 16]. Soparat and Nanakorn [14] and Coetzee et al. [17] used meshfree methods for pullout problems. Though many discussions on the advantages of the meshfree methods have been reported [15–18], it should still be noted that the treatment of essential boundary conditions is not straightforward as the conventional FEM and the use of shape functions of any desired order of continuity may lead to computational difficulties and complexity in deriving the coefficients of the stiffness matrix [18, 19]. Other contributions can be found in the literature regarding anchors’ modeling using the microplane model [20, 21].
Silling and Askari [22, 23] proposed peridynamics as an alternative technique for the solution of crack related problems where direct displacement is used in the formulations instead of displacement derivatives, and deformation continuity is not based on assumptions. Mainly, bonds containing constitutive information of the materials are used to reproduce the nonlocal interacting forces between particles over a certain distance. Moreover, in contrast to the partial differential equations used in the classical formulation, this theory uses spatial integral equations which permit spontaneous crack occurring at multiple sites and freely extend along an arbitrary path without extra remedial techniques [24]. Therefore, the peridynamic theory has obvious advantages in handling cracks problems and predicting the progressive failure process in solid mechanics. Numerous achievements have been obtained using this model during the past few years including deformation of onedimensional bar [25, 26], progressive damage of composite laminates [27, 28], damage and fracture of membranes and fibers [29], static and dynamic fracture of plain and reinforced concrete structure [30–34], dynamic fracture in functionally graded materials [35], and coupling with classical continuum mechanics [36].
In this paper, an extended bondlevel model for concrete is proposed and the interaction between the anchor bolt and concrete is represented by a modified shortrange force. The analysis of anchor bolt pullout problem is carried out using the bondbased peridynamic (PD) model. Moreover, the capabilities of the improved numerical method to capture the progressive damage process and the extreme load of the anchor bolt are validated. Finally, a parametric study was performed to investigate the influence of the size of the horizon and the embedded length. Comparison of the experiments and simulations with those in the literature is also carried out.
2. Introduction of the Peridynamic Theory (PD)
The PD theory may be viewed as a special version of particle method or meshfree method. It is based on assumptions that an object possesses a spatial domain modeled as a discrete set of particles and each particle owns a subregion within a certain radius called the material horizon as shown in Figure 1. The peridynamic equation of motion at any time is given by Silling [22] as follows:where denotes the mass density; is the prescribed external body force density; and are the displacement vector and acceleration vector, respectively. Also, is a pairwise force vector in the peridynamic bond that represents the nonlocal interactions between the particle and the rest of the particles. Unlike MD, in peridynamics, the magnitude of depends upon the initial reference configuration and the state of the bond, that is, the relative position vector and the relative displacement vector , as follows:
The direction of the pairwise force vector can be expressed as ; thus, a general expression for the peridynamic force function can be written aswhere is a scalarvalued force function and . For a microelastic material [22], another expression form of the pairwise force vector can be derived from a pairwise potential function such thatThe general form of the linear microelastic potential is obtained as follows:where represents the bond elastic stiffness known as “micromodulus” function, is the scalarvalued Boolean influence function, is a constant that depends on , and denotes the relative elongation of a bond and it is given asIf the relative elongation , then the pairwise force between the particles does not exist. According to the above expressions, the corresponding pairwise force becomeswhere is a historydependent function which depends on the value of and it can be written as denotes the critical stretch of the bond and can be obtained by mathematical derivation and by processing the classical fracture parameters. The value of local damage at point within a peridynamic material can be denoted by the percent of broken bonds aswhere can be viewed as a matrix representing the damage of material point. Silling [22] introduced an original constitutive model for quasibrittle materials; the influence function is assumed to be 1.0, and the micromodulus function can be correlated with the classical elastic constants through the equivalent of elastic strain energy in elasticity and PD theory. For 2D plane stress problems, as illustrated in Figure 2, from the conventional theory of linear elasticity, the strain energy density due to a uniform principal strain state and uniform shear strain state can be calculated aswhere is the elastic modulus, is the Poisson ratio, and denotes the strain. On the other hand, from the twodimensional peridynamic theory, , the strain energy density can be calculated as [34]where denotes the thickness. By solving the equations and , the constant can be obtained asThe critical stretch is associated with the fracture energy and the bond would break when the elongation goes beyond the critical value . In 2D plane stress conditions, the work required to break all the bonds per unit fracture area can be derived as in [23]:where is the distance between the point and the crack surface as shown in Figure 3 and can be obtained by substituting into (13) as follows:
3. Peridynamic Model for Anchor Bolts
3.1. Peridynamic Model for Concrete and Steel
As peridynamics do not need the continuous displacement field, there are no concepts of stress or strain required in the model. Thus, the constitutive model is defined through the relationship between the bond stretch and the pairwise force among material particles or, in other words, the material damage is introduced at the bond level. Analogous to the softening function proposed by Gerstle et al. [37], a simplified softening function is introduced in the present work and the bond force and bond stretch relationship are illustrated in Figure 4.
(a) Concrete
(b) Steel
As shown in Figure 4(a), the yield compressive and tensile stretch limits of concrete are denoted as and , respectively. For small bond stretches, , the bond remains in the linear elastic range. However, if is less than the compressive stretch limit , the pairwise force of concrete remains constant. On the other hand, an abrupt drop will happen in pairwise force if reaches the tensile stretch limit and the pairwise force will become until reaches the critical stretch . A brittle fracture will occur at this point and the pairwise force will drop to zero. As steel may yield in tension as shown in Figure 4(b), the pairwise force of steel remains constant if reaches the tensile stretch limit of steel . , , and can be given asTo provide required fracture energy in tension, according to (13),Substituting (11a) and (11b) into (13), the relationship between and can be computed asThe value of here is 0.5 as advised by Cheng et al. [35], so value can be obtained as
3.2. Interaction between Steel and Concrete
Stress is transferred mainly by adhesion, mechanical interaction, and friction between steel and concrete. Adhesion comes from chemical bonding and stresses are generated during curing of concrete and experimentally it is very tough to measure them. Due to the limited experimental information, it is difficult to determine the friction coefficient. As stresses due to adhesion and friction are relatively small, therefore, only mechanical interaction is considered in the present simulation. In this paper, we introduce a shortrange force [23, 24, 38] to reproduce the interaction between the concrete and anchor bolt aswhere denotes the distance between two types of particles and the constant denotes the average micromodulus which is assumed to be
3.3. Discretization and Numerical Implementation
Peridynamic equations of motion (see (1)) could be solved by utilizing a numerical approximation method which involves the discretization of the reference configuration into particles with a certain volume. So, the integrals in (1) can be replaced by the finite sums:where denotes the acceleration of the point at time step , is the total number of particles within the horizon of the point , and are the pairwise forces and the body force density at time step , and denotes the equivalent calculation volume (see [38]). To solve quasistatic problems by applying peridynamic equations of motion, Kilic and Madenci [39] introduced artificial damping to attain a steadystate solution. The Adaptive Dynamic Relaxation (ADR) scheme proposed by Underwood [40] is used to determine the most effective damping coefficient at each time step , as follows: where is the resultant force density vector; is the damping coefficient at the th iteration and is the modified density at the point and they can be given as
Here, is the time step size, is a unit vector along , , or , and is the diagonal stiffness matrix of the system, given as
Finally, with the assumptions that and , velocities and displacement at point for the next time step can be obtained by central difference explicit integration as follows:
Both the numerical algorithm and constitutive modeling are implemented in Fortran90 language based on Visual Studio using an inhouse peridynamic code.
4. Peridynamic Results for Anchor Bolt Pullout in Concrete
4.1. Problem Setup
Because of the distinctive advantages in solving crack propagation problems, the peridynamic method is adopted in the present study to model the anchor bolt pullout in concrete. The anchor bolt pullout experiment conducted by Vervuurt et al. [6] is considered here. Table 1 lists the geometrical parameters of the experiment. Figure 5 shows the corresponding geometrical features of the prediction model.

and are the width and length of the concrete block as shown in Figure 5. is the support span and it is the distance from the edge of the anchor to the support and is the embedded depth. The breadth and thickness of the anchor bolt are and , respectively. All the concrete blocks and anchor bolts are 0.1 m thick, which is very small as compared to other dimensions, so the model is treated as a plane stress problem in the analysis. In most of the past research works [9, 10, 14], steel bolt was also not directly modeled and there were no experimental values of friction coefficients. Pullout load is modeled by controlling the vertical displacement of the top surface of the bolt head and for the convenience of simulations only the anchor head is assumed to be in contact with the concrete. As a result, both the friction resistance caused by the pullout movement and the vertical deformation of the anchor bar were ignored. In the present study, both the plain concrete and the steel anchor are modeled as discrete particles of corresponding volumes. In addition, a shortrange force is used to model the frictionless contact between concrete and the anchor bolt. To reduce the computation cost and modeling difficulty, a constant inplane grid number of 240 is used for each prediction model so that the grid spacing of three groups is set to be 1.25 mm, 2.5 mm, and 3.75 mm, respectively. Representative of quasistatic loading, a very low velocity of m/s was applied on the top particles of the anchor bolt. The material properties are depicted in Table 2 and the corresponding numerical models are shown in Figure 6.

4.2. Results and Discussions
To investigate the influence of the size of material horizon on the loaddisplacement response of the anchor bolt, three different material horizon sizes (i.e., , , and ) were used. Figure 7 shows the relationship between the vertical displacement of the top particles and the pullout force in Group 2. It is observed that the results obtained with different material horizons approximately give the same peak load. For the larger material horizon, a flatter initial slope and greater peak displacement are witnessed. As the material horizon is related to critical stretch for the proposed concrete model (see (18)), a lower critical stretch shows that the bond is more brittle than the material. As shown in Figure 8, the curve is more close to the experiment results by Vervuurt et al. [6] when is adopted; therefore, the default material horizon is set to be three times the grid spacing in the following simulation.
Figure 8 also shows the lattice model, FE, and EFG results by Vervuurt et al. [6], Feenstra [11], and Soparat and Nanakorn [14], respectively. The peridynamic model shows very good results in the upper section of the experimental curve for loaddisplacement behavior. Figure 9 shows loaddisplacement curves for all groups. The predicted peak loads and experimental values comparison is given in Table 3. Embedded depth has an influence on the pullout mechanism; however, the peak load does not increase proportionally with the embedded depth which shows that size has effect on pullout loads [6].

The values of the predicted peak loads in the present study are a little higher than the experimental values which may be because of the higher Poisson’s ratio; 0.33 instead of normal 0.3, based on the peridynamic theory [22, 23] which can be seen in Figure 8 and Table 3.
Although the precision obtained by the current model is not very high as compared to EFGM in peak load prediction, it was based on many simplifications [14]. (a) Half of the model was considered due to symmetry. (b) Steel bolt was not directly modeled. (c) Supports were treated as points. (d) Crack propagation was assumed from the upper edge corner of the bolt head, among others. These simplifications may affect the reliability of simulation results. In peridynamics, both support and specimen are modeled as particles, so simulations are more reliable than in [14]. Cracks can appear spontaneously as the material response includes damage without the need for any special criterion. Crack trajectories of the present model are compared with the EFGM results by Soparat and Nanakorn [14] for different groups as shown in Figure 10. The crack branching can be easily captured in PD theory while there is no localized branching in EFGM.
(a) Group 1 (mm, displacement × 100)
(b) Group 2 (mm, displacement × 100)
(c) Group 3 (mm, displacement × 100)
The crack pattern from experimental results is different from most of the previous numerical methods results [6, 11, 14]. The reasons can be the assumption of a 2D model instead of the real 3D complex fracture model for simulations, heterogeneous nature of concrete [6, 10], and experimental errors in the tests [10]. Figure 11 compares the failure mode of the present study with the experimental observations [6]. It is worth noting that both the crack direction and the crack branching obtained in the present study show good resemblance with the experimental result by Vervuurt et al. [6].
(a) Group 1 (mm, displacement × 50)
(b) Group 2 (mm, displacement × 50)
(c) Group 3 (mm, displacement × 50)
5. Conclusions
In this paper, the problem of anchor bolt pullout in plain concrete was investigated with an improved peridynamic model. An extended constitutive model was used to refine the behaviors of concrete material simulation. A shortrange force was introduced to simulate anchor bolt and concrete interaction in 2D. The numerical discretization and iteration algorithms were implemented with an inhouse peridynamic FORTRAN code. It was observed that the crack propagation of concrete was exposed in more detail by the proposed approach as compared to conventional FEM or EFGM. Compared with the results from the literatures and experiments, it can be concluded that the extreme failure load and the final failure mode of the anchor bolt by analysis of the peridynamic approach match well those of the experimental observations. From all the results and comparisons, this approach was proved to be a promising method for solving the problem of anchor bolt pullout in plain concrete.
Nomenclature
:  The mass density 
:  The material horizon 
:  The prescribed external body force density 
:  The displacement vector 
:  The relative position vector 
:  The relative displacement vector 
:  The pairwise force function 
:  The pairwise potential function 
:  The influence function 
:  The micromodulus function 
:  The fracture energy 
:  The damage of material point at time 
:  The tensile stretch limit of steel 
:  The compressive stretch limit of concrete 
:  The tensile stretch limit of concrete 
:  The adjustment coefficients 
:  The equivalent calculation volume 
:  The resultant force density vector 
:  The damping coefficient at the th iteration 
:  The modified density at the point 
:  The time step size 
:  The diagonal stiffness matrix. 
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
References
 W. C. Stone and N. J. Carino, “Deformation and Failure in LargeScale Pullout Tests,” ACI Journal Proceedings, vol. 80, pp. 501–513, 1983. View at: Publisher Site  Google Scholar
 P. Bocca, “The application of pullout test to high strength concrete estimation,” Matériaux et Constructions, vol. 17, no. 3, pp. 211–216, 1984. View at: Publisher Site  Google Scholar
 R. Ballarini, S. P. Shah, and L. M. Keer, “Failure Characteristics of Short Anchor Bolts Embedded in a Brittle Material,” Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences, vol. 404, no. 35, 1986. View at: Publisher Site  Google Scholar
 J. Hashimoto and K. Takiguchi, “Experimental study on pullout strength of anchor bolt with an embedment depth of 30 mm in concrete under high temperature,” Nuclear Engineering and Design, vol. 229, no. 23, pp. 151–163, 2004. View at: Publisher Site  Google Scholar
 S.Y. Kim, C.S. Yu, and Y.S. Yoon, “Sleevetype expansion anchor behavior in cracked and uncracked concrete,” Nuclear Engineering and Design, vol. 228, no. 13, pp. 273–281, 2004. View at: Publisher Site  Google Scholar
 A. Vervuurt, J. G. M. Van Mier, and E. Schlangen, “Analyses of anchor pullout in concrete,” Materials and Structures, vol. 27, no. 5, pp. 251–259, 1994. View at: Publisher Site  Google Scholar
 A. Ali, “FEM analysis of concrete structures subjected to modeI and mixedmode loading conditions,” Computers & Structures, vol. 61, no. 6, pp. 1043–1055, 1996. View at: Publisher Site  Google Scholar
 J. Alfaiate, E. B. Pires, and J. A. C. Martins, “A finite element analysis of nonprescribed crack propagation in concrete,” Computers & Structures, vol. 63, no. 1, pp. 17–26, 1997. View at: Publisher Site  Google Scholar
 G. Etse, “Finite element analysis of failure response behavior of anchor bolts in concrete,” Nuclear Engineering and Design, vol. 179, no. 2, pp. 245–252, 1998. View at: Publisher Site  Google Scholar
 C. Xu, C. Heping, L. Bin, and Z. Fangfang, “Modeling of anchor bolt pullout in concrete based on a heterogeneous assumption,” Nuclear Engineering and Design, vol. 241, no. 5, pp. 1345–1351, 2011. View at: Publisher Site  Google Scholar
 P. H. Feenstra, Computational Aspects of Biaxial Stress in Plain and Reinforced Concrete, Delft University of Technology, Delft, The Netherlands, 1993.
 P. O. Bouchard, F. Bay, Y. Chastel, and I. Tovena, “Crack propagation modelling using an advanced remeshing technique,” Computer Methods Applied Mechanics and Engineering, vol. 189, no. 3, pp. 723–742, 2000. View at: Publisher Site  Google Scholar
 M. Jirásek, “Comparative study on finite elements with embedded discontinuities,” Computer Methods Applied Mechanics and Engineering, vol. 188, no. 1, pp. 307–330, 2000. View at: Publisher Site  Google Scholar
 P. Soparat and P. Nanakorn, “Analysis of anchor bolt pullout in concrete by the elementfree Galerkin method,” Engineering Structures, vol. 30, no. 12, pp. 3574–3586, 2008. View at: Publisher Site  Google Scholar
 T. Belytschko, Y. Krongauz, D. Organ, M. Fleming, and P. Krysl, “Meshless methods: An overview and recent developments,” Computer Methods Applied Mechanics and Engineering, vol. 139, no. 14, pp. 3–47, 1996. View at: Publisher Site  Google Scholar
 S. Li and W. K. Liu, “Meshfree and particle methods and their applications,” Applied Mechanics Reviews, vol. 55, no. 1, pp. 1–34, 2002. View at: Publisher Site  Google Scholar
 C. J. Coetzee, P. A. Vermeer, and A. H. Basson, “The modelling of anchors using the material point method,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 29, no. 9, pp. 879–895, 2005. View at: Publisher Site  Google Scholar
 A. Charles and H. Claire, “The use of meshless methods in geotechnics,” in Proceedings of the 1st International Symposium on Comptational Geomechnics (COMGEO I), JuanLesPins, France, 2009. View at: Google Scholar
 S. S. Pandey, P. K. Kasundra, and S. D. Daxini, “Introduction of meshfree methods and implementation of Element Free Galerkin (EFG) method to beam problem,” International Journal on Theoretical and Applied Research in Mechanical Engineering 2, pp. 85–89, 2013. View at: Google Scholar
 M. Marcon, J. Vorel, K. Ninčević, and R. WanWendner, “Modeling adhesive anchors in a discrete element framework,” Materials , vol. 10, no. 8, article no. 917, 2017. View at: Publisher Site  Google Scholar
 J. Ožbolt, Y. Li, and I. Kožar, “Microplane model for concentre with relaxed kinematic constraint,” International Journal of Solids and Structures, vol. 38, no. 16, pp. 2683–2711, 2001. View at: Publisher Site  Google Scholar
 S. A. Silling, “Reformulation of elasticity theory for discontinuities and longrange forces,” Journal of the Mechanics and Physics of Solids, vol. 48, no. 1, pp. 175–209, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 S. A. Silling and E. Askari, “A meshfree method based on the peridynamic model of solid mechanics,” Computers & Structures, vol. 83, no. 1718, pp. 1526–1535, 2005. View at: Publisher Site  Google Scholar
 E. Madenci and E. Oterkus, Peridynamic Theory and Its Applications, SpringerVerlag, New York, NY, USA, 2014. View at: Publisher Site
 O. Weckner and R. Abeyaratne, “The effect of longrange forces on the dynamics of a bar,” Journal of the Mechanics and Physics of Solids, vol. 53, no. 3, pp. 705–728, 2005. View at: Publisher Site  Google Scholar  MathSciNet
 F. Bobaru, M. Yang, L. F. Alves, S. A. Silling, E. Askari, and J. Xu, “Convergence, adaptive refinement, and scaling in 1D peridynamics,” International Journal for Numerical Methods in Engineering, vol. 77, no. 6, pp. 852–877, 2009. View at: Publisher Site  Google Scholar
 B. Kilic, A. Agwai, and E. Madenci, “Peridynamic theory for progressive damage prediction in centercracked composite laminates,” Composite Structures, vol. 90, no. 2, pp. 141–151, 2009. View at: Publisher Site  Google Scholar
 Y.L. Hu, Y. Yu, and H. Wang, “Peridynamic analytical method for progressive damage in notched composite laminates,” Composite Structures, vol. 108, no. 1, pp. 801–810, 2014. View at: Publisher Site  Google Scholar
 S. A. Silling and F. Bobaru, “Peridynamic modeling of membranes and fibers,” International Journal of NonLinear Mechanics, vol. 40, no. 23, pp. 395–409, 2005. View at: Publisher Site  Google Scholar
 W. Gerstle, “Peridynamic modeling of plain and reinforced concrete structures,” in Proceedings of the 18th International Conference on Structural Mechanics in Reactor Technology, Beijing, China, 2005. View at: Google Scholar
 D. Huang, Q. Zhang, and P. Z. Qiao, “Damage and progressive failure of concrete structures using nonlocal peridynamic modeling,” Science China Technological Sciences, vol. 54, no. 3, pp. 591–596, 2011. View at: Publisher Site  Google Scholar
 F. Shen, Q. Zhang, and D. Huang, “Damage and failure process of concrete structure under uniaxial compression based on peridynamics modeling,” Mathematical Problems in Engineering, vol. 2013, Article ID 631074, 5 pages, 2013. View at: Publisher Site  Google Scholar
 D. Huang, G. Lu, C. Wang, and P. Qiao, “An extended peridynamic approach for deformation and fracture analysis,” Engineering Fracture Mechanics, vol. 141, pp. 196–211, 2015. View at: Publisher Site  Google Scholar
 D. Huang, G. Lu, and P. Qiao, “An improved peridynamic approach for quasistatic elastic deformation and brittle fracture analysis,” International Journal of Mechanical Sciences, vol. 9495, article no. 2937, pp. 111–122, 2015. View at: Publisher Site  Google Scholar
 Z. Cheng, G. Zhang, Y. Wang, and F. Bobaru, “A peridynamic model for dynamic fracture in functionally graded materials,” Composite Structures, vol. 133, pp. 529–546, 2015. View at: Publisher Site  Google Scholar
 F. Han, G. Lubineau, Y. Azdoud, and A. Askari, “A morphing approach to couple statebased peridynamics with classical continuum mechanics,” Computer Methods Applied Mechanics and Engineering, vol. 301, pp. 336–358, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 W. Gerstle, N. Sau, and E. Aguilera, “Micropolar peridynamic constitutive model for concrete,” in Proceedings of the 6th International Conference on Fracture Mechanics of Concrete and Concrete Structures (SMiRT 19 '07), Toronto, Canada, June 2007. View at: Google Scholar
 M. L. Parks, R. B. Lehoucq, S. J. Plimpton, and S. A. Silling, “Implementing peridynamics within a molecular dynamics code,” Computer Physics Communications, vol. 179, no. 11, pp. 777–783, 2008. View at: Publisher Site  Google Scholar
 B. Kilic and E. Madenci, “An adaptive dynamic relaxation method for quasistatic simulations using the peridynamic theory,” Theoretical and Applied Fracture Mechanics, vol. 53, no. 3, pp. 194–204, 2010. View at: Publisher Site  Google Scholar
 P. Underwood, “Dynamic relaxation,” in Computational Methods for Transient Analysis, T. Belytschko and T. J. R. Hughes, Eds., pp. 245–265, Elsevier Science Publishers, Amsterdam, 1983. View at: Google Scholar
Copyright
Copyright © 2018 Jiezhi Lu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.