1. Introduction
Near-surface soil compaction in earthworks is mainly performed by vibratory rollers, as they have a greater depth effect than static and oscillatory rollers and therefore allow the installation of greater layer thicknesses. The rotation of an unbalance mounted concentrically in the drum axis of the roller provides for its dynamic excitation and causes a predominantly vertically directed load transfer into the soil.
During the entire compaction process, the dynamically excited drum and the soil form an oscillating, interacting system with changing contact conditions. Depending on the tuning of operating or process parameters (travel speed, excitation frequency, amplitude, etc.), design parameters (ratio of static axle load to vibrating mass, eccentricity, mass of unbalance, etc.) and soil stiffness (grain size distribution, grain shape, mineral structure, water content, bedding density, loading history, etc.), the roller thereby operates in different modes of operation [
1]. The resulting characteristic motion behavior of the compaction device can be used to derive and calculate measurement values that allow for an assessment of the compaction state of the soil or its load-bearing capacity. Systems for continuous compaction control (ccc) and intelligent compaction (IC) [
1,
2,
3,
4,
5,
6] make use of this basic idea.
The dynamic motion behavior of the roller drum as a result of the vibratory excitation is a complex and highly nonlinear process. The occurring nonlinearities result primarily from the curved shape of the drum (geometric nonlinearity), from the plastic behavior of the soil (material nonlinearity) and in particular from the aforementioned interaction between the roller and soil (system nonlinearity due to changing contact conditions).
Compaction work usually starts in “Continuous Contact”, where there is permanent contact between the drum of the roller and the soil. During the loading phase, the contact area between the drum and soil surface increases steadily due to the curved drum shape, which in turn influences the reaction forces of the soil and changes the motion behavior of the drum. In the course of the unloading phase, the drum may be lifted off the ground, causing the two subsystems to move separately for a short time. This loss of contact occurs only after the dynamic force between the drum and the soil has assumed a maximum soil reaction force. The motion behavior changes significantly when the two subsystems are separated (change in mode of operation). Thus, the system reacts to the increasing compaction and the resulting increasing soil stiffness with additional nonlinearity (loss of contact). If a loss of contact occurs between the drum and the soil, it is well known that several modes of operation can occur, in addition to “partial uplift”, especially “double jump”, and possibly also “tumbling” and “chaos” [
1]. In this paper, the loss of contact is considered exclusively in the “partial uplift” mode of operation. However, all considerations can also be applied analogously to the “double jump” mode of operation.
In the last decades, numerous research studies have been conducted to better understand the interaction between the vibrating drum of a vibratory roller [
2,
4,
7,
8,
9] or an oscillating roller [
10,
11] and the soil to be compacted.
Lumped parameter models [
7,
8,
12] and cone models [
1,
3] have already been used in the literature to analyze the drum–soil mechanics but only under the assumption of a constant contact geometry between the drum of the vibrating roller and the soil surface.
Dynamic elasto-plastic finite element models have been developed (e.g., [
13,
14]) to investigate the interaction between a vibrating drum and the soil. Kenneally et al. [
15] already highlight the importance of the consideration of the variable contact geometry between the drum and soil in their FE analysis. Most recent studies primarily focus on the application of machine learning or artificial neural networks (ANNs) on intelligent compaction (e.g., [
16]).
This paper presents a comparatively simple semi-analytical vibratory roller–soil simulation model to describe the interaction system, characterized by its challenging contact conditions. Previous models consider, if at all, only a constant geometry of the contact area between the drum and soil. However, the width of the contact area increases during the loading phase when the drum penetrates the soil, making the drum contact width a variable parameter. The presented model is able to account for this variability of the drum contact width in the loading phase of vibratory roller compaction.
The semi-analytical model is part of an ongoing research project with the goal to develop a novel ICMV for vibratory rollers that is largely independent of machine (drum geometry and mass ratios) and process parameters (excitation frequency, excitation amplitude, and travel speed) and, at best, influenced solely by the stiffness of the soil.
First, the equations of motion for modeling the drum of the roller and the soil are derived separately by means of the subsystem technique. This is followed by the coupling of the two interacting subsystems via the formulation of contact conditions (contact problem). In this way, the three possible operating phases per revolution of the exciter—loading, unloading and loss of contact—can be mapped.
2. Modeling of the Roller Subsystem
The actual compaction tool is the cylindrically curved drum of the roller, in the axis of which an unbalance (exciter) is mounted on a drive shaft which rotates at a specific exciter frequency f. Due to the quadratic excitation, the drum, which is considered rigid, undergoes a circular translatory motion, whereby the soil is compacted by a predominantly vertically directed loading.
The drum is connected to the frame of the roller via rubber buffers, modeled as a Kelvin–Voigt element (a linear elastic spring
and a viscous dashpot
in parallel). The roller subsystem under consideration thus corresponds to a two-mass oscillator with the absolute displacement of the frame
and the absolute displacement of the drum
as position coordinates. The initial position of the entire system in static equilibrium is shown in
Figure 1 (left).
If the rubber buffers and frame are suitably designed, the frame and drum can be regarded as largely decoupled from vibration so that the energy of the dynamic excitation can theoretically be fully utilized as compaction energy [
17]. Under this condition, the frame can be considered to be at rest (
), which means that all time-dependent terms of the frame displacement
as well as the forces in the rubber buffers
and
can be neglected [
1]. The influence of the frame on the motion behavior of the drum is thus limited to its dead weight (
), where
is the mass of the frame. At the same time, the roller subsystem is now reduced to a single degree of freedom with the position coordinate
(see
Figure 1, right). From the application of the principle of linear momentum [
18],
where
and
are the masses of the drum and the eccentric, respectively, follows the equation of motion for the roller subsystem by means of the resulting soil contact force
between the drum and soil:
The soil contact force
in Equations (
1) and (
2) is defined positive for compression and comprises three parts: the inertia force
of the vibrating mass, the static load
of the dead weight and the force
of the dynamic excitation.
e denotes the eccentricity of the unbalance, and
is the angular frequency of the excitation. The product of the angular frequency
and time
t gives the rotation angle of the eccentric mass at time
t (see
Figure 1, right).
3. Modeling of the Soil Subsystem
The soil model presented in the following is not able to represent the processes that occur during dynamic compaction with the vibratory rollers, but it is very well able to describe the compaction state of the homogeneous soil after compaction, which is ultimately also used to assess the compaction quality. It is pointed out that for modeling the dynamic processes in the soil, the use of higher-order material models (e.g., hypoplasticity [
19]) is inevitable. However, the main focus of the presented soil model is on the determination of the resulting soil contact force as well as the governing motion quantities (accelerations, velocities, and displacements). As an essential simplification, only vertical movements on homogeneous soil are taken into account.
The modeling of the soil subsystem is based on considerations by
Wolf [
20]. An arbitrary load application area, in this case, the approximately rectangular contact area of the drum and soil (half side lengths are denoted by
a and
b), is converted into a circle with an equivalent radius, which represents the top surface of a downward unlimited truncated cone and is approximated as a Kelvin–Voigt element. The Poisson’s ratio
of the soil defines the opening angle of the cone. Vibratory roller compaction is primarily used for cohesionless, compressible soils with Poisson’s ratios of
, for which
Wolf [
20] does not consider any resonant mass. The associated dynamic parameters of the elastic isotropic half-space, the elastic spring stiffness
k and the dashpot coefficient
c, respectively, are defined as [
20]
where
G is the shear modulus of the soil, and
is the density of the soil. The equation for the calculation of the dashpot coefficient presented by
Wolf [
20] is adapted in Equation (
4) by introducing the factor
. Previous studies [
3] have shown that increased damping is required in the one-dimensional model to represent spatial effects and to obtain realistic results.
The essential simplification of the cone model according to
Wolf [
20] can be summarized in the consideration of only one-dimensional waves, which cause a vertical radiation into the soil with constant velocity. The predominantly vertical loading of the soil by a vibratory roller thus justifies the use of
Wolf’s model [
20] for the mechanical model presented in this paper. However, the Kelvin–Voigt element is only able to describe the elastic material behavior like it is observed in the far-field of dynamic compaction work. Therefore, an additional elastic spring with stiffness
is introduced, which is connected in series with the Kelvin–Voigt element. The spring deforms proportionally during loading (compression) but remains locked otherwise (during unloading and loss of contact) and, therefore, simulates the near-field of the load application area and enables the consideration of plastic deformations during the loading phase. The elastic half-space remains unaffected. The geometric damping (radiation damping) in the far-field is captured by the damping parameter
c. Material damping in the near field is not taken into account since it plays only a minor role [
21]. Two approaches for the estimation of the spring stiffness
are presented in [
1]:
Option 1: The stiffness of the spring
is estimated by determining the permanent plastic deformation of the soil due to a static roller pass in the field [
1]:
where
is the accumulative plastic deformation gained during the loading phase, which equals the plastic deformation of the soil after the roller pass. This means that the determined stiffness may well be greater than the stiffness simulating the elastic soil behavior. For well-compacted soils, this will generally be the case.
Option 2: A dimensionless condition factor
(
) is defined that describes the ratio of elastic to elasto-plastic deformations and provides a relationship between
and
k [
1]:
where
and
are the absolute elastic and plastic deformations of the soil, and
is a compensation variable to account for the periodic increase in vertical displacements of the soil under cyclic loading, each with respect to the coordinate system of the initial state of the mechanical model. The better the soil is compacted, the lower is the plastic part of the deformations, and the condition factor
tends towards a limit value close to 1.
A predefined condition factor thus dictates a certain spring stiffness as a function of the elastic spring stiffness from the cone model. The series connection of the single spring and the elastic Kelvin–Voigt element results in an additional degree of freedom in the soil model (position coordinate of the absolute displacement ), which is, however, coupled to the absolute displacement of the soil .
The soil model presented thus takes into account the fact that deviating conditions from more distant soil zones prevail in the load application area and thus also enable the comparability and separation of relative plastic () and elastic () deformations. For the calculations and results in this paper, a condition factor of is defined to simulate the“partial uplift” mode of operation since this is the desirable mode of operation for vibratory roller compaction.
5. Variable Drum Contact Width
The contact area between the two subsystems decisively determines the dynamic interaction between the drum and soil. Due to the cylindrically curved drum geometry, the contact area is a variable quantity during the loading phase. The drum increasingly penetrates the soil, and the contact width between the subsystems steadily increases, while the length of the load area remains largely constant and equals the length of the drum.
This circumstance causes significant nonlinearity in the interaction system and thus significantly influences the reaction forces in the soil, which in turn have a decisive influence on the motion behavior of the drum. In contrast to previous research [
1,
8], the presented model enables the consideration of the variable contact width
during the loading phase. The instantaneous drum contact width is determined using the circular arc segment shown in
Figure 3, which is defined by the sagitta
and the circular chord
s.
The sagitta
is calculated with respect to the relative elasto-plastic deformations of the soil. Therefore, the compensation variable
is subtracted from the absolute soil displacement
:
With the known drum radius
r and the calculated elasto-plastic deformations
, the circular chord
s and, subsequently, the circular arc length
can be determined according to Equations (
12) and (
13):
During the unloading phase, the width is assumed to be constant and corresponds to the value at the end of the loading phase since only elastic deformations occur.
If the drum changes directly from the unloading phase to a loading phase without losing contact with the soil beforehand, the contact width has a finite value at the beginning of the loading phase. However, if the drum loses contact and hits the soil again in the subsequent loading phase, an infinitely small contact width would result at the time of contact initialization. This circumstance would inevitably lead to singularities in the further course of the calculation and consequently to numerical problems in the calculation (see Equations (
3), (
4) and (
11) to (
13)). Therefore, a finitely small value for the width of the contact width at the beginning of the loading phase is defined by extending Equation (
11) by the variable
in the case of contact loss:
The variable drum contact width results in a variable spring stiffness
k and dashpot coefficient
c of the Kelvin–Voigt element during the loading phase. A definition of the stiffness of spring
via a condition factor
(see Equation (
6)) also makes
dependent on the variable drum contact width. Therefore, an incremental calculation is required.
It is noted that the settlement trough that forms during the loading phase is not taken into account in the presented method for the determination of the contact width and the calculated arc length
is, therefore, slightly overestimated (see
Figure 3). However, in non-cohesive (well-compacted) soils, the influence of the settlement trough plays only a minor role with regard to the drum contact width.
6. Calculations and Results
In this chapter, selected results from simulations with the presented mechanical roller–soil model are presented.
Table 1 gives a summary of the roller (HAMM H13i [
22]) and soil parameters used for the calculations. The soil parameters were selected to ensure the “partial uplift” mode of operation for the roller. To guarantee a steady state of the simulation, 90 revolutions of the eccentric mass were simulated (simulation duration = 3 s). The following graphs show a section of four excitation periods.
Figure 4 illustrates the development of the variable contact width of the drum during the loading phase (marked “L”, highlighted in green). Since the roller operates in the “partial uplift” mode of operation due to the simulation parameters listed in
Table 1, the quantity
according to Equation (
14) must be taken into account when calculating
for the loading phase so that the contact area between the drum and the soil assumes a finite value in the first loading increment and singularities are avoided.
In the subsequent unloading phase (marked “U”, highlighted in red), the drum contact width is kept constant and thus corresponds to the value at the end of the loading phase. The soil follows the upward movement of the drum until the two subsystems finally separate. During the time of contact loss (marked “T”, highlighted in blue), the contact width
assumes the value zero. In the modeling, however, it must be kept constant (as well as
) in order to determine the soil parameters (see Equations (
3) and (
4), and
Section 5, respectively) via the cone model and to be able to calculate the motion quantities of the soil in the creep-back phase.
Figure 5 shows the development of the soil parameters
k,
c and
in the individual operating phases (loading, unloading, and loss of contact) resulting from the course of the drum contact width shown in
Figure 4. Due to the specification of the condition factor
(see
Table 1), the spring stiffness
also changes during the loading phase (see Equation (
6)), which means that all stiffness and damping parameters increase during the loading phase.
For the subsequent operating phases unloading and loss of contact, a constant course for the parameters of the Kelvin–Voigt element results due to the given contact width (see
Figure 4). According to Equation (
6), the state factor
in the purely elastically modeled unloading phase approaches its limit value close to 1 and the plastic parameter
, thus tending to infinity. In the loss of contact phase, on the other hand, the soil provides no resistance. Therefore, the state factor
assumes its second limit value and becomes zero, as does the plastic parameter. However, since the spring stiffness
is not included in the calculations during unloading and loss of contact anyway, it is kept constant in the modeling, illustrated by the dotted line in magenta in
Figure 5.
The determined displacement quantities are shown in
Figure 6. According to the described modeling, the roller “moves towards the center of the earth” with constant speed due to the increase in the plastic deformations in each loading phase. In reality, due to its travel motion, the roller constantly moves “uphill” from the impressed deformation to the soil surface that has not yet been passed over and is still higher at this point in time. The work for this constant “uphill travel” must be conducted by the drive of the roller. Vibration displacement and total deformation are not identical due to the plasticity taken into account. The variable
compensates for this circumstance and allows the calculation of the relative displacements
from the calculated absolute value
. During the operational phase change from unloading to loss of contact, the drum and the soil separate from each other. In this phase, the soil shows elastic creep back towards its initial position until the drum changes to the next loading cycle and there is again a steady increase in plastic deformations. In the subsequent unloading phase, the soil reacts purely elastically because the single spring is locked (
) and thus cannot expand. The maximum of the plastic and elastic deformations does not occur at the same time due to the dynamics in the system (mass inertia and damping).
Figure 7 shows the temporal course of the determined soil contact force
(black) and its components for four excitation periods. The contact force is calculated from its three components, the inertia force
of the vibrating mass shifted upwards by the static axle load
(blue), and reduced by the force
(green) resulting from the rotation of the eccentric mass. If the blue and green lines are superimposed, the soil contact force does not contribute to the inertia force (see Equation (
2)). During this period, the drum and the soil subsystems move separately, and the soil contact force becomes zero. A change of the operating phase from loss of contact to loading, or from unloading to loss of contact results in a pronounced kink in the course of the inertia force since the drum acceleration changes abruptly at these points. The inertia force shows recurring movement during each unbalance revolution, amplitudes and contact times are identical, which, in turn, means that the roller operates in the “partial uplift” mode. The maximum contact force coincides with the change from the loading to the unloading phase.
The soil contact force may also be calculated by means of the soil subsystem. The sum of the elastic spring force and the elastic damping force provides the soil reaction force, which becomes zero in the loss of contact phase since the spring and damping forces have the same magnitude. At the moment when the drum reaches the lower reversal point of its motion, i.e., the maximum relative vibration displacement per unbalance revolution, the vibration velocity of the drum and, consequently, also the velocity-proportional damper force become zero.
Figure 8 shows the course of the plastic spring force
as a function of the absolute plastic deformations
. These plastic deformations are the cumulative sum of the permanent deformations of each loading phase of the soil (compare
Figure 6). The depicted relationship thus illustrates, to a certain extent, the geometric nonlinearity of the system resulting from the consideration of the variable drum contact width and the resulting variable spring stiffness
in the loading phase. If the calculations were performed with a constant contact area between the drum and the soil—and therefore a constant stiffness
—the course of
would be a linear saw-tooth.
The force–displacement diagram in
Figure 9 shows the relation between the soil contact force
and the corrected relative vibration displacement of the drum
and illustrates the loading and unloading of the soil by the drum. The corrected relative vibration displacement (“folded back” to the initial position of the coordinate
) is used on the abscissa instead of the total or absolute amount of vertical displacements of the drum in order to ensure the comparability of the diagrams from the simulation with measured force-displacement diagrams.
The force-displacement line of the elastic spring force (red) illustrates the main nonlinearities observed during compaction with vibratory rollers: the variable drum contact area—caused by the cylindrical geometry of the drum—and the plastic deformations occurring in the loading phase result in a curved force-displacement line of the elastic spring force during loading. The kink in the soil contact force at the operating phase change from loading to unloading () also indicates plastic deformations in the loading phase. Furthermore, the inherent nonlinearity of the system due to the changing contact conditions at the operating phase transitions from unloading to loss of contact, or from loss of contact to loading leads to a kink in the velocity-proportional damper force (magenta).
As can be seen from the determined displacement quantities in
Figure 6, the unloading phase represents a purely linear elastic process. Due to the chosen representation of the forces over the corrected vibration displacement of the drum
, the force-displacement line of the elastic spring force in the corresponding operating phase results in a “loop”. This loop disappears for a plot of the elastic spring force over the absolute vertical displacement of the drum. The unloading phase is terminated by the loss of contact (
); thus, the essential nonlinearity of the system, according to the present contact problem, is captured accordingly. The comparatively simple modeling of the complex interaction system is therefore able to adequately represent the essential nonlinearities and practice-relevant characteristics that occur during dynamic compaction with vibratory rollers.
7. Discussion
To further illustrate the relevance of a consideration of the variable contact geometry between the drum and soil, additional simulations were performed with a constant drum contact width. The contact width between the drum and soil was assumed to be
(see
Figure 10) and
(see
Figure 11), respectively, to ensure comparability with the results in
Figure 9. The essential nonlinearity resulting from the curved shape of the drum is thus not taken into account in
Figure 10 and
Figure 11. The soil parameters
k,
c, and
therefore assume constant values during the entire simulation.
Compared with
Figure 9, the velocity-proportional damper force in
Figure 10 shows a much steeper curve at the beginning of the loading phase, while the elastic spring force shows a flatter curve. A comparison of the force-displacement diagrams in
Figure 10 and
Figure 11, which only differ in the choice of the constant drum contact width, highlights the influence of the predefined contact width on the shape of the force-displacement diagram. The contact area chosen in
Figure 11, which is twice as large, already causes an almost abrupt increase in the damper force (and thus also in the soil contact force) at the beginning of the loading phase and has a significant influence on the fullness and shape of the force-displacement diagram. If the plastic parameter
were infinitely large (purely elastic behavior of the soil), or if the spring with stiffness
were also connected in parallel with the dashpot
c, a vertical jump in the damping force
(or soil contact force
) would also result at the transition from the loss of contact phase to the loading phase.
In reality, however, the area between the drum and the soil at the moment of contact initiation is very small due to the cylindrical curvature of the drum, and thus the contact stress is significant. Therefore, the soil beneath the drum fails in a manner similar to a static ground failure. The impulse due to the drum impact is thus weakened accordingly, and the contact force
does not increase abruptly but represents a rather continuous function starting from zero [
1].
Force-displacement diagrams with such a shape as shown in
Figure 10 and
Figure 11 can be observed with good approximation for LWD tests (constant, circular contact area) [
21] but do actually not occur for roller compaction in reality [
1,
17]. Consequently, a consideration of the variable contact width is not only justified but also necessary in order to be able to perform realistic simulations.