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

    誘導(dǎo)分布對(duì)高速列車多胞結(jié)構(gòu)耐撞性的影響

    2022-06-26 00:40:10梁習(xí)鋒陳劼昊李治祥馬聞許平
    關(guān)鍵詞:褶皺數(shù)量有限元

    梁習(xí)鋒,陳劼昊,李治祥,馬聞,許平

    (1.中南大學(xué)交通運(yùn)輸工程學(xué)院,軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室,湖南長沙,410075;2.中南大學(xué)軌道交通安全關(guān)鍵技術(shù)國際聯(lián)合研究實(shí)驗(yàn)室,湖南長沙,410075;3.中南大學(xué)鐵道車輛安全技術(shù)國家與地方聯(lián)合工程研究中心,湖南長沙,410075)

    高速列車碰撞事故容易造成大量人員傷亡和重大經(jīng)濟(jì)損失,因此,需要提高被動(dòng)安全防護(hù)性能來降低事故損失,而被動(dòng)安全防護(hù)主要依靠列車端部的吸能結(jié)構(gòu)將沖擊動(dòng)能轉(zhuǎn)化為變形吸能[1]。薄壁多胞結(jié)構(gòu)以質(zhì)量小、吸能效率高以及制作成本低而被廣泛應(yīng)用于汽車、列車和航空等交通運(yùn)輸設(shè)備的吸能結(jié)構(gòu)中[2-4]。薄壁結(jié)構(gòu)性能可以通過耐撞性指標(biāo)進(jìn)行評(píng)估。研究人員對(duì)如何增強(qiáng)吸能結(jié)構(gòu)的耐撞性進(jìn)行了大量研究[5-7],以確保設(shè)計(jì)的薄壁結(jié)構(gòu)滿足被動(dòng)安全防護(hù)需求。比較常用的耐撞性評(píng)價(jià)指標(biāo)有吸能量EA、峰值力FPC、比吸能ESA、壓潰力效率ECL和平均破碎力FMC,其中,ESA和FPC分別作為結(jié)構(gòu)使用效率和被動(dòng)安全防護(hù)性能的重要參考指標(biāo)而被重點(diǎn)研究[8]。降低吸能結(jié)構(gòu)的FPC能夠有效減小碰撞中的減速度,在結(jié)構(gòu)中加入開孔能有效降低吸能結(jié)構(gòu)的

    GUPTA 等[10]在圓管中加入圓形誘導(dǎo)孔,發(fā)現(xiàn)能夠有效避免結(jié)構(gòu)的整體彎曲,并且使結(jié)構(gòu)的FPC降低。LI等[11]對(duì)泡沫填充方管的4個(gè)角和4條邊上的誘導(dǎo)孔進(jìn)行研究,通過有限元計(jì)算對(duì)比了不同形狀的誘導(dǎo)孔對(duì)于結(jié)構(gòu)耐撞性的影響。NIKKHAH等[12]采用在方管和圓管的管壁中加入方形、圓形和六邊形的誘導(dǎo)設(shè)計(jì)方案,研究了不同形狀的開孔設(shè)計(jì)對(duì)于單胞管的影響,此外,還將不同數(shù)量的誘導(dǎo)孔均勻分布在結(jié)構(gòu)的管壁中,并通過優(yōu)化設(shè)計(jì)確定了最佳數(shù)量配置。姚如洋等[13]對(duì)開孔圓管在軸向載荷下的變形和吸能量進(jìn)行了試驗(yàn)和理論研究,得到了單層孔數(shù)量和多層孔數(shù)量對(duì)于結(jié)構(gòu)耐撞性的影響。而對(duì)于多胞結(jié)構(gòu),其研究思路與單胞開孔結(jié)構(gòu)的研究思路基本相同。PIRMOHAMMAD 等[14]研究了誘導(dǎo)孔對(duì)方形和八邊形雙滾筒結(jié)構(gòu)的影響,其中誘導(dǎo)孔的形狀和大小對(duì)雙滾筒結(jié)構(gòu)的耐撞性的影響明顯。TRAN等[15]在多層方形嵌套管中加入開孔設(shè)計(jì),在每一層的方管端部都增加圓形孔來提高結(jié)構(gòu)的耐撞性,研究了嵌套管對(duì)結(jié)構(gòu)耐撞性的影響。ESA 等[16]在圓形嵌套管的每一層圓管添加誘導(dǎo)孔,提高了吸能結(jié)構(gòu)的耐撞性。

    目前,人們對(duì)于誘導(dǎo)孔的研究多集中于簡單管結(jié)構(gòu)方面,對(duì)于多胞管中的開孔方案有待進(jìn)一步研究。高速動(dòng)車組前端多胞吸能結(jié)構(gòu)作為主要的吸能部件,需要吸收列車碰撞過程中大部分動(dòng)能,提高其耐撞性是主要的設(shè)計(jì)目標(biāo)。本文以文獻(xiàn)[17-18]中的多胞結(jié)構(gòu)為研究對(duì)象,根據(jù)其研究的多胞吸能結(jié)構(gòu)形狀,將誘導(dǎo)孔加入其中,分析在不同胞壁上開孔對(duì)于結(jié)構(gòu)耐撞性的影響。

    1 開孔孔結(jié)構(gòu)的幾何結(jié)構(gòu)和有限元模型

    1.1 幾何結(jié)構(gòu)

    本文研究的多胞吸能結(jié)構(gòu)橫截面形狀如圖1所示。該結(jié)構(gòu)是由4個(gè)正六邊形和1個(gè)八邊形組成的五胞元結(jié)構(gòu),正六邊形的邊長為56 mm,八邊形上下2 條邊長度為56 mm,而左右2 條邊長度為 51 mm,整個(gè)結(jié)構(gòu)長度為550 mm,壁厚為5 mm。

    圖1 多胞開孔結(jié)構(gòu)橫截面示意圖Fig.1 Cross section view of multi-cell windowed structure

    通過在多胞結(jié)構(gòu)不同的胞壁上添加圓形誘導(dǎo)孔研究誘導(dǎo)孔對(duì)于多胞結(jié)構(gòu)耐撞性的影響。在結(jié)構(gòu)端部開一層孔,具體位置位于結(jié)構(gòu)端部第1個(gè)折疊位置,如圖2(a)所示。在本文研究的結(jié)構(gòu)中,開孔位置設(shè)置在第1個(gè)褶皺的折彎處,根據(jù)無孔結(jié)構(gòu)的變形模式,將開孔的位置設(shè)置在距離端部40 mm處,在整個(gè)結(jié)構(gòu)中,開孔位置如圖2(b)所示,誘導(dǎo)孔的半徑為15 mm。

    圖2 多胞結(jié)構(gòu)開孔設(shè)計(jì)Fig.2 Windowed design of multi-cell structure

    在確定好開孔位置后,在不同的胞壁上引入誘導(dǎo)孔,形成如圖3 所示的6 種誘導(dǎo)結(jié)構(gòu),其中,紅線為加入圓形誘導(dǎo)孔的邊,多胞結(jié)構(gòu)中八邊形胞元位于內(nèi)側(cè),六邊形胞元位于外側(cè),通過結(jié)構(gòu)中心點(diǎn)到開孔邊的距離定義誘導(dǎo)孔分布規(guī)律,根據(jù)距離由近到遠(yuǎn)變化將誘導(dǎo)孔位置分成內(nèi)側(cè)誘導(dǎo)孔、中間誘導(dǎo)孔和外側(cè)誘導(dǎo)孔,而根據(jù)誘導(dǎo)孔對(duì)稱分布原則,采用4個(gè)誘導(dǎo)孔和8個(gè)誘導(dǎo)孔能夠囊括所有邊上的誘導(dǎo)孔分布情況。根據(jù)不同的開孔數(shù)量和位置,分布將這6 種結(jié)構(gòu)命名為S4-1,S4-2,S4-3,S8-1,S8-2和S8-3。

    圖3 多胞結(jié)構(gòu)開孔方案說明Fig.3 Description of multi-cell structure window scheme

    1.2 有限元模型和驗(yàn)證

    采用非線性有限元軟件LS-DYNA對(duì)列車端部多胞結(jié)構(gòu)在沖擊載荷下的力學(xué)行為進(jìn)行分析,有限元模型見圖4。

    圖4 多胞結(jié)構(gòu)沖擊行為有限元模型Fig.4 Finite element model of multi-cell structure in impact behavior

    1)多胞吸能結(jié)構(gòu)。為了平衡好計(jì)算時(shí)間和結(jié)果的精確度,采用邊長為5 mm[17]的Beltschko-Tsay四節(jié)點(diǎn)方形殼單元對(duì)多胞吸能結(jié)構(gòu)過程進(jìn)行模擬。多胞結(jié)構(gòu)的前端板和后端板采用邊長為5 mm的正六面體實(shí)體單元進(jìn)行建模。吸能結(jié)構(gòu)由6008 鋁合金構(gòu)成,采用彈塑性本構(gòu)模型建模,材料參數(shù)見文獻(xiàn)[18],彈塑性材料使用MAT_24來定義。

    2)試驗(yàn)臺(tái)車模型采用邊長為40 mm 的方形殼單元進(jìn)行建模,以節(jié)省計(jì)算時(shí)間。

    3)軌道模型采用邊長為40 mm 的方形殼單元進(jìn)行建模。

    4)固定剛性墻。在整個(gè)沖擊過程中,使用自動(dòng)單面接觸算法模擬多胞結(jié)構(gòu)自身變形產(chǎn)生的接觸,采用自動(dòng)面對(duì)面接觸算法模擬多胞結(jié)構(gòu)、鐵軌、臺(tái)車和剛性墻之間的接觸。在接觸過程中,靜摩擦因數(shù)設(shè)為0.3,動(dòng)摩擦因數(shù)設(shè)為0.1[19]。軌道采用完全約束,沖擊過程中的重力設(shè)置為9.8 N/kg。為了確保所建立的有限元模型能準(zhǔn)確計(jì)算多胞結(jié)構(gòu)在沖擊載荷下的行為,采用動(dòng)態(tài)沖擊試驗(yàn)對(duì)有限元模型進(jìn)行驗(yàn)證。試驗(yàn)在中南大學(xué)高速列車碰撞測試中心進(jìn)行,通過4 個(gè)20 號(hào)螺栓將吸能結(jié)構(gòu)固定在質(zhì)量為2 t 的試驗(yàn)臺(tái)車上,通過空氣炮賦予臺(tái)車速度,在列車撞擊位置放置1臺(tái)激光測速儀記錄結(jié)構(gòu)壓縮時(shí)的速度,為15.3 m/s。吸能結(jié)構(gòu)的壓縮過程采用高速攝影儀以3 000 幀/s 的幀率拍攝,最終得到開孔多胞結(jié)構(gòu)的力-時(shí)間對(duì)比曲線如圖5所示。從圖5 可見:多胞結(jié)構(gòu)試驗(yàn)得到的力-時(shí)間曲線和仿真的力-時(shí)間曲線整體上較吻合,在整個(gè)試驗(yàn)過程中,6 個(gè)波峰位置的力和時(shí)間基本一致,其中,仿真結(jié)果的谷值要稍高于試驗(yàn)結(jié)果的谷值,這是因?yàn)樵嚰苤谱鞴に嚨南拗?,?nèi)部有部分裂痕,這使得2條曲線存在一定差異。試驗(yàn)和仿真中峰值力、持續(xù)時(shí)間和吸能量的對(duì)比見表1。從表1可見結(jié)構(gòu)的峰值力、持續(xù)時(shí)間和吸能量這3個(gè)關(guān)鍵參數(shù)相對(duì)誤差分別為-2.67%,-3.02%和0.21%,這些差異對(duì)于后續(xù)的分析影響較小,因此,可采用有限元模型計(jì)算開孔多胞結(jié)構(gòu)的耐撞性。

    表1 關(guān)鍵參數(shù)的試驗(yàn)結(jié)果和仿真結(jié)果對(duì)比Table 1 Comparion of verification of optimal solution of key parameters

    圖5 力與時(shí)間的關(guān)系Fig.5 Relationship between crash force and time

    1.3 耐撞性評(píng)價(jià)指標(biāo)

    采用EA,F(xiàn)MC,F(xiàn)PC,ESA和ECL(EA為吸能量,F(xiàn)MC為平均壓潰力,F(xiàn)PC為最大壓潰力,ESA為結(jié)構(gòu)單位質(zhì)量所吸收的能量,ECL為結(jié)構(gòu)破碎過程中的載荷波動(dòng))這5 個(gè)耐撞性指標(biāo)對(duì)薄壁結(jié)構(gòu)的耐撞性進(jìn)行分析,所使用的計(jì)算方法見文獻(xiàn)[4]。其中,EA,F(xiàn)MC和ESA越高,說明結(jié)構(gòu)的吸能效果越好;ECL越接近于1,則結(jié)構(gòu)設(shè)計(jì)越合理;FPC越低,則安全性越好。

    在一般情況下,EA通過沖擊力-位移曲線中積分得到,但由于不同構(gòu)型的結(jié)構(gòu)在壓實(shí)時(shí)位移不一致,選擇相同的壓縮值來計(jì)算結(jié)構(gòu)的吸能量不準(zhǔn)確。本文引入吸能效率f[20]來計(jì)算結(jié)構(gòu)最大的位移,并將此時(shí)壓縮量對(duì)應(yīng)的吸能量作為結(jié)構(gòu)最終吸能量,計(jì)算式為:

    其中:s為結(jié)構(gòu)壓縮位移;Fmax為結(jié)構(gòu)折疊過程中所產(chǎn)生的除初始峰值力外的最大力。結(jié)構(gòu)在折疊過程中,EA和位移近似呈線性增加,此時(shí)f也隨著位移增加而增加。當(dāng)壓縮到密實(shí)階段時(shí),壓縮力急劇增加,而位移變化很小,此時(shí),f出現(xiàn)極大值,隨后f開始下降。因此,當(dāng)f達(dá)到最大值時(shí),認(rèn)為此時(shí)結(jié)構(gòu)到達(dá)壓實(shí)階段,吸能結(jié)構(gòu)的吸能量達(dá)到最大,可得到結(jié)構(gòu)的最大位移Dmax。在1次壓縮吸能過程中,結(jié)構(gòu)的EA,f和Dmax之間的關(guān)系如圖6所示。

    圖6 同一折疊過程中力、吸能量和吸能效率與位移的關(guān)系Fig.6 Relationships among force,EA,f and displacement in the same folding process

    2 開孔結(jié)構(gòu)變形模式分析

    2.1 4孔結(jié)構(gòu)變形模式分析

    通過有限元對(duì)誘導(dǎo)結(jié)構(gòu)進(jìn)行數(shù)值計(jì)算,其中有4個(gè)誘導(dǎo)孔的3種結(jié)構(gòu)在壓實(shí)時(shí)的變形結(jié)果如圖7 所示。從圖7 可見:在壓實(shí)時(shí)刻,這3 種結(jié)構(gòu)的內(nèi)側(cè)和外側(cè)的胞壁分別產(chǎn)生了7個(gè)和5個(gè)完整的褶皺,此外,S4-1 和S4-2 的內(nèi)側(cè)胞壁比S4-3 多半個(gè)褶皺。在最外側(cè)的胞壁中,S4-2 所形成的變形模式最穩(wěn)定,每個(gè)褶皺的折疊波長相近,而S4-1 和S4-3 的折疊波長出現(xiàn)了先增加后減小的情況,從而導(dǎo)致與最外側(cè)胞壁相鄰的結(jié)構(gòu)壁也受到了影響。對(duì)于內(nèi)側(cè)胞壁的變形模式,S4-2 和S4-3 以同時(shí)折疊的方式吸能,而S4-1 以交叉折疊的方式吸能。此外,S4-1 和S4-2 內(nèi)側(cè)的變形模式與S4-3 的變形模式不同,在第1次折疊時(shí),S4-1和S4-2的折疊方向與S4-3 內(nèi)側(cè)的折疊方向相反,從而使其內(nèi)側(cè)產(chǎn)生的褶皺比S4-3產(chǎn)生的褶皺多半個(gè)。

    圖7 4個(gè)誘導(dǎo)孔結(jié)構(gòu)的變形結(jié)果Fig.7 Deformation results of four windowed structures

    2.2 8孔結(jié)構(gòu)變形模式分析

    8 個(gè)誘導(dǎo)孔的3 種結(jié)構(gòu)在壓實(shí)時(shí)的變形結(jié)果如圖8所示。從圖8可見:S8-1在壓實(shí)時(shí),其外側(cè)胞壁產(chǎn)生6 個(gè)完整褶皺,內(nèi)側(cè)胞壁產(chǎn)生8 個(gè)完整褶皺;S8-2外側(cè)和內(nèi)側(cè)胞壁的褶皺數(shù)量分別為5個(gè)和7個(gè),而S8-3外側(cè)的褶皺數(shù)量為6個(gè),內(nèi)側(cè)胞壁的變形受到了誘導(dǎo)孔的影響,出現(xiàn)了5個(gè)和6個(gè)不同數(shù)量的褶皺;S8-1的變形較穩(wěn)定,變形模式可控,同一胞壁產(chǎn)生的褶皺波長較一致;S8-2 中,最外側(cè)胞壁的褶皺數(shù)量較少,內(nèi)側(cè)胞壁變形穩(wěn)定,產(chǎn)生的褶皺數(shù)量較多;而S8-3 中,雖然最外側(cè)胞壁褶皺較多,但內(nèi)側(cè)胞壁褶皺數(shù)量在所有結(jié)構(gòu)中最少,其內(nèi)側(cè)胞壁的變形模式穩(wěn)定性最低。

    圖8 8個(gè)誘導(dǎo)孔結(jié)構(gòu)的變形結(jié)果Fig.8 Deformation results of eight windowed structures

    2.3 開孔位置對(duì)于結(jié)構(gòu)變形的影響

    將在八邊形胞元上進(jìn)行開孔的結(jié)構(gòu)定義為內(nèi)側(cè)開孔結(jié)構(gòu),在六邊形胞元上進(jìn)行開孔的結(jié)構(gòu)定義為外側(cè)開孔結(jié)構(gòu)。不同位置下開孔對(duì)結(jié)構(gòu)褶皺數(shù)的影響見表2。從表2 可以看到:在結(jié)構(gòu)內(nèi)側(cè)的胞壁上進(jìn)行開孔,能夠增加結(jié)構(gòu)內(nèi)側(cè)胞壁的褶皺數(shù)量,而在內(nèi)側(cè)的每一條邊上都添加開孔時(shí),對(duì)多胞結(jié)構(gòu)最外側(cè)的褶皺數(shù)量產(chǎn)生影響,使外側(cè)胞壁的變形模式向內(nèi)胞壁的變形模式轉(zhuǎn)變,整體提高結(jié)構(gòu)的折疊次數(shù);與在結(jié)構(gòu)內(nèi)側(cè)增加誘導(dǎo)孔相比,單獨(dú)在最外側(cè)增加誘導(dǎo)孔對(duì)外側(cè)和內(nèi)側(cè)胞壁的變形模式影響較小。同樣地,在與內(nèi)側(cè)胞壁連接的外側(cè)胞壁上增加誘導(dǎo)孔對(duì)于結(jié)構(gòu)變形模式影響不大。作為連接外側(cè)和內(nèi)測的連接壁,其兩側(cè)胞壁的變形模式?jīng)]有受到影響,而在次外測的胞壁上加入誘導(dǎo)孔,雖然提高了最外側(cè)胞壁的褶皺數(shù)量,但同樣將結(jié)構(gòu)外側(cè)不規(guī)則變形的特性傳遞到了內(nèi)側(cè)胞壁上,使內(nèi)側(cè)胞壁中產(chǎn)生的褶皺數(shù)量降低。通過結(jié)構(gòu)內(nèi)側(cè)開孔方案和外側(cè)開孔方案,結(jié)構(gòu)的變形模式呈現(xiàn)出內(nèi)側(cè)胞壁變形較有序折疊,每一次折疊波長相近,而最外側(cè)胞壁的折疊波長不規(guī)律,每一次折疊波長都不一致,這是導(dǎo)致外側(cè)胞壁的褶皺數(shù)量少于內(nèi)側(cè)胞壁褶皺數(shù)量的主要原因。受結(jié)構(gòu)自身變形的影響,誘導(dǎo)孔會(huì)增大相鄰胞壁變形模式對(duì)其自身變形的影響,導(dǎo)致外側(cè)的不規(guī)則變形影響到內(nèi)側(cè)胞壁,因此,在內(nèi)側(cè)加入誘導(dǎo)孔能夠更好地提高結(jié)構(gòu)整體的褶皺數(shù)量,有利于結(jié)構(gòu)產(chǎn)生有序的折疊變形。

    表2 不同開孔位置對(duì)于胞壁褶皺數(shù)量的影響Table 2 Effects of different window positions on number of cell folds

    3 開孔結(jié)構(gòu)耐撞性分析

    將結(jié)構(gòu)根據(jù)不同的誘導(dǎo)孔數(shù)量進(jìn)行對(duì)比,耐撞性指標(biāo)如圖9 所示。從圖9(a),(b)和(c)可以看到:誘導(dǎo)孔的位置在由結(jié)構(gòu)內(nèi)側(cè)向外側(cè)變化時(shí),多胞結(jié)構(gòu)的FPC,EA和ESA都出現(xiàn)下降的情況,說明在多胞結(jié)構(gòu)外側(cè)增設(shè)誘導(dǎo)孔,能夠取得較好的降低FPC的效果,但對(duì)于提高結(jié)構(gòu)的EA和ESA會(huì)起到負(fù)面效果。從圖9(d)可以看到4 個(gè)誘導(dǎo)孔和8 個(gè)誘導(dǎo)孔結(jié)構(gòu)的變化趨勢不一致,這是由于S8-2 和S8-3 在吸能過程中,變形出現(xiàn)了一定的整體屈曲現(xiàn)象,造成在一定時(shí)間內(nèi)力均值上升,影響了吸能結(jié)構(gòu)的變形模式。從結(jié)構(gòu)的EA可以看出,盡管結(jié)構(gòu)S8-2 和S8-3 使FMC提高,但對(duì)于結(jié)構(gòu)的總吸能量是不利的。從圖9(e)可見,相對(duì)于內(nèi)側(cè)胞壁,在結(jié)構(gòu)外側(cè)胞壁加入誘導(dǎo)孔更利于提高結(jié)構(gòu)的ECL。

    圖9 不同開孔方案的耐撞性指標(biāo)Fig.9 Crashworthiness indexes of different window schemes

    4 開孔結(jié)構(gòu)選擇

    表3所示為上述6個(gè)開孔結(jié)構(gòu)的有限元仿真結(jié)果,可以看到誘導(dǎo)結(jié)構(gòu)無法同時(shí)在多個(gè)耐撞性指標(biāo)中達(dá)到最優(yōu)。為了更進(jìn)一步評(píng)估誘導(dǎo)結(jié)構(gòu)的耐撞性并找出最佳的誘導(dǎo)孔分布方案,使用復(fù)雜比例評(píng)估方法(complex proportional assessment,COPRAS)通過FPC,ECL和ESA對(duì)開孔結(jié)構(gòu)的性能進(jìn)行排序,選出較好的結(jié)構(gòu)。

    表3 開孔結(jié)構(gòu)的有限元仿真結(jié)果Table 3 Finite element results of windowed structures

    COPRAS 是一種用來解決對(duì)多準(zhǔn)則問題的決策方法,能夠從多個(gè)相沖突的目標(biāo)中,選擇出合適的方案[8]。該方法的步驟如下。

    步驟1:使用備選方案的標(biāo)準(zhǔn)值生成決策矩陣并進(jìn)行量綱一化。

    其中:xij為第j個(gè)評(píng)價(jià)指標(biāo)對(duì)于第i個(gè)備選方案的性能值;i≤m,j≤n,m和n分別為備選方案和評(píng)價(jià)指標(biāo)的數(shù)量。

    步驟2:確定每一個(gè)評(píng)價(jià)指標(biāo)的權(quán)重。

    對(duì)于要評(píng)估的對(duì)象,不同標(biāo)準(zhǔn)的重要性可能不同,采用權(quán)重表征不同標(biāo)準(zhǔn)的重要性。COPRAS中,權(quán)重的計(jì)算過程如下。

    1)將不同評(píng)價(jià)指標(biāo)之間進(jìn)行兩兩比較,總計(jì)進(jìn)行N=n(n-1)/2次比較。這里選擇其中比較重要的FPC,ESA和ECL作為評(píng)價(jià)指標(biāo)。

    2)對(duì)所有指標(biāo)的重要性進(jìn)行評(píng)分,不同指標(biāo)在進(jìn)行比較時(shí),指標(biāo)的重要性會(huì)出現(xiàn)相等或不相等共2 種情況。當(dāng)2 個(gè)指標(biāo)重要性相等時(shí),2 個(gè)指標(biāo)得分都為2;當(dāng)2 個(gè)指標(biāo)重要性不相等時(shí),其中更重要的指標(biāo)得分為3,另一個(gè)指標(biāo)得分為1。在結(jié)構(gòu)耐撞性評(píng)估中,一般認(rèn)為ESA和FPC是比較重要的參數(shù),兩者同等重要,而ECL則作為次要考慮目標(biāo)。

    3)通過下式得到不同指標(biāo)的權(quán)重ωj;

    其中:Sij為每一次2個(gè)指標(biāo)比較的評(píng)分;G為評(píng)價(jià)指標(biāo)的總分。

    步驟3:計(jì)算加權(quán)歸一化決策矩陣D。

    其中:dij為第i個(gè)備選方案的第j個(gè)標(biāo)準(zhǔn)的加權(quán)歸一化值。

    步驟4:對(duì)所有備選方案的有益屬性和無益屬性指標(biāo)的加權(quán)歸一化值求和。

    其中:S-為無益屬性指標(biāo);S+為有益屬性指標(biāo);S-i為第i個(gè)無益屬性指標(biāo);S+i為第i個(gè)有益屬性指標(biāo)。

    在本文的評(píng)價(jià)指標(biāo)中,F(xiàn)PC是無益屬性指標(biāo)。

    步驟5:計(jì)算相對(duì)優(yōu)先級(jí)結(jié)果和性能的定量效用Ui。

    其中:S-min為無益屬性指標(biāo)中的最小值。

    最終得到的備選方案中,其Ui越高則代表該備選方案的排名越高,最終得到的結(jié)果排名見表4。從表4 可以看到:綜合排名中,無孔初始結(jié)構(gòu)的評(píng)分最低。在3個(gè)指標(biāo)的權(quán)重中,S8-3由于綜合考慮了多個(gè)耐撞性指標(biāo),其排名最高,使結(jié)構(gòu)的FPC降低。此外,結(jié)構(gòu)的ECL提高,而ESA降低,這是由于使用COPRAS 方法中的權(quán)重分配原則,對(duì)耐撞性指標(biāo)的選擇有了偏向性,最終選擇的結(jié)果是一個(gè)較均衡的方案,因此,會(huì)忽略部分指標(biāo)的性能。該評(píng)分結(jié)果可供誘導(dǎo)孔結(jié)構(gòu)設(shè)計(jì)時(shí)參考,在實(shí)際應(yīng)用中,需要根據(jù)具體的設(shè)計(jì)要求進(jìn)行選擇。COPRAS 得到的最優(yōu)結(jié)構(gòu)和初始結(jié)構(gòu)耐撞性對(duì)比結(jié)果如表5 所示,可見結(jié)構(gòu)的FPC降低了14.20%,ECL提高了20.63%。

    表4 開孔結(jié)構(gòu)性能排名Table 4 Performance ranking of windowed structures

    表5 最優(yōu)開孔結(jié)構(gòu)和無孔結(jié)構(gòu)耐撞性對(duì)比Table 5 Comparison of crashworthiness between optimal windowed structure and initial structure

    5 結(jié)論

    1)建立了沖擊載荷下多胞結(jié)構(gòu)的有限元模型,并通過試驗(yàn)驗(yàn)證了模型的有效性。

    2)在結(jié)構(gòu)內(nèi)側(cè)胞壁加入誘導(dǎo)孔能夠提高結(jié)構(gòu)內(nèi)側(cè)胞壁的褶皺數(shù)量,同時(shí),對(duì)于結(jié)構(gòu)外側(cè)胞壁的褶皺數(shù)量有一定的影響。而在結(jié)構(gòu)外側(cè)胞壁加入誘導(dǎo)孔,能夠提升結(jié)構(gòu)外側(cè)胞壁的褶皺數(shù)量,但可能減少結(jié)構(gòu)內(nèi)側(cè)胞壁的褶皺數(shù)量。

    3)在結(jié)構(gòu)中由內(nèi)往外布置誘導(dǎo)孔,在多胞結(jié)構(gòu)內(nèi)側(cè)加入誘導(dǎo)孔與在結(jié)構(gòu)外側(cè)加入誘導(dǎo)孔相比,結(jié)構(gòu)的EA和ESA更高,而在結(jié)構(gòu)外側(cè)加入誘導(dǎo)孔,F(xiàn)PC更低。結(jié)構(gòu)內(nèi)側(cè)和外側(cè)誘導(dǎo)孔的布置側(cè)重點(diǎn)不同,兩者對(duì)于耐撞性指標(biāo)的影響是相互沖突的,如何平衡好兩者之間的關(guān)系,得到一個(gè)優(yōu)質(zhì)的多胞結(jié)構(gòu),有待進(jìn)一步研究。

    4)使用COPRAS 選擇出S8-3 作為綜合性能最均衡的開孔方案,結(jié)構(gòu)的FPC降低了14.20%,ECL提高了20.63%。

    猜你喜歡
    褶皺數(shù)量有限元
    動(dòng)漫人物衣服褶皺的畫法(1)
    統(tǒng)一數(shù)量再比較
    一點(diǎn)褶皺
    Coco薇(2017年6期)2017-06-24 23:08:56
    頭發(fā)的數(shù)量
    復(fù)雜褶皺構(gòu)造分析與找煤
    夢幻褶皺
    Coco薇(2015年7期)2015-08-13 22:33:38
    我國博物館數(shù)量達(dá)4510家
    磨削淬硬殘余應(yīng)力的有限元分析
    基于SolidWorks的吸嘴支撐臂有限元分析
    箱形孔軋制的有限元模擬
    上海金屬(2013年4期)2013-12-20 07:57:18
    亚洲欧美日韩高清在线视频| 999精品在线视频| 嫁个100分男人电影在线观看| 老熟妇仑乱视频hdxx| 国产精品亚洲美女久久久| 成人18禁高潮啪啪吃奶动态图| 热99re8久久精品国产| 国产亚洲av嫩草精品影院| tocl精华| 婷婷丁香在线五月| 久久国产亚洲av麻豆专区| 国产视频一区二区在线看| 18禁观看日本| 久久香蕉精品热| 午夜日韩欧美国产| 一区福利在线观看| 欧美av亚洲av综合av国产av| 精品一区二区三区视频在线观看免费| 看黄色毛片网站| 国产精品亚洲美女久久久| 成人欧美大片| 午夜日韩欧美国产| 九色国产91popny在线| 久久这里只有精品19| 麻豆国产av国片精品| 丁香欧美五月| 久久人人精品亚洲av| 亚洲国产欧美日韩在线播放| 91老司机精品| 两个人看的免费小视频| 欧洲精品卡2卡3卡4卡5卡区| 中亚洲国语对白在线视频| 88av欧美| 丝袜在线中文字幕| 成年版毛片免费区| 9191精品国产免费久久| 啦啦啦韩国在线观看视频| 可以在线观看毛片的网站| 欧美精品啪啪一区二区三区| 亚洲全国av大片| 男女下面进入的视频免费午夜 | 精品欧美国产一区二区三| 欧美黄色淫秽网站| 亚洲欧美日韩高清在线视频| 亚洲色图av天堂| 午夜精品国产一区二区电影| 国产精品一区二区三区四区久久 | 中文字幕最新亚洲高清| 女性被躁到高潮视频| 在线观看66精品国产| 一进一出抽搐gif免费好疼| 久久天堂一区二区三区四区| 亚洲狠狠婷婷综合久久图片| bbb黄色大片| 色尼玛亚洲综合影院| 在线观看舔阴道视频| 侵犯人妻中文字幕一二三四区| 日本五十路高清| 久久天堂一区二区三区四区| 母亲3免费完整高清在线观看| 久久这里只有精品19| 91老司机精品| 欧美成人午夜精品| 搡老岳熟女国产| 在线天堂中文资源库| 高清黄色对白视频在线免费看| 日本免费a在线| 精品欧美国产一区二区三| 国产激情久久老熟女| 精品国产国语对白av| 亚洲欧美激情在线| 午夜福利欧美成人| 老汉色av国产亚洲站长工具| 久久精品成人免费网站| 午夜久久久在线观看| 热re99久久国产66热| 欧美成人午夜精品| av视频免费观看在线观看| 黄色视频,在线免费观看| 高清毛片免费观看视频网站| 久久伊人香网站| 国产亚洲精品综合一区在线观看 | 18禁黄网站禁片午夜丰满| 亚洲avbb在线观看| 一区二区三区国产精品乱码| 在线观看免费视频日本深夜| 午夜福利一区二区在线看| 亚洲性夜色夜夜综合| 日韩大码丰满熟妇| 中文字幕色久视频| 一卡2卡三卡四卡精品乱码亚洲| 一区二区日韩欧美中文字幕| 亚洲欧美激情在线| 国产精品99久久99久久久不卡| 久久香蕉精品热| 亚洲少妇的诱惑av| 一级a爱片免费观看的视频| 国产精品影院久久| 精品人妻在线不人妻| 国产av一区二区精品久久| 国产精品久久久av美女十八| 欧美在线黄色| 黄色 视频免费看| 手机成人av网站| 精品高清国产在线一区| 国产激情欧美一区二区| 国产精品野战在线观看| 日本三级黄在线观看| 欧美日韩亚洲国产一区二区在线观看| 老司机午夜十八禁免费视频| 久久久国产成人免费| 一区在线观看完整版| 又大又爽又粗| www.www免费av| 黄色片一级片一级黄色片| 国产精品永久免费网站| 黄色成人免费大全| 日本欧美视频一区| 午夜福利免费观看在线| tocl精华| 91老司机精品| 午夜视频精品福利| 一二三四社区在线视频社区8| a级毛片在线看网站| 99久久综合精品五月天人人| 欧美精品啪啪一区二区三区| 国产又爽黄色视频| 欧美成人性av电影在线观看| 两性夫妻黄色片| 亚洲国产中文字幕在线视频| 国产免费男女视频| 好男人电影高清在线观看| 国产精品免费一区二区三区在线| 亚洲国产精品成人综合色| 免费不卡黄色视频| 中文字幕久久专区| 精品一区二区三区视频在线观看免费| 亚洲视频免费观看视频| 给我免费播放毛片高清在线观看| 欧美久久黑人一区二区| 国产午夜福利久久久久久| 男女做爰动态图高潮gif福利片 | 日韩欧美一区二区三区在线观看| 亚洲国产精品成人综合色| 亚洲第一欧美日韩一区二区三区| 9191精品国产免费久久| 天天躁夜夜躁狠狠躁躁| 亚洲精品美女久久av网站| 18禁裸乳无遮挡免费网站照片 | 中文字幕另类日韩欧美亚洲嫩草| 制服人妻中文乱码| 男女午夜视频在线观看| 后天国语完整版免费观看| 欧美日本视频| 他把我摸到了高潮在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 又黄又爽又免费观看的视频| 亚洲av日韩精品久久久久久密| 成人国语在线视频| 很黄的视频免费| 无限看片的www在线观看| 国产主播在线观看一区二区| 嫁个100分男人电影在线观看| 天堂影院成人在线观看| 午夜精品国产一区二区电影| 女性生殖器流出的白浆| 老司机午夜十八禁免费视频| 曰老女人黄片| 在线观看免费视频网站a站| 久久久久久久久久久久大奶| 国产日韩一区二区三区精品不卡| 亚洲三区欧美一区| 美女高潮到喷水免费观看| 欧美日韩一级在线毛片| 精品久久久久久久人妻蜜臀av | 欧美中文日本在线观看视频| 精品卡一卡二卡四卡免费| 黑人巨大精品欧美一区二区蜜桃| 免费看美女性在线毛片视频| 老司机福利观看| 精品福利观看| 亚洲,欧美精品.| 日韩大码丰满熟妇| 一级毛片高清免费大全| 国产精品影院久久| 啪啪无遮挡十八禁网站| 高潮久久久久久久久久久不卡| 中文亚洲av片在线观看爽| 最好的美女福利视频网| 欧美绝顶高潮抽搐喷水| 国产亚洲欧美精品永久| 成人亚洲精品av一区二区| 亚洲自偷自拍图片 自拍| 麻豆一二三区av精品| 久久久久久久久免费视频了| 久久青草综合色| 精品第一国产精品| 国产成人啪精品午夜网站| 午夜成年电影在线免费观看| 亚洲精品av麻豆狂野| 国产乱人伦免费视频| 亚洲午夜精品一区,二区,三区| 国产一级毛片七仙女欲春2 | 亚洲国产精品久久男人天堂| 中文字幕av电影在线播放| 日韩国内少妇激情av| 老司机靠b影院| 久久国产精品男人的天堂亚洲| 长腿黑丝高跟| 免费人成视频x8x8入口观看| 操出白浆在线播放| 91麻豆av在线| 啦啦啦韩国在线观看视频| 老司机在亚洲福利影院| 女人高潮潮喷娇喘18禁视频| 亚洲第一青青草原| 两个人视频免费观看高清| 亚洲色图 男人天堂 中文字幕| 99久久综合精品五月天人人| 麻豆一二三区av精品| 亚洲电影在线观看av| 美女免费视频网站| 欧美黑人欧美精品刺激| 亚洲性夜色夜夜综合| 久久九九热精品免费| 精品久久久久久成人av| 99国产精品99久久久久| 夜夜看夜夜爽夜夜摸| 日韩 欧美 亚洲 中文字幕| 90打野战视频偷拍视频| 国产精品爽爽va在线观看网站 | 日本免费a在线| 99香蕉大伊视频| 精品久久久久久久久久免费视频| 国产成+人综合+亚洲专区| 精品国产超薄肉色丝袜足j| 天天躁夜夜躁狠狠躁躁| 一二三四在线观看免费中文在| 97碰自拍视频| 亚洲精品美女久久久久99蜜臀| av网站免费在线观看视频| 色精品久久人妻99蜜桃| 国产av又大| 黑人欧美特级aaaaaa片| 久久久久久久久中文| 国产av一区在线观看免费| 午夜福利影视在线免费观看| 精品电影一区二区在线| 亚洲av美国av| 18禁美女被吸乳视频| www.www免费av| 99riav亚洲国产免费| 色av中文字幕| 亚洲国产精品合色在线| 国产又爽黄色视频| 亚洲精品中文字幕一二三四区| 亚洲人成电影免费在线| 人人妻,人人澡人人爽秒播| 亚洲av日韩精品久久久久久密| 欧美黄色淫秽网站| 麻豆成人av在线观看| 亚洲精品粉嫩美女一区| 别揉我奶头~嗯~啊~动态视频| 精品卡一卡二卡四卡免费| 国产精品久久久av美女十八| 最近最新中文字幕大全免费视频| 欧美日韩亚洲国产一区二区在线观看| 美女国产高潮福利片在线看| 久久精品国产亚洲av香蕉五月| 这个男人来自地球电影免费观看| 级片在线观看| 亚洲 欧美 日韩 在线 免费| 国产精品九九99| 操美女的视频在线观看| 99国产精品一区二区三区| 男人舔女人下体高潮全视频| 免费观看人在逋| 99久久精品国产亚洲精品| 超碰成人久久| 久久精品国产清高在天天线| 久久香蕉激情| 亚洲色图综合在线观看| 日本免费一区二区三区高清不卡 | 亚洲一区高清亚洲精品| 精品电影一区二区在线| 中文字幕人妻熟女乱码| 亚洲成人久久性| 色精品久久人妻99蜜桃| 欧美精品亚洲一区二区| 高清黄色对白视频在线免费看| 99国产精品99久久久久| 亚洲中文av在线| 亚洲情色 制服丝袜| 日韩精品免费视频一区二区三区| 精品一区二区三区四区五区乱码| 午夜福利高清视频| 在线观看免费午夜福利视频| 一进一出抽搐gif免费好疼| 久久久久九九精品影院| 久久国产精品影院| 母亲3免费完整高清在线观看| 成人永久免费在线观看视频| 久久久精品欧美日韩精品| 日日夜夜操网爽| 精品国产美女av久久久久小说| 精品欧美一区二区三区在线| 无遮挡黄片免费观看| 人人妻人人澡人人看| 激情视频va一区二区三区| 亚洲国产欧美一区二区综合| 精品福利观看| 中亚洲国语对白在线视频| 欧美精品啪啪一区二区三区| 国产伦人伦偷精品视频| 欧美性长视频在线观看| 久久中文字幕人妻熟女| 成人永久免费在线观看视频| 免费女性裸体啪啪无遮挡网站| 久久久久久久久中文| 变态另类丝袜制服| 国产精品野战在线观看| 成人免费观看视频高清| 色播在线永久视频| 国产精品爽爽va在线观看网站 | 熟妇人妻久久中文字幕3abv| 制服丝袜大香蕉在线| 亚洲中文日韩欧美视频| 国产一区二区三区在线臀色熟女| 真人做人爱边吃奶动态| 露出奶头的视频| 91精品国产国语对白视频| 国产欧美日韩精品亚洲av| 黑人欧美特级aaaaaa片| 国产一卡二卡三卡精品| 色尼玛亚洲综合影院| av在线播放免费不卡| 国产单亲对白刺激| 日韩 欧美 亚洲 中文字幕| www国产在线视频色| 婷婷精品国产亚洲av在线| 99在线视频只有这里精品首页| 在线观看免费视频日本深夜| 大香蕉久久成人网| 亚洲三区欧美一区| 国产av又大| 久久久久久免费高清国产稀缺| 电影成人av| 成年版毛片免费区| 一边摸一边抽搐一进一小说| 日本黄色视频三级网站网址| 午夜福利视频1000在线观看 | 午夜精品久久久久久毛片777| 一区二区三区激情视频| 欧美成人性av电影在线观看| 无限看片的www在线观看| 丝袜美腿诱惑在线| 亚洲最大成人中文| 亚洲aⅴ乱码一区二区在线播放 | 亚洲在线自拍视频| 亚洲,欧美精品.| 九色国产91popny在线| 久久中文看片网| 成人18禁高潮啪啪吃奶动态图| 国产蜜桃级精品一区二区三区| 国产激情久久老熟女| av欧美777| 国产成人啪精品午夜网站| 免费少妇av软件| 自拍欧美九色日韩亚洲蝌蚪91| 一进一出抽搐gif免费好疼| 色精品久久人妻99蜜桃| 国产aⅴ精品一区二区三区波| 亚洲免费av在线视频| 欧美精品啪啪一区二区三区| 亚洲久久久国产精品| 18禁裸乳无遮挡免费网站照片 | 亚洲精品久久国产高清桃花| 9热在线视频观看99| 欧美激情高清一区二区三区| 日韩大尺度精品在线看网址 | 性少妇av在线| 欧美亚洲日本最大视频资源| 日韩av在线大香蕉| 亚洲国产欧美网| 午夜免费观看网址| av片东京热男人的天堂| 一级黄色大片毛片| 十八禁网站免费在线| 可以免费在线观看a视频的电影网站| 九色亚洲精品在线播放| 国产精品精品国产色婷婷| 久久久久久久久免费视频了| 日本欧美视频一区| 国产精品美女特级片免费视频播放器 | 欧美另类亚洲清纯唯美| 精品电影一区二区在线| a在线观看视频网站| 国产精品亚洲美女久久久| 国产精品二区激情视频| 国产精品亚洲av一区麻豆| ponron亚洲| 亚洲人成电影观看| 色老头精品视频在线观看| 满18在线观看网站| 亚洲午夜精品一区,二区,三区| 男女下面插进去视频免费观看| 欧美性长视频在线观看| a级毛片在线看网站| 亚洲人成网站在线播放欧美日韩| 午夜免费激情av| 妹子高潮喷水视频| 国产成人一区二区三区免费视频网站| 男女午夜视频在线观看| 免费久久久久久久精品成人欧美视频| 国产高清激情床上av| av天堂久久9| 十八禁网站免费在线| 免费一级毛片在线播放高清视频 | 免费在线观看视频国产中文字幕亚洲| cao死你这个sao货| 女同久久另类99精品国产91| 一级作爱视频免费观看| 村上凉子中文字幕在线| 色尼玛亚洲综合影院| 午夜福利,免费看| 中文字幕av电影在线播放| 国产精品自产拍在线观看55亚洲| 少妇熟女aⅴ在线视频| 国产精品 国内视频| 一区二区三区激情视频| 三级毛片av免费| 女人高潮潮喷娇喘18禁视频| 国产精品久久电影中文字幕| 欧美老熟妇乱子伦牲交| 真人做人爱边吃奶动态| 久久精品国产亚洲av高清一级| 欧美日韩瑟瑟在线播放| 视频在线观看一区二区三区| 一进一出抽搐gif免费好疼| av片东京热男人的天堂| 91在线观看av| 一个人免费在线观看的高清视频| 不卡av一区二区三区| 母亲3免费完整高清在线观看| 亚洲成人免费电影在线观看| 无人区码免费观看不卡| 亚洲国产精品999在线| 多毛熟女@视频| 亚洲成av片中文字幕在线观看| 国产99久久九九免费精品| 日韩视频一区二区在线观看| 女性被躁到高潮视频| 国产成人影院久久av| 午夜视频精品福利| 中文字幕色久视频| 亚洲成人国产一区在线观看| 亚洲中文日韩欧美视频| 日韩av在线大香蕉| 色哟哟哟哟哟哟| 性少妇av在线| 男女做爰动态图高潮gif福利片 | 久久精品国产亚洲av香蕉五月| 9色porny在线观看| 黄色视频,在线免费观看| 一边摸一边抽搐一进一出视频| 久久久久九九精品影院| 波多野结衣高清无吗| 色综合欧美亚洲国产小说| 他把我摸到了高潮在线观看| 麻豆一二三区av精品| 亚洲欧美一区二区三区黑人| 在线观看66精品国产| 在线观看免费日韩欧美大片| 伊人久久大香线蕉亚洲五| 亚洲成a人片在线一区二区| 美女午夜性视频免费| 亚洲天堂国产精品一区在线| 国产三级黄色录像| 日韩国内少妇激情av| 亚洲黑人精品在线| 国产高清视频在线播放一区| 大码成人一级视频| 亚洲久久久国产精品| 欧美乱码精品一区二区三区| 久久精品成人免费网站| 嫩草影视91久久| 亚洲av电影不卡..在线观看| 天堂动漫精品| 搞女人的毛片| 又大又爽又粗| 久久久久亚洲av毛片大全| 亚洲av电影不卡..在线观看| 国产欧美日韩精品亚洲av| 国产精品亚洲美女久久久| 亚洲精品久久国产高清桃花| 亚洲av成人不卡在线观看播放网| 亚洲国产中文字幕在线视频| 亚洲av成人av| 欧美性长视频在线观看| 久久午夜亚洲精品久久| 亚洲国产精品sss在线观看| 久久影院123| 久久久精品欧美日韩精品| 欧美精品亚洲一区二区| 国产欧美日韩综合在线一区二区| 亚洲国产精品999在线| 两个人看的免费小视频| 丝袜美腿诱惑在线| 一区二区三区激情视频| 亚洲欧美激情综合另类| 久久天堂一区二区三区四区| netflix在线观看网站| 午夜福利在线观看吧| 久久人人97超碰香蕉20202| 欧美另类亚洲清纯唯美| 国产一区二区在线av高清观看| 一本大道久久a久久精品| 亚洲五月婷婷丁香| av中文乱码字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 天堂√8在线中文| 国产精品亚洲美女久久久| 一区二区三区精品91| 精品电影一区二区在线| 国产成人啪精品午夜网站| 国产一区二区在线av高清观看| 少妇熟女aⅴ在线视频| 一区在线观看完整版| 亚洲欧美日韩无卡精品| 午夜免费鲁丝| 国产成人精品无人区| 一区二区三区高清视频在线| 午夜福利成人在线免费观看| 国产三级黄色录像| 国产精品免费视频内射| 性少妇av在线| 久久天堂一区二区三区四区| 精品久久久精品久久久| 十八禁人妻一区二区| 久久久久亚洲av毛片大全| 级片在线观看| 欧美激情高清一区二区三区| 亚洲精品中文字幕在线视频| 老司机福利观看| 亚洲精品国产精品久久久不卡| 天天一区二区日本电影三级 | 国产精品影院久久| 久久午夜综合久久蜜桃| av在线天堂中文字幕| 精品国产亚洲在线| 精品午夜福利视频在线观看一区| 99在线人妻在线中文字幕| 亚洲国产精品合色在线| 操出白浆在线播放| 国产免费av片在线观看野外av| 在线国产一区二区在线| 亚洲人成电影观看| 国产成人系列免费观看| 可以在线观看毛片的网站| 国产成人精品久久二区二区91| 一进一出抽搐gif免费好疼| 法律面前人人平等表现在哪些方面| 久久久久久国产a免费观看| 一级a爱视频在线免费观看| 久久久精品国产亚洲av高清涩受| 国产欧美日韩一区二区三| 国产伦一二天堂av在线观看| 好男人在线观看高清免费视频 | 久久久久久久久久久久大奶| 波多野结衣高清无吗| 久久人妻av系列| 少妇裸体淫交视频免费看高清 | 日日摸夜夜添夜夜添小说| 人成视频在线观看免费观看| 神马国产精品三级电影在线观看 | 亚洲avbb在线观看| or卡值多少钱| 999久久久精品免费观看国产| 欧美成人免费av一区二区三区| 国产一区二区在线av高清观看| 九色亚洲精品在线播放| 欧美色欧美亚洲另类二区 | 在线国产一区二区在线| 在线观看免费日韩欧美大片| 欧美在线黄色| 亚洲av片天天在线观看| 欧美乱妇无乱码| 一区福利在线观看| 韩国av一区二区三区四区| 国产av又大| 精品人妻在线不人妻| 美国免费a级毛片| 欧美人与性动交α欧美精品济南到| 巨乳人妻的诱惑在线观看| 一级毛片高清免费大全| 男人的好看免费观看在线视频 | 好男人电影高清在线观看| 欧美大码av| 国产伦人伦偷精品视频| 亚洲专区中文字幕在线| 黄色丝袜av网址大全| 丝袜美腿诱惑在线| 久久久久久人人人人人| 在线播放国产精品三级| 制服丝袜大香蕉在线| 亚洲专区中文字幕在线| 精品第一国产精品| 亚洲黑人精品在线| 午夜影院日韩av| 日日摸夜夜添夜夜添小说| 亚洲人成电影观看| 国产不卡一卡二| 99久久99久久久精品蜜桃| 青草久久国产| 久热这里只有精品99|