The adaptive element techniques of contact problem are studied by means of penalty method, and the error estimators are discussed. Based on error estimators, algorithm of the adaptive element techniques is developed, ...The adaptive element techniques of contact problem are studied by means of penalty method, and the error estimators are discussed. Based on error estimators, algorithm of the adaptive element techniques is developed, then the Gauss - Newton iterations are used which allow the nonlinear problem to be transformed into a sequence of linear sub- problems then easily solved. In addition, the algorithm can be applied into the simulation of de -bonding of fiber - reinforced composites.展开更多
Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed...Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed finite ele- ment equations,thus resulting in a reduction by half in the dimension of final governing equations.Second,an algorithm composed of contact condition probes and elasto-plastic iterations is utilized to solve the governing equation,which distinguishes two kinds of nonlinearities,and makes the solution unique.In addition,Positive-Negative Sequence Modifica- tion Method is used to condense the finite element equations of each substructure and an analytical integration is intro- duced to determine the elasto-plastic status after each time step or each iteration,hence the computational efficiency is en- hanced to a great extent.Finally,several test and practical examples are pressented showing the validity and versatility of these methods and algorithms.展开更多
Based on the analysis of [7] and [10], we present the mixed finite element approximation of the variational inequality resulting from the contact problem in elasticity. The convergence rate of the stress and displacem...Based on the analysis of [7] and [10], we present the mixed finite element approximation of the variational inequality resulting from the contact problem in elasticity. The convergence rate of the stress and displacement field are both improved from O（h3/4） to quasi-optimal O（h│logh│^1/4）. If stronger but reasonable regularity is available, the convergence rate can be optimal O（h）.展开更多
In this paper, we provide a new mixed finite element approximation of the variational inequality resulting from the unilateral contact problem in elasticity. We use the continuous piecewise P2-P1 finite element to app...In this paper, we provide a new mixed finite element approximation of the variational inequality resulting from the unilateral contact problem in elasticity. We use the continuous piecewise P2-P1 finite element to approximate the displacement field and the normal stress component on the contact region. Optimal convergence rates are obtained under the reasonable regularity hypotheses. Numerical example verifies our results.展开更多
The possibilities of the particle finite element method(PFEM)for modeling geotechnical problems are increasingly evident.PFEM is a numerical approach to solve large displacement and large strain continuum problems tha...The possibilities of the particle finite element method(PFEM)for modeling geotechnical problems are increasingly evident.PFEM is a numerical approach to solve large displacement and large strain continuum problems that are beyond the capabilities of classical finite element method(FEM).In PFEM,the computational domain is reconfigured for optimal solution by frequent remeshing and boundary updating.PFEM inherits many concepts,such as a Lagrangian description of continuum,from classic geomechanical FEM.This familiarity with more popular numerical methods facilitates learning and application.This work focuses on G-PFEM,a code specifically developed for the use of PFEM in geotechnical problems.The article has two purposes.The first is to give the reader an overview of the capabilities and main features of the current version of the G-PFEM and the second is to illustrate some of the newer developments of the code.G-PFEM can solve coupled hydro-mechanical static and dynamic problems involving the interaction of solid and/or deformable bodies.Realistic constitutive models for geomaterials are available,including features,such as structure and destructuration,which result in brittle response.The solutions are robust,solidly underpinned by numerical technology including mixedfield formulations,robust and mesh-independent integration of elastoplastic constitutive models and a rigorous and flexible treatment of contact interactions.The novel features presented in this work include the contact domain technique,a natural way to capture contact interactions and impose contact constraints between different continuum bodies,as well as a new simplified formulation for dynamic impact problems.The code performance is showcased by the simulation of several soil-structure interaction problems selected to highlight the novel code features:a rigid footing insertion in soft rock,pipeline insertion and subsequent lateral displacement on over-consolidated clay,screw-pile pull-out and the dynamic impact of a free-falling spherical penetrometer into clay.展开更多
We consider a mathematical model which describes the static frictional contact between a piezoelectric body and a conductive foundation.A non linear electro-elastic constitutive law is used to model the piezoelectric ...We consider a mathematical model which describes the static frictional contact between a piezoelectric body and a conductive foundation.A non linear electro-elastic constitutive law is used to model the piezoelectric material.The unilateral contact is modelled using the Signorini condition,nonlocal Coulomb friction law with slip dependent friction coefficient and a regularized electrical conductivity condition.Existence and uniqueness of a weak solution is established.A finite elements approximation of the problem is presented,a priori error estimates of the solutions are derived and a convergent successive iteration technique is proposed.展开更多
Thermal damage caused by frictional heat of rolling-sliding contact is one of the most important failure forms of wheel and rail. Many studies of wheel-rail frictional heating have been devoted to the temperature fiel...Thermal damage caused by frictional heat of rolling-sliding contact is one of the most important failure forms of wheel and rail. Many studies of wheel-rail frictional heating have been devoted to the temperature field, but few literatures focus on wheel-rail thermal stress caused by frictional heating. However, the wheel-rail creepage is one of important influencing factors of the thermal stress In this paper, a thermo-mechanical coupling model of wheel-rail rolling-sliding contact is developed using thermo-elasto-plastic finite element method. The effect of the wheel-rail elastic creepage on the distribution of heat flux is investigated using the numerical model in which the temperature-dependent material properties are taken into consideration. The moving wheel-rail contact force and the frictional heating are used to simulate the wheel rolling on the rail. The effect of the creepage on the temperature rise, thermal strain, residual stress and residual strain under wheel-rail sliding-rolling contact are investigated. The investigation results show that the thermally affected zone exists mainly in a very thin layer of material near the rail contact surface during the rolling-sliding contact. Both the temperature and thermal strain of rail increase with increasing creepage. The residual stresses induced by the frictional heat in the surface layer of rail appear to be tensile. When the creepage is large, the frictional heat has a significant influence on the residual stresses and residual strains of rail. This paper develops a thermo-meehanical coupling model of wheel-rail rolling-sliding contact, and the obtained results can help to understand the mechanism of wheel/rail frictional thermal fatigue.展开更多
This study provides a comprehensive analysis of collision and impact problems’ numerical solutions, focusing ongeometric, contact, and material nonlinearities, all essential in solving large deformation problems duri...This study provides a comprehensive analysis of collision and impact problems’ numerical solutions, focusing ongeometric, contact, and material nonlinearities, all essential in solving large deformation problems during a collision.The initial discussion revolves around the stress and strain of large deformation during a collision, followedby explanations of the fundamental finite element solution method for addressing such issues. The hourglassmode’s control methods, such as single-point reduced integration and contact-collision algorithms are detailedand implemented within the finite element framework. The paper further investigates the dynamic responseand failure modes of Reinforced Concrete (RC) members under asymmetrical impact using a 3D discrete modelin ABAQUS that treats steel bars and concrete connections as bond slips. The model’s validity was confirmedthrough comparisons with the node-sharing algorithm and system energy relations. Experimental parameterswere varied, including the rigid hammer’s mass and initial velocity, concrete strength, and longitudinal and stirrupreinforcement ratios. Findings indicated that increased hammer mass and velocity escalated RC member damage,while increased reinforcement ratios improved impact resistance. Contrarily, increased concrete strength did notsignificantly reduce lateral displacement when considering strain rate effects. The study also explores materialnonlinearity, examining different materials’ responses to collision-induced forces and stresses, demonstratedthrough an elastic rod impact case study. The paper proposes a damage criterion based on the residual axialload-bearing capacity for assessing damage under the asymmetrical impact, showing a correlation betweendamage degree hammer mass and initial velocity. The results, validated through comparison with theoreticaland analytical solutions, verify the ABAQUS program’s accuracy and reliability in analyzing impact problems,offering valuable insights into collision and impact problems’ nonlinearities and practical strategies for enhancingRC structures’ resilience under dynamic stress.展开更多
文摘The adaptive element techniques of contact problem are studied by means of penalty method, and the error estimators are discussed. Based on error estimators, algorithm of the adaptive element techniques is developed, then the Gauss - Newton iterations are used which allow the nonlinear problem to be transformed into a sequence of linear sub- problems then easily solved. In addition, the algorithm can be applied into the simulation of de -bonding of fiber - reinforced composites.
基金The Project Supported by National Natural Science Foundation of China
文摘Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed finite ele- ment equations,thus resulting in a reduction by half in the dimension of final governing equations.Second,an algorithm composed of contact condition probes and elasto-plastic iterations is utilized to solve the governing equation,which distinguishes two kinds of nonlinearities,and makes the solution unique.In addition,Positive-Negative Sequence Modifica- tion Method is used to condense the finite element equations of each substructure and an analytical integration is intro- duced to determine the elasto-plastic status after each time step or each iteration,hence the computational efficiency is en- hanced to a great extent.Finally,several test and practical examples are pressented showing the validity and versatility of these methods and algorithms.
文摘Based on the analysis of [7] and [10], we present the mixed finite element approximation of the variational inequality resulting from the contact problem in elasticity. The convergence rate of the stress and displacement field are both improved from O（h3/4） to quasi-optimal O（h│logh│^1/4）. If stronger but reasonable regularity is available, the convergence rate can be optimal O（h）.
文摘In this paper, we provide a new mixed finite element approximation of the variational inequality resulting from the unilateral contact problem in elasticity. We use the continuous piecewise P2-P1 finite element to approximate the displacement field and the normal stress component on the contact region. Optimal convergence rates are obtained under the reasonable regularity hypotheses. Numerical example verifies our results.
基金financial support by Severo Ochoa Centre of Excellence (2019-2023) Grant No. CEX2018-000797-Sfunded by MCIN/AEI/10.13039/501100011033+1 种基金research projects BIA2017-84752-RPID2020-119598RB-I00
文摘The possibilities of the particle finite element method(PFEM)for modeling geotechnical problems are increasingly evident.PFEM is a numerical approach to solve large displacement and large strain continuum problems that are beyond the capabilities of classical finite element method(FEM).In PFEM,the computational domain is reconfigured for optimal solution by frequent remeshing and boundary updating.PFEM inherits many concepts,such as a Lagrangian description of continuum,from classic geomechanical FEM.This familiarity with more popular numerical methods facilitates learning and application.This work focuses on G-PFEM,a code specifically developed for the use of PFEM in geotechnical problems.The article has two purposes.The first is to give the reader an overview of the capabilities and main features of the current version of the G-PFEM and the second is to illustrate some of the newer developments of the code.G-PFEM can solve coupled hydro-mechanical static and dynamic problems involving the interaction of solid and/or deformable bodies.Realistic constitutive models for geomaterials are available,including features,such as structure and destructuration,which result in brittle response.The solutions are robust,solidly underpinned by numerical technology including mixedfield formulations,robust and mesh-independent integration of elastoplastic constitutive models and a rigorous and flexible treatment of contact interactions.The novel features presented in this work include the contact domain technique,a natural way to capture contact interactions and impose contact constraints between different continuum bodies,as well as a new simplified formulation for dynamic impact problems.The code performance is showcased by the simulation of several soil-structure interaction problems selected to highlight the novel code features:a rigid footing insertion in soft rock,pipeline insertion and subsequent lateral displacement on over-consolidated clay,screw-pile pull-out and the dynamic impact of a free-falling spherical penetrometer into clay.
文摘We consider a mathematical model which describes the static frictional contact between a piezoelectric body and a conductive foundation.A non linear electro-elastic constitutive law is used to model the piezoelectric material.The unilateral contact is modelled using the Signorini condition,nonlocal Coulomb friction law with slip dependent friction coefficient and a regularized electrical conductivity condition.Existence and uniqueness of a weak solution is established.A finite elements approximation of the problem is presented,a priori error estimates of the solutions are derived and a convergent successive iteration technique is proposed.
基金supported by National Natural Science Foundation of China(Grant Nos.51175438,U1134202)National Science and Technology Support Program of China(Grant No.2009BAG12A01)Program for New Century Excellent Talents in University of China(Grant No.NCET-08-0824)
文摘Thermal damage caused by frictional heat of rolling-sliding contact is one of the most important failure forms of wheel and rail. Many studies of wheel-rail frictional heating have been devoted to the temperature field, but few literatures focus on wheel-rail thermal stress caused by frictional heating. However, the wheel-rail creepage is one of important influencing factors of the thermal stress In this paper, a thermo-mechanical coupling model of wheel-rail rolling-sliding contact is developed using thermo-elasto-plastic finite element method. The effect of the wheel-rail elastic creepage on the distribution of heat flux is investigated using the numerical model in which the temperature-dependent material properties are taken into consideration. The moving wheel-rail contact force and the frictional heating are used to simulate the wheel rolling on the rail. The effect of the creepage on the temperature rise, thermal strain, residual stress and residual strain under wheel-rail sliding-rolling contact are investigated. The investigation results show that the thermally affected zone exists mainly in a very thin layer of material near the rail contact surface during the rolling-sliding contact. Both the temperature and thermal strain of rail increase with increasing creepage. The residual stresses induced by the frictional heat in the surface layer of rail appear to be tensile. When the creepage is large, the frictional heat has a significant influence on the residual stresses and residual strains of rail. This paper develops a thermo-meehanical coupling model of wheel-rail rolling-sliding contact, and the obtained results can help to understand the mechanism of wheel/rail frictional thermal fatigue.
基金the authority of the National Natural Science Foundation of China(Grant Nos.52178168 and 51378427)for financing this research work and several ongoing research projects related to structural impact performance.
文摘This study provides a comprehensive analysis of collision and impact problems’ numerical solutions, focusing ongeometric, contact, and material nonlinearities, all essential in solving large deformation problems during a collision.The initial discussion revolves around the stress and strain of large deformation during a collision, followedby explanations of the fundamental finite element solution method for addressing such issues. The hourglassmode’s control methods, such as single-point reduced integration and contact-collision algorithms are detailedand implemented within the finite element framework. The paper further investigates the dynamic responseand failure modes of Reinforced Concrete (RC) members under asymmetrical impact using a 3D discrete modelin ABAQUS that treats steel bars and concrete connections as bond slips. The model’s validity was confirmedthrough comparisons with the node-sharing algorithm and system energy relations. Experimental parameterswere varied, including the rigid hammer’s mass and initial velocity, concrete strength, and longitudinal and stirrupreinforcement ratios. Findings indicated that increased hammer mass and velocity escalated RC member damage,while increased reinforcement ratios improved impact resistance. Contrarily, increased concrete strength did notsignificantly reduce lateral displacement when considering strain rate effects. The study also explores materialnonlinearity, examining different materials’ responses to collision-induced forces and stresses, demonstratedthrough an elastic rod impact case study. The paper proposes a damage criterion based on the residual axialload-bearing capacity for assessing damage under the asymmetrical impact, showing a correlation betweendamage degree hammer mass and initial velocity. The results, validated through comparison with theoreticaland analytical solutions, verify the ABAQUS program’s accuracy and reliability in analyzing impact problems,offering valuable insights into collision and impact problems’ nonlinearities and practical strategies for enhancingRC structures’ resilience under dynamic stress.
