pith. sign in

arxiv: 1907.02165 · v1 · pith:QA7KPJEZnew · submitted 2019-07-04 · 🧮 math.NA · cs.NA

Analysis and numerical simulation of the nonlinear beam equation with moving ends

Natanael Quintino , Mauro Rincon This is my paper

Pith reviewed 2026-05-25 09:45 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords beam equationmoving boundariesHermite polynomialsfinite element methodstability analysisconvergenceNewmark methodenergy decay
0
0 comments X

The pith

Hermite polynomial finite elements achieve quadratic convergence for the beam equation with moving ends.

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

The paper develops a numerical method for the small-amplitude motion of an elastic beam with internal damping when the ends move over time. It applies Hermite polynomials to discretize the spatial operator on the time-dependent domain and combines this with finite-difference time stepping. Stability and error estimates are derived for the resulting semi-discrete and fully-discrete schemes, establishing quadratic convergence in both space and time. Numerical tests in one and two dimensions confirm the predicted rates and demonstrate uniform energy decay.

Core claim

An efficient numerical method is constructed to solve this moving boundary problem. The stability and convergence of the method is studied, and the errors of both the semi-discrete and fully-discrete schemes are derived, using Hermite's polynomials as a base function, proving that the method has order of quadratic convergence in space and time. Numerical simulations using the finite element method associated with the finite difference method are employed for one-dimensional and two-dimensional cases, validating the theoretical results through comparison tables and showing the uniform decay rate for energy.

What carries the argument

Hermite finite-element space applied directly to the moving-boundary beam problem, combined with Newmark time integration.

If this is right

  • The semi-discrete and fully-discrete schemes remain stable for the moving-boundary formulation.
  • Quadratic convergence holds for both one- and two-dimensional problems.
  • Uniform energy decay rates are preserved under the discrete scheme.
  • Approximate solutions match exact solutions at the predicted rates in validation tests.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar direct discretizations may apply to other hyperbolic problems whose domains change with time.
  • The method could be extended to test more general nonlinearities or different damping terms.
  • Engineering models of beams with time-varying supports become directly computable without coordinate transformations.

Load-bearing premise

The moving-boundary problem can be discretized directly with a standard Hermite finite-element space while preserving the stability and approximation properties that hold for fixed domains.

What would settle it

A computed convergence rate below quadratic order in a test with rapid boundary motion would contradict the derived error estimates.

Figures

Figures reproduced from arXiv: 1907.02165 by Mauro Rincon, Natanael Quintino.

Figure 1
Figure 1. Figure 1: Asymptotic behavior of energy T ∗ = 770, in example S1, and T ∗ = 560, in example S2 , for both moving borders. In [PITH_FULL_IMAGE:figures/full_fig_p022_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Moving boundaries 5.2 Homogeneos Solution Now we will present the evolution of the numerical solution of the problem (6), obtained from the (18), by considering f(y, t) = 0 and the initial conditions from taking t = 0 in exact functions of [PITH_FULL_IMAGE:figures/full_fig_p022_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Homogeneous solutions for example S1 (a) Under the effect of the initial velocity (b) Post-effect of initial velocity [PITH_FULL_IMAGE:figures/full_fig_p023_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Homogeneous solutions for example S2 [PITH_FULL_IMAGE:figures/full_fig_p023_4.png] view at source ↗
read the original abstract

The numerical analysis for the small amplitude motion of an elastic beam with internal damping is investigated in domain with moving ends. An efficient numerical method is constructed to solve this moving boundary problem. The stability and convergence of the method is studied, and the errors of both the semi-discrete and fully-discrete schemes are derived, using Hermite's polynomials as a base function, proving that the method has order of quadratic convergence in space and time. Numerical simulations using the finite element method associated with the finite difference method (Newmark's method) are employed, for one-dimensional and two-dimensional cases. To validate the theoretical results, tables are shown comparing approximate and exact solutions. In addition, numerically the uniform decay rate for energy and the order of convergence of the approximate solution are also shown.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper studies the small-amplitude motion of a damped nonlinear elastic beam whose ends move in time. It constructs a Hermite finite-element semi-discretization in space together with Newmark time stepping, claims to derive stability and a priori error estimates establishing quadratic convergence in both space and time for the moving-boundary problem, and presents numerical experiments in one and two dimensions that compare computed solutions against exact solutions and illustrate energy decay.

Significance. If the error analysis correctly controls the additional transport terms induced by the time-dependent domain, the work would supply a useful rigorous justification for applying standard Hermite elements to variable-domain beam models. The explicit derivation of both semi-discrete and fully-discrete estimates together with numerical verification of the predicted rates would be a concrete contribution to the numerical analysis of moving-boundary hyperbolic problems.

major comments (2)
  1. [Abstract and numerical-method construction] The abstract asserts that error estimates yielding quadratic convergence are derived for the moving-boundary problem, yet the provided text supplies neither the precise weak formulation after application of the Reynolds transport theorem nor the function spaces in which the estimates are stated. Without these, it is impossible to verify whether the convective terms arising from d/dt ∫_{Ω(t)} are absorbed without degrading the O(h²) rate claimed for the Hermite interpolant.
  2. [Error estimates for semi-discrete scheme] Standard Hermite interpolation error on a fixed interval yields O(h²) in H¹ and O(h³) in L², but the moving ends introduce an additional consistency term in the error equation. The manuscript must show explicitly (presumably in the section deriving the error estimates) that this term is controlled by the same quadratic order; otherwise the fixed-domain theory invoked does not transfer directly.
minor comments (2)
  1. [Introduction] The abstract mentions both one- and two-dimensional simulations but does not clarify whether the two-dimensional case is a genuine 2-D plate or a collection of 1-D beams; this should be stated explicitly in the introduction.
  2. [Numerical results] Tables comparing approximate and exact solutions are referenced but their captions should include the precise norms (L², H¹, energy) and the mesh sizes or time steps used.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive suggestions. We address the two major comments below and will revise the manuscript to improve explicitness of the weak formulation and the control of moving-boundary consistency terms.

read point-by-point responses

  1. Referee: [Abstract and numerical-method construction] The abstract asserts that error estimates yielding quadratic convergence are derived for the moving-boundary problem, yet the provided text supplies neither the precise weak formulation after application of the Reynolds transport theorem nor the function spaces in which the estimates are stated. Without these, it is impossible to verify whether the convective terms arising from d/dt ∫_{Ω(t)} are absorbed without degrading the O(h²) rate claimed for the Hermite interpolant.

    Authors: The weak formulation obtained via the Reynolds transport theorem appears in Section 2, together with the natural energy space (H² ∩ H¹₀ on the instantaneous domain) in which the estimates are stated. The convective terms are absorbed into the stability estimate without order reduction, as they are controlled by the L² norm of the velocity and the bounded domain speed. To address the referee’s concern we will add an explicit statement of the weak form and the function spaces to the abstract and a short clarifying paragraph at the beginning of Section 3. revision: yes

  2. Referee: [Error estimates for semi-discrete scheme] Standard Hermite interpolation error on a fixed interval yields O(h²) in H¹ and O(h³) in L², but the moving ends introduce an additional consistency term in the error equation. The manuscript must show explicitly (presumably in the section deriving the error estimates) that this term is controlled by the same quadratic order; otherwise the fixed-domain theory invoked does not transfer directly.

    Authors: Section 4 derives the error equation for the semi-discrete scheme and isolates the consistency term generated by the moving boundary. Under the standing smoothness assumptions on the domain motion, this term is bounded by C h² (‖u‖_{H³} + ‖u_t‖_{H²}) via the standard Hermite interpolation estimates and the trace theorem on the moving ends; the resulting contribution remains O(h²) in the energy norm. We will insert a dedicated remark or short lemma immediately after the error equation to make this bounding step fully explicit. revision: yes

Circularity Check

0 steps flagged

No circularity; derivation of error estimates is self-contained

full rationale

The paper constructs semi-discrete and fully-discrete schemes with Hermite basis functions for the moving-boundary beam problem and derives stability plus error bounds directly, claiming quadratic convergence from those estimates. No quoted equations or steps reduce the claimed O(h²) rate to a fitted parameter, self-definition, or load-bearing self-citation chain. The abstract states the errors are derived and the order proven, without any reduction of the central result to its own inputs by construction. This is the normal case of an independent numerical analysis.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only the abstract is available; the ledger is therefore minimal and reflects the domain assumptions stated there.

axioms (1)
  • domain assumption Small-amplitude motion permits the use of the stated nonlinear beam model with internal damping.
    Explicitly referenced in the first sentence of the abstract.

pith-pipeline@v0.9.0 · 5658 in / 1179 out tokens · 38974 ms · 2026-05-25T09:45:31.853733+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

20 extracted references · 20 canonical work pages

  1. [1]

    Study of nonlinear vibration of Euler-Bernoulli beams by using analytical approximate techniques

    S. Bagheri, A. Nikkar e H. Ghaffarzadeh, “Study of nonlinear vibration of Euler-Bernoulli beams by using analytical approximate techniques”, Latin American Journal of Solids and Structures, 11 (2014) 157–168

    work page 2014

  2. [2]

    Spectral element methods for nonlinear spatio- temporal dynamics of an Euler-Bernoulli beam

    P. Z. Bar-Yoseph, D. Fisher e O. Gottlieb, “Spectral element methods for nonlinear spatio- temporal dynamics of an Euler-Bernoulli beam”, Computational Mechanics, 19 (1996) 136– 15

    work page 1996

  3. [3]

    Non-linear vibration of Euler-Bernoulli beams

    A. Barari, H. D. Kaliji, M.Ghadimi e G. Domairry, “Non-linear vibration of Euler-Bernoulli beams”, Latin American Journal of Solids and Structures , 8 (2011) [S.p.]

    work page 2011

  4. [4]

    Remark on the decay for damped string and beam equations

    P. Biler, “Remark on the decay for damped string and beam equations”,Nonlinear Analisys, 10 (1986) 839–842

    work page 1986

  5. [5]

    Remarks on the beam evolution equations in noncylindrical domains

    C. S. Q. Caldas, J. Limaco e R. K. Barreto, “Remarks on the beam evolution equations in noncylindrical domains”, Nonlinear Analysis, 70 (2009) 693–710

    work page 2009

  6. [6]

    Error estimates of finite element approximations for problems in linear elasticity - Part 3. Problems in elastodynamics discrete time approximations

    S.-I. Chou e C.-C.Wang, “Error estimates of finite element approximations for problems in linear elasticity - Part 3. Problems in elastodynamics discrete time approximations”, Archive for Rational Mechanics and Analysis , 73 (1980) 159–182

    work page 1980

  7. [7]

    The Finite Element Method for Elliptic Problems

    P.G. Ciarlet, “The Finite Element Method for Elliptic Problems”, North-Holland, Amster- dam, 1978

    work page 1978

  8. [8]

    Beam Equation with weak-internal damping in domain with moving boudary

    H. R. Clark, M. A. Rincon e R. D. Rodrigues, “Beam Equation with weak-internal damping in domain with moving boudary”, Applied Numerical Mathematics, 47 (2003) 139–157

    work page 2003

  9. [9]

    The initial value problem for a nonlinear semi-infinite string

    R. W. Dickey, “The initial value problem for a nonlinear semi-infinite string”, Proceedings of the Royal Society Edinburgh Section A , 82 (1978) 19–26

    work page 1978

  10. [10]

    A Finite Element Solution of the Beam Equation via Matlab

    S. R. Gunakata, D. M. G. Comissiong, K. Jordan e A. Sankar, “A Finite Element Solution of the Beam Equation via Matlab”, International Journal of Applied Science and Technology , 2 (2012) 80–88

    work page 2012

  11. [11]

    Dynamics of Transversely Vibrating Beams Using Four Engineering Theories

    S. M. Han, H. Benaroya e T. Wei, “Dynamics of Transversely Vibrating Beams Using Four Engineering Theories”, Journal of Sound and Vibration , 225 (1999) 935–988

    work page 1999

  12. [12]

    Nonlinear equations of motion for the elastic bending and torsion of twisted nonuniform rotor blades

    D. H. Hodgese e E. H. Dowell, “Nonlinear equations of motion for the elastic bending and torsion of twisted nonuniform rotor blades”, [S.l.]; National Aeronautics and Space Adiministration (NASA), 1974 (NASA TN D-7818)

    work page 1974

  13. [13]

    On the Nonlinear Deformation Geometry of Euler- Bernoulli

    D. H. Hodgese e R. A. Ormiston, “On the Nonlinear Deformation Geometry of Euler- Bernoulli”, [S.l.]; National Aeronautics and Space Adiministration (NASA) , 1960 (NASA A-7985)

    work page 1960

  14. [14]

    Vorlesungen ¨ uber Mechanik

    G. Kirchhoff, “Vorlesungen ¨ uber Mechanik”, Tauber, Leipzig, 1883

    work page

  15. [15]

    Vibrations of elastic strings: Mathematical aspects I and II

    L. A. Medeiros, J. L. Ferrel e S. B. Menezes, “Vibrations of elastic strings: Mathematical aspects I and II”, Journal of Computational Analysis and Applications , 4(2) (2002) 91–127; 4(3) (2002) 211-263

    work page 2002

  16. [16]

    Wave Propagation in Single-Walled Car- bon Nanotube Under Longitudinal Magnetic Field Using Nonlocal Euler-Bernoulli Beam Theory

    S. Narendar, S. S. Gupta e S. Gopalakrishnan, “Wave Propagation in Single-Walled Car- bon Nanotube Under Longitudinal Magnetic Field Using Nonlocal Euler-Bernoulli Beam Theory”, Applied Mathematical Modelling, 36 (2012) 4529–4538

    work page 2012

  17. [17]

    Nonlinear quasi-static finite element formulations for vis- coelastic Euler-Bernoulli and Timoshenko beams

    G. S. Payette e J. N. Reddy, “Nonlinear quasi-static finite element formulations for vis- coelastic Euler-Bernoulli and Timoshenko beams”, International Journal for Numerical Methods in Biomedical Engineering, 26 (2010) 1736–1755

    work page 2010

  18. [18]

    Existence, uniqueness and asymptotic behavior for solutions of the nonlinear bema equation

    D. C. Pereira, “Existence, uniqueness and asymptotic behavior for solutions of the nonlinear bema equation”, Nonlinear Analisys, 8 (1990) 613–623

    work page 1990

  19. [19]

    Semi-internal Stabilization for a Nonlinear Bernoulli-Euler Equation

    M. Tucsnak, “Semi-internal Stabilization for a Nonlinear Bernoulli-Euler Equation”, Math- ematical Methods in the Applied Sciences , 19 (1996) 897-907

    work page 1996

  20. [20]

    Isogeometric analysis of nonlinear Euler-Bernoulli beam vibrations

    O. Weeger, U. Wever and B. Simeon, “Isogeometric analysis of nonlinear Euler-Bernoulli beam vibrations”, Nonlinear Dynamics, 72 (2013) 813–835

    work page 2013