Yade

Pfacet in periodic boundary conditions

Asked by Alice Ezzeddine

Hello!
I am trying to simulate a plate being installed in sand. For that, I imported a 3d mesh of the plate with "gtsPFacet". I then set walls for the sample in x and y directions while periodic boundary conditions on z. I saw on the forum that it can be done by setting O.cell.setBox on x and y directions to large values.
The code for yade and the mesh file are in the following link:
https://filesender.renater.fr/?s=download&token=6328dd93-0240-4c43-bd2e-6cc02f35fe44

My questions are:

1. If I would like to have a plate that is not hollow inside, should I provide the mesh of the interior of the plate? I mean not just
just the external geometry or am I wrong?

2. Is it okay how I simulated periodic boundary conditions on the z dimension or am I missing something?

If you have any other suggestions or comments please let me know.
Thank you a lot!
Alice

Question information

Language:
English Edit question
Status:
Solved
For:
Yade Edit question
Assignee:
No assignee Edit question
Solved by:
Bruno Chareyre
Solved:
Last query:
Last reply:
Revision history for this message
Bruno Chareyre (bruno-chareyre) said : 		#1

Hello, I don't understand the problem very well.
What is x,y,z? Is there gravity (which axis?)?
Why do you care about the plate being hollow? Isn't it a rigid object described by it's external shape?
Regards
Bruno

Revision history for this message
Alice Ezzeddine (aliceezzeddine) said : 		#2

Hi,
x, y, and z are the 3 directions in a coordinate system.
My point is I would like to set walls as boundaries on two directions ( by setting wallMask = 15) while boundary conditions on the third direction. Would it work by setting O.cell.setBox on the directions with walls to large values?

There is an acceleration of 100*gravity on the y-axis as I am trying to simulate an experiment in the geotechnical centrifuge.

Whether it is hollow or not plays a role in the mass applied on the sample.

Thank you,
Alice

Revision history for this message
Best Bruno Chareyre (bruno-chareyre) said : 		#3

For the uniaxial periodicity, I think yes, what you do is correct.

The hollowness is really not a problem. If you know the mass you can just assign it: plate.phys.mass = whatever.

I hope it helps.

Bruno

Revision history for this message
Alice Ezzeddine (aliceezzeddine) said : 		#4

It does help, thank you!

Alice

Revision history for this message
Alice Ezzeddine (aliceezzeddine) said : 		#5

Thanks Bruno Chareyre, that solved my question.