yade-users team mailing list archive

[Luis Barbosa <question695592@xxxxxxxxxxxxxxxxxxxxx>]

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

Hi all,

I am used to create pakings from gts surfaces. However, I am having problem with this right now.
I am trying to create a dense pack in gts predicates, but it seems it is operating in an inverse way. Instead of the pack be presented inside the predicate, it appears outside. The predicate is the void region inside the pack. It seems a boolean operation, but has nothing to do with my code. I already tried creating gts binary or ascii, or even changing the language ( -*- coding: utf-8 -*-).

Here is the script I am using

#!/usr/bin/python
# -*- coding: utf-8 -*-
# F.T. Silva grupodemfeagri@xxxxxxxxx

from yade.pack import *
from yade import pack,export,qt
from math import *
from yade import plot
from yade import ymport
import gts

######################################Material#################################################
O.materials.append(JCFpmMat(type=1,young=70e5,poisson=0.3,frictionAngle=radians(30),density=0.0025,tensileStrength=194,cohesion=194,jointNormalStiffness=10,jointShearStiffness=10,jointCohesion=194,jointFrictionAngle=radians(30),jointDilationAngle=0.0,label='spheres'))
######################################Selector############################################
sel=22
######################################Surface#############################################
surf=gts.read(open('agg'+str(sel)+'.gts'))

if surf.is_closed():
	pred=pack.inGtsSurface(surf)
	radius=0.03 #(cm)
	O.bodies.append(pack.randomDensePack(pred,radius=radius))

O.bodies.append(pack.gtsSurface2Facets(surf,wire=True))


######################################Engine#################################################
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, Key="Wei", recordCracks=True)]
	),

  GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=500,timestepSafetyCoefficient=0.5),
  NewtonIntegrator(damping=0.1,gravity=[0,0,-981]),#cm/s2

]

from yade import qt
qt.View()
qt.Controller()



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