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

    基于邊界元法的CAE軟件開發(fā)

    2025-02-20 00:00:00任永真劉思凡李子昊
    科技創(chuàng)新與應(yīng)用 2025年5期
    關(guān)鍵詞:軟件開發(fā)可視化

    摘" 要:該文基于Python和Qt框架研發(fā)一款具有前后處理模塊的二維邊界元法CAE軟件,為用戶打造一個(gè)直觀的可視化操作界面,通過該界面可實(shí)現(xiàn)幾何模型構(gòu)建、單元結(jié)點(diǎn)信息獲取、邊界條件加載和材料屬性施加等前處理功能,以及顯示位移、應(yīng)力云圖的后處理功能,可直觀展現(xiàn)邊界元法計(jì)算結(jié)構(gòu)獲得的物理參量。最后以含切口薄板結(jié)構(gòu)受力分析為例,利用所開發(fā)的軟件進(jìn)行計(jì)算,展現(xiàn)所開發(fā)軟件的功效。

    關(guān)鍵詞:邊界元法;可視化;軟件開發(fā);Python;CAE軟件

    中圖分類號(hào):TP391.7" " " 文獻(xiàn)標(biāo)志碼:A" " " " " 文章編號(hào):2095-2945(2025)05-0054-04

    Abstract: In this paper, a two-dimensional boundary element method CAE software with pre-and post-processing modules is developed based on Python and Qt framework, which creates an intuitive visual operation interface for users. Through this interface, pre-processing functions such as geometric model construction, element node information acquisition, and boundary conditions can be realized. Pre-processing functions such as loading and material property application, as well as post-processing functions to display displacement and stress nephograms, can intuitively display the physical parameters obtained by the boundary element method calculation structure. Finally, this paper takes the stress analysis of thin plate structures with notches as an example, and uses the developed software to calculate, which demonstrates the effectiveness of the developed software.

    Keywords: boundary element method; visualization; software development; Python; CAE software

    相對(duì)有限元法商業(yè)化程度高、軟件種類繁多而言,邊界元法的通用軟件較少[1]。現(xiàn)有較多的是邊界元法應(yīng)用程序,數(shù)據(jù)要靠手動(dòng)輸入或文件導(dǎo)入,計(jì)算完成之后,結(jié)果也是以文件的形式給出,需要用戶在文件中尋找需要的信息。邊界元法程序計(jì)算得到的結(jié)果并不能直觀展示結(jié)構(gòu)相應(yīng)的物理參量,僅有的數(shù)值結(jié)果對(duì)于未經(jīng)過培訓(xùn)學(xué)習(xí)的工程技術(shù)人員和研究人員來說難以閱讀。因此邊界元法CAE軟件的開發(fā)還處于發(fā)展的起步階段。

    高宇等[2]基于UG軟件開發(fā)了針對(duì)二維邊界元法的前處理模塊。邊煜凱[3]基于Qt開發(fā)了采用間接邊界元法的地震模擬分析軟件。衛(wèi)軍朝等[4]基于MATLAB環(huán)境介紹了一種CAE軟件的前置處理、解析計(jì)算及后置處理程序?qū)崿F(xiàn)的方法和編程技巧。秦曉晨等[5]基于Python和Abaqus研發(fā)了一款應(yīng)用于無加強(qiáng)U形波紋管的力學(xué)性能分析軟件。這些軟件多數(shù)是為了解決某一特定問題而編制,缺乏通用性。

    本文基于Python開發(fā)了一款二維邊界元力學(xué)分析軟件,實(shí)現(xiàn)在UI界面輸入結(jié)構(gòu)計(jì)算信息、展示力學(xué)分析結(jié)果。軟件開發(fā)完后使用PyInstaller庫將其打包為Python應(yīng)用程序[6],方便用戶通過.exe文件運(yùn)行,界面友好,操作方便。

    1" 軟件概述

    1.1" 軟件開發(fā)環(huán)境及應(yīng)用工具

    本軟件基于二維常規(guī)邊界元法,在Windows10系統(tǒng)下PyCharm環(huán)境中開發(fā),針對(duì)結(jié)構(gòu)的彈性力學(xué)問題進(jìn)行幾何建模、賦予材料屬性、施加邊界條件,再進(jìn)行邊界元法計(jì)算,并將計(jì)算得到的位移和應(yīng)力以云圖的形式展示。

    在界面開發(fā)方面,采用Python語言,使用PyQt工具包創(chuàng)建桌面GUI應(yīng)用程序,構(gòu)建了前處理界面和后處理界面。前處理界面借助OpenGL技術(shù),能夠流暢地渲染2D圖形,為用戶提供清晰的視覺展示。后處理界面則充分利用numpy庫進(jìn)行高效的矩陣運(yùn)算,確保數(shù)據(jù)處理的準(zhǔn)確性和速度。此外軟件還集成matplotlib庫,并采用scipy庫中的griddata函數(shù)進(jìn)行插值處理,用于繪制位移、應(yīng)力云圖。

    軟件的核心邊界元法使用Fortran語言編寫,使用Python的subprocess模塊以及os庫對(duì)數(shù)據(jù)和算法程序進(jìn)行文件操作。

    1.2" 軟件設(shè)計(jì)流程圖

    2" 軟件設(shè)計(jì)

    2.1" 軟件界面設(shè)計(jì)

    軟件的界面先由QTDesigner工具進(jìn)行布局設(shè)置,窗體、按鈕、文字框等界面元素的大小、位置設(shè)置均在工具中完成,然后通過PyUIC將QTDesigner生成的UI文件轉(zhuǎn)成Python源碼,之后將Python源碼添加在main文件中進(jìn)行編輯,用以修改參數(shù)以及添加交互功能。

    在main文件中創(chuàng)建Ui_MainWindow類,該類繼承自QMainWindow,負(fù)責(zé)構(gòu)建主窗體,設(shè)計(jì)窗體布局及按鈕等元素,同時(shí)可以添加如setText方法(設(shè)置文本信息)、clicked方法(點(diǎn)擊交互)等參數(shù)和交互功能。

    在窗體界面中使用OpenGL Widget控件來顯示繪圖過程以及后處理云圖。首先在openglwidgettest文件中創(chuàng)建MyGLWidget1類,該類繼承自QOpenGLWidget,以利用OpenGL功能進(jìn)行圖形渲染。之后在Ui_MainWindow類中,實(shí)例化MyGLWidget1類的對(duì)象,以實(shí)現(xiàn)其圖形渲染函數(shù)的調(diào)用。

    2.2" 前處理模塊設(shè)計(jì)

    2.2.1" 幾何模型構(gòu)建

    該模塊用于生成結(jié)構(gòu)的幾何模型。首先在Ui_MainWindow類中的lineEdit控件中獲取一條邊的參數(shù),將這些參數(shù)信息保存至datasave.Py文件中。通過MyGLWidget1類的對(duì)象,調(diào)取該類中的繪圖函數(shù)drawfunc,函數(shù)讀取datasave文件以渲染二維圖形,并在主窗體的右側(cè)區(qū)域進(jìn)行顯示。

    用戶可通過左側(cè)文本框輸入結(jié)點(diǎn)編號(hào),以逆時(shí)針方向繪制二維任意多邊形,系統(tǒng)會(huì)根據(jù)用戶輸入的2個(gè)坐標(biāo)來定義一條邊,并將這條邊均勻劃分為N個(gè)單元,同時(shí)生成N+1個(gè)邊界結(jié)點(diǎn)。生成的邊界結(jié)點(diǎn)顯示如圖3所示。

    2.2.2" 單元結(jié)點(diǎn)信息構(gòu)建

    該模塊針對(duì)前一個(gè)模塊輸入的幾何參數(shù)信息自動(dòng)生成邊界結(jié)點(diǎn)信息和單元信息。Ui_MainWindow類中的pushButton控件被點(diǎn)擊后會(huì)調(diào)用datasave文件中的結(jié)點(diǎn)和單元信息保存函數(shù),用以按照特定數(shù)據(jù)形式生成和保存單元結(jié)點(diǎn)信息。此外,lineEdit控件會(huì)獲取用戶輸入的內(nèi)結(jié)點(diǎn)坐標(biāo)和編號(hào)信息,并將其保存至datasave文件中,進(jìn)行對(duì)內(nèi)結(jié)點(diǎn)的繪制和存儲(chǔ)。

    邊界結(jié)點(diǎn)信息和單元信息程序會(huì)按照邊界元法要求的數(shù)據(jù)格式自動(dòng)生成并保存。邊界結(jié)點(diǎn)信息獲取如圖4所示。內(nèi)點(diǎn)信息則需要用戶根據(jù)需要手動(dòng)輸入并保存。

    2.2.3" 邊界條件信息構(gòu)建

    該模塊用以獲取結(jié)構(gòu)的位移和應(yīng)力邊界條件。本軟件使用鼠標(biāo)來選擇單元以施加邊界條件,在鼠標(biāo)移動(dòng)過程中,通過Ui_MainWindow類中mouseMoveEvent方法(鼠標(biāo)移動(dòng)事件)捕獲坐標(biāo)信息。當(dāng)mousePressEvent方法(鼠標(biāo)按下事件)被觸發(fā)時(shí),系統(tǒng)會(huì)獲取當(dāng)前的坐標(biāo)信息,并將其與datasave文件中的單元結(jié)點(diǎn)坐標(biāo)進(jìn)行對(duì)比,判定該單元是否被鼠標(biāo)點(diǎn)擊,并對(duì)已點(diǎn)擊單元及其結(jié)點(diǎn)進(jìn)行加粗顯示。當(dāng)mouseReleaseEvent方法(鼠標(biāo)釋放事件)發(fā)生,即鼠標(biāo)抬起時(shí),系統(tǒng)會(huì)進(jìn)行加粗單元和結(jié)點(diǎn)的繪制,此時(shí)點(diǎn)擊保存按鈕即pushButton控件,相應(yīng)函數(shù)會(huì)記錄這些加粗結(jié)點(diǎn)的坐標(biāo)和編號(hào)等信息。此外,lineEdit控件會(huì)讀取用戶輸入的數(shù)據(jù)信息,并與鼠標(biāo)選定的結(jié)點(diǎn)信息相結(jié)合,從而構(gòu)成結(jié)構(gòu)的完整邊界條件。位移邊界條件的施加如圖5所示。

    2.2.4" 材料參數(shù)信息設(shè)置

    該模塊用來獲取結(jié)構(gòu)的基本參數(shù)信息,例如結(jié)構(gòu)對(duì)稱性、彈性模量、泊松比等。選擇類數(shù)據(jù)在Ui_MainWindow類中的comboBox控件獲取,輸入類數(shù)據(jù)在lineEdit控件中獲取。點(diǎn)擊保存參數(shù)的pushButton控件會(huì)調(diào)用數(shù)據(jù)保存函數(shù)對(duì)用戶在前處理模塊輸入的數(shù)據(jù)信息進(jìn)行保存,點(diǎn)擊輸出信息文檔按鈕則會(huì)調(diào)用datasave文件中的write函數(shù),通過Python的文件操作根據(jù)相應(yīng)格式生成.dat文件。材料參數(shù)信息設(shè)置如圖6所示。

    2.3" 邊界元法計(jì)算模塊設(shè)計(jì)

    邊界元法程序使用Fortran語言編寫完成,主要分為3個(gè)步驟:

    1)根據(jù)格式要求將幾何模型前處理信息輸入到.dat文件中。

    2)運(yùn)行邊界元Fortran算法程序,經(jīng)過編譯構(gòu)建.exe文件,執(zhí)行.exe文件,生成結(jié)果.dat文件。

    3)根據(jù)需求從結(jié)果.dat文件中提取數(shù)據(jù)。

    本軟件針對(duì)邊界元法程序的計(jì)算執(zhí)行特點(diǎn),先將用戶在前處理模塊的數(shù)據(jù)輸入板塊輸入的信息進(jìn)行保存,并按照相應(yīng)格式生成.dat文件。通過Python的subprocess模塊對(duì)邊界元法程序生成的.exe文件進(jìn)行調(diào)用執(zhí)行,此時(shí)會(huì)生成結(jié)果.dat文件。

    用戶可以查看當(dāng)前生成的結(jié)果.dat文件,對(duì)感興趣的結(jié)果進(jìn)行查閱,也可以在后處理模塊進(jìn)行交互操作,軟件會(huì)在全局生成多個(gè)內(nèi)點(diǎn),之后調(diào)用.exe文件重新計(jì)算內(nèi)點(diǎn)位移和應(yīng)力并讀取結(jié)果.dat文件,將結(jié)果經(jīng)過插值計(jì)算,呈現(xiàn)位移、應(yīng)力云圖。

    2.4" 后處理模塊設(shè)計(jì)

    軟件的后處理部分主要以云圖的形式展示邊界元法計(jì)算后各邊界點(diǎn)和內(nèi)點(diǎn)的位移或應(yīng)力大小。

    云圖繪制部分使用numpy庫,在10×10的矩形區(qū)域內(nèi)生成500×500個(gè)點(diǎn),并將這些點(diǎn)的X坐標(biāo)和Y坐標(biāo)通過numpy.meshgrid函數(shù)分別存儲(chǔ)起來。系統(tǒng)會(huì)對(duì)全局生成的所有點(diǎn)進(jìn)行篩選,每隔20個(gè)點(diǎn)進(jìn)行一次判斷。這里采用射線法[7]來確定點(diǎn)是否位于多邊形內(nèi)部。如果點(diǎn)位于多邊形內(nèi),系統(tǒng)會(huì)記錄該點(diǎn)的編號(hào)和坐標(biāo)信息,這些信息將用于后續(xù)的結(jié)點(diǎn)計(jì)算。借助scipy庫的griddata函數(shù),能夠?qū)Y(jié)構(gòu)其余邊界結(jié)點(diǎn)和生成的內(nèi)點(diǎn)(這些點(diǎn)都具有位移和應(yīng)力數(shù)據(jù))進(jìn)行插值處理,從而生成一個(gè)包含多邊形內(nèi)部所有點(diǎn)數(shù)據(jù)的矩陣。最后,根據(jù)該矩陣和邊界元法計(jì)算結(jié)果,利用matplotlib在QOpenGLWidget中繪制出相應(yīng)的云圖。

    3" 應(yīng)用舉例

    使用本軟件對(duì)如圖7所示的含切口薄板結(jié)構(gòu)進(jìn)行力學(xué)分析。該模型下側(cè)邊界均為固定鉸支座約束,上側(cè)邊界為?滓=2.0 MPa。材料的彈性模量E=3.9 GPa、泊松比v=0.25。

    4" 結(jié)束語

    本文介紹了基于Python和Qt開發(fā)的二維邊界元法CAE軟件。該軟件的前處理模塊允許用戶通過直觀的可視化界面構(gòu)建幾何模型、獲取單元結(jié)點(diǎn)信息、施加邊界條件和設(shè)置材料參數(shù)信息,并導(dǎo)出數(shù)據(jù)供邊界元法程序計(jì)算使用,經(jīng)邊界元法計(jì)算后,結(jié)果可在后處理界面中以位移和應(yīng)力云圖進(jìn)行全局可視化展示,也可在.dat文件中對(duì)感興趣的數(shù)據(jù)進(jìn)行讀取。文章還詳細(xì)描述了如何構(gòu)建軟件界面和構(gòu)建前后處理模塊。軟件界面友好,操作便捷,為用戶提供了一種直觀、高效的力學(xué)分析方法。

    參考文獻(xiàn):

    [1] 姚振漢,段小華,尹欣,等.邊界元法軟件及其在工程與教學(xué)中的應(yīng)用[J].重慶建筑大學(xué)學(xué)報(bào),2000(6):84-87.

    [2] 高宇,曹悅,劉彪,等.基于邊界元法的CAE軟件開發(fā)及應(yīng)用[J].武漢大學(xué)學(xué)報(bào)(工學(xué)版),2023,56(9):1066-1072.

    [3] 邊煜凱.三維復(fù)雜場(chǎng)地地震動(dòng)放大效應(yīng)并行邊界元模擬及軟件開發(fā)[D].天津:天津城建大學(xué),2022.

    [4] 衛(wèi)軍朝,張廣鵬,方英武,等.基于Matlab的CAE軟件實(shí)現(xiàn)方法[J].制造業(yè)自動(dòng)化,2004(6):51-53.

    [5] 秦曉晨,劉靜,何進(jìn),等.基于Python和Abaqus的波紋管力學(xué)性能分析軟件開發(fā)[J].重型機(jī)械,2023(5):77-84.

    [6] 王瑞文,王圣輝,邊潤根,等.Python程序打包成exe可執(zhí)行文件的方法探究[J].無線互聯(lián)科技,2017(12):52-53.

    [7] 章磊,何芬,李鴻赟.一種基于奇異射線法檢測(cè)點(diǎn)在多邊形內(nèi)的方法[J].計(jì)算機(jī)應(yīng)用研究,2020,37(S2):133-135.

    猜你喜歡
    軟件開發(fā)可視化
    基于CiteSpace的足三里穴研究可視化分析
    基于Power BI的油田注水運(yùn)行動(dòng)態(tài)分析與可視化展示
    云南化工(2021年8期)2021-12-21 06:37:54
    基于CGAL和OpenGL的海底地形三維可視化
    “融評(píng)”:黨媒評(píng)論的可視化創(chuàng)新
    信息安全環(huán)境下的計(jì)算機(jī)軟件開發(fā)研究
    建模技術(shù)在軟件開發(fā)中的應(yīng)用
    核電廠定期試驗(yàn)管理軟件開發(fā)
    科技視界(2016年21期)2016-10-17 20:50:03
    FPGA軟件開發(fā)過程中編碼規(guī)范的作用
    科技視界(2016年20期)2016-09-29 13:35:10
    日韩中文字幕欧美一区二区 | 激情视频va一区二区三区| 婷婷色av中文字幕| 在线亚洲精品国产二区图片欧美| av国产久精品久网站免费入址| 亚洲精品国产色婷婷电影| 国产亚洲欧美精品永久| 精品久久久久久电影网| 亚洲精品自拍成人| 亚洲成人av在线免费| 亚洲中文av在线| 午夜福利在线观看免费完整高清在| 亚洲天堂av无毛| 精品一区在线观看国产| 一边摸一边做爽爽视频免费| 亚洲国产精品999| 一级片免费观看大全| 国产精品人妻久久久影院| 亚洲成人手机| 午夜精品国产一区二区电影| 大码成人一级视频| 久久热在线av| 免费av中文字幕在线| 国产爽快片一区二区三区| 国产免费一区二区三区四区乱码| 两个人看的免费小视频| 免费不卡的大黄色大毛片视频在线观看| 久久久久久伊人网av| 欧美日韩精品成人综合77777| 又大又黄又爽视频免费| 卡戴珊不雅视频在线播放| 老汉色∧v一级毛片| 免费看不卡的av| 老鸭窝网址在线观看| xxx大片免费视频| 男男h啪啪无遮挡| 久久久久久久国产电影| 亚洲成人av在线免费| 日本wwww免费看| 国产综合精华液| 午夜激情久久久久久久| 美女午夜性视频免费| 自线自在国产av| 韩国高清视频一区二区三区| 日韩一区二区三区影片| 久久久久国产精品人妻一区二区| 亚洲精品成人av观看孕妇| 国产片特级美女逼逼视频| 蜜桃国产av成人99| 久久久久久久亚洲中文字幕| 午夜影院在线不卡| 国产男人的电影天堂91| 国产精品人妻久久久影院| 亚洲国产欧美在线一区| 91久久精品国产一区二区三区| 国产成人精品婷婷| 9色porny在线观看| 国产精品久久久久久久久免| 亚洲久久久国产精品| 亚洲精品自拍成人| 国产av国产精品国产| 如何舔出高潮| 美女脱内裤让男人舔精品视频| 国产国语露脸激情在线看| 免费观看无遮挡的男女| 在线观看人妻少妇| 天天躁夜夜躁狠狠久久av| 少妇的丰满在线观看| 夜夜骑夜夜射夜夜干| 日本91视频免费播放| 一级毛片黄色毛片免费观看视频| 2018国产大陆天天弄谢| 青春草国产在线视频| 午夜免费男女啪啪视频观看| 欧美亚洲 丝袜 人妻 在线| 国产一区亚洲一区在线观看| 久久精品国产自在天天线| 国产精品久久久久久久久免| 国产精品免费视频内射| 最近的中文字幕免费完整| 亚洲精品在线美女| 国产精品久久久av美女十八| 国产成人午夜福利电影在线观看| 黄片无遮挡物在线观看| 好男人视频免费观看在线| 极品少妇高潮喷水抽搐| 欧美精品一区二区大全| 欧美少妇被猛烈插入视频| www日本在线高清视频| 超碰97精品在线观看| 在线观看一区二区三区激情| 欧美av亚洲av综合av国产av | av在线老鸭窝| 精品少妇黑人巨大在线播放| 边亲边吃奶的免费视频| 久久久精品国产亚洲av高清涩受| 美女xxoo啪啪120秒动态图| 久久影院123| 99九九在线精品视频| 老司机影院成人| 亚洲欧美清纯卡通| 麻豆精品久久久久久蜜桃| 一边摸一边做爽爽视频免费| 亚洲国产色片| 精品国产一区二区三区四区第35| av在线观看视频网站免费| 我的亚洲天堂| 中文字幕亚洲精品专区| 色94色欧美一区二区| 亚洲av男天堂| 最近的中文字幕免费完整| 久久av网站| 69精品国产乱码久久久| 欧美人与性动交α欧美软件| 青春草视频在线免费观看| 一级片免费观看大全| 成人18禁高潮啪啪吃奶动态图| 久热久热在线精品观看| 伦理电影大哥的女人| 日韩一卡2卡3卡4卡2021年| 午夜福利影视在线免费观看| 亚洲三级黄色毛片| 91aial.com中文字幕在线观看| av在线播放精品| 欧美成人午夜免费资源| 边亲边吃奶的免费视频| 观看美女的网站| 亚洲精品国产av蜜桃| 国产av精品麻豆| xxx大片免费视频| 久久久久国产精品人妻一区二区| 中文字幕人妻丝袜一区二区 | 精品少妇久久久久久888优播| 国产精品久久久久久久久免| 欧美激情高清一区二区三区 | 1024视频免费在线观看| 美女xxoo啪啪120秒动态图| 老汉色∧v一级毛片| 欧美最新免费一区二区三区| www.熟女人妻精品国产| 亚洲图色成人| 久久人妻熟女aⅴ| 亚洲精品国产一区二区精华液| 青春草国产在线视频| 美女高潮到喷水免费观看| 国产国语露脸激情在线看| 亚洲欧美清纯卡通| 99国产精品免费福利视频| 在线看a的网站| 男女午夜视频在线观看| 桃花免费在线播放| 中文字幕色久视频| 精品一区在线观看国产| 人妻人人澡人人爽人人| 亚洲欧美一区二区三区久久| 亚洲人成电影观看| 青草久久国产| 一区在线观看完整版| 中文欧美无线码| 97人妻天天添夜夜摸| 侵犯人妻中文字幕一二三四区| 999精品在线视频| 亚洲少妇的诱惑av| 在线精品无人区一区二区三| av网站免费在线观看视频| 精品一区二区三区四区五区乱码 | 街头女战士在线观看网站| av卡一久久| 国产色婷婷99| av免费观看日本| 日韩免费高清中文字幕av| 在线观看免费视频网站a站| 香蕉国产在线看| 国产免费一区二区三区四区乱码| 国产毛片在线视频| 一本色道久久久久久精品综合| 亚洲欧美精品综合一区二区三区 | 久久ye,这里只有精品| 好男人视频免费观看在线| 午夜91福利影院| 亚洲熟女精品中文字幕| 99热国产这里只有精品6| 国产国语露脸激情在线看| 国产成人精品婷婷| 99久久综合免费| 欧美成人午夜精品| 精品视频人人做人人爽| 国产福利在线免费观看视频| 在线观看免费高清a一片| 欧美在线黄色| 亚洲一区二区三区欧美精品| 久久韩国三级中文字幕| h视频一区二区三区| 欧美日韩精品成人综合77777| 亚洲四区av| 欧美激情 高清一区二区三区| 成人国产麻豆网| 下体分泌物呈黄色| 国产欧美亚洲国产| 日韩av免费高清视频| 美女视频免费永久观看网站| 精品人妻偷拍中文字幕| 丝袜美足系列| 久久av网站| 十分钟在线观看高清视频www| 久久久久国产精品人妻一区二区| 黄频高清免费视频| 亚洲精品av麻豆狂野| 国产淫语在线视频| 99热网站在线观看| 丝袜在线中文字幕| 91午夜精品亚洲一区二区三区| 卡戴珊不雅视频在线播放| 在线观看免费视频网站a站| 精品国产乱码久久久久久小说| 老司机亚洲免费影院| 亚洲成av片中文字幕在线观看 | 国产精品久久久久久精品电影小说| 欧美日本中文国产一区发布| 久久ye,这里只有精品| 亚洲少妇的诱惑av| 最黄视频免费看| 男人爽女人下面视频在线观看| 在线天堂中文资源库| 高清av免费在线| 天天躁夜夜躁狠狠躁躁| 久久99精品国语久久久| 亚洲欧美中文字幕日韩二区| 黑丝袜美女国产一区| 成人漫画全彩无遮挡| 国产一区有黄有色的免费视频| 成人黄色视频免费在线看| 性色avwww在线观看| 亚洲视频免费观看视频| 国产一区亚洲一区在线观看| 亚洲精品国产av蜜桃| a级毛片在线看网站| 国产免费一区二区三区四区乱码| 欧美+日韩+精品| 亚洲av福利一区| 18禁动态无遮挡网站| 美女大奶头黄色视频| 国产1区2区3区精品| 一区二区日韩欧美中文字幕| 亚洲一级一片aⅴ在线观看| 在线观看www视频免费| 在线看a的网站| 久久韩国三级中文字幕| 老汉色av国产亚洲站长工具| 亚洲国产精品一区三区| 欧美变态另类bdsm刘玥| 2018国产大陆天天弄谢| 少妇精品久久久久久久| 亚洲综合色惰| 丝瓜视频免费看黄片| 亚洲精品,欧美精品| 2021少妇久久久久久久久久久| 国产一区二区三区综合在线观看| 国产精品免费大片| av在线老鸭窝| 桃花免费在线播放| 建设人人有责人人尽责人人享有的| 国产有黄有色有爽视频| 黄色毛片三级朝国网站| 国产又色又爽无遮挡免| 各种免费的搞黄视频| 亚洲精品国产av蜜桃| 亚洲av国产av综合av卡| av女优亚洲男人天堂| 久久国产精品大桥未久av| 久久综合国产亚洲精品| 宅男免费午夜| 免费av中文字幕在线| 日本91视频免费播放| 欧美日韩精品成人综合77777| 少妇人妻精品综合一区二区| 边亲边吃奶的免费视频| 9色porny在线观看| 日韩电影二区| 99热网站在线观看| 久久精品国产鲁丝片午夜精品| 制服诱惑二区| 亚洲精品第二区| 亚洲欧美一区二区三区黑人 | 看免费成人av毛片| 日韩大片免费观看网站| 人人澡人人妻人| 久久免费观看电影| 69精品国产乱码久久久| 日本av手机在线免费观看| 999精品在线视频| 亚洲精品久久成人aⅴ小说| 久久99一区二区三区| 精品久久蜜臀av无| 中文乱码字字幕精品一区二区三区| 黑人猛操日本美女一级片| 国产av国产精品国产| 黄网站色视频无遮挡免费观看| 两个人免费观看高清视频| 亚洲激情五月婷婷啪啪| 精品国产一区二区三区久久久樱花| 丝袜美腿诱惑在线| 国产一区二区在线观看av| 久久 成人 亚洲| 久久久久久人妻| 久久精品国产亚洲av高清一级| 国产精品免费大片| 青春草国产在线视频| 色94色欧美一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 在现免费观看毛片| 十八禁高潮呻吟视频| 伦理电影大哥的女人| 日韩一卡2卡3卡4卡2021年| 99久国产av精品国产电影| 蜜桃在线观看..| 久久精品国产自在天天线| 亚洲欧美一区二区三区久久| 中文字幕人妻丝袜一区二区 | 国产成人a∨麻豆精品| 亚洲精品美女久久久久99蜜臀 | 欧美精品国产亚洲| 新久久久久国产一级毛片| 最新中文字幕久久久久| 免费黄网站久久成人精品| 一区在线观看完整版| 亚洲,一卡二卡三卡| 亚洲伊人色综图| 人体艺术视频欧美日本| 性少妇av在线| 久久久久精品性色| 最近最新中文字幕大全免费视频 | 久久99热这里只频精品6学生| 久久亚洲国产成人精品v| 熟女少妇亚洲综合色aaa.| 国产一区有黄有色的免费视频| 国产成人91sexporn| 日韩制服丝袜自拍偷拍| 性少妇av在线| 久久久a久久爽久久v久久| 亚洲综合色惰| 波多野结衣av一区二区av| 日本wwww免费看| 久久人人97超碰香蕉20202| 国产精品一国产av| 欧美日韩亚洲高清精品| 中文字幕人妻熟女乱码| 熟妇人妻不卡中文字幕| 美女主播在线视频| 三上悠亚av全集在线观看| 69精品国产乱码久久久| 18禁观看日本| 极品人妻少妇av视频| www.自偷自拍.com| 日韩制服骚丝袜av| 女人高潮潮喷娇喘18禁视频| 18禁裸乳无遮挡动漫免费视频| 男女午夜视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 极品少妇高潮喷水抽搐| 成年女人毛片免费观看观看9 | 在线 av 中文字幕| 母亲3免费完整高清在线观看 | 精品国产露脸久久av麻豆| 丰满少妇做爰视频| 精品国产露脸久久av麻豆| 精品久久久精品久久久| 国产成人精品久久久久久| 高清欧美精品videossex| 欧美日韩一区二区视频在线观看视频在线| 99久国产av精品国产电影| 天堂中文最新版在线下载| 人妻人人澡人人爽人人| 在线看a的网站| 国产片特级美女逼逼视频| 永久免费av网站大全| 国产老妇伦熟女老妇高清| 免费观看在线日韩| 欧美日韩av久久| 交换朋友夫妻互换小说| 九色亚洲精品在线播放| 亚洲精品,欧美精品| 在线观看www视频免费| 91aial.com中文字幕在线观看| 亚洲成人一二三区av| 亚洲熟女精品中文字幕| 熟妇人妻不卡中文字幕| 色吧在线观看| 日韩精品免费视频一区二区三区| 91久久精品国产一区二区三区| kizo精华| 观看美女的网站| 日韩,欧美,国产一区二区三区| 亚洲色图综合在线观看| 国产高清不卡午夜福利| 超碰成人久久| 久久免费观看电影| 妹子高潮喷水视频| av福利片在线| 亚洲精品一区蜜桃| 久久久久国产一级毛片高清牌| 亚洲第一av免费看| 国产精品免费大片| 久久亚洲国产成人精品v| 人人澡人人妻人| 国产一区二区三区综合在线观看| 色播在线永久视频| 国产成人精品久久久久久| 亚洲国产色片| 亚洲精品,欧美精品| 国产成人欧美| 菩萨蛮人人尽说江南好唐韦庄| videos熟女内射| 各种免费的搞黄视频| 青春草视频在线免费观看| 夫妻午夜视频| 亚洲久久久国产精品| 啦啦啦视频在线资源免费观看| 日韩熟女老妇一区二区性免费视频| 久久午夜福利片| 免费女性裸体啪啪无遮挡网站| 男人爽女人下面视频在线观看| 日韩av在线免费看完整版不卡| 一区二区三区激情视频| 国产国语露脸激情在线看| 国产日韩欧美亚洲二区| 免费久久久久久久精品成人欧美视频| 午夜免费男女啪啪视频观看| 日本爱情动作片www.在线观看| videosex国产| 日本欧美国产在线视频| 男女啪啪激烈高潮av片| 久热久热在线精品观看| 七月丁香在线播放| 精品人妻在线不人妻| 黑人欧美特级aaaaaa片| videosex国产| 亚洲精品日本国产第一区| 亚洲精品视频女| 天堂俺去俺来也www色官网| av福利片在线| 亚洲国产毛片av蜜桃av| 国产精品 欧美亚洲| 18禁国产床啪视频网站| 少妇猛男粗大的猛烈进出视频| 夫妻性生交免费视频一级片| 五月天丁香电影| 国产黄色免费在线视频| 亚洲成国产人片在线观看| 青春草视频在线免费观看| 亚洲欧洲精品一区二区精品久久久 | 日本欧美国产在线视频| 免费在线观看完整版高清| 欧美精品高潮呻吟av久久| 超色免费av| 综合色丁香网| 国产亚洲av片在线观看秒播厂| 亚洲av中文av极速乱| 9191精品国产免费久久| 亚洲男人天堂网一区| 国产一区亚洲一区在线观看| 中文字幕制服av| 午夜日本视频在线| 午夜日韩欧美国产| 亚洲欧美成人精品一区二区| 99热网站在线观看| 飞空精品影院首页| 这个男人来自地球电影免费观看 | 国产精品亚洲av一区麻豆 | 亚洲av免费高清在线观看| 国产欧美亚洲国产| 精品国产国语对白av| 中文字幕人妻丝袜制服| 97人妻天天添夜夜摸| 国产黄频视频在线观看| 三上悠亚av全集在线观看| 日韩制服骚丝袜av| 欧美精品高潮呻吟av久久| 亚洲内射少妇av| 国产不卡av网站在线观看| 精品第一国产精品| 成年人免费黄色播放视频| av卡一久久| 精品国产乱码久久久久久男人| 五月天丁香电影| 1024视频免费在线观看| 美女视频免费永久观看网站| 国产一区二区在线观看av| 欧美在线黄色| 久久久久人妻精品一区果冻| 欧美日本中文国产一区发布| 中文字幕最新亚洲高清| 美女午夜性视频免费| 熟妇人妻不卡中文字幕| 国产成人精品久久久久久| 亚洲,一卡二卡三卡| 少妇被粗大的猛进出69影院| 国产成人免费观看mmmm| 久久久久久人人人人人| 亚洲人成电影观看| 在线免费观看不下载黄p国产| 精品人妻偷拍中文字幕| 欧美精品亚洲一区二区| 国产成人精品无人区| 亚洲国产精品成人久久小说| 纵有疾风起免费观看全集完整版| 我要看黄色一级片免费的| 免费观看av网站的网址| 国产精品久久久久成人av| 在线天堂最新版资源| 我要看黄色一级片免费的| 制服丝袜香蕉在线| 日本欧美国产在线视频| 国产成人精品久久二区二区91 | 久久久久久伊人网av| 黄网站色视频无遮挡免费观看| 在线看a的网站| 国产 精品1| 精品国产一区二区三区久久久樱花| 激情五月婷婷亚洲| av视频免费观看在线观看| 国产精品欧美亚洲77777| 久久这里有精品视频免费| 国产男女内射视频| 久热这里只有精品99| 大片电影免费在线观看免费| 日本91视频免费播放| 女人精品久久久久毛片| 亚洲精品一区蜜桃| av片东京热男人的天堂| 久久精品国产a三级三级三级| 欧美日韩精品成人综合77777| 蜜桃国产av成人99| 哪个播放器可以免费观看大片| 啦啦啦啦在线视频资源| 国产精品偷伦视频观看了| 妹子高潮喷水视频| 欧美日韩精品成人综合77777| 精品99又大又爽又粗少妇毛片| 亚洲欧洲国产日韩| 亚洲国产精品国产精品| 又大又黄又爽视频免费| av线在线观看网站| 一本—道久久a久久精品蜜桃钙片| 久久久久久久亚洲中文字幕| 午夜激情久久久久久久| 欧美变态另类bdsm刘玥| 两个人看的免费小视频| 日韩,欧美,国产一区二区三区| 卡戴珊不雅视频在线播放| 免费高清在线观看日韩| 中文字幕精品免费在线观看视频| 在线看a的网站| 免费黄色在线免费观看| 欧美精品一区二区免费开放| 国产一级毛片在线| 国产亚洲欧美精品永久| 日韩熟女老妇一区二区性免费视频| 国产精品99久久99久久久不卡 | 91在线精品国自产拍蜜月| 天天操日日干夜夜撸| 熟妇人妻不卡中文字幕| 亚洲精品国产色婷婷电影| 叶爱在线成人免费视频播放| 性色av一级| 欧美精品国产亚洲| 国产白丝娇喘喷水9色精品| 97精品久久久久久久久久精品| 国产成人91sexporn| 欧美黄色片欧美黄色片| 五月天丁香电影| 国产极品粉嫩免费观看在线| 午夜免费男女啪啪视频观看| 最近手机中文字幕大全| av免费观看日本| 一区二区日韩欧美中文字幕| av国产久精品久网站免费入址| 在线观看免费日韩欧美大片| 国产日韩欧美在线精品| 一区二区三区四区激情视频| 精品福利永久在线观看| 老熟女久久久| 亚洲情色 制服丝袜| 国产福利在线免费观看视频| 丝瓜视频免费看黄片| 色哟哟·www| 99re6热这里在线精品视频| 好男人视频免费观看在线| 男女下面插进去视频免费观看| 国产男人的电影天堂91| 国产深夜福利视频在线观看| 熟女少妇亚洲综合色aaa.| 十八禁网站网址无遮挡| 欧美 亚洲 国产 日韩一| 只有这里有精品99| 中文字幕制服av| 国产精品人妻久久久影院| 超碰97精品在线观看| 久久久久久伊人网av| 男人添女人高潮全过程视频| 成人18禁高潮啪啪吃奶动态图| 午夜福利在线免费观看网站| 美女国产高潮福利片在线看| 香蕉国产在线看| 精品亚洲成国产av| 精品第一国产精品| 美女脱内裤让男人舔精品视频| 亚洲欧美一区二区三区国产| 亚洲精品aⅴ在线观看| 亚洲精品av麻豆狂野| 亚洲人成77777在线视频| 超碰成人久久| 你懂的网址亚洲精品在线观看| 亚洲五月色婷婷综合| 精品一区二区三区四区五区乱码 |