Box/Wall - Sphere contact with periodic boundaries
Hi,
I'm lost about the contact between box or wall with spheres with periodic boundaries (this is not related to the allowBiggerThan
Consider 9 spheres lying on a box, in a periodic cell. The cell includes exactly the 9 spheres, and the box may change in size. Its surface may approach the one of the cell, or may be bigger (see script below)
- if the surface of the box is 90% the one of the cell (set howBig = 0.45 in the script below), I get my 9 sphere-box interactions, fine !
- if the box fills exactly (in surface) the cell (howBig = 0.5) => only 6 spheres are detected as in contact with the box.. (or 4 if I put radius = 0.5, so that there is normally exact geometrical contact)
- if howBig = 0.6 (the box extends beyond the limits of the cell) => 6 spheres in contact, for both radii tested (0.5 - exact contact - or 0.55 - overlap)
- if howBig = 1000 => there are 3 spheres in contact with the box.....
If I use a wall I got 3 spheres-wall contact for radius = 0.5, and 6 for radius = 0.55, but not 9.... May someone explain to me these results ? (or confirm that there is a bug ?)
Jérôme
PS : the script :
# encoding: utf-8
O.periodic=1
O.cell.
for i in range(0,3):
for j in range (0,3):
howBig=0.45
O.bodies.
O.engines=[
ForceResetter()
,InsertionSort
,InteractionLoop(
[Ig2_
[Ip2_
[Law2_
]
O.dt=1.e-5
O.step()
print '\nhowBig =',howBig,':', O.interactions.
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Jérôme Duriez for more information if necessary.