謝濤濤, 羅 祎, 肖大為

水下非正交聲學(xué)角反射器聲散射特性研究

謝濤濤, 羅 祎, 肖大為

(海軍工程大學(xué) 兵器工程學(xué)院, 湖北 武漢, 430033)

水下角反射器具有較強(qiáng)的遠(yuǎn)場聲散射特性, 由于存在生產(chǎn)裝配誤差, 標(biāo)準(zhǔn)正交金屬角反射器無法滿足各板面絕對垂直, 因此生產(chǎn)誤差是對角反射器聲散射特性產(chǎn)生影響的重要因素。文中研究非正交角反射器的聲散射特性, 利用夾逼定理思想改變標(biāo)準(zhǔn)角反射器平板之間的夾角為90°、85°和80°, 并進(jìn)行定性分析; 采用結(jié)構(gòu)建模軟件ANSYS依次對非正交二面角、三面角反射器進(jìn)行模型構(gòu)建及網(wǎng)格劃分; 再利用聲學(xué)仿真軟件SYSNOISE根據(jù)結(jié)構(gòu)有限元+流體邊界元的方法對其遠(yuǎn)場聲散射特性進(jìn)行仿真, 并與標(biāo)準(zhǔn)正交角反射器在相同的載荷條件與聲場條件下的聲散射特性進(jìn)行對比, 得到聲波不同入射角度下的目標(biāo)聲散射強(qiáng)度以及目標(biāo)強(qiáng)度隨入射聲波頻率的變化規(guī)律。仿真結(jié)果表明, 正交角反射器聲散射特性優(yōu)于非正交角反射器, 正交二面角反射器優(yōu)于三面角反射器, 二面角反射器生產(chǎn)誤差在2°以內(nèi)不會對聲散射特性造成實質(zhì)性影響。

水下角反射器; 生產(chǎn)誤差; 結(jié)構(gòu)有限元; 流體邊界元; 聲散射; 目標(biāo)強(qiáng)度

0 引言

角反射器由金屬板材或者覆蓋金屬的其他材料與非金屬組成, 具有結(jié)構(gòu)簡單、回波能力強(qiáng)、造價低等優(yōu)點, 最初用于雷達(dá)電磁波對抗領(lǐng)域[1-3]。水下角反射器作為一種無源聲反射裝置, 在聲波入射其表面時, 不僅會產(chǎn)生鏡面反射, 還會產(chǎn)生繞射、透射和遠(yuǎn)場散射, 因其較好的聲散射特性, 可作為聲學(xué)假目標(biāo)對主動聲吶進(jìn)行誘騙[4-6]。

金屬平板是構(gòu)成角反射器的基礎(chǔ), 相互垂直的金屬板構(gòu)成的角反射器反聲性能較好, 但水下角反射器在實際加工條件下并非相互垂直, 聲波入射到角反射器時經(jīng)過多次反射回到入射聲源處, 而角反射器平板間角度的變化會直接影響到角反射器的聲回溯特性, 使得聲波反射方向無序化, 從而減小了聲源處所能接收到的目標(biāo)強(qiáng)度。由于生產(chǎn)工藝的限制, 滿足各平板之間相互垂直的條件, 無疑會極大地增加水下角反射器的生產(chǎn)成本。因此, 研究水下非正交角反射器的聲散射特性具有重要意義。

畢亞峰等[7]結(jié)合Biot流體理論進(jìn)行了基于變換聲學(xué)的角反射器分析; 張小鳳等[8]結(jié)合“柱的方法”研究了水下充氣圓柱的聲散射特性; 魏克難等[9]通過耦合聲學(xué)邊界元法(boundary element method, BEM)分析了水中彈性球殼的聲散射特性; Zeidan[10]采用解析法分析了有限長加肋柱狀結(jié)構(gòu)的聲散射情況; 陳鑫等[11-12]通過流固耦合的方法對水下標(biāo)準(zhǔn)二面角反射器的聲散射特性進(jìn)行了研究, 結(jié)合聲束彈跳法驗證了SYSNOISE軟件仿真數(shù)據(jù)的準(zhǔn)確性, 并對產(chǎn)生的誤差進(jìn)行分析, 認(rèn)為加工角度誤差使得角反射器平面不垂直會對角反射器聲散射特性造成較大影響。

文中以正交角反射器為基礎(chǔ), 通過結(jié)構(gòu)有限元+流體邊界元的方法, 分別對不正交彈性二面角和三面角的聲散射特性進(jìn)行研究, 分析角度誤差對角反射器聲散射特性的影響, 得出不影響角反射器聲散射特性的加工誤差范圍。

1 有限元-邊界元基本理論

1.1 有限元法簡介

有限元法的基本概念是用較簡單的問題代替復(fù)雜問題后再求解。它將求解域看成是由許多稱為有限元的小的互連子域組成, 對每一單元假定1個合適的(較簡單的)近似解, 然后推導(dǎo)求解這個域總的滿足條件(如結(jié)構(gòu)的平衡條件), 從而得到問題的解。

有限元法實際上就是求解下面的方程

另一方面如果能知道模態(tài), 通過模態(tài)疊加, 就可以得到聲腔的反應(yīng)

這樣就轉(zhuǎn)化為求解系數(shù)向量的問題, 通過求解下列方程得到特征值向量

1.2 BEM簡介

BEM又稱邊界積分方程法, 是將待求解問題的微分方程變成邊界上的積分方程, 再引入邊界上的有限個單元將積分方程離散, 得到只含有邊界上節(jié)點未知量的方程組, 并進(jìn)行數(shù)值求解。由于BEM是將區(qū)域上的控制方程轉(zhuǎn)化為沿著區(qū)域邊界的積分方程, 因此只需要定義邊界上的單元, 結(jié)合邊界條件求解, 這樣就使處理問題的維數(shù)降低一維, 可以有效地將復(fù)雜的三維幾何模型簡化為二維圖形。與同樣復(fù)雜的完整三維有限元模型相比, 該模型更小、更易于創(chuàng)建和檢驗, 并在更短的時間內(nèi)得出結(jié)果。但BEM所建立方程組的系數(shù)矩陣稠密, 而且一般非對稱, 矩陣元素分量的計算量很大, 增加了計算時間。由于BEM所利用的微分算子基本解能自動滿足無限遠(yuǎn)處的條件, 因而BEM特別便于無限域以及半無限域問題。BEM又可分為直接BEM和間接BEM。

1) 直接BEM的系統(tǒng)方程為

2) 間接BEM的系統(tǒng)方程為

1.3 有限元耦合間接邊界元

在求解空氣中結(jié)構(gòu)的聲散射問題時, 由于空氣的特性阻抗較之于結(jié)構(gòu)阻抗十分微小, 基本可以忽略空氣對于結(jié)構(gòu)振動的影響。但對于水中結(jié)構(gòu), 由于水的聲特性阻抗是空氣的100多倍, 且固體的結(jié)構(gòu)阻抗和水的特性阻抗相差1個數(shù)量級, 會產(chǎn)生流體與結(jié)構(gòu)之間的耦合振動, 因此不能忽略水的阻抗作用。因此水中結(jié)構(gòu)的聲散射問題實質(zhì)上是對流體結(jié)構(gòu)振動的復(fù)合問題, 需進(jìn)行聲振耦合計算, 使用結(jié)構(gòu)有限元+流體邊界元的方法對問題進(jìn)行分析解決。

在物理坐標(biāo)系中, 直接邊界元流體模型和有限元結(jié)構(gòu)模型的耦合系統(tǒng)方程如下

2 模型構(gòu)建及分析流程

文中采用ANSYS和SYSNOISE軟件相結(jié)合的方式對不標(biāo)準(zhǔn)水下聲學(xué)角反射器的聲散射特性進(jìn)行了仿真分析。SYSNOISE善于處理復(fù)雜結(jié)構(gòu)的振動聲學(xué)問題, 但沒有有限元及邊界元模型網(wǎng)格生成的前處理功能。因此, 先采用ANSYS進(jìn)行模型構(gòu)建及網(wǎng)格劃分, 而后將網(wǎng)格模型以*.cdb的格式輸出。SYSNOISE則負(fù)責(zé)導(dǎo)入此網(wǎng)格文件, 繼而進(jìn)行結(jié)構(gòu)-聲耦合后處理, 具體流程見圖1。

圖1 仿真分析流程圖

2.1 ANSYS模型前處理

在ANSYS模型前處理中, 為了減少計算誤差, 在建模時需通過布爾運(yùn)算實現(xiàn)對復(fù)雜結(jié)構(gòu)的精確建模。水下角反射器屬于薄壁結(jié)構(gòu), 采用SHELL181類型單元對模型進(jìn)行網(wǎng)格劃分, 為了保證計算精度, 網(wǎng)格大小劃分小于入射聲波波長的1/6, 網(wǎng)格法線方向一致。二面角反射器及三面角反射器結(jié)構(gòu)網(wǎng)格模型如圖2和圖3所示。

圖2 二面角反射器網(wǎng)格劃分圖

圖3 三面角反射器網(wǎng)格劃分圖

彈性金屬板的聲散射特性與金屬板厚度和入射波頻率有關(guān), 利用MATLAB軟件對水下彈性板聲反射系數(shù)進(jìn)行研究, 考慮到金屬密度及實驗所用信號發(fā)生器最高可輸出15 kHz聲波, 因此分析0～20 kHz內(nèi)0～20 mm厚彈性金屬板的反射系數(shù), 如圖4和圖5所示。

圖4 彈性金屬板反射系數(shù)隨入射波角度及頻率變化

圖5 反射系數(shù)隨金屬板厚度變化曲線

由圖可知, 彈性金屬板反射系數(shù)隨入射波頻率增大而增大, 入射波頻率在12 kHz時, 反射系數(shù)趨近于1。入射波頻率為15 kHz時, 反射系數(shù)隨金屬板厚度增大而增大, 在板厚超過10 mm時, 反射器系數(shù)大于0.85, 在板厚為18 mm時, 反射系數(shù)趨近于1; 考慮到金屬的密度及角反射器體積, 取角反射器厚度為10 mm。

完成模型構(gòu)建與網(wǎng)格劃分后, 將網(wǎng)格數(shù)據(jù)輸出為*.cdb文件, 便于導(dǎo)入SYSNOISE軟件中進(jìn)行聲學(xué)后處理分析。

2.2 SYSNOISE聲學(xué)后處理

在SYSNOISE中進(jìn)行聲-結(jié)構(gòu)耦合處理時, 需要先將結(jié)構(gòu)模型的分析類型定義為FEM Structure, 導(dǎo)入模型網(wǎng)格后進(jìn)行檢查及材料屬性定義; 流體邊界元模型需定義為BEM Indirect, 同理在模型網(wǎng)格導(dǎo)入后仍需定義流體材料屬性及耦合面, 最后完成有限元結(jié)構(gòu)模型與流體邊界元模型的耦合。

為了使計算滿足遠(yuǎn)場條件, 定義聲源距離聲場中心100 m, 聲源入射波為平面波, 入射聲波幅值為1 Pa, 頻率設(shè)定為0～15 kHz, 方向為聲場中心方向。場點與聲源的位置相同, 滿足收發(fā)合置的條件。

定義求解頻率范圍和步長, 對模型網(wǎng)格和場點上的聲學(xué)量, 如聲壓、聲強(qiáng)和振速等參數(shù)進(jìn)行計算。場點計算完成后重置聲源, 再對場點進(jìn)行一次運(yùn)算, 此時得到的計算結(jié)果是散射聲場的聲學(xué)量, 不含入射聲場。根據(jù)需要提取指定場點在某頻率范圍內(nèi)的計算結(jié)果, 繪制該場點的頻率響應(yīng)函數(shù)曲線。

3 仿真算例分析

3.1 彈性非正交二面角聲散射特性分析

標(biāo)準(zhǔn)二面角反射器是由2塊矩形金屬板相互垂直構(gòu)成(如圖6所示), 當(dāng)聲波入射時, 由于方向回溯特性, 可產(chǎn)生較為理想的后向反射。

圖6 標(biāo)準(zhǔn)二面角反射器示意圖

平板互不垂直時, 構(gòu)成了非正交二面角反射器, 通過研究非正交角反射器的聲散射特性, 分析加工誤差對角反射器帶來的影響。

圖7 非正交二面角反射器示意圖

在入射聲波屬性一致的前提下, 分析圖8～圖10可以得出如下結(jié)論。

圖8 不同入射頻率下90°二面角反射器目標(biāo)強(qiáng)度隨入射角變化曲線

1) 二面角反射器對入射聲波具有較強(qiáng)的頻率特性, 目標(biāo)強(qiáng)度隨入射聲波頻率的增大而增大;

2) 二面角反射器目標(biāo)強(qiáng)度隨平板間誤差角變化明顯, 誤差越大, 目標(biāo)強(qiáng)度越低, 且波動性越強(qiáng);

圖9 不同入射頻率下85°二面角反射器目標(biāo)強(qiáng)度隨入射角變化曲線

圖10 不同入射頻率下80°二面角反射器目標(biāo)強(qiáng)度隨入射角變化曲線

3)目標(biāo)強(qiáng)度隨聲波入射角變化規(guī)律基本一致, 高頻聲波入射時, 始終在45°處取得最大目標(biāo)強(qiáng)度。

基于二面角反射器對入射波的頻率特性, 將15 kHz高頻下二面角反射器目標(biāo)強(qiáng)度進(jìn)行提取對比, 如圖11所示, 可知:

1) 正交角二面反射器目標(biāo)強(qiáng)度在相同入射條件下大于非正交二面角反射器, 90°正交角反射器目標(biāo)強(qiáng)度比85°二面角反射器大6 dB;

2) 正交角反射器散射寬度優(yōu)于非正交角反射器, 在入射角為30°～60°范圍內(nèi)具有較強(qiáng)的目標(biāo)強(qiáng)度值;

3) 90°二面角反射器隨入射波角度改變目標(biāo)強(qiáng)度變化較小, 穩(wěn)定性強(qiáng), 85°和80°非正交角反射器則波動性較大。

圖11 不同夾角下15 kHz二面角反射器目標(biāo)強(qiáng)度隨入射角變化曲線

綜上所述, 二面正交角反射器聲散射特性優(yōu)于非正交角反射器, 在相同入射條件下具有更強(qiáng)的目標(biāo)強(qiáng)度與散射寬度。

3.2 彈性非正交三面角反射器散射特性分析

標(biāo)準(zhǔn)三面角反射器是由3塊三角形板相互垂直構(gòu)成的(如圖12所示), 聲波在入射后可在一定象限內(nèi)形成較強(qiáng)的散射聲場, 因此工程上常用多個三面角反射器組合來制作具有較大散射寬度的角反射器。

圖12 標(biāo)準(zhǔn)三面正交角反射器

非正交三面角反射器(如圖13所示), 是改變軸方向2個反射面Ⅰ、Ⅱ之間的夾角, 使其兩面夾角依次也為90°、85°和80°。

圖13 非正交三面角反射器

三面角反射器結(jié)構(gòu)參數(shù)和分析參數(shù)與3.1節(jié)相同, 仿真結(jié)果如圖14～圖16所示。

圖14 不同入射頻率下90°三面角反射器目標(biāo)強(qiáng)度隨入射角變化曲線

圖15 不同入射頻率下85°三面角反射器目標(biāo)強(qiáng)度隨入射角變化曲線

圖16 不同入射頻率下80°三面角反射器目標(biāo)強(qiáng)度隨入射角變化曲線

在相同的入射條件及載荷條件下, 可知:

2) 非正交三面角反射器最大目標(biāo)強(qiáng)度值無固定點, 頻率曲線波動性較強(qiáng), 而正交角反射器則在45°處取得最大值。

同理, 將三面角反射器高頻入射波下頻率與誤差角數(shù)據(jù)進(jìn)行提煉可以得到如下結(jié)論, 如圖17所示。

圖17 不同夾角下15 kHz三面角反射器目標(biāo)強(qiáng)度隨入射角變化曲線

1) 正交角三面反射器目標(biāo)強(qiáng)度值在相同入射條件下大于非正交三面角反射器, 90°正交角反射器目標(biāo)強(qiáng)度值比85°二面角反射器大3 dB;

2) 正交三面角反射器散射寬度優(yōu)于非正交角反射器, 在入射角為25°～65°范圍內(nèi)具有較強(qiáng)的目標(biāo)強(qiáng)度;

3) 正交三面角反射器隨入射波角度改變目標(biāo)強(qiáng)度變化較小, 穩(wěn)定性強(qiáng), 80°非正交角反射器則波動性極大, 無明顯變化特征。

由上述分析可知, 在相同結(jié)構(gòu)和流體材料屬性及入射聲源條件下:

1) 正交角反射器的聲散射特性優(yōu)于非正交角反射器, 角反射器板間夾角增大, 目標(biāo)強(qiáng)度與散射寬度急劇減小;

2) 正交二面角反射器目標(biāo)強(qiáng)度大于正交三面角反射器, 但散射寬度略小。

圖18 86°～90°二面角反射器目標(biāo)強(qiáng)度隨入射角變化曲線

4 結(jié)論

針對水下彈性非正交角反射器的特點, 運(yùn)用結(jié)構(gòu)有限元+流體邊界元的方法, 結(jié)合仿真軟件ANSYS與SYSNOISE對非正交角反射器的聲散射特性進(jìn)行了仿真, 對非正交彈性二面角反射器與非正交彈性三面角反射器的散射聲場進(jìn)行仿真分析, 得到以下結(jié)論:

1) 正交角反射器聲散射特性強(qiáng)于非正交角反射器;

2) 水下角反射器目標(biāo)強(qiáng)度變化隨著入射波頻率的增加而增大, 且最大值出現(xiàn)在入射角度為45°附近;

3) 相同幾何尺寸下二面角反射器聲散射特性更強(qiáng), 目標(biāo)強(qiáng)度值大于三面角反射器, 散射寬度略小于三面角反射器;

4) 二面角反射器加工誤差控制在2°以內(nèi)時, 角反射器聲散射特性并不會發(fā)生較大變化, 且2°的誤差相比于絕對垂直, 更容易實現(xiàn)生產(chǎn)加工, 具有較為重要的實際意義。

受篇幅所限, 文中僅對彈性二面角與三面角反射器的非正交形態(tài)進(jìn)行了分析, 對于剛性情況下非正交夾角的角反射器, 異型面構(gòu)成角反射器仍需進(jìn)一步研究。

[1] 湯廣富, 李華, 甘榮兵, 等.海戰(zhàn)場環(huán)境下角反射器干擾分析[J].電子信息對抗技術(shù), 2015, 30(5): 39-45, 84.

Tang Guang-fu, Li Hua, Gan Rong-bing, et al.Analysis of Corner Reflector Under Naval Battle Field[J].Electronic Warfare Technology, 2015, 30(5): 39-45, 84.

[2] 張志遠(yuǎn), 趙原源.新型二十面體三角形角反射器的電磁散射特性分析[J].指揮控制與仿真, 2018, 274(4): 138- 142.

Zhang Zhi-yuan, Zhao Yuan-yuan.Analysis of Electromagnetic Scattering Characteristic for New TypeIcosahedrons Triangular Trihedral Corner[J].Command Control & Simulation, 2018, 274(4): 138-142.

[3] 趙虎辰.一種X波段全向雷達(dá)角反射器陣列設(shè)計[J].河北省科學(xué)院學(xué)報, 2019, 36(1): 30-33.

Zhao Hu-chen.Design of an X-Band Omnidirectional Radar Angular Reflector Array[J].Journal of the Hebei Academy of Sciences, 2019, 36(1): 30-33.

[4] 徐海珠, 袁延藝, 劉雄厚, 等.線陣聲誘餌對抗智能魚雷尺度識別效果分析[J].艦船科學(xué)技術(shù), 2017, 39(5): 135-138.

Xu Hai-zhu, Yuan Yan-yi, Liu Xiong-hou, et al.On Performance Analysis of Linear Array Decoy in Confronting Smart Torpedo[J].Ship Science and Technology, 2017, 39(5): 135-138.

[5] 陳建青, 趙俊杰.一種尺度目標(biāo)模擬器的模擬逼真度分析[J].魚雷技術(shù), 2014, 22(6): 442-446.

Chen Jian-qing, Zhao Jun-jie.Fidelity Analysis of a Scale Target Simulator[J].Torpedo Technology, 2014, 22(6): 442-446.

[6] 羅祎, 文無敵, 陳鑫．一種水下目標(biāo)尺度特征無源模擬優(yōu)化方法[J]．海軍工程大學(xué)學(xué)報, 2018, 30(1): 98-102．

Luo Yi, Wen Wu-di, Chen Xin.A Passive Simulating Optimization Method of Measure Characteristic of Underwa- ter Object[J].Journal of Naval University of Engineering, 2018, 30(1): 98-102．

[7] 畢亞峰, 賈晗, 楊軍.基于變換聲學(xué)的角反射器的設(shè)計[J].應(yīng)用聲學(xué), 2019, 38(1): 52-57.

Bi Ya-feng, Jia Han, Yang Jun.The Design of Acoustical Corner Reflector Based on Transformation Acoustics[J].Applied Acoustics, 2019, 38(1): 52-57.

[8] 張小鳳, 趙俊渭, 王尚斌, 等.水下物理目標(biāo)聲反射器研究[J].兵工學(xué)報, 2004, 25(5): 600-603.

Zhang Xiao-feng, Zhao Jun-wei, Wang Shang-bin, et al.A Study on Underwater Passive Sound Reflectors[J].Acta Armamentarii, 2004, 25(5): 600-603.

[9] 魏克難, 李威, 雷明, 等.基于耦合邊界元法的水下目標(biāo)低頻聲散射特性[J].艦船科學(xué)技術(shù), 2014, 36(10): 32- 36.

Wei Ke-nan, Li Wei, Lei Ming, et al.The Low-frequency Acoustic Scattering Characteristics Study on Underwater Targets by the Coupled Boundary Element Method[J].Ship Science and Technology, 2014, 36(10): 32-36.

[10] Zeidan N I.GNSS Receivers for Weak Signals[A].The GNSS Technology and Applications Series[M].1rd ed.London, UK: Artech House, 2006.

[11] 陳鑫, 羅祎.水下剛性角反射器散射特性[J].聲學(xué)技術(shù), 2019, 38(3): 278-283.

Chen Xin, Luo Yi.Acoustic Scattering Characteristics of Underwater Rigid Corner Reflector[J].Technical Acoustics, 2019, 38(3): 278-283.

[12] 陳鑫, 羅祎, 李愛華.水下彈性角反射器聲散射特性[J].兵工學(xué)報, 2018, 39(11): 2236-2242.

Chen Xin, Luo Yi, Li Ai-hua.Acoustic Scattering Characteristics of Underwater Elastic Corner Reflectors[J].Acta Armamentarii, 2018, 39(11): 2236-2242.

Analysis of Acoustic Scattering Characteristics of Underwater Non Orthogonal Acoustic Corner Reflectors

XIE Tao-tao, LUO Yi, XIAO Da-wei

(College of Weaponry Engineering, Naval University of Engineering, Wuhan 430033, China)

Underwater corner reflectors possess strong far-field acoustic scattering characteristics, and because of the assembly error, the standard orthogonal metal corner reflector cannot meet the absolute vertical of each plate.The production error is an important factor that affects the sound scattering characteristics of corner reflectors.The acoustic scattering characteristics of non-orthogonal corner reflectors were studied.The pinch theorem idea was used to change the angle between standard corner reflectors (90°, 85°, and 80°), and a qualitative analysis was performed.The structural modeling software ANSYS was used to model and mesh the non-orthogonal dihedral and trihedral corner reflectors.The acoustic simulation software SYSNOISE was used to analyze the acoustic scattering characteristics of non-orthogonal dihedral and trihedral corner reflectors.Furthermore, far-field acoustic scattering characteristics were simulated and analyzed by the structural finite element and fluid boundary element methods.The acoustic scattering characteristics were compared with standard orthogonal corner reflectors under similar loading and acoustic field conditions.The acoustic scattering intensity of the target under different incident angles of sound waves and variation law of the target intensity with the frequency of the incident sound waves were obtained.The results showed that the acoustic scattering characteristics of orthogonal corner and orthogonal dihedral corner reflectors were better than those of non-orthogonal corner and trihedral corner reflectors, respectively.The production error of the dihedral corner reflector was less than 2° and will not have a substantial impact on the acoustic scattering characteristics.

underwater corner reflector; production error; structural finite element method; fluid boundary element method; acoustic scattering; target strength

謝濤濤, 羅祎, 肖大為.水下非正交聲學(xué)角反射器聲散射特性研究[J].水下無人系統(tǒng)學(xué)報, 2022, 30(1): 94-101.

TJ630; TB566

A

2096-3920(2022)01-0094-08

10.11993/j.issn.2096-3920.2022.01.012

2020-08-14;

2021-03-29.

謝濤濤(1995-),男, 碩士, 主要研究方向為兵器工程.

(責(zé)任編輯: 楊力軍)