In this work, the complex coupled deformation and flow processes occurring in the soil around a CPTu penetrometer during a test in structured natural clays are investigated by means of the Geotechnical Particle Finite Element Method (GPFEM) code. The GPFEM implementation adopted incorporates a fully coupled hydro-mechanical formulation, based on regularized, mixed low-order linear strain triangles. To capture the relevant features of the mechanical response of natural structured clays, the soil behaviour is described using a finite deformation, non-associative elastic-plastic model for this class of geomaterials, referred as FD_Milan model. The model formulation is based on a multiplicative decomposition of the deformation gradient and on t...