Lecture 32 - Application of Numerical Modelling
Lecture 32 - Application of Numerical Modelling
Suspension effect: Bolts pre-tensioned to less than the design value may allow the roof to sag locally. This
causes a disturbance in the intact structure of the roof which may act as a zone where roof failure can
initiate.
Friction effect: It relies on the generation of a stress normal to bedding planes in the roof. For this
mechanism to be most effective, the distribution of normal stress on these planes must be as even as possible.
After the installation of pre-tensioned rock bolts, the roof often moves, and as a consequence, the load in the
rock bolts increases by an amount dependent upon their stiffness and the amount of movement.
The equation for pre-tension expressed as a function of the last four parameters above
is well known and widely used.
TORQUE-TENSION EQUATION
If a rock bolt is anchored at one end and a torque T applied to a nut bearing on a washer
(collar) at the other end, the tensile force F developed in the bolt is estimated by following
equation (Shigley, 1963):
𝑻
𝑭 = (𝟏 )
𝑲𝒅
where d is the bolt diameter K
(2)
= mean bolt diameter (average of major and minor diameters at the thread)
𝜽=𝒕𝒂𝒏
−𝟏
( 𝒍
𝝅 𝒅𝒎 )
= lead or pitch of thread
𝛂= 𝐩𝐫𝐨𝐣𝐞𝐜𝐭𝐞𝐝 𝐭𝐡𝐫𝐞𝐚𝐝𝐟𝐥𝐚𝐧𝐤 𝐚𝐧𝐠𝐥𝐞
= Coefficient of thread friction (nut/thread)
= Coefficient of collar friction (nut/washer)
THREAD FRICTION
Resistance to Applied Torque
• Friction between the thread surfaces in the nut and bolt provides an opposing moment to the applied torque which
generates a shear (torsional) stress in the bolt.
• As the total resistance to the applied torque comes from friction at both the thread and collar, it is possible to define
the effect of thread friction as a percentage of the total resistance.
• It is observed that if the coefficients of friction at the thread and collar are equal, then each provides approximately
equal resistance to the torque.
• Therefore an attempt to further reduce the torque resistance would be most effective if it was concentrated on
methods of reducing thread friction
Stress in the Bolt
The maximum shear stress a rock bolt developed during pre-tensioning is at the thread minor diameter and is given by:
( ) {( )}
′
𝟏𝟔 𝑲 𝒅𝑭 𝒅𝒎 𝒕𝒂𝒏 𝜽 + 𝝁𝒕 𝒔𝒆𝒄 𝜶
𝝉= 𝟑 𝑲 ′=
𝝅 ( 𝒅𝒊 ) 𝟐𝒅 𝟏 − 𝝁𝒕 𝒕𝒂𝒏 𝜽 𝒔𝒆𝒄 𝜶
is the minor diameter of the thread. The shear stress is a direct result of friction on the thread. During
tensioning of the bolt, this shear stress combines with the axial stress to give the total stress. Axial stress is
given by
= (4)
where is the core area (section at the minor diameter)
Based on the octahedral shear stress theory of failure, the maximum equivalent
normal stress is given by:
(5)
and must not exceed the yield stress of the bolt material. An assumption in the formulation of equations (4) and (5)
is that the bolt does not bend during installation and pre-tensioning.
To do this the objective of the test is to defeat the resin bond and to
characterise the bond failure in terms of :
BOND STRENGTH
and
SYSTEM STIFFNESS
BOND STRENGTH
BOND STRENGTH
APPLIED LOAD (kN)
and
SYSTEM STIFFNESS
SS
FN E
STIF
So why do we use
“SHORT ENCAPSULATION”
pull testing ?
SEPT-Principles and Objectives
If we do pull test a fully encapsulated bolt, installed using
a modern high strength/high stiffness polyester resin, we
would
BOLT
FACE PLATE
HYDRAULIC RAM
NUT
PULL BAR
DIAL GAUGE
HYDRAULIC PUMP WITH
FORCE INDICATOR
MONOPOD
SEPT: Equipment with Pull Bar
BOLT
SHIMS
BEARING PLATE
HYDRAULIC RAM
NUT
DIAL GAUGE
HYDRAULIC PUMP WITH
FORCE INDICATOR
MONOPOD
SEPT: Equipment with Pull Bar
HYDRAULIC JACK
• 30 tonne capacity, hollow cylinder
• centre hole diameter > bolt (35mm ID)
• single acting, spring return
• 50mm stroke
• effective area of cylinder 47cm^2
• closed height 165mm
• weight 10 kg
PUMP
• hand pump
• 700 bar working pressure
• oil capacity 1 litre
• length (typically) 560mm
• weight 5.6 kg
HOSE
• 3m reinforced hydraulic hose quick release couplers
SEPT: Equipment with Pull Bar
Steel shims
Bearing Plate
Tape bolt 300mm from bolt end to ensure required bond length, taping at least
300mm of the bolt
Ensure surface condition of bolt is free from rust, grease, paint or dirt or any
other surface contamination
300mm
SEPT: Bolt Preparation
SEPT: Bolt and Resin Capsule Diameter
3
Resin Bond Length
4
The purpose of calculating capsule length is to restrict the excess resin beyond
300mm
Prepare capsule using tie wraps to give correct amount of resin for 300mm
encapsulation
Required
length
1. Eye protection should be worn at all times. Twisting these capsules puts the contents under pressure
and pills could burst. Also, when cutting the centre section resin may burst out under pressure.
2. Wear gloves. Resin is a skin irritant.
3. If resin comes into contact with skin, wash off immediately.
SEPT: Resin capsule resizing
SEPT: Installing bolt
Rockbolt
Rockbolt
Rockbolt
Plate
Jack
Bolt Length Bolt rib Bolt core Hole Resin Resin Spare
900 O O O capsule O length thread
mm (mm) (mm) (mm) (mm) (mm) (mm)
130 185
140 210
BOND STRENGTH
By plotting
Applied load
Against
APPLIED LOAD (kN)
Bond displacement
we can determine
E SS
FFN
BOND STRENGTH
ST I
And
SYSTEM STIFFNESS
BOND DISPLACEMENT (mm)
SEPT: Analysis-Bond displacement
Bond Displacement
Measured Displacement
Applied Load
SEPT: Analysis-Bond displacement
Bond displacement (mm) = Measured displacement - Stretch in free length steel
We need to calculate the stretch in the steel for the bolt free length
PULL TEST RESULTS COLLIERY WHITE SECTION 35 - Split 218 Road L1-L2 DATE 14/3/02
130 185
140 210
Example:
100 kN
Imposed load (force) = 100kN (10 tonne)
APPLIED LOAD (kN)
.x…mm
BOND DISPLACEMENT (mm)
and
then,
LAnch = Y / GF
Obviously, a bolt with a larger Grip Factor will
have more available resistance over a longer length,
as shown in figure B.
In fact, the “Full Resistance Length” (L FR) of a
fully grouted bolt, which is the zone in which the
force available to resist rock movement is at least Effect of the Grip Factor on the resistance available from
equal to the yield strength of the bolt, is the total 10-ton roof bolts to react against roof loads. (A) Grip
bolt length L minus the length of the anchor: Factor = 0.5 tons/in; (B) Grip Factor = 1.0 ton/in.
LFR = L - LAnch
CAUSES OF POOR RESIN BOLT ANCHORAGE
1. Weak rock
It requires a longer grouted length to achieve the same anchorage as strong rock
Low grip factor
2. Poor installation
Defective grout can result from improper storage (too hot, too cold, too wet, or shelf life
exceeded), or (rarely) from manufacturing problems.
Improper mixing can occur if the proper spin time is not followed. Underspinning can
result in inadequate mixing, while overspinning can destroy the partially cured resin.
Improper mixing can also occur with long bolts where the top of the hole has less time to
mix before the bottom sets up. The temperature of the resin at the time of installation can
also affect the cure time.
Improper holes can be too long, too short, too large, or too smooth. The proper grout
cartridge must also be matched to the hole and the bolt being installed.
Finger gloving occurs when the plastic cartridge wrapper remains intact around the
hardened resin. It is more likely if the bolt is not rotated as it is inserted in the hole.
IMPROVING THE PERFORMANCE OF FULLY GROUTED BOLTS
Bridge element
Hydraulic
Hydraulic hand
jack pump
Circular
ring
Conical wedge
1.5m
Testing after 1st cycle of roof blasting Testing after 2nd cycle of roof blasting
Arrangement of in situ pull-out test
Arrangement of in situ pull-out test
Rock mass
Grout
Bolt
Confinement
Force
MN4021:MPD
Failure mechanism of rock bolt/cable bolt
1. Failure of system
The reinforcement system failed
• Shear or yield strength of individual components of the system
exceeded by the shear, compressive or tensile stresses generated with
bed separation
• Induced shear stress exceeds the shear strength between the
interfaces.
2. Failure of surrounding rock mass
Having weak/laminated strata/geological discontinuities
Considerable roof sagging with time under influence of induce
stresses.
Failure modes of rock bolts
Bond strength is mainly frictional and hence depends on the pressure build-up at
the interface which in turn depends on the dilational movement against the
confining grout or rock.
1. The first quadrant shows the variation of bond strength with axial displacement. It represents the
pullout test graph;
2. The second quadrant relates the confining pressure at the bolt-grout interface to the bond strength
using bond strength Equation;
3. The third quadrant shows the relationship between axial and lateral displacements. Since the apparent
dilation angle is not constant, the relation is non-linear and asymptotically approaches an ultimate
lateral displacement;
4. The dilation acts outward on the grout column and creates the interface pressure as illustrated by the
fourth quadrant. The straight lines show that the grout may split under the dilational pressure
In the fourth quadrant of the BSM, the dilatational behavior of grout is: (1) elastic; (2) fully split or (3) a transition
zone of partially split with an elastic portion.
Combining the above three equations, the displacement at the bolt–
grout interface can be expressed in terms of the internal pressure in
the form
Displacement at the bolt–grout interface for the totally split grout cylinder is:
Successive stages in the failure during a pull test
The mechanical interlock and frictional resistance is
related with the initial linear response during a pull test
From stage 2
• The failure mechanism is dependent on the radial
confining pressure.
• The stress drop may correspond to radial fracturing of
the grout annulus and/or shear failure through the grout
flutes.
• Cable displacement increases.
• Grout annulus can be pushed aside is determined by the
radial stiffness of the confining medium.
• When the radial stiffness is low, the favorable failure
mechanism is lateral displacement of the wedges;
• when it is higher, dilation is suppressed and failure is
more likely to occur by shear of the grout flutes and pull
out along a cylindrical frictional surface.
Application of Numerical Modelling
Simulation for estimation of angle of draw based on
excavation width and depth of cover
Coal seam
Panel width,120 m
Floor
70
Numerical simulation for caving angle with excavation width
1. Generate block
2. Assign constitutive model and properties to strata
3. Fix boundary
4. Initiate In-situ stress
5. Applied truncated load (for unmodelled depth of cover)
6. Solve (insitu model)
7. Excavate (model null range)
8. Solve (Quarter symmetry model)
9. Initiate Z velocity
10.Define stress
11.Define vertical displacement
12.Plot history
13.Solve
New Numerical modelling procedure- A Simple Example
;;floor=5
gen zone brick p0 (0, 0,-20) p1 (40,0,-20) p2 (0,40,-20) p3 (0,0,0) size (40,40,10) ratio (1,1,0.9) group floor
;;coal=4 ini xdisp=0
gen zone brick p0 (0,0,0) p1 (40,0,0) p2 (0,40,0) p3 (0,0,4) size (40,40,8) ratio (1,1,1) group coal ini ydisp=0
;;roof=5 ini zdisp=0
gen zone brick p0 (0,0,4) p1 (40,0,4) p2 (0,40,4) p3 (0,0,24) size (40,40,40) ratio (1,1,1) group roof model null range x(17.5,22.5) y(10,30) z(0,4)
;; floor=5 step 50
model mech elastic range group floor sel cable id 1 begin (20,20,4) end (20,20,5.8) nseg 9
prop dens 2500 bulk 4.67e9 shear 2.80e9 range group floor sel cable pretension 2.94e4
;;coal=4 sel cable prop emod=2.06e11 ytens=2.7e5 xcarea=3.8e-4 gr_coh=2.5e5
model mech elastic range group coal gr_k=2.0e9 gr_per=8.48e-2
prop den 1440 bulk 1.67e9 shear 1.00e9 range group coal step 2000
;;roof=5 save dev.sav
model mech elastic range group roof
prop dens 2500 bulk 4.67e9 shear 2.80e9 range group roof
;boundary condition
fix x range x -0.1 0.1
fix x range x 39.9 40.1
fix y range y -0.1 0.1
fix y range y 39.9 40.1
fix z range z -20.1 -19.9
p1
;******insitu stress
set gravity 9.81
p3
;sxx=syy=2.4+0.01h
;h= 100m
;szz=0.025h
ini sxx= -3.40e6 grad 0 0 0.010e6 range z -20 24
ini syy= -3.400e6 grad 0 0 0.010e6 range z -20 24
ini szz= -2.500e6 grad 0 0 0.025e6 range z -20 24 p2 p0
;truncated overburden load for 76m
apply szz= -1.9e6 range z 23.9 24.1
Solve
save insitu.sav
CALIBRATION OF MODEL
Design applied support using numerical modelling
Ribs (under stress)
Goaf edge
rock bolt supports in the models.
Line of extraction
It is examined in the middle portion of the panels when maximum diagonal length of line of extraction is attained
to develop maximum induced stress during depillaring.
A section of working height including immediate roof strata is cut along red dotted line as shown in Figure.
Application of numerical modelling in support design and strata control
during underground mining
Original gallery,
Two rows of 1.8m long bolts for RLH= 1.75m
RBGES
1.6m
Goaf edge, RLH=1-1.2m
Rib, 3 m
Roof of split
The values of RLH measured at goaf edge without RBGES is shown in Figure (left side).
Significant improvement in safety factor at the goaf edge is observed after application of RBGES (Right side).
Formulation for roof bolt-based goaf edge support during depillaring with stowing
where, rock load height at goaf edge in meter during depillaring with stowing, H is depth of cover in meter, W is
width of gallery in meter, RMR is Rock Mass Rating
Developed empirical formula for design of goaf edge support for semi-mechanised depillaring with stowing method of
strata control management in goaf is valid for depth of cover from 76-192m, rock mass rating from 42-58.5 and
gallery width from 4 to 5.5m.
It is also found that 1.8-2.0m bolt length for the RBGES is suitable for the studied range of variation in geo-mining
conditions.
Further, it would improve the mining cycle time and reduce the consumption of timber with the application of
RBGES during depillaring with stowing method of strata control management.
Assignment for Review
Thank you
Email: [email protected]
WhatsApp: 9470313115
MN4021:MPD