Hi Kate. I'm a little confused here, and something seems circular:
"Cr = (e1-e2)/log(sigma'p/sigma'1)
In above equation I have everything except sigma'p, but I found out that I can calculate it if I have OCR from equation:
OCR = sigma p/gamma*z
Do you have a consolidation test with rebound? If so, you have Cr and sigma'p as directly as possible.
By definition, OCR is sigma'p/[sigma' at present]. The denominator is (assuming a single layer) gamma*depth - u. The OCR equation you have seems to be neglecting the pore pressure u. If there is no pore pressure, then the present value of sigma' is gamma*z (assuming z is the depth), or total stress = effective stress.
Have never used Plaxis, so I don't know what its soft soil model is like. Does Plaxis understand pore water pressure? Our people use FLAC, and we go for the simplest models we can justify for any particular case, often just elasto-plastic Mohr-Coulomb, but our MC strength is always based on effective stress. (FLAC does understand water tables, gridded pore pressure, etc.)
Also, what is MES?