October 5, 2018 at 4:49 amRetroBJJSubscriber
I just want to ask about my problem, which is related to excavations, open pit mine stability, etc. Geomechanics in general. There's something called in situ stress, and below ground it exists. These are applied at opposite faces of our model, and have the same magnitude at opposite faces of the outer boundary of the model. (See pic attached)
So basically, I need to apply a compressive pressure in a face x with magnitude X, and in the opposite face I need to apply THE SAME pressure. Same procedure needs to be followed for the rest of the faces of the box.
The issue is that, they seem to cancel to each other, so the final model is like it was never under the action of in situ stress, and only under the effect of gravity... Is there a way to fix this? it's really frustrating, because I am not getting the results I am supposed to get.
Cheers and thanks in advance. Hope you get my problem.
October 5, 2018 at 4:30 pmpeteroznewmanSubscriber
The good setup is to take a cube, and on three faces, apply a displacement constraint with zero displacement in X on the face normal to the X axis, a zero displacement in Y on the face normal to the Y axis, and a zero displacement in Z on the face normal to the Z axis. Now you have a cube that isn't going anywhere.
On the opposite face to the one that had the displacement constraint, apply a pressure. The reaction force on the face with the displacement constraint will have an equal and opposite force to the face with the pressure applied.
Let's say Z axis points up and gravity is working in the -Z direction. If you are doing it right, the pressure on the side walls should be a function of Z.
October 6, 2018 at 2:33 amRetroBJJSubscriber
Peter, thanks so much for your answer.
When I first read your comment, it made me a lot of sense, in terms of physics. But as you may know, in other softwares this situation is a bit different (in ABAQUS for example, stresses don't cancel each other out when applied at faces).
I recently have made other simulations, but the results I'm getting aren't the expected ones.
It's supposed to have increased displacements in the inner part of the excavation, but the image attached doesn't show that sadly. I know that in this case the model is symmetric and I could have the desired results by cutting it in half, but that's not the idea... I want to know the general procedure to model this kind of situation.
I've tried with fixed boundaries, zero disp BC, etc. Any other advice dear Peter?
Have a nice weekend.
ps: Another question... is there a way to change the angle from radian to degrees in the engineering data? I have friction angle in radians and that's weird to me.
October 8, 2018 at 6:40 pmRetroBJJSubscriber
Is there a different way to proceed in this case? please, this problem is truly taking my sleep away.
October 8, 2018 at 10:08 pmpeteroznewmanSubscriber
I don't know what you were expecting with the attached image, or what you meant by "It's supposed to have increased displacements in the inner part of the excavation", maybe you could explain it in more detail.
If you want to upload a Workbench Project Archive, I will look more closely at your model.
October 8, 2018 at 10:52 pmRetroBJJSubscriber
Peter, thanks again for your time and efforts in helping me.
When I said "increased displacements in the inner part of the excavation" I meant that the right simulation would have given me a displacement contour where vectors have direction towards the free faces of the excavation, so it's like the excavation shrunk (the technical term is convergence). The rest of faces have to remain without displacement (outer part of the box). In consequence, the highest displacements occur near the excavation.
Also, the min and max principal stresses should always be concentrated around the excavation (induced stresses), like in this example (keep in mind that this is not the same project, just an example to illustrate the idea):
And well... that's it. As a remainder, this is a cube (box) with an excavation (mining stope) inside and is under the action of s1, s2 and s3 compressive stresses just like this image below. And, the boundaries of the box should be constrained to have zero displacements.
The issue is that these stress cancel each other out if I model this situation as in other softwares. I hope now you can understand what I meant.
My project is attached, dear Peter. Thanks again for everything.
October 9, 2018 at 1:26 ampeteroznewmanSubscriber
Here is what I wanted you to do. For the Y face, set the displacement in Y to 0 while leaving X and Z free.
What you did was set all three to zero like is shown below.
If you make the same change for dx and dz, setting just one component to zero while leaving the other two free, you will get the following result.
If you plot the Y displacement, you can see the Average is -93 mm.
So if instead of Y = 0 you set Y = 93, then the center of the cave will have zero displacement and the cave walls will collapse in. Do the same for the X and Z axis.
I hope you like this.
Attached is a zip file with an AIM 19.2 project
October 9, 2018 at 2:13 amRetroBJJSubscriber
Hello again Peter, that looks a lot better! thanks for your kind support.
But... the displacement doesn't seem right to me. Is there a way to change this?
Ideally, higher displacements should be near the excavation, while the outer part of the box should be zero. I understand that what we got is an obvious result of the boundary conditions applied...so, by following this procedure should I renounce at looking at the deformation? (1)
A second question Peter, to finally close this topic... in the 'engineering data' in the material library, I have friction angle in radians and I can't convert it to degrees. Is there a way to change it? because I have that cell in gray, and in the unit systems I don't see the option :S (2)
Also, I don't understand why this value is also so low... 0.01 rad is like 0.57 deg, which is non realistic for the materials that I am aware of.
Thanks again, kind regards.
October 9, 2018 at 3:08 amRetroBJJSubscriber
October 9, 2018 at 3:28 ampeteroznewmanSubscriber
1) Plotting strain is the best way to see the local deformation in the material. Don't focus on the displacement because that just increases linearly from the point that didn't move, so the magnitude is not that useful.
2) You are stuck with radians for that field. Other places in ANSYS you can use degrees. I don't know much about Mohr-Coulomb material model so I can't comment on the magnitude of that parameter.
Good Luck and Cheers,
October 9, 2018 at 3:30 amRetroBJJSubscriber
Peter, thanks so much for all your time and kind support.
Have a fabulous week!
- You must be logged in to reply to this topic.
Boost Ansys Fluent Simulations with AWS
Computational Fluid Dynamics (CFD) helps engineers design products in which the flow of fluid components is a significant challenge. These different use cases often require large complex models to solve on a traditional workstation. Click here to join this event to learn how to leverage Ansys Fluids on the cloud, thanks to Ansys Gateway powered by AWS.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
- Saving & sharing of Working project files in .wbpz format
- Understanding Force Convergence Solution Output
- An Unknown error occurred during solution. Check the Solver Output…..
- Solver Pivot Warning in Beam Element Model
- Colors and Mesh Display
- whether have the difference between using contact and target bodies
- What is the difference between bonded contact region and fixed joint
- How to calculate the residual stress on a coating by Vickers indentation?
- The solver engine was unable to converge on a solution for the nonlinear problem as constrained.
- User manual
© 2023 Copyright ANSYS, Inc. All rights reserved.