← Back to team overview

yade-users team mailing list archive

Re: [Question #404135]: intRadius(aabbenlargefactor) doesn't work

 

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

weijy gave more information on the question:
Forgot to paste the script. The testing script mentioned in #4 is like
this:

**********************************
O.materials.append(CohFrictMat(young = 1e6, poisson = 0.5,normalCohesion = 500000, shearCohesion = 500000, frictionAngle = 0,density = 2700,label='spheres'))

O.bodies.append([
    utils.sphere(center=(0.01,0,0.01),radius=0.008,material = 'spheres'),
    utils.sphere(center=(-0.01,0,0.01),radius=0.008,material = 'spheres'),
    utils.sphere(center=(0.01,0,0.03),radius=0.008,material = 'spheres'),
    utils.sphere(center=(-0.01,0,0.03),radius=0.008,material = 'spheres'),
    utils.sphere(center=(0.01,0,-0.01),radius=0.01,fixed=True,material = 'spheres'),
    utils.sphere(center=(-0.01,0,-0.01),radius=0.01,fixed=True,material = 'spheres'),
])

O.engines=[
    ForceResetter(),
    InsertionSortCollider([Bo1_Sphere_Aabb(aabbEnlargeFactor = 1.5,label = 'bo1s')]),
    InteractionLoop(
        [Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor = 1.5,label = 'ig2ss')],
        [Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(label = 'ip2cc')],
        [Law2_ScGeom6D_CohFrictPhys_CohesionMoment()]
    ),
    NewtonIntegrator(damping=0.1),
    PyRunner(command = 'stick()',iterPeriod = 10000,label = 'sticker')
]

ip2cc.setCohesionNow = True
O.step()
num_intr = 0
print "current step is ", O.iter
for i in O.interactions:
    num_intr = num_intr + 1
print 'interaction number is ', num_intr
num_intr = 0
        
def stick():
    global num_intr
    print "current step is ", O.iter
    for i in O.interactions:
        num_intr = num_intr+1
    print "interaction number is ", num_intr
    num_intr = 0
    sticker. iterPeriod = 100
    O.pause()

O.dt=0.005*utils.PWaveTimeStep()
*****************************************

-- 
You received this question notification because your team yade-users is
an answer contact for Yade.