# Behaviour of an interface between pavement layers obtained using Digital Image Correlation

T. Attia · H. Di Benedetto · C. Sauzéat · S. Pouget

Received: 5 June 2020 / Accepted: 7 January 2021

© The Author(s) 2021, corrected publication 2022

**Abstract** To improve pavement design methods, the behaviour of the interfaces between pavement layers must be precisely taken into account. This paper describes the thermomechanical behaviour of an interface in the small strain domain. Thanks to the innovative 2T3C Hollow Cylinder Apparatus (2T3C HCA), the behaviour of an interface was observed directly with 3D Digital Image Correlation (3D DIC). The studied samples were constituted of two layers of bituminous mixtures with a bitumen emulsion at the interface. From cyclic tests performed at four different frequencies (from 0.01 to 0.3 Hz) and four different temperatures (from 10 to 40 °C), the complex interface stiffnesses in the tension–compression mode and in the shear mode were obtained. The time–temperature superposition principle was verified for the abovementioned stiffnesses. Knowing a 2S2P1D (2 Springs, 2 Parabolic elements, 1 Dashpot) model of the bitumen of the tack coat, the norms of the complex interface stiffnesses were modelled with 2S2P1D models by using a SHStS (Shift, Homothety, Shift and time Shift) transformation.

**Keywords** Pavement · Interface · Bituminous mixture · Tack coat · 2S2P1D · 3D digital image correlation

## 1 Introduction

When designing pavements, layers made of bituminous mixtures are usually assumed to be perfectly bonded in the mechanical calculations (i.e. continuity of displacements at the interface). In this case, the maximum strain is obtained at the bottom of the structure and strain in the top layer is reduced [1] which is advantageous for pavement durability [2]. However, premature degradations of pavements linked with bond failure have been observed. Understanding the mechanical behaviour of interfaces between pavement layers is necessary for designing more durable roads.

The experimental research on interfaces behaviour has focused on bonding strength measurements by testing the interface failure, for instance with the guillotine test [3–5]. The failure tests can be monotonic [6] or cyclic for fatigue characterisation [7–11]. Multiple factors such as temperature [12], dosage of tack coat or type of bituminous mixtures layers were investigated [13]. These tests were either done at high strain levels for monotonic failure tests or with an important number of cycles for fatigue tests. But it is also important to know the interface behaviour for a

<footnote>T. Attia (✉) · H. Di Benedetto · C. Sauzéat
LTDS (CNRS 5513), University of Lyon/ENTPE, 3 rue
Maurice Audin, 69518 Vaulx-en-Velin, France
e-mail: thomas.attia@gmail.com</footnote>

<footnote>S. Pouget
Eiffage Infrastructures, 8 rue du Dauphiné, 69960 Corbas,
France</footnote>
few number of cycles at small strain amplitude. It allows computing more accurately the strain at the bottom of the structure due to a single vehicle which is an important piece of information in pavement design. Up to the knowledge of the author, very few authors attempted such a characterisation [14–16].

Concerning the mechanical modelling of interfaces, instead of the simple assumption of a perfect bond at the interface, an interface stiffness can be defined [1]. It links the stress applied to the relative displacement between the boundaries of the layers. This approach is easy to implement in a numerical simulation of a pavement. Interfaces can also be modelled as layers. Usually, a tack coat made of bitumen emulsion is applied to guarantee the bond at the interface between the layers. For a small number of cycles and for small strain amplitudes, bitumen presents a linear viscoelastic behaviour. It has been shown from in situ measurements that interfaces with a tack coat are better modelled by a homogeneous viscoelastic layer rather than an elastic layer [17]. Laboratory experiments also suggest that the behaviour of interfaces for small strain, assuming they have a thickness or not, is viscoelastic. Freire et al. [15] modelled the normal stiffness of the interface using the 2S2P1D model, which is a linear viscoelastic model, when Hristov et al. [16] used a sigmoid function to model the norm of the shear interface stiffness. However, it is not clear whether a linear viscoelastic domain exists for the interface since very few studies focus on the influence of the amplitude of loading. Isailović et al. [10] did not observe a linear viscoelastic domain since the shear interface stiffness changed with the amplitude of loading. Yet, their data was limited to one frequency and two temperatures, and the tested amplitudes were quite high.

The objective of this paper is to characterise the thermomechanical behaviour of an interface with a tack coat in the small strain domain using the innovative apparatus 2T3C Hollow Cylinder Apparatus (2T3C HCA). This work focuses on highlighting the relation between viscoelastic properties of the tack coat bitumen and the interface behaviour. It is part of a research project (chair between University of Lyon/ ENTPE and the company Eiffage Infrastructures) dedicated to pavement materials and structures.

## 2 2T3C hollow cylinder apparatus (2T3C HCA)

The 2T3C HCA was developed at the University of Lyon/ENTPE to study interfaces between bituminous mixtures [18]. The tested samples are hollow cylinders composed of two layers of bituminous mixtures as illustrated in Fig. 1. Their dimensions are: external radius of 86 mm, internal radius of 61 mm and height of 125 mm. The dimensions were chosen to ensure rather homogeneous stress and strain fields in the hollow cylinder. Axial and rotation loadings can be applied to the sample independently using a servohydraulic press. Various loadings (cyclic or monotonic, from small deformation to failure) can be applied. The samples are placed in a climate chamber allowing to control temperature during the tests. PT100 sensors measure sample temperature on the inside surface and on the outside surface of the sample.

Load cells embedded in the servohydraulic press measure the axial force and the torque, up to 100 kN for axial force and up to 2 kN.m for torque. Four non-contact sensors (eddy currents sensors) are attached to the top of the samples (Fig. 2). Their targets are fixed to the bottom of the sample. The non-contact sensors measure the displacement between the top and the bottom of the sample in the axial and rotation directions with an accuracy of 0.1 µm on a 1000 µm range. Two pairs on opposite sides of the sample are used to calculate the average vertical displacement and the average rotation displacement. The press can use the sensors measurements to control precisely the applied loadings.

In addition to these global displacement measurements, the 3D Digital Image Correlation (3D DIC) is used to obtain the three-dimensional displacement field on the surface of the sample. The principle of DIC is to compare the pictures of a surface with a random pattern before and after transformation. A subset of the reference picture (captured before the transformation)

<figure><img src="image_1.png" /><figcaption>Fig. 1 Hollow cylinder sample of the 2T3C HCA</figcaption></figure>
<figure><img src="image_2.png" /><figcaption>**Fig. 2** Picture and scheme of the 2T3C HCA with measuring systems</figcaption></figure>

is chosen around a point of the surface where displacements are searched. Then a correlation algorithm finds the deformed subset after transformation in the following pictures and deduces the displacement of the point [19]. The subset is a window of 25 × 25 pixels in this study. For 3D DIC two pairs of cameras are used on opposite sides of the sample as shown in Fig. 2, with an angle of 30° between cameras in each pair. The cameras have a maximum capture frequency of 10 Hz. They are placed outside the climate chamber (which has glass windows on its sides) at 600 mm from the sample surface. Considering the digital sensors characteristics, one pixel on a picture represents about 150 µm of material. The accuracy of the displacement measurements with DIC is about 1/100th of a pixel which corresponds to 1.5 µm in our setup. The application of DIC to study interfaces is scarce [20, 21] and for 3D DIC it is even more rare.

From the displacement field, a method developed at the University of Lyon/ENTPE [18] is used to obtain the strain tensor components $\varepsilon_{zz}$ and $\varepsilon_{0z}$ in the layers and the displacement gaps at the interface in the vertical direction $\Delta u_z$ and for the rotation $\Delta u_\theta$, ($r, \theta, z$) being the cylindrical coordinates. The example of calculation of $\varepsilon_{zz}$ and $\Delta u_z$ is illustrated in Fig. 3. Horizontal strips (100 × 4 mm) are drawn on the digital images. The number of strips depends on the sample dimensions. For the samples tested in the present study, there are 12 strips in each layer. The displacement component $u_z$ is averaged on each strip and is affected to the coordinate Z of the centre of the strip. By plotting the linear regression in the $Z-u_z$ diagram, the slopes of the two lines (see Fig. 3) give $\varepsilon_{zz}$ in each layer. The assumption of a homogeneous

strain field is verified experimentally. The regression lines do not intersect at the interface and the gap between them at the level of the interface gives the vertical displacement gap, $\Delta u_z$. A similar procedure is used to get $\varepsilon_{0z}$ and $\Delta u_\theta$, based on the displacement related to rotation $u_\theta$ obtained with 3D DIC. The results from the two opposite sides of the samples are then averaged to have one value at each state of deformation.

# 3 Materials and procedure

## 3.1 Studied materials

In this paper, the studied double-layered specimens are composed of two bituminous mixtures layers and a tack coat at the interface. First, half of a slab of 600 × 400 × 75 mm³ of EME 20/14 is compacted using a French wheel compactor [22]. The EME 2 mixture is commonly used in France as a base course. Here it has a nominal maximum aggregate size (NMAS) of 14 mm, it contains 30% of Reclaimed Asphalt Pavement (RAP) and its binder is a bitumen with a 15/25 grade. The binder represents 5.6% of the total mix weight.

24 h after the compaction of the first layer, an emulsion of bitumen graded 160/220 is applied with a brush to act as a tack coat. The residual bitumen content is 450 g/m². The linear viscoelastic behaviour of the bitumen 160/220 used for the tack coat has been characterised with a Dynamic Shear Rheometer (DSR) at the Eiffage Infrastructures laboratory in Corbas. The results are used and discussed in Sect. 4.
<figure><img src="image_3.png" /><figcaption>**Fig. 3** Vertical displacement $u_z$ from 3D DIC (picture with a global extension of the sample of 12.5 µm) and illustration of the computation of $\varepsilon_{zz}$ in each asphalt layers and of the</figcaption></figure>

displacement gap $\Delta u_z$ at the interface from average values of $u_z$ obtained in each of the 24 strips of 10 cm × 4 mm

After another 24 h, the upper layer of the slab is compacted with the wheel compactor. It is a BBSG 3 0/10 mixture with a NMAS of 10 mm, a RAP content of 20% and a 50/70 bitumen. The binder content is 5.4% of the mix weight. The aggregate grading curves of both bituminous mixtures are presented in Fig. 4. The sample coring takes place a few days later starting with the inside coring of the hollow cylinder and finishing with the outside. Three hollow cylinder samples are cored from the slab and named as G8, D8 and M8.

<figure><img src="image_4.png" /><figcaption>**Fig. 4** Aggregate grading curves of bituminous mixtures</figcaption></figure>

## 3.2 Advanced complex modulus test with 2T3C HCA

To characterise the viscoelastic behaviour of the interface in the small strain domain, sinusoidal cycles are applied to the samples. For each of the 3 samples, a temperature/frequency sweep was conducted at 4 different temperatures (10, 20, 30, 40 °C) and 4 different frequencies (0.01, 0.03, 0.1, 0.3 Hz) as described in Fig. 5. For each couple of temperature and frequency, 5 axial cycles and 5 rotation cycles were performed.

For axial cycles, a global vertical displacement amplitude of 25 µm (measured with the non-contact sensors) was imposed. This would correspond to a global vertical strain of 200 µm/m for a hollow cylinder sample with the same dimensions but constituted of one homogeneous material. The vertical strain in the layers and the displacement gaps could only be obtained from the DIC analysis after the end of the mechanical testing. For axial cycles, the measured vertical strain amplitudes $\varepsilon_{zz\_0}$ were between 150 µm/m and 220 µm/m in the upper layer and between 100 and 170 µm/m in the lower layer. The vertical displacement gaps amplitudes at the interface were between 1 and 2 µm except for the test at 40 °C for the D8 sample where vertical displacement gaps were below 1 µm and could not be analysed. The latter
<figure><img src="image_5.png" /><figcaption>**Fig. 5** Temperature/frequency sweep for the advanced complex modulus test with the 2T3C HCA, tests are controlled in displacement mode with an amplitude calculated to get 200 µm/</figcaption></figure>

temperature was removed from the analysis. The amplitude values for strain and displacement gaps vary for each couple of temperature and frequency because they depend on the ratios between the layers' moduli and the interface stiffness.

For rotation cycles, a global rotation displacement amplitude of 90 µm at the level of the non-contact sensors was imposed. It would correspond to a global shear strain of 200 µm/m in an homogeneous sample of the same size. The measured shear strain amplitudes $\epsilon_{0z-0}$ in the bituminous mixtures were between 200 and 275 µm/m in the upper layer and between 100 and 150 µm/m in the lower layer. The rotation displacement gaps amplitudes were between 5 and 10 µm.

During the mechanical loading, 50 pictures per cycle were taken with each camera except for tests at 0.3 Hz where only 30 pictures per cycle could be taken. Cameras can take pictures at a maximum rate of 10 Hz, preventing from testing at higher frequencies than 0.3 Hz. The capture time of a picture is about 6 ms. An acquisition system allowed obtaining the load cells measurements at the exact same time than a picture was taken. At the end of the mechanical testing, the pictures were analysed with DIC as presented above. For each picture the displacement gap at the interface is obtained, $\Delta u_z$ for axial tests and

m in global strain ($\epsilon_{zz}$ for axial tests or $\epsilon_{0z}$ for rotation tests), as if the sample was constituted with one homogeneous material

$\Delta u_0$ for rotation tests. The stresses are calculated using the assumption that they are homogeneous in the sample. An example of displacement gap and stress measurement is presented on Fig. 6, where each dot represents a state of deformation, i.e. an instant when pictures are taken. For one sample, about 28,800 pictures are captured during an advanced complex modulus test.

### 3.3 Overview of the performed tests

The behaviour of one interface configuration (see Sect. 3.1) was evaluated using the advanced complex modulus test performed with the 2T3C HCA. Three 2T3C HCA samples were tested for this configuration. In addition, the bitumen used for the tack coat was tested using a shear complex modulus test performed with a DSR. One sample of bitumen was tested. All the tests performed in this study are summed up in Table 1.

### 3.4 2S2P1D (2 springs, 2 parabolic elements, 1 dashpot) model

For small strain amplitudes and few cycles applied, bitumen and bituminous mixtures present a linear
<figure><img src="image_6.png" /><figcaption>**Table 1** Overview of the performed tests</figcaption></figure>

<table><thead><tr><th>Tested material</th><th>Performed test</th><th>Experimental device</th><th>Samples</th></tr></thead><tbody><tr><td>Bilayered material: BBSG 3 on EME 2 with emulsion of pure bitumen 160/ 220 at the interface</td><td>Advanced complex modulus test</td><td>2T3C HCA</td><td>3</td></tr><tr><td>Pure bitumen 160/220</td><td>Shear complex modulus test</td><td>DSR</td><td>1</td></tr></tbody></table>

viscoelastic behaviour. The strain limit of the linear domain depends on the temperature and on the frequency. Orders of magnitude for these limits are 100 µm/m for mixtures [23, 24] and 10,000 µm/m for bitumen [25]. The fact that the interface behaviour is linear in the small strain domain has not been proved and it is not the purpose of this paper. The complex modulus is considered as a linear viscoelastic property obtained by conducting a sinusoidal loading. In the case of small non-linearity, an equivalent complex modulus can be considered as both strain and stress signals are close to sinuses.

For axial tests, when a sinusoidal normal stress with an amplitude of $\sigma_{zz\_0}$ is applied to the interface the vertical displacement gap is a sinus with the same frequency, an amplitude $\Delta u_{z\_0}$ and with a phase lag $\phi_{Kzz}$. The normal complex interface stiffness $K_{zz}^*$ is defined in Eq. 1.

$$K_{zz}^* = \frac{\sigma_{zz\_0}}{\Delta u_{z\_0}} e^{i\phi_{Kzz}} = |K_{zz}^*| e^{i\phi_{Kzz}} \quad (1)$$

For rotation tests, if a sinusoidal shear stress loading with an amplitude of $\tau_{\theta z\_0}$ is applied to the

interface the rotation displacement gap response is sinusoidal with the same frequency, an amplitude $\Delta u_{\theta\_0}$ and with a phase lag $\phi_{K\theta z}$. The shear complex interface stiffness $K_{\theta z}^*$ is introduced in Eq. 2.

$$K_{\theta z}^* = \frac{\tau_{\theta z\_0}}{\Delta u_{\theta\_0}} e^{i\phi_{K\theta z}} = |K_{\theta z}^*| e^{i\phi_{K\theta z}} \quad (2)$$

The 2S2P1D model developed at the University of Lyon/ENTPE is used to model the complex modulus dependency on the frequency for bituminous mixtures and binders [26]. The complex modulus expression for 2S2P1D model for the shear complex modulus $G_{2S2P1D}$, the normal complex interface stiffness $K_{zz\_2S2P1D}$ and the shear complex interface stiffness $K_{\theta z\_2S2P1D}$ are given in Eqs. 3, 4 and 5, respectively. There are 7 constants to the model: $G_{00}$ (respectively, $K_{zz\_00}$ and $K_{\theta z\_00}$ for the normal and shear complex interface stiffnesses) the static modulus which is the asymptotic value for low frequencies, $G_0$ (resp. $K_{zz\_0}$ and $K_{\theta z\_0}$) the glassy modulus which is the asymptotic value for high frequencies, $k$, $h$ and $\delta$ that are calibration parameters, $\tau$ a characteristic time depending on the temperature and $\beta$ a viscosity parameter.
This model is used in this paper to model the bitumen shear complex modulus $G^*$ and the complex interface stiffnesses $K^*_{zz}$ and $K^*_{\theta z}$. In Eqs. 3, 4 and 5 $\omega$ is the angular frequency.

$$G^*_{2S2P1D}(\omega) = G_{00} + \frac{G_0 - G_{00}}{1 + \delta(i\omega\tau)^{-k} + (i\omega\tau)^{-h} + (i\omega\beta\tau)^{-1}} \quad (3)$$

$$K^*_{zz\_2S2P1D}(\omega) = K_{zz\_00} + \frac{K_{zz\_0} - K_{zz\_00}}{1 + \delta(i\omega\tau)^{-k} + (i\omega\tau)^{-h} + (i\omega\beta\tau)^{-1}} \quad (4)$$

$$K^*_{\theta z\_2S2P1D}(\omega) = K_{\theta z\_00} + \frac{K_{\theta z\_0} - K_{\theta z\_00}}{1 + \delta(i\omega\tau)^{-k} + (i\omega\tau)^{-h} + (i\omega\beta\tau)^{-1}} \quad (5)$$

The dependency of the complex modulus on the frequency can be related to the dependency on the temperature. When a bituminous material respects the Time-Temperature Superposition Principle (TTSP), its behaviour is equivalent for low frequencies and high temperatures (and for high frequencies and low temperatures). At each temperature $T$, it is possible to define a shift factor $a_T$ in order to superimpose isothermal curves of the norm (or of the phase angle) of the complex modulus, forming a unique curve called the master curve. The master curve represents the material behaviour at an arbitrarily chosen reference temperature. The Williams-Landel-Ferry (WLF) formula presented in Eq. 6 links the shift factors $a_T$ to the temperature $T$ with regards to the reference temperature $T_{ref}$ using two constants $C_1$ and $C_2$.

$$\log(a_T) = \frac{-C_1(T - T_{ref})}{C_2 + T - T_{ref}} \quad (6)$$

The shift factors are determined for a reference temperature chosen amongst the experimental temperatures. In order to compare the results with other tests, it is necessary to obtain the shift factors and thus the constants $C_1$ and $C_2$ of the WLF formula for other reference temperatures. Equations 7 and 8 express the constants $C'_1$ and $C'_2$ for the WLF formula at a reference temperature $T'_{ref}$ from the WLF constants $C_1$ and $C_2$ at the reference temperature $T_{ref}$.

$$C'_1 = \frac{C_1 C_2}{C'_2} \quad (7)$$

$$C'_2 = C_2 + T'_{ref} - T_{ref} \quad (8)$$

### 3.5 SHStS (Shift, Homothety, Shift and time Shift) transformation

It has been shown that the viscoelastic properties of a bituminous mixture depend directly on the binder behaviour. The mixture complex modulus $G^*_{mix}$ can be derived from the bitumen modulus $G^*_{binder}$ using the SHStS transformation according to Eq. 9 for given angular frequency $\omega$ and temperature $T$ [27]. The glassy moduli $G_{0,mix}$ and $G_{0,binder}$ as wells as the static moduli $G_{00,mix}$ and $G_{00,binder}$ of both the mixture and the binder appear in Eq. 9, independently of any rheological model, along with the parameter $\alpha$ introduced for the time shift.

$$\frac{G^*_{mix}(\omega, T) - G_{00,mix}}{G_{0,mix} - G_{00,mix}} = \frac{G^*_{binder}(10^\alpha \omega, T) - G_{00,binder}}{G_{0,binder} - G_{00,binder}} \quad (9)$$

If the linear viscoelastic behaviour of both mixture and binder are modelled using 2S2P1D, then the SHStS transformation implies that the constants $k, h, \delta, \beta$ are the same in the mixture model and in the binder model. If the TTSP is verified for the binder, then it is also verified for the mixture and the shift factors $a_T$ are the same. In this case, the relation between the characteristic times $\tau$ of mixture and binder for any temperature $T$ is the Eq. 10.

$$\tau_{mix}(T) = 10^{\alpha} \tau_{binder}(T) \quad (10)$$

The SHStS transformation is used in Sect. 4 to obtain the norm of the complex interface stiffnesses from the shear complex modulus of the bitumen of the tack coat.

## 4 Results and discussion

### 4.1 Normal and shear complex interface stiffnesses

Shifting the isothermal curves, it was possible to build the master curves of the norm of the normal complex interface stiffness, Fig. 7, and of the norm of the shear
<figure><img src="image_7.png" /><figcaption>**Fig. 7** Master curves of the norm of the normal complex interface stiffness $|K_{*zz}|$ for the three samples (G8, D8, M8) at the reference temperature of 15 °C with 2S2P1D model</figcaption></figure>

complex interface stiffness, Fig. 8. In both cases, the TTSP is valid. The shift factors used to plot the master curves of the norm of the complex interface stiffnesses at the reference temperature 15 °C are presented in Fig. 9. They were obtained from the WLF formula, first obtained at the experimental reference temperature (one for each sample) and then expressed for the reference temperature 15 °C using Eqs. 7 and 8. The WLF constants of all samples were averaged to plot the WLF formula on Fig. 9. The constants at a reference temperature of 15 °C are $C_1 = 7.22$ and $C_2 = 112.00$.

A good repeatability was observed between the samples, especially for shear tests. Interfaces are consistently more rigid in the tension–compression

mode than in the shear mode at the same temperature and frequency.

<figure><img src="image_8.png" /><figcaption>**Fig. 9** Shift factors for interface and bitumen</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 8** Master curves of the norm of the shear complex interface stiffness $|K_{*\theta z}|$ for the three samples (G8, D8, M8) at the reference temperature of 15 °C with 2S2P1D model</figcaption></figure>
## 4.2 Interface behaviour modelling

The objective of this section is to model the interface behaviour with a 2S2P1D model from the bitumen behaviour.

From the results of the DSR testing, the shear complex modulus of the 160/220 bitumen was modelled with a 2S2P1D model. The parameters of this model are presented in Table 2, G00 was found to be nil as it is commonly seen for bitumen. To compare the bitumen behaviour and the interface behaviour, the norm of the bitumen shear complex modulus and the norm of the complex interface stiffnesses are normalised. The phase angle is not taken into account in this study. For the bitumen, the expression of the normalisation depends on the static modulus G00 and on the glassy modulus G0 that were estimated when calibrating the 2S2P1D model. This expression of the normalised norm of the shear complex modulus |G*|norm is presented in Eq. 11, taking into account that G00 is nil.

$$|G^*|_{\text{norm}} = \frac{|G^*|}{G_0} \qquad (11)$$

For the interface, the normalised norms |K*<sub>zz</sub>|norm and |K*<sub>θz</sub>|norm of the interface stiffnesses were obtained using Eqs. 12 and 13.

$$|K_{zz}^*|_{\text{norm}} = \frac{|K_{zz}^*| - K_{zz\_0}}{K_{zz\_0} - K_{zz\_0}} \qquad (12)$$

$$|K_{\theta z}^*|_{\text{norm}} = \frac{|K_{\theta z}^*| - K_{\theta z\_0}}{K_{\theta z\_0} - K_{\theta z\_0}} \qquad (13)$$

$K_{zz\_0}$ and $K_{zz\_00}$ for the interface are found with the master curves of the three samples on Fig. 7, whereas $K_{\theta z\_0}$ and $K_{\theta z\_00}$ are found with the master curves on Fig. 8. These constants are presented in Table 3. The equivalent frequencies of the master curves of the interface stiffnesses were then multiplied by 50,000 to superimpose the master curves of interface and bitumen on Fig. 10. This corresponds to a time shift

with a coefficient $\alpha = 4.69$ as defined in the SHStS transformation. The result of the transformation is plotted on Fig. 10.

On Fig. 10 there is a good match between the normalised shear modulus of the bitumen |G*|norm and the normalised norms |K*<sub>zz</sub>|norm and |K*<sub>θz</sub>|norm of the interface stiffnesses. The interface behaviour dependency on time and temperature is thus the same than for the bitumen of the tack coat. Based on this observation, the behaviour of the interface was modelled with a 2S2P1D model using $K_{zz\_0}$, $K_{zz\_00}$, $K_{\theta z\_0}$ and $K_{\theta z\_00}$ from Table 3, the parameters k, h, $\delta$, $\beta$ of 2S2P1D model of the bitumen (Table 2) and $\tau$ computed with Eq. 10 using $\alpha = 4.69$. These operations correspond to a SHStS transformation to obtain the interface behaviour from the bitumen behaviour. The parameters of the 2S2P1D models for the interface (one for the normal stiffness, one for the shear stiffness) are presented in Table 2.

2S2P1D models were plotted in Fig. 7 for the axial modulus and in Fig. 8 for the shear modulus. The models fit well the experimental data.

## 5 Conclusions

Cyclic tension–compression tests and cyclic rotation tests at small strain amplitude were performed on hollow cylinder samples. Three samples of the same bilayered material composed of two bituminous mixtures layers and a bitumen emulsion at the interface were tested. The tests were done at four different frequencies (0.01 to 0.3 Hz) and four different temperatures (from 10 to 40 °C). The viscoelastic behaviour of an interface between bituminous mixtures layers has been observed thanks to the Digital Image Correlation. This optical measurement technique made it possible to evaluate precisely the displacement gap at the interface.

The normal complex interface stiffness $K^*_{zz}$ and the shear complex interface stiffness $K^*_{\theta z}$ of the studied interface were obtained. The TTSP was verified for the

<table><caption>**Table 2** 2S2P1D constants for the bitumen at the reference temperature 15 °C</caption><thead><tr><th></th><th>G<sub>00</sub> [MPa]</th><th>G<sub>0</sub> [MPa]</th><th>k [-]</th><th>h [-]</th><th>δ [-]</th><th>τ [s]</th><th>β [-]</th></tr></thead><tbody><tr><td>BP Lavera 160/220</td><td>0</td><td>850</td><td>0.20</td><td>0.53</td><td>2.30</td><td>2.10<sup>-6</sup></td><td>300</td></tr></tbody></table>
<table><caption>**Table 3** 2S2P1D constants for the interface at the reference temperature 15 °C</caption><thead><tr><th></th><th>K<sub>XX_oo</sub> [MPa/mm]</th><th>K<sub>XX_0</sub> [MPa/mm]</th><th>k [-]</th><th>h [-]</th><th>δ [-]</th><th>τ [s]</th><th>β [-]</th></tr></thead><tbody><tr><td>Normal modulus “XX = zz”</td><td>22</td><td>3500</td><td>0.20</td><td>0.53</td><td>2.30</td><td>0.1</td><td>300</td></tr><tr><td>Shear modulus “XX = θz”</td><td>2.2</td><td>420</td><td>0.20</td><td>0.53</td><td>2.30</td><td>0.1</td><td>300</td></tr></tbody></table>

<figure><img src="image_10.png" /><figcaption>**Fig. 10** Master curves of the normalised norm of the normal interface stiffness (all samples), of the shear interface stiffness (all samples) and of the bitumen modulus along with the normalised 2S2P1D model of the bitumen in continuous line</figcaption></figure>

norm of both these moduli and it was quantified with WLF formula. Knowing the 2S2P1D model of the bitumen of the tack coat, the norm of the shear and of the normal interface stiffnesses were also modelled with the 2S2P1D model by using a SHStS transformation. This approach allows modelling the interface viscoelastic behaviour from the viscoelastic behaviour of the bitumen used for the tack coat.

This study shows that the interface between two bituminous mixtures with a tack coat presents a thermo-viscoelastic behaviour for small strain. The norm of the normal and of the shear interface stiffnesses can be modelled using the 2S2P1D model. The interface behaviour is related to the viscoelastic behaviour of the bitumen of the tack coat.

**Acknowledgements** This work is part of the industrial chair between the University of Lyon/ENTPE and the company EIFFAGE dedicated to pavement materials and structures.

**Compliance with ethical standards**

**Conflict of interest** The authors declare that they have no conflict of interest.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

**References**

1. Uzan J, Livneh H, Eshed Y (1978) Investigation of adhesion properties between asphaltic-concrete layers. Annu Meet AAPT 47:495-521

2. Di Benedetto H, Corté JF (2005) *Matériaux routiers bitumeux 2 : constitution et propriétés thermomécaniques des mélanges*.

3. Leutner R (1979) Untersuchung des schichtenverbundes beim bituminosen oberbau. Bitumen 41:84-91

4. Raab C, Partl MN (2004) Effect of tack coats on interlayer shear bond of pavements. In: Proceedings of the 8th Conference on Asphalt Pavements for Southern Africa (CAPSA'04).

5. Collop AC, Sutanto MH, Airey GD, Elliott RC (2009) Shear bond strength between asphalt layers for laboratory prepared samples and field cores. Constr Build Mater 23:2251-2258. https://doi.org/10.1016/j.conbuildmat.2008.11.017

6. Canestrari F, Ferrotti G, Lu X, Millien A, Partl MN, Petit C, Phelipot-Mardelé A, Piber H, Raab C (2013) Mechanical testing of interlayer bonding in asphalt pavements. Adv Interlab Test Eval Bitum Mater. https://doi.org/10.1007/978-94-007-5104-0_6

7. Romanoschi A S (1999) Characterization of Pavement Layer Interfaces. Ph.D. thesis. Louisiana State University.

8. Diakhate M, Phelipot A, Millien A, Petit C (2006) Shear Fatigue Behaviour of Tack Coats in Pavements. Road Mater Pavement Des 7:201-222. https://doi.org/10.1080/14680629.2006.9690033

9. Tozzo C, D'Andrea A, Cozzani D, Meo A (2014) Fatigue Investigation of the Interface Shear Performance in Asphalt Pavement. Mod Appl Sci 8:1-11. https://doi.org/10.5539/ma.v8n2p1

10. Isailović I, Falchetto AC, Wistuba M (2017) Fatigue investigation on asphalt mixture layers' interface. Road
Mater Pavement Des 18:514-534. https://doi.org/10.1080/14680629.2017.1389087

11. Hristov B (2018) Influence of Different Interface Properties on the Interlayer Bond Shear Stiffness. IOP Conf Ser Mater Sci Eng. https://doi.org/10.1088/1757-899X/365/3/032056

12. Canestrari F, Santagata E (2005) Temperature Effects on the Shear Behaviour of Tack Coat Emulsions Used in Flexible Pavements. Int J Pavement Eng 6:39-46

13. Raposeiras AC, Vega-Zamanillo A, Calzada-Pérez MA, Castro-Fresno D (2012) Influence of Surface Macro-texture and Binder Dosage on the Adhesion between Pavement Layers. Constr Build Mater 28:187-192

14. Raab C, Partl M N, Fourquet E, Abd El Halim AO (2017) Static and Cyclic Evaluation of Interlayer Bonding. In: Proceedings of the 10th international conference on the bearing capacity of roads, railways and airfields (BCRRA 2017) 1511-16. https://doi.org/10.1201/9781315100333-216

15. Freire R, Di Benedetto H, Sauzéat C, Pouget S, Lesueur D (2018) Linear Viscoelastic Behaviour of Geocells Interface within Bituminous Mixtures. KSCE J Civ Eng 22(6):2082-2088. https://doi.org/10.1007/s12205-018-1696-9

16. Hristov B (2018) Master Function for Analytical Determination of the Interlayer Bond Shear Stiffness and Fatigue Functions for Asphalt Pavements. MATEC Web of Conf 251:01001. https://doi.org/10.1051/matecconf/201825101001

17. Grellet D, Doré G, Chupin O, Piau J-M (2018) Highlighting of the Viscoelastic Behaviour of Interfaces in Asphalt Pavements – a Possible Origin to Top-down Cracking. Road Mater Pavement Des 19(3):581-590. https://doi.org/10.1080/14680629.2018.1418721

18. Attia T, Di Benedetto H, Sauzeat C, Pouget S (2020) 2T3C HCA, a new hollow cylinder device using Digital Image Correlation to measure properties of interfaces between asphalt layers. Constr Build Mater 247:118499

19. Sutton MA, McNeill SR, Helm JD, Chao YJ (2000) Advances in Two-Dimensional and Three-Dimensional Computer Vision. Photomechanics 77:323-372. https://doi.org/10.1007/3-540-48800-6_10

20. Tkari R, Millien A, Fouchal F, Pop I-O, Petit C (2016) Pavement Interface Damage Behavior in Tension Monotonie Loading. Constr Build Mater 106(March):430-442. https://doi.org/10.1016/j.conbuildmat.2015.12.020

21. Cho SH, Safavizadeh SA, Kim YR (2017) Verification of the Applicability of the Time-temperature Superposition Principle to Interface Shear Stiffness and Strength of GlasGrid-Reinforced Asphalt Mixtures. Road Mater Pavement Des 18(4):766-784. https://doi.org/10.1080/14680629.2016.1189350

22. NF EN 12697-33+A1 (2007) *Mélange Bitumineux - Méthodes D'essai Pour Mélange Hydrocarboné À Chaud - Partie 33: Confection D'éprouvettes Au Compacteur de Plaque.*

23. Nguyen QT, Di Benedetto H, Sauzeat C (2015) Linear and nonlinear viscoelastic behaviour of bituminous mixtures. Mater Struct 48(7):2339-2351. https://doi.org/10.1617/s11527-014-0316-5

24. Mangiafico S, Babadopulos LFAL, Sauzeat C, Di Benedetto H (2018) Nonlinearity of bituminous mixtures. Mech Time-Depend Mater 22(1):29-49. https://doi.org/10.1007/s11043-017-9350-3

25. Airey GD, Rahimzadeh B, Collop AC (2003) Viscoelastic Linearity Limits for Bituminous Materials. Mater Struct 36(10):643-647

26. Olard F, Di Benedetto H (2003) General '2S2P1D' Model and Relation Between the Linear Viscoelastic Behaviours of Bituminous Binders and Mixes. Road Mater Pavement Des 4(2):185-224. https://doi.org/10.1080/14680629.2003.9689946

27. Di Benedetto H, Olard F, Sauzéat C, Delaporte B (2004) Linear Viscoelastic Behaviour of Bituminous Materials: From Binders to Mixes. Road Mater Pavement Des 5(sup1):163-202. https://doi.org/10.1080/14680629.2004.9689992

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.