Analyzing static three-dimensional elastic domains with a new infinite boundary element formulation

A new two-dimensionally mapped infinite boundary element (IBE) is presented. The formulation is based on a triangular boundary element (BE) with linear shape functions instead of the quadrilateral IBEs usually found in the literature. The infinite solids analyzed are assumed to be three-dimensional, linear-elastic and isotropic, and Kelvin fundamental solutions are employed. One advantage of the proposed formulation over quadratic or higher order elements is that no additional degrees of freedom are added to the original BE mesh by the presence of the IBEs. Thus, the IBEs allow the mesh to be reduced without compromising the accuracy of the result. Two examples are presented, in which the numerical results show good agreement with authors using quadrilateral IBEs and analytical solutions. (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.

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.

Non-linear boundary element formulation with tangent operator to analyse crack propagation in quasi-brittle materials

This work deals with analysis of cracked structures using BEM. Two formulations to analyse the crack growth process in quasi-brittle materials are discussed. They are based on the dual formulation of BEM where two different integral equations are employed along the opposite sides of the crack surface. The first presented formulation uses the concept of constant operator, in which the corrections of the nonlinear process are made only by applying appropriate tractions along the crack surfaces. The second presented BEM formulation to analyse crack growth problems is an implicit technique based on the use of a consistent tangent operator. This formulation is accurate, stable and always requires much less iterations to reach the equilibrium within a given load increment in comparison with the classical approach. Comparison examples of classical problem of crack growth are shown to illustrate the performance of the two formulations. (C) 2009 Elsevier Ltd. All rights reserved.

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.

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.

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.

Effect of impeller type and stirring frequency on the behavior of an AnSBBR in the treatment of low-strength wastewater

The influence of impeller type and stirring frequency on the performance of a mechanically stirred anaerobic sequencing batch reactor containing immobilized biomass on an inert support (AnSBBR - Anaerobic Sequencing Batch Biofilm Reactor) was evaluated. The biomass was immobilized on polyurethane foam cubes placed in a stainless-steel basket inside a glass cylinder. Each 8-h batch run consisted of three stages: feed (10 min), reaction (460 min) and discharge (10 min) at 30 degrees C. Experiments were performed with four impeller types, i.e., helical, flat-blade, inclined-blade and curved-blade turbines, at stirring frequencies ranging from 100 to 1100 rpm. Synthetic wastewater was used in all experiments with an organic-matter concentration of 530 +/- 37 mg/L measured as chemical oxygen demand (COD). The reactor achieved an organic-matter removal efficiency of around 87% under all investigated conditions. Analysis of the four impeller types and the investigated stirring frequencies showed that mass transfer in the liquid phase was affected not only by the applied stirring frequency but also by the agitation mode imposed by each impeller type. The best reactor performance at all stirring frequencies was obtained when agitation was provided by the flat-blade turbine impeller. (C) 2010 Elsevier Ltd. All rights reserved.

Estimation of Mass Transfer Velocity Based on Measured Turbulence Parameters

The aim of this study is to quantify the mass transfer velocity using turbulence parameters from simultaneous measurements of oxygen concentration fields and velocity fields. The surface divergence model was considered in more detail, using data obtained for the lower range of beta (surface divergence). It is shown that the existing models that use the divergence concept furnish good predictions for the transfer velocity also for low values of beta, in the range of this study. Additionally, traditional conceptual models, such as the film model, the penetration-renewal model, and the large eddy model, were tested using the simultaneous information of concentration and velocity fields. It is shown that the film and the surface divergence models predicted the mass transfer velocity for all the range of the equipment Reynolds number used here. The velocity measurements showed viscosity effects close to the surface, which indicates that the surface was contaminated with some surfactant. Considering the results, this contamination can be considered slight for the mass transfer predictions. (C) 2009 American Institute of Chemical Engineers AIChE J, 56: 2005-2017; 2010

Enhanced aeroelastic energy harvesting by exploiting combined nonlinearities: theory and experiment

Converting aeroelastic vibrations into electricity for low power generation has received growing attention over the past few years. In addition to potential applications for aerospace structures, the goal is to develop alternative and scalable configurations for wind energy harvesting to use in wireless electronic systems. This paper presents modeling and experiments of aeroelastic energy harvesting using piezoelectric transduction with a focus on exploiting combined nonlinearities. An airfoil with plunge and pitch degrees of freedom (DOF) is investigated. Piezoelectric coupling is introduced to the plunge DOF while nonlinearities are introduced through the pitch DOF. A state-space model is presented and employed for the simulations of the piezoaeroelastic generator. A two-state approximation to Theodorsen aerodynamics is used in order to determine the unsteady aerodynamic loads. Three case studies are presented. First the interaction between piezoelectric power generation and linear aeroelastic behavior of a typical section is investigated for a set of resistive loads. Model predictions are compared to experimental data obtained from the wind tunnel tests at the flutter boundary. In the second case study, free play nonlinearity is added to the pitch DOF and it is shown that nonlinear limit-cycle oscillations can be obtained not only above but also below the linear flutter speed. The experimental results are successfully predicted by the model simulations. Finally, the combination of cubic hardening stiffness and free play nonlinearities is considered in the pitch DOF. The nonlinear piezoaeroelastic response is investigated for different values of the nonlinear-to-linear stiffness ratio. The free play nonlinearity reduces the cut-in speed while the hardening stiffness helps in obtaining persistent oscillations of acceptable amplitude over a wider range of airflow speeds. Such nonlinearities can be introduced to aeroelastic energy harvesters (exploiting piezoelectric or other transduction mechanisms) for performance enhancement.

Piezoceramic Composite Actuators for Flow Control in Low Reynolds Number Airflow

The research presented here employs solid-state actuators for flow separation delay or for forced attachment of separated flow seen in airfoils at low Reynolds numbers. To reduce separation, periodic excitation to the flow around the leading edge of the airfoil is induced by Macro-Fiber Composite actuated clamped-free unimorph benders. An electromechanical model of the unimorph is briefly presented and parametric study is conducted to aid the design of a unimorph to output high deformation at a desired frequency. The optimum frequency and amplitude for lift improvement at post-stall angles are identified experimentally. Along with aerodynamic force and structural displacement measurements, helium bubble flow visualization is used to verify existing separated flow, and the attached flow induced by flow control. The lift enhancement induced by several flow control techniques is compared. A symmetric and non-uniform (3D) flow excitation results in the maximum lift enhancement at post-stall region at the lowest power consumption level. A maximum lift coefficient increase of 27.5% (in the post-stall region) is achieved at 125 Hz periodic excitation, with the 3D symmetric actuation mode at 5 m/s and the reduced frequency of 3.78. C(l,max) is increased 7.6% from the baseline.

Is Nitric Oxide a Mediator of the Effects of Low-Intensity Electrical Stimulation on Bone in Ovariectomized Rats?

Low-intensity electrical stimulation (LIES) may counteract the effects of ovariectomy (OVX) on nitric oxide synthase (NOS) expression, osteocyte viability, bone structure, and microarchitecture in rats (Lirani-Galvo et al., Calcif Tissue Int 84:502-509, 2009). The aim of the present study was to investigate if these effects of LIES could be mediated by NO. We analyzed the effects of NO blockage (by l-NAME) in the response to LIES on osteocyte viability, bone structure, and microarchitecture in OVX rats. Sixty rats (200-220 g) were divided into six groups: sham, sham-l-NAME (6 mg/kg/day), OVX, OVX-l-NAME, OVX-LIES, and OVX-LIES-l-NAME. After 12 weeks, rats were killed and tibiae collected for histomorphometric analysis and immunohistochemical detection of endothelial NOS (eNOS), inducible NOS (iNOS), and osteocyte apoptosis (caspase-3 and TUNEL). In the presence of l-NAME, LIES did not counteract the OVX-induced effects on bone volume and trabecular number (as on OVX-LIES). l-NAME blocked the stimulatory effects of LIES on iNOS and eNOS expression of OVX rats. Both l-NAME and LIES decreased osteocyte apoptosis. Our results showed that in OVX rats l-NAME partially blocks the effects of LIES on bone structure, turnover, and expression of iNOS and eNOS, suggesting that NO may be a mediator of some positive effects of LIES on bone.

Low-Intensity Electrical Stimulation Counteracts the Effects of Ovariectomy on Bone Tissue of Rats: Effects on Bone Microarchitecture, Viability of Osteocytes, and Nitric Oxide Expression

Low Intensity Electrical Stimulation (LIES) has been used for bone repair, but little is known about its effects on bone after menopause. Osteocytes probably play a role in mediating this physical stimulus and they could act as transducers through the release of biochemical signals, such as nitric oxide (NO). The aim of the present study was to investigate the effects of LIES on bone structure and remodeling, NOS expression and osteocyte viability in ovariectomized (OVX) rats. Thirty rats (200-220 g) were divided into 3 groups: SHAM, OVX, and OVX subjected to LIES (OVX + LIES) for 12 weeks. Following the protocol, rats were sacrificed and tibias were collected for histomorphometric analysis and immunohistochemical detection of endothelial NO synthase (eNOS), inducible NOS (iNOS), and osteocyte apoptosis (caspase-3 and TUNEL). OVX rats showed significant (p < 0.05 vs. SHAM) decreased bone volume (10% vs. 25%) and trabecular number (1.7 vs. 3.9), and increased eroded surfaces (4.7% vs. 3.2%) and mineralization surfaces (15.9% vs. 7.7%). In contrast, after LIES, all these parameters were significantly different from OVX but not different from SHAM. eNOS and iNOS were similarly expressed in subperiosteal regions of tibiae cortices of SHAM, not expressed in OVX, and similarly expressed in OVX + LIES when compared to SHAM. In OVX, the percentage of apoptotic osteocytes (24%) was significantly increased when compared to SHAM (11%) and OVX + LIES (8%). Our results suggest that LIES counteracts some effects of OVX on bone tissue preserving bone structure and microarchitecture, iNOS and eNOS expression, and osteocyte viability.

Acoustic Simulation of a Special Switched Reluctance Drive by Means of Field-Circuit Coupling and Multiphysics Simulation

The approach presented in this paper consists of an energy-based field-circuit coupling in combination with multi-physics simulation of the acoustic radiation of electrical machines. The proposed method is applied to a special switched reluctance motor with asymmetric pole geometry to improve the start-up torque. The pole shape has been optimized, subject to low torque ripple, in a previous study. The proposed approach here is used to analyze the impact of the optimization on the overall acoustic behavior. The field-circuit coupling is based on a temporary lumped-parameter model of the magnetic part incorporated into a circuit simulation based on the modified nodal analysis. The harmonic force excitation is calculated by means of stress tensor computation, and it is transformed to a mechanical mesh by mapping techniques. The structural dynamic problem is solved in the frequency domain using a finite-element modal analysis and superposition. The radiation characteristic is obtained from boundary element acoustic simulation. Simulation results of both rotor types are compared, and measurements of the drive are presented.