基于壓電分流技術(shù)的PT 對稱梁散射特性研究

王剛，羅彩明，張曉東

（湖南大學(xué)汽車車身先進(jìn)設(shè)計(jì)制造國家重點(diǎn)實(shí)驗(yàn)室，湖南長沙 410082）

宇稱-時間（Parity-Time，PT）對稱由Bender 和Boettcher 于1998 年提出［1］，他們發(fā)現(xiàn)在滿足PT 對稱的條件下，系統(tǒng)的哈密頓量即使是非厄米的也可以具有實(shí)數(shù)的本征值.系統(tǒng)具有PT 對稱性的條件是勢函數(shù)共軛對稱，即勢函數(shù)的實(shí)部偶對稱、虛部奇對稱.PT 對稱系統(tǒng)的一個重要性質(zhì)為PT 對稱性的自發(fā)性破缺，具體表現(xiàn)是：當(dāng)勢函數(shù)虛部在某一閾值以下時，哈密頓量的本征值全部為實(shí)數(shù)，系統(tǒng)處于PT對稱相；當(dāng)勢函數(shù)虛部越過該閾值時，哈密頓量的本征值開始有復(fù)數(shù)產(chǎn)生，此時對應(yīng)于PT 破缺相.這一閾值也被稱為相變點(diǎn)或奇異點(diǎn)（Exceptional Point）.PT 對稱將量子力學(xué)拓展到了非厄米范圍，此后它便成了量子力學(xué)中的一個重要研究方向.現(xiàn)階段非厄米PT 對稱的研究工作已經(jīng)拓展到光學(xué)和聲學(xué)領(lǐng)域中，表現(xiàn)出許多奇特的物理現(xiàn)象.在光學(xué)領(lǐng)域的研究中，調(diào)節(jié)折射率的分布，使折射率在傳播方向上滿足實(shí)部偶對稱、虛部奇對稱，即可構(gòu)造PT 對稱光學(xué)系統(tǒng)［2-6］；而在聲學(xué)領(lǐng)域的研究中，當(dāng)聲學(xué)系統(tǒng)滿足一定條件時，它也會有PT 對稱性，例如Fleury 等人［7］通過調(diào)節(jié)兩個麥克風(fēng)連接的阻抗使得麥克風(fēng)的等效質(zhì)量密度共軛對稱，從而利用PT 對稱性實(shí)現(xiàn)了一種聲學(xué)隱身傳感器.

PT 對稱性的關(guān)鍵點(diǎn)在于引入均衡的增益和損耗，而目前，固體介質(zhì)中針對彎曲波的PT 對稱性的研究相對較少，其中一個重要的原因是自然介質(zhì)難以實(shí)現(xiàn)能量損耗和能量增益之間的平衡.壓電材料由于能夠?qū)崿F(xiàn)機(jī)械能與電能之間的轉(zhuǎn)化，并且具有易于調(diào)控的特點(diǎn)，有望解決這一問題.Christensen 等人［4］利用壓電半導(dǎo)體中的聲電效應(yīng)構(gòu)造了聲子PT對稱系統(tǒng)，并分析了其中的單向無反射現(xiàn)象.Hou 等人［8］基于壓電分流單元，提出了一種可調(diào)PT 對稱系統(tǒng)，通過改變分流電路的阻抗便可以調(diào)節(jié)PT 對稱系統(tǒng)奇異點(diǎn)出現(xiàn)的頻率.上述研究工作所提出的PT 對稱系統(tǒng)均是針對縱波而言的，所獲得的研究成果為PT 對稱系統(tǒng)的應(yīng)用奠定了理論基礎(chǔ).而與縱波相比，彎曲波的理論模型更為復(fù)雜，因此針對彎曲波的PT對稱系統(tǒng)同樣值得研究.

本文首先基于歐拉梁的假設(shè)，推導(dǎo)了針對彎曲波的PT 對稱條件；然后基于該條件，利用壓電分流單元設(shè)計(jì)一種PT 對稱梁，并采用等效介質(zhì)法和有限元仿真，得到增益單元、損耗單元的等效質(zhì)量密度和等效彎曲剛度；最后通過傳遞矩陣法和有限元仿真證明了PT 對稱梁中的單向無反射現(xiàn)象，并通過計(jì)算散射矩陣的特征值和特征向量說明了單向無反射現(xiàn)象源于奇異點(diǎn)的存在.本文所提出的PT 對稱梁對于彎曲波有許多潛在的應(yīng)用，包括增強(qiáng)傳感、彎曲波放大和非對稱控制等.

1 PT 對稱梁的設(shè)計(jì)與分析

1.1 彎曲波PT 對稱條件

根據(jù)歐拉梁假設(shè)，彎曲波沿x軸傳播時滿足控制方程［9］：

式中：w(x)、Deff(x)、ρeff(x)和hb分別表示梁的橫向位移、等效彎曲剛度、等效質(zhì)量密度和厚度.若梁滿足PT對稱性，則根據(jù)式（1）可以得到：

公式（3）和公式（4）表明，若滿足PT 對稱性，則要求等效質(zhì)量密度ρeff(x)和等效彎曲剛度Deff(x)關(guān)于原點(diǎn)共軛對稱.

1.2 PT 對稱梁設(shè)計(jì)

壓電分流單元已經(jīng)被證明可以有效地改變基體材料的等效彎曲剛度而幾乎不改變其等效質(zhì)量密度［10］，因此，可以通過設(shè)計(jì)合理的壓電分流單元使其滿足公式（3）和（4）所提出的等效質(zhì)量密度和等效彎曲剛度共軛對稱條件.圖1（a）為所設(shè)計(jì)的PT對稱梁示意圖，基體梁上表面貼有兩塊壓電片，兩塊壓電片分別連接由正、負(fù)電阻與電感并聯(lián)組成的分流電路，從而與基體梁部分構(gòu)成損耗單元、增益單元.分流電路中，正、負(fù)電阻的作用分別是衰減和放大彎曲波.正電阻Rsh和電感Lsh并聯(lián)后的分流阻抗ZL、負(fù)電阻-Rsh與電感Lsh并聯(lián)后的分流阻抗ZG分別為：

式中：ω為角頻率.

PT 對稱梁中的負(fù)電阻可以通過圖1（b）所示的Non-Foster電路來實(shí)現(xiàn)［11-12］，大電感可以用圖1（c）所示的Antoniou’s 模擬電感電路實(shí)現(xiàn)［13-14］.此外，分流電路中的電感會與壓電片的固有電容產(chǎn)生諧振，從而使得增益單元、損耗單元的等效彎曲剛度Deff(x)與彎曲波頻率相關(guān)，該現(xiàn)象在第2節(jié)進(jìn)一步說明.

圖1 PT對稱梁示意圖、Non-Foster電路及Antoniou’s模擬電感電路Fig.1 Schematic illustrations of PT symmetric beam，the Non-Foster circuit and the Antoniou’s circuit

1.3 PT對稱梁的特性分析

本節(jié)基于等效介質(zhì)理論計(jì)算增益單元、損耗單元的等效參數(shù)，并通過傳遞矩陣法，針對圖1（a）所示PT 對稱梁，計(jì)算彎曲波分別從左端和右端兩種入射情況下的透射系數(shù)和反射系數(shù).同時，使用有限元軟件COMSOL Multiphysics 進(jìn)行有限元仿真，以驗(yàn)證理論計(jì)算結(jié)果的準(zhǔn)確性.

如圖2 所示，令基體梁的厚度及彈性模量分別為hb、Eb，壓電片的厚度和長度分別為hp、Lp.假設(shè)壓電片除了垂直于x軸的端面外，其余表面均為自由狀態(tài)，且壓電片僅在z軸存在極化，因此壓電片的本構(gòu)方程可簡化為［15-16］：

圖2 等效參數(shù)求解原理示意圖Fig.2 Configurations to calculate the effective parameters

式中：T1和S1分別為壓電片x軸方向的應(yīng)力和應(yīng)變；D3和E3分別為壓電片的電位移和電場強(qiáng)度d31分別對應(yīng)于壓電片柔度系數(shù)、恒應(yīng)力介電常數(shù)、壓電應(yīng)力常數(shù).

另外，圖2 所示單元內(nèi)任意一點(diǎn)x方向的位移u(x，z) 為縱波位移u0(x) 與彎曲波位移wp(x) 的組合：

由于壓電片的厚度與基體梁的厚度處于同一個數(shù)量級，因此電場強(qiáng)度E3在z方向不能被認(rèn)為是恒定的.本文采用文獻(xiàn)［17］提出的假設(shè)，認(rèn)為電場強(qiáng)度E3在z方向呈線性分布，即電場強(qiáng)度E3和電勢V(x，z)可以用如下公式表述：

電位移D3滿足控制方程［18］：

聯(lián)立公式（6）～公式（10），便可以求出系數(shù)a(x)、b(x)、c(x)和壓電片上表面的電壓Vup［17］：

以基體梁的中性面為基準(zhǔn)面，對應(yīng)力求積分，可以得到單元x-y截面上的彎矩M、法向力N和剪切力T：

其中系數(shù)I、J、K、F和G通過下式計(jì)算：

基于公式（15）～公式（17），圖2所示單元的運(yùn)動控制方程為：

令γn表示特征方程（21）的特征根，則縱波和彎曲波wp(x)的表達(dá)式分別為：

對于沒有貼壓電片的基體梁，縱波wb(x)和彎曲波解耦，運(yùn)動控制方程分別為［19］：

對于基體梁中未貼壓電片的部分及貼有壓電片的部分，本文分別定義向量Yb和Yp：

根據(jù)公式（22）～公式（28），公式（15）～公式（17）可以用以下矩陣形式表述：

其中：

根據(jù)連續(xù)性條件，結(jié)合公式（31）和公式（32），可以得到增益單元或損耗單元的傳遞矩陣Mi，其中i=G表示增益單元，i=L表示損耗單元：

式中：Pp為對角矩陣.

1.3.1 等效參數(shù)

本文僅對波長遠(yuǎn)大于單元長度Lp的低頻彎曲波進(jìn)行分析，采用等效介質(zhì)理論計(jì)算結(jié)構(gòu)的等效參數(shù)［20］.因?yàn)樵鲆鎲卧蛽p耗單元的等效參數(shù)求解原理相同，對應(yīng)的等效參數(shù)計(jì)算過程相同，所以本文僅以損耗單元為例介紹等效參數(shù)求解過程.

令x=0，根據(jù)公式（36），Yb(Lp)與Yb(0)之間滿足：

求解等效質(zhì)量密度ρeff時，需要獲取單元的質(zhì)量，這要求單元任意一點(diǎn)的橫向加速度均相同.如圖2（a）所示，損耗單元整體相對于虛線方框所示的原始位置作幅值為A1的橫向簡諧振動，并且限制其左、右邊界的縱向位移和旋轉(zhuǎn)角度為0.依據(jù)上述邊界條件，基于公式（38）可求解向量中Yb(0)、Yb(Lp)的未知元素，從而得到損耗單元的等效質(zhì)量密度ρeff：

如圖2（b）所示，求解等效彎曲剛度Deff時，在單元左、右兩側(cè)相對于虛線方框所示的原始位置分別施加大小為-α1、α1的轉(zhuǎn)角，令wp(0)=wp(Lp)=0.彎曲波和縱波相互耦合，若限制單元兩端縱向位移為零，則需要確定增益單元、損耗單元的中性面.但是該中性面對于本文中的增益單元、損耗單元來說很難確定，因此，本文只限制單元左邊界的縱向位移(0)=0，而令右邊界的縱向力為零.依據(jù)上述邊界條件便可以基于公式（38）求解向量中Yb(0)、Yb(Lp)的未知元素，從而得到等效彎曲剛度Deff：

為了驗(yàn)證等效介質(zhì)法計(jì)算結(jié)果的正確性，本文針對損耗單元和增益單元，使用COMSOL 分別建立了如圖2所示2種情況下的二維有限元模型.有限元模型中，認(rèn)為單元滿足平面引力假設(shè)，施加的邊界條件與等效介質(zhì)法中的邊界條件一致.

求解等效質(zhì)量密度ρeff時，分別在左、右邊界上，對單元所受的z方向力求線積分即可得到對應(yīng)邊界上的剪切力T(0)、T(Lp)；而求解等效彎曲剛度Deff時，在右邊界上對單元所受的y方向力矩求線積分可得到力矩M(Lp).得到剪切力和力矩之后，便可以通過公式（39）和公式（40）分別計(jì)算得到單元的等效質(zhì)量密度ρeff及等效彎曲剛度Deff.

等效參數(shù)法和有限元仿真中，損耗單元和增益單元的幾何參數(shù)及材料參數(shù)如表1 所示.另外，電感與壓電片的固有電容會產(chǎn)生諧振，因此取諧振頻率f0=500 Hz，分流電阻Rsh=1 000 Ω.等效參數(shù)法和有限元仿真的計(jì)算結(jié)果在2.1節(jié)進(jìn)行討論.

表1 損耗單元和增益單元的幾何參數(shù)及材料參數(shù)Tab.1 Geometrical and material parameters of loss and gain units

1.3.2 透射系數(shù)及反射系數(shù)

為了表述方便，本文用下標(biāo)l表示彎曲波從左端入射，用下標(biāo)r 表示彎曲波從右端入射，將彎曲波從左端入射時的透射系數(shù)和反射系數(shù)分別定義為左透射系數(shù)tl和左反射系數(shù)rl.同理，定義彎曲波從右端入射時的右透射系數(shù)tr和右反射系數(shù)rr.

在PT 對稱梁的左端或右端施加一個單位位移，即

式中：L表示兩塊壓電片之間的距離.求解公式（41）便可以分別得到左透射系數(shù)tl、左反射系數(shù)rl、右反射系數(shù)rr、右透射系數(shù)tr，進(jìn)而可以得到彎曲波的散射矩陣S［21］：

為驗(yàn)證傳遞矩陣法計(jì)算結(jié)果的正確性，本文使用COMSOL 軟件建立如圖3 所示的有限元仿真模型，基體梁上含有增益單元、損耗單元.分別在左端激勵點(diǎn)和右端激勵點(diǎn)施加力載荷，產(chǎn)生彎曲波.為了抑制邊界反射，基體梁的兩端各連接一個完美匹配層，其他邊界保持自由狀態(tài).完美匹配層的長度至少為研究波長的1∕2［22］，在本文中，波長的最大值約為0.384 4 m，則本文中完美匹配層的長度設(shè)為0.2 m.由于激勵點(diǎn)關(guān)于原點(diǎn)O對稱，所以觀測點(diǎn)相對原點(diǎn)也要對稱分布.為了避免受到激勵點(diǎn)的影響，觀測點(diǎn)與激勵點(diǎn)之間要保持一定距離.觀測點(diǎn)A和觀測點(diǎn)B的坐標(biāo)分別為（-0.35 m，0）和（0.35 m，0），左端激勵點(diǎn)和右端激勵點(diǎn)的坐標(biāo)分別為（-0.5 m，8.0×10-4m）和（0.5 m，8.0×10-4m）.圖4 為試驗(yàn)裝置圖，在圖4 中，梁的兩端附著有藍(lán)丁膠，用來減少反射.信號發(fā)生器產(chǎn)生的信號經(jīng)過功率放大器后作用在壓電激勵器上，產(chǎn)生彎曲激勵.激光測振儀用來測量指定點(diǎn)的位移響應(yīng)，它與信號調(diào)理器、數(shù)據(jù)采集卡、電腦共同組成信號采集處理系統(tǒng).試驗(yàn)中所需要的負(fù)電阻和電感可以分別通過圖1（b）和圖1（c）所示電路來實(shí)現(xiàn).

圖3 有限元仿真示意圖Fig.3 Configuration for finite elements simulation

圖4 試驗(yàn)裝置圖Fig.4 Diagram of experimental setup

傳遞矩陣法和有限元仿真中，壓電片長度Lp=45 mm，兩塊壓電片之間的距離L=300 mm，分流電阻Rsh=1 900 Ω，分流電路的諧振頻率f0=500 Hz，其他參數(shù)與表1 中所列參數(shù)相同.有限元仿真模型的長度為2 005 mm.

2 結(jié)果和討論

2.1 等效參數(shù)

基于上述等效介質(zhì)法和有限元仿真，計(jì)算得到增益單元和損耗單元的等效質(zhì)量密度ρeff和等效彎曲剛度Deff.為了表達(dá)方便，本文采用歸一化等效質(zhì)量密度ρeff∕ρb和歸一化等效彎曲剛度Deff∕Db來表征增益單元、損耗單元的等效參數(shù).等效介質(zhì)法和有限元仿真計(jì)算得到的歸一化等效參數(shù)分別如圖5 和圖6 所示.由圖5（a）可知，通過2 種方法得到的增益單元和損耗單元的歸一化等效質(zhì)量密度實(shí)部Re(ρeff∕ρb)在計(jì)算頻率范圍內(nèi)幾乎保持不變，并且與分流阻抗無關(guān).由圖5（b）可知，歸一化等效質(zhì)量密度虛部Im(ρeff∕ρb)數(shù)值很小，可以認(rèn)為，等效質(zhì)量密度ρeff在計(jì)算的頻率范圍內(nèi)為實(shí)數(shù).由圖6（a）和圖6（b）可以發(fā)現(xiàn)，增益單元和損耗單元的等效彎曲剛度Deff共軛.因此，本文提出的PT 對稱梁滿足公式（3）和公式（4）所提出的等效質(zhì)量密度ρeff和等效彎曲剛度Deff共軛對稱條件.由于電感和分流電阻Rsh的存在，增益單元和損耗單元的等效彎曲剛度Deff與彎曲波頻率相關(guān).等效介質(zhì)法和有限元仿真獲取的結(jié)果總體上相吻合，存在部分差異的原因主要在于等效介質(zhì)法存在的局限性，如壓電片z方向上電場線性分布的假設(shè)等.

圖5 歸一化等效質(zhì)量密度Fig.5 The normalized effective mass density

圖6 歸一化等效彎曲剛度Fig.6 The normalized effective bending stiffness

為了進(jìn)一步表征歸一化等效彎曲剛度Deff∕Db與分流電阻Rsh、彎曲波頻率的關(guān)系，本文采用等效介質(zhì)法計(jì)算不同分流電阻Rsh及不同頻率的彎曲波激勵下?lián)p耗單元的歸一化等效彎曲剛度Deff∕Db，計(jì)算結(jié)果分別如圖7 和圖8 所示.由圖7 可知，當(dāng)彎曲波的頻率遠(yuǎn)離530 Hz時，分流電阻Rsh的取值不會影響歸一化等效彎曲剛度實(shí)部Re(Deff∕Db) 和虛部Im(Deff∕Db)，并且虛部Im(Deff∕Db)接近0，表明此時損耗單元幾乎沒有消耗彎曲波能量；當(dāng)彎曲波頻率為530 Hz 附近時，歸一化等效彎曲剛度實(shí)部Re(Deff∕Db)出現(xiàn)極值，并且增大分流電阻Rsh會使得歸一化等效彎曲剛度虛部Im(Deff∕Db)數(shù)值減小，從而加強(qiáng)損耗作用.當(dāng)分流電路的諧振頻率f0=800 Hz 時，可以發(fā)現(xiàn)圖8（a）和圖8（b）所示結(jié)果分別與圖7（a）和圖7（b）所示結(jié)果相近，只是歸一化等效彎曲剛度Deff∕Db的梯度變大，并且隨著分流電路諧振頻率f0的改變，極值出現(xiàn)的頻率變?yōu)?50 Hz附近.

圖7 諧振頻率f0=500 Hz時損耗單元?dú)w一化等效彎曲剛度與分流電阻Rsh、彎曲波頻率的關(guān)系Fig.7 The normalized effective bending stiffness of loss element in relationship to shunting resistance and flexural wave frequency when resonant frequency f0=500 Hz

圖8 諧振頻率f0=800 Hz時損耗單元?dú)w一化等效彎曲剛度與分流電阻Rsh、彎曲波頻率的關(guān)系Fig.8 The normalized effective bending stiffness of loss element in relationship to shunting resistance and flexural wave frequency when resonant frequency f0=800 Hz

由圖5～圖8 可知，本文提出的PT 對稱梁滿足公式（3）和公式（4）所提出的等效質(zhì)量密度ρeff和等效彎曲剛度Deff共軛對稱條件.通過改變分流電阻Rsh和諧振頻率f0，增益單元、損耗單元的彎曲剛度Deff會產(chǎn)生較大變化.因此，本文所提出的PT 對稱梁具有易于調(diào)控的特點(diǎn).

2.2 透射系數(shù)、反射系數(shù)及散射矩陣

基于傳遞矩陣法和有限元仿真得到的透射系數(shù)和反射系數(shù)分別如圖9 和圖10 所示.圖9 表示彎曲波分別從左端和右端入射的情況下透射系數(shù)的幅值，由圖9 可知，2 種方法得到的透射系數(shù)幅值一致，這是因?yàn)镻T 對稱梁具有互易性，即左透射系數(shù)與右透射系數(shù)相等.在圖9 中的陰影部分內(nèi)，透射系數(shù)的幅值大于1，這是由于PT 對稱梁是非保守系統(tǒng)，與外界存在能量交換.圖10（a）和圖10（b）分別表示左反射系數(shù)幅值和右反射系數(shù)幅值.由圖10 可以發(fā)現(xiàn)，PT 對稱梁對彎曲波的非對稱散射特性，即當(dāng)彎曲波分別從左、右兩端入射時，2 種情況下的反射系數(shù)幅值不相等，2 種情況下的反射系數(shù)都存在零點(diǎn)，該零點(diǎn)對應(yīng)的彎曲波頻率即為單向無反射點(diǎn)，如圖10（a）中的520.5 Hz 及圖10（b）中的511 Hz.由圖9 可以發(fā)現(xiàn)，在單向無反射點(diǎn)511 Hz 和520.5 Hz，透射系數(shù)的幅值接近1.

圖9 透射系數(shù)幅值圖Fig.9 The amplitude of transmission coefficients

圖10 反射系數(shù)幅值圖Fig.10 The amplitude of reflection coefficients

根據(jù)PT對稱性，左反射系數(shù)rl、右反射系數(shù)rr、透射系數(shù)t之間滿足：

根據(jù)公式（44），可以得到廣義守恒定律［23］：

式中：T=|t|2表示透射率；Rl=|rl|2表示彎曲波左端入射時的反射率；Rr=|rr|2表示彎曲波右端入射時的反射率.

由公式（45）可知，當(dāng)T＜1時，T+=1，表明彎曲波從左、右兩端入射時，PT 對稱梁中不存在耗散或者放大.當(dāng)T＞1 時，T-=1，PT 對稱梁對彎曲波進(jìn)行放大.在單向無反射點(diǎn)，由于反射系數(shù)為0，根據(jù)公式（45）可推知T=1，表明此時PT 對稱梁對彎曲波完全透射，這與圖9所示結(jié)果相符.

為了更好地闡述單向無反射特性，本文從有限元仿真的結(jié)果中，提取PT 對稱梁在2 個單向無反射點(diǎn)511 Hz 和520.5 Hz 下的歸一化彎曲波位移場幅值，結(jié)果如圖11所示.由圖11可知，區(qū)域Ⅰ和區(qū)域Ⅱ能直觀地體現(xiàn)PT 對稱梁的散射特性，區(qū)域Ⅰ代表左端激勵點(diǎn)到損耗單元左邊界的范圍；區(qū)域Ⅱ代表增益單元右邊界與右端激勵點(diǎn)之間的范圍.當(dāng)511 Hz的彎曲波從左端入射時，區(qū)域Ⅰ中有明顯的干涉條紋，表明區(qū)域Ⅰ中入射彎曲波與反射彎曲波形成了駐波，導(dǎo)致區(qū)域Ⅰ中各點(diǎn)的彎曲位移幅值明顯不一致.當(dāng)511 Hz 的彎曲波從右端入射時，區(qū)域Ⅱ中各點(diǎn)的彎曲位移幅值近似相等，表明區(qū)域Ⅱ中只有向右傳播的彎曲波而幾乎不存在反射彎曲波，并且區(qū)域Ⅰ和區(qū)域Ⅱ的彎曲位移幅值相近，表明此時PT 對稱梁對右端入射的彎曲波完全透射.當(dāng)520.5 Hz 的彎曲波從左端入射時，區(qū)域Ⅰ中幾乎無反射彎曲波存在，并且彎曲波幾乎完全透射至區(qū)域Ⅱ；而從右端入射時，區(qū)域Ⅱ中反射彎曲波與入射彎曲波形成了駐波，導(dǎo)致區(qū)域Ⅱ中各點(diǎn)之間的彎曲位移幅值有較大差別.

圖11 兩個單向無反射點(diǎn)下的歸一化彎曲位移場幅值Fig.11 The normalized amplitude of flexural wave displacement fields at two unidirectional reflectionless points

通過散射矩陣S也可以表征PT對稱梁的散射性質(zhì).由公式（43）定義的散射矩陣，可以得到散射矩陣的特征值λ1、λ2和特征向量V1、V2，結(jié)果分別如圖12和圖13 所示.圖12 和圖13 中的實(shí)心圓點(diǎn)代表單向無反射點(diǎn)520.5 Hz，陰影部分表示透射率T＞1 的彎曲波頻率范圍.將圖12（a）與圖9所示結(jié)果對比可以發(fā)現(xiàn)：當(dāng)T≤1 時，兩個特征值的模均為1，散射矩陣的兩個特征向量滿足，此時對應(yīng)于PT對稱相；而當(dāng)T＞1 時，其中一個特征值的模大于1，另一個特征值的模小于1，散射矩陣的兩個特征向量滿足此時對應(yīng)于PT 破缺相［24］.在單向無反射點(diǎn)520.5 Hz，同時也是PT 對稱相與PT 破缺相之間的臨界點(diǎn)（相變點(diǎn)，也即奇異點(diǎn)），兩個特征值發(fā)生簡并.由圖13 可知，在單向無反射點(diǎn)520.5 Hz，兩個特征向量也會發(fā)生簡并，這說明PT 對稱梁的奇異點(diǎn)即為單向無反射點(diǎn).此外，圖10（a）所示的單向無反射點(diǎn)也可以通過特征向量來解釋.例如，在520.5 Hz 的奇異點(diǎn)，兩個特征向量的第二分量為零，表明從左端入射時反射率為0，這與圖10（a）中的結(jié)果一致.

圖12 特征值幅值及相角Fig.12 The amplitude and phase angle of eigenvalues

圖13 特征向量實(shí)部和虛部Fig.13 The real and imaginary parts of the eigenvectors

3 結(jié)論

1）利用壓電分流單元設(shè)計(jì)一種針對彎曲波的PT 對稱梁，其中壓電分流單元連接了正∕負(fù)電阻并聯(lián)電感組成的分流電路，構(gòu)成損耗單元、增益單元；采用等效介質(zhì)法計(jì)算損耗單元、增益單元的等效彎曲剛度和等效質(zhì)量密度，并通過有限元仿真對等效介質(zhì)法進(jìn)行驗(yàn)證.結(jié)果表明，本文所提出的設(shè)計(jì)滿足等效質(zhì)量密度和等效彎曲剛度共軛對稱.

2）采用傳遞矩陣法計(jì)算PT 對稱梁的透射系數(shù)和反射系數(shù)，并與有限元仿真的結(jié)果進(jìn)行比較，驗(yàn)證傳遞矩陣法的正確性.傳遞矩陣法和有限元仿真的結(jié)果表明，所提出的PT 梁針對彎曲波具有非對稱散射特性，并且存在單向無反射點(diǎn).

3）通過計(jì)算散射矩陣的特征值與特征向量，分析PT 對稱梁的奇異點(diǎn)與單向無反射點(diǎn)之間的關(guān)系，結(jié)果表明，PT對稱梁的奇異點(diǎn)即為單向無反射點(diǎn).