Seismic and Fire Assessment and Upgrading Process for Historical Buildings: The Case Study of Palazzo Colonna in Caggiano

The assessment and retrofit of existing masonry structures with historical and cultural value in highly seismic zones are challenging issues in earthquake engineering. In fact, the historic and recent earthquakes have shown the problem of the seismic vulnerability of existing masonry constructions. A historical masonry palace located in Caggiano (Salerno, Italy) is used herein as a case study, showing the vulnerability assessment and the seismic upgrading process. The case study building has a masonry structural type at the first two floors while there is a third floor realized in reinforced concrete and a fourth floor realized with a wood structure. The building was characterized by a remarkable seismic vulnerability and needed seismic upgrading operations. After the vulnerability assessment process, some design suggestions are proposed for the seismic upgrading of the building. The structure before and after the upgrading operations has been checked through nonlinear static and dynamic analyses. Then, coherently with the “Sismabonus” approach, the attribution of the seismic risk class, performed through numerical analyses, is founded on two parameters, namely, the expected annual mean losses (PAM), related to economic factors, and the Life Safety Index (IS-V), related to the structure seismic safety. Finally, the overcoming of the different classes of risk is shown and compared with the amount of the retrofit operations, their costs, and the impact on the existing space. Moreover, fire assessment has been investigated. In fact, in many cases, the buildings such as the case study structure are intended for public activities such as museums, so specific fire requirements, like fire resistance, are necessary. This topic became relevant especially if the structure is equipped with particular structural retrofit interventions which can be altered and modified in case of a fire. The paper presents the results of advanced thermomechanical analyses on the historical masonry palace under investigation. Since the case study building has a masonry structural type at the first two floors while there is a third floor realized in reinforced concrete, the fire analyses were conducted on the third and fourth floors, which may be more vulnerable to fire.

The assessment and retrofit of existing masonry structures with historical and cultural value in highly seismic zones are challenging issues in earthquake engineering. In fact, the historic and recent earthquakes have shown the problem of the seismic vulnerability of existing masonry constructions. A historical masonry palace located in Caggiano (Salerno, Italy) is used herein as a case study, showing the vulnerability assessment and the seismic upgrading process. The case study building has a masonry structural type at the first two floors while there is a third floor realized in reinforced concrete and a fourth floor realized with a wood structure. The building was characterized by a remarkable seismic vulnerability and needed seismic upgrading operations. After the vulnerability assessment process, some design suggestions are proposed for the seismic upgrading of the building. The structure before and after the upgrading operations has been checked through nonlinear static and dynamic analyses. Then, coherently with the "Sismabonus" approach, the attribution of the seismic risk class, performed through numerical analyses, is founded on two parameters, namely, the expected annual mean losses (PAM), related to economic factors, and the Life Safety Index (IS-V), related to the structure seismic safety. Finally, the overcoming of the different classes of risk is shown and compared with the amount of the retrofit operations, their costs, and the impact on the existing space. Moreover, fire assessment has been investigated. In fact, in many cases, the buildings such as the case study structure are intended for public activities such as museums, so specific fire requirements, like fire resistance, are necessary. This topic became relevant especially if the structure is equipped with particular structural retrofit interventions which can be altered and modified in case of a fire. The paper presents the results of advanced thermomechanical analyses on the historical masonry palace under investigation. Since the case study building has a masonry structural type at the first two floors while there is a third floor realized in reinforced concrete, the fire analyses were conducted on the third and fourth floors, which may be more vulnerable to fire.

INTRODUCTION
The presented project involves the recovery of Palazzo Colonna and the surrounding areas, of which the Palazzo is the central fulcrum. The project involves a real urban redevelopment aimed at both the recomposition of this piece of the ancient center with the adjacent areas already previously requalified and the creation of safe spaces with a view to recovering the historical image of the site and its safety in case of calamitous events. Within the urban redevelopment project, the Palazzo Colonna is a true dynamic fulcrum, becoming a sort of incubator for cultural activities, strictly connected to the socioeconomic relaunch of the entire community. The surveys showed situations of extreme degradation of some vertical structures, whose roofs and internal horizontals are largely collapsed. The opportunity to make upgrading interventions allows a better use of the entire building sector in conditions of maximum safety. Therefore, the project includes the structural and functional upgrading of the building aimed at creating exhibition spaces, conference rooms, rooms for coworking, and services for the community in general. This work presents the structural assessment and upgrading design for the case study building, considered together with the goal of the functional upgrading of the whole structure. The assessment and upgrading of historical existing buildings in seismic regions are fundamental to reduce the associated risk. The vulnerability assessment of historical masonry structures has been investigated in the literature through different methods and focusing on different modeling issues, issues such as presented in Giordano et al. (2019), Chieffo et al. (2019a), Clementi et al. (2020), Ferrante et al. (2019), and de Silva (2020), to name a few.
Then, after the understanding of the failure mechanics, the best upgrading strategy can be chosen for the structure as a compromise among different needs (structural, architectural, economical, etc.). The criteria for the choice of the retrofit methods have been highly investigated in the literature (Moehle, 2000;Thermou and Elnashai, 2006;Formisano et al., 2008Formisano et al., , 2016aFormisano et al., , 2017aCalvi, 2013;Formisano and Sahoo, 2015;Terracciano et al., 2015;Miano et al., 2017b;Chieffo et al., 2019b). These techniques are usually implemented to evaluate the improving of the structural performance and also to evaluate and minimize the costs and the losses (Aslani and Miranda, 2005;Jalayer et al., 2012Jalayer et al., , 2015Liel and Deierlein, 2013;Miano et al., 2018aMiano et al., , 2019a. Herein, the Pushover Analysis has been used to assess the performance of the building with respect to vertical and seismic loads and to provide an estimation of the vulnerability index (Frascadore et al., 2015;NTC, 2018). The nonlinear dynamic analysis procedures have also been implemented to implement the seismic assessment. Cloud Analysis is selected herein. It is based on a simple regression in the logarithmic space of structural response vs. seismic intensity using a non-scaled ground motions set. The mathematical formulation is very simple, and the application is very quick and efficient for fragility assessment (Jalayer et al., 2007(Jalayer et al., , 2015Celik and Ellingwood, 2010;Miano et al., 2017a). The seismic assessment of the case structure has shown deficiencies in seismic and vertical loads conditions. The upgrading techniques have been thought to reduce the criticisms by solving the vertical loads deficiencies and reaching a higher level of the vulnerability index. The technique used to improve the structural performance of the masonry is the "Cuciture attive per la muratura" (CAM) technique. Moreover, in the third floor, RC jacketing is used for the columns, while steel jacketing is used for the beams. The vulnerability index values for the limit state of Life Safety (LS) before and after upgrading are estimated. The fragility curve for the limit state of LS is calculated before and after upgrading. Then, a comparison between the results from pushover and Cloud analyses is proposed. The seismic risk class is calculated before and after the upgrading (Guidelines, 2017;Cosenza et al., 2018). The assigned class is a function of two parameters, e.g., the expected annual mean losses (Perdite Annue medie, PAM) and the Life Safety Index (Indice di sicurezza della vita, IS-V).
Moreover, fire assessment has been investigated. In fact, in many cases, the buildings such as the case study structure are intended for public activities such as museums, so specific fire requirements, like fire resistance, are necessary. This topic became relevant especially if the structure is equipped with particular structural retrofit interventions which can be altered and modified in case of a fire. The paper presents the results of advanced thermomechanical analyses of the historical masonry palace, located in Caggiano (Salerno, Italy), hosting a museum. Since the case study building has a masonry structural type at the first two floors while there is a third floor realized in reinforced concrete and a fourth floor realized with a wood structure, the fire analyses were conducted on the third and fourth floors, which may be more vulnerable to fire.

Building Description
As previously highlighted, the project involves the recovery of Palazzo Colonna and the surrounding areas. The structure is composed of two masonry floors and a reinforced concrete third floor. The building is composed of two blocks placed around a central courtyard characterized by the presence of a central monumental staircase which leads to the main floor. Figure 1 presents the building south view.
The Palace is today morphologically altered by the restoration interventions that have turned it into the town house, after the 1980 Earthquake. These operations have irreversibly compromised the original 18th-century morphology of which the only remains are concentrated in the inner courtyard where there are still visible remains of the decorative apparatuses for which restoration is planned. For the restoration of the decorations, a careful mapping of the identified and classified degradation has been drawn up according to the categories codified by the UNI 11182 (2006) standard. The surveys showed situations of extreme degradation of some vertical structures, whose roofs and internal horizontals are largely collapsed. Figure 2 shows  the degradation analysis results for one of the frontal views of the building.
Moreover, particular attention was dedicated to the plant system upgrading. All the solutions respected the fundamental principles of containment of consumption, modularity, maintainability, and reversibility. Then, regarding the respect of the legislation for the overcoming of architectural barriers, all the provisions of the current legislation have been satisfied. Finally, Figure 3 shows a typical floor view of the building. It is possible to see that different thicknesses of the masonry walls are present in the building, varying from about 0.5-1.0 m. The stairs have been realized in masonry. Finally, the slabs are done with a composite system of steel joists and thin hollow bricks and cast in place of concrete.
With respect to the characterization of the mechanical material properties, Table C8.5.1 of the Circolare (2019) has been used. The mechanical properties of the masonry related to "irregular masonry of soft stone (tuff, calcarenite, etc.)" have been chosen. Conservatively, while waiting for specific destructive/non-destructive tests on materials, the minimum value in the range proposed by the code has been taken. The used mechanical material properties are described as follows: • fm = average compressive strength of the masonry: 140 N/cm 2 • τ 0 = average shear strength of the masonry: 2.8 N/cm 2 • E = average value of the Young modulus: 900 N/mm 2 • G = average value of the tangential elasticity modulus 300 N/mm 2 • w = average specific weight of the masonry: 16 kN/m 3 Finally, it is to note that in the upgrading design, the mechanical properties values have been improved based on the designed operations on the masonry, such as cracks injections, joints restillation, and the scuci-cuci method of patching.

Introduction to the Assessment and Upgrading Operations
The technique used to improve the structural performance of the masonry is the CAM technique. The CAM consists of a "packing" of the masonry with stainless steel strips arranged in the horizontal and vertical directions, passing through the masonry thickness and closed on themselves after applying a pretension. The adoption of a system of tensioning diffused in the three orthogonal directions, in particular also in the transversal direction, improves the monolithic and the mechanical behaviors of the masonry body, increasing its shear strength and bending in plane and out of plane. Through the special connection elements (stainless steel plates, equipped with a funnel-shaped hole), the stainless steel strips allow to create a continuous system of trusses, able to retrace all the irregularities of the masonry, both horizontally and vertically. Moreover, in the third floor, RC jacketing is used for the columns, while steel jacketing is used for the beams.

Nonlinear Static Analysis Procedure
Nonlinear static analysis (e.g., Pushover Analysis) is proposed herein. The CDS (CDS, 2018) software has been used. The modeling of masonry structures is strictly related to the choice of modeling single units or more complex aggregates (Chiumiento and Formisano, 2019). This issue is out of the scope of this paper. Moreover, the identification and calibration of the structural model of masonry buildings are fundamental issues (for example, see Catinari et al., 2017;Clementi et al., 2017, to name a few), but are outside the goals of this work. The structure has been modeled using the finite element method (FEM), inserting two types of elements: (1) one-dimensional line element for RC members and (2) the two-dimensional shell (quad) element for the masonry structure. With respect to masonry structure, initially, a first model schematizes the masonry walls with rectangular two-dimensional elements representative of singlemembrane stresses with deformability in their own plane, as plate, therefore, with shear, extensional, and flexion deformability typologies. Then, a second model known as the "equivalent frame, " schematizes the panels and the floor strips with beamtype finite elements with two nodes and rigid offsets to model the panels or the floor band. With respect to RC line members, their force-deformation behavior is evaluated considering that, for each structural member, inelastic action is considered at the FIGURE 4 | Pushover curves before upgrading operations in terms of (A) base shear/displacement and (B) acceleration/displacement. Pushover curves after upgrading operations in terms of (C) base shear/displacement and (D) acceleration/displacement. ends while the middle remains elastic. The nonlinear momentcurvature relationship in the element ends is obtained through section analysis. The safety verifications are proposed using the limit state method with the partial coefficients (NTC, 2018). LS limit state results are presented herein, while in the design phase, also serviceability limit states have been verified. The LS limit state is achieved when the maximum shear strength (brittle safety verifications) or the maximum rotational capacity (ductile safety verifications) is reached for a structural element. The safety checks for the joints are also considered. Pushover Analysis is implemented based on the code recommendations (NTC, 2018). Vertical loads analysis is done before Pushover Analysis. Then, 32 pushovers are implemented, using two force distributions with an eccentricity of 5% in each direction. The first distribution of the forces is proportional to the floor masses. The second one is proportional to the shape of the first vibration mode (NTC, 2018;Circolare, 2019). Figures 4A,B presents the pushover results (base shear or acceleration vs. top displacement). The 32 pushovers are presented by color in groups of four. Each group corresponds to a certain direction of the forces application ±30% of the orthogonal force. The four pushovers for each direction refer to the two force distributions with the application of ±5% eccentricity associated to the structure geometrical barycenter. The red curves refer to forces parallel to the horizontal positive X-Direction (Figure 3); the blue curves refer to forces parallel to the horizontal negative X-Direction; the black curves refer to forces parallel to horizontal positive Y-Direction; the green curves refer to forces parallel to horizontal negative Y-Direction. For the LS limit state, the verifications are not satisfied. Figure 5A presents for the LS limit state with red color the elements for which the safety checks are not verified and with green color the ones for which the safety checks are verified. A consisting portion of masonry panels, columns, and beams are not verified. The failures of the masonry panels are mainly axial-bending crises. RC columns and beams suffer premature flexural and shear failures. The choice and the investigation of the right retrofit methods are very discussed topics inside the scientific community (Moehle, 2000; Thermou and Elnashai, 2006;Calvi, 2013;Formisano and Mazzolani, 2015;Formisano et al., 2016bFormisano et al., , 2017bMiano et al., 2019b). Herein, the choice has been as a compromise between architectural and structural recommendations, remaining in the economical funding established for this design project. Then, the same procedure has been repeated after the implementation on the upgrading operations. The same modeling approach presented before is used for the upgraded structure. Pushover Analysis has been implemented. A set of 32 pushovers has been implemented using two different distributions of forces and using a 5% eccentricity in each direction. Figures 4C,D present the pushover results (base shear or acceleration vs. top displacement). The red curves refer to forces parallel to horizontal positive X-Direction (Figure 3); the blue curves refer to forces parallel to horizontal negative X-Direction; the black curves refer to forces parallel to horizontal positive Y-Direction; the green curves refer to forces parallel to horizontal negative Y-Direction. Figure 5B presents the LS limit state with red color and the elements for which the safety checks are not verified and with green color the ones for which the safety checks are verified. The number of unverified elements strongly decreases after the upgrading operations. It is to note that some panels are still unsafe, while all the premature failures for the RC members are solved.

Vulnerability Index Calculation
The Italian code (NTC, 2018) recommends to verify the seismic vulnerability before and after upgrading with respect to LS limit state. This process has been followed for the case study structure.
The calculation of the vulnerability index (ζ E ) is defined and discussed in the literature (Frascadore et al., 2015;Cosenza et al., 2018). The ζ E is an efficient parameter to evaluate the vulnerability of the building. In this section, a brief summary of the procedure to calculate ζ E is presented. For more details, see Frascadore et al. (2015). Based on the pushover response for the Multi Degree of Freedom (MDOF) system, it is possible to have the response for the corresponding Single Degree of Freedom (SDOF). The SDOF is characterized by the period T * , yield strength F * y , and ultimate displacement d * u . The calculation of the capacity in terms of return period and then of peak ground acceleration (PGA), for which the failure is achieved, is based on these parameters. The quantification of the ζ E is obtained in the Acceleration Displacement Response Spectrum (ADRS) space (Fajfar, 1999(Fajfar, , 2000, in which the abscissas are the spectral displacements and the ordinates are the spectral accelerations. The elastic spectrum has to be scaled until the spectrum containing the point performance (S ae ; S de ) of the SDOF is intercepted. It is identified by the inclination line T * and the displacement d * maxSLV . The ζ E is the ratio between the demand PGA (based on the code recommendations for the LS limit state) and the capacity PGA of the structure: where PGA LS_Capacity is the PGA corresponding to the achievement of the first failure at the LS limit state inside the structure, while PGA LS_Demand is obtained from the elastic code spectrum for the specific site for the LS limit state. The

Cloud Analysis and Fragility Estimation
The estimation of the analytical fragility curves through nonlinear dynamic analyses can be very useful to have a more complete measure of the vulnerability. In this section, the fragility curves are presented for the investigated structure both before and after the upgrading operations. A group of 30 ground motion records has been applied to the SDOF previously described (e.g., among the 32 SDOF systems, the most penalized SDOF, that minimizes ζ E , has been used). A brief summary of Cloud Analysis is proposed. More details about the application of the Cloud Analysis for SDOF systems are provided in Jalayer et al. (2019). The first step of the procedure is to identify the structural response parameter. As presented in Jalayer et al. (2007), the critical demand to capacity ratio for the LS limit state (DCR LS ), equal to the mechanism that brings the structure closest to the onset of the limit state, is chosen as the structural response parameter. It can be calculated for each nonlinear time-history analysis. The DCR LS parameter, equal to one at the achieving of the LS limit state, is obtained for the SDOF system as the ratio between the maximum demand in terms of displacement for each record among all the steps of the nonlinear dynamic analysis and the limit state capacity, calculated based on Italian code (NTC, 2018). The LS capacity corresponds to the ultimate displacement of the bilinear system of the SDOF. The ratio between demand and capacity for SDOF and MDOF is equal because both the demand and the capacity are multiplied by the modal participation factor to pass from SDOF to MDOF. Herein, LS limit state is considered, but the procedure can be used for other limit states. This choice allows to have a comparison with the trend of ζ E before and after the upgrading operations. The record selection is important to have a proper implementation of Cloud Analysis. A set of 30 records is chosen from both the NGA-West2 and RESORCE database (see Miano et al., 2018b, andMiano et al., 2019a, about the record selection recommendations). This records set has a range of magnitudes between 5.2 and 7.2, and closest distanceto-ruptured area (R RUP ) up to around 50 km. The shear wave velocity in the upper 30 m of soil, Vs30, at the structure's site is in higher than 800 m/s (e.g., soil class A, accordingly to Italian code, NTC, 2018). Then, the selected records are mainly on NTC 2018 site classes A and B. Only one of the two horizontal components of each record is considered. The records are free field or measured on the ground level. Table 1 presents the records set.
The structural fragility is estimated through Cloud Analysis. It considers a structure subjected to a group of unscaled records of specific values of the first-mode spectral acceleration Sa(T), adopted herein as the intensity measure (IM) parameter. The records are applied to the structure, and the DCR LS is estimated.
This gives a set of values that form the basis for the Cloud Analysis calculations. The statistical properties of the Cloud response are calculated through a logarithmic linear regression applied to the response. This is equivalent to fitting a power-law curve to the cloud response in the original (arithmetic) scale. The result is a curve that finds the median drift demand for a given level of Sa: where ln(a) and b are the constants of the regression. The logarithmic standard deviation β DCR|Sa is the root mean sum of the square of the residuals: where DCR i and Sa, i are the DCR values and the corresponding Sa for record number i within the cloud response set, and N is total number of records. In conclusion, the fragility curve based on the Cloud Analysis is: Figures 6A,B shows the Cloud regression results for the LS limit state with reference to the pre and post upgrading structure. The scatter plots for Cloud Analysis data D = {(S a,i , DCR LS , i ), i = 1:30, for the group of records outlined in Table 1, are shown in Figures 6A,B. For each data point (cyan-colored squares), the corresponding record number is shown. The figures illustrate also the Cloud Analysis regression prediction model (i.e., regression line and the estimated parameters a and b) fitted to the data. The line DCR LS =1, for which LS limit state is achieved, is proposed with a red dashed line. Figure 6C shows the resulting Cloud Analysis-based fragility curve with reference to the bare and the upgraded building. The values of the median and the logarithmic standard deviation of fragility curves are, respectively, 0.25 and 0.32 before upgrading and 0.40 and 0.39 after the upgrading.

CALCULATION OF THE CLASS OF RISK BEFORE AND AFTER THE UPGRADING
The class of risk (as suggested in the Italian Guidelines on the seismic classification, Guidelines, 2017 andL.205/17, 2017) is evaluated for the investigated structure before and after the upgrading based on two parameters derived after the assessment process of the building. The economical parameter is the expected annual mean losses (PAM). The structural safety parameter is the Life Safety Index (IS-V). The 2017 Italian law (Balance Law 2017) proposed a campaign for the seismic upgrading of existing buildings (e.g., "Sismabonus" plan). This plan stimulated the evaluation of the existing building and the prevention of the seismic risk. In relation to this, the Italian guidelines provided operational tools for the seismic risk classification of the structures. Eight risk classes, with increasing risk from letter A+ to G, are considered. Two methods for evaluating the class of risk for a structure are considered: the conventional and the simplified methods (Guidelines, 2017). The conventional method is presented for this case study. The steps for determining the class of risk are the following: • PGA evaluation: the PGAs at the building site, PGA D , for achieving the various limit states are estimated; • Analysis of the structure: the capacity PGA values, PGA C , for which LS and damage limitation (DL) limit states are reached, are estimated; • IS-V class definition, as relation between PGA C (for LS limit state) and PGA D of the site. The percentage value, calculated based on Table 2.1 defined by the guidelines, defines the seismic risk class as a consequence of vulnerability index; • Vulnerability analysis: it defines the damages (structural and not structural) based on the achievement of predefined levels by the response parameters, calculating return periods, T rC , and its inverse annual average frequency of exceeding λ; reconstruction cost percentage (CR%) associated with λ values for the limit states considered; • PAM class definition as the area under the curve of the direct economic losses, calculated as function of the annual average frequency of exceeding the events that produce the achievement of a predefined limit state. The PAM Class definition is based on Table 2.2 of the guidelines. Then, the risk class is defined as the worst class between the PAM and the IS-V classes. Figure 7 shows risk class evaluation pre and post upgrading.

FIRE RESISTANCE ASSESSMENT
Structural fire safety is one of the primary aspects that has to be considered in the design of buildings especially where the intended of use involves particularly fire-risky activities (e.g., elevated fire loads) or the structure materials are vulnerable to fire (e.g., steel or wood). Indeed, given its nature, the most vulnerable material to fire is steel, where sometimes the passive fire protections are necessary (de Silva et al., 2017(de Silva et al., , 2019, while the fire performance of the masonry and concrete has often been taken for granted considering its noncombustible nature and ability to function as a thermal barrier, preventing heat and fire spread. The low concrete thermal conductivity, compared to steel, cause inside the fire-exposed concrete members, strong temperature gradients and the core region may take a long time to heat up. Thus, while the compressive strength of concrete is rapidly lost beyond a critical temperature, which is not too dissimilar to the equivalent temperature for loss of steel strength, structural effectiveness is not affected until the bulk of the material reaches the same temperature. This requires an analysis of the thermal response of the entire structural element (Fletcher et al., 2007). At present, structural fire safety (fire resistance) of members is generally achieved through prescriptive approaches, which are based on standard fire resistance tests or empirical calculation methods. These prescriptive based approaches (PBAs) are design approaches that will allow the building to meet a preset standard when it comes to fire safety. Applying already deemed to satisfy solutions to the design of the building will inherently be safe, it is a code-compliant building. Therefore, a building that is designed in accordance with the code is considered an acceptable building from a fire risk perspective.
Sometimes, the building deviates from a "standard" building, and it is very difficult to impossible to apply the standard set of rules set out in the code. This of course does not necessarily mean that the building is unsafe; it just means that it cannot comply with the deemed to satisfy solutions. In these cases, the Fire Engineering Approach (FEA) could be useful for the evaluation of the real level of building fire safety. Indeed, the recent move toward performancebased fire design, by using the Fire Safety Engineering (FSE) criteria, has increased the focus on the use of advanced computer simulations for evaluating fire resistance of structure. In particular, the FEA consists of detailed analysis of the fire, considering natural fire curves, which combine more sophisticated calculation (advanced methods) for structural models and the temperature distribution in the elements, the mechanical and geometric nonlinear structural responses are taken into account.
This requirement, implemented in the National Codes of European member countries, is explained by achieving the following five performance levels: -the load-bearing capacity of the construction can be assumed for a specific period of time -the generation and spread of both fire and smoke within the works are limited -the spread of fire to neighboring construction works must be limited -occupants have to be able to leave the works or be rescued by other means -the safety of rescue teams must be taken into consideration.
The results of each application of the FEA to the fire safety should be evaluated through the analysis of the achievement of these objectives. Figure 8 explains the steps for applying the FEA.
In the following, the PBA and the FEA were applied to the case study of Palazzo Colonna in Caggiano. In particular, since the main structure of Palazzo Colonna in Caggiano is in masonry, the concrete part of the structure is analyzed, with the aim of observing all the effects that may arise in the structural elements during fire.

Fire Modeling
The fire event represents a very complex phenomenon involving different interacting physical and chemical processes, such as combustion, radiation, and transport phenomena under turbulent conditions. Fire evolution and the severity of the consequences depend on the boundary conditions, such as the confined or unconfined environment. Given the complexity of the scenarios, the simulation tools are becoming fundamental for design purposes of buildings in case of fire. Indeed, numerical modeling approaches allow to obtain cost-effective and reproducible information, with a degree of confidence dependent on the quality of the input information and the assumptions made by the model itself. For a building-confined fire, like in this case, the most accredited models are the zone models. Zone models are numerical tools used for the evaluation of the gases temperature within a compartment during a fire, and they are based on several hypothesis. The main hypothesis in zone models is that the compartment is vertically divided into zones (Figure 9), and each zone has homogeneous properties in terms of temperature and composition: a hot layer with combustion products, located near the ceiling, and a cold layer with fresh clean air at the bottom, separated by a moving interface, which is considered an adiabatic surface. The zone properties (and the layer height) can vary over time according to global conservation equations (mass and energy). In particular, in one zone model (Figure 9B), the temperature is considered to be uniform within the whole compartment. This type of model is thus valid in case of generalized fires, contrary to two zones models, which are valid in case of localized fires.
These models are implemented in the software (Cadorin and Franssen, 2003), which is used for the simplified fluid-dynamic analyses in this case study. In the following, all the details are described.
Performance Level, Fire Scenarios, and Fire Models The main objective of fire safety checks concerns the mechanical resistance and stability of the building in case of fire. In the case of PBA, the required time of fire resistance is 90 min of standard fire curve, while for the FEA, the safety performance level required for the structure is maintaining the fire resistance requirements for the duration of the natural fires. For the definition of the natural fire curves, the description of the structure fire scenarios is necessary. The design fire scenario is a qualitative description of the fire development during the time, identifying key events that characterize the fire (Del Prete et al., 2016). Since each floor has a different occupancy like museum and offices, the scenario is chosen to maximize the fire risk based on engineering judgment. In particular, the analyzed fire scenario consists of the library located at the third floor of the structure. Furthermore, this part of the structure is in concrete, while the other parts are in masonry.
The library was considered a fire compartment of 36 m 2 , and the post-flashover fire is modeled by two-zone model, which assumes homogeneous temperature, density, internal energy, and pressure of the gas in the compartment. The combustion model was considered as "extended fire duration, " which supposes that the release of mass may be limited by the quantity of oxygen available in the compartment. The total mass of fuel is burnt inside the compartment (safe procedure) then the fire duration is increased compared to the input one (Cadorin et al., 2001). In this model, no external combustion is assumed, the entire fire load delivers its energy inside the compartment. If the fire is ventilation controlled, the pyrolysis rate is proportional to the oxygen coming in the library and is not a physical model because pyrolysis is not directly dependent on oxygen concentration. It has been established for design procedures in order to avoid uncertainties on the maximum pyrolysis rate per unit floor area and therefore to be on the safe side concerning the fire duration. In this case, the fire load was estimated equal to 500 MJ/m 2 of paper. Figure 10A shows the data RHR, obtained in accordance with EN 1991EN -1-2 (2002 and the RHR computed obtained as a result of Ozone. Figure 10A shows that the fire is ventilationcontrolled, indeed, the fire is limited by the availability of oxygen in the compartment and the RHR computed has a lower stationary phase than the RHR data. This is confirmed also by the comparison between the ventilation factor (O comp ) of the compartment and the critic ventilation factor (O cr ) (Compagnone et al., 2017): O comp = 0.021 m 1/2 < O cr = 0.063 m 1/2 (5) Equation (5) shows that O comp < O cr ; this is the case in which the ventilation controls the fire. Figure 10B shows the natural fire curve of the library compared with the standard fire curve, which gives higher temperatures and therefore more demanding for the structure.

Thermomechanical Analyses
As said before, the thermomechanical analyses purpose concerns the mechanical resistance and stability of the building in fire. In the case of PBA, the required time of fire resistance is 90 min of standard fire, while for the FEA, the safety performance level required for the structure is maintaining the structural stability for the duration of the library natural fire. In the case of PBA, the simplified tabular method is considered, while applying the FEA criteria, advanced thermomechanical analyses were conducted on the concrete part of the building, selected as significant substructure (Figure 11).

Prescriptive Based Analyses Results
The PBA approach was considered applying the simplified tabular data according with Table 5.2a and Table 5.5 of the Eurocode (EN 1992-1-2, section 5). Table 2 shows that the beam "B1" does not satisfy the required resistance time of 90 min, while for the other elements, both the b min and a dimensions, allow to reach 90 min.
Since the tabular data are often very conservative, an analysis of the substructure subjected to the ISO 834 curve was carried out in order to check the eventual differences and to take into account the structural thermal effects, which are neglected in the tabular method. The results are shown in Figure 12. Figure 12 shows the variation of the bending moment during the thermal transitory due to the statically redundancies and restraints to thermal expansions, which causes the emergence of additional stresses in structural elements. As Table 2 shows, the only element which does not satisfy the 90 min of fire resistance is the beam "B1" with tabular data; furthermore, the same beam satisfies the request if modeled within the entire structure always with ISO834 curve, with a sufficient margin of safety ( Figure 13A); this is the demonstration of the conservative of the simplified tabular data. Furthermore, Figure 13 shows that, on the one hand, the resistance of the beam decreases with the increase of time (because the temperature is increasing), and on the other hand, the stress bending moment increases due to the structure redundancies.

Performance-Based Analyses Results
In order to evaluate the real level of fire safety of the concrete substructure, an advanced thermomechanical analysis was conducted with the library fire curve, checking the stability of the substructure for the entire duration of the library fire (about 300 min).
In addition, in this case, there is a variation of the bending moment during the thermal transitory due to the structural redundancy, but without causing the structural collapse. The beam "B1, " which is the weakest element in the tabular data analyses, is verified also in this case, with a greater safety margin than the case of the ISO834 fire curve (Figure 13B).

CONCLUSIONS
The assessment and retrofit of existing masonry structures with historical and cultural value in highly seismic zones are challenging issues in earthquake engineering. A historical masonry palace, located in Caggiano (Salerno, Italy), currently hosting a museum of the country culture of the town and previously a location of the town municipality, is used herein as a case study to show the vulnerability assessment and the seismic upgrading process for the existing masonry structures with historical and cultural value. The case study building has a masonry structural type at the first two floors, while there is a third floor realized in reinforced concrete and a fourth floor realized with a wood structure. The masonry palace has a configuration of the structural system similar to the big part of the masonry buildings of the historical center of different Italian towns. In fact, most of these Italian masonry buildings were built in a time period where the knowledge about the seismic behavior of structures, performance of materials, and earthquake actions was very limited. As a consequence, these are characterized by a remarkable seismic vulnerability and need seismic upgrading operations compatibility with existing spaces and the "Italian authority for the preservation of cultural heritage" recommendations. After the vulnerability assessment process, some design suggestions are proposed for the seismic upgrading of the building. The structure before and after the upgrading operations has been checked through nonlinear static and dynamic analyses. Then, coherently with the "Sismabonus" approach, the attribution of the seismic risk class, performed through numerical analyses, is founded on two parameters, namely, the expected annual mean losses (PAM), related to economic factors, and the Life Safety Index (IS-V), related to the structure seismic safety. Finally, the overcoming of the different classes of risk is shown and compared with the amount of the retrofit operations, their costs, and the impact on the existing space. Moreover, fire assessment has been investigated. In fact, in many cases, the buildings such as the case study structure are intended for public activities such as museums, so specific fire requirements, like fire resistance, are necessary. This topic became relevant especially if the structure is equipped with particular structural retrofit interventions which can be altered and modified in case of a fire. The paper presents the results of advanced thermomechanical analyses of the historical masonry palace, located in Caggiano (Salerno, Italy), hosting a museum. Since the case study building has a masonry structural type at the first two floors while there is a third floor realized in reinforced concrete, where a library is located, the fire analyses were conducted on the third floor, which may be more vulnerable to fire, making the fire load considerable. The structural fire analyses show that the simplified tabular analyses are conservative, indeed, the beam, which does not satisfy the required R90, is verified if an advanced analysis of a substructure is conducted, referring to the ISO834 fire curve.
The safety check is even more precise in the case of a natural fire curve analysis, where the library fire curve is considered. In this case, the safety margin is greater than the ISO834, since the input temperatures and therefore the temperatures in the structural elements are, on average, lower. This result underlines the efficacy of using the PBA, which certainly returns more reliable results.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.