yade-users team mailing list archive

[William Chevremont <question676053@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Open => Answered

William Chevremont proposed the following answer:
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).

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