Engineering Design, Kinematic and Dynamic Analysis of High Lugs Rigid Driving Wheel, a Traction Device for Conventional Agricultural Wheeled Tractors

: Traction failure on loose terrain is common in conventional agricultural wheeled tractors due to poor traction ability and lower power transfer efﬁciency of drive wheels, which leads to excessive energy consumption and soil compaction in agriculture. To overcome the problem, this paper presents a new design of a rigid lugged wheel for use in ﬁeld tillage operations. This wheel was designed according to ﬁeld operational requirements and also provided with rubber pads for smooth on-road drives. Kinematic and dynamic analysis of new wheel designs were carried out to study how they move, how they interact with the soil, and how they generate drive force in loose terrain soil. The relationship of wheel lug motion trajectories, displacement, and velocity of the wheel relative to ﬁeld conditions, different travel reduction rates, and lug penetration/wheel sinkage were analyzed. Wheel-terrain interaction and shear stress-shear displacement relationships when the wheel is driven in soft, deformable terrain were studied using classic soil mechanics principles. It is found that the component of thrust in the direction of driving, i.e., driving force, is ranged between 81.52% and 86.17%, while the vertical component is reported to be less than 30% and further decreases to 9%, which is the compaction avoiding factor. The relationships, thus developed, of wheel parameters, soil stress and thrust characteristics, and wheel drive force were derived and revealed that the traction performance, power transfer efﬁciency, and trafﬁcability of tractors in loose terrain can be improved by using the newly proposed wheel. A ﬁnite element method was used to analyze the designed wheel model for structural stability and optimization. The theoretical analysis results of the new drive wheel are convincing, so further tests and ﬁeld operation research are recommended for sustainable adoption.


Introduction
Conventional agricultural tractors are prominent off-road ground vehicles and the main power source in field tillage operations for crop cultivation.Tillage and ploughing are the major cultivation operations, so the traction demand for tillage is increasing with the development of mechanized agriculture [1].Consequently, the size of tractors has increased in recent years, which requires maintaining the power-to-weight ratio of tractors for efficient operations, soil management, and energy conservation [2].Weight control for all agricultural vehicles, particularly tractors, is strongly recommended for arable soil management and for reducing the energy demand of crop production [1,[3][4][5].
Tractor traction or drawbar force development through wheel-soil interaction is a complex process that adversely affects the terrain, soil, and energy demand of the trac-tor [1].Researchers around the world are working to find solutions for efficient traction development to prevent soil compaction and promote energy conservation.All the attempts that were made for the purpose are adhesion-based and focused on increasing the wheel-terrain contact area and normal axle loads at contact points [3,[6][7][8].All-wheel drives, dual and multiple wheels, ballasting techniques, inflation pressure adjustments, and split and continuous type rubber tracks are examples of practiced methods.Due to their size, weight, and resistance to movement, these methods cannot be used in agricultural field operations [4].
Lugged-type rigid wheels have some advantages over tire wheels for traction intensification on loose deformable terrains.Many forms of lugged wheels are in use in agriculture, especially in low-land farming and paddy fields, on conventional tractors and boat-type tractors that have limited applicability in some specific conditions [9].Some new types of moveable lugged wheels and transformable walking wheels were developed for ground vehicles that are specific to small machines and not applicable to conventional agricultural tractors [10,11].To further enhance the applicability and advantage of lug-type rigid wheels, it is required to improve their structure and design characteristics, with a special focus on agricultural field conditions.
With regard to the agricultural environment and field conditions, energy crises, agricultural sustainability issues, soil degradation caused by mechanical compaction issues (mainly by machine weight), and environmental quality problems, it is of great importance to design highly efficient traction wheels for agricultural tractors [9,[11][12][13].For designing and studying agricultural vehicles, especially agricultural wheeled tractors, it is important to consider machine weight, field mobility, traction ability (pull ability), soil trafficability, and soil behavior under applied wheel forces.It is more effective to use terramechanics principles, soil mechanics theories, and vehicle engineering principles to make the operations energy-efficient, economical, and environment-friendly [14][15][16][17].
A new design of the traction device "rigid lugged wheel" presented in this paper was designed for the high traction demand of intensive tillage operations with a similar mounting mechanism to tire wheels [4,18].It is easily replaceable with tire wheels and also provides rubber pad accessories to use for on-road driving.Rubber pad assemblies can be mounted and dismounted to meet field and road drive requirements.It is used in the field without rubber pads, where sharp edge high lugs (150 mm) penetrate to soil and apply exertion pressure to the terrain soil, causing thrust in response to the lug working surface.The horizontal component of this thrust force is the driving force of the wheel that moves the tractor in the forward direction and also provides drawbar force to draw the tillage implements/ploughs through the soil.The lugged wheel is intended to generate traction force through insertion pressure rather than shearing pressure, and the majority of vertical pressure is applied in the lateral direction.The thrust or drive force of a rigid lugged wheel is independent of normal axle load and only depends on the soil thrust (soil strength) and the area of the working faceplate of the lug.Theoretical analysis is conducted for the newly designed rigid lugged wheel that provides in-depth theoretical modeling to assess and confirm its working principles.It will elaborate on how this wheel would help to maximize the traction ability of agricultural wheeled tractors and how it interacts with the terrain.The FEM analysis will help optimize the model structure of the wheel for lighter weight and longer life for feasible adoption.

Design Concept
A wheel prototype with a new design structure for a high-lug rigid wheel was proposed.An easily mountable and replaceable lugged wheel, which is a novel mechanism for agricultural wheeled tractors, ensures full-time four-wheel drive in complex off-road terrain and diversified agricultural environments.Additionally, it is convenient to control the traction mechanism of wheeled tractors.As the rigid lugged wheel was designed to drive tractors on soft soil surfaces, due to the penetration of lugs into the soil, it may face some resistance force.The rigidly lugged wheel derives its driving force from the lugs' penetration into the soil, resulting in a pre-tillage effect on the terrain soil beneath the tracks.
When the Agricultural wheeled tractor is driven on soft terrain, the wide surface rim of the rigid lugged wheel contacts the soft surface to support most of the weight and avoid soil compaction.The trafficability of Agricultural wheeled tractors can be improved due to the larger contact area of the newly designed rigid lugged wheel as rear drive wheels.When Agricultural wheeled tractors are driven on hard terrain or paved roads, the road pad accessory components are mounted to the wheel and the rear wheels are adjusted to contact the ground for a smoother drive.

Design Requirements
The rigid lugged wheel as a driving wheel for the agricultural tractor can adjust its status based on the driving area and the different characteristics of the agricultural field soil.On soft terrain, due to the height of the lugs, the wheel will improve the tractive force performance, field operation energy consumption, and trafficability of the field soil.The rigid lugged wheel's main requirement is field operations and traction performance on soft, relatively weaker soils with less trafficability.
The height of the lugs of the rigid lugged wheel is relative to the normal tillage depth of 120 to 180 mm, so the lugs will get supported by the subsurface, and with a higher contact area, the ground pressure on soft soil will get reduced, which will help to avoid soil compaction.The working surface or working faceplate of the lug is the object of study in the design process.The driving force and traction ability of the wheel depend on the working surface's contact with the soil; as much as the lugs penetrate the soil, they deliver the driving force through soil contact.
The design parameters of the rigid lugged wheel should be determined in such a way that it can produce maximum driving force with optimum energy consumption and minimum damage to the working environment.Rim width, lug length, lug height, and lug angle, all should be in optimum working performance position.The rim setting should be adjusted to withstand wet soil working and avoid soil sticking and wheel blocking during wet soil field operations.
Based on the above considerations, three significant problems must be solved.First, a mechanism should be designed to transform the wheel into a lightweight and highly efficient to generate driving force, whereas by implementing tracks the weight of the tractor and internal motion resistance increased by a significant amount.Second, the power transmission or power transfer mechanism to pull the tillage implement must be solved.Third, these two mechanisms can be integrated to reduce the amount of energy and space required and also avoid soil compaction in agriculture.

Structure and Working Principle
The model of the proposed rigid lugged wheel was developed in Pro-E, and MAT-LAB was employed for the solution of the design equation and design calculations.This computerized model of a rigid lugged wheel is held to make further design calculations, justifications and confirmations, kinematic and dynamic analysis, and drive force calculations.Pro-E software is an element-based tool that first develops the parts individually and then unites them by joining the parts in patterns to make the full assembly.The preliminary design of the rigid lugged wheel is presented in Figure 1, where: 1.The rim of the wheel, 2. Lug (150 mm) and lug working faceplate (the main research object), 3. Rubber pad (for on-road drive) and fixing assembly, 4. Wheel axis, 5. Lug angle (α = 30 • ), 6.The direction of motion or rolling of the wheel.
calculations.Pro-E software is an element-based tool that first develops the parts individually and then unites them by joining the parts in patterns to make the full assembly.The preliminary design of the rigid lugged wheel is presented in Figure 1, where: 1.The rim of the wheel, 2. Lug (150 mm) and lug working faceplate (the main research object), 3. Rubber pad (for on-road drive) and fixing assembly, 4. Wheel axis, 5. Lug angle (α = 30 °), 6.The direction of motion or rolling of the wheel.

Wheel Structure
The whole wheel structure of the rigid lugged wheel consists of three main functional components.(1) Rim, including internal strengthening supports and mounting accessories.(2) Lugs that are arranged on the rim in herringbone arrangement at a 30 ° mounting angle with the axis of the wheel and lug height of 150 mm were set for optimum traction and driveability performance.(3) Road pad/rubber pad assemblies with internal supports and outer rubber pad that can be mounted to the wheel for on-road drive.
The main function and purpose of a rigid lugged wheel are field drive and adequate traction force provision for high-draught implements in diverse agricultural field conditions.Therefore, the road pad assemblies are only accessory items for minor road drives, while the rim and lugs are two major components for the functional requirements.The main design considerations were given to the two components, and some details of the designs are as follows: The rim of the rigid lugged wheel is designed according to mounting provisions available for the tire wheel, as the same is used for rigid lugged wheel mounting, and according to available space and the drive mechanism.The width of the rigid lugged wheel rim was kept as wide as possible according to the drive criteria of the tractor and optimum working performance.Rim setting was adjusted to withstand all terrain/soil conditions and all working environments, such as dry and hard-untilled surfaces, soft loose tilled soil, and wet or saturated soil working conditions, and to avoid soil sticking and wheel blocking during wet soil field operations.
The main research object of this study was the lug of the rigid lugged wheel, so the design parameters of the lug are of much importance.On the lug, the working faceplate was the more specific object to design, arrange, and analyze for optimum wheel performance.The driving force of the lug is dependent on working surface parameters such as length and width, the mounting angle of the lug with the wheel axis, the angle of the lug's working faceplate according to the reference axis, and its interaction with soil (soil penetration).Therefore, consideration was given to lug design, and analysis was performed for lug design and lug arrangements.

Wheel Structure
The whole wheel structure of the rigid lugged wheel consists of three main functional components.(1) Rim, including internal strengthening supports and mounting accessories.
(2) Lugs that are arranged on the rim in herringbone arrangement at a 30 • mounting angle with the axis of the wheel and lug height of 150 mm were set for optimum traction and driveability performance.(3) Road pad/rubber pad assemblies with internal supports and outer rubber pad that can be mounted to the wheel for on-road drive.
The main function and purpose of a rigid lugged wheel are field drive and adequate traction force provision for high-draught implements in diverse agricultural field conditions.Therefore, the road pad assemblies are only accessory items for minor road drives, while the rim and lugs are two major components for the functional requirements.The main design considerations were given to the two components, and some details of the designs are as follows: The rim of the rigid lugged wheel is designed according to mounting provisions available for the tire wheel, as the same is used for rigid lugged wheel mounting, and according to available space and the drive mechanism.The width of the rigid lugged wheel rim was kept as wide as possible according to the drive criteria of the tractor and optimum working performance.Rim setting was adjusted to withstand all terrain/soil conditions and all working environments, such as dry and hard-untilled surfaces, soft loose tilled soil, and wet or saturated soil working conditions, and to avoid soil sticking and wheel blocking during wet soil field operations.
The main research object of this study was the lug of the rigid lugged wheel, so the design parameters of the lug are of much importance.On the lug, the working faceplate was the more specific object to design, arrange, and analyze for optimum wheel performance.The driving force of the lug is dependent on working surface parameters such as length and width, the mounting angle of the lug with the wheel axis, the angle of the lug's working faceplate according to the reference axis, and its interaction with soil (soil penetration).Therefore, consideration was given to lug design, and analysis was performed for lug design and lug arrangements.

Design of Lug Working Faceplate
Figure 2a shows the fixed coordinate system "xyz", where: the x-axis is the iron wheel axis, the y-axis points to the tractor's forward direction, and the z-axis is vertically upward.In the "x o y o z o " coordinate system, that is obtained by rotating "xyz" about 30 degrees in the positive direction (counter-clockwise) around the z-axis.

Design of Lug Working Faceplate
Figure 2a shows the fixed coordinate system "xyz", where: the x-axis is the iron whee axis, the y-axis points to the tractor's forward direction, and the z-axis is vertically up ward.In the "xoyozo" coordinate system, that is obtained by rotating "xyz" about 30 de grees in the positive direction (counter-clockwise) around the z-axis.In this way, the z-axis remains the same, and the transformation relationship betwee the coordinate systems "xyz" and "xoyozo" is as follows:  =  cos  −  sin   =  cos  +  sin   =  (1 In the "   xoyozo" coordinate system, the small ellipse at the lug base of the iro wheel (as shown in Figure 2b) corresponds to the following equation: where: R = 1276/2 = 638 mm is the radius of the iron wheel rim; α = 30°; the small ellipt equation ( 2), now written as a parametric equation:  = .sin   = 0  = .cos  (3 In Equation (3), θ is a small ellipse parameter, and the equation of a small ellipse i the xyz coordinate system is: The range of x in the equation (according to Figure 2b) is 345  140.

Mounting Angle of the Lug Working Faceplate
As shown in Figure 3a,b, line 1: iron wheel axis; plane 4: is a plane perpendicular t line 1 wheel axis; point p1: is the intersection of plane 4 and line 1 wheel axis; point p2: the intersection of plane 4 and smaller elliptic curve; line 2: the intersection of point p and point p2.In this way, the z-axis remains the same, and the transformation relationship between the coordinate systems "xyz" and "x o y o z o " is as follows: In the "x o y o z o " coordinate system, the small ellipse at the lug base of the iron wheel (as shown in Figure 2b) corresponds to the following equation: where: R = 1276/2 = 638 mm is the radius of the iron wheel rim; α = 30 • ; the small elliptic equation ( 2), now written as a parametric equation: In Equation (3), θ is a small ellipse parameter, and the equation of a small ellipse in the xyz coordinate system is: The range of x in the equation (according to Figure 2b) is 345 ≥ x ≥ 140.

Mounting Angle of the Lug Working Faceplate
As shown in Figure 3a,b, line 1: iron wheel axis; plane 4: is a plane perpendicular to line 1 wheel axis; point p1: is the intersection of plane 4 and line 1 wheel axis; point p2: is the intersection of plane 4 and smaller elliptic curve; line 2: the intersection of point p1 and point p2.line 3: intersect line 2 at point p2 and go along the wheel lug surface and plane 4; angle gamma (γ): mounting angle of iron wheel lug faceplate; that is the angle between line 2 and line 3 at vertex point p2.
In the coordinate system xyz, line 3 is parallel to the z-axis, so the angle between line 2 and line 3 is the angle between the vector and the z-axis.
That is, the installation angle of the lug is equal to the parameter theta (θ) in the small elliptic parametric equation.From Equation ( 4): In the above equation, the value range of x is [140,345]; α = 30°; R = 638 mm; The width of the lug along the axis of the iron wheel is 205 mm (see Figure 2b).The variable x in Equation ( 7) is replaced by the lug width parameter w (range: 0 to 205): According to Equation ( 8), the relationship curve between the installation angle of the lug and the width "w" along the axis of the wheel can be drawn, as shown in Figure 4.In Figure 3, the obliquely upward direction vector of line 2 obtained from Equation ( 4) is as follows: In the coordinate system xyz, line 3 is parallel to the z-axis, so the angle between line 2 and line 3 is the angle between the vector and the z-axis.
That is, the installation angle of the lug is equal to the parameter theta (θ) in the small elliptic parametric equation.
From Equation ( 4): In the above equation, the value range of x is [140,345]; α = 30 • ; R = 638 mm; The width of the lug along the axis of the iron wheel is 205 mm (see Figure 2b).The variable x in Equation ( 7) is replaced by the lug width parameter w (range: 0 to 205): According to Equation ( 8), the relationship curve between the installation angle of the lug and the width "w" along the axis of the wheel can be drawn, as shown in Figure 4.
Agriculture 2023, 13, x FOR PEER REVIEW 6 of 21 plane 4; angle gamma (γ): mounting angle of iron wheel lug faceplate; that is the angle between line 2 and line 3 at vertex point p2.
In Figure 3, the obliquely upward direction vector of line 2 obtained from Equation ( 4) is as follows: In the coordinate system xyz, line 3 is parallel to the z-axis, so the angle between line 2 and line 3 is the angle between the vector and the z-axis.
That is, the installation angle of the lug is equal to the parameter theta (θ) in the small elliptic parametric equation.From Equation ( 4): In the above equation, the value range of x is [140,345]; α = 30°; R = 638 mm; The width of the lug along the axis of the iron wheel is 205 mm (see Figure 2b).The variable x in Equation ( 7) is replaced by the lug width parameter w (range: 0 to 205): According to Equation ( 8), the relationship curve between the installation angle of the lug and the width "w" along the axis of the wheel can be drawn, as shown in Figure 4. To verify the correctness of Equation ( 8), the values of the installation angle "γ" at different values of w are obtained by a ProE three-dimensional modeling test, and the results are compared with those calculated by Equation ( 8).The results are shown in Table 1.
Comparisons between measured values of the three-dimensional model of installation angle gamma and calculated values of theoretical Equation ( 8) at different values of w.Because of the ProE 5.0 version of the three-dimensional software, its measurement accuracy is four digits after the decimal point.From Table 1, we can see that Equation (8) for calculating the installation angle γ is correct.

Driving Wheel Rotational and Translational Motion Relationship
It is assumed that when the tractor is driven on hard non-deformable or soft deformable terrain, the vehicle body maintains its movement in the horizontal plane.The vertical position is not affected by the rotation of the driving wheel, which is the rigid lugged wheel.Therefore, the driving wheel axis maintains a horizontal plane of movement.The radius R o represents the rolling radius of the wheel on hard un-deformable terrain; on soft, deformable terrain, when lugs penetrate the soil and the rim is not submerged, the rolling radius of the wheel is R, which is the radius of the rim.If the rim is submerged, R 1 represents the rolling radius of the wheel.
In Figure 5, the point "O" represents the wheel center or axle point of the rear axle of the tractor that moves along the horizontal direction with velocity "v" when the driving wheel rotates with an angular velocity "ω".This is pure rotational movement along imaginary plane 1, which is the ground surface.A perpendicular line passes through the axle point O, and the point of intersection with plane 1 is the instantaneous axis of rotation.
Agriculture 2023, 13, x FOR PEER REVIEW 8 of 21 over the terrain surface during the driving process is called travel reduction  and can be calculated by the following equation: The rolling radius of the driving wheel varies with field operational conditions, terrain properties, wheel structural parameters, applied normal and drawbar loads, and several other factors.Therefore, in practice, it is inconvenient to use the actual travel reduction for theoretical computations and analysis.Therefore, the maximum radius  of the rigid lugged wheel is used for the calculation of travel reduction by using Equation (12).
The rolling radius of the tractor wheel changes with the applied drawbar or traction load, which can be determined during the field-driving test.If the tractor operated at wheel rpm  under a certain applied load and traveled a certain distance s, the rotating radius of the wheel could be calculated by the following equation: The tractor wheel radius,  may vary according to the terrain and operational conditions; it is used for kinematic and dynamic analysis and determination and calculation of tractor traction performance parameters in field operations with the applied drawbar.

Motion Trajectory of the Lug Tip of the Rigid Lugged Wheel
While the driving wheel "rigid lugged wheel" rotates on the loose terrain/ field soil surface, the lugs of the wheel continuously pass through the soil.This process forms the cut marks in the soil that reflect the interaction between the wheel and the terrain soil.Therefore, it is important to analyze the wheel-soil interaction relative to the motion of the wheel to determine the stresses of the wheel lug on soil and the soil cutting and shearing process.The theoretical shape and size of the lug cut hole are determined by the motion trajectory, and they change with the structural parameters of the wheel lug, slip rate, i.e., travel reduction of the driving wheel, and height of lug penetration/wheel submergence.
When the rigid lugged wheel is driven on a flat, non-deformable surface, the terrain To produce the driving force through wheel-terrain interaction to move the tractor on the soft, deformable soil, the wheel rotates under the action of applied torque and normal load on the wheel.In this process, some slippage happens between the wheel and the terrain due to the shearing and sliding of surface soil.Ideally, the driving wheel rotates without slipping or sliding.
Figure 5 shows that the actual rolling motion of the driving wheel is in accordance with plan 3. The velocity v t is the theoretical velocity of the point "O", and "R o " is the maximum rolling radius, with "ω" as the angular velocity.Equation ( 9) can be used to calculate the theoretical velocity.
During the driving of a tractor for practical applications such as soil tillage, a certain amount of slip is always produced by the driving wheel, even without towed load conditions.Therefore, the actual velocity "v" of the tractor is less than the theoretical velocity "v t ".Equation ( 10) is used to calculate the actual velocity.
where, R 1 is the equivalent radius or actual rolling radius of the wheel in operational conditions.The loss of velocity caused by the slipping or shearing of the driving wheel over the terrain surface during the driving process is called travel reduction T r and can be calculated by the following equation: The rolling radius of the driving wheel varies with field operational conditions, terrain properties, wheel structural parameters, applied normal and drawbar loads, and several other factors.Therefore, in practice, it is inconvenient to use the actual travel reduction for theoretical computations and analysis.Therefore, the maximum radius R o of the rigid lugged wheel is used for the calculation of travel reduction by using Equation (12).
The rolling radius of the tractor wheel changes with the applied drawbar or traction load, which can be determined during the field-driving test.If the tractor operated at wheel rpm n under a certain applied load and traveled a certain distance s, the rotating radius of the wheel could be calculated by the following equation: The tractor wheel radius, R 1 may vary according to the terrain and operational conditions; it is used for kinematic and dynamic analysis and determination and calculation of tractor traction performance parameters in field operations with the applied drawbar.

Motion Trajectory of the Lug Tip of the Rigid Lugged Wheel
While the driving wheel "rigid lugged wheel" rotates on the loose terrain/ field soil surface, the lugs of the wheel continuously pass through the soil.This process forms the cut marks in the soil that reflect the interaction between the wheel and the terrain soil.Therefore, it is important to analyze the wheel-soil interaction relative to the motion of the wheel to determine the stresses of the wheel lug on soil and the soil cutting and shearing process.The theoretical shape and size of the lug cut hole are determined by the motion trajectory, and they change with the structural parameters of the wheel lug, slip rate, i.e., travel reduction of the driving wheel, and height of lug penetration/wheel submergence.
When the rigid lugged wheel is driven on a flat, non-deformable surface, the terrain surface contact points are located on the tips of the lugs, which are the highest points on the wheel's outer radius along the vertical axis.Therefore, kinematic analysis of the rigid lugged wheel is equivalent to finding the rotational and translational conditions of the contact points.A diagram representing the rotational conditions of the rigid lugged wheel with axes and plans of motion is shown in Figure 5, where point O represents the center of the wheel or wheel rotation axis that is at a height from the plan of motion (line 3) equivalent to the radius R o .A point M located on the tip of the lug in rest condition at time t = 0 has coordinates X M = 0, which means the horizontal displacement is 0, and Y M = 2R, which is the height of point M from the ground surface equal to the wheel diameter.
Assuming the wheel begins to rotate and moves with a uniform angular velocity ω and a linear velocity v.After a time t, the point M turns into M .By subtending an angle θ (Figure 5), the coordinates of the point M with respect to the ground-surface contact point reference frame can be determined by the following equations: The above equations are for ideal motion without slipping on a hard plain surface.However, in real conditions, in which wheels slip and slide on the terrain surface, the slip factor should be included, and the maximum radius, R o , should be replaced with the equivalent radius, R 1 , that is defined in Equation (13).Therefore, the equations to determine the position of the point are as follows when (θ = ωt): From the above equations, it is found that the motion trajectory of point M is associated with the rolling radius of the rigid lugged wheel and the wheel travel reduction rate.The motion trails of point M with different travel reduction rates are shown in Figure 6, and with different lug penetration or wheel submergence in Figure 7.The vertical axis (y-axis) represents the vertical position of the point, and the horizontal axis (x-axis) represents the horizontal position of the point at any instant of time while the wheel is rolling.The curves 1, 2, 3, and 4 in Figure 6 represent the trajectories of point M with a travel reduction rate of 0.0, 0.1, 0.25, and 0.4, respectively.Curve 1 is the trajectory of a point with ideal motion without slip, which means no travel reduction, and this pattern of motion is called cycloidal motion.Curves 2, 3, and 4 show the reduction in horizontal distance traveled per rotation of the wheel to calculate the travel reduction rate.As the travel reduction rate increases due to slipping and shearing between the surfaces, it produces more soil breakage and damage and reduces the distance traveled.Curves 1, 2, 3, and 4 in Figure 7 represent the trajectories of point M with lug penetration and wheel submergence of 0, 75, 150, and 165 mm, respectively.The effect of lug penetration and wheel submergence on rolling motion and horizontal distance traveled is obvious.
Agriculture 2023, 13, x FOR PEER REVIEW 9 of 21 of the wheel or wheel rotation axis that is at a height from the plan of motion (line 3) equivalent to the radius  .A point  located on the tip of the lug in rest condition at time  = 0 has coordinates  = 0, which means the horizontal displacement is 0, and  = 2, which is the height of point  from the ground surface equal to the wheel diameter.
Assuming the wheel begins to rotate and moves with a uniform angular velocity  and a linear velocity .After a time , the point  turns into  .By subtending an angle  (Figure 5), the coordinates of the point  with respect to the ground-surface contact point reference frame can be determined by the following equations: The above equations are for ideal motion without slipping on a hard plain surface.However, in real conditions, in which wheels slip and slide on the terrain surface, the slip factor should be included, and the maximum radius,  , should be replaced with the equivalent radius,  , that is defined in Equation (13).Therefore, the equations to determine the position of the point are as follows when ( = ): ) From the above equations, it is found that the motion trajectory of point  is associated with the rolling radius of the rigid lugged wheel and the wheel travel reduction rate.The motion trails of point  with different travel reduction rates are shown in Figure 6, and with different lug penetration or wheel submergence in Figure 7.The vertical axis (yaxis) represents the vertical position of the point, and the horizontal axis (x-axis) represents the horizontal position of the point at any instant of time while the wheel is rolling.The curves 1, 2, 3, and 4 in Figure 6 represent the trajectories of point  with a travel reduction rate of 0.0, 0.1, 0.25, and 0.4, respectively.Curve 1 is the trajectory of a point with ideal motion without slip, which means no travel reduction, and this pattern of motion is called cycloidal motion.Curves 2, 3, and 4 show the reduction in horizontal distance traveled per rotation of the wheel to calculate the travel reduction rate.As the travel reduction rate increases due to slipping and shearing between the surfaces, it produces more soil breakage and damage and reduces the distance traveled.Curves 1, 2, 3, and 4 in Figure 7 represent the trajectories of point  with lug penetration and wheel submergence of 0, 75, 150, and 165 mm, respectively.The effect of lug penetration and wheel submergence on rolling motion and horizontal distance traveled is obvious.The instantaneous velocity  of point  at any instant of the time can be determined by the following equation for the center of the wheel reference frame when the wheel is rotating with uniform angular velocity  and translating with a uniform linear velocity : When the rigid lugged wheel undergoes pure rotation at a certain velocity, the relationship between the vertical displacement and velocity of the edge points of the wheel lugs is shown in Figure 8a.The vertical axes (y-axis) denote the vertical displacement and velocity of the point.The horizontal axis (x-axis) denotes the horizontal displacement.Curves 1 and 2 represent the vertical displacement and velocity of the point  respectively.When the rigid lugged wheel goes through pure rotation at a slip rate of 0.4 on the soft deformable terrain surface, the relationship between the vertical displacement and velocity of the edge points of the wheel lugs is shown in Figure 8b.Curves 1 and 2 represent the vertical displacement and velocity, respectively.The instantaneous velocity v M of point M at any instant of the time can be determined by the following equation for the center of the wheel reference frame when the wheel is rotating with uniform angular velocity ω and translating with a uniform linear velocity v: When the rigid lugged wheel undergoes pure rotation at a certain velocity, the relationship between the vertical displacement and velocity of the edge points of the wheel lugs is shown in Figure 8a.The vertical axes (y-axis) denote the vertical displacement and velocity of the point.The horizontal axis (x-axis) denotes the horizontal displacement.Curves 1 and 2 represent the vertical displacement and velocity of the point M respectively.When the rigid lugged wheel goes through pure rotation at a slip rate of 0.4 on the soft deformable terrain surface, the relationship between the vertical displacement and velocity of the edge points of the wheel lugs is shown in Figure 8b.Curves 1 and 2 represent the vertical displacement and velocity, respectively.The instantaneous velocity  of point  at any instant of the time can be determined by the following equation for the center of the wheel reference frame when the wheel is rotating with uniform angular velocity  and translating with a uniform linear velocity :  =   sin() +   () When the rigid lugged wheel undergoes pure rotation at a certain velocity, the relationship between the vertical displacement and velocity of the edge points of the wheel lugs is shown in Figure 8a.The vertical axes (y-axis) denote the vertical displacement and velocity of the point.The horizontal axis (x-axis) denotes the horizontal displacement.Curves 1 and 2 represent the vertical displacement and velocity of the point  respectively.When the rigid lugged wheel goes through pure rotation at a slip rate of 0.4 on the soft deformable terrain surface, the relationship between the vertical displacement and velocity of the edge points of the wheel lugs is shown in Figure 8b.Curves 1 and 2 represent the vertical displacement and velocity, respectively.A comparative analysis of Figure 8a,b indicates that when the drive wheel slips, the displacement and velocity of the wheel lug tips (contact points) are repeated in the space.Thus, the vertical displacement and velocity of the wheel occur at lug tips or edge points, and the motion condition that is produced is more complex than pure rotation.
The impacts of the soil on the wheel during the wheel-soil interaction process produce power losses.The motion or rolling resistance is produced by the wheel penetration/sinkage and the soil reaction force, which is generated as the wheel lug penetrates and breaks/excavates the soil.The wheel lug impacting the soil backward also produces power losses, but it is necessary to obtain a wheel thrust or drive force from the soil.

Dynamic Analysis
In the process of cultivation or tillage, the tractor's forward speed is v m , along the positive direction of the y-axis.The angular speed of the iron wheel is ω, along the negative direction of the x-axis.Then, there is a rotation and translation transformation relationship between the coordinate system x 1 y 1 z 1 fixed on the iron wheel at the time 't' and the fixed coordinate system xyz at the beginning, and the rotation transformation relationship is shown in Figure 9.
Agriculture 2023, 13, x FOR PEER REVIEW 11 of 21 A comparative analysis of Figure 8a,b indicates that when the drive wheel slips, the displacement and velocity of the wheel lug tips (contact points) are repeated in the space.Thus, the vertical displacement and velocity of the wheel occur at lug tips or edge points, and the motion condition that is produced is more complex than pure rotation.
The impacts of the soil on the wheel during the wheel-soil interaction process produce power losses.The motion or rolling resistance is produced by the wheel penetration/sinkage and the soil reaction force, which is generated as the wheel lug penetrates and breaks/excavates the soil.The wheel lug impacting the soil backward also produces power losses, but it is necessary to obtain a wheel thrust or drive force from the soil.

The Normal Vector of the Lug Working Surface
In the xoyozo coordinate system, the normal vector of the working surface of the lug is ( ⃗ =  ⃗), by combining Equation (1), the normal vector of the lug working surface in the xyz coordinate system is: where: where: In the process of cultivation or tillage, the tractor's forward speed is  , along the positive direction of the y-axis.The angular speed of the iron wheel is ω, along the negative direction of the x-axis.Then, there is a rotation and translation transformation relationship between the coordinate system x1y1z1 fixed on the iron wheel at the time 't' and the fixed coordinate system xyz at the beginning, and the rotation transformation relationship is shown in Figure 9.  From Figure 9, following equation is derived: where Considering the rotation and translation relationship between the coordinate systems x 1 y 1 z 1 and xyz, there is: By adding Equation (19) to Equation ( 20), the expression of the normal direction vector of the lug working surface in the fixed coordinate system xyz is obtained as follows:

Driving Force Analysis of the Lug
Most of the upland farming areas in China cultivated two crops a year, so farmers tilled the soil twice a year with rotary or vertical tillage implements.Considering the physical characteristics of farmland soil, we can characterize it as an isotropic and completely plastic material.This study assumes that the depth of surface soil compacted by the iron wheel due to submergence of the iron wheel rim of radius R = 638 mm during cultivation is h = 15 mm, as shown in Figure 10.According to this figure: Agriculture 2023, 13, x FOR PEER REVIEW 13 of 21  By the above equation, the value of β is: As the width of the iron wheel B = 440 mm, the wheel soil contact area of the iron wheel rim is: It results in: S o = 60992.75144mm 2 .
Equation ( 4) was used to determine the corresponding position of two points A and B in Figure 2b, where x A = 345 and x B = 140 were on the lug in the fixed coordinate system xyz: By taking x A = 345 and x B = 140 into the above equation, it results in θ A = 18.19 • and θ B = 7.28 • , respectively, and by putting them in Equation ( 21), the trajectory equations of points A and B are obtained in return as: That resulted in the following equation: and Based on ProE three-dimensional software modeling and measurement, the area of the lug working surface shown in Figure 2b, is S = 34,343.2mm 2 .The force or soil thrust on the lug's working surface is related to the reaction force of the lug on the soil.Regardless of the friction force between the soil and the lug surface, it is assumed that the normal pressure of the soil on the driving or working surface of the lug is p l , the unit of which is N•m −2 .By using Equation (22), the positive pressure of soil, i.e., soil thrust, on the working surface of the lug is as follows: In the equation: S = area of lug working surface (34,343.2mm 2 ) it varies with the height of penetration of the lug, p l to be determined represents the compressive or normal stress of soil on the lug working surface, ω =1.28 rad s −1 , and 2.37 s ≥ t ≥ 2.18 s.Equation (27) shows that the axial component of the soil thrust acting on the lug, i.e., the x-axis component, is a constant, i.e., F x = S.p l .sin α = 0.5Sp l Therefore, the axial component accounts for 50% of the total force.The component of soil thrust acting on the lug working surface that helps to move the whole machine in the forward direction is the y-axis component, which is to say, the driving force provided by the interaction of the wheel lug and soil is: It can be seen that the lug's driving force is proportional to the angle α the smaller the angle, the greater the driving force.In addition, the driving force of the lug is related to the relative alignment of the working surface in the soil.In the beginning, the driving force is the smallest when the lug starts to penetrate and gradually increases until it is completely submerged in the soil.The fluctuation range of the driving force is 81.52% to 86.17%.
The vertical component of soil thrust acting on the lug is as follows: It can be seen that the vertical force decreases from 29.23% to 8.67%.The smaller vertical force helps to reduce the compaction effect of the whole machine on the soil [19,20].

Soil Thrust-Driving Force Relationship Analysis of Lugged Wheel
Predicting the normal pressure p of soil applied to the wheel contact surface or (p l = p) applied to the lug working surface is critical for analyzing the tractor drive wheel's traction performance.The principles used to predict the soil mechanical behavior and response under wheel loading are called the pressure-sinkage relationship and shear failure criteria.The pressure-sinkage relationship suggested by scientists who worked on soil mechanics is presented in Equation (30) [21].
This pressure sinkage relationship states that the average pressure applied to a plate on a uniform terrain is dependent on an exponential function of sinkage 'z', in which the exponent 'n' represents a soil property reliant on the size of the wheel and 'k' is the proportionality constant.A more explanatory and descriptive form of the pressure sinkage relationship of vehicle wheels is presented in Equation (31), which is a more widely relied-upon expression and is used in the practice of terramechanics and vehicle traction mechanics [21].
The proportionality constant 'k' has been revised to justify the dimensions of the contact area.The letter b represents the wheel width, and k 1 and k 2 remain parameters of soil sinkage.The soil parameters n, k 1 , and k 2 are empirical and have no intrinsic physical meaning.The parameters can be determined by pressure-sinkage data from field experiments.Although the pressure sinkage relationship is a semi-empirical expression, it is most important in terramechanics and tractor traction mechanics as the shear strength of the soil is mainly related to normal load (ground pressure).Additionally, the soil thrust developed under the drive wheel or wheel drive force and motion resistance is relative to sinkage (i.e., the relationship between pressure and sinkage).
The shear stress-displacement relationship, or share failure expression for maximum share stress and applied normal pressure, can be expressed by the Mohr-Coulomb equation that adequately describes the relation in many cases [17].
To describe the shared stress-displacement relationship for a deformable terrain, an expression was suggested by [22], in which the researcher [23] proposed a modification for loose plastic soil that is presented in the following equation and widely adopted in practice.
Therefore, Equation (32) has the following form: In the above expression, τ is the shear stress, j is the soil shearing displacement, K is the shear modulus, which can be determined by fitting the shear stress-shear displacement equation [23] to each data set, c is the soil cohesion, and φ is internal friction angle, which can be determined by the Mohr-Coulomb failure criteria.
The share displacement that is developed along the contact patch of the rigid wheel can be determined by analysis of the slip velocity v j on the wheel-terrain interface, as shown in Figure 11a [17,24,25].The tangential component of the absolute velocity v of a certain point on the circumference of the wheel relative to the terrain is described as slip velocity v j [17].The magnitude of the slip velocity v j of a point on the rigid wheel circumference defined by the angle β can, therefore, be expressed as below [24].
Agriculture 2023, 13, x FOR PEER REVIEW 16 of 21 is the travel reduction due to the wheel slip, and p is the normal pressure on the wheelterrain contact surface that is equal to the normal load of the wheel on the contact surface.As the front contact angle of the rigid lugged wheel is determined (i.e., 0 −  ) (Figure 11) [17], the rigid lugged wheel thrust force or drive force,  , can be determined by integrating the stresses over the wheel-terrain contact surface.
))/  (38) From the above equation, it is found that the driving force of a rigid lugged wheel is related directly to wheel width, radius, and sinkage and inversely to travel reduction rate.
When the tractor is operating on soft terrain and the majority of the rigid wheel's lugs are submerged in the terrain, the terrain material/soil fills the spaces between the lugs.In this case, the rigid lugged wheel acts like a smooth rigid wheel (Figure 11), so Equation (38), can be used for the thrust/drive force analysis of a rigid lugged wheel.However, over medium to hard terrain, the entire wheel load may be supported on the tips or faces of lugs, and there is minimum or no contact of the carcass (rim surface) with the terrain.In this case, the thrust-sinkage analysis, already done for the single lug, Equation (28), can be used by integrating the stresses over the whole contact area of lugs or the number of lugs by considering ( = ) as follows: The parameter j, which is the shared displacement of a point along the contact surface for rigid wheels and loose soil, can be determined by the following equation: According to soil shear theory and wheel-terrain interaction, the rigid lugged wheel has two main soil reaction forces along the y-axis: (1) F s due to the shear stress τ between the wheel bottom surface and soil surface, and (2) F l due to the normal terrain pressure p on the lug working surface.
The thrust-slip relationship for a rigid wheel, Equation (37), can be determined by using the shear stress-shear displacement relationship, Equation (34), by considering uniformly distributed normal pressure and the flat wheel-terrain contact patch of contact length l and contact area A [17].
where c, φ, and K are the cohesion, angle of shearing resistance or internal friction, and shear deformation parameters of the terrain material, "the surface soil", respectively.T r is the travel reduction due to the wheel slip, and p is the normal pressure on the wheel-terrain contact surface that is equal to the normal load of the wheel on the contact surface.
As the front contact angle of the rigid lugged wheel is determined (i.e., 0 − β o ) (Figure 11) [17], the rigid lugged wheel thrust force or drive force, F s , can be determined by integrating the stresses over the wheel-terrain contact surface.
From the above equation, it is found that the driving force of a rigid lugged wheel is related directly to wheel width, radius, and sinkage and inversely to travel reduction rate.
When the tractor is operating on soft terrain and the majority of the rigid wheel's lugs are submerged in the terrain, the terrain material/soil fills the spaces between the lugs.In this case, the rigid lugged wheel acts like a smooth rigid wheel (Figure 11), so Equation (38), can be used for the thrust/drive force analysis of a rigid lugged wheel.However, over medium to hard terrain, the entire wheel load may be supported on the tips or faces of lugs, and there is minimum or no contact of the carcass (rim surface) with the terrain.In this case, the thrust-sinkage analysis, already done for the single lug, Equation (28), can be used by integrating the stresses over the whole contact area of lugs or the number of lugs by considering (ωt = β) as follows: The net driving force, F d , produced by wheel-terrain interaction can be found by: Rolling resistance or motion resistance M r to the wheel of a moving tractor on deformable terrain could be produced due to the penetration of lugs and sinkage of the wheel.The power consumed in the motion resistance of the tire wheel is equal to the work conducted on soil compaction and rut formation, but in the case of high lugs rigid wheel, work is expended on the breaking and loosening of soil by the penetration of sharp-edged lugs.
The radial pressure σ applied on the wheel due to wheel-terrain interaction (Figure 11b) is equal to the normal pressure p below a sinkage plate at the same sinkage depth z in a pressure sinkage test; therefore, the following expression is developed.
By combining the pressure-sinkage relationship with the above equation and including the effect of travel reduction rate, the motion resistance of a rigid wheel can be determined as follows: where z o represents the wheel's maximum sinkage, as shown in Figure 11b.The above analysis shows that the motion resistance of a rigid lugged wheel depends on the wheel width, wheel radius, and sinkage depth; a higher travel reduction rate also causes higher motion resistance.
From the dynamic analysis of the rigid lugged wheel, it is found that the forces on the wheel include shear and normal stresses on the wheel rim and normal pressure on the wheel lug working surface.The normal pressure on the lug working surface results in a higher driving force from the rigid lugged wheel that is independent of the wheel axle's normal load.The motion resistance of a rigid lugged wheel could be higher due to the penetration of lugs into the soil, but this will help to break and loosen the soil and get compensated by the lower soil tillage draft.In this way, the rigid lugged wheel will help to consume tractor power efficiently and avoid soil compaction in soil tillage operations.

Model Evaluation (FEM Analysis) for Structural Stability and Optimization
The designed model of the rigid lugged wheel was analyzed in the ANSYS workbench platform for the structural optimization of a new wheel.The behavior of the structural components under applied forces on the designed model was analyzed on assessing stresses, strains, and deformations resulting from the analysis of applied loads.The purpose of this simulation analysis was to check the material suitability and design parameters to select the best-fit material and optimum parameters, e.g., the thickness of the steel sheet used for the wheel structure with the minimum wheel structure weight.

Material Properties (Engineering Data)
The material was selected based on design criteria, working conditions, and operational requirements.It was a Manganese Iron alloy called high Carbon Manganese steel (Steel No. 15 Mn) sheets that contain Carbon (0.12-0.18),Silicon (0.17-0.37), and Manganese (0.70-1.00) percent as chemical mass fractions.The engineering properties data of the material are as follows: the density of 15Mn steel is 7.8 × 10 3 kg m −3 , the modulus of elasticity is 200-220 G Pa, and the Poisson's ratio is 0.25-0.30.In the field, tractor wheels are constantly stressed and, in extreme cases, face high-intensity impact loads.The material used for the fabrication of the newly designed modern type rigid lugged wheel was chosen considering all possible loading conditions so that it can safely withstand applied stresses and all types of loads.

Static Structural Analysis
The structure of the rigid lugged wheel under static loads was analyzed for structural stability and optimization using ANSYS Workbench mechanical simulation software.ProE's (Creo Elements) 2D and 3D engineering design platform was used to generate 3D wheel model geometry, and the entire wheel structure assembly was saved in "igs" file format, which is compatible with and importable into ANSYS Workbench.In ANSYS R17 Workbench, the static structural project solver was used for analysis, engineering data was modified according to test specifications and analysis requirements, and the '.igs' geometry file was imported into the project schematic.After making some corrections and modifications, the final structure model was prepared, and the mesh was applied in Figure 12a,b.The meshing, loads, supports, and displacement setup were defined, and resulting behavior parameters, e.g., stress-strain and deformations were selected.The project was analyzed by the problem solver option in the design modeler of ANSYS.Results were generated in graphical form (color schematic), as presented in Figure 12c-h which explained the effect of applied loads on the wheel body.
the tractor and applied torque that acted on the lug working surface as the lugs of the wheel were in direct contact with the soil.The normal component of reaction forces to the lugworking surface provides the drive and traction force to the tractor.Therefore, the total force applied to the lugs would be equal to the maximum drawbar force of the tractor.To avoid the failure risk of the wheel structure in adverse field conditions, 20% of the total force was included as a safety factor.The reaction force, or soil thrust, developed on the wheel against the normal weight of the tractor and applied torque that acted on the lug working surface as the lugs of the wheel were in direct contact with the soil.The normal component of reaction forces to the lug-working surface provides the drive and traction force to the tractor.Therefore, the total force applied to the lugs would be equal to the maximum drawbar force of the tractor.To avoid the failure risk of the wheel structure in adverse field conditions, 20% of the total force was included as a safety factor.
In terms of the contact surface with the wheel, two different field conditions were considered.The first one was soft, easily penetrable soil where the whole lug can penetrate easily; in this case, a total of six lugs were considered to be in contact with the soil.In another case, the soil was hard enough that it allowed very little penetration, and only two lugs were in contact with the soil surface.The total force was divided accordingly in both cases; when force was divided among six lugs, there was less force per lug, so the stresses and deformations were less (Figure 12c-e) compared to another case where total force was applied only on two lugs, where stresses and deformations were higher (Figure 12f-h).
In both cases, the overall effects of loads and resulting stresses and deformations were very low and there was no situation was observed that affect the structural stability.Therefore, it is recommendable to suggest a reduced thickness of the steel sheet that might help to reduce the weight of the wheel, and ultimately, the total weight of the tractor helps to get optimum traction with minimum damage to the soil and the environment.

Wheel Prototype Development and Driving Tests
The designed rigid lugged wheel model was converted to a prototype after the structure stability analysis in the ANSYS workbench and verification of the designed parameters.Two rigid lugged wheel full-body structures were fabricated and developed according to the size and operational requirements of the specified tractor for this study.At a local industrial manufacturing factory, wheels were developed using a specific material (15Mn steel).
On-road drive and field tests were carried out to verify the design parameters, theoretical analysis, and traction performance (Figure 13), which represented the significant improvement in the wheeled tractor's traction performance with rigid lugged wheels as drive wheels compared to conventional tire wheels.Traction and wheel-soil interaction performance of a rigid lugged wheel is already reported [2,4], and further studies and analyses are under process to verify the effectiveness of the rigid lugged wheel for off-road drive operations of the conventional wheeled tractor.In terms of the contact surface with the wheel, two different field conditions were considered.The first one was soft, easily penetrable soil where the whole lug can penetrate easily; in this case, a total of six lugs were considered to be in contact with the soil.In another case, the soil was hard enough that it allowed very little penetration, and only two lugs were in contact with the soil surface.The total force was divided accordingly in both cases; when force was divided among six lugs, there was less force per lug, so the stresses and deformations were less (Figure 12c-e) compared to another case where total force was applied only on two lugs, where stresses and deformations were higher (Figure 12f-h).
In both cases, the overall effects of loads and resulting stresses and deformations were very low and there was no situation was observed that affect the structural stability.Therefore, it is recommendable to suggest a reduced thickness of the steel sheet that might help to reduce the weight of the wheel, and ultimately, the total weight of the tractor helps to get optimum traction with minimum damage to the soil and the environment.

Wheel Prototype Development and Driving Tests
The designed rigid lugged wheel model was converted to a prototype after the structure stability analysis in the ANSYS workbench and verification of the designed parameters.Two rigid lugged wheel full-body structures were fabricated and developed according to the size and operational requirements of the specified tractor for this study.At a local industrial manufacturing factory, wheels were developed using a specific material (15Mn steel).On-road drive and field tests were carried out to verify the design parameters, theoretical analysis, and traction performance (Figure 13), which represented the significant improvement in the wheeled tractor's traction performance with rigid lugged wheels as drive wheels compared to conventional tire wheels.Traction and wheel-soil interaction performance of a rigid lugged wheel is already reported [2,4], and further studies and analyses are under process to verify the effectiveness of the rigid lugged wheel for offroad drive operations of the conventional wheeled tractor.

Conclusions
Traction demand from conventional agricultural wheeled tractors has increased in recent years for cultivation operations according to the requirement of intensive traction, improved trafficability, and drivability without increasing machine weight to avoid soil compaction and excessive energy consumption.A new design of a high lug rigid wheel is developed and tested as a traction device.The main component of the research object is the design of the lug working faceplate.The working principle of the rigid lugged wheel

Conclusions
Traction demand from conventional agricultural wheeled tractors has increased in recent years for cultivation operations according to the requirement of intensive traction, improved trafficability, and drivability without increasing machine weight to avoid soil compaction and excessive energy consumption.A new design of a high lug rigid wheel is developed and tested as a traction device.The main component of the research object is the design of the lug working faceplate.The working principle of the rigid lugged wheel is theoretically confirmed by kinematic and dynamic analysis.The relationships and analysis are elucidated for wheel design parameters, terrain soil characteristics, rolling motion, and drive force generation behavior.A computer simulation of the designed model by the finite element method is presented that verifies the structure's stability and provides a reference for optimization.The complete wheel structure is developed and tested on field conditions while the theoretical analysis/modeling and design parameters optimization are in process.The results of theoretical analysis and field tests are consistent, indicating that the rigid lugged wheel is well suited to soft deformable terrain, with improved traction and drivability.Theoretical analysis is performed for the newly designed rigid lugged wheel that provides in-depth theoretical modeling to assess and confirm its working principles.Insertion pressure, rather than shearing pressure, generates driving power in a rigid lugged wheel with a sharp edge and flat contact surface.In this way, it can produce higher traction, give A pre-tillage effect to the top soil, and reduces the tillage draught requirement.As confirmed by theoretical analysis, the component of thrust along the direction of driving ranged from 81.52% to 86.17%, which is actually the driving force and relates to lug faceplate alignment and lug penetration depth.Moreover, the vertical component of thrust reported was less than 30% and decreases further to less than 9% as the wheel displaces in the soil, which will avoid soil compaction under the wheel.To test lugged wheels' field performance and conduct a comparative analysis with tire wheels, traction, and wheel-soil interaction performance are already been reported [2,4], which provide the validation base for these theoretical results.Further research on the rigid lugged wheel is required to improve the wheel structure parameters and adaptability in various adverse field conditions.The results of this study could provide a trusted technical reference for the future design and construction of off-road agricultural vehicles, agricultural tractors, and particularly traction devices/drive wheels for field operations.

Patents
High-power tractor drive wheel.

Figure 1 .
Figure 1.An illustration of the complete modeled structure of the rigid lugged wheel.

Figure 1 .
Figure 1.An illustration of the complete modeled structure of the rigid lugged wheel.

Figure 2 .
Figure 2. Lug face plate design.(a) Design model illustration, (b) Ellipses of lug working faceplate

Figure 3 .
Figure 3. Schematic explanation diagram of working faceplate mounting angle (a) Design model, (b) Prototype "labeling is explained in paragraphs where cited".

Figure 4 .
Figure 4.The relationship between the installation angle of the lug and its axial width.

Figure 3 .
Figure 3. Schematic explanation diagram of working faceplate mounting angle (a) Design model, (b) Prototype "labeling is explained in paragraphs where cited".

Figure 3 .
Figure 3. Schematic explanation diagram of working faceplate mounting angle (a) Design model, (b) Prototype "labeling is explained in paragraphs where cited".

Figure 4 .
Figure 4.The relationship between the installation angle of the lug and its axial width.

Figure 4 .
Figure 4.The relationship between the installation angle of the lug and its axial width.

Figure 5 .
Figure 5.The rolling condition diagram of the rigid lugged wheel.

Figure 5 .
Figure 5.The rolling condition diagram of the rigid lugged wheel.

Figure 6 .
Figure 6.Motion trail or trajectories of point M at different travel reduction rates, Corves 1, 2, 3 and 4 represent travel reduction rate of 0.0, 0.1, 0.25 and 0.4 respectively.

Figure 6 .
Figure 6.Motion trail or trajectories of point M at different travel reduction rates, Corves 1, 2, 3 and 4 represent travel reduction rate of 0.0, 0.1, 0.25 and 0.4 respectively.

Figure 8 .
Figure 8. Relationships of vertical displacement and velocity of lug edge point M at different travel reduction rates (a) 0, and (b) 0.4, and corves 1 and 2, represent vertical displacement and velocity respectively.

Figure 8 .
Figure 8. Relationships of vertical displacement and velocity of lug edge point M at different travel reduction rates (a) 0, and (b) 0.4, and corves 1 and 2, represent vertical displacement and velocity respectively.

Figure 8 .
Figure 8. Relationships of vertical displacement and velocity of lug edge point M at different travel reduction rates (a) 0, and (b) 0.4, and corves 1 and 2, represent vertical displacement and velocity respectively.

4 . 1 .
Determination of the Lugged Wheel Driving Force 4.1.1.The Normal Vector of the Lug Working Surface In the x o y o z o coordinate system, the normal vector of the working surface of the lug is ( → n o = → j ), by combining Equation (1), the normal vector of the lug working surface in the xyz coordinate system is: → n = A. → n o , where :

Figure 9 .
Figure 9. Coordinate rotation transformation in the tillage process.From Figure9, following equation is derived:

Figure 9 .
Figure 9. Coordinate rotation transformation in the tillage process.

Figure 10 ..
Figure 10.Schematic of the rigid lugged wheel in working condition: 1-Front untilled soil surface line; 2-Back compressed soil surface line; h: compressed depth of cultivated soil = 15 mm.Equation (4) was used to determine the corresponding position of two points A and B in Figure 2 (b), where x = 345 and x = 140 were on the lug in the fixed coordinate system xyz:    =  tan  .sin   sin   cos

Figure 10 .
Figure 10.Schematic of the rigid lugged wheel in working condition: 1-Front untilled soil surface line; 2-Back compressed soil surface line; h: compressed depth of cultivated soil = 15 mm.

Figure 11 .
Figure 11.The schematic force diagram of shear-displacement under the wheel terrain interface (a) slip velocity  and (b) with applied radial pressure .

Figure 11 .
Figure 11.The schematic force diagram of shear-displacement under the wheel terrain interface (a) slip velocity v and (b) with applied radial pressure σ.

Figure 12 .
Figure 12.Static structural simulation analysis process in ANSYS under applied loads and results of the analysis in different surface soil conditions.Figure 12. Static structural simulation analysis process in ANSYS under applied loads and results of the analysis in different surface soil conditions.

Figure 12 .
Figure 12.Static structural simulation analysis process in ANSYS under applied loads and results of the analysis in different surface soil conditions.Figure 12. Static structural simulation analysis process in ANSYS under applied loads and results of the analysis in different surface soil conditions.

Figure 13 .
Figure 13.Rigid lugged wheel on-road drive tests, (a) concrete, (b) asphalt, and (c) field drive and traction test.

Figure 13 .
Figure 13.Rigid lugged wheel on-road drive tests, (a) concrete, (b) asphalt, and (c) field drive and traction test.

Table 1 .
Measured and calculated values of the installation angle (γ) of the working faceplate in the 3D model.