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

    基于Monte Carlo的風(fēng)力機(jī)葉片強(qiáng)度可靠性分析

    2024-11-05 00:00:00王國(guó)富畢俊喜王妍
    科技創(chuàng)新與應(yīng)用 2024年30期

    摘 要:隨著風(fēng)力機(jī)應(yīng)用數(shù)量增多,以及應(yīng)用規(guī)模擴(kuò)大,風(fēng)力機(jī)葉片可靠度的保證尤為關(guān)鍵。利用Monte Carlo客觀分析葉片強(qiáng)度可靠性,為葉片安全狀態(tài)長(zhǎng)效維持提供依據(jù),實(shí)現(xiàn)風(fēng)力機(jī)的平穩(wěn)運(yùn)轉(zhuǎn)。該文在介紹Monte Carlo的基礎(chǔ)上,分析風(fēng)力機(jī)類型及其葉片狀態(tài)和葉片參數(shù),并以有限元模型構(gòu)建的方式分析風(fēng)力機(jī)葉片載荷及強(qiáng)度,最后探究基于Monte Carlo的風(fēng)力機(jī)葉片可靠性。

    關(guān)鍵詞:Monte Carlo;風(fēng)力機(jī);葉片強(qiáng)度;可靠性分析;有限元模型

    中圖分類號(hào):TK83 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):2095-2945(2024)30-0060-04

    Abstract: As the number of wind turbine applications increases and the scale of applications expands, ensuring the reliability of wind turbine blades is particularly critical. Monte Carlo is used to objectively analyze blade strength reliability to provide a basis for long-term maintenance of blade safety conditions and achieve smooth operation of wind turbines. Based on the introduction by Monte Carlo, this paper analyzes the types of wind turbines, their blade states and blade parameters, analyzes the load and strength of wind turbines blades by means of finite element model construction, and finally explores the reliability of wind turbines blades based on Monte Carlo.

    Keywords: Monte Carlo; wind turbine; blade strength; reliability analysis; finite element model

    風(fēng)力機(jī)葉片物理故障客觀存在,為減少葉片彎矩失效事件,務(wù)必對(duì)葉片強(qiáng)度可靠性分析。當(dāng)Monte Carlo用于分析風(fēng)力機(jī)葉片強(qiáng)度及可靠性,基于有限元模型掌握額定風(fēng)速、切出風(fēng)速2種工況下葉片的失效概率和變異系數(shù)。根據(jù)0.000 28(0.000 35)較低失效概率和14.2%(18.7%)較低變異系數(shù),選定精確度高、適用性強(qiáng)的子集模擬法展開分析。

    1 Monte Carlo概述

    Monte Carlo中文名稱為蒙特卡洛模擬法,是結(jié)構(gòu)可靠度分析法的一種。Monte Carlo程序結(jié)構(gòu)簡(jiǎn)單,受幾何條件影響甚微,易于實(shí)現(xiàn)[1]。其應(yīng)用原理,即圍繞隨機(jī)事物特點(diǎn)進(jìn)行接近真實(shí)性的描述,并逼真模擬物理實(shí)驗(yàn)過程,誤差、精度易確定。這一方法計(jì)算公式如公式(1),以隨機(jī)試驗(yàn)的形式計(jì)算積分,分布密度函數(shù)用f(r)表示,隨機(jī)變量用g(r)表示,視所計(jì)算積分為函數(shù)中自變量的數(shù)學(xué)期望。

    <g>=g(r)f(r)dr 。 (1)

    試驗(yàn)中,若干觀察值r1、r2……rN指函數(shù)中f(r)的N個(gè)子樣;對(duì)應(yīng)觀察值的隨機(jī)變量值,用g(r1)、g(r2)、g(rN)表示,隨機(jī)變量值的算數(shù)平均值(公式(2))又可稱為積分近似值。

    gN=∑g(ri)。 (2)

    Monte Carlo在可靠性分析中能夠提供精確解,但要以模擬次數(shù)為依據(jù),圍繞失效率值分析模擬數(shù)量,當(dāng)失效率小于0.001,說明模擬數(shù)量之多。一般來說, Monte Carlo抽樣次數(shù)在102~104 pf之間,但耗時(shí)較長(zhǎng),不易實(shí)現(xiàn)。因此,要想獲得Monte Carlo高精度解,要借助適當(dāng)方法選擇適合樣本。當(dāng)Monte Carlo用于風(fēng)力機(jī)葉片強(qiáng)度可靠性分析,考慮到風(fēng)力機(jī)葉片彎矩失效事件的失效率偏低,模擬次數(shù)約100萬次,為減少工程量,可在Monte Carlo方法上開發(fā)減方差方法——子集模擬法[2]。

    2 風(fēng)力機(jī)類型及其葉片

    2.1 風(fēng)力機(jī)類型

    以結(jié)構(gòu)為依據(jù)對(duì)風(fēng)力機(jī)分類,一般而言,風(fēng)力機(jī)由葉片、轉(zhuǎn)子、齒輪箱、定子、變槳系統(tǒng)和發(fā)電機(jī)等組成,在組成機(jī)構(gòu)的協(xié)作下實(shí)現(xiàn)空氣動(dòng)能向機(jī)械能的轉(zhuǎn)化,其中,葉片獲取風(fēng)能,并向轉(zhuǎn)子機(jī)構(gòu)傳力[3]。齒輪箱連接低速軸和高速軸,通過轉(zhuǎn)速切換驅(qū)動(dòng)發(fā)電機(jī)工作。

    2.2 葉片狀態(tài)

    風(fēng)力機(jī)葉片是風(fēng)輪的重要部件,葉片工作強(qiáng)度、工作穩(wěn)定性影響風(fēng)力機(jī)性能,當(dāng)葉片可靠性得到保證,則風(fēng)力機(jī)故障率會(huì)大大降低,進(jìn)而實(shí)現(xiàn)風(fēng)力機(jī)常態(tài)化運(yùn)轉(zhuǎn)[4]。由于風(fēng)力機(jī)大多處于野外區(qū)域、沿海區(qū)域,風(fēng)力較大、風(fēng)速較快,極易導(dǎo)致葉片受損、失效,所以要增強(qiáng)葉片抗失效能力,確保葉片強(qiáng)度,使其達(dá)到安全狀態(tài)。葉片強(qiáng)度、葉片應(yīng)力屬于葉片結(jié)構(gòu)分析變量,用Xi表示,則Z是Xi函數(shù),通過特定概率模型進(jìn)行變量表達(dá),葉片結(jié)構(gòu)失效概率用Pf表示,公式如下

    Pf=P(Z<0) 。 (3)

    葉片可靠度指既定時(shí)間內(nèi)和已有條件下,完成預(yù)定功能的概率,用1-Pf表示。針對(duì)風(fēng)力機(jī)葉片狀態(tài)判斷的過程,即葉片強(qiáng)度可靠性分析的過程。

    2.3 葉片參數(shù)

    風(fēng)力機(jī)葉片數(shù)量受風(fēng)機(jī)轉(zhuǎn)速、風(fēng)輪尖速比2個(gè)因素影響,一般情況下,水平軸向風(fēng)力機(jī)處于高速狀態(tài),尖速比在5~15之間,葉片數(shù)目2~4或1~2片。

    葉輪平面直徑(D)用公式(4)計(jì)算

    P=·η1·η2=0.49VDCp η η,(4)

    式中:P為輸出功率,V1為設(shè)計(jì)風(fēng)速,ρ為設(shè)計(jì)風(fēng)速,η、η分別代表發(fā)電機(jī)效率和傳動(dòng)效率,Cp表示風(fēng)能利用系數(shù)。

    葉尖速比(λ0)用公式(5)表示

    λ0= , (5)

    式中:ω、R分別表示風(fēng)力機(jī)葉輪旋轉(zhuǎn)角速度和葉輪圓平面半徑,V1代表葉輪設(shè)計(jì)風(fēng)速。

    3 風(fēng)力機(jī)葉片建模分析

    借助ANSYS有限元分析軟件對(duì)風(fēng)力機(jī)葉片構(gòu)建有限元模型,假設(shè)某風(fēng)力機(jī)額定功率為600 kW,葉片設(shè)計(jì)風(fēng)速13 m/s,風(fēng)能利用系數(shù)0.4,風(fēng)力機(jī)效率0.79,空氣密度1.225 kg/m3。由公式(1)—(4)轉(zhuǎn)化求葉輪直徑(D),結(jié)果為43 m。根據(jù)風(fēng)力機(jī)類型、葉片大小等參數(shù),最終確定尖速比為6。風(fēng)力機(jī)葉片翼型整體式設(shè)計(jì)保證(圖1),保證葉片結(jié)構(gòu)簡(jiǎn)單。

    建模過程中,以上述參數(shù)為依據(jù),并借助SHELL99單元?jiǎng)澐志W(wǎng)格。其中,所選定單元的節(jié)點(diǎn)數(shù)量有8個(gè),單個(gè)節(jié)點(diǎn)自由度有6個(gè)??紤]到葉片結(jié)構(gòu)分層性,通過構(gòu)建有限元幾何模型直觀呈現(xiàn)葉片應(yīng)力信息。以葉片簡(jiǎn)化設(shè)計(jì)為原則,保證風(fēng)力機(jī)減負(fù)運(yùn)行,針對(duì)葉片主梁結(jié)構(gòu)夾芯處理,其中芯料采用輕木和PVC,實(shí)現(xiàn)降本、提質(zhì)、減重、高強(qiáng)。

    4 風(fēng)力機(jī)葉片載荷及強(qiáng)度分析

    4.1 葉片載荷

    風(fēng)力機(jī)葉片承受空氣動(dòng)力、離心力及重力,分別在3種力的作用下,使葉片承受彎曲與扭轉(zhuǎn)、拉伸、拉壓。對(duì)于空氣動(dòng)力載荷,需掌握葉片平面內(nèi)(外)單位長(zhǎng)度受力與葉片半徑的關(guān)系情況:輪轂周圍風(fēng)力機(jī)葉片單位長(zhǎng)度的平面外空氣載荷與葉片半徑大小呈正比例相關(guān),考慮到葉尖損耗的存在,隨著葉片運(yùn)行時(shí)間的延長(zhǎng),則葉尖損耗加劇,一般而言,20 m半徑單位長(zhǎng)度的平面外空氣載荷會(huì)逐漸減小。風(fēng)輪平面內(nèi)空氣載荷在葉片半徑周圍分布較均勻,不會(huì)因半徑增加而加大。對(duì)于離心力載荷,葉片離心力拉力影響因素主要有2個(gè),分別是輪轂附近葉素微元質(zhì)量、微元距離輪轂,確切而言,受微元線速度影響,拉力隨著葉片徑向延長(zhǎng)而增大;受葉素微元質(zhì)量影響下,拉力隨著葉片徑向延長(zhǎng)而減小。對(duì)于重力載荷,來自重力的剪力、壓力受葉素微元質(zhì)量影響,主要呈遞減趨勢(shì)。

    4.2 葉片強(qiáng)度

    借助有限元模型分析風(fēng)力機(jī)葉片強(qiáng)度,對(duì)此明確邊界條件,將葉片簡(jiǎn)化為懸臂梁,葉片根部指圍繞根部截面完全約束。接下來分析有限元數(shù)值,圍繞單個(gè)葉片模擬,針對(duì)微元的空氣動(dòng)力、離心力、重力映射到模型節(jié)點(diǎn),為葉片應(yīng)力分析做準(zhǔn)備,真正掌握葉片強(qiáng)度。隨著葉片運(yùn)行,記錄運(yùn)轉(zhuǎn)方向的變形位移情況,并重點(diǎn)標(biāo)記,相對(duì)而言,葉尖位置的變形、位置幅度較大,葉根處變形、位移幅度相對(duì)較小,運(yùn)轉(zhuǎn)及承載保持同向,變形及位移最大。風(fēng)力機(jī)葉片的材料為玻璃纖維/環(huán)氧,基于各向異性表現(xiàn),所以對(duì)葉片強(qiáng)度校核,全面掌握層合板結(jié)構(gòu)狀態(tài),真實(shí)得知葉片應(yīng)力情況。風(fēng)力機(jī)葉片蒙皮結(jié)構(gòu)為3層合板結(jié)構(gòu),基于Hill-蔡強(qiáng)度理論獲取材料的強(qiáng)度校核數(shù)值,將其與1對(duì)比,如果數(shù)值小于1,說明葉片強(qiáng)度處于安全狀態(tài)。強(qiáng)度校核結(jié)果顯示,第一n8RUFP5Kc5tqP57AIvHZ7bHzdzHk7YNxBbZTFKD+4dA=層合板結(jié)構(gòu)的強(qiáng)度校核數(shù)值為0.368 4,第二層和第三層的強(qiáng)度校核數(shù)值分別是0.425 6和0.826 5,說明三層合板結(jié)構(gòu)安全,風(fēng)機(jī)葉片未破壞。

    從圖1內(nèi)容可知,層板鋪設(shè)纖維材料的顏色發(fā)生變化,顏色漸變的過程,即層板強(qiáng)度增強(qiáng)的過程。

    5 基于Monte Carlo的風(fēng)力機(jī)葉片可靠性分析

    5.1 需求分析

    葉片在風(fēng)力機(jī)建造中的占比達(dá)到18%左右,葉片質(zhì)量關(guān)乎到風(fēng)力機(jī)運(yùn)行效果。由于風(fēng)力機(jī)運(yùn)行環(huán)境惡劣,意味著葉片運(yùn)行期間承受較多載荷,一旦載荷超過葉片承受范圍,極易導(dǎo)致風(fēng)力機(jī)出現(xiàn)故障,因風(fēng)力機(jī)停運(yùn)造成嚴(yán)重的經(jīng)濟(jì)損失[5]?;诖耍獙?duì)葉片可靠性展開分析。在需求導(dǎo)向下,合理選擇葉片可靠性分析方法——Monte Carlo法,這一方法在傳統(tǒng)計(jì)算方法的基礎(chǔ)上升級(jí)而來,滿足高維計(jì)算、便捷計(jì)算的需求。如今Monte Carlo法細(xì)分多種,其中,子集模擬法具有誤差小、效率高、精度高、適用性強(qiáng)的優(yōu)點(diǎn),將其用于風(fēng)力機(jī)葉片可靠性分析具有可行性。

    5.2 翼型參數(shù)

    風(fēng)力機(jī)葉片翼型的參數(shù)見表1,取1—8葉素序號(hào)。

    葉片距根部約18%處達(dá)到最大弦長(zhǎng)4.5 m,轉(zhuǎn)子高度處平均額定風(fēng)速為11.14 m/s,額定轉(zhuǎn)速為12.09 rpm,風(fēng)機(jī)的切入風(fēng)速2.8 m/s,切出風(fēng)速為23.9 m/s。

    風(fēng)力機(jī)翼型葉片在工作中一直處于轉(zhuǎn)動(dòng)狀態(tài),考慮到失速現(xiàn)象的存在,針對(duì)翼型升力系數(shù)、阻力系數(shù)修正,獲得相對(duì)真實(shí)的翼型啟動(dòng)特性,為葉片可靠性分析做準(zhǔn)備。建模分析過程中,針對(duì)葉片根部彎矩重點(diǎn)分析,目的是減少計(jì)算量、縮短計(jì)算時(shí)間,基于梁?jiǎn)卧邢拊P瞳@知葉片葉素的氣動(dòng)力,并顯現(xiàn)于模型節(jié)點(diǎn),具體掌握彎矩情況,針對(duì)葉片根部進(jìn)行加固處理,保證葉片可靠性。

    5.3 環(huán)境模擬

    要想掌握風(fēng)力機(jī)葉片的工作狀態(tài)和強(qiáng)度情況,應(yīng)首先模擬風(fēng)場(chǎng)環(huán)境,根據(jù)風(fēng)的形成原理,以及脈動(dòng)風(fēng)特征、參數(shù),計(jì)算得知風(fēng)載荷。風(fēng)場(chǎng)環(huán)境模擬的過程中,選擇諧波合成法,充分利用計(jì)算機(jī)快速計(jì)算、準(zhǔn)確計(jì)算的優(yōu)勢(shì)完成模擬,順風(fēng)向風(fēng)速與時(shí)間的關(guān)系情況如圖2所示。模擬過程中,每間隔10 min對(duì)不同葉片方向進(jìn)行風(fēng)速記錄,如果實(shí)際運(yùn)行時(shí)間為580 s,前1 min排除瞬態(tài)響應(yīng)的干擾,后520 s為采樣時(shí)間,時(shí)間步長(zhǎng)0.1 s。

    5.4 Monte Carlo——子集模擬

    Monte Carlo用于分析風(fēng)力機(jī)葉片強(qiáng)度可靠性,隨著葉片強(qiáng)度要求、可靠度要求提高,加之,風(fēng)力機(jī)葉片失效概率低,Monte Carlo中子集模擬方法更具適用性。下面圍繞子集模擬的應(yīng)用原理、應(yīng)用過程、計(jì)算分析具體探究。

    5.4.1 應(yīng)用原理

    子集模擬是圍繞每一階段條件失效事件來代表最終小概率失效事件。實(shí)踐過程中,基于馬爾科夫鏈原理形成達(dá)標(biāo)的新樣本,并構(gòu)建真實(shí)樣本空間,面向失效概率空間高效計(jì)算、高效模擬。

    5.4.2 應(yīng)用過程

    第一,確定建議概率密度函數(shù)?;诟咚垢怕史植己瘮?shù),以及0~2π隨機(jī)相位角序列形成新隨機(jī)相位角序列,根據(jù)平穩(wěn)隨機(jī)過程模擬結(jié)果得知高斯分布標(biāo)準(zhǔn)差值。第二,先后進(jìn)行Monte Carlo計(jì)算、生成N個(gè)馬爾科夫鏈樣本、所要求失效概率達(dá)標(biāo)。上述是通過馬爾科夫鏈蒙特卡羅法產(chǎn)生符合條件的樣本的過程。

    5.4.3 計(jì)算分析

    經(jīng)工況模擬求得計(jì)算結(jié)果,工況一:風(fēng)速15 m/s(額定風(fēng)速);工況二:風(fēng)速25 m/s(切出風(fēng)速)。分別進(jìn)行Monte Carlo模擬(數(shù)據(jù)集10萬個(gè))和子集模擬(數(shù)據(jù)集1 000個(gè))。Monte Carlo模擬需CJpPMAE9YPPb86SMslLRuA==要的樣本數(shù)目較多,不同風(fēng)速下的葉片失效概率、變異系數(shù)存在差異,風(fēng)速15 m/s時(shí),失效概率為0.000 30,變異系數(shù)為17.8%;風(fēng)速25 m/s時(shí),失效概率為0.000 41,變異系數(shù)為15.5%。子集模擬所需的樣本數(shù)目較少,風(fēng)速15 m/s時(shí),失效概率為0.000 28,變異系數(shù)為14.2%;風(fēng)速25 m/s時(shí),失效概率為0.000 35,變異系數(shù)為18.7%。見表2。

    比較發(fā)現(xiàn),子集模擬的計(jì)算結(jié)果相對(duì)較小,風(fēng)速15m/s時(shí),失效概率為0.000 28,小于Monte Carlo模擬失效概率的0.000 30;風(fēng)速25 m/s時(shí),失效概率為0.000 35,小于Monte Carlo模擬失效概率的0.000 41。子集模擬在額定風(fēng)速下變異系數(shù)14.2%小于Monte Carlo模擬變異系數(shù)17.8%,說明子集模擬計(jì)算精度較高。綜合來看,子集模擬的樣本量較少,且失效概率偏低,計(jì)算精度偏高,適用于風(fēng)力機(jī)葉片強(qiáng)度可靠性分析,真實(shí)反映葉片工作狀態(tài),為葉片高強(qiáng)度、可靠性工作提供依據(jù),實(shí)現(xiàn)風(fēng)力機(jī)的常態(tài)運(yùn)行。

    6 結(jié)論

    綜上所述,風(fēng)力機(jī)在風(fēng)能向機(jī)械能轉(zhuǎn)化方面發(fā)揮重要作用,風(fēng)力機(jī)葉片是重要部件,其強(qiáng)度及可靠性影響風(fēng)力機(jī)工作效率。通過Monte Carlo中的子集模擬法展開分析,發(fā)現(xiàn)額定風(fēng)速工況下葉片失效概率偏低,且精度較高,能夠以安全狀態(tài)服務(wù)于風(fēng)力機(jī)運(yùn)行。

    參考文獻(xiàn):

    [1] 鞠浩,王旭東,陸佳紅.神經(jīng)網(wǎng)絡(luò)正逆預(yù)測(cè)結(jié)合的風(fēng)力機(jī)葉片強(qiáng)度可靠性研究[J].太陽(yáng)能學(xué)報(bào),2024,45(1):291-298.

    [2] 潘月月,李正農(nóng),張雨坤,等.基于實(shí)測(cè)風(fēng)場(chǎng)的沿海風(fēng)力機(jī)葉片流固耦合特性分析[J].太陽(yáng)能學(xué)報(bào),2023,44(12):330-340.

    [3] 趙清鑫,張?zhí)m挺.基于徑向基神經(jīng)網(wǎng)絡(luò)的風(fēng)力機(jī)葉片鋪層優(yōu)化[J].太陽(yáng)能學(xué)報(bào),2020,41(4):229-234.

    [4] 刁曉航,孫鵬文,馬志坤,等.基于相變量的風(fēng)力機(jī)葉片宏觀拓?fù)鋬?yōu)化設(shè)計(jì)[J].太陽(yáng)能學(xué)報(bào),2023,44(3):198-203.

    [5] 馬航,畢俊喜,葛新宇,等.基于優(yōu)化等壽命疲勞極限模型的風(fēng)力機(jī)葉片壽命預(yù)測(cè)[J].太陽(yáng)能學(xué)報(bào),2023(10):362-369.

    基金項(xiàng)目:自治區(qū)直屬高校基本科研業(yè)務(wù)費(fèi)項(xiàng)目(JY20220305);內(nèi)蒙古自治區(qū)研究生科研創(chuàng)新項(xiàng)目(S20231130Z);內(nèi)蒙古科技計(jì)劃項(xiàng)目(2022YFHH0055)

    第一作者簡(jiǎn)介:王國(guó)富(1996-),男,碩士研究生。研究方向?yàn)閺?fù)雜系統(tǒng)的可靠性。

    *通信作者:畢俊喜(1974-),男,工學(xué)博士,教授,博士研究生導(dǎo)師。研究方向?yàn)閺?fù)雜機(jī)電裝備可靠性設(shè)計(jì),制造、優(yōu)化與控制。

    国产精品爽爽va在线观看网站| 国产精品亚洲美女久久久| 国产熟女午夜一区二区三区| 日本免费a在线| 精品乱码久久久久久99久播| 国内少妇人妻偷人精品xxx网站 | 久久久国产成人免费| 国产精品av久久久久免费| 一区福利在线观看| 国产一级毛片七仙女欲春2| 天堂av国产一区二区熟女人妻 | 看免费av毛片| 精品日产1卡2卡| tocl精华| 欧美日韩乱码在线| 黄片小视频在线播放| 91字幕亚洲| 99国产精品一区二区三区| 在线播放国产精品三级| 午夜激情福利司机影院| 久久精品国产亚洲av高清一级| 国产成年人精品一区二区| 亚洲国产精品sss在线观看| 亚洲熟女毛片儿| 国产熟女xx| 久久午夜综合久久蜜桃| 色av中文字幕| 国产亚洲av嫩草精品影院| 舔av片在线| 国产一区在线观看成人免费| 国产午夜福利久久久久久| 一区二区三区激情视频| 成人永久免费在线观看视频| 精品久久久久久久人妻蜜臀av| 一本综合久久免费| 亚洲中文字幕日韩| 亚洲乱码一区二区免费版| 天堂动漫精品| 亚洲一码二码三码区别大吗| 国产成人aa在线观看| 好男人在线观看高清免费视频| 亚洲第一欧美日韩一区二区三区| 最新美女视频免费是黄的| 看免费av毛片| 激情在线观看视频在线高清| 人妻夜夜爽99麻豆av| 午夜福利18| 亚洲精品在线观看二区| 亚洲av成人不卡在线观看播放网| 国产爱豆传媒在线观看 | 中文字幕熟女人妻在线| 波多野结衣高清作品| 久久精品国产亚洲av香蕉五月| 免费在线观看影片大全网站| 久热爱精品视频在线9| 日本熟妇午夜| 亚洲午夜理论影院| 国产成人av激情在线播放| 欧美黄色淫秽网站| 国产熟女午夜一区二区三区| 51午夜福利影视在线观看| av有码第一页| 五月玫瑰六月丁香| 18禁国产床啪视频网站| 久久亚洲精品不卡| 毛片女人毛片| 欧美乱妇无乱码| 久9热在线精品视频| 国产熟女午夜一区二区三区| 两性夫妻黄色片| 久久久久九九精品影院| www.www免费av| 国产高清有码在线观看视频 | 亚洲精品在线美女| 不卡av一区二区三区| 99国产精品99久久久久| 黄色a级毛片大全视频| 欧洲精品卡2卡3卡4卡5卡区| 99久久精品国产亚洲精品| 欧美色视频一区免费| 中国美女看黄片| 欧美黑人巨大hd| 久久午夜亚洲精品久久| 久久午夜亚洲精品久久| 久久精品成人免费网站| a级毛片在线看网站| 成人永久免费在线观看视频| 国产久久久一区二区三区| 欧美一级毛片孕妇| 很黄的视频免费| 在线观看日韩欧美| 精品熟女少妇八av免费久了| 亚洲人与动物交配视频| 18禁裸乳无遮挡免费网站照片| 亚洲av熟女| 久久人妻福利社区极品人妻图片| 亚洲国产欧美人成| 亚洲国产欧洲综合997久久,| 国产私拍福利视频在线观看| av在线播放免费不卡| 人人妻,人人澡人人爽秒播| 日韩欧美 国产精品| 不卡一级毛片| 变态另类成人亚洲欧美熟女| 亚洲国产精品合色在线| netflix在线观看网站| 精品久久久久久久人妻蜜臀av| 99精品在免费线老司机午夜| 超碰成人久久| 欧美高清成人免费视频www| 欧美乱色亚洲激情| 99国产精品99久久久久| 黄色a级毛片大全视频| 一级毛片精品| 啦啦啦观看免费观看视频高清| 脱女人内裤的视频| 真人做人爱边吃奶动态| av福利片在线观看| netflix在线观看网站| 日日爽夜夜爽网站| 高潮久久久久久久久久久不卡| 午夜久久久久精精品| 成人18禁在线播放| 国产精品永久免费网站| 亚洲欧美日韩高清在线视频| 精品国产乱子伦一区二区三区| 男女那种视频在线观看| 免费观看精品视频网站| 亚洲午夜理论影院| 亚洲午夜理论影院| 免费在线观看视频国产中文字幕亚洲| 成人国产综合亚洲| 手机成人av网站| 亚洲av五月六月丁香网| 精品少妇一区二区三区视频日本电影| 国产精品乱码一区二三区的特点| 欧美日本视频| 亚洲成av人片在线播放无| 中文字幕高清在线视频| 黄色毛片三级朝国网站| 麻豆一二三区av精品| 日日爽夜夜爽网站| 成人欧美大片| 成人一区二区视频在线观看| 韩国av一区二区三区四区| 国产精品亚洲一级av第二区| 少妇粗大呻吟视频| 国产精品免费视频内射| 精品不卡国产一区二区三区| 神马国产精品三级电影在线观看 | 成在线人永久免费视频| 色播亚洲综合网| 99riav亚洲国产免费| 天堂√8在线中文| 美女大奶头视频| 国产免费av片在线观看野外av| 天天添夜夜摸| 琪琪午夜伦伦电影理论片6080| 女生性感内裤真人,穿戴方法视频| 欧美日韩中文字幕国产精品一区二区三区| 欧美性猛交黑人性爽| 欧美性猛交╳xxx乱大交人| 欧美日韩乱码在线| or卡值多少钱| 色播亚洲综合网| 操出白浆在线播放| 日韩有码中文字幕| 亚洲欧美激情综合另类| 久久这里只有精品19| 两人在一起打扑克的视频| 伦理电影免费视频| 99国产精品99久久久久| 欧美日韩瑟瑟在线播放| 久久精品91蜜桃| 又粗又爽又猛毛片免费看| 99久久精品热视频| 九九热线精品视视频播放| 狠狠狠狠99中文字幕| 禁无遮挡网站| 国产日本99.免费观看| www.精华液| 黄色毛片三级朝国网站| 又黄又粗又硬又大视频| 小说图片视频综合网站| ponron亚洲| 一a级毛片在线观看| 1024香蕉在线观看| 男女之事视频高清在线观看| 久久精品国产亚洲av香蕉五月| 男女下面进入的视频免费午夜| 亚洲一区二区三区色噜噜| 国产精品影院久久| 18禁黄网站禁片免费观看直播| 国产不卡一卡二| 国产精品免费一区二区三区在线| 欧美日韩中文字幕国产精品一区二区三区| 性色av乱码一区二区三区2| 久久亚洲真实| 制服诱惑二区| 国内精品久久久久久久电影| 首页视频小说图片口味搜索| 国产精品 欧美亚洲| 脱女人内裤的视频| 中文字幕高清在线视频| 免费在线观看视频国产中文字幕亚洲| 白带黄色成豆腐渣| 12—13女人毛片做爰片一| 最好的美女福利视频网| 午夜福利欧美成人| 国产成年人精品一区二区| 免费观看人在逋| 在线观看免费日韩欧美大片| 99热6这里只有精品| 国内久久婷婷六月综合欲色啪| 观看免费一级毛片| 欧美日韩国产亚洲二区| 欧美+亚洲+日韩+国产| 国产精品久久久人人做人人爽| 中文字幕熟女人妻在线| 久久精品成人免费网站| 国内精品久久久久久久电影| 制服人妻中文乱码| 手机成人av网站| 人人妻人人澡欧美一区二区| 黄片小视频在线播放| 黄色视频,在线免费观看| 国产亚洲av嫩草精品影院| 夜夜夜夜夜久久久久| 蜜桃久久精品国产亚洲av| 国产激情欧美一区二区| 欧美日韩瑟瑟在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕熟女人妻在线| 一个人观看的视频www高清免费观看 | 老鸭窝网址在线观看| 欧美黄色淫秽网站| 99精品在免费线老司机午夜| 欧美 亚洲 国产 日韩一| 国产视频一区二区在线看| 好男人电影高清在线观看| 国产精品香港三级国产av潘金莲| 国产精品一区二区三区四区久久| 淫秽高清视频在线观看| 久久精品成人免费网站| 三级毛片av免费| 黄色丝袜av网址大全| 91在线观看av| 国产亚洲av高清不卡| 久久久久亚洲av毛片大全| 国产免费av片在线观看野外av| 淫秽高清视频在线观看| 中国美女看黄片| 夜夜夜夜夜久久久久| 亚洲人成网站高清观看| 搞女人的毛片| 我的老师免费观看完整版| 国产一区在线观看成人免费| 免费高清视频大片| 国产成人精品无人区| 亚洲人与动物交配视频| 国产久久久一区二区三区| 亚洲av成人一区二区三| 亚洲中文字幕日韩| 久久久国产欧美日韩av| 久久精品综合一区二区三区| 母亲3免费完整高清在线观看| 很黄的视频免费| www.熟女人妻精品国产| 看片在线看免费视频| 757午夜福利合集在线观看| 欧美+亚洲+日韩+国产| 国产蜜桃级精品一区二区三区| 国产精品av久久久久免费| 国产午夜福利久久久久久| 成人三级黄色视频| 人成视频在线观看免费观看| 成人18禁高潮啪啪吃奶动态图| 婷婷精品国产亚洲av| 在线观看一区二区三区| 国产黄a三级三级三级人| 最近最新免费中文字幕在线| 97超级碰碰碰精品色视频在线观看| 久久久久久久久免费视频了| 白带黄色成豆腐渣| 一区福利在线观看| 亚洲九九香蕉| 18禁黄网站禁片午夜丰满| 精品第一国产精品| 亚洲va日本ⅴa欧美va伊人久久| 大型黄色视频在线免费观看| 亚洲精品美女久久久久99蜜臀| tocl精华| 精品久久久久久久久久久久久| 国产伦一二天堂av在线观看| 12—13女人毛片做爰片一| 久久香蕉国产精品| 高清毛片免费观看视频网站| 国产精品久久电影中文字幕| 国产免费男女视频| 村上凉子中文字幕在线| 国产精品av久久久久免费| 精品熟女少妇八av免费久了| 精品一区二区三区av网在线观看| 久热爱精品视频在线9| 国产真人三级小视频在线观看| 亚洲狠狠婷婷综合久久图片| 久久久国产成人免费| 国产精品 欧美亚洲| 国产精品乱码一区二三区的特点| 丁香欧美五月| 久久中文字幕一级| x7x7x7水蜜桃| 久久久久免费精品人妻一区二区| www.自偷自拍.com| 男女床上黄色一级片免费看| 久久久久九九精品影院| 成人亚洲精品av一区二区| 天天一区二区日本电影三级| 2021天堂中文幕一二区在线观| 嫩草影院精品99| 黄片大片在线免费观看| 在线播放国产精品三级| 超碰成人久久| 免费看a级黄色片| 国产一区二区在线av高清观看| 99riav亚洲国产免费| 99热这里只有精品一区 | 亚洲精品美女久久久久99蜜臀| 国产成人影院久久av| www.精华液| 在线a可以看的网站| 一a级毛片在线观看| 日本a在线网址| 欧美丝袜亚洲另类 | 1024香蕉在线观看| 男女视频在线观看网站免费 | 婷婷精品国产亚洲av在线| 九九热线精品视视频播放| 免费在线观看成人毛片| 国产亚洲精品一区二区www| 18禁黄网站禁片免费观看直播| 性色av乱码一区二区三区2| 老司机午夜十八禁免费视频| 欧美高清成人免费视频www| 不卡一级毛片| 欧美黑人欧美精品刺激| 最近在线观看免费完整版| 亚洲第一电影网av| 亚洲电影在线观看av| 亚洲片人在线观看| 日韩大码丰满熟妇| 黑人欧美特级aaaaaa片| 成人特级黄色片久久久久久久| 一进一出好大好爽视频| 十八禁人妻一区二区| 99热这里只有是精品50| 亚洲精品一区av在线观看| 亚洲人成77777在线视频| 无人区码免费观看不卡| 精品第一国产精品| 日本黄大片高清| 波多野结衣高清作品| 99在线视频只有这里精品首页| 精品国产乱子伦一区二区三区| 欧美成人午夜精品| 国产欧美日韩精品亚洲av| 国产熟女午夜一区二区三区| 亚洲 欧美一区二区三区| 成人特级黄色片久久久久久久| 丰满的人妻完整版| 1024视频免费在线观看| 午夜久久久久精精品| 亚洲国产高清在线一区二区三| 精品乱码久久久久久99久播| 1024手机看黄色片| 中文字幕久久专区| 伊人久久大香线蕉亚洲五| 91成年电影在线观看| av片东京热男人的天堂| 香蕉久久夜色| www日本黄色视频网| 免费搜索国产男女视频| 日日夜夜操网爽| 国产三级在线视频| 精品久久久久久久末码| 日本一本二区三区精品| 免费在线观看影片大全网站| 一个人免费在线观看电影 | 色av中文字幕| 中文亚洲av片在线观看爽| 国产午夜精品久久久久久| 首页视频小说图片口味搜索| 一级黄色大片毛片| 琪琪午夜伦伦电影理论片6080| 成年免费大片在线观看| av免费在线观看网站| 宅男免费午夜| 女警被强在线播放| 麻豆一二三区av精品| 色综合亚洲欧美另类图片| 狂野欧美白嫩少妇大欣赏| 天天一区二区日本电影三级| 香蕉av资源在线| 操出白浆在线播放| 激情在线观看视频在线高清| 国产精品,欧美在线| 国产成人aa在线观看| 无限看片的www在线观看| 一a级毛片在线观看| 国产欧美日韩一区二区精品| 亚洲一区二区三区不卡视频| 国产精品av视频在线免费观看| 最近最新免费中文字幕在线| 脱女人内裤的视频| 国产野战对白在线观看| 啪啪无遮挡十八禁网站| 两个人看的免费小视频| 日韩免费av在线播放| 日本a在线网址| 欧美久久黑人一区二区| ponron亚洲| а√天堂www在线а√下载| 亚洲第一电影网av| 99久久久亚洲精品蜜臀av| 欧美性猛交黑人性爽| 色av中文字幕| 精品国产超薄肉色丝袜足j| 欧美一区二区国产精品久久精品 | 日本一本二区三区精品| 嫩草影视91久久| 性色av乱码一区二区三区2| 精品久久久久久久久久久久久| 12—13女人毛片做爰片一| 久久人妻福利社区极品人妻图片| 亚洲国产精品久久男人天堂| 久久久久久免费高清国产稀缺| 嫩草影院精品99| 国产片内射在线| 两个人的视频大全免费| 一夜夜www| 99热这里只有是精品50| 亚洲精品一区av在线观看| 色精品久久人妻99蜜桃| 国产激情欧美一区二区| 啦啦啦免费观看视频1| 一本久久中文字幕| 久久久久久久久中文| 欧美在线黄色| 毛片女人毛片| 国产免费男女视频| 宅男免费午夜| 国产97色在线日韩免费| 日本黄色视频三级网站网址| 午夜福利免费观看在线| 国产久久久一区二区三区| 精品电影一区二区在线| 午夜福利在线观看吧| 欧美日韩瑟瑟在线播放| 黑人巨大精品欧美一区二区mp4| 久久久久久大精品| xxx96com| 天天躁狠狠躁夜夜躁狠狠躁| 午夜免费激情av| 久久精品国产综合久久久| 两个人免费观看高清视频| e午夜精品久久久久久久| 一二三四在线观看免费中文在| 露出奶头的视频| 亚洲第一电影网av| 免费看a级黄色片| 亚洲色图 男人天堂 中文字幕| 国产精品亚洲美女久久久| 久久久水蜜桃国产精品网| 亚洲精品国产一区二区精华液| 男人舔奶头视频| 桃红色精品国产亚洲av| 国产伦一二天堂av在线观看| 长腿黑丝高跟| 男女那种视频在线观看| 色噜噜av男人的天堂激情| 午夜老司机福利片| 国产精品免费视频内射| 国产黄色小视频在线观看| 91老司机精品| 一区福利在线观看| 国产成人欧美在线观看| 国产av在哪里看| 亚洲欧美日韩高清专用| 免费无遮挡裸体视频| av在线天堂中文字幕| 高清在线国产一区| 精品久久久久久,| 久99久视频精品免费| 亚洲欧美精品综合一区二区三区| 九九热线精品视视频播放| 伦理电影免费视频| 午夜影院日韩av| 一级作爱视频免费观看| 三级男女做爰猛烈吃奶摸视频| 极品教师在线免费播放| 日日爽夜夜爽网站| 午夜免费观看网址| www国产在线视频色| 日韩高清综合在线| 精品国内亚洲2022精品成人| 亚洲成人国产一区在线观看| 99久久综合精品五月天人人| 成人国产一区最新在线观看| 99精品久久久久人妻精品| 麻豆成人av在线观看| 色精品久久人妻99蜜桃| www日本在线高清视频| 欧美一区二区国产精品久久精品 | 黑人操中国人逼视频| 女生性感内裤真人,穿戴方法视频| 精品国内亚洲2022精品成人| 免费在线观看完整版高清| 色综合欧美亚洲国产小说| av中文乱码字幕在线| 好男人在线观看高清免费视频| 国产精品免费视频内射| 日韩欧美免费精品| 精品久久久久久久人妻蜜臀av| 91麻豆精品激情在线观看国产| 色av中文字幕| 久久国产乱子伦精品免费另类| 操出白浆在线播放| 久久精品91蜜桃| 亚洲一码二码三码区别大吗| 99久久99久久久精品蜜桃| 国产99白浆流出| 久久九九热精品免费| 欧美三级亚洲精品| 禁无遮挡网站| 两个人的视频大全免费| 久久香蕉激情| 又紧又爽又黄一区二区| 一夜夜www| 9191精品国产免费久久| 久久久久久久久久黄片| 成人18禁在线播放| 啦啦啦免费观看视频1| 欧美色欧美亚洲另类二区| 国产爱豆传媒在线观看 | 欧美日韩国产亚洲二区| 欧美性猛交黑人性爽| 欧美最黄视频在线播放免费| 久久中文看片网| 国产片内射在线| 欧美黑人精品巨大| 色综合欧美亚洲国产小说| 国产野战对白在线观看| 国产成人一区二区三区免费视频网站| 久久久久久免费高清国产稀缺| 免费看十八禁软件| 老司机在亚洲福利影院| а√天堂www在线а√下载| 亚洲一区二区三区不卡视频| 精品久久久久久久毛片微露脸| 成人欧美大片| 变态另类丝袜制服| 两个人看的免费小视频| 香蕉av资源在线| 国产精品永久免费网站| 欧美黑人精品巨大| 国产激情久久老熟女| 精华霜和精华液先用哪个| 少妇的丰满在线观看| 欧美日本亚洲视频在线播放| 精品无人区乱码1区二区| 亚洲专区中文字幕在线| 亚洲一卡2卡3卡4卡5卡精品中文| 动漫黄色视频在线观看| 男人的好看免费观看在线视频 | 一夜夜www| 欧美 亚洲 国产 日韩一| 久久人妻福利社区极品人妻图片| 韩国av一区二区三区四区| 最近视频中文字幕2019在线8| 亚洲av中文字字幕乱码综合| 日韩三级视频一区二区三区| 看免费av毛片| 哪里可以看免费的av片| 99精品在免费线老司机午夜| 国产av不卡久久| 国产探花在线观看一区二区| 一二三四社区在线视频社区8| 亚洲欧美日韩高清在线视频| 五月玫瑰六月丁香| 50天的宝宝边吃奶边哭怎么回事| 亚洲av美国av| 啦啦啦观看免费观看视频高清| 五月玫瑰六月丁香| 精品国内亚洲2022精品成人| 变态另类丝袜制服| 免费一级毛片在线播放高清视频| 19禁男女啪啪无遮挡网站| 波多野结衣巨乳人妻| 亚洲av日韩精品久久久久久密| 在线国产一区二区在线| 视频区欧美日本亚洲| 久久精品国产综合久久久| 亚洲九九香蕉| 亚洲va日本ⅴa欧美va伊人久久| 亚洲七黄色美女视频| 国产伦在线观看视频一区| 精品免费久久久久久久清纯| 亚洲国产精品合色在线| 99久久无色码亚洲精品果冻| 婷婷精品国产亚洲av| 这个男人来自地球电影免费观看| 全区人妻精品视频| 两性夫妻黄色片| ponron亚洲| 国产区一区二久久|