Post-Earthquake Damage Assessment—Case Study of the Educational Building after the Zagreb Earthquake

: In the wake of recent strong earthquakes in Croatia, there is a need for a detailed and more comprehensive post-earthquake damage assessment. Given that masonry structures are highly vulnerable to horizontal actions caused by earthquakes and a majority of the Croatian building stock is made of masonry, this ﬁeld is particularly important for Croatia. In this paper, a complete assessment of an educational building in Zagreb Lower Town is reported. An extensive program of visual inspection and geometrical surveys has been planned and performed. Additionally, an in situ shear strength test is presented. After extensive ﬁeldwork, collected data and results were input in 3Muri software for structural modeling. Moreover, a non-linear static (pushover) analysis was performed to individuate the possible failure mechanisms and to compare real-life damage to software results.


Introduction
On 22 March 2020, at 6 h 22 min, Zagreb Metropolitan area was hit by an earthquake of medium magnitude M L = 5.5, and intensity of VII, in the epicenter, according to the EMS-98 scale [1]. At 7 h 1 min followed the strongest subsequent earthquake of magnitude M L = 5.0 and intensity of VI. The main earthquake damaged most of the buildings in the Lower Town, including residential buildings, universities, schools, kindergartens, hospitals and public buildings. The vast majority of buildings built after the first mandatory earthquake regulations in former Yugoslavia (1964) [2,3] either remained intact or suffered small damage. Nonetheless, the larger part of the city's historical center (Upper and Lower Town) was severely damaged because the buildings in the center were built before any seismic regulations. The damage to historical buildings is enormous. Numerous museums, churches and university buildings have been severely damaged (Figure 1). At the end of the year, Croatia was hit by another devastating earthquake with an epicenter in Petrinja, located approx. 50 km from Zagreb (M L = 6.3). The quake caused subsequent damage to already damaged buildings, but to a lesser extent.
As well as most parts of the European region, many existing buildings in Croatia are built in masonry. Given that most of the so-called "strategic" buildings of cultural significance and high historical importance are built using masonry, such a condition is suggesting that the assessment and rehabilitation of existing masonry structures must be conducted on a very high level [4][5][6][7][8]. An important part of the structural assessment is numerical analysis. When it comes to existing buildings, a more refined non-linear analysis should be adopted. Non-linear static analysis or pushover analysis is important and is recommended in Eurocode 8-3 as a reference method for such situations. As well as most parts of the European region, many existing buildings in Croatia are built in masonry. Given that most of the so-called "strategic" buildings of cultural significance and high historical importance are built using masonry, such a condition is suggesting that the assessment and rehabilitation of existing masonry structures must be conducted on a very high level [4][5][6][7][8]. An important part of the structural assessment is numerical analysis. When it comes to existing buildings, a more refined non-linear analysis should be adopted. Non-linear static analysis or pushover analysis is important and is recommended in Eurocode 8-3 as a reference method for such situations.
After the earthquake, the first to respond to a disaster were civil engineers who led and coordinated the entire organization of building assessment and damage detection. Various similar post-earthquake assessment procedures are used worldwide [9][10][11].
In the first week, a large number of buildings were inspected, with a rapid post-earthquake assessment. The most endangered buildings in the city's center were the ones under cultural heritage protection. The aim of a rapid assessment of buildings is to determine the degree of damage to buildings concerning the protection of life and property, that is, to determine if the buildings are usable, temporarily unusable or unusable. Emerging technological advances allow the usage of artificial intelligence in the post-earthquake assessment process in the form of machine learning methods for more efficient and precise results [6,[12][13][14][15].
Zagreb's historic urban complex is a protected area regulated by the Law on the Protection and Preservation of Cultural Heritage. The area is divided into two zones, zone A and zone B (Figure 2) [1]. Zone A includes the oldest and most architecturally valuable After the earthquake, the first to respond to a disaster were civil engineers who led and coordinated the entire organization of building assessment and damage detection. Various similar post-earthquake assessment procedures are used worldwide [9][10][11].
In the first week, a large number of buildings were inspected, with a rapid postearthquake assessment. The most endangered buildings in the city's center were the ones under cultural heritage protection. The aim of a rapid assessment of buildings is to determine the degree of damage to buildings concerning the protection of life and property, that is, to determine if the buildings are usable, temporarily unusable or unusable. Emerging technological advances allow the usage of artificial intelligence in the postearthquake assessment process in the form of machine learning methods for more efficient and precise results [6,[12][13][14][15].
Zagreb's historic urban complex is a protected area regulated by the Law on the Protection and Preservation of Cultural Heritage. The area is divided into two zones, zone A and zone B ( Figure 2) [1]. Zone A includes the oldest and most architecturally valuable parts of Zagreb and is characterized by densely-built blocks of buildings made of stone, brick or a combination of materials. Most buildings consist of massive longitudinal and orthogonal walls and masonry ceiling vaults or wooden ceiling beams and wooden roofs ( Figure 3) [1]. Many hospitals, schools, business premises, residential and government buildings, cultural institutions, monuments, churches and chapels are located in zone A and are protected either as part of a historic urban complex or as individual heritage buildings per se. A total of 72% [1] of buildings in zone A suffered major damage due to the earthquake; to compare the damage suffered by this area is almost proportional to the Sustainability 2021, 13,6353 3 of 25 value of its cultural heritage. Zone B consists of a variety of urban patterns and a large number of immensely valuable buildings [16]. According to the World Bank report, in the educational sector, 106 buildings intended for preschool education, 214 primary and secondary school buildings and education centers, and 12 pupils' dorms were damaged. In the higher education subsector, the damage was reported to 152 buildings. In addition, the buildings of 29 research institutes were also affected. The total value of damage and losses to the education sector is estimated at EUR 1.8 billion at pre-disaster prices, with 97.9% affecting the City of Zagreb [1].
parts of Zagreb and is characterized by densely-built blocks of buildings made of stone, brick or a combination of materials. Most buildings consist of massive longitudinal and orthogonal walls and masonry ceiling vaults or wooden ceiling beams and wooden roofs ( Figure 3) [1]. Many hospitals, schools, business premises, residential and government buildings, cultural institutions, monuments, churches and chapels are located in zone A and are protected either as part of a historic urban complex or as individual heritage buildings per se. A total of 72% [1] of buildings in zone A suffered major damage due to the earthquake; to compare the damage suffered by this area is almost proportional to the value of its cultural heritage. Zone B consists of a variety of urban patterns and a large number of immensely valuable buildings [16]. According to the World Bank report, in the educational sector, 106 buildings intended for preschool education, 214 primary and secondary school buildings and education centers, and 12 pupils' dorms were damaged. In the higher education subsector, the damage was reported to 152 buildings. In addition, the buildings of 29 research institutes were also affected. The total value of damage and losses to the education sector is estimated at EUR 1.8 billion at pre-disaster prices, with 97.9% affecting the City of Zagreb [1].  In 2020, there were destructive earthquakes all around the world, causing loss of lives, building collapses, and severe structural and non-structural damage and economic losses (e.g., M7.0 Aegean Sea (Turkey-Greece) [17], M6.7 Elazig (East Turkey) [18], M5.5 parts of Zagreb and is characterized by densely-built blocks of buildings made of stone, brick or a combination of materials. Most buildings consist of massive longitudinal and orthogonal walls and masonry ceiling vaults or wooden ceiling beams and wooden roofs ( Figure 3) [1]. Many hospitals, schools, business premises, residential and government buildings, cultural institutions, monuments, churches and chapels are located in zone A and are protected either as part of a historic urban complex or as individual heritage buildings per se. A total of 72% [1] of buildings in zone A suffered major damage due to the earthquake; to compare the damage suffered by this area is almost proportional to the value of its cultural heritage. Zone B consists of a variety of urban patterns and a large number of immensely valuable buildings [16]. According to the World Bank report, in the educational sector, 106 buildings intended for preschool education, 214 primary and secondary school buildings and education centers, and 12 pupils' dorms were damaged. In the higher education subsector, the damage was reported to 152 buildings. In addition, the buildings of 29 research institutes were also affected. The total value of damage and losses to the education sector is estimated at EUR 1.8 billion at pre-disaster prices, with 97.9% affecting the City of Zagreb [1].  In 2020, there were destructive earthquakes all around the world, causing loss of lives, building collapses, and severe structural and non-structural damage and economic losses (e.g., M7.0 Aegean Sea (Turkey-Greece) [17], M6.7 Elazig (East Turkey) [18], M5.5 In 2020, there were destructive earthquakes all around the world, causing loss of lives, building collapses, and severe structural and non-structural damage and economic losses (e.g., M7.0 Aegean Sea (Turkey-Greece) [17], M6.7 Elazig (East Turkey) [18], M5.5 and M6.4 Croatia [1]). Identification of vulnerability characteristics and earthquake performance assessment of existing structures are essential steps in reducing earthquake losses, and the topic of seismic assessment of existing masonry structures is actual worldwide. Based on the available state-of-the-art literature on assessment and rehabilitation of existing masonry structures (e.g., [19][20][21][22][23]), this paper presents the Croatian perspective and shows it on an actual case study.

The Case Study
The subject of this paper is a building (Figures 4 and 5) located at Vlaška Street 38 as an attached building inside a block. Today's building was built in 1895 by adapting and upgrading two one-story buildings that were built in the early 19th century. The building was upgraded in 1906, while the building's current shape established complete reconstruction in 1997. The building was retrofitted in 1997 for educational purposes and seismic strengthening was not implemented. The building has a rectangular ground plan with the main orientation, which is the longer side of the building in the east-west direction. The building's external dimensions are 12 × 53 m, with two wings: one, 4.4 × 7.6 m and the other, 4.2 × 5.4 m located at the south side of the building. The total floor area of the building is approximately 685 m 2 . The building consists of a basement, first, second, third floor and attic. According to the Croatian seismic hazard map [24], the building is located in the area of peak ground acceleration intensity of a gR = 0.255 g for a return period of 475 years. The building serves as an educational-scientific institution. The condition of the building before the earthquake, regarding the vertical loads, was satisfactory, and the building was regularly maintained. and M6.4 Croatia [1]). Identification of vulnerability characteristics and earthquake performance assessment of existing structures are essential steps in reducing earthquake losses, and the topic of seismic assessment of existing masonry structures is actual worldwide. Based on the available state-of-the-art literature on assessment and rehabilitation of existing masonry structures (e.g., [19][20][21][22][23]), this paper presents the Croatian perspective and shows it on an actual case study.

The Case Study
The subject of this paper is a building (Figures 4 and 5) located at Vlaška Street 38 as an attached building inside a block. Today's building was built in 1895 by adapting and upgrading two one-story buildings that were built in the early 19th century. The building was upgraded in 1906, while the building's current shape established complete reconstruction in 1997. The building was retrofitted in 1997 for educational purposes and seismic strengthening was not implemented. The building has a rectangular ground plan with the main orientation, which is the longer side of the building in the east-west direction. The building's external dimensions are 12 × 53 m, with two wings: one, 4.4 × 7.6 m and the other, 4.2 × 5.4 m located at the south side of the building. The total floor area of the building is approximately 685 m 2 . The building consists of a basement, first, second, third floor and attic. According to the Croatian seismic hazard map [24], the building is located in the area of peak ground acceleration intensity of agR = 0.255 g for a return period of 475 years. The building serves as an educational-scientific institution. The condition of the building before the earthquake, regarding the vertical loads, was satisfactory, and the building was regularly maintained.    The original drawings show foundations that are approximately 1.0 m wide and 1.15 m deep in relation to the surrounding terrain. They were probably built in brick or stone, which is in line with the construction technology of the time. The building is built of a solid brick of a standard format 30 × 15 × 6.5 (unusual for today's standards) used in the late 19th century. The load-bearing wall thicknesses vary throughout the building, reducing with height, and are 51, 43, 28 cm (Figures 6 and 7). Plaster thickness also varies throughout the building from 3 to 6 cm. The ceiling structures before the reconstruction in 1997 were wooden beams, except for the basement ceiling and the first floor where the masonry vaults are located. After the 1997 reconstruction, the attic and 2nd floor ceilings (Figures 8 and 9) were replaced with reinforced concrete slabs, 12 cm and 16 cm, respectively. The 1st floor ceiling is a semi-precast masonry/concrete floor system (Fert ceiling) inserted between the existing wooden beams. In contrast, the ceilings on the ground floor and basement remained masonry vaults. The building also has two auxiliary staircases at the ends of the building made during the 1997 reconstruction and are made of reinforced concrete, and the main central staircase is older and is made of prefabricated stone stairs supported by a wall on one and a beam on the other edge. The original drawings show foundations that are approximately 1.0 m wide and 1.15 m deep in relation to the surrounding terrain. They were probably built in brick or stone, which is in line with the construction technology of the time. The building is built of a solid brick of a standard format 30 × 15 × 6.5 (unusual for today's standards) used in the late 19th century. The load-bearing wall thicknesses vary throughout the building, reducing with height, and are 51, 43, 28 cm (Figures 6 and 7). Plaster thickness also varies throughout the building from 3 to 6 cm. The ceiling structures before the reconstruction in 1997 were wooden beams, except for the basement ceiling and the first floor where the masonry vaults are located. After the 1997 reconstruction, the attic and 2nd floor ceilings (Figures 8 and 9) were replaced with reinforced concrete slabs, 12 cm and 16 cm, respectively. The 1st floor ceiling is a semi-precast masonry/concrete floor system (Fert ceiling) inserted between the existing wooden beams. In contrast, the ceilings on the ground floor and basement remained masonry vaults. The building also has two auxiliary staircases at the ends of the building made during the 1997 reconstruction and are made of reinforced concrete, and the main central staircase is older and is made of prefabricated stone stairs supported by a wall on one and a beam on the other edge.

Assessment Procedure
The first step in a complete post-earthquake building assessment is a rapid, preliminary assessment of the usability [25,26] of all buildings damaged in the earthquake. Additionally, it is of great importance to preserve the three-dimensional data of the facades of culturally-protected goods in the form of point clouds obtained by laser scanning or drone imaging. The mentioned data can also be used to assess existing structures for the creation of a 3D numerical model. Similar technology was used in the following papers [6,27,28]. In cases where it is needed, detailed assessment and available Non-Destructive Testing (NDT) assessment methods [29] are favorably used. A rapid preliminary assessment is conducted as early as possible after the earthquake, bearing in mind the safety of civil engineers in the field. In Croatia, this type of assessment consisted of a quick visual inspection of individual elements of the load-bearing structure, stating the appropriate degree of damage and deciding on the classification of the building into one of six possible categories ( Figure 10): U1 Usable without limitations (Green label), U2 Usable with recommendations Sustainability 2021, 13, 6353 7 of 25 (Green label), PN1 Temporary unusable-detailed inspection needed (Yellow label), PN2 Temporary unusable-emergency interventions needed (Yellow label), N1 Unusable due to external impacts (Red label) and N2 Unusable due to damage (Red label). [6,27,28]. In cases where it is needed, detailed assessment and available Non-Destructive Testing (NDT) assessment methods [29] are favorably used. A rapid preliminary assessment is conducted as early as possible after the earthquake, bearing in mind the safety of civil engineers in the field. In Croatia, this type of assessment consisted of a quick visual inspection of individual elements of the load-bearing structure, stating the appropriate degree of damage and deciding on the classification of the building into one of six possible categories ( Figure 10): U1 Usable without limitations (Green label), U2 Usable with recommendations (Green label), PN1 Temporary unusable-detailed inspection needed (Yellow label), PN2 Temporary unusable-emergency interventions needed (Yellow label), N1 Unusable due to external impacts (Red label) and N2 Unusable due to damage (Red label). Figure 10. Six categories of usability divided into three original labels (in Croatian) [30]. Figure 10. Six categories of usability divided into three original labels (in Croatian) [30].

Rapid Preliminary Assessment Results
A rapid assessment of the building in question was conducted on 23 March 2020. After a quick visual inspection of individual elements of the load-bearing structure, a decision was made to classify the building as temporarily unusable (Yellow label) with a recommendation for a detailed assessment (PN1). Basic conclusions of the preliminary assessment are:

Rapid Preliminary Assessment Results
A rapid assessment of the building in question was conducted on 23 March 2020. After a quick visual inspection of individual elements of the load-bearing structure, a decision was made to classify the building as temporarily unusable (Yellow label) with a recommendation for a detailed assessment (PN1). Basic conclusions of the preliminary assessment are:

•
There is visible damage in the form of cracks on the wall coverings, arches ( Figure  11a), vaults and ceilings ( Figure 11b) on all floors; • Separation and local decay of plaster; • Minor local damage is visible on structural elements (walls, columns, arches);

•
In the eastern part of the building, diagonal cracks are visible on the load-bearing walls.
The second floor and attic suffered minor damage, while the most severe damage is found on the eastern (Figure 12a,b) and central staircase wings (Figure 13a,b). Recommendations were given that the building can be used with a restriction in the zones where there is a danger of plaster falling. Additionally, the eastern and western staircase can be used with a restriction in the number of people, while the central staircase is not to be used until a detailed assessment is conducted.  The second floor and attic suffered minor damage, while the most severe damage is found on the eastern (Figure 12a,b) and central staircase wings (Figure 13a,b). Recommendations were given that the building can be used with a restriction in the zones where there is a danger of plaster falling. Additionally, the eastern and western staircase can be used with a restriction in the number of people, while the central staircase is not to be used until a detailed assessment is conducted.

Detailed Assessment Results
According to the current standards for the design of structures-a series of Eurocodes, HRN EN 1990-1998 and the relevant national annexes, the building that is the subject of this study is in the range of peak ground acceleration of 0.255 g; that is, the expected earthquake intensity is IX according to EMS-98 scale for a return period of 475 years. No geotechnical tests have been performed for the site in question for this article, but based on empirical data, a category B foundation soil (deposits of very compacted sand, gravel or hard clay, at least several tens of meters deep) or category C (deep deposits of compacted or medium-compacted sand, gravel or hard clay with a thickness of several tens of meters to several hundred meters) can be assumed. Moreover, based on the latest findings obtained from the research of the Croatian Geological Institute in cooperation with

Detailed Assessment Results
According to the current standards for the design of structures-a series of Eurocodes, HRN EN 1990-1998 and the relevant national annexes, the building that is the subject of this study is in the range of peak ground acceleration of 0.255 g; that is, the expected earthquake intensity is IX according to EMS-98 scale for a return period of 475 years. No geotechnical tests have been performed for the site in question for this article, but based on empirical data, a category B foundation soil (deposits of very compacted sand, gravel or hard clay, at least several tens of meters deep) or category C (deep deposits of compacted or medium-compacted sand, gravel or hard clay with a thickness of several tens of meters to several hundred meters) can be assumed. Moreover, based on the latest findings obtained from the research of the Croatian Geological Institute in cooperation with

Detailed Assessment Results
According to the current standards for the design of structures-a series of Eurocodes, HRN EN 1990-1998 and the relevant national annexes, the building that is the subject of this study is in the range of peak ground acceleration of 0.255 g; that is, the expected earthquake intensity is IX according to EMS-98 scale for a return period of 475 years. No geotechnical tests have been performed for the site in question for this article, but based on empirical data, a category B foundation soil (deposits of very compacted sand, gravel or hard clay, at least several tens of meters deep) or category C (deep deposits of compacted or medium-compacted sand, gravel or hard clay with a thickness of several tens of meters to several hundred meters) can be assumed. Moreover, based on the latest findings obtained from the research of the Croatian Geological Institute in cooperation with the University of Zagreb, a seismic microzonation map was prepared according to Eurocode 8 standards for the Zagreb area [31]. According to the mentioned seismic microzonation (2017-2019), the soil in the immediate vicinity of the assessed building belongs to the category of soil type C. Soil type C causes a certain amplification of the soil shaking, which must be taken into account when assessing the condition of the structure.
All damage, structural and non-structural, is photographed and described. They are plotted in the floor plans of the building (Figures 14-16). The building was inspected from the air by an unmanned aerial vehicle, and no damage was observed to the main load-bearing structure or the building's roof structure. Decorative crosses, statues and reliefs were also inspected. For the purposes of digital preservation, the 3D model of the building was made on the basis of photogrammetric images. repaired and reinforced. Before that, it is necessary to do all the research work to determine the characteristics of the masonry and other necessary data for the structural analysis.

In Situ Masonry Shear Strength Tests
In order to assess the condition of the structure after the earthquake and corresponding analysis of the existing condition of the load-bearing structure, in situ tests were carried out. Determination of shear strength (mortar in the composition of load-bearing ma- repaired and reinforced. Before that, it is necessary to do all the research work to determine the characteristics of the masonry and other necessary data for the structural analysis.

In Situ Masonry Shear Strength Tests
In order to assess the condition of the structure after the earthquake and corresponding analysis of the existing condition of the load-bearing structure, in situ tests were carried out. Determination of shear strength (mortar in the composition of load-bearing ma- repaired and reinforced. Before that, it is necessary to do all the research work to determine the characteristics of the masonry and other necessary data for the structural analysis.

In Situ Masonry Shear Strength Tests
In order to assess the condition of the structure after the earthquake and corresponding analysis of the existing condition of the load-bearing structure, in situ tests were carried out. Determination of shear strength (mortar in the composition of load-bearing ma- A detailed inspection of the building revealed the following damage: at the ground floor, the damage is visible in the form of cracks on the wall coverings, arches, vaults and ceilings, as well as separation and local decay of plaster. Cracks in the barrel vaults are mostly parallel to the supporting joints, probably due to lateral movements during the earthquake. They are the result of the occurrence of tensile stresses perpendicular to the supporting joint. Such cracks can cause hinge formation and consequent loss of stability if they propagate deep enough. Fortunately, the cracks in the assessed building are narrow and mostly found in the plaster. Minor local damage to structural elements (walls, columns, arches) is also visible. In the central core of the building where the main staircase is located and in the eastern part of the building, diagonal cracks are visible on the load-bearing walls, which can also be seen on the north side of the building. Damage is visible on all floors in the form of cracks and falling plaster on the walls. Minor local damage to the walls on the 1st floor is also visible, and cracks at the joints of partition walls and ceilings are locally visible. In the central wing of the building where the main staircase is located and in the eastern wing of the building, diagonal cracks are visible on the load-bearing walls, which can also be seen on the north side of the building.
The 2nd floor and attic suffered minor damage. Particular attention should be paid to the central part of the building, occupying the wing with the staircase. The formation of cracks on the central wing transverse walls is clearly visible, indicating a possible failure mechanism out-of-plane of the entire wing. A wedge was made, and the cracks were interconnected and propagated inside the building (they also appear in the stair beams). There was no displacement of the walls out-of-plane, but the preconditions for its failure were met. The central wing needs to be strengthened as soon as possible as part of the entire building's renovation.
The east wing was also damaged at the ground floor and 1st floor level. The cracks that appeared propagated were through the entire wall of the south facade of the wing. It is unfavorable that the cracks are joined and continue to the transversely-joined walls and lintels. The cause of such cracks can be the slight contribution of a torsional response of the building as a whole, where the boundary elements are the most loaded ones, and their failure occurs. Additionally, that part with the building is connected to the neighboring building. Although this can generally have a positive effect on the whole building, in the case of walls on the east wing, such a boundary condition can cause additional forces. If the walls are not well connected to the diaphragms by a tensile compression connection, this can cause them to fail. Since there has been no displacement of the wall out-of-plane, it is not in danger of collapsing, but it should be strengthened soon, and further damage propagation should be prevented. Minor damage in the form of cracks on the walls can be seen on the west and east staircases.
The building can be used in its entirety except for the main staircase. Depending on the possibilities in the future, a static and seismic analysis of the existing condition of the building should be made, and regardless of whether the entire building will be reinforced, the main staircase and other walls with cracks along the entire width of the walls must be repaired and reinforced. Before that, it is necessary to do all the research work to determine the characteristics of the masonry and other necessary data for the structural analysis.

In Situ Masonry Shear Strength Tests
In order to assess the condition of the structure after the earthquake and corresponding analysis of the existing condition of the load-bearing structure, in situ tests were carried out. Determination of shear strength (mortar in the composition of load-bearing masonry) of solid brick masonry [32,33] was performed "in situ" using a small hydraulic press "Holmatro" with a load capacity of 200 kN. The mortar was moved horizontally in the vicinity of one brick in order to determine the shear strength. At the same time, the structure of the existing wall was minimally damaged. A total of eleven positions on the load-bearing walls were selected for testing the shear strength of the masonry (Figures 14-16 The shear strength test of the masonry was carried out in nine places. After removing the plaster, it was found that due to the method of masonry (no bricks were found in the longitudinal direction-the inner part of the wall is built of bricks "on edge"), conditions to perform the test were not met for positions PS-PR-3 and PS-PR-4.

Results of the Shear Strength Tests
The shear strength of the masonry was obtained based on the registered horizontal force H umax acting on one brick at the time of reaching the shear strength in that brick and the corresponding mortar area on both sides of the shear is transmitted (Ag + Ad). The test method can be seen in Figure 17 and photographs (Figures 18 and 19). The shear strength test of the masonry was carried out in nine places. After removing the plaster, it was found that due to the method of masonry (no bricks were found in the longitudinal direction-the inner part of the wall is built of bricks "on edge"), conditions to perform the test were not met for positions PS-PR-3 and PS-PR-4.

Results of the Shear Strength Tests
The shear strength of the masonry was obtained based on the registered horizontal force Humax acting on one brick at the time of reaching the shear strength in that brick and the corresponding mortar area on both sides of the shear is transmitted (Ag + Ad). The test method can be seen in Figure 17 and photographs (Figures 18 and 19).  The shear strength test of the masonry was carried out in nine places. After removing the plaster, it was found that due to the method of masonry (no bricks were found in the longitudinal direction-the inner part of the wall is built of bricks "on edge"), conditions to perform the test were not met for positions PS-PR-3 and PS-PR-4.

Results of the Shear Strength Tests
The shear strength of the masonry was obtained based on the registered horizontal force Humax acting on one brick at the time of reaching the shear strength in that brick and the corresponding mortar area on both sides of the shear is transmitted (Ag + Ad). The test method can be seen in Figure 17 and photographs (Figures 18 and 19).          It is the shear strength f v with the contribution of σ 0 -vertical stress. During the test, each measuring position is precisely located to calculate the vertical load (G 0 ), that is, vertical stress (σ 0 ) from the numerical model. Therefore, for each test site, in addition to the floor plan position, the height of the measuring position from the upper edge of the ceiling structure (h) is recorded. When analyzing the shear strength of mortar, the vertical constant load is taken into account, that is, vertical stresses at a particular test position. Shear strength according to Mohr-Coulomb law is calculated by Equation (1).
It can be seen from the diagram that the shear strength without vertical pressure is 0.316 MPa and that coefficient of friction is 1.303. Due to the high value in comparison to EC standard recommendation, the friction coefficient was taken as µ = 0.40 (according to [34]). The obtained shear strength of masonry without vertical pressure, i.e., cohesion, is higher than the one provided by the regulations for the case when there are no tests (f vm0 = 0.32 MPa > 0.10 MPa). The results show that the quality of masonry is good in contrast to similar buildings from that period.

Numerical Modeling
The 3D numerical model of the assessed building is obtained in 3Muri software. The macro-element approach is adopted due to computational efficiency and high precision [35]. Its versatility in modeling (implementing elements of various materials, realistic floor stiffnesses, strengthening and many more) makes it highly valuable in a region where a vast majority of building stock is made of masonry. Similar case studies in 3Muri software were used as a base for our research [36][37][38].
Macro-element approach implies equivalent-frame method which uses non-linear beam elements. Macro-elements (or non-linear beam elements) are divided into three categories which are piers, spandrels and rigid nodes. In piers and spandrels, all deformation is concentrated, and they are connected with rigid nodes. Figure 21 shows an equivalent frame model made of mentioned macro-elements. It is the shear strength fv with the contribution of -vertical stress. During the test, each measuring position is precisely located to calculate the vertical load (G0), that is, vertical stress ( ) from the numerical model. Therefore, for each test site, in addition to the floor plan position, the height of the measuring position from the upper edge of the ceiling structure (h) is recorded. When analyzing the shear strength of mortar, the vertical constant load is taken into account, that is, vertical stresses at a particular test position. Shear strength according to Mohr-Coulomb law is calculated by Equation (1).
It can be seen from the diagram that the shear strength without vertical pressure is 0.316 MPa and that coefficient of friction is 1.303. Due to the high value in comparison to EC standard recommendation, the friction coefficient was taken as = 0.40 (according to [34]). The obtained shear strength of masonry without vertical pressure, i.e., cohesion, is higher than the one provided by the regulations for the case when there are no tests (fvm0 = 0.32 MPa > 0.10 MPa). The results show that the quality of masonry is good in contrast to similar buildings from that period.

Numerical Modeling
The 3D numerical model of the assessed building is obtained in 3Muri software. The macro-element approach is adopted due to computational efficiency and high precision [35]. Its versatility in modeling (implementing elements of various materials, realistic floor stiffnesses, strengthening and many more) makes it highly valuable in a region where a vast majority of building stock is made of masonry. Similar case studies in 3Muri software were used as a base for our research [36][37][38].
Macro-element approach implies equivalent-frame method which uses non-linear beam elements. Macro-elements (or non-linear beam elements) are divided into three categories which are piers, spandrels and rigid nodes. In piers and spandrels, all deformation is concentrated, and they are connected with rigid nodes. Figure 21 shows an equivalent frame model made of mentioned macro-elements. Non-linear static pushover analysis [39,40] allows us to check the overstrength ratio used in linear analysis and it gives us more detailed insight into critical elements, possible failure mechanisms, and global behavior of the building as a whole. Pushover analysis is performed with constant gravity loads and monotonically-increasing horizontal loads. Two different distributions of the horizontal loads along the structure's height are used Non-linear static pushover analysis [39,40] allows us to check the overstrength ratio used in linear analysis and it gives us more detailed insight into critical elements, possible failure mechanisms, and global behavior of the building as a whole. Pushover analysis is performed with constant gravity loads and monotonically-increasing horizontal loads. Two different distributions of the horizontal loads along the structure's height are used for the pushover analysis. The first distribution has a uniform pattern where the horizontal load is proportional to the mass of the building. The second distribution has a modal pattern where the horizontal load is distributed along with the building's height proportionally to the first vibration mode shape of the building determined through elastic analysis (Figures 22 and 23). These horizontal loads are applied at the location of the masses in the model, i.e., at each floor level in the center of masses. Moreover, accidental eccentricity is taken into account to cover uncertainty in the calculation of the center of masses of the building. The 5% of the building's length perpendicular to seismic load direction is taken into account on each side for both x (longitudinal) and y (transverse) directions.
for the pushover analysis. The first distribution has a uniform pattern where the horizontal load is proportional to the mass of the building. The second distribution has a modal pattern where the horizontal load is distributed along with the building's height proportionally to the first vibration mode shape of the building determined through elastic analysis (Figures 22 and 23). These horizontal loads are applied at the location of the masses in the model, i.e., at each floor level in the center of masses. Moreover, accidental eccentricity is taken into account to cover uncertainty in the calculation of the center of masses of the building. The 5% of the building's length perpendicular to seismic load direction is taken into account on each side for both x (longitudinal) and y (transverse) directions.  In each incremental step, internal forces are redistributed according to the element equilibrium, and stiffness is degraded in plastic range. Additionally, ductility is controlled with maximum drift which is 0.004 for shear failure and 0.008 for bending failure [34,41]. Generally, masonry walls have three main failure modes as described in [41][42][43][44]. Turnsek-Cacovic constitutive law is used as diagonal cracking is usually the dominant failure mode for existing unreinforced masonry structures [45,46]. In the work of [47], diagonal failure strength is correlated with shear strength used for Turnsek-Cacovic constitutive law. Shear strength obtained by in situ testing is compared with shear strength approximation from the visual MQI method which is explained in more detail in [48] and further developed in [49,50]. Compared shear strengths are very close for this case study which implies good precision of the MQI method. Cracked stiffness of vertical elements is used in a model as recommended in [34], so that cracking that occurs during lifetime because of expected earthquakes of a smaller magnitude is taken into account. Shear and flexural stiffness are taken as half value of initial stiffness. In a 3D model, floors are modeled as horizontally rigid diaphragms, which is precise enough due to the real in-plane stiffness of the horizontal floor structures. Axial in-plane stiffness of rigid diaphragms in software is infinite and the mass of the real slab is taken into account. Many similar old masonry buildings have an unfavorable distribution of seismic forces due to traditional flexible timber floors [22]. In seismic analysis, the roof is In each incremental step, internal forces are redistributed according to the element equilibrium, and stiffness is degraded in plastic range. Additionally, ductility is controlled with maximum drift which is 0.004 for shear failure and 0.008 for bending failure [34,41]. Generally, masonry walls have three main failure modes as described in [41][42][43][44]. Turnsek-Cacovic constitutive law is used as diagonal cracking is usually the dominant failure mode for existing unreinforced masonry structures [45,46]. In the work of [47], diagonal failure strength is correlated with shear strength used for Turnsek-Cacovic constitutive law. Shear strength obtained by in situ testing is compared with shear strength approximation from the visual MQI method which is explained in more detail in [48] and further developed in [49,50]. Compared shear strengths are very close for this case study which implies good precision of the MQI method. Cracked stiffness of vertical elements is used in a model as recommended in [34], so that cracking that occurs during lifetime because of expected earthquakes of a smaller magnitude is taken into account. Shear and flexural stiffness are taken as half value of initial stiffness. In a 3D model, floors are modeled as horizontally rigid diaphragms, which is precise enough due to the real in-plane stiffness of the horizontal floor structures. Axial in-plane stiffness of rigid diaphragms in software is infinite and the mass of the real slab is taken into account. Many similar old masonry buildings have an unfavorable distribution of seismic forces due to traditional flexible timber floors [22]. In seismic analysis, the roof is excluded from the load-bearing structure because it does not significantly affect the response of the structure and does not contribute to the global resistance of the structure. Although it was left out of the structural part, its contribution in the form of load on the structure itself was not neglected.
The mean values of material characteristics used in the numerical model (Table 2) are a combination of the literature review [34,51] and on-site testing. Regarding experimental in situ tests and detailed inspection of the structure knowledge level 2 (normal knowledge) can be defined. Based on the achieved knowledge level, confidence factor was taken as 1.2.  In each incremental step, internal forces are redistributed according to the element equilibrium, and stiffness is degraded in plastic range. Additionally, ductility is controlled with maximum drift which is 0.004 for shear failure and 0.008 for bending failure [34,41]. Generally, masonry walls have three main failure modes as described in [41][42][43][44]. Turnsek-Cacovic constitutive law is used as diagonal cracking is usually the dominant failure mode for existing unreinforced masonry structures [45,46]. In the work of [47], diagonal failure strength is correlated with shear strength used for Turnsek-Cacovic constitutive law. Shear strength obtained by in situ testing is compared with shear strength approximation from the visual MQI method which is explained in more detail in [48] and further developed in [49,50]. Compared shear strengths are very close for this case study which implies good precision of the MQI method. Cracked stiffness of vertical elements is used in a model as recommended in [34], so that cracking that occurs during lifetime because of expected earthquakes of a smaller magnitude is taken into account. Shear and flexural stiffness are taken as half value of initial stiffness. In a 3D model, floors are modeled as horizontally rigid diaphragms, which is precise enough due to the real in-plane stiffness of the horizontal floor structures. Axial in-plane stiffness of rigid diaphragms in software is infinite and the mass of the real slab is taken into account. Many similar old masonry buildings have an unfavorable distribution of seismic forces due to traditional flexible timber floors [22]. In seismic analysis, the roof is excluded from the load-bearing structure because it does not significantly affect the response of the structure and does not contribute to the global resistance of the structure.

Material Characteristic Value
Although it was left out of the structural part, its contribution in the form of load on the structure itself was not neglected.
The mean values of material characteristics used in the numerical model (Table 2) are a combination of the literature review [34,51] and on-site testing. Regarding experimental in situ tests and detailed inspection of the structure knowledge level 2 (normal knowledge) can be defined. Based on the achieved knowledge level, confidence factor was taken as 1.2. According to the work in [52], the building is classified as regular in height but irregular in floor plan, requiring 3D modeling. The building is classified as a torsional stiff system.
First, static analysis is performed according to [53]. Next, the seismic analysis is done. The educational building belongs to importance class III because its seismic resistance is of great importance given the consequences associated with a collapse. Hence, importance factor is γ I = 1.2. Three PGA values are used for two limit states. According to the new law "Law on the Reconstruction of Earthquake-Damaged Buildings in the City of Zagreb, Krapina-Zagorje County and Zagreb County (NN 102/2020)" [54], ultimate limit state return period can be different depending on the level of strengthening for old masonry buildings damaged in the recent earthquakes. Limit state of significant damage with a return period of 475 and limit state of damage limitation with a return period of 95 years were checked [55]. In the new law [54], the return period of 225 years which corresponds to a probability of exceedance of 20% in 50 years is introduced for a limit state of significant damage.
Elastic response spectrums for acceleration are calculated for all three return periods taking into account parameters for soil type C, which is found on the location of the building. Altogether, 24 pushover analyses are performed; for x-and y-direction in both orientation, with two load distributions, without and with −/+ 5% of accidental eccentricity.
The result of the performed seismic analysis is a capacity curve that shows the ratio of the shear force in the base of the structure and the displacement of the control node. The control node was selected in the immediate vicinity of the center of mass and is located on the top floor of the building. Obtained capacity curves for all 24 analyses can be seen in Figure 24. Bilinearized pushover curves for the x-and y-direction are shown in Figures 25 and 26. Total base shear in kN is plotted on the y-axis and the displacement of the control nodes in mm is plotted on the x-axis.             After the response of the structure, the capacity of the structure is obtained and checks are carried out according to the basic requirements relating to the state of structural damage, defined by limit states. Parameters for equivalent SDOF systems from Figures 27 and 28 are shown in Table 3. These parameters are obtained during bilinearization based on the equivalent energy principal and are used for target displacement determination. Results are also given in the form of a parameter α (Table 4), where α is the ratio between the limit capacity acceleration of the building and reference peak ground acceleration on type A ground. Parameter α is given for all limit states. A problem arises with After the response of the structure, the capacity of the structure is obtained and checks are carried out according to the basic requirements relating to the state of structural damage, defined by limit states. Parameters for equivalent SDOF systems from Figures 27 and 28 are shown in Table 3. These parameters are obtained during bilinearization based on the equivalent energy principal and are used for target displacement determination. Results are also given in the form of a parameter α (Table 4), where α is the ratio between the limit capacity acceleration of the building and reference peak ground acceleration on type A ground. Parameter α is given for all limit states. A problem arises with old masonry buildings that often cannot be strengthened to such an extent that today's building codes regarding seismic resistance are satisfied. Hence, a new legal document [54] that followed recent earthquakes allows for different levels of earthquake resistance after reconstruction. For the return period of 475 years, none of the 24 analyses satisfies the limit state of significant damage, while for the return period of 225 years, only 8 of the 24 analyses satisfies the same limit state. For a return period of 95 years, none of the 24 analyses satisfies the limit state of limited damage.
The target displacement is determined in accordance with Eurocode 8 (Appendix B) [52]. Since the building's natural period of vibration is less than the period of T C , the target displacement is determined by the procedure for short periods shown in Figure 29. Buildings with short natural periods of vibration do not comply with the equal displacement rule [56] as buildings with medium-long and long natural periods of vibration.

95
0.560 0.363 For the return period of 475 years, none of the 24 analyses satisfies the limit state of significant damage, while for the return period of 225 years, only 8 of the 24 analyses satisfies the same limit state. For a return period of 95 years, none of the 24 analyses satisfies the limit state of limited damage.
The target displacement is determined in accordance with Eurocode 8 (Appendix B). [52], Since the building's natural period of vibration is less than the period of TC, the target displacement is determined by the procedure for short periods shown in Figure 29. Buildings with short natural periods of vibration do not comply with the equal displacement rule [56] as buildings with medium-long and long natural periods of vibration. Since 3Muri does not consider the bending failure of walls out-of-plane during the global non-linear static analysis, it is necessary to subsequently perform an analysis of the bending of the walls out-of-plane. This takes into account the seismic action perpendicular to the walls. The analysis is performed for the boundary condition near collapse. A return period of 475 years is used. Walls are considered to have passed the out-of-plane bending check if their MRd/MEd ratio is greater than or equal to 1.0. In Figure 30, walls that did not satisfy the check are colored red. Since 3Muri does not consider the bending failure of walls out-of-plane during the global non-linear static analysis, it is necessary to subsequently perform an analysis of the bending of the walls out-of-plane. This takes into account the seismic action perpendicular to the walls. The analysis is performed for the boundary condition near collapse. A return period of 475 years is used. Walls are considered to have passed the out-of-plane bending check if their M Rd /M Ed ratio is greater than or equal to 1.0. In Figure 30, walls that did not satisfy the check are colored red. The 3Muri program does not consider out-of-plane loss of stability of local mechanisms in the global analysis either. It assumes that proper connection is established between walls and between walls and diaphragms. That way, out-of-plane local mechanisms are prevented so that the global in-plane response of the building can be evaluated [43]. Therefore, the resistance to local losses of stability is checked in a special program module; certain parts of a single wall are checked and the interaction of parts of several individual walls that can together form different local mechanisms. Figure 31 shows an example of a local mechanism where the kinematic block is colored red. The 3Muri program does not consider out-of-plane loss of stability of local mechanisms in the global analysis either. It assumes that proper connection is established between walls and between walls and diaphragms. That way, out-of-plane local mechanisms are prevented so that the global in-plane response of the building can be evaluated [43]. Therefore, the resistance to local losses of stability is checked in a special program module; certain parts of a single wall are checked and the interaction of parts of several individual walls that can together form different local mechanisms. Figure 31 shows an example of a local mechanism where the kinematic block is colored red.
Local mechanisms are arbitrarily defined according to the structure's shape, common failure mechanisms and earthquake damage. The emergence of local mechanisms is often due to the poor interconnection of walls and walls with floor structures. Linear kinematic analysis is used. Defining a local mechanism consists of three steps. To begin with, it is necessary to define a kinematic block that is a part of a wall that is considered absolutely rigid, and that is subject to movement or tilting relative to another block or the rest of the wall. Then, the boundary conditions are defined and finally, the load needs to be set. Some of the possible local failure mechanisms of the observed structure and its resistance to the shown forms of failure are presented below ( Figure 32 and Table 5). Parameter α represents the ratio between the spectral acceleration of the activation of the mechanisms and spectral acceleration of the seismic demand. The 3Muri program does not consider out-of-plane loss of stability of local mechanisms in the global analysis either. It assumes that proper connection is established between walls and between walls and diaphragms. That way, out-of-plane local mechanisms are prevented so that the global in-plane response of the building can be evaluated [43]. Therefore, the resistance to local losses of stability is checked in a special program module; certain parts of a single wall are checked and the interaction of parts of several individual walls that can together form different local mechanisms. Figure 31 shows an example of a local mechanism where the kinematic block is colored red. Local mechanisms are arbitrarily defined according to the structure's shape, common failure mechanisms and earthquake damage. The emergence of local mechanisms is often due to the poor interconnection of walls and walls with floor structures. Linear kinematic analysis is used. Defining a local mechanism consists of three steps. To begin with, it is necessary to define a kinematic block that is a part of a wall that is considered absolutely rigid, and that is subject to movement or tilting relative to another block or the rest of the wall. Then, the boundary conditions are defined and finally, the load needs to be set. Some of the possible local failure mechanisms of the observed structure and its resistance to the shown forms of failure are presented below ( Figure 32 and Table 5). Parameter α represents the ratio between the spectral acceleration of the activation of the mechanisms and spectral acceleration of the seismic demand.  After the performed analyses, the real damage was compared with the damage distribution previewed in 3Muri (Figures 33-36). Description of building's capacity through displacement rather than forces allows us to better understand and accurately predict a building's response in the form of a damage initiation and propagation throughout all phases of analysis all the way until the formation of failure mechanism and collapse. That is possible, thanks to incremental non-linear static pushover analysis that shows us damage patterns for every macro-element in every step of the analysis. Compared results are very similar to real damage, which implies good accuracy of the used analysis and software.  After the performed analyses, the real damage was compared with the damage distribution previewed in 3Muri (Figures 33-36). Description of building's capacity through displacement rather than forces allows us to better understand and accurately predict a building's response in the form of a damage initiation and propagation throughout all phases of analysis all the way until the formation of failure mechanism and collapse. That is possible, thanks to incremental non-linear static pushover analysis that shows us damage patterns for every macro-element in every step of the analysis. Compared results are very similar to real damage, which implies good accuracy of the used analysis and software. displacement rather than forces allows us to better understand and accurately predict a building's response in the form of a damage initiation and propagation throughout all phases of analysis all the way until the formation of failure mechanism and collapse. That is possible, thanks to incremental non-linear static pushover analysis that shows us damage patterns for every macro-element in every step of the analysis. Compared results are very similar to real damage, which implies good accuracy of the used analysis and software.      Another important piece of information is that the studied building is a part of a masonry building aggregate, which is not rare in Zagreb and other European cities. The behavior of buildings in the aggregate is extremely complex, where many parameters affect their response during an earthquake. Moreover, buildings at the ends of aggregates are often more severely damaged than those inside aggregates during earthquakes [57,58]. Modeling of the entire aggregate or, at least, adjacent buildings is recommended. Adjacent buildings can be modeled as mutually independent, mutually interconnected by the same walls, or mutually interconnected by various connections that accurately simulate real coherence. Those connections transfer axial and/or shear load between adjacent buildings in aggregate and they can be modeled as linear or non-linear elements [58][59][60]. Due to the lack of data on neighboring buildings, modeling of only the observed building is often resorted to. This is on the safe side in the case of a unit inside the aggregate when deformations are observed, but there is a risk of misinterpretation of the failure mechanism [61].
The observed building was modeled as an isolated unit due to insufficient data on neighboring buildings. In order to roughly assess the effect of neighboring buildings on the observed building, the vulnerability index proposed in the articles [62,63] is used. The index is based on five categories depending on the height of the surrounding buildings, the position of the building in the unit, the mismatch of floor heights, differences in material characteristics between buildings and differences in the areas of openings on the facade walls. This method was developed as an upgrade to the vulnerability index based on the ten parameters described in [64]. According to the mentioned five parameters, the buildings in the aggregate have a minimal impact on the observed building, and it can be assumed that the modeling of the building is as isolated as on the safe side.

Discussion and Conclusions
Based on similar case studies from the recent earthquakes in Italy [37,38,[65][66][67][68][69], the assessment of the existing educational building is made in order to determine its seismic resistance. The building was damaged in the recent earthquakes in Croatia. The building was not constructed according to the principles of seismic design, but the reconstruction in 1997 partially improved the condition of the structure by adding transverse walls and replacing old wooden beams with reinforced concrete floors. Such rigid diaphragms enable good connection of all walls and thus, better behavior of the building in an earthquake. That is why increasing the stiffness of the traditional timber floors in old masonry buildings is usually one of the first measures in seismic retrofitting. On the other hand, new research suggests that replacing traditional wooden floors with rigid diaphragms, i.e., RC floors, can induce some unwanted consequences such as cracks on the edges of the two materials or, in the worst scenario, disintegration and collapse of the masonry walls. However, for earthquakes with expected smaller magnitudes as the ones in Zagreb, this measure is convenient and serves to reinforce the existing structure to horizontal actions.
The assessment of the existing unconfined masonry structure was performed using 3Muri with an equivalent frame method. Rather than linear analysis, non-linear static (pushover) seismic analysis is used due to its various benefits. Limit state checks were carried out for the return periods of 95, 225 and 475 years. The results are in line with the expected behavior with respect to the existing wall distribution in the structure. The structure is less rigid and has greater displacements in the y-direction. Additionally, the capacity of the structure is smaller for the y-direction. There is also a small eccentricity between the center of rigidity and the center of mass that causes a slight but negative impact of torsion on the global behavior of the structure. The critical elements are the walls of the central stairwell and the transverse walls on the west side of the building. Additionally, bending failure out-of-plane was checked. Local failure mechanisms were also analyzed using linear kinematic analysis. A comparison of the actual damage with the damage obtained as a result of the conducted non-linear static seismic analysis in the 3Muri program was also conducted.
After assessing the building's behavior in an earthquake, it is obvious that strengthening is needed to raise the seismic resistance level of the building. Damage caused by an earthquake must be repaired to prevent the progression of damage and the possible threat to the global resistance and stability of the building.
As Croatia was hit by two severe earthquakes last year, the knowledge of 'build back better' is fully appreciated. That means that sustainable materials and innovative concepts [70][71][72][73] should be used and energy efficiency ensured [74,75]. Various methods are used for strengthening of masonry buildings [76][77][78] and materials such as FRP and TRM will probably be adopted due to their compatibility and reversibility.
This study provides a detailed insight into the behavior of the building during the earthquake, as well as the extent and distribution of damage and critical elements for earthquakes of different intensities.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy reasons.

Conflicts of Interest:
The authors declare no conflict of interest.