pith. sign in

arxiv: 2512.21483 · v2 · submitted 2025-12-25 · ⚛️ physics.class-ph

Asymptotically exact dimension reduction of functionally graded anisotropic rods

Khanh Chau Le This is my paper

Pith reviewed 2026-05-16 20:07 UTC · model grok-4.3

classification ⚛️ physics.class-ph
keywords dimension reductionvariational-asymptotic methodfunctionally graded rodsanisotropic elasticityasymptotic exactnessPrager-Synge identitylong-wave approximationone-dimensional model
0
0 comments X

The pith

The variational-asymptotic method produces an asymptotically exact one-dimensional model for functionally graded anisotropic rods directly from three-dimensional elasticity.

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

This paper develops a one-dimensional theory for slender rods whose material properties vary along their length and exhibit general anisotropy by applying the variational-asymptotic method to the full three-dimensional elastic equations. Dual cross-sectional problems are solved numerically to obtain rigorous upper and lower bounds on the average transverse energy density, after which the Prager-Synge identity supplies an error estimate in the energetic norm that proves the reduction is asymptotically exact. The same error bound is shown to hold for low-frequency vibrations, and the resulting dispersion relations agree with exact three-dimensional solutions for waves in composite rods. A sympathetic reader would care because conventional rod theories can produce deflection errors as large as 20 percent, while this framework reduces the discrepancy below 3 percent and demonstrates O(h/L) convergence in numerical benchmarks.

Core claim

The developed one-dimensional model captures the long-wave asymptotic behavior of the three-dimensional elastic body with high fidelity. Numerical benchmarks indicate that while the naive rod theory incurs errors up to 20 percent in deflection predictions, the current VAM framework reduces this discrepancy to below 3 percent, with log-log convergence studies confirming the theoretical O(h/L) accuracy. The estimate is extended to the dynamic regime for low-frequency vibrations, where one-dimensional dispersion relations match exact analytical three-dimensional solutions for wave propagation in composite rods.

What carries the argument

Numerical solution of dual cross-sectional problems, which furnish rigorous upper and lower bounds on the average transverse energy density and are combined with the Prager-Synge identity to produce the asymptotic error estimate.

If this is right

  • Rigorous upper and lower bounds are obtained for the average transverse energy density through the dual cross-sectional problems.
  • The error estimate in the energetic norm extends directly to low-frequency dynamic problems.
  • Dispersion relations of the one-dimensional model coincide with exact three-dimensional solutions in the long-wave limit for composite rods.
  • Deflection predictions improve from errors of 20 percent in naive rod theory to below 3 percent with confirmed O(h/L) convergence.

Where Pith is reading between the lines

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

  • The same dual-problem technique could be adapted to derive reduced models for functionally graded plates or shells under similar slenderness assumptions.
  • Engineering software could incorporate the numerical cross-sectional solver to obtain accurate effective stiffnesses for arbitrary anisotropic grading without ad-hoc assumptions.
  • The approach supplies a systematic route to test whether a given material variation is slow enough for the one-dimensional description to remain accurate.

Load-bearing premise

The rod must be slender, with small thickness-to-length ratio, and the material properties must vary slowly enough along the length that the cross-sectional problems remain well-posed.

What would settle it

A sequence of three-dimensional finite-element computations on rods with progressively smaller h/L ratios in which the deflection error of the one-dimensional model fails to decrease proportionally to h/L.

Figures

Figures reproduced from arXiv: 2512.21483 by Khanh Chau Le.

Figure 1
Figure 1. Figure 1: A rod. The coordinates x α span the connected 2D domain S, such that the first moment of area vanishes, R S x α da = 0. For an anisotropic FG-rod, the dynamic behavior is governed by the ac￾tion functional expressed in curvilinear coordinates. Following the general formulation for slender rods [7, 20], the action functional reads: Z t1 t0 Z V L(x α , w˙ , ∇w) dv dt = Z t1 t0 Z L 0 Z S [T(x α , w˙ ) − W(x α… view at source ↗
Figure 2
Figure 2. Figure 2: Normalized torsional stiffness gˆ as a function of the aspect ratio a: (i) at fixed δ = 1, γU = 1 and four different γL = 0.2, 0.4, 0.6, 0.8 (left), and (ii) at fixed γL = 0.1, γU = 0.4 and four different δ = 1, 2, 3, 4 (right). As shown in [PITH_FULL_IMAGE:figures/full_fig_p016_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Numerical solution of the stress function [PITH_FULL_IMAGE:figures/full_fig_p017_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Percentage contribution of the transverse bending stiffness [PITH_FULL_IMAGE:figures/full_fig_p018_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Percentage contribution of the transverse cross stiffness [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Contour distributions of the plane strain stress functions [PITH_FULL_IMAGE:figures/full_fig_p019_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Normalized torsional stiffness gˆ as a function of the aspect ratio a for FG-rod Barium Titanate/Corundum with δ = 4. significant deviation of the local 3D stress state from the simplified predic￾tions of classical rod theories. These results underscore the necessity of the variational-asymptotic approach for accurately capturing the energetic land￾scape of low-symmetry functionally graded materials [PITH… view at source ↗
Figure 8
Figure 8. Figure 8: Contour plots of the coupled stress functions [PITH_FULL_IMAGE:figures/full_fig_p022_8.png] view at source ↗
read the original abstract

This study utilizes the variational-asymptotic method to establish a one-dimensional theory for functionally graded rods characterized by general anisotropy from the three-dimensional elasticity theory. A distinctive feature of this dimension reduction procedure is the numerical solution of dual cross-sectional problems, which provide rigorous upper and lower bounds for the average transverse energy density. By employing the Prager-Synge identity, we derive an error estimate in the energetic norm to establish the asymptotic exactness of the model. This estimate is extended to the dynamic regime for low-frequency vibrations. Furthermore, the dynamic validity of the theory is confirmed by comparing the one-dimensional dispersion relations with exact analytical three-dimensional solutions for wave propagation in composite rods. The results show that the developed one-dimensional model captures the long-wave asymptotic behavior of the three-dimensional elastic body with high fidelity. Numerical benchmarks indicate that while the naive rod theory incurs errors up to $20\%$ in deflection predictions, the current VAM framework reduces this discrepancy to below $3\%$, with log-log convergence studies confirming the theoretical $O(h/L)$ accuracy.

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 manuscript develops a one-dimensional theory for functionally graded anisotropic rods using the variational-asymptotic method starting from three-dimensional elasticity. Key features include numerical solution of dual cross-sectional problems to obtain rigorous bounds on the average transverse energy density, application of the Prager-Synge identity to derive an error estimate establishing asymptotic exactness, extension to dynamic low-frequency vibrations, and validation through comparison with exact three-dimensional dispersion relations for composite rods. Numerical results show reduction of deflection errors from up to 20% to below 3% with O(h/L) convergence.

Significance. If the numerical aspects are solidified, the work offers a valuable contribution to dimension reduction techniques by providing asymptotically exact models with error bounds for generally anisotropic and functionally graded materials, which are challenging for analytical approaches. The use of dual problems for bounds and Prager-Synge for error control is a strength, and the dynamic validation adds credibility.

major comments (2)
  1. [Numerical benchmarks] Numerical benchmarks section: the log-log plots and global deflection error reductions (from 20% to <3%) are reported for the slenderness parameter h/L, but no mesh-refinement studies, residual estimates, or convergence data are provided for the dual cross-sectional problems. Since the Prager-Synge identity in the error estimate relies on these numerical bounds approaching the exact infima, the absence of such verification weakens the claim of rigorous asymptotic exactness for arbitrary anisotropy and functional grading.
  2. [Dynamic extension] Dynamic regime extension: the manuscript states that the error estimate is extended to low-frequency vibrations, but does not detail how the Prager-Synge identity is modified for the time-dependent case or whether the same static dual-problem bounds suffice without additional dynamic cross-sectional corrections.
minor comments (2)
  1. [Abstract] Abstract: the phrase 'log-log convergence studies confirming the theoretical O(h/L) accuracy' should specify whether the studies address only the rod slenderness or also the discretization of the cross-section dual problems.
  2. [Notation and definitions] Notation: ensure consistent use of symbols for the average transverse energy density and the dual energies throughout; a table summarizing the bound quantities would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and positive overall assessment. We address each major comment point by point below and will incorporate the suggested clarifications and additional data into the revised manuscript.

read point-by-point responses

  1. Referee: Numerical benchmarks section: the log-log plots and global deflection error reductions (from 20% to <3%) are reported for the slenderness parameter h/L, but no mesh-refinement studies, residual estimates, or convergence data are provided for the dual cross-sectional problems. Since the Prager-Synge identity in the error estimate relies on these numerical bounds approaching the exact infima, the absence of such verification weakens the claim of rigorous asymptotic exactness for arbitrary anisotropy and functional grading.

    Authors: We agree that explicit verification of the numerical accuracy of the dual cross-sectional problems is necessary to fully support the rigor of the error bounds. In the revised manuscript we will add mesh-refinement studies, residual estimates, and convergence tables for representative anisotropic and graded cross-sections, confirming that the computed upper and lower bounds approach the exact infima to within the tolerance required by the Prager-Synge estimate. revision: yes

  2. Referee: Dynamic regime extension: the manuscript states that the error estimate is extended to low-frequency vibrations, but does not detail how the Prager-Synge identity is modified for the time-dependent case or whether the same static dual-problem bounds suffice without additional dynamic cross-sectional corrections.

    Authors: For low-frequency vibrations the inertial contributions enter as higher-order perturbations in the variational-asymptotic expansion. Consequently the leading-order transverse energy density remains governed by the same static dual problems, and the Prager-Synge identity carries over directly with only an additional O((h/L)^2) remainder that is already absorbed in the existing error estimate. We will insert a short derivation of this reduction in the revised text to make the extension explicit. revision: yes

Circularity Check

0 steps flagged

Derivation proceeds from 3D elasticity via independent numerical cross-section solves; no reduction to inputs by construction

full rationale

The central chain begins with the three-dimensional linear elasticity equations, applies the variational-asymptotic method to obtain a one-dimensional model, and solves dual cross-sectional problems numerically to produce rigorous upper and lower bounds on average transverse energy density. The Prager-Synge identity is then invoked to convert those bounds into an a-posteriori error estimate in the energetic norm. Dispersion relations are validated against exact analytical three-dimensional solutions for composite rods. None of these steps equates a derived quantity to a fitted parameter or to a self-citation that itself assumes the target result. The numerical convergence of the dual problems is a separate question of discretization accuracy and does not render the derivation circular.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The reduction rests on standard assumptions of slender-body asymptotics and the existence of well-posed cross-sectional problems; no new free parameters or invented entities are introduced beyond the numerical solution procedure.

axioms (2)
  • domain assumption The rod is slender with slow axial variations in material properties
    Required for the separation of scales in the variational-asymptotic procedure
  • domain assumption The dual cross-sectional problems admit unique solutions that furnish rigorous energy bounds
    Invoked to justify the Prager-Synge error estimate

pith-pipeline@v0.9.0 · 5472 in / 1292 out tokens · 19345 ms · 2026-05-16T20:07:42.004876+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

43 extracted references · 43 canonical work pages

  1. [1]

    Propagation of harmonic waves in composite circular cylindrical shells

    Armenàkas, A.E., 1967. Propagation of harmonic waves in composite circular cylindrical shells. I. Theoretical investigation. AIAA Journal 5(4), 740–744

    work page 1967

  2. [2]

    Free vibration analysis of functionally graded beams with simply supported edges

    Aydogdu, M., Taskin, V., 2007. Free vibration analysis of functionally graded beams with simply supported edges. Materials & Design 28(5), 1651–1656

    work page 2007

  3. [3]

    On the proof of the Saint-Venant principle for bodies of arbitrary shape

    Berdichevsky, V.L., 1974. On the proof of the Saint-Venant principle for bodies of arbitrary shape. Journal of Applied Mathematics and Mechan- ics 38(5), 799–813

    work page 1974

  4. [4]

    Variational-asymptotic method of construct- ing a theory of shells

    Berdichevsky, V.L., 1979. Variational-asymptotic method of construct- ing a theory of shells. Journal of Applied Mathematics and Mechanics 43(4), 711–736

    work page 1979

  5. [5]

    High frequency, long wave shell vibration

    Berdichevsky, V.L., Le, K.C., 1980. High frequency, long wave shell vibration. Journal of Applied Mathematics and Mechanics 44(3), 520– 525

    work page 1980

  6. [6]

    On the energy of an elastic rod

    Berdichevsky, V.L., 1981. On the energy of an elastic rod. Journal of Applied Mathematics and Mechanics 45(4), 518–529. 34

    work page 1981

  7. [7]

    Variational Principles of Continuum Mechan- ics, Vol

    Berdichevsky, V.L., 2009. Variational Principles of Continuum Mechan- ics, Vol. 2. Springer Verlag, Berlin

    work page 2009

  8. [8]

    A new simple shear and normal deformations theory for functionally graded beams

    Bourada, M., Kaci, A., Houari, M.S.A., Tounsi, A., 2015. A new simple shear and normal deformations theory for functionally graded beams. Steel and Composite Structures 18(2), 409–423

    work page 2015

  9. [9]

    International Journal of Mechanical Sciences 45(3), 519–539

    Chakraborty, A., Gopalakrishnan, S., Reddy, J., 2003.Anewbeamfinite element for the analysis of functionally graded materials. International Journal of Mechanical Sciences 45(3), 519–539

    work page 2003

  10. [10]

    Two-dimensional elasticity solution of elastic strips and beams made of functionally graded mate- rials under tension and bending

    Chu, P., Li, X.F., Wu, J.X., Lee, K., 2015. Two-dimensional elasticity solution of elastic strips and beams made of functionally graded mate- rials under tension and bending. Acta Mechanica 226, 2235–2253

    work page 2015

  11. [11]

    Concepts and Applications of Finite Element Analy- sis

    Cook, R.D., 2007. Concepts and Applications of Finite Element Analy- sis. John Wiley & Sons

    work page 2007

  12. [12]

    Decaying states of plane strain in a semi-infinite strip and boundary conditions for plate theory

    Gregory, R.D., Wan, F.Y., 1984. Decaying states of plane strain in a semi-infinite strip and boundary conditions for plate theory. Journal of Elasticity 14(1), 27–64

    work page 1984

  13. [13]

    The pressurized hollow cylinder or disk problem for functionally graded isotropic linearly elastic materials

    Horgan, C.O., Chan, A.M., 1999. The pressurized hollow cylinder or disk problem for functionally graded isotropic linearly elastic materials. Journal of Elasticity 55, 43–59

    work page 1999

  14. [14]

    Analytical solution for functionally graded anisotropic cantilever beam subjected to linearly distributed load

    Huang, D.J., Ding, H.J., Chen, W.Q., 2007. Analytical solution for functionally graded anisotropic cantilever beam subjected to linearly distributed load. Applied Mathematics and Mechanics 28(7), 855–860

    work page 2007

  15. [15]

    Three-dimensional elasticity solution for bend- ing of functionally graded rectangular plates

    Kashtalyan, M., 2004. Three-dimensional elasticity solution for bend- ing of functionally graded rectangular plates. European Journal of Mechanics-A/Solids 23(5), 853–864

    work page 2004

  16. [16]

    Recent progress of functionally gradient materials in Japan

    Koizumi, M., 1992. Recent progress of functionally gradient materials in Japan. In: 16th Annual Conference on Composites and Advanced Ceramic Materials I. Vol. 13, p. 333

    work page 1992

  17. [17]

    New concepts for linear beam theory with arbitrary geometry and loading

    Ladevèze, P., Simmonds, J., 1998. New concepts for linear beam theory with arbitrary geometry and loading. European Journal of Mechanics- A/Solids 17(3), 377–402. 35

    work page 1998

  18. [18]

    The theory of piezoelectric shells

    Le, K.C., 1986. The theory of piezoelectric shells. Journal of Applied Mathematics and Mechanics 50(1), 98–105

    work page 1986

  19. [19]

    High frequency vibrations and wave propagation in elastic shells: variational-asymptotic approach

    Le, K.C., 1997. High frequency vibrations and wave propagation in elastic shells: variational-asymptotic approach. International Journal of Solids and Structures 34(30), 3923–3939

    work page 1997

  20. [20]

    Vibrations of Shells and Rods

    Le, K.C., 1999. Vibrations of Shells and Rods. Springer Verlag, Berlin

    work page 1999

  21. [21]

    An asymptotically exact theory of smart sand- wich shells

    Le, K.C., Yi, J.H., 2016. An asymptotically exact theory of smart sand- wich shells. International Journal of Engineering Science 106, 179–198

    work page 2016

  22. [22]

    An asymptotically exact theory of functionally graded piezoelectric shells

    Le, K.C., 2017. An asymptotically exact theory of functionally graded piezoelectric shells. International Journal of Engineering Science 112, 42–62

    work page 2017

  23. [23]

    Asymptotically exact theory of functionally graded elastic beams

    Le, K.C., Tran, T.M., 2025. Asymptotically exact theory of functionally graded elastic beams. International Journal of Engineering Science 209, 104214

    work page 2025

  24. [24]

    Dynamic variational-asymptotic proce- dure for laminated composite shells—part II: High-frequency vibration analysis

    Lee, C.Y., Hodges, D.H., 2009. Dynamic variational-asymptotic proce- dure for laminated composite shells—part II: High-frequency vibration analysis. Journal of Applied Mechanics 76(1), 011003

    work page 2009

  25. [25]

    A unified approach for analyzing static and dynamic be- haviors of functionally graded Timoshenko and Euler–Bernoulli beams

    Li, X.F., 2008. A unified approach for analyzing static and dynamic be- haviors of functionally graded Timoshenko and Euler–Bernoulli beams. Journal of Sound and Vibration 318(4-5), 1210–1229

    work page 2008

  26. [26]

    Partial Differential Equation Toolbox User’s Guide

    Mathworks, 1995. Partial Differential Equation Toolbox User’s Guide. The MathWorks, Inc., Natick, MA

    work page 1995

  27. [27]

    Exact solutions for functionally graded and laminated elastic materials

    Mian, M.A., Spencer, A., 1998. Exact solutions for functionally graded and laminated elastic materials. Journal of the Mechanics and Physics of Solids 46(12), 2283–2295

    work page 1998

  28. [28]

    High frequency vibrations of crystal plates

    Mindlin, R.D., 1961. High frequency vibrations of crystal plates. Quar- terly of Applied Mathematics 19(1), 51–61

    work page 1961

  29. [29]

    General study of axisymmetric waves in layered anisotropic fibers and their composites

    Nayfeh, A.H., Nagy, P.B., 1996. General study of axisymmetric waves in layered anisotropic fibers and their composites. The Journal of the Acoustical Society of America 99(2), 931–941. 36

    work page 1996

  30. [30]

    Recent development status of functionally gradient materials

    Niino, M., Maeda, S., 1990. Recent development status of functionally gradient materials. ISIJ International 30(9), 699–703

    work page 1990

  31. [31]

    Exact solution for functionally graded anisotropic elastic compositelaminates.JournalofCompositeMaterials37(21), 1903–1920

    Pan, E., 2003. Exact solution for functionally graded anisotropic elastic compositelaminates.JournalofCompositeMaterials37(21), 1903–1920

    work page 2003

  32. [32]

    Approximations in elasticity based on the concept of function space

    Prager, W., Synge, J.L., 1947. Approximations in elasticity based on the concept of function space. Quarterly of Applied Mathematics 5(3), 241–269

    work page 1947

  33. [33]

    Axisymmetric bending of functionally graded circular and annular plates

    Reddy, J.N., Wang, C.M., Kitipornchai, S., 1999. Axisymmetric bending of functionally graded circular and annular plates. European Journal of Mechanics-A/Solids 18(2), 185–199

    work page 1999

  34. [34]

    Theories and analysis of functionally graded beams

    Reddy, J.N., Ruocco, E., Loya, J.A., Neves, A.M., 2021. Theories and analysis of functionally graded beams. Applied Sciences 11(15), 7159

    work page 2021

  35. [35]

    30 Years of functionally graded materials: An overview of manufacturing methods, applications and future challenges

    Saleh, B., Jiang, J., Fathi, R., Al-Hababi, T., Xu, Q., Wang, L., Song, D., Ma, A., 2020. 30 Years of functionally graded materials: An overview of manufacturing methods, applications and future challenges. Compos- ites Part B: Engineering 201, 108376

    work page 2020

  36. [36]

    An elasticity solution for functionally graded beams

    Sankar, B.V., 2001. An elasticity solution for functionally graded beams. Composites Science and Technology 61(5), 689–696

    work page 2001

  37. [37]

    Free and forced vibration of a function- ally graded beam subjected to a concentrated moving harmonic load

    Şimşek, M., Kocatürk, T., 2009. Free and forced vibration of a function- ally graded beam subjected to a concentrated moving harmonic load. Composite Structures 90(4), 465–473

    work page 2009

  38. [38]

    Bending and free vibration of functionally graded beams using various higher-order shear deformation beam theo- ries

    Thai, H.T., Vo, T.P., 2012. Bending and free vibration of functionally graded beams using various higher-order shear deformation beam theo- ries. International Journal of Mechanical Sciences 62(1), 57–66

    work page 2012

  39. [39]

    Saint-Venant’s principle

    Toupin, R.A., 1965. Saint-Venant’s principle. Archive for Rational Me- chanics and Analysis 18(2), 83–96

    work page 1965

  40. [40]

    Asymptotic theory for static behavior of elastic anisotropic I-beams

    Volovoi, V.V., Hodges, D.H., Berdichevsky, V.L., Sutyrin, V.G., 1999. Asymptotic theory for static behavior of elastic anisotropic I-beams. International Journal of Solids and Structures 36(7), 1017–1043. 37

    work page 1999

  41. [41]

    On Timoshenko- like modeling of initially curved and twisted composite beams

    Yu, W., Hodges, D.H., Volovoi, V., Cesnik, C.E., 2002. On Timoshenko- like modeling of initially curved and twisted composite beams. Interna- tional Journal of Solids and Structures 39(19), 5101–5121

    work page 2002

  42. [42]

    Elasticity solutions versus asymptotic sec- tional analysis of homogeneous, isotropic, prismatic beams

    Yu, W., Hodges, D.H., 2004. Elasticity solutions versus asymptotic sec- tional analysis of homogeneous, isotropic, prismatic beams. Journal of Applied Mechanics 71(1), 15–23

    work page 2004

  43. [43]

    Generalized Timoshenko theory of the variational asymptotic beam sectional analysis

    Yu, W., Hodges, D.H., 2005. Generalized Timoshenko theory of the variational asymptotic beam sectional analysis. Journal of the American Helicopter Society 50(1), 46–55. 38

    work page 2005