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

    基于多核CPU的復(fù)雜液壓產(chǎn)品快速并行優(yōu)化方法

    2022-05-12 09:32:12寧志強(qiáng)衛(wèi)立新化建輝高有山
    關(guān)鍵詞:配流柱塞泵對稱軸

    寧志強(qiáng) 衛(wèi)立新 張 瑞 權(quán) 龍 化建輝 高有山

    (1.太原科技大學(xué)機(jī)械工程學(xué)院,太原 030024;2.山西工程技術(shù)學(xué)院機(jī)械電子工程系,陽泉 045000;3.太原理工大學(xué)機(jī)械與運(yùn)載工程學(xué)院,太原 030002)

    0 引言

    復(fù)雜機(jī)械產(chǎn)品的參數(shù)如尺寸、零部件選型等將對產(chǎn)品性能產(chǎn)生影響,同時產(chǎn)品性能也會有各種限制條件,復(fù)雜機(jī)械產(chǎn)品優(yōu)化設(shè)計往往具有多維非線性特征,因此獲取全局最優(yōu)解需要較大計算量。復(fù)雜產(chǎn)品關(guān)鍵參數(shù)的優(yōu)化問題是機(jī)械優(yōu)化設(shè)計的一個重要研究方向[1-2]。液壓技術(shù)在工業(yè)設(shè)備領(lǐng)域應(yīng)用廣泛,對液壓產(chǎn)品性能要求逐漸向高壓、高速、大功率的方向發(fā)展。隨著計算機(jī)硬件性能的極大提升,液壓仿真技術(shù)也得到飛速發(fā)展,逐漸成為液壓產(chǎn)品設(shè)計人員必不可少的開發(fā)工具[3]。

    智能優(yōu)化算法與液壓仿真軟件結(jié)合是一種發(fā)展趨勢,可極大減輕設(shè)計人員重復(fù)設(shè)計過程[4-5]。德國亞深工業(yè)大學(xué)早期開發(fā)的DSH是一款液壓系統(tǒng)仿真軟件,具有面向原理圖建模的功能,但不具備優(yōu)化模塊。浙江大學(xué)對DSH軟件進(jìn)行了二次開發(fā),擴(kuò)充了一些液壓仿真模塊和優(yōu)化模塊[6]。國內(nèi)學(xué)者運(yùn)用優(yōu)化工具對液壓產(chǎn)品進(jìn)行了參數(shù)設(shè)計[7-10]。

    當(dāng)前液壓仿真軟件進(jìn)行優(yōu)化設(shè)計和人機(jī)交互方面存在一些問題。液壓產(chǎn)品進(jìn)行參數(shù)優(yōu)化時,需提前安裝專業(yè)仿真軟件,無法脫離仿真平臺獨立運(yùn)行在Windows操作系統(tǒng)。在仿真軟件平臺進(jìn)行參數(shù)優(yōu)化時,主要是以串行迭代進(jìn)行優(yōu)化,專業(yè)仿真平臺對內(nèi)存和CPU需求很高,一臺計算機(jī)同時運(yùn)行兩個仿真平臺非常緩慢,因此基于仿真平臺軟件的并行仿真不易實現(xiàn),且高性能工作站攜帶不便。此外,仿真軟件的人機(jī)交互界面不友好,更改參數(shù)初始值需手動輸入,模型參數(shù)較多時,手動修改參數(shù)操作煩瑣,重復(fù)工作量大。

    為解決以上問題,多核 CPU加速技術(shù)是一個可行的途徑。多核CPU技術(shù)的出現(xiàn)使得并行計算的廣泛應(yīng)用成為可能,計算機(jī)上多核CPU技術(shù)使并行計算技術(shù)的應(yīng)用推向微型計算機(jī)[11]。多核 CPU 相比早期單核處理器有了較大的計算性能提升,且處理器的核心數(shù)有增多的趨勢,主要代表性企業(yè)有Intel和 AMD 公司[12]。國內(nèi)外文獻(xiàn)中關(guān)于液壓產(chǎn)品多核并行優(yōu)化方法的相關(guān)報道較少。本文提出一種基于多核CPU的復(fù)雜液壓產(chǎn)品快速多進(jìn)程并行優(yōu)化方法。以三角槽參數(shù)優(yōu)化為研究對象,利用泵出口壓力試驗驗證COVDE求解器仿真模型的準(zhǔn)確性,在此基礎(chǔ)上通過粒子群算法尋優(yōu)來降低輸出流量脈動。

    1 多核并行優(yōu)化方法框架

    多核CPU的復(fù)雜液壓產(chǎn)品快速并行優(yōu)化方法的主要原理為基于群體智能優(yōu)化算法,對CVODE仿真程序動態(tài)調(diào)整仿真參數(shù)。CVODE求解器能提高仿真速度,是本文提出的第一級加速策略。COVDE可執(zhí)行仿真程序看作群體智能算法的群體粒子,每個粒子(COVDE可執(zhí)行程序)被分配給不同的設(shè)計參數(shù),參數(shù)化過程通過對參數(shù)文件的寫入來實現(xiàn)。從COVDE可執(zhí)行程序生成的輸出文件中可提取產(chǎn)品性能參數(shù),并計算出群體粒子的適應(yīng)值,同時判斷是否滿足產(chǎn)品性能約束條件。區(qū)別于傳統(tǒng)的串行優(yōu)化仿真方法,每個COVDE可執(zhí)行程序是一個dos程序,通過ShellExecute函數(shù)即可啟動運(yùn)行,自動生成一個進(jìn)程,同時啟動多個dos程序則進(jìn)程分配權(quán)歸操作系統(tǒng)所有,不需要人為干預(yù)。以多個dos進(jìn)程同時進(jìn)行仿真可充分利用多核處理器,與群體優(yōu)化算法模塊進(jìn)行數(shù)據(jù)交互,直到滿足迭代終止條件?;谌后w智能優(yōu)化算法的參數(shù)分配模塊通過群體進(jìn)化策略來動態(tài)分配設(shè)計參數(shù)。鑒于專業(yè)液壓仿真軟件人機(jī)交互性差,可視化交互模塊也是必要的組成部分。多核CPU并行優(yōu)化方法主要框架結(jié)構(gòu)如圖1所示。

    圖1 多核并行優(yōu)化方法框圖Fig.1 Frame of multi-core parallelism optimization

    2 非對稱軸向柱塞泵配流盤設(shè)計

    變排量非對稱軸向柱塞泵(VAPP)可用于挖掘機(jī)動臂能量回收[13-16]。VAPP與普通泵的最大區(qū)別在于把普通配流盤的兩個配流窗口設(shè)計成為3個配流窗口,配流窗口A和B分別接差動缸的無桿腔和有桿腔。通過改變斜盤傾角使配流窗口A和B獲得不同的流量達(dá)到驅(qū)動液壓缸運(yùn)動,而不需要其他輔助元件。通過改變斜盤角度可使VAPP處于馬達(dá)工況或泵工況,處于馬達(dá)狀態(tài)時差動缸無桿腔的部分液壓油通過配流窗口T存在蓄能器中。當(dāng)VAPP處于泵工況時,蓄能器中的液壓能被釋放出來,減少電機(jī)做功,從而達(dá)到勢能回收和利用的目的。非對稱軸向柱塞泵結(jié)構(gòu)如圖2所示。

    圖2 非對稱軸向柱塞泵結(jié)構(gòu)圖Fig.2 VAPP structure diagram

    景健[17]對非對稱軸向柱塞泵配流盤建模時,將三角槽結(jié)構(gòu)配流面積的變化過程劃分成腰型孔變化區(qū)域和三角槽變化區(qū)域,并對非對稱泵物理樣機(jī)進(jìn)行了試驗。試驗結(jié)果表明該方法能夠合理模擬非對稱軸向柱塞泵配流面積。因此,本文采用該方法對非對稱軸向柱塞泵建模。

    三角槽不同階段的劃分如圖3所示,具體配流面積計算方法如下:

    圖3 三角槽配流窗口過程劃分Fig.3 Dividing of flow procedure of triangle groove

    (1)加速增大階段(0<φ≤φ1max)

    在加速增大階段,柱塞腔底孔僅與三角槽部分連通。配流面積定義為位于柱塞腔底孔與三角槽相交處的三角槽頂面所包圍的面積。圖4為配流盤三角槽的結(jié)構(gòu),柱塞腔底孔與三角槽之間的配流面積為

    圖4 配流盤三角槽結(jié)構(gòu)圖Fig.4 Triangle groove structure of valve plate

    (1)

    式中l(wèi)——三角槽寬度(底邊長)

    L——三角槽長度φ——轉(zhuǎn)角

    R——柱塞腔底孔分布圓半徑

    (2)穩(wěn)定階段(φ1max<φ≤φ2max)

    在穩(wěn)定階段柱塞腔底孔與配流盤三角槽完全接通,此時配流面積為

    (2)

    (3)加速減小階段(φ2max<φ≤φ3max)

    在加速減小階段,柱塞腔底孔與三角槽部分接通,其配流面積為

    (3)

    腰型孔變化區(qū)域內(nèi)的配流面積分為弓形增大、線性增大、穩(wěn)定、線性減小和弓形減小,如圖5所示配流面積計算方法為:

    圖5 腰型孔配流窗口過程劃分Fig.5 Dividing of flow procedure of waist window

    (1)弓形增大階段(φ1max<φ≤φ4max)

    柱塞腔底孔與腰型孔接通,其配流面積形狀為弓形,如圖5所示。配流面積為

    (4)

    其中

    φ′=φ4max-φ

    (5)

    式中Ssec——柱塞腔底孔與腰型孔形成的扇形面積

    SΔ——扇形面積中包含的三角形面積

    r——腰形槽半徑

    (2)線性增大階段(φ4max<φ≤φ5max)

    此時柱塞腔底孔與配流盤腰型孔重疊面積輪廓為腰型輪廓。其配流面積為

    S3=πR2+2Rr(φ-φ4max)

    (6)

    (3)穩(wěn)定階段(φ5max<φ≤φ6max)

    柱塞腔底孔與配流盤腰型孔完全連通,如圖6所示,其配流面積為

    圖6 腰型孔配流窗口過程劃分Fig.6 Dividing of flow procedure of waist window

    S3=πr2+2Rr(φ5max-φ4max)

    (7)

    (4)線性減小階段(φ6max<φ≤φ7max)

    配流面積形狀仍為腰型孔輪廓,但其變化趨勢為線性減小,如圖6所示,配流面積為

    S3=πr2+2Rr(φ7max-φ)

    (8)

    (5)弓形減小階段(φ7max<φ≤φ8max)

    柱塞腔底孔與配流盤重疊面積輪廓變?yōu)楣危鐖D6所示,配流面積為

    (9)

    其中

    φ′=φ-φ7max

    (10)

    3 CVODE仿真模型試驗驗證

    3.1 非對稱泵CVODE求解器可執(zhí)行程序

    在SimulationX軟件中建立非對稱泵變量機(jī)構(gòu)仿真模型,如圖7所示。SimulationX 軟件采用物理對象建模方法,且能提供模型導(dǎo)出代碼功能,SimulationX 模型能以 C 代碼形式導(dǎo)出,導(dǎo)出時可選擇定步長求解器和COVDE求解器。導(dǎo)出的C代碼可生成獨立的后綴名為exe的CVODE可執(zhí)行程序、同時生成具有固定文件名及文件類型的模型參數(shù)文件(parameters.txt)、仿真控制參數(shù)文件(solversettings.txt)和仿真結(jié)果文件(Outputs1.txt)。所生成的dos可執(zhí)行文件能夠在Windows操作系統(tǒng)下運(yùn)行,不需要安裝SimulationX 軟件平臺就可以運(yùn)行。通過VC++自動執(zhí)行CVODE仿真程序的調(diào)用函數(shù)為ShellExecute(NULL,"open","CVODE.exe",NULL,NULL,SW_SHOWNORMAL)。

    圖7 非對稱軸向柱塞泵SimulationX仿真模型Fig.7 VAPP model of SimulationX

    根據(jù)SimulationX 軟件中非對稱軸向柱塞泵模型,導(dǎo)出C代碼。VC++6.0通過Open workspace打開C代碼中后綴名為dsw的文件,編譯生成CVODE求解器的dos可執(zhí)行文件,其運(yùn)行界面如圖8所示。

    圖8 CVODE仿真程序Fig.8 CVODE executable program

    3.2 試驗驗證CVODE求解器的準(zhǔn)確性

    通過樣機(jī)試驗驗證CVODE求解器的dos可執(zhí)行程序仿真的準(zhǔn)確性。試驗過程所使用的設(shè)備參數(shù)見表1,圖9為樣機(jī)實物圖和試驗方案[18-20]。試驗和仿真條件設(shè)為:非對稱軸向柱塞泵轉(zhuǎn)速為600 r/min,斜盤傾角設(shè)定為18°,B口和T口加載壓力分別取5、10、15、20 MPa。當(dāng)對B口進(jìn)行加載,且T口不加載時,B口在4個不同壓力等級下,樣機(jī)試驗和仿真結(jié)果如圖10a所示。從圖10a中可知,仿真與試驗結(jié)果基本吻合,B口的壓力脈動隨著加載壓力的提高而增大,總體脈動幅值仍然較小。當(dāng)對B口和T口同時加載壓力時,B口和T口分別在4種壓力等級下的壓力情況如圖10b、10c所示。仿真結(jié)果和試驗結(jié)果基本一致,B口和T口的工作壓力呈現(xiàn)一定周期性波動,且T口壓力脈動幅值略大于窗口B。通過試驗驗證,結(jié)果表明壓力實測和CVODE可執(zhí)行程序得到的仿真曲線吻合度較高。

    表1 試驗設(shè)備參數(shù)Tab.1 Proposed system parameters

    圖9 變排量非對稱泵樣機(jī)和試驗方案Fig.9 VAPP prototype and experiment method1.電動機(jī) 2.非對稱泵 3.比例伺服閥 4.液壓泵 5.變量機(jī)構(gòu)電動機(jī) 6、7.溢流閥 8.DSpace控制器

    4 配流盤并行優(yōu)化的實現(xiàn)過程

    針對液壓產(chǎn)品開發(fā)特點,設(shè)計了一種多核并行優(yōu)化機(jī)制。非對稱軸向柱塞泵的模型復(fù)雜度相對較高,因此可以在仿真方案中設(shè)計多個可獨立運(yùn)行的仿真程序,分別賦不同的參數(shù)值,在經(jīng)過一次迭代后,微粒群算法優(yōu)化器能從輸出文件output.txt中獲取計算結(jié)果,得到粒子群算法所需的的目標(biāo)函數(shù)適應(yīng)值和約束條件判斷結(jié)果,粒子群算法更新每個粒子的參數(shù)。粒子群算法具有并行性,每個粒子的適應(yīng)值計算相互獨立,因此能夠同時多進(jìn)程并行運(yùn)行多個可執(zhí)行程序,以實現(xiàn)提高計算效率的目的。多進(jìn)程運(yùn)行的啟動方式借助ShellExecute函數(shù)實現(xiàn)。

    4.1 三角槽優(yōu)化參數(shù)選取

    三角槽長度和寬度對流量脈動和柱塞腔壓力具有較大的影響,鑒于三配流口的配流槽尺寸限定,三角槽長度優(yōu)化的空間較小,因此選取三角槽寬度作為設(shè)計變量。為研究三角槽寬度對B口和T口流量脈動影響,選取不同的三角槽寬度進(jìn)行仿真,分別取三角槽寬度為1、3 mm進(jìn)行仿真。由圖11a可知,隨著B口三角槽寬度變小,B口流量脈動變小;由圖11b可知,隨著T口三角槽寬度變小,T口流量脈動變??;由圖11c可知,隨著B口和T口三角槽寬度變小,柱塞腔壓力變大。在非對稱泵設(shè)計中,要求泵出口流量脈動較小,同時也要求柱塞腔壓力較小,但流量脈動較小和柱塞腔壓力較小的要求存在沖突,因此,本文提出通過改變B口和T口三角槽寬度,在限制最大柱塞腔壓力的同時,達(dá)到B口和T口的流量脈動值最優(yōu)化。

    圖11 不同寬度的三角槽對流量脈動和柱塞腔壓力的影響Fig.11 Effect of different triangle groove on flow pulsation and chamber pressure

    4.2 目標(biāo)函數(shù)和約束條件

    粒子群算法的目標(biāo)函數(shù)設(shè)置為B口和T口的流量脈動率之和,即目標(biāo)函數(shù)

    σall=σB+σT

    (11)

    其中

    (12)

    式中qBmax、qTmax——B口、T口最大流量

    qBmin、qTmin——B口、T口最小流量

    優(yōu)化設(shè)計變量分別為三角槽B口和T口的寬度lB和lT,如圖4b所示。約束條件為柱塞腔最大壓力不得大于限制值。約束條件的施加方法,最常見的有懲罰函數(shù)法等方法,但無法將約束條件公式化,不適合采用懲罰函數(shù)法,因此提出一種約束條件處理策略,通過3個規(guī)則處理約束條件。基于粒子群算法的非對稱泵多核并行優(yōu)化流程圖如圖12所示。

    圖12 非對稱泵多核并行優(yōu)化流程圖Fig.12 Flowchart of VAPP parallelism optimization

    4.3 三角槽多核并行優(yōu)化對比

    利用VC++軟件進(jìn)行SimulationX軟件的二次開發(fā),編制出友好的人機(jī)交互界面,如圖13a所示,多進(jìn)程并行優(yōu)化界面如圖13b所示。

    圖13 二次開發(fā)界面Fig.13 Secondary development interface

    優(yōu)化前配流盤B口和T口三角槽寬度均為2 mm,流量脈動率之和為0.763 268,柱塞腔壓力為2.3 MPa。在20輪迭代過程中,約束條件設(shè)為柱塞腔壓力為2.3 MPa,優(yōu)化后B口和T口三角槽寬度分別為3 mm和1.46 mm,B口和T口脈動率之和下降為0.491 952 78,即柱塞腔壓力不變的情況下,通過優(yōu)化三角槽寬度使出口脈動率降低了36%。排量為45 mL/r時20次迭代過程脈動率變化趨勢如圖14所示,優(yōu)化后流量和壓力如圖15所示。

    圖14 20次迭代過程出口脈動率變化曲線Fig.14 Variation trend of pulse rate in 20 iterative process

    圖15 三角槽參數(shù)優(yōu)化后的輸出流量和柱塞腔壓力Fig.15 Output flow and plunger chamber pressure of optimized triangle groove

    對于非對稱軸向柱塞泵的設(shè)計,輸入設(shè)定排量為55 mL/r,約束條件設(shè)定柱塞腔壓力不超過2.3 MPa,優(yōu)化后B口和T口三角槽寬度分別為1.24、2.62 mm,無法到達(dá)搜索范圍下限,因為三角槽寬度下限值將導(dǎo)致柱塞腔壓力大于許可值2.3 MPa。排量55 mL/r時,在5 MPa柱塞腔壓力限制條件下,三角槽寬度達(dá)到搜索范圍下限值,即0.5 mm,此時柱塞腔壓力約束條件不起作用。排量55 mL/r相比于排量45 mL/r,脈動率下降,即隨著排量的增加,泵出口脈動率變小。不同排量和壓力限制條件下三角槽尺寸尋優(yōu)結(jié)果見表2。圖16為不同排量和壓力限制條件下三角槽優(yōu)化后輸出流量脈動率。

    表2 三角槽寬度優(yōu)化結(jié)果Tab.2 Optimized result of triangle groove

    圖16 不同排量和壓力限制條件下三角槽優(yōu)化后輸出流量脈動率Fig.16 Output flow and plunger chamber pressure of optimized triangle groove with different displacements and pressure limits

    5 不同核心數(shù)目處理器對運(yùn)行速度的影響

    設(shè)置并行優(yōu)化進(jìn)行20次迭代,每輪迭代同時并行執(zhí)行10個仿真程序,共計執(zhí)行200次dos仿真程序。對比試驗的仿真處理器配置見表3,仿真耗時對比見表4。對于配置A,Intel(R)Core(TM)i7-10700F是具有8核16線程的處理器,理論上該處理器相比單個exe至少可提高4倍效率,實際提高效率為7倍左右,基本符合理論上提升效率;對于配置B,具有2核處理器,理論上該處理器相比單個exe至少可提高2倍效率,實際提高效率約為2.5倍,基本符合理論提升效率。8核處理器i7-10700F對比雙核處理器i5-7200U提升約5倍效率。并行優(yōu)化性能與核心數(shù)基本成正比,但也與仿真時所用計算機(jī)的實際性能相關(guān)。

    表3 工作站主要性能參數(shù)Tab.3 Parameters of workstation

    表4 仿真耗時對比Tab.4 Comparison of simulation time s

    6 結(jié)論

    (1)提出了一種基于多核CPU的復(fù)雜液壓產(chǎn)品快速多進(jìn)程并行優(yōu)化方法,該方法可解決液壓產(chǎn)品仿真對專業(yè)軟件依賴的問題。CVODE求解器能加快仿真速度,物理樣機(jī)試驗結(jié)果與CVODE求解器仿真程序結(jié)果吻合度較高,證明CVODE仿真程序的準(zhǔn)確性。在基于粒子群算法的并行優(yōu)化過程的運(yùn)行速率與多核CPU的核心數(shù)目呈正比關(guān)系。

    (2)利用粒子群算法對非對稱軸向柱塞泵的三角槽寬度參數(shù)進(jìn)行優(yōu)化,以降低泵輸出流量脈動。研究表明減小三角槽寬度能夠降低流量脈動,但會引起柱塞腔壓力上升。在不增加柱塞腔壓力的條件下,運(yùn)用多核并行優(yōu)化方法對非對稱軸向柱塞泵優(yōu)化后,泵出口流量脈動相比優(yōu)化前降低了36%。

    猜你喜歡
    配流柱塞泵對稱軸
    重載柱塞泵球面配流副承載特性研究*
    基于HP濾波與ARIMA-GARCH模型的柱塞泵泄漏量預(yù)測
    某型飛機(jī)液壓柱塞泵彈簧斷裂分析
    微觀織構(gòu)配流副熱-流-固耦合潤滑特性
    PSO-BP網(wǎng)絡(luò)模型在柱塞泵故障診斷中的應(yīng)用
    電子測試(2018年11期)2018-06-26 05:56:02
    先找對稱軸!
    抓牢對稱軸突破二次函數(shù)
    有幾條對稱軸
    基于改進(jìn)最小二乘支持向量機(jī)的柱塞泵故障診斷模型
    鐵路編組站動態(tài)配流分層模型
    亚洲中文字幕日韩| 久久久国产成人免费| 伊人亚洲综合成人网| 伊人亚洲综合成人网| 亚洲天堂av无毛| 国产精品秋霞免费鲁丝片| 欧美亚洲 丝袜 人妻 在线| 汤姆久久久久久久影院中文字幕| 麻豆av在线久日| 99香蕉大伊视频| 老司机午夜十八禁免费视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩亚洲国产一区二区在线观看 | 久久国产精品男人的天堂亚洲| 久久天堂一区二区三区四区| 19禁男女啪啪无遮挡网站| 久久久国产精品麻豆| 亚洲一码二码三码区别大吗| 老鸭窝网址在线观看| 母亲3免费完整高清在线观看| 大香蕉久久成人网| 超色免费av| 99久久精品国产亚洲精品| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品一区三区| 一级毛片女人18水好多| 久久久久久久精品精品| 精品国产国语对白av| 一级a爱视频在线免费观看| 久久久久久久国产电影| 精品一区二区三卡| 亚洲成av片中文字幕在线观看| 色94色欧美一区二区| 最近最新免费中文字幕在线| 欧美精品人与动牲交sv欧美| 精品熟女少妇八av免费久了| 爱豆传媒免费全集在线观看| 三级毛片av免费| 法律面前人人平等表现在哪些方面 | 亚洲欧美精品综合一区二区三区| 欧美日韩黄片免| 亚洲欧美色中文字幕在线| 人人妻人人爽人人添夜夜欢视频| 国产有黄有色有爽视频| 亚洲精品日韩在线中文字幕| 亚洲天堂av无毛| 亚洲国产欧美在线一区| 午夜免费鲁丝| 国产精品熟女久久久久浪| 欧美亚洲日本最大视频资源| 亚洲专区中文字幕在线| 欧美变态另类bdsm刘玥| 国产在线免费精品| 免费观看人在逋| 国产精品偷伦视频观看了| 老司机午夜十八禁免费视频| 黑人欧美特级aaaaaa片| 搡老熟女国产l中国老女人| 欧美日韩av久久| 老司机靠b影院| 日韩视频在线欧美| 国产不卡av网站在线观看| av有码第一页| 老鸭窝网址在线观看| 视频区图区小说| 丰满人妻熟妇乱又伦精品不卡| 精品卡一卡二卡四卡免费| 久久久久国产精品人妻一区二区| 国产精品影院久久| 久久天堂一区二区三区四区| 人人妻人人澡人人爽人人夜夜| 超色免费av| 777米奇影视久久| 高潮久久久久久久久久久不卡| 老熟妇乱子伦视频在线观看 | 啦啦啦 在线观看视频| 高清黄色对白视频在线免费看| 99香蕉大伊视频| av一本久久久久| 久久久久久久精品精品| 老熟妇乱子伦视频在线观看 | 日本a在线网址| 青草久久国产| 99精品久久久久人妻精品| 免费久久久久久久精品成人欧美视频| 国产高清videossex| 国产在视频线精品| 在线亚洲精品国产二区图片欧美| 日韩欧美一区视频在线观看| 国产一区二区三区综合在线观看| 大型av网站在线播放| 人人妻人人澡人人看| 法律面前人人平等表现在哪些方面 | 老熟女久久久| 美女高潮喷水抽搐中文字幕| 国产伦人伦偷精品视频| 麻豆国产av国片精品| 老司机福利观看| 国精品久久久久久国模美| 色婷婷久久久亚洲欧美| 久久久久久人人人人人| 国产真人三级小视频在线观看| www.av在线官网国产| 久久久久久久久久久久大奶| 一区二区三区激情视频| 无遮挡黄片免费观看| 久久亚洲国产成人精品v| 国产男女超爽视频在线观看| 亚洲av日韩在线播放| tocl精华| 女警被强在线播放| 9色porny在线观看| 欧美国产精品va在线观看不卡| 悠悠久久av| 久久久精品免费免费高清| 啦啦啦 在线观看视频| 免费看十八禁软件| 日韩中文字幕视频在线看片| 啦啦啦中文免费视频观看日本| 一二三四在线观看免费中文在| 极品少妇高潮喷水抽搐| 亚洲,欧美精品.| 亚洲七黄色美女视频| 视频区图区小说| 黄色怎么调成土黄色| 国产精品一区二区精品视频观看| 日韩制服丝袜自拍偷拍| 亚洲男人天堂网一区| 最近最新免费中文字幕在线| 国产97色在线日韩免费| 丁香六月欧美| 欧美成人午夜精品| 国产在线观看jvid| 精品少妇黑人巨大在线播放| 亚洲激情五月婷婷啪啪| 老司机福利观看| 国产免费一区二区三区四区乱码| 老熟妇仑乱视频hdxx| 人妻一区二区av| 又大又爽又粗| 亚洲精品国产一区二区精华液| 我要看黄色一级片免费的| 成年av动漫网址| 亚洲久久久国产精品| 狠狠婷婷综合久久久久久88av| 热99re8久久精品国产| 两个人看的免费小视频| 曰老女人黄片| 波多野结衣av一区二区av| 女人久久www免费人成看片| 亚洲熟女精品中文字幕| 欧美另类一区| 欧美激情久久久久久爽电影 | 欧美成人午夜精品| 久久天堂一区二区三区四区| 1024视频免费在线观看| 中文精品一卡2卡3卡4更新| 久久午夜综合久久蜜桃| 水蜜桃什么品种好| 国产一区有黄有色的免费视频| 国产精品一区二区在线观看99| 如日韩欧美国产精品一区二区三区| 成人18禁高潮啪啪吃奶动态图| 亚洲va日本ⅴa欧美va伊人久久 | 欧美精品一区二区大全| 99热网站在线观看| 在线观看免费高清a一片| 国产无遮挡羞羞视频在线观看| 热99久久久久精品小说推荐| av网站免费在线观看视频| 操美女的视频在线观看| 91精品伊人久久大香线蕉| 天天躁夜夜躁狠狠躁躁| 日本精品一区二区三区蜜桃| 国产色视频综合| 老司机午夜十八禁免费视频| av天堂久久9| 亚洲熟女毛片儿| 18禁国产床啪视频网站| 国产精品偷伦视频观看了| 搡老乐熟女国产| 久久国产精品影院| 极品少妇高潮喷水抽搐| 19禁男女啪啪无遮挡网站| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久久免费视频了| 一级a爱视频在线免费观看| 亚洲欧美色中文字幕在线| 久久人妻熟女aⅴ| 欧美一级毛片孕妇| av电影中文网址| 日本五十路高清| 大香蕉久久网| 国产欧美日韩综合在线一区二区| 少妇的丰满在线观看| 精品第一国产精品| 亚洲va日本ⅴa欧美va伊人久久 | 另类亚洲欧美激情| 久久久久久久久久久久大奶| 国产亚洲午夜精品一区二区久久| 精品福利观看| 欧美在线黄色| 国产精品久久久久久人妻精品电影 | 国产精品 国内视频| 国产激情久久老熟女| 午夜激情av网站| 久久天堂一区二区三区四区| 久久精品成人免费网站| svipshipincom国产片| 亚洲精品成人av观看孕妇| 美女中出高潮动态图| 亚洲国产精品一区二区三区在线| 大码成人一级视频| 老司机在亚洲福利影院| 久热这里只有精品99| 久久天躁狠狠躁夜夜2o2o| 欧美国产精品一级二级三级| 国产精品久久久久成人av| 激情视频va一区二区三区| 男女床上黄色一级片免费看| 两性夫妻黄色片| 国产av国产精品国产| 少妇精品久久久久久久| 色94色欧美一区二区| 人人妻人人澡人人看| 99九九在线精品视频| 欧美精品高潮呻吟av久久| 亚洲国产精品一区二区三区在线| 成人三级做爰电影| 欧美精品亚洲一区二区| 国产不卡av网站在线观看| 免费人妻精品一区二区三区视频| 久久影院123| 亚洲五月色婷婷综合| 人妻 亚洲 视频| a级毛片在线看网站| 国产片内射在线| 国产亚洲一区二区精品| 精品少妇黑人巨大在线播放| 国产熟女午夜一区二区三区| 国产日韩一区二区三区精品不卡| 亚洲熟女毛片儿| 国产在线视频一区二区| 国产精品秋霞免费鲁丝片| 美国免费a级毛片| 天天影视国产精品| 老鸭窝网址在线观看| 亚洲国产精品一区三区| 午夜视频精品福利| 亚洲视频免费观看视频| 如日韩欧美国产精品一区二区三区| 一级毛片精品| 嫩草影视91久久| 中文字幕人妻丝袜制服| 亚洲国产欧美在线一区| 日韩一区二区三区影片| 日韩欧美一区二区三区在线观看 | 韩国高清视频一区二区三区| 宅男免费午夜| 如日韩欧美国产精品一区二区三区| 精品第一国产精品| 国产欧美日韩一区二区三区在线| 亚洲av日韩精品久久久久久密| 可以免费在线观看a视频的电影网站| 欧美精品亚洲一区二区| 黄片小视频在线播放| 操出白浆在线播放| 美国免费a级毛片| 精品国产一区二区三区四区第35| 国产黄色免费在线视频| 人人澡人人妻人| 丝瓜视频免费看黄片| 久久精品国产亚洲av高清一级| 欧美中文综合在线视频| 国产成人欧美| 亚洲中文av在线| 欧美一级毛片孕妇| 久久国产精品影院| 欧美日韩中文字幕国产精品一区二区三区 | 每晚都被弄得嗷嗷叫到高潮| 色94色欧美一区二区| 婷婷丁香在线五月| 国产精品99久久99久久久不卡| 亚洲第一青青草原| 国产精品秋霞免费鲁丝片| 五月天丁香电影| 国产91精品成人一区二区三区 | 一级片'在线观看视频| 国产精品国产三级国产专区5o| 亚洲精品美女久久av网站| 精品一区二区三区四区五区乱码| 亚洲欧美一区二区三区久久| 国产一区有黄有色的免费视频| 久久久久国内视频| 一级毛片女人18水好多| 男女午夜视频在线观看| 交换朋友夫妻互换小说| 欧美另类亚洲清纯唯美| 热re99久久精品国产66热6| 91成年电影在线观看| 亚洲av男天堂| 侵犯人妻中文字幕一二三四区| 国产精品九九99| 中文字幕色久视频| 亚洲熟女毛片儿| 成在线人永久免费视频| 亚洲欧洲日产国产| 女人被躁到高潮嗷嗷叫费观| 亚洲七黄色美女视频| 国产高清videossex| 狠狠精品人妻久久久久久综合| 丝袜美腿诱惑在线| 91大片在线观看| 丝瓜视频免费看黄片| 亚洲一码二码三码区别大吗| 黄色怎么调成土黄色| 丝袜喷水一区| 最新的欧美精品一区二区| 一级毛片精品| 精品视频人人做人人爽| 欧美人与性动交α欧美精品济南到| 一区二区三区乱码不卡18| 另类亚洲欧美激情| 99热国产这里只有精品6| av不卡在线播放| 捣出白浆h1v1| 亚洲精品乱久久久久久| 建设人人有责人人尽责人人享有的| 人人妻人人澡人人爽人人夜夜| 中文精品一卡2卡3卡4更新| 国产精品久久久久久精品古装| 性色av乱码一区二区三区2| av天堂久久9| av天堂在线播放| 18在线观看网站| 夜夜骑夜夜射夜夜干| 热99久久久久精品小说推荐| 视频区图区小说| 欧美精品啪啪一区二区三区 | 人人妻人人澡人人看| 五月天丁香电影| 国产欧美亚洲国产| 欧美日韩亚洲国产一区二区在线观看 | 国产不卡av网站在线观看| 三级毛片av免费| 侵犯人妻中文字幕一二三四区| 国产男女超爽视频在线观看| 成人影院久久| 日本五十路高清| 欧美亚洲 丝袜 人妻 在线| av电影中文网址| 亚洲av美国av| 色综合欧美亚洲国产小说| 午夜视频精品福利| 啪啪无遮挡十八禁网站| 中文字幕av电影在线播放| 国产xxxxx性猛交| 日韩有码中文字幕| 国产精品久久久久久精品电影小说| 亚洲av美国av| 99久久人妻综合| 伊人亚洲综合成人网| 真人做人爱边吃奶动态| 国产又爽黄色视频| 在线观看免费午夜福利视频| a级毛片黄视频| 人人妻人人澡人人爽人人夜夜| 正在播放国产对白刺激| 成人国产一区最新在线观看| 一级,二级,三级黄色视频| 一区二区三区四区激情视频| 久久精品熟女亚洲av麻豆精品| 国产国语露脸激情在线看| 国产精品一区二区精品视频观看| 精品久久久久久久毛片微露脸 | 蜜桃国产av成人99| 永久免费av网站大全| 在线观看免费高清a一片| 久久毛片免费看一区二区三区| 亚洲欧美成人综合另类久久久| 久热爱精品视频在线9| 国内毛片毛片毛片毛片毛片| 青春草亚洲视频在线观看| 久久久久久久大尺度免费视频| 久久久水蜜桃国产精品网| 美女扒开内裤让男人捅视频| 国产精品影院久久| 国产色视频综合| 美女国产高潮福利片在线看| 国产亚洲欧美在线一区二区| 高清在线国产一区| 亚洲国产毛片av蜜桃av| 制服诱惑二区| 曰老女人黄片| 俄罗斯特黄特色一大片| 亚洲精品久久午夜乱码| 久热这里只有精品99| 日本wwww免费看| 久久久久国产精品人妻一区二区| 国产免费福利视频在线观看| 国产成人免费无遮挡视频| 欧美激情 高清一区二区三区| 中文精品一卡2卡3卡4更新| 91字幕亚洲| 欧美精品一区二区免费开放| 一边摸一边抽搐一进一出视频| 亚洲av电影在线进入| 欧美久久黑人一区二区| 波多野结衣av一区二区av| 国产精品秋霞免费鲁丝片| 黑人巨大精品欧美一区二区mp4| 51午夜福利影视在线观看| 黄色怎么调成土黄色| 国产又色又爽无遮挡免| av网站在线播放免费| 亚洲少妇的诱惑av| 一区在线观看完整版| 秋霞在线观看毛片| 亚洲av日韩在线播放| 色播在线永久视频| 亚洲男人天堂网一区| 国产成人啪精品午夜网站| 亚洲国产毛片av蜜桃av| 亚洲欧美成人综合另类久久久| 国产精品久久久久久精品古装| 中国美女看黄片| 国产老妇伦熟女老妇高清| 国产精品自产拍在线观看55亚洲 | 精品免费久久久久久久清纯 | 国产国语露脸激情在线看| 国产精品久久久久久精品古装| 色精品久久人妻99蜜桃| 欧美国产精品va在线观看不卡| 国产精品自产拍在线观看55亚洲 | 成人黄色视频免费在线看| 亚洲欧美日韩另类电影网站| 99国产精品一区二区蜜桃av | 国产极品粉嫩免费观看在线| 超碰97精品在线观看| 亚洲av国产av综合av卡| 777米奇影视久久| 男人爽女人下面视频在线观看| 免费高清在线观看日韩| 丁香六月欧美| 另类亚洲欧美激情| 国产精品免费大片| 日本91视频免费播放| 精品久久久久久电影网| 欧美亚洲日本最大视频资源| 亚洲精品久久久久久婷婷小说| 99九九在线精品视频| 日韩熟女老妇一区二区性免费视频| 国产成人av激情在线播放| 久久亚洲国产成人精品v| 午夜成年电影在线免费观看| 美女中出高潮动态图| 老司机午夜福利在线观看视频 | 久久久久久久大尺度免费视频| 久久九九热精品免费| 在线 av 中文字幕| 国产免费福利视频在线观看| 一本综合久久免费| 99精国产麻豆久久婷婷| 91精品伊人久久大香线蕉| 亚洲精品粉嫩美女一区| 久久久久久久久免费视频了| 蜜桃在线观看..| www.熟女人妻精品国产| 美女午夜性视频免费| 国产成人欧美在线观看 | 亚洲中文日韩欧美视频| 19禁男女啪啪无遮挡网站| 国产99久久九九免费精品| 久久亚洲国产成人精品v| 男女国产视频网站| 午夜福利乱码中文字幕| 美女扒开内裤让男人捅视频| 国产精品香港三级国产av潘金莲| 成人国产av品久久久| 少妇粗大呻吟视频| 亚洲精品自拍成人| 黑人巨大精品欧美一区二区mp4| 视频在线观看一区二区三区| 99热全是精品| 老鸭窝网址在线观看| 国产精品1区2区在线观看. | 手机成人av网站| 母亲3免费完整高清在线观看| 成年人午夜在线观看视频| 99国产精品一区二区蜜桃av | 丰满少妇做爰视频| 岛国毛片在线播放| 国产99久久九九免费精品| 啦啦啦 在线观看视频| 午夜福利视频精品| 亚洲精品粉嫩美女一区| 国产成人精品久久二区二区免费| 日韩熟女老妇一区二区性免费视频| 丝袜美足系列| 黑丝袜美女国产一区| 亚洲va日本ⅴa欧美va伊人久久 | 国产精品一区二区免费欧美 | 国产av国产精品国产| 国产免费现黄频在线看| 女人爽到高潮嗷嗷叫在线视频| 天堂中文最新版在线下载| 中文字幕精品免费在线观看视频| 男女免费视频国产| 精品福利观看| 麻豆av在线久日| 久热爱精品视频在线9| 国产成人精品久久二区二区91| 中文字幕精品免费在线观看视频| 亚洲性夜色夜夜综合| 国产免费视频播放在线视频| 黄色视频不卡| 真人做人爱边吃奶动态| 亚洲人成77777在线视频| 在线观看舔阴道视频| 中文字幕另类日韩欧美亚洲嫩草| 日韩免费高清中文字幕av| 日韩视频一区二区在线观看| 丰满少妇做爰视频| 黄频高清免费视频| 美女国产高潮福利片在线看| 久久免费观看电影| 久久人人爽人人片av| 久久九九热精品免费| 天堂俺去俺来也www色官网| 亚洲熟女精品中文字幕| 黑人巨大精品欧美一区二区mp4| 久久性视频一级片| 亚洲成人免费电影在线观看| 中文欧美无线码| 国产日韩欧美视频二区| 亚洲美女黄色视频免费看| 午夜视频精品福利| 亚洲伊人久久精品综合| 无遮挡黄片免费观看| 大片电影免费在线观看免费| 国产av一区二区精品久久| 久久午夜综合久久蜜桃| 国产一卡二卡三卡精品| 久久午夜综合久久蜜桃| 亚洲精品美女久久av网站| 亚洲av片天天在线观看| 精品免费久久久久久久清纯 | 操美女的视频在线观看| 亚洲成人免费av在线播放| 欧美另类亚洲清纯唯美| 久久天堂一区二区三区四区| 亚洲欧美日韩高清在线视频 | av网站免费在线观看视频| 亚洲五月色婷婷综合| 久久久国产成人免费| 久久精品国产a三级三级三级| av片东京热男人的天堂| 欧美av亚洲av综合av国产av| 欧美中文综合在线视频| 国产又色又爽无遮挡免| 淫妇啪啪啪对白视频 | 精品久久蜜臀av无| 久久亚洲精品不卡| 黄片播放在线免费| 女人久久www免费人成看片| 欧美人与性动交α欧美软件| 69av精品久久久久久 | av线在线观看网站| 欧美黄色淫秽网站| 一区福利在线观看| 丰满迷人的少妇在线观看| 美女高潮到喷水免费观看| 在线观看舔阴道视频| 国产福利在线免费观看视频| 国产精品国产三级国产专区5o| 久久精品久久久久久噜噜老黄| av在线老鸭窝| 国产亚洲欧美精品永久| 久久天堂一区二区三区四区| av在线播放精品| 国产伦人伦偷精品视频| 99热全是精品| 一区在线观看完整版| videos熟女内射| 多毛熟女@视频| 欧美日本中文国产一区发布| 欧美日韩国产mv在线观看视频| 高清黄色对白视频在线免费看| 在线观看www视频免费| 色婷婷av一区二区三区视频| 日韩欧美一区视频在线观看| 脱女人内裤的视频| 丝袜美腿诱惑在线| 少妇 在线观看| 别揉我奶头~嗯~啊~动态视频 | 欧美日韩亚洲国产一区二区在线观看 | 日本av免费视频播放| 极品少妇高潮喷水抽搐| 丰满人妻熟妇乱又伦精品不卡| 黄色视频不卡| 欧美精品啪啪一区二区三区 | 丰满少妇做爰视频| 久久天躁狠狠躁夜夜2o2o| 午夜精品国产一区二区电影| 精品欧美一区二区三区在线| a级片在线免费高清观看视频| 精品高清国产在线一区| 69av精品久久久久久 | 美女中出高潮动态图| 高潮久久久久久久久久久不卡| 香蕉丝袜av| 国产亚洲一区二区精品| 桃花免费在线播放| 欧美老熟妇乱子伦牲交| 亚洲av片天天在线观看|