Simulation on Sequential Construction Process and Structure of the Pisa Tower

The leaning of structures happens all around the world and generates impacts on different extents; thus, it is important to learn about the causes behind. In this report, the sequential construction of a typical leaning structure, the Tower of Pisa, is discussed and simulated by using a finite element code, PLAXIS. The simulation is performed on a two-dimensional plane with simplifications taken into consideration in making modeling feasible under limi-tations. Three distinct models are built with one as a control variable, while the other two models are set up with exact eccentricity. Data are obtained from the analysis and are plotted in a graph to clearly show the relationship between the tilting angle and construction phases. With reasonable and completed simulation, the study is able to show the significant role compressible subsoil plays in impacting the tilting performance of a tall building.


Introduction
The leaning tower of Pisa is one of the most remarkable architectural structures in medieval Europe. Constructed in 1173, the tower is round and is constructed throughout of white marble, inlaid on the exterior with colored marbles. The uneven settling of the campanile's foundations during its construction gave the structure a marked inclination that is now about 17 feet (5. Currently, the tower is approximately 60 meters tall; the highest side is 56.67 m, and the lowest side is 55.86 m. And the outer diameter is 15.484 m. The overall form of the tower is a hollow cylinder subdivided into eight orders, with spiral staircase leading to the top. Colonnade is engraved on the wall of the base floor; the middle six stories are surrounded by marble columns and loggias; a belfry, which is always denoted as the Bell Chamber, is constructed at the top. The significant events in the construction history of the leaning tower of Pisa were shown in Table 1.
Began in 1173, the tower's construction process was interrupted after 5 years when three and one third floors were completed, and the tower was leaning to the north about 0.25 degrees. The tower continued to move towards a more upright position when the work finished [2].
However, when the 7 th floor was finished in 1278, the tower tilted back to south by about 0.6 degrees and increased to 1.6 degrees in about 90 years. In 1817, the first measurement of the tilting angle was recorded, showing that the angle is about 5 degrees toward the south. After centuries of repairs and adjustments, the inclination decreased approximately to about 4 degrees. As Pisans say, the tower is "banana-shaped", with the bell tower 1.5 degrees closer to the vertical than the base [3]. The inclining angle of the tower over time could be seen from Figure 1.
The objective of this research is to simulate the sequential construction of the leaning tower of Pisa. One of the important reasons is that the tower of Pisa is currently the world's most ancient and famous leaning structure which means that there are more historical recordings and articles published about this leaning tower. Those can give us more freedom in choosing data for simulation and the tilting of the Tower of Pisa itself is also a very representative architectural problem.
To this end, the finite element software PLAXIS 8.2 [5] is utilized to capture the process of staged construction. This is because PLAXIS has the functions to simulate the consolidation problems and step construction. Those functions are required for this simulation to be achieved. Even though the version of the code restraint the study to two-dimensional plane strain modeling, and the premises  do not accurately reflect the geometrical features of the Pisa tower, the proposed study is useful in gaining insight into how subsoil foundations impact the performance of a tall structure. The study also provides an opportunity to utilize an advanced finite element code in modeling complex boundary-value problems, such as those arising from the sequential construction of a world heritage.
Soil Structure Interaction (SSI) is the process of mutual interaction between soil and structure, that is, under external forces, the response of the soil influences the motion of the structure and vice versa. Conventional structural design methods usually neglect the SSI effects, which is feasible under light structures in relatively stiff soil, but the effects become prominent when heavy building is resting on soft soils; the leaning tower of Pisa, which is mostly comprised of heavy stones and was built on compressible soil, can be a typical example of the latter situation. Therefore, this research is a simple exploration into the soil structure interaction between the leaning tower of Pisa and the soil beneath.

Subsoil Condition
The tower of Pisa was built on highly compressible soils and started leaning at the very beginning. The ground profile beneath the tower consists of three separate layers as shown in Figure 2. Horizon A is about 10 m thick, consisting of a  ground water elevation is 3 meters above the mean sea level, while aquifers in Horizon A, located at a depth between 1 m and 2 m, create challenging problems on protecting the tower. This information is quite essential since the depth of soil layer strongly affects the structural response [6].

Assumptions on the Tilting Structure
Based on the aforementioned background of the tilting structure, following assumptions can be made. The first one is that the mass and weight of the building are uniformly distributed. Then the building is rectangular in shape and the building's deformation is small compared to its rigid-body motion. There is no relative movement between the soil and the building at their interface. The material of the building is elastic and the building is tilting in the plane defined along width and height; therefore, it will be a two-dimension simulation. Finally, the ground elevation is 0 m.

Numerical Model
Some simplifications of the problem are conducted to make the simulation feasible under limitations. The structure was concerned to be an 8-storied perfectly vertical building in height of 60 m, in length of 100 m, in width of 19 m, with 5 m embedded underground as shown in Figure 3.
When simulating, the building was assumed to tilt in a certain direction. Additional information is shown in Table 2.
The subsoil profile was assumed to consist of four layers in Mohr-Coulomb model to simulate their constitutive behavior. The profile is given in Table 3 and   Table 4.

Finite Element Model
In order to simplify the loading condition, the weight of each floor is taken as concentrated load at the centroid as shown in Figure 4, which has the same effect as distributed load.
The initial simulation model (model 1) is a vertical ghost building with vertical point loads on the centroid of each floor as shown in Figure 5. The structure of sand layers will cause uneven sinking and make the building tilt without other setting. However, the 2D model in PLAXIS is unable to simulate non-linear conditions so that an angle of tilt in the initial setting is required to generate an eccentricity which can make simulation closer to realistic condition.  In order to locate the load with exact eccentricity, it is assumed that the entire building starts to tilt at 0.6 degrees as the construction begins. The geometry model (model 2) could be built in 2D as shown in Figure 6.
The building could be defined as crossing plates with high stiffness with deformation ignored. The weight of each plate is included in the overall weight of the floor; therefore, the input magnitude each member weighs should be zero. Table 5 shows the properties of each member.
The finite elements are produced by meshing the entire geometry with medium element distribution and refining the cluster around the box foundation for more precise calculation. The meshing result is shown in Figure 7.

Simulation Results
The equivalent point load could be defined by: Therefore, the simulation results could be obtained from PLAXIS as shown in Table 6 and Figure 10.

Advantages and Weaknesses
There are several advantages about the current modeling process. The numerical model can be set up with ease, since the shape of the horizontal section of the Journal of Building Construction and Planning Research    Figure 10. Simulation results of all models.
Pisa tower is a symmetrical circle and the weight of each floor is uniformly distributed, it is reasonable to regard the weight as concentrated load at the centroid. Similar leaning angle could be easily simulated in two-dimension model.
Also the model has its limitations. The first one is that PLAXIS 2D is more suitable in simulating buildings with the same sectional surfaces. However, the tower of Pisa is a hollow cylinder with certain degrees of leaning, meaning that the sections are different; therefore, the results derived cannot be applied to real situations. Soil structure of the 2D model is over-simplified, for the Pisa tower, in reality, first tilted to the north and then to the south during the construction process. In addition, the presence of a gap between the soil and foundation of tower in cohesive soil is ignored, which has a significant effect on the capacity performance [8]. The simulation also lacks precise information, including soil condition and water levels, making it much more constrained.

Conclusions
In For model 1 and 2, the uprising tendency of the tilting angle corresponds to the initial assumption that the angle of tilt has a positive relationship with the weight.
In the simulation above, the external force comes from the leaning tower of Pisa itself. Since the tower began tilting when the fourth order was constructing, the gravity from the above upper orders can provide the simulation with external force. Moreover, in the models above, certain areas around the building have been circled and meshed with finer elements, showing that the simulation focuses on the direct interaction. The SSI effect is reflected by the simulated tilting angle between the Pisa Tower and the plumb line, displaying the significant role that compressible subsoil has played in the performance of a tall tilting building.