Natural clays exhibit inherent anisotropy, a result of their stress history and the size and arrangement of particles during sedimentation. As a consequence, when identical soils undergo the same loading program, their strengths differ depending on the testing angle used. Cone penetration testing (CPTu) stands as one of the most established in-situ tests. During this test, the soil undergoes significant deformations and experiences a rotation of principal stresses. While the mobilization of undrained shear strength is uncertain for isotropic materials, with ambiguity about whether which undrained shear strength is mobilized (triaxial compression, triaxial extension…), the uncertainty is even greater for anisotropic materials. As numerical methods advance, more realistic CPTu simulations become feasible, employing a combination of hydromechanical formulations and appropriate constitutive models. These simulations provide valuable insights into the stress paths induced to the soil during CPTu testing. In this study, we use the Particle Finite Element [1] method to simulate CPTu in undrained anisotropic clays. The constitutive response is represented by SCLAY-1 [2], a modified version of Modified Cam-Clay that integrates mixed isotropic and rotational hardening with an associated flow rule. The model assumes isotropic elasticity inside the yield surface, which is a rotated ellipsoid in the p’- q plane (triaxial conditions), and the model is equipped with two hardening laws: one describes how the size of the yield curve changes with plastic volumetric strains, while the other relates to the creation or removal of fabric due to plastic straining. This work reports two sets of simulations. The first set aims to explore which type of shearing is mobilized during CPTu testing in isotropic materials; this is achieved by simulating CPTu testing in soils with identical undrained triaxial compression strength but different triaxial extension strength. The second set of simulations illustrates the effect of initial fabric and rate of destructuration on cone metrics and cone factors.