• <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
    黄片小视频在线播放| 日韩视频一区二区在线观看| 亚洲中文av在线| kizo精华| 国产三级黄色录像| 成年人黄色毛片网站| 久久九九热精品免费| 后天国语完整版免费观看| 国产麻豆69| 黄色a级毛片大全视频| av免费在线观看网站| 韩国高清视频一区二区三区| 亚洲,欧美精品.| 性高湖久久久久久久久免费观看| 欧美日韩一级在线毛片| 精品少妇久久久久久888优播| 久久久久国产精品人妻一区二区| 在线观看免费视频网站a站| 国产一卡二卡三卡精品| 青草久久国产| 天堂8中文在线网| 亚洲av片天天在线观看| 成人国产一区最新在线观看| 十八禁网站免费在线| 欧美亚洲 丝袜 人妻 在线| 欧美日韩亚洲国产一区二区在线观看 | 色婷婷久久久亚洲欧美| 三上悠亚av全集在线观看| 1024香蕉在线观看| 欧美日韩视频精品一区| 脱女人内裤的视频| 国产精品国产三级国产专区5o| 精品一区二区三卡| 亚洲欧美清纯卡通| 亚洲国产看品久久| 日韩人妻精品一区2区三区| 欧美精品av麻豆av| 80岁老熟妇乱子伦牲交| 一区福利在线观看| 欧美老熟妇乱子伦牲交| 午夜福利,免费看| 99热国产这里只有精品6| 久久九九热精品免费| 18禁黄网站禁片午夜丰满| 女警被强在线播放| 两人在一起打扑克的视频| 久久久国产一区二区| 正在播放国产对白刺激| 我的亚洲天堂| 搡老熟女国产l中国老女人| 免费在线观看日本一区| 日韩电影二区| 亚洲av成人不卡在线观看播放网 | 丝袜美腿诱惑在线| 大香蕉久久网| 少妇被粗大的猛进出69影院| 国产精品久久久av美女十八| 好男人电影高清在线观看| 在线观看一区二区三区激情| 国产不卡av网站在线观看| 十八禁网站免费在线| 国产福利在线免费观看视频| 国产精品欧美亚洲77777| 成年人午夜在线观看视频| 搡老岳熟女国产| 亚洲欧美一区二区三区黑人| 午夜久久久在线观看| 国产老妇伦熟女老妇高清| 国产又爽黄色视频| 国产福利在线免费观看视频| 亚洲精品国产一区二区精华液| 中文字幕精品免费在线观看视频| 久久狼人影院| 这个男人来自地球电影免费观看| 色视频在线一区二区三区| 国产成人欧美| 成在线人永久免费视频| 老司机午夜福利在线观看视频 | 又紧又爽又黄一区二区| 欧美黑人欧美精品刺激| 免费在线观看黄色视频的| 日日摸夜夜添夜夜添小说| 黄色 视频免费看| 99国产精品免费福利视频| 欧美成狂野欧美在线观看| 91老司机精品| 好男人电影高清在线观看| 国产免费现黄频在线看| 操美女的视频在线观看| 久久久久网色| 宅男免费午夜| 别揉我奶头~嗯~啊~动态视频 | 久久人妻熟女aⅴ| 成人国产av品久久久| 另类精品久久| 91成年电影在线观看| 午夜老司机福利片| 一本久久精品| 欧美黑人精品巨大| 国产av又大| 久久99热这里只频精品6学生| 老汉色∧v一级毛片| www.999成人在线观看| 美女主播在线视频| 一区二区三区四区激情视频| 新久久久久国产一级毛片| 欧美日韩亚洲综合一区二区三区_| 亚洲免费av在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 性色av乱码一区二区三区2| 久久久久久亚洲精品国产蜜桃av| 老熟妇乱子伦视频在线观看 | 国产精品亚洲av一区麻豆| 一本大道久久a久久精品| 亚洲欧美色中文字幕在线| 亚洲av日韩在线播放| 日本猛色少妇xxxxx猛交久久| 久久国产精品影院| 91麻豆精品激情在线观看国产 | 女人高潮潮喷娇喘18禁视频| 日韩人妻精品一区2区三区| av天堂久久9| 啦啦啦啦在线视频资源| 国产免费现黄频在线看| 亚洲国产中文字幕在线视频| 97精品久久久久久久久久精品| 午夜福利,免费看| 午夜精品久久久久久毛片777| 日日摸夜夜添夜夜添小说| 国产在线一区二区三区精| 精品国产一区二区三区久久久樱花| 在线av久久热| 动漫黄色视频在线观看| 亚洲九九香蕉| 精品一区在线观看国产| 欧美亚洲日本最大视频资源| 满18在线观看网站| av天堂久久9| 老司机靠b影院| videosex国产| 国产成人啪精品午夜网站| 欧美成人午夜精品| 精品国产一区二区久久| 一边摸一边抽搐一进一出视频| 国产又色又爽无遮挡免| 亚洲欧美成人综合另类久久久| 女人爽到高潮嗷嗷叫在线视频| 搡老岳熟女国产| 欧美亚洲 丝袜 人妻 在线| 一区二区三区四区激情视频| 男人添女人高潮全过程视频| 69精品国产乱码久久久| 亚洲成人手机| 精品人妻在线不人妻| 成年人黄色毛片网站| 精品高清国产在线一区| 在线永久观看黄色视频| 国产成人av教育| 夜夜骑夜夜射夜夜干| 亚洲精品久久成人aⅴ小说| 两性夫妻黄色片| 免费一级毛片在线播放高清视频 | www日本在线高清视频| 欧美日韩黄片免| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成av片中文字幕在线观看| 欧美精品一区二区大全| 日韩制服骚丝袜av| 人人妻人人添人人爽欧美一区卜| 一边摸一边抽搐一进一出视频| 极品少妇高潮喷水抽搐| 日韩,欧美,国产一区二区三区| 久久免费观看电影| 性少妇av在线| 免费在线观看日本一区| 97精品久久久久久久久久精品| 精品一区二区三卡| 亚洲一码二码三码区别大吗| 99国产精品99久久久久| 国产在线观看jvid| 香蕉国产在线看| 黑人操中国人逼视频| 久久国产精品男人的天堂亚洲| 十八禁人妻一区二区| 窝窝影院91人妻| 日本91视频免费播放| 免费看十八禁软件| 久久精品国产综合久久久| 考比视频在线观看| 日韩有码中文字幕| 欧美性长视频在线观看| 日韩,欧美,国产一区二区三区| 日韩大片免费观看网站| 亚洲精品第二区| 国产精品.久久久| 国产97色在线日韩免费| 蜜桃国产av成人99| 日本黄色日本黄色录像| 国产熟女午夜一区二区三区| 亚洲av男天堂| 女警被强在线播放| 免费一级毛片在线播放高清视频 | 国产片内射在线| 99久久综合免费| 国产在线一区二区三区精| 久久人妻福利社区极品人妻图片| 精品久久久久久电影网| 欧美日韩国产mv在线观看视频| 亚洲国产毛片av蜜桃av| 男女之事视频高清在线观看| 亚洲avbb在线观看| 久久国产精品影院| kizo精华| 欧美人与性动交α欧美精品济南到| 高潮久久久久久久久久久不卡| 老汉色∧v一级毛片| 久久国产精品男人的天堂亚洲| 国内毛片毛片毛片毛片毛片| 国产麻豆69| 国产精品99久久99久久久不卡| 日本黄色日本黄色录像| 一级毛片电影观看| 丝袜喷水一区| 国产熟女午夜一区二区三区| 又紧又爽又黄一区二区| 老司机影院毛片| 精品第一国产精品| 精品一区二区三卡| 免费不卡黄色视频| 18禁观看日本| 成人手机av| 欧美乱码精品一区二区三区| 亚洲中文av在线| 欧美中文综合在线视频| 新久久久久国产一级毛片| 亚洲精品日韩在线中文字幕| 老司机深夜福利视频在线观看 | 亚洲欧美清纯卡通| 亚洲精品美女久久久久99蜜臀| 国产亚洲欧美精品永久| 国产日韩欧美亚洲二区| 中国美女看黄片| 男男h啪啪无遮挡| 国产91精品成人一区二区三区 | 超碰成人久久| 国产成人精品在线电影| 国精品久久久久久国模美| av欧美777| 亚洲欧美成人综合另类久久久| 久热爱精品视频在线9| 国产有黄有色有爽视频| av国产精品久久久久影院| 国产真人三级小视频在线观看| 国产熟女午夜一区二区三区| 侵犯人妻中文字幕一二三四区| 俄罗斯特黄特色一大片| 啦啦啦在线免费观看视频4| 男女无遮挡免费网站观看| 女警被强在线播放| 久久人人爽人人片av| 欧美激情高清一区二区三区| 国产成人a∨麻豆精品| 亚洲,欧美精品.| 可以免费在线观看a视频的电影网站| 不卡一级毛片| 亚洲激情五月婷婷啪啪| 欧美性长视频在线观看| 岛国毛片在线播放| 欧美精品av麻豆av| h视频一区二区三区| 99久久精品国产亚洲精品| 黄网站色视频无遮挡免费观看| 男女之事视频高清在线观看| 午夜91福利影院| 精品人妻熟女毛片av久久网站| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久大尺度免费视频| 亚洲av欧美aⅴ国产| 在线观看舔阴道视频| 美女高潮喷水抽搐中文字幕| 夫妻午夜视频| 精品人妻一区二区三区麻豆| 国产精品秋霞免费鲁丝片| 天堂俺去俺来也www色官网| 免费观看人在逋| 精品久久久久久电影网| 欧美国产精品va在线观看不卡| 国产精品九九99| 正在播放国产对白刺激| 精品一区二区三卡| 爱豆传媒免费全集在线观看| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区三区av在线| 欧美人与性动交α欧美软件| 老司机影院毛片| 精品国产乱子伦一区二区三区 | 色视频在线一区二区三区| 一区福利在线观看| 成人三级做爰电影| 热re99久久精品国产66热6| 亚洲av电影在线进入| 久久99热这里只频精品6学生| 国产av一区二区精品久久| 99热国产这里只有精品6| 亚洲欧美清纯卡通| 在线亚洲精品国产二区图片欧美| av电影中文网址| 精品亚洲成a人片在线观看| 亚洲国产日韩一区二区| 午夜老司机福利片| 99九九在线精品视频| 天天影视国产精品| 电影成人av| 黑人欧美特级aaaaaa片| 亚洲国产成人一精品久久久| 90打野战视频偷拍视频| 99热网站在线观看| 曰老女人黄片| 亚洲第一青青草原| 亚洲熟女精品中文字幕| 黄色毛片三级朝国网站| 男女午夜视频在线观看| 男女下面插进去视频免费观看| 多毛熟女@视频| 久久久久久人人人人人| 国产男女超爽视频在线观看| 伦理电影免费视频| 青春草视频在线免费观看| av片东京热男人的天堂| 97在线人人人人妻| 99久久精品国产亚洲精品| 午夜视频精品福利| 中文字幕人妻熟女乱码| 老汉色∧v一级毛片| 嫁个100分男人电影在线观看| 精品福利观看| 纵有疾风起免费观看全集完整版| 成年人午夜在线观看视频| 天堂中文最新版在线下载| 精品国产国语对白av| 国产在线视频一区二区| 脱女人内裤的视频| 女人久久www免费人成看片| 人人妻人人添人人爽欧美一区卜| 欧美黑人欧美精品刺激| av国产精品久久久久影院| 午夜福利,免费看| 国产高清videossex| 久久精品国产亚洲av香蕉五月 | 欧美xxⅹ黑人| 99国产精品一区二区三区| 少妇裸体淫交视频免费看高清 | 日韩熟女老妇一区二区性免费视频| 50天的宝宝边吃奶边哭怎么回事| 老汉色av国产亚洲站长工具| 高潮久久久久久久久久久不卡| 亚洲第一青青草原| 人人妻人人添人人爽欧美一区卜| 十八禁网站免费在线| 精品第一国产精品| 高清视频免费观看一区二区| 精品国产乱子伦一区二区三区 | 亚洲欧美一区二区三区久久| 免费高清在线观看日韩| 91九色精品人成在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 肉色欧美久久久久久久蜜桃| 可以免费在线观看a视频的电影网站| 丝袜脚勾引网站| 搡老岳熟女国产| 如日韩欧美国产精品一区二区三区| 97精品久久久久久久久久精品| 亚洲av欧美aⅴ国产| 9热在线视频观看99| 亚洲欧美一区二区三区黑人| 国产精品熟女久久久久浪| 黄频高清免费视频| 女性被躁到高潮视频| 精品免费久久久久久久清纯 | 黄色a级毛片大全视频| 亚洲精品一二三| 国产高清国产精品国产三级| 欧美激情 高清一区二区三区| 精品一区在线观看国产| 极品少妇高潮喷水抽搐| 18禁裸乳无遮挡动漫免费视频| 欧美精品高潮呻吟av久久| 一本—道久久a久久精品蜜桃钙片| 乱人伦中国视频| 精品一区二区三区四区五区乱码| 十八禁高潮呻吟视频| 国产区一区二久久| av在线老鸭窝| 高清av免费在线| 巨乳人妻的诱惑在线观看| 高清av免费在线| 一区福利在线观看| 正在播放国产对白刺激| 一本一本久久a久久精品综合妖精| 精品一区二区三区av网在线观看 | 中文字幕另类日韩欧美亚洲嫩草| 欧美 亚洲 国产 日韩一| 亚洲欧美精品自产自拍| 亚洲一区中文字幕在线| 这个男人来自地球电影免费观看| 动漫黄色视频在线观看| 最黄视频免费看| 亚洲专区中文字幕在线| 色婷婷av一区二区三区视频| 99国产极品粉嫩在线观看| 精品国产一区二区久久| 丁香六月天网| 国产精品免费大片| av欧美777| 国产精品欧美亚洲77777| 国产91精品成人一区二区三区 | 亚洲国产欧美日韩在线播放| 日韩一卡2卡3卡4卡2021年| 丝袜喷水一区| 欧美少妇被猛烈插入视频| 亚洲国产精品999| 日本av手机在线免费观看| 日韩三级视频一区二区三区| 亚洲国产欧美网| 俄罗斯特黄特色一大片| 欧美激情久久久久久爽电影 | 久久综合国产亚洲精品| 十八禁网站网址无遮挡| 国产成人精品久久二区二区免费| 少妇裸体淫交视频免费看高清 | 亚洲精品在线美女| 亚洲人成电影观看| 12—13女人毛片做爰片一| 日本一区二区免费在线视频| av天堂在线播放| av又黄又爽大尺度在线免费看| 精品国产一区二区三区四区第35| 欧美国产精品一级二级三级| 欧美黑人欧美精品刺激| 9191精品国产免费久久| 黄频高清免费视频| 国产精品免费视频内射| 国产欧美日韩一区二区精品| 成年av动漫网址| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美+亚洲+日韩+国产| 色婷婷av一区二区三区视频| 久久国产精品男人的天堂亚洲| 高清在线国产一区| 国产欧美日韩精品亚洲av| 国产1区2区3区精品| 99精国产麻豆久久婷婷| 丝袜脚勾引网站| 免费观看av网站的网址| 黄色视频在线播放观看不卡| 嫁个100分男人电影在线观看| 老熟妇仑乱视频hdxx| 日本wwww免费看| 精品亚洲成a人片在线观看| 中文精品一卡2卡3卡4更新| 午夜福利一区二区在线看| 中国国产av一级| 久久久水蜜桃国产精品网| 99国产精品99久久久久| 51午夜福利影视在线观看| 少妇 在线观看| 如日韩欧美国产精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 成年女人毛片免费观看观看9 | 国产日韩欧美视频二区| 免费在线观看黄色视频的| 人人妻,人人澡人人爽秒播| 亚洲成国产人片在线观看| 成年动漫av网址| 久久国产精品人妻蜜桃| 女性生殖器流出的白浆| 99国产精品一区二区三区| 又大又爽又粗| 亚洲专区中文字幕在线| 精品久久蜜臀av无| 亚洲国产精品999| 午夜福利,免费看| 精品少妇久久久久久888优播| 亚洲国产日韩一区二区| 欧美国产精品va在线观看不卡| 大陆偷拍与自拍| 亚洲精品av麻豆狂野| 18禁观看日本| 久久久久精品人妻al黑| 正在播放国产对白刺激| 在线永久观看黄色视频| 国产精品欧美亚洲77777| 日韩欧美一区二区三区在线观看 | 国产精品99久久99久久久不卡| 中文字幕人妻丝袜制服| 99国产精品99久久久久| 性高湖久久久久久久久免费观看| 青春草亚洲视频在线观看| 人妻 亚洲 视频| 免费久久久久久久精品成人欧美视频| 国产日韩欧美在线精品| 女人精品久久久久毛片| 日韩制服丝袜自拍偷拍| 亚洲国产av新网站| 在线亚洲精品国产二区图片欧美| 亚洲国产av新网站| 一区二区日韩欧美中文字幕| 一级毛片女人18水好多| 免费观看av网站的网址| 黄网站色视频无遮挡免费观看| 国产一区二区在线观看av| 青春草视频在线免费观看| 久久综合国产亚洲精品| 肉色欧美久久久久久久蜜桃| 国产日韩欧美视频二区| 日韩制服骚丝袜av| 两个人免费观看高清视频| 久久 成人 亚洲| 亚洲欧美一区二区三区黑人| 热re99久久精品国产66热6| 日韩制服骚丝袜av| 久久天堂一区二区三区四区| 亚洲中文日韩欧美视频| 老鸭窝网址在线观看| 999久久久国产精品视频| 中文字幕色久视频| 亚洲成人国产一区在线观看| 国产在视频线精品| 桃红色精品国产亚洲av| 精品亚洲乱码少妇综合久久| 男女国产视频网站| 大陆偷拍与自拍| 一级黄色大片毛片| 18禁黄网站禁片午夜丰满| 美女中出高潮动态图| 一级毛片女人18水好多| 狠狠狠狠99中文字幕| 999久久久精品免费观看国产| 色婷婷av一区二区三区视频| a级片在线免费高清观看视频| 国产成人a∨麻豆精品| 欧美精品人与动牲交sv欧美| 亚洲精品国产av蜜桃| 国产成人啪精品午夜网站| 19禁男女啪啪无遮挡网站| 老司机深夜福利视频在线观看 | 涩涩av久久男人的天堂| av天堂久久9| 天天添夜夜摸| 国产视频一区二区在线看| 老司机午夜十八禁免费视频| 午夜福利乱码中文字幕| 窝窝影院91人妻| 伊人亚洲综合成人网| 亚洲精品国产av蜜桃| 国产精品免费大片| 国产精品.久久久| 97精品久久久久久久久久精品| 一区二区av电影网| 少妇被粗大的猛进出69影院| 王馨瑶露胸无遮挡在线观看| 桃红色精品国产亚洲av| 无遮挡黄片免费观看| 久久久久久久久免费视频了| 在线看a的网站| 女人被躁到高潮嗷嗷叫费观| 亚洲avbb在线观看| 国产精品一区二区在线观看99| 国产精品免费视频内射| 亚洲国产av影院在线观看| 欧美变态另类bdsm刘玥| av电影中文网址| 久久九九热精品免费| 国产在线观看jvid| 久久青草综合色| 日韩人妻精品一区2区三区| 亚洲av成人不卡在线观看播放网 | 国产三级黄色录像| 九色亚洲精品在线播放| 精品一区二区三区四区五区乱码| 国产又色又爽无遮挡免| 一本久久精品| 亚洲一卡2卡3卡4卡5卡精品中文| 成年动漫av网址| 精品少妇内射三级| 久久性视频一级片| 国产成人系列免费观看| 亚洲欧美日韩高清在线视频 | 久久久精品94久久精品| 黑丝袜美女国产一区| 妹子高潮喷水视频| 久久久久国产一级毛片高清牌| 国产欧美日韩精品亚洲av| a级毛片黄视频| 日韩一卡2卡3卡4卡2021年| 亚洲第一青青草原| 热99国产精品久久久久久7| 99久久国产精品久久久| 亚洲精品一二三| 亚洲自偷自拍图片 自拍| 天天操日日干夜夜撸| 国产视频一区二区在线看| 久久久国产成人免费| 乱人伦中国视频| 美女扒开内裤让男人捅视频| 欧美黑人欧美精品刺激| 一本综合久久免费| 高清黄色对白视频在线免费看| 丰满饥渴人妻一区二区三| 欧美日韩一级在线毛片| 亚洲第一av免费看| 亚洲欧洲精品一区二区精品久久久| 欧美日本中文国产一区发布| 精品国产国语对白av| 天天操日日干夜夜撸|