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

    滑雪板夾角對跳臺滑雪飛行階段氣動(dòng)特性的影響

    2018-08-02 01:32:16張文毅
    體育科學(xué) 2018年7期

    胡 齊,陳 騏,張文毅

    ?

    滑雪板夾角對跳臺滑雪飛行階段氣動(dòng)特性的影響

    胡 齊1,陳 騏1,張文毅2

    1.國家體育總局體育科學(xué)研究所, 北京 100061; 2.北京體育大學(xué), 北京 100084

    目的:分析研究滑雪板夾角對跳臺滑雪飛行階段氣動(dòng)特性的影響。方法:建立精細(xì)化三維幾何模型與網(wǎng)格模型,采用部分時(shí)均(partially averaged Navier-Stokes,PANS)湍流模型對氣動(dòng)特性進(jìn)行計(jì)算流體力學(xué)(computational fluid dynamics,CFD)數(shù)值模擬,獲取不同滑雪板夾角下(24°、28°、32°、36°和40°)力和力矩、壓力分布以及流場形態(tài)。結(jié)果:1)總升力與總阻力、升力系數(shù)與阻力系數(shù)、總力矩與滑雪板力矩均隨滑雪板夾角增大先減小后增大最后小幅減小,但總升阻比基本保持不變,總體呈略微下降趨勢;2)滑雪板升力隨滑雪板夾角增大呈現(xiàn)單調(diào)遞增的趨勢,而滑雪板阻力隨滑雪板夾角增大先減小后增大;3)滑雪板正面壓力分布規(guī)律基本相同且比較均勻,而滑雪板背面下半部分壓力隨滑雪板夾角變化比較明顯;4)運(yùn)動(dòng)員身后以及滑雪板身后的流場形態(tài)主要以回流渦形式出現(xiàn),滑雪板背面回流渦結(jié)構(gòu)明顯比運(yùn)動(dòng)員身后復(fù)雜得多,而且回流渦強(qiáng)度也明顯大得多。結(jié)論:1)滑雪板夾角變化對總升阻比影響很小,但對總升力、總阻力、升力系數(shù)以及阻力系數(shù)產(chǎn)生一定影響,同時(shí)對運(yùn)動(dòng)員氣動(dòng)特性也產(chǎn)生一定影響。2)滑雪板氣動(dòng)特性在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)中起著更為重要作用,不僅要關(guān)注運(yùn)動(dòng)員身體姿態(tài),而且應(yīng)更加關(guān)注滑雪板姿態(tài);3)總力矩大小和滑雪板力矩大小隨滑雪板夾角變化趨勢一致,而且滑雪板力矩在總力矩中始終占據(jù)主導(dǎo)地位,選取合適的滑雪板夾角對飛行階段穩(wěn)定性至關(guān)重要,建議優(yōu)選的滑雪板夾角范圍為24°~32°。

    氣動(dòng)特性;跳臺滑雪飛行階段;滑雪板夾角;計(jì)算流體力學(xué);V型

    1 前言

    整個(gè)跳臺滑雪過程通常依次分為4個(gè)不同的階段:助滑、起跳、飛行和著陸,空氣動(dòng)力學(xué)在這4個(gè)階段中均起著重要的作用。在起跳后,運(yùn)動(dòng)員應(yīng)在飛行階段早期盡快達(dá)到一個(gè)穩(wěn)定的飛行位置,并確保受到小阻力[13],同時(shí)應(yīng)完成向后旋轉(zhuǎn)和向前旋轉(zhuǎn)角動(dòng)量的完全平衡[16]。跳臺滑雪運(yùn)動(dòng)員在0.5s內(nèi)完成飛行姿態(tài)的穩(wěn)定[21],并且在這個(gè)短時(shí)間內(nèi)(約為15 m的飛行距離),很有可能運(yùn)動(dòng)員不能進(jìn)行任何有意的姿態(tài)調(diào)整,可以說起跳階段決定隨后的飛行姿態(tài)。同樣,這也意味著起跳時(shí)所發(fā)生的錯(cuò)誤不能在飛行階段得到有效糾正,但成功起跳的優(yōu)勢可能會(huì)因飛行過程中錯(cuò)誤而喪失。一般來說,飛行階段被認(rèn)為是整個(gè)跳臺滑雪過程中氣動(dòng)特性體現(xiàn)最為明顯的階段。諸多對跳臺滑雪飛行階段氣動(dòng)特性研究文獻(xiàn)[2,5,6,8-10,12,17,18,20]采用了風(fēng)洞試驗(yàn)測量或計(jì)算機(jī)模擬與現(xiàn)場測量相結(jié)合等方法。普遍認(rèn)為,滑雪器材(滑雪板和滑雪服)在跳臺滑雪飛行階段的作用不可忽略。與以往的滑雪板平行飛行方式相比,V型飛行方式具有更好的氣動(dòng)性能,使運(yùn)動(dòng)員能夠更多地向前傾斜[19,20]。近年來,為了使氣流經(jīng)過滑雪板更加平緩,在飛行過程中滑雪板姿態(tài)越來越受到人們的關(guān)注。Virmavirta等[22]通過風(fēng)洞試驗(yàn)研究單個(gè)滑雪板空氣動(dòng)力特性,將滑雪板懸掛在六分量風(fēng)洞平衡測量裝置上,這種試驗(yàn)研究方法或多或少偏理想化,顯然并不是所有的結(jié)果都適用于真實(shí)飛行階段。

    上述這些方法均不能直觀地顯示或分析運(yùn)動(dòng)員周圍的流場。能夠在運(yùn)動(dòng)過程中分析氣動(dòng)力、顯示壓力分布和詳細(xì)流場信息的另一種技術(shù)是計(jì)算流體力學(xué)(CFD),該技術(shù)最近被應(yīng)用于自行車和游泳等運(yùn)動(dòng)中。CFD不需要物理裝置,成本低,可以模擬大多數(shù)物理現(xiàn)象,只需一次模擬就可以收集大量的數(shù)據(jù)。盡管CFD已經(jīng)被認(rèn)為是未來跳臺滑雪空氣動(dòng)力學(xué)研究的重要工具,但迄今為止發(fā)表的CFD相關(guān)研究較少[1,3,4,7,14,15],尤其是針對飛行階段就更少[1,4,7,14]。Keizo等[3]人采用CFD技術(shù)研究了跳雪運(yùn)動(dòng)員在起跳過程中姿態(tài)對氣動(dòng)特性的影響,著重比較兩種不同跳躍者的起跳情況,結(jié)論認(rèn)為,起跳過程中氣動(dòng)力特性在較短的時(shí)間內(nèi)發(fā)生了動(dòng)態(tài)變化,兩種起跳方式的氣動(dòng)力特性差異較大。Gardan等[1]人采用CFD方法來探討迎風(fēng)角和速度對氣動(dòng)力的影響,數(shù)值結(jié)果表明,在早期飛行階段速度對升力系數(shù)和阻力系數(shù)的影響很小,相反,迎風(fēng)角變化對作用于運(yùn)動(dòng)員身體上的升力和阻力有很大的影響。

    綜上所述,國內(nèi)、外報(bào)道姿態(tài)與跳臺滑雪飛行階段氣動(dòng)特性的CFD研究主要集中于運(yùn)動(dòng)員身體姿態(tài),鮮有報(bào)道滑雪板姿態(tài)的CFD研究[7,14],即使有也僅是涉及而已。本研究擬建立運(yùn)動(dòng)員/滑雪板整體系統(tǒng)的精細(xì)化三維幾何模型與網(wǎng)格模型,采用PANS湍流模型進(jìn)行CFD數(shù)值模擬,獲取不同滑雪板夾角下力和力矩、壓力分布以及流場形態(tài),分析研究V型飛行方式滑雪板夾角對跳臺滑雪飛行階段氣動(dòng)特性的影響。

    2 研究對象與方法

    2.1 研究對象

    研究對象為跳臺滑雪運(yùn)動(dòng)員/滑雪板整體系統(tǒng)。根據(jù)Müller等[9]人統(tǒng)計(jì)分析結(jié)果,跳臺滑雪運(yùn)動(dòng)員的身體形態(tài)特征選取平均值,即身高為177 cm,身體質(zhì)量指數(shù)(BMI)為19.5,軀干(坐高)與身高比值0.532,滑雪板長度258 cm,寬11.5 cm。

    2.2 研究方法

    2.2.1 控制方程

    本研究采用大渦模擬技術(shù),該技術(shù)對預(yù)測鈍體周圍的流動(dòng)分離有優(yōu)勢。為了獲得更精確的結(jié)果,采用了PANS湍流模型,其控制方程表達(dá)式如下:

    式中i為已分解流場速度,為時(shí)間,為流體密度,為介質(zhì)粘性系數(shù),為湍流粘性系數(shù),f為未分解湍動(dòng)能比率,為未分解湍動(dòng)能耗散率比率,k為未分解局部時(shí)均化湍動(dòng)能,為未分解局部時(shí)均化湍動(dòng)能耗散率。

    其中:

    采用有限體積法對控制方程進(jìn)行離散計(jì)算,壓力和速度的耦合采用SIMPLEC算法,時(shí)間采用二階差分格式進(jìn)行離散,湍動(dòng)能以及速度項(xiàng)采用二階迎風(fēng)格式,時(shí)間步長為0.000 1 s。

    2.2.2 模型建立與網(wǎng)格無關(guān)性驗(yàn)證

    根據(jù)研究對象特征,建立跳臺滑雪運(yùn)動(dòng)員/滑雪板整體系統(tǒng)的三維實(shí)體模型,對運(yùn)動(dòng)員的身體特征進(jìn)行精細(xì)化建模,手指、耳朵、眼罩、面部特征可以從圖1中清晰分辨。

    圖1 運(yùn)動(dòng)員/滑雪板整體系統(tǒng)精細(xì)化實(shí)體模型圖

    Figure 1. The Refined Solid Model of the Whole Athlete/Skis System

    注:左圖為正視圖,右圖為側(cè)視圖。

    跳臺滑雪運(yùn)動(dòng)員/滑雪板整體系統(tǒng)計(jì)算區(qū)域的尺寸大小為長13.5 m、寬5 m、高6 m,如圖2所示。采用商用前處理軟件ICEM對計(jì)算區(qū)域進(jìn)行四面體網(wǎng)格離散,在計(jì)算區(qū)域中心對稱面上體積網(wǎng)格分布和運(yùn)動(dòng)員表面網(wǎng)格分布如圖3所示。運(yùn)動(dòng)員模型的表面網(wǎng)格分辨率最小達(dá)到0.5 mm。對于體網(wǎng)格,在運(yùn)動(dòng)員周圍區(qū)域配置了較細(xì)的網(wǎng)格單元,以高精度捕獲尾流結(jié)構(gòu)。另一方面,為了降低計(jì)算量,較粗的網(wǎng)格單元被分配到遠(yuǎn)離運(yùn)動(dòng)員的區(qū)域。

    圖2 計(jì)算區(qū)域示意圖

    Figure. 2. Computational Domain

    圖3 網(wǎng)格分布示意圖

    Figure 3. Mesh Distribution

    注:左圖為體網(wǎng)格,右圖為運(yùn)動(dòng)員表面網(wǎng)格。

    為了滿足PANS模型計(jì)算要求,運(yùn)動(dòng)員周邊采用網(wǎng)格加密,并保證y plus控制在100之內(nèi),此控制要求在先前的研究中獲得了證實(shí)[3]。針對上述網(wǎng)格模型,通過均勻加密選定了4種網(wǎng)格密度,網(wǎng)格點(diǎn)數(shù)從800萬~2 357萬,進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證,驗(yàn)證結(jié)果如表1所示。4種網(wǎng)格驗(yàn)證計(jì)算獲得的升阻比均為1.8。由此可見,即使選定800萬網(wǎng)格節(jié)點(diǎn)的計(jì)算域離散方案,也能夠準(zhǔn)確預(yù)測跳臺滑雪運(yùn)動(dòng)員/滑雪板整體系統(tǒng)的氣動(dòng)力學(xué)性能。

    表1 網(wǎng)格無關(guān)性驗(yàn)證結(jié)果

    2.2.3 邊界條件與計(jì)算工況

    邊界條件的設(shè)置如下:

    1)進(jìn)口設(shè)置為速度進(jìn)口,進(jìn)口速度根據(jù)飛行速度設(shè)定;2)出口設(shè)置為壓力出口,壓力為大氣壓101 325 Pa;3)中間截面為周期性條件;4)其他壁面為無滑移邊界條件;5)氣體為不可壓縮空氣;6)在常重力下,重力加速度設(shè)置為g0=9.807 m/s2。

    跳臺滑雪運(yùn)動(dòng)員在飛行階段的姿態(tài)參數(shù)包括迎風(fēng)角、滑雪板與速度方向夾角α、滑雪板與身體夾角θ、上半身彎曲角度β、滑雪板夾角λ,如圖4所示。Müller等[10,17]人的研究成果給出這些姿態(tài)參數(shù)通常的取值范圍。同時(shí),圖5顯示了在飛行階段整個(gè)系統(tǒng)受力情況。阻力系數(shù)(C)和升力系數(shù)(C)的計(jì)算方法見公式(11)和(12),為模型的正面投影區(qū)域面積。

    在本研究中,開展了5種不同滑雪板夾角下跳臺滑雪飛行階段動(dòng)特性和流動(dòng)結(jié)構(gòu)特性CFD研究。其中,滑雪板夾角λ分別取值為24°、28°、32°、36°和40°,同時(shí)迎風(fēng)角φ=30°、滑雪板與速度方向夾角α=30°、滑雪板與身體夾角θ=16°、上半身彎曲角度β=18°,速度V=29 m/s。據(jù)此,為每個(gè)滑雪板夾角λ的運(yùn)動(dòng)員/滑雪板整體系統(tǒng)生成計(jì)算網(wǎng)格,分別進(jìn)行了CFD數(shù)值模擬。

    圖4 飛行階段姿態(tài)參數(shù)示意圖

    Figure 4. Attitude Parameters in Flight

    圖5 飛行階段整體系統(tǒng)受力情況圖

    Figure. 5. All the Forces Acting on the Athlete/Skis System during Flight Phase

    3 結(jié)果

    3.1 力和力矩

    作用于運(yùn)動(dòng)員/滑雪板整體系統(tǒng)上的氣動(dòng)力可分為升力和阻力,而且這些力的作用點(diǎn)絕大多數(shù)不會(huì)在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)質(zhì)心上,因此極有可能會(huì)產(chǎn)生相應(yīng)的力矩。表2列出了不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)以及滑雪板的力學(xué)特性結(jié)果。圖6~10顯示了運(yùn)動(dòng)員/滑雪板整體系統(tǒng)以及滑雪板的力學(xué)特性隨滑雪板夾角變化曲線。圖6和圖7中各項(xiàng)力均為作用在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)或滑雪板上同一性質(zhì)的合力。圖8中阻力系數(shù)與升力系數(shù)分別由公式(11)和(12)計(jì)算得到。圖9中力矩是相對于運(yùn)動(dòng)員/滑雪板整體系統(tǒng)質(zhì)心的力矩,“+”代表力矩作用后使整體系統(tǒng)逆時(shí)針旋轉(zhuǎn),“-”代表力矩作用后使整體系統(tǒng)順時(shí)針旋轉(zhuǎn)。圖10中的升阻比由升力除以阻力計(jì)算得到的。

    不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)總升力以及滑雪板升力隨滑雪板夾角變化曲線如圖6所示??梢钥闯觯偵﹄S滑雪板夾角增大先減小后增大最后小幅減小,即當(dāng)夾角不大于28°時(shí),總升力隨夾角增大而減?。划?dāng)夾角大于28°且小于36°時(shí),總升力隨夾角增大而增大;當(dāng)夾角大于36°時(shí),總升力隨夾角增大小幅減小。滑雪板升力隨滑雪板夾角增大呈現(xiàn)單調(diào)遞增的趨勢,即當(dāng)夾角不大于28°時(shí),升力隨夾角增大而小幅增加;當(dāng)夾角大于28°且小于36°時(shí),升力隨夾角增大而大幅增大;當(dāng)夾角大于36°時(shí),升力隨夾角增大小幅增加。總升力的最大增幅是滑雪板升力的兩倍左右。

    表2 力學(xué)特性結(jié)果

    圖6 升力變化曲線圖

    Figure 6. Lift Force Changes with Ski Opening Angle

    不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)總阻力以及滑雪板阻力隨滑雪板夾角變化曲線如圖7所示??梢钥闯?,總阻力隨滑雪板夾角增大變化趨勢與總升力一致,即先減小后增大最后小幅減小。滑雪板阻力隨滑雪板夾角增大先減小后增大,即當(dāng)夾角不大于28°時(shí),阻力隨夾角增大而減?。划?dāng)夾角大于28°且小于36°時(shí),阻力隨夾角增大而大幅增大;當(dāng)滑雪板夾角大于36°時(shí),阻力隨夾角增大小幅增加??傋枇Φ淖畲笤龇腔┌遄枇Φ膬杀蹲笥摇?/p>

    圖7 阻力變化曲線圖

    Figure 7. Drag Force Changes with Ski Opening Angle

    不同雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)阻力系數(shù)與升力系數(shù)隨滑雪板夾角變化曲線如圖8所示。不難看出,阻力系數(shù)隨滑雪板夾角增大變化趨勢與升力系數(shù)一致,即先減小后增大最后小幅減小,同樣與總升力以及總阻力隨滑雪板夾角增大變化趨勢一致。

    圖8 阻力系數(shù)、升力系數(shù)變化曲線圖

    Figure 8. Changes of Drag Coefficient and Lift Coefficient with Ski Opening Angle

    不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)總力矩以及滑雪板力矩隨滑雪板夾角變化曲線如圖9所示??梢钥闯?,總力矩與滑雪板力矩均是產(chǎn)生逆時(shí)針旋轉(zhuǎn)的效果,其總力矩大小和滑雪板力矩大小隨滑雪板夾角增大變化趨勢一致,即先減小后增大最后小幅減小,同樣與總升力以及總阻力隨滑雪板夾角增大變化趨勢一致。

    圖9 力矩變化曲線圖

    Figure 9. Torque Changes with Ski Opening Angle

    圖10 升阻比變化曲線圖

    Figure 10. Lift-to-drag Ratio Changes with Ski Opening Angle

    不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)總升阻比以及滑雪板升阻比隨滑雪板夾角變化曲線如圖10所示??梢钥闯?,總升阻比隨滑雪板夾角增大變化基本保持不變,總體呈略微下降趨勢。另一方面,滑雪板升阻比隨滑雪板夾角增大先大幅增大(在滑雪板夾角不大于28°之前),之后略微下降。

    3.2 壓力分布

    不同滑雪板夾角λ下運(yùn)動(dòng)員/滑雪板整體系統(tǒng)的正面壓力分布如圖11所示。對于運(yùn)動(dòng)員來說,正面壓力在不同夾角下分布規(guī)律基本一致,高壓區(qū)主要集中于頭部、上半身、胯關(guān)節(jié)、膝關(guān)節(jié)等部位,低壓區(qū)主要集中于身體側(cè)部。同時(shí)不難看出,夾角為28°和32°時(shí)正面高壓區(qū)面積偏小,尤其是滑雪板夾角為28°時(shí)正面高壓區(qū)面積明顯最小。對于滑雪板來說,正面壓力在不同夾角下分布規(guī)律基本相同,而且分布比較均勻。

    圖11 正面壓力分布圖

    Figure 11. Frontal Pressure Distribution

    不同滑雪板夾角λ下滑雪板背面壓力分布如圖12所示??梢钥闯觯哓?fù)壓區(qū)主要集中在滑雪板頂端和中部。進(jìn)一步可看出,滑雪板背面上半部分的壓力隨滑雪板夾角變化比較細(xì)微,然而滑雪板背面下半部分壓力變化比較明顯。

    圖12 滑雪板背面壓力分布圖

    Figure 12. Pressures Distribution of the Ski Back

    3.3 流場形態(tài)

    不同滑雪板夾角λ下渦流形態(tài)和氣流速度流線如圖13所示,左邊為背面正視圖,右邊為側(cè)視圖。運(yùn)動(dòng)員身后以及滑雪板后面的渦流形態(tài)均主要以回流渦結(jié)構(gòu)形式出現(xiàn),此結(jié)構(gòu)會(huì)加速能量耗散。對于運(yùn)動(dòng)員來說,其身后主要產(chǎn)生2對渦結(jié)構(gòu),分別在肩與上臂結(jié)合處以及腰與胯結(jié)合處的后方,當(dāng)夾角為28°時(shí),身后的回流渦結(jié)構(gòu)強(qiáng)度最低。對于滑雪板來說,無論角度是否改變,其后側(cè)有6對渦結(jié)構(gòu)產(chǎn)生,當(dāng)夾角為28°時(shí),后側(cè)的回流渦結(jié)構(gòu)強(qiáng)度也相對最低,之后隨著夾角增大,其后側(cè)的回流渦結(jié)果變得更加無序,尤其是滑雪板下半部分的后側(cè)。不難看出,當(dāng)夾角為24°時(shí),運(yùn)動(dòng)員腰與胯結(jié)合處后方的回流渦結(jié)構(gòu)強(qiáng)度較小,渦流尺度較大的區(qū)域主要集中在滑雪板后側(cè);當(dāng)夾角大于28°時(shí),運(yùn)動(dòng)員身后的回流渦強(qiáng)度高于滑雪板后側(cè),在腰與胯結(jié)合處后方形成尺度較大的回流渦。同時(shí),觀察到在負(fù)壓力梯度作用下,氣流速度在運(yùn)動(dòng)員以及滑雪板上均發(fā)生了流動(dòng)分離,尤其是在滑雪板下半部分。當(dāng)夾角為24°時(shí)流動(dòng)分離程度以及尾跡尺寸明顯最大,當(dāng)夾角為28°時(shí)流動(dòng)分離程度以及尾跡尺寸明顯最小,當(dāng)夾角大于28°時(shí)流動(dòng)分離程度以及尾跡尺寸逐漸變大。

    4 討論

    本研究中,所獲得的氣動(dòng)特性部分結(jié)果與Gardan等[1]人的部分研究結(jié)果相似。通過總結(jié)風(fēng)洞試驗(yàn)測量結(jié)果以及數(shù)值計(jì)算結(jié)果,他們公布了迎風(fēng)角為29°、滑雪板與身體夾角θ為17.4°、上半身彎曲角度β為16.6°、滑雪板夾角λ為37.9°的姿態(tài)條件下總升阻比為1.89。在本研究中,在迎風(fēng)角為30°、滑雪板與身體夾角θ為16°、上半身彎曲角度β為18°、滑雪板夾角λ為40°的姿態(tài)條件下總升阻比為1.818。兩者的結(jié)果非常接近,驗(yàn)證了本研究CFD計(jì)算結(jié)果的正確性。

    從圖6~8可以看出,總升力與總阻力、阻力系數(shù)與升力系數(shù)隨滑雪板夾角增大變化趨勢一致。同時(shí),圖10顯示出總升阻比隨滑雪板夾角變化基本保持不變,總體呈略微下降趨勢,也進(jìn)一步說明總升力和總阻力變化趨勢的一致性。因此,可以認(rèn)為,滑雪板夾角變化對運(yùn)動(dòng)員/滑雪板整體系統(tǒng)升阻比影響很小,但對阻力系數(shù)與升力系數(shù)有一定影響。

    圖13 流場形態(tài)示意圖

    Figure 13. Flow Field Form

    注:左圖為后視圖,右圖為側(cè)視圖。

    N?rstrud等[14]通過數(shù)值計(jì)算研究認(rèn)為,V型飛行方式代表了更為先進(jìn)的空氣動(dòng)力,滑雪板在產(chǎn)生良好的升力方面起著重要作用。從表2、圖6和圖7可以看出,不僅滑雪板升力一直占據(jù)總升力63%以上,而且滑雪板阻力也一直占據(jù)總阻力65%以上,可以認(rèn)為,滑雪板氣動(dòng)特性在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)中起著更為重要作用,滑雪板氣動(dòng)特性的改善將會(huì)獲得更為有利的效果,在分析飛行階段姿態(tài)氣動(dòng)特性時(shí),不僅要關(guān)注運(yùn)動(dòng)員身體姿態(tài),而且應(yīng)更加關(guān)注滑雪板姿態(tài)。

    另外,圖6展示了滑雪板升力隨滑雪板夾角增大呈現(xiàn)單調(diào)遞增的趨勢,同時(shí),圖8顯示了在滑雪板夾角大于28°后升力系數(shù)總體呈增大趨勢,即增大夾角能夠提高總氣動(dòng)升力系數(shù)。同時(shí),總力矩大小和滑雪板力矩大小隨滑雪板夾角變化趨勢一致,而且滑雪板力矩值是總力矩值的1.17倍以上,如圖9和表2所示。當(dāng)夾角大于28°時(shí),隨著夾角增大,總力矩值和滑雪板力矩值越大,產(chǎn)生向前傾的效果越明顯,可能會(huì)對運(yùn)動(dòng)員的穩(wěn)定飛行越不利。本研究結(jié)果與Virmavirta[22]的風(fēng)洞試驗(yàn)研究結(jié)論一致。他們重點(diǎn)研究了單個(gè)滑雪板的空氣動(dòng)力特性,結(jié)論認(rèn)為,大的滑雪板夾角改善了氣動(dòng)升力系數(shù),但是,張開如此寬度的滑雪板夾角的飛行姿態(tài)可能會(huì)使運(yùn)動(dòng)員飛行變得不穩(wěn)定。因此,運(yùn)動(dòng)員應(yīng)根據(jù)自身技術(shù)特點(diǎn)以及訓(xùn)練效果來選擇滑雪板夾角,切忌一味為了改善氣動(dòng)升力系數(shù)而選擇大的滑雪板夾角??偠灾?,滑雪板力矩在總力矩中始終占據(jù)主導(dǎo)地位,選取合適的滑雪板夾角對飛行穩(wěn)定性至關(guān)重要,優(yōu)選的滑雪板夾角范圍為24°~32°。

    隨著滑雪板夾角變化,滑雪板背面下半部分壓力變化比較明顯,同時(shí)滑雪板正面壓力基本不變且滑雪板背面上半部分的壓力變化也比較細(xì)微(圖11、12),這也是導(dǎo)致滑雪板升力、阻力以及力矩變化的主要原因。當(dāng)夾角為28°時(shí),滑雪板背面下半部分的次高負(fù)壓區(qū)域明顯偏少,低負(fù)壓區(qū)域明顯偏多,尤其在滑雪板底端部位,故此時(shí)滑雪板正面背面壓差偏小,從而滑雪板升力、阻力以及向前傾的力矩值偏小,尤其是力矩值,這也從表2以及圖6~9中不難看出。

    運(yùn)動(dòng)員身后以及滑雪板后側(cè)的流場形態(tài)主要以回流渦形式出現(xiàn),回流渦結(jié)構(gòu)會(huì)使能量耗散速度加快。如圖13所示,滑雪板后側(cè)回流渦結(jié)構(gòu)明顯比運(yùn)動(dòng)員身后復(fù)雜得多,而且綜合尺度及強(qiáng)度也明顯大得多,說明滑雪板背面能量損失更大,導(dǎo)致滑雪板的升力和阻力數(shù)值較大,并且在總升力和總阻力中占據(jù)較大的比重,從而進(jìn)一步驗(yàn)證了滑雪板氣動(dòng)特性在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)中的確起著更為重要作用。同時(shí),對于運(yùn)動(dòng)員來說,當(dāng)夾角為24°時(shí)運(yùn)動(dòng)員身后回流渦結(jié)構(gòu)尺度以及流動(dòng)分離程度明顯偏大,即此時(shí)運(yùn)動(dòng)員的升力與阻力偏大;當(dāng)夾角為28°時(shí)運(yùn)動(dòng)員身后回流渦結(jié)構(gòu)尺度以及流動(dòng)分離程度明顯最小,即此時(shí)運(yùn)動(dòng)員的升力與阻力最小,當(dāng)夾角大于28°時(shí)運(yùn)動(dòng)員身后回流渦結(jié)構(gòu)尺度以及流動(dòng)分離程度逐漸變大,直至夾角為40°時(shí)運(yùn)動(dòng)員身后回流渦結(jié)構(gòu)尺度以及流動(dòng)分離程度又略微變小一些,即運(yùn)動(dòng)員的升力與阻力先變大后略微變小。從表2以及圖6~9中不難看出,運(yùn)動(dòng)員的升力和阻力隨滑雪板夾角變化趨勢也與上述變化趨勢一致,這兩種結(jié)果之間也起到良好的相互驗(yàn)證效果。這些研究結(jié)果表明,滑雪板夾角變化對運(yùn)動(dòng)員氣動(dòng)特性產(chǎn)生一定影響,運(yùn)動(dòng)員的升力以及阻力的變化趨勢與總升力以及總阻力變化趨勢一致。

    5 結(jié)論

    1.滑雪板夾角變化對運(yùn)動(dòng)員/滑雪板整體系統(tǒng)的升阻比影響很小,但對總升力、總阻力、升力系數(shù)以及阻力系數(shù)均產(chǎn)生一定影響,同時(shí)對運(yùn)動(dòng)員本身的氣動(dòng)特性也產(chǎn)生一定影響??偵涂傋枇﹄S滑雪板夾角的變化趨勢與運(yùn)動(dòng)員的升力和阻力變化趨勢一致。

    2.滑雪板氣動(dòng)特性在運(yùn)動(dòng)員/滑雪板整體系統(tǒng)中起著更為重要作用,滑雪板氣動(dòng)特性的改善將會(huì)獲得更為有利的效果,在分析飛行階段姿態(tài)氣動(dòng)特性時(shí),不僅要關(guān)注運(yùn)動(dòng)員身體姿態(tài),更加要關(guān)注滑雪板姿態(tài)。

    3.總力矩大小和滑雪板力矩大小隨滑雪板夾角變化趨勢一致,而且滑雪板力矩在總力矩中始終占據(jù)主導(dǎo)地位,選取合適的滑雪板夾角對飛行階段穩(wěn)定性至關(guān)重要,運(yùn)動(dòng)員應(yīng)根據(jù)自身技術(shù)特點(diǎn)以及訓(xùn)練效果來選擇滑雪板夾角,切記一味為了改善氣動(dòng)升力系數(shù)而選擇大的滑雪板夾角,建議優(yōu)選的滑雪板夾角范圍為24°~32°。

    [1] GARDAN N, SCHNEIDER A, POLIDORI G,. Numerical investigation of the early flight phase in ski-jumping[J]. J Biomech, 2017, 59(1): 29-34.

    [2] JUNG A, STAAT M, MULLER W. Flight style optimization in ski jumping on normal, large, and ski flying hills[J]. J Biomech, 2014, 47(3): 716-722.

    [3] KEIZO Y, MAKOTO T, JUN I,. Effect of posture on the aerodynamic characteristics during take-off in ski jumping [J]. J Biomech, 2016, 49(15): 3688-3696.

    [4] LEE K D, PARK M J, KIM K Y. Optimization of ski jumper’s posture considering lift-to-drag ratio and stability[J]. J Biomech, 2012, 45(12): 2125–2132.

    [5] MARQUES-BRUNA P, GRIMSHAW P. Mechanics of flight in ski jumping: aero-dynamic stability in pitch[J]. Sports Technol, 2009, 2(1–2): 24–31.

    [6] MARQUES-BRUNA P, GRIMSHAW P. Mechanics of flight in ski jumping: aero-dynamic stability in roll and yaw[J]. Sports Technol, 2009, 2(3-4): 111–120.

    [7] MEILE W, REISENBERGER E, MAYER M,. Aerodynam-ics of ski jumping: experiments and CFD simulations[J]. Exp Fluids, 2006, 41(6): 949–964.

    [8] MULLER W.Performance factors in ski jumping[J]. J Biomech, 2006, 39(1): 192–213.

    [9] MULLER W, GROSCHL W, MULLER R,. Underweight in ski jumping: The solution of the problem[J]. Int J Sports Med, 2006, 27(11): 926-934.

    [10] MULLER W, PLATZER D, SCHMOLZER B. Dynamics of human flight on skis: improvements in safety and fairness in ski jumping[J]. J Biomech, 1996, 29(8): 1061–1068 .

    [11] MULLER W, PLATZER D, SCHMOLZER B. Scientific approach to ski safety[J]. Nature, 1995, 375(6531): 455-466.

    [12] MURAKAMI M, IWASE M, SEO K,.High-speed video image analysis of ski jumping flight posture[J]. Sports Eng, 2014, 17(4): 217–225.

    [13] MURAKAMI M, IWASE M, SEO K,. Ski jumping flight skill analysis based on high-speed video image[J]. Procedia Eng, 2010, 2(2): 2381–2386.

    [14] N?RSTRUD H, ?YE I J. On CFD simulation of ski jumping[M]//Computational Fluid Dynamics for Sport Simulation, 2009, 63–82.

    [15] RYU M, CHO L, CHO J. Aerodynamic analysis on postures of ski jumpers during flight using computational fluid dynamics[J]. TJp Soc Aeronaut S, 2015, 58(4): 204–212.

    [16] SCHWAMEDER H. Biomechanics research in ski jumping: 1991–2006[J]. Sport Biomech, 2008, 7(1): 114–136.

    [17] SCHMOLZER B, MULLER W. Individual flight styles in ski jumping: results obtained during Olympic games competitions[J]. J Biomech, 2005, 38(5): 1055–1065.

    [18] SCHMOLZER B, MULLER W. The importance of being light: aerodynamic forces and weight in ski jumping[J]. J Biomech, 2002, 35(8): 1059–1069.

    [19] SEO K, MURAKAMI M, YOSHIDA K. Optimal flight technique for V-style ski jumping[J]. Sports Eng, 2004, 7(2): 97–104.

    [20] SEO K, WATANABE I, MURAKAMI M. Aerodynamic force data for a V-style ski jumping flight[J]. Sports Eng, 2004, 7(1): 31–39.

    [21] VIRMAVIRTA M, ISOLEHTO J, KOMI P,. Characteristics of the early flight phase in the Olympic ski jumping competition[J]. J Biomech, 2005, 38(11): 2157–2163.

    [22] VIRMAVIRTA M, KIVEKAS J. Aerodynamics of isolated ski jumping ski[C]. XXV ISB Congress, Glasgow, 2015: 215-223.

    Effect of the Ski Opening Angle on the Aerodynamic Characteristics during Flight in Ski-jumping

    HU Qi1, CHEN Qi1, ZHANG Wen-yi2

    1. China Institute of Sport Science, Beijing 100061, China; 2. Beijing Sport University, Beijing 100084, China.

    Objective: This study aimed to investigate the effects of the ski opening angle on aerodynamic characteristics during the flight in ski-jumping. Methods: Refined 3D geometric models and mesh models were established, and aerodynamic characteristics were numerically simulated by using computational fluid dynamics (CFD) with partially averaged Navier-Stokes (PANS) turbulence model. The force and moment, pressure distribution and flow field form of different ski opening angles (24°, 28°, 32°, 36° and 40°) were obtained. Results: 1) As ski opening angle increases ,the total lift force and that drag force, the coefficients of lift and drag, the total torque and the skis’ torque all decrease firstly, then increase, and finally decrease slightly, but the total lift-to-drag ratio almost remains unchanged, and shows a slight downward trend as a whole. 2) The skis’ lift force increases monotonously with the increase of ski opening angle, while the skis’ drag force decreases firstly and then increases. 3) The front pressure distribution of skis is approximately the same and fairly uniform, but the back pressures distribution of the lower part of skis varies obviously with ski opening angle. 4) The flow field form behind athlete and skis mainly appears in the form of reflux vortex. The vortex structure behind skis is much more complicated than that behind the athlete, and the vortex intensity is also much greater. Conclusion: 1) The ski opening angle has little effect on the total lift-to-drag ratio, but has some influence on the total lift and total drag and the coefficients of lift and drag, as well as on the aerodynamic characteristics of the athlete. 2) The aerodynamic characteristics of skis play a more important role in the whole athlete/skis system, so more attention should be paid to skis’ attitude while paying attention to athletes' body posture. 3) The variation trend of total torque and skis’ torque with ski opening angle is consistent, and skis’ torque has always been dominant in the total torque, so it is very important to select the proper ski opening angle for the stability of flight phase, and it is suggested that the optimum range of ski opening angle is 24°~32°.

    1000-677X(2018)07-0042-08

    10.16469/j.css.201807017

    G804.6

    A

    2018-06-09;

    2018-07-12

    國家體育總局體育科學(xué)研究所基本科研業(yè)務(wù)費(fèi)資助項(xiàng)目(18-29)、中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助課題(2018RC016)。

    胡齊,男,助理研究員,碩士,主要研究方向?yàn)轶w育工程、流體與機(jī)械工程、運(yùn)動(dòng)生物力學(xué),E-mail:huqi@ciss.cn; 陳騏,男,副研究員,學(xué)士,主要研究方向?yàn)轶w育工程,E-mail:chenqi@ciss.cn; 張文毅,女,助理研究員,碩士,主要研究方向?yàn)檫\(yùn)動(dòng)生物力學(xué)、高等教育管理, E-mail: zhangwy1228@126.com。

    国产高清国产精品国产三级| a在线观看视频网站| 人妻一区二区av| 久久天堂一区二区三区四区| 国产免费av片在线观看野外av| 亚洲精品国产精品久久久不卡| av超薄肉色丝袜交足视频| 亚洲一区高清亚洲精品| 丁香六月欧美| 亚洲 国产 在线| 日本精品一区二区三区蜜桃| 国产成人影院久久av| 精品久久蜜臀av无| 欧美精品啪啪一区二区三区| 亚洲成人免费av在线播放| 国产一卡二卡三卡精品| 午夜亚洲福利在线播放| 在线看a的网站| 欧美色视频一区免费| 欧美精品av麻豆av| 亚洲成人手机| 久久国产精品影院| 国产免费av片在线观看野外av| 亚洲一码二码三码区别大吗| √禁漫天堂资源中文www| 国产免费男女视频| 午夜福利欧美成人| 亚洲国产欧美一区二区综合| 欧美中文综合在线视频| 国产亚洲av高清不卡| 在线观看免费视频网站a站| 精品免费久久久久久久清纯 | 国产免费男女视频| 一a级毛片在线观看| 91av网站免费观看| 欧美成人免费av一区二区三区 | 99re6热这里在线精品视频| 99riav亚洲国产免费| 日韩欧美在线二视频 | 最新的欧美精品一区二区| 欧美中文综合在线视频| 777久久人妻少妇嫩草av网站| 美女午夜性视频免费| 亚洲欧美一区二区三区黑人| 色综合欧美亚洲国产小说| 看黄色毛片网站| 婷婷成人精品国产| 国产免费现黄频在线看| 最新在线观看一区二区三区| 天天躁日日躁夜夜躁夜夜| 黑人巨大精品欧美一区二区蜜桃| 国产有黄有色有爽视频| av网站在线播放免费| 亚洲va日本ⅴa欧美va伊人久久| 精品久久久久久,| 国产精品自产拍在线观看55亚洲 | 无遮挡黄片免费观看| 欧美日韩亚洲综合一区二区三区_| 热99re8久久精品国产| 黄色毛片三级朝国网站| 亚洲av成人不卡在线观看播放网| 国内毛片毛片毛片毛片毛片| 777米奇影视久久| av不卡在线播放| 精品久久久久久,| 国产伦人伦偷精品视频| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久免费高清国产稀缺| 免费在线观看完整版高清| 欧美日韩成人在线一区二区| 国产成人免费无遮挡视频| 三上悠亚av全集在线观看| 99国产精品免费福利视频| 久久性视频一级片| 中文字幕高清在线视频| 天天影视国产精品| 高清毛片免费观看视频网站 | 免费在线观看日本一区| 精品电影一区二区在线| 久久中文字幕一级| svipshipincom国产片| 在线观看66精品国产| 欧美老熟妇乱子伦牲交| 操出白浆在线播放| 一级毛片高清免费大全| 亚洲国产精品sss在线观看 | 国产精品美女特级片免费视频播放器 | 中文字幕另类日韩欧美亚洲嫩草| 老司机福利观看| 手机成人av网站| 欧美黑人精品巨大| 男人的好看免费观看在线视频 | 好看av亚洲va欧美ⅴa在| 无限看片的www在线观看| 日韩中文字幕欧美一区二区| 国产熟女午夜一区二区三区| 久久久精品国产亚洲av高清涩受| 精品久久久久久电影网| 国产精品免费大片| 中出人妻视频一区二区| 中文字幕最新亚洲高清| 亚洲在线自拍视频| 久久精品国产综合久久久| 女性被躁到高潮视频| 精品国产超薄肉色丝袜足j| 精品亚洲成国产av| 丰满的人妻完整版| 天天添夜夜摸| 最新的欧美精品一区二区| e午夜精品久久久久久久| av网站在线播放免费| 啦啦啦视频在线资源免费观看| 侵犯人妻中文字幕一二三四区| 99re在线观看精品视频| 色精品久久人妻99蜜桃| 国产1区2区3区精品| 精品人妻1区二区| 欧美日韩乱码在线| xxx96com| 18禁观看日本| 啪啪无遮挡十八禁网站| 亚洲 欧美一区二区三区| 免费在线观看影片大全网站| 一区二区日韩欧美中文字幕| 成在线人永久免费视频| av网站在线播放免费| 又紧又爽又黄一区二区| 欧美日本中文国产一区发布| 亚洲一区高清亚洲精品| 国产亚洲精品久久久久久毛片 | 久久人妻av系列| 超碰97精品在线观看| 色94色欧美一区二区| 国产97色在线日韩免费| 亚洲一区中文字幕在线| 99国产精品免费福利视频| x7x7x7水蜜桃| 在线观看免费视频日本深夜| 国产成人精品在线电影| 老司机福利观看| 国产精品一区二区精品视频观看| 亚洲精品粉嫩美女一区| 法律面前人人平等表现在哪些方面| 午夜两性在线视频| 精品久久久久久久毛片微露脸| 很黄的视频免费| 亚洲国产精品合色在线| 老汉色av国产亚洲站长工具| 婷婷丁香在线五月| 欧美 日韩 精品 国产| 精品国产一区二区三区久久久樱花| 91麻豆av在线| 少妇粗大呻吟视频| 69av精品久久久久久| 丰满迷人的少妇在线观看| 久久人人97超碰香蕉20202| 免费观看人在逋| 国产成人啪精品午夜网站| 国产亚洲av高清不卡| 欧美日本中文国产一区发布| 成熟少妇高潮喷水视频| 午夜精品国产一区二区电影| 香蕉丝袜av| 欧美国产精品一级二级三级| 精品人妻1区二区| 啪啪无遮挡十八禁网站| 黑人操中国人逼视频| 午夜福利视频在线观看免费| tocl精华| 国产精品国产高清国产av | 中文字幕高清在线视频| 1024香蕉在线观看| 亚洲成a人片在线一区二区| 手机成人av网站| www.999成人在线观看| 亚洲免费av在线视频| 777久久人妻少妇嫩草av网站| 亚洲精品美女久久久久99蜜臀| 日本撒尿小便嘘嘘汇集6| 下体分泌物呈黄色| 国产男女超爽视频在线观看| 欧美中文综合在线视频| 极品少妇高潮喷水抽搐| 欧美成狂野欧美在线观看| 咕卡用的链子| 欧美黑人欧美精品刺激| 91麻豆av在线| 老熟女久久久| 国产精品成人在线| 9热在线视频观看99| 老汉色av国产亚洲站长工具| 成人18禁在线播放| 亚洲美女黄片视频| 后天国语完整版免费观看| 美国免费a级毛片| av视频免费观看在线观看| 亚洲人成电影免费在线| 变态另类成人亚洲欧美熟女 | 99re6热这里在线精品视频| 十分钟在线观看高清视频www| 国产亚洲精品久久久久久毛片 | 久热爱精品视频在线9| 捣出白浆h1v1| 日韩欧美在线二视频 | 久久久久视频综合| 久久人人97超碰香蕉20202| 久久亚洲真实| 黑人操中国人逼视频| 欧美午夜高清在线| 国产单亲对白刺激| 久久人人97超碰香蕉20202| 亚洲五月天丁香| 啦啦啦视频在线资源免费观看| 亚洲片人在线观看| 日本vs欧美在线观看视频| 亚洲专区国产一区二区| 男女免费视频国产| 久久 成人 亚洲| 两个人免费观看高清视频| 亚洲人成电影免费在线| 在线观看免费视频日本深夜| 69av精品久久久久久| 99国产精品99久久久久| 99re6热这里在线精品视频| 精品午夜福利视频在线观看一区| av中文乱码字幕在线| a级毛片在线看网站| 热re99久久国产66热| 久久久精品免费免费高清| 18禁黄网站禁片午夜丰满| 在线观看66精品国产| 大陆偷拍与自拍| 成人18禁高潮啪啪吃奶动态图| 国产亚洲精品第一综合不卡| 国产欧美日韩综合在线一区二区| 天堂俺去俺来也www色官网| 免费女性裸体啪啪无遮挡网站| 9色porny在线观看| 午夜福利乱码中文字幕| 老司机午夜福利在线观看视频| 欧美日韩瑟瑟在线播放| 久久久水蜜桃国产精品网| 午夜成年电影在线免费观看| 免费黄频网站在线观看国产| 亚洲人成伊人成综合网2020| 999久久久国产精品视频| 久久精品成人免费网站| 色婷婷久久久亚洲欧美| 国产男女超爽视频在线观看| 91精品国产国语对白视频| 黄色女人牲交| 777米奇影视久久| 亚洲精品国产区一区二| 在线观看午夜福利视频| 97人妻天天添夜夜摸| 国产成人精品在线电影| 免费观看精品视频网站| 啦啦啦免费观看视频1| 欧美成狂野欧美在线观看| 色婷婷久久久亚洲欧美| 精品国产一区二区三区久久久樱花| 最近最新中文字幕大全电影3 | 在线天堂中文资源库| 精品欧美一区二区三区在线| 悠悠久久av| 亚洲黑人精品在线| 女人爽到高潮嗷嗷叫在线视频| 国产精品免费一区二区三区在线 | 黄网站色视频无遮挡免费观看| 成在线人永久免费视频| 夜夜躁狠狠躁天天躁| 色在线成人网| 国产极品粉嫩免费观看在线| 麻豆乱淫一区二区| 亚洲aⅴ乱码一区二区在线播放 | 王馨瑶露胸无遮挡在线观看| 精品国产美女av久久久久小说| 很黄的视频免费| 动漫黄色视频在线观看| 中国美女看黄片| 欧美乱色亚洲激情| 亚洲精品美女久久久久99蜜臀| 亚洲avbb在线观看| 国产成人系列免费观看| 日韩有码中文字幕| 婷婷成人精品国产| 最近最新免费中文字幕在线| 成人永久免费在线观看视频| 国产免费现黄频在线看| 国产成人免费无遮挡视频| 18禁国产床啪视频网站| 国产高清videossex| 无遮挡黄片免费观看| 欧美日韩一级在线毛片| 手机成人av网站| 亚洲性夜色夜夜综合| 久久热在线av| 他把我摸到了高潮在线观看| 欧美中文综合在线视频| 国产亚洲精品第一综合不卡| 精品久久久久久电影网| 纯流量卡能插随身wifi吗| 18在线观看网站| 欧美激情久久久久久爽电影 | 国产一区二区激情短视频| 国产日韩一区二区三区精品不卡| 18禁国产床啪视频网站| 最近最新中文字幕大全电影3 | 婷婷丁香在线五月| 国产成人系列免费观看| 老司机午夜福利在线观看视频| 99国产极品粉嫩在线观看| 后天国语完整版免费观看| 午夜两性在线视频| 桃红色精品国产亚洲av| 天堂俺去俺来也www色官网| 免费人成视频x8x8入口观看| 亚洲一区中文字幕在线| 亚洲av欧美aⅴ国产| 又黄又粗又硬又大视频| 日本a在线网址| 咕卡用的链子| 中文字幕人妻丝袜制服| av电影中文网址| 免费高清在线观看日韩| 天天躁日日躁夜夜躁夜夜| 久久久国产精品麻豆| 怎么达到女性高潮| 久久久水蜜桃国产精品网| 我的亚洲天堂| 捣出白浆h1v1| 丁香六月欧美| 露出奶头的视频| 一级毛片精品| 亚洲av成人av| 50天的宝宝边吃奶边哭怎么回事| 99国产精品免费福利视频| 精品第一国产精品| 欧美黑人欧美精品刺激| 亚洲 欧美一区二区三区| 亚洲成av片中文字幕在线观看| 正在播放国产对白刺激| 亚洲,欧美精品.| av一本久久久久| 成年动漫av网址| 校园春色视频在线观看| 国产淫语在线视频| 欧美中文综合在线视频| 国产欧美日韩一区二区三| 日韩欧美在线二视频 | 丁香六月欧美| 91成年电影在线观看| 国产激情欧美一区二区| 视频区欧美日本亚洲| 日日夜夜操网爽| 久久香蕉国产精品| 丝袜美足系列| 久久国产精品大桥未久av| 丝袜美足系列| 免费观看人在逋| 国产精品久久久久成人av| 91av网站免费观看| 99国产精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 老司机亚洲免费影院| 狠狠婷婷综合久久久久久88av| 午夜福利影视在线免费观看| 男女之事视频高清在线观看| 99国产精品一区二区三区| 男女之事视频高清在线观看| 欧美不卡视频在线免费观看 | 久久草成人影院| 免费黄频网站在线观看国产| 中文字幕人妻熟女乱码| 一区二区三区国产精品乱码| 精品第一国产精品| 国产日韩欧美亚洲二区| 国产精品 国内视频| 免费在线观看日本一区| 欧美成人午夜精品| 国产91精品成人一区二区三区| 亚洲欧美激情综合另类| 麻豆乱淫一区二区| 久久中文字幕一级| 亚洲精品久久午夜乱码| 精品国产一区二区三区四区第35| 69精品国产乱码久久久| 久热爱精品视频在线9| 99精国产麻豆久久婷婷| 久久影院123| 在线视频色国产色| 69精品国产乱码久久久| 九色亚洲精品在线播放| 人人妻人人澡人人爽人人夜夜| 久久ye,这里只有精品| 色在线成人网| 久9热在线精品视频| 精品福利观看| 亚洲九九香蕉| 久久 成人 亚洲| 1024香蕉在线观看| 国产亚洲欧美精品永久| ponron亚洲| 老鸭窝网址在线观看| 国产精品98久久久久久宅男小说| 天天躁日日躁夜夜躁夜夜| 欧美性长视频在线观看| 他把我摸到了高潮在线观看| 91麻豆av在线| 日韩免费av在线播放| 桃红色精品国产亚洲av| 少妇 在线观看| 国产男女超爽视频在线观看| 熟女少妇亚洲综合色aaa.| 悠悠久久av| 久久久久久久久久久久大奶| 久久精品国产综合久久久| 99精品欧美一区二区三区四区| 少妇的丰满在线观看| 国产成人欧美| 在线观看免费日韩欧美大片| 男人舔女人的私密视频| 国产精品亚洲av一区麻豆| 一级a爱片免费观看的视频| 午夜成年电影在线免费观看| 99精品欧美一区二区三区四区| 国精品久久久久久国模美| 90打野战视频偷拍视频| 欧美日韩瑟瑟在线播放| 大香蕉久久成人网| 大型黄色视频在线免费观看| 可以免费在线观看a视频的电影网站| 在线视频色国产色| 国产高清videossex| 美女高潮喷水抽搐中文字幕| 国产1区2区3区精品| 欧美日韩亚洲综合一区二区三区_| 亚洲国产精品合色在线| 男人的好看免费观看在线视频 | 老司机靠b影院| 国产欧美日韩综合在线一区二区| 亚洲国产毛片av蜜桃av| 黑人猛操日本美女一级片| 国产成人免费观看mmmm| 久久久水蜜桃国产精品网| 成人精品一区二区免费| 亚洲av日韩在线播放| av免费在线观看网站| 制服诱惑二区| 韩国精品一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 大陆偷拍与自拍| 亚洲国产精品合色在线| 亚洲五月婷婷丁香| av线在线观看网站| 精品午夜福利视频在线观看一区| 精品乱码久久久久久99久播| 在线国产一区二区在线| 男女下面插进去视频免费观看| 国产亚洲av高清不卡| 精品熟女少妇八av免费久了| 无遮挡黄片免费观看| 欧美一级毛片孕妇| 亚洲成人国产一区在线观看| 国产成人精品久久二区二区免费| 天天操日日干夜夜撸| 精品国产超薄肉色丝袜足j| 国产精品秋霞免费鲁丝片| 国产xxxxx性猛交| 最近最新中文字幕大全电影3 | 黄色a级毛片大全视频| 国产成人精品在线电影| 午夜激情av网站| 欧美色视频一区免费| 国产亚洲精品第一综合不卡| 正在播放国产对白刺激| 午夜福利免费观看在线| 色老头精品视频在线观看| 国产精品成人在线| 精品亚洲成a人片在线观看| 法律面前人人平等表现在哪些方面| 亚洲久久久国产精品| 9191精品国产免费久久| 日韩人妻精品一区2区三区| 成人av一区二区三区在线看| 成年女人毛片免费观看观看9 | 欧美另类亚洲清纯唯美| 久久久久精品人妻al黑| 欧美 亚洲 国产 日韩一| 麻豆av在线久日| 亚洲全国av大片| 18禁裸乳无遮挡动漫免费视频| 久久久精品免费免费高清| 脱女人内裤的视频| 免费女性裸体啪啪无遮挡网站| 国产精品一区二区精品视频观看| 别揉我奶头~嗯~啊~动态视频| av天堂在线播放| 国产精品永久免费网站| 亚洲精品国产精品久久久不卡| 国产精品久久久久成人av| 国产黄色免费在线视频| av线在线观看网站| 在线播放国产精品三级| 久久久久久久午夜电影 | 日韩制服丝袜自拍偷拍| 亚洲avbb在线观看| 久久香蕉国产精品| 老司机在亚洲福利影院| 一级作爱视频免费观看| 久久狼人影院| 黄色毛片三级朝国网站| 久久亚洲真实| 自线自在国产av| 午夜成年电影在线免费观看| 国产精华一区二区三区| 午夜免费成人在线视频| 国产精品一区二区在线不卡| 极品教师在线免费播放| 天堂俺去俺来也www色官网| e午夜精品久久久久久久| 久久亚洲精品不卡| 伦理电影免费视频| 精品人妻熟女毛片av久久网站| 午夜福利影视在线免费观看| 国产精品成人在线| 亚洲国产毛片av蜜桃av| 高潮久久久久久久久久久不卡| 国产伦人伦偷精品视频| 超碰97精品在线观看| 中文字幕色久视频| 亚洲成av片中文字幕在线观看| 国产精品久久久久成人av| 两性夫妻黄色片| 亚洲精品在线观看二区| 国内久久婷婷六月综合欲色啪| av有码第一页| 久久久久国产精品人妻aⅴ院 | 在线观看一区二区三区激情| 我的亚洲天堂| 国产成人精品无人区| 久久久久久人人人人人| 校园春色视频在线观看| 欧美激情极品国产一区二区三区| 欧美日韩av久久| 18禁黄网站禁片午夜丰满| 国产av一区二区精品久久| 亚洲中文日韩欧美视频| 久久精品91无色码中文字幕| 婷婷精品国产亚洲av在线 | 老司机影院毛片| 亚洲欧美激情在线| 国产单亲对白刺激| 99在线人妻在线中文字幕 | 99精品久久久久人妻精品| 亚洲国产看品久久| 日日摸夜夜添夜夜添小说| 一级片免费观看大全| 国产在线一区二区三区精| 欧美大码av| 国产有黄有色有爽视频| 欧美亚洲日本最大视频资源| 色综合欧美亚洲国产小说| 久久香蕉国产精品| 又紧又爽又黄一区二区| 首页视频小说图片口味搜索| 在线免费观看的www视频| 国精品久久久久久国模美| av超薄肉色丝袜交足视频| 精品人妻1区二区| 91老司机精品| 亚洲国产精品一区二区三区在线| 每晚都被弄得嗷嗷叫到高潮| 亚洲avbb在线观看| 国产一区二区三区在线臀色熟女 | 捣出白浆h1v1| 国产精品99久久99久久久不卡| 亚洲欧美激情综合另类| 深夜精品福利| 9热在线视频观看99| 老司机午夜十八禁免费视频| 深夜精品福利| 下体分泌物呈黄色| 日本撒尿小便嘘嘘汇集6| 99精品欧美一区二区三区四区| 亚洲一区中文字幕在线| 99久久国产精品久久久| 欧美激情 高清一区二区三区| 国产精品乱码一区二三区的特点 | 精品亚洲成a人片在线观看| 一二三四社区在线视频社区8| 怎么达到女性高潮| 国产精品av久久久久免费| 侵犯人妻中文字幕一二三四区| 18禁裸乳无遮挡免费网站照片 | 国产高清视频在线播放一区| 18禁裸乳无遮挡动漫免费视频| 热99久久久久精品小说推荐| 91精品国产国语对白视频| 99国产精品一区二区蜜桃av | 视频区图区小说| 国产不卡一卡二| 亚洲欧美日韩高清在线视频| 午夜亚洲福利在线播放| 成人特级黄色片久久久久久久| 日本撒尿小便嘘嘘汇集6| 在线永久观看黄色视频| 免费一级毛片在线播放高清视频 | videos熟女内射| 在线播放国产精品三级| 视频区欧美日本亚洲| a级毛片黄视频| 欧美日韩视频精品一区| 天天躁狠狠躁夜夜躁狠狠躁| 久久ye,这里只有精品|