楊金花，張鵬君

（長(zhǎng)沙理工大學(xué)土木與建筑學(xué)院，長(zhǎng)沙 410114）

電－熱－力載下BNNTs增強(qiáng)壓電板的非線性動(dòng)力響應(yīng)

楊金花，張鵬君

（長(zhǎng)沙理工大學(xué)土木與建筑學(xué)院，長(zhǎng)沙 410114）

研究了電－熱－力載荷作用下硼氮納米管（BNNTs）增強(qiáng)壓電板的非線性動(dòng)力響應(yīng)。運(yùn)用板的非線性應(yīng)變幾何關(guān)系及考慮溫度效應(yīng)的壓電理論，建立了BNNTs增強(qiáng)壓電板的本構(gòu)關(guān)系。通過(guò)變分推導(dǎo)出了結(jié)構(gòu)的非線性動(dòng)力控制方程。在空間上采用差分法離散，在時(shí)間上采用Newton-Newmark方法離散，整個(gè)問(wèn)題采用迭代法進(jìn)行求解。算例中，詳細(xì)討論了電壓、溫度、載荷、體積比以及外激勵(lì)頻率等因素對(duì)BNNTs增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響。

非線性動(dòng)力響應(yīng)；壓電；板；BNNT；電－熱－力載

硼氮納米管（BNNT）具有類似于碳納米管（CNT）的結(jié)構(gòu)，有很好的力學(xué)性能，并且與CNT相比，BNNT具有更好的抗氧化能力及強(qiáng)的壓電特性。另外，BNNT還具有穩(wěn)定的半導(dǎo)體性能［1－2］。硼氮納米管的這些性質(zhì)使其似乎可以作為各種納米級(jí)的電子和光子設(shè)備的替代材料。因此，BNNT越來(lái)越廣泛地用作復(fù)合結(jié)構(gòu)的加固組分。隨著科技的發(fā)展，一種具有壓電性質(zhì)的新型智能納米復(fù)合材料，即用壓電材料作基體而B(niǎo)NNT作增強(qiáng)材料，吸引了科學(xué)界的興趣?，F(xiàn)有文獻(xiàn)表明，對(duì)BNNTs增強(qiáng)壓電結(jié)構(gòu)的研究還很少，且大多討論的是線性靜力問(wèn)題。因此，有必要對(duì)這種結(jié)構(gòu)的非線性動(dòng)力行為進(jìn)行更廣泛的研究。

目前大多數(shù)研究都局限于討論無(wú)BNNT增強(qiáng)的壓電結(jié)構(gòu)的動(dòng)力行為。隨著微分求積（DQ）逐層建模技術(shù)的開(kāi)發(fā)和實(shí)現(xiàn)，Zhang等［7］研究了多層壓電復(fù)合材料板的自由振動(dòng)?；诟唠A剪切板彈性理論和壓電理論，Mao等［8］分析了壓電功能梯度板的非線性動(dòng)態(tài)響應(yīng)和主動(dòng)振動(dòng)控制。基于高階剪切變形板理論，F(xiàn)akhari等［9］使用有限元分析具有表面粘貼壓電層的功能梯度板在熱、電和力載條件下的非線性固有頻率和頻率響應(yīng)。為了分析簡(jiǎn)支壓電夾層板的自由振動(dòng)，Benjeddou等［10］提出了一種二維封閉式的解法。Behjat等［11］使用不同的機(jī)電載荷下的有限元方法，進(jìn)行了功能梯度壓電板的非線性靜態(tài)分析和自由振動(dòng)分析。Larbi等［12］給出了壓電復(fù)合結(jié)構(gòu)振聲問(wèn)題的理論公式及有限元解。Xia等［1］研究了熱環(huán)境下表面粘貼壓電纖維增強(qiáng)復(fù)合材料板的非線性振動(dòng)及動(dòng)態(tài)響應(yīng)。Dash等［14］分析了內(nèi)埋或表面粘結(jié)壓電層的復(fù)合材料層合板的非線性自由振動(dòng)特性。Wang等［15］使用有限元模型，研究了壓電復(fù)合材料板的負(fù)速度反饋控制的動(dòng)力穩(wěn)定。Liew等［16］使用一種有效的無(wú)網(wǎng)格方法，討論了帶壓電傳感器或驅(qū)動(dòng)器的復(fù)合材料層合板的主動(dòng)控制。Zhang等［17］提出了由1－3個(gè)纖維增強(qiáng)的壓電復(fù)合層合板及正交復(fù)合層合板的結(jié)構(gòu)振動(dòng)控制分析方程。Yang等［18］研究了復(fù)合壓電殼的動(dòng)力不穩(wěn)定性。Saviz等［19］給出了簡(jiǎn)支壓電層合圓柱殼層的動(dòng)力分析。Ying等［20］分析了壓電厚殼在邊界隨機(jī)激發(fā)平面應(yīng)變條件下的隨機(jī)響應(yīng)。D'Ottavio等［21］研究了內(nèi)埋壓電材料的多層殼的自由振動(dòng)問(wèn)題。以上研究中均沒(méi)有考慮BNNT增強(qiáng)壓電結(jié)構(gòu)。

近年來(lái)出現(xiàn)了一些對(duì)于BNNTs增強(qiáng)壓電結(jié)構(gòu)靜力屈曲問(wèn)題的研究。Salehi-Khojin等［22］研究了電－熱－力載下PVDF彈性介質(zhì)中BNNT的屈曲問(wèn)題，Mosallaie Barzoki等［23］研究了具有彈性核心的BNNTs增強(qiáng)的PVDF圓柱殼的扭轉(zhuǎn)屈曲。在另一項(xiàng)研究中，Mosallaie Barzoki等［24］研究了BNNT增強(qiáng)內(nèi)埋壓電殼的非線性屈曲。使用基于非局部圓柱壓電連續(xù)殼理論的虛位移原理，Ghorbanpour Arani等［25］討論了電－熱－力負(fù)荷條件下內(nèi)埋彈性介質(zhì)中雙壁硼氮納米管的軸向屈曲。到目前為止，BNNTs增強(qiáng)壓電板的非線性動(dòng)力研究尚未見(jiàn)報(bào)道。

出于這種考慮，本文我們研究電－熱－力載下BNNTs增強(qiáng)壓電板的非線性動(dòng)力響應(yīng)?？臻g上使用有限差分法，時(shí)間上采用Newton-Newmark方法，整個(gè)問(wèn)題通過(guò)迭代求解。討論了電壓，溫度，力載及體積比等對(duì)BNNTs增強(qiáng)壓電板的非線性動(dòng)力響應(yīng)的影響。

1 基本方程

圖1為一BNNTs增強(qiáng)壓電板，板的長(zhǎng)度為a，寬度為b，厚度為h，板的密度為ρ0。坐標(biāo)系o－xyz置于板的中面（z＝0），壓電板受有橫向動(dòng)載荷q，施加電壓為V及溫度升高ΔT。

圖1 BNNTs增強(qiáng)壓電板幾何示意圖Fig.1 Geometry of piezoelectric plate reinforced with BNNTs

下標(biāo)（，）表示對(duì)坐標(biāo)變量求偏導(dǎo)。

其非線性應(yīng)變－位移關(guān)系可表示為

式中，εx，εy，εxy為殼中面上相應(yīng)點(diǎn)的應(yīng)變，κx，κy，κxy分別為殼中面上相應(yīng)點(diǎn)的曲率和扭率改變量，且

對(duì)電－熱－力載荷作用下的壓電結(jié)構(gòu)而言，本構(gòu)關(guān)系可寫(xiě)為［23，29－30］

一般而言，鵬氮納米管（BNNTs）有兩種對(duì)稱結(jié)構(gòu)，即Z字形的和扶手椅形狀的。Z字形管具有縱向壓電響應(yīng)，而扶手椅形管存在與扭轉(zhuǎn)耦合的電偶極矩，因此本文采用前者。智能壓電板由兩種材料組成，基體材料為PVDF，增強(qiáng)材料為BNNTs，結(jié)構(gòu)整體的材料常數(shù)可用'XY（或YX）［27－28］矩形模型計(jì)算得到，其表達(dá)式為

上標(biāo)r，m分別指增強(qiáng)及基體材料，Vf為增強(qiáng)材料（BNNTs）在基體材料（PVDF）中所占的體積比。BNNTs增強(qiáng)壓電板的總能量Π可寫(xiě)為

式中，K代表動(dòng)能，U為應(yīng)變能，Γ為外載所做的功。應(yīng)變能為

考慮采用的是z字形結(jié)構(gòu)的BNNTs，則有

運(yùn)用瑞茲變分原理δΠ＝0，可得BNNTs增強(qiáng)壓電板的非線性動(dòng)力控制方程為

則電－熱－力載下BNNTs增強(qiáng)壓電板的無(wú)量綱動(dòng)力控制方程可寫(xiě)為

2 求解方法

為了尋求滿足初始條件（21）及邊界條件（20）的非線性運(yùn)動(dòng)微分方程組（19）的解，我們將無(wú)量綱位移函數(shù)W，V，U在空間和時(shí)域內(nèi)離散，以求其近似解。

在空間上采用差分法離散。對(duì)于線性項(xiàng)的處理，型如U，ξ，V，η，U，ξξ，W，ξ，W，ξξ，W，ξξξξ等各階偏導(dǎo)數(shù)項(xiàng)，以W，ξξ為例，采用如下的差分格式：

對(duì)照差分格式，可容易的求得控制方程中其余各線性項(xiàng)的各階偏導(dǎo)數(shù)的差分表達(dá)式。

在時(shí)間上將外載荷的作用時(shí)間τ等分為n等份，記每一等份為Δτ，整個(gè)問(wèn)題采用迭代法求解。在任一迭代步j(luò)中，控制方程和定解條件中的非線性項(xiàng)被線性化，且寫(xiě)為如下的形式［31］

式中，（y）jp是前二次迭代值的均值。對(duì)于初始迭代步，可以應(yīng)用二次外推法求解，即

對(duì)于不同的迭代步，系數(shù)Δ1，Δ2及Δ3可以分別取如下的數(shù)值：

對(duì)于三次非線性項(xiàng)的處理方法與二次非線性項(xiàng)的處理方法相同。

同時(shí)，采用Newmark方法，將控制方程中對(duì)時(shí)間變量的偏導(dǎo)數(shù)項(xiàng)如W，ττ在時(shí)間上離散。將外載荷的作用時(shí)間τ等分為n等份，記每一等份為Δτ＝τ／n，則在網(wǎng)格點(diǎn)j的加速度項(xiàng)W，ττ可表示為

經(jīng)過(guò)差分處理和線性化以后，非線性微分控制方程由最初的偏微分方程組轉(zhuǎn)化為差分格式的線性代數(shù)方程組。聯(lián)立求解此方程組，對(duì)每一時(shí)間步采用迭代法求解，迭代過(guò)程持續(xù)到相鄰兩次所得計(jì)算結(jié)果的誤差小于某一給定精度（本文為10－12），在第j步取得收斂解后，該迭代步終止，再進(jìn)行第j+1步的計(jì)算。

3 數(shù)值計(jì)算與討論

下面算例中討論了電－熱－力載下BNNTs增強(qiáng)壓電板的非線性動(dòng)力響應(yīng)。壓電板的四邊簡(jiǎn)支，幾何尺寸為a／h＝10，b／h＝10，壓電板基體材料為PVDF，增強(qiáng)材料為BNNT，材料常數(shù)如表1所示。橫向動(dòng)載荷取為

式中，

為外激勵(lì)幅值，ω為外激勵(lì)頻率（圖6除外，各算例中ω＝1）。各圖中垂直坐標(biāo)W

0

為板的無(wú)量綱中心撓度，水平坐標(biāo)τ為無(wú)量綱時(shí)間。

表1 材料常數(shù)Tab.1 Mechanical，electrical and thermal properties of PVDF and BNNT

電－熱－力載荷下BNNT所占體積比Vf對(duì)板中點(diǎn)非線性動(dòng)力響應(yīng)的影響如圖2（a）所示。本算例中，電壓V＝25 v，溫度ΔT＝50 K，外激勵(lì)幅值＝0.1。從圖中可看出隨著B(niǎo)NNT在基體中所占體積比Vf的增加，板中點(diǎn)的動(dòng)力響應(yīng)幅值減小。這是因?yàn)轶w積比Vf的增加會(huì)使得結(jié)構(gòu)的整體剛度變大，因而動(dòng)力響應(yīng)幅值會(huì)減小。另外，圖2（b）還給出了當(dāng)Vf＝0.6時(shí)本文的計(jì)算結(jié)果和用Ansys軟件模擬分析的結(jié)果，此圖中外激勵(lì)幅值、電壓、溫度與圖2（a）中的相同。從圖2（b）可見(jiàn)，本文的計(jì)算結(jié)果與用Ansys模擬分析所得的結(jié)果吻合較好，說(shuō)明本文的力學(xué)模型和分析方法是合理的。

圖3給出了正負(fù)電壓對(duì)不同體積比下BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響。圖中溫度ΔT＝0，外激勵(lì)幅值＝0.1。從圖中可看出在硼氮納米管增強(qiáng)壓電板兩端施加負(fù)電壓會(huì)使結(jié)構(gòu)的響應(yīng)幅值減小。這是由于施加負(fù)電壓會(huì)使硼氮納米管沿x方向極化，產(chǎn)生收縮，因而使得結(jié)構(gòu)更加緊湊堅(jiān)固，結(jié)構(gòu)的剛度增大，故結(jié)構(gòu)的非線性動(dòng)力響應(yīng)幅值減小。圖中同時(shí)也給出了正電壓對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響，從圖中可看出，正電壓會(huì)使BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)幅值增大。另外，比較圖3（a）、（b）可知，隨著硼氮納米管所占體積比Vf的增大，正電壓對(duì)壓電板非線性動(dòng)力響應(yīng)的影響要減弱，而負(fù)電壓的影響要增強(qiáng)，圖中難以看出此變化規(guī)律，為此將其列于表2所示。

圖2 體積比Vf對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響Fig.2 Effect of volume fraction Vfon the nonlinear dynamic response of BNNT-reinforced piezoelectric plate

圖3 電壓對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響Fig.3 Effect of applied voltage on the nonlinear dynamic response of BNNT-reinforced piezoelectric plate

表2 電壓對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響Tab.2 Effect of app lied voltage on the nonlinear dynam icresponse of BNNT-reinforced piezoelectric p late

溫度對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響如圖4所示。圖中硼氮納米管所占體積比Vf分別為0.6和0.8，電壓V＝0，外激勵(lì)幅值＝0.1。從圖中可看出，壓電板非線性動(dòng)力響應(yīng)幅值隨著溫度的升高而變大。另外，比較圖3（a）與圖3（b）可知，隨著硼氮納米管所占體積比Vf的增大，溫度對(duì)壓電板非線性動(dòng)力響應(yīng)的影響要減?。▍⒁?jiàn)表3）。

圖4 溫度對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響Fig.4 Effect of temperature on the nonlinear dynamic response of BNNT-reinforced piezoelectric plate

表3 溫度對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響Tab.3 Effect of tem perature on the nonlinear dynam ic response of BNNT-reinforced piezoelectric plate

圖5給出了橫向動(dòng)載荷幅值對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響。圖中體積比Vf＝0.6，溫度ΔT＝50 K，電壓V＝25 v。由圖5可知，當(dāng)壓電板承受的外載荷幅值越大，其在振動(dòng)過(guò)程中達(dá)到的幅值也越大。同樣，比較圖3（a）與圖3（b）可知，荷載對(duì)壓電板非線性動(dòng)力響應(yīng)的影響隨著Vf的增大而減?。▍⒁?jiàn)表4）。

表4 載荷對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響Tab.4 Effect ofmechanical load on the nonlinear dynam ic response of BNNT-reinforced piezoelectric plate

圖5 外載荷對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響Fig.5 Effect ofmechanical load on the nonlinear dynamic response of BNNT-reinforced piezoelectric plate

圖6給出了外激勵(lì)頻率對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響。圖中體積比Vf＝0.6，溫度ΔT＝50 K，電壓V＝25 v，外激勵(lì)幅值＝0.1。從圖中可看出，壓電板非線性動(dòng)力響應(yīng)幅值隨著外激勵(lì)頻率ω的增加而減小。比較圖3（a）與圖3（b）可知，隨著體積比Vf的增大，外激勵(lì)頻率對(duì)壓電板非線性動(dòng)力響應(yīng)的影響也減?。▍⒁?jiàn)表5）。

表5 外激勵(lì)頻率對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響Tab.5 Effect of external excitation frequency on the nonlinear dynam ic response of BNNT-reinforced piezoelectric p late

圖6 外激勵(lì)頻率對(duì)BNNT增強(qiáng)壓電板非線性動(dòng)力響應(yīng)的影響Fig.6 Effect of external excitation frequency on the nonlinear dynamic response of BNNT-reinforced piezoelectric plate

4 結(jié) 論

本文建立了電－熱－力載荷作用下BNNTs增強(qiáng)壓電板的非線性動(dòng)力響應(yīng)控制方程。研究結(jié)果表明電壓、溫度、荷載、體積比以及外激勵(lì)頻率等因素對(duì)BNNTs增強(qiáng)壓電板非線性動(dòng)力響應(yīng)有著重要的影響。結(jié)果如下：

（1）在BNNTs兩端施加正或負(fù)電壓可使壓電板非線性動(dòng)力響應(yīng)幅值增加或減??；

（2）BNNTs增強(qiáng)壓電板非線性動(dòng)力響應(yīng)幅值隨著外載荷及溫度的增加而增大；隨著外激勵(lì)頻率及硼氮納米管所占體積比的增加而減小。

（3）隨著硼氮納米管所占體積比Vf的增大，正電壓、溫度、荷載以及外激勵(lì)頻率等因素對(duì)壓電板非線性動(dòng)力響應(yīng)的影響會(huì)減小。

［1］Chen CW，Lee M H，Clark S J，Band gap modification of single-walled carbon nanotube and boron nitride nanotube under a transverse electric field［J］.Nanotechnology，2004，15：1837.

［2］Bansal N P，Hurst JB，Choi SR，Boron nitride nanotubereinforced glass composites［J］.J Am Ceram Soc，2006，89：388.

［3］Chen Y，Zou J，Campbell S J，et al，Boron Nitride nanotubes：Pronounced Resistance to Oxidation，Appl.Phys.Lett.84（2004）2430－2432.

［4］Suryavanshi A P，Yu M F，Wen J，etal.Elasticmodulus and resonance behavior of boron nitride nanotubes［J］.App. Phys.Lett.，2004，84：2527－2529.

［5］Terrones M，Romo-Herrera JM，Cruz-Silva E，etal.Mu?oz-Sandoval，J.J.Velázquez-Salazar，H.Terrones，Y.Bando，D.Golberg，Pure and Doped Boron Nitride Nanotubes，Mater.Today，2007，10：30－38.

［6］Zhi C，Bando Y，Tang C，et al.Characteristics of boron nitride nanotube-polyaniline composites［J］.Angew Chem Int Ed，2005，44：7929.

［7］Zhang Zan，F(xiàn)eng C，Liew K M.Three-dimensional vibration analysis of multilayered piezoelectric composite plate［J］.International Journal of Engineering Science，2006，44（7）：397－408.

［8］Mao Yi-qi，F(xiàn)u Yi-ming，Nonlinear dynamic response and active vibration control for piezoelectric functionally graded plate［J］.Journal of Sound and Vibration，2010，329（11）：2015－2028.

［9］Fakhari V，Ohadi A，Yousefian P.Nonlinear free and forced vibration behavior of functionally graded plate with piezoelectric layers in thermal environment［J］.Composite Structures，2011，93（9）：2310－2321.

［10］Benjeddou A，DeüJ F.A two-dimensional closed-form solution for the free-vibrations analysis of piezoelectric sandwich plates［J］.International Journal of Solids and Structures，2002，39（6）：1463－1486.

［11］Behjat B，Khoshravan M R.Geometrically nonlinear static and free vibration analysis of functionally graded piezoelectric plates［J］.Composite Structures，2012，94（3）：874－882.

［12］LarbiW，DeüJF，Ohayon R，F(xiàn)inite element formulation of smart piezoelectric composite plates coupled with acoustic fluid［J］Composite Structures，2012，94（2）：501－509.

［13］Xia Xian-kun，Shen Hui-shen.Nonlinear vibration and dynamic response of FGM plates with piezoelectric fiber reinforced composite actuators［J］.Composite Structures，2009，90（2）：254－262.

［14］Dash P，Singh B N.Nonlinear free vibration of piezoelectric laminated composite plate［J］.Finite Elements in Analysis and Design，2009，45（10）：686－694.

［15］Wang SY，Quek S T，Ang K K，Dynamic stability analysis of finite element modeling of piezoelectric composite plates ［J］.International Journal of Solids and Structures，2004，41 （3－4）：745－764.

［16］Liew K M，He X Q，Tan M J，et al.Dynamic analysis of laminated composite plates with piezoelectric sensor／actuator patches using the FSDTmesh-freemethod［J］.International Journal of Mechanical Sciences，2004，46（3）：411－431.

［17］Zhang H Y，Shen Y P.Vibration suppression of laminated plates with 1－3 piezoelectric fiber-reinforced composite layers equipped with interdigitated electrodes［J］.Composite Structures，2007，79（2）：220－228.

［18］Yang X M，Shen Y P.Dynamic instability of laminated piezoelectric shells［J］International Journal of Solids and Structures，2001，38：2291－2303.

［19］Saviz M R，Mohammadpourfard M.Dynamic analysis of a laminated cylindrical shell with piezoelectric layers under dynamic loads［J］.Finite Elements in Analysis and Design，2010，46：770－781.

［20］Ying ZG，Zhu X Q.Response analysis of piezoelectric shells in plane strain under random excitations［J］.Acta Mechanica Solida Sinica，2009，22：152－160.

［21］D'Ottavio M，Ballhause D，Kr?plin B，et al.Closed-form solutions for the free-vibration problem of multilayered piezoelectric shells［J］.Computers＆Structures，2006，84：1506－1518.

［22］Salehi-Khojin A，Jalili N.Buckling ofboronnitridenanotube reinforced piezoelectric polymeric composites subject to combined electro-thermo-mechanical loadings［J］.Compos Sci Technol，2008，68：1489－501.

［23］Mosallaie Barzoki A A，Ghorbanpour Arani A，Kolahchi R，et al.Electro-thermo-mechanical torsional buckling of a piezoelectric polymeric cylindrical shell reinforced by DWBNNTswith an elastic core［J］.Appl.Math.Model，2012，36：2977－2989.

［24］Mosallaie Barzoki A A，Ghorbanpour Arani A，Kolahchi R，et al.Nonlinear buckling response of embedded piezoelectric cylindrical shell reinforced with BNNT under electro-thermomechanical loadings using HDQM［J］.Composites Part B：Engineering，2013，44（1）：722－727.

［25］Ghorbanpour Arani A，Amir S，Shajari A R，et al.Electrothermo-mechanical buckling of DWBNNTs embedded in bundle of CNTs using nonlocal piezoelasticity cylindrical shell theory［J］.Composites Part B：Engineering，2012，43：195 －203.

［26］Sai N，Mele E J.Microscopic theory for nanotube piezoelectricity［J］.Phys Rev B，2003，68：1405.

［27］Tan P，Tong L，Steven G P.A 3D modeling technique for predicting the linear elastic properties of opened-packing woven fabric unit cell［J］.Composite Structures，1997，38：261－271.

［28］Tan P，Tong L.Micro-electromechanics models for piezoelectric-fiber-reinforced composite materials［J］.Compos.Sci.Technol，2001，61：759－769.

［29］Yang J，Xiang H J.Thermo-electro-mechanical characteristics of functionally graded piezoelectric actuators［J］.Smart Mater.Struct，2007，16：784－797.

［30］Yang J，Kitipornchai S，Liew K M.Nonlinear analysis of thermo-electro-mechanical behavior of shear deformable FGM plates with piezoelectric actuators［J］.Inter.J.Numeri.Methods Eng，2004，59：1605－1632.

［31］傅衣銘，結(jié)構(gòu)非線性動(dòng)力分析［M］.廣州：暨南大學(xué)出版社，1997.

Nonlinear dynam ic response of a piezoelectric p late reinforced w ith BNNTs under electro-thermo-mechanical loadings

YANG Jin-hua，ZHANG Peng-jun
（School of Civil Engineering and Architecture，Changsha University of Science＆Technology，Changsha 410114，China）

The nonlinear dynamic response of a piezoelectric plate reinforced with boron nitride nanotubes （BNNTs）under electro-thermo-mechanical loadings was studied.Using nonlinear strains based on the plate theory and the piezoelectric theory considering thermal effects，the constitutive relations of the piezoelectric plate reinforced with BNNTs were established.Then，the dynamic governing equations of the structure were derived with the variational principle.Using the differencemethod in space and Newton-Newmark method in time，thewhole problem was solved with the iteration method.In numerical examples，the effects of voltage，temperature，mechanical load，volume ratio and external excitation frequency on the nonlinear dynamic response of the piezoelectric plate reinforced with BNNTs were discussed.

nonlinear dynamic response；piezoelectric；plate；BNNT；electro-thermo-mechanical loadings

O341

A

10.13465／j.cnki.jvs.2015.21.026

國(guó)家自然科學(xué)基金（11102028）

2015－04－27 修改稿收到日期：2015－07－15

楊金花女，博士，副教授，1974年生