An Attempt to Predict the Failure Time of Abandoned Mine Pillars

Summary  Left in place pillars of abandoned mines are subject to weathering (e.g action of water, bacteria) that degrades their mechanical strength and eventually leads to collapse. A simple weathering model is proposed, that is governed by two parameters: the rate of progression of weathered front and the rate of degradation of the compressive strength with time. Both plane strain and axisymmetric analyses are performed and closed form solutions of the variation with time of the bearing capacity of the pillar are given. Experimental data of the tests conducted on gypsum and anhydrite specimens attacked by water are presented. It is shown that in order to fit the experimental data a third parameter must be introduced. New closed form solutions are given and the data are used for estimating the time to failure of abandoned gypsum mines in Northern Italy. Authors’ address: Dr. Riccardo Castellanza, Research Assistant, Department of Structural Engineering, Milan University of Technology (Politecnico), Piazza Leonardo da Vinci 32, 20133 Milan, Italy     

An elastoplastic strainhardening model for soil allowing for hydraulic bonding–debonding effects

The paper presents a strainhardening constitutive model for unsaturated soil behaviour based on energy conjugated stress variables in the framework of superposed continua. The proposed constitutive law deals with hydro‐mechanical coupling phenomena. The main purpose is to develop within a consistent framework a model that can deal with possible mechanical instabilities occurring in partially saturated materials. The loss of capillary effects during wetting processes can, in fact, play a central role in unstable processes. Therefore, it will be shown that the bonding effects due to surface tensions can be described in a mathematical framework similar to that employed for bonded geomaterials to model weathering or diagenesis effects, either mechanically or chemically induced. The results of several simulations of common laboratory tests on partially saturated soil specimens are shown. The calculated behaviour appears to be in good qualitative agreement with that observed in the laboratory. In particular it is shown that volumetric collapse phenomena due to hydraulic debonding effects can be successfully described by the model. Finally, it will be highlighted the ability of the model to naturally capture the transition to a fully saturated condition and to deal with possible mechanical instabilities in the unsaturated regime. Copyright © 2008 John Wiley & Sons, Ltd.     

A Generalized Backward Euler algorithm for the numerical integration of an isotropic hardening elastoplastic model for mechanical and chemical degradation of bonded geomaterials

An extended version of the classical Generalized Backward Euler (GBE) algorithm is proposed for the numerical integration of a three‐invariant isotropic‐hardening elastoplastic model for cemented soils or weak rocks undergoing mechanical and non‐mechanical degradation processes. The restriction to isotropy allows to formulate the return mapping algorithm in the space of principal elastic strains. In this way, an efficient and robust integration scheme is developed which can be applied to relatively complex yield surface and plastic potential functions. Moreover, the proposed algorithm can be linearized in closed form, thus allowing for quadratic convergence in the global Newton iteration. A series of numerical experiments are performed to illustrate the accuracy and convergence properties of the algorithm. Selected results from a finite element analysis of a circular footing on a soft rock layer undergoing chemical weathering are then presented to illustrate the algorithm performance at the boundary value problem level. Copyright © 2002 John Wiley & Sons, Ltd.     

A constitutive model for sand in triaxial compression

Although there is a large number of constitutive models for sand available in the literature it is believed that a fresh approach, striking a balance between complexity and theoretical rigour, is desirable. The approach here has certain conceptual links with the Cam Clay series of elastic–plastic models, but includes the more general starting assumption that the yield function, plastic potential and failure locus should be given quite distinct mathematical expressions. Possible physical bases for the proposed forms are discussed. Ways in which the parameters required to define the model may be determined are suggested and the use of the model is then demonstrated. Firstly, it is shown that, where a limited set of experimental data is available, the model is flexible enough to be able to match the test results. Secondly, it is shown that, where a wide range of test results has been produced, it is possible to determine the model constitutive parameters from a small number of tests and proceed to make satisfactory predictionsfor other, quite different, types of test. The model is developed for sand at a single initial density, but the way in which the constitutive parameters might be expected to vary with density is discussed. The model is described for conditions of triaxial compression, and extension to more general stress states will be needed before it can be put to the test of incorporation in, for example, a finite element program.     

Do mixing models with different input requirement yield similar streamflow source contributions? Case study: A tropical montane catchment

Hydrogeochemical based mixing models have been successfully used to investigate the composition and source identification of streamflow. The applicability of these models is limited due to the high costs associated with data collection and the hydrogeochemical analysis of water samples. Fortunately, a variety of mixing models exist, requiting different amount of data as input, and in data scarce regions it is likely that preference will be given to models with the lowest requirement of input data. An unanswered question is if models with high or low input requirement are equally accurate. To this end, the performance of two mixing models with different input requirement, the mixing model analysis (MMA) and the end-member mixing analysis (EMMA), were verified on a tropical montane headwater catchment (21.7 km²) in the Ecuadorian Andes. Nineteen hydrogeochemical tracers were measured on water samples collected weekly during 3 years in streamflow and eight potential water sources or end-members (precipitation, lake water, soil water from different horizons and springs). Results based on 6 conservative tracers, revealed that EMMA (using all tracers) and MMA (using pair-combinations out of the 6 conservative ones), identified the same end-members: rainfall, soil water and spring water., as well as, similar contribution fractions to streamflow from rainfall 21.9% and 21.4%, soil water 52.7% and 52.3%, and spring water 26.1% and 28.7%, respectively. Our findings show that a hydrogeochemical mixing model requiring a few tracers can provide similar outcomes than models demanding more tracers as input data. This underlines the value of a preliminary detailed hydrogeochemical characterization as basis to derive the most cost-efficient monitoring strategy.     

Water transport and tracer mixing in volcanic ash soils at a tropical hillslope: A wet layered sloping sponge

Andosol soils formed in volcanic ash provide key hydrological services in montane environments. To unravel the subsurface water transport and tracer mixing in these soils we conducted a detailed characterization of soil properties and analyzed a 3-year data set of sub-hourly hydrometric and weekly stable isotope data collected at three locations along a steep hillslope. A weakly developed (52–61 cm depth), highly organic andic (Ah) horizon overlaying a mineral (C) horizon was identified, both showing relatively similar properties and subsurface flow dynamics along the hillslope. Soil moisture observations in the Ah horizon showed a fast responding (few hours) “rooted” layer to a depth of 15 cm, overlying a “perched” layer that remained near saturated year-round. The formation of the latter results from the high organic matter (33–42%) and clay (29–31%) content of the Ah horizon and an abrupt hydraulic conductivity reduction in this layer with respect to the rooted layer above. Isotopic signatures revealed that water resides within this soil horizon for short periods, both at the rooted (2 weeks) and perched (4 weeks) layer. A fast soil moisture reaction during rainfall events was also observed in the C horizon, with response times similar to those in the rooted layer. These results indicate that despite the perched layer, which helps sustain the water storage of the soil, a fast vertical mobilization of water through the entire soil profile occurs during rainfall events. The latter being the result of the fast transmissivity of hydraulic potentials through the porous matrix of the Andosols, as evidenced by the exponential shape of the water retention curves of the subsequent horizons. These findings demonstrate that the hydrological behavior of volcanic ash soils resembles that of a “layered sponge,” in which vertical flow paths dominate.     

Wavelet analyses of neural networks based river discharge decomposition

The problem of discharge forecasting using precipitation as input is still very active in Hydrology, and has a plethora of approaches to its solution. But, when the objective is to simulate discharge values without considering the phenomenology behind the processes involved, Artificial Neural Networks, ANN give good results. However, the question of how the black box internally solve this problem remains open. In this research, the classical rainfall-runoff problem is approached considering that the total discharge is a sum of components of the hydrological system, which from the ANN perspective is translated to the sum of three signals related to the fast, middle and slow flow. Thus, the present study has two aims (a) to study the time-frequency representation of discharge by an ANN hydrologic model and (b) to study the capabilities of ANN to additively decompose total river discharge. This study adds knowledge to the open problem of the physical interpretability of black-box models, which remains very limited. The results show that total discharge is adequately simulated in the time frequency domain, although less power spectrum is evident during the rainy seasons in the ANN model, due to fast flow underestimation. The wavelet spectrum of discharge represents well the slow, middle and fast flow components of the system with transit times of 256, 12–64 and 2–12 days, respectively. Interestingly, these transit times are remarkably similar to those of the soil water reservoirs of the studied system, a small headwater catchment in the tropical Andes. This result needs further research because it opens the possibility of determining MMT on a fraction of the cost of isotopic based methods. The cross-power spectrum indicates that the error in the simulated discharge is more related to the misrepresentation of the fast and the middle flow components, despite limitations in the recharge period of the slow flow component. With respect to the representation of individual signals of the slow, middle and fast flows components, the three neurons were uncapable to individually represent such flows. However, the combination of pairs of these signals resemble the dynamics and the spectral content of the aforementioned flows signals. These results show some evidence that signal processing techniques may be used to infer information about the hydrological functioning of a basin.     

Nonlinear numerical assessment of the seismic response of hillside RC buildings

Major damage has been reported in hilly areas after major earthquakes,primarily because of two special conditions:the variation in the seismic ground motion due to the inclined ground surface and the irregularities caused by a stepped base level in the structure.The aim of this study is to evaluate possible differences in the responses of Chilean hillside buildings through numerical linear-elastic and nonlinear analyses.In the first step,a set of response-spectrum analyses were performed on four simplified 2D structures with mean base inclination angles of 0°,15°,30°,and 45°.The structures were designed to comply with Chilean seismic codes and standards,and the primary response parameters were compared.To assess the seismic performance of the buildings,nonlinear static(pushover)and dynamic(time-history)analyses were performed with SeismoStruct software.Pushover analyses were used to compare the nonlinear response at the maximum roof displacement and the damage patterns.Time-history analyses were performed to assess the nonlinear dynamic response of the structures subjected to seismic ground motions modified by topographic effects.To consider the topographic modification,acceleration records were obtained from numerical models of soil,which were calculated using the rock acceleration record of the Mw 8.01985 Chilean earthquake.Minor differences in the structure responses(roof displacements and maximum element forces and moments)were caused by the topographic effects in the seismic input motion,with the highly predominant ones being the differences caused by the step-back configuration at the base of the structures.High concentrations of shear forces in short walls were observed,corresponding to the walls located in the upper zone of the foundation system.The response of the structures with higher angles was observed to be more prone to fragile failures due to the accumulation of shear forces.Even though hillside buildings gain stiffness in the lower stories,resulting in lower design roof displacement,maximum roof displacements for nonlinear time-history analyses remained very close for all the models that were primarily affected by the drifts of the lower stories.Additionally,vertical parasitic accelerations were considered for half the time-history analyses performed here.The vertical component seems to considerably modify the axial load levels in the shear walls on all stories.     

On the Measurement of the Tensile Strength of Soft Rocks

Summary. This paper reports on a comparative study of various types of experimental tests for measuring the tensile strength of rocks and rock-like materials. A critical assessment is presented of some widely used laboratory techniques on the basis of experimental data from the literature and from the laboratory investigation performed in this study.Tests were carried out using a triaxial apparatus recently set up at Milan University of Technology. This was designed to reduce random misalignments between specimen and loading frame, which are typical of conventional triaxial cells with external tie bars. The apparatus was then modified to perform various types of tests for determining the tensile strength of materials. An artificial building stone and a natural calcarenite of the Gravina di Puglia geological formation, sampled at Montescaglioso (Matera-Italy) were tested in this research. The experimental investigation included various types of tests, namely the uni-axial, the Brazilian, the ring, the three and four points bending and the Luong test. Specimens of both materials were also compressed to failure in unconfined conditions and loaded cyclically in unconfined tension and compression, Young’s moduli being measured by means of local instrumentation.