Estimation of Isotropic Hyperelasticity Constitutive Models to Approximate the Atomistic Simulation Data for Aluminium and Tungsten Monocrystals

In this paper,the choice and parametrisation of finite deformation polyconvex isotropic hyperelastic models to describe the behaviour of a class of defect-free monocrystalline metal materials at the molecular level is examined.The article discusses some physical,mathematical and numerical demands which in our opinion should be fulfilled by elasticity models to be useful.A set of molecular numerical tests for aluminium and tungsten providing data for the fitting of a hyperelastic model was performed,and an algorithm for parametrisation is discussed.The proposed models with optimised parameters are superior to those used in non-linear mechanics of crystals.

NEW CONSTITUTIVE RELATIONSHIP OF INCOMPRESSIBLE HYPERELASTICITY

Starting from an incremental stress strain relationship of the form of Hooke's law, the paper derives expressions of a new constitutive relation for incompressible hyperelasticity. Assuming the stiffness coefficients being a power series of strain invari

A MESH-FREE ALGORITHM FOR DYNAMIC IMPACT ANALYSIS OF HYPERELASTICITY

A mesh-free method based on local Petrov-Galerkin formulation is presented to solve dynamic impact problems of hyperelastic material.In the present method,a simple Heaviside test function is chosen for simplifying domain integrals.Trial function is constructed by using a radial basis function （RBF） coupled with a polynomial basis function,in which the shape function possesses the kronecker delta function property.So,additional treatment is not required for imposing essential boundary conditions.Governing equations of impact problems are established and solved node by node by using an explicit time integration algorithm in a local domain,which is very similar to that of the collocation method except that numerical integration can be implemented over local domain in the present method.Numerical results for several examples show that the present method performs well in dealing with the dynamic impact problem of hyperelastic material.

Isogeometric Analysis of Hyperelastic Material Characteristics for Calcified Aortic Valve

This study explores the implementation of computed tomography(CT)reconstruction and simulation techniques for patient-specific valves,aiming to dissect the mechanical attributes of calcified valves within transcatheter heart valve replacement(TAVR)procedures.In order to facilitate this exploration,it derives pertinent formulas for 3D multi-material isogeometric hyperelastic analysis based on Hounsfield unit(HU)values,thereby unlocking foundational capabilities for isogeometric analysis in calcified aortic valves.A series of uniaxial and biaxial tensile tests is executed to obtain an accurate constitutive model for calcified active valves.To mitigate discretization errors,methodologies for reconstructing volumetric parametric models,integrating both geometric and material attributes,are introduced.Applying these analytical formulas,constitutive models,and precise analytical models to isogeometric analyses of calcified valves,the research ascertains their close alignment with experimental results through the close fit in displacement-stress curves,compellingly validating the accuracy and reliability of the method.This study presents a step-by-step approach to analyzing themechanical characteristics of patient-specific valves obtained fromCT images,holding significant clinical implications and assisting in the selection of treatment strategies and surgical intervention approaches in TAVR procedures.

Determination of Material Parameters of EVA Foam under Uniaxial Compressive Testing Using Hyperelastic Models

The objective of this research was to determine the mechanical parameter from EVA foam and also investigate its behavior by using Blatz-Ko,Neo-Hookean,Mooney model and experimental test.The physical characteristic of EVA foam was also evaluated by scanning electron microscopy(SEM).The results show that Blatz-Ko and Neo-Hookean model can fit the curve at 5%and 8%strain,respectively.The Mooney model can fit the curve at 50%strain.The modulus of rigidity evaluated from Mooney model is 0.0814±0.0027 MPa.The structure of EVA foam from SEM image shows that EVA structure is a closed cell with homogeneous porous structure.From the result,it is found that Mooney model can adjust the data better than other models.This model can be applied for mechanical response prediction of EVA foam and also for reference value in engineering application.

A review on material models for isotropic hyperelasticity

Dozens of hyperelastic models have been formulated and have been extremely handy in understanding the complex mechanical behavior of materials that exhibit hyperelastic behavior(characterized by large nonlinear elastic deformations that are completely recoverable)such as elastomers,polymers,and even biological tis-sues.These models are indispensable in the design of complex engineering com-ponents such as engine mounts and structural bearings in the automotive and aerospace industries and vibration isolators and shock absorbers in mechanical systems.Particularly,the problem of vibration control in mechanical system dy-namics is extremely important and,therefore,knowledge of accurate hyperelastic models facilitates optimum designs and the development of three‐dimensional finite element system dynamics for studying the large and nonlinear deformation beha-vior.This review work intends to enhance the knowledge of 15 of the most com-monly used hyperelastic models and consequently help design engineers and scientists make informed decisions on the right ones to use.For each of the models,expressions for the strain‐energy function and the Cauchy stress for both arbitrary loading assuming compressibility and each of the three loading modes(uniaxial tension,equibiaxial tension,and pure shear)assuming incompressibility are pro-vided.Furthermore,the stress–strain or stress–stretch plots of the model's pre-dictions in each of the loading modes are compared with that of the classical experimental data of Treloar and the coefficient of determination is utilized as a measure of the model's predictive ability.Lastly,a ranking scheme is proposed based on the model's ability to predict each of the loading modes with minimum deviations and the overall coefficient of determination.

Active Truss Metamaterials: Modelling and Tuning of Band Gaps

Periodic composite structures, like acoustic metamaterials (AMMs) and phononic crystals (PCs), are able to prevent the propagation of sound and elastic waves for some specific frequency ranges, leading to the emergence of so-called band gaps. So far, the optimization of the metamaterial properties and therefore of the band gaps has been typically performed on passive PCs and AMMs. Hence, the band gap properties cannot be tuned anymore after the production process of the metamaterials;this problem can be overcome thanks to the use of active materials. In this work, material and geometric nonlinearities are exploited to actively tune the frequency ranges of the band gaps of an architected AMM characterized by a three-dimensional periodicity. Specifically, a hyperelastic piezoelectric composite, that can be obtained by embedding piezo nanoparticles in a soft polymeric matrix, is considered to assess the effects of the nonlinearities on the behavior of sculptured microstructures, taking advantage of instability-induced pattern transformation and piezoelectricity to actively tune the band gaps. .

Active Truss Metamaterials: Modelling and Tuning of Band Gaps

Periodic composite structures, like acoustic metamaterials (AMMs) and phononic crystals (PCs), are able to prevent the propagation of sound and elastic waves for some specific frequency ranges, leading to the emergence of so-called band gaps. So far, the optimization of the metamaterial properties and therefore of the band gaps has been typically performed on passive PCs and AMMs. Hence, the band gap properties cannot be tuned anymore after the production process of the metamaterials;this problem can be overcome thanks to the use of active materials. In this work, material and geometric nonlinearities are exploited to actively tune the frequency ranges of the band gaps of an architected AMM characterized by a three-dimensional periodicity. Specifically, a hyperelastic piezoelectric composite, that can be obtained by embedding piezo nanoparticles in a soft polymeric matrix, is considered to assess the effects of the nonlinearities on the behavior of sculptured microstructures, taking advantage of instability-induced pattern transformation and piezoelectricity to actively tune the band gaps. .

FFT-Based Numerical Method for Nonlinear Elastic

In theoretical research pertaining to sealing, a contact model must be used to obtain the leakage channel. However, for elastoplastic contact, current numerical methods require a long calculation time. Hyperelastic contact is typically simplifed to a linear elastic contact problem, which must be improved in terms of calculation accuracy. Based on the fast Fourier transform, a numerical method suitable for elastoplastic and hyperelastic frictionless contact that can be used for solving two-dimensional and three-dimensional (3D) contact problems is proposed herein. The nonlinear elastic contact problem is converted into a linear elastic contact problem considering residual deformation (or the equivalent residual deformation). Results from numerical simulations for elastic, elastoplastic, and hyperelastic contact between a hemisphere and a rigid plane are compared with those obtained using the fnite element method to verify the accuracy of the numerical method. Compared with the existing elastoplastic contact numerical methods, the proposed method achieves a higher calculation efciency while ensuring a certain calculation accuracy (i.e., the pressure error does not exceed 15%, whereas the calculation time does not exceed 10 min in a 64 × 64 grid). For hyperelastic contact, the proposed method reduces the dependence of the approximation result on the load, as in a linear elastic approximation. Finally, using the sealing application as an example, the contact and leakage rates between complicated 3D rough surfaces are calculated. Despite a certain error, the simplifed numerical method yields a better approximation result than the linear elastic contact approximation. Additionally, the result can be used as fast solutions in engineering applications.

Application of MLPG in large deformation analysis

Two-dimensional large deformation analysis of hyperelastic and elasto-plastic solids based on the Meshless Local Petrov-Galerkin method （MLPG） is presented. A material configuration based the nonlinear MLPG formulation is introduced for the large deformation analysis of both path-dependent and path-independent materials. The supports of the MLS approximation functions cover the same sets of nodes during material deformation, thus the shape function needs to be computed only in the initial stage. The multiplicative hyperelasto-plastic constitutive model is adopted to avoid objective time integration for stress update in large rota- tion. With this constitutive model, the computational formulations for path-dependent and path-independent materials become identical. Computational efficiency of the nonlinear MLPG method is discussed and optimized in several aspects to make the MLPG an O（N） algorithm. The numerical examples indicate that the MLPG method can solve large deformation problems accurately. Moreover, the MLPG computations enjoy better convergence rate than the FEM under very large particle distortion.

T-Splines for Isogeometric Analysis of Two-Dimensional Nonlinear Problems

Nonlinear behaviors are commonplace in many complex engineering applications,e.g.,metal forming,vehicle crash test and so on.This paper focuses on the T-spline based isogeometric analysis of two-dimensional nonlinear problems including general large deformation hyperelastic problems and small deformation elastoplastic problems,to reveal the advantages of local refinement property of T-splines in describing nonlinear behavior of materials.By applying the adaptive refinement capability of T-splines during the iteration process of analysis,the numerical simulation accuracy of the nonlinear model could be increased dramatically.The Bézier extraction of the T-splines provides an element structure for isogeometric analysis that can be easily incorporated into existing nonlinear finite element codes.In addition,T-splines show great superiority of modeling complex geometries especially when the model is irregular and with hole features.Several numerical examples have been tested to validate the accuracy and convergence of the proposed method.The obtained results are compared with those from NURBS-based isogeometric analysis and commercial software ABAQUS.

Constitutive Modelling of Elastomeric Seal Material under Compressive Loading

Elastomers are used in numerous engineering applications such as sealing components, it is therefore important to devise a method that can accurately predict elastomers’ response to load. Many applications that employ the use of these materials subject them to a nonlinear large strain;therefore the simple Hooke’s law is not sufficient to describe their material behaviour. This paper presents an approach to obtain material properties of elastomer under compression loading, based on hyperelastic strain formulation, through experimental test and finite element modelling. The paper focuses on the isotropic incompressible behaviour exhibited by elastomers, and obtains strain energy functions that satisfy the characteristic properties of a hyperelastic model. Data obtained from compression test on a nitrile rubber (NBR) specimen were used as material input into ABAQUS®—a finite element analysis software. A least square fitting technique was used to determine the coefficients of various stable hyperelastic models, based on Drucker’s stability criteria within the software. The strain energy functions obtained concentrate on material parameters which are related to physical quantities of the material molecular network they are subjected to in practical application. The approach benefits from mathematical simplicity, and possesses the property of the deformation mode dependency. Furthermore, a model validation procedure using a step-by-step method for parameters estimation is explained. The work herein is a nonlinear finite element modelling process that leads to an optimal solution and can be employed not only for elastomeric seals, but also for similar engineering assets.

A 3-D Visco-Hyperelastic Constitutive Model for Rubber with Damage for Finite Element Simulation

A constitutive model to describe the behavior of rubber from low to high strain rates is presented.For loading,the primary hyperelastic behavior is characterized by the six parameter Ogden’s strain-energy potential of the third order.The rate-dependence is captured by the nonlinear second order BKZ model using another five parameters,having two relaxation times.For unloading,a single parameter model has been presented to define Hysteresis or continuous damage,while Ogden’s two term model has been used to capture Mullin’s effect or discontinuous damage.Lastly,the Feng-Hallquist failure surface dictates the ultimate failure for element deletion.The proposed model can accurately predict the response of rubber using a limited set of experimental data.The model has been validated here for the case of rubber but can be extended to a wide range of polymers.

Influence of random uncertainties of anisotropic fibrous model parameters on arterial pressure estimation

This paper deals with a stochastic approach based on the principle of the maximum entropy to investigate the effect of the parameter random uncertainties on the arterial pressure. Motivated by a hyperelastic, anisotropic, and incompressible constitutive law with fiber families, the uncertain parameters describing the mechanical behavior are considered. Based on the available information, the probability density functions are attributed to every random variable to describe the dispersion of the model parameters. Numerous realizations are carried out, and the corresponding arterial pressure results are compared with the human non-invasive clinical data recorded over a mean cardiac cycle. Furthermore, the Monte Carlo simulations are performed, the convergence of the probabilistic model is proven. The different realizations are useful to define a reliable confidence region, in which the probability to have a realization is equM to 95%. It is shown through the obtained results that the error in the estimation of the arterial pressure can reach 35% when the estimation of the model parameters is subjected to an uncertainty ratio of 5%. Finally, a sensitivity analysis is performed to identify the constitutive law relevant parameters for better understanding and characterization of the arterial wall mechanical behaviors.

Bubble-Enriched Smoothed Finite Element Methods for Nearly-Incompressible Solids

This work presents a locking-free smoothed finite element method(S-FEM)for the simulation of soft matter modelled by the equations of quasi-incompressible hyperelasticity.The proposed method overcomes well-known issues of standard finite element methods(FEM)in the incompressible limit:the over-estimation of stiffness and sensitivity to severely distorted meshes.The concepts of cell-based,edge-based and node-based S-FEMs are extended in this paper to three-dimensions.Additionally,a cubic bubble function is utilized to improve accuracy and stability.For the bubble function,an additional displacement degree of freedom is added at the centroid of the element.Several numerical studies are performed demonstrating the stability and validity of the proposed approach.The obtained results are compared with standard FEM and with analytical solutions to show the effectiveness of the method.

Modeling of micropipette aspiration and optical tweezers stretching of erythrocytes with or without Malaria parasite

The erythrocytes play an important role in the human body. The healthy erythrocytes can undergo extremely large deformation while passing through small capillaries. Their infection by Malaria Plasmodium falcipurum （P.f.） will lead to capillary blockage and blood flow obstruction. Many experimental and computational methods have been applied to study the increase in stickiness and decrease in deformability of the Malaria （P.f.） infected erythrocytes. The novelty of this paper lies in the establishment of an multi-component model for investigating mechanical properties of Malaria （P.f.） infected erythrocytes, especially of their enclosed parasites. Finite element method was applied to simulate the erythrocytes＇ deformation in micropipette aspiration and optical tweezers stretching using the computational software ABAQUS. The comparisons between simulations and experiments were able to quantitatively conclude the effects of stiffness and stickiness of the parasitophorous vacuole membrane on the cells＇ deformation, which could not be obtained from experiments directly.

An improved dynamic model for a silicone material beam with large deformation

Dynamic modeling for incompressible hyperelastic materials with large deformation is an important issue in biomimetic applications. The previously proposed lower-order fully parameterized absolute nodal coordinate formulation（ANCF） beam element employs cubic interpolation in the longitudinal direction and linear interpolation in the transverse direction, whereas it cannot accurately describe the large bending deformation. On this account, a novel modeling method for studying the dynamic behavior of nonlinear materials is proposed in this paper. In this formulation, a higher-order beam element characterized by quadratic interpolation in the transverse directions is used in this investigation. Based on the Yeoh model and volumetric energy penalty function, the nonlinear elastic force matrices are derived within the ANCF framework. The feasibility and availability of the Yeoh model are verified through static experiment of nonlinear incompressible materials. Furthermore,dynamic simulation of a silicone cantilever beam under the gravity force is implemented to validate the superiority of the higher-order beam element. The simulation results obtained based on the Yeoh model by employing three different ANCF beam elements are compared with the result achieved from a commercial finite element package as the reference result. It is found that the results acquired utilizing a higher-order beam element are in good agreement with the reference results,while the results obtained using a lower-order beam element are different from the reference results. In addition, the stiffening problem caused by volumetric locking can be resolved effectively by applying a higher-order beam element. It is concluded that the proposed higher-order beam element formulation has satisfying accuracy in simulating dynamic motion process of the silicone beam.

NONLINEAR TRAVELING WAVES IN A COMPRESSIBLE MOONEY-RIVLIN ROD I.LONG FINITE-AMPLITUDE WAVES

In literature,nonlinear traveling waves in elastic circular rods have only been studied based on single partial differential equation(pde)models,and here we consider such a problem by using a more accurate coupled-pde model.We derive the Hamiltonian from the model equations for the long finite-amplitude wave approximation,analyze how the number of singular points of the system changes with the parameters,and study the features of these singular points qualitatively.Various physically acceptable nonlinear traveling waves are also discussed,and corresponding examples are given.In particular,we find that certain waves,which cannot be counted by the single-equation model,can arise.

CAVITATED BIFURCATION FOR INCOMPRESSIBLE HYPERELASTIC MATERIAL

The spherical cavitated bifurcation for a hyperelastic solid sphere made of the incompressible Valanis-Landel material under boundary dead-loading is examined. The analytic solution for the bifurcation problem is obtained. The catastrophe and concentration of stresses are discussed. The stability of solutions is discussed through the energy comparison. And the growth of a pre-existing micro-void is also observed.

A comparative study of 85 hyperelastic constitutive models for both unfilled rubber and highly filled rubber nanocomposite material

Nonlinear finite element analysis is widely used for structural optimization of the design and the reliability analysis of complex elastomeric components.However,high-precision numerical results cannot be achieved without reliable strain energy functions(SEFs)of the rubber or rubber nanocomposite material.Although hyperelastic constitutive models have been studied for nearly 80 years,selecting one that accurately describes rubber's mechanical response is still a challenge.This work reviews 85 isotropic SEFs based on both the phenomenological theory and the micromechanical network theory proposed from the 1940s to 2019.A fitting algorithm which can realize the automatic fitting optimization and determination of the parameters of all SEFs reviewed is developed.The ability of each SEF to reproduce the experimental data of both the unfilled and highly filled rubber nanocomposite is quantitatively assessed based on a new proposed evaluation index.The top 30 SEFs for the unfilled rubber and the top 14 SEFs for the highly filled rubber nanocomposite are presented in the ranking lists.Finally,some suggestions on how to select an appropriate hyperelastic constitutive model are given,and the perspective on the future progress of constitutive models is summarized.