Ramms Rock Manual
Ramms Rock Manual
RAMMS
rapid mass movements simulation
Manuscript update
September 2022
2
Table of Content
1 Introduction ........................................................................................................................ 6
3 Theory ............................................................................................................................... 19
3.3 Free Flight Motion with Gravity and Gyroscopic Forces ........................................... 22
5
CHAPTER 1: INTRODUCTION
1 Introduction
1.1 Motivation
Mitigation of natural hazards relies increasingly on numerical process models to predict the area
inundated by rapid geophysical mass movements. These movements include
• snow avalanches,
• torrent based debris flows and hillslope debris flows,
• mudslides,
• ice avalanches and glacier lake outbreaks
• rockfalls and rock avalanches.
Process models are used by engineers to predict the speed and reach of these hazardous movements
in complex terrain. The preparation of hazard maps is a primary application. The models are especially
helpful when proposing technical mitigation measures, such as dams and embankments or rockfall
protection barriers. The models allow hazard engineers to optimize limited financial resources by
studying the influence of different hazard scenarios on defense options.
1.2 RAMMS
The RAMMS (RApid Mass Movements Simulation) software system contains three process modules:
• RAMMS::AVALANCHE
• RAMMS::DEBRISFLOW
• RAMMS::ROCKFALL
The RAMMS::AVALANCHE and RAMMS::DEBRISFLOW modules are designed for flow phenomena
containing fast moving particulate debris of snow and rocks. In the avalanche module, the interstitial
fluid is air, whereas in the debris flow module the interstitial fluid is mud. The RAMMS::AVALANCHE
and RAMMS::DEBRISFLOW models are used to calculate the motion of the movement from initiation
to runout in three-dimensional terrain. The models use depth-averaged equations and predict the
slope-parallel velocities and flow heights. This information is sufficient for most engineering
applications. Information in the slope-perpendicular direction (e.g. mass and velocity distribution) is
lost; however, this is seldom of practical interest. Both models require an accurate digital
representation of the terrain. Engineers specify initial conditions (location and size of the release mass)
and friction parameters, depending on terrain (e.g. roughness, vegetation) and material (e.g. snow, ice
or mud content of the debris flow).
The RAMMS::ROCKFALL module is used to study the rigid body motion of falling rocks. The model
predicts rock trajectories in general three-dimensional terrain. Rock trajectories are governed by the
interaction between the rock and ground. The model contains six primary state variables: three
translational speeds and three rotational velocities of the falling rock. From these, kinetic energy,
runout distance and jump heights can be derived. Generalized rock shapes are modeled. Rock
orientation and rotational speed are included in the rock/ground interaction. The RAMMS::ROCKFALL
module is therefore fundamentally different from the RAMMS::AVALANCHE and
RAMMS::DEBRISFLOW modules because it is based on hard-contact, rigid-body Lagrangian mechanics,
6
CHAPTER 1: INTRODUCTION
not Eulerian flow mechanics. It also differs from existing rockfall modules because the rock/ground
interaction is not governed entirely by simple rebound mechanics, but frictional (dissipative)
rock/ground interactions. These govern the onset of rock jumping. The RAMMS::ROCKFALL module
predicts all rigid-body motions – rock sliding, rolling, jumping and skipping.
The RAMMS::ROCKFALL module was coupled to the same user-friendly visualization tool used in the
RAMMS::AVALANCHE and RAMMS::DEBRISFLOW modules. The visualization tool allows easy
preparation, execution, visualization and interpretation of simulations.
In all RAMMS modules, new constitutive models have been developed and implemented, thanks to
calibration and verification at full scale test sites such as St. Léonard/Walenstadt (rockfall, mitigation
measures), Vallée de la Sionne (snow avalanches) and Illgraben (debris flow). At present, two new
scientific RAMMS modules are under development: RAMMS::AVAL_EXTENDED and
RAMMS::DBF_EXTENDED.
The RAMMS web page [Link] provides useful information such as a moderated discussion
forum, frequently asked questions (FAQ) or recent software updates. Please visit this web page
frequently to stay up to date.
7
CHAPTER 1: INTRODUCTION
DISCLAIMER
RAMMS is intended to be used as a tool to support experienced users. The interpretation of the
simulation results has to be done by a rockfall expert who is familiar with the local as well as with
the topographic and geological situation of the investigation area. In no event shall SLF/WSL be liable
for any damage or lost profits arising, directly or indirectly, from the use of RAMMS. Swiss law
applies. Court of jurisdiction is Davos. If you encounter problems, please contact ramms@[Link].
8
CHAPTER 2: INSTALLATION AND SETUP
2.2 Installation
Please download the RAMMS::ROCKFALL setup file “ramms_rock_user_setup_64.zip” from
[Link] (Downloads section). Please make sure that you have a 64-bit Windows system.
• Click on the path given above or copy the path to any browser. A window pops up and the
automatic download of the file ramms_rock_user_setup_64.zip starts after clicking Yes.
• Unzip the file to a temporary location.
• You must have Administrator privileges on the target machine. If you do not have such privileges,
the installer cannot modify the system configuration of the machine and the installation will
fail. Note that you do not need Administrator privileges to run RAMMS afterwards.
• Read first, install afterwards! Please read the whole installation process once, before you begin
the installation.
• Start the file “ramms<version>_rock_user_setup_64.exe”.
9
CHAPTER 2: INSTALLATION AND SETUP
Step 1: Welcome
The welcome dialog introduces you to the English setup program and will guide you through the
installation process. Click Next to continue.
Step 2: Readme
10
CHAPTER 2: INSTALLATION AND SETUP
Read the license agreement carefully and accept it by activating the check box in the lower left
corner. If you do not accept the license agreement, you are not able to proceed with the
installation. After accepting the license agreement, click Next to continue the installation.
Choose your destination directory. This dialog shows the amount of space available on your hard disk
and required for the installation. Click Next to start the installation process.
11
CHAPTER 2: INSTALLATION AND SETUP
RAMMS is copying the files to the destination location. The window shows the installation progress.
RAMMS finished copying the files. Click Next to finish the installation process.
12
CHAPTER 2: INSTALLATION AND SETUP
To ensure that all important system libraries are installed on your target machine follow the
instructions below:
The welcome dialog introduces you to the English setup program and will guide you through the
installation process of the IDL Visual Studio Merge Modules. Click Next to continue.
Figure 2-8: IDL Visual Studio Merge Modules - welcome dialog window.
13
CHAPTER 2: INSTALLATION AND SETUP
Figure 2-9: IDL Visual Studio Merge Modules - ready to install the program.
14
CHAPTER 2: INSTALLATION AND SETUP
After having successfully installed RAMMS and the necessary files on your personal computer, you
will notice the RAMMS icon on your desktop (for all users):
Additionally, a new application folder is created in Start → Programs (for all users):
15
CHAPTER 2: INSTALLATION AND SETUP
2.3 Licensing
Access to RAMMS is controlled by a personal use license. Personal use licenses are time limited licenses
tied to a single personal computer. This method of licensing requires a machine’s unique host ID
to be incorporated into a license request file. After the license request file is sent to SLF/WSL, you
will receive a license key. Entering the license key on a personal computer enables full RAMMS
functionality for the specific personal computer. For more information please visit
[Link] Alternatively, the license can also be installed on a Windows Server and
accessed by different users (only one at a time) by RDC (Remote Desktop Connection).
Double-click the RAMMS icon or use Start → Programs → RAMMS Rockfall → Run RAMMS Rockfall
to start RAMMS for the first time. Whenever you start RAMMS, the splash screen below will pop
up:
Click on the image. It will disappear and RAMMS will start up. The following dialog window appears:
16
CHAPTER 2: INSTALLATION AND SETUP
Click the button to create your personal license request file. In Figure 2-16 enter your full
name and the name of your company.
In the next dialog window, choose the destination directory of your personal license request file and
save it to your target machine. Your personal license request file should look similar to Figure 2-17.
An order confirmation email is sent to your email address. We then process your order and send
you an invoice. As soon as we received your payment, we will send you your personal license key.
Your personal license key is named similar to ROCK_Muster_Test_RAMMS.txt. Open the file in a
text editor. It should look similar to Figure 2-18.
17
CHAPTER 2: INSTALLATION AND SETUP
Now, restart RAMMS (as explained before). The IDL splash screen appears (Figure 2-14) and then the
dialog window of Figure 2-15 shows up (RAMMS - Licensing). Copy the license key (in this example:
ROCKFALL ebei-flhl-ilkq-behe-1i5m) and paste it at the field LICENCE KEY (see Figure 2-15). Notice
that there is the prefix ROCKFALL. This prefix is part of the license key and has to be inserted as
well! If RAMMS accepts your installation key, you successfully finished the installation.
2.4 Update
When you start RAMMS it will automatically check for updates on the internet. This can lead to an
error message, if your firewall blocks the executable [Link] (this file starts the IDL-Virtual Machine
you need to run RAMMS). Please unblock this file for your firewall. You can also disable the
AutoWebUpdate-function by unchecking Help → Advanced... → AutoWebUpdate. In the same way
you can enable the AutoWebUpdate-function by checking Help → Advanced... → AutoWebUpdate.
18
CHAPTER 3: THEORY
3 Theory
3.1 Overview
The RAMMS::ROCKFALL model utilizes a hard- contact, rigid-body approach to model rockfall
trajectories in general three-dimensional terrain (Leine et al., 2013). The program is designed to be
used by hazard engineers to predict rockfall velocity and runout for hazard mapping and planning of
rockfall mitigation measures. The calculation engine and user interface were developed as part of a
joint research project between the WSL Institute for Snow and Avalanche SLF and the Institute of
Mechanics, Swiss Federal Institute of Technology (ETHZ) between the years 2010-2013. The rockfall
model is the third RAMMS module, following the RAMMS::AVALANCHE and RAMMS::DEBRISFLOW
modules and offers many of the same user-friendly features. The RAMMS::ROCKFALL model was
officially released after a period of calibration and application testing in April, 2015.
To date most rockfall models utilize simple rebound mechanics to describe the complex interaction
between the rock and the ground (Bourrier et al. 2012; Dorren 2003; Dorren and Seijmonsbergen 2003;
Schweizer 2015, Volkwein et al. 2011). Rock geometries consisted of simplified shapes, mostly spheres
or ellipsoids. The rock-ground interaction was parameterized using apparent restitution coefficients to
model the rock jumping. To account for the wide variation of possible jump distances and heights (even
in homogenous terrain, see Glover 2015), random, stochastic methods were used to define the
bandwidth of possible restitution coefficients. Rockfall modeling was therefore both quasi-
deterministic and quasi-stochastic.
In RAMMS the rock-ground interaction is parameterized by frictional operators that act at the rock
surface. Compared to rebound models (that employ apparent restitution coefficients to model
entire ground-rock interaction), the hard-contact, rigid-body approach applies contact forces to the
rock’s edges and corner points. The primary advantage of using hard-contact approach is that the
role of rock shape is accounted for in the ground-rock interaction. This facilitates a natural modeling
of the four primary modes of rock motion: sliding, rolling, skipping and jumping – without the use of
random, stochastic methods to define the rebound parameters. All four modes of rock propagation
are modeled in RAMMS. Long and widespread rock runout is generally associated with the jumping
mode; however, rock stopping requires a transition from jumping to a rolling/sliding mode. Modeling
all four modes is essential for a realistic, self-consistent and risk-based rockfall hazard analysis. The
natural variation of jumps is defined automatically by the rock shape and orientation at impact. The
statistical spread of rockfall runout and dispersion is generated only by changing the initial conditions.
Ground parameters are not random: they are deterministic in the sense that one material type is
assigned to describe hardness and the general tendency of the terrain to react to a rock impact. In
RAMMS the clear separation between stochastic initial conditions and deterministic boundary
conditions simplifies and enhances the construction of engineering-based hazard scenarios and the
interpretation of model results. At present, the RAMMS rockfall model contains six ground categories:
soft, medium soft, medium, medium hard and hard. A seventh category snow has been introduced to
model rockfall motion on snow in high mountain terrain.
Applying RAMMS to a rockfall problem necessitates that rocks of different shapes and sizes can be
easily specified. In RAMMS::ROCKFALL the rock-body is modelled as a convex hull polyhedron. The
19
CHAPTER 3: THEORY
shape of a rock body is user defined by providing a point cloud file, defining the surface geometry of
the rock. Shapes can be simple geometric forms, such as equant, platy or columnar. A unique feature
of RAMMS is that real rock geometries obtained from laser scans during field investigations can be
used in a modeling application. Over time, the user can build-up and manage a rockfall library
containing rock shapes representative of different geologic settings. At present, the rocks are
considered indestructible; that is, they do not fragment or change form during the analysis.
Over and above weathering processes, the geometric relationships of rock-mass discontinuities
(joints, fractures, contacts, bedding, asperities and schistosity) govern block shape, size and release
mechanism (Jaboyedoff, 2011). With RAMMS::ROCKFALL preconditioning the shape and size and
number of possible release orientations of detachable rocks is an essential part of the analysis. The
observation that different basic geological settings produce characteristic rock shapes has been well
documented by Fityus et al. 2013, among others. Some commonly encountered rock shapes and
the associated geological setting are given below (Figure 3-1).
Figure 3-1: Photographs of rock masses and their aggregate forms. Top left: An example of equant
cubic rock forms generated in a sequence of sandstones exposed to an extensional deformation
regime, the primary joint sets are near equally spaced and orthogonal to one another. Topright:
The complex joint of this granodioritic rock mass results in highly irregular and angular rock block
forms. Bottom left: The uplifted and folded limestone sequence is well bedded producing
distinguished slabs which detach as pronounced platy rock forms. Bottom right: Distinguished
columnar jointed basalt sequence produces the characteristic elongate rock forms (Glover 2015).
The specification of general rock geometries will be discussed in the next section 3.2.
Another feature of the RAMMS::ROCKFALL model is the inclusion of rock rotations in both the airborne
phase and during the interaction with the ground. The RAMMS::ROCKFALL model includes gyroscopic
20
CHAPTER 3: THEORY
forces induced by rock rotations. These forces are necessary to model wheel-like rock skipping and
jumping modes that are often responsible for extreme runout. To model ground interaction
considering rocks with arbitrary geometry and rotational speed requires methods to accurately track
the rock orientation relative to the ground. RAMMS::ROCKFALL employs quaternion algebra for this
purpose. This method tracks rotation sequences even when non-linear contact forces change the
translational and rotational direction of the rock. Modeling the rock-ground contacts in this way
permits the entire mechanics of an impact to be simulated deterministically. The moment arms and
torques responsible for how different rock-shapes convert translational movement into angular
momentum and influence rebound heights are computed, and therefore allows an accurate
modeling of rolling, skipping, sliding and jumping.
The three-dimensional motion equations including rock rotations and gyroscopic forces will be
presented in section 3.3.
Complex mountain terrain is modelled using a high- resolution digital elevation model (DEM). The
specification of the DEM will be discussed in more detail in section 3.10.
21
CHAPTER 3: THEORY
Figure 3-2: Laser scans of real rocks are captured in the field. The point cloud representing the
rocks geometry are then used by the rockfall model to create a convex-hull polyhedron
representative of the rock-body.
Figure 3-3: Rock is generated from a point cloud and converted into a rigid-body polyhedral.
where M is the constant and diagonal mass matrix (containing the mass and three moments
of inertia I ). The vector 𝑢̇ contains the rock’s three translational and three rotational velocities.
The rock-body’s motion is governed by a number of forces which determine its trajectory.
Gravitational force (Fg) acts globally; a drag force (D) is implemented to represent the effects of
trees, undergrowth and soil deformation. Along with gyroscopic forces G which can cause rocks
of irregular shape to become upright and rotate about a rolling axis. All force terms h are a
function of the rock’s position q and velocity u forming the force vector h:
F + D
h ( q, u ) = g (3.2)
G
22
CHAPTER 3: THEORY
The contact of the rigid rock-body is detected by continually measuring the vertical gap length
g N between the rock-body’s corner points (P) and the terrain projections (Q) (Figure 3-4). The
gap length is defined as
g N ( X , Y , Z ) = Z − Z m ( X m. , Ym ) (3.3)
Then, when g N 0 there is no contact; when g N 0 there is contact and the contact forces , acting
at the contact point P, are computed. (The contact forces are denoted using the Greek letter lambda
because the contact forces are Lagrangian multipliers that enforce the non-penetration constraint).
Minimal penetration with the terrain is permitted to allow the assessment of the contact condition
(Eq. 3.3).This is a non-physical penetration and purely for numerical purposes.
Contact forces are modeled as hard unilateral constraints with Coulomb friction using non-
smooth contact dynamics approaches (see Acary and Brogliato, 2008, Glocker, 2001 and Moreau,
1988). For the case of contact, the governing equations of motion now become
where the direction of the contact forces is given by W (q) . There can be a number of active
contact forces depending on the rock-body’s configuration at the point of contact. Ultimately it is
the combination of these forces λ (and force directions W (q) ) that allows the complex rotations
and trajectory deviations that are inherent to rockfall to be simulated.
23
CHAPTER 3: THEORY
The advantage of this hard-contact rigid-body approach is that the contact forces are applied
directly about these contact points, respecting the configuration (orientation and kinetics) of the
impact. This is achieved by considering the contact pair (Q, P) within the contact frame C = (n,t1 ,
t2 ) which is attached to the terrain surface at contact point Q (Figure 3-5).
24
CHAPTER 3: THEORY
Figure 3-5: Contact frame C at point Q detected with the gap function gN.
The normal force component λN is resolved with a contact cone differential inclusion, in which the
transient normal force vector over the finite contact period can be computed. Over the contact
period this is a set-valued normal force considering all periods of contact identified with the gap
function gN.
The tangential force component λT is assumed to obey spatial Coulomb’s friction law (see Figure 3-6).
Stiction of the contact γT = 0 occurs as long as the magnitude of the tangential force ‖𝑘𝜆𝑇‖ is less
than µλN, in which λN is the applied normal force and µ the friction coefficient. The direction is also
resolved with a normal cone inclusion projecting a friction disc on to the surface (Figure 3-6). The
formulation covers both sticking and sliding cases.
25
CHAPTER 3: THEORY
N+ + N− = 0 (3.5)
Impulsive normal forces can also induce impulsive tangential forces. While this is mainly seen in the
elastic impacts of superballs (Cross, 1999), and therefore in the rockfall model εT is set at εT = 0
since these effects are absent.
To determine the resultant force direction acting on the rock-body the configuration of the impact
must be computed. This requires finding the relative velocity between the contact points P and the
terrain Q. Importantly, the velocity of contact point P is composed of the translational velocity with
respect to the body’s center of mass vS and its angular velocity KΩ in the fixed body frame (K); for
which P also has a fixed position vector relative to the center of mass S. That is, the contact algorithm
in the rigid-body approach considers the rotational speed of the rock at contact. Because the forces
are then applied at points away from the center of mass, and with a direction respecting the impact
configuration to a body with three degrees of translational and rotational freedom, torques and
moment, arms can act generating rotations and rebounds that represent the true mechanics of an
impact.
26
CHAPTER 3: THEORY
friction force is of great importance because it controls when the rock slides, rolls or jumps. Drag,
on the other hand, acts at the rock’s center of mass and therefore creates no additional
rotational moments. It acts in the direction opposite to the rock’s movement (velocity). There are
two drag forces in the RAMMS::ROCKFALL model. The first accounts for vegetation drag; the
second accounts for the viscoplastic drag due to terrain deformation during ground contact.
Figure 3-7: Rock impact scar on soft soil; the scar morphology is tapered widening towards the
accumulation of earth at the scar end, where an earth ramp structure is formed. This is
modeled as a climbing friction from the beginning of the scar s = 0 at first contact which tends
towards high friction at the end of the scar.
To simulate ground deformation within the framework of a hard contact model requires introducing
a slip (s) dependent friction that acts during sliding and accounts for the increase in friction due to
material accumulation behind the rock- body as it slides through the impact (Figure 3-7 and Figure
3-8). The slip dependent friction is an extension of the Coulomb friction model in which the friction
value µ is made dependent on the slip distance (s) travelled by the center of mass µ(s) (Figure 3-8).
27
CHAPTER 3: THEORY
Figure 3-8: Sliding friction in RAMMS is governed by a slip-dependent material law. At rock impact
slip is s=0 and sliding friction is given by min; with s>0, friction increases according to the
coefficient . At some point s the maximum friction max is reached. After contact, the friction
exponentially decreases with coefficient . Therefore describes the duration of the friction as the
rock is leaving the scar (ramping).
Moreover,
T = (s) N (3.6)
force λN enforces the non-penetrability constraint; the force λT acts tangentially on the terrain
surface (see Figure 3-6). The dependence of the friction coefficient on the slip distance (s) is:
2
( s) = min + ( max − min ) arctan( s ) (3.7)
where µmin, µmax and κ are parameters of the friction model. The initial friction encountered at the
contact where s = 0 is µmin. Over the slip period, µ tends toward µmax for large slip values, see
Figure 3-8. The parameter κ controls how quickly the friction increases from µmin to µmax. Typically
µmin < µmax meaning that the friction increases the longer the rock is in contact with the ground. It
is entirely possible that there are brittle ground materials where the opposite behavior (µmax >
µmin) is encountered.
28
CHAPTER 3: THEORY
The slip distance (s) is a transition state variable having a time evolution which is described by a
simple differential equation:
v if g N 0
s = s (3.8)
− s gN 0
The parameter β controls how quickly the friction is released as the rock departs the ground scar.
If β is large, friction is immediately removed as the rock moves away from the ground. Conversely,
when β is small, sliding friction can act, even after the rock is no longer in contact with the ground.
The parameter β is linked to the penetration depth of the rock into the ground. Larger penetration
depths (softer materials) are associated with smaller β values.
m
Fv = − C v v s2 (3.9)
2
Ground drag acts when the rock is in contact with the ground (gN < 0) as the rock is sliding on the
2
terrain surface (s > 0). The drag force Fv is proportional to the square of the rock velocity v s as
well as the mass of the rock m . That is, heavier and faster moving rocks will experience more drag
than smaller, slower moving rocks, as they penetrate the ground surface. The drag force is
proportional to the rock’s total kinetic energy. The drag coefficient C v varies between 0.0 m-1 (hard)
and 1.0 m-1 (soft).
29
CHAPTER 3: THEORY
3.8 Forest/Vegetation
Forest drag is given by (Figure 3-9):
if Z Z h
Fdf = −C f vs Cf = f (3.10)
0 Z Zh
The idea behind forest drag is that a resisting force acts on the rock’s center of mass when it is located
below the drag layer height Zh. This force is linearly proportional to the rock velocity v s . The forest is
parameterized by the effective height of the vegetation layer Zh as well as the drag coefficient f . The
effective height Zh roughly corresponds to the height of the forest but in some cases, for example old
forests, the drag force in the tree crowns might be negligible and therefore the effective height could
be smaller than the real tree height. The model does not account for a Z-dependency in forest structure
as it assumes a homogeneous layer with mean drag properties. Typical values for Zh are between 5 m
and 30 m (default value is 30m); typical values for f range between 100 kg/s and 1’000 kg/s. Three
different forest types are implemented in RAMMS::ROCKFALL for now:
This model description is a simplified summary of a more complete and detailed explanation provided
by Leine et al. (2013) in which the time stepping methods are also explained. While it is possible with
the rigid-body approach to model multi-body interaction of many particles, and it would be possible
to include a fragmentation law in the model, these features are not yet implemented in the model and
remain as future additions to the model.
30
CHAPTER 3: THEORY
Figure 3-9: Forest drag Fdf is implemented to act on the center of gravity of the rock-body at
height Z.
• Extra Soft
• Soft
• Medium Soft
• Medium
• Medium Hard
• Hard
• Extra Hard
• Snow
31
CHAPTER 3: THEORY
Selection of the appropriate terrain material model is the primary task of the hazard engineer when
using RAMMS. If there is uncertainty about the specific terrain material we recommend the user to
compare the results of different terrain scenarios. Please consider the examples in Table 3.1 as a
preliminary suggestion.
Rocks penetrate
meadow surface
leaving impact
Medium Soft scars. Soil is deep, Meadow
few rock
fragments. Rank
vegetation.
32
CHAPTER 3: THEORY
Meadow is deep,
but contains rock
fragments. The
meadow can be
covered with
Medium vegetation. Soil Meadow
structure of a
medium
deepness. Rank
vegetation.
Penetration
depths are small.
Ground is flat. Non-paved
Medium Rocky debris is mountain roads,
Hard present. Shallow mountain
surface soil. meadow, pebble
Usually little
(initial) vegetation.
33
CHAPTER 3: THEORY
Ground is very
hard and is
marginally
Extra Hard deformed by Bedrock, cliff
rocks. No
vegetation and no
surface soil.
Rocks slide on
Snow Snow
snow surface.
The used parameters for every terrain material are explained in the table below:
34
CHAPTER 3: THEORY
With this current version, we urge the user to go for softer categories when the rocks get large (see
table below, considering the bigger scarring effects and larger energy dissipation of larger rocks).
Moving towards softer soil for larger rocks is always the preferred direction. The block sizes have to be
compared with the assumed grain size distributions of the terrain material on a qualitative basis. In
Table 3.3, suggestions of terrain categories for typical alpine geomorphological terrain units are listed
for different rock sizes. Please note, bare rock surfaces should be chosen as Hard to Extra Hard,
regardless of rock size.
Table 3.3: Suggestions of terrain categories based on empirical test assuming a 2m DEM resolution
Typical morphological Thin, fine Thick, coarse Thick, medium- Cultural land
alpine units grained Alpine grained grained low- (vegetated)
talus Alpine talus altitude talus
Block size scenario
0.1-0.5 m3 Medium hard - Medium hard Medium Medium
Medium - Medium
1-5 m3 Medium Medium – Medium soft Medium soft -
Medium soft Soft
5-50 m3 Medium Medium soft Soft Soft – Extra
Soft
With this deterministic modeling approach, the influences of rock-shape on the dynamics and
persistence of runout can be simulated. This is important because the model is highly sensitive to
rock shape, which in the case of rebound approaches has to be treated with stochastic methods
(Bourrier et al., 2009). The role of rock-shape in runout dynamics is crucial in determining the
rotational and rebound behavior. For specific rock-shapes, characteristic of geological zones,
distinctive runout behavior such as extreme jump heights and runout distances are observed.
Dynamics of this kind are decisive for hazard mapping and rockfall protection structures; and with
full three-dimensional data of rock position, velocities, rotations and energies (see Table 3.4 for a
full list), rockfall management and the design of protection structures can be optimized. The listed
data (Table 3.4) are available as *.log files which can be generated from simulations. The application
of rigid-body theory to rockfall modeling has advanced the capacity to include detailed and hazard
specific information on rock-shapes and sizes. This allows the inclusion of lithology and geological
setting to establish realistic initial conditions for a hazard simulation.
36
CHAPTER 3: THEORY
Figure 3-10: High resolution three dimensional terrain model forms simulation frame O in which
the four sided planes form the tessellated terrain surface with which the rock-body can come into
contact.
DEM resolution
A minimal cell size of 2-5m should be used to obtain realistic results. The cell size can be increased
in case of long runout distances of > 2 km and scenarios with big block sizes (>10-50 m3). If the
project perimeter is rather small, corresponding runout distances short (<50m) and the
variability of the terrain is high within a few meters, then the cell size of the DEM should be
lower, e.g. 1.0 or even 0.5m.
37
CHAPTER 4: SETTING UP A SIMULATION
4 Setting up a simulation
4.1 Preparations
To successfully start a new RAMMS project, a few important preparations are necessary. Topographic
input data (DEM in ASCII- or GEOTIFF-format), project boundary coordinates and georeferenced
maps or remote sensing imagery should be prepared in advance (.tif format and .tfw-file, maps and
imagery are not mandatory, but nice to have). Georeferenced datasets have to be in the same
Cartesian coordinate system (e.g. Swiss CH1903 LV03) as the DEM. Polar coordinate systems in
degree (e.g. WGS84 Longitude Latitude) are not supported. For more information about specific
national coordinate systems please contact the national topographic agency in your country.
ASCII X, Y, Z data (regular and irregular) can be converted within RAMMS into an ASCII or GEOTIFF grid.
A wizard will guide you through the conversion process. The following interpolation methods are
available: LINEAR or INVERSE DISTANCE.
The header of an ESRI ASCII grid must contain the information shown in Figure 4-1.
Figure 4-1: Example ESRI ASCII grid. Figure 4-2: Example ASCII X, Y, Z single space
data.
38
CHAPTER 4: SETTING UP A SIMULATION
Figure 4-3: The same project extent (area of interest) can be used to calculate different scenarios
with different input parameters.
4.2 Preferences
To ease the file handling, we recommend setting the preferences prior to start with simulations. The
preferences set the path to the working directory and the necessary files such as DEM, maps and ortho-
imagery. If the path to the maps and the imagery files is set correctly in the preferences, RAMMS will
automatically open the georeferenced data when you generate a project.
Use Track → Preferences to open the RAMMS preferences window or click the button .
For resetting the general preferences use Help → Advanced… → Reset General Preferences.
Figure 4-4: General tab of RAMMS Figure 4-5: Rockfall tab of RAMMS
preferences. preferences.
39
CHAPTER 4: SETTING UP A SIMULATION
General Tab
Setting Purpose
Working Directory Set your working directory. VERY IMPORTANT: DO NOT USE BLANKS
in the working directory path!
Map Directory Set the folder where you place your georeferenced digital maps
(consists of a .tif file and a corresponding .twf file (world-file).
Orthophoto Directory Set the folder where you place your digital georeferenced
orthophotos (aerial picture, consists of a .tif file and a corresponding
.twf file (world-file).
DEM Directory Set the folder where you place the Digital Elevation Models (format
ASCII grid)
Rockfall Tab
Setting Purpose
Nr of Colorbar Colors Set default number of colorbar colors.
Rock Magnification *X Set values between 1 and 100 for magnification of the rock size in
the visualization.
GIF-Animation Interval [s] Set interval for GIF animation images in seconds.
Background Color Set background color.
Animation Delay [s] Set animation delay to the animation speed.
The following exercise Working directory shows how to choose a new working directory. All further
settings can be changed in a similar manner. The settings are saved, until they are changed again
manually.
40
CHAPTER 4: SETTING UP A SIMULATION
Choosing the right working directory is very useful and saves a lot of time searching for files and
folders.
VERY IMPORTANT: Do NOT use blanks or special characters in the path names!
• Click (or use Track → Preferences or Ctrl+P) to open the RAMMS preferences window.
• Click into the field Working directory. A window pops up where you can choose your new
working directory. Click OK in both windows. Do this also for other directories if necessary.
Figure 4-6: RAMMS preferences Figure 4-7: Browse for the correct folder.
41
CHAPTER 4: SETTING UP A SIMULATION
• Click or Track → New... → Project Wizard to open the RAMMS Project Wizard.
• The following window pops up.
42
CHAPTER 4: SETTING UP A SIMULATION
Step 1:
• Enter a project name (1)
• Add some project details (2)
• The project location (3) suggested is the current working directory. To change the location
click into the Location field. A second window appears and you can browse for a different
folder (see figure below, VERY IMPORTANT: Do NOT use BLANCS or special characters in
the project location path!).
• Click Next (4)
Figure 4-9: Step 1 of the RAMMS Project Figure 4-10: Window to browse for a new
Wizard Project Information. project location.
Step 2:
• Click into the Select DEM-file field to
browse for the DEM file. Locate your
DEM file in the folder set in the RAMMS
preferences.
• The Grid Resolution field shows you the
dimension of a single grid cell.
• Click Next.
43
CHAPTER 4: SETTING UP A SIMULATION
Step 3:
Enter the X- and Y-coordinates of the lower left and upper right corner of your project area,
using the Swiss Coordinate System CH1903 LV03 (or another Cartesian coordinate system), as it is
shown below for the Vallée de la Sionne area.
Step 4:
• Check the project summary, especially if
a DEM-file was found.
• To make changes click Previous, to
create the project click Create Project.
• If several matching .tif-files exist,
RAMMS shows a list with all these files.
Project creation:
The creation process can take a while. Different status bars will pop up and show the progress
of the project creation process.
44
CHAPTER 4: SETTING UP A SIMULATION
Figure 4-15: Files and directories created with a new RAMMS::ROCKFALL project.
Table 4.1: Listing of files and directories created with a new RAMMS::ROCKFALL project.
File / Folder Purpose
doc (folder) Folder containing input and output log files
logfiles (folder) Project creation and calculation log files
output (folder) Folder containing calculated scenarios
rocks (folder) Folder to save rock files (.pts)
[Link] GEOTIFF raster file with altitude values
[Link] Internal binary file containing DEM information
[Link] Calculated slope angles of DEM
[Link] Calculated planar curvatures of DEM
_.xml Input file
_.xyz Topographic data used in RAMMS
45
CHAPTER 4: SETTING UP A SIMULATION
Figure 4-16: Active project with lines and corners for resizing.
• To move the model without changing size or aspect ratio, move the cursor to the
model and check if the cursor turns to . Then click and hold the left mouse
button and drag the model to the desired position.
• To resize the model without changing the aspect ratio, use the mouse wheel to zoom
in or out. Alternatively, you can resize the model by changing the percentage value in
the horizontal toolbar .
46
CHAPTER 4: SETTING UP A SIMULATION
After activating the rotation button , the model can be rotated along the rotation axis’,
by moving the cursor directly on one of the axis until the cursor changes from to . Otherwise
a freehand rotation in any direction is possible.
Click to switch from 3D to 2D view. This button then changes to and by clicking again, you
will return to 3D view.
Figure 4-18: 3D view of example model. Figure 4-19: 2D view of example model.
In 2D mode you have all possibilities that work for the 3D mode. It works for input files as well as
for simulations. For the following functions of RAMMS it is necessary to switch from 3D to 2D
view:
INPUT: OUTPUT:
47
CHAPTER 4: SETTING UP A SIMULATION
4.4.2 Colorbar
As soon as a parameter is shown in the project, the colorbar appears on the right side of the
main window. It can be turned on and off by clicking on . The colorbar can be moved
anywhere in the screen (and can get lost). Use Project → Get Colorbar to find a lost colorbar.
Changing the minimum and maximum values of the colorbar as well as changing the number
of colors used is done in the panel ROCKFALL (right of the map window) in the tab Display.
48
CHAPTER 4: SETTING UP A SIMULATION
It is possible to change the map or imagery of a project anytime. Take into account, that the
corresponding .tfw-file (world-file) has to be in the same folder as the actual map (.tif ). If this is not
the case, the map will not be found!
• Information on the image dimensions (x-Dim and y-Dim, pixel) and size (in MB) are
provided and might be a selection criterion.
• Select the map you wish to add and click Load selected map.
49
CHAPTER 4: SETTING UP A SIMULATION
Exercise 4.3f : How to save input files and program settings manually
a. Input file:
• In case you want to save the input file manually before running a simulation, go to Track
→ Save. This is helpful, when for example a release line was loaded but you wish to close
the project before doing the simulation.
• If you wish to save a copy of your file under a new name, go to Track → Save
Copy As or click .
• A window pops up to choose an old file which should be overwritten or to type in a new
name, then click Save.
• Continue working on the original file, not the just saved one.
b. Position settings
• If you have moved and/or rotated your project for a better view, you can save this
position by going on Extras → Save Active Position.
• You can now get back to this position anytime by choosing Extras → Reload Position.
50
CHAPTER 4: SETTING UP A SIMULATION
51
CHAPTER 4: SETTING UP A SIMULATION
Release Point
1. Click or use Input → Release... → Set New Release Point, move the cursor to the
desired position and click with the left mouse button at the position of the release
point.
2. Use Input → Release... → Enter Coordinates (X/Y) to manually enter the coordinates
(X/Y) of a release point.
Save your release point location with Input → Release... → Save Point Location. See exercise 4.4a
below for more details.
Release Line
Draw a new release line by clicking or use Input → Release... → Draw New Release Line. Draw
the line by clicking points with your left mouse button and finish the line by a click on your right mouse
button.
52
CHAPTER 4: SETTING UP A SIMULATION
Release Area
Draw a new release area by clicking or use Input → Polygon Shapefile... → Draw New Polygon
Shapefile. Finish the release polygon with right-click. It is possible to define several release areas
within the same shapefile. There are two possibilities to load existing release areas:
• Click .
• Click into the project where you want to define your rockfall release point (1).
• The lower right status bar then displays the position of the release point within the
terrain (2).
• To save the coordinates of the release point go to Input → Release → Save Point
Location and enter a file name.
53
CHAPTER 4: SETTING UP A SIMULATION
• Click .
• Click into the project where you want to start drawing the outline of the release
line.
• Continue drawing the release line by moving the cursor and clicking the left mouse
button. If you would like to delete one step of the drawing, click the middle mouse
button.
• To end the release line, click the right mouse button.
Before the release line is created, you have to answer the following questions:
• Add more polylines? You can either answer with Yes and create a second release line
as explained above or answer with No and continue with the next step.
• Choose a new release line file name: e.g. line_rel, the addition _rel helps the user to
run a simulation because the model recognizes the shapefile as a release line.
54
CHAPTER 4: SETTING UP A SIMULATION
55
CHAPTER 4: SETTING UP A SIMULATION
• Select a predefined rock shape from the rock library (1) or select a .pts-file of a
rock shape from another source (2).
• By pressing the left mouse button and moving the mouse you can move the 3D
visualization of the rock interactively and look at it from any direction (3).
56
CHAPTER 4: SETTING UP A SIMULATION
• Enter a file name or use the suggested name and click to save the new .pts-file (6).
RAMMS automatically creates a rocks folder in your project directory. It’s strongly
recommended to save your new rocks in this directory.
57
CHAPTER 4: SETTING UP A SIMULATION
• Terrain Material
• Forest areas
Exercise 4.6a below demonstrates how to draw such a terrain shapefile. Figure 4-33 then explains how
to assign terrain types to selected shapefiles. Each terrain category defines the parameters of the
rockfall slippage friction law as well as a ground drag value. The ground drag value accounts for the
viscoplastic drag due to terrain deformation during ground contact. The terrain classes are described
in detail in Table 3.1. It is important to note, that you also have to consider rock size when choosing
different terrain types, see Table 3.3 for more information about rock size vs. terrain categories.
The friction parameters and drag forces should be defined for every terrain material shapefile. The
list below gives an overview on some possible terrain materials. You can choose between the
categories: Snow, Extra Soft, Soft, Medium Soft, Medium, Medium Hard, Hard and Extra Hard. The
used values have a large impact on the simulation results and therefore it is important to critically
think about the chosen parameters. The Exercise 4.6a shows how to define these terrain material
parameters. Fehler! Verweisquelle konnte nicht gefunden werden. gives an overview on the different
predefined parameters that are used for the categories. Please use these parameters with the
awareness that they are currently based on case studies. It remains an ongoing research task to
reassure these parameters rigorously.
ASCII-files exist for every friction parameter µmin, µmax, κ (average gradient from µmin to µmax after
ground contact of a rock), β (average gradient from µmax to µmin after rock loses contact with ground),
and forest drag, if the shapefiles in the friction tab or forest tab are specified. If only the overall terrain
material is specified (and no friction and forest shapefiles) no ASCII-files will be created and used.
Before the polygon shapefile is created, you have to answer the following questions:
• Add more polygons?
You can either answer with Yes and create a second polygon as explained above or answer
with No and continue with the next step.
• Choose a new polygon file name:
Enter a new name according to the terrain material represented by the polygon(s) (e.g.
bedrock). The ending *.shp is added automatically. The polygon shapefile will now be
created and opened directly.
• Create polygon shapefiles as described above for all the different and important
terrain materials inside the area of interest. These polygons can then be loaded in
the Terrain Tab in the Run Simulation dialogue window (see Figure 4-33).
• Alternatively you can draw the polygon shapfiles in a GIS software and load them
directly in the Run a simulation dialogue.
59
CHAPTER 4: SETTING UP A SIMULATION
RAMMS::ROCKFALL will apply a linear viscous drag force, which is acting only within the forest areas.
Detailed information about the forest drag forces can be found in Chapter 3.8.
The difference between the forest types considers the stem wood area per hectare (basal area). The
drag force is calibrated according to the NAIS document (“Nachhaltigkeit und Erfolgskontrolle im
Schutzwald, Frehner et al. 2005). This document describes the silvicultural management of protection
forests for different forest habitats and specific natural hazards. The revised requirement profile of a
rockfall protection forest includes numbers of trees and other information. Please consider that the
rock mass is essential for the protective effect of a forested area.
This simple but efficient forest representation applied in RAMMS::ROCKFALL will be further developed
in the future at SLF/WSL. User feedback or comments are very welcome. Please contact ramms@[Link]
for any ideas or suggestions.
60
CHAPTER 4: SETTING UP A SIMULATION
Figure 4-31: Project with emerging polygon shapefiles which represent forested areas.
• Create polygon shapefiles as described above for all the areas with different and
important surface covers inside the area of interest.
4.7.1 General
In this tab you specify the name of the output folder and the dump step (time interval between
different dumps saved to disk).
4.7.2 Terrain
Here you specify the overall terrain type and load the additional terrain shapefiles, specify their
terrain type and load them to the list. The available terrain types are:
61
CHAPTER 4: SETTING UP A SIMULATION
• Extra soft
• Soft
• Medium soft
• Medium
• Medium hard
• Hard
• Extra hard
• Snow
A detailed description of the terrain types is given in Table 3.1. Please make sure that the selected
shapes are depicted in the list otherwise they will not be considered in the simulation.
4.7.3 Forest
Here you load the forest shapfiles and specify the forest type and/or the type Lake/River/Moor to
stop rocks immediately. The available types are:
A detailed description of the terrain types is given in Table 3.1. Please make sure that the selected
shapes are depicted in the list otherwise they will not be considered in the simulation.
4.7.4 Release
Here you specify the random start orientations of the rocks and the release type.
To introduce variability into the rockfall simulation you have to specify the number of random initial
orientations of the rocks. This number multiplied with the number of rocks to calculate equals the
number of total simulations to perform. This number is automatically updated and displayed in the
Run ‘XXX’ Simulations button.
• Point
• Line/Multipoint
• Area
In the number of points field you can choose how many individual release points are nestled along
the line or polygon.
You can optionally set the rocks initial velocities and initial rotation velocities along all three axis (x,
y, z).
62
CHAPTER 4: SETTING UP A SIMULATION
To set the rock offset (initial fall height of the rocks measured from the center of mass) you have two
options. RAMMS can automatically calculate the minimal offset that is necessary to start the rock
(Automatic). Alternatively, you can set the rock offset manually. This offset should be high enough such
that the rock is not sticking in the terrain and cannot start.
4.7.5 Rock
In the Rock tab you can choose between three different types of rocks:
You can also choose a folder containing different rocks. You have to fill the folder with rocks generated
with the Rock Builder in advance. All rocks in the specified folder will be calculated.
63
CHAPTER 4: SETTING UP A SIMULATION
General Tab:
(4) Digital Elevation Model Information shows you which DEM is used for the simulation.
(5) Calculation Domain (optional): Narrow down your calculation area by using a polygon
shapefile.
(6) Activate the box Stop at first contact only if you wish to stop your simulation as soon as
the rock reaches the terrain for the first time.
64
CHAPTER 4: SETTING UP A SIMULATION
Terrain Tab
65
CHAPTER 4: SETTING UP A SIMULATION
1
Rock Tab → Cuboid / Sphere
2
(1) Choose rock type: Click Cuboid /
Sphere to run a RAMMS simulation
with a cuboid / sphere.
66
CHAPTER 4: SETTING UP A SIMULATION
67
CHAPTER 4: SETTING UP A SIMULATION
(3) The Nr of Grid Points show how many Figure 4-40: Area Release Tab
rocks will be simulated in minimum
(within one orientation). You can vary
the number of used grid points and
with it the number of simulated rocks
in minimum. How to select the
number of random orientations is
described in the Release tab for Point
release section.
Output Tab
In order to achieve statistical relevant results at the object of interest (e.g. village, road, endangered
object), after a first test with low numbers of released rocks and short calculation times, some main
model runs should be undertaken. The first main model runs should start with
and compare the new barrier plot for above polygon shapefile with the one from the first model run
with 10-15 Nr Random Orientations.
If the statistical results differ significantly, then increase your Nr Random Orientation and/or
Release Locations, until the statistical results for the barrier plot do not differ significantly anymore.
69
CHAPTER 4: SETTING UP A SIMULATION
If the scenario already exists, RAMMS will ask you if you want to overwrite the scenario:
If you click No or Cancel, you will be able to rename your scenario in the General tab. Click Yes, and the
old scenario is deleted (all the files and subdirectories within the scenario folder).
Scenario Preparation
If you specified terrain and/or forest shapefiles, RAMMS will gather the data and saving the
appropriate friction parameters ASCII files, showing messages like the one below in the lower left
status bar:
After that, RAMMS creates all the input files for this scenario…
70
CHAPTER 4: SETTING UP A SIMULATION
For every rock trajectory, RAMMS generates an output file (.rts). After the completion of all simulations
or if the user clicks Cancel, RAMMS will open the simulation files in Statistic Mode, see next chapter,
and present the Scenario Logfile:
71
CHAPTER 4: SETTING UP A SIMULATION
Additionally, a trajectory shapefile is created for every scenario, containing all the trajectories as
polyline shapefiles. Attributes of these polyline shapefiles are corresponding trajectory filename as
well as the rock name. The name of this shapefile is <scenario_name>_Trajectories.shp.
If some of the rocks could not release, because the specified Z-offset was not big enough, RAMMS will
show a message like this:
5 Results
RAMMS offers a variety of tools and visualizations to interpret the simulation results. There are two
different modes to open simulation results: a) Statistic Mode to visualize a large number of simulations
(> 100) and b) Trajectory Mode to analyze single trajectories in detail. The two modes are described in
this chapter.
Example:
We assume that RAMMS has saved 63 values (63 values of jump height, 63 values of velocity, etc.) in
a given cell, see figure below (showing only part of the 63 values):
73
CHAPTER 5: RESULTS
RAMMS then calculates the following statistic values for this given cell out of the 63 values (e.g. Jump
Height):
Use Help → Advanced… → Additional Preferences… → Edit or use the button (Additional
Preferences) in the lower left toolbar.
74
CHAPTER 5: RESULTS
……….
Find the keyword QUANTILE_VALUES, change the values accordingly and press the Save button.
Changing these values results in adjusted statistical analyses and dropdown menus, according to the
values the user entered (and saved!).
Figure 5-4: Adjusted quantile values in Statistic Mode; Left: Statistic Summary Plot, Right: Quantile
dropdown menu in upper right toolbar of GUI.
Wikipedia: “The probability density function (PDF), or density of a continuous random variable, is a
function that describes the relative likelihood for this random variable to take on a given value. The
probability of the random variable falling within a particular range of values is given by the integral of
this variable’s density over that range—that is, it is given by the area under the density function but
above the horizontal axis and between the lowest and greatest values of the range. The probability
density function is non-negative everywhere, and its integral over the entire space is equal to one.”
75
CHAPTER 5: RESULTS
Wikipedia: “In probability theory and statistics, the cumulative distribution function (CDF), or just
distribution function, describes the probability that a real-valued random variable X with a given
probability distribution will be found to have a value less than or equal to x.”
Wikipedia: “In statistics, the empirical distribution function is the distribution function associated with
the empirical measure of the sample. This cumulative distribution function is a step function that jumps
up by 1/n at each of the n data points. The empirical distribution function estimates the cumulative
distribution function underlying of the points in the sample and converges with probability 1 according
to the Glivenko–Cantelli theorem.”
Boxplot:
Wikipedia: “In descriptive statistics, a box plot or boxplot is a convenient way of graphically depicting
groups of numerical data through their quartiles. Box plots are non-parametric: they display variation
in samples of a statistical population without making any assumptions of the underlying statistical
distribution. The spacings between the different parts of the box indicate the degree of dispersion
(spread) and skewness in the data, and show outliers.”
Figure 5-5: Boxplot Explanations (Q1: lower quartile = 25 %, Q3: upper quartile = 75 %, IQR:
interquartile range).
76
CHAPTER 5: RESULTS
Open a simulation scenario with Track → Open… → Rockfall Scenario or click in the toolbar and
then choose a scenario in a project’s output folder.
Figure 5-7: The file name filter shows information about the scenario that you are opening.
You can enter a file string or click OK to open all files.
Enter a file name filter for the specific results or click OK, if you are interested in all the simulations
in the scenario folder. Click OK to proceed.
77
CHAPTER 5: RESULTS
If you would like to open several scenarios, click Yes to choose another scenario from the output
folder. It is important that the scenarios are saved in the same project folder. Click No if you want
to analyze only the results from one scenario.
RAMMS will open the scenario and show the 95%-Quantile of the kinetic Rock Energy (kJ). The 95%-
Quantile is the default selection of the quantile dropdown menu in the upper right toolbar. You can
change the default selection in the Additional Preferences (Keyword: QUANTILE, values between 0 - 5
correspond to the position of the quantile in the dropdown menu, see Figure 5-4).
The ROCKFALL tab on the right shows useful Statistic Mode information. The trajectory mode should
be OFF – you are in the statistic mode. You find general information about the selected scenario:
In Statistic Mode, only Jump Height (H), Rock Velocity (V), Resultant Rotational Rock Velocity and
78
CHAPTER 5: RESULTS
All other results are available in Trajectory Mode only. For the visualization of the parameters you can
choose between the values Mean, Median, 90%, 95%, 99% or Max values. The dropdown is located in
• Nr of Rocks
• Nr of Deposited Rocks
79
CHAPTER 5: RESULTS
• Barrier Plot
• Line Profile
• Quantile dropdown
menu
Choose the result parameter that you are interested in (e.g. Jump Height, Rock Velocity, Kinetic Rock
Energy or Resultant Rotational Rock Velocity). The results appear in the main window. Exercise 5.1
shows how to produce and analyze a Summary Plot.
80
CHAPTER 5: RESULTS
A Barrier Plot contains the same statistical information as a Summary Plot, but for a certain region of
interest (line or polygon). Additionally, for smaller datasets (< 10000 data values), an empirical
distribution function plot (EDF) is added to the CDF plot. For small datasets, the EDF plot (red line in
lower plot, see Figure 5-14) can deviate from the CDF plot. For large datasets, the EDF plot is not
shown, as the two plots are then identical.
It is important to note, that the analysis and interpretation of the barrier plot should be carefully
undertaken. The number of trajectories passing through certain grid cells (from a polyline or polygon
region) is crucial to establish a confidential statistic with reliable results. For this, a third axis has been
added to the line profile plot, the axis “Data Values”. This parameter indicates how the data values are
distributed along the line profile. For a rule of thumb, if data values drop far below numbers of 100,
then the statistical analysis becomes strongly dependent on the peaks (trajectories with very high
81
CHAPTER 5: RESULTS
kinetic energies and jumping heights). Deriving for example the Q95-value out of such strongly
asymmetrical distribution might result in an overrated, non-representative Q95 value. In the following
Figure 5-13, the line profile plot along a 87m long section is shown. The peak regions with jumping
heights of ~3m coincide to very low data values of 10 and less (right red circles). In contrast, where
data values exceed 200 or more (left red circles), the jumping heights are lower and are more
representative for this scenario.
Barrier plots are only possible for Jump Height, Velocity, Kin. Energy and Rot. Velocity results.
Draw a line profile by clicking (e.g. a dam or places that you are specifically interested in e.g. a place
where several trajectories pass through). RAMMS opens both a Line Profile Plot and the Barrier Plot
for the line:
Figure 5-13: Line Profile Plot depicting the Topography (green line, right y-axis), the jump height
(rock position) on top of the topography (red dots, Rocks), the Data Values distribution (blue
dashed line, left outer y-axis) and the Jump Height values (red line, left inner y-axis).
82
CHAPTER 5: RESULTS
Figure 5-14: Barrier Plot of a line profile. Scenario and line profile name are marked in the upper
right corner (red box). An empirical distribution plot (EDF, red line) is added to the CDF plot for
smaller datasets.
Use Statistics → Barrier Plot or the horizontal toolbar button to select a polyline or polygon
shapefile you wish to create a Barrier Plot for.
A window opens, displaying the Statistics Barrier Plot (see figure above). In the Data information part
of the Barrier Plot you find the information about the selected scenario and the name of either the
polyline or the polygon shapefile. Additionally you will find a statistical summary with the most
important statistic values, as well as information about the total number of trajectories passing the
region (3134) and stopping in the region (1), see Figure below (Traj./Stopped: 3134/1).
83
CHAPTER 5: RESULTS
Source:
The probability that a rock arrives in a given cell. The release cells (source) feeding a given cell are
taken into account when calculating the Source reach probability.
Total:
Same as above, but release cells are not considered; Total reach probability is calculated from total
released rocks.
Normal reach probabilities (Total and Source Reach Probability) are calculated on a grid cell basis. This
can lead to lower values than expected. With this feature it is possible to calculate a “Source Reach
Probability” for a polygon region.
84
CHAPTER 5: RESULTS
This feature can be used to analyze rock impacts on polygon areas, e.g. galleries, roads, dams, etc.
85
CHAPTER 5: RESULTS
• impact_data.txt and
• impact_trajectories.txt
impact_data.txt:
This file contains details of every trajectory impact (impact time and location, kin. Energy, velocities
and impact angle). The first entry “i” is the index-value into “impact_trajectories.txt”.
impact_trajectories.txt:
This file shows all trajectories passing your shapefile, even if they do not impact.
Use Extras → Open Trajectories from Shapefile to select a shapefile (polygon or polyline), and open all
the trajectories passing or stopping in the specified shapefile.
86
CHAPTER 5: RESULTS
Figure 5-19: Rockfall trajectories; different colors depict different kinetic rock energies (KJ).
87
CHAPTER 5: RESULTS
Open trajectories with the button , the keyboard combination Ctrl+T or go to Track → Open →
Rockfall Trajectories. Choose the trajectories that you are interested in from the output scenario
folder and click Open.
In Figure 5-21 you could also filter your selected trajectories. Otherwise click OK.
It is also possible to add trajectories from other scenarios (but within the same project). Click Yes in
Figure 5-22 to select more trajectories, click No otherwise.
88
CHAPTER 5: RESULTS
It is suggested to open only up to 100 single trajectories at one time (e.g. from a specific scenario).
Although it is possible to open more trajectories, it is not recommended, because the memory usage
will increase strongly and the handling of your visualization will get very slow. Control if the trajectory
mode is ON (in the general tab of the right field) and the toolbar shows the number of trajectories.
The buttons provide quick access to Jump Height, Velocity and Kinetic Rock Energy of the
simulation.
The drop-do wn menu Results offers the following functions in the trajectory mode:
(the unit of rotation can be adjusted in Additional Preferences, ROT_UNIT, default rotation
unit is rot/s)
• Total Rock Energy (kin + pot) (kJ)
• Kinetic Rock Energy (kJ)
• Kinetic Rock Energy (translational) (kJ)
• Kinetic Rock Energy (rotational) (kJ)
• Slippage (m)
• Friction Value ()
Figure 5-24: Results Velocity Figure 5-25: Results Jump Figure 5-26: Results Kinetic
Height Rock Energy
The values of Jump Height, Rock Velocity and Kinetic Rock Energy give a good overview of the
dimension of the rockfall event in the trajectory mode. You can choose the parameters in the
horizontal toolbar under Results.
90
CHAPTER 5: RESULTS
Information is refreshed
Figure 5-27: Trajectory information – General tab. The information (filename, start position of
trajectory, etc.) is refreshed in the General tab on the right side. Additionally, the Rock tab indicates
the rock used for the selected trajectory (see next figure).
91
CHAPTER 5: RESULTS
Rock Viewer
• Line Profile
92
CHAPTER 5: RESULTS
Click on a rock trajectory. You can find the function View Trajectory XY Plot in the horizontal toolbar
(button ) or in the menu Trajectory → View Trajectory XY Plot. If you want to analyze a trajectory
in detail, the XY plot of a trajectory is the way to go. The graph shows the currently active parameter
(you can change it by clicking one of the buttons in the upper horizontal toolbar or by
selecting another result parameter via the Results menu). You can keep the Trajectory XY Plot by
saving it, see Exercise 5.2b below ”How to create a Trajectory XY Plot”.
93
CHAPTER 5: RESULTS
▪ Red dots: points of contact (consider the time step of your simulation – missing contact
points are possible if time steps are too large).
▪ Bottom scale: projected profile distance (m).
Go to Extras → Profile... → Draw Line Profile or click to draw a line profile. The line profile
function provides a graph of the currently active parameter along a specific line through the rockfall
area. This is helpful when it is of interest to know the values and maximum values at these places
(e.g. close to a road or dam). Line profiles are saved in the file [Link] in the project directory. If
you want to keep a line profile, you have to save it, see Exercise 5.2c ”How to draw a line profile”.
• Activate the project by clicking on it once, then click or choose Extras → Profile…
→ Draw New Line Profile.
• Define the line profile in the same way you specify a new release line (see Exercise 4.4c
“How to create a new release line”). Finish the line profile with a right-click on the mouse
button.
94
CHAPTER 5: RESULTS
Plot explanations:
▪ Green line: track profile (altitude, y-axis on the right side).
▪ Red dots: active parameter (rock position) on top of the topography (altitude,
y-axis on the right side).
▪ Red line: active parameter (inner y-axis on the left side).
▪ Blue dashed line: the Data Values distribution.
▪ Bottom x-axis: projected profile distance (m).
If you change the active parameter or the Min and Max values in the Display tab
in RAMMS, the plot will be directly updated.
• To save the coordinates of the points belonging to the line profile, select Extras →
Profile… → Save Line Profile Points and enter a file name.
• To save the line profile parameter’s data (distance (m) and the active parameter,
e.g. the jump height (m)) at the current dump step, select Extras → Profile… → Export
Line Profile Plot Data and enter a file name.
• Activate the project by clicking on it once, then click or choose Extras → Profile…
→ Draw New Line Profile.
• Click the middle mouse button once.
• A window pops up and you can browse for the line profile you wish to open.
95
CHAPTER 5: RESULTS
To see the exact values of the simulation results, check the Trajectory Data Log File which shows
the results for every dump step of a single trajectory. After running a simulation click on a rock
trajectory in the trajectory mode and go to Trajectory → View Trajectory Data Log File to open it.
96
CHAPTER 5: RESULTS
Further information of a simulation run is available under Trajectory → Standard Output Log File
and
97
CHAPTER 5: RESULTS
The input file is a simple xml file and can be viewed from within RAMMS. Use Trajectory → View Input
File (xml) to open a selected trajectory’s input file.
98
CHAPTER 5: RESULTS
It is possible to start an animation of all trajectories. Switch to 3D mode by clicking (2D mode is
working as well).
Change the speed of the animation in the Display tab on the right side (speed slider from fast to slow).
It is possible to export your results as an image in different formats (e.g. .png, .jpg, .gif, .tif etc.). Click
or choose Track → Export… → Image File and define a file name with the corresponding extension.
An image of the visible part in the viewer will then be exported.
GIF animation
Creating a GIF animation is only possible in output mode. Click or choose Track → Export... → GIF
Animation. Enter a file name and location and wait until the simulation stopped. As soon as the
simulation finished, the GIF animation file is saved. In the Preferences you can define the interval for
the GIF animation (GIF animation interval (s)) in the Rockfall tab.
99
6 References and further reading
6.1 References
Maps and aerial images
➔ All topographic base maps and aerial images from Source: Swiss Federal Office of Topography
Manual Literature
Acary, V. and Brogliato, B., 2008: Numerical Methods for nonsmooth Dynamical Systems;
Applications in Mechanics and Electronics. In: Dynamics of Non-Smooth Systems. Lecture
Notes in Applied and Computational Mechanics, 35, Springer, Berlin.
Bourrier, F., Dorren, L., Nicot, F., Berger, F., and Darve, F., (2009): Toward objective rockfall
trajectory simulation using a stochastic impact model. Geomorphology, 110(3-4): 68-79.
Bourrier, F., Berger, F., Tardif, P., Dorren, L., and Hungr, O., 2012: Rockfall rebound: comparison
of detailed field experiments and alternative modeling approaches. Earth Surface Processes
and Landforms, 37(6): 656-665.
Cross, R., (1999): The bounce of a ball. American Journal of Physics, 67(3):222–227. [8] Fityus,
S., Giacomini, A., and Buzzi, O., (2013): The significance of geology for the morphology of
potentially unstable rocks. Engineering Geology, 162(0): 43-52.
Dorren, L., 2003: A review of rockfall mechanics and modeling approaches. Progress in Physical
Geography, 27(1):6987.
Fityus, S., Giacomini, A., and Buzzi, O. 2013: The significance of geology for the morphology of
potentially unstable rocks. Engineering Geology, 162:43 - 52.
Frehner, M., Wasser, B., Schwitter, R. 2005: Nachhaltigkeit und Erfolgskontrolle im Schutzwald.
Wegleitung für Pflegemassnahmen in Wäldern mit Schutzfunktion. Bern: Bundesamt für
Umwelt. 564p.
Glocker, Ch., 2001: Set-Valued Forces Laws. In: Dynamics of Non-Smooth Systems. Lecture
Notes in Applied and Computational Mechanics, 1, Springer, Berlin.
Glover, J., (2015): Rock-shape and its role in rockfall dynamics. Doctoral thesis, Durham
University.
100
Latham, J.-P.; Munjiza, A.; Garcia, X.; Xiang, J. and Guises, R., 2008: Three-dimensional particle
shape acquisition and use of shape library for DEM and FEM/DEM simulation. In: Minerals
Engineerings, Elsevier, 797-805.
Leine, R. I.; Schweizer, A.; Christen, M.; Glover, J.; Bartelt, P. and Gerber, W., 2014: Simulation
of rockfall trajectories with consideration of rock shape. Multibody System Dynamics 32,
2: 241 - 271.
Moreau, J.J., 1988: Unilateral contact and dry friction in finite freedom dynamics. In: Non-
Smooth Mechanics and Applications, Springer, Berlin, 1-82.
Schweizer, A. 2015: Ein nichtglattes mechanisches Modell für Steinschlag. Dissertation ETH
Zürich Nr. 22559.
Volkwein, A., Schellenberg, K., Labiouse, V., Agliardi, F., Berger, F., Bourrier, F., Dorren, L. K. A.,
Gerber, W., and Jaboyedoff, M. 2011: Rockfall characterization and structural protection - a
review. Natural Hazards and Earth System Science, 11(9):2617 - 2651.
6.2 Publications
The development of RAMMS::ROCKFALL is based on scientific findings published in international
scientific journals. A list of the most important scientific publications about RAMMS::ROCKFALL
and its applications is given below (chronological order). This publication-list is also available on
our webpage: [Link]
Leine, R.I.; Capobianco, G.; Bartelt, P.; Christen, M.; Caviezel, A., 2021: Stability of rigid body
motion through an extended intermediate axis theorem: application to rockfall simulation.
Multibody System Dynamics, 52: 431-455. doi: 10.1007/s11044-021-09792-y
Caviezel, A.; Ringenbach, A.; Demmel, S.E.; Dinneen, C.E.; Krebs, N.; Bühler, Y.; Christen, M.;
Meyrat, G.; Stoffel, A.; Hafner, E.; Eberhard, L.A.; Von Rickenbach, D.; Simmler, K.; Mayer, P.;
Niklaus, P.S.; Birchler, T.; Aebi, T.; Cavigelli, L.; Schaffner, M.; ... Bartelt, P., 2021: The relevance of
rock shape over mass - implications for rockfall hazard assessments. Nature Communications, 12:
5546 (9 pp.). doi: 10.1038/s41467-021-25794-y
Lu, G.; Ringenbach, A.; Caviezel, A.; Sanchez, M.; Christen, M.; Bartelt, P., 2021: Mitigation effects
of trees on rockfall hazards: does rock shape matter?. Landslides, 18: 59-77. doi: 10.1007/s10346-
020-01418-2
101
Lu, G.; Caviezel, A.; Christen, M.; Demmel, S.E.; Ringenbach, A.; Bühler, Y.; Dinneen, C.E.; Gerber,
W.; Bartelt, P., 2019: Modelling rockfall impact with scarring in compactable soils. Landslides, 16:
2353-2367. doi: 10.1007/s10346-019-01238-z
Caviezel, A.; Demmel, S.E.; Ringenbach, A.; Bühler, Y.; Lu, G.; Christen, M.; Dinneen, C.E.;
Eberhard, L.A.; Von Rickenbach, D.; Bartelt, P., 2019: Reconstruction of four-dimensional rockfall
trajectories using remote sensing and rock-based accelerometers and gyroscopes. Earth Surface
Dynamics, 7, 1: 199-210. doi: 10.5194/esurf-7-199-2019
Caviezel, A.; Schaffner, M.; Cavigelli, L.; Niklaus, P.; Bühler, Y.; Bartelt, P.; Magno, M.; Benini, L.,
2018: Design and evaluation of a low-power sensor device for induced rockfall experiments. IEEE
Transactions on Instrumentation and Measurement, 67, 4: 767-779. doi:
10.1109/TIM.2017.2770799
Leine, R.I.; Schweizer, A.; Christen, M.; Glover, J.; Bartelt, P.; Gerber, W., 2014: Simulation of
rockfall trajectories with consideration of rock shape. Multibody System Dynamics, 32, 2: 241-
271. doi: 10.1007/s11044-013-9393-4
Leine, R.I.; Capobianco, G.; Bartelt, P.; Christen, M.; Caviezel, A., 2021: Stability of rigid body
motion through an extended intermediate axis theorem: application to rockfall simulation.
Multibody System Dynamics, 52: 431-455. doi: 10.1007/s11044-021-09792-y
Caviezel, A.; Ringenbach, A.; Demmel, S.E.; Dinneen, C.E.; Krebs, N.; Bühler, Y.; Christen, M.;
Meyrat, G.; Stoffel, A.; Hafner, E.; Eberhard, L.A.; Von Rickenbach, D.; Simmler, K.; Mayer, P.;
Niklaus, P.S.; Birchler, T.; Aebi, T.; Cavigelli, L.; Schaffner, M.; ... Bartelt, P., 2021: The relevance of
rock shape over mass - implications for rockfall hazard assessments. Nature Communications, 12:
5546 (9 pp.). doi: 10.1038/s41467-021-25794-y
Lu, G.; Ringenbach, A.; Caviezel, A.; Sanchez, M.; Christen, M.; Bartelt, P., 2021: Mitigation effects
of trees on rockfall hazards: does rock shape matter?. Landslides, 18: 59-77. doi: 10.1007/s10346-
020-01418-2
Sanchez, M.A.; Caviezel, A., 2020: Full-scale testing of rockfall nets in real terrain. Results of tests
at Chant Sura: 13th September and 4th October, 2019. WSL Berichte, 97. 81 p.
Ringenbach, A.; Caviezel, A.; Demmel, S.E.; Lu, G.; Bühler, Y.; Christen, M.; Bartelt, P.; Meier, L.,
2020: Three-dimensional trajectory reconstruction of induced single block rockfall experiments.
In: Da Fontoura, S.A.B.; Rocca, R.J.; Pavón Mendoza, J.F. (eds), 2020: Rock Mechanics for Natural
Resources and Infrastructure Development - Full Papers. Proceedings. 14th international concress
on rock mechanics and rock engineering (ISRM 2019), Foz do Iguassu, Brazil. 2887-2894.
Wendeler, C.; Lanter, A.; Lu, G.; Caviezel, A.; Ringenbach, A.; Bartelt, P., 2020: New rockfall testing
method of flexible rockfall barriers. In: Da Fontoura, S.A.B.; Rocca, R.J.; Pavón Mendoza, J.F. (eds),
2020: Rock Mechanics for Natural Resources and Infrastructure Development - Full Papers.
Proceedings. 14th international concress on rock mechanics and rock engineering (ISRM 2019),
Foz do Iguassu, Brazil. 3416-3423.
Caviezel, A.; Lu, G.; Demmel, S.E.; Ringenbach, A.; Bühler, Y.; Christen, M.; Bartelt, P., 2019:
RAMMS::ROCKFALL - a modern 3-dimensional simulation tool calibrated on real world data. In:
102
2019: 53rd US rock mechanics/geomechanics symposium. 53rd US rock mechanics/geomechanics
symposium, New York, USA. ARMA 19-223 (8 pp.).
Caviezel, A.; Demmel, S.E.; Ringenbach, A.; Bühler, Y.; Lu, G.; Christen, M.; Dinneen, C.E.;
Eberhard, L.A.; Von Rickenbach, D.; Bartelt, P., 2019: Reconstruction of four-dimensional rockfall
trajectories using remote sensing and rock-based accelerometers and gyroscopes. Earth Surface
Dynamics, 7, 1: 199-210. doi: 10.5194/esurf-7-199-2019
Lu, G.; Caviezel, A.; Christen, M.; Demmel, S.E.; Ringenbach, A.; Bühler, Y.; Dinneen, C.E.; Gerber,
W.; Bartelt, P., 2019: Modelling rockfall impact with scarring in compactable soils. Landslides, 16:
2353-2367. doi: 10.1007/s10346-019-01238-z
Caviezel, A.; Gerber, W., 2018: Brief communication: measuring rock decelerations and rotation
changes during short-duration ground impacts. Natural Hazards and Earth System Science, 18, 11:
3145-3151. doi: 10.5194/nhess-18-3145-2018
Lu, G.; Caviezel, A.; Christen, M.; Ringenbach, A.; Meyrat, G.; Bartelt, P., 2019: Non-smooth
mechanics modelling of rock-tree and rock-forest interactions. In: Papadrakakis, M.; Fragiadakis,
M. (eds), 2019: COMPDYN 2019. 7th international conference on computational methods in
structural dynamics and earthquake engineering. Proceedings. COMPDYN 2019. 7th ECCOMAS
thematic conference on computational methods in structural dynamics and earthquake
engineering, Crete, Greece. 1612-1629.
Caviezel, A.; Schaffner, M.; Cavigelli, L.; Niklaus, P.; Bühler, Y.; Bartelt, P.; Magno, M.; Benini, L.,
2018: Design and evaluation of a low-power sensor device for induced rockfall experiments. IEEE
Transactions on Instrumentation and Measurement, 67, 4: 767-779. doi:
10.1109/TIM.2017.2770799
Gerber, W.; Caviezel, A., 2018: Measurement and analysis of ground contacts during rockfall
events. In: Cardoso, A.S.; Borges, J.L.; Costa, P.A.; Gomes, A.T.; Marques, J.C.; Vieira, C.S. (eds),
2018: Numerical methods in geotechnical engineering IX. 9th European conference on numerical
methods in geotechnical engineering (NUMGE 2018), Porto. 149-154.
Caviezel, A.; Bühler, Y.; Lu, G.; Christen, M.; Bartelt, P., 2018: Experimental validation of numerical
rockfall trajectory models. In: Cardoso, A.S.; Borges, J.L.; Costa, P.A.; Gomes, A.T.; Marques, J.C.;
Vieira, C.S. (eds), 2018: Numerical methods in geotechnical engineering IX. 9th European
conference on numerical methods in geotechnical engineering (NUMGE 2018), Porto. 875-883.
Eichenberger, V.; McArdell, B.; Christen, M.; Trappmann, D.; Stoffel, M., 2017: Wenn
Baumwunden dazu beitragen, Steinschlagmodelle weiterzuentwickeln. Schweizerische Zeitschrift
für Forstwesen, 168, 2: 84-91. doi: 10.3188/szf.2017.0084
Niklaus, P.; Birchler, T.; Aebi, T.; Schaffner, M.; Cavigelli, L.; Caviezel, A.; Magno, M.; Benini, L.,
2017: StoneNode: a low-power sensor device for induced rockfall experiments. In: 2017: 2017
IEEE sensors applications symposium (SAS). 2017 IEEE sensors applications symposium (SAS),
Glassboro, NJ, USA. 16807856 (6 pp.). doi: 10.1109/SAS.2017.7894081
Gerber, W.; Caviezel, A., 2017: Diversity of the results from drop weight tests. In: Corominas, J.;
Moya, J.; Janeras, M. (eds), 2017: RocExs 2017. 6th interdisciplinary workshop on rockfall
protection. 6th interdisciplinary workshop on rockfall protection (RocExs 2017), Barcelona,
Catalonia, Spain. 63-66.
103
Caviezel, A.; Christen, M.; Bühler, Y.; Bartelt, P., 2017: Calibration methods for numerical rockfall
models based on experimental data. In: Corominas, J.; Moya, J.; Janeras, M. (eds), 2017: RocExs
2017. 6th interdisciplinary workshop on rockfall protection. 6th interdisciplinary workshop on
rockfall protection (RocExs 2017), Barcelona, Catalonia, Spain. 59-62.
Leine, R.I.; Schweizer, A.; Christen, M.; Glover, J.; Bartelt, P.; Gerber, W., 2014: Simulation of
rockfall trajectories with consideration of rock shape. Multibody System Dynamics, 32, 2: 241-
271. doi: 10.1007/s11044-013-9393-4
Glover, J.; Denk, M.; Bourrier, F.; Volkwein, A.; Gerber, W., 2012: Measuring the kinetic energy
dissipation effects of rock fall attenuating systems with video analysis. In: Koboltschnig, G.; Hübl,
J.; Braun, J. (eds), 2012: 12th congress INTERPRAEVENT 2012. 23 to 26 April 2012. Grenoble,
France. Conference proceedings "Protection of living space from natural hazards". 12th congress
INTERPRAEVENT 2012, Grenoble, France, April 23-26, 2012. 151-160.
104
List of Figures
107
List of tables
108