simSigmaP function

Simulating Preconsolidation Stress