Currently, there is a lack of in-situ or model test results for cone penetration tests conducted in deep dense sand layers under high overburden stresses, restricting the development of empirical relationships between cone penetration test results and the characteristics of such deep dense sand layers. This study addresses this gap by proposing an empirical relationship to predict the relative density of dense silica sand based on stress level and cone tip resistance. The relationship is developed through cone penetration tests performed in a calibration chamber using dense sand specimens subjected to high stresses, along with numerical simulations employing the large deformation finite element method. The Arbitrary Lagrangian Eulerian method was employed to regularly regenerate the mesh to prevent soil element distortion around the cone tip. Additionally, the modified Mohr-Coulomb model was integrated to capture the stress-strain behavior of dense silica sand under high stresses. The empirical relationship provides predictions with acceptable accuracy, as the discrepancies between the predicted and measured relative density values fall within ±30%.