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.
_______________________________________________
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
