Hybrid Method to Assess Sensitive Process Interruption Costs Due to Faults in Electric Power Distribution Networks

This paper shows a new hybrid method for risk assessment regarding interruptions in sensitive processes due to faults in electric power distribution systems. This method determines indices related to long duration interruptions and short duration voltage variations (SDVV), such as voltage sags and swells in each customer supplied by the distribution network. Frequency of such occurrences and their impact on customer processes are determined for each bus and classified according to their corresponding magnitude and duration. The method is based on information regarding network configuration, system parameters and protective devices. It randomly generates a number of fault scenarios in order to assess risk areas regarding long duration interruptions and voltage sags and swells in an especially inventive way, including frequency of events according to their magnitude and duration. Based on sensitivity curves, the method determines frequency indices regarding disruption in customer processes that represent equipment malfunction and possible process interruptions due to voltage sags and swells. Such approach allows for the assessment of the annual costs associated with each one of the evaluated power quality indices.

An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 1: Rods

A fully conserving algorithm is developed in this paper for the integration of the equations of motion in nonlinear rod dynamics. The starting point is a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, which results in an extremely simple update of the rotational variables. The weak form is constructed with a non-orthogonal projection corresponding to the application of the virtual power theorem. Together with an appropriate time-collocation, it ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that nonlinear hyperelastic materials (and not only materials with quadratic potentials) are permitted without any prejudice on the conservation properties. Spatial discretization is performed via the finite element method and the performance of the scheme is assessed by means of several numerical simulations.

Hybrid and multi-field variational principles for geometrically exact three-dimensional beams

This paper addresses the development of several alternative novel hybrid/multi-field variational formulations of the geometrically exact three-dimensional elastostatic beam boundary-value problem. In the framework of the complementary energy-based formulations, a Legendre transformation is used to introduce the complementary energy density in the variational statements as a function of stresses only. The corresponding variational principles are shown to feature stationarity within the framework of the boundary-value problem. Both weak and linearized weak forms of the principles are presented. The main features of the principles are highlighted, giving special emphasis to their relationships from both theoretical and computational standpoints. (C) 2010 Elsevier Ltd. All rights reserved.

Effects of elastic indenter deformation on spherical instrumented indentation tests: the reduced elastic modulus

Although the Hertz theory is not applicable in the analysis of the indentation of elastic-plastic materials, it is common practice to incorporate the concept of indenter/specimen combined modulus to consider indenter deformation. The appropriateness was assessed of the use of reduced modulus to incorporate the effect of indenter deformation in the analysis of the indentation with spherical indenters. The analysis based on finite element simulations considered four values of the ratio of the indented material elastic modulus to that of the diamond indenter, E/E(i) (0, 0.04, 0.19, 0.39), four values of the ratio of the elastic reduced modulus to the initial yield strength, E(r)/Y (0, 10, 20, 100), and two values of the ratio of the indenter radius to maximum total displacement, R/delta(max) (3, 10). Indenter deformation effects are better accounted for by the reduced modulus if the indented material behaves entirely elastically. In this case, identical load-displacement (P - delta) curves are obtained with rigid and elastic spherical indenters for the same elastic reduced modulus. Changes in the ratio E/E(i), from 0 to 0.39, resulted in variations lower than 5% for the load dimensionless functions, lower than 3% in the contact area, A(c), and lower than 5% in the ratio H/E(r). However, deformations of the elastic indenter made the actual radius of contact change, even in the indentation of elastic materials. Even though the load dimensionless functions showed only a little increase with the ratio E/E(i), the hardening coefficient and the yield strength could be slightly overestimated when algorithms based on rigid indenters are used. For the unloading curves, the ratio delta(e)/delta(max), where delta(e) is the point corresponding to zero load of a straight line with slope S from the point (P(max), delta(max)), varied less than 5% with the ratio E/E(i). Similarly, the relationship between reduced modulus and the unloading indentation curve, expressed by Sneddon`s equation, did not reveal the necessity of correction with the ratio E/E(i). The most affected parameter in the indentation curve, as a consequence of the indentation deformation, was the ratio between the residual indentation depth after complete unloading and the maximum indenter displacement, delta(r)/delta(max) (up to 26%), but this variation did not significantly decrease the capability to estimate hardness and elastic modulus based on the ratio of the residual indentation depth to maximum indentation depth, h(r)/h(max). In general, the results confirm the convenience of the use of the reduced modulus in the spherical instrumented indentation tests.

Modeling and simulation of severe slugging in air-water pipeline-riser systems

A mathematical model, numerical simulations and stability and flow regime maps corresponding to severe slugging in pipeline riser systems, are presented. In the simulations air and water were used as flowing fluids. The mathematical model considers continuity equations for liquid and gas phases, with a simplified momentum equation for the mixture, neglecting inertia. A drift-flux model, evaluated for the local conditions in the riser, is used as a closure law. The developed model predicts the location of the liquid accumulation front in the pipeline and the liquid level in the riser, so it is possible to determine which type of severe slugging occurs in the system. The numerical procedure is convergent for different nodalizations. A comparison is made with experimental results corresponding to a catenary riser, showing very good results for slugging cycle and stability and flow regime maps. (c) 2010 Elsevier Ltd. All rights reserved.

Tempering temperature effects on abrasive wear of mottled cast iron

The effects of different tempering temperatures (300-600 degrees C) on abrasive wear resistance of mottled cast iron were studied. Abrasive wear tests were carried out using the rubber-wheel test on quartz sand and the pin test on Al(2)O(3) abrasive cloths. The retained austenite content of the matrix was determined by X-ray diffraction. The wear surface of the specimens was examined by scanning electron microscopy for identifying the wear micromechanism. Bulk hardness and matrix hardness before and after the tests were measured. The results showed that in the two-body (pin-on-disc test) system, the main wear mechanism was microcutting and high matrix hardening was presented. The wear rates presented higher correlation with the retained austenite than with the bulk and matrix hardness. In the three-body system (sand-rubber wheel), the wear surfaces presented indentations due to abrasive rolling. The wear rates had better correlation with both the bulk and matrix hardness (before and after the wear test) than with the retained austenite content. There are two groups of results, high and low wear rates corresponding to each tribosystem, two-body abrasive wear and three-body abrasive wear, respectively. (C) 2009 Elsevier B.V. All rights reserved.

Characterization of the elastic-plastic region in AISI/SAE 1070 steel by the magnetic barkhausen noise

The magnetic Barkhausen energy in the rolling and transversal directions of AISI/SAE 1070 annealed surfaces is studied. The measurements were made in the samples under applied tension in the elastic-plastic region for different angular directions. The outcomes evidence that the magnetic anisotropy coefficient can be used to characterize the linear and nonlinear elastic limits of the material tinder tensile tresses. The results also show that the area of the curve corresponding to the angular dependence of the number of Barkhausen jumps with average energy presents a maximum value that corresponds to the elastic limit of the sample. (C) 2008 Elsevier Ltd. All rights reserved.

A procedure to evaluate environmental rehabilitation in limestone quarries

A procedure to evaluate mine rehabilitation practices during the operational phase was developed and validated. It is based on a comparison of actually observed or documented practices with internationally recommended best practices (BP). A set of 150 BP statements was derived from international guides in order to establish the benchmark. The statements are arranged in six rehabilitation programs under three categories: (1) planning (2) operational and (3) management, corresponding to the adoption of the plan-do-check-act management systems model to mine rehabilitation. The procedure consists of (i) performing technical inspections guided by a series of field forms containing BP statements; (ii) classifying evidences in five categories; and (iii) calculating conformity indexes and levels. For testing and calibration purposes, the procedure was applied to nine limestone quarries and conformity indexes were calculated for the rehabilitation programs in each quarry. Most quarries featured poor planning practices, operational practices reached high conformity levels in 50% of the cases and management practices scored moderate conformity. Despite all quarries being ISO 14001 certified, their management systems pay low attention to issues pertaining to land rehabilitation and biodiversity. The best results were achieved by a quarry whose expansion was recently submitted to the environmental impact assessment process, suggesting that public scrutiny may play a positive role in enhancing rehabilitation practices. Conformity indexes and levels can be used to chart the evolution of rehabilitation practices at regular intervals, to establish corporate goals and for communication with stakeholders. (C) 2010 Elsevier Ltd. All rights reserved.

Dynamic Design of Piezoelectric Laminated Sensors and Actuators using Topology Optimization

Sensors and actuators based on piezoelectric plates have shown increasing demand in the field of smart structures, including the development of actuators for cooling and fluid-pumping applications and transducers for novel energy-harvesting devices. This project involves the development of a topology optimization formulation for dynamic design of piezoelectric laminated plates aiming at piezoelectric sensors, actuators and energy-harvesting applications. It distributes piezoelectric material over a metallic plate in order to achieve a desired dynamic behavior with specified resonance frequencies, modes, and enhanced electromechanical coupling factor (EMCC). The finite element employs a piezoelectric plate based on the MITC formulation, which is reliable, efficient and avoids the shear locking problem. The topology optimization formulation is based on the PEMAP-P model combined with the RAMP model, where the design variables are the pseudo-densities that describe the amount of piezoelectric material at each finite element and its polarization sign. The design problem formulated aims at designing simultaneously an eigenshape, i.e., maximizing and minimizing vibration amplitudes at certain points of the structure in a given eigenmode, while tuning the eigenvalue to a desired value and also maximizing its EMCC, so that the energy conversion is maximized for that mode. The optimization problem is solved by using sequential linear programming. Through this formulation, a design with enhancing energy conversion in the low-frequency spectrum is obtained, by minimizing a set of first eigenvalues, enhancing their corresponding eigenshapes while maximizing their EMCCs, which can be considered an approach to the design of energy-harvesting devices. The implementation of the topology optimization algorithm and some results are presented to illustrate the method.

Design and Characterization of a Biomimetic Piezoelectric Pump Inspired on Group Fish Swimming Effect

Flow pumps are important tools in several engineering areas, such as in the fields of bioengineering and thermal management solutions for electronic devices. Nowadays, many of the new flow pump principles are based on the use of piezoelectric actuators, which present some advantages such as miniaturization potential and lower noise generation. In previous work, authors presented a study of a novel pump configuration based on placing an oscillating bimorph piezoelectric actuator in water to generate flow. It was concluded that this oscillatory behavior (such as fish swimming) yields vortex interaction, generating flow rate due to the action and reaction principle. Thus, following this idea the objective of this work is to explore this oscillatory principle by studying the interaction among generated vortex from two bimorph piezoelectric actuators oscillating inside the same pump channel, which is similar to the interaction of vortex generated by frontal fish and posterior ones when they swim together in a group formation. It is shown that parallel-series configurations of bimorph piezoelectric actuators inside the same pump channel provide higher flow rates and pressure for liquid pumping than simple parallel-series arrangements of corresponding single piezoelectric pumps, respectively. The scope of this work includes structural simulations of bimorph piezoelectric actuators, fluid flow simulations, and prototype construction for result validation.

Recycling Krylov subspaces for efficient large-scale electrical impedance tomography

Electrical impedance tomography (EIT) captures images of internal features of a body. Electrodes are attached to the boundary of the body, low intensity alternating currents are applied, and the resulting electric potentials are measured. Then, based on the measurements, an estimation algorithm obtains the three-dimensional internal admittivity distribution that corresponds to the image. One of the main goals of medical EIT is to achieve high resolution and an accurate result at low computational cost. However, when the finite element method (FEM) is employed and the corresponding mesh is refined to increase resolution and accuracy, the computational cost increases substantially, especially in the estimation of absolute admittivity distributions. Therefore, we consider in this work a fast iterative solver for the forward problem, which was previously reported in the context of structural optimization. We propose several improvements to this solver to increase its performance in the EIT context. The solver is based on the recycling of approximate invariant subspaces, and it is applied to reduce the EIT computation time for a constant and high resolution finite element mesh. In addition, we consider a powerful preconditioner and provide a detailed pseudocode for the improved iterative solver. The numerical results show the effectiveness of our approach: the proposed algorithm is faster than the preconditioned conjugate gradient (CG) algorithm. The results also show that even on a standard PC without parallelization, a high mesh resolution (more than 150,000 degrees of freedom) can be used for image estimation at a relatively low computational cost. (C) 2010 Elsevier B.V. All rights reserved.

Piezoresistive sensor design using topology optimization

Piezoresistive materials, materials whose resistivity properties change when subjected to mechanical stresses, are widely utilized in many industries as sensors, including pressure sensors, accelerometers, inclinometers, and load cells. Basic piezoresistive sensors consist of piezoresistive devices bonded to a flexible structure, such as a cantilever or a membrane, where the flexible structure transmits pressure, force, or inertial force due to acceleration, thereby causing a stress that changes the resistivity of the piezoresistive devices. By applying a voltage to a piezoresistive device, its resistivity can be measured and correlated with the amplitude of an applied pressure or force. The performance of a piezoresistive sensor is closely related to the design of its flexible structure. In this research, we propose a generic topology optimization formulation for the design of piezoresistive sensors where the primary aim is high response. First, the concept of topology optimization is briefly discussed. Next, design requirements are clarified, and corresponding objective functions and the optimization problem are formulated. An optimization algorithm is constructed based on these formulations. Finally, several design examples of piezoresistive sensors are presented to confirm the usefulness of the proposed method.

A ductile damage criterion at various stress triaxialities

The paper discusses the effect of stress triaxiality on the onset and evolution of damage in ductile metals. A series of tests including shear tests and experiments oil smooth and pre-notched tension specimens wits carried Out for it wide range of stress triaxialities. The underlying continuum damage model is based oil kinematic definition of damage tensors. The modular structure of the approach is accomplished by the decomposition of strain rates into elastic, plastic and damage parts. Free energy functions with respect to fictitious undamaged configurations as well as damaged ones are introduced separately leading to elastic material laws which are affected by increasing damage. In addition, a macroscopic yield condition and a flow rule are used to adequately describe the plastic behavior. Numerical simulations of the experiments are performed and good correlation of tests and numerical results is achieved. Based oil experimental and numerical data the damage criterion formulated in stress space is quantified. Different branches of this function are taken into account corresponding to different damage modes depending oil stress triaxiality and Lode parameter. In addition, identification of material parameters is discussed ill detail. (C) 2007 Elsevier Ltd. All rights reserved.

Emissions from the premixed combustion of gasified polyethylene

An investigation was conducted on pollutants emitted from steady-state, steady-flow gasification and combustion of polyethylene (PE) in a two-stage furnace. The polymer, in pulverized form, was first pyrolyzed at 1000 degrees C, and subsequently, its gaseous pyrolyzates were burned, upon mixing with air at high temperatures (900-1100 degrees C). The motivation for this indirect type of burning PE was to attain nominally premixed combustion of the pyrolyzate gases with air, thereby achieving lower pollutant emissions than those emanating from the direct burning of the solid PE polymer. This work assessed the effluents of the two-stage furnace and examined the effects of the combustion temperature, as well as the polymer feed rate and the associated fuel/air equivalence ratio (0.3 < phi < 1.4). It was found that, whereas the yield of pyrolysis gas decreased with an increasing polymer feed rate, its composition was nearly independent of the feed rate. CO2 emissions peaked at an equivalence ratio near unity, while the CO emissions increased with an increasing equivalence ratio. The total light volatile hydrocarbon and semivolatile polycyclic aromatic hydrocarbon (PAH) emissions of combustion increased with an increasing equivalence ratio. The generated particulates were mostly submicrometer in size. Overall, PAH and soot emissions from this indirect burning of PE were an order of magnitude lower than corresponding emissions from the direct burning of the solid polymer, obtained previously in this laboratory using identical sampling and analytical techniques. Because pyrolysis of this polymer requires a nominal heat input that amounts to only a diminutive fraction of the heat released during its combustion, implementation of this technique is deemed advantageous.

Influence of Microstructures on Polarization Resistance and Hydrogen-Induced Cracking of a Micro-Alloyed Steel Submitted to Different Cooling Rates

The effect of different microstructures on the polarization resistance (Rp) and the hydrogen-induced cracking (HIC) of a micro-alloyed steel austenitized and submitted to different cooling rates was studied. Samples 19.1 x 6 x 2 mm, containing the whole thickness of the plate were extracted from a 20 mm plate and heat treated on a quenching dilatometer, were submitted to Rp and HIC corrosion tests. Both Rp and HIC tests followed as close as possible ASTM G59 and NACE standard TM0284-2003, in this case, modified only with regard to the size of the samples. Steel samples transformed from austenite by a slow cooling (cooling rate of 0.5 degrees C.s(-1)) showed higher susceptibility to hydrogen-induced cracking, with large cracks in the middle of the sample propagating along segregation bands, corresponding to the centerline of the plate thickness. For cooling rates of 10 degrees C.s(-1), only small cracks were found in the matrix and micro cracks nucleated at non-metallic inclusions. For higher cooling rates (40 degrees C.s(-1)) very few small cracks were detected, linked to non-metallic inclusions. This result suggests that structures formed by polygonal structures and segregation bands (were cutectoid microconstituents predominate) have higher susceptibility to HIC. Structures predominantly formed by acicular ferrite make it difficult to propagate the cracks among non-oriented and interlaced acicular ferrite crystals. Smaller segregation bands containing eutectoid products also help inhibit cracking and crack propagation; segregation bands can function as pipelines for hydrogen diffusion and offer a path of stress concentration for the propagation of cracks, frequently associated to non-metallic inclusions. Polarization resistance essays performed on the steel in theas received condition, prior to any heat treatment, showed larger differences between the regions of the plate, with a considerably lower Rp in the centerline. The austenitization heat treatments followed by cooling rates of 0.5 e 10 degrees C.s(-1) made more uniform the corrosion resistance along the thickness of the plate. The effects of heat treatments on the corrosion resistance are probably related to the microconstituent formed, allied to the chemical homogenization of the impurities concentrated on the centerline of the plate.