Research Article | Open Access
H. Rodrigues, R. Sousa, H. Vitorino, N. Batalha, H. Varum, P. Fernandes, "Characterisation of Portuguese RC Precast Industrial Building Stock", Advances in Civil Engineering, vol. 2020, Article ID 7517205, 19 pages, 2020. https://doi.org/10.1155/2020/7517205
Characterisation of Portuguese RC Precast Industrial Building Stock
The construction of a vulnerability model requires reliable information on the features of the buildings in the study. The purpose of this work is the characterisation of the precast industrial buildings in Portuguese industrial park, based on the survey of 73 design projects of existing buildings. The collected data are based on a previous study on the features that influence the seismic response of this type of buildings. The parameters collected are associated with the global geometry and specific elements characteristics (e.g., column dimensions, reinforcement ratios, and connections details), to the mechanical properties of the materials and other parameters that can give some important information in the characterisation of the buildings (e.g., construction year and localization). In the end, a comparison with other available databases, namely, from Italy and Turkey, is done in order to conclude about the similarity. This information is important to define representative experimental specimens and numerical simulation to conduce seismic risk analysis.
The development of seismic risk studies requires the knowledge of the building characteristics of a given typology under analysis. The information of geometric and mechanical parameters allows the definition of reliable numerical models that can be used to derive fragility functions capable of describing the relation between the seismic intensity and building limit states.
Contrarily to what is observed for residential buildings, limited information is available in regard to the properties of precast buildings. To the authors’ knowledge, in Portugal, only the study developed for the European Commission  presents a general description of the main typologies of the industrial building stock, but its information is insufficient to enable the development of representative numerical models.
To overcome the previous limitation, a survey was carried out analysing the structural design project of dozens of industrial reinforced concrete (RC) buildings built in the Portuguese continental territory over the last 50 years. The statistical analysis of the information gathered provides indication on the material properties as well as the geometry typically used in industrial precast RC buildings.
The present work provides the statistics associated with the properties that influence the seismic performance of these buildings. The information collected enabled the characterisation of the global geometry and the mechanical properties of the materials, as well as other local systems that may influence the seismic response of these buildings, such as columns size and reinforcement ratios, beam-to-column connection, and cladding systems.
The results obtained were also compared against previous studies conducted in Italy  and Turkey , allowing the identification of common structural characteristics. The analysis of previous studies is also important to assess the possibility to derive common properties for the Mediterranean countries and assess the possibility to use common database that can be used to define representative experimental specimens and numerical simulation to conduce seismic risk analysis.
2. Portuguese Industrial Building Stock
This report attempts to identify and quantify the main structural properties of the RC precast buildings in Portugal. Considering that these buildings represent only a fraction of the total industrial facilities in Portugal, it appears important to put this typology in perspective within the Portuguese industrial reality.
The study carried out for the European Commission  presents a general description of the main characteristics of the industrial building stock in Europe. Despite being mainly focused on prefabricated steel structures, this report provides also important indications on the overall characteristics of the industrial building stock, which includes RC precast buildings. As regards the buildings typologies, Figure 1 indicates that RC precast structures represent about 1/3 of the total industrial park in Portugal.
Despite not being directly linked to the buildings characteristics, the type of activity developed (Figure 2) provides an indication on the size of the buildings and type of load admitted in the design process, which may condition the size of the structural elements. For example, in heavy industries, the presence of large capacity cranes is expected, as opposed to what is expected, for example, in warehousing buildings.
Following the previous comment, it is not surprising to observe that the majority of the industrial buildings have only one bay along the main direction (assumed as longitudinal direction) and that in only 9% of the cases the buildings feature 3 or more spans (Figure 3).
Regarding the main geometrical properties, Figure 4 shows that the majority of buildings have a main span length between 20 and 30 m, whilst, along the transverse direction, the distance between two consecutive frames is essentially between 5 and 7 m. Finally, in terms of columns height, the study reveals that typically this varies between 6 and 8 m (Figure 5).
3. Characteristics of RC Precast Buildings in Portugal
The database presented in this chapter was built based on the information collected after analysing 73 design projects of existing RC precast buildings in the Portuguese mainland. The identification of the buildings sought to reflect an adequate geographical and temporal representation. Regarding the geographical distribution, Figure 6 shows a comparison between the location of the collected projects and the actual manufacturing industry according to the data available in Pordata . The figure shows that, in both images, the buildings tend to be concentrated at the centre and northern regions of the Portuguese coast side.
Regarding the construction period, the buildings analysed were built over the last 50 years showing a clear concentration after 1990 (Figure 7(a)). A first reason for this concentration is related with the very limited, or absent, information in the design project of older buildings. Moreover, until the 60s, the precast industry in Portugal was mainly focused on the construction of elements for slabs and cladding panels. It was only during the 70s that an important growth in the precast systems could be observed . Figure 7(b) shows also the correlation between the year of construction and the seismic zones, showing that, for the buildings consulted, there is a tendency for a concentration of newer structures in the higher seismic zones (lower seismic hazard).
In addition to the location and period of construction, the characterisation of the RC precast buildings involved the collection of several structural parameters regarding the global geometry, material properties, columns section details, and beam-to-column connections, as described in Table 1. The data, and corresponding statistics, associated with each property are presented in the following sections, together with additional, indirect, structural parameters that could provide important indications on the design principles adopted for different periods or regions.
In the following sections, the evolution of the properties defined in Table 1 is evaluated considering a total of 73 buildings. However, it should be mentioned that in some design projects it was not possible to identify all properties, namely, the type of claddings and the beam-to-column connections.
3.2. Global Geometry
The first geometric property addressed regards the identification of the main typologies of precast RC buildings. Based on the survey carried out, it was possible to identify 4 main typologies that are illustrated and described in Figure 8 and Table 2, respectively.
The most traditional configuration consists of prestressed beams, with up to 45 m spans, simply supported on precast columns fixed at the base . Yet, this configuration may present deviations regarding the properties of the longitudinal beams, which can have variable (T1), constant (T2), or truss (T3) configuration, or at the column properties which can feature a “Y” shape at the region of the connection with the longitudinal beams (T4).
Based on the collected data, it is apparent that the majority of the buildings in Portugal correspond to the Typology 1 (86%) and Typology 3 (11%), while only 2% of the cases are associated with Typology 2 and Typology 4 (Figure 9).
As regards the number of storeys, it is apparent that the majority of the precast buildings are single-storey buildings (Figure 10(a)) with a total height that is below 10 m for the majority of the cases (Figure 10(b)).
The values of building height are in line with the study conducted in  for the general industrial facilities in Portugal, whose results are illustrated in Figure 5. Regarding the number of storeys, it is important to note that, except in one case, the upper storeys are only representative of a fraction of the total area of the building (Figure 11).
In terms of plan geometry, the properties vary significantly depending on the building direction analysed. In the direction along the longer beams (assumed as longitudinal hereinafter), the number of spans is generally low (1 or 2), and the beams length can reach values up to 50 m (Figure 12). In addition, when analysed with respect to the construction period, the spans length seems to increase with the year of construction, which is probability related with improvements on the manufacturing and construction processes (Figure 13).
On the other hand, in the transverse direction, the number of spans is typically higher and features smaller lengths, up to 15 m (Figure 14).
Finally, based on the data collected, three main types of cladding systems were identified: infill masonry and horizontal and vertical RC precast panels. Among these, the use of horizontal panels seems to be most common, being observed in 50% of the buildings, while infill masonry and vertical panels were identified in 33% and 17% of the cases, respectively. The increase in the use of vertical cladding panels in the last decades (Figure 15) indicates that this solution appears to become more appealing, as a consequence of improved precast construction processes as well as transport and installation methods. On the other hand, it is still somehow surprising that infill masonry represents an important percentage of the buildings in the last decade and about one-third of the total number of buildings.
3.3. Material Properties
The material properties presented in the following plots refer to the classes indicated in the design projects and, therefore, do not necessarily correspond to the actual values found in the existing buildings.
Regarding the concrete class, Figure 16 shows a large dispersion of the concrete compressive strength (corresponding to the cylinder test) considered in the design process, despite the apparent larger employment of concrete strengths between 20 MPa and 30 MPa. Notwithstanding the great variability, when examined with respect to the year of construction, it seems clear that the concrete classes increase with the period of construction.
Regarding the reinforcement, the number of classes is much lower than the ones found for concrete, and it is apparent that most of the RC members were built with A400 and A500 steel grade (Figure 17(a)). As for the reinforcing steel, it is possible to observe an evolution of the steel grade with the year of construction (Figure 17(b)).
3.4. Column Properties
In this section, the properties of the columns, namely, the section dimensions and reinforcement ratios, are described. The emphasis given to the properties of the vertical elements stems from their relative importance in the seismic behaviour of these structures, with respect to the beam members that, by virtue of their properties and structural arrangement, are expected to remain essentially undamaged. In terms of columns cross-section, Figure 18 shows the frequency and associated normal and log-normal distributions of the columns width and length (defined as the dimension along the longitudinal direction).
The analysis of the results revealed that in only about one-third of the cases the columns are square. At the same time, the collected data indicated that, for the majority of the rectangular columns, the length-to-width ratio remains below 2 (Figure 19(a)).
Another relevant geometric parameter regards the height-to-length ratio, as it provided an indication on the target aspect ratios generally considered in these buildings. Figure 19(b) indicates that, for the analysed cases, this ratio varies essentially between 12 and 25, corresponding to large slenderness and, is therefore indicative that these columns are more prone to develop flexural type of failures, at least in cases where the cladding system is continuous along the height of the column.
In order to assess a possible reason for the variability observed in the relation between the height and length of the columns, this ratio, was plotted against the year of the project and the seismic zone associated with the building location (Figure 20). Considering the above considerations, it would be expected that the ratio could decrease (indicative of less slender columns) with the increase of the level of seismic loads, being assessed in terms of seismic regulation (year of project) or seismic intensity (seismic zone). Nonetheless, the results presented in Figure 20 indicate that the ratio is essentially independent of the year of construction and, contrary to what would be expected, it appears to increase with the increase in the seismic loads. It is noted that, in the Portuguese seismic zonation, the seismic intensity decreases with the increase in the seismic zone. These conclusions seem to indicate that the seismic combination of loads was not the conditioning action or that, at the time of design, no seismic regulation was available. The latter, however, seems less plausible given that most of the collected data refers to buildings designed after 1990 (see Figure 7) and, therefore, after the introduction of the first reference seismic regulation, the Regulamento de Segurança e Ações para Estruturas de Edifícios e Pontes . Furthermore, the increase in buildings height over time (for higher versatility and to satisfy different industry requirements) was certainly somehow compensated with the optimization of the structural design, construction processes, and the use of higher quality materials, allowing reduction of the elements cross-section requirements.
Finally, the properties of the columns are evaluated in terms of longitudinal and transverse reinforcement ratios (see Figure 21). The results indicate that the mean value of the longitudinal reinforcement ratio is approximately 1.5%, while for the transverse case, the average value approaches 0.3%. These values are somehow higher that the ones observed in conventional RC structures without seismic design that, according to Furtado, Costa, Arêde, and Rodrigues  and Sousa, Costa, Costa, Romão, and Candeias , are about 50% of the ones measured in the precast buildings. This low level of transverse reinforcement results from the fact that, up to 2019, the concrete design codes (REBA  and REBAP ) define the spacing of the stirrups based on the diameter of the longitudinal bars and size of the elements and never related with the seismic zone.
As for the height-to-column ratio, it is also difficult to find a trend between the reinforcement ratios and the year or seismic zone where the buildings are located. Nonetheless, it is apparent there is a slight increase in both ratios with the year of the project (see Figure 22). The apparent increase in both ratios with the seismic zones is somehow justified by the fact that the buildings collected in higher seismic zones (lower seismic hazard) have been built later and, therefore, reflect the tendency for an increase in reinforcement ratios with the year of construction (see Figure 7(b)).
3.5. Beam-to-Column Connection
Past earthquakes showed that beam-to-column connections represent one of the main sources of damage in precast structures [11–14]. The adequate characterisation of these elements is therefore critical to study the structural behaviour of these structures.
In addition to the friction force developed at the surface between the beam and the columns corbel, the lateral strength of these systems depends on the number and size of connecting dowels. Among the analysed design projects, in nearly 60% of the cases, it was possible to access the details about the dowel connection. In the remaining cases, however, no reference to these elements was found, which may indicate that, in a reasonable amount of buildings, the beam-to-column connections could be ensured simply by friction. Considering the limited information often available in older projects, it is also possible that in some cases a standard dowel connection could be assumed during the construction phase, without specifying the details in the design project.
Regarding the cases in which the dowels were detailed, the variability is significant in terms of both the number and diameter of the dowels (Figure 23). It is, thus, not surprising that the total dowel area found in each beam-to-column connection follows a distribution that is essentially uniform (Figure 24(a)).
Looking for a reliable relation between the dowel properties with structural or nonstructural parameters, it was found that the area of the dowels appears to be correlated with the beam span, notwithstanding the large dispersion observed (Figure 24(b)).
Regarding the evolution of the dowel properties with the year of construction, the results presented in Figure 25(a)) show that, only in the recent projects (essentially after 2000), the dowel connections are described in detail. As noted before, these results regard the description found in the design projects, which does not necessarily mean that no dowels were considered in the construction phase.
On the other hand, it is apparent that the dowels properties are not dictated by seismic provisions, since the area of dowels is not proportional to the design seismic action (Figure 25(b)); it is recalled that the seismic loads decrease with the increase in the seismic zone number.
Regarding the corbel properties, the results obtained from the data collected reveal values that vary between 10 and 50 cm (Figure 26). This figure also shows that, contrary to the dowel area, the corbel span seems to be independent of the longitudinal beam span.
3.6. Summary and Statistical Indicators
This section presents an overview of the data collected from 73 precast RC buildings, summarized in Table 3, highlighting the main findings together with the statistics derived for some of the properties.
Regarding the representativeness of the actual building stock, it was showed that the geographical location of the buildings analysed follows closely the actual distribution of industrial buildings. From a temporal point of view, most of the information collected refers to building built after 1990 and, therefore, after the introduction of the seismic regulation . In theory, this could introduce a bias in the information gathered. When looking to the Italian example, around 20% of the precast buildings were constructed before 1980 and around 60% before 1996, without seismic concerns . However, considering that in Portugal the prefabrication industry started to gain a relevant scale in the construction of industrial facilities especially during the 80s and 90s—the National Association of Prefabricated Concrete Industry (ANIPB) was founded in 1975—the asymmetry in the year of the project of the buildings (Figure 7) is certainly a reflection of the evolution of the prefabricated construction and, by that, the bias is not so significant. Moreover, the regressions carried out appear to indicate that the structural properties are independent of the seismic zone and, hence, the design seismic action.
4. Comparison of Portuguese Database with Other European Countries
In this section, the properties identified in the Portuguese building stock are compared with similar data collected in Italy and Turkey. Considering the damage observed in precast structures after recent earthquakes in these countries, the comparison of the properties of the collected buildings in view of the properties found in the other countries can provide indications of the expected behaviour of these structures if subjected to an earthquake in Portugal.
The comparisons made are naturally conditioned by the properties analysed in different studies. Yet, it was possible to establish comparisons regarding the structural typology, main dimensions, and columns properties, including longitudinal and transverse reinforcement ratios.
Regarding the data corresponding to Italy, it comprises the contribution of 3 different studies, namely, one carried out in Emilia-Romagna and Toscana regions and one comprising information collected from different regions promoted by the consortium ReLUIS and the Department of Civil Protection in Italy, hereinafter referred to as DPC/ReLUIS. Additional information regarding the studies consulted can be found in the work carried out in  for Italy and by Senel and Kayhan (2010), for Turkey.
The first parameter addresses the date of construction of the buildings. This information may be relevant in the sense that the properties of the buildings can change as a function of code requirements, material properties, construction process, etc. The results presented in Figure 27 show that the buildings analysed in Portugal are more recent that the ones in Italy and Turkey. The organization of the data around the year 1998 is related with the information available in the Turkish database, which adopted this threshold associated with the implementation of a new code that establishes higher seismic loads .
Regarding the main building typology, according to Casotto, Silva, Crowley, Nascimbene, and Pinho  and , the most common typology found in Italy and Turkey is consistent with the one found in Portugal that essentially consists of 1-storey portal frames with beams simply supported on cantilever column fixed at the base. Regarding the beams cross-section, the results shown in Figure 28 indicate that those are also consistent, at least in regard to the Italian ones (data relative to this property was available only in the study carried out in Emilia-Romagna).
Regarding the beam span in the main direction, Figure 29 shows that the difference between the countries is not significant. In fact, the apparent larger values measured in the Portuguese buildings can be, at least, partially related to younger age of these buildings with respect to the other databases. As noted in Figure 13 and in , the beam spans in the main direction tend to increase with the year of construction.
As regards the columns properties, Figure 30 indicates that the columns height in Portugal appears to be higher than the ones measured in Turkey. On the other hand, the ratio between the height and the length of the columns, which is a measure of their slenderness, is in line with the data collected in Italy (Figure 31).
Finally, regarding the reinforcement ratios, different trends were observed whether related with the longitudinal or transverse reinforcement ratios. In regard to the longitudinal ratio, Figure 32 shows small differences between the Portuguese and Turkish buildings. In one side, values lower than 1% were observed only in Portugal; it is also apparent that a larger percentage of buildings feature large reinforcement ratios (e.g., higher than 2%). On the other hand, based on the results present in Figure 33, it seems clear that the Portuguese buildings feature lower levels of transverse reinforcement, notwithstanding the significantly younger age of the buildings included in the Portuguese database with respect to the Turkish ones (see Figure 27). This lower ratio may anticipate lower ductility capacity (lower confinement level) and higher shear vulnerability of Portuguese buildings when subjected to seismic loads.
5. Final Comments
The determination of the seismic vulnerability of buildings requires the detailed knowledge of their structural properties. For this purpose, a survey was carried out considering a total of 73 design projects of existing buildings. Based on the information collected, it was possible to characterise, from a statistical point of view, the most relevant structural properties of the buildings in Portugal.
Regarding the representativeness of the collected information, the localization of the buildings consulted follows, in a reasonable manner, the actual distribution of the industrial buildings stock. In terms of period of construction, to the authors’ knowledge, there is no information available in the literature to assess if the observed concentration in the period after 1990 is consistent with the age of the existing building stock.
Based on the analysis of the information gathered, the most common typology consists of 1-storey framed buildings with beams simply supported on cantilever columns fixed at the base. Generally, in one of the directions, the number of spans is lower and features longer lengths, achieved through beams with variable cross-section. In the other direction, the number of spans is typically higher but features lower lengths; in this case, the beams present constant cross-section, often used to collect the water from the roof. In a small number of cases (approximately 25%), the buildings feature an intermediate storey that occupies up to 40% of the total area.
In terms of material properties, an increase in both concrete and steel strength was observed with the year of construction. Overall, the steel grades A400 and A500 appear as the most used reinforcement, whilst for the concrete a concentration of buildings employing concrete with compressive strength (associated with cylinder tests) was observed between 20 and 30 MPa.
Regarding the properties of the columns, in addition to the height, cross-section dimensions, and reinforcement ratios (longitudinal and transverse), parameters indicative of the structural robustness such as the height-to-length and length-to-width ratios are also analysed. Additional characteristics relevant from a seismic point of view are also discussed in this document, namely, the type of cladding systems and beam-to-column connection details, including corbel spans and dowel properties.
In summary, the outcome of this study identifies the most conventional typologies of precast RC buildings in Portugal as well as the statistics associated with the most relevant structural parameters. The information provided is essential to estimate the seismic vulnerability of precast RC buildings and consequently estimate the seismic risk associated with a relevant portion of the industrial building park. Finally, the properties of the Portuguese buildings were compared with other studies carried out in other Mediterranean countries, namely, in Italy and Turkey. This comparison showed some consistency regarding the most conventional typologies. However, it was observed that the Portuguese buildings appear to be relatively taller with longer beam spans and, perhaps more importantly, present a relatively lower transverse reinforcement ratio in the columns.
Despite being merely indicative and considering the poor behaviour of precast buildings in recent earthquakes in Italy and Turkey, this comparison raises some questions regarding the seismic safety of the RC precast building stock in Portugal. Additional studies are thus required to assess in more detail the seismic vulnerability of these structures. The statistical data made available in the present work is of paramount importance to realistically characterise the Portuguese building stock and conduct more detailed seismic risk studies.
The data obtained are based on real projects. The statistical analysis reported can be used; however, the project details are not available due to data protection issues.
Conflicts of Interest
The authors declare that they have no conflicts of interest that could have appeared to influence the work reported in this paper.
This work was financially supported by Project POCI-01-0145-FEDER-028439-“SeismicPRECAST: Seismic Performance ASSessment of existing Precast Industrial Buildings and Development of Innovative Retrofitting Sustainable Solutions” funded by FEDER funds through COMPETE2020-ProgramaOperacionalCompetitividade e Internacionalização (POCI) and by national funds (PIDDAC) through FCT/MCTES. The fourth author acknowledges FCT (Fundação para a Ciência e Tecnologia) for the Ph.D. grant with reference SFRH/BD/139723/2018.
- A. Braconi, “Prefabricated steel structures for low-rise buildings in seismic areas (Precasteel),” Brussels, 2013, https://op.europa.eu/en/publication-detail/-/publication/fd7658c8-192e-4eea-ad6d-da916abd4af0.
- D. Bellotti, “Capannoni monopiano prefabbricati: distribuzione proba- bilistica dei sistemi e sottosistemi strutturali dagli anni sessanta ad oggi single-storey precast buildings: probabilistic distribution of structural systems and subsystems from the sixties,” Progett Sismica, vol. 5, no. 3, pp. 41–70, 2014.
- S. M. Senel and A. H. Kayhan, “Fragility based damage assessment in existing precast industrial buildings: a case study for Turkey,” Structural Engineering and Mechanics, vol. 34, no. 1, pp. 39–60, 2010.
- FFMS, “Número de empresas de industrias transformadoras,” Pordata, 2017, https://www.pordata.pt/.
- C. Chastre and V. Lúcio, “Estruturas pré-moldadas no mundo-aplicações e Comportamento Estrutural,” Fundação da Faculdade de Ciência e Tecnologias da Universidade NOVA de Lisboa, Lisbon, Portugal, 2012.
- RSA, “Regulamento de segurança e acções para estruturas de edifícios e pontes-decreto-lei n.o 235/83,” 1983.
- A. Furtado, C. Costa, A. Arêde, and H. Rodrigues, “Geometric characterization of Portuguese RC buildings with masonry infill walls,” European Journal of Environmental and Civil Engineering Impact Factor, vol. 20, no. 4, pp. 396–411, 2015.
- R. Sousa, A. Costa, A. C. Campos Costa, X. Romão, and P. Candeias, “Caracterização do comportamento sísmico de edifícios de betão armado representativos do edificado português sem dimensionamento sismorresistente,” Revista Portuguesa de Engenharia de Estruturas, pp. 105–114, 2017.
- REBA, Regulamento de Estruturas de Betão Armado. Decreto-Lei N.° 47723, REBA, Lisboa, Portugal, 1967.
- REBAP, Regulamento de Estruturas de Betão Armado e Pré-Esforçado. Decreto-Lei N.° 349-C/83, REBAP, Lisboa, Portugal, 1983.
- A. Belleri, E. Brunesi, R. Nascimbene, M. Pagani, and P. Riva, “Seismic performance of precast industrial facilities following major earthquakes in the Italian territory,” Journal of Performance of Constructed Facilities, vol. 29, no. 5, p. 04014135, 2015.
- G. Magliulo, M. Ercolino, C. Petrone, O. Coppola, and G. Manfredi, “The Emilia earthquake: seismic performance of precast reinforced concrete buildings,” Earthquake Spectra, vol. 30, no. 2, pp. 891–912, 2014.
- D. A. Bournas, P. Negro, and F. F. Taucer, “Performance of industrial buildings during the Emilia earthquakes in Northern Italy and recommendations for their strengthening,” Bulletin of Earthquake Engineering, vol. 12, no. 5, pp. 2383–2404, 2014.
- D. Beilic, C. Casotto, R. Nascimbene, D. Cicola, and D. Rodrigues, “Seismic fragility curves of single storey RC precast structures by comparing different Italian codes,” Earthquakes and Structures, vol. 12, no. 3, pp. 359–374, 2017.
- M. Deyanova, S. Pampanin, and R. Nascimbene, “Assessment of silgle-storey precast concrete industrial buildings with hinged beam-column connections with and without dowels,” in Proceedings of the Second European Conference on Earthquake Engineering and Seismology, pp. 1–14, Istanbul, Turkey, August 2014.
- C. Casotto, V. Silva, H. Crowley, R. Nascimbene, and R. Pinho, “Seismic fragility of Italian RC precast industrial structures,” Engineering Structures, vol. 94, pp. 122–136, 2015.
Copyright © 2020 H. Rodrigues et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.