Continuous Softening as a State of Hyperelasticity: Examples of Application to the Softening Behavior of the Brain Tissue.

The continuous softening behavior of the brain tissue, i.e., the softening in the primary loading path with an increase in deformation, is modeled in this work as a state of hyperelasticity up to the onset of failure. That is, the softening behavior is captured via a core hyperelastic model without the addition of damage variables and/or functions. Examples of the application of the model will be provided to extant datasets of uniaxial tension and simple shear deformations, demonstrating the capability of the model to capture the whole-range deformation of the brain tissue specimens, including their softening behavior. Quantitative and qualitative comparisons with other models within the brain biomechanics literature will also be presented, showing the clear advantages of the current approach. The application of the model is then extended to capturing the rate-dependent softening behavior of the tissue by allowing the parameters of the core hyperelastic model to evolve, i.e., vary, with the deformation rate. It is shown that the model captures the rate-dependent and softening behaviors of the specimens favorably and also predicts the behavior at other rates. These results offer a clear set of advantages in favor of the considered modeling approach here for capturing the quasi-static and rate-dependent mechanical properties of the brain tissue, including its softening behavior, over the existing models in the literature, which at best may purport to capture only a reduced set of the foregoing behaviors, and with ill-posed effects.

Modelling the rate-dependent mechanical behaviour of the brain tissue.

A new modelling approach is employed in this work for application to the rate-dependent mechanical behaviour of the brain tissue, as an incompressible isotropic material. Extant datasets encompassing single- and multi-mode compression, tension and simple shear deformation(s) are considered, across a wide range of deformation rates from quasi-static to rates akin to blast loading conditions, in the order of 1000 s-1 . With a simple functional form and a reduced number of parameters, the model is shown to capture the considered rate-dependent behaviours favourably, including in both single- and multi-mode deformation fits, and over all range of deformation rates. The provided modelling results here are obtained from either first fitting the model to the quasi-static data, or/and predicting the behaviour at a different rate than those used for calibrating the model parameters. Given its simplicity, versatility, predictive capability and accuracy, the application of the utilised modelling framework in this work to the rate-dependent mechanical behaviour of the brain tissue is proposed.

Modelling brain tissue elasticity with the Ogden model and an alternative family of constitutive models†.

The Ogden model is often considered as a standard model in the literature for application to the deformation of brain tissue. Here, we show that, in some of those applications, the use of the Ogden model leads to the non-convexity of the strain-energy function and mis-prediction of the correct concavity of the experimental stress-stretch curves over a range of the deformation domain. By contrast, we propose a family of models which provides a favourable fit to the considered datasets while remaining free from the highlighted shortcomings of the Ogden model. While, as we discuss, those shortcomings might be due to the artefacts of the testing protocols, the proposed family of models proves impervious to such artefacts. This article is part of the theme issue 'The Ogden model of rubber mechanics: Fifty years of impact on nonlinear elasticity'.

The Ogden model of rubber mechanics: 50 years of impact on nonlinear elasticity.

We place the Ogden model of rubber elasticity, published in Proceedings of the Royal Society 50 years ago, in the wider context of the theory of nonlinear elasticity. We then follow with a short interview of Ray Ogden FRS and introduce the papers collected for this Theme Issue. This article is part of the theme issue 'The Ogden model of rubber mechanics: Fifty years of impact on nonlinear elasticity'.

Shear waves in a nonlinear relaxing media: A three-dimensional perspective.

An important one-dimensional rheological model for the propagation of a linearly polarized shear wave was recently obtained and proposed by Cormack and Hamilton [(2018). J. Acoust. Soc. Am. 143(2), 1035-1048]. We show that it is possible to embed such a result within a wider and complete set of general three-dimensional models derived within the theoretical framework of rigorous continuum mechanics. We show that, following this approach, we are able to derive in a simple and straightforward way the equations that govern the propagation of circularly polarized shear waves. When the phase of such waves is constant, we find the same equation for linearly polarized shear waves already proposed elsewhere. Moreover, we show that, under appropriate asymptotic assumptions, our results are indifferent with respect to the choice of the objective time derivative used in the constitutive class.

Managing awareness can avoid hysteresis in disease spread: an application to coronavirus Covid-19.

A SEIR-type model is investigated to evaluate the effects of awareness campaigns in the presence of factors that can induce overexposure to disease. We find that high levels of overexposure can drive system dynamics towards a backward phenomenology and that increasing people awareness through balanced and aware information can be crucial to avoid dangerous dynamical transitions as hysteresis or transient oscillations before disease eradication. Investigations in the time dependent regimes are provided to support the results. Google Trends data in the context of Covid19 are also used to stress how low levels of awareness, combined with high overexposure, can be related to recent episodes of epidemic resurgence in Europe. Our results suggest that the interplay between overexposure and awareness is a point that should not be underestimated both in the current and future management of the Covid19 emergency.

Generalization of the Zabolotskaya equation to all incompressible isotropic elastic solids.

We study elastic shear waves of small but finite amplitude, composed of an anti-plane shear motion and a general in-plane motion. We use a multiple scales expansion to derive an asymptotic system of coupled nonlinear equations describing their propagation in all isotropic incompressible nonlinear elastic solids, generalizing the scalar Zabolotskaya equation of compressible nonlinear elasticity. We show that for a general isotropic incompressible solid, the coupling between anti-plane and in-plane motions cannot be undone and thus conclude that linear polarization is impossible for general nonlinear two-dimensional shear waves. We then use the equations to study the evolution of a nonlinear Gaussian beam in a soft solid: we show that a pure (linearly polarized) shear beam source generates only odd harmonics, but that introducing a slight in-plane noise in the source signal leads to a second harmonic, of the same magnitude as the fifth harmonic, a phenomenon recently observed experimentally. Finally, we present examples of some special shear motions with linear polarization.

Linearly polarized waves of finite amplitude in pre-strained elastic materials.

We study the propagation of linearly polarized transverse waves in a pre-strained incompressible isotropic elastic solid. Both finite and small-but-finite amplitude waves are examined. Irrespective of the magnitude of the wave amplitude, these waves may propagate only if the (unit) normal to the plane spanned by the directions of propagation and polarization is a principal direction of the left Cauchy-Green deformation tensor associated with the pre-strained state. A rigorous asymptotic analysis of the equations governing the propagation of waves of small but finite amplitude reveals that the time scale over which the nonlinear effects become significant depends on the direction along which the wave travels. Moreover, we design theoretically an experimental procedure to determine the Landau constants of the fourth-order weakly nonlinear theory of elasticity.

Parametric resonance in DNA.

We consider a simple mesoscopic model of DNA in which the binding of the RNA polymerase enzyme molecule to the promoter sequence of the DNA is included through a substrate energy term modeling the enzymatic interaction with the DNA strands. We focus on the differential system for solitary waves and derive conditions--in terms of the model parameters--for the occurrence of the parametric resonance phenomenon. We find that what truly matters for parametric resonance is not the ratio between the strength of the stacking and the inter-strand forces but the ratio between the substrate and the inter-strands. On the basis of these results, the standard objection that longitudinal motion is negligible because of the second order seems to fail, suggesting that all the studies involving the longitudinal degree of freedom in DNA should be reconsidered when the interaction of the RNA polymerase with the DNA macromolecule is not neglected.

Proper formulation of viscous dissipation for nonlinear waves in solids.

To model nonlinear viscous dissipative motions in solids, acoustical physicists usually add terms linear in E, the material time derivative of the Lagrangian strain tensor E, to the elastic stress tensor σ derived from the expansion to the third (sometimes fourth) order of the strain energy density E=E(tr E,tr E(2),tr E(3)). Here it is shown that this practice, which has been widely used in the past three decades or so, is physically wrong for at least two reasons and that it should be corrected. One reason is that the elastic stress tensor σ is not symmetric while E is symmetric, so that motions for which σ+σ(T)≠0 will give rise to elastic stresses that have no viscous pendant. Another reason is that E is frame-invariant, while σ is not, so that an observer transformation would alter the elastic part of the total stress differently than it would alter the dissipative part, thereby violating the fundamental principle of material frame indifference. These problems can have serious consequences for nonlinear shear wave propagation in soft solids as seen here with an example of a kink in almost incompressible soft solids.

Automated estimation of collagen fibre dispersion in the dermis and its contribution to the anisotropic behaviour of skin.

Collagen fibres play an important role in the mechanical behaviour of many soft tissues. Modelling of such tissues now often incorporates a collagen fibre distribution. However, the availability of accurate structural data has so far lagged behind the progress of anisotropic constitutive modelling. Here, an automated process is developed to identify the orientation of collagen fibres using inexpensive and relatively simple techniques. The method uses established histological techniques and an algorithm implemented in the MATLAB image processing toolbox. It takes an average of 15 s to evaluate one image, compared to several hours if assessed visually. The technique was applied to histological sections of human skin with different Langer line orientations and a definite correlation between the orientation of Langer lines and the preferred orientation of collagen fibres in the dermis (p < 0.001, R(2) = 0.95) was observed. The structural parameters of the Gasser-Ogden-Holzapfel (GOH) model were all successfully evaluated. The mean dispersion factor for the dermis was κ = 0.1404±0.0028. The constitutive parameters µ, k(1) and k(2) were evaluated through physically-based, least squares curve-fitting of experimental test data. The values found for µ, k(1) and k(2) were 0.2014 MPa, 243.6 and 0.1327, respectively. Finally, the above model was implemented in ABAQUS/Standard and a finite element (FE) computation was performed of uniaxial extension tests on human skin. It is expected that the results of this study will assist those wishing to model skin, and that the algorithm described will be of benefit to those who wish to evaluate the collagen dispersion of other soft tissues.

Third- and fourth-order constants of incompressible soft solids and the acousto-elastic effect.

Acousto-elasticity is concerned with the propagation of small-amplitude waves in deformed solids. Results previously established for the incremental elastodynamics of exact non-linear elasticity are useful for the determination of third- and fourth-order elastic constants, especially in the case of incompressible isotropic soft solids, where the expressions are particularly simple. Specifically, it is simply a matter of expanding the expression for ρv(2), where ρ is the mass density and v the wave speed, in terms of the elongation e of a block subject to a uniaxial tension. The analysis shows that in the resulting expression: ρv(2) = a+be+ce(2), say, a depends linearly on µ; b on µ and A; and c on µ, A, and D, the respective second-, third, and fourth-order constants of incompressible elasticity, for bulk shear waves and for surface waves.

Weierstrass's criterion and compact solitary waves.

Weierstrass's theory is a standard qualitative tool for single degree of freedom equations, used in classical mechanics and in many textbooks. In this Brief Report we show how a simple generalization of this tool makes it possible to identify some differential equations for which compact and even semicompact traveling solitary waves exist. In the framework of continuum mechanics, these differential equations correspond to bulk shear waves for a special class of constitutive laws.

Introducing mesoscopic information into constitutive equations for arterial walls.

We propose a new elastic constitutive law for arterial tissue in which the limiting polymeric chain extensibility of both collagen and elastin fibres is accounted for. The elastic strain-energy function is separated additively into two parts: an isotropic contribution associated with the matrix (incorporating the elastin fibre network) and an anisotropic one associated with the collagen fibres. Information on the limiting extensibility in each case provides some mesoscopic input into the model. The (logarithm-based) model is compared with the Fung-Demiray exponential model and certain other recently proposed models. Some aspects of the elastic response under extension and inflation of a thin-walled circular cylindrical tube (the artery) are then examined and compared with the corresponding response of a rubber-like tube. We point out that our model, when both isotropic and anisotropic terms are included, can be developed to accommodate changing mechanical properties associated with degradation of the elastin and collagen by considering the material constants that define the limit of chain extensibility to evolve in time.

The relevance of nonlinear stacking interactions in simple models of double-stranded DNA.

Single molecule DNA experiments provide interesting data that allow a better understanding of the mechanical interactions between the strands and the nucleotides of this molecule. In some sense, these experiments complement the classical ones about DNA thermal denaturation. It is well known that the original Peyrard-Bishop (PB) model by means of a harmonic stacking potential and a nonlinear substrate potential has been able to predict the existence of a critical temperature of full denaturation of the molecule. In the present paper, driven by the findings of single molecule experiments, we substitute the original harmonic intra-strand stacking potential with a Duffing type potential. By elementary and analytical arguments, we show that with this choice it is possible to obtain a sharp transition in the classical domain wall solution of the PB model and the compactification of the classical solitary wave solutions of other models for the dynamics of DNA. We discuss why these solutions may improve our knowledge of the DNA dynamics in several directions.

Solitary and compactlike shear waves in the bulk of solids.

We show that a model proposed by Rubin, Rosenau, and Gottlieb [J. Appl. Phys. 77, 4054 (1995)] for dispersion caused by an inherent material characteristic length belongs to the class of simple materials. Therefore, it is possible to generalize the idea of Rubin, Rosenau, and Gottlieb to include a wide range of material models, from nonlinear elasticity to turbulence. Using this insight, we are able to fine-tune nonlinear and dispersive effects in the theory of nonlinear elasticity in order to generate pulse solitary waves and also bulk traveling waves with compact support.

Finite amplitude elastic waves propagating in compressible solids.

The paper studies the interaction of a longitudinal wave with transverse waves in general isotropic and unconstrained hyperelastic materials, including the possibility of dissipation. The dissipative term chosen is similar to the classical stress tensor describing a Stokesian fluid and is commonly used in nonlinear acoustics. The aim of this research is to derive the corresponding general equations of motion, valid for any possible form of the strain energy function and to investigate the possibility of obtaining some general and exact solutions to these equations by reducing them to a set of ordinary differential equations. Then the reductions can lead to some exact closed-form solutions for special classes of materials (here the examples of the Hadamard, Blatz-Ko, and power-law strain energy densities are considered, as well as fourth-order elasticity). The solutions derived are in a time-space separable form and may be interpreted as generalized oscillatory shearing motions and generalized sinusoidal standing waves. By means of standard methods of dynamical systems theory, some peculiar properties of waves propagating in compressible materials are uncovered, such as for example, the emergence of destabilizing effects. These latter features exist for highly nonlinear strain energy functions such as the relatively simple power-law strain energy, but they cannot exist in the framework of fourth-order elasticity.