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.展开更多
This work presents a numerical methodology for modeling the Winkler supports and nonlinear conditions by proposing new boundary conditions. For the boundary conditions of Winkler support model, the surface tractions a...This work presents a numerical methodology for modeling the Winkler supports and nonlinear conditions by proposing new boundary conditions. For the boundary conditions of Winkler support model, the surface tractions and the displacements normal to the surface of the solid are unknown, but their relationship is known by means of the ballast coefficient, whereas for nonlinear boundary conditions, the displacements normal to the boundary of the solid are zero in the positive direction but are allowed in the negative direction. In those zones, detachments of nodes might appear, leading to a nonlinearity, because the number of nodes that remain fixed or of the detached ones (under tensile tractions) is unknown. The proposed methodology is applied to the 3D elastic receding contact problem using the boundary element method. The surface t r actions and the displacements of the common int erface bet ween the two solids in contac t under the influence of different supports are calculated as well as the boundary zone of the solid where the new boundary conditions are applied. The problem is solved by a double-iterative met hod, so in the final solut ion, t here are no t r act ions or pene trations between the two solids or at the boundary of the solid where the nonlinear boundary conditions are Simula ted. The effectiveness of the proposed method is verified by examples.展开更多
A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the co...A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.展开更多
Heat transfer and deformation of initial solidification shell in soft contact continuous casting mold under high frequency electromagnetic field were analyzed using numerical simulation method; the relative electromag...Heat transfer and deformation of initial solidification shell in soft contact continuous casting mold under high frequency electromagnetic field were analyzed using numerical simulation method; the relative electromagnetic parameters were obtained from the previous studies. Owing to the induction heating of a high frequency electromagnetic field （20 kHz）, the thickness of initial solidification shell decreases, and the temperature of strand surface and slit copper mold increases when compared with the case without the electromagnetic filed. The viscosity of flux de- creases because of the induction heating of the high frequency electromagnetic field, and the dimension of the flux channel increases with electromagnetic pressure; thus, the deformation behavior of initial solidification shell was different before and after the action of high frequency electromagnetic field. Furthermore, the abatement mechanism of oscillation marks under high frequency electromagnetic field was explained.展开更多
The focus of the current contribution is on the development of the unified geometrical formulation of contact algorithms in a covariant form for various geometrical situations of contacting bodies leading to contact p...The focus of the current contribution is on the development of the unified geometrical formulation of contact algorithms in a covariant form for various geometrical situations of contacting bodies leading to contact pairs: surface-to-surface, line-to-surface, point-to-surface, line-to-line, point-to-line, point-to-point. The construction of the corresponding computational contact algorithms are considered in accordance with the geometry of contact bodies in a covariant form. These forms can be easily discredited within finite element methods independently of order of approximation and, therefore, the result is straightforwardly applied within iso-geometric finite element methods. This approach is recently became known as geometrically exact theory of contact interaction [10]. Application for contact between bodies with iso- and anisotropic surface, for contact between cables and curvilinear beams as well as recent development for contact between cables and bodies is straightforward. Recent developments include the improvement of the curve-to-surface (deformable) contact algorithm.展开更多
To improve the application of discrete element models(DEM)to the design of agricultural crushers,in this study a new highly accurate model is elaborated.The model takes into account the fiber structure,porous nature o...To improve the application of discrete element models(DEM)to the design of agricultural crushers,in this study a new highly accurate model is elaborated.The model takes into account the fiber structure,porous nature of the material and the leaf sheath coating structure.Dedicated experimental tests are conducted to determine the required“intrinsic”and basic contact parameters of the considered banana straw materials.A large number of bonding parameters are examined in relation to the particle aggregation model in order to characterize different actual banana straws.Using the particle surface energy contact model,the viscosity characteristics of the crushed material are determined together with the related stacking angle(considered as the main response factor).Through single factor experiment analysis,it is found that when the surface energy is 0.9 J·m-2,the relative error between simulations and physical experiments is 5.288%.展开更多
Discontinuous deformation problems are common in rock engineering. Numerical analysis methods based on system models of the discrete body can better solve these problems. One of the most effective solutions is discont...Discontinuous deformation problems are common in rock engineering. Numerical analysis methods based on system models of the discrete body can better solve these problems. One of the most effective solutions is discontinuous deformation analysis （DDA） method, but the DDA method brings about rock embedding problems when it uses the strain assumption in elastic deformation and adopts virtual springs to simulate the contact problems. The multi-body finite element method （FEM） proposed in this paper can solve the problems of contact and deformation of blocks very well because it integrates the FEM and multi-body system dynamics theory. It is therefore a complete method for solving discontinuous deformation problems through balance equations of the contact surface and for simulating the displacement of whole blocks. In this study, this method was successfully used for deformation analysis of underground caverns in stratified rock. The simulation results indicate that the multi-body FEM can show contact forces and the stress states on contact surfaces better than DDA, and that the results calculated with the multi-body FEM are more consistent with engineering practice than those calculated with DDA method.展开更多
In the numerical study of rough surfaces in contact problem, the flexible body beneath the roughness is commonly assumed as a half-space or a half-plane. The surface displacement on the boundary, the displacement comp...In the numerical study of rough surfaces in contact problem, the flexible body beneath the roughness is commonly assumed as a half-space or a half-plane. The surface displacement on the boundary, the displacement components and state of stress inside the half-space can be determined through the convolution of the traction and the corresponding influence function in a closed-form. The influence function is often represented by the Boussinesq-Cerruti solution and the Flamant solution for three-dimensional elasticity and plane strain/stress, respectively. In this study, we rigorously show that any numerical model using the above mentioned half-space solution is a special form of the boundary element method(BEM). The boundary integral equations(BIEs) in the BEM is simplified to the Flamant solution when the domain is strictly a half-plane for the plane strain/stress condition. Similarly, the BIE is degraded to the Boussinesq-Cerruti solution if the domain is strictly a half-space. Therefore, the numerical models utilizing these closed-form influence functions are the special BEM where the domain is a half-space(or a half-plane). This analytical work sheds some light on how to accurately simulate the non-half-space contact problem using the BEM.展开更多
The adhesion between the mining machine and the deep-sea sediments will significantly affect the driving performance of the mining machine in the deep-sea environment.When the mining machine and the deep-sea sediment ...The adhesion between the mining machine and the deep-sea sediments will significantly affect the driving performance of the mining machine in the deep-sea environment.When the mining machine and the deep-sea sediment interaction simulation was carried out,the accuracy of the particle interaction parameters will directly affect the simulation results.This study proposed a method to systematically calibrate the interaction parameters between deep-sea sediment and grouser through the combination of experiment and simulation.The uniaxial compression test and macro adhesion test and corresponding discrete element numerical simulation were carried out,modifying the contact parameters until the simulation results are close to the experimental results.Then the micro-parameters of the JKR adhesion contact model were back calibrated with the test results,and the contact parameters between soil particle-soil particle and soil particle-metal are calibrated.Besides,the adhesion test shows that the adhesion forces were ranked in the order of 5052<STi80<TA2<TC4 under the same surface roughness,which indicates the aluminum alloy 5052 has the best anti-adhesion performance.The relationship between surface adhesion force and microscopic contact parameters was studied by discrete element numerical simulation,and the result shows that the coefficient of static friction and the coefficient of rolling friction has little effect on adhesion force.While it is mainly affected by the coefficient of restitution and surface energy,the surface adhesion force tends to decrease with the increase of the coefficient of restitution and increase with the growth of surface energy.The obtained parameters of soil particle to soil particle and soil particle to metal affecting the adhesion can contribute to the design optimization for the grouser of mining machines to decrease surface adhesion and enhance its movability and mining efficiency.展开更多
Recently, gears of high strength, reliability, and surface-damage-resistant under severe service conditions are required to achieve the weight saving and downsizing of a product. For the high-speed condition in partic...Recently, gears of high strength, reliability, and surface-damage-resistant under severe service conditions are required to achieve the weight saving and downsizing of a product. For the high-speed condition in particular, it is important to understand the influence of the surface properties on the scuffing resistance. If the effective surface profile to improve the lubrication property was found, the metal surfaces could be obtained with both surface strength and surface lubricity. Herein, the influence of surface properties modified with fine shot peening, which can form the arbitrary surface profile, on the scuffing resistance in the rolling-sliding contact machine element, was investigated. The scuffing test was performed using a two-cylinder rolling contact test machine. In a specific sliding, a faster roller of 60% and a sliding velocity of 1.75 m/s were utilized. The scuffing test results with shot-peened test rollers and those with non-shot-peened test roller were compared. The influence of the surface roughness of the shot-peened test roller was also discussed. We found that the shot-peened roller had a better scuffing resistance compared with the roller without the shot-peening process.展开更多
文摘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.
文摘This work presents a numerical methodology for modeling the Winkler supports and nonlinear conditions by proposing new boundary conditions. For the boundary conditions of Winkler support model, the surface tractions and the displacements normal to the surface of the solid are unknown, but their relationship is known by means of the ballast coefficient, whereas for nonlinear boundary conditions, the displacements normal to the boundary of the solid are zero in the positive direction but are allowed in the negative direction. In those zones, detachments of nodes might appear, leading to a nonlinearity, because the number of nodes that remain fixed or of the detached ones (under tensile tractions) is unknown. The proposed methodology is applied to the 3D elastic receding contact problem using the boundary element method. The surface t r actions and the displacements of the common int erface bet ween the two solids in contac t under the influence of different supports are calculated as well as the boundary zone of the solid where the new boundary conditions are applied. The problem is solved by a double-iterative met hod, so in the final solut ion, t here are no t r act ions or pene trations between the two solids or at the boundary of the solid where the nonlinear boundary conditions are Simula ted. The effectiveness of the proposed method is verified by examples.
基金The project supported by the National Natural Science Foundation of China (19772025)
文摘A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.
基金Item Sponsored by National Natural Science Foundation of China (59734080)
文摘Heat transfer and deformation of initial solidification shell in soft contact continuous casting mold under high frequency electromagnetic field were analyzed using numerical simulation method; the relative electromagnetic parameters were obtained from the previous studies. Owing to the induction heating of a high frequency electromagnetic field （20 kHz）, the thickness of initial solidification shell decreases, and the temperature of strand surface and slit copper mold increases when compared with the case without the electromagnetic filed. The viscosity of flux de- creases because of the induction heating of the high frequency electromagnetic field, and the dimension of the flux channel increases with electromagnetic pressure; thus, the deformation behavior of initial solidification shell was different before and after the action of high frequency electromagnetic field. Furthermore, the abatement mechanism of oscillation marks under high frequency electromagnetic field was explained.
文摘The focus of the current contribution is on the development of the unified geometrical formulation of contact algorithms in a covariant form for various geometrical situations of contacting bodies leading to contact pairs: surface-to-surface, line-to-surface, point-to-surface, line-to-line, point-to-line, point-to-point. The construction of the corresponding computational contact algorithms are considered in accordance with the geometry of contact bodies in a covariant form. These forms can be easily discredited within finite element methods independently of order of approximation and, therefore, the result is straightforwardly applied within iso-geometric finite element methods. This approach is recently became known as geometrically exact theory of contact interaction [10]. Application for contact between bodies with iso- and anisotropic surface, for contact between cables and curvilinear beams as well as recent development for contact between cables and bodies is straightforward. Recent developments include the improvement of the curve-to-surface (deformable) contact algorithm.
基金Youth Fund of National Natural Science Foundation of China,61801175,JJ,https://www.nsfc.gov.cn/.
文摘To improve the application of discrete element models(DEM)to the design of agricultural crushers,in this study a new highly accurate model is elaborated.The model takes into account the fiber structure,porous nature of the material and the leaf sheath coating structure.Dedicated experimental tests are conducted to determine the required“intrinsic”and basic contact parameters of the considered banana straw materials.A large number of bonding parameters are examined in relation to the particle aggregation model in order to characterize different actual banana straws.Using the particle surface energy contact model,the viscosity characteristics of the crushed material are determined together with the related stacking angle(considered as the main response factor).Through single factor experiment analysis,it is found that when the surface energy is 0.9 J·m-2,the relative error between simulations and physical experiments is 5.288%.
文摘Discontinuous deformation problems are common in rock engineering. Numerical analysis methods based on system models of the discrete body can better solve these problems. One of the most effective solutions is discontinuous deformation analysis （DDA） method, but the DDA method brings about rock embedding problems when it uses the strain assumption in elastic deformation and adopts virtual springs to simulate the contact problems. The multi-body finite element method （FEM） proposed in this paper can solve the problems of contact and deformation of blocks very well because it integrates the FEM and multi-body system dynamics theory. It is therefore a complete method for solving discontinuous deformation problems through balance equations of the contact surface and for simulating the displacement of whole blocks. In this study, this method was successfully used for deformation analysis of underground caverns in stratified rock. The simulation results indicate that the multi-body FEM can show contact forces and the stress states on contact surfaces better than DDA, and that the results calculated with the multi-body FEM are more consistent with engineering practice than those calculated with DDA method.
文摘In the numerical study of rough surfaces in contact problem, the flexible body beneath the roughness is commonly assumed as a half-space or a half-plane. The surface displacement on the boundary, the displacement components and state of stress inside the half-space can be determined through the convolution of the traction and the corresponding influence function in a closed-form. The influence function is often represented by the Boussinesq-Cerruti solution and the Flamant solution for three-dimensional elasticity and plane strain/stress, respectively. In this study, we rigorously show that any numerical model using the above mentioned half-space solution is a special form of the boundary element method(BEM). The boundary integral equations(BIEs) in the BEM is simplified to the Flamant solution when the domain is strictly a half-plane for the plane strain/stress condition. Similarly, the BIE is degraded to the Boussinesq-Cerruti solution if the domain is strictly a half-space. Therefore, the numerical models utilizing these closed-form influence functions are the special BEM where the domain is a half-space(or a half-plane). This analytical work sheds some light on how to accurately simulate the non-half-space contact problem using the BEM.
基金Project(12072309)supported by the National Natural Science Foundation of ChinaProject(19B546)supported by the Education Department Foundation of Hunan Province,ChinaProject(2019RS1059)supported by the Hunan Innovative Province Construction Project,China。
文摘The adhesion between the mining machine and the deep-sea sediments will significantly affect the driving performance of the mining machine in the deep-sea environment.When the mining machine and the deep-sea sediment interaction simulation was carried out,the accuracy of the particle interaction parameters will directly affect the simulation results.This study proposed a method to systematically calibrate the interaction parameters between deep-sea sediment and grouser through the combination of experiment and simulation.The uniaxial compression test and macro adhesion test and corresponding discrete element numerical simulation were carried out,modifying the contact parameters until the simulation results are close to the experimental results.Then the micro-parameters of the JKR adhesion contact model were back calibrated with the test results,and the contact parameters between soil particle-soil particle and soil particle-metal are calibrated.Besides,the adhesion test shows that the adhesion forces were ranked in the order of 5052<STi80<TA2<TC4 under the same surface roughness,which indicates the aluminum alloy 5052 has the best anti-adhesion performance.The relationship between surface adhesion force and microscopic contact parameters was studied by discrete element numerical simulation,and the result shows that the coefficient of static friction and the coefficient of rolling friction has little effect on adhesion force.While it is mainly affected by the coefficient of restitution and surface energy,the surface adhesion force tends to decrease with the increase of the coefficient of restitution and increase with the growth of surface energy.The obtained parameters of soil particle to soil particle and soil particle to metal affecting the adhesion can contribute to the design optimization for the grouser of mining machines to decrease surface adhesion and enhance its movability and mining efficiency.
文摘Recently, gears of high strength, reliability, and surface-damage-resistant under severe service conditions are required to achieve the weight saving and downsizing of a product. For the high-speed condition in particular, it is important to understand the influence of the surface properties on the scuffing resistance. If the effective surface profile to improve the lubrication property was found, the metal surfaces could be obtained with both surface strength and surface lubricity. Herein, the influence of surface properties modified with fine shot peening, which can form the arbitrary surface profile, on the scuffing resistance in the rolling-sliding contact machine element, was investigated. The scuffing test was performed using a two-cylinder rolling contact test machine. In a specific sliding, a faster roller of 60% and a sliding velocity of 1.75 m/s were utilized. The scuffing test results with shot-peened test rollers and those with non-shot-peened test roller were compared. The influence of the surface roughness of the shot-peened test roller was also discussed. We found that the shot-peened roller had a better scuffing resistance compared with the roller without the shot-peening process.
