Optimal Design of Block Quay Walls

Block walls, consisting of stacked unreinforced prefabricated concrete blocks, are commonly used for the construction of quay walls in the presence of rocky subgrades. A traditional design of block quay walls is based on manual design iterations, envisaging sufficient safety against ultimate limit states (ULS) such as sliding, overturning or loss of bearing capacity of the foundation soil. In addition, the designer should consider stability during the different construction stages of the block wall, referred to as construction constraints. This design process can be laborious, while the resulting designs comprise a large volume of concrete. In order to optimize block quay walls, we developed an automated design procedure in the framework of gradient-based optimization, accounting for the various ULS and construction constraints encountered in engineering practice. The design checks for a block quay wall are first explained in detail. This includes global ULS requirements that apply to the block wall as a whole, and internal ULS requirements to consider sliding and overturning of separate blocks. During all construction stages, the block wall has to be stable, which imposes additional construction constraints. Block walls consisting of rectangular blocks and chamfered blocks are optimized. The resulting designs obtained with the automated design procedure satisfy all design requirements, and have a realistic layout. Furthermore, the influence of the different construction stages is studied, demonstrating the practicality of the proposed automated design procedure.


INTRODUCTION
Block walls are gravity retaining walls that consist of unreinforced prefabricated concrete blocks. While block walls use a lot of material, they are commonly used as a retaining wall in the case of rocky subgrades that are characterized by a large bearing capacity, and where soil penetrating constructions such as combi-walls, sheet piling, or slurry walls are difficult or impossible to install (de Gijt and Broeken, 2013). Block quay walls reach typical retaining heights of up to 30 m and can be both constructed in dry conditions, or by placing blocks under water. The blocks are commonly prefabricated near the construction site, and stacked using construction cranes. Block quay walls are characterized by permeable joints, providing drainage that results in little or no water pressure differences over the wall.
A traditional design of block quay walls is based on manual iterations, governed by ultimate limit states (ULS) that require sufficient safety against failure mechanisms such as sliding, overturning or loss of bearing capacity of the foundation soil (de Gijt and Broeken, 2013). This traditional design approach, with manual iterations based on engineering judgement, is laborious and time consuming. Therefore, the question arises whether numerical optimization techniques could yield similar or better, more economic designs, at the same time reduce engineering costs.
Numerical design optimization has already been applied to a variety of geotechnical applications, including the design of piled foundations (Chan et al., 2009), underground excavations (Ren et al., 2014), and rock bolt design (Nguyen et al., 2015), mostly based on genetic algorithms and simplified design methods. For example, for cantilever retaining walls, Camp and Akin (2012) and Kaveh et al. (2011) employ a Big Bang-Big Crunch optimization algorithm, Khajehzadeh et al. (2013), Moayyeri et al. (2019), and Talatahari et al. (2012) use a gravitational search in the context of particle swarm optimization, whereas Gandomi et al. (2016) and Nandha Kumar and Suribabu (2017) apply evolutionary algorithms. While such genetic algorithms are easy to understand and implement, they can only cope with a limited number of design variables, such as the geometry of a cantilever wall that is described by a limited number of parameters. Furthermore, genetic algorithms require many iterations to converge. As the number of design variables increases, a gradient-based optimization algorithm becomes more attractive, while requiring careful treatment of the objective function and constraints. For the particular case of block quay walls, such a numerical optimization has been applied by Shafieefar and Mirjalili (2013) and Shafieefar and Mirjalili (2014), using a gradient-based sequential quadratic programming algorithm. This algorithm is particularly suited to cope with the large design freedom, varying the number, size and position of the blocks.
The typical loads for block quay walls include the weight of the blocks, the lateral earth pressures on the blocks, water pressures, surcharge loading, and loading by mooring ships. As design requirements, sufficient safety against sliding, overturning or loss of bearing capacity of the foundation soil are imposed. In addition, seismic loads are to be considered in seismic prone areas (Richards Jr. and Elms, 1979). A classical quasi-static seismic design commonly used in practice is based on the use of seismic earth pressure coefficients obtained from a Mononobe-Okabe analysis. Such a quasi-static analysis represents a strong simplification, and a lot of research is devoted to provide improvements, such as Newmark's sliding block method (Li et al., 2010;Caltabiano et al., 2012;Pain et al., 2017) or the use of dedicated models to account for the behavior of the interfaces between the blocks under complex loading scenarios, such as torsional effects (Casapulla and Maione, 2018;Casapulla et al., 2019). In the present paper, however, seismic loads are not considered.
Apart from these essential ULS stability requirements, the design of block walls also includes considerations on constructability: the block wall should be stable at every stage of the construction of the wall, and engineering a construction sequence forms an integral part of the design. Block walls are typically gradually backfilled a couple of times during construction for logistical reasons. The number of backfill operations and the corresponding stacking sequence hence provides additional design freedom, where the stability of the freestanding stack of blocks is considered as an additional constructability constraint. Other constructability constraints include e.g. a maximum block weight, governed by the cost and availability of construction cranes. This paper develops a gradient-based shape optimization algorithm for block walls, aiming at reducing the material use while accounting for realistic boundary conditions, load cases, and construction constraints. The width, height, and position of the blocks are adjusted in order to obtain an optimal layout. Accounting for ULS envisages sufficient safety against failure mechanisms such as sliding, overturning, or loss of bearing capacity of the foundation soil. Additional construction constraints include the maximum weight of the blocks and the stability of the block wall during construction. This constructability constrained optimization represents a first step toward a holistic design approach to quay walls (de Gijt, 2015), aiming at an optimal technical and economical design, and considering the total lifecycle costs that include the costs of engineering, construction, maintenance, demolition and environmental impacts.
This paper is organized as follows. In section 2, the design requirements are outlined by means of an example block wall, for which the different elements of the analysis are explained in detail. In section 3, the shape optimization problem is defined for a wall consisting of rectangular blocks, where a two-block wall allows to visualize the optimization process. Furthermore, the results of an optimization of a realistic wall with more blocks are discussed. Since chamfered blocks are commonly used in practice and allow for more design freedom, section 4 extends the optimization problem to such chamfered blocks, and the optimization results are compared with the case of rectangular blocks. Section 5 summarizes the main findings of the paper.

Problem Outline
In order to demonstrate the design principles of a block wall, a simple layout with four rectangular blocks is considered, as shown on Figure 1. The base of the wall is located at z = −6 m, and every block i = 1, . . . , n block , n block = 4 has a height t i = 2 m so that the total wall height H tot = 8 m. Furthermore, the water level w l = −0.5 m on the left side is lower than the water level w r = 0 m on the right side, resulting in an unfavorable water pressure on the wall. In the following sections, the material properties of the blocks, the backfill, and the foundation soil are listed. The forces acting on the structure are subsequently discussed, such as variable loads, effective soil pressures, and water pressures. This allows to check ULS requirements, including sliding failure, overturning, and loss of bearing capacity of the soil. Global ULS requirements apply for the block wall as a whole, whereas internal ULS requirements consider sliding and overturning for each block separately.

Material Properties
The concrete blocks are considered to be rigid and have a unit weight γ c = 24 kN/m 3 . The blocks have a compressive strength f cu = 45 MPa, which limits the stresses at the interface between the blocks. The foundation soil is characterized by an internal friction angle φ = 28 • , a cohesion c = 83 kPa, and dry and saturated unit weights γ dry = 20 kN/m 3 and γ sat = 22 kN/m 3 , respectively. These soil properties are representative  for rocky soils. The backfill is characterized by an internal friction angle φ = 40 • , a zero cohesion c = 0 kPa, and dry and wet unit weights γ dry = 20 kN/m 3 and γ sat = 22 kN/m 3 , respectively. Furthermore, the strength of the soil-block interface is characterized by an interface friction angle δ = 2/3φ.

Variable Loads and Load Combinations
The variable loads include the point loads Q x and Q z , distributed vertical loads q z on top of the wall, and surcharge loading q 0 next to the wall (Figure 1). Four load cases (LCs) are considered, denoted as LC 1 to LC 4 , as listed in Table 1. Sliding failure of the wall dominates LC 2 , in absence of the vertically distributed load and vertical point load on top of the wall. Overturning of the wall governs LC 4 , whereas the loss of load bearing capacity of the soil is dominant for LC 3 , with all loads considered. A negative Q x and q 0 are always unfavorable, and are considered in every LC.

Soil Pressures
The backfill generates soil stresses that represent the predominant loading on the block wall. Since the block wall has a denticulated back, the calculation of soil pressures is more complex than for straight retaining walls such as diaphragm walls. The soil pressures on block walls are computed by considering a vertical fictitious back plane (Figure 1), as described in the British Standard (BSi, 2015). In a manual design calculation, the fictitious back plane is determined on the basis of engineering judgement, roughly following the back of the block wall. In order to compute soil stresses in a gradient-based optimization scheme, an automated calculation is required. Therefore, the convex hull of the ensemble of blocks is considered as the fictitious back plane, since its sound mathematical definition allows for an algorithmic treatment. The fictitious back planes for the example block wall, computed from the convex hull, are shown in Figure 2. A different fictitious back plane is considered for the equilibrium of each block. On the fictitious back plane, the coefficient of active soil pressure is computed using the method of Müller-Breslau (1906), assuming a straight slip surface of the failure wedge behind the wall: where α is the slope of the fictitious back plane, β the slope of the ground surface (considered zero in this paper), and δ is the wall friction. The active soil pressure coefficients are listed for each block in Table 2.
While the wall can have a small embedment depth in which (stabilizing) passive soil pressures could be generated, these passive soil pressures are neglected since the soil near the front of the block wall is often eroded or stirred as a result of the construction work.

Weight of the Enclosed Soil
Apart from the weight of the blocks, the weight of the soil above protruding blocks (Figure 3) provides a stabilizing effect on the sliding and overturning failure mechanisms. The net weight of this enclosed soil is accounted for in order to avoid overconservative designs. Furthermore, buoyancy results in a weight difference of the soil above and below the water table.

Water Pressures
The water pressure acting on the example wall is shown in Figure 4, where a distinction is made between the water pressure acting on the block wall as a whole (for global equilibrium) and the water pressure acting on individual blocks (for internal equilibrium). The water pressure is linearly interpolated along the bottom and top surfaces, accounting for water flow through the permeable joints.

Global ULS Requirements
The global ULS requirements include sliding and overturning of the block wall, as well as loss of bearing capacity of the soil. The definition of the corresponding safety factors for each of these failure mechanisms is discussed next.   The global safety factor of the sliding failure mechanism is defined as: where V res and H res are the resulting vertical and horizontal forces, respectively, and f cg = 0.6 is the Coulomb friction coefficient between the bottom block and the foundation. The global safety factor for overturning: accounts for the resulting stabilizing moments M stab and the resulting destabilizing moments M destab around the toe of the block wall. The foundation has a width B, corresponding to the width of the first block. The resulting vertical force V res and horizontal force H res acting on the foundation are calculated from the sliding equilibrium. In order to compute the eccentricity of the force on the foundation, the moment relative to the center of the foundation M center is required. The moment around the toe of the wall M ′ = M stab + M destab . The moment around the center of the foundation M center is computed as: The eccentricity e of the load is computed as: The effective foundation width equals B ′ = B − 2e and the pressure on the foundation q eff is calculated as: The bearing capacity q ult is computed using Terzaghi's equation (Tomlinson and Boorman, 2001): where N c , N q , and N γ are bearing capacity coefficients, represented by non-linear functions of the angle of internal friction φ. In the present paper, the equations of Meyerhof are applied (Tomlinson and Boorman, 2001): The parameters a and b are 0.1054 and 9.6, respectively, for a rough foundation. The depth factors d c , d q and d γ are equal to 1 because a surface foundation is considered. The form factors s c , s q , and s γ are also equal to 1 for the plane strain case at hand. The slope factors i c , i q , and i γ take into account the angle at which the resulting loads occur, and are calculated as: where q is equal to the surcharge next to the foundation. The effective width B ′ , the horizontal and vertical resulting forces H res and V res are already known. φ and c correspond to the internal friction angle and the cohesion of the foundation soil and the effective bulk weight equals γ ′ = γ sat −γ w below the groundwater level and γ ′ = γ dry − γ w above the groundwater level. Applying Terzaghi's formula to the example block wall, the second term in Equation (7) is omitted because there is no effective vertical stress q ′ on the foundation. The global safety factor for the bearing capacity of the soil is defined as follows: The safety factors for the example block wall for all four load cases are shown in Table 3. It should be noted that sliding is the most critical failure mechanism. For sliding, LC 2 is critical since there are no vertical loads on top of the quay wall. Overturning is critical for LC 2 and LC 4 . The bearing capacity controls LC 3 , since all vertical loads are present. The vertical point load at the top of the wall has a major influence on the foundation pressure q eff for two reasons: (a) the eccentricity of the load considerably reduces the effective width B eff , and (b) the additional vertical forces cause higher stresses on the foundation. The increase in foundation pressure q eff is the main reason for the reduction of the safety factor for bearing capacity rather than the decreasing slope factors i c , i q , and i γ , as a consequence of the large eccentricity of the external load.

Internal ULS Requirements
The eccentricity of inter-block forces must be restricted so that a minimum contact area is subjected to pressure, avoiding decompression (loss of contact stresses) on one side of the contact surface, at the same time preventing excessive compressive stresses on the (small) remaining contact area. The compression ratio C compression,k is the ratio between the compressed part A ′ k and the total contact surface A k between two blocks: The index k refers to the upper block and varies from 2 to n block . In addition, the maximum contact stress σ max,k between the blocks is also restricted. The compression ratio and the maximum contact stress depend on the eccentricity e k of the vertical force V k on the contact surface k. If e k /B k < 1/6, the vertical contact force is located in the central core (middle third) of the contact surface, and decompression does not occur, i.e., C compression,k = 1. Correspondingly, the maximum contact stress equals: If e k /B k > 1/6, the vertical force is located outside the central core, resulting in decompression. The compression ratio is then calculated as: in which case the maximum contact stress σ max,k equals: where B k is the width of the contact surface. The eccentricity e k of the vertical force is calculated as e k = M center,k V k . The maximum contact stress of σ max,k is compared to the ultimate compressive strength of the blocks f cu by means of the ratio of F σ max, k : The safety factors for the internal ULS requirements are shown in Table 3 for each load combination. It should be noted that the safety factors for sliding and overturning of block 1 are equal to those of the global equilibrium (Table 3). Furthermore, the upper block 4 is the most sensitive to sliding due to a large horizontal point force on the top of the wall. Therefore, in practice, the upper block is typically casted continuously in-situ as a capping beam in order to distribute the forces over larger portions of the block wall. The overturning failure mechanism is critical for block 2. Decompression does not occur and the maximum stresses are limited. Therefore, concrete C35/45 with a compressive strength f cu = 45 MPa is more than sufficient to resist the resulting stresses.

Constraints on Protrusion
Blocks that excessively protrude above the underlying block are prone to break under their own weight or may overturn during construction. Therefore, two conditions are imposed for all but the first block (k = 2, . . . , n block ). The first condition prevents breakage: blocks may not protrude over a distance more than their height. This rule is also referred to as the 45 • rule because cracks in concrete are considered at an angle of 45 • . The condition can be expressed with a safety factor SF protrusion1,k : where t k is the thickness and U k is the protrusion length of the considered block. The safety factor is calculated for protrusions at the front as well as at the back of the quay wall. For the example block wall under consideration, only protrusion at the back is relevant. The second condition prevents the blocks from overturning during construction. The center of gravity of block k is enforced to be located within the contact surface with the supporting block k − 1. The corresponding safety factor SF protrusion2,k is defined as: where l k is the width of the supporting and L k is the width of the block.

Construction Constraints
The above safety factors guarantee sufficient safety with respect to ULS failure mechanisms. However, the global and internal ULS requirements assume that the block wall is supported by the backfill. For logistical and financial reasons, the block wall is only backfilled once or a couple of times during construction. As a result, the stability of the freestanding stack of blocks has to be guaranteed during the construction of the wall. If a number of blocks is stacked, overturning of the stack should be avoided. The most unfavorable situation occurs when the backfill is completed after stacking of the entire block wall. In order to account for constructability, safety factors are defined to allow to stack a number of blocks in each building stage without the need for gradual backfilling. If two blocks have been stacked, the center of gravity of the stack must be within the contact surface with the lower block, if three blocks have been stacked, the center of gravity of the stack must meet these conditions, and so on. In the following, we consider n stage construction stages where n blocki represents the number of blocks placed in stage i = 1 . . . n stage , and consider stability of the freestanding stacks of blocks. The safety factors SF constructability,m,n for the constructability of stage m are expressed as follows: SF constructability,m,n = L fic,k,n + 2z k,n L fic,k,n where z k,n is the center of gravity of the group of blocks of the considered check n. L fic,k,n is the fictitious length of the group of blocks of which the center of gravity z k,n is calculated. This fictitious length is obtained by dividing the area of the group of blocks under consideration by the total height of the group. The index n varies from 1 to (n BlockinStage,m × (n BlockinStage,m − 1)). The checks are performed for overturning to both sides. For the 4-block example, this corresponds to the six safety factors listed in Table 4.

Optimization Problem
The geometry of a block wall with rectangular blocks is uniquely defined by the number of blocks n block in the wall, the width L i and height t i , and the x-coordinate of the bottom left corner of each block x i . Therefore, the complete quay wall layout is described by the vector of design variables x: x = {L; t; T] = [L 1 L 2 L 3 ... L n ; x 2 x 3 ... x n−1 ; t 1 t 2 t 3 ... t n−1 } T (23) where x 1 and x n are fixed at x 1 = 0 and x n = 0. In addition, the height of the highest block t n = H tot − n−1 k=1 t k is not a design variable since the total height H tot of the wall is fixed a priori. As a result, the number of design variables equals n var = 3n Blocks − 3.
The total cost of a block wall depends on many factors including logistic considerations, location of the construction site, and availability of equipment and is therefore difficult to objectively determine. Therefore, the total volume of concrete is considered as the objective function f 0 to be minimized.
Most of the constraints to which the design has to comply correspond to the ULS safety of the block wall. The safety factors SF sliding,k , SF overturning,k , and SF bearing capacity are restricted in order to provide sufficient stability of the block wall, both considering global and local equilibrium. Furthermore, the compression ratio C compression,k is constrained, limiting the stress F σ max, k . In addition, cantilever equilibrium correspond to the factors SF protrusion1,k and SF protrusion2,k and the construction constraints in correspondence to the safety factor SF constructability,m,n should apply.
Furthermore, a weight limit M max per block is imposed. This limit is case specific and corresponds to the maximum capacity of the cranes available at the construction site. Furthermore, the contact area between the different blocks must have a minimum width B min in order to prevent local stress concentrations and failure. The minimum width of the contact surface is a constraint, in addition to the sliding constraint: sliding is only verified based on the ratio of forces and does not take into account the size of the contact surface. Finally, the relative weight difference between blocks is limited, avoiding large variations in block weight and reducing overall construction costs.

Optimization Algorithm
An appropriate algorithm is selected for the solution of the constrained optimization problem. Since the number of blocks in a block wall is rarely larger than about 15, and the forward calculation takes very little computation time, the sensitivities can easily be determined numerically by means of finite differences. Therefore, a gradient-based optimization algorithm is preferred. In this paper, the active set method, as implemented in the Matlab function fmincon, is employed. The active set method turned out to be the most performant algorithm provided by the fmincon function. The active set method typically works well for highly constrained problems, as in the present case.
The following convergence requirements have been considered: iterations are performed until a step tolerance f 0 = 10 −8 of the objective function is reached. A design is considered feasible if the constraints max(f ) is less than the assumed value 10 −4 , and the maximum number of iterations of the search algorithm is set to 70. The gradients of both the objective function as well as the constraints are computed using a finite difference approach.

Optimization of a Two-Block Wall
The optimization algorithm is first applied to a block wall consisting of 2 blocks, only considering the lengths L 1 and L 2 of the blocks as design variables. The objective function and the constraints can thus be visualized using a contour plot in (L 1 ,L 2 )-space. This visual representation allows to check whether the search algorithm actually converges to the desired optimum.
For various start values of the algorithm, the same optimum is found. In the optimum (L = {4.9615 m, 6.6154 m} T ), all design requirements are met. Figure 5 shows a contour plot of the objective function and all constraints. For sliding failure, LC2 is controlling the design. For the overturning failure mechanism, LC2 and LC4 are controlling. LC1 or LC4 control the loss of bearing capacity as well as the decompression, depending on the value of the design variables L 1 and L 2 . This example demonstrates the importance of taking all load cases into account.
The constraints for sliding, overturning, and maximum weight of block 2 are linear or constant and only depend on L 2 . Similarly, the protrusion, construction, and weight regularization constraints are linear in L 1 and L 2 . However, the constraints on sliding and overturning of block 1 depend on both L 1 and L 2 in non-linear way. This can be explained as follows: at a value of L 1 = 7.7 m for the sliding failure and L 1 = 5 m for overturning, the weight of block 2 is not required for stabilizing block 1. As a result, the constraint does not depend on L 2 . However, if L 1 becomes smaller, the constraints require an increasing length L 2 of block 2 to provide sufficient safety against sliding and overturning of block 1.
The contour lines of the constraints on decompression and bearing capacity are non-linear, due to the different conditions in calculating those safety factors. It should be noted that L 2 should not be less than 3.5 m to comply with the decompression requirements. If the length is greater than L 2 , the length L 1 must also be greater for the same reason. The bearing capacity requirement imposes a minimum length L 1 = 5 m. From length L 2 = 4.5 m, an increasing L 2 involves increasing L 1 in order to resist the extra weight of block 2.
At the bisector of the contour plot in Figure 5, a discontinuity occurs in constraints of the (a) sliding failure of block 1, (b) tilting of block 1, (c) loss of bearing capacity of the wall and (d) decompression of the contact area between block 1 and block 2. These discontinuities occur when L 1 = L 2 , and the friction angle on the fictitious back plane is replaced by an interface friction angle δ.

Optimization of a Wall With 12 Blocks
The algorithm is also used to optimize the design of a wall with 12 blocks and a retaining height of 30 m. The base level is located at z = −28 m. The soil properties, the water levels, the forces on the wall and at ground level, the load cases, the concrete strength and the friction coefficients are the same as for the example block wall.
The design variables are described in Equation (23), where the x-coordinates of the left corner of the lowest and highest block are fixed to 0. The thickness of the highest block is fixed in order to obtain a retaining height of 30 m. The resulting number of design variables is 33. To be realistic, the blocks are limited to a maximum weight of 80 tons. Since the blocks are typically 2 m long, the weight limit is set at 40 tons/m. For the lower block, however, the weight is increased to 80 tons/m in order to have a sufficient bearing capacity. The factor F weight is equal to 0.5 and the minimum width of the contact surface between the blocks B min is fixed at 4 m. The design variables are box constrained: the block width varies between 1 and 15 m, the x-coordinate between 0 and 5 m and the heights of the blocks range from 0.5 to 4 m. The height of the lower block, however, is limited to a minimum of 2 m because that block must be sufficiently thick so it will not break and provide sufficient load-bearing capacity.
The outcome of the active-set optimization algorithm strongly depends on the initial values of the design variables. Different initial values can lead to convergence to (sub)optimal solutions. Therefore, series of 100 simulations are performed with random initial values and the five best designs are retained. If the total volume of the five best designs are very close to each other, it is assumed that the optimum has been found. The initial values are randomly selected but (1) must meet the box constraints, (2) must have a positive initial value for the height of the upper block and (3) the contact surface between blocks must have a minimal width of 1 m.
Six different building scenarios for the wall are considered. The first scenario simulates the design of a 12-block wall in a single stage. Scenario 2 to 6 simulate series in which construction takes place in 2 stages of 6 blocks, 3 stages of 4 blocks, 4 stages of 3 blocks, 6 stages of 2 blocks, and 12 stages of a single block, respectively. The optimization problem has 33 design variables and 531 constraints (at maximum, when n Stage = n block ). Simulation times vary between 20 and 40 s on a regular PC. Since the optimization problem is non-convex, a series of 100 designs is obtained from random starting values, and the best designs are retained. Computing 100 designs takes 1 h. The total volume and wall layout of the best designs of the hundred simulations is shown in Figure 6. Figure 6 also shows how the volume evolves with the number of construction stages, ranging from 170.0 m 3 /m for a single stage to 153.9 m 3 /m for construction in 12 stages. Additional construction stages allow to further reduce the volume of the design. The results can be used to perform a cost-benefit analysis to determine an optimal number of construction stages: if replenishment of the backfill is relatively inexpensive, then a choice for 3, 4, or 6 construction stages is an excellent choice. If the cost of additional construction stages is larger than the reduced material cost, then building in 1 or 2 stages is recommended. A maximum of approximately 15.2 m 3 /m of concrete can be saved by building in 6 stages compared to building in 1 stage. The difference between building in 3 or 6 stages is small.
In the optimized designs, as shown in Figure 6, a relief platform can be observed at the level of blocks 8, 9 or 10. The relief platform is located at a position low enough to allow the stabilizing action of the superimposed mass of soil and high enough so that it is beneficial for the internal equilibrium of 8 to 10 of the 12 blocks. In the case of series 1, the platform is at a lower position near block 5. The reason for this is that with a higher platform, as in series 4, the structure is not stable during the construction of the wall in 1 stage. Therefore, the platform is located at a lower position. The blocks above that location are almost all large blocks that keep themselves in balance and provide the counterweight to balance the earth pressures. Series 2 is of an intermediate form, with a platform at the height of block 8. At the front of the quay wall, an arch shape can be seen in each case. The arch at the front can be explained by the effect of force around the toe of the wall: the more the blocks are located to the rear, the greater the lever arm of their weight and therefore the greater the stabilizing moment that prevents tilting.

Optimization Problem
In practice, the corners of the rectangular blocks are commonly chamfered, so that a smooth face is obtained. This section investigates to what extent the use of chamfered blocks can further optimize the design. Chamfered blocks are characterized by the lengths of the blocks L, the x-coordinates of the bottom left corner of the blocks x and the heights of the blocks t. The chamfers, considered at an angle of 45 • as shown in Figure 7, are 4 additional parameters, where the height of the chamfers is expressed as a percentage (d1 i to d4 i ) of the block height. The vectors d1, d2, d3, and d4 are added to the set of design variables: x = {L; t; T; d1; d2; d3; d4} = {L 1 , L 2 , L 3 , ..., L n ; x 2 , x 3 , ..., x n−1 ; t 1 , t 2 , t 3 , ..., t n−1 ; ... ..., d1 2 , d1 3 , ..., d1 n ; d2 1 , d2 2 , ..., d2 n−1 ; d3 1 , d3 2 , ..., d3 n−1 ; The chamfers d1 and d4 of the lowest block and d2 and d3 of the highest block are not considered, in accordance to actual engineering practice. The number of design variables n var increases to 7 × n block − 7 with respect to that of the rectangular block wall.
The constraints f 12 and f 13 impose that the two lower and the two upper chamfers of a block do not exceed the block length. The constraint f 14 ensures that the width of the top of the quay wall always remains greater than the minimum width B min , considered 4 m in this paper. The conditional constraints f 15 to f 18 have been introduced to avoid chamfers as shown in Figure 8A, since such a situation is never considered in practice where a layout as shown in Figure 8B is commonly preferred. Figure 8A shows that both parameters UH 1 and UH 2 , and the bevel d1 2 of the upper block are positive and therefore the limitation is not met. For all previous conditional constraints, a term F cond has been added to improve the convergence of the algorithm. The term F cond takes a small value and allows the constraints to be violated in a controlled way, where the tolerance is determined by trial-and-error. Constraint f 19 imposes that either bevel d1 or bevel d2 of a particular block must be equal to zero and limitation f 20 that either bevel d3 k or bevel d4 k of a particular block k must be equal to zero. This corresponds to a bevel either at the top or at the bottom corner. Correspondingly, the blocks can only take a trapezoidal shape, avoiding octagonal shapes.

Optimization of a Wall With 12 Blocks
In this section, a block wall consisting of 12 chamfered blocks is optimized. Apart from the chamfering, all input parameters and algorithmic parameters correspond to the 12 rectangular block case. Box constraints are also applied, where the length of the blocks varies between 1 and 15 m, the x-coordinate varies between 0 and 5 m and the height of the blocks varies from 0.5 to 4 m. The thickness of the lower block is limited to at least 2 m. The parameters relating to the chamfers may vary from 0 to a maximum of 70% in order to avoid sharp, fragile corners.
A series of 50 designs is simulated using random start values and the best designs are retained. The initial values meet the same conditions as specified in section 3. In addition, the initial values must meet the following condition to prevent chamfers from exceeding the length of that side: Per series, 50 simulations have been conducted. Each series contains different construction constraints, i.e., building the wall in 2 stages of 6 blocks, 3 stages of 4 blocks, 4 stages of 3 blocks, 6 stages of 2 blocks, and 12 stages of a single block, respectively. The factor F cond has been determined by trial-and-error to 0.025 for series 1 and series 3, and to 0.010 for series 2. The time to optimize the design is about 2 min, which is significantly longer than that for rectangular block walls. The results of different series are summarized in Figure 6. The resulting designs are smooth at the front and the back. In each design there is a platform, which has also been observed in the design of block walls from rectangular blocks, that exploits the stabilizing effect of the soil. The resulting designs herein obtained are similar to the block walls that are being designed in practice.
The volume of the block walls does not differ much with respect to the number of construction stages n Stage . Building in 2 stages is the most advantageous in terms of volume. The additional constraints, ensure that the construction constraints are only active to a limited extent and therefore have little influence on the total quay wall volume. The construction constraints therefore are less important in the general shape of the wall, in contrast to what was observed for the rectangular block walls.

CONCLUSION
In this paper, an automated design procedure for block quay walls taking into account construction constraints has been proposed in the framework of gradient-based optimization. Block walls consisting of rectangular or chamfered blocks have been considered, accounting for the various ULS and construction constraints encountered in engineering practice, aiming at minimal material use while maintaining a target factor of safety.
The design checks for a block quay wall have been explained in detail. This includes global ULS requirements that apply to the block wall as a whole, and internal ULS requirements that consider sliding and overturning of separate blocks. During construction, the stability of the block wall has to be guaranteed during all construction stages, which imposes additional construction constraints.
Block walls consisting of rectangular blocks and chamfered blocks have been optimized. The resulting designs obtained with the automated design procedure satisfy all design requirements, and have a realistic layout. The importance of the construction constraints has been demonstrated, as well as the corresponding additional volume of concrete in the wall. A number of typologies are identified, which depend on the construction constraints. In addition, the optimization of a block wall with chamfered blocks, as used in practice, results in very realistic designs, and allows for a fully automated design and optimization of the block walls. As such, the use of the proposed automated design procedure allows for a strong reduction of the engineering time, typically spent to manually optimize the design.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.

AUTHOR CONTRIBUTIONS
LL performed the implementation and analysis of the cases, where SF, MS, and HV were involved in the planning and supervision of the work. SF and LL prepared a first draft of the manuscript. All authors aided in interpreting the results and worked on the manuscript, discussed the results and commented on the manuscript.