New question #647591 on Yade:
https://answers.launchpad.net/yade/+question/647591
Here is a part of my codes:
__________________________________________________________________________________________________________________
contact = Law2_ScGeom_FrictViscoPhys_CundallStrackVisco(traceEnergy=True)
newton=NewtonIntegrator(damping=damp)
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Box_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(),Ig2_Box_Sphere_ScGeom()],
[Ip2_FrictViscoMat_FrictViscoMat_FrictViscoPhys()],
[contact]
),
## We will use the global stiffness of each body to determine an
optimal timestep (see
https://yade-dem.org/w/images/1/1b/Chareyre&Villard2005_licensed.pdf)
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
triax,
newton,
# save data for Paraview
VTKRecorder(fileName='post/3d-vtk-',recorders=['all'],iterPeriod =
50000)
]
O.dt=.5*utils.PWaveTimeStep()
__________________________________________________________________________________________________________________
During my calculation, I print the O.dt and utils.PWaveTimeStep and find that
O.dt is larger than utils.PWaveTimeStep.
I really wonder whether the O.dt is determined by the codes
"O.dt=.5*utils.PWaveTimeStep()" or timestepSafetyCoefficient?
The particles I ues are 2-sphere clumps.
Thank you very much.
--
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 : [email protected]
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
