Revised identification of strain gradient elastic parameters

Abstract: The work reported in ``Granular micromechanics-based identification of isotropic strain gradient parameters for elastic geometrically nonlinear deformations" misidentified key terms in the grain-pair objective relative displacement when accounting for the second gradient of placement. In this paper, we correct that oversight by deriving a revised expression for the grain-pair objective relative displacement within the granular micromechanics framework. The amended terms, which resemble Christoffel symbols expressed in terms of strain gradients, modify the contributions of both the normal and tangential components to the strain energy and, consequently, alter the identified strain gradient elastic parameters. Importantly, the identification of the standard (first gradient) elastic tensor remains unchanged. This brief paper presents the corrected derivation, the resulting stiffness tensors for anisotropic strain gradient elasticity, and updated analytical expressions for the material parameters in both 2D and 3D isotropic settings.

• The paper corrects a key projection error in granular micromechanics, leading to revised expressions for normal and tangential strain energy contributions.
• It derives explicit closed-form formulas for first-gradient and second-gradient stiffness tensors in both 2D and 3D isotropic settings.
• The revised framework improves modeling accuracy for size-dependent phenomena and underpins advanced simulation techniques in higher-order elasticity.

Revised Identification of Strain Gradient Elastic Parameters

Overview and Motivation

The paper "Revised identification of strain gradient elastic parameters" (2506.14932) addresses a critical deficiency in the granular micromechanical framework proposed for the identification of strain gradient elastic coefficients, as originally formulated in prior work (notably, Barchiesi et al., ZAMM 2021). The accurate computation of higher-order elastic parameters is indispensable for modeling and engineering nano-structured materials and metamaterials that exhibit intrinsic length-scale effects. Most continuum approaches struggle with parameter identification for second-gradient elasticity due to ambiguities in representing the micromechanical behavior of discrete microstructures.

This study exposes and corrects a key misidentification in the prior granular micromechanics formalism: the projection of the grain-pair objective relative displacement that accounts for the second gradient of placement. The correction leads to revised expressions for the normal and tangential contributions to the strain energy, which in turn yield amended formulas and explicit closed-form expressions for the higher-order stiffness tensors in both 2D and 3D isotropic cases.

Core Mathematical Correction

The primary contribution is the corrected expression for the grain-pair objective relative displacement, $u^{np}$, which is foundational for the computation of internal energy in granular micromechanics-based second-gradient continua. The previous literature erroneously replaced the second-gradient term by a symmetrized gradient of the Green-St. Venant strain. This work shows that the correct term, structurally analogous to a generalized Christoffel symbol, is

$H_{ibc} = G_{ib,c} + G_{ic,b} - G_{bc,i}$

where $G_{ij}$ is the Green-St. Venant tensor and $F_{ij}$ the deformation gradient. The treatment rigorously derives the effect of this correction on both the normal and tangential displacement contributions.

Consequently, the density of the elastic energy is derived by integration over the orientation space (unit circle in 2D, unit sphere in 3D), properly scaling the contributions by symmetries and providing exhaustive index symmetrizations.

Explicit Stiffness Tensor Expressions

From the corrected energy density, new expressions are derived for the first-gradient and second-gradient stiffness tensors ($\mathbb{C}$, $\mathbb{M}$, and $\mathbb{D}$). For isotropic settings, the paper presents closed-form analytic expressions for all nonzero components as functions of the normal and tangential integrated stiffnesses ($\bar{k}_\eta$, $\bar{k}_\tau$), the distance $L$ between grain-pairs, and—through explicit formulas—standard elastic parameters ($Y_{2D/3D}$, $\nu_{2D/3D}$):

• In both 2D and 3D cases, the first-gradient tensor identifications ($\mathbb{C}$) remain unchanged, preserving correspondence with classical elasticity (Lamé constants, Young’s modulus, Poisson ratio).
• For the second-gradient tensor ($\mathbb{D}$), the amendments lead to modified structure and coefficients. The explicit identification for the independent components ensures the proper physical interpretation and enables immediate use in higher-order continuum formulations.

In 2D, six independent parameters $d_1$ through $d_6$ are given as explicit combinations of $\bar{k}_\eta$, $\bar{k}_\tau$. In 3D, the formulation admits seven groups, with five independent isotropic constants ($c_3$ through $c_7$), directly related to Mindlin’s coefficients, and all expressed as rational functions of the underlying elastic moduli and geometric parameters.

Implications and Theoretical Significance

The rigorous correction of the micromechanical-to-continuum translation, particularly for higher-gradient elasticity, is essential for both the interpretation of experiments (e.g., identification from nano-indentation, acoustic wave dispersion) and the reliable prediction of size-dependent mechanical phenomena. Notably:

• The corrected framework accurately describes the emergence and influence of multiple characteristic lengths associated with different deformation modes, which are central for capturing size effects and boundary layer phenomena in nano- and metamaterials.
• The revised parameter identification reinforces the link between discrete microstructural mechanics and continuum second-gradient theories, supporting ongoing efforts in multiscale modeling, homogenization, and materials-by-design pipelines.
• The separation and explicitness in 2D versus 3D cases facilitate direct application in computational schemes (e.g., C1 FE formulations, meshless methods) and in the analytical treatment of boundary value problems involving higher-order elasticity.

The revisited structure clarifies that only the second-gradient (and not the first-gradient) elastic parameters are altered by the correction, which anchors the physical meaningfulness of prior standard elasticity results even as higher-order theory is robustly completed.

Strong Results and Claims

• The identification of higher-order (second-gradient) stiffness tensors is fundamentally altered; prior analytic and numerical results based on the misidentified term must be revisited with the amended formulas.
• Analytical expressions for all nonzero components of second-gradient stiffness tensors in isotropic 2D and 3D settings are given explicitly, parametrized by both microstructural (grain-pair) and macroscopic elastic quantities.
• The framework generalizes directly to anisotropic media with appropriate orientation-dependent stiffnesses, supporting advanced applications in architected metamaterials with nontrivial symmetry.

Future Developments

These corrected identifications open pathways for:

• More reliable parameter extraction procedures from atomistic and discrete simulations, and from experimental mechanics at sub-micrometer scales.
• The robust modeling of materials and metamaterials whose response is governed by size effects and higher-order interactions, especially for the development of programmable and adaptive microstructures.
• Extensions to coupled multi-physics generalizations (e.g., flexoelectricity, gradient thermoelasticity) where the structure of higher-order moduli critically depends on the accurate micromechanical derivation.

Subsequent work can focus on systematic recalibration of previously published numerical results, experimental fits, and optimization studies using the correct coefficients.

Conclusion

This work establishes a corrected and internally consistent micromechanics-based method for identifying strain gradient elastic parameters, thereby resolving a critical foundation for higher-order continuum elasticity modeling. The explicit and rigorously derived formulas are immediately relevant for researchers working on the analytical, computational, and experimental fronts of gradient elasticity, ensuring the fidelity and predictive capacity of second-gradient theories in capturing real material behavior at small scales.