Temperature variations in masonry façades can induce expansion and contraction movements. When these movements are restrained, cracking and material degradation may occur, especially in older buildings lacking movement or expansion joints. Such temperature variations arise from factors as solar radiation, shading, material color, reflectivity, and environmental conditions. This study investigates the magnitude and spatial distribution of surface temperature variations (ΔT) on exterior masonry wall surfaces using outdoor infrared (IR) thermography. A better understanding of the magnitude and distribution of ΔT is essential for accurate damage assessment and for improving the attributability of observed damage to temperature effects rather than to other causes. Field data were collected in Delft, the Netherlands. Thermal images were captured with an IR camera to identify temperature differences across various points on exterior wall surfaces under direct solar radiation and varying shading conditions. The acquired imagery was analyzed using temperature histograms and profiles to quantify thermal gradients over the surface area of the façades. Results revealed significant spatial temperature variations, with measured ΔT values reaching up to 13 °C between the warmest and coolest zones on individual façades. Even where façades showed no pronounced surface gradients, temperature differences of up to 6 °C occurred between different, contiguous exterior walls of the same building. The study demonstrates that outdoor thermography, combined with targeted image processing, effectively identifies thermal gradients on masonry façades. These gradients reflect uneven thermal responses under real environmental conditions, which can accelerate moisture-related damage, cracking, and material fatigue. The findings emphasize the need to account for surface temperature heterogeneity in damage assessment of existing structures.

This paper investigates settlement-induced damages in unreinforced masonry (URM) walls using a high-fidelity block-based numerical modeling approach. The research aims to address gaps in the understanding of settlement effects on URM walls with flanges, particularly with respect to their seismic out-of-plane (OOP) behavior. A parametric study is conducted on four wall specimens with varying geometries, boundary conditions, and settlement scenarios, including symmetric and asymmetric patterns. The numerical models are developed via a high-fidelity block-based finite element method that simulates masonry using expanded blocks connected by zero-thickness joints, allowing for detailed analysis of cracking patterns and damage mechanisms. Different damage states, from no visible cracks to near-collapse conditions, are identified in the response of the walls and are used as initial conditions for subsequent monotonic static pushover OOP loading. The results highlight the significant influence of settlement-induced pre-damages on the OOP response of URM walls, with varying degrees of impact observed across different specimen configurations. The findings underscore the importance of considering even “light” settlement-induced pre-damages when assessing the seismic performance of URM structures, particularly in subsidence-prone regions. Under symmetric hogging, such pre-damage level can reduce OOP stiffness and peak strength by up to 41% and 20%, respectively. This study lays the groundwork for future investigations into the seismic behavior of pre-damaged masonry structures under dynamic loading and offers valuable insights for the development of more accurate assessment and mitigation strategies for buildings subjected to settlement deformations.

This paper presents a framework for developing fragility curves for masonry buildings on strip foundations exposed to subsidence using non-linear finite element (NLFE) analyses. A 2D plane-stress model of a masonry façade is used to evaluate the probability of cracking damage resulting from settlements. The model simulates the behaviour of typical Dutch two-storey historical buildings, using an established modelling approach to represent the non-linear behaviour of the façade, the transversal walls and the strip foundation, supported by a base interface for soil-foundation interaction. Settlements are imposed at the bottom of the interface, characterizing their intensity with the angular distortion. The damage severity is objectively quantified using the scalar parameter Ψ, computed considering the number, length, and width of the cracks. Cumulative probability functions are derived from 864 numerical analyses that account for realistic variations in building and soil features, including 3 masonry materials, 2 strip foundation systems, 2 interface soil materials and the 72 possible settlement patterns. The effect of each selected variation is evaluated individually. The proposed curves reveal a probability of over 25% for cracks up to 5 mm in width when the angular distortion equal to 0.2% (or 1/500), the threshold deemed acceptable by international codes, is applied to the models. Doubling the applied angular distortion results in an approximate doubling of the probability of damage. While the proposed curves are specific to the selected geometry, the framework can be adapted to accommodate different façade geometries, enabling the development of more comprehensive fragility functions.

Temperature effects are frequently cited as the cause of light cracking in masonry façades, yet most modelling studies idealise thermal loading as uniform steps and represent restraint as fully fixed, assumptions that tend to exaggerate damage. This work evaluates whether realistic, non-uniform temperature gradients, like those produced by shading and insolation, together with soil–structure interaction as the dominant restraint mechanism, can generate cracking patterns consistent with field observations. A coupled thermo-mechanical FEM model with a homogenised masonry continuum and tensile softening is employed; the façade–foundation–soil system is modelled explicitly, and damage is quantified using a crack-based index Ψ. A parametric campaign (1200 simulations) spans two façade typologies (clay masonry on unreinforced masonry foundations; calcium-silicate on reinforced concrete strips), three layered soils, 33 geometries, and multiple vertical and two-dimensional gradient shapes. The results indicate that gradient shape is decisive: widely distributed vertical gradients trigger visible damage (Ψ≥1) at roughly half the temperature differential required by more localised gradients, with visible damage becoming likely around ΔT≈20 °C (warming) and ≈25 °C (cooling) for the most adverse shapes. Restraint stiffness governs severity: stiffer sandy profiles increase tensile stresses and cracking, whereas softer profiles accommodate thermal movement; relative to uniform, fully restrained models, crack initiation is delayed by ∼15–20 °C and cracking is less distributed. Geometric discontinuities also dominate sensitivity: larger/more openings and low vertical-masonry ratios promote earlier localisation, while overall length/height is secondary. Fragility-like curves provide thresholds useful for assessment and mitigation.

Soil heterogeneity, due to variations in the subsurface stratigraphy or properties within a layer, can trigger or amplify differential settlements that affect buildings and infrastructure and can thus lead to (increase in) damage. The state-of-the-art mainly focuses on the effect of heterogeneous properties within a layer on engineering problems. From this, it is known that the variation in properties can increase the vulnerability of a structure. However, nearly always variations in the soil lithological conditions are disregarded, while they can influence subsidence potentially even more. Lithological variations are relevant both at the scale of individual buildings as well as different scales (city, regional, country), for which often detailed soil information is not available. Thus, for a better prediction of potential building damage related to subsidence, knowledge about the scale and influence of lithological variations is needed. This paper describes an approach to quantify and investigate the influence of lithological heterogeneity at the scale of a single building. Moreover, this exploratory study evaluates the influence of lithological heterogeneity on the spatial variability of settlements, intending to upscale the approach to regional application. Two independent datasets at high resolution (site-specific) and low resolution (national level) are used to retrieve the stratigraphic conditions for the area selected for the analyses. One-, Two- and Three-dimensional numerical models, based on the collected information are used to simulate the consolidation process and settlement due to a uniform load imposed on the surface level of the study area. Additional analyses investigate the influence of loading conditions and groundwater table. The parameter “correlation length” is used to quantify the spatial variability of the soil layer thickness and then of the computed settlements. The analyses reveal that the spatial variability of the soil strata thickness matches that of the computed settlements, ranging from 2 to 10 meters. In other words, the lithological variability of the soil leads to differential settlements occurring at the scale of man-made structures such as houses, roads, and embankments. Thus, the results encourage including the contribution of lithological heterogeneity in models and predictions of differential settlement at the scale of individual structures. Moreover, the statistical properties, in terms of mean, spread and distribution shape, of the settlement computed through in-situ specific models, match with those derived at the national scale. These results are expected to support the identification of areas potentially influenced by lithological soil heterogeneity, thus showing potential for upscaling to regional or national levels.

Subsidence caused by natural or human-induced factors can occur unevenly, resulting in differential settlements. Existing unreinforced masonry (URM) buildings are susceptible to damage from differential settlements. However, the extent of the damage varies between structures, depending on factors such as the magnitude and pattern of the settlements, along with the features of the building and the properties of the underlying soil. Non-linear finite element analyses (NLFEA) are often used for studying the damage response, accounting for variability in soil and structural features. This study uses 6912 NLFEA, including 8 variations in façade geometry, 3 masonry materials, 2 soils, 2 shallow foundation systems, and 72 settlement patterns, to develop fragility curves for URM buildings undergoing subsidence. Old Dutch URM buildings with strip foundations are modelled using 2D plane-stress façade models, accounting for non-linear smeared shearing, cracking and crushing of masonry and 3D effects of transverse walls. Settlement troughs are applied at a non-linear soil-foundation interface, with angular distortion (β) progressively increasing to quantify settlement intensity and building deformation. As β increases, the NLFE models exhibit progressive cracking damage, with severity objectively assessed through the parameter Ψ considering crack width, length, and number. Then, the distortion β is used as the demand parameter to develop the fragility curves. The analysis shows that long façades are twice as likely to experience 5 mm cracks from settlement damage compared to short façades under an applied β of 2 ‰ (1/500). For this applied β, proposed as an acceptable limit for many structures in the Eurocode, half of the models exhibit cracks up to 5 mm wide. Therefore, while 1/500 may be considered safe for structural integrity, it can still lead to noticeable damage. Light damage occurs even at angular distortion values below 0.5 ‰ (1/2000), with 10 % of models showing cracks up to 1 mm wide.

In the Netherlands, subsidence due to different causes is linked to damage to the ubiquitous masonry structures. Finite element (FE) analyses can be used to assess the response of the structures subjected to settlements. This paper presents the comparison between three-dimensional FE modelling strategies to investigate the response of an unreinforced masonry building on a strip foundation. The aim is to investigate whether different modelling approaches demonstrate consistent results. The soil-structure system is modelled employing two strategies: a coupled model, in which the structure is tied to the soil volume, and an uncoupled approach that divides the soil and structure into two sub-systems. Two displacement fields, imposed at the bottom of the soil volume, idealize various shapes of the subsidence troughs, with increasing intensity measured by their distortion. Non-linear interfaces are used to simulate the soil-foundation interaction, and their stiffness values vary based on the type of model. The displacements, interface stresses and crack patterns of the selected modelling strategies are consistent. The interface types do not influence the response of the façade, whereas the shape of the settlement does play a key role. The uncoupled models exhibit, on average, slightly higher values of damage than coupled models for a given imposed distortion. The two modelling strategies require almost the same computational time and show similar convergence. Because of the limited contribution of small soil volumes in uncoupled models, the superstructure sub-system can be directly utilized to assess the response of structures undergoing vertical displacements, thereby reducing the modelling burden.

In this study, 2D and 3D modelling strategies are used to represent the behaviour of historical masonry buildings on strip foundations undergoing settlements. The application focuses on a two-story building, typical of the Dutch architectural heritage. An improved 2D modelling is presented: It includes the effect of the lateral walls to replicate the response of the detailed 3D models. The masonry strip foundation is modelled and supported by a no-tension interface, which represents the soil-foundation interaction. Two settlement configurations, hogging and sagging, are applied to the models, and their intensity is characterized using their angular distortion. The improved 2D model that includes the stiffness and weight of the lateral walls agrees in terms of displacements, stress and damage with the detailed 3D models. Conversely, the simplified 2D façade models without lateral walls exhibit different cracking, and damage from 2 to 7 times lower at an applied angular distortion of 2‰ (1/500). The improved 2D model requires less computational and modelling burden, resulting in analyses from 9 to 40 times faster than the 3D models. The results prove the importance of identifying and including the 3D effects that affect the response of structures subjected to settlements.

Historical masonry façades are sensitive to various damaging processes. A recent study, looking at the initiation and progression of cracks in masonry, in the range of 0.1 to 5 mm in width and thus corresponding to light damage [1], has allowed for the calibration of finite-element models that include a material model capable of accurately replicating this damage and which is populated with material properties corresponding to existing structures [6]. The models, which also include a soil-structure interaction boundary designed to account for the effect of the soil during earthquake vibrations [7], have been used to determine the fragility of masonry buildings via the proxy of 2D walls [2]. In the study presented herein, the finite element models are employed to replicate the geometry of (historical) masonry facades to determine their sensitivity to light damage as a consequence of the two damaging processes observed to be most common for this type of façade, namely (differential) settlements and (earthquake) vibrations [3]. The masonry façades were first pre-damaged via settlement distortions which generate just-visible cracks in the order of 0.1 mm to 1 mm in width. Then, an acceleration time history corresponding to two different Dutch earthquake events and two recordings of traffic-induced building vibrations [9] were separately applied at the base of the models. In this manner, the effect of existing damage could be assessed in regards to the aggravation generated by vibrations. The settlement part of the study revealed that long façades were more vulnerable to applied soil distortions, for instance. Then, subsequent vibrations further increased damage for intensities measured with a peak ground velocity (PGV) larger than 2 mm/s while the control set of virgin or uncracked façades remained undamaged at this PGV. At 32 mm/s, many pre-damaged façades also exceeded the light damage range. At equal PGV, the traffic vibrations, with a larger number of effective cycles, resulted in increased damage aggravation in comparison to the earthquake recordings.

This study aims to investigate the damage response of unreinforced masonry (URM) façades resting on strip foundations and subjected to ground settlements via numerical models. The models depict the non-linear constitutive behaviour of both the masonry, via smeared cracking, and of the soil-foundation interaction, via nonlinear interface elements. The influence of building features, such as the masonry material, the length over height (L/H) ratio of the geometry, the wall thickness, the number and size of openings and different types of strip foundations (i.e. reinforced concrete and unreinforced) is examined. A sensitivity study additionally investigates the influence of the interface stiffness and its constitutive model. A Gaussian curve is used to replicate the shape of the ground settlements; These simulate the loss of support underneath the foundation due to urban subsidence. Eight settlement shapes are applied in the FE models, including both symmetric and asymmetric profiles, while the angular distortion is used to measure their intensity. A new aspect is that the extent of the induced damage to the façade is assessed objectively using a damage parameter that represents the number, length and width of cracks in a single scalar value. The method distinguishes between the applied settlement profile at the bottom of the interface and the retrieved settlement profile measured on the façade. The analyses indicate that for a value of the angular distortion equal to 2 ‰ (or 1/500), computed from the resulting deformations of the façades, 60% of the models exhibit serviceability damage associated with cracks of about 5 mm width. Accordingly, the limit values available in the literature are observed to be too optimistic and not conservative in relation to the analyses presented in this study. A key outcome is that facades with an L/H smaller or equal to 1 do not exhibit cracks wider than 1 mm. Façades on reinforced concrete foundations were observed to be less susceptible to settlement damage, compared to unreinforced ones.

This study presents a novel framework in which numerical modelling contributes to the performance of district-scale, subsidence-induced damage assessment in cities where ground settlements affect entire quarters. Therein, the implementation of expeditious procedures offers geotechnical engineers the possibility of contributing beyond the typical site scale. For this purpose, several “typified” hydro-geomechanical-loading (HGL) models, which represent (simplified) scenarios of masonry buildings undergoing settlements, were set up to account for different predisposing or triggering factors (i.e., soil heterogeneity, loading conditions, and groundwater variations) of settlement occurrence in built-up areas. These models exploit multi-source, wide-area input datasets encompassing the hydro-mechanical properties of geomaterials, in situ investigations and measurements (e.g., groundwater levels in wells), and innovative remote sensing data (i.e., DInSAR techniques). With reference to a district in Rotterdam City (the Netherlands), which was built on “soft soils”, the numerical simulations of different scenarios (i) provide an overview of the comparative role of predisposing or triggering factors on settlement occurrence and (ii) allow assessments of the expected induced damage to masonry buildings over 30 years with the exploitation of fragility curves. Considering the widespread diffusion of such geohazards, the proposed approach could help prioritise (rather expensive) maintenance work to the built heritage within sustainable strategies for subsidence risk mitigation.