Respected everyone,
I am trying to simulate an UCS sample with crack in PFC. The ball is getting deleted with in the specified region but further compaction the ball from is coming in between the region.
The code is attached with the message and problem of the image.
; fname: make_sample.p2dat
;
; Create unbonded sample
;=========================================================================
model new
model largestrain on
model title ‘Testing Bonded Particle Model’
; Set the domain extent
model domain extent -0.05 0.05 -0.1 0.1 condition destroy
contact cmat default model linear method deform emod 1.0e9 kratio 0.0
contact cmat default property dp_nratio 0.5
; create walls that extend past the edges of the sample
wall create vertices -0.025,0.038 0.025,0.038 id 1
wall create vertices -0.025,-0.038 0.025,-0.038 id 2
wall create vertices -0.019,-0.040 -0.019,0.040 id 3
wall create vertices 0.019,-0.040 0.019,0.040 id 4
wall create vertices -0.01,0.0075 0.01,0.0075 id 6
wall create vertices -0.01,-0.0075 0.01,-0.0075 id 7
wall create vertices -0.01,-0.0075 -0.01,0.0075 id 8
wall create vertices 0.01,-0.0075 0.01,0.0075 id 9
wall group ‘w1’ range id 6 9
model random 10001
ball distribute poros 0.35 …
resolution 0.04 …
radius 75e-6 4.75e-3 …
box -0.019 0.019 -0.038 0.038
ball attribute density 2500 damp 0.7
; Calm the system
;model cycle 1000 calm 10
measure create id 1 rad [0.4*(math.min(0.038,0.076))]
[porosity = measure.porosity(measure.find(1))]
; Solve the system to a target limit (here the average force ratio)
; Use density scaling to quickly reach equilibrium
ball delete range rectangle position-x -0.01 0.01 position-y -0.0075 0.0075
model mechanical timestep scale
model solve cycle 5000;ratio-average 1e-4
model mechanical timestep auto
model calm
; delete side walls
wall delete range id 3
wall delete range id 4
wall delete range id 6
wall delete range id 7
wall delete range id 8
wall delete range id 9
model save ‘unbonded’
;=========================================================================
; eof: make_sample.p2dat