I am analyzing different lateral earth pressure conditions (K₀), and I need to obtain an accurate initial stress distribution for values up to K₀ = 1.5. My current approach and the challenges are summarized below:

• Initially, I relied on the stress state generated by gravity loading using a Mohr–Coulomb plasticity model. With a Poisson’s ratio of 0.33, this produced a reasonable stress distribution corresponding to K₀ ≈ 0.5.
• To achieve higher K₀ values (e.g., 1.5), it is not physically valid to increase Poisson’s ratio above 0.5. Therefore, I introduced a predefined initial stress field to impose the desired lateral stress.
• Because Abaqus/Explicit (CEL) is sensitive to dynamic effects, I applied non-reflecting boundary conditions and used a smoothed amplitude curve to minimize wave propagation during gravity loading.

• Despite these measures, the stress distribution still does not match the stable, near-equilibrium field that can be achieved in Abaqus/Standard using the geostatic step. The Explicit model continues to show significant transient variations, and the resulting K₀ distribution is not sufficiently accurate.

  I would greatly appreciate guidance on any of the following:

• Best practices for generating and transferring a geostatic or K₀-consistent stress state into an Abaqus/Explicit CEL analysis—for example, through Abaqus/Standard → Explicit workflows, restarts, or keyword-based stress initialization.
• Recommended boundary-condition setups in CEL that maintain lateral confinement and reduce wave reflection while establishing the initial stress field.
• Any example input files, documentation, or published references demonstrating successful K₀ initialization for slope problems in CEL.
• Any alternative methods you would recommend for achieving a realistic initial K₀ stress state in Abaqus/Explicit (CEL).