No Arabic abstract
The overarching goal of this work is to develop an accurate, robust, and stable methodology for finite deformation modeling using strong-form peridynamics (PD) and the correspondence modeling framework. We adopt recently developed methods that make use of higher-order corrections to improve the computation of integrals in the correspondence formulation. A unified approach is presented that incorporates the reproducing kernel (RK) and generalized moving least square (GMLS) approximations in PD to obtain higher-order gradients. We show, however, that the improved quadrature rule does not suffice to handle correspondence-modeling instability issues. In Part I of this paper, a bond-associative, higher-order core formulation is developed that naturally provides stability. Numerical examples are provided to study the convergence of RK-PD, GMLS-PD, and their bond-associat
The overarching goal of this work is to develop an accurate, robust, and stable methodology for finite deformation modeling using strong-form peridynamics (PD) and the correspondence modeling framework. We adopt recently developed methods that make use of higher-order corrections to improve the computation of integrals in the correspondence formulation. A unified approach is presented that incorporates the reproducing kernel (RK) and generalized moving least square (GMLS) approximations in PD to obtain higher-order gradients. We show, however, that the improved quadrature rule does not suffice to handle correspondence-modeling instability issues. In Part I of this paper, a bond-associative, higher-order core formulation is developed that naturally provides stability. Numerical examples are provided to study the convergence of RK-PD, GMLS-PD, and their bond-associat
In this work, we study the reproducing kernel (RK) collocation method for the peridynamic Navier equation. We first apply a linear RK approximation on both displacements and dilatation, then back-substitute dilatation, and solve the peridynamic Navier equation in a pure displacement form. The RK collocation scheme converges to the nonlocal limit and also to the local limit as nonlocal interactions vanish. The stability is shown by comparing the collocation scheme with the standard Galerkin scheme using Fourier analysis. We then apply the RK collocation to the quasi-discrete peridynamic Navier equation and show its convergence to the correct local limit when the ratio between the nonlocal length scale and the discretization parameter is fixed. The analysis is carried out on a special family of rectilinear Cartesian grids for the RK collocation method with a designated kernel with finite support. We assume the Lam{e} parameters satisfy $lambda geq mu$ to avoid adding extra constraints on the nonlocal kernel. Finally, numerical experiments are conducted to validate the theoretical results.
As a nonlocal extension of continuum mechanics, peridynamics has been widely and effectively applied in different fields where discontinuities in the field variables arise from an initially continuous body. An important component of the constitutive model in peridynamics is the influence function which weights the contribution of all the interactions over a nonlocal region surrounding a point of interest. Recent work has shown that in solid mechanics the influence function has a strong relationship with the heterogeneity of a materials micro-structure. However, determining an accurate influence function analytically from a given micro-structure typically requires lengthy derivations and complex mathematical models. To avoid these complexities, the goal of this paper is to develop a data-driven regression algorithm to find the optimal bond-based peridynamic model to describe the macro-scale deformation of linear elastic medium with periodic heterogeneity. We generate macro-scale deformation training data by averaging over periodic micro-structure unit cells and add a physical energy constraint representing the homogenized elastic modulus of the micro-structure to the regression algorithm. We demonstrate this scheme for examples of one- and two-dimensional linear elastodynamics and show that the energy constraint improves the accuracy of the resulting peridynamic model.
Stable and accurate modeling of thin shells requires proper enforcement of all types of boundary conditions. Unfortunately, for Kirchhoff-Love shells, strong enforcement of Dirichlet boundary conditions is difficult because both functional and derivative boundary conditions must be applied. A popular alternative is to employ Nitsches method to weakly enforce all boundary conditions. However, while many Nitsche-based formulations have been proposed in the literature, they lack comprehensive error analyses and verifications. In fact, existing formulations are variationally inconsistent and yield sub-optimal convergence rates when used with common boundary condition specifications. In this paper, we present a novel Nitsche-based formulation for the linear Kirchhoff-Love shell that is provably stable and optimally convergent for general sets of admissible boundary conditions. To arrive at our formulation, we first present a framework for constructing Nitsches method for any abstract variational constrained minimization problem. We then apply this framework to the linear Kirchhoff-Love shell and, for the particular case of NURBS-based isogeometric analysis, we prove that the resulting formulation yields optimal convergence rates in both the shell energy norm and the standard $L^2$-norm. In the process, we derive the Euler-Lagrange equations for general sets of admissible boundary conditions and show that the Euler-Lagrange boundary conditions typically presented in the literature is incorrect. We verify our formulation by manufacturing solutions for a new shell obstacle course that encompasses flat, parabolic, hyperbolic, and elliptic geometric configurations. These manufactured solutions allow us to robustly measure the error across the entire shell in contrast with current best practices where displacement and stress errors are only measured at specific locations.
We introduce a technique to automatically convert local boundary conditions into nonlocal volume constraints for nonlocal Poissons and peridynamic models. The proposed strategy is based on the approximation of nonlocal Dirichlet or Neumann data with a local solution obtained by using available boundary, local data. The corresponding nonlocal solution converges quadratically to the local solution as the nonlocal horizon vanishes, making the proposed technique asymptotically compatible. The proposed conversion method does not have any geometry or dimensionality constraints and its computational cost is negligible, compared to the numerical solution of the nonlocal equation. The consistency of the method and its quadratic convergence with respect to the horizon is illustrated by several two-dimensional numerical experiments conducted by meshfree discretization for both the Poissons problem and the linear peridynamic solid model.