Influence of an extensive green roof in the non-linear structural behavior of self – construction households Influencia de un techo verde extensivo en el comportamiento estructural no lineal de viviendas de autoconstrucción

In recent years, ideas such as green roofs gained importance as an alternative for environmental adaptation of society to new ambient conditions. In Bogotá city, the roofs at popular neighborhood represent a potential area for their implementation. This present research discusses how the weight of a green roof impacts a self-construction dwelling unit from a structural capacity demand point of view. Green roof weight is a random variable whose normality was investigated through the Shapiro-Wilk goodness of fit. The magnification factors for this dead-like load come from the Nowak-Collins method using the confidence intervals on the mean weight value. The ATC-40 pushover methodology helped in finding the non-linear behavior of a controlled but similar structure as those found in the popular neighborhoods of Bogotá city. According to this methodology, damage caused on these structures due to seismic demand increases by the green roof incorporation. However, damage increase is a function of local soil conditions.


Introduction
Over the last few decades, the effects of climate change and environmental degradation left consequences throughout the world. Environmentally sustainable ideas represent a way to mitigate the effects. From an environmental point of view, green roofs are defined as a tool that reduces the thermal sensation in urban areas while retaining runoff by alleviating the flow of rainwater in the sewer system, by improving air quality (Aboelata, 2021).
From a structural perspective on one hand, green roofs can be defined as just one added load on the roof of a house. On the other hand, this load tends to be variable depending on the climate conditions of the area, due to changes in weight of the substrate of plants as saturation can happen. With the increase in the mass of the last story of a house, the potential inertial force of this story also increases as a result of an external force for example seismic forces. As consequence, the horizontal displacement of the structure also increases. Thus, increasing the upper story mass in a structure, increases the probability of developing a new failure mechanism. For the present research a Self-Construction-Households (SCH) define an unregulated popular masonry dwelling unit, that in most cases lacks structural design. However, these SCH constitute also a large area in the city for green roof potential implementation. In the municipality of Soacha there was a productive green roof alternative using plastic bottles containers in a house built with modular concrete panels embedded in thin sheet profiles made of galvanized steel.
The roof of this house got four uniaxial seismic accelerometers and through computational modeling, the researchers reported an increase of 241% in the tensile stresses associated with the load combinations and a 56% increase in natural structural period (Olaya et al., 2014). A study using numerical methods characterized a reinforced concrete moment resistance frame, with a 3, 6 and 8 story configurations. This structure was designed in accordance with current regulatory standards without accounting the weight of green roof. Structural natural periods increased 5.75%, 3.03% and 2.29% for the 3, 6 and 8 story frames respectively, compared to a normal gravel roof finish (Bianchini et al., 2016). Within a life cycle investigation for green buildings carried out in 2017, non -linear plus fragility analysis helped in analyzing a 4-story-high commercial building designed in reinforced concrete. Intensive and extensive green roof models experienced a slightly earlier loss of lateral load capacity at a roof drift ratio of around 4% due to significant P-Delta effects. This is evidence of an increase in the structural period of 5.34% and 8.15% respectively with respect to the structure without green roof _bookmark48 (Welsh-Huggins & Liel, 2017).
Current knowledge of the influence of a green roof from a structural point of view focuses on properly design of reinforced concrete and precast concrete structures. The present research studies the non-linear behavior of typical SCH in Bogotá based on moment-resistant concrete frames neither with structural design nor confined masonry under the presence of a green roof.

Green roofs
Green roofs are a sustainable means that help to mitigate negative effects on the environment. Amongst the main advantages they represent there are tax benefits for the owner of said green roof, better air quality, rainwater retention and an increase in the commercial value of the property (Manso et al., 2021). (Figure 1) shows an alternative construction option. This demonstrates the community's interest in this sustainable mechanism and allows the development of research that contributes to the analysis of both environmental benefits and structural effects.
Green roofs are globally classified as extensive (thickness<15cm) or intensive (thickness>15cm) (Nektarios et al., 2021). These are generally made of a waterproofing layer to avoid contact between the support structure and biological elements, a drainage layer which allows the passage of water but not solid elements, a filtration layer where the substrate is contained, a layer of substrate and finally the vegetation layer. The growth medium is made up of materials that increase in weight when a humid condition is present. In such a way green roof mass depends directly on the climate characteristics of the area where it is, the type of plants in the vegetal layer, the retention capacity of the roof in its different layers and the evaporation rate (A'saf et al., 2020).   (Apuntes, 2020) 0.02 and 0.01, let to a Normal probability distribution function. Using different thicknesses for substrate is a common practice allowing different vegetal species during a green roof solution design. Having the information of 11 projects, helped in characterizing these green roof solutions accounting local solutions based on the catalog of vegetation in domestic gardens in Bogotá (Sierra Guerrero and Amarillo Suárez, 2014) and the Bogotá green roof guide. Subsequently, the probable thickness that each type of plant should have to survive local weather conditions is a function of natural conditions. A satellite image analysis helped in determining the approximate area covered in green roofs. A safety factor calibration methodology helped to find the average weight approximately at 208.05 / ! . The following steps show the methodology followed to determine said factors.

Green roof weight
• Step-1. Formulate the limit state function and the design equation. • Step-2. Obtain the initial design point for the number of random variables ( − 1) name it ! * . • Step-3. Compute the normal equivalent for the non-normal probability distribution functions of the random variables ! ! ! , ( ! ! ! ). • Step-4. Determine the partial derivatives of the limit state function. • Step-5. Calculate the column vector using (Equation 1): • Step-6. Determine the new design point in reduced variants (Equation 2). * = Where !"#$%! is the reliability of the target to be found. Where is the resistance, is the green roof dead load and is the live load. The design equation is (Equation 7): Where , ! and ! are safety factors for resistance, dead load and live load respectively. In step 2, the design points correspond to the mean value of the random variables (identified with an asterisk as a superscript). For live load, a Gumbel type I probability distribution function is used (Arroyo Amell and Ramos Castillo, 2015), and for structural resistance a LogNormal distribution (Nowak and Collins, 1976). A ratio between live load and twice the dead load is used as suggested by Nowak and Collins, and from the following equations (Equation 8): The design point of the resistance is, then (Equation 9): * = (9) In step 3, the equivalent for resistance as a normal probability function has the following parameters (Equation 10) and (Equation 11): where ! is the coefficient of variation for resistance, which for reinforced concrete elements is considered as 0.12 (Arroyo Amell and Ramos Castillo, 2015). The normal equivalent for the probability function of the live load uses the following parameters (Equation 12) (Equation 13) and (Equation 14): where ! is the coefficient of variation of the live load, which is considered as 0.25 (Galambos et al., 1982) (Equation 15) (Equation 16): The partial derivatives obtained from step 4 are (Equation 17): Where ! is the coefficient of variation of dead load, which according to results previously obtained corresponds to 0.23. For step 5 the matrix indicates the relationship between the random variables under consideration. Each variable is considered to be independent (Nowak and Collins, 1976 Convergence rapidly reports a safety factor for green roof weight of 1.1, according to (Equation 26): So, then for weight, a normal distribution function works with a mean value of 228.86 / ! and a standard deviation of 52.65 / ! as shown in (Figure 2).
A confidence interval on the mean value for a green roof seems a better tool for all probability purposes. Thus, the following equations (Equation 27)

Self-Construction Households (SCH)
Urban planning failures from the conception of Bogotá city, led to an unplanned, disorganized and unqualified construction activity know as SCH, which currently populate vast areas of the city. Due to budget constraints, a formal construction is not a solution for a large portion of Bogotá's population. Thus, SCH are a solution to shelter families within the city. This phenomenon started to grow as early as the 1960s (López Borbón, 2018). Even though SCH have basic services such as water and electricity, they lack structural design, representing a liability for owners. Ancient construction knowledge is part of Colombia architectural heritage. However, it doesn't account the need for a rigorous study and seismic capacity, which nowadays, is mandatory. In the following part of this document, there is a brief description of a typical SCH in terms of geometry, materials and their constitutive models used to perform non-linear analysis.

SCH scope geometries
A total of three SCH models of 1, 2 and 3 stories helped to understand the structural response when a green roof is present, after construction. The distribution of spaces is such that 35% of the length of the walls are in a transverse direction and 65% in a longitudinal direction. (Figure 3) depicts the floor plants for each one of the stories for a typical SCH.

Materials
Structural analysis used OpenSees including three main materials for structural definition: concrete, masonry and steel. The degree of concrete confinement is a function of transverse reinforcement separation. According to the model proposed by (Mander et al., 1989), with a separation of 0.2m between stirrups, the effective confinement stress is zero. The constitutive model for concrete follows the parameters described in (Mohd Yassin, 1994) (Table 1) (Figure 4).  (López, 2016)  For steel hysteresis, the constitutive model presented in (Scott and Filippou, 2016) was part of the present analysis, with the following parameters: (Table 2 Table 2. Parameters for steel constitutive material Figure 5. Constitutive model for steel (Scott & Filippou, 2016) Where 1 is the yield stress in tension, 1 is the yielding strain, 2 is the ultimate tensile strength (UTS), 2 is the strain at the UTS, 3 is considered as 20% of the last stress in tension and 3 is the corresponding strain. (Marques et al., 2020). In order to model the impact of masonry on the structure, the equivalent diagonal strut model is used. Figure 6 presents this model, which consists of a diagonal element whose width depends on the relative flexural stiffness of the masonry with respect to the columns of the surrounding frame (Al-Chaar and Lamb, 2002). Said element is modeled to receive only axial forces and since it is considered a constitutive model of zero tensile strength, in the pushover analysis only one of the two elements contribute to the stiffness of the structure.

The same magnitudes are considered in the compressive direction of stresses and strains due to the isotropic properties of steel. Confined masonry is widely used in Latin America, however neither structural analysis nor safety controls are mandatory for its design, mainly due to limited local resources and the lack of an effective management and control policy in the construction sector. The diffusion of this type of structures towards a more efficient application requires a better understanding of their structural behavior
A set of two typical masonry units available locally helped in modeling the compressive response of such composite element. The ANSYS model is available in (Figure 7).

In this model, the assumed compressive strength for mortar !"
! was of 7.5 MPa (Valbuena Porras & Mena Serna, 2011) having a thickness of 17mm. This value varies from 13 mm (Rincon, 2010) to 20 mm (Fundación Suiza de Cooperación para el Desarrollo Técnico, 2016). Since this thickness is variable and depends on the builder's skills, a mean value of 17 mm is chosen with a lime-to-sand ratio of 1:4 commonly used for simple masonry (Sánchez de Guzmán, 2001). The Kent-Scott-Park model (Filippou and Mazzoni, 2014) for masonry helped in setting the parameters as masonry constitutive model. These parameters are available in (Table 3) ( Figure  8).
Where is the masonry compressive strength, 0 is the maximum associated strain, is the crushing strength of the material and is the strain at the crushing strength.

Load pattern applied for the non-linear analysis
The application of the Pushover methodology involves the gradual pushing of the structure in one direction to observe its behavior as it goes into the non-linear range. However, the applied monotonic load pattern can follow several options. There are reports of lateral loads following an increasing triangular or uniform shape along the structure's height (Abhilash et al., 2009). However, inertial loading patterns arise with a debate about its applicability to seismic demands (Rahmani et al., 2019). In order to consider this effect, a lateral load pattern could be proportional to the main mode of vibration of the structure (R & S, 2000). This means neglecting the contribution of the other modes of vibration. For this research, the triangular load pattern was appropriate, since it provides an excellent approximation of the seismic capacity curve for medium and low rise structures (Sun et al., 2003). At the same time even though the load distribution is not accurate, the errors are on the security side as it underestimates the structural capacity (Khoshnoudian et al., 2011). The triangular shape, and the load distribution is available in (Figure 9): Table 3. Parameters of masonry constitutive material Figure 8. Kent-Scott-Park constitutive model for masonry units (Filippou and Mazzoni, 2014)

Demand and capacity curves
To calculate the performance of a buildings with a specific target displacement or until the structure reaches a collapse limit which allows to show the weaknesses of the original design (Bhandari, 2020). In other words, this analysis finds a correlation between applied lateral load and roof sway. The structure is subjected to gravity loads and static lateral loads that increase monotonically. (Figure 10) presents the results obtained from said analysis applied to the scope SCH. In all cases there is evidence of a decrease in the capacity of the structure with the presence of an extensive green roof. It is interesting to see that the smaller the mass of the structure, the higher the force-to-base shear ratio.
The seismic demand in Bogotá uses on the seismic microzoning of the city. (Figure 11) shows the distribution of the 16 seismic response classifications available as a function of soil type.
Where ! is the shear at the base, is the weight of the structure, ! is the mass participation factor, is the magnitude for earth's gravity acceleration, (9.81 / ! ), !""# is the displacement at the control node, ! is the modal participation factor. And !,!""# is the modal amplitude with respect to the − ℎ vibration mode. The subscript indicates the vibration mode corresponding to the analysis. (Equation 31) For the demand curve: Where is the period of the structure in seconds. Figure 11. Seismic microzoning of Bogotá (IDIGER, 2020)

Performance point
The performance point is known as the intersection between the capacity curve and the demand curve in ADRS coordinates. (Figure 12) shows the intersection of the capacity curve with the demand curve, taking place in early values of ! . According to the displacement-based design approach, it is necessary to define seismic performance levels, which are associated with a level of damage exhibited by structures that are commonly identified through deformations or drift limits (Aşıkoğlu et al., 2020). For this study, the performance levels base their criteria on experimental investigation performed in Peru where a house made of perforated bricks, studied the effects of horizontal forces acting on the structure. As a result, the following ranges of levels of performance are available: at 13% of the curve measured horizontally corresponds to an immediate occupancy ( ) level where the damage is minor, at 32% the performance level is life safety ( ) with moderate damage, and finally with a 39% collapse prevention ( ) is the corresponding level (Zavala Toledo, 2019). The performance limits are plotted on the capacity curve in order to show the level at which the performance point happens. For example, for a two-story house without a green roof, placed at the "Hills" soil type, the performance is at the life safety level ( ). This procedure repeated the procedure to study the 1, 2 and 3-story house models in each of the 16 types of soil condition. (Figure 13) shows the performance status for the analyzed models considering the action of the green roof load.
It is evident that the use of a green roof on these structures can cause in certain areas of the Bogotá a noticeable increase in the performance level expected from the structure in terms of seismic activity. The hills soil type seems to be the only one that does demand lower forces at the SCH, remaining the performance level in ( ).  (Welsh-Huggins and Liel, 2017) presents an increase in the capacity curve for a reference reinforced concrete structure. This, due to the use of a shallow and a deep green roof. The results obtained in the present research are similar even though the present study is made of concrete and masonry. In contrast, in the work done by (Olaya et al., 2014) and (Bianchini et al., 2016) there is an increase in floor drift. It is not feasible to make a numerical comparison since the structures are different; while in the first reference a one-story prefabricated structure was the scope of study, in the second one a reinforced concrete story structure was the scope of research.
In both, an increase in the structural period is evident. For the present work there are also increases in the period of the three scope structures (1-story, 2-story and 3-story SCH). (Table 4) presents the natural periods obtained, in each case.

Conclusions and future work
The conclusions of this work focus on two parts; the first part regarding the weight characterization of the green roof, and the second one regarding the structural impact of the implementation of this load on the structure of a SCH.
The mean weight of an extensive green roof in saturated condition with a 95% probability can be found between 1.925 / ! and 2.565 / ! . It is also clear that as the SCH gets higher, the structural vulnerability gets higher, with an increase in this condition by the green roof presence. For one-story structures, it is possible to implement a green roof in the northwestern part of the city. For structures with more than one story, it is not advisable because the structure, when subjected to an earthquake, can go from suffering moderate damage to severe damage and even reaching a stage close to collapse.   Table 4. Natural periods of SCH with and without green roof