A non

Scientific Reports volume 13, Article number: 6562 (2023) Cite this article

312 Accesses

1 Citations

Metrics details

In the past, to model fibre stiffness of finite-radius fibres, previous finite-strain (nonlinear) models were mainly based on the theory of non-linear strain-gradient (second-gradient) theory or Kirchhoff rod theory. We note that these models characterize the mechanical behaviour of polar transversely isotropic solids with infinitely many purely flexible fibres with zero radius. To introduce the effect of fibre bending stiffness on purely flexible fibres with zero radius, these models assumed the existence of couple stresses (contact torques) and non-symmetric Cauchy stresses. However, these stresses are not present on deformations of actual non-polar elastic solids reinforced by finite-radius fibres. In addition to this, the implementation of boundary conditions for second gradient models is not straightforward and discussion on the effectiveness of strain gradient elasticity models to mechanically describe continuum solids is still ongoing. In this paper, we develop a constitutive equation for a non-linear non-polar elastic solid, reinforced by embedded fibers, in which elastic resistance of the fibers to bending is modelled via the classical branches of continuum mechanics, where the development of the theory of stresses is based on non-polar materials; that is, without using the second gradient theory, which is associated with couple stresses and non-symmetric Cauchy stresses. In view of this, the proposed model is simple and somewhat more realistic compared to previous second gradient models.

Fibre-reinforced composite materials have often been used in recent engineering applications. The rapid growth in manufacturing industries has led to the need for the improvement of materials in terms of strength, stiffness, density, and lower cost with improved sustainability. Fibre-reinforced composite materials have emerged as one of the materials possessing such improvement in properties serving their potential in a variety of applications1,2,3,4. The infusion of natural synthetic or natural fibers in the fabrication of composite materials has revealed significant applications in a variety of fields such as biomedical, automobile, mechanical, construction, marine and aerospace5,6,7,8. In biomechanics, some soft tissues can be modelled as fibre-reinforced composite materials9,10. In modern heavy engineering, the heavy traditional materials are gradually being replaced by fibre-reinforced polymer composite structures of lower weight and higher strength. These structures, such as railroads and bridges, are always under the action of dynamic moving loads caused by the moving vehicular traffic. Hence, in view of the above, a rigourous construction of a mechanical constitutive model, based on the sound theory of continuum mechanics, for non-polar fibre-reinforced solids, is paramount, and is of valuable interest in engineering designs and would find many practical applications.

The long history11,12,13 of mechanics of non-polar fiber-reinforced solids has, in general, significantly enriched and advanced the knowledge of solid mechanics. A boundary value problem for a non-polar elastic solid reinforced by (finite radius) fibres can be solved using the Finite Element Method (FEM), if small elements are permittable to mesh the fibres. If we treat the fibres to be an isotropic solid but have a different material properties from the matrix (material that is not attributable to the fibers) properties, we can use an inhomogeneous strain energy function

in solving the FEM problem, where \(\lambda _1,\lambda _2\) and \(\lambda _3\) are the pricipal stretches. We note that, due to the finite radius of the fibres, bending resistance due to changes in the curvature for the fibres, is observed. However, if the fibre radius is significantly small, meshing the fibres and the matrix can be troublesome and hence it may not be possible to seek a boundary value solution via the FEM. To overcome this significantly small radius problem, a FEM solution can be obtained using a transversely elastic strain energy function13

where \({\varvec{U}}\) is the right-stretch tensor and \({\varvec{a}}\) is the unit preferred vector in the reference configuration. We note that this transversely isotropic model contains infinitely many purely flexible fibres with zero radius; hence this model cannot model elastic resistance due to changes in the curvature for the fibres. We emphasize that the Cauchy stress in both isotropic and transversely isotropic models is symmetric and this is actually observed in a non-polar solid in the absence of a couple stress. To model the effect of elastic resistance due to changes in the curvature for the fibres, recent models14,15,16,17 that are framed in the setting of the non-linear strain-gradient theory or Kirchhoff rod theory18, were developed. We note that these second-gradient models characterize the mechanical behaviour of (polar) transversely isotropic solids with infinitely many purely flexible fibres with zero radius. But, in order to simulate the effect of fibre bending stiffness on purely flexible fibres with zero radius, the second-gradient models introduce the existence of a couple stress and a non-symmetric Cauchy stress in the constitutive equations; we must emphasize that both of these stresses are not present on deformations of actual non-polar elastic solids reinforced by finite-radius fibres. In general, higher gradient elasticity models are used to describe mechanical structures at the micro-and nano-scale or to regularize certain ill-posed problems by means of these higher gradient contributions. Discussion on the effectiveness of higher gradient elasticity models to mechanically describe continuum solids is still ongoing19,20,21.

Hence, the objective of this paper is to propose an approximate model to simulate the mechanical behaviour of actual non-polar elastic solids reinforced by finite- radius fibres, where the Cauchy stress is symmetric and fibre bending resistance is caused by changes in curvature of the fibres. We focus on changes in fibre curvature, since in composite solids, these changes play an important role in the mechanical behaviour of solids. Since our model contains infinitely many fibres with zero radius, we exclude the effects due to fibre 'twist’. In fact Spencer and Soldatos17 stated that

"In doing this, we exclude effects due to fibre 'splay’ and fibre 'twist’, both of which feature in liquid crystal theory, but it is plausible that in fibre composite solids the major factor is fibre curvature."

Our proposed model does not require couple stresses (which are not observed in actual non-polar elastic solids reinforced by finite-radius fibres) to describe elastic resistance of the fibers to bending.

Spectral approach14,22 is used in the modelling and this is preliminary described in Sects. "Preliminaries" and "Strain energy function", where in Sect. "Strain energy function" a strain energy function contains a vector that governs the changes in the fibre curvature. A prototype of the strain energy is given in Sect. "Strain energy prototype" and boundary value problems to study the effect of fibre bending resistance are presented in Sect. "Boundary value problem".

Deformation due to application of boundary displacement and boundary traction. \(B_r\) is the reference (undeformed) configuration, \(B_t\) is the current configuration, \({\varvec{x}}\) and \({\varvec{y}}\) are, respectively, the position vectors of X in the reference and current configurations, where X represents a generic particle of the solid body.

In this communication, all subscripts i, j and k take the values 1,2,3, unless stated otherwise. In terms of spectral invariants, the deformation gradient \({\varvec{F}}\) is described via

where \({\varvec{y}}\) and \({\varvec{x}}\) denote, respectively, the position vectors of a solid body particle in the current and reference configurations (see Fig. 1 ) ; \(\lambda _i\) is a principal stretch, \({\varvec{v}}_i\) is an eigenvector of the left stretch tensor \({\varvec{V}}= {\varvec{F}}(\lambda _i,{\varvec{v}}_i,{\varvec{v}}_i)\) and \({\varvec{u}}_i\) is an eigenvector of the right- stretch tensor \({\varvec{U}}= {\varvec{F}}(\lambda _i,{\varvec{u}}_i,{\varvec{u}}_i)\). Note that the right Cauchy-Green tensor \({\varvec{C}}= {\varvec{F}}(\lambda _i^2,{\varvec{u}}_i,{\varvec{u}}_i)\) and the rotation tensor \({\varvec{R}}= {\varvec{F}}(\lambda _i=1,{\varvec{v}}_i,{\varvec{u}}_i)\), where \({\varvec{F}}={\varvec{R}}{\varvec{U}}\). We only consider incompressible elastic solids, where \(\det {\varvec{F}}=1\), \(\det\) indicates the determinant of a tensor and the effect of body forces is assumed negligible. The summation convention is not used here.

To model the kinematics of the embedded fibers, we assume the body, regarded as a homogenized continuum consisting of matrix material and fibers together. We model this material by considering a transversely elastic solid with the preferred unit directions \({\varvec{a}}({\varvec{x}})\) in the reference configuration and these preferred directions becomes the vector

in the current configuration, where \({\varvec{f}}\) is a unit vector. In our proposed model, the directional derivative of the fibre unit vector in the fibre direction, i.e.,

plays an important role in modelling elastic resistance due to changes in curvature for the fibres. In view of this we endow a vector \({\varvec{d}}\) associated with \({\varvec{c}}\) (we will make the association clear later) in (5), which is independent of \({\varvec{F}}\), i.e.14,15

where

\({\bar{{\varvec{F}}}}({\varvec{x}})\) is the deformation tensor independent of \({\varvec{F}}\), i.e., \({\varvec{d}}\) is not embedded in the matrix, and so in general its image \({\bar{{\varvec{F}}}}^{-T}{\varvec{d}}\) in the current configuration is not directly connected to the deformation of the matrix. Clearly from (6), we have \({\varvec{d}}\cdot {\varvec{a}}=0\) (See Fig. 2 for geometric interpretation). If we let \({\bar{{\varvec{F}}}}={\varvec{F}}\), we then have the association \({\varvec{c}}= {\varvec{F}}^{-T}{\varvec{d}}\)14,15. To facilitate the process of modelling, we express the vector

where \({\varvec{k}}\) is a unit vector with the property \({\varvec{a}}\cdot {\varvec{k}}=0\). To model elastic resistance due to changes in curvature of the fibres, we assume the objective strain energy

for every rotation tensor \({\varvec{Q}}\). Following, the work of Shariff22,23, W can be characterised by the spectral invariants

and the scaler \(\rho\), where \(\lambda _i\) and \({\varvec{u}}_i\) are, repectively, the eigenvalues and eigenvectors of \({\varvec{U}}\). Hence, we can express

taking note the \({W}_{(a)}\) must satisfy the P-property described in24 associated with the coalescence of principal stretches \(\lambda _i\). In view that W should be independent of the sign of \({\varvec{a}}\) and \({\varvec{d}}\), we express

Geometric meaning of directional derivative vectors: \({\varvec{F}}\ne {\bar{{\varvec{F}}}}\), \({\varvec{a}}\cdot {\varvec{d}}={\varvec{b}}\cdot {\varvec{c}}= {\bar{{\varvec{b}}}}\cdot {\bar{{\varvec{c}}}} = 0\), \({\bar{{\varvec{b}}}} = {\bar{{\varvec{F}}}}{\varvec{a}}=\iota {\bar{{\varvec{f}}}}\), and \({\bar{{\varvec{c}}}} = {\displaystyle \frac{\partial {\bar{{\varvec{f}}}}}{\partial {\varvec{x}}}}{\varvec{a}}\).

The evaluation of stress tensors requires spectral the Lagrangian spectral tensor components of \({\displaystyle \frac{\partial W}{\partial {\varvec{C}}}}\) i.e.,

The Eulerian spectral components of the Cauchy stress \({\varvec{T}}\) for an incompressible body with respect to the spectral Eulerian basis \(\{ {\varvec{v}}_1,{\varvec{v}}_2,{\varvec{v}}_3\}\) are

In this paper, following the work of Shariff22, we nominate a prototype strain energy function that satisfy the P-property. We emphasize that our proposed non-linear strain energy function is consistent with the theory of infinitesimal elasticity. To ensure this consistency, we start our non-linear prototype construction via developing its infinitesimal strain energy counterpart.

Before we construct strain energy prototypes for finite strain deformation, we give a brief description on infinitesimal elasticity. When the gradient of the displacement field \({\varvec{u}}\) is very small

where \(\Vert \bullet \Vert\) is an appropriate norm and the magnitude of e is much less than unity. Up to O(e),

where \({\varvec{E}}\) is the infinitesimal strain. The most general quadratic form of the strain energy function is

where

where \(\mu , \mu _1, \mu _2, \kappa _1, \kappa _2,\kappa _3\) are ground-state material constants and their limitations are given in Online Appendix A.

We propose finite-strain energy function that is consistent with its infinitesimal counterpart. This can be easily done, following the work of Shariff22, by extending the above infinitesimal strain energy function using spectral generalized strains for finite deformations. The proposed strain energy function is

where

with the properties22

We could also include the following property, when appropriate, \(r_\alpha\) to represent physical strain measures with the extreme deformation values

We could easily extend (21) to (23) to construct a more general strain energy function (see for example22), but the strain energy function proposed in the Section should suffice to illustrate our model.

To illustrate our theory, we consider two simple deformations, pure bending and finite torsion of a right circular cylinder, where their displacements are known. For boundary value problems, where the displacements are unknown, the construction of solutions are described in Online Appendix B.

To plot the results in this section, for simplicity, we use

and the ground-state values

are those associated with skeletal muscle tissue10,25. Since our model is new and there are no experimental values for the following bending stiffness ground-state constants, we use the ad hoc values

to plot the graphs. Take note that the above values satisfy the restrictions given in Appendix A.

Bending of a rectangular block into a sector of a cylindrical tube.

Consider the problem of pure bending in plane strain, depicted in Fig. 3, in which a rectangular slab of incompressible material is bent into a sector of a circular annulus defined by

where \((r,\theta ,z)\) is the cylindrical polar coordinate for the current configuration and \((x_1,x_2,x_3)\) is the Cartesian referential coordinate with the basis \(\{ {\varvec{g}}_1 , {\varvec{g}}_2, {\varvec{g}}_3 ={\varvec{e}}_z \}\).

The formula employed here could be used to compare our theory with experiment (for example, a three point bending test experiment described in reference26).

The deformation tensor has the form

From the incompressibility condition \(\det {\varvec{F}}=1\) and the boundary conditions \(\theta (0)=0\) and \(r(A)=a\) we obtain

where \(r(B) = b\). Hence, in view of (3), (30) and (31), we have

and the spectral basis vectors are \({\varvec{u}}_i={\varvec{g}}_i\), \({\varvec{v}}_1={\varvec{e}}_r\), \({\varvec{v}}_2={\varvec{e}}_\theta\) and \({\varvec{v}}_3={\varvec{e}}_z\).

In this section we study the case \({\varvec{a}}={\varvec{g}}_2\), hence, \(a_1=a_3=0\) and \(a_2=1\). If we let \({\bar{{\varvec{F}}}}={\varvec{F}}\), we get

The strain energy function is simplified, i.e.

The non-zero Cauchy stress components simply becomes

where \(\sigma _1=\sigma _{rr}\), \(\sigma _2=\sigma _{\theta \theta }\) and \(\sigma _3=\sigma _{zz}\) are cylindrical components of the Cauchy stress. Since \(\sigma _i\) depends only on r, the equilibrium equation simply becomes

If we assume that \(\sigma _{rr} = 0\) at \(r = b\), we then have

Hence, we can evaluate

and with the above expression for p we obtain the stress-strain relations for \(\sigma _{\theta \theta }\) and \(\sigma _{zz}\). The bending moment \({\mathcal {M}}\), and the normal force \({\mathcal {N}}\), per unit length in the \(x_3\) direction, and applied to a section of constant \(\theta\), are

In Figs. 4 and 5, the behaviours of, respectively, the radial and hoop stresses are depicted using \({\displaystyle \frac{\chi }{B}}=1\) and the material is deformed to \({\displaystyle \frac{a}{B}}=1\). It is clear from these figures the magnitude of the stresses is higher for an elastic solid with fibre bending resistance than for a solid with perfectly flexible fibres.

Radial behaviour of stress \(\sigma _{rr}\). (a) Elastic solid with fibre bending resistance. (b) Elastic solid with no fibre bending resistance.

Radial behaviour of stress \(\sigma _{\theta \theta }\). (a) Elastic solid with fibre bending resistance. (b) Elastic solid with no fibre bending resistance.

The \({\mathcal {M}}\) values for a material with and without fibre bending resistance are, respectively, 46.44514245 kPaM\(^2\) and 35.55851694 kPaM\(^2\). The \({\mathcal {N}}\) values for a material with and without fibre bending resistance are, respectively, 30.58637503 kPaM and 23.29228593 kPaM. Hence, bending stiffness increases the magnitude of \({\mathcal {M}}\) and \({\mathcal {N}}\).

In this section we consider an incompressible thick-walled circular cylindrical annulus with the initial geometry

where R, \(\Theta\) and Z are reference polar coordinates with the corresponding basis \(B_R=\{ {\varvec{E}}_R,{\varvec{E}}_\Theta ,{\varvec{E}}_Z \}\). The boundary value problem illustrated here could be used in an experiment (see, for example, reference27) to verify our theoretical predictions.

The deformation is depicted in Fig. 6 and is described by

where \(\tau\) is the amount of torsional twist per unit deformed length and \(\lambda _z\) is the axial stretch. In the above formulation, r, \(\theta\) and z are cylindrical polar coordinates in the deformed configuration with the corresponding basis \(B_C=\{ {\varvec{e}}_r,{\varvec{e}}_\theta ,{\varvec{e}}_z \}\). Here, we have allowed \({\varvec{e}}_r={\varvec{E}}_R\), \({\varvec{e}}_\theta ={\varvec{E}}_\Theta\) and \({\varvec{e}}_z={\varvec{E}}_Z\). The deformation gradient is

where \(\gamma =r\tau\) and in this paper, we only consider \(\lambda _z \ge 1\). The Lagrangian principal directions are:

where

with

In the case of pure torsion, \(\lambda _z=1\) and we have \({\hat{\gamma }}=\gamma\). The principal stretches for a combined extension and torsion deformation are

Torsion and extension of a cylinder.

In this section we consider the case when \({\varvec{a}}={\varvec{E}}_z\), hence, \(a_1=0\), \(a_2=s\) and \(a_3=c\). If we let \({\bar{{\varvec{F}}}}={\varvec{F}}\) and using

we obtain

The strain energy function then takes the form

The Cauchy stress

In view of \({\varvec{a}}\equiv [0,0,1]^T\), we have \(a_1=0\), \(a_2=s\) and \(a_3=c\) and

where

The normal force \({\mathcal {N}}\) and the torque per unit deformed area \({\mathcal {M}}\) applied at the ends of the cylinder are as follows:

To remove the hydrostatic pressure term in (54)\(_1\), we use the equilibrium relation

and reformulate (54)\(_1\) in the form

It is clear from Fig. 7 that, for an axial stretch \(\lambda _z=1.5\), we require more torque to twist an elastic solid cylinder with fibre bending stiffness.

Torque, \({\mathcal {M}}\) versus \(\tau\). (a) Elastic solid with fibre bending stiffness. (b) Elastic solid with no fibre bending stiffness. \(\lambda _z=1.5\).

We have modelled elastic resistance due to changes in the curvature of the fibres without using the second gradient theory. In view of this, the proposed hyperelastic model is simple and does not contain couple stresses (which is required in a second gradient model). Hence, the proposed model is more realistic in the sense that a carbon fiber reinforced polymer is a non-polar material, where couple stresses do not exist. In the near future, FEM solutions of the proposed model will be obtained and we will extend this model to polymers that are reinforced with a family of two fibres.

All data generated or analysed during this study are included in this published article [and its supplementary information files].

Chukov, D. et al. Structure, mechanical and thermal properties of polyphenylene sulfide and polysulfone impregnated carbon fiber composites. Polymers 11, 684 (2019).

Article CAS PubMed PubMed Central Google Scholar

Linul, E., Lell, D., Movahedi, N., Codrean, C. & Fiedler, T. Compressive properties of Zinc Syntactic Foams at elevated temperatures. Compos. Part B Eng. 167, 122–134 (2019).

Article CAS Google Scholar

Sherif, G., Chukov, D., Tcherdyntsev, V. & Torokhov, V. Effect of formation route on the mechanical properties of the polyethersulfone composites reinforced with glass fibers. Polymers 11, 1364 (2019).

Article CAS PubMed PubMed Central Google Scholar

Yashas Gowda, T. G. et al. Polymer matrix-natural fiber composites: an overview. Cogent. Eng. 5, 1446667 (2018).

Article Google Scholar

Clyne, T. W. & Hull, D. An Introduction to Composite Materials 3rd edn. (Cambridge University Press, 2019).

Book Google Scholar

Monteiro, S. N. et al. Fique fabric: a promising reinforcement for polymer composites. Polymers 10, 246 (2018).

Article Google Scholar

Movahedi, N. & Linul, E. Quasi-static compressive behavior of the ex-situ aluminum-alloy foam-filled tubes under elevated temperature conditions. Mater. Lett. 206, 182–184 (2017).

Article CAS Google Scholar

Zagho, M. M., Hussein, E. A. & Elzatahry, A. A. Recent overviews in functional polymer composites for biomedical applications. Polymers 10, 739 (2018).

Article PubMed PubMed Central Google Scholar

Shariff, M. H. B. M. Physical invariant strain energy function for passive myocardium. J. Biomech. Model. Mechanobiol. 12(2), 215–223 (2012).

Article Google Scholar

Shariff, M. H. B. M. On the spectral constitutive modelling of transversely isotropic soft tissue: physical invariants. Int. J. Eng. Sci. 120, 199–219 (2017).

Article MathSciNet MATH Google Scholar

Pipkin, A. C. Stress analysis for fiber-reinforced materials. Adv. Appl. Mech. 19, 1–51 (1979).

Article MATH Google Scholar

Shariff, M. H. B. M. & Merodio, J. Residually stressed fiber reinforced solids: a spectral approach. Materials 13(18), 4076 (2020).

Article ADS PubMed PubMed Central MATH Google Scholar

Spencer, A. J. M. Deformations of Fibre-Reinforced Materials (Oxford University Press, 1972).

MATH Google Scholar

Shariff, M. H. B. M., Merodio, J. & Bustamante, R. Finite deformations of fibre bending stiffness: a spectral approach. J. Appl. Comput. Mech. 8(4), 1332–1342 (2022).

Google Scholar

Shariff, M. H. B. M., Merodio, J. & Bustamante, R. Nonlinear elastic constitutive relations of residually stressed composites with stiff curved fibres. Appl. Math. Mech. 43(10), 1515–1530 (2022).

Article MathSciNet MATH Google Scholar

Soldatos, K. P., Shariff, M. H. B. M. & Merodio, J. On the constitution of polar fibre-reinforced materials. Mech. Adv. Mater. Struct. 28(21), 2255–2266 (2020).

Article Google Scholar

Spencer, A. J. M. & Soldatos, K. P. Finite deformations of fibre-reinforced elastic solids with fibre bending stiffness. Int. J. Non-Linear Mech. 42, 355–368 (2007).

Article ADS Google Scholar

Steigmann, D. Theory of elastic solids reinforced with fibers resistant to extension, flexure and twist. Int. J. Non-Linear Mech. 47, 734–742 (2012).

Article ADS Google Scholar

Hadjesfandiari, A. R. & Dargush, G. F. Couple stress theory for solids. Int. J. Solids Struct. 48, 2496–2510 (2011).

Article Google Scholar

Hadjesfandiari, A. R. & Dargush, G. F. Evolution of Generalized Couple-stress Continuum Theories:acriticalanalysis. Preprint arXiv:1501.03112. (2015).

Neff, P., Münch, I., Ghiba, I. & Madeo, A. On some fundamental misunderstandings in the indeterminate couple stress model. A comment on recent papers of A.R. Hadjesfandiari and G.F. Dargush. Int. J. Solids Struct. 81, 233–243 (2016).

Article Google Scholar

Shariff, M. H. B. M. A generalized strain approach to anisotropic elasticity. Sci. Rep. 12(1), 1–22 (2022).

Article Google Scholar

Shariff, M. H. B. M. On the smallest number of functions representing isotropic functions of scalars, vectors and tensors. Q. J. Mech. Appl. Math.https://doi.org/10.1093/qjmam/hbac022 (2023).

Article Google Scholar

Shariff, M. H. B. M. Anisotropic separable free energy functions for elastic and non-elastic solids. Acta Mech. 227(11), 3213–3237 (2016).

Article MathSciNet MATH Google Scholar

Morrow, D. A., Donahue, T. L. H., Odegard, G. M. & Kaufman, K. R. Transversely isotropic tensile material properties of skeletal muscle tissue. J. Mech. Behav. Biomed. Mater. 3, 124–129 (2010).

Article PubMed Google Scholar

Mathieu, S., Hamila, N., Bouillon, F. & Boisse, P. Enhanced modeling of 3D composite preform deformations taking into account local fiber bending stiffness. Compos. Sci. Technol. 117, 322–333 (2015).

Article CAS Google Scholar

Lu, W.-Y. et al. Solid cylinder torsion for large shear deformation and failure of engineering materials. Exp. Mech. 61, 307–320 (2021).

Article CAS Google Scholar

Download references

Department of Mathematics, Khalifa University of Science and Technology, Abu Dhabi, UAE

M. H. B. M. Shariff

Departamento de Matemática Aplicada a las TIC, ETS de Ingeniería de Sistemas Informáticos, Universidad Politécnica de Madrid, 28031, Madrid, Spain

J. Merodio

Departamento de Ingeniería Mecánica, Universidad de Chile, Beauchef 851, Santiago Centro, Santiago, Chile

R. Bustamante

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

M.H.B.M.S. Writing-original draft, Writing-review, and editing. All authors reviewed the manuscript.

Correspondence to M. H. B. M. Shariff.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Shariff, M.H.B.M., Merodio, J. & Bustamante, R. A non-second-gradient model for nonlinear elastic bodies with fibre stiffness. Sci Rep 13, 6562 (2023). https://doi.org/10.1038/s41598-023-33670-6

Download citation

Received: 03 February 2023

Accepted: 17 April 2023

Published: 21 April 2023

DOI: https://doi.org/10.1038/s41598-023-33670-6

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.