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

    “天河一號”工業(yè)設(shè)計仿真云平臺中間件研究

    2015-10-29 03:32:20鄧子云劉楊兵肖久如
    中國機(jī)械工程 2015年6期
    關(guān)鍵詞:天河中間件超級計算機(jī)

    鄧子云 章 兢 劉楊兵 肖久如

    1.湖南大學(xué),長沙,410082  2.國家超級計算長沙中心,長沙,410082

    “天河一號”工業(yè)設(shè)計仿真云平臺中間件研究

    鄧子云1,2章兢1,2劉楊兵1,2肖久如1

    1.湖南大學(xué),長沙,4100822.國家超級計算長沙中心,長沙,410082

    在“天河一號”超級計算機(jī)上搭建了工業(yè)設(shè)計仿真云平臺,研發(fā)了該平臺的核心部件——工業(yè)設(shè)計仿真云平臺中間件,該中間件軟件集成了多款大型CAE軟件、SLURM作業(yè)調(diào)度軟件、License管理軟件、計算資源管理功能軟件。提出了工業(yè)設(shè)計仿真云平臺的總體架構(gòu),以及工業(yè)設(shè)計仿真云平臺中間件的設(shè)計思想、體系架構(gòu),解決了作業(yè)狀態(tài)轉(zhuǎn)換、作業(yè)提交等關(guān)鍵技術(shù)問題。使用工業(yè)設(shè)計仿真云平臺對火星著陸氣囊的碰撞分析進(jìn)行了仿真實(shí)驗(yàn),結(jié)果表明利用工業(yè)設(shè)計仿真云平臺中間件軟件可以簡便地使用“天河一號”超級計算機(jī)的計算資源。

    “天河一號”超級計算機(jī);工業(yè)設(shè)計仿真云中間件;體系架構(gòu);設(shè)計思想;關(guān)鍵技術(shù);仿真實(shí)驗(yàn)

    0 引言

    工業(yè)設(shè)計領(lǐng)域主要使用CAD、CAE軟件來進(jìn)行模擬仿真,以縮短產(chǎn)品研發(fā)周期,提升產(chǎn)品質(zhì)量。CAD軟件用于設(shè)計建模,CAE軟件用于對設(shè)計進(jìn)行驗(yàn)證、優(yōu)化、性能分析等。然而機(jī)械設(shè)備的設(shè)計越來越復(fù)雜,CAE軟件又以有限元求解計算最為通用,模型越復(fù)雜、網(wǎng)格劃分越細(xì)(求解結(jié)果越精確),對計算機(jī)性能的要求就越高,普通PC機(jī)應(yīng)對大型CAE軟件的大規(guī)模并行計算需求顯得力不從心,這就需要用超級計算機(jī)來支持大規(guī)模并行計算。

    超級計算機(jī)擁有大量的計算資源,比如上萬個核的CPU(中央處理單元)、PB(Peta Byte,拍字節(jié))級的內(nèi)存、高速的局域網(wǎng)絡(luò)??赏ㄟ^大型CAE軟件的并行模塊與超級計算機(jī)的調(diào)度軟件配合來對有限元模型進(jìn)行大規(guī)模并行計算求解,合并計算結(jié)果后返回給用戶。

    大型CAE軟件十分昂貴,在超級計算機(jī)上適配各種大型CAE軟件后,在超級計算機(jī)上搭建服務(wù)于工業(yè)設(shè)計仿真需求的云計算服務(wù)平臺(以下簡稱“超算集成平臺”)成為繼超級計算機(jī)硬件研發(fā)的后續(xù)軟件研發(fā)的熱點(diǎn)[1]。本文的研究焦點(diǎn)就是超級計算機(jī)(“天河一號”)環(huán)境下的工業(yè)設(shè)計仿真云平臺中的中間件軟件。

    1 國內(nèi)外研究現(xiàn)狀與趨勢

    綜合國內(nèi)外情況來看,工業(yè)設(shè)計仿真云平臺中間件軟件的研究現(xiàn)狀與趨勢如下:

    (1)Web服務(wù)和代理技術(shù)已滲透進(jìn)入云計算中間件領(lǐng)域。王海偉等[2]提出了基于面向服務(wù)的體系架構(gòu)(service-oriented architecture,SOA)的協(xié)同仿真環(huán)境框架,并最終建了一個簡單的原型CAE平臺,但只是構(gòu)建了一個簡單的框架及原型系統(tǒng)。于加晴等[3]提出了一種在分布式環(huán)境下使用Web方式進(jìn)行協(xié)同設(shè)計的方法,對CAE集成的協(xié)同設(shè)計系統(tǒng)的體系結(jié)構(gòu)和實(shí)現(xiàn)技術(shù)進(jìn)行了介紹,服務(wù)資源通過Jini TM和SOA框架來實(shí)現(xiàn)。

    (2)工業(yè)設(shè)計仿真云平臺都需要利用中間件軟件的支撐使用戶可以通過Web方式來使用HPC(High Performance Computing,高性能計算)的計算資源。Rius等[4]提出在超級計算機(jī)中構(gòu)建一個名為“T-Coffee”的門戶系統(tǒng)來進(jìn)行基因組學(xué)的研究,該系統(tǒng)提供了低成本和高效的專門工具和接口來開發(fā)用戶的交互功能,并支持對遠(yuǎn)程計算資源的訪問。高潭等[5]提出在“銀河”超級計算機(jī)上構(gòu)建一個高性能計算門戶系統(tǒng),系統(tǒng)中的中間件應(yīng)用訪問資源的適配器技術(shù),直接采用網(wǎng)格中間件GT4(Globus Toolkit 4),利用GT4的內(nèi)置適配器對PBS、LSF、SLURM等任務(wù)調(diào)度器進(jìn)行適配。

    (3)已搭建的工業(yè)設(shè)計仿真云平臺結(jié)構(gòu)復(fù)雜,均以層次結(jié)構(gòu)分塊實(shí)現(xiàn)。Kim等[6]提出了一個具有大自由度的線性結(jié)構(gòu)分析系統(tǒng),系統(tǒng)具有三層結(jié)構(gòu),并以“tachyon”超級計算機(jī)作為有限元計算的并行求解平臺。遲學(xué)斌等[7]描述實(shí)現(xiàn)了一種網(wǎng)格中間件SCE(super computing environment,超級計算環(huán)境),該軟件的體系架構(gòu)分為前端服務(wù)、系統(tǒng)服務(wù)和HPC資源三層。張仙偉等[8]根據(jù)網(wǎng)格高性能計算環(huán)境下并行計算的需求,設(shè)計了一個可運(yùn)行在網(wǎng)格上的GPCS(grid HPC-based parallel computing system,高性能計算平臺上的并行計算系統(tǒng))。

    (4)國內(nèi)正努力基于超級計算機(jī)自主研發(fā)工業(yè)設(shè)計仿真云平臺中間件。張洪武等[9]已經(jīng)研發(fā)自主的面向CAE工程與科學(xué)計算集成化軟件平臺SiPESC。國家超算濟(jì)南中心基于“神威藍(lán)光”超級計算機(jī)正在研發(fā)基于超級計算的工業(yè)設(shè)計資源中心技術(shù)開發(fā)與集成技術(shù)平臺,其中就包括該平臺的中間件軟件。

    為充分利用超級計算機(jī)的計算資源、降低用戶使用的難度、促進(jìn)在各種行業(yè)的應(yīng)用,國內(nèi)外都在研究構(gòu)建基于Web的面向CAE或其他行業(yè)應(yīng)用的門戶系統(tǒng),并且都使用中間件軟件來屏蔽超級計算機(jī)任務(wù)調(diào)度器、操作系統(tǒng)的差異,采用的中間件技術(shù)各有差異,有的是采用開源軟件組合,有的是采用網(wǎng)格中間件,有的提出了模型并進(jìn)行了少量的試驗(yàn)。但普遍由于CAE軟件品種繁多、接口復(fù)雜,集成的CAE軟件品種不多,并且在國內(nèi)還未見基于“天河一號”超級計算機(jī)的工業(yè)設(shè)計仿真云平臺中間件出現(xiàn)。

    現(xiàn)階段,需要先行研發(fā)出具有自主知識產(chǎn)權(quán)、能適配“天河一號”甚至更多超級計算機(jī)的工業(yè)設(shè)計仿真云平臺中間件軟件,再考慮如何以網(wǎng)格形式進(jìn)行資源整合。

    2 中間件架構(gòu)設(shè)計

    2.1仿真云平臺總體設(shè)計

    超算集成平臺的總體架構(gòu)設(shè)計如圖1所示,該平臺分為前端系統(tǒng)、中間件系統(tǒng)、集群端系統(tǒng)三個部分。前端系統(tǒng)負(fù)責(zé)與用戶交互。中間件系統(tǒng)進(jìn)行CAE軟件集成,屏蔽異種集群、CAE軟件的差異,為云計算環(huán)境下的“超算集成平臺”前端系統(tǒng)提供統(tǒng)一的接口,支持各種工業(yè)設(shè)計仿真的大規(guī)模并行計算需求,進(jìn)行作業(yè)的提交和監(jiān)控,并將計算結(jié)果返回給前端系統(tǒng)顯示,支持對計算資源(包括CAE軟件License)的管控,支持協(xié)同仿真的業(yè)務(wù)流程組合。

    圖1 平臺總體架構(gòu)設(shè)計

    2.2中間件系統(tǒng)設(shè)計思想

    中間件系統(tǒng)采用SOA架構(gòu)設(shè)計,將各種功能組件封裝為Web Service供前端系統(tǒng)使用,形成系統(tǒng)集成與通信的接口;采用XML作為通信報文,使用XML Schema進(jìn)行數(shù)據(jù)校驗(yàn);在報文內(nèi)部標(biāo)簽中約定某個標(biāo)簽來描述目標(biāo)功能模塊及安全認(rèn)證要求,從而進(jìn)行路由配置,在代碼中實(shí)現(xiàn)報文加解密及路由分發(fā)功能;使用Spring作為Web Service組件對應(yīng)用的JavaBean管理容器,借助AOP(aspect-oriented programming,面向方面的程序設(shè)計)編程思想,形成簡單易配的業(yè)務(wù)流程處理機(jī)制。

    我們面對的計算機(jī)集群主要是Windows集群和Linux集群(“天河一號”超級計算機(jī))2種,針對不同的目標(biāo)集群需要由中間件系統(tǒng)來生成不同的作業(yè)與資源管理指令批處理腳本,再由中間件向目標(biāo)集群提交執(zhí)行,在執(zhí)行過程中監(jiān)控并取得中間和最終計算結(jié)果。為及時調(diào)度干預(yù)計算作業(yè),可通過中間件獲取集群中各種資源的狀態(tài)及使用情況,資源主要包括計算節(jié)點(diǎn)、CPU/GPU、CAE軟件License等。

    2.3中間件體系架構(gòu)

    超算集成平臺的中間件就像一條ESB(enterprise service bus,企業(yè)服務(wù)總線)總線,其他的系統(tǒng)以可插拔的方式掛接在上面,通過Web Service進(jìn)行交互式處理,有同步和異步兩種數(shù)據(jù)交換方式。在ESB總線上,以開源的Tuscany作為SOA架構(gòu)基礎(chǔ)支撐軟件。Tuscany的SCA(service component architecture,服務(wù)組件架構(gòu))容器與Spring的IoC容器相互配合,進(jìn)行構(gòu)件的裝配,如圖2所示。

    圖2 中間件容器配合

    SCA容器可進(jìn)行組合構(gòu)件(composite)和構(gòu)件(component)的裝配,但構(gòu)件的最終實(shí)現(xiàn)由Spring Bean來實(shí)現(xiàn),因此Spring是實(shí)際上的構(gòu)件工廠,并以單例的方式提供多線程的構(gòu)件服務(wù),從而擴(kuò)展Tuscany的構(gòu)件裝配能力和構(gòu)件功能范圍。

    ESB總線的架構(gòu)如圖3所示。中間件通過Spring集成的Quartz作為任務(wù)調(diào)度器,以很小的時間片(不大于5 s)觸發(fā)啟動一個進(jìn)程的線程,分別有作業(yè)監(jiān)控、作業(yè)提交、資源監(jiān)控、License管理4個進(jìn)程,Spring作為組件的容器成熟穩(wěn)定,使進(jìn)程之間相對獨(dú)立,線程之間互不影響,以提高中間件系統(tǒng)的并發(fā)處理能力。

    圖3 ESB總線的總體架構(gòu)

    Spring可使用@AspectJ、Schema等方式來配置AOP的情況,在配置中明確目標(biāo)對象是以切面形成的連接點(diǎn),以及針對這些連接點(diǎn)的增強(qiáng),中間件各連接點(diǎn)配置的增強(qiáng)為前置增強(qiáng)和后置增強(qiáng),以作出一些可供共享的通用處理功能,如日志記錄、權(quán)限檢查、報文加密、報文解密等。

    3 關(guān)鍵技術(shù)

    在超級計算云環(huán)境下,研制工業(yè)設(shè)計仿真云平臺中間件的主要目的是讓用戶在互聯(lián)網(wǎng)環(huán)境下可以簡便地使用海量的計算資源來進(jìn)行工業(yè)設(shè)計工程應(yīng)用,而需要大規(guī)模并行計算的情形主要就是針對CAE軟件的求解器計算應(yīng)用,在“天河一號”超級計算機(jī)上是以一個作業(yè)來對應(yīng)一個求解器計算任務(wù)的,因此首先需要考慮集群環(huán)境下的作業(yè)狀態(tài)是如何轉(zhuǎn)換的,在引入本文所研制出的中間件后又是如何轉(zhuǎn)換的;其次就是要考慮如何來提交作業(yè),以及如何監(jiān)控作業(yè)執(zhí)行的情況并指出在什么情況下作出什么樣的處理。下面以工業(yè)設(shè)計仿真云平臺中間件軟件中的作業(yè)狀態(tài)轉(zhuǎn)換功能為例進(jìn)行分析。

    3.1SLURM作業(yè)狀態(tài)控制

    作業(yè)是超級計算機(jī)進(jìn)行調(diào)度的基本單位,常用的調(diào)度器軟件主要有LSF、SLURM等。“天河一號”超級計算機(jī)主要使用SLURM。作業(yè)在提交給SLURM后,將經(jīng)歷一系列的狀態(tài),如圖4所示。

    圖4 作業(yè)狀態(tài)轉(zhuǎn)換圖

    作業(yè)提交給SLURM后,即進(jìn)入Pending狀態(tài),Pending表示作業(yè)在隊(duì)列中排隊(duì),等待分配資源;作業(yè)被SLURM分配資源后進(jìn)入Running狀態(tài);用戶將Running狀態(tài)下的作業(yè)掛起,則釋放資源,讓出CPU,進(jìn)入Suspended狀態(tài),但內(nèi)存并不釋放;作業(yè)運(yùn)行結(jié)束轉(zhuǎn)為Completed狀態(tài);作業(yè)運(yùn)行失敗轉(zhuǎn)入Failed狀態(tài);作業(yè)在Pending、Suspended或Running狀態(tài)時被取消則進(jìn)入Canncelled狀態(tài);如果作業(yè)在運(yùn)行超出了請求的運(yùn)行時間則進(jìn)入Timeout狀態(tài);如果分配給作業(yè)的節(jié)點(diǎn)出現(xiàn)故障,且作業(yè)提交時沒有指定節(jié)點(diǎn)故障容忍參數(shù),則轉(zhuǎn)入Nodefail狀態(tài)。

    3.2中間件作業(yè)狀態(tài)控制

    “天河一號”超級計算機(jī)通過SLURM調(diào)度器調(diào)度求解作業(yè),提交作業(yè)之后的狀態(tài)由SLURM調(diào)度器控制,在此基礎(chǔ)之上需要對狀態(tài)進(jìn)行重新封裝和控制,要考慮解決以下問題:

    (1)用戶以Web方式提交工業(yè)設(shè)計仿真計算作業(yè),中間件要根據(jù)用戶Web方式提交的作業(yè)需求來自動向“天河一號”超級計算機(jī)提交作業(yè)。

    (2)作業(yè)提交后作業(yè)狀態(tài)交由SLURM調(diào)度器控制,中間件要以較小的時間片(不大于5 s)不間斷地輪詢SLURM以獲取作業(yè)提交后的狀態(tài),并作出相應(yīng)的處理。

    (3)考慮計費(fèi)功能、網(wǎng)絡(luò)與軟件故障處理,需要適當(dāng)增加作業(yè)的控制狀態(tài)。

    (4)由于SLURM調(diào)度器內(nèi)的作業(yè)狀態(tài)有可能跳轉(zhuǎn)很快而在中間件的時間片中監(jiān)控不到,在中間件的作業(yè)狀態(tài)上允許非人為干預(yù)下的跨狀態(tài)跳轉(zhuǎn)。

    根據(jù)以上情況,中間件軟件設(shè)計的作業(yè)狀態(tài)轉(zhuǎn)換情況如圖5所示。

    圖5 中間件軟件設(shè)計的作業(yè)狀態(tài)轉(zhuǎn)換情況

    中間件并不修改SLURM的源代碼來控制作業(yè)的狀態(tài),而是集成SLURM后在SLURM原有狀態(tài)控制功能的基礎(chǔ)上進(jìn)行狀態(tài)轉(zhuǎn)換的改進(jìn)。根據(jù)SLURM調(diào)度器中的作業(yè)狀態(tài)模型和中間件調(diào)整后的作業(yè)狀態(tài)模型的對比情況,對狀態(tài)轉(zhuǎn)換作了以下調(diào)整:

    (1)為讓用戶通過前端系統(tǒng)可將新建作業(yè)刪除而增加了Deleted(被刪除)狀態(tài),但作業(yè)一旦通過前端系統(tǒng)提交則不允許再行刪除。

    (2)增加了Waitcommitted(等待提交)狀態(tài)。在前端系統(tǒng)中提交作業(yè)后就進(jìn)入Waitcommitted狀態(tài);如果通過中間件的作業(yè)提交線程提交作業(yè)成功,則轉(zhuǎn)為Pending狀態(tài);如果作業(yè)提交失敗,則重試,重試3次仍提交不成功置為Failed狀態(tài),在失敗原因中注明原因是作業(yè)提交失敗。

    (3)作業(yè)提交后交由SLURM管控作業(yè)狀態(tài),但中間件的監(jiān)控進(jìn)程通過Linux命令在每個時間片的線程中得到作業(yè)的狀態(tài)并更新,因此會引入非人為干預(yù)情況下的跨狀態(tài)跳轉(zhuǎn)情況,增加了從Pending跳轉(zhuǎn)到Failed、Completed狀態(tài)的轉(zhuǎn)換關(guān)系,以及從Suspended跳轉(zhuǎn)到Failed、Completed狀態(tài)的轉(zhuǎn)換關(guān)系。

    (4)取消了原來的Timeout和Nodefail狀態(tài),合并到Failed狀態(tài),但在Failed狀態(tài)的失敗原因中注明是什么原因。

    3.3作業(yè)提交功能設(shè)計

    作業(yè)提交功能的設(shè)計思想如圖6所示,提交的過程如下:

    圖6 作業(yè)提交功能設(shè)計思想

    (1)以Spring集成的Quartz作為時鐘觸發(fā)器,時間片到時觸發(fā)作業(yè)監(jiān)控提線程查詢作業(yè)隊(duì)列中狀態(tài)為“Waitcommitted”的作業(yè)。

    (2)根據(jù)作業(yè)所需要的求解器類型與版本、模型文件和指令文件的情況、所需計算資源的情況等來綜合生成腳本文件。

    (3)調(diào)用腳本文件生成器的接口來通過SFTP程序?qū)⒛P臀募椭噶钗募?、腳本文件傳到“天河一號”超級計算機(jī)上。

    (4)、(5)生成執(zhí)行腳本文件的指令通過指令執(zhí)行通道執(zhí)行提交批處理作業(yè)指令。如果提交成功得到作業(yè)的ID號,將作業(yè)狀態(tài)變更為“Pending”;如果提交失敗,將作業(yè)狀態(tài)變更為“Failed”。

    以上過程中的狀態(tài)轉(zhuǎn)換工作均由作業(yè)狀態(tài)轉(zhuǎn)換器來完成,如果發(fā)生異常則置作業(yè)狀態(tài)為“F”,并記錄失敗原因。

    以作業(yè)提交進(jìn)程時間片到時線程執(zhí)行的程序設(shè)計思想為例,算法思想描述如下:

    //得到一個新的作業(yè)

    NewJob newJob=SelectANewJob();

    //如果作業(yè)的執(zhí)行集群為“天河一號”

    if(newJob.object equal‘天河一號’){

    //生成作業(yè)提交腳本文件

    GenerateCommitScriptFile(newJob);

    //置作業(yè)上傳次數(shù)為0

    newJob.uploadCount=0;

    //用戶作業(yè)預(yù)扣費(fèi)

    boolean moneySuccess=PreMoney(newJob,user);

    //如果預(yù)扣費(fèi)成功

    if(moneySuccess){

    //上傳并提交作業(yè)

    UploadAndCommit(newJob);

    }else{

    //修改作業(yè)狀態(tài)為“Failed” ,并設(shè)置失敗原因

    ModeJobStatus(newJob,'Failed','用戶賬上余額不足');

    }

    }

    //上傳并提交作業(yè)功能

    void UploadAndCommit(newJob){

    //記錄作業(yè)上傳次數(shù)

    newJob.uploadCount++;

    //如果已經(jīng)上傳3次

    if(newJob.uploadCount>3){

    //作業(yè)狀態(tài)為'Failed',并設(shè)置失敗原因

    ModeJobStatus(newJob,'Failed','作業(yè)文件上傳失敗');

    //返回用戶作業(yè)預(yù)扣費(fèi)

    ReturnMoney(newJob,user);

    return;

    }

    //通過文件傳輸通道上傳作業(yè)指令文件、模型文件及其它文件,返回是否成功標(biāo)志

    booleanuploadJobFilesFlag=UploadJobFiles(newJob);

    //通過文件傳輸通道上傳作業(yè)提交腳本文件

    Boolean uploadScriptFileFlag

    =UploadScriptFile(newJob);

    //如果上傳成功

    if(uploadJobFilesFlag and uploadScriptFileFlag){

    //通過以批處理方式提交作業(yè),返回是否成功

    boolean commitFlag=CommitJob(newJob);

    //如果提交成功

    if(commitFlag){

    //修改作業(yè)狀態(tài)為'Pending'

    ModiJobStatus(newJob,'Pending');

    }else{

    //修改作業(yè)狀態(tài)為“Failed”,并設(shè)置原因

    ModiJobStatus(newJob,'Failed','作業(yè)提交失敗');

    //返回用戶作業(yè)預(yù)扣費(fèi)

    ReturnMoney(newJob,user);

    }

    }else{

    //利用遞歸方法再次上傳和提交作業(yè)

    UploadAndCommit(newJob);

    }

    }

    4 工業(yè)設(shè)計仿真實(shí)例

    下面以火星著陸多室連通氣囊ANSYS/LS-DYNA有限元分析計算的前處理、求解計算、后處理[10]為例進(jìn)一步說明中間件的應(yīng)用。

    4.1火星著陸多室連通氣囊?guī)缀文P图坝邢拊P偷慕?/p>

    火星著陸多室連通氣囊緩沖裝置示意圖見圖7a。在該氣囊緩沖裝置中,總共有A、B、C、D 4個完全對稱的氣囊子系統(tǒng),四面體探測器的每個面通過法蘭與4個氣囊子系統(tǒng)相連,這樣探測器就被包裹在這4個氣囊子系統(tǒng)中。氣囊子系統(tǒng)A示意圖見圖7b,由A1~A6總共6個直徑相同的球囊組合而成,其中3個角上的球囊A1、A3和A5被稱為角囊,中間位置球囊A2、A4、A6被稱為中間囊,角囊和中間囊之間的隔膜上有通氣孔進(jìn)行連接,可以使得氣體在氣囊子系統(tǒng)內(nèi)部自由流通。在該多室連通氣囊系統(tǒng)中,氣囊子系統(tǒng)C的中間囊上有3個通氣孔分別與相鄰的氣囊子系統(tǒng)A、B、D的中間囊相連。在探測器著陸沖擊過程中,氣囊子系統(tǒng)C在減速傘相對的下端,所以其最先受到?jīng)_擊,其內(nèi)部氣體受到探測器壓力作用經(jīng)過通氣孔被擠壓到氣囊子系統(tǒng)A、B、D中[10]。

    (a)整體結(jié)構(gòu)示意圖 (b)氣囊子系統(tǒng)A圖7 火星著陸連通多室氣囊緩沖裝置

    根據(jù)上述幾何結(jié)構(gòu)及圖8所示結(jié)構(gòu)尺寸,建立了火星著陸多室連通氣囊緩沖系統(tǒng)的有限元模型,如圖9所示(Dbag= 1058 mm)。其中,采用LS-DYNA的*AIRBAG_INTERACTION定義氣囊之間的連通,氣囊子系統(tǒng)與火星探測器之間的連接法蘭一側(cè)與氣囊節(jié)點(diǎn)重合,另一側(cè)與剛性探測器節(jié)點(diǎn)從屬。每個球囊均采用控制體積法來計算其內(nèi)部氣體的熱力學(xué)參數(shù)。球囊的單元為Belytschko-Tsay膜單元,單元在厚度方向的積分點(diǎn)個數(shù)為3,在面內(nèi)采用單點(diǎn)積分?;鹦翘綔y器和地面的單元均為Belytschko-Tsay殼單元,單元在厚度方向的積分點(diǎn)個數(shù)為2,在面內(nèi)采用單點(diǎn)積分。球囊的材料本構(gòu)模型為LS-DYNA的MAT_FABRIC,火星探測器和剛性地面的材料本構(gòu)模型均采用LS-DYNA的MAT_RIGID。在接觸計算中,采用自動面面接觸算法來計算球囊與火星探測器以及球囊與剛性地面之間的接觸。

    圖8 火星著陸多室連通氣囊緩沖裝置尺寸圖

    圖9 火星著陸多室連通氣囊緩沖系統(tǒng)有限元模型(Dbag =1058 mm)

    在緩沖著陸過程中,探測器與多室氣囊一起以20 m/s的初速度沖擊到火星表面上,著陸初速度與豎直方向的夾角為α,在有限元模型中采用LS-DYNA的INITIAL_VELOCITY來實(shí)現(xiàn)。為模擬火星重力環(huán)境,在有限元分析中設(shè)置環(huán)境重力場g=3.76 m/s2,火星的表面大氣壓力為650 Pa。

    4.2求解計算

    用戶使用工業(yè)設(shè)計仿真云平臺進(jìn)行火星著陸氣囊的有限元碰撞分析計算,從提交作業(yè)到取回計算結(jié)果的完整流程如下。

    (1)用戶通過超算集成平臺門戶新建一個LS-DYNA作業(yè),設(shè)置好相應(yīng)的參數(shù)和輸入文件,如請求的CPU數(shù)量、內(nèi)存數(shù)量等,保存并提交作業(yè)。

    (2)中間件的作業(yè)提交線程獲取到狀態(tài)為“Waitcommitted”狀態(tài)的作業(yè)列表,其中就包括名為“火星著陸氣囊”的作業(yè),生成作業(yè)提交腳本文件shield.sh,這個文件實(shí)際上是一個腳本文件,其核心的算法思想如下:

    #得出需要的結(jié)點(diǎn)數(shù)量

    int nodenum=0;

    #“天河一號”每個結(jié)點(diǎn)2個CPU,每個CPU擁有8個核

    # processnum為用戶請求的并行數(shù)

    if(processnum%16==0)

    nodenum= processnum/16;

    else

    nodenum= processnum/16+1;

    #得出空間的結(jié)點(diǎn)列表

    Node idlenodelist[]=new Node[nodenum];

    int j=0;

    for(int i=0;i

    if(j== nodenum) break;

    if(cn[i] is idle) idlenodelist[j]=cn[i];

    }

    If(j< nodenum) return 0;#資源數(shù)量不足

    #執(zhí)行并行計算

    lsdyna145 -dis -machinesidlenodelist[0]:2…

    idlenodelist[idlenodelist.length-1]:2

    i=kfilename pr=dyna

    接下來通過文件傳輸通道把腳本文件、待求解的火星著陸氣囊模型文件上傳到“天河一號”的計算指定目錄下,再通過指令通道執(zhí)行批處理指令,如下所示:

    sbatch shield.sh

    執(zhí)行后如果成功則將作業(yè)交由SLURM進(jìn)行調(diào)度,待所有資源需求滿足后,SLURM調(diào)度器會調(diào)度相應(yīng)LS-DYNA軟件進(jìn)行加載求解。

    (3)在作業(yè)運(yùn)行過程中,用戶可以通過門戶系統(tǒng)的作業(yè)監(jiān)控模塊隨時了解作業(yè)當(dāng)時的運(yùn)行狀態(tài),監(jiān)控的狀態(tài)來源于中間件的作業(yè)監(jiān)控線程采集的作業(yè)數(shù)據(jù)。在作業(yè)狀態(tài)發(fā)生改變時,用戶可通過門戶查看和下載本次求解計算生成的中間或結(jié)果文件,這些文件由中間件的作業(yè)監(jiān)控進(jìn)程下載得到。

    (4)下載得到火星著陸氣囊的有限元碰撞分析結(jié)果后,將下載得到的文件導(dǎo)入到Windows下的LS-PREPOST軟件中即可得到所需的計算結(jié)果。

    5 下一步研究工作

    (1)集成更多的主流CAE軟件,如Fluent、ADINA、Nastran等;集成其他調(diào)度器軟件,如LSF。

    (2)解決在線后處理問題。雖然現(xiàn)有平臺可以Web方式在線下載計算結(jié)果文件,但由于計算結(jié)果文件尺寸較大,常達(dá)GB數(shù)量級,用戶下載時間較長,可考慮采用VNC(virtual network computer,虛擬網(wǎng)絡(luò)計算機(jī))方式使用戶可直接以Web方式在線查看和分析計算結(jié)果,從而免去下載計算結(jié)果的步驟。

    (3)基于工作流模型進(jìn)行多學(xué)科聯(lián)合仿真及其在工業(yè)設(shè)計仿真云平臺上的實(shí)現(xiàn),可先從單學(xué)科CAD、CAE軟件的集成,單CAE軟件多學(xué)科聯(lián)合仿真與優(yōu)化分析開始,再進(jìn)一步研究多學(xué)科多軟件集成聯(lián)合仿真。

    (4)多超算中心工業(yè)設(shè)計仿真云平臺的網(wǎng)格化處理及中間件的網(wǎng)格化研究?,F(xiàn)有研發(fā)的平臺體現(xiàn)了用戶以簡單的方式使用云資源的優(yōu)點(diǎn),但尚不能利用網(wǎng)格資源,可考慮多超算中心聯(lián)合研發(fā)、構(gòu)建計算資源網(wǎng)格中間件,從而通過平臺可以使用異種超級計算機(jī)的計算資源。

    [1]Perumalla K S.Tutorial:Parallel Simulation on Supercomputers[C]//2012 Winter Simulation Conference,WSC 2012.Berlin,2012:00919620785775.

    [2]Wang Haiwei,Liu Geng,Han Bing,et al.Collaborative Simulation Environment Framework Based on SOA[C]//2008 12th International Conference on Computer Supported Cooperative Work in Design,CSCWD 2008.Xi’an,2008:416-419.

    [3]Yu Jiaqing,Cha Jianzhong,Lu Yiping,et al.A CAE-integrated Distributed Collaborative Design System for Finite Element Analysis of Complex Product Based on SOOA[J].Advances in Engineering Software,2010,41(4):590-603.

    [4]Rius J,Cores F,Solsona F,et al.A User-friendly Web Portal for T-Coffee on Supercomputers[J].BMCBioinformatics England,2011,12(5):150-154.

    [5]高潭.高性能計算門戶的研究與實(shí)現(xiàn)[D].長沙:國防科學(xué)技術(shù)大學(xué),2009.

    [6]Kim J,Lee S,Jung H,et al.Development of Large-scale Structural Analysis System on a Supercomputer[C]//2011 International Conference on High Performance Computing and Simulation,HPCS 2011.United States,2011:825-829.

    [7]遲學(xué)斌,肖海力,王小寧,等.面向科學(xué)計算的網(wǎng)格環(huán)境[J].集成技術(shù),2012,1(1):68-76.

    Chi Xuebin,Xiao Haili,Wang Xiaoning,et al.Scientific Computing Grid and SCE Middleware[J].Journal of Integration Technology,2012,1(1):68-76.

    [8]張仙偉 張璟.基于網(wǎng)格計算平臺的并行計算系統(tǒng)研究與實(shí)現(xiàn)[J].計算機(jī)工程與應(yīng)用,2012,48(7):5-8,11.

    Zhang Xianwei,Zhang Jing.Research and Implementation of Grid Computing-based Parallel Computing System[J].Computer Engineering and Applications,2012,48(7):5-8.

    [9]張洪武,陳飆松,李云鵬,等.面向集成化CAE軟件開發(fā)的SiPESC 研發(fā)工作進(jìn)展[J].計算機(jī)輔助工程,2011,20(2):39-49.

    Zhang Hongwu,Chen Biaosong,Li Yunpeng,et al.Advancement of Design and Implementation of SiPESC for Development of Integrated CAE Software Systems[J].Computer Aided Engineering,2011,20(2):39-49.[10]文桂林,雷志華,尹建武,等.全向式多室連通氣囊的緩沖特性研究[J].振動與沖擊,2013,32(8):13-17.

    Wen Guilin,Lei Zhihua,Yin Jianwu,et al.Cushion Characteristics of An Omni-directional and Multi-chamber Airbag[J].Journal of Vibration and Shock,2013,32(8):13-17.

    (編輯袁興玲)

    Research on Middleware of Cloud Platform for Industrial Design and Simulation on “Tianhe No.1” Super Computer

    Deng Ziyun1,2Zhang Jing1,2Liu Yangbing1,2Xiao Jiuru1

    1.Hunan University,Changsha,410082 2.National Supercomputing Center in Changsha,Changsha,410082

    A cloud middleware was developed as a core component of the cloud platform for industrial design and simulation which was built on the “Tianhe No.1” super computer.The middleware integrated a large number of CAE software,including SLURM job scheduling software,License management software,and computing resource management function,etc.The whole framework of the cloud platform for industrial design and simulation,and the architecture of the middleware software were designed.Key technical problems including operational state transition,job submission etc,were investigated and solved,the cloud platform and the middleware software were developed for applications on the statics analysis of the shield machine.An experiment for the collision analysis of the Mars landing airbag was carried out on the cloud platform by using the middleware.The results show that,it is convenient to use “Tianhe No.1” super computing resources by way of the middleware in the cloud platform for industrial design and simulation.

    “Tianhe No.1” super computer;middleware in cloud platform for industrial design and simulation; architecture;design thought;key technology;simulation experiment

    2012-02-04

    國家科技支撐計劃資助項(xiàng)目(2012BAH09B02);長沙市重大科技專項(xiàng)資助項(xiàng)目(K1204006-11-1)

    TP311.52< class="emphasis_italic">DOI

    :10.3969/j.issn.1004-132X.2015.06.011

    鄧子云,男,1979年生。湖南大學(xué)電氣與信息工程學(xué)院博士研究生、教授、高級工程師。主要研究方向?yàn)镠PC、物流信息技術(shù)。章兢,男,1959年生。湖南大學(xué)電氣與信息工程學(xué)院教授、博士研究生導(dǎo)師。劉楊兵,男,1985年生。湖南大學(xué)信息科學(xué)與工程學(xué)院碩士研究生。肖久如,男,1989年生。湖南大學(xué)機(jī)械與運(yùn)載工程學(xué)院碩士研究生。

    猜你喜歡
    天河中間件超級計算機(jī)
    超級計算機(jī)
    湖南省瀟水涔天河水庫擴(kuò)建工程通過竣工驗(yàn)收
    中國水利(2022年1期)2022-02-13 07:37:00
    超級計算機(jī)及其在航空航天領(lǐng)域中的應(yīng)用
    科技傳播(2019年22期)2020-01-14 03:06:36
    一條天河走運(yùn)來
    北方音樂(2019年13期)2019-08-21 02:14:32
    RFID中間件技術(shù)及其應(yīng)用研究
    電子制作(2018年14期)2018-08-21 01:38:10
    美國制造出全球最快超級計算機(jī)
    基于VanConnect中間件的設(shè)計與開發(fā)
    電子測試(2018年10期)2018-06-26 05:54:02
    每秒100億億次 中國超級計算機(jī)
    天河CBD:集聚創(chuàng)新,遇見城市未來
    空中之家(2017年11期)2017-11-28 05:27:45
    “天河二號”獲全球超算五連冠等
    日本爱情动作片www.在线观看| 欧美精品人与动牲交sv欧美| 日韩一本色道免费dvd| 777米奇影视久久| 日韩精品有码人妻一区| 蜜桃久久精品国产亚洲av| 草草在线视频免费看| 国产色婷婷99| 色哟哟·www| 欧美97在线视频| 国产成人精品婷婷| 免费看光身美女| a 毛片基地| www.色视频.com| 色5月婷婷丁香| 狂野欧美白嫩少妇大欣赏| 啦啦啦视频在线资源免费观看| 国产亚洲91精品色在线| 丰满人妻一区二区三区视频av| av在线观看视频网站免费| 在线亚洲精品国产二区图片欧美 | 超碰97精品在线观看| 亚洲精品国产色婷婷电影| 丝袜在线中文字幕| 多毛熟女@视频| 国产精品不卡视频一区二区| 男女国产视频网站| 五月开心婷婷网| 少妇熟女欧美另类| 日本av免费视频播放| 99久久精品一区二区三区| 热99国产精品久久久久久7| 日韩中字成人| 亚洲国产日韩一区二区| 国产在视频线精品| 99久久中文字幕三级久久日本| 国产片特级美女逼逼视频| 高清欧美精品videossex| 男女边摸边吃奶| 涩涩av久久男人的天堂| 久久韩国三级中文字幕| 免费高清在线观看视频在线观看| 精品一品国产午夜福利视频| 丝袜脚勾引网站| 国产精品人妻久久久影院| 女性生殖器流出的白浆| 亚洲国产欧美日韩在线播放 | a级一级毛片免费在线观看| 99九九在线精品视频 | 99热国产这里只有精品6| a级片在线免费高清观看视频| 亚洲欧美中文字幕日韩二区| 美女国产视频在线观看| 伊人久久国产一区二区| 久久99热这里只频精品6学生| av在线观看视频网站免费| 亚洲不卡免费看| 精品久久久久久久久亚洲| 国产成人精品福利久久| 亚洲精品中文字幕在线视频 | 亚洲精品一区蜜桃| 国产成人精品婷婷| 色5月婷婷丁香| av线在线观看网站| 2021少妇久久久久久久久久久| 夫妻午夜视频| 九色成人免费人妻av| 日韩不卡一区二区三区视频在线| 最近中文字幕2019免费版| 伦精品一区二区三区| 亚洲国产精品国产精品| 看十八女毛片水多多多| 超碰97精品在线观看| 校园人妻丝袜中文字幕| 如日韩欧美国产精品一区二区三区 | 亚洲成色77777| 99久国产av精品国产电影| 又粗又硬又长又爽又黄的视频| 亚洲,欧美,日韩| 亚洲欧美日韩另类电影网站| 久久久精品94久久精品| 日日啪夜夜撸| 久久精品国产亚洲av天美| 最近中文字幕高清免费大全6| 亚洲天堂av无毛| 国产精品人妻久久久久久| 日韩免费高清中文字幕av| 日日撸夜夜添| 免费av不卡在线播放| 女的被弄到高潮叫床怎么办| 草草在线视频免费看| 欧美人与善性xxx| 少妇猛男粗大的猛烈进出视频| 亚洲精品久久久久久婷婷小说| 少妇精品久久久久久久| av国产久精品久网站免费入址| 国产极品粉嫩免费观看在线 | 精品熟女少妇av免费看| 2018国产大陆天天弄谢| 亚洲av中文av极速乱| 国产在线一区二区三区精| 亚洲精品aⅴ在线观看| 亚洲国产精品成人久久小说| 色哟哟·www| 我的老师免费观看完整版| 久久亚洲国产成人精品v| 高清欧美精品videossex| 日日撸夜夜添| a级片在线免费高清观看视频| 国产精品一区www在线观看| 国产男人的电影天堂91| 男女啪啪激烈高潮av片| 国产成人freesex在线| 韩国高清视频一区二区三区| 精品国产乱码久久久久久小说| 精品视频人人做人人爽| 国产精品久久久久成人av| 99热网站在线观看| 免费av不卡在线播放| videos熟女内射| 免费观看的影片在线观看| 色婷婷av一区二区三区视频| 少妇高潮的动态图| av黄色大香蕉| 99久国产av精品国产电影| 免费av中文字幕在线| 如日韩欧美国产精品一区二区三区 | 亚洲,欧美,日韩| 久久久a久久爽久久v久久| 中文字幕免费在线视频6| 秋霞在线观看毛片| 国产精品一二三区在线看| 秋霞在线观看毛片| 青春草国产在线视频| 97精品久久久久久久久久精品| 大香蕉97超碰在线| 国语对白做爰xxxⅹ性视频网站| 亚洲,一卡二卡三卡| 在线观看美女被高潮喷水网站| 精品久久久久久久久av| 多毛熟女@视频| 黑人猛操日本美女一级片| 九色成人免费人妻av| 2018国产大陆天天弄谢| 午夜久久久在线观看| 大片电影免费在线观看免费| 少妇丰满av| 国产亚洲5aaaaa淫片| a级毛片在线看网站| 天堂中文最新版在线下载| 国产精品久久久久久久电影| 人体艺术视频欧美日本| 国产黄片美女视频| 成人亚洲欧美一区二区av| 久久国产精品大桥未久av | av福利片在线| 国产精品成人在线| 亚洲激情五月婷婷啪啪| 国产亚洲5aaaaa淫片| 日本与韩国留学比较| 国产熟女午夜一区二区三区 | 不卡视频在线观看欧美| 十八禁网站网址无遮挡 | 日韩大片免费观看网站| 亚洲精品乱久久久久久| 亚洲一级一片aⅴ在线观看| 精品午夜福利在线看| 中文欧美无线码| 日韩伦理黄色片| 少妇丰满av| 亚州av有码| 欧美另类一区| 精品一区在线观看国产| 9色porny在线观看| 黄色毛片三级朝国网站 | 涩涩av久久男人的天堂| 午夜免费观看性视频| h日本视频在线播放| 全区人妻精品视频| 亚洲国产精品999| 国产午夜精品一二区理论片| 久久久久久久久久人人人人人人| 这个男人来自地球电影免费观看 | 一本大道久久a久久精品| av卡一久久| 国产片特级美女逼逼视频| 国产亚洲一区二区精品| 久久女婷五月综合色啪小说| 亚洲熟女精品中文字幕| 亚洲欧美日韩卡通动漫| 欧美 日韩 精品 国产| 三上悠亚av全集在线观看 | 日韩精品有码人妻一区| 人人妻人人澡人人看| 爱豆传媒免费全集在线观看| 久久 成人 亚洲| 人妻少妇偷人精品九色| 久久精品国产鲁丝片午夜精品| 九草在线视频观看| 天天躁夜夜躁狠狠久久av| 日韩精品有码人妻一区| 亚洲自偷自拍三级| 国产av国产精品国产| 国产美女午夜福利| 观看免费一级毛片| 国产成人freesex在线| a级毛片在线看网站| 亚洲av男天堂| 久久综合国产亚洲精品| 一本一本综合久久| 国产91av在线免费观看| 精品一品国产午夜福利视频| .国产精品久久| 婷婷色麻豆天堂久久| 色婷婷av一区二区三区视频| 9色porny在线观看| 老司机影院成人| 在线亚洲精品国产二区图片欧美 | 久久久久久久大尺度免费视频| 在线观看人妻少妇| kizo精华| 亚洲成人av在线免费| av网站免费在线观看视频| 亚洲欧美成人综合另类久久久| 我要看日韩黄色一级片| 99久国产av精品国产电影| 三级国产精品片| 国产色爽女视频免费观看| 久久久久久久久久成人| 成人毛片60女人毛片免费| 久久精品国产自在天天线| 国产亚洲午夜精品一区二区久久| 国产免费一区二区三区四区乱码| 最近2019中文字幕mv第一页| 久久精品国产亚洲av天美| 在线观看免费视频网站a站| 久久久久国产网址| 久久久久久久久久人人人人人人| 亚洲精品久久久久久婷婷小说| 极品人妻少妇av视频| 国产视频首页在线观看| 亚洲,一卡二卡三卡| 一级毛片我不卡| 国产精品99久久99久久久不卡 | 肉色欧美久久久久久久蜜桃| 日日啪夜夜撸| 最近最新中文字幕免费大全7| 亚洲av电影在线观看一区二区三区| 亚洲成人手机| 欧美日韩av久久| 国产视频内射| 中文字幕人妻熟人妻熟丝袜美| 国产精品国产av在线观看| 久久国产精品大桥未久av | videossex国产| 亚洲精品成人av观看孕妇| 久久久国产欧美日韩av| 亚洲无线观看免费| 高清黄色对白视频在线免费看 | 国产精品一二三区在线看| 在线观看人妻少妇| 国产在线视频一区二区| 欧美最新免费一区二区三区| 在线看a的网站| 午夜av观看不卡| 有码 亚洲区| 日韩欧美精品免费久久| 国产精品久久久久久精品电影小说| 久久久久久久大尺度免费视频| 久久久精品免费免费高清| 免费av不卡在线播放| 国产 一区精品| 一级毛片电影观看| 卡戴珊不雅视频在线播放| 中文乱码字字幕精品一区二区三区| 特大巨黑吊av在线直播| 国产av一区二区精品久久| 夫妻性生交免费视频一级片| 亚洲av中文av极速乱| 少妇丰满av| av在线播放精品| 亚洲av不卡在线观看| 久久狼人影院| 婷婷色麻豆天堂久久| 青春草亚洲视频在线观看| 精品久久久精品久久久| 亚洲国产欧美日韩在线播放 | 亚洲伊人久久精品综合| 黄色怎么调成土黄色| 少妇 在线观看| 国产成人精品一,二区| 三级经典国产精品| 国产黄频视频在线观看| 国产国拍精品亚洲av在线观看| 毛片一级片免费看久久久久| 如何舔出高潮| 欧美一级a爱片免费观看看| 中国美白少妇内射xxxbb| 纵有疾风起免费观看全集完整版| 亚洲怡红院男人天堂| 日本av手机在线免费观看| 秋霞在线观看毛片| 日韩欧美一区视频在线观看 | 观看美女的网站| 国产精品99久久99久久久不卡 | a级毛片在线看网站| 七月丁香在线播放| 国产高清有码在线观看视频| 午夜激情福利司机影院| 韩国高清视频一区二区三区| 久久精品熟女亚洲av麻豆精品| 久久国产精品大桥未久av | 黄色视频在线播放观看不卡| 99热这里只有是精品50| .国产精品久久| 免费观看a级毛片全部| √禁漫天堂资源中文www| 亚洲国产精品成人久久小说| av在线老鸭窝| 国产伦理片在线播放av一区| 日韩一本色道免费dvd| 色婷婷av一区二区三区视频| 国产免费视频播放在线视频| 精品99又大又爽又粗少妇毛片| 亚洲久久久国产精品| 亚洲av中文av极速乱| 大片电影免费在线观看免费| 大香蕉97超碰在线| av黄色大香蕉| 免费av不卡在线播放| 国内揄拍国产精品人妻在线| 丰满少妇做爰视频| 亚洲国产毛片av蜜桃av| 亚洲精品亚洲一区二区| 99九九在线精品视频 | 一边亲一边摸免费视频| 妹子高潮喷水视频| av.在线天堂| 免费大片18禁| 九九爱精品视频在线观看| 大话2 男鬼变身卡| 国产亚洲欧美精品永久| 免费看日本二区| 性高湖久久久久久久久免费观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品第二区| av.在线天堂| 少妇猛男粗大的猛烈进出视频| 伊人久久国产一区二区| 在线观看一区二区三区激情| 国产极品粉嫩免费观看在线 | 亚洲精品国产色婷婷电影| 中文字幕人妻熟人妻熟丝袜美| 2022亚洲国产成人精品| 嫩草影院入口| 午夜免费男女啪啪视频观看| 大话2 男鬼变身卡| 精品午夜福利在线看| 性色avwww在线观看| 国产乱人偷精品视频| 欧美日韩综合久久久久久| 99久久精品热视频| 黄色毛片三级朝国网站 | 午夜福利在线观看免费完整高清在| 国产极品粉嫩免费观看在线 | 欧美精品一区二区免费开放| 亚洲色图综合在线观看| 中文字幕免费在线视频6| 少妇高潮的动态图| 亚洲美女视频黄频| 精品久久久精品久久久| 亚洲美女视频黄频| 日韩av不卡免费在线播放| 欧美日韩在线观看h| 黄色毛片三级朝国网站 | 一区二区三区精品91| 久久久久久久精品精品| 亚洲人成网站在线播| 最近中文字幕2019免费版| 又粗又硬又长又爽又黄的视频| 午夜激情福利司机影院| 高清视频免费观看一区二区| 欧美日韩精品成人综合77777| 男男h啪啪无遮挡| 十八禁高潮呻吟视频 | av在线播放精品| 乱人伦中国视频| 我的老师免费观看完整版| av卡一久久| 亚洲丝袜综合中文字幕| 久久热精品热| 亚洲情色 制服丝袜| 亚洲精品456在线播放app| 波野结衣二区三区在线| 亚洲国产欧美在线一区| 国产中年淑女户外野战色| 久久鲁丝午夜福利片| 亚洲一级一片aⅴ在线观看| 97精品久久久久久久久久精品| 尾随美女入室| 日本vs欧美在线观看视频 | 人体艺术视频欧美日本| 国产亚洲91精品色在线| 婷婷色av中文字幕| av天堂中文字幕网| 亚洲国产欧美在线一区| 亚洲国产精品国产精品| 欧美xxⅹ黑人| 亚洲色图综合在线观看| 又黄又爽又刺激的免费视频.| 国产精品福利在线免费观看| 丰满少妇做爰视频| 蜜桃在线观看..| 看非洲黑人一级黄片| 热99国产精品久久久久久7| 午夜视频国产福利| 国产极品粉嫩免费观看在线 | 国产一区亚洲一区在线观看| 亚洲精品成人av观看孕妇| 色94色欧美一区二区| 久久久精品94久久精品| 欧美3d第一页| 久久人妻熟女aⅴ| 嫩草影院新地址| 亚洲国产毛片av蜜桃av| 日本色播在线视频| 国内精品宾馆在线| 国产伦理片在线播放av一区| 最近最新中文字幕免费大全7| 少妇熟女欧美另类| 波野结衣二区三区在线| 国国产精品蜜臀av免费| 久久久久久久大尺度免费视频| 成人毛片a级毛片在线播放| 啦啦啦中文免费视频观看日本| 国产日韩欧美亚洲二区| 蜜臀久久99精品久久宅男| 精品少妇内射三级| 成年女人在线观看亚洲视频| 久久国产乱子免费精品| 亚洲四区av| 伦精品一区二区三区| 国产成人免费观看mmmm| 国产精品国产三级国产专区5o| 黄色日韩在线| av一本久久久久| 日日啪夜夜爽| 久久热精品热| 一二三四中文在线观看免费高清| 精品视频人人做人人爽| 国产欧美日韩综合在线一区二区 | 亚洲av在线观看美女高潮| 国内少妇人妻偷人精品xxx网站| 国产高清有码在线观看视频| 丝袜喷水一区| 国产午夜精品久久久久久一区二区三区| 王馨瑶露胸无遮挡在线观看| 蜜桃在线观看..| 久久久久久伊人网av| 99热6这里只有精品| 成人亚洲欧美一区二区av| 亚洲一级一片aⅴ在线观看| 丰满饥渴人妻一区二区三| 国产男女超爽视频在线观看| 夫妻午夜视频| 精品99又大又爽又粗少妇毛片| 又大又黄又爽视频免费| 国产成人免费无遮挡视频| 亚洲三级黄色毛片| 中文欧美无线码| 亚洲欧美清纯卡通| 麻豆成人av视频| 亚洲一级一片aⅴ在线观看| 久久99热6这里只有精品| 亚洲av二区三区四区| 国国产精品蜜臀av免费| 亚洲美女搞黄在线观看| 一区在线观看完整版| 丝袜在线中文字幕| 日本午夜av视频| 免费看日本二区| 日日啪夜夜撸| 国产美女午夜福利| 男女啪啪激烈高潮av片| 国产老妇伦熟女老妇高清| 99精国产麻豆久久婷婷| 18禁在线无遮挡免费观看视频| 五月开心婷婷网| 青青草视频在线视频观看| 午夜影院在线不卡| 少妇被粗大猛烈的视频| 国产精品嫩草影院av在线观看| 精品酒店卫生间| 大码成人一级视频| 国产精品蜜桃在线观看| 我的老师免费观看完整版| 人人澡人人妻人| 日本欧美视频一区| 3wmmmm亚洲av在线观看| 亚洲人与动物交配视频| 欧美日韩一区二区视频在线观看视频在线| 精品少妇内射三级| 国产在线免费精品| 国产综合精华液| 久久这里有精品视频免费| 免费黄频网站在线观看国产| 91aial.com中文字幕在线观看| 在现免费观看毛片| 成人免费观看视频高清| 久久久久人妻精品一区果冻| 久久国产精品男人的天堂亚洲 | 色婷婷久久久亚洲欧美| 九草在线视频观看| 亚洲av电影在线观看一区二区三区| 日本91视频免费播放| 哪个播放器可以免费观看大片| 中文乱码字字幕精品一区二区三区| 欧美激情极品国产一区二区三区 | 一个人免费看片子| 欧美日韩视频精品一区| 高清欧美精品videossex| 一级毛片黄色毛片免费观看视频| 午夜久久久在线观看| 在线观看人妻少妇| 视频区图区小说| 热re99久久精品国产66热6| 亚洲电影在线观看av| 91久久精品国产一区二区成人| 免费看不卡的av| 亚洲国产欧美在线一区| 日韩一区二区三区影片| 全区人妻精品视频| 久久 成人 亚洲| 亚洲,欧美,日韩| 亚洲情色 制服丝袜| 久久婷婷青草| 精品久久国产蜜桃| 男人添女人高潮全过程视频| 日本vs欧美在线观看视频 | 欧美三级亚洲精品| 欧美日本中文国产一区发布| 精品国产乱码久久久久久小说| 青春草国产在线视频| 国产黄色免费在线视频| 简卡轻食公司| 制服丝袜香蕉在线| 人妻人人澡人人爽人人| 久久精品熟女亚洲av麻豆精品| 97超视频在线观看视频| 亚洲av在线观看美女高潮| 国内少妇人妻偷人精品xxx网站| 亚洲四区av| av.在线天堂| 久久97久久精品| 欧美+日韩+精品| 国产日韩一区二区三区精品不卡 | 国产片特级美女逼逼视频| 人体艺术视频欧美日本| 热re99久久国产66热| 日本黄色日本黄色录像| 在线观看美女被高潮喷水网站| 日日摸夜夜添夜夜添av毛片| 多毛熟女@视频| 国产在线男女| 色婷婷av一区二区三区视频| 99热这里只有精品一区| 观看av在线不卡| 极品教师在线视频| 国产亚洲5aaaaa淫片| 亚洲精品成人av观看孕妇| 久久热精品热| 99热全是精品| 久久久午夜欧美精品| 少妇熟女欧美另类| freevideosex欧美| 免费在线观看成人毛片| 美女大奶头黄色视频| 亚洲人与动物交配视频| 纵有疾风起免费观看全集完整版| 麻豆成人午夜福利视频| av.在线天堂| 国产在视频线精品| 国产探花极品一区二区| 日本爱情动作片www.在线观看| 大又大粗又爽又黄少妇毛片口| 国产男女内射视频| 精品一品国产午夜福利视频| 国产亚洲午夜精品一区二区久久| 国产免费又黄又爽又色| 一边亲一边摸免费视频| 精品一区二区三卡| 色哟哟·www| 久久婷婷青草| 成人特级av手机在线观看| 亚洲三级黄色毛片| 一级a做视频免费观看| 中文在线观看免费www的网站| 中国三级夫妇交换| 在线观看人妻少妇| 男女边吃奶边做爰视频| 国产一区有黄有色的免费视频| 国产亚洲91精品色在线| 色5月婷婷丁香| av卡一久久| 日本黄大片高清| 黄色视频在线播放观看不卡| 欧美日韩在线观看h| 国产精品久久久久久精品古装| 久久久久久伊人网av| 亚洲丝袜综合中文字幕| 国产精品免费大片| 十八禁网站网址无遮挡 | 国产成人午夜福利电影在线观看| 久久ye,这里只有精品| 十八禁网站网址无遮挡 |