Preprint
Article

Α Rapid Seismic Damage Assessment (RASDA) tool for RC Buildings based on an Artificial Intelligence Algorithms

Altmetrics

Downloads

122

Views

71

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

25 March 2023

Posted:

30 March 2023

You are already at the latest version

Alerts
Abstract
In the current manuscript, a novel software application for Rapid Damage Assessment of RC buildings subjected to earthquake excitation is presented based on Artificial Neural Networks. The software integrates the use of a novel deep learning methodology for Rapid Damage Assessment into modern software development platforms, while the developed graphical user interface promotes the ease of use even from non-experts. The aim is to foster actions both in the pre- and post- earthquake phase. The structure of the source code permits the usage of the application either autonomously as a software tool for Rapid Visual Inspections of buildings prior to or after a strong seismic event or as a component of Building Information Modelling systems in the framework of digitizing building data and properties. The methodology implemented for the estimation of the RC buildings’ damage states is based on the theory and algorithms of Pattern Recognition problems. The effectiveness of the developed software is successfully tested using an extended, numerically generated database of RC buildings subjected to recorded seismic events
Keywords: 
Subject: Computer Science and Mathematics  -   Artificial Intelligence and Machine Learning

1. Introduction

High seismicity observed in many regions globally has always been a critical concern for the safety and the prosperity of their communities. To avoid losses caused by the occurrence of strong earthquake events researchers have been stimulated to develop solutions, such as seismic damage assessment methods, that provide increased and justified knowledge on the structural safety of the built environment exposed to seismic hazard. The application of seismic damage assessment methodologies mainly results in the evaluation of the expected damage level for a range of seismic intensities. This outcome could be either used in the pre- or post- earthquake assessment phase. In the pre-earthquake assessment phase such methods are used to verify the design effectiveness and evaluate the resistance level against the expected seismic intensity for existing structures contributing to retrofit prioritization and planning optimum retrofit schemes. In the post-earthquake assessment phase, such results could be used to rapidly estimate the actual damage level of buildings after earthquake events in aid to first responders and civil protection authorities.
Seismic damage assessment methods can be classified according to the anticipated complexity and the required computational effort. Several methods for Rapid Visual Inspection (RVI) (or Rapid Visual Screening (RVS)) are proposed [1,2,3], based on survey inspection and collection of data, thus not involving analyses and material strength justification. RVI methods are mainly based on visual observation and identification; therefore, they can be widely used for large-scale assessment (i.e., for entire cities or municipalities). However, in this manner, damage evaluation and assessment results strongly depend on expert judgment, thus affecting their effectiveness and reliability. The RVI methods are applicable for both pre-and post-earthquake assessment of structures contributing to the prioritization of interventions and development of emergency plans. It should be outlined that despite the inherent simplifications and assumptions, the RVI methods are valuable since they can be applied in a short time after an earthquake event, providing a first-stage evaluation of the damage extent and highlighting the need for detailed assessment and retrofit measures. The reliability of the assessment results is increased when RVI methods are used in conjunction with results of non-destructive tests and simplified calculations, while the highest level of reliability is achieved when the assessment is based on analysis of the detailed structural model highlighting the retrofit needs, also providing an optimum retrofit scheme and prioritization of interventions.
In line with the above, it is challenging to develop a framework, and a reliable tool for RVS of structures that could overcome the inherent disadvantages (mainly related to expert judgment), will be based on easily (in-situ) collected data, and require low computational effort. To this end, the current paper proposes a novel framework and a software application that exploits modern techniques and Machine Learning (ML) algorithms and is capable of extracting, in real-time, an estimation of the seismic damage assessment of individuals or groups of RC buildings either in the pre- or post- seismic phase. More specifically, the developed software tool is based on the expression of the seismic damage assessment as a Pattern Recognition (PR) problem [4,5]. The first step of RVI methods, either in the pre- or post-seismic phase, is the classification of structures into classes accounting for their seismic damage extend. This approach is compatible with the scope of the PR problems, which are dedicated to the classification of objects of the same type into pre-defined classes. The basic idea for the framework and software developed is to relate the macroscopic structural parameters (collected in the framework of RVI) and the seismic input parameters (based on selected accelerograms) to the expected seismic damage extent that is qualitatively and quantitively defined through Seismic Damage Classes (SDC). Therefore, a preliminary rating of a structure is possible, consisting a valuable tool for rapid assessment and retrofit prioritization for both pre- and post-earthquake assessment. Regarding the computational framework, the Multilayered Feedforward Perceptron Neural Networks (MFPNN) were selected to solve the PR problem. Several research studies have highlighted the effectiveness of this type of Artificial Neural Networks (ANN) within the last three decades [6,7,8,9,10]. In addition, it should be noted that the use of ANNs as a computational tool in software development for the estimation of seismic damage assessment of structures has been proposed in the last years [11,12].
The software developed and presented herein includes a user-friendly and powerful Graphic User Interface (GUI), enabling the quick and easy insertion of the required structural and seismic input data. The features available include the graphical representation of useful structural and seismic input data parameters and the calculation report of relevant indexes defining the RC buildings’ Seismic Damage State (SDC), thus enabling rapid damage assessment. The development of the software using modern development platforms [13,14] and the structure of the corresponding source code enables its wide implementation for pre or post -earthquake seismic damage assessment of existing RC buildings. The software developed may form an independent module used locally (i.e. in desktops, laptops, or portable devices), or it can be included in integrated systems, such as “Building Information Modelling” (BIM) systems, which are the modern trend for the optimum design and management of structures, based on digitization of building data and properties. Due to the fact that BIM [15,16,17,18] is based on integrated systems covering all the stages of buildings’ life cycle, seismic damage prediction indices could be also incorporated. The inclusion of seismic damage assessment procedures in BIM systems has been proposed and investigated in several research studies and state-of-the-art papers [19,20,21,22]. The developed software includes a series of modern techniques which characterize and represent the next generation of the software application, dedicated to the rapid yet reliable and consistent seismic damage assessment in the pre- or post- earthquake phase of individuals or groups of RC buildings. The proposed framework using already trained MFPNN [23,24] and the relevant software developed is successfully applied to selected case studies with either full or limited availability of the required structural property values. The results are presented and discussed herein, highlighting the effectiveness of the tool developed for the rapid seismic damage assessment of RC buildings using available RVI data.

2. Methodologies for Rapid Visual Inspection (RVI) of Buildings

Many methodologies are available worldwide for either pre- or post-earthquake rapid visual inspection of structures, while the most popular and widely used is the FEMA method for Rapid Visual Screening (RVS) of Buildings [1,2]. The New Zealand National Society for Earthquake Engineering [25], the Japan Building Disaster Prevention Association [26], and the GNTD [27] are also widely applied methodologies for Rapid Visual Screening. The RVS procedure prescribed by FEMA uses a methodology based on the building survey inspection, and it requires the filling of a data collection form based on visual observation of the building from the exterior and, if possible, the interior. The two-page data collection form includes building identification information (i.e., usage, area, floor number, etc.), a photograph of the building, sketches, and documentation of pertinent data related to seismic performance.
Simple survey procedures for Seismic Risk Assessment are proposed and applied to urban building stocks to provide damage statistics [28]. In most cases, the expected damage extent is classified as slight, moderate, major, and collapse according to the indicative qualitative definitions of Table 1, based on FEMA.

3. The Proposed Method for the Rapid Seismic Damage Assessment of RC Buildings

The rapid seismic damage assessment method proposed is based on the formulation and solution of a Pattern Recognition (PR) problem. The PR is one of the problems which can be solved using Machine Learning (ML) algorithms [4,5]. For this reason, Artificial Neural Networks (ANN) and, more specifically, the Multilayered Feedforward Perceptron Neural Networks (MFPΝN) are implemented. The theory and the applications of the MFPΝN are described in detail in several books (e.g. [29,30,31]).
The first step towards the development of a methodology for the damage assessment of RC buildings using PR and MFPNN is the definition of the relevant input and output parameters, considering the basic principles of the PR problems and the solution procedure using MFPNN. By definition, PR is the procedure for detecting and classifying objects of the same type into specific classes. A PR problem can be defined using three different approaches: the supervised learning approach, the unsupervised learning approach, or the reinforcement learning approach [31]. The software application developed herein applies the supervised learning approach. Therefore, the target of the applied algorithm is the consistent classification of objects into pre-defined classes. For this reason, the creation of an appropriate data-set and the training of a properly configured MFPNN using this data-set are required. The general configuration of a MFPNN, that can be used to solve a PR problem with n pre-defined classes implementing the supervised learning approach, is depicted in Figure 1. The population of the input parameters, i.e., the parameters that describe the objects (elements of the input vectors x) is m. Thus, in this case, the target is the classification of objects described using the x vectors (having m parameters) into n predefined classes.
In order to use PR techniques and ML algorithms for the rapid seismic damage assessment, the following should be applied:
(a)
The type and number of the characteristics of the objects (i.e., the input x vectors) should be selected and classified into pre-defined classes. Furthermore, the corresponding classes' number and type should also be defined (i.e., the output o vectors). An important issue, which is also briefly described in Figure 1, is the mapping of output vectors (oj elements) to the classes of the problem. More specifically, the configuration of o vectors for each one of the pre-defined classes must be defined. By definition, when the MFPNN extracts an output vector o with oj=1 and all the other elements of o are equal to 0, then the corresponding object (x vector) is classified to the class j. For the seismic damage assessment, the input vectors x should contain parameters that are crucial for the seismic performance of the RC buildings (structural parameters), as well as parameters that describe the seismic excitation (seismic parameters) [32]. To this end, the input vectors x consists of two sub-vectors, namely the sub-vector xstruct and the sub-vector xseism (Equation (1)):
Preprints 70171 i001
(b)
The structural and the seismic parameters should be selected. Since the target of the current paper is the development of the software application based on the PR approach for the rapid seismic damage assessment of RC buildings, no further investigation regarding the optimum selection of the structural and the seismic parameters was performed. Besides, as explained in the next section, the structure of the software code allows for modifications in the case of MFPNNs with different input parameters (structural and seismic). At this stage of development, already trained and successfully tested MFPNNs were used, while description of the parametric investigation for the optimum configuration of these MFPNNs is given in [24]. A brief description of the selected input and output parameters of the introduced MFPNNs as well as their configuration parameters, are presented herein. Regarding the selection of the structural parameters, four parameters were selected with a view to considering parameters that are critical for the seismic performance and the ones are also considered in the framework of taxonomy and classification systems proposed for fragility assessment. The parameters selected are the total height of the building (Htot), the ratio of the base shear that is received by RC walls along two perpendicular between them axes x and y (ratio nvx and ratio nvy), and the structural eccentricity e0 (i.e. the distance between the mass center and the stiffness center of storeys). However, it must be noted that the estimation of the parameter’s values could be difficult in case of RVI. For this reason, as it will be presented in the section 4.4, the developed software allows the input of user-defined parameters, provided that the values are known from previous studies or measurements. In case the values are unknown, a parametric investigation is automatically performed (considering a realistic range of values for specific input parameters) to account for the effect of their variation on the classification of the examined building. This feature renders the software applicable to buildings with both reliable known and unknown (or non-reliable known) structural properties; however, the parametric investigation introduces an inherent uncertainty that is generally acceptable in the framework of the rapid seismic damage assessment methods. The selected seismic parameters are well-documented parameters for the description of the seismic excitations [33,34], widely used in several research studies. The seismic parameters used herein are summarized in Table 2.
The final form of the input vectors x is given in Equation 2.
x = [ x s e i s m x s t r u c t ] T x s e i s m = [ P G A P G V P G D I a S E D C A V A S I H I E P A P G V / P G A P P U D B D S D ] T x s t r u c t = [ H t o t e 0 n v x n v y ] T
(c)
The Seismic Damage Classes (SDC) should be qualitatively and quantitively defined considering appropriate Engineering Demand Parameters (EDPs) and relevant threshold values. The EDPs (which in the present case are also defined as Seismic Damage Indices, (SDI)) could be either global or local [35]. Threshold values should be defined to highlight damage initiation for the limit state considered. The Maximum Interstorey Drift Ratio (MIDR) which is a SDI that refers to buildings’ global performance is selected herein as EDP. Several (5 to 3) Damage States (DS) – which are mapped to SDC – are proposed in the literature for RC buildings, as well as the relevant threshold values, as presented in Table 3 [36]. The MFPNNs used within the software code were trained using the Maximum Interstorey Drift Ratio (MIDR). The SDC considered in the framework of the proposed approach are three, having the threshold values shown in Table 3.
As already stated, the MFPNNs used in the current version of the developed software were trained using three SDC. This approach is more compatible to RVI methods, where the characterization of the seismic damage of buildings is not detailed. The MFPNNs were trained using three different training data-sets created for three RC building classes, classified according to masonry infills existence. In particular, the classes considered are (a) buildings without masonry infills or with light masonry infills (Bare Buildings, BB), (b) buildings with masonry infills at all stories (Masonry Buildings MB), and (c) buildings with masonry infills at all stories except the ground storey (Buildings with Pilotis, PB). More details about the training procedure of these networks are given in [24]. Based on the above, the current version of the software has built-in three trained MFPNNs (see also section 4.2.4). The optimum configuration of the used MFPNNs is presented in Figure 2.

4. Description of the Developed Software

4.1. General Description

A detailed description of the software application is presented in the current section. The software was developed on the MATLAB platform [13,14] and is currently provided as an executable file, while the development of a relevant online tool in Python is currently underway. The main advantages of the MATLAB and the Python software development platforms are their tools for the programming of effective and functional GUIs, as well as the inclusion of Machine Learning (ML) algorithm libraries required for the software proposed. An overview of the software application is presented in Figure 3. A comprehensive and detailed GUI is developed, including all the necessary features for a prompt and efficient operation (described in section 4.2.). The software was developed on the basis of event-driven programming; therefore, the user controls all the procedures through appropriate functions that are activated using the GUI’s control components. More details about these control components are presented in sections 4.2 and 4.3. A detailed description of the procedures is provided in Figure 3.
Τhe developed software provides (a) prediction of the expected DS (SDC) of RC buildings subjected to seismic ground excitations (considering two horizontal components of the corresponding seismic records) using properly trained ANNs which can be alternatively added by the user and (b) estimation the buildings’ seismic damage rating using the First Stage Rapid Visual Screening (FSRVS) procedure, which is practically an RVI methodology proposed by Earthquake Planning and Protection Organization of Greece (E.P.P.O.) for pre-earthquake assessment [37]. Furthermore, the software provides (as intermediate results) the parameters that describe the selected seismic excitation (spectra and seismic parameters) and properly selected structural parameters that are critical for the seismic damage assessment of a studied RC building.
The software is outlined in the sections below. In particular :
-
Description of the GUI and its components used for the insertion of the required data and the presentation of the results.
-
Description of the source code, i.e., the functions used and the interaction between them, along with the required data processing and the corresponding flowcharts.

4.2. Description of the GUI’s components

The proposed application is controlled through a basic/main window organized in four panel containers (Figure 4). These panels correspond to four different procedure categories:
  • Panel 1 entitled “INPUT DATA” contains the GUI tools/components for the structural and seismic data input.
  • Panel 2 entitled “SEISMIC PARAMETERS and SPECTRUMS” includes the presentation of the spectra and the seismic parameters of the selected input seismic excitation.
  • Panel 3 entitled “STRUCTURAL PARAMETERS” includes the illustration of the elastic spectrum of the RC building studied, and the structural parameters that are crucial for the seismic damage assessment in the framework of RVI.
  • Panel 4 entitled “RESULTS” presents the final results of the software, i.e., DS (SDC) prediction (extracted by the selected ANN) for the RC building subjected to the input seismic excitation, as well as its seismic damage assessment and prioritization based on the rating system proposed within the FSRVS procedure.

4.2.1. Panel “INPUT DATA”

The Panel “INPUT DATA” contains the GUI tools for the seismic and structural data insertion. Figure 5 presents the procedures performed by the GUI tools of the panel.
The two push buttons entitled “FILE” activate a window with a File Selector for the selection of the accelerograms of the earthquake excitation (in text file format) in directions 1 and 2 (Figure 6).
The push button entitled “INPUT OF THE STRUCTURAL DATA” activates a window that contains the appropriate GUI control tools/components (i.e., “static text”, “edit text” and “popup menu”) for the insertion of the required structural data (Figure 7(a)). It should be outlined that the inserted structural data are used for the formation of the input vector for the trained ANNs and the DS prediction as well as for the seismic damage assessment and prioritization of the RC building studied, in line with the FSRVS rating system proposed by E.P.P.O. [37].
Using the push button entitled “IMPORT and SHOW DATA”, the window of Figure 7(b) is activated. This window contains the structural data inserted by the user. Thus, an additional check for erroneous or omitted data can be performed. Towards this check the values of inserted data are presented using blue fonts and the omitted data using red hyphens. If the user confirms that all the inserted data are correct, the push button “OK” finishes the procedure. Then the two windows of Figure 7 are closed, and the inserted data are stored in a matrix for the further calculations. In case that errors are detected, the push button “CORRECTIONS” is available for the relevant corrections. Then, the user should go back to the window of Figure 7(a) and proceed with the corrections.

4.2.2. Panel “SEISMIC PARAMETERS and SPECTRUMS”

This panel is the first among the two panels used to present the intermediate results. More specifically, the panel “SEISMIC PARAMETERS and SPECTRUMS” is related to the presentation of the seismic and spectral parameters of the selected seismic excitations. As shown in Figure 8, the panel contains the push button entitled “CALCULATIONS” which activates the connected external function (see section 4.3.) to calculate the acceleration, velocity, and displacement spectra for the two horizontal components of the seismic excitation (directions 1 and 2) and the corresponding (20) selected seismic parameters [33,34]. Diagrams with the calculated spectra are provided within the panel, along with the seismic parameters presented in “static text” GUI components (Figure 8).
It should be noted that the values of the seismic parameters presented in the panel are the geometric means of their values extracted from the input accelerograms along the two horizontal directions. However, the seismic parameters corresponding to the individual accelerograms are also printed in the output file (in .txt or .html format) along with the input data and the other results (Figure 3).

4.2.3. Panel “STRUCTURAL PARAMETERS”

This panel (Figure 9) is the second among the two panels used to present the intermediate results, i.e., the selected structural parameters. In order to run the corresponding external function (see section 4.3.) and depict these parameters, the push button entitled “CALCULATIONS” should be used.
The most indicative representation of this panel is the depiction of (a) the elastic acceleration spectrum of the selected excitation, (b) the elastic design spectrum of the studied RC building, and (c) the predicted range of the building’s uncoupled fundamental elastic eigenperiods, in the same figure (Figure 9(b)). The parameters above are estimated for the two horizontal directions 1 and 2 and are presented in separate figures. Finally, it should be outlined that the equations used for the estimation of the values of the uncoupled fundamental elastic eigenperiods of the studied buildings are well-documented formulas for RC buildings with or without RC shear walls [38,39,40].

4.2.4. Panel “RESULTS”

The panel presented in Figure 10 is the last panel depicting the final results of the application and extracting the output file.
When the push button entitled “PREDICTION OF DAMAGE LEVEL” is used, the window of Figure 11 is activated. This window enables the selection of the trained MFPNN that will be used to estimate the DS (SDC) of the studied RC building for the input earthquake ground motion selected.
Two general options are available:
(a)
Selection of one of the three, embedded in the software application, trained MFPNNs using a pulldown menu. Each one of these MFPNNs was trained considering three general classes of RC buildings regarding the masonry infills. Buildings without masonry infills or with light masonry infills, buildings with masonry infills at all stories, and buildings with masonry infills at all stories except of the ground storey (buildings with pilotis) were used. These embedded MFPNNs were trained and successfully applied as presented in [23,24].
(b)
Selection of an ANN trained by the user. For this purpose, the push button “File” is available, activating a window with a File Selector. The selected file should be compatible with the MATLAB objects containing trained ANNs (“.mat” files). In addition, the imported ANNs should be compatible with the shape of ANNs used in the current version of the software application (section 3).
In Figure 11, detailed instructions for selecting the most suitable ANN for predicting the DS of the studied RC buildings are provided. These instructions are in line with the structural data inserted by the user (Figure 7). For example, suppose that structural parameters inserted by the user for the RC building studied are related with the existence of masonry infills at all stories except from the ground storey (i.e., RC building with pilotis), a notation will then appear in the window, suggesting the implementation of a compatible ANN (i.e., the implementation of an ANN trained using data-set created for RC buildings with pilotis).
After selecting the trained ANN, the software activates the function that simulates this network and extracts the prediction regarding the SDC (i.e. the DS) of the studied RC building (see section 4.3.). This prediction is presented in the properly designed field of the panel (Figure 10). In addition, a short description of the extracted DS is presented in the corresponding static text area (Figure 10). The software also activates the function that performs all the required calculations for the rating of the studied RC building according to the FSRVS procedure. This rating as well as the classification of the studied building in categories with a view to prioritize interventions (i.e., low, medium, or high), are also presented in the panel (Figure 10).
Finally, the panel contains two push buttons, entitled “SAVE and NEW RUN” and “SAVE and EXIT”. Each one activates a different function (see section 4.3.). The first one activates a function that creates an output file (in .txt or .html format) which includes all input data and the results, cleaning the memory for a new run. In contrast, the second one activates a function that also creates an output file with the data provided and the results produced and closes the application. The output could be incorporated even in Building Information Modelling systems.

4.3. Description of the Structure of the Source Code

The source code of the application was developed on the basis of a central/main function that controls the other components/functions of the program (Figure 3). The basic operation of the main function is the programming of the main window and all the other windows of the GUI. In addition, the main function controls all functions of the program, which perform the required calculations, and manages the transferring of the data between them. In other words, the main function controls the insertion of the user data, calls the functions that (using the input data) calculate the intermediate and the final results, and manages their representation. The code of the main function is organized through different code “blocks” which perform specific operations. These “blocks” are summarized in Figure 12.
The functions are separated into two main categories:
(a)
The functions embedded in the main function (internal functions). These functions are used to develop the user interface push buttons, i.e., the operations performed when the user calls them (callback functions). A brief description of these functions is presented in Table 4.
(b)
The external functions, used from the callback functions. These functions perform specific procedures which are mainly computational. The external functions contain the program code, leading to the intermediate and final application results (Figure 3). A brief description of these functions is presented in Table 5.
Figure 13 briefly describes the interconnection between the internal (Table 4) and the external (Table 5) functions. Τhe transmission of data between the functions starting from the data input to the final presentation of the results on the screen and the creation of the final output file are also presented.

4.4. Methodology for the Estimation of Unknown input Structural Parameters

As mentioned in section 4.3, the input vectors of the MFPNN used are formed utilizing the external function “FORM_ANN_INPUT” (Table 5). Moreover, as stated in section 3 (Equation (2)), the input vectors of the MFPNNs embedded in the current version of the software application include 14 seismic (Table 2) and 4 structural parameters (Equation 2).
The seismic parameters are estimated by the processing of the input records (accelerograms) that could be obtained from freely available databases (e.g., PEER strong motion database [41]) or in real time (immediately after an earthquake event) by the stakeholders of accelerometric stations using well-documented formulas [33,34]. These formulas are inserted in the external function “SEISMIC_PARS” (Table 5).
On the contrary, the structural parameters required as input are, in general, difficult to estimate instantly since they are not limited to geometrical parameters (like the total height of buildings Htot, Figure 14) that could be easily measured or defined during RVI. Input parameters like the base shear ratio that is received by RC walls (if exist) along the two principal directions (ratios nv1 and nv2) and the structural eccentricity e0 require time-consuming calculations (no compatible to the RVI methods) in order to be defined.
As already mentioned in section 3 the software application proposed herein can be implemented for the seismic damage assessment of RC buildings independently of reliable data availability. Therefore, the structural parameters of the MFPNNs’ input vectors (Equation (2)) could be reliably known (measured in situ, obtained by technical reports and drawings, or calculated) or completely unknown. In the first case, the user can directly insert the values of the four required structural parameters. In the second case, the user defines the missing parameters as unknown. Reasonable assumptions (based on the available data (Figure 7)) are necessary for the use of the software developed, since it is based on MFPNNs that require the input of the aforementioned structural parameters (Equation (2)). The parameters that are difficult to estimate, using the available structural input data, are the ratios nv1 and nv2 and the value of the structural eccentricity e0.
Based on the above, parametric analyses are proposed in order to apply the software developed consistently and reliably, using the external function “FORΜ_ANN _INPUT” (Table 5). The first approximation of nv1, nv2 and e0 values, is a reasonable, user defined value range (i.e. minimum and maximum values). This value range is estimated accounting for the structural data inserted by the users to the corresponding GUI window (Figure 7).
Regarding the structural eccentricity e0 the minimum (mine0) and the maximum (maxe0) are defined on the basis of the following two structural data inserted by the user (Figure 7): (a) Regularity (or not) of the building in plan, and (b) Regular (or not) distribution of masonry infills in plans. According to the input of the aforementioned structural parameters, three levels of structural eccentricity are defined, namely low, medium and high eccentricity. These levels are related to specific value ranges of e0, as illustrated in Table 6. All values of the structural eccentricities are expressed as percentage of the buildings’ plan dimensions Li (i=1,2, Figure 14), following the rationale of seismic codes such as EN1998-1 [42].
As regards the assumptions for the estimation of the ratios nv1 and nv2, the following three structural data inserted by the user (Figure 7) are considered: (a) The existence (or not) of RC shear walls along the two horizontal directions (axes) x(≡1), y(≡2) if this information is reliably known, (b) the year or the period of the building’s construction, and (c) the application (or not) of seismic code provisions for the design of the building. All the structural data described above are considered to define the maximum and the minimum values of nv1 and nv2. In the current version of the software application, the limit values of nv1 and nv2 were estimated considering the construction practices followed in Greece. However, an expansion of the software source code in order to cover different limit values for nv1 and nv2, is relatively easy. The limit values of nv1 and nv2 inserted in the current version of the software application are presented in Table 7. Different value range is considered according to the period of construction and the seismic code considered for the design (no code, limited requirements, current code provisions). The latter is related to the assumption that relevant code provisions regarding the base shear ratios nv1 and nv2 are applied.
Using the ranges of possible minimum and maximum values of structural parameters nv1, nv2 and e0 (Table 6, Table 7) the eight input vectors for the MFPNN of Equation (3) are automatically formed by means of the external function “FORM_ANN_INPUT”. The combination of all the possible nv1, nv2 and eo parameter values (maximum and minimum), results in eight different input vectors used to define eight different classifications of the examined RC building to the three pre-defined SDC (i.e. to the corresponding DS) of Table 3. The worst classification is presented in the panel “RESULTS” (Figure 10).
x = [ x s e i s m x s t r u c t ] [ x s t r u c t ] 1 = [ H t o t min e 0 min n v x min n v y ] , [ x s t r u c t ] 3 = [ H t o t min e 0 min n v x max n v y ] , [ x s t r u c t ] 5 = [ H t o t min e 0 max n v x min n v y ] , [ x s t r u c t ] 7 = [ H t o t min e 0 max n v x max n v y ] [ x s t r u c t ] 2 = [ H t o t max e 0 min n v x min n v y ] , [ x s t r u c t ] 4 = [ H t o t max e 0 min n v x max n v y ] , [ x s t r u c t ] 6 = [ H t o t max e 0 max n v x min n v y ] , [ x s t r u c t ] 8 = [ H t o t max e 0 max n v x max n v y ]
Finally, it should be outlined that the parametric analysis described above is not mandatory, since user defined input values could also be applied. More specifically, if the user can estimate more reliably the minimum and the maximum values of the unknown structural parameters it is possible to insert these values using the GUI and run separately the software application for each one of them. In other words, in this case the user can perform the required parametric analyses manually.

5. Numerical Applications

The current section applies the software developed for both pre-and post-earthquake seismic damage assessment of RC buildings. Two different cases of data availability are considered, namely full data availability for RC buildings used for the pre-earthquake assessment and limited data availability for the case of post-earthquake assessment. The latter is frequently the case in real-time applications, where the need for rapid estimation of the seismic damage level is evident, along with the difficulties for structural data collection. In this case, the parametric investigation considering several scenarios regarding the values of the required structural parameters of the examined RC buildings that are not known is performed (section 4.4.). The ability of the developed software to instantly provide the classification of an RC building to the pre-defined DS, based on the use of the embedded MFPNN, makes the parametric investigation possible, even in real-time analyses.

5.1. Software Application for the Pre-Earthquake Assessment of RC Buildings with Full Data Availability

In this case, the software was applied for 30 RC buildings with reliably known structural parameters presented in Table 8 and explained in Figure 14. These buildings were also analysed through Non-linear Time History Analyses (NTHA), which is the most effective numerical method for estimating the seismic DS (in the present paper using the MIDR seismic damage index). At first, the selected RC buildings were modeled and designed considering the EN1992-1-1 [46] and EN1998-1 [42] recommendations. The buildings, designed with current code provisions, were subsequently assessed by means of NTHA for 65 known seismic excitations extracted from the European Strong-Motion Database [47] and PEER [41]. The ranges of the seismic parameters’ values of the 65 selected seismic excitations are given in [24].
All the selected buildings are rectangular in plan and regular in elevation according to the provisions of EN1998-1. It must be noted that for each one of the 30 selected buildings three different versions were considered regarding their masonry infills: (a) without masonry infills or with masonry infills with weak contribution to the buildings’ seismic response (“Bare Buildings”- “BB”), (b) with strong masonry infills at all storeys (“Fully Infilled Buildings”- “FIB”) and (c) with strong masonry infills at all storeys except of the ground storey (“Pilotis Buildings” – “PB”). The contribution of the masonry infills to the buildings’ seismic response was considered in the framework of NTHA, using the model proposed in [48]. At the design stage, the contribution of the masonry infills was ignored. Thus, the masonry infills were considered only as vertical loads at that stage.
The NTHA performed led to 1950(=30·65) MIDR values. Thus, for the three versions of the 30 selected RC buildings 5850(=1950·3) MIDR values were calculated. Using these values, 5850 classifications to the three pre-defined SDC (i.e. DS) of Table 3 were performed.
The classification of the 5850 generated samples to the three pre-defined SDC using the results of NTHA constitutes the data-set for the examination of the efficiency of the developed software application. A MATLAB compatible script was developed to run the application automatically 5850 times using the data of the 90 selected RC buildings and the 65 selected seismic excitations. The classification results based on NTHA and the relevant derived from the software application were subsequently compared. The comparisons are illustrated in terms of Confusion Matrices (CM) [4,49], which are the most effective tools for the evaluation of the effectiveness of trained MFPNN used for the solution of PR problems.
The CMs for the 5850 samples, corresponding to the 90 known RC buildings (with full data availability, namely reliable knowledge of their structural data), are given in Figure 15.
The performance of the MFPNNs, embedded in the current version of the software application is significantly high. However, the most crucial point in the present case is that the implementation of the MFPNNs can be achieved in almost real time, exploiting the developed GUI capabilities for the input of data for a large number of buildings with reliably known structural parameters.

5.2. Software Application for the Post-Earthquake Real-Time Assessment of RC Buildings with Limited Data Availability

In this case, the software was applied to three case study RC buildings with structural parameters that are considered unknown in order to check the reliability of the parametric procedure based on the proposed assumptions described in section 4.4. In fact, the structural parameters of these buildings are known and presented in Table 9. However, in order to proceed with the reliability check, the structural data of this Table were not directly inserted into the application to investigate the effectiveness of the parametric analyses performed by the external function “FORM_ANN_INPUT” when the structural parameters of the MFPNNs’ input vectors are considered unknown (section 4.4).
The case study buildings were designed according to the provisions of EN1992-1-1 and EN1998-1. The information about the application of these codes for the design of the buildings of Table 9 was considered as known. Three different versions of each one of the three case study buildings were considered, differentiated regarding their masonry infills; i.e., “Bare Buildings”- “BB”, “Fully Infilled Buildings”- “FIB” and “Pilotis Buildings”- “PB” (section 5.1). After the design procedure, the case study buildings were assessed using NTHA for the 65 seismic excitations, also used for the analyses of the 30 known RC buildings presented in the previous section. Based on the analyses results, the calculation of the seismic damage indices (i.e. the MIDR values) and consequently the classifications to the three pre-defined DS of Table 3 were extracted. Since the structural parameters nv1, nv2, and e0 were considered as unknown, the assumptions described in section 4.4. were followed. To this end, the data required to define the range of the possible values of unknown structural parameters for the three case study buildings were estimated and inserted into the application (Table 10).
Considering the structural data of Table 10 and under the assumptions of section 4.4, the unknown values of the structural parameters nv1, nv2, and e0 of the three case study buildings fluctuate into the ranges presented in the Table 11.
The eight input vectors of Equation (3) were formed for each one of the case study buildings and the 65 seismic excitations, using the values of Table 11 as input. Thus, for each one of the three case study buildings, 520(=65 earthquakes·8 input vector types of Equation (3)) input vectors were created, and 1560(=520·3 versions of buildings) input vectors were introduced to the MFPNNs, embedded in the current version of the developed software application. Consequently, the software was automatically applied 4680(=1560·3 case study buildings) times using an appropriate script. In order to check the reliability of the results, the software was also applied, considering the real values of the structural parameters nv1, nv2, and e0 (Table 9).
The Overall Accuracy (OA) index values, i.e., the percentages of the correct classifications of the case study buildings subjected to the 65 selected seismic excitations into the three pre-defined SDC (DS), as calculated by the MFPNNs, are presented in Figure 16. These OA index values correspond to each one of the eight input vectors Xi (i=1-8) of Equation (3) [OA(Xi)] as well as to the input vectors XR which include the actual values of the structural parameters nv1, nv2, and e0 ([OA(XR)] – black bars).
The main conclusion that arises from the study of Figure 16 is that the parametric procedure, based on the consideration of eight input vectors Xi (Equation (3)), can effectively approach the percentages of the correct classifications extracted using the input vectors XR. More specifically, in the case of BB the differences between the OA(Xi) index values and the corresponding OA(XR) index values (i.e., the Δ(OA) values in Figure 16) fluctuate between -1.8% and -18.3%, while the corresponding lower and upper values of differences in the cases of FIB and PB are -0.6%/-9.4% and -1.9%/-19.9% respectively. Bearing in mind that the feasible level of accuracy of the RVI methods is generally low, the abovementioned deviations between the OA(Xi) and the OA(XR) index values can be considered acceptable. However, it must be outlined that the considered value ranges of the unknown structural parameters nv1, nv2, and e0 strongly affect the deviation between the OA(Xi) and the OA(XR) index values. In the present paper, the ranges summarized in Table 11 were used. The software proposed enables the selection of different ranges for the values on the unknown structural parameters on the basis of the user’s expertise. Further investigation of these ranges is possible, but out of the scope of the present paper.
Finally, the CMs corresponding to the classifications performed using the input vectors XR and the input vectors Xi that provides the min[Δ(OA)] and the max[Δ(OA)] of Figure 16 are presented in Figure 17. The study of the CMs extracted using these Xi lead to the conclusion that the vast majority of the Recall (R-index) and the Precision (P-index) indices are close to the corresponding values of the CMs extracted using XR as input vectors. Thus, the classifications of the three case study buildings to the three pre-defined DS extracted using the input vectors Xi are not significantly different from the corresponding classifications extracted using the input vectors XR with the actual values of the structural parameters nv1, nv2, and e0. However, as mentioned above, the approach based on the use of input vectors Xi (if the input vectors XR are unknown) can be significantly improved if the users’ expertise can lead to more reliable estimation of the unknown structural parameters nv1, nv2, and e0 than the one presented in Table 11. To this end, the software is developed open to include user-defined values.

6. Conclusions

The scope of the current paper is to emerge and highlight the functions, the applicability, and the advantages of a new software application developed for the rapid seismic damage assessment of RC buildings. The developed software is intended to instantly provide assessment results and is used for both pre-and post-earthquake assessments. Whereas the required seismic data can be instantly available either from strong-motion databases or in real time after a seismic event, the level of availability of the structural data (parameters) is low due to the limited time and limited resources.
The selected ANN-based method for the seismic damage assessment of RC buildings is a method based on the definition and the solution of the Pattern Recognition (PR) problems. More specifically, by formulating the structures’ seismic damage assessment problem (compatible with “first stage simplified methodologies”) in terms of PR problems and of ANNs’ function, the software application extracts predictions regarding the seismic Damage State (DS) of RC buildings in real-time. This DS is quantified considering three Seismic Damage Classes (SDC) defined using specific threshold values of a Seismic Damage Index (SDI). In the current paper, the Maximum Interstorey Drift Ratio (MIDR), which is extensively used in research studies, was selected as SDI. The selection of three DS was in line with the DS proposed in RVIs, in order to be directly applicable. The seismic and structural parameters were selected as input parameters for the ANNs in the current version of the proposed application; however, they can be easily changed in future versions due to the flexibility of the source code developed.
The current version of the developed software application was evaluated by predicting the DS of two general categories of RC buildings, namely buildings with reliable known values of the selected structural parameters and buildings with unknown (or no-reliable known) values of the selected structural parameters. In both cases, the effectiveness of the proposed application regarding its capability to extract reliable results in real time was highlighted. However, it must be noted that the accuracy of the results extracted by the application depends on the effectiveness of the ANNs used. The latter is not the subject of the present paper since already trained ANNs were used. In line with the above, the usage of user-defined ANNs trained with different databases is also possible.
In any case, the capability for direct (in real time) and reliable estimation of the DS of a building or building inventories in the pre-earthquake phase or after a strong earthquake, avoiding time-consuming modeling and analyses, without exclusively requiring detailed RVI input data (implementing parametric analyses), renders the software developed a valuable calculational tool available to the authorities for either for retrofit prioritization in the pre-earthquake phase or emergency planning after an strong earthquake event.

Author Contributions

Conceptualization, Konstantinos Morfidis; methodology, Konstantinos Morfidis and Sotiria Stefanidou; software, Konstantinos Morfidis and Olga Markogiannaki; validation, Sotiria Stefanidou and Olga Markogiannaki; investigation, Konstantinos Morfidis, Sotiria Stefanidou and Olga Markogiannaki; writing—original draft preparation, Konstantinos Morfidis; writing—review and editing, Sotiria Stefanidou and Olga Markogiannaki; project administration, Konstantinos Morfidis; funding acquisition, Sotiria Stefanidou and Olga Markogiannaki. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Regional Development Fund of the European Union and Greek national funds through the Operational Program Competitiveness, Entrepreneurship and Innovation, under the call RESEARCH CREATE INNOVATE – Second Round (project code: T2EDK-03412).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are only available after contacting the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. FEMA P-154. (2015). Rapid visual screening of buildings for potential seismic hazards: A Handbook (Issue January). [CrossRef]
  2. FEMA 356. (2000). Prestandard and Commentary for the Seismic Rehabilitation of Buildings FEMA-356 (Issue November).
  3. GNDT. Detection of seismic vulnerability of masonry buildings - Instructions for filling in the 2th level form (in Italian). CNR, Rome, Italy (1993).
  4. S. Theodoridis, K. Koutroumbas, Pattern Recognition, fourth ed. Elsevier, 2008. [CrossRef]
  5. C.M. Bishop, Pattern Recognition and Machine Learning, Springer, 2006.
  6. H. Adeli, Neural networks in civil engineering: 1989-2001, Comput Aid Civ Infrastruct Eng 16 (2001) 126–42. [CrossRef]
  7. Shahin, M.A.; Jaksa, M.B.; Maier, H.R. State of the art of artificial neural networks in geotechnical engineering. Electr J Geotech Eng EJGE 2008, 8, 1–26. [Google Scholar]
  8. Jegadesh, S.J.S.; Jayalekshmi, S. A review on artificial neural network concepts in structural engineering applications. Int J Appl Civ Env Eng 2015, 1, 6–11. [Google Scholar]
  9. H. Sun, H.V. Burton, H. Huang, Machine learning applications for building structural design and performance assessment: State - of - the - art review, Journal of Building Engineering 33 (2021) Article 101816. [CrossRef]
  10. E. Harirchian, S.E.A. Hosseini, K. Jadhav, V. Kumari, S. Rasulzade, E. Işık, M. Wasif, T. Lahmer, A review on application of soft computing techniques for the rapid visual safety evaluation and damage classification of existing buildings, Journal of Building Engineering 43 (2021) Article 102536. [CrossRef]
  11. Y. Xu, X. Lu, B. Cetiner, & E. Taciroglu, Real-time regional seismic damage assessment framework based on long short-term memory neural network, Computer-Aided Civil and Infrastructure Engineering, (2020). [CrossRef]
  12. S. Ruggieri, A. Cardellicchio, V. Leggieri, G. Uva, Machine-learning based vulnerability analysis of existing buildings, Autom. Constr. 132 (2021) Article 103936. [CrossRef]
  13. Matlab R2022a, Deep Learning toolbox User’s guide, (2022). https://www.mathworks.com/help/deeplearning/.
  14. Matlab R2022a, App Building, (2022). https://www.mathworks.com/products/matlab/app-designer.html.
  15. X. Wang, P.E. Love, M.J. Kim, C.-S. Park, C.-P. Sing, L. Hou, A conceptual framework for integrating building information modeling with augmented reality, Autom. Constr. 34 (2013) 37-44. [CrossRef]
  16. X. Liu, X. Wang, G. Wright, J.C.P. Cheng, X. Li, R. Liu, A state-of-the-art review on the integration of building information modeling (BIM) and geographic information system (GIS), ISPRS International Journal of Geo-Information 6(2) (2017). [CrossRef]
  17. R. Crotty, The Impact of Building Information Modelling Transforming Construction, first ed. Taylor and Francis, London, UK, 2011. [CrossRef]
  18. R. Garber, BIM Design: Realising the Creative Potential of Building Information Modelling, Wiley, 2012 ISBN: 1322229473. ISBN: 978-1-118-71979-4.
  19. M. Alirezaei, M. Noori, O. Tatari, K.R. Mackie, A. Elgamal, BIM-based damage estimation of buildings under earthquake loading condition, Procedia Eng. 145 (2016) 1051-1058. [CrossRef]
  20. Ζ. Xu, H. Zhang, Χ. Lu, Υ. Xu, Ζ. Zhang, Υ. Li, A prediction method of building seismic loss based on BIM and FEMA P-58, Autom. Constr. 102 (2019) 245–257. [CrossRef]
  21. S. Christodoulou, D. Vamvatsikos, C. Georgiou, A BIM-based framework for forecasting and visualizing seismic damage, cost and time to repair, In: Proceedings of the European Conference on Product and Process Modelling, Cork, Ireland, 14–16 September 2011.
  22. C. Georgiou, D. Vamvatsikos, Damage Assessment, Cost Estimating, and Scheduling for Post-Earthquake Building Rehabilitation Using BIM, In: Proceedings of the 31st International Conference of CIB W78, Orlando, Florida, USA, 23-25 June 2014 398-405.
  23. K. Morfidis, K. Kostinakis, Approaches to the rapid seismic damage prediction of r/c buildings using artificial neural networks, Eng Struct 165 (2018) 120-141. [CrossRef]
  24. K. Morfidis, K. Kostinakis, Comparative evaluation of MFP and RBF neural networks’ ability for instant estimation of r/c buildings’ seismic damage level, Eng Struct 197 (2019) 1–19. [CrossRef]
  25. New Zealand Society for Earthquake Engineering (NZSEE). Assessment and Improvement of the Structural Performance of Buildings in Earthquakes; Recommendations of a NZSEE Study Group on Earthquake Risk Buildings, June 2006; NZSEE: Wellington, New Zealand, 2006. [Google Scholar]
  26. Japan Building Disaster Prevention Association 1990. Standard for evaluation of seismic capacity and guidelines for seismic retrofit design of existing reinforced concrete buildings. 1977 (revised 1990) (in Japanese).
  27. GNDT (1993). Seismic Risk of public buildings – Part 1 – Methodology Aspects (in Italian). CNR, Rome, Italy.
  28. H. Sucuoğlu, U. Yazgan, Simple survey procedures for seismic risk assessment in urban building stocks. In: S. T. Wasti, G. Ozcebe editors. Seismic Assessment and Rehabilitation of Existing Buildings, Vol 29, Earth and Environmental Sciences: Kluwer Academic Publishers, London, (2003) 97-118.
  29. K. Gurney, An Introduction to Neural Networks, UCL Press, 1997. [CrossRef]
  30. L. Fausett, Fundamentals of Neural Networks: Architectures, Algorithms and Applications, Pearson, 1994.
  31. S. Haykin, Neural networks and learning machines, third ed. Prentice Hall, 2009.
  32. de Lautour, P. Omenzetter, Prediction of seismic-induced structural damage using artificial neural networks, Eng Struct 31(2) (2009) 600–606. [CrossRef]
  33. S.L. Kramer, Geotechnical earthquake engineering, Prentice-Hall, 1996.
  34. SeismoSoft. SeismoSignal v.5.1.0; 2014. www.seismosoft.com.
  35. A.J. Kappos, Seismic damage indices for RC buildings: evaluation of concepts and procedures, Construction Research Communications Limited ISSN 1365-0556, (1997) 78-87. [CrossRef]
  36. Masi, M. Vona, M. Mucciarelli, Selection of natural and synthetic accelerograms for seismic vulnerability studies on reinforced concrete frames, J Struct Eng 137 (2011) 367–378. [CrossRef]
  37. First Stage Rapid Visual Screening – Greek Rapid Visual Investigation Methodology (TOE) 5th edition, Earthquake Planning and Protection Organization of Greece (E.P.P.O), 2020. Available at. https://www.oasp.gr/node/74.
  38. H. Crowley, R. Pinho, Revisiting Eurocode 8 formulae for periods of vibration and their employment in linear seismic analysis, Earthq. Eng. Struct. Dyn. 39(2) (2010) 223-235. [CrossRef]
  39. R.K. Goel, A.K. Chopra, Period formulas for moment resisting frame buildings, J. Struct. Eng. 123(11) (1997) 1454–1461. [CrossRef]
  40. R.K. Goel, A.K. Chopra, Period formulas for concrete shear wall buildings, J. Struct. Eng. 124(4) (1998) 426–433. [CrossRef]
  41. PEER (Pacific Earthquake Engineering Research Centre). Strong motion database (2003): https://ngawest2.berkeley.
  42. EN1998-1 (Eurocode 8). Design of structures for earthquake resistance - part 1: general rules, seismic actions and rules for buildings, European Committee for Standardization, (2005).
  43. Royal Decree on the Seismic Code for Building Structures, Government Gazette Issue A No. 36 (1959), (in Greek).
  44. NEAK: The New Greek Antiseismic Regulations, Bulletin of the Technical Chamber of Greece, Νο. 1757, Earthquake Planning and Protection Organization, Athens, Greece (1993), (In Greek).
  45. ΕAΚ/2000: Greek Seismic Code, Earthquake Planning and Protection Organization, Athens, Greece (1999), (In Greek).
  46. EN1992-1-1 (Eurocode 2). Design of concrete structures, Part 1-1: General rules and rules for buildings. European Committee for Standardization, (2005).
  47. European Strong-Motion Database, (2003). http://isesd.hi.is/ESD_Local/frameset.htm.
  48. F.J. Crisafulli, Seismic behaviour of reinforced concrete structures with masonry infills, Ph.D. Thesis, University of Canterbury, Christchurch, New Zealand 1997. [CrossRef]
  49. T. Fawcett, An introduction to ROC analysis, Pattern Recognition Letters 27(8) (2006) 861–874. [CrossRef]
Figure 1. General description of the configuration of a MFPNN used for the solution of PR problems
Figure 1. General description of the configuration of a MFPNN used for the solution of PR problems
Preprints 70171 g001
Figure 2. Configuration parameters of the optimum MFPNNs used in the current version of the software.
Figure 2. Configuration parameters of the optimum MFPNNs used in the current version of the software.
Preprints 70171 g002
Figure 3. Overview of the developed software application.
Figure 3. Overview of the developed software application.
Preprints 70171 g003
Figure 4. Main window of the proposed application (before the insertion of the input data).
Figure 4. Main window of the proposed application (before the insertion of the input data).
Preprints 70171 g004
Figure 5. Brief description of the GUI control tools/components of the panel “INPUT DATA”.
Figure 5. Brief description of the GUI control tools/components of the panel “INPUT DATA”.
Preprints 70171 g005
Figure 6. The file selector activated when the two push buttons “FILE” are pushed by the user.
Figure 6. The file selector activated when the two push buttons “FILE” are pushed by the user.
Preprints 70171 g006
Figure 7. The windows used for the insertion and checking of the structural data.
Figure 7. The windows used for the insertion and checking of the structural data.
Preprints 70171 g007
Figure 8. The window used for the presentation of the seismic parameters (a) before and (b) after running.
Figure 8. The window used for the presentation of the seismic parameters (a) before and (b) after running.
Preprints 70171 g008
Figure 9. The window used for the presentation of the structural parameters (a) before and (b) after running.
Figure 9. The window used for the presentation of the structural parameters (a) before and (b) after running.
Preprints 70171 g009
Figure 10. The window used for the presentation of the final results (a) before and (b) after running.
Figure 10. The window used for the presentation of the final results (a) before and (b) after running.
Preprints 70171 g010
Figure 11. The file selector activated for the selection of the trained MFPNN.
Figure 11. The file selector activated for the selection of the trained MFPNN.
Preprints 70171 g011
Figure 12. General description of the structure of the software application’s main function.
Figure 12. General description of the structure of the software application’s main function.
Preprints 70171 g012
Figure 13. Brief flow chart of the source code of the developed software application.
Figure 13. Brief flow chart of the source code of the developed software application.
Preprints 70171 g013
Figure 14. Description of the structural parameters imported to the input vectors of the trained MFPNN.
Figure 14. Description of the structural parameters imported to the input vectors of the trained MFPNN.
Preprints 70171 g014
Figure 15. CMs of the classifications performed by the MFPPNs embedded in the current version of the software application (a) for BB, (b) for FIB, and (c) for PB.
Figure 15. CMs of the classifications performed by the MFPPNs embedded in the current version of the software application (a) for BB, (b) for FIB, and (c) for PB.
Preprints 70171 g015
Figure 16. OA index values of the classifications extracted by the used MFPNNs considering input vectors with assumed values for the unknown structural parameters: (a) for BB, (b) for FIB, and (c) PB.
Figure 16. OA index values of the classifications extracted by the used MFPNNs considering input vectors with assumed values for the unknown structural parameters: (a) for BB, (b) for FIB, and (c) PB.
Preprints 70171 g016aPreprints 70171 g016b
Figure 17. CMs correspond to the classifications extracted by the used MFPNNs considering input vectors with assumed values for the unknown structural parameters: (a) for BB, (b) for FIB, and (c) PB.
Figure 17. CMs correspond to the classifications extracted by the used MFPNNs considering input vectors with assumed values for the unknown structural parameters: (a) for BB, (b) for FIB, and (c) PB.
Preprints 70171 g017
Table 1. The seismic parameters of the input vectors for the trained MFPNNs of the developed software.
Table 1. The seismic parameters of the input vectors for the trained MFPNNs of the developed software.
Damage State Qualitative description
Slight damage No permanent drift. Structure substantially retains original strength and stiffness. Minor cracking of facades, partitions, and ceilings as well as structural elements. All systems important to normal operation are functional.
Moderate damage Wider cracks at non-structural elements, in-plane or out-of-plane. Some residual strength and stiffness left in all stories. Gravity-load bearing elements function. No out-of-plane failure of walls or tipping of parapets. Some permanent drift. Damage to partitions. Cracking of facades, partitions, and ceilings as well as structural elements. Flexural and shear cracks at structural elements, concrete spalling.
Major damage Local in-plane and out-of-plane failure of nonstructural walls, infills, etc. Wide flexural and shear cracks at structural elements, hoop fracture, buckling of longitudinal reinforcement, initiation of concrete core crushing. Little residual stiffness and strength, but loadbearing columns and walls function. Large permanent drifts. Some exits blocked. Infills and unbraced parapets failed or at incipient failure.
Collapse Loss of load-carrying capacity, locally or globally.
Table 2. The seismic parameters of the input vectors for the trained MFPNNs of the developed software.
Table 2. The seismic parameters of the input vectors for the trained MFPNNs of the developed software.
Seismic Parameter
1 Peak Ground Acceleration - PGA 8 Housner Intensity - HI
2 Peak Ground Velocity -PGV 9 Effective Peak Acceleration - EPA
3 Peak Ground Displacement - PGD 10 Vmax/Amax (PGV/PGA)
4 Arias Intensity - Ia 11 Predominant Period - PP
5 Specific Energy Density - SED 12 Uniform Duration - UD
6 Cumulative Absolute Velocity - CAV 13 Bracketed Duration - BD
7 Acceleration Spectrum Intensity - ASI 14 Significant Duration - SD
Table 3. Definition of the SDC (DS) using the MIDR as EDP.
Table 3. Definition of the SDC (DS) using the MIDR as EDP.
MIDR [%] <0.25 0.25-0.5 0.5-1.0 1.0-1.5 >1.5
SDC (5 classes) Null Slight Moderate Heavy Destruction
SDC (3 classes) Slight (“S”) Moderate (“M”) Heavy (“H”)
Description No (or repairable) structural damages Significant but repairable structural damages Non-repairable structural damages
Table 4. Description of the internal functions of application (built-in in the file of the main function).
Table 4. Description of the internal functions of application (built-in in the file of the main function).
Function / [GUI component – Panel] Procedure
Acc1_Callback / [Pushbutton “FILE” - Panel 1], (Figure 5)
{Acc2_Callback / [Pushbutton “FILE” - Panel 1]}, (Figure 5)
Insertion (through the activation of a window with file selector), saving, and plotting of the accelerogram of the selected excitation in the direction 1 {2}
StrInpData_Callback / [“INPUT OF THE STRUCTURAL DATA” – Panel 1], (Figure 5) Creates and activates a GUI window for the insertion of the structural data (this window activates a second window for the checking of the imported data using the function ChStrDat Callback activated by the push button “IMPORT and SHOW DATA”)
SePar_Callback / [“CALCULATIONS” – Panel 2], (Figure 8) Calculates (using the external function “SEISMIC_PARS”) and presents of the excitation’s acceleration, velocity and displacement spectra and the seismic parameters
StrPar_Callback / [“CALCULATIONS” – Panel 3], (Figure 9)
  • Calculates and plots the excitation’s acceleration spectrum and the design elastic acceleration spectrum (external function “GR_CODE_EL_SPECTR”)
  • Calculates the structural parameters used for the rapid estimation of the building’s seismic damage assessment
Results_Callback / [“PREDICTION OF DAMAGE LEVEL” – Panel 4], (Figure 10) Classifies the building to one of the three pre-defined DS (SDC) and its rating according to the FSRVS procedure proposed by E.P.P.O. [by calling the external functions: “TRAINED_ANN SELECTOR”, “ANN_CALCS” and “FSRVS_CALCS”]
SavEx_Callback / [“SAVE & EXIT” - Panel 4], (Figure 10)
  • Calls the external function “FINAL_OUTPUT” that creates the output file which includes all input data and results,
  • Closes the application
SavNRun_Callback / [“SAVE & NEW RUN” - Panel 4], (Figure 10)
  • Calls the external function “FINAL_OUTPUT” that creates the output file which includes all input data and results,
  • Clears the memory (external function “CLEAR_MEM”) and restarts the application
Table 5. Description of the external functions of the software.
Table 5. Description of the external functions of the software.
Function / Called from callback function Procedure
“SEISMIC_PARS” / SePar_Callback Calculates the seismic parameters of the selected excitation in directions 1 and 2
“GR_CODE_EL_SPECTR” / StrPar_Callback Calculates the design elastic acceleration spectrum (according to the valid seismic codes at the year of the building’s design)
“FORΜ_ANN _INPUT” / Results_Callback Configuration of the vectors used as input of the selected trained ANNs
“TRAINED_ANN_SELECTOR” / Results_Callback Creates and activates a window with GUI tools for the selection and insertion of a trained ANN (Figure 11)
“ANN_CALCS” / Results_Callback Simulates the imported ANN for the rapid prediction of the DS of building
“FSRVS_CALCS” / Results_Callback Required calculations for the rating of building and its classification into the priority classes, defined in the framework of FSRVS method
“FINAL_OUTPUT” / SavEx_Callback and SavNRun_Callback
(a)
Creates and activates a window for the saving of the final output file
(b)
Creates the output file that includes all input data and (intermediate, final) results
“CLEAR_MEM” / SavEx_Callback and SavNRun_Callback Clears the memory from the parameters which are defined during the running of the application
Table 6. Relation assumed between the structural data inserted by the user and the structural eccentricity e0.
Table 6. Relation assumed between the structural data inserted by the user and the structural eccentricity e0.
Case 1 Case 2A Case 2B Case 3
Regularity in plan Yes Yes No No
Regular distribution of masonries Yes No Yes No
Eccentricity Low Medium High
mine0i / maxe0i 0.05Li / 0.075Li 0.075Li / 0.15Li 0.15Li / 0.175Li
Table 7. Relation between the structural data inserted by the user and the ratios nv1 and nv2 (minnv1/maxnv2).
Table 7. Relation between the structural data inserted by the user and the ratios nv1 and nv2 (minnv1/maxnv2).
Period (year) of construction
Existence of RC shear walls Usage of seismic code <1959 1959-1984 1984-1992 1992-2000 2000-2010 >2010
- R/D 1959(1) Expansion of R/D 1959 NEAK(2) EAK/2000(3) EAK/2000 and Eurocodes
No Yes minnv=maxnv=0.0
No
Yes Yes 0.05/0.15 0.1/0.30 0.2/0.40 0.25/0.45 0.275/0.45 0.35/0.65
No 0.025/0.05 0.025/0.05 0.05/0.075 0.075/0.10 0.075/0.10 0.10/0.20
Unknown Yes 0.0/0.15 0.0/0.30 0.0/0.40 0.0/0.45 0.0/0.45 0.0/0.65
No 0.0/0.05 0.0/0.05 0.0/0.075 0.0/0.10 0.0/0.10 0.0/0.20
(1) R/D: Royal Decree on the Seismic Code for Building Structures [43]. (2) NEAK: The New Greek Antiseismic Regulations [44].(3) ΕAΚ/2000: Greek Seismic Code [45].
Table 8. Structural data of the 30 examined known RC buildings used for the prediction of their DS.
Table 8. Structural data of the 30 examined known RC buildings used for the prediction of their DS.
Name nvx nvy Htot (m) Lx (m) Ly (m) ex (m) ey (m)
1 SFxy_3 0.00 0.00 9.60 13.50 10.00 0.00 0.00
2 SFxy_5 0.00 0.00 16.00 20.00 14.00 0.00 0.00
3 SFxy_7 0.00 0.00 22.40 20.00 14.00 0.00 0.00
4 SWxy_3 73.00 76.00 9.60 15.00 10.00 0.00 0.00
5 SWxy_5 77.00 80.00 16.00 19.00 16.40 0.00 0.00
6 SWxy_7 57.00 64.00 22.40 19.00 16.40 0.00 0.00
7 SFExy_3 41.00 41.00 9.60 15.00 15.00 0.00 0.00
8 SFExy_5 46.00 50.00 16.00 21.00 18.50 0.00 0.00
9 SFExy_7 43.00 46.00 22.40 21.00 18.50 0.00 0.00
10 SFExFy_3 43.00 0.00 9.60 17.00 12.50 0.00 0.00
11 SFExFy_5 41.00 0.00 16.00 20.00 15.00 0.00 0.00
12 SFExFy_7 38.00 0.00 22.40 20.00 15.00 0.00 0.00
13 SWxFy_3 77.00 0.00 9.60 15.00 10.00 0.00 0.00
14 SWxFy_5 68.00 0.00 16.00 20.00 15.00 0.00 0.00
15 SWxFy_7 51.00 0.00 22.40 20.00 15.00 0.00 0.00
16 AFxy_3 0.00 0.00 9.60 13.00 9.00 0.942 0.272
17 AFxy_5 0.00 0.00 16.00 17.50 10.00 2.545 0.395
18 AFxy_7 0.00 0.00 22.40 17.50 10.00 2.35 0.420
19 AFExy_3 52.00 46.00 9.60 13.50 9.00 4.12 2.14
20 AFExy_5 43.00 42.00 16.00 16.00 14.50 3.28 2.61
21 AFExy_7 37.00 36.00 22.40 16.00 14.50 2.98 2.35
22 AFExFy_3 47.00 0.00 9.60 13.50 9.00 0.71 2.11
23 AFExFy_5 38.00 0.00 16.00 16.00 14.50 0.45 2.61
24 AFExFy_7 35.00 0.00 22.40 16.00 14.50 0.45 2.45
25 AWxFy_3 64.00 0.00 9.60 14.50 9.00 0.30 3.51
26 AWxFy_5 69.00 0.00 16.00 14.00 16.00 2.80 1.11
27 AWxFy_7 65.00 0.00 22.40 14.00 16.00 2.76 1.20
28 AWxy_3 64.00 58.00 9.60 13.50 10.00 5.55 3.81
29 AWxy_5 65.00 72.00 16.00 16.25 16.25 3.11 5.46
30 AWxy_7 59.00 67.00 22.40 16.25 16.25 2.79 5.27
Notation: The names of the buildings code their structural characteristics as follows:. x, y= axes of buildings (referring to directions 1 and 2 respectively); S=Symmetric structural system along directions 1 (axis x) and 2 (axis y); A=Asymmetric structural system; W=Wall system according to EN1998-1 (the R/C shear walls receive more than 65% of base shear force); F=Frame system according to EN1998-1 (the R/C shear walls receive less than 35% of base shear force); FE=Dual system equivalent to Frame system according to EN1998-1 (the R/C shear walls receive 35%-50% of the base shear force); 3, 5, 7 are the number of storeys
Table 9. Structural data of the 3 examined buildings with unknown data used for the prediction of their DS.
Table 9. Structural data of the 3 examined buildings with unknown data used for the prediction of their DS.
Name nvx=nv1 nvy=nv2 Htot (m) Lx (m) Ly (m) ex (m) ey (m)
1 SFExFy_3 0.62 0.0 9.6 10.0 15.0 0.0 0.0
2 SFExFy_5 0.60 0.0 16.0 10.0 15.0 0.0 0.0
3 SFExFy_8 0.58 0.0 25.6 10.0 15.0 0.0 0.0
Table 10. Structural data of the three versions of the three case study buildings.
Table 10. Structural data of the three versions of the three case study buildings.
Building Year of construction Regularity in plan Existence of RC shear walls Regular distribution of masonries Strong masonry infills Pilotis
Dir 1 Dir 2
SFExFy_3B >2010 Yes Unknown Unknown Yes No No
SFExFy_5B >2010 Yes Yes Unknown No
SFExFy_8B >2010 Yes Unknown No No
SFExFy_3F >2010 Yes Unknown Unknown Yes Yes No
SFExFy_5F >2010 Yes Yes Unknown No
SFExFy_8F >2010 Yes Unknown No No
SFExFy_3P >2010 Yes Unknown Unknown Yes Yes Yes
SFExFy_5P >2010 Yes Yes Unknown No
SFExFy_8P >2010 Yes Unknown No No
Notation: The names of the three case study buildings are coded using the symbols used for the 30 known RC buildings in section 5.1. (Table 8). The characters B, F, P denote the “Bare Buildings” the “Fully Infilled Buildings” and the “Pilotis Buildings”.
Table 11. Assumed ranges of the values of the case study buildings’ unknown structural parameters.
Table 11. Assumed ranges of the values of the case study buildings’ unknown structural parameters.
Building nv ratios Eccentricity e0
Dir 1 (nv1) Dir 2 (nv2)
minnv1 maxnv2 minnv1 maxnv2 mine0 maxe0
SFExFy_3B 0.00 0.65 0.00 0.65 0.9014 1.352
SFExFy_5B 0.35 0.65 0.00 0.65 1.3521 2.704
SFExFy_8B 0.00 0.65 0.00 0.00 1.3521 2.704
SFExFy_3F 0.00 0.65 0.00 0.65 0.9014 1.352
SFExFy_5F 0.35 0.65 0.00 0.65 1.3521 2.704
SFExFy_8F 0.00 0.65 0.00 0.00 1.3521 2.704
SFExFy_3P 0.00 0.65 0.00 0.65 0.9014 1.352
SFExFy_5P 0.35 0.65 0.00 0.65 1.3521 2.704
SFExFy_8P 0.00 0.65 0.00 0.00 1.3521 2.704
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated