Pith. sign in

REVIEW 3 major objections 4 minor 31 references

Series solutions for clamped peridynamic beams using fourth-order eigenfunctions

T0 review · 3 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Peridynamic beam deflections converge ten times faster with fourth-order beam eigenfunctions.

desk verdict A solid Galerkin extension for peridynamic beams with a real but fixable gap in the statement of nonlocal boundary-layer conditions. read the letter →

arxiv 2412.09702 v1 pith:RHUUWVPU submitted 2024-12-12 physics.class-ph

Ziyu Wang , Ivan C. Christov This is my paper · ORCID
classification physics.class-ph
keywords peridynamicsEuler-Bernoullibeameigenfunctionsfourth-orderSturm-LiouvilleproblemGalerkinseriessolutionnonlocalelasticityconvergencerateclamped-clamped
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper shows that the deflection of a clamped-clamped peridynamic beam can be written as a series in the natural fourth-order eigenfunctions of the clamped beam problem, not as a Fourier sine series. For the static loads considered, these beam-function series converge so rapidly that a given error tolerance is reached with roughly ten times fewer terms than the previously used sine series. If correct, accurate analytical solutions of the peridynamic beam equation reduce to inverting a small linear system, making the nonlocal beam theory much easier to use in practice.

What carries the argument

The fourth-order beam eigenfunctions, defined by the clamped eigenvalue problem with transcendental relations $\coth\lambda^s_m-\cot\lambda^s_m=0$ and $\tanh\lambda^c_m+\tan\lambda^c_m=0$, form a complete orthonormal basis in $L^2[-1,1]$. Their orthonormality and completeness let the paper project the integrodifferential peridynamic equation onto this basis, producing a square, diagonally dominant but non-symmetric linear system for the series coefficients. The matrix entries are evaluated analytically in terms of sine and hyperbolic-sine integrals, and the whole procedure avoids the Castigliano-theorem trick previously used to replace clamped conditions with moment conditions.

What would settle it

Compute a numerical solution of the peridynamic beam equation for the Heaviside step load with horizon $\Delta=1$, using a consistent nonlocal boundary treatment (such as periodic extension or explicit boundary-layer corrections) outside $[-1,1]$, and compare it with the five-term beam-function series plotted in Fig. 7(b); a disagreement larger than the plot resolution would invalidate the series solution's reliance on whole-line extension.

Watch

Extended reading notes

Core claim

The central claim is that expanding the peridynamic beam equation's solution in the orthonormal beam eigenfunctions — the sine-like and cosine-like solutions of $d^4\psi/dX^4=\lambda^4\psi$ with $\psi(\pm 1)=\psi'(\pm 1)=0$ — yields series whose coefficients decay as $O(m^{-4})$ for odd modes and $O(m^{-5})$ for even modes under the loads examined. This is one order faster than the $O(m^{-3})$ decay of Fourier sine coefficients, and the resulting partial sums converge as roughly $m^{-3.5}$ in root-mean-square error versus $m^{-2.5}$ for the sine series. For the offset point-load example, ten terms of the beam-function series match a thousand terms of the sine series, so the higher-order basis is not only natural but also practically superior.

Load-bearing premise

The solution is represented by the eigenfunction series on the whole real line, so the finite-horizon integrals are evaluated without imposing additional nonlocal boundary conditions outside $[-1,1]$; if those boundary conditions change the physical solution materially, the computed deflections may not match the actual clamped-clamped peridynamic beam.

Editorial extensions

If this is right

  • Static peridynamic beam deflections under arbitrary loads can be obtained analytically by inverting a small, diagonally dominant matrix rather than summing a long trigonometric series.
  • The same basis could handle time-dependent loads by converting the problem into a system of ordinary differential equations for the time-varying coefficients.
  • The method yields a sequence of closed-form approximate solutions that converge uniformly to the classical beam solution as the horizon $\Delta\to 0^+$.
  • The approach sidesteps the need to construct a separate solution for supported beams and then impose clamping via Castigliano's theorem, because the eigenfunctions already satisfy the clamped boundary conditions.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The observed one-order-faster coefficient decay suggests the beam functions are better adapted to the fourth-order nonlocal operator; a similar advantage may hold for other high-order nonlocal boundary-value problems, such as peridynamic plate theories.
  • The paper implicitly extends the displacement beyond $[-1,1]$ as the analytic continuation of the eigenfunction series, which for large horizons (e.g., $\Delta=1$) is a physical assumption about nonlocal boundary conditions that could be checked against a fully discretized peridynamic simulation.
  • A directly testable extension is to apply the same beam-function basis to simply-supported or clamped-free beams, predicting a comparable convergence advantage over the standard Fourier basis.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper develops Galerkin-type series solutions for the static peridynamic beam equation derived by Yang et al. [28], using the clamped-clamped fourth-order ('beam') eigenfunctions (5). For odd, even, and step loads, the deflection is expanded as a series in these eigenfunctions, substituted into the integrodifferential equation (11), and projected to obtain linear systems for the coefficients with explicitly evaluated matrix entries (Appendix). The paper then compares the beam-function series with the Fourier sine series of Yang et al. [30] for an offset point load, reporting that coefficients decay one order faster and that a given error tolerance is reached with roughly ten times fewer terms. A final section studies the effect of the horizon size and notes convergence to classical beam theory as the horizon tends to zero.

Significance. If the mathematical setup is made fully precise, the paper offers a practically useful analytical tool: the Galerkin matrix entries are given in closed form, the Δ→0 limit correctly recovers the classical fourth-order operator, and the numerical illustrations suggest rapid coefficient decay. The comparison with the sine series is a genuine contribution to the literature on analytical peridynamic beam solutions. However, the significance is currently limited by an unstated nonlocal boundary condition, and the headline convergence claims rest on empirical coefficient fits rather than on a demonstrated asymptotic result. With those points addressed, the paper would be a solid addition to the analytical peridynamics literature.

major comments (3)
  1. [§3, Eq. (11); §4 and Appendix Eqs. (34)–(39)] Equation (11) is not a closed boundary-value problem on [−1,1]. For X within 2Δ of either boundary, the double integral samples U(X+H) and U(X+Ξ+H) at arguments outside [−1,1], so a nonlocal boundary condition (volume constraint) on a boundary layer of width at least 2Δ must be specified. The manuscript never states one. Instead, the derivations in §4 and the Appendix implicitly extend the series (7) to all real X using the analytic formulas (5); those eigenfunctions do not vanish for |X|>1, so this is not equivalent to a rigidly clamped boundary layer. This matters quantitatively: for Δ=0.1 the affected layer is 20% of the half-length, and Fig. 7 uses Δ=1, where the entire domain is affected. The abstract's claim of solving the 'clamped–clamped peridynamic beam' is therefore conditional on this unstated convention. Please state the nonlocal boundary condition explicitly, or reframe the claims as solutions of Eq. (11) with the specific eigenfunction-extension convention (5); ideally, validate against an independent numerical solution for a stated volume constraint.
  2. [§4.4, Fig. 6(c)] The RMSE is computed between the m-term truncation and the 1000-term truncation of the same series, for each method separately. This measures self-convergence of each series, not the error relative to the solution of a stipulated peridynamic boundary-value problem. The claim that the fourth-order eigenfunction series 'achieve a comparable precision with far fewer terms' should either be clearly labeled as a statement about self-convergence of the two truncations, or be backed by a comparison to a common reference solution, such as a high-resolution numerical solution of Eq. (11) with the same nonlocal boundary condition. Because the two series may converge to different functions if different extensions are used, the comparison in Fig. 6(a) alone is not sufficient.
  3. [§4.1–§4.3, Figs. 3–5] The central efficiency claim rests on the asserted coefficient decay rates O(m^{−4}) and O(m^{−5}). These rates are inferred by drawing reference lines of slope −4.2, −5.6, and −4.8 through the first ten computed coefficients; no asymptotic analysis is given. The fitted slopes do not exactly equal the claimed integer rates, so the 'one order faster' statement is not precisely quantified. Please provide a derivation or, failing that, a much longer coefficient sequence and an explicit statement of how the reference slopes were obtained (fitting range, least-squares procedure). The same concern applies to the convergence-rate slopes in Fig. 6(c).
minor comments (4)
  1. [Eq. (3a)] After nondimensionalization, the condition 't > 0' should read 'T > 0' to be consistent with the dimensionless variables introduced in Eq. (2).
  2. [Appendix Eq. (34b)] In the second term of Eq. (34b), the kernel is written with I_s_{1,n}, but the summation index is m; this should presumably be I_s_{1,m}.
  3. [Fig. 7 caption] The caption appears to use the placeholder 'ε' for the horizon, while the text and equations use Δ; please align the notation.
  4. [Appendix heading text] There is a typo: 'definining' should be 'defining' in the sentence introducing the integral evaluations.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: series coefficients and convergence rates are computed afresh from a published peridynamic beam model; the self-citations are background only.

full rationale

The derivation is self-contained against external inputs. Eq. (11) is adopted from Yang et al. [28], and the eigenfunctions (5) come from standard clamped-beam theory (Rayleigh [24]; Chandrasekhar [4]; Papanicolaou [17,19]). The paper does not fit any parameter to a target solution: the coefficients a_s_m and a_c_m are obtained by solving the Galerkin systems (14)/(22), i.e. Eqs. (17)/(25), whose matrix entries are evaluated analytically in the Appendix, and the plotted coefficient decay rates are outputs of that calculation. The claimed speed-up is a comparison against the independent Fourier-sine benchmark from Yang et al. [30], re-dimensionalized and truncated at 1000 terms; the RMSE curves are computed directly from those partial sums, not imposed. The only author-overlapping citations, [18] and [19], supply background on beam-function orthonormality and are also supported by non-self references, so they are not load-bearing. A modeling caveat does exist: for X within 2Delta of the endpoints, the integro-differential operator in Eq. (11) samples U outside [-1,1], and the paper never states a nonlocal volume constraint; its series substitution silently imposes a particular extension of the eigenfunctions. This is a well-posedness/modeling concern, not circularity, because the extension is an extra assumption rather than a reduction of the output to the input by construction.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

The central claim rests on standard spectral theory, the adopted peridynamic beam model, and an implicit treatment of the nonlocal boundary region; no free parameters are fitted to data.

assumptions (3)
  • standard math The fourth-order beam eigenfunctions form a complete orthonormal basis of L^2[-1,1], and the eigenfunction expansion (7) converges to the solution of the peridynamic boundary-value problem.
    Invoked in Section 2 via standard ODE theorems [5] and used throughout Section 4 to justify substituting series (7) into Eq. (11) and truncating at M terms.
  • domain assumption The peridynamic beam equation (8) from Yang et al. [28] is the correct governing equation for the nonlocal Euler-Bernoulli beam, and its dimensionless static form (11) is well-posed with clamped boundaries.
    The paper adopts Eq. (8) as given and solves it; it does not derive or validate the model.
  • domain assumption The displacement U is defined on the whole real line by the analytic continuation of the eigenfunction series, so finite-horizon integrals in Eq. (11) can be evaluated without additional nonlocal boundary-layer conditions outside [-1,1].
    This is unstated but necessary for evaluating the peridynamic operator near boundaries and for the Delta=1 case in Fig. 7.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Series solutions for clamped peridynamic beams using fourth-order eigenfunctions." pith.science (2026). https://pith.science/paper/RHUUWVPU

@misc{pith2026241209702,
  author       = {Pith},
  title        = {Pith review of: Series solutions for clamped peridynamic beams using fourth-order eigenfunctions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RHUUWVPU}},
  note         = {Machine review of arXiv:2412.09702}
}
read the original abstract

We propose an analytical approach to solving nonlocal generalizations of the Euler--Bernoulli beam. Specifically, we consider a version of the governing equation recently derived under the theory of peridynamics. We focus on the clamped--clamped case, employing the natural eigenfunctions of the fourth derivative subject to these boundary conditions. Static solutions under different loading conditions are obtained as series in these eigenfunctions. To demonstrate the utility of our proposed approach, we contrast the series solution in terms of fourth-order eigenfunctions to the previously obtained Fourier sine series solution. Our findings reveal that the series in fourth-order eigenfunctions achieve a given error tolerance (with respect to a reference solution) with ten times fewer terms than the sine series. The high level of accuracy of the fourth-order eigenfunction expansion is due to the fact that its expansion coefficients decay rapidly with the number of terms of the series, one order faster than the Fourier series in our examples.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

31 extracted references · 24 canonical work pages

  1. [28]

    Z Angew Math Mech (ZAMM) 102, e202200132

    Yang Z, Naumenko K, Altenbach H, Ma CC, Oterkus E, Oterkus S (2022) Some analytical solutions to peridynamic beam equations. Z Angew Math Mech (ZAMM) 102, e202200132. DOI 10.1002/zamm.202200132

    doi:10.1002/zamm.202200132 2022

  2. [30]

    Appl Sci 13, 10025

    Yang Z, Naumenko K, Ma CC, Chen Y (2023) Closed-form analytical solutions for the deflection of elastic beams in a peridynamic framework. Appl Sci 13, 10025. DOI 10.3390/app131810025

    work page doi:10.3390/app131810025 2023

  3. [1]

    In: Kienzler R, Ott I, Altenbach H (eds) Theories of Plates and Shells: Critical Review and New Applications, Springer, Berlin/Heidelberg, pp 1–12

    Altenbach H, Zhilin PA (2004) The theory of simple elastic shells. In: Kienzler R, Ott I, Altenbach H (eds) Theories of Plates and Shells: Critical Review and New Applications, Springer, Berlin/Heidelberg, pp 1–12. DOI 10.1007/978-3-540-39905-6 1

    work page doi:10.1007/978-3-540-39905-6 2004

  4. [2]

    Springer-Verlag, Berlin/Heidelberg

    Altenbach H, Maugin GA, Erofeev V (eds) (2011) Mechanics of Generalized Continua, Advanced Structured Materials, vol 7. Springer-Verlag, Berlin/Heidelberg. DOI 10.1007/978- 3-642-19219-7

    doi:10.1007/978- 2011

  5. [3]

    Comptes Rendus – Mecanique 346, 320–335

    Challamel N (2018) Static and dynamic behaviour of nonlocal elastic bar using integral strain-based and peridynamic models. Comptes Rendus – Mecanique 346, 320–335. DOI 10.1016/j.crme.2017.12.014

    work page doi:10.1016/j.crme.2017.12.014 2018

  6. [4]

    Dover Publications, New York, NY

    Chandrasekhar S (1981) Hydrodynamic and Hydromagnetic Stability. Dover Publications, New York, NY

    work page 1981

  7. [5]

    McGraw-Hill, New York, NY

    Coddington EA, Levinson N (1955) Theory of Ordinary Differential Equations. McGraw-Hill, New York, NY

    work page 1955

  8. [6]

    Int J Solids Struct 69-70, 152–168

    Diyaroglu C, Oterkus E, Oterkus S, Madenci E (2015) Peridynamics for bending of beams and plates with transverse shear deformation. Int J Solids Struct 69-70, 152–168. DOI 10.1016/j.ijsolstr.2015.04.040

    work page doi:10.1016/j.ijsolstr.2015.04.040 2015

Show all 31 references
  1. [7]

    Release 1.2.0 of 2024-03-15

    DLMF (2024) NIST Digital Library of Mathematical Functions. Release 1.2.0 of 2024-03-15. URL https://dlmf.nist.gov/, F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller, B. V. Saunders, H. S. Cohl, and M. A. McClain, eds

    2024

  2. [8]

    Theor Appl Mech Lett 8, 351–354

    Huang Z (2018) The singularity in the state-based peridynamic solution of uniaxial tension. Theor Appl Mech Lett 8, 351–354. DOI 10.1016/j.taml.2018.05.008

    2018 doi

  3. [9]

    J Comput Phys 493, 112466

    Kim KH, Bhalla APS, Griffith BE (2023) An immersed peridynamics model of fluid-structure interaction accounting for material damage and failure. J Comput Phys 493, 112466. DOI 10.1016/j.jcp.2023.112466

    2023

  4. [10]

    John Wiley & Sons, New York, NY

    Kraus H (1967) Thin Elastic Shells. John Wiley & Sons, New York, NY

    1967

  5. [11]

    Int J Solids Struct 49, 2887–2897

    Mikata Y (2012) Analytical solutions of peristatic and peridynamic problems for a 1D infinite rod. Int J Solids Struct 49, 2887–2897. DOI 10.1016/j.ijsolstr.2012.02.012

    2012 doi

  6. [12]

    Eur J Mech A/Solids 101, 104978

    Mikata Y (2023) Analytical solutions of peristatics and peridynamics for 3D isotropic mate- rials. Eur J Mech A/Solids 101, 104978. DOI 10.1016/j.euromechsol.2023.104978

    2023

  7. [13]

    Compos Struct 279, 114728

    Naumenko K, Eremeyev V A (2022) A non-linear direct peridynamics plate theory. Compos Struct 279, 114728. DOI 10.1016/j.compstruct.2021.114728

    2022

  8. [14]

    Math Mech Solids 22, 1639–1653

    Nishawala VV, Ostoja-Starzewski M (2017) Peristatic solutions for finite one- and two- dimensional systems. Math Mech Solids 22, 1639–1653. DOI 10.1177/1081286516641180

    2017 doi

  9. [15]

    Int J Solids Struct 51, 3177–3183

    O’Grady J, Foster J (2014) Peridynamic beams: A non-ordinary, state-based model. Int J Solids Struct 51, 3177–3183. DOI 10.1016/j.ijsolstr.2014.05.014

    2014 doi

  10. [16]

    Int J Solids Struct 51, 4572–4579

    O’Grady J, Foster J (2014) Peridynamic plates and flat shells: A non-ordinary, state-based model. Int J Solids Struct 51, 4572–4579. DOI 10.1016/j.ijsolstr.2014.09.003

    2014 doi

  11. [17]

    PhD thesis, University of Louisiana at Lafayette

    Papanicolaou NC (2003) A Galerkin Spectral Method for Fourth-Order Boundary Value Prob- lems. PhD thesis, University of Louisiana at Lafayette. URLhttps://www.proquest.com/ dissertations-theses/galerkin-spectral-method-fourth-order-boundary/ docview/305212891/se-2

    2003

  12. [18]

    J Phys: Conf Ser 2675, 012016

    Papanicolaou NC, Christov IC (2023) Orthonormal eigenfunction expansions for sixth- order boundary value problems. J Phys: Conf Ser 2675, 012016. DOI 10.1088/1742- 6596/2675/1/012016, 2308.00673

    2023 arXiv

  13. [19]

    Int J Num Meth Fluids 59, 945–965

    Papanicolaou NC, Christov CI, Homsy GM (2009) Galerkin technique based on beam func- tions in application to the parametric instability of thermal convection in a vertical slot. Int J Num Meth Fluids 59, 945–965. DOI 10.1002/fld.1845

    2009 doi

  14. [20]

    In: Altenbach H, Pouget J, Rousseau M, Collet B, Michelitsch T (eds) Generalized 18 Z

    Porubov A V, Osokina AE, Michelitsch TM (2018) Nonlocal approach to square lattice dy- namics. In: Altenbach H, Pouget J, Rousseau M, Collet B, Michelitsch T (eds) Generalized 18 Z. Wang and I. C. Christov Models and Non-classical Approaches in Complex Materials 1, Springer In...

    2018 doi

  15. [21]

    J Mech Phys Solids 48, 175–209

    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

    2000 doi

  16. [22]

    Int J Non-Linear Mech 40, 395–409

    Silling SA, Bobaru F (2005) Peridynamic modeling of membranes and fibers. Int J Non-Linear Mech 40, 395–409. DOI 10.1016/j.ijnonlinmec.2004.08.004

    2005 doi

  17. [23]

    Adv Appl Mech 44, 73–168

    Silling SA, Lehoucq RB (2010) Peridynamic theory of solid mechanics. Adv Appl Mech 44, 73–168. DOI 10.1016/S0065-2156(10)44002-8

    2010 doi

  18. [24]

    Macmillan and Co., London

    Strutt JW (1877) The Theory of Sound, vol 1. Macmillan and Co., London. URL http: //books.google.com/books?id=GyI5AAAAMAAJ&oe=UTF-8

  19. [25]

    Math Mech Solids 20, 998–1010

    Taylor M, Steigmann DJ (2015) A two-dimensional peridynamic model for thin plates. Math Mech Solids 20, 998–1010. DOI 10.1177/1081286513512925

    2015 doi

  20. [26]

    McGraw- Hill, New York, NY

    Timoshenko S, Woinowsky-Krieger S (1959) Theory of Plates and Shells, 2nd edn. McGraw- Hill, New York, NY

    1959

  21. [27]

    Comput Model Eng Sci 124, 527–544

    Yang Z, Oterkus E, Oterkus S (2020) A state-based peridynamic formulation for func- tionally graded Euler-Bernoulli beams. Comput Model Eng Sci 124, 527–544. DOI 10.32604/cmes.2020.010804

    2020

  22. [29]

    Mech Res Commun126, 104000

    Yang Z, Naumenko K, Ma CC, Altenbach H, Oterkus E, Oterkus S (2022) Some closed form series solutions to peridynamic plate equations. Mech Res Commun126, 104000. DOI j.mechrescom.2022.104000

    2022

  23. [31]

    Eur J Mech A/Solids 101, 105075

    Yang Z, Naumenko K, Ma CC, Oterkus E, Oterkus S (2023) Peridynamic analysis of curved elastic beams. Eur J Mech A/Solids 101, 105075. DOI 10.1016/j.euromechsol.2023.105075

    2023

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.