Discrete element modelling for wheel-soil interaction and the analysis of the effect of gravity

The discrete element method (DEM) is widely seen as one of the more accurate, albeit more computationally demanding approaches for terramechanics modelling. Part of its appeal is its explicit consideration of gravity in the formulation, making it easily applicable to the study of soil in reduced gravity environments. The parallel particles (P²) approach to terramechanics modelling is an alternate approach to traditional DEM that is computationally more efficient at the cost of some assumptions. Thus far, this method has mostly been applied to soil excavation maneuvers. The goal of this work is to implement and validate the P² approach on a single wheel driving over soil in order to evaluate the applicability of the method to the study of wheel-soil interaction. In particular, the work studies how well the method captures the effect of gravity on wheel-soil behaviour. This was done by building a model and first tuning numerical simulation parameters to determine the critical simulation frequency required for stable simulation behaviour and then tuning the physical simulation parameters to obtain physically accurate results. The former were tuned via the convergence of particle settling energy plots for various frequencies. The latter were tuned via comparison to drawbar pull and wheel sinkage data collected from experiments carried out on a single wheel testbed with a martian soil simulant in a reduced gravity environment. Sensitivity of the simulation to model parameters was also analyzed. Simulations produced promising data when compared to experiments as far as predicting experimentally observable trends in drawbar pull and sinkage, but also showed limitations in predicting the exact numerical values of the measured forces.     

Dynamic behaviour of an open lugged wheel under paddy soil conditions

Many experimental studies of open lugged wheel-soil interaction have been conducted, mainly based on the condition of constant slip and sinkage. As a result the reaction force to lugs seemed to be equal to the soil cutting resistance to a metal surface. However, analyses based on such methods do not appear to represent the actual behaviour of lugged wheel-soil interaction, especially when the lugs are spaced widely. The actual motion the wheel axle. In this study, an experimental device for a model lugged wheel was constructed to investigate the characteristics of the interaction between a lugged wheel and soil. Experiments were conducted under several test conditions of soil including paddy soil with a hard pan. The result of both theoretical and experimental data indicated that slip and sinkage of a lugged wheel showed a fluctuation with rotation angle of which the period is equal to the angular lug spacing. In each test soil condition used, the motion of the lugged wheel and the reaction forces acting on each lug from the soil for a free sinking wheel were different from that of the condition of constant slip and sinkage. It was found that the results obtained from this study could clarify the actual behaviour of lugged wheel-soil interaction.     

Experimental study on wheel-soil interaction mechanics using in-wheel sensor and particle image velocimetry Part I: Analysis and modeling of normal stress of lightweight wheeled vehicles

This study aims to develop a wheel-soil interaction model for a lightweight wheeled vehicle by measuring the normal stress distribution beneath the wheel. The main contribution of this work is to clarify the wheel-soil interaction using a wheel testbed that equips multiple sensory systems. An in-wheel sensor accurately measures the normal stress distribution as well as the contact angles of the wheel. Particle image velocimetry with a standard off-the-shelf camera analyzes soil flow beneath the wheel. The proposed model for the normal stress distribution is formulated based on these experimental data and takes into account the following phenomena for the lightweight vehicles that have not been considered in the classical model: (1) the normal stress distribution takes the form of a Gaussian curve; (2) the normal stress distribution concentrates in the front region of the wheel contact patch; (3) the distribution is divided into two areas with the boundary determined by the maximum normal stress angle; and (4) the maximum normal stress exponentially decreases as the slip ratio increases. Then, the proposed model is experimentally validated. Furthermore, a simulation study for the wheel driving characteristics using the proposed model confirms the accuracy of the proposed model.     

Prediction of wheel performance by analysis of normal and tengential stress distributions under the wheel-soil interface

Prediction of wheel performance by analysis of normal stress distribution under the wheel-soil interface was reported by one of our research members. In this study analysis of both normal and tangential stress distributions are included for the prediction of wheel performance. A visco-elastic soil model based on a three-element Maxwell model is used to evaluate normal stress distribution under a wheel running on soft ground. The values of the parameters characterizing the visco-elastic behavior of the soil can be derived from plate penetration tests. A rigid wheel-soil interface model is used to evaluate the tangential stress distribution under the wheel-soil interface. Shear deformation modulus, cohesion and angle of internal shearing resistance of the soil are derived from shear-displacement tests. Test results indicate that both maximum normal and shear stress occur in front of the wheel axle, and the location of peak normal stress shifts backwards towards the wheel axle while that of tangential stress shifts forwards when slippage is increased from a low value. Increasing slippage also causes a decrease in normal stress and an increase in tangential stress. Coefficients of traction and tractive efficiency are low at low slippage, increase with an increase in slippage, and level off at higher slippage.     

Some recent developments in vehicle—terrain interaction studies

This report presents a review of recent developments in the study of vehicle3-terrain interaction, based primarily on papers included in Session 2A of the 10th International Conference of ISTVS. Analytical studies of wheel-soil interaction using visco-elastic soil models and elasto-plastic and rigid-plastic finite element methods are reviewed. Results of experimental investigations of stress distributions on the wheel-soil interface as well as in the soil under various loading conditions are summarized. Various methods for predicting off-road vehicle performance are examined.     

Experimental and DEM analyses on wheel-soil interaction

In this paper, the wheel-soil interaction for a future lunar exploration mission is investigated by physical model tests and numerical simulations. Firstly, a series of physical model tests was conducted using the TJ-1 lunar soil simulant with various driving conditions, wheel configurations and ground void ratios. Then the corresponding numerical simulations were performed in a terrestrial environment using the Distinct Element Method (DEM) with a new contact model for lunar soil, where the rolling resistance and van der Waals force were implemented. In addition, DEM simulations in an extraterrestrial (lunar) environment were performed. The results indicate that tractive efficiency does not depend on wheel rotational velocity, but decreases with increasing extra vertical load on the wheel and ground void ratio. Rover performance improves when wheels are equipped with lugs. The DEM simulations in terrestrial environment can qualitatively reproduce the soil deformation pattern as observed in the physical model tests. The variations of traction efficiency against the driving condition, wheel configuration and ground void ratio attained in the DEM simulations match the experimental observations qualitatively. Moreover, the wheel track is found to be less evident and the tractive efficiency is higher in the extraterrestrial environment compared to the performance on Earth.     

Experimental study and analysis on driving wheels’ performance for planetary exploration rovers moving in deformable soil

Planetary rovers are different from conventional terrestrial vehicles in many respects, making it necessary to investigate the terramechanics with a particular focus on them, which is a hot research topic at the budding stage. Predicting the wheel-soil interaction performance from the knowledge of terramechanics is of great importance to the mechanical design/evaluation/optimization, dynamics simulation, soil parameter identification, and control of planetary rovers. In this study, experiments were performed using a single-wheel testbed for wheels with different radii (135 and 157.35 mm), widths (110 and 165 mm), lug heights (0, 5, 10, and 15 mm), numbers of lugs (30, 24, 15, and 8), and lug inclination angles (0°, 5°, 10°, and 20°) under different slip ratios (0, 0.1, 0.2, 0.3, 0.4, 0.6, etc.). The influences of the vertical load (30 N, 80 N, and 150 N), moving velocity (10, 25, 40, and 55 mm/s), and repetitive passing (four times) were also studied. Experimental results shown with figures and tables and are analyzed to evaluate the wheels’ driving performance in deformable soil and to draw conclusions. The driving performance of wheels is analyzed using absolute performance indices such as drawbar pull, driving torque, and wheel sinkage and also using relative indices such as the drawbar pull coefficient, tractive efficiency, and entrance angle. The experimental results and conclusions are useful for optimal wheel design and improvement/verification of wheel-soil interaction mechanics model. The analysis methods used in this paper, such as those considering the relationships among the relative indices, can be referred to for analyzing the performance of wheels of other vehicles.     

Longitudinal skid model for wheels of planetary rovers based on improved wheel sinkage considering soil bulldozing effect

To successfully deploy a wheeled mobile robot on deformable rough terrains, the wheel-terrain interaction mechanics should be considered. Skid terramechanics is an essential part of the wheel terramechanics and has been studied by the authors based on the wheel sinkage obtained using a linear displacement sensor that does not consider soil bulldozing effect. The sinkage measured by a newly developed wheel via detecting the entrance angle is about 2 times of that measured by the linear displacement sensor. On the basis of the wheel sinkage that takes the soil bulldozing effect into account, a linear function is proposed to the sinkage exponent. Soil flow in the rear region of wheel-soil interface is considered in the calculation of soil shear displacement, and its average velocity is assumed to be equal to the tangential velocity component of the transition point of shear stress. To compute the normal stress in the rear region directly, the connection of the entrance and leaving points is supposed as the reference of wheel sinkage. The wheel performance can be accurately estimated using the proposed model by comparing the simulation results against the experimental data obtained using two wheels and on two types of sands.     

Analysis and prediction of tyre-soil interaction and performance using finite elements

The finite element method [FEM] of analysis previously developed for prediction of rigid wheel-soil interaction is improved and extended to take into account (a) the effect of flexibility of tyre carcass where energy losses now occur in development of mobility, (b) a simpler requirement for specification of boundary condition using input loading, and (c) normal and tangential load stress from the tyre distributed across the tyre-soil interface and varying with slip. The comparisons of analytically computed (predicted) drawbar pull with actual experimentally obtained drawbar pull results for tests in three types of tyres show good correlations. The effect of inflation pressure on development of tyre deformation energy losses can be seen from the analytically computed values.     

Prediction of soil vertical stress under off-road tire using smoothed-particle hydrodynamics

Evaluation of vertical stress distribution in clay-loam soil using Smoothed-Particle Hydrodynamics-Finite Element Analysis (SPH-FEA) technique is presented in this research. The moist soil is modelled using the hydrodynamic elastic–plastic material and Murnaghan equation of state, while the tire is modelled using FEA in Visual Environment’s Pam-Crash software. Soil-tire interaction is performed using the node symmetric node to segment with edge treatment method. A single-wheel tester in a soil bin environment was utilized to provide experimental data. The objectives of the experimental test were to (1) calculate maximum subsoil stresses in the subsoil at 1 to 15 passes of a wheel with loads of 1, 2, 4, and 5 kN on soil with moisture levels of 0, 10, 17, and 24%.; (2) calibrate soil with different levels of moisture (3) compare predicted soil stresses with experiments. The maximum stress at 20 cm depth increased with increasing soil moisture and also with high levels of tire load. In contrast, successive traffic showed a decreasing effect on soil stress. The coefficient of determination 0.97 shows the predictions agreed very well with experiments. The moist soil-tire interaction model will be further used to analyze the soil stress in different soil depths and different forward velocity.     

Visualization and analysis of wheel camber angle effect for slope traversability using an in-wheel camera

This paper visualizes and analyzes an effect of a wheel camber angle for the slope traversability in sandy terrain. An in-wheel camera developed in this work captures the wheel-soil contact phenomenon generated beneath the wheel through a transparent section of the wheel surface. The images taken by the camera are then analyzed using the particle image velocimetry. The soil flows with various wheel camber angles are analyzed with regard to the soil failure observed on the slope surface. The analysis reveals that the slope failure and soil accumulation in front of the wheel significantly affect the wheel forces and distributions of the wheel sinkage in the wheel width direction. Further, the side force of the wheel in traversing a slope decreases as the slip ratio increases because the shear stress in the slope downward direction decreases owing to the slope failure.     

2D FE–DEM analysis of tractive performance of an elastic wheel for planetary rovers

We have been developing a simulation program for use with soil–wheel interaction problems by coupling Finite Element Method (FEM) and Discrete Element Method (DEM) for which a wheel is modeled by FEM and soil is expressed by DEM. Previous two-dimensional FE–DEM was updated to analyze the tractive performance of a flexible elastic wheel by introducing a new algorithm learned from the PID-controller model. In an elastic wheel model, four structural parts were defined using FEM: the wheel rim, intermediate part, surface layer, and wheel lugs. The wheel rigidity was controlled by varying the Young’s Modulus of the intermediate part. The tractive performance of two elastic wheels with lugs for planetary rovers of the European Space Agency was analyzed. Numerical results were compared with experimentally obtained results collected at DLR Bremen, Germany. The FE–DEM result was confirmed to depict similar behaviors of tractive performance such as gross tractive effort, net traction, running resistance, and wheel sinkage, as in the results of experiments. Moreover, the tractive performance of elastic wheels on Mars was predicted using FE–DEM. Results clarified that no significant difference of net traction exists between the two wheels.     

Experimental validation of distinct element simulation for dynamic wheel–soil interaction

The deformation behaviour of the soil during dynamic wheel–soil interaction was studied by using the discontinuum modelling technique, distinct or discrete element method (DEM). The simulation model was developed using DEM for two types of soil, soil-A (coarse sand) and soil-B (medium sand). A transparent sided soil bin was used to observe the soil deformation. Three CCD video camera photographic images of the validation experiments were analyzed and compared with the simulation program results.This paper presents the simulation and validation results for two types of soil at three different vertical loadings of 4.9, 9.8 and 14.7 N. Wheel sinkage, vertical and horizontal draft force acting on the rigid wheel and the soil deformation images from the validation experiments were some of the data used to compare the simulation program results with the validation experiments. The simulation program was helpful to understand the complex deformation behaviour of the soils. The simulated results for the deformation behaviour of soil-B showed better correlation with the validation experiments than soil-A. The results obtained have also been compared with the previous work on DEM to explain phenomena such as the high simulated sinkage of the rigid wheel.     

Linear normal stress under a wheel in skid for wheeled mobile robots running on sandy terrain

Wheeled mobile robots are often used on high risk rough terrain. Sandy terrains are widely distributed and tough to traverse. To successfully deploy a robot in sandy environment, wheel-terrain interaction mechanics in skid should be considered. The normal and shear stress is the basis of wheel-soil interaction modeling, but the normal stress in the rear region on the contact surface is computed through symmetry in classical terramechanics equations. To calculate that directly, a new reference of wheel sinkage is proposed. Based on the new reference, both the wheel sinakge and the normal stress can be given using a quadratic equation as the function of wheel-soil contact angle. Moreover, the normal stress can be expressed as a linear function of the wheel sinkage by introducing a constant coefficient named as sand stiffness in this paper. The linearity is demonstrated by the experimental data obtained using two wheels and on two types of sands. The sand stiffness can be estimated with high accuracy and it decreases with the increase of skid ratio due to the skid-sinkage phenomenon, but increases with the increase of vertical load. Furthermore, the sand stiffness can be utilized directly to compare the stiffness of various sandy terrains.     

Discrete element method-based studies on dynamic interactions of a lugged wheel with granular media

Terramechanics plays an important role in determining the design and control of autonomous robots and other vehicles that move on granular surfaces. Traction capabilities, slippage, and sinkage of a robot are governed by the interaction of a robot’s appendage with the operating terrain. It is important to understand how the terrain flows under this appendage during such an interaction. In this work, dynamics of soil performance and locomotion performance of a lugged wheel travelling on soft soil are numerically investigated. Studies are conducted with a two-dimensional model by using the discrete element method to analyze the interactions between a lugged wheel and the soil. The soil performance is studied by examining the force distribution and evolution of force networks during the course of the wheel travel. For two different control modes, namely, slip-based wheel control and angular velocity-based wheel control, the performance parameters such as, sinkage, traction, traction efficiency, and power consumption of the wheel are compared for various wheel configurations. The findings of this work are expected to be useful for optimal design and control of the lugged wheel travelling on deformable surfaces.     

Effect of wheel loads on the elevated roads between retaining structures

Deformations and stresses in a saturated soil between two retaining structures under wheel loads are analyzed. The problem has a practical importance especially in case of elevated roads. Expressions are given in graphical forms for one, two and four point loads which are representation of one and two vehicles with different axle lengths. Results provide information of predict contacting stresses at the wheel-soil interface and lateral pressures for design of retaining walls.     

Development of in-wheel sensor system for accurate measurement of wheel terrain interaction characteristics

Planetary rovers need high mobility on a rough terrain such as sandy soil, because such a terrain often impedes the rover mobility and causes significant wheel slip. Therefore, the accurate estimation of wheel soil interaction characteristics is an important issue. Recent studies related to wheel soil interaction mechanics have revealed that the classical wheel model has not adequately addressed the actual interaction characteristics observed through experiments. This article proposes an in-wheel sensor system equipped with two sensory devices on the wheel surface: force sensors that directly measure the force distribution between the wheel and soil and light sensors that accurately detect the wheel soil surface boundary line. This sensor design enables the accurate measurement of wheel terrain interaction characteristics such as wheel force distribution, wheel–soil contact angles, and wheel sinkage when the powered wheel runs on loose sand. In this article, the development of the in-wheel sensor system is introduced along with its system diagram and sensor modules. The usefulness of the in-wheel sensor system is then experimentally evaluated via a single wheel test bench. The experimental results confirm that explicit differences can be observed between the classical wheel model and practical data measured by the in-wheel sensor system.     

Wheels' performance of Mars exploration rovers: Experimental study from the perspective of terramechanics and structural mechanics

The soft ground and stones on the surface of Mars may cause sinkage of and damage to the wheels of a Mars rover. Therefore, we analyzed the performance of a wheel-step Mars rover from the perspective of terramechanics and structural mechanics. Using China’s Mars rover wheel prototype and wheel-soil interaction testbed, we obtained the driving performance of the wheels of a Mars rover under various conditions, including full skidding conditions. The vertical load set in the tests was determined using the gravity of Mars and the mass of the wheel-step Mars rover. The results indicate that the driving torque, sinkage, and resistance coefficient have a linear relationship with the wheel vertical load. An analysis of the structural mechanical characteristics of the wheel of the Mars rover was conducted by testing the radial, axial, torsional, and deflection stiffness. We found that the wheel ribs can improve the stiffness of the wheel but may reduce its driving performance. The analysis methods and evaluation indices can be used to analyze the performance of the wheels of other Mars rovers. Furthermore, the findings of this study can be used to optimize wheel design and motion control of wheel-step Mars exploration rovers.     

Analysis of a resistance force for the locked-wheel of push-pull locomotion rovers using large subsidence

This paper provides a quantitative analysis of the resistance force of the locked-wheel for push-pull locomotion rovers using intentional sinkage. Our previous study has confirmed that push-pull locomotion using intentional subsidence at an initial position can contribute to improving the traveling performance. The key factor of this scheme is the resistance force of the locked-wheel. However, the resistance force at different sinkage conditions and wheel sizes (e.g., mass, width, and diameter) remains unclear, especially for the individual locked-wheel. The detailed investigation of this interaction can contribute to the accurate estimation of rover mobility. This paper, therefore, investigates the locked-wheel and soil interaction at different sinkage conditions experimentally, especially focusing on the intentional sinkage condition. Additionally, the resistance force is considered theoretically through the knowledge based on the soil flow patterns beneath the locked-wheel. The experimental results confirmed that the resistance force of the locked-wheel rose as the initial sinkage, wheel size, and weight increases. Furthermore, the theoretical calculation suggested the resistance force increased with a similar tendency of the experimental data.     

2D dynamics of a stratified elastic subsoil layer

Summary The aim of this contribution is to propose a 2D model for an elastic stratified subsoil layer, which takes into account the effect of inhomogeneity on the dynamic behaviour of a medium. The obtained equations describe the subsoil response to the time-dependent foundation soil loadings, and can be applied to the analysis of dynamic subsoil-structure interactions. It is shown under which conditions the effect of microproperties of a medium on its dynamic response cannot be neglected.