yade-users team mailing list archive

[xuelong14 <question647591@xxxxxxxxxxxxxxxxxxxxx>]

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.