The cavitated bifurcation problem in a solid sphere composed oftwo compressible hyper-elas- tic materials is examined. Thebifurcation solution for the composed sphere under a uniform radialtensile boundary dead-load i...The cavitated bifurcation problem in a solid sphere composed oftwo compressible hyper-elas- tic materials is examined. Thebifurcation solution for the composed sphere under a uniform radialtensile boundary dead-load is obtained. The bifurcation curves andthe stress contributions subsequent to the cavita- tion are given.The right and left bifurcation as well as the catastrophe andconcentration of stresses are ana- lyzed. The stability of solutionsis discussed through an energy comparison.展开更多
In this paper, the strain energy function proposed by Shang and Cheng was generalized by introducing a nonlinear term. Void formation and growth in the interior of a sphere composed of compressible hyper-elastic mater...In this paper, the strain energy function proposed by Shang and Cheng was generalized by introducing a nonlinear term. Void formation and growth in the interior of a sphere composed of compressible hyper-elastic material, subjected to a prescribed uniform displacement, was examined. A parametric cavitated bifurcation solution for the radial deformed function was obtained. Stability of the solution of the cavitated bifurcation equation was discussed. With the appearance of a cavity, an interesting feature of the radial deformation near the deformed cavity wall is the transition from extension to compression.展开更多
Dynamical responses, such as motion and destruction of hyper-elastic cylindrical shells subject to periodic or suddenly applied constant load on the inner surface, are studied within a framework of finite elasto-dynam...Dynamical responses, such as motion and destruction of hyper-elastic cylindrical shells subject to periodic or suddenly applied constant load on the inner surface, are studied within a framework of finite elasto-dynamics. By numerical computation and dynamic qualitative analysis of the nonlinear differential equation, it is shown that there exists a certain critical value for the internal load describing motion of the inner surface of the shell. Motion of the shell is nonlinear periodic or quasi-periodic oscillation when the average load of the periodic load or the constant load is less than its critical value. However, the shell will be destroyed when the load exceeds the critical value. Solution to the static equilibrium problem is a fixed point for the dynamical response of the corresponding system under a suddenly applied constant load. The property of fixed point is related to the property of the dynamical solution and motion of the shell. The effects of thickness and load parameters on the critical value and oscillation of the shell are discussed.展开更多
An efficient hyper-elastic model that can reflect the primary mechanical behaviors of Carbon-Kevlar hybrid woven reinforcement was developed and implemented with VUMAT constitutive code for preforming simulation.The m...An efficient hyper-elastic model that can reflect the primary mechanical behaviors of Carbon-Kevlar hybrid woven reinforcement was developed and implemented with VUMAT constitutive code for preforming simulation.The model parameters were accurately determined through the uniaxial and bias-extension tests.To calibrate the simulation code,preforming experiments of hybrid woven reinforcement over the hemisphere mold and tetrahedron mold were respectively conducted to validate the proposed hyper-elastic model.The comparison between the simulations and experiments shows that the model can not only accurately capture shear angle distribution and geometry shape after deformation,but also accurately predict the force–displacement curve and potential fiber tensile failure during the preforming process.This result indicates that the proposed model can be used to predict the preforming behavior of Carbon-Kevlar hybrid woven reinforcement,and simulate its manufacturing process of complicated geometry.展开更多
The problems on the non-uniqueness and stability of a two-family fiber- reinforced anisotropic incompressible hyper-elastic square sheet under equibiaxial tensile dead loading are examined within the framework of fini...The problems on the non-uniqueness and stability of a two-family fiber- reinforced anisotropic incompressible hyper-elastic square sheet under equibiaxial tensile dead loading are examined within the framework of finite elasticity. For a two-family fiber-reinforced square sheet, which is in-plane symmetric and subjected to the in-plane symmetric tension in dead loading on the edges, three symmetrically deformed configu- rations and six asymmetrically deformed configurations are possible for any values of the loading. Moreover, another four bifurcated asymmetrically deformed configurations are possible for the loading beyond a certain criticM value. The stability of all the solutions is discussed in comparison with the energy of the sheet. It is shown that only one of the symmetric solutions is stable when the loading is less than the critical value. However, this symmetric solution will become unstable when the loading is larger than the critical value, while one of the four bifurcated asymmetric solutions will be stable.展开更多
The deformations and stresses of a rotating cylindrical hollow disk made of incompressible functionally-graded hyper-elastic material are theoretically analyzed based on the finite elasticity theory.The hyper-elastic ...The deformations and stresses of a rotating cylindrical hollow disk made of incompressible functionally-graded hyper-elastic material are theoretically analyzed based on the finite elasticity theory.The hyper-elastic material is described by a new micro-macro transition model.Specially,the material shear modulus and density are assumed to be a function with a power law form through the radial direction,while the material inhomogeneity is thus reflected on the power index m.The integral forms of the stretches and stress components are obtained.With the obtained complicated integral forms,the composite trapezoidal rule is utilized to derive the analytical solutions,and the explicit solutions for both the stretches and the stress components are numerically obtained.By comparing the results with two classic models,the superiority of the model in our work is demonstrated.Then,the distributions of the stretches and normalized stress components are discussed in detail under the effects of m.The results indicate that the material inhomogeneity and the rotating angular velocity have significant effects on the distributions of the normalized radial and hoop stress components and the stretches.We believe that by appropriately choosing the material inhomogeneity and configuration parameters,the functionally-graded material(FGM)hyper-elastic hollow cylindrical disk can be designed to meet some unique requirements in the application fields,e.g.,soft robotics,medical devices,and conventional aerospace and mechanical industries.展开更多
The dynamical formation of cavity in a hyper_elastic sphere composed of two materials with the incompressible strain energy function, subjected to a suddenly applied uniform radial tensile boundary dead_load, was stud...The dynamical formation of cavity in a hyper_elastic sphere composed of two materials with the incompressible strain energy function, subjected to a suddenly applied uniform radial tensile boundary dead_load, was studied following the theory of finite deformation dynamics. Besides a trivial solution corresponding to the homogeneous static state, a cavity forms at the center of the sphere when the tensile load is larger than its critical value. An exact differential relation between the cavity radius and the tensile land was obtained. It is proved that the evolution of cavity radius with time displays nonlinear periodic oscillations. The phase diagram for oscillation, the maximum amplitude, the approximate period and the critical load were all discussed.展开更多
The growth of a prolate or oblate elliptic micro-void in a fiber reinforced anisotropic incompressible hyper-elastic rectangular thin plate subjected to uniaxial extensions is studied within the framework of finite el...The growth of a prolate or oblate elliptic micro-void in a fiber reinforced anisotropic incompressible hyper-elastic rectangular thin plate subjected to uniaxial extensions is studied within the framework of finite elasticity. Coupling effects of void shape and void size on the growth of the void are paid special attention to. The deformation function of the plate with an isolated elliptic void is given, which is expressed by two parameters to solve the differential equation. The solution is approximately obtained from the minimum potential energy principle. Deformation curves for the void with a wide range of void aspect ratios and the stress distributions on the surface of the void have been obtained by numerical computation. The growth behavior of the void and the characteristics of stress distributions on the surface of the void are captured. The combined effects of void size and void shape on the growth of the void in the thin plate are discussed. The maximum stresses for the void with different sizes and different void aspect ratios are compared.展开更多
We construct an unconventional divergence preserving discretization of updated Lagrangian ideal magnetohydrodynamics(MHD)over simplicial grids.The cell-centered finite-volume(FV)method employed to discretize the conse...We construct an unconventional divergence preserving discretization of updated Lagrangian ideal magnetohydrodynamics(MHD)over simplicial grids.The cell-centered finite-volume(FV)method employed to discretize the conservation laws of volume,momentum,and total energy is rigorously the same as the one developed to simulate hyperelasticity equations.By construction this moving mesh method ensures the compatibility between the mesh displacement and the approximation of the volume flux by means of the nodal velocity and the attached unit corner normal vector which is nothing but the partial derivative of the cell volume with respect to the node coordinate under consideration.This is precisely the definition of the compatibility with the Geometrical Conservation Law which is the cornerstone of any proper multi-dimensional moving mesh FV discretization.The momentum and the total energy fluxes are approximated utilizing the partition of cell faces into sub-faces and the concept of sub-face force which is the traction force attached to each sub-face impinging at a node.We observe that the time evolution of the magnetic field might be simply expressed in terms of the deformation gradient which characterizes the Lagrange-to-Euler mapping.In this framework,the divergence of the magnetic field is conserved with respect to time thanks to the Piola formula.Therefore,we solve the fully compatible updated Lagrangian discretization of the deformation gradient tensor for updating in a simple manner the cell-centered value of the magnetic field.Finally,the sub-face traction force is expressed in terms of the nodal velocity to ensure a semi-discrete entropy inequality within each cell.The conservation of momentum and total energy is recovered prescribing the balance of all the sub-face forces attached to the sub-faces impinging at a given node.This balance corresponds to a vectorial system satisfied by the nodal velocity.It always admits a unique solution which provides the nodal velocity.The robustness and the accuracy of this unconventional FV scheme have been demonstrated by employing various representative test cases.Finally,it is worth emphasizing that once you have an updated Lagrangian code for solving hyperelasticity you also get an almost free updated Lagrangian code for solving ideal MHD ensuring exactly the compatibility with the involution constraint for the magnetic field at the discrete level.展开更多
An explicit integration scheme for rate-dependent crystal plasticity （CP） was developed. Additive decomposition of the velocity gradient tensor into lattice and plastic parts is adopted for describing the kinematics...An explicit integration scheme for rate-dependent crystal plasticity （CP） was developed. Additive decomposition of the velocity gradient tensor into lattice and plastic parts is adopted for describing the kinematics; the Cauchy stress is calculated by using a hypo-elastic formulation, applying the Jaumann stress rate. This CP scheme has been implemented into a commercial finite element code （CPFEM）. Uniaxial compression and roiling processes were simulated. The results show good accuracy and reliability of the integration scheme. The results were compared with simulations using one hyper-elastic CPFEM implementation which involves multiplicative decomposition of the deformation gradient tensor. It is found that the hypo-elastic implementation is only slightly faster and has a similar accuracy as the hyper-elastic formulation.展开更多
Various mixed formulations of the finite element method (FEM) yield matrix equations involving zero diagonal entries. They are then dealt with by a penaltymethod so that they become non-zero but near zero terms. Howev...Various mixed formulations of the finite element method (FEM) yield matrix equations involving zero diagonal entries. They are then dealt with by a penaltymethod so that they become non-zero but near zero terms. However, the penalty has tobe chosen properly. If it is too large, the matrix equation may become ill-conditioned. Onthe other hand, the matrix equation may give incorrect answer if the penalty is too small.In non-linear regime, the difficulty is more serious because the magnitude order of the matrix varies considerably in the entire loading history. The paper suggests an iteration solution and applies it to non-linear FEM of rubber-like hyper-elasticity. This type of analysisis highly non-linear both in physics and in geometry as well as the strong constraint of incompressibility. The iteration solution is demonstrated to possess super precision and excellent convergence characteristics.展开更多
Because of the limited space of the launch rockets, deployable mechanisms are always used to solve the phenomenon. One dimensional deployable mast can deploy and support antenna, solar sail and space optical camera. T...Because of the limited space of the launch rockets, deployable mechanisms are always used to solve the phenomenon. One dimensional deployable mast can deploy and support antenna, solar sail and space optical camera. Tape-spring hyperelastic hinges can be folded and extended into a rod like configuration. It utilizes the strain energy to realize self-deploying and drive the other structures. One kind of triangular prism mast with tape-spring hyperelastic hinges is proposed and developed. Stretching and compression stiffness theoretical model are established with considering the tape-spring hyperelastic hinges based on static theory. The finite element model of ten-module triangular prism mast is set up by ABAQUS with the tape-spring hyperelastic hinge and parameter study is performed to investigate the influence of thickness, section angle and radius. Two-module TPM is processed and tested the compression stiffness by the laser displacement sensor, deploying repeat accuracy by the high speed camera, modal shape and fundamental frequency at cantilever position by LMS multi-channel vibration test and analysis system, which are used to verify precision of the theoretical and finite element models of ten-module triangular prism mast with the tape-spring hyperelastic hinges. This research proposes an innovative one dimensional triangular prism with tape-spring hyperelastic hinge which has great application value to the space deployable mechanisms.展开更多
Mechanical properties, such as the deformation and stress distributions for venous walls under the combined load of transmural pressure and axial stretch, are examined within the framework of nonlinear elasticity with...Mechanical properties, such as the deformation and stress distributions for venous walls under the combined load of transmural pressure and axial stretch, are examined within the framework of nonlinear elasticity with one kind of hyper-elastic strain energy functions. The negative pressure instability problem of the venous wall is explained through energy comparison. First, the deformation equation of the venous wall under the combined loads is obtained with a thin-walled circular cylindrical tube. The deformation curves and the stress distributions for the venous wall are given under the normal transmural pressure, and the regulations are discussed. Then, the deformation curves of the venous wall under the negative transmural pressure or the internal pressure less than the external pressure are given. Finally, the negative pressure instability problem is discussed through energy comparison.展开更多
In the present work we describe a new expression of total stress taking into account the passive and active contributions but especially the pore level stress.Special attention is paid to the effective stress and osmo...In the present work we describe a new expression of total stress taking into account the passive and active contributions but especially the pore level stress.Special attention is paid to the effective stress and osmotic pressure gradient in numerical simulation to understand the mechanical behavior of the human arterial wall.The new model aims to predict the rigidity of the artery,by using the theoretical model of hyper-elastic,anisotropic and dynamical behavior of human common carotid artery.The principal obtained result showed that:the osmosis phenomenon is the best parameter to explain the loss water in arterial tissue.This loss of water causes the rigidity of the artery which thus can be controlled by the osmosis phenomenon.All the results are in good agreement with the expected results of the literature and could play the important role in the diagnosis of the patients with the CVD(Cardiovascular Disease).展开更多
Constitutive modeling is crucial for engineering design and simulations to accurately describe material behavior.However,traditional phenomenological models often struggle to capture the complexities of real materials...Constitutive modeling is crucial for engineering design and simulations to accurately describe material behavior.However,traditional phenomenological models often struggle to capture the complexities of real materials under varying stress conditions due to their fixed forms and limited parameters.While recent advances in deep learning have addressed some limitations of classical models,purely data-driven methods tend to require large data sets,lack interpretability,and struggle to generalize beyond their training data.To tackle these issues,we introduce“Fusion-based Constitutive model(FuCe):Toward model-data augmentation in constitutive modeling.”This approach combines established phenomenological models with an Input Convex Neural Network architecture,designed to train on the limited and noisy force-displacement data typically available in practical applications.The hybrid model inherently adheres to necessary constitutive conditions.During inference,Monte Carlo dropout is employed to generate Bayesian predictions,providing mean values and confidence intervals that quantify uncertainty.We demonstrate the model's effectiveness by learning two isotropic constitutive models and one anisotropic model with a single fiber direction,across six different stress states.The framework's applicability is also showcased in finite element simulations across three geometries of varying complexities.Our results highlight the framework's superior extrapolation capabilities,even when trained on limited and noisy data,delivering accurate and physically meaningful predictions across all numerical examples.展开更多
Dynamical cavitation and oscillation of an anisotropic two-family fiber-reinforced incompressible hyper-elastic sphere subjected to a suddenly applied constant boundary dead load are examined within the framework of f...Dynamical cavitation and oscillation of an anisotropic two-family fiber-reinforced incompressible hyper-elastic sphere subjected to a suddenly applied constant boundary dead load are examined within the framework of finite elasto-dynamics.An exact differential equation between the radius of the cavity and the applied load is obtained.The curves for the variation of the maximum radius of the cavity with the load and the phase diagrams are obtained by vibration theories and numerical computation.It is shown that there exists a critical value for the applied load.When the applied load is larger than the critical value,a spherical cavity will suddenly form at the center of the sphere.It is proved that the evolution of the cavity radius with time follows that of nonlinear periodic oscillation,and oscillation of the anisotropic sphere is not the same as that of the isotropic sphere.展开更多
基金the National Natttral Science Foundation of China(No.19802012)
文摘The cavitated bifurcation problem in a solid sphere composed oftwo compressible hyper-elas- tic materials is examined. Thebifurcation solution for the composed sphere under a uniform radialtensile boundary dead-load is obtained. The bifurcation curves andthe stress contributions subsequent to the cavita- tion are given.The right and left bifurcation as well as the catastrophe andconcentration of stresses are ana- lyzed. The stability of solutionsis discussed through an energy comparison.
文摘In this paper, the strain energy function proposed by Shang and Cheng was generalized by introducing a nonlinear term. Void formation and growth in the interior of a sphere composed of compressible hyper-elastic material, subjected to a prescribed uniform displacement, was examined. A parametric cavitated bifurcation solution for the radial deformed function was obtained. Stability of the solution of the cavitated bifurcation equation was discussed. With the appearance of a cavity, an interesting feature of the radial deformation near the deformed cavity wall is the transition from extension to compression.
基金the National Natural Science Foundation of China(Nos.10772104 and10402018)the Shanghai Leading Academic Discipline Project(No.Y0103)
文摘Dynamical responses, such as motion and destruction of hyper-elastic cylindrical shells subject to periodic or suddenly applied constant load on the inner surface, are studied within a framework of finite elasto-dynamics. By numerical computation and dynamic qualitative analysis of the nonlinear differential equation, it is shown that there exists a certain critical value for the internal load describing motion of the inner surface of the shell. Motion of the shell is nonlinear periodic or quasi-periodic oscillation when the average load of the periodic load or the constant load is less than its critical value. However, the shell will be destroyed when the load exceeds the critical value. Solution to the static equilibrium problem is a fixed point for the dynamical response of the corresponding system under a suddenly applied constant load. The property of fixed point is related to the property of the dynamical solution and motion of the shell. The effects of thickness and load parameters on the critical value and oscillation of the shell are discussed.
基金support from the Young Fund of Natural Science Foundation of Shaanxi Province,China(Nos.2020JQ-701 and 2020JQ-121)the Open Fund of State Key Laboratory of Applied Optics,China(No.SKLAO2020001A09).
文摘An efficient hyper-elastic model that can reflect the primary mechanical behaviors of Carbon-Kevlar hybrid woven reinforcement was developed and implemented with VUMAT constitutive code for preforming simulation.The model parameters were accurately determined through the uniaxial and bias-extension tests.To calibrate the simulation code,preforming experiments of hybrid woven reinforcement over the hemisphere mold and tetrahedron mold were respectively conducted to validate the proposed hyper-elastic model.The comparison between the simulations and experiments shows that the model can not only accurately capture shear angle distribution and geometry shape after deformation,but also accurately predict the force–displacement curve and potential fiber tensile failure during the preforming process.This result indicates that the proposed model can be used to predict the preforming behavior of Carbon-Kevlar hybrid woven reinforcement,and simulate its manufacturing process of complicated geometry.
基金supported by the National Natural Science Foundation of China(No.10772104)the Shanghai Leading Academic Discipline Project(No.S30106)
文摘The problems on the non-uniqueness and stability of a two-family fiber- reinforced anisotropic incompressible hyper-elastic square sheet under equibiaxial tensile dead loading are examined within the framework of finite elasticity. For a two-family fiber-reinforced square sheet, which is in-plane symmetric and subjected to the in-plane symmetric tension in dead loading on the edges, three symmetrically deformed configu- rations and six asymmetrically deformed configurations are possible for any values of the loading. Moreover, another four bifurcated asymmetrically deformed configurations are possible for the loading beyond a certain criticM value. The stability of all the solutions is discussed in comparison with the energy of the sheet. It is shown that only one of the symmetric solutions is stable when the loading is less than the critical value. However, this symmetric solution will become unstable when the loading is larger than the critical value, while one of the four bifurcated asymmetric solutions will be stable.
基金supported by the National Natural Science Foundation of China(No.11972144)the Shanxi Province Specialized Research and Development Breakthrough in Key Core and Generic Technologies(Key Research and Development Program)(No.2020XXX017)the Fundamental Research Program of Shanxi Province of China(No.202203021211134)。
文摘The deformations and stresses of a rotating cylindrical hollow disk made of incompressible functionally-graded hyper-elastic material are theoretically analyzed based on the finite elasticity theory.The hyper-elastic material is described by a new micro-macro transition model.Specially,the material shear modulus and density are assumed to be a function with a power law form through the radial direction,while the material inhomogeneity is thus reflected on the power index m.The integral forms of the stretches and stress components are obtained.With the obtained complicated integral forms,the composite trapezoidal rule is utilized to derive the analytical solutions,and the explicit solutions for both the stretches and the stress components are numerically obtained.By comparing the results with two classic models,the superiority of the model in our work is demonstrated.Then,the distributions of the stretches and normalized stress components are discussed in detail under the effects of m.The results indicate that the material inhomogeneity and the rotating angular velocity have significant effects on the distributions of the normalized radial and hoop stress components and the stretches.We believe that by appropriately choosing the material inhomogeneity and configuration parameters,the functionally-graded material(FGM)hyper-elastic hollow cylindrical disk can be designed to meet some unique requirements in the application fields,e.g.,soft robotics,medical devices,and conventional aerospace and mechanical industries.
文摘The dynamical formation of cavity in a hyper_elastic sphere composed of two materials with the incompressible strain energy function, subjected to a suddenly applied uniform radial tensile boundary dead_load, was studied following the theory of finite deformation dynamics. Besides a trivial solution corresponding to the homogeneous static state, a cavity forms at the center of the sphere when the tensile load is larger than its critical value. An exact differential relation between the cavity radius and the tensile land was obtained. It is proved that the evolution of cavity radius with time displays nonlinear periodic oscillations. The phase diagram for oscillation, the maximum amplitude, the approximate period and the critical load were all discussed.
基金supported by the National Natural Science Foundation of China (Nos. 10772104 and 10872045)the Innovation Project of Shanghai Municipal Education Commission (No. 09YZ12)the Shanghai Leading Academic Discipline Project (No. S30106)
文摘The growth of a prolate or oblate elliptic micro-void in a fiber reinforced anisotropic incompressible hyper-elastic rectangular thin plate subjected to uniaxial extensions is studied within the framework of finite elasticity. Coupling effects of void shape and void size on the growth of the void are paid special attention to. The deformation function of the plate with an isolated elliptic void is given, which is expressed by two parameters to solve the differential equation. The solution is approximately obtained from the minimum potential energy principle. Deformation curves for the void with a wide range of void aspect ratios and the stress distributions on the surface of the void have been obtained by numerical computation. The growth behavior of the void and the characteristics of stress distributions on the surface of the void are captured. The combined effects of void size and void shape on the growth of the void in the thin plate are discussed. The maximum stresses for the void with different sizes and different void aspect ratios are compared.
基金support by Fondazione Cariplo and Fondazione CDP(Italy)under the project No.2022-1895.
文摘We construct an unconventional divergence preserving discretization of updated Lagrangian ideal magnetohydrodynamics(MHD)over simplicial grids.The cell-centered finite-volume(FV)method employed to discretize the conservation laws of volume,momentum,and total energy is rigorously the same as the one developed to simulate hyperelasticity equations.By construction this moving mesh method ensures the compatibility between the mesh displacement and the approximation of the volume flux by means of the nodal velocity and the attached unit corner normal vector which is nothing but the partial derivative of the cell volume with respect to the node coordinate under consideration.This is precisely the definition of the compatibility with the Geometrical Conservation Law which is the cornerstone of any proper multi-dimensional moving mesh FV discretization.The momentum and the total energy fluxes are approximated utilizing the partition of cell faces into sub-faces and the concept of sub-face force which is the traction force attached to each sub-face impinging at a node.We observe that the time evolution of the magnetic field might be simply expressed in terms of the deformation gradient which characterizes the Lagrange-to-Euler mapping.In this framework,the divergence of the magnetic field is conserved with respect to time thanks to the Piola formula.Therefore,we solve the fully compatible updated Lagrangian discretization of the deformation gradient tensor for updating in a simple manner the cell-centered value of the magnetic field.Finally,the sub-face traction force is expressed in terms of the nodal velocity to ensure a semi-discrete entropy inequality within each cell.The conservation of momentum and total energy is recovered prescribing the balance of all the sub-face forces attached to the sub-faces impinging at a given node.This balance corresponds to a vectorial system satisfied by the nodal velocity.It always admits a unique solution which provides the nodal velocity.The robustness and the accuracy of this unconventional FV scheme have been demonstrated by employing various representative test cases.Finally,it is worth emphasizing that once you have an updated Lagrangian code for solving hyperelasticity you also get an almost free updated Lagrangian code for solving ideal MHD ensuring exactly the compatibility with the involution constraint for the magnetic field at the discrete level.
文摘An explicit integration scheme for rate-dependent crystal plasticity （CP） was developed. Additive decomposition of the velocity gradient tensor into lattice and plastic parts is adopted for describing the kinematics; the Cauchy stress is calculated by using a hypo-elastic formulation, applying the Jaumann stress rate. This CP scheme has been implemented into a commercial finite element code （CPFEM）. Uniaxial compression and roiling processes were simulated. The results show good accuracy and reliability of the integration scheme. The results were compared with simulations using one hyper-elastic CPFEM implementation which involves multiplicative decomposition of the deformation gradient tensor. It is found that the hypo-elastic implementation is only slightly faster and has a similar accuracy as the hyper-elastic formulation.
文摘Various mixed formulations of the finite element method (FEM) yield matrix equations involving zero diagonal entries. They are then dealt with by a penaltymethod so that they become non-zero but near zero terms. However, the penalty has tobe chosen properly. If it is too large, the matrix equation may become ill-conditioned. Onthe other hand, the matrix equation may give incorrect answer if the penalty is too small.In non-linear regime, the difficulty is more serious because the magnitude order of the matrix varies considerably in the entire loading history. The paper suggests an iteration solution and applies it to non-linear FEM of rubber-like hyper-elasticity. This type of analysisis highly non-linear both in physics and in geometry as well as the strong constraint of incompressibility. The iteration solution is demonstrated to possess super precision and excellent convergence characteristics.
基金Supported by National Natural Science Foundation of China(Grant No.51605001)Joint Funds of the National Natural Science Foundation of China(Grant No.U1637207)Anhui University Research Foundation for Doctor(Grant No.J01003222)
文摘Because of the limited space of the launch rockets, deployable mechanisms are always used to solve the phenomenon. One dimensional deployable mast can deploy and support antenna, solar sail and space optical camera. Tape-spring hyperelastic hinges can be folded and extended into a rod like configuration. It utilizes the strain energy to realize self-deploying and drive the other structures. One kind of triangular prism mast with tape-spring hyperelastic hinges is proposed and developed. Stretching and compression stiffness theoretical model are established with considering the tape-spring hyperelastic hinges based on static theory. The finite element model of ten-module triangular prism mast is set up by ABAQUS with the tape-spring hyperelastic hinge and parameter study is performed to investigate the influence of thickness, section angle and radius. Two-module TPM is processed and tested the compression stiffness by the laser displacement sensor, deploying repeat accuracy by the high speed camera, modal shape and fundamental frequency at cantilever position by LMS multi-channel vibration test and analysis system, which are used to verify precision of the theoretical and finite element models of ten-module triangular prism mast with the tape-spring hyperelastic hinges. This research proposes an innovative one dimensional triangular prism with tape-spring hyperelastic hinge which has great application value to the space deployable mechanisms.
基金Project supported by the National Natural Science Foundation of China (Nos. 10772104 and 10872045)the Innovation Project of Shanghai Municipal Education Commission (No. 09YZ12)the Shanghai Leading Academic Discipline Project (No. S30106)
文摘Mechanical properties, such as the deformation and stress distributions for venous walls under the combined load of transmural pressure and axial stretch, are examined within the framework of nonlinear elasticity with one kind of hyper-elastic strain energy functions. The negative pressure instability problem of the venous wall is explained through energy comparison. First, the deformation equation of the venous wall under the combined loads is obtained with a thin-walled circular cylindrical tube. The deformation curves and the stress distributions for the venous wall are given under the normal transmural pressure, and the regulations are discussed. Then, the deformation curves of the venous wall under the negative transmural pressure or the internal pressure less than the external pressure are given. Finally, the negative pressure instability problem is discussed through energy comparison.
文摘In the present work we describe a new expression of total stress taking into account the passive and active contributions but especially the pore level stress.Special attention is paid to the effective stress and osmotic pressure gradient in numerical simulation to understand the mechanical behavior of the human arterial wall.The new model aims to predict the rigidity of the artery,by using the theoretical model of hyper-elastic,anisotropic and dynamical behavior of human common carotid artery.The principal obtained result showed that:the osmosis phenomenon is the best parameter to explain the loss water in arterial tissue.This loss of water causes the rigidity of the artery which thus can be controlled by the osmosis phenomenon.All the results are in good agreement with the expected results of the literature and could play the important role in the diagnosis of the patients with the CVD(Cardiovascular Disease).
基金Anusandhan National Research Foundation(ANRF)via grant no.CRG/2023/007667 and from the Ministry of Port,Shipping,and Waterways via letter no.ST-14011/74/MT(356529).
文摘Constitutive modeling is crucial for engineering design and simulations to accurately describe material behavior.However,traditional phenomenological models often struggle to capture the complexities of real materials under varying stress conditions due to their fixed forms and limited parameters.While recent advances in deep learning have addressed some limitations of classical models,purely data-driven methods tend to require large data sets,lack interpretability,and struggle to generalize beyond their training data.To tackle these issues,we introduce“Fusion-based Constitutive model(FuCe):Toward model-data augmentation in constitutive modeling.”This approach combines established phenomenological models with an Input Convex Neural Network architecture,designed to train on the limited and noisy force-displacement data typically available in practical applications.The hybrid model inherently adheres to necessary constitutive conditions.During inference,Monte Carlo dropout is employed to generate Bayesian predictions,providing mean values and confidence intervals that quantify uncertainty.We demonstrate the model's effectiveness by learning two isotropic constitutive models and one anisotropic model with a single fiber direction,across six different stress states.The framework's applicability is also showcased in finite element simulations across three geometries of varying complexities.Our results highlight the framework's superior extrapolation capabilities,even when trained on limited and noisy data,delivering accurate and physically meaningful predictions across all numerical examples.
基金supported by the National Natural Science Foundation of China (Grant Nos.10772104 and 10872045)the innovation project of Shanghai Municipal Education Commission (Grant No.09YZ12)Shanghai Leading Academic Discipline Project (Grant No.S30106)
文摘Dynamical cavitation and oscillation of an anisotropic two-family fiber-reinforced incompressible hyper-elastic sphere subjected to a suddenly applied constant boundary dead load are examined within the framework of finite elasto-dynamics.An exact differential equation between the radius of the cavity and the applied load is obtained.The curves for the variation of the maximum radius of the cavity with the load and the phase diagrams are obtained by vibration theories and numerical computation.It is shown that there exists a critical value for the applied load.When the applied load is larger than the critical value,a spherical cavity will suddenly form at the center of the sphere.It is proved that the evolution of the cavity radius with time follows that of nonlinear periodic oscillation,and oscillation of the anisotropic sphere is not the same as that of the isotropic sphere.
