Investigation of biogenic materials and ferroelectrets for energy harvesting on vibrating aircraft structures

In this publication the application of novel piezoelectric materials for energy harvesting on vibrating aircraft structures is investigated. These materials have significant advantages over conventional piezoelectric transducer materials like piezoceramics. In particular, biogenic materials in the form of wood-based materials and ferroelectrets in the form of irradiation cross-linked polypropylene are the subject of the investigation. The material characterization in terms of mechanical and electromechanical properties is shown for both material types. For the wood materials a compression test is used as the material has load-bearing properties. The ferroelectrets provide high compliances and are therefore investigated in a tensile test for material characterisation as well as in a four-point flexural test regarding its behaviour when glued to a dynamically bending surface. Additionally an FE-model of the material model for ferroelectrets is presented, which is validated by experimental results. An estimation of the power output is given for different concepts with both kinds of materials.


Introduction
Climate change undoubtedly represents one of the main challenges of human kind of the current and future centuries. The still growing world population with a simultaneous striving of wealth including mobility for more and more people requires new technological innovations and a fundamental change in infrastructure and human habits. In 2016 air traffic caused 2.8% of carbon dioxide emissions [1], a number that could rise in case of a very likely increase of flight passengers around the globe. More ecologically friendly air traffic will be a part of the solution. In every phase of 1 3 aircraft operation, lightweight structures are prone to vibrations, whose mechanical energy remains unused. Within the framework of the project EnerVib (Innovative Energy Generation Concepts from Vibrating Aircraft Structures) the research question is clarified how and to what extent the vibration energy can be converted into electrical energy effectively. The aim is to power microelectronics like sensors in difficultly accessible areas of aircraft. Thus cables with their associated weights can be reduced and eventually fuel can be saved. The investigation focuses in particular on biogenic materials and ferroelectrets, which are to be integrated into the aircraft as transducer materials for flat or freely oscillating Energy-Harvesters (EH). With a high ecoefficiency of wood materials and a high structural flexibility of ferroelectrets, these materials have significant advantages over conventional piezoelectric materials. As an example of an aircraft-component a wingbox (cf. Fig. 1) is chosen providing cyclic deformations in the form of vibrations for Energy Harvesting with both material kinds. In Holzmann et al. [2] it is shown that realistic strain magnitudes can be calculated to around 4.5 ⋅ 10 −5 m ⋅ m −1 at a frequency of around 1.5 Hz.
Due to its load-bearing properties piezoelectric wood materials can be integrated directly into the flux of forces, replacing material of the wing structure e.g. of ribs, stringers, struts or cross beams. To ensure stability of the mechanical structure, sandwich-arrangements of wood and metals are a possible solution. Ferroelectrets on the other hand are highly compliant materials that can be attached to large surfaces of the structure using local surface strains for energy conversion. Therefore ferroelectrets are passively applied to the structure in parallel to the force flux. Because of the difference in application, different kinds of experiments are conducted for the two material types in this work.

Piezoelectric constitutive equations
The material characterisation and numerical simulation in the following investigations of piezoelectric wood-based materials and ferroelectrets are based on physical fundamentals. The piezoelectric properties are coupled with the mechanical properties via the piezoelectric constitutive equations and with compliances s ij , piezoelectric charge constants ij and permittivities ij . S is the strain tensor, T the stress tensor, E the electric field and D the electric charge density. For sufficiently small surface elements, the electric field can be considered constant and the piezoelectric material can be approximated through a plate capacitor. For a constant electric field E = 0 Eq. 2 is simplified to By including with the charge q and the assumption that the specimen is regarded as a plate capacitor with the electrode and force application surfaces A E and A f the equation can be transformed to with the force f as axial load. For an alternating current as the output of a harmonically excited piezoelectric material with the frequency , the power output can be expressed as with the electrical voltage u el , optimal load resistance R opt = ⋅ C , the frequency and the capacitance of the piezoelectric transducer C. Using the piezoelectric constitutive equations with open electrodes ( D = 0 ) and by eliminating u el the power output subject to a strain excitation S can be expressed as with the Young's modulus Y, the surface Area A piezo of the transducer electrode and the thickness t piezo .

Wood materials
Due to its chemical and mechanical properties, wood is applied in the most diverse areas. As a natural and renewable raw material it is used for energetic purposes as well as a construction material. Because of its specific strength and rigidity the fibre composite material comes into effect in lightweight design for load-bearing structures. In addition, wood has good fatigue strength and a fail-safe character. Furthermore, wood has piezoelectric properties, which were first discovered in the 1940s and are still the subject of current research activities [4]. Figure 2 shows the schematic cell structure of cellulose. The piezoelectric effect is made possible by a crystalline structure, which consists of different cell layers. If shear stresses are induced into the structure, a displacement of charge centres occurs, which results in electric dipole moments and a polarization. Its magnitude depends on the applied load and on the parameters of the wood itself, which are mainly limited to the type of wood, humidity and temperature. The piezoelectric charge constant of Eq. 5 can be used as a characteristic value for the conversion of mechanical force into electric charge. For wood, the charge constant is in the order of 5% of that of quartz crystals [5]. In a direct comparison of the transducer characteristics of wood to piezo crystals and especially to modern piezoceramics (Table 4), piezoceramics exceed wood materials. In addition, the low piezoelectric charge constant of wood complicates energy harvesting, because the low currents must first cover the power loss of the energy processing before any electrical energy can be stored and used by consumers. But to compensate for this problem, the wood offers a significant advantage over conventional piezo materials. Due to its mechanical properties, large volumes of the transducer material can be integrated. The aim is to substitute the existing structures with piezoelectric load-bearing structures. This substitution can occur without significant weight increase, as shown in Fig. 3 for aluminium.
A beam cross-section made of aluminium, which is used as a reference is compared to three selected cross-sections with the same bending stiffness. In option one, a cross-section is shown that must have the same height as the reference due to the construction space, but is made wider. The crosssection of Wood Beam 1 has a mass increase of 20 percent compared to the aluminium cross-section. In the case where there are no limits on size, the aluminium cross-section can be substituted without an increase in mass by using Wood Beam 2. The third option with enormous lightweight potential is the composite beam. If there are construction space limitations, a cross-section can be designed from a variable composite of several materials. In this example, 25 percent of the width of the aluminium of the reference cross-section was substituted by wood. The increase in mass is only five percent. This offers a realistic potential for integrating the piezoelectric properties of the wood.
The production process of wood-based energy harvesters brings further advantages. While piezoceramics are timeconsuming and thus cost-intensive to produce, certain production steps, such as the polarization imprinting after the sintering process, are not required for wood-based materials. In addition to economic advantages, wood is a renewable and eco-friendly material. Its environmental friendliness results from its sustainable regeneration and good recyclability. At the Institute of Structural Mechanics and Lightweight Design (SLA) at RWTH Aachen University multifunctional structures and mechanisms are developed, which combine the advantages of wood-based materials and thus create a pioneering concept of energy harvesting.  Figure 4 illustrates two basic conceptual approaches. On the one hand, flat, integrative structures are developed that combine the strength, stiffness and piezoelectric properties of wood. The focus is on ribs, skin panels, stringer profiles, struts and cross beams. Hybrid materials have a laminate structure of fibre composite materials and the transducer layers of wood incorporated into the supporting structure, which are idealized as a disk. A further side effect is the increase in intrinsic material damping due to the possible conversion of deformation energy into thermal energy. In addition to multifunctional structures, the development of absorber mechanisms is a major motivation. These are frequency-tuned, cantilever-resonators which are applied to vibrating structures. Vibrations are induced into the absorber masses by base excitations, which lead to excessive amplitudes at the resonance frequencies. The spring elements, which form the connection of the absorber mass and the base point, are made of piezoelectric material. With these two concepts, a wide range of vibrations can be used to convert mechanical energy into electrical energy in aircraftstructures. The possibility of passive vibration suppression with the multifunctional structures and mechanisms creates a double-use. Moreover, the increase in structural durability and the reduction of acoustic vibrations have further synergistic effects that qualify the use of wood-based materials as a transducer material. As mentioned above, the tests conducted in this work focus the first of the two concepts.

Ferroelectrets
In addition modern ferroelectrets represented by irradiation cross-linked polypropylene (IXPP) are investigated. Due to high piezoelectric transducer coefficients and high structural flexibility, these electromechanical transducer materials combine the advantages of piezoceramics and electroactive polymers. Electrets provide a permanent electric field caused by a forced charge separation in the material itself. This charge separation can be induced by a so-called corona charging, leading to an electrical breakdown in the dielectric. Through metallic electrodes vaporized onto the material, compensation charges appear when connected to an electrical circuit. Through a cyclic deformation of the material, the electric field changes periodically and hence the compensation charges can be used as an alternating electric current. The electret-material has gained piezoelectric properties and is therefore called piezo-of ferroelectret. In Fig. 5 two types of ferroelectrets are shown, that are investigated in the Research group System Reliability, Adaptive Structures, and Machine Acoustics (SAM).
One type is irradiated cross-linked polypropylene (IXPP) as shown on the left. It provides an increased temperature stability of the piezoelectric transducer coefficients compared to regular polypropylene ferroelectrets achieved by a radiation process of the foam material. A hot-pressing process reduces the thickness and increases the compliance of the material in thickness-direction at the same time. After the corona charging process electric surface-charges emerge on the voids of the foam, that cause a quasi-permanent electric field [8]. IXPP is mostly used as a 33 -transducer. However, the material also works in the 31-mode. IXPP usually provides 33 -constants in the range from 300 up to 700 pC/N being in the same order of magnitude as piezoceramics. The charge constants perpendicular to the polarization axis are about two orders of magnitude lower [9,10]. Another type is parallel tunnel fluorethylenepropylene (FEP), having a macroscopic tunnel-structure, that is formed by a fusion of two waved layers of FEP. The structure is shown on the right-hand-side of Fig. 5. The tunnels work similar to the air-voids of the foams as macroscopic cavities separating the charges after corona charging. When the material is exposed to an external load perpendicular to the tunnels length axis, the tunnels deform and an electric potential can be measured between the electrodes. Parallel-tunnel FEP is designed as a d 31 -transducer despite having high 33 -constants of about 1300 pC/N as well [9]. Ferroelectrets are a subject of research for Energy Harvesting only since 2012 and provide a huge potential for the use in lightweightstructures. Analytical models for the electrical power-output of ferroelectret-EH have been developed by Pondrom for harvesters working in both the 33 and 31 -mode [11]. These models describe systems with one degree of freedom including a seismic mass. The systems are optimized for a certain frequency. Regarding ferroelectrets the project addresses the question of research, if and to what extend distributed EH based on ferroelectrets can be used in a wide frequency range. This work focuses on the ferroelectret IXPP, due to the fact that it can be produced in larger sheets and that it is more robust than parallel-tunnel FEP.

Investigation of piezoelectric wood-based materials
There are various investigations on the piezoelectric behaviour of wood-based materials. These include the most diverse types of wood, which have, however, been insufficiently investigated with regard to major influential parameters such as humidity and temperature. The results show significant deviations and also contradictions. Therefore, only few conclusions can be derived to select a suitable wood type as transducer material. Specific values characterizing the performance must therefore be determined. In the course of this problematic, an extensive material characterization is carried out to determine the electro-mechanical properties of different biogenic transducer materials. To measure the piezoelectric effect in wood-based materials, a test setup is designed and verified. As a fundamental for the characterization of the electro-mechanical properties of different biogenic transducer materials, the piezoelectric effect is first metrologically registered on a wooden cuboid and the results are compared with the literature references. Thus, the general test setup including the measurement chain are verified. Figure 6 shows the general test setup with specimen, insulation, testing machine and charging cable. Currently there are no normative specifications defined. Hence the wooden cuboid is selected as a simple specimen geometry which can be regarded as a plate capacitor. First, the specimen is milled out of a solid wood block perpendicular to the LR plane at a grain angle of = 45 • (see Fig. 2). With axial load in the X or Y direction, the greatest possible shear is thus obtained and the piezoelectric effect is at its maximum. A conductive bronze coating is applied to the two surfaces in the LR plane, forming the electrodes. According to theory, they are polarized and thus referred to as polarization surfaces in the following. Electric plugs are attached to this coating with silver conductive adhesive, to which the highly insulated charge cable is connected. Investigations of various insulation materials have shown that plastics such as the highly insulating Teflon are particularly unsuitable for insulating samples and specimen holders, since static charging of these materials is accompanied by changes in the electrical field, so that the piezoelectric effect can no longer be differentiated from these disturbances. Glass ceramic plates (Macor) do not show this effect and are therefore chosen for the insulation between specimen and the specimen holder of the hydraulic testing machine type SCHENK POZ 160. This machine applies a cyclical pressure load to the specimen.  [6,7], c schematic structure in thickness direction with charge distribution

Material characterization
To be able to measure the charges applied to the electrodes, a charge amplifier type KISTLER 5015 is connected to the cables.
To determine material parameters compression tests were performed on wood samples with a grain angle of 45 • to the direction of loading. From these tests, a Young's modulus of 1334 MPa and a yield strength of 12.2 MPa were determined as mechanical parameters. The stress-strain graphs are shown in Fig. 7

Comparison of theoretical polarization surface and reference surface
For this comparison, a cuboid of European beech with a wood moisture content of 8% with four identical electrode surfaces is investigated. Two electrodes were applied to the polarization surfaces and two to reference surfaces perpendicular to them. The two samples were subjected to a swelling pressure load of the form with f 0 = 1 kN, = 2 Hz and f = 0.25 kN. The test setup with the specimen is shown in Fig. 8. Figure 9 shows the graphs of charge and force for the polarization surface A Pol and the reference surface A Ref . Both surfaces have a sinusoidal charge curve that follows the force curve. There are significant differences between the maximum amplitudes of the charge lines. As described theoretically, the polarizations occur at the expected surfaces. At the same time, there are also charge changes at the reference surfaces that can be measured.

Comparison of 0 • and 45 • grain to load direction
For this comparison, two cubes made of European beech with a wood moisture content of 8% and edge lengths of 30mm are investigated. The grain directions were differentiated between = 0 • and = 45 • . The cuboids were set under  the same pressure load as the cube in Sect. 2.2. Under this load, no piezoelectric effect is expected to be registered in the 0 • due to the absence of shear stress, which is a premise for the polarization of the electrodes. On the other hand, a sinusoidal load curve is expected for the 45 • -specimen with a sinusoidal force curve. Figure 10 shows the curves of charge and force for the 0 • and 45 • samples. As expected, the 0 • specimen does not feature a piezoelectric effect. The 45 • -specimen shows a sinusoidal charge curve that follows the force curve. According to Eq. 5 the piezoelectric charge constant is calculated to be 0.044 pC/N and is thus in the range from 0.037 to 2.3 pC/N given in the literature [5,12]. Compared to typical piezoelectric materials (Table 4) the piezoelectric properties of wood are significantly lower. Despite this, the following example demonstrates the potential for using this material as an energy harvester.Integrating wood with a mass of 1 kg and a thickness of 5mm, under ideal load conditions with the assumptions of a frequency of 1.5 Hz at a strain of 4.5 ⋅ 10 −5 m ⋅ m −1 an electrical energy of 2.37 nW can be achieved according to Eq. 7. Through applied lightweight design, this quantity can be expanded without significant weight increase and thus wood gains more importance as an energy harvester.

Investigation of ferroelectrets
The application and experimental investigation of ferroelectrets on real aircraft-structures is costly and can be problematic considering safety standards. Therefore the focus is the experimental investigation of simplified structures with EH and especially the simulation of these. Due to the fact that aircraft-structures have curved surfaces, analytical models as the ones by Pondrom [11] are neglected for the simulation of the system. Instead a numerical modelling approach is pursued. For the latter one, first a material model of the ferroelectret IXPP is defined and in the next step the system investigated in this work is outlined. Moreover the experimental setup is described and the validation of the model with experimental data is presented.

IXPP material model
For the simulation of the material a linear model is used, because small local strains of aircraft base structures are assumed to affect the material. This assumption seems reasonable, as polymer foams show a linear elastic behaviour for small deflections [13]. In Fig. 11 a unit cell of the material is shown with a local coordinate system. The coordinates 1 and 2 represent in-plane directions, the number 3 the direction of polarisation perpendicular to the 1-and 2-direction. Based on the literature [14] and the manufacturing process it is assumed, that the parameters in 1-and 2-direction are in the same order of magnitude so they are assumed identical. The parameters in 3-direction have different values. Hence, a linear transverse isotropic material model is chosen in this work. The parameters are included in Eq. 1 by means of compliances s ij , piezoelectric charge constants ij and permittivities ij . Additionally to the massdensity , the equations are imported into a commercial numerical simulation software. Viscoelastic and non-linear effects are neglected in this work.

Material characterization
To determine all relevant material properties from Sect. 3.1, measurements on several IXPP-samples are used aside manufacturer data. With a capacitance measurement the static capacitance of the samples is determined. For samples with a size of 60 by 60 mm 2 an average capacitance of 302.4 pF results with a standard deviation of 27.69%. Since a planar attachment of the ferroelectrets without a seismic mass to a vibrating structure primarily leads to a deformation of the ferroelectret in the directions 1 and 2, the derivation of the Young's modulus in 3-direction is not described at this point. Typical values are 0.54 [14] to 1.1 MPa [10]. By a tensile test the Young's modulus in plane direction is measured, as shown in Fig. 12. For this purpose four samples without electrode are tested with a given strain amplitude.
Using the slope of the graphs, the Young's modulus can be calculated to 133.3 N∕m 2 in average as shown in Table 2. Furthermore a tensile test of an electroded specimen is conducted. After prestressing the specimen a force amplitude of f = 5 N and a frequency of = 2 Hz are applied. For the charge measurements, a laboratory charge amplifier was used that ensures the absence of any load-resistance by very  high input impedance. According to the charge amplifier principle, the generated charge is defined by the negative feedback capacitor (This also applies to the charge measurements for the investigated wood materials). The corresponding graph is shown in Fig. 13. The average charge amplitude is determined as 372.3 pC. Using Eq. 5 the piezoelectric charge constant 31 is further calculated as 0.252 pC/N. As mentioned in Sect. 1.3, the investigated IXPP-ferroelectret is designed for the use in its 33-mode.

Four-point flexural test
In the case of a planar attachment of the ferroelectrets onto aircraft-components, surface strains of the vibrating structure are transferred directly to the ferroelectret-EH. Due to the much lower stiffness and density of the polymer compared to materials in aircraft like Aluminium or CFRP, the mechanical stiffness-and mass-properties of the aircraftstructure stay approximately the same. In a four-point flexural test the behaviour of IXPP-ferroelectrets glued to Aluminium-plates is investigated like outlined in Fig. 14.
The force f acts as a symmetric excitation on both ends of the plate. The aluminium plate is pressed symmetrically at two rotatable bearings with its upper side.
By design a constant momentum is caused in the sample holder by the application of the force f between the bearings and therefore a constant curvature and surface strain is achieved. When the plate and herewith the IXPP-sample is exposed to a load, the resulting electric voltage at the electrodes is measured.
The dimensions of the sample holder are adapted to the tested sample size with an edge length of 60 mm. The constants for the experimental setup are shown in Table 3. For a preliminary design the maximum force is calculated with a bending line. The parameters in Fig. 14     The constants can be derived from the boundary conditions w(b) = w(c) = 0 as and By means of the bending line, the maximum deflection of the sample holder exposed to the force f can be estimated. Taking a yield strength of 100 MPa the maximum force not leading to a plastic deformation of the aluminium structure of 1 mm thickness is calculated as 244.5 N. The corresponding maximum deflection and surface strain of the plate are −3.522 mm and 0.141% respectively. The emerging surface strain through a force of 1 N is calculated to 5.767 ⋅ 10 −4 %.

Finite element model
With the system shown in Fig. 14, a finite element (FE) model is developed based on the geometry data in Table 3 discretization the aluminium plate and the IXPP-ferroelectret by means of hexahedric elements. The equation of motion for the mechanical system can be expressed in its most general form as Reliability (LBF) in Darmstadt. In MATLAB a model reduction of the piezo-mechanical structure and a conversion to a state-space representation is conducted as described in [15, p. 77].

Experimental validation
For the validation of the created model, an experimental setup is used as shown in Fig. 16 representing the system in Fig. 14. Aluminium-rolls are both used for the rotatable bearings and for the contact point for the induction of the shaker force. At the connection point between shaker and base plate, both force and acceleration can be measured.
In a first test, the electrical response of the electret in terms of voltage with open electrodes is measured when applying a static load to the setup shown in Fig. 16. For different forces the generated voltage is shown in Fig. 17. The simulated behaviour fits the experimental results in the investigated load range.
For the system identification, the transfer function of acceleration over excitation force at the connection point is measured. Additionally the transfer function of the electrode voltage over excitation force is identified. Both graphs are shown in Fig. 18 with their corresponding simulated behaviour.
A dominant resonance can be observed in both the mechanical and electrical transfer function, representing the first bending mode of the plate. A model fit is achieved by choosing a mass of 0.75 kg for the base plate of the setup and a structural damping of 3% . The damping matrix is implemented by a Rayleigh damping model. A good agreement of the measured and simulated curves can be achieved however in the higher frequency range the measured mechanical-electrical transfer function shows a lower amplitude than the simulated one. An interesting effect is furthermore perceived by the emerging peaks between 200 and 150 Hz only in the measured mechanical-electrical transfer function. These can be explained by local asymmetric vibration modes of the plate, that cannot be observed in the mechanical transfer function because of the centred position of the sensor. Furthermore the aluminiumplates cannot be placed in the setup perfectly symmetrical like in the simulation. In the frequency range up to 200 Hz the simulations based on the developed material model show a good agreement with experimental results. Therefore numerical methods will be further used in future works to determine the electrical power-output when applying ferroelectrets to more complex models of aircraft-components. For the case of the assumed strain amplitude in the winbox (Fig. 1) a power output of around 3.6 ⋅ 10 −13 W is calculated using Eq. 7, which is very low since the material is used in its 31-mode. However by using a simple transmission as shown in Fig. 19 to convert the in-plane strain of the host structure into a compression of IXPP-sheets in 33-mode the power output is increased tremendously. This is firstly due to the strain increase in the material and secondly due to the much higher 33 -coefficients. The transmission of in-plane strain S 1 to the strain in the ferroelectret S 3 is calculated by Using a realistic transmission ratio of 1/20 (i. e. an angle of 2.9 • ) with a projected mechanism length of 0.1 m a strain S 3 of 0.018 mm −1 results. Mechanical losses are neglected due to the high compliance of IXPP compared to the metallic structures. With a stack of 15 IXPP-sheets with a respective surface area of 1 cm 2 the power output is increased to 4.2 μW using Eq. 7. Although IXPP-ferroelectrets have much lower 31 -constants compared to typical piezoelectric materials (Table 4), they can effectively be used for energy harvesting in bending host structures.

Conclusion
Alternative piezoelectric materials for Energy Harvesting on vibrating aircraft structures were presented in the form of biogenic wood materials with load-bearing properties and a high eco-efficiency and ferroelectrets with remarkable structural flexibility and high piezoelectric transducer coefficients. For both materials the piezoelectric effect was registered with a respective experimental setup. For wood materials the assumptions from the literature were confirmed with a determined charge constant of = 0.044 pC/N. The low piezoelectric charge constant of wood complicates energy harvesting. The lightweight potential of wood, which enables the substitution of existing loadbearing structures, can compensate for this disadvantage. The test setup and the measurement chain were verified for wood materials with regard to environmental disturbances. An initial specimen geometry was defined for the material characterization. Compression tests were performed on wood samples with a grain angle of 45 • to the direction of loading. Young's Modulus and the yield strength were determined. Investigations on different wood species under variation of the parameters will follow in future works. From this, specific values will be derived to describe the dynamic electro-mechanical performance of the transducer materials. Approaches for normative test methods will be developed to standardize the procedure. After selecting a potential material and the optimal geometry, structural integration will follow. For ferroelectrets it was shown, that the static and dynamic modelling for small structural strains of host structures is possible with linear material models, that are also used for piezoceramics. The material parameters for the model were derived from manufacturer data, capacitance measurements and by means of a tensile test. Furthermore a preliminary design of a dynamic four-point flexural test for the investigation of IXPP-ferroelectrets glued on aluminium plates was used to create an electromechanical numerical model. Finally the model was validated with measurements from an experimental test setup regarding its mechanical and electromechanical transfer functions. In the frequency range up to 200 Hz the simulations based on the developed material model showed a good agreement with experimental results. Despite not being designed for the use in 31-mode the investigated IXPP-ferroelectret showed a charge constant of 31 = 0.252 pC/N. This work represents the first step for the dynamic modelling of ferroelectrets attached to surfaces, without the aim of maximizing the electrical power output. In future works, optimizations regarding the ferroelectret material and the geometry of the EH will be in focus of research, as well as their experimental validation. A concept of a lightweight energy harvester design was introduced to compensate for the disadvantage of low d 31 -coefficients, by converting the in-plane strain energy to a perpendicular compression of a ferroelectret stack. First power estimations of a specific mechanism nad a stack with a surface area of 1 cm 2 showed a power output of 4.2 μW . Furthermore the modelling approach will also be used for wood materials. For both material kinds the electrical power-output subject to more realistic load conditions in aircraft structures will be calculated and optimized using more complex numerical models. The simulations will further be validated on simplified demonstrator components.