← Back to team overview

yade-users team mailing list archive

Re: [Question #204014]: Not sure about utils.PWaveTimeStep()

 

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

Christian Sommerfeld posted a new comment:
Hi,

thank you guys!

@Jan: I found it in the DEM-Formulation at the point "Estimation of
\Dtcr by wave propagation speed"

My programm is still working but i'm not sure if i get the correct
results.

So now i understand that i have to use both (PWaveTimeStep and PWaveTimeStep) in my case. Because my facet is dynamic. 
I can give a short summary of my simulation to make it more understandable:

I created a ball and a box which are located in a cylinder and the box
and the box rotates around the z axis enclosed by the spheres. That's
what i said before.

oriBody = Quaternion(Vector3(0,0,1),(math.pi))
box=O.bodies.append(geom.facetBox((.5,0,.5),(.3,.3,.3),oriBody, wallMask=63, color=(0,0,1), dynamic=True, material='mat_box'))

cylinder=O.bodies.append(utils.geom.facetCylinder((0,0,1),1.5,2,orientation=Quaternion((1,
0, 0),0),segmentsNumber=20,wallMask=6, color=(0,1,0), fixed=True))

sphereRadius = 0.06
nbSpheres = (10,10,30)
#nbSpheres = (50,50,50)
for i in xrange(nbSpheres[0]):
    for j in xrange(nbSpheres[1]):
        for k in xrange(nbSpheres[2]):
            x = (i*2 - nbSpheres[0])*sphereRadius*1.1
            y = (j*2 - nbSpheres[1])*sphereRadius*1.1
            z = (k*2 - nbSpheres[2])*sphereRadius*1.1
            sp1=utils.sphere([x,y,z+10],sphereRadius,color=(1,0,0),material='mat_spheres')
            O.bodies.append(sp1)

O.engines=[
   ForceResetter(),
   InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
   InteractionLoop(
      [Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
      [Ip2_FrictMat_FrictMat_FrictPhys()],
      [Law2_ScGeom_FrictPhys_CundallStrack()]

   ),
  NewtonIntegrator(damping=0.4, gravity=(0,0,-9.81)),
	RotationEngine(rotateAroundZero=True,zeroPoint=(0,0,0),rotationAxis=(0,0,1),angularVelocity=100*(2*pi/60),ids=box,label='rotor')
]

O.dt=0.2*utils.PWaveTimeStep()

Should i add the GlobalStiffnessTimeStepper to my code?

Christian

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


Follow ups