1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
# -*- encoding=utf-8 -*-
"""
Show basic wall functionality (infinite axis-aligned planes).
"""
O.materials.append(FrictMat(young=30e9, density=1000, poisson=.2, frictionAngle=.5))
O.bodies.append([wall(1, axis=2, sense=-1), wall(-5, axis=0, sense=1), wall(1, axis=1), wall((1, 0, 0), 0), sphere([0, 0, 0], .5), sphere([-4, -4, -3], .5)])
Gl1_Wall(div=10)
from yade import qt
qt.Controller()
qt.View()
O.engines = [
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb(), Bo1_Facet_Aabb(), Bo1_Wall_Aabb()]),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom(), Ig2_Facet_Sphere_ScGeom(), Ig2_Wall_Sphere_ScGeom()],
[Ip2_FrictMat_FrictMat_FrictPhys()],
[Law2_ScGeom_FrictPhys_CundallStrack()],
),
NewtonIntegrator(damping=0.01, gravity=[1e2, 1e2, 1e2]),
]
O.dt = PWaveTimeStep()
O.save('/tmp/a.xml')
O.saveTmp()
O.run()
#O.bodies.append([
# facet([[-1,-1,0],[1,-1,0],[0,1,0]],dynamic=False,color=[1,0,0],young=1e3),
# facet([[1,-1,0],[0,1,0,],[1,.5,.5]],dynamic=False,young=1e3)
#])
#import random
|
