A BEM model applied to failure analysis of multi-fractured structures

Due to manufacturing or damage process, brittle materials present a large number of micro-cracks which are randomly distributed. The lifetime of these materials is governed by crack propagation under the applied mechanical and thermal loadings. In order to deal with these kinds of materials, the present work develops a boundary element method (BEM) model allowing for the analysis of multiple random crack propagation in plane structures. The adopted formulation is based on the dual BEM, for which singular and hyper-singular integral equations are used. An iterative scheme to predict the crack growth path and crack length increment is proposed. This scheme enables us to simulate the localization and coalescence phenomena, which are the main contribution of this paper. Considering the fracture mechanics approach, the displacement correlation technique is applied to evaluate the stress intensity factors. The propagation angle and the equivalent stress intensity factor are calculated using the theory of maximum circumferential stress. Examples of multi-fractured domains, loaded up to rupture, are considered to illustrate the applicability of the proposed method. (C) 2011 Elsevier Ltd. All rights reserved.

A numerical model for the description of the nonlinear behaviour of multi-leaf masonry walls

A nonlinear finite element model was developed to simulate the nonlinear response of three-leaf masonry specimens, which were subjected to laboratory tests with the aim of investigating the mechanical behaviour of multiple-leaf stone masonry walls up to failure. The specimens consisted of two external leaves made of stone bricks and mortar joints, and an internal leaf in mortar and stone aggregate. Different loading conditions, typologies of the collar joints, and stone types were taken into account. The constitutive law implemented in the model is characterized by a damage tensor, which allows the damage-induced anisotropy accompanying the cracking process to be described. To follow the post-peak behaviour of the specimens with sufficient accuracy it was necessary to make the damage model non-local, to avoid mesh-dependency effects related to the strain-softening behaviour of the material. Comparisons between the predicted and measured failure loads are quite satisfactory in most of the studied cases. (c) 2007 Elsevier Ltd. All rights reserved.

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.

Constitutive model for long term municipal solid waste mechanical behavior

This paper presents some improvements in the model proposed by Machado et al. [Machado SL, Carvalho MF, Vilar OM. Constitutive model for municipal solid waste. J Geotech Geoenviron Eng ASCE 2002; 128(11):940-51] now considering the influence of biodegradation of organic matter in the mechanical behavior of municipal solid waste. The original framework considers waste as composed of two component groups; fibers and organic paste. The particular laws of behavior are assessed for each component group and then coupled to represent waste behavior. The improvements introduced in this paper take into account the changes in the properties of fibers and mass loss due to organic matter depletion over time. Mass loss is indirectly calculated considering the MSW gas generation potential through a first order decay model. It is shown that as the biodegradation process occurs the proportion of fibers increases, however, they also undergo a degradation process which tends to reduce their ultimate tensile stress and Young modulus. The way these changes influence the behavior of MSW is incorporated in the final framework which captures the main features of the MSW stress-strain behavior under different loading conditions. (C) 2007 Elsevier Ltd. All rights reserved.

Defining Nutrient and Biochemical Oxygen Demand Baselines for Tropical Rivers and Streams in So Paulo State (Brazil): A Comparison Between Reference and Impacted Sites

Determining reference concentrations in rivers and streams is an important tool for environmental management. Reference conditions for eutrophication-related water variables are unavailable for Brazilian freshwaters. We aimed to establish reference baselines for So Paulo State tropical rivers and streams for total phosphorus (TP) and nitrogen (TN), nitrogen-ammonia (NH(4) (+)) and Biochemical Oxygen Demand (BOD) through the best professional judgment and the trisection methods. Data from 319 sites monitored by the So Paulo State Environmental Company (2005 to 2009) and from the 22 Water Resources Management Units in So Paulo State were assessed (N = 27,131). We verified that data from different management units dominated by similar land cover could be analyzed together (Analysis of Variance, P = 0.504). Cumulative frequency diagrams showed that industrialized management units were characterized by the worst water quality (e.g. average TP of 0.51 mg/L), followed by agricultural watersheds. TN and NH(4) (+) were associated with urban percentages and population density (Spearman Rank Correlation Test, P < 0.05). Best professional judgment and trisection (median of lower third of all sites) methods for determining reference concentrations showed agreement: 0.03 & 0.04 mg/L (TP), 0.31 & 0.34 mg/L (TN), 0.06 & 0.10 mg-N/L (NH(4) (+)) and 2 & 2 mg/L (BOD), respectively. Our reference concentrations were similar to TP and TN reference values proposed for temperate water bodies. These baselines can help with water management in So Paulo State, as well as providing some of the first such information for tropical ecosystems.

Anaerobic sequencing batch biofilm reactor applied to automobile industry wastewater treatment: Volumetric loading rate and feed strategy effects

This paper presents a technological viability study of wastewater treatment in an automobile industry by an anaerobic sequencing batch biofilm reactor containing immobilized biomass (AnSBBR) with a draft tube. The reactor was operated in 8-h cycles, with agitation of 400 rpm, at 30 degrees C and treating 2.0 L wastewater per cycle. Initially the efficiency and stability of the reactor were studied when supplied with nutrients and alkalinity. Removal efficiency of 88% was obtained at volumetric loading rate (VLR) of 3.09 mg COD/L day. When VLR was increased to 6.19 mg COD/L day the system presented stable operation with reduction in efficiency of 71%. In a second stage the AnSBBR was operated treating wastewater in natura, i.e., without nutrients supplementation, only with alkalinity, thereby changing feed strategy. The first strategy consisted in feeding 2.0 L batch wise (10 min), the second in feeding 1.0 L of influent batch wise (10 min) and an additional 1.0 L fed-batch wise (4 h), both dewatering 2.0 L of the effluent in 10 min. The third one maintained 1.0 L of treated effluent in the reactor, without discharging, and 1.0 L of influent was fed fed-batch wise (4 h) with dewatering 1.0 L of the effluent in 10 min. For all implemented strategies (VLR of 1.40, 2.57 and 2.61 mg COD/L day) the system presented stability and removal efficiency of approximately 80%. These results show that the AnSBBR presents operational flexibility, as the influent can be fed according to industry availability. In industrial processes this is a considerable advantage, as the influent may be prone to variations. Moreover, for all the investigated conditions the kinetic parameters were obtained from fitting a first-order model to the profiles of organic matter, total volatile acids and methane concentrations. Analysis of the kinetic parameters showed that the best strategy is feeding 1.0 L of influent batchwise (10 min) and 1.0 L fed-batch wise (4 h) in 8-h cycle. (c) 2007 Elsevier B.V. All rights reserved.

An electromechanical finite element model for piezoelectric energy harvester plates

Vibration-based energy harvesting has been investigated by several researchers over the last decade. The goal in this research field is to power small electronic components by converting the waste vibration energy available in their environment into electrical energy. Recent literature shows that piezoelectric transduction has received the most attention for vibration-to-electricity conversion. In practice, cantilevered beams and plates with piezoceramic layers are employed as piezoelectric energy harvesters. The existing piezoelectric energy harvester models are beam-type lumped parameter, approximate distributed parameter and analytical distributed parameter solutions. However, aspect ratios of piezoelectric energy harvesters in several cases are plate-like and predicting the power output to general (symmetric and asymmetric) excitations requires a plate-type formulation which has not been covered in the energy harvesting literature. In this paper. an electromechanically coupled finite element (FE) plate model is presented for predicting the electrical power output of piezoelectric energy harvester plates. Generalized Hamilton`s principle for electroelastic bodies is reviewed and the FE model is derived based on the Kirchhoff plate assumptions as typical piezoelectric energy harvesters are thin structures. Presence of conductive electrodes is taken into account in the FE model. The predictions of the FE model are verified against the analytical solution for a unimorph cantilever and then against the experimental and analytical results of a bimorph cantilever with a tip mass reported in the literature. Finally, an optimization problem is solved where the aluminum wing spar of an unmanned air vehicle (UAV) is modified to obtain a generator spar by embedding piezoceramics for the maximum electrical power without exceeding a prescribed mass addition limit. (C) 2009 Elsevier Ltd. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part I - A two-phase flow pattern map and a flow pattern based phenomenological model for two-phase flow frictional pressure drops

An updated flow pattern map was developed for CO2 on the basis of the previous Cheng-Ribatski-Wojtan-Thome CO2 flow pattern map [1,2] to extend the flow pattern map to a wider range of conditions. A new annular flow to dryout transition (A-D) and a new dryout to mist flow transition (D-M) were proposed here. In addition, a bubbly flow region which generally occurs at high mass velocities and low vapor qualities was added to the updated flow pattern map. The updated flow pattern map is applicable to a much wider range of conditions: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to +25 degrees C (reduced pressures from 0.21 to 0.87). The updated flow pattern map was compared to independent experimental data of flow patterns for CO2 in the literature and it predicts the flow patterns well. Then, a database of CO2 two-phase flow pressure drop results from the literature was set up and the database was compared to the leading empirical pressure drop models: the correlations by Chisholm [3], Friedel [4], Gronnerud [5] and Muller-Steinhagen and Heck [6], a modified Chisholm correlation by Yoon et al. [7] and the flow pattern based model of Moreno Quiben and Thome [8-10]. None of these models was able to predict the CO2 pressure drop data well. Therefore, a new flow pattern based phenomenological model of two-phase flow frictional pressure drop for CO2 was developed by modifying the model of Moreno Quiben and Thome using the updated flow pattern map in this study and it predicts the CO2 pressure drop database quite well overall. (C) 2007 Elsevier Ltd. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part II - An updated general flow boiling heat transfer model based on flow patterns

Corresponding to the updated flow pattern map presented in Part I of this study, an updated general flow pattern based flow boiling heat transfer model was developed for CO2 using the Cheng-Ribatski-Wojtan-Thome [L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside horizontal tubes, Int. J. Heat Mass Transfer 49 (2006) 4082-4094; L. Cheng, G. Ribatski, L. Wojtan, J.R. Thome, Erratum to: ""New flow boiling heat transfer model and flow pattern map for carbon dioxide evaporating inside tubes"" [Heat Mass Transfer 49 (21-22) (2006) 4082-4094], Int. J. Heat Mass Transfer 50 (2007) 391] flow boiling heat transfer model as the starting basis. The flow boiling heat transfer correlation in the dryout region was updated. In addition, a new mist flow heat transfer correlation for CO2 was developed based on the CO2 data and a heat transfer method for bubbly flow was proposed for completeness sake. The updated general flow boiling heat transfer model for CO2 covers all flow regimes and is applicable to a wider range of conditions for horizontal tubes: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to 25 degrees C (reduced pressures from 0.21 to 0.87). The updated general flow boiling heat transfer model was compared to a new experimental database which contains 1124 data points (790 more than that in the previous model [Cheng et al., 2006, 2007]) in this study. Good agreement between the predicted and experimental data was found in general with 71.4% of the entire database and 83.2% of the database without the dryout and mist flow data predicted within +/-30%. However, the predictions for the dryout and mist flow regions were less satisfactory due to the limited number of data points, the higher inaccuracy in such data, scatter in some data sets ranging up to 40%, significant discrepancies from one experimental study to another and the difficulties associated with predicting the inception and completion of dryout around the perimeter of the horizontal tubes. (C) 2007 Elsevier Ltd. All rights reserved.

Model to predict the chloride penetration in concrete with profile forming peak: verification of its precision

Chloride attack in marine environments or in structures where deicing salts are used will not always show profiles with concentrations that decrease from the external surface to the interior of the concrete. Some profiles show an increase in chloride concentrations from when a peak is formed. This type of profile must be analyzed in a different way from the traditional model of Fick`s second law to generate more precise service life models. A model for forecasting the penetration of chloride ions as a function of time for profiles having formed a peak. To confirm the efficiency of this model, it is necessary to observe the behavior of a chloride profile with peak in a specific structure over a period of time. To achieve this, two chloride profiles with different ages (22 and 27 years) were extracted from the same structure. The profile obtained from the 22-year sample was used to estimate the chloride profile at 27 years using three models: a) the traditional model using Fick`s second law and extrapolating the value of C(S)-external surface chloride concentration; b) the traditional model using Fick`s second law and shifting the x-axis to the peak depth; c) the previously proposed model. The results from these models were compared with the actual profile measured in the 27-year sample and the results were analyzed. The model was presented with good precision for this study of case, requiring to be tested with other structures in use.

A reference business process for ecological niche modelling

Ecological niche modelling combines species occurrence points with environmental raster layers in order to obtain models for describing the probabilistic distribution of species. The process to generate an ecological niche model is complex. It requires dealing with a large amount of data, use of different software packages for data conversion, for model generation and for different types of processing and analyses, among other functionalities. A software platform that integrates all requirements under a single and seamless interface would be very helpful for users. Furthermore, since biodiversity modelling is constantly evolving, new requirements are constantly being added in terms of functions, algorithms and data formats. This evolution must be accompanied by any software intended to be used in this area. In this scenario, a Service-Oriented Architecture (SOA) is an appropriate choice for designing such systems. According to SOA best practices and methodologies, the design of a reference business process must be performed prior to the architecture definition. The purpose is to understand the complexities of the process (business process in this context refers to the ecological niche modelling problem) and to design an architecture able to offer a comprehensive solution, called a reference architecture, that can be further detailed when implementing specific systems. This paper presents a reference business process for ecological niche modelling, as part of a major work focused on the definition of a reference architecture based on SOA concepts that will be used to evolve the openModeller software package for species modelling. The basic steps that are performed while developing a model are described, highlighting important aspects, based on the knowledge of modelling experts. In order to illustrate the steps defined for the process, an experiment was developed, modelling the distribution of Ouratea spectabilis (Mart.) Engl. (Ochnaceae) using openModeller. As a consequence of the knowledge gained with this work, many desirable improvements on the modelling software packages have been identified and are presented. Also, a discussion on the potential for large-scale experimentation in ecological niche modelling is provided, highlighting opportunities for research. The results obtained are very important for those involved in the development of modelling tools and systems, for requirement analysis and to provide insight on new features and trends for this category of systems. They can also be very helpful for beginners in modelling research, who can use the process and the experiment example as a guide to this complex activity. (c) 2008 Elsevier B.V. All rights reserved.

Adaptable Strut-and-Tie Model for Design and Verification of Four-Pile Caps

A large percentage of pile caps support only one column, and the pile caps in turn are supported by only a few piles. These are typically short and deep members with overall span-depth ratios of less than 1.5. Codes of practice do not provide uniform treatment for the design of these types of pile caps. These members have traditionally been designed as beams spanning between piles with the depth selected to avoid shear failures and the amount of longitudinal reinforcement selected to provide sufficient flexural capacity as calculated by the engineering beam theory. More recently, the strut-and-tie method has been used for the design of pile caps (disturbed or D-region) in which the load path is envisaged to be a three-dimensional truss, with compressive forces being supported by concrete compressive struts between the column and piles and tensile forces being carried by reinforcing steel located between piles. Both of these models have not provided uniform factors of safety against failure or been able to predict whether failure will occur by flexure (ductile mode) or shear (fragile mode). In this paper, an analytical model based on the strut-and-tie approach is presented. The proposed model has been calibrated using an extensive experimental database of pile caps subjected to compression and evaluated analytically for more complex loading conditions. It has been proven to be applicable across a broad range of test data and can predict the failures modes, cracking, yielding, and failure loads of four-pile caps with reasonable accuracy.

PERFORMANCE OF THE MODEL CALLED ACQUANETXL IN THE ALLOCATION OF WATER IN COMPLEX WATER RESOURCE SYSTEMS

This article presents a tool for the allocation analysis of complex systems of water resources, called AcquaNetXL, developed in the form of spreadsheet in which a model of linear optimization and another nonlinear were incorporated. The AcquaNetXL keeps the concepts and attributes of a decision support system. In other words, it straightens out the communication between the user and the computer, facilitates the understanding and the formulation of the problem, the interpretation of the results and it also gives a support in the process of decision making, turning it into a clear and organized process. The performance of the algorithms used for solving the problems of water allocation was satisfactory especially for the linear model.

A transient three-dimensional heat transfer model of the human body

The objective of this work is to develop an improved model of the human thermal system. The features included are important to solve real problems: 3D heat conduction, the use of elliptical cylinders to adequately approximate body geometry, the careful representation of tissues and important organs, and the flexibility of the computational implementation. Focus is on the passive system, which is composed by 15 cylindrical elements and it includes heat transfer between large arteries and veins. The results of thermal neutrality and transient simulations are in excellent agreement with experimental data, indicating that the model represents adequately the behavior of the human thermal system. (C) 2009 Elsevier Ltd. All rights reserved.