Inverse analysis FOR two-dimensional structures using the boundary element method

Inverse analysis is currently an important subject of study in several fields of science and engineering. The identification of physical and geometric parameters using experimental measurements is required in many applications. In this work a boundary element formulation to identify boundary and interface values as well as material properties is proposed. In particular the proposed formulation is dedicated to identifying material parameters when a cohesive crack model is assumed for 2D problems. A computer code is developed and implemented using the BEM multi-region technique and regularisation methods to perform the inverse analysis. Several examples are shown to demonstrate the efficiency of the proposed model. (C) 2010 Elsevier Ltd. All rights reserved,

A methodological approach for the performance optimization of transient seepage models through inverse analysis

The topic of the Ph.D project focuses on the modelling of the soil-water dynamics inside an instrumented embankment section along Secchia River (Cavezzo (MO)) in the period from 2017 to 2018 and the quantification of the performance of the direct and indirect simulations . The commercial code Hydrus2D by Pc-Progress has been chosen to run the direct simulations. Different soil-hydraulic models have been adopted and compared. The parameters of the different hydraulic models are calibrated using a local optimization method based on the Levenberg - Marquardt algorithm implemented in the Hydrus package. The calibration program is carried out using different types of dataset of observation points, different weighting distributions, different combinations of optimized parameters and different initial sets of parameters. The final goal is an in-depth study of the potentialities and limits of the inverse analysis when applied to a complex geotechnical problem as the case study. The second part of the research focuses on the effects of plant roots and soil-vegetation-atmosphere interaction on the spatial and temporal distribution of pore water pressure in soil. The investigated soil belongs to the West Charlestown Bypass embankment, Newcastle, Australia, that showed in the past years shallow instabilities and the use of long stem planting is intended to stabilize the slope. The chosen plant species is the Malaleuca Styphelioides, native of eastern Australia. The research activity included the design and realization of a specific large scale apparatus for laboratory experiments. Local suction measurements at certain intervals of depth and radial distances from the root bulb are recorded within the vegetated soil mass under controlled boundary conditions. The experiments are then reproduced numerically using the commercial code Hydrus 2D. Laboratory data are used to calibrate the RWU parameters and the parameters of the hydraulic model.

On the analysis of notched concrete beams: From measurement with digital image correlation to identification with boundary element method of a cohesive model

A way of coupling digital image correlation (to measure displacement fields) and boundary element method (to compute displacements and tractions along a crack surface) is presented herein. It allows for the identification of Young`s modulus and fracture parameters associated with a cohesive model. This procedure is illustrated to analyze the latter for an ordinary concrete in a three-point bend test on a notched beam. In view of measurement uncertainties, the results are deemed trustworthy thanks to the fact that numerous measurement points are accessible and used as entries to the identification procedure. (C) 2010 Elsevier Ltd. All rights reserved.

Analysis of PVC pipe-wall viscoelasticity during water hammer

This research work focuses on the analysis of hydraulic transients in polyvinyl chloride (PVC) pipes, which are characterized by a viscoelastic rheological behavior. Transient pressure data were collected in a pipe rig consisting of a set of PVC pipes. The creep function of the PVC pipes was determined by using an inverse transient model based on collected transient pressure data and compared with that obtained by carrying out mechanical tensile tests of PVC pipe specimens. The numerical results obtained from the transient solver have shown that the attenuation, dispersion, and shape of transient pressures were well described. The incorporation of the viscoelastic mechanical behavior in the hydraulic transient model has provided an excellent fitting between numerical results and observed data. Calibrated creep function based on inverse analysis fit the one determined by mechanical tests well, which emphasized the importance of pipe-wall viscoelasticity in hydraulic transients in PVC pipes.

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.

Inverse convection problem of simultaneous estimation of two boundary heat fluxes in parallel plate channels

This paper deals with the use of the conjugate gradient method of function estimation for the simultaneous identification of two unknown boundary heat fluxes in parallel plate channels. The fluid flow is assumed to be laminar and hydrodynamically developed. Temperature measurements taken inside the channel are used in the inverse analysis. The accuracy of the present solution approach is examined by using simulated measurements containing random errors, for strict cases involving functional forms with discontinuities and sharp-corners for the unknown functions. Three different types of inverse problems are addressed in the paper, involving the estimation of: (i) Spatially dependent heat fluxes; (ii) Time-dependent heat fluxes; and (iii) Time and spatially dependent heat fluxes.

Avoiding Overfitting in Inverse Modeling of Chloride Migration in Concrete

Inverse analysis for reactive transport of chlorides through concrete in the presence of electric field is presented. The model is solved using MATLAB’s built-in solvers “pdepe.m” and “ode15s.m”. The results from the model are compared with experimental measurements from accelerated migration test and a function representing the lack of fit is formed. This function is optimised with respect to varying amount of key parameters defining the model. Levenberg-Marquardt trust-region optimisation approach is employed. The paper presents a method by which the degree of inter-dependency between parameters and sensitivity (significance) of each parameter towards model predictions can be studied on models with or without clearly defined governing equations. Eigen value analysis of the Hessian matrix was employed to investigate and avoid over-parametrisation in inverse analysis. We investigated simultaneous fitting of parameters for diffusivity, chloride binding as defined by Freundlich isotherm (thermodynamic) and binding rate (kinetic parameter). Fitting of more than 2 parameters, simultaneously, demonstrates a high degree of parameter inter-dependency. This finding is significant as mathematical models for representing chloride transport rely on several parameters for each mode of transport (i.e., diffusivity, binding, etc.), which combined may lead to unreliable simultaneous estimation of parameters.

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.

Thin panels of cement composites reinforced with recycled fibres for the shear strengthening of reinforced concrete elements

A new technique was developed for producing thin panels of a cement based material reinforced with relatively high content of steel fibres originated from the industry of tyre recycling. Flexural tests with notched and un-notched specimens were carried out to characterize the mechanical properties of this Fibre Reinforced Cement Composite (FRCC) and the results are presented and discussed. The values of the fracture mode I parameters of the developed FRCC were determined by performing inverse analysis with test results obtained in three point notched beam bending tests. To appraise the potentialities of these FRCC panels for the increase of the shear capacity of reinforced (RC) beams, numerical research was performed on the use of developed FRCC panel for shear reinforcement by applying the panels in the lateral faces of RC beams deficiently reinforced in shear.

Tensile stress–crack width law for steel fibre reinforced self-compacting concrete obtained from indirect (splitting) tensile tests

In this work, the fracture mode I parameters of steel fibre reinforced self-compacting concrete (SFRSCC) were derived from the numerical simulation of indirect splitting tensile tests. The combined experimental and numerical research allowed a comparison between the stress-crack width (σ - w) relationship acquired straightforwardly from direct tensile tests, and the σ - w response derived from inverse analysis of the splitting tensile tests results. For this purpose a comprehensive nonlinear 3D finite element (FE) modeling strategy was developed. A comparison between the experimental results obtained from splitting tensile tests and the corresponding FE simulations confirmed the good accuracy of the proposed strategy to derive the σ – w for these composites. It is concluded that the post-cracking tensile laws obtained from inverse analysis provided a close relationship with the ones obtained from the experimental uniaxial tensile tests.

Recycled steel fibre reinforced concrete failing in bending and in shear

Recent research is showing that the addition of Recycled Steel Fibres (RSF) from wasted tyres can decrease significantly the brittle behaviour of cement based materials, by improving its toughness and post-cracking resistance. In this sense, Recycled Steel Fibre Reinforced Concrete (RSFRC) seems to have the potential to constitute a sustainable material for structural and non-structural applications. To assess this potential, experimental and numerical research was performed on the use of RSFRC in elements failing in bending and in beams failing in shear. The values of the fracture mode I parameters of the developed RSFRC were determined by performing inverse analysis with test results obtained in three point notched beam bending tests. To assess the possibility of using RSF as shear reinforcement in Reinforced Concrete (RC) beams, three point bending tests were executed with three series of RSFRC beams flexurally reinforced with a relatively high reinforcement ratio of longitudinal steel bars in order to assure shear failure for all the tested beams. By performing material nonlinear simulations with a computer program based on the finite element method (FEM), the applicability of the fracture mode I crack constitutive law derived from the inverse analysis is assessed for the prediction of the behaviour of these beams. The performance of the formulation proposed by RILEM TC 162 TDF and CEB-FIP 2010 for the prediction of the shear resistance of fibre reinforced concrete elements was also evaluated.

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.

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.

Numerical evaluation of the modulus of longitudinal elasticity in structural round timber elements of the Eucalyptus genus

Currently, the standards that deal with the determination of the properties of rigidity and strength for structural round timber elements do not take in consideration in their calculations and mathematical models the influence of the existing irregularities in the geometry of these elements. This study has as objective to determine the effective value of the modulus of longitudinal elasticity for structural round timber pieces of the Eucalyptus citriodora genus by a technique of optimization allied to the Inverse Analysis Method, to the Finite Element Method and the Least Square Method.

Developing a reliable strategy to infer the effective soil hydraulic properties from field evaporation experiments for agro-hydrological models

The Richards equation has been widely used for simulating soil water movement. However, the take-up of agro-hydrological models using the basic theory of soil water flow for optimizing irrigation, fertilizer and pesticide practices is still low. This is partly due to the difficulties in obtaining accurate values for soil hydraulic properties at a field scale. Here, we use an inverse technique to deduce the effective soil hydraulic properties, based on measuring the changes in the distribution of soil water with depth in a fallow field over a long period, subject to natural rainfall and evaporation using a robust micro Genetic Algorithm. A new optimized function was constructed from the soil water contents at different depths, and the soil water at field capacity. The deduced soil water retention curve was approximately parallel but higher than that derived from published pedo-tranfer functions for a given soil pressure head. The water contents calculated from the deduced soil hydraulic properties were in good agreement with the measured values. The reliability of the deduced soil hydraulic properties was tested in reproducing data measured from an independent experiment on the same soil cropped with leek. The calculation of root water uptake took account for both soil water potential and root density distribution. Results show that the predictions of soil water contents at various depths agree fairly well with the measurements, indicating that the inverse analysis is an effective and reliable approach to estimate soil hydraulic properties, and thus permits the simulation of soil water dynamics in both cropped and fallow soils in the field accurately. (C) 2009 Elsevier B.V. All rights reserved.