Nonlinear Joint Model In 3DEC
The Nonlinear model is the same as the MohrCoulomb Joint Model In 3DEC except that the normal and shear stiffness may depend on the normal stress as in the Continuously Yielding Joint Model In 3DEC. The tiffness behavior differs from the ContinuouslyYielding Joint Model in that during unloading, the stiffness is constant and is equal to the stiffness at the point of unloading.
Formulations
The response to normal loading is expressed incrementally as
where the normal stiffness, \(k^{current}_n\), is given by
representing the observed increase of stiffness with normal stress, where \(k_n\) and \(e_n\) are input model parameters.
If the normal stress descreases (becomes less compressive), then the normal stiffness is constant and equal to the normal stiffness at the start of unloading. The stiffness will remain constant as the joint is reloaded, until the main loading path is rejoined. This is illustrated in Figure 1.
A maximim and minimum value of normal stiffness may also be supplied by the user. If not provided, these will default to \(k_n\).
The shear stiffness can also vary with normal stress. The relationship is essentially the same:
where \(k_s\) and \(e_s\) are input model parameters. A maximim and minimum value of shear stiffness may also be supplied by the user. If not provided, these will default to \(k_s\).
Summary of Nonlinear Parameters
The model parameters associated with the nonlinear model are summarized in Table 1. The model is accessed in 3DEC with the block contact jmodel assign
nonlinear
command.
Parameter 
Description 
Keyword 

\(k_n\) 
initial joint normal stiffness (STRESS/LENGTH) 

\(k^{max}_n\) 
maximum joint normal stiffness (STRESS/LENGTH). Defaults to k_n. 

\(k^{min}_n\) 
minimum joint normal stiffness (STRESS/LENGTH). Defaults to k_n. 

\(k^{current}_n\) 
current calculated joint normal stiffness (read only) 

\(e_n\) 
joint normal stiffness exponent 

\(k_s\) 
initial joint shear stiffness (STRESS/LENGTH) 

\(k^{max}_s\) 
maximum joint shear stiffness (STRESS/LENGTH). Defaults to k_s. 

\(k^{min}_s\) 
minimum joint shear stiffness (STRESS/LENGTH). Defaults to k_s. 

\(k^{current}_s\) 
current calculated joint shear stiffness (read only) 

\(e_s\) 
joint shear stiffness exponent 

\(\phi\) 
friction angle (DEGREES) 

\(\phi_{res}\) 
residual friction angle (DEGREES). Defaults to phi. 

\(c\) 
cohesion (STRESS) 

\(c_{res}\) 
residual cohesion (STRESS). Defaults to 0. 

\(\psi\) 
dilation angle (DEGREES) 

\(u_{cs}\) 
shear displacement at which dilation stops (LENGTH). Default is infinity 

\(T_f\) 
tensile strength (STRESS) 

\(T_f^{res}\) 
residual tensile strength (STRESS). Default is 0. 

Example: Normal Loading
The following example illustrates the behavior with a simple test model. The problem applies a varying normal load to the top block of a twoblock system.
To view this example in 3DEC, use the menu command . Choose “3DEC/ Constitutive_Models/ Nonlinear” and select “nonlinear.prj” to load. The data file used is shown at the end of this example.
The modeled joint is first subjected to an increasing normal load. This is followed to by an unloadreload excursion to test the functionality of the nonlinear joint normal stiffness calculation. The shear stiffness is kept constant. Figure 2 shows the model; the joint is defined by one contact that is composed of 10 subcontacts. The properties used in this test are:
Property Keyword 
Param. 
Description 
Value 


\(\rho\) 
block mass density 
2600 kg/m^{3} 

\(K\) 
bulk modulus of block 
4 GPa 

\(G\) 
shear modulus of block 
3 GPa 

\(k_n\) 
joint initial normal stiffness 
10 GPa/m 

\(k_s\) 
joint shear stiffness 
10 GPa/m 

\(e_n\) 
joint normal stiffness exponent 
1.1 

\(k^{max}_n\) 
maximum joint normal stiffness 
1000 GPa/m 

\(\phi\) 
joint intrinsic friction angle 
30º 

\(R\) 
joint roughness parameter 
0.1 mm 

\(\psi\) 
dilation angle 
15º 
A graph of normal stress verus normal dispacement during the first loading step is shown in Figure 3. Initially the line is linear with a slope equal to the input value of k_n. After about 0.15 mm of displacement, the stiffness increases (the line becomes steeper).
The joint is then unloaded and the response is shown in Figure 4. The stiffness remains constant and equal to the tangent stiffness at the start of unloading.
Finally, the joint is reloaded (Figure 5). The normal stress increases back up along the unload line until rejoining the main loading curve, at which time the stiffness starts to increase again.
The change in normal stiffness throughout the test is shown in Figure 6.
Data File
compression.dat
model new
model largestrain on
;
block create brick 0.15,0.15 0.10,0.10 0.10,0
block create brick 0.10,0.10 0.10,0.10 0,0.10
block zone generate edgelength 0.2
block zone cmodel assign elastic
block zone prop dens=0.0026 bulk=4000 shear=3000
;
; Nonlinear slip model
block contact jmodel assign nonlinear
block contact prop stiffnessnormal=1e4 stiffnessshear=1e4 ...
friction=30.0 exponentnormal 1.1 knmaximum 1e6
;
block contact mattab default jmodel nonlinear
block contact materialtable default prop stiffnessnormal=1e4 ...
stiffnessshear=1e4 ...
friction=30.0 ...
exponentnormal 1.1 ...
knmaximum 1e6 ...
dilation 15
block hide range posz 0 1
block gridpoint apply vel 0 0 0 range posz 1 0.1
block hide off
;; test insitu
;block insitu stress 0 0 10 0 0 0
;
; normal load
block gridpoint group 'top' range posz 0.1
block gridpoint apply velz 0.01 range group 'top'
fish def kn
cx = block.subcontact.near(0.1,0.1,0)
kn = block.subcontact.prop(cx,'kncurrent')
end
fish his kn
block contact his stressnormal pos 0.1 0.1 0
block contact his dispnorm pos 0.1 0.1 0
;
model cycle 5000
model save 'load'
; now unload
block gridpoint apply velz 0.01 range group 'top'
model cyc 2000
; test save/restore
model save 'unload'
model rest 'unload'
; now reload
block gridpoint apply velz 0.01 range group 'top'
model cyc 3000
model save 'reload'
Example: Shear Loading
The same model as above is loaded in shear. An insitu normal stress of 10 MPa is first applied. This results in an initial normal stiffness of 31.6 GPa/m. The top of the model is fixed during shear loading. This causes the normal stress on the joint to increase once dilation starts. This increase in normal stress causes the normal stiffness to increase. A plot of shear stress verusus shear displacement is shown in Figure 7 and the variation in normal stiffness is shown in Figure 8
Data File
shear.dat
model new
model largestrain on
;
block create brick 0.15,0.15 0.10,0.10 0.10,0
block create brick 0.10,0.10 0.10,0.10 0,0.10
block zone generate edgelength 0.2
block zone cmodel assign elastic
block zone prop dens=0.0026 bulk=4000 shear=3000
;
; Nonlinear slip model
block contact jmodel assign nonlinear
block contact prop stiffnessnormal=1e4 stiffnessshear=1e4 ...
exponentnormal 1.5 knmaximum 1e6 ...
friction=30.0 dilation 15
block contact mattab default jmodel nonlinear
block contact materialtable default prop stiffnessnormal=1e4 ...
stiffnessshear=1e4 ...
exponentnormal 1.5 knmaximum 1e6 ...
friction=30.0 dilation 15
block hide range posz 0 1
block gridpoint apply vel 0 0 0 range posz 1 0.1
block hide off
;
; normal load
block insitu stress 0 0 10 0 0 0
block gridpoint apply velz 0 range posz 0.1
;
model solve
;
block contact reset disp
; shear load
block hide range posz .1 0.0
block gridpoint apply velx=0.005 range posz .1 1.1
block gridpoint apply vely 0 range posz 1 1
block hide off
fish def kn
cx = block.subcontact.near(0.1,0.1,0)
kn = block.subcontact.prop(cx,'kncurrent')
end
fish his kn
block contact his stressnormal pos 0.1 0.1 0
block contact his dispnorm pos 0.1 0.1 0
block contact hist dispshear pos 0.1 0.1 0
block contact hist stressshear pos 0.1 0.1 0
model cyc 15000
model save 'shear'
program return
Example: Shear Loading with Nonlinear Shear Stiffness
The model above is used but a changing normal stress is applied to the top boundary to observe how the shear stiffness is affected. An insitu normal stress of 1 MPa is first applied. This results in an initial shear stiffness of 10 GPa/m (the minimum). The joint is then sheared for 1000 steps. The normal stress is then increased in increments of 1 MPa with shearing after each normal stress change. This increase in normal stress causes the shear stiffness to increase. After 4000 steps, the normal stress is decreased and then increased again. The variation in shear stiffness is shown in Figure 9. You can see how the shear stiffness increases with each increase in normal stress, but between 4000 and 6000 steps, the shear stiffness doesn’t change as the joint expeirneces unloading and reloading. Finally, normal stiffness is increased further and the shear stiffness reaces its maximum value (100 GPa/m). A plot of shear stress verusus shear displacement is shown in Figure 10.
Data File
shear2.dat
model new
model largestrain off
;
block create brick 0.15,0.15 0.10,0.10 0.10,0
block create brick 0.10,0.10 0.10,0.10 0,0.10
block zone generate edgelength 0.2
block zone cmodel assign elastic
block zone prop dens=0.0026 bulk=4000 shear=3000
;
; Nonlinear slip model
block contact jmodel assign nonlinear
block contact prop stiffnessnormal=1e4 stiffnessshear=1e4 ...
exponentshear 1.5 ksmaximum 1e5 ...
friction=30.0 dilation 15
block hide range posz 0 1
block gridpoint apply vel 0 0 0 range posz 1 0.1
block hide off
;
; normal load
block insitu stress 0 0 1 0 0 0
block face apply stresszz 1 range posz 0.1
;
model solve
;
block contact reset disp
; shear load
block hide range posz .1 0.0
block gridpoint apply velx=0.0005 range posz .1 1.1
block gridpoint apply vely 0 range posz 1 1
block hide off
fish def kn
cx = block.subcontact.near(0.1,0.1,0)
kn = block.subcontact.prop(cx,'kncurrent')
ks = block.subcontact.prop(cx,'kscurrent')
end
fish his kn
block contact his stressnormal pos 0.1 0.1 0
block contact his dispnorm pos 0.1 0.1 0
block contact hist dispshear pos 0.1 0.1 0
block contact hist stressshear pos 0.1 0.1 0
fish his ks
model cyc 1000
block face apply stresszz 1 range posz 0.1 ; total = 2
model cyc 1000
block face apply stresszz 1 range posz 0.1 ; total = 3
model cyc 1000
block face apply stresszz 1 range posz 0.1 ; total = 4
model cyc 1000
; unload
block face apply stresszz 2 range posz 0.1
model cyc 1000
; reload
block face apply stresszz 2 range posz 0.1
model cyc 1000
; get to the max
block face apply stresszz 5 range posz 0.1 ; total = 5
model cyc 1000
model save 'shear'
program return
Properties
 cohesion f
Joint cohesion.
 cohesionresidual f
Residual joint cohesion.
 dilation f
Joint dilation angle in degrees.
 dilationzero f
Shear displacement at which joint will cease to dilate.
 exponentnormal f
Joint normal stiffness exponent
 exponentshear f
Joint shear stiffness exponent
 friction f
Friction angle in degrees.
 frictionresidual f
Residual joint friction angle in degrees.
 stiffnessnormal f
Initial joint normal stiffness (stress / distance)
 stiffnessshear f
Joint shear stiffness (stress / distance)
 kncurrent f
Current calculated joint normal stiffness (stress / distance) (read only)
 knmaximum f
Maximum joint normal stiffness (stress / distance)
 knminimum f
Minimum joint normal stiffness (stress / distance)
 kscurrent f
Current calculated joint shear stiffness (stress / distance) (read only)
 ksmaximum f
Maximum joint shear stiffness (stress / distance)
 ksminimum f
Minimum joint shear stiffness (stress / distance)
 tension f
Tensile strength.
 tensionresidual f
Residual joint tensile strength.
Was this helpful? ...  Itasca Software © 2022, Itasca  Updated: Mar 09, 2023 