<?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">WJET</journal-id><journal-title-group><journal-title>World Journal of Engineering and Technology</journal-title></journal-title-group><issn pub-type="epub">2331-4222</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/wjet.2015.34016</article-id><article-id pub-id-type="publisher-id">WJET-60874</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Chemistry&amp;Materials Science</subject><subject> Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Parameters Optimization for Piezoelectric Harvesting Energy from Pavement Based on Taguchi’s Orthogonal Experiment Design
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>hunhua</surname><given-names>Sun</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>Hongbing</surname><given-names>Wang</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>Guangqing</surname><given-names>Shang</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>Jianhong</surname><given-names>Du</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Mechanic and Electronic Engineering, Suzhou Vocational University, Suzhou, China</addr-line></aff><pub-date pub-type="epub"><day>03</day><month>11</month><year>2015</year></pub-date><volume>03</volume><issue>04</issue><fpage>149</fpage><lpage>157</lpage><history><date date-type="received"><day>5</day>	<month>August</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>31</month>	<year>October</year>	</date><date date-type="accepted"><day>3</day>	<month>November</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; 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><p>
 
 
  To effectively harvest vibration energy from pavement without affecting driving comfort and safety, parameter optimization was done with the orthogonal experiment design and the finite element analysis. L
  <sub>16(</sub>4
  <sup>4</sup>) Taguchi’s orthogonal experiments were carried out with planted depth, PZT material, PZT diameter and thickness as optimization parameters and with open voltage and pavement displacement as optimization objectives. The experiment results were obtained via the finite element method. By using range analysis method, the dominance degree of the influencing factors and the optimum condition was obtained for the two objectives, respectively. Further, the multi-objective optimization was performed based on a weight grade method. The combined optimum conditions in order of their dominance degree are PZT diameter 35 mm, PZT thickness 6 mm, planted depth 50 mm and material PZT4. The validity of optimization scheme was confirmed.
 
</p></abstract><kwd-group><kwd>Piezoelectric Energy Harvesting</kwd><kwd> Pavement</kwd><kwd> Taguchi’s Orthogonal Experiment</kwd><kwd> Finite Element  Method</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Piezoelectric harvesting energy from pavement is a new process of capturing the wasted pavement vibration and transferring it into electricity to store it for later use or to power electronic sensors [<xref ref-type="bibr" rid="scirp.60874-ref1">1</xref>] . Since the 1990s, the technology has attracted more and more attention because of its unique advantages such as high energy convert efficiency, simple piezoelectric device, and applicability for miniaturization [<xref ref-type="bibr" rid="scirp.60874-ref2">2</xref>] . Some studies have been carried out on the effect of cymbal piezoelectric transducers of structural parameters on potential electrical energy and coupling with pavement [<xref ref-type="bibr" rid="scirp.60874-ref3">3</xref>] . Electromechanical analysis of piezoelectric harvesting unit from road vibration was also done with the finite element method (FEM) [<xref ref-type="bibr" rid="scirp.60874-ref4">4</xref>] . The effect of traffic flow on characteristics of piezoelectric harvesting unit from pavement was analyzed [<xref ref-type="bibr" rid="scirp.60874-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.60874-ref6">6</xref>] . However, little application of the technology has been realized until now [<xref ref-type="bibr" rid="scirp.60874-ref7">7</xref>] . The main reason consists in the relative small current and low power compared to some other harvesting energy techniques such as electromagnetic and photovoltaic. Therefore, it is of significance for the practical application of piezoelectric energy harvesting to improve process efficiency through parameter optimization.</p><p>It has been well known that Taguchi’s approach is an efficient way to correlate various factors for parameter optimization and to identify the most influencing factors based on the experimental design and analysis method with good consistency and reproducibility. To our knowledge, the studies on the application of Taguchi experimental design for the purpose of high efficiency of piezoelectric harvesting energy are few [<xref ref-type="bibr" rid="scirp.60874-ref8">8</xref>] . And the application of the method to parameter optimization of piezoelectric harvesting device from pavement is not yet reported.</p><p>In this paper, the optimization of piezoelectric harvesting device from pavement was studied based on Taguchi experimental design for high efficiency of harvesting energy and good coupling with pavement. The parameters, including embedded depth, thickness and diameter of PZT disk, and material type of PZT disk, were involved. Firstly, the influence of factors on piezoelectric harvesting device from pavement was investigated. Then, the orthogonal experiments were carried out and analyzed with the output voltage and relative displacement ratio as optimization objectives. Dominance degree of the factors and the optimum conditions for piezoelectric harvesting energy from pavement were studied.</p></sec><sec id="s2"><title>2. Optimization Objectives and Parameters</title><p><xref ref-type="fig" rid="fig1">Figure 1</xref> shows the harvester-pavement coupling system, on which a piezoelectric harvester is used to scavenge pavement vibration caused by vehicle traveling and transfer it into electricity. A successful harvester used in the system should have the abilities that it can convert mechanical energy into electrical power as more as possible and maintain the safety and comfort performance of the pavement. Thus, two basic performances are required for the harvester as following: 1) high efficiency of harvesting energy; 2) good coupling with pavement.</p><p>The efficiency of a harvester is always represented by the electromechanical coupling factor k and the energy transmission coefficient<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x6.png" xlink:type="simple"/></inline-formula>. However, it is difficult to calculate k and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x7.png" xlink:type="simple"/></inline-formula> in the complicated harvester-pavement coupling system. Instead of them, the voltage at open circuit, which is directly related to the energy convert efficiency, is adopted to evaluate the energy harvesting ability in this paper.</p><p>The output voltage of PZT disk under the applied stress can be obtained from the following equation in [<xref ref-type="bibr" rid="scirp.60874-ref9">9</xref>] :</p><disp-formula id="scirp.60874-formula67"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1560199x8.png"  xlink:type="simple"/></disp-formula><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Harvester-pavement coupling system</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1560199x9.png"/></fig><disp-formula id="scirp.60874-formula68"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1560199x10.png"  xlink:type="simple"/></disp-formula><p>where, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x11.png" xlink:type="simple"/></inline-formula>is piezoelectric strain constant tensor; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x12.png" xlink:type="simple"/></inline-formula>is the piezoelectric voltage constant of PZT; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x13.png" xlink:type="simple"/></inline-formula>is the relative dielectric constant of PZT; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x14.png" xlink:type="simple"/></inline-formula>is the dielectric constant of vacuum; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x15.png" xlink:type="simple"/></inline-formula>is the stress tensor; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x16.png" xlink:type="simple"/></inline-formula>is thickness of PZT disk; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x17.png" xlink:type="simple"/></inline-formula>is the electric potential(the difference of the electric potential is referred as voltage) along the direction of 3 caused by polarization changing, which can be obtained via FEM.</p><p>Meanwhile, the good coupling with pavement reflects that there is no much change of pavement displacement with or without the harvester under pavement while a vehicle passes through. Here, the relative deformation ratio p is adopted to express the coupling property.</p><disp-formula id="scirp.60874-formula69"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1560199x18.png"  xlink:type="simple"/></disp-formula><p>where, d<sub>0</sub> and d<sub>w</sub> represent the pavement displacements without and with the harvester under the pavement, respectively.</p><p>If the relative deformation ratio p is less, it indicates that the harvester has good coupling with pavement. The pavement displacement d<sub>0</sub> and d<sub>w</sub> can be obtained according to the maximum displacement value of the pavement via FEM. Therefore, the output voltage <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x19.png" xlink:type="simple"/></inline-formula> and the relative deformation ratio p are taken to be the optimization objectives.</p><p>Parameters, which affect the harvesting energy efficiency and coupling with pavement, deal with three types. One is related to the harvester’s ceramic material, another relates to the harvester structure and the other is the embedded depth under pavement. If the cylindrical harvester is adopted in the harvester-pavement piezoelectric system, the structural parameters of the harvester include the diameter and thickness of PZT disk．Therefore, the detailed parameters including the embedded depth of the harvester under pavement, PZT material, the diameter and thickness of PZT disk, are chosen to optimize.</p></sec><sec id="s3"><title>3. Parameter Optimization by Taguchi’s Orthogonal Experiment Design</title><sec id="s3_1"><title>3.1. The Taguchi’s Orthogonal Experiment Design</title><p>According to the Taguchi parameter design methodology, a <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x20.png" xlink:type="simple"/></inline-formula> standard orthogonal array with four factors in four levels each in 16 runs is employed. The <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x21.png" xlink:type="simple"/></inline-formula> considering the four above-discussed parameters is shown in <xref ref-type="table" rid="table1">Table 1</xref>. Each row of the orthogonal array represents a specified set of factor levels to be tested.</p><p>The results of open voltage<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x22.png" xlink:type="simple"/></inline-formula>, pavement displacement d<sub>w</sub> and d<sub>0</sub> are obtained via FEM, the relative deformation ratio p is then calculated according to Equation (3). All are shown in <xref ref-type="table" rid="table1">Table 1</xref>. It can be seen that open voltage<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x23.png" xlink:type="simple"/></inline-formula>, pavement displacement d<sub>w</sub> and relative deformation ratio p vary in a range of 53.593 - 140.638 V 0.0648 - 0.0663 mm, and 3.68% - 6.08%, respectively. The open voltage <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x24.png" xlink:type="simple"/></inline-formula> varies in wide range and the relative deformation ratio p has less change under the conditions in orthogonal experiment design. Here, the pavement displacement d<sub>0 </sub>without the harvester is d<sub>0</sub> = 0.0625 mm through the finite element analysis.</p><p>The orthogonal experiment was analyzed by range analysis method for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x25.png" xlink:type="simple"/></inline-formula> and p as shown in <xref ref-type="table" rid="table2">Table 2</xref>.</p></sec><sec id="s3_2"><title>3.2. Single-Objective Optimization</title><p>In <xref ref-type="table" rid="table2">Table 2</xref>, the arithmetical parameter k<sub>i</sub> is the average value of the objective variable, where i denotes the level of the certain factor. R is the range calculated as the difference between the maximum and the minimum values of k<sub>i</sub> for a certain factor. According to the value of R, the dominance degree of each factor can be determined.</p><p><xref ref-type="table" rid="table3">Table 3</xref> shows the analysis results according to single objective of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x26.png" xlink:type="simple"/></inline-formula> and p. It shows that the dominance degree of various parameters as well as the optimum conditions is different for each optimization objective. PZT thickness and diameter are the key factors to the two objectives. Meanwhile, the influence of PZT material is relatively small for all of the objectives. Planted depth is the most important factor for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x27.png" xlink:type="simple"/></inline-formula> while it is not so important for p. Obviously, this is because of the different requirement for the levels of factors aiming at the different optimization objectives.</p><p>Experiments were run according to the chosen optimum conditions. The optimized values are shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>. The voltage objective <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x28.png" xlink:type="simple"/></inline-formula> achieved under optimum conditions is 151.8 V. The optimized value is com-</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Schedule and Results of Taguchi’s orthogonal experiments</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Run</th><th align="center" valign="middle"  colspan="4"  >Factors</th><th align="center" valign="middle"  colspan="3"  >Results</th></tr></thead><tr><td align="center" valign="middle" >1(Planted depth h/mm)</td><td align="center" valign="middle" >2(PZT materials)</td><td align="center" valign="middle" >3(PZT thickness/mm)</td><td align="center" valign="middle" >4(PZT diameter/mm)</td><td align="center" valign="middle" >Open voltage U<sub>3</sub>/V</td><td align="center" valign="middle" >Pavement displacement d<sub>w</sub>/mm</td><td align="center" valign="middle" >Relative deformation ratio p</td></tr><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1(20)</td><td align="center" valign="middle" >1(PZT4)</td><td align="center" valign="middle" >1(2)</td><td align="center" valign="middle" >1(20)</td><td align="center" valign="middle" >53.593</td><td align="center" valign="middle" >0.0663</td><td align="center" valign="middle" >6.08</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >1(20)</td><td align="center" valign="middle" >2(PZT5A)</td><td align="center" valign="middle" >2(4)</td><td align="center" valign="middle" >2(25)</td><td align="center" valign="middle" >87.267</td><td align="center" valign="middle" >0.0655</td><td align="center" valign="middle" >4.80</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >1(20)</td><td align="center" valign="middle" >3(PZT5H)</td><td align="center" valign="middle" >3(5)</td><td align="center" valign="middle" >3(30)</td><td align="center" valign="middle" >62.627</td><td align="center" valign="middle" >0.0654</td><td align="center" valign="middle" >4.64</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >1(20)</td><td align="center" valign="middle" >4(PZT8)</td><td align="center" valign="middle" >4(6)</td><td align="center" valign="middle" >4(35)</td><td align="center" valign="middle" >98.81</td><td align="center" valign="middle" >0.0648</td><td align="center" valign="middle" >3.68</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >2(30)</td><td align="center" valign="middle" >1(PZT4)</td><td align="center" valign="middle" >2(4)</td><td align="center" valign="middle" >3(30)</td><td align="center" valign="middle" >134.121</td><td align="center" valign="middle" >0.0654</td><td align="center" valign="middle" >4.64</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >2(30)</td><td align="center" valign="middle" >2(PZT5A)</td><td align="center" valign="middle" >1(2)</td><td align="center" valign="middle" >4(35)</td><td align="center" valign="middle" >86.819</td><td align="center" valign="middle" >0.0660</td><td align="center" valign="middle" >5.60</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >2(30)</td><td align="center" valign="middle" >3(PZT5H)</td><td align="center" valign="middle" >4(6)</td><td align="center" valign="middle" >1(20)</td><td align="center" valign="middle" >78.95</td><td align="center" valign="middle" >0.0651</td><td align="center" valign="middle" >4.16</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >2(30)</td><td align="center" valign="middle" >4(PZT8)</td><td align="center" valign="middle" >3(5)</td><td align="center" valign="middle" >2(25)</td><td align="center" valign="middle" >86.676</td><td align="center" valign="middle" >0.0654</td><td align="center" valign="middle" >4.64</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >3(40)</td><td align="center" valign="middle" >1(PZT4)</td><td align="center" valign="middle" >3(5)</td><td align="center" valign="middle" >4(35)</td><td align="center" valign="middle" >133.342</td><td align="center" valign="middle" >0.0651</td><td align="center" valign="middle" >4.16</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >3(40)</td><td align="center" valign="middle" >2(PZT5A)</td><td align="center" valign="middle" >4(6)</td><td align="center" valign="middle" >3(30)</td><td align="center" valign="middle" >139.898</td><td align="center" valign="middle" >0.0649</td><td align="center" valign="middle" >3.84</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >3(40)</td><td align="center" valign="middle" >3(PZT5H)</td><td align="center" valign="middle" >1(2)</td><td align="center" valign="middle" >2(25)</td><td align="center" valign="middle" >57.358</td><td align="center" valign="middle" >0.0662</td><td align="center" valign="middle" >5.92</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >3(40)</td><td align="center" valign="middle" >4(PZT8)</td><td align="center" valign="middle" >2(4)</td><td align="center" valign="middle" >1(20)</td><td align="center" valign="middle" >79.207</td><td align="center" valign="middle" >0.0660</td><td align="center" valign="middle" >5.60</td></tr><tr><td align="center" valign="middle" >13</td><td align="center" valign="middle" >4(50)</td><td align="center" valign="middle" >1(PZT4)</td><td align="center" valign="middle" >4(6)</td><td align="center" valign="middle" >2(25)</td><td align="center" valign="middle" >126.082</td><td align="center" valign="middle" >0.0650</td><td align="center" valign="middle" >4.00</td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" >4(50)</td><td align="center" valign="middle" >2(PZT5A)</td><td align="center" valign="middle" >3(5)</td><td align="center" valign="middle" >1(20)</td><td align="center" valign="middle" >120.888</td><td align="center" valign="middle" >0.0658</td><td align="center" valign="middle" >5.28</td></tr><tr><td align="center" valign="middle" >15</td><td align="center" valign="middle" >4(50)</td><td align="center" valign="middle" >3(PZT5H)</td><td align="center" valign="middle" >2(4)</td><td align="center" valign="middle" >4(35)</td><td align="center" valign="middle" >140.638</td><td align="center" valign="middle" >0.0656</td><td align="center" valign="middle" >4.96</td></tr><tr><td align="center" valign="middle" >16</td><td align="center" valign="middle" >4(50)</td><td align="center" valign="middle" >4(PZT8)</td><td align="center" valign="middle" >1(2)</td><td align="center" valign="middle" >3(30)</td><td align="center" valign="middle" >103.561</td><td align="center" valign="middle" >0.0662</td><td align="center" valign="middle" >5.92</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Analysis of orthogonal experiments for two single objectives</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Optimization objective</th><th align="center" valign="middle" >Analysis type</th><th align="center" valign="middle" >Planted depth/mm</th><th align="center" valign="middle" >PZT materials</th><th align="center" valign="middle" >PZT thickness/mm</th><th align="center" valign="middle" >PZT diameter/mm</th></tr></thead><tr><td align="center" valign="middle"  rowspan="5"  >Open voltage U<sub>3</sub>/V</td><td align="center" valign="middle" >k<sub>1</sub></td><td align="center" valign="middle" >75.574</td><td align="center" valign="middle" >111.785</td><td align="center" valign="middle" >75.333</td><td align="center" valign="middle" >83.160</td></tr><tr><td align="center" valign="middle" >k<sub>2</sub></td><td align="center" valign="middle" >96.642</td><td align="center" valign="middle" >108.718</td><td align="center" valign="middle" >110.308</td><td align="center" valign="middle" >89.346</td></tr><tr><td align="center" valign="middle" >k<sub>3</sub></td><td align="center" valign="middle" >102.4512</td><td align="center" valign="middle" >84.893</td><td align="center" valign="middle" >100.883</td><td align="center" valign="middle" >110.052</td></tr><tr><td align="center" valign="middle" >k<sub>4</sub></td><td align="center" valign="middle" >122.7922</td><td align="center" valign="middle" >92.064</td><td align="center" valign="middle" >110.935</td><td align="center" valign="middle" >114.902</td></tr><tr><td align="center" valign="middle" >R</td><td align="center" valign="middle" >47.218</td><td align="center" valign="middle" >26.891</td><td align="center" valign="middle" >35.602</td><td align="center" valign="middle" >31.742</td></tr><tr><td align="center" valign="middle"  rowspan="5"  >Relative deformation ratio p</td><td align="center" valign="middle" >k<sub>1</sub></td><td align="center" valign="middle" >4.8</td><td align="center" valign="middle" >4.72</td><td align="center" valign="middle" >5.88</td><td align="center" valign="middle" >5.28</td></tr><tr><td align="center" valign="middle" >k<sub>2</sub></td><td align="center" valign="middle" >4.76</td><td align="center" valign="middle" >4.88</td><td align="center" valign="middle" >5.00</td><td align="center" valign="middle" >4.84</td></tr><tr><td align="center" valign="middle" >k<sub>3</sub></td><td align="center" valign="middle" >4.88</td><td align="center" valign="middle" >4.92</td><td align="center" valign="middle" >4.68</td><td align="center" valign="middle" >4.76</td></tr><tr><td align="center" valign="middle" >k<sub>4</sub></td><td align="center" valign="middle" >5.04</td><td align="center" valign="middle" >4.96</td><td align="center" valign="middle" >3.92</td><td align="center" valign="middle" >4.60</td></tr><tr><td align="center" valign="middle" >R</td><td align="center" valign="middle" >0.28</td><td align="center" valign="middle" >0.24</td><td align="center" valign="middle" >1.96</td><td align="center" valign="middle" >0.68</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Optimization result for two single objectives</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Optimization objective</th><th align="center" valign="middle" >Optimization conditions (in order of dominance degree)</th></tr></thead><tr><td align="center" valign="middle" >Open voltage U<sub>3</sub>/V</td><td align="center" valign="middle" >Planted depth 50 mm, PZT thickness 6 mm, PZT diameter 35 mm, PZT material PZT4</td></tr><tr><td align="center" valign="middle" >Relative deformation ratio p</td><td align="center" valign="middle" >PZT thickness 6 mm, PZT diameter 35 mm, planted depth 20 mm, PZT material PZT4</td></tr></tbody></table></table-wrap><fig-group id="fig2"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Comparison of experimental results and optimum conditions.</title></caption><fig id ="fig2_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1560199x29.png"/></fig></fig-group><p>pared to the several highest values of the orthogonal experimental results in <xref ref-type="table" rid="table1">Table 1</xref>, as shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>(a). The optimization objective of open voltage under the optimum conditions is higher than the results from the orthogonal experiment. While, in <xref ref-type="fig" rid="fig2">Figure 2</xref>(b), the relative deformation ratio p of 1.92% under the selected optimum conditions is listed in contrast to the smallest values of the orthogonal experimental results in <xref ref-type="table" rid="table1">Table 1</xref>. From this figure, the optimization objective of p under the optimum conditions is the smallest one compared to the ones in <xref ref-type="table" rid="table1">Table 1</xref>.</p></sec><sec id="s3_3"><title>3.3. Multi-Objectives Optimization</title><p>Considering the practical application of harvester-pavement piezoelectric system, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x30.png" xlink:type="simple"/></inline-formula>and p have to be considered simultaneously. It is necessary to optimize parameters for requirements of all objectives. Here, the weight grade method is used to solve the multi-objectives optimization problem. According to the method, each objective variable is converted to a normalized value in a range of 0 - 1. Then, the multi-objectives optimization problem is conversed to a single-objective problem by introducing a comprehensive indicator. The new evaluation indicator, Y, of the system is computed as a mean product according to the following equation:</p><disp-formula id="scirp.60874-formula70"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/1-1560199x31.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x32.png" xlink:type="simple"/></inline-formula> is the weight factor of the two objective variables and is often determined according to the requirement of the users. Provided that U<sub>3</sub> is the same as important as p for the application of the harvester-pavement piezoelectric system, both <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x33.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x34.png" xlink:type="simple"/></inline-formula> take value 1 for the objective variables.</p><p>The following discussion is on the multi-objectives optimization. <xref ref-type="table" rid="table4">Table 4</xref> shows the range analysis of the orthogonal experiments for the objective variable, Y. The analysis results show that the optimum condition for the maximum value of Y in order of dominance degree is PZT thickness 6 mm, PZT diameter 35 mm, planted depth 50 mm and material PZT4. This is different from our above single-objective analysis results. It can be seen that the harvester structural parameters including diameter and thickness are the two most important factors for the overall objective while the PZT material has the most minor influence on it.</p><p>In the optimum conditions, the open voltage <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x35.png" xlink:type="simple"/></inline-formula> and relative displacement ratio p are 151.826 V and 3.52% (Corresponding to the pavement displacement 0.0647 mm), respectively. Obviously, the predicted value of open voltage is higher than these obtained from the orthogonal experiments, while the predicted value of relative deformation ratio is smaller than those in <xref ref-type="table" rid="table1">Table 1</xref>. It indicates the validity of the applied optimization procedure.</p></sec></sec><sec id="s4"><title>4. Effect of Each Factor on Objectives</title><sec id="s4_1"><title>4.1. Effect of Planted Depth</title><p><xref ref-type="fig" rid="fig3">Figure 3</xref> shows the effect of planted depth on open voltage <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x36.png" xlink:type="simple"/></inline-formula> and relative deformation ratio p. As can be observed from the figure, the open voltage increases with the increase in the planted depth. The relative deforma-</p><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Effect of Planted depth</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1560199x37.png"/></fig><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Analysis for multi-objectives optimization</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Analysis type</th><th align="center" valign="middle" >Planted depth/mm</th><th align="center" valign="middle" >PZT materials</th><th align="center" valign="middle" >PZT thickness/mm</th><th align="center" valign="middle" >PZT diameter/mm</th></tr></thead><tr><td align="center" valign="middle" >k<sub>1</sub></td><td align="center" valign="middle" >74.789</td><td align="center" valign="middle" >101.852</td><td align="center" valign="middle" >56.855</td><td align="center" valign="middle" >72.288</td></tr><tr><td align="center" valign="middle" >k<sub>2</sub></td><td align="center" valign="middle" >90.427</td><td align="center" valign="middle" >97.040</td><td align="center" valign="middle" >96.197</td><td align="center" valign="middle" >56.209</td></tr><tr><td align="center" valign="middle" >k<sub>3</sub></td><td align="center" valign="middle" >92.584</td><td align="center" valign="middle" >79.442</td><td align="center" valign="middle" >94.759</td><td align="center" valign="middle" >99.962</td></tr><tr><td align="center" valign="middle" >k<sub>4</sub></td><td align="center" valign="middle" >104.416</td><td align="center" valign="middle" >83.882</td><td align="center" valign="middle" >114.406</td><td align="center" valign="middle" >106.043</td></tr><tr><td align="center" valign="middle" >R</td><td align="center" valign="middle" >29.627</td><td align="center" valign="middle" >22.410</td><td align="center" valign="middle" >57.552</td><td align="center" valign="middle" >49.834</td></tr></tbody></table></table-wrap><p>tion ratio also increases with the planted depth except of 30 mm. It shows that the output voltage benefits from the bigger of the planted depth.</p></sec><sec id="s4_2"><title>4.2. Effect of PZT Materials</title><p>The effect of different PZT materials on open voltage and relative deformation ratio is shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>. It can be seen that the material of PZT4 can be harvested the maximum voltage and caused the minimum pavement deformation. The material of PZT5A also possesses better performance for harvesting pavement vibration with consideration of the harvested electricity and pavement deformation.</p></sec><sec id="s4_3"><title>4.3. Effect of PZT Thickness</title><p><xref ref-type="fig" rid="fig5">Figure 5</xref> shows the effect of PZT thicknesses on open voltage and relative deformation ratio. The relative deformation ratio decreases with increasing thickness of PZT disks as shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>. However, the open voltage shows a irregular trend.</p></sec><sec id="s4_4"><title>4.4. Effect of PZT diameter</title><p><xref ref-type="fig" rid="fig6">Figure 6</xref> shows the effect of PZT diameters on open voltage and relative deformation ratio. From the figure, open voltage increases and relative deformation ratio decreases while PZT diameter increases. But both vary in nonlinearity.</p></sec></sec><sec id="s5"><title>5. Finite Element Analysis Model</title><p>In order to save manpower, material and financial resources, the finite element method was used in the paper to</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Effect of PZT materials</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1560199x38.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Effect of PZT thickness</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1560199x39.png"/></fig><p>perform the electromechanical analysis. <xref ref-type="fig" rid="fig7">Figure 7</xref> shows the finite element model via ANSYS. Here, only one tire load is considered because the size of the harvester is small comparing to the contact area. Provided that the diameter of tire load area is 213 mm and the contact stress is 0.7 MPa. The diameter of asphalt pavement is 1000 mm. The height of pavement is set to be 150 mm according to the standard structure of pavement.</p><p>The mesh of Plane 42, which is suitable for elastic material properties, is used for asphalt pavement and end cap. Mean while, Plane 13 is adopted to mesh the piezoelectric material. The material properties can be found in [<xref ref-type="bibr" rid="scirp.60874-ref10">10</xref>] . The simulation results are filled in <xref ref-type="table" rid="table1">Table 1</xref> based on the parameter combination of the orthogonal experimental design.</p></sec><sec id="s6"><title>6. Conclusion</title><p>To optimize the harvester parameters and the embedded position in piezoelectric harvesting from pavement vi-</p><p>bration, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/1-1560199x40.png" xlink:type="simple"/></inline-formula>Taguchi’s orthogonal experiment schedule was designed and carried out considering the</p><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Effect of PZT diameter</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1560199x41.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Finite element model</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/1-1560199x42.png"/></fig><p>planted depth, PZT material, PZT thickness and diameter as optimization parameters. The finite element method was adopted to do experiment. It is demonstrated that the open voltage varies in wide range and the relative deformation ratio has less change under the conditions in orthogonal experiment design. Range analysis results show that the planted depth is the most important factor for open voltage while it is not so important for relative deformation ratio. The PZT diameter and thickness are the key important factors for all the two optimization objectives. Different PZT materials have fewer influences on two objectives. By using the weight grade method, the optimum conditions obtained from the two-objectives optimization are PZT thickness 6 mm, PZT diameter 35 mm, planted depth 50 mm and PZT4 material, which are in order of their dominance degrees. The validity of the optimum conditions is confirmed and is of practical value for the application of piezoelectric energy harvesting technology in pavement.</p></sec><sec id="s7"><title>Acknowledgements</title><p>This research was supported by the National Natural Science Foundation of China (No. 51175359) and the 4th “333 Engineering” Research Funding Project of Jiangsu Province (BRA2014086).</p></sec><sec id="s8"><title>Cite this paper</title><p>ChunhuaSun,HongbingWang,GuangqingShang,JianhongDu, (2015) Parameters Optimization for Piezoelectric Harvesting Energy from Pavement Based on Taguchi’s Orthogonal Experiment Design. World Journal of Engineering and Technology,03,149-157. doi: 10.4236/wjet.2015.34016</p></sec></body><back><ref-list><title>References</title><ref id="scirp.60874-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Huang, R.-B., Niu, Y.-L., Zhao, H.-D. and Chang, H.-B. (2012) Technical Approach and Research Prospect of Piezoelectric Energy Harvest from Highway. China Journal of Highway and Transport, 25, 1-8. (In Chinese)</mixed-citation></ref><ref id="scirp.60874-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Priya, S. and Inman, D.J. (2009) Energy Harvesting Technologies. Springer, New York. 
http://dx.doi.org/10.1007/978-0-387-76464-1</mixed-citation></ref><ref id="scirp.60874-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Zhao, H.D., Yu, J. and Ling, J.M. (2010) Finite Element Analysis of Cymbal Piezoelectric Transducers for Harvesting Energy from Asphalt Pavement. Journal of the Ceramic Society of Japan, 118, 909-915. 
http://dx.doi.org/10.2109/jcersj2.118.909</mixed-citation></ref><ref id="scirp.60874-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Sun, C.-H., Zhang, Y.-K., Du, J.-H. and Shang, G.-Q. (2013) Electromechanical Analysis of Piezoelectric Harvesting Unit from Road Vibration with FEA. Advanced Materials Research, 726-731, 3144-3147.</mixed-citation></ref><ref id="scirp.60874-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Lee, J., Lee, M., OH, J., et al. (2010) Study on the Energy Conversion from the Dynamic Load of Vehicles on the Road Using Piezoelectric Material. Materials Science Forum, 658, 57-60. 
http://dx.doi.org/10.4028/www.scientific.net/MSF.658.57</mixed-citation></ref><ref id="scirp.60874-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Sun, C.-H. and Shang, G.-Q. (2014) Effect of Traffic Flow on Characteristics of Piezoelectric Harvesting Unit. Applied Mechanics and Materials, 672-674, 902-905. http://dx.doi.org/10.4028/www.scientific.net/AMM.672-674.902</mixed-citation></ref><ref id="scirp.60874-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Innowattech Ltd. (2008) Innowattech’s Solution—The Innowattech Piezo Electric Generator (IPEG&amp;trade;) [EB/OL].</mixed-citation></ref><ref id="scirp.60874-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Wu, S.-L., Zhang, Q. and Huang, Z.P. (2010) Optimization of Resonance Frequency of Spherical-Cymbal Transducer Based on Orthogonal Experiment. Chinese Journal of Scientific Instrument, 31, 683-691.</mixed-citation></ref><ref id="scirp.60874-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Vives, A.A. (2008) Piezoelectric Transducers and Applications. 2nd Edition, Springer, New York. 
http://dx.doi.org/10.1007/978-3-540-77508-9</mixed-citation></ref><ref id="scirp.60874-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Zhao, X.M. (2006) Experimental Design Methods. Science Press, Beijing. (In Chinese)</mixed-citation></ref></ref-list></back></article>