{"id":"dba5c024-21b3-47f3-9958-941f0a7dcd74","arxiv_id":"2507.13098","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A single hierarchy of beam models, from Euler-Bernoulli to a defect-capable micromorphic beam, is derived from one 3D higher-order elasticity framework, with simpler models recovered as singular limits.","lead":"This paper derives a family of beam models from a three-dimensional theory of materials with internal structure, ranging from classical beams to a model that includes material defects. The authors show that the simpler classical beam models emerge as limits of the most general one when penalty coefficients grow.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 39's e→∞ limit proof does not control the N_{jα} affine mode; under the stated one-point anchoring the limiting ODE system may not be exactly the semi-holonomic one.","rationale":"The reader's CONDITIONAL verdict is appropriate: the paper constructs a plausible hierarchy and the main proofs are sketched rather than rigorous. I agree that the Taylor ansatz in eqs. (9)–(10) and the boundary-anchoring assumptions are genuine limitations. However, the most load-bearing weakness for the central claim is specifically in the proof of the singular-limit theorems: the parameter-continuity argument contains a false Hilbert-space statement, and it does not verify that the e→∞ (or d→∞) limiting ODE system coincides with the semi-holonomic (or holonomic) system for all components, particularly N_{jα}. This is a distinct concern from the reader's Taylor-ansatz concern, so my agreement is partial. The concern is concrete and testable: a direct symbolic elimination in the Appendix B system, or a finite-e numerical sweep, would settle whether the limiting identification holds. Because the gap is fillable and the theorem is plausible, the verdict should remain CONDITIONAL rather than moving to REJECT or ACCEPT.","tokens_in":21570,"tokens_out":10885,"duration_ms":133769,"concrete_test":"Take the pure bending subsystem in Appendix B, eqs. (51). Formally eliminate eQ by differentiating the N-equations and adding them to the P-equations, then let e→∞. Write the resulting limiting ODE system and check whether, under the stated one-point anchoring of N_{21} and N_{12} (with all other data identical to eq. (52)), it forces N_{21} and N_{12} to be the same constants as in the semi-holonomic problem and reduces exactly to eq. (52). Alternatively, solve eqs. (51) numerically for e = 10, 10², 10³, …, with one-point anchoring, and measure ‖u_non−holo − u_semi−holo‖_L∞ and ‖N_{jα} − const‖_L∞; if these do not tend to zero, Theorem 39 as stated is false.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the singular-limit hierarchy: Theorem 39 (e→∞ sends non-holonomic to semi-holonomic) and Theorem 46 (d→∞ sends semi-holonomic to holonomic). The proof in §4.2 argues that after eliminating the penalty terms the equations become a first-order linear ODE system whose coefficients depend continuously on 1/e, and then asserts that the solution depends continuously on this parameter because C0(R*+,R) is a Hilbert space. That assertion is false: C0 with the sup norm is not a Hilbert space, and continuous parameter dependence of boundary-value problems also requires uniform well-posedness, which is not shown. More importantly, the limiting system is not fully identified. For the N_{jα} components, dividing eq. (31.4) by e and letting e→∞ gives only N_{jα,11}=0, i.e. an affine mode. The semi-holonomic model, by contrast, has N_{jα} constant, forced by N=∇P and eq. (18). The theorem's hypothesis that N_{jα} is 'fully specified on at least one point' does not determine the affine slope, and the assertion that the boundary data are 'consistent with ∇P=N' is not checked. If the natural boundary condition on N_{jα,1} is not compatible, the e→∞ solution may retain a non-holonomic affine component, so convergence to the semi-holonomic solution would fail. Theorem 46 has the same structure: the proof is deferred to the same parameter-continuity argument, and the limiting identification of P_{α} is not verified. The claim may be repairable via a direct energy estimate e‖N−∇P‖²≤C plus a careful boundary-layer analysis, but the proof as written does not establish it.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a family of one-dimensional beam models obtained by a formal dimensional reduction of a three-dimensional higher-order elasticity energy with three kinematic fields: the displacement u, the micro-distortion P, and the third-order tensor N. Three regimes are considered: holonomic (N=∇P, P=∇u), semi-holonomic (N=∇P), and non-holonomic (no constraint). The authors derive the Euler–Lagrange systems (31), (33), (34), discuss pure traction and bending sub-systems, and claim that the semi-holonomic and holonomic models emerge as singular limits of the non-holonomic model when the penalty coefficients e and d tend to infinity (Theorem 39 and Theorem 46). They also relate the holonomic and semi-holonomic models to higher-order Euler–Bernoulli and Timoshenko beams and interpret the fields N_i jk − N_i kj and N_i jk,l as dislocation and disclination densities.","tokens_in":21979,"tokens_out":15243,"duration_ms":161933,"significance":"If Theorems 39 and 46 are correct, the paper offers a valuable hierarchical unification: one non-holonomic one-dimensional continuum contains the classical beam models as singular limits, and dislocation and disclination fields are encoded in the same framework. The variational derivation from the three-dimensional energy is explicit, and the reduction to concrete ODE systems in the appendices is a useful reference. The identification of the classical Euler–Bernoulli and Timoshenko energies as special cases of (25) and (21) is clearly presented. The main mathematical novelty is the claimed singular-limit hierarchy, and the two central limit theorems are precisely where the proof is incomplete. The value of the paper will depend on whether the gaps identified below can be closed.","major_comments":[{"comment":"The proof of Theorem 39 is built on an invalid functional-analytic argument. The authors claim that the system (37), whose coefficients depend continuously on 1/e, can be viewed as a first-order linear ODE \"valued in C0(R*_+, R)\" and that, since this space is a Hilbert space, the solutions depend continuously on the parameter and are well-defined at 1/e=0. Two things are wrong: C0 with the sup norm is not a Hilbert space, and 0 is not in R*_+. More importantly, continuity of coefficients alone does not imply continuous dependence of solutions of a boundary-value problem on the parameter; one needs uniform well-posedness estimates, for example uniform bounds on the fundamental matrix and on the inverse of the boundary operator. No such estimates are provided. Since the convergence claims (40) rest on this step, Theorem 39 is not established.","section":"§4.2, eqs. (36)–(38) and Theorem 39"},{"comment":"The identification of the limiting ODE system is incomplete for the components N_i jα. Dividing eq. (31.4) by e and passing to the limit gives only N_i jα,11 = O(1/e), so the limit is affine in X1, whereas in the semi-holonomic target N_i jα is constant by eq. (18). The hypothesis that N_i jα is fully specified at one boundary point fixes the intercept but not the slope N_i jα,1; the natural boundary condition in eq. (31) involves (c + eℓ^4/12)N_i jα,1, whose limit is not controlled. Consequently the claimed convergence ‖N − ∇P‖_{L∞} → 0 in the proof, and the assertion that the limiting system is exactly eq. (33), are not justified.","section":"§4.2, Theorem 39 (limit e→∞, transverse components of N)"},{"comment":"The proof sketch for Theorem 39 does not derive the limiting system. Eq. (36) is asserted without showing the algebra that eliminates P; the reduction to the first-order form (37) is not displayed; and eq. (42) is obtained by an unspecified differentiation-and-subtraction procedure. The passage from (42) to (44) assumes N → ∇P, which is essentially the conclusion to be proved. A rigorous proof would require a direct energy estimate controlling e‖N − ∇P‖ in a suitable norm, followed by passage to the limit in the weak formulation; no such estimate is given.","section":"§4.2, eqs. (41)–(44)"},{"comment":"Theorem 46 is not a consequence of \"the exact same arguments\". Eq. (45) gives O(1/d) relations for u_i,1 − P_i 1 and for N_i 1α − P_i α,1, but the holonomic target (34) is a sixth-order scalar problem for u, and the convergence of the remaining components P_i α to u_i,α is not shown. In particular, the dℓ^4/12 terms in eq. (33) involve P_i α,11; without uniform a priori bounds on the relevant derivatives as d → ∞, the limiting identification of P_i α is open. The boundary conditions also need to be shown to pass to the limit, which is not done.","section":"§4.2, Theorem 46 (limit d→∞)"},{"comment":"The reduction from the three-dimensional energy to the one-dimensional beam energies is a formal Taylor truncation. No error estimate is provided for the approximations u_i,α ≈ P_i α, P_i j,α ≈ N_i jα, and N_i jk,α ≈ 0, so the energies (14), (20), and (24) are a modeling ansatz rather than a proved asymptotic limit. This should be stated explicitly; it does not affect the internal consistency of the one-dimensional ODE systems, but it limits the claim that the beam models are derived from the three-dimensional theory.","section":"§2, eqs. (9)–(10)"}],"minor_comments":[{"comment":"The definition \"2 [sym P]_i^j := P_i^j − P_j^i\" is the negative of the symmetrization used in eq. (6); in the linearisation of (Id+P)^T(Id+P) the sign must be plus.","section":"Eq. (4)"},{"comment":"The surjectivity of the homogenised boundary/bulk force correspondence is asserted without proof or reference; it should be stated as an assumption, since the treatment of the one-dimensional model as independent depends on it.","section":"After eq. (17)"},{"comment":"There are several typos that should be corrected, including \"freODEm\" in the second bullet of Section 1.2, \"remainder fo this article\" at the start of Section 3, and \"nono-holonomic\" in the proof of Theorem 46.","section":"§1.2 and §3"},{"comment":"The sentence \"This shows that ‖N − ∇P‖_{L∞} → 0\" is stronger than what eq. (38) establishes: the displayed bound concerns only the i,j1 components N_i j1 − P_i j,1. The transverse components N_i jα − P_i j,α are not controlled by this argument.","section":"§4.2, eq. (38)"},{"comment":"The notation C0(R*_+,R) is ambiguous: if the parameter 1/e is meant to include 0, the domain should be a compact interval such as [0,ε], and the space should be specified with the sup norm and identified as a Banach space, not a Hilbert space.","section":"§4.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript appears to be a conference-proceedings contribution, and the two central theorems are advertised as proofs but are not rigorous in the current form. The Hilbert-space sentence is a clear red flag. The affine-mode issue for N_i jα is the most serious technical gap: if it cannot be closed, the e → ∞ hierarchy claim should be weakened. The modeling ansatz in Section 2 is acceptable for a mechanics paper, but for a math.AP venue the formal nature of the reduction should be made explicit. I recommend major revision rather than rejection because the framework and the variational derivations are coherent and the flaws are potentially repairable with direct energy estimates."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe thing to know: this paper is worth a referee's time, but the singular-limit theorems it advertises are not proved as written. The genuinely new piece is the non-holonomic 1D model, where N carries both torsion (dislocations) and curvature (disclinations) in a beam, and the idea that the semi-holonomic and holonomic beam models are e→∞ and d→∞ limits of that model. That hierarchy, if it works, is a useful unification for the generalized-continuum beam literature. The holonomic and semi-holonomic special cases reduce honestly to known gradient beam models (Lurie–Solyaev, Wang et al., Papargyri-Beskou et al.), and the paper says so, which I respect.\n\nWhat it does well: the 3D-to-1D reduction via the Taylor ansatz is clearly laid out, the connected-component analysis in §4.1 is a nice way to organize the 39 fields, and the identification of dislocations with N^j_k - N^k_j and disclinations with N^j_{kl} - N^j_{lk} is clean. The simplified pure-traction and pure-bending ODE systems in the appendices are useful anchors.\n\nWhere it is soft:\n\n1. The proof of Theorem 39 is not a proof. The step 'C0(R*+,R) is a Hilbert space' is false, and continuous parameter dependence of solutions on 1/e needs uniform well-posedness, which is not shown. More substantively, the limiting identification is incomplete: from eq. (31.4), dividing by e and sending e→∞ gives only N^j_{α,11}=0, i.e. an affine mode. The semi-holonomic model has N^j_α constant, forced by N=∇P and eq. (18). The condition that N^j_α is specified at one point does not kill the affine slope, and the compatibility of the boundary data with ∇P=N is asserted, not checked. If the natural boundary condition on N^j_{α,1} is incompatible, the limit may retain a non-holonomic affine component. Theorem 46 has the same structure; the limiting identification of P_α is deferred to the same continuity argument and is not verified.\n\n2. The 1D reduction rests on the transversal Taylor ansatz (9)–(10), with N^j_{k,α}≈0. If N varies across the section, the energies and limit theorems do not follow. The paper says 'safely neglect' without a small-parameter estimate.\n\n3. The energy simplifications in (21) and (25) need boundary anchoring of N (or P). That is stated, but it restricts the claimed reductions to problems with those anchoring conditions.\n\nAre these fixable? Probably. The e→∞ claim could be repaired by a direct energy estimate e‖N−∇P‖²≤C plus a boundary-layer analysis, and the C0 error is a side issue. But as written, the central theorems are not established.\n\nVerdict: conditional accept for a workshop proceedings; for a journal, major revision. The paper deserves a serious referee, not a desk reject, because the framework and the hierarchy claim are original and the gaps look fillable. I'd read a revised version, and I'd cite the non-holonomic model once it is on firmer ground.\n\n[Signature]","headline":"A genuinely new non-holonomic beam model with a plausible but unproved singular-limit hierarchy; the limit theorems need real work before the central claim is established.","tokens_in":22523,"tokens_out":4137,"would_cite":false,"duration_ms":39287,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["74K10","74A35","74B05","35B25"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that Euler–Bernoulli and Timoshenko beam models are singular limits of one non-holonomic continuum model carrying dislocations and disclinations.","keywords":["beams","generalised continua","higher-order elasticity","micromorphic elasticity","dislocations","disclinations","defects","singular limits"],"falsifier":"A concrete check: solve the non-holonomic system (31) with a transverse-dependent $N$, for instance $N^1_{12}(X^1,X^2)=N^0 + \\varepsilon X^2$, so that assumption (10) is violated, and test whether $\\|u_{\\mathrm{non-holo}}-u_{\\mathrm{semi-holo}}\\|_{L^\\infty}$ still tends to zero as $e\\to\\infty$; Theorem 39 predicts it should, and a failure would show the Taylor ansatz is load-bearing.","tokens_in":21379,"feed_emoji":"📐","tokens_out":9764,"duration_ms":97223,"temperature":0.7,"pith_summary":"This paper tries to establish that the classical beam hierarchy is one family: a higher-order Euler–Bernoulli model, a generalised Timoshenko model, and a fully relaxed model that carries both dislocations and disclinations all descend from a single three-dimensional higher-order elasticity framework. Starting from three kinematic fields — the displacement $u$, the micro-distortion $P$, and a third-order tensor $N$ — the authors impose a truncated Taylor ansatz across the cross-section to obtain one-dimensional beam energies and static equilibrium systems. The load-bearing results are Theorem 39, stating that as the penalty coefficient $e$ enforcing $\\nabla P = N$ grows, the non-holonomic solution converges in $L^\\infty$ to the semi-holonomic solution, and Theorem 46, stating that as $d$ enforcing $\\nabla u = P$ grows, the semi-holonomic solution converges to the holonomic solution. If correct, classical beam theories gain a defect-theoretic parent model and a precise statement of how they are recovered as limits.","feed_headline":"A single beam model contains Euler-Bernoulli and Timoshenko as limits","feed_subtitle":"Raising two penalty coefficients turns a defect-rich continuum model first into Timoshenko, then Euler-Bernoulli.","key_machinery":"The central object is the third-order tensor field $N$, added to the classical displacement $u$ and the micromorphic micro-distortion $P$, with energy (6) in which the terms $d\\|\\nabla u - P\\|^2$ and $e\\|\\nabla P - N\\|^2$ act as tunable penalties. The argument runs through a 1D-in-3D reduction: assumptions (9)–(10), namely $u_{i,\\alpha}\\approx P_{i\\alpha}$, $P_{ij,\\alpha}\\approx N_{ij\\alpha}$ and $N_{ijk,\\alpha}\\approx 0$, replace the cross-sectional dependence by a Taylor ansatz and turn the three-dimensional energy into a unidimensional system of ordinary differential equations with only longitudinal derivatives. The proof of the hierarchy works by rewriting the non-holonomic Euler–Lagrange system as a first-order ODE system continuously parametrised by $1/e$ and $1/d$, so the limit $e\\to\\infty$ or $d\\to\\infty$ is a regular point of the parametrisation; the limiting system is then read off as the semi-holonomic or holonomic equations.","core_discovery":"The central claim is that the holonomic, semi-holonomic and non-holonomic beam models are not independent theories but three levels of one framework, connected by singular limits: Theorem 39 shows $\\lim_{e\\to\\infty}\\|u_{\\mathrm{non-holo}}-u_{\\mathrm{semi-holo}}\\|_{L^\\infty}=0$ (and similarly for $P$ and $N$), while Theorem 46 shows $\\lim_{d\\to\\infty}\\|u_{\\mathrm{semi-holo}}-u_{\\mathrm{holo}}\\|_{L^\\infty}=0$. The coefficient $e$ multiplies $\\|\\nabla P - N\\|^2$ in the energy, and $d$ multiplies $\\|\\nabla u - P\\|^2$, so increasing each coefficient enforces the corresponding constraint. In the holonomic regime the model reduces to a higher-order Euler–Bernoulli beam with energy $a u_{,1}u_{,1}+b u_{,11}u_{,11}+c u_{,111}u_{,111}$; in the semi-holonomic regime it generalises the Timoshenko beam; and in the non-holonomic regime the unconstrained $N$ and $P$ fields allow both dislocation densities and disclination densities to appear from the kinematics alone.","pith_inferences":["If the same penalty structure is carried into the nonlinear regime, the limit theorems suggest a physical interpretation: the measured response of a real beam should interpolate between Timoshenko and Euler–Bernoulli behaviour as micro-structural stiffness grows, which would give a protocol for calibrating $d$ and $e$ from experiments.","The Taylor ansatz (9)–(10) should admit a two-dimensional analogue, producing a hierarchy of plate theories with dislocations and disclinations; testing that extension would show whether the beam result is a special case of a general dimensional-reduction principle.","The consistency conditions in Theorems 39 and 46 require the boundary data to match $\\nabla P=N$ or $\\nabla u=P$; without that matching, boundary layers are likely to appear, so the $L^\\infty$ convergence as $e$ or $d$ grows may require an additional surface-energy correction that the paper does not treat."],"forward_implications":["When $e\\to\\infty$, the non-holonomic beam equations reduce to the semi-holonomic system; when $d\\to\\infty$, that system reduces to the holonomic system, so engineering beam models form a single chain rather than separate theories.","Because the non-holonomic model leaves $N$ free, a purely kinematic bending problem can excite both a dislocation density $N^1_{12}-N^1_{21}$ and a disclination density $N^1_{12}$ through the $d\\ell^4/12$ curl-type coupling, so defects appear without extra constitutive assumptions.","Choosing $c=0$ in the holonomic energy yields the standard Euler–Bernoulli beam equation, and choosing $c=0$, $\\ell=0$, $f_2=0$ in the semi-holonomic energy yields the Timoshenko beam energy, placing the textbook models inside the hierarchy as parameter limits.","The simplified pure-traction and pure-bending subsystems give explicit ODE systems that can be solved numerically, so the predicted hierarchy is testable in closed form for basic load cases."],"supporting_citations":[{"why":"Supplies the three-dimensional higher-order elasticity model with $u$, $P$ and $N$ from which all beam energies are derived.","marker":"[7]"},{"why":"Provides the relaxed micromorphic model and the curl-based dislocation density that the homogenised $d\\ell^4/12$ term is compared to.","marker":"[29]"},{"why":"Defines the classical micromorphic theory whose kinematics the semi-holonomic model matches when $\\nabla P=N$.","marker":"[25]"},{"why":"Introduces the Timoshenko beam that the semi-holonomic model generalises.","marker":"[36]"},{"why":"Offers a strain-gradient Timoshenko beam model that reduces to the semi-holonomic energy for a specific parameter choice.","marker":"[41]"},{"why":"Provides a gradient-elasticity extension of Euler–Bernoulli beams with which the holonomic higher-order model is compared.","marker":"[24]"}],"fun_headline_variants":["One beam model embeds Euler-Bernoulli and Timoshenko","Euler-Bernoulli and Timoshenko emerge as limits of one beam theory","Unified beam framework links defect-rich model to classical ones","Three beam regimes collapse into one via penalty limits","Single continuum yields Euler-Bernoulli, Timoshenko, and defect models"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire reduction depends on assuming that the kinematic fields vary exactly linearly across the thin cross-section, with the third-order field $N$ constant across the section; if $N$ actually varies across the section, the simplified beam energies and both limit theorems no longer follow.","fun_headline_variants_meta":{"raw":{"variants":["One beam model embeds Euler-Bernoulli and Timoshenko","Euler-Bernoulli and Timoshenko emerge as limits of one beam theory","Unified beam framework links defect-rich model to classical ones","Three beam regimes collapse into one via penalty limits","Single continuum yields Euler-Bernoulli, Timoshenko, and defect models"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000858,"raw_usage":{"total_tokens":3750,"prompt_tokens":998,"completion_tokens":2752,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":614,"completion_tokens_details":{"reasoning_tokens":2662}},"tokens_in":614,"tokens_out":2752,"duration_ms":24907,"temperature":1.0,"reasoning_tokens":2662,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:32:38.126197+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete check: solve the non-holonomic system (31) with a transverse-dependent $N$, for instance $N^1_{12}(X^1,X^2)=N^0 + \\varepsilon X^2$, so that assumption (10) is violated, and test whether $\\|u_{\\mathrm{non-holo}}-u_{\\mathrm{semi-holo}}\\|_{L^\\infty}$ still tends to zero as $e\\to\\infty$; Theorem 39 predicts it should, and a failure would show the Taylor ansatz is load-bearing.","supporting_citations":[{"cited_title":"T wo-Scale Geomet ric Modelling for Defective Media","cited_arxiv_id":null,"evidence_quote":"Supplies the three-dimensional higher-order elasticity model with $u$, $P$ and $N$ from which all beam energies are derived."},{"cited_title":"A Unifying Perspective: The Relaxed Linear Micromorphic Continuum","cited_arxiv_id":null,"evidence_quote":"Provides the relaxed micromorphic model and the curl-based dislocation density that the homogenised $d\\ell^4/12$ term is compared to."},{"cited_title":"Micro-Structure in Linear Elasticity","cited_arxiv_id":null,"evidence_quote":"Defines the classical micromorphic theory whose kinematics the semi-holonomic model matches when $\\nabla P=N$."},{"cited_title":"On the Correction for Shear of the Diff erential Equation for Transverse Vibra- tions of Prismatic Bars","cited_arxiv_id":null,"evidence_quote":"Introduces the Timoshenko beam that the semi-holonomic model generalises."},{"cited_title":"A Micro Scale Timoshenko B eam Model Based on Strain Gra- dient Elasticity Theory","cited_arxiv_id":null,"evidence_quote":"Offers a strain-gradient Timoshenko beam model that reduces to the semi-holonomic energy for a specific parameter choice."},{"cited_title":"Revisiting Bending Theories f or Gradient Elastic Beams","cited_arxiv_id":null,"evidence_quote":"Provides a gradient-elasticity extension of Euler–Bernoulli beams with which the holonomic higher-order model is compared."}],"review_version":1}