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

    航空發(fā)動機數值仿真智能綜合集成平臺架構研究

    2022-08-11 00:41:50王佳川
    計算機集成制造系統(tǒng) 2022年7期
    關鍵詞:優(yōu)化

    黃 博,姚 燁,王佳川

    (1.西北工業(yè)大學 機電學院,陜西 西安 710072;2.中國航發(fā)商用航空發(fā)動機有限責任公司,上海 200241; 3.上海市航空發(fā)動機數字孿生重點實驗室,上海 200241)

    1 問題的提出

    航空發(fā)動機研發(fā)數字化是集成了設計、仿真分析、優(yōu)化驗證于一體的系統(tǒng)工程,體現了一個國家的工業(yè)基礎、裝備先進性和科技水平。數字化仿真技術作為支撐航空發(fā)動機自主研發(fā)的重要手段,可大幅提高研發(fā)效率、減少反復實物試驗、縮短研制周期、降低研制成本。近年來,航空發(fā)動機研發(fā)過程與核心要素日益復雜[1-3]、研發(fā)輔助工具種類繁多、跨部門研發(fā)協同頻繁、多學科知識融合交叉等給傳統(tǒng)的數字化研發(fā)體系帶來了諸多挑戰(zhàn)。為了突破傳統(tǒng)體系的局限性,大型高端裝備數字化研發(fā)設計環(huán)境正向著集成化、平臺化、智能化的方向發(fā)展。

    國外大型工業(yè)軟件及綜合集成平臺領域起步較早,形成了以美國、德國、法國為代表的巨頭工業(yè)軟件企業(yè)和產品[1]。美國Ansys推出的WorkBench[2]是一種設計仿真集成環(huán)境,提供開放的框架結構,將豐富的分析工具整合在一起,實現了分析過程與分析工具的高度統(tǒng)一和緊密結合,由于耦合度較高,目前只適配Ansys系列軟件。Isight[3]是一套過程集成、優(yōu)化設計和穩(wěn)健性設計的軟件平臺,可以將數字技術、推理技術和設計探索技術有效融合,實現各軟件間的數據流傳遞,完成多學科多領域聯合仿真和優(yōu)化流程。另外,德國西門子Simcenter[4]、法國達索SIMULIA Suite等平臺也實現了面向多學科分析軟件的集成與流程分析。典型的科學計算工作流系統(tǒng)Kepler[5],是一套開源的科學數據分析與建模軟件應用平臺,通過可視化模式設計科學工作流并支持基于網格的分布式計算,但其針對具體的仿真流程建模過程需要花費大量時間[6]。

    隨著高端裝備數字研發(fā)能力提升及工業(yè)軟件自主化需求迫切,國內科研院所和企業(yè)開始投入到自主工業(yè)軟件集成平臺的研發(fā)中。大連理工大學陳飆松教授團隊[7-8]的SiPESC平臺,實現了面向核心任務管理、系統(tǒng)集成、大規(guī)模計算工程數據庫的集成開發(fā)環(huán)境,并完成了大量計算分析軟件集成案例。安世亞太研發(fā)的多學科異構集成平臺軟件SimCube[9],支持自主及商業(yè)計算分析軟件的聯合仿真、性能優(yōu)化、界面集成、結果可視化等功能。中國工程物理研究院針對多領域數值模擬軟件的封裝、計算任務的集成、異構資源屏蔽、圖形化交互等問題設計了數值模擬工作流管理平臺HSWAP[10],并將其應用在材料、光學、力學等領域,主要針對文件兼任務的松耦合。另外數值風洞、數值水池等研發(fā)團隊也開發(fā)了類似的集成平臺。

    綜上可見,國外仿真集成平臺產品在軟件成熟度和穩(wěn)定性上具有較強的優(yōu)勢,但大部分屬于通用軟件,缺少面向航空發(fā)動機專業(yè)特征的知識庫。國內的仿真集成平臺多數將重點放在自研求解器集成和流程集成方面,對廣泛商業(yè)軟件集成不夠深入。因此,現急需一套覆蓋航空發(fā)動機數字化研發(fā)設計全流程,以仿真數據為核心,流程靈活定義,經驗知識模板固化,支持多學科軟件、多物理場耦合的自動化智能化的航空發(fā)動機數值仿真智能綜合集成平臺(簡稱數值仿真集成平臺)。本文將重點開展航空發(fā)動機數值仿真智能綜合集成平臺架構的研究,為平臺的開發(fā)提供先進的方法理念與架構解決方案。研究思路如圖1所示。

    2 集成平臺需求分析

    數值仿真集成平臺貫通航空發(fā)動機的設計、計算分析和數據管理等環(huán)節(jié),其研制目標包括:實現各部門仿真工具的統(tǒng)一管理入口,系統(tǒng)化集成自研程序和商用軟件;跨部門仿真流程的自動化,不同學科仿真軟件間輸入輸出數據聯通流轉及流程自定義;仿真專家經驗數據固化積累,異構仿真過程數據有序管理以及面向型號/部件的仿真模板定義。

    圍繞數值仿真集成平臺建設目標,提升跨部門、多學科協同仿真效率,總結集成平臺架構需要滿足的6點需求分別是:①按照仿真場景定制仿真流程需求,支持可視化建模方式構建仿真流程和仿真模板參數化配置,形成自主的仿真流程規(guī)則表達語言;②靈活擴展仿真應用軟件需求,將常用的仿真軟件封裝成組件,同時支持用戶方便擴展新增商業(yè)、自研的仿真軟件;③全流程仿真數據管理需求,基于平臺的數據管理既要滿足仿真數據結構化管理,又要滿足仿真流程及模板的有組織管理;④仿真數據挖掘分析的需求,面向仿真優(yōu)化需求提供常用的多目標優(yōu)化算法庫,同時集成常用的機器學習引擎,為后續(xù)數據挖掘分析提供平臺基礎;⑤充分運用本地與集群資源需求,平臺的部署方案支持本地部署運行,同時也要支持基于集群/超算的遠程部署;⑥自研/商業(yè)軟件多物理場耦合需求,支持多學科軟件弱耦合的多物理場耦合需求,實現耦合界面數據映射算法,多求解器耦合策略,耦合節(jié)點間通信和基于適配器求解器封裝。

    3 平臺架構方法研究

    3.1 基于科學工作流仿真流程自定義方法

    航空發(fā)動機研發(fā)設計過程涉及CAD、CAE、分析優(yōu)化等相關軟件,包括整機與部件仿真涵蓋結構、流體、燃燒等多學科復雜分析,不同分析優(yōu)化目標對流程與關鍵參數的要求均有不同,本文采用科學工作流引擎實現靈活、高可擴展、可視化交互的仿真流程定義?;诳茖W工作流的仿真流程自定義方法主要包括兩方面:

    (1)可視化建模語言體系 是科學計算工作流中要素的規(guī)則體系,需要定義語素(語言符號)、語法和語義,并實現相應的解析器(相當于常規(guī)編程語言的編譯)。集成軟件(計算模型)在該建模語言中抽象為一個可配置參數的圖形組件,圖形組件之間按規(guī)則(計算步驟)組織起來,便構成科學計算工作流,如圖2所示。其中,圖形組件間是否能夠組合與連接,由圖形組件的參數決定,而這些參數都是標準化的接口。每一個圖形元素即為工作流中的節(jié)點,每個節(jié)點可依據其功能設置不同的屬性,通過輸入輸出參數來體現,參數遵循統(tǒng)一的標準,支持與其他圖元節(jié)點的數據通訊。

    (2)可視化流程建模引擎 基于統(tǒng)一的建模語言開展工作流建模引擎開發(fā),相當于C++、Java等高級編程語言的集成開發(fā)環(huán)境(Integrated Development Environment, IDE),它由工作流建模環(huán)境、工作流模型庫、模型管理器3個主要部分構成。工作流建模環(huán)境為用戶提供建模過程的支持,模型管理器提供管理模型庫與工作流模型模板的管理功能,工作流模型庫保存與管理計算模型和工作流模型模板的元數據與其對應的模型文件。模塊組織關系如圖3所示。

    3.2 插件式集成軟件擴展機制

    在工作流引擎的基礎上,采用組件封裝和插件擴展的軟件設計技術,可提高航空發(fā)動機研發(fā)CAx工具軟件的復用,集成平臺中核心模塊均采用插件開發(fā)實現,使得軟件體系結構具有較強的開放性和集成性,有利于系統(tǒng)應用功能的動態(tài)擴展和系統(tǒng)核心功能的維護。

    插件的本質是動態(tài)鏈接庫,從軟件的結構層級來說,高于所有的軟件功能模塊,因此插件可以調用任意的軟件接口,復雜的功能拓展主要通過多態(tài)實現。插件式集成軟件擴展機制主要包括以下3個方面:

    (1)插件開放接口規(guī)則設計 充分考慮異構CAx軟件系統(tǒng)在仿真流程中的作用、功能特性、資源需求、數據需求等,研究并制定動態(tài)插拔的異構插件規(guī)則體系,并編制相應的應用程序接口(Applicaition Programming Interface, API)。

    (2)軟件集成模板設計器 仿真軟件的插件集成需要考慮軟件的輸入參數配置和軟件驅動運行。即使同一個軟件在不同流程中也需要不同的參數,因此需要一套簡單易用的模板設計器,以可視化方式配置所集成軟件的參數界面,并給出參數化模板生成對應的軟件運行腳本,支持目標軟件靜默執(zhí)行。

    (3)功能插件配置管理器 插件支持動態(tài)的加載與卸載,支持插件屬性配置,以及運行環(huán)境配置,方便實現具有基礎共性功能的插件共享。插件支持通過Web進行維護和升級,插件更新以及新增時,支持插件相關文件的上傳/下載功能。插件管理機制如圖4所示。

    3.3 跨部件與多學科耦合封裝技術

    目前多物理場耦合方法分為強耦合(直接耦合式解法)和弱耦合(分離解法),強耦合方程的理論推導過程相當復雜,由于不同的物理場中的物質屬性不同,強耦合方法求解時矩陣通常是病態(tài)的,這導致無法或很難得到問題正確解,不適用于成熟的商業(yè)軟件之間進行耦合。然而弱耦合可以有效地避免這些問題,可以利用現有的單一物理場求解器/軟件,通過耦合中間件模塊組合在一起,求解多物理場問題。平臺集成的商業(yè)和自研仿真軟件需要按照一定的封裝規(guī)則開發(fā)適配器,接入到耦合中間件模塊。耦合中間件架構如圖5所示。

    集成平臺借助耦合中間件實現多物理場弱耦合,其所需研究的關鍵技術如下:

    (1)耦合界面數據映射算法 求解器中的數據通常是以離散化網格來表征的,網絡通常由節(jié)點和面組成。位移、速度或力數據值通常與節(jié)點相關聯,而壓力與面相關聯。弱耦合涉及到兩個求解器在耦合界面上的數據匹配問題,即兩個不同的求解器進行耦合計算時,由于劃分網格的不同,會導致數據不能直接傳遞到另一個求解器,這就需要對數據進行映射處理??梢愿鶕詈锨蠼馄骶W格的特征,采取最近鄰映射或者最近投影映射法。

    (2)多求解器耦合策略 耦合策略分為串行顯示、串行隱式、并行顯式及并行隱式方案。串行是指一個參與者一個接一個地交錯執(zhí)行,并行是指兩個求解器同時執(zhí)行。使用顯式方案,兩個參與者每個時間步只運行一次。使用隱式方案,求解器會不斷迭代,直到滿足收斂標準。迭代方法一般在高斯—塞德爾(Gauss-Seidel)迭代的基礎上,采用各種加速方法進行迭代算法。

    (3)耦合節(jié)點間通信 將需要耦合的求解器劃分為一個整體的模塊,將對應求解器的適配器與集成軟件連接,通過映射模塊,在不同求解器之間基于Socket或信息傳遞接口(Message Passing Interface, MPI)進行點對點通訊。模塊化通訊使用耦合策略控制不同求解器的運行順序,在求解器之間不斷進行網格數據的通訊。將要耦合的物理場的耦合交界面網格劃分為不同的域,分配給通訊模塊的不同進程組,可以直接在進程組內不同進程之間進行通訊,也可以在不同進程組的進程之間進行通訊,將耦合交界面上的網格數據通過MPI和Socket通訊功能進行傳輸。通過映射模塊,使用插值方法,對不匹配的網格數據進行匹配。

    (4)基于適配器求解器封裝 因為不同求解器提供的接口不同,難以直接使用已有的集成平臺中的耦合功能進行分離耦合運算,所以需要開發(fā)一個連接集成平臺以及求解器的中間層,即適配器。適配器連接求解器和集成軟件,其中封裝了一些核心功能去調用耦合策略、映射方法以及通訊方式,不同求解器開發(fā)需要相應的適配器。通過適配器讀取配置文件中的信息來控制耦合,配置文件包括:耦合求解器想要交換的數據值;耦合交界面網格ID(可以通過適配器的功能定義網格ID);參與耦合的求解器以及物理場網格的名字;時間步長;映射方法;通訊方式;耦合策略等。

    3.4 異構仿真工程數據管理技術

    航空發(fā)動機仿真數據具有大規(guī)模、數據密集、類型多、動態(tài)存儲、高并發(fā)、格式不透明等特點。大部分商業(yè)軟件采用專有的工程項目數據庫管理過程數據,數值仿真集成平臺管理的數據包含商業(yè)軟件工程項目數據、仿真流程模板參數數據、仿真人員信息、仿真過程關鍵數據等。數據涉及業(yè)務邏輯、參與人員、多種軟件等,具有更高的復雜度。因此,需要統(tǒng)一的接口和管理方式將這些數據組織起來,方便后續(xù)復用與分析挖掘。

    (1)異構仿真數據的組織與管理 仿真模型、仿真結果等數據采用文件的形式進行管理,但是需要結構化數據將仿真對象、仿真人員、配置參數(環(huán)境參數、載荷信息、邊界條件等)這些數據聯系起來,形成清晰的組織結構。采用文件數據庫和結構化數據庫混合模式,設計樹狀組織管理結構,可以與產品物料清單(Bill of Material, BOM)聯系起來,便于擴展和查找。結構化數據庫可以保障仿真數據各要素一致完整,并支持不同粒度的數據追溯。文件數據庫通過索引方便高效地管理大規(guī)模數據/模型(TB級),支持跨文件、跨分區(qū)、跨磁盤和跨節(jié)點等擴展方案。

    (2)仿真任務分析過程數據管理 仿真過程數據不同于仿真數據,具有輸入動態(tài)性、模板依賴性和強干預性等特征?;谄脚_的仿真過程需要維護仿真流程定義數據和引用的仿真模板數據。如模板參數無法滿足仿真節(jié)點設置需求,需要啟動軟件手動配置參數,管理該節(jié)點軟件對應的工程數據。采用元數據標簽法,實現同一型號不同人員、不同工況的多維關系網,最大化經驗知識、模板的復用。

    3.5 本地與遠程混合環(huán)境部署架構

    設計分析人員在開展仿真分析工作時會有很多場景,可能是本地運行小的算例開展一些驗證工作,也可能是采用某個軟件開展單學科性能分析任務,亦或是多人協同完成面向具體型號的多學科聯合仿真性能分析。數值仿真集成平臺是各種仿真軟件的集中入口,應該提供靈活的部署方案,滿足本地、遠程、多人協同等應用場景仿真任務需求。如何在數值仿真智能綜合集成平臺有效結合和使用遠程超算的計算資源,并實現操作人員在本地辦公電腦中便捷、靈活的進行人機交互,是對兼顧網絡環(huán)境、應用模式、數據傳輸保障等因素的巨大挑戰(zhàn)。下面給出綜合集成平臺在混合環(huán)境下的部署方案:

    (1)基于超算遠程部署方案 需要將綜合平臺所有功能部件、集成的仿真軟件和仿真數據等完全部署到超算環(huán)境中,在超算環(huán)境中實現綜合集成平臺數據的組織管理、程序的運行管理和人機交互界面的渲染等。對于部署工作本身而言,這種部署方案簡單直接沒有額外的通信工作。在進行操作時,可以配套遠程作業(yè)管理系統(tǒng)進行使用,在本地辦公電腦中,通過作業(yè)管理系統(tǒng)提交帶用戶界面(User Interface, UI)的作業(yè),作業(yè)管理系統(tǒng)負責從超算中啟動綜合集成平臺,并能夠在本地使用諸如虛擬網絡控制臺(Virtual Network Console, VNC)客戶端等工具軟件,渲染并顯示綜合集成平臺UI界面,最終可以通過該界面實現人機交互,實現使用綜合集成平臺的目的。

    (2)本地與超算混合部署方案 要求綜合集成平臺客戶端(UI人機交互部分)和后臺功能部件分開部署,能夠實現兩者的數據通信。這種部署方式,對于仿真數據而言,需要兼顧本地和遠程兩部分,需要提供額外的數據同步操作,仿真組件庫、流程模板庫等完成數據同步后才能在多個用戶之間共享。而仿真軟件則可以將負責計算的部分部署在超算,以充分使用超算算力資源,建模和前處理則可以部署在本地辦公電腦,實現高效渲染和顯示。另外,對于超算集群而言,作為服務端需要開辟專門的對外網絡端口,這在具體實施中增加了政策風險。混合部署方案拓撲如圖6所示。

    3.6 仿真設計/數據智能優(yōu)化分析方法

    仿真分析的目的是為了支撐優(yōu)化設計,多目標優(yōu)化是在復雜裝備設計分析過程中常用的技術手段,是在多輸入條件中尋找一組非劣解集。實際工程問題的優(yōu)化主要是修改目標幾何模型,或修改仿真分析輸入參數條件,來得到滿足工程設計的最優(yōu)解。在集成平臺中集成常用的優(yōu)化方法,如BFGS算法,序列線性規(guī)劃(Sequential Linear Programming, SLP)算法,序列二次規(guī)劃(Sequence Quadratic Program, SQP)算法和遺傳算法(Genetic Algorithm, GA)等,方便在仿真流程中加入多目標優(yōu)化組件。另外,隨著仿真集成平臺數據的積累,大規(guī)模仿真過程和結果數據可支撐AI與CAE結合,在一定程度上減少了計算資源浪費,基于工作流規(guī)則和接口集成TensorFlow等AI引擎,為各專業(yè)學科數據分析提供軟硬件環(huán)境。集成平臺中要實現仿真設計/數據智能優(yōu)化分析方法,架構設計要考慮以下兩個方面:

    (1)多目標優(yōu)化算法集成 構建優(yōu)化算法庫,為集成平臺提供豐富的抽樣方法、優(yōu)化算法和代理模型等優(yōu)化資源,將算法封裝成組件。優(yōu)化算法組件作為仿真工作流中的邏輯控制節(jié)點,為CAD/CAE組件提供設計變量,優(yōu)化組件使用代理模型擬合采樣數據的變量和響應值之間的關系,然后使用優(yōu)化算法,基于構建的代理模型進行迭代尋優(yōu),可自動執(zhí)行迭代尋優(yōu)過程,直到迭代結束,最后輸出優(yōu)化結果。

    (2)機器學習引擎集成 基于集成平臺統(tǒng)一規(guī)則集成飛槳、TensorFlow、PyTorch、NumPy等機器學習引擎,基于插件模式開發(fā)機器學習組件,借助異構仿真工程數據管理模塊打通與平臺底層數據的接口。為上層業(yè)務智能分析提供建模環(huán)境和數據基礎,降低用戶構建AI環(huán)境和數據鏈接的成本,專注業(yè)務模型的開發(fā)。

    4 應用實例分析

    在平臺架構方法的指導下,下面給出航空發(fā)動機數值仿真智能綜合集成平臺設計實例。平臺采用分層架構模式及模塊化設計理念,有利于降低模塊之間的邏輯耦合度,便于協同研發(fā)與功能的迭代升級,為基于工作流的仿真集成環(huán)境后續(xù)功能擴展提供了良好的技術架構。

    本實例采用四層架構,分別是基礎設施層、數據層、核心層和交互層。按照軟件工程架構設計思想,將系統(tǒng)的功能和業(yè)務邏輯進行劃分,并在不同層的功能設計上落實本文提出的平臺架構方法,設計實例架構如圖7所示。

    (1)交互層 交互層是直接與用戶打交道的圖形界面。航空發(fā)動機數值仿真流程中涉及到的交互界面,包括仿真軟件組件封裝、仿真流程搭建、仿真模板管理、組件庫管理、仿真數據管理、仿真任務管理等功能模塊。

    (2)核心層 核心層是平臺的邏輯和控制核心,將基礎功能的底層控制機制、基礎功能封裝成核心層的中間件模塊,通過統(tǒng)一接口有機組裝起來。包括插件擴展封裝機制、仿真工作流引擎、多物理場耦合引擎、本地遠程部署機制、作業(yè)調度與管理、智能優(yōu)化引擎等。

    (3)數據層 數據是軟件平臺的血液,貫通所有模塊與業(yè)務,將各類數據統(tǒng)籌管理,方便功能模塊的共享,包括優(yōu)化算法庫、仿真案例庫、仿真過程數據庫、型號模型數據庫等。

    (4)基礎設施層 基礎設施層是軟件平臺運行的底層硬件基礎,同時在軟件架構設計時也要兼顧網絡拓撲、通訊帶寬的限制、硬件計算能力等。本實例考慮集群與本地的混合部署模式,因此,考慮集群、磁盤陣列、本地工作站等硬件資源。

    表1給出了本文應用實例平臺與商業(yè)軟件/開源軟件在架構功能方面的對比。由表1可知,Ansys公司的WorkBench和西門子公司的Simcenter平臺針對各自廠家的工具軟件接口支持性較好,但是在底層功能擴展和復雜計算環(huán)境支持方面還不完善;Isight和Kepler平臺偏向仿真流程管理,驅動軟件形式單一粒度較粗,針對異構仿真工程數據的管理缺少功能模塊和接口;而本文應用實例平臺在設計過程中完全兼顧和覆蓋了這些先進架構功能。

    表1 架構功能對比

    綜上可見,通過該實例充分落實了本文提出的6大平臺架構指導方法,按照軟件架構思想在不同層設計了相應的功能模塊,從平臺的先進性、完整性方面驗證了平臺架構方法的合理性和有效性。

    5 結束語

    為了滿足航空發(fā)動機數值化研發(fā)體系向集成化、智能化、平臺化方向發(fā)展,圍繞數值仿真集成平臺建設目標,提升跨部門、多學科協同仿真效率,本文總結了集成平臺架構需要滿足的六大需求。針對六大需求開展了平臺架構方法研究,并提出基于科學工作流仿真流程自定義方法、插件式集成軟件擴展機制、跨部件與多學科耦合封裝技術、異構仿真工程數據管理技術、本地與遠程混合環(huán)境部署架構、仿真設計/數據智能優(yōu)化分析方法六大方法指導平臺設計開發(fā)。最后,通過航空發(fā)動機數值仿真智能綜合集成平臺設計實例與主流軟件平臺在先進功能方面的對比分析,驗證了提出方法的先進性、合理性和有效性。本文重點開展了航空發(fā)動機數值化研發(fā)集成平臺架構方法的研究,未來將進一步在集成平臺實施關鍵技術和體系化研發(fā)方法方面開展研究。

    猜你喜歡
    優(yōu)化
    超限高層建筑結構設計與優(yōu)化思考
    房地產導刊(2022年5期)2022-06-01 06:20:14
    PEMFC流道的多目標優(yōu)化
    能源工程(2022年1期)2022-03-29 01:06:28
    民用建筑防煙排煙設計優(yōu)化探討
    關于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數”優(yōu)化運算——以2021年解析幾何高考題為例
    圍繞“地、業(yè)、人”優(yōu)化產業(yè)扶貧
    事業(yè)單位中固定資產會計處理的優(yōu)化
    消費導刊(2018年8期)2018-05-25 13:20:08
    4K HDR性能大幅度優(yōu)化 JVC DLA-X8 18 BC
    幾種常見的負載均衡算法的優(yōu)化
    電子制作(2017年20期)2017-04-26 06:57:45
    激情视频va一区二区三区| 亚洲国产高清在线一区二区三 | 美女大奶头视频| 99精品欧美一区二区三区四区| 精品免费久久久久久久清纯| 麻豆国产av国片精品| 国产精品永久免费网站| 免费在线观看视频国产中文字幕亚洲| 国产高清视频在线播放一区| 美女高潮喷水抽搐中文字幕| av视频在线观看入口| 十分钟在线观看高清视频www| 国产成人免费无遮挡视频| 精品国产亚洲在线| 日韩欧美一区二区三区在线观看| 中文字幕另类日韩欧美亚洲嫩草| 美女 人体艺术 gogo| 免费在线观看日本一区| 久久精品91蜜桃| 国产蜜桃级精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 国产麻豆69| 一级毛片精品| 每晚都被弄得嗷嗷叫到高潮| 国产精品免费一区二区三区在线| 国产精品久久视频播放| 免费久久久久久久精品成人欧美视频| 亚洲va日本ⅴa欧美va伊人久久| 涩涩av久久男人的天堂| 亚洲狠狠婷婷综合久久图片| 欧美色视频一区免费| 精品第一国产精品| 亚洲三区欧美一区| 一进一出好大好爽视频| 国产精品久久视频播放| 亚洲国产精品久久男人天堂| 岛国在线观看网站| 国产又爽黄色视频| 国产亚洲精品久久久久5区| 亚洲欧美激情综合另类| 亚洲国产日韩欧美精品在线观看 | 色综合站精品国产| av福利片在线| 脱女人内裤的视频| 老汉色∧v一级毛片| 国产成年人精品一区二区| 十八禁人妻一区二区| 多毛熟女@视频| 亚洲午夜精品一区,二区,三区| 成人三级黄色视频| 亚洲精品美女久久av网站| 久久精品成人免费网站| 日本在线视频免费播放| 成人国产综合亚洲| 99re在线观看精品视频| 他把我摸到了高潮在线观看| 欧美日韩亚洲综合一区二区三区_| 久久久久久久精品吃奶| 午夜老司机福利片| 亚洲av片天天在线观看| 精品人妻1区二区| 亚洲av电影在线进入| 久久久国产成人免费| 中文字幕久久专区| 少妇粗大呻吟视频| 国产欧美日韩精品亚洲av| 真人一进一出gif抽搐免费| 国内精品久久久久精免费| 久久久久久免费高清国产稀缺| 中文字幕人妻丝袜一区二区| 香蕉国产在线看| 午夜免费鲁丝| 18禁观看日本| av有码第一页| 国产成人av教育| 一区二区三区国产精品乱码| 12—13女人毛片做爰片一| 午夜福利,免费看| 午夜福利视频1000在线观看 | 99在线人妻在线中文字幕| 99在线视频只有这里精品首页| 免费观看精品视频网站| 国产精品综合久久久久久久免费 | 国产精品一区二区在线不卡| 精品少妇一区二区三区视频日本电影| avwww免费| 久久天躁狠狠躁夜夜2o2o| 正在播放国产对白刺激| 中文亚洲av片在线观看爽| 88av欧美| 欧美成狂野欧美在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲中文字幕日韩| 免费高清视频大片| 麻豆国产av国片精品| 真人做人爱边吃奶动态| 黑人巨大精品欧美一区二区蜜桃| 国内毛片毛片毛片毛片毛片| 成年女人毛片免费观看观看9| 岛国在线观看网站| 日韩欧美一区二区三区在线观看| 在线观看免费日韩欧美大片| 99精品在免费线老司机午夜| 日韩欧美国产一区二区入口| 少妇裸体淫交视频免费看高清 | 啪啪无遮挡十八禁网站| 国产精品一区二区在线不卡| 国产成人欧美| 女生性感内裤真人,穿戴方法视频| 一区在线观看完整版| 免费一级毛片在线播放高清视频 | 欧美成人免费av一区二区三区| 激情在线观看视频在线高清| 午夜福利一区二区在线看| 一本久久中文字幕| 可以在线观看毛片的网站| 此物有八面人人有两片| 精品国产乱码久久久久久男人| 一卡2卡三卡四卡精品乱码亚洲| 国产精品98久久久久久宅男小说| 91麻豆精品激情在线观看国产| 久久午夜亚洲精品久久| 黄色毛片三级朝国网站| 久久久久亚洲av毛片大全| 免费女性裸体啪啪无遮挡网站| 午夜久久久在线观看| 国产亚洲精品久久久久久毛片| 深夜精品福利| 国产成人av教育| 又黄又粗又硬又大视频| 精品人妻在线不人妻| 亚洲五月婷婷丁香| 日日夜夜操网爽| 最好的美女福利视频网| 久久久久久久精品吃奶| 成人18禁高潮啪啪吃奶动态图| 999精品在线视频| 91成年电影在线观看| 首页视频小说图片口味搜索| 午夜激情av网站| 亚洲av美国av| 女警被强在线播放| 免费在线观看视频国产中文字幕亚洲| 在线永久观看黄色视频| 禁无遮挡网站| 91精品国产国语对白视频| 亚洲专区中文字幕在线| 黄频高清免费视频| 黄色丝袜av网址大全| 12—13女人毛片做爰片一| 丝袜美腿诱惑在线| 一卡2卡三卡四卡精品乱码亚洲| 日韩中文字幕欧美一区二区| 琪琪午夜伦伦电影理论片6080| 女人被狂操c到高潮| 国产精品一区二区精品视频观看| 宅男免费午夜| 国产国语露脸激情在线看| 欧美绝顶高潮抽搐喷水| 精品国产美女av久久久久小说| 69精品国产乱码久久久| 久久精品亚洲精品国产色婷小说| 丰满的人妻完整版| 午夜精品久久久久久毛片777| 国产亚洲精品综合一区在线观看 | 亚洲少妇的诱惑av| 久久亚洲真实| 国产精品电影一区二区三区| 美女免费视频网站| 久久性视频一级片| 91在线观看av| 丰满的人妻完整版| 国产成人av教育| 美女国产高潮福利片在线看| 亚洲av美国av| 最新在线观看一区二区三区| 一边摸一边做爽爽视频免费| 美女高潮到喷水免费观看| 首页视频小说图片口味搜索| 中文亚洲av片在线观看爽| 国产精品爽爽va在线观看网站 | 99精品在免费线老司机午夜| 搡老岳熟女国产| 国产色视频综合| 精品国产国语对白av| 日韩精品免费视频一区二区三区| 人妻久久中文字幕网| 欧美日韩黄片免| 叶爱在线成人免费视频播放| 国产激情久久老熟女| 麻豆久久精品国产亚洲av| 久久人妻福利社区极品人妻图片| 91大片在线观看| 黑人巨大精品欧美一区二区蜜桃| 最近最新免费中文字幕在线| 夜夜躁狠狠躁天天躁| 欧美+亚洲+日韩+国产| 亚洲aⅴ乱码一区二区在线播放 | av天堂久久9| 999久久久国产精品视频| 亚洲七黄色美女视频| 人人澡人人妻人| 无遮挡黄片免费观看| 一本久久中文字幕| 国产精品亚洲美女久久久| 亚洲,欧美精品.| 免费在线观看黄色视频的| 岛国在线观看网站| 欧美激情极品国产一区二区三区| 亚洲avbb在线观看| 国产高清视频在线播放一区| 啪啪无遮挡十八禁网站| 真人做人爱边吃奶动态| 国产av一区二区精品久久| 人人妻人人澡人人看| www日本在线高清视频| 男女床上黄色一级片免费看| 可以免费在线观看a视频的电影网站| 亚洲 欧美 日韩 在线 免费| 久久中文字幕一级| 我的亚洲天堂| 搡老岳熟女国产| 亚洲一卡2卡3卡4卡5卡精品中文| 成人国产综合亚洲| 露出奶头的视频| 正在播放国产对白刺激| 色在线成人网| 欧美色视频一区免费| 色老头精品视频在线观看| 日韩三级视频一区二区三区| 国产精品久久电影中文字幕| 成人国产一区最新在线观看| 视频在线观看一区二区三区| 中文字幕高清在线视频| 一进一出好大好爽视频| 欧美丝袜亚洲另类 | 极品人妻少妇av视频| 国产视频一区二区在线看| 国产精品av久久久久免费| 国产日韩一区二区三区精品不卡| 亚洲国产精品999在线| 欧美色视频一区免费| 国产精品国产高清国产av| 久久香蕉国产精品| 日本在线视频免费播放| 91成人精品电影| 两性夫妻黄色片| 亚洲va日本ⅴa欧美va伊人久久| 在线天堂中文资源库| svipshipincom国产片| 精品人妻1区二区| 两性夫妻黄色片| 又黄又爽又免费观看的视频| 91成年电影在线观看| 国产日韩一区二区三区精品不卡| 中文字幕人成人乱码亚洲影| 好男人在线观看高清免费视频 | 精品国产美女av久久久久小说| 两人在一起打扑克的视频| 亚洲精品在线美女| 日韩av在线大香蕉| 国产亚洲精品av在线| 日韩欧美三级三区| 国产亚洲欧美精品永久| 欧美+亚洲+日韩+国产| 男人的好看免费观看在线视频 | 久久久久久久久中文| 亚洲激情在线av| 中文字幕色久视频| 自线自在国产av| 亚洲五月婷婷丁香| 免费无遮挡裸体视频| 黑丝袜美女国产一区| 国内精品久久久久精免费| 男女做爰动态图高潮gif福利片 | 丝袜在线中文字幕| 麻豆av在线久日| 熟妇人妻久久中文字幕3abv| 女生性感内裤真人,穿戴方法视频| 正在播放国产对白刺激| 久久久国产欧美日韩av| bbb黄色大片| 夜夜看夜夜爽夜夜摸| 国产成人精品久久二区二区91| 中文字幕高清在线视频| 91精品三级在线观看| 久久伊人香网站| 国产成人系列免费观看| 精品国产一区二区久久| 妹子高潮喷水视频| 国产精品久久久久久人妻精品电影| 色播在线永久视频| 日韩三级视频一区二区三区| 九色国产91popny在线| 男女下面进入的视频免费午夜 | 亚洲天堂国产精品一区在线| 村上凉子中文字幕在线| 精品一区二区三区四区五区乱码| 亚洲国产欧美日韩在线播放| 又大又爽又粗| 国产乱人伦免费视频| 免费看美女性在线毛片视频| 国产成人精品久久二区二区免费| 欧美日本中文国产一区发布| 日韩欧美国产在线观看| 一区福利在线观看| 欧美丝袜亚洲另类 | 美女扒开内裤让男人捅视频| 久久国产乱子伦精品免费另类| 国产成人av教育| 色哟哟哟哟哟哟| 免费看十八禁软件| 午夜激情av网站| 久久精品国产亚洲av高清一级| 精品卡一卡二卡四卡免费| 在线天堂中文资源库| xxx96com| 1024香蕉在线观看| 午夜福利成人在线免费观看| 国产精品影院久久| 热99re8久久精品国产| 身体一侧抽搐| a在线观看视频网站| 97人妻天天添夜夜摸| 黄频高清免费视频| 国产伦一二天堂av在线观看| 亚洲狠狠婷婷综合久久图片| 成人精品一区二区免费| 亚洲精品国产精品久久久不卡| 久久人妻av系列| 日韩欧美一区视频在线观看| 亚洲男人的天堂狠狠| 美女大奶头视频| 久久久久久亚洲精品国产蜜桃av| 亚洲av成人不卡在线观看播放网| 精品熟女少妇八av免费久了| 久久婷婷成人综合色麻豆| 视频在线观看一区二区三区| 韩国av一区二区三区四区| 国产成人系列免费观看| 日韩av在线大香蕉| 免费一级毛片在线播放高清视频 | 亚洲第一欧美日韩一区二区三区| 搞女人的毛片| 色综合亚洲欧美另类图片| 在线观看日韩欧美| 欧美精品啪啪一区二区三区| 亚洲av美国av| 久久这里只有精品19| 午夜激情av网站| 亚洲av日韩精品久久久久久密| 日本三级黄在线观看| 国产一级毛片七仙女欲春2 | 18美女黄网站色大片免费观看| 成人欧美大片| 久热这里只有精品99| av片东京热男人的天堂| 亚洲人成网站在线播放欧美日韩| or卡值多少钱| av天堂久久9| 午夜影院日韩av| 精品欧美国产一区二区三| 国产免费男女视频| 一二三四社区在线视频社区8| 精品欧美一区二区三区在线| 免费在线观看影片大全网站| 日韩av在线大香蕉| 99国产精品99久久久久| 色综合亚洲欧美另类图片| 国内精品久久久久久久电影| 97碰自拍视频| 99久久精品国产亚洲精品| 婷婷精品国产亚洲av在线| 国产麻豆成人av免费视频| 日韩免费av在线播放| 国产精品美女特级片免费视频播放器 | 美女大奶头视频| www国产在线视频色| 成人三级黄色视频| 国产亚洲欧美精品永久| 一a级毛片在线观看| 欧美另类亚洲清纯唯美| 男女下面进入的视频免费午夜 | 国产激情欧美一区二区| 国产精品九九99| 99国产精品99久久久久| 亚洲狠狠婷婷综合久久图片| 一本综合久久免费| а√天堂www在线а√下载| 女性被躁到高潮视频| 亚洲七黄色美女视频| 又黄又爽又免费观看的视频| 日本 欧美在线| 少妇 在线观看| 亚洲五月天丁香| 亚洲欧美日韩另类电影网站| 国产色视频综合| 丝袜在线中文字幕| 国产成人av教育| 精品人妻1区二区| 成人国语在线视频| 久久国产精品男人的天堂亚洲| 亚洲av电影不卡..在线观看| 黄色女人牲交| 人人澡人人妻人| 久久久水蜜桃国产精品网| 看片在线看免费视频| 亚洲av电影在线进入| 美国免费a级毛片| а√天堂www在线а√下载| 手机成人av网站| 国产精品 国内视频| 国产精品亚洲一级av第二区| 如日韩欧美国产精品一区二区三区| 国产午夜福利久久久久久| 别揉我奶头~嗯~啊~动态视频| 久久人人97超碰香蕉20202| aaaaa片日本免费| 国产成人系列免费观看| 久久伊人香网站| 欧美午夜高清在线| 亚洲欧美日韩高清在线视频| 美女高潮到喷水免费观看| 91在线观看av| 大型av网站在线播放| 侵犯人妻中文字幕一二三四区| 91国产中文字幕| 久久久久久人人人人人| 夜夜夜夜夜久久久久| 搡老岳熟女国产| 久久伊人香网站| 亚洲av美国av| 热re99久久国产66热| 91av网站免费观看| 欧美日韩精品网址| 久久久久九九精品影院| 久久久久久人人人人人| 最新在线观看一区二区三区| 亚洲精品久久成人aⅴ小说| 午夜久久久久精精品| 国产精品一区二区三区四区久久 | 两性午夜刺激爽爽歪歪视频在线观看 | 成人国产一区最新在线观看| av片东京热男人的天堂| 亚洲精品久久成人aⅴ小说| 国产av一区在线观看免费| 在线观看www视频免费| 国产精品 国内视频| 欧美成人午夜精品| 国产国语露脸激情在线看| 日韩成人在线观看一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 欧美午夜高清在线| 91精品国产国语对白视频| 少妇裸体淫交视频免费看高清 | 法律面前人人平等表现在哪些方面| 国产精品免费一区二区三区在线| 男男h啪啪无遮挡| 免费看美女性在线毛片视频| 在线观看66精品国产| 麻豆国产av国片精品| 国产男靠女视频免费网站| 亚洲精品国产一区二区精华液| 黄片小视频在线播放| 中文字幕av电影在线播放| 免费观看人在逋| 国语自产精品视频在线第100页| 欧美乱妇无乱码| 日韩欧美国产一区二区入口| 熟女少妇亚洲综合色aaa.| 美女高潮喷水抽搐中文字幕| 国产亚洲精品av在线| 亚洲国产精品成人综合色| 日韩中文字幕欧美一区二区| 极品人妻少妇av视频| 国产精品香港三级国产av潘金莲| 国产成人av激情在线播放| 日韩视频一区二区在线观看| 老司机在亚洲福利影院| 69精品国产乱码久久久| 丁香欧美五月| 久久国产精品人妻蜜桃| 免费高清在线观看日韩| 亚洲欧美精品综合一区二区三区| 欧美色视频一区免费| 亚洲成人精品中文字幕电影| 久久午夜综合久久蜜桃| 亚洲狠狠婷婷综合久久图片| 一二三四在线观看免费中文在| 国产欧美日韩一区二区精品| 亚洲国产精品sss在线观看| 日韩欧美国产在线观看| 国产一区二区三区在线臀色熟女| 国产精品98久久久久久宅男小说| 久99久视频精品免费| 国产又爽黄色视频| 黄色a级毛片大全视频| 成在线人永久免费视频| 国产精品1区2区在线观看.| 久9热在线精品视频| 日韩成人在线观看一区二区三区| 精品第一国产精品| 欧美日韩亚洲综合一区二区三区_| 在线观看66精品国产| 国产国语露脸激情在线看| 久久狼人影院| 国产99久久九九免费精品| 日本免费a在线| 国产亚洲精品久久久久久毛片| 久久这里只有精品19| 精品人妻在线不人妻| 一级毛片精品| 久久久久久久精品吃奶| 国产精品综合久久久久久久免费 | 久久精品亚洲熟妇少妇任你| 91成人精品电影| 久久久久久国产a免费观看| 国产1区2区3区精品| 一二三四社区在线视频社区8| 亚洲精品粉嫩美女一区| 国产日韩一区二区三区精品不卡| 亚洲第一青青草原| 亚洲在线自拍视频| 19禁男女啪啪无遮挡网站| 男男h啪啪无遮挡| 国产三级黄色录像| 免费观看精品视频网站| 欧美日本亚洲视频在线播放| 操美女的视频在线观看| av中文乱码字幕在线| 亚洲一区高清亚洲精品| 欧美在线黄色| 亚洲国产中文字幕在线视频| 在线视频色国产色| 法律面前人人平等表现在哪些方面| 欧美黄色片欧美黄色片| 亚洲国产欧美网| 叶爱在线成人免费视频播放| 国产精品国产高清国产av| 中国美女看黄片| 757午夜福利合集在线观看| 日韩av在线大香蕉| 97人妻天天添夜夜摸| 国产精品久久视频播放| 欧美成人午夜精品| 青草久久国产| 一区二区三区高清视频在线| 黄色片一级片一级黄色片| 人人妻,人人澡人人爽秒播| 在线av久久热| 国产精品香港三级国产av潘金莲| 成人18禁在线播放| 久久中文看片网| 看免费av毛片| av有码第一页| 国产亚洲精品一区二区www| 午夜亚洲福利在线播放| 亚洲色图av天堂| 亚洲五月婷婷丁香| 女人被躁到高潮嗷嗷叫费观| 精品国产一区二区久久| 波多野结衣av一区二区av| 精品午夜福利视频在线观看一区| 午夜福利一区二区在线看| 精品第一国产精品| 99国产精品免费福利视频| 99久久99久久久精品蜜桃| 一本久久中文字幕| 国产精品久久久av美女十八| 免费观看精品视频网站| 久久亚洲真实| 国产男靠女视频免费网站| 久久久久久久久免费视频了| 满18在线观看网站| 韩国精品一区二区三区| 久久人妻av系列| 亚洲七黄色美女视频| 亚洲专区中文字幕在线| 午夜久久久久精精品| 午夜亚洲福利在线播放| 欧美一区二区精品小视频在线| 成人18禁在线播放| 香蕉丝袜av| 黄色 视频免费看| 亚洲国产中文字幕在线视频| 嫩草影院精品99| 不卡一级毛片| 国产成人啪精品午夜网站| 亚洲狠狠婷婷综合久久图片| 亚洲中文字幕一区二区三区有码在线看 | 黄色毛片三级朝国网站| 亚洲精品久久国产高清桃花| 成人三级黄色视频| 麻豆一二三区av精品| a在线观看视频网站| 人人妻人人澡欧美一区二区 | 久久亚洲精品不卡| 亚洲在线自拍视频| 国产精品精品国产色婷婷| 欧美激情久久久久久爽电影 | 999久久久精品免费观看国产| 精品国产一区二区久久| 日韩 欧美 亚洲 中文字幕| 国产aⅴ精品一区二区三区波| 久久久水蜜桃国产精品网| 欧美乱色亚洲激情| 欧美成人一区二区免费高清观看 | 婷婷六月久久综合丁香| 国产日韩一区二区三区精品不卡| av有码第一页| 岛国在线观看网站| 男人的好看免费观看在线视频 | 精品久久久久久久久久免费视频| 亚洲av日韩精品久久久久久密| 欧美精品啪啪一区二区三区|