Design methodology of piezoelectric energy-harvesting skin using topology optimization

Abstract This paper describes a design methodology for piezoelectric energy harvester s that thinly encapsulate the mechanical devices and expl oit resonances from higher- order vibrational modes. The direction of polarization determines the sign of the pi ezoelectric tensor to avoid cancellations of electric fields from opposite polarizations in the same circuit. The resultant modified equations of state are solved by finite element method (FEM). Com- bining this method with the solid isotropic material with penalization (SIMP) method for piezoelectric material, we have developed an optimization methodology that optimizes the piezoelectric material layout and polarization direc- tion. Updating the density function of the SIMP method is performed based on sensitivity analysis, the sequen- tial linear programming on the early stage of the opti- mization, and the phase field method on the latter stage

Analysis of elastic functionally graded materials under large displacements via high-order tetrahedral elements

The purpose of this study is to present a position based tetrahedral finite element method of any order to accurately predict the mechanical behavior of solids constituted by functionally graded elastic materials and subjected to large displacements. The application of high-order elements makes it possible to overcome the volumetric and shear locking that appears in usual homogeneous isotropic situations or even in non-homogeneous cases developing small or large displacements. The use of parallel processing to improve the computational efficiency, allows employing high-order elements instead of low-order ones with reduced integration techniques or strain enhancements. The Green-Lagrange strain is adopted and the constitutive relation is the functionally graded Saint Venant-Kirchhoff law. The equilibrium is achieved by the minimum total potential energy principle. Examples of large displacement problems are presented and results confirm the locking free behavior of high-order elements for non-homogeneous materials. (C) 2011 Elsevier B.V. All rights reserved.

Performance of polymeric reinforcements in vehicle structures submitted to frontal impact

Preformed structural reinforcements have shown good performance in crash tests, where the great advantage is their weight. These reinforcements are designed with the aim of increasing the rigidity of regions with large deformations, thus stabilising sections of the vehicle that work as load path during impact. The objective of this work is to show the application of structural reinforcements made of polymeric material PA66 in the field of vehicle safety, through finite element simulations. Simulations of frontal impact at 50 km/h and in ODB (offset deformable barrier) at 57 km/h configurations (standards such as ECE R-94 and ECE R-12) were performed in the software LS-DYNA R (R) and MADYMO (R). The simulations showed that the use of polymeric reinforcements leads to a 70% reduction in A-pillar intrusion, a 65% reduction in the displacement of the steering column and a 59% reduction in the deformation in the region of the occupant legs and feet. The level of occupant injuries was analysed by MADYMO (R) software, and a reduction of 23.5% in the chest compression and 80% in the tibia compression were verified. According to the standard, such conditions lead to an improvement in the occupant safety in a vehicle collision event.

Abrasividade pendular e a resistência mecânica das rochas

Esse trabalho emprega o método para avaliar a abrasividade proposto por Golovanevskiy e Bearman (2008). Esse método, ensaio de abrasão por impacto deslizante (Gouging Abrasion Test), é realizado em condições de alta tensão/alto impacto de desgaste. O método consiste de uma ponteira cilíndrica com uma ponta cônica de 90º, que, em trajetória pendular, atinge uma amostra de rocha com energia de impacto de 300 J e velocidade da ordem de 5,2 m/s. O Gouging Abrasion Index (Gi) é calculado como sendo a média do diâmetro da ponta cônica, após desgaste, em milímetros e o resultado é multiplicado por 10. Esse trabalho verificou a adequabilidade do Gouging Abrasion Test, para um pequeno número de amostras de rocha, que representam, qualitativamente, os principais tipos de rocha encontrados em trabalhos de corte, perfuração e britagem no Brasil, e a sua correlação com outros ensaios consagrados como a resistência à compressão, o desgaste Amsler e a dureza Knoop. Essa análise mostrou alta correlação entre Gi e a dureza Knoop (R² = 0,94), baixa correlação com o desgaste

Making the case for grid-connected photovoltaics in Brazil

In the developed world, grid-connected photovoltaics (PVs) are the fastest-growing segment of the energy market. From 1999 to 2009, this industry had a 42% compound annual growth-rate. From 2009 to 2013, it is expected to grow to 45%, and in 2013 the achievement of grid parity - when the cost of solar electricity becomes competitive with conventional retail (including taxes and charges) grid-supplied electricity - is expected in many places worldwide. Grid-connected PV is usually perceived as an energy technology for developed countries, whereas isolated, stand-alone PV is considered as more suited for applications in developing nations, where so many individuals still lack access to electricity. This rationale is based on the still high costs of PV when compared with conventional electricity. We make the case for grid-connected PV generation in Brazil, showing that with the declining costs of PV and the rising prices of conventional electricity, urban populations in Brazil will also enjoy grid parity in the present decade. We argue that governments in developing nations should act promptly and establish the mandates and necessary conditions for their energy industry to accumulate experience in grid-connected PV, and make the most of this benign technology in the near future. (C) 2010 Elsevier Ltd. All rights reserved.

A Comparative Wet Collapse Buckling Study for the Carcass Layer of Flexible Pipes

When there is a failure on the external sheath of a flexible pipe, a high value of hydrostatic pressure is transferred to its internal plastic layer and consequently to its interlocked carcass, leading to the possibility of collapse. The design of a flexible pipe must predict the maximum value of external pressure the carcass layer can be subjected to without collapse. This value depends on the initial ovalization due to manufacturing tolerances. To study that problem, two numerical finite element models were developed to simulate the behavior of the carcass subjected to external pressure, including the plastic behavior of the materials. The first one is a full 3D model and the second one is a 3D ring model, both composed by solid elements. An interesting conclusion is that both the models provide the same results. An analytical model using an equivalent thickness approach for the carcass layer was also constructed. A good correlation between analytical and numerical models was achieved for pre-collapse behavior but the collapse pressure value and post-collapse behavior were not well predicted by the analytical model. [DOI: 10.1115/1.4005185]

Non-linear magnetohydrodynamic simulations of density evolution in Tore Supra sawtoothing plasmas

The plasma density evolution in sawtooth regime on the Tore Supra tokamak is analyzed. The density is measured using fast-sweeping X-mode reflectometry which allows tomographic reconstructions. There is evidence that density is governed by the perpendicular electric flows, while temperature evolution is dominated by parallel diffusion. Postcursor oscillations sometimes lead to the formation of a density plateau, which is explained in terms of convection cells associated with the kink mode. A crescent-shaped density structure located inside q = 1 is often visible just after the crash and indicates that some part of the density withstands the crash. 3D full MHD nonlinear simulations with the code XTOR-2F recover this structure and show that it arises from the perpendicular flows emerging from the reconnection layer. The proportion of density reinjected inside the q = 1 surface is determined, and the implications in terms of helium ash transport are discussed. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4766893]

Numerical simulation of drop impact and jet buckling problems using the eXtended Pom-Pom model

This work presents numerical simulations of two fluid flow problems involving moving free surfaces: the impacting drop and fluid jet buckling. The viscoelastic model used in these simulations is the eXtended Pom-Pom (XPP) model. To validate the code, numerical predictions of the drop impact problem for Newtonian and Oldroyd-B fluids are presented and compared with other methods. In particular, a benchmark on numerical simulations for a XPP drop impacting on a rigid plate is performed for a wide range of the relevant parameters. Finally, to provide an additional application of free surface flows of XPP fluids, the viscous jet buckling problem is simulated and discussed. (C) 2011 Elsevier B.V. All rights reserved.

Influence of pattern gradation on the design of piezocomposite energy harvesting devices using topology optimization

Piezoelectric materials can be used to convert oscillatory mechanical energy into electrical energy. Energy harvesting devices are designed to capture the ambient energy surrounding the electronics and convert it into usable electrical energy. The design of energy harvesting devices is not obvious, requiring optimization procedures. This paper investigates the influence of pattern gradation using topology optimization on the design of piezocomposite energy harvesting devices based on bending behavior. The objective function consists of maximizing the electric power generated in a load resistor. A projection scheme is employed to compute the element densities from design variables and control the length scale of the material density. Examples of two-dimensional piezocomposite energy harvesting devices are presented and discussed using the proposed method. The numerical results illustrate that pattern gradation constraints help to increase the electric power generated in a load resistor and guides the problem toward a more stable solution. (C) 2012 Elsevier Ltd. All rights reserved.

Theoretical models to predict the mechanical behavior of thick composite tubes

This paper shows theoretical models (analytical formulations) to predict the mechanical behavior of thick composite tubes and how some parameters can influence this behavior. Thus, firstly, it was developed the analytical formulations for a pressurized tube made of composite material with a single thick ply and only one lamination angle. For this case, the stress distribution and the displacement fields are investigated as function of different lamination angles and reinforcement volume fractions. The results obtained by the theoretical model are physic consistent and coherent with the literature information. After that, the previous formulations are extended in order to predict the mechanical behavior of a thick laminated tube. Both analytical formulations are implemented as a computational tool via Matlab code. The results obtained by the computational tool are compared to the finite element analyses, and the stress distribution is considered coherent. Moreover, the engineering computational tool is used to perform failure analysis, using different types of failure criteria, which identifies the damaged ply and the mode of failure.

Effect of parametric uncertainties on the performance of a piezoelectric energy harvesting device

The use of piezoelectric materials for the development of electromechanical devices for the harvesting or scavenging of ambient vibrations has been extensively studied over the last decade. The energy conversion from mechanical (vibratory) to electrical energy is provided by the electromechanical coupling between mechanical strains/stresses and electric charges/voltages in the piezoelectric material. The majority of the studies found in the open literature present a tip-mass cantilever piezoelectric device tuned on the operating frequency. Although recent results show that these devices can be quite effective for harvesting small amounts of electrical energy, little has been published on the robustness of these devices or on the effect of parametric uncertainties on the energy harvested. This work focuses on a cantilever plate with bonded piezoelectric patches and a tip-mass serving as an energy harvesting device. The rectifier and storage electric circuit was replaced by a resistive circuit (R). In addition, an alternative to improve the harvesting performance by adding an inductance in series to the harvesting circuit, thus leading to a resonant circuit (RL), is considered. A coupled finite element model leading to mechanical (displacements) and electrical (charges at electrodes) degrees of freedom is considered. An analysis of the effect of parametric uncertainties of the device on the electric output is performed. Piezoelectric and dielectric constants of the piezoelectric active layers and electric circuit equivalent inductance are considered as stochastic parameters. Mean and confidence intervals of the electric output are evaluated.

The reduced modulus in the analysis of sharp instrumented indentation tests

In the analysis of instrumented indentation data, it is common practice to incorporate the combined moduli of the indenter (E-i) and the specimen (E) in the so-called reduced modulus (E-r) to account for indenter deformation. Although indenter systems with rigid or elastic tips are considered as equivalent if E-r is the same, the validity of this practice has been questioned over the years. The present work uses systematic finite element simulations to examine the role of the elastic deformation of the indenter tip in instrumented indentation measurements and the validity of the concept of the reduced modulus in conical and pyramidal (Berkovich) indentations. It is found that the apical angle increases as a result of the indenter deformation, which influences in the analysis of the results. Based upon the inaccuracies introduced by the reduced modulus approximation in the analysis of the unloading segment of instrumented indentation applied load (P)-penetration depth (delta) curves, a detailed examination is then conducted on the role of indenter deformation upon the dimensionless functions describing the loading stages of such curves. Consequences of the present results in the extraction of the uniaxial stress-strain characteristics of the indented material through such dimensional analyses are finally illustrated. It is found that large overestimations in the assessment of the strain hardening behavior result by neglecting tip compliance. Guidelines are given in the paper to reduce such overestimations.

Analysis, manufacture and characterization of Ni/Cu functionally graded structures

In this work, an experimental and numerical analysis and characterization of functionally graded structures (FGSs) is developed. Nickel (Ni) and copper (Cu) materials are used as basic materials in the numerical modeling and experimental characterization. For modeling, a MATLAB finite element code is developed, which allows simulation of harmonic and modal analysis considering the graded finite element formulation. For experimental characterization, Ni-Cu FGSs are manufactured by using spark plasma sintering technique. Hardness and Young's modulus are found by using microindentation and ultrasonic measurements, respectively. The effective gradation of Ni/Cu FGS is addressed by means of optical microscopy, energy dispersive spectrometry, scanning electron microscopy and hardness testing. For the purpose of comparing modeling and experimental results, the hardness curve, along the gradation direction, is used for identifying the gradation profile; accordingly, the experimental hardness curve is used for approximating the Young's modulus variation and the graded finite element modeling is used for verification. For the first two resonance frequency values, a difference smaller than 1% between simulated and experimental results is obtained. (C) 2012 Elsevier Ltd. All rights reserved.

Assessment of nose protector for sport activities: finite element analysis

There has been a significant increase in the number of facial fractures stemming from sport activities in recent years, with the nasal bone one of the most affected structures. Researchers recommend the use of a nose protector, but there is no standardization regarding the material employed. Clinical experience has demonstrated that a combination of a flexible and rigid layer of ethylene vinyl acetate (EVA) offers both comfort and safety to practitioners of sports. The aim of the present study was the investigation into the stresses generated by the impact of a rigid body on the nasal bone on models with and without an EVA protector. For such, finite element analysis was employed. A craniofacial model was constructed from images obtained through computed tomography. The nose protector was modeled with two layers of EVA (1 mm of rigid EVA over 2 mm of flexible EVA), following the geometry of the soft tissue. Finite element analysis was performed using the LS Dyna program. The bone and rigid EVA were represented as elastic linear material, whereas the soft tissues and flexible EVA were represented as hyperelastic material. The impact from a rigid sphere on the frontal region of the face was simulated with a constant velocity of 20 m s-1 for 9.1 mu s. The model without the protector served as the control. The distribution of maximal stress of the facial bones was recorded. The maximal stress on the nasal bone surpassed the breaking limit of 0.130.34 MPa on the model without a protector, while remaining below this limit on the model with the protector. Thus, the nose protector made from both flexible and rigid EVA proved effective at protecting the nasal bones under high-impact conditions.

Limiting the influence of friction on the split Hopkinson pressure bar tests by using a ring specimen

The deformation of a ring under axial compression is analyzed in order to estimate a favorable ring specimen geometry capable of limiting the influence of friction on the stress-strain curve obtained from SHPB tests. The analysis shows that the use of a ring specimen with a large inner diameter and a small radial thickness offers some advantages comparing with the traditional disk sample. In particular, it can improve the reliability of the test results for ductile materials in the presence of friction. Based on the deformation analysis of a ductile ring under compression, a correction coefficient is proposed to relate the actual material stress strain curve with the reading from the SHPB. It is shown using finite element simulation that the proposed correction can be used for a wide range of conventional ductile materials. Experimental results with steel alloys indicate that the correction procedure is an effective technique for an accurate measurement of the dynamic material strength response. (C) 2012 Elsevier Ltd. All rights reserved.