王聞航, 儲景瑞, 盧 超, 陳 偉, 楊 鋼
(1.華中科技大學(xué)機(jī)械科學(xué)與工程學(xué)院FESTO 氣動中心, 湖北武漢 430074;2.中國艦船研究設(shè)計中心, 湖北武漢 430064)
防波板位于潛航器表面,航行時可保持航行器流線外形,以減小阻力[1];而工作時,需開啟防波板,以保證工作空間。由于潛航器內(nèi)部空間十分緊張,需對某現(xiàn)有防波板啟閉機(jī)構(gòu)進(jìn)行優(yōu)化[2],減小其開啟過程中的最大驅(qū)動力矩[3],從而降低防波板機(jī)構(gòu)驅(qū)動器尺寸,同時降低開啟功率,減少啟閉過程的能源消耗。
防波板啟閉機(jī)構(gòu)由電機(jī)驅(qū)動,通過連桿機(jī)構(gòu)傳遞動力,驅(qū)動防波板開啟和關(guān)閉。對于連桿機(jī)構(gòu)中力矩優(yōu)化的研究[4-7]較多,劉杰等[8]對載具連桿機(jī)構(gòu)展開狀態(tài)下的力矩進(jìn)行了優(yōu)化,優(yōu)化后力矩下降60%,取得較好優(yōu)化效果。但目前對于連桿機(jī)構(gòu)運動全過程的優(yōu)化研究較少,全過程動態(tài)優(yōu)化[9]即綜合考慮整個運動過程的受力、力矩變化,能夠在一定程度上提高優(yōu)化效果。
本研究對防波板的啟閉連桿機(jī)構(gòu)進(jìn)行全過程動態(tài)優(yōu)化設(shè)計,以啟閉過程中所需的最大驅(qū)動力矩為優(yōu)化目標(biāo),結(jié)合Fluent流體仿真[10]獲得啟閉過程防波板所受負(fù)載力的開度-水動力函數(shù),通過ADAMS建立動力學(xué)仿真模型以獲取驅(qū)動力矩變化特性[11],最后通過控制變量法實現(xiàn)了目標(biāo)力矩優(yōu)化。
防波板啟閉裝置三維結(jié)構(gòu)如圖1所示,通過四桿機(jī)構(gòu)帶動防波板完成啟閉動作。
圖1 防波板啟閉裝置三維模型Fig.1 3D model of shutter opening and closing device
機(jī)構(gòu)的平面示意圖如圖2所示。啟閉連桿機(jī)構(gòu)曲柄OA桿初末位置固定,連桿BC與防波板固連,防波板繞C點旋轉(zhuǎn)40°,即∠BCB1。優(yōu)化參數(shù)為鉸接點B與鉸接點C橫縱坐標(biāo)。機(jī)構(gòu)的約束條件為:
圖2 平面連桿機(jī)構(gòu)示意圖Fig.2 Schematic diagram of planar linkage mechanism
(1) 防波板開啟后需滿足工作空間;
(2) 防波板啟閉過程中無死點。
由上述約束條件分別得B,C兩點坐標(biāo)的約束方程如式(1)、式(2)所示,式中單位為mm。
(1)
(2)
由該方程組,可得出C點和B點的約束范圍,分別如圖2三角形陰影和四邊形陰影所示。由于運動過程中,桿長應(yīng)保持不變,以C點為坐標(biāo)原點推導(dǎo)得到B點橫縱坐標(biāo)限制關(guān)系,如式(3)所示:
(3)
由式(3)可知,當(dāng)A,A1,C點位置確定時,B點橫縱坐標(biāo)滿足線性關(guān)系,如圖2所示。
為確定驅(qū)動力大小,需首先計算防波板所受負(fù)載力,而在航行過程中,防波板所受的負(fù)載力近似認(rèn)為等于航行過程中的水動力。
本研究以防波板為研究對象,基于Fluent開展水動力仿真分析。為真實模擬海水流過防波板的狀況,對防波板及附近的部分潛航器結(jié)構(gòu)進(jìn)行數(shù)值計算,計算模型如圖3所示。劃分網(wǎng)格約為160萬個并驗證了網(wǎng)格無關(guān)性[12],計算模型設(shè)定為k-ε中的RNG湍流模型[13]。
圖3 流體計算模型Fig.3 Fluid computational model
為綜合考慮防波板啟閉全過程的負(fù)載力,將防波板繞旋轉(zhuǎn)點C開啟角度從0°~50°均分為11種工況進(jìn)行數(shù)值計算。
針對0°~50°間的11種工況,分別開展了仿真計算,得出了防波板所受水動力及表面應(yīng)力。其中,防波板全閉(∠BCB1=0°)、半開(∠BCB1=25°)、全開(∠BCB1=50°)工況壓力云圖,如圖4所示。
圖4 防波板壓力云圖Fig.4 Pressure contour map of shutter
防波板所受水動力可分解為x,y,z方向,而z方向受力被連桿結(jié)構(gòu)平衡,防波板負(fù)載力主要為x,y方向所受水動力。防波板所受x和y方向水動力,如表1所示。
表1 防波板水動力Tab.1 Hydrodynamics of shutter
基于上述仿真結(jié)果擬合得出防波板開度-水動力函數(shù),如式(4)所示:
(4)
式中,Fx—— 防波板x方向受力
Fy—— 防波板y方向受力,N
x—— 防波板繞下旋轉(zhuǎn)點旋轉(zhuǎn)角度,(°)
防波板x方向受力隨開啟角度變大先變小和增大,在開啟角度為40°,即防波板與航行方向平行時防波板水動力約為0。y方向受力隨開啟角度變大而變大,符合實際情況。
因防波板開啟特定角度后,其位置與旋轉(zhuǎn)點C位置相關(guān),故需驗證防波板開度-水動力函數(shù)在C點位置不同時的通用性??紤]到旋轉(zhuǎn)點C可變范圍較小,在C點約束范圍內(nèi),對3個邊界頂點位置分別進(jìn)行半開、全開兩種工況仿真,仿真結(jié)果如表2所示。
表2 三極限位置仿真結(jié)果Tab.2 Simulation results of three extreme positions
由表可見,旋轉(zhuǎn)點C位置對防波板水動力無明顯影響,可認(rèn)為不同的旋轉(zhuǎn)點位置下的開度-水動力函數(shù)相同。
防波板啟閉連桿機(jī)構(gòu)為平面連桿機(jī)構(gòu),基于ADAMS[14-15]搭建動力學(xué)模型并開展仿真分析,得出啟閉連桿機(jī)構(gòu)驅(qū)動防波板啟閉過程所需驅(qū)動力矩的變化特性。
動力學(xué)模型如圖5所示,從左至右分別為點O,A,B,C。OA桿為驅(qū)動桿,在O鉸接點處施加位移驅(qū)動,使得防波板經(jīng)歷開啟(0~5 s)、暫位(5~7 s)、關(guān)閉(7~12 s)過程,開啟階段驅(qū)動桿OA勻速運動至∠ACA1=80°,暫位階段驅(qū)動桿OA保持靜止,關(guān)閉階段驅(qū)動桿OA勻速運動回原位;在BC桿對應(yīng)防波板受力中心位置施加水動力,該力作用點分別添加x,y方向集中力,力函數(shù)為開度-水動力函數(shù)。由于防波板開啟和關(guān)閉過程為勻速、慢速啟閉過程,可認(rèn)為驅(qū)動桿OA所受力矩等于克服水動力負(fù)載所需力矩。
圖5 ADAMS動力學(xué)模型圖Fig.5 ADAMS dynamic model diagram
得到驅(qū)動力矩T以及功率P變化圖,如圖6所示。由圖6可知,啟閉連桿機(jī)構(gòu)所需驅(qū)動力矩在開啟過程中逐漸增大,在接近防波板呈開啟狀態(tài)時,所需力矩急劇增大后達(dá)到最大值,關(guān)閉過程與開啟過程變化規(guī)律對稱,啟閉連桿機(jī)構(gòu)在開啟狀態(tài)時的傳力性能對最大驅(qū)動力矩有重要影響。
圖6 ADAMS仿真力矩圖Fig.6 ADAMS dynamic model diagram
優(yōu)化參數(shù)為B點坐標(biāo)與C點坐標(biāo),為得出B,C點坐標(biāo)的優(yōu)化解,基于控制變量法探究兩鉸接點位置對優(yōu)化目標(biāo)的影響。
控制C點坐標(biāo)不變,O,A,C點坐標(biāo)分別設(shè)置為O(0,0),A(464.363,24.336),C(2473.803,-399.479),坐標(biāo)單位為mm。在B點限制范圍內(nèi)更改B點坐標(biāo)進(jìn)行動力學(xué)仿真,不同位置連桿機(jī)構(gòu)示意圖如圖7所示,從左至右依次為位置1~5,重復(fù)仿真過程,得到仿真結(jié)果曲線圖如圖8所示,其中,由于開啟后,連桿機(jī)構(gòu)存在死點,B點位于位置1時,仿真發(fā)散。
圖7 5組不同B點位置示意圖Fig.7 Schematic diagram of five sets of different B point positions
圖8 不同B點位置仿真結(jié)果圖Fig.8 Simulation result plot for different B point positions
圖9 3組不同C點位置示意圖Fig.9 Schematic diagram of three sets of different C point positions
由圖8可知,驅(qū)動力矩整體變化趨勢相同,驅(qū)動力矩隨防波板開啟角度增加而增大,不同B點坐標(biāo)對應(yīng)的最大力矩及功率值見表3。
表3 B點坐標(biāo)對最大力矩影響結(jié)果Tab.3 Results of influence of B point coordinates on maximum torque
由表3可知,B點位置越靠近防波板,即橫縱坐標(biāo)越大,啟閉全過程中最大力矩越小,最大力矩的最小值為11040 N·m。此時,B點位置在防波板開啟狀態(tài)下∠A1B1C最小,即在現(xiàn)有約束下傳動角最小,傳力性能較好,仿真結(jié)果符合機(jī)械傳動相關(guān)理論。
因B點坐標(biāo)與C點坐標(biāo)相關(guān),無法控制B點坐標(biāo)不變,由上述計算過程可知,B點在其可行域右邊界時,最大驅(qū)動力矩最小。保持O點與A點坐標(biāo)不變,C點坐標(biāo)不同時,對應(yīng)的B點均選取其可行域右邊界上的點,不同位置連桿機(jī)構(gòu)示意圖如9所示,仿真結(jié)果曲線圖,如圖10所示。
圖10 不同C點位置仿真結(jié)果圖Fig.10 Simulation result plot for different C point positions
由仿真結(jié)果圖可知,驅(qū)動力矩隨防波板開啟角度增加而增大,不同C點坐標(biāo)對應(yīng)的最大力矩值見表4。
表4 C點坐標(biāo)對最大力矩影響結(jié)果Tab.4 Results of influence of C point coordinates on maximum torque
由表4可知,C點位于可行域三角形左下角時,最大驅(qū)動力矩最小,最大力矩的最小值為8811.6 N·m。此時,在防波板開啟狀態(tài)下,∠A1B1C最小,即在現(xiàn)有約束條件下,傳動角最小,傳力性能較好,仿真結(jié)果符合機(jī)械傳動相關(guān)理論。
綜上可知,在C點坐標(biāo)不變情況下,B點位置越靠近防波板,即約束范圍內(nèi)橫坐標(biāo)越大,最大驅(qū)動力矩越小;在保持B點位置處于右邊界條件下,C點橫、縱坐標(biāo)越小,對應(yīng)的最大驅(qū)動力矩越小,且對最大力矩的影響弱于B點位置對最大驅(qū)動力矩影響。
因此在現(xiàn)有約束下,可認(rèn)為C點處于約束范圍左下位,B點處于約束范圍右邊界上時,啟閉全過程最大驅(qū)動力矩最小。優(yōu)化前后啟閉全過程驅(qū)動力矩及功率變化特性如圖11所示,最大驅(qū)動力矩及功率值如表5所示。優(yōu)化后的最大驅(qū)動力矩由17934 N·m降為8811.6 N·m,最大驅(qū)動力矩從3.9468 kW下降為2.4366 kW;驅(qū)動力矩和驅(qū)動功率分別減少了50.86%和38.26%。
表5 防波板參數(shù)優(yōu)化對比表Tab.5 Comparison table of shutter parameter optimization
圖11 優(yōu)化前后對比圖Fig.11 Comparison graph before and after optimization
本研究以減小防波板啟閉機(jī)構(gòu)最大驅(qū)動力矩為目標(biāo),開展了防波板啟閉連桿機(jī)構(gòu)全過程動態(tài)優(yōu)化,得出如下結(jié)論:
(1) 通過Fluent仿真計算求得防波板在潛航器運動時開啟過程中所受負(fù)載力的變化,得到開度-水動力函數(shù);同時通過流體仿真驗證了開度-水動力函數(shù)的通用性;
(2) 基于ADAMS搭建了防波板啟閉機(jī)構(gòu)的動力學(xué)模型,仿真結(jié)果表明驅(qū)動力矩隨著防波板開度增加而增大;
(3) 基于控制變量法探究了優(yōu)化參數(shù)B,C點坐標(biāo)對最大驅(qū)動力矩影響,其中B點越靠近右邊界,C點越靠近左下邊界,對應(yīng)的最大驅(qū)動力矩越小;
(4) 優(yōu)化后最大驅(qū)動力矩從17934 N·m減少到8811.6 N·m,同比下降50.86%,與此同時,最大驅(qū)動功率從3.9468 kW下降為2.4366 kW,同比下降38.26%。
本研究優(yōu)化工作考慮了啟閉連桿機(jī)構(gòu)開啟關(guān)閉全過程的受力與驅(qū)動力矩變化,大幅減小了防波板啟閉機(jī)構(gòu)最大驅(qū)動力矩和驅(qū)動功率。這對減小驅(qū)動設(shè)備尺寸、節(jié)約設(shè)備能源有著重要意義。