In this paper, with the use of the friction problem in elasticity as the background, the existence and uniqueness for the solution of the nonlinear, indifferentiable mixed variational inequality are discussed. Its cor...In this paper, with the use of the friction problem in elasticity as the background, the existence and uniqueness for the solution of the nonlinear, indifferentiable mixed variational inequality are discussed. Its corresponding boundary variational inequality and the existence and uniqueness, of solution are given. This provides the theoretical basis for using boundary element method to solve the mixed variational inequality.展开更多
A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle ...A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.展开更多
Three dimensional frictional contact problems are formulated as linear complementarity problems based on the parametric variational principle. Two aggregate-functionbased algorithms for solving complementarity problem...Three dimensional frictional contact problems are formulated as linear complementarity problems based on the parametric variational principle. Two aggregate-functionbased algorithms for solving complementarity problems are proposed. One is called the self-adjusting interior point algorithm, the other is called the aggregate function smoothing algorithm. Numerical experiment shows the efficiency of the proposed two algorithms.展开更多
Based on the theory and technique of nonlinear geometric field theory of continuum, a more general incremental variational equation for elastic and plastic large deformation in co-moving coordinate is established in t...Based on the theory and technique of nonlinear geometric field theory of continuum, a more general incremental variational equation for elastic and plastic large deformation in co-moving coordinate is established in this paper. An expression for two and three-ditnensicnal continua is derived, and the incremental variational equation for large deformation of changing boundary contact and the variational inequality in rate form tire obtained, which provides the theoretical basis for the computation of elastic-plastic large deformation contact problem with friction.展开更多
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.展开更多
The formulation of boundary element method for handling contact problems with friction and the technique for high speed contact analysis are presented. This formulation is based on the idea of modifying the length of...The formulation of boundary element method for handling contact problems with friction and the technique for high speed contact analysis are presented. This formulation is based on the idea of modifying the length of contact elements without altering the total number of elements. The high precision of solution and high speed analysis are verified according to the results of conventional method and analysis method.展开更多
An approach is proposed for modeling and anal- yses of rigid multibody systems with frictional translation joints and driving constraints. The geometric constraints of translational joints with small clearance are tre...An approach is proposed for modeling and anal- yses of rigid multibody systems with frictional translation joints and driving constraints. The geometric constraints of translational joints with small clearance are treated as bilat- eral constraints by neglecting the impact between sliders and guides. Firstly, the normal forces acting on sliders, the driv- ing constraint forces （or moments） and the constraint forces of smooth revolute joints are all described by complementary conditions. The frictional contacts are characterized by a set- valued force law of Coulomb＇s dry friction. Combined with the theory of the horizontal linear complementarity problem （HLCP）, an event-driven scheme is used to detect the transi- tions of the contact situation between sliders and guides, and the stick-slip transitions of sliders, respectively. And then, all constraint forces in the system can be computed easily. Secondly, the dynamic equations of multibody systems are written at the acceleration-force level by the Lagrange multiplier technique, and the Baumgarte stabilization method is used to reduce the constraint drift. Finally, a numerical example is given to show some non-smooth dynamical behaviors of the studied system. The obtained results validate the feasibility of algorithm and the effect of constraint stabilization.展开更多
Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction ...Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints is presented. In particular, a wheeled multi-body system is considered. Here, the state transition of stick-slip between wheel and ground is transformed into a nonlinear complementarity problem （NCP）. An iterative algorithm for solving the NCP is then presented using an event-driven method. Dynamical equations of the multi-body system with holonomic and nonholonomic constraints are given using Routh equations and a con- straint stabilization method. Finally, an example is used to test the proposed numerical method. The results show some dynamical behaviors of the wheeled multi-body system and its constraint stabilization effects.展开更多
The objective of this study is to investigate the effects of the Coulomb dry friction model versus the modified Coulomb friction model on the dynamic behavior of the slider-crank mechanism with a revolute clearance jo...The objective of this study is to investigate the effects of the Coulomb dry friction model versus the modified Coulomb friction model on the dynamic behavior of the slider-crank mechanism with a revolute clearance joint. The normal and tangential forces acting on the contact points between the journal and the bearing are described by using a Hertzian-based contact force model and the Coulomb friction models, respectively.The dynamic equations of the mechanism are derived based on the Lagrange equations of the first kind and the Baumgarte stabilization method. The frictional force is solved via the linear complementarity problem(LCP) algorithm and the trial-and-error algorithm.Finally, three numerical examples are given to show the influence of the two Coulomb friction models on the dynamic behavior of the mechanism. Numerical results show that due to the stick friction, the slider-crank mechanism may exhibit stick-slip motion and can balance at some special positions, while the mechanism with ideal joints cannot.展开更多
The Lagrange-I equations and measure differential equations for multibody systems with unilateral and bilateral constraints are constructed. For bilateral constraints, frictional forces and their impulses contain the ...The Lagrange-I equations and measure differential equations for multibody systems with unilateral and bilateral constraints are constructed. For bilateral constraints, frictional forces and their impulses contain the products of the filled-in relay function induced by Coulomb friction and the absolute values of normal constraint reactions. With the time-stepping impulse-velocity scheme, the measure differential equations are discretized. The equations of horizontal linear complementarity problems （HLCPs）, which are used to compute the impulses, are constructed by decomposing the absolute function and the filled-in relay function. These HLCP equations degenerate into equations of LCPs for frictional unilateral constraints, or HLCPs for frictional bilateral constraints. Finally, a numerical simulation for multibody systems with both unilateral and bilateral constraints is presented.展开更多
Our aim is to investigate the solutions to the boundary layer problem of a power-law non-Newtonian fluid along an impermeable sheet moving with a constant velocity in an otherwise quiescent fluid environment. In the a...Our aim is to investigate the solutions to the boundary layer problem of a power-law non-Newtonian fluid along an impermeable sheet moving with a constant velocity in an otherwise quiescent fluid environment. In the absence of an exact solution in closed form, numerical solutions for the velocity distribution in the boundary layer for different power exponents will be presented. Our goal is to give an iterative transformation method for the determination of the skin friction parameter and the boundary layer thickness for different parameter values and the dependence of the skin friction parameter and the boundary layer thickness on the power exponent are examined.展开更多
文摘In this paper, with the use of the friction problem in elasticity as the background, the existence and uniqueness for the solution of the nonlinear, indifferentiable mixed variational inequality are discussed. Its corresponding boundary variational inequality and the existence and uniqueness, of solution are given. This provides the theoretical basis for using boundary element method to solve the mixed variational inequality.
文摘A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.
基金The project supported by the National Natural Science foundation of china(10225212,50178016.10302007)the National Kev Basic Research Special Foundation and the Ministry of Education of China
文摘Three dimensional frictional contact problems are formulated as linear complementarity problems based on the parametric variational principle. Two aggregate-functionbased algorithms for solving complementarity problems are proposed. One is called the self-adjusting interior point algorithm, the other is called the aggregate function smoothing algorithm. Numerical experiment shows the efficiency of the proposed two algorithms.
文摘Based on the theory and technique of nonlinear geometric field theory of continuum, a more general incremental variational equation for elastic and plastic large deformation in co-moving coordinate is established in this paper. An expression for two and three-ditnensicnal continua is derived, and the incremental variational equation for large deformation of changing boundary contact and the variational inequality in rate form tire obtained, which provides the theoretical basis for the computation of elastic-plastic large deformation contact problem with friction.
基金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.
文摘The formulation of boundary element method for handling contact problems with friction and the technique for high speed contact analysis are presented. This formulation is based on the idea of modifying the length of contact elements without altering the total number of elements. The high precision of solution and high speed analysis are verified according to the results of conventional method and analysis method.
基金supported by the National Natural Science Foundation of China(11372018 and 11172019)
文摘An approach is proposed for modeling and anal- yses of rigid multibody systems with frictional translation joints and driving constraints. The geometric constraints of translational joints with small clearance are treated as bilat- eral constraints by neglecting the impact between sliders and guides. Firstly, the normal forces acting on sliders, the driv- ing constraint forces （or moments） and the constraint forces of smooth revolute joints are all described by complementary conditions. The frictional contacts are characterized by a set- valued force law of Coulomb＇s dry friction. Combined with the theory of the horizontal linear complementarity problem （HLCP）, an event-driven scheme is used to detect the transi- tions of the contact situation between sliders and guides, and the stick-slip transitions of sliders, respectively. And then, all constraint forces in the system can be computed easily. Secondly, the dynamic equations of multibody systems are written at the acceleration-force level by the Lagrange multiplier technique, and the Baumgarte stabilization method is used to reduce the constraint drift. Finally, a numerical example is given to show some non-smooth dynamical behaviors of the studied system. The obtained results validate the feasibility of algorithm and the effect of constraint stabilization.
基金Project supported by the National Natural Science Foundation of China(Nos.11372018 and 11572018)
文摘Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints is presented. In particular, a wheeled multi-body system is considered. Here, the state transition of stick-slip between wheel and ground is transformed into a nonlinear complementarity problem （NCP）. An iterative algorithm for solving the NCP is then presented using an event-driven method. Dynamical equations of the multi-body system with holonomic and nonholonomic constraints are given using Routh equations and a con- straint stabilization method. Finally, an example is used to test the proposed numerical method. The results show some dynamical behaviors of the wheeled multi-body system and its constraint stabilization effects.
基金Project supported by the National Natural Science Foundation of China(No.11772021)
文摘The objective of this study is to investigate the effects of the Coulomb dry friction model versus the modified Coulomb friction model on the dynamic behavior of the slider-crank mechanism with a revolute clearance joint. The normal and tangential forces acting on the contact points between the journal and the bearing are described by using a Hertzian-based contact force model and the Coulomb friction models, respectively.The dynamic equations of the mechanism are derived based on the Lagrange equations of the first kind and the Baumgarte stabilization method. The frictional force is solved via the linear complementarity problem(LCP) algorithm and the trial-and-error algorithm.Finally, three numerical examples are given to show the influence of the two Coulomb friction models on the dynamic behavior of the mechanism. Numerical results show that due to the stick friction, the slider-crank mechanism may exhibit stick-slip motion and can balance at some special positions, while the mechanism with ideal joints cannot.
基金supported by the National Natural Science Foundation of China (10672007)
文摘The Lagrange-I equations and measure differential equations for multibody systems with unilateral and bilateral constraints are constructed. For bilateral constraints, frictional forces and their impulses contain the products of the filled-in relay function induced by Coulomb friction and the absolute values of normal constraint reactions. With the time-stepping impulse-velocity scheme, the measure differential equations are discretized. The equations of horizontal linear complementarity problems （HLCPs）, which are used to compute the impulses, are constructed by decomposing the absolute function and the filled-in relay function. These HLCP equations degenerate into equations of LCPs for frictional unilateral constraints, or HLCPs for frictional bilateral constraints. Finally, a numerical simulation for multibody systems with both unilateral and bilateral constraints is presented.
文摘Our aim is to investigate the solutions to the boundary layer problem of a power-law non-Newtonian fluid along an impermeable sheet moving with a constant velocity in an otherwise quiescent fluid environment. In the absence of an exact solution in closed form, numerical solutions for the velocity distribution in the boundary layer for different power exponents will be presented. Our goal is to give an iterative transformation method for the determination of the skin friction parameter and the boundary layer thickness for different parameter values and the dependence of the skin friction parameter and the boundary layer thickness on the power exponent are examined.
