Export stress contour of interaction between polyhedra and spheres

Hi,
Very recently, I introduce PDFEngine, that seems to make what you want, but I tested it only for Sphere-Sphere interaction (maybe need extra development to handle polyhedra-sphere interaction). It give you global stress repartition, distance between objects and relative velocity.

To use it, you just have to add

PDFEngine(filename="PDF.txt", virtPeriod=.1)

to your engines. It'll write "PDF.txt" with column like normalStress_xx(theta, phi). You can change angle discretization (see options [1])

Be sure to update yade before using it (committed only 3 days ago)

[1] https://yade-dem.org/doc/yade.wrapper.html?highlight=pdfengine#yade.wrapper.PDFEngine

Thank you so much dear William :)

your advices helped me a lot but, is it possible for you to send me your full script about this problem that you exported stress contour about Sphere-Sphere interaction which you mentioned above .

Hi,

Full script is useless because of very specific stuff done inside, this will confuse more than explain.

The only modification to be done is in the O.engines part. My looks like:

engines = [ForceResetter(),
    InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor=1.5, label="aabb")],verletDist=-0.1,allowBiggerThanPeriod=False),
    InteractionLoop(
        [Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.5,label="Ig2")],
        [Ip2_FrictMat_FrictMat_LubricationPhys(eta=m_viscosity,eps=m_epsilon)],
        [Law2_ScGeom_ImplicitLubricationPhys(activateNormalLubrication=True,
                                                                activateTangencialLubrication=True,
                                                                activateTwistLubrication=True,
                                                                activateRollLubrication=True,
                                                                debug=True,
                                                                resolution = 3,
                                                                theta = 1,
                                                                SolutionTol = 1.e-8,
                                                                MaxIter = 50)]
    ),
    NewtonIntegrator(damping=0.),
    GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8, defaultDt=1e-6,label="TimeStepper",viscEl=False),
    PDFEngine(filename="PDF.txt", virtPeriod=.1)
    ];

The output file (PDF.txt) is then post-proceeded with MatLab in order to plot the stress tensor component magnitude as a function of the position (theta, phi).