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

    結(jié)構(gòu)可靠性的概率設(shè)計(jì)法綜述

    2017-01-04 08:22:08張學(xué)軍林延鑫吳颯王也
    裝備環(huán)境工程 2016年3期
    關(guān)鍵詞:不確定性可靠性概率

    張學(xué)軍,林延鑫,吳颯,王也

    (1.海軍裝備部飛機(jī)辦公室,北京 100071;2.北京航空航天大學(xué),北京 100191)

    結(jié)構(gòu)可靠性的概率設(shè)計(jì)法綜述

    張學(xué)軍1,林延鑫2,吳颯2,王也2

    (1.海軍裝備部飛機(jī)辦公室,北京 100071;2.北京航空航天大學(xué),北京 100191)

    綜述了目前國際上結(jié)構(gòu)可靠性設(shè)計(jì)法領(lǐng)域現(xiàn)有的概率設(shè)計(jì)法,如隨機(jī)有限元法、一階及二階可靠性設(shè)計(jì)法、蒙特卡羅法及拉丁超立方采樣法、隨機(jī)過程及隨機(jī)場。討論了具有代表性的概率設(shè)計(jì)法,如多項(xiàng)式混沌擴(kuò)展法及 KL變換?;陔S機(jī)過程概念的隨機(jī)擴(kuò)展法具有良好的分析收斂特性,直接采用隨機(jī)擴(kuò)展法進(jìn)行隨機(jī)分析是一種高效的可靠性解算方法。Karhunen和Loeve基于由協(xié)方差函數(shù)導(dǎo)出的正交坐標(biāo)函數(shù),分別對連續(xù)時(shí)間過程進(jìn)行了描述,KL變換能有效減少相關(guān)數(shù)據(jù)集的維度。當(dāng)協(xié)方差函數(shù)已知時(shí),可采用 KL擴(kuò)展法表征不確定性系統(tǒng)的隨機(jī)特征;而當(dāng)結(jié)構(gòu)響應(yīng)的協(xié)方差函數(shù)未知時(shí),可采用PCE代替KL擴(kuò)展法來表示這種類型的不確定性。出于對工程系統(tǒng)性能的嚴(yán)格要求,安全及可靠性余量的進(jìn)一步縮小,加之嚴(yán)酷的市場競爭等因素,結(jié)構(gòu)可靠性的概率設(shè)計(jì)法將迅速在多學(xué)科設(shè)計(jì)領(lǐng)域中得到廣泛應(yīng)用并產(chǎn)生顯著的工程及經(jīng)濟(jì)效益,有力地推動我國裝備質(zhì)量的跨越式發(fā)展。

    結(jié)構(gòu)可靠性設(shè)計(jì);概率設(shè)計(jì)法;可靠度;不確定性

    隨著現(xiàn)代結(jié)構(gòu)對關(guān)鍵和復(fù)雜設(shè)計(jì)需求的增加,越來越需要一種能夠?qū)V泛存在于結(jié)構(gòu)系統(tǒng)載荷、材料性能、結(jié)構(gòu)尺寸、加工過程以及使用環(huán)境中的不確定性進(jìn)行精確度量的解算方法[1—2]。當(dāng)前,科研人員廣泛認(rèn)可的不確定度定義如下:根據(jù)所用到的信息,表征賦予被測量量值分散性的非負(fù)參數(shù)[3]。何謂可靠性,就是結(jié)構(gòu)設(shè)計(jì)達(dá)到設(shè)計(jì)指標(biāo),能夠保證結(jié)構(gòu)有效穩(wěn)定的工作,而不發(fā)生破壞。重新提出結(jié)構(gòu)可靠性設(shè)計(jì),并不意味著傳統(tǒng)設(shè)計(jì)不考慮可靠性,相反由于傳統(tǒng)設(shè)計(jì)忽略了影響結(jié)構(gòu)可靠性因素的隨機(jī)性,而是選取較為保守的極值與較大的安全系數(shù)作為設(shè)計(jì)依據(jù),使得結(jié)構(gòu)可靠性過高,反而影響了結(jié)構(gòu)性能的發(fā)揮,與現(xiàn)今高效率承載結(jié)構(gòu)的設(shè)計(jì)思想相違背,也造成了物力上的浪費(fèi)。因此基于概率論和數(shù)理統(tǒng)計(jì)的新結(jié)構(gòu)可靠性設(shè)計(jì)方法正越來越受到設(shè)計(jì)人員的歡迎。這種方法認(rèn)為載荷、材料性能、結(jié)構(gòu)尺寸服從一定分布規(guī)律,依據(jù)此規(guī)律計(jì)算出結(jié)構(gòu)的可靠度與設(shè)計(jì)要求的可靠度進(jìn)行比較,從而定量地表達(dá)結(jié)構(gòu)的可靠性指標(biāo)是否滿足設(shè)計(jì)要求。結(jié)構(gòu)可靠性的概率設(shè)計(jì)法不僅可以為穩(wěn)健設(shè)計(jì)提供初始設(shè)計(jì)指南,而且還可以用來辨識不確定性對結(jié)構(gòu)系統(tǒng)的哪些性能有重要影響,或者可以為增加結(jié)構(gòu)系統(tǒng)安全性和系統(tǒng)效能提出有效的研究、測試及質(zhì)量控制方法[4—14]。

    1 結(jié)構(gòu)可靠性設(shè)計(jì)法

    結(jié)構(gòu)可靠性設(shè)計(jì)法按其提出時(shí)基于的表述方式可以分為兩類:非概率設(shè)計(jì)法和概率設(shè)計(jì)法,如圖1所示。

    圖1 結(jié)構(gòu)可靠性設(shè)計(jì)的分類Fig.1 Classification of structural reliability design

    非概率方法是用于處理關(guān)于參數(shù)真值的非精確知識。一種典型的非概率方法是區(qū)間分析,該方法把區(qū)間定義為一種新類型的數(shù),即“區(qū)間數(shù)”,任一不確定參數(shù)可能的變動范圍通過一區(qū)間表示,只需知道參數(shù)的上、下界[15]。另一種典型的是模糊理論,它描述了對象外延的不確定性和模糊性,并建立了有效處理不精確問題的工具[16]??赡苄岳碚撌菑囊粋€(gè)命題抽象出變量的分布,即自然語言命題的形式化[17]。證據(jù)理論則將不同證據(jù)上的互補(bǔ)和冗余信息依據(jù)某種優(yōu)化準(zhǔn)則組合起來,產(chǎn)生對所研究事件的一致性描述[18]。

    概率方法的理論基礎(chǔ)是概率密度函數(shù)信息,主要采用隨機(jī)變量、隨機(jī)過程和隨機(jī)場來描述系統(tǒng)的不確定性。隨機(jī)有限元方法能夠考慮材料性能的分散性、載荷的隨機(jī)性以及其他不確定性特征,適合于不確定性和隨機(jī)性模型的求解[19]。一階、二階可靠性方法在設(shè)計(jì)驗(yàn)算點(diǎn)處對極限狀態(tài)函數(shù)進(jìn)行一、二階泰勒級數(shù)展開(Taylor Series Expansion,TSE),以實(shí)現(xiàn)可靠度的解算[20]。蒙特卡羅模擬方法(Monte Carlo Simulation,MCS)通過隨機(jī)模擬和統(tǒng)計(jì)試驗(yàn)來求解結(jié)構(gòu)可靠度數(shù)值。拉丁超立方采樣法(Latin Hypercube Sampling,LHS)是一種多維分層采樣方法(Stratified Sampling,SS),其最大優(yōu)勢就在于任何大小的抽樣數(shù)目都能很容易地產(chǎn)生[21]。隨機(jī)場和隨機(jī)過程用于準(zhǔn)確估計(jì)和描述隨機(jī)空間和時(shí)間變化數(shù)據(jù)的隨機(jī)特性,對于更加真實(shí)地進(jìn)行不確定性系統(tǒng)的仿真分析顯得極為重要,其主要根據(jù)不同位置之間的相關(guān)性建立變量空間變異性的數(shù)學(xué)模型[22]。結(jié)構(gòu)可靠性設(shè)計(jì)方法主要的設(shè)計(jì)流程如圖2所示。

    圖2 結(jié)構(gòu)可靠性設(shè)計(jì)流程Fig.2 Structural reliability design flow

    2 概率設(shè)計(jì)法

    如果整個(gè)結(jié)構(gòu)或結(jié)構(gòu)的一部分超過了某一特定狀態(tài),則這個(gè)結(jié)構(gòu)或結(jié)構(gòu)的部分構(gòu)件就無法正常運(yùn)行,把這個(gè)特定的狀態(tài)稱為極限狀態(tài)。如果結(jié)構(gòu)的極限狀態(tài)超過了所要求的值,則我們認(rèn)為這個(gè)結(jié)構(gòu)是不可靠的[23—27]。對于大多數(shù)結(jié)構(gòu)來說,極限狀態(tài)可以分為兩類:

    1)最終極限狀態(tài)。最終極限狀態(tài)與部分結(jié)構(gòu)或整體結(jié)構(gòu)的失效相關(guān)。最常見的最終極限狀態(tài)有:腐蝕、疲勞、老化、起火、塑形變形、結(jié)構(gòu)連續(xù)倒塌、斷裂等。當(dāng)超過極限狀態(tài)時(shí),會導(dǎo)致人員傷亡及重大經(jīng)濟(jì)損失等風(fēng)險(xiǎn),因此極限狀態(tài)應(yīng)該具有非常低的發(fā)生概率。

    2)使用極限狀態(tài)。使用極限狀態(tài)與結(jié)構(gòu)的正常使用功能性故障相關(guān)。使用極限狀態(tài)包括過度變形、過度振動、排水和泄露、局部變形等。由于使用極限狀態(tài)的風(fēng)險(xiǎn)要小于最終極限狀態(tài),因此允許使用極限狀態(tài)有更高的失效概率[28—29]。當(dāng)然,人們也不會使用那些出現(xiàn)過度變形或振動等問題的結(jié)構(gòu)。

    一般來說,極限狀態(tài)表示結(jié)構(gòu)抗力和結(jié)構(gòu)載荷之間的安全邊界??煽慷戎笖?shù)表示為可靠度邊界的均值到失效曲面g(·)=0的距離。一維問題可靠度指數(shù)的幾何示意如圖 3所示??煽慷戎笖?shù)作為評價(jià)結(jié)構(gòu)可靠性的一個(gè)良好技術(shù)指標(biāo),其基本原理表示為位置測度μg到極限狀態(tài)曲面的距離??煽慷戎笖?shù)以不確定性比例參數(shù)σg作為距離度量的單位。圖3中的陰影部分的面積表示為結(jié)構(gòu)的失效概率。

    圖3 極限狀態(tài)g(·)的概率密度Fig.3 Probability density for limit-stateg(·)

    在某些情況下,若結(jié)構(gòu)抗力R和結(jié)構(gòu)載荷S均服從正態(tài)分布,且彼此不相關(guān),由于g(·)=0為R和S的線性函數(shù),因此極限狀態(tài)函數(shù)也為正態(tài)分布。

    3 典型的概率設(shè)計(jì)法

    解算多維問題的失效概率時(shí),往往會導(dǎo)致維度災(zāi)難問題,為此研究人員提出了許多可靠性積分問題的數(shù)值解法[30—32]。TSE 常用于極限狀態(tài)函數(shù)g(x)=0的線性化,在這種方法中,采用一階或二階TSE計(jì)算結(jié)構(gòu)系統(tǒng)的可靠性,這兩種方法分別稱為一階二次矩法(First-order Second Moment,F(xiàn)OSM)和二階二次矩法(Second-order Second Moment,SOSM)。FOSM又稱為均值一階二次矩法(Mean Value First Order Second Moment Method,MVFOSM),是一種在均值點(diǎn)處的點(diǎn)擴(kuò)展法,在 FOSM中,二階矩是最高階數(shù)的統(tǒng)計(jì)結(jié)果。盡管 FOSM比較簡單,易于實(shí)現(xiàn),但已經(jīng)證明在解算低失效概率問題(Pf<10-5)或高度非線性問題的可靠性時(shí),其解算精度并不高。相較而言,在 SOSM中,由于增加了二階項(xiàng),使得解算規(guī)模顯著增加。然而,其解算精度的提高卻是非常有限的。

    可靠度指數(shù)法從本質(zhì)上說,等價(jià)于從標(biāo)準(zhǔn)正態(tài)空間中的結(jié)構(gòu)響應(yīng)面(極限狀態(tài)近似曲面)上搜索一個(gè)最優(yōu)點(diǎn),使得從原點(diǎn)到這個(gè)點(diǎn)的距離為最短的一個(gè)數(shù)學(xué)優(yōu)化問題。Hasofer和Lind從地理學(xué)的角度對可靠度指數(shù)進(jìn)行了解釋,并通過引入Hasofer-Lind(HL)變換對FOSM進(jìn)行了改進(jìn)。在變換過程中,設(shè)計(jì)向量X被變換為標(biāo)準(zhǔn)獨(dú)立高斯向量U。由于采用了HL變換,并且設(shè)計(jì)向量在U空間具有旋轉(zhuǎn)對稱性,U空間中的設(shè)計(jì)點(diǎn)表示為最大概率密度點(diǎn)或最大似然點(diǎn)。由于這個(gè)設(shè)計(jì)點(diǎn)對名義失效概率Pf= Φ(-β)來說,是最重要的影響因子,因此稱為設(shè)計(jì)驗(yàn)算點(diǎn)(Most Probable Point,MPP)。

    不同的近似響應(yīng)面模型g(U)=0對應(yīng)于不同的失效概率解算方法。如果響應(yīng)面模型是由于 MPP處的一階函數(shù)逼近來表示,則稱為一階可靠性方法(First-order Reliability Method,F(xiàn)ORM),如果響應(yīng)面模型是由MPP處的二階函數(shù)逼近來表示,則稱為二階可靠性方法(Second-order Reliability Method,SORM)。此外,如果響應(yīng)面模型是由MPP處的更高階函數(shù)逼近來表示,則稱為高階可靠性方法(Higher-order Reliability Method,HORM)。從發(fā)展歷程來看,HL變換法通常稱為FORM,或者也稱為高級一階二次矩法(advanced FOSM)或擴(kuò)展一次二階矩法(extended FOSM)。由于概率分布不再由一階和二階矩來共同近似表示[33],因此,稱其為FORM要比HL方法更加普遍。

    在FORM中,極限狀態(tài)函數(shù)用MPP處的正切平面函數(shù)來逼近表示。近似 FORM主要用于確定失效概率的邊界,如果最大失效點(diǎn)處極限狀態(tài)函數(shù)的近似模型精確,則求得失效概率的邊界是準(zhǔn)確的,否則這種方法會產(chǎn)生很大的誤差。當(dāng)失效曲面高度非線性時(shí),F(xiàn)ORM的解算結(jié)果很不準(zhǔn)確,此時(shí),采用 FORM解算失效概率時(shí),會導(dǎo)致解算結(jié)果出現(xiàn)震蕩或?qū)е率諗康藉e(cuò)誤的解算結(jié)果[34]。

    4 擴(kuò)展的概率設(shè)計(jì)法

    基于隨機(jī)過程概念的隨機(jī)擴(kuò)展法具有良好的分析收斂特性,因而,直接采用隨機(jī)擴(kuò)展法進(jìn)行隨機(jī)分析是一種高效的可靠性解算方法[35]。隨機(jī)擴(kuò)展的目的是通過引入一系列旨在表征隨機(jī)系統(tǒng)特性的多項(xiàng)式來更好的描述系統(tǒng)的不確定性。

    4.1 多項(xiàng)式混沌擴(kuò)展法

    自從Ghanem和Spanos[36]提出譜隨機(jī)有限元法(Spectral Stochastic Finite Element Method,SSFEM)以來,多項(xiàng)式混沌擴(kuò)展法(Polynomial Chaos Expansion,PCE)已經(jīng)在包括結(jié)構(gòu)響應(yīng)等各個(gè)工程領(lǐng)域的不確定性分析中得到應(yīng)用。PCE表示為隨機(jī)變量的正交多項(xiàng)式。在SSFEM中,所采用的隨機(jī)變量通常稱為標(biāo)準(zhǔn)正態(tài)隨機(jī)變量,而所用的多項(xiàng)式一般為Hermite多項(xiàng)式。PCE依均方值收斂,任意階數(shù)的 PCE由正交多項(xiàng)式組成,這一性質(zhì)可以簡化統(tǒng)計(jì)矩的解算過程。

    Tatang[37]提出了一種概率配點(diǎn)法(Probabilistic Collocation Method,PCM),主要用于將隨機(jī)系統(tǒng)響應(yīng)映射到PCE中。Tatang等人將每個(gè)配點(diǎn)的Delta函數(shù)作為Galerkin法中的試函數(shù),并根據(jù)所選配點(diǎn)位置的模型輸出響應(yīng),確定PCE的系數(shù),即PCE的根。Isukapalli[38]指出 PCM在求解大規(guī)模問題時(shí)存在局限性,并基于模型輸出響應(yīng)對模型輸入變量的偏導(dǎo)數(shù),提出了一種隨機(jī)響應(yīng)面法,在隨機(jī)響應(yīng)面法中采用FORTRAN程序庫中的ADIFOR函數(shù)解算模型輸出響應(yīng)的偏導(dǎo)數(shù)。

    最近,研究人員應(yīng)用 PCE解算屈曲特征值問題,并基于MCS確定了PCE的系數(shù)。Xiu[39—40]等采用Askey策略擴(kuò)展了PCE,使其能夠用于表示各種不同的分布函數(shù)。Askey策略主要用于對超幾何正交多項(xiàng)式進(jìn)行分類,并表明他們之間的極限轉(zhuǎn)換關(guān)系。例如:Laguerre多項(xiàng)式可以由Jacobi多項(xiàng)式得到,也可以用來產(chǎn)生Hermite多項(xiàng)式。

    每種方法都有它的局限性,在PCM中,當(dāng)PCE具有多個(gè)自由度時(shí),所需的配點(diǎn)數(shù)就會呈指數(shù)增加,因而,當(dāng)有大量配點(diǎn)沒有被采樣時(shí),用于確定PCE系數(shù)的配點(diǎn)不能保證是如圖4a所示的空間填充試驗(yàn)設(shè)計(jì)??臻g填充試驗(yàn)設(shè)計(jì)是指根據(jù)合理定義的設(shè)計(jì)準(zhǔn)則(如最大化最小距離準(zhǔn)則)選擇樣本點(diǎn),從而使得樣本點(diǎn)能夠充滿整個(gè)設(shè)計(jì)空間。

    圖4 概率配點(diǎn)法與分層采樣法設(shè)計(jì)點(diǎn)對比Fig.4 Comparison of design points of probabilistic collocation method and stratified sampling

    如果我們對圖 5中所示的概率密度函數(shù)的截尾部分感興趣,就需要重新考慮是否能夠用 PCM選擇樣本點(diǎn)數(shù)據(jù),因?yàn)镻CM所選擇的設(shè)計(jì)點(diǎn)主要集中在高概率區(qū)域。

    圖5 概率密度函數(shù)中的感興趣區(qū)域Fig.5 Regions of interest in probability density function

    根據(jù)前面的方法,可以將所使用的隨機(jī)擴(kuò)展法,包括PCE和KL擴(kuò)展法區(qū)分為如圖6所示的兩種類型,即非嵌入式隨機(jī)擴(kuò)展法和嵌入式隨機(jī)擴(kuò)展法。嵌入式隨機(jī)擴(kuò)展法在進(jìn)行系統(tǒng)分析時(shí),無法用顯式表達(dá)式描述系統(tǒng)不確定性,而是用隱式表達(dá)式來表征系統(tǒng)的隨機(jī)特性。實(shí)際上,這就意味著采用嵌入式方法時(shí),需要對結(jié)構(gòu)分析代碼進(jìn)行修改。采用非嵌入式隨機(jī)擴(kuò)展法時(shí),只需將分析代碼當(dāng)作“黑箱”進(jìn)行處理。PCE作為非嵌入式隨機(jī)擴(kuò)展法的一種類型,能夠在不干涉結(jié)構(gòu)系統(tǒng)有限元分析模型的基礎(chǔ)上構(gòu)建系統(tǒng)的響應(yīng)近似模型。因此,PCE有時(shí)也稱為隨機(jī)響應(yīng)面法。另一種類型的非嵌入式隨機(jī)擴(kuò)展法是PCM。與之相比,嵌入式隨機(jī)擴(kuò)展法采用PCE和KL擴(kuò)展直接修改有限元分析模型的剛度矩陣。SSFEM 和Galerkin有限元法均為嵌入式隨機(jī)擴(kuò)展法[41]。

    圖6 嵌入式及非嵌入式隨機(jī)擴(kuò)展法Fig.6 Intrusive and non-intrusive formulation

    當(dāng)協(xié)方差函數(shù)已知時(shí),可采用KL擴(kuò)展法表征不確定性系統(tǒng)的隨機(jī)特征;而當(dāng)結(jié)構(gòu)響應(yīng)的協(xié)方差函數(shù)未知時(shí),可采用PCE代替KL擴(kuò)展法來表示這種類型的不確定性。

    4.2 KL變換

    正交變換能有效減少相關(guān)數(shù)據(jù)集的維度。對于隨機(jī)參數(shù)可變性的正交變換又稱為 KL變換、Hotelling變換或主成分分析。Karhunen和 Loeve基于由協(xié)方差函數(shù)導(dǎo)出的正交坐標(biāo)函數(shù),分別對連續(xù)時(shí)間過程進(jìn)行了描述。Hotelling提出采用主成分法分析多個(gè)隨機(jī)變量之間的結(jié)構(gòu)相關(guān)性[42]。當(dāng)對一個(gè)連續(xù)過程進(jìn)行有限個(gè)數(shù)據(jù)點(diǎn)的周期采樣時(shí),其自協(xié)方差矩陣的主成分分析結(jié)果與KL變換結(jié)果相同。從本質(zhì)上說,對于均值為0的向量,其主成分分析和KL變換的解算結(jié)果是相同的。

    5 工程實(shí)例

    在航空發(fā)動機(jī)適航要求中,將原發(fā)失效能夠引起發(fā)動機(jī)危害性影響的部件定義為發(fā)動機(jī)壽命限制件(Engine Life Limited Part,LLP),如旋轉(zhuǎn)輪盤、大型旋轉(zhuǎn)封嚴(yán)裝置等,設(shè)計(jì)中主要通過降低 LLP的失效概率來提高整機(jī)的安全性[43—44]。目前國內(nèi)普遍采用基于安全壽命理念的部件壽命管理方法來防止LLP的失效在一定程度上提高了航空發(fā)動機(jī)安全性水平[45]。發(fā)動機(jī)結(jié)構(gòu)完整性大綱MIL- HDBK-1783B中要求除了考慮結(jié)構(gòu)安全性和可靠性,還需要考慮壽命周期費(fèi)用問題,傳統(tǒng)的安全壽命管理方法將某一存活率下(通常為 99.87%或 99.99%)零部件萌生規(guī)定長度裂紋所需的循環(huán)數(shù)定義為部件的批準(zhǔn)壽命,通過限制使用壽命的方法來保障LLP的安全,取得了良好的效果。在處理載荷、工作環(huán)境和材料等參數(shù)的隨機(jī)性對LLP安全性的影響方面,傳統(tǒng)方法存在一些不足。為此,國外發(fā)動機(jī)制造商在型號研制中提出結(jié)構(gòu)概率可靠度計(jì)算的部件壽命管理方法進(jìn)一步降低 LLP失效導(dǎo)致的危險(xiǎn),在民用領(lǐng)域美國聯(lián)邦航空局(Federal Aviation Administration,F(xiàn)AA)在其適航規(guī)章中提出在發(fā)動機(jī)的聯(lián)合定義階段通過系統(tǒng)安全性分析確定LLP后,必須通過可靠度計(jì)算表明LLP在預(yù)期使用壽命期內(nèi)的失效概率風(fēng)險(xiǎn)每飛行小時(shí)低于 10-8,發(fā)動機(jī)才能獲得最終的型號合格證[46]。在軍用領(lǐng)域中,2002年美國頒布的 MIL-HDBK-1783B[47]中引入了概率裕度設(shè)計(jì)概念,并要求發(fā)動機(jī)內(nèi)部和安裝于外部的任何部件由于高循環(huán)疲勞導(dǎo)致失效的概率應(yīng)每飛行小時(shí)低于10-7。另外,美國在“綜合高性能渦輪發(fā)動機(jī)技術(shù)計(jì)劃(Integrated High Performance Turbine Engine Technology,IHPTET)”研究中,提出了概率設(shè)計(jì)系統(tǒng)(PDS)的新概念[48]。綜上所述,對LLP在使用壽命期內(nèi)的失效進(jìn)行概率可靠度計(jì)算成為航空發(fā)動機(jī)型號適航取證過程以及保證其結(jié)構(gòu)完整性的關(guān)鍵技術(shù)和重要實(shí)施步驟之一。

    目前開展發(fā)動機(jī) LLP結(jié)構(gòu)概率可靠性計(jì)算難點(diǎn)在于高維、小失效概率事件(FAR33.75條款要求導(dǎo)致危險(xiǎn)影響事件的 LLP發(fā)生概率每飛行小時(shí)小于 10-8)的定量要求驗(yàn)證問題。概率設(shè)計(jì)法的飛速發(fā)展及其在關(guān)鍵領(lǐng)域的廣泛應(yīng)用,為破解上述技術(shù)難點(diǎn)奠定了堅(jiān)實(shí)的根基。

    6 結(jié)語

    綜述了目前國際上在結(jié)構(gòu)可靠性設(shè)計(jì)法領(lǐng)域先進(jìn)高效的概率設(shè)計(jì)法,討論了當(dāng)前在各個(gè)關(guān)鍵核心領(lǐng)域所廣泛應(yīng)用的結(jié)構(gòu)可靠性設(shè)計(jì)法。出于對工程系統(tǒng)性能的嚴(yán)格要求,安全及可靠性余量的進(jìn)一步縮小,以及嚴(yán)酷的市場競爭等因素,可靠性分析方法迅速在多學(xué)科設(shè)計(jì)領(lǐng)域得到應(yīng)用。包含不確定性的結(jié)構(gòu)設(shè)計(jì)中,采用確定性方法可能要比考慮不確定性的概率方法導(dǎo)致出現(xiàn)更大的失效概率。這主要是由于在確定性方法中要求精確滿足設(shè)計(jì)需求,但是參數(shù)的任意擾動均可能引起系統(tǒng)響應(yīng)違反約束。在進(jìn)行非傳統(tǒng)結(jié)構(gòu)的設(shè)計(jì)時(shí),往往缺乏相關(guān)數(shù)據(jù)和足夠的先驗(yàn)知識,而對不確定性的正確認(rèn)識對于安全、高效的決策至關(guān)重要。

    通過分析可以看出,任何方法都是在一定假設(shè)的基礎(chǔ)上提出的,當(dāng)某種方法的適用條件成立時(shí),該方法可以給出較好的計(jì)算結(jié)果;當(dāng)該方法的適用條件違背時(shí),該方法的計(jì)算結(jié)果可能會有較大誤差。因此,如何有效地改進(jìn)現(xiàn)有方法,提出新的方法,使其具有更廣泛的適用范圍是未來結(jié)構(gòu)可靠性設(shè)計(jì)法的主要研究方向。

    中國制造 2025的基本方針之一——“質(zhì)量為先”明確指出:“堅(jiān)持把質(zhì)量作為建設(shè)制造強(qiáng)國的生命線,強(qiáng)化企業(yè)質(zhì)量主體責(zé)任,加強(qiáng)質(zhì)量技術(shù)攻關(guān)、自主品牌培育。建設(shè)法規(guī)標(biāo)準(zhǔn)體系、質(zhì)量監(jiān)管體系、先進(jìn)質(zhì)量文化,營造誠信經(jīng)營的市場環(huán)境,走以質(zhì)取勝的發(fā)展道路”。結(jié)構(gòu)可靠性設(shè)計(jì)法及其應(yīng)用性研究,尤其是概率設(shè)計(jì)法正是提升產(chǎn)品質(zhì)量、工程質(zhì)量的強(qiáng)大引擎和關(guān)鍵保證。

    [1] 程家軍, 李春枝, 陳穎. 艦載環(huán)境振動搖擺特性測試分析[J]. 裝備環(huán)境工程, 2015, 12(1): 114—119. CHENG Jia-jun, LI Chun-zhi, CHEN Ying. Testing Analysis of Vibration and Sway Characteristics in Shipboard Environment[J]. Equipment Environmental Engineering, 2015, 12 (1): 114—119.

    [2] 李勘. 裝備研制項(xiàng)目可靠性風(fēng)險(xiǎn)預(yù)警技術(shù)研究[J]. 裝備環(huán)境工程, 2014, 11(4): 125—130. LI Kan. Research on Pre-warning Technology for Reliability Risk of Equipment Research Project[J]. Equipment Environmental Engineering, 2014, 11(4): 125—130.

    [3] 倪育才. 實(shí)用測量不確定度評定[M]. 第4版. 北京: 中國質(zhì)檢出版社, 2014: 7—50. NI Yu-cai. Evaluation of the Uncertainty of Practical Measurement[M]. Fourth Edition. Beijing: China Quality Inspection Press, 2014: 7—50.

    [4] 范義龍,吳颯,李琦. 基于PLC的冷卻系統(tǒng)自整定模糊控制研究[J]. 環(huán)境技術(shù), 2014(6): 66—71. FAN Yi-long, WU Sa, LI Qi. Research of the PLC Based Parameter Self-tuning Fuzzy Control for Cooling System [J]. Environmental Technology, 2014(6): 66—71.

    [5] HORáLEK V. Analysis of Basic Probability Distributions, Their Properties and Use in Determining Type B Evaluation of Measurement Uncertainties[J]. Measurement, 2013, 46 (1): 16—23.

    [6] CHU-AGOR M L, MU?OZ-CARPENA R, KIKER G A, et al. Simulating the Fate of Florida Snowy Plovers with Sea-level rise: Exploring Research and Management Priorities with a Global Uncertainty and Sensitivity Analysis Perspective[J]. Ecological Modelling, 2012, 224(1): 33—47.

    [7] SCHWABE O, SHEHAB E, ERKOYUNCU J. Uncertainty Quantification Metrics for Whole Product Life Cycle Cost Estimates in Aerospace Innovation[J]. Progress in Aerospace Sciences, 2015, 77: 1—24.

    [8] DANOWSKY B P, CHRSTOS J R, KLYDE D H, et al. Evaluation of Aeroelastic Uncertainty Analysis Methods [J]. Journal of Aircraft, 2010, 47(4): 1266—1273.

    [9] MIRZAEI M, GALAVI H, FAGHIH M, et al. Model Calibration and Uncertainty Analysis of Runoff in the Zayanderood River Basin Using Generalized Likelihood Uncertainty Estimation (GLUE) Method[J]. Aqua, 2013, 62(5): 309—320.

    [10] WU Y, PAN X, KANG R, et al. Multi-parameters Uncertainty Analysis of Logistic Support Process Based on GERT[J]. Journal of Systems Engineering & Electronics, 2014, 25(6): 1011—1019.

    [11] ZHAO Y, FAN F, WANG J, et al. Uncertainty Analysis for Bulk Power Systems Reliability Evaluation Using Taylor Series and Nonparametric Probability Density Estimation[J]. International Journal of Electrical Power & Energy Systems, 2015, 64: 804—814.

    [12] MOON H J. Uncertainty Analysis in Mould Spore Transportation and Its Application in an Existing Building[J]. Indoor & Built Environment, 2010, 19(19): 355—365.

    [13] CAMACHO R A, MARTIN J L, DIAZ-RAMIREZ J, et al. Uncertainty Analysis of Estuarine Hydrodynamic Models: An Evaluation of Input Data Uncertainty in the Weeks Bay Estuary, Alabama[J]. Applied Ocean Research, 2014, 47 (2): 138—153.

    [14] SHI F, LONG H, ZHAN M, et al. Uncertainty Analysis on Process Responses of Conventional Spinning Using Finite Element Method[J]. Structural & Multidisciplinary Optimization, 2014, 49(5): 839—850.

    [15] 朱慶鵬, 孫宇, 陳浩, 等. 基于區(qū)間層次分析法的高速送料機(jī)可靠性綜合分配方法[J]. 鍛壓技術(shù), 2015(4): 131—136. ZHU Qing-peng, SUN Yu, CHEN Hao, et. Comprehensive Reliability Allocation Method for High-speed Feeder Based on Interval Analytic Hierarchy Process[J]. Forging & Stamping Technology, 2015(4): 131—136.

    [16] 趙遠(yuǎn), 焦健, 趙廷弟. 基于模糊理論的風(fēng)險(xiǎn)評價(jià)方法[J].系統(tǒng)工程與電子技術(shù), 2015(8): 1825—1831. ZHAO Yuan, JIAO Jian, ZHAO Ting-di. Risk AssessmentMethod Based on Fuzzy Logic[J]. Systems Engineering and Electronics, 2015 (08): 1825—1831.

    [17] 嚴(yán)巍, 趙川. 可能性理論研究[J]. 數(shù)字通信, 2014(5): 24—27. YAN Wei, ZHAO Chuan. Brief Analysis of the Possibility Theory[J]. Digital Communication, 2014(5): 24—27.

    [18] 唐和生, 鄧立新, 胡長遠(yuǎn), 等. 基于證據(jù)理論和微分演化的結(jié)構(gòu)不確定分析[J]. 同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014(3): 325—330. TANG He-sheng, DENG Li-xin, HU Chang-yuan, et al. Evidence Theory and Differential Evolution for Uncertainty Quantification of Structures[J]. Journal of Tongji University(Natural Science), 2014(3): 325—330.

    [19] 洪興福, 胡祥濤, 鄒滿玲. 基于隨機(jī)有限元方法的地面雷達(dá)支耳連接結(jié)構(gòu)可靠性評估[J]. 機(jī)械設(shè)計(jì)與研究, 2014(6): 59—61. HONG Xing-fu, HU Xiang-tao, ZOU Man-ling. Reliability Assessment for Radar Supporting Structure Based on Stochastic Finite Element Method[J]. Mechanical Design and Research, 2014(6): 59—61.

    [20] 郭彪, 顧德華, 董玉革, 等. 可靠性分析的 FORM 和SORM 組合法[J]. 合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014(1): 30—33. GUO Biao, GU De-hua, DONG Yu-ji, et al. Combination of FORM and SORM for Reliability Analysis[J]. Journal of Hefei University of Technology (Natural Science), 2014(1): 30—33.

    [21] 伍國軍, 陳衛(wèi)忠, 譚賢君, 等. 基于拉丁超立方抽樣的有限元可靠度程序開發(fā)及應(yīng)用[J]. 巖土力學(xué), 2015 (2): 550—554. WU Guo-jun, CHEN Wei-zhong, TAN Xian-jun, et al. Program Development of Finite Element Reliability Method and Its Application Based on Latin Hypercube Sampling[J]. Rock and Soil Mechanics, 2015(2): 550—554.

    [22] 王濤, 周國慶, 陰琪翔, 等. 巖土參數(shù)隨機(jī)場離散的三角形單元局部平均法[J]. 巖土力學(xué), 2014(5): 1482—1488. WANG Tao, ZHOU Guo-qing, YIN Qi-xiang, et al. Local Average Method of Triangular Elements for Discretization of Random Field of Geotechnical Parameters[J]. Rock and Soil Mechanics, 2014(5): 1482—1488.

    [23] DAERYONG P, ROESNER L A. Evaluation of Pollutant Loads from Stormwater BMPs to Receiving Water Using Load Frequency Curves with Uncertainty Analysis[J]. Water Research, 2012, 46(20): 6881—6890.

    [24] MANDAL T K, GHOSH S, PUJARI N N. Seismic Fragility Analysis of a Typical Indian PHWR Containment: Comparison of Fragility Models[J]. Structural Safety, 2015, 58: 11—19.

    [25] BERNUZZI C, GABBIANELLI G, GOBETTI A, et al. Beam Design for Steel Storage Racks[J]. Journal of Constructional Steel Research, 2016, 116(1): 156—172.

    [26] BECK A T, áVILA DA SILVA C R J. Strategies for Finding the Design Point under Bounded Random Variables[J]. Structural Safety, 2016, 58: 79—93.

    [27] FREEMAN R R. Simplified Uncertainty Analysis of Layer of Protection Analysis Results[J]. Process Safety Progress, 2013, 32(4): 351—360.

    [28] 李鵬, 李傳日, 李桃. 基于可靠性分配的整機(jī)加速因子計(jì)算方法研究[J]. 裝備環(huán)境工程, 2015, 12(3): 53—57. LI Peng, LI Chuan-ri, LI Tao. Research on the Calculation Method of Overall Acceleration Factor Based on Reliability Allocation[J]. Equipment Environmental Engineering, 2015, 12(3): 53—57.

    [29] 向樹紅. “航天器環(huán)境新技術(shù)研究”專題序言[J]. 裝備環(huán)境工程, 2015, 12(3): 3—4. XIANG Shu-hong. The Introduction of the New Technology Research of Spacecraft Environment[J]. Equipment Environmental Engineering, 2015, 12(3): 3—4.

    [30] 楊云峰, 湯鵬杰. 服務(wù)器集群網(wǎng)絡(luò)中糾纏疑似攻擊檢測技術(shù)研究[J]. 計(jì)算機(jī)仿真, 2015, 32(5): 315—317. TANG Peng-jie, YANG Yun-feng. Research on Detection Technology for Entanglement Suspected Attack in Server Cluster Network[J]. Computer Simulation, 2015, 32(5): 315—317.

    [31] 王紅, 王希誠, 李克秋. 混合Kriging代理模型的高維參數(shù)估計(jì)優(yōu)化算法[J]. 大連理工大學(xué)學(xué)報(bào). 2015(2): 215—222. WANG Hong, WANG Xi-cheng, LI Ke-qiu. Hybrid Kriging Surrogate Model Optimization Algorithm for High-dimension Parameter Estimation[J]. Journal of Dalian University of Technology, 2015(2): 215—222.

    [32] 蘇盈盈, 李太福, 易軍, 等. 基于KICA子空間虛假鄰點(diǎn)判別的軟傳感器變量選擇方法[J]. 機(jī)械工程學(xué)報(bào), 2015(4): 15—21. SU Ying-ying, LI Tai-fu, YI Jun, et al. Variable Selection for Nonlinear Soft Sensor Based on False Nearest Neighbors in KICA Space [J]. Journal of Mechanical Engineering, 2015(4): 15—21.

    [33] MELCHERS R E. Structural Reliability Analysis and Prediction[M]. UK: Ellis Horwood Limited, 1987.

    [34] 仲偉君, 李德勝, 崔亮. 一種制導(dǎo)彈藥電子部件可靠性預(yù)計(jì)方法[J]. 裝備環(huán)境工程, 2014, 11(2): 103—105. ZHONG Wei-jun, LI De-sheng, CUI Liang. A Method for the Reliability Prediction of Electronic Components in Guided Munitions[J]. Equipment Environmental Engineering, 2014, 11 (2): 103—105.

    [35] CAMERON R H, MARTIN W T. The Orthogonal Development of Nonlinear Functionals in Series of Fourier-Hermite Functionals[J]. Annals of Mathematics, 1947, 48: 385—392.

    [36] GHANEM R G, SPANOS P D. Stochastic Finite Elements: A Spectral Approach [M]. NY: Springer, 1991.

    [37] TATANG M A. Direct Incorporation of Uncertainty in Chemical and Environmental Engineering Systems [D]. Cambrigde, MA: Massachusetts Institute of Technology, 1995.

    [38] ISUKAPALLI S S. Uncertainty Analysis of Transport-Transformation Models[D]. New Brunswick, NJ: the State University of New Jersey, 1999.

    [39] XIU D, KARNIADAKIS G. The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations[J]. SIAM Journal of Scientific Computing, 2002, 24(2): 619—644.

    [40] XIU D, LUCOR D, SU C, et al. Stochastic Modeling of Flow-structure Interactions Using Generalized Polynomial Chaos[J]. Journal of Fluids Engineering, Transactions of the ASME, 2002, 24(51): 51—59.

    [41] BABUSKA I, TEMPONE R, ZOURARIS G E. Galerkin Finite Element Approximations of Stochastic Elliptic Differential Equations[J]. SIAM Journal on Numerical Analysis, 2004, 42: 800—825.

    [42] HOTELLING H. Analysis of a Complex of Statistical Variables into Principal Components[J]. Journal of Educational Psychology, 1933, 24: 498—520.

    [43] U S Department of Transportation, Federal Aviation Administration. Advisory Circular 33. 70-1: Guidance Material for Aircraft Engine-Life-Limited Parts Requirements [R]. Washington DC: FAA, AC33 70-1, 2009: 30—50.

    [44] VITTAL S, HAJELA P, JOSHI A. Review of Approaches to Gas Turbine Life Management[R]. AIAA-2004-4372, 2004.

    [45] 王衛(wèi)國. 輪盤低循環(huán)疲勞壽命預(yù)測模型和試驗(yàn)評估方法研究[D]. 南京: 南京航空航天大學(xué), 2006. WANG Wei-guo. Disc LCF Life Prediction Models and Experiment Assessment Methodologies[D]. Nanjing: Nanjing University of Aeronautics & Astronautics, 2006.

    [46] 丁水汀, 張弓, 蔚奪魁, 等. 航空發(fā)動機(jī)適航概率風(fēng)險(xiǎn)評估方法研究綜述[J]. 航空動力學(xué)報(bào), 2011, 26(7): 1142—1143. DING Shui-ting, ZHANG Gong, WEI Duo-kui, et al. Review of Probabilistic Risk Assessment on Aero-engine Airworthiness[J]. Journal of Aerospace Power, 2011, 26 (7): 1142—1143.

    [47] MIL-HDBK-1783B, Engine Structural Integrity Program [S].

    [48] MCCLUNG R C, LEVERANT G R, ENRIGHT M P. Turbine Rotor Material Design-Phase II[R]. Washington, DC: FAA Grant 99-G-016, 2008: 50—55.

    Review on Probabilistic Design Method for Structural Reliability

    ZHANG Xue-jun1, LIN Yan-xin2, WU Sa2, WANG Ye2

    (1. Office of the Navy Equipment Department, Beijing 100071, China; 2.Beihang University, Beijing 100191, China)

    The existing probabilistic design methods in the field of current international structure reliability design method are reviewed, such as stochastic finite element method, first order and second order reliability design method, Monte Carlo method and Latin hypercube sampling method, random processes and random fields. The representative probabilistic design methods are discussed, such as polynomial chaos expansion method and KL transform. Stochastic expansion method based on the concept of stochastic process has a good analytic convergence property, and stochastic analysis is an efficient method for reliability calculation. Karhunen and Loeve are derived from the covariance function based on the orthogonal coordinate function, the continuous time process is described respectively, and KL transform can effectively reduce the dimension of the data set. When the covariance function is known, the KL expansion of the stochastic characteristics of the method can be used to characterize the uncertainty system, and when the structural response of the covariance function is unknown, the PCE can instead of KL expansion method to represent this type of uncertainty. Considering the strict requirements of the performance of engineering system, safety and reliability margins are further narrowed, and coupled with the harsh market competition and other factors, the probability design method of structure reliability will be widely applied in the field of the multidisciplinary design and pro-duce significant engineering and economic benefits, effectively promoting the great leap forward development of China's equipment quality.

    structural reliability design; probability design method; reliability; uncertainty

    10.7643/ issn.1672-9242.2016.03.027

    TJ01;V416

    A

    1672-9242(2016)03-0161-08

    2016-02-02;

    2016-04-02

    Received:2016-02-02;Revised:2016-04-02

    張學(xué)軍(1964—),男,黑龍江集賢人,高級工程師,主要研究方向?yàn)楹娇昭b備管理。

    Biography:ZHANG Xue-jun (1964—), Male, from Jixian, Heilongjiang, Senior engineer, Research focus: aeronautical equipment management.

    猜你喜歡
    不確定性可靠性概率
    法律的兩種不確定性
    法律方法(2022年2期)2022-10-20 06:41:56
    第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
    第6講 “統(tǒng)計(jì)與概率”復(fù)習(xí)精講
    概率與統(tǒng)計(jì)(一)
    概率與統(tǒng)計(jì)(二)
    可靠性管理體系創(chuàng)建與實(shí)踐
    英鎊或繼續(xù)面臨不確定性風(fēng)險(xiǎn)
    中國外匯(2019年7期)2019-07-13 05:45:04
    電子制作(2017年2期)2017-05-17 03:55:06
    具有不可測動態(tài)不確定性非線性系統(tǒng)的控制
    基于可靠性跟蹤的薄弱環(huán)節(jié)辨識方法在省級電網(wǎng)可靠性改善中的應(yīng)用研究
    電測與儀表(2015年6期)2015-04-09 12:01:18
    97热精品久久久久久| 深夜a级毛片| 亚洲欧美日韩卡通动漫| 欧美三级亚洲精品| 国产伦一二天堂av在线观看| 亚洲av免费高清在线观看| 变态另类丝袜制服| 日本-黄色视频高清免费观看| 少妇的逼水好多| 麻豆精品久久久久久蜜桃| 欧美性感艳星| 免费av毛片视频| 在线天堂最新版资源| 麻豆乱淫一区二区| 性色avwww在线观看| 高清日韩中文字幕在线| 97人妻精品一区二区三区麻豆| 国产成人freesex在线 | 精品久久久久久久人妻蜜臀av| а√天堂www在线а√下载| 精华霜和精华液先用哪个| 久久久国产成人精品二区| 长腿黑丝高跟| 啦啦啦观看免费观看视频高清| 97超视频在线观看视频| 国产高清激情床上av| 男女啪啪激烈高潮av片| 插逼视频在线观看| 永久网站在线| 尾随美女入室| 精品熟女少妇av免费看| 能在线免费观看的黄片| 久久久久久大精品| 国产成人a区在线观看| 久久久色成人| 我的女老师完整版在线观看| 搡女人真爽免费视频火全软件 | 久久久久久久久久成人| 免费黄网站久久成人精品| 51国产日韩欧美| 精品人妻熟女av久视频| 国产乱人视频| 精品久久久久久久久亚洲| 一进一出抽搐gif免费好疼| 久久久久久伊人网av| 九九久久精品国产亚洲av麻豆| 日韩强制内射视频| 精品无人区乱码1区二区| 桃色一区二区三区在线观看| 99热这里只有是精品50| 午夜福利在线观看免费完整高清在 | 最近手机中文字幕大全| 最新在线观看一区二区三区| 十八禁国产超污无遮挡网站| 全区人妻精品视频| 九九久久精品国产亚洲av麻豆| 夜夜爽天天搞| 中文字幕免费在线视频6| 色5月婷婷丁香| 亚洲国产高清在线一区二区三| 变态另类成人亚洲欧美熟女| 午夜影院日韩av| 国产一区二区在线av高清观看| 亚洲av第一区精品v没综合| 免费看av在线观看网站| 特级一级黄色大片| 日本成人三级电影网站| 菩萨蛮人人尽说江南好唐韦庄 | 身体一侧抽搐| 亚洲中文日韩欧美视频| 精品久久久久久久人妻蜜臀av| 欧美最黄视频在线播放免费| 麻豆一二三区av精品| 日韩大尺度精品在线看网址| 熟女人妻精品中文字幕| а√天堂www在线а√下载| 国产精品久久视频播放| 亚洲精品影视一区二区三区av| 国产一区二区在线av高清观看| 日韩欧美免费精品| 日韩欧美免费精品| 夜夜爽天天搞| 最近最新中文字幕大全电影3| 小说图片视频综合网站| 婷婷六月久久综合丁香| 日本a在线网址| 熟女人妻精品中文字幕| 伦精品一区二区三区| 日本免费一区二区三区高清不卡| 在线观看美女被高潮喷水网站| 亚洲一区高清亚洲精品| 欧美绝顶高潮抽搐喷水| 国产蜜桃级精品一区二区三区| 成年女人永久免费观看视频| 中文字幕精品亚洲无线码一区| 国产在视频线在精品| 国产激情偷乱视频一区二区| 国产精品久久久久久亚洲av鲁大| 亚洲精品粉嫩美女一区| 麻豆av噜噜一区二区三区| 国产色婷婷99| 人妻制服诱惑在线中文字幕| 亚洲成人中文字幕在线播放| 18禁在线播放成人免费| 老司机福利观看| 欧美人与善性xxx| 免费看a级黄色片| 久久久久久久久久黄片| 综合色丁香网| 亚洲第一电影网av| 欧美最新免费一区二区三区| 深夜精品福利| 欧美色欧美亚洲另类二区| 国产伦在线观看视频一区| 日本与韩国留学比较| 国产精品人妻久久久影院| 欧美色视频一区免费| 黄色欧美视频在线观看| 小说图片视频综合网站| 女生性感内裤真人,穿戴方法视频| 18禁在线无遮挡免费观看视频 | 亚洲精品在线观看二区| 天美传媒精品一区二区| 美女cb高潮喷水在线观看| 在线看三级毛片| 特级一级黄色大片| 3wmmmm亚洲av在线观看| 18禁裸乳无遮挡免费网站照片| 日韩大尺度精品在线看网址| 国产大屁股一区二区在线视频| 九九在线视频观看精品| 欧美bdsm另类| 欧美成人一区二区免费高清观看| 精品人妻一区二区三区麻豆 | 日韩欧美一区二区三区在线观看| 淫秽高清视频在线观看| 日本撒尿小便嘘嘘汇集6| 午夜精品国产一区二区电影 | 午夜福利视频1000在线观看| 亚洲国产精品国产精品| 久久热精品热| 国产一区亚洲一区在线观看| 乱人视频在线观看| 亚洲一区高清亚洲精品| 特级一级黄色大片| 日韩欧美一区二区三区在线观看| 亚洲五月天丁香| 亚洲av成人精品一区久久| 在线观看免费视频日本深夜| 国国产精品蜜臀av免费| 99在线人妻在线中文字幕| 99久国产av精品国产电影| 亚洲欧美日韩高清在线视频| 我的老师免费观看完整版| 床上黄色一级片| 嫩草影院入口| 亚洲不卡免费看| 夜夜看夜夜爽夜夜摸| 色综合亚洲欧美另类图片| 给我免费播放毛片高清在线观看| 日本在线视频免费播放| 波多野结衣巨乳人妻| 国产午夜福利久久久久久| 女生性感内裤真人,穿戴方法视频| 午夜精品在线福利| 国产探花在线观看一区二区| 精华霜和精华液先用哪个| 欧美日韩乱码在线| 性欧美人与动物交配| 精品福利观看| www.色视频.com| 可以在线观看的亚洲视频| 99精品在免费线老司机午夜| 精品一区二区三区人妻视频| 我要看日韩黄色一级片| 91精品国产九色| 欧美日韩综合久久久久久| 免费电影在线观看免费观看| 能在线免费观看的黄片| 欧美+日韩+精品| 男人舔女人下体高潮全视频| 久久午夜福利片| 亚洲人成网站高清观看| av在线蜜桃| 天堂√8在线中文| 国产精品一二三区在线看| 国产伦一二天堂av在线观看| 婷婷色综合大香蕉| 乱系列少妇在线播放| 亚洲无线观看免费| 久久精品国产自在天天线| 美女cb高潮喷水在线观看| 丝袜喷水一区| av福利片在线观看| 亚洲最大成人av| 国产精品一区www在线观看| 久久精品国产亚洲av涩爱 | 寂寞人妻少妇视频99o| 黄色视频,在线免费观看| 美女 人体艺术 gogo| 一卡2卡三卡四卡精品乱码亚洲| 成人特级黄色片久久久久久久| 久久久久久伊人网av| 超碰av人人做人人爽久久| aaaaa片日本免费| 亚洲av中文字字幕乱码综合| 国产午夜精品论理片| 嫩草影视91久久| 国产国拍精品亚洲av在线观看| 精品久久久久久久末码| 精品国内亚洲2022精品成人| 久久久久久久久大av| 日韩欧美精品免费久久| 国产精品99久久久久久久久| 久久国产乱子免费精品| 最近中文字幕高清免费大全6| 又黄又爽又刺激的免费视频.| av.在线天堂| 亚洲av一区综合| 国产视频内射| 人人妻人人澡人人爽人人夜夜 | 久久久午夜欧美精品| 99国产精品一区二区蜜桃av| 中文在线观看免费www的网站| 又粗又爽又猛毛片免费看| 嫩草影院入口| 亚洲人成网站在线观看播放| 搡女人真爽免费视频火全软件 | 免费人成在线观看视频色| 国产女主播在线喷水免费视频网站 | 五月玫瑰六月丁香| av女优亚洲男人天堂| 搞女人的毛片| 亚洲自偷自拍三级| 午夜精品国产一区二区电影 | 日本免费一区二区三区高清不卡| 久久久久国产网址| 亚洲第一区二区三区不卡| 小蜜桃在线观看免费完整版高清| 18禁在线播放成人免费| 国产一区二区三区在线臀色熟女| 日韩欧美三级三区| 国产一区亚洲一区在线观看| 免费看美女性在线毛片视频| 国产黄片美女视频| 婷婷六月久久综合丁香| 亚洲国产精品合色在线| 亚洲国产欧美人成| 晚上一个人看的免费电影| 两个人的视频大全免费| 久久久久久大精品| 淫妇啪啪啪对白视频| 一级a爱片免费观看的视频| 精品欧美国产一区二区三| 18+在线观看网站| 国产三级中文精品| 男插女下体视频免费在线播放| 成人一区二区视频在线观看| 成人特级黄色片久久久久久久| 欧美一区二区亚洲| av黄色大香蕉| 精华霜和精华液先用哪个| 亚洲欧美日韩高清专用| 成人亚洲欧美一区二区av| 黄色视频,在线免费观看| 免费在线观看成人毛片| 12—13女人毛片做爰片一| 乱人视频在线观看| 黄色欧美视频在线观看| 自拍偷自拍亚洲精品老妇| 美女内射精品一级片tv| 欧美潮喷喷水| 国产精品一及| 欧洲精品卡2卡3卡4卡5卡区| 国产男靠女视频免费网站| 毛片一级片免费看久久久久| 非洲黑人性xxxx精品又粗又长| 女同久久另类99精品国产91| 国产精品人妻久久久影院| 女的被弄到高潮叫床怎么办| 日本五十路高清| 午夜免费激情av| 国产精品乱码一区二三区的特点| 亚洲精品一区av在线观看| 亚洲va在线va天堂va国产| 国产色婷婷99| 色吧在线观看| 99热这里只有是精品在线观看| 日本黄色视频三级网站网址| 亚洲av熟女| 一级av片app| 久久精品人妻少妇| 毛片女人毛片| 国产亚洲精品久久久久久毛片| 女生性感内裤真人,穿戴方法视频| 久久久久精品国产欧美久久久| 97人妻精品一区二区三区麻豆| 别揉我奶头~嗯~啊~动态视频| 精品福利观看| 久久久成人免费电影| 狂野欧美白嫩少妇大欣赏| 全区人妻精品视频| 欧美成人一区二区免费高清观看| 国产欧美日韩精品一区二区| 日本 av在线| 国产69精品久久久久777片| 欧美不卡视频在线免费观看| 亚洲最大成人av| 久久久久久久久久久丰满| 五月伊人婷婷丁香| 国产视频内射| 亚洲精品日韩在线中文字幕 | 国产精品电影一区二区三区| 亚洲综合色惰| 久久精品国产亚洲av天美| 亚洲成人久久性| 插阴视频在线观看视频| 大香蕉久久网| 国产男人的电影天堂91| 男女那种视频在线观看| 亚洲成人精品中文字幕电影| 国产美女午夜福利| 免费观看在线日韩| 久久午夜福利片| 看片在线看免费视频| 国产精品一区二区三区四区久久| 色视频www国产| 精品人妻偷拍中文字幕| 亚洲七黄色美女视频| 欧美日本视频| 亚洲一级一片aⅴ在线观看| 欧美另类亚洲清纯唯美| 天天躁夜夜躁狠狠久久av| 国产亚洲精品久久久com| 波野结衣二区三区在线| 亚洲精品影视一区二区三区av| av专区在线播放| 十八禁网站免费在线| 国产69精品久久久久777片| 久久精品人妻少妇| 99热6这里只有精品| 日韩人妻高清精品专区| 在线免费观看不下载黄p国产| 国产精品综合久久久久久久免费| 国产高清有码在线观看视频| 搡老熟女国产l中国老女人| 久久精品国产亚洲av香蕉五月| 中出人妻视频一区二区| 97人妻精品一区二区三区麻豆| 热99re8久久精品国产| 人妻夜夜爽99麻豆av| 春色校园在线视频观看| av福利片在线观看| 黄色视频,在线免费观看| 午夜视频国产福利| 久久久成人免费电影| 校园春色视频在线观看| 久久人人精品亚洲av| 三级毛片av免费| h日本视频在线播放| 国产精品三级大全| 又黄又爽又刺激的免费视频.| 欧美日韩在线观看h| 一边摸一边抽搐一进一小说| 伦理电影大哥的女人| 女人十人毛片免费观看3o分钟| 一个人看视频在线观看www免费| 久久中文看片网| 老师上课跳d突然被开到最大视频| 国产av在哪里看| 热99re8久久精品国产| 99视频精品全部免费 在线| 一边摸一边抽搐一进一小说| 联通29元200g的流量卡| 天堂动漫精品| 男女那种视频在线观看| 国产av不卡久久| 久久久国产成人免费| 成人特级黄色片久久久久久久| 九九久久精品国产亚洲av麻豆| 午夜精品国产一区二区电影 | 欧美一级a爱片免费观看看| 九九久久精品国产亚洲av麻豆| 日本欧美国产在线视频| 成人一区二区视频在线观看| 亚州av有码| 欧美成人免费av一区二区三区| 国产一区二区三区av在线 | 免费观看人在逋| 国产伦在线观看视频一区| 免费电影在线观看免费观看| 国产在视频线在精品| 一进一出抽搐gif免费好疼| 午夜爱爱视频在线播放| 午夜福利高清视频| 亚洲国产色片| 嫩草影视91久久| 乱系列少妇在线播放| 国产高清视频在线播放一区| videossex国产| av天堂中文字幕网| 成人午夜高清在线视频| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久久久久黄片| 亚洲av成人精品一区久久| 成人三级黄色视频| 国产精品久久视频播放| 精品熟女少妇av免费看| 最近2019中文字幕mv第一页| 蜜臀久久99精品久久宅男| 麻豆精品久久久久久蜜桃| 高清午夜精品一区二区三区 | 亚洲欧美日韩卡通动漫| 久久久久久久久中文| 中国国产av一级| av在线老鸭窝| 久久久久久久午夜电影| 久久久久久久久中文| 美女xxoo啪啪120秒动态图| 国产精品久久久久久久久免| 国内少妇人妻偷人精品xxx网站| 一本一本综合久久| 男人舔奶头视频| 欧美xxxx黑人xx丫x性爽| 成人特级黄色片久久久久久久| 免费大片18禁| h日本视频在线播放| 国产av不卡久久| 日日摸夜夜添夜夜添小说| 色综合站精品国产| 最近最新中文字幕大全电影3| 波多野结衣高清作品| 真人做人爱边吃奶动态| 免费大片18禁| 赤兔流量卡办理| 超碰av人人做人人爽久久| 亚洲国产精品合色在线| 免费观看人在逋| 白带黄色成豆腐渣| 久久九九热精品免费| 国产爱豆传媒在线观看| 久久久国产成人精品二区| 一区福利在线观看| 欧美中文日本在线观看视频| 搡女人真爽免费视频火全软件 | 日韩一区二区视频免费看| 日韩强制内射视频| 中文资源天堂在线| 狠狠狠狠99中文字幕| 亚洲丝袜综合中文字幕| 男女做爰动态图高潮gif福利片| 精品无人区乱码1区二区| 久久热精品热| 最后的刺客免费高清国语| 高清午夜精品一区二区三区 | 麻豆精品久久久久久蜜桃| 亚洲国产精品国产精品| 亚洲av免费在线观看| 一卡2卡三卡四卡精品乱码亚洲| 精品国产三级普通话版| 午夜激情福利司机影院| www日本黄色视频网| 国产乱人偷精品视频| 国产男人的电影天堂91| 99热精品在线国产| 听说在线观看完整版免费高清| 不卡视频在线观看欧美| 亚洲va在线va天堂va国产| 12—13女人毛片做爰片一| 淫秽高清视频在线观看| 少妇熟女aⅴ在线视频| 亚洲天堂国产精品一区在线| av国产免费在线观看| 又爽又黄a免费视频| 亚洲av成人av| 亚洲av中文字字幕乱码综合| 成人午夜高清在线视频| 国产一级毛片七仙女欲春2| 1024手机看黄色片| 日韩欧美精品v在线| 国产乱人偷精品视频| 国产男人的电影天堂91| 国产伦在线观看视频一区| 色播亚洲综合网| 国产成人精品久久久久久| 美女 人体艺术 gogo| 美女免费视频网站| 国产精品一区二区免费欧美| 狂野欧美激情性xxxx在线观看| 又黄又爽又免费观看的视频| 久久国内精品自在自线图片| 久久久午夜欧美精品| 亚洲久久久久久中文字幕| 久久久久久伊人网av| 天堂av国产一区二区熟女人妻| 亚洲av一区综合| 波野结衣二区三区在线| 最近2019中文字幕mv第一页| www日本黄色视频网| 又黄又爽又刺激的免费视频.| 欧美一区二区国产精品久久精品| 成人高潮视频无遮挡免费网站| 午夜福利视频1000在线观看| 久久久午夜欧美精品| 亚洲av电影不卡..在线观看| 久久精品影院6| 欧美日韩国产亚洲二区| 一本久久中文字幕| 男人舔奶头视频| 欧洲精品卡2卡3卡4卡5卡区| 99在线人妻在线中文字幕| 国产伦在线观看视频一区| 最后的刺客免费高清国语| 日本免费一区二区三区高清不卡| 欧美极品一区二区三区四区| 国产男人的电影天堂91| 欧美+日韩+精品| 男人的好看免费观看在线视频| 亚洲激情五月婷婷啪啪| 国产久久久一区二区三区| 哪里可以看免费的av片| 伊人久久精品亚洲午夜| 国模一区二区三区四区视频| 久久精品人妻少妇| 亚洲人成网站在线播| 变态另类丝袜制服| 日韩,欧美,国产一区二区三区 | 国产乱人偷精品视频| 亚洲美女搞黄在线观看 | 日韩精品有码人妻一区| 我的女老师完整版在线观看| 人妻夜夜爽99麻豆av| 我要看日韩黄色一级片| 麻豆精品久久久久久蜜桃| 成人永久免费在线观看视频| 欧美不卡视频在线免费观看| 国产白丝娇喘喷水9色精品| 国产亚洲精品久久久久久毛片| av天堂在线播放| 日韩av在线大香蕉| 伊人久久精品亚洲午夜| 男人舔奶头视频| 一夜夜www| 国产欧美日韩一区二区精品| 亚洲无线观看免费| 亚洲在线自拍视频| 你懂的网址亚洲精品在线观看 | 特大巨黑吊av在线直播| 久久久久国产网址| 亚州av有码| 欧美xxxx性猛交bbbb| 身体一侧抽搐| 在线观看av片永久免费下载| 男插女下体视频免费在线播放| 99久国产av精品| 91狼人影院| 插逼视频在线观看| 91久久精品国产一区二区成人| 久久久久国内视频| h日本视频在线播放| av在线天堂中文字幕| 国产精品永久免费网站| 精品人妻视频免费看| 国产爱豆传媒在线观看| 午夜精品在线福利| 亚洲久久久久久中文字幕| av天堂在线播放| 午夜激情欧美在线| 亚洲人成网站高清观看| 久久中文看片网| 少妇熟女欧美另类| 国产精品综合久久久久久久免费| 亚洲国产色片| 天美传媒精品一区二区| 精品欧美国产一区二区三| 男女那种视频在线观看| 看免费成人av毛片| av免费在线看不卡| 精品一区二区三区视频在线| 精品一区二区三区视频在线观看免费| 成人无遮挡网站| 成人特级av手机在线观看| 啦啦啦韩国在线观看视频| 搡老岳熟女国产| 一本一本综合久久| 国内精品美女久久久久久| 成年女人永久免费观看视频| 男插女下体视频免费在线播放| 狂野欧美白嫩少妇大欣赏| 精品福利观看| 久久久久久久午夜电影| 国产伦一二天堂av在线观看| 日本免费a在线| 欧美成人a在线观看| 成人高潮视频无遮挡免费网站| 亚洲电影在线观看av| 99riav亚洲国产免费| 亚洲中文字幕日韩| 久久久久国内视频| 亚洲成av人片在线播放无| av视频在线观看入口| 国产日本99.免费观看| 在线观看av片永久免费下载| 色噜噜av男人的天堂激情| 极品教师在线视频| 12—13女人毛片做爰片一| 成人国产麻豆网| 99热这里只有精品一区| 日韩av不卡免费在线播放| 亚洲国产精品成人综合色| 在线国产一区二区在线| 国产一区二区激情短视频| 免费观看的影片在线观看| 国产亚洲精品久久久久久毛片| 亚洲电影在线观看av|