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 wit...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.展开更多
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 obta...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.展开更多
The nonlinear vibration problem is studied for a thin-walled rubber cylindrical shell composed of the classical incompressible Mooney-Rivlin material and subjected to a radial harmonic excitation. With the KirchhofF-L...The nonlinear vibration problem is studied for a thin-walled rubber cylindrical shell composed of the classical incompressible Mooney-Rivlin material and subjected to a radial harmonic excitation. With the KirchhofF-Love hypothesis, DonnelFs nonlinear shallow shell theory, hyperelastic constitutive relation, Lagrange equations and small strain hypothesis, a system of nonlinear differential equations describing the large-deflection vibration of the shell is derived. First, the natural frequencies of radial, circumferential and axial vibrations axe studied. Then, based on the bifurcation diagrams and the Poincare sections, the nonlinear behaviors describing the radial vibration of the shell are illustrated. Examining the influences of structural and material parameters on radial vibration of the shell shows that the vibration modes are highly sensitive to the thickness-radius ratio when the ratio is less than a certain critical value. Moreover, in terms of the results of multimodal expansion, it is found that the response of the shell to radial motion is more regular than that without considering the coupling between modes, while there are more phenomena for the uncoupled case.展开更多
The partial differential equation for isotropic hyperelastic constitutive models has been postulated and derived from the balance between stored energy and stress work done. The partial differential equation as a func...The partial differential equation for isotropic hyperelastic constitutive models has been postulated and derived from the balance between stored energy and stress work done. The partial differential equation as a function of three invariants has then been solved by Lie group methods. With geometric meanings of deformations, the general solution boils down to a particular three-term solution. The particular solution has been applied for several isotropic hyperelastic materials. For incompressible materials, vulcanized rubber containing 8% sulfur and Entec Enflex S4035A thermoplastic elastomer, three coefficients have been determined from uniaxial tension data and applied to predict the pure shear and equibiaxial tension modes. For a slightly compressible rubber material, the coefficients have also been extracted from the confined volumetric test data.展开更多
Primary objective of automobile seats is to offer adequate level of safety and comfort to the seated human occupant, primarily against vibration. Ideally, any sort of automotive seat is constructed by mechanical frame...Primary objective of automobile seats is to offer adequate level of safety and comfort to the seated human occupant, primarily against vibration. Ideally, any sort of automotive seat is constructed by mechanical framework, cushion, backrest and headrest. The frame structures are made of metallic alloys, while the cushion, backrest and headrest are made of polyurethane foam material. During the design phase of automotive seat, the greatest challenge is to assign realistic material properties to foam material;as it is non-linear in nature and exhibit hysteresis at low level stress. In this research paper, a car seat has been modelled in finite element environment by implementing both hyperelastic and viscoelastic material properties to polyurethane foam. The car seat has been excited with the loads due to car acceleration and human object and the effects of vibration in terms of vertical acceleration at different locations have been measured. The aims of this simulation study are to establish a car seat with the foam material properties as accurately as possible and provide a finite element set up of car seat to monitor the vertical acceleration responses in a reasonable way. The RMS acceleration values for headrest, backrest and cushion have been found to be 0.91 mm/sec2, 0.54 mm/sec2 and 0.47 mm/sec2, respectively, which showed that the car seat foam can effectively be modelled through combined hyperelastic and viscoelastic material formulations. The simulation outputs have been validated through real life testing data, which clearly indicates that this computerized simulation technique is capable of anticipating the acceleration responses at different car seat segments in a justified way.展开更多
The e?ect of material anisotropy on the bifurcation for void formation in anisotropic compressible hyperelastic materials is examined. Numerical solutions are obtained in an anisotropic sphere, whose material is tra...The e?ect of material anisotropy on the bifurcation for void formation in anisotropic compressible hyperelastic materials is examined. Numerical solutions are obtained in an anisotropic sphere, whose material is transversely isotropic in the radial direction. It is shown that the bifur- cation may occur either to the right or to the left, depending on the degree of material anisotropy. The deformation and stress contribution in the sphere before cavitation are di?erent from those after cavitation. The stability of solutions is discussed through a comparison of energy.展开更多
Bifurcation problems both for static and dynamical cavitation in a solid sphere composed of the incompressible hyperelasticplastic material, under uniformly distributed tensile boundary dead load were studied. For eac...Bifurcation problems both for static and dynamical cavitation in a solid sphere composed of the incompressible hyperelasticplastic material, under uniformly distributed tensile boundary dead load were studied. For each problem, cavity forms at the center of the sphere when the tensile load is larger than its critical value. Bifurcation curves and the growth curves for the plastic deformation region were given. For static cavitation, the deformation displays three stages, namely, fully elastic, elasto-plastic and fully plastic stages. For dynamical cavitation, the cavity grows without bound and the sphere displays plastic flow.展开更多
Rubber-like materials that are commonly used in structural applications are modelled using hyperelastic material models.Most of the hyperelastic materials are nearly incompressible,which poses challenges,i.e.,volumetr...Rubber-like materials that are commonly used in structural applications are modelled using hyperelastic material models.Most of the hyperelastic materials are nearly incompressible,which poses challenges,i.e.,volumetric locking during numerical modelling.There exist many formulations in the context of the finite element method,among which the mixed displacementpressure formulation is robust.However,such a displacement-pressure formulation is less explored in meshfree methods,which mitigates the problem associated with mesh distortion during large deformation.This work addresses this issue of alleviating volumetric locking in the element-free Galerkin method(EFGM),which is one of the popular meshfree methods.A two-field mixed variational formulation using the perturbed Lagrangian approach within the EFGM framework is proposed for modelling nearly incompressible hyperelastic material models,such as Neo-Hookean and Mooney-Rivlin.Taking advantage of the meshless nature of the EFGM,this work introduces a unique approach by randomly distributing pressure nodes across the geometry,following specific guidelines.A wide spectrum of problems involving bending,tension,compression,and contact is solved using two approaches of the proposed displacement-pressure node formulation involving regular and irregular pressure node distribution.It is observed that both approaches give accurate results compared to the reference results,though the latter offers flexibility in the pressure nodal distribution.展开更多
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 for...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.展开更多
In this paper, the dynamic characteristics are examined for a cylindrical membrane composed of a transversely isotropic incompressible hyperelastic material under an applied uniform radial constant pressure at its inn...In this paper, the dynamic characteristics are examined for a cylindrical membrane composed of a transversely isotropic incompressible hyperelastic material under an applied uniform radial constant pressure at its inner surface. A second-order nonlinear ordinary differential equation that approximately describes the radial oscillation of the inner surface of the membrane with respect to time is obtained. Some interesting conclusions are proposed for different materials, such as the neo-Hookean material, the Mooney-Rivlin material and the Rivlin-Saunders material. Firstly, the bifurcation conditions depending on the material parameters and the pressure loads are determined. Secondly, the conditions of periodic motion are presented in detail for membranes composed of different materials. Meanwhile, numerical simulations are also provided.展开更多
Torsional instability of an incompressible thermo-hyperelastic cylindrical rod, subjected to axial stretching and large torsions, is examined within the framework of finite elasticity. When the cylinder is stretched a...Torsional instability of an incompressible thermo-hyperelastic cylindrical rod, subjected to axial stretching and large torsions, is examined within the framework of finite elasticity. When the cylinder is stretched and twisted by a sufficiently large degree, a knot may form suddenly at one point. This inherent elastic instability is analyzed with the minimum potential energy principle and the critical values of torsion are obtained. The distribution of stresses as well as the tensile force and the torque are studied. Effect of tem- perature change is specifically discussed.展开更多
The present study is focused on the constitutive modeling for the mechanical behavior of rubber reinforced with filler particles.A filler-dependent energy density function is proposed with all the continuum mechanics-...The present study is focused on the constitutive modeling for the mechanical behavior of rubber reinforced with filler particles.A filler-dependent energy density function is proposed with all the continuum mechanics-based necessities of an effective hyperelastic material model.The proposed invariant-based energy function comprises a single set of material parameters for a material subjected to several modes of loading conditions.The model solution agrees well with existing experimental results.Later,the effect of varying concentrations of filler particles in the rubber matrix is also studied.展开更多
A molecular network constitutive model is proposed in this study.Based on the concept of molecular network decomposition,the molecular network was decomposed into a James-Guth 3-chain network model and an Arruda-Boyce...A molecular network constitutive model is proposed in this study.Based on the concept of molecular network decomposition,the molecular network was decomposed into a James-Guth 3-chain network model and an Arruda-Boyce 8-chain network model.Considering that the single molecular chain of rubber is constrained by the surrounding molecular networks, the 3-chain model and the 8-chain model were consequently modified using the tube theory.The proposed model contains four material parameters that were determined by fitting the uniaxial tensile test data from rubber-like materials,and the values of these parameters were utilized to determine the pure shear and the biaxial extension.The proposed model agrees well with the experimental data and can predict the test results of experiments performed under the conditions of pure shear and equi-biaxial deformation with acceptable accuracy.展开更多
基金the National Key Research and Development Program of China(2018YFB1502501)the National Natural Science Foundation of China(52003024).
文摘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.
文摘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.
基金supported by the National Natural Science Foundation of China (Nos.11672069,11702059,11872145).
文摘The nonlinear vibration problem is studied for a thin-walled rubber cylindrical shell composed of the classical incompressible Mooney-Rivlin material and subjected to a radial harmonic excitation. With the KirchhofF-Love hypothesis, DonnelFs nonlinear shallow shell theory, hyperelastic constitutive relation, Lagrange equations and small strain hypothesis, a system of nonlinear differential equations describing the large-deflection vibration of the shell is derived. First, the natural frequencies of radial, circumferential and axial vibrations axe studied. Then, based on the bifurcation diagrams and the Poincare sections, the nonlinear behaviors describing the radial vibration of the shell are illustrated. Examining the influences of structural and material parameters on radial vibration of the shell shows that the vibration modes are highly sensitive to the thickness-radius ratio when the ratio is less than a certain critical value. Moreover, in terms of the results of multimodal expansion, it is found that the response of the shell to radial motion is more regular than that without considering the coupling between modes, while there are more phenomena for the uncoupled case.
文摘The partial differential equation for isotropic hyperelastic constitutive models has been postulated and derived from the balance between stored energy and stress work done. The partial differential equation as a function of three invariants has then been solved by Lie group methods. With geometric meanings of deformations, the general solution boils down to a particular three-term solution. The particular solution has been applied for several isotropic hyperelastic materials. For incompressible materials, vulcanized rubber containing 8% sulfur and Entec Enflex S4035A thermoplastic elastomer, three coefficients have been determined from uniaxial tension data and applied to predict the pure shear and equibiaxial tension modes. For a slightly compressible rubber material, the coefficients have also been extracted from the confined volumetric test data.
文摘Primary objective of automobile seats is to offer adequate level of safety and comfort to the seated human occupant, primarily against vibration. Ideally, any sort of automotive seat is constructed by mechanical framework, cushion, backrest and headrest. The frame structures are made of metallic alloys, while the cushion, backrest and headrest are made of polyurethane foam material. During the design phase of automotive seat, the greatest challenge is to assign realistic material properties to foam material;as it is non-linear in nature and exhibit hysteresis at low level stress. In this research paper, a car seat has been modelled in finite element environment by implementing both hyperelastic and viscoelastic material properties to polyurethane foam. The car seat has been excited with the loads due to car acceleration and human object and the effects of vibration in terms of vertical acceleration at different locations have been measured. The aims of this simulation study are to establish a car seat with the foam material properties as accurately as possible and provide a finite element set up of car seat to monitor the vertical acceleration responses in a reasonable way. The RMS acceleration values for headrest, backrest and cushion have been found to be 0.91 mm/sec2, 0.54 mm/sec2 and 0.47 mm/sec2, respectively, which showed that the car seat foam can effectively be modelled through combined hyperelastic and viscoelastic material formulations. The simulation outputs have been validated through real life testing data, which clearly indicates that this computerized simulation technique is capable of anticipating the acceleration responses at different car seat segments in a justified way.
基金Project supported by the National Natural Science Foundation of China (No.10272069) and the Shanghai Key SubjectProgram.
文摘The e?ect of material anisotropy on the bifurcation for void formation in anisotropic compressible hyperelastic materials is examined. Numerical solutions are obtained in an anisotropic sphere, whose material is transversely isotropic in the radial direction. It is shown that the bifur- cation may occur either to the right or to the left, depending on the degree of material anisotropy. The deformation and stress contribution in the sphere before cavitation are di?erent from those after cavitation. The stability of solutions is discussed through a comparison of energy.
基金Project supported by National Natural Science Foundation of China （Grant Nos. 10402018, 10272069）
文摘Bifurcation problems both for static and dynamical cavitation in a solid sphere composed of the incompressible hyperelasticplastic material, under uniformly distributed tensile boundary dead load were studied. For each problem, cavity forms at the center of the sphere when the tensile load is larger than its critical value. Bifurcation curves and the growth curves for the plastic deformation region were given. For static cavitation, the deformation displays three stages, namely, fully elastic, elasto-plastic and fully plastic stages. For dynamical cavitation, the cavity grows without bound and the sphere displays plastic flow.
基金supported by the DST-SERB and VSSC,ISRO of the project titled“Functionality Enhancement through Design and Development of Advanced Finite Element Algorithms for STR tools”under IMPRINT.IIC(IMP/2019/000276)scheme.
文摘Rubber-like materials that are commonly used in structural applications are modelled using hyperelastic material models.Most of the hyperelastic materials are nearly incompressible,which poses challenges,i.e.,volumetric locking during numerical modelling.There exist many formulations in the context of the finite element method,among which the mixed displacementpressure formulation is robust.However,such a displacement-pressure formulation is less explored in meshfree methods,which mitigates the problem associated with mesh distortion during large deformation.This work addresses this issue of alleviating volumetric locking in the element-free Galerkin method(EFGM),which is one of the popular meshfree methods.A two-field mixed variational formulation using the perturbed Lagrangian approach within the EFGM framework is proposed for modelling nearly incompressible hyperelastic material models,such as Neo-Hookean and Mooney-Rivlin.Taking advantage of the meshless nature of the EFGM,this work introduces a unique approach by randomly distributing pressure nodes across the geometry,following specific guidelines.A wide spectrum of problems involving bending,tension,compression,and contact is solved using two approaches of the proposed displacement-pressure node formulation involving regular and irregular pressure node distribution.It is observed that both approaches give accurate results compared to the reference results,though the latter offers flexibility in the pressure nodal distribution.
基金supported by the National Natural Science Foundation of China (11772186 and 11272203)
文摘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.
基金Project supported by the National Natural Science Foundation of China (Nos. 10872045 and 10772104)the Program for New Century Excellent Talents in University (No. NCET-09-0096)the Fundamental Research Funds for the Central Universities (No. DC10030104)
文摘In this paper, the dynamic characteristics are examined for a cylindrical membrane composed of a transversely isotropic incompressible hyperelastic material under an applied uniform radial constant pressure at its inner surface. A second-order nonlinear ordinary differential equation that approximately describes the radial oscillation of the inner surface of the membrane with respect to time is obtained. Some interesting conclusions are proposed for different materials, such as the neo-Hookean material, the Mooney-Rivlin material and the Rivlin-Saunders material. Firstly, the bifurcation conditions depending on the material parameters and the pressure loads are determined. Secondly, the conditions of periodic motion are presented in detail for membranes composed of different materials. Meanwhile, numerical simulations are also provided.
基金The project supported by the National Natural Science Foundation of China （10402018, 10272069） and Shanghai Key Project Program （Y0103） The English text was polished by Keren Wang.
文摘Torsional instability of an incompressible thermo-hyperelastic cylindrical rod, subjected to axial stretching and large torsions, is examined within the framework of finite elasticity. When the cylinder is stretched and twisted by a sufficiently large degree, a knot may form suddenly at one point. This inherent elastic instability is analyzed with the minimum potential energy principle and the critical values of torsion are obtained. The distribution of stresses as well as the tensile force and the torque are studied. Effect of tem- perature change is specifically discussed.
文摘The present study is focused on the constitutive modeling for the mechanical behavior of rubber reinforced with filler particles.A filler-dependent energy density function is proposed with all the continuum mechanics-based necessities of an effective hyperelastic material model.The proposed invariant-based energy function comprises a single set of material parameters for a material subjected to several modes of loading conditions.The model solution agrees well with existing experimental results.Later,the effect of varying concentrations of filler particles in the rubber matrix is also studied.
基金National Natural Science Foundation of China under Grant No. 11372074.
文摘A molecular network constitutive model is proposed in this study.Based on the concept of molecular network decomposition,the molecular network was decomposed into a James-Guth 3-chain network model and an Arruda-Boyce 8-chain network model.Considering that the single molecular chain of rubber is constrained by the surrounding molecular networks, the 3-chain model and the 8-chain model were consequently modified using the tube theory.The proposed model contains four material parameters that were determined by fitting the uniaxial tensile test data from rubber-like materials,and the values of these parameters were utilized to determine the pure shear and the biaxial extension.The proposed model agrees well with the experimental data and can predict the test results of experiments performed under the conditions of pure shear and equi-biaxial deformation with acceptable accuracy.
