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.
_______________________________________________
Mailing list: https://launchpad.net/~yade-users
Post to : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
