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

    雙層超彈性鉸鏈展開沖擊分析與優(yōu)化

    2016-06-23 09:27:30郭宏偉劉榮強(qiáng)陶建國
    振動(dòng)與沖擊 2016年9期
    關(guān)鍵詞:沖擊優(yōu)化

    楊 慧, 郭宏偉, 王 巖, 劉榮強(qiáng), 陶建國

    (哈爾濱工業(yè)大學(xué) 機(jī)電學(xué)院 機(jī)器人技術(shù)與系統(tǒng)國家實(shí)驗(yàn)室,哈爾濱 150080)

    雙層超彈性鉸鏈展開沖擊分析與優(yōu)化

    楊慧, 郭宏偉, 王巖, 劉榮強(qiáng), 陶建國

    (哈爾濱工業(yè)大學(xué) 機(jī)電學(xué)院 機(jī)器人技術(shù)與系統(tǒng)國家實(shí)驗(yàn)室,哈爾濱150080)

    摘要:超彈性鉸鏈依靠自身大撓度彈性折疊存儲的彈性勢能實(shí)現(xiàn)彈性展開,集驅(qū)動(dòng)、旋轉(zhuǎn)和鎖定于一體。研究了雙層超彈性鉸鏈展開沖擊性能,提出了降低雙層超彈性鉸鏈展開沖擊優(yōu)化方法。為了確保有限元模型的準(zhǔn)確性,搭建實(shí)驗(yàn)平臺對雙層超彈性鉸鏈進(jìn)行動(dòng)力學(xué)展開性能測試,基于實(shí)驗(yàn)結(jié)果對有限元模型進(jìn)行修正。采用正交法試驗(yàn)設(shè)計(jì)建立展開沖擊特性代理模型。以展開沖擊角度和鎖定時(shí)間最小為目標(biāo)函數(shù),采用改進(jìn)的非支配遺傳算法對雙層超彈性鉸鏈進(jìn)行參數(shù)化優(yōu)化設(shè)計(jì),得到最優(yōu)結(jié)構(gòu)尺寸。對最優(yōu)結(jié)構(gòu)尺寸進(jìn)行有限元展開分析,誤差不大于5.008%,表明代理模型準(zhǔn)確性。最優(yōu)非等厚度與等厚度雙層超彈性鉸鏈對比,發(fā)現(xiàn)前者沖擊角度和鎖定時(shí)間分別降低52.154%和29.104%,優(yōu)化后結(jié)構(gòu)的沖擊得到明顯改善。通過參數(shù)研究發(fā)現(xiàn),展開沖擊對外側(cè)帶簧厚度較為敏感,鎖定時(shí)間隨帶簧厚度增加表現(xiàn)出先減少后增大的趨勢。

    關(guān)鍵詞:沖擊;超彈性鉸鏈;動(dòng)力學(xué)分析;優(yōu)化;星載天線

    航天器受發(fā)射能力和體積的約束,如太陽能陣列、星載天線和固體反射器等受運(yùn)載器有效載荷和承載能力的限制。超彈性鉸鏈?zhǔn)且环N薄壁圓柱殼體結(jié)構(gòu),類似于帶彈簧(簡稱“帶簧”),依靠自身大撓度彈性變形實(shí)現(xiàn)折疊和彈性展開,集驅(qū)動(dòng)、旋轉(zhuǎn)和鎖定于一體,不含移動(dòng)部件、不存在摩擦,并具有輕質(zhì)特點(diǎn),能夠解決航天器發(fā)射體積有效性問題。并且超彈性鉸鏈能夠有效降低星載可展開機(jī)構(gòu)復(fù)雜度和發(fā)射成本,在航天領(lǐng)域具有極大的應(yīng)用前景。近年來出現(xiàn)了一些自鎖定帶簧超彈性鉸鏈形式[1]替代傳統(tǒng)鉸鏈,并成功應(yīng)用于波音回彈反射器[2]和火星快速航天器[3]。

    在航天器達(dá)到預(yù)定姿態(tài)時(shí),通過釋放約束使超彈性鉸鏈驅(qū)動(dòng)可展開機(jī)構(gòu)展開。Seffen等[4]分別基于能量法和沖量-動(dòng)量方程建立細(xì)長超彈性鉸鏈可移動(dòng)點(diǎn)鉸鏈展開動(dòng)力學(xué)模型,但未對展開之后超彈性鉸鏈的沖擊進(jìn)行分析。Soykasap[5]對三縫圓管超彈性鉸鏈的展開動(dòng)力學(xué)性能進(jìn)行分析,將展開過程分為展開、鎖定過程和鎖定三個(gè)階段。如果動(dòng)能過大超彈性鉸鏈將發(fā)生屈曲并繼續(xù)旋轉(zhuǎn),產(chǎn)生的振蕩會(huì)對其它結(jié)構(gòu)產(chǎn)生不利影響。Yao等[6]采用有限元方法對復(fù)合材料雙縫圓管超彈性鉸鏈展開過程中阻尼因素的影響進(jìn)行了研究。Mallikarachchi等[7]通過實(shí)驗(yàn)研究了雙縫圓管超彈性鉸鏈展開過程的沖擊,并對由兩個(gè)縱向縫超彈性鉸鏈構(gòu)成1 m長輕質(zhì)自展開圓桿進(jìn)行展開動(dòng)力學(xué)分析,通過實(shí)驗(yàn)驗(yàn)證桿輕質(zhì)桿繞一個(gè)小型航天器兩次折疊和展開過程中不會(huì)出現(xiàn)損壞。

    本文研究的雙層超彈性鉸鏈由彈性金屬合金構(gòu)成,剛度較大,在展開末端時(shí)殼體突然翻轉(zhuǎn)可能對可展開機(jī)構(gòu)產(chǎn)生過沖,會(huì)對其他結(jié)構(gòu)的穩(wěn)定性產(chǎn)生影響。為了降低雙層超彈性鉸鏈展開末端沖擊,研究了雙層超彈性鉸鏈動(dòng)力學(xué)展開沖擊性能,基于代理模型法提出了一種快速有效的降低雙層超彈性鉸鏈展開沖擊的優(yōu)化方法。搭建動(dòng)力學(xué)展開實(shí)驗(yàn)平臺對雙層超彈性鉸鏈進(jìn)行展開性能測試,基于實(shí)驗(yàn)結(jié)果對有限元模型進(jìn)行修正。由于雙層超彈性鉸鏈折疊和展開具有高度非線性,采用兩因子五水平全因子正交法進(jìn)行試驗(yàn)設(shè)計(jì),借助于ABAQUS/Explicit軟件進(jìn)行數(shù)值仿真建立25組實(shí)驗(yàn)樣本點(diǎn),利用四次多項(xiàng)式建立表征雙層超彈性鉸鏈展開沖擊特性的響應(yīng)面數(shù)學(xué)模型。以展開沖擊角度和鎖定時(shí)間最小為目標(biāo)函數(shù),采用改進(jìn)的非支配遺傳算法對雙層超彈性鉸鏈進(jìn)行參數(shù)化優(yōu)化設(shè)計(jì),得到最優(yōu)結(jié)構(gòu)尺寸?;趯?shí)驗(yàn)設(shè)計(jì)的樣本點(diǎn)進(jìn)行參數(shù)化研究,分析位于折疊外側(cè)的兩層帶簧厚度對雙層超彈性鉸鏈展開沖擊的影響。

    1超彈性鉸鏈大撓度折展特性

    1.1靜力學(xué)彎曲特性

    超彈性鉸鏈?zhǔn)且活愌乜v向開縫薄壁圓柱殼體結(jié)構(gòu),依靠自身大撓度彈性變形存儲彈性勢能實(shí)現(xiàn)彈性展開,折疊和展開過沖中力矩具有高度非線性,如圖1所示。定義帶簧彎曲方向和橫截面曲率方向相反時(shí)為背向彎曲,相同時(shí)為對向彎曲。圖 1中實(shí)心箭頭表示對向和背向彎曲折疊力矩路徑,空心箭頭表示相應(yīng)的展開路徑。帶簧超彈性鉸鏈對向折展時(shí)力矩路徑相同,但是背向折展時(shí)的力矩路徑不同,展開時(shí)的峰值力矩遠(yuǎn)小于折疊時(shí)的峰值力矩。

    圖1 帶簧超彈性鉸鏈背向和對向彎曲力矩曲線Fig.1 Moment curves for opposite and equal bend of tape-spring

    1.2動(dòng)力學(xué)展開特性

    不同超彈性鉸鏈在對向和背向折疊過程中力矩因結(jié)構(gòu)屈曲出現(xiàn)峰值,隨著折疊角度增加力矩趨于穩(wěn)定,超彈性鉸鏈動(dòng)力學(xué)展開過程恢復(fù)曲線如圖 2所示。展開時(shí),角度隨著時(shí)間增加而降低,該過程中力矩保持穩(wěn)定。在展開接近平衡位置附近發(fā)生突然翻轉(zhuǎn),展開力矩出現(xiàn)峰值能夠?qū)φ麄€(gè)機(jī)構(gòu)具有鎖定作用,超彈性鉸鏈嘗試鎖定,在t0時(shí)刻實(shí)現(xiàn)鎖定。但是當(dāng)展開峰值力矩過大時(shí),展開過程中會(huì)出現(xiàn)一定量沖擊角θ,會(huì)對可展開機(jī)構(gòu)中的其它部件產(chǎn)生影響。

    圖2 超彈性鉸鏈展開恢復(fù)曲線Fig.2 Deployment angle-time for double-layer tape-spring hinge

    雙層超彈性鉸鏈由彎曲面相對的雙層帶簧疊加而成,一端固定,另一端固連長度為l1的剛性桿構(gòu)成折展機(jī)構(gòu),如圖3所示,按照圖示方向折疊時(shí)內(nèi)側(cè)、外側(cè)兩層帶簧將分別發(fā)生背向、對向彎曲,展開時(shí)內(nèi)側(cè)、外側(cè)兩層帶簧分別發(fā)生對向、背向彎曲。雙層超彈性鉸鏈折疊時(shí)帶簧彎曲變形存儲彈性勢能,展開時(shí)彈性勢能轉(zhuǎn)換為剛性桿和鉸鏈的動(dòng)能,在展開接近直線位置時(shí)若動(dòng)能大于外側(cè)帶簧產(chǎn)生突然翻轉(zhuǎn)所需要的能量,則外側(cè)帶簧將會(huì)發(fā)生背向彎曲。若能量不變,則帶簧變形量隨厚度增加而降低。圖3中內(nèi)側(cè)兩層帶簧厚度均為w0,外側(cè)兩層帶簧厚度分別為w1和w2。

    圖3 含超彈性鉸鏈的剛性桿幾何示意圖Fig.3 Geometry for the double-layer tape-spring hinge

    在之前的文獻(xiàn)[8-9]已經(jīng)對影響沖擊性能的其它因素進(jìn)行了優(yōu)化,如帶簧橫截面半徑、中心角、分離距離等因素。根據(jù)Seffen對超彈性鉸鏈折展時(shí)能量的分析,并考慮材料阻尼對能量耗散和固連剛性桿慣性力影響,通過改變外側(cè)帶簧厚度w1和w2研究雙層超彈性鉸鏈展開沖擊特性。

    2展開實(shí)驗(yàn)研究

    圖4為雙層超彈性鉸鏈展開實(shí)驗(yàn)裝置示意圖。測試的超彈性鉸鏈橫截面半徑為17.8 mm,中心角為76°,橫截面間距為16 mm,厚度均為0.12 mm。在精密光學(xué)平臺上雙層超彈性鉸鏈一端固定,另一端連接著外徑為22 mm和厚度為2.5 mm不銹鋼剛性圓桿,在剛性圓桿中央連接一個(gè)球形滾輪,用于補(bǔ)償重力對整個(gè)機(jī)構(gòu)展開沖擊的影響。

    圖4 展開實(shí)驗(yàn)裝置示意圖Fig.4 Dynamic deploying experimental setup

    美國VRI公司生產(chǎn)的Phantom V12.1高速相機(jī)固定在雙層超彈性鉸鏈折疊區(qū)域正上方,用以記錄展開過程中雙層超彈性鉸鏈變形位置和角度。雙層超彈性鉸鏈初始折疊角度為150°,高速相機(jī)采集速率設(shè)定為600 幀/秒。圖5為雙層超彈性鉸鏈動(dòng)力學(xué)展開過程中的位置圖。

    圖5 雙層超彈性鉸鏈初始折疊150°時(shí)展開過程圖Fig.5 Photos taken during deployment of double-layer tape-spring hinge folded 150 deg

    3雙層超彈性鉸鏈動(dòng)力學(xué)展開數(shù)值仿真

    為了對超彈性鉸鏈進(jìn)行展開沖擊分析,必須先把其折疊再進(jìn)行展開動(dòng)力學(xué)分析,而超彈性鉸鏈折疊和展開過程是高度非線性的過程,分析耗時(shí)較長。為此,先采用仿真方法建立實(shí)驗(yàn)樣本點(diǎn),利用多項(xiàng)式響應(yīng)面法建立雙層超彈性鉸鏈展開沖擊特性的代理模型,再采用改進(jìn)的非支配遺傳算法進(jìn)行多目標(biāo)優(yōu)化分析。

    3.1雙層超彈性鉸鏈折展有限元模型

    由于薄壁圓柱殼體類帶簧結(jié)構(gòu)大撓度彎曲是高度非線性過程,采用隱式求解器求解,將會(huì)遇到不收斂、數(shù)值穩(wěn)定和奇異性等問題,而顯示求解器可以避免這類問題的出現(xiàn)。所以,在ABAQUS中采用4節(jié)點(diǎn)縮減積分單元(S4R)建立有限元模型,如圖 6所示,利用顯示求解器(Explicit)進(jìn)行求解。

    圖6 雙層超彈性鉸鏈有限元模型及邊界條件Fig.6 Finte element model and boundaries for the DLTS hinge

    為了簡化模型,在雙層超彈性鉸鏈兩端分別建立參考點(diǎn),按照轉(zhuǎn)動(dòng)慣量相等的原則把夾持端質(zhì)量等效到參考點(diǎn)上,通過施加運(yùn)動(dòng)耦合約束把參考點(diǎn)與周圍區(qū)域連接,模擬夾持端進(jìn)行加載。參考點(diǎn)2(RP-2)釋放沿z軸移動(dòng)的自由度和饒y軸旋轉(zhuǎn)自由度,繞y軸旋轉(zhuǎn)角度為166°。在超彈性鉸鏈固定端通過shell to solid鏈接一個(gè)尺寸為30 mm×30 mm×12 mm固體,在靠近超彈性鉸鏈通過網(wǎng)格劃分出3層固體單元(C3D8),在另一端布置一層無限單元(CIN3D8),以此模擬實(shí)驗(yàn)中固定端中阻尼耗散對動(dòng)力學(xué)展開能量的耗散作用。在分析中添加自接觸(self-contact)以模擬折展過程中帶簧之間的接觸[10],同時(shí),為了提高計(jì)算效率,摩擦因數(shù)設(shè)置為零。

    在ABAQUS中使用結(jié)構(gòu)阻尼假設(shè)的動(dòng)力學(xué)分析包括穩(wěn)態(tài)響應(yīng)分析和隨機(jī)響應(yīng)分析,瞬態(tài)動(dòng)力學(xué)分析不能直接使用結(jié)構(gòu)阻尼。而且,系統(tǒng)結(jié)構(gòu)阻尼特性與結(jié)構(gòu)內(nèi)摩擦機(jī)理相關(guān),在雙層超彈性鉸鏈有限元模型中帶簧之間設(shè)置為無摩擦接觸,整個(gè)模型中未考慮摩擦。鉸鏈展開動(dòng)力學(xué)沖擊分析屬于瞬態(tài)動(dòng)力學(xué)分析,所以在沖擊分析中與結(jié)構(gòu)阻尼相關(guān)量均設(shè)為零,不會(huì)對結(jié)論產(chǎn)生影響。

    3.2材料特性

    雙層超彈性鉸鏈帶簧材料為Ni36CrTiAl,是奧氏體高彈性合金,具有較高的強(qiáng)度和彈性模量、良好的耐腐蝕性。通過拉伸試驗(yàn)[8]測量出材料Ni36CrTiAl性能參數(shù),見表1。

    表1 超彈性鉸鏈Ni36CrTiAl材料參數(shù)

    3.3仿真模型驗(yàn)證

    雙層超彈性鉸鏈準(zhǔn)靜態(tài)折展有限元模型在之前研究中已通過實(shí)驗(yàn)驗(yàn)證了準(zhǔn)確性[8]。為了保證雙層超彈性鉸鏈展開時(shí)動(dòng)能從零狀態(tài)開始,在準(zhǔn)靜態(tài)折疊步之后增加能量耗散步,在展開步中的分析時(shí)間具有真實(shí)意義。仿真模型尺寸與實(shí)驗(yàn)測試件保持一致。圖7為修正的有限元模型仿真值與實(shí)驗(yàn)值對比曲線。仿真與實(shí)驗(yàn)時(shí)超彈性鉸鏈?zhǔn)状蔚竭_(dá)平衡位置的時(shí)間分別是0.275 s、0.278 s,過沖角度分別是82.21°、77.96°,完全鎖定時(shí)間分別是0.752 s、 0.735 s,與實(shí)驗(yàn)值的相對誤差分別為1.08%、-5.45%、-2.26%,表明修正后模型的準(zhǔn)確性。

    圖7 仿真與實(shí)驗(yàn)恢復(fù)曲線對比Fig.7 Compared between simulation and experiment for rotation-time curves

    4雙層超彈性鉸鏈展開沖擊特性代理模型

    響應(yīng)面法(RSM)是統(tǒng)計(jì)方法和數(shù)學(xué)方法的結(jié)合,通過實(shí)驗(yàn)設(shè)計(jì)選取設(shè)計(jì)變量進(jìn)行有限次試驗(yàn)得到樣本點(diǎn),然后對樣本點(diǎn)進(jìn)行數(shù)據(jù)擬合。通過響應(yīng)面法建立代理模型,并對代理模型進(jìn)行優(yōu)化具有高效性,并且便于實(shí)現(xiàn),在眾多優(yōu)化問題中得到應(yīng)用[11-12]。

    該研究中雙層超彈性鉸鏈展開沖擊特性參數(shù)為最大應(yīng)力Sf、沖擊角θ和鎖定時(shí)間t0,這些特性參數(shù)可以用如下的基函數(shù)表示:

    (1)

    (2)

    式中,β0、βij、βiij、βiiij、βiijj分別為相應(yīng)階次基函數(shù)系數(shù),i,j(j=1, 2, …,n)是第i,j個(gè)獨(dú)立變量。

    利用最小二乘法可以計(jì)算四次多項(xiàng)式函數(shù)的系數(shù),

    bf=(ΦTΦ)-1ΦTy

    (3)

    式中,bf是四次多項(xiàng)式函數(shù)的系數(shù)bf= (β1,β2, …,βN),N是基函數(shù)φi序號。矩陣Φ由如下M個(gè)樣本點(diǎn)構(gòu)成,

    (4)

    式中,φN(x)M是第M個(gè)樣本點(diǎn)的第N個(gè)基函數(shù)。

    在以下分析中利用之前對雙層超彈性鉸鏈結(jié)構(gòu)優(yōu)化[8]的結(jié)果,雙層超彈性鉸鏈基本幾何尺寸為:帶簧橫截面半徑R=17.046 mm,橫截面中心角φ=84.698 °,帶簧分離距離s=19.844 mm,縱向長度L=126 mm,內(nèi)側(cè)兩層帶簧厚度固定為w0=0.12 mm。與雙層超彈性鉸鏈固連剛性桿的長度、外徑和厚度分別為l1=165 mm、D=22 mm和t=2.5 mm,材料為鋁合金,彈性模量E=70 GPa,泊松比σ=0.33,密度ρ=2.7×103kg/m3。

    為了建立雙層超彈性鉸鏈展開沖擊性能參數(shù)的代理模型,需要在設(shè)計(jì)空間中有足夠的樣本點(diǎn)?;诙蛩匚逅秸粚?shí)驗(yàn)設(shè)計(jì)法對樣本點(diǎn)進(jìn)行設(shè)計(jì),如表2所示。利用25組樣本點(diǎn)將式(4)代入到式(3)可以計(jì)算出多項(xiàng)式函數(shù)系數(shù),建立雙層超彈性鉸鏈展開沖擊特性最大應(yīng)力Sf、沖擊角度θ和鎖定時(shí)間t0關(guān)于外側(cè)的內(nèi)、外兩層帶簧厚度w1和w2的代理模型如下

    (5)

    表2 25組設(shè)計(jì)樣本點(diǎn)的有限元結(jié)果

    θ=-50 806.9+(1.46×106w2+

    (6)

    (7)

    (8)

    (9)

    (10)

    (11)

    (12)

    表3 設(shè)計(jì)樣本點(diǎn)代理模型的精度評價(jià)

    分析發(fā)現(xiàn)相關(guān)系數(shù)和復(fù)相關(guān)系數(shù)均大于0.98,相對誤差不大于7.5%,表明代理模型具有足夠的精度。

    5超彈性鉸鏈優(yōu)化設(shè)計(jì)

    雙層超彈性鉸鏈展開末端存在的沖擊和鎖定時(shí)間的長度會(huì)對其他結(jié)構(gòu)產(chǎn)生影響,而折疊時(shí)應(yīng)力集中現(xiàn)象會(huì)影響其使用性能和次數(shù)。超彈性鉸鏈的帶簧背向彎曲時(shí)折展過程力矩路徑不同,與折展屈曲和突然翻轉(zhuǎn)對應(yīng)的彈性勢能不同,帶簧的厚度變化也改變屈曲和突然翻轉(zhuǎn)的彈性勢能。所以,選取雙層超彈性鉸鏈展開沖擊角θ和鎖定時(shí)間t0為目標(biāo)函數(shù),折疊最大應(yīng)力Sf為約束變量,外層的內(nèi)、外層帶簧厚度w1、w2為自變量,建立如下優(yōu)化模型:

    (13)

    目標(biāo)函數(shù)fo,是各個(gè)目標(biāo)分量fk及相應(yīng)的權(quán)重因子wk與比例因子sk的和:

    (14)

    式中,k(k=1, 2, …,q)為目標(biāo)函數(shù)序號,q(q=2)是目標(biāo)函數(shù)分量的總數(shù)。如果對不同數(shù)量級的目標(biāo)函數(shù)分量均選用相同的權(quán)重因子和比例因子,那么數(shù)量級較小對目標(biāo)函數(shù)的影響會(huì)被削弱。所以,最大應(yīng)力Sf、沖擊角θ和鎖定時(shí)間t0的比例因子分別選取為s1=0.5、s2=40.0、s3=0.25,權(quán)重因子分別為w1=w2=w3=1.0。在iSIGHT中采用改進(jìn)的非支配遺傳算法[13]對式(13)中模型進(jìn)行多目標(biāo)優(yōu)化設(shè)計(jì),設(shè)置種群代數(shù)為100,每代種群數(shù)120。圖 8為雙層超彈性鉸鏈展開沖擊優(yōu)化的Pareto。

    圖8 Pareto最優(yōu)解Fig.8 Pareto optimal results

    由圖8可知沖擊角θ和鎖定時(shí)間t0呈現(xiàn)負(fù)相關(guān)特性,沖擊角θ增大時(shí)鎖定時(shí)間t0變短。最優(yōu)解為w1=0.138 mm、w2=0.129 mm,并計(jì)算該最優(yōu)結(jié)構(gòu)的有限元解與代理模型解、優(yōu)化模型約束條件進(jìn)行對比見表4。最優(yōu)結(jié)構(gòu)與代理模型的相對誤差不大于5.008%,表明所選取的最優(yōu)結(jié)構(gòu)滿足優(yōu)化模型約束條件,且驗(yàn)證了該最優(yōu)解的可用性和代理模型的準(zhǔn)確性。

    表4 最優(yōu)解的代理模型和有限元模型結(jié)果對比

    表4中列出最優(yōu)解與等厚度雙層超彈性鉸鏈(表2中第13號樣本點(diǎn))有限元仿真得到的展開性能參數(shù)見表 5,等厚度結(jié)構(gòu)和最優(yōu)非等厚結(jié)構(gòu)的展開角度-時(shí)間曲線進(jìn)行對比如圖9所示。

    由表5和圖9對比可知,最優(yōu)非等厚結(jié)構(gòu)使得雙層超彈性鉸鏈展開沖擊角度降低52.154%,鎖定時(shí)間縮短了29.104%,明顯提高了雙層超彈性鉸鏈結(jié)構(gòu)的展開性能。

    表5 優(yōu)化結(jié)構(gòu)和初始等厚度結(jié)構(gòu)展開性能參數(shù)對比

    圖9 優(yōu)化結(jié)構(gòu)和初始等厚度結(jié)構(gòu)展開角度-時(shí)間對比Fig.9 Compared deployment angle-time curves between original and optimal configurations

    6參數(shù)研究

    根據(jù)表 2研究內(nèi)層、外層帶簧厚度w1、w2對雙層超彈性鉸鏈折疊最大應(yīng)力Sf、展開沖擊角θ和鎖定時(shí)間t0的影響。進(jìn)行參數(shù)研究的雙層超彈性鉸鏈基本幾何尺寸為:帶簧橫截面半徑R=17.046 mm,橫截面中心角=84.698°,帶簧分離距離s=19.844 mm,縱向長度L=126 mm,鉸鏈折疊時(shí)發(fā)生背向彎曲的內(nèi)側(cè)帶簧厚度均固定為w0=0.12 mm,剛性桿長度和外徑分別為l1=165 mm和D=20 mm。超彈性鉸鏈外側(cè)兩層帶簧的厚度w1和w2變化范圍均是從0.10 mm~0.14 mm。

    圖10和圖 11分別為內(nèi)層和外層帶簧厚度w1、w2對超彈性鉸鏈折疊最大應(yīng)力Sf的影響曲線??梢钥闯觯S著w1和w2的增加,Sf增加范圍分別為6.333%~12.906%和12.044%~18.969%。

    圖10 外側(cè)內(nèi)層帶簧厚度對最大應(yīng)力的影響Fig.10 Effect of inner layer thickness on the maximum stress

    圖11 外側(cè)外層帶簧厚度對最大應(yīng)力的影響Fig.11 Effect of outer layer thickness on the maximum stress

    圖12和圖13分別為w1、w2對超彈性鉸鏈展開沖擊角θ的影響曲線。從圖中可以看出,θ隨著w1和w2的增加而降低的范圍分別為47.211%~66.925%和50.154%~68.770=%。

    圖12 外側(cè)內(nèi)層帶簧厚度對沖擊角度的影響Fig.12 Effect of inner layer thickness on the shock angle

    圖13 外側(cè)外層帶簧厚度對沖擊角度的影響Fig.13 Effect of inner layer thickness on the shock angle

    圖14和圖15分別為w1、w2對超彈性鉸鏈展開鎖定時(shí)間t0的影響曲線。從圖中可以看出,t0隨著w1的增加先降低后增加,t0先降低的范圍為26.866%~35.484%,之后t0增加的范圍為3.571%~7.843%;t0隨著w2的增加先降低后增加,t0先降低的范圍為17.460%~27.957%,t0后增加的范圍為3.175%~5.769%。

    圖14 外層內(nèi)層帶簧厚度對鎖定時(shí)間的影響Fig.14 Effect of outer layer thickness on the locked time

    圖15 外側(cè)外層帶簧厚度對鎖定時(shí)間的影響Fig.15 Effect of outer layer thickness on the locked time

    通過對雙層超彈性鉸鏈展開沖擊特性進(jìn)行參數(shù)研究發(fā)現(xiàn),外側(cè)兩層帶簧片厚度w1、w2對其折疊最大應(yīng)力Sf、展開沖擊角θ和鎖定時(shí)間t0的影響都很顯著:隨著厚度w1、w2增加,Sf增大,θ降低,而t0則是先降低后增大。主要是由于折疊時(shí)內(nèi)側(cè)產(chǎn)生背向彎曲的帶簧厚度不變,相當(dāng)于折疊存貯的彈性勢能不變,而外側(cè)帶簧厚度越大在展開發(fā)生突然翻轉(zhuǎn)時(shí)產(chǎn)生相同變形需要的能量越大,所以沖擊角度隨著外側(cè)帶簧厚度增大而降低。但是外側(cè)帶簧厚度過大時(shí)產(chǎn)生突然翻轉(zhuǎn)產(chǎn)生相同變形量存儲過多的彈性勢能,會(huì)導(dǎo)致雙層超彈性鉸鏈二次沖擊角度增大,進(jìn)而使其鎖定時(shí)間增加。

    7結(jié)論

    基于ABAQUS/Explicit有限元軟件,以降低雙層超彈性鉸鏈展開沖擊角度和縮短鎖定時(shí)間為目標(biāo),以外側(cè)兩層帶簧厚度為設(shè)計(jì)變量進(jìn)行優(yōu)化設(shè)計(jì),并對展開沖擊性能進(jìn)行參數(shù)研究,得到如下結(jié)論:

    (1) 搭建實(shí)驗(yàn)平臺對雙層超彈性鉸鏈進(jìn)行動(dòng)力學(xué)展開性能測試,基于實(shí)驗(yàn)結(jié)果對有限元模型進(jìn)行了修正,保證了仿真分析的準(zhǔn)確性。

    (2) 基于有限元仿真結(jié)果,利用四次多項(xiàng)式響應(yīng)面法建立雙層超彈性鉸鏈展開沖擊性能代理模型,誤差不大于7.543%,表明代理模型能夠精確表征雙層超彈性鉸鏈外側(cè)帶簧厚度和展開沖擊性能之間的關(guān)系。

    (3) 以雙層超彈性鉸鏈展開沖擊角度最小和鎖定時(shí)間最短為目標(biāo)函數(shù),以外側(cè)兩層帶簧厚度為設(shè)計(jì)變量進(jìn)行展開沖擊性能優(yōu)化,得到了最優(yōu)非等厚雙層超彈性鉸鏈結(jié)構(gòu)參數(shù),與等厚度雙層超彈性鉸鏈相比最優(yōu)結(jié)構(gòu)展開沖擊角度減小了52.154%、鎖定時(shí)間縮短了29.104%。

    (4) 基于正交試驗(yàn)設(shè)計(jì)對雙層超彈性鉸鏈展開沖擊特性進(jìn)行參數(shù)研究,發(fā)現(xiàn)外側(cè)兩層帶簧片厚度增加時(shí),最大應(yīng)力增大、展開沖擊角θ降低,而鎖定時(shí)間t0則是先降低后增大。

    參 考 文 獻(xiàn)

    [1] Lane S A,Murphey T W,Zatman M. Overview of the innovative space-based radar antenna technology program[J]. Journal of Spacecraft and Rockets,2011,48(1):135-145.

    [2] Anon. Hughes graphite antennas installed on Mast-2 craft[N]. Space News,Nov. 1994.

    [3] Adams D S,Mobrem M. Lenticular jointed antenna deployment anomaly and resolution onboard the Mars Express Spacecraft[J]. Journal of Spacecraft and Rockets,2009,46(2):403-410.

    [4] Seffen K A,Pellegrino S. Deployment dynamics of tape springs[J]. Proceedings of the Royal Society of London,Series A:Mathematical and Physical Sciences,1999,455(1983):1003-1048.

    [5] Soykasap O. Deployment analysis of a self-deployable composite boom[J]. Composite Structures,2009,89:374-381.

    [6] Yao X F,Ma Y J,Yin Y J,et al. Design theory and dynaimic mechanical characterization of the deployable composite tube hinge[J]. Science China:Physics,Mechanics and Astronomy,2011,54(4):633-639.

    [7] Mallikarachchi H M Y C,Pellegrino S. Deployment dynamics of ultrathin composite booms with tape-spring hinges[J]. Journal of Spacecraft and Rockets,2014,51(2):604-613.

    [8] Yang H,Liu R Q,Wang Y,et al. Experiment and multiobjective optimization design of tape-spring hinges[J]. Structural and Multidisciplinary Optimization,2015,51(6):1373-1384.

    [9] Yang H,Deng Z Q,Liu R Q,et al. Optimizing the quasi-static folding and deploying of thin-walled tube flexure hinges with double slots[J]. Chinese Journal of Mechanical Engineering,2014,27(2):279-286.

    [10] Mallikarachchi H M Y C,Pellegrino S. Quasi-static folding and deployment of ultrathin composite tape-spring hinges[J]. Journal of Spacecraft and Rockets,2011,48(1):187-198.

    [11] Hou S J,Liu T Y,Dong D,et al. Factor screening and multivariable crashworthiness optimization for vehicle side impact by factorial design[J]. Structural and Multidisciplinary Optimization,2014,49:147-167.

    [12] Li M,Deng Z Q,Guo H W,et al. Crashworthiness analysis on alternative square honeycomb structure under axial loading[J]. Chinese Journal of Mechanical Engineering,2013,26(4):784-792.

    [13] Deb K,Pratap A,Agarwal S,et al. A fast and elitist multiobjective genetic algorithm:NSGA-Ⅱ[J]. IEEE Transactions on Evolutionary Computation,2002,6(2):182-197.

    Dynamic analysis and optimization for double-layer tape-spring hinges unfolding

    YANG Hui, GUO Hong-wei, WANG Yan, LIU Rong-qiang, TAO Jian-guo

    (State Key Laboratory of Robotics and System, School of Mechatronics Engineering, Harbin Institute of Technology, Harbin 150080, China)

    Abstract:Flexible tape-spring hinges can be folded elastically and can self-unfold by releasing stored strain energy with fewer component parts and slight weights. The unfolding dynamic behaviors of double-layer tape-spring (DLTS) hinges were studied and optimized to reduce unfolding shock. A platform for unfolding test of DLTS hinges was built and tests were performed to verify the correctness of the FE numerical model for double-layer tape-spring hinges. The unfolding shock angle and the locked time were taken as objectives, the thicknesses of outer layer tape-spring and inner layer one in outer bending direction of a DLTS hinge were taken as design variables. The modified non-dominated sorting genetic algorithm was employed to conduct the optimal design of DLTS hinges. The finite element models for the non-equal thickness optimal design based on numerical methods were established and validated. The optimization results had relative errors no more than 5.008%. Comparing unfolding behaviors of the non-equal thickness optimal design and those of the equal thickness one for DLTS hinges, it was shown that the unfolding shock angle and locked time for the former reduce 52.154% and 29.104%, respectively; after optimization, the shock of the structure is reduced obviously. The parameter study showed that the unfolding shock angle is more sensitive to the outer layer thickness of tape-spring, and the locked time decreases firstly and then increases with increase in the thickness of tape-spring.

    Key words:unfolding shock; double-layer tape-spring hinge; dynamic analysis; optimization; spaceborne antenna

    基金項(xiàng)目:國家高校111計(jì)劃(B07018);機(jī)器人技術(shù)與系統(tǒng)國家實(shí)驗(yàn)室自主課題(SKLRS201401A02);中央高?;A(chǔ)研究基金(HIT.NSRIF. 2015050)

    收稿日期:2015-01-27修改稿收到日期:2015-05-22

    通信作者劉榮強(qiáng) 男,博士,教授,1965年生

    中圖分類號:V214.3+6

    文獻(xiàn)標(biāo)志碼:A

    DOI:10.13465/j.cnki.jvs.2016.09.004

    第一作者 楊慧 女,博士生,1986年6月生

    猜你喜歡
    沖擊優(yōu)化
    超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
    民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
    關(guān)于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
    一種新型沖擊器拆卸及預(yù)緊裝置
    奧迪Q5換擋沖擊
    奧迪A8L換擋沖擊
    基于低碳物流的公路運(yùn)輸優(yōu)化
    一汽奔騰CA7165AT4尊貴型車換擋沖擊
    无遮挡黄片免费观看| 久久亚洲精品不卡| 国产精品久久久久久精品电影| 天堂动漫精品| 一本一本综合久久| 50天的宝宝边吃奶边哭怎么回事| 一二三四在线观看免费中文在| 女人高潮潮喷娇喘18禁视频| АⅤ资源中文在线天堂| 欧美日韩一级在线毛片| 国产成人av激情在线播放| 国产精品1区2区在线观看.| 色在线成人网| 亚洲精品美女久久av网站| 久久精品91无色码中文字幕| 色尼玛亚洲综合影院| 天天一区二区日本电影三级| 国产精品av视频在线免费观看| 九九热线精品视视频播放| 午夜免费激情av| √禁漫天堂资源中文www| 欧美日韩国产亚洲二区| 在线视频色国产色| 亚洲av日韩精品久久久久久密| 国模一区二区三区四区视频 | 国产伦一二天堂av在线观看| 亚洲欧美日韩无卡精品| 日日干狠狠操夜夜爽| 国产69精品久久久久777片 | 在线视频色国产色| 成人午夜高清在线视频| 老司机福利观看| 国产伦一二天堂av在线观看| 色噜噜av男人的天堂激情| 精品国产乱码久久久久久男人| 亚洲色图 男人天堂 中文字幕| 国产人伦9x9x在线观看| 亚洲精品在线美女| 亚洲精品色激情综合| 麻豆av在线久日| 精品久久久久久成人av| 一级黄色大片毛片| 亚洲片人在线观看| 1024视频免费在线观看| 国产又黄又爽又无遮挡在线| 国产一区二区三区视频了| 18禁国产床啪视频网站| 真人一进一出gif抽搐免费| 色尼玛亚洲综合影院| 欧美日韩精品网址| 国产精品av视频在线免费观看| 久久亚洲真实| 天堂影院成人在线观看| 身体一侧抽搐| 淫秽高清视频在线观看| 老司机午夜福利在线观看视频| 免费在线观看日本一区| 一级黄色大片毛片| 亚洲精品久久成人aⅴ小说| 我要搜黄色片| 久久久久九九精品影院| 午夜老司机福利片| 特大巨黑吊av在线直播| 色综合欧美亚洲国产小说| 又爽又黄无遮挡网站| 欧美中文日本在线观看视频| 在线播放国产精品三级| 欧美日韩亚洲国产一区二区在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲成a人片在线一区二区| 欧美3d第一页| 亚洲国产中文字幕在线视频| 久久久久精品国产欧美久久久| 精品第一国产精品| 国产片内射在线| 欧美国产日韩亚洲一区| 国产aⅴ精品一区二区三区波| 最新美女视频免费是黄的| 欧美日韩国产亚洲二区| 午夜福利欧美成人| 亚洲国产欧美一区二区综合| 老司机福利观看| 亚洲av五月六月丁香网| 国产精品,欧美在线| 中亚洲国语对白在线视频| 亚洲男人的天堂狠狠| 岛国在线观看网站| 婷婷精品国产亚洲av在线| 91成年电影在线观看| av片东京热男人的天堂| 精品高清国产在线一区| 女人高潮潮喷娇喘18禁视频| 日韩av在线大香蕉| 欧美中文综合在线视频| 一区二区三区激情视频| 麻豆国产97在线/欧美 | 老汉色∧v一级毛片| 波多野结衣高清作品| 后天国语完整版免费观看| 国内精品久久久久精免费| 亚洲中文字幕一区二区三区有码在线看 | 国产成人啪精品午夜网站| 9191精品国产免费久久| 好男人电影高清在线观看| 免费看a级黄色片| 亚洲狠狠婷婷综合久久图片| aaaaa片日本免费| 亚洲av成人一区二区三| 亚洲成人中文字幕在线播放| 69av精品久久久久久| 美女扒开内裤让男人捅视频| 精品国产亚洲在线| 亚洲av电影在线进入| 免费在线观看黄色视频的| 中文字幕av在线有码专区| 99国产极品粉嫩在线观看| 国产精品av久久久久免费| 又爽又黄无遮挡网站| 婷婷亚洲欧美| 亚洲,欧美精品.| 最近在线观看免费完整版| 亚洲 欧美一区二区三区| 日韩精品青青久久久久久| 国产精品精品国产色婷婷| 欧美丝袜亚洲另类 | 日本成人三级电影网站| www.999成人在线观看| 久久精品成人免费网站| 99热这里只有是精品50| 91字幕亚洲| 亚洲性夜色夜夜综合| 成人三级做爰电影| 成人国产综合亚洲| 九色成人免费人妻av| 亚洲人与动物交配视频| 亚洲最大成人中文| 国产av一区二区精品久久| 高清毛片免费观看视频网站| 国产成人精品久久二区二区免费| 巨乳人妻的诱惑在线观看| 欧美黑人巨大hd| 老司机午夜福利在线观看视频| 黄色成人免费大全| 国产亚洲欧美在线一区二区| 欧美另类亚洲清纯唯美| 欧美日韩黄片免| 天堂√8在线中文| 精品一区二区三区四区五区乱码| 国产乱人伦免费视频| 深夜精品福利| 女人被狂操c到高潮| 手机成人av网站| 久久中文看片网| 久久香蕉激情| 亚洲午夜精品一区,二区,三区| 99久久精品热视频| 国产麻豆成人av免费视频| 欧美 亚洲 国产 日韩一| 免费av毛片视频| tocl精华| 日韩国内少妇激情av| 大型黄色视频在线免费观看| 长腿黑丝高跟| 亚洲性夜色夜夜综合| 女人爽到高潮嗷嗷叫在线视频| 伊人久久大香线蕉亚洲五| 欧美av亚洲av综合av国产av| 国产精品一区二区三区四区久久| 欧美在线一区亚洲| 国产成人啪精品午夜网站| 99国产精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 19禁男女啪啪无遮挡网站| 久久这里只有精品19| 日日夜夜操网爽| 日本黄色视频三级网站网址| 他把我摸到了高潮在线观看| 欧美成狂野欧美在线观看| 精品国内亚洲2022精品成人| 91字幕亚洲| 叶爱在线成人免费视频播放| 亚洲人成电影免费在线| 黄色丝袜av网址大全| 两性午夜刺激爽爽歪歪视频在线观看 | 九色国产91popny在线| 亚洲一区中文字幕在线| 日韩免费av在线播放| 草草在线视频免费看| 俄罗斯特黄特色一大片| 国产av又大| 伦理电影免费视频| 在线a可以看的网站| 男人舔女人的私密视频| 久久婷婷人人爽人人干人人爱| 丝袜美腿诱惑在线| 国产成人影院久久av| 蜜桃久久精品国产亚洲av| 中文资源天堂在线| 日本五十路高清| 最近视频中文字幕2019在线8| 久久性视频一级片| 精品一区二区三区av网在线观看| 99国产精品一区二区三区| 国产精品久久久久久精品电影| 国产成人aa在线观看| 老熟妇乱子伦视频在线观看| 精品久久久久久久毛片微露脸| 精品日产1卡2卡| 国产视频一区二区在线看| 国产精品日韩av在线免费观看| 精品午夜福利视频在线观看一区| 久久婷婷成人综合色麻豆| 我要搜黄色片| 亚洲av中文字字幕乱码综合| 久久精品国产清高在天天线| 国产精品亚洲av一区麻豆| 全区人妻精品视频| 亚洲aⅴ乱码一区二区在线播放 | 久久久久久久精品吃奶| 色综合婷婷激情| 久久久国产欧美日韩av| 后天国语完整版免费观看| 不卡av一区二区三区| 免费高清视频大片| 午夜久久久久精精品| 亚洲精品在线观看二区| 又大又爽又粗| 在线十欧美十亚洲十日本专区| 日韩成人在线观看一区二区三区| 香蕉丝袜av| 好男人电影高清在线观看| 亚洲成人精品中文字幕电影| 91成年电影在线观看| 欧美丝袜亚洲另类 | 国产精品一区二区精品视频观看| 国产三级在线视频| 亚洲国产欧美人成| 热99re8久久精品国产| 91老司机精品| 久久国产精品影院| 国产精品av视频在线免费观看| 欧美成人性av电影在线观看| 中出人妻视频一区二区| 国产伦一二天堂av在线观看| 午夜福利在线在线| 一个人免费在线观看的高清视频| 日韩大码丰满熟妇| 精品第一国产精品| 中文亚洲av片在线观看爽| 久久精品aⅴ一区二区三区四区| 亚洲一码二码三码区别大吗| 国产精品亚洲av一区麻豆| 国产伦在线观看视频一区| 男女午夜视频在线观看| 琪琪午夜伦伦电影理论片6080| 国产精品美女特级片免费视频播放器 | 亚洲欧美激情综合另类| 亚洲av中文字字幕乱码综合| 女同久久另类99精品国产91| 一本综合久久免费| 在线观看www视频免费| 国产精品久久久久久精品电影| 免费看美女性在线毛片视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲成av人片在线播放无| 欧美精品亚洲一区二区| 一边摸一边做爽爽视频免费| 亚洲中文av在线| 日本一二三区视频观看| 亚洲人成电影免费在线| 真人一进一出gif抽搐免费| 18禁黄网站禁片午夜丰满| 亚洲人与动物交配视频| 欧美色视频一区免费| 国内精品久久久久精免费| 黄色 视频免费看| 日韩欧美一区二区三区在线观看| 国产黄a三级三级三级人| a在线观看视频网站| www.熟女人妻精品国产| 久久婷婷成人综合色麻豆| 久久久精品欧美日韩精品| 精品日产1卡2卡| 黑人巨大精品欧美一区二区mp4| 日本精品一区二区三区蜜桃| 非洲黑人性xxxx精品又粗又长| 欧美绝顶高潮抽搐喷水| 欧美又色又爽又黄视频| 淫秽高清视频在线观看| 婷婷丁香在线五月| 熟女电影av网| 天堂动漫精品| 国产日本99.免费观看| 搡老熟女国产l中国老女人| 国产又色又爽无遮挡免费看| 一本久久中文字幕| 91在线观看av| 久久精品夜夜夜夜夜久久蜜豆 | 丁香六月欧美| 亚洲熟妇中文字幕五十中出| 天堂影院成人在线观看| 国产欧美日韩一区二区三| 精品午夜福利视频在线观看一区| 国产野战对白在线观看| 很黄的视频免费| 亚洲 欧美 日韩 在线 免费| 久久精品成人免费网站| 小说图片视频综合网站| 1024视频免费在线观看| 变态另类成人亚洲欧美熟女| 男女视频在线观看网站免费 | 一个人免费在线观看电影 | 日韩免费av在线播放| 一本久久中文字幕| 国产亚洲av高清不卡| 一个人观看的视频www高清免费观看 | 欧美久久黑人一区二区| 亚洲,欧美精品.| 国产精品 欧美亚洲| 国产真人三级小视频在线观看| 成年免费大片在线观看| 精品电影一区二区在线| 午夜精品久久久久久毛片777| 国产精品99久久99久久久不卡| 精品免费久久久久久久清纯| 亚洲专区字幕在线| 国产伦人伦偷精品视频| 50天的宝宝边吃奶边哭怎么回事| 中文字幕人妻丝袜一区二区| 精品久久久久久久久久免费视频| 亚洲18禁久久av| 国产精品永久免费网站| 欧美午夜高清在线| 国产精品 欧美亚洲| 在线观看一区二区三区| 国产精品久久电影中文字幕| 啦啦啦免费观看视频1| 亚洲一区高清亚洲精品| 婷婷精品国产亚洲av在线| 日韩有码中文字幕| 又粗又爽又猛毛片免费看| 香蕉国产在线看| 激情在线观看视频在线高清| 男女视频在线观看网站免费 | 在线观看免费日韩欧美大片| 亚洲中文字幕日韩| 青草久久国产| 亚洲电影在线观看av| 18禁裸乳无遮挡免费网站照片| 老熟妇乱子伦视频在线观看| 人成视频在线观看免费观看| 国产精品98久久久久久宅男小说| 麻豆av在线久日| 婷婷精品国产亚洲av| 免费人成视频x8x8入口观看| 亚洲专区国产一区二区| 高清毛片免费观看视频网站| 国产一级毛片七仙女欲春2| 黄片大片在线免费观看| 麻豆成人av在线观看| 欧美zozozo另类| 欧美 亚洲 国产 日韩一| 日韩欧美免费精品| 欧美精品啪啪一区二区三区| 一个人观看的视频www高清免费观看 | 国产精品野战在线观看| a级毛片a级免费在线| 麻豆一二三区av精品| 国产三级黄色录像| 真人做人爱边吃奶动态| 听说在线观看完整版免费高清| 亚洲国产欧美网| 露出奶头的视频| 亚洲av成人精品一区久久| 国产成年人精品一区二区| 正在播放国产对白刺激| 深夜精品福利| 丁香六月欧美| 国产精品自产拍在线观看55亚洲| 人人妻人人看人人澡| 亚洲18禁久久av| 人妻丰满熟妇av一区二区三区| 一二三四社区在线视频社区8| svipshipincom国产片| 久久久国产欧美日韩av| 少妇的丰满在线观看| 一级a爱片免费观看的视频| 欧美中文日本在线观看视频| 午夜福利欧美成人| 亚洲男人的天堂狠狠| 日本 av在线| 欧美日韩亚洲国产一区二区在线观看| 在线永久观看黄色视频| svipshipincom国产片| 丰满的人妻完整版| 亚洲在线自拍视频| 老汉色av国产亚洲站长工具| 日本一本二区三区精品| 长腿黑丝高跟| 亚洲午夜理论影院| 极品教师在线免费播放| 亚洲国产精品合色在线| 在线观看午夜福利视频| 妹子高潮喷水视频| 欧美日韩福利视频一区二区| 免费看a级黄色片| 久久精品国产99精品国产亚洲性色| 日韩欧美精品v在线| 性欧美人与动物交配| 五月伊人婷婷丁香| 法律面前人人平等表现在哪些方面| 色噜噜av男人的天堂激情| 亚洲av电影不卡..在线观看| 精品国产美女av久久久久小说| 久久久久精品国产欧美久久久| 亚洲,欧美精品.| 人妻丰满熟妇av一区二区三区| 亚洲精品粉嫩美女一区| 色综合婷婷激情| 国内毛片毛片毛片毛片毛片| 国产精品乱码一区二三区的特点| 老鸭窝网址在线观看| 婷婷精品国产亚洲av| 午夜免费激情av| 我的老师免费观看完整版| 久9热在线精品视频| 日韩欧美 国产精品| 99热6这里只有精品| 亚洲成人精品中文字幕电影| √禁漫天堂资源中文www| 丰满人妻一区二区三区视频av | 好男人在线观看高清免费视频| 老汉色∧v一级毛片| 日本一区二区免费在线视频| 亚洲欧美日韩无卡精品| 亚洲第一电影网av| 男女床上黄色一级片免费看| 午夜日韩欧美国产| a在线观看视频网站| 亚洲 欧美 日韩 在线 免费| 99久久国产精品久久久| 97人妻精品一区二区三区麻豆| 国产精品久久电影中文字幕| 日本精品一区二区三区蜜桃| 亚洲av成人一区二区三| 亚洲国产精品成人综合色| 国产精品免费视频内射| 两性夫妻黄色片| 好男人电影高清在线观看| 欧美黄色淫秽网站| 黄色女人牲交| 窝窝影院91人妻| 叶爱在线成人免费视频播放| 国产精品98久久久久久宅男小说| 最近在线观看免费完整版| 国产精品影院久久| 精品久久久久久成人av| 中文字幕人妻丝袜一区二区| 韩国av一区二区三区四区| 69av精品久久久久久| 亚洲熟女毛片儿| tocl精华| 久久人人精品亚洲av| 日韩欧美在线二视频| 1024香蕉在线观看| 麻豆一二三区av精品| 精品免费久久久久久久清纯| 亚洲熟女毛片儿| 国产麻豆成人av免费视频| 成人三级黄色视频| 国产aⅴ精品一区二区三区波| 女生性感内裤真人,穿戴方法视频| 亚洲一区二区三区不卡视频| 中国美女看黄片| 一进一出抽搐gif免费好疼| 国产69精品久久久久777片 | 精品久久久久久久末码| 亚洲中文字幕一区二区三区有码在线看 | 桃色一区二区三区在线观看| 午夜a级毛片| 淫妇啪啪啪对白视频| av有码第一页| 国产精品国产高清国产av| 国产亚洲av嫩草精品影院| 亚洲av成人不卡在线观看播放网| 精品日产1卡2卡| 免费在线观看完整版高清| 国产1区2区3区精品| а√天堂www在线а√下载| 久久国产精品影院| 国产成人aa在线观看| 美女免费视频网站| 亚洲av成人精品一区久久| 最近视频中文字幕2019在线8| 精品福利观看| 国产熟女午夜一区二区三区| 亚洲国产精品sss在线观看| 中文字幕最新亚洲高清| 久久久久久国产a免费观看| 国产免费男女视频| 午夜福利视频1000在线观看| 欧美又色又爽又黄视频| 国内久久婷婷六月综合欲色啪| 亚洲熟妇熟女久久| 国产激情欧美一区二区| av视频在线观看入口| 国产午夜福利久久久久久| av片东京热男人的天堂| 午夜免费观看网址| 国产亚洲精品av在线| 99久久无色码亚洲精品果冻| 亚洲精品中文字幕在线视频| 老司机午夜福利在线观看视频| 久久中文看片网| 亚洲成av人片在线播放无| 色精品久久人妻99蜜桃| 嫩草影视91久久| 欧美中文综合在线视频| 久久久久久国产a免费观看| 日韩精品中文字幕看吧| 午夜精品在线福利| 午夜影院日韩av| 精品欧美一区二区三区在线| av欧美777| 伦理电影免费视频| 五月玫瑰六月丁香| 一卡2卡三卡四卡精品乱码亚洲| 欧美另类亚洲清纯唯美| 国产一区在线观看成人免费| 女人被狂操c到高潮| 亚洲成人精品中文字幕电影| 国产亚洲av嫩草精品影院| 国产精品一区二区三区四区免费观看 | 亚洲av成人av| 国产av麻豆久久久久久久| 久久精品综合一区二区三区| 日韩欧美在线乱码| 黄色视频,在线免费观看| 天天一区二区日本电影三级| 免费搜索国产男女视频| 一级黄色大片毛片| 禁无遮挡网站| 成人av一区二区三区在线看| bbb黄色大片| 日韩免费av在线播放| 老司机午夜十八禁免费视频| 两个人的视频大全免费| 中文字幕精品亚洲无线码一区| 男人舔女人的私密视频| 亚洲成a人片在线一区二区| 97人妻精品一区二区三区麻豆| 欧美日韩亚洲国产一区二区在线观看| 国产精品av视频在线免费观看| 最新美女视频免费是黄的| 亚洲精品美女久久av网站| 国产成人av激情在线播放| 亚洲国产欧美人成| 亚洲一区二区三区不卡视频| 国产成人精品久久二区二区91| 又紧又爽又黄一区二区| 制服诱惑二区| 国产三级中文精品| 免费看日本二区| 亚洲 欧美 日韩 在线 免费| 一本一本综合久久| 久久精品人妻少妇| 91麻豆av在线| 人人妻,人人澡人人爽秒播| 午夜老司机福利片| 日韩精品免费视频一区二区三区| 一区二区三区国产精品乱码| 黄频高清免费视频| 激情在线观看视频在线高清| 国产成人av教育| 丰满的人妻完整版| 一本综合久久免费| 1024手机看黄色片| 日本 av在线| 国产亚洲精品综合一区在线观看 | 久久精品人妻少妇| 午夜久久久久精精品| 在线观看美女被高潮喷水网站 | 亚洲国产精品成人综合色| 伊人久久大香线蕉亚洲五| 国产人伦9x9x在线观看| 日本免费a在线| 麻豆久久精品国产亚洲av| 欧美成人免费av一区二区三区| 日本精品一区二区三区蜜桃| 婷婷亚洲欧美| 午夜激情福利司机影院| 午夜福利欧美成人| 在线国产一区二区在线| 午夜福利在线观看吧| 亚洲自拍偷在线| 免费观看精品视频网站| 99久久无色码亚洲精品果冻| 日本五十路高清| 国内精品久久久久精免费| 男女视频在线观看网站免费 | 亚洲第一电影网av| 精品国产乱码久久久久久男人| 成人av一区二区三区在线看| 高清在线国产一区| 日本免费a在线| 精品电影一区二区在线| 国产aⅴ精品一区二区三区波| 亚洲欧美一区二区三区黑人| 午夜视频精品福利| 级片在线观看| 亚洲精品中文字幕一二三四区| www.自偷自拍.com| 久久久久久久久免费视频了| 亚洲精品中文字幕一二三四区| 国产三级中文精品|