yade-users team mailing list archive

[Luis Barbosa <question291326@xxxxxxxxxxxxxxxxxxxxx>]

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

    Status: Answered => Open

Luis Barbosa is still having a problem:
Hi Jan,

Sure I can.

#!/usr/bin/python
# -*- coding: utf-8 -*-
#======================================
from yade import pack
import gts
from yade import utils
surf=gts.read(open('tetrahedron.gts'))
pred=pack.inGtsSurface(surf)
vol = surf.volume()
#======================================
# Spheres
O.materials.append(JCFpmMat(type=1,young=1e8,poisson=0.3,frictionAngle=radians(30),density=3000,tensileStrength=1e6,cohesion=1e6,jointNormalStiffness=1e7,jointShearStiffness=1e7,jointCohesion=1e6,jointFrictionAngle=radians(20),jointDilationAngle=0.0,label='sphereso'))
#======================================
O.bodies.append(pack.gtsSurface2Facets(surf,color=(0,0,1)))
sp = pack.randomDensePack(pred,radius=0.08,rRelFuzz=0.01, material = 'sphereo',useOBB=False,returnSpherePack=True)
sp.toSimulation(color=(0.9,0.8,0.6))
#======================================
O.engines=[
	ForceResetter(),
	InsertionSortCollider([Bo1_Sphere_Aabb(),Bo1_Facet_Aabb()]),
	InteractionLoop(
		[Ig2_Sphere_Sphere_ScGeom(),Ig2_Facet_Sphere_ScGeom()],
		[Ip2_JCFpmMat_JCFpmMat_JCFpmPhys(cohesiveTresholdIteration=1)],
		[Law2_ScGeom_JCFpmPhys_JointedCohesiveFrictionalPM(smoothJoint=False)]
	),
  GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.2),
  NewtonIntegrator(damping=0.9,gravity=[0,0,-9.81]),
]
ve=utils.getSpheresVolume()
print vol,ve
#=======================================
from yade import qt
qt.View()
qt.Controller()
#########################################################

And here follows the gts surface...
https://drive.google.com/open?id=0B9LhrwAf_vasanZ4bVc4b0FHVUE

Thanks,
Luis

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