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.展开更多
Through defining slide yield function and floating potential function of thermo-contact surface, the complementary equation of thermo-contact boundary has been reached, the fundamental equations to solve 3D thermo-con...Through defining slide yield function and floating potential function of thermo-contact surface, the complementary equation of thermo-contact boundary has been reached, the fundamental equations to solve 3D thermo-contact coupled problem have been listed. On this foundation, the finite element equation and definite solution condition of contact heat transfer have been given out. Based on virtual work principle and contact element technology, the finite element equation of 3D elastic contact system has been deduced under the effect of thermal stress. The pseudo load brought by contact gap have been introduced into this equation in order to reflect the contact state change. During iteration, once contact rigidity matrix is formed, it won’t change, which will make calculation reduce greatly.展开更多
In this paper, the linear finite element approximation to the elastic contact problem with curved contact boundary is considered. The error bound O(h[sup ?]) is obtained with requirements of two times continuously dif...In this paper, the linear finite element approximation to the elastic contact problem with curved contact boundary is considered. The error bound O(h[sup ?]) is obtained with requirements of two times continuously differentiable for contact boundary and the usual regular triangulation, while I.Hlavacek et. al. Obtained the error bound O(h[sup ?]) with requirements of three times continuously differentiable for contact boundary and extra regularities of triangulation (c.f. [2]). [ABSTRACT FROM AUTHOR]展开更多
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）.展开更多
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.展开更多
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.展开更多
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.展开更多
The paper presents the formulation and approximation of a static thermoelasticity problem that describes bilateral frictional contact between a deformable body and a rigid foundation. The friction is in the form of a ...The paper presents the formulation and approximation of a static thermoelasticity problem that describes bilateral frictional contact between a deformable body and a rigid foundation. The friction is in the form of a nonmonotone and multivalued law. The coupling effect of the problem is neglected. Therefore, the thermic part of the problem is considered independently on the elasticity problem. For the displacement vector, we formulate one substationary problem for a non-convex, locally Lipschitz continuous functional representing the total potential energy of the body. All problems formulated in the paper are approximated with the finite element method.展开更多
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.展开更多
Based on the theory of reliability-based structural shape optimization, exact expressions of the sensibility using the stochastic finite element method for contact problems were derived in detail, and the basic steps ...Based on the theory of reliability-based structural shape optimization, exact expressions of the sensibility using the stochastic finite element method for contact problems were derived in detail, and the basic steps of structural optimization were given. A coattail-type tenon/mortise of an aero-engine was optimized. In this model, the maximum equivalent stress of the nodes on the boundary of the tenon was the objective function; the width of tooth’s neck and the side surface’s slope angle of a tenon were design variables, with constraints of tension stress, extrusion stress and reliability index. The result showed that the distributions of the contact pressure between tenon and mortise, the equivalence stress and reliability index were more reasonable. It validates the correctness of the optimization model and the reliability-based structural shape optimization, and provides valuable references for structural design of the tenon/mortise.展开更多
A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations ...A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations are derived in detail for the interface element.All the involved matrices are of the same form as those of a solid element,which makes the incorporation of the model into a finite element program straightforward.Three examples are then numerically simulated using the interface element.Reasonable results confirm the correctness of the proposed model and motivate its application in hydromechanical contact problems in the future.展开更多
With rollers as elastic bodies and workpieces as elastoplastic bodies,the rolling problem can be viewed as a friction elastic-plastic contact problem.With fewer assumptions in the simulation of strip-rolling process,a...With rollers as elastic bodies and workpieces as elastoplastic bodies,the rolling problem can be viewed as a friction elastic-plastic contact problem.With fewer assumptions in the simulation of strip-rolling process,a boundary element method（BEM）for two-dimensional elastoplastic finite strain and finite deformation analysis of contact problems with friction was presented.All the equations for contact problems,which include multi-nonlinearities,were obtained.Incremental and iterative procedures were used to find contact pressure and friction stress.Moreover,initial strain rate algorithm and work-hardening material behavior can be assumed in the plastic analysis.Several examples were presented,and the results of contact pressure and friction stress were in excellent agreement with those of analysis.展开更多
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.
文摘Through defining slide yield function and floating potential function of thermo-contact surface, the complementary equation of thermo-contact boundary has been reached, the fundamental equations to solve 3D thermo-contact coupled problem have been listed. On this foundation, the finite element equation and definite solution condition of contact heat transfer have been given out. Based on virtual work principle and contact element technology, the finite element equation of 3D elastic contact system has been deduced under the effect of thermal stress. The pseudo load brought by contact gap have been introduced into this equation in order to reflect the contact state change. During iteration, once contact rigidity matrix is formed, it won’t change, which will make calculation reduce greatly.
文摘In this paper, the linear finite element approximation to the elastic contact problem with curved contact boundary is considered. The error bound O(h[sup ?]) is obtained with requirements of two times continuously differentiable for contact boundary and the usual regular triangulation, while I.Hlavacek et. al. Obtained the error bound O(h[sup ?]) with requirements of three times continuously differentiable for contact boundary and extra regularities of triangulation (c.f. [2]). [ABSTRACT FROM AUTHOR]
文摘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）.
基金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.
文摘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 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.
基金supported by the Minisitry of Science of the Republic of Serbia (No. 144005)
文摘The paper presents the formulation and approximation of a static thermoelasticity problem that describes bilateral frictional contact between a deformable body and a rigid foundation. The friction is in the form of a nonmonotone and multivalued law. The coupling effect of the problem is neglected. Therefore, the thermic part of the problem is considered independently on the elasticity problem. For the displacement vector, we formulate one substationary problem for a non-convex, locally Lipschitz continuous functional representing the total potential energy of the body. All problems formulated in the paper are approximated with the finite element method.
文摘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.
基金National Natural Science Foundation of China (5 9875 0 3 7)
文摘Based on the theory of reliability-based structural shape optimization, exact expressions of the sensibility using the stochastic finite element method for contact problems were derived in detail, and the basic steps of structural optimization were given. A coattail-type tenon/mortise of an aero-engine was optimized. In this model, the maximum equivalent stress of the nodes on the boundary of the tenon was the objective function; the width of tooth’s neck and the side surface’s slope angle of a tenon were design variables, with constraints of tension stress, extrusion stress and reliability index. The result showed that the distributions of the contact pressure between tenon and mortise, the equivalence stress and reliability index were more reasonable. It validates the correctness of the optimization model and the reliability-based structural shape optimization, and provides valuable references for structural design of the tenon/mortise.
基金supported by the Innovation Plan for Postgraduate Students sponsored by the Education Department of Jiangsu Province,China (CX08B 107Z)
文摘A hydromechanical interface element is proposed for the consideration of the hydraulic-mechanical coupling effect along the interface.The fully coupled governing equations and the relevant finite element formulations are derived in detail for the interface element.All the involved matrices are of the same form as those of a solid element,which makes the incorporation of the model into a finite element program straightforward.Three examples are then numerically simulated using the interface element.Reasonable results confirm the correctness of the proposed model and motivate its application in hydromechanical contact problems in the future.
基金Item Sponsored by National Natural Science Foundation of China(50575155)
文摘With rollers as elastic bodies and workpieces as elastoplastic bodies,the rolling problem can be viewed as a friction elastic-plastic contact problem.With fewer assumptions in the simulation of strip-rolling process,a boundary element method（BEM）for two-dimensional elastoplastic finite strain and finite deformation analysis of contact problems with friction was presented.All the equations for contact problems,which include multi-nonlinearities,were obtained.Incremental and iterative procedures were used to find contact pressure and friction stress.Moreover,initial strain rate algorithm and work-hardening material behavior can be assumed in the plastic analysis.Several examples were presented,and the results of contact pressure and friction stress were in excellent agreement with those of analysis.
基金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.
