Tunnel engineering – influence of the type and the quantity of measurements in the back analysis of geomechanical parameters

The monitoring data collected during tunnel excavation can be used in inverse analysis procedures in order to identify more realistic geomechanical parameters that can increase the knowledge about the interested formations. These more realistic parameters can be used in real time to adapt the project to the real structure in situ behaviour. However, monitoring plans are normally designed for safety assessment and not especially for the purpose of inverse analysis. In fact, there is a lack of knowledge about what types and quantity of measurements are needed to succeed in identifying the parameters of interest. Also, the optimisation algorithm chosen for the identification procedure may be important for this matter. In this work, this problem is addressed using a theoretical case with which a thorough parametric study was carried out using two optimisation algorithms based on different calculation paradigms, namely a conventional gradient-based algorithm and an evolution strategy algorithm. Calculations were carried for different sets of parameters to identify several combinations of types and amount of monitoring data. The results clearly show the high importance of the available monitoring data and the chosen algorithm for the success rate of the inverse analysis process.

Joint modelling of longitudinal and survival data on breast cancer

Tese de Doutoramento em Ciências (Especialidade em Matemática)

The influence of fibre orientation on the post-cracking tensile behaviour of steel fibre reinforced self-compacting concrete

Adding fibres to concrete provides several advantages, especially in terms of controlling the crack opening width and propagation after the cracking onset. However, distribution and orientation of the fibres toward the active crack plane are significantly important in order to maximize its benefits. Therefore, in this study, the effect of the fibre distribution and orientation on the post-cracking tensile behaviour of the steel fibre reinforced self-compacting concrete (SFRSCC) specimens is investigated. For this purpose, several cores were extracted from distinct locations of a panel and were subjected to indirect (splitting) and direct tensile tests. The local stress-crack opening relationship (σ-w) was obtained by modelling the splitting tensile test under the finite element framework and by performing an Inverse Analysis (IA) procedure. Afterwards the σ-w law obtained from IA is then compared with the one ascertained directly from the uniaxial tensile tests. Finally, the fibre distribution/orientation parameters were determined adopting an image analysis technique.

A model to simulate the moment-rotation and crack width of FRC members reinforced with longitudinal bars

The present work describes a model for the determination of the moment–rotation relationship of a cross section of fiber reinforced concrete (FRC) elements that also include longitudinal bars for the flexural reinforcement (R/FRC). Since a stress–crack width relationship (σ–w)(σ–w) is used to model the post-cracking behavior of a FRC, the σ–w directly obtained from tensile tests, or derived from inverse analysis applied to the results obtained in three-point notched beam bending tests, can be adopted in this approach. For a more realistic assessment of the crack opening, a bond stress versus slip relationship is assumed to simulate the bond between longitudinal bars and surrounding FRC. To simulate the compression behavior of the FRC, a shear friction model is adopted based on the physical interpretation of the post-peak compression softening behavior registered in experimental tests. By allowing the formation of a compressive FRC wedge delimited by shear band zones, the concept of concrete crushing failure mode in beams failing in bending is reinterpreted. By using the moment–rotation relationship, an algorithm was developed to determine the force–deflection response of statically determinate R/FRC elements. The model is described in detail and its good predictive performance is demonstrated by using available experimental data. Parametric studies were executed to evidence the influence of relevant parameters of the model on the serviceability and ultimate design conditions of R/FRC elements failing in bending.

Seismic vulnerability of existing masonry buildings: Nonlinear parametric analysis

Existing masonry structures are usually associated to a high seismic vulnerability, mainly due to the properties of the materials, weak connections between floors and load-bearing walls, high mass of the masonry walls and flexibility of the floors. For these reasons, the seismic performance of existing masonry structures has received much attention in the last decades. This study presents the parametric analysis taking into account the deviations on features of the gaioleiro buildings - Portuguese building typology. The main objective of the parametric analysis is to compare the seismic performance of the structure as a function of the variations of its properties with respect to the response of a reference model. The parametric analysis was carried out for two types of structural analysis, namely for the non-linear dynamic analysis with time integration and for the pushover analysis with distribution of forces proportional to the inertial forces of the structure. The Young's modulus of the masonry walls, Young's modulus of the timber floors, the compressive and tensile non-linear properties (strength and fracture energy) were the properties considered in both type of analysis. Additionally, in the dynamic analysis, the influences of the vis-cous damping and of the vertical component of the earthquake were evaluated. A pushover analysis proportional to the modal displacement of the first mode in each direction was also carried out. The results shows that the Young's modulus of the masonry walls, the Young's modulus of the timber floors and the compressive non-linear properties are the pa-rameters that most influence the seismic performance of this type of tall and weak existing masonry structures. Furthermore, it is concluded that that the stiffness of the floors influences significantly the strength capacity and the collapse mecha-nism of the numerical model. Thus, a study on the strengthening of the floors was also carried out. The increase of the thickness of the timber floors was the strengthening technique that presented the best seismic performance, in which the reduction of the out-of-plane displacements of the masonry walls is highlighted.

Thermal conditions affecting heat transfer in FDM/FFE: a contribution towards the numerical modelling of the process

The performance of parts produced by Free Form Extrusion (FFE), an increasingly popular additive manufacturing technique, depends mainly on their dimensional accuracy, surface quality and mechanical performance. These attributes are strongly influenced by the evolution of the filament temperature and deformation during deposition and solidification. Consequently, the availability of adequate process modelling software would offer a powerful tool to support efficient process set-up and optimisation. This work examines the contribution to the overall heat transfer of various thermal phenomena developing during the manufacturing sequence, including convection and radiation with the environment, conduction with support and between adjacent filaments, radiation between adjacent filaments and convection with entrapped air. The magnitude of the mechanical deformation is also studied. Once this exercise is completed, it is possible to select the material properties, process variables and thermal phenomena that should be taken in for effective numerical modelling of FFE.

Methodology for real-time adaptation of tunnels support using the observational method

The observational method in tunnel engineering allows the evaluation in real time of the actual conditions of the ground and to take measures if its behavior deviates considerably from predictions. However, it lacks a consistent and structured methodology to use the monitoring data to adapt the support system in real time. The definition of limit criteria above which adaptation is required are not defined and complex inverse analysis procedures (Rechea et al. 2008, Levasseur et al. 2010, Zentar et al. 2001, Lecampion et al. 2002, Finno and Calvello 2005, Goh 1999, Cui and Pan 2012, Deng et al. 2010, Mathew and Lehane 2013, Sharifzadeh et al. 2012, 2013) may be needed to consistently analyze the problem. In this paper a methodology for the real time adaptation of the support systems during tunneling is presented. In a first step limit criteria for displacements and stresses are proposed. The methodology uses graphics that are constructed during the project stage based on parametric calculations to assist in the process and when these graphics are not available, since it is not possible to predict every possible scenario, inverse analysis calculations are carried out. The methodology is applied to the “Bois de Peu” tunnel which is composed by two tubes with over 500 m long. High uncertainty levels existed concerning the heterogeneity of the soil and consequently in the geomechanical design parameters. The methodology was applied in four sections and the results focus on two of them. It is shown that the methodology has potential to be applied in real cases contributing for a consistent approach of a real time adaptation of the support system and highlight the importance of the existence of good quality and specific monitoring data to improve the inverse analysis procedure.

Modelling integral viscoelastic flows in OpenFOAM

[Excerpt] A large number of constitutive equations were developed for viscoelastic fluids, some empirical and other with strong physical foundations. The currently available macroscopic constitutive equations can be divided in two main types: differential and integral. Some of the constitutive equations, e.g. Maxwell are available both in differential and integral types. However, relevant in tegral models, like K - BKZ, just possesses the integral form. (...)

Spatio-temporal modelling of environmental data

Programa Doutoral em Matemática e Aplicações.

The group inverse of a product

In this paper, we characterize the existence and give an expression of the group inverse of a product of two regular elements by means of a ring unit.

Further results on the inverse along an element in semigroups and rings

In this paper, we introduce a new notion in a semigroup $S$ as an extension of Mary's inverse. Let $a,d\in S$. An element $a$ is called left (resp. right) invertible along $d$ if there exists $b\in S$ such that $bad=d$ (resp. $dab=b$) and $b\leq_\mathcal{L}d$ (resp. $b\leq_\mathcal{R}d$). An existence criterion of this type inverse is derived. Moreover, several characterizations of left (right) regularity, left (right) $\pi$-regularity and left (right) $*$-regularity are given in a semigroup. Further, another existence criterion of this type inverse is given by means of a left (right) invertibility of certain elements in a ring. Finally we study the (left, right) inverse along a product in a ring, and, as an application, Mary's inverse along a matrix is expressed.

The inverse along a product and its applications

In this paper, we study the recently defined notion of the inverse along an element. An existence criterion for the inverse along a product is given in a ring. As applications, we present the equivalent conditions for the existence and expressions of the inverse along a matrix.

Dynamic modelling and analysis of hydraulic forces in radial blood pumps

Dissertação de mestrado integrado em Biomedical Engineering Biomaterials, Biomechanics and Rehabilitation

Modelling tyre-road noise with data mining techniques

The research aimed to establish tyre-road noise models by using a Data Mining approach that allowed to build a predictive model and assess the importance of the tested input variables. The data modelling took into account three learning algorithms and three metrics to define the best predictive model. The variables tested included basic properties of pavement surfaces, macrotexture, megatexture, and uneven- ness and, for the first time, damping. Also, the importance of those variables was measured by using a sensitivity analysis procedure. Two types of models were set: one with basic variables and another with complex variables, such as megatexture and damping, all as a function of vehicles speed. More detailed models were additionally set by the speed level. As a result, several models with very good tyre-road noise predictive capacity were achieved. The most relevant variables were Speed, Temperature, Aggregate size, Mean Profile Depth, and Damping, which had the highest importance, even though influenced by speed. Megatexture and IRI had the lowest importance. The applicability of the models developed in this work is relevant for trucks tyre-noise prediction, represented by the AVON V4 test tyre, at the early stage of road pavements use. Therefore, the obtained models are highly useful for the design of pavements and for noise prediction by road authorities and contractors.

Modelling Z→ττ processes in ATLAS with τ-embedded Z→μμ data

This paper describes the concept, technical realisation and validation of a largely data-driven method to model events with Z→ττ decays. In Z→μμ events selected from proton-proton collision data recorded at s√=8 TeV with the ATLAS experiment at the LHC in 2012, the Z decay muons are replaced by τ leptons from simulated Z→ττ decays at the level of reconstructed tracks and calorimeter cells. The τ lepton kinematics are derived from the kinematics of the original muons. Thus, only the well-understood decays of the Z boson and τ leptons as well as the detector response to the τ decay products are obtained from simulation. All other aspects of the event, such as the Z boson and jet kinematics as well as effects from multiple interactions, are given by the actual data. This so-called τ-embedding method is particularly relevant for Higgs boson searches and analyses in ττ final states, where Z→ττ decays constitute a large irreducible background that cannot be obtained directly from data control samples.