A virtual calibration chamber was built using a three-dimensional model based on the discrete-element method. The chamber was then filled with a scaled granular equivalent of Ticino sand, the material properties of which were selected by curve-fitting triaxial tests. Cone penetration tests were then performed under different initial densities and isotropic stresses. Penetration resistance in the virtual calibration chamber was affected by the same cone/chamber size effect that affects physical calibration chambers and was corrected accordingly. The corrected cone resistance obtained from the virtual calibration chamber cone penetration tests shows good quantitative agreement with correlations that summarise previous physical results.

A three-dimensional discrete element model is used to investigate the effect of grain crushing on the tip resistance measured by cone penetration tests (CPT) in calibration chambers. To do that a discrete analogue of pumice sand, a very crushable microporous granular material, is created. The particles of the discrete model are endowed with size-dependent internal porosity and crushing resistance. A simpli-\ud fied Hertz–Mindlin elasto-frictional model is used for contact interaction. The model has 6 material parameters that are calibrated using one oedometer test and analogies with similar geomaterials. The calibration is validated reproducing other element tests. To fill a calibration chamber capable of containing a realistic sized CPT the discrete analogue is up-scaled by a factor of 25. CPT is then performed at two different densities and three different confinement pressures. Cone tip resistance in the crushable material is practically insensitive to initial density, as had been observed in previous physical experiments. The same CPT series is repeated but now particle crushing is disabled. The ratios of cone tip resistance between the two types of simulation are in good agreement with previous experimental comparisons of hard and crushable soils. Microscale exploration of the models indicates that crushing disrupts the buttressing effect of chamber walls on the cone. DEM modelling of cone penetration tests in a double-porosity crushable granular material.Peer ReviewedPostprint (author's final draft

A virtual calibration chamber was developed using a three-dimensional (3D) discrete element method (DEM) to perform cone penetration tests (CPTs) on a discrete analogue of Ticino sand. The macroscale response of the DEM model was previously shown to be in good quantitative agreement with that of analogous physical models. In the current study the performance of the model at meso and microscale levels of resolution is examined. The microscale response is examined using particle displacements and contact force distributions. The mesoscale behaviour is examined using stress and strain fields obtained through appropriate averaging and interpolating procedures. Four CPTs are examined at the steady-state penetration stage. The effects of radial boundary conditions, initial stress state, initial average density, and particle rotational inertia are examined. The ability of the micro and mesoscale data to identify and explain the relevant mechanisms underlying the significant differences in the macroscale response of the models is discussed. Comparisons with similar phenomena observed in physical tests are also highlighted.

A three—dimensional numerical model was implemented to PFC3D Code to simulate cone penetration test in Ticino sand in a calibration chamber. The model is calibrated using laboratory test results. The full/half/quarter calibration chamber was used to examine the effect of symmetry on the results. Examination of specimen homogeneity was done by (i) visual observation of the network of the contact forces developing between particles, (ii) examination of the porosity, d50 and Cu distributions inside the specimen by using a representative elementary volume (REV). Some overall results from these simulations are also presented here and compared with the experimental results from a calibration chamber test database.

In this paper the development and experimental validation of a numerical model of two-dimensional unsaturated flow in a double-porosity medium is presented. The model is based on the coupled formulation for flow in macro- and micropores obtained by homogenization. It was applied to simulate the axisymmetrical tension disk infiltration experiments that were carried out in a double-porosity medium. The physical model was a three-dimensional periodic structure, composed of porous spheres made of sintered clay and embedded in Hostun fine sand HN38. The hydraulic parameters of both porous materials were determined by inverse analysis of independent infiltration experiments performed on sand and sintered clay. The effective parameters of the double-porosity medium were calculated from the solution of the local boundary value problem, obtained from the homogenization procedure. The cumulative infiltration curve and the global dimensions of the humidified zone obtained from the numerical solution are in good agreement with the observations. Moreover, numerical simulations showed the existence of a narrow zone of local nonequilibrium that moves with the infiltration front. Upstream of this zone, the infiltration bulb is in the local equilibrium conditions.

Calibration chambers are frequently used to verify, adapt, or both verify and adapt empirical relations between different state variables and in situ test results. Virtual calibration chambers (VCC) built with 3D discrete element models may be used to extend and partially substitute costly physical testing series. VCC are used here to explore the mechanics of flat dilatometer penetration and expansion. Results obtained for a simulation of physical tests in Ticino sand are presented. Blade tip resistance during penetration is in good agreement with the experiments. A piston-like design is used for the blade so that larger displacements may be applied than it is possible with a membrane. Initial piston pressures in the expansion curves are very low, strongly affected by the scaled-up grain sizes. Despite that difficulty, expansion curves may be easily interpreted to recover dilatometer moduli ED close to those observed in the physical experiments. Particle-scale examination of the results allows a firmer understanding of the current limitations and future potential of the technique.