Formulas For The Calculation of Breakwaters: January 2001
Formulas For The Calculation of Breakwaters: January 2001
net/publication/281274780
CITATIONS READS
0 3,407
1 author:
Miguel Losada
University of Granada
389 PUBLICATIONS 4,496 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
Asistencia Técnica para la Realización del Estudio de la Evolución de la Duna de Valdevaqueros View project
All content following this page was uploaded by Miguel Losada on 26 August 2015.
Introduction
This article analyzes the research done by Iribarren and his colleagues on
rubble mound breakwaters. The first part is a description of the impact that his initial
1938 formula had on the calculation of these structures, and the following sections
examine the experiments that he carried out in the 50’s and 60’s, which were based
on this formula. Apart from his important contribution to our knowledge about the
calculation of rubble mound breakwaters, Iribarren, whose work was afterwards
continued by Suarez Bores, also laid the scientific and technical foundations for
maritime engineering in Spain. This article studies the validity of Iribarren’s
proposals, and evaluates his technical and scientific legacy .
Throughout history, maritime structures, such as harbors and ports, have
played a vital role in transportation and communication systems. In ancient times,
harbors were built in places that were naturally protected from the sea, such as bays,
coves, inlets, estuaries, etc. From the end of the 19th century on into the 20th, the
progressive increase in maritime traffic, as well as ship size meant that breakwaters
had to be built in deeper waters. This was an even greater necessity on the Iberian
Peninsula because of its orography, and the scarcity of rivers and navigable estuaries.
Initially, such maritime structures were always built in protected areas, and
their shape depended on the configuration of embayments and existing rocky areas.
In the majority of cases, the typology of the breakwater section was that of a vertical
wall constructed by means of masonry. With the passage of time, new construction
materials, such as concrete, permitted the construction of artificial ashlars, e. g. the
floating caisson. The calculation methods of this section were very rudimentary with
little or no theoretical basis. In the first half of the 20th Century several failures
occurred in vertical structures, which brought to light the need for an immediate
solution, one that would be both technically viable and economically feasible.
1
Iribarren’s first period as a maritime engineer
In 1929 Iribarren was appointed Director of the recently created Grupo de
Puertos de Guipúzcoa, (Marín Balda 1999). From the very beginning, there were three
elements that would have a profound influence on his professional life: (1) his
responsibility for the protection of human life, property, ships, and maritime
structures; (2) the vulnerability of vertical breakwaters against the action of the
waves; (3) the lack of available methods and technical devices with which to evaluate
the effect of the action of the sea on breakwaters.
Iribarren thus embarked upon his professional activity as a maritime engineer
with dedication, enthusiasm, intelligence, and realism. He made full use of his talent
for observation, as well as his knowledge of the principles of rational mechanics
(among other things). After studying structures designed by other engineers and
observing the devastating effect of the action of the sea on the harbors of Guipúzcoa,
the region of his birth, he realized that vertical breakwaters were not capable of
withstanding the action of the waves breaking against their walls. He thus conceived
the idea of constructing sloping breakwaters with granular elements, which would
force the waves to break. This was the premise upon which he based all of his work.
Initially, Iribarren protected the vertical breakwaters with a layer of natural
elements, such as quarry stones, or artificial ones such as parallelpipedical concrete
blocks. When he saw that this solution worked, he decided to focus his activity on
the design of rubble mound breakwaters, which when built with successive layers of
granular material, caused the waves to break, thus gradually dissipating the energy
of the waves. He devoted almost thirty years of his life to this activity. The resulting
procedure for the design and calculation of such breakwaters was used very
successfully in the construction of more than 10 kms of breakwaters in Spain as well
as many other countries.
However, Iribarren’s path was fraught with obstacles. The sea, of course, was
the principal obstacle involved, but there were others, such as the lack of
comprehension of certain colleagues, which did not facilitate his work. Nonetheless,
his publications and the structures that he designed are testimony to his success in
overcoming all of these difficulties.
2
The first formulas
In 1938 Iribarren wrote the article (or paper, according to the author), titled
“A formula for the calculation of rubble mound breakwaters”. In it he proposes a
procedure for the deduction of a formula that combines the resistance characteristics
of the breakwater, weight [P], and density of the pieces [ γ s , ] and angle of the
protective level [ α , ] with the characteristics and height of the largest possible wave
N ∗H 3 ∗d
P= (1)
(cosα − senα ) 3 ∗ (d − 1) 3
In this formula, d = γ s / γ w , y, γ w is the specific weight of the water, and the practical
3
become discouraged by the many failed attempts and errors, which are an inherent
part of any research.
It is evident that Iribarren possessed all of these qualities in abundance since
he was obliged to endure the harsh and even hostile context prevalent in post-Civil
War Spain for those who wished to engage in research.
4
Research: 1954-1964
The positive response of the breakwaters in Orio and San Juan de Luz, the
accurate evaluation of the damage to the Mustafá Breakwater in Algeria, as well as
the enthusiastic reception of the formula in Navigation Conferences, not to mention
his evident vocation for research, encouraged Iribarren to continue with his work.
However, he perceived that only by carrying out experiments based on his own
criteria would he be able to arrive at a definitive formula and truly evaluate its
coefficients. Finally, Iribarren with the help of Nogales, undertook the construction
of the Laboratorio de Puertos of the Centro de Estudios y Experimentación de Obras
Públicas. He built a wave tank and a wave basin near the Retiro1, in a building
contiguous to the former Civil Engineering School in Madrid. This became his
headquarters, so to speak, or the place where he carried out his research with the
help of Suárez Bores, Tejedor and Sánchez-Naverac, all of whom are cited in
Iribarren’s 1965 article.
5
• direction of approach, θ
3. Parameters of the dike
• angle of the slope, α
• weight of the pieces of the main layer, P
• specific weight of the pieces of the main layer, γ s
Iribarren number, I r
Most of the experiments were carried out with waves breaking on the slope.
In order to know beforehand the wave height and period, Iribarren adopted a
“threshold value” for the quotient of the slope of the breakwater and the wave
H
steepness, ( tan α / ). This criterion was presented in a research paper given at
L0
the Navigation Conference held in Lisbon in 1949. The quotient is known today as
the Iribarren number (often termed the “surf similarity parameter”), and is generally
6
The research carried out in recent years had shown that the distribution of the
energy flow of a wave train reaching a slope in The research carried out in recent
years had shown that the distribution of the incoming energy flux reaching a slope
into (1) dissipated energy (2) transmitted energy flux, and (3) reflected energy flux
depends on the value of I r . Iribarren affirmed that for a certain “threshold value” of
the parameter, the wave does not break on the rubble mound. Years later, the type of
breaker associated with the “threshold value” became known as collapsing type of
wave breaking. The stability of the pieces of the slope is critical for this type of
breaker, and under these conditions, the reflection coefficient as defined by the
reflection quotient of the wave heights of the reflected and incident trains is superior
to 0.40.
Given that the flow regime is turbulent and according to dimensional
analysis, the weight of the pieces of the main layer can be expressed by the formula:
H d
P = φ (α , , failure criterion, type of piece) ∗ ( )∗ H3 (3)
L0 (d − 1)3
The dependence of the weight, P, of the pieces of the main layer on its relative
density appears with the original mathematical structure of the 1938 formula. Years
later, Losada and Giménez Curto (1982) showed that with that structure it is possible
to obtain an accurate adjustment of the experimental data for pieces with a relative
density in the range, d>2.00, which is normal in practical engineering. Furthermore,
the formula establishes the dependence of P on the cube of the wave height. This
dependence is derived from the dimensional analysis, and Iribarren very wisely did
not try to modify it.
The main objective of the research carried out by Iribarren and his colleagues
in 1954-1964 was the study of the influence of the four factors included in the
functional, φ , on the stability of the pieces of the main layer. For this reason, he
N ∗H 3 ∗d
P= (4)
( f ∗ cosα − senα ) 3 ∗ (d − 1) 3
In his work during that period the following elements were studied:
• the type of wave breaking and the ascending and descending equilibrium
• the friction and the interlocking between pieces
• the evolution of the failure and the temporal sequence of wave heights
7
• the type of failure and its dependence on the connection between pieces
This research was so complicated because all four factors intervened in each
of the four elements. The solution to the problem, expressed as a formula with
known coefficients, necessarily required data that would permit the effects to be
separated from the causes. Iribarren knew that the only way that he could do this
was to carry out experiments with a very specific type of design. The desire to obtain
a wave flume became almost an obsession for him, a need which he fortunately was
able to satisfy.
In the 1938 formula, the difference between the sine and cosine of the angle of
the slope was used to evaluate the dependence of the weight of the piece of the angle
of the slope, supposing that the friction coefficient and the interlocking between the
pieces is approximately equal to the unit. In posterior derivations, Iribarren and
Nogales (1965) kept a friction coefficient, f, between the pieces which affected the
cosine function in the following way: f ∗ cosα − senα .
In the American formula, published in 1953, 1959 and 1961, along with the
experimental results, this dependence was expressed by means of the function cotan
α , without the friction coefficient and interlocking between the pieces. The
variability of the experimental results was expressed by a sole coefficient that
Hudson called K D . Iribarren, who was convinced that the totally theoretical
deduction of his formula was correct, tried to discover which of the two formulas
was the most accurate. In order to do this, he first proceeded to determine the friction
coefficients, f, corresponding to different forms of elements, quarry stones, artificial
blocks, and tetrapods. The formula used by Iribarren is both intelligent and precise
since it reflects the mentality of an engineer who asks himself which formula leads to
fewer errors and provides more information about the task to be carried out.
To evaluate the friction coefficient, Iribarren devised the slopemeter with
which he could measure the angle at which a slope breaks in function with the
number and shape of the pieces. Thanks to a series of carefully planned experiments,
Iribarren discovered that the friction coefficient increases as the number of pieces of
the slope decreases, and that only when the slope consists of an important number of
edges, does the friction coefficient approach the unit. Moreover, his experiments
confirmed that the pieces with a tendency to rub against and interlock with each
8
other (i.e. tetrapods), have a larger coefficient than the pieces that only rub against
each other (i.e. quarry blocks).
The experimental results obtained with the slopemeter explained and
quantified why tetrapods are more resistant to waves than quarry blocks of the same
weight. However, with these results there arose the question of the number of pieces
that made up the “unit of response” of the slope in the face of the action of the waves
(in the words of Iribarren, the active or resistant slope). As will be discussed in one of
the following sections, Iribarren considered that this active slope was made up of six
equal sides of a cube, independently of the type of piece involved.
Some years after Iribarren first proposed this basic aspect of the stability of
the pieces, very slender pieces began being used, such as the dolo, whose hydraulic
stability was analyzed in the laboratory by applying the coefficient K D and ignoring
the mechanism by which the piece resisted the action of the waves.
The N coefficient
Once the value of the friction coefficient was resolved, the results of the
experiments in the wave flume were used to evaluate the N coefficient, which in
Iribarren’s proposal should only depend on the type of piece, and whether the
dominant water flow over the slope is ascending or descending.
Firstly, he observed that after various series of waves of increasing height
during which the structure remained intact when a certain wave height was reached
(which Iribarren called the breaking point), the accumulated curve of the number of
edges separated from the layer tended to increase with the wave height. The slope
thus would become stabilized after a certain number of series of waves of the height
that had been tested.
When the wave would reach a certain height, known as total breakage, the
damage would extend with increased velocity, finally affecting the layer at the depth
of an edge, or more specifically, of the equivalent side of the cube. In this way the
first layer became totally damaged, and the stability of the breakwater, seriously
affected. According to Iribarren, the average accumulated curve of loosened pieces
under conditions of partial stability was practically identical for quarry blocks,
artificial blocks, and tetrapods.
The relocation in the slope of the pieces detached from the damaged area
depended on the characteristics of the water flow over the slope. Accordingly, the
9
deposit for steeper slopes forms a berm underneath the damaged area. Just the
opposite occurs in the case of mild slopes, in which the berm builds upwards the
failure area. By analyzing the profiles of mild slopes under total damage, it was
possible to obtain the approximate location of the transition point between both
types of behavior. As it turned out, this depended on the type of piece involved. The
behavior associated with the formation of the berm below the damaged zone was
called descending equilibrium, since the pieces became detached when the water
mass was descending the slope. Iribarren was able to define with great precision the
worst conditions of descending equilibrium. According to him, such conditions occur
when the water flow falls freely, and is not slowed by the lowering of the contact
level with the slope, originated by the oscillatory movement.
Iribarren and Nogales (1964) considered that the dikes should be designed so
that for the calculation height, no damage would be produced. Taking into account
that the relation of the wave heights of total damage (100%) and the initiation of
damage (0%) is around 1.6, the relation between the weights of the elements or the
coefficients N 100% and N 0% of the formula should be 4.1. Therefore, these values of N
In order to reduce the volume of the quarry block, Iribarren proposed the
construction of a leveled parapet “...at low tide or a little above in order to be able to
work above the water (in dry conditions)...” and crowned “...at the height of 1.5*H
over the highest level of the sea at rest, which is that reached by the pressure
diagram of the breaker, although the almost vertical foam, which scarcely pushes,
can reach greater heights.” Furthermore, Iribarren recommended that the quarry
layer, “...ought to reach a height of 0.75*H, the same as the crest of the breaker...”.
To calculate the crown, Iribarren and Nogales (1964) proposed, “..to only use
the pressure diagram of the breaker very approximately and conveniently
modified…. and admit that the presence of the rubble mound reduces this pressure
10
by half”. They add, “It should be pointed out that the pressure diagram of the
breaker is only provisional and must be made more precise in the future..”.
Unfortunately, Iribarren died before he could finish the huge task that he had set
himself. He would have liked to have carried out new experiments to quantify the
behavior of the head, the effects of oblique wave incidence, the laws of pressure on
the crown, the influence of the randomness of the wave heights and periods, etc.
Years later after the 1965 publication, Bruun and Johannesson (1978) Bruun
and Gunbak (1978) for crowns, and Losada and Giménez Curto (1979) based on data
11
from Iribarren y Nogales (1965) for rip-rap, rubble mounds, blocks, and tetrapods,
confirmed that the worst stability conditions of a slope of these pieces is produced
when the wave breaks in collapsing. In this type of wave breaking, the time that
elapses between the two waves that consecutively arrive at the slope is practically
equal to the time taken by the water flow to ascend and descend the slope, there
being no interference between the two movements. This situation corresponds to the
worst case of free fall. When the wave breaks in plunging, the variation of the
average sea level affects the descent of the water down the slope, whereas when the
wave breaks in surging, the movement on the slope is basically governed by the
partially stationary oscillation that is produced on it.
The experiments always used the same core material, shape of the section
type and position in the canal, by which the reflected energy flow only depended on
the type of damage in the slope, which in turn depends on I r . Under those
López and Losada (1998) showed that for breakers with I r > 2.5 full plunging,
collapsing and surging and taking into account the reflection produced in the slope
to define the calculated wave height, the N coefficient of Iribarren’s formula is
h
constant for total damage and weakly depends on the relative depth , , at the foot
L
of the breakwater for the initiation of damage.
In 1953 and afterwards in 1960, Hedar, in the light of the 1938 formula,
deduced expressions to analyze ascending and descending equilibriums. On the
basis of this research, Iribarren established the worst conditions of descending
equilibrium of the pieces when the slope is rigid and the ebb of the water is free. He
also defined the conditions for which the worst change of equilibrium, whether
ascending or descending, is produced for each type of piece.
In recent years in certain countries where there are a lot of available quarry
stones close to the structure, deep water, and very high waves, breakwaters called
berms have been designed in which the stones can move with the action of the waves,
and the behavior of the section is similar to that of a beach with coarse-grained sand.
These breakwaters are built with gently inclined slopes in the section where the wave
12
breaks. In Spain such structures can be found in Tarragon and Alicante, among other
places. Iribarren was critical of the construction of breakwaters under conditions of
ascending equilibrium since the coefficient of the formula is greater than that of
descending equilibrium. In other words, all other parameters being equal, the piece
needs more weight.
13
their immediate objectives. It has often been discussed whether or not the calculation
wave proposed is adequate, in the light of what we now know about waves.
Iribarren determined the coefficient of the formula for the calculation of
wave height at undefined depths in function with the fetch, after collecting data
pertaining to the Cantabrian Sea. After all of these years and considering the fact that
Iribarren was a person accustomed to observing the sea, it can be assumed that the
calculation wave determined by the fetch formula, is an approximation of the
significant wave height, H s of a sea state in zones of depth limited wave heights. It is
thus possible to calculate breakwaters by applying the formula, as long as they are
found in depth limited wave heights, and the expected damage will be equal to or
less than the so -called total damage.
However, in areas without limitation of wave height by the depth or by the
generation process, the application of Iribarren’s formula gives optimistic results
since the arrival of larger waves is possible, whose relation with the significant wave
height is substantially greater than 1.8. With such waves arriving at the breakwater,
damage can occur that is greater than the so-called total damage. This important
topic was the object of a series of experimental and theoretical studies by Vidal,
Losada and Mansard (1995), who proposed as the calculation wave height a statistics
of order of the sea state.
Iribarren was aware of the occurrence of groups of waves, but to our
knowledge, never analyzed their influence on the stability of the pieces of the main
layer. Medina, Fassardi and Hudspeth (1990) and Medina, Hudspeth and Fassardi
(1994) focus on this problem. Their work has opened a line of research that would
surely have been of great interest to Iribarren.
A global analysis of all of three factors that affect breakwater stability was
proposed in the 1980’s by Suárez Bores, who like Iribarren, showed himself to be
ahead of his time by elaborating a multivariate method for the calculation of
maritime structures. Suárez Bores continues to work in this line of research to this
day. The results of his work can be found in the publication, Recommendations for
Maritime Structures, R.O.M. 0.0 (Losada 2001), as well as in research on the
application of neural networks on the calculation of maritime structures (Medina,
1999).
14
Oblique incidence and head stability
15
layer, which has a surprising level of accuracy. Losada and Giménez Curto (1981)
studied the run-down of the water flow of the slope, R d for different types of
Rd
damage and pieces. For collapsing breakers, the quotient < 1 , this limit
H
depending on the type of piece. Evidently, the maximum flow occurs slightly above
the point of maximum run-down. For this reason, when the result is lower than this
figure, it is not necessary to continue with the main layer.
Most of these criteria and standards are generally included in the technical
prescriptions for maritime structures built in Spain. However, certain engineers,
dazzled by the light coming from foreign shores, have occasionally built maritime
structures without following Iribarren’s criteria and standards. The result has been at
the very least, a loss of durability, and at the most, important failures in the structure.
Following Iribarren and Nogales’ book, it is possible to design a breakwater,
based on a geometry that previously defines its principal components, such as the
crown (ibid, p. 275). Paraphrasing the authors, this book is “... an important source
regarding such an important question as the design of breakwaters with a certain
technical basis, and without recurring to a deceptive comparison with already
existing structures...”. Although this contribution is mentioned somewhat less often,
it is doubtlessly one of the most important legacies that Iribarren left to maritime
engineering.
16
certain colleagues, who from their administrative ivory towers could not or would
not lend their support and assistance to the research efforts of two brilliant and
irreplaceable engineers.
Suárez Bores, recently retired from the Escuela de Ingenieros de Caminos,
Canales y Puertos in Madrid, has left an important group of followers, who now are
professors in Civil Engineering Schools and Faculties of Oceanography all over
Spain. In many of these academic institutions research has been and is presently
being carried out on the behavior of breakwaters (e.g. references in text).
17
us all, and with respect and admiration for the teacher as well as the man, we give
him our most heartfelt thanks.
References
Bruun, P. and Johannesson, P., 1976. “Parameters affecting stability of rubble
mounds”. J. Waterways, harbors, Coastal Eng. Div, ASCE 102 (WW2): 141-164
Bruun, P. and Gunbak, A.R.., 1978. “Stability of sloping structures in relation to
18
Iribarren, R. and C. Nogales, 1964. Obras Marítimas: Oleaje y Diques. Editorial Dossat,
S.A.
Iribarren, R., 1965. “Formule pour le calcul des diques en enrochements naturels ou
élements artificiels”. XXI Int. Navigation Congr., Stockholm, Section II-1.
Iribarren, R. and C. Nogales, 1965. “Fórmula para el cálculo de los diques de
escollera y elementos artificiales”. Revista de Obras Públicas, Madrid.
Larras, M.J., 1952. “L`équilibrie sous-marin d`unmasif de matériaux soumis á la
houle”. Le Genie Civil, September.
Losada, M.A. and Giménez Curto, 1979. “The joint effect of the wave height and
period on the stability of rubble mound breakwaters using Iribarren´s
Number”. Coastal Engineering, 3: 77-96.
Losada, M.A. and Giménez Curto, 1981. “Flow characteristics on rough, permeable
slopes under wave action”. Coastal Engineering, 4: 187-206.
Losada, M.A. and Giménez Curto, 1982. “Mound breakwaters under oblique wave
attack: a working hipótesis”. Coastal Engineering, 4: 187-206.
Marín Balda, J.L., 1999. La costa guipuzcoana, punto de partida de la personalidad
ingenieril de D. Ramón Iribarren Cavanilles. Real Sociedad Vascongada de los
Amigos del País. Publisher: Delegación en Corte, Dpt. de Publicaciones.
Martín, F.L., Losada, M.A. and R. Medina, 1999. “Wave loads on rubble mound
breakwater crown walls”. Coastal Engineering, (37), Nº 2: 147-174
Medina, J.R., Fassardi, C. and Hudspeth, R.T., 1990. “Effects of wave groups on the
stability of rubble mound breakwaters”. Proc. 22 nd ICCE, ASCE: 1552-1563.
Medina, J.R., Hudspeth, R.T., and Fassardi, C., 1994. “Breakwater armor damage due
to wave groups”. J. Waterway, Ports, Coastal and Ocean Eng. Vol 120, Nº.2:
179-198.
Medina, J,R., 1999. “Neural network modeling of runup and overtopping”. Proc. Of
the Int. Conf. Coastal Structures´99, Santander, Spain, Ed. I. J. Losada,
Balkema: 421-429.
Recomendaciones de Obras Marítimas, 1999. Ponencia de la ROM-0.1, (author Miguel A.
Losada). Puertos del Estado. Madrid
Suárez Bores, P., 1985. “Método multivariado para el cálculo de las Obras
Marítimas”. Class notes. E.T.S. de Ingenieros de Caminos, Canales y
Puertos. U. Politécnica de Madrid.
19
Vidal, C., Losada, M.A. and R. Medina, 1991. “Stability of mound breakwater´s head
and trunk”. J. of Waterway, Ports and Coastal Eng. ASCE 117, Nº 6: 570-586
Vidal, C., Losada, M.A. and Mansard, P.D., 1995. “Suitable wave-height parameter
for characterizing breakwater stability”. J. of Waterway, Port, Coastal and
Ocean Eng. ASCE, Vol. 121 (2): 88-97
20