yade-users team mailing list archive

[Henry <question253483@xxxxxxxxxxxxxxxxxxxxx>]

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

Henry posted a new comment:
I modified the  "triax.py" file in "example/concrete/", according to the
example file of densityscale (timeStepperUsage.py).

The followin script is the modified part:

.....

O.engines=[
	ForceResetter(),
	InsertionSortCollider([
		Bo1_Sphere_Aabb(aabbEnlargeFactor=enlargeFactor,label='bo1s'),
		Bo1_Facet_Aabb()
	]),
	InteractionLoop(
		[
			Ig2_Sphere_Sphere_ScGeom(interactionDetectionFactor=enlargeFactor,label='ss2d3dg'),
			Ig2_Facet_Sphere_ScGeom(),
		],
		[
			Ip2_CpmMat_CpmMat_CpmPhys(cohesiveThresholdIter=O.iter+5),
			Ip2_FrictMat_CpmMat_FrictPhys(),
			Ip2_FrictMat_FrictMat_FrictPhys(),
		],
		[
			Law2_ScGeom_CpmPhys_Cpm(),
			Law2_ScGeom_FrictPhys_CundallStrack(),
		],
	),
	PyRunner(iterPeriod=1,command="addForces()"),
	NewtonIntegrator(damping=.3,label="newton"),
	CpmStateUpdater(iterPeriod=50,label='cpmStateUpdater'),
	PyRunner(command='plotAddData()',iterPeriod=10),
	PyRunner(iterPeriod=50,command='stopIfDamaged()'),
      GlobalStiffnessTimeStepper(label='ts'),
]

# run one step
O.step()

# reset interaction detection enlargement
bo1s.aabbEnlargeFactor=ss2d3dg.interactionDetectionFactor=1.0

O.timingEnabled=1
ts.active=True
newton.densityScaling=True


#... but not that of cell, hence we have to adjust it or the problem becomes unstable
for b in O.bodies:
    if isinstance(b.shape,Sphere):
        b.state.mass=1.0e20*b.state.mass
        


# initialize auto-updated plot
if runInGui:
  .........


Thanks a lot!

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.