@itascan
If I have several layers of soft soil with two constitutive models mohr-coulomb and modified cam clay, arranged in an irregular pattern (cam clay between mohr coulombs and vice versa), is it okay to initialize initial stress by simple zone initialize stress command and assigning gravity for any layers irrespective of their positions? I have seen examples where the pressure-effective for cam clay is calculated separately, so how to initialize stress in soil layers arranged in irregular pattern of constitutive model?
1 Like
You can use zone initialize stress command to initialize stress, however, for the modified cam clay model, the mean effective pressure must be initialized before stepping, this can be done via a fish function:
fish define camclay_ini_p
loop foreach local p_z zone.list
local mean_p = -(zone.stress.xx(p_z) + zone.stress.yy(p_z) …
+ zone.stress.zz(p_z))/3.0 - z.pp(p_z)
zone.prop(p_z,‘pressure-effective’) = mean_p
endloop
end
1 Like