Modelling of settlement induced building damage

This thesis focuses on the modelling of settlement induced damage to masonry buildings. In densely populated areas, the need for new space is nowadays producing a rapid increment of underground excavations. Due to the construction of new metro lines, tunnelling activity in urban areas is growing. One of the consequences is a greater attention to the risk of damage on existing structures. Thus, the assessment of potential damage of surface buildings has become an essential stage in the excavation projects in urban areas (Chapter 1). The current damage risk assessment procedure is based on strong simplifications, which not always lead to conservative results. Object of this thesis is the development of an improved damage classification system, which takes into account the parameters influencing the structural response to settlement, like the non-linear behaviour of masonry and the soil-structure interaction. The methodology used in this research is based on experimental and numerical modelling. The design and execution of an experimental benchmark test representative of the problem allows to identify the principal factors and mechanisms involved. The numerical simulations enable to generalize the results to a broader range of physical scenarios. The methodological choice is based on a critical review of the currently available procedures for the assessment of settlement-induced building damage (Chapter 2). A new experimental test on a 1/10th masonry façade with a rubber base interface is specifically designed to investigate the effect of soil-structure interaction on the tunnelling-induced damage (Chapter 3). The experimental results are used to validate a 2D semi-coupled finite element model for the simulation of the structural response (Chapter 4). The numerical approach, which includes a continuum cracking model for the masonry and a non-linear interface to simulate the soil-structure interaction, is then used to perform a sensitivity study on the effect of openings, material properties, initial damage, initial conditions, normal and shear behaviour of the base interface and applied settlement profile (Chapter 5). The results assess quantitatively the major role played by the normal stiffness of the soil-structure interaction and by the material parameters defining the quasi-brittle masonry behaviour. The limitation of the 2D modelling approach in simulating the progressive 3D displacement field induced by the excavation and the consequent torsional response of the building are overcome by the development of a 3D coupled model of building, foundation, soil and tunnel (Chapter 6). Following the same method applied to the 2D semi-coupled approach, the 3D model is validated through comparison with the monitoring data of a literature case study. The model is then used to carry out a series of parametric analyses on geometrical factors: the aspect ratio of horizontal building dimensions with respect to the tunnel axis direction, the presence of adjacent structures and the position and alignment of the building with respect to the excavation (Chapter 7). The results show the governing effect of the 3D building response, proving the relevance of 3D modelling. Finally, the results from the 2D and 3D parametric analyses are used to set the framework of an overall damage model which correlates the analysed structural features with the risk for the building of being damaged by a certain settlement (Chapter 8). This research therefore provides an increased experimental and numerical understanding of the building response to excavation-induced settlements, and sets the basis for an operational tool for the risk assessment of structural damage (Chapter 9).

Computational modelling of building damage due to tunnelling

Underground constructions in soft ground may lead to settlement damage to existing buildings. In The Netherlands the situation is particularly complex, because of the combination of soft soil, fragile pile foundations and brittle, unreinforced masonry façades. The tunnelling design process in urban areas requires a reliable risk damage assessment. In the engineering practice the current preliminary damage assessment is based on the limiting tensile strain method (LTSM). Essentially this is an uncoupled analysis, in which the building is modelled as an elastic beam subject to imposed Greenfield settlements and the induced tensile strains are compared with a limit value for the material. The soil-structure interaction is included only as a ratio between the soil and the building stiffness. In this paper, a coupled approach is evaluated. The soil-structure interaction in terms of normal and shear behaviour is represented by interface elements and a cracking model for masonry is included. This project aims to improve the existing damage classification system for masonry buildings subjected to tunnel-induced settlement, in order to evaluate the necessity of strengthening techniques or mitigation measures.

Raninfall and Slope Characteristics Affecting Soil Erosion on Hillsopes

A soil erosion model including interrill erosion and rill erosion was developed for hillslope. The rainfall and slope characteristics affecting soil erosion on hillslopes was analyzed by simulation. The results show that the slope length and gradient, rainfall intensity and process have varying influence on runoff and soil erosion. The unit discharge of runoff and erosion rate increases with increases in the rainfall intensity and the slope length. The same precipitation but with different rainfall processes may cause different runoff and erosion results. The effect of the slope gradient on runoff and soil erosion can be both positive and negative. There exists a critical slope gradient for runoff and soil erosion.

Estimation of wind erosion rates by using Cs-137 tracing technique: A case study in Tariat-Xilin Gol transect, Mongolian Plateau

Wind erosion is one of the major environmental problems in semi-arid and arid regions. Here we established the Tariat-Xilin Gol transect from northwest to southeast across the Mongolian Plateau, and selected seven sampling sites along the transect. We then estimated the soil wind erosion rates by using the Cs-137 tracing technique and examined their spatial dynamics. Our results showed that the Cs-137 inventories of sampling sites ranged from 265.63 +/- 44.91 to 1279.54 +/- 166.53 Bq.m(-2), and the wind erosion rates varied from 64.58 to 419.63 t.km(-2).a(-1) accordingly. In the Mongolia section of the transect (from Tariat to Sainshand), the wind erosion rate increased gradually with vegetation type and climatic regimes; the wind erosion process was controlled by physical factors such as annual precipitation and vegetation coverage, etc., and the impact of human activities was negligible. While in the China section of the transect (Inner Mongolia), the wind erosion rates of Xilin Hot and Zhengxiangbai Banner were thrice as much as those of Bayannur of Mongolia, although these three sites were all dominated by typical steppe. Besides the physical factors, higher population density and livestock carrying level should be responsible for the higher wind erosion rates in these two regions of Inner Mongolia.

Modelling of industrial multi-physics processes—a key role for computational mechanics.

In the analysis of industrial processes, there is an increasing emphasis on systems governed by interacting continuum phenomena. Mathematical models of such multi-physics processes can only be achieved for practical simulations through computational solution procedures—computational mechanics. Examples of such multi-physics systems in the context of metals processing are used to explore some of the key issues. Finite-volume methods on unstructured meshes are proposed as a means to achieve efficient rapid solutions to such systems. Issues associated with the software design, the exploitation of high performance computers, and the concept of the virtual computational-mechanics modelling laboratory are also addressed in this context.

Computer modelling of human behaviour in aircraft fire accidents

The mathematical simulation of the evacuation process has a wide and largely untapped scope of application within the aircraft industry. The function of the mathematical model is to provide insight into complex behaviour by allowing designers, legislators, and investigators to ask ‘what if’ questions. Such a model, EXODUS, is currently under development, and this paper describes its evolution and potential applications. EXODUS is an egress model designed to simulate the evacuation of large numbers of individuals from an enclosure, such as an aircraft. The model tracks the trajectory of each individual as they make their way out of the enclosure or are overcome by fire hazards, such as heat and toxic gases. The software is expert system-based, the progressive motion and behaviour of each individual being determined by a set of heuristics or rules. EXODUS comprises five core interacting components: (i) the Movement Submodel — controls the physical movement of individual passengers from their current position to the most suitable neighbouring location; (ii) the Behaviour Submodel — determines an individual's response to the current prevailing situation; (iii) the Passenger Submodel — describes an individual as a collection of 22 defining attributes and variables; (iv) the Hazard Submodel — controls the atmospheric and physical environment; and (v) the Toxicity Submodel — determines the effects on an individual exposed to the fire products, heat, and narcotic gases through the Fractional Effective Dose calculations. These components are briefly described and their capabilities and limitations are demonstrated through comparison with experimental data and several hypothetical evacuation scenarios.

Temporal modelling of short-term rainfall using Cox processes

We study two marked point process models based on the Cox process. These models are used to describe the probabilistic structure of the rainfall intensity process. Mathematical formulation of the models is described and some second-moment characteristics of the rainfall depth, and aggregated processes are considered. The derived second-order properties of the accumulated rainfall amounts at different levels of aggregation are used in order to examine the model fit. A brief data analysis is presented. Copyright © 1998 John Wiley & Sons, Ltd.

Finite volume methods applied to the computational modelling of welding phenomena

This paper presents the computational modelling of welding phenomena within a versatile numerical framework. The framework embraces models from both the fields of computational fluid dynamics (CFD) and computational solid mechanics (CSM). With regard to the CFD modelling of the weld pool fluid dynamics, heat transfer and phase change, cell-centred finite volume (FV) methods are employed. Additionally, novel vertex-based FV methods are employed with regard to the elasto-plastic deformation associated with the CSM. The FV methods are included within an integrated modelling framework, PHYSICA, which can be readily applied to unstructured meshes. The modelling techniques are validated against a variety of reference solutions.

A model to predict the life of pneumatic conveyor bends

A new approach to the prediction of bend lifetime in pneumatic conveyors, subject to erosive wear is described. Mathematical modelling is exploited. Commercial Computational Fluid Dynamics (CFD) software is used for the prediction of air flow and particle tracks, and custom code for the modelling of bend erosion and lifetime prediction. The custom code uses a toroidal geometry, and employs a range of empirical data rather than trying to fit classical erosion models to a particular circumstance. The data used was obtained relatively quickly and easily from a gas-blast erosion tester. A full-scale pneumatic conveying rig was used to validate a sample of the bend lifetime predictions, and the results suggest accuracy of within ±65%, using calibration methods. Finally, the work is distilled into user-friendly interactive software that will make erosion lifetime predictions for a wide range of bends under varying conveying conditions. This could be a valuable tool for the pneumatic conveyor design or maintenance engineer.

Computational modelling of particle degradation in dilute phase pneumatic conveyors

The aim of this paper is to develop a mathematical model with the ability to predict particle degradation during dilute phase pneumatic conveying. A numerical procedure, based on a matrix representation of degradation processes, is presented to determine the particle impact degradation propensity from a small number of particle single impact tests carried out in a new designed laboratory scale degradation tester. A complete model of particle degradation during dilute phase pneumatic conveying is then described, where the calculation of degradation propensity is coupled with a flow model of the solids and gas phases in the pipeline. Numerical results are presented for degradation of granulated sugar in an industrial scale pneumatic conveyor.

Mathematical modelling: a laser soldering process for an optoelectronics butterfly package

For sensitive optoelectronic components, traditional soldering techniques cannot be used because of their inherent sensitivity to thermal stresses. One such component is the Optoelectronic Butterfly Package which houses a laser diode chip aligned to a fibre-optic cable. Even sub-micron misalignment of the fibre optic and laser diode chip can significantly reduce the performance of the device. The high cost of each unit requires that the number of damaged components, via the laser soldering process, are kept to a minimum. Mathematical modelling is undertaken to better understand the laser soldering process and to optimize operational parameters such as solder paste volume, copper pad dimensions, laser solder times for each joint, laser intensity and absorption coefficient. Validation of the model against experimental data will be completed, and will lead to an optimization of the assembly process, through an iterative modelling cycle. This will ultimately reduce costs, improve the process development time and increase consistency in the laser soldering process.

Computational modelling of bubbles, droplets and particles in metals reduction and refining

A multi-phase framework is typically required for the CFD modelling of metals reduction processes. Such processes typically involve the interaction of liquid metals, a gas (often air) top space, liquid droplets in the top space and injection of both solid particles and gaseous bubbles into the bath. The exchange of mass, momentum and energy between the phases is fundamental to these processes. Multi-phase algorithms are complex and can be unreliable in terms of either or both convergence behaviour or in the extent to which the physics is captured. In this contribution, we discuss these multi-phase flow issues and describe an example of each of the main “single phase” approaches to modelling this class of problems (i.e., Eulerian–Lagrangian and Eulerian–Eulerian). Their utility is illustrated in the context of two problems – one involving the injection of sparging gases into a steel continuous slab caster and the other based on the development of a novel process for aluminium electrolysis. In the steel caster, the coupling of the Lagrangian tracking of the gas phase with the continuum enables the simulation of the transient motion of the metal–flux interface. The model of the electrolysis process employs a novel method for the calculation of slip velocities of oxygen bubbles, resulting from the dissolution of alumina, which allows the efficiency of the process to be predicted.

Computational modelling of multi-physics and multi-scale processes in parallel

This paper provides an overview of the developing needs for simulation software technologies for the computational modelling of problems that involve combinations of interactions amongst varying physical phenomena over a variety of time and space scales. Computational modelling of such problems requires software tech1nologies that enable the mathematical description of the interacting physical phenomena together with the solution of the resulting suites of equations in a numerically consistent and compatible manner. This functionality requires the structuring of simulation modules for specific physical phenomena so that the coupling can be effectively represented. These multi-physics and multi-scale computations are very compute intensive and the simulation software must operate effectively in parallel if it is to be used in this context. An approach to these classes of multi-disciplinary simulation in parallel is described, with some key examples of application to2 challenging engineering problems.