Processing math: 100%
Research article Special Issues

Two-dimensional implementation of the coarsening method for linear peridynamics

  • Peridynamic theory was introduced to overcome the limitations of classical continuum mechanics (CCM) in handling discontinuous material response. However, for certain problems, it is computationally expensive with respect to CCM based approaches. To reduce the computational time, a coarsening method was developed and its capabilities were demonstrated for one-dimensional structures by substituting a detailed model with a surrogate model with fewer degrees of freedom. The objective of this study is to extend the application of coarsening method for linear peridynamics for two-dimensional analysis. Moreover, the existing one-dimensional coarsening method was further explored by considering various different micromodulus functions. The numerical results demonstrated that coarsening approach has a potential to reduce the computational time with high accuracy for both one-dimensional and two-dimensional problems.

    Citation: Yakubu Galadima, Erkan Oterkus, Selda Oterkus. Two-dimensional implementation of the coarsening method for linear peridynamics[J]. AIMS Materials Science, 2019, 6(2): 252-275. doi: 10.3934/matersci.2019.2.252

    Related Papers:

    [1] Holm Altenbach, Oleksiy Larin, Konstantin Naumenko, Olha Sukhanova, Mathias Würkner . Elastic plate under low velocity impact: Classical continuum mechanics vs peridynamics analysis. AIMS Materials Science, 2022, 9(5): 702-718. doi: 10.3934/matersci.2022043
    [2] Kai Friebertshäuser, Christian Wieners, Kerstin Weinberg . Dynamic fracture with continuum-kinematics-based peridynamics. AIMS Materials Science, 2022, 9(6): 791-807. doi: 10.3934/matersci.2022049
    [3] Erkan Oterkus, Selda Oterkus . Recent advances in peridynamic theory: A review. AIMS Materials Science, 2024, 11(3): 515-546. doi: 10.3934/matersci.2024026
    [4] Aylin Ahadi, Jakob Krochmal . Anisotropic peridynamic model—Formulation and implementation. AIMS Materials Science, 2018, 5(4): 742-755. doi: 10.3934/matersci.2018.4.742
    [5] Zoubir Chaieb, Ould Mohamed Ouarda, Azzeddine Abderrahmane Raho, Mouhyddine Kadi-Hanifi . Effect of Fe and Si impurities on the precipitation kinetics of the GPB zones in the Al-3wt%Cu-1wt%Mg alloy. AIMS Materials Science, 2016, 3(4): 1443-1455. doi: 10.3934/matersci.2016.4.1443
    [6] Shaohua Chen, Yaopengxiao Xu, Yang Jiao . Modeling solid-state sintering with externally applied pressure: a geometric force approach. AIMS Materials Science, 2017, 4(1): 75-88. doi: 10.3934/matersci.2017.1.75
    [7] K.A. Lazopoulos, E. Sideridis, A.K. Lazopoulos . On Λ-Fractional peridynamic mechanics. AIMS Materials Science, 2022, 9(5): 684-701. doi: 10.3934/matersci.2022042
    [8] Hugo F. Lopez . Microstructural features associated with the effect of temperature on the dimensional stability of an automotive Al-A319 alloy. AIMS Materials Science, 2016, 3(2): 634-644. doi: 10.3934/matersci.2016.2.634
    [9] Erkan Oterkus, Timon Rabczuk, Selda Oterkus . Special Issue: Peridynamics and its applications. AIMS Materials Science, 2024, 11(3): 602-604. doi: 10.3934/matersci.2024030
    [10] Bozo Vazic, Hanlin Wang, Cagan Diyaroglu, Selda Oterkus, Erkan Oterkus . Dynamic propagation of a macrocrack interacting with parallel small cracks. AIMS Materials Science, 2017, 4(1): 118-136. doi: 10.3934/matersci.2017.1.118
  • Peridynamic theory was introduced to overcome the limitations of classical continuum mechanics (CCM) in handling discontinuous material response. However, for certain problems, it is computationally expensive with respect to CCM based approaches. To reduce the computational time, a coarsening method was developed and its capabilities were demonstrated for one-dimensional structures by substituting a detailed model with a surrogate model with fewer degrees of freedom. The objective of this study is to extend the application of coarsening method for linear peridynamics for two-dimensional analysis. Moreover, the existing one-dimensional coarsening method was further explored by considering various different micromodulus functions. The numerical results demonstrated that coarsening approach has a potential to reduce the computational time with high accuracy for both one-dimensional and two-dimensional problems.


    Peridynamics is a nonlocal continuum theory developed by Silling [1] with the objective of overcoming the shortcomings associated with the Classical Continuum theory particularly as it relates to the ability to model fracture and damage in materials. To achieve this objective, the equation of motion of the classical theory was cast in an integral form. This paved the way for the unification of the mathematical modelling of continuous media, fracture and particles in a single modelling framework. Thus, peridynamics provided a robust structure that can find application in situations involving evolution and propagation of discontinuities such as crack nucleation and growth using the same field equation as in the continuous case. Another feature of the Peridynamic formulation is that it is a nonlocal theory that incorporates the concept of long-range forces by allowing interaction of particles located at a finite distance from each other through a pairwise force field.

    The peridynamic theory has extensively been investigated and is shown to have capabilities of modelling a wide range of engineering problems [2,3,4,5,6,7,8,9,10]. However, for certain problems, it suffers from the disadvantage of being computationally more expensive compared to the classical continuum theory which is attributable to its nonlocal formulation. In an effort to mitigate against this issue, researchers have proposed a number of multiscale modelling techniques. An adaptive refinement method was proposed in [11] for bond-based peridynamic problems involving a variable horizon. A peridynamic-based framework for hierarchical multiscale modelling was proposed in [12] that couples molecular dynamics and peridynamics. These techniques allow coupling of models at different scales for effective inter-scale information transfer. A coarsening method was also proposed in [13]. The key idea of the coarsening method is to obtain accurate results by solving a reduced order model. The reduction in the order of model is achieved by substituting a detailed model with a surrogate model with fewer degrees of freedom. The objective of this study is to extend the application of a coarsening method for linear peridynamics proposed in [13] for two-dimensional analysis. Moreover, the existing one-dimensional coarsening method was further explored by considering different micromodulus functions other than that considered in [14].

    Recalling from the classical continuum theory, the equation of motion of a medium arising from conservation of momentum is as follows:

    ρ(x)¨u(x,t)=σ +b(x,t) (1)

    where ρ is the mass density of the medium, σ is the stress tensor, b is a vector of body force density, is the acceleration vector field of the material point x within the medium at time t. The key challenge in using Eq 1 to model discontinuous system behaviour is the presence of the divergence operator, which implies the existence of the spatial derivative of the stress field and consequently the displacement field within the domain of interest. However, since these field variables are not continuous over features such as crack tip and crack surfaces, the derivatives in such instance are undefined. In the peridynamic formulation, the equation of motion was casted such that the integral operators replaced the derivatives on the right hand side of Eq 1. A "bond-based" peridynamic equation of motion was originally proposed by Silling [1] as

    ρ(x)¨u(x,t)=Hxf(u(q,t)u(x,t),qx)dVq+b(x,t) (2)

    where u is the displacement vector field, Hx is the neighbourhood of the particle located at point x. f is the pairwise force function with a unit of force per unit volume squared that a particle located at point q exerts on a particle located at point x.

    If we assume a linear material behaviour, then the pairwise force function [1] takes the form

    f(η,ξ)=C(ξ)ηξ,η (3)

    where C(ξ) is a tensor valued function called the micromodulus function given by

    C(ξ)=f(0,ξ)η            ξ (4)

    ξ = q-x is the relative position of the particles in the reference configuration, and η = u(q, t)-u(x, t) is their relative displacement. The peridynamic equation of motion (Eq 2) therefore specialises to

    ρ(x)¨u(x,t)=HxC(x,q)(u(q,t)u(x,t))dVq+b(x,t) (5)

    The micromodulus function C(ξ) is required to satisfy certain properties [1]. Moreover, the range of admissible functions that can be used is extensive [14]. However, only a few of them will be discussed here.

    The first sets of examples of the micromodulus functions to be given belong to a class of materials called microelastic materials. A peridynamic material is said to be microelastic if the pairwise force function, f, is derived from a scalar micropotential, w:

    f(η,ξ)=w(η,ξ)η (6)

    The micropotential, w is the energy of a single bond and represents a local strain energy density [2] which is given by

    w(η,ξ)=c(|ξ|)s2|ξ|2 (7)

    where the stretch, s can be defined as

    s=|ξ+η||ξ||ξ| (8)

    From Eqs 7 and 8, Eq 6 becomes

    f(η,ξ)=c(|ξ|)s(ξ+η)|ξ+η|         ξ,η (9)

    where c(|ξ|) is a function that is characteristic of the material and is called the bond constant. Differentiating Eq 9 with respect to η and assuming small displacements yields

    C(ξ)= c(|ξ|)|ξ|3[ξ2xξxξyξxξyξ2y]           ξ (10)

    for two-dimensional configurations. For one-dimensional configurations, Eq 10 specialises to

    C(ξ)= c(|ξ|)|ξ|           ξ (11)

    In order to determine the bond constants for microelastic materials, the strain energy density of a given point can be defined [2] as

    W=12HwdVξ (12)

    Let δ be the radius of the influence domain or horizon that is characteristic of the material. Then, using Eq 7 and assuming a one-dimensional problem, the peridynamic strain energy density of a material point becomes

    W=A4δδc(|ξ|)s2|ξ|dξ (13)

    If we assume c(|ξ|) to be independent of bond length, i.e., c(|ξ|)=c0, then Eq 13 can be written for constant stretch for all bonds as

    W=A4δδc0s2|ξ|dξ=Ac0s2δ24 (14)

    For the same condition, the strain energy density from classical continuum theory can be written as

    W=Es22 (15)

    Equating Eqs 14 and 15 yields

    c0=2EAδ2 (16)

    Hence, the micromodulus function for this material is

    C(ξ)= {2EAδ2|ξ|,if   |ξ|δ0,     if   |ξ|>δ (17)

    In the case where c(|ξ|) is dependent on the bond length and assumed to have a triangular profile, i.e., c(|ξ|)=c1(1|ξ|δ), then Eq 13 becomes

    W=14δδc1(1|ξ|δ)s2|ξ|dξ=Ac1s2δ212 (18)

    Equating Eq 18 with Eq 15 gives

    c1=6EAδ2 (19)

    The micromodulus function in this case is

    C(ξ)= {6EAδ2|ξ|(1|ξ|δ),if   |ξ|δ0,              if   |ξ|>δ (20)

    In the case where the dependency of c(|ξ|) is in the form of an inverted triangular profile [3], in which case, bonds closer to the principal point are softer than bonds further away from it, such that, c(|ξ|)=c2(|ξ|/δ) for bonds within the horizon. In this case, Eq 13 becomes

    W=A4δδc2(|ξ|δ)s2|ξ|dξ=Ac2s2δ26 (21)

    Equating Eq 21 with Eq 15 gives

    c2=3EAδ2 (22)

    The micromodulus function for this scenario then becomes

    C(ξ)= {3EAδ3,if   |ξ|δ0,  if   |ξ|>δ (23)

    Notice that the micromodulus function in this case is a constant function. Figure 1a shows a plot of the three different bond constants as functions of α while Figure 1b is a plot of the corresponding micromodulus functions as functions of α, where α=E/6A.

    Figure 1.  (a) Bond elastic functions, (b) Micromodulus functions as ξ changes.

    In the two-dimensional case, let δ be the radius of sphere H describing the horizon. Then, the two-dimensional equivalence of the micromodulus Eqs 17, 20 and 23 are respectively given as

    C(ξ)= {9Eπδ3h|ξ|3[ξ2xξxξyξxξyξ2y],if   |ξ|δ0,                   if   |ξ|>δ (24)
    C(ξ)= {36Eπδ3h|ξ|3(1|ξ|δ)[ξ2xξxξyξxξyξ2y],if   |ξ|δ0,                           if   |ξ|>δ (25)
    C(ξ)= {8Eπδ4h|ξ|2[ξ2xξxξyξxξyξ2y],if   |ξ|δ0,                   if   |ξ|>δ (26)

    Notice that Eq 26 unlike its one-dimensional counterpart is not a constant function.

    The goal of the coarsening method is to provide a multiscale framework based on peridynamic theory that will predict macroscale response of a medium based on the underlying evolving microstructure, without necessarily resolving all microstructural details. This is achieved by a process of successive elimination of points from the medium called coarsening. Each successive coarsening results in a medium with a reduced material points as well as reduced level of geometrical detail. However, the properties of the coarsened medium are determined such that the effects of the excluded material points are implicitly included in the coarsened level simulation. The coarsening formulation presented here is adapted from [14] for the completeness of the paper.

    To coarsen a detailed model as proposed in [13], let B0 be a linear elastic peridynamic body as shown in Figure 2. Let A0 be the set of linear admissible displacement fields on B0, and let C0:B0×B0 be the micromodulus tensor associated with the material of the body where L is the set of all second order tensors. Assume r0 is a positive number that represents the maximum interaction distance for all points in B0, such that if

    |qx|>r0C0(x,q)=0x,qB0 (27)
    Figure 2.  A Peridynamic body showing model levels 0, 1 and 2.

    Let B1B0, and let A1 represent the set of admissible displacement fields on B1. B0 and B1 are called level 0 and level 1 body, respectively. The objective of this coarsening process is to express the internal forces acting on B1 in terms of its own displacements only while implicitly accounting for forces points in B0-B1 exert on B1. To achieve this objective, let x be an arbitrary point in B1 and let δ1 be a positive number. Let H1x be the closed neighbourhood of x in B0 with radius r1

    H1x={qB0||qx|r1} (28)

    Let R1x = H1xB1. Suppose u1A1 is given and let u0A0 satisfy the compatibility condition

    u0(p)=u1(p)pR1x (29)

    Outside of R1x, assume that u0 satisfies the equilibrium condition, neglecting interactions between H1x and its exterior:

    L0(z)+b(z)=0              zH1xR1x (30)

    where

    L0(z)=H1xC0(z,p)(u0(p)u0(z))dVpzH1x (31)

    Also, assume that there is no body force density applied outside R1x

    b(z)=0              zH1xR1x (32)

    Further assume that for a given displacement field u1 Eqs 29 and 30 have a unique solution u0 on H1x, and let Sx0, 1 be the resolvent kernel that generates this solution:

    u0(p)=R1xS0,1x(p,q)(u1(q))dVqpH1x (33)

    From Eqs 29 and 33, we can infer that

    S0,1x(p,q)=IΔ(pq)      pR1x,qH1x (34)

    where I is the isotropic tensor and ∆ is the three dimensional Dirac delta function. For the special case where u0 represents a linear rigid translation of points in R1x through an arbitrary displacement u, then all points in H1x-R1x will also translate by u. Therefore, from Eq 33, we see that

    ˉu=[R1xS0,1x(p,q)dVq]ˉupH1x (35)

    from which we can obtain the following identity

    R1xS0,1x(p,q)dVq=1pH1x (36)

    Subtracting u0(z) from both sides of Eq 33 and using the identity given in Eq 36, we can obtain

    u0(p)u0(z)=R1xS0,1x(p,q)(u1(q)u1(z))dVqp,zH1x (37)

    Substituting this result in Eq 31 yields

    L0(z)=H1xC0(z,p)[R1xS0,1x(p,q)(u1(q)u1(z))dVq]dVpzH1x (38)

    Reversing the order of integration and rearranging gives

    L0(z)=R1x[H1xC0(z,p)S0,1x(p,q)dVp](u1(q)u1(z))dVqzH1x (39)

    Recalling that x is an arbitrary point in B1, if we denote the force density in any such choice of x by

    L0(x)=L1(x)      xB1 (40)

    Then, from Eqs 39 and 40 we can obtain

    L1(x)=R1xC1(x,q)(u1(q)u1(x))dVqxB1 (41)

    where the coarsened micromodulus in level 1 C1: B1×B1 is defined by

    C1(z,q)=H1xC0(x,p)S0,1x(p,q)dVpx,qB1 (42)

    Similarly, for any level m in the coarsening process, the force density can be obtained from

    Lm(x)=RmxCm(x,q)(um(q)um(x))dVqxBm (43)

    where the coarsened micromodulus in level m, Cm: Bm×Bm is defined by

    Cm(z,q)=HmxCm1(x,p)Sm1,mx(p,q)dVpx,qBm (44)

    To carry out the numerical implementation of the coarsening method described in the preceding section, B0 is discretised into nodes, which for simplicity are taken to have equal volume v. Let xi be the position of node i in level mi. For any nodes i and j, let

    C0i,j=vC0(xi,xj) (45)

    The coarsened micromodulus C1i,j is obtained by discretizing Eq 42 as

    C1i,j=vNk=1C0i,kS0,1k,j (46)

    To evaluate Eq 46, the resolvent kernel S0,1k,j needs to be determined. A convenient means of obtaining the kernel function was proposed in [13]. Successive coarsening to higher levels can be achieved by following the same procedure as above, so that for any m≥1,

    Cmi,j=vNmk=1Cm1i,kSm1,mk,j (47)

    To compute the coarsened micromodulus at level m for node i, the resolvent kernel Sm1,mk,j must be determined. A method was suggested in [13] for determining the resolvent kernel. This method proceeds by creating and merging two sub-systems of equations into a global system of equations. To create such a system, two vectors {u0} representing the displacement field in the level 0 body and {u1} representing the displacement field in level 1 body are defined. The relationship between these two displacement fields is formulated as follows:

    The displacement fields ui0 and ui1 for nodes in level 1 body are constrained to satisfy Eq 29:

    u0i=u1iiB1 (48)

    The second system of equations is obtained by ensuring that the nodes complementary to nodes in level 1 body H1x-R1x satisfy Eqs 30, 31 and 32. Concatenation of these two systems of equations results in a global system of equations that has the following form in one dimension:

    [10000100C0i,i1PiC0i,i+1C0N,N1PN]{u01u0Ru0iu0N}={u11u1R00} (49)

    where R is the number of nodes in R1x, N is the number of nodes in H1x and the diagonal elements Pi are given by

    Pi=jiC0i,j (50)

    In shorthand notation, Eq 49 can be written as

    [A]{u0}={b} (51)

    where [A] is N×N square matrix in one-dimension and 2N×2N in two-dimensions. Let [A]-1 be the inverse of [A], [S0, 1] in one-dimension will be the leftmost R columns of [A]-1. In two-dimensions, [S0, 1] is the leftmost 2R columns of [A]-1.

    In this section, the coarsened form of various micromodulus functions will be presented. One-dimensional cases will be followed by a two dimensional case.

    The one-dimensional problem is a homogeneous bar with a length of 1 m. The elastic modulus of the material is assumed to be 200 GPa. Three scenarios will be investigated by considering the following micromodulus functions:

    i.          C(ξ)= {2EAδ2|ξ|,if   |ξ|δ0,    if   |ξ|>δ (52)
    ii.          C(ξ)= {6EAδ2|ξ|(1|ξ|δ),if   |ξ|δ0,           if   |ξ|>δ (53)
    iii.         C(ξ)= {3EAδ3,if   |ξ|δ0,  if   |ξ|>δ (54)

    The level 0 horizon is specified as δ0 = 25 mm. Note that this same horizon will be used in the coarsened models. The bar is discretized into 400 nodes with spacing ∆x = 2.5 mm. Coarsening of the detailed model to the level 1 body is carried out by retaining every fourth node of the level 0 body. Similarly, coarsening of the level 1 body to the level 2 body proceeds by retaining every second node in the level 1 body. Figure 3 shows a schematic representation of the coarsening process. The micromodulus of the detailed model (level 0) C0 as well as the coarsened micromodulus C1 and C2 for the three scenarios are shown in Figures 4 to 6 These curves are characterized by sharp peaks consistent with the fact that the coarsened micromodulus functions are defined only at their respective coarsened regions.

    Figure 3.  Schematic representation of the one-dimensional coarsening process.
    Figure 4.  Micromodulus function associated with a uniform bond constant.
    Figure 5.  Micromodulus function associated with a bond constant function having triangular profile.
    Figure 6.  Micromodulus function associated with a bond constant function having inverted triangular profile.

    The two-dimensional problem is a 500×500 mm square plate with a thickness of 50 mm. Let the elastic modulus of the plate to be 200 GPa and its micromodulus function is given by

    C(ξ)= {36Eπδ3h|ξ|3(1|ξ|δ)[ξ2xξxξyξxξyξ2y],if   |qx|δ0,                           if   |qx|>δ (55)

    Coarsening is carried out as shown schematically in Figure 7. Every second node in every second row of the level 0 model is retained in the coarsened level 1 model. Similarly, every second node in the second row of level 1 is retained in level 2 body. Figure 8 shows the micromodulus function C0 of the level 0 model as well as the coarsened micromodulus C1 of the level 1 model. It is noticed that the coarsened micromodulus C1 is characterised by sharp peaks which is a reflection of the fact that it is only defined at the coarsened region. Figure 8a shows the response of the i-th material point in the x-direction for a displacement of the j-th point in the x-direction. Figure 8b shows the response of the i-th material point in the x-direction if the j-th point is displaced in the y-direction. Figure 8c shows the response in the y-direction of the i-th point when the j-th point displaces in the x-direction. Finally, Figure 8d shows the response of the i-th point in the y-direction when the j-th point is displaced in the y-direction.

    Figure 7.  Detailed level 0 and coarsened level 1 for coarsening of two-dimensional micromodulus function.
    Figure 8.  Coarsening of two-dimensional micromodulus function: (a) x-x interaction (b) x-y interaction (c) y-x interaction, and (d) y-y interaction.

    This section comprises five numerical experiments to illustrate the application of the coarsening method in solving solid mechanics problems. The systems are assumed to be in a state of static equilibrium and thus an implicit static solution scheme has been utilized. The first three examples focus a one-dimensional bar under tension loading cases for homogeneous and composite materials with and/or without defect. The last two examples demonstrate the capability of the coarsening approach for two dimensional plate problems for isotropic and composite materials.

    This example is a numerical experiment aimed at illustrating the capability of the coarsening method in capturing the correct bulk behaviour of a system even with reduced degrees of freedom as opposed to if the degrees of freedom were to be reduced without the coarsening process. Consider a bar of length 1 m of material having a micromodulus function of the form:

    C(ξ)={6EAδ2|ξ|(1|ξ|δ),if   |qx|δ0,             if   |qx|>δ (56)

    The maximum interaction distance δ = 10∆x, where ∆x is the length of a material point. The detailed model (level 0) is discretised into 1000 points, so that ∆x = 1 mm. The elastic modulus of the material is E = 200 GPa. Coarsening of the level 0 model to level 1 proceeds by retaining every fourth node in level 0. Level 2 model retains every second node in level 1 model as shown in Figure 3. A body force density of b = 200 N/mm3 is applied to the rightmost level 2 node, while the third leftmost level 2 node is constrained to have zero displacement. The displacement field associated with the level 0 and coarsened levels 1 and 2 are shown in Figure 9. It is observed that even with reduced degrees of freedom as the detailed model is coarsened into level 1 and subsequently level 2, the displacement fields obtained from all these models are very similar to each other. This would not have been the case if the reduction in degrees of freedom were not achieved through the coarsening process.

    Figure 9.  Displacement fields for a one-dimensional homogeneous bar under tension loading.

    The analytical solution based on the classical theory to the axial bar problem is of the form

    u(x)=(FEA)x (57)

    where x is the distance from the support, u(x) is the displacement at point x, E is the elastic modulus, A is the cross sectional area, and F is the force applied at the free end of the bar. The result from the analytical model and those from the detail level 0 model as well as the results from coarsened levels 1 and 2 model are almost identical.

    To illustrate the need for the coarsening process in achieving model order reduction, the number of nodes was reduced from 1000 to 250 and further to 125 without following the coarsening procedure. As shown in Figure 9, although the results from the coarsened models and the results from the models with reduced number of nodes are very close to each other, coarsened model results agree better with the analytical solution.

    In this example, a one-dimensional composite bar with a length of 1m and a cross-sectional area of 100 mm2 is considered. The composite bar has a periodic microstructure consisting of alternating strips Shard and Ssoft of the hard and the soft materials, respectively. Each strip is 50 mm in length. The horizon for both hard and soft section is 50 mm. The micromodulus of the composite system has the form:

    C(ξ)={6E1Aδ2|ξ|(1|ξ|δ),if|qx|δand(xShardandqShard)2E2Aδ2|ξ|,if|qx|δand(xSsoftorqSsoft)0,otherwise (58)

    In other words, bonds with both ends in hard strip have hard properties while bonds with either ends in soft strip have soft properties. The elastic moduli of the hard and soft materials are E1 = 200 GPa and E2 = 9 GPa, respectively. The level 0 horizon is specified as δ0 = 10∆x. Coarsening from the level 0 to the level 1 body is carried out as shown in Figure 3 by retaining every fourth node of the level 0 body. Similarly, the level 1 body is coarsened to level 2 body by retaining every second node in the level 1 body.

    A body force density of b = 0.001 N/mm3 is applied to the three rightmost level 2 nodes, while the three leftmost level 2 nodes are constrained to have zero displacement. The computed displacements of coarsened level 0, 1 and 2 are shown in Figure 10. As expected, computation from the level 0 revealed the most detailed microstructural information. As the computational regions get coarsened, these microstructural details get smoothed out. However, both detailed and coarsened models have similar global stretch. This demonstrates the fact that the effective properties produced by the coarsening process accurately reflect the bulk properties of the composite.

    Figure 10.  Displacement fields for levels 0, 1, and 2 for a one-dimensional composite bar under tension loading.

    This example considers a homogeneous bar with a defect. The bar has a length of 1 m. The constitutive model of the bar is given by

    C(ξ)=φ(x,q){2EAδ2|ξ|,if|qx|δ0,if|qx|>δ (59)

    where φ is a degrading term given by

    φ={0.1,ifxxdqorqxdx1.0otherwise (60)

    The horizon is specified as δ0 = 50 mm in the level 0 and the bar is discretised into 200 nodes. The defect is located at the center of the bar as shown in Figure 11. The level 1 contains every third node in the level 0 and similarly, the level 2 contains every third node of the level 1 as shown in Figure 12. Prescribed displacement boundary conditions are applied to the three leftmost and three rightmost level 2 nodes. The values of the prescribed displacements are given by u0i=u1i=u2i=0.01xi, where xi is the coordinate of the node.

    Figure 11.  One-dimensional homogeneous bar with a defect.
    Figure 12.  Schematic representation of the coarsening process for example 2.

    The displacement fields for the level 0 model as well as the coarsened level 1 and level 2 models near the defect are shown in Figure 13. The three levels give identical displacement profile except for close to the site of the defect. This reflects the fact that due to fewer nodes in the coarsened models, wider spacing exits between nodes.

    Figure 13.  Coarsened displacement fields in a one-dimensional homogeneous bar with a defect.

    This example concerns a two-dimensional homogeneous plate under tension loading. The plate is a 500 mm×500 mm square plate. The micromodulus of the plate is chosen as

    C(ξ)= {8Eπδ4h|ξ|2[ξ2xξxξyξxξyξ2y],if   |qx|δ0,                   if   |qx|>δ (61)

    The horizon for the level 0 body is specied as δ0 = 3.015∆x. The plate is discretised into 200 points in both x and y-directions with every point having ∆x = ∆y = 2.5 mm. The thickness of the plate is ∆x. Coarsening of the detailed model is achieved by retaining every second point in the second row as shown schematically in Figure 14.

    Figure 14.  Discretisation and coarsening of the two-dimensional homogeneous plate.

    Two boundary layers each consisting of one grid of 100 points are created; one along the top edge and the other along the bottom edge. At the top boundary layer, every second point of the detailed model is constrained to have zero displacement. A body force density of 2000 N/mm3 is applied to every second point in the bottom boundary layer in the detailed model.

    The displacement field obtained from solving both a detailed (level 0) and the coarsened model (level 1) are similar as can be seen from Figure 15 and Figure 16. The displacement profile of a specific grid of points along the y-axis and x-axis are plotted and shown in Figure 17 and Figure 18 respectively. It is observed that both models result in very similar displacement fields.

    Figure 15.  Profile of x-axis displacement field: (a) Detailed (Level 0) model and (b) Coarsened (Level 1) model.
    Figure 16.  Profile of y-axis displacement field: (a) Detailed (Level 0) model and (b) Coarsened (Level 1) model.
    Figure 17.  Profile of displacement field along a grid line of points in the y-direction.
    Figure 18.  Profile of displacement field along a grid line of points in the x-direction.

    This example will consider a two-dimensional 500 mm×500 mm square plate with periodic microstructure consisting of alternative strips Shard and Ssoft of hard and soft materials, respectively as shown in Figure 19. Each strip is 50 mm in length. The horizon for both hard and soft sections is 50 mm. The micromodulus of the composite system is given by

    C(ξ)={36Eπδ3h|ξ|3(1|ξ|δ)[ξ2xξxξyξxξyξ2y],if  |qx|δ and (xShard and qShard)9Eπδ3h|ξ|3[ξ2xξxξyξxξyξ2y],if   |qx|δ   and (xSsoft or qSsoft)0,otherwise    (62)
    Figure 19.  Discretisation and coarsening of the two-dimensional composite plate.

    The horizon is specified as δ0 = 3.015∆x. The plate is discretized into 200 points in both x and y-directions with every point having ∆x = ∆y = 2.5 mm. Coarsening of the detailed model and the applied boundary condition is same as that of previous example.

    The computed displacement profile along the y-direction for the identical boundary condition for the detailed model (level 0) and the coarsened model (level 1) is shown in Figure 20. It is observed that because the detailed model included more material points, it has been able to resolve more microstructural details than the coarsened models. However, both models have the same global stretch thereby demonstrating that the effective properties obtained using this coarsening method is an accurate representation of the bulk properties of the composite.

    Figure 20.  Displacement fields for levels 0, 1, and 2 for a two-dimensional plate under tension loadin.

    Determination of computational cost savings arising from coarsening a numerical model will follow from the idea suggested in [1]. Consider a body that is discretized into N0 nodes in the level 0 model. Let J=aNpM be the number of arithmetic operations used by a linear solver to solve the fully coarsened level M model, where a and p are constants (n = 3 for Gaussian elimination [1]), and NM is the total number of nodes in level M model. Suppose each level has 1/α as many nodes as the previous level, where α is a constant, such that

    NM=N0αDM (63)

    where D is the dimension of the problem (In the example on two dimensional homogeneous plate, α = 2 and D = 2).

    The computational effort in determining material properties at level m+1 from level m is proportional to Nm since in doing so, the inverse of the matrix ANm×Nm has to be computed. Therefore, for some positive constant b, the total arithmetic operations that will be used to coarsen a numerical model and solve the fully coarsened level M model can be written as

    J=a(N0αDM)n+M1m=0bN0αDMJ=a(N0αDM)n+bN011αDM11αJ<a(N0αDM)n+bN0αα1 (64)

    It can be concluded from Eq 64 that the computational effort of coarsening up to level M in a linear solver for a boundary value problem is reduced by a factor of αnDM over the computational effort needed to solve the detailed level 0 model. For a two-dimensional problem such as the example on two-dimensional homogeneous plate, this factor will be 23×2×M = 26M. The computational price paid to determine the coarsened properties of the model is less than bN0α/(α1) which is independent of M.

    The coarsening method developed in [13] extended for two-dimensional application in this work has demonstrated its robustness in modelling a range of problems. Moreover, several one-dimensional problems were considered to investigate the effect of different micromodulus functions other than considered in [14]. The examples solved in this work have demonstrated that this method is able to relate the bulk behaviour of materials to the details of their microstructure in two-dimensional problems. This allows us to reduce the order of the problem without losing accuracy in predicting the bulk response of the system. This translates to a reduction in computational time and storage required. Another advantage offered by the coarsening method is that the model retains its attributes such as boundary conditions and applied force density after the coarsening process. This is not the case if the model is reduced without going through the coarsening process. In this case, the boundary conditions and applied force density can no longer be identically applied as in the original model.

    It is instructive to mention a number of caveats related to this method. This coarsening method is currently valid only for linear and static problems. This limitation is a consequence of assumption made in obtaining Eq 30.

    The first author is supported by the government of the federal republic of Nigeria through the Petroleum Technology Development Fund (PTDF).

    The authors declare that there is no conflict of interest regarding the publication of this manuscript.



    [1] Silling SA (2000) Reformulation of elasticity theory for discontinuities and long-range forces. J Mech Phys Solids 48: 175–209. doi: 10.1016/S0022-5096(99)00029-0
    [2] Silling SA, Askari E (2005) A meshfree method based on the peridynamic model of solid mechanics. Comput Struct 83: 1526–1535. doi: 10.1016/j.compstruc.2004.11.026
    [3] Silling SA, Zimmermann M, Abeyaratne R (2003) Deformation of a Peridynamic Bar. J Elasticity 73: 173–190. doi: 10.1023/B:ELAS.0000029931.03844.4f
    [4] Ha YD, Bobaru F (2011) Characteristics of dynamic brittle fracture captured with peridynamics. Eng Fract Mech 78: 1156–1168. doi: 10.1016/j.engfracmech.2010.11.020
    [5] Parks ML, Lehoucq RB, Plimpton SJ, et al. (2008) Implementing peridynamics within a molecular dynamics code. Comput Phys Commun 179: 777–783. doi: 10.1016/j.cpc.2008.06.011
    [6] Chen X, Gunzburger M (2011) Continuous and discontinuous finite element methods for a peridynamics model of mechanics. Comput Method Appl M 200: 1237–1250. doi: 10.1016/j.cma.2010.10.014
    [7] Ghajari M, Iannucci L, Curtis P (2014) A peridynamic material model for the analysis of dynamic crack propagation in orthotropic media. Comput Method Appl M 276: 431–452.
    [8] Huang D, Lu GD, Wang CW, et al. (2015) An extended peridynamic approach for deformation and fracture analysis. Eng Fract Mech 141: 196–211. doi: 10.1016/j.engfracmech.2015.04.036
    [9] Ha YD, Bobaru F (2010) Studies of dynamic crack propagation and crack branching with peridynamics. Int J Fracture 162: 229–244. doi: 10.1007/s10704-010-9442-4
    [10] Agwai A, Guven I, Madenci E (2011) Predicting crack propagation with peridynamics: a comparative study. Int J Fracture 171: 65. doi: 10.1007/s10704-011-9628-4
    [11] Bobaru F, Ha YD (2011) Adaptive refinement and multiscale modeling in 2D peridynamics. J Multiscale Com 9: 635–659. doi: 10.1615/IntJMultCompEng.2011002793
    [12] Rahman R, Foster JT, Haque A (2014) A multiscale modeling scheme based on peridynamic theory. Int J Multiscale Com 12: 223–248. doi: 10.1615/IntJMultCompEng.2014007954
    [13] Silling SA (2011) A coarsening method for linear peridynamics. Int J Multiscale Com 9: 609–622. doi: 10.1615/IntJMultCompEng.2011002674
    [14] Bobaru F, Yang M, Alves LF, et al. (2009) Convergence, adaptive refinement, and scaling in 1D peridynamics. Int J Numer Meth Eng 77: 852–877. doi: 10.1002/nme.2439
  • This article has been cited by:

    1. Cong Tien Nguyen, Selda Oterkus, Erkan Oterkus, An energy-based peridynamic model for fatigue cracking, 2021, 241, 00137944, 107373, 10.1016/j.engfracmech.2020.107373
    2. N. Zhu, E. Oterkus, Calculation of Stress Intensity Factor using Displacement Extrapolation Method in Peridynamic Framework, 2020, 36, 1727-7191, 235, 10.1017/jmech.2019.62
    3. Yakubu Kasimu Galadima, Erkan Oterkus, Selda Oterkus, Model order reduction of linear peridynamic systems using static condensation, 2020, 1081-2865, 108128652093704, 10.1177/1081286520937045
    4. Yakubu Kasimu Galadima, Wenxuan Xia, Erkan Oterkus, Selda Oterkus, 2021, 9780128200698, 371, 10.1016/B978-0-12-820069-8.00018-4
    5. Yakubu Kasimu Galadima, Wenxuan Xia, Erkan Oterkus, Selda Oterkus, A computational homogenization framework for non-ordinary state-based peridynamics, 2022, 0177-0667, 10.1007/s00366-021-01582-6
    6. Yakubu Kasimu Galadima, Erkan Oterkus, Selda Oterkus, Static condensation of peridynamic heat conduction model, 2022, 27, 1081-2865, 2689, 10.1177/10812865221081160
    7. Erkan Oterkus, Mazdak Ghajari, 2022, 9780128035818, 10.1016/B978-0-12-822944-6.00023-2
    8. Yakubu Kasimu Galadima, Wenxuan Xia, Erkan Oterkus, Selda Oterkus, Peridynamic computational homogenization theory for materials with evolving microstructure and damage, 2022, 0177-0667, 10.1007/s00366-022-01696-5
    9. Yakubu Kasimu Galadima, Selda Oterkus, Erkan Oterkus, Islam Amin, Abdel-Hameed El-Aassar, Hosam Shawky, Modelling of viscoelastic materials using non-ordinary state-based peridynamics, 2023, 0177-0667, 10.1007/s00366-023-01808-9
    10. Yakubu Kasimu Galadima, Selda Oterkus, Erkan Oterkus, Islam Amin, Abdel-Hameed El-Aassar, Hosam Shawky, Effect of phase contrast and inclusion shape on the effective response of viscoelastic composites using peridynamic computational homogenization theory, 2024, 31, 1537-6494, 155, 10.1080/15376494.2023.2218364
    11. Erkan Oterkus, Selda Oterkus, Recent advances in peridynamic theory: A review, 2024, 11, 2372-0484, 515, 10.3934/matersci.2024026
    12. Valeriy A. Buryachenko, Transformation Field Analysis and Clustering Discretization Method in Peridynamic Micromechanics of Composites, 2024, 6, 2522-896X, 626, 10.1007/s42102-023-00113-9
    13. Valeriy A. Buryachenko, Peridynamic Micromechanics of Composites: A Review, 2024, 6, 2522-896X, 531, 10.1007/s42102-024-00122-2
    14. Yakubu Kasimu Galadima, Selda Oterkus, Erkan Oterkus, Islam Amin, Abdel-Hameed El-Aassar, Hosam Shawky, A nonlocal method to compute effective properties of viscoelastic composite materials based on peridynamic computational homogenization theory, 2023, 319, 02638223, 117147, 10.1016/j.compstruct.2023.117147
  • Reader Comments
  • © 2019 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4832) PDF downloads(741) Cited by(14)

Figures and Tables

Figures(20)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog