01 CRSP-3D Users Manual Entire Report
01 CRSP-3D Users Manual Entire Report
01 CRSP-3D Users Manual Entire Report
February 2012
FHWA-CFL/TD-12-007
4. Title and Subtitle
5. Report Date
February 2012
6. Performing Organization Code
J2009-05
7. Author(s)
DTFH68-07-D-00001
Final Report
April 2008 February 2012
14. Sponsoring Agency Code
HFTS-16.4
COTR: Matthew DeMarco, FHWA-CFLHD. Advisory Panel Members: Khamis Haramy and Roger
Surdahl, FHWA-CFLHD; Khalid Mohamed, FHWA-EFLHD; Rich Barrows, FHWA-FLH, and Barry Siel,
FHWA-OTS. This project was funded under the FHWA Federal Lands Highway Coordinated Technology
Implementation Program (CTIP).
16. Abstract
Rockfall is a natural result of weathering, erosion, and gravity on steep natural rock slopes and rock
cuts. When rocks become loosened from a steep slope, they may travel down the slope in a combination of
rolling, bouncing, sliding, launching, and free fall. When rockfall comes in contact with roads, structures or
railroads, it can cause destruction to property and becomes a safety issue. As traffic and building
developments increase in volume in mountainous areas, the need to protect people and structures from
rockfall increases in importance.
The original Colorado Rockfall Simulation Program (CRSP) was a 2D model developed for modeling
rocks rolling and bouncing down slopes. It has been used to predict the behavior of rockfall, determine the
need for rockfall mitigation, and aid in the design of mitigation. New numerical techniques that can be
applied to modeling of rockfall movement have become available since CRSP was first written in the mid
1980s. CRSP-3D uses the Discrete Element Method and provides a theoretically more representative
dynamic simulation of rock and slope interaction using the equations of motion than the approach employed
by previous versions of CRSP. In addition, CRSP-3D includes an option where the slope profile can be
loaded from scanned slope data, such as LIDAR. CRSP-3D estimates preferred rockfall paths, distribution
on the slope, rock bounce height, rock velocity, rock kinetic energy and rock rollout distances.
This users manual presents the theoretical basis for the computer program, discusses program
execution including inputs and outputs, and includes examples to demonstrate use of the program.
17. Key Words
Unclassified
Form DOT F 1700.7 (8-72)
Unclassified
22. Price
164
Reproduction of completed page authorized
LENGTH
in
ft
yd
mi
Multiply By
To Find
Symbol
inches
feet
yards
miles
25.4
0.305
0.914
1.61
millimeters
meters
meters
kilometers
mm
m
m
km
square inches
square feet
square yard
acres
square miles
645.2
0.093
0.836
0.405
2.59
square millimeters
square meters
square meters
hectares
square kilometers
mm2
m2
m2
ha
km2
milliliters
liters
cubic meters
cubic meters
mL
L
m3
m3
grams
kilograms
megagrams (or "metric ton")
g
kg
Mg (or "t")
Celsius
lux
candela/m2
lx
cd/m2
newtons
kilopascals
N
kPa
AREA
in2
ft2
yd2
ac
mi2
VOLUME
fl oz
gal
ft3
yd3
fluid ounces
gallons
cubic feet
cubic yards
29.57
3.785
0.028
0.765
NOTE: volumes greater than 1000 L shall be shown in m3
MASS
oz
ounces
28.35
lb
pounds
0.454
T
short tons (2000 lb)
0.907
TEMPERATURE (exact degrees)
F
Fahrenheit
5 (F-32)/9
or (F-32)/1.8
ILLUMINATION
fc
foot-candles
10.76
fl
foot-Lamberts
3.426
FORCE and PRESSURE or STRESS
lbf
poundforce
4.45
lbf/in2
poundforce per square inch
6.89
APPROXIMATE CONVERSIONS FROM SI UNITS
Symbol
Multiply By
To Find
Symbol
millimeters
meters
meters
kilometers
0.039
3.28
1.09
0.621
inches
feet
yards
miles
in
ft
yd
mi
square millimeters
square meters
square meters
hectares
square kilometers
0.0016
10.764
1.195
2.47
0.386
square inches
square feet
square yards
acres
square miles
in2
ft2
yd2
ac
mi2
milliliters
liters
cubic meters
cubic meters
0.034
0.264
35.314
1.307
fluid ounces
gallons
cubic feet
cubic yards
fl oz
gal
ft3
yd3
grams
kilograms
megagrams (or "metric ton")
0.035
2.202
1.103
ounces
pounds
short tons (2000 lb)
oz
lb
T
1.8C+32
Fahrenheit
0.0929
0.2919
foot-candles
foot-Lamberts
fc
fl
0.225
poundforce
lbf
LENGTH
mm
m
m
km
AREA
mm2
m2
m2
ha
km2
VOLUME
mL
L
m3
m3
MASS
g
kg
Mg (or "t")
Celsius
ILLUMINATION
lx
cd/m2
lux
candela/m2
newtons
kilopascals
2
0 145should be made to comply with
dfSection 4 of ASTM E380.
i h (Revised March
lbf/i 2003)
*SI is the symbol for the International System of Units. Appropriate rounding
______________________________________________________________________________
ii
TABLE OF CONTENTS
CHAPTER 1 BACKGROUND AND METHODOLOGY ......................................................1
INTRODUCTION ...................................................................................................................... 1
PURPOSE .................................................................................................................................. 1
LIMITATIONS .......................................................................................................................... 2
PREVIOUS CRSP VERSIONS ................................................................................................. 2
Development of CRSP............................................................................................................. 3
CRSP-3D .................................................................................................................................... 3
Discrete Element Method ........................................................................................................ 3
Contact Mechanics Between Falling Rock and a Slope .......................................................... 6
Stiffness and Damping (Slope Hardness) ................................................................................ 8
SYSTEM REQUIREMENTS .................................................................................................. 10
INSTALLATION ..................................................................................................................... 10
CD Installation....................................................................................................................... 10
Internet Installation ................................................................................................................ 11
FUTURE VERSIONS OF CRSP-3D ....................................................................................... 11
CHAPTER 2 - GENERAL PROGRAM EXECUTION ...........................................................13
STARTING CRSP-3D ............................................................................................................. 13
MAIN MENU........................................................................................................................... 14
File Menu Item ...................................................................................................................... 14
Units Menu Item .................................................................................................................... 15
Help Menu Item ..................................................................................................................... 15
TOOLBAR ............................................................................................................................... 16
TABS ........................................................................................................................................ 17
Project Tab ............................................................................................................................. 17
Slope Geometry Tab .............................................................................................................. 18
Slope Mechanical Properties Tab .......................................................................................... 27
Analysis Partitions Tab.......................................................................................................... 30
Falling Rock Parameters Tab ................................................................................................ 33
______________________________________________________________________________
iii
______________________________________________________________________________
v
LIST OF FIGURES
Figure 1. Schematic. Blocks in contact. .................................................................................... 5
Figure 2. Schematic. Visco-elastic contact model in normal direction. ................................... 7
Figure 3. Schematic. Elastic contact model in the tangential direction. ................................... 8
Figure 4. Screenshot. CRSP-3D start-up splash screen. ......................................................... 13
Figure 5. Screenshot. Main Menu. .......................................................................................... 14
Figure 6. Screenshot. File Menu Item. .................................................................................... 14
Figure 7. Screenshot. Units Menu........................................................................................... 15
Figure 8. Screenshot. Help Menu. .......................................................................................... 16
Figure 9. Screenshot. Example of Build version. ................................................................... 16
Figure 10. Screenshot. Toolbar. .............................................................................................. 16
Figure 11. Screenshot. Tabs. ................................................................................................... 17
Figure 12. Screenshot. Project Tab. ........................................................................................ 17
Figure 13. Screenshot. Slope Geometry Tab. ......................................................................... 18
Figure 14. Screenshot. Surface Points Tab. ............................................................................ 19
Figure 15. Screenshot. Select Excel Spreadsheet X, Y, Z Coordinate Data. ......................... 19
Figure 16. Screenshot. Highlight Paste Zone in CRSP-3D. ................................................... 19
Figure 17. Screenshot. Pasted X, Y, Z Coordinate Data......................................................... 20
Figure 18. Screenshot. Initialize scrollbars. ............................................................................ 20
Figure 19. Screenshot. Moving and deleting points on the screen. ........................................ 21
Figure 20. Screenshot. Graphic rotation and translation controls. .......................................... 21
Figure 21. Screenshot. Image Rotation. .................................................................................. 22
Figure 22. Screenshot. LIDAR Preprocessor. ......................................................................... 22
Figure 23. Diagram. The maximum surface roughness of 6 feet generated by spheres with
centers 12 feet apart. ................................................................................................... 24
Figure 24. Screenshot. Surface Resolution Selection. ............................................................ 25
Figure 25. Screenshot. Specifying Surface Resolution, Offsets and Spin value when
loading a LIDAR file. ................................................................................................. 25
Figure 26. Screenshot. Use of the Spin control to place the view in the X-Z plane. .............. 26
Figure 27. Screenshot. Select Refresh on the Slope Geometry Tab to create a slope
surface. ........................................................................................................................ 27
Figure 28. Screenshot. Select a slope material. ...................................................................... 27
Figure 29. Screenshot. Enter roughness and hardness values associated with a slope
material type................................................................................................................ 28
Figure 30. Screenshot. Identify the material type under the Slope Mechanical Properties
tab. ............................................................................................................................... 29
Figure 31. Screenshot. Place points around the perimeter of the identified material. ............ 29
Figure 32. Screenshot. Place points around the perimeter of the first material. ..................... 29
Figure 33. Screenshot. Boulder field zone (magenta) drawn on the screen using the mouse. 30
Figure 34. Screenshot. Enter partition name, direction, and location. .................................... 31
Figure 35. Screenshot. X and Z coordinates of the Partition. ................................................. 32
Figure 36. Screenshot. Dragging a Partition. .......................................................................... 32
Figure 37. Screenshot. The spherical elements and tubular connectors that make up
various rock shapes. .................................................................................................... 33
Figure 38. Screenshot. Falling Rock Parameters Input........................................................... 34
Figure 39. Screenshot. Magnified Simulation Tab Controls. ................................................ 35
______________________________________________________________________________
vi
Figure 79. Screenshot. Histogram showing kinetic energy frequency for simulation. ........... 68
Figure 80. Schematic. Slope profile identifying launch points and source areas (after
Brawner, 1994). .......................................................................................................... 70
Figure 81. Schematic. Volumetric Joint Count (JV). Measurements are made normal to the
direction of joint set. Line 1 intersects 5 joints, line 2 intersects 6 joints, and line 3
intersects 4 joints......................................................................................................... 73
Figure 82. Photograph. Moist Soft Soil. ................................................................................. 79
Figure 83. Photograph. Moist Soft Soil. ................................................................................. 80
Figure 84. Photograph. Soft Soil Indentation. ........................................................................ 81
Figure 85. Photograph. Dry Soft Soil. .................................................................................... 82
Figure 86. Photograph. Intermediate Soil. .............................................................................. 83
Figure 87. Photograph. Moderately Soft Soil. ........................................................................ 84
Figure 88. Photograph. Firm Soil. .......................................................................................... 85
Figure 89. Photograph. Firm Soil. .......................................................................................... 86
Figure 90. Photograph. Firm Soil. .......................................................................................... 87
Figure 91. Photograph. Gravel to Cobble Talus. .................................................................... 88
Figure 92. Photograph. Gravel to Cobble Talus. .................................................................... 89
Figure 93. Photograph. Gravel to Cobble Talus. .................................................................... 90
Figure 94. Photograph. Boulder Talus. ................................................................................... 91
Figure 95. Photograph. Boulder Talus. ................................................................................... 92
Figure 96. Photograph. Boulder Talus. ................................................................................... 93
Figure 97. Photograph. Hard Bedrock. ................................................................................... 94
Figure 98. Photograph. Hard Bedrock. ................................................................................... 95
______________________________________________________________________________
viii
LIST OF TABLES
______________________________________________________________________________
ix
two dimensional
three dimensional
Caltrans
CFLHD
CDOT
COR
COTR
CRSP
CSM
DEM
DR
EFLHD
FHWA
FLH
kip
kJ
lbm
lbf
lb
pounds mass
pounds force, pounds weight
pounds force, pounds weight
LIDAR
ODOT
OTS
psi
s
SI units
second
International System of units (e.g. m, N, kPa, etc.)
US units
WFLHD
X,Y,Z
______________________________________________________________________________
x
ACKNOWLEDGMENTS
The contributions and cooperation of the FHWA Advisory Panel; Khamis Haramy and Roger
Surdahl, FHWA-CFLHD; Khalid Mohamed, FHWA-EFLHD; Rich Barrows, FHWA-FLH, and ,
Barry Siel, FHWA-OTS are gratefully acknowledged. The efforts of Matthew DeMarco, in the
role of FHWA-CFLHDs Contracting Officers Technical Representative (COTR), Ty Ortiz and
C.K. Su of the Colorado Department of Transportation (CDOT) and Dr. Jerry Higgins of the
Colorado School of Mines (CSM), are appreciated.
______________________________________________________________________________
xi
(blank page)
______________________________________________________________________________
xii
INTRODUCTION
When designing rockfall protective structures such as ditches, berms, fences, and walls, it has
become common practice to first model the rockfall trajectory using computer simulations and
then assign barrier location(s) and design guidelines based on program output. Computer
modeling of rockfall trajectory is helpful in estimating how rockfall parameters such as bounce
height, velocity, kinetic energy, and rollout distance change along the slope length.
Current programs used for modeling rockfall trajectory require input of a 2D slope profile, slope
characteristics, and typical rockfall size. 2D modeling causes problems when simulating the
rotation and slope interaction of non-spherical rocks because obscure rock shapes have different
modes of rotation depending on velocity. For example, at high speeds a cylindrical shaped rock
will roll end-over-end, and act more like a sphere, while at slower speeds these rocks will tumble
and roll along the long axis. This may result in inaccurate estimates of rockfall properties,
especially rock rollout and bounce height. Additionally, current programs use a stochastic semiempirical approach which does not directly address the fundamental physics of falling rock
interaction with a slope face.
CRSP-3D uses the Discrete Element Method for dynamic model simulation using the equations
of motion, which provide a more accurate approximation of rock and slope interaction than the
methods employed by previous versions of CRSP. CRSP-3D uses 3D slope geometry to achieve
a more accurate approach in modeling the interaction and movement of rockfall on a slope
surface. The use of a 3D slope provides the ability to model different possible rockfall paths on
a section of slope, and has the capability to model the rotational movement of non-spherical
rocks.
This manual is intended to provide the user with the background required to apply CRSP-3D to
rockfall hazards studies.
PURPOSE
The CRSP-3D program is intended to provide an up-to-date analytical tool for evaluating
existing or potential rock fall events and to provide data to be used in the design of rockfall
fences, rockfall attenuators, catch ditches, catch berms and other rockfall protection structures.
The goal has been to improve the accuracy of rockfall simulations, in particular to improve
predictions of rock bounce heights and rock runout distance at the base of a slope.
Some specific features:
CRSP-3D uses a new approach which applies the Discrete Element Method to rockslope interaction and directly models interactions such as impact, rolling, launching,
sliding, and damping during a rockfall event.
______________________________________________________________________________
1
The slope hardness is represented using one numerical input value (hardness coefficient)
instead of two (tangential and normal coefficients of restitution) as in previous versions.
The slope profile can be de displayed either as apparent 2D or full 3D. The code is
actually 3D only; the 2D effect is obtained by looking at a section through the slope.
The user can build rocks of several shapes: spherical, cylindrical, discoidal, tetrahedral
and prismatic (blocky).
Analysis Partition parameters (e.g. equivalent to barrier height and position, if barrier
position is being modeled) can be input into the slope profile and kinetic energy of the
falling rock can be obtained.
LIMITATIONS
This program has been tested and is believed to be a reliable engineering tool; however, it should
be recognized that this is first version of an entirely new code. No responsibility is assumed by
the authors; Summit Peak Technologies, Yeh & Associates, the Federal Highway Administration
(FHWA), The Colorado Department of Transportation (CDOT), the Colorado School of Mines
(CSM) or any employees of the above for any errors, mistakes or misrepresentations that may
occur from any use of this program.
PREVIOUS CRSP VERSIONS
The original CRSP version 1.0 was written by Timothy J. Pfeiffer for a Master of Engineering
thesis in Geological Engineering at CSM in 1988. It was originally developed for CDOT for use
in conjunction with the construction of I-70 through Glenwood Canyon, Colorado. A reasonable
estimate of probable bounce height and velocity of rockfall events was needed in order to design
rockfall fences and alternative catchment ditches. CRSP was developed to simulate rockfall
events at a site based on estimates of slope profile, slope irregularities, slope material properties
and rockfall size. The project was supervised
by Dr. Jerry Higgins and the research was partially funded by CDOT. Later versions included:
Version 2.1 (1991) manual revised by Richard D. Andrew and program revised by Robert
Beck (CDOT).
Version 3.0 (1993) manual revised by Richard D. Andrew and Robert Beck and program
revised by Richard J. Schultz, Timothy J. Pfeiffer, Richard D. Andrew, and Robert B.
Beck (CDOT).
Version 3.0a (1995) program revised by Richard D. Andrew (CDOT).
Version 4.0 (2000) manual and program revised by Christopher L. Jones (CSM), Jerry D.
Higgins, (CSM) and Richard D. Andrew (CDOT).
______________________________________________________________________________
2
Development of CRSP
From version 1.0 through version 4.0, the CRSP algorithm did not change. Version 4.0 was
produced to be compatible with Windows 95 and Windows NT, it was reprogrammed using
Visual Basic, and the input/text editor section was modified to increase user-friendliness. The
input coefficients (normal coefficient of restitution and tangential coefficient of frictional
resistance) were recalibrated through the accumulation of rock rolling data provided by Caltrans
and the original CRSP calibration effort. The Caltrans data consisted of seven sets of data. Six
of the data sets are from Monterey County, CA (Duffy and Hoon, 1996; Duffy, 1996; Smith and
Duffy, 1990). The other set was collected in Switzerland (Duffy, 1992).
CRSP-3D
The CRSP-3D algorithm is a departure from previous versions of CRSP which were semiempirical in nature to a more fundamental physical computational model approach. The new
algorithm uses a combination of the general equations of motion and the Discrete Element
Method (DEM) to simulate rock velocity and the contact forces between rock and slope. DEM
assumes a linear elastic contact force-displacement relationship between two impact bodies. The
normal contact formulation is linear elastic with a viscous damping coefficient characterized by
two parameters; normal stiffness and viscosity. The tangential force depends on the coefficient
of restitution (COR) and the relative tangential velocity of the two contact objects. An
explanation of the program algorithm is included below.
In CRSP-3D, the definition of surface roughness has changed so that only one value is needed to
model the surface roughness for all rock sizes. The surface roughness is no longer dependent on
falling rock size. Additionally, CRSP-3D uses one input value, the hardness coefficient, to
model slope material properties. The hardness coefficient is a function of the coefficient of
restitution, which is a measure of the elasticity of the collision between a falling rock and the
slope, and a damping coefficient which models how much tangential resistance the rock is
subjected to upon sliding on the slope. A discussion of the equations behind CRSP-3D rockfall
modeling is included in the following section.
Discrete Element Method
The Discrete Element Method (DEM) is a numerical technique designed to solve problems in
applied mechanics that exhibit gross discontinuous material and geometrical behavior. DEM is
used to analyze multiple interacting rigid or deformable bodies undergoing large dynamic or
pseudo static, absolute or relative motion, governed by complex constitutive behavior.
DEM essentially is based on the numerical solution of the equation of motion and the principle
of dynamic relaxation. Kinematics equations are established for each discrete body. The
velocities, accelerations, and positions of the bodies are updated by calculating the contact forces
between them. Depending on different physical problems, DEM programs should at least
include the following three aspects:
______________________________________________________________________________
3
The following section discusses the discrete element method specifically related to CRSP-3D,
which discretizes the particles as three-dimensional spheres that contact each other at their
surfaces. Some general features of DEM are also included in this section.
Equation of Motion
Figure 1 shows two blocks I and II in contact. Their positions are defined by vectors R1 and R2.
The blocks have masses m1 and m2, linear velocity vectors v1 and v2, and angular velocity vectors
1 and 2 , respectively. The equation of motion for element i at discretized time step n is:
Ma ni Cv ni Pni x ni f ni
(1)
where:
x in , v in and a in are the position, velocity and acceleration vectors of the ith
element at the nth time step. The 3D case results in:
x x
i T
n
v x
i T
n
i
n
a x
i
n
i
n
, y in , z in , ni
, y in , z in , ni
i
n
, y in , z in , ni
(2)
where:
M mass matrix, (kg, lbm)
C damping matrix, (kg/s, lbm/s-gc), where gc is the gravitational constant.
Pni resultant contact force vector applied on the ith element at the nth step (N,
lbf)
i
fn the boundary force or body force vector applied on the ith element at the nth
step (N, lbf)
______________________________________________________________________________
4
II
Pt
I
Pn
v1
Fn2
v2
Fn1
R1
m1g
R2
m2g
x
Figure 1. Schematic. Blocks in contact.
Numerically solving equation (1) in the time domain gives accelerations, velocities,
displacements and resultant forces. The stress/strain relationship inside of the discrete assembly
is obtained by an averaging method. The average stress tensor of the volume V of the
representative of volume element (RVE) can be obtained by:
ij
1
V
mp
x F
p 1 c 1
c
i
c
j
(3)
where:
x ic is position vector at contact point c
ij
1
Fij F ji
2
(4)
______________________________________________________________________________
5
where
Fij is contact force.
where:
fn P
M / h C / 2 v
n 1 / 2
M / h C / 2
M / h C / 2
(8)
(9)
The explicit integration algorithm used in DEM analysis is simple and straightforward compared
to implicit schemes. It achieves fast and accurate solutions for dynamic problems such as
rockfall, impact, and granular flow problems. However, this algorithm is only conditionally
stable. The time step must be adequately small to maintain stability conditions. In CRSP-3D, the
time step is automatically set and adjusted by the program. When the algorithm is used to solve
static (or pseudo static) problems, dynamic relaxation procedures (DR) must be performed in
order to achieve rapid convergence.
Contact Mechanics Between Falling Rock and a Slope
The visco-elastic contact model, illustrated in Figure 2 as a spring and dashpot, is the most
popular constitutive contact model used in DEM analysis. Because of its simplicity, the
calculations are very efficient. The assumption of a linear elastic contact force-displacement
relationship between two impact bodies is a good approximation, and is still widely used in
engineering.
______________________________________________________________________________
6
The normal contact formulation is linear elastic with a viscous damper characterized by two
parameters: normal stiffness kn and viscosity C. The model works for both compression and
tension forces based on the relative distance between the two contact points. The normal force is
defined by:
1
k n Cv rn
N 2 n n
0 || max
max
(10)
where:
ks f | N |
1
sign(v rs ) k s t s 2
T
2
1
sign(v rs ) f N s 2 k s f | N |
(11)
where:
ks = shear contact stiffness
f
______________________________________________________________________________
7
Normal force
tangential unit vector
The direction of the tangential force is the reverse of the tangential relative velocity. The
magnitude of the tangential force is equal to the relative displacement of the falling rock and the
slope ( t ) times the shear contact stiffness, which is the first term of Equation 11. The second
term of Equation 11 refers to the case when the maximum static friction force is larger than the
Coulomb friction force (threshold or maximum friction at which movement between surfaces can
commence). This mechanical model is shown in Figure 3.
There may be three cases influencing the tangential force and rotation of the falling rock. When
the tangential force is less than the Coulomb frictional force, the rock is either in a static state or
rolling without sliding on the slope. When the tangential force is equal to or greater than the
Coulomb frictional force, the rock may slide without rolling, or slide while rolling on the slope.
Since the DEM numerical scheme is based on integrating the equation of motion, the velocity,
position, and rock rotation angle is updated based on the calculation of resultant forces on the
rock in each case.
Stiffness and Damping (Slope Hardness)
The key to successful modeling using DEM is proper selection of the stiffness and damping
coefficients. Theoretically, the damping coefficient C can be derived from material properties
such as the coefficient of restitution:
______________________________________________________________________________
8
k n mi m j
( mi m j )
1
C 2 ln 2
2
e ln(1 / e)
(12)
where:
mi and mj = the masses of particles i and j, respectively.
e = restitution coefficient of the material.
kn = the normal contact stiffness.
For the falling rock problem, mj is the mass of the slope element, which can be assumed to be
infinite. Equation 12 may be rewritten as:
knm
1
C 2 ln
2
2
e ln 1 / e
(13)
where:
m =
mi mj
mi + mj
Since the damping coefficient and the COR are directly related to the hardness coefficient, they
are used in calculating effects of rock and slope interaction. The damping and COR have an
inverse linear relationship. For example, a hard slope will have a relatively high COR and a low
damping value while a soft soil slope will have a low COR and a high damping value.
Slope and roadway materials have been grouped under general descriptions (Soft Clay, Hard
Clay, Sand, Gravel, Soft Bedrock, Hard Bedrock, Boulder Field, Talus, Concrete Road and
Asphalt Road). Because of the large range in hardness coefficient values observed in some
materials, the hardness coefficient overlap between different slopes under the same material
description was used to construct a table of recommended ranges of hardness coefficient values.
This table, shown in Material Descriptions in the CRSP-3D Help, is to be used as a starting point
for choosing hardness coefficient values.
During CRSP-3D rockfall simulations, the hardness coefficient values, for example, may need to
be calibrated to produce accurate output. This calibration process entails running rockfall
simulations of existing rockfalls (calibration sites) and comparing the output to observed rockfall
parameters at the existing sites. In this manner input parameters such as slope hardness,
roughness and boulder size at a required location can be estimated more accurately. An
automatic calibration mode capability is included in the advanced upgrade version of CRSP-3D.
______________________________________________________________________________
9
SYSTEM REQUIREMENTS
Minimum System Requirements:
Internet Installation
This public version of CRSP-3D can be found at
www.cflhd.gov/programs/techDevelopment/geotech/CRSP-3D.
Individually right click on "Entire Document - CRSP-3D Users Manual, Colorado Rockfall
Simulation Program and Self-Extracting Executable Program CRSP-3D, Colorado Rockfall
Simulation Program then save them onto your hard drive. Once copied, follow the CD
Installation directions shown above.
FUTURE VERSIONS OF CRSP-3D
Currently no additional revisions of CRSP-3D are under development though discussion is in
progress at FHWA with respect to the future of the program. It is possible that versions of the
code will be developed that incorporate changes and any corrections. Improvement to the
capabilities of the code may also be considered.
______________________________________________________________________________
11
(blank page)
______________________________________________________________________________
12
Thereafter, a Browse for Folder screen will open up requesting the user to Select a Project
Folder (an existing saved run) or Make a New Folder (for a new run). Once this has been
accomplished the splash screen will disappear for as long as 2 minutes before another screen
opens.
Note: though there is no physical indication on the screen that CRSP-3D is running at this point, it is processing
in the background and a startup screen (possibly minimized) will appear in a minute or so.
______________________________________________________________________________
13
MAIN MENU
The main menu in CRSP-3D (Figure 5) consists of three options: File, Units, and Help.
______________________________________________________________________________
14
The Units menu (Figure 7) allows the user to enter data in either customary US units (i.e. ft, lb,
etc.) or SI units (i.e. m, kN, etc.). The default setting for a new project is in customary US units.
If the user has input data in SI units, and accidentally enters those values while the units are set
to US units, then the user can simply select SI from the units menu. The input values will
remain the same, but the program will perform calculations in SI units.
Convert to US and Convert to SI are used when a user has entered data in one unit and then
wishes to see the equivalent problem in the other unit. The SI units designation in the top part
of the menu will automatically be checked, indicating the current unit in use.
Units->US: Use English US units, with no unit conversions.
Units->SI: Use metric SI units, with no unit conversions.
Units->Convert to US: Convert all current SI units to US units.
Units->Convert to SI: Convert all current US units to SI units.
Slight round-off errors may occur during conversions. Multiple conversions between units will
increase the magnitude round-off errors. The use of these menu items is discouraged for design
work. All US and SI units are displayed in the Report.
Help Menu Item
Selecting Help will open the Help files Table of Contents (Figure 8). The help function is
navigated by clicking on links in the Table of Contents listing, and clicking on Return to Table of
Contents and other links within each Help section. To hide Help, drag the frame all the way to
the right side of the screen:
______________________________________________________________________________
15
TOOLBAR
The toolbar is located to the right of the Menu, with icons corresponding to File menu items.
The toolbar (Figure 10) at the top of the screen gives you shortcut access to common menu
options that also appear under the File menu. From left to right:
New
Open
Save
Save As
Help
______________________________________________________________________________
16
TABS
Input parameters and output results in CRSP-3D are organized into the following tabs (Figure
11):
Project Tab
Slope Geometry Tab
Slope Mechanical Properties Tab
Analysis Partitions Tab
Falling Rock Parameters Tab
Simulation Tab
Output Tab
Project Tab
Basic project information such as Project Name, Project Number, Company Name, Location,
Map Coordinates, etc., is input under the Project Tab (Figure 12). The Project Name field will
be displayed over the slope image on the screen graphic as a title to identify the project. The
information entered into the Project Tab appears in the Report in order to identify the project,
data set and the analysis.
______________________________________________________________________________
17
______________________________________________________________________________
18
______________________________________________________________________________
19
Next, left-click the mouse and move points to the desired location (Figure 19).
Points may be deleted by right-clicking the point, or selecting the row of the grid and
pressing the Delete key.
______________________________________________________________________________
20
The slope graphic may be rotated and translated using the screen slider controls identified as
follows (Figures 20 and 21):
Tilt
Move Z
Zoom
Move Y
Spin
Move X
______________________________________________________________________________
21
Data Format:
CRSP-3D currently supports two ASCII text file formats for slope geometry input (file extension
.asc or .txt). These two file formats are current standards used by the Federal Highway
Administration (FHWA). CRSP-3D allows the user to convert the data from SI units to US units,
or from US units to SI units, after importing the data file.
The first file format contains a header section of any content and length (which is ignored),
followed by the line "X Y Z" (uppercase separated by spaces only), followed by data points in
the form "Easting Northing Elevation" (separated by spaces) as shown below.
SLIDE AREA 1
US State Plane 1983, Colorado Central 0502, NAD 83
US Survey Feet
XYZ
______________________________________________________________________________
22
The second file format contains a header section of any content and length (which is ignored),
followed by the line "X Y Z R G B" (uppercase separated by spaces only), followed by data
points in the form "Easting Northing Elevation Red Green Blue" (separated by spaces only), with
Red, Green, and Blue referring to color components from 0-255) as shown below:
SLIDE AREA 1
US State Plane 1983, Colorado Central 0502, NAD 83
US Survey Feet
XYZRGB
2938198.64 1679321.49 9114.45 0 0 4
2938199.59 1679321.55 9114.40 131 59 21
2938200.75 1679321.17 9114.25 5 0 4
A LIDAR text file with a large number of X Y Z records will run slowly and may be prone to
crashes. A solution is to "thin out" the raw LIDAR data to a fewer number of records using
software specifically designed for this. Another option to thin out the number of records is to
increase the surface resolution value (see Figure 24 and the discussion under Surface
Resolution). This will reduce the number of data records to a more manageable size. Input
reading time with a file size of five thousand records can be very slow (CRSP-3D was intended
for sizes between 500 and 2000 records). A thousand records seem to be a good compromise
between accuracy and execution time though the user should experiment to determine what filesize versus execution-time compromise is acceptable.
Once a LIDAR file is received it can initially be reviewed, edited and saved using a suitable text
editor (e.g. a freeware text editor like Crimson Editor works well). The edited file can then be
manipulated in spreadsheet program like Microsoft Excel which provides the capability to
perform data operations that an editor is not intended for. This pre-processed file can then be
saved as a text file (.txt, space delimited). Before proceeding to load the file into CRSP-3D it is
imperative that the file be reviewed in a text editor like Crimson Editor to make certain that the
file contains only blank space field delimiters. At this point, heading and X Y Z lines can be
added to the beginning of the file. Once saved this file is ready to be loaded into CRSP-3D.
Slope Surface Roughness:
The slope surface roughness is simulated using virtual spheres to create the surface. In order to
create the slope surface, CRSP-3D assumes a fixed 12 feet lateral spacing between sphere
centers regardless of the diameter of the spheres. These virtual spheres are constrained to vary in
radius between a minimum of 6 feet and a maximum of 60 feet which, respectively, correspond
to a surface roughness or amplitude of 6 feet over a base length of 12 feet (the maximum
roughness possible) and an approximately smooth surface. It is this variation in radius that
creates the variation in slope surface roughness. Roughness is defined as the perpendicular
______________________________________________________________________________
23
distance between the intersection of two adjacent virtual spheres and a tangent line to the two
spheres as shown below in Figure 23.
6 ft
12 ft
Figure 23. Diagram. The maximum surface roughness of 6 feet generated by spheres with centers 12 feet
apart.
Surface Resolution:
CRSP-3D uses a fixed 12 feet center-to-center distance for modeling a slope surface when
surface coordinates are cut and paste from a spreadsheet, entered manually or inserted
graphically. This is also the case for a LIDAR file regardless of the spacing of data points
contained in the file. This 12 feet spacing controls the roughness of the generated slope surface
and is referred to as the resolution when processing a LIDAR file. Therefore, using a surface
resolution of 12 ft when importing a LIDAR file will create a slope surface with the same
roughness as that generated by CRSP-3D when creating a slope surface created using manually
inserted coordinates. Loading a LIDAR file and specifying a resolution value less than 12 ft (a
higher resolution) will only degrade performance with no other improvements. Using a surface
resolution value greater than 12 ft (a lower resolution) will reduce rockfall simulation time and
will reduce the uploaded file size (the number of LIDAR records will be "thinned out"), but will
cause CRSP-3D to interpolate between points to achieve a virtual 12 ft surface resolution
resulting in a smoothing of the model slope surface.
CRSP-3D will maintain a constant surface point spacing (points that define the surface of the
slope being analyzed) of 12 feet - either by extrapolating new point data points from data
originally provided at a greater spacing than 12 feet or by thinning the number of data points
provided at a closer spacing than 12 feet. For example, on a 1000 foot by 1000 foot slope there
will be approximately 7000 surface points spaced 12 feet apart. CRSP-3D will reduce the
number of surface data points to 7000 for large input data files or increase the number of surface
data points to 7000 for small input data files.
By increasing the Surface Resolution value to greater than 12 (see Figure 24), one is increasing
the spacing between points that define the surface and therefore decreasing the surface
resolution. For example, by increasing the surface resolution value to 48, the total number of
surface points defining a 1000 foot by 1000 foot slope will reduce to about 450.
______________________________________________________________________________
24
Users should use their judgment to define the surface resolution when processing LIDAR files.
For example, for smooth planar slopes, selecting a surface resolution of 48 feet will have little
impact on simulated slope topography compared to selecting a 12 feet surface resolution.
Conversely, on very rough slopes, a 12 feet surface resolution will show more representative
slope topography compared to choosing a 48 feet surface resolution. Surface resolution is set at
the time of upload, and cannot be adjusted after the file has been uploaded in CRSP-3D.
Note: Use multiples of 12 when changing surface resolution values e.g. 12, 24, 36, 48, etc. (shown in order of
decreasing surface resolution).
Figure 25. Screenshot. Specifying Surface Resolution, Offsets and Spin value when loading a LIDAR file.
Offsets:
The leading digits that are common to each of the coordinate values can be eliminated for
purposes of analysis of a local feature like a slope. These are referred to as offsets and can be
removed in two different ways; (1) manually (Figure 25) as described in the following paragraph
or (2) via the Autocalc Transform option (see Figure 24), which will perform the appropriate
offset calculations automatically without the need to select offsets.
______________________________________________________________________________
25
Enter X, Y and Z offsets to reduce the number of digits that need to be stored in the geometry
file. When the file is loaded, these values will be subtracted from the X (Easting), Y (Northing),
and Z (Elevation) coordinates in the file, thereby defining position in local coordinates.
Note: The screen graphics rendering does not support coordinate values with more than 5 digits, so removing offsets may be
necessary in order to reduce the number of coordinate digits to no more than 5 to avoid potential display problems.
Spin:
A rotation angle may be entered for Spin to rotate the image about the Z-axis in order to align the
fall line of slope along the X-Axis (Figure 26). Increments of 90 degrees should be used to
achieve this.
Load File:
Clicking the Load File button imports the LIDAR data into the Surface Points Tab (Figure 25).
Refresh:
Select Refresh (Figure 27) on the Slope Geometry Tab to create or to update the slope surface if
slope coordinate data has been newly input or modified.
Figure 26. Screenshot. Use of the Spin control to place the view in the X-Z plane.
______________________________________________________________________________
26
Figure 27. Screenshot. Select Refresh on the Slope Geometry Tab to create a slope surface.
A slope material may be selected by clicking the down-arrow icon in the ID cell (Figure
28):
Note: There are only nine colors that can be used to represent zones and correspondingly there are only nine
material IDs that can be used on a single slope. The first nine colors that appear in the drop down box are
available - the remaining gradient color options that show are disabled.
______________________________________________________________________________
27
Enter Roughness values or ranges (Figure 29). Roughness values between 0.01 (smooth)
and 6 (rough) can be entered. Some judgment as to the selection of this value is
recommended since the numbers represent the amplitude of the roughness measured in
feet with respect to a flat surface over a reference distance of 12 feet. See the discussion
on slope Surface Roughness presented in previously.
Enter Hardness values or Ranges (Figure 29). Hardness varies between 0.01 and 1 and
represents the amount of energy that is transmitted between bodies during impact.
Hardness also correlates to the degree of damping i.e. a high hardness will have low
damping ratio and vice versa since the damping coefficient and the coefficient of
restitution (COR) are directly related to the hardness coefficient as discussed previously
in Chapter 1.
Note: The surface roughness, when entered into the program as a single number, such as 2.0, will result in a
constant surface roughness being modeled in the program. If a range of surface roughness is input, such as 0.1 2.0, the program will randomly choose roughness values between 0.1 and 2.0 for each element location on the
slope. The user can also use a range for the hardness coefficient, which is entered the same way as the roughness.
A range of hardness values can be used to simulate dispersed trees and large boulders on the slope.
Figure 29. Screenshot. Enter roughness and hardness values associated with a slope material type.
The mouse location is shown along the bottom of the program window. Several coordinates are
listed:
Mouse: (mouse coordinates)
Surface: (coordinates of mouse projected on the slope surface). Use surface coordinates
for placing points on the slope surface.
Center: (coordinates of viewing center)
Span: (horizontal width of view, similar to zoom)
______________________________________________________________________________
28
Note: When placing points near the edges or top of the slope using a mouse, the points are sometimes set far
behind the slope surface. To alleviate this problem, rotate the slope image to look straight down on the slope
Identify
theanmaterial
type
(similar to
aerial view)
and then place the points.
Select the material from the drop down menu (Figure 30).
Figure 30. Screenshot. Identify the material type under the Slope Mechanical Properties tab.
Click on the slope in regions where this material may be found (Figure 31).
Figure 31. Screenshot. Place points around the perimeter of the identified material.
Select the next material and define its extent (Figure 32). Click on the appropriate
material and then click the mouse on the slope at points corresponding to locations of this
material on the slope.
Figure 32. Screenshot. Place points around the perimeter of the first material.
______________________________________________________________________________
29
To draw a zone on the display, simply use the left mouse button to create a series of
points along the perimeter of the zone. Points can be deleted by placing the cursor on the
point to be deleted and clicking the right mouse button. Points can be dragged by holding
the left mouse button on the point and dragging it to the desired location. The zone
within two perimeters or between a perimeter and the boundary of the model (depending
on which is applicable) will be assigned the input roughness and hardness values.
After entering a set of mechanical properties data click Build to update the slope (Figure
33).
Figure 33. Screenshot. Boulder field zone (magenta) drawn on the screen using the mouse.
______________________________________________________________________________
30
Note: The partition does not act as a physical barrier. Also note that the number of partitions generally should
not exceed four because this may increase execution time significantly and create high graphics memory demand
which may lead to program instability.
Partitions can be located where observed rockfall bounce heights (tree scars) or rollouts were
measured or where rock fall experimental data are available. This information can then be used
to calibrate the model to the observed rockfall parameters (see Calibration later in this
document).
Select the Analysis Partition tab to enter locations of analysis partitions on the slope.
______________________________________________________________________________
31
231ft
70ft
Alternatively, clicking the right mouse button on the desired location on the slope will
add points representing the position of a partition.
CRSP-3D will calculate Kinetic Energy, Velocity and Rock Trajectory Height (bounce
height) for falling rocks at the location of each partition (rocks must actually pass the
partition for the energy and height calculations to be performed). Partitions can be
deleted by clicking the right mouse button on the base of the Partition to be deleted.
They can be dragged up and down the slope by holding the left mouse button on the base
of the Partition and dragging it to the desired location (see Figure 36).
______________________________________________________________________________
32
the zone where rocks are released (the program will randomly choose rockfall starting
locations within the material zone designated as the release zone),
the number of rocks to release,
the unit weight of the rock,
the shape of the rocks defined by the number of spherical elements combined to create
each rock shape. See Figure 37. Number of spherical elements to create a falling rock
shape: one spherical element = sphere, two spherical elements = cylinder, three spherical
elements = triangular disc, four spherical elements = tetrahedron, eight spherical
elements = cube), and
the size of the falling rock defined in terms the characteristic dimension for each shape. In
the case of:
(1) a single element sphere, size is the diameter;
(2) a two-element cylinder, size is the length of the cylinder;
(3) a three element triangular disc, size is the triangular side length;
(4) a four element tetrahedron, size is the height; and
(5) an eight element cube, size is the length of the side of the cube.
Figure 37. Screenshot. The spherical elements and tubular connectors that make up various rock shapes.
Rock size can be entered as a single number if the simulated rock(s) will be equidimensional. If
a range of rock sizes (e.g. 1- 4.5 feet) are to be simulated, a range of rock sizes are randomly
selected within the designated range. This rock size corresponds to the characteristic dimension
of the rock shape specified.
______________________________________________________________________________
33
Rock size and shape can have an effect on output parameters such as runout distance, bounce
height, velocity, and kinetic energy. CRSP-3D does not output the simulated rock size or shape
in the Stats tab (see Figure 41); however, this information is displayed in the Report if the user
chooses to generate it. CRSP-3D simulates rocks and outputs the results in the Stats tab in the
order that they are entered in the Falling Rock Parameters tab. The user can model several
different rock shapes and sizes in one rockfall simulation and observe the difference in output by
using the numbers listed on the left hand column of the Stats tab.
Simulation Tab
The Simulation tab performs the actual execution of the simulation. The rockfall simulation
model is run from the Simulation tab. See Figures 39 and 40.
______________________________________________________________________________
34
Output Tab
The Output tab contains a Stats and a Report sub-tab.
______________________________________________________________________________
35
Stats Sub-Tab
Output of the rockfall simulation can be viewed under the Stats sub-tab (Figure 41). Output can
only be viewed at barrier Partitions, so at least one partition must be included to view output.
The Stats sub-tab displays various rockfall properties for each rock at each Analysis
Partition.
The numbers down the left-hand column of the Stats sub-tab represent the simulation
sequence numbers of each rock rolled in the analysis. CRSP-3D first simulates the
rockfall identified by data input in the first row under the Falling Rock Parameters tab and
thereafter simulates rockfalls sequentially as identified by the data input in each
succeeding row.
Figure 41. Screenshot. Rock Trajectory Height simulation output for each rock passing Analysis
Partition BP-1.
Clicking on a cell in the Stats grid displays a histogram of the selected rockfall property at
the selected Analysis Partition (Figure 42).
Figure 42. Screenshot. Rock Trajectory Height Frequency output for Partition BP-1 as presented in the
Report option.
Checkboxes allow various image views and rockfall properties to be specified to show in
the Stats and Report (Figure 43).
Note: It is advisable to check only those boxes of direct interest (generally 5), otherwise the volume of data to
display can overwhelm the amount of graphics memory available leading to program instability.
______________________________________________________________________________
36
Figure 43. Screenshot. Checkboxes that select various image views in the simulation report.
Report Tab
The Report tab provides the option to automatically generate a report of a simulation run
complete with graphics. This report can be cut-and-paste into another document using the
options provided.
Note: The model display should be rotated and translated to the desired viewing position before generating the
Report under the Report Sub-Tab.
Note: The report may take 5 or more minutes to process depending on computer hardware characteristics and on
occasion may appear to have stop working. It is important not to click on any other tabs or buttons during this
time. The Report tab will be automatically selected when the report completes.
______________________________________________________________________________
37
To Print the report use the print buttons at the top of the Report tab (Figure 45).
To cut and paste the entire report into another document, click in the Report body and
then:
(1) to select all
type Ctrl-A,
(2) to copy
type Ctrl-C, and
(3) to paste into a MS-Word document
type Ctrl-V.
______________________________________________________________________________
38
DELETING FILES
To delete unwanted CRSP-3D rockfall simulation files, go to the drive and folder location where
the files are saved and delete them. CRSP-3D does not have an option to delete files within the
program.
CALIBRATION
CRSP-3D rockfall simulations may need to be calibrated to produce acceptable output. The
calibration process entails running rockfall simulations of existing monitored/measured rockfalls
(calibration sites) and comparing the simulation output parameters to observed rockfall
parameters obtained at these sites. In this manner, input parameters, such as slope hardness and
roughness, at a required location can be adjusted to simulate a rockfall more accurately.
The hardness coefficient and roughness values can be calibrated by entering candidate roughness
and hardness values and comparing the output of the simulation to the observed bounce height
or runout distance at a specific location on the slope (with the Analysis Partition positioned as
appropriate for simulation). The input values are adjusted until the resulting simulated bounce
height and runout values correspond to actual measured values.
Rock shape and size variables can be added to bounce height and runout distance to improve the
calibration, if desired. Rock size and shape can have an effect on output parameters such as
runout distance, bounce height, velocity, and kinetic energy. CRSP-3D does not output the
simulated rock size or shape in the Stats tab; this information is displayed in the Report.
CRSP-3D simulates rocks and outputs the results in the Stats tab in the order that they are
entered in the Falling Rock Parameters tab. As indicated previously the user can model several
different rock shapes and sizes in one rockfall simulation and observe the difference in output by
using the numbers listed on the left hand column of the Stats tab
Actual velocity and kinetic energy values may require a rock rolling experiment, possibly with
the use of high speed cameras and reference lines. Bounce height and runout distance can be
measured on a natural slope with rockfall accumulation.
______________________________________________________________________________
39
(blank page)
______________________________________________________________________________
40
TUTORIAL EXAMPLE
This example can be loaded at program startup by selecting the Tutorial folder when the Select
Project Folder window opens. A complete copy of the CRSP-3D report output for this example
can be found in Appendix B.
Creating a New Project
To create a new, blank project:
You can browse to the location where you would like your new project to be stored.
Select Make New Folder to make a new project. A new project folder will appear.
______________________________________________________________________________
41
A new window opens up at a blank Project tab which can then be filled with the
appropriate identifying project information (Figure 48).
Slope Geometry
In this example, slope geometry data are cut from an Excel spreadsheet and pasted directly into
Surface Points, as shown in Figure 49.
Note: That the graphic display routine is not always stable with coordinate values expressed in more than 5 digits.
This should be confirmed by trial and error.
______________________________________________________________________________
42
Figure 49. Screenshot. Cut and paste data from Excel spreadsheet.
Once the geometry data has been pasted, click Refresh to create the solid geometry. A
graphic model will be generated, as shown in Figure 49.
______________________________________________________________________________
43
Figure 50. Screenshot. Graphic rendition of slope geometry.
Table 1. Procedures for editing slope points.
Function
Procedure
Add point
Delete Point
Place the mouse cursor on the point to be deleted and right click the mouse
Move Point
material ID, or clicking a grid cell in the desired row. There is a drop-down menu of material
types (which can easily be typed over with an alternative user-defined description, if desired), as
well as an adjacent drop-down color palette that associates a material type with a colored area on
the screen graphic.
In this example (Figure 51) two material types are being modeled, Boulder Field and Talus:
Boulder Field has arbitrarily been assigned a hardness of 0.7-0.9 and the Talus has been
assigned a hardness of 0.3-0.6. Talus has arbitrarily been assigned a roughness of 0.3-0.5
and the Talus has been assigned a roughness of 0.3-0.6 (see Figure 52).
______________________________________________________________________________
45
Note: Hardness and Roughness values can be manually calibrated from known rockfall information in order to
obtain realistic rockfall response data, as previously described in this manual.
Note: If a range of values are input (e.g., 0.1-2.0 for surface roughness), the program will randomly choose
roughness values between 0 and 2.0 for each element location on the slope within the zone occupied by that
material.
The area within the zone will be assigned the input roughness and hardness values.
______________________________________________________________________________
46
Press Build to update the graphic to show the distribution of the material types on the
slope (Figure 53).
Analysis Partitions
Select the Analysis Partition tab to enter locations of analysis partitions on the slope. The user
enters a partition name under ID (Figure 54). Thereafter (Figure 55):
Enter the location of the analysis partition by either clicking the left mouse button on the
profile or manually entering the coordinates. The location of the partition can be entered
by manually entering a position across the slope face as represented by a fixed X
coordinate.
Note: Partitions can be deleted by clicking the right mouse button on the base of the partition to be deleted.
Partitions can be dragged up and down the slope by holding the left mouse button on the base of the partition and
dragging it to the desired location.
______________________________________________________________________________
47
or:
Enter a fixed height Z of the top of the barrier (the elevation) coordinate above the slope.
Alternatively, clicking the left mouse button on the desired location on the slope will add a
point representing the position of a partition. Give the partition a name (such as BP-1) and
a Z coordinate height of the top of the partition (such as 70 ft). Add another partition
called BP-2 and select a horizontal X coordinate position on the slope face of 231 ft.
Note: CRSP-3D requires only the coordinate of the partition to be entered (either X for location on slope face or
Z for the fixed elevation of the top of the barrier). CRSP-3D will calculate rockfall kinetic energy, velocity and
bounce height for rocks as they fall past the partition locations.
______________________________________________________________________________
48
______________________________________________________________________________
49
Simulation
Click Run Simulation to execute (Figure 58). A disclaimer precedes each run (Figure 57).
______________________________________________________________________________
50
Figure 58. Screenshot. Rock trajectory height traces generated by the rockfall simulation.
______________________________________________________________________________
51
3D LIDAR EXAMPLE
The following example is a 3D simulation of a rockfall using LIDAR data. The LIDAR
information is read into CRSP-3D using built-in program features; however, some preprocessing of the LIDAR data is required prior to downloading the data to the program. The
steps involved in creating a new project folder are not repeated here since they are common to all
simulations and have been presented previously.
This example can be loaded at program startup by selecting the Canyon LIDAR 3D Sphere folder
when the Select Project Folder window opens. A complete copy of the CRSP-3D report output
for this example can be found in Appendix B.
Slope Geometry
LIDAR data requirements are specific and should be followed carefully. The LIDAR data need
to be formatted as described in Chapter 2, LIDAR Preprocessor Sub-Tab.
The following provides a step-by-step procedure that was used to manipulate a LIDAR file to be
uploaded by CRSP-3D. Not all LIDAR files may be created the same, so the procedure may
differ on occasion. The LIDAR file is initially in .txt format, but had commas used as delimiters
(see Figure 59).
1. To remove the comma delimiters, upload the .txt file into Microsoft Excel by choosing the
Data tab and then choosing From Text (Figure 60). This will upload the .txt file into
Microsoft Excel with the ability to remove the comma delimiters.
______________________________________________________________________________
52
2. Once the LIDAR points are showing on the Excel screen, make sure two significant digits
or less are used for the LIDAR points. Also, add X Y Z on the top of each column, and
add a title for project tracking (Figure 61). Save the file as a .txt file.
3. Now use the text editor (in this example Crimson Editor is used) to upload the .txt file and
view the delimiters in the file (Figure 62).
______________________________________________________________________________
53
4. Make sure all delimiters are set to spaces (Figure 63). Save the file as a .txt file. The
saved file will now be ready for upload to CRSP-3D.
______________________________________________________________________________
54
By checking the AutoCalc Transform box (Figure 64), the program automatically calculates the
X, Y and Z offsets, which is usually the best approach to follow, though the user may choose to
leave it unchecked and provide the offsets manually.
Figure 64. Screenshot. AutoCalc Transform box and Load File button.
To load the file, the Load File button is clicked and the program will request the location and the
name of the file to load (Figure 65).
______________________________________________________________________________
55
Once the file loads successfully, the calculated offsets will be displayed (if the AutoCalc
Transform box was checked), the X, Y, Z coordinates of the surface will be displayed, and a
graphic rendition of the slope geometry will appear (see Figure 66).
Note: This may take several minutes to appear depending on the size of the file being processed.
______________________________________________________________________________
56
Figure 67. Screenshot. Defining materials zone and adding mechanical properties information.
Analysis Partitions
For this example four locations (in blue) on the slope were analyzed as indicated in Figure 68.
These analytical partitions were defined using the X-coordinate option to locate them across the
slope face.
______________________________________________________________________________
57
______________________________________________________________________________
58
Figure 69. Screenshot. Falling rock parameters (release zone shown in green).
Simulation
Simulation of twenty round rocks, each of 15 feet in diameter, falling from random locations
within the identified green release zone was performed. The combined velocity traces shwing
both velocities and rockfall path for all of the 20 rockfalls are presented in Figure 70.
______________________________________________________________________________
59
Figure 70. Screenshot. Velocity traces of 15 ft diameter round rocks falling down slope.
Output
Statistics
Statistical information gathered during the simulation is presented under the Stats tab.
Histograms of all the values of a parameter of interest in a particular analysis partition can be
displayed by highlighting (see in blue) the specific parameter as indicated out in Figure 71.
______________________________________________________________________________
60
Figure 71. Screenshot. Frequency information for all velocities calculated at Analysis Partition 1.
______________________________________________________________________________
61
2D EXAMPLE
The following example is a 2D simulation of a rockfall using manually input slope geometry
data. The steps involved in creating a new project folder are not repeated here since they are
common to all simulations and have been presented previously.
This example can be loaded at program startup by selecting the Canyon 2D Sphere folder when
the Select Project Folder window opens. A complete copy of the CRSP-3D report output for this
example can be found in Appendix B.
Slope Geometry
In this example the slope geometry being considered is two-dimensional. Since CRSP-3D is
inherently three-dimensional, two-dimensional analysis can be accommodated by rotating the
profile of the slope into the X-Z plane.
In a CRSP-3D graphic slope representation, the X-axis (slope length projected onto the
horizontal plane) represents the horizontal component of the slope cross-section and the Z-axis
(slope height) represents the vertical component of the slope cross-section. In 2D the Y-axis
(slope face width) is normal to the screen or printout. A 2D geometry, therefore is entered in X,
Z coordinates. In order to maintain compatibility with the 3D program requirements, a constant
Y coordinate is required. A fixed slope face width is calculated if a constant Y coordinate is
entered as part of each X Y Z set of coordinates (Figure 72).
In the current example, slope geometry data was input manually using slope coordinates from an
X, Z coordinate table.
______________________________________________________________________________
62
Figure 72. Screenshot. Rotated view of slope geometry with slope coordinates showing how a virtual 2D
geometry is created from an inherently 3D model.
Coordinates are entered for the crest of the slope i.e. the highest elevation (largest Z
coordinate) and the corresponding X coordinate (lowest X coordinate).
The next pair of X, Z coordinates corresponds to the lower edge of the uppermost material
zone.
The following pair of X, Z coordinates represents the next lowermost boundary and so on
for each consecutive material zone down the slope. This process is repeated until the
bottom of the slope is reached. See the generated profile in Figure 74.
______________________________________________________________________________
63
Figure 73. Screenshot. View showing slope mechanical properties and graphic rendering.
Figure 74. Plot. Material profile down slope (ft). Each point represents a change in slope material.
______________________________________________________________________________
64
Analysis Partitions
The selected partitions are identified in Figure 75 (blue lines)
______________________________________________________________________________
65
Simulation
Kinetic Energy simulation results are shown graphically in Figure 77.
Figure 77. Screenshot. Rotated simulation view showing rockfall kinetic energy.
______________________________________________________________________________
66
______________________________________________________________________________
67
Statistics
Figure 79 shows a sample of the statistical output from CRSP-3D consisting of a histogram
showing the distribution of kinetic energies for each rock that was simulated falling down the
slope.
Figure 79. Screenshot. Histogram showing kinetic energy frequency for simulation.
______________________________________________________________________________
68
______________________________________________________________________________
69
Figure 80. Schematic. Slope profile identifying launch points and source areas (after Brawner, 1994).
As already indicated,a critical part of rock fall modeling is the generation of a slope profile
defined topographically by X and Y coordinates in the case of 2D profiles or by X, Y and Z
coordinates in the case of 3D profiles. Generally a slope can be divided into three areas, each
displaying differing behavior during a rockfall event. Ideally, these zones will be separate, but
on many slopes, boundaries of the zones can overlap. The first zone is the rockfall source area
typically located near the top of the slope composed of bedrock outcrops, cut rock slopes, or
coarse surficial deposits. The middle zone is referred to as the travel path and represents the
course that the rock will following after separating from the source zone. The runout zone is the
area of deceleration and stopping for rockfall and can be located on slope or at the base of the
slope.
Since CRSP-3D only performs 3D analyses, X, Y and Z coordinates are mandatory though the Y
coordinate (laterally across the face of the slope) can be set to a fixed value, equal to roughly half
to a third the height of the slope, to allow what is effectively a 2D analysis. To perform what is
effectively a 2D analysis, only a 2D (X-Z) geometry profile needs to be surveyed/measured
since the third dimension (Y) is provided as a user-estimated constant value. The slope may then
be viewed in 2D so that the slope cross-section has only X and Z Cartesian coordinates visible.
Zone Boundary Selection
Slopes are divided into different material zones that identify areas of differing properties.
Properties that are generally of significance to rockfall sites are the relative hardness of both the
source rocks and various slope zones as well as the surface roughnesses of each of these zones.
Surface materials along the rockfall profile should be identified in the field in order to define
property characteristics required for modeling. Some indication of the relative hardness of the
surface and description of the surface material is required including amount of exposure as a
percentage the slope surface being surveyed.
In the field, the location zone boundaries can effectively be established by measuring the slope
angle and length and width of a zone. Various methods can be used to achieve this. CRSP-3D,
______________________________________________________________________________
70
however, has limits to the number of zones that can be employed during a single simulation. The
number of possible zones is controlled by the number of unique material ID values representing
these zones that are available in CRSP-3D.
Rock Size and Shape
Rock size and shape influence the degree to which slope surface roughness affects rockfall
behavior. Rockfall data show that spherical or near spherical rocks move faster downslope than
other irregular shapes. Large rocks tend to have longer runout because of their greater kinetic
energy.
Rock shapes can be generalized as spherical, tabular, or discoidal etc. to facilitate analysis.
These shapes can have a large effect on the rock trajectory during a rockfall. Spheres generally
represent a worst case scenario since they comprise the greatest mass for a given size and tend to
accelerate faster and decelerate slower than angular rocks. Elongated shapes, such as cylinders,
normally take a longer slope distance to reach high velocities and have shorter runout distances
compared to spherical rocks. At high velocities, elongated rocks act much like spherical rocks
because they are rotating end-over-end so the rock contacts the slope at just two points. During
acceleration and deceleration, angular rocks tend to tumble and fall while spherical rocks roll
without needing to overcome any undulation in their rock surface, resulting in a longer runout
distance.
Any disintegration of rocks as they fall tends to dissipate energy, so disintegration can influence
modeling results. Furthermore rock sizes that accumulate at the toe of rockfall zones in the
runout area may not be representative of the initial rock sizes at launch or early in their
trajectories.
Although falling rocks may break up during descent, the largest rocks found at the base of the
rockfall path that can be identified as having fallen from the source area may be considered a
reasonable representation of a maximum rock size case for CRSP-3D modeling. If no rocks are
available at the base of the path, then rock size can be established by measuring joint spacing in
the source area and estimating a probable size for falling rock.
CRSP-3D allows spherical, cylindrical, discoidal, tetrahedral and cube-shaped rocks. The
investigator should choose the appropriate rock shape for CRSP-3D rockfall simulation. To
represent the probable worst case on a slope, a spherical rock shape can be used.
Block Size
The discontinuity orientation and number of sets determine the shape of the blocks, which may
resemble cubes, rhombohedrons, tetrahedrons, sheets, or irregular shapes. The block size and
inter-block shear strength determine the mechanical behavior of a rock mass. ISRM Working
Groups (1978, 1981) suggest standard terminology for field descriptions of rock masses to
enable evaluations of block size and shape (see Table 2).
______________________________________________________________________________
71
Blocky
Approximately equidimensional
Tabular
Columnar
Irregular
Crushed
Block size (Ib) can be estimated by selecting several typical in-situ blocks, measuring them and
establishing an average size (ISRM Working Group, 1978). An alternate method of describing
block size is use of the Volumetric Joint Count method (ISRM Working Group, 1978, 1981). A
volume joint count (Jv) is defined as the sum of the number of joints per meter for each joint set
(by definition Jv is metric). The number of joints of each set should be counted along a 5-10 m
(16-32 ft) sampling interval along a line that runs perpendicular to each set.
The procedure for calculating Jv is illustrated in Figure 81. The rock mass has three
discontinuity sets, and thus the joints belonging to each set are counted along three sampling
lines; two lines(lines 1 and 2) are 10 m (32 ft) long, and one line (Line 3) is 7 m (23 ft) long.
As shown in Figure 81, the joint counts were: Line 1 = 5, Line 2 = 6, and Line 3 = 4. Then the
calculation for Jv is performed as shown:
Jv = 5/10 + 6/10 + 4/7
Jv = 0.5 + 0.6 + 0.6 = 1.7 joints /meter3
Using Table 2, the example rock mass is composed of large blocks based on a Jv value of 1.7.
______________________________________________________________________________
72
Figure 81. Schematic. Volumetric Joint Count (JV). Measurements are made normal to the direction of
joint set. Line 1 intersects 5 joints, line 2 intersects 6 joints, and line 3 intersects 4 joints.
Rock Density
The rock density was found to have little effect on CRSP-3D output. A value of 165 lb/ft3 is a
good average that produces acceptable CRSP-3D output and is the default in this program. If a
more precise estimate of density is desired, Table 3 can be used.
Table 3. Typical rock unit weights and densities (After: Hoek and Bray, 1981).
ROCK TYPE
U.S UNITS:
UNIT
WEIGHT
(lb/ft3)
160 to 190
2550 to 3060
Metamorphic rocks
quartzite, gneiss, slate
160 to 180
2550 to 2850
150 to 180
2340 to 2850
110 to 150
1730 to 2340
METRIC UNITS:
DENSITY
(kg/m3)
found to be appropriate for modeling the proper amount of rock and slope interaction, though
other experience may differ. Table 4 provides some suggestions on roughness values to be input
in CRSP-3D.
Note: It is emphasized that these values are, at best, preliminary starting points and should not be relied on - there
are no substitutes for calibration analyses and judgment.
It should be recognized that the surface roughness, when entered into the program as a single
number such as 2.0, will result in a constant surface roughness being modeled on the slope
surface. If a range of surface roughness is input, such as 0.5 - 2.0, the program will randomly
choose roughness values between 0.5 and 2.0. Single roughness values were found to produce
accurate results for all slope types, and should be used in the initial rockfall simulations to
establish suitable roughness values that result in acceptable rockfall simulation output values.
Note: The program can be sensitive to the roughness value used, and starting with a range of roughness values
without previously establishing acceptable output can result in erroneous results.
Roughness Value
3.0 and above
1.0 - 2.0
0.1 - 1.0
Comments
Generally only used for very rough slope surfaces
where high bounce heights are predicted or have
been observed.
Generally for use on most slope surfaces with most
falling rock geometries, even if the slope is
relatively smooth. This compensates for the nonuniformity in most rock shapes and slope surfaces.
Use judiciously. Roughness values below 2.0 may
produce higher than expected values when modeling
spherical rocks.
______________________________________________________________________________
74
The current program allows modeling of a hardness range or a single hardness value. A single
hardness value is used the majority of the time, and should always be used for the initial rockfall
simulations in order to get the model to calibrate satisfactorily. Only after an acceptable hardness
coefficient resulting in accurate calibrated program output has been established should a range of
hardness values be input to perform sensitivity analyses.
A hardness range could be used to model trees or boulders on a slope. For example, for a slope
with an estimated hardness of 0.5 that also has trees, instead of using a single hardness value of
0.2 to account for the cumulative effect of trees, use a hardness range of 0.1 - 0.5 to more
accurately model the random affect of trees. Each time a rock is rolled down the slope, the
hardness of the slope will vary randomly from 0.1 - 0.5, simulating a rock occasionally glancing
off of trees or hitting them directly as it bounces down the slope.
CRSP-3D maps the input hardness value to the tangential damping coefficient and the coefficient
of restitution (COR) values. The damping and COR have an inverse linear relationship. For
example, a hard slope will have a relatively high COR and a low damping value while a soft soil
slope will have a low COR and a high damping value. Slope and roadway materials have been
grouped under few general descriptions (Soft Clay, Hard Clay, Sand, Gravel, Soft Bedrock, Hard
Bedrock, Boulder field, Talus, Concrete Road and Asphalt Road) which were then used to
construct a table of preliminary roughness and hardness coefficient values (Table 5). This table is
provided only as a starting point for choosing hardness coefficient values.
Table 5. Preliminary hardness values and associated hardness measure.
SLOPE
MATERIAL
TYPE
SOFT
MATERIAL
DESCRIPTION
Soft clay
Loose sand
0.1-0.3
0.5-0.8
0.7-0.9
0.3-0.5
FIRM
0.4-0.7
Hard clay
Soft bedrock
GRAVEL TO
Gravel
COBBLE TALUS Cobbles
BOULDER
Boulder field
TALUS
HARD
Fresh hard rock
BEDROCK
0.2-0.6
CALIBRATION
The calibration process entails running rockfall simulations of existing monitored/measured
rockfalls (calibration sites) and comparing the simulated output parameters (e.g. rollout, bounce
______________________________________________________________________________
75
height, velocity) to observed rockfall parameters obtained at these sites. In this manner, input
parameters, such as slope hardness and roughness can be adjusted to simulate a rockfall more
accurately.
Two main approaches can be used to calibrate numerical models:
______________________________________________________________________________
76
REFERENCES
REFERENCES
Bartingale, R. J., 2009, Calibration of Colorado Rockfall Simulation Program (CRSP) Version
5.0. ME Thesis, Geological Engineering, Colorado School of Mines, Golden, CO, 88p.
Brawner, C.O. 1994. Rockfall Hazard Mitigation Methods Participant Workbook. National
Highway Institute Course No. 13219, FHWA SA-93-085, U.S. Department of Transportation,
Federal Highway Administration, Washington, DC.
Hoek, E., and Bray, J.W., 1981, Rock Slope Engineering, 3rd Edition. Institution of Mining and
Metallurgy, London, England, 358p
International Society for Rock Mechanics (ISRM) Working Group. 1978. Suggested Methods for
the Quantitative Description of Discontinuities in Rock Masses. Int. J. Rock Mech. Min. Sci. &
Geomech. Abst., Vol. 15, pp. 319-368.
International Society for Rock Mechanics (ISRM) Working Group. 1981. Suggested Methods for
the Quantitative Description of Discontinuities in Rock Masses. In: Rock Characterization
Testing and Monitoring: ISRM Suggested Methods (E.T. Brown, editor), Pergamon Press, pp. 352.
Jones, C.L., Higgins, J.D., Andrew, R.D., 2000, Colorado Rockfall Simulation Program Version
4.0. MI 66, Colorado Geological Survey, Denver, CO, 127p.
Pfeiffer, T.J., Higgins, J.D., Andrew, R.D., Schultz, R., and Beck, R.B,. 1995, Colorado Rockfall
Simulation Program. (Version 2.1) Colorado Department of Transportation, Denver, CO, 127p.
Pfeiffer, T.J., Higgins, J.D., Schultz, R., and Andrew, R.D., 1991, Colorado Rockfall Simulation
Program. (Version 2.1) Colorado Department of Transportation, Denver, CO, 127p.
Pfeiffer, T.J. and Higgins, J.D., 1991, Colorado Rockfall Simulation Program. (Version 1.0)
Misc. Report No. 39, Colorado Geological Survey, Denver, CO, 71p.
Pfeiffer, T.J., 1989, Rockfall Hazard Analysis Using Computer Simulation of Rockfalls. ME
Thesis, Geological Engineering, Colorado School of Mines, Golden, CO, 103p.
Ritchie, A.M. 1963. Evaluation of Rockfall and its Control. Highway Research Record 17,
Highway Research Board, National Research Council, Washington, DC, pp.13-28.
______________________________________________________________________________
77
REFERENCES
(blank page)
______________________________________________________________________________
78
HARDNESS RANGE
MATERIAL DESCRIPTION
SOFT
0.1-0.3
______________________________________________________________________________
79
HARDNESS RANGE
MATERIAL DESCRIPTION
SOFT
0.1-0.3
______________________________________________________________________________
80
HARDNESS RANGE
MATERIAL DESCRIPTION
SOFT
0.1-0.3
______________________________________________________________________________
81
HARDNESS RANGE
MATERIAL DESCRIPTION
SOFT
0.1-0.3
______________________________________________________________________________
82
HARDNESS RANGE
MATERIAL DESCRIPTION
INTERMEDIATE
0.3-0.5
______________________________________________________________________________
83
HARDNESS RANGE
MATERIAL DESCRIPTION
MODERATELY SOFT
0.3-0.5
______________________________________________________________________________
84
HARDNESS RANGE
MATERIAL DESCRIPTION
FIRM
0.4-0.6
______________________________________________________________________________
85
HARDNESS RANGE
MATERIAL DESCRIPTION
FIRM
0.4-0.6
______________________________________________________________________________
86
HARDNESS RANGE
MATERIAL DESCRIPTION
FIRM
0.4-0.6
______________________________________________________________________________
87
HARDNESS RANGE
MATERIAL DESCRIPTION
GRAVEL TO COBBLE
TALUS
0.2-0.6
______________________________________________________________________________
88
HARDNESS RANGE
MATERIAL DESCRIPTION
GRAVEL TO COBBLE
TALUS
0.2-0.6
______________________________________________________________________________
89
HARDNESS RANGE
MATERIAL DESCRIPTION
GRAVEL TO COBBLE
TALUS
0.2-0.6
______________________________________________________________________________
90
HARDNESS RANGE
MATERIAL DESCRIPTION
BOULDER TALUS
0.5-0.8
______________________________________________________________________________
91
HARDNESS RANGE
MATERIAL DESCRIPTION
BOULDER TALUS
0.5-0.8
______________________________________________________________________________
92
HARDNESS RANGE
MATERIAL DESCRIPTION
BOULDER TALUS
0.5-0.8
______________________________________________________________________________
93
HARDNESS RANGE
MATERIAL DESCRIPTION
HARD BEDROCK
0.7-0.9
Hard rock
______________________________________________________________________________
94
HARDNESS RANGE
MATERIAL DESCRIPTION
HARD BEDROCK
0.7-0.9
Hard rock
______________________________________________________________________________
95
(blank page)
______________________________________________________________________________
96
________________
PROJECT INFORMATION
Project Name
Tutorial
Material Properties
Material 1
Material 2
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
English
Boulder Field
Magenta
3-5
0.7-0.9
Talus
Gray
0.5-2
0.3-.6
Units
ft
ft
SI
Boulder Field
Magenta
0.91-1.52
0.7-0.9
Talus
Gray
0.15-0.61
0.3-.6
Units
m
Analysis Partition 2
ID
Color
Dir
Location
ID
Color
Dir
Location
English
BP-1
Cyan
Z
70
BP-2
Green
X
231
Units
ft
ft
SI
BP-1
Cyan
Z
21
BP-2
Green
X
70
Units
______________________________________________________________________________
97
Rockfall Parameters
Rockfall Parameter
1
English
Boulder Field
Units
SI
Boulder Field
Units
Release Zone
Rocks to Release
Shape
Size
165
3
One_Element_Sphere
3
pcf
25.9
3
One_Element_Sphere
0.9144
kN/m3
ft
______________________________________________________________________________
98
______________________________________________________________________________
99
______________________________________________________________________________
100
______________________________________________________________________________
101
______________________________________________________________________________
102
______________________________________________________________________________
103
______________________________________________________________________________
104
Max Energy, kJ
4
Percent Passing
67%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
0.0
0.1
0.1
0.1
0.1
Energy, ft-kips
3
3
3
3
4
Energy, kJ
4
4
4
4
5
Rollout X, ft
282
290
298
302
307
Rollout Z, ft
9
12
14
14
15
______________________________________________________________________________
105
______________________________________________________________________________
106
______________________________________________________________________________
108
Max Energy, kJ
4
Percent Passing
67%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
0.1
0.2
0.2
0.2
0.2
Energy, ft-kips
2
3
3
3
3
Energy, kJ
3
4
4
4
4
Rollout X, ft
282
290
298
302
307
Rollout Z, ft
9
12
14
14
15
______________________________________________________________________________
109
______________________________________________________________________________
110
______________________________________________________________________________
111
(blank page)
______________________________________________________________________________
112
_______________
PROJECT INFORMATION
Project Name
Canyon Lidar 3D
Sphere
Material Properties
Material 1
Material 2
Material 3
Material 4
Material 5
Material 6
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
English
Firm Soil
Yellow
1
0.3
Rock Cliff
Cyan
1
0.9
Vegetated Soil
Red
0.5
0.2
Release Zone
Green
0.5
0.9
Small Talus
Gray
1
0.4
Boulder Field
Units
ft
ft
ft
ft
ft
SI
Firm Soil
Yellow
0.3
0.3
Rock Cliff
Cyan
0.3
0.9
Vegetated Soil
Red
0.15
0.2
Release Zone
Green
0.15
0.9
Small Talus
Gray
0.3
0.4
Boulder Field
Units
m
______________________________________________________________________________
113
Material 7
Material 8
Material 9
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
DarkGray
4
0.6
Hard Bedrock
Magenta
1
0.9
Firm Soil
Yellow
1
0.3
Pavement
Blue
0.1
0.8
DarkGray
1.22
0.6
Hard Bedrock
Magenta
0.3
0.9
Firm Soil
Yellow
0.3
0.3
Pavement
Blue
0.03
0.8
ft
ft
ft
ft
Analysis Partition 2
Analysis Partition 3
Analysis Partition 4
ID
Color
Dir
Location
ID
Color
Dir
Location
ID
Color
Dir
Location
ID
Color
Dir
Location
English
AP-1
Cyan
X
809
AP-2
Cyan
X
986
AP-3
Cyan
X
1137
AP-4
Cyan
X
1273
Units
ft
ft
ft
ft
SI
AP-1
Cyan
X
247
AP-2
Cyan
X
301
AP-3
Cyan
X
347
AP-4
Cyan
X
388
Units
Rockfall Parameters
Rockfall Parameter 1
Release Zone
Rocks to Release
Shape
Size
English
Release Zone
165
20
One_Element_Sphere
15
Units
pcf
ft
SI
Release Zone
25.9
20
One_Element_Sphere
4.572
Units
kN/m3
m
______________________________________________________________________________
114
______________________________________________________________________________
115
______________________________________________________________________________
116
Max Energy, kJ
165092
Percent Passing
100%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
26.5
37.0
52.6
66.3
69.3
Velocity, ft/s
106
119
132
140
145
Energy, ft-kips
68374
85361
104479
117386
123231
Energy, kJ
92704
115735
141656
159074
167081
Rollout X, ft
1412
1479
1545
1582
1621
Rollout Z, ft
68
81
94
237
304
______________________________________________________________________________
117
Velocity (ft/s)
______________________________________________________________________________
118
Max Energy, kJ
173683
Percent Passing
100%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
11.9
15.4
25.9
41.3
48.0
Velocity, ft/s
88
99
115
132
147
Energy, ft-kips
46888
58354
78852
115491
129455
Energy, kJ
63571
79116
106906
156582
175514
Rollout X, ft
1412
1479
1545
1582
1621
Rollout Z, ft
68
81
94
237
304
______________________________________________________________________________
119
Velocity (ft/s)
______________________________________________________________________________
120
______________________________________________________________________________
121
Max Energy, kJ
105795
Percent Passing
95%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
8.6
10.9
12.3
18.9
20.5
Velocity, ft/s
71
85
97
107
116
Energy, ft-kips
29577
43156
54940
72928
79016
Energy, kJ
40098
58507
74483
98870
107123
Rollout X, ft
1418
1483
1547
1584
1623
Rollout Z, ft
67
78
88
94
101
______________________________________________________________________________
122
Velocity (ft/s)
______________________________________________________________________________
123
Max Energy, kJ
96760
Percent Passing
95%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
20.2
59.0
64.7
66.6
68.3
Velocity, ft/s
77
88
97
104
111
Energy, ft-kips
36220
44999
55133
66695
72262
Energy, kJ
49106
61007
74747
90422
97970
Rollout X, ft
1418
1483
1547
1584
1623
Rollout Z, ft
67
78
88
94
101
______________________________________________________________________________
124
Velocity (ft/s)
______________________________________________________________________________
125
______________________________________________________________________________
126
_____________
PROJECT INFORMATION
Project Name
Project Number
Company Name
Location
Designer
Reviewer
Date
Canyon 2D Sphere
1235
Acme Slope
Consulting
Colorado
D. Ziner
R. Viewer
2-25-11
Material Properties
Material 1
Material 2
Material 3
Material 4
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
English
FirmSoil
Yellow
1
0.3
Rock Cliff
Cyan
1
0.9
Vegetated Soil
Red
0.5
0.2
Vegetated Soil
Red
0.5
Units
ft
ft
ft
ft
SI
FirmSoil
Yellow
0.3
0.3
Rock Cliff
Cyan
0.3
0.9
Vegetated Soil
Red
0.15
0.2
Vegetated Soil
Red
0.15
Units
m
______________________________________________________________________________
127
Material 6
Material 7
Material 8
Material 9
Material 10
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
ID
Color
Roughness
Hardness
0.2
Release Zone
Green
0.5
0.9
Small Talus
Gray
1
0.4
Boulder Field
DarkGray
4
0.6
Hard Bedrock
Magenta
1
0.9
Firm Soil Again
Yellow
1
0.3
Pavement
Blue
0.1
0.8
0.2
Release Zone
Green
0.15
0.9
Small Talus
Gray
0.3
0.4
Boulder Field
DarkGray
1.22
0.6
Hard Bedrock
Magenta
0.3
0.9
Firm Soil Again
Yellow
0.3
0.3
Pavement
Blue
0.03
0.8
ft
ft
ft
ft
ft
ft
Analysis Partition 2
Analysis Partition 3
Analysis Partition 4
Analysis Partition 5
ID
Color
Dir
Location
ID
Color
Dir
Location
ID
Color
Dir
Location
ID
Color
Dir
Location
ID
Color
Dir
Location
English
AP-1
Cyan
X
457
AP-2
Cyan
X
626
AP-3
Cyan
X
938
AP-4
Cyan
X
989
AP-5
Cyan
X
1126
Units
ft
ft
ft
ft
ft
SI
AP-1
Cyan
X
139
AP-2
Cyan
X
191
AP-3
Cyan
X
286
AP-4
Cyan
X
301
AP-5
Cyan
X
343
Units
Rockfall Parameters
Rockfall Parameter 1
Release Zone
Rocks to Release
Shape
Size
English
Release Zone
165
20
One_Element_Sphere
15
Units
pcf
ft
SI
Release Zone
25.9
20
One_Element_Sphere
4.572
Units
kN/m3
m
______________________________________________________________________________
128
RockVelocity
RockKineticEnergyFtKips
______________________________________________________________________________
129
RockKineticEnergyKJ
______________________________________________________________________________
130
______________________________________________________________________________
131
Max Energy, kJ
188319
Percent Passing
100%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
48.1
51.8
54.2
55.5
56.9
Velocity, ft/s
137
146
153
156
160
Energy, ft-kips
116607
126677
133205
137025
141331
Energy, kJ
158095
171748
180599
185778
191616
Rollout X, ft
1129
1163
1193
1212
1233
Rollout Z, ft
68
74
81
86
92
______________________________________________________________________________
132
______________________________________________________________________________
133
______________________________________________________________________________
134
______________________________________________________________________________
135
Max Energy, kJ
153745
Percent Passing
100%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
8.4
8.8
9.2
9.4
9.6
Velocity, ft/s
119
130
137
140
144
Energy, ft-kips
84256
102457
110678
113854
116859
Energy, kJ
114236
138913
150060
154365
158439
Rollout X, ft
1129
1163
1193
1212
1233
Rollout Z, ft
68
74
81
86
92
______________________________________________________________________________
136
______________________________________________________________________________
137
______________________________________________________________________________
138
______________________________________________________________________________
139
Max Energy, kJ
212758
Percent Passing
100%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
13.7
14.6
15.3
15.8
16.3
Velocity, ft/s
149
156
162
166
170
Energy, ft-kips
133469
143669
151673
156145
160696
Energy, kJ
180952
194782
205633
211696
217866
Rollout X, ft
1129
1163
1193
1212
1233
Rollout Z, ft
68
74
81
86
92
______________________________________________________________________________
140
______________________________________________________________________________
141
______________________________________________________________________________
142
______________________________________________________________________________
143
Max Energy, kJ
268230
Percent Passing
100%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
33.4
36.5
38.9
40.7
42.7
Velocity, ft/s
171
177
183
187
191
Energy, ft-kips
175484
185376
193883
198631
203676
Energy, kJ
237920
251332
262865
269303
276143
Rollout X, ft
1129
1163
1193
1212
1233
Rollout Z, ft
68
74
81
86
92
______________________________________________________________________________
144
______________________________________________________________________________
145
______________________________________________________________________________
146
Max Energy, kJ
85122
Percent Passing
100%
Cumulative Probability
Cumulative Probability
50%
75%
90%
95%
98%
Bounce Ht, ft
8.4
24.4
31.9
35.7
39.1
Velocity, ft/s
32
88
95
99
104
Energy, ft-kips
39275
46779
50577
59745
63732
Energy, kJ
53249
63422
68572
81001
86407
Rollout X, ft
1129
1163
1193
1212
1233
Rollout Z, ft
68
74
81
86
92
______________________________________________________________________________
147
______________________________________________________________________________
148
______________________________________________________________________________
149