单轴压缩实验5.0版

; fname: fragment.p3dat
;
; Illustrate how the Fragment, DFN and Model Result logic can be used
; to study the fragmentation of a Bonded Particle Model.
;
;==============================================================================
new
set precision 15
title 'Fragmentation of a Bonded Particle Model'

;------------------------------------------------------------------------------
; --------------- INITIAL PHASE: CREATION OF A BONDED SPECIMEN ----------------
;------------------------------------------------------------------------------

domain extent -5 5 -5 5 -7.5 7.5 condition destroy
cmat default model linear property kn 1e7
cmat default type ball-ball model linearpbond property kn 5e6 proximity 0.1

wall generate box -2.5 2.5 -2.5 2.5 -5 5
set random 10001
ball distribute porosity 0.3 radius 0.25 0.3 box -2.5 2.5 -2.5 2.5 -5 5
ball attribute density 1.0e3 damp 0.7

cycle 2000 calm 10
solve aratio 5e-3 calm 1000

contact method bond gap 0.1
contact property pb_kn 1e7 pb_ks 1e7 pb_ten 4e5 pb_coh 4e5 pb_fa 20.0

wall delete
cmat default type ball-facet model linear property kn 1e7 ks 1e7 fric 0.5

wall generate id 1 name top polygon -4.0 -4.0 5.0 ...
-4.0 4.0 5.0 ...
4.0 4.0 5.0 ...
4.0 -4.0 5.0


wall generate id 2 name bottom polygon -4.0 -4.0 -5.0 ...
-4.0 4.0 -5.0 ...
4.0 4.0 -5.0 ...
4.0 -4.0 -5.0

wall attribute centrotation 0.0 0.0 5.0 range z 5.0
wall attribute centrotation 0.0 0.0 -5.0 range z -5.0
clean
cycle 10
solve arat 5e-3 elastic
wall attribute zvel 0.0 gradient 0.0 0.0 -0.001
history id 1 mechanical age
history id 2 wall contactforce id 1
history id 3 wall contactforce id 2
save initial

;------------------------------------------------------------------------------
; ------------ CASE 0: NO TRACKING OF FRAGMENTS - NO BALL RESULTS -------------
;------------------------------------------------------------------------------

solve time 50.0
save final-case0

;------------------------------------------------------------------------------
; ------------ CASE 1: CONTINUOUS UPDATE OF FRAGMENTS WITH BALL RESULTS --------
;------------------------------------------------------------------------------
restore initial

fragment register ball-ball
fragment activate time 0.1
ball result addattribute fragment time 0.1 activate
wall result addattribute velocity time 0.1 activate

solve time 50.0
save final-case1

;call post_treat.p3fis
;@make_movie_case1

;------------------------------------------------------------------------------
; ------------ CASE 2: CUSTOM UPDATE OF FRAGMENTS AND BALL RESULTS ------------
;------------------------------------------------------------------------------
restore initial

call fracture.p3fis

@track_init
solve time 50.0
@track_end
save final-case2

;call post_treat.p3fis
;@make_movie_case2

return
;==============================================================================
自己看看 我的是PFC5.0