; fname: cmat3.dat (2D)
;
;  Illustrate advanced usage of the CMAT
;  (Default and Optional slots filled with different data)
;=========================================================================
model new
model large-strain on
model random 10001

[k_set = 1.0e8]
[c_set = 1.0e8]
[f_set = 0.0  ]
[rad   = 0.05 ]
[un    = 0.01 ]

model domain extent -1.5e-2 1.5e-2 condition destroy
wall generate box -1.0e-2 1.0e-2 ;one-wall
ball generate number 100 radius 0.8e-3 0.8e-3 box -1.0e-2 1.0e-2 group 'glass'
ball generate number 100 radius 0.8e-3 0.8e-3 box -1.0e-2 1.0e-2 group 'steel'
wall group 'Bot' facet range position-z -0.01 -0.01
wall group 'South' facet range position-y -0.01 -0.01
wall group 'North' facet range position-y 0.01 0.01
wall group 'East' facet range position-x -0.01 -0.01
wall group 'West' facet range position-x 0.01 0.01

ball attribute density 2500.0 range group 'glass'
ball attribute density 7500.0 range group 'steel'

; excerpt-oohc-start
contact cmat add 1 model hertz                          ...
                 property hz_shear 30e9 hz_poiss 0.3    ...
                          fric 0.25                     ...
                          dp_nratio 0.2                 ...
             range group 'glass' matches 2

contact cmat add 2 model hertz                          ...
                 property hz_shear 70e9 hz_poiss 0.3    ...
                          fric 0.05                     ...
                          dp_nratio 0.2                 ...
             range group 'steel' matches 2

contact cmat default type ball-ball                     ...
                   model hertz                          ...
                   property hz_shear 50e9 hz_poiss 0.3  ...
                           fric 0.05                    ...
                           dp_nratio 0.2

contact cmat default type ball-facet                    ...
                   model linear                         ...
                   property kn 1e6 dp_nratio 0.2

;This block was addded to define different properties for the bottom wall 


contact cmat add 3 model linear                        ...
                property kn 3e6 dp_nratio 0.2    ...
                range group 'North' matches 1
                
contact cmat add 4 model linear                        ...
                property kn 4e6 dp_nratio 0.2    ...
                range group 'South' matches 1
                
contact cmat add 5 model linear                        ...
                property kn 5e6 dp_nratio 0.2    ...
                range group 'East' matches 1

contact cmat add 6 model linear                        ...
                property kn 6e6 dp_nratio 0.2    ...
                range group 'West' matches 1
                 
contact cmat add 7 model burger property ...
                        'bur_knk'  [k_set] 'bur_cnk' [c_set] ...
                        'bur_knm'  [k_set] 'bur_cnm' [c_set] ...
                        'bur_ksk'  [k_set] 'bur_csk' [c_set] ...
                        'bur_ksm'  [k_set] 'bur_csm' [c_set] ...
                        'bur_fric' [f_set] ...
                        range group 'Bot' matches 1




; excerpt-oohc-end

model gravity 9.81
model solve
model save 'cmat3'

program return
;=========================================================================
; eof: cmat3.dat (2D)