yade-users team mailing list archive

[Karol Brzezinski <question700199@xxxxxxxxxxxxxxxxxxxxx>]

Question #700199 on Yade changed:
https://answers.launchpad.net/yade/+question/700199

    Status: Open => Answered

Karol Brzezinski proposed the following answer:
Hi,

you can iterate over interactions and find read penetration depth. Based
on this, you can compute overlapping volume manually. Please note it is
not a ready-for-use code, but just schematic pseudocode:

###
totalOverlapVolume = 0

for ii in O.interactions: #iterate over all the interactions
    b1 = O.bodies[ii.id1]
    b2 = O.bodies[ii.id2]
    penetrationDepth = ii.geom.penetrationDepth # Should work for Mindlin law. You can inspect interactions manually to find out whether .penetrationDepth is a right property to read.
    overlapVolume = yourFunctionComputingOverlap(b1,b2,penetrationDepth)
    totalOverlapVolume += overlapVolume

###
yourFunctionComputingOverlap() should check the shape of the particles b1 and b2 (whether it is a sphere, wall, wall, etc), extract some shape properties (e.g. radii) and compute volume based on geometrical relationships.

Best wishes,
Karol

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.