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

    堆貨荷載下高樁碼頭結(jié)構(gòu)整體安全度分析

    2013-04-08 20:23:01王元戰(zhàn)
    水道港口 2013年5期
    關(guān)鍵詞:概率分布正態(tài)分布碼頭

    王元戰(zhàn),李 妲

    (天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室天津市港口與海洋工程重點(diǎn)實(shí)驗(yàn)室,天津300072)

    堆貨荷載下高樁碼頭結(jié)構(gòu)整體安全度分析

    王元戰(zhàn),李 妲

    (天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室天津市港口與海洋工程重點(diǎn)實(shí)驗(yàn)室,天津300072)

    以大型通用有限元軟件ANSYS為平臺(tái),建立了高樁碼頭結(jié)構(gòu)整體極限承載力分析的非線性有限元數(shù)值模型。將蒙特卡洛方法與有限元數(shù)值模型相結(jié)合,求取出堆載下高樁碼頭結(jié)構(gòu)整體極限承載力的概率分布及其統(tǒng)計(jì)參數(shù)。以結(jié)構(gòu)整體極限承載力(結(jié)構(gòu)抗力)和作用荷載(作用效應(yīng))為隨機(jī)變量,構(gòu)建反映高樁碼頭結(jié)構(gòu)整體安全度的功能函數(shù),采用JC法計(jì)算結(jié)構(gòu)的安全指標(biāo),建立了高樁碼頭結(jié)構(gòu)整體安全度分析的有效方法。研究了堆貨荷載下高樁碼頭結(jié)構(gòu)損傷與尺寸型式改變對(duì)極限承載力樣本概率分布與統(tǒng)計(jì)參數(shù)的影響,得出了堆載下高樁碼頭結(jié)構(gòu)極限承載力概率分布為正態(tài)分布,結(jié)構(gòu)損傷與尺寸型式改變對(duì)其統(tǒng)計(jì)參數(shù)影響不大的結(jié)論。計(jì)算表明,利用此法分析堆貨荷載下高樁碼頭結(jié)構(gòu)整體安全度簡(jiǎn)便且可行。

    高樁碼頭;堆貨荷載;整體安全度;極限承載力;概率統(tǒng)計(jì)

    可靠度理論是分析和度量結(jié)構(gòu)安全性的一種先進(jìn)手段。國(guó)內(nèi)外工程界對(duì)結(jié)構(gòu)的可靠度分析可分為構(gòu)件層次可靠度分析與結(jié)構(gòu)系統(tǒng)整體可靠度分析兩個(gè)層次。目前,對(duì)于高樁碼頭結(jié)構(gòu)體系可靠度與安全度研究工作多集中在構(gòu)件層次上,李鑫[1]對(duì)在役高樁碼頭的橫梁可靠度評(píng)估方法進(jìn)行了研究,考慮了包括氯離子、混凝土碳化及強(qiáng)度等因素隨時(shí)間變化對(duì)高樁碼頭橫梁結(jié)構(gòu)抗力的影響,對(duì)高樁碼頭的橫梁可靠度進(jìn)行計(jì)算;史青芬[2]把碼頭結(jié)構(gòu)分為結(jié)構(gòu)構(gòu)件、各分單元與整體結(jié)構(gòu)三個(gè)層次,并把碼頭結(jié)構(gòu)劃分為多個(gè)組成單元,采用分層綜合結(jié)構(gòu)模型,對(duì)高樁碼頭結(jié)構(gòu)進(jìn)行了可靠度分析,并給出安全等級(jí)評(píng)估。在高樁碼頭結(jié)構(gòu)整體可靠度分析的層次上,所開展的分析研究工作較少,仍處于探索階段。20世紀(jì)90年代初,有學(xué)者對(duì)高樁碼頭應(yīng)用模糊數(shù)學(xué)理論進(jìn)行過可靠度評(píng)估[3],隨著有限元法的廣泛應(yīng)用和可靠度計(jì)算方法的發(fā)展,亦有學(xué)者利用ANSYS對(duì)高樁墩式碼頭的可靠性指標(biāo)進(jìn)行了初步研究[4]。文獻(xiàn)[5]研究了水平荷載作用下高樁碼頭整體極限承載力的概率分布及統(tǒng)計(jì)參數(shù),將結(jié)構(gòu)的極限狀態(tài)方程顯示表達(dá)并計(jì)算水平荷載下結(jié)構(gòu)的安全指標(biāo)。但對(duì)于高樁碼頭結(jié)構(gòu),在其實(shí)際服役過程中并非僅承受水平向荷載的作用,結(jié)構(gòu)安全性更多是受豎向荷載作用的影響,更需探討結(jié)構(gòu)在豎向外部荷載(堆貨荷載等)的作用下極限承載力的統(tǒng)計(jì)分布規(guī)律,進(jìn)一步完善高樁碼頭結(jié)構(gòu)的整體安全度計(jì)算分析方法。

    本文以大型通用有限元軟件ANSYS為平臺(tái),建立了高樁碼頭結(jié)構(gòu)整體極限承載力分析的非線性有限元數(shù)值模型。將蒙特卡洛方法與有限元數(shù)值模型相結(jié)合,求取出堆載下高樁碼頭結(jié)構(gòu)整體極限承載力的概率分布及其統(tǒng)計(jì)參數(shù)。以結(jié)構(gòu)整體極限承載力(結(jié)構(gòu)抗力)和作用荷載(作用效應(yīng))為隨機(jī)變量,構(gòu)建反映高樁碼頭結(jié)構(gòu)整體安全度的功能函數(shù),采用JC法計(jì)算結(jié)構(gòu)的安全指標(biāo),建立了高樁碼頭結(jié)構(gòu)整體安全度分析的有效方法。研究了堆貨荷載下高樁碼頭結(jié)構(gòu)損傷與尺寸型式改變對(duì)極限承載力樣本概率分布與統(tǒng)計(jì)參數(shù)的影響,得出了堆載下高樁碼頭結(jié)構(gòu)極限承載力概率分布為正態(tài)分布,結(jié)構(gòu)損傷與尺寸型式改變對(duì)其統(tǒng)計(jì)參數(shù)影響不大的結(jié)論。計(jì)算表明,利用此法分析堆貨荷載下高樁碼頭結(jié)構(gòu)整體安全度簡(jiǎn)便且可行。

    1 基本分析思路

    高樁碼頭整體極限承載力代表了高樁碼頭結(jié)構(gòu)整體的承載性能和安全狀況。本文以高樁碼頭結(jié)構(gòu)整體極限承載力及作用荷載作為基本的隨機(jī)變量,可建立反映高樁碼頭結(jié)構(gòu)整體安全度的功能函數(shù)如下

    式中:R代表高樁碼頭整體極限承載力;S代表作用于高樁碼頭結(jié)構(gòu)上的荷載。

    在式(1)中,作用于高樁碼頭結(jié)構(gòu)上的荷載作為隨機(jī)變量考慮,其概率分布與統(tǒng)計(jì)參數(shù)已在文獻(xiàn)[6]中說明,而高樁碼頭整體極限承載力的概率分布與統(tǒng)計(jì)參數(shù)現(xiàn)僅有水平向荷載作用下的參考數(shù)據(jù),尚無(wú)堆貨荷載下的相關(guān)統(tǒng)計(jì)數(shù)據(jù)作為計(jì)算參考。為解決這一問題,本文利用蒙特卡洛模擬與非線性有限元結(jié)合,模擬出堆貨荷載作用下高樁碼頭整體極限承載力的概率分布與統(tǒng)計(jì)參數(shù),根據(jù)所得參數(shù)使用JC法將非正態(tài)分布的設(shè)計(jì)變量當(dāng)量化為正態(tài)分布的設(shè)計(jì)變量,隨后即可采用工程中常用的計(jì)算方法對(duì)式(1)計(jì)算得出可靠指標(biāo)。

    2 極限承載力概率分布與統(tǒng)計(jì)參數(shù)

    2.1 極限承載力樣本求取

    2.1.1 模型的輸入變量樣本

    影響結(jié)構(gòu)極限承載力的因素有很多,如荷載種類、結(jié)構(gòu)尺寸、材料性能、計(jì)算方法等。本文考慮的荷載為堆貨荷載,并考慮結(jié)構(gòu)幾何尺寸、材料性能等影響因素如表1所示。利用ANSYS中PDS模塊提供的具有較高抽樣效率的拉丁超立方抽樣方法對(duì)表1中的影響因素進(jìn)行抽樣,形成求取極限承載力樣本模型的輸入變量樣本,樣本容量為400。

    2.1.2 有限元模型概況

    選取典型梁板式高樁碼頭標(biāo)準(zhǔn)結(jié)構(gòu)段為研究對(duì)象,其結(jié)構(gòu)斷面如圖1所示。本碼頭模型中鋼筋混凝土材料采用整體式的SOLID65單元模擬?;炷敛牧锨?zhǔn)則采用多線性等向強(qiáng)化模型(MISO),其本構(gòu)關(guān)系曲線按現(xiàn)行規(guī)范[7]確定,破壞準(zhǔn)則為改進(jìn)的William-Warnke五參數(shù)破壞準(zhǔn)則,鋼筋材料屈服準(zhǔn)則采用雙線性等向強(qiáng)化模型(BISO)。樁側(cè)土體采用SOLID45單元模擬,采用Drucker-Prager屈服準(zhǔn)則。為更好的模擬樁-土相互作用,本模型在樁與土體之間建立接觸單元(TARGE170和CONTA174)來(lái)模擬其粘結(jié)、滑移和脫離等現(xiàn)象。

    2.1.3 極限承載力求解方法

    對(duì)高樁碼頭模型進(jìn)行有限元數(shù)值分析,通過逐步施加荷載增量來(lái)求解結(jié)構(gòu)的繼續(xù)承載力。考慮幾何非線性和材料非線性的高樁碼頭結(jié)構(gòu)增量平衡方程為

    在每個(gè)荷載增量步內(nèi),對(duì)式(2)按照增量-Newton Raphson迭代法進(jìn)行求解,并輔以線性搜索及自適應(yīng)下降功能,來(lái)保證求解速度和收斂穩(wěn)定性[8]。隨著荷載的逐步施加,高樁碼頭結(jié)構(gòu)的剛度隨之發(fā)生變化,當(dāng)荷載產(chǎn)生的應(yīng)力使切線剛度矩陣趨于奇異時(shí),結(jié)構(gòu)無(wú)法繼續(xù)承載,此時(shí)的荷載即為極限承載力[9]。在碼頭上部結(jié)構(gòu)上選取一特征點(diǎn)(位置在圖1中標(biāo)出),可繪制豎向荷載—位移曲線如圖2。從圖2中可看出,加載后期,荷載—位移曲線的斜率已趨于零,即荷載增加很小,結(jié)構(gòu)也會(huì)產(chǎn)生較大的位移,可以判定結(jié)構(gòu)已經(jīng)達(dá)到極限承載力。由此可見,本文選取的求解結(jié)構(gòu)極限承載力的方法是可行且較準(zhǔn)確的。

    2.1.4 極限承載力樣本求取流程

    本文利用ANSYS參數(shù)化設(shè)計(jì)語(yǔ)言(APDL)編寫程序,并運(yùn)用宏技術(shù)將求解極限承載力樣本的各功能進(jìn)行分塊,將影響因素抽樣、有限元分析及樣本輸出模塊形成的獨(dú)立宏文件與主程序連接,具有較好的穩(wěn)定性。整個(gè)求解程序的基本流程如圖3所示。

    2.2 樣本概率統(tǒng)計(jì)分析

    2.2.1 工況設(shè)置

    目前,實(shí)際在服役的高樁碼頭結(jié)構(gòu)普遍存在不同程度的損傷,且其實(shí)際結(jié)構(gòu)型式或設(shè)計(jì)尺寸會(huì)與前文介紹的典型結(jié)構(gòu)段有一定區(qū)別。為研究結(jié)構(gòu)損傷以及結(jié)構(gòu)型式與尺寸變化對(duì)極限承載力概率分布的影響,本文除設(shè)置了典型尺寸與型式的無(wú)損傷結(jié)構(gòu)工況外,還按照結(jié)構(gòu)損傷情況、結(jié)構(gòu)尺寸改變情況和結(jié)構(gòu)型式改變(圖4)設(shè)計(jì)組合工況。工況設(shè)置的具體情況詳見表2。根據(jù)文獻(xiàn)[10-11],高樁碼頭多于樁頂與樁帽連接處、橫梁跨中底部出現(xiàn)裂縫等損傷,面板處多現(xiàn)大面積剝蝕破壞,本文采用直接減小截面尺寸的方法來(lái)模擬面板大面積的剝蝕破壞,采用在結(jié)構(gòu)易出現(xiàn)裂縫的位置處設(shè)置缺陷來(lái)模擬損傷,損傷模擬示意如圖4。

    2.2.2 概率分布檢驗(yàn)與確定

    得到不同工況下的極限承載力樣本值后,需將樣本值進(jìn)行無(wú)量綱化處理,使其更具一般性。處理方法為用樣本值除以其所對(duì)應(yīng)工況下的極限承載力標(biāo)準(zhǔn)值,這里的極限承載力標(biāo)準(zhǔn)值為各輸入隨機(jī)變量取標(biāo)準(zhǔn)值時(shí)計(jì)算所得的極限承載力值。為進(jìn)一步減少確定樣本分布規(guī)律的干擾因素,本文采用拉依達(dá)準(zhǔn)則[12]對(duì)無(wú)量綱化后的樣本值進(jìn)行異常值的剔除。

    為確定所得樣本的概率分布,以工況(一)為例,首先繪制無(wú)量綱化后的樣本分組直方圖如圖5所示,觀察其概率密度曲線與工程中常見的正態(tài)分布、對(duì)數(shù)正態(tài)分布和極值I型分布相似,故采用概率紙檢驗(yàn)的方法來(lái)觀察樣本點(diǎn)在指定分布概率紙上的分布情況,工況(一)的檢驗(yàn)結(jié)果如圖6所示,可以直觀的看出樣本點(diǎn)與正態(tài)累積概率曲線、對(duì)數(shù)正態(tài)累積概率曲線較接近,基本呈直線分布,與極值I型累積概率曲線擬合較差,尤其在尾部樣本點(diǎn)曲線較遠(yuǎn)。工況(二)到工況(十四)的概率紙檢驗(yàn)結(jié)果與工況(一)類似,限于篇幅原因不在此一一列出?;谝陨咸攸c(diǎn),可以認(rèn)為各工況下樣本服從正態(tài)分布和對(duì)數(shù)正態(tài)分布的可能性要大于其服從極值I型分布。

    為進(jìn)一步選擇最優(yōu)分布,分別采用χ2擬合優(yōu)度檢驗(yàn)、Kolmogorov-Smirnov檢驗(yàn)、Lilliefors檢驗(yàn)、Cramer-Von Mises檢驗(yàn)[13]及A-D檢驗(yàn)[14]共五種假設(shè)檢驗(yàn)方法對(duì)各工況下所得極限承載力樣本進(jìn)行分析,假設(shè)檢驗(yàn)結(jié)果匯總于表3,其中,A代表正態(tài)分布,B代表對(duì)數(shù)正態(tài)分布;0代表接受假設(shè)檢驗(yàn),1代表拒絕假設(shè)檢驗(yàn);p值為在顯著性水平為0.05的情況下不拒絕接受該假設(shè)檢驗(yàn)的檢驗(yàn)值,p值越大表示不拒絕接受假設(shè)檢驗(yàn)的可能性越大;m值為在顯著性水平為0.05的情況下不拒絕接受該假設(shè)檢驗(yàn)的統(tǒng)計(jì)值,與相應(yīng)的臨界統(tǒng)計(jì)值比較,m值越小表示不拒絕接受假設(shè)檢驗(yàn)的可能性越大。

    從表3中可看出,有六種工況正態(tài)分布接受假設(shè)檢驗(yàn)方法的次數(shù)較對(duì)數(shù)正態(tài)分布多一些,剩余八種工況下,兩種分布接受假設(shè)檢驗(yàn)方法的次數(shù)一樣多,但進(jìn)一步比較正態(tài)分布與對(duì)數(shù)正態(tài)分布的檢驗(yàn)p值與檢驗(yàn)統(tǒng)計(jì)m值可知,正態(tài)分布的檢驗(yàn)p值更大,檢驗(yàn)統(tǒng)計(jì)m值更小。由此可認(rèn)為,正態(tài)分布為各工況下極限承載力樣本的最優(yōu)分布。為進(jìn)一步驗(yàn)證上述結(jié)果的可信度,采用專門用于檢驗(yàn)樣本是否服從正態(tài)分布的d′Agostino檢驗(yàn)法與Jarque-Bera檢驗(yàn)法[15]對(duì)上述十四種工況極限承載力樣本進(jìn)行檢驗(yàn),所得結(jié)果均不拒絕接受原假設(shè)。

    綜合以上假設(shè)檢驗(yàn)結(jié)果可知,在堆貨荷載作用下,高樁碼頭結(jié)構(gòu)損傷、尺寸與型式改變不影響極限承載力的概率分布類型,可確定堆載下高樁碼頭極限承載力概率分布為正態(tài)分布。

    2.2.3 統(tǒng)計(jì)參數(shù)確定

    利用極大似然估計(jì)法統(tǒng)計(jì)計(jì)算,并剔除異常樣本值后,得到各工況堆貨荷載作用下高樁碼頭結(jié)構(gòu)極限承載力樣本統(tǒng)計(jì)參數(shù)如表4所示。比較工況(一)~(九)可知,結(jié)構(gòu)損傷對(duì)極限承載力統(tǒng)計(jì)參數(shù)略有影響,從變異系數(shù)的比較中可看出,有損結(jié)構(gòu)極限承載力統(tǒng)計(jì)樣本值的離散程度較無(wú)損結(jié)構(gòu)略大,且損傷程度越高、損傷數(shù)量越多,結(jié)構(gòu)的極限承載力統(tǒng)計(jì)樣本值離散程度越大。分別將工況(十)~(十四)與工況(一)比較,總體上各工況下樣本統(tǒng)計(jì)參數(shù)比較接近,相對(duì)誤差均在5%以內(nèi),可知結(jié)構(gòu)尺寸與型式改變對(duì)極限承載力統(tǒng)計(jì)參數(shù)影響不大。

    3 堆載下高樁碼頭結(jié)構(gòu)整體安全度分析

    通過計(jì)算高樁碼頭整體可靠指標(biāo)來(lái)分析結(jié)構(gòu)的安全度,功能函數(shù)用式(1)表達(dá),式中的兩個(gè)隨機(jī)變量的平均值和標(biāo)準(zhǔn)差()分別為

    港口工程堆貨荷載的概率分布為極值I型分布[6],首先使用JC法將其當(dāng)量化為正態(tài)分布的設(shè)計(jì)變量,隨后采用工程中常用的改進(jìn)一次二階矩法[16]對(duì)式(1)計(jì)算各工況的實(shí)際可靠指標(biāo)如表5所示,結(jié)果顯示結(jié)構(gòu)極限承載能力隨結(jié)構(gòu)損傷程度與損傷數(shù)量的增加而下降,高度損傷結(jié)構(gòu)的實(shí)際可靠指標(biāo)與無(wú)損結(jié)構(gòu)相比有明顯降低,可認(rèn)為結(jié)構(gòu)的整體安全度明顯下降,這與工程實(shí)際中的認(rèn)識(shí)是一致的。

    若采用無(wú)損典型尺寸結(jié)構(gòu)極限承載力概率分布的統(tǒng)計(jì)參數(shù)(工況一)計(jì)算其他工況下結(jié)構(gòu)的可靠指標(biāo),并將其與實(shí)際可靠指標(biāo)進(jìn)行對(duì)比,計(jì)算其相對(duì)誤差均小于5%,詳見表5。故在堆載下對(duì)高樁碼頭結(jié)構(gòu)進(jìn)行整體安全度分析時(shí),可采用無(wú)損結(jié)構(gòu)極限承載力的概率分布及統(tǒng)計(jì)參數(shù)對(duì)實(shí)際工程結(jié)構(gòu)進(jìn)行分析計(jì)算,計(jì)算結(jié)果精度可以滿足工程要求。

    4 結(jié)語(yǔ)

    本文以可靠度理論為基礎(chǔ),結(jié)合結(jié)構(gòu)極限承載力的相關(guān)內(nèi)容,將高樁碼頭結(jié)構(gòu)整體極限承載力及作用荷載作為基本的隨機(jī)變量,建立了反映高樁碼頭結(jié)構(gòu)整體安全度的功能函數(shù)。利用ANSYS有限元軟件對(duì)堆載下高樁碼頭結(jié)構(gòu)極限承載力進(jìn)行了大量樣本統(tǒng)計(jì)分析,得出了堆載下高樁碼頭結(jié)構(gòu)極限承載力概率分布為正態(tài)分布,結(jié)構(gòu)損傷與尺寸型式改變對(duì)其統(tǒng)計(jì)參數(shù)影響不大的科學(xué)結(jié)論。在此基礎(chǔ)上,對(duì)堆貨荷載作用下高樁碼頭結(jié)構(gòu)進(jìn)行了整體安全度的分析,結(jié)果表明利用此法分析堆載下高樁碼頭結(jié)構(gòu)整體安全度簡(jiǎn)便且可行。

    [1]李鑫.在役高樁碼頭橫梁可靠度評(píng)估方法研究[D].天津:天津大學(xué),2010.

    [2]史青芬.高樁碼頭結(jié)構(gòu)安全性評(píng)估[D].重慶:重慶交通大學(xué),2010.

    [3]竺存宏,韓玉,劉衛(wèi)民.海港高樁碼頭結(jié)構(gòu)可靠度的模糊數(shù)學(xué)評(píng)估[J].水道港口,1992(4):1-11. ZHU C H,HAN Y,LIUW M.Fuzzy Mathematical Evaluation for the Reliability of the Structure of Wharf in Seaport[J].Journal of Waterway and Harbor,1992(4):1-11.

    [4]李聲文,張日向,韓麗華.基于ANSYS的高樁墩式碼頭可靠度分析[J].水道港口,2011,32(4):276-280. LISW,ZHANGR X,HAN L H.Reliability analysis of pile pier based on ANSYS[J].Journal of Waterway and Harbor,2011,32(4):276-280.

    [5]王元戰(zhàn),高樹奇,張智凱,等.梁板式高樁碼頭極限承載力概率分布[J].水道港口,2012(6):515-522. WANG Y Z,GAO SQ,ZHANG Z K.Probability distribution of ultimate bearing capacity of beam-slab high-piled wharf[J]. Journal of Waterway and Harbor,2012(6):515-522.

    [6]港口工程結(jié)構(gòu)可靠度設(shè)計(jì)統(tǒng)一標(biāo)準(zhǔn)編制組.港口工程結(jié)構(gòu)可靠度[M].北京:人民交通出版社,1992.

    [7]GB50010-2002,混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范[S].

    [8]程進(jìn),江見鯨,肖汝誠(chéng),等.大跨度鋼拱橋結(jié)構(gòu)極限承載力分析[J].工程力學(xué),2003,20(2):7-10. CHENG J,JIANG J J,XIAO R C,et al.ULTIMATE LOAD-CARRYING CAPACITY OF LONG-SPAN STEEL ARCH BRIDGES[J].ENGINEERING MECHANICS,2003,20(2):7-10.

    [9]劉占省.空間結(jié)構(gòu)體系可靠度分析及關(guān)鍵構(gòu)件的判定與應(yīng)用[D].哈爾濱:哈爾濱工程大學(xué),2007.

    [10]王能貝.某高樁碼頭損傷原因及承載能力研究[D].長(zhǎng)沙:長(zhǎng)沙理工大學(xué),2010.

    [11]盧秀榮.天津港高樁碼頭損壞原因分析及修復(fù)方法[J].水道港口,1999(3):39-42. LU X R.On causes of failure of open type pier on piles of Tianjin port and its repairing method[J].Journal of waterway and harbor,1999(3):39-42.

    [12]肖明耀.誤差理論與應(yīng)用[M].北京:計(jì)量出版社,1985:170-171.

    [13]張思俊,吳世偉,甘峰,等.結(jié)構(gòu)可靠度分析中隨機(jī)變量分布函數(shù)的擬合檢驗(yàn)[C]//吳世偉.結(jié)構(gòu)安全度與可靠度分析論文集.南京:河海大學(xué)出版社,1988.

    [14]P A W Lewis.Distribution of the Anderson-Darling Statistic[J].The Annals of Mathematical Statistics,1961,32(4):1 118-1 124.

    [15]梁馮珍,宋占杰,張玉環(huán).應(yīng)用概率統(tǒng)計(jì)[M].天津:天津大學(xué)出版社,2004.

    [16]GB50158-2010,港口工程結(jié)構(gòu)可靠性設(shè)計(jì)統(tǒng)一標(biāo)準(zhǔn)[S].

    Integral safety degree analysis of high-piled wharf structure under heaped load

    WANG Yuan-zhan,LIDa
    (National Key Laboratory of Water Conservancy Engineering Simulation and Security,Tianjin Key Laboratory of Port and Ocean Engineering,Tianjin University,Tianjin 300072,China)

    Using the large universal finite element software ANSYS as a platform,a nonlinear finite element model of high-piled wharf was established to analyze its ultimate bearing capacity in this paper.The probability distribution of ultimate bearing capacity of high-piled wharf under heaped load and its statistical parameter were obtained by combination of the Monte Carlo method and the finite element numerical model.A effective method was also established to analyze the integral safety degree of high-piled wharf,which could calculate the safety index of structure by using JC-method.This method proposed a function reflected the integral safety degree of high-piled wharf by taking the ultimate bearing capacity of the structure(structural resistance)and the load(load action effect)as random variables.The probability distribution of ultimate bearing capacity of high-piled wharf was analyzed under heaped load.The impact of probability distribution and statistical parameters of ultimate bearing capacity of the damage,the changes of the structure′s sizes and type were considered.The results show that the probability distribution of the ultimate bearing capacity of high-piled wharf under heaped load are normal distribution and its statistical parameter is not affected by the damage,the size,or the type change of the structure.The calculation shows that it is simple and feasible to analyze the integral safety of high-piled wharf by using this method.

    high-piled wharf;heaped load;integral safety degree;ultimate bearing capacity;probability statistics

    U 656.1+13

    A

    1005-8443(2013)05-0430-07

    海南將打造東南亞航運(yùn)樞紐

    2013-04-10;

    2013-05-10

    國(guó)家自然科學(xué)基金(51279128)

    王元戰(zhàn)(1958-),男,天津市人,教授,博士生導(dǎo)師,主要從事港口海岸及近海工程結(jié)構(gòu)設(shè)計(jì)理論和方法、巖土力學(xué)及土與結(jié)構(gòu)相互作用、結(jié)構(gòu)振動(dòng)分析理論和方法等方面的研究工作。

    Biography:WANG Yuan-zhan(1958-),male,professor.

    本刊從海南省交通運(yùn)輸廳獲悉,到2020年,海南將建成港口設(shè)施完善、航運(yùn)資源集聚、航運(yùn)服務(wù)主要功能完備、航運(yùn)秩序規(guī)范的,以??诟?、洋浦港為雙核的,背靠華南腹地,面向東南亞,具有資源配置能力的區(qū)域性航運(yùn)樞紐。據(jù)介紹,國(guó)際旅游島建設(shè)上升為國(guó)家戰(zhàn)略以來(lái),海南交通物流業(yè)有了初步發(fā)展。目前,全省擁有萬(wàn)噸級(jí)以上深水泊位34個(gè),即將建成的深水泊位9個(gè),港口吞吐能力明顯增強(qiáng)。(殷缶,梅深)

    猜你喜歡
    概率分布正態(tài)分布碼頭
    全自動(dòng)化碼頭來(lái)了
    離散型概率分布的ORB圖像特征點(diǎn)誤匹配剔除算法
    關(guān)于概率分布函數(shù)定義的辨析
    科技視界(2016年19期)2017-05-18 10:18:46
    基于對(duì)數(shù)正態(tài)分布的出行時(shí)長(zhǎng)可靠性計(jì)算
    前往碼頭
    基于概率分布的PPP項(xiàng)目風(fēng)險(xiǎn)承擔(dān)支出測(cè)算
    正態(tài)分布及其應(yīng)用
    在碼頭上釣魚
    正態(tài)分布題型剖析
    χ2分布、t 分布、F 分布與正態(tài)分布間的關(guān)系
    久久天堂一区二区三区四区| 国产一区二区三区综合在线观看| 欧美日本亚洲视频在线播放| 国产人伦9x9x在线观看| 久久久国产精品麻豆| 亚洲专区国产一区二区| 免费久久久久久久精品成人欧美视频| 亚洲精品中文字幕一二三四区| 在线永久观看黄色视频| 亚洲自拍偷在线| 悠悠久久av| 悠悠久久av| 女同久久另类99精品国产91| 又紧又爽又黄一区二区| 热99国产精品久久久久久7| 欧美日韩国产mv在线观看视频| 99国产精品99久久久久| 日日摸夜夜添夜夜添小说| 亚洲精品国产一区二区精华液| 亚洲性夜色夜夜综合| 亚洲欧美精品综合一区二区三区| 两人在一起打扑克的视频| 国产免费现黄频在线看| 国产真人三级小视频在线观看| 激情视频va一区二区三区| 好看av亚洲va欧美ⅴa在| 久久久久精品国产欧美久久久| 极品人妻少妇av视频| 18美女黄网站色大片免费观看| 成人18禁在线播放| 1024香蕉在线观看| 日本五十路高清| 久久精品91蜜桃| 一本综合久久免费| 美女扒开内裤让男人捅视频| 国产熟女午夜一区二区三区| 黄频高清免费视频| www.999成人在线观看| 午夜福利欧美成人| 精品电影一区二区在线| 国产精品98久久久久久宅男小说| 久久国产亚洲av麻豆专区| 久久久久精品国产欧美久久久| 99久久精品国产亚洲精品| 日日干狠狠操夜夜爽| 午夜福利欧美成人| 午夜两性在线视频| 国产精品99久久99久久久不卡| 99国产精品一区二区三区| 久久久久国产精品人妻aⅴ院| 国产精品国产高清国产av| 日韩免费高清中文字幕av| 亚洲av片天天在线观看| 亚洲精品久久午夜乱码| 夜夜爽天天搞| 久久精品aⅴ一区二区三区四区| av在线播放免费不卡| 欧美av亚洲av综合av国产av| 天天躁狠狠躁夜夜躁狠狠躁| 自线自在国产av| 99久久久亚洲精品蜜臀av| 精品电影一区二区在线| 波多野结衣高清无吗| 国产亚洲精品第一综合不卡| 亚洲国产精品sss在线观看 | 欧美日韩亚洲高清精品| 亚洲国产欧美日韩在线播放| 麻豆久久精品国产亚洲av | 香蕉丝袜av| 超碰97精品在线观看| 天堂影院成人在线观看| 桃红色精品国产亚洲av| 国产一区二区三区在线臀色熟女 | 国产国语露脸激情在线看| 99国产精品一区二区蜜桃av| 久久香蕉国产精品| 黄色视频,在线免费观看| 久久性视频一级片| 亚洲五月婷婷丁香| 真人一进一出gif抽搐免费| 午夜精品在线福利| aaaaa片日本免费| 国产区一区二久久| 色婷婷av一区二区三区视频| 亚洲成国产人片在线观看| 亚洲免费av在线视频| 国产高清国产精品国产三级| 久久香蕉激情| 国产在线观看jvid| 国产aⅴ精品一区二区三区波| 亚洲欧美一区二区三区黑人| 国产精品二区激情视频| 国产av一区二区精品久久| 搡老岳熟女国产| 村上凉子中文字幕在线| 97碰自拍视频| 欧美+亚洲+日韩+国产| 丝袜在线中文字幕| 不卡一级毛片| 精品福利观看| a级毛片黄视频| 亚洲精品av麻豆狂野| 精品一区二区三区av网在线观看| 欧美中文综合在线视频| 女人精品久久久久毛片| 黄色片一级片一级黄色片| 欧美日韩一级在线毛片| 99riav亚洲国产免费| 亚洲黑人精品在线| 狠狠狠狠99中文字幕| 精品久久久久久久毛片微露脸| 丰满的人妻完整版| 久久亚洲精品不卡| 久久精品国产亚洲av香蕉五月| 色在线成人网| 国产三级在线视频| 亚洲色图av天堂| 亚洲精品av麻豆狂野| 日韩大尺度精品在线看网址 | 高清在线国产一区| 一级毛片女人18水好多| 亚洲在线自拍视频| 脱女人内裤的视频| 99久久99久久久精品蜜桃| 脱女人内裤的视频| 一区二区三区国产精品乱码| 脱女人内裤的视频| 老司机深夜福利视频在线观看| 精品国产亚洲在线| 日本欧美视频一区| 亚洲第一青青草原| 深夜精品福利| 嫩草影院精品99| 午夜福利在线观看吧| 极品人妻少妇av视频| 免费久久久久久久精品成人欧美视频| 一级毛片高清免费大全| 99国产精品一区二区三区| 不卡一级毛片| 欧美在线黄色| a级片在线免费高清观看视频| 啦啦啦 在线观看视频| 亚洲黑人精品在线| 国产熟女午夜一区二区三区| 一级片'在线观看视频| 高清黄色对白视频在线免费看| 国产一区二区激情短视频| 19禁男女啪啪无遮挡网站| 波多野结衣av一区二区av| 亚洲精品在线美女| 国产精品一区二区精品视频观看| 亚洲精品成人av观看孕妇| 欧美日韩av久久| 超碰97精品在线观看| 国产野战对白在线观看| 日韩欧美一区视频在线观看| 亚洲成av片中文字幕在线观看| 狠狠狠狠99中文字幕| 成人黄色视频免费在线看| 欧美黄色淫秽网站| 日韩大尺度精品在线看网址 | 国产免费男女视频| 亚洲午夜精品一区,二区,三区| 一二三四在线观看免费中文在| 香蕉久久夜色| 香蕉久久夜色| 女人爽到高潮嗷嗷叫在线视频| 国产成人免费无遮挡视频| 亚洲av成人一区二区三| 一级毛片高清免费大全| 欧美人与性动交α欧美软件| 老司机亚洲免费影院| 亚洲精品成人av观看孕妇| 免费在线观看影片大全网站| 神马国产精品三级电影在线观看 | 欧美乱色亚洲激情| 九色亚洲精品在线播放| 精品少妇一区二区三区视频日本电影| 欧美亚洲日本最大视频资源| 涩涩av久久男人的天堂| 精品国产一区二区三区四区第35| 91麻豆精品激情在线观看国产 | 人妻久久中文字幕网| 法律面前人人平等表现在哪些方面| 国产精品永久免费网站| 国产aⅴ精品一区二区三区波| 人妻丰满熟妇av一区二区三区| 精品高清国产在线一区| 99国产精品一区二区三区| 国产精品1区2区在线观看.| 黄网站色视频无遮挡免费观看| 老熟妇仑乱视频hdxx| av天堂在线播放| 巨乳人妻的诱惑在线观看| 欧美午夜高清在线| 国产黄a三级三级三级人| 黄色视频,在线免费观看| 久久精品91无色码中文字幕| av福利片在线| 精品午夜福利视频在线观看一区| 黄色丝袜av网址大全| 91大片在线观看| 国产激情欧美一区二区| 国产乱人伦免费视频| 亚洲av成人不卡在线观看播放网| 岛国视频午夜一区免费看| 日韩三级视频一区二区三区| 免费看十八禁软件| 91麻豆av在线| 亚洲国产欧美网| 国产伦一二天堂av在线观看| 久久久久国产精品人妻aⅴ院| 999久久久精品免费观看国产| 精品久久久久久,| 亚洲精品国产精品久久久不卡| 久久久国产成人免费| 亚洲精品久久午夜乱码| 在线观看www视频免费| 欧美+亚洲+日韩+国产| 美女高潮喷水抽搐中文字幕| 久久亚洲真实| 欧美一级毛片孕妇| 超碰成人久久| 日本免费一区二区三区高清不卡 | 亚洲精品一区av在线观看| 欧美日本中文国产一区发布| 午夜免费激情av| 91大片在线观看| 曰老女人黄片| 日韩成人在线观看一区二区三区| 五月开心婷婷网| 高清毛片免费观看视频网站 | 一级a爱视频在线免费观看| 成年版毛片免费区| 精品电影一区二区在线| 99久久久亚洲精品蜜臀av| 一区二区日韩欧美中文字幕| 啦啦啦免费观看视频1| 亚洲国产看品久久| 男人舔女人下体高潮全视频| 老汉色av国产亚洲站长工具| 叶爱在线成人免费视频播放| 别揉我奶头~嗯~啊~动态视频| 黄色怎么调成土黄色| 黄网站色视频无遮挡免费观看| 日本精品一区二区三区蜜桃| 久久久久国内视频| 法律面前人人平等表现在哪些方面| 欧美日韩瑟瑟在线播放| 国产又爽黄色视频| 国产黄a三级三级三级人| 十八禁人妻一区二区| 正在播放国产对白刺激| 亚洲精品成人av观看孕妇| 两个人免费观看高清视频| 亚洲黑人精品在线| 国产成人欧美在线观看| 亚洲av成人一区二区三| 在线观看免费视频日本深夜| 国产精品亚洲一级av第二区| 欧美日韩亚洲高清精品| 久久精品国产亚洲av香蕉五月| 天堂动漫精品| 精品免费久久久久久久清纯| 69精品国产乱码久久久| √禁漫天堂资源中文www| 看黄色毛片网站| 欧美成人午夜精品| 精品国产一区二区三区四区第35| 美女高潮到喷水免费观看| 日韩高清综合在线| e午夜精品久久久久久久| 亚洲欧美日韩高清在线视频| 无限看片的www在线观看| 波多野结衣高清无吗| 在线永久观看黄色视频| 叶爱在线成人免费视频播放| 成人黄色视频免费在线看| 丰满人妻熟妇乱又伦精品不卡| 美女福利国产在线| 侵犯人妻中文字幕一二三四区| 在线视频色国产色| 在线天堂中文资源库| 啪啪无遮挡十八禁网站| 97超级碰碰碰精品色视频在线观看| 精品第一国产精品| 精品福利观看| 国产在线观看jvid| 成在线人永久免费视频| 桃色一区二区三区在线观看| 老司机亚洲免费影院| 久久亚洲精品不卡| 这个男人来自地球电影免费观看| 色尼玛亚洲综合影院| 国产精品乱码一区二三区的特点 | 久久人妻av系列| 又黄又爽又免费观看的视频| 搡老岳熟女国产| 亚洲精华国产精华精| 国产高清videossex| 久热爱精品视频在线9| 亚洲男人天堂网一区| www.www免费av| 两个人看的免费小视频| 欧美成人午夜精品| 中文字幕人妻丝袜一区二区| 亚洲精华国产精华精| 欧美中文日本在线观看视频| xxxhd国产人妻xxx| 亚洲精品中文字幕在线视频| 18禁国产床啪视频网站| 久久精品亚洲av国产电影网| av天堂久久9| 日韩中文字幕欧美一区二区| 乱人伦中国视频| 色综合站精品国产| 中文字幕人妻熟女乱码| 精品国产国语对白av| 亚洲美女黄片视频| 欧美在线一区亚洲| 国产片内射在线| 久久久久久人人人人人| 国产成人精品在线电影| 国产区一区二久久| 色老头精品视频在线观看| 十八禁网站免费在线| 久久久精品国产亚洲av高清涩受| 黑人操中国人逼视频| 国产真人三级小视频在线观看| 又大又爽又粗| 91精品国产国语对白视频| 日本五十路高清| cao死你这个sao货| 久久热在线av| 亚洲专区字幕在线| 在线观看一区二区三区| 日本黄色日本黄色录像| 欧美色视频一区免费| 欧美日韩亚洲国产一区二区在线观看| 亚洲一区二区三区不卡视频| 天堂俺去俺来也www色官网| 久久亚洲精品不卡| 日韩视频一区二区在线观看| 亚洲第一av免费看| 亚洲熟妇中文字幕五十中出 | 又黄又粗又硬又大视频| 成人三级黄色视频| av网站免费在线观看视频| 国产成人啪精品午夜网站| 欧美日韩av久久| 宅男免费午夜| 精品熟女少妇八av免费久了| 黑人操中国人逼视频| 久久狼人影院| 久久天躁狠狠躁夜夜2o2o| 50天的宝宝边吃奶边哭怎么回事| 又紧又爽又黄一区二区| 最好的美女福利视频网| 亚洲国产精品999在线| 亚洲精品久久成人aⅴ小说| 国产一区二区在线av高清观看| 国产国语露脸激情在线看| 在线观看免费高清a一片| 日韩精品青青久久久久久| 久久精品人人爽人人爽视色| 久久这里只有精品19| 黑人操中国人逼视频| 亚洲成人免费av在线播放| 中文字幕人妻丝袜制服| 久久草成人影院| 日韩精品中文字幕看吧| 搡老岳熟女国产| 最近最新免费中文字幕在线| 国内久久婷婷六月综合欲色啪| 亚洲精品一卡2卡三卡4卡5卡| av中文乱码字幕在线| 国产亚洲精品久久久久5区| 国产97色在线日韩免费| 欧美乱色亚洲激情| 日韩精品免费视频一区二区三区| 男女下面插进去视频免费观看| 色尼玛亚洲综合影院| 九色亚洲精品在线播放| 国产精品久久久人人做人人爽| 国产黄色免费在线视频| 成人影院久久| 天堂俺去俺来也www色官网| 欧美黑人欧美精品刺激| 在线观看免费午夜福利视频| 一级毛片高清免费大全| av中文乱码字幕在线| 国产伦人伦偷精品视频| 成人亚洲精品av一区二区 | 黄色视频不卡| 国产精品秋霞免费鲁丝片| 欧美乱码精品一区二区三区| 丁香欧美五月| 精品福利永久在线观看| 久久伊人香网站| 欧美日韩黄片免| 90打野战视频偷拍视频| 成人国语在线视频| 不卡av一区二区三区| 天天影视国产精品| 在线观看免费日韩欧美大片| 两性午夜刺激爽爽歪歪视频在线观看 | 看片在线看免费视频| 美女大奶头视频| 欧美精品一区二区免费开放| 好看av亚洲va欧美ⅴa在| 男人舔女人下体高潮全视频| 99国产精品一区二区蜜桃av| 亚洲avbb在线观看| 一二三四在线观看免费中文在| 88av欧美| 女警被强在线播放| 成年版毛片免费区| 成年人免费黄色播放视频| 日本三级黄在线观看| 欧美日韩黄片免| 免费人成视频x8x8入口观看| 欧美乱码精品一区二区三区| 国产欧美日韩一区二区三| 欧美日韩视频精品一区| 精品高清国产在线一区| 国产av一区二区精品久久| 免费不卡黄色视频| 精品国产超薄肉色丝袜足j| 久久中文字幕一级| 最近最新中文字幕大全电影3 | 19禁男女啪啪无遮挡网站| 一级片'在线观看视频| 悠悠久久av| 国产精品久久久久久人妻精品电影| 一级毛片精品| 在线观看一区二区三区| 国产成人一区二区三区免费视频网站| 午夜激情av网站| 天天添夜夜摸| 国产精品一区二区免费欧美| 欧美av亚洲av综合av国产av| 亚洲一码二码三码区别大吗| 精品午夜福利视频在线观看一区| 一进一出抽搐gif免费好疼 | 日韩欧美在线二视频| 午夜两性在线视频| 91精品三级在线观看| 一级片'在线观看视频| 熟女少妇亚洲综合色aaa.| 国产熟女xx| 91精品三级在线观看| 亚洲国产中文字幕在线视频| av天堂在线播放| 欧美日韩视频精品一区| 一级a爱视频在线免费观看| 欧美日韩乱码在线| 丁香六月欧美| 另类亚洲欧美激情| 国产亚洲欧美98| 黑人巨大精品欧美一区二区mp4| 成人免费观看视频高清| 久久久国产一区二区| 欧美亚洲日本最大视频资源| 中文字幕人妻丝袜制服| 淫妇啪啪啪对白视频| 色在线成人网| 欧美乱妇无乱码| 国产精品乱码一区二三区的特点 | 亚洲专区字幕在线| 午夜福利,免费看| 欧美中文日本在线观看视频| 妹子高潮喷水视频| 国产成人精品无人区| 窝窝影院91人妻| 日韩 欧美 亚洲 中文字幕| 亚洲人成电影免费在线| 亚洲精品一卡2卡三卡4卡5卡| 久久香蕉国产精品| 欧美日韩亚洲国产一区二区在线观看| 妹子高潮喷水视频| 久久亚洲真实| 一级毛片高清免费大全| 国产主播在线观看一区二区| 免费在线观看黄色视频的| 国产乱人伦免费视频| 亚洲狠狠婷婷综合久久图片| 午夜久久久在线观看| 97碰自拍视频| 一边摸一边抽搐一进一小说| 91在线观看av| 男女床上黄色一级片免费看| 成年女人毛片免费观看观看9| 老汉色∧v一级毛片| 亚洲七黄色美女视频| 窝窝影院91人妻| 黄色女人牲交| 性欧美人与动物交配| 悠悠久久av| 一个人观看的视频www高清免费观看 | 高清av免费在线| 99久久人妻综合| 国产极品粉嫩免费观看在线| 国产无遮挡羞羞视频在线观看| 精品福利永久在线观看| 久久精品亚洲精品国产色婷小说| 天堂√8在线中文| 最近最新中文字幕大全电影3 | 人成视频在线观看免费观看| 好看av亚洲va欧美ⅴa在| 午夜福利免费观看在线| 亚洲在线自拍视频| 久热这里只有精品99| 国产97色在线日韩免费| 亚洲成人久久性| 人成视频在线观看免费观看| 欧美乱色亚洲激情| 天天添夜夜摸| 美女福利国产在线| 琪琪午夜伦伦电影理论片6080| 99精品在免费线老司机午夜| 日本精品一区二区三区蜜桃| 国产av一区在线观看免费| 免费少妇av软件| 桃色一区二区三区在线观看| 波多野结衣高清无吗| 欧美日韩一级在线毛片| 色哟哟哟哟哟哟| 制服诱惑二区| 999久久久精品免费观看国产| 这个男人来自地球电影免费观看| 亚洲精品一卡2卡三卡4卡5卡| 久久精品国产亚洲av高清一级| 国产熟女xx| 黄色视频,在线免费观看| 免费一级毛片在线播放高清视频 | 亚洲成人国产一区在线观看| 成熟少妇高潮喷水视频| 国产精品久久久av美女十八| 日本撒尿小便嘘嘘汇集6| 好男人电影高清在线观看| 国产极品粉嫩免费观看在线| 成人av一区二区三区在线看| 欧美日韩视频精品一区| 丰满饥渴人妻一区二区三| 免费观看人在逋| 国产激情欧美一区二区| 日本黄色日本黄色录像| 欧美中文综合在线视频| 国产男靠女视频免费网站| 亚洲人成伊人成综合网2020| 日本黄色视频三级网站网址| 亚洲一码二码三码区别大吗| 成人永久免费在线观看视频| 日本黄色日本黄色录像| www.熟女人妻精品国产| 欧美av亚洲av综合av国产av| e午夜精品久久久久久久| 欧美+亚洲+日韩+国产| 99久久国产精品久久久| 精品人妻1区二区| 成年女人毛片免费观看观看9| 亚洲第一av免费看| 美女扒开内裤让男人捅视频| 亚洲精品久久成人aⅴ小说| 欧美另类亚洲清纯唯美| 一区福利在线观看| 日韩欧美在线二视频| 午夜福利一区二区在线看| √禁漫天堂资源中文www| 在线观看免费视频网站a站| 国产精品久久电影中文字幕| 国产免费av片在线观看野外av| 亚洲精品国产色婷婷电影| 免费av毛片视频| 国产精品一区二区在线不卡| 中文字幕人妻丝袜一区二区| 亚洲欧美日韩高清在线视频| 亚洲熟妇熟女久久| 久久 成人 亚洲| 国产精品爽爽va在线观看网站 | 国产三级在线视频| 精品欧美一区二区三区在线| 男女下面进入的视频免费午夜 | 久久久久国产精品人妻aⅴ院| 极品教师在线免费播放| 国产亚洲精品一区二区www| 久久人妻福利社区极品人妻图片| 黑人操中国人逼视频| 91成年电影在线观看| 久久久久久久久中文| 国产精品自产拍在线观看55亚洲| 欧美日韩国产mv在线观看视频| 成在线人永久免费视频| 91九色精品人成在线观看| 久久久精品欧美日韩精品| 欧美日韩国产mv在线观看视频| 国产av精品麻豆| 欧美性长视频在线观看| cao死你这个sao货| 变态另类成人亚洲欧美熟女 | 国产成人精品无人区| 欧美国产精品va在线观看不卡| 日本免费a在线| www日本在线高清视频| 交换朋友夫妻互换小说| 激情在线观看视频在线高清| 高清黄色对白视频在线免费看| 亚洲国产精品合色在线| 交换朋友夫妻互换小说| 国产精品亚洲一级av第二区| 大香蕉久久成人网| 热re99久久精品国产66热6| 国产av又大| ponron亚洲| 久久久精品国产亚洲av高清涩受|