Question #680307 on Yade changed:
https://answers.launchpad.net/yade/+question/680307
Status: Answered => Open
Gavin He is still having a problem:
Hello Jérôme,
With your responses, I have been able to generate aabbwalls surrounding
polyhedra specimen.
However, after I modify the code of triax-basic.py, trying to run fast
triaxial test on polyhedrons, I found the aabbwalls cannot run as the
same way as the example does. Eventually I found I need an IPhysFunctor
like "Ig2_Box_Polyhedra_ScGeom" (like "Ig2_Box_Sphere_ScGeom" when
simulating with spheres or clumps) to make the aabbwalls able to
interact with polyhedras, which I cannot find.
So now the question looks to be "how to make aabbwalls able to interact
with polyhedras"?
Could you help me with this problem?
_______________________________________________________________________________________________________________________
CODE:
from yade import pack,polyhedra_utils
gravel = PolyhedraMat()
gravel.density = 2600 #kg/m^3
gravel.young = 1E7 #Pa
gravel.poisson = 20000/1E7
gravel.frictionAngle = 0.5 #rad
mn,mx=Vector3(0,0,0),Vector3(10,10,10)
O.materials.append(PolyhedraMat(young=15e6,poisson=.4,frictionAngle=0,density=0,label='frictionless'))
walls=aabbWalls(extrema=(mn,mx),thickness=1e-10,material='frictionless')
wallIds=O.bodies.append(walls)
polyhedra_utils.fillBox((0,0,0),
(10,10,10),gravel,sizemin=[1,1,1],sizemax=[0.1,0.1,0.1],seed=4)
triax=TriaxialCompressionEngine(
wall_bottom_id=wallIds[2],
wall_top_id=wallIds[3],
wall_left_id=wallIds[0],
wall_right_id=wallIds[1],
wall_back_id=wallIds[4],
wall_front_id=wallIds[5],
internalCompaction=False,
## define the rest of triax params here
## see in pkg/dem/PreProcessor/TriaxialTest.cpp:524 etc
## which are assigned in the c++ preprocessor actually
sigmaIsoCompaction=-50e3,
sigmaLateralConfinement=-50e3,
max_vel=10,
strainRate=0.01,
label="triax"
)
O.engines=[
ForceResetter(),
InsertionSortCollider([Bo1_Polyhedra_Aabb(),Bo1_Wall_Aabb(),Bo1_Facet_Aabb()],verletDist=5e-6),
InteractionLoop(
[Ig2_Wall_Polyhedra_PolyhedraGeom(),
Ig2_Polyhedra_Polyhedra_PolyhedraGeom(), Ig2_Facet_Polyhedra_PolyhedraGeom()],
[Ip2_PolyhedraMat_PolyhedraMat_PolyhedraPhys()], # collision "physics"
[Law2_PolyhedraGeom_PolyhedraPhys_Volumetric()] # contact law --
apply forces
),
GlobalStiffnessTimeStepper(active=1,timeStepUpdateInterval=100,timestepSafetyCoefficient=0.8),
triax,
# you can add TriaxialStateRecorder and such here…
NewtonIntegrator(damping=.4)
]
O.dt=polyhedra_utils.PWaveTimeStep()
--
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 : yade-users@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-users
More help : https://help.launchpad.net/ListHelp
