BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Michel Destrade

Michel Destrade contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.softphysics.class-phBiological Physicsnlin.PSphysics.app-phcond-mat.mtrl-sciTissues and Organs

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
39works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

39 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&#39;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