Determination of the soil pressure distribution around a cone penetrometer

The objective of this paper was to investigate the pressure distribution around a cone penetrometer using a pressure sensing mat under laboratory conditions. The investigation was conducted under (1) constrained conditions using cylindrical split pipe molds and (2) unconstrained conditions using a soil box. These tests were conducted in Capay clay and Yolo loam soil containing two different moisture conditions and two compaction levels.In the constrained tests, a maximum radial pressure of 111 kPa was observed in the Capay clay soil with 3.4–4.3% d.b. moisture content and three blows of compaction (cone index value of 2040 kPa) when using the 41 mm diameter split pipe mold. These pressure levels decreased to 82 and 22 kPa, respectively, when 65 and 88 mm diameter molds were used. In both the Capay clay and Yolo loam tests, the average radial pressure and average cone index values showed similar trends.In the unconstrained tests, a maximum pressure of 9.0 kPa was observed in the Capay clay with 4.5% d.b. moisture content and three blows of compaction (cone index value of 550 kPa) at a horizontal distance of 25.4 mm from the vertical axis of the cone penetrometer and minimum pressure levels in the range of 0.2–0.3 kPa when the horizontal distance of the penetrometer was in the range of 56.8–66 mm. The pressure levels are much smaller than the ones obtained in the constrained tests and may suggest that the pressure distribution under field conditions is small at a distance of 25.4 mm or higher from the tip of the cone.The experimental data were statistically analyzed to identify significant factors. The results of the analysis for the constrained test indicated that the mold diameter and number of blows significantly increased the pressure readings within the soil mass. Increasing the mold diameter led to a decrease in the average radial pressure and increasing the number of blows contributed to an increase in the average radial pressure. In the unconstrained test, the average radial pressure distribution at a given point were significantly influenced by the horizontal distance of the point from the vertical axis passing through the center of the penetrometer shaft, soil type, and soil moisture content. Higher pressure values were obtained in the Capay clay tests compared to the Yolo loam tests. In all cases, the pressure levels were greater for the drier soil than for the moist soil.     

The use of soil compaction levels in the selection of field-safe sugarcane transport vehicles

Data was collected for single bundle and nucleus estate trailers aimed at selecting the trailer units that could safely travel in the sugarcane fields without causing detrimental soil compaction. The proportion of trailers carrying loads in excess of established safe axle loads was assessed. Over 60% single bundle trailers traveling in sandy loam and sandy clay loam fields were found not to induce detrimental soil compaction. Nucleus estate trailers, however, were sufficiently loaded to cause significant soil compaction. Working under soil moisture contents of 21.4–27.1% (dry basis), safe loads were found to be payloads of 64.9 and 46.1 kN carried by single bundle and nucleus estate trailers (respectively) on a single axle having two 10 ply 18.4×30 tires with an inflation pressure of 207 kPa.     

Modeling compaction in agricultural soils

Research was conducted to quantify the effect of tire variables (section width, diameter, inflation pressure); soil variables (soil moisture content, initial cone index, initial bulk density); and external variables (travel speed, axle load, number of tire passes) on soil compaction and to develop models to assess compaction in agricultural soils. Experiments were conducted in a laboratory soil bin at the Asian Institute of Technology in three soils, namely: clay soil (CS), silty clay loam soil (SCLS), and silty loam soil (SLS). A dimensional analysis technique was used to develop the compaction models. The axle load and the number of tire passes proved to be the most dominant factors which influenced compaction. Up to 13% increase in bulk density and cone index were observed when working at 3 kN axle load in a single pass using a 8.0–16 tire. Most of the compaction occurred during the first three passes of the tire. It was also found that the aspect ratio, tire inflation pressure and soil moisture content have significant effect on soil compaction. The initial cone index did not show significant effect. The compaction models provided good predictions even when tested with actual field data from previous studies. Thus, using the models, a decision support system could be developed which may be able to provide useful recommendations for appropriate soil management practices and solutions to site-specific compaction problems.     

Non-linear finite element analysis of cone penetration in layered sandy loam soil – Considering precompression stress state

Axisymmetric finite element (FE) method was developed to simulate cone penetration process in layered granular soil. The FE was modeled using ABAQUS/Explicit, a commercially available package. Soil was considered as a non-linear elastic plastic material which was modeled using variable elastic parameters of Young’s Modulus and Poisson’s ratio and Drucker–Prager criterion with yield stress dependent material hardening property. The material hardening parameters of the model were estimated from the USDA-ARS National Soil Dynamics Laboratory – Auburn University (NSDL-AU) soil compaction model. The stress–strain relationship in the NSDLAU compaction model was modified to account for the different soil moisture conditions and the influence of precompression stress states of the soil layers. A surface contact pair (‘slave-master’) algorithm in ABAQUS/Explicit was used to simulate the insertion of a rigid cone (RAX2 ABAQUS element) into deformable and layered soil medium (CAX4R ABAQUS element). The FE formulation was verified using cone penetration data collected on a soil chamber of Norfolk sandy loam soil which was prepared in two compaction treatments that varied in bulk density in the hardpan layer of (1) 1.64 Mg m⁻³ and (2) 1.71 Mg m⁻³. The FE model successfully simulated the trend of cone penetration in layered soils indicating the location of the sub-soil compacted (hardpan) layer and peak cone penetration resistance. Modification of the NSDL-AU model to account for the actual soil moisture content and inclusion of the influence of precompression stress into the strain behavior of the NSDL-AU model improved the performance of FE in predicting the peak cone penetration resistance. Modification of the NSDL-AU model resulted in an improvement of about 42% in the finite element-predicted soil cone penetration forces compared with the FE results that used the NSDL-AU ‘virgin’ model.     

Comparisons of different procedures of pre-compaction stress determination on weakly structured soils

Compaction is an important component of soil degradation. In this regard, the pre-compaction stress (σ_pc) concept is considered useful in mechanized agriculture nowadays. When the external forces exceed the internal strength (σ_pc) of soil, soil structure and soil physical quality will deteriorate. This concept was introduced at first for confined consolidation of non-structured, homogenized and saturated subsoils in civil engineering, though it is also suitable for agricultural conditions where the topsoil and subsoil are considered and both are often structured, heterogeneous and unsaturated. The best method for predicting σ_pc is by the plate sinkage test (PST) in the field, but it is expensive and time-consuming. This study was conducted to find an alternative laboratory method besides the confined compaction test (CCT) for predicting σ_pc. The CCT may not be a good method, especially at higher water contents, and for soils with low organic matter content, because of low sharpness of the critical region on the stress–strain curves. The study was performed on five soil types with a range of soil textures and organic matter content from central Iran using three loading types and three pF (i.e. Log [soil matric suction in cm]) values of 2.3, 2.7 and 2.9 with two replicates. Loading types consisted of CCT, the semi-confined compaction test (SCCT) and the kneading compaction test (KCT) at three maximum (or pre-compaction) stresses of 200, 400 and 600 kPa. The experiment was a completely randomized factorial design. The aim was to determine how accurately each loading type test could predict σ_pc of soil pre-compacted by one of the other methods. The applied combinations of CCT–SCCT, SCCT–CCT and KCT–CCT mean that the soil was pre-compacted by the first loading type and evaluated by the second one. The results showed that σ_pc and the sharpness of the σ_pc region were significantly affected by loading types as well as the soil conditions. The sharpest σ_pc region was observed in SCCT and the least sharp in CCT with the overall order being CCT–SCCT > SCCT–CCT > KCT–CCT. The sharpness of the σ_pc region was reduced for the soil samples with higher water content and coarser texture. Regardless of the soil and loading conditions, the prediction by SCCT was consistently more accurate than by CCT. The prediction of σ_pc by SCCT was more precise in comparison with CCT especially at higher stresses and soil water contents. However, the prediction of σ_pc by SCCT was very accurate at pF values of 2.7 and 2.9, and with low σ_pc values, when compared with the actual values of the σ_pc. For the clay soil at a pF value of 2.3, no pre-compaction region (i.e. zero σ_pc) could be determined by CCT at a maximum axial stress of 600 kPa. This was because of the incompressibility of soil water at this near-saturated soil condition at high stress. However, the sharpness of the critical region in SCCT was high enough to predict σ_pc satisfactorily. There was no significant difference between the combinations of SCCT–CCT and KCT–CCT in predicting σ_pc. The SCCT is a compromise method that lies between CCT and PST. SCCT has the advantage of using a limited and definite soil volume that can be modeled as a soil element. Marginal effects of disturbance caused by coring/sampling as well as pre-test sample preparation seem to have minor effects on the stress–strain curve determined by SCCT in comparison with CCT. Moreover, the soil volume needed for this test is the same as for CCT.     

Establishment of safe axle loads for sugarcane soils under varying moisture content

Studies were conducted for the establishment of safe axle loads for sugarcane hauling vehicles beyond which detrimental soil compaction would be induced. The treatments involved running a loaded test vehicle in field strips previously chosen at random. Safe loads were established by testing the level of significance of the difference in induced soil compaction between treated and non-treated sections. Working under soil moisture contents of 21.4–27.1% (dry basis), safe axle loads for two 18.4 × 30 tires were found to be 55.6 and 60.0 kN for sandy clay loam and sandy loam soils with initial dry bulk density about 1.434 g/cm³. These corresponded to ground contact pressures of 111 and 120 kPa, respectively.     

Measuring soil compaction and soil behavior under the tractor tire using strain transducer

Soil compaction can occur due to machine traffic and is an indicator of soil physical structure degradation. For this study 3 strain transducers with a maximum displacement of 5 cm were used to measure soil compaction under the rear tire of MF285 tractor. In first series of experiments, the effect of tractor traffic was investigated using displacement transducers and cylindrical cores. For the second series, only strain transducers were used to evaluate the effect of moisture levels of 11%, 16% and 22%, tractor velocities of 1, 3 and 5 km/h, and three depths of 20, 30 and 40 cm on soil compaction, and soil behavior during the compaction process was investigated. Results showed that no significant difference was found between the two methods of measuring the bulk density. The three main factors were significant on soil compaction at a probability level of 1%. The mutual binary effect of moisture and depth was significant at 1%, and the interaction of moisture, velocity, and depth were significant at 5%. The soil was compressed in the vertical direction and elongated in the lateral direction. In the longitudinal direction, the soil was initially compressed by the approaching tractor, then elongated, and ultimately compressed again.     

Effects of tillage systems on compaction and crop yield of Albic Luvisol in Croatia

This four-year experiment was conducted in north-west Slavonia (agricultural area of Croatia) to evaluate the effects of different tillage systems on compaction of silty loam soil (Albic Luvisol). The compared tillage systems were: (1) conventional tillage (CT), (2) conservation tillage (CM), (3) no-tillage system (NT), and the crop rotation was corn (Zea mays L.) – winter wheat (Triticum aestivum L.) – corn – winter wheat. For detecting the soil compaction, bulk density and penetration resistance were measured during the growing seasons. In all seasons and tillage systems, the bulk density and penetration resistance increased with depth and the greatest increase from surface to the deepest layer in average was observed at CT system. The bulk density and penetration resistance increased at all tillage systems during the experiment, but the greatest increase was also observed at CT system. The greatest bulk density (1.66 Mg m⁻³) and the greatest increase of 6.4% were observed at CT system in the layer 30–35 cm. In the first season, the bulk density was the greatest at NT system, but during the experiment the lowest average increase of 1.9% was observed at this system. The greatest penetration resistance of all measurements (5.9 MPa) was observed in the last season at CT system in depth of 40 cm. The lowest average increase of penetration resistance 11.4% was also observed at NT system. The highest yield of corn in the first season was achieved with CT system while in other seasons the highest yield of winter wheat and corn was achieved with CM system.     

Estimation of a three-dimensional tyre footprint using dynamic soil–tyre contact pressures

A method for estimating the three-dimensional (3D) footprint of a 16.9R38 pneumatic tyre was developed. The method was based on measured values of contact pressure at the soil–tyre interface and wheel contact length determined from the contact pressures and the depths and widths of ruts formed in the soil. The 3D footprint was investigated in an area of the field where the pressure sensors of the tyre passed in a soft clay soil. The tyre was instrumented with six miniature pressure sensors, three on the lug face and the remaining three on the under-tread region between two lugs. The instrumented tyre was run at a constant forward speed of 0.27 m/s and 23% slip on a soft soil, 0.48 MPa cone index, 25.6% d.b. moisture content for four wheel load and tyre pressure combination treatments. The 3D footprint assessment derived from soil–tyre interface stress used in this research is a unique methodology, which could precisely relate the trend profile of the 3D footprint to the measured rut depth. The tyre–soil interface contact pressure distributions results showed that as inflation pressure increased the soil strength increased significantly near the centre of the tyre as a compaction increase sensed with the cone penetrometer.     

Vehicular traffic effects on development and yield of corn (maize)

A 52 plot experiment was performed during the growing season of 1977 in a Ste. Rosalie clay soil, using a randomized complete-block design with 13 treatments of machinery traffic within each of four blocks. Three vehicle contact pressures, four numbers of tractor passes, and a control of zero traffic were used to relate the growth and yield variables to wheel traffic and the resulting soil compaction. The number of days required for (a) plant emergence, (b) tasselling, and (c) silking were minimum in moderately compacted plots. The plant growth rate monitored at 41, 54, 68, and 105 days from the seeding time was different from plot to plot. Growth models at different times of the season were derived in terms of the wheel traffic variables, and plant and ear moisture content variation with traffic treatment was investigated. Yield and ear yield increased with increases in machine contact pressure and passes reaching a peak around 500 kPa and dropped off for further increases. The reduction in yield was over 35% in some cases, suggesting that careful traffic planning is essential to obtain better production in agricultural fields. Prediction models were obtained for all the plant growth characteristics in terms of traffic variables. A relation for yield in terms of soil bulk density was established for a dry season.     

Cone penetration resistance equation as a function of the clay ratio, soil moisture content and specific weight

The cone penetrometer is a simple versatile device which is widely used to monitor the strength of a soil in terms of its resistance to the penetration of a standard cone. The soil penetration resistance is a function of soil moisture content, soil specific weight and soil type. The soil type is characterised by means of a clay ratio which is the ratio of the clay content of the soil to the content of silt and sand.Based on the classical bearing capacity theories for strip foundations, a general cone penetration resistance equation is developed to represent the variability of cohesion and friction angle by means of soil type and moisture content. The empirical relationship is shown to give an accurate prediction of the cone penetration resistance for a wide range of soils from a loamy sand to a heavy clay (clay ratios 0.10–1.60) and over a wide spectrum of soil moisture contents from 10 to 65% w/w.     

Effects of timber skidding on chemical characteristics of herbaceous cover, forest floor and topsoil on skidroad in an oak (Quercus petrea L.) forest

This paper examines the effects of timber harvesting by skidding on some soil properties (sand, silt, clay, pH, organic carbon, bulk density and compaction), herbaceous cover (unit mass) and forest floor (unit weight) properties. Also N (%), P, K, Na, Ca, Mg, Fe, Zn, Cu and Mn (ppm) were determined in all herbaceous cover, forest floor and two soil depth (0–5 cm and 5–10 cm) on skidroad of an oak (Quercus petrea L.) stand in Istanbul Belgrad Forest – Turkey. In this study, obtained results are; the forest floor and the herbaceous cover amount on the skidroad have been found considerably lower than undisturbed area. There were some crucial changes in the characteristics of the soil which has been investigated down to 10 cm depth. Soil bulk density was found quite high in the samples taken from the skidroad subject to compaction compared to the ones on the undisturbed area. Nevertheless, no important difference had been detected between the skidroad and the undisturbed area at both soil depths in terms of organic carbon contents. Moreover, the soil acidity (pH) values showed noteworthy differences in the analysis of soil samples taken from both soil depths on the skidroad and on the undisturbed area. Fe and Cu contents of herbaceous samples on skidroad were significantly higher than undisturbed area. Forest floor on skidroad had significantly higher K content, and significantly lower Zn, Mn and N content compared to undisturbed area. P, Fe, Zn and Mn contents were found significantly lower in 0–5 cm soil depth on skidroad than undisturbed area. In 5–10 cm soil depth, concentrations of N, P, Fe, Zn and Mn were significantly lower, while Mg and Cu contents were significantly higher than undisturbed area. Results indicate that long-term harvest using skidding techniques on these sites had adversely affected soil cation concentrations, physical soil conditions and mass of herbaceous cover and forest floor.     

Soil stress state orientation beneath a tire at various loads and inflation pressures

Stress state transducers (SSTs) were used to determine the orientation of the major principal stress, σ₁, in soil beneath the centeline of an 18.4R38 radial-ply R-1 drive tire operated at 10% slip. Two soils, a sandy loam and a clay loam, were each prepared twice to obtain two density profiles. One profile of each soil had a hardpan and the soil above the hardpan was loose. The soil in the second profile was loosely tilled. The stress state was determined at a depth of 358 mm in the sandy loam and 241 mm in the clay loam soil. The tire was operated at two dynamic loads (13.2 and 25.3 kN), each at two levels of inflation pressure (41 and 124 kPa). When the orientation of σ₁ was determined directly beneath the axle, the mean angles of tilt in the direction of travel ranged from 6 to 23 degrees from vertical. Inflation pressure did not significantly affect the angle when the dynamic load was 13.2 kN in the sandy loam soil, and neither inflation pressure nor dynamic load significantly affected the angle in the clay loam soil. When the dynamic load was 25.3 kN in the sandy loam soil, the orientation of the major principal stress determined directly beneath the axle was tilted significantly more in the direction of travel when the tire was at 41 kPa inflation pressure than when at 124 kPa. These changes in stress orientation demonstrate the importance of measuring the complete stress state in soil, rather than stresses along only one line of action. The changing orientation of σ₁ as the tire passes over the soil indicates the soil undergoes kneading and supports future investigation of the contribution of changes in stress orientation to soil compaction.     

Collapse failure in dry clay soils caused by tine implements

The relation between forces applied to the soil and the resultant soil reaction was studied in dry clay soils under a quasi-static condition. As a tine advanced in dry compact clay soils at 5.2% dry basis moisture content, masses of soil collapsed one by one in front of the tine. The horizontal and vertical components of soil resistance measured were cyclic and in phase, with distinct peak and trough values. The peak values and trough values indicated the soil stress conditions before and immediately after each failure occurred. The frequency of failure depended on the size of the tine. The magnitude of the peak values depended on level of compaction and trough values on density of collapsed mass. The paper presents the details of observations.     

Prediction of clay soil compaction

Field measurements were made of soil density and moisture patterns under different vehicle tire paths with varying external pressures and number of passes. In addition, laboratory index tests were performed to determine the compaction behaviour of the same soil. Using these results, a prediction equation of dry density in terms of applied pressure and moisture content was obtained for the clay soil. A previously developed equation for sandy soil was modified for the complete range of moisture contents encountered. Estimation of shear strength for the clay soil was made using plastic and liquid limits.     

Effect of wetting and drying on soil physical properties

Agricultural soils are subject to seasonal wetting and drying cycles. Effect of drying stress, as influenced by one cycle of wetting and drying, on physical properties of a clay–loam soil was investigated in the laboratory. The physical properties studied were soil bulk density, cone penetration resistance, shear strength, adhesion and aggregate size and stability. Three drying stress treatments were made by wetting air-dried soil of initial moisture content of 12% (on dry weight basis) to three different higher moisture contents, namely 27, 33 and 40%, and then drying each of them back to their original moisture content of 12%. Thus, the soil was subjected to three different degrees of drying stress. The results showed that the soil strength indicated by cone penetration resistance and cohesion, and soil aggregate size, increased with the degree of drying stress. However, the soil bulk density did not change significantly with the drying stress.     

Soil compaction and tires for harvesting and transporting sugarcane

Four tire types (A, block-shape tread; B, rib-shape tread; C, low-lug tread; D, high-lug tread) used to harvest and transport sugarcane were compared regarding the compaction induced to the soil. Tires were tested at three inflation pressures (207, 276, 345 kPa) and six loads ranging from 20 to 60 kN/tire. Track impressions were traced, and 576 areas were measured to find equations relating inflation pressure, load, contact surface and pressure. Contact surface increased with increasing load and decreasing inflation pressure; however, the contact pressure presented no defined pattern of variation, with tire types A and B generating lower contact pressure. The vertical stresses under the tires were measured and simulated with sensors and software developed at the Colombian Sugarcane Research Center (Cenicaña). Sensors were placed at 10, 30, 50 and 70 cm depth. Tire types A and B registered vertical stresses below 250 kPa at the surface. These two tires were better options to reduce soil compaction. The equations characterizing the tires were introduced into a program to simulate the vertical stress. Simulated and measured stresses were adjusted in an 87–92% range. Results indicate a good correlation between the tire equations, the vertical stress simulation and the vertical stress measurement.     

Traffic-soil-plant (maize) relations

Twenty-five treatments consisting of three vehicle contact pressures, 62, 41 and 31 kPa (0.63, 0.42, 0.32 kg/cm²), four numbers of tractor passes (1, 5, 10, 15,) before and after seeding groups, and a control of zero traffic were used to study the effect of soil compaction on corn plant root growth and distribution in a Ste. Rodalie clay soil. The average dry bulk density values for 0–20 cm depths measured during the season varied from a minimum of 0.89 g/cm³ to a maximum of 1.12 g/cm³ depending on the severity of the treatment. Root distribution maps were obtained for all the treatments by field measurements coupled with root washing methods. An average root density of 5.7 mg/g of soil in an uncompacted control plot was reduced to less than 2 mg/g in a plot with 15 passes of 0.63 kg/cm² contact pressure. Soil penetration resistance values in various plots were compared, and a statistical model was obtained in terms of the traffic treatments, soil moisture content and depth. Yield reductions and penetration resistance were compared to root distrubution density results.     

Performance prediction of animal drawn vehicle tyres in sand

Four animal drawn vehicle (ADV) tyres of 5.00–19, 6.00–19, 7.00–19 and 8.00–19 sizes were tested in sand under various but controlled conditions in an indoor soil bin. A tyre test carriage with four-bar parallel linkage was developed for accommodating a single wheel of different sizes. Performance tests were conducted at five levels of inflation pressure and load. The sand compaction level was varied in the range of 3.4–4.5 MPa/m and forward speed of the test carriage was maintained at 3.1 km/h. Performance of the tyres 7.00–19 and 8.00–19 was identical and offered less rolling resistance as compared to other tyres. However, their use in camel carts may not be recommended beyond the payload of 6 kN on single wheel with inflation pressure and sand compaction range of 172–379 kPa and 3.4 –4.5 MPa/m, respectively. Based on the experimental results, empirical models were developed to predict the performance of tyres. The accuracy of prediction of the developed empirical models was compared with that of existing semi-empirical approaches. Model with sand mobility number considered relatively simple and convenient to use in the field and yields reasonably good prediction for rolling resistance and sinkage.     

Modification of a mouldboard plough surface using arrays of polyethylene protuberances

Surface design modifications have recently exhibited a means of reducing soil-tool adhesion. The tribological characteristics of soil-burrowing animals were employed on tillage machinery to study the effect on adhesion. Considering the characteristics of dung beetles, ultra high molecular weight polyethylene (UHMW-PE) protuberances were mounted as embossed arrays on a mouldboard plough. To investigate a suitable geometry of such protuberances, five shapes were studied (flat, semi-spherical, semi-oblate, semi short-prolate and semi long-prolate) using a combination of base diameters and protrusion heights. The dimensionless height to diameter ratio (HDR) was used to characterize the geometry. Experiments were conducted to evaluate the resultant influence of various geometrical shapes and sizes of the protuberances (base diameter: 20–50 mm; protrusion height: 0–50 mm) on lowering the ploughing resistance of the mouldboard plough in Bangkok clay soil. A comparison was made between the modified and the conventional plough in dry (21.8% (d.b.)), sticky (37.2% (d.b.)), wet (49.1% (d.b.)) and flooded (64.3% (d.b.)) Bangkok clay soil at 1, 3 and 5 km/h forward speeds. Percent reduction in ploughing resistance of bionic mouldboard plough in these soil conditions with HDR = 0 was 1–6% in dry soil, 16–22% in sticky soil, 14–20% in wet soil and 8–12% in flooded soil. With HDR = 0.25 the ploughing resistance was reduced by 2–7% in dry soil, 18–36% in sticky soil, 17–33% in wet soil and 15–28% in flooded soil. Similarly with HDR = 0.5, it reduced by 10–16% in sticky soil, 6–17% in wet soil and 12–26% in flooded soil. Whereas, HDR > 0.5 increased the ploughing resistance by 7–29%.