賈文超，陳美霞，楊 丹

(華中科技大學(xué) 船舶與海洋工程學(xué)院，湖北 武漢 430074)

二維彈性平板繞流鎖定發(fā)生機(jī)理研究

賈文超，陳美霞，楊 丹

(華中科技大學(xué) 船舶與海洋工程學(xué)院，湖北 武漢 430074)

本文采用雙向流固耦合來(lái)研究一端固定二維平板渦激振動(dòng)特性，通過(guò)動(dòng)網(wǎng)格技術(shù)來(lái)實(shí)現(xiàn)流場(chǎng)和結(jié)構(gòu)之間力和位移的相互傳遞。對(duì)不同雷諾數(shù)下平板繞流渦分離進(jìn)行計(jì)算，成功捕捉到了鎖定現(xiàn)象，并對(duì)其發(fā)生機(jī)理進(jìn)行研究。計(jì)算結(jié)果表明，當(dāng)渦脫落頻率接近結(jié)構(gòu)固有頻率、結(jié)構(gòu)表面壓力分布與結(jié)構(gòu)模態(tài)一致、尾渦強(qiáng)度達(dá)到一定水平，即可產(chǎn)生鎖定現(xiàn)象。

鎖定；雙向流固耦合；渦激共振

0 引 言

艦船螺旋槳發(fā)出的“唱音”是螺旋槳噪聲中一種非常特別的噪聲源，它的產(chǎn)生與流固之間強(qiáng)耦合作用密切相關(guān)，這種強(qiáng)耦合作用即為“鎖定”現(xiàn)象。由于鎖定現(xiàn)象涉及到流場(chǎng)和結(jié)構(gòu) 2 個(gè)方面，其發(fā)生機(jī)理非常復(fù)雜。Philippe Ausoni等[1]以三維等直翼型為研究對(duì)象，進(jìn)行不同來(lái)流速度下水動(dòng)力性能和振動(dòng)特性試驗(yàn)，試驗(yàn)成功捕捉到鎖定現(xiàn)象，當(dāng)渦脫落頻率接近結(jié)構(gòu)扭轉(zhuǎn)固有頻率時(shí)，結(jié)構(gòu)振級(jí)達(dá)到最大，且在鎖定雷諾數(shù)區(qū)間內(nèi)，渦脫落頻率均與結(jié)構(gòu)固有頻率保持一致。A.Zobeiri等[2]對(duì)不同隨邊形狀三維等直翼在不同來(lái)流速度下的渦激振動(dòng)特性進(jìn)行了試驗(yàn)研究，試驗(yàn)結(jié)果表明，尖銳隨邊比更鈍的隨邊更能降低鎖定雷諾數(shù)區(qū)間內(nèi)的振級(jí)，但對(duì)鎖定頻率大小及區(qū)間寬度影響較小。Eliasson R等[3]給出了相同來(lái)流速度不同隨邊形狀翼型結(jié)構(gòu)的振級(jí)，對(duì)稱隨邊更有利于上下緣漩渦交替發(fā)放，從而激發(fā)較強(qiáng)的渦激振動(dòng)；非對(duì)稱隨邊則不利于漩渦的產(chǎn)生于發(fā)放，渦激振動(dòng)振級(jí)較低。

彈性結(jié)構(gòu)繞流鎖定研究的難點(diǎn)在于計(jì)算方法，需要采用雙向流固耦合方法進(jìn)行計(jì)算。目前針對(duì)雙向流固耦合算法的研究，主要方法有 2 種：一種是采用彈簧振子系統(tǒng)來(lái)模擬，結(jié)構(gòu)被視為剛體；另一種是流場(chǎng)和結(jié)構(gòu)同步求解，每一時(shí)間步交換耦合面上的力和位移，考慮了結(jié)構(gòu)彈性對(duì)流場(chǎng)的反饋效應(yīng)。

B.S. Carmo等[4]對(duì)不同來(lái)流速度下彈簧-圓柱系統(tǒng)振動(dòng)特性進(jìn)行了數(shù)值模擬，將得到的頻率鎖定區(qū)間劃分為初始區(qū)域和穩(wěn)定區(qū)域，在鎖定區(qū)間內(nèi)結(jié)構(gòu)響應(yīng)并未

保持恒定，在初始區(qū)域內(nèi)結(jié)構(gòu)振級(jí)急劇上升達(dá)到峰值；而在穩(wěn)定區(qū)域，振級(jí)逐漸下降。

隨著大型計(jì)算機(jī)的發(fā)展以及數(shù)值計(jì)算軟件的開發(fā)，雙向流固耦合計(jì)算成為了可能。Gluck等[5]采用CSD 和 CFD 數(shù)值算法軟件包來(lái)同時(shí)進(jìn)行結(jié)構(gòu)和流場(chǎng) 2個(gè)域的計(jì)算，通過(guò) MPCCI 插值算法進(jìn)行結(jié)構(gòu)和流場(chǎng)兩個(gè)獨(dú)立耦合面節(jié)點(diǎn)之間位移和力的相互傳遞，從而將結(jié)構(gòu)和流場(chǎng) 2 個(gè)域耦合在一起。Herwig Peters等[6]采用Workbench 中的多物理場(chǎng)耦合模塊對(duì)一端固支懸臂平板渦激振動(dòng)進(jìn)行了計(jì)算，研究了結(jié)構(gòu)楊氏模量對(duì)懸臂板位移響應(yīng)的影響規(guī)律。結(jié)果表明，當(dāng)流體介質(zhì)從空氣變?yōu)樗?，結(jié)構(gòu)共振楊氏模量增加；結(jié)構(gòu)在動(dòng)水中共振楊氏模量相比于靜水中楊氏模量并無(wú)顯著差異。

文獻(xiàn)[6]通過(guò)改變楊氏模量來(lái)研究鎖定現(xiàn)象，工程應(yīng)用性不強(qiáng)。本文在其算法的基礎(chǔ)上，針對(duì)艦船螺旋槳經(jīng)常出現(xiàn)的“唱音”現(xiàn)象，對(duì)一端固支柔性平板繞流進(jìn)行了變來(lái)流速度的二維數(shù)值模擬，研究了其振動(dòng)特性和水動(dòng)力特性，對(duì)平板及翼型這類狹長(zhǎng)結(jié)構(gòu)繞流問(wèn)題中常見的“鎖定”現(xiàn)象發(fā)生機(jī)理進(jìn)行研究，為螺旋槳“唱音”問(wèn)題的解決提供參考。

1 雙向流固耦合計(jì)算方法

1.1 初始流場(chǎng)控制方程

本文流場(chǎng)介質(zhì)為不可壓黏性流體，且不考慮熱交換，在當(dāng)前雷諾數(shù)下，選用基于時(shí)間平均的湍流模型。具體方法是將湍流運(yùn)動(dòng)看做是時(shí)間平均流動(dòng)和湍流脈動(dòng) 2 種運(yùn)動(dòng)的線性疊加，因此，可將速度和壓力進(jìn)行如下處理：式中：u，v，w 為平均速度；u’，v’，w’ 為脈動(dòng)速度；p 為平均壓力；p’ 為脈動(dòng)壓力。

將其代入瞬態(tài)的 N-S 方程，可得時(shí)均化的動(dòng)量方程[7]：

顯然上式多了一個(gè)應(yīng)力項(xiàng)，需要引入湍流模型以使其封閉。本文采用 SST k-ω 湍流模型，本模型是一個(gè)結(jié)合了 k-ε 的分區(qū)域湍流模型，在靠近壁面的邊界層內(nèi)使用 k-ω 模型，能更有效模擬較大的分離流動(dòng)，而在遠(yuǎn)離邊界層的自由流場(chǎng)中則使用 k-ε 模型，克服了 kω 對(duì)自由來(lái)流條件不太敏感的缺點(diǎn)，本模型是模擬平板繞流湍流運(yùn)動(dòng)的理想模型。

1.2 懸臂板振動(dòng)控制方程

懸臂板受迫振動(dòng)有限元控制方程為：

式中：M 為質(zhì)量矩陣；C 為阻尼矩陣；K 為剛度矩陣；w 為節(jié)點(diǎn)位移；F 為平板表面分布激勵(lì)力，即第1.1 小節(jié)計(jì)算得到平板表面脈動(dòng)壓力。

1.3 雙向流固耦合數(shù)據(jù)交換

雙向流固耦合的關(guān)鍵在于流場(chǎng)與平板之間的力和位移的相互傳遞，本文通過(guò) Ansys Workbench 平臺(tái)中的System Coupling 模塊來(lái)實(shí)現(xiàn)這種傳遞。首先求解 1.1 小節(jié)中初始流場(chǎng)（平板視為剛體），得到剛性平板表面脈動(dòng)壓力然后通過(guò)求解 1.2 小節(jié)中彈性平板控制方程，可得彈性平板表面位移響應(yīng)。將作為流場(chǎng)內(nèi)邊界的速度邊界條件，采用動(dòng)網(wǎng)格技術(shù)，重新求解 N-S 方程，得到更新之后壁面壓力，再將其重新加載到結(jié)構(gòu)表面，求解結(jié)構(gòu)動(dòng)響應(yīng)，如此反復(fù)迭代，直至壓力和位移滿足收斂條件。圖 1 為懸臂平板雙向流固耦合求解流程圖。

圖1 懸臂板雙向流固耦合求解流程圖Fig.1 Flow char of two-way fluid-structure interaction algorithm

2 數(shù)值計(jì)算模型

2.1 流場(chǎng) CFD 模型及結(jié)構(gòu)有限元模型

本文結(jié)構(gòu)模型為二維彈性懸臂平板，其基本尺寸為：平板沿來(lái)流方向長(zhǎng) L = 0.41 m，平板厚度 t = 0.02 m，密度 ρ = 7 850 kg/m3，楊氏模量 E = 32 921 Pa，泊松比 v = 0.3。平板有限元模型采用 solid186 單元，如圖 2 所示，沿來(lái)流方向劃分 216 份，沿厚度方向劃分 16 份，沿展向劃分 1 份，以模擬展向無(wú)限長(zhǎng)平板，結(jié)構(gòu)單元總數(shù)為 3 586。

流體計(jì)算域如圖 3 和圖 4 所示，流動(dòng)方向?yàn)?x 軸

正向，y 軸方向?yàn)榇怪庇趤?lái)流方向，左端流場(chǎng)入口邊界到平板左端的距離為 75 倍板厚，出口邊界到平板右端距離為 140 倍板厚，計(jì)算域上下邊界到平板上下壁面的距離為 40 倍板厚。網(wǎng)格劃分采用分塊劃分方法，結(jié)構(gòu)近壁處采用加密的結(jié)構(gòu)化網(wǎng)格，以更準(zhǔn)確地模擬近壁面復(fù)雜的流動(dòng)情況，遠(yuǎn)場(chǎng)則采用稀疏的非結(jié)構(gòu)化網(wǎng)格，確保獲得較快的計(jì)算速度。

圖2 懸臂板有限元模型Fig.2 Finite element model of plate

圖3 全局流場(chǎng)網(wǎng)格Fig.3 Fluid domain

圖4 近壁面網(wǎng)格Fig.4 Boundary layer

2.2 邊界條件

平板表面及流域上下為無(wú)滑移固壁：u = 0，v = 0。

2.3 數(shù)值求解方法

本文采用有限體積法來(lái)離散二維不可壓縮非定常粘性流體控制方程（2），動(dòng)量采用二階迎風(fēng)格式，時(shí)間項(xiàng)采用一階隱式積分法，速度和壓力耦合項(xiàng)采用Simple 算法來(lái)求解。

3 彈性平板振動(dòng)特性分析

由于雙向流固耦合計(jì)算十分緩慢，為了捕捉到平板鎖定現(xiàn)象，在進(jìn)行雙向流固耦合計(jì)算之前，需要確定計(jì)算雷諾數(shù)區(qū)間。靜水中懸臂平板第 2、3 階固有頻率分別為 0.184 5 Hz 和 0.567 5 Hz，通過(guò)對(duì)剛性平板繞流進(jìn)行試算，雷諾數(shù)為 500 和 1 280 時(shí)，渦脫落頻率分別為 0.185 7 Hz 和 0.561 8 Hz，與平板第 2、3 階固有頻率比較接近，因此，在 500 和 1 280 這 2 個(gè)雷諾數(shù)附近進(jìn)行變速度雙向流固耦合計(jì)算。

圖 5 給出了雙向流固耦合算法下不同雷諾數(shù)懸臂板最大位移響應(yīng)計(jì)算結(jié)果。圖 5（a）為結(jié)構(gòu)第 2、3 階共振峰，圖 5（b）為結(jié)構(gòu)第二階共振峰局部放大圖。從圖 5（a）可以看出，當(dāng)雷諾數(shù)在 500 和 1 280 附近，平板發(fā)生極大的位移響應(yīng)，這 2 個(gè)峰值工況對(duì)應(yīng)的彈性平板尾部渦發(fā)放頻率分別為 0.186 4 Hz 和 0.557 3 Hz，與表 2 中得到的靜水中平板第 2 階和第 3 階固有頻率基本吻合，在這 2 個(gè)雷諾數(shù)下，懸臂板被漩渦激發(fā)而產(chǎn)生共振。觀察圖 5（b），當(dāng)雷諾數(shù)接近鎖定區(qū)間時(shí)，懸臂板自由端位移響應(yīng)緩慢增大到“位移平臺(tái)”區(qū)，而當(dāng)雷諾數(shù)超過(guò)鎖定區(qū)間上限 532 時(shí)，結(jié)構(gòu)位移響應(yīng)迅速下降。

圖5 不同雷諾數(shù)下 2 種算法得到自由端總響應(yīng)的對(duì)比Fig.5 Bending vibration amplitude of the free edge with different Reynolds number

圖 6（a）給出了剛性和彈性平板渦脫落頻率隨雷諾數(shù)的變化關(guān)系。圖 6（b）給出了第 2 階共振峰處的頻率鎖定區(qū)間。雙向耦合算法在第 2 階渦激共振區(qū)域捕捉到鎖定現(xiàn)象，鎖定雷諾數(shù)區(qū)間為 520～532。在該區(qū)間內(nèi)，渦脫落頻率大體保持恒定，而剛性平板則不

能得到這樣的鎖定區(qū)間。

圖 7 為文獻(xiàn)[14]中給出的某翼型繞流表面一點(diǎn)速度響應(yīng)隨雷諾數(shù)的變化曲線。本文得到的響應(yīng)曲線與文獻(xiàn)[14]翼型實(shí)驗(yàn)得到的曲線非常相似，規(guī)律基本一致。

圖6 剛性及彈性平板的渦脫落頻率隨雷諾數(shù)的變化規(guī)律的對(duì)比Fig.6 Vortex shedding frequency of the rigid and elastic plate at different Reynolds number

圖7 某固定翼型繞流位移響應(yīng)[8]Fig.7 Vortex-induced vibration amplitude of a hydrofoil

4 鎖定形成條件研究

4.1 頻率條件

渦脫落頻率與結(jié)構(gòu)固有頻率相一致為鎖定發(fā)生的前提條件。表 1 給出了平板固有頻率及圖 5（a）中 2個(gè)峰值工況下渦脫落頻率的對(duì)比。由表 1 可知，當(dāng)渦脫落頻率接近結(jié)構(gòu)某一階固有頻率時(shí)，漩渦將激發(fā)結(jié)構(gòu)產(chǎn)生共振，此時(shí)流-固之間產(chǎn)生較強(qiáng)的耦合作用。

表1 平板固有頻率及對(duì)應(yīng)峰值工況下渦脫落頻率的對(duì)比Tab.1 The second and third order natural frequency of the plate under water and votex shedding frequency of the elastic plate

圖8 靜水中平板第 2，3 階模態(tài)云圖Fig.8 The second and third order modes of the cantilever plate under water

4.2 脈動(dòng)壓力分布條件

平板和翼型這類狹長(zhǎng)結(jié)構(gòu)更有利于渦激共振的產(chǎn)生，本質(zhì)上是由于這種線型使得脈動(dòng)壓力的分布特性與結(jié)構(gòu)固有模態(tài)發(fā)生高度吻合，從而激發(fā)結(jié)構(gòu)與流場(chǎng)之間的強(qiáng)耦合作用，最終產(chǎn)生鎖定。圖 9 所示為 2 個(gè)共振峰工況對(duì)應(yīng)平板上下表面脈動(dòng)壓力幅值分布，從導(dǎo)邊到隨邊，其幅值開始是緩慢增大，接近隨邊時(shí)則急劇增大，激勵(lì)力主要集中在平板自由端，而其他區(qū)域則非常小，這是因?yàn)闇u激振蕩脈動(dòng)壓力完全由尾渦脫落而產(chǎn)生，漩渦主要在尾部形成，因此激勵(lì)力也主要集中中尾端。圖 10 給出了平板法向脈動(dòng)壓力合力沿流向的時(shí)域分布，由圖可知平板升力方向合力（激勵(lì)力）幅值自導(dǎo)邊至隨邊相位基本相同，而幅值則逐漸

增大，與圖 8（a）給出的結(jié)構(gòu)第 2 階固有模態(tài)分布更加匹配。

圖9 平板上下表面壓力幅值分布Fig.9 The pressure distribution along the upper surface and the lower surface

圖10 平板法向脈動(dòng)壓力合力時(shí)域分布Fig.10 Normal pressure in time domain

從渦激振動(dòng)幅值來(lái)看，圖 5（a）中第 2 階共振峰值明顯小于第 1 階的，而第 2 階共振峰工況對(duì)應(yīng)平板表面脈動(dòng)壓力明顯大于第 1 階，顯然這是由于脈動(dòng)壓力分布特性與 1 階共振峰處的固有模態(tài)更加吻合。在第 2 階共振峰處，流-固之間耦合更加強(qiáng)，結(jié)構(gòu)振級(jí)更高；而第 3 階共振峰處，激勵(lì)力分布與模態(tài)分布不再匹配，流固之間為弱耦合作用，結(jié)構(gòu)振級(jí)較低。

4.3 隨邊形狀及尾渦強(qiáng)度條件

渦激振動(dòng)強(qiáng)度與脈動(dòng)壓力的幅值、相位與分布都密切相關(guān)，脈動(dòng)壓力則是由結(jié)構(gòu)尾部上下緣漩渦周期性交替發(fā)放產(chǎn)生，而不同的尾部形狀對(duì)漩渦的產(chǎn)生與脫落會(huì)產(chǎn)生重大影響。圖 11 所示為一剛性翼，弦長(zhǎng) 1 m，攻角為 0，來(lái)流速度為 0.5 m/s，分別對(duì)隨邊厚度為0.02 m，0.03 m 和 0.04 m 三個(gè)剛性翼型進(jìn)行繞流數(shù)值模擬。為了方便對(duì)比，各模型尾流場(chǎng)中渦量測(cè)點(diǎn)位置使用無(wú)量綱位置，3 種模型渦量監(jiān)測(cè)測(cè)點(diǎn)分別取在距離隨邊 1 ～ 15 倍隨邊厚度處。

圖11 翼型形狀及尾流場(chǎng)測(cè)點(diǎn)分布Fig.11 Geometrical model and location of the measurement point in the wake

圖12 不同隨邊厚度模型對(duì)應(yīng)渦強(qiáng)度的對(duì)比Fig.12 Vorticity magnitude of the foil with different thickness

圖 12 給出了不同隨邊厚度翼型尾流場(chǎng)中渦脫落與耗散過(guò)程中，漩渦中心強(qiáng)度的對(duì)比。由圖可知，相同來(lái)流速度下，較大的隨邊厚度翼型在近場(chǎng)將會(huì)形成強(qiáng)度較大的渦。鈍翼尾渦強(qiáng)度明顯強(qiáng)于尖翼，而隨邊處脈動(dòng)壓力幅值也更大，更容易激起結(jié)構(gòu)極強(qiáng)的振動(dòng)，如圖 13 所示。圖 14 所示為不同隨邊厚度模型渦量

場(chǎng)，尾部越“尖銳”，越不易產(chǎn)生較強(qiáng)的漩渦發(fā)放；尾部越“鈍”，則越有利于漩渦的產(chǎn)生與發(fā)放。

圖13 隨邊表面一點(diǎn)壓力時(shí)間歷程的對(duì)比Fig.13 The history of the pressure on the trailing edges with different thickness

圖14 不同隨邊厚度模型對(duì)應(yīng)渦脫落渦量云圖Fig.14 Vortex arrangement along the wake of the trailing edges with different thickness

5 結(jié) 語(yǔ)

本文采用計(jì)算流體力學(xué)方法和有限元方法來(lái)實(shí)現(xiàn)雙向流固耦合，并以懸臂平板為研究對(duì)象，計(jì)算在不同雷諾數(shù)下，彈性平板的振動(dòng)響應(yīng)特性，成功捕捉到了頻率鎖定現(xiàn)象，并給出了鎖定的形成條件，為螺旋槳唱音問(wèn)題提供了指導(dǎo)方向，最終得到以下結(jié)論：

1）當(dāng)懸臂平板尾部渦發(fā)放頻率與靜水中平板的固有頻率相一致的時(shí)候，尾部漩渦的周期性發(fā)放將會(huì)誘導(dǎo)平板產(chǎn)生極大的振動(dòng)響應(yīng)。

2）渦與結(jié)構(gòu)發(fā)生“鎖定”現(xiàn)象，表現(xiàn)為結(jié)構(gòu)位移響應(yīng)峰值區(qū)域加寬，這是由于在共振頻率附近區(qū)域，結(jié)構(gòu)表面脈動(dòng)壓力幅值大小完全按結(jié)構(gòu)模態(tài)分布，使得結(jié)構(gòu)共振區(qū)域加寬。

3）結(jié)構(gòu)和流場(chǎng)同時(shí)滿足頻率、壓力分布和隨邊形狀幾個(gè)條件時(shí)，流-固之間將產(chǎn)生強(qiáng)耦合作用，結(jié)構(gòu)振幅達(dá)到最大，流場(chǎng)渦脫落頻率出現(xiàn)鎖定現(xiàn)象，在此鎖定區(qū)內(nèi)，渦脫落頻率保持不變，結(jié)構(gòu)振幅也保持在一個(gè)很高的水平。而在其他固有頻率處，雖能發(fā)生渦激共振，但振級(jí)大大低于鎖定時(shí)的振級(jí)，此時(shí)流-固之間的耦合相對(duì)較弱。

[1]AUSONI P, FARHAT M, AIT BOUZIAD Y, et al. Kármán vortex shedding in the wake of a 2D hydrofoil: Measurement and numerical simulation[C]//IAHR int. Meeting of WG on Cavitation and Dynamic Problems in Hydraulic Machinery and Systems. Barcelone, Spain, 2006 .

[2]ZOBEIRI A, AUSONI P, AVELLAN F, et al. How oblique trailing edge of a hydrofoil reduces the vortex-induced vibration[J]. Journal of Fluids and Structures, 2012, 32: 78–89.

[3]ELIASSON R, LARSSON L, ORYCH M. Principles of yacht design[M]. A&C Black, 2014.

[4]CARMO B S, SHERWIN S J, BEARMAN P W, et al. Flow-induced vibration of a circular cylinder subjected to wake interference at low Reynolds number[J]. Journal of Fluids and Structures, 2011, 27(4): 503–522.

[5]GLüCK M, BREUER M, DURST F, et al. Computation of fluid–structure interaction on lightweight structures[J]. Journal of Wind Engineering and Industrial Aerodynamics, 2001, 89(14/15): 1351–1368.

[6]PETERS H, CHEN L, KESSISSOGLOU N. The effect of flow on the natural frequencies of a flexible plate[C]//INTER-NOISE and NOISE-CON Congress and Conference Proceedings. Melbourne Australia: Institute of Noise Control Engineering, 2014: 610–616.

[7]ANDERSON J D, WENDT J F. Computational fluid dynamics[M]. New York: McGraw-Hill, 1995.

[8]AUSONI P. Turbulent vortex shedding from a blunt trailing edge hydrofoil[D]. Lausanne: école Polytechnique Fédérale de Lausanne, 2009.

Mechanism research of lock-in on two-dimensional flow around elastic plate

JIA Wen-chao, CHEN Mei-xia, YANG Dan
(School of Naval Architecture and Ocean Engineering, Huazhong University of Science and Technology, Wuhan 430074, China)

Vortex-induced vibration of a cantilever plate is investigated numerically based on two-way fluid-structure interaction, which is achieved through the exchange of force and displacement data between the fluid and the structure field. In this work, vortex-induced vibration of a cantilever plate at different Reynolds number is simulated, simultaneously successfully captures the lock-in phenomena, and the mechanism was studied. The results suggest that when the vortex shedding frequency is close to the natural frequency of the structure, the surface pressure distribution is consistent with the structure modal, vortex intensity reaches a certain level, the lock-in phenomenon will occur.

lock-in；two-way fluid-structure interaction；vorticity induced vibration

U661.1

A

1672 – 7619(2016)11 – 0028 – 06

10.3404/j.issn.1672 – 7619.2016.11.005

2016 – 03 – 01；

2016 – 07 – 29

國(guó)家自然科學(xué)基金資助項(xiàng)目（51179071）

賈文超(1986 – )，男，碩士研究生，主要研究方向?yàn)榱骷ふ駝?dòng)。