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 [email protected] 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. _______________________________________________ Mailing list: https://launchpad.net/~yade-users Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-users More help : https://help.launchpad.net/ListHelp
