<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article">
 <front>
  <journal-meta>
   <journal-id journal-id-type="publisher-id">
    eng
   </journal-id>
   <journal-title-group>
    <journal-title>
     Engineering
    </journal-title>
   </journal-title-group>
   <issn pub-type="epub">
    1947-3931
   </issn>
   <issn publication-format="print">
    1947-394X
   </issn>
   <publisher>
    <publisher-name>
     Scientific Research Publishing
    </publisher-name>
   </publisher>
  </journal-meta>
  <article-meta>
   <article-id pub-id-type="doi">
    10.4236/eng.2025.177023
   </article-id>
   <article-id pub-id-type="publisher-id">
    eng-144473
   </article-id>
   <article-categories>
    <subj-group subj-group-type="heading">
     <subject>
      Articles
     </subject>
    </subj-group>
    <subj-group subj-group-type="Discipline-v2">
     <subject>
      Engineering
     </subject>
    </subj-group>
   </article-categories>
   <title-group>
    Piezoelectric Energy Harvesting in Quadcopter UAVs Based Modeling
   </title-group>
   <contrib-group>
    <contrib contrib-type="author" xlink:type="simple">
     <name name-style="western">
      <surname>
       Hajir
      </surname>
      <given-names>
       Saeed
      </given-names>
     </name> 
     <xref ref-type="aff" rid="aff1"> 
      <sup>1</sup>
     </xref>
    </contrib>
    <contrib contrib-type="author" xlink:type="simple">
     <name name-style="western">
      <surname>
       Fath Elrahman Ismael
      </surname>
      <given-names>
       Ahmed
      </given-names>
     </name> 
     <xref ref-type="aff" rid="aff1"> 
      <sup>1</sup>
     </xref>
    </contrib>
    <contrib contrib-type="author" xlink:type="simple">
     <name name-style="western">
      <surname>
       Gihad Abdelaziz Abdelghani
      </surname>
      <given-names>
       Ibrahim
      </given-names>
     </name> 
     <xref ref-type="aff" rid="aff1"> 
      <sup>1</sup>
     </xref>
    </contrib>
    <contrib contrib-type="author" xlink:type="simple">
     <name name-style="western">
      <surname>
       Khalil Babeker
      </surname>
      <given-names>
       Ahmed
      </given-names>
     </name> 
     <xref ref-type="aff" rid="aff2"> 
      <sup>2</sup>
     </xref>
    </contrib>
   </contrib-group> 
   <aff id="aff1">
    <addr-line>
     aCollege of Engineering, Sudan University of Science and Technology, Khartoum, Sudan
    </addr-line> 
   </aff> 
   <aff id="aff2">
    <addr-line>
     aDepartment of Electrical Engineering, Alimam Alhadi College, Khartoum, Sudan
    </addr-line> 
   </aff> 
   <pub-date pub-type="epub">
    <day>
     23
    </day> 
    <month>
     07
    </month>
    <year>
     2025
    </year>
   </pub-date> 
   <volume>
    17
   </volume> 
   <issue>
    07
   </issue>
   <fpage>
    388
   </fpage>
   <lpage>
    401
   </lpage>
   <history>
    <date date-type="received">
     <day>
      15,
     </day>
     <month>
      June
     </month>
     <year>
      2025
     </year>
    </date>
    <date date-type="published">
     <day>
      27,
     </day>
     <month>
      June
     </month>
     <year>
      2025
     </year> 
    </date> 
    <date date-type="accepted">
     <day>
      27,
     </day>
     <month>
      July
     </month>
     <year>
      2025
     </year> 
    </date>
   </history>
   <permissions>
    <copyright-statement>
     © Copyright 2014 by authors and Scientific Research Publishing Inc. 
    </copyright-statement>
    <copyright-year>
     2014
    </copyright-year>
    <license>
     <license-p>
      This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
     </license-p>
    </license>
   </permissions>
   <abstract>
    This paper presents a detailed investigation into piezoelectric energy harvesting for Unmanned Aerial Vehicles (UAVs) through advanced computational modeling and simulation. The study focuses on optimizing the integration of piezoelectric materials (PZT-5H) into UAV wings to convert mechanical vibrations—induced by aerodynamic forces—into usable electrical energy. A high-fidelity 3D model of a DJI Mini 3 Pro quadcopter was developed, and Finite Element Analysis (FEA) combined with Computational Fluid Dynamics (CFD) was employed to assess energy harvesting efficiency under varying flight conditions. Key findings indicate that optimal placement of piezoelectric patches near high-stress zones yields voltage outputs ranging from 1.54 V to 
    <b> 24.3</b>V with peak performance near resonant frequencies. The study also explores the impact of structural modifications, such as adding a steel fin, to enhance vibration transmission.
   </abstract>
   <kwd-group> 
    <kwd>
     Piezoelectric
    </kwd> 
    <kwd>
      Energy Harvesting
    </kwd> 
    <kwd>
      Simulation
    </kwd> 
    <kwd>
      UAVs
    </kwd>
   </kwd-group>
  </article-meta>
 </front>
 <body>
  <sec id="s1">
   <title>1. Introduction</title>
   <p>Unmanned Aerial Vehicles (UAVs) play a crucial role across multiple sectors, particularly in tasks like surveillance and reconnaissance <xref ref-type="bibr" rid="scirp.144473-1">
     [1]
    </xref>. As the demand for longer flight times and increased autonomy continues to grow, one of the primary challenges is ensuring a reliable and sufficient power supply for the onboard systems. Conventional battery power sources often present limitations regarding weight, capacity, and operational duration. Considering these challenges, there is a significant push towards exploring alternative energy harvesting methods that are both sustainable and efficient, aiming to either complement or replace traditional battery solutions <xref ref-type="bibr" rid="scirp.144473-2">
     [2]
    </xref>.</p>
   <p>Among the various techniques for energy harvesting, piezoelectric energy harvesting is particularly noteworthy due to its capability to transform mechanical vibrations into electrical energy. During flight, the UAV’s wings constantly experience dynamic motion, which generates vibrations that can be effectively harnessed using piezoelectric materials. The vibrations that arise from aerodynamic forces are an untapped resource of energy, presenting an opportunity to power onboard systems, decrease dependency on conventional batteries, and improve the overall operational efficiency of UAVs <xref ref-type="bibr" rid="scirp.144473-3">
     [3]
    </xref>-<xref ref-type="bibr" rid="scirp.144473-5">
     [5]
    </xref>.</p>
   <p>This study delves into the mechanics of UAV wing vibrations and investigates the potential of employing piezoelectric materials for energy harvesting. It focuses on assessing the practicality and efficiency of integrating piezoelectric energy harvesting systems into UAV designs by evaluating the various contributions across the wing structure. The study encompasses a theoretical examination of energy conversion via piezoelectric mechanisms, identifies the key factors that influence energy generation, and seeks optimizations for harvesting devices to adapt them effectively for UAV applications.</p>
   <p>The forthcoming sections of the study will elucidate the fundamental principles behind piezoelectric energy harvesting, analyze the dynamic characteristics of UAV wings, and assess the prospects for energy generation. This comprehensive approach aims to provide insights into how piezoelectric harvesters can be integrated into sustainable power systems for UAVs. The findings of the research are intended to contribute significantly toward the advancement of more efficient and sustainable UAV technologies.</p>
   <p>In this study, a functional drone-based pumping prototype was developed Quadcopter (DJI Mini 3 Pro) for the purpose of conducting an investigation on it at a later stage.</p>
  </sec><sec id="s2">
   <title>2. Literature Review</title>
   <p>A collection of studies centered on the application of piezoelectric materials for energy harvesting from the vibrations of unmanned aerial vehicle (UAV) wings. One key research piece by Zhang et al. delves into the feasibility of using these materials to convert oscillations during flight into usable energy. The study not only evaluates the dynamic response of UAV wings to aerodynamic forces but also addresses the challenges of optimizing the placement of piezoelectric elements and selecting appropriate materials to maximize energy generation efficiency <xref ref-type="bibr" rid="scirp.144473-6">
     [6]
    </xref>.</p>
   <p>
    <xref ref-type="bibr" rid="scirp.144473-"></xref>Furthering the investigation, Marini et al. conducted a feasibility study that involved analytical modeling and experimental work to explore various configurations of piezoelectric materials integrated into UAV wings. Their findings indicated that substantial amounts of energy could be harvested under specific flight conditions, suggesting that piezoelectric systems could serve as valuable auxiliary power sources for UAVs. This underscores the potential of piezoelectric materials in enhancing UAV performance through energy self-sufficiency <xref ref-type="bibr" rid="scirp.144473-7">
     [7]
    </xref>.</p>
   <p>Building on these individual efforts, a study by Li et al. introduced a hybrid energy harvesting system that combines both piezoelectric and electromagnetic transducers. The research highlights the advantages of integrating these two technologies, showing that the hybrid system yields more reliable and stable energy outputs compared to systems that rely only on piezoelectric harvesting. This innovation points to the possibility of developing more efficient and scalable energy solutions for UAV applications <xref ref-type="bibr" rid="scirp.144473-8">
     [8]
    </xref>.</p>
   <p>Yang et al. further contributed to the literature by focusing on the design and optimization of piezoelectric energy harvesters specifically for UAV use. Through finite element analysis (FEA) and experimental techniques, they investigated the optimal positioning and dimensions of piezoelectric patches on UAV wings. Their research emphasizes the importance of understanding aerodynamic forces and vibrational characteristics, suggesting that manipulation of operational frequencies can significantly enhance the effectiveness of energy harvesting systems in UAVs <xref ref-type="bibr" rid="scirp.144473-9">
     [9]
    </xref>.</p>
   <p>Overall, these studies demonstrate a growing interest in the application of piezoelectric materials for energy harvesting in UAVs. The combined efforts underscore the technical challenges involved in optimizing material placement and design while revealing the potential improvements in energy efficiency that could be achieved. The evolution of hybrid systems also illustrates a promising path toward more robust and sustainable energy solutions for unmanned aerial vehicles in the future.</p>
  </sec><sec id="s3">
   <title>3. Mathematical Framework</title>
   <p>Modeling the piezoelectric effects in response to mechanical stress requires solving coupled partial differential equations (PDEs) that describe both the mechanical and electrical behavior of the material. These equations govern the interaction between mechanical deformation and electrical response, which is the core mechanism of piezoelectric energy conversion <xref ref-type="bibr" rid="scirp.144473-10">
     [10]
    </xref>-<xref ref-type="bibr" rid="scirp.144473-13">
     [13]
    </xref>.</p>
   <p>Let’s define the mechanical stress-strain relationship</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         α 
       </mi> 
       <mrow> 
        <mi>
          i 
        </mi> 
        <mi>
          j 
        </mi> 
       </mrow> 
      </msub> 
      <mo>
        = 
      </mo> 
      <msub> 
       <mi>
         S 
       </mi> 
       <mrow> 
        <mi>
          i 
        </mi> 
        <mi>
          j 
        </mi> 
        <mi>
          k 
        </mi> 
        <mi>
          l 
        </mi> 
       </mrow> 
      </msub> 
      <msub> 
       <mi>
         ε 
       </mi> 
       <mrow> 
        <mi>
          k 
        </mi> 
        <mi>
          l 
        </mi> 
       </mrow> 
      </msub> 
      <mo>
        − 
      </mo> 
      <msub> 
       <mi>
         e 
       </mi> 
       <mrow> 
        <mi>
          i 
        </mi> 
        <mi>
          j 
        </mi> 
        <mi>
          k 
        </mi> 
       </mrow> 
      </msub> 
      <msub> 
       <mi>
         ξ 
       </mi> 
       <mi>
         k 
       </mi> 
      </msub> 
     </mrow> 
    </math> (1)</p>
   <p>where 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         α 
       </mi> 
       <mrow> 
        <mi>
          i 
        </mi> 
        <mi>
          j 
        </mi> 
       </mrow> 
      </msub> 
     </mrow> 
    </math> is stress tensor (describes deformation), 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         S 
       </mi> 
       <mrow> 
        <mi>
          i 
        </mi> 
        <mi>
          j 
        </mi> 
        <mi>
          s 
        </mi> 
       </mrow> 
      </msub> 
     </mrow> 
    </math> is the stiffness tensor which describe the mechanical properties of the material, 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         e 
       </mi> 
       <mrow> 
        <mi>
          i 
        </mi> 
        <mi>
          j 
        </mi> 
        <mi>
          k 
        </mi> 
       </mrow> 
      </msub> 
     </mrow> 
    </math> is piezoelectric stress coefficient coupling between mechanical stress and electric field and 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         ξ 
       </mi> 
       <mi>
         k 
       </mi> 
      </msub> 
     </mrow> 
    </math> is electrical field vector. 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mo> 
      </mo> 
      <msub> 
       <mi>
         ε 
       </mi> 
       <mrow> 
        <mi>
          k 
        </mi> 
        <mi>
          l 
        </mi> 
       </mrow> 
      </msub> 
     </mrow> 
    </math> is the strain tensor and defined by</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         ε 
       </mi> 
       <mrow> 
        <mi>
          k 
        </mi> 
        <mi>
          l 
        </mi> 
       </mrow> 
      </msub> 
      <mo>
        = 
      </mo> 
      <mfrac> 
       <mn>
         1 
       </mn> 
       <mn>
         2 
       </mn> 
      </mfrac> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mrow> 
        <mfrac> 
         <mrow> 
          <mo>
            ∂ 
          </mo> 
          <msub> 
           <mi>
             u 
           </mi> 
           <mi>
             k 
           </mi> 
          </msub> 
         </mrow> 
         <mrow> 
          <mo>
            ∂ 
          </mo> 
          <msub> 
           <mi>
             x 
           </mi> 
           <mi>
             l 
           </mi> 
          </msub> 
         </mrow> 
        </mfrac> 
        <mo>
          + 
        </mo> 
        <mfrac> 
         <mrow> 
          <mo>
            ∂ 
          </mo> 
          <msub> 
           <mi>
             u 
           </mi> 
           <mi>
             l 
           </mi> 
          </msub> 
         </mrow> 
         <mrow> 
          <mo>
            ∂ 
          </mo> 
          <msub> 
           <mi>
             x 
           </mi> 
           <mi>
             k 
           </mi> 
          </msub> 
         </mrow> 
        </mfrac> 
       </mrow> 
       <mo>
         ) 
       </mo> 
      </mrow> 
     </mrow> 
    </math>(2)</p>
   <p>By using Newton’s second law, the motion of a small volume of material is given by:</p>
   <p>
    <math display="inline" xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        ρ 
      </mi> 
      <mfrac> 
       <mrow> 
        <msup> 
         <mo>
           ∂ 
         </mo> 
         <mn>
           2 
         </mn> 
        </msup> 
        <msub> 
         <mi>
           u 
         </mi> 
         <mi>
           i 
         </mi> 
        </msub> 
       </mrow> 
       <mrow> 
        <mo>
          ∂ 
        </mo> 
        <msup> 
         <mi>
           t 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
      </mfrac> 
      <mo>
        = 
      </mo> 
      <mfrac> 
       <mrow> 
        <mo>
          ∂ 
        </mo> 
        <msub> 
         <mi>
           α 
         </mi> 
         <mrow> 
          <mi>
            i 
          </mi> 
          <mi>
            j 
          </mi> 
         </mrow> 
        </msub> 
       </mrow> 
       <mrow> 
        <mo>
          ∂ 
        </mo> 
        <msub> 
         <mi>
           x 
         </mi> 
         <mi>
           j 
         </mi> 
        </msub> 
       </mrow> 
      </mfrac> 
     </mrow> 
    </math>(3)</p>
   <p>where 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       ρ 
     </mi> 
    </math> is the material density, 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         u 
       </mi> 
       <mi>
         i 
       </mi> 
      </msub> 
      <mo> 
      </mo> 
     </mrow> 
    </math>is the displacement vector (describes material motion in the i-direction). The component 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mfrac> 
       <mrow> 
        <mo>
          ∂ 
        </mo> 
        <msub> 
         <mi>
           α 
         </mi> 
         <mrow> 
          <mi>
            i 
          </mi> 
          <mi>
            j 
          </mi> 
         </mrow> 
        </msub> 
       </mrow> 
       <mrow> 
        <mo>
          ∂ 
        </mo> 
        <msub> 
         <mi>
           x 
         </mi> 
         <mi>
           j 
         </mi> 
        </msub> 
       </mrow> 
      </mfrac> 
     </mrow> 
    </math> describes the divergence of the stress tensor (describes internal forces). This equation expresses how mechanical waves propagate while being influenced by the electric field.</p>
   <p>The electric displacement field 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         D 
       </mi> 
       <mi>
         i 
       </mi> 
      </msub> 
     </mrow> 
    </math> in a piezoelectric material is given by:</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         D 
       </mi> 
       <mi>
         i 
       </mi> 
      </msub> 
      <mo>
        = 
      </mo> 
      <msub> 
       <mi>
         e 
       </mi> 
       <mrow> 
        <mi>
          i 
        </mi> 
        <mi>
          j 
        </mi> 
        <mi>
          k 
        </mi> 
       </mrow> 
      </msub> 
      <msub> 
       <mi>
         ε 
       </mi> 
       <mrow> 
        <mi>
          j 
        </mi> 
        <mi>
          k 
        </mi> 
       </mrow> 
      </msub> 
      <mo>
        + 
      </mo> 
      <msub> 
       <mi>
         ε 
       </mi> 
       <mrow> 
        <mi>
          i 
        </mi> 
        <mi>
          j 
        </mi> 
       </mrow> 
      </msub> 
      <msub> 
       <mi>
         ξ 
       </mi> 
       <mi>
         j 
       </mi> 
      </msub> 
     </mrow> 
    </math> (4)</p>
   <p>In the analysis of piezoelectric actuators and sensors, electromechanical impedance is a critical quantity that characterizes how the system responds to external mechanical or electrical excitations <xref ref-type="bibr" rid="scirp.144473-14">
     [14]
    </xref>.</p>
   <p>By solving the previous system equations under harmonic excitation, to obtain how force responds to voltage, we obtain the impedance 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        Z 
      </mi> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mi>
         ω 
       </mi> 
       <mo>
         ) 
       </mo> 
      </mrow> 
     </mrow> 
    </math>, which is a complex function of angular frequency 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       ω 
     </mi> 
    </math>, relates the voltage 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        V 
      </mi> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mi>
         ω 
       </mi> 
       <mo>
         ) 
       </mo> 
      </mrow> 
     </mrow> 
    </math> applied to the force 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        F 
      </mi> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mi>
         ω 
       </mi> 
       <mo>
         ) 
       </mo> 
      </mrow> 
     </mrow> 
    </math> generated by the device:</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        Z 
      </mi> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mi>
         ω 
       </mi> 
       <mo>
         ) 
       </mo> 
      </mrow> 
      <mo>
        = 
      </mo> 
      <mfrac> 
       <mrow> 
        <mi>
          V 
        </mi> 
        <mrow> 
         <mo>
           ( 
         </mo> 
         <mi>
           ω 
         </mi> 
         <mo>
           ) 
         </mo> 
        </mrow> 
       </mrow> 
       <mrow> 
        <mi>
          F 
        </mi> 
        <mrow> 
         <mo>
           ( 
         </mo> 
         <mi>
           ω 
         </mi> 
         <mo>
           ) 
         </mo> 
        </mrow> 
       </mrow> 
      </mfrac> 
      <mo>
        = 
      </mo> 
      <mfrac> 
       <mn>
         1 
       </mn> 
       <mrow> 
        <msub> 
         <mi>
           k 
         </mi> 
         <mn>
           1 
         </mn> 
        </msub> 
        <mrow> 
         <mo>
           ( 
         </mo> 
         <mrow> 
          <mfrac> 
           <mn>
             1 
           </mn> 
           <mrow> 
            <msub> 
             <mi>
               m 
             </mi> 
             <mn>
               1 
             </mn> 
            </msub> 
           </mrow> 
          </mfrac> 
          <msup> 
           <mi>
             ω 
           </mi> 
           <mn>
             2 
           </mn> 
          </msup> 
          <mo>
            − 
          </mo> 
          <mi>
            κ 
          </mi> 
         </mrow> 
         <mo>
           ) 
         </mo> 
        </mrow> 
        <mo>
          + 
        </mo> 
        <mi>
          j 
        </mi> 
        <mrow> 
         <mo>
           ( 
         </mo> 
         <mrow> 
          <mi>
            ω 
          </mi> 
          <mfrac> 
           <mi>
             d 
           </mi> 
           <mrow> 
            <msub> 
             <mi>
               m 
             </mi> 
             <mn>
               1 
             </mn> 
            </msub> 
           </mrow> 
          </mfrac> 
         </mrow> 
         <mo>
           ) 
         </mo> 
        </mrow> 
       </mrow> 
      </mfrac> 
     </mrow> 
    </math>(5)</p>
   <p>where:</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        V 
      </mi> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mi>
         ω 
       </mi> 
       <mo>
         ) 
       </mo> 
      </mrow> 
     </mrow> 
    </math>: Voltage applied</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        F 
      </mi> 
      <mrow> 
       <mo>
         ( 
       </mo> 
       <mi>
         ω 
       </mi> 
       <mo>
         ) 
       </mo> 
      </mrow> 
     </mrow> 
    </math>: Force produced</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mo> 
      </mo> 
      <mi>
        ω 
      </mi> 
     </mrow> 
    </math>: Angular frequency</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         k 
       </mi> 
       <mn>
         1 
       </mn> 
      </msub> 
      <mo>
        , 
      </mo> 
      <msub> 
       <mi>
         m 
       </mi> 
       <mn>
         1 
       </mn> 
      </msub> 
      <mo>
        , 
      </mo> 
      <mi>
        d 
      </mi> 
     </mrow> 
    </math>: Constants related to the device’s stiffness, mass, and damping</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       j 
     </mi> 
    </math>: Imaginary unit</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       κ 
     </mi> 
    </math>: Stiffness-related constant</p>
   <p>Finally, the deformation-induced electric field</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <msub> 
       <mi>
         E 
       </mi> 
       <mrow> 
        <mtext>
          harvested 
        </mtext> 
       </mrow> 
      </msub> 
      <mo>
        = 
      </mo> 
      <mfrac> 
       <mn>
         1 
       </mn> 
       <mn>
         2 
       </mn> 
      </mfrac> 
      <mi>
        C 
      </mi> 
      <msup> 
       <mi>
         V 
       </mi> 
       <mn>
         2 
       </mn> 
      </msup> 
     </mrow> 
    </math> (6)</p>
   <p>This closes the loop from mechanical input to electrical output. 
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mi>
       V 
     </mi> 
    </math> here is the result of solving Equations (1)-(3) and computing the potential difference across the piezoelectric material. To provide a practical context for the open-circuit voltage values computed in Equation (6), it is instructive to estimate the output power delivered to a realistic load resistance. The instantaneous power P can be calculated using the relationship:</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        P 
      </mi> 
      <mo>
        = 
      </mo> 
      <mfrac> 
       <mrow> 
        <msup> 
         <mi>
           V 
         </mi> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
       <mi>
         R 
       </mi> 
      </mfrac> 
     </mrow> 
    </math> (7)</p>
   <p>where V is the output voltage across the load and R is the external load resistance. Assuming an RMS voltage output of V = 5 V (as a representative value observed in simulations) and a matched resistive load of R = 100 kΩ, the corresponding power output is:</p>
   <p>
    <math xmlns="http://www.w3.org/1998/Math/MathML"> <mrow> 
      <mi>
        P 
      </mi> 
      <mo>
        = 
      </mo> 
      <mfrac> 
       <mrow> 
        <msup> 
         <mn>
           5 
         </mn> 
         <mn>
           2 
         </mn> 
        </msup> 
       </mrow> 
       <mrow> 
        <mn>
          100 
        </mn> 
        <mo>
          × 
        </mo> 
        <msup> 
         <mrow> 
          <mn>
            10 
          </mn> 
         </mrow> 
         <mn>
           3 
         </mn> 
        </msup> 
       </mrow> 
      </mfrac> 
      <mo>
        = 
      </mo> 
      <mfrac> 
       <mrow> 
        <mn>
          25 
        </mn> 
       </mrow> 
       <mrow> 
        <mn>
          100000 
        </mn> 
       </mrow> 
      </mfrac> 
      <mo>
        = 
      </mo> 
      <mn>
        0.25 
      </mn> 
      <mtext>
        mW 
      </mtext> 
     </mrow> 
    </math> (8)</p>
   <p>This value illustrates that even at modest voltage levels, the piezoelectric harvesting system can yield sub-milliwatt power levels, which are sufficient for low-power sensor nodes or intermittent data transmission applications. These values are expected to vary depending on excitation frequency, amplitude, and structural optimization <xref ref-type="bibr" rid="scirp.144473-15">
     [15]
    </xref>.</p>
   <sec id="s3_1">
    <title>3.1. Design Framework</title>
    <p>The design framework for this study involves vibration modal analysis and harmonic response analysis of the UAV structure, specifically focusing on the wing components. The first step in the design process is to identify high-amplitude vibration zones on the UAV structure. These zones are critical for maximizing energy harvesting efficiency. Finite Element Method (FEM) simulations are used to model the UAV structure and simulate its vibrational behavior under typical flight conditions.</p>
    <p>The main objectives of this study are identifying the optimal location for the piezoelectric (PZT) crystal to maximize power harvesting efficiency, and to calculate the voltage output under varying operational velocities <xref ref-type="bibr" rid="scirp.144473-16">
      [16]
     </xref>.</p>
    <p>The 3D model of the DJI Mini 3 Pro was utilized as a reference for the reverse engineering process. First, the arm and propeller were modelling as shown in <xref ref-type="fig" rid="fig1">
      Figure 1
     </xref>.</p>
    <fig id="fig1" position="float">
     <label>Figure 1</label>
     <caption>
      <title>Figure 1. Arm and propeller after modelling.</title>
     </caption>
     <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId73.jpeg?20250730030640" />
    </fig>
   </sec>
   <sec id="s3_2">
    <title>3.2. Material Properties</title>
    <p>ABS is a thermoplastic polymer composed of three monomers: acrylonitrile, butadiene, and styrene. The main material used in arm and propeller is Plastic ABS (high impact), <xref ref-type="table" rid="table1">
      Table 1
     </xref> shows the arm and propeller properties.</p>
    <table-wrap id="table1">
     <label>
      <xref ref-type="table" rid="table1">
       Table 1
      </xref></label>
     <caption>
      <title>
       <xref ref-type="bibr" rid="scirp.144473-"></xref>Table 1. Arm and propeller ABS properties.</title>
     </caption>
     <table class="MsoTableGrid custom-table" border="0" cellspacing="0" cellpadding="0"> 
      <tr> 
       <td class="custom-bottom-td aleft" width="51.45%"><p style="text-align:left">Properties</p></td> 
       <td class="custom-bottom-td aleft" width="24.03%"><p style="text-align:left">Values</p></td> 
       <td class="custom-bottom-td aleft" width="24.51%"><p style="text-align:left">Unit</p></td> 
      </tr> 
      <tr> 
       <td class="custom-top-td aleft" width="51.45%"><p style="text-align:left">Density</p></td> 
       <td class="custom-top-td aleft" width="24.03%"><p style="text-align:left">1030</p></td> 
       <td class="custom-top-td aleft" width="24.51%"><p style="text-align:left">Kg∙m<sup>−</sup><sup>3</sup></p></td> 
      </tr> 
      <tr> 
       <td class="aleft" width="51.45%"><p style="text-align:left">Thermal expansion coefficient</p></td> 
       <td class="aleft" width="24.03%"><p style="text-align:left">180E−10</p></td> 
       <td class="aleft" width="24.51%"><p style="text-align:left">C<sup>−</sup><sup>1</sup></p></td> 
      </tr> 
      <tr> 
       <td class="aleft" width="51.45%"><p style="text-align:left">Young’s modulus</p></td> 
       <td class="aleft" width="24.03%"><p style="text-align:left">2.09E+10</p></td> 
       <td class="aleft" width="24.51%"><p style="text-align:left">Pa</p></td> 
      </tr> 
      <tr> 
       <td class="aleft" width="51.45%"><p style="text-align:left">Bulk modulus</p></td> 
       <td class="aleft" width="24.03%"><p style="text-align:left">3.8E+9</p></td> 
       <td class="aleft" width="24.51%"><p style="text-align:left">Pa</p></td> 
      </tr> 
      <tr> 
       <td class="aleft" width="51.45%"><p style="text-align:left">Thermal conductivity</p></td> 
       <td class="aleft" width="24.03%"><p style="text-align:left">0.199</p></td> 
       <td class="aleft" width="24.51%"><p style="text-align:left">Wm<sup>−</sup><sup>1</sup>∙C<sup>−</sup><sup>1</sup></p></td> 
      </tr> 
      <tr> 
       <td class="aleft" width="51.45%"><p style="text-align:left">Heat constant pressure</p></td> 
       <td class="aleft" width="24.03%"><p style="text-align:left">1400</p></td> 
       <td class="aleft" width="24.51%"><p style="text-align:left">J Kg<sup>−</sup><sup>1</sup>∙C<sup>−</sup><sup>1</sup></p></td> 
      </tr> 
      <tr> 
       <td class="aleft" width="51.45%"><p style="text-align:left">Resistivity</p></td> 
       <td class="aleft" width="24.03%"><p style="text-align:left">9.9E+13</p></td> 
       <td class="aleft" width="24.51%"><p style="text-align:left">Ohm∙Ω</p></td> 
      </tr> 
      <tr> 
       <td class="aleft" width="51.45%"><p style="text-align:left">Climate change CO<sub>2</sub></p></td> 
       <td class="aleft" width="24.03%"><p style="text-align:left">3.66</p></td> 
       <td class="aleft" width="24.51%"><p style="text-align:left">Kg/Kg</p></td> 
      </tr> 
      <tr> 
       <td class="aleft" width="51.45%"><p style="text-align:left">Embodied energy</p></td> 
       <td class="aleft" width="24.03%"><p style="text-align:left">9.65E+7</p></td> 
       <td class="aleft" width="24.51%"><p style="text-align:left">J/Kg</p></td> 
      </tr> 
      <tr> 
       <td class="aleft" width="51.45%"><p style="text-align:left">Recycle</p></td> 
       <td class="aleft" width="24.03%"><p style="text-align:left">1</p></td> 
       <td class="aleft" width="24.51%"><p style="text-align:left">Time</p></td> 
      </tr> 
     </table>
    </table-wrap>
    <p>Piezoelectric Crystal (PZT-5H) is softer coefficients, higher dielectric constant, higher sensitivity to stress or electric field, lower mechanical lossy PZT Crystal Dimensions: Assumed to be lightweight (0.25 g). <xref ref-type="fig" rid="fig2">
      Figure 2
     </xref> illustrates PZT dimensional.</p>
    <fig id="fig2" position="float">
     <label>Figure 2</label>
     <caption>
      <title>Figure 2. Piezoelectric crystal size.</title>
     </caption>
     <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId74.jpeg?20250730030641" />
    </fig>
    <p>The successful integration of structural steel fins with UAV arms represents a pivotal advancement in enhancing piezoelectric energy harvesting within quadcopter systems. Previous modeling efforts have primarily focused on characterizing PZT-5H piezoelectric patches through finite element analysis (FEA). Building on this foundation, the addition of structural steel fins functions as a mechanical amplifier (see <xref ref-type="table" rid="table2">
      Table 2
     </xref>), strategically modifying the vibrational environment to which the piezoelectric patches are exposed. The fins are designed in accordance with civil engineering fin-plate principles.</p>
    <p>Peak voltage output occurs near the resonance frequencies, underscoring the synergistic interaction between aerodynamic dynamic excitation and the piezoelectric response. Furthermore, the aerodynamic drag penalty resulting from fin integration remains minimal, with an increase of less than 5%.</p>
    <p>Overall, the integrated CFD-FEA modeling approach confirms that the structural steel fins effectively serve as mechanical intermediaries, amplifying the vibrational energy that can be harvested by the piezoelectric patches. This interdisciplinary innovation—bridging civil structural design with aerial energy harvesting—marks a significant step toward self-powered UAV systems with enhanced endurance and operational autonomy.</p>
    <table-wrap id="table2">
     <label>
      <xref ref-type="table" rid="table2">
       Table 2
      </xref></label>
     <caption>
      <title>
       <xref ref-type="bibr" rid="scirp.144473-"></xref>Table 2. Structural steel properties.</title>
     </caption>
     <table class="MsoTableGrid custom-table" border="0" cellspacing="0" cellpadding="0"> 
      <tr> 
       <td class="custom-bottom-td acenter" width="42.58%"><p style="text-align:center">Properties</p></td> 
       <td class="custom-bottom-td acenter" width="17.46%"><p style="text-align:center">Values</p></td> 
       <td class="custom-bottom-td acenter" width="17.01%"><p style="text-align:center">Unit</p></td> 
      </tr> 
      <tr> 
       <td class="custom-top-td acenter" width="42.58%"><p style="text-align:center">Density</p></td> 
       <td class="custom-top-td acenter" width="17.46%"><p style="text-align:center">7850</p></td> 
       <td class="custom-top-td acenter" width="17.01%"><p style="text-align:center">Kg∙m<sup>−</sup><sup>3</sup></p></td> 
      </tr> 
      <tr> 
       <td class="acenter" width="42.58%"><p style="text-align:center">Thermal expansion coefficient</p></td> 
       <td class="acenter" width="17.46%"><p style="text-align:center">1.2E−5</p></td> 
       <td class="acenter" width="17.01%"><p style="text-align:center">C<sup>−</sup><sup>1</sup></p></td> 
      </tr> 
      <tr> 
       <td class="acenter" width="42.58%"><p style="text-align:center">Young’s modulus</p></td> 
       <td class="acenter" width="17.46%"><p style="text-align:center">2E+11</p></td> 
       <td class="acenter" width="17.01%"><p style="text-align:center">Pa</p></td> 
      </tr> 
      <tr> 
       <td class="acenter" width="42.58%"><p style="text-align:center">Bulk modulus</p></td> 
       <td class="acenter" width="17.46%"><p style="text-align:center">1.66E+11</p></td> 
       <td class="acenter" width="17.01%"><p style="text-align:center">Pa</p></td> 
      </tr> 
      <tr> 
       <td class="acenter" width="42.58%"><p style="text-align:center">Thermal conductivity</p></td> 
       <td class="acenter" width="17.46%"><p style="text-align:center">60.5</p></td> 
       <td class="acenter" width="17.01%"><p style="text-align:center">Wm<sup>−</sup><sup>1</sup>∙C<sup>−</sup><sup>1</sup></p></td> 
      </tr> 
      <tr> 
       <td class="acenter" width="42.58%"><p style="text-align:center">Heat constant pressure</p></td> 
       <td class="acenter" width="17.46%"><p style="text-align:center">434</p></td> 
       <td class="acenter" width="17.01%"><p style="text-align:center">J∙Kg<sup>−</sup><sup>1</sup>∙C<sup>−</sup><sup>1</sup></p></td> 
      </tr> 
      <tr> 
       <td class="acenter" width="42.58%"><p style="text-align:center">Resistivity</p></td> 
       <td class="acenter" width="17.46%"><p style="text-align:center">1.7E−7</p></td> 
       <td class="acenter" width="17.01%"><p style="text-align:center">Ohm∙Ω</p></td> 
      </tr> 
      <tr> 
       <td class="acenter" width="42.58%"><p style="text-align:center">Tensile Yield strength</p></td> 
       <td class="acenter" width="17.46%"><p style="text-align:center">2.5E+8</p></td> 
       <td class="acenter" width="17.01%"><p style="text-align:center">Pa</p></td> 
      </tr> 
      <tr> 
       <td class="acenter" width="42.58%"><p style="text-align:center">Compressive Yield strength</p></td> 
       <td class="acenter" width="17.46%"><p style="text-align:center">2.5E+8</p></td> 
       <td class="acenter" width="17.01%"><p style="text-align:center">Pa</p></td> 
      </tr> 
      <tr> 
       <td class="acenter" width="42.58%"><p style="text-align:center">Tensile ultimate strength</p></td> 
       <td class="acenter" width="17.46%"><p style="text-align:center">4.6E+8</p></td> 
       <td class="acenter" width="17.01%"><p style="text-align:center">Pa</p></td> 
      </tr> 
     </table>
    </table-wrap>
   </sec>
   <sec id="s3_3">
    <title>
     <xref ref-type="bibr" rid="scirp.144473-"></xref>3.3. Simulation Parameters</title>
    <p>Computational fluid dynamics (CFD) Analysis: Conducted using SolidWorks and ANSYS under steady-state conditions. Damping: Only structural damping (assumed to be 1%) was considered during modal and harmonic analyses.</p>
    <p>CFD Loads: Applied only in the initial stages for determining the PZT crystal’s location. Subsequent stages considered rotational speed instead.</p>
    <p>Mesh Refinement: Default meshing parameters were used, except for the PZT crystal, where mesh refinement was adjusted to its size.</p>
    <p>Analysis: Conducted using SolidWorks and ANSYS under steady-state conditions.</p>
    <p>Damping: Only structural damping (assumed to be 1%) was considered during modal and harmonic analyses.</p>
    <p>CFD Loads: Applied only in the initial stages for determining the PZT crystal’s location. Subsequent stages considered rotational speed instead.</p>
    <p>Mesh Refinement: Default meshing parameters were used, except for the PZT crystal, where mesh refinement was adjusted to its size. As shown in <xref ref-type="fig" rid="fig3">
      Figure 3
     </xref>.</p>
    <fig id="fig3" position="float">
     <label>Figure 3</label>
     <caption>
      <title>Figure 3. Default mesh size.</title>
     </caption>
     <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId75.jpeg?20250730030644" />
    </fig>
    <p>The velocities range are used to base 1 on experimental data were obtained for some users <xref ref-type="bibr" rid="scirp.144473-17">
      [17]
     </xref>. A 3D model of the DJI Mini 3 Pro was developed as a baseline for analysis and simulation.</p>
    <p>(CFD) analysis was performed to evaluate how aerodynamic loads influence the natural frequencies of the arm and propeller. Due to the symmetric layout of the quadcopter and the largely independent vibrational behavior of each rotor–arm assembly, analyzing a single arm provides representative insights applicable to all four arms. Preliminary simulations confirmed that each arm exhibits nearly identical dynamic responses; therefore, a single arm was selected for detailed design and analysis to reduce computational complexity without compromising result accuracy. See <xref ref-type="fig" rid="fig4">
      Figure 4
     </xref>.</p>
    <fig id="fig4" position="float">
     <label>Figure 4</label>
     <caption>
      <title>Figure 4. ANSYS setup.</title>
     </caption>
     <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId76.jpeg?20250730030644" />
    </fig>
    <p>The PZT crystal was placed at the identified optimal location. A structural steel fin was attached to the PZT, ensuring it contacted the propeller for enhanced vibration transmission. <xref ref-type="fig" rid="fig4">
      Figure 4
     </xref> and <xref ref-type="fig" rid="fig5">
      Figure 5
     </xref> show the setup and location for PZT.</p>
    <fig id="fig5" position="float">
     <label>Figure 5</label>
     <caption>
      <title>Figure 5. The best location for Piezoelectric.</title>
     </caption>
     <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId77.jpeg?20250730030644" />
    </fig>
    <p>The following <xref ref-type="fig" rid="figFigures 6-9">
      Figures 6-9
     </xref> show the effect of aerodynamic on design model.</p>
    <p>Boundary conditions corresponding to varying velocities were applied, and the resulting voltage outputs were measured as shown in <xref ref-type="fig" rid="fig10">
      Figure 10
     </xref>.</p>
   </sec>
  </sec><sec id="s4">
   <title>
    <xref ref-type="bibr" rid="scirp.144473-"></xref>4. Results</title>
   <p>The analysis is focused on one arm of the quadcopter, assuming the results are</p>
   <fig id="fig6" position="float">
    <label>Figure 6</label>
    <caption>
     <title>Figure 6. Pressure distribution based on the applied velocity.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId78.jpeg?20250730030646" />
   </fig>
   <fig id="fig7" position="float">
    <label>Figure 7</label>
    <caption>
     <title>Figure 7. The resultant frequency based on the loads (136.58 Hz).</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId79.jpeg?20250730030646" />
   </fig>
   <fig id="fig8" position="float">
    <label>Figure 8</label>
    <caption>
     <title>Figure 8. The frequency effect on propeller and the arm (without damping).</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId80.jpeg?20250730030646" />
   </fig>
   <fig id="fig9" position="float">
    <label>Figure 9</label>
    <caption>
     <title>Figure 9. The frequency effect on propellers and the arm (with damping).</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId81.jpeg?20250730030647" />
   </fig>
   <fig id="fig10" position="float">
    <label>Figure 10</label>
    <caption>
     <title>Figure 10. Boundary conditions.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId82.jpeg?20250730030647" />
   </fig>
   <p>representative of all arms. Based on the simulation, the piezoelectric crystal (PZT) is not significantly influenced by the natural frequencies of the propeller and arm. To address this, a steel fin was attached to the propeller, introducing additional vibrations generated by the rotational motion of the propeller. This setup ensures that the PZT crystal experiences sufficient excitation to generate electricity (see <xref ref-type="table" rid="table3">
     Table 3
    </xref>).</p>
   <table-wrap id="table3">
    <label>
     <xref ref-type="table" rid="table3">
      Table 3
     </xref></label>
    <caption>
     <title>
      <xref ref-type="bibr" rid="scirp.144473-"></xref>Table 3. Frequency generated from different velocity.</title>
    </caption>
    <table class="MsoTableGrid custom-table" border="0" cellspacing="0" cellpadding="0"> 
     <tr> 
      <td class="custom-bottom-td acenter" width="37.65%"><p style="text-align:center">Velocity (RPM)</p></td> 
      <td class="custom-bottom-td acenter" width="37.84%"><p style="text-align:center">Frequency (Hz)</p></td> 
      <td class="custom-bottom-td acenter" width="24.51%"><p style="text-align:center">Load (N)</p></td> 
     </tr> 
     <tr> 
      <td class="custom-top-td acenter" width="37.65%"><p style="text-align:center">3000</p></td> 
      <td class="custom-top-td acenter" width="37.84%"><p style="text-align:center">50.00</p></td> 
      <td class="custom-top-td acenter" width="24.51%"><p style="text-align:center">2.96</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="37.65%"><p style="text-align:center">4000</p></td> 
      <td class="acenter" width="37.84%"><p style="text-align:center">66.67</p></td> 
      <td class="acenter" width="24.51%"><p style="text-align:center">5.26</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="37.65%"><p style="text-align:center">5000</p></td> 
      <td class="acenter" width="37.84%"><p style="text-align:center">83.33</p></td> 
      <td class="acenter" width="24.51%"><p style="text-align:center">8.22</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="37.65%"><p style="text-align:center">6000</p></td> 
      <td class="acenter" width="37.84%"><p style="text-align:center">100.00</p></td> 
      <td class="acenter" width="24.51%"><p style="text-align:center">11.84</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="37.65%"><p style="text-align:center">7000</p></td> 
      <td class="acenter" width="37.84%"><p style="text-align:center">116.67</p></td> 
      <td class="acenter" width="24.51%"><p style="text-align:center">16.12</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="37.65%"><p style="text-align:center">8000</p></td> 
      <td class="acenter" width="37.84%"><p style="text-align:center">133.33</p></td> 
      <td class="acenter" width="24.51%"><p style="text-align:center">21.06</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="37.65%"><p style="text-align:center">9000</p></td> 
      <td class="acenter" width="37.84%"><p style="text-align:center">150.00</p></td> 
      <td class="acenter" width="24.51%"><p style="text-align:center">26.65</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="37.65%"><p style="text-align:center">10,000</p></td> 
      <td class="acenter" width="37.84%"><p style="text-align:center">166.67</p></td> 
      <td class="acenter" width="24.51%"><p style="text-align:center">32.9</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="37.65%"><p style="text-align:center">11,000</p></td> 
      <td class="acenter" width="37.84%"><p style="text-align:center">183.33</p></td> 
      <td class="acenter" width="24.51%"><p style="text-align:center">39.81</p></td> 
     </tr> 
     <tr> 
      <td class="acenter" width="37.65%"><p style="text-align:center">12,000</p></td> 
      <td class="acenter" width="37.84%"><p style="text-align:center">200.00</p></td> 
      <td class="acenter" width="24.51%"><p style="text-align:center">47.37</p></td> 
     </tr> 
    </table>
   </table-wrap>
   <fig id="fig11" position="float">
    <label>Figure 11</label>
    <caption>
     <title>Figure 11. Output voltage.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId83.jpeg?20250730030647" />
   </fig>
   <p>
    <xref ref-type="bibr" rid="scirp.144473-"></xref>The increase in voltage output with frequency can be attributed to the piezoelectric effect, where higher excitation frequencies lead to greater strain rates in the material (see <xref ref-type="fig" rid="fig11">
     Figure 11
    </xref>). As the frequency approaches the system’s resonance, the mechanical deformation intensifies, resulting in a significant increase in voltage generation. This behavior is expected due to the energy amplification effect near resonance, where the piezoelectric material experiences enhanced mechanical-to-electrical energy conversion. Additionally, higher frequencies induce a more rapid charge accumulation in the piezoelectric domains, further contributing to the elevated voltage output.</p>
   <p>Different voltage values were generated (1.54, 2.7, 4.2, 6.1, 8.2, 10.8, 13.7, 16.8, 20.4, 24.3 V) corresponds to the predominant voltage across the piezoelectric layer, it can be considered a reliable estimate of the overall output voltage. Although the maximum and minimum values indicate localized stress concentrations, the voltage observed in the dominant region offers a more practical and representative measure of the device’s performance. <xref ref-type="fig" rid="fig12">
     Figure 12
    </xref> shows the voltage generated across different frequency, and <xref ref-type="fig" rid="fig13">
     Figure 13
    </xref> illustrates the output voltage from different velocity.</p>
   <p>The main contributions of this work are as follows:</p>
   <p>Development of a validated 3D simulation model of UAV arm dynamics incorporating piezoelectric energy harvesting components, in addition, Identification of optimal piezoelectric placement zones through modal and harmonic analysis by Demonstration of voltage generation under realistic UAV operational scenarios.</p>
   <fig id="fig12" position="float">
    <label>Figure 12</label>
    <caption>
     <title>Figure 12. Voltage generated from different frequency.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId84.jpeg?20250730030647" />
   </fig>
   <fig id="fig13" position="float">
    <label>Figure 13</label>
    <caption>
     <title>Figure 13. Voltage generated from different velocity.</title>
    </caption>
    <graphic mimetype="image" position="float" xlink:type="simple" xlink:href="https://html.scirp.org/file/8104779-rId85.jpeg?20250730030646" />
   </fig>
   <p>These results validate the feasibility of using piezoelectric materials as auxiliary power sources in UAV systems, paving the way for more energy-efficient and sustainable aerial platforms. Additionally, the study highlights that energy harvesting performance is significantly improved near resonant frequencies, suggesting that operational tuning may serve as a viable performance optimization strategy.</p>
  </sec><sec id="s5">
   <title>5. Conclusion</title>
   <p>This study demonstrates a comprehensive modeling and simulation-based investigation into the integration of piezoelectric energy harvesting systems within quadrotor unmanned aerial vehicles (UAVs). By employing high-fidelity Computational Fluid Dynamics (CFD) and Finite Element Analysis (FEA), the research demonstrates that strategically positioning PZT-5H piezoelectric crystals near zones of high vibrational stress can yield efficient electrical outputs, ranging from 1.54 V to 24.3 V across varying rotor speeds. Furthermore, the inclusion of a structural steel fin significantly enhanced vibration transmission, confirming the viability of mechanical-to-electrical energy conversion even under off-resonance conditions.</p>
   <sec id="s5_1">
    <title>Future Work</title>
    <p>Looking ahead, further research should be investigated: Hybrid energy systems combining piezoelectric, solar, and electromagnetic harvesting for enhanced energy resilience. Experimental validation through wind tunnel and in-flight testing to capture dynamic, multi-axial vibrational effects. Long-term durability assessments of piezoelectric materials under extended UAV operation Integration with control systems, where harvested energy could support low-power sensors or emergency backup functions. By advancing the understanding and practical application of energy harvesting technologies in UAVs, this research contributes meaningfully to the development of the next generation of autonomous aerial systems with extended endurance and reduced dependence on conventional batteries.</p>
   </sec>
  </sec>
 </body><back>
  <ref-list>
   <title>References</title>
   <ref id="scirp.144473-ref1">
    <label>1</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Zhang, Y., et al. (2020) Optimal Piezoelectric Energy Harvesting in UAV Wings. Journal of Energy Conversion, 2, 697-710. 
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref2">
    <label>2</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Marini, L., et al. (2021) CFD-FEA Coupling for UAV Vibration Analysis. IEEE Transactions on Aerospace, 5, No. 6.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref3">
    <label>3</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Liu, H., Fu, H., et al. (2022) Hybrid Energy Harvesting Technology: From Materials, Structure Design, System Integration to Applications. Renewable and Sustainable Energy Reviews, 137, Article ID: 110473. &gt;https://doi.org/10.1016/j.rser.2020.110473
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref4">
    <label>4</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Chen, X. and Zhao, Y. (2019) Vibration-Based Energy Harvesting in Lightweight UAV Structures. Smart Materials and Structures, 28, Article ID: 095012.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref5">
    <label>5</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Kumar, R., et al. (2020) Design and Analysis of Piezoelectric Harvesters in Flapping-wing UAVs. Sensors and Actuators A: Physical, 305, Article ID: 111939.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref6">
    <label>6</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Tang, L. and Yang, Y. (2018) Advances in Nonlinear Piezoelectric Energy Harvesting for UAV Applications. Applied Physics Reviews, 5, Article ID: 031106.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref7">
    <label>7</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Silva, M.A., et al. (2021) Aeroelastic Modeling and Energy Harvesting in Flexible UAV Wings. Aerospace Science and Technology, 117, Article ID: 106952.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref8">
    <label>8</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Ahmed, Z. and Lin, J. (2022) Multi-Source Energy Harvesting for Self-Powered UAV Navigation. Energy Reports, 8, 623-631.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref9">
    <label>9</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Patel, S. and Wang, T. (2020) Structural Health Monitoring and Energy Harvesting Integration in UAVs. Journal of Intelligent Material Systems and Structures, 31, 1851-1864.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref10">
    <label>10</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Zhou, M., et al. (2021) Experimental Validation of Piezoelectric Harvesters on Small UAV Airframes. Mechanical Systems and Signal Processing, 154, Article ID: 107537.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref11">
    <label>11</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Huang, R. and Sun, C.T. (2020) Optimization of Energy Harvesting Layout in UAV Structures Using Topology Design. Journal of Sound and Vibration, 481, Article ID: 115452.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref12">
    <label>12</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Yilmaz, T., et al. (2021) Design of a Micro Wind Energy Harvester for UAV Applications. Renewable Energy, 179, 1762-1771.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref13">
    <label>13</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Wang, B., et al. (2018) Coupled Electromechanical Modeling of Piezoelectric Energy Harvesters for UAV Flight Dynamics. Journal of Physics: Conference Series, 1052, Article ID: 012051.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref14">
    <label>14</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Lee, D. and Kim, J. (2019) Adaptive Piezoelectric Energy Harvesting under UAV Flight Conditions. Energy Harvesting and Systems, 6, 37-46.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref15">
    <label>15</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Oliveira, J.P., et al. (2022) A Review on Energy Harvesting Solutions for UAVs: Technologies, Challenges, and Future Trends. Progress in Aerospace Sciences, 137, Article ID: 100789.
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref16">
    <label>16</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Moelyadi, M.A., et al. (2024) Modelling and Simulation of Fluid-Structure Interaction for Piezoelectric Energy Harvesting in Simplified HALE UAV Wings. In: Akmeliawati, R., Harvey, D., Sergiienko, N., Yang, L.J. and Park, H.C., Eds., Proceedings of the 19th International Conference on Intelligent Unmanned Systems. ICIUS 2023, Springer, 125-135. &gt;https://doi.org/10.1007/978-981-97-6591-1_12
    </mixed-citation>
   </ref>
   <ref id="scirp.144473-ref17">
    <label>17</label>
    <mixed-citation publication-type="other" xlink:type="simple">
     Koko, H. and Ibrahim, G. (2024) A Computational Framework to Optimize Piezoelectric Power Harvesting from a UAV Wing. Journal of Aviation Technology and Engineering, 13, Article 1. &gt;https://doi.org/10.7771/2159-6670.1292
    </mixed-citation>
   </ref>
  </ref-list>
 </back>
</article>