;fname:flat_joint.p3dat
;
;Create sample with flat joint contacts
model restore 'unbonded sample'

;here we want all existing contacts to assume flat_joint properties
contact model flatjoint
contact property fj_track on

;bond ball-ball contacts
contact method bond gap 0 range contact type 'ball-ball'

;set stiffness of bonds using methods
contact method deform emod 60e9 krat 1.5

;set bond strengths
contact property fj_ten 9.40e6 fj_coh 73e6 fj_fa 53

;set ball-ball friction to non-zero value
contact property fj_fric 0.577 range contact type 'ball-ball'

;Set new ball-ball contacts to be flatjoint contacts (unbonded)
contact cmat default model flatjoint...
           property fj_track on...
           method deform emod 60e9 krat 1.5
           type 'ball-ball'

;Reset ball displacement
ball attribute displacement multiply 0.0
;Set the linear force to 0.0 and a force a reset of the linear contact forces.
contact property lin_force 0.0 0.0 0.0 lin_mode 1
ball attribute force-contact multiply 0.0 moment-contact multiply 0.0
model cycle 1
model solve ratio-average 1e-5

model save 'flatjoint'

geometry import 'DiscCutter.stl' positon 0 0 0 

;================================================================================
;eof: flat_joint.p3dat