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
|
# -*- encoding=utf-8 -*-
#
# https://yade-dem.org/wiki/ScriptParametricStudy
#
## read parameters from table here
from yade import plot
readParamsFromTable(gravity=-9.81, density=2400, initialSpeed=10, noTableOk=True)
from yade.params.table import *
print(gravity, density, initialSpeed)
O.materials.append(FrictMat(young=30e9, density=density, poisson=.3)) ## use the 'density' parameter here
O.engines = [
ForceResetter(),
InsertionSortCollider([
Bo1_Sphere_Aabb(),
Bo1_Box_Aabb(),
]),
InteractionLoop([Ig2_Sphere_Sphere_ScGeom(), Ig2_Box_Sphere_ScGeom()], [
Ip2_FrictMat_FrictMat_FrictPhys(),
], [
Law2_ScGeom_FrictPhys_CundallStrack(),
]),
PyRunner(iterPeriod=100, command='myAddPlotData()', label='plotDataCollector'),
NewtonIntegrator(damping=0.4, gravity=(0, 0, gravity)) ## here we use the 'gravity' parameter
]
O.bodies.append([box([0, 50, 0], extents=[1, 50, 1], fixed=True, color=[1, 0, 0]), sphere([0, 0, 10], 1, color=[0, 1, 0])])
O.bodies[1].state.vel = (0, initialSpeed, 0) ## assign initial velocity
O.dt = .8 * PWaveTimeStep()
## o.saveTmp('initial')
def myAddPlotData():
s = O.bodies[1]
plot.addData({'t': O.time, 'y_sph': s.state.pos[1], 'z_sph': s.state.pos[2]})
plot.plots = {'y_sph': ('z_sph',)}
# run 30000 iterations
O.run(20000, True)
# write some results to a common file
# (we rely on the fact that 2 processes will not write results at exactly the same time)
# 'a' means to open for appending
#file('multi.out','a').write('%s %g %g %g %g\n'%(O.tags['description'],gravity,density,initialSpeed,O.bodies[1].state.pos[1]))
file('multi.out', 'a').write('%g %g %g %g\n' % (gravity, density, initialSpeed, O.bodies[1].state.pos[1]))
print('gnuplot', plot.saveGnuplot(O.tags['id']))
|
