Contact Model Research Articles

Establishment of Whole-Rice-Plant Model and Calibration of Characteristic Parameters Based on Segmented Hollow Stalks

To address the limitations of the current discrete element model of rice plants in terms of accurately reflecting structural differences and threshing characteristics, this study proposes a whole-rice-plant modeling method based on segmented hollow stalks and establishes a whole-rice-plant model that accurately represents the bending and threshing characteristics of the actual rice plant. Initially, based on the characteristics of the rice plant, the rice stalk was segmented into three sections of hollow stalks with distinct structures—namely, the primary stalk, the secondary stalk, and the tertiary stalk—ensuring that the model closely resembles actual rice plants. Secondly, the mechanical and contact parameters for each structure of the rice plant were measured and calibrated through mechanical and contact tests. Subsequently, utilizing the Hertz–Mindlin contact model, a multi-dimensional element particle arrangement method was employed to establish a discrete element model of the entire rice plant. The bending characteristics of the stalk and the threshing characteristics of the rice were calibrated using three-point bending tests and impact threshing tests. The results indicated calibration errors in the bending resistance force of 4.46%, 3.95%, and 2.51% for the three-section stalk model, and the calibration error for the rice model’s threshing rate was 1.86%, which can accurately simulate the bending characteristics of the stalk and the threshing characteristics of the rice plant. Finally, the contact characteristics of the model were validated through a stack angle verification test, which revealed that the relative error of the stacking angle did not exceed 7.52%, confirming the accuracy of the contact characteristics of the rice plant model. The findings of this study provide foundational models and a theoretical basis for the simulation of and analytical applications related to rice threshing and cleaning.

Analysis of Short-range Contact Forces in Peridynamics Endowed with an Improved Nonlocal Contact Model

The motion and deformation laws of multi object interactions in computational models depend on contact algorithms. However, research on the peridynamic contact problem is limited. In this paper, in order to effectively prevent non-physical intrusion during contact, the inherent problems of the contact algorithm in peridynamic discrete models are analyzed, and a force boundary contact method using nonlinear contact stiffness is proposed. Through numerical calculations and geometric analysis of the peridynamic discrete model, the characteristics and reasons for the variation of contact force in the peridynamic contact model under fixed contact stiffness are discovered. To address the nonlinear reduction of contact force and potential non-physical intrusion in the peridynamic contact model, a force boundary peridynamic contact method is proposed by introducing a nonlinear changing contact stiffness function. Then the characteristics of contact force variation in different kinds of contact functions are studied and the parameter setting of contact functions is discussed. The results show that this method can effectively prevent non-physical intrusion and the calculation error is acceptable. This paper lays a foundation for further research on contact constitutive model based on Peridynamic framework.

A novel elastoplastic impact contact model for thin orthotropic layer

A complex stress state is often an obstacle in obtaining analytical solutions to elastoplastic contact problems of orthotropic structural materials. In this study, an analytical model is presented for investigating the impact contact between a rigid body and a thin orthotropic layer situated on a rigid foundation. By assuming that the local indentation during impact contact is due to the elastoplastic deformation, a theoretical study is carried out to predict the contact response of the orthotropic layer, which obeys an elastic-perfectly plastic stress-strain law. A relationship between contact force and indentation is derived, and the coefficient governing the rebound response is determined. The presented results show generally good agreement with the experimental and numerical results available in the literature. The impact contact model can also be utilized in the impact response analysis of coated structures. Parametric analysis results indicate that the elastic model tends to overestimate the impact resistance of thin layers. The elastoplastic contact law can accurately account for the decrease in contact force due to plastic indentation and permanent deformation. Moreover, the yield strength significantly influences the impact contact time and the permanent indentation deformation of the thin-layer structure.

Dynamic wheel–rail adhesion characteristics on wet curved tracks considering surface roughness

Wheel–rail rolling contact is a complex tribological issue, especially the existence of a thin water film on the rail surface, which leads to low adhesion problems for the traction and braking performance of the railway trains. The objective of this paper is to propose a non-Hertzian dynamic wheel–rail adhesion model for trains passing through wet curved tracks considering wheel–rail surface roughness. A vehicle-track coupling dynamics model was first developed to output the wheel–rail dynamic contact parameters under wet conditions. These parameters were subsequently utilized as dynamic inputs for the present wheel–rail adhesion model to simulate dynamic curve passing. The present numerical model includes the normal contact model and tangential model. For the normal contact issue, an elastohydrodynamic lubrication model was used to accurately determine the normal contact stress within the contact patch along the contact traces. Based on the normal contact characteristics, the advanced extended creep force model was used to compute the tangential stress. Finally, the dynamic wheel–rail adhesion coefficient can be obtained. This model considers the effects of time variation and changes in contact geometry, resulting in sudden increases or decreases in the adhesion coefficient at certain locations. At a speed of 400 km/h, a higher spin creepage occurs, leading to an increased displacement of the third-body layer and causing the maximum tangential stress to shift from one side of the contact patch to the other. Under the combined influence of the creepages, the maximum tangential stress occurs at the trailing edge of the contact patch.

Physics-based numerical implementation framework towardsmuti-scale contact problem

This paper establishes a general computational framework to solve the muti-scale contact problem by integrating the statistical contact model with the finite element format. Compared to existing models, the proposed method is applicable to most geometric configurations and can effectively evaluate the pressure distribution. In this work, a modified Karush-Kuhn-Tucker (KKT) condition is proposed by the assumption that asperity height obeys the Gaussian distribution. Therefore, in the variational formula, the contact contribution is decomposed into body contribution and asperity contribution, corresponding to the nominal smooth surface and roughness, respectively. Then the linearization and constraint enforcement of these two components are derived, followed by a nonlinear Newton-Raphson-based iterative algorithm. The contact patch test and Hertz contact test are conducted, and the predicted results are consistent with the theoretical and experimental values, confirming the effectiveness and accuracy of the proposed approach. It is worth noting that in the Hertz contact test, the contact pressure distribution varies progressively with the roughness level and external force, tending to the Hertz limit or Gaussian limit. This means that the proposed method can be applied to any roughness and load. Finally, the contact behaviors of the transmission interface of a piezoelectric actuator, i.e., a typical multi-scale contact problem, are studied as an engineering application case.

Research on an equivalent calculation method of the sealing structure of subsea pipeline connectors based on dimensionality reduction method

For the sealing structure of subsea pipeline connectors, it is crucial to understand the maximum equivalent stress of the sealing ring, the normal force generated by the pressure effect of the internal fluid acting on the contact area of the flanges and the sealing ring, and the displacement of the flange in the pre-tightening state and the working state. This paper proposes an equivalent calculation method (ECM) to quickly compute the values of these three variables. A theoretical model for the parallel contact of the axially symmetric object and an elastic foundation is proposed based on the dimensionality reduction method and Hertzian theory. Considering the contact and structural characteristics of the flanges and the sealing ring of the sealing structure, an equivalent hollow cylinder model is introduced. By combining these two models, along with the long hollow cylinder theory, theoretical formulas for the maximum equivalent stress, the normal force, and the flange displacement in the pre-tightening state and the working state of the sealing structure are derived. In addition, finite element simulations were conducted on the sealing structure of different sizes and working states. A comparative analysis was performed between the simulation results and the theoretical formulas results. The results show that, within the yield strength of the sealing ring, the maximum relative deviation of the three theoretical formulas is less than 10%. The theoretical formulas can be used for the design and engineering application of the sealing structure of subsea pipeline connectors.

The Impact of the Yeoh Model’s Variability in Contact on Knee Joint Mechanics

The aim of this study was to assess the impact of the variability of the Yeoh model when modeling the contact of bones through cartilage in the knee in compression and flexion–extension within a hybrid knee model. Firstly, a Sobol sequence of 64 samples and four variables representing the Yeoh parameters of the cartilage of the femur and tibia was generated. Based on these samples, 2 × 64 finite element contact models of the geometry of the sphere plane were generated and solved for healthy tissue affected by osteoarthritis. The resulting indentation curves were incorporated into a multibody knee joint model. The obtained results suggested that cartilage variability severely affected the knee in compression by up to 32%. However, the same variability also affected the flexion–extension motion, although to a lesser extent, with a relative change to the range of angular displacements of almost 7%. Osteoarthritic tissue was consistently more affected by this variability, suggesting that when modeling degenerated tissue, complex joint models are necessary.

Fiber random contact model for evaluating electrical resistance variation of carbon fiber composites under various damage modes

AbstractThe variations in electrical resistance of carbon fiber composites under different damage modes are critical for the design of electrically conductive components. Here, we developed a fiber random contact model to analyze the changes in electrical resistance of carbon fiber composites subjected to various damage scenarios. An ideal rectangular grid was constructed in matrix form utilizing the nodal voltage and current method. The electrical resistance of each fiber branch was calculated based on the normal distribution of fiber orientation angles. To simulate the random distribution of contact branches among fiber contact points within the resistor network, contact branches were randomly removed from the ideal rectangular grid. The electric conductive properties of carbon fiber tows, unidirectional carbon fiber composites, and 3D woven carbon fiber composites were measured to validate the fiber random contact model. The findings are anticipated to contribute to the optimization of carbon fiber composites for electrical applications.Highlights Fiber random lap model for evaluating electrical behavior was developed. The conductive behaviors of the composites were measured. Electrical resistance variation under various damages was investigated.

Fretting wear modeling of bolted joint interface with microscopic roughness using the 1D microslip friction model and equivalent thin layer

In this work, an improved 1D microslip model at the frictional joint interface considering the effect of surface topography is developed, which is further used to evaluate the fretting wear performance of the joint structure. The rough surface contact of the joint interface is characterized by the micro-asperity Greenwood and Williamson statistical contact model (GW model). A thin layer with elastic tangential stiffness, exhibiting equivalent properties to accord with the micro-asperity contact of the rough interface, is supposed to connect with the two components comprising the joint interface. The nonuniform clamping pressure distribution is considered and the resulted stick-slip transitions along the interface are derived using the continuum elastic beam model. The tangential stiffness of micro-asperities for rough surface contact is derived and incorporated into the continuum beam model to obtain the relative slip displacement along the interface. Fretting wear depth resulted from the cyclic microslip at the interface is evaluated based on the relative sliding distance and the Archard wear model. Results show that the minimum tangential force causing microslip motion increases with surface roughness and preload. The slip length increases with the tangential force whilst decreases with the increased surface roughness. The wear depth progressively increases from stuck region to slip region, and decreases as the surface roughness and preload increase. Experimental validation is performed and excellent agreement is observed between model predictions and experimental results. The developed model provides important new insights into the prediction of fretting wear and loosening behavior of the bolted joints.

Particulate behaviour of soft granular materials: a case study on lentils

In this article, lentils are used as a case study to characterise the particulate behaviour of soft granular materials. Experiments were carried out on a single lentil particle or a pair of lentils. The single particle or a pair of particles in contact were compressed vertically to crushing or to a fixed vertical load. Then, in inter-particle tests, pairs of particles were slid over each other at a constant vertical load, and the tangential stiffness and coefficient of friction were estimated. The pairs of lentil particles in contact were also subjected to repeated normal and tangential loading. The presence or absence of the cover of lentil particles (shell) was found to affect their behaviour significantly under these loading conditions. The lentil particles have a very compliant shell and stiff core in normal loading; the stiffness of the shell is constant, and the core follows Hertz contact law. The lentil particles show less variability in their crushing strength, with high Weibull modulus (~ 10), in comparison to other natural granular materials like sand. With repeated cycles of vertical loading, the contact between a lentil pairs of particles becomes more stiff and less damp. In tangential loading, the coefficient of friction between lentil particles decreases with normal force while the contact stiffness increases. Further, in cyclic tangential loading, the coefficient of friction decreases and the contact stiffness increases with cycles. A simple contact model is also proposed to use in discrete element simulations.Graphical abstract

Non-Linear Increase of the Real Contact Area of PMMA Blocks and the Related Contact Model

Many experiments have supported the contact models, such as the GW and MB models, but the majority of previous validations have been performed under light loads, resulting in a linear relationship between normal force and contact area. However, the real contact area fraction should never equal one; there must be a limit smaller than the apparent area, implying that the real contact area cannot increase linearly indefinitely. In this paper, the real contact area between two polymethylmethacrylate (PMMA) blocks under heavy load is measured using the total reflection method, and the contact area is analyzed using the image processing method. The results show that the real contact area increases with normal load linearly in light loads but non-linearly in heavy loads; the number of contact spots increases with load linearly in light loads but also non-linearly in heavy loads, synchronous with the change in the real contact area. The GW, MB, and Zhao, Maietta, and Chang (ZMC) models were used to predict the experiment results, but none of them predicted the non-linear stage. A revised GW model based on the bulk deformation hypothesis performs better in predicting the non-linear stage. The study’s findings can be applied to PMMA or other similar materials, and they can serve as a useful reference for future research on the contact mechanisms of other materials.

Quantitative reservoir evaluation and hydrocarbon volumetrics :an integrated petrophysical and 3-D static modeling approach in ‘Hamphidex’ field, Niger-Delta, Nigeria

The viability and hydrocarbon volumetrics of the Hamphidex field reservoirs were evaluated using petrophysical analysis and three-dimensional (3D) static modeling of three key reservoir sands (X, Y, and Z). The petrophysical analysis involved detailed characterization of lithologies, net-to-gross ratios, porosity, hydrocarbon saturation, water saturation, and permeability. The volumetric attributes of the field were derived from the constructed 3D reservoir models. Two geostatistical methods—Sequential Gaussian Simulation (SGS) and Sequential Indicator Simulation (SIS)—were utilized for facies and property modeling, focusing on porosity and permeability. Comprehensive facies models and petrophysical property models were developed for each reservoir sand (X, Y, and Z), integrating 3D visualizations of facies distributions, porosity, permeability, water saturation, and fluid contact models. These models also included cross-sectional views, providing a detailed spatial representation of the reservoir's characteristics. Two main facies, sand and shale, were identified, with sand acting as the hydrocarbon reservoir. Well-log correlations in six wells (Hamdex-02, 06, 01, 05, 04 and 07) of the field show that sands X, Y, and Z have lateral continuity and their hydrocarbon potential varies between wells. Specifically, Sand ‘X’ contains hydrocarbons in Hamdex-05 and Hamdex-07, Sand ‘Y’ is hydrocarbon-rich in Hamdex-05, and Sand ‘Z’ shows hydrocarbon presence in Hamdex-06, Hamdex-05, and Hamdex-04. These hydrocarbons are confined within an anticlinal structural trap that closes on a fault. For Sand X, Y, and Z, porosity values range between 16 and 25%, permeability varies from 10 to 1600 mD, water saturation lies within 7 to 50%, and hydrocarbon saturation spans from 50 to 93%. The volumetric assessments from the build models showed that Sand X, Y and Z have Stock-Tank-Oil-Initially-In-Place (STOIIP) of 75.864, 8.566 and 80.177 Million Barrels (MMbbl) respectively. In addition to Sand Y STOIIP it also has Gas-Initially-In-Place (GIIP) of 14.870 Billion Standard Cubic Feet (Bscf). The analysis indicates that the field contains substantial hydrocarbon reserves, with the petrophysical properties confirming the reservoirs' high quality. These findings demonstrate that the field is economically viable and suitable for development, presenting a strong potential for successful exploitation.

Modeling the Stress Field in MSLA-Fabricated Photosensitive Resin Components: A Combined Experimental and Numerical Approach

This study presents an experimental and numerical investigation into the stress field in cylinders manufactured from photosensitive resin using the Masked Stereolithography (MSLA) technique. For material characterization, tensile and bending test data from resin specimens were utilized. The stress field in resin disks was experimentally analyzed using photoelasticity and Digital Image Correlation (DIC) methods, subjected to compressive loads, according to the cylinder–plane contact model. Images were captured during the experiments using polarizing film and a low-cost CPL lens, coupled to a smartphone. The experimental results were compared with numerical and analytical simulations, where the formation of fringes and regions indicating the direction and magnitude of normal and shear stresses were observed, with variations ranging from 0.6% to 8.2%. The convergence of the results demonstrates the feasibility of using parts produced with commercially available photosensitive resin on non-professional printers for studying contact theory and stress fields. In the future, this methodology is intended to be applied to studies on stress in gears.

Coefficient of restitution in a low-velocity normal impact using elastoplastic contact stiffness

Abstract Elastoplastic deformation during particle impact occurs widely in many engineering applications. The material properties characterising both the elastic and plastic behaviour play an important role in particle impact. A non-linear contact stiffness-based model representing the elastic and plastic deformation of the material is used to obtain the coefficient of restitution (COR) during the impact of a sphere on a deformable substrate. The model consists of the Maxwell combination of perfectly plastic component and a non-linear elastic component. The proposed model is used to estimate the plastic energy dissipation during the impact. An analytical solution is obtained for residual contact radius and coefficient of restitution expressed in terms of experimentally determinable parameters. Our approach yields a single dimensionless parameter referred to as the ‘indentation parameter’, Λ, and it is shown that the impact response and coefficient of restitution for various impact situations can be determined based on this indentation parameter. The proposed model accurately predicts the residual contact radius and coefficient of restitution, validated through experimental results of low-velocity impacts (1 m/s - 4 m/s) over a flat sample of aluminium alloy (Al6061) impacted by steel and zirconia balls. The present model is further compared with other existing theoretical contact models for the elastoplastic impact and the extension of the present model for other dissipative systems is also discussed.

Contact force calculation and evolution analysis of granular systems based on micro-CT experiment

The calculation of inter-granule contact force in three-dimensional (3D) granular systems is a key and challenging aspect of granular mechanics research. Two elastic rubber balls are used as research objects for in-situ flat pressing Micro-CT experiments. Based on the Hertzian contact theory and Tatara large deformation contact theory, the contact model of elastic balls is verified, and the theoretical formula of the contact force of elastic balls based on the experiment is obtained. Taking the 3D granular systems as research object, in-situ probe loading experiment of micro-CT is carried out to obtain the 2D image sequence of the granules, after a series of digital transformations, the digital body images emerge, the contact force networks of the 3D granular systems under different loading conditions are obtained by constructing pore network models. The contact force distribution and evolution law of the granular systems are analyzed. The relation among the number of strong contacts, the distribution evolution, and the stability of the granular system is explored. The results show that the two elastic ball contact model conforms to the Hertzian contact theory and Tatara large deformation contact theory, and the contact force fitting formula based on experiment can characterize the contact force between two granules reasonably and effectively. The contact force of granules under probe loading is distributed in a net-like pattern starting from the contact point of the indenter and gradually transmitted to the lower and the surrounding area. The trend of average contact force is consistent with the trend of the contact times, showing a significant phase transition. With the increase of contact times, the frequency of particle compression increases, resulting in a greater contact force between granules, ultimately stabilizing at about 10.5 N. The number of strong contacts accounts for 45% to 50% of the total number of contacts, distributed throughout the whole granular system and supporting the network structure of the granular system. The larger values are concentrated below the indenter and exhibit a branching distribution. In the loading process, an equilibrium point is established at <i>z</i> = 14 mm, where the number of strong contacts reaches the peak. The network structure of strong contact force is spread throughout the entire 3D granular system, establishing the main skeleton that can withstand external loads. As the loading continues, the total value of strong contact forces increases, and their distribution within the granular system becomes more uniform.

Continuous Contact Model of Metamorphic Groove Pin Pair with Clearance

Continuous Contact Model of Metamorphic Groove Pin Pair with Clearance

Interfacial stress distribution analysis of natural fiberreinforced epoxy composites: a finite element approach

The strength of fiber-reinforced composites is greatly influenced by the bonding at the fiber-matrix interface. Experimental methods to study this interface are often challenging, making numerical approaches essential for evaluating the interfacial behavior in fiber-reinforced composites. This study investigates the stress and strain distribution in the fiber, matrix, and fibermatrix interface regions of natural fiber-reinforced single-fiber composites under tensile loading using the finite element method. Interface conditions were modeled using cohesive elements, with the composites represented in two dimensions through ABAQUS 6.14 software. The tie constrains contact model was employed to define binding interactions between the cohesive element, the fiber, and the matrix. The maximum stress value resulting from the simulation process is 202 MPa and a strain of 0.0449 mm. The stress is effectively distributed to the fiber, demonstrating that the cohesive element used in composite analysis under tensile loading serves as a reliable link between the fiber and the matrix. The simulation results revealed a maximum stress value of 202 MPa and a corresponding strain of 0.0449 mm. The stress distribution effectively transferred to the fiber, demonstrating the capability of cohesive elements to represent the interfacial bond in composites under tensile loading. These findings confirm that cohesive element modeling is reliable method for analyzing fibermatrix interactions in natural fiber reinforced composites, providing insights for optimizing composite performance.

Polyaxial Stress-Dependent Tensorial Permeability Variations of a Columnar Jointed Rock Mass: Insights from 3D Distinct Element Method

AbstractStress-induced permeability changes in geological formations have implications for various geo-engineering applications. Fractures are the predominant fluid flow pathways within tight geological formations such as jointed basaltic rock mass. Three-dimensional permeability variations of a columnar jointed rock mass using the three-dimensional block distinct element method are investigated. First, a permeability tensor is constructed by simulating fluid flow through joint sets using cubic law and determining its principal components through eigenvector analysis. Subsequently, numerical simulations are conducted to gain insights into the impact of principal stresses on joint deformation and fluid flow. Joint deformation is controlled by the joint’s mechanical properties, contact models, and stress state. Different isotropic and anisotropic stress loading scenarios on the permeability tensor evolution are investigated. The research findings indicate that isotropic stress loading led to uniform normal compression on fractures, resulting in a monotonic reduction of all permeability tensor components and increased permeability anisotropy. Anisotropic stress loading, particularly the intermediate stress, often neglected, significantly influenced the rock’s overall permeability. While the directions of the principal components remained consistent during isotropic loading, they exhibited changes in both direction and magnitude during anisotropic loading. The relative angle between the joint plane and the principal stresses was crucial in controlling joint behavior, including normal compression, shear failure, and dilation. Additionally, data analysis revealed that the exponential empirical model provided an excellent fit for permeability–stress data.

Calibration and Testing of Discrete Element Simulation Parameters for Ultrasonic Vibration-Cutter-Soil Interaction Model

This paper established an accurate discrete element for ultrasonic vibration-cutter-soil interaction model to study the interaction mechanism between the soil-engaging component and the soil. In order to reduce the interaction between calibration parameters and improve the calibration accuracy, it is proposed that the soil constitutive, contact parameters, and bonding parameters be calibrated by combining the soil repose angle experiment and the soil resistance experiment of ultrasonic vibration cutting. The study adopts the Hertz-Mindlin (no slip) contact model used in EDEM, to explore soil particle interactions. The central composite design is used to achieve systematic investigation. 3-factor 3-level orthogonal design experiment was employed using the coefficient of restitution, the coefficient of static friction, and the coefficient of rolling friction as key test factors and soil’s repose angle as the response index. Based on the Hertz-Mindlin with bonding contact model, Design-Expert 13.0 software was used to design the Plackett-Burman experiment, the steepest ascent, and the Box-Behnken experiment. With the maximum soil cutting resistance in ultrasonic vibration cutting experiment used as the response value, the adhesion parameters were optimized, and the optimal solution combination was obtained as: Normal Stiffness = 4.635 × 106 N/m, Shear Stiffness = 3.401 × 106 N/m, and Bonded Disk Radius = 2.57 mm. The optimal parameter combinations obtained from the calibration experiments were verified in two ways: ultrasonic vibration cutting and non-ultrasonic vibration cutting. The results showed that the errors between the simulation values and the actual values of the two comparative experiments were less than 5%, and the model calibrated for the three parameters can be used to study the drag reduction mechanism of ultrasonic vibration cutting in soil.

Calibration of Discrete Element Simulation Parameters and Model Construction for the Interaction Between Coastal Saline Alkali Soil and Soil-Engaging Components

There are approximately 36.7 million hectares of saline alkali land available in China. To enhance the comprehensive utilization value of coastal saline alkali land and boost crop yields in such areas, it is essential to conduct research on optimizing the operational performance of high-performance soil contact components in light of the soil characteristics of coastal saline alkali land. Discrete element simulation can be employed to investigate the operational mechanisms of various key components. Nevertheless, at present, there is a dearth of discrete element models for the key physical parameters and soil structure of coastal saline alkali land soil. In this article, typical coastal saline alkali field soil was sampled, and the physical properties of the saline alkali soil, including salt content, moisture content, particle size distribution, and particle size, as well as intrinsic parameters such as soil compaction, density, Poisson’s ratio, and shear modulus, were measured. The Hertz Mindlin with Bonding contact model was employed. Physical experiments on soil accumulation angles at different depths were carried out using the cylindrical lifting method. Subsequently, by means of the discrete element method and the BBD experimental design method, a response surface model was established, and an optimization analysis was performed on the optimal parameters for the soil–soil collision recovery coefficient, static friction coefficient, and dynamic friction coefficient at each depth. Test benches for measuring the collision recovery coefficient, static friction coefficient, and rolling friction coefficient of saline alkali soil at -65Mn were set up, calculation formulas for each parameter were derived, and the contact parameters between soil at different depths and 65Mn were obtained. The results of the sliding friction angle test on different depths of saline alkali soil at -65Mn were further verified using the discrete element method, with a maximum error of 3.11%, which falls within the allowable range. This suggests that the calibration results of the discrete element simulation parameters for the interaction between soil and contact components are reliable, providing data and model support for future research on enhancing the operational performance of high-performance contact components.