This paper presents a nonlinear micropolar nonclassical continuum theory (MPNCCT) for finite deformation, finite strain deformation physics of thermosviscoelastic solid medium with memory (polymeric micropolar solids)...This paper presents a nonlinear micropolar nonclassical continuum theory (MPNCCT) for finite deformation, finite strain deformation physics of thermosviscoelastic solid medium with memory (polymeric micropolar solids) based on classical rotations cΘand their rates. Contravariant second Piola-Kirchhoff stress and moment tensors, in conjunction with finite deformation measures derived by the authors in recent paper, are utilized in deriving the conservation and balance laws and the constitutive theories based on conjugate pairs in entropy inequality and the representation theorem. This nonlinear MPNCCT for TVES with rheology: 1) incorporates nonlinear ordered rate dissipation mechanism based on Green’s strain rates up to order n;2) also incorporates an additional ordered rate dissipation mechanism due to microconstituents, the viscosity of the medium and the rates of the symmetric part of the rotation gradient (of cΘ) tensor up to order n, referred to as micropolar dissipation or micropolar viscous dissipation mechanism;3) incorporates the primary mechanism of memory or rheology due to long chain molecules of the polymer and the viscosity of the medium by using the contravaraint second Piola-Kirchhoff stress tensor and its rates up to order m, resulting in a relaxation spectrum;4) incorporates second mechanism of memory or rheology due to nonclassical physics, interaction of microconstituents with the viscous medium and long chain molecules by considering rates of the contravariant second Piola-Kirchhoff moment tensor up to order m, resulting in relaxation of second Piola-Kirchhoff moment tensor. This results in another relaxation spectrum for the second Piola-Kirchhoff moment tensor due to microconstituents, referred to as micropolar relaxation spectrum consisting of micropolar relaxation time constants of the material. This nonlinear MPNCCT for TVES with memory is thermodynamically and mathematically consistent, and the mathematical model consisting of conservation and balance laws and the constitutive theories has closure and naturally reduces to linear MPNCCT based on infinitesimal deformation assumption. BMM is the essential balance law for all MPNCCT and is used in the present work as well. In the absence of this balance law, a valid thermodynamically and mathematically consistent nonlinear MPNCCT is not possible. The nonlinear MPNCCT based on rotations (cΘ+αΘ) and αΘ(ignoring cΘ) is not considered due to the fact that even the linear MPNCCT based on these rotations is invalid and is thermodynamically and mathematically inconsistent MPNCCT.展开更多
In the present paper, some additional new definitions on the kinematics and dynamics are introduced, and the dynamical equations of Boussinesq type, Kirchhoff type, Signorini type and Nowozilov type for finite deforma...In the present paper, some additional new definitions on the kinematics and dynamics are introduced, and the dynamical equations of Boussinesq type, Kirchhoff type, Signorini type and Nowozilov type for finite deformable polar elastic media are systematically derived from the consideration of Euler angles as angular coordinates and the dynamical equations of Cauchy type presented by Dluzewski.展开更多
This paper proposes a geometrically nonlinear total Lagrangian Galerkin meshfree formulation based on the stabilized conforming nodal integration for efficient analysis of shear deformable beam.The present nonlinear a...This paper proposes a geometrically nonlinear total Lagrangian Galerkin meshfree formulation based on the stabilized conforming nodal integration for efficient analysis of shear deformable beam.The present nonlinear analysis encompasses the fully geometric nonlinearities due to large deflection,large deformation as well as finite rotation.The incremental equilibrium equation is obtained by the consistent linearization of the nonlinear variational equation.The Lagrangian meshfree shape function is utilized to discretize the variational equation.Subsequently to resolve the shear and membrane locking issues and accelerate the computation,the method of stabilized conforming nodal integration is systematically implemented through the Lagrangian gradient smoothing operation.Numerical results reveal that the present formulation is very effective.展开更多
Recently,numerous studies have demonstrated that the physics-informed neural network(PINN)can effectively and accurately resolve hyperelastic finite deformation problems.In this paper,a PINN framework for tackling hyp...Recently,numerous studies have demonstrated that the physics-informed neural network(PINN)can effectively and accurately resolve hyperelastic finite deformation problems.In this paper,a PINN framework for tackling hyperelastic-magnetic coupling problems is proposed.Since the solution space consists of two-phase domains,two separate networks are constructed to independently predict the solution for each phase region.In addition,a conscious point allocation strategy is incorporated to enhance the prediction precision of the PINN in regions characterized by sharp gradients.With the developed framework,the magnetic fields and deformation fields of magnetorheological elastomers(MREs)are solved under the control of hyperelastic-magnetic coupling equations.Illustrative examples are provided and contrasted with the reference results to validate the predictive accuracy of the proposed framework.Moreover,the advantages of the proposed framework in solving hyperelastic-magnetic coupling problems are validated,particularly in handling small data sets,as well as its ability in swiftly and precisely forecasting magnetostrictive motion.展开更多
This paper presents a mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and ordered rate constitutive theories in Lagrangian description derived using entropy ...This paper presents a mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and ordered rate constitutive theories in Lagrangian description derived using entropy inequality and the representation theorem for thermoviscoelastic solids (TVES) with rheology. The CBL and the constitutive theories take into account finite deformation and finite strain deformation physics and are based on contravariant deviatoric second Piola-Kirchhoff stress tensor and its work conjugate covariant Green’s strain tensor and their material derivatives of up to order m and n respectively. All published works on nonlinear dynamics of TVES with rheology are mostly based on phenomenological mathematical models. In rare instances, some aspects of CBL are used but are incorrectly altered to obtain mass, stiffness and damping matrices using space-time decoupled approaches. In the work presented in this paper, we show that this is not possible using CBL of CCM for TVES with rheology. Thus, the mathematical models used currently in the published works are not the correct description of the physics of nonlinear dynamics of TVES with rheology. The mathematical model used in the present work is strictly based on the CBL of CCM and is thermodynamically and mathematically consistent and the space-time coupled finite element methodology used in this work is unconditionally stable and provides solutions with desired accuracy and is ideally suited for nonlinear dynamics of TVES with memory. The work in this paper is the first presentation of a mathematical model strictly based on CBL of CCM and the solution of the mathematical model is obtained using unconditionally stable space-time coupled computational methodology that provides control over the errors in the evolution. Both space-time coupled and space-time decoupled finite element formulations are considered for obtaining solutions of the IVPs described by the mathematical model and are presented in the paper. Factors or the physics influencing dynamic response and dynamic bifurcation for TVES with rheology are identified and are also demonstrated through model problem studies. A simple model problem consisting of a rod (1D) of TVES material with memory fixed at one end and subjected to harmonic excitation at the other end is considered to study nonlinear dynamics of TVES with rheology, frequency response as well as dynamic bifurcation phenomenon.展开更多
The inflation tests of rubbery membranes have been widely employed as an efficient method to characterize the stress response as biaxial loading states.However,most of the previous theoretical works have employed clas...The inflation tests of rubbery membranes have been widely employed as an efficient method to characterize the stress response as biaxial loading states.However,most of the previous theoretical works have employed classic hyperelastic models to analyze the deformation behaviors of inflated membranes.The classic models have been demonstrated to lack the ability to capturing the biaxial deformation of rubbers.To address this issue,we have combined the analytical method and the finite element simulation to investigate the deformation response of soft membranes with different constitutive relationships.For the analytical method,the governing ordinary differential equations have been set up for the boundary value problem of inflation tests and further solved using the shooting method.The analytical results are consistent with those obtained from finite element simulation.The results show that the deformation belongs to the unequal biaxial condition rather than the equi-biaxial state unless a neo-Hookean model is adopted.We also perform a parameter study using the extended eight-chain model,which shows that a change in different parameters affects the mechanical response of inflation tests variously.This work may shed light on the future experimental characterization of soft materials using inflation experiments.展开更多
The theory of plasticity as a special field of continuum mechanics deals with the irreversible,i.e.permanent,deformation of solids.Under the action of given loads or deformations,the state of the stresses and strains ...The theory of plasticity as a special field of continuum mechanics deals with the irreversible,i.e.permanent,deformation of solids.Under the action of given loads or deformations,the state of the stresses and strains or the strain rates in these bodies is described.In this way,it complements the theory of elasticity for the reversible behavior of solids.In practice,it has been observed that many materials behave elastically up to a certain load(yield point),beyond that load,however,increasingly plastic or liquid-like.The combination of these two material properties is known as elastoplasticity.The classical elastoplastic material behavior is assumed to be time-independent or rate-independent.In contrast,we call a time-or rate-dependent behavior visco-elastoplastic and visco-plastic—if the elastic part of the deformation is neglected.In plasticity theory,because of the given loads the states of the state variables stress,strain and temperature as well as their changes are described.For this purpose,the observed phenomena are introduced and put into mathematical relationships.The constitutive relations describing the specific material behavior are finally embedded in the fundamental relations of continuum theory and physics.Historically,the theory of plasticity was introduced in order to better estimate the strength of constructions.An analysis based purely on elastic codes is not in a position to do this,and can occasionally even lead to incorrect interpretations.On the other hand,the entire field of forming techniques requires a theory for the description of plastic behavior.Starting from the classical description of plastic behavior with small deformations,the present review is intended to provide an insight into the state of the art when taking into account finite deformations.展开更多
The overbroken rock mass of gob areas is made up of broken and accumulated rock blocks compressed to some extent by the overlying strata. The beating pressure of the gob can directly affect the safety of mining fields...The overbroken rock mass of gob areas is made up of broken and accumulated rock blocks compressed to some extent by the overlying strata. The beating pressure of the gob can directly affect the safety of mining fields, formarion of road retained along the next goaf and seepage of water and methane through the gob. In this paper, the software RFPA＇2000 is used to construct numerical models. Especially the Euler method of control volume is proposed to solve the simulation difficulty arising from plastically finite deformations. The results show that three characteristic regions occurred in the gob area： （1） a naturally accumulated region, 0-10 m away from unbroken surrounding rock walls, where the beating pressure is nearly zero; （2） an overcompacted region, 10-20 m away from unbroken walls, where the beating pressure results in the maximum value of the gob area; （3） a stable compaction region, more than 20 m away from unbroken walls and occupying absolutely most of the gob area, where the beating pressures show basically no differences. Such a characteristic can exolain the easy-seeoaged “O”-ring phenomena around mining fields very well.展开更多
A nonlinear wave equation of elastic rod taking account of finite deformation, transverse inertia and shearing strain is derived by means of the Hamilton principle in this paper. Nonlinear wave equation and truncated ...A nonlinear wave equation of elastic rod taking account of finite deformation, transverse inertia and shearing strain is derived by means of the Hamilton principle in this paper. Nonlinear wave equation and truncated nonlinear wave equation are solved by the Jacobi elliptic sine function expansion and the third kind of Jacobi elliptic function expansion method. The exact periodic solutions of these nonlinear equations are obtained, including the shock wave solution and the solitary wave solution. The necessary condition of exact periodic solutions, shock solution and solitary solution existence is discussed.展开更多
On the basis of classical linear theory on longitudinal, torsional and flexural waves in thin elastic rods, and taking finite deformation and dispersive effects into consideration, three kinds of nonlinear evolution e...On the basis of classical linear theory on longitudinal, torsional and flexural waves in thin elastic rods, and taking finite deformation and dispersive effects into consideration, three kinds of nonlinear evolution equations are derived. Qualitative analysis of three kinds of nonlinear equations are presented. It is shown that these equations have homoclinic or heteroclinic orbits on the phase plane, corresponding to solitary wave or shock wave solutions, respectively. Based on the principle of homogeneous balance, these equations are solved with the Jacobi elliptic function expansion method. Results show that existence of solitary wave solution and shock wave solution is possible under certain conditions. These conclusions are consistent with qualitative analysis.展开更多
It is noted that the behavior of most piezoelectric materials is temperaturedependent and such piezo-thermo-elastic coupling phenomenon has become even more pronounced in thecase of finite deformation. On the other ha...It is noted that the behavior of most piezoelectric materials is temperaturedependent and such piezo-thermo-elastic coupling phenomenon has become even more pronounced in thecase of finite deformation. On the other hand, for the purpose of precise shape and vibrationcontrol of piezoelectric smart structures, their deformation under external excitation must beideally modeled. This demands a thorough study of the coupled piezo-thermo-elastic response underfinite deformation. In this study, the governing equations of piezoelectric structures areformulated through the theory of virtual displacement principle and a finite element method isdeveloped. It should be emphasized that in the finite element method the fully coupledpiezo-thermo-elastic behavior and the geometric non-linearity are considered. The method developedis then applied to simulate the dynamic and steady response of a clamped plate to heat flux actingon one side of the plate to mimic the behavior of a battery plate of satellite irradiated under thesun. The results obtained are compared against classical solutions, whereby the thermal conductivityis assumed to be independent of deformation. It is found that the full-coupled theory predicts lesstransient response of the temperature compared to the classic analysis. In the steady state limit,the predicted temperature distribution within the plate for small heat flux is almost the same forboth analyses. However, it is noted that increasing the heat flux will increase the deviationbetween the predictions of the temperature distribution by the full coupled theory and by theclassic analysis. It is concluded from the present study that, in order to precisely predict thedeformation of smart structures, the piezo-thermo-elastic coupling, geometric non-linearity and thedeformation dependent thermal conductivity should be taken into account.展开更多
A series of monotonic tensile experiments of thermo-induced shape memory polyurethane （TSMPU） at different loading rates were carried out to investigate the interaction between the internal heat production and the m...A series of monotonic tensile experiments of thermo-induced shape memory polyurethane （TSMPU） at different loading rates were carried out to investigate the interaction between the internal heat production and the mechanical deformation. It is shown that the tem- perature variation on the surfaces of the specimens due to the internal heat production affects the mechanical properties of TSMPU remarkably. Then, based on irreversible thermodynamics, the Helmholtz free energy was decomposed into three parts, i.e., the instantaneous elastic free energy, visco-plastic free energy and heat free energy. The total deformation gradient was decomposed into the mechanical and thermal parts, and the mechanical deformation gradient was further divided into the elastic and visco-plastic components. The Hencky＇s logarithmic strain was used in the current configuration. The heat equilibrium equation of internal heat production and heat exchange was derived in accordance with the first and second thermodynamics laws. The temperature of specimens was contributed by the internal heat production and the ambient temperature simultaneously, and a thermo-mechanically coupled thermo-elasto-visco-plastie model was established. The effect of temperature variation of specimens on the mechanical properties of the material was considered in this work. Finally, the capability of the proposed model was validated by comparing the simulated results with the corresponding experimental data of TSMPU.展开更多
The problem of finite deformation of an incompressible rectangular rubber ring with an internal rigid body, where the ring is subjected to equal axial loads at its two ends, is examined. A reasonable mathematical mode...The problem of finite deformation of an incompressible rectangular rubber ring with an internal rigid body, where the ring is subjected to equal axial loads at its two ends, is examined. A reasonable mathematical model is formulated by using the nonlinear field theory and the implicit analytical solutions are derived. Then numerical simulations are implemented to further illustrate the results and obtain some meaningful conclusions. The deformation of the lateral surface of the ring becomes larger with the increasing axial loads, the decreasing ratio of the inner and outer radii and the increasing height of the ring.展开更多
A new finite strain elatoplastic J2-flow model with coupling effects of both isotropic and anisotropic hardening is proposed with the co-rotational logarithmic rate.In terms of certain single-variable shape functions ...A new finite strain elatoplastic J2-flow model with coupling effects of both isotropic and anisotropic hardening is proposed with the co-rotational logarithmic rate.In terms of certain single-variable shape functions representing uniaxial loading and unloading curves,explicit multi-axial expressions for the three hardening quantities incorporated in the new model proposed are derived in unified forms for the purpose of automatically and accurately simulating complex pseudoelastic-to-plastic transition effects of shape memory alloys(SMAs)under multiple loading-unloading cycles.Numerical examples show that with only a single parameter of direct physical meaning for each cycle,accurate and explicit simulations may be achieved for extensive data from multiple cycle tests.展开更多
In this paper, a finite element method （FEM）-based multi-phase problem based on a newly proposed thermal elastoplastic constitutive model for saturated/unsaturated geomaterial is discussed. A program of FEM named as...In this paper, a finite element method （FEM）-based multi-phase problem based on a newly proposed thermal elastoplastic constitutive model for saturated/unsaturated geomaterial is discussed. A program of FEM named as SOFT, adopting unified field equations for thermo-hydro-mechanical-air （THMA） behavior of geomaterial and using finite element-finite difference （FE-FD） scheme for so/l-water-air three-phase coupling problem, is used in the numerical simulation. As an application of the newly proposed numerical method, two engineering problems, one for slope failure in unsaturated model ground and another for in situ heating test related to deep geological repository of high-level radioactive waste （HLRW）, are simulated. The model tests on slope failure in unsaturated Shirasu ground, carried out by Kitamura et al. （2007）, is simulated in the framework of soil-water-air three-phase coupling under the condition of constant temperature. While the in situ heating test reported by Munoz （2006） is simulated in the same framework under the conditions of variable temperature hut constant air pressure.展开更多
In engineering component design,material models are increasingly used in finite element simulations for an expeditious and less costly analysis of the design prototypes.As such,researchers strive to formulate models t...In engineering component design,material models are increasingly used in finite element simulations for an expeditious and less costly analysis of the design prototypes.As such,researchers strive to formulate models that are less complex,robust,and accurate.In the realm of hyperelastic materials,phenomenological-based Carroll’s model is highly promising due to its simplicity and accuracy.This work proposes its further improvement by modifying the strain energy density function to comply with the restriction that it should vanish at reference configuration and adding a compressible term to capture the practical behavior of elastomeric materials and to avoid numerical problems during finite element implementation.The model constants for both the original and the modified versions were obtained by fitting their respective expressions to the classical Treloar’s experimental data using the Levenberg-Marquardt algorithm.The modified model was implemented using Abaqus CAE 2016 via a vectorized user material(VUMAT)subroutine.Comparisons of the model predictions with Treloar’s experimental data demonstrated the superiority of the modified version particularly in the equibiaxial loading mode.Moreover,the simulated and the experimentally observed behavior agreed to a great accuracy thus making the modified model suitable for simulating the loading response of components fabricated of elastomeric materials.展开更多
By using Hamilton-type variation principle in non-conservation system, the nonlinear equation of wave motion of a elastic thin rod was derived according to Lagrange description of finite deformation theory. The dissip...By using Hamilton-type variation principle in non-conservation system, the nonlinear equation of wave motion of a elastic thin rod was derived according to Lagrange description of finite deformation theory. The dissipation caused due to viscous effect and the dispersion introduced by transverse inertia were taken into consideration so that steady traveling wave solution can be obtained. Using multi-scale method the nonlinear equation is reduced to a KdV-Burgers equation which corresponds with saddle-spiral heteroclinic orbit on phase plane. Its solution is called the oscillating-solitary wave or saddle-spiral shock wave. If viscous effect or transverse inertia is neglected, the equation is degraded to classical KdV or Burgers equation. The former implies a propagating solitary wave with homoclinic on phase plane, the latter means shock wave and heteroclinic orbit.展开更多
Based on the Reddy's theory of plates with the effect of higher-order shear deformations, the governing equations for bending of orthotropic plates with finite deformations were established. The differential quadr...Based on the Reddy's theory of plates with the effect of higher-order shear deformations, the governing equations for bending of orthotropic plates with finite deformations were established. The differential quadrature (DQ) method of nonlinear analysis to the problem was presented. New DQ approach, presented by Wang and Bert (DQWB), is extended to handle the multiple boundary conditions of plates. The techniques were also further extended to simplify nonlinear computations. The numerical convergence and comparison of solutions were studied. The results show that the DQ method presented is very reliable and valid. Moreover, the influences of geometric and material parameters as well as the transverse shear deformations on nonlinear bending were investigated. Numerical results show the influence of the shear deformation on the static bending of orthotropic moderately thick plate is significant.展开更多
In this paper, based on the finite deformation S-R decomposition theorem, the definition of the body moment is renewed as the stem of its internal and external. The expression of the increment rate of the deformation ...In this paper, based on the finite deformation S-R decomposition theorem, the definition of the body moment is renewed as the stem of its internal and external. The expression of the increment rate of the deformation energy is derived and the physical meaning is clarified. The power variational principle and the complementary power variational principle for finite deformation mechanics are supplemented and perfected.展开更多
A systematic study has been conducted aiming to attain an insight into the influence of coefficient of roll speed asymmetry, crystal orientation and structure on the deformation behavior, and crystallographic orientat...A systematic study has been conducted aiming to attain an insight into the influence of coefficient of roll speed asymmetry, crystal orientation and structure on the deformation behavior, and crystallographic orientation development during foil rolling. Simulations were successfully carried out by using crystal plasticity finite element method（CPFEM）,and a novel computational framework is presented for the representation of virtual polycrystalline grain structures. It has been found that asymmetric rolling（ASR） is more efficient in producing plastic deformation since it develops additional shear strain and activity of slip system compared with symmetric rolling（SR）. For ASR, increase in the length of the shear zone, and decrease in the amount of the pressure and roll force would lead to further reduction. The shear strain path in SR and ASR is strictly influenced by the misorientation of neighbor grains, and corresponding {1 1 1} pole figures offer direct evidence of the spread of crystallographic orientation around the normal direction. The activity of slip systems was examined in detail and found that the predicted results are consistent with the surface layer model. The accuracy of the developed CPFEM model is verified by the fact that the simulated results of roll force coincide well with the experimental results.展开更多
文摘This paper presents a nonlinear micropolar nonclassical continuum theory (MPNCCT) for finite deformation, finite strain deformation physics of thermosviscoelastic solid medium with memory (polymeric micropolar solids) based on classical rotations cΘand their rates. Contravariant second Piola-Kirchhoff stress and moment tensors, in conjunction with finite deformation measures derived by the authors in recent paper, are utilized in deriving the conservation and balance laws and the constitutive theories based on conjugate pairs in entropy inequality and the representation theorem. This nonlinear MPNCCT for TVES with rheology: 1) incorporates nonlinear ordered rate dissipation mechanism based on Green’s strain rates up to order n;2) also incorporates an additional ordered rate dissipation mechanism due to microconstituents, the viscosity of the medium and the rates of the symmetric part of the rotation gradient (of cΘ) tensor up to order n, referred to as micropolar dissipation or micropolar viscous dissipation mechanism;3) incorporates the primary mechanism of memory or rheology due to long chain molecules of the polymer and the viscosity of the medium by using the contravaraint second Piola-Kirchhoff stress tensor and its rates up to order m, resulting in a relaxation spectrum;4) incorporates second mechanism of memory or rheology due to nonclassical physics, interaction of microconstituents with the viscous medium and long chain molecules by considering rates of the contravariant second Piola-Kirchhoff moment tensor up to order m, resulting in relaxation of second Piola-Kirchhoff moment tensor. This results in another relaxation spectrum for the second Piola-Kirchhoff moment tensor due to microconstituents, referred to as micropolar relaxation spectrum consisting of micropolar relaxation time constants of the material. This nonlinear MPNCCT for TVES with memory is thermodynamically and mathematically consistent, and the mathematical model consisting of conservation and balance laws and the constitutive theories has closure and naturally reduces to linear MPNCCT based on infinitesimal deformation assumption. BMM is the essential balance law for all MPNCCT and is used in the present work as well. In the absence of this balance law, a valid thermodynamically and mathematically consistent nonlinear MPNCCT is not possible. The nonlinear MPNCCT based on rotations (cΘ+αΘ) and αΘ(ignoring cΘ) is not considered due to the fact that even the linear MPNCCT based on these rotations is invalid and is thermodynamically and mathematically inconsistent MPNCCT.
文摘In the present paper, some additional new definitions on the kinematics and dynamics are introduced, and the dynamical equations of Boussinesq type, Kirchhoff type, Signorini type and Nowozilov type for finite deformable polar elastic media are systematically derived from the consideration of Euler angles as angular coordinates and the dynamical equations of Cauchy type presented by Dluzewski.
基金supported by the National Natural Science Foundation of China (10972188)the Program for New Century Excellent Talents in University from China Education Ministry (NCET-09-0678)
文摘This paper proposes a geometrically nonlinear total Lagrangian Galerkin meshfree formulation based on the stabilized conforming nodal integration for efficient analysis of shear deformable beam.The present nonlinear analysis encompasses the fully geometric nonlinearities due to large deflection,large deformation as well as finite rotation.The incremental equilibrium equation is obtained by the consistent linearization of the nonlinear variational equation.The Lagrangian meshfree shape function is utilized to discretize the variational equation.Subsequently to resolve the shear and membrane locking issues and accelerate the computation,the method of stabilized conforming nodal integration is systematically implemented through the Lagrangian gradient smoothing operation.Numerical results reveal that the present formulation is very effective.
基金supported by the National Natural Science Foundation of China(Nos.12072105 and 11932006)。
文摘Recently,numerous studies have demonstrated that the physics-informed neural network(PINN)can effectively and accurately resolve hyperelastic finite deformation problems.In this paper,a PINN framework for tackling hyperelastic-magnetic coupling problems is proposed.Since the solution space consists of two-phase domains,two separate networks are constructed to independently predict the solution for each phase region.In addition,a conscious point allocation strategy is incorporated to enhance the prediction precision of the PINN in regions characterized by sharp gradients.With the developed framework,the magnetic fields and deformation fields of magnetorheological elastomers(MREs)are solved under the control of hyperelastic-magnetic coupling equations.Illustrative examples are provided and contrasted with the reference results to validate the predictive accuracy of the proposed framework.Moreover,the advantages of the proposed framework in solving hyperelastic-magnetic coupling problems are validated,particularly in handling small data sets,as well as its ability in swiftly and precisely forecasting magnetostrictive motion.
文摘This paper presents a mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and ordered rate constitutive theories in Lagrangian description derived using entropy inequality and the representation theorem for thermoviscoelastic solids (TVES) with rheology. The CBL and the constitutive theories take into account finite deformation and finite strain deformation physics and are based on contravariant deviatoric second Piola-Kirchhoff stress tensor and its work conjugate covariant Green’s strain tensor and their material derivatives of up to order m and n respectively. All published works on nonlinear dynamics of TVES with rheology are mostly based on phenomenological mathematical models. In rare instances, some aspects of CBL are used but are incorrectly altered to obtain mass, stiffness and damping matrices using space-time decoupled approaches. In the work presented in this paper, we show that this is not possible using CBL of CCM for TVES with rheology. Thus, the mathematical models used currently in the published works are not the correct description of the physics of nonlinear dynamics of TVES with rheology. The mathematical model used in the present work is strictly based on the CBL of CCM and is thermodynamically and mathematically consistent and the space-time coupled finite element methodology used in this work is unconditionally stable and provides solutions with desired accuracy and is ideally suited for nonlinear dynamics of TVES with memory. The work in this paper is the first presentation of a mathematical model strictly based on CBL of CCM and the solution of the mathematical model is obtained using unconditionally stable space-time coupled computational methodology that provides control over the errors in the evolution. Both space-time coupled and space-time decoupled finite element formulations are considered for obtaining solutions of the IVPs described by the mathematical model and are presented in the paper. Factors or the physics influencing dynamic response and dynamic bifurcation for TVES with rheology are identified and are also demonstrated through model problem studies. A simple model problem consisting of a rod (1D) of TVES material with memory fixed at one end and subjected to harmonic excitation at the other end is considered to study nonlinear dynamics of TVES with rheology, frequency response as well as dynamic bifurcation phenomenon.
基金supported by the National Natural Science Foundation of China(Grant Nos.12211530061 and 12321002)the Zhejiang Provincial Natural Science Foundation of China(Grant No.LD22A020001)the 111 Project(Grant No.B21034).
文摘The inflation tests of rubbery membranes have been widely employed as an efficient method to characterize the stress response as biaxial loading states.However,most of the previous theoretical works have employed classic hyperelastic models to analyze the deformation behaviors of inflated membranes.The classic models have been demonstrated to lack the ability to capturing the biaxial deformation of rubbers.To address this issue,we have combined the analytical method and the finite element simulation to investigate the deformation response of soft membranes with different constitutive relationships.For the analytical method,the governing ordinary differential equations have been set up for the boundary value problem of inflation tests and further solved using the shooting method.The analytical results are consistent with those obtained from finite element simulation.The results show that the deformation belongs to the unequal biaxial condition rather than the equi-biaxial state unless a neo-Hookean model is adopted.We also perform a parameter study using the extended eight-chain model,which shows that a change in different parameters affects the mechanical response of inflation tests variously.This work may shed light on the future experimental characterization of soft materials using inflation experiments.
文摘The theory of plasticity as a special field of continuum mechanics deals with the irreversible,i.e.permanent,deformation of solids.Under the action of given loads or deformations,the state of the stresses and strains or the strain rates in these bodies is described.In this way,it complements the theory of elasticity for the reversible behavior of solids.In practice,it has been observed that many materials behave elastically up to a certain load(yield point),beyond that load,however,increasingly plastic or liquid-like.The combination of these two material properties is known as elastoplasticity.The classical elastoplastic material behavior is assumed to be time-independent or rate-independent.In contrast,we call a time-or rate-dependent behavior visco-elastoplastic and visco-plastic—if the elastic part of the deformation is neglected.In plasticity theory,because of the given loads the states of the state variables stress,strain and temperature as well as their changes are described.For this purpose,the observed phenomena are introduced and put into mathematical relationships.The constitutive relations describing the specific material behavior are finally embedded in the fundamental relations of continuum theory and physics.Historically,the theory of plasticity was introduced in order to better estimate the strength of constructions.An analysis based purely on elastic codes is not in a position to do this,and can occasionally even lead to incorrect interpretations.On the other hand,the entire field of forming techniques requires a theory for the description of plastic behavior.Starting from the classical description of plastic behavior with small deformations,the present review is intended to provide an insight into the state of the art when taking into account finite deformations.
基金Projects 2005CB221502 supported by the Vital Foundational 973 Program of China, 50225414 by the National Outstanding Youth Foundation,20040350222 by China Postdoctoral Science FoundationBK 2004033 by Jiangsu Natural Science Foundation
文摘The overbroken rock mass of gob areas is made up of broken and accumulated rock blocks compressed to some extent by the overlying strata. The beating pressure of the gob can directly affect the safety of mining fields, formarion of road retained along the next goaf and seepage of water and methane through the gob. In this paper, the software RFPA＇2000 is used to construct numerical models. Especially the Euler method of control volume is proposed to solve the simulation difficulty arising from plastically finite deformations. The results show that three characteristic regions occurred in the gob area： （1） a naturally accumulated region, 0-10 m away from unbroken surrounding rock walls, where the beating pressure is nearly zero; （2） an overcompacted region, 10-20 m away from unbroken walls, where the beating pressure results in the maximum value of the gob area; （3） a stable compaction region, more than 20 m away from unbroken walls and occupying absolutely most of the gob area, where the beating pressures show basically no differences. Such a characteristic can exolain the easy-seeoaged “O”-ring phenomena around mining fields very well.
基金Project supported by the National Natural Science Foundation of China (No. 10472076).
文摘A nonlinear wave equation of elastic rod taking account of finite deformation, transverse inertia and shearing strain is derived by means of the Hamilton principle in this paper. Nonlinear wave equation and truncated nonlinear wave equation are solved by the Jacobi elliptic sine function expansion and the third kind of Jacobi elliptic function expansion method. The exact periodic solutions of these nonlinear equations are obtained, including the shock wave solution and the solitary wave solution. The necessary condition of exact periodic solutions, shock solution and solitary solution existence is discussed.
基金Project supported by the National Natural Science Foundation of China (No.10772129)the Youth Science Foundation of Shanxi Province of China (No.2006021005)
文摘On the basis of classical linear theory on longitudinal, torsional and flexural waves in thin elastic rods, and taking finite deformation and dispersive effects into consideration, three kinds of nonlinear evolution equations are derived. Qualitative analysis of three kinds of nonlinear equations are presented. It is shown that these equations have homoclinic or heteroclinic orbits on the phase plane, corresponding to solitary wave or shock wave solutions, respectively. Based on the principle of homogeneous balance, these equations are solved with the Jacobi elliptic function expansion method. Results show that existence of solitary wave solution and shock wave solution is possible under certain conditions. These conclusions are consistent with qualitative analysis.
基金the National Natural Science Foundation of China (Nos.10132010 and 50135030)the Foundation of In-service Doctors of Xi'an Jiaotong University
文摘It is noted that the behavior of most piezoelectric materials is temperaturedependent and such piezo-thermo-elastic coupling phenomenon has become even more pronounced in thecase of finite deformation. On the other hand, for the purpose of precise shape and vibrationcontrol of piezoelectric smart structures, their deformation under external excitation must beideally modeled. This demands a thorough study of the coupled piezo-thermo-elastic response underfinite deformation. In this study, the governing equations of piezoelectric structures areformulated through the theory of virtual displacement principle and a finite element method isdeveloped. It should be emphasized that in the finite element method the fully coupledpiezo-thermo-elastic behavior and the geometric non-linearity are considered. The method developedis then applied to simulate the dynamic and steady response of a clamped plate to heat flux actingon one side of the plate to mimic the behavior of a battery plate of satellite irradiated under thesun. The results obtained are compared against classical solutions, whereby the thermal conductivityis assumed to be independent of deformation. It is found that the full-coupled theory predicts lesstransient response of the temperature compared to the classic analysis. In the steady state limit,the predicted temperature distribution within the plate for small heat flux is almost the same forboth analyses. However, it is noted that increasing the heat flux will increase the deviationbetween the predictions of the temperature distribution by the full coupled theory and by theclassic analysis. It is concluded from the present study that, in order to precisely predict thedeformation of smart structures, the piezo-thermo-elastic coupling, geometric non-linearity and thedeformation dependent thermal conductivity should be taken into account.
基金Financial supports by National Natural Science Foundation of China (11572265,11202171)Excellent Youth Found of Sichuan Province (2017JQ0019)+1 种基金Open Project of Traction Power State Key Laboratory(TPL1606)Exploration Project of Traction Power State Key Laboratory (2017TPL_T04)
文摘A series of monotonic tensile experiments of thermo-induced shape memory polyurethane （TSMPU） at different loading rates were carried out to investigate the interaction between the internal heat production and the mechanical deformation. It is shown that the tem- perature variation on the surfaces of the specimens due to the internal heat production affects the mechanical properties of TSMPU remarkably. Then, based on irreversible thermodynamics, the Helmholtz free energy was decomposed into three parts, i.e., the instantaneous elastic free energy, visco-plastic free energy and heat free energy. The total deformation gradient was decomposed into the mechanical and thermal parts, and the mechanical deformation gradient was further divided into the elastic and visco-plastic components. The Hencky＇s logarithmic strain was used in the current configuration. The heat equilibrium equation of internal heat production and heat exchange was derived in accordance with the first and second thermodynamics laws. The temperature of specimens was contributed by the internal heat production and the ambient temperature simultaneously, and a thermo-mechanically coupled thermo-elasto-visco-plastie model was established. The effect of temperature variation of specimens on the mechanical properties of the material was considered in this work. Finally, the capability of the proposed model was validated by comparing the simulated results with the corresponding experimental data of TSMPU.
基金supported by the National Natural Science Foundation of China (Nos. 10872045, 10721062 and 10772104)the Program for New Century Excellent Talents in University (No. NCET-09-0096)the Fundamental Research Funds for the Central Universities
文摘The problem of finite deformation of an incompressible rectangular rubber ring with an internal rigid body, where the ring is subjected to equal axial loads at its two ends, is examined. A reasonable mathematical model is formulated by using the nonlinear field theory and the implicit analytical solutions are derived. Then numerical simulations are implemented to further illustrate the results and obtain some meaningful conclusions. The deformation of the lateral surface of the ring becomes larger with the increasing axial loads, the decreasing ratio of the inner and outer radii and the increasing height of the ring.
基金Project supported by the National Natural Science Foundation of China(No.11372172)and the Start-up Fund from Jinan University in Guangzhou of China。
文摘A new finite strain elatoplastic J2-flow model with coupling effects of both isotropic and anisotropic hardening is proposed with the co-rotational logarithmic rate.In terms of certain single-variable shape functions representing uniaxial loading and unloading curves,explicit multi-axial expressions for the three hardening quantities incorporated in the new model proposed are derived in unified forms for the purpose of automatically and accurately simulating complex pseudoelastic-to-plastic transition effects of shape memory alloys(SMAs)under multiple loading-unloading cycles.Numerical examples show that with only a single parameter of direct physical meaning for each cycle,accurate and explicit simulations may be achieved for extensive data from multiple cycle tests.
文摘In this paper, a finite element method （FEM）-based multi-phase problem based on a newly proposed thermal elastoplastic constitutive model for saturated/unsaturated geomaterial is discussed. A program of FEM named as SOFT, adopting unified field equations for thermo-hydro-mechanical-air （THMA） behavior of geomaterial and using finite element-finite difference （FE-FD） scheme for so/l-water-air three-phase coupling problem, is used in the numerical simulation. As an application of the newly proposed numerical method, two engineering problems, one for slope failure in unsaturated model ground and another for in situ heating test related to deep geological repository of high-level radioactive waste （HLRW）, are simulated. The model tests on slope failure in unsaturated Shirasu ground, carried out by Kitamura et al. （2007）, is simulated in the framework of soil-water-air three-phase coupling under the condition of constant temperature. While the in situ heating test reported by Munoz （2006） is simulated in the same framework under the conditions of variable temperature hut constant air pressure.
基金supported by the National Natural Science Foundation of China(Grants 11632005 and 11672086).
文摘In engineering component design,material models are increasingly used in finite element simulations for an expeditious and less costly analysis of the design prototypes.As such,researchers strive to formulate models that are less complex,robust,and accurate.In the realm of hyperelastic materials,phenomenological-based Carroll’s model is highly promising due to its simplicity and accuracy.This work proposes its further improvement by modifying the strain energy density function to comply with the restriction that it should vanish at reference configuration and adding a compressible term to capture the practical behavior of elastomeric materials and to avoid numerical problems during finite element implementation.The model constants for both the original and the modified versions were obtained by fitting their respective expressions to the classical Treloar’s experimental data using the Levenberg-Marquardt algorithm.The modified model was implemented using Abaqus CAE 2016 via a vectorized user material(VUMAT)subroutine.Comparisons of the model predictions with Treloar’s experimental data demonstrated the superiority of the modified version particularly in the equibiaxial loading mode.Moreover,the simulated and the experimentally observed behavior agreed to a great accuracy thus making the modified model suitable for simulating the loading response of components fabricated of elastomeric materials.
文摘By using Hamilton-type variation principle in non-conservation system, the nonlinear equation of wave motion of a elastic thin rod was derived according to Lagrange description of finite deformation theory. The dissipation caused due to viscous effect and the dispersion introduced by transverse inertia were taken into consideration so that steady traveling wave solution can be obtained. Using multi-scale method the nonlinear equation is reduced to a KdV-Burgers equation which corresponds with saddle-spiral heteroclinic orbit on phase plane. Its solution is called the oscillating-solitary wave or saddle-spiral shock wave. If viscous effect or transverse inertia is neglected, the equation is degraded to classical KdV or Burgers equation. The former implies a propagating solitary wave with homoclinic on phase plane, the latter means shock wave and heteroclinic orbit.
文摘Based on the Reddy's theory of plates with the effect of higher-order shear deformations, the governing equations for bending of orthotropic plates with finite deformations were established. The differential quadrature (DQ) method of nonlinear analysis to the problem was presented. New DQ approach, presented by Wang and Bert (DQWB), is extended to handle the multiple boundary conditions of plates. The techniques were also further extended to simplify nonlinear computations. The numerical convergence and comparison of solutions were studied. The results show that the DQ method presented is very reliable and valid. Moreover, the influences of geometric and material parameters as well as the transverse shear deformations on nonlinear bending were investigated. Numerical results show the influence of the shear deformation on the static bending of orthotropic moderately thick plate is significant.
文摘In this paper, based on the finite deformation S-R decomposition theorem, the definition of the body moment is renewed as the stem of its internal and external. The expression of the increment rate of the deformation energy is derived and the physical meaning is clarified. The power variational principle and the complementary power variational principle for finite deformation mechanics are supplemented and perfected.
基金financially supported by the National Natural Science Foundation of China (Nos. 51374069 and U1460107)
文摘A systematic study has been conducted aiming to attain an insight into the influence of coefficient of roll speed asymmetry, crystal orientation and structure on the deformation behavior, and crystallographic orientation development during foil rolling. Simulations were successfully carried out by using crystal plasticity finite element method（CPFEM）,and a novel computational framework is presented for the representation of virtual polycrystalline grain structures. It has been found that asymmetric rolling（ASR） is more efficient in producing plastic deformation since it develops additional shear strain and activity of slip system compared with symmetric rolling（SR）. For ASR, increase in the length of the shear zone, and decrease in the amount of the pressure and roll force would lead to further reduction. The shear strain path in SR and ASR is strictly influenced by the misorientation of neighbor grains, and corresponding {1 1 1} pole figures offer direct evidence of the spread of crystallographic orientation around the normal direction. The activity of slip systems was examined in detail and found that the predicted results are consistent with the surface layer model. The accuracy of the developed CPFEM model is verified by the fact that the simulated results of roll force coincide well with the experimental results.
