Cpm model
Hello everyone,
I'm using the cpm model to do a simple 2D unconfined compression test. After generating the packing, I load the packing into the compression test in which I erase the other walls except the upper and the lower walls but it seems that there is no cohesion since the packing explodes after loading. I wonder why there is no cohesion within the particles.
Here is my srcipts for making the packing:
from yade import pack
O.materials.
O.materials.
sp = pack.SpherePack()
size = .3
sp.makeCloud(
O.bodies.
walls=aabbWalls
wallIds=
#O.bodies.
#O.bodies.
#O.cell.hSize = Matrix3(size,0,0, 0,size,0, 0,0,.1)
print len(O.bodies)
r = numpy.zeros(400)
zdim = 5.e-2
for i in range(400):
p=O.bodies[i]
p.state.
p.state.mass = 2650 * zdim * pi * p.shape.radius**2 # 0.1 = thickness of cylindrical particle
inertia = 0.5 * p.state.mass * p.shape.radius**2
p.state.inertia = (.5*inertia,
r[i] = p.shape.radius
O.dt = .5*utils.
print O.dt
triax=TriaxialS
wall_
wall_
wall_
wall_
wall_
wall_
internalComp
stressMask = 3, #1 is x, 2 is y, 4 is z, Bitmask, if imposed goal`s are stresses (0 for none, 7 for all, 1 for direction 1, 5 for directions 1 and 3, etc. :ydefault:`7)
goal1 = -1.e7*(
goal2 = -1.e7*(
#goal1 = -1.e7,
#goal2 = -1.e7,
goal3 = 0.,
max_vel=0.01,
label = 'triax'
)
O.engines = [
ForceResetter(),
InsertionSor
InteractionLoop(
[
[
[
),
triax,
NewtonIntegr
]
while 1:
O.run(1000,True)
unb=
area=0.
for i in range(400):
p = O.bodies[i]
area += pi*p.shape.
areatotal = triax.width*
void=
voidratio=
print 'unbalanced force:'
if unb<0.001:
O.
break
any kind of help is appreciated.
Question information
- Language:
- English Edit question
- Status:
- Solved
- For:
- Yade Edit question
- Assignee:
- No assignee Edit question
- Solved by:
- Jan Stránský
- Solved:
- Last query:
- Last reply: