Note to testers

VPI software is currently in beta testing stage.

If you wish to test the software and helps us move forward, please leave a comment on this page or contact us at one of the following addresses:

Dave Cebon - dc29@hermes.cam.ac.uk
Riccardo Isola - riccardo.isola@nottingham.ac.uk

Thanks,

The VPI team.

Read More

permanent behaviour

Permanent deformation in UGMs is mainly related to resilient stresses or strains and number of load cycles. The effects of these two factors are often separated and the models for predicting permanent deformations can assume the form of:
                                                                                                               (1)
Amongst the numerous models that are available in the literature, Hornych (2004) shows that the Gidel model is particularly suitable to describe the variation of permanent deformations with stress level:
                                                                                              (2)
The link between permanent deformation and number of load cycles can instead be expressed by the simple Sweere model:
                                                                                                                                                (3)
or by the slightly more elaborated Hornych model:
                                                                                                                          (4)
All these models are derived empirically by fitting curves to laboratory data and do not take into account the effect that moisture content has on the permanent behaviour of the granular material. In theory, a different set of fitting parameters would need to be calculated for each value of moisture content considered. The ME-PDG, instead, uses a modified version of the Tseng and Lytton model that relates permanent strains to resilient strains, number of cycles and moisture content:
                                                                                                                      (5)
where β, ρ and ε₀/ε_r are functions of water content.

Read More

Resilient Behaviour

In dealing with granular materials, the elastic properties are often defined by the resilient modulus (which replaces the modulus of elasticity to indicate the nonlinearity of the behaviour) and Poisson’s ratio. A common approach in dealing with the stress dependency of granular materials stiffness is the K-θ model, which expresses resilient modulus as a function of the sum of the principal stresses, or bulk stress.
                                                                                                                                (1)
Thanks to its simplicity, this model has been widely adopted for analysis of stress dependence of material stiffness. However, it has the drawbacks of assuming a constant Poisson’s ratio and of ignoring the effect of other stress parameters apart from the bulk stress. For this reason, several more elaborated versions of this model have been suggested such as, for instance, the so called “Universal Model”, first suggested by Uzan and then adopted by the ME-PDG:
                                                                                             (2)
This model includes the 3D effect of shear, or deviator, stress by means of the octahedral stress. Most of the models of this family have empirical origins and are mainly used in the US.
A more theoretical approach is taken by those models that characterise the non-linear behaviour by decomposing stresses and strains into volumetric and shear components. In this case, resilient modulus and Poisson’s ratio are replaced by bulk and shear moduli:
                                                                                                                                                  (3)
                                                                                                                                                (4)
where
                                                                                                                                 (5)
                                                                                                                                            (6)
                                                                                                                                 (7)
                                                                                                                            (8)
In 1980 Boyce used this type of approach to develop a theoretical nonlinear elastic model for the stress-strain relationship of granular materials:
                                                                                                                         (9)
                                                                                                                         (10)
this model has then been employed by many other researchers, particularly in mainland Europe. It is important to note that it is possible to express resilient modulus and Poisson’s ratio in terms of bulk and shear moduli:
                                                                                                                                   (11)
                                                                                                                                          (12)
This can be useful considering that our current approach to including the effect of moisture content in granular materials makes use of the resilient modulus.

Read More

UGMs - Introduction

Over the years, many researchers have studied the complex behaviour of granular materials, using laboratory and in situ testing techniques. An extensive literature review was carried out by Lekarp et al. (2000) to collect findings from previous research and summarize in two companion papers the state of knowledge on resilient and plastic properties of granular materials. In this work the authors point out that these properties are affected by numerous factors such as stress, density, grading, moisture content, stress history, particle shape and load frequency, nonetheless the effect of stress parameters is certainly dominant.
If granular materials are simulated by means of Layered Elastic algorithms or using the Method of Equivalent Thickness, this stress dependency is usually dealt with through iterative processes (Figure 1):
Figure 1: Iterative determination of resilient modulus
The complexity of the problem meant that a large number of models can be found in the literature. An interesting review of available models for the prediction of permanent deformations in UGMs has also been conducted as part of the SAMARIS project by Hornych et al. in 2004.
Hereafter we present a selection of the main approaches to resilient and permanent behaviour of UGMs that might be implemented in the WLPPS.

Read More

Results 2

Normalised Aggregate Force and Fourth-Power Forces
Figure 1 is a comparative plot of the fleet normalised aggregate forces generated by the four simulation methods.
Figure 1: Comparison of fleet normalised aggregate force histories for each of the four simulated vehicle fleets
The ‘reference’ and ‘target’ time histories agree very well, while the time histories from the phase shifted and randomised QCM fleets appear to follow a similar pattern. This comes down to the inherent dynamics of each model; the ‘reference’ and ‘target’ histories were generated by pitch-plane models, the randomised QCMs and phase-shifted models used QCMs that lack pitch interaction and wheelbase filtering.
SRI Curves
Figure 2 and Figure 3 are comparative plots of the SRI curves generated by each of the four simulated fleets and a re-plotting of the region 0.8<SRI<1, respectively. From Figure 2 and Figure 3 it can be seen that the ‘target’ and the phase-shifted SRI curves agree well, as expected.
Due to the fact that vehicles with similar suspension characteristics correlate well with each other, all curves show a sharp peak close to the maximum value of one and a secondary peak in the region of SRI=0.85, corresponding respectively to those vehicles that are entirely air sprung and those that have leaf-sprung trailers.
Figure 2: Comparison plot of SRI curves for each of the four simulation methods
Figure 3: Comparison plot of SRI curves re-plotted for the range 0.8<SRI<1
Computation Time and Summary of Results
Table 1 provides a summary of computation time and accuracy metrics for each simulation method.

Method	Simulation Time (One Run) [sec]	Simulation Time (20 years, @ 1 per week) [days]	R2 of SRI Curves Relative to ‘Reference’	R2 of SRI Curves Relative to ‘Target’	Correlation of ATFs to ‘Reference’ Fleet	Correlation of ATFs to ‘Target’ Fleet
Reference	147600	1777	1	-	1	-
Target	26000	313	0.91	1	0.99	1
Random QCM	24000	290	0.76	0.76	0.75	0.74
Phase Shifted QCM	90	1.5	0.53	0.89	0.52	0.48

Table 1: Summary of simulation time, goodness-of-fit, and correlation of normalised aggregate forces for all four simulated fleets
Although the phase-shifted method includes the overhead time of generating the ‘target’ SRI curve first (i.e., from measured data or more realistic simulations), this is a ‘one-time’ cost and the phase-shifted models still represent a significant decrease in computation time. For example, if one was simulating 20 years worth of traffic in weekly intervals, 1040 separate traffic calculations would be required. Using the randomised pitch-plane models of the ‘target fleet’ would take approximately one month of CPU time. Conversely, the phase-shifted method would require only 1.5 days of CPU time, accounting for the overhead of generating the ‘target’ SRI distribution; a 99.5% reduction in computing time.
Given the substantial computational benefit of the phase-shifting method, and the excellent agreement of the SRI statistics, it is believed that the phase-shifted QCMs are still the best available method for simulating dynamic tyre forces for whole-life pavement performance calculations in which the effects of millions of axle loads need to be simulated over the lifetime of the road surface.

Read More

Spatial repeatability 2

Spatial repeatability arises because trucks are similar in weights, dimensions, and dynamic characteristics and travel at similar speeds. As a result, each vehicle will apply its peak forces at approximately the same places along the pavement surface (Cole, Cebon 1992).
Cole et al. defined the spatial repeatability index (SRI) as the correlation coefficient between a dynamic tyre force histories,
                                                                                                         (1)
where x and y are dynamic tyre forces histories,
m_x and m_y are the mean forces of x and y, respectively,
and σ_x and σ_y are the standard deviations of x and y, respectively (Cole, Cebon 1992; Cole et al. 1996).
An alternate measure of spatial repeatability, suggested by Cole, is the fleet normalised aggregate tyre force,
                                                                                                                              (2)
where
F_i,jk is the tyre force at the i-th road point due to the k-th axle of the j-th vehicle,
m is the number of axles on each vehicle,
N_V is the total number of vehicles,
and is the mean of the double sum in the numerator (Cole et al. 1996).
The fleet normalised aggregate force gives a spatial (time-domain) picture of the cumulative pattern of traffic loading.

Read More

Examples 2

Example 1 – Full depth crack
A full-depth crack allows the water to reach the top of the granular layer, saturating it. In order to model how the moisture spreads from this point throughout the layer, we keep the top left element of the structure saturated (Figure 1 to Figure 3).
Figure 1: Moisture infiltration through a full-depth crack – t=0
Figure 2: Moisture infiltration through a full-depth crack – t=4
Figure 3: Moisture infiltration through a full-depth crack – t=8
Example 2 – Pipe leak
Similarly to the case of a full depth crack, the effect of a leaking pipe can be simulated by keeping a node internal to the structure constantly saturated (Figure 4 to Figure 6).
Figure 4: Moisture infiltration from a leaking pipe – t=0
Figure 5: Moisture infiltration from a leaking pipe – t=4
Figure 6: Moisture infiltration from a leaking pipe – t=8
Example 3 – Micro-cracks
The presence of micro-cracks over a certain length of pavement could be modelled by keeping the degree of saturation of some nodes underneath the asphalt layer constant but lower than 1 (Figure 7 to Figure 9).
Figure 7: Moisture infiltration through micro-cracks– t=0
Figure 8: Moisture infiltration through micro-cracks– t=4
Figure 9: Moisture infiltration through micro-cracks– t=8

Read More