Static and Vibration Analysis of Two-phase Peridynamic Rods or Strings under General Boundary Conditions — Exact, FDM and FEM Approaches

Static and dynamic analysis of two-phase peridynamic one-dimensional systems is thoroughly investigated in this paper for various boundary conditions by exact and numerical approaches. Such a nonlocal study concerns the axial or torsional responses of two-phase peridynamic rods or the behavior of a two-phase peridynamic string. For a normalized exponential kernel, the integro-differential equation of the two-phase peridynamic problem can be converted as a linear higher-order differential equation with constant coefficients. Exact static solutions of fixed–fixed and fixed–free two-phase peridynamic axial rod (TPPAR) under uniform distributed loading and concentrated force are analytically derived. It is shown that the pure peridynamic rod problem can be ill-posed due to a conflict between the natural boundary condition and the peridynamic constitutive boundary condition. The two-phase peridynamic problem may be justified from physical arguments, and is also seen as a regularized peridynamic problem. The exact eigen-frequencies of the two-phase peridynamic rods are also analytically obtained for various boundary conditions, including fixed–free, fixed–fixed or free–free boundary conditions. These theoretical approaches are corroborated by numerical investigations based on Finite Difference Method (FDM), Finite Element Method (FEM), exact FEM and variational iterative methods. Two FDM schemes are presented for the statics and the vibration of two-phase peridynamic rods, with convergence orders of 1 and 2, respectively. Moreover, three different FE procedures are proposed as: (1) FEM for direct solving the integro-differential equations of two-phase peridynamic one-dimensional system, (2) FEM for discretizing the pure differential equations, and (3) Exact FEM to obtain both the exact stiffness matrix and exact dynamic stiffness matrix of two-phase peridynamic one-dimensional system, based on Eisenberger’s methodology. The order of convergence for the first and the second proposed FEM is equal to 2 and 4, according to the type of selected finite elements. Moreover, the third proposed FEM based on the exact shape function of the peridynamic problem gives the exact solution from one single FEM element. In addition, the variational iterative method of Stodola–Vianello is successfully applied to the two-phase peridynamic eigenvalue problems.