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

    可展薄膜的Miura彈性折痕建模與分析

    2023-01-11 10:23:32唐愉真郭宏偉王治易劉榮強(qiáng)
    關(guān)鍵詞:折痕構(gòu)型傾角

    唐愉真,劉 超,肖 洪,郭宏偉,王治易,謝 超,劉榮強(qiáng)

    (1.機(jī)器人技術(shù)與系統(tǒng)國(guó)家重點(diǎn)實(shí)驗(yàn)室(哈爾濱工業(yè)大學(xué)),哈爾濱 150001;2.中國(guó)航天科技集團(tuán)有限公司 空間結(jié)構(gòu)與機(jī)構(gòu)技術(shù)實(shí)驗(yàn)室,上海 201108)

    隨著航天事業(yè)的發(fā)展,深空探測(cè)、地球觀測(cè)、射電天文學(xué)等航天任務(wù)對(duì)航天器的結(jié)構(gòu)與可展機(jī)構(gòu)的尺寸、質(zhì)量、可靠性等提出了更高的要求??烧归_(kāi)薄膜具有質(zhì)量輕、尺寸大等一系列優(yōu)點(diǎn),其廣泛應(yīng)用于太陽(yáng)帆、大孔徑天線、遮陽(yáng)罩中,如NASA 研制的NanoSail-D太陽(yáng)帆[1],R3D2衛(wèi)星中的P-DAHGR 薄膜反射陣面天線[2]等。隨著空間任務(wù)對(duì)大尺寸薄膜的需求越來(lái)越迫切,如何實(shí)現(xiàn)可展薄膜的大折展比折疊,折疊后如何實(shí)現(xiàn)高效有序展開(kāi),以及展開(kāi)后如何保證薄膜表面的高精度成為研究熱點(diǎn)之一。

    Miura-ori是由Miura[3-4]提出的一種折紙方法,該研究的初衷是致力于解決地圖的折疊問(wèn)題,目前這種折紙方法廣泛應(yīng)用于航空航天、建筑等眾多領(lǐng)域。如Wright等[5]開(kāi)發(fā)了一種使用Miura-ori折疊方法的裝配過(guò)程,以及一個(gè)用于制造3 m× 3 m太陽(yáng)帆的原型裝配臺(tái);Velez等[6]提出了一種基于擴(kuò)展Miura-ori折紙圖形的可折疊可重構(gòu)深空通信天線;Yu等[7]設(shè)計(jì)了一種基于Miura-ori的氣動(dòng)可折疊驅(qū)動(dòng)器驅(qū)動(dòng)的爬行軟體機(jī)器人;Li 等[8]根據(jù)Miura-Ori原理提出了一種設(shè)計(jì)和制造可伸縮電路板的方法;Horner等[9]詳細(xì)討論了Miura-Ori太陽(yáng)帆模型的構(gòu)造技術(shù),并概述了各種可能的展開(kāi)方法;肖洪等[10]基于 Miura-ori折紙?jiān)硖岢霰∧し雷o(hù)罩折疊方法,并進(jìn)行了參數(shù)分析。

    薄膜在折疊過(guò)程中局部發(fā)生變形,由此產(chǎn)生折痕[11]。折痕改變了膜面材料的狀態(tài)與行為,且在薄膜完全展開(kāi)狀態(tài)下具有一定的不可恢復(fù)性[12-13],對(duì)折疊態(tài)薄膜展開(kāi)過(guò)程與展開(kāi)結(jié)果產(chǎn)生了一定程度的影響。在研究折痕的彈塑性變形方面,國(guó)內(nèi)外許多學(xué)者將薄膜折痕處假設(shè)為完全塑性變形,以便對(duì)折痕區(qū)域展開(kāi)研究。Woo等[14]采用有限元法模擬薄膜折疊和單軸拉伸試驗(yàn),以計(jì)算薄膜的有效模量;Furuya等[15-16]建立了由彈性鉸鏈和板組成的一維可展開(kāi)薄膜分析模型,通過(guò)數(shù)值分析研究了薄膜的抗彎剛度對(duì)結(jié)構(gòu)的影響,還對(duì)一維Z型折疊薄膜進(jìn)行了數(shù)值模擬和實(shí)驗(yàn)驗(yàn)證,將大變形理論作為彈性問(wèn)題進(jìn)行理論分析,研究了折疊膜結(jié)構(gòu)太陽(yáng)帆的展開(kāi)性能;Satou等[17-18]采用有限元方法研究可展開(kāi)薄膜的折疊特性,估計(jì)了折痕中可能出現(xiàn)的塑性變形,還對(duì)一維Z型折疊薄膜進(jìn)行了彈塑性有限元分析,檢驗(yàn)折痕的力學(xué)性能;Reynolds等[19]對(duì)光學(xué)應(yīng)用的薄膜提出了一種彈性折疊方法,使用特定的最小半徑折痕而不是剛性折痕;侯德高等[20]對(duì)薄硅橡膠板的交叉折疊進(jìn)行了數(shù)值分析,研究了薄膜的變形和應(yīng)力;謝志民等[21]對(duì)柔性復(fù)合薄膜材料的折疊損傷進(jìn)行了實(shí)驗(yàn)研究,并討論了由于壓縮屈曲導(dǎo)致的復(fù)合薄膜在拉伸中出現(xiàn)的伸長(zhǎng)破壞的現(xiàn)象。上述研究的薄膜折痕多數(shù)以一維折痕為研究對(duì)象,很少涉及二維折痕的研究,而且,目前對(duì)薄膜展開(kāi)過(guò)程的分析大部分停留在理論分析階段,沒(méi)有實(shí)現(xiàn)相應(yīng)的工程應(yīng)用。深入研究折痕處的力學(xué)行為對(duì)薄膜展開(kāi)過(guò)程是否順暢以及展開(kāi)后的表面精度至關(guān)重要,二維彈性折痕的可行性分析具有重大的研究意義。

    本文提出基于Miura彈性折痕的可展薄膜構(gòu)型,可以兼顧大尺度、輕量化、高展收可靠性等性能要求,并對(duì)折痕參數(shù)作出合理優(yōu)化以實(shí)現(xiàn)大折展比及小折痕總長(zhǎng),同時(shí)建立數(shù)值模型,分析折痕在折疊過(guò)程中力學(xué)行為變化,驗(yàn)證薄膜二維彈性折痕的可行性,使得展開(kāi)后薄膜的平整度及表面精度達(dá)到最優(yōu)。

    1 空間薄膜可展機(jī)構(gòu)構(gòu)型設(shè)計(jì)

    空間薄膜可展開(kāi)機(jī)構(gòu)通常由展開(kāi)支撐機(jī)構(gòu)、薄膜和張拉系統(tǒng)組成?;诹呅蔚膸缀螛?gòu)型,根據(jù)空間薄膜可展機(jī)構(gòu)具有薄膜結(jié)構(gòu)可折疊、輕量化、高折展比、高剛度和大尺寸等性能要求,本文提出了中心展開(kāi)空間薄膜可展機(jī)構(gòu),如圖1所示。中心展開(kāi)的展開(kāi)方式具有良好的同步展開(kāi)性并且易于控制,同時(shí)在面密度、收攏展開(kāi)尺寸和動(dòng)力學(xué)特性等方面也具有優(yōu)勢(shì)。

    (a)收攏狀態(tài)

    (b)展開(kāi)狀態(tài)

    中心展開(kāi)空間薄膜可展機(jī)構(gòu)選用開(kāi)縫圓管彈性伸桿作為展開(kāi)支撐機(jī)構(gòu),卷曲后能夠發(fā)生較大的彈性變形,從而存儲(chǔ)一定的彈性變形能;當(dāng)解除約束后彈性變形自動(dòng)恢復(fù),實(shí)現(xiàn)自驅(qū)動(dòng)展開(kāi),且桿件展開(kāi)狀態(tài)具有一定的彎曲剛度和強(qiáng)度。張拉系統(tǒng)采用網(wǎng)狀懸線張拉方案,通過(guò)內(nèi)外懸線以及連接線的合理配置,有效減輕了繩索總質(zhì)量,并極大程度地消除了褶皺,提高了薄膜的表面精度和平整度,有效改善了中心展開(kāi)空間薄膜可展機(jī)構(gòu)整體系統(tǒng)的動(dòng)力學(xué)特性。

    中心展開(kāi)空間薄膜可展機(jī)構(gòu)將薄膜整體分解成了6個(gè)三角形薄膜,借鑒折紙方法Miura-ori來(lái)實(shí)現(xiàn)三角形薄膜的折疊展開(kāi)。參照Miura-ori的谷線和脊線,合理配置薄膜的折痕,展開(kāi)時(shí)靠近中心的薄膜角點(diǎn)固定,彈性伸桿帶動(dòng)邊緣的薄膜角點(diǎn)進(jìn)行展開(kāi),當(dāng)彈性伸桿完全展開(kāi)時(shí)完成整個(gè)薄膜的展開(kāi)。而針對(duì)中心展開(kāi)空間薄膜可展機(jī)構(gòu)中的三角形薄膜,如何設(shè)計(jì)優(yōu)化Miura-ori折痕參數(shù),來(lái)滿足機(jī)構(gòu)性能要求,實(shí)現(xiàn)高折展比和對(duì)薄膜較小的塑性損傷,需要進(jìn)一步深入研究。

    塑性折痕會(huì)使薄膜產(chǎn)生塑性變形,甚至造成薄膜的損傷,不僅會(huì)影響薄膜的展開(kāi)效率及展開(kāi)穩(wěn)定可靠度,而且可能存在殘余應(yīng)力以及應(yīng)力集中的現(xiàn)象,導(dǎo)致薄膜產(chǎn)生褶皺乃至撕裂,影響薄膜的展開(kāi)形態(tài)及表面精度。在完成優(yōu)化的Miura-ori三角形薄膜折疊方法的基礎(chǔ)上,研究其彈性折痕的可行性。而建立折痕的解析模型來(lái)分析其彈塑性,模型很復(fù)雜且很難可視化直觀分析薄膜各處的應(yīng)力分布情況。可以通過(guò)有限元來(lái)觀察應(yīng)力分布是否存在應(yīng)力集中以及最大應(yīng)力是否超過(guò)屈服強(qiáng)度來(lái)驗(yàn)證Miura-ori彈性折痕的可行性。在得到的優(yōu)化折痕參數(shù)基礎(chǔ)上建立數(shù)值分析模型,利用有限元軟件ABAQUS/Explicit顯示分析模塊來(lái)模擬折疊過(guò)程,分析折疊后薄膜的應(yīng)力分布并驗(yàn)證Miura-ori彈性折痕的可行性。

    2 Miura-ori折痕參數(shù)分析和優(yōu)化

    Miura-ori折疊方法中的兩個(gè)主要折痕參數(shù)分別是折痕間距和折痕傾角,基于折痕參數(shù)的幾何和數(shù)值關(guān)系建立分析模型,對(duì)折痕參數(shù)進(jìn)行分析和優(yōu)化,以提高薄膜的折展比,同時(shí)通過(guò)減小折痕總長(zhǎng)來(lái)降低對(duì)薄膜的損傷。

    2.1 Miura-ori 折疊薄膜折展比計(jì)算

    Miura-ori折疊因其快速簡(jiǎn)單的展開(kāi)方式而廣泛應(yīng)用于衛(wèi)星天線和太陽(yáng)能聚光器等空間任務(wù)中,將Miura-ori折疊應(yīng)用到中心展開(kāi)空間薄膜機(jī)構(gòu)的薄膜折疊設(shè)計(jì)中。圖2(a)中繪制了Miura-ori基本單元的折痕參數(shù)及折展構(gòu)型,其構(gòu)型由折痕長(zhǎng)度a,b和折痕傾角θ這3個(gè)參數(shù)決定,并由4個(gè)二面角來(lái)描述其折展構(gòu)型狀態(tài)。由此基本單元可知,Miura-ori只需要一次折展即可完成折疊及展開(kāi),折展比大,折展效率高。中心展開(kāi)空間薄膜可展機(jī)構(gòu)由6個(gè)三角形薄膜組成,將Miura-ori應(yīng)用于三角形構(gòu)型,得到如圖2(b)所示的基于Miura-ori的三角形薄膜折痕分布圖,圖中d為折痕間距,L為薄膜邊長(zhǎng),θ為折痕傾角,s為薄膜厚度。

    (a)Miura-ori基本單元

    (b)基于Miura-ori的三角形薄膜

    為了計(jì)算基于Miura-ori的三角形薄膜的折展比,需要先計(jì)算Miura-ori折疊的展開(kāi)體積Vd和折疊體積Vf。根據(jù)幾何關(guān)系,可以計(jì)算薄膜的展開(kāi)體積為

    (1)

    基于Miura-ori折疊的三角形薄膜折疊過(guò)程和收攏狀態(tài)如圖3所示,其收攏狀態(tài)的包絡(luò)體積可被視為具有平行四邊形底部的棱柱。反之,要展開(kāi)基于Miura-ori折疊的三角形薄膜,只需要對(duì)3個(gè)角點(diǎn)進(jìn)行張拉。

    圖3 基于Miura-ori折疊的薄膜折疊過(guò)程和收攏狀態(tài)

    棱柱的高度與薄膜的厚度和折疊后的薄膜層數(shù)有關(guān),根據(jù)幾何關(guān)系和棱柱體積計(jì)算公式,推導(dǎo)基于Miura-ori折疊的三角形薄膜的折疊體積為

    (2)

    根據(jù)折展比的定義,進(jìn)一步推導(dǎo)了基于Miura-ori的三角形薄膜的折展比表達(dá)式:

    (3)

    由式(3)可知,基于Miura-ori的三角形薄膜的折展比只與薄膜的折痕傾角θ和折痕間距d有關(guān)。

    當(dāng)L=25 m,折痕間距不應(yīng)超過(guò)薄膜邊長(zhǎng),所以0 m≤d≤25 m,折痕傾角0°≤θ≤90°,繪制基于Miura-ori的三角形薄膜的折展比關(guān)于折痕參數(shù)的曲線,如圖4所示。由圖4可知,在折痕間距一定時(shí),基于Miura-ori折疊的薄膜的的折展比隨折痕傾角的增大而增大。在折痕傾角一定時(shí),折展比隨折痕間距的增大而減小。

    圖4 Miura-ori折疊的折展比關(guān)于折痕參數(shù)的曲線

    2.2 基于Miura-ori折疊薄膜折痕長(zhǎng)度計(jì)算

    (4)

    (5)

    折痕總長(zhǎng)為豎直折痕和水平傾斜折痕長(zhǎng)度總和,即

    (6)

    式中:θ為折痕傾角,L為薄膜邊長(zhǎng),d為折痕間距,ΣL為折痕總長(zhǎng)。由式(6)易知,Miura-ori折疊的折痕總長(zhǎng)與折痕間距、折痕傾角和薄膜邊長(zhǎng)有關(guān),合理設(shè)計(jì)這些參數(shù)可以有效縮短折痕長(zhǎng)度,降低對(duì)薄膜造成的損傷程度。

    2.3 薄膜折痕參數(shù)優(yōu)化

    根據(jù)薄膜可展機(jī)構(gòu)的設(shè)計(jì)要求,對(duì)折痕參數(shù)進(jìn)行優(yōu)化以實(shí)現(xiàn)高折展比、短折痕長(zhǎng)度,綜合折展比和折痕長(zhǎng)度,定義折展比與折痕長(zhǎng)度的比Θ/ΣL為優(yōu)化目標(biāo),基于Miura-ori折疊的Θ/ΣL表達(dá)式為

    (7)

    將三角形薄膜邊長(zhǎng)L=25 m帶入式(7)中進(jìn)行簡(jiǎn)化,可得最終的優(yōu)化目標(biāo)函數(shù):

    (8)

    式中其他折痕參數(shù)均為已知量,只需求解d和θ兩個(gè)設(shè)計(jì)變量的最優(yōu)解。計(jì)算后求得基于Miura-ori折疊的最優(yōu)折痕參數(shù)為θ=80°,d=0.25 m。

    3 二維彈性折痕數(shù)值建模分析

    一維彈性折痕較為容易實(shí)現(xiàn),而二維彈性折痕必然會(huì)存在折痕交點(diǎn),此處會(huì)產(chǎn)生應(yīng)力集中,容易產(chǎn)生塑性變形。為了分析二維彈性折痕,選擇最為簡(jiǎn)單的四折痕基本單元進(jìn)行分析,分別對(duì)折痕傾角θ=90°及θ<90°的四折痕基本單元建立數(shù)值仿真模型,分析二維彈性折痕的可行性。圖2(a)中所示的Miura-ori基本單元,也可看作折痕傾角θ<90°的四折痕基本單元,先對(duì)基本單元進(jìn)行數(shù)值建模分析,為后續(xù)建立多折痕的Miura-ori彈性折痕分析模型奠定基礎(chǔ)。

    3.1 折痕傾角θ=90°彈性折痕有限元仿真分析

    為更好地分析薄膜的二維彈性折痕,對(duì)折痕傾角θ=90°的四折痕基本單元進(jìn)行分析,建立如圖5所示的幾何模型,薄膜基本幾何參數(shù)為:厚度t=25 μm,有效邊長(zhǎng)L=100 mm。根據(jù)Euler-Bernoulli梁理論可知,當(dāng)薄膜的材料參數(shù)確定,其彈性折痕的曲率半徑與薄膜厚度成正比,數(shù)值模型以厚度為25 μm的聚酰亞胺薄膜為例來(lái)進(jìn)行分析。為實(shí)現(xiàn)交叉折疊,根據(jù)折痕的峰線和谷線分布,第一步要將AD、BC兩邊沿著折痕EF對(duì)折,第二步則將AB、CD兩邊沿著折痕GH對(duì)折。選用ABAQUS/Explicit顯示分析模塊來(lái)進(jìn)行薄膜折疊過(guò)程的數(shù)值分析。

    單元類型和網(wǎng)格尺寸直接決定了創(chuàng)建有限元模型的質(zhì)量,進(jìn)而影響仿真分析的計(jì)算速度和準(zhǔn)確性。為了對(duì)薄膜折疊和展開(kāi)過(guò)程進(jìn)行比較準(zhǔn)確的數(shù)值模擬,要選擇正確的單元類型以及合適的網(wǎng)格尺寸。在確定網(wǎng)格尺寸時(shí),不僅要保證薄膜折痕附近的網(wǎng)格質(zhì)量和計(jì)算精度,而且要保證合適的計(jì)算時(shí)間。根據(jù)收斂性測(cè)試,最終網(wǎng)格劃分均采用四邊形進(jìn)階算法,網(wǎng)格大小設(shè)定為1 mm,薄膜選用M3D4R單元。此外,在進(jìn)行顯示分析步時(shí),對(duì)計(jì)算時(shí)長(zhǎng)的選取也會(huì)影響收斂性及計(jì)算速度。針對(duì)薄膜二維彈性折痕,每一個(gè)折疊過(guò)程劃分為一個(gè)計(jì)算步,綜合考慮仿真的準(zhǔn)確性、收斂性以及計(jì)算速度,設(shè)定每一個(gè)分析步時(shí)間為0.1 s。綜上,在ABAQUS中建立折痕傾角θ=90°薄膜的有限元模型,如圖5 (c)所示。薄膜添加的材料為聚酰亞胺,其常溫(22 ℃)狀態(tài)下的主要特性如表1所示。

    (a) 薄膜的具體尺寸 (b)導(dǎo)入ABAQUS的幾何模型 (c)有限元模型

    表1 聚酰亞胺材料特性

    折痕傾角θ=90°的四折痕基本單元需要兩步折疊得到最終的折疊構(gòu)型,所以分別建立Step-1和Step-2兩個(gè)計(jì)算步來(lái)進(jìn)行邊界條件的添加和折疊結(jié)果的計(jì)算。為了更好地觀察折疊后的應(yīng)力分布情況,邊界條件添加到特定的節(jié)點(diǎn)上,通過(guò)控制節(jié)點(diǎn)來(lái)實(shí)現(xiàn)折疊,相較于將剛體與薄膜綁定之后控制剛體的邊界條件添加方法,節(jié)點(diǎn)控制法對(duì)薄膜的構(gòu)型影響更小,所得變形及應(yīng)力結(jié)果只跟彈性折痕有關(guān),不受剛體的影響,結(jié)果更加準(zhǔn)確。如圖6(a)所示,在薄膜的中間定義一條路徑EF以便于后處理提取路徑上各節(jié)點(diǎn)應(yīng)力值,Step-1中邊界條件為如圖6所示的節(jié)點(diǎn)集向著中心點(diǎn)O沿著X軸旋轉(zhuǎn)90°,得到如圖6(b)所示的折疊構(gòu)型。

    對(duì)Step-1的仿真結(jié)果進(jìn)行后處理,提取所定義的路徑EF上各節(jié)點(diǎn)應(yīng)力值,繪制如圖7所示的Step-1后Von Mises應(yīng)力沿路徑分布圖。由圖7可知,路徑上各節(jié)點(diǎn)應(yīng)力呈對(duì)稱分布,在路徑25 mm及75 mm處出現(xiàn)應(yīng)力最大值,為0.67 MPa;而路徑50 mm處應(yīng)力值為0.48 MPa,造成這一應(yīng)力分布的原因是,只在靠近邊緣的特定節(jié)點(diǎn)處添加了邊界條件,靠近中間處的節(jié)點(diǎn)處于自由狀態(tài),所以彈性折痕附近處的路徑EF的應(yīng)力分布應(yīng)呈對(duì)稱性,但應(yīng)力最大處不一定出現(xiàn)在路徑中點(diǎn)處,且各節(jié)點(diǎn)應(yīng)力值都很小,非常接近。經(jīng)歷Step-1折疊后薄膜的應(yīng)力值仍處于彈性范圍內(nèi),未發(fā)生塑性變形,初步驗(yàn)證了一維彈性折痕的可行性。

    (a)Step-1邊界條件

    (b)Step-1后薄膜構(gòu)型

    圖7 Step-1后Von Mises應(yīng)力沿路徑分布

    在Step-1基礎(chǔ)上,進(jìn)行第二步折疊Step-2的仿真計(jì)算。如圖8(a)所示,根據(jù)第二步折疊的特性,定義兩條路徑EF及GH以便于后處理提取路徑上各節(jié)點(diǎn)應(yīng)力值。Step-2中邊界條件為如圖8(a)所示的節(jié)點(diǎn)集向著中心點(diǎn)O沿著Y軸旋轉(zhuǎn)90°,得到如圖8(b)所示的折疊構(gòu)型,即折痕傾角θ=90°的四折痕基本單元的最終折疊構(gòu)型。

    并對(duì)Step-2的仿真結(jié)果進(jìn)行后處理,提取所定義的路徑EF及GH上各節(jié)點(diǎn)應(yīng)力值,繪制如圖9所示的Step-2后Von Mises應(yīng)力沿路徑分布圖。由圖9可知,兩條路徑上各節(jié)點(diǎn)應(yīng)力總體呈對(duì)稱分布,薄膜邊界附近應(yīng)力并不大,當(dāng)靠近兩條折痕交點(diǎn)時(shí),應(yīng)力急劇增大,并且在兩條折痕交點(diǎn)處應(yīng)力達(dá)到最大值,為2.794 MPa。通過(guò)和材料屬性對(duì)比分析可知,經(jīng)歷Step-2折疊后薄膜的應(yīng)力值仍處于彈性范圍內(nèi),未發(fā)生塑性變形,初步驗(yàn)證了二維彈性折痕的可行性。路徑EF的應(yīng)力分布對(duì)稱性較好,而路徑GH的應(yīng)力分布總體趨勢(shì)呈現(xiàn)對(duì)稱性,少量區(qū)域不完全對(duì)稱,造成這一應(yīng)力分布的原因是,OE和OF完全對(duì)稱且受力情況相同,而OG會(huì)被折疊到OH的內(nèi)側(cè),造成應(yīng)力值存在一定的差異,但是兩條路徑的最大應(yīng)力值都出現(xiàn)在點(diǎn)O處,因?yàn)辄c(diǎn)O為四折痕的交點(diǎn),符合彈塑性力學(xué)規(guī)律。

    (a)Step-2邊界條件

    (b)Step-2后薄膜構(gòu)型

    圖9 Step-2后Von Mises應(yīng)力沿路徑的應(yīng)力分布

    3.2 折痕傾角θ<90°彈性折痕有限元仿真分析

    為了進(jìn)一步分析驗(yàn)證折痕傾角θ<90°的四折痕基本單元折疊過(guò)程,建立如圖10(a)所示的幾何模型,薄膜基本幾何參數(shù)為:厚度t=25 μm,有效邊長(zhǎng)L=100 mm。為了實(shí)現(xiàn)折痕傾角θ<90°四折痕基本單元的折疊,實(shí)際上是根據(jù)折痕的脊線和谷線分布,將相鄰表面沿著他們的連接折痕折疊實(shí)現(xiàn)的,而且僅需一個(gè)步驟就能完成所有折痕的折疊。同樣選用ABAQUS/Explicit顯示分析模塊來(lái)進(jìn)行折痕傾角θ<90°四折痕基本單元折疊過(guò)程的數(shù)值分析。

    在ABAQUS中建立折痕傾角θ<90°四折痕基本單元的有限元模型,并在如圖10(b)所示的折痕處定義兩條路徑AOA′及BOB′,其中薄膜選用M3D4R單元,網(wǎng)格劃分均采用四邊形進(jìn)階算法,網(wǎng)格大小為1 mm,同樣薄膜添加的材料為聚酰亞胺。利用建立的有限元模型,對(duì)薄膜折痕傾角θ<90°四折痕基本單元的折疊過(guò)程進(jìn)行仿真分析,如上所述其折疊過(guò)程只需要一個(gè)步驟,所以建立一個(gè)分析步Step-1來(lái)進(jìn)行邊界條件的添加和折疊結(jié)果的計(jì)算即可。如圖10(c)所示,兩組相鄰的表面分別沿著傾斜折痕旋轉(zhuǎn)90°得到如圖11所示的折疊構(gòu)型。

    (a)幾何模型 (b)有限元模型 (c)邊界條件

    圖11 折痕傾角θ<90°的折疊過(guò)程

    對(duì)仿真結(jié)果進(jìn)行后處理,提取所定義的路徑AOA′及BOB′上各節(jié)點(diǎn)應(yīng)力值,分別繪制如圖12所示的Von Mises應(yīng)力分布曲線。由圖12可知,在路徑AOA′上,最大應(yīng)力值發(fā)生在40 mm處,在路徑BOB′上,最大應(yīng)力值發(fā)生在70 mm處,這兩處坐標(biāo)對(duì)應(yīng)位置均為O點(diǎn),應(yīng)力值為2.83 MPa,應(yīng)力值仍處于彈性范圍內(nèi),未發(fā)生塑性變形,初步驗(yàn)證了折痕傾角θ<90°四折痕基本單元的二維彈性折痕的可行性。

    (a)路徑AOA′

    (b)路徑BOB′

    4 Miura-ori彈性折痕仿真和樣機(jī)研制

    4.1 基于Miura-ori折疊過(guò)程有限元仿真分析

    參照上一章中的數(shù)值模型建立方法,建立了如圖13(a)所示的基于Miura-ori的三角形薄膜的幾何模型,薄膜基本幾何參數(shù)為:厚度t=25μm,有效邊長(zhǎng)L=100 mm。同樣的,Miura-ori折疊的實(shí)現(xiàn)是根據(jù)折痕的峰線和谷線分布,將相鄰表面沿著他們的連接折痕折疊,并且僅需一個(gè)步驟就能完成所有折痕的折疊。在折痕交點(diǎn)處的5個(gè)節(jié)點(diǎn)分別定義為節(jié)點(diǎn)1、2、3、4、5,以便于分析折痕交點(diǎn)處的應(yīng)力情況。同樣選用ABAQUS/Explicit顯示分析模塊來(lái)進(jìn)行基于Miura-ori薄膜折疊過(guò)程的數(shù)值分析。單元的類型、材料的屬性以及網(wǎng)格尺寸均與上一章相同,進(jìn)而得到如圖13(b)所示的基于Miura-ori薄膜彈性折疊的有限元模型。

    (a)Miura-ori幾何模型 (b) Miura-ori有限元模型 (c) Miura-ori邊界條件

    利用建立的有限元模型,對(duì)薄膜基于Miura-ori的折疊過(guò)程進(jìn)行仿真分析,其折疊過(guò)程同樣只需要一個(gè)步驟,所以建立一個(gè)分析步Step-1來(lái)進(jìn)行邊界條件的添加和折疊結(jié)果的計(jì)算即可。如圖13(c)所示,各組相鄰的表面分別沿著傾斜折痕旋轉(zhuǎn)90°,得到如圖14所示的折疊構(gòu)型。

    在上述仿真中,不難看出在折疊過(guò)程中多條折痕交點(diǎn)處的應(yīng)力相比于薄膜其他位置應(yīng)力更大,為此,提取上述5個(gè)節(jié)點(diǎn)的應(yīng)力變化曲線,如圖15所示??梢杂^察出3個(gè)特征,首先,在折疊過(guò)程進(jìn)行到一半左右,各個(gè)節(jié)點(diǎn)的應(yīng)力達(dá)到最大值,其中節(jié)點(diǎn)2的應(yīng)力峰值最大,為3.6 MPa,但折疊后薄膜的應(yīng)力值仍處于彈性范圍內(nèi),未發(fā)生塑性變形,驗(yàn)證了薄膜基于Miura-ori折疊多折痕彈性折痕的可行性;其次,可以發(fā)現(xiàn)節(jié)點(diǎn)1處的應(yīng)力峰值最小,節(jié)點(diǎn)4、5的應(yīng)力峰值與節(jié)點(diǎn)1相近,而節(jié)點(diǎn)1、4、5均為薄膜臨近邊界上的節(jié)點(diǎn),相應(yīng)可以得出結(jié)論,臨近邊界上的折痕交點(diǎn)相較于薄膜內(nèi)部的折痕交點(diǎn)在折疊過(guò)程中產(chǎn)生的應(yīng)力較??;最后,在幾何模型中,節(jié)點(diǎn)3和節(jié)點(diǎn)5是成幾何對(duì)稱的,對(duì)于所得到的應(yīng)力變化曲線中,節(jié)點(diǎn)3和節(jié)點(diǎn)5的應(yīng)力變化在折疊前0.05 s內(nèi)幾乎完全一致,在后0.05 s內(nèi)有一點(diǎn)偏差,符合應(yīng)力分布變化規(guī)律。

    圖14 Miura-ori折疊過(guò)程

    圖15 薄膜各折痕交點(diǎn)處應(yīng)力隨時(shí)間變化曲線

    4.2 空間可展開(kāi)薄膜機(jī)構(gòu)樣機(jī)研制

    根據(jù)所提出的空間薄膜可展機(jī)構(gòu)構(gòu)型設(shè)計(jì)方案,研制薄膜邊長(zhǎng)為2.5 m,厚度為25 μm的原理樣機(jī),其中薄膜基于Miura-ori的折痕參數(shù)設(shè)計(jì)優(yōu)化方法,對(duì)三角形薄膜進(jìn)行折痕參數(shù)優(yōu)化及二維彈性折痕的研制。桿件采用碳纖彈性伸桿,懸線采用凱夫拉繩索,接頭采用尼龍基玻璃纖維復(fù)合材料3D打印制作??臻g可展開(kāi)薄膜機(jī)構(gòu)樣機(jī)的收攏態(tài)及展開(kāi)態(tài)構(gòu)型如圖16所示。

    (a) 收攏態(tài)

    (b) 展開(kāi)態(tài)

    整個(gè)機(jī)構(gòu)在彈性伸桿的作用下沿著徑向方向同步展開(kāi),并通過(guò)繩索驅(qū)動(dòng)薄膜的展開(kāi),整個(gè)展開(kāi)過(guò)程流暢無(wú)卡滯。利用樣機(jī)進(jìn)行多次重復(fù)展收實(shí)驗(yàn),樣機(jī)均順利地實(shí)現(xiàn)了軸向的重復(fù)展收,且展收過(guò)程穩(wěn)定可靠,驗(yàn)證了所提出的構(gòu)型設(shè)計(jì)的合理性。樣機(jī)中的薄膜展開(kāi)過(guò)程如圖17所示,其中圖17(a)為Miura-ori薄膜的收攏態(tài)構(gòu)型,收攏態(tài)體積符合上述計(jì)算公式,所占空間小,滿足大折展比要求,且折痕總長(zhǎng)數(shù)值較小,驗(yàn)證了薄膜折痕參數(shù)優(yōu)化的有效性。根據(jù)Miura-ori的展開(kāi)方法,在繩索牽引下張拉3個(gè)角點(diǎn),實(shí)現(xiàn)薄膜的高效有序展開(kāi),得到圖17(d)所示的薄膜展開(kāi)態(tài)構(gòu)型,其展開(kāi)形態(tài)平整,具有較高的表面精度,驗(yàn)證了可展薄膜的Miura彈性折痕的可行性及優(yōu)越性。

    圖17 薄膜展開(kāi)過(guò)程

    5 結(jié) 論

    1)本文設(shè)計(jì)了一種空間薄膜可展機(jī)構(gòu),滿足薄膜結(jié)構(gòu)可折疊、輕量化、高折展比、高剛度和大尺寸等性能要求。建立幾何模型,計(jì)算了基于Miura-ori薄膜折疊方式的折展比和折痕總長(zhǎng),發(fā)現(xiàn)折展比和折痕總長(zhǎng)與折痕傾角和折痕間距有關(guān),并優(yōu)化得到最優(yōu)折痕參數(shù)。

    2)分別建立折痕傾角θ=90°及θ<90°的四折痕基本單元以及基于Miura彈性折痕的三角形薄膜數(shù)值仿真模型,發(fā)現(xiàn)其關(guān)鍵路徑上各節(jié)點(diǎn)應(yīng)力值均處于彈性范圍內(nèi),未發(fā)生塑性變形,驗(yàn)證了二維彈性折痕的可行性。

    3) 研制了空間薄膜可展機(jī)構(gòu)原理樣機(jī)并進(jìn)行了重復(fù)展收試驗(yàn),展開(kāi)過(guò)程流暢無(wú)卡滯,驗(yàn)證了機(jī)構(gòu)構(gòu)型設(shè)計(jì)的合理性,薄膜的大折展比和小折痕總長(zhǎng)表明了折痕參數(shù)優(yōu)化的有效性,薄膜展開(kāi)后表面平整說(shuō)明了二維彈性折痕的可行性及優(yōu)越性。

    通過(guò)對(duì)基于Miura彈性折痕的空間薄膜可展機(jī)構(gòu)的構(gòu)型設(shè)計(jì)、建模分析及樣機(jī)研制,為空間薄膜可展開(kāi)結(jié)構(gòu)的彈性折痕設(shè)計(jì)分析和進(jìn)一步研發(fā)提供了有益參考。后續(xù)將采用非接觸測(cè)量的方法來(lái)測(cè)試可展薄膜Miura彈性折痕的應(yīng)力分布情況,量化驗(yàn)證數(shù)值模擬工作的準(zhǔn)確性。

    猜你喜歡
    折痕構(gòu)型傾角
    《紡織品織物折痕回復(fù)角的測(cè)定》正式發(fā)布
    地球軸傾角的改斜歸正
    激光傾角儀在CT引導(dǎo)下經(jīng)皮肺穿刺活檢中的應(yīng)用
    車輪外傾角和前束角匹配研究
    北京汽車(2021年3期)2021-07-17 08:25:42
    分子和離子立體構(gòu)型的判定
    折痕
    青春(2017年5期)2017-05-22 11:57:33
    航天器受迫繞飛構(gòu)型設(shè)計(jì)與控制
    雙艙船
    遙感衛(wèi)星平臺(tái)與載荷一體化構(gòu)型
    午夜影院在线不卡| 日韩一区二区视频免费看| 只有这里有精品99| 日日啪夜夜爽| av在线播放精品| 有码 亚洲区| 久久久久国产网址| 赤兔流量卡办理| 热re99久久国产66热| 婷婷色麻豆天堂久久| 日本欧美视频一区| 狠狠精品人妻久久久久久综合| 男人爽女人下面视频在线观看| 日日摸夜夜添夜夜爱| 又大又黄又爽视频免费| 免费观看av网站的网址| 热99国产精品久久久久久7| 久久97久久精品| 日韩中文字幕视频在线看片| 亚洲激情五月婷婷啪啪| 高清黄色对白视频在线免费看| 高清视频免费观看一区二区| 午夜福利乱码中文字幕| 如何舔出高潮| 免费黄色在线免费观看| 一区二区av电影网| 最黄视频免费看| 日韩人妻精品一区2区三区| 亚洲性久久影院| 乱码一卡2卡4卡精品| 亚洲欧美日韩另类电影网站| av天堂久久9| 国产激情久久老熟女| 十分钟在线观看高清视频www| 中文欧美无线码| 一级,二级,三级黄色视频| 久久久久久久亚洲中文字幕| 亚洲国产成人一精品久久久| 水蜜桃什么品种好| 国产男女内射视频| 国产一区二区三区av在线| 亚洲中文av在线| 欧美变态另类bdsm刘玥| 国产高清三级在线| av线在线观看网站| a级毛色黄片| 在现免费观看毛片| 亚洲av免费高清在线观看| 美女福利国产在线| 国产1区2区3区精品| 日韩成人av中文字幕在线观看| 啦啦啦在线观看免费高清www| 日本av免费视频播放| 1024视频免费在线观看| 国产女主播在线喷水免费视频网站| 看十八女毛片水多多多| 另类精品久久| 欧美xxxx性猛交bbbb| 王馨瑶露胸无遮挡在线观看| 少妇猛男粗大的猛烈进出视频| 国产成人av激情在线播放| 激情五月婷婷亚洲| 成人国语在线视频| 亚洲欧美成人综合另类久久久| 亚洲精品乱码久久久久久按摩| 国产熟女欧美一区二区| 不卡视频在线观看欧美| 99re6热这里在线精品视频| 亚洲成av片中文字幕在线观看 | 超色免费av| 这个男人来自地球电影免费观看 | 亚洲,一卡二卡三卡| 欧美另类一区| 2018国产大陆天天弄谢| 十分钟在线观看高清视频www| 亚洲色图 男人天堂 中文字幕 | 国产精品成人在线| 日韩av不卡免费在线播放| 亚洲av欧美aⅴ国产| 黑丝袜美女国产一区| 日韩熟女老妇一区二区性免费视频| 国产欧美亚洲国产| 国产av一区二区精品久久| 美女中出高潮动态图| 最近手机中文字幕大全| 国产亚洲午夜精品一区二区久久| av免费在线看不卡| 国产极品粉嫩免费观看在线| 哪个播放器可以免费观看大片| 亚洲精品日韩在线中文字幕| 男女啪啪激烈高潮av片| freevideosex欧美| 欧美精品一区二区大全| 人人妻人人澡人人爽人人夜夜| 色网站视频免费| 久久久久国产网址| 女人久久www免费人成看片| 国产白丝娇喘喷水9色精品| 国产成人精品在线电影| 国产探花极品一区二区| 亚洲精品美女久久久久99蜜臀 | 999精品在线视频| 亚洲欧美中文字幕日韩二区| 咕卡用的链子| 亚洲av免费高清在线观看| 熟女人妻精品中文字幕| 侵犯人妻中文字幕一二三四区| 最近的中文字幕免费完整| a 毛片基地| 99热6这里只有精品| 亚洲欧美日韩另类电影网站| 久久精品久久久久久久性| 麻豆精品久久久久久蜜桃| 黄网站色视频无遮挡免费观看| 美女国产高潮福利片在线看| 亚洲国产色片| 国产精品嫩草影院av在线观看| 亚洲高清免费不卡视频| 春色校园在线视频观看| 丝袜美足系列| 一区二区三区乱码不卡18| 男男h啪啪无遮挡| 免费观看无遮挡的男女| 国产精品一区二区在线不卡| 亚洲内射少妇av| av又黄又爽大尺度在线免费看| 中国三级夫妇交换| 成人国语在线视频| 成人黄色视频免费在线看| 26uuu在线亚洲综合色| 国产麻豆69| 精品一区二区免费观看| 人成视频在线观看免费观看| 日本黄大片高清| 欧美xxⅹ黑人| 免费女性裸体啪啪无遮挡网站| 亚洲第一区二区三区不卡| 久久人妻熟女aⅴ| 1024视频免费在线观看| 精品国产露脸久久av麻豆| 肉色欧美久久久久久久蜜桃| 国产精品久久久久久久电影| 久久人妻熟女aⅴ| av.在线天堂| 亚洲国产精品一区二区三区在线| 午夜日本视频在线| a级毛片在线看网站| 纵有疾风起免费观看全集完整版| 国产男女超爽视频在线观看| 亚洲欧美精品自产自拍| 欧美日韩视频精品一区| 日韩三级伦理在线观看| 午夜福利视频精品| 高清不卡的av网站| 伦理电影免费视频| 一边亲一边摸免费视频| 多毛熟女@视频| 久久ye,这里只有精品| 久久午夜综合久久蜜桃| 久久 成人 亚洲| 日本欧美国产在线视频| av不卡在线播放| 色婷婷久久久亚洲欧美| a级毛片黄视频| 99热这里只有是精品在线观看| 久久久久久人人人人人| 日日啪夜夜爽| 午夜福利乱码中文字幕| 99国产精品免费福利视频| 中文字幕另类日韩欧美亚洲嫩草| 制服丝袜香蕉在线| 日本黄大片高清| 精品久久久久久电影网| 午夜久久久在线观看| 黑人欧美特级aaaaaa片| 久久99热这里只频精品6学生| 亚洲色图综合在线观看| 丰满少妇做爰视频| 精品少妇内射三级| 黄色视频在线播放观看不卡| 一级毛片我不卡| 亚洲四区av| 妹子高潮喷水视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩伦理黄色片| 在线观看三级黄色| 久久热在线av| 成年人午夜在线观看视频| 亚洲国产成人一精品久久久| 男人添女人高潮全过程视频| 黄色配什么色好看| 国产精品 国内视频| 国产国语露脸激情在线看| 少妇人妻 视频| 国产片内射在线| 精品卡一卡二卡四卡免费| 九色亚洲精品在线播放| 十八禁高潮呻吟视频| 久久av网站| 国产女主播在线喷水免费视频网站| 波多野结衣一区麻豆| 大片电影免费在线观看免费| 国产精品久久久久久av不卡| 国产精品国产三级国产专区5o| 亚洲少妇的诱惑av| 国产成人午夜福利电影在线观看| 精品国产一区二区三区久久久樱花| 街头女战士在线观看网站| 亚洲熟女精品中文字幕| 欧美日韩一区二区视频在线观看视频在线| 啦啦啦视频在线资源免费观看| 国产一区二区激情短视频 | 午夜久久久在线观看| 精品午夜福利在线看| 啦啦啦中文免费视频观看日本| av女优亚洲男人天堂| 亚洲精品456在线播放app| 国产深夜福利视频在线观看| 18禁国产床啪视频网站| 亚洲第一区二区三区不卡| 丝袜人妻中文字幕| 日韩伦理黄色片| 五月玫瑰六月丁香| 欧美人与性动交α欧美精品济南到 | 亚洲精品乱码久久久久久按摩| 人成视频在线观看免费观看| 日本-黄色视频高清免费观看| 久久久国产精品麻豆| 国产精品偷伦视频观看了| 99热网站在线观看| 久久影院123| 色婷婷久久久亚洲欧美| 丝瓜视频免费看黄片| 欧美变态另类bdsm刘玥| 黄色一级大片看看| 伊人亚洲综合成人网| 日日撸夜夜添| 丰满迷人的少妇在线观看| 女性被躁到高潮视频| 久久人人97超碰香蕉20202| 亚洲国产精品999| 亚洲第一区二区三区不卡| 亚洲国产av新网站| 黄色怎么调成土黄色| 日韩一区二区三区影片| 日韩视频在线欧美| 十八禁网站网址无遮挡| 国产女主播在线喷水免费视频网站| 日韩伦理黄色片| 激情五月婷婷亚洲| www日本在线高清视频| 自线自在国产av| 亚洲,欧美,日韩| 伊人亚洲综合成人网| 日韩大片免费观看网站| 亚洲 欧美一区二区三区| 美女国产高潮福利片在线看| 久久久久久人人人人人| 久久久欧美国产精品| 纵有疾风起免费观看全集完整版| 美女主播在线视频| 免费av中文字幕在线| av不卡在线播放| 乱人伦中国视频| 伦精品一区二区三区| 最近最新中文字幕大全免费视频 | 国产极品粉嫩免费观看在线| 日本猛色少妇xxxxx猛交久久| 亚洲国产最新在线播放| 内地一区二区视频在线| 午夜视频国产福利| 美女大奶头黄色视频| 少妇人妻 视频| 久久久久久久精品精品| 韩国av在线不卡| 亚洲精品一二三| 久久久久久久久久久久大奶| 亚洲精品色激情综合| 亚洲av电影在线观看一区二区三区| 97精品久久久久久久久久精品| 男女无遮挡免费网站观看| 人体艺术视频欧美日本| 国产麻豆69| 精品少妇内射三级| 国产日韩一区二区三区精品不卡| 2021少妇久久久久久久久久久| 免费黄网站久久成人精品| 色5月婷婷丁香| freevideosex欧美| 制服人妻中文乱码| 亚洲精品日本国产第一区| 日本91视频免费播放| 美女脱内裤让男人舔精品视频| 精品国产一区二区久久| 欧美日韩国产mv在线观看视频| 亚洲图色成人| 国产精品.久久久| 精品一区二区免费观看| 啦啦啦中文免费视频观看日本| 欧美成人精品欧美一级黄| 在线天堂中文资源库| 哪个播放器可以免费观看大片| 久久综合国产亚洲精品| 人妻系列 视频| 国产毛片在线视频| 欧美亚洲 丝袜 人妻 在线| 久久久久久久精品精品| 如何舔出高潮| 在线天堂中文资源库| 婷婷色综合大香蕉| 日本av手机在线免费观看| 黄色 视频免费看| 十八禁高潮呻吟视频| 国产免费现黄频在线看| 欧美bdsm另类| 久久97久久精品| 自线自在国产av| av有码第一页| 亚洲精品,欧美精品| 制服人妻中文乱码| 春色校园在线视频观看| 亚洲精品第二区| 免费在线观看完整版高清| 69精品国产乱码久久久| 少妇猛男粗大的猛烈进出视频| 亚洲图色成人| 99精国产麻豆久久婷婷| 一本—道久久a久久精品蜜桃钙片| 国产精品熟女久久久久浪| videos熟女内射| 精品福利永久在线观看| 国产精品三级大全| 一级黄片播放器| 亚洲在久久综合| 国产国语露脸激情在线看| av国产久精品久网站免费入址| 亚洲国产精品国产精品| 精品国产一区二区三区久久久樱花| 久久久久视频综合| 久久久久精品性色| 亚洲三级黄色毛片| 成年美女黄网站色视频大全免费| 日本猛色少妇xxxxx猛交久久| 亚洲欧美日韩另类电影网站| 亚洲国产精品999| 精品国产乱码久久久久久小说| 国语对白做爰xxxⅹ性视频网站| 精品一区二区三卡| 午夜福利网站1000一区二区三区| 国国产精品蜜臀av免费| 国产精品久久久久久精品古装| 伦理电影大哥的女人| 女人精品久久久久毛片| 精品一区二区三区视频在线| 乱码一卡2卡4卡精品| 91午夜精品亚洲一区二区三区| 青青草视频在线视频观看| 自线自在国产av| 国产女主播在线喷水免费视频网站| 男女国产视频网站| 青春草国产在线视频| 国产一区二区激情短视频 | 丝瓜视频免费看黄片| 亚洲国产精品一区二区三区在线| 少妇的逼好多水| 建设人人有责人人尽责人人享有的| 十八禁网站网址无遮挡| 亚洲国产av新网站| 99九九在线精品视频| 亚洲精品久久成人aⅴ小说| av在线观看视频网站免费| 丰满迷人的少妇在线观看| 18禁观看日本| 精品熟女少妇av免费看| 国产视频首页在线观看| 韩国高清视频一区二区三区| 五月开心婷婷网| √禁漫天堂资源中文www| 999精品在线视频| 国产福利在线免费观看视频| 久久影院123| 黑人巨大精品欧美一区二区蜜桃 | 一二三四中文在线观看免费高清| 精品亚洲乱码少妇综合久久| 欧美精品人与动牲交sv欧美| 人人妻人人添人人爽欧美一区卜| 高清在线视频一区二区三区| 亚洲av中文av极速乱| 王馨瑶露胸无遮挡在线观看| 亚洲,欧美精品.| 自拍欧美九色日韩亚洲蝌蚪91| 欧美亚洲日本最大视频资源| 老司机影院毛片| 丰满乱子伦码专区| 免费黄网站久久成人精品| 亚洲精品aⅴ在线观看| 成人漫画全彩无遮挡| 精品亚洲乱码少妇综合久久| 亚洲激情五月婷婷啪啪| 欧美性感艳星| 免费黄频网站在线观看国产| 国产黄色免费在线视频| 国产精品久久久久久av不卡| 亚洲国产看品久久| 制服诱惑二区| av电影中文网址| a级片在线免费高清观看视频| 国产麻豆69| 亚洲欧美日韩另类电影网站| 久久久久久久久久人人人人人人| 激情五月婷婷亚洲| 高清在线视频一区二区三区| 日产精品乱码卡一卡2卡三| 精品卡一卡二卡四卡免费| 午夜免费观看性视频| 最新中文字幕久久久久| 99国产精品免费福利视频| 菩萨蛮人人尽说江南好唐韦庄| 一二三四在线观看免费中文在 | 日本与韩国留学比较| 国内精品宾馆在线| 精品人妻熟女毛片av久久网站| 高清av免费在线| 香蕉精品网在线| 精品国产一区二区三区久久久樱花| 久久热在线av| 一级爰片在线观看| 精品国产一区二区三区久久久樱花| 成人影院久久| 少妇精品久久久久久久| 免费高清在线观看视频在线观看| 丰满少妇做爰视频| 国产成人午夜福利电影在线观看| 18禁裸乳无遮挡动漫免费视频| 成人二区视频| 国产免费一级a男人的天堂| 美女福利国产在线| 日韩成人伦理影院| 大码成人一级视频| 欧美成人精品欧美一级黄| 欧美精品一区二区大全| 91精品伊人久久大香线蕉| 亚洲人成网站在线观看播放| 日韩制服骚丝袜av| 一级黄片播放器| 久久精品国产自在天天线| 99热国产这里只有精品6| av女优亚洲男人天堂| 国产高清不卡午夜福利| 久久影院123| 巨乳人妻的诱惑在线观看| 国产一区二区在线观看日韩| 亚洲人与动物交配视频| 精品熟女少妇av免费看| 亚洲精品av麻豆狂野| 国产av一区二区精品久久| 亚洲一级一片aⅴ在线观看| 久热久热在线精品观看| 男女啪啪激烈高潮av片| 欧美日韩国产mv在线观看视频| 色5月婷婷丁香| 2018国产大陆天天弄谢| 丝瓜视频免费看黄片| 蜜臀久久99精品久久宅男| av不卡在线播放| 国产免费现黄频在线看| 国产精品女同一区二区软件| 天天影视国产精品| 亚洲欧洲国产日韩| 国产成人aa在线观看| 中文字幕人妻丝袜制服| 免费观看无遮挡的男女| 久久毛片免费看一区二区三区| 精品午夜福利在线看| 欧美成人午夜免费资源| 两性夫妻黄色片 | 涩涩av久久男人的天堂| 国产精品无大码| 亚洲精品第二区| 亚洲av欧美aⅴ国产| 精品卡一卡二卡四卡免费| 亚洲精品美女久久久久99蜜臀 | 日韩制服骚丝袜av| 国产无遮挡羞羞视频在线观看| 一级,二级,三级黄色视频| 国产精品久久久久成人av| 不卡视频在线观看欧美| 国产一区亚洲一区在线观看| 中文字幕免费在线视频6| 亚洲av中文av极速乱| 国内精品宾馆在线| 成人亚洲精品一区在线观看| 国产精品女同一区二区软件| 国产免费又黄又爽又色| 搡女人真爽免费视频火全软件| 国产69精品久久久久777片| 最近最新中文字幕大全免费视频 | 女性被躁到高潮视频| 欧美最新免费一区二区三区| 中文字幕av电影在线播放| 国产伦理片在线播放av一区| 国产激情久久老熟女| 久久久a久久爽久久v久久| 一区二区av电影网| 亚洲国产成人一精品久久久| 大香蕉久久成人网| 少妇 在线观看| 美女大奶头黄色视频| 高清欧美精品videossex| 国产男女超爽视频在线观看| 亚洲精品av麻豆狂野| 久久av网站| 久久久精品94久久精品| 精品一区二区三卡| 精品久久蜜臀av无| 在线观看人妻少妇| 国产成人91sexporn| 九九爱精品视频在线观看| 日本av免费视频播放| 亚洲五月色婷婷综合| 一级毛片黄色毛片免费观看视频| 丝袜美足系列| 少妇高潮的动态图| 久久精品久久精品一区二区三区| 高清视频免费观看一区二区| 欧美精品av麻豆av| 精品人妻在线不人妻| 久久久国产一区二区| 丁香六月天网| 久久久久久伊人网av| av黄色大香蕉| 亚洲国产av新网站| 老熟女久久久| 亚洲美女黄色视频免费看| 另类精品久久| 午夜日本视频在线| 最近的中文字幕免费完整| av网站免费在线观看视频| 久久久久网色| 久久99一区二区三区| 亚洲国产色片| 久久人妻熟女aⅴ| 春色校园在线视频观看| 精品国产露脸久久av麻豆| 黄色毛片三级朝国网站| 国产在线视频一区二区| 久久青草综合色| 黄色配什么色好看| 中文字幕av电影在线播放| 女人久久www免费人成看片| 亚洲 欧美一区二区三区| 亚洲精品美女久久av网站| 精品少妇久久久久久888优播| 十八禁网站网址无遮挡| 亚洲精品久久午夜乱码| 91精品三级在线观看| 综合色丁香网| 看免费成人av毛片| 日韩中文字幕视频在线看片| 国产日韩欧美在线精品| 中国三级夫妇交换| 又黄又爽又刺激的免费视频.| 中国三级夫妇交换| 亚洲国产毛片av蜜桃av| 亚洲精品乱久久久久久| 巨乳人妻的诱惑在线观看| 国产激情久久老熟女| 国产午夜精品一二区理论片| 久久久久网色| 国产免费福利视频在线观看| av不卡在线播放| 在线看a的网站| 国产 精品1| 18在线观看网站| 免费看av在线观看网站| 成人黄色视频免费在线看| 欧美日韩亚洲高清精品| 亚洲国产av影院在线观看| 高清在线视频一区二区三区| 99热网站在线观看| 欧美精品高潮呻吟av久久| 久久人人爽人人爽人人片va| 午夜视频国产福利| 中文字幕最新亚洲高清| 亚洲久久久国产精品| 久久97久久精品| 久久久久久久久久人人人人人人| 日韩免费高清中文字幕av| 成人毛片60女人毛片免费| 97精品久久久久久久久久精品| 国产一区二区三区av在线| 中文字幕免费在线视频6| 女人精品久久久久毛片| 国产免费又黄又爽又色| 97在线视频观看| 欧美+日韩+精品| 国产一级毛片在线| 视频区图区小说| 久久久国产欧美日韩av| 亚洲成国产人片在线观看| 日韩电影二区| 欧美国产精品一级二级三级| 高清视频免费观看一区二区| 国产免费一区二区三区四区乱码| 狠狠婷婷综合久久久久久88av| 老司机影院成人| 午夜日本视频在线| 飞空精品影院首页| 精品国产一区二区久久| 免费黄频网站在线观看国产| 色哟哟·www| 色5月婷婷丁香| 国产成人a∨麻豆精品| 18禁观看日本| 又大又黄又爽视频免费|