• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    碼頭系泊系統(tǒng)快速分析軟件開(kāi)發(fā)

    2022-12-22 03:58:46時(shí)永鵬閆鳳超
    船舶 2022年6期
    關(guān)鍵詞:纜繩系泊拉力

    時(shí)永鵬 閆鳳超

    (中國(guó)船舶及海洋工程設(shè)計(jì)研究院 上海 200001)

    0 引 言

    船舶系泊事故時(shí)有發(fā)生,系泊事故可能會(huì)造成結(jié)構(gòu)損壞、纜繩斷裂以及船舶狀態(tài)的突變等。尤其發(fā)生纜繩斷裂時(shí),纜繩內(nèi)部能量突然釋放會(huì)造成纜繩快速回彈,極易造成人員傷亡[1-2]。因此系泊安全受到船東越來(lái)越多的關(guān)注,采取的措施包括標(biāo)識(shí)纜繩反彈區(qū)、設(shè)置護(hù)欄以及進(jìn)行系泊分析等,其中系泊分析是在詳細(xì)設(shè)計(jì)階段和系泊實(shí)施前可以完成的主要的系泊風(fēng)險(xiǎn)評(píng)估方式[3]?;谙挡从?jì)算、分析,可以對(duì)系泊設(shè)計(jì)進(jìn)行評(píng)估、優(yōu)化,也可以對(duì)船舶可承受的極限環(huán)境載荷進(jìn)行預(yù)先判斷。

    規(guī)范計(jì)算是實(shí)船設(shè)計(jì)中最常用的系泊計(jì)算方法,常用的規(guī)范包括IACS REC.10[4-5]和各船級(jí)社規(guī)范。但規(guī)范計(jì)算方法過(guò)于簡(jiǎn)化,許多重要的系泊因素(包括船舶類型、碼頭條件、系泊布置以及纜繩材料等)都無(wú)法考慮在內(nèi),計(jì)算精度較低。

    石油公司組織OCIMF 的出版物“MOORING EQUIPMENT GUIDELINES”(下文稱為MEG,最新版為MEG4)中為液貨船提供了更詳盡的系泊分析指導(dǎo)和相關(guān)數(shù)據(jù)[6-7],并被很多系泊分析方法借鑒,然而其并沒(méi)有提供直接可用的系泊分析工具。

    使用成熟的商業(yè)軟件是目前最可靠的系泊分析方法,但不同的軟件有不同的特點(diǎn)或適用性,如由法國(guó)船級(jí)社開(kāi)發(fā)的Ariane 系泊分析軟件,其計(jì)算方法可考慮風(fēng)、浪、流等復(fù)雜的環(huán)境條件,適用于單點(diǎn)、多點(diǎn)和碼頭系泊等不同的系泊場(chǎng)景。然而該軟件有較高的使用門(mén)檻,需要輸入繁復(fù)的數(shù)據(jù)參數(shù),易出錯(cuò)且計(jì)算過(guò)程也很耗時(shí)。如要完成一次多參數(shù)的方案優(yōu)化,一般需要一周或更長(zhǎng)時(shí)間,并且其計(jì)算原理也更適用于考慮波浪載荷的海工產(chǎn)品。

    本文借鑒已有的系泊分析方法和軟件,對(duì)碼頭系泊系統(tǒng)的組成及特性進(jìn)行分析,實(shí)現(xiàn)1 套具有更低用戶使用門(mén)檻、更高計(jì)算效率,且兼具準(zhǔn)確性的碼頭系泊分析軟件的開(kāi)發(fā)。

    1 碼頭系泊系統(tǒng)的組成

    本文以碼頭系泊作為研究對(duì)象,系泊系統(tǒng)的形成過(guò)程是船舶、碼頭與環(huán)境的相互作用過(guò)程,因此本文將碼頭系泊系統(tǒng)分為環(huán)境系統(tǒng)、船舶系統(tǒng)和碼頭系統(tǒng)3 部分。為提高準(zhǔn)確性,計(jì)算中需將各個(gè)系統(tǒng)中的重要因素納入計(jì)算模型,但也需要結(jié)合開(kāi)發(fā)目標(biāo)進(jìn)行合理簡(jiǎn)化。

    1.1 環(huán)境系統(tǒng)

    船舶可能受到的載荷種類非常豐富,包括常見(jiàn)的風(fēng)、浪、流,以及潮汐、涌浪和船行波等,而船舶在不同系泊場(chǎng)景下需側(cè)重考慮的環(huán)境因素也不同。

    碼頭系泊一般位于遮蔽水域,環(huán)境因素較簡(jiǎn)單,其中規(guī)范和OCIMF 推薦的碼頭系泊環(huán)境因素如表1所示。

    表1中,風(fēng)速為高于地面10 m 處的30 s 內(nèi)的平均風(fēng)速;30 s 被視為系泊系統(tǒng)的平均響應(yīng)速度。環(huán)境載荷從船尾指向船首時(shí),其角度定為0°。

    結(jié)合表1及實(shí)船設(shè)計(jì)經(jīng)驗(yàn),環(huán)境系統(tǒng)中僅考慮風(fēng)和流,可以滿足大部分船舶的碼頭系泊分析需求。此外,僅考慮定常的風(fēng)、流載荷,簡(jiǎn)化環(huán)境系統(tǒng),可提高系泊分析效率。

    1.2 船舶系統(tǒng)

    船舶系統(tǒng)主要包括船舶參數(shù)和船舶系泊,前者直接影響環(huán)境載荷的大??;后者影響系泊力的內(nèi)部分配。

    船舶系統(tǒng)的詳細(xì)組成如下頁(yè)圖1所示,本軟件將盡可能準(zhǔn)確地對(duì)船舶系統(tǒng)進(jìn)行描述。

    圖1 船舶系統(tǒng)的主要組成

    1.3 碼頭系統(tǒng)

    碼頭系統(tǒng)包括碼頭高度、碼頭纜樁布置、碼頭水深、碼頭碰墊的布置和物理性能等(見(jiàn)圖2),其中碼頭布置影響系泊纜繩的實(shí)際長(zhǎng)度和角度、碼頭水深影響環(huán)境載荷的大小、碰墊直接影響系泊力及船舶位移。本軟件將盡可能真實(shí)地表達(dá)碼頭系統(tǒng)信息。

    圖2 碼頭系統(tǒng)示意圖

    2 系泊系統(tǒng)的數(shù)學(xué)模型

    在定常環(huán)境載荷下,系泊力與環(huán)境力將達(dá)到平衡,系泊系統(tǒng)達(dá)到穩(wěn)定狀態(tài);在穩(wěn)定的載荷作用下,系泊系統(tǒng)可視為一個(gè)定常系統(tǒng)。本文將定常定義為一個(gè)相對(duì)的時(shí)間長(zhǎng)度,是船舶在一定的環(huán)境載荷下,系泊系統(tǒng)與之達(dá)到平衡所需的時(shí)間。參照IACS REC.10 與MEG4 中的描述,該時(shí)間長(zhǎng)度定為30 s;30 s 是一個(gè)足夠短的時(shí)間段,允許我們將變化的外界環(huán)境離散為不同大小的定常環(huán)境的疊加,從而將碼頭系泊分析轉(zhuǎn)變?yōu)槎ǔ7治觥?/p>

    當(dāng)碼頭系泊系統(tǒng)中出現(xiàn)涌浪、波浪和船行波等不穩(wěn)定因素時(shí),系泊系統(tǒng)無(wú)法在短時(shí)間內(nèi)達(dá)到穩(wěn)定狀態(tài),一般需要針對(duì)特定的計(jì)算案例進(jìn)行CFD 水動(dòng)力分析并進(jìn)行時(shí)域仿真計(jì)算,本文暫不考慮這種情況。

    2.1 坐標(biāo)系

    本文中坐標(biāo)系的定義如圖3所示。假設(shè)船舶右舷靠泊碼頭,x軸正向?yàn)榭v向向艏、y軸正向?yàn)闄M向指向左舷、z軸正向?yàn)榇瓜蛳蛏希睾徒嵌榷家阅鏁r(shí)針為正。

    圖3 坐標(biāo)系定義

    本文中,全局坐標(biāo)為基于碼頭的固定坐標(biāo)系,局部坐標(biāo)為基于船體的運(yùn)動(dòng)坐標(biāo)系。在初始狀態(tài)下,全局坐標(biāo)與局部坐標(biāo)有相同且以上文描述為基準(zhǔn)的x、y、z的坐標(biāo)軸方向;船體運(yùn)動(dòng)即局部坐標(biāo)系在全局坐標(biāo)系中的運(yùn)動(dòng)。

    2.2 平衡方程

    將船舶視為剛體,船舶運(yùn)動(dòng)存在6 向自由度。系泊系統(tǒng)達(dá)到平衡狀態(tài)時(shí),船體所受的合力為0。將船舶受力分解到6 個(gè)方向,分別為Fx、Fy、Fz、Mx、My及Mz,可形成6 個(gè)平衡方程。基于后文第2.6節(jié)論述,對(duì)船舶運(yùn)動(dòng)進(jìn)行簡(jiǎn)化,簡(jiǎn)化后僅需3 個(gè)力平衡方程,即Fx=0、Fy=0 及Mz=0。

    船體受到的力可以歸類為纜繩力、碰墊力和環(huán)境力。由于每根纜繩長(zhǎng)度、角度及位置不同,需單獨(dú)計(jì)算每根纜繩的拉力及其在3 個(gè)方向的分力;碰墊力同樣需要根據(jù)位置分別計(jì)算。將這些力疊加,形成平衡方程如下:

    式中:n為纜繩的數(shù)量;m為碰墊的數(shù)量;FxiLine、FyiLine、MziLine分別為第i根纜繩在x、y以及z向的分力和力矩;FjFender、MzjFender代表第j個(gè)碰墊在y向和z向的力與力矩;同理,F(xiàn)xWind、FyWind、MzWind和FxCurr、FyCurr、MzCurr代表風(fēng)和流在3 個(gè)方向的力和力矩。下文將對(duì)上式中的各個(gè)力分別進(jìn)行分析計(jì)算。

    2.3 環(huán)境力

    環(huán)境力是環(huán)境與船體的作用結(jié)果。本文中環(huán)境載荷的計(jì)算使用MEG4 中的環(huán)境系數(shù)法。MEG4 基于大量的模型實(shí)驗(yàn)和計(jì)算流體力學(xué)軟件分析、得到無(wú)因次化的風(fēng)載和流載系數(shù);該系數(shù)是環(huán)境載荷角度、水深吃水比、船型特征、船舶載況等的函數(shù);環(huán)境系數(shù)再結(jié)合實(shí)際計(jì)算的風(fēng)和流的速度、船舶主尺度、受風(fēng)面積等參數(shù),可計(jì)算得到環(huán)境載荷的 大小。

    如以x方向的流載荷FXc的計(jì)算為例,其計(jì)算公式如下:

    式中:ρc為水的密度,kg/m3;Vc為水流速度,m/s;LBP為船長(zhǎng),m;T為計(jì)算吃水,m。CXc為MEG4中提供的無(wú)因次化的流載系數(shù),計(jì)算方法如下:

    MEG4 中,函數(shù)fc通過(guò)圖表的方式呈現(xiàn)。本軟件通過(guò)讀取圖表數(shù)據(jù)、使用數(shù)值擬合和插值方法實(shí)現(xiàn)該函數(shù)的計(jì)算,可適應(yīng)任意合理的用戶輸入。其他環(huán)境系數(shù)和環(huán)境載荷的計(jì)算同理,更多的計(jì)算方法和公式請(qǐng)參照MEG4 中的描述。

    2.4 纜繩力

    纜繩的位置、長(zhǎng)度、角度及材質(zhì)等因素都會(huì)影響其拉力大小。本軟件中建立纜繩真實(shí)的三維空間模型,并定義其物理特性。實(shí)際系泊狀態(tài)下,纜繩可能存在多個(gè)拐點(diǎn),如典型的纜繩軌跡為“絞車-羊角滾輪-導(dǎo)纜孔-碼頭纜樁”,該4 點(diǎn)可確定纜繩的三維形狀。計(jì)算模型中,纜繩的三維模型也通過(guò)一系列控制點(diǎn)描述,如:P0 表示纜繩在碼頭的固定點(diǎn)(碼頭纜樁);P1 表示纜繩在船體的出口點(diǎn)(導(dǎo)纜孔);P2、P3 代表纜繩在船上的控制點(diǎn)(羊角、帶纜樁、絞車等)。可以發(fā)現(xiàn):P0 在全局坐標(biāo)系中的坐標(biāo)不會(huì)變化,而P1、P2、P3 等在船體局部坐標(biāo)系中的坐標(biāo)不會(huì)變化;船體運(yùn)動(dòng)會(huì)導(dǎo)致P0 與P1 相對(duì)位置的變化,而船體上的點(diǎn)(P2、P3 等)的相對(duì)位置不會(huì)發(fā)生變化。

    船體位移引起纜繩長(zhǎng)度變化,在船體坐標(biāo)系發(fā)生角位移θ和線位移(?x,?y)后,纜繩控制點(diǎn)P1(x,y)的新全局坐標(biāo)為(x′,y′),兩者有以下關(guān)系:

    根據(jù)P1 的新坐標(biāo),可計(jì)算纜繩的伸長(zhǎng)量?l及拉力FLine。伸長(zhǎng)量與拉力的函數(shù)關(guān)系通過(guò)用戶自定義的拉力-變形率函數(shù)體現(xiàn),如下列公式中的函數(shù)f。

    式中:MBL為纜繩的最小破斷負(fù)荷;el為纜繩伸長(zhǎng)率,即伸長(zhǎng)量Δl與纜繩無(wú)拉力狀態(tài)下的原始長(zhǎng)度lLine的比值。由于纜繩僅能提供單向力,拉力-變形率函數(shù)在0 點(diǎn)一階不連續(xù)。

    此外,纜繩的初始拉力與尾索配置也會(huì)影響纜繩最終的拉力大小。

    (1)初始拉力的影響

    實(shí)船系泊時(shí),會(huì)在系泊初始狀態(tài)下,通過(guò)絞車或其他方式給予纜繩初始拉力。初始拉力可減小船舶在環(huán)境載荷下的運(yùn)動(dòng),也可能降低纜繩的最大拉力,是一個(gè)重要的系泊參數(shù)。對(duì)應(yīng)到計(jì)算模型中,船舶的初始狀態(tài)應(yīng)是環(huán)境載荷為0、纜繩的初始拉力與碰墊力達(dá)到平衡的狀態(tài)。系泊計(jì)算的第一步必須基于用戶給出的初始拉力,確定系泊的初始狀態(tài),包括船舶的初始位移及纜繩的初始變形。

    (2)尾索的影響

    當(dāng)纜繩的剛度太大時(shí),一些船舶可能會(huì)設(shè)置尾索以優(yōu)化纜繩受力。由于主纜繩和尾索有不同的物理特性(即不同的拉力-變形率函數(shù)),在相同的軸向拉力下,兩者有不同的變形率,可通過(guò)下列公式計(jì)算得到平衡狀態(tài)下兩者各自的變形。

    式中:f代表纜繩的拉力-變形率函數(shù);下標(biāo)Line代表主纜繩,Tail 代表尾索。

    2.5 碰墊力

    碰墊在被壓縮(與船體的接觸點(diǎn)產(chǎn)生y向負(fù)位移)時(shí),會(huì)對(duì)船體產(chǎn)生反彈力。對(duì)于碼頭系泊系統(tǒng),風(fēng)和流為離岸方向時(shí),纜繩受力較大;風(fēng)和流為向岸方向時(shí),主要靠碰墊抵抗環(huán)境載荷。

    碰墊反彈力與變形率也存在類似于纜繩變形的函數(shù)關(guān)系,本文不再詳述;且該函數(shù)同樣在0 點(diǎn)一階不連續(xù)。碰墊僅會(huì)產(chǎn)生y向分力FFender,且始終為正;該力乘以碰墊的x向坐標(biāo)可得到碰墊產(chǎn)生的z向力矩MzFender。

    2.6 船舶運(yùn)動(dòng)的簡(jiǎn)化

    船舶運(yùn)動(dòng)引起導(dǎo)纜孔位置變化,從而導(dǎo)致纜繩拉力變化。導(dǎo)纜孔作為船體上的一個(gè)點(diǎn),具有三向運(yùn)動(dòng)自由度??芍庇^看出,導(dǎo)纜孔的x向和y向位移對(duì)于纜繩變形有顯著的影響。本節(jié)通過(guò)分析導(dǎo)纜孔的z向位移引起的纜繩長(zhǎng)度變化,反向推導(dǎo)船舶運(yùn)動(dòng)的影響。

    船舶的升沉、橫傾和縱傾會(huì)導(dǎo)致導(dǎo)纜孔的z向位移;同時(shí),該三向運(yùn)動(dòng)也主要是使導(dǎo)纜孔產(chǎn)生z向位移。

    (1)船舶的升沉

    船舶所受的纜繩力、碰墊力和環(huán)境力中,只有纜繩力會(huì)產(chǎn)生顯著的z向分力,造成船舶升沉。

    MEG4 中推薦的纜繩水平夾角為15°,在不利的輕載載況下,建議不超過(guò)25°。以某5 萬(wàn)載重噸油輪為例,其配置纜繩16 根,其中倒纜4 根。假設(shè)纜繩的水平夾角為不利的25°,在MEG4 基準(zhǔn)環(huán)境條件下,其纜繩拉力極值約為27 t,均值約15 t,16 根纜繩產(chǎn)生z向分力約為100 t。

    該船輕載吃水的水線面積約為4 500 m2,在100 t 的z向壓力下,升沉約20 mm;在25°水平夾角下,會(huì)導(dǎo)致纜繩約8.5 mm 的長(zhǎng)度變化。纜繩長(zhǎng)度平均約為30 m,其長(zhǎng)度變化率約為2.8×10-4。對(duì)于MBL約為60 t 的常規(guī)纖維纜,其產(chǎn)生的拉力變化不足0.1 t,故船舶的升沉運(yùn)動(dòng)可被合理忽略。

    (2)船舶的橫傾與縱傾

    風(fēng)和流為90°時(shí),船舶所受的橫向(y向)環(huán)境力最大,作為考慮橫傾影響的極端工況。

    船舶輕載時(shí),風(fēng)載是主要環(huán)境載荷,流載約是其5%,且流載的橫傾力臂很小。纜繩的y向拉力與風(fēng)載基本相等,風(fēng)載的力臂從上建到水面分布,可以認(rèn)為風(fēng)載與纜繩產(chǎn)生的x向力矩相互抵消,產(chǎn)生的橫傾力矩可忽略,如下頁(yè)圖4所示。

    圖4 不同吃水下的船舶受力

    滿載吃水時(shí),流載增大到風(fēng)載的75%左右,且流載力臂增大、風(fēng)載力臂減小,可認(rèn)為兩者產(chǎn)生的x向力矩大致抵消,橫傾力矩主要由纜繩產(chǎn)生。假設(shè)纜繩的水平夾角為10°,同樣以上述油輪為例,纜繩平均拉力為15 t,除去倒纜的纜繩數(shù)量為12根,導(dǎo)纜孔距水線約6 m,纜繩產(chǎn)生的x向力矩約為1×103量級(jí)。該船滿載時(shí),橫傾1°所需力矩約為8×103量級(jí),計(jì)算可得橫傾引起的導(dǎo)纜孔z向位移約35 mm、纜繩長(zhǎng)度變化約6 mm,由此表明橫傾導(dǎo)致的纜繩拉力變化可忽略。

    艏艉纜繩對(duì)船舶產(chǎn)生的縱傾力矩本身會(huì)相互抵消,同時(shí)船舶縱傾所需力矩顯著大于橫傾,類似方法分析得到,船舶縱傾的影響可忽略。

    上述分析表明:導(dǎo)致導(dǎo)纜孔產(chǎn)生z向位移的船舶升沉、橫傾和縱傾三向位移可以被合理忽略。值得注意的是,本節(jié)所述的船舶運(yùn)動(dòng)簡(jiǎn)化是在一定數(shù)量級(jí)上的分析,并非嚴(yán)謹(jǐn)推導(dǎo)。

    3 軟件開(kāi)發(fā)及用戶界面

    本軟件使用Visual Studio 環(huán)境和C#語(yǔ)言,完成系泊系統(tǒng)分析的算法實(shí)現(xiàn)和用戶界面開(kāi)發(fā)。算法實(shí)現(xiàn)保證計(jì)算效率和準(zhǔn)確性,用戶界面則應(yīng)保證以較低的門(mén)檻引導(dǎo)用戶快速、準(zhǔn)確地輸入必要的數(shù)據(jù),并直觀地輸出計(jì)算結(jié)果。

    按照上文所述的系泊系統(tǒng)劃分方式,軟件引導(dǎo)用戶依次輸入船舶、環(huán)境、碼頭、纜繩及系泊布置信息。

    3.1 基本信息輸入

    (1)輸入船型信息和主尺度信息,見(jiàn)圖5。

    圖5 用戶輸入-船舶參數(shù)

    (2)輸入環(huán)境載荷見(jiàn)下頁(yè)圖6,用戶可以批量設(shè)置需要計(jì)算的風(fēng)和流的方向、以及在不同方向上的速度,軟件會(huì)自動(dòng)生成所有環(huán)境組合。

    圖6 用戶輸入-環(huán)境參數(shù)

    (3)輸入碼頭信息(見(jiàn)圖7),含水深、碼頭高度和碰墊的物理特性。軟件通過(guò)收集常用的碰墊數(shù)據(jù)信息,對(duì)不同尺寸碰墊的變形-彈力曲線進(jìn)行回歸,得到無(wú)因次曲線;用戶輸入典型變形下的碰墊彈力值,可自動(dòng)生成碰墊的變形-彈力曲線。此外,軟件也允許用戶輸入自定義的碰墊變形-彈力曲線。

    圖7 用戶輸入-碼頭參數(shù)

    (4)輸入纜繩信息(見(jiàn)圖8),包括主纜繩與尾索的物理特性、纜繩MBL以及初始拉力。軟件內(nèi)置多種類型纜繩的物理特性數(shù)據(jù)庫(kù),用戶可快速選擇,也可通過(guò)輸入?yún)?shù)、自由定義纜繩物理特性。軟件會(huì)根據(jù)用戶輸入,繪制纜繩的拉力-變形率曲線圖。

    圖8 用戶輸入-纜繩參數(shù)

    此外,軟件可將用戶的所有輸入另存為獨(dú)立的系泊計(jì)算模型文件,方便保存與傳輸。

    3.2 系泊布置信息

    系泊系統(tǒng)的布置信息包括所有纜繩的三維控制點(diǎn)信息、碼頭碰墊和碼頭纜樁的布置信息,數(shù)據(jù)量較大,且在方案優(yōu)化時(shí)需要頻繁修改。

    在常見(jiàn)的系泊軟件中,需要用戶讀取系泊圖紙并在軟件中輸入所有控制點(diǎn)的三維坐標(biāo),效率很低,尤其發(fā)生頻繁的方案修改時(shí),需要反復(fù)進(jìn)行“圖紙修改-讀取坐標(biāo)-修改軟件輸入”的流程,極易發(fā)生數(shù)據(jù)不同步的錯(cuò)誤。

    通過(guò)定義與CAD 繪圖軟件的標(biāo)準(zhǔn)接口,可將存儲(chǔ)在CAD 圖紙中的系泊布置信息直接讀取至本軟件,節(jié)省大量用戶時(shí)間;系泊圖紙發(fā)生修改時(shí),其信息也會(huì)同步到本軟件,保證數(shù)據(jù)一致性,并非常高效、準(zhǔn)確地完成系泊設(shè)計(jì)的迭代優(yōu)化。如下頁(yè)圖9展示的用于本軟件輸入的標(biāo)準(zhǔn)化系泊布置圖紙,其形式與用于送審的系泊布置圖紙基本一致。

    圖9 用戶輸入-系泊布置信息

    3.3 軟件輸出

    軟件會(huì)直接輸出所有環(huán)境和吃水組合工況的纜繩拉力極值,將其作為系泊設(shè)計(jì)的“Working Load Limit”,并根據(jù)纜繩材質(zhì)計(jì)算“Ship design MBL” (“Working Load Limit”和“Ship design MBL”基于MEG4 中的定義)和纜繩最小破斷負(fù)荷;同時(shí)還會(huì)輸出計(jì)算的中間數(shù)據(jù),用戶可進(jìn)一步進(jìn)行數(shù)據(jù)分析,或挖掘規(guī)律進(jìn)行針對(duì)性的系泊優(yōu)化等,中間數(shù)據(jù)包括所有環(huán)境和吃水組合下,所有纜繩的受力大小、船體的位移和作用在船體上的環(huán)境載荷等,如圖10所示。

    圖10 軟件輸出-部分輸出數(shù)據(jù)

    軟件計(jì)算結(jié)果對(duì)于任何的輸入修改都會(huì)即時(shí)響應(yīng),即本軟件可實(shí)現(xiàn)從“設(shè)計(jì)輸入 ”到“設(shè)計(jì)結(jié)果”的即時(shí)數(shù)據(jù)關(guān)聯(lián)。

    3.4 軟件計(jì)算的準(zhǔn)確性校核

    本軟件的計(jì)算結(jié)果經(jīng)過(guò)多型船和多種系泊方案的計(jì)算驗(yàn)證,其計(jì)算結(jié)果與計(jì)算變量之間的影響關(guān)系符合一般的系泊規(guī)律,且經(jīng)過(guò)與OCIMF 計(jì)算方法和Ariane 軟件的計(jì)算結(jié)果進(jìn)行對(duì)比,其計(jì)算精度可以滿足使用要求。表2為多種計(jì)算方法的計(jì)算結(jié)果對(duì)比(纜繩拉力極值)。

    表2 多種計(jì)算方法的計(jì)算結(jié)果對(duì)比(纜繩拉力極值)

    表2中,計(jì)算所用環(huán)境條件為MEG4 中的基準(zhǔn)環(huán)境條件。OCIMF 計(jì)算方法中,纜繩受力采用等分法并乘以相應(yīng)系數(shù)得到,詳細(xì)可參考MEG 提供的算法示例。

    4 結(jié) 語(yǔ)

    本文通過(guò)對(duì)碼頭系泊系統(tǒng)的分析和合理簡(jiǎn)化,建立算法和用戶界面,完成了一套簡(jiǎn)單易用、快速且準(zhǔn)確的系泊分析軟件的開(kāi)發(fā)。本軟件的使用可以達(dá)到以下目的:

    (1)軟件較低的使用門(mén)檻可以適用于更多的場(chǎng)景和用戶;

    (2)計(jì)算的簡(jiǎn)便、快速性能夠讓用戶對(duì)于系泊優(yōu)化從被動(dòng)轉(zhuǎn)變?yōu)橹鲃?dòng);

    (3)軟件的目的是建立數(shù)據(jù)間的關(guān)系,系泊系統(tǒng)只是船舶系統(tǒng)的一部分,計(jì)算的即時(shí)性讓用戶能夠直觀看到“設(shè)計(jì)”與“結(jié)果”的數(shù)據(jù)關(guān)聯(lián),保證船舶設(shè)計(jì)的連續(xù)性和可靠性;

    (4)自主可控的軟件有非常好的擴(kuò)展性,可以直接嵌入到船舶設(shè)計(jì)大系統(tǒng)中,實(shí)現(xiàn)更大范圍的系統(tǒng)優(yōu)化。

    在本文軟件開(kāi)發(fā)的基礎(chǔ)上,下一步將在軟件中集成更多系泊設(shè)計(jì)的知識(shí)和經(jīng)驗(yàn),逐步實(shí)現(xiàn)從“被動(dòng)”給出計(jì)算結(jié)果到“主動(dòng)”提供設(shè)計(jì)優(yōu)化建議的改變。

    猜你喜歡
    纜繩系泊拉力
    基于CFD方法的系泊狀態(tài)FSO舵載荷計(jì)算
    高韌性海洋用三級(jí)系泊鏈鋼的研制
    山東冶金(2022年2期)2022-08-08 01:50:50
    MEG4 規(guī)范對(duì)導(dǎo)纜孔的新要求及 相應(yīng)的設(shè)計(jì)改進(jìn)
    基于CFD的螺旋槳拉力確定方法
    船舶靠泊過(guò)程中系纜作業(yè)仿真
    自不量力
    深水工程船舶纜繩測(cè)力裝置設(shè)計(jì)與應(yīng)用
    跟蹤導(dǎo)練(三)(3)
    水文測(cè)船調(diào)速電動(dòng)纜繩收放機(jī)設(shè)計(jì)與應(yīng)用
    等動(dòng)拉力和重力拉力蝶泳劃臂動(dòng)作生物力學(xué)分析
    在线观看免费视频日本深夜| 黄片小视频在线播放| 久久精品影院6| 啪啪无遮挡十八禁网站| 精品久久蜜臀av无| 中文资源天堂在线| 12—13女人毛片做爰片一| 日本 av在线| 嫩草影院精品99| 久久午夜亚洲精品久久| 最新美女视频免费是黄的| 搡老熟女国产l中国老女人| 亚洲人成网站在线播放欧美日韩| 国产高清激情床上av| 成在线人永久免费视频| 精品一区二区三区四区五区乱码| 麻豆一二三区av精品| 精品不卡国产一区二区三区| 久久久成人免费电影| 国产伦人伦偷精品视频| xxxwww97欧美| 真人做人爱边吃奶动态| 亚洲专区中文字幕在线| 国产精品九九99| 两个人视频免费观看高清| 久9热在线精品视频| www.熟女人妻精品国产| 亚洲精品一卡2卡三卡4卡5卡| 国产男靠女视频免费网站| 免费高清视频大片| 欧美色视频一区免费| www.熟女人妻精品国产| 女同久久另类99精品国产91| 一区福利在线观看| 91久久精品国产一区二区成人 | 午夜福利在线观看吧| 色视频www国产| 一级毛片精品| 亚洲专区中文字幕在线| 90打野战视频偷拍视频| 免费在线观看成人毛片| 搡老熟女国产l中国老女人| 午夜福利高清视频| 成人特级av手机在线观看| 两性夫妻黄色片| 国产精品 欧美亚洲| 国产精品99久久久久久久久| 无遮挡黄片免费观看| 欧美日韩中文字幕国产精品一区二区三区| 又大又爽又粗| 日韩欧美国产在线观看| 亚洲成人免费电影在线观看| 久久久成人免费电影| 欧美在线一区亚洲| 在线观看免费午夜福利视频| 99国产精品99久久久久| 亚洲真实伦在线观看| 19禁男女啪啪无遮挡网站| 桃红色精品国产亚洲av| 免费观看精品视频网站| 国产午夜福利久久久久久| 国产蜜桃级精品一区二区三区| 亚洲av第一区精品v没综合| 桃色一区二区三区在线观看| 国产伦一二天堂av在线观看| 村上凉子中文字幕在线| 国产黄色小视频在线观看| 日本撒尿小便嘘嘘汇集6| 日日夜夜操网爽| 国产午夜福利久久久久久| 久久天堂一区二区三区四区| 黄片大片在线免费观看| 久久久久久久久免费视频了| 国产一区二区三区视频了| 1024香蕉在线观看| a在线观看视频网站| 蜜桃久久精品国产亚洲av| 天天一区二区日本电影三级| 美女午夜性视频免费| 精品国产乱码久久久久久男人| 亚洲人与动物交配视频| 在线国产一区二区在线| 91字幕亚洲| 99在线人妻在线中文字幕| 久久九九热精品免费| 最新美女视频免费是黄的| 日本成人三级电影网站| 男女床上黄色一级片免费看| 欧美色视频一区免费| 淫妇啪啪啪对白视频| 三级男女做爰猛烈吃奶摸视频| 国产真人三级小视频在线观看| 噜噜噜噜噜久久久久久91| 久久午夜亚洲精品久久| 亚洲 欧美一区二区三区| 亚洲av五月六月丁香网| 男女下面进入的视频免费午夜| 色综合站精品国产| 国产麻豆成人av免费视频| 99热这里只有是精品50| 在线观看66精品国产| 国产乱人伦免费视频| 成人三级黄色视频| 亚洲av第一区精品v没综合| 老司机午夜十八禁免费视频| 免费在线观看日本一区| 国产极品精品免费视频能看的| 亚洲美女视频黄频| 99久久无色码亚洲精品果冻| 婷婷丁香在线五月| 国内精品一区二区在线观看| 久久欧美精品欧美久久欧美| 久久精品影院6| 免费观看的影片在线观看| 日韩欧美在线二视频| 国产精品乱码一区二三区的特点| 一a级毛片在线观看| 日韩欧美免费精品| 很黄的视频免费| 美女被艹到高潮喷水动态| 熟女电影av网| 深夜精品福利| 麻豆国产av国片精品| 国产又黄又爽又无遮挡在线| 亚洲一区二区三区色噜噜| 宅男免费午夜| 成年人黄色毛片网站| 精品久久久久久久久久久久久| 久久久久国产精品人妻aⅴ院| 人人妻人人澡欧美一区二区| 中文字幕熟女人妻在线| 色视频www国产| 99久久99久久久精品蜜桃| 一区二区三区高清视频在线| 视频区欧美日本亚洲| 一级作爱视频免费观看| 美女被艹到高潮喷水动态| 黄色 视频免费看| 久久中文看片网| 中亚洲国语对白在线视频| 欧美另类亚洲清纯唯美| 小说图片视频综合网站| 岛国视频午夜一区免费看| 欧美黑人巨大hd| 亚洲男人的天堂狠狠| 亚洲va日本ⅴa欧美va伊人久久| 精品人妻1区二区| 欧美性猛交黑人性爽| 亚洲人与动物交配视频| 岛国在线观看网站| 亚洲av免费在线观看| 久久久色成人| 韩国av一区二区三区四区| 无遮挡黄片免费观看| 99国产精品一区二区蜜桃av| 人妻夜夜爽99麻豆av| 亚洲18禁久久av| 久久久久国内视频| 99精品在免费线老司机午夜| av天堂在线播放| 久久中文字幕人妻熟女| 美女扒开内裤让男人捅视频| 日本成人三级电影网站| 免费高清视频大片| 亚洲精品456在线播放app | 亚洲五月天丁香| av黄色大香蕉| 老汉色∧v一级毛片| 欧美黑人欧美精品刺激| 可以在线观看的亚洲视频| 国产激情偷乱视频一区二区| 又粗又爽又猛毛片免费看| 最近视频中文字幕2019在线8| 国产亚洲欧美98| 岛国在线观看网站| 免费高清视频大片| 欧美中文日本在线观看视频| 国产视频一区二区在线看| 国产精品,欧美在线| 国产亚洲精品久久久com| 日韩大尺度精品在线看网址| av在线蜜桃| 在线永久观看黄色视频| 毛片女人毛片| 成人精品一区二区免费| 一a级毛片在线观看| 99riav亚洲国产免费| 中文字幕久久专区| 美女大奶头视频| 一边摸一边抽搐一进一小说| 中文字幕高清在线视频| 黄色女人牲交| 美女 人体艺术 gogo| 一本综合久久免费| 宅男免费午夜| 欧美中文日本在线观看视频| 一级黄色大片毛片| 中文字幕熟女人妻在线| 法律面前人人平等表现在哪些方面| 国产极品精品免费视频能看的| 亚洲精品美女久久av网站| 日本撒尿小便嘘嘘汇集6| 亚洲av五月六月丁香网| 午夜福利免费观看在线| 国内久久婷婷六月综合欲色啪| 国产99白浆流出| 亚洲av日韩精品久久久久久密| 99久国产av精品| av天堂中文字幕网| 淫秽高清视频在线观看| 久久久久精品国产欧美久久久| 1000部很黄的大片| 日本a在线网址| 午夜久久久久精精品| 久久久久久人人人人人| 老鸭窝网址在线观看| 欧美一级a爱片免费观看看| 国产av在哪里看| 亚洲黑人精品在线| 美女被艹到高潮喷水动态| 少妇人妻一区二区三区视频| 精品一区二区三区视频在线 | 身体一侧抽搐| 亚洲av电影在线进入| 成年女人毛片免费观看观看9| 亚洲人成网站在线播放欧美日韩| 少妇裸体淫交视频免费看高清| 国产人伦9x9x在线观看| 国产三级在线视频| 成年人黄色毛片网站| 韩国av一区二区三区四区| 国产精品av视频在线免费观看| 真人一进一出gif抽搐免费| 欧美日本视频| 国产精品免费一区二区三区在线| 国产视频内射| 全区人妻精品视频| 成人午夜高清在线视频| 亚洲av中文字字幕乱码综合| 免费av毛片视频| 女警被强在线播放| 在线观看66精品国产| 久久久久久国产a免费观看| 特大巨黑吊av在线直播| 好男人电影高清在线观看| 最近最新中文字幕大全免费视频| 久久热在线av| 美女午夜性视频免费| 久久久色成人| 在线十欧美十亚洲十日本专区| 久久久久久久久久黄片| 色综合欧美亚洲国产小说| 亚洲精品在线美女| 欧美激情在线99| 精品国内亚洲2022精品成人| 久久这里只有精品中国| 国产97色在线日韩免费| 每晚都被弄得嗷嗷叫到高潮| 91在线精品国自产拍蜜月 | 精品一区二区三区视频在线观看免费| 亚洲av熟女| 怎么达到女性高潮| 日韩欧美一区二区三区在线观看| 性欧美人与动物交配| 国产亚洲精品久久久久久毛片| 听说在线观看完整版免费高清| 伦理电影免费视频| 亚洲国产精品999在线| 久久九九热精品免费| 搞女人的毛片| 琪琪午夜伦伦电影理论片6080| 亚洲第一欧美日韩一区二区三区| 欧美激情久久久久久爽电影| 一本久久中文字幕| 国产精品一区二区精品视频观看| 狂野欧美激情性xxxx| 白带黄色成豆腐渣| 91麻豆精品激情在线观看国产| 国产精华一区二区三区| 日本三级黄在线观看| 不卡一级毛片| 亚洲精品美女久久久久99蜜臀| 午夜精品久久久久久毛片777| 精品一区二区三区视频在线观看免费| 国产极品精品免费视频能看的| 哪里可以看免费的av片| 日韩高清综合在线| 午夜福利成人在线免费观看| 在线观看舔阴道视频| 亚洲真实伦在线观看| 久久久久久大精品| 国产成人一区二区三区免费视频网站| 亚洲国产日韩欧美精品在线观看 | 中文字幕高清在线视频| 伊人久久大香线蕉亚洲五| 久久久国产精品麻豆| 精品无人区乱码1区二区| 久久久久免费精品人妻一区二区| 欧美日本亚洲视频在线播放| 欧美成人性av电影在线观看| 色综合欧美亚洲国产小说| 亚洲精品一卡2卡三卡4卡5卡| 久久天堂一区二区三区四区| 毛片女人毛片| 国内精品久久久久精免费| 欧美+亚洲+日韩+国产| av在线天堂中文字幕| 成人鲁丝片一二三区免费| 成年版毛片免费区| 亚洲精品色激情综合| 精品熟女少妇八av免费久了| 激情在线观看视频在线高清| 琪琪午夜伦伦电影理论片6080| 午夜久久久久精精品| 91av网一区二区| 国产美女午夜福利| 成年女人看的毛片在线观看| 欧美成人一区二区免费高清观看 | 麻豆av在线久日| 欧美乱码精品一区二区三区| 免费av毛片视频| 一级作爱视频免费观看| 男人和女人高潮做爰伦理| 嫩草影院入口| 婷婷丁香在线五月| 欧美不卡视频在线免费观看| 国产精品香港三级国产av潘金莲| 真人一进一出gif抽搐免费| 小蜜桃在线观看免费完整版高清| 日韩欧美 国产精品| 成年女人永久免费观看视频| 国产亚洲av嫩草精品影院| 国产成人欧美在线观看| 日本熟妇午夜| 久久精品aⅴ一区二区三区四区| 制服丝袜大香蕉在线| 久久人妻av系列| av天堂中文字幕网| 18禁国产床啪视频网站| 国产淫片久久久久久久久 | 丰满人妻一区二区三区视频av | 成人国产一区最新在线观看| 91在线观看av| 人人妻,人人澡人人爽秒播| 18美女黄网站色大片免费观看| 亚洲天堂国产精品一区在线| 亚洲精品在线美女| 美女高潮喷水抽搐中文字幕| 国产日本99.免费观看| 亚洲在线观看片| 免费在线观看影片大全网站| 国产av不卡久久| e午夜精品久久久久久久| 国产亚洲精品综合一区在线观看| 亚洲自拍偷在线| 免费看美女性在线毛片视频| 级片在线观看| 国产99白浆流出| 九九热线精品视视频播放| 可以在线观看的亚洲视频| 熟妇人妻久久中文字幕3abv| 国产精品香港三级国产av潘金莲| 久久久久国内视频| 最新美女视频免费是黄的| 亚洲自拍偷在线| 国产一区二区在线av高清观看| 国产精品美女特级片免费视频播放器 | 国产久久久一区二区三区| 亚洲成人免费电影在线观看| 亚洲国产欧洲综合997久久,| 精品久久蜜臀av无| 成年版毛片免费区| 国产成年人精品一区二区| 99久久精品国产亚洲精品| 变态另类成人亚洲欧美熟女| 精品久久久久久久久久久久久| 99国产精品99久久久久| 丝袜人妻中文字幕| 久久性视频一级片| 免费在线观看亚洲国产| 最近在线观看免费完整版| 桃色一区二区三区在线观看| 久久午夜综合久久蜜桃| 免费无遮挡裸体视频| 一二三四社区在线视频社区8| 91久久精品国产一区二区成人 | 成人午夜高清在线视频| 中文字幕精品亚洲无线码一区| xxx96com| 麻豆久久精品国产亚洲av| 亚洲18禁久久av| 给我免费播放毛片高清在线观看| 久久国产精品人妻蜜桃| 久久久久国产精品人妻aⅴ院| av天堂在线播放| or卡值多少钱| 国产麻豆成人av免费视频| 日韩欧美国产在线观看| 在线观看免费视频日本深夜| 99久久精品热视频| 99riav亚洲国产免费| 亚洲欧美日韩高清在线视频| 国产极品精品免费视频能看的| 黄色女人牲交| 午夜福利在线观看吧| 日本在线视频免费播放| 在线观看66精品国产| 一级a爱片免费观看的视频| 日本免费a在线| 99久国产av精品| 婷婷亚洲欧美| 黄色女人牲交| 一本一本综合久久| 亚洲欧美日韩高清在线视频| 国产爱豆传媒在线观看| 国产精品爽爽va在线观看网站| 美女 人体艺术 gogo| 特级一级黄色大片| 啦啦啦免费观看视频1| 久久亚洲真实| 久久久久精品国产欧美久久久| 国产伦一二天堂av在线观看| 91av网一区二区| 国产高清有码在线观看视频| 午夜福利成人在线免费观看| 少妇熟女aⅴ在线视频| 美女午夜性视频免费| 久久亚洲真实| 国产成人av激情在线播放| 精品国产超薄肉色丝袜足j| 国产成年人精品一区二区| 日本一二三区视频观看| 99久久成人亚洲精品观看| 亚洲精品一卡2卡三卡4卡5卡| 男女下面进入的视频免费午夜| 九九热线精品视视频播放| 欧美成狂野欧美在线观看| 看片在线看免费视频| 国内精品久久久久精免费| 欧美精品啪啪一区二区三区| 91av网一区二区| 欧美黑人巨大hd| 色尼玛亚洲综合影院| 欧美日韩瑟瑟在线播放| 最好的美女福利视频网| 日本免费a在线| 国产人伦9x9x在线观看| 麻豆久久精品国产亚洲av| 90打野战视频偷拍视频| 婷婷丁香在线五月| 久久香蕉精品热| 19禁男女啪啪无遮挡网站| 日本撒尿小便嘘嘘汇集6| 日本熟妇午夜| 18禁国产床啪视频网站| av女优亚洲男人天堂 | 日韩 欧美 亚洲 中文字幕| 美女cb高潮喷水在线观看 | 免费在线观看亚洲国产| 亚洲欧美日韩无卡精品| 校园春色视频在线观看| 国产高清激情床上av| www.熟女人妻精品国产| 久久这里只有精品中国| 精品久久久久久久久久久久久| 国产私拍福利视频在线观看| 中文字幕人妻丝袜一区二区| 窝窝影院91人妻| 九九热线精品视视频播放| 精品国产乱子伦一区二区三区| 99视频精品全部免费 在线 | bbb黄色大片| 9191精品国产免费久久| 日本与韩国留学比较| 极品教师在线免费播放| 首页视频小说图片口味搜索| www日本在线高清视频| 精华霜和精华液先用哪个| 18禁裸乳无遮挡免费网站照片| 黑人操中国人逼视频| 亚洲,欧美精品.| 美女 人体艺术 gogo| 国产午夜精品论理片| 亚洲真实伦在线观看| 露出奶头的视频| 亚洲一区高清亚洲精品| 亚洲va日本ⅴa欧美va伊人久久| 国产伦精品一区二区三区视频9 | 又黄又粗又硬又大视频| 日韩 欧美 亚洲 中文字幕| 亚洲五月天丁香| 国产高清videossex| 老汉色∧v一级毛片| 精品久久久久久久毛片微露脸| 亚洲avbb在线观看| 免费在线观看影片大全网站| 51午夜福利影视在线观看| 精品久久久久久成人av| 国产精品一区二区三区四区免费观看 | h日本视频在线播放| 男女做爰动态图高潮gif福利片| cao死你这个sao货| 欧美中文综合在线视频| 亚洲 欧美 日韩 在线 免费| 欧美黑人欧美精品刺激| 最近在线观看免费完整版| 一级毛片高清免费大全| 久久人人精品亚洲av| 久久亚洲真实| 看免费av毛片| 又粗又爽又猛毛片免费看| 偷拍熟女少妇极品色| 91麻豆av在线| 国产亚洲精品一区二区www| 2021天堂中文幕一二区在线观| 岛国在线免费视频观看| 制服人妻中文乱码| 老汉色∧v一级毛片| 淫秽高清视频在线观看| 国内精品久久久久久久电影| 亚洲专区字幕在线| 久久精品夜夜夜夜夜久久蜜豆| 亚洲五月天丁香| 欧美一级a爱片免费观看看| 亚洲精品一区av在线观看| 亚洲国产精品成人综合色| 亚洲第一电影网av| 欧美三级亚洲精品| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕精品亚洲无线码一区| 老司机深夜福利视频在线观看| 国产一区二区在线观看日韩 | 人人妻人人澡欧美一区二区| 香蕉丝袜av| 国产激情久久老熟女| 午夜福利在线观看免费完整高清在 | 欧美黑人欧美精品刺激| 国产精品 欧美亚洲| 亚洲真实伦在线观看| 在线观看日韩欧美| 99视频精品全部免费 在线 | 白带黄色成豆腐渣| 免费看美女性在线毛片视频| 国产精品电影一区二区三区| 欧美黑人欧美精品刺激| 九色国产91popny在线| 九九久久精品国产亚洲av麻豆 | 老熟妇乱子伦视频在线观看| 99国产综合亚洲精品| 一级毛片高清免费大全| 久久久久久国产a免费观看| 久久热在线av| 久久亚洲精品不卡| 狂野欧美激情性xxxx| 午夜视频精品福利| 中亚洲国语对白在线视频| 一区二区三区国产精品乱码| 国产91精品成人一区二区三区| 女人被狂操c到高潮| 日韩欧美国产一区二区入口| 1024手机看黄色片| 九九在线视频观看精品| 天天躁日日操中文字幕| 国产免费男女视频| 一个人免费在线观看的高清视频| 国产探花在线观看一区二区| 国产精品一及| 日本三级黄在线观看| 亚洲成人精品中文字幕电影| 久久久久性生活片| 久久久久亚洲av毛片大全| 国产淫片久久久久久久久 | 欧美日韩乱码在线| 少妇的逼水好多| 久久久国产精品麻豆| 国产1区2区3区精品| 长腿黑丝高跟| 亚洲天堂国产精品一区在线| 国内揄拍国产精品人妻在线| 两人在一起打扑克的视频| av在线蜜桃| 精品久久久久久久末码| 亚洲最大成人中文| 91av网站免费观看| 国产精品电影一区二区三区| 美女黄网站色视频| 亚洲天堂国产精品一区在线| 日本 av在线| 色视频www国产| 最好的美女福利视频网| 桃色一区二区三区在线观看| 波多野结衣高清无吗| 国产麻豆成人av免费视频| 999久久久国产精品视频| 在线免费观看不下载黄p国产 | 人妻久久中文字幕网| 日本免费一区二区三区高清不卡| bbb黄色大片| 村上凉子中文字幕在线| 99久久国产精品久久久| 亚洲真实伦在线观看| 精品久久蜜臀av无| 啦啦啦免费观看视频1| 成熟少妇高潮喷水视频| 一a级毛片在线观看| 欧美日韩亚洲国产一区二区在线观看| 日韩 欧美 亚洲 中文字幕| 国产精品永久免费网站| 1024香蕉在线观看| 欧美最黄视频在线播放免费| 国产精品爽爽va在线观看网站| 日韩国内少妇激情av| 免费搜索国产男女视频| 国产午夜福利久久久久久| 亚洲人与动物交配视频|