How do you normally go about initializing the stresses in the domain? I know it’s possible to input initial particle stresses using a text file, but how do you go about doing this if you’re specifying the particles to be created at the Gauss points?
Seems like you’d calculate the initial stresses by calculating the vertical stress and assuming some K_0, finding the stresses acting at each particle level. If I don’t specify particle locations and just let the software generate them at the Gauss points, is there a way to provide these initial stresses in a way that they can be added at the correct locations?
I was considering using time-varying forces to first get the particles in equilibrium based on gravity and then impose some additional loading. However, doing it this way, there are massive settlements initially such that the material points aren’t anywhere close to the Gauss points so the stresses calculated are not reliable. Right now, it seems like I’d have to manually determine where the Gauss points would be for the mesh and create an input file that contains this information. Then, calculate the stresses that match the ordering of these point and put that in the input file for the initial stresses.
When the initial stresses and strains are important to capture correctly, I’m guessing this is the way that folks have been dealing with it. Just want to make sure I’m not making things harder than they need to be or that I’m not reinventing the wheel here…