薛?鋼，劉延俊，薛祎凡，劉大輝

內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的動力學(xué)研究

薛?鋼1, 2，劉延俊1, 2，薛祎凡1，劉大輝3

(1. 山東大學(xué)海洋研究院，青島 266237；2. 山東大學(xué)機(jī)械工程學(xué)院，濟(jì)南 250061；3. 中集海洋工程研究院有限公司，煙臺 264670)

內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置是一種新型的波浪能發(fā)電系統(tǒng)，采用全封閉式結(jié)構(gòu)，具有可靠性高、使用壽命長等特點．本文提出了內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的海況自適應(yīng)能量輸出方案，能夠根據(jù)海況特點自適應(yīng)地調(diào)整偏心轉(zhuǎn)子的偏心距，從而提高能量輸出效率和穩(wěn)定性．針對波浪能發(fā)電裝置的整體水動力學(xué)特性和內(nèi)部偏心轉(zhuǎn)子的機(jī)構(gòu)動力學(xué)特性，建立了一種基于ABAQUS-XFlow聯(lián)合仿真的動力學(xué)耦合分析方法，研究了偏心轉(zhuǎn)子質(zhì)量、摩擦系數(shù)、波浪參數(shù)等因素對偏心轉(zhuǎn)子動力學(xué)特性的影響，并利用原理樣機(jī)實驗分析了隨機(jī)波浪激勵下偏心轉(zhuǎn)子的動力學(xué)特性．其中，偏心轉(zhuǎn)子的運(yùn)動性能間接反映了波浪能發(fā)電系統(tǒng)的功率輸出特性．結(jié)果表明：偏心轉(zhuǎn)子質(zhì)量是影響動力學(xué)特性的敏感因素，對其運(yùn)動特征的影響呈現(xiàn)出非單調(diào)性規(guī)律；而偏心轉(zhuǎn)子的運(yùn)動對偏心轉(zhuǎn)子與中心軸之間的摩擦系數(shù)的敏感度較低，且摩擦系數(shù)對偏心轉(zhuǎn)子運(yùn)動特征的影響呈現(xiàn)出單調(diào)性規(guī)律；波浪能發(fā)電裝置受到隨機(jī)波浪激勵后，在垂蕩和縱搖方向的運(yùn)動幅值較大，在其他方向上的運(yùn)動幅值相對較小，偏心轉(zhuǎn)子有較大概率保持沿單一方向轉(zhuǎn)動，且轉(zhuǎn)動較為平穩(wěn)，預(yù)期能夠穩(wěn)定捕獲波浪能．相關(guān)研究結(jié)論對內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的研發(fā)具有重要的借鑒意義．

偏心轉(zhuǎn)子；波浪能發(fā)電裝置；動力學(xué)；耦合分析

海洋中蘊(yùn)藏著大量清潔可再生能源，具有重要的開發(fā)價值[1]．波浪能是海洋可再生能源的一種，具有儲量大、分布范圍廣等特點，受到學(xué)者的廣泛關(guān)注[2].據(jù)統(tǒng)計，全球波浪能的儲量為1～10TW，可開發(fā)量巨大[3]．

為合理開發(fā)波浪能，學(xué)者研究了形式各樣的波浪能發(fā)電裝置(wave energy converter，WEC)，如振蕩水柱式、振蕩浮子式、擺式等[4-5]．盛松偉等[6]研制了一種鷹式大型漂浮式波浪能發(fā)電裝置，在設(shè)計波況下系統(tǒng)總效率可達(dá)20%以上．Xie等[7]提出了一種基于空間雙X形機(jī)構(gòu)的振蕩浮子式波浪能發(fā)電裝置，實現(xiàn)了應(yīng)用于跨海大橋的傳感器的自供電．Zhang等[8]提出了一種采用雙作用液壓缸的倒立擺式波浪能發(fā)電裝置，并實驗驗證了設(shè)計的合理性．

可靠性是波浪能發(fā)電裝置的重要考核指標(biāo). Hinck[9]最早提出了采用全封閉結(jié)構(gòu)的內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置，將能量轉(zhuǎn)換單元與海水完全隔離，避免運(yùn)動部件受到海水腐蝕，具有可靠性高、使用壽命長等特點，備受關(guān)注．基于此，Chen等[10]分析了浮體運(yùn)動阻尼和發(fā)電機(jī)阻尼對內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置發(fā)電性能的影響，Yurchenko等[11]提出了一種基于參數(shù)擺的新型波浪能發(fā)電裝置，能夠改善小幅值、低頻率海況下的發(fā)電性能．本文的研究對象即為內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置．

波浪能發(fā)電裝置的動力學(xué)研究是分析能量轉(zhuǎn)換過程、掌握發(fā)電規(guī)律、優(yōu)化控制策略的關(guān)鍵．王冬姣等[12]研究了含慣性擺波浪能發(fā)電裝置的浮標(biāo)的動力學(xué)響應(yīng)，分析了負(fù)載阻尼系數(shù)對浮標(biāo)運(yùn)動及波浪能轉(zhuǎn)換效率的影響．Wei等[13]將光滑粒子流體動力學(xué)模型應(yīng)用于擺式波浪能發(fā)電裝置的水動力學(xué)分析和錨泊系統(tǒng)的動力學(xué)分析，發(fā)現(xiàn)將裝置的固有頻率設(shè)計為當(dāng)?shù)夭ɡ祟l率，能夠提高能量輸出效果．Arzaghi等[14]基于點吸收式波浪能發(fā)電裝置的流體動力學(xué)分析結(jié)果，構(gòu)造了馬爾可夫鏈模型，估計輸出功率的穩(wěn)態(tài)分布，實現(xiàn)了對波浪能發(fā)電裝置發(fā)電功率的長期準(zhǔn)確預(yù)測．Zhan等[15]提出了一種適用于波浪能發(fā)電裝置的自適應(yīng)分層模型預(yù)測控制策略，能夠解決由于海況變化而引起的動力學(xué)建模的不準(zhǔn)確性問題．

目前，進(jìn)行波浪能發(fā)電裝置的動力學(xué)研究時，一般先計算漂浮式結(jié)構(gòu)物的水動力學(xué)響應(yīng)，再將水動力學(xué)響應(yīng)施加于能量轉(zhuǎn)換系統(tǒng)，從而獲得能夠輸出的電能數(shù)據(jù)．但是，漂浮式結(jié)構(gòu)物的水動力學(xué)與能量轉(zhuǎn)換系統(tǒng)的機(jī)構(gòu)動力學(xué)存在耦合關(guān)系，現(xiàn)有的研究思路無法獲得水動力學(xué)與機(jī)構(gòu)動力學(xué)的耦合特性，導(dǎo)致動力學(xué)分析結(jié)果不精確．特別是研發(fā)新型內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置時，掌握其動力學(xué)特性至關(guān)重要．

本文首先介紹了內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的結(jié)構(gòu)及工作原理，然后提出了一種基于ABAQUS軟件和XFlow軟件的動力學(xué)耦合分析新方法，最后分析了偏心轉(zhuǎn)子質(zhì)量、摩擦系數(shù)、波浪參數(shù)等因素對偏心轉(zhuǎn)子運(yùn)動特性的影響，并利用原理樣機(jī)模型開展了相關(guān)實驗．其中，本文的主要創(chuàng)新性為提出了海況自適應(yīng)液壓PTO方案，建立了動力學(xué)耦合分析方法，能夠同步開展內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的水動力學(xué)分析和機(jī)構(gòu)動力學(xué)分析，提高了動力學(xué)分析結(jié)果的精度，有利于新型內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的研發(fā)及性能評估．

1?結(jié)構(gòu)及工作原理

1.1?機(jī)械結(jié)構(gòu)

內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的機(jī)械結(jié)構(gòu)如圖1所示．

圖1?內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的機(jī)械結(jié)構(gòu)

上殼體和下殼體形成封閉的內(nèi)部空間，將安裝于內(nèi)部的機(jī)械元件與外部海水完全隔離，能夠有效避免機(jī)械元件受到海水腐蝕．液壓能量輸出系統(tǒng)(power-take-off，PTO)固定在支架和下殼體構(gòu)成的框架結(jié)構(gòu)內(nèi)部，并盡量降低液壓PTO的安裝高度，從而降低波浪能發(fā)電裝置的重心．偏心轉(zhuǎn)子為捕能元件，可以利用液壓缸主動調(diào)節(jié)偏心轉(zhuǎn)子的偏心距，以適應(yīng)不同海況條件．在波浪激勵作用下，偏心轉(zhuǎn)子驅(qū)動中心軸旋轉(zhuǎn)，中心軸獲得的轉(zhuǎn)速經(jīng)過增速箱增速后，輸入到液壓PTO．

1.2?液壓能量輸出系統(tǒng)

內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的液壓能量輸出系統(tǒng)如圖2所示．

1—油箱；2—過濾器；3—普通單向閥；4—雙向液壓泵；5—增速箱；6—帶彈簧的單向閥；7—壓力表；8—蓄能器；9—可變節(jié)流閥；10—液壓馬達(dá)；11—發(fā)電機(jī)；12—流量計；13—溢流閥；14—帶彈簧的液壓缸；15—偏心轉(zhuǎn)子

增速箱的輸出軸驅(qū)動雙向液壓泵旋轉(zhuǎn)，為液壓PTO提供動力輸入．在雙向液壓泵的吸油口使用普通單向閥，減少吸油口處的壓力損失．在雙向液壓泵的排油口使用帶彈簧的單向閥，確保閥芯能夠快速復(fù)位．可變節(jié)流閥能夠調(diào)節(jié)進(jìn)入液壓馬達(dá)的液壓油的流量，溢流閥能夠控制液壓PTO的最大工作壓力.蓄能器起到削峰填谷的作用，使作用于液壓馬達(dá)的液壓油保持壓力平穩(wěn)．帶彈簧的液壓缸能夠根據(jù)液壓系統(tǒng)的壓力，自適應(yīng)調(diào)節(jié)偏心轉(zhuǎn)子的偏心距．

1.3?海況自適應(yīng)發(fā)電原理

受波浪的激勵，內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置發(fā)生晃動，使得偏心轉(zhuǎn)子的重心位置偏離重力勢能最低點．在重力的作用下，偏心轉(zhuǎn)子向重力勢能最低點運(yùn)動，從而帶動中心軸轉(zhuǎn)動．增速箱能夠提升中心軸的轉(zhuǎn)動速度，并驅(qū)動液壓泵轉(zhuǎn)動．

中心軸既可以沿順時針方向旋轉(zhuǎn)，也可以沿逆時針方向旋轉(zhuǎn)，因此，液壓泵既可以正向旋轉(zhuǎn)，也可以反向旋轉(zhuǎn)．液壓泵輸出的高壓液壓油進(jìn)入液壓馬達(dá)，使液壓馬達(dá)驅(qū)動發(fā)電機(jī)發(fā)電．

當(dāng)波浪較大時，中心軸的旋轉(zhuǎn)速度較快，使得液壓泵的旋轉(zhuǎn)速度也較快，液壓泵輸出液壓油的流量較大，在不改變可變節(jié)流閥的節(jié)流面積的條件下，可變節(jié)流閥入口處的壓力較大，從而使進(jìn)入液壓缸有桿腔的液壓油壓力增大，高壓液壓油使彈簧壓縮，驅(qū)動液壓缸的活塞桿回縮，從而減小了偏心轉(zhuǎn)子的偏心距，進(jìn)而能夠降低中心軸受波浪激勵后的旋轉(zhuǎn)速度．

當(dāng)波浪較小時，中心軸的旋轉(zhuǎn)速度較慢，使得液壓泵的旋轉(zhuǎn)速度也較慢，液壓泵輸出液壓油的流量較小，在不改變可變節(jié)流閥的節(jié)流面積的條件下，可變節(jié)流閥入口處的壓力較小，從而使進(jìn)入液壓缸有桿腔的液壓油壓力減小，在彈簧力的作用下，使液壓缸的活塞桿伸出，從而增大了偏心轉(zhuǎn)子的偏心距，進(jìn)而能夠提高中心軸受波浪激勵后的旋轉(zhuǎn)速度．

此外，通過調(diào)節(jié)可變節(jié)流閥的節(jié)流面積，能夠獲得不同的液壓傳動性能．

因此，內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置能夠根據(jù)波浪狀況自動調(diào)節(jié)偏心轉(zhuǎn)子的偏心距，適應(yīng)復(fù)雜多變的海況條件，有效提高發(fā)電效率和發(fā)電穩(wěn)定性．

2?聯(lián)合仿真模型

2.1?ABAQUS的參數(shù)設(shè)置

為便于數(shù)值模擬分析，對內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的模型進(jìn)行簡化，如圖3所示．

圖3?幾何模型簡化

將上殼體和下殼體合并為殼體，將殼體、中心軸、偏心轉(zhuǎn)子分別導(dǎo)入到ABAQUS軟件，并在ABAQUS軟件中繪制系泊纜模型．對各零件分別進(jìn)行網(wǎng)格劃分，將各零件的網(wǎng)格模型組合為裝配體，設(shè)置殼體與流場域的流固耦合交界面以及偏心轉(zhuǎn)子與中心軸之間的摩擦接觸條件．

仿真模型中各零件的參數(shù)設(shè)置和初步劃分的網(wǎng)格單元數(shù)量如表1所示．

表1?幾何模型的參數(shù)設(shè)置

Tab.1?Parameter setting of geometric model

按照最大時間步長不大于最小網(wǎng)格尺寸與最大運(yùn)動速度之間比值的原則，同時考慮計算成本，設(shè)置時間步長為0.0002s，并導(dǎo)出幾何模型的.stl格式文件和仿真模型的.inp格式文件．

2.2?XFlow的參數(shù)設(shè)置

XFlow采用了基于粒子法的格子波爾茲曼技術(shù)，不需要對流場域進(jìn)行網(wǎng)格劃分．在XFlow中，采用線性波模型，流場域環(huán)境變量和波浪參數(shù)的設(shè)置如表2所示．

表2?流場域環(huán)境變量和波浪參數(shù)

Tab.2 Environmental variables and wave parameters of flow field

將ABAQUS中導(dǎo)出的幾何模型的.stl格式文件導(dǎo)入XFlow，設(shè)置Structural coupling為two way，設(shè)置時間步長為0.0002s，進(jìn)行雙向耦合分析．

2.3?建立耦合模型

首先運(yùn)行XFlow，實現(xiàn)流場域的初始化；然后采用命令行方式運(yùn)行ABAQUS，實現(xiàn)機(jī)構(gòu)動力學(xué)分析的初始化；使用相同的仿真時間和時間步長，實現(xiàn)ABAQUS和XFlow的同步迭代求解；最后定義數(shù)據(jù)傳輸接口，ABAQUS向XFlow傳遞機(jī)構(gòu)動力學(xué)參數(shù)，XFlow向ABAQUS傳遞水動力學(xué)激勵參數(shù)，實現(xiàn)ABAQUS與XFlow的聯(lián)合仿真，如圖4所示．

聯(lián)合仿真為瞬態(tài)模擬過程，設(shè)置總模擬時間為5s，達(dá)到設(shè)定的模擬時間后，仿真過程停止．

在ABAQUS中設(shè)置不同的全局種子間距，分別獲得網(wǎng)格單元總數(shù)為12455、15248、22695、36109、43989的網(wǎng)格模型．在ABAQUS-XFlow聯(lián)合仿真模型中，設(shè)置入射波速度為1m/s、波高為0.2m、波浪頻率為0.25Hz、偏心轉(zhuǎn)子的密度為1000kg/m3、偏心轉(zhuǎn)子與中心軸之間的摩擦系數(shù)為0.1，進(jìn)行網(wǎng)格無關(guān)性驗證．不同網(wǎng)格模型條件下偏心轉(zhuǎn)子轉(zhuǎn)角的時歷曲線如圖5所示．

圖4?ABAQUS與XFlow的聯(lián)合仿真

圖5?網(wǎng)格質(zhì)量對轉(zhuǎn)動角度的影響

當(dāng)網(wǎng)格單元總數(shù)分別為12455、15248和22695時，偏心轉(zhuǎn)子轉(zhuǎn)角的時歷曲線變化規(guī)律有較大差別.而當(dāng)網(wǎng)格單元總數(shù)分別為22695、36109和43989時，偏心轉(zhuǎn)子轉(zhuǎn)角的時歷曲線變化規(guī)律相近．說明網(wǎng)格單元總數(shù)超過22695后，繼續(xù)加密網(wǎng)格對數(shù)值模擬結(jié)果的影響較小．因此，為確保較高的數(shù)值模擬精度，選取網(wǎng)格單元總數(shù)為36109的網(wǎng)格模型，各零件的網(wǎng)格單元數(shù)量參見表1．

3?結(jié)果與討論

3.1?偏心轉(zhuǎn)子質(zhì)量的影響

在XFlow中，設(shè)置入射波速度為1m/s，波高為0.2m，波浪頻率為0.25Hz．在ABAQUS中，設(shè)置偏心轉(zhuǎn)子與中心軸之間的摩擦系數(shù)為0.1．為便于分析偏心轉(zhuǎn)子質(zhì)量對其運(yùn)動特性的影響，在保持體積不變的情況下，分別設(shè)置偏心轉(zhuǎn)子的密度為800kg/m3、1200kg/m3、2000kg/m3、4000kg/m3，對應(yīng)的偏心轉(zhuǎn)子質(zhì)量分別為3.83kg、5.74kg、9.57kg、19.13kg，分析不同偏心轉(zhuǎn)子質(zhì)量對其運(yùn)動性能的影響，如圖6和圖7所示．在工程實踐中，可以通過不同材料的組合，調(diào)配出對應(yīng)的材料密度．

圖6?偏心轉(zhuǎn)子質(zhì)量對轉(zhuǎn)動角度的影響

圖7?偏心轉(zhuǎn)子質(zhì)量對轉(zhuǎn)動速度的影響

圖6和圖7中，轉(zhuǎn)角和轉(zhuǎn)動速度的正值表明偏心轉(zhuǎn)子沿順時針方向旋轉(zhuǎn)，負(fù)值表明偏心轉(zhuǎn)子沿逆時針方向旋轉(zhuǎn)．偏心轉(zhuǎn)子的轉(zhuǎn)角越大、轉(zhuǎn)動速度越快，說明波浪能發(fā)電系統(tǒng)的功率輸出特性越好．

3.2?摩擦系數(shù)的影響

偏心轉(zhuǎn)子處的機(jī)械摩擦、液壓傳動系統(tǒng)的阻尼、發(fā)電機(jī)的輸出負(fù)載等均是影響波浪能捕獲效率的重要因素．為便于數(shù)值模擬研究，將上述因素統(tǒng)一轉(zhuǎn)化為偏心轉(zhuǎn)子處的摩擦力．

在XFlow中，設(shè)置入射波速度為1m/s，波高為0.2m，波浪頻率為0.25Hz．在ABAQUS中，設(shè)置偏心轉(zhuǎn)子的密度為1000kg/m3，對應(yīng)的質(zhì)量為4.78kg．同時，分別設(shè)置偏心轉(zhuǎn)子與中心軸之間的摩擦系數(shù)為0.05、0.10、0.15、0.20，分析不同摩擦系數(shù)對偏心轉(zhuǎn)子運(yùn)動性能的影響，如圖8和圖9所示．

圖8?摩擦系數(shù)對轉(zhuǎn)動角度的影響

圖9?摩擦系數(shù)對轉(zhuǎn)動速度的影響

3.3?波浪參數(shù)的影響

在ABAQUS中，設(shè)置偏心轉(zhuǎn)子與中心軸之間的摩擦系數(shù)為0.1，設(shè)置偏心轉(zhuǎn)子的密度為1000kg/m3，對應(yīng)的質(zhì)量為4.78kg．在XFlow中，設(shè)置入射波速度為0.5m/s，設(shè)置不同的波高和頻率fr，分析波浪參數(shù)對偏心轉(zhuǎn)子運(yùn)動性能的影響，如圖10和圖11所示．

圖10?波浪參數(shù)對轉(zhuǎn)動角度的影響

圖11?波浪參數(shù)對轉(zhuǎn)動速度的影響

在相同的摩擦系數(shù)和偏心轉(zhuǎn)子質(zhì)量條件下，波浪參數(shù)對偏心轉(zhuǎn)子運(yùn)動特性的影響規(guī)律呈現(xiàn)出一定的復(fù)雜性．保持波高恒定，逐漸加大波浪頻率，偏心轉(zhuǎn)子轉(zhuǎn)動角度的上限值先增大，后減??；而下限值有減小的趨勢；但總幅值有逐漸增大的趨勢．偏心轉(zhuǎn)子的轉(zhuǎn)動速度也有相似的變化規(guī)律．保持波浪頻率恒定，增大波高，偏心轉(zhuǎn)子轉(zhuǎn)動角度的上限值和下限值均有增大的趨勢，但偏心轉(zhuǎn)子的轉(zhuǎn)動速度幅值變化不明顯．此外，波高0.50m、頻率0.5Hz和波高0.75m、頻率0.1Hz兩種波浪狀況下，單位波前寬度上的波浪能功率不同，但能夠獲得相似的偏心轉(zhuǎn)子運(yùn)動特性．

4?原理樣機(jī)實驗

為進(jìn)一步探索內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的偏心轉(zhuǎn)子動力學(xué)特性，建立了其原理實驗?zāi)Ｐ?，如圖12所示．

圖12?內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的原理實驗?zāi)Ｐ?/p>

與ABAQUS-XFlow聯(lián)合仿真中的簡化模型相似，按照1∶5的縮小比例，采用相同的封閉式殼體結(jié)構(gòu)將捕能機(jī)構(gòu)與流體介質(zhì)隔離．捕能機(jī)構(gòu)的偏心轉(zhuǎn)子與中心軸固定連接，中心軸又通過聯(lián)軸器與編碼器的輸入軸固定連接．因此，編碼器能夠記錄偏心轉(zhuǎn)子的轉(zhuǎn)動角度．同時，將姿態(tài)傳感器和數(shù)據(jù)采集卡固定在波浪能發(fā)電裝置殼體內(nèi)部，從而能夠在線記錄波浪能發(fā)電裝置的姿態(tài)數(shù)據(jù)和偏心轉(zhuǎn)子的運(yùn)動特性．

所選用的數(shù)據(jù)采集卡的采樣頻率為10Hz，編碼器與姿態(tài)傳感器的性能如表3所示．

表3?編碼器與姿態(tài)傳感器的性能參數(shù)

Tab.3 Performance parameters of encoder and attitude sensor

與ABAQUS-XFlow聯(lián)合仿真中的簡化模型不同，受實驗條件限制，中心軸轉(zhuǎn)動時的摩擦系數(shù)未能定量測量，造波器也無法精確模擬XFlow中的線性波模型，所產(chǎn)生的波浪激勵條件存在較大的隨機(jī)性，因而實驗數(shù)據(jù)僅可用于定性分析偏心轉(zhuǎn)子的運(yùn)動特性．

圖13?波浪能發(fā)電裝置沿著x軸、y軸、z軸加速度

圖14?波浪能發(fā)電裝置沿著x軸、y軸、z軸角加速度

圖15?偏心轉(zhuǎn)子的運(yùn)動情況

5?結(jié)?語

本文介紹了一種內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置，提出了海況自適應(yīng)液壓PTO方案，能夠根據(jù)海況特點自適應(yīng)調(diào)節(jié)偏心轉(zhuǎn)子的偏心距，維持較高的能量捕獲效率．建立了適用于內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置的ABAQUS-XFlow聯(lián)合仿真模型，能夠同步開展水動力學(xué)分析和機(jī)構(gòu)動力學(xué)分析，提高了動力學(xué)分析結(jié)果的精度．

通過數(shù)值模擬和實驗分析發(fā)現(xiàn)：偏心轉(zhuǎn)子的質(zhì)量對其運(yùn)動特性的影響呈現(xiàn)出非單調(diào)性規(guī)律，是影響運(yùn)動特性的敏感因素；偏心轉(zhuǎn)子與中心軸之間的摩擦系數(shù)對偏心轉(zhuǎn)子運(yùn)動特性的影響呈現(xiàn)出單調(diào)性規(guī)律，但偏心轉(zhuǎn)子的運(yùn)動特性對偏心轉(zhuǎn)子與中心軸之間的摩擦系數(shù)的敏感度較低；波浪能發(fā)電裝置受到隨機(jī)波浪激勵后，在垂蕩和縱搖方向的運(yùn)動幅值較大，在其他方向上的運(yùn)動幅值相對較小；偏心轉(zhuǎn)子有較大概率保持沿單一方向轉(zhuǎn)動；偏心轉(zhuǎn)子轉(zhuǎn)動較為平穩(wěn)，能夠穩(wěn)定地捕獲波浪能．

本文為波浪能發(fā)電裝置的水動力學(xué)和機(jī)構(gòu)動力學(xué)耦合分析提供了新方法，也為提高波浪能發(fā)電裝置的海況自適應(yīng)發(fā)電能力提供了新思路．

受數(shù)值模擬方法和所研制的波浪能發(fā)電裝置原理樣機(jī)功能的限制，無法直接獲得發(fā)電裝置的輸出功率特性，僅能夠利用偏心轉(zhuǎn)子的運(yùn)動特性間接反映輸出功率特性．在未來的研究工作中，將研制具有海況自適應(yīng)液壓PTO的內(nèi)置偏心轉(zhuǎn)子式波浪能發(fā)電裝置，利用實驗分析輸出功率特性，并驗證其海況自適應(yīng)發(fā)電能力．

［1］ Jeff T. Power from the oceans：Blue energy[J]. Nature，2014，508(7496)：302-304.

［2］ 鄭崇偉，李崇銀. 海洋強(qiáng)國視野下的“海上絲綢之路”海洋新能源評估[J]. 哈爾濱工程大學(xué)學(xué)報，2020，41(2)：175-183.

Zheng Chongwei，Li Chongyin. Evaluation of new marine energy for the Maritime Silk Road from the perspective of maritime power[J]. Journal of Harbin Engineering University，2020，41(2)：175-183(in Chinese).

［3］ Esteban M，Leary D. Current developments and future prospects of offshore wind and ocean energy[J]. Applied Energy，2012，90(1)：128-136.

［4］ Aderinto T，Li H. Review on power performance and efficiency of wave energy converters[J]. Energies，2019，12(22)：4329-1-4329-24.

［5］ 唐友剛，曲志森，趙?青，等. 新型浮式浮力擺波浪能裝置運(yùn)動特性試驗研究[J]. 天津大學(xué)學(xué)報(自然科學(xué)與工程技術(shù)版)，2017，50(10)：1029-1036.

Tang Yougang，Qu Zhisen，Zhao Qing，et al. Test study of kinetic characteristic of a new floating flap wave energy converter[J]. Journal of Tianjin University(Sci-ence and Technology)，2017，50(10)：1029-1036(in Chinese).

［6］ 盛松偉，王坤林，吝紅軍，等. 100kW鷹式波浪能發(fā)電裝置“萬山號”實海況試驗[J]. 太陽能學(xué)報，2019(3)：709-714.

Sheng Songwei，Wang Kunlin，Lin Hongjun，et al. Open sea tests of 100kW wave energy convertor sharp eagle Wanshan[J]. Acta Energiae Solaris Sinica，2019(3)：709-714(in Chinese).

［7］ Xie Q，Zhang T，Pan Y，et al. A novel oscillating buoy wave energy harvester based on a spatial double X-shaped mechanism for self-powered sensors in sea-crossing bridges[J]. Energy Conversion and Management，2020，204：112286-1-112286-14.

［8］ Zhang D，Li W，Zhao H，et al. Design of a hydraulic power take-off system for the wave energy device with an inverse pendulum[J]. China Ocean Engineering，2014，28(2)：283-292.

［9］ Hinck Ⅲ. Wave Power Generator：US，3231749[P]. 1963-04-12.

［10］ Chen H M，DelBalzo D R. Eccentric rotating wave energy converter[C]//2019 IEEE Underwater Technology(UT). Kaohsiung，Taiwan，China，2019：1-5.

［11］ Yurchenko D，Alevras P. Parametric pendulum based wave energy converter[J]. Mechanical Systems and Signal Processing，2018，99：504-515.

［12］ 王冬姣，王英毅，邱守強(qiáng)，等. 含慣性擺發(fā)電裝置之浮標(biāo)的水動力分析[J]. 船舶力學(xué)，2019，23(4)：46-55.

Wang Dongjiao，Wang Yingyi，Qiu Shouqiang，et al. Hydrodynamics of a buoy with inertial pendulum power generation device[J]. Journal of Ship Mechanics，2019，23(4)：46-55(in Chinese).

［13］ Wei Z，Edge B L，Dalrymple R A，et al. Modeling of wave energy converters by GPUSPH and project chrono[J]. Ocean Engineering，2019，183：332-349.

［14］ Arzaghi E，Abaei M M，Abbassi R，et al. A Markovian approach to power generation capacity assessment of floating wave energy converters[J]. Renewable Energy，2020，146：2736-2743.

［15］ Zhan S，Na J，Li G，et al. Adaptive model predictive control of wave energy converters[J]. IEEE Transactions on Sustainable Energy，2018，11(1)：229-238.

Dynamic Study of the Wave Energy Converter with Inner Eccentric Rotor

Xue Gang1, 2，Liu Yanjun1, 2，Xue Yifan1，Liu Dahui3

(1. Institute of Marine Science and Technology，Shandong University，Qingdao 266237，China；2. School of Mechanical Engineering，Shandong University，Jinan 250061，China；3. CIMC Offshore Engineering Institute Co.，Ltd.，Yantai 264670，China)

The novel wave energy converter with an inner eccentric rotor is a totally enclosed structure，which is characterized by high reliability and long working life. The wave-condition-adaptive power-take-off scheme is proposed for the wave energy converter with the inner eccentric rotor. The scheme will change the eccentric distance adaptively according to the wave conditions and improve the energy output efficiency and stability. The research on the dynamics of wave energy converter with inner eccentric rotor involves two aspects，namely，the hydrodynamics of the entire equipment and the mechanism dynamics of the inner eccentric rotor. A coupling analysis method for hydrodynamics and mechanism dynamics based on ABAQUS and XFlow joint simulation model is established. The factors that influence the motion characteristics of the eccentric rotor，such as eccentric rotor mass，friction coefficient，and wave parameters，are analyzed. Moreover，the motion characteristics of the eccentric rotor under random wave conditions are experimentally analyzed using the principle prototype. The motion characteristics of the eccentric rotor could reflect the power output characteristics indirectly. Notably，the eccentric rotor mass is a sensitive factor that affects the motion characteristics of the eccentric rotor，and its influence is nonmonotonic. The motion characteristics of the eccentric rotor are less sensitive to the friction coefficient between the eccentric rotor and the central shaft，and the influence of the friction coefficient is monotonic. When excited by a random wave，the wave energy converter obtains a larger motion amplitude in the heaving and pitching directions than that in other directions. The eccentric rotor has a greater probability of keeping rotating in a single direction，and the rotational motion of eccentric rotor is relatively stable. Thus，the eccentric rotor can stably capture wave energy. This paper can be an important reference for developing the wave energy converter with the inner eccentric rotor.

eccentric rotor；wave energy converter；dynamics；coupling analysis

10.11784/tdxbz202010018

TK79

A

0493-2137(2022)02-0191-08

2020-10-12；

2021-01-04.

薛?鋼（1990—??），男，博士，副研究員，xuegangzb@163.com.Email：m_bigm@tju.edu.cn

劉延俊，lyj111ky@163.com.

國家重點研發(fā)計劃資助項目(2017YFE0115000，2016YFE0205700)；國家自然科學(xué)基金資助項目(U1706230).

the National Key Research and Development Program of China(No.2017YFE0115000，No.2016YFE0205700)，the National Natural Science Foundation of China(No.U1706230).

(責(zé)任編輯：王曉燕)