model rest 'TBM_high_stress' ; fish define elas_props _bulk = _young/3.0/(1-2*_poiss) _shear = _young/2.0/(1+_poiss) end ;Contitinum material paramters ; SS-I Transversely Isotropic [global _density_SSI = 2400] [global _dip_SSI = 0] [global _dipdir_SSI = 0] [global _youngnormal_SSI = 1.5e9] [global _shearnormal_SSI = 0.3e9] [global _youngplane_SSI = 4.5e9] [global _poisson_SSI = 0.2] ; ; ;Discontitinum material paramters ;BSS-I - STRAINSOFTEN CWFS FOR SANDSTONE [global _density_BSSI = 2400] [global _young_BSSI = 5.5e9] [global _poisson_BSSI = 0.2] [global _coh_peak_BSSI = 4500e3] [global _chesion_res_BSSI = 100e3] [global _friction_peak_BSSI = 20] [global _friction_res_BSSI = 55] [global _ten_peak_BSSI = 2250e3] ; lab scale tensile strength [global _ten_res_BSSI = 20e3] [global _dilation_peak_BSSI = 12] [global _dilation_res_BSSI = 0] ; ; ; ASSIGN MODELS - contitinum zones ; ;Transversely Isotropic Zones ;SS-I block zone cmodel assign anisotropic range group 'SS-I' block zone property density @_density_SSI dip @_dip_SSI dip-direction @_dipdir_SSI & young-normal @_youngnormal_SSI poisson-normal @_poisson_SSI shear-normal @_shearnormal_SSI & young-plane @_youngplane_SSI poisson-plane @_poisson_SSI range group 'SS-I' ; ; ASSIGN MODELS - discontitinum zones ; ;B-SS-I fish set @_young = @_young_BSSI fish set @_poiss = @_poisson_BSSI @elas_props ; block zone cmodel assign strain-softening range group 'B-SS-I' block zone prop dens @_density_BSSI shear @_shear bulk @_bulk cohesion @_coh_peak_BSSI friction @_friction_peak_BSSI ... tension @_ten_peak_BSSI dilation @_dilation_peak_BSSI table-cohesion 1 table-friction 2 table-tension 3 table-dilation 4 range group 'B-SS-I' ; ;DEFINE TABLES FOR STRAIN SOFTENING MODELS - BSS-II ;cohesion tables table 1 add 0 @_coh_peak_BSSI 0.002 @_chesion_res_BSSI ;friction tables table 2 add 0 @_friction_peak_BSSI 0.005 @_friction_res_BSSI ;tension tables table 3 add 0 @_ten_peak_BSSI 0.01 @_ten_res_BSSI ;dilation table for all units table 4 add 0 @_dilation_peak_BSSI 0.01 @_dilation_res_BSSI ; ; ; ;Discontinuities ; ; default contact props block contact cmodel default 2 material 3 block contact property material 3 stiff-norm 4e9 stiff-shear 4e8 friction 30 ; ; make construction joints elastic block contact cmodel assign area range jointset-id 50 block contact prop stiff-norm 1e11 stiff-shear 1e11 range jointset-id 50 ;return ; ; ;Sandstone joints ;bedding block contact cmodel assign residual range jointset-id 20 block contact prop stiff-norm 6e9 stiff-shear 6e8 friction 32 friction-residual 30 cohesion 10e3 cohesion-residual 0 tension 0 range jointset-id 20 ;id-list ;Joints block contact cmodel assign residual range jointset-id 30 block contact prop stiff-norm 10e9 stiff-shear 10e8 friction 40 friction-residual 30 cohesion 0e3 tension 0 range jointset-id 30 ; ;Initial boundary conditions block gridpoint apply vel-x 0 range pos-x -60 pos-y -100 28.2 block gridpoint apply vel-x 0 range pos-x 60 pos-y -100 28.2 block gridpoint apply vel-y 0 range pos-x -60 60 pos-y -100 ; ; model gravity 0 -9.81 ;boundary stress ;block edge apply stress 0, 0, -25e3 range position-x -95,95 position-y 28.2 ; ;in-situ stress block insitu stress -3771.8e3 0 -676.8e3 gradient-y 91.2e3 0 24.0e3 ... stress-ZZ -2300.8e3 gradient-z 0 55.6e3 range position-y -100 28.2 ; ; model solve model save 'initial1.sav' block gridpoint initialize displacement-x 0 block gridpoint initialize displacement-y 0 block contact reset displacement model save 'initial.sav'