亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于樣本船的興波阻力智能預(yù)報(bào)方法

        2022-08-17 11:24:04周廣利蔣彩霞
        船舶力學(xué) 2022年8期
        關(guān)鍵詞:興波船型船體

        陳 帥,周廣利,王 超,蔣彩霞

        (1.中國(guó)船舶科學(xué)研究中心,江蘇 無錫 214082;2.哈爾濱工程大學(xué) 船舶工程學(xué)院,哈爾濱 150001)

        0 引 言

        為了提高船型優(yōu)化設(shè)計(jì)系統(tǒng)的計(jì)算效率,近似模型的引入至關(guān)重要。目前以減小興波阻力性能為優(yōu)化目標(biāo)的船型設(shè)計(jì)問題較多,而船體曲面變形技術(shù)和興波阻力預(yù)報(bào)技術(shù)為近似模型訓(xùn)練提供樣本數(shù)據(jù),隨著計(jì)算機(jī)和人工智能技術(shù)的迅速發(fā)展,構(gòu)建基于樣本船型的興波阻力高精度近似計(jì)算模型來解決優(yōu)化過程中的海量計(jì)算具有重要研究意義。

        船體曲面變形技術(shù)需要解決的是如何以較少的設(shè)計(jì)變量對(duì)船體曲面進(jìn)行變換,本文根據(jù)自主開發(fā)的NURBS技術(shù)[1],采用改進(jìn)平移法和自由變形技術(shù)對(duì)船體曲面進(jìn)行修改生成海量的船型樣本方案。船體興波阻力快速預(yù)報(bào)可以通過求解定常興波問題得到,在上兩個(gè)世紀(jì)中,出現(xiàn)了一些經(jīng)典的興波勢(shì)流理論,如Michell 理論、Kelvin 源(或Havelock 源)、Neumann-Kelvin(NK)理論和新細(xì)長(zhǎng)體理論等,這些理論為后續(xù)的興波問題研究奠定了基礎(chǔ)。Dawson[2]在自由面和船體表面布置Rankine 源來求解定常興波問題,使得Rankine 源法在船舶興波阻力預(yù)報(bào)上取得了巨大的成功。本文開發(fā)了基于貼體網(wǎng)格[3]的Dawson方法計(jì)算程序,只需數(shù)十秒就可以計(jì)算一次興波,適用于樣本船型預(yù)報(bào)。

        在近似模型方面,馮佰威[4-5]綜述了近似技術(shù)在船型自動(dòng)優(yōu)化中的意義并進(jìn)行應(yīng)用;萬德成等[6]討論了近似模型在船型優(yōu)化中的應(yīng)用現(xiàn)狀以及研宄難點(diǎn);劉志強(qiáng)[7]以O(shè)PTShip-SJTU[12]求解器為基礎(chǔ),分析了不同近似模型的計(jì)算精度;Huang[8]提出了一種基于徑向基函數(shù)的近似模型降低了貨船優(yōu)化的計(jì)算成本;Vesting[9]采用近似模型與局部搜索方法相結(jié)合大幅度提高了船用螺旋槳的設(shè)計(jì)效率;Priftis[10]采用近似模型研究了不確定條件下船舶的參數(shù)化設(shè)計(jì)和多目標(biāo)優(yōu)化問題,在設(shè)計(jì)初期對(duì)總阻力、運(yùn)費(fèi)率和成本進(jìn)行優(yōu)化;Kim[11]采用響應(yīng)面Kriging模型優(yōu)化Series-60船型,結(jié)果表明響應(yīng)面法大大降低計(jì)算成本。為了實(shí)現(xiàn)興波阻力的智能預(yù)報(bào),本文分別對(duì)樣本船型生成、興波阻力預(yù)報(bào)和近似模型構(gòu)建展開研究,首先提出曲面變形技術(shù)的理論方法,以ITM 和FFD 兩種方法對(duì)船體變形進(jìn)行應(yīng)用。其次,給出興波理論和數(shù)值計(jì)算方法,驗(yàn)證本文興波阻力計(jì)算程序的可靠性。最后,建立基于SLE-BP 神經(jīng)網(wǎng)絡(luò)的興波阻力預(yù)報(bào)近似模型,利用KCS船型進(jìn)行智能預(yù)報(bào)方法的應(yīng)用和驗(yàn)證。

        1 曲面變形技術(shù)

        曲面變形技術(shù)是實(shí)現(xiàn)船型優(yōu)化設(shè)計(jì)的基礎(chǔ),同時(shí)為近似模型訓(xùn)練提供樣本船型。本文根據(jù)NURBS 技術(shù)[1],通過修改曲面控制網(wǎng)格可以實(shí)現(xiàn)曲面變形,該方法使得變形后的曲面仍為NURBS 曲面。如圖1所示,在船體曲面生成的過程中,采用ITM和FFD法對(duì)曲線控制頂點(diǎn)進(jìn)行修改,通過變形后的控制網(wǎng)格和基函數(shù)計(jì)算出變形后的船體曲面。

        圖1 曲面變形技術(shù)的程序?qū)崿F(xiàn)流程Fig.1 Program implementation flow of surface deformation technology

        1.1 改進(jìn)平移法

        根據(jù)修改函數(shù)式(1),以曲面v方向某曲線的變形為例,ITM 法原理示意圖如圖2 所示,f1(x)和f2(x)為變形前后的函數(shù)曲線,g(x)為修改函數(shù)曲線。針對(duì)船體曲面變形,以控制網(wǎng)格點(diǎn)作為輸入即可得到變形中的位移大小和方向,再通過曲面重構(gòu)就可以得到變形船體。以KCS船體曲面變形為例,從圖3可以看出,式(1)中x1、x2、α1和α2四個(gè)參數(shù)含義均有體現(xiàn),給定參數(shù)擾動(dòng)區(qū)間即可實(shí)現(xiàn)船體曲面的整體變形。ITM 法的特點(diǎn)是設(shè)計(jì)變量涉及較少,非常利于船體幾何(半船、整船)的修改,并且基于NURBS技術(shù)的船體曲面控制網(wǎng)格變換可以保證新船型的光順性和連續(xù)性。

        圖2 ITM法原理示意圖Fig.2 Schematic diagram of ITM method

        圖3 KCS船體曲面變形示意圖Fig.3 Deformation diagram of KCS hull surface

        1.2 自由變形方法

        自由變形(free-form deformation,FFD)方法[3]是Sederberg 和Parry 提出的一種網(wǎng)格自由變形方法。它作為常用的幾何造型技術(shù),目前被廣泛應(yīng)用于計(jì)算機(jī)動(dòng)畫、機(jī)器人導(dǎo)航和船體幾何重構(gòu)等領(lǐng)域。

        在一個(gè)含有待變形物體的長(zhǎng)方體中構(gòu)造局部坐標(biāo)系O'-STU,在笛卡爾坐標(biāo)系O-XYZ下的任意點(diǎn)X在局部坐標(biāo)系中的坐標(biāo)為(s,t,u),則有

        式中,B代表Bernstein多項(xiàng)式基函數(shù),定義為

        由式(5)可知,初始船體曲面控制網(wǎng)格點(diǎn)是長(zhǎng)方體內(nèi)控制頂點(diǎn)的線性函數(shù)。將部分控制頂點(diǎn)的位置作為設(shè)計(jì)變量,通過控制長(zhǎng)方體框架的變形來達(dá)到船體曲面控制網(wǎng)格變形的目的。長(zhǎng)方體框架內(nèi)的船體曲面控制網(wǎng)格點(diǎn)的局部坐標(biāo)為(s,t,u),部分控制頂點(diǎn)的位置變化將導(dǎo)致X移動(dòng)到X':

        如圖4 所示,以KCS 球艏變形為例,在船艏處布置控制框架確定曲面變形區(qū)域,利用等分?jǐn)?shù)設(shè)置框架內(nèi)控制點(diǎn)位置,控制點(diǎn)在x、y、z方向上的位移控制球艏曲面在三個(gè)方向上的變形,如圖4(a)所示,在控制框架內(nèi)的曲面控制網(wǎng)格隨著部分控制點(diǎn)的位移而變形,通過曲面正算可計(jì)算出曲面節(jié)點(diǎn)的變形,如圖4(b)所示,船艏縱剖線變形前后對(duì)比顯示,整個(gè)球艏區(qū)域變形節(jié)點(diǎn)保持協(xié)調(diào)性。

        圖4 KCS船體球艏變形Fig.4 KCS hull bulbous bow deformation

        2 興波阻力預(yù)報(bào)

        2.1 興波理論

        本文采用C++語言開發(fā),基于Dawson 法的興波阻力計(jì)算步驟[2],首先在船體表面上用面元法計(jì)算疊模擾動(dòng)流,然后在船體表面和自由面上利用物面邊界條件、自由邊界條件和輻射條件計(jì)算興波擾動(dòng)流。計(jì)算模型離散如圖5所示,把求解的疊模繞流和興波擾動(dòng)流問題轉(zhuǎn)化為求解船體表面和自由面分布源點(diǎn)強(qiáng)度問題,有

        圖5 計(jì)算模型離散示意圖Fig.5 Schematic diagram of discrete calculation model

        式中:Φ為總速度勢(shì),滿足拉普拉斯方程;φr為疊模擾動(dòng)速度勢(shì);φw為興波擾動(dòng)速度勢(shì);為疊模表面源強(qiáng);σB和σF分別為船體表面和自由面源強(qiáng);SB為船體表面;SF為自由面;為場(chǎng)點(diǎn)p到船體表面上點(diǎn)qb和自由面上點(diǎn)qw的距離。

        2.2 數(shù)值計(jì)算驗(yàn)證

        以Wigley 和KCS船型作為算例,圖6~7為自由面在各方向的加密增長(zhǎng)率示意圖,自由面采用貼體網(wǎng)格。通過計(jì)算域大小變化對(duì)結(jié)果的影響分析,確定自由面上游擴(kuò)展0.5L,下游擴(kuò)展2.5L,y方向延伸1.5L。如圖所示對(duì)船體首尾部及船寬方向進(jìn)行了加密。

        圖6 Wigley船型自由面網(wǎng)格示意圖Fig.6 Schematic diagram of Wigley ship free surface grid

        圖7 KCS船體及自由面網(wǎng)格示意圖Fig.7 Schematic diagram of KCS hull and free surface grid

        圖8(a)和圖9(a)為以自由面橫向網(wǎng)格數(shù)為變量的興波阻力系數(shù)收斂圖,從圖中可以看出,隨著面元數(shù)的增加興波阻力系數(shù)趨近一個(gè)值,因此Dawson方法具有較好的收斂性。圖8(b)和圖9(b)為兩種船型興波阻力系數(shù)隨不同F(xiàn)r數(shù)下航行的變化曲線,通過對(duì)比國(guó)內(nèi)外不同文獻(xiàn)的興波理論計(jì)算值,計(jì)算結(jié)果雖有些不同,但隨Fr數(shù)的變化較為一致,說明了本文程序計(jì)算興波阻力的可靠性。

        圖8 Wigley船型興波阻力系數(shù)Fig.8 Wave-making resistance coefficient of Wigley ship form

        圖9 KCS船型興波阻力系數(shù)Fig.9 Wave-making resistance coefficient of KCS ship form

        3 近似模型

        3.1 SLE方法

        圖10 SLE方法樣本生成過程Fig.10 Sample generation process of SLE method

        圖11 試驗(yàn)樣本空間填充性能對(duì)比Fig.11 Comparison of filling performance of test sample space

        3.2 BP神經(jīng)網(wǎng)絡(luò)

        圖12 為BP 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu),包含輸入層、隱含層和輸出層。輸入層對(duì)當(dāng)前輸入的數(shù)據(jù)進(jìn)行處理,隱含層把輸入層的數(shù)據(jù)處理好后傳向輸出層,輸出層給出當(dāng)此正向傳播的輸出值和訓(xùn)練誤差。步驟如下:

        圖12 BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)Fig.12 BP neural network structure

        步驟1:初始化BP 網(wǎng)絡(luò)參數(shù),包括輸入層、隱含層和輸出層神經(jīng)元個(gè)數(shù),初始化權(quán)值w和閾值b,設(shè)置最大訓(xùn)練次數(shù)、最小精度和學(xué)習(xí)率。

        步驟2:正向過程計(jì)算,通過輸入試驗(yàn)樣本值,計(jì)算各層輸出和訓(xùn)練損失函數(shù),如果滿足終止訓(xùn)練的要求,則近似模型訓(xùn)練完成,此時(shí)神經(jīng)網(wǎng)絡(luò)達(dá)到預(yù)期目標(biāo)。

        步驟3:誤差計(jì)算,判斷是否滿足收斂條件。

        步驟4:反向過程計(jì)算,如果不滿足終止條件,將誤差反向傳輸并更新權(quán)值和偏置,繼續(xù)訓(xùn)練。

        3.3 近似模型程序設(shè)計(jì)

        如圖13 所示,本文通過SLE 試驗(yàn)設(shè)計(jì)方法和BP 神經(jīng)網(wǎng)絡(luò)相結(jié)合,建立的近似模型可代替Dawson方法預(yù)報(bào)船體興波阻力。

        圖13 近似模型流程圖Fig.13 Flow chart of approximate model

        4 KCS船型計(jì)算分析

        4.1 船型參數(shù)化設(shè)計(jì)

        采用ITM 和FFD 方法對(duì)KCS 前半部分變形,設(shè)計(jì)變量如圖14 所示,α1和α2控制前半船在船長(zhǎng)和船寬方向的變形,dx、dy、dz為控制框架里部分控制點(diǎn)沿著三個(gè)方向的位移,表1顯示了約束范圍均轉(zhuǎn)換為相對(duì)于船長(zhǎng)的無因次量。采用SLE 方法生成400 個(gè)含有5 個(gè)特征變量的樣本船型,取Fr=0.26 并計(jì)算興波阻力值,從圖15中可以看出,較多樣本船在興波阻力方面表現(xiàn)優(yōu)于母型船。

        圖14 KCS船體變形示意圖Fig.14 Diagram of KCS hull deformation

        表1 變形方法參數(shù)設(shè)置Tab.1 Parameter setting of deformation method

        圖15 船型樣本集Fig.15 Sample set of ship type

        4.2 網(wǎng)絡(luò)參數(shù)對(duì)性能的影響分析

        為了分析不同網(wǎng)絡(luò)參數(shù)對(duì)近似模型預(yù)測(cè)精度的影響,本文選擇均方誤差Emse和相關(guān)系數(shù)R2為評(píng)價(jià)指標(biāo):

        式中,yi、y?i分別為樣本的理想值和預(yù)測(cè)值,為理想值的平均值,n為樣本數(shù)。

        設(shè)網(wǎng)絡(luò)模型訓(xùn)練樣本數(shù)為400,最大訓(xùn)練次數(shù)為50 000次,訓(xùn)練精度為0.000 01。圖16顯示的是訓(xùn)練模型評(píng)價(jià)指標(biāo)隨學(xué)習(xí)率的變化,可以看出相關(guān)系數(shù)越大,均方誤差越小。如圖17所示,隱含層的神經(jīng)元節(jié)點(diǎn)數(shù)影響近似模型的預(yù)測(cè)精度和訓(xùn)練產(chǎn)生的時(shí)間成本。根據(jù)數(shù)據(jù)對(duì)比結(jié)果可看出,本文學(xué)習(xí)率為0.99時(shí),網(wǎng)絡(luò)模型的訓(xùn)練性能較好,當(dāng)隱含層節(jié)點(diǎn)數(shù)為6層時(shí),相關(guān)系數(shù)較大,時(shí)間成本相對(duì)較低。

        圖16 學(xué)習(xí)率對(duì)性能評(píng)價(jià)指標(biāo)的影響Fig.16 Influence of learning rate on performance indicators

        圖17 隱含層節(jié)點(diǎn)數(shù)對(duì)性能評(píng)價(jià)指標(biāo)的影響Fig.17 Influence of number of hidden layer nodes on performance

        4.3 訓(xùn)練實(shí)驗(yàn)及結(jié)果分析

        隱含層節(jié)點(diǎn)數(shù)和學(xué)習(xí)率參數(shù)確定之后,按照表2 網(wǎng)絡(luò)參數(shù)開始訓(xùn)練近似模型。為了驗(yàn)證近似模型預(yù)測(cè)的準(zhǔn)確性,在設(shè)計(jì)空間內(nèi),隨機(jī)生成60 個(gè)測(cè)試樣本船型并計(jì)算興波阻力值。近似模型總誤差隨訓(xùn)練次數(shù)的變化曲線如圖18所示,可以看出總誤差E逐漸收斂為一個(gè)極小值。然后,使用訓(xùn)練之后的近似模型對(duì)測(cè)試樣本船型集的60組樣本進(jìn)行預(yù)測(cè),圖19為測(cè)試樣本預(yù)測(cè)值的相對(duì)誤差統(tǒng)計(jì)圖。

        表2 近似模型結(jié)構(gòu)參數(shù)Tab.2 Structure parameters of approximate model

        圖18 訓(xùn)練次數(shù)與總誤差的關(guān)系Fig.18 Relationship between training times and total error

        圖19 測(cè)試樣本的相對(duì)誤差變化Fig.19 Relative error variation of the test sample

        利用訓(xùn)練好的模型對(duì)60 組測(cè)試樣本Cws進(jìn)行預(yù)測(cè),平均相對(duì)誤差為0.282%,相關(guān)系數(shù)R2為0.9998,均方誤差為3.14E-04,運(yùn)行時(shí)間為2 ms。其中,部分測(cè)試樣本的網(wǎng)絡(luò)輸出如表3 所示,表中輸入數(shù)據(jù)為歸一化后的值,通過分析表明,本文構(gòu)建的智能預(yù)報(bào)方法預(yù)測(cè)樣本船型的興波阻力的精度高且速度快。

        表3 測(cè)試樣本的網(wǎng)絡(luò)輸出Tab.3 Network output of test samples

        5 結(jié) 論

        本文建立了基于樣本船的興波阻力智能預(yù)報(bào)方法,結(jié)合曲面變形技術(shù)和興波阻力預(yù)報(bào)技術(shù),實(shí)現(xiàn)了在設(shè)計(jì)空間內(nèi)船體興波阻力的快速預(yù)報(bào),并對(duì)近似模型的相關(guān)參數(shù)和訓(xùn)練展開分析,得出如下結(jié)論:

        (1)提出了基于曲面控制網(wǎng)格的曲面變形方法,編寫了ITM 和FFD 方法的曲面變形程序,對(duì)KCS船型進(jìn)行參數(shù)化設(shè)計(jì)并生成樣本船型集。

        (2)在曲面變形技術(shù)和興波阻力預(yù)報(bào)技術(shù)基礎(chǔ)上,通過SLE試驗(yàn)設(shè)計(jì)方法和BP神經(jīng)網(wǎng)絡(luò)相結(jié)合,建立的近似模型可快速預(yù)報(bào)船體興波阻力,耗時(shí)成本為毫秒量級(jí)。

        (3)通過網(wǎng)絡(luò)模型參數(shù)分析和訓(xùn)練,測(cè)試樣本預(yù)測(cè)值相關(guān)系數(shù)為0.9998,均方誤差為3.14E-04,本文近似模型的預(yù)測(cè)精度較高。

        猜你喜歡
        興波船型船體
        船體行駛過程中的壓力監(jiān)測(cè)方法
        船舶興波運(yùn)動(dòng)的非線性現(xiàn)象算法研究與仿真
        中國(guó)古代著名的海船船型
        ——福船
        焊接殘余應(yīng)力對(duì)船體結(jié)構(gòu)疲勞強(qiáng)度的影響分析
        焊接(2015年9期)2015-07-18 11:03:51
        中國(guó)古代著名的海船船型—廣船
        赴美軍“仁慈”號(hào)醫(yī)院船駐船體會(huì)
        Anti-ageing effects of a new Dimethylaminoethanol-based formulation on DGalactose induced skin ageing model of rat
        水下爆炸氣泡作用下船體總縱強(qiáng)度估算方法
        內(nèi)傾船型的非線性橫搖仿真分析
        影響潛航體水面興波因素分析
        中文字幕人妻被公上司喝醉| 日日噜噜噜夜夜狠狠久久蜜桃 | 久久久久亚洲AV片无码乐播| 日韩一区二区av伦理| 九九影院理论片私人影院| 亚洲男人的天堂在线播放| 精品无人区无码乱码大片国产| 国产三级av在线精品| 人妻精品久久久久中文字幕| 欧美mv日韩mv国产网站| 国产成人亚洲综合小说区| av免费一区二区久久| 日韩人妻无码一区二区三区久久| 国内少妇偷人精品视频免费| 极品美女尤物嫩模啪啪| 东北熟妇露脸25分钟| 国产亚洲真人做受在线观看| 亚洲永久精品ww47永久入口| 亚洲综合有码中文字幕| 刺激一区仑乱| 日韩精品无码一区二区三区视频| 国产西西裸体一级黄色大片| 国产成人高清在线观看视频| 久久亚洲私人国产精品va| 日韩欧美在线播放视频| 人妻少妇中文字幕专区| 一色桃子中文字幕人妻熟女作品| 护士奶头又白又大又好摸视频| 男女在线免费视频网站| 中文字幕乱码亚洲无限码| 中文成人无码精品久久久不卡 | 黄片视频免费在线播放观看 | 无码不卡高清毛片免费| 精品女同一区二区三区免费播放 | 丝袜国产高跟亚洲精品91| 久久国产女同一区二区| 高清毛茸茸的中国少妇| 国产精品久久久久久久久鸭 | h视频在线观看视频在线| 久久久亚洲欧洲日产国码aⅴ| 亚洲精品无码久久毛片|