Integrated Cost-Analysis Approach for Seismic and Thermal Improvement of Masonry Building Façades

The combination of structural and thermal efficiency is a new frontier in civil engineering. Indeed, the retrofitting strategies should optimize costs and technical solutions from these two points of view. If a technical solution is able to provide an improvement of both structural and energetic behavior, then the utility of the intervention can better justify the economic investment. In this paper, a meso-scale approach (i.e., façade-scale) for integrated interventions applied on masonry façades is proposed. The structural performance of the façade is evaluated by considering base shear and ductility of the structural element through non-linear static analyses. Moreover, the thermal indicator, that is the thermal transmittance, is computed with a simplified approach in terms of an equivalent wall, taking into account the role of the windows and doors of the façade. As proof of concept, the procedure is applied to a façade of an existing masonry building. Economic and environmental iso-cost curves are obtained to tune the interventions conceived for a real case study, analyzing the benefit offered by different retrofitting solutions.


Introduction
Nowadays, existing masonry buildings do not fulfill standards requirements in terms of seismic and energetic performances according to current regulations. In the past, masonry buildings were generally not designed to withstand seismic actions, due to improper vertical connections between walls and horizontal connections between more or less flexible diaphragms [1] or roofs and walls [2,3]. These aspects worsen the structural behavior in case of irregular in-plan buildings [4,5]. In addition, in some cases, the masonry texture is such that localized damages can occur, making this aspect extremely difficult to evaluate in a seismic analysis. For these reasons, it is necessary to adopt structural techniques, such as ferro-cement, reinforced plaster, grout and epoxy injections, steel plates, stainless reticulatus grids, FRP, GFRP [6], bionatural aggregates and others. These techniques improve the in-plane and out-of-plane behavior of masonry walls [7][8][9][10], either undamaged or damaged, giving a more monolithic behavior, greater strength and stiffness. Moreover, the level of knowledge of the material is a key aspect, and for that, only experimental and in situ tests can really be reliable [11][12][13]. As for the thermal aspect, as is well-known, masonry does not have high performance levels in terms of thermal insulation, in case of absence of insulating layers. Although insulating from the exterior side would provide optimum performance and would ensure the integrity of the masonry walls, in the case of historic stone buildings, the addition of new insulation is generally restricted to the interior side of the building for aesthetic reasons [14]. Nevertheless, stone masonry was introduced in [27,28] and applies the procedure to an existing masonry building discussing the different implications of each assumption. It also gives practical recommendations on how to evaluate the level of improvement in an integrated way.
The three levels of investigation (single panel, single façade, entire building) can provide increased levels of information to decide the most appropriate retrofitting strategy.
Section 2 illustrates the integrated approach at a mesoscale level (i.e., the masonry façade), which is intermediate between the microscale level (i.e., the single masonry panel) and the macroscale level (i.e., the entire masonry building). Section 3 applies the procedure to a façade of a case study of an existing masonry building located in Italy. Section 4 defines the demand curves, discusses the results and generalizes them in a larger perspective

Procedure
The mesoscale approach proposed here consists in performing structural and thermal analysis on external façades of buildings, to evaluate the economic and environmental iso-cost curves resulting from the adopted interventions. The latter aspect has the purpose of understanding which type of intervention is optimal from a structural, sustainable and economic point of view.
The procedure consists of the following steps [ Figure 1]: 1. Acquisition of the structural and thermal parameters needed for the analysis (mechanical parameters, such as tensile strength, shear strength, etc. and thermal parameters, such as thermal conductivity, thicknesses, etc.) 2.
Definition of a set of integrated interventions, namely retrofitting strategies that have positive effects on either the seismic response, or the thermal performance or both of them.

3.
Identification of a performance indicator at mesoscale level: as for the structural behavior, the variation of base shear capacity ∆V and the corresponding variation of ductility capacity ∆µ are considered, defined by well-known methods (non-linear static analysis); as for the thermal side, the variation of thermal transmittance ∆U is taken into account.

4.
Economic and environmental iso-cost curves representing the relationships between the thermal capacity indicator (∆U) and the structural capacity indicators (∆V or ∆µ). For each integrated intervention, after an economic budget (investment) or an environmental impact in terms of CO 2eq are fixed, one can calculate, as explained in Sections 2.2 and 2.3, the corresponding pair of capacity indicators. That represents a single point in the graphs ∆U−∆V or ∆U−∆µ. By varying the economic budget or the environmental impact, several points are obtained. Finally, the curves fitting these points can be regarded as iso-cost curves. Moreover, the iso-performance curves, that will be the subject of future work, will express, for the same seismic or energetic performance, the economic investment needed or the environmental impact caused by each integrated intervention.

5.
Definition of dimensionless parameters c U and c R , defined in [28], to identify the demand for thermal and seismic performances respectively, with the expressions: where PGA i is the seismic peak ground acceleration at the site, PGA max is the maximum value in all the Italian regions (or any other reference zone in the world); DD i is the degree day value of the site and DD max the maximum Italian value. They are given for the different construction sites respectively by [30,31].
Buildings 2020, 10, 143 4 of 18 6. The hypothesis of a correlation between energy efficiency demand and seismic demand, through demand curves identified by these analytical expressions: They are based on a simple criterion of proportionality of increment of capacity indicators (∆ , ∆ or ∆ ) with the demand indicators ( , ). As stated in [28], is a corrective and tuning parameter that can be defined by the decision-maker according to political, social or stakeholder reasons. The values of this parameter α can be defined by using a short number of case studies that are assumed as emblematic to the political, social or stakeholder goals. The tuning parameter α should assume higher values in low seismic and high thermal demand areas. By contrast, in high seismic and low thermal demand areas, the α values may be lower. One can see α as an additional parameter that modifies the relationship between the thermal and seismic indicator according to necessities out of technical requirements. Without any political need, the tuning parameter can be assumed equal to 1.
Target points for each site can then be obtained by intersecting the capacity curves (obtained with step (4)) and the demand curves (obtained with step (5) and (6)).

Thermal Performance Indicator
From a thermal point of view, all the structural and non-structural elements are subjected to a heat flux that causes thermal dispersions. The main parameter to evaluate these dispersions is the thermal transmittance , defined as the rate of transfer of heat through a unitary surface with a difference of temperature of 1 °C. For its calculation, the UNI EN ISO 6946:2008, [32] gives the following expression: is the thermal resistance of the internal surface, m K/W ; is the resistance of the i-th layer, m K/W , where s is the thickness of the i-th wall layer m ; is the resistance of the external surface, m K/W .

The term
represents the thermal conductivity of the element and measures the tendency of transmitting heat, which can be found in UNI 10351:1994 [33].
At the same time, the variation of thermal transmittance Δ was adopted as a reference parameter to evaluate the thermal variation due to the intervention. Therefore Δ has been expressed by the ratio between the performance variation (the difference between the value after ( ) and before the integrated intervention ( )) to the initial value : They are based on a simple criterion of proportionality of increment of capacity indicators (∆U, ∆V or ∆µ) with the demand indicators (c U , c R ). As stated in [28], α is a corrective and tuning parameter that can be defined by the decision-maker according to political, social or stakeholder reasons. The values of this parameter α can be defined by using a short number of case studies that are assumed as emblematic to the political, social or stakeholder goals. The tuning parameter α should assume higher values in low seismic and high thermal demand areas. By contrast, in high seismic and low thermal demand areas, the α values may be lower. One can see α as an additional parameter that modifies the relationship between the thermal and seismic indicator according to necessities out of technical requirements. Without any political need, the tuning parameter can be assumed equal to 1.
Target points for each site can then be obtained by intersecting the capacity curves (obtained with step (4)) and the demand curves (obtained with step (5) and (6)).

Thermal Performance Indicator
From a thermal point of view, all the structural and non-structural elements are subjected to a heat flux that causes thermal dispersions. The main parameter to evaluate these dispersions is the thermal transmittance U, defined as the rate of transfer of heat through a unitary surface with a difference of temperature of 1 • C. For its calculation, the UNI EN ISO 6946:2008, [32] gives the following expression: (3) R si is the thermal resistance of the internal surface, m 2 K/W ; The term λ i represents the thermal conductivity of the element and measures the tendency of transmitting heat, which can be found in UNI 10351:1994 [33].
At the same time, the variation of thermal transmittance ∆U was adopted as a reference parameter to evaluate the thermal variation due to the intervention. Therefore ∆U has been expressed by the ratio between the performance variation (the difference between the value after (U 1 ) and before the integrated intervention (U 0 )) to the initial value U 0 : Each intervention entails a reduction of the thermal transmittance U and the parameter ∆U allows us to measure its entity.
The mesoscale approach, based on considering an entire façade with openings, requires the evaluation of the thermal transmittance considering these elements.
The variables that influence the calculation of the heat transmission of a transparent layer are glass and support typologies including, for instance, spacers. For that calculation, one can combine the heat transmission weighting them with respect to the area and adding to this contribution to the effect of the thermal bridging determined in the interface glass-support. The following expression, taken from the UNI EN ISO 10077-1:2007 standard [34], defines the thermal transmittance of the window: where: This expression is valid for new buildings. For existing buildings, one can refer to the tables reported in Attachment F of UNI EN ISO 10077-1:2007. In particular, the following procedure can be adopted: (a) select the thermal characteristics of the frame support U f ; (b) select the thermal characteristics of the glass U g ; (c) cross the values of U f and U g by selecting the percentage of the support with respect to the entire opening and find the value of U w of the opening with the chosen characteristics.
Once that U w is defined, a unique value of thermal transmittance U has to be defined, representative of the whole façade and capable to take into account the thermal transmittances of the walls U m and of the windows U w . To define such a parameter, an expression analogous to Equation (5), considering that the masonry portions take the place of the support frame of the windows, is: in which: The linear thermal transmittance, Ψ, is the heat flux in steady-state divided by the thermal bridging length and by the difference of temperature between the elements located on each side of the thermal bridge. According to UNI EN ISO 10211:2008 [35], the thermal bridging is a part of the building where the thermal resistance is not uniform but significantly varies due to many situations, e.g., a variation of thickness, connections between wall and floors or wall and roof, etc.
The thermal bridging implies an increase of the specific heat transfer with respect to that occurring in the surrounding materials, creating a path of minimum resistance for heat transfer.
The standard UNI EN ISO 14683:2008 [36] indicates how to determine the linear transmittance of thermal bridging. Some reference values are given in Attachment A of the mentioned standard. Some common types of linear thermal bridging are displayed in Table 1 and Figure 2. The thermal bridging implies an increase of the specific heat transfer with respect to that occurring in the surrounding materials, creating a path of minimum resistance for heat transfer.
The standard UNI EN ISO 14683:2008 [36] indicates how to determine the linear transmittance of thermal bridging. Some reference values are given in Attachment A of the mentioned standard. Some common types of linear thermal bridging are displayed in Table 1 and Figure 2.  In such a way, the thermal transmittance can be easily calculated for the as-built and for the retrofitted state to obtain the variation needed for the iso-cost curves.

Structural Performance Indicators
With the purpose of determining the relative variation of the wall structural resistance, the variation of the maximum base shear Δ and the variation of the ductility Δ are selected as reference indicators.
The relative variation of the structural parameter Δ is assessed by the ratio of the performance variation between its value after ( ) and before the retrofitting ( ) to the initial value ( ): Similarly, the variation of the ductility Δ reads: After each intervention, the structural performance of the façade improves, therefore, the parameters Δ and Δ measure the structural improvement in terms of base shear and ductility, respectively. In such a way, the thermal transmittance can be easily calculated for the as-built and for the retrofitted state to obtain the variation needed for the iso-cost curves.

Structural Performance Indicators
With the purpose of determining the relative variation of the wall structural resistance, the variation of the maximum base shear ∆V and the variation of the ductility ∆µ are selected as reference indicators.
The relative variation of the structural parameter ∆V is assessed by the ratio of the performance variation between its value after (V 1 ) and before the retrofitting (V 0 ) to the initial value (V 0 ): Buildings 2020, 10, 143 Similarly, the variation of the ductility ∆µ reads: After each intervention, the structural performance of the façade improves, therefore, the parameters ∆V and ∆µ measure the structural improvement in terms of base shear and ductility, respectively.
These two parameters are calculated through non-linear static (pushover) analyses on the entire façade. The non-linear analysis allows then obtaining the capacity curves from which it is possible to identify the parameters' values necessary for the study. The maximum base shear V represents the maximum value of the base shear shown in the analysis' steps, instead, the ductility µ is calculated by the ratio between the maximum displacement and the displacement at the onset of yielding (identified on the equivalent bilinear curve by which the capacity curve is approximated).

Application of the Integrated Approach at Mesoscale Level
In this section, the integrated approach applies at the mesoscale level to a case study consisting of a façade of an existing stone masonry building.

Case of Study
The building under examination is the seat of a Municipality located in Northern Italy ( Figure 3). The two-storey building is made of stone masonry arranged in good texture, with a structural thickness of 60 cm at the ground floor and 52 cm on the first floor. The numerical analysis is performed for the case under study with 3D MACRO software [37], which models masonry buildings through a macro-element approach. This method can be very useful for existing buildings, especially historic buildings for which the macro-element approach reveals to be reliable if compared with common finite element models [38]. The assumptions of nonlinear material according to the mechanical parameters listed in Table 2

Application of the Integrated Approach at Mesoscale Level
In this section, the integrated approach applies at the mesoscale level to a case study consisting of a façade of an existing stone masonry building.

Case of Study
The building under examination is the seat of a Municipality located in Northern Italy ( Figure  3). The two-storey building is made of stone masonry arranged in good texture, with a structural thickness of 60 cm at the ground floor and 52 cm on the first floor. The numerical analysis is performed for the case under study with 3D MACRO software [37], which models masonry buildings through a macro-element approach. This method can be very useful for existing buildings, especially historic buildings for which the macro-element approach reveals to be reliable if compared with common finite element models [38]. The assumptions of nonlinear material according to the mechanical parameters listed in Table 2-and estimated from [39]-are made. The considered façade, displayed in Figure 3a, is quite regular with five openings (one door and four windows) for each storey. The macro-element model with geometrical dimensions is shown in Figure 3b.

Assumptions on the Integrated Interventions
As for masonry, the assumption is that the properties of blocks and mortar joints of each panel are homogenized.
The considered types of interventions are similar to those assumed in [28] and are indicated in Figure 4.  The considered façade, displayed in Figure 3a, is quite regular with five openings (one door and four windows) for each storey. The macro-element model with geometrical dimensions is shown in Figure 3b.

Assumptions on the Integrated Interventions
As for masonry, the assumption is that the properties of blocks and mortar joints of each panel are homogenized.
The considered types of interventions are similar to those assumed in [28] and are indicated in Figure 4. The two latter techniques have a significant improvement in mechanical behavior and a low thermal conductivity. Nevertheless, they imply a reduction of the wall ductility, which has to be considered as detrimental in the seismic zone. Solution (f) is an application of a GFRP net on both sides of the wall with a mesh of approximate dimensions 100 mm × 100 mm. It has a tensile strength equal to 3500 N and axial stiffness of 230 kN/m. This solution does not increase the thermal resistance due to the absence of an insulating layer, but permits to improve the structural response.

The Procedure of Estimation of the Varied Mechanical and Thermal Properties
The variation of mechanical properties in the presence of integrated interventions can be estimated, with the purpose of large scale analysis, through the indications reported in [39]. For existing masonry buildings, corrective coefficients are given to modify the strength and elastic modulus values. For instance, when transverse connections are considered, the structural improvement is evaluated by incrementing the masonry mechanical properties by 30%. Analogously, the increase of mechanical properties of the stone masonry when ferro-cement is added is 50%. As for the fiber-reinforced polymer-based interventions, an elasto-plastic constitutive law is considered, with an assumed strain limit. Beyond the strength value, the maximum delamination force is considered for the fiber composite system.
An FRP-based solution increases both flexural and shear strength: the increment of tensile Solution (a) consists in a polystyrene insulating panel with a thermal conductivity λ = 0.04 W/mK. This is a typical intervention for which the structural strength is not increased, but the energy efficiency is. The same insulating panel is adopted in solution (b), with the addition of transverse connectors to get an additional improvement in the structural sense. Solution (c) consists in the classical ferro-cement made of two layers of concrete (λ = 1.60 W/mK) on the two sides of the masonry wall, applied over an armature of metal mesh and closely spaced steel rods. This type of intervention allows improving the mechanical properties of masonry, but does not offer good thermal insulation. In solution (d) Carbon Fiber Reinforced Polymers (CFRP) strips are installed on both sides of the masonry panel. Their thickness is 0.133 mm and the thermal conductivity is λ = 0.08 W/mK. The tensile mechanical strength of CFRP is f CFRP = 4.83 GPa and the normal elastic modulus is E CFRP = 252 GPa. Analogously, solution (e) consists in Glass Fiber Reinforced Polymers (GFRP) of thickness 0.48 mm and thermal conductivity of λ = 0.04 W/mK. GFRP tensile strength is f GFRP = 2.56 GPa and the elastic modulus E CFRP = 80.7 GPa. The width of the strips and their disposition vary depending on the extension of the intervention. Only horizontal and vertical strips are considered. The two latter techniques have a significant improvement in mechanical behavior and a low thermal conductivity. Nevertheless, they imply a reduction of the wall ductility, which has to be considered as detrimental in the seismic zone. Solution (f) is an application of a GFRP net on both sides of the wall with a mesh of approximate dimensions 100 mm × 100 mm. It has a tensile strength equal to 3500 N and axial stiffness of 230 kN/m. This solution does not increase the thermal resistance due to the absence of an insulating layer, but permits to improve the structural response.

The Procedure of Estimation of the Varied Mechanical and Thermal Properties
The variation of mechanical properties in the presence of integrated interventions can be estimated, with the purpose of large scale analysis, through the indications reported in [39]. For existing masonry buildings, corrective coefficients are given to modify the strength and elastic modulus values. For instance, when transverse connections are considered, the structural improvement is evaluated by incrementing the masonry mechanical properties by 30%. Analogously, the increase of mechanical properties of the stone masonry when ferro-cement is added is 50%. As for the fiber-reinforced polymer-based interventions, an elasto-plastic constitutive law is considered, with an assumed strain limit. Beyond the strength value, the maximum delamination force is considered for the fiber composite system.
An FRP-based solution increases both flexural and shear strength: the increment of tensile strength is evaluated by considering the ultimate strength of the reinforcement f f R given by the minimum between the design strength of the fibers and the delamination strength. The increases of strength in vertical (∆σ t,v ) and horizontal (∆σ t,v ) direction respectively are [40]: where: s is the wall thickness; A h , A v are the reinforcement area placed in horizontal and vertical direction per unit of length; p h , p v are the spacings between horizontal or vertical strips.
Analogously, it is possible to define the increment of shear strength: where V R is the shear strength of reinforced masonry that is determined as the sum of the strength of the unreinforced wall (V Rm ) and the shear strength offered by the fibers V R f , namely: in which: The thermal transmittance U, and its variation ∆U needed for the iso-cost curves, is calculated for the façade with the indications reported in UNI EN ISO 6946:2008, as stated in Section 2.2 In particular, the façade under examination has eight windows and two doors with aluminum frame support with double glazed glass 4-9-4. The transmittance of the frame support is U f = 3.0 W/m 2 K, and of the glass, U g = 3.1 W/m 2 K. By crossing these values, the window has a value of thermal transmittance of U w = 3.3 W/m 2 K. The wood door is assumed to have U d = 1.5 W/m 2 K (wood thermal conductivity of λ = 0.22 W/mK). As for the linear thermal transmittance, it has been assumed Ψ = 0.60 W/mK. Finally, the thermal transmittance of the entire façade in the as-built configuration is U = 2.937 W/m 2 K. In case of the retrofitting interventions described in Section 3.2., the values of the thermal transmittance are analogously calculated (Tables 3 and 4). All the interventions cause a reduction of the thermal transmission; the best improvement is obtained with the polystyrene panel. The corresponding base shear capacity and ductility values are reported in Tables 5-8.  Table 3.
Thermal transmittance U of the entire façade with the different integrated interventions-economic analysis.

Economic Cost Analysis
Analogously to the study performed in [28], an economical relationship has been found between ∆V (or ∆µ), that is the seismic performance indicator, and ∆U, the energetic performance indicator. Six different solutions have been considered (Figure 2), with costs ranging from 100 €/m 2 to 350 €/m 2 . For the sake of comparison, the same unitary costs (supply and manpower) as those considered in [28] have been used for each individual intervention (Section 3.2.) as indicated in Table 9. For each intervention, the layer thickness has been tuned to meet a target cost between 100 €/m 2 to 350 €/m 2 . It is then possible to calculate by means of the expressions reported in Sections 2.2 and 2.3. the capacity curves ∆U − ∆V and ∆U − ∆µ, each one representing the performance offered by each intervention among the six selected. Each curve is obtained by curve fitting of the points, each one representative of a single integrated intervention. In particular, the points in the graph (displayed in Figure 5a) have been interpolated with a proper curve fitting through hyperbolic regression curves: Buildings 2020, 10, x FOR PEER REVIEW 12 of 18 (a) (b) Figure 5. Iso-cost ∆ − ∆ (a) and ∆U−∆ (b) regression curves. Figure 5 shows that the maximum structural improvement is attained with FRP stripes (solutions d and e) in terms of increase of shear base strength ∆ . Nevertheless, these interventions significantly reduce the increase in ductility. Table 11 reports the corresponding regression coefficients for the curves ∆ − ∆ . Table 11. Economic cost regression coefficients for ∆ − ∆ regression curves.

Environmental Cost Analysis
The same analysis performed in economic terms is proposed here from an environmental point of view. The aim of this analysis is to limit the CO2eq emissions that are normally produced when retrofitting interventions are made. The CO2eq emission is one of the main indicators chosen to measure building sustainability.
As any other commercial product, each building product implies a certain level of emissions: extraction of raw materials, manufacturing process, transport of materials from the source to the destination, the emissions during the construction, use, maintenance works and demolition, in a perspective of life cycle assessment. Provided that the aim of the paper is the choice of the optimal consolidation technique, only the CO2eq emissions related to the phase of production of the materials for intervention are considered. Although this analysis can be extended to the entire life of the building, this paper neglects that aspect, which will be considered in future works. The assumed parameter of CO2eq emission varies from 5 kgCO2eq/m 3 to 35 kgCO2eq/m 3 . The specific environmental costs for each material are reported in Table 12. The parameters α 0 and α 1 are obtained by means of a minimum square criterion and are listed in Table 10. Hyperbolic type curves provided the best estimation by using a reduced number of parameters to describe this relationship. It is worthy to notice that not always the fitting capacity curves are of hyperbolic shape. Indeed, for this set of retrofitting strategies, this is the case, but for other retrofitting techniques the response will generally be different. If, for example, integrated approaches resulting in similar variations of thermal transmittance but significantly increase structural strength, the fitting capacity curve would be flatter and tending to a horizontal line. The analogous approach has been adopted for the ∆U−∆µ regression curves (Figure 5b).  Figure 5 shows that the maximum structural improvement is attained with FRP stripes (solutions d and e) in terms of increase of shear base strength ∆V. Nevertheless, these interventions significantly reduce the increase in ductility. Table 11 reports the corresponding regression coefficients for the curves ∆U − ∆µ.

Environmental Cost Analysis
The same analysis performed in economic terms is proposed here from an environmental point of view. The aim of this analysis is to limit the CO 2eq emissions that are normally produced when retrofitting interventions are made. The CO 2eq emission is one of the main indicators chosen to measure building sustainability.
As any other commercial product, each building product implies a certain level of emissions: extraction of raw materials, manufacturing process, transport of materials from the source to the destination, the emissions during the construction, use, maintenance works and demolition, in a perspective of life cycle assessment. Provided that the aim of the paper is the choice of the optimal consolidation technique, only the CO 2eq emissions related to the phase of production of the materials for intervention are considered. Although this analysis can be extended to the entire life of the building, this paper neglects that aspect, which will be considered in future works. The assumed parameter of CO 2eq emission varies from 5 kgCO 2eq /m 3 to 35 kgCO 2eq /m 3 . The specific environmental costs for each material are reported in Table 12. The regression parameters α 0 and α 1 are reported in Table 13 and in Table 14, whilst Figure 6 shows the environmental regression curves for ∆U−∆µ (∆U−∆V).  The regression parameters α0 and α1 are reported in Table 13 and in Table 14, whilst Figure 6 shows the environmental regression curves for ∆U−∆ (∆U−∆ ).  It should be noticed that, for the same motivation mentioned above, solutions (a) and (b) are the best choices to obtain a strong improvement of the thermal performance of the wall, whereas solutions (d) and (e) guarantee an intermediate increment of the base shear, but, if the quantity of the FRP stripes increases, there is a reduction of ductility which is not optimal from a seismic point of view. It should be noticed that, for the same motivation mentioned above, solutions (a) and (b) are the best choices to obtain a strong improvement of the thermal performance of the wall, whereas solutions (d) and (e) guarantee an intermediate increment of the base shear, but, if the quantity of the FRP stripes increases, there is a reduction of ductility which is not optimal from a seismic point of view.

Definition of Demand Curves and Discussion of Results
Steps (5) and (6) of the procedure described in Section 2.1 lead to graphs reporting demand and capacity curves, which synthetically describe the structural, economic and environmental impact of the interventions adopted. Some sites are chosen to compare the demand of the site construction of the case study (Rocchetta di Vara, Italy) to that of other significant sites in Italy (Torino, L'Aquila, Catania, Cagliari, Italy). The demand lines are obtained, as already discussed, according to the equation: with symbols defined in Section 2.1. As visible in Figures 7 and 8, the demand line of Rocchetta di Vara stands between those chosen as reference locations. The hyperbolic graphs in the same Figures represent the trend of capacity obtained by the retrofitting techniques, in terms of seismic and thermal improvements. The graphs in Figure 7 are plotted in terms of economic costs (investment in €), and in Figure 8, in terms of environmental costs (CO 2 ). The optimal point for each location is given by the intersection of the lines with the hyperboles. The decision-makers can act through two aspects: which hyperbole should be considered (i.e., which level of economic or environmental cost) and which "tuning parameter" α has to be assumed (i.e., the slope of the line graphs) to give more importance to the thermal or to the seismic aspect, provided that the location influences the inclination of the line through the coefficients of Equation (1

Conclusions
This paper presented a mesoscale approach for masonry buildings to optimize retrofitting interventions capable of improving the structural and thermal performances of an existing building. These integrated interventions can be read in terms of seismic indicators, such as variation of base shear strength and variation of ductility, and in terms of a thermal indicator, such as variation of thermal transmittance. The performance indicators are calculated for external façades of masonry buildings (mesoscale level, different from the macroscale involving the whole building), considering their openings (windows, doors). Economic and environmental iso-cost curves have been obtained from cost analysis to measure the level of financial commitment and sustainable impact on the environment associated with each type of intervention. The graphs reporting demand and capacity curves were revealed to be effective and of immediately readable since, for a fixed maximum budget

Conclusions
This paper presented a mesoscale approach for masonry buildings to optimize retrofitting interventions capable of improving the structural and thermal performances of an existing building. These integrated interventions can be read in terms of seismic indicators, such as variation of base shear strength and variation of ductility, and in terms of a thermal indicator, such as variation of thermal transmittance. The performance indicators are calculated for external façades of masonry buildings (mesoscale level, different from the macroscale involving the whole building), considering their openings (windows, doors). Economic and environmental iso-cost curves have been obtained from cost analysis to measure the level of financial commitment and sustainable impact on the environment associated with each type of intervention. The graphs reporting demand and capacity curves were revealed to be effective and of immediately readable since, for a fixed maximum budget As for the demand curves, the tuning parameter α has been here set equal to 1. Table 15 lists the parameters needed for the definition of the demand curves. These graphs are of immediate use: for a fixed maximum budget that a public administration or a private owner has, it is possible to estimate the level of performance that can be gained in terms of thermal and structural behavior. Therefore, for instance with a budget of 100 €/m 2 one should expect in Rocchetta di Vara an increase of base shear strength by 7%, of ductility by 10% and of thermal transmittance by 8%, whereas with a budget of 300 €/m 2 the expected increase of base shear and thermal transmittance are respectively by 14% and 17%. It should be noticed that the information gathered by this procedure is able to give an overall look at the improvement obtained by integrated approaches. An alternative way to proceed is to define the capacity curves by keeping constant the type of intervention and tuning its thickness, or coupling two interventions, one with more impact on the structural performance, and the other more effective from a thermal point of view. This procedure, which will be investigated in the future, will allow us to refine the decisional practice about the intervention to adopt.
It can be observed that there is not a monotonic trend of improvement of the performance indicators by increasing the environmental costs: indeed, for instance, the curve corresponding to 350 €/m 2 stands between that obtained for 100 €/m 2 and that for 300 €/m 2 (Figure 8). This aspect is due to the fact that for the FRP interventions, beyond a threshold value, the increase of the cost (and therefore of the entity of the intervention) does not imply an increment of ductility, but instead a decrement of it. In other words, a greater thickness of FRP strips worsens the structural performance. This influence is stronger with respect to the other interventions and therefore the fitted curves do not follow a monotonic trend. A way to avoid this apparently counter-intuitive response is, as mentioned above, to plot the fitted curves for interventions with similar and monotonic effects.

Conclusions
This paper presented a mesoscale approach for masonry buildings to optimize retrofitting interventions capable of improving the structural and thermal performances of an existing building. These integrated interventions can be read in terms of seismic indicators, such as variation of base shear strength and variation of ductility, and in terms of a thermal indicator, such as variation of thermal transmittance. The performance indicators are calculated for external façades of masonry buildings (mesoscale level, different from the macroscale involving the whole building), considering their openings (windows, doors). Economic and environmental iso-cost curves have been obtained from cost analysis to measure the level of financial commitment and sustainable impact on the environment associated with each type of intervention. The graphs reporting demand and capacity curves were revealed to be effective and of immediately readable since, for a fixed maximum budget that a public administration or a private owner has at their disposal, it is possible to find the level of improvement that can be gained in terms of thermal and structural behavior. An alternative way to proceed is to define the capacity curves by keeping constant the type of intervention and tuning its thickness, or combining two interventions. This procedure, which will be investigated in the future, will allow us to refine the decisional practice about the intervention to adopt. Finally, another future development will extend the procedure to a macroscale level, by considering seismic and thermal performance indicators able to properly take into account the performance of the whole building.
Author Contributions: L.G. and M.S. conceived the methodology and its application to real cases; L.G. analyzed the case study and wrote the paper. S.P. made calculations and produced tables and figures. All authors have read and agreed to the published version of the manuscript.