yade-users team mailing list archive

[Christopher Stanbridge <question270352@xxxxxxxxxxxxxxxxxxxxx>]

New question #270352 on Yade:
https://answers.launchpad.net/yade/+question/270352

Hello all,

The Yade script I've posted below should, according to the aabbEnlargeFactor and interactionDetectionFactor values, generate bonds between adjacent spheres. However, it doesn't -- the number of contacts formed by it is zero. Does anyone here know why this might be? All help is appreciated.

- Christopher Stanbridge

import sys, time, random, os, gts, math
from yade import ymport
from yade import pack
from yade import qt
from yade import export

mat1=O.materials.append(CohFrictMat(young=1e8, poisson=1/0.25, density=2000, frictionAngle=radians(30), alphaKr=0.89, alphaKtw=0.67, etaRoll=-1, isCohesive=True, momentRotationLaw=True,normalCohesion=20000,shearCohesion=30000))
sphereIce=O.materials[mat1]

mat2=O.materials.append(CohFrictMat(young=1e8, poisson=1/0.25, density=2000, frictionAngle=radians(30), alphaKr=0.89, alphaKtw=0.67, etaRoll=-1, isCohesive=False, momentRotationLaw=True,normalCohesion=0,shearCohesion=0))
wallIce=O.materials[mat2]

O.bodies.append([
	utils.wall(-1,0,material=wallIce,color=(0.8,0.8,0.8)),
	utils.wall(1,0,material=wallIce,color=(0.8,0.8,0.8)),
	utils.wall(-1,1,material=wallIce,color=(0.8,0.8,0.8)),
	utils.wall(1,1,material=wallIce,color=(0.8,0.8,0.8)),
	utils.wall(-1,2,material=wallIce,color=(0.8,0.8,0.8)),
	utils.wall(1,2,material=wallIce,color=(0.8,0.8,0.8))
])

for x in range(-9, 9):
	for y in range(-9, 9):
		for z in range(-9, 9):
			r = 0.05
			if (x + y + z) % 2 == 1:
				r = 0.045
			O.bodies.append(utils.sphere(center=(x/10.0,y/10.0,z/10.0),radius=r,material=sphereIce,color=(0.8,0.2,0.2),fixed=False))


O.engines=[
	ForceResetter(),
	InsertionSortCollider([	Bo1_Sphere_Aabb(aabbEnlargeFactor=1.3), Bo1_Wall_Aabb()]),
	InteractionLoop([Ig2_Sphere_Sphere_ScGeom6D(interactionDetectionFactor=1.3), Ig2_Wall_Sphere_L3Geom()],
		[Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=False)],
		[Law2_ScGeom6D_CohFrictPhys_CohesionMoment(always_use_moment_law=False, creep_viscosity=1, neverErase=False,shear_creep=False,twist_creep=False,useIncrementalForm=False)]
	),

	NewtonIntegrator(damping=0.2, gravity=[0, 0, 0], label='integrator')
]



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