Search

International audience; We propose a chemical-mechanical model of myosin heads in sarcomeres, within the classical description of rigid sliding filaments. In our case, myosin heads have two mechanical degrees-of-freedom (dofs) - one of which associated with the so-called power stroke - and two possible chemical states, i.e. bound to an actin site or not. Our major motivations are twofold: (1) to derive a multiscale coupled chemical-mechanical model, and (2) to thus account - at the macroscopic scale - for mechanical phenomena that are out of reach for classical muscle models. This model is first written in the form of Langevin stochastic equations, and we are then able to obtain the corresponding Fokker-Planck partial differential equations governing the probability density functions associated with the mechanical dofs and chemical states. This second form is important, as it allows to monitor muscle energetics, and also to compare our model with classical ones, such as the Huxley'57 model to which our equations are shown to reduce under two different types of simplifying assumptions. This provides insight, and gives a Langevin form for Huxley'57. We then show how we can calibrate our model based on experimental data - taken here for skeletal muscles - and numerical simulations demonstrate the adequacy of the model to represent complex physiological phenomena, in particular the fast isometric transients in which the power stroke is known to have a crucial role, thus circumventing a limitation of many classical models.

Despite recent advances, the myocardial microstructure remains imperfectly understood. In particular, bundles of cardiomyocytes have been observed but their three-dimensional organisation remains debated and the associated mechanical consequences unknown. One of the major challenges remains to perform multiscale observations of the mechanical response of the heart wall. For this purpose, in this study, a full-field Mueller polarimetric imager (MPI) was combined, for the first time, with an in-situ traction device. The full-field MPI enables to obtain a macroscopic image of the explored tissue, while providing detailed information about its structure on a microscopic scale. Specifically it exploits the polarization of the light to determine various biophysical quantities related to the tissue scattering or anisotropy properties. Combined with a mechanical traction device, the full-field MPI allows to measure the evolution of such biophysical quantities during tissue stretch. We observe separation lines on the tissue, which are associated with a fast variation of the fiber orientation, and have the size of cardiomyocyte bundles. Thus, we hypothesize that these lines are the perimysium, the collagen layer surrounding these bundles. During the mechanical traction, we observe two mechanisms simultaneously. On one hand, the azimuth shows an affine behavior, meaning the orientation changes according to the tissue deformation, and showing coherence in the tissue. On the other hand, the separation lines appear to be resistant in shear and compression but weak against traction, with a forming of gaps in the tissue.

Of interest in this work are nematic continua that exhibit electromechanical coupling. The first part of this paper presents a novel variational formulation with a potential energy depending on four independent variables (the displacement, director, specific polarization and electric displacement perturbation). Variations of the potential energy with respect to each one of these variables lead to the governing mechanical equilibrium and constitutive relations plus Maxwell’s equations. The proposed variational formulation is next applied to the study of bifurcation of an infinite layer of a nematic liquid crystal confined between two parallel plates and subjected to a uniform electric field perpendicular to these plates under full anchoring boundary conditions. As the electric field exceeds a critical value, the nematic directors which are initially parallel to the plates, rotate and tend to align with the electric field orientation. This phenomenon, termed in the literature as Freedericksz transition, is treated here as a bifurcation problem using a fully 2D formulation. It is shown that the solution corresponding to the lowest applied electric field, also termed the critical load, is uniform in the direction parallel to the plates and that the corresponding bifurcated path is stable near this critical load. This result holds for arbitrary positive constants of the Frank-Oseen energy (and values of electric susceptibility constants that allow bifurcation) and justifies the 1D treatment of the Freedericksz transition in 2D settings that is widely adopted in the liquid crystal literature. An asymptotic analysis of the supercritical, stable bifurcated equilibrium path about the critical load is also presented and compared with the exact bifurcated solution.

International audience; An unresolved issue in patient-specific models of cardiac mechanics is the choice of an appropriate constitutive law, able to accurately capture the passive behavior of the myocardium, while still having uniquely identifiable parameters tunable from available clinical data. In this paper, we aim to facilitate this choice by examining the practical identifiability and model fidelity of constitutive laws often used in cardiac mechanics. Our analysis focuses on the use of novel 3D tagged MRI, providing detailed displacement information in three dimensions. The practical identifiability of each law is examined by generating synthetic 3D tags from in silico simulations, allowing mapping of the objective function landscape over parameter space and comparison of minimizing parameter values with original ground truth values. Model fidelity was tested by comparing these laws with the more complex transversely isotropic Guccione law, by characterizing their passive end-diastolic pressure–volume relation behavior, as well as by considering the in vivo case of a healthy volunteer. These results show that a reduced form of the Holzapfel–Ogden law provides the best balance between identifiability and model fidelity across the tests considered.

International audience; This paper investigates the effects of plasticity on the effective fracture toughness. A layered material is considered as a modelling system. An elastic-plastic phase-field model and a surfing boundary condition are used to study how the crack propagates throughout the material and the evolution of the effective toughness as a function of the layer angle. We first study three idealized situations, where only one property among fracture toughness, Young's modulus and yield strength is heterogeneous whereas the others are uniform. We observe that in the case of toughness and strength heterogeneity, the material exhibits anomalous isotropy: the effective toughness is equal to the largest of the point-wise values for any layer angle except when the layers are parallel to the macroscopic direction of propagation. As the layer angle decreases, the crack propagates along the brittle-to-tough interfaces, whereas it goes straight when the layers have different yield strength but uniform toughness. We find that smooth deflections in the crack path do not induce any overall toughening and that the effective toughness is not proportional to either the cumulated fracture energy or the cumulated plastic work. In the case of elastic heterogeneity, the material is anisotropic in the sense of the effective toughness, as the latter varies as a function of the layer angle. Four toughening mechanisms are active: stress fluctuations, crack renucleation, plastic dissipation and plastic blunting. Finally, we consider a layered medium comprised of compliant-tough-weak and stiff-brittle-strong phases, as it is the case for many structural composites. We observe a transition from an interface-dominated to a plasticity-dominated failure regime, as the phase constituents become more ductile. The material is anisotropic in the sense of the effective toughness.

International audience; The viscoplastic behavior of coarse and bimodal fine/ultrafine grained (F/UFG) Al5083 alloywas investigated between 20 °C and 200 °C through tensile tests at various strain rates, andstress relaxation tests to deduce the strain rate sensitivity (SRS). The plastic strain fields weremeasured by correlation of SEM images. In the F/UFG material at high temperature, very highstrains were measured in shear bands which sometimes crossed the whole gage width andexhibited intensive grain boundary sliding (GBS). Both the SRS and ductility rose with thetemperature, and as the strain rate decreased, mainly due to a rising contribution of GBS, whichaccommodated a much larger fraction of the global strain in the F/UFG material. Theboundary between the temperature–strain rate domains where grain refinement led either tostrengthening or to softening was determined. Finite element simulations of tension andrelaxation tests with viscoplastic grains and sliding grain boundaries captured the macro-scalebehavior of the F/UFG material. It also provided some insight into the mechanisms ofcorrelated and cooperative GBS and grain rotation along percolation paths (both inter andintragranular), probably, responsible for macro shear banding.

International audience; The resistance to mode I failure of rubbers is studied by submitting single edge notch samples to uniaxial tension. Reproducing the seminal work of Rivlin and Thomas (J Polym Sci 10:291–318, 1953), single edge notch tension specimens, presenting notches of various lengths, are stretched until break. A styrene butadiene rubber, unfilled and filled with carbon-black, and an unfilled rubber from the latter mentioned work, were considered. When the notch is smaller than one fifth of the sample width, mode I crack opening is observed, leading to catastrophic failure that creates smooth mirror-like crack surfaces. Nonetheless, the experimental force-elongation responses show that the mode I critical energy release rate cannot be calculated by a classical Griffith elastic failure analysis. When notches are longer, the SENT samples are not submitted to pure uniaxial tension only. Structural bending leads to uncontrolled mixed mode crack propagation. The surfaces created when the long notches propagate are rough and bifurcations are witnessed for the filled rubbers.

International audience; This paper focuses on the development of an appropriate Digital Image Correlation (DIC) methodology based on Image Registration and dedicated for characterizing the plastic deformation in single crystals. A pure nickel single crystal specimen is plastically deformed in tension and investigated by DIC technique. Based on the measured kinematic fields, the proposed method enables to identify the slip activity on the crystal surface and to locate precisely the slip band interfaces at microscale which behave as kinematic discontinuities. The computed displacement data are projected on a well-defined physical basis containing slip details, then the strain fields can be derived directly from a set of analytical functions. The possible errors in displacement induced by this projection approach are evaluated. Finally, some results of the evaluated strain fields are presented. It demonstrates that the developed DIC methodology allows quantitative characterization of a heterogeneous deformation process and promotes further relationships to be established between slip activity and strain field evolution in single crystals.

International audience; This paper presents a rate-independent ana- lytical model for porous Tresca ( J3 -dependent) mate- rials containing general ellipsoidal voids. The model is based on the nonlinear variational homogenization method which uses a linear comparison material to esti- mate the response of the nonlinear porous solid and is denoted as “MVAR”. Specifically, the model is derived by an original approach starting from a novel porous single crystal model (Mbiakop et al. in Int J Solids Struct 64–65:100–119, 2015b, J Mech Phys Solids 84:436–467, 2015c) by considering the limiting case of infinite slip systems which leads readily to the cor- responding Tresca criterion. The MVAR yield surface is then validated using FEM on different unit-cells and various parameters including several porosity levels, several stress triaxiality ratios, different Lode angle and general void shapes and orientations. The MVAR model is found to be in good agreement with the finite element results for all cases considered in this study. Both the MVAR and the FEM computations reveal a strong sensitivity upon the microstructure anisotropy (void shape and orientation), and a dependence of the effective behavior on the third invariant of the applied stress. To the best knowledge of the authors, this is the first model in the literature that is able to deal with porous Tresca material and general void shapes and orientations. Moreover, the MVAR is used in a pre- dictive manner to investigate the complex response of porous Tresca cases with strong coupling between the J3-dependent matrix behavior and the (morphologi- cal) anisotropy induced by the shape and orientation of the voids. The simplicity of the present analytical study opens the possibility to adapt the present model to experimental results for various materials.

International audience; The great success of stents in treating cardiovas-cular disease is actually undermined by their long-term fatigue failure. The high variability of stent failure incidence suggests that it is due to several correlated aspects, such as loading conditions, material properties, component design, surgical procedure, and patient functional anatomy. Numerical and experimental non-clinical assessments are included in the recommendations and requirements of several regulatory bodies and they are thus exploited in the analysis of stent fatigue performance. Optimization-based simulation methodologies have been developed as well, to improve the fatigue endurance of novel designs. This paper presents a review on the fatigue issue in metallic stents, starting from a description of clinical evidence about stent fracture up to the analysis of computational approaches available from the literature. The reported discussion on both the experimental and numerical framework aims at providing a general insight into stent lifetime prediction as well as at understanding the factors which affect stent fatigue performance for the design of novel components.

CSF leakages constitute a major complication of intradural procedures, especially for posterior fossa and skull base surgery. Dural suture watertightness is a decisive issue, and neurosurgeons routinely use different products to reinforce their dural closure. We have designed an experimental system capable of testing CSF leak pressure levels in order to compare two types of sutures in vitro and particularly four different sealants.Twenty-five fresh human cadaveric dural samples were removed and prepared for testing in a pressure chamber system connected to a hydraulic pressure motor. CSF leak levels were objectively registered. First, simple interrupted stitches were compared to running simple closure on 50-mm linear suture. Secondly, four sealants (two sealants/glues, Bioglue®, Duraseal®; two haemostatics, Tachosil®, Tissucol®) were tested. Statistical analysis was performed with paired Student's t-test.No significant difference between interrupted closure and running suture was observed (p = 0.079). All sealants increased the watertightness of the suture significatively. However, comparison of the means of the differences for each product revealed large variations. In the conditions of our experiment, one sealant (Duraseal®) and one haemostatic (Tachosil®) seemed to show better results. We observed two different types of leakage: at the dura-sealant interface and through the sealant itself.We have developed an experimental device capable of testing dural closure watertightness. Interrupted stitch suturing seemed no different from running simple closure. On the contrary, the sealants tested show different watertightness capacities.

International audience; In the case of Lipschitz extremals of vectorial variational problems, an important class of strong variations originates from smooth deformations of the corresponding non-smooth graphs. These seemingly singular variations, which can be viewed as combinations of weak inner and outer variations, produce directions of differentiability of the functional and lead to singularity-centered necessary conditions on strong local minima: an equality, arising from stationarity, and an inequality, implying configurational stability of the singularity set. To illustrate the underlying coupling between inner and outer variations, we study in detail the case of smooth surfaces of gradient discontinuity representing, for instance, martensitic phase boundaries in non-linear elasticity.

International audience; We derive a new general jump condition on a broken Weierstrass–Erdmann extremal of a vectorial variational problem. Such extremals, containing surfaces of gradient discontinuity, are ubiquitous in shape optimization and in the theory of elastic phase transformations. The new condition, which does not have a one dimensional analog, reflects the stationarity of the singular surface with respect to two-scale variations that are nontrivial generalizations of Weierstrass needles. The over-determinacy of the ensuing free boundary problem suggests that typical stable solutions must involve microstructures or chattering controls.

International audience; We study Hamiltonian dynamics of a Fermi–Pasta–Ulam (FPU) chain with bi-stable elements. We show, that a quasi-static driving of a ‘cold’ chain beyond the spinodal threshold leads to complex dynamical behavior involving equipartition which suggests thermalization. The subsequent quasi-static cycling between the two energy wells produces reversible temperature oscillations which we link to the release (or absorbtion) of the latent heat. By adopting canonical distribution we obtain a thermodynamical description of the chain which agrees well with numerically computed time-averaged behavior of the corresponding dynamical system.

International audience; Within continuum dislocation theory, one-dimensional energy functional of a bent beam, made of a single crystal, is derived. By relaxing the continuously differentiable minimizer of this energy functional, we construct a sequence of piecewise smooth deflections and piecewise constant plastic distortions reducing the energy and exhibiting polygonization. The number of polygons can be estimated by comparing the surface energy of small angle tilt boundaries with the contribution of the gradient terms from the weak minimizer in the bulk energy.

International audience; The paper addresses the identification of a planar crack for Zener type linear viscoelastic solids. Under the condition of low frequency, the Zener model of viscoelasticity establishes the equivalence between viscoelasticity and elasticity and the equations are reduced to a Helmholtz type problem for time harmonic loadings. The solutions to the crack identification problems are then obtained from the corresponding solutions in elasticity, using only one frequency.

International audience; The influence of hydrostatic and uniaxial stress states on the porosity and permeability of sandstones has been investigated. The experimental procedure uses a special triaxial cell which allows permeability measurements in the axial and radial directions. The core sleeve is equipped with two pressure samplers placed distant from the ends. They provide mid-length axial permeability measure as opposed to the overall permeability measure, which is based on the flow imposed through the pistons of the triaxial cell. The core sleeve is also equipped to perform flows in two directions transverse to the axis of the sample. Two independent measures of axial and complementary radial permeability are thus obtained. Both Fontainebleau sandstone specimens with a porosity of about 5.8% to 8% and low permeability ranging from 2.5 mD to 30 mD and Bentheimer sandstone with a porosity of 24% and a high permeability of 3D have been tested. The initial axial permeability values obtained by each method are in good agreement for the Fontainebleau sandstone. The Bentheimer sandstone samples present an axial mid-length permeability 1.6 times higher than the overall permeability. A similar discrepancy is also observed in the radial direction, also it relates essentially to the shape of flow lines induced by the radial flow. All the tested samples have shown a higher stress dependency of overall and radial permeability than mid-length permeability. The effect of compaction damage at the pistons/sample and radial ports/sample interfaces is discussed. The relevance of directional permeability measurements during continuous uniaxial compression loadings has been shown on the Bentheimer sandstone until the failure of the sample. We can efficiently measure the influence of brittle failure associated to dilatant regime on the permeability: It tends to increase in the failure propagation direction and to decrease strongly in the transverse direction.

International audience; Optical full-field measurement methods such as Digital Image Correlation (DIC) are increasingly used in the field of experimental mechanics, but they still suffer from a lack of information about their metrological performances. To assess the performance of DIC techniques and give some practical rules for users, a collaborative work has been carried out by the Workgroup "Metrology" of the French CNRS research network 2519 "MCIMS (Mesures de Champs et Identification en Mécanique des Solides / Full-field measurement and identification in solid mechanics, http://www.ifma.fr/lami/gdr2519)". A methodology is proposed to assess the metrological performances of the image processing algorithms that constitute their main component, the knowledge of which being required for a global assessment of the whole measurement system. The study is based on displacement error assessment from synthetic speckle images. Series of synthetic reference and deformed images with random patterns have been generated, assuming a sinusoidal displacement field with various frequencies and amplitudes. Displacements are evaluated by several DIC packages based on various formulations and used in the French community. Evaluated displacements are compared with the exact imposed values and errors are statistically analyzed. Results show general trends rather independent of the implementations but strongly correlated with the assumptions of the underlying algorithms. Various error regimes are identified, for which the dependence of the uncertainty with the parameters of the algorithms, such as subset size, gray level interpolation or shape functions, is discussed.

The purpose of this article is to revise some concepts on defects nucleation based on bifurcation of equilibrium solutions. Equilibrium solutions are obtained on a homogeneous body and on a body with an infinitesimal defect such as cavity under the same prescribed dead load. First void formation and growth in non linear mechanics are examined. A branch of radial transformation bifurcates from the undeformed configuration in presence of a small cavity. Two cases of behaviour are examined. One case is the growth of the cavity by only the deformation of the shell. In another modelling the cavity evolves like a damaged zone, the transition between the sound part and the damaged one is governed by a local criterium. Each configuration leads to the definition of a nucleation criterion based on a presence of a bifurcation state, common state of the homogeneous body and a body with an infinitesimal defect.

International audience; An aluminum alloy(1) was tested at quasi-static to dynamic strain-rates (from 10(-1) to 5 10(3) s(-1)), using a single measuring device, a modified Split Hopkinson Bar. A wave separation technique [Bussac et al., J Mech Phys Solids 50:321-350, 2002] based on the maximum likelihood method was applied to process the strain and velocity measurements recorded at various points on each bar. With this method, it is possible to compute the stress, strain, displacement and velocity at any point on the bar. Since the measurement time is unlimited, the maximum strain measured in a given specimen no longer decreases with the strain-rate, as occurs with the classical Split Hopkinson Bar method.

International audience; We propose a series of quasicontinuum approximations for the simplest lattice model of a fully dynamic martensitic phase transition in one dimension. The approximations are dispersive and include various non-classical corrections to both kinetic and potential energies. We show that the well-posed quasicontinuum theory can be constructed in such a way that the associated closed-form kinetic relation is in an excellent agreement with the predictions of the discrete theory.

International audience; The determination of the chloride diffusion coefficient of a concrete is needed to help the prediction of the service life of concrete structure. In this paper, we propose first a critical review of models for chloride diffusion coefficients already used in literature at different scales and then we develop an analytical model, which takes into account the characteristics of the different phases of concrete. These materials are treated as a three-phase composite, consisting of a cement continuous phase, of an aggregates dispersed phase and of an interface transition zone. Chloride diffusion coefficient using an n-layered inclusion-based micromechanical modeling is predicted. The details of calculations are summarized hereafter and experimental measurements obtained on mortars are compared with predicted results.

A boundary-domain integral equation for a coated half-space (elastically isotropic homogeneous substratum, possibly anisotropic coating layer) is developed. The half-space fundamental solution is used, so that the discretization is limited to the potential contact zone (boundary elements), the potentially plastic part of the substratum and the coating layer (domain integration cells). Steady-state elastoplastic analysis is implemented within this framework, for plane-strain conditions, for solving rolling and/or sliding contact problems, where at the moment the contact load comes from either a purely elastic contact analysis or is of Hertz type. The constitutive integration is of implicit type. In order to improve accuracy and computational efficiency, infinite elements are used. Comparison of numerical results with other sources, when available, is satisfactory. The present formulation is also used to compute the contact pressure for an isotropic (or anisotropic) coating on an isotropic homogeneous half-space indented by an elastic punch.