Interaction Finite Element Model Research Articles

Study on mechanical properties and vibration reduction effect of a novel mesh-type elastic cushion plate with orifice: Theory, simulation and test

Elastic cushion plates in fastener systems play a critical role in controlling vibrations and noise, yet enhancing their damping capacity while maintaining cost-efficiency remains a persistent challenge. This paper presents a Novel Mesh-Type Elastic Cushion Plate with Orifices (NMTECPO), featuring innovative structural modifications that leverage air-damping mechanisms to achieve superior performance. Unlike conventional elastic cushions, NMTECPO incorporates a unique mesh structure with variable orifices and chamber volumes, significantly boosting its damping properties. Key structural parameters influencing the stiffness and damping characteristics of the NMTECPO such as orifice aperture, chamber volume, loading frequency, and loading amplitude, are identified through theoretical derivation and validated by extensive comparative tests. A fluid–structure interaction finite element model further demonstrates enhanced damping performance by optimizing the orifice aperture and chamber volume. Dynamic simulations also reveal that NMTECPO exhibits superior vibration attenuation compared to traditional cushion plates. The proposed NMTECPO offers a novel approach by effectively integrating structural innovations and parameter optimization to improve vibration damping, providing a cost-efficient solution with potential for future engineering applications and structural refinement.

On the deformation of CRTS-Ⅲ ballastless track during the casting process of self-compacting concrete: Numerical simulations and random forest-based prediction models

CRTS (China railway track system)-Ⅲ slab track is a main structural type of ballastless track in Chinse high-speed railway (HSR), which requires precise positioning of prefabricated concrete track slabs. However, the cast-in-place process of the self-compacting concrete (SCC) filling layer induces a large deformation as the flowing SCC generates an up-lift force on the bottom of the track slab. Currently, the deformation mechanism and influencing factors of the track slab during the SCC casting process have not been thoroughly studied. In this work, a fluid-solid interaction finite element model (FEM) is created to calculate the upward deformation of the CRTS-Ⅲ track slab caused by the flowing SCC. The main influencing factors, including the SCC’s casting height and speed, SCC’s viscosity, constraint stiffness of the limiting beams, and vertical temperature gradient, have been systematically analyzed and discussed. The relationship between the influencing factors and the track slab’s deformation is obtained. Furthermore, machine learning prediction models using random forest algorithms are developed to evaluate the deformation and up-lift force. A dataset comprising 603 calculation cases obtained from FEM simulations is used to train the prediction models. The random forest-based models developed in this study give an accurate estimation of the deformation and up-lift force based on the SCC properties and construction parameters. In addition, the feature importance of each influencing factor is obtained and studied. The findings of this study can provide guidance for the control of the deformation in CRTS-Ⅲ track slab during the SCC construction process.

Numerical modeling techniques for noise emission of free railway wheels

AbstractIn this article, we consider the numerical prediction of the noise emission from a wheelset in laboratory conditions. We focus on the fluid–structure interaction leading to sound emission in the fluid domain by analyzing three different methods to account for acoustic sources. These are a discretized baffled piston using the discrete calculation method (DCM), a closed cylindrical volume using the boundary element method (BEM) and radiating elastic disks in a cubic enclosure solved with the finite element method (FEM). We provide the validation of the baffled piston and the BEM using measurements of the noise emission of a railway wheel by considering ground reflections in the numerical models. Selected space-resolved waveforms are compared with experimental results as well as with a fluid–structure interaction finite element model. The computational advantage of a discretized disk mounted on a baffle and BEM compared to FEM is highlighted, and the baffled pistons limitations caused by a lack of edge radiation effects are investigated.

SEISMIC MITIGATION EFFECT FOR LARGE-SPACE UNDERGROUND STRUCTURES CONSIDERING SPATIALLY VARYING SOIL PROPERTIES

The seismic response of the large-space underground structure (LSUS) is significantly influenced by the physical properties of the surrounding soil media, while the soil owns a strong spatial variability. This study proposes a seismic response analysis process of the soil-LSUS interaction system is proposed, which can consider the characteristic of the spatially distributed soil properties. The proposed process begins with establishing the spatially random field model of the soil properties using the improved latent space method. Then, the model is calibrated based on the real data and Bayesian approach, and the realization of the random field is accomplished. Further, the soil-LSUS interaction finite element (FE) model is established, which incorporating the soil physical properties generated from the random field. Finally, the nonlinear time-history analysis of the soil-LSUS interaction FE model is conducted. As an illustration of the proposed process, a typical LSUS located in Guangzhou is selected as an example, and the seismic mitigation measure which the lead-filled steel tube damper (LFSTD) is installed between the intermediate column and the top beam is adopted for the LSUS. The influence of the spatial variability of soil properties on the seismic mitigation effect of the LSUS is investigated. Results indicate that the spatial variability of the soil properties can cause a minor influence on the force and deformation of the intermediate column and the energy dissipation ratio between the LFSTD and structure, while it can bring a significant influence on the maximum deformation and force and the shape of the hysteresis loop of the LFSTD.

Finite element modeling and analysis of effect of preexisting cervical degenerative disease on the spinal cord during flexion and extension.

Recent studies have emphasized the importance of dynamic activity in the development of myelopathy. However, current knowledge of how degenerative factors affect the spinal cord during motion is still limited. This study aimed to investigate the effect of various types of preexisting herniated cervical disc and the ligamentum flavum ossification on the spinal cord during cervical flexion and extension. A detailed dynamic fluid-structure interaction finite element model of the cervical spine with the spinal cord was developed and validated. The changes of von Mises stress and maximum principal strain within the spinal cord in the period of normal, hyperflexion, and hyperextension were investigated, considering various types and grades of disc herniation and ossification of the ligamentum flavum. The flexion and extension of the cervical spine with spinal canal encroachment induced high stress and strain inside the spinal cord, and this effect was also amplified by increased canal encroachments and cervical hypermobility. The spinal cord might evade lateral encroachment, leading to a reduction in the maximum stress and principal strain within the spinal cord in local-type herniation. Although the impact was limited in the case of diffuse type, the maximum stress tended to appear in the white matter near the encroachment site while compression from both ventral and dorsal was essential to make maximum stress appear in the grey matter. The existence of canal encroachment can reduce the safe range for spinal cord activities, and hypermobility activities may induce spinal cord injury. Besides, the ligamentum flavum plays an important role in the development of central canal syndrome.Significance.This model will enable researchers to have a better understanding of the influence of cervical degenerative diseases on the spinal cord during extension and flexion.

Study on Temperature Distribution Law of Tunnel Portal Section in Cold Region Considering Fluid–Structure Interaction

The design of tunnels in cold regions contributes greatly to the feasibility and sustainability of highways. Based on the heat transfer mechanism of the tunnel surrounding rock–lining–air, this paper uses FEPG software to carry out secondary excavation and development, then the air heat convection calculation model is established by using a three-dimensional extension of the characteristic-based operator-splitting (CBOS) finite-element method and the explicit characteristic–Galerkin method. By coupling with the heat conduction model of the tunnel lining and surrounding rock, the heat conduction-thermal convection fluid–structure interaction finite-element calculation model of tunnels in cold regions is established. Relying on the Qinghai Hekashan tunnel project, the temperature field of the tunnel portal section is calculated and studied by employing the fluid–structure interaction finite-element model and then compared with the field monitoring results. It is found that the calculated values are basically consistent with the measured values over time, which proves the reliability of the model. The calculation results are threefold: (1) The temperature of the air, lining, and surrounding rock in the tunnel changes sinusoidally with the ambient temperature. (2) The temperature of each layer gradually lags behind, and the temperature variation amplitude of the extreme value of the layer temperature gradually decreases with the increase in the radial distance of the lining. (3) In the vicinity of the tunnel entrance, the lining temperature of each layer remains unchanged, and the temperature gradually decreases or increases with the increase in the depth. The model can be used to study and analyze the temperature field distribution law of the lining and surrounding rock under different boundary conditions, and then provide a calculation model with both research and practical value for the study of the temperature distribution law of tunnels in cold regions in the future.

Study on the Seismic Response Characteristics of Shield Tunnels with Different General Segment Assembly Methods and Widths

Shield tunnels assembled with general ring segments are widely used in urban areas. Segment assembly methods and widths cause changes in the mechanical properties of the structure and influence the seismic response of shield tunnels. To investigate the influence of the assembly method and width of the general ring segment on the seismic performance of a shield tunnel, a three-dimensional refined soil–structure dynamic interaction finite element model of the shield tunnel was established based on ABAQUS, and the mechanical response and joint deformation of the general ring lining under seismic loads were studied. The simulation results show the following: (i) The overall deformation of the tunnel lining is not significantly affected by the assembly method, and the difference is only 5.24% under a 0.4 g earthquake. (ii) The seismic responses of general ring tunnels with different assembly methods are quite different, and the mechanical properties of the shield tunnel assembled with the straight assembly method are better than those of the shield tunnel assembled with staggered joints, but the deformation of the structure is larger. Under the action of a 0.1 g earthquake, the radial force, circumferential force, and bending moment of the staggered 90° assembly tunnel are respectively reduced by 13.6%, 11.1%, and 17.8% compared with the staggered 45° assembly structure, but the maximum intra-opening deformation increases by 0.19, 0.58, and 2.4 mm, respectively. (iii) The internal force distribution of the bolts is controlled by the deformation of the joint; compared with the CF90 and TF assembled tunnels, the mechanical properties and deformation characteristics of the CF45 and CF90 assembled tunnels are more reasonable. (iv) The extrados and intrados joint opening deformation and shear dislocation of the 1.2 m wide general ring segment under the staggered assembly increase by 1.2 mm and 1.03 mm, respectively, compared with the 1.5 m wide segment, while the radial force, circumferential force, and bending moment are reduced by 24.4%, 36.5%, and 41.7%, respectively, indicating that the seismic performance of the shield tunnel with a segment width of 1.5 m is better than that of the shield tunnel with a width of 1.2 m.

A compact DEA-based soft peristaltic pump for power and control of fluidic robots.

Fluid-driven robotic systems typically use bulky and rigid power supplies, considerably limiting their mobility and flexibility. Although various forms of low-profile soft pumps have been demonstrated, they either are limited to specific working fluids or generate limited flow rates or pressures, making them ill-suited for widespread robotics applications. In this work, we introduce a class of centimeter-scale soft peristaltic pumps for power and control of fluidic robots. An array of high power density robust dielectric elastomer actuators (DEAs) (each weighing 1.7 grams) were adopted as soft motors, operated in a programmed pattern to produce pressure waves in a fluidic channel. We investigated and optimized the dynamic performance of the pump by analyzing the interaction between the DEAs and the fluidic channel with a fluid-structure interaction finite element model. Our soft pump achieved a maximum blocked pressure of 12.5 kilopascals and a run-out flow rate of 39 milliliters per minute with a response time of less than 0.1 second. The pump can generate bidirectional flow and adjustable pressure through control of drive parameters such as voltage and phase shift. Furthermore, the use of peristalsis makes the pump compatible with various liquids. To illustrate the versatility of the pump, we demonstrate mixing a cocktail, powering custom actuators for haptic devices, and performing closed-loop control of a soft fluidic actuator. This compact soft peristaltic pump opens up possibilities for future on-board power sources for fluid-driven robots in a variety of applications, including food handling, manufacturing, and biomedical therapeutics.

Influence analysis of overlying soil layer to seismic behavior of inclined liquefiable soil and pile interaction system

The statistical data of earthquake hazard indicate that the liquefaction and large lateral flow deformation of ground is one of the main reasons of pile foundation failure. It is significant to focus on the seismic dynamic behavior of pile foundation in inclined liquefiable soil. Based on the Opensees finite element platform, an integrated finite numerical model of inclined liquefiable soil-pile group-structure interaction is established in this study, according to the shaking table model test of inclined liquefiable soil. The nonlinearity of pile-soil contact and the shear localization of soil layer are considered in the numerical simulation. The comparison of the numerical simulation result with the shaking table test results are verified. In order to discuss the influence of different overlying soil layer on the seismic behavior of soil and pile system, a typical inclined liquefiable soil-pile group-structure interaction finite element model is established by the proposed numerical simulation method. The results show that the overlying soil layer has an obvious influence on the seismic dynamic behavior of the pore pressure, the horizontal displacement of soil and pile. The increase of the thickness of overlying soil can reduce the liquefaction degree of sand and improve the mechanical performance of pile foundation.

Multiscale Finite-Element Analysis of Damage Behavior of Curved Ramp Bridge Deck Pavement Considering Tire–Bridge Interaction Effect

The present study developed a three-dimensional (3D) multiscale modeling approach to investigate the mesoscale damage behavior of curved ramp bridge deck pavement subjected to tire loading. First, a full-scale tire–bridge interaction finite element (FE) model was established to solve the macroscopic response of deck pavement during vehicle turning. Subsequently, the 3D mesostructure of the asphalt concrete layer was reconstructed from X-ray computer tomography (CT) images through a digital image processing (DIP) technology. Finally, a homogenization method was adopted to link the mechanical properties of deck pavement at two different scales, and a mapping procedure was employed to transfer the boundary displacements from macroscale pavement target element to the mesoscale representative volume element (RVE) model. Also, the bilinear cohesive elements were applied to simulate damage initiation in the RVEs. The results showed that smaller curvature radius and higher travel speed can promote unbalanced stress and local damage of the outer loading zone. From the mesoscale viewpoint, reducing curvature radius or increasing travel speed could result in crack direction shifts and increase the irregularity of microcrack distribution. In addition, transverse microcracks are commonly found in the center of the tire load, while longitudinal microcracks are distributed on the deck pavement surface around the tire edges and the dual-tire clearance center.

Modified p-y curves for monopile foundation with different length-to-diameter ratio

The soil reaction of the monopile foundation subjected to lateral loading in offshore wind turbines is typically assessed relying on p-y curves advocated by API. However, this method is inadequate for gradually increasing monopile diameters and significantly underestimates the lateral soil confinement. In the present works, a 3D pile-soil interaction finite element model was first established, considering the soil suction and strain hardening characteristics for the normally consolidated clay in China’s sea. Modifications to the p-y curves in API were accomplished in the comparative process between the lateral soil resistance-displacement curves retrieved from the finite element model and the representative expression. Furthermore, the prediction accuracy for the corrected p-y curves has been proved by forecasting the monopile lateral bearing capacity with varying length-to-diameter ratios, which also demonstrates that the modified p-y curves could successfully reflect the lateral soil confinement of the normally consolidated clay and flexible piles. It also provides an approach to assess the deformation response and horizontal ultimate bearing capacity of monopiles with different length-to-diameter ratios.

Development of Empirical Relations to Predict Ground Vibrations due to Underground Metro Trains

The vibration caused due to underground metro trains in densely populated urban areas and its adverse effect on the nearby structures and occupants is a growing concern for the policymakers and government bodies, demanding its rapid and correct assessment. Metro bodies in India use Research Designs and Standards Organisation guidelines similar to Federal Transit Administration guidelines for vibration assessment. However, in the guidelines, there is limited discussion on soil geotechnical properties, tunnel depth and axle load, which considerably affect the magnitude of train-induced ground vibration. This study focuses on developing easily comprehendible empirical relations to predict metro train-induced ground vibrations considering the effect of these parameters. Multiple non-linear regression is used to establish the empirical relations utilising the datasets generated from a two-dimensional train-track-tunnel-soil dynamic interaction (TTTSDI) finite element model. The TTTSDI model is based on the two-step methodology available in the literature and is validated with field measurement results of the Delhi metro sites. The developed empirical relations predicted the ground vibrations with maximum and minimum errors of 2.66% and 0.84%, respectively.

Mitigating hydrodynamic ram effect by a novel perforated lattice filled tank

High-velocity impact on fluid-filled structures has become an increasingly critical issue for the vulnerability of aircraft fuel tanks. Hydrodynamic Ram (HRAM) is a typical phenomenon that occurs when a flying object impacts against a fluid-filled tank with a high velocity, leading to excessive structural damage and catastrophic failure. Thus, an effective HRAM mitigation method is highly desired to provide better protection for the aircraft fuel tanks. In this paper, a novel perforated lattice filled tank was proposed by filling the empty space inside a fuel tank to mitigate HRAM effect and reduce the shock wave pressure. Perforated lattices were utilized to reduce the cavity expansion and structural damage while the holes on the lattice walls could guarantee the connectivity for each unit cells. Experimental tests were conducted to investigate the impact response and cavity evolution of the perforated lattice filled tanks with 1 × 1, 3 × 3 and 5 × 5 unit cells. The deformations of the front, back and lateral walls of the lattice-filled tanks were recorded by using a high-speed video camera and strain gauges (SG). Besides, A fluid–structure interaction finite element model of the fuel tank based on the SPH method was established to study the mitigation performance of HRAM phenomenon. A parametric analysis was conducted to explore the influence of unit-cell number, impact velocity, volume fraction of fluid and liquid level pressure on the impact response. The results indicated that the structural deformation and damage can be alleviated by the lattice-filled design, and it could be considered as an effective solution to mitigate the HRAM effect.

Numerical investigation of load-induced fatigue cracking in curved ramp bridge deck pavement considering tire-bridge interaction

Load-induced fatigue cracking is one of the main distresses of bridge deck pavement; however, due to the complexity of curved ramp bridge structures, it is difficult to effectively consider this damage behavior in the deck pavement design procedure. To deal with this problem, a three-dimensional (3D) tire-bridge interaction finite element (FE) model was developed in the present study to investigate the load-induced fatigue damage behavior within the curved ramp bridge deck pavement. In this FE model, an inflated tire model was developed to simulate the rolling tire load, and the Yeoh model was used to define the constitutive relationship of the hyperelastic material. In addition, the extended finite element method (XFEM) was utilized to describe the crack initiation and evolution, and a direct cycle algorithm was used to calculate the fatigue crack propagation law of deck pavement under different structural design parameters. The results indicated that the developed tire-bridge interaction FE model furnishes an in-depth insight into the fatigue crack propagation within curved ramp bridge deck pavement. The consideration of the tire-bridge interaction and structural design parameters has a significant influence on the accurate prediction of the fatigue crack evolution law within different deck pavement locations. Because the developed full-scale FE model is capable of simultaneously considering multiple factors, including the curvature radius, transverse and longitudinal slope, traffic speeds and tire-bridge interaction loads, it can be expected to be a mechanistic tool to facilitate and enhance fatigue damage analysis of curved ramp bridge deck pavement.

An Accurate and Efficient Fitness-For-Service Assessment Method of Pipes with Defects under Surface Load

With continued urbanization in China, the construction of urban gas pipelines is increasing, and the safety of gas pipelines are also increasingly affected by urban development and the increased scope of buildings and roads. Pipes with defects are more likely to fail under the surface loads. In this study, uniaxial tensile tests of high-density polyethylene (HDPE) pipes were carried out to obtain the real material parameters of pipe. A pipeline-soil interaction finite element model of HDPE pipeline with defects under surface load was established. The failure mechanism of the urban gas pipeline was studied and the influence of parameters such as internal pressure, defect position, defect depth on the mechanical behavior, and failure of pipelines were analyzed. A failure criterion for HDPE pipes with defects under surface load was proposed based on the limit-state curves obtained under different working conditions. Furthermore, an accurate and efficient fitness-for-service assessment procedure of pipes with defects under surface load was proposed. The results showed that maximum Mises stress of the pipeline gradually increased with increasing surface load and the position of maximum stress changed from the top and bottom of the pipe to the defect position and both sides of the pipe. Finally, when Mises stress of the HDPE pipe exceeds the yield limit, failure will occur. Internal pressure, defect location, and defect depth were found to influence the failure process and critical surface load of the pipeline. Safety evaluation curves of the gas pipeline with defects under surface load were obtained by calculating the critical failure load of the pipeline under various working conditions. Finally, a nonlinear fitting method was used to derive a formula for calculating the critical surface load under different defect parameters. The proposed method provides a useful reference for urban gas pipeline safety management.

Development and Validation of TSFEM to Blunt Impact Based On Artificial Intelligence and Information Technology

The numerical simulation method is highly repeatable and has a short research cycle, which has been widely popular in predicting blunt ballistic impact to the thorax. A Thorax Simplified Finite Element Model (TSFEM) of projectile-structure interaction was developed by using ABAQUS software. The projectiles made of plastic stick which were used in Bir’s whole body impact testing. The force dynamic response, deflection dynamic response and force - deflection dynamic response were obtained, and validated with the test response corridors of the human body (PMHS) announced by the researchers of Wayne State University, we confirmed the finite element model can be used to predict the ballistic blunt thorax injury risk, in order to provide certain reference to further formulate and assess the design criterion and injury criterion of non-lethal kinetic energy weapons.

Nonlinear multiscale analysis of coronary atherosclerotic vulnerable plaque artery: fluid-structural modeling with micromechanics.

A unique three-dimensional (3D) computational multiscale modeling approach is proposed to investigate the influence of presence of microcalcification particles on the stress field distribution in the thin cap layer of a coronary atherosclerotic vulnerable plaque system. A nested 3D modeling analysis framework spanning the multiscale nature of a coronary atherosclerotic vulnerable plaque is presented. At the microscale level, a micromechanical modeling approach, which is based on computational finite-element (FE) representative unit cell, is applied to obtain the homogenized nonlinear response of the calcified tissue. This equivalent response effectively allows the integration of extremely small microcalcification inclusions in a global biomechanical FE model. Next, at the macroscale level, a 3D patient-based fluid-structure interaction FE model, reconstructing a refined coronary artery geometry with calcified plaque lesion, is generated to study the mechanical behavior of such multi-component biomechanical system. It is shown that the proposed multiscale modeling approach can generate a higher resolution of stress and strain field distributions within the coronary atherosclerotic vulnerable plaque system and allow the assessment of the local concentration stress around the microcalcifications in plaque cap layers. A comparison of stress field distributions within cap layers with and without inclusion of microcalcifications is also presented.

High-Fidelity Finite Element Modeling and Analysis of Adaptive Gas Turbine Stator-Rotor Flow Interaction at Off-Design Conditions

ABSTRACTThe objective of this work is to computationally investigate the impact of an incidence-tolerant rotor blade concept on gas turbine engine performance under off-design conditions. When a gas turbine operates at an off-design condition such as hover flight or takeoff, large-scale flow separation can occur around turbine blades, which causes performance degradation, excessive noise, and critical loss of operability. To alleviate this shortcoming, a novel concept which articulates the rotating turbine blades simultaneous with the stator vanes is explored. We use a finite-element-based moving-domain computational fluid dynamics (CFD) framework to model a single high-pressure turbine stage. The rotor speeds investigated range from 100% down to 50% of the designed condition of 44,700 rpm. This study explores the limits of rotor blade articulation angles and determines the maximal performance benefits in terms of turbine output power and adiabatic efficiency. The results show articulating rotor blades can achieve an efficiency gain of 10% at off-design conditions thereby providing critical leap-ahead design capabilities for the U.S. Army Future Vertical Lift (FVL) program.

Finite element analysis on mechanical behavior of semi-exposed pipeline subjected to debris flows

Before the debris flow occurs, the continuous rainfall may cause the pipeline to be partially exposed (semi-exposed pipeline). The semi-exposed pipeline will be under great threat when the debris flow comes. To study the mechanical behavior of semi-exposed pipeline subjected to the debris flow, a finite element model of interaction between soil and semi-exposed pipeline under the impact of debris flow was established for the practical working conditions of semi-exposed pipeline. Effects of relevant parameters were analyzed, including the velocity of debris flow, impact angle of debris flow, massive stone, and parameters of corrosion pit (i.e. the depth, length, and width of corrosion pit). The results show that the velocity and the impact angle of debris flow and massive stone in the debris flow have a significant influence on the stress and deformation of the pipeline, but the exposed length of pipeline has little impact. Non-pressure pipeline can withstand the greater velocity of debris flow. In addition, parameters of corrosion pit (depth, length, and width) all have some impact on pipeline stress. Based on the result of multiple regression analysis, the effects of relevant parameters of corrosion pit on the maximum Von Mises stress of pipeline are ranked as follows: corrosion depth (A) > corrosion width(B) > corrosion length (L). People can propose an effective protection strategy for the pipeline subjected to debris flow based on the research results.

Finite element modeling of interaction between non-pneumatic mechanical elastic wheel and soil

The tire/soil interaction is an important and complex research topic in vehicle-terrain mechanics, which has a great influence on vehicle mobility and soil compaction. The purpose of this work was to develop a detailed nonlinear finite element model for parametric analysis of the interaction between mechanical elastic wheel and the soil. The three-dimensional finite element model of the mechanical elastic wheel, which considers material nonlinearity, geometric nonlinearity, as well as large contact deformation between the mechanical elastic wheel and soil, was developed based on the physical model. The reliability and accuracy of the finite element model were validated by comparing the predicted wheel static loading characteristics in rigid plane with those of experiments. The finite element model of the soil is modeled as a two-layer system and the Extended Drucker–Prager model was used to describe the nature soil properties. The deformation and stress of the soil and mechanical elastic wheel under the static loading condition were studied in detail. Moreover, the effects of the axial load on mechanical elastic wheel/soil interaction were also analyzed. The research results could offer reference for an optimum structural design of a mechanical elastic wheel and a prediction of soil compaction caused by non-pneumatic tires.