BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michel Destrade

Michel Destrade appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softBiological Physicsphysics.med-phnlin.PSTissues and Organsphysics.class-phcond-mat.mtrl-scimath-phmath.MPphysics.app-ph

Catalog footprint

What is connected

98works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

98 published item(s)

preprint2022arXiv

Non-destructive mapping of stress, strain and stiffness of thin elastically deformed materials

Knowing the stress within a soft material is of fundamental interest to basic research and practical applications, such as soft matter devices, biomaterial engineering, and medical sciences. However, it is challenging to measure stress fields in situ in a non-invasive way. It becomes even more difficult if the mechanical properties of the material are unknown or altered by the stress. Here we present a robust non-destructive technique capable of measuring in situ stress and strain in elastically deformed thin films without the need to know their material properties. The technique is based on measuring elastic wave speeds, and then using a universal dispersion curve we derived for Lamb wave to predict the local stress and strain. Using optical coherence tomography, we experimentally verified the method for a rubber sheet, a cling film, and the leather skin of a musical instrument.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Representing the stress and strain energy of elastic solids with initial stress and transverse texture anisotropy

Real-world solids, such as rocks, soft tissues, and engineering materials, are often under some form of stress. Most real materials are also, to some degree, anisotropic due to their microstructure, a characteristic often called the `texture anisotropy'. This anisotropy can stem from preferential grain alignment in polycrystalline materials, aligned micro-cracks, or structural reinforcement, such as collagen bundles in biological tissues, steel rods in prestressed concrete and reinforcing fibres in composites. Here we establish a framework for initially stressed solids with transverse texture anisotropy. We consider that the strain energy per unit mass of the reference is an explicit function of the elastic deformation gradient, the initial stress tensor, and the texture anisotropy. We determine the corresponding constitutive relations and develop examples of nonlinear strain energies which depend explicitly on the initial stress and direction of texture anisotropy. As an application, we then employ these models to analyse the stress distribution of an inflated initially stressed cylinder with texture anisotropy, and the tension of a welded metal plate. We also deduce the elastic moduli needed to describe linear elasticity from stress reference with transverse texture anisotropy. As an example we show how to measure the stress with small-amplitude shear waves.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A combined experimental and numerical study of stab-penetration forces

The magnitude of force used in a stabbing incident can be difficult to quantify, although the estimate given by forensic pathologists is often seen as `critical' evidence in medico-legal situations. The main objective of this study is to develop a quantitative measure of the force associated with a knife stabbing biological tissue, using a combined experimental and numerical technique. A series of stab-penetration tests were performed to quantify the force required for a blade to penetrate skin at various speeds and using different `sharp' instruments. A computational model of blade penetration was developed using ABAQUS/EXPLICIT, a non-linear finite element analysis (FEA) commercial package. This model, which incorporated element deletion along with a suitable failure criterion, is capable of systematically quantifying the effect of the many variables affecting a stab event. This quantitative data could, in time, lead to the development of a predictive model that could help indicate the level of force used in a particular stabbing incident.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A robust anisotropic hyperelastic formulation for the modelling of soft tissue

The Holzapfel-Gasser-Ogden (HGO) model for anisotropic hyperelastic behaviour of collagen fibre reinforced materials was initially developed to describe the elastic properties of arterial tissue, but is now used extensively for modelling a variety of soft biological tissues. Such materials can be regarded as incompressible, and when the incompressibility condition is adopted the strain energy Ψof the HGO model is a function of one isotropic and two anisotropic deformation invariants. A compressible form (HGO-C model) is widely used in finite element simulations whereby the isotropic part of Ψis decoupled into volumetric and isochoric parts and the anisotropic part of Ψis expressed in terms of isochoric invariants. Here, by using three simple deformations (pure dilatation, pure shear and uniaxial stretch), we demonstrate that the compressible HGO-C formulation does not correctly model compressible anisotropic material behaviour, because the anisotropic component of the model is insensitive to volumetric deformation due to the use of isochoric anisotropic invariants. In order to correctly model compressible anisotropic behaviour we present a modified anisotropic (MA) model, whereby the full anisotropic invariants are used, so that a volumetric anisotropic contribution is represented. The MA model correctly predicts an anisotropic response to hydrostatic tensile loading, whereby a sphere deforms into an ellipsoid. It also computes the correct anisotropic stress state for pure shear and uniaxial deformation. To look at more practical applications, we developed a finite element user-defined material subroutine for the simulation of stent deployment in a slightly compressible artery. Significantly higher stress triaxiality and arterial compliance are computed when the full anisotropic invariants are used (MA model) instead of the isochoric form (HGO-C model).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An ultrasonic measurement of stress in steel without calibration: the angled shear wave identity

Measuring stress levels in loaded structures is crucial to assess and monitor their health, and to predict the length of their remaining structural life. However, measuring stress non-destructively has proved quite challenging. Many ultrasonic methods are able to accurately predict in-plane stresses in a controlled laboratory environment, but struggle to be robust outside, in a real world setting. That is because they rely either on knowing beforehand the material constants (which are difficult to acquire) or they require significant calibration for each specimen. Here we present a simple ultrasonic method to evaluate the in-plane stress in situ directly, without knowing any material constants. This method only requires measuring the speed of two angled shear waves. It is based on a formula which is exact for incompressible solids, such as soft gels or tissues, and is approximately true for compressible "hard" solids, such as steel and other metals. We validate the formula against virtual experiments using Finite Element simulations, and find it displays excellent accuracy, with a small error of the order of 1%.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

At least three invariants are necessary to model the mechanical response of incompressible, transversely isotropic materials

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Counter-intuitive results in acousto-elasticity

We present examples of body wave and surface wave propagation in deformed solids where the slowest and the fastest waves do not travel along the directions of least and greatest stretch, respectively. These results run counter to commonly accepted theory, practice, and implementation of the principles of acousto-elasticity in initially isotropic solids. For instance we find that in nickel and steel, the fastest waves are along the direction of greatest compression, not greatest extension (and vice-versa for the slowest waves), as soon as those solids are deformed. Further, we find that when some materials are subject to a small-but-finite deformations, other extrema of wave speeds appear in non-principal directions. Examples include nickel, steel, polystyrene, and a certain hydrogel. The existence of these ``oblique'', non-principal extremal waves complicates the protocols for the non-destructive determination of the directions of extreme strains.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deficiencies in numerical models of anisotropic nonlinearly elastic materials

Incompressible nonlinearly hyperelastic materials are rarely simulated in Finite Element numerical experiments as being perfectly incompressible because of the numerical difficulties associated with globally satisfying this constraint. Most commercial Finite Element packages therefore assume that the material is slightly compressible. It is then further assumed that the corresponding strain-energy function can be decomposed additively into volumetric and deviatoric parts. We show that this decomposition is not physically realistic, especially for anisotropic materials, which are of particular interest for simulating the mechanical response of biological soft tissue. The most striking illustration of the shortcoming is that with this decomposition, an anisotropic cube under hydrostatic tension deforms into another cube instead of a hexahedron with non-parallel faces. Furthermore, commercial numerical codes require the specification of a `compressibility parameter' (or `penalty factor'), which arises naturally from the flawed additive decomposition of the strain-energy function. This parameter is often linked to a `bulk modulus', although this notion makes no sense for anisotropic solids; we show that it is essentially an arbitrary parameter and that infinitesimal changes to it result in significant changes in the predicted stress response. This is illustrated with numerical simulations for biaxial tension experiments of arteries, where the magnitude of the stress response is found to change by several orders of magnitude when infinitesimal changes in `Poisson's ratio' close to the perfect incompressibility limit of 1/2 are made.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dominant negative Poynting effect for soft tissues

We identify three distinct shearing modes for simple shear deformations of transversely isotropic soft tissue which allow for both positive and negative Poynting effects (that is, they require compressive and tensile lateral normal stresses, respectively, in order to maintain simple shear). The positive Poynting effect is that usually found for isotropic rubber. Here, specialisation of the general results to three strain-energy functions which are quadratic in the anisotropic invariants, linear in the isotropic strain invariants and consistent with the linear theory suggests that there are two Poynting effects which can accompany the shearing of soft tissue: a dominant negative effect in one mode of shear and a relatively small positive effect in the other two modes. We propose that the relative inextensibility of the fibres relative to the matrix is the primary mechanism behind this large negative Poynting effect.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effects of nerve bundle geometry on neurotrauma evaluation

Objective: We confirm that alteration of a neuron structure can induce abnormalities in signal propagation for nervous systems, as observed in brain damage. Here, we investigate the effects of geometrical changes and damage of a neuron structure in 2 scaled nerve bundle models, made of myelinated nerve fibers or unmyelinated nerve fibers. Methods: We propose a 3D finite element model of nerve bundles, combining a real-time full electromechanical coupling, a modulated threshold for spiking activation, and independent alteration of the electrical properties for each fiber. We then simulate mechanical compression and tension to induce damage at the membrane of a nerve bundle made of 4 fibers. We examine the resulting changes in strain and neural activity by considering in turn the cases of intact and traumatized nerve membranes. Results: Our results show lower strain and lower electrophysiological impairments in unmyelinated fibers than in myelinated fibers, higher deformation levels in larger bundles, and higher electrophysiological impairments in smaller bundles. Conclusion: We conclude that the insulation sheath of myelin constricts the membrane deformation and scatters plastic strains within the bundle, that larger bundles deform more than small bundles, and that small fibers tolerate a higher level of elongation before mechanical failure.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Elastic Cherenkov effects in transversely isotropic soft materials-I: Theoretical analysis, simulations and inverse method

A body force concentrated at a point and moving at a high speed can induce shear-wave Mach cones in dusty-plasma crystals or soft materials, as observed experimentally and named the elastic Cherenkov effect (ECE). The ECE in soft materials forms the basis of the supersonic shear imaging (SSI) technique, an ultrasound-based dynamic elastography method applied in clinics in recent years. Previous studies on the ECE in soft materials have focused on isotropic material models. In this paper, we investigate the existence and key features of the ECE in anisotropic soft media, by using both theoretical analysis and finite element (FE) simulations, and we apply the results to the non=invasive and non-destructive characterization of biological soft tissues. We also theoretically study the characteristics of the shear waves induced in a deformed hyperelastic anisotropic soft material by a source moving with high speed, considering that contact between the ultrasound probe and the soft tissue may lead to finite deformation. On the basis of our theoretical analysis and numerical simulations, we propose an inverse approach to infer both the anisotropic and hyperelastic parameters of incompressible transversely isotropic (TI) soft materials. Finally, we investigate the properties of the solutions to the inverse problem by deriving the condition numbers in analytical form and performing numerical experiments. In Part II of the paper, both ex vivo and in vivo experiments are conducted to demonstrate the applicability of the inverse method in practical use.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electro-mechanical response of a 3D nerve bundle model to mechanical loads leading to axonal injury

Objective: Traumatic brain injuries and damage are major causes of death and disability. We propose a 3D fully coupled electro-mechanical model of a nerve bundle to investigate the electrophysiological impairments due to trauma at the cellular level. Methods: The coupling is based on a thermal analogy of the neural electrical activity by using the finite element software Abaqus CAE 6.13-3. The model includes a real-time coupling, modulated threshold for spiking activation and independent alteration of the electrical properties for each 3-layer fibre within a nerve bundle as a function of strain. Results: Results of the coupled electro-mechanical model are validated with previously published experimental results of damaged axons. Here, the cases of compression and tension are simulated to induce (mild, moderate and severe) damage at the nerve membrane of a nerve bundle, made of four fibres. Changes in strain, stress distribution, and neural activity are investigated for myelinated and unmyelinated nerve fibres, by considering the cases of an intact and of a traumatized nerve membrane. Conclusion: A fully coupled electro-mechanical modelling approach is established to provide insights into crucial aspects of neural activity at the cellular level due to traumatic brain injury. Significance: One of the key findings is the 3D distribution of residual stresses and strains at the membrane of each fibre due to mechanically induced electrophysiological impairments, and its impact on signal transmission.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electro-mechanically guided growth and patterns

Several experiments have demonstrated the existence of an electro-mechanical effect in many biological tissues and hydrogels, and its actual influence on growth, migration, and pattern formation. Here, to model these interactions and capture some growth phenomena found in Nature, we extend volume growth theory to account for an electro-elasticity coupling. Based on the multiplicative decomposition, we present a general analysis of isotropic growth and pattern formation of electro-elastic solids under external mechanical and electrical fields. As an example, we treat the case of a tubular structure to illustrate an electro-mechanically guided growth affected by axial strain and radial voltage. Our numerical results show that a high voltage can enhance the non-uniformity of the residual stress distribution and induce extensional buckling, while a low voltage can delay the onset of wrinkling shapes and can also generate more complex morphologies. Within a controllable range, axial tensile stretching shows the ability to stabilise the tube and help form more complex 3D patterns, while compressive stretching promotes instability. Both the applied voltage and external axial strain have a significant impact on guiding growth and pattern formation. Our modelling provides a basic tool for analysing the growth of electro-elastic materials, which can be useful for designing a pattern prescription strategy or growth self-assembly in Engineering.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electrostatically tunable axisymmetric vibrations of soft electro-active tubes

Due to their unique electromechanical coupling properties, soft electro-active (SEA) resonators are actively tunable, extremely suitable, and practically important for designing the next-generation acoustic and vibration treatment devices. In this paper, we investigate the electrostatically tunable axisymmetric vibrations of SEA tubes with different geometric sizes. We consider both axisymmetric torsional and longitudinal vibrations for an incompressible SEA cylindrical tube under inhomogeneous biasing fields induced by radial electric voltage and axial pre-stretch. We then employ the state-space method, which combines the state-space formalism in cylindrical coordinates with the approximate laminate technique, to derive the frequency equations for two separate classes of axisymmetric vibration of the tube subjected to appropriate boundary conditions. We perform numerical calculations to validate the convergence and accuracy of the state-space method and to illuminate that the axisymmetric vibration characteristics of SEA tubes may be tuned significantly by adjusting the electromechanical biasing fields as well as altering the tube geometry. The reported results provide a solid guidance for the proper design of tunable resonant devices composed of SEA tubes

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electrothermal equivalent three-dimensional Finite Element Model of a single neuron

Objective: We propose a novel approach for modelling the inter-dependence of electrical and mechanical phenomena in nervous cells, by using electro-thermal equivalences in finite element (FE) analysis so that existing thermo-mechanical tools can be applied. Methods: First, the equivalence between electrical and thermal properties of the nerve materials is established, and results of a pure heat conduction analysis performed in Abaqus CAE Software 6.13-3 are validated with analytical solutions for a range of steady and transient conditions. This validation includes the definition of equivalent active membrane properties that enable prediction of the action potential. Then, as a step towards fully coupled models, electromechanical coupling is implemented through the definition of equivalent piezoelectric properties of the nerve membrane using the thermal expansion coefficient, enabling prediction of the mechanical response of the nerve to the action potential. Results: Results of the coupled electro-mechanical model are validated with previously published experimental results of deformation for the squid giant axon, crab nerve fibre and garfish olfactory nerve fibre. Conclusion: A simplified coupled electro-mechanical modelling approach is established through an electro-thermal equivalent FE model of a nervous cell for biomedical applications. Significance: One of the key findings is the mechanical characterization of the neural activity in a coupled electro-mechanical domain, which provides insights into the electromechanical behaviour of nervous cells, such as thinning of the membrane. This is a first step towards modelling 3D electromechanical alteration induced by trauma at nerve bundle, tissue and organ levels.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Extreme softness of brain matter in simple shear

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gent models for the inflation of spherical balloons

We revisit an iconic deformation of non-linear elasticity: the inflation of a rubber spherical thin shell. We use the 3-parameter Mooney and Gent-Gent (GG) phenomenological models to explain the stretch-strain curve of a typical inflation, as these two models cover a wide spectrum of known models for rubber, including the Varga, Mooney-Rivlin, one-term Ogden, Gent-Thomas and Gent models. We find that the basic physics of inflation exclude the Varga, one-term Ogden and Gent-Thomas models. We find the link between the exact solution of non-linear elasticity and the membrane and Young-Laplace theories often used a priori in the literature. We compare the performance of both models on fitting the data for experiments on rubber balloons and animal bladder. We conclude that the GG model is the most accurate and versatile model on offer for the modelling of rubber balloon inflation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Guided waves in pre-stressed hyperelastic plates and tubes: Application to the ultrasound elastography of thin-walled soft materials

In vivo measurement of the mechanical properties of thin-walled soft tissues (e.g., mitral valve, artery and bladder) and in situ mechanical characterization of thin-walled artificial soft biomaterials in service are of great challenge and difficult to address via commonly used testing methods. Here we investigate the properties of guided waves generated by focused acoustic radiation force in immersed pre-stressed plates and tubes, and show that they can address this challenge. To this end, we carry out both (i) a theoretical analysis based on incremental wave motion in finite deformation theory and (ii) finite element simulations. Our analysis leads to a novel method based on the ultrasound elastography to image the elastic properties of pre-stressed thin-walled soft tissues and artificial soft materials in a non-destructive and non-invasive manner. To validate the theoretical and numerical solutions and demonstrate the usefulness of the corresponding method in practical measurements, we perform (iii) experiments on polyvinyl alcohol cryogel phantoms immersed in water, using the Verasonics V1 System equipped with a L10-5 transducer. Finally, potential clinical applications of the method have been discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Head-to-nerve analysis of electromechanical impairments of diffuse axonal injury

The aim was to investigate mechanical and functional failure of diffuse axonal injury (DAI) in nerve bundles following frontal head impacts, by finite element simulations. Anatomical changes following traumatic brain injury are simulated at the macroscale by using a 3D head model. Frontal head impacts at speeds of 2.5-7.5 m/s induce mild-to-moderate DAI in the white matter of the brain. Investigation of the changes in induced electromechanical responses at the cellular level is carried out in two scaled nerve bundle models, one with myelinated nerve fibres, the other with unmyelinated nerve fibres. DAI occurrence is simulated by using a real-time fully coupled electromechanical framework, which combines a modulated threshold for spiking activation and independent alteration of the electrical properties for each three-layer fibre in the nerve bundle models. The magnitudes of simulated strains in the white matter of the brain model are used to determine the displacement boundary conditions in elongation simulations using the 3D nerve bundle models. At high impact speed, mechanical failure occurs at lower strain values in large unmyelinated bundles than in myelinated bundles or small unmyelinated bundles; signal propagation continues in large myelinated bundles during and after loading, although there is a large shift in baseline voltage during loading; a linear relationship is observed between the generated plastic strain in the nerve bundle models and the impact speed and nominal strains of the head model. The myelin layer protects the fibre from mechanical damage, preserving its functionalities.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Influence of initial residual stress on growth and pattern creation for a layered aorta

Residual stress is ubiquitous and indispensable in most biological and artificial materials, where it sustains and optimizes many biological and functional mechanisms. The theory of volume growth, starting from a stress-free initial state, is widely used to explain the creation and evolution of growth-induced residual stress and the resulting changes in shape, and to model how growing bio-tissues such as arteries and solid tumors develop a strategy of pattern creation according to geometrical and material parameters. This modelling provides promising avenues for designing and directing some appropriate morphology of a given tissue or organ and achieve some targeted biomedical function. In this paper, we rely on a modified, augmented theory to reveal how we can obtain growth-induced residual stress and pattern evolution of a layered artery by starting from an existing, non-zero initial residual stress state. We use experimentally determined residual stress distributions of aged bi-layered human aortas and quantify their influence by a magnitude factor. Our results show that initial residual stress has a more significant impact on residual stress accumulation and the subsequent evolution of patterns than geometry and material parameters. Additionally, we provide an essential explanation for growth-induced patterns driven by differential growth coupled to an initial residual stress. Finally, we show that initial residual stress is a readily available way to control growth-induced pattern creation for tissues and thus may provide a promising inspiration for biomedical engineering.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring the linear and nonlinear elastic properties of brain tissue with shear waves and inverse analysis

We use supersonic shear wave imaging (SSI) technique to measure not only the linear but also the nonlinear elastic properties of brain matter. Here, we tested six porcine brains ex vivo and measured the velocities of the plane shear waves induced by acoustic radiation force at different states of pre-deformation when the ultrasonic probe is pushed into the soft tissue. We relied on an inverse method based on the theory governing the propagation of small-amplitude acoustic waves in deformed solids to interpret the experimental data. We found that, depending on the subjects, the resulting initial shear modulus mu0 varies from 1.8 to 3.2 kPa, the stiffening parameter b of the hyperelastic Demiray-Fung model from 0.13 to 0.73, and the third- (A) and fourth-order (D) constants of weakly nonlinear elasticity from -1.3 to -20.6 kPa and from 3.1 to 8.7 kPa, respectively. Paired t-test performed on the experimental results of the left and right lobes of the brain shows no significant difference. These values are in line with those reported in the literature on brain tissue, indicating that the SSI method, combined to the inverse analysis, is an efficient and powerful tool for the mechanical characterization of brain tissue, which is of great importance for computer simulation of traumatic brain injury and virtual neurosurgery.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mechanical characterization of brain tissue in simple shear at dynamic strain rates

During severe impact conditions, brain tissue experiences a rapid and complex deformation, which can be seen as a mixture of compression, tension and shear. Moreover, diffuse axonal injury (DAI) occurs in animals and humans when both the strains and strain rates exceed 10% and 10/s, respectively. Knowing the mechanical properties of brain tissue in shear at these strains and strain rates is thus of particular importance, as they can be used in finite element simulations to predict the occurrence of brain injuries under different impact conditions. In this research, an experimental setup was developed to perform simple shear tests on porcine brain tissue at strain rates < 120/s. The maximum measured shear stress at strain rates of 30, 60, 90 and 120/s was 1.15 +/- 0.25 kPa, 1.34 +/- 0.19 kPa, 2.19 +/- 0.225 kPa and 2.52 +/- 0.27 kPa, (mean +/- SD), respectively, at the maximum amount of shear, K = 1. Good agreement of experimental, theoretical (Ogden and MooneyRivlin models) and numerical shear stresses was achieved (p = 0.7866 - 0.9935). Specimen thickness effects (2.0 - 10.0 mm thick specimens) were also analyzed numerically and we found that there is no significant difference (p = 0.9954) in the shear stress magnitudes, indicating a homogeneous deformation of the specimens during simple shear tests. Stress relaxation tests in simple shear were also conducted at different strain magnitudes (10% - 60% strain) with the average rise time of 14 ms. This allowed us to estimate elastic and viscoelastic parameters (mu = 4942.0 Pa and Prony parameters: g1 = 0.520, g2 = 0.3057, tau1 = 0.0264 s, tau2 = 0.011 s) that can be used in FE software to analyze the hyperviscoelastic behavior of brain tissue.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mechanical characterization of brain tissue in tension at dynamic strain rates

Mechanical characterization of brain tissue at high loading velocities is crucial for modeling Traumatic Brain Injury (TBI). During severe impact conditions, brain tissue experiences compression, tension and shear. Limited experimental data is available for brain tissue in extension at dynamic strain rates. In this research, a High Rate Tension Device (HRTD) was developed to obtain dynamic properties of brain tissue in extension at strain rates of < 90/s. In vitro tensile tests were performed to obtain properties of brain tissue at strain rates of 30, 60 and 90/s up to 30% strain. The brain tissue showed a stiffer response with increasing strain rates, showing that hyperelastic models are not adequate. Specifically, the tensile engineering stress at 30% strain was 3.1 +/- 0.49 kPa, 4.3 +/- 0.86 kPa, 6.5 +/- 0.76 kPa (mean +/- SD) at strain rates of 30, 60 and 90/s, respectively. Force relaxation tests in tension were also conducted at different strain magnitudes (10-60% strain) with the average rise time of 24 ms, which were used to derive time dependent parameters. One-term Ogden, Fung and Gent models were used to obtain material parameters from the experimental data. Numerical simulations were performed using a one-term Ogden model to analyze hyperelastic behavior of brain tissue up to 30% strain. The material parameters obtained in this study will help to develop biofidelic human brain finite element models, which can subsequently be used to predict brain injuries under impact conditions and as a reconstruction and simulation tool for forensic investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modified multiplicative decomposition model for tissue growth: Beyond the initial stress-free state

The multiplicative decomposition model is widely employed for predicting residual stresses and morphologies of biological tissues due to growth. However, it relies on the assumption that the tissue is initially in a stress-free state, which conflicts with the observations that any growth state of tissue is under a significant level of residual stresses that helps to maintain its ideal mechanical conditions. Here, we propose a modified multiplicative decomposition model where the initial state (or reference configuration) of biological tissue is endowed with residual stress instead of being stress-free. Releasing theoretically the initial residual stress, the initially stressed state is first transmitted into a virtual stress-free state, resulting in an initial elastic deformation. The initial virtual stress-free state subsequently grows to another counterpart with a growth deformation and the latter is further integrated into its natural configuration with an excessive elastic deformation that ensures tissue compatibility. With this decomposition, the total deformation may be expressed as the product of elastic deformation, growth deformation and initial elastic deformation, while the corresponding free energy density depends on the initial residual stress and the total deformation. We address three key issues: explicit expression of the free energy density,predetermination of the initial elastic deformation, and initial residual stress. Finally, we consider a tubular organ to demonstrate the effects of the proposed initial residual stress on stress distribution and on shape formation through an incremental stability analysis. Our results suggest that the initial residual stress exerts a major influence on the growth stress and the morphology of tissues.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Morphology of residually-stressed tubular tissues: Beyond the elastic multiplicative decomposition

Many interesting shapes appearing in the biological world are formed by the onset of mechanical instability. In this work we consider how the build-up of residual stress can cause a solid to buckle. In all past studies a fictitious (virtual) stress-free state was required to calculate the residual stress. In contrast, we use a model which is simple and allows the prescription of any residual stress field. We specialize the analysis to an elastic tube subject to a two-dimensional residual stress, and find that incremental wrinkles can appear on its inner or its outer face, depending on the location of the highest value of the residual hoop stress. We further validate the predictions of the incremental theory with finite element simulations, which allow us to go beyond this threshold and predict the shape, number and amplitude of the resulting creases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On anisotropic elasticity and questions concerning its Finite Element implementation

We give conditions on the strain-energy function of nonlinear anisotropic hyperelastic materials that ensure compatibility with the classical linear theories of anisotropic elasticity. We uncover the limitations associated with the volumetric deviatoric separation of the strain energy used, for example, in many Finite Element (FE) codes in that it does not fully represent the behavior of anisotropic materials in the linear regime. This limitation has important consequences. We show that, in the small deformation regime, a FE code based on the volumetric-deviatoric separation assumption predicts that a sphere made of a compressible anisotropic material deforms into another sphere under hydrostatic pressure loading, instead of the expected ellipsoid. For finite deformations, the commonly adopted assumption that fibres cannot support compression is incorrectly implemented in current FE codes and leads to the unphysical result that under hydrostatic tension a sphere of compressible anisotropic material deforms into a larger sphere.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Prescribing patterns in growing tubular soft matter by initial residual stress

Initial residual stress is omnipresent in biological tissues and soft matter, and can affect growth-induced pattern selection significantly. Here we demonstrate this effect experimentally by letting soft tubes grow in the presence or absence of initial residual stress and by observing different growth pattern evolutions. These experiments motivate us to model the mechanisms at play when a growing bilayer tubular organ spontaneously displays buckling patterns on its inner surface. We demonstrate that not only differential growth, geometry and elasticity, but also initial residual stress distribution, exert a notable influence on these pattern phenomena. Prescribing an initial residual stress distribution offers an alternative or a more effective way to implement pattern selection for growable bio-tissues or soft matter. The results also show promise for the design of 4D bio-mimic printing protocols or for controlling hydrogel actuators.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stability analysis of charge-controlled soft dielectric plates

We examine the stability of a soft dielectric plate deformed by the coupled effects of a mechanical pre-stress applied on its lateral faces and an electric field applied through its thickness under charge control. The electric field is created by spraying charges on the major faces of the plate: although in practice this mode of actuation is harder to achieve than a voltage-driven deformation, here we find that it turns out to be much more stable in theory and in simulations. First we show that the electromechanical instability based on the Hessian criterion associated with the free energy of the system does not occur at all for charge-driven dielectrics for which the electric displacement is linear in the electric field. Then we show that the geometric instability associated with the formation of small-amplitude wrinkles on the faces of the plate that arises under voltage control does not occur either under charge control. This is in complete contrast to voltage-control actuation, where Hessian and wrinkling instabilities can occur once certain critical voltages are reached. For the mechanical pre-stresses, two modes that can be implemented in practice are used: equi-biaxial and uni-axial. We confirm the analytical and numerical stability results of homogeneous deformation modes with Finite Element simulations of real actuations, where inhomogeneous fields may develop. We find complete agreement in the equi-biaxial case, and very close agreement in the uni-axial case, when the pre-stress is due to a dead-load weight. In the latter case, the simulations show that small inhomogeneous effects develop near the clamps, and eventually a compressive lateral stress emerges, leading to a breakdown of the numerics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Straightening wrinkles

We consider the elastic deformation of a circular cylindrical sector composed of an incompressible isotropic soft solid when it is straightened into a rectangular block. In this process, the circumferential line elements on the original inner face of the sector are stretched while those on the original outer face are contracted. We investigate the geometrical and physical conditions under which the latter line elements can be contracted to the point where a localized incremental instability develops. We provide a robust algorithm to solve the corresponding two-point boundary value problem, which is stiff numerically. We illustrate the results with full incremental displacement fields in the case of Mooney-Rivlin materials and also perform an asymptotic analysis for thin sectors

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Straightening: Existence, uniqueness and stability

One of the least studied universal deformations of incompressible nonlinear elasticity, namely the straightening of a sector of a circular cylinder into a rectangular block, is revisited here and, in particular, issues of existence and stability are addressed. Particular attention is paid to the system of forces required to sustain the large static deformation, including by the application of end couples. The influence of geometric parameters and constitutive models on the appearance of wrinkles on the compressed face of the block is also studied. Different numerical methods for solving the incremental stability problem are compared and it is found that the impedance matrix method, based on the resolution of a matrix Riccati differential equation, is the more precise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Strain energy function for isotropic non-linear elastic incompressible solids with linear finite strain response in shear and torsion

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The effect of an exterior electric field on the instability of dielectric plates

We investigate the theoretical nonlinear response, Hessian stability, and possible wrinkling behaviour of a voltage-activated dielectric plate immersed in a tank filled with silicone oil. Fixed rigid electrodes are placed on the top and bottom of the tank, and an electric field is generated by a potential difference between the electrodes. We solve the associated incremental boundary value problem of superimposed, inhomogeneous small-amplitude wrinkles, signalling the onset of instability. We decouple the resulting bifurcation equation into symmetric and antisymmetric modes. For a neo-Hookean dielectric plate, we show that a potential difference between the electrodes can induce a thinning of the plate and thus an increase of its planar area, similar to the scenarios encountered when there is no silicone oil. However, we also find that, depending on the material and geometric parameters, an increasing applied voltage can also lead to a thickening of the plate, and thus a shrinking of its area. In that scenario, Hessian instability and wrinkling bifurcation may then occur spontaneously once some critical voltages are reached.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Torsion instability of soft solid cylinders

The application of pure torsion to a long and thin cylindrical rod is known to provoke a twisting instability, evolving from an initial kink to a knot. In the torsional parallel-plate rheometry of stubby cylinders, the geometrical constraints impose zero displacement of the axis of the cylinder, preventing the occurrence of such twisting instability. Under these experimental conditions, wrinkles occur on the cylinder's surface at a given critical angle of torsion. Here we investigate this subclass of elastic instability--which we call torsion instability--of soft cylinders subject to a combined finite axial stretch and torsion, by applying the theory of incremental elastic deformation superimposed on finite strains. We formulate the incremental boundary elastic problem in the Stroh differential form, and use the surface impedance method to build a robust numerical procedure for deriving the marginal stability curves. We present the results for a Mooney-Rivlin material and study the influence of the material parameters on the elastic bifurcation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tuning the pull-in instability of soft dielectric elastomers through loading protocols

Pull-in (or electro-mechanical) instability occurs when a drastic decrease in the thickness of a dielectric elastomer results in electrical breakdown, which limits the applications of dielectric devices. Here we derive the criteria for determining the pull-in instability of dielectrics actuated by different loading methods: voltage-control, charge-control, fixed pre-stress and fixed pre-stretch, by analyzing the free energy of the actuated systems. The Hessian criterion identifies a maximum in the loading curve beyond which the elastomer will stretch rapidly and lose stability, and can be seen as a path to failure. We present numerical calculations for neo-Hookean ideal dielectrics, and obtain the maximum allowable actuation stretch of a dielectric before failure by electrical breakdown. We find that applying a fixed pre-stress or a fixed pre-stretch to a charge-driven dielectric may decrease the stretchability of the elastomer, a scenario which is the opposite of what happens in the case of a voltage-driven dielectric. Results show that a reversible large actuation of a dielectric elastomer, free of the pull-in instability, can be achieved by tuning the actuation method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Wrinkles in the opening angle method

We investigate the stability of the deformation modeled by the opening angle method, often used to give a measure of residual stresses in arteries and other biological soft tubular structures. Specifically, we study the influence of stiffness contrast, dimensions and inner pressure on the onset of wrinkles when an open sector of a soft tube, coated with a stiffer film, is bent into a full cylinder. The tube and its coating are made of isotropic, incompressible, hyperelastic materials. We provide a full analytical exposition of the governing equations and the associated boundary value problem for the large deformation and for the superimposed small-amplitude wrinkles. For illustration, we solve them numerically with a robust algorithm in the case of Mooney-Rivlin materials. We confront the results to experimental data that we collected for soft silicone sectors. We study the influence of axial stretch and inner pressure on the stability of closed-up coated tubes with material parameters comparable with those of soft biological tubes such as arteries and veins, although we do not account for anisotropy. We find that the large deformation described in the opening angle method does not always exist, as it can become unstable for certain combinations of dimensions and material parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Poynting effect of brain matter in torsion

We investigate experimentally and model theoretically the mechanical behaviour of brain matter in torsion. Using a strain-controlled rheometer we perform torsion tests on fresh porcine brain samples. We quantify the torque and the normal force required to twist a cylindrical sample at constant twist rate. Data fitting gives a mean value for the shear modulus μ = 900 $\pm$ 312 Pa and for the second Mooney-Rivlin parameter $c_2$ = 297 $\pm$ 189 Pa, indicative of extreme softness. Our results show that brain always displays a positive Poynting effect; in other words, it expands in the direction perpendicular to the plane of twisting. We validate the experiments with Finite Element simulations and show that when a human head experiences a twisting motion in the horizontal plane, the brain can experience large forces in the axial direction.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

The mechanics of human brain organoids

Organoids are prototypes of human organs derived from cultured human stem cells. They provide a reliable and accurate experimental model to study the physical mechanisms underlying the early developmental stages of human organs morphogenesis and, in particular, the early morphogenesis of the cortex. Here, we propose a mathematical model to elucidate the role played by two mechanisms which have been experimentally proven to be crucial in shaping human brain organoids: the contraction of the inner core of the organoid and the microstructural remodeling of the outer cortex. Our results show that both mechanisms are crucial for the final shape of the organoid and can explain the origin of brain pathologies such as lissencephaly (smooth brain).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Non-invasive evaluation of skin tension lines with elastic waves

Background: Since their discovery by Karl Langer in the 19th Century, Skin Tension Lines (STLs) have been used by surgeons to decide the location and orientation of an incision. Although these lines are patient-specific, most surgeons rely on generic maps to determine their orientation. Beyond the imprecise pinch test, there remains no accepted method for determining STLs in vivo. Methods: (i) The speed of an elastic motion travelling radially on the skin of canine cadavers was measured with a commercial device called the Reviscometer. (ii) Similar to the original experiments conducted by Karl Langer, circular excisions were made on the skin and the geometric changes to the resulting wounds and excised samples were used to determine the orientation of STLs. Results: A marked anisotropy in the speed of the elastic wave travelling radially was observed. The orientation of the fastest wave was found to correlate with the orientation of the elongated wound (P<0.001, R^2 = 74%). Similarly, the orientation of fastest wave was the same for both in vivo and excised isolated samples, indicating that the STLs have a structural basis. Resulting wounds expanded by an average area of 9% (+16% along STL and -10% across) while excised skin shrunk by an average area of 33% (23% along STL and 10% across). Conclusion: Elastic surface wave propagation has been validated experimentally as a robust method for determining the orientation of STLs non-destructively and non-invasively. This study has implications for the identification of STLs and for the prediction of skin tension levels, both important factors in both human and veterinary reconstructive surgery.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Initial Stress Symmetry and Applications in Elasticity

An initial stress within a solid can arise to support external loads or from processes such as thermal expansion in inert matter or growth and remodelling in living materials. For this reason it is useful to develop a mechanical framework of initially stressed solids irrespective of how this stress formed. An ideal way to do this is to write the free energy density $Ψ= Ψ(\boldsymbol F, \boldsymbol τ)$ in terms of initial stress $\boldsymbol τ$ and the elastic deformation gradient $\boldsymbol F$. In this paper we present a new constitutive condition for initially stressed materials, which we call the initial stress symmetry (ISS). We focus on two consequences of this symmetry. First we examine how ISS restricts the free energy density $Ψ= Ψ(\boldsymbol F, \boldsymbol τ) $ and present two examples of $Ψ(\boldsymbol F, \boldsymbol τ)$ that satisfy ISS. Second we show that the initial stress can be derived from the Cauchy stress and the elastic deformation gradient. To illustrate we take an example from biomechanics and calculate the optimal Cauchy stress within an artery subjected to internal pressure. We then use ISS to derive the optimal target residual stress for the material to achieve after remodelling.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A High Rate Tension Device for Characterizing Brain Tissue

The mechanical characterization of brain tissue at high loading velocities is vital for understanding and modeling Traumatic Brain Injury (TBI). The most severe form of TBI is diffuse axonal injury (DAI) which involves damage to individual nerve cells (neurons). DAI in animals and humans occurs at strains > 10% and strain rates > 10/s. The mechanical properties of brain tissues at these strains and strain rates are of particular significance, as they can be used in finite element human head models to accurately predict brain injuries under different impact conditions. Existing conventional tensile testing machines can only achieve maximum loading velocities of 500 mm/min, whereas the Kolsky bar apparatus is more suitable for strain rates > 100/s. In this study, a custom-designed high rate tension device is developed and calibrated to estimate the mechanical properties of brain tissue in tension at strain rates < 90/s, while maintaining a uniform velocity. The range of strain can also be extended to 100% depending on the thickness of a sample. The same apparatus can be used to characterize the dynamic behavior of skin and other soft biological tissues by using appropriately sized load cells with a capacity of 10 N and above.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bending instabilities of soft tissues

Rubber components and soft tissues are often subjected to large bending deformations "in service". The circumferential line elements on the inner face of a bent block can contract up to a certain critical stretch ratio lambda_cr (say) before bifurcation occurs and axial creases appear. For several models used to describe rubber, it is found that lambda_cr = 0.56, allowing for a 44% contraction. For models used to describe arteries it is found, somewhat surprisingly, that the strain-stiffening effect promotes instability. For example, the models used for the artery of a 70 year old human predict that lambda_cr = 0.73, allowing only for a 27% contraction. Tensile experiments conducted on pig skin indicate that bending instabilities should occur even earlier there.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bimodular rubber buckles early in bending

A block of rubber eventually buckles under severe flexure, and several axial wrinkles appear on the inner curved face of the bent block. Experimental measurements reveal that the buckling occurs earlier ---at lower compressive strains--- than expected from theoretical predictions. This paper shows that if rubber is modeled as being bimodular, and specifically, as being stiffer in compression than in tension, then flexure bifurcation happens indeed at lower levels of compressive strain than predicted by previous investigations (these included taking into account finite size effects, compressibility effects, and strain-stiffening effects.) Here the effect of bimodularity is investigated within the theory of incremental buckling, and bifurcation equations, numerical methods, dispersion curves, and field variations are presented and discussed. It is also seen that Finite Element Analysis software seems to be unable to encompass in a realistic manner the phenomenon of bending instability for rubber blocks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Characterising the Anisotropic Mechanical Properties of Excised Human Skin

The mechanical properties of skin are important for a number of applications including surgery, dermatology, impact biomechanics and forensic science. In this study we have investigated the influence of location and orientation on the deformation characteristics of 56 samples of excised human skin. Uniaxial tensile tests were carried out at a strain rate of 0.012s$^{-1}$ on skin from the back. Digital Image Correlation was used for 2D strain measurement and a histological examination of the dermis was also performed. The mean ultimate tensile strength (UTS) was 21.6$\pm$8.4MPa, the mean failure strain 54$\pm$17%, the mean initial slope 1.18$\pm$0.88MPa, the mean elastic modulus 83.3$\pm$34.9MPa and the mean strain energy was 3.6$\pm$1.6MJ/m$^3$. A multivariate analysis of variance has shown that these mechanical properties of skin are dependent upon the orientation of Langer lines (P$<$0.0001-P=0.046). The location of specimens on the back was also found to have a significant effect on the UTS (P =0.0002), the elastic modulus (P=0.001) and the strain energy (P=0.0052). The histological investigation concluded that there is a definite correlation between the orientation of Langer Lines and the preferred orientation of collagen fibres in the dermis (P$<$0.001). The data obtained in this study will provide essential information for those wishing to model the skin using a structural constitutive model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Circularly-polarized plane waves in a deformed Hadamard material

Small amplitude inhomogeneous plane waves propagating in any direction in a homogeneously deformed Hadamard material are considered. Conditions for circular polarization are established. The analysis relies on the use of complex vectors (or bivectors) to describe the slowness and the polarization of the waves. Generally, homogeneous circularly-polarized plane waves may propagate in only two directions, the directions of the acoustic axes, in a homogeneously deformed Hadamard material. For inhomogeneous circularly-polarized plane waves, the number of possibilities is far greater. They include an infinity of "transverse waves", as well as "longitudinal waves", and the superposition of "transverse waves" and "longitudinal waves", where "transverse" and "longitudinal" are used in the bivector sense. Each and every possibility of circular polarization is examined in turn, and explicit examples of solutions are given in every case.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Compression Instabilities of Tissues with Localized Strain Softening

The stress-strain relationship of biological soft tissues affected by Marfan's syndrome is believed to be non-convex. More specifically, Haughton and Merodio recently proposed a strain-energy density leading to localized strain softening, in order to model the unusual mechanical behavior of these isotropic, incompressible tissues. Here we investigate how this choice of strain energy affects the results of some instabilities studies, such as those concerned with the compression of infinite and semi-infinite solids, slabs, and cylinders, or with the bending of blocks, and draw comparisons with known results established previously for the case of a classical neo-Hookean solid. We find that the localized strain softening effect leads to early instability only when instability occurs at severe compression ratios for neo-Hookean solids, as is the case for bulk, surface, and bending instabilities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Determination of Friction Coefficient in Unconfined Compression of Brain Tissue

Unconfined compression tests are more convenient to perform on cylindrical samples of brain tissue than tensile tests in order to estimate mechanical properties of the brain tissue because they allow for homogeneous deformations. The reliability of these tests depends significantly on the amount of friction generated at the specimen/platen interface. Thus, there is a crucial need to find an approximate value of the friction coefficient in order to predict a possible overestimation of stresses during unconfined compression tests. In this study, a combined experimental-computational approach was adopted to estimate the dynamic friction coefficient mu of porcine brain matter against metal platens in compressive tests. Cylindrical samples of porcine brain tissue were tested up to 30% strain at variable strain rates, both under bonded and lubricated conditions in the same controlled environment. It was established that mu was equal to 0.09 +/- 0.03, 0.18 +/- 0.04, 0.18 +/- 0.04 and 0.20 +/- 0.02 at strain rates of 1, 30, 60 and 90/s, respectively. Additional tests were also performed to analyze brain tissue under lubricated and bonded conditions, with and without initial contact of the top platen with the brain tissue, with different specimen aspect ratios and with different lubricants (Phosphate Buffer Saline (PBS), Polytetrafluoroethylene (PTFE) and Silicon). The test conditions (lubricant used, biological tissue, loading velocity) adopted in this study were similar to the studies conducted by other research groups. This study will help to understand the amount of friction generated during unconfined compression of brain tissue for strain rates of up to 90/s.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Differences in tension and compression in the nonlinearly elastic bending of beams

The classical flexure problem of non-linear incompressible elasticity is revisited for elastic materials whose mechanical response is different in tension and compression---the so-called bimodular materials. The flexure problem is chosen to investigate this response since the two regions, one of tension and one of compression, can be identified easily using simple intuition. Two distinct problems are considered: the first is where the stress is assumed continuous across the boundary of the two regions, which assumption has a sound physical basis. The second problem considered is more speculative: it is where discontinuities of stress are allowed. It is shown that such discontinuities are necessarily small for many applications, but might nonetheless provide an explanation for the damage incurred by repeated flexure. Some experimental evidence of the possibility of bimodularity in elastomers is also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Elastic interface acoustic waves in twinned crystals

A new type of Interface Acoustic Waves (IAW) is presented, for single-crystal orthotropic twins bonded symmetrically along a plane containing only one common crystallographic axis. The effective boundary conditions show that the waves are linearly polarized at the interface, either transversally or longitudinally. Then the secular equation is obtained in full analytical form using new relationships for the displacement-traction quadrivector at the interface. For Gallium Arsenide and for Silicon, it is found that the IAWs with transverse (resp. longitudinal) polarization at the interface are of the Stoneley (resp. leaky) type.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Explicit secular equations for piezoacoustic surface waves: Rayleigh modes

The existence of a two-partial Rayleigh wave coupled to an electrical field in 2mm piezoelectric crystals is known but has rarely been investigated analytically. It turns out that the Z-cut, X-propagation problem can be fully solved, up to the derivation of the secular equation as a polynomial in the squared wave speed. For the metallized (unmetallized) boundary condition, the polynomial is of degree 10 (48). The relevant root is readily identified and the full description of the mechanical and electrical fields follows. The results are illustrated in the case of the superstrong piezoelectric crystal, Potassium niobate, for which the effective piezoelectric coupling coefficient is calculated to be about 0.1

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Explicit secular equations for piezoacoustic surface waves: Shear-Horizontal modes

Attention is given to surface waves of shear-horizontal modes in piezoelectric crystals permitting the decoupling between an elastic in-plane Rayleigh wave and a piezoacoustic anti-plane Bleustein-Gulyaev wave. Specifically, the crystals possess $\bar{4}$ symmetry (inclusive of $\bar{4}2$m, $\bar{4}3$m, and 23 classes) and the boundary is any plane containing the normal to a symmetry plane (rotated $Y$-cuts about the $Z$ axis). The secular equation is obtained explicitly as a polynomial not only for the metallized boundary condition but, in contrast to previous studies on the subject, also for other types of boundary conditions. For the metallized surface problem, the secular equation is a quadratic in the squared wave speed; for the un-metallized surface problem, it is a sextic in the squared wave speed; for the thin conducting boundary problem, it is of degree 16 in the speed. The relevant root of the secular equation can be identified and the complete solution is then found (attenuation factors, field profiles, etc.). The influences of the cut angle and of the conductance of the adjoining medium are illustrated numerically for GaAs ($\bar{4}3$m), BaLaGa$_3$O$_7$ ($\bar{4}2$m) and Bi$_{12}$GeO$_{20}$ (23). Indications are given on how to apply the method to crystals with 222 symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finite amplitude inhomogeneous waves in Mooney-Rivlin viscoelastic solids

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finite-amplitude inhomogeneous plane waves of exponential type in incompressible elastic materials

It is proved that elliptically-polarized finite-amplitude inhomogeneous plane waves may not propagate in an isotropic elastic material subject to the constraint of incompressibility. The waves considered are harmonic in time and exponentially attenuated in a direction distinct from the direction of propagation. The result holds whether the material is stress-free or homogeneously deformed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Influence of Preservation Temperature on the Measured Mechanical Properties of Brain Tissue

The large variability in experimentally measured mechanical properties of brain tissue is due to many factors including heterogeneity, anisotropy, age dependence and post-mortem time. Moreover, differences in test protocols also influence these measured properties. This paper shows that the temperature at which porcine brain tissue is stored or preserved prior to testing has a significant effect on the mechanical properties of brain tissue, even when tests are conducted at the same temperatures. Three groups of brain tissue were stored separately for at least one hour at three different preservation temperatures, i.e., ice cold, room temperature (22C) and body temperature (37C), prior to them all being tested at room temperature (approx. 22C). Significant differences in the corresponding initial elastic shear modulus mu (Pa) (at various amounts of shear, K, i.e., 0-0.2) were observed. The initial elastic moduli were 1043 +/- 271 Pa, 714 +/- 210 Pa and 497 +/- 156 Pa (mean +/- SD) at preservation temperatures of ice cold, 22C and 37C, respectively. Based on this investigation, it is strongly recommended that brain tissue samples must be preserved at an ice-cold temperature prior to testing in order to minimize the difference between the measured in vitro test results and the in vivo properties. A by-product of the study is that simple shear tests allow for large, almost perfectly homogeneous deformation of brain matter.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inhomogeneous "longitudinal" circularly-polarized plane waves in anisotropic elastic crystals

Conditions on the elastic stiffnesses of anisotropic crystals are derived such that circularly polarized longitudinal inhomogeneous plane waves with an isotropic slowness bivector may propagate for any given direction of the normal to the sagittal plane. Once this direction is chosen, then the wave speed, the direction of propagation, and the direction of attenuation are expressed in terms of the mass density, the elastic stiffnesses, and the angle between the normal to the sagittal plane and the normals (also called "optic axes") to the planes of central circular section of a certain ellipsoid. In the special case where this angle is zero, and in this special case only, such waves cannot propagate.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inhomogeneous "longitudinal" plane waves in a deformed elastic material

A homogeneous isotropic compressible Hadamard material has the property that an infinitesimal longitudinal homogeneous plane wave may propagate in every direction when the material is maintained in a state of arbitrary finite static homogeneous deformation. Here, as regards the wave, 'homogeneous' means that the direction of propagation of the wave is parallel to the direction of eventual attenuation and 'longitudinal' means that the wave is linearly polarized in a direction parallel to the direction of propagation. In other words, the displacement is of the form u = n cos k(n.x - ct), where n is a real vector. It is seen that the Hadamard material is the most general one for which a 'longitudinal' inhomogeneous plane wave may also propagate in any direction of a predeformed body. Here, 'inhomogeneous' means that the wave is attenuated in a direction distinct from the direction of propagation and 'longitudinal' means that the wave is elliptically polarized in the plane containing these two directions. In other words, the displacement is of the form u = Re {S exp[i (omega S.x - ct)]}, where S is a complex vector (or bivector). Then a Generalized Hadamard material is introduced. It is the most general homogeneous isotropic compressible material which allows the propagation of infinitesimal 'longitudinal' inhomogeneous plane circularly polarized waves for all choices of the isotropic directional bivector. Finally, the most general forms of response functions are found for homogeneously deformed isotropic elastic materials in which 'longitudinal' inhomogeneous plane waves may propagate with a circular polarization in each of the two planes of central circular section of the B^n-ellipsoid, where B is the left Cauchy-Green strain tensor corresponding to the primary pure homogeneous deformation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inhomogeneous Deformation of Brain Tissue During Tension Tests

Mechanical characterization of brain tissue has been investigated extensively by various research groups over the past fifty years. These properties are particularly important for modelling Traumatic Brain Injury (TBI). In this research, we present the design and calibration of a High Rate Tension Device (HRTD) capable of performing tests up to a maximum strain rate of 90/s. We use experimental and numerical methods to investigate the effects of inhomogeneous deformation of porcine brain tissue during tension at different specimen thicknesses (4.0-14.0 mm), by performing tension tests at a strain rate of 30/s. One-term Ogden material parameters (mu = 4395.0 Pa, alpha = -2.8) were derived by performing an inverse finite element analysis to model all experimental data. A similar procedure was adopted to determine Young's modulus (E= 11200 Pa) of the linear elastic regime. Based on this analysis, brain specimens of aspect ratio (diameter/thickness) S < 1.0 are required to minimise the effects of inhomogeneous deformation during tension tests.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inhomogeneous shear of orthotropic incompressible non-linearly elastic solids: singular solutions and biomechanical interpretation

We present a detailed study of rectilinear shear deformation in the framework of orthotropic nonlinear elasticity, under Dirichlet and mixed-boundary conditions. We take a slab made of a soft matrix, reinforced with two families of extensible fibers. We consider the case where the shear occurs along the bissectrix of the angle between the two privileged directions aligned with the fibers. We show that if the two families of parallel fibers are mechanically equivalent, then only smooth solutions are possible, whereas if the mechanical differences among the two families of fibers is pronounced, then strain singularities may develop. We determine the precise conditions for the existence of such singular solutions for the standard reinforcing orthotropic model. We then extend our findings to some orthotropic models of interest in biomechanical applications, and we discuss the possible relevance of the singular solutions to biomechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Initial stresses in elastic solids: Constitutive laws and acoustoelasticity

On the basis of the nonlinear theory of elasticity, the general constitutive equation for an isotropic hyperelastic solid in the presence of initial stress is derived. This derivation involves invariants that couple the deformation with the initial stress and in general, for a compressible material, it requires 10 invariants, reducing to 9 for an incompressible material. Expressions for the Cauchy and nominal stress tensors in a finitely deformed configuration are given along with the elasticity tensor and its specialization to the initially stressed undeformed configuration. The equations governing infinitesimal motions superimposed on a finite deformation are then used to study the combined effects of initial stress and finite deformation on the propagation of homogeneous plane waves in a homogeneously deformed and initially stressed solid of infinite extent. This general framework allows for various different specializations, which make contact with earlier works. In particular, connections with results derived within Biot's classical theory are highlighted. The general results are also specialized to the case of a small initial stress and a small pre-deformation, i.e. to the evaluation of the acoustoelastic effect. Here the formulas derived for the wave speeds cover the case of a second-order elastic solid without initial stress and subject to a uniaxial tension [Hughes and Kelly, Phys. Rev. 92 (1953) 1145] and are consistent with results for an undeformed solid subject to a residual stress [Man and Lu, J. Elasticity 17 (1987) 159]. These formulas provide a basis for acoustic evaluation of the second- and third-order elasticity constants and of the residual stresses. The results are further illustrated in respect of a prototype model of nonlinear elasticity with initial stress, allowing for both finite deformation and nonlinear dependence on the initial stress.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large acoustoelastic effect

Classical acoustoelasticity couples small-amplitude elastic wave propagation to an infinitesimal pre-deformation, in order to reveal and evaluate non-destructively third-order elasticity constants. Here, we see that acoustoelasticity can be also be used to determine fourth-order constants, simply by coupling a small-amplitude wave with a small-but-finite pre-deformation. We present results for compressible weakly nonlinear elasticity, we make a link with the historical results of Bridgman on the physics of high pressures, and we show how to determine "D", the so-called fourth-order elasticity constant of soft (incompressible, isotropic) solids by using infinitesimal waves.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Mechanical Characterization of Brain Tissue in Compression at Dynamic Strain Rates

Traumatic brain injury (TBI) occurs when local mechanical load exceeds certain tolerance levels for brain tissue. Extensive research has been done previously for brain matter experiencing compression at quasistatic loading; however, limited data is available to model TBI under dynamic impact conditions. In this research, an experimental setup was developed to perform unconfined compression tests and stress relaxation tests at strain rates < 90/s. The brain tissue showed a stiffer response with increasing strain rates, showing that hyperelastic models are not adequate. Specifically, the compressive nominal stress at 30% strain was 8.83 +/- 1.94, 12.8 +/- 3.10 and 16.0 +/- 1.41 kPa (mean +/- SD) at strain rates of 30, 60 and 90/s, respectively. Relaxation tests were also conducted at 10%-50% strain with the average rise time of 10 ms, which can be used to derive time dependent parameters. Numerical simulations were performed using one-term Ogden model with initial shear modulus mu_0 = 6.06 +/- 1.44, 9.44 +/- 2.427 and 12.64 +/- 1.227 kPa (mean +/- SD) at strain rates of 30, 60 and 90/s, respectively. A separate set of bonded and lubricated tests were also performed under the same test conditions to estimate the friction coefficient mu, by adopting combined experimental - computational approach. The values of mu were 0.1 +/- 0.03 and 0.15 +/- 0.07 (mean +/- SD) at 30 and 90/s strain rates, respectively, indicating that pure slip conditions cannot be achieved in unconfined compression tests even under fully lubricated test conditions. The material parameters obtained in this study will help to develop biofidelic human brain finite element models, which can subsequently be used to predict brain injuries under impact conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Non-principal surface waves in deformed incompressible materials

The Stroh formalism is applied to the analysis of infinitesimal surface wave propagation in a statically, finitely and homogeneously deformed isotropic half-space. The free surface is assumed to coincide with one of the principal planes of the primary strain, but a propagating surface wave is not restricted to a principal direction. A variant of Taziev's technique [Sov. Phys. Acoust. 35 (1989) 535] is used to obtain an explicit expression of the secular equation for the surface wave speed, which possesses no restrictions on the form of the strain energy function. Albeit powerful, this method does not produce a unique solution and additional checks are necessary. However, a class of materials is presented for which an exact secular equation for the surface wave speed can be formulated. This class includes the well-known Mooney-Rivlin model. The main results are illustrated with several numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear Correction to the Euler Buckling Formula for Compressed Cylinders with Guided-Guided End Conditions

Euler's celebrated buckling formula gives the critical load $N$ for the buckling of a slender cylindrical column with radius $B$ and length $L$ as \[ N / (π^3 B^2) = (E/4)(B/L)^2, \] where $E$ is Young's modulus. Its derivation relies on the assumptions that linear elasticity applies to this problem, and that the slenderness $(B/L)$ is an infinitesimal quantity. Here we ask the following question: What is the first nonlinear correction in the right hand-side of this equation when terms up to $(B/L)^4$ are kept? To answer this question, we specialize the exact solution of incremental non-linear elasticity for the homogeneous compression of a thick compressible cylinder with lubricated ends to the theory of third-order elasticity. In particular, we highlight the way second- and third-order constants ---including Poisson's ratio--- all appear in the coefficient of $(B/L)^4$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear transverse waves in deformed dispersive solids

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On interface waves in misoriented pre-stressed incompressible elastic solids

Some relationships, fundamental to the resolution of interface wave problems, are presented. These equations allow for the derivation of explicit secular equations for problems involving waves localized near the plane boundary of anisotropic elastic half-spaces, such as Rayleigh, Sholte, or Stoneley waves. They are obtained rapidly, without recourse to the Stroh formalism. As an application, the problems of Stoneley wave propagation and of interface stability for misaligned predeformed incompressible half-spaces are treated. The upper and lower half-spaces are made of the same material, subject to the same prestress, and are rigidly bonded along a common principal plane. The principal axes in this plane do not however coincide, and the wave propagation is studied in the direction of the bisectrix of the angle between a principal axis of the upper half-space and a principal axis of the lower half-space.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On Magnetoacoustic Waves in Finitely Deformed Elastic Solids

In this paper, in the context of the quasi-magnetostatic approximation, we examine incremental motions superimposed on a static finite deformation of a magnetoelastic material in the presence of an applied magnetic field. Explicit expressions are obtained for the associated magnetoacoustic (or magnetoelastic moduli) tensors in the case of an incompressible isotropic magnetoelastic material, and these are then used to study the propagation of incremental plane waves. The propagation condition is derived in terms of a generalized acoustic tensor and the results are illustrated by obtaining explicit formulas in two special cases: first, when the material is undeformed but subject to a uniform bias field and second for a prototype model of magnetoelastic interactions in the finite deformation regime. The results provide a basis for the experimental determination of the material parameters of a magneto-sensitive elastomer from measurements of the speed of incremental waves for different pre-strains, bias magnetic fields, and directions of propagation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On stress-dependent elastic moduli and wave speeds

On the basis of the general nonlinear theory of a hyperelastic material with initial stress, initially without consideration of the origin of the initial stress, we determine explicit expressions for the stress-dependent tensor of incremental elastic moduli. In considering three special cases of initial stress within the general framework, namely hydrostatic stress, uniaxial stress and planar shear stress, we then elucidate in general form the dependence of various elastic moduli on the initial stress. In each case the effect of initial stress on the wave speed of homogeneous plane waves is studied and it is shown how various special theories from the earlier literature fit within the general framework. We then consider the situation in which the initial stress is a pre-stress associated with a finite deformation and, in particular, we discuss the specialization to the second-order theory of elasticity and highlight connections between several classical approaches to the topic, again with special reference to the influence of higher-order terms on the speed of homogeneous plane waves. Some discrepancies arising in the earlier literature are noted.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Rectilinear Shear of Compressible and Incompressible Elastic Slabs

We review some pseudo-planar deformations for the equations of incompressible isotropic nonlinear elasticity first introduced in 1985 by Rajagopal and Wineman. We extend this class of deformations to compressible isotropic and transverse isotropic materials, and also consider the influence of gravity. We consider some new approximate solutions and we discuss the possible relevance of such solutions to the understanding of the complex structure of the fields equations of nonlinear elasticity, using weakly nonlinear theories.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the third- and fourth-order constants of incompressible isotropic elasticity

Consider the constitutive law for an isotropic elastic solid with the strain-energy function expanded up to the fourth order in the strain, and the stress up to the third order in the strain. The stress-strain relation can then be inverted to give the strain in terms of the stress with a view to considering the incompressible limit. For this purpose, use of the logarithmic strain tensor is of particular value. It enables the limiting values of all nine fourth-order elastic constants in the incompressible limit to be evaluated precisely and rigorously. In particular, it is explained why the three constants of fourth-order incompressible elasticity $μ$, $\bar{A}$, and $\bar D$ are of the same order of magnitude. Several examples of application of the results follow, including determination of the acoustoelastic coefficients in incompressible solids and the limiting values of the coefficients of nonlinearity for elastic wave propagation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Onset of Non-Linearity in the Elastic Bending of Blocks

The classical flexure problem of non-linear incompressible elasticity is revisited assuming that the bending angle suffered by the block is specified instead of the usual applied moment. The general moment-bending angle relationship is then obtained and is shown to be dependent on only one non-dimensional parameter: the product of the aspect ratio of the block and the bending angle. A Maclaurin series expansion in this parameter is then found. The first-order term is proportional to $μ$, the shear modulus of linear elasticity; the second-order term is identically zero, because the moment is an odd function of the angle; and the third-order term is proportional to $μ(4β-1)$, where $β$ is the non-linear shear coefficient, involving third-order and fourth-order elasticity constants. It follows that bending experiments provide an alternative way of estimating this coefficient, and the results of one such experiment are presented. In passing, the coefficients of Rivlin's expansion in exact non-linear elasticity are connected to those of Landau in weakly (fourth-order) non-linear elasticity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Proper formulation of viscous dissipation for nonlinear waves in solids

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rayleigh waves and surface stability for Bell materials in compression; comparison with rubber

The stability of a Bell-constrained half-space in compression is studied. To this end, the propagation of Rayleigh waves on the surface of the material when it is maintained in a static state of triaxial prestrain is considered. The prestrain is such that the free surface of the half-space is a principal plane of deformation. The exact secular equation is established for surface waves traveling in a principal direction of strain with attenuation along the principal direction normal to the free plane. As the half-space is put under increasing compressive loads, the speed of the wave eventually tends to zero and the bifurcation criterion, or stability equation, is reached. Then the analysis is specialized to specific forms of strain energy functions and prestrain, and comparisons are made with results previously obtained in the case of incompressible neo-Hookean or Mooney-Rivlin materials. It is found that these rubber-like incompressible materials may be compressed more than "Bell empirical model" materials, but not as much as "Bell simple hyperelastic" materials, before the critical stretches, solutions to the bifurcation criterion, are reached. In passing, some classes of incompressible materials which possess a relative-universal bifurcation criterion are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rayleigh waves in anisotropic crystals rotating about the normal to a symmetry plane

The propagation of surface acoustic waves in a rotating anisotropic crystal is studied. The crystal is monoclinic and cut along a plane containing the normal to the symmetry plane; this normal is also the axis of rotation. The secular equation is obtained explicitly using the "method of the polarization vector", and it shows that the wave is dispersive and decelerates with increasing rotation rate. The case of orthorhombic symmetry is also treated. The surface wave speed is computed for 12 monoclinic and 8 rhombic crystals, and for a large range of the rotation rate/wave frequency ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rayleigh waves in symmetry planes of crystals: explicit secular equations and some explicit wave speeds

Rayleigh waves are considered for crystals possessing at least one plane of symmetry. The secular equation is established explicitly for surface waves propagating in any direction of the plane of symmetry, using two different methods. This equation is a quartic for the squared wave speed in general, and a biquadratic for certain directions in certain crystals, where it may itself be solved explicitly. Examples of such materials and directions are found in the case of monoclinic crystals with the plane of symmetry at $x_3=0$. The cases of orthorhombic materials and of incompressible materials are also treated.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scalar evolution equations for shear waves in incompressible solids: A simple derivation of the Z, ZK, KZK, and KP equations

We study the propagation of two-dimensional finite-amplitude shear waves in a nonlinear pre-strained incompressible solid, and derive several asymptotic amplitude equations in a simple, consistent, and rigorous manner. The scalar Zabolotskaya (Z) equation is shown to be the asymptotic limit of the equations of motion for all elastic generalized neo-Hookean solids (with strain energy depending only on the first principal invariant of Cauchy-Green strain). However, we show that the Z equation cannot be a scalar equation for the propagation of two-dimensional shear waves in general elastic materials (with strain energy depending on the first and second principal invariants of strain). Then we introduce dispersive and dissipative terms to deduce the scalar Kadomtsev-Petviashvili (KP), Zabolotskaya-Khokhlov (ZK) and Khokhlov-Zabolotskaya-Kuznetsov (KZK) equations of incompressible solid mechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Shear instability in skin tissue

We propose two toy-models to describe, predict, and interpret the wrinkles appearing on the surface of skin when it is sheared. With the first model, we account for the lines of greatest tension present in human skin by subjecting a layer of soft tissue to a pre-stretch, and for the epidermis by endowing one of the layer's faces with a surface tension. For the second model, we consider an anisotropic model for the skin, to reflect the presence of stiff collagen fibres in a softer elastic matrix. In both cases, we find an explicit bifurcation criterion, linking geometrical and material parameters to a critical shear deformation accompanied by small static wrinkles, with decaying amplitudes normal to the free surface of skin.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Simple shear is not so simple

For homogeneous, isotropic, nonlinearly elastic materials, the form of the homogeneous deformation consistent with the application of a Cauchy shear stress is derived here for both compressible and incompressible materials. It is shown that this deformation is not simple shear, in contrast to the situation in linear elasticity. Instead, it consists of a triaxial stretch superposed on a classical simple shear deformation, for which the amount of shear cannot be greater than 1. In other words, the faces of a cubic block cannot be slanted by an angle greater than 45deg by the application of a pure shear stress alone. The results are illustrated for those materials for which the strain energy function does not depend on the principal second invariant of strain. For the case of a block deformed into a parallelepiped, the tractions on the inclined faces necessary to maintain the derived deformation are calculated.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Slight compressibility and sensitivity to changes in Poisson's ratio

Finite Element simulations of rubbers and biological soft tissue usually assume that the material being deformed is slightly compressible. It is shown here that in shearing deformations the corresponding normal stress distribution can exhibit extreme sensitivity to changes in Poisson's ratio. These changes can even lead to a reversal of the usual Poynting effect. Therefore the usual practice of arbitrarily choosing a value of Poisson's ratio when numerically modelling rubbers and soft tissue will, almost certainly, lead to a significant difference between the simulated and actual normal stresses in a sheared block because of the difference between the assumed and actual value of Poisson's ratio. The worrying conclusion is that simulations based on arbitrarily specifying Poisson's ratio close to 1/2 cannot accurately predict the normal stress distribution even for the simplest of shearing deformations. It is shown analytically that this sensitivity is due to the small volume changes which inevitably accompany all deformations of rubber-like materials. To minimise these effects, great care should be exercised to accurately determine Poisson's ratio before simulations begin.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Small-amplitude inhomogeneous plane waves in a deformed Mooney-Rivlin material

The propagation of small-amplitude inhomogeneous plane waves in an isotropic homogeneous incompressible Mooney--Rivlin material is considered when the material is maintained in a state of finite static homogeneous deformation. Disturbances of complex exponential type are sought and all propagating inhomogeneous solutions to the equations of motion are given, as well as the conditions for linear, elliptical, or circular polarization. % It is seen that a great variety of solutions arises. These include some original solutions, such as circularly-polarized plane waves which propagate with an arbitrary complex scalar slowness, or linearly-polarized waves for which the direction of propagation is not necessarily orthogonal to the direction of attenuation. % Throughout the paper, geometrical interpretations and explicit examples are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Some results on finite amplitude elastic waves propagating in rotating media

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stoneley waves and interface stability of Bell materials in compression; Comparison with rubber

Two semi-infinite bodies made of prestressed, homogeneous, Bell-constrained, hyperelastic materials are perfectly bonded along a plane interface. The half-spaces have been subjected to finite pure homogeneous predeformations, with distinct stretch ratios but common principal axes, and such that the interface is a common principal plane of strain. Constant loads are applied at infinity to maintain the deformations and the influence of these loads on the propagation of small-amplitude interface (Stoneley) waves is examined. In particular, the secular equation is found and necessary and sufficient conditions to be satisfied by the stretch ratios to ensure the existence of such waves are given. As the loads vary, the Stoneley wave speed varies accordingly: the upper bound is the `limiting speed' (given explicitly), beyond which the wave amplitude cannot decay away from the interface; the lower bound is zero, where the interface might become unstable. The treatment parallels the one followed for the incompressible case and the differences due to the Bell constraint are highlighted. Finally, the analysis is specialized to specific strain energy densities and to the case where the bimaterial is uniformly deformed (that is when the stretch ratios for the upper half-space are equal to those for the lower half-space.) Numerical results are given for `simple hyperelastic Bell' materials and for `Bell's empirical model' materials, and compared to the results for neo-Hookean incompressible materials.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Surface acoustic waves in rotating orthorhombic crystals

The propagation of surface (Rayleigh) waves over a rotating orthorhombic crystal is studied. The crystal possesses three crystallographic axes, normal to the symmetry planes: the half-space is cut along a plane normal to one of these axes, the wave travels in the direction of another, and the rotation occurs at a uniform rate about any of the three axes. The secular equation for the surface wave speed is found explicitly; in contrast to the non-rotating case, it is dispersive (frequency-dependent). Both Coriolis and centrifugal accelerations appear in the equations of motion: none can be neglected in favor of the other, even at small rotation rates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Surface waves in a deformed isotropic hyperelastic material subject to an isotropic internal constraint

An isotropic elastic half space is prestrained so that two of the principal axes of strain lie in the bounding plane, which itself remains free of traction. The material is subject to an isotropic constraint of arbitrary nature. A surface wave is propagated sinusoidally along the bounding surface in the direction of a principal axis of strain and decays away from the surface. The exact secular equation is derived by a direct method for such a principal surface wave; it is cubic in a quantity whose square is linearly related to the squared wave speed. For the prestrained material, replacing the squared wave speed by zero gives an explicit bifurcation, or stability, criterion. Conditions on the existence and uniqueness of surface waves are given. The bifurcation criterion is derived for specific strain energies in the case of four isotropic constraints: those of incompressibility, Bell, constant area, and Ericksen. In each case investigated, the bifurcation criterion is found to be of a universal nature in that it depends only on the principal stretches, not on the material constants. Some results related to the surface stability of arterial wall mechanics are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Surface waves in a stretched and sheared incompressible elastic material

In this paper we analyze the effect of a combined pure homogeneous strain and simple shear in a principal plane of the latter on the propagation of surface waves for an incompressible isotropic elastic half-space whose boundary is normal to the glide planes of the shear. This generalizes previous work in which, separately, pure homogeneous strain and simple shear were considered. For a special class of materials the secular equation is obtained in explicit form and then specialized to recover results obtained previously for the two cases mentioned above. A method for obtaining the secular equation for a general form of strain-energy function is then outlined. In general this is very lengthy and the result is not listed, but, for the case in which there is no normal stress on the half-space boundary, the result is given, for illustration, in respect of the so-called generalized Varga material. Numerical results are given to show how the surface wave speed depends on both the underlying pure homogeneous strain and the superimposed simple shear. Further numerical results are provided for the Gent model of limiting chain extensibility.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Surface waves in deformed Bell materials

Small amplitude inhomogeneous plane waves are studied as they propagate on the free surface of a predeformed semi-infinite body made of Bell constrained material. The predeformation corresponds to a finite static pure homogeneous strain. The surface wave propagates in a principal direction of strain and is attenuated in another principal direction, orthogonal to the free surface. For these waves, the secular equation giving the speed of propagation is established by the method of first integrals. This equation is not the same as the secular equation for incompressible half-spaces, even though the Bell constraint and the incompressibility constraint coincide in the isotropic infinitesimal limit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Temperature Effects on Brain Tissue in Compression

Extensive research has been carried out for at least 50 years to understand the mechanical properties of brain tissue in order to understand the mechanisms of traumatic brain injury (TBI). The observed large variability in experimental results may be due to the inhomogeneous nature of brain tissue and to the broad range of test conditions. However, test temperature is also considered as one of the factors influencing the properties of brain tissue. In this research, the mechanical properties of porcine brain have been investigated at 22C (room temperature) and at 37C (body temperature) while maintaining a constant preservation temperature of approximately 4-5C. Unconfined compression tests were performed at dynamic strain rates of 30 and 50/s using a custom made test apparatus. There was no significant difference (p = 0.8559 - 0.9290) between the average engineering stresses of the brain tissue at the two different temperature conditions. The results of this study should help to understand the behavior of brain tissue at different temperature conditions, particularly in unconfined compression tests.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The explicit secular equation for surface acoustic waves in monoclinic elastic crystals

The secular equation for surface acoustic waves propagating on a monoclinic elastic half-space is derived in a direct manner, using the method of first integrals. Although the motion is at first assumed to correspond to generalized plane strain, the analysis shows that only two components of the mechanical displacement and of the tractions on planes parallel to the free surface are nonzero. Using the Stroh formalism, a system of two second order differential equations is found for the remaining tractions. The secular equation is then obtained as a quartic for the squared wave speed. This explicit equation is consistent with that found in the orthorhombic case. The speed of subsonic surface waves is then computed for twelve specific monoclinic crystals.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The incompressible limit in linear anisotropic elasticity, with application to surface waves and elastostatics

Incompressibility is established for three-dimensional and two-dimensional deformations of an anisotropic linearly elastic material, as conditions to be satisfied by the elastic compliances. These conditions make it straightforward to derive results for incompressible materials from those established for the compressible materials. As an illustration, the explicit secular equation is obtained for surface waves in incompressible monoclinic materials with the symmetry plane at x_3=0. This equation also covers the case of incompressible orthotropic materials. The displacements and stresses for surface waves are often expressed in terms of the elastic stiffnesses, which can be unbounded in the incompressible limit. An alternative formalism in terms of the elastic compliances presented recently by Ting is employed so that surface wave solutions in the incompressible limit can be obtained. A different formalism, also by Ting, is employed to study the solutions to two-dimensional elastostatic problems. In the special case of incompressible monoclinic material with the symmetry plane at x_3=0, one of the three Barnett-Lothe tensors S vanishes while the other two tensors H and L are the inverse of each other. Moreover, H and L are diagonal with the first two diagonal elements being identical. An interesting physical phenomenon deduced from this property is that there is no interpenetration of the interface crack surface in an incompressible bimaterial. When only the inplane deformation is considered, it is shown that the image force due to a line dislocation in a half-space or in a bimaterial depends only on the magnitude, not on the direction, of the Burgers vector.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The speed of interfacial waves polarized in a symmetry plane

The surface-impedance matrix method is used to study interfacial waves polarized in a plane of symmetry of anisotropic elastic materials. Although the corresponding Stroh polynomial is a quartic, it turns out to be analytically solvable in quite a simple manner. A specific application of the result concerns the calculation of the speed of a Stoneley wave, polarized in the common symmetry plane of two rigidly bonded anisotropic solids. The corresponding algorithm is robust, easy to implement, and gives directly the speed (when the wave exists) for any orientation of the interface plane, normal to the common symmetry plane. Through the examples of the couples (Aluminum)-(Tungsten) and (Carbon/epoxy)-(Douglas pine), some general features of a Stoneley wave speed are verified: the wave does not always exist; it is faster than the slowest Rayleigh wave associated with the separated half-spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Third- and fourth-order elasticity of biological soft tissues

In the theory of weakly non-linear elasticity, Hamilton et al. [J. Acoust. Soc. Am. \textbf{116} (2004) 41] identified $W = μI_2 + (A/3)I_3 + D I_2^2$ as the fourth-order expansion of the strain-energy density for incompressible isotropic solids. Subsequently, much effort focused on theoretical and experimental developments linked to this expression in order to inform the modeling of gels and soft biological tissues. However, while many soft tissues can be treated as incompressible, they are not in general isotropic, and their anisotropy is associated with the presence of oriented collagen fiber bundles. Here the expansion of $W$ is carried up to fourth-order in the case where there exists one family of parallel fibers in the tissue. The results are then applied to acoustoelasticity, with a view to determining the second- and third-order nonlinear constants by employing small-amplitude transverse waves propagating in a deformed soft tissue.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Transverse Waves in Nonlinearly Elastic Solids and the Milne-Pinney (or Ermakov) Equation

We establish a connection between the general equations of nonlinear elastodynamics and the nonlinear ordinary differential equation of Pinney [Proc. Amer. Math. Soc. 1 (1950) 681]. As a starting point, we use the exact travelling wave solutions of nonlinear elasticity discovered by Carroll [Acta Mechanica 3 (1967) 167]. The connection provides a method for finding new exact and approximate dynamic solutions for neo-Hookean and Mooney-Rivlin solids, and for the general third- and fourth-order elasticity models of incompressible solids.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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 was observed. The structural parameters of the Gasser-Ogden-Holzapfel (GOH) model were all successfully evaluated. 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.

0 citations0 reviewsNo code linkedOpen access