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

    筋的截面形狀對(duì)薄壁結(jié)構(gòu)振動(dòng)疲勞性能的影響

    2017-06-13 10:44:14劉雙燕涂玉倩苗應(yīng)剛李玉龍
    航空工程進(jìn)展 2017年2期
    關(guān)鍵詞:筋條固支簡(jiǎn)支

    劉雙燕,涂玉倩,苗應(yīng)剛,李玉龍

    (西北工業(yè)大學(xué) 航空學(xué)院,西安 710072)

    筋的截面形狀對(duì)薄壁結(jié)構(gòu)振動(dòng)疲勞性能的影響

    劉雙燕,涂玉倩,苗應(yīng)剛,李玉龍

    (西北工業(yè)大學(xué) 航空學(xué)院,西安 710072)

    筋條的截面形狀對(duì)薄壁結(jié)構(gòu)的振動(dòng)疲勞性能具有重要影響。對(duì)L型和T型加筋薄板進(jìn)行振動(dòng)試驗(yàn),研究固支條件下L型筋和T型筋對(duì)薄板頻率和振動(dòng)疲勞壽命的影響,并從理論角度分析其原因;利用MSC有限元軟件研究在正弦分布載荷下六種不同截面形狀(矩形、L型、T型、I型、Z型、∏型)的筋條在給定邊界條件(固支、簡(jiǎn)支)下,對(duì)薄板的頻率、應(yīng)力和振動(dòng)疲勞壽命的影響,并對(duì)加筋板的壽命進(jìn)行預(yù)測(cè)。結(jié)果表明:六種加筋構(gòu)型中,Z型加筋板的第一階固有頻率最大,加筋板Von Mises應(yīng)力最大點(diǎn)出現(xiàn)在同樣的位置;固支條件下,加筋板的破壞位置在板的左右兩側(cè),Z型加筋板壽命最長(zhǎng);簡(jiǎn)支條件下,破壞位置出現(xiàn)在板的中央,L型加筋板的壽命最長(zhǎng)。

    薄壁加筋結(jié)構(gòu);截面形狀;固有頻率;Von Mises;振動(dòng)疲勞壽命

    0 引 言

    薄壁結(jié)構(gòu)廣泛應(yīng)用于航空、航天、汽車、船舶等領(lǐng)域。由于該類結(jié)構(gòu)剛度、強(qiáng)度、屈曲、振動(dòng)等力學(xué)性能薄弱,在實(shí)際應(yīng)用中一般都采用加筋構(gòu)型,例如飛機(jī)機(jī)翼整體壁板、發(fā)動(dòng)機(jī)機(jī)匣、導(dǎo)彈殼體、航天飛機(jī)的外部燃料箱等都屬于典型的加筋結(jié)構(gòu)。筋條的形狀、尺寸與布局直接影響著結(jié)構(gòu)的重量與性能。

    國(guó)內(nèi)外很多學(xué)者研究了筋條的布局、尺寸對(duì)結(jié)構(gòu)強(qiáng)度的影響。Y.C.Lam等[1]采用變厚度法系統(tǒng)地研究了板結(jié)構(gòu)的加筋布局問題。J.Luo等[2]提出了一種考慮筋條位置與方向的加筋設(shè)計(jì)方法,并用于靜載和特征值問題。L.A.Krog等[3]研究了靜載和自由振動(dòng)下含加強(qiáng)肋板結(jié)構(gòu)的拓?fù)鋬?yōu)化問題。D.Bojczuk等[4]采用加筋板單元研究了平面內(nèi)載荷和彎曲載荷下Kirchhof板的加筋布局和形狀優(yōu)化。A.Kallassy等[5]和J.Fatemi等[6]以筋條位置為變量、剛度為目標(biāo)研究了平板結(jié)構(gòu)的加筋布局。

    S.H.Crandall等[7]首次提出振動(dòng)疲勞的概念,并將其定義為結(jié)構(gòu)在振動(dòng)疲勞下產(chǎn)生的不可逆累積性損傷或破壞。H.Y.Liou等[8]基于隨機(jī)振動(dòng)理論提出了一個(gè)用于估計(jì)結(jié)構(gòu)疲勞壽命的修正模型。王明珠等[9]提出了一種結(jié)構(gòu)隨機(jī)振動(dòng)疲勞壽命估算的樣本法,通過(guò)該樣本法能夠處理在頻域內(nèi)用譜密度描述的寬帶隨機(jī)振動(dòng)載荷的情況。張積亭等[10]在飛機(jī)典型構(gòu)件振動(dòng)疲勞壽命估計(jì)中提出,利用隨機(jī)響應(yīng)功率譜密度求出的特征頻率作為平均頻率,在數(shù)據(jù)處理中是一種較方便的方法。郭建平等[11]針對(duì)機(jī)載電子設(shè)備隨機(jī)振動(dòng)問題,采用MSC.FATIGUE軟件進(jìn)行隨機(jī)振動(dòng)疲勞分析,計(jì)算了疲勞壽命大小及分布;通過(guò)對(duì)比計(jì)算結(jié)果與試驗(yàn)結(jié)果,證明了利用隨機(jī)振動(dòng)疲勞分析方法預(yù)測(cè)此類產(chǎn)品的疲勞壽命是可行的。陳安等[12]采用了一系列試驗(yàn)方法研究鋁鋰合金鉚接結(jié)構(gòu)疲勞性能。

    加筋板和振動(dòng)疲勞已有很多研究成果,但是將兩者結(jié)合的研究仍鮮有報(bào)道。本文主要內(nèi)容是將二者結(jié)合起來(lái),研究在固支和簡(jiǎn)支條件下筋條的截面形狀對(duì)加筋板振動(dòng)疲勞性能的影響,為航空航天等領(lǐng)域的應(yīng)用提供理論和試驗(yàn)依據(jù)。

    1 振動(dòng)疲勞試驗(yàn)

    1.1 試驗(yàn)件設(shè)計(jì)

    為了模擬航空結(jié)構(gòu)中典型的加筋薄壁結(jié)構(gòu)特征,試驗(yàn)件設(shè)計(jì)為雙加筋板,平板厚度為1.6 mm,加筋板平面尺寸及筋條布局如圖1所示。筋條的截面形狀分別選用L型和T型,如圖2所示(筋條長(zhǎng)度均為200 mm)。

    薄壁結(jié)構(gòu)的平板和筋條分別為2024-T3和LY12CZ鋁合金,兩種材料的化學(xué)成分如表1所示,性能參數(shù)如下:彈性模量E分別為72和70 GPa,泊松比ν均為0.3,密度ρ為2 770和2 800 kg/m3,屈服極限σs為275和270 MPa。為了保證試驗(yàn)數(shù)據(jù)的可靠性,L型和T型試樣各取6件,并采用6件不加筋的平板作為參照試樣組。

    表1 鋁合金成分表

    試驗(yàn)中設(shè)計(jì)的夾具由45號(hào)鋼制成,分為a、b兩部分,如圖3所示。a部分通過(guò)8枚M12螺栓安裝固定在振動(dòng)臺(tái)臺(tái)面上,將薄壁結(jié)構(gòu)試樣水平放置在夾具a部分的上表面,之后將夾具b部分置于試樣上方,用20個(gè)M8沉頭六角螺栓將夾具a部分、試樣和夾具b部分連接,模擬四邊固支的邊界條件。

    1.2 試驗(yàn)方案

    試驗(yàn)所使用的設(shè)備包括m+p振動(dòng)控制儀及配套的Vibcontrol控制軟件、東菱電動(dòng)振動(dòng)臺(tái)、BE120-3AA應(yīng)變片(規(guī)格為120 Ω)、東華DH3841應(yīng)變放大器和DEWE-801動(dòng)態(tài)信號(hào)采集系統(tǒng)。

    薄壁結(jié)構(gòu)的振動(dòng)疲勞試驗(yàn)原理如圖4所示。貼有應(yīng)變片的試樣被夾具以四邊固支的方式安裝夾持在振動(dòng)臺(tái)上如圖5所示,應(yīng)變片位于圖中數(shù)字標(biāo)號(hào)附近。

    振動(dòng)疲勞激勵(lì)系統(tǒng)由m+p Vibcontrol振動(dòng)控制儀和東菱電動(dòng)振動(dòng)臺(tái)組成。振動(dòng)激勵(lì)的加載譜等振動(dòng)參數(shù)由振動(dòng)控制儀傳給振動(dòng)臺(tái),從而驅(qū)動(dòng)振動(dòng)臺(tái)以產(chǎn)生相應(yīng)的振動(dòng),加速度計(jì)將采集到的振動(dòng)信號(hào)反饋給控制儀,控制儀根據(jù)控制譜的要求對(duì)激勵(lì)信號(hào)進(jìn)行實(shí)時(shí)修正以保證振動(dòng)疲勞加載的持續(xù)性。試驗(yàn)過(guò)程中試樣的固有頻率和應(yīng)變均方根(Root Mean Square,簡(jiǎn)稱RMS)值變化歷程由應(yīng)變片采集到的動(dòng)態(tài)應(yīng)變信號(hào)經(jīng)動(dòng)態(tài)應(yīng)變放大器放大后輸入DEWE-801動(dòng)態(tài)信號(hào)采集系統(tǒng)進(jìn)行存儲(chǔ)和后續(xù)處理分析。

    試驗(yàn)加載歷程如下:首先對(duì)各試樣施加1 g正弦掃頻載荷,得到相應(yīng)的頻率響應(yīng)曲線和第一階固有頻率;然后,對(duì)各試樣加載頻帶范圍在其固有頻率0.8~1.2倍之間的幅值為10 g2/Hz的隨機(jī)振動(dòng)載荷譜,以固有頻率降低5%[13]作為試樣疲勞破壞的判據(jù),對(duì)各試樣進(jìn)行足夠長(zhǎng)時(shí)間的振動(dòng),從而得到其振動(dòng)疲勞壽命。對(duì)薄壁結(jié)構(gòu)采用隨機(jī)激勵(lì),而不是正弦或正弦駐留等激勵(lì)方式,是因?yàn)楸”诮Y(jié)構(gòu)比梁結(jié)構(gòu)復(fù)雜,用正弦激勵(lì)難以保證加載的穩(wěn)定性,而且隨機(jī)激勵(lì)方式更接近薄壁結(jié)構(gòu)服役時(shí)的真實(shí)載荷環(huán)境。

    對(duì)三種板各進(jìn)行六次重復(fù)試驗(yàn),將每種板得到的壽命值進(jìn)行統(tǒng)計(jì)分析,得到其壽命在置信度為95%時(shí)的上下置信區(qū)間和中值疲勞壽命,結(jié)果如表2所示。

    表2 試樣的振動(dòng)疲勞試驗(yàn)結(jié)果

    2 振動(dòng)疲勞試驗(yàn)的模擬驗(yàn)證

    2.1 有限元建模

    利用有限元分析軟件MSC.PATRAN 建立試驗(yàn)件的有限元模型。平板和筋條的尺寸與試驗(yàn)保持一致,平板的四邊采用固支的邊界條件,使得其有效平面尺寸為220 mm×220 mm,板結(jié)構(gòu)為殼單元,網(wǎng)格細(xì)化到5 mm,筋條采用實(shí)體單元,筋條與平板之間為綁定連接,T型加筋板的有限元模型如圖6所示。

    在有限元模擬中,為了得到準(zhǔn)確的振動(dòng)疲勞壽命,需要輸入結(jié)構(gòu)的品質(zhì)數(shù)Q[14],

    (1)

    式中:f0為固有頻率;Δf為半功率帶寬,Δf=f2-f1。f0、f1和f2可以從試驗(yàn)所得對(duì)應(yīng)加筋板的頻率響應(yīng)曲線得到。

    利用MSC.PATRAN及MSC.NASTRAN,結(jié)合品質(zhì)數(shù)模擬得到各板的頻率響應(yīng)關(guān)系,再利用MSC.FATIGUE有限元軟件對(duì)各板施加試驗(yàn)過(guò)程中所加載的載荷譜,并進(jìn)行振動(dòng)疲勞壽命分析,得到平板、T型和L型加筋板的損傷圖,從而得到其振動(dòng)疲勞壽命。

    2.2 模擬驗(yàn)證結(jié)果

    通過(guò)試驗(yàn)和模擬得到平板、T型和L型加筋板的固有頻率和振動(dòng)疲勞壽命,結(jié)果如表3所示。

    表3 薄壁結(jié)構(gòu)的固有頻率和振動(dòng)疲勞壽命比較

    從表3可以看出:

    (1) 試驗(yàn)測(cè)得的頻率比模擬的頻率略高,主要原因是因?yàn)樵囼?yàn)過(guò)程中板的第一階固有頻率是通過(guò)施加1 g的正弦掃頻載荷得到頻率響應(yīng)曲線,并且試驗(yàn)中的結(jié)構(gòu)較薄,在振動(dòng)過(guò)程中板的中央位置的振幅大于板的厚度,從而使板的幾何剛度不容忽視,而有限元模擬采用殼單元,忽略了該項(xiàng),試驗(yàn)中剛度項(xiàng)需附加上幾何剛度,所以試驗(yàn)結(jié)果偏高。

    (2) 不管是試驗(yàn)結(jié)果還是數(shù)值結(jié)果,加筋板的第一階固有頻率均高于不加筋板,這主要是由于筋的引入增加了結(jié)構(gòu)質(zhì)量,同時(shí)也提高了結(jié)構(gòu)的剛度,并且剛度的增加程度遠(yuǎn)大于結(jié)構(gòu)質(zhì)量的增加。因此工程中常使用加筋板代替不加筋板,以獲得高固有頻率,延遲達(dá)到共振的時(shí)間。

    (3) 對(duì)于振動(dòng)疲勞壽命,試驗(yàn)得到的結(jié)果比數(shù)值模擬的結(jié)果偏低,主要原因是板在加工過(guò)程中存在內(nèi)部缺陷或表面微裂紋,導(dǎo)致試驗(yàn)結(jié)果比模擬結(jié)果提前破壞;試驗(yàn)中三種板的疲勞壽命相差不大,考慮到三種板的加載頻率均是在其本身共振帶內(nèi),而兩種加筋板的第一階固有頻率遠(yuǎn)高于不加筋板,所以從振動(dòng)疲勞循環(huán)次數(shù)上看,兩種加筋板的振動(dòng)疲勞壽命明顯高于不加筋板。

    綜上所述,利用有限元軟件模擬得到試樣的固有頻率和試驗(yàn)結(jié)果誤差在5%以內(nèi),驗(yàn)證了數(shù)值模擬的可靠性。由于振動(dòng)疲勞試樣結(jié)果分散性較大,業(yè)內(nèi)認(rèn)為模擬結(jié)果和試驗(yàn)結(jié)果在一個(gè)數(shù)量級(jí)即認(rèn)為具有很好的吻合性,本文振動(dòng)疲勞壽命試驗(yàn)和模擬誤差在50%以內(nèi),說(shuō)明本文采用的邊界條件、數(shù)值模型和材料參數(shù)是合理的。

    對(duì)于試驗(yàn)和模擬中的加筋板而言,T型加筋板的固有頻率和振動(dòng)疲勞壽命均高于L型加筋板。因此,有必要深入研究筋的不同截面形狀對(duì)薄板的振動(dòng)疲勞性能的影響,本文從有限元模擬角度來(lái)探討筋的截面形狀對(duì)薄壁結(jié)構(gòu)的頻率、應(yīng)力和壽命的影響。

    3 筋的不同截面形狀對(duì)薄板振動(dòng)疲勞性能的影響

    3.1 有限元模擬

    為了研究筋條的截面形狀對(duì)薄板的振動(dòng)疲勞性能的影響,對(duì)焊接了不同截面形狀筋條的加筋板施加單頻正弦壓強(qiáng)載荷,載荷幅值為8 500 Pa,加載頻率為加筋板的第一階固有頻率,邊界條件分為四邊固支與四邊簡(jiǎn)支兩種。加筋結(jié)構(gòu)的平板尺寸為220 mm×220 mm×1.6 mm,板和筋的材料性能與1.1小節(jié)一致。筋的截面形狀分為六種:矩形、L型、T型、I型、Z型、∏型。焊接于板上的兩條筋間距為150 mm,且在板上對(duì)稱分布。為保證所焊筋條的重量相等,六種筋的橫截面面積相等且筋與板的接觸面積相等,具體形狀和尺寸如圖7所示。首先對(duì)含有六種截面形狀的加筋板進(jìn)行模態(tài)分析,得到其模態(tài)頻率;然后對(duì)各試樣加載頻帶范圍在其第一階固有頻率的0.8~1.2倍之間的幅值為10 g2/Hz的隨機(jī)振動(dòng)載荷譜,計(jì)算其應(yīng)力及壽命值。

    3.2 結(jié)果與討論

    3.2.1 對(duì)固有頻率的影響

    對(duì)于具有n自由度的線性阻尼振動(dòng)系統(tǒng),其運(yùn)動(dòng)微分方程為[13]

    (2)

    若忽略系統(tǒng)阻尼的影響,自由振動(dòng)的運(yùn)動(dòng)微分方程表達(dá)為

    (3)

    設(shè)該方程的解為x=Φejwt,其中Φ為n階振幅列陣,將其帶入式(3)得到

    K-ω2M=0

    (4)

    當(dāng)系統(tǒng)發(fā)生振動(dòng)時(shí),式(4)為齊次線性方程組,具有非零解的條件為其系數(shù)矩陣的行列式等于零,即

    |K-ω2M|=0

    (5)

    式(5)稱為系統(tǒng)的特征方程,將其展開可獲得以ω2為未知量的n次代數(shù)方程,求解得到n個(gè)ω2值,即系統(tǒng)的n個(gè)固有頻率,由小到大排列,最小的稱之為第一階固有頻率(基本頻率),其余依次為第二階固有頻率、第三階固有頻率等。通過(guò)式(5)可以求出任意結(jié)構(gòu)的固有頻率,對(duì)于任意給定結(jié)構(gòu),其固有頻率的平方值與其質(zhì)量成反比,與其剛度成反比。

    對(duì)于質(zhì)量相同、筋條截面形狀不同的加筋板來(lái)說(shuō),在給定的邊界條件下,結(jié)構(gòu)剛度不同,其固有頻率必然不同。四邊固支和簡(jiǎn)支條件下,六種不同截面形狀的加筋板結(jié)構(gòu)的前三階固有頻率值如圖8所示。

    從圖8可以看出:在固支條件下加筋板的第一階固有頻率值普遍高于在簡(jiǎn)支條件下的值,其根本原因在于固支條件下相應(yīng)結(jié)構(gòu)的邊界約束較大,使得其結(jié)構(gòu)剛度顯著提高;在任何同一邊界條件下,其第一階固有頻率值的關(guān)系為Z型>L型>T型>I型>∏型>矩形,這是因?yàn)樵谝浑A模態(tài)下,結(jié)構(gòu)承受抗彎剛度,而抗彎剛度取決于其在撓度方向的慣性矩[15],經(jīng)過(guò)計(jì)算六種截面形狀的筋的慣性矩,得到其大小關(guān)系為Z型>L型>T型>I型>∏型>矩形,這與加筋板第一階固有頻率的大小關(guān)系一致。

    3.2.2 對(duì)結(jié)構(gòu)VonMises應(yīng)力分布的影響

    隨機(jī)振動(dòng)條件下,筋條在薄板的布局相同,因此在同樣的邊界條件,其VonMises應(yīng)力分布云圖基本不變。固支邊界條件作用下,I型加筋板的VonMises應(yīng)力分布云圖如圖9(a)所示,位置1、位置2和位置3分別是第一階固有頻率作用下出現(xiàn)局部應(yīng)力最大值的位置。簡(jiǎn)支條件下,在第一階固有頻率時(shí),加筋板的VonMises應(yīng)力分布云圖如圖9(b)所示,局部應(yīng)力最大值的位置分布出現(xiàn)在位置1、位置2和位置4。

    從圖1和圖9可以看出:位置1在加筋板的中心位置,由于所加載荷頻帶處于第一階固有頻率的0.8~1.2倍的共振區(qū)域,中心區(qū)域的應(yīng)變最大,而所研究的結(jié)構(gòu)在線彈性范圍內(nèi),因而在加筋板的中心區(qū)域出現(xiàn)了局部應(yīng)力最大值;位置2在筋條連接在薄板位置的區(qū)域,局部質(zhì)量增加易出現(xiàn)應(yīng)力集中;固支條件下,位置3出現(xiàn)局部應(yīng)力集中,而在簡(jiǎn)支條件下,局部應(yīng)力集中現(xiàn)象沒有出現(xiàn)在相應(yīng)位置,而是出現(xiàn)在位置4,原因在于簡(jiǎn)支條件下,加筋板結(jié)構(gòu)四周邊界處可以轉(zhuǎn)動(dòng),而筋條增加了局部剛度,阻礙了板的轉(zhuǎn)動(dòng),因而在筋條與邊界最近處位置4的位置出現(xiàn)應(yīng)力集中;在簡(jiǎn)支條件下應(yīng)力集中的區(qū)域位置1和位置2的范圍都比在固支條件下的大,究其原因在于簡(jiǎn)支條件加筋板的四周邊界處可以轉(zhuǎn)動(dòng),使得其在振動(dòng)條件下相應(yīng)區(qū)域內(nèi)發(fā)生形變更容易,沿面外方向的應(yīng)變方向更大;特別在中心區(qū)域的應(yīng)力值也更大,在位置1處,簡(jiǎn)支條件下加筋板的最大VonMises應(yīng)力為174MPa,而在固支條件下其最大值僅為142MPa。

    事實(shí)上,同樣的加筋構(gòu)型下,簡(jiǎn)支條件下的局部VonMises應(yīng)力峰值一般大于固支條件下的峰值,四邊固支/四邊簡(jiǎn)支條件下加筋板的局部VonMises應(yīng)力峰值如表4所示。

    表4 四邊固支和四邊簡(jiǎn)支條件下加筋板的局部Von Mises應(yīng)力

    從表4可以看出:①在固支條件下,比較點(diǎn)焊不同截面形狀筋條的加筋板局部應(yīng)力峰值大小發(fā)現(xiàn),位置1處的VonMises峰值最大,位置3處的VonMises峰值最小;②在簡(jiǎn)支條件下,位置4處的VonMises峰值最大,位置1處的VonMises峰值最小,這兩種現(xiàn)象的出現(xiàn),是由于邊界約束條件和筋的耦合作用產(chǎn)生的;③在固支條件下,截面形狀不同,其在中心區(qū)域位置1的局部VonMises值不同,其大小排序?yàn)榫匦?∏型>I型>T型>Z型>L型,因?yàn)榻Y(jié)構(gòu)處于一階共振頻帶內(nèi),其振動(dòng)形式為沿著薄板法線方向上下振動(dòng),并且研究的結(jié)構(gòu)是線彈性結(jié)構(gòu),其應(yīng)力值與應(yīng)變/位移成線性關(guān)系,故結(jié)構(gòu)的抗彎剛度越大,其變形越小,其局部VonMises值越小,對(duì)于具有同樣質(zhì)量的加筋板來(lái)說(shuō),其固有頻率越大,因而其中心區(qū)域的VonMises峰值大小排序正好與其第一階固有頻率大小關(guān)系相反;④含矩形截面筋條的加筋板結(jié)構(gòu)在簡(jiǎn)支條件已達(dá)到屈服,這是因?yàn)槠湓诤?jiǎn)支條件下剛度太小。

    3.2.3 對(duì)疲勞壽命的影響

    大多數(shù)工程結(jié)構(gòu)的振動(dòng)疲勞失效往往是在一系列變幅循環(huán)載荷作用下產(chǎn)生疲勞損傷并累積而造成的。線性疲勞累積損傷理論是指材料或結(jié)構(gòu)在循環(huán)載荷作用下,其疲勞損傷可線性累加,各個(gè)量級(jí)應(yīng)力之間相互獨(dú)立并且互不相關(guān)。當(dāng)損傷累積到某一數(shù)值時(shí),結(jié)構(gòu)件就發(fā)生疲勞破壞。線性累積損傷理論中最經(jīng)典和常用的是Palmgren-Miner理論,簡(jiǎn)稱Miner理論[14]。

    其表達(dá)式為

    D=∑(ni/Ni)

    (6)

    式中:D為疲勞損傷;ni為某應(yīng)力范圍和平均應(yīng)力下的循環(huán)次數(shù);i為范圍變量;Ni為該應(yīng)力范圍和平均應(yīng)力下試樣完全疲勞破壞所對(duì)應(yīng)循環(huán)次數(shù)。

    按照Miner線性疲勞損傷累積準(zhǔn)則,當(dāng)D=1時(shí),試樣完全疲勞失效破壞。根據(jù)Miner理論計(jì)算得到在不同邊界和加筋構(gòu)型條件下振動(dòng)疲勞破壞區(qū)域和壽命值。四邊簡(jiǎn)支的加筋板結(jié)構(gòu)振動(dòng)破壞位置在板的中央位置位置1,而四邊簡(jiǎn)支的加筋板結(jié)構(gòu)的振動(dòng)疲勞破壞位置在板的左右兩邊位置,即位置3,如圖9所示。具體不同構(gòu)型加筋板的振動(dòng)疲勞壽命如圖10所示。

    從圖10可以看出:除了矩形筋外的五種筋條構(gòu)型對(duì)板振動(dòng)疲勞壽命的影響不隨邊界約束條件的不同而改變,其關(guān)系為L(zhǎng)型>Z型>T型>I型>∏型,這與加筋板在破壞位置的應(yīng)力大小關(guān)系相反,即破壞位置應(yīng)力越大,壽命越短。

    4 結(jié) 論

    (1) 通過(guò)平板及加筋板的振動(dòng)疲勞試驗(yàn)和有限元數(shù)值模擬結(jié)果的對(duì)比,驗(yàn)證了所采取模擬過(guò)程方法的準(zhǔn)確性和可行性。

    (2) 在筋條質(zhì)量相等且與平板的接觸面積相同的情況下,采用同樣的邊界條件,當(dāng)筋條的截面撓度方向的慣性矩較大時(shí),筋的抗彎剛度較大,則加筋結(jié)構(gòu)的第一階固有頻率較大;相應(yīng)地,若不同截面形狀的加筋板結(jié)構(gòu)在同樣的外載荷作用下,筋條的截面慣性矩較大時(shí),對(duì)應(yīng)的加筋結(jié)構(gòu)危險(xiǎn)點(diǎn)處的應(yīng)力會(huì)較小,則此種加筋結(jié)構(gòu)的振動(dòng)疲勞壽命會(huì)較長(zhǎng)。

    [1]LamYC,SanthikumarS.Automatedriblocationandoptimizationforplatestructures[J].Structural&MultidisciplinaryOptimization, 2003, 25(1): 35-45.

    [2]LuoJ,GeaHC.Asystematictopologyoptimizationapproachforoptimalstiffenerdesign[J].StructuralOptimization, 1998, 16(4): 280-288.

    [3]KrogLA,OlhoffN.Optimumtopologyandreinforcementdesignofdiskandplatestructureswithmultiplestiffnessandeigenfrequencyobjectives[J].Computers&Structures, 1999, 72(4-5): 535-563.

    [4]BojczukD,SzteleblakW.Optimizationoflayoutandshapeofstiffenersin2Dstructures[J].Computers&Structures, 2008, 86(13-14): 1436-1446.

    [5]KallassyA,MarcelinJL.Optimizationofstiffenedplatesbygeneticsearch[J].StructuralOptimization, 1997, 13(2): 134-141.

    [6]FatemiJ,TrompetteP.Optimaldesignofstiffenedplatestructures[R].AIAA-2002-1672, 2002.

    [7]CrandallSH,MarkWD.Randomvibrationinmechanicalsystems[M].NewYork:AcedemicPress,1963: 90-92.

    [8]LiouHY,WuWF,ShinCS.Amodifiedmodelfortheestimationoffatiguelifederivedfromrandomvibrationtheory[J].ProbabilisticEngineeringMechanics, 1999, 14(3): 281-288.

    [9] 王明珠, 姚衛(wèi)星, 孫偉. 結(jié)構(gòu)隨機(jī)振動(dòng)疲勞壽命估算的樣本法[J]. 中國(guó)機(jī)械工程, 2008, 19(8): 972-975.WangMingzhu,YaoWeixing,SunWei.Sampleapproachforfatiguelifepredictionofstructuresunderrandomvibration[J].ChinaMechanicalEngineering, 2008, 19(8): 972-975.(inChinese)

    [10] 張積亭, 周蘇楓. 飛機(jī)典型構(gòu)件振動(dòng)疲勞壽命分析[J]. 機(jī)械科學(xué)與技術(shù), 2002, 21(增刊): 3,16.ZhangJiting,ZhouSufeng.Analysisofvibrationfratiguelifeoftypicalaircraftcomponents[J].MechanicalScienceandTechnology, 2002, 21(S1): 3,16.(inChinese)

    [11] 郭建平, 任康, 楊龍, 等. 基于MSC.Fatigue的電子設(shè)備隨機(jī)振動(dòng)疲勞分析[J]. 航空計(jì)算技術(shù), 2008, 3(4): 48-50.GuoJianping,RenKang,YangLong,etal.FatigueanalysisonrandomvibrationofelectronicequipmentbasedonMSC.Fatiguesoftware[J].AeronauticalComputingTechnique, 2008, 3(4): 48-50.(inChinese)

    [12] 陳安, 許飛, 閆文偉, 等. 基于DFR法的Al-Li-S4鋁鋰合金鉚接結(jié)構(gòu)[J]. 航空工程進(jìn)展, 2016, 7(1): 101-105.ChenAn,XuFei,YanWenwei,et.al.FatiguereliabilityanalysisofAl-Li-S4rivetedjointsbasedondetailfatiguerating[J].AdvancesinAeronauticalScienceandEngineering, 2016, 7(1): 101-105.(inChinese)

    [13] 方同, 薛璞. 振動(dòng)理論及應(yīng)用[M]. 西安: 西北工業(yè)大學(xué)出版社, 1998: 69-70.FangTong,XuePu.Theoryandapplicationofvibration[M].Xi’an:NorthwesternPolytechnicalUniversityPress, 1998: 69-70.(inChinese)

    [14]MinerMA.Cumulativedamageinfatigue[J].AmericanSocietyofMechanicalEngineers-JournalofAppliedMechanics, 1945, 12(3): 159-164.

    [15] 茍文選. 材料力學(xué)(Ⅱ)[M]. 北京: 科學(xué)出版社, 2010: 346-347.GouWenxuan.Strengthofmaterials(Ⅱ)[M].Beijing:SciencePress, 2010: 346-347.(inChinese)

    (編輯:趙毓梅)

    Effects of Stiffeners’ Cross Section Shape on Thin Walled Structures’Vibration Fatigue Property

    Liu Shuangyan, Tu Yuqian, Miao Yinggang, Li Yulong

    (School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China)

    Stiffeners’ cross section shape plays an important role on the vibration fatigue property of stiffened panels. The vibration tests of two kinds of stiffened panels are carried out. The stiffened panels with L-shaped stiffener and T-shaped stiffener are excited by base motion at the resonance frequency, respectively. The effects of different stiffeners’ shape on the frequency and vibration fatigue life of thin-walled structures are studied experimentally. The possible mechanism is analyzed theoretically and the test process is validated by MSC software. Based on the reasonable FEM simulation, the thin-walled stiffened panels are excited by a sinusoidal distributed load under two boundary conditions(fully clamped and simply supported) in MSC simulation. The effects of stiffeners with six shapes (rectangle, L-shape, T-shape, I-shape, Z-shape and ∏-shape) on the frequency, stress distribution and vibration fatigue life of thin walled structures are discussed. In addition, the predicted study on the life of stiffened panel are carried out. The results show that the first natural frequency of panel with Z-shaped stiffener is the maximum among all the configurations. The peak value of the stress occurred at the same position. When the edges of the panel are fully clamped, the damage position occurred at the two sides of the stiffener and the fatigue life of Z-shaped stiffened panels is the longest. When the edges of the panels are simply supported, the damage position occurred at the center region of the panels and the L-shaped stiffened panels had the longest life.

    thin-walled stiffened structure; cross shape; natural frequency; Von Mises; vibration fatigue life

    2016-12-12;

    2016-12-23

    陜西省科技研究發(fā)展攻關(guān)計(jì)劃 (2009K01-39,2012GY2-26)

    李玉龍,liyulong@nwpu.edu.cn

    1674-8190(2017)02-190-09

    V214.4+2; O346.2

    A

    10.16615/j.cnki.1674-8190.2017.02.011

    劉雙燕(1984-),女,博士研究生。主要研究方向:金屬薄壁結(jié)構(gòu)抑振修補(bǔ)。

    涂玉倩(1992-),女,碩士研究生。主要研究方向:復(fù)合材料薄壁結(jié)構(gòu)的屈曲分析。

    苗應(yīng)剛(1983-),男,博士,講師。主要研究方向:Hopkinson桿技術(shù)、聚合物力學(xué)行為及變形機(jī)理、本構(gòu)方程。

    李玉龍(1961-),男,博士,教授。主要研究方向:飛行器結(jié)構(gòu)抗墜毀設(shè)計(jì)、飛行器結(jié)構(gòu)抗離散源撞擊設(shè)計(jì)、極端環(huán)境下先進(jìn)材料及結(jié)構(gòu)的力學(xué)行為及其優(yōu)化設(shè)計(jì)。

    猜你喜歡
    筋條固支簡(jiǎn)支
    直升機(jī)設(shè)備安裝平臺(tái)的強(qiáng)度優(yōu)化設(shè)計(jì)
    固支結(jié)構(gòu)對(duì)壓阻式壓力傳感器的輸出影響研究
    熱壓彎成形工藝下的高筋蒙皮結(jié)構(gòu)優(yōu)化研究
    復(fù)合材料加筋板高速?zèng)_擊的損傷研究
    基于復(fù)合材料的山體滑坡固支模型及配套傳感器布置方法申請(qǐng)?zhí)? 202010514582.0
    傳感器世界(2021年2期)2021-03-27 19:36:46
    一角點(diǎn)支撐對(duì)面兩邊固支正交各向異性矩形薄板彎曲問題的辛疊加解
    簡(jiǎn)支U型梁橋的抗彎性能分析與試驗(yàn)
    簡(jiǎn)支超寬鋼箱梁設(shè)計(jì)與頂推受力分析
    贛龍鐵路簡(jiǎn)支T梁預(yù)制場(chǎng)建設(shè)
    約束分布對(duì)混合邊界薄板固有頻率影響的研究★
    山西建筑(2018年6期)2018-03-22 07:14:15
    999精品在线视频| 亚洲情色 制服丝袜| 91精品三级在线观看| 午夜免费男女啪啪视频观看| 亚洲欧美一区二区三区国产| 久久精品国产鲁丝片午夜精品| 中文乱码字字幕精品一区二区三区| a 毛片基地| 国产亚洲一区二区精品| 精品一区二区三卡| 黄色一级大片看看| 视频在线观看一区二区三区| 久久久久精品久久久久真实原创| 亚洲欧洲日产国产| 国内精品宾馆在线| 久久久久久久大尺度免费视频| 亚洲国产成人一精品久久久| 在线精品无人区一区二区三| 亚洲欧美精品自产自拍| 亚洲综合色惰| 亚洲精品久久午夜乱码| 最近的中文字幕免费完整| 91成人精品电影| 伊人久久精品亚洲午夜| 亚洲av中文av极速乱| 97超视频在线观看视频| 久久国产亚洲av麻豆专区| 亚洲精品成人av观看孕妇| 久久精品国产亚洲av天美| 在线播放无遮挡| 成人亚洲欧美一区二区av| 国产黄频视频在线观看| 久久午夜综合久久蜜桃| 少妇熟女欧美另类| 午夜免费男女啪啪视频观看| 18在线观看网站| 日韩中文字幕视频在线看片| 一级毛片我不卡| 特大巨黑吊av在线直播| 日韩,欧美,国产一区二区三区| 日本黄大片高清| 国产伦精品一区二区三区视频9| 国产视频首页在线观看| 九色亚洲精品在线播放| 精品人妻在线不人妻| 国产有黄有色有爽视频| 色吧在线观看| 久久午夜福利片| 青春草国产在线视频| 国产淫语在线视频| 伊人亚洲综合成人网| 日韩大片免费观看网站| 一级毛片我不卡| 国产深夜福利视频在线观看| 波野结衣二区三区在线| 精品国产乱码久久久久久小说| 夫妻午夜视频| 亚洲精品第二区| 成人综合一区亚洲| 成年美女黄网站色视频大全免费 | 国产高清有码在线观看视频| 国产精品一二三区在线看| 精品一品国产午夜福利视频| 这个男人来自地球电影免费观看 | 国产色婷婷99| 日韩,欧美,国产一区二区三区| 一级,二级,三级黄色视频| 精品久久久久久电影网| 欧美性感艳星| 国产黄色视频一区二区在线观看| 秋霞伦理黄片| 免费人妻精品一区二区三区视频| 中文字幕最新亚洲高清| 国产69精品久久久久777片| 美女中出高潮动态图| 国产熟女欧美一区二区| 国产精品秋霞免费鲁丝片| 高清不卡的av网站| 国产乱人偷精品视频| 看十八女毛片水多多多| 亚洲欧美成人综合另类久久久| 丰满迷人的少妇在线观看| 最近的中文字幕免费完整| av视频免费观看在线观看| 91久久精品电影网| 国产精品三级大全| 人妻少妇偷人精品九色| 两个人免费观看高清视频| 男女边摸边吃奶| 有码 亚洲区| 精品久久久精品久久久| 中文字幕久久专区| 亚洲精品自拍成人| 视频中文字幕在线观看| 亚洲不卡免费看| 成年人免费黄色播放视频| 日本av免费视频播放| 9色porny在线观看| 在线亚洲精品国产二区图片欧美 | 国产免费视频播放在线视频| 精品人妻一区二区三区麻豆| 亚洲av日韩在线播放| 成人国语在线视频| 久久久久精品性色| kizo精华| 哪个播放器可以免费观看大片| 一级二级三级毛片免费看| 亚洲精品第二区| 日韩成人伦理影院| 亚洲精品,欧美精品| 国产高清国产精品国产三级| 香蕉精品网在线| 欧美+日韩+精品| 一级,二级,三级黄色视频| 中文字幕最新亚洲高清| 亚洲国产欧美日韩在线播放| 免费高清在线观看日韩| 又粗又硬又长又爽又黄的视频| 成人无遮挡网站| 午夜精品国产一区二区电影| 全区人妻精品视频| 免费av中文字幕在线| 最近的中文字幕免费完整| 久久亚洲国产成人精品v| 制服诱惑二区| av国产久精品久网站免费入址| 午夜精品国产一区二区电影| 制服丝袜香蕉在线| 久久ye,这里只有精品| 国产日韩欧美视频二区| 大香蕉久久网| 黑人高潮一二区| 国产精品嫩草影院av在线观看| 91成人精品电影| 满18在线观看网站| 美女xxoo啪啪120秒动态图| 在线免费观看不下载黄p国产| 少妇的逼水好多| 交换朋友夫妻互换小说| 啦啦啦中文免费视频观看日本| 国产无遮挡羞羞视频在线观看| 麻豆成人av视频| 欧美人与性动交α欧美精品济南到 | 另类精品久久| 国产男女超爽视频在线观看| 日韩av免费高清视频| 狂野欧美白嫩少妇大欣赏| 亚洲精华国产精华液的使用体验| 欧美日韩在线观看h| 一区二区三区四区激情视频| 一个人免费看片子| 国产精品国产三级专区第一集| 欧美国产精品一级二级三级| 久久久久久久大尺度免费视频| 女人久久www免费人成看片| 成人国产av品久久久| 在线观看人妻少妇| 亚洲av国产av综合av卡| 毛片一级片免费看久久久久| 秋霞伦理黄片| 26uuu在线亚洲综合色| 王馨瑶露胸无遮挡在线观看| 少妇被粗大的猛进出69影院 | 婷婷色综合大香蕉| 一级毛片我不卡| 精品一品国产午夜福利视频| 26uuu在线亚洲综合色| 午夜福利影视在线免费观看| 国产精品人妻久久久久久| 国产午夜精品久久久久久一区二区三区| 中国国产av一级| 18禁观看日本| 男的添女的下面高潮视频| 狂野欧美激情性xxxx在线观看| 51国产日韩欧美| 97超碰精品成人国产| 69精品国产乱码久久久| 国产精品嫩草影院av在线观看| 男人添女人高潮全过程视频| 男人操女人黄网站| 热99国产精品久久久久久7| 成人午夜精彩视频在线观看| 女性生殖器流出的白浆| 天天操日日干夜夜撸| 18禁动态无遮挡网站| 久久精品久久精品一区二区三区| 中文字幕最新亚洲高清| 精品一区二区三区视频在线| 免费看av在线观看网站| 亚洲av不卡在线观看| 乱人伦中国视频| 赤兔流量卡办理| 亚洲性久久影院| 中文字幕人妻丝袜制服| 汤姆久久久久久久影院中文字幕| 久久久久久久大尺度免费视频| 国产淫语在线视频| 亚洲av国产av综合av卡| 18禁在线播放成人免费| av不卡在线播放| 高清毛片免费看| 18+在线观看网站| 黄色配什么色好看| 国产精品.久久久| 秋霞伦理黄片| 九九在线视频观看精品| 寂寞人妻少妇视频99o| 国产女主播在线喷水免费视频网站| 精品国产露脸久久av麻豆| 亚洲精品视频女| 亚洲精品日韩av片在线观看| 啦啦啦在线观看免费高清www| xxxhd国产人妻xxx| 最近2019中文字幕mv第一页| 蜜桃国产av成人99| kizo精华| 天天躁夜夜躁狠狠久久av| 国产伦精品一区二区三区视频9| 国产av国产精品国产| 男人爽女人下面视频在线观看| 自线自在国产av| 人妻制服诱惑在线中文字幕| 色吧在线观看| 日本黄色片子视频| 久久人妻熟女aⅴ| 青青草视频在线视频观看| 午夜影院在线不卡| 成人国产麻豆网| 99久久精品国产国产毛片| 国产精品不卡视频一区二区| 亚州av有码| 久久久精品区二区三区| 久久人妻熟女aⅴ| 啦啦啦中文免费视频观看日本| 亚洲欧美色中文字幕在线| 日本与韩国留学比较| 大片免费播放器 马上看| 考比视频在线观看| 最近最新中文字幕免费大全7| 久久精品夜色国产| 我的老师免费观看完整版| 久久久久久久精品精品| 久久这里有精品视频免费| 日韩精品有码人妻一区| 国产精品久久久久久精品电影小说| 欧美+日韩+精品| 亚州av有码| 成人二区视频| 久久免费观看电影| 国产在视频线精品| 久久久久久久久久人人人人人人| 国产片内射在线| 伊人久久国产一区二区| 欧美bdsm另类| 狂野欧美白嫩少妇大欣赏| 99视频精品全部免费 在线| 亚洲av电影在线观看一区二区三区| 亚洲不卡免费看| 国产精品 国内视频| 日韩欧美精品免费久久| 亚洲,欧美,日韩| 成人国产av品久久久| 久久婷婷青草| 亚洲av男天堂| 狠狠精品人妻久久久久久综合| 制服丝袜香蕉在线| 啦啦啦视频在线资源免费观看| 国产色婷婷99| 精品国产乱码久久久久久小说| 简卡轻食公司| 日本免费在线观看一区| 亚洲成人手机| 成年人免费黄色播放视频| 伦理电影免费视频| 天天影视国产精品| 七月丁香在线播放| 久久久久久久国产电影| 日本av免费视频播放| 欧美精品亚洲一区二区| 男女免费视频国产| 午夜福利影视在线免费观看| 少妇高潮的动态图| 国产一区二区在线观看日韩| av天堂久久9| 汤姆久久久久久久影院中文字幕| 在线观看免费视频网站a站| 国产成人精品在线电影| 3wmmmm亚洲av在线观看| 91精品国产国语对白视频| 老女人水多毛片| 亚洲成人av在线免费| 亚洲怡红院男人天堂| av女优亚洲男人天堂| av在线老鸭窝| 国产免费福利视频在线观看| 日本爱情动作片www.在线观看| 自线自在国产av| 国产女主播在线喷水免费视频网站| 男女国产视频网站| 久久精品国产亚洲网站| 亚洲精华国产精华液的使用体验| 视频在线观看一区二区三区| 啦啦啦视频在线资源免费观看| 天天影视国产精品| 日本色播在线视频| 日韩中文字幕视频在线看片| 97精品久久久久久久久久精品| 国产精品国产av在线观看| 亚洲怡红院男人天堂| 少妇熟女欧美另类| www.av在线官网国产| 一本—道久久a久久精品蜜桃钙片| 高清不卡的av网站| 黑人巨大精品欧美一区二区蜜桃 | 在线亚洲精品国产二区图片欧美 | av一本久久久久| 一级黄片播放器| 午夜免费观看性视频| 亚洲av中文av极速乱| 青春草亚洲视频在线观看| 飞空精品影院首页| 亚洲精品一二三| 色网站视频免费| 97超视频在线观看视频| 国产日韩一区二区三区精品不卡 | 亚洲无线观看免费| av一本久久久久| 秋霞伦理黄片| a级毛片免费高清观看在线播放| 在线精品无人区一区二区三| 97在线视频观看| 日韩视频在线欧美| 亚洲精品色激情综合| 美女视频免费永久观看网站| 婷婷色综合大香蕉| 成人二区视频| 国产午夜精品久久久久久一区二区三区| 曰老女人黄片| 熟女av电影| 一级片'在线观看视频| 亚洲精品自拍成人| 国产黄色视频一区二区在线观看| 看免费成人av毛片| 久久综合国产亚洲精品| 人妻一区二区av| 激情五月婷婷亚洲| 久久亚洲国产成人精品v| 91久久精品国产一区二区三区| 青春草视频在线免费观看| 亚洲一级一片aⅴ在线观看| 极品少妇高潮喷水抽搐| 免费不卡的大黄色大毛片视频在线观看| 高清在线视频一区二区三区| 草草在线视频免费看| 少妇熟女欧美另类| 亚洲美女黄色视频免费看| 在线免费观看不下载黄p国产| 日韩人妻高清精品专区| 最近的中文字幕免费完整| 日韩人妻高清精品专区| 麻豆乱淫一区二区| 欧美日韩视频精品一区| 亚洲精品国产av蜜桃| 国产精品一国产av| 飞空精品影院首页| 免费播放大片免费观看视频在线观看| 久久精品国产鲁丝片午夜精品| 亚洲国产av影院在线观看| 亚洲av不卡在线观看| 精品久久久久久久久亚洲| 特大巨黑吊av在线直播| 亚洲美女视频黄频| 99热这里只有精品一区| 亚洲怡红院男人天堂| 亚洲av在线观看美女高潮| 日韩制服骚丝袜av| 亚洲国产精品国产精品| 免费不卡的大黄色大毛片视频在线观看| 在线观看人妻少妇| 高清av免费在线| 成人综合一区亚洲| 男女边吃奶边做爰视频| 尾随美女入室| 亚洲,一卡二卡三卡| 国产在线免费精品| 亚洲高清免费不卡视频| 狠狠婷婷综合久久久久久88av| 能在线免费看毛片的网站| 欧美最新免费一区二区三区| av电影中文网址| 熟女电影av网| 国产精品嫩草影院av在线观看| 三级国产精品欧美在线观看| 午夜激情av网站| 精品人妻偷拍中文字幕| 人妻少妇偷人精品九色| 在线观看免费日韩欧美大片 | 欧美日韩精品成人综合77777| 亚洲精品国产av蜜桃| 日本wwww免费看| 久久久久久久久久久丰满| 人人妻人人澡人人看| 国产爽快片一区二区三区| 国产极品粉嫩免费观看在线 | 免费高清在线观看视频在线观看| 国产一级毛片在线| 一区二区三区精品91| 韩国av在线不卡| 国产精品无大码| 亚洲精品中文字幕在线视频| 欧美bdsm另类| 丝袜喷水一区| 黑人猛操日本美女一级片| 国产精品国产av在线观看| 成人18禁高潮啪啪吃奶动态图 | 丰满乱子伦码专区| 精品少妇黑人巨大在线播放| 亚洲婷婷狠狠爱综合网| 大话2 男鬼变身卡| 麻豆乱淫一区二区| 国产一区二区在线观看日韩| 丁香六月天网| 国产精品欧美亚洲77777| 久久毛片免费看一区二区三区| av天堂久久9| 国产成人freesex在线| 久久国产亚洲av麻豆专区| 精品国产一区二区三区久久久樱花| 久久国内精品自在自线图片| 久久久久久久久久久久大奶| 国产日韩欧美视频二区| 大香蕉久久成人网| 欧美日韩亚洲高清精品| 国产精品一国产av| 99热这里只有精品一区| 成人国产av品久久久| 人人妻人人澡人人看| 免费黄网站久久成人精品| 色吧在线观看| 考比视频在线观看| av在线播放精品| 黑丝袜美女国产一区| 最后的刺客免费高清国语| 我的老师免费观看完整版| 韩国av在线不卡| 观看av在线不卡| 三上悠亚av全集在线观看| 一级毛片 在线播放| 日韩视频在线欧美| 狂野欧美激情性xxxx在线观看| 国产又色又爽无遮挡免| 久久99热这里只频精品6学生| 成人免费观看视频高清| 91精品国产国语对白视频| 中文字幕人妻丝袜制服| 最近2019中文字幕mv第一页| 国产精品一国产av| 乱码一卡2卡4卡精品| 国产欧美另类精品又又久久亚洲欧美| 国产又色又爽无遮挡免| 国产成人免费无遮挡视频| 精品国产露脸久久av麻豆| 下体分泌物呈黄色| av免费观看日本| 亚洲av在线观看美女高潮| 亚洲综合色网址| 丝袜喷水一区| 久久久国产一区二区| 少妇被粗大的猛进出69影院 | 成人无遮挡网站| 成人毛片a级毛片在线播放| 成人18禁高潮啪啪吃奶动态图 | 蜜臀久久99精品久久宅男| 婷婷色av中文字幕| 少妇丰满av| 热re99久久精品国产66热6| 亚洲精品,欧美精品| 欧美bdsm另类| av黄色大香蕉| 男女高潮啪啪啪动态图| 亚洲精品美女久久av网站| 一区二区三区免费毛片| 日韩人妻高清精品专区| 香蕉精品网在线| 十分钟在线观看高清视频www| 国产av一区二区精品久久| 啦啦啦中文免费视频观看日本| 在线天堂最新版资源| 国产 精品1| 69精品国产乱码久久久| 日本欧美视频一区| 色视频在线一区二区三区| 人体艺术视频欧美日本| 免费黄色在线免费观看| av福利片在线| 国产日韩欧美在线精品| 国产日韩一区二区三区精品不卡 | 色视频在线一区二区三区| 人体艺术视频欧美日本| 韩国av在线不卡| 国产伦精品一区二区三区视频9| 久久精品熟女亚洲av麻豆精品| 日本av手机在线免费观看| 一级毛片黄色毛片免费观看视频| 国产精品国产三级国产专区5o| 国产精品秋霞免费鲁丝片| 精品一品国产午夜福利视频| 久久97久久精品| 男女免费视频国产| 香蕉精品网在线| 五月玫瑰六月丁香| 一个人看视频在线观看www免费| 好男人视频免费观看在线| 久久人妻熟女aⅴ| 性色av一级| 亚洲欧美成人综合另类久久久| 在线观看免费视频网站a站| 国产高清有码在线观看视频| 中文字幕人妻丝袜制服| 热re99久久精品国产66热6| 我的女老师完整版在线观看| 日日摸夜夜添夜夜添av毛片| 好男人视频免费观看在线| 久久热精品热| 亚洲精品456在线播放app| 亚洲精品色激情综合| 高清视频免费观看一区二区| 免费观看的影片在线观看| 2022亚洲国产成人精品| 中文字幕人妻熟人妻熟丝袜美| 成人影院久久| 国产av码专区亚洲av| 人妻系列 视频| 桃花免费在线播放| 亚洲精品一二三| 午夜老司机福利剧场| 超色免费av| 丝袜喷水一区| 日本av手机在线免费观看| 国产免费福利视频在线观看| 如何舔出高潮| 精品久久国产蜜桃| 波野结衣二区三区在线| 日韩一区二区三区影片| 成人国语在线视频| 一个人看视频在线观看www免费| 国产精品嫩草影院av在线观看| 国产熟女午夜一区二区三区 | 天堂8中文在线网| 全区人妻精品视频| 欧美亚洲日本最大视频资源| 久久精品国产亚洲网站| 超色免费av| 纯流量卡能插随身wifi吗| 夜夜爽夜夜爽视频| 欧美日韩成人在线一区二区| 欧美精品亚洲一区二区| 汤姆久久久久久久影院中文字幕| 另类精品久久| 日本wwww免费看| 丁香六月天网| 日韩av免费高清视频| 亚洲不卡免费看| 精品人妻熟女毛片av久久网站| 久久韩国三级中文字幕| 国产成人精品福利久久| av不卡在线播放| av在线观看视频网站免费| 高清毛片免费看| 亚洲精品色激情综合| 亚洲精华国产精华液的使用体验| 久久 成人 亚洲| 国产av精品麻豆| 色婷婷久久久亚洲欧美| a级毛片在线看网站| 国产69精品久久久久777片| 中文欧美无线码| 久久久久久人妻| 国产一区亚洲一区在线观看| 大又大粗又爽又黄少妇毛片口| 亚洲欧洲日产国产| 国产色婷婷99| 观看美女的网站| 精品国产一区二区三区久久久樱花| 精品人妻熟女av久视频| 免费观看a级毛片全部| 丝瓜视频免费看黄片| 777米奇影视久久| 国产一区有黄有色的免费视频| 男女国产视频网站| 久久鲁丝午夜福利片| 国产免费一区二区三区四区乱码| 男女免费视频国产| 午夜激情福利司机影院| 亚洲美女黄色视频免费看| 久热久热在线精品观看| 亚洲欧美成人综合另类久久久| 亚洲国产精品国产精品| 777米奇影视久久| 日韩成人伦理影院| 久久 成人 亚洲| 在线天堂最新版资源| 中文字幕av电影在线播放| 精品熟女少妇av免费看| 三级国产精品片| 色婷婷av一区二区三区视频| 精品人妻在线不人妻| 丰满乱子伦码专区| 国产精品一区二区三区四区免费观看| 人妻人人澡人人爽人人| 91午夜精品亚洲一区二区三区| 婷婷色av中文字幕| 日韩不卡一区二区三区视频在线| 日本91视频免费播放| 狂野欧美激情性xxxx在线观看|