Search

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 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.

We investigate how two finite-amplitude, transverse, plane body waves may be superposed to propagate in a deformed hyperelastic incompressible solid. We find that the equations of motion reduce to a well-determined system of partial differential equations, making the motion controllable for all solids. We find that in deformed Mooney-Rivlin materials, they may travel along any direction and be polarised along any transverse direction, an extension of a result by Boulanger and Hayes [Quart. J. Mech. Appl. Math. 45 (1992) 575]. Furthermore, their motion is governed by a linear system of partial differential equations, making the Mooney-Rivlin special in that respect. We select another model to show that for other materials, the equations are nonlinear. We use asymptotic equations to reveal the onset of nonlinearity for the waves, paying particular attention to how close the propagation direction is to the principal axes of pre-deformation.

Soft materials such as rubbers, silicones, gels and biological tissues have a nonlinear response to large deformations, a phenomenon which in principle can be captured by hyperelastic models. The suitability of a candidate hyperelastic strain energy function is then determined by comparing its predicted response to the data gleaned from tests and adjusting the material parameters to get a good fit, an exercise which can be deceptive because of nonlinearity. Here we propose to generalise the approach of Rivlin and Saunders [Phil Trans A 243 (1951) 251-288] who, instead of reporting the data as stress against stretch, manipulated these measures to create the 'Mooney plot', where the Mooney-Rivlin model is expected to produce a linear fit. We show that extending this idea to other models and modes of deformation (tension, shear, torsion, etc.) is advantageous, not only (a) for the fitting procedure, but also to (b) delineate trends in the deformation which are not obvious from the raw data (and may be interpreted in terms of micro-, meso-, and macro-structures) and (c) obtain a bounded condition number \k{appa} over the whole range of deformation; a robustness which is lacking in other plots and spaces.

Background: Needle-based procedures such as fine needle aspiration (FNA) and thermal ablation, are often applied for thyroid nodule diagnosis and therapeutic purposes, respectively. With blood vessels and nerves nearby, these procedures can pose risks in damaging surrounding critical structures. Purpose: The development and validation of innovative strategies to manage these risks require a test object with well-characterized physical properties. For this work, we focus on the application of ultrasound-guided thermal radio-frequency ablation (RFA). Methods: We have developed an anthropomorphic phantom mimicking the thyroid and surrounding anatomical and physiological structures that are relevant to ultrasound-guided thermal ablation. The phantom was composed of a mixture of polyacrylamide, water, and egg white extract and was cast using molds in multiple steps. The thermal, acoustical, and electrical characteristics were experimentally validated. The ablation zones were analyzed via non-destructive T2-weighted MRI scans utilizing the relaxometry changes of coagulated egg albumen, and the temperature distribution was monitored using an array of fiber Bragg sensors. Results: The physical properties of the phantom were verified both on ultrasound as well as its response to thermal ablation. The final temperature achieved (92{\deg}C), the median percentage of the nodule ablated (82.1%), the median volume ablated outside the nodule (0.8 mL), and the median number of critical structures affected (0) were quantified. Conclusion: An anthropomorphic phantom that can provide a realistic model for development and training in ultrasound-guided needle-based thermal interventions for thyroid nodules has been presented. In the future, this model can also be extended to novel needle-based diagnostic procedures., Comment: 19 pages, 10 figures, 3 tables

We consider shear wave propagation in soft viscoelastic solids of rate type. Based on objective stress rates, the constitutive model accounts for finite strain, incompressibility, as well as stress- and strain-rate viscoelasticity. The theory generalises the standard linear solid model to three-dimensional volume-preserving motions of large amplitude in a physically-consistent way. The nonlinear equations governing shear motion take the form of a one-dimensional hyperbolic system with relaxation. For specific objective rates of Cauchy stress (lower- and upper-convected derivatives), we study the propagation of acceleration waves and shock waves. Then we show that both smooth and discontinuous travelling wave solutions can be obtained analytically. We observe that the amplitude and velocity of steady shocks are very sensitive to variations of the stress relaxation time. Furthermore, the existence of steady shocks is conditional. Extension of these results to the case of multiple relaxation mechanisms and of the Jaumann stress rate is attempted. The analysis of simple shearing motions is more involved in these cases.

The protocol of quantum reading refers to the quantum enhanced retrieval of information from an optical memory, whose generic cell stores a bit of information in two possible lossy channels. In the following we analyze the case of a particular class of optical receiver, based on photon counting measurement, since they can be particularly simple in view of real applications. We show that a quantum advantage is achievable when a transmitter based on two-mode squeezed vacuum (TMSV) states is combined with a photon counting receiver, and we experimentally confirm it. In this paper, after introducing some theoretical background, we focus on the experimental realisation, describing the data collection and the data analysis in detail.

We show that porcine brain matter can be modelled accurately as a very soft rubber-like material using the Mooney-Rivlin strain energy function, up to strains as high as 60\%. This result followed from simple shear experiments performed on small rectangular fresh samples ($2.5$ cm$^3$ and $1.1$ cm$^3$) at quasi-static strain rates. They revealed a linear shear stress--shear strain relationship ($R^2> 0.97 $), characteristic of Mooney-Rivlin materials at large strains. We found that porcine brain matter is about 30 times less resistant to shear forces than a silicone gel. We also verified experimentally that brain matter exhibits the positive Poynting effect of nonlinear elasticity, and numerically that the stress and strain fields remain mostly homogeneous throughout the thickness of the samples in simple shear.

We find the strain energy function for isotropic incompressible solids exhibiting a linear relationship between shear stress and amount of shear, and between torque and amount of twist, when subject to large simple shear or torsion deformations. It is inclusive of the well-known neo-Hookean and the Mooney-Rivlin models, but also can accommodate other terms, as certain arbitrary functions of the principal strain invariants. Effectively, the extra terms can be used to account for several non-linear effects observed experimentally but not captured by the neo-Hookean and Mooney-Rivlin models, such as strain stiffening effects due to limiting chain extensibility.

The modelling of off-axis simple tension experiments on transversely isotropic nonlinearly elastic materials is considered. A testing protocol is proposed where normal force is applied to one edge of a rectangular specimen with the opposite edge allowed to move laterally but constrained so that no vertical displacement is allowed. Numerical simulations suggest that this deformation is likely to remain substantially homogeneous throughout the specimen for moderate deformations. It is therefore further proposed that such tests can be modelled adequately as a homogenous deformation consisting of a triaxial stretch accompanied by a simple shear. Thus the proposed test should be a viable alternative to the standard biaxial tests currently used as material characterisation tests for transversely isotropic materials in general and, in particular, for soft, biological tissue. A consequence of the analysis is a kinematical universal relation for off-axis testing that results when the strain-energy function is assumed to be a function of only one isotropic and one anisotropic invariant, as is typically the case. The universal relation provides a simple test of this assumption, which is usually made for mathematical convenience. Numerical simulations also suggest that this universal relation is unlikely to agree with experimental data and therefore that at least three invariants are necessary to fully capture the mechanical response of transversely isotropic materials.

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 non-linear 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 polarised) 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 polarisation.

We propose a model for the mechanical behavior of protein materials. Based on a limited number of experimental macromolecular parameters (persistence and contour lengths, rate of unfolding dissipation) we obtain the macroscopic behavior of keratin fibers (human, cow, and rabbit hair), taking into account the damage and residual stretches effects which are fundamental in many functions of life. We support our theoretical results by showing that our model is robust and able to reproduce with high quantitive accuracy the cyclic experimental behavior of different keratinous protein materials we tested. We also show the capability of describing, even if with lower precision, the dissipation and permanent strain effects in spider silks.

In mean-field theory, the non-local state of fluid molecules can be taken into account using a statistical method. The molecular model combined with a density expansion in Taylor series of the fourth order yields an internal energy value relevant to the fourth-gradient model, and the equation of isother-mal motions takes then density's spatial derivatives into account for waves travelling in both liquid and vapour phases. At equilibrium, the equation of the density profile across interfaces is more precise than the Cahn and Hilliard equation, and near the fluid's critical-point, the density profile verifies an Extended Fisher-Kolmogorov equation, allowing kinks, which converges towards the Cahn-Hillard equation when approaching the critical point. Nonetheless, we also get pulse waves oscillating and generating critical opalescence., Comment: Some travelling waves of density near a fluid critical-point

A great variety of models can describe the non-linear response of rubber to uni-axial tension. Yet an in-depth understanding of the successive stages of large extension is still lacking. We show that the response can be broken down in three steps, which we delineate by relying on a simple formatting of the data, the so-called Mooney transform. First, the small-to-moderate regime, where the polymeric chains unfold easily and the Mooney plot is almost linear. Second, the strain-hardening regime, where blobs of bundled chains unfold to stiffen the response in correspondence to the "upturn" of the Mooney plot. Third, the limiting-chain regime, with a sharp stiffening occurring as the chains extend towards their limit. We provide strain-energy functions with terms accounting for each stage, that (i) give an accurate local and then global fitting of the data; (ii) are consistent with weak non-linear elasticity theory; and (iii) can be interpreted in the framework of statistical mechanics. We apply our method to Treloar's classical experimental data and also to some more recent data. Our method not only provides models that describe the experimental data with a very low quantitative relative error, but also shows that the theory of non-linear elasticity is much more robust that seemed at first sight.

The theory of Lie remarkable equations, i.e. differential equations characterized by their Lie point symmetries, is reviewed and applied to ordinary differential equations. In particular, we consider some relevant Lie algebras of vector fields on $\mathbb{R}^k$ and characterize Lie remarkable equations admitted by the considered Lie algebras., Comment: 17 pages

Using the theory of $1+1$ hyperbolic systems we put in perspective the mathematical and geometrical structure of the celebrated circularly polarized waves solutions for isotropic hyperelastic materials determined by Carroll in Acta Mechanica 3 (1967) 167--181. We show that a natural generalization of this class of solutions yields an infinite family of \emph{linear} solutions for the equations of isotropic elastodynamics. Moreover, we determine a huge class of hyperbolic partial differential equations having the same property of the shear wave system. Restricting the attention to the usual first order asymptotic approximation of the equations determining transverse waves we provide the complete integration of this system using generalized symmetries., Comment: 19 pages

We propose a general multiscale approach for the mechanical behavior of three-dimensional networks of macromolecules undergoing strain-induced unfolding. Starting from a (statistically based) energetic analysis of the macromolecule unfolding strategy, we obtain a three-dimensional continuum model with variable natural configuration and an energy function analytically deduced from the microscale material parameters. The comparison with the experiments shows the ability of the model to describe the complex behavior, with residual stretches and unfolding effects, observed in different biological materials.

On one hand, classical Monte Carlo and molecular dynamics (MD) simulations have been very useful in the study of liquids in nanotubes, enabling a wide variety of properties to be calculated in intuitive agreement with experiments. On the other hand, recent studies indicate that the theory of continuum breaks down only at the nanometer level; consequently flows through nanotubes still can be investigated with Navier-Stokes equations if we take suitable boundary conditions into account. The aim of this paper is to study the statics and dynamics of liquids in nanotubes by using methods of non-linear continuum mechanics. We assume that the nanotube is filled with only a liquid phase; by using a second gradient theory the static profile of the liquid density in the tube is analytically obtained and compared with the profile issued from molecular dynamics simulation. Inside the tube there are two domains: a thin layer near the solid wall where the liquid density is non-uniform and a central core where the liquid density is uniform. In the dynamic case a closed form analytic solution seems to be no more possible, but by a scaling argument it is shown that, in the tube, two distinct domains connected at their frontiers still exist. The thin inhomogeneous layer near the solid wall can be interpreted in relation with the Navier length when the liquid slips on the boundary as it is expected by experiments and molecular dynamics calculations., Comment: 27 pages

Two questions related to elastic motions are raised and addressed. First: in which theoretical framework can the equations of motion be written for an elastic half-space put into uniform rotation? It is seen that nonlinear finite elasticity provides such a framework for incompressible solids. Second: how can finite amplitude exact solutions be generated? It is seen that for some finite amplitude transverse waves in rotating incompressible elastic solids with general shear response, the solutions are obtained by reduction of the equations of motion to a system of ordinary differential equations equivalent to the system governing the central motion problem of classical mechanics. In the special case of circularly-polarized harmonic progressive waves, the dispersion equation is solved in closed form for a variety of shear responses, including nonlinear models for rubberlike and soft biological tissues. A fruitful analogy with the motion of a nonlinear string is pointed out.

New exact solutions are exhibited within the framework of finite viscoelasticity. More precisely, the solutions correspond to finite-amplitude, transverse, linearly-polarized, inhomogeneous motions superposed upon a finite homogeneous static deformation. The viscoelastic body is composed of a Mooney-Rivlin viscoelastic solid, whose constitutive equation consists in the sum of an elastic part (Mooney-Rivlin hyperelastic model) and a viscous part (Newtonian viscous fluid model). The analysis shows that the results are similar to those obtained for the purely elastic case; inter alia, the normals to the planes of constant phase and to the planes of constant amplitude must be orthogonal and conjugate with respect to the B-ellipsoid, where B is the left Cauchy-Green strain tensor associated with the initial large static deformation. However, when the constitutive equation is specialized either to the case of a neo-Hookean viscoelastic solid or to the case of a Newtonian viscous fluid, a greater variety of solutions arises, with no counterpart in the purely elastic case. These solutions include travelling inhomogeneous finite-amplitude damped waves and standing damped waves.

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.

We propose a simple approach, based on the minimization of the total (entropic plus unfolding) energy of a two-state system, describing the stretch-induced unfolding of macromolecules (proteins, silks, nanopolymers, DNA/RNA). The model is fully analytical and enlightens the role of the different energetic components regulating the unfolding evolution. As an explicit example of application we compare the analytical results with the titin Atomic Force Microscopy experiments showing the ability of the model to quantitatively reproduce the mechanical behavior of macromolecules unfolding.

We present a phenomenological approach to dispersion in nonlinear elasticity. A simple, thermomechanically sound, constitutive model is proposed to describe the (non-dissipative) properties of a hyperelastic dispersive solid, without recourse to a microstructure or a special geometry. As a result, nonlinear and dispersive waves can travel in the bulk of such solids, and special waves emerge, some classic (periodic waves or pulse solitary waves of infinite extend), some exotic (kink or pulse waves of compact support). We show that for incompressible dispersive power-law solids and forth-order elasticity solids, solitary waves can however only exist in the case of linear transverse polarization. We also study the influence of pre-stretch and hardening. We provide links with other (quasi-continuum, asymptotic) theories; in particular, an appropriate asymptotic multiscale expansion specializes our exact equations of motion to the vectorial MKdV equation, for any hyperelastic material.

In order to model nonlinear viscous dissipative motions in solids, acoustical physicists usually add terms linear in dot{E}, the material time derivative of the Lagrangian strain tensor E, to the elastic stress tensor sigma derived from the expansion to the third- (sometimes fourth-) order of the strain energy density e=e(trace(E), trace(E^2), trace(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 sigma is not symmetric while dot{E) is symmetric, so that motions for which sigma + sigma^T <> 0 will give rise to elastic stresses which have no viscous pendant. Another reason is that dot{E} is frame-invariant, while sigma 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.