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 coupled boundary element/differential equation method formulation for plate-beam interaction analysis

In this work, a new boundary element formulation for the analysis of plate-beam interaction is presented. This formulation uses a three nodal value boundary elements and each beam element is replaced by its actions on the plate, i.e., a distributed load and end of element forces. From the solution of the differential equation of a beam with linearly distributed load the plate-beam interaction tractions can be written as a function of the nodal values of the beam. With this transformation a final system of equation in the nodal values of displacements of plate boundary and beam nodes is obtained and from it, all unknowns of the plate-beam system are obtained. Many examples are analyzed and the results show an excellent agreement with those from the analytical solution and other numerical methods. (C) 2009 Elsevier Ltd. All rights reserved.

Simplified models for cross-section stress demands on C-section purlins in uplift

The objective of this paper is to provide and verify simplified models that predict the longitudinal stresses that develop in C-section purlins in uplift. The paper begins with the simple case of flexural stress: where the force has to be applied at the shear center, or the section braced in both flanges. Restrictions on load application point and restraint of the flanges are removed until arriving at the more complex problem of bending when movement of the tension flange alone is restricted, as commonly found in purlin-sheeting systems. Winter`s model for predicting the longitudinal stresses developed due to direct torsion is reviewed, verified, and then extended to cover the case of a bending member with tension flange restraint. The developed longitudinal stresses from flexure and restrained torsion are used to assess the elastic stability behavior of typical purlin-sheeting systems. Finally, strength predictions of typical C-section purlins are provided for existing AISI methods and a newly proposed extension to the direct strength method that employs the predicted longitudinal stress distributions within the strength prediction. (C) 2009 Elsevier Ltd. All rights reserved.

A hybrid-Trefftz formulation for plane elasticity with selective enrichment of the approximations

This work is related to the so-called non-conventional finite element formulations. Essentially, a methodology for the enrichment of the initial approximation which is typical of the meshless methods and based on the clouds concept is introduced in the hybrid-Trefftz formulation for plane elasticity. The formulation presented allows for the approximation and direct enrichment of two independent fields: stresses in the domains and displacements on the boundaries of the elements. Defined by a set of elements and interior boundaries sharing a common node, the cloud notion is employed to select the enrichment support for the approximation fields. The numerical analysis performed reveals an excellent performance of the resulting formulation, characterized by the good approximation ability and a reduced computational effort. Copyright (C) 2009 John Wiley & Sons, Ltd.

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.

Analysis of stiffened plates by the boundary element method

In this paper, a formulation for representation of stiffeners in plane stress by the boundary elements method (BEM) in linear analysis is presented. The strategy is to adopt approximations for the displacements in the central line of the stiffener. With this simplification the Spurious oscillations in the stress along stiffeners with small thickness is prevented. Worked examples are analyzed to show the efficiency of these techniques, especially in the insertion of very narrow sub-regions, in which quasi-singular integrals are calculated, with stiffeners that are much stiffer than the main domain. The results obtained with this formulation are very close to those obtained with other formulations. (C) 2007 Elsevier Ltd. All rights reserved.

Consolidation of elastic-plastic saturated porous media by the boundary element method

This paper presents a domain boundary element formulation for inelastic saturated porous media with rate-independent behavior for the solid skeleton. The formulation is then applied to elastic-plastic behavior for the solid. Biot`s consolidation theory, extended to include irreversible phenomena is considered and the direct boundary element technique is used for the numerical solution after time discretization by the implicit Euler backward algorithm. The associated nonlinear algebraic problem is solved by the Newton-Raphson procedure whereby the loading/unloading conditions are fully taken into account and the consistent tangent operator defined. Only domain nodes (nodes defined inside the domain) are used to represent all domain values and the corresponding integrals are computed by using an accurate sub-elementation scheme. The developments are illustrated through the Drucker-Prager elastic-plastic model for the solid skeleton and various examples are analyzed with the proposed algorithms. (c) 2008 Elsevier B.V. All rights reserved.

A STUDY OF PENALTY FORMULATIONS USED IN THE NUMERICAL APPROXIMATION OF A RADIALLY SYMMETRIC ELASTICITY PROBLEM

We consider a class of two-dimensional problems in classical linear elasticity for which material overlapping occurs in the absence of singularities. Of course, material overlapping is not physically realistic, and one possible way to prevent it uses a constrained minimization theory. In this theory, a minimization problem consists of minimizing the total potential energy of a linear elastic body subject to the constraint that the deformation field must be locally invertible. Here, we use an interior and an exterior penalty formulation of the minimization problem together with both a standard finite element method and classical nonlinear programming techniques to compute the minimizers. We compare both formulations by solving a plane problem numerically in the context of the constrained minimization theory. The problem has a closed-form solution, which is used to validate the numerical results. This solution is regular everywhere, including the boundary. In particular, we show numerical results which indicate that, for a fixed finite element mesh, the sequences of numerical solutions obtained with both the interior and the exterior penalty formulations converge to the same limit function as the penalization is enforced. This limit function yields an approximate deformation field to the plane problem that is locally invertible at all points in the domain. As the mesh is refined, this field converges to the exact solution of the plane problem.

Lifecycle assessment of fuel ethanol from sugarcane in Brazil

This paper presents the lifecycle assessment (LCA) of fuel ethanol, as 100% of the vehicle fuel, from sugarcane in Brazil. The functional unit is 10,000 km run in an urban area by a car with a 1,600-cm(3) engine running on fuel hydrated ethanol, and the resulting reference flow is 1,000 kg of ethanol. The product system includes agricultural and industrial activities, distribution, cogeneration of electricity and steam, ethanol use during car driving, and industrial by-products recycling to irrigate sugarcane fields. The use of sugarcane by the ethanol agribusiness is one of the foremost financial resources for the economy of the Brazilian rural area, which occupies extensive areas and provides far-reaching potentials for renewable fuel production. But, there are environmental impacts during the fuel ethanol lifecycle, which this paper intents to analyze, including addressing the main activities responsible for such impacts and indicating some suggestions to minimize the impacts. This study is classified as an applied quantitative research, and the technical procedure to achieve the exploratory goal is based on bibliographic revision, documental research, primary data collection, and study cases at sugarcane farms and fuel ethanol industries in the northeast of SA o pound Paulo State, Brazil. The methodological structure for this LCA study is in agreement with the International Standardization Organization, and the method used is the Environmental Design of Industrial Products. The lifecycle impact assessment (LCIA) covers the following emission-related impact categories: global warming, ozone formation, acidification, nutrient enrichment, ecotoxicity, and human toxicity. The results of the fuel ethanol LCI demonstrate that even though alcohol is considered a renewable fuel because it comes from biomass (sugarcane), it uses a high quantity and diversity of nonrenewable resources over its lifecycle. The input of renewable resources is also high mainly because of the water consumption in the industrial phases, due to the sugarcane washing process. During the lifecycle of alcohol, there is a surplus of electric energy due to the cogeneration activity. Another focus point is the quantity of emissions to the atmosphere and the diversity of the substances emitted. Harvesting is the unit process that contributes most to global warming. For photochemical ozone formation, harvesting is also the activity with the strongest contributions due to the burning in harvesting and the emissions from using diesel fuel. The acidification impact potential is mostly due to the NOx emitted by the combustion of ethanol during use, on account of the sulfuric acid use in the industrial process and because of the NOx emitted by the burning in harvesting. The main consequence of the intensive use of fertilizers to the field is the high nutrient enrichment impact potential associated with this activity. The main contributions to the ecotoxicity impact potential come from chemical applications during crop growth. The activity that presents the highest impact potential for human toxicity (HT) via air and via soil is harvesting. Via water, HT potential is high in harvesting due to lubricant use on the machines. The normalization results indicate that nutrient enrichment, acidification, and human toxicity via air and via water are the most significant impact potentials for the lifecycle of fuel ethanol. The fuel ethanol lifecycle contributes negatively to all the impact potentials analyzed: global warming, ozone formation, acidification, nutrient enrichment, ecotoxicity, and human toxicity. Concerning energy consumption, it consumes less energy than its own production largely because of the electricity cogeneration system, but this process is highly dependent on water. The main causes for the biggest impact potential indicated by the normalization is the nutrient application, the burning in harvesting and the use of diesel fuel. The recommendations for the ethanol lifecycle are: harvesting the sugarcane without burning; more environmentally benign agricultural practices; renewable fuel rather than diesel; not washing sugarcane and implementing water recycling systems during the industrial processing; and improving the system of gases emissions control during the use of ethanol in cars, mainly for NOx. Other studies on the fuel ethanol from sugarcane may analyze in more details the social aspects, the biodiversity, and the land use impact.

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.

Modeling and Analysis of Piezoelectric Energy Harvesting From Aeroelastic Vibrations Using the Doublet-Lattice Method

Multifunctional structures are pointed out as an important technology for the design of aircraft with volume, mass, and energy source limitations such as unmanned air vehicles (UAVs) and micro air vehicles (MAVs). In addition to its primary function of bearing aerodynamic loads, the wing/spar structure of an UAV or a MAV with embedded piezoceramics can provide an extra electrical energy source based on the concept of vibration energy harvesting to power small and wireless electronic components. Aeroelastic vibrations of a lifting surface can be converted into electricity using piezoelectric transduction. In this paper, frequency-domain piezoaeroelastic modeling and analysis of a canti-levered platelike wing with embedded piezoceramics is presented for energy harvesting. The electromechanical finite-element plate model is based on the thin-plate (Kirchhoff) assumptions while the unsteady aerodynamic model uses the doublet-lattice method. The electromechanical and aerodynamic models are combined to obtain the piezoaeroelastic equations, which are solved using a p-k scheme that accounts for the electromechanical coupling. The evolution of the aerodynamic damping and the frequency of each mode are obtained with changing airflow speed for a given electrical circuit. Expressions for piezoaeroelastically coupled frequency response functions (voltage, current, and electrical power as well the vibratory motion) are also defined by combining flow excitation with harmonic base excitation. Hence, piezoaeroelastic evolution can be investigated in frequency domain for different airflow speeds and electrical boundary conditions. [DOI:10.1115/1.4002785]

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

No relationship found between point heteroplasmy in mitochondrial DNA control region and age range, sex and haplogroup in human hairs

The analysis of heteroplasmy (presence of more than one type of mitochondrial DNA in an individual) is used as a tool in human identification studies, anthropology, and most currently in studies that relate heteroplasmy with longevity. The frequency of heteroplasmy and its correlation with age has been analyzed using different tissues such as blood, muscle, heart, bone and brain and in different regions of mitochondrial DNA, but this analysis had never been performed using hair samples. In this study, samples of hair were sequenced in order to ascertain whether the presence or not of heteroplasmy varied according to age, sex and origin of haplogroup individuals. The samples were grouped by age (3 groups), gender (male and female) and haplogroup of origin (European, African and Native American), and analyzed using the chi-square statistical test (chi(2)). Based in statistical results obtained, we conclude that there is no relationship between heteroplasmy and sex, age and haplogroup origin using hair samples.

Multiple views of sustainable urban mobility: The case of Brazil

The aim of this work is to identify key factors of a sustainable urban mobility concept in a particular context. A multiple criteria decision analysis method was developed to identify the main variables associated to the concept. Looking at the results obtained in 11 cities of the five Brazilian regions, we conclude that the method is able to capture the different views and approaches discussed in the formulation of the mobility concept. Therefore, it can be used as a starting point for the formulation of public policies and also in the development of tools designed for monitoring the mobility conditions. (C) 2008 Elsevier Ltd. All rights reserved.

Wavelet-Galerkin method for one-dimensional elastoplasticity and damage problems: Constitutive modeling and computational aspects

This work presents an analysis of the wavelet-Galerkin method for one-dimensional elastoplastic-damage problems. Time-stepping algorithm for non-linear dynamics is presented. Numerical treatment of the constitutive models is developed by the use of return-mapping algorithm. For spacial discretization we can use wavelet-Galerkin method instead of standard finite element method. This approach allows to locate singularities. The discrete formulation developed can be applied to the simulation of one-dimensional problems for elastic-plastic-damage models. (C) 2007 Elsevier Inc. All rights reserved.