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

    仿生梯度圓環(huán)防護(hù)系統(tǒng)的耐撞性設(shè)計

    2022-08-01 07:31:00邢運張橋楊先鋒劉華楊嘉陵
    航空學(xué)報 2022年6期
    關(guān)鍵詞:抗沖擊圓環(huán)梯度

    邢運,張橋,楊先鋒,*,劉華,楊嘉陵

    1. 北京航空航天大學(xué) 航空科學(xué)與工程學(xué)院,北京 100083 2. 北京航空航天大學(xué) 先進(jìn)結(jié)構(gòu)沖擊與仿生力學(xué)實驗室,北京 100083

    隨著飛行器結(jié)構(gòu)在強(qiáng)動載荷下的抗沖擊力學(xué)性能要求越來越高,相應(yīng)的結(jié)構(gòu)與人體抗沖擊安全防護(hù)問題也受到了廣泛的關(guān)注,成為探究學(xué)術(shù)前沿和突破關(guān)鍵技術(shù)過程中亟待解決的核心課題。針對飛行器結(jié)構(gòu)的抗沖擊防護(hù)問題,國內(nèi)外學(xué)者已經(jīng)開展了大量的研究工作,設(shè)計開發(fā)了一系列傳統(tǒng)的抗沖擊防護(hù)結(jié)構(gòu),主要包括由金屬或復(fù)合材料制成的薄壁構(gòu)件(如波紋梁、方管、多邊形管、圓管、錐管等)、多胞材料(如蜂窩、泡沫、點陣材料)、填充薄壁管件、三明治夾芯結(jié)構(gòu)等。傳統(tǒng)結(jié)構(gòu)的抗沖擊設(shè)計涉及多種復(fù)雜的非線性力學(xué)現(xiàn)象,除了彈塑性動力大變形外,還伴有材料疲勞損傷、斷裂、結(jié)構(gòu)屈曲和坍塌,以及摩擦和接觸等復(fù)雜物理因素,其中還有許多與之相關(guān)的固體力學(xué)和沖擊動力學(xué)重大基礎(chǔ)問題沒有解決。此外,目前的傳統(tǒng)抗沖擊防護(hù)結(jié)構(gòu)并沒有全部參與到能量的吸收與耗散過程中,不能滿足較寬廣的適用環(huán)境,因此防護(hù)效率在一定程度上受到了限制,抗沖擊性能很難達(dá)到最優(yōu)化,同時也造成了很大資源浪費。

    大自然為設(shè)計和制造具有優(yōu)異力學(xué)性能的先進(jìn)材料和結(jié)構(gòu)提供了豐富的仿生靈感。通過對天然生物抗沖擊策略的研究與探索,對生物身體防護(hù)機(jī)制的模仿與學(xué)習(xí),國內(nèi)外學(xué)者在仿生抗沖擊防護(hù)研究領(lǐng)域開展了大量工作,研制出多種多樣的新型抗沖擊仿生材料與結(jié)構(gòu)。這些仿生結(jié)構(gòu)具有出色的防護(hù)性能,例如有限的初始峰值力、足夠長的工作行程、較高的能量吸收能力,因為它們被設(shè)計為采用塑性變形或破碎的形式將大部分沖擊動能轉(zhuǎn)化為非彈性能耗散。但是,采用塑性變形或破碎形式能量耗散策略的能量吸收結(jié)構(gòu)只能使用一次。當(dāng)受到外部沖擊載荷頻繁作用時,以上防護(hù)結(jié)構(gòu)無法滿足提高可重復(fù)使用耐沖擊性的特定要求。盡管某些可重用防護(hù)結(jié)構(gòu)的抗沖擊性能目前無法與一次性防護(hù)結(jié)構(gòu)相提并論,但它們還具有許多其他優(yōu)異的性能,例如可重構(gòu)性、可編程性和廣泛的適應(yīng)性。特別是當(dāng)抗沖擊結(jié)構(gòu)需要承受極端和復(fù)雜的外部載荷時,結(jié)構(gòu)的可重復(fù)使用性、可編程性以及形狀可重構(gòu)性在防護(hù)系統(tǒng)中就顯得非常重要。

    為了解決以上問題,本文嘗試從自然進(jìn)化的動物抗沖擊策略中尋求有效的解決方案。在自然界中,動物進(jìn)化出的大多數(shù)身體防護(hù)機(jī)制都是可重用且高效的,動物利用其生物組織的彈性或黏彈性來抵抗沖擊載荷,這些機(jī)制很少引起生物組織的塑性變形,從而防止了沖擊載荷對身體的不可逆轉(zhuǎn)傷害,這對研究沖擊防護(hù)具有重要的生物啟示。近期,Xing和Yang在天然昆蟲角質(zhì)層的空間層次結(jié)構(gòu)中發(fā)現(xiàn)了一種抗沖擊策略。這一發(fā)現(xiàn)表明,昆蟲甲殼角質(zhì)層沿厚度方向具有不連續(xù)的指數(shù)剛度梯度(DC-EXP),可使結(jié)構(gòu)內(nèi)部應(yīng)力達(dá)到最小,從而達(dá)到最佳的抗沖擊性能和防御效果。由此推理,剛度分布對人造防護(hù)結(jié)構(gòu)的抗沖擊性能也應(yīng)產(chǎn)生重大影響。但是,前人在設(shè)計可重復(fù)使用防護(hù)結(jié)構(gòu)時尚未考慮剛度分布對其抗沖擊能力和防護(hù)性能的確切影響。比如,不同的剛度梯度在沖擊載荷下對防御結(jié)構(gòu)的防護(hù)性能的作用、不同的剛度梯度設(shè)計方案對防御結(jié)構(gòu)的耐沖擊性能的提高程度、最適合結(jié)構(gòu)防護(hù)的剛度梯度分布等問題仍然不清楚。

    為了解決上述問題并設(shè)計出易于制造的防護(hù)結(jié)構(gòu),本文受昆蟲角質(zhì)層剛度分布的生物啟發(fā),提出一種新型可重復(fù)使用彈性梯度圓環(huán)陣列結(jié)構(gòu),旨在通過調(diào)節(jié)圓環(huán)間彈性模量、半徑和厚度的分布,改變圓環(huán)陣列的相對剛度梯度,從而減小與內(nèi)部結(jié)構(gòu)的碰撞力,提高圓環(huán)的沖擊防護(hù)特性。

    1 仿生剛度梯度圓環(huán)系統(tǒng)

    1.1 幾何構(gòu)造和有限元模型

    仿生設(shè)計思路如圖1所示,使用一系列彈性環(huán)陣列結(jié)構(gòu)代替昆蟲角質(zhì)層,類似于集成不同剛度梯度分布的空間分層結(jié)構(gòu),對內(nèi)部結(jié)構(gòu)起到抗沖擊防護(hù)作用。

    圖1 仿生圓環(huán)防護(hù)系統(tǒng)設(shè)計示意圖Fig.1 Schematic diagram of bionic ring protection system design

    根據(jù)Timoshenko提出的彈性解,彈性環(huán)的等效剛度的計算公式為

    (1)

    式中:為圓環(huán)基本材料的彈性模量;為圓環(huán)的半徑;為圓環(huán)的厚度;為圓環(huán)的寬度。根據(jù)式(1), 可以通過調(diào)整3個參數(shù)(彈性模量、半徑、厚度)對彈性環(huán)系統(tǒng)的剛度梯度編程,從而獲得最佳的抗沖擊性能設(shè)計方案(見圖2)。

    圖2 仿生剛度梯度圓環(huán)防護(hù)系統(tǒng)Fig.2 Bionic stiffness gradient ring protection system

    為研究仿生彈性梯度圓環(huán)系統(tǒng)的防護(hù)性能,使用HyperMesh軟件構(gòu)建了具有不同等效剛度梯度分布的單列圓環(huán)系統(tǒng)有限元模型,模型包括3個部分:剛性質(zhì)量塊、彈性環(huán)系統(tǒng)、剛性底座(見圖3),并用非線性顯式有限元算法工具ANSYS/LS-DYNA來計算此仿生結(jié)構(gòu)在沖擊載荷下的動力響應(yīng)。約束剛性質(zhì)量塊其他方向的自由度,使其只能沿方向(沖擊方向,見圖3)進(jìn)行平移運動,并以的初速度垂直向下撞擊彈性圓環(huán)系統(tǒng)。將剛性底座最下層單元節(jié)點固支約束。彈性圓環(huán)系統(tǒng)采用線彈性材料模型,可以分別調(diào)整單個圓環(huán)的彈性模量、半徑、厚度來設(shè)計圓環(huán)系統(tǒng)的等效剛度梯度分布。有限元模型采用8節(jié)點實體單元劃分網(wǎng)格,在圓環(huán)厚度方向上至少劃分4層,這可以提供良好的計算精度并同時節(jié)省計算成本(見圖3)。為避免彈性環(huán)之間發(fā)生滲透或相交現(xiàn)象,采用自動單面接觸來模擬圓環(huán)系統(tǒng)自身的接觸行為;采用自動面-面接觸來模擬撞擊質(zhì)量塊、剛性底座與圓環(huán)系統(tǒng)之間的接觸行為,各接觸面之間的靜動摩擦系數(shù)均設(shè)置為0.2。

    圖3 仿生圓環(huán)防護(hù)系統(tǒng)有限元模擬示意圖Fig.3 Schematic diagram of finite element simulation of bionic ring protection system

    為模仿甲蟲角質(zhì)層防止應(yīng)力波傳播及其防護(hù)性能的功能特點,通過參考之前研究昆蟲角質(zhì)層剛度分布的方法,沿沖擊方向同樣提出7種等效剛度梯度(見圖4),分別為指數(shù)函數(shù)分布(EXP)、線性函數(shù)分布(LIN)、對數(shù)函數(shù)分布(LG)、三角函數(shù)分布(SIN-Ⅰ、SIN-Ⅱ)、二次函數(shù)分布(SQ-Ⅰ、SQ-Ⅱ),并且對具有這7種不同剛度梯度分布的仿生圓環(huán)系統(tǒng)的防護(hù)性能進(jìn)行比較。以上7種等效剛度梯度分布曲線具有凹形或凸形輪廓,本文中沿沖擊方向(方向)斜率增大的剛度曲線定義為凹曲線,斜率減小的剛度曲線定義為凸曲線(見圖4)。因此,凹形剛度梯度有EXP、SIN-Ⅰ、SQ-Ⅰ,凸型剛度梯度有LG、SQ-Ⅱ、SQ-Ⅱ。

    圖4 沿梯度彈性圓環(huán)受沖擊方向(x方向)不同等效剛度分布形式Fig.4 Distribution forms of different equivalent stiffness along impact direction (x direction) of gradient elastic ring

    Xing和Yang的研究表明,甲蟲中表皮層對彈性應(yīng)力波的衰減作用以及對整個角質(zhì)層的抗沖擊性能影響最大。中表皮最外層的最大剛度約為其最內(nèi)層的8倍。本文以天然甲蟲角質(zhì)層的中表皮層為生物靈感,以半徑為11.6 mm、厚度為3. 2 mm、寬度為12.5 mm的彈性鋼環(huán)作為參考基本環(huán),發(fā)展出圓環(huán)陣列防護(hù)系統(tǒng),從而模仿中表皮層的抗沖擊特性。如果參考基本鋼環(huán)的彈性模量為210 GPa,其等效剛度約為48.179 kN/mm。根據(jù)參考基本鋼環(huán),分別通過調(diào)整環(huán)的彈性模量、半徑、厚度,可以得到圖4中所示的具有7種等效剛度梯度分布的仿生圓環(huán)防護(hù)系統(tǒng)。

    1.2 理論分析模型

    基于多胞結(jié)構(gòu)受撞擊的理論分析模型,提出一種改進(jìn)的一維質(zhì)量-彈簧模型,研究1.1節(jié)提出的仿生剛度梯度圓環(huán)防護(hù)系統(tǒng)在沖擊載荷下的動態(tài)響應(yīng)。為便于理解,首先考慮了單個環(huán)在質(zhì)量塊的撞擊作用下的改進(jìn)質(zhì)量-彈簧模型,如圖5(a)所示。如圖5(b)所示,單個環(huán)可以簡化為1個集中質(zhì)量()和1個無質(zhì)量的代表單個圓環(huán)壓縮特性的彈簧(),可將圓環(huán)壓縮特性近似看作線彈性。

    質(zhì)量塊與單環(huán)上邊緣之間的接觸由接觸彈簧(圖5(b)所示的彈簧)表示。在前人研究中,對于理想彈塑性的多胞結(jié)構(gòu)受撞擊問題,常常將與受壓材料的屈服應(yīng)力相關(guān)的恒力彈簧表示為接觸彈簧,它可以準(zhǔn)確地預(yù)測質(zhì)量塊與圓環(huán)碰撞過程中的能量耗散。然而,在該研究中,恒力彈簧不再適用于撞擊質(zhì)量塊和線彈性圓環(huán)之間的碰撞,因為它們之間接觸力的改變不能被忽略。

    圖5 仿生圓環(huán)防護(hù)系統(tǒng)理論模型示意圖Fig.5 Schematic diagram of theoretical model of bionic ring protection system

    基于接觸力學(xué),可以引入空心率確定撞擊質(zhì)量塊與線彈性薄壁圓環(huán)之間的碰撞接觸力

    =(1-)

    (2)

    式中:為2個實體之間的接觸力;為空心率,計算公式為

    (3)

    對于撞擊質(zhì)量塊和線彈性薄壁圓環(huán)的接觸問題,非赫茲接觸力學(xué)問題的解可提供圓環(huán)的徑向變形(接觸彈簧的變形),計算公式為

    (4)

    (5)

    式中:為等效模量,計算公式為

    (6)

    式中:為撞擊質(zhì)量塊的彈性模量;為圓環(huán)的彈性模量;為撞擊質(zhì)量塊的泊松比;為圓環(huán)的泊松比。

    對于表示薄壁圓環(huán)壓縮力學(xué)行為的線彈性彈簧(),Timoshenko推導(dǎo)出了彈性圓環(huán)的等效剛度,如式(1)所示。則薄壁彈性圓環(huán)的彈簧力可以表示為

    =

    (7)

    式中:為圓環(huán)徑向減小的距離。

    對于圖5(c)所示的單列梯度彈性圓環(huán)系統(tǒng),幾何可重復(fù)性可以將整個結(jié)構(gòu)簡化成一系列由1個無質(zhì)量線彈性彈簧連接1個集中質(zhì)量塊凝聚的質(zhì)量-彈簧模型(見圖5(d))。假定單列梯度彈性圓環(huán)沿方向受到?jīng)_擊,為在撞擊瞬間圓環(huán)系統(tǒng)第個集中質(zhì)量塊的位置坐標(biāo)(距受沖擊端的距離),為第個集中質(zhì)量塊距其初始位置的位移,為第個彈簧的壓縮距離,計算公式為

    =-+1

    (8)

    通過對以上質(zhì)量-彈簧模型的分析,給出了梯度彈性環(huán)系統(tǒng)在沖擊載荷作用下的動態(tài)響應(yīng)控制方程

    (9)

    (10)

    (11)

    式中:為撞擊質(zhì)量塊的質(zhì)量;為撞擊質(zhì)量塊與圓環(huán)系統(tǒng)中第1個圓環(huán)碰撞時的接觸力;為第個彈性環(huán)的質(zhì)量;為第個彈性環(huán)的等效剛度,可以通過將第個彈性環(huán)的彈性模量、半徑、厚度、寬度代入式(1)計算得到。

    問題的初始條件為

    (12)

    (13)

    式中:為撞擊質(zhì)量塊的初始沖擊速度。模型中彈簧不能承受拉力,因此,如果彈簧的壓縮距離<0,此時在彈簧中沒有任何力存在,彈簧力等于零(=0),直到它受到進(jìn)一步壓縮。本文理論模型通過MATLAB編程求解,采用四階Runge-Kutta算法計算預(yù)測撞擊質(zhì)量塊和第1個圓環(huán)之間的碰撞接觸力以及圓環(huán)系統(tǒng)末端與剛性壁之間的相互作用力。

    1.3 有限元模型驗證

    為了檢查1.1節(jié)中有限元分析的有效性,基于先前的實驗對仿生剛度梯度圓環(huán)系統(tǒng)受沖擊載荷問題進(jìn)行了比較研究。在文獻(xiàn)[31]中,實驗測試了半徑為19.05 mm、厚度為1 mm、寬度為10 mm 的彈性黃銅環(huán)受沖擊載荷時的動態(tài)響應(yīng),如圖6所示,銅環(huán)材料的密度和彈性模量分別為8 199 kg/m、120.4 GPa。建立了對應(yīng)于實驗的理論模型和有限元模型,計算了撞擊質(zhì)量塊和銅環(huán)之間的接觸力(見圖6(a))以及金屬環(huán)和剛性壁之間的碰撞力(見圖6(b)),可以發(fā)現(xiàn)有限元數(shù)值結(jié)果與實驗測量及理論預(yù)測結(jié)果分別吻合良好。

    圖6 單個金屬圓環(huán)受撞擊時的典型碰撞接觸力和反作用力Fig.6 Typical collision contact force and reaction force of single metal ring being impacted

    此外,將圖5(d)中的九環(huán)質(zhì)量-彈簧理論模型對1、10、20 m/s 3種初始碰撞速度的沖擊預(yù)測結(jié)果與有限元數(shù)值模擬結(jié)果進(jìn)行對比,結(jié)果如圖7所示。從圖7中可以看出,在低速沖擊情況下,分別通過理論模型預(yù)測和有限元模擬獲得的接觸力和碰撞反作用力彼此保持一致。以上驗證表明有限元模型構(gòu)建的精確性是完全可接受的,可作為可靠手段對仿生梯度圓環(huán)防護(hù)系統(tǒng)進(jìn)、一步進(jìn)行參數(shù)分析。

    圖7 圓環(huán)系統(tǒng)受初速度1 m/s的質(zhì)量塊撞擊時的數(shù)值模擬與理論預(yù)測結(jié)果比較Fig.7 Comparison between numerical simulation and theoretical prediction results of ring system impacted by 1 m/s initial velocity mass

    2 單列仿生剛度梯度圓環(huán)的參數(shù)化分析

    第1節(jié)通過理論模型和實驗對比驗證了有限元分析的準(zhǔn)確性,本節(jié)采用有限元模型來研究具有不同等效剛度梯度分布的單列仿生圓環(huán)系統(tǒng)的抗沖擊能力;同時評估剛度曲線輪廓的凹凸形對防護(hù)性能的影響。評判仿生圓環(huán)系統(tǒng)防護(hù)性能的重要參考指標(biāo)之一就是沖擊遠(yuǎn)端圓環(huán)和內(nèi)部基底(剛性壁)之間的碰撞力,此碰撞力減小有助于降低內(nèi)部被保護(hù)結(jié)構(gòu)的損傷風(fēng)險,這意味著仿生剛度梯度圓環(huán)系統(tǒng)具有良好的抗沖擊能力。

    根據(jù)式(1),圓環(huán)系統(tǒng)等效剛度梯度的分布取決于基體材料的彈性模量、圓環(huán)半徑、圓環(huán)壁厚、寬度。本文圓環(huán)的寬度保持不變(=12.5 mm), 等效剛度分布僅由基體材料(彈性模量)和圓環(huán)半徑與壁厚之比(/)決定。當(dāng)圓環(huán)材料相同時(彈性模量保持不變),影響等效剛度的核心因素是徑厚比(/),徑厚比越大,剛度越?。粡胶癖仍叫?,剛度越大。通過控制徑厚比(/)可以得到相同的等效剛度分布形式,但這并不意味著在下面的情況下,比值/可以被視為一個獨立的參數(shù)分析。因為半徑或厚度的單獨變化將導(dǎo)致2種不同的質(zhì)量(或慣性)分布,并最終導(dǎo)致圓環(huán)系統(tǒng)的抗沖擊響應(yīng)和防護(hù)特性不同。因此,需要討論在獲得相同的等效剛度分布時,彈性模量分布、半徑分布和厚度分布對圓環(huán)系統(tǒng)防護(hù)性能的影響。

    2.1 彈性模量分布對防護(hù)特性的影響

    根據(jù)式(1),保持其他參數(shù)不變,調(diào)整彈性模量以滿足圖4所示的7種不同剛度梯度分布形式。選擇半徑為11.6 mm、厚度為3.2 mm、寬度為12.5 mm的不銹鋼環(huán)作為參考圓環(huán)。參考環(huán)材料的彈性模量為=210 GPa、泊松比=03、密度=7 800 kg/m,進(jìn)一步可計算出參考環(huán)的等效剛度為48.179 kN/mm。圖8(a)給出了調(diào)整彈性模量分布達(dá)到不同等效剛度分布的圓環(huán)系統(tǒng)在,圖8(b)、圖8(c)則給出了此梯度圓環(huán)系統(tǒng)受到質(zhì)量為0.2 kg、初速度為1 m/s的質(zhì)量塊撞擊作用時,沖擊遠(yuǎn)端圓環(huán)與剛性壁之間的碰撞力曲線。

    圖8 通過控制彈性模量設(shè)計出的具有不同等效剛度梯度的圓環(huán)系統(tǒng)受沖擊時的反作用力Fig.8 Reaction forces of ring system with different equivalent stiffness gradients under impact designed by controlling elastic modulus

    結(jié)果表明,所有碰撞力-時間曲線只有一個波峰。為了便于描述,凹形剛度梯度(EXP、SIN-Ⅰ、SQ-Ⅰ)定義為Ⅰ型梯度,凸型剛度梯度(LG、SQ-Ⅱ、SQ-Ⅱ)和LIN定義為Ⅱ型梯度。研究發(fā)現(xiàn),Ⅰ型梯度的峰值力小于Ⅱ型梯度的峰值力。與Ⅱ型梯度中最大的峰值力(由LG得到)相比,Ⅰ型梯度中最小的峰值力(由EXP得到)降低了13.6%。此外,與其他類型相比,Ⅰ型梯度中的峰值力有所延遲。以上結(jié)果說明,凹形剛度梯度會減緩沖擊力的傳播,并對內(nèi)部結(jié)構(gòu)具有更好的防護(hù)性能。造成這種現(xiàn)象的一個重要原因是凹型曲線(Ⅰ型梯度)可以更有效地衰減彈性應(yīng)力波的傳播并同時在復(fù)合結(jié)構(gòu)中保持更低的應(yīng)力狀態(tài)。

    在以上模型中,沿仿生梯度圓環(huán)系統(tǒng)受沖擊方向(方向)具有指數(shù)剛度梯度(EXP)的模型對內(nèi)部結(jié)構(gòu)產(chǎn)生的碰撞力最小(圖8(c))。這是因為EXP剛度梯度在衰減彈性應(yīng)力波傳播方面具有優(yōu)勢,在撞擊過程中可有效減小作用在內(nèi)部結(jié)構(gòu)上的碰撞力。此外,還有一個重要事實可以解釋以上發(fā)現(xiàn),相比于其他梯度分布形式,EXP剛度梯度可以通過將剛性、柔順性2個相互排斥的特性很好地結(jié)合起來,在材料的承載力和彈性之間達(dá)到更好的折衷策略。

    2.2 半徑分布對防護(hù)特性的影響

    保持其他參數(shù)不變,通過調(diào)整圓環(huán)半徑獲得具有圖4中所示的不同剛度梯度分布的圓環(huán)系統(tǒng),在相同沖擊條件下的防護(hù)性能。計算發(fā)現(xiàn),可得到與2.1節(jié)類似的結(jié)論。

    圖9給出了具有不同圓環(huán)半徑分布形式的梯度圓環(huán)系統(tǒng)在沖擊載荷作用下,圓環(huán)系統(tǒng)遠(yuǎn)端與內(nèi)部剛性壁之間的碰撞力曲線。具有Ⅰ型梯度(EXP、SIN-Ⅰ、SQ-Ⅰ)的圓環(huán)系統(tǒng)產(chǎn)生的碰撞力曲線只有一個波峰,而具有Ⅱ型梯度(LG、SQ-Ⅱ、SQ-Ⅱ、LIN)的圓環(huán)系統(tǒng)產(chǎn)生的碰撞力曲線有2個波峰(圖9)。Ⅱ型梯度產(chǎn)生的第1個碰撞力波峰出現(xiàn)在Ⅰ型梯度的波峰之前,第2個波峰出現(xiàn)在其他波峰之后。研究還可以發(fā)現(xiàn),Ⅱ型梯度產(chǎn)生的第2個波峰值小于第1個波峰值。這是因為此時的圓環(huán)系統(tǒng)質(zhì)量分布不再像圖8所示的那樣均勻。圓環(huán)的質(zhì)量會隨著半徑的增加而增加,因此,圓環(huán)系統(tǒng)遠(yuǎn)端的慣性也隨之增加。在相同的碰撞條件下,圓環(huán)的壓縮深度隨著慣性的增加而減小,從而導(dǎo)致較小的碰撞反力。與Ⅱ型梯度相比(LG、SQ-Ⅱ、SQ-Ⅱ、LIN),I型梯度(EXP、SIN-Ⅰ、SQ-Ⅰ)仍在遠(yuǎn)端圓環(huán)與內(nèi)部結(jié)構(gòu)之間的碰撞力曲線中只產(chǎn)生一個峰值(見圖9(b)),這不僅延遲了內(nèi)部結(jié)構(gòu)受到碰撞的時間,而且同時減少了沖擊初期防護(hù)系統(tǒng)與被保護(hù)結(jié)構(gòu)之間的碰撞次數(shù)。與第2.1節(jié)結(jié)論相似,與具有2個I型剛度梯度(SIN-Ⅰ、SQ-Ⅰ)的圓環(huán)系統(tǒng)相比,具有指數(shù)剛度梯度(EXP)的圓環(huán)系統(tǒng)對內(nèi)部結(jié)構(gòu)產(chǎn)生更小的碰撞力峰值,相比SIN-Ⅰ梯度下降了2.32%,相比SQ-Ⅰ梯度下降了1.42%(圖9)。

    圖9 通過控制半徑設(shè)計出的具有不同等效剛度梯度的圓環(huán)系統(tǒng)受沖擊時的反作用力Fig.9 Reaction forces of ring system with different equivalent stiffness gradients under impact designed by controlling radius

    2.3 壁厚分布對防護(hù)特性的影響

    剛度梯度分布也取決于圓環(huán)的壁厚,討論了不同壁厚對梯度圓環(huán)系統(tǒng)防護(hù)性能的影響,其中/與2.2節(jié)保持相同。圖10繪制了在相同的沖擊條件下,具有不同壁厚分布的仿生剛度梯度圓環(huán)系統(tǒng)與內(nèi)部結(jié)構(gòu)碰撞的作用力曲線。I型梯度和Ⅱ型梯度的碰撞力-時間曲線基本上保持相同的趨勢。與2.2節(jié)結(jié)果一致,Ⅱ型梯度(LG、SQ-Ⅱ、SQ-Ⅱ、LIN)產(chǎn)生的碰撞力曲線具有2個峰值(見圖10),但是,第2峰值力大大增加,并且當(dāng)其峰值過高時,可能會對內(nèi)部結(jié)構(gòu)造成嚴(yán)重破壞。以上現(xiàn)象也可以通過梯度圓環(huán)系統(tǒng)中的質(zhì)量分布來解釋。圓環(huán)系統(tǒng)遠(yuǎn)端的慣性隨著壁厚的減小而不斷減小,因此,即使圓環(huán)具有相同的剛度分布,壁更薄的圓環(huán)更容易被壓縮更大的深度并導(dǎo)致更大的反作用力峰值。

    圖10 通過控制壁厚設(shè)計出的具有不同等效剛度梯度的圓環(huán)系統(tǒng)受沖擊時的反作用力Fig.10 Reaction forces of ring system with different equivalent stiffness gradients under impact designed by controlling wall thickness

    與具有不同彈性模量分布和不同半徑分布的梯度圓環(huán)系統(tǒng)在沖擊載荷下的動態(tài)響應(yīng)一致,通過調(diào)整厚度分布得到具有Ⅰ型剛度梯度分布的圓環(huán)系統(tǒng)會對內(nèi)部結(jié)構(gòu)產(chǎn)生更小的峰值力。結(jié)果表明,具有EXP剛度梯度的圓環(huán)系統(tǒng)對內(nèi)部結(jié)構(gòu)產(chǎn)生的碰撞力峰值相比SIN-Ⅰ梯度下降了4.43%,相比SQ-Ⅰ梯度下降了2.58%,可以得到與前面各節(jié)相同的結(jié)論。

    2.4 不同單列仿生剛度梯度圓環(huán)的防護(hù)特性比較

    為便于描述,將具有不同彈性模量分布的梯度圓環(huán)系統(tǒng)稱為E型,將具有不同半徑分布的梯度圓環(huán)系統(tǒng)稱為R型,將具有不同壁厚分布的梯度圓環(huán)系統(tǒng)稱為T型。

    如2.1~2.3節(jié)中所述,具有EXP剛度分布的圓環(huán)系統(tǒng)在相同沖擊條件下對內(nèi)部防護(hù)結(jié)構(gòu)產(chǎn)生的碰撞力最小,這與天然甲蟲角質(zhì)層結(jié)構(gòu)演化的結(jié)果一致??紤]到結(jié)構(gòu)的制造成本和重量,仿生剛度梯度圓環(huán)系統(tǒng)的質(zhì)量是一個重要的設(shè)計指標(biāo)。圖11分別給出了具有不同彈性模量,半徑和厚度分布的每種梯度圓環(huán)系統(tǒng)對內(nèi)部防護(hù)結(jié)構(gòu)產(chǎn)生的碰撞力峰值及其自身總質(zhì)量。圖11(a)~圖11(c)顯示,E型、R型、T型梯度圓環(huán)系統(tǒng)對內(nèi)部結(jié)構(gòu)產(chǎn)生的碰撞力峰值隨分布形式改變具有相同的變化趨勢,其中EXP剛度梯度導(dǎo)致最小的碰撞反力。E型梯度圓環(huán)系統(tǒng)的質(zhì)量不變(見圖12(a)),R型梯度圓環(huán)系統(tǒng)的質(zhì)量隨梯度分布形式改變的趨勢與T型梯度圓環(huán)系統(tǒng)相反(見圖11(e)、圖11(f))。這是因為半徑與厚度分別為式(1)的分母和分子,當(dāng)分別調(diào)節(jié)半徑和厚度這2個參數(shù)得到相同的等效剛度時,半徑和厚度的增大和減小趨勢相反,這使得圓環(huán)系統(tǒng)最終獲得的質(zhì)量正好相反。同樣對于EXP剛度梯度分布,T型梯度圓環(huán)系統(tǒng)對內(nèi)部結(jié)構(gòu)產(chǎn)生的碰撞力峰值比E型梯度圓環(huán)系統(tǒng)產(chǎn)生的大1.86%,而質(zhì)量卻比E型梯度圓環(huán)系統(tǒng)降低了22.4%;比R型對內(nèi)部結(jié)構(gòu)產(chǎn)生的碰撞力峰值大3.91%,而質(zhì)量下降了更多,達(dá)到46.4%。因此,當(dāng)工程師設(shè)計人工合成先進(jìn)防御結(jié)構(gòu)時,應(yīng)同時考慮其防護(hù)性能,質(zhì)量以及制造成本來選擇一種最佳折衷策略。

    圖11 具有不同等效剛度梯度分布的圓環(huán)系統(tǒng)在相同沖擊載荷作用下產(chǎn)生的反作用力峰值和自身質(zhì)量Fig.11 Reaction force peak and its own mass of ring system with different equivalent stiffness gradient distributions under the same impact load

    本文還分別計算了E型、R型、T型梯度圓環(huán)系統(tǒng)對內(nèi)部結(jié)構(gòu)產(chǎn)生的碰撞力沖量(見表1)。研究發(fā)現(xiàn),具有指數(shù)剛度分布(EXP)的T型梯度圓環(huán)系統(tǒng)對內(nèi)部防護(hù)結(jié)構(gòu)產(chǎn)生的碰撞力沖量最小。以上結(jié)果表明,具有指數(shù)剛度分布(EXP)的T型梯度圓環(huán)系統(tǒng)可能會成為工程應(yīng)用中最佳的選擇,因為它以最輕的質(zhì)量使得對內(nèi)部防護(hù)結(jié)構(gòu)產(chǎn)生的碰撞沖量達(dá)到最小。

    表1 E型、R型、T型梯度圓環(huán)系統(tǒng)對內(nèi)部結(jié)構(gòu)產(chǎn)生的碰撞力沖量Table 1 Impulses of impact force generated by E-type, R-type and T-type gradient ring systems on internal structure

    2.5 碰撞能量傳遞關(guān)系

    以具有指數(shù)剛度分布(EXP)的T型梯度圓環(huán)系統(tǒng)為例,計算在碰撞過程中的能量傳遞關(guān)系,來說明抗沖擊結(jié)構(gòu)和碰撞物的能量是如何轉(zhuǎn)換的。本算例中,圓環(huán)系統(tǒng)的參數(shù)與2.3節(jié)相同。質(zhì)量為0.2 kg、初速度為1 m/s的質(zhì)量塊撞擊圓環(huán)系統(tǒng)后,質(zhì)量塊變?yōu)榉聪蜻\動,動能下降(見圖12(a)),耗散的動能轉(zhuǎn)化為圓環(huán)系統(tǒng)的變形能和振動動能,這2部分能量組成圓環(huán)系統(tǒng)的總能量。當(dāng)質(zhì)量塊撞擊最外部圓環(huán)(圓環(huán)1)時,質(zhì)量塊動能開始減小,圓環(huán)1能量增加。隨后,碰撞能量依次向內(nèi)部傳遞,位置越靠內(nèi)的圓環(huán)吸收能量越高(見圖12(b));當(dāng)能量首次傳遞到最內(nèi)部圓環(huán)時,由于最內(nèi)部圓環(huán)與內(nèi)部結(jié)構(gòu)發(fā)生碰撞,從而反向壓縮較外部圓環(huán),使得能量由內(nèi)向外傳遞。與由外向內(nèi)傳遞時趨勢相反,由內(nèi)向外傳遞時各圓環(huán)吸收的能量逐次降低,最終各圓環(huán)的能量趨于穩(wěn)定。在碰撞過程中,質(zhì)量塊的動能先逐漸降低至零,之后反向加速,動能逐漸增加,等碰撞結(jié)束系統(tǒng)穩(wěn)定后的質(zhì)量塊動能大約減少20%。

    圖12 碰撞過程中各部分能量Fig.12 Energy of each part during collision

    3 仿生圓環(huán)陣列的防護(hù)特性分析

    根據(jù)以上研究,可以知道具有EXP剛度梯度的單列圓環(huán)系統(tǒng)具有最佳的防護(hù)性能。在此基礎(chǔ)上,進(jìn)一步研究了圓環(huán)陣列結(jié)構(gòu)(見圖13)的抗沖擊特性。以上設(shè)計策略可以應(yīng)用于各種尺寸比例和更多的工程防護(hù)結(jié)構(gòu)中。例如,在撞擊速度為1 m/s的質(zhì)量塊沖擊作用下,研究了仿生圓環(huán)陣列結(jié)構(gòu)的防護(hù)性能,撞擊圓環(huán)陣列結(jié)構(gòu)的質(zhì)量塊質(zhì)量為0.2 kg。圓環(huán)陣列中環(huán)的參數(shù)(彈性模量=210 GPa、半徑=10 mm、寬度=12.5 mm),均保持不變。壁厚的分布使每個圓環(huán)陣列的剛度梯度滿足指數(shù)(EXP)分布。

    圖13 仿生梯度圓環(huán)陣列防護(hù)系統(tǒng)Fig.13 Bionic gradient ring array protection system

    圖14顯示了在具有指數(shù)剛度梯度(EXP)的圓環(huán)陣列結(jié)構(gòu)中,彈性應(yīng)力傳播過程的數(shù)值模擬結(jié)果。在撞擊開始時,沖擊力主要作用在撞擊質(zhì)量塊正下方的3列圓環(huán)系統(tǒng)上。隨著沖擊過程的進(jìn)行,質(zhì)量塊正下方兩側(cè)的圓環(huán)系統(tǒng)也開始承受沖擊載荷。但是,數(shù)值結(jié)果表明,單列圓環(huán)系統(tǒng)離撞擊質(zhì)量塊的距離越遠(yuǎn),受沖擊載荷的影響越小。

    在沖擊應(yīng)力傳播過程中,最大應(yīng)力值沿沖擊方向逐漸減小(見圖14)。當(dāng)沖擊載荷作用于圓環(huán)陣列結(jié)構(gòu)時,最外層圓環(huán)產(chǎn)生的最大應(yīng)力為145.8 MPa,而當(dāng)應(yīng)力波傳輸?shù)皆撽嚵凶顑?nèi)層的圓環(huán)時,圓環(huán)陣列中最大應(yīng)力僅為71.6 MPa,降低了50.9%。這表明該圓環(huán)陣列結(jié)構(gòu)對彈性應(yīng)力波的傳播具有良好的衰減效果,可以減小結(jié)構(gòu)內(nèi)環(huán)的應(yīng)力水平防止該陣列結(jié)構(gòu)發(fā)生破壞。另外,圖14顯示在沖擊方向上具有指數(shù)剛度梯度(EXP)的圓環(huán)陣列結(jié)構(gòu)可以使彈性應(yīng)力波沿著沖擊方向的傳播減速。為了清晰地解釋這個問題,研究沿沖擊方向?qū)鞑ゾ嚯x分為3個相等的部為了比較沿沖擊方向具有不同剛度梯度分布的仿生圓環(huán)陣列結(jié)構(gòu)的防護(hù)性能,本文還分別計算了圓環(huán)陣列結(jié)構(gòu)與內(nèi)部基板之間的碰撞力(見圖15)。分析表明,這些受生物啟發(fā)的圓環(huán)陣列結(jié)構(gòu)可以幫助降低作用在被保護(hù)結(jié)構(gòu)(基板)上的碰撞力。研究證明凹形剛度梯度(EXP、SIN-Ⅰ、SQ-Ⅰ)具有更好的防護(hù)效果,在凹形剛度梯度中,EXP剛度分布導(dǎo)致產(chǎn)生更小的碰撞力峰值。在相同沖擊條件下,與沒有受任何保護(hù)的基板相比,作用在受仿生剛度梯度(EXP)圓環(huán)陣列結(jié)構(gòu)保護(hù)的基板上的碰撞力峰值降低了96%。盡管被保護(hù)的基板受沖擊的時間增加,但作用在基板上的沖量總量減少了68%。

    圖14 仿生梯度圓環(huán)陣列防護(hù)系統(tǒng)經(jīng)受撞擊不同時刻應(yīng)力圖Fig.14 Stress diagram of bionic gradient ring array protection system at different moments of impact

    圖15 仿生圓環(huán)陣列結(jié)構(gòu)與底部基座之間的碰撞力曲線Fig.15 Collision force curves between bionic ring array structure and bottom base

    分:第1部分包括圖13所示的1、2、3號圓環(huán),第2部分包括圖13所示的4、5、6號圓環(huán),第3部分包括圖13所示的7、8、9號圓環(huán)。彈性波經(jīng)過第1~3部分傳播路徑分別花費0.04、0.07、0.11 ms(見圖14),證明應(yīng)力波速沿沖擊方向降低。

    為說明仿生梯度圓環(huán)的抗斜撞擊特性,本文計算了質(zhì)量塊與豎直方向成45°角、以1 m/s的初速度撞擊梯度圓環(huán)陣列的工況。結(jié)果表明,在同時期圓環(huán)產(chǎn)生的應(yīng)力有所降低(見圖16(a)),這是因為產(chǎn)生斜撞擊后,縱向的撞擊速度減小造成的。但是由于橫向載荷的作用,在橫向產(chǎn)生的變形更大(見圖16(b))。由于同縱列各圓環(huán)相連但各縱列圓環(huán)之間存在間隙,因此橫向變形沒有傳遞更遠(yuǎn),這對結(jié)構(gòu)的防護(hù)能力是有利的。

    圖16 仿生梯度圓環(huán)陣列防護(hù)系統(tǒng)經(jīng)受斜撞擊時應(yīng)力和變形圖Fig.16 Stress and deformation diagram of bionic gradient ring array protection system subjected to oblique impact

    此外,考慮相鄰縱列圓環(huán)之間的橫向擠壓非線性特點,建立縱列無間隔圓環(huán)陣列,并與有間隔圓環(huán)陣列的抗沖擊特性進(jìn)行對比分析。結(jié)果表明,考慮了縱列圓環(huán)之間的橫向擠壓非線性特點的圓環(huán)陣列,受到相同撞擊載荷時,應(yīng)力場分布更廣,但是最大應(yīng)力減小(見圖17)。這是因為,圓環(huán)之間的非線性橫向擠壓作用耗散了更多沖擊能量,使圓環(huán)內(nèi)部應(yīng)力水平降低。此外,圓環(huán)之間的非線性橫向擠壓作用還降低了圓環(huán)陣列與內(nèi)部結(jié)構(gòu)的碰撞力(見圖18),提升了防護(hù)效果。

    圖17 考慮橫向擠壓和非線性變形特點的應(yīng)力圖(T=0.119 95 ms)Fig.17 Strain diagram considering transverse extrusion and nonlinear deformation (T=0.119 95 ms)

    圖18 考慮橫向擠壓和非線性變形特點的碰撞反力曲線Fig.18 Collision reaction curve considering transverse extrusion and nonlinear deformation characteristics

    4 結(jié) 論

    1) 分析了具有不同剛度梯度的單列圓環(huán)系統(tǒng)在沖擊載荷下的動態(tài)響應(yīng),結(jié)果表明,凹形剛度梯度可以顯著改善仿生圓環(huán)系統(tǒng)的防護(hù)性能,使圓環(huán)系統(tǒng)內(nèi)部始終保持較低的應(yīng)力狀態(tài)并且降低了施加到內(nèi)部受保護(hù)結(jié)構(gòu)上的碰撞力和沖擊沖量。

    2) 通過研究圓環(huán)彈性模量、半徑和厚度分布對仿生剛度梯度圓環(huán)系統(tǒng)防護(hù)特性的影響,結(jié)果表明,通過調(diào)整圓環(huán)厚度分布達(dá)到相同等效剛度梯度的形式可以達(dá)到最佳的防護(hù)效果,它可以使整個圓環(huán)系統(tǒng)總質(zhì)量最輕,同時對內(nèi)部防護(hù)結(jié)構(gòu)產(chǎn)生的碰撞沖量最小。

    3) 通過調(diào)整圓環(huán)厚度分布進(jìn)一步研究了在沖擊載荷作用下仿生圓環(huán)陣列結(jié)構(gòu)的防護(hù)性能,結(jié)果表明,具有指數(shù)剛度梯度(EXP)的仿生圓環(huán)陣列結(jié)構(gòu)可以使自身的應(yīng)力水平達(dá)到最低,并且對內(nèi)部防護(hù)結(jié)構(gòu)的沖擊影響減到最小。提出的仿生策略顯著提高了防護(hù)結(jié)構(gòu)的抗沖擊性能和防護(hù)效果,與沒有防護(hù)的情況相比,可大大減小內(nèi)部防護(hù)結(jié)構(gòu)受到損傷的風(fēng)險。

    猜你喜歡
    抗沖擊圓環(huán)梯度
    加權(quán)全能量最小的圓環(huán)形變
    一個改進(jìn)的WYL型三項共軛梯度法
    豬圓環(huán)病毒病的發(fā)生、診斷和防治
    一例鴨圓環(huán)病毒病的診斷
    內(nèi)置加勁環(huán)T型管節(jié)點抗沖擊承載力計算
    一種自適應(yīng)Dai-Liao共軛梯度法
    一類扭積形式的梯度近Ricci孤立子
    圓環(huán)上的覆蓋曲面不等式及其應(yīng)用
    ACR抗沖擊改性劑的合成及其在聚乳酸改性中的應(yīng)用研究
    中國塑料(2015年6期)2015-11-13 03:03:09
    一種抗沖擊減震橡膠
    国产熟女午夜一区二区三区 | 国产永久视频网站| 国产又色又爽无遮挡免| 妹子高潮喷水视频| 视频中文字幕在线观看| 成人国产麻豆网| 国产熟女欧美一区二区| 成人午夜精彩视频在线观看| 少妇人妻精品综合一区二区| 久久综合国产亚洲精品| 搡女人真爽免费视频火全软件| 三级国产精品欧美在线观看| 久热这里只有精品99| 女的被弄到高潮叫床怎么办| 亚洲精品国产av成人精品| 91久久精品国产一区二区三区| 亚洲美女视频黄频| 亚洲av男天堂| 久久99精品国语久久久| 国产熟女午夜一区二区三区 | 日韩成人伦理影院| av女优亚洲男人天堂| 国产黄色视频一区二区在线观看| 久久国内精品自在自线图片| 国产黄片视频在线免费观看| 在线观看免费高清a一片| 一区二区三区免费毛片| 久久精品夜色国产| tube8黄色片| 免费看光身美女| av在线app专区| 美女cb高潮喷水在线观看| 久久影院123| 男女边吃奶边做爰视频| 国产一区二区在线观看日韩| 2022亚洲国产成人精品| 精品人妻偷拍中文字幕| 精品国产一区二区三区久久久樱花| 精品国产一区二区三区久久久樱花| 精品一品国产午夜福利视频| 国产成人精品一,二区| 2021少妇久久久久久久久久久| 丰满少妇做爰视频| av国产久精品久网站免费入址| 日本av手机在线免费观看| 久久午夜综合久久蜜桃| 久久久久久久久大av| 大香蕉久久成人网| 欧美精品一区二区大全| 99热这里只有精品一区| 丁香六月天网| 欧美亚洲日本最大视频资源| 在线观看www视频免费| 丁香六月天网| 精品人妻在线不人妻| 国产精品.久久久| 99国产综合亚洲精品| 亚洲性久久影院| 亚洲国产精品999| 99久久精品一区二区三区| 男人添女人高潮全过程视频| 精品久久久噜噜| 精品人妻在线不人妻| 尾随美女入室| 成人亚洲欧美一区二区av| 欧美老熟妇乱子伦牲交| 三级国产精品欧美在线观看| 亚洲av综合色区一区| 最近中文字幕高清免费大全6| 国产亚洲午夜精品一区二区久久| 精品国产一区二区久久| 自拍欧美九色日韩亚洲蝌蚪91| 卡戴珊不雅视频在线播放| 久久综合国产亚洲精品| 国产亚洲午夜精品一区二区久久| 亚洲精品av麻豆狂野| 亚洲精品乱久久久久久| 天堂8中文在线网| 久久久久久伊人网av| 丰满乱子伦码专区| 精品久久久久久电影网| 男女啪啪激烈高潮av片| 国产精品一国产av| www.色视频.com| 一级毛片电影观看| 国产一区二区在线观看日韩| 大码成人一级视频| 91久久精品电影网| 成年女人在线观看亚洲视频| 精品人妻在线不人妻| 日日啪夜夜爽| 精品一区二区三卡| 久热这里只有精品99| 在线播放无遮挡| 国产在线视频一区二区| videosex国产| 国产精品.久久久| 免费av不卡在线播放| a级毛片在线看网站| 欧美xxⅹ黑人| 草草在线视频免费看| 日韩av免费高清视频| 2018国产大陆天天弄谢| 精品久久蜜臀av无| 欧美国产精品一级二级三级| 欧美另类一区| 色吧在线观看| 国产一区二区三区av在线| 伊人亚洲综合成人网| 亚洲欧美成人精品一区二区| 国产av码专区亚洲av| 亚洲精品久久久久久婷婷小说| 色吧在线观看| 亚洲色图综合在线观看| av国产久精品久网站免费入址| 青春草国产在线视频| 日本91视频免费播放| 高清在线视频一区二区三区| 亚洲五月色婷婷综合| 最后的刺客免费高清国语| 国产国拍精品亚洲av在线观看| 中文字幕最新亚洲高清| 97超视频在线观看视频| 女人久久www免费人成看片| 老司机影院毛片| 成人毛片a级毛片在线播放| 国产精品99久久久久久久久| 亚洲丝袜综合中文字幕| 成人漫画全彩无遮挡| 五月伊人婷婷丁香| 插逼视频在线观看| kizo精华| 91成人精品电影| 99热6这里只有精品| 内地一区二区视频在线| 国产 精品1| 国产黄片视频在线免费观看| 国产高清有码在线观看视频| 美女视频免费永久观看网站| 熟女av电影| 亚洲国产日韩一区二区| 只有这里有精品99| 草草在线视频免费看| 国产欧美另类精品又又久久亚洲欧美| 日韩强制内射视频| 午夜91福利影院| 纯流量卡能插随身wifi吗| 亚洲不卡免费看| 交换朋友夫妻互换小说| 国产精品人妻久久久影院| 一本大道久久a久久精品| 又黄又爽又刺激的免费视频.| 国产毛片在线视频| 久久久亚洲精品成人影院| 男女高潮啪啪啪动态图| 女性被躁到高潮视频| 在线观看人妻少妇| 成人亚洲欧美一区二区av| 十分钟在线观看高清视频www| 91在线精品国自产拍蜜月| 国产精品久久久久久精品古装| 在线观看人妻少妇| 久久99一区二区三区| 国产国拍精品亚洲av在线观看| 九色亚洲精品在线播放| 大又大粗又爽又黄少妇毛片口| 国产成人精品在线电影| 亚洲精品日韩av片在线观看| 久久ye,这里只有精品| 亚洲精品成人av观看孕妇| 亚洲三级黄色毛片| 国产无遮挡羞羞视频在线观看| 黄色欧美视频在线观看| 国产成人精品在线电影| 国产精品偷伦视频观看了| 亚洲av福利一区| 大片免费播放器 马上看| 另类精品久久| 少妇的逼水好多| 插逼视频在线观看| 久久久久精品久久久久真实原创| 一级二级三级毛片免费看| 色婷婷久久久亚洲欧美| 新久久久久国产一级毛片| 九草在线视频观看| 国产成人精品无人区| 日韩大片免费观看网站| 大片免费播放器 马上看| 欧美97在线视频| 一区二区三区精品91| 热re99久久国产66热| 男人添女人高潮全过程视频| 视频区图区小说| 久久久久网色| 国产免费现黄频在线看| 日韩强制内射视频| 亚洲欧美色中文字幕在线| 大香蕉久久成人网| 亚洲四区av| 国产精品不卡视频一区二区| 国产一区二区在线观看av| a级毛片黄视频| 在线精品无人区一区二区三| 日本猛色少妇xxxxx猛交久久| 欧美成人精品欧美一级黄| 亚洲第一区二区三区不卡| 久久影院123| 亚洲少妇的诱惑av| 哪个播放器可以免费观看大片| 我要看黄色一级片免费的| 国产伦理片在线播放av一区| 一级片'在线观看视频| 少妇人妻精品综合一区二区| 午夜福利在线观看免费完整高清在| 成人影院久久| 久久精品国产a三级三级三级| 天天躁夜夜躁狠狠久久av| 黄色毛片三级朝国网站| 在线播放无遮挡| 亚洲av.av天堂| 国产日韩欧美视频二区| 精品国产露脸久久av麻豆| av一本久久久久| 国产深夜福利视频在线观看| 91aial.com中文字幕在线观看| 春色校园在线视频观看| 九草在线视频观看| 丰满迷人的少妇在线观看| 日日啪夜夜爽| av福利片在线| 最近的中文字幕免费完整| 亚洲经典国产精华液单| 国国产精品蜜臀av免费| 亚洲美女搞黄在线观看| 国产极品天堂在线| 女性被躁到高潮视频| 久久久久国产网址| 久久久久人妻精品一区果冻| 亚洲三级黄色毛片| 99久久中文字幕三级久久日本| 麻豆乱淫一区二区| 永久网站在线| 91久久精品国产一区二区成人| 两个人免费观看高清视频| 亚洲怡红院男人天堂| 韩国av在线不卡| 大陆偷拍与自拍| 国产精品女同一区二区软件| 两个人免费观看高清视频| 亚洲国产精品专区欧美| 亚洲成人一二三区av| 亚洲欧美精品自产自拍| 热re99久久国产66热| 日韩在线高清观看一区二区三区| 日韩av免费高清视频| 高清毛片免费看| 视频在线观看一区二区三区| 最近中文字幕高清免费大全6| 国产精品久久久久久久久免| 日韩av在线免费看完整版不卡| 午夜福利视频精品| 国产精品人妻久久久久久| 免费播放大片免费观看视频在线观看| 最近手机中文字幕大全| 亚洲久久久国产精品| 97超视频在线观看视频| 久久久久国产精品人妻一区二区| 伦理电影大哥的女人| 国产视频首页在线观看| 男女无遮挡免费网站观看| 久久精品国产自在天天线| a 毛片基地| 夫妻性生交免费视频一级片| 国精品久久久久久国模美| 赤兔流量卡办理| 色婷婷av一区二区三区视频| 三上悠亚av全集在线观看| 男女边摸边吃奶| av国产久精品久网站免费入址| 亚洲,欧美,日韩| 少妇的逼好多水| 精品亚洲乱码少妇综合久久| 99热国产这里只有精品6| 在线免费观看不下载黄p国产| 日本免费在线观看一区| 精品国产国语对白av| 国产精品99久久久久久久久| 久久午夜综合久久蜜桃| 欧美97在线视频| 国产成人91sexporn| 亚洲成人一二三区av| 22中文网久久字幕| 大香蕉97超碰在线| 大陆偷拍与自拍| av国产精品久久久久影院| 少妇被粗大的猛进出69影院 | 免费高清在线观看视频在线观看| 欧美xxxx性猛交bbbb| 免费观看av网站的网址| 最后的刺客免费高清国语| 国产欧美亚洲国产| 国产日韩一区二区三区精品不卡 | 看非洲黑人一级黄片| 韩国av在线不卡| 两个人的视频大全免费| 一本一本综合久久| 成年美女黄网站色视频大全免费 | 国产免费福利视频在线观看| 我要看黄色一级片免费的| 女人久久www免费人成看片| 免费大片18禁| 日日摸夜夜添夜夜添av毛片| 精品久久久久久久久av| 男女无遮挡免费网站观看| 亚洲av.av天堂| 久久久久久人妻| 欧美精品一区二区免费开放| 欧美丝袜亚洲另类| 国产午夜精品久久久久久一区二区三区| 熟女人妻精品中文字幕| 午夜福利网站1000一区二区三区| 日本欧美国产在线视频| 天堂中文最新版在线下载| 卡戴珊不雅视频在线播放| 国产探花极品一区二区| 精品少妇久久久久久888优播| 日韩在线高清观看一区二区三区| 成人影院久久| 国产精品不卡视频一区二区| 天天操日日干夜夜撸| 国语对白做爰xxxⅹ性视频网站| 国产成人freesex在线| 在线观看国产h片| 99久久精品一区二区三区| 国产男人的电影天堂91| 亚洲av成人精品一区久久| 久久久久久久久大av| 成年av动漫网址| 亚洲天堂av无毛| 人妻系列 视频| 男女边摸边吃奶| 亚洲欧美中文字幕日韩二区| av网站免费在线观看视频| 亚洲四区av| 国产av码专区亚洲av| 丝袜美足系列| 国产男女超爽视频在线观看| 最后的刺客免费高清国语| 国产成人精品在线电影| 精品卡一卡二卡四卡免费| 国产欧美亚洲国产| 91久久精品国产一区二区成人| 免费看光身美女| 人人妻人人添人人爽欧美一区卜| 精品少妇黑人巨大在线播放| xxxhd国产人妻xxx| 久久精品国产亚洲av涩爱| av天堂久久9| 少妇被粗大的猛进出69影院 | 综合色丁香网| 国产成人一区二区在线| 男女国产视频网站| 熟女av电影| 日韩强制内射视频| 亚洲无线观看免费| xxxhd国产人妻xxx| 大香蕉久久网| 亚洲色图综合在线观看| 看十八女毛片水多多多| 亚洲成人手机| 国产有黄有色有爽视频| 大香蕉97超碰在线| 成人毛片a级毛片在线播放| 精品一区二区三区视频在线| 久久免费观看电影| a级毛色黄片| 国产永久视频网站| 麻豆乱淫一区二区| 欧美日韩视频精品一区| 亚洲国产成人一精品久久久| 最新中文字幕久久久久| 日日啪夜夜爽| 国产欧美另类精品又又久久亚洲欧美| 日韩电影二区| 狂野欧美激情性xxxx在线观看| 欧美变态另类bdsm刘玥| 欧美性感艳星| 一个人看视频在线观看www免费| 欧美日韩视频高清一区二区三区二| 久久av网站| 成人漫画全彩无遮挡| 午夜福利,免费看| 久久综合国产亚洲精品| 亚洲在久久综合| 久久青草综合色| 精品少妇黑人巨大在线播放| 国产免费视频播放在线视频| 最近中文字幕2019免费版| 插阴视频在线观看视频| 男女免费视频国产| 国产 精品1| 亚洲激情五月婷婷啪啪| 丰满饥渴人妻一区二区三| 91国产中文字幕| 国产精品久久久久久精品古装| 成人黄色视频免费在线看| 国产精品久久久久久久久免| 国产精品秋霞免费鲁丝片| 国产乱人偷精品视频| 桃花免费在线播放| 成人综合一区亚洲| 97在线视频观看| 赤兔流量卡办理| 成人18禁高潮啪啪吃奶动态图 | 精品久久国产蜜桃| 亚洲欧美成人精品一区二区| 免费观看性生交大片5| 欧美日本中文国产一区发布| 国产成人a∨麻豆精品| 99精国产麻豆久久婷婷| 亚洲精品日本国产第一区| 美女大奶头黄色视频| 亚洲欧美日韩另类电影网站| 国产午夜精品久久久久久一区二区三区| 嘟嘟电影网在线观看| 26uuu在线亚洲综合色| 国产爽快片一区二区三区| 久久精品久久久久久噜噜老黄| 久久精品夜色国产| 天天操日日干夜夜撸| 久久精品夜色国产| 在线观看免费日韩欧美大片 | 三级国产精品欧美在线观看| 亚洲欧美日韩另类电影网站| av免费观看日本| 熟妇人妻不卡中文字幕| 18禁裸乳无遮挡动漫免费视频| 青春草视频在线免费观看| 国产片内射在线| 久久久久久久久久成人| 大片免费播放器 马上看| 你懂的网址亚洲精品在线观看| 国内精品宾馆在线| 日韩人妻高清精品专区| 美女视频免费永久观看网站| 国产成人91sexporn| 免费少妇av软件| a级毛片免费高清观看在线播放| 亚洲欧美成人精品一区二区| 久久 成人 亚洲| 插逼视频在线观看| 国产精品99久久99久久久不卡 | 欧美成人午夜免费资源| 免费不卡的大黄色大毛片视频在线观看| 天天影视国产精品| 在线观看国产h片| 免费大片黄手机在线观看| 精品卡一卡二卡四卡免费| 成年美女黄网站色视频大全免费 | 少妇被粗大猛烈的视频| 国产色爽女视频免费观看| 亚洲欧洲日产国产| 伊人久久精品亚洲午夜| 丰满乱子伦码专区| a级片在线免费高清观看视频| 国产精品欧美亚洲77777| 亚洲无线观看免费| 美女内射精品一级片tv| 两个人的视频大全免费| 国产无遮挡羞羞视频在线观看| 亚洲av成人精品一二三区| 少妇精品久久久久久久| 黑人猛操日本美女一级片| 久久av网站| 国产精品久久久久久av不卡| 国产女主播在线喷水免费视频网站| 久久人人爽人人片av| av免费观看日本| 中文字幕人妻熟人妻熟丝袜美| 国产精品免费大片| 伊人久久精品亚洲午夜| 高清欧美精品videossex| 一级,二级,三级黄色视频| 嘟嘟电影网在线观看| 久久久久久久久大av| av专区在线播放| 久久人人爽av亚洲精品天堂| 欧美最新免费一区二区三区| 午夜激情av网站| 在线观看免费日韩欧美大片 | 日本猛色少妇xxxxx猛交久久| 亚洲性久久影院| 国产免费一区二区三区四区乱码| 伊人亚洲综合成人网| 蜜桃国产av成人99| 一本—道久久a久久精品蜜桃钙片| 成人午夜精彩视频在线观看| 免费观看无遮挡的男女| 满18在线观看网站| 热re99久久精品国产66热6| 亚洲怡红院男人天堂| 美女中出高潮动态图| 久久综合国产亚洲精品| 精品国产露脸久久av麻豆| 狂野欧美白嫩少妇大欣赏| av福利片在线| 男人添女人高潮全过程视频| 精品少妇久久久久久888优播| 熟女电影av网| 99热这里只有是精品在线观看| 精品国产一区二区久久| 男女无遮挡免费网站观看| 人人妻人人添人人爽欧美一区卜| 中文字幕人妻丝袜制服| 日本黄色片子视频| 七月丁香在线播放| 久久久欧美国产精品| 国产成人免费无遮挡视频| 国产免费一级a男人的天堂| 亚洲精品中文字幕在线视频| 不卡视频在线观看欧美| 插阴视频在线观看视频| 免费高清在线观看视频在线观看| 老熟女久久久| 国产片内射在线| 久久久久久久久久久免费av| 水蜜桃什么品种好| 成人国产麻豆网| 在线观看免费视频网站a站| 亚洲色图综合在线观看| 纵有疾风起免费观看全集完整版| 晚上一个人看的免费电影| 亚州av有码| 亚洲国产精品999| 一级毛片我不卡| 亚洲精品日韩在线中文字幕| 熟妇人妻不卡中文字幕| 免费av中文字幕在线| 免费大片黄手机在线观看| 久久久久网色| 欧美bdsm另类| 秋霞伦理黄片| 男男h啪啪无遮挡| 交换朋友夫妻互换小说| 中文字幕免费在线视频6| 国产男女内射视频| 国产精品.久久久| 欧美精品亚洲一区二区| 草草在线视频免费看| 欧美变态另类bdsm刘玥| videossex国产| 亚洲精品久久久久久婷婷小说| 亚洲少妇的诱惑av| 丰满迷人的少妇在线观看| 全区人妻精品视频| 久久久精品区二区三区| 黑人高潮一二区| 99国产精品免费福利视频| 免费日韩欧美在线观看| 久久精品国产自在天天线| 在线观看免费视频网站a站| 91精品国产国语对白视频| 日韩人妻高清精品专区| 亚洲精品一区蜜桃| 成人漫画全彩无遮挡| 狠狠婷婷综合久久久久久88av| 国产无遮挡羞羞视频在线观看| 精品久久蜜臀av无| 久久久久国产网址| 亚洲欧美一区二区三区国产| www.色视频.com| 久久av网站| 亚洲av免费高清在线观看| 日韩中字成人| av视频免费观看在线观看| 欧美日韩国产mv在线观看视频| 只有这里有精品99| 午夜影院在线不卡| 久久97久久精品| 日韩一本色道免费dvd| 极品少妇高潮喷水抽搐| 综合色丁香网| 中文天堂在线官网| 亚洲精品一区蜜桃| 黄片播放在线免费| 国产免费一级a男人的天堂| 成年av动漫网址| 久久午夜综合久久蜜桃| 成年美女黄网站色视频大全免费 | av又黄又爽大尺度在线免费看| 国产成人freesex在线| 国产 精品1| 久久久久久伊人网av| 在线精品无人区一区二区三| 久久久久精品性色| 九九久久精品国产亚洲av麻豆| kizo精华| 一本一本综合久久| 日韩不卡一区二区三区视频在线| 亚洲精华国产精华液的使用体验| 一本大道久久a久久精品| 美女主播在线视频| 国产午夜精品一二区理论片| 日韩欧美精品免费久久| 中文欧美无线码| 一级毛片aaaaaa免费看小| 亚洲高清免费不卡视频| 国产精品偷伦视频观看了| 美女中出高潮动态图| 久久午夜福利片| 国产一区二区三区av在线| 久久国产精品男人的天堂亚洲 | 久久精品国产a三级三级三级| 9色porny在线观看| 91久久精品电影网|