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
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
|
# -*- encoding=utf-8 -*-
# The script implements the show case in this article [Mani2013]
from yade import utils, plot
o = Omega()
fr = 0.5
rho = 2000
tc = 0.001
en = 0.7
et = 0.7
o.dt = 1.0
r1 = 1.0
r2 = 1.0
Gamma = 20.6 * 1e-3
Theta = 0
VB3 = 74.2 * 1e-12
CapillarType = "Lambert"
mat1 = O.materials.append(
ViscElCapMat(frictionAngle=fr, density=rho, Vb=VB3, gamma=Gamma, theta=Theta, Capillar=True, CapillarType=CapillarType, tc=tc, en=en, et=et)
)
d = 0.9999
id1 = O.bodies.append(sphere(center=[0, 0, 0], radius=r2, material=mat1, fixed=True, color=[0, 1, 0]))
id2 = O.bodies.append(sphere(center=[0, (r1 + r2) * d, 0], radius=r2, material=mat1, fixed=True, color=[0, 1, 0]))
id3 = O.bodies.append(sphere(center=[(r1 + r2) * d, (r1 + r2) * d, 0], radius=r2, material=mat1, fixed=True, color=[0, 1, 0]))
id4 = O.bodies.append(sphere(center=[(r1 + r2) * d, (r1 + r2) * d * 2, 0], radius=r2, material=mat1, fixed=True, color=[0, 1, 0]))
id5 = O.bodies.append(sphere(center=[(r1 + r2) * d * 2, (r1 + r2) * d, 0], radius=r2, material=mat1, fixed=True, color=[0, 1, 0]))
Vf = 0.0e-1
Vfmin = 0.0e-1
O.bodies[id1].state.Vf = Vf
O.bodies[id1].state.Vmin = Vfmin
O.bodies[id2].state.Vf = Vf
O.bodies[id2].state.Vmin = Vfmin
O.bodies[id3].state.Vf = Vf
O.bodies[id3].state.Vmin = Vfmin
O.bodies[id4].state.Vf = Vf
O.bodies[id4].state.Vmin = Vfmin
O.bodies[id5].state.Vf = Vf
O.bodies[id5].state.Vmin = Vfmin
vel = 0.0
O.bodies[id1].state.vel = [0, 0, vel]
O.bodies[id2].state.vel = [0, 0, -vel]
O.bodies[id3].state.vel = [vel, 0, 0]
O.bodies[id4].state.vel = [-vel, 0, 0]
o.engines = [
ForceResetter(),
InsertionSortCollider([Bo1_Sphere_Aabb()], verletDist=(r1 + r2) * 5.0),
InteractionLoop(
[Ig2_Sphere_Sphere_ScGeom()],
[Ip2_ViscElCapMat_ViscElCapMat_ViscElCapPhys()],
[Law2_ScGeom_ViscElCapPhys_Basic()],
),
LiqControl(label='lqc'),
NewtonIntegrator(damping=0, gravity=[0, 0, 0]),
PyRunner(command='showData()', iterPeriod=1),
]
def showData():
print("Step %d" % O.iter)
print("idB=%d, Vf=%s, Vmin=%s;" % (id1, O.bodies[id1].state.Vf, O.bodies[id1].state.Vmin))
print("idB=%d, Vf=%s, Vmin=%s;" % (id2, O.bodies[id2].state.Vf, O.bodies[id2].state.Vmin))
print("idB=%d, Vf=%s, Vmin=%s;" % (id3, O.bodies[id3].state.Vf, O.bodies[id3].state.Vmin))
print("idB=%d, Vf=%s, Vmin=%s;" % (id4, O.bodies[id4].state.Vf, O.bodies[id4].state.Vmin))
print("idB=%d, Vf=%s, Vmin=%s;" % (id5, O.bodies[id5].state.Vf, O.bodies[id5].state.Vmin))
try:
print("Interaction[1, 2].Vb=%s" % (O.interactions[id1, id2].phys.Vb))
except:
pass
try:
print("Interaction[2, 3].Vb=%s" % (O.interactions[id2, id3].phys.Vb))
except:
pass
try:
print("Interaction[3, 4].Vb=%s" % (O.interactions[id3, id4].phys.Vb))
except:
pass
try:
print("Interaction[3, 5].Vb=%s" % (O.interactions[id3, id5].phys.Vb))
except:
pass
print()
showData()
O.run(1, True)
for i in range(5):
O.bodies[i].state.Vf = 0
O.bodies[i].state.Vmin = 0
O.interactions[id1, id2].phys.Vmax = 5.0
lqc.addLiqInter(id1, id2, 1.0)
O.interactions[id2, id3].phys.Vmax = 5.0
lqc.addLiqInter(id2, id3, 1.0)
O.interactions[id3, id4].phys.Vmax = 5.0
lqc.addLiqInter(id3, id4, 1.0)
O.interactions[id3, id5].phys.Vmax = 5.0
lqc.addLiqInter(id3, id5, 1.0)
O.run(1, True)
vel = 1.0
O.bodies[id3].state.vel = [vel, 0, 0]
O.bodies[id4].state.vel = [vel, 0, 0]
O.bodies[id5].state.vel = [vel, 0, 0]
from yade import qt
qt.View()
|
