Alize Mu v130 GB
Alize Mu v130 GB
Alize Mu v130 GB
version 1.3
User manual
January 2011
© LCPC - 2011
Official distributor: itech 12-16 rue de Vincennes 93100 Montreuil - FRANCE
Tel.: +33 1 48 70 47 41 - Fax: +33 1 48 59 12 24 – E-mail: [email protected]
www.itech-soft.com
Contents
Contents
6. Performing computation 25
In this manual, the "Alize-LCPC Routes" software will be also designated by both the terms "Alize-LCPC" and
"Alize". It is assumed that the user has been adequately trained in running the Microsoft Windows operating
system and is capable of mastering the various associated interfaces and peripherals.
Alize-LCPC Routes has been set up to implement the rational mechanical design method for pavement
structures, as developed by the LCPC and SETRA French organizations; this method constitutes the regulatory
approach for designing pavements throughout the French national road network and moreover has been
adopted by many other road project development agencies. The basis of the rational design technique is
presented in Appendix A1 for what concerns the pavement mechanical design, and Appendix A2 for pavement
frost-thaw verification, which is also considered by the program.
Alize-LCPC may be distributed in the complete integral version including the three modules Alize-mechanical,
Alize-frost thaw and Alize-back calculation, or in one of the three following limited versions: Alize-mechanical
alone, or Alize-frost thaw alone, or Alize-mechanical + Alize-back calculation. This User’s manual concerns
both the integral version and the three partial versions of the software:
The users of the integral version are concerned by the whole manual (11 parts) and the two
appendixes A1 and A2.
The users of the limited version Alize-mechanical are concerned by the parts 1 to 8 and Appendix A1,
The users of the limited version Alize-frost thaw are concerned by the parts 1 to 3, 9 to 10 and
Appendix A2,
The users of the limited version Alize-mechanical + back calculation are concerned by the parts 1 to 8
and 11, and Appendix A1.
The French Design guide for pavement structures, LCPC-SETRA 1994 which present in greater detail
the principles of the rational pavement design
The Catalogue of new standard pavement structures, LCPC-SETRA 1998
The Technical Guide of Variant Specifications, SETRA 2003.
These two later documents describe the application conditions of rational design to the pavement
structures of the French national road network, accordingly to the French official specifications.
The computation algorithm of the Alize-frost thaw module is derived from the Gel1d software, a MS-Dos
application developed by LCPC, not yet distributed since 2002. For more details about the Gel1d software,
refer to «Gel1d, Modélisation de la congélation des structures de chaussées multicouches», LCPC 1999. In
addition to this computation algorithm, the Alize-frost thaw module is built around a new graphical user
interface (GUI) including help facilities, which are not comprised in the original Gel1d Dos program.
Alize-mechanical module :
Definition of the pavement structure: thickness, elasticity parameters of the different layers and
interface condition between layers;
Definition of the loading applied at the surface (reference load or other loading, named “special
loading”);
Determination of the stresses and/or strains allowable by the different materials, according to their
damage law parameters and the traffic condition;
Computation of the stresses and strains created in the pavement materials by the traffic loads;
Graphical display of the mechanical computation results;
Assistance and support for the practical choice of both hypotheses and numerical values of the various
parameters need by the mechanical computation, in accordance with the specifications of the
Technical Guide LCPC-SETRA and/or the New Structures Catalogue;
Management of a library (the mechanicam library) including both standard materials which mechanical
properties are defined by the LCPC-SETRA documents referred above, and customized materials
defined by the user.
The protection of Alize-LCPC against hacking is ensured by means of a single-workstation key or network key,
which comes delivered with the software. This protection key must imperatively be installed onto either the
parallel port of the microcomputer (individual installation) or the server workstation (network installation),
preliminary to any use of the software.
In the case of a network key installation, it is also necessary to perform the network server addressing
procedure and install the electronic license manager. For this step, reference is directed to the key installation
notice (contained within the "install_key.pdf" file).
2.3 - Comments
It is recommended to accept all options proposed by the installation program, whether following the automatic
or manual installation steps.
In the Version 1.3 release, the configuration is limited to the following items:
Choice of the language used for the whole input and output (screen displays, files and printer) and the
help facilities: French or English (figure 3.3).
Definition of the reference load for the mechanical computation (cf. Part 5 below).
Choice of personal options for the Alize-mechanical module (figure 3.4 and 3.5):
– The default value of the Poisson coefficient,
– The default interface condition,
– The printer layout in case of computation including variants (cf $4.3), and the printer layout
for allowable values printing,
– The default value of the equivalent temperature for the bituminous materials.
Choice of personal options for the Alize-frost thaw module: Display or not display of the results of the
simplified frost-thaw method Iatm = f(Qpf) , cf. figure 3.6 and $10.2.
Personal option “Wide Windows” (figure 3.5) : if selected, the width and height of some windows are
automatically increased in order to solve readability problems, mainly occurring with the Operating
system Windows Vista.
Figure 3.3: Menu “Customize” of the main menu bar, choice of the language
Figure 3.4: Menu “Customize” of the main menu bar, personal options - default interface condition
and Poisson value, Printer layout for the printing of mechanical results and allowable values
Figure 3.5: Menu “Customize” of the main menu bar, personal options - default equivalent
temperature value for bituminous materials and Wide-windows option
Figure 3.6: Menu “Customize” of the main menu bar, personal options – Display or not display
the results of the simplified frost-thaw method Iatm = f(Qpf) (Alize-frost thaw module)
The results of the computation can be saved – if requested –in output Ascii files as follows :
Allowable values computation files : extension .adm, cf. Part 8 below of the present manual;
Results of mechanical computation : extension .res, cf. Part 7 below;
Results of frost-thaw computation : extension .res, cf. §10.2 below;
Results of backcalculation : extension -retro.res, cf. Part11 beloow.
The various layers of material constituting the structure have constant thickness; moreover, their expansion
within the horizontal plane XoY is infinite. Expansion along the vertical direction ZZ of the lower layer of the
multilayer foundation, which generally represents the substratum or supporting soil, is also assumed infinite.
The following set of parameters are needed for each layer :
thickness H;
Young's modulus E of the material;
Poisson's ratio of the material (denoted "Nu" in Alize-LCPC); and
interface conditions at the layer's upper and lower extremes, with the adjoining layers.
Three types of interface are available in order to characterize how the interface functions between adjacent
layers: bonded, sliding, or semi-bonded. The semi-bonded interface condition is specified by the Technical
guide French pavement design LCPC-SETRA, for modeling the contact between some materials. In case of semi-
bonded interface, two successive computations are automatically performed, the first one considering bonded
interface and the second one considering sliding interface. The semi-bonded condition is then considered as
the mean value of the stresses and strains resulting from these two computations.
A base structure made of 3 layers is then initialized within the main window. Depending on the pavement
structure to design, the user must then modify this basic model structure and complete all the data boxes so
that define the structure mechanical characteristics, as reproduced by the example in Figure 4.2.
For this purpose, the following set of commands, functions and help features have been made available:
Table 4.1: List of the commands available for the creation of a new pavement
structure for mechanical computation
The decimal separator used for entering numerical values must be that for which the Windows
operating system has been configured on the host PC (see Settings/Regional and Language
Options/Decimal Symbol).
It is strongly advised to fill the structure title boxe in order to facilitate its subsequent identification,
e.g. during future use of the data or printout of computation results.
The maximum number of layers of materials constituting the pavement structure is limited to 15, with
a minimum of 1.
The default interface condition may be defined using the “Customize” command of the main menu
bar, which open the window “Alize-Lcpc – Personal options” (see Part 3 above and figure 4.4). The
default value of the Poisson coefficient, and the default value of equivalent temperature of
bituminous materials, may also be modified using the same command.
The computation of internal stresses and strains within the materials is performed at two levels in
each material layer. These levels implicitly consist of the interfaces between layers, i.e. those levels
at which maximum internal stresses are generally obtained. The "Computation levels" command
enables modifying this implicit choice.
For defining one layer of the pavement structure as a predefined library material, click the yellow box
“Other”, which automatically opens the window “Material library – value of E and Nu” (figure 4.3).
Then it is possible to select either a standardized material as specified by the Pavement design guide
LCPC-SETRA, or a personal material of the user’s library. For more information about the material
library, see $8.2.
Figure 4.3: Use of the Material library for defining the mechanical parameter
of a given layer, by clicking the boxe “Other” of the MAerial type column
The help button "Nature of interfaces" recalls the set of hypotheses stipulated in both the French
pavement design guide LCPC-SETRA and the 1998 Catalogue (figure 4.4).
The help button "Minimum technical thicknesses" recalls both the minimum and maximum thickness
values for material implementation, which had been incorporated when establishing the 1998
Catalogue (figure 4.5).
Figure 4.4: The help screen “Interfaces”, activated by the command “Nature of interfaces”
Figure 4.5: The help screen “Minimal and maximal thicknesses of materials according
to 1998 Catalogue”, activated by the command “Mini-maxi thicknesses”
For a given structure, it is possible to define up to 25 computation variants. The command “See-
manage variants” opens the window “Alize-LCPC – Manage-definition of alternative data” (figure 4.6).
With respect to the base structure data, these computation variants pertain to either:
– variations in thickness of one or several layers, or
– variations in Young's modulus of one or several layers.
Defining computation variants simplifies the search for a solution by sequencing automatically the
computations specific to each predefined variant.
Consulting the forms of 1998 Catalogue of new pavements may also be done using the command ”1998
Structures catalogue” of the window ”Computation of allowable values” (cf. Part 8).
The "File/Open" and "File/Save" commands are also executed according to a standard Windows procedure; they
are used to open and save the files containing previously-defined structures.
Figure 4.7 : Commande for consulting and/or importing structure from the 1998 French Catalogue
The reference load is defined in the Alize-LCPC program by use of the command "Configure Alize/Reference
load" of the main menu bar. The reference load defined in this manner will be recorded by the software for
future use, until a different reference load is eventually set.
The definition of a reference load is not mandatory however. Should one not be predefined, the computations
would thus have to be conducted, for a loading applied at the pavement surface, with a loading named "special
load", which is to be defined as indicated in Section 5.2.
The vertical computation profiles, associated with the various computation levels (see Section 4.2), are used to
establish those points at which the mechanical calculation of internal stresses will be performed by the Alize
computation . The predefined vertical computation profiles for the reference load are as follows (conventional
profiles as guided by the rational design approach):
Case of the reference load composed of an isolated wheel: a vertical profile of a single computation,
corresponding with the axis of the circle that defines the load;
Case of the reference load composed of a dual-wheel: two vertical computation profiles, one
corresponding to the axis of the circle defining one of the two loads, and the other corresponding to
the central axis of symmetry of the dual-wheel.
The computation profiles connected with the reference load cannot not be modified. If such a modification is
requested, use the Special load definition procedure as described in Section 5.2.
The input of these two data sets (loading characteristics and vertical computation profiles) is activated by the
"File/New/Special loads" command of the main menu bar.
Saving to a file
Once a special loading has been completely defined and before continuing with use of this application, it is
recommended to save the entire set of data relative to the special loading, by means of the "File/Save as"
command on the main menu bar (see Section 4.1). The general comments expressed in Section 4.1, with
respect to both saving and opening data files via the standard Windows dialog boxes, are also applicable
herein.
Figure 5.3: Definition of a special load: Example of a 16-wheel convoy, definition of loads and
vertical computation profiles
This method serves to reproduce those loadings for which assimilation of the tire-pavement contact surface by
a circle or a small number of circles would prove to be overly simplistic.
5.4 - Navigating between the Structure window and Special load window
It is possible that both the "Pavement structure definition" and the "Special loading definition"
worksheets are simultaneously active. Only one of these two windows may be displayed on the
screen, while the other must remain invisible. The Window command on the main menu bar enables
switching back and forth from one window display to the other.
6 6. Performing computation
General case:
The Alize computation is launched through the "Compute/Alize" command on the main menu bar (figure 6.1).
The "Execute mechanical computations" window (see figure 6.2) allows defining both the pavement structure
and the loading for which the computation has been requested. All of the following combinations are indeed
possible:
computation for the pavement structure that could either be displayed in the "Structure" window or
retrieved from the "Structure" file;
reference loading or special load;
in the case of loading by means of a special load, the computation would eventually be displayed in
the "Loading" window or retrieved from the "Loading" file;
a "Standard" or "Grid-seca" type of computation.
Figure 6.2: Initiating the Alize computation: Choice of data and type of computation
Figure 6.3: Execution of the Alize computation using the "Fast computation (ref. load)" command
The Grid-seca computation actually relies upon the supplemental data that have been defined using the
"Definition of a Grid-seca computation" window (cf. figures 6.4, 6.5 and 6.6).
The supplemental data necessary for running the computations are as follows (see figures 6.4 and 6.5):
grid parameterization: in a grid-seca computation situation, the points on the XoY plane that outline
the computation profiles create a cluster of rectangular points, as characterized by grid length, width
and the dimensions of its elementary mesh, which serve to designate the grid "interval". The Alize
program automatically proposes a computation grid on the basis of the loading geometry (either
reference load or special load). The user can then modify, if so desired, this initial grid by means of
altering the following parameters:
– Grid symmetries: It would be possible to deactivate the eventual symmetries inherent in the
grid;
– Border: This element represents the distance between the grid boundaries and the center of
the closest load. The overhang is identical everywhere along the XX and YY directions;
– Interval: the grid interval directly influences the number of computation profiles, the amount
of memory required and the Alize computation time. For purposes of illustration, the
maximum grid size is approximately 850,000 computation profiles on a P3 type of PC, with
128 megabytes of RAM.
computation criteria: In order to limit both the computation time and volume of results, it is advised
to reduce the number of parameters to be computed, by specifying these paramters. They may vary
depending on both the material layer and computation point within each layer. In order to simplify
declaration of these parameters, 6 options will be proposed, each containing a pre-established list of
parameter to be computed (table 6.1). Option selection is then carried out by one or several
successive clicks on the box corresponding to the chosen computation level.
Displayed Number of
Designation of the parameters to be computed
indication parameters
SigmaT 7 SigmaXX, SigmaYY,Sigma1, Sigma2, Teta (XoY), p and q
EpsilonT 7 EpsiXX, EpsiYY, Epsi1, Epsi2, Teta (XoY), EpsiV and EpsiD
SigmaZ 3 Sigma ZZ, p and q
EpsilonZ 3 EpsiZZ, EpsiV and EpsiD
SigXX...ZX 10 SigmaXX, SigmaYY, SigmaZZ, SigmaXY, SigmaYZ, SigmaZX, Sig1, Sig2, Teta
(YoZ) and W
EpsXX...ZX 10 EpsilonXX, EpsYY, EpsZZ, EpsilonXY, EpsYZ, EpsZX, Eps1, Eps2, Teta (YoZ) and
W
Note 1:
In the above list, Epsi or Eps designates Epsilon (strain), Sig designates Sigma (stress), Teta the rotation angle
of stresses or strains, and W the vertical displacement.
Note 2:
In a situation with "Structures" data that contain several variants, the “Grid-seca type of computation is only
carried out for the first variant (i.e. base data). The computation of another variant requires to return to the
main "Structure" window and then transposing the data from the targeted variant, within the base data.
Figure 6.5: Execution of the Grid-seca computation: View of the grid, example
7 7. Results of mechanical
computation
"Results" window:
This screen primarily features the following:
a review of the structure serving as the object of the computation; the loading identifier used for the
computation is also recalled herein;
a total of 8 tables, displayed one by one, presenting the following set of results:
Table 1:
Stresses and strains at the computation points. At each computation level, the minimum strain
value (EpsT) and minor primary stress value (SigmaT) within the horizontal plane XoY, along with the
maximum values of both strain (EpsZ) and stress (SigmaZ) in the ZZ direction.
In more common design applications, this table serves to summarize the set of results used for directly
designing the structure.
Table 2:
Stresses and strains at the computation points. Localization and orientation of the minimum EpsT
and SigmaT values in the XoY plane, as well as the maximum EpsZ and SigmaZ values in the ZZ direction,
as listed in Table 1. The notations used here are the following:
– - in the case of computations with the reference load: R = wheel axis, J = twinning axis, m =
non-oriented XX or YY direction (stress rotation), see Table 8;
– - in the case of a special loading: Pk = vertical computation profile n°k, m = non-oriented XX
or YY direction (stress rotation), see Table 8.
Table 3:
Strains in the XX, YY and ZZ directions, at each computation level and at the selected vertical computation
profile (see the choice of profile displayed using the vertical profile cursor).
Table 4:
Stresses in the XX, YY and ZZ directions, at each computation level and at the selected vertical
computation profile.
Table 5:
Shear strains XY, YZ and ZX at each computation level and at the selected vertical computation profile.
Table 6:
Shear stresses XY, YZ and ZX at each computation level and at the selected vertical computation profile.
Table 7:
Primary major and minor strains within the XoY plane, and angle of rotation with the XX axis. Values
calculated at each computation level and at the selected vertical computation profile.
Table 8:
Primary major and minor stresses within the XoY plane, and angle of rotation with the XX axis. Values
calculated at each computation level and at the selected vertical computation profile.
– The deflection values calculated at the pavement surface computation points are also
presented. The vertical profile cursor enables selecting the profile to be displayed.
– In the case of inputting a reference load, the value of the radius of curvature is also
calculated within either the wheel axis or the twinning axis, depending on the nature of the
reference load.
Figure 7.1: Main window for the display of the mechanical computation results
The command “Drawing” allows to display the evolution curves of the main results with the no of the variant or
thickness (in case of thickness variants) or modules (in case of module variants), as presented in figures 7.2 and
7.3. Printing the evolution curves is possible using the command “Print”, or saving them in Ascii files (extension
***-var.res) using the command “Save”.
Note:
In a computation situation containing variants, the printout of results is only terminated once the "Computation
results" screen has been closed, so as to enable, if and when applicable, transcribing the results of several
variants onto a single print page.
The results of a Grid-seca computation may also be presented in terms of damage created by the load in the
pavement materials. In this case the parameters needed for the damage calculation are defined by mean of
the window “Fatigue and wandering parameters for damage integration”, which is activated by mean of the
buton “Damage parameters” of the Results window (figure 7.5). The procedure for implementation of damage
computation is detailed below, cf. $7.3.
Representation in the p, q or EpsilonD, EpsilonV diagram of the paths of stresses or strains in the various
materials submitted to the loading effect is also possible herein. It would thus be necessary to have indicated,
for these specific materials and computation levels, one of the SigmaT, EpsilonT, SigmaZ or EpsilonZ criteria
(see Section 6.3).
The "Min-max values" command establishes a table that contains, for each calculated parameter (see
Section 6.3), the minimum and maximum values obtained (figure 7.6).
An integrated 2D-3D graphics routine enables visualizing computation results in the form of color isovalue maps
(figure 7.8). The reference plane for establishing such maps must be the horizontal XoY plane. In order to
automatically switch from the 2D representation to the 3D representation, the user will simply have to press
the appropriate option buttons. A right mouse click opens a dialog box that enables defining the personalized
visualization options (only available in English language).
The capture of screen-display of result curves and surfaces, for insertion as an image within Windows-based
word and image processing applications, is generated by simultaneously pressing the "Alt + Prt Sc" keys
followed by the "Edit/Paste special" command, which is generally proposed in the main application menu.
Figure 7.5: Results from a grid-seca computation: Choice of results visualization mode
Figure 7.6: Results from a Grid-seca computation: Example of results visualization in the "2D profiles" mode
Horizontal profiles of results along the longitudinal XX direction
Figure 7.7: Results from a Grid-seca computation: Example of results visualization in the "2D profiles" mode
Vertical ZZ profiles of results
Figure 7.9: Results of a Grid-seca computation: Dialog box for the saving of results in Ascii output files
Loading:
– The moving direction of the load at the surface of the pavement is necessarily the
longitudinal XX axis.
– The longitudinal XX axis is necessarily a global symmetry axis of the load (figure 7.10).
– A damage criteria related to flexural-tensile fatigue (treated materials) has been defined, as
detailed in Appendix A1. Then the parameter governing the damage is the minor flexural-
tensile strain Epsi2 (bituminous materials), or the minor flexural-tensile stress Sigma2
(hydraulic bounded materials and concrete).
Moreover, the damage calculation should be performed for any other component of the stress or
strain tensor. However, it should be pointed out that only the damages associated to the minor
tensile stresses and strains and the vertical compressive strains are considered by the the LCPC-SETRA
pavement design model.
The damage is calculated according to the continuous integration of the Miner law associated to the
selected parameter Epsi2, Sig2 or EpsiZZ. Integration is done along each longitudinal line of the grid.
Figure 7.11 details the Miner integration implemented in the Alize software. The damage due to the
passage of the load will be displayed as the transversal profile D = f(Y) resulting from the damage
calculated along each longitudinal line.
The transversal wandering of the load is taken into account, if requested, by mean of a standard
normal Gaussian cumulative distribution representative of the lateral position of the load trajectories
from the central axis. The damage resulting from this lateral wandering is computed by the Miner
law.
Input data for damage calculation : in addition to data relative to the pavement structure and the
loading, data for damage computation have to be defined in the window «Damage parameters» which
is activated by the command “Damage parameters” of the window “Results” (figure 7.12) :
– Fatigue parameters: allowable value for the considered layer and fatigue slope -1/b. The
determination of the allowable value may be done by mean of the « Allowable values »
window (cf. Part 8), using the aggressiveness coefficient value CAM = 1.
– Lateral wandering of the load: it is defined by the over-width of the wheel path (possibly
zero in case of completely canalized traffic), with corresponds to 2 standard deviations of the
normal Gaussian distribution.
The results of the damage computation are display by mean of the “2D Profiles” command of the “Results”
window, presented before in $ 7.2. The drawing of the transversal damage profiles D = f(Y) is activated using
the tick box “Damage D=f(Y)” as illustrated by figure 7.13. Once the parameter(s) concerned by the damage
has (have) been selected, the display of the corresponding damage profile is launched by mean of the button
“Drawing”. For each selected parameter, 2 transversal damage profiles are drawn, eg. without and with lateral
wandering.
Figure 7.10 : Load data for Grid-seca computation with damage calculation – Example.
The longitudinal XX axis is necessarily the load moving direction, and a symmetry axis of the load
Unloading
Déchargement
inter-pic
inter-peaks ul12 ul23
Tensile
Pics de
peaks
traction t1
t2 t3
Example : Airfield pavement structure loaded by a 6-wheel airfield landing gear – minor
strain profiles Epsi2 at the bottom of the base layer (GB3) along the longitudinal PX1 axis
b
Damage law (treated or untreated materials and soils) : s = K x N
with s = resilient stress or strain (= major Sigma2, Epsi2 or EpsiZZ, or other parameter …)
1 s
Elementary damage (Miner) : De = =[ ] with = -1/b
N K
-1
D = s ds
K
Integral expression of the Miner law (continuous integration along the moving
axis of the load) :
+
-1 ds(x) x si x 0
D=
K -
s pos[
dx
] dx with pos(x) =
0 si x < 0
Example : damage created by the passing of a 6-wheel airfield landing gear (triple-axle) :
Miner integration of the minor strain profile Epsi2 = f(x) shown above
1pas
1 passage : Dtridem = K1/b [ t1-1/b + t2-1/b + t3-1/b – ul12-1/b – ul23-1/b]
Npas 1pas
N passages (without wandering) : Dtridem = N Dtridem
Npas 1/b
Dtridem = tadm [ t1-1/b + t2-1/b + t3-1/b – ul12-1/b – ul23-1/b]
1/b
with : tadm = allowable strain associated to N passages according to the fatigue law
Figure 7.11: Expression of the damage according to the Miner continuous integration along the load moving axis
Figure 7.12 : Grid-seca computation, complementary data needed for damage calculation
8 8. Computation of allowable
values
The approach undertaken, along with the abbreviations and notations used and the expression of the various
allowable value laws, have all been drawn from this Guide, which should be consulted when seeking greater
precision.
The "Computation of allowable values" window is activated by means of the "Allowable values" command on the
main menu bar of the "Defining a structure" or "Defining a special load" window.
The following computational possibilities and features are offered as part of the "Allowable values" worksheet:
Computation of cumulative truck traffic (NPL), based on the typical set of data: average annual daily
traffic, rate of annual increase (tg - geometric) or (at - arithmetic), and service period (p);
Inverse computation of two of the five parameters mentioned above, based on data input from three
of these parameters;
Direct computation of allowable values from data on traffic, material choices, and the set of
associated parameters and coefficients that characterize material fatigue behavior (for treated
materials) or rutting (for untreated materials and soils);
Indirect computation of the allowable traffic for a given material, based on input data of the stress
borne by the material and the computation risk parameter, by means of inverting the fatigue law;
Indirect computation of the computation risk parameter, from input data on the stress borne by the
material and the traffic applied (figure 8.2);
Figure 8.2: Computation of allowable values: Example of inverse risk = f(Sigmat, traffic)
Consultation of the materials library that recalls the behavioral parameters of both standard materials
("system" materials) and customized materials defined by the user ("user" materials);
Consultation of the help functions extracted from both the Guide French pavement design and the
1998 Catalogue as a review of the following parameters:
– traffic classification,
– coefficients of traffic stress,
– computation risk for establishing allowable values;
Lastly, consultation of the entire set of structures established by the 1998 Catalogues of New
Structures.
The mechanical library includes 4 categories: bituminous material, materials treated with hydraulic binder,
concrete and untreated materials and soil. Each material category includes:
The standardized materials as defined by the technical uide Pavement design LCPC-SETRA 1994 and
the 1998 Catalogue. The values of the mechanical parameters of such standardized materials cannot
be modified (status = “system”).
The personal materials which have been introduced by the user in the library (staus = user). For
introducing, modifyning or removing a personal material in the library, use the command “Libraries”
of the main menu bar and follow the procedure launched by the butons “Add one material” or
“Remove material”.
Figure 8.3 shows how to access to the mechanical library from the Alize main window. Figure 8.4 shows the
mechanical library form for hydraulic bounded materials as example, including both standardized and personal
materials.
Figure 8.4 : Consulting the Mechanical library, example including personal materials
In common practice, the total thickness of the model is always more than twenty or thirty meters, so that the
temperature at the bottom of the model may be considered as constant; not depending on the time. For
instance, the Lcpc-Setra rational design method specifies that the bottom of the model is located 40 meters
above the top of the capping layer (also called the foundation layer in this manual).
As for the mechanical computation, the thickness of each layer of the Frost-thaw model is assumed to be
constant, and its lateral extension in the horizontal plane XoY is infinite (mono-directional model assumptions).
Each layer is defined by the following parameters:
Thickness H ;
Voluminal mass ;
Water content of the material ;
Thermal conductibility of the material in unfrosted situation ;
Thermal conductibility of the material in frost situation.
The Alize-frost computation kernel calculates the evolution with the time, of the temperature at the surface of
the model and at each interface level between adjacent layers. Theses results are also presented in the terms
of the evolution curves in the time of the frost quantity at the same levels, and frost index at the surface of
the pavement. The surface frost index is converted into atmospheric index. In the majority of cases, the
finalized result of frost computation is presented in term of the frost index allowable by a given pavement
structure (main unknown of the problem to solve), depending on the frost quantity transmitted at the top of
the capping layer. This frost quantity transmitted at the top of the capping layer is an input of the problem. It
can be evaluated by mean of a specific help-sheet included in the Alize-frost module, accordingly to the Lcpc-
Setra ‘s design Guide.
The evolution with the time of the frost front depth (s) is (are) also calculated and edited.
Alize-frost is based on a finite differences algorithm, which unknown parameters are the temperature and the
frost front depth speed. Three different thermal behaviour areas are considered:
The unfrosted area, which the thermal equilibrium follows from the heat equation without source
term ;
The frost area, which is characterized by thermal behaviour parameters different from the unfrosted
area ;
The frost front(s), which is (are) separating surface(s) between the unfrosted and the frost areas.
Latent heat absorption/dissipation phenomena develop on the frost front(s), which therefore present
thermal gradient discontinuity.
Alize-frost is usually used for the frost-thaw verification of pavements, which is based on the solving of the uni-
directional propagation heat equation in a multilayered structure. More generally, Alize-frost can also be used
for the computation of the temperature evolution curves with the time in a multilayered structure with or
without development of negative temperature. The loading of the model is defined both by the initial vertical
and the surface and bottom temperature profiles. For instance, the program can be used for the evaluation of
hot temperature in bituminous layers in connection with the evaluation of rutting in theses materials.
The building of a new pavement structure for the Alize-frost module is initialized by the command
Files/New/Alize-Gel of the main menu bar (figure 9.1).
Figure 9.1 : Command for the building of a new pavement structure for Alize-frost
Then two modes for building a new Alize-frost pavement structure are possible:
2-Case 2 : a pavement structure for mechanical Alize-Lcpc computation has previously been defined :
The direct exportation of pre-defined mechanical structure towards the Alize-frost module is proposed (cf.
figure 9.3).
The exported structure includes the same number of layers as the initial structure used for mechanical
computation. The thicknesses of layers are kept, with the exception of the lowest layer (see §9.2). In case the
initial mechanical structure includes one or several standard materials, then an equivalency is automatically
set between the standard materials of the Alize-mechanical library, and the standard materials of the Alize-
frost library. This equivalency is defined by the Pavement design Guide Lcpc-Setra (1994), which set for
standard materials both the mechanical and thermal parameters values.
Depending on the studied project, the user has to modify the pavement structure shown on the screen Frost-
thaw module: definition of the structure. He also has to fill and/or modify the different numerical value
defining the thermal behaviour of each material. To that effect, the following command, functions and help
routines may be used (cf. table 9.1) :
Table 9.1 : Commands for building and modifying a structure for Alize-frost computation
The decimal separator used for numerical value input is obligatorily the separator defined by the
Windows system as specified by the Windows configuration panel (cf. Configuration panel/Regional
and linguistic options/Decimal symbol).
It is recommended to give a title to the data (cf. Title line), in order to later facilitate its
identification, for instance as the data are again read on file in further Alize sessions, or as the
computation results are edited.
The lowest layer thickness is not directly given by the user. It is automatically calculated by the
program, from:
– The upper level Zpf of the pavement foundation in case of standard initial and boundary
conditions (cf. Lcpc-Setra ‘s rational design method, frost-thaw verification §9.4) ;
– The total thickness of the model defined by the vertical initial temperature profile To=f(z), in
case of special initial and/or boundary condition different from the Lcpc-Setra ones (cf.
§9.4).
The number of layers making up the whole pavement structure (including soil) is limited to 15. The
minimal number is 2.
Each layer may be considered as composed of a standard material defined by the Alize-Gel library. For
allocating a standard material to a layer, click the corresponding cell of the column “material type”
associated to this layer. Then the list of the standard materials composing the Alize-gel library is
shown on the screen (cf. figure 9.5). The name of the standard material selected for a given layer is
shown in the cell “material type”. In case on not standard library material, the material is named
“other”.
Alize-frost library material includes two types of material:
– Standard materials as defined by the Lcpc-Setra’s Design guide (1994). These materials are
marked by their alphanumeric acronym as defined by the Lcpc-Setra’s Design guide. Their
status is “system”.
– User personal material, which has been previously introduced by the user in the Alize-frost
material library. These materials are marked by the alphanumeric acronym defined by the
user. Their status is “user”.
For complementary instructions about the use of the Alize-frost material library, consult §9.6.
Figure 9.5: Alize-Gel material library main screen, activated by the command «Material type»
of the main form Alize-Frost-thaw module, definition of the structure
The initial values of the parameters Simulated duration et Output deltaT (see the main form Alize-
Frost-thaw module: definition of the structure) are fixed to 60 days and 12 hours. The simulated
duration is the time period which will be taken into account for the computation of the temperature,
frost indexes and frost quantities evolution curves (ie from t=0 to t=simulated duration). The data
Deltat is the time interval used for the definition of these evolution curves. The values 60 days and 12
hours are appropriate for performing almost all the standard Frost-thaw verifications according to the
Lcpc-Setra rational method. It possible to modify these values, depending on the problem to solve.
The commands Files/Open and Files/Save are used for reading and saving data in existing data files. They are
also operated by mean of the standard Windows dialog boxes.
The choice Lcpc-Sétra method lead to select the initial and boundary conditions of the Lcpc-Setra Frost-thaw
verification method, as defined by the Lcpc-Setra Design guide (1994). If this first option is ticked, then the
command Display assumption of the frame Complementary activates the visualization on screen of these
conditions (cf. figure 9.6).
Figure 9.6: Visualization of the initial vertical temperature profile and the upper and lower boundary condition
temperature evolution curves, as defined by the Lcpc-Setra Design Guide (1994)
Special assumptions:
If the option Special assumptions is ticked, then the command Define-change-display of the frame
Complementary data open the sheet Alize-Lcpc Frost-thaw, special assumptions.
This sheet is used for defining initial and boundary conditions different from the conditions specified by the
Lcpc-Setra Frost-thaw verification method. An example of such special temperature initial and boundary
conditions is shown on the figure 9.7.
The Special condition data files for the storage of not standard initial temperature profile and
boundary conditions have the extension .hag. These files have been created:
By saving the data defined by the user in the 3 tables To=f(z), Ts=f(t) and Tb=f(t) of the sheet Alize-
Lcpc Frost-thaw, special assumptions. This saving is operated by the command Save on of this sheet.
Or by mean of any Ascii text editor, apart from the Alize-Lcpc program. The command Model of file of
sheet Alize-Lcpc Frost-thaw, special assumptions indicated the format and the arrangement of the
Special condition data files wearing the extension .hag.
The command Drawing of the sheet Alize-Lcpc Frost-thaw, special assumptions activates the drawing
of the initial and boundary temperature conditions defined by the tables To=f(t), Ts=f(t) and Tb=f(t).
An example is presented in the figure 9.8.
The sheet Alize-Lcpc Frost-thaw, special assumptions makes it also possible to define the values of the
coefficients A and B of the linear relationship between the surface frost index Isurf and the
atmospheric frost index Iatm, namely Isurf = A x (Iatm-B). The Lcpc-Setra Frost-thaw standard
verification method specified the values A = 0.7 et B = 10°Cxday (cf. Annex A2, §A2.5).
9.5 - Calculation of the frost quantity Qpf allowable by the pavement foundation
Nota :
The value of the frost quantity Qpf allowable by pavement foundation is not in itself an input data necessary to
the computation by mean of the Alize-Lcpc Frost-thaw module. Only the input data for the definition of the
geometry and the thermal characteristics of the pavement structure, and the initial and boundary temperature
conditions are strictly necessary.
Yet the data of the frost quantity Qpf is recommended if possible, when usual (standard) frost-thaw
verification is performed.
Then the Qpf data will be used by the program to automatically calculate the surface frost index IS and
especially the atmospheric frost index Iatm allowable by the pavement structure, associated to the given value
of Qpf.
In almost the frost-thaw pavement verifications, the Iatm value is compared to the reference atmospheric frost
index of the project (depending on its geographic localization and the owner policy), leading to express a
decision about the admissibility of the modelized pavement structure for what concerns frost-thaw behaviour.
Figure 9.9 : Determination of the frost quantity Qpf allowable by the pavement foundation,
according to the frost-thaw verification method specified by the Lcpc-Setra Design guide (1994)
The command Qpf calculation of the main sheet Alize-Frost-thaw module: definition of the structure open the
form Qpf computation. This form enables the determination of Qpf according to the Qpf determination way
specified by the Lcpc-Setra Design guide. The Help menus Frost sensibility and Qm computation: Help
reproduce the recommendations of the Lcpc-Setra Design guide. One will refer to this Design guide for
complementary information. An example of Qpf computation is shown in the figure 9.9.
It also possible to directly define the value of the frost quantity Qpf without using the form Qpf computation.
The Qpf value may be directly given in the cell Frost quantity allowable by the pavement foundation Qpf of the
main sheet Alize-Frost-thaw module: definition of the structure.
The material library of the Alize-Frost thaw module is used in order to allocate to one or some layers of the
pavement structure, the thermal parameters of material included in the library.
The command Libraries/Alize-Gel library of the main menu may be launched from the main menu bar of both
the main sheets Alize-Lcpc mechanical computation and Alize-Lcpc-Frost thaw computation. This command
opens the Frost-thaw material library, and makes it possible to consult the library and to add/remove user
materials (cf. figure 4.3).
The form Alize-Lcpc Frost-thaw module, material library shows the list of the standard and user materials and
their thermal parameter values (cf. figure 9.5).
A personal user material may be added by mean of the command Add 1 material. Then the user has to give the
acronym and the value of the thermal characteristics of this new material. The modification of the library due
to this new material introduction will be validated or cancelled when leaving the Library form.
It is possible to modify the name and/or the parameter values of a personal user material. For selecting the
user material to modify, use a left mouse-click on the material to be modified, then the content of the
different cells may be changed.
The command Remove 1 material may be used for removing a personnel user material from the Frost-thaw
material library.
10 Launching computation
and editing results
In the first case (command Alize-Gel computation) the data used for the computation are obligatorily the data
shown on the main screen (pavement structure data) and specified by the selected option of the frame
Complementary data (initial and boundary temperature conditions). This first launching computation mode is
recommended in the most cases. In fact this first launching mode is a short-cut making easier and faster the
launching of the computation and the editing of its results.
In the second case (command Computation/Alize-gel of the main menu bar) the sheet Alize-Frost-thaw module,
launching the computation is opened (figure 10.2). Then the following choices are proposed:
Option 1 Current data (Frost-thaw main form) : then the data used for the computation are those
shown on the main screen (pavement structure data) and specified by the selected option of the
frame Complementary data (initial and boundary temperature conditions). In fact this launching
computation way (option 1) is quite equivalent to the previous one (short-cut command Alize-Gel
computation of the main Alize-Frost thaw sheet)
Option 2 Input file Alize-Gel (…dag) : then the data are read in pre-existing data file (extension .dag).
In this second case, the data relating to the initial and boundary temperature conditions are
obligatorily those specified by the Lcpc-Setra Design Guide (1994).
Figure 10.2 : Launching the Alize-Frost thaw computation from the command
Computation/Alize-Gel of the main menu bar
As soon as the computation are finished, the sheet Alize-Frost-thaw module - Computation results is
automatically open. The results are presented in 5 graphs, which can be individually shown on screen according
to choice done by the user.
Graph 1 :
Graph 1 represents the variations of surface pavement frost index (Is) and atmospheric frost index Iatm, as a
function of the frost quantity Qpf transmitted at the top Zpf of the pavement foundation. This first graph also
presents the variation curve Iatm=f(Qpf), resulting from the simplified method of calculation proposed by the
Lcpc-Setra Design guide (1994). This last curve is presented on graph 1, providing that the two following
condition are simultaneously verified:
All the materials making up the simulated pavement (including soil) are standard materials of the
Frost-thaw material library (materials having the status « user », cf. §9.6.1) In the opposite case, the
simplified Iatm=f(Qpf) cannot be applied, the linearization coefficients of the simplified relationship
being not defined for the not-standard materials.
The option Recall the results of the simplified method must have been previously selected. This option
is accessible through the Customize/Personal options of the main menu bar (cf. figure 10.3).
In case the frost quantity allowable by the pavement foundation (Qpf-adm) has been previously determined
(cf. §9.5), the atmospheric frost index Iatm associated to Qpf is automatically calculated and displayed on
screen. The comparison between this Iatm value and the atmospheric reference index IR is usually considered
by the user, in order to express a decision about the admissibility of the modelized pavement structure for
what concerns frost-thaw behaviour. Figure 10.4 presents an example of computation result curves IS = f(Qpf)
and Iatm = f(Qpf).
Graph 2 :
Graph 2 represents the evolution curves of frost quantity at interface between layers interfaces and pavement
surface level, versus the time (cf. example shown on the figure 10.5).
Graph 3 :
Graph 3 represents the vertical temperature profiles calculated at 6 value of the time, uniformly distributed
over the whole computation period (cf. example shown on figure 10.6).
Graph 4 :
Graph 4 represents the evolution curves of temperature at interface between layers and at pavement surface
level, versus the time (cf. example shown on the figure 10.7).
Graph 5 :
Graph 5 represents the evolution curve(s) of the frost front depth(s), versus the time (cf. example shown on
the figure 10.8).
The whole results of the current computation (graphs 1 to 5) may be stored in a temporary internal buffer,
using the button Memo Cj of the Alize-Frost results screen (j is the number of the concerned computation,
which is incremented by 1 at each launching of a new computation). Then it will possible to automatically
recall these stored results and superpose them to the results of further computations, using the button Remind
Cj of the Alize-Frost results screen. By this way, a direct comparison between the results of two computations
using different data may be directly displayed on screen, as shown on figure 10.9. The contain of the internal
buffer (results Cj) may be replaced at any moment by the results of the current computation (Ck), using the
button Memo Ck of the Alize-Frost results screen.
The sheet French frost indexes may be opened from the sheet Alize-Frost computation results. It displays the
frost reference indexes IR of the main French meteorological stations. This consultation may also be done from
the main form Alize Frost-thaw structure definition (cf. figure 10.10).
Lastly, the commands Print and Save of the sheet Alize-Frost computation results activate the edition of the
computation results, or their saving on Ascii file. The tables to print or save and the printing or output format
are defined thought different option boxes.
The figure 10.11 presents an example of printing of Alize Frost-thaw computation results.
Figure 10.3: Choice of the option Recall the results of the simplified method – or not – for the computation
of the relation Iatm = f(Qtpf), accessible from the Customize/Personal option of the main menu bar
Figure 10.4: Results of Alize-Frost module computation – Variation of the surface frost indexes
with the frost quantity transmitted at the top of the subgrade (example)
Figure 10.5: Results of Alize-Frost module computation – Evolution curves with time of the frost
quantities at interfaces and at the surface of the pavement (example)t
Figure 10.7: Results of Alize-Frost module computation – Evolution curves with time of the temperature at
interfaces and at the surface of the pavement (example)
Figure 10.9: Automatic superposition of the results curves of two different computations
by mean of the functions Memo Cj and Remind Cj
Figure 10.10: Help sheet Reference indexes of the main French meteorological stations,
according to the 1998 Catalogue
The pavement structure is simulated by a linear multilayer structure, according to the principles and the
hypothesis of the rational pavement design method implemented in the Alize software. The thickness, the
Poisson’s coefficient of each layer and the interface condition between layers are known data. The unknown
parameters of the problem to solve are the Young moduli of the materials. They are determined by a manual or
numerical procedure, base on the optimization of the adjustment between the measured deflections and the
deflections calculated by the Alize elastic-multilayer model. In common practice, the back-calculated moduli
may contribute to the evaluation of the actual damage of an in-service pavement. They can also be used for
determining the overlay of a damaged pavement structure.
Data concerning the pavement structure: number of layers, thicknesses and Poisson’s coefficients of
each layer, interface condition between layers. Initial values of the unknown elasticity moduli are
also needed. Depending on the type of back-calculation to performed, these initial module values are
used for initializing the numerical adjustment algorithm (Newton Raphson type), or to built the
computation grid, as detailed in $11.2
Data concerning the pavement loading: type of load (single or dual), radius of the print-foot, vertical
force or uniform pressure applied, and dual-wheel spacing in the case of dual-wheel.
Deflection bowl defined by the deflection measured at several points (up to 13) located at different
distances from the center of the load.
Back-calculation based on the automatic optimization of the adjustment between the measured and
computed deflection bowl. Then this problem of over-dimensioned system of linear equations is
solved by mean of an optimized Newton-Raphson algorithm implemented in the Alize software.
Back-calculation according to a grid of moduli values defined by the user. For each layer, the series
of module values to considered are predefined (arithmetical suite) by mean of its minimal value, its
arithmetic ratio, and the number of modulus values to consider. In this case, the program performs
the direct computation of the deflection basin associated to each module combination. The results
are finally classified in the order of decreasing adjustment precision.
Computation of the surface moduli. For a given point of deflection measure, the surface module is
defined as the Young module of the semi-infinite linear-elastic model (Boussinesq one’s) leading to
the same deflection as the measured one. The curve of surface moduli may bring some contribution
for the evaluation of the linear or not-linear behavior of the pavement, or even the substratum
depth.
In operational application, the most frequently mode of modulus back-computation is the Newton-Raphson
automatical procedure. The other two modes are generally reserved for specific application.
The access to the Alize-Back calculation module is done by mean of the command ”Calculation/Alize-Back
calculation” of the main menu bar (figure 11.1).
Pavement structure
Exportation of the pavement structure defined in the main window of the program Alize (cf. Part 4 of
the present manual). This structure is automatically transferred in the Back-calculation window as
soon as the Back-calculation module is activated.
Or using the command “New data” of the menu bar of the window Data for Back-calculation (cf.
figure 11.2). By mean of the command “Open data” of the same menu bar, it is also possible to load a
predefined structure saved in back-calculation data files (extension .mwd).
Figure 11.2: Building a pavement structure for back-calculation, from the window “Data for back-calculation”
The number of layers, thicknesses and Poisson’s coefficients, interface conditions are the basic data given by
the user. The Young modulus of a given layer may be a fixed data (no back-calculation adjustment for this
layer), or unknown parameter to be determined by mean of back-calculation.
Figure 11.3: Pavement structure data for back-calculation – option Surface moduli (example)
Figure 11.4: Pavement structure data for back-calculation – option Newton-Raphson algorithm (example)
Figure 11.5: Pavement structure data for back-calculation – option Grid (example)
In the first case (fixed module values), the material library may help to defined the fixed Young modulus of
standardized of personal materials.
In the second case (unknown Young modulus to be back-calculated), initial modulus values are requested:
If the Newton-Raphson option is selected: the initial values of Young modulus may have a significant
effect on the results of the back-calculation. Moreover, the numerical solution to the optimization
problem is not unique. Therefore, it is recommended to introduce initial values as close as possible to
the results of the back-calculation, or close lower values, using a manual iterative approach if
necessary.
If the grid option is selected, the minimal value of the modulus ranges for each material as to be
entered.
Figure 11.3 presents an example of data in the case of Surface moduli back-calulation. Figure 11.4 presents an
example in the case of Newton-Raphson back-calculation. Figure 11.5 presents an example in the case of grid
back-calculation.
In the cases of Newton-Raphson and Grid computation, the tick-box associated to each layer is used in order to
specify whether the modulus of this layer has to be back-calulated or not. If not (un-ticked box), the Young
modulus remains equal to the initial value given in the Structure data tab.
In the case of Newton-Raphson bac-kcalculation, the tick-box “Free/dependant” associated to each layer
permits to fix a dependency relationship between the modulus of the layer, and the modulus of the layer
directly below. If “Dependant” condition is selected, then the ration of these two moduli will be remain
constant, equal to the ration of their initial values defined in the Structure data tab.
Loading data:
The load used for back-calculation is defined by mean of the tab ”Load data” (figure 11.6). It may consist in a
single circular uniform vertical loading applied by a wheel or the plate of the FWD or HWD test device. It also
may consist in the dual wheel of the Deflectometer or the Curviameter.
Deflection data:
The measured deflections are entered by mean of the “Measured deflection” tab (figure 11.7).
Multi-deflection bowls data. This second mode for defining the deflection basins for back-calculation
makes easier the sequential treatment of a more or less great number of deflection measurements
along a survey itinerary. The opening of a “Multimeasure” file is done by mean of the “Multimeasure”
button in the tab “Measures data”. The arrangement of the multi-measure deflection Ascii files is
detailed in figure 11.8, and described bemow.
The multimeasure data file may include up to 512 columns corresponding to a maximum number of 511
deflection basins. Such a file includes the following parameters:
Title and identifier of the deflection series concerned (facultative), without limitation of the
number of lines.
Reminder about the key-words used by Multimeasure data files (facultative).
Key-word EpaisseurJ= (facultative): It allows to take into account the thicknesses of the
layer no J along the whole itinerary, which may be different from one deflection station to
another, and also different from the single thickness defined by the “Structure data” tab.
Unit = m.
Key-word ModuleJ= (facultative): It allows to take into account the values of fixed modulus
of the layer no J (not concerned by the backcalcultion adjustment) which may be different
from one deflection station to another, and also different from single modulus value defined
by the “Structure data” tab. Unit = MPa.
Key-word Poids= (facultative): It allows to take into account the values of the load applied
to the pavement, which may be different from one deflection station to another, and also
different from the single load value defined by the “Load data” tab. Unit = MN.
Key-word Unit= (facultative): It concerns the unit used for the measured deflection values.
Default unity is mm/100. The key-word Unit=micron specifies that the micrometer (µm) is
used instead of mm/100.
Positioning of the measure deflection points (obligatory), eg. Geophone distances from the
center of the load in the case of FWD measurement. These data are placed in the first
column named “X”. Unit = m.
Deflection basin in arranged in vertical columns (obligatory), up to 511 basins. Unit=mm/100,
unless the key-word Unit=micron is declared before.
Key-word Fin fichier for marking the end of the record.
The vertical scroll-cursor may be used in order to select one of the deflection bowl, and display the
corresponding curve on the graph. The tick-box “Echel” allows keeping unchanged the scales of the
graph, whatever the selected deflection bowl in the series of deflection basins.
Figure 11.9: Automatic importation of Dynatest or Carl-Bro FWD/HWD measure files - Example
Remark :
If a FWD or HWD Dynatest measure file or a FWD or HWD Carl-Bro measure file (Ascii format) is opened by
mean of the “Multimeasure” button (instead of a standard multimeasure file), then the configuration of the
FWD or HWD test device and the components of the measure are automatically detected and analyzed by the
program (cf. figure 11.9).
According to the choices selected by the user (force levels and drop nos), the corresponding Alize standard
multimeasure deflection file is automatically built.
Figure 11.10 presents an example of the back-calculation results, in the case of moduli surface computation.
Figure 11.11 presents another example, in the case of Newton-Raphson computation, and figure 11.12 in the
case of grid computation.
Figure 11.10: Results of Alize-Back calculation – Option surface moduli computation (example)
Back-calculated moduli : they are presented in the table showing the pavement structure, in the
column named “Esurf” in the case of surface moduli computation, or “Current E” in the case of
Newton-Raphson and Grid computation. The modulus values displayed in this table correspond to the
options and parameters selected by the user. In the case of Newton-Raphson and Grid computation,
these displayed values are depending on the scroll-cursor setting, which defines the no of iteration to
display. The default display of results is the moduli set giving the best adjustment between measured
and computed deflections, which is the first in the order of decreasing adjustment precision.
Option Surface moduli : the graph shows the surface moduli as a function of the distance between the
corresponding deflection sensor and the center of the load.
Option NewtonR and Grid computation: the graph shows both the measured deflection basin, and the
calculated one depending on the user’s setting by mean of the scroll-cursor.
Saving on file of the results of back-calculation is done in Ascii format, by mean of the command
“Save”, according to the display parameters fixed by the user.
The command “Filter” is illustrated by figure 11.13. It allows to automatically eject from the
complete retro-calculated results, the sets of moduli not included within the range specified by the
minimal-maximal limits given for the concerned layers.
Figure 11.13: Window « Results of back-calculation – Definition of the filter parameters (facultative)
The role of the buffer buttons M1, M2 and M3 is illustrated by figure 11 .14. Clicking the button Mj
launches the storage of the selected deflection profile (computed) in the internal buffer Bj. Then
ticking the tick-box Rj displays the drawing of the memorized basin on the graph. This function may
be useful in order to directly evaluate the effect of variable computation hypothesis – as layer
thicknesses, interface conditions, fixed modulus values … - on the results of the back-calculation.
Moreover, the general rational design approach has been the focus of a variety of transpositions, aimed at its
adaptation to other types of design: differing traffic loads (e.g. industrial carriers, airport loads), other design
problem areas (e.g. reinforcement of damaged pavements). The Alize software program is also applicable to
these specific uses, without its implementation conditions being necessarily, in this case, entirely contained
within an official method.
In any event, this software facilitates the practical running of numerical computations required for designing
pavement structures. The design of pavement structures extends far beyond these numerical aspects and still
lies within the domain of pavement mechanics specialists, whose expertise cannot be replaced despite the
level of effectiveness this software is able to offer.
This choice depends upon the framework and specificities of the project under examination. The pavement
structure may be of the flexible type, the thick bituminous type, semi-rigid, rigid, mixed or inverse. Each type
of structure is associated with the possible use of materials, themselves divided into five families: hydrocarbon
materials, materials treated with hydraulic binders, concretes, untreated materials, and soils.
2.2.1 - for treated materials (hydrocarbon, treated with hydraulic binders, and
concretes): The model adopted is that of fatigue-induced damage submitted to repeated
tensile stresses as a result of bending exerted by traffic (the Wöhler model). The allowable
stress is then a decreasing function of the number of loadings applied by the traffic and
may be expressed in the following general form:
Stadm = AxPxNb
where:
A and b = descriptive parameters of the material's fatigue behavior, as generated from
fatigue tests on laboratory specimens (b < 0);
P= additional coefficients for incorporating various factors, which have not been
reproduced by raw results from the laboratory fatigue tests;
N= number of rolling load cycles (i.e. cumulative traffic) over the entire service
life cycle ascribed to the pavement.
2.2.2 - for untreated materials (untreated gravels and soils): The model adopted is that
of damage caused by the accumulation of irreversible plastic strains (rutting), as generated
from the vertical compressive stresses exerted by the traffic. The expression of allowable
vertical compressive stresses does not stem directly from laboratory tests, but rather from
more empirical considerations:
Szadm = AxNb
where:
A and b = parameters of an empirical origin, independent of the specific untreated
material being examined, yet in practice varying as a function of both the type
of pavement and the intensity of traffic (new vs. reinforced pavement, light vs.
heavy traffic, etc.).
N= number of rolling load cycles – same as for Section 2.2.1 above.
The stresses created by rolling loads within the various material layers are calculated using the Burmister
linear elastic multilayer model (see Section 3.1), with Smaxj designating the highest stress from the standpoint
of the damage criteria defined above and with support provided by material layer n°j.
The thickness determination of each material layer then relies upon verifying the non-breaking criterion of the
material under investigation throughout the pavement service life, i.e.:
Smaxj ≤ Sadmj
where:
Smaxj: maximum stress generated by the traffic within material layer n°j, as stipulated by
the damage mode associated with the type of material being considered (bending-
induced horizontal tension or vertical compression, see Section 2.2);
The ultimate determination of thicknesses Hj of the various layers of materials constituting the pavement
structure does not, in general, lead to a direct solution (absence of explicit solution {Hj = f(N)}). In most cases,
this determination is performed by means of an inverse iterative approach, as diagrammed in Figure A1 below.
The numerical solution to the problem presented is moreover not at all unique. The final choice also depends
on both technological considerations (e.g. minimum and maximum layer thicknesses as a function of the
material, traffic classification, etc.) and empirical considerations (e.g. thickness of the surface layers).
▼
Output:
Continuation of
the project, other
verifications: e.g.
freezing conditions
This model relies upon a mechanical description of the structure using a semi-infinite foundation, generated by
superimposing constant thickness material layers, featuring isotropic linear elastic behavior. The number of
descriptive parameters for the mechanical behavior of each material can thus be reduced to two: strain
modulus E, and Poisson's ratio .
Each interface between adjacent layers may be set as bonded, sliding or semi-bonded.
The loading applied by rolling loads at the pavement surface is represented by a set of circular disks, each
loaded by a static and uniform vertical pressure.
Any combination of circular elementary loads of this type is indeed possible, a feature that enables
reproducing the profiles of axles, bogie units, a vehicle or of the widest array of vehicles.
Within more commonplace design operations, a unique loading, designated as the reference loading, is
typically applied at the model surface; this reference loading is specific to the context in which the rational
design method is being applied.
For example, a 65-kN twinning (3.25 kN per wheel, with a 6.62 MPa contact pressure) constitutes, except in
special cases, the reference design load for pavement structures on France's road and highway networks.
In some airport pavement design applications, the reference load will be taken as the isolated 30-kN wheel
load, etc.
P = 3.25 kN
Q LCPC-SETRA Method
3r Reference load:
= r R = 0,125 m
P = 3,25 kN
0 Q = 0,662 MPa
.
6
6
2 E1, 1 h1: Surface layer
E2, 2 h2: Base course
M
P E3, 3 h3: Foundation layer
a
E4, 4 h4: Subgrade
The computation model may potentially provide insight into the full stress and strain tensors, as well as the
three displacement components, at any point of the structure. According to common practice, the model
results strictly necessary for achieving the design objective are reduced to the maximum stress Smaxj
supported by each material layer of index j (from j = 1 for the surface layer to j = nc for the soil support).
Depending on the type of material involved, Smaxj serves to represent the following magnitudes:
On a large part of the French territory the taking into account of frost-thaw periods is a determining factor in
the design of pavements for low and mean traffic, so long as the maintaining of the traffic of heavy trucks
during the thaw periods is pursued.
Until the end of the 70s, the frost-thaw verification of pavement design has been done on the basis of the
empirical relationship Hng > k x H. Hng was the total thickness of frost resistant material of the structure, and
H was the maximal frost depth against which protection was requested, with k close to 80%. In the middle of
the 70s, the Lcpc and the Setra have undertaken the development of a new method for the frost-thaw
verification of pavements, based on a rational approach in continuity with the rational mechanical design of
pavement thicknesses.
This new method can be distinguished from the empirical previous one, by the following points:
A better taking into account of the frost susceptibility of soils, based on the results of the swelling
test in laboratory ;
A better characterization of the winter rigour against which protection is searched (reference winter),
by mean of the winter frost index instead of the frost depth ;
A better taking into account of the respective effects of frost resistance of soil and the mechanical
resistance of the pavement ;
A verification method applicable to road contexts different from those of the national French
network, which was impossible by the previous empirical method.
This rational frost-thaw verification of pavement structure has been implemented in the module Alize-Lcpc Gel
(also referred as Alize-Lcpc Frost-thaw module).
As the thaw stage occurs the frost resistance soils immediately recover their initial water content, with a very
few exceptions. Therefore the normal mechanical properties are also immediately recovered. On the contrary,
the water content of frost susceptible soils is much higher than initial water content (before the beginning of
the frost period), because of the thawing of ice cells. At this moment, the water content may exceed the
liquidity limit of the soil. Therefore, the bearing capacity of a frost susceptible soil is considerably reduced,
indeed even neutralized, during this thaw period.
The figure A2-1 shows an example of severe damages created by the thaw after an intense frost period.
Appendix A2: Principles of the Frost-thaw verification
To opt for another structure (soil treatment, thicker pavement structure or foundation layers, used of
other materials, …)
To accept possible traffic interruptions during the thaw period.
Frost Quantity :
It is the square root of the frost index (unity [°C.day]1/2). The following frost quantities are used:
- Qatm atmospheric frost quantity;
- Qs pavement surface frost quantity;
- Qt (or Qpf) : frost quantity transmitted at the bottom of the pavement structure.
The figure A2-3 represents the general progress of frost-thaw verification studies. The verification of the frost-
thaw behaviour involves comparing:
the atmospheric frost index chosen as reference, IR, that characterizes the severity of winter against
which it has been chosen to protect the pavement,
and the atmospheric frost index that the pavement is able to withstand, which is called the allowed
frost index IA. This index is evaluated according to the frost susceptibility of the subgrade, the
thermal protection and mechanical function fulfilled by the pavement.
The pavement structure is designed so that the allowed frost index IA of the pavement is higher than the
reference frost index IR. If the frost index happens to be larger than IA, the structure has to be modified (other
higher thicknesses, or materials, or pavement foundation …), or thaw barriers may have to be installed to
protect the pavement during thaw periods.
The 1998 Catalogue also indicates the frost indexes of the exceptional winter and the rigorous not exceptional
winters recorded by the main meteorological stations on the 1951-1997 period.
Appendix A2: Principles of the Frost-thaw verification
Logiciel
Pavement owner’s specifications towards frost-thaw
verification ?
Determination of the
allowable atmospheric frost
index IA = f(Qpf)
Figure A2-4 : Laboratory equipment for the frost swelling test 0/5mm of the
Nancy’s Public works regional laboratory
Stages for the computation of the allowable frost index of the pavement
The determination of the allowable frost index IA includes 5 successive stages, from the subgrade to the
pavement surface, as shown by the figure A2-5.
IA = IS/0.70 + 10 (°Cxj)
step 5 (IA determination)
Stage 1 : Analysis of the frost susceptibility of the soil supporting the pavement
The frost quantity Qng, which represents the thermal protection due to the frost resistant materials of
the pavement foundation layer and soil. Qng is evaluated from the nature and the thickness of these
frost resistant materials.
It is accepted that the frost can penetrate into the frost-susceptible subgrade. This frost penetration (over the
thickness e) corresponding to the additional mechanical frost quantity Qm. The decrease in bearing capacity
due to this penetration creates an increase in the maximal tensile stress S2 in the pavement structure induced
by the traffic, in comparison with the maximal stress S1 over the normal situation other than the frost-thaw
season. The increase of the stress S2 is limited to about 5%. During the thawing phase, the Young modulus of
the subgrade subject to thaw is taken as equal to a tenth of the normal modulus in the situation other than the
frost/thaw period, ie E(Pf)/10. Iterative computations by mean of the Alize-mechanical module lead to the
determination of the thickness e (penetration thickness of the frost front in the frost susceptible part of the
subgrade). Then the quantity Qm is calculated using the relationship: Qm= e/10 (unities: sqare(°C.day) and
cm).
Stage 3 : Determination of the frost quantity allowable by the pavement foundation layer
The quantity of frost considered as allowable at the top of the capping layer is obtained as follows:
Qpf = Qg + Qng +Qm
Two approaches are proposed by the Pavement design technical guide (1994) for the computation of the
thermal protection due to the pavement structure:
either a thermal calculation to work out the propagation of frost in the pavement. Then the frost
propagation is calculated using the Alize – Frost-thaw module. This first approach is the more accurate
and is recommended in the normal application of the French pavement design method.
Or using simplified equation (manual application).
The frost propagation through the pavement leads to determine the relationship Qpf = f(IS) by mean
of the Alize – Frost-thaw module. It implements the Fourier’s model in order to solve the heat
equation propagation in multilayered mixed unfrosted-frosted structure by mean of a finite
differences algorithm. All the hypotheses taken into account by the Alize – Frost-thaw module are
defined by the Lcpc-Setra Pavement design guide Lcpc-Setra (1994). This includes the geometrical
definition of the pavement structure, its thermal characteristics, and the definition of the vertical
initial temperature profile and the temperature boundary condition at the top and the bottom of the
pavement structure + subgrade.
The simplified method is less accurate than the previous frost propagation method. It is based on
manual calculation using a linearized relationship between the frost quantities Qpf and Qs. For a given
layer, the two coefficients of this linear relationship depend on the nature of the material and the
thickness of the layer, as described by the Lcpc-Setra Pavement design technical guide.
The relationship between Qpf and IS established in the previous stage 4 leads to determine the value of IS =
finverse(Qpf). The atmospheric frost index corresponding to the surface frost index IS is calculated by mean of a
empirical relationship between the indexes IA and IS. This relation accounts for the convection and radiation
phenomena at the pavement surface. At average altitude, with little to average sunshine and an atmospheric
frost index not exceeding 210 °C.day, the Lcpc-Setra Pavement design guide recommends the following
equation: IA = IS/0.7 + 10 (unity: °C.day).
For all the other cases (very severe frost, much sunshine) an individual study will be need.