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
|
# -*- encoding=utf-8 -*-
"""
Same example as CohesiveCylinderSphere.py but using gridConnections instead of chainedCylinder.
"""
from yade import qt
from yade.gridpfacet import *
from numpy import linspace
### Engines need to be defined first since the function gridConnection creates the interaction
O.engines = [
ForceResetter(),
InsertionSortCollider([
Bo1_Sphere_Aabb(),
Bo1_GridConnection_Aabb(),
]),
InteractionLoop(
# Geometric interactions
[
Ig2_GridNode_GridNode_GridNodeGeom6D(),
Ig2_Sphere_GridConnection_ScGridCoGeom(), # used for the cohesive sphere-cylinder interaction
],
[
# Interaction phusics
Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True, setCohesionOnNewContacts=False),
],
# Interaction law
[
Law2_ScGeom6D_CohFrictPhys_CohesionMoment(),
Law2_ScGridCoGeom_CohFrictPhys_CundallStrack(), # used for the cohesive sphere-cylinder interaction
]
),
NewtonIntegrator(gravity=(0, 0, 0), damping=0.3, label='newton'),
PyRunner(command='main()', iterPeriod=1000),
]
O.dt = 5e-07
O.materials.append(
CohFrictMat(
young=8e5,
poisson=0.3,
density=4e3,
frictionAngle=radians(30),
normalCohesion=1e5,
shearCohesion=1e5,
momentRotationLaw=True,
label='gridNodeMat'
)
)
O.materials.append(
CohFrictMat(
young=8e5,
poisson=0.3,
density=4e3,
frictionAngle=radians(30),
normalCohesion=1e5,
shearCohesion=1e5,
momentRotationLaw=True,
label='gridCoMat'
)
)
O.materials.append(
CohFrictMat(
young=8e5,
poisson=0.3,
density=4e3,
frictionAngle=radians(30),
normalCohesion=4e4,
shearCohesion=1e5,
momentRotationLaw=False,
label='spheremat'
)
)
rCyl = 0.006
nL = 15
L = 0.3
### Create all nodes first :
nodesIds = []
for i in linspace(0, L, nL):
nodesIds.append(O.bodies.append(gridNode([i, 0, 0], rCyl, wire=False, fixed=False, material='gridNodeMat')))
### Now create connection between the nodes
for i, j in zip(nodesIds[:-1], nodesIds[1:]):
O.bodies.append(gridConnection(i, j, rCyl, material='gridCoMat'))
### Set fixed nodes
O.bodies[nodesIds[0]].state.blockedDOFs = 'xyzXYZ'
O.bodies[nodesIds[-1]].state.blockedDOFs = 'xyzXYZ'
IdSphere = O.bodies.append(sphere([0.15, 0, 2. * rCyl], rCyl, wire=False, fixed=False, material='spheremat'))
def main():
global Fn, Ft
if O.iter > 50000:
O.bodies[IdSphere].dynamic = False
O.bodies[IdSphere].state.vel[2] = 0.1
qt.View()
O.saveTmp()
|
