The discrete element method(DEM) has been extensively adopted to investigate many complex geotechnical related problems due to its capability to incorporate the discontinuous nature of granular materials. In particula...The discrete element method(DEM) has been extensively adopted to investigate many complex geotechnical related problems due to its capability to incorporate the discontinuous nature of granular materials. In particular, when simulating large deformations or distortion of soil(e.g. cavity expansion),DEM can be very effective as other numerical solutions may experience convergence problems. Cavity expansion theory has widespread applications in geotechnical engineering, particularly to the problems concerning in situ testing, pile installation and so forth. In addition, the behaviour of geomaterials in a macro-level is utterly determined by microscopic properties, highlighting the importance of contact models. Despite the fact that there are numerous contact models proposed to mimic the realistic behaviour of granular materials, there are lack of studies on the effects of these contact models on the soil response.Hence, in this study, a series of three-dimensional numerical simulations with different contact constitutive models was conducted to simulate the response of sandy soils during cylindrical cavity expansion. In this numerical investigation, three contact models, i.e. linear contact model, rolling resistance contact model,and Hertz contact model, are considered. It should be noted that the former two models are linear based models, providing linearly elastic and frictional plasticity behaviours, whereas the latter one consists of nonlinear formulation based on an approximation of the theory of Mindlin and Deresiewicz. To examine the effects of these contact models, several cylindrical cavities were created and expanded gradually from an initial radius of 0.055 m to a final radius of 0.1 m. The numerical predictions confirm that the calibrated contact models produced similar results regarding the variations of cavity pressure, radial stress, deviatoric stress, volumetric strain, as well as the soil radial displacement. However, the linear contact model may result in inaccurate predictions when highly angular soil particles are involved. In addition, considering the excessive soil displacement induced by the pile installation(i.e. cavity expansion), a minimum distance of11 a(a is the cavity radius) is recommend for practicing engineers to avoid the potential damages to the existing piles and adjacent structures.展开更多
The models constructed by particle flow simulation method can effectively simulate the heterogeneous substance characteristics and failure behaviors of rocks.However,existing contact models overlook the rock cracks,an...The models constructed by particle flow simulation method can effectively simulate the heterogeneous substance characteristics and failure behaviors of rocks.However,existing contact models overlook the rock cracks,and the various simulation methods that do consider cracks still exhibit certain limitations.In this paper,based on Flat-Joint model and Linear Parallel Bond model,a crack contact model considering linked substance in the crack is proposed by splitting the crack contact into two portions:linked portion and unlinked portion for calculation.The new contact model considers the influence of crack closure on the contact force-displacement law.And a better compressive tensile strength ratio(UCS/T)was obtained by limiting the failure of the contact bond to be solely controlled by the contact force and moment of the linked portion.Then,by employing the FISH Model tool within the Particle Flow Code,the contact model was constructed and verified through contact force–displacement experiments and loading-unloading tests with cracked model.Finally,the contact model was tested through simulations of rock mechanics experiments.The results indicate that the contact model can effectively simulate the axial and lateral strain laws of rocks simultaneously and has a relatively good reproduction of the bi-modularity of rocks.展开更多
This study proposes a general imperfect thermal contact model to predict the thermal contact resistance at the interface among multi-layered composite structures.Based on the Green-Lindsay(GL)thermoelastic theory,semi...This study proposes a general imperfect thermal contact model to predict the thermal contact resistance at the interface among multi-layered composite structures.Based on the Green-Lindsay(GL)thermoelastic theory,semi analytical solutions of temperature increment and displacement of multi-layered composite structures are obtained by using the Laplace transform method,upon which the effects of thermal resistance coefficient,partition coefficient,thermal conductivity ratio and heat capacity ratio on the responses are studied.The results show that the generalized imperfect thermal contact model can realistically describe the imperfect thermal contact problem.Accordingly,it may degenerate into other thermal contact models by adjusting the thermal resistance coefficient and partition coefficient.展开更多
In this paper,an incremental contact model is developed for the elastic self-affine fractal rough surfaces under plane strain condition.The contact between a rough surface and a rigid plane is simplified by the accumu...In this paper,an incremental contact model is developed for the elastic self-affine fractal rough surfaces under plane strain condition.The contact between a rough surface and a rigid plane is simplified by the accumulation of identical line contacts with half-width given by the truncated area divided by the contact patch number at varying heights.Based on the contact stiffness of two-dimensional flat punch,the total stiffness of rough surface is estimated,and then the normal load is calculated by an incremental method.For various rough surfaces,the approximately linear load-area relationships predicted by the proposed model agree well with the results of finite element simulations.It is found that the real average contact pressure depends significantly on profile properties.展开更多
The complex geometrical features of mechanical components significantly influence contact interactions and system dynamics.However,directly modeling contact forces on surfaces with intricate geometries presents consid...The complex geometrical features of mechanical components significantly influence contact interactions and system dynamics.However,directly modeling contact forces on surfaces with intricate geometries presents considerable challenges.This study focuses on the helically twisted wire rope-sheave contact and proposes a contact force model that incorporates complex geometric features through a parameter identification approach.The model's impact on contact forces and system dynamics is thoroughly investigated.Leveraging a point contact model and an elliptic integral approximation,a loss function is formulated using the finite element(FE)contact model results as the reference data.Geometric parameters are subsequently determined by optimizing this loss function via a genetic algorithm(GA).The findings reveal that the contact stiffness increases with the wire rope pitch length,the radius of principal curvature,and the elliptic eccentricity of the contact zone.The proposed contact force model is integrated into a rigid-flexible coupled dynamics model,developed by the absolute node coordinate formulation,to examine the effects of contact geometry on system dynamics.The results demonstrate that the variations in wire rope geometry alter the contact stiffness,which in turn affects dynamic rope tension through frictional energy dissipation.The enhanced model's predictions exhibit superior alignment with the experimental data,thereby validating the methodology.This approach provides new insights for deducing the contact geometry from kinetic parameters and monitoring the performance degradation of mechanical components.展开更多
Aero engine seal coatings can effectively improve the air tightness of aircraft engines and increase fuel efficiency.However,due to the frictional forces between the blades and the coating,the coating often flakes off...Aero engine seal coatings can effectively improve the air tightness of aircraft engines and increase fuel efficiency.However,due to the frictional forces between the blades and the coating,the coating often flakes off,resulting in damage to the blades and causing eco-nomic losses.Therefore,it is necessary to analyze the friction between the blades and the coating.In this paper,three ceramic-based high-temperature seal coatings with different polyphenylene ester contents and a pure Yttria-stabilised zirconia coating were prepared by atmo-spheric plasma spraying(APS).The hardness and modulus of elasticity of the coated surfaces were obtained by hardness and modulus of elasticity tests,and the coatings were subjected to high-speed touch abrasion tests.The Hertzian contact model was used to calculate the maximum normal contact load on the coating during the process.The test values were compared with the theoretical values and it was found that the calculated values were always greater than the test values.In order to make the Hertzian contact model more accurate in calculating the touching and abrasion forces,the contact coefficients in the Hertzian contact model were optimized.Substituting the optimized coeffi-cients into the Hertzian contact model,the results show that the calculated results after optimizing the coefficients are much closer to the test values,with deviations from the test values ranging from 1%to 38%.展开更多
A comprehensive numerical investigation into mixed⁃mode delamination is presented in this study.It aims to assess the impact of thermal and friction effects through mixed⁃mode flexure crack propagation testing.Finite ...A comprehensive numerical investigation into mixed⁃mode delamination is presented in this study.It aims to assess the impact of thermal and friction effects through mixed⁃mode flexure crack propagation testing.Finite element analysis was employed to model the delamination process,incorporating a contact cohesive zone model.This model couples the traction⁃separation law,the contact law,and the Coulomb friction law simultaneously.The thermomechanical analysis in this study is performed using a sequentially coupled approach,implemented with the finite element software ABAQUS.The findings underscore the importance of this study.展开更多
A kind of improved contact frictional model on basis of traditional Coulomb Friction model is adopted.Corresponding contact element is also given.The contact algorithm on basis of augmented Lagrange method is introduc...A kind of improved contact frictional model on basis of traditional Coulomb Friction model is adopted.Corresponding contact element is also given.The contact algorithm on basis of augmented Lagrange method is introduced and successfully applied to complex contact friction problem.Test example and actual engineering case all show that the algorithm of the model is efficient and computation results agree well with general rules.展开更多
Because the result of the MB fractal model contradicts with the classical contact mechanics, a revised elastoplastic contact model of a single asperity is developed based on fractal theory. The critical areas of a sin...Because the result of the MB fractal model contradicts with the classical contact mechanics, a revised elastoplastic contact model of a single asperity is developed based on fractal theory. The critical areas of a single asperity are scale dependent, with an increase in the contact load and contact area, a transition from elastic, elastoplastic to full plastic deformation takes place in this order. In considering the size distribution function, analytic expression between the total contact load and the real contact area on the contact surface is obtained. The elastic, elastoplastic and full plastic contact load are obtained by the critical elastic contact area of the biggest asperity and maximun contact area of a single asperity. The results show that a rough surface is firstly in elastic deformation. As the load increases, elastoplastic or full plastic deformation takes place. For constant characteristic length scale G, the slope of load-area relation is proportional to fractal dimension D. For constant fractal dimension D, the slope of load-area relation is inversely proportional to G. For constant D and G, the slope of load-area relation is inversely proportional to property of the material ~b, namely with the same load, the material of rough surface is softer, and the total contact area is larger. The contact mechanics model provides a foundation for study of the friction, wear and seal performance of rough surfaces.展开更多
A contact model for describing the contact mechanics between the stator and slider of the standing wave linear ultrasonic motor was presented. The proposed model starts from the assumption that the vibration character...A contact model for describing the contact mechanics between the stator and slider of the standing wave linear ultrasonic motor was presented. The proposed model starts from the assumption that the vibration characteristics of the stator is not affected by the contact process. A modified friction models was used to analyze the contact problems. Firstly, the dynamic normal contact force, interface friction force, and steady-state characteristics were analyzed. Secondly, the influences of the contact layer material, the dynamic characteristics of the stator, and the pre-load on motor performance were simulated. Finally, to validate the contact model, a linear ultrasonic motor based on in-plane modes was used as an example. The corresponding results show that a set of simulation of motor performances based on the proposed contact mechanism is in good agreement with experimental results. This model is helpful to understanding the operation principle of the standing wave linear motor and thus contributes to the design of these types of motor.展开更多
A new elastic-plastic impact-contact model is proposed in this paper. By adopting the principle of minimum acceleration for elastic-plastic continue at finite deformation, and with the aid of finite difference method,...A new elastic-plastic impact-contact model is proposed in this paper. By adopting the principle of minimum acceleration for elastic-plastic continue at finite deformation, and with the aid of finite difference method, the proposed model is applied in the problem of dynamic response of a clamped thin circular plate subjected to a projectile impact centrally. The impact force history and response characteristics of the target plate is studied in detail. The theoretical predictions of the impact force and plate deflection are in good agreements with those of LDA experimental data. Linear expressions of the maximum impact force/transverse deflection versus impact velocity are given on the basis of the theoretical results.展开更多
The spindle barrel finishing is commonly used to improve the surface integrity of the important parts of the high-end equipment while it is difficult to provide enough test artifacts for the traditional trial and erro...The spindle barrel finishing is commonly used to improve the surface integrity of the important parts of the high-end equipment while it is difficult to provide enough test artifacts for the traditional trial and error experiment to obtain the desirable processing technology.The EDEM simulation of the spindle barrel finishing can provide effective help for the process design,however,the difference between the simulation and experiment is closely related to the selection of the contact model during simulation.In this paper,simulations and experiments are conducted based on the identical apparatus and conditions to facilitate the comparison and validation between each other.Based on the Hertz contact theory,the effect of the material properties of contact objects and the relative position of the workpiece on the contact force is qualified.The expression of the correlation coefficient of the contact model is deduced.Then the formula for calculating the contact force between the barrel finishing abrasive and the workpiece that includes influence coefficient of the material properties and the relative positions is established.Finally,the contact force calculation formula is verified by changing the rotating speed.The result shows that the material correction coefficient ranges from 1.41 to 2.38,which is inversely related to the equivalent modulus E.The position correction coefficient ranges from 2.0 to 2.3.The relative error value between the calculation result and the experimental test result is from 0.58%to 14.07%.This research lay a theoretical foundation for the correction theory of the core elements of the spindle barrel finishing process.展开更多
The contact sti ness of a mechanical bonding surface is an important parameter in determining the normal and radial contact force. To improve the calculation accuracy of the contact force model, the surface roughness ...The contact sti ness of a mechanical bonding surface is an important parameter in determining the normal and radial contact force. To improve the calculation accuracy of the contact force model, the surface roughness of the bonding surface and the energy loss that necessarily occurs during the impact process should be considered com?prehensively. To study the normal contact force of a revolute joint with clearance more accurately in the case of dry friction, a nonlinear sti ness coe cient model considering the surface roughness was established based on fractal theory, which considers the elastic, elastic?plastic, and plastic deformations of the asperities of the contact surface during the contact process. On this basis, a modified nonlinear spring damping model was established based on the Lankarani–Nikravesh contact force model. The laws influencing the surface roughness, recovery coe cient, initial velocity, and clearance size on the impact force were revealed, and were compared with the Lankarani–Nikravesh model and a hybrid model using MATLAB. The maximum impact force was obtained using a modified contact force model under di erent initial velocities, di erent clearances, and di erent degrees of surface roughness, and the calculated results were then compared with the experiment results. This study indicates that the modified model can be used more widely than other models, and is suitable for both large and small clearances. In particular, the modified model is more accurate when calculating the contact force of a revolute joint with a small clearance.展开更多
The development of a general discrete element method for irregularly shaped particles is the core issue of the simulation of the dynamic behavior of granular materials.The general energy-conserving contact theory is u...The development of a general discrete element method for irregularly shaped particles is the core issue of the simulation of the dynamic behavior of granular materials.The general energy-conserving contact theory is used to establish a universal discrete element method suitable for particle contact of arbitrary shape.In this study,three dimentional(3D)modeling and scanning techniques are used to obtain a triangular mesh representation of the true particles containing typical concave particles.The contact volumebased energy-conserving model is used to realize the contact detection between irregularly shaped particles,and the contact force model is refined and modified to describe the contact under real conditions.The inelastic collision processes between the particles and boundaries are simulated to verify the robustness of the modified contact force model and its applicability to the multi-point contact mode.In addition,the packing process and the flow process of a large number of irregular particles are simulated with the modified discrete element method(DEM)to illustrate the applicability of the method of complex problems.展开更多
A new contact model of traveling wave ultrasonic motor （TWUSM） with a visco-elastic stator frictional layer was presented. In this model, the initial boundaries were revised, and the rotor revolution speed could he ...A new contact model of traveling wave ultrasonic motor （TWUSM） with a visco-elastic stator frictional layer was presented. In this model, the initial boundaries were revised, and the rotor revolution speed could he calculated iteratively. This model was compared with compliant slider and rigid stator model. The results of motor characteristics simulations showed that the motors based on this model would gain bigger stall torque. Then the friction and wear characteristics of two models were analyzed. The motors based on this model had lower coefficient of friction and better wear resistance.展开更多
The distinct element method（DEM） incorporated with a novel bond contact model was applied in this paper to shed light on the microscopic physical origin of macroscopic behaviors of weathered rock, and to achieve the...The distinct element method（DEM） incorporated with a novel bond contact model was applied in this paper to shed light on the microscopic physical origin of macroscopic behaviors of weathered rock, and to achieve the changing laws of microscopic parameters from observed decaying properties of rocks during weathering. The changing laws of macroscopic mechanical properties of typical rocks were summarized based on the existing research achievements. Parametric simulations were then conducted to analyze the relationships between macroscopic and microscopic parameters, and to derive the changing laws of microscopic parameters for the DEM model. Equipped with the microscopic weathering laws, a series of DEM simulations of basic laboratory tests on weathered rock samples was performed in comparison with analytical solutions. The results reveal that the relationships between macroscopic and microscopic parameters of rocks against the weathering period can be successfully attained by parametric simulations. In addition, weathering has a significant impact on both stressestrain relationship and failure pattern of rocks.展开更多
The forming performance of sheet metals in the deep-drawing process with ultrasonic vibrations can be improved by the surface effect between the sheet metal and the die.A sheet metal friction test with ultrasonic vibr...The forming performance of sheet metals in the deep-drawing process with ultrasonic vibrations can be improved by the surface effect between the sheet metal and the die.A sheet metal friction test with ultrasonic vibrations is performed to explore the cause of the surface effect.The frictional characteristics are investigated,and the corresponding friction expressions are established based on the contact mechanics and the elastic–plastic contact model for rough surfaces.Friction is caused by the elastic–plastic deformation of contacting asperities under normal loads.The actual contacting region between two surfaces increases with normal loads,whereas the normal distance decreases.The normal distance between the contacting surfaces is changed,asperities generate a tangential deformation with ultrasonic vibrations,and the friction coefficient is eventually altered.Ultrasonic vibrations are applied on a 40Cr steel punch at the frequency of 20 kHz and the amplitude of 4.2μm.In the friction tests,the punch is perpendicular to the surface of the magnesium alloy AZ31B sheet metals and is sliding with a relative velocity of 1 mm/s.The test results show that the friction coefficient is decreased by approximately 40%and the theoretical values are in accordance with the test values;Ultrasonic vibrations can clearly reduce wear and improve the surface quality of parts.展开更多
Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft...Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.展开更多
The existence of clearance in the joints of mechanisms system is inevitable.The movements of the real mechanism are deflection from the ideal mechanism due to the clearances and the motion accuracy is decreased.The ef...The existence of clearance in the joints of mechanisms system is inevitable.The movements of the real mechanism are deflection from the ideal mechanism due to the clearances and the motion accuracy is decreased.The effects of the hinge clearance on the crank and rocker mechanism system are studied.The system dynamics equation with clearance is presented.The contact dynamics model is established using the nonlinear equivalent spring-damp model and the friction effect is considered by using Coulomb friction model.Then the models are incorporated into ADAMS,and based on the model,large numbers numeric simulations are made.The regularity of contact forces in clearance are studied in detail.And the effects of clearance size,clearance friction on the mechanism dynamics characteristic are analyzed.The simulation results can predict the effects of clearance on the mechanism dynamics characteristic preferably.展开更多
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.展开更多
文摘The discrete element method(DEM) has been extensively adopted to investigate many complex geotechnical related problems due to its capability to incorporate the discontinuous nature of granular materials. In particular, when simulating large deformations or distortion of soil(e.g. cavity expansion),DEM can be very effective as other numerical solutions may experience convergence problems. Cavity expansion theory has widespread applications in geotechnical engineering, particularly to the problems concerning in situ testing, pile installation and so forth. In addition, the behaviour of geomaterials in a macro-level is utterly determined by microscopic properties, highlighting the importance of contact models. Despite the fact that there are numerous contact models proposed to mimic the realistic behaviour of granular materials, there are lack of studies on the effects of these contact models on the soil response.Hence, in this study, a series of three-dimensional numerical simulations with different contact constitutive models was conducted to simulate the response of sandy soils during cylindrical cavity expansion. In this numerical investigation, three contact models, i.e. linear contact model, rolling resistance contact model,and Hertz contact model, are considered. It should be noted that the former two models are linear based models, providing linearly elastic and frictional plasticity behaviours, whereas the latter one consists of nonlinear formulation based on an approximation of the theory of Mindlin and Deresiewicz. To examine the effects of these contact models, several cylindrical cavities were created and expanded gradually from an initial radius of 0.055 m to a final radius of 0.1 m. The numerical predictions confirm that the calibrated contact models produced similar results regarding the variations of cavity pressure, radial stress, deviatoric stress, volumetric strain, as well as the soil radial displacement. However, the linear contact model may result in inaccurate predictions when highly angular soil particles are involved. In addition, considering the excessive soil displacement induced by the pile installation(i.e. cavity expansion), a minimum distance of11 a(a is the cavity radius) is recommend for practicing engineers to avoid the potential damages to the existing piles and adjacent structures.
基金supported by the Natural Science Foundation of Heilongjiang Province(No.ZD2021E006)the National Natural Science Foundation of China(Nos.52174075 and 52074110).
文摘The models constructed by particle flow simulation method can effectively simulate the heterogeneous substance characteristics and failure behaviors of rocks.However,existing contact models overlook the rock cracks,and the various simulation methods that do consider cracks still exhibit certain limitations.In this paper,based on Flat-Joint model and Linear Parallel Bond model,a crack contact model considering linked substance in the crack is proposed by splitting the crack contact into two portions:linked portion and unlinked portion for calculation.The new contact model considers the influence of crack closure on the contact force-displacement law.And a better compressive tensile strength ratio(UCS/T)was obtained by limiting the failure of the contact bond to be solely controlled by the contact force and moment of the linked portion.Then,by employing the FISH Model tool within the Particle Flow Code,the contact model was constructed and verified through contact force–displacement experiments and loading-unloading tests with cracked model.Finally,the contact model was tested through simulations of rock mechanics experiments.The results indicate that the contact model can effectively simulate the axial and lateral strain laws of rocks simultaneously and has a relatively good reproduction of the bi-modularity of rocks.
基金Projects(42477162,52108347,52178371,52168046,52178321,52308383)supported by the National Natural Science Foundation of ChinaProjects(2023C03143,2022C01099,2024C01219,2022C03151)supported by the Zhejiang Key Research and Development Plan,China+6 种基金Project(LQ22E080010)supported by the Exploring Youth Project of Zhejiang Natural Science Foundation,ChinaProject(LR21E080005)supported by the Outstanding Youth Project of Natural Science Foundation of Zhejiang Province,ChinaProject(2022M712964)supported by the Postdoctoral Science Foundation of ChinaProject(2023AFB008)supported by the Natural Science Foundation of Hubei Province for Youth,ChinaProject(202203)supported by Engineering Research Centre of Rock-Soil Drilling&Excavation and Protection,Ministry of Education,ChinaProject(202305-2)supported by the Science and Technology Project of Zhejiang Provincial Communication Department,ChinaProject(2021K256)supported by the Construction Research Founds of Department of Housing and Urban-Rural Development of Zhejiang Province,China。
文摘This study proposes a general imperfect thermal contact model to predict the thermal contact resistance at the interface among multi-layered composite structures.Based on the Green-Lindsay(GL)thermoelastic theory,semi analytical solutions of temperature increment and displacement of multi-layered composite structures are obtained by using the Laplace transform method,upon which the effects of thermal resistance coefficient,partition coefficient,thermal conductivity ratio and heat capacity ratio on the responses are studied.The results show that the generalized imperfect thermal contact model can realistically describe the imperfect thermal contact problem.Accordingly,it may degenerate into other thermal contact models by adjusting the thermal resistance coefficient and partition coefficient.
基金supported by the National Natural Science Foundation of China(Grant Nos.12372100,12302126,and 12302141)the China Postdoctoral Science Foundation(Grant No.2023M732799)+1 种基金the Fundamental Research Funds for the Central Universities(Grant No.xzy012024020)Sihe Wang also thanks the support from the China Scholarship Council(CSC).
文摘In this paper,an incremental contact model is developed for the elastic self-affine fractal rough surfaces under plane strain condition.The contact between a rough surface and a rigid plane is simplified by the accumulation of identical line contacts with half-width given by the truncated area divided by the contact patch number at varying heights.Based on the contact stiffness of two-dimensional flat punch,the total stiffness of rough surface is estimated,and then the normal load is calculated by an incremental method.For various rough surfaces,the approximately linear load-area relationships predicted by the proposed model agree well with the results of finite element simulations.It is found that the real average contact pressure depends significantly on profile properties.
基金supported by the National Key Research and Development Program of China(No.2023YFC3010400)。
文摘The complex geometrical features of mechanical components significantly influence contact interactions and system dynamics.However,directly modeling contact forces on surfaces with intricate geometries presents considerable challenges.This study focuses on the helically twisted wire rope-sheave contact and proposes a contact force model that incorporates complex geometric features through a parameter identification approach.The model's impact on contact forces and system dynamics is thoroughly investigated.Leveraging a point contact model and an elliptic integral approximation,a loss function is formulated using the finite element(FE)contact model results as the reference data.Geometric parameters are subsequently determined by optimizing this loss function via a genetic algorithm(GA).The findings reveal that the contact stiffness increases with the wire rope pitch length,the radius of principal curvature,and the elliptic eccentricity of the contact zone.The proposed contact force model is integrated into a rigid-flexible coupled dynamics model,developed by the absolute node coordinate formulation,to examine the effects of contact geometry on system dynamics.The results demonstrate that the variations in wire rope geometry alter the contact stiffness,which in turn affects dynamic rope tension through frictional energy dissipation.The enhanced model's predictions exhibit superior alignment with the experimental data,thereby validating the methodology.This approach provides new insights for deducing the contact geometry from kinetic parameters and monitoring the performance degradation of mechanical components.
基金supported by Basic Research Funds for Central Universities(3122019189).
文摘Aero engine seal coatings can effectively improve the air tightness of aircraft engines and increase fuel efficiency.However,due to the frictional forces between the blades and the coating,the coating often flakes off,resulting in damage to the blades and causing eco-nomic losses.Therefore,it is necessary to analyze the friction between the blades and the coating.In this paper,three ceramic-based high-temperature seal coatings with different polyphenylene ester contents and a pure Yttria-stabilised zirconia coating were prepared by atmo-spheric plasma spraying(APS).The hardness and modulus of elasticity of the coated surfaces were obtained by hardness and modulus of elasticity tests,and the coatings were subjected to high-speed touch abrasion tests.The Hertzian contact model was used to calculate the maximum normal contact load on the coating during the process.The test values were compared with the theoretical values and it was found that the calculated values were always greater than the test values.In order to make the Hertzian contact model more accurate in calculating the touching and abrasion forces,the contact coefficients in the Hertzian contact model were optimized.Substituting the optimized coeffi-cients into the Hertzian contact model,the results show that the calculated results after optimizing the coefficients are much closer to the test values,with deviations from the test values ranging from 1%to 38%.
文摘A comprehensive numerical investigation into mixed⁃mode delamination is presented in this study.It aims to assess the impact of thermal and friction effects through mixed⁃mode flexure crack propagation testing.Finite element analysis was employed to model the delamination process,incorporating a contact cohesive zone model.This model couples the traction⁃separation law,the contact law,and the Coulomb friction law simultaneously.The thermomechanical analysis in this study is performed using a sequentially coupled approach,implemented with the finite element software ABAQUS.The findings underscore the importance of this study.
文摘A kind of improved contact frictional model on basis of traditional Coulomb Friction model is adopted.Corresponding contact element is also given.The contact algorithm on basis of augmented Lagrange method is introduced and successfully applied to complex contact friction problem.Test example and actual engineering case all show that the algorithm of the model is efficient and computation results agree well with general rules.
基金Supported by National Natural Science Foundation of China(Grant Nos.51105304,51475364)Shaanxi Provincial Natural Science Basic Research Plan of China(Grant No.2015JM5212)
文摘Because the result of the MB fractal model contradicts with the classical contact mechanics, a revised elastoplastic contact model of a single asperity is developed based on fractal theory. The critical areas of a single asperity are scale dependent, with an increase in the contact load and contact area, a transition from elastic, elastoplastic to full plastic deformation takes place in this order. In considering the size distribution function, analytic expression between the total contact load and the real contact area on the contact surface is obtained. The elastic, elastoplastic and full plastic contact load are obtained by the critical elastic contact area of the biggest asperity and maximun contact area of a single asperity. The results show that a rough surface is firstly in elastic deformation. As the load increases, elastoplastic or full plastic deformation takes place. For constant characteristic length scale G, the slope of load-area relation is proportional to fractal dimension D. For constant fractal dimension D, the slope of load-area relation is inversely proportional to G. For constant D and G, the slope of load-area relation is inversely proportional to property of the material ~b, namely with the same load, the material of rough surface is softer, and the total contact area is larger. The contact mechanics model provides a foundation for study of the friction, wear and seal performance of rough surfaces.
基金Funded by the National Basic Research Program (973 program) (No. 2011CB707602)the Digital Manufacturing Equipment and Technology National Key Laboratory,Huazhong University of Science and Technology (No. DMETKF2009002)National Sciences Foundation-Guangdong Natural Science Foundation,China (No.U0934004)
文摘A contact model for describing the contact mechanics between the stator and slider of the standing wave linear ultrasonic motor was presented. The proposed model starts from the assumption that the vibration characteristics of the stator is not affected by the contact process. A modified friction models was used to analyze the contact problems. Firstly, the dynamic normal contact force, interface friction force, and steady-state characteristics were analyzed. Secondly, the influences of the contact layer material, the dynamic characteristics of the stator, and the pre-load on motor performance were simulated. Finally, to validate the contact model, a linear ultrasonic motor based on in-plane modes was used as an example. The corresponding results show that a set of simulation of motor performances based on the proposed contact mechanism is in good agreement with experimental results. This model is helpful to understanding the operation principle of the standing wave linear motor and thus contributes to the design of these types of motor.
基金The project supported by the National Natural Science Foundation of China(10532020)
文摘A new elastic-plastic impact-contact model is proposed in this paper. By adopting the principle of minimum acceleration for elastic-plastic continue at finite deformation, and with the aid of finite difference method, the proposed model is applied in the problem of dynamic response of a clamped thin circular plate subjected to a projectile impact centrally. The impact force history and response characteristics of the target plate is studied in detail. The theoretical predictions of the impact force and plate deflection are in good agreements with those of LDA experimental data. Linear expressions of the maximum impact force/transverse deflection versus impact velocity are given on the basis of the theoretical results.
基金Supported by Program National Natural Science Foundation of China(Grant Nos.51875389,51975399,52075362)Key Program of Natural Science Foundation of Shanxi Province of China(Grant No.201801D111002)Scientific and Technological Innovation Project for Excellent Talents in Shanxi Province of China(Grant No.201805D211031).
文摘The spindle barrel finishing is commonly used to improve the surface integrity of the important parts of the high-end equipment while it is difficult to provide enough test artifacts for the traditional trial and error experiment to obtain the desirable processing technology.The EDEM simulation of the spindle barrel finishing can provide effective help for the process design,however,the difference between the simulation and experiment is closely related to the selection of the contact model during simulation.In this paper,simulations and experiments are conducted based on the identical apparatus and conditions to facilitate the comparison and validation between each other.Based on the Hertz contact theory,the effect of the material properties of contact objects and the relative position of the workpiece on the contact force is qualified.The expression of the correlation coefficient of the contact model is deduced.Then the formula for calculating the contact force between the barrel finishing abrasive and the workpiece that includes influence coefficient of the material properties and the relative positions is established.Finally,the contact force calculation formula is verified by changing the rotating speed.The result shows that the material correction coefficient ranges from 1.41 to 2.38,which is inversely related to the equivalent modulus E.The position correction coefficient ranges from 2.0 to 2.3.The relative error value between the calculation result and the experimental test result is from 0.58%to 14.07%.This research lay a theoretical foundation for the correction theory of the core elements of the spindle barrel finishing process.
基金Supported by National Natural Science Foundation of China(Grant No.51775475)Hebei Provincial Natural Science Foundation of China(Grant No.E2016203463)
文摘The contact sti ness of a mechanical bonding surface is an important parameter in determining the normal and radial contact force. To improve the calculation accuracy of the contact force model, the surface roughness of the bonding surface and the energy loss that necessarily occurs during the impact process should be considered com?prehensively. To study the normal contact force of a revolute joint with clearance more accurately in the case of dry friction, a nonlinear sti ness coe cient model considering the surface roughness was established based on fractal theory, which considers the elastic, elastic?plastic, and plastic deformations of the asperities of the contact surface during the contact process. On this basis, a modified nonlinear spring damping model was established based on the Lankarani–Nikravesh contact force model. The laws influencing the surface roughness, recovery coe cient, initial velocity, and clearance size on the impact force were revealed, and were compared with the Lankarani–Nikravesh model and a hybrid model using MATLAB. The maximum impact force was obtained using a modified contact force model under di erent initial velocities, di erent clearances, and di erent degrees of surface roughness, and the calculated results were then compared with the experiment results. This study indicates that the modified model can be used more widely than other models, and is suitable for both large and small clearances. In particular, the modified model is more accurate when calculating the contact force of a revolute joint with a small clearance.
基金supported by the National Key Research and Development Program of China(2018YFA0605902)the National Natural Science Foundation of China(42176241 and 11872136).
文摘The development of a general discrete element method for irregularly shaped particles is the core issue of the simulation of the dynamic behavior of granular materials.The general energy-conserving contact theory is used to establish a universal discrete element method suitable for particle contact of arbitrary shape.In this study,three dimentional(3D)modeling and scanning techniques are used to obtain a triangular mesh representation of the true particles containing typical concave particles.The contact volumebased energy-conserving model is used to realize the contact detection between irregularly shaped particles,and the contact force model is refined and modified to describe the contact under real conditions.The inelastic collision processes between the particles and boundaries are simulated to verify the robustness of the modified contact force model and its applicability to the multi-point contact mode.In addition,the packing process and the flow process of a large number of irregular particles are simulated with the modified discrete element method(DEM)to illustrate the applicability of the method of complex problems.
基金Sponsored by the Natural Sciences Foundation of China (Grant No 50575055)the National 973 Key Development Project in Fundamental Science Research (Grant No2007CB607602)
文摘A new contact model of traveling wave ultrasonic motor （TWUSM） with a visco-elastic stator frictional layer was presented. In this model, the initial boundaries were revised, and the rotor revolution speed could he calculated iteratively. This model was compared with compliant slider and rigid stator model. The results of motor characteristics simulations showed that the motors based on this model would gain bigger stall torque. Then the friction and wear characteristics of two models were analyzed. The motors based on this model had lower coefficient of friction and better wear resistance.
基金funded by the National Basic Research Programs of China(Grant Nos.2011CB013504 and 2014CB046901)the National Funds for Distinguished Young Scientists of China(Grant No.51025932)the National Nature Science Foundation of China(Grant No.41372272)
文摘The distinct element method（DEM） incorporated with a novel bond contact model was applied in this paper to shed light on the microscopic physical origin of macroscopic behaviors of weathered rock, and to achieve the changing laws of microscopic parameters from observed decaying properties of rocks during weathering. The changing laws of macroscopic mechanical properties of typical rocks were summarized based on the existing research achievements. Parametric simulations were then conducted to analyze the relationships between macroscopic and microscopic parameters, and to derive the changing laws of microscopic parameters for the DEM model. Equipped with the microscopic weathering laws, a series of DEM simulations of basic laboratory tests on weathered rock samples was performed in comparison with analytical solutions. The results reveal that the relationships between macroscopic and microscopic parameters of rocks against the weathering period can be successfully attained by parametric simulations. In addition, weathering has a significant impact on both stressestrain relationship and failure pattern of rocks.
基金Projects(51775480,51305385)supported by the National Natural Science Foundation of ChinaProject(E2018203143)supported by the Natural Science Foundation of Hebei Province,China
文摘The forming performance of sheet metals in the deep-drawing process with ultrasonic vibrations can be improved by the surface effect between the sheet metal and the die.A sheet metal friction test with ultrasonic vibrations is performed to explore the cause of the surface effect.The frictional characteristics are investigated,and the corresponding friction expressions are established based on the contact mechanics and the elastic–plastic contact model for rough surfaces.Friction is caused by the elastic–plastic deformation of contacting asperities under normal loads.The actual contacting region between two surfaces increases with normal loads,whereas the normal distance decreases.The normal distance between the contacting surfaces is changed,asperities generate a tangential deformation with ultrasonic vibrations,and the friction coefficient is eventually altered.Ultrasonic vibrations are applied on a 40Cr steel punch at the frequency of 20 kHz and the amplitude of 4.2μm.In the friction tests,the punch is perpendicular to the surface of the magnesium alloy AZ31B sheet metals and is sliding with a relative velocity of 1 mm/s.The test results show that the friction coefficient is decreased by approximately 40%and the theoretical values are in accordance with the test values;Ultrasonic vibrations can clearly reduce wear and improve the surface quality of parts.
基金Special Fund of Marine Commonweal Industry under contact Nos 201105016 and 201205007supported by National Marine Environment Forecasting Centrethe National Natural Science Foundation of China under contact No.41176012
文摘Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.
基金Sponsored by Program for Changjiang Scholars and Innovative Research Team in University(Grant No.IRT0520)
文摘The existence of clearance in the joints of mechanisms system is inevitable.The movements of the real mechanism are deflection from the ideal mechanism due to the clearances and the motion accuracy is decreased.The effects of the hinge clearance on the crank and rocker mechanism system are studied.The system dynamics equation with clearance is presented.The contact dynamics model is established using the nonlinear equivalent spring-damp model and the friction effect is considered by using Coulomb friction model.Then the models are incorporated into ADAMS,and based on the model,large numbers numeric simulations are made.The regularity of contact forces in clearance are studied in detail.And the effects of clearance size,clearance friction on the mechanism dynamics characteristic are analyzed.The simulation results can predict the effects of clearance on the mechanism dynamics characteristic preferably.
基金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.
