4.2.1. Laboratory Drop Weight Impact
The multi-layer pavement slab was subjected to an 1181-kg drop weight impact. The cross-section of the multi-layer pavement is the same as that given in
Figure 7. The basic mechanical properties of the, HSC, ECC and AC are determined according to ASTM standards, and the results are summarized in
Table 9. The drop weight was a cylindrical projectile with a hemispheric head (100 mm in diameter and made with high strength steel), and the pavement slab was subjected to two times of impact from the same drop height of 1.5 m. During the test, the multi-layer pavement slab was placed on the top of compacted soil/sand in a steel strong box. Directly below the slab was the geocell (MiraCell MC-100) which was filled with compacted soil/sand. This was to enhance the strength of the soil/sand layer and to provide a solid sub-base as to simulate the practical condition. The setup for the multi-layer pavement slab is given in
Figure 13. Various instruments were also installed to monitor the response of the pavement during the drop weight test.
Figure 14 shows the positioning of the potentiometers. A photodiode system was used to trigger the data acquisition system during the test. It consists of two photodiodes and two laser sources placed 100 mm vertically apart. The data acquisition system would be triggered when the falling projectile crosses the top laser emitter. Impact velocity could be determined using the time interval that the projectile took to cross the second laser emitter.
In the numerical model, the AUTOMATIC_SURFACE_TO_SURFACE contact algorithm is employed to model the interaction between pavement slab and soil. The contact algorithm TIEBREAK is also used to simulate the interface behaviour between the AC and HSC layer, which is the same as that in the numerical model for blast loading as mentioned in
Section 4.1.2. As the impactor might penetrate the AC layer, the erosion method is employed in the simulation, in which the maximum plastic failure strain of 0.2 is used to delete the distorted element once the actual strain exceeds this predefined failure strain. The multi-layer pavement slab, the drop weight head and soil mass are discretised in space with one point Gauss integration eight-node hexahedron elements. The geogrid is simulated with the four-node Belytschko–Tsay shell element that allows no bending resistance. In addition, the PLASTIC_KINEMATIC material model is employed to describe the bi-linear behaviour of the geogrid under tensile loading [
23]. The parameters used for the geogrid are the same as given in
Table 4.
The Drucker–Prager model is used to simulate soil mass. In the laboratory drop weight test, the upper soil layers were compacted and reinforced with geocell material, which would enhance the strength of the soil, and the lower layer has no reinforcement. Hence, in the numerical model, it is necessary to consider the function of the geocell material. From the laboratory test on the geocell-encased sand [
26], it is observed that the geocell confinement did not change the friction angle of soil while significant cohesion occurred in the granular soil, which indicated that for the geocell-reinforced sand layer, the strength and stiffness behaviour of the soil would be enhanced. However, in the numerical model, it is difficult to model and mesh the geocell material due to its complex geometry. Hence, it would be preferable to use the composite model to consider the enhancement of the shear strength and stiffness of the geocell-reinforced sand layer. Madhavi et al. [
27] purposed an empirical equation to calculate Young’s modulus of the geocell-reinforced sand using the secant tensile modulus of the geocell material and Young’s modulus parameter of the unreinforced sand, which could be expressed as:
in which
is the Young’s modulus of the geocell-reinforced sand,
M is the secant modulus of the geocell material at axial strain 2.5% in kN/m and
is the confining pressure from the geocell in kPa.
is the dimensionless modulus parameter of the unreinforced sand, which is a modulus number in the hyperbolic model developed by Duncan and Chang [
28]. The confining pressure
could be calculated as:
where
is the initial diameter of the geocell and
is the axial strain of the geocell at failure; the induced cohesion in the geocell-reinforced sand is then related to the increase in the confining pressure
:
in which
is the enhanced cohesion and
is the coefficient of passive earth pressure.
In the current study, the geocell MC-100 from Polyfelt is used. The geocell dimension is a rhombus with two diagonal lengths of 203 mm and 244 mm. Thus, the equivalent diameter is calculated as about 177.5 mm. The secant modulus
M for the geocell is obtained as 278 kN/m from the tensile test [
20]. Additionally, the value of
is taken as 4.8%. The modulus parameter
for unreinforced sand in the current study is taken as 727 MPa according to the curve fitting from the triaxial test, and hence, the confining pressure, enhanced cohesion and Young’s modulus for geocell-reinforced sand could be calculated based on Equations (21)–(23). The parameters for the unreinforced sand and geocell reinforced sand are summarized in
Table 10.
During the impact test, the deformation of the drop weight head was negligible compared to the deformation of the pavement slab. Hence, the drop weight head is modelled with a rigid body in the current study. For the configuration of the drop weight head, the simple cylindrical shape is modelled instead of modelling the head with weight mass (as shown in
Figure 15). The simple cylindrical head has a diameter of 100 mm with a length of 1292 mm. The total mass for the simple cylindrical head is about 1181 kg, from which the density of the drop head would be obtained. The properties for the drop head are listed in
Table 11. The convergence study is conducted, and it was found that a 5-mm element size gave a stable response, which is therefore applied for the simulation. The numerical model of the multi-layer pavement slab under drop weight impact load is given in
Figure 16.
For the 1.5-m drop weight impact, the drop weight head is assigned with 5.02 m/s for the first impact. For the second impact, due to the penetration of the AC layer in the first impact, the distance between the laser diode system and the face of the pavement slab is increased, while in the numerical model, the impact head is just placed at the position right before reaching the surface. Hence, the initial velocity in numerical model for the second impact is determined as the sum of the experimental recorded velocity and velocity caused by gravity acceleration. Thus, the velocity is calculated as 5.06 m/s for the second impact with the gravity acceleration of 9.8 m/s2.
After the first impact, the fully restarted method in LSDYNA is used to conduct the second impact simulation. At the beginning of the second impact, the stress and residual velocity within the pavement slab and the velocity from the first impact are set to be zero. The damage factor and plastic strain is retained in order to check the accumulated damaged behaviour after the second impact. The simulation of the second impact is carried out when the downward velocity of impactor reached zero, as it is very time consuming to continue to simulate the vibration of pavement slab after impact. It should be noticed that the pavement slab stopped rebounding, which might bring numerical errors (i.e., energy unbalance).
4.2.2. Numerical Result
The damaged situations for the multi-layer slab under the first drop weight impact for the experiment and numerical model are presented in
Figure 17. It is found from the numerical simulation that the AC layer is penetrated through, and the drop weight head is impeded by the HSC layer due to the high compressive strength. Such findings from the numerical model are consistent with the observation in the physical test. It is observed that no severe and moderate damage happened in the ECC layer. Additionally, the integrity of both HSC and ECC layers is kept after the first impact.
The damage situations for the multi-layer pavement slab under the second drop weight impact and the numerical model are given in
Figure 18. The simulation results indicate that the HSC layer stopped the impactor; however, severe damage happens at the top surface and in the middle at the side of the HSC layer. Severe damage occurs at the rear surface of the HSC layer, as well. For the ECC layer, it is shown that the cracking occurs and propagates from the centre at the top surface. However, the damage area is smaller than that of the HSC layer. The cracks are also found in the middle at the side of the ECC layer, and similar to the experiment, the crack does not propagate through the thickness of the layer.
Figure 19 compares the numerical results with the experimental results about the rear surface of ECC layer after the second impact for the multi-layer slab. It is found that the damaged contour and cracking pattern were similar. Such results further indicate the reasonable close agreement between the numerical analysis and the experimental results using the KCC model for AC material.
Table 12 compares the numerical results with the experimental data for the multi-layer slab under two times of the drop weight impact. After the first impact, the settlement from the numerical results is close to that obtained from the experiment at P2. A big deviation is found in between the numerical results and experiment results on the settlement at P3. This might be attributed to the dislocation of the potentiometer at P3 upon impact, which could lead to a misreporting of the settlement at P3. In addition, the settlement at P1 is lower in the numerical results compared with the experimental results. This could be attributed to the erosion process in the numerical model or, more specifically, the energy released phenomena. At the second impact, since no erosion technology is employed for the HSC layer, the numerical results are more reasonable. In the impact simulation, the HSC layer suffered deformation instead of being penetrated, and hence, no erosion technique was adopted to trigger the energy released phenomena. Thus, the settlement profile of the pavement slab decreased with the increase of the radial distance from the centre, which indicated that the whole pavement slab bent during the impact. In general, the numerical models provide reasonable estimations on the settlements of the multi-layer slab under impacts, in particular under the second impact load, though discrepancies exist at certain regions.