A boundary element formulation for analysis of elastoplastic plates with geometrical nonlinearity

In this paper a new boundary element method formulation for elastoplastic analysis of plates with geometrical nonlinearities is presented. The von Mises criterion with linear isotropic hardening is considered to evaluate the plastic zone. Large deflections are assumed but within the context of small strain. To derive the boundary integral equations the von Karman`s hypothesis is taken into account. An initial stress field is applied to correct the true stresses according to the adopted criterion. Isoparametric linear elements are used to approximate the boundary unknown values while triangular internal cells with linear shape function are adopted to evaluate the domain value influences. The nonlinear system of equations is solved by using an implicit scheme together with the consistent tangent operator derived along the paper. Numerical examples are presented to demonstrate the accuracy and the validity of the proposed formulation.

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.

BEM formulation for von Karman plates

This work deals with nonlinear geometric plates in the context of von Karman`s theory. The formulation is written such that only the boundary in-plane displacement and deflection integral equations for boundary collocations are required. At internal points, only out-of-plane rotation, curvature and in-plane internal force representations are used. Thus, only integral representations of these values are derived. The nonlinear system of equations is derived by approximating all densities in the domain integrals as single values, which therefore reduces the computational effort needed to evaluate the domain value influences. Hyper-singular equations are avoided by approximating the domain values using only internal nodes. The solution is obtained using a Newton scheme for which a consistent tangent operator was derived. (C) 2009 Elsevier Ltd. All rights reserved.

GENERALIZED FINITE ELEMENT METHOD FOR NONLINEAR THREE-DIMENSIONAL ANALYSIS OF SOLIDS

The Generalized Finite Element Method (GFEM) is employed in this paper for the numerical analysis of three-dimensional solids tinder nonlinear behavior. A brief summary of the GFEM as well as a description of the formulation of the hexahedral element based oil the proposed enrichment strategy are initially presented. Next, in order to introduce the nonlinear analysis of solids, two constitutive models are briefly reviewed: Lemaitre`s model, in which damage and plasticity are coupled, and Mazars`s damage model suitable for concrete tinder increased loading. Both models are employed in the framework of a nonlocal approach to ensure solution objectivity. In the numerical analyses carried out, a selective enrichment of approximation at regions of concern in the domain (mainly those with high strain and damage gradients) is exploited. Such a possibility makes the three-dimensional analysis less expensive and practicable since re-meshing resources, characteristic of h-adaptivity, can be minimized. Moreover, a combination of three-dimensional analysis and the selective enrichment presents a valuable good tool for a better description of both damage and plastic strain scatterings.

A layered finite element for reinforced concrete beams with bond-slip effects

The behaviour of reinforced concrete members is affected by the slipping of steel bars inserted in the concrete matrix. A tension-stiffening effect and crack evolution occur from the beginning of slipping; thus, the assessment of those phenomena requires the introduction of a bond-slip interaction model. This work presents a beam-layered model, including the constitutive relationships of materials and their interaction, according to the CEB-FIP Model Code 1990. To eliminate the finite element sub-division procedure, a continuous slip function is imposed into the element domain. The results are continuous descriptions of bond stress in the steel-concrete interface, as well as concrete and steel stresses along the element. (C) 2007 Elsevier Ltd. All rights reserved.

Multiple random crack propagation using a boundary element formulation

This paper proposes a boundary element method (BEM) model that is used for the analysis of multiple random crack growth by considering linear elastic fracture mechanics problems and structures subjected to fatigue. The formulation presented in this paper is based on the dual boundary element method, in which singular and hyper-singular integral equations are used. This technique avoids singularities of the resulting algebraic system of equations, despite the fact that the collocation points coincide for the two opposite crack faces. In fracture mechanics analyses, the displacement correlation technique is applied to evaluate stress intensity factors. The maximum circumferential stress theory is used to evaluate the propagation angle and the effective stress intensity factor. The fatigue model uses Paris` law to predict structural life. Examples of simple and multi-fractured structures loaded until rupture are considered. These analyses demonstrate the robustness of the proposed model. In addition, the results indicate that this formulation is accurate and can model localisation and coalescence phenomena. (C) 2010 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.

Solution of a heterogeneous modeling of a horizontal-flow anaerobic immobilized biomass (HAIB) reactor by the sequencing method

A modeling study was completed to develop a methodology that combines the sequencing and finite difference methods for the simulation of a heterogeneous model of a tubular reactor applied in the treatment of wastewater. The system included a liquid phase (convection diffusion transport) and a solid phase (diffusion reaction) that was obtained by completing a mass balance in the reactor and in the particle, respectively. The model was solved using a pilot-scale horizontal-flow anaerobic immobilized biomass (HAIB) reactor to treat domestic sewage, with the concentration results compared with the experimental data. A comparison of the behavior of the liquid phase concentration profile and the experimental results indicated that both the numerical methods offer a good description of the behavior of the concentration along the reactor. The advantage of the sequencing method over the finite difference method is that it is easier to apply and requires less computational time to model the dynamic simulation of outlet response of HAIB.

Dendritic Arm Spacing Affecting Mechanical Properties and Wear Behavior of Al-Sn and Al-Si Alloys Directionally Solidified under Unsteady-State Conditions

Alloys of Al-Sn and Al-Si are widely used in tribological applications such as cylinder liners and journal bearings. Studies of the influence of the as-cast microstructures of these alloys on the final mechanical properties and wear resistance can be very useful for planning solidification conditions in order to permit a desired level of final properties to be achieved. The aim of the present study was to contribute to a better understanding about the relationship between the scale of the dendritic network and the corresponding mechanical properties and wear behavior. The Al-Sn (15 and 20 wt pct Sn) and Al-Si (3 and 5 wt pct Si) alloys were directionally solidified under unsteady-state heat flow conditions in water-cooled molds in order to permit samples with a wide range of dendritic spacings to be obtained. These samples were subjected to tensile and wear tests, and experimental quantitative expressions correlating the ultimate tensile strength (UTS), yield tensile strength, elongation, and wear volume to the primary dendritic arm spacing (DAS) have been determined. The wear resistance was shown to be significantly affected by the scale of primary dendrite arm spacing. For Al-Si alloys, the refinement of the dendritic array improved the wear resistance, while for the Al-Sn alloys, an opposite effect was observed, i.e., the increase in primary dendrite arm spacing improved the wear resistance. The effect of inverse segregation, which is observed for Al-Sn alloys, on the wear resistance is also discussed.

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.

Multimodal passive vibration control of sandwich beams with shunted shear piezoelectric materials

This work presents a performance analysis of multimodal passive vibration control of a sandwich beam using shear piezoelectric materials, embedded in a sandwich beam core, connected to independent resistive shunt circuits. Shear piezoelectric actuators were recently shown to be more interesting for higher frequencies and stiffer structures. In particular, for shunted damping, it was shown that equivalent material loss factors of up to 31% can be achieved by optimizing the shunt circuit. In the present work, special attention is given to the design of multimodal vibration control through independent shunted shear piezoelectric sensors. In particular, a parametric analysis is performed to evaluate optimal configurations for a set of modes to be damped. Then, a methodology to evaluate the modal damping resulting from each shunted piezoelectric sensor is presented using the modal strain energy method. Results show that modal damping factors of 1%-2% can be obtained for three selected vibration modes.

Chemical composition of mixed construction and demolition recycled aggregates from the State of Sao Paulo

Construction and Demolition Waste (CDW) represents. about 50% of the total Brazilian municipal solid waste: thus, recycling represents huge benefits both in environmental and economic perspectives. Herein, the chemical characterization results of three samples from two different recycling plants from the State of Sao Paulo is prevented. The results demonstrated that the visual classification into grey and red is not related to the chemical composition but mostly to the grain size fraction. The chemical composition of the CDW varies according to the content of cement paste, natural aggregates (quartz sand or granite), red ceramic and clay. Furthermore, the production of recycled concrete aggregates requires two crushing stages to meet the technical standards. The sand fraction (below 4.8 mm) presents high grades of SiO(2), which indicates the liberation of cement paste to fines (< 0.15 mm). The fines have a great potential to be used in the cement industry.

Whirlwind: a new cryptographic hash function

A new cryptographic hash function Whirlwind is presented. We give the full specification and explain the design rationale. We show how the hash function can be implemented efficiently in software and give first performance numbers. A detailed analysis of the security against state-of-the-art cryptanalysis methods is also provided. In comparison to the algorithms submitted to the SHA-3 competition, Whirlwind takes recent developments in cryptanalysis into account by design. Even though software performance is not outstanding, it compares favourably with the 512-bit versions of SHA-3 candidates such as LANE or the original CubeHash proposal and is about on par with ECHO and MD6.

Three-phase harmonic distortion state estimation algorithm based on evolutionary strategies

This paper presents a new methodology to estimate unbalanced harmonic distortions in a power system, based on measurements of a limited number of given sites. The algorithm utilizes evolutionary strategies (ES), a development branch of evolutionary algorithms. The problem solving algorithm herein proposed makes use of data from various power quality meters, which can either be synchronized by high technology GPS devices or by using information from a fundamental frequency load flow, what makes the overall power quality monitoring system much less costly. The ES based harmonic estimation model is applied to a 14 bus network to compare its performance to a conventional Monte Carlo approach. It is also applied to a 50 bus subtransmission network in order to compare the three-phase and single-phase approaches as well as the robustness of the proposed method. (C) 2010 Elsevier B.V. All rights reserved.

Power Factor Calculation by the Finite Element Method

The use of finite element analysis (FEA) to design electrical motors has increased significantly in the past few years due the increasingly better performance of modern computers. Even though the analytical software remains the most used tool, the FEA is widely used to refine the analysis and gives the final design to be prototyped. The power factor, a standard data of motor manufactures data sheet is important because it shows how much reactive power is consumed by the motor. This data becomes important when the motor is connected to network. However, the calculation of power factor is not an easy task. Due to the saturation phenomena the input motor current has a high level of harmonics that cannot be neglected. In this work the FEA is used to evaluate a proposed (not limitative) methodology to estimate the power factor or displacement factor of a small single-phase induction motor. Results of simulations and test are compared.