{"id":"6d2cbb2d-4577-4e68-827a-2af06917044e","arxiv_id":"2412.09702","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Clamped-clamped peridynamic beam deflections are expanded in fourth-order beam eigenfunctions, giving series solutions that converge about ten times faster than the prior sine-series approach.","lead":"The paper introduces an analytical series method for the clamped-clamped peridynamic beam, using beam vibration eigenfunctions that naturally satisfy the clamped ends. The new series converges roughly ten times faster than the existing Fourier sine series, so accurate solutions need far fewer terms.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"For finite horizons the clamped condition is imposed only at the endpoints; the solution depends on an unspecified extension of the eigenfunctions outside [-1,1], so the headlined comparison may not describe the clamped peridynamic beam.","rationale":"The reader's weakest assumption identifies exactly the load-bearing gap: the nonlocal boundary layer is not specified, and the series calculation implicitly defines an extension of the displacement outside [-1,1] through the eigenfunctions themselves. I agree that this is the weakest point. The rest of the paper is internally consistent: the matrix formulas in the Appendix, including the factors of Δ, reduce correctly to the classical diagonal limit as Δ→0, and the numerical comparisons are internally coherent. The faster decay of beam-function coefficients relative to sine-series coefficients for clamped problems is plausible and is demonstrated consistently in Figs. 3-6, although it is empirical rather than proven. However, all of this evidence addresses the solution of Eq. (11) under the paper's implicit extension, not the clamped-clamped peridynamic beam as a physical boundary-value problem. Since the reader already made the verdict CONDITIONAL on this assumption, my read does not change the verdict. The concrete test would settle whether the implicit extension is harmless for the reported examples; if the difference is negligible, the conditional can be relaxed, and if it is large, the paper should be revised to state the extension explicitly and re-run the comparisons under a properly posed nonlocal clamped condition.","tokens_in":14202,"tokens_out":15433,"duration_ms":167038,"concrete_test":"Discretize Eq. (11) on the extended interval [-1-2Δ,1+2Δ] using composite quadrature (or finite differences) with the peridynamic clamped volume constraint U=0 on the appended layers, for loads P(X)=-X and P(X)=H(X) at Δ=0.1 and Δ=1. Compare the interior deflection with the paper's ten-term beam-function series and compute the L2 difference; if it exceeds the apparent error between the one-term and ten-term series shown in Figs. 3 and 5, the implicit-extension convention is quantitatively important and the 'clamped-clamped' identification is not yet established. Then recompute the beam-function coefficients of the corrected solution to test whether the claimed O(m^-4)/O(m^-5) decay survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"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) 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, Section 4 and the Appendix substitute the series (7) and evaluate shifted terms via the analytic formulas (5) for all real X, thereby imposing a particular extension of the clamped eigenfunctions beyond the beam. That extension is a modeling choice, not the usual peridynamic clamped condition of a rigidly fixed boundary layer: ψ_s^m and ψ_c^m vanish with zero slope at X=±1, but they do not vanish for |X|>1, and their hyperbolic parts continue to vary. For Δ=0.1 the affected layer is 20% of the half-length, and Fig. 7 uses Δ=1, where the extension affects the whole domain. No independent numerical solution with a stated volume constraint, and no comparison to a different series with a different extension, is provided. Because the abstract reports solutions for the clamped-clamped peridynamic beam and quantifies a ten-fold speedup, the central claim is conditional on this unstated convention; a different, equally plausible nonlocal clamping rule could change both the deflections and the coefficient-decay comparison.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":14407,"tokens_out":6336,"duration_ms":65236,"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":[{"comment":"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.","section":"§3, Eq. (11); §4 and Appendix Eqs. (34)–(39)"},{"comment":"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.","section":"§4.4, Fig. 6(c)"},{"comment":"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).","section":"§4.1–§4.3, Figs. 3–5"}],"minor_comments":[{"comment":"After nondimensionalization, the condition 't > 0' should read 'T > 0' to be consistent with the dimensionless variables introduced in Eq. (2).","section":"Eq. (3a)"},{"comment":"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}.","section":"Appendix Eq. (34b)"},{"comment":"The caption appears to use the placeholder 'ε' for the horizon, while the text and equations use Δ; please align the notation.","section":"Fig. 7 caption"},{"comment":"There is a typo: 'definining' should be 'defining' in the sentence introducing the integral evaluations.","section":"Appendix heading text"}],"recommendation":"major_revision","confidential_remarks":"The algebraic core of the paper is sound and the authors have done the community a service by providing explicit formulas for the Galerkin matrices. The main reservation is that the paper claims to solve a clamped peridynamic beam without specifying the nonlocal volume constraint; this is a modeling issue, not a mere presentation problem. I would not reject because the construct is well defined once the extension is stated, but the authors should be required to either adopt a standard clamped boundary layer and recompute the reported results, or explicitly limit the claims to the extension implied by Eq. (5). The convergence-rate claims also need more support than a visual fit to ten coefficients."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper does something concrete: it expands the clamped peridynamic beam deflection in the natural fourth-order beam eigenfunctions, works out analytic Galerkin matrix formulas, and shows faster convergence than the prior sine-series approach. That is a genuine, reproducible improvement. The derivation checks out—substitution, orthogonality, and the Δ→0 limit all behave as expected, and the explicit formulas in Eqs. (36) and (39) are useful for benchmarking and truncation analysis.\n\nThe soft spot the reader flagged is real: the boundary-value problem is not fully closed. Equation (11) involves integrals that sample U at arguments beyond [-1,1] whenever X is within 2Δ of the boundary, but the manuscript never states a nonlocal volume constraint on that boundary layer. The series (7) implicitly provides an extension through the analytic eigenfunction formulas, which is a modeling choice. For Δ=0.1 the affected layer is 20% of the half-length; for Δ=1 the extension affects the whole domain. The paper acknowledges that large Δ is \"a mathematical exercise,\" but the abstract's \"clamped-clamped peridynamic beam\" language is stronger than what is actually solved. This is fixable by stating the extension explicitly or by comparing to a rigidly fixed boundary layer, and it would strengthen the paper.\n\nOther concerns are minor. The claimed decay rates (O(m^-4), O(m^-5)) are inferred from reference lines fit to the first ten coefficients, not proven; and the RMSE is measured against a 1000-term truncation of the same series, not an independent numerical solution. Neither threatens the core message, especially since the comparison to Yang et al.'s sine series is direct and fair.\n\nWho this is for: researchers working on analytical solutions in peridynamics or nonlocal beam models, and anyone who needs benchmark solutions for nonlocal mechanics. It deserves a serious referee. With a moderate revision clarifying the boundary-layer convention and a more careful statement of the decay-rate claims, it would be a solid contribution.\n\nI would send it to peer review.","headline":"A solid Galerkin extension for peridynamic beams with a real but fixable gap in the statement of nonlocal boundary-layer conditions.","tokens_in":14968,"tokens_out":3182,"would_cite":false,"duration_ms":35102,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Peridynamic beam deflections converge ten times faster with fourth-order beam eigenfunctions.","keywords":["peridynamics","Euler-Bernoulli beam","beam eigenfunctions","fourth-order Sturm-Liouville problem","Galerkin series solution","nonlocal elasticity","convergence rate","clamped-clamped beam"],"falsifier":"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.","tokens_in":13960,"feed_emoji":"📐","tokens_out":4001,"duration_ms":38174,"temperature":0.7,"pith_summary":"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.","feed_headline":"Peridynamic beam series converge ten times faster","feed_subtitle":"Fourth-order beam eigenfunctions cut the terms needed to match a 1000-term sine series to about ten.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the peridynamic beam equation (8) that the paper solves.","marker":"[28]"},{"why":"Provides the Fourier sine series solution and the point-load setup used as the comparison baseline.","marker":"[30]"},{"why":"Defines the fourth-order beam eigenfunctions and establishes their orthonormality.","marker":"[17]"},{"why":"Also defines and uses these eigenfunctions, serving as the source for the spectral properties the paper relies on.","marker":"[19]"},{"why":"Provides the standard theorem that the eigenfunctions form a complete set in $L^2$, justifying the series expansion.","marker":"[5]"}],"fun_headline_variants":["Beam eigenfunctions cut series terms tenfold","Peridynamic beams solved with 10x fewer terms","Tenfold faster convergence for peridynamic beams","Fourth-order eigenfunctions beat sine series in beams","Series for beams: 10x less terms with better basis"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Beam eigenfunctions cut series terms tenfold","Peridynamic beams solved with 10x fewer terms","Tenfold faster convergence for peridynamic beams","Fourth-order eigenfunctions beat sine series in beams","Series for beams: 10x less terms with better basis"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00021,"raw_usage":{"total_tokens":1376,"prompt_tokens":877,"completion_tokens":499,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":493,"completion_tokens_details":{"reasoning_tokens":423}},"tokens_in":493,"tokens_out":499,"duration_ms":4741,"temperature":1.0,"reasoning_tokens":423,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T16:51:00.315083+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Appl Sci 13, 10025","cited_arxiv_id":null,"evidence_quote":"Provides the Fourier sine series solution and the point-load setup used as the comparison baseline."},{"cited_title":"PhD thesis, University of Louisiana at Lafayette","cited_arxiv_id":null,"evidence_quote":"Defines the fourth-order beam eigenfunctions and establishes their orthonormality."},{"cited_title":"Int J Num Meth Fluids 59, 945–965","cited_arxiv_id":null,"evidence_quote":"Also defines and uses these eigenfunctions, serving as the source for the spectral properties the paper relies on."},{"cited_title":"McGraw-Hill, New York, NY","cited_arxiv_id":null,"evidence_quote":"Provides the standard theorem that the eigenfunctions form a complete set in $L^2$, justifying the series expansion."}],"review_version":1}