Open access peer-reviewed chapter

Earthscraper: A Smart Solution for Developing Future Underground Cities

Written By

Faham Tahmasebinia, Kevin Yu, Jiachen Bao, George Chammoun, Edwin Chang, Samad Sepasgozar and Fernando Alonso Marroquin

Submitted: April 5th, 2019 Reviewed: June 3rd, 2019 Published: May 13th, 2020

DOI: 10.5772/intechopen.87217

Chapter metrics overview

625 Chapter Downloads

View Full Metrics


This chapter reports on the finite element analysis of the “earthscraper,” proposed by BNKR Arquitectura. It was proposed as an alternative building method for the future, as it requires less surface area and lower operating costs than an equivalent aboveground structure. A 2D model of the cross section of the structure was created using Strand7 for steady-state thermal analysis. This solver gave internal temperature ranging between 20 and 38°C between the bottom apex and the surface, respectively. This provides a comfortable temperature by default, displaying the lesser dependency on heating and cooling costs. A 3D model was also created to analyze the effect of lateral earth pressure by the use of the linear static solver. Results give a maximum lateral displacement of 527 mm and 19.8 mm on the exterior and interior walls, respectively. The model was used for earthquake analysis in accordance with AS/NZS1170.4, requiring the natural frequency and spectral response solvers. Twenty-five modal frequencies were found, with 99.6% of the mass of the structure contributing to the direction under analysis. The maximum horizontal displacement of the structure under the designed earthquake loads was found to be 19.2 mm.


  • earthscraper
  • future cities
  • underground development
  • numerical modelling
  • finite element analysis

1. Introduction

The United Nations has projected that the world population will reach 9.7 billion by 2050, with the increase of 2.5 billion as to date. The increased pressure on housing solutions to accommodate for the growing population while reducing the carbon footprint has led BNKR Arquitectura [1], based in Mexico City, to formulate a radical solution called the “earthscraper.”

Such a building is an inverted pyramid with the ability to accommodate 100,000 while using a small surface area, allowing easy implementation into developed cities. In addition, its founding into the soil eliminates any wind loading while also ensuring a constant warm temperature imposed by the surrounding soil, reducing the energy consumption for heating and cooling.

The projected populated rise in Sydney is forecast as between 2 million [2], while urban sprawl is limited by the surrounding Blue Mountains, countryside, and ocean [3] (Clarke, 2016). As such, the main objective of this project will investigate the feasibility of implementing the earthscraper in Sydney’s geology and climatology. The theoretical loads imposed by earth pressures, thermal loads, and earthquakes loads will be determined, and the structures response will be modeled using Strand7, a finite element method processor.

Post-processing of the results will give the displacement of the building under each respective load, providing an insight into the expected structural performance. The possibility of such a structure as an addition in the future cities of 2050 shall be determined in this chapter (Figure 1).

Figure 1.

Conceptual design of earthscraper.


2. Structural members

Below is the description of the members of the structure:

Floor system: reinforced concrete one-way floor slabs, 0.3 m thick with N16 reinforcement, sitting on steel beams and acting compositely.

Beams: steel beams of varying length extending from bottom of the exterior wall to the top of the interior wall, interconnected to form a truss system. Steel 360UB44.7 beams supporting the concrete slabs.

Columns: steel UB360.44.7 column spacings at 6.5 m, with floor-to-floor height of 5 m. Steel columns replaced by truss system at locations of intersection.

Walls: reinforced concrete with 1.5–2 m thickness and N20 reinforcement, acting compositely.

Foundation: concrete on rock-concrete pads for columns extruding into the concrete.

Cables: rock bolts to tie back structure.

Ground: water table at 60 m depth, soft soil of unit weight 17 kN/m3 from depths 0–10 m, and then Hawkesbury sandstone with unit weight 24 kN/m3 from depths 10–300 m.

Ceiling: perlucor glass ceiling with improved strength and heat insulation compared to normal glass.

Soil: spring dampeners on the exterior of the structure to simulate the soil and the structure’s deformations within the soil.

Table 1 shows the general element sizes used in the design of the structure. Element sizes were chosen based on design standards and structural strength, combined with numerical analysis to determine the most appropriate member size for the expected loads.

Details of the structural elementsSuggested structural element sizeSuggested design standards
Steel columns460UB82.1AS4100 [4]
Reinforced concrete floor system40 MPa, thickness 0.3 m, N20 rebarsAS3600 [5]
Reinforced concrete interior wall40 MPa, thickness 2 m, N20 rebarsAS3600 [5]
Reinforced concrete exterior wall50 MPa, thickness 2 m, N20 rebarsAS3600 [5]
Truss beams360UB50.7AS4100 [4]
Perlucor ceilingThickness 0.1 mAS1288 [6]

Table 1.

Structural elements.


3. Structural system

Figure 2a shows the dimensions of the structure, with a Perlucor glass ceiling and reinforced concrete walls and a reinforced concrete core running through the center of the structure. The green lines represent the outline of the interior concrete core, the black lines are the outline of the whole structure, the dark blue lines represent the floors, and the light blue lines on the top represent the Perlucor glass.

Figure 2.

(a) CAD rotated section view of entire structure and (b) CAD plan view from ground level.

Figure 3a and b shows the interior design of the structure for the first level, with the steel cables forming a truss system between the exterior wall and the interior core to provide buckling and bending resistance from the lateral earth pressure loads and earthquake loads.

Figure 3.

(a) CAD section view of first level, (b) CAD plan view of first level and (c) spring dampeners.

Figure 3c shows the spring dampeners with fixed restraints at the end away from the wall and free restraints at the connection with the wall to simulate the behavior of the structure within the soil mass.

The concrete floor slabs use N16 reinforcement bars at a spacing of 100 mm, while the reinforced concrete wall and core use N20 reinforcement bars at a spacing of 100 mm. The reinforcement helps the structure resist tensile loads due to the large bending moments caused by the lateral earth pressure.

The system resists vertical loads through load transfer from the roof and floor slabs into the columns, through the beams and truss system. The vertical loads within the columns are then transferred into the surrounding sandstone. The beams and columns of each individual floor have been omitted from the design, and the 60 floors of the 300 m deep structure have been simplified to 6 floors, as the vertical loads are not the focus of this structural design, since the surrounding Hawkesbury sandstone will be able to withstand the dead and live loading of the structure.

The system resists lateral loads with the truss system and floor slabs acting as struts, reducing the overall deflection of the structure and distributing lateral loads throughout the floor and into the surrounding bedrock. The floors resist loads by transferring vertical live and dead loads into the support steel beams, and the slabs resist lateral loads by transferring the loads either into the interior wall or exterior wall. The beams resist loads by transferring horizontal and vertical loads into nearby columns or walls. In addition, the beam truss elements provide different load paths and brace the structure to further reduce the stress concentrations on members. The columns resist load by transferring horizontal and vertical loads into the base pads and then into the surrounding bedrock.

Since the structure and its 50 floors are very complex to design for using Australian design codes, it was simplified into 6 larger floors to consider overall trends of stresses and deflections. This was done to locate areas of significant stresses and deflection and where the structure is most likely to fail. Strategies were then implemented to counteract these problem areas, such as increasing wall thickness or concrete strength.


4. Loads

4.1 Lateral earth loads (earth pressure)

The loads imposed onto the earthscraper by the surrounding soil and water were calculated using effective earth pressures, which are shown below and expanded in Appendix A:


Eq. (1) shows the relationship between vertical effective pressure and lateral earth pressure through the at-rest coefficient of earth pressure, K0. Eq. (2) gives an expression for K0, given by Jaky [7] as a function of the internal friction angle ϕ. These lateral earth pressures act on the exterior wall, causing horizontal deflection of external and internal wall after load transfer.

Figure 4 shows the lateral pressures varying linearly, with a change in gradient at 10 m depth with the change from soft soil to Hawkesbury sandstone. At a depth of 60 m is the water table, introducing pore water pressures and increasing the lateral loading on the structure. This representation of earth pressures.

Figure 4.

Vertical profile of lateral earth pressure—a combination of soil pressure and water pressure.

4.2 Thermal loads

The thermal loads imposed onto the structure include the solar heat flux, the convection currents between the structure and outside fluid, and the external ambient pressure of the surface. The ambient temperature surrounding the earthscraper was determined through research of past studies. The topsoil temperature was averaged across 1980–2006 for the top 20 cm, determined as 19.65°C [8]. The temperature of the soil is taken to linearly increase with depth, with a determined thermal gradient of 0.07°C/m.

The ambient temperature of the surface was chosen as the mean temperature at the Sydney Observatory Hill Bureau Weather Station since 1859, determined as 22.1°C [8]. However, due to the incoming solar radiation, the heat flux also contributes to the heat transfer to the Perlucor ceiling. Calculations for an absorbing blackbody give a heat flux of 1.37kW/m2; however, due to cloud albedo, atmospheric absorption, and radiation reflection, the solar exposure is taken as 69.44W/m2[8].

4.3 Earthquake loads

The response of the structure under earthquake loads requires the use of dynamic analysis, as the earthscraper falls into the earthquake design category III (EDCIII), as determined in Appendix C. This method shall comply with AS/NZS1170.4 [9] Section 7 with the use of a modal response spectrum analysis.


5. Structural design

The initial analysis was undertaken with all concrete wall elements and floors as the same plain concrete with 40 MPa strength and 0.5 m wall thickness. Upon linear static and spectral analysis, it was found that horizontal deflections along the second, third, and fourth floors were much greater than the first, fifth, and sixth floors. As a response, concrete reinforcement was added throughout the structure to better resist the lateral loads, and the concrete wall thickness was increased to 2 m for the second, third, and fourth floors. In addition, the concrete strength used was increased to 50 MPa for the second, third, and fourth floors. The first, fifth, and sixth floors were also increased in wall thickness to 1.5 m.

To further combat the deflections and stresses found within the structure, a steel beam truss system was created to connect the exterior wall to the interior wall, providing additional strutting and load transfer throughout the structure. 360UB44.7 beams were chosen for their relative light weight and strength. Perlucor glass was chosen for the ceiling design for its superior strength and thermal properties compared to standard glass and to fulfill the need to have a transparent ceiling material for light transfer downwards into the structure.


6. Strand7 models

6.1 2D model

A two-dimensional model was created to model the temperature within the structure. A quad 4 plate was created with the thermal properties of air and the dimensions of a vertical cut through the middle of the structure. The structure was then subdivided, and the top layer of plates was changed to the structural properties of glass with the thermal properties of Perlucor. The sidewalls were then converted to the structural and thermal properties of concrete (Figure 5).

Figure 5.

2D model of the earthscraper created using Strand7 cross section showing the perlucor ceiling (red), green concrete exterior walls (green), and the air inside (blue).

6.2 3D model

The model was formed by creating tri-3 elements as the exterior and interior walls. The floors and ceiling were then created using quad 4 elements. The truss beams were then made using beam 2 elements. Each element type was then separated into six groups for their respective floors. The plate elements were then subdivided. Finally, spring dampener beam 2 elements were created by selecting the nodes on each face of the structure and then extruding by increment of 1 m in the negative Z direction or the face’s equivalent axis perpendicular to the vertical axis pointing outward. Fixed conditions were then applied to the base node of the structure (tip of the pyramid), as well as the node of the spring dampeners not attached to the structure (Figure 6).

Figure 6.

Stages of 3D model creation—exterior walls, interior walls and floors, and beam truss (from left to right).


7. Numerical analysis

7.1 Thermal analysis

The transfer of thermal energy occurs in two thermal paths. The Perlucor ceiling at the surface ambient temperature absorbs the solar radiation, transporting the heat internally through conduction and heating the internal air through convection. In addition, the walls of the structure also at the ambient temperature of the surrounding soil conduct this energy through the wall system and expel this energy through convection within the structure.

An ambient temperature of 22.1°C was applied at the top surface, and the linearly increasing ambient temperature of the soil was applied along the sidewalls. A convection coefficient of 11 and 23 W/mK was applied along the Perlucor ceiling and sidewall, respectively, and heat flux from the sun of 69.44 W/m2 was applied along the top surface. The results of the static thermal solver are shown below, alongside the model used.

Figure 7 shows a maximum temperature of 38°C at the bottom vertex, a region designed to harness geothermal energy. This, however, will not be sufficient for energy production according to the thermal analysis. Otherwise, the building reaches thermal equilibrium at comfortable temperatures. Ranging between 19 and 25°C for the upper third of the structure almost eliminates the need for cooling or heating costs. Although the lower section exceeds standard room temperature (25°C), the cooling costs would be lower compared to a standard aboveground structures.

Figure 7.

Temperature contour of steady-state heat solver.

7.2 Earth pressure analysis

7.2.1 Linear static solver

The 3D model shown in Figure 6 was used to analyze the horizontal displacement of the earthscraper under lateral earth pressures. To achieve this, the model was run under the linear static solver, with the following results produced. Figure 8 portrays the horizontal displacements under lateral earth pressures. The internal walls experience a maximum deflection of 28 mm, whereas the external walls experience 527 mm of deflection.

Figure 8.

(a) Interior deflection contour and (b) exterior deflection contour.

7.3 Earthquake analysis

Dynamic analysis can be used to calculate the response under earthquake loading. The acceleration vector,U¨, and displacement vector, U, of a structure under the action of a time-dependent force, F(t), are given by the matrix element equation:


where M is the mass matrix and K is the stiffness matrix of the structure in Eq. (3). This can be solved using the transient solver; however, it is quite computationally demanding. Alternatively, the Strand7’s spectral response solver coupled with the natural frequency solver can be used to determine the structure response under an external loading.

7.3.1 Natural frequency solver

The natural frequencies of the structure can be found by removing the external and damping forces and introducing an oscillatory function, U(t), shown below in Eq. (4). This yields the characteristic equation showed in Eq. (5):


where w is the natural frequency and the solution gives a polynomial function of w2 with order equal to the degrees of freedom of the structure. The number of degrees of freedom equals to the number of modes, and the modal shapes are found by normalizing the vibrational modes, shown below in Eq. (6):


The natural frequency solver uses this method to find the modal shapes and the respective frequency, wi, with the modes and respective frequencies found tabulated below. This is as specified by Clause 7.4, with 25 modes found to have converged. Table 2 shows the different modes and their respective frequencies.

ModeFrequency (Hz)ModeFrequency (Hz)ModeFrequency (Hz)ModeFrequency (Hz)ModeFrequency (Hz)

Table 2.

Frequency and modes of the natural frequency solver.

7.3.2 Spectral response solver

The method of spectral response involves determining the peak response of the structure due to an applied acceleration, which is given by the acceleration response shown in Eq. (7):


The spectral response acceleration was input into Strand7 in compliance with AS/NZS1170.4 Clause 7.2, in which the horizontal design spectrum, Cd(T), was determined and is plotted in Figure 9.

Figure 9.

Design response spectrum.

The maximum displacement is then found by method of square root of the sum of the squares (SRSS) satisfying Clause 7.4.3, in which the contribution of each mode is superposed such that


where Γiis the mass participation factor, a measure of the contribution of mode i. A detailed explanation is given in Appendix C. The displacement contour was plotted using the SRSS method. It is shown in Figure 10 that the maximum displacement due to earthquake loads is 19.8 mm, insignificant relative to the scale of the structure. It should be noted, however, that the design standard AS1170.4 [6] used to generate the numerical analysis does not account for liquefaction, settlement, or fault rupture effects on the soil and as such may be underestimating the effects of earthquake loads on underground structures.

Figure 10.

(a) Exterior deflection under earthquake loads and (b) internal deflection under earthquake loads.


8. Conclusions

This chapter provides a 3D linear static and spectral analysis, as well as a 2D thermal analysis of an inverted pyramid structure modeled off the architectural design of the firm BNKR Arquitectura [8].

Both the 2D and 3D models have been successful in determining the loading conditions of the structure, as well as identifying key structural areas of increased stress, deflection, and heat. Using this analysis, modifications to the structure were made to reduce the effect of these loads and bring the structure within an acceptable level. A final design with a simplified six-floor system arrived upon to allow a sufficiently accurate model while maintaining the viability of 3D modeling within the finite element analysis software Strand7.

While the results found are encouraging, further research should be undertaken to investigate the liquefaction effect of earthquake loads on underground structures and soil. In addition, more detailed structural modeling could be undertaken to investigate the minutia of detail for each individual story, rather than simplifying into six larger general floors.

Therefore, the structure could feasibly be considered in 2050, with a rapidly increasing sophistication of software and material engineering. The demands of the human population in 30 years may necessitate such an ambitious design.


Figure 11 portrays the soil pressures experienced by the external walls of the structure, with the relevant values kept parametrized for the determination of Eq. (10).

Figure 11.

Earth pressure profile.

For substitution in Eq. (10), the parameters and their respective values are shown in Table 3. These values were found through literature review and govern the lateral earth pressures experienced by the structure.

Internal friction angle, ϕ°47.2
At-rest coefficient of earth pressure, K00.26627
Soft soil unit weight, γsskN/m317
Water unit weight, γwkN/m310
Hawkesbury sandstone unit weight, γhskN/m324
Depth of soft soil, dssm10
Depth of water table, dwm60

Table 3.

Earth pressure parameters.

The lateral pressure is given as


where u(z) is the pore water pressure, given as the hydrostatic water pressure shown below


and σhzis the effective lateral pressure, expressed as


where K0is the coefficient of earth pressure at rest, given by the following expression given by Jaky [7]


and σvzis the effective vertical earth pressure. This is the contribution of the different soil layers less the effect of pore water pressure, given as


The Macaulay brackets seen above operate as follows:


This allows one expression for earth pressure at any depth. Replacing all the relations together gives


At the bottom of the structure, i.e., z = 300 m, the lateral earth pressure is calculated as 3.66 MPa, shown below in Eq. (10) using the parameters in Table 3:


The thermal properties of the materials used in the structure are given in Table 4. These parameters determine the conductivity and transfer of heat throughout the structure. The reinforced concrete conducts a relatively higher amount of heat compared to the Perlucor ceiling, both absorbing external heat from the soil and surface, respectively.

MaterialReinforced concretePerlucorAir
Coefficient of thermal energy W/mK1.370.780.0257
Specific heat capacity J/kgK8808401005

Table 4.

Material thermal properties.

The heat flux on the Perlucor ceiling can be estimated for a blackbody using Stefan-Boltzmann law. The heat flux generated at the sun’s surface, qsun,is given by


Hence, for radius, R, of the sun and distance D from the sun to the earth, the heat flux experienced by earth is


The parameters needed for Eq. (17) are shown in Table 5, giving the value of

StefanBoltzmann constant,σ(W/m2K4)5.67×108
Temperature ofSun,TsunK5780
Temperature of space,TspaceK3
Distance fromsun,Dkm1.5×108
Radius ofsun,Rkm6.957×105
Radiation coefficient ofsun,ε1

Table 5.

Heat flux parameters.


This however assumes the earth acts as a blackbody, absorbing all incoming radiation. Measurements by the Bureau of Meteorology have shown that the incoming solar radiation is of the value of 69.44 W/m2, significantly lesser than the blackbody value. This is due to solar irradiance, cloud reflection, atmospheric absorbance, and the reflection of the earth surface. The thermal boundary conditions are shown in Figure 12.

Figure 12.

Thermal conditions.

Pells (2004) determined the bearing capacity of Hawkesbury sandstone as 20–103 MPa and 6–14 MPa for laboratory tests and field tests, respectively. Taking the field tests as true, this classifies the rock as subsoil Class Be (AS/ZS1170.4 Clause 4.2.2).

Using AS/NZS1170.4, the hazard factor, Z, in Sydney is 0.08. For an important factor of 3 (BCA 2016), the structure has an annual probability of exceedance of 1/1000 (BCA 2016).

AS/NZS1170.4 gives the probability factor of 1.3 for an annual probability of exceedance of 1/1000.

AS/NZS1170.4, for a structure height of 300 m, the earthquake design category is EDCIII. Structural ductility factor (μ) and structural performance factor (Sp) are found as 2 and 0.77, respectively (AS/NZS1170.4).

The spectral shape factor, ChT, is found by using AS/NZS1170.4 Clause 6.4. This requires the period of vibration of the structure, which was found using the natural frequency solver. The horizontal design response spectrum,CdT, is found using AS/NZS1170.4 Clause 7.2, given as


The scaling factor is multiplied by the gravitational acceleration to give the horizontal acceleration in terms of m/s2. This process is summarized in Table 6, showing the modes, the respective frequency and period, then the calculated and design spectral response factor, and then the mass participation factor. The total mass participation of the 25 nodes is 99.461%, satisfying AS/NZS1170.4 Clause 7.4.2.

ModeFrequency (Hz)Period (s)Ch(T)Cd(T) (g)Cd(T) (m/s2)Γi%

Table 6.

Natural frequency results and spectral response.


  1. 1. Sol Robles J. La Redención del Ornamento: Recuperando la Dimensión Simbólica de la Arquitectura. RevistArquis. 2012;1(2):26-35
  2. 2. McDonald P, Kippen R. Projecting future migration levels: Should rates or numbers by used? People and Place. 2002;10:82-83
  3. 3. Morse JM et al. Developing Grounded Theory: The Second Generation. New York: Routledge; 2016
  4. 4. Trahair N, Bradford MA. Behaviour and Design of Steel Structures to AS4100: Australian. New York: CRC Press; 2014
  5. 5. Gilbert R. Creep and shrinkage models for high strength concrete–Proposals for inclusion in AS3600. Australian Journal of Structural Engineering. 2002;4(2):95-106
  6. 6. AS AS. Glass in Buildings-Selection and Installation. Standards Association of Australia. 1989
  7. 7. Jaky J. Pressure in silos. In: Proceedings of 2nd ICSM. Japan: ICSM; 1948
  8. 8. Hudson D et al. ACCESS-S1: The new Bureau of Meteorology multi-week to seasonal prediction system. Journal of Southern Hemisphere Earth Systems Science. 2017;67(3):132-159
  9. 9. Sharma R, Richards P. The influence of Helmholtz resonance on internal pressures in a low-rise building. Journal of Wind Engineering and Industrial Aerodynamics. 2003;91(6):807-828

Written By

Faham Tahmasebinia, Kevin Yu, Jiachen Bao, George Chammoun, Edwin Chang, Samad Sepasgozar and Fernando Alonso Marroquin

Submitted: April 5th, 2019 Reviewed: June 3rd, 2019 Published: May 13th, 2020