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

    具有等距密排微小孔穿刺模板的力學(xué)分析

    2016-08-05 02:24:49趙世海王婷慧董九志蔣秀明天津工業(yè)大學(xué)機(jī)械工程學(xué)院天津300387天津工業(yè)大學(xué)現(xiàn)代機(jī)電裝備技術(shù)天津市重點(diǎn)實(shí)驗(yàn)室天津300387
    關(guān)鍵詞:單胞等距細(xì)觀

    趙世海,王婷慧,董九志,蔣秀明(1.天津工業(yè)大學(xué) 機(jī)械工程學(xué)院,天津 300387;2.天津工業(yè)大學(xué) 現(xiàn)代機(jī)電裝備技術(shù)天津市重點(diǎn)實(shí)驗(yàn)室,天津 300387)

    ?

    具有等距密排微小孔穿刺模板的力學(xué)分析

    趙世海1,2,王婷慧1,2,董九志1,2,蔣秀明1,2
    (1.天津工業(yè)大學(xué) 機(jī)械工程學(xué)院,天津 300387;2.天津工業(yè)大學(xué) 現(xiàn)代機(jī)電裝備技術(shù)天津市重點(diǎn)實(shí)驗(yàn)室,天津 300387)

    穿刺模板用于正交疊層碳布織機(jī)中,其上開(kāi)有近千個(gè)緊密排列的微孔使得力學(xué)分析較為困難.首先根據(jù)擬板理論用2種方法求出了等效彈性模量:一種是從細(xì)觀層次出發(fā)由漸進(jìn)均勻化理論得出;另一種是從宏觀角度應(yīng)用薄板理論,分別得出了剩余彈性模量和孔密度的關(guān)系圖.然后根據(jù)工程需要利用擬板求出了穿刺模板的等效板厚和自振頻率.最后設(shè)計(jì)了實(shí)驗(yàn)以測(cè)試?yán)碚摲治鏊玫年P(guān)系圖,綜合2種理論方法和實(shí)驗(yàn)所得的3個(gè)結(jié)果找出了在不同孔密度區(qū)間里求等效彈性模量的最優(yōu)方法.

    穿刺模板;等距密排微小孔;均勻化理論;等效彈性模量;擬板理論;力學(xué)分析

    modulus;pseudo-plates theory;mechanical analysis

    等距密排穿刺模板在碳布的整體穿刺工藝中用于制造碳纖維立體織物[1].碳纖維立體織物是高性能碳纖維在三維空間多方向連續(xù)排布形成的纖維增強(qiáng)體.整體穿刺的過(guò)程是將鋼針?lè)胖贸蒢向鋼針陣列,固定于工作臺(tái)上,碳布置于穿刺模板下方,在穿刺模板帶動(dòng)下穿過(guò)Z向鋼針陣列,重復(fù)上述過(guò)程直到滿足所需厚度,使用Z向纖維置換Z向鋼針,制成立體織物.

    在整體穿刺過(guò)程中穿刺模板向下移動(dòng),接觸靜止的碳布后,帶動(dòng)碳布一同穿過(guò)靜置的Z向鋼針陣列,當(dāng)碳布與鋼針陣列接觸時(shí)開(kāi)始穿刺過(guò)程,鋼針對(duì)碳布有作用力,導(dǎo)致碳布變形,與此同時(shí),碳布對(duì)穿刺模板的微小孔邊產(chǎn)生力的作用.

    在此過(guò)程中穿刺模板起到為鋼針定位和給碳布加壓的作用,隨著碳布層數(shù)的增加,碳纖維立體織物的厚度加大,鋼針在穿刺中承受較大的反作用力后會(huì)傾斜、彎曲甚至折斷,穿刺模板的微小孔約束了鋼針的上述運(yùn)動(dòng);另外,由于各層碳布是否密實(shí)直接影響立體織物的力學(xué)性能,這時(shí)穿刺模板起到了加壓的作用.因此,穿刺模板必須有足夠的剛度,以保證在整個(gè)穿刺過(guò)程中發(fā)生盡可能小的彎曲變形.由于開(kāi)孔對(duì)板的彎曲剛度影響很大,含有近千個(gè)微小孔的穿刺模板和加工板所用的材料相比,其彎曲剛度大為削弱.

    穿刺模板的外部尺寸符合薄板彎曲理論的條件,但是利用ANSYS直接建模分析時(shí)會(huì)出現(xiàn)2個(gè)障礙導(dǎo)致分析中止:一是由于孔很小,應(yīng)力和撓度以及其它變形分量等一些物理量隨位置變化呈現(xiàn)震蕩特征,使得力學(xué)分析過(guò)程復(fù)雜且結(jié)果不可靠;二是如果直接應(yīng)用有限元進(jìn)行數(shù)值模擬,劃分網(wǎng)格時(shí)由于每個(gè)微孔的四周都有孔,每個(gè)孔的圓周上至少要分4個(gè)節(jié)點(diǎn)。穿刺模板上孔的數(shù)目多,描述該類結(jié)構(gòu)需要很密的網(wǎng)格從而形成很多單元,穿刺模板上有近千個(gè)孔使得上述計(jì)算量巨大,超出微型計(jì)算機(jī)的計(jì)算能力,導(dǎo)致無(wú)法求解或者求解精度過(guò)低.

    穿刺模板的結(jié)構(gòu)特點(diǎn)是在板的中心區(qū)域按照一定規(guī)律開(kāi)了很多孔,孔徑比板的輪廓尺寸小得多.本文采用擬板理論,把穿刺模板的微小孔抹去,用同尺寸的實(shí)心板代替穿刺模板,由此可知代替穿刺模板完成力學(xué)分析的等效實(shí)心板與穿刺模板應(yīng)該有同樣的彈性常數(shù),以保證等效實(shí)心板和穿刺模板在受到相同的激勵(lì)載荷時(shí)響應(yīng)相同.

    1 穿刺模板的力學(xué)分析

    圖1表示了穿刺模板的細(xì)部結(jié)構(gòu).由圖1可以看出,穿刺模板的結(jié)構(gòu)特點(diǎn)是緊密排列的微孔中心距相等,排布顯周期性特征.

    圖1 穿刺模板結(jié)構(gòu)模型Fig.1 Structure model of piercing template

    本文先根據(jù)擬板理論把密排孔板的孔抹去使其等價(jià)為另一種材料制成的實(shí)心板,也就與是密排孔板等價(jià)的擬板,在不改變板的形狀和尺寸的條件下求出擬板的彈性模量定義為密排孔板的等效彈性模量[2].借鑒復(fù)合材料板的分析方法,用細(xì)觀和宏觀2種方法求出2組等效彈性模量[3-4].

    細(xì)觀方法應(yīng)用基于小參數(shù)攝動(dòng)法的漸進(jìn)均勻化理論和有限元分析的理論[5].先根據(jù)細(xì)觀力學(xué)中的代表單元法劃分出單胞,然后將單胞的應(yīng)力應(yīng)變以及應(yīng)變位移的關(guān)系式代入虛功方程,將位移展開(kāi)為含小參數(shù)的漸進(jìn)表達(dá)式后代入上述方程,利用周期性性質(zhì)和散度定理化簡(jiǎn)方程得出等效彈性模量的張量表達(dá)式;最后利用通用有限元分析的格式將張量式改寫(xiě)為矩陣方程,求出剛度矩陣從而得到等效彈性模量[6].

    根據(jù)結(jié)構(gòu)的周期性劃分出的代表單元稱為單胞,理想單胞可以通過(guò)空間上無(wú)重疊、無(wú)縫隙的延拓即可重現(xiàn)真實(shí)的密排孔板結(jié)構(gòu)并且要有規(guī)則的幾何外形.此外如果單胞外廓閉合,則在二維坐標(biāo)系中是單連通域,有利于物理方程的化簡(jiǎn)中直接應(yīng)用高斯散度定理.基于以上原因穿刺模板的單胞模型如圖2所示[7].

    圖2 單胞模型Fig.2 Model of base cell

    建立2個(gè)尺度的坐標(biāo)系,宏觀用X表示,細(xì)觀用Y表示.引入小量,其物理意義是把宏觀坐標(biāo)系X中的一個(gè)點(diǎn)放大倍放入細(xì)觀坐標(biāo)系Y中,建立了一個(gè)模型在宏觀和細(xì)觀上對(duì)應(yīng)關(guān)系.穿刺模板在垂直于板中面的方向孔的形狀和尺寸沿厚度不變,在彈性力學(xué)中可歸結(jié)為的平面問(wèn)題處理并且密排孔板的彈性矩陣設(shè)為正交各向異性材料的形式,可知待求量為E11、E12、E22和E66.

    經(jīng)數(shù)學(xué)推導(dǎo)得出等效彈性模量的表達(dá)式

    宏觀方法是基于經(jīng)典薄板理論,以密排孔板和擬板受相同的外載荷時(shí)的最大撓度相等為條件建立方程,從中提出抗彎剛度的關(guān)系式從而得到等效彈性模量[9].又由于等效前后板的外形尺寸不能變,所以微孔數(shù)目越多,板的實(shí)體部分的體積相應(yīng)就越小,每個(gè)小孔周邊承受的載荷就越大,這就建立起孔的數(shù)量和彈性模量的關(guān)系.綜上可得Timoshenko理論,其中正方形薄板的中心撓度的精確值為

    2 案例分析

    密排孔板的外形尺寸為20 cm×20 cm,等距密排孔區(qū)域面積為256 cm2,厚度為1 cm.有976個(gè)微孔,孔徑為4 mm,孔邊間距1 mm.密排孔板的彈性模量為210×106Pa,泊松比μ=0.3,密度為7.85×103kg/m3,表面受壓力載荷為4×104N/m2.用MAPLE解有限元矩陣的結(jié)果為:

    由觀察結(jié)果可知,在板上的微孔是圓形且等距排列時(shí),密排孔板的彈性矩陣顯各向同性材料的特征.

    另一種解法是由上文推出的邊界條件直接用ANSYS建單胞模型[10].為了保證單胞的變形協(xié)調(diào)還需再加上周期性邊界條件,即保證單胞外邊界相對(duì)的2個(gè)面的所有點(diǎn)的位移變化一致[11].具體實(shí)施的方法是在WORKBENCH中用“強(qiáng)迫位移”命令添加位移約束. 圖3為單胞水平方向有一個(gè)單位初應(yīng)變的變形,剩余2種情況類似.

    圖3 單胞變形圖Fig.3 Deformation of base cell

    ANSYS與MAPLE給出的結(jié)果相差無(wú)幾,原因是2種方法本質(zhì)是一樣的,只是同一個(gè)矩陣方程用2種計(jì)算器求解而已.用MAPLE的解是先借助ANSYS APDL劃分網(wǎng)格,導(dǎo)出此時(shí)節(jié)點(diǎn)和單元信息然后導(dǎo)入MAPLE,用ANSYS求解是建模、劃分網(wǎng)格、加載邊界條件全程借助ANSYS的模塊,優(yōu)點(diǎn)是可以借助ANSYS強(qiáng)大的后處理器得到直觀的變形圖像.

    從宏觀角度來(lái)看板的中心撓度就是最大撓度故取密排孔板中心區(qū)域的小樣[12].令密排孔板的中心撓度為w,彈性模量E,擬板的中心撓度為wH彈性模量EH.借助ANSYS APDL建立小樣模型,求出w和w′再代入(3)求出.用同樣的方法能求出不同開(kāi)孔直徑下的[13].

    由于最終要得出開(kāi)孔數(shù)量對(duì)板的彈性模量的削弱程度,故定義剩余彈性模量為EH/E,其中分母為加工密排孔板所用材料的彈性模量.開(kāi)孔數(shù)量用孔密度表示,定義為α等于孔板的體積和當(dāng)量板體積比值.2種方法求出的剩余彈性模量和孔密度的關(guān)系如圖4所示.

    圖4 剩余彈性模量-孔密度圖Fig.4 Residual modulus of elasticity-density of hole

    觀察圖4得出2條結(jié)論:一是當(dāng)板上的微孔是圓形、等距排列且在孔密度小于0.4時(shí),細(xì)觀方法和宏觀方法數(shù)值相近,可以推測(cè)出按此方法分析時(shí)最佳孔密度為0.3到0.4,但是隨著孔的數(shù)量增多,計(jì)算所得剩余彈性模量數(shù)值的離散度增加,到底哪種更接近實(shí)際由后續(xù)試驗(yàn)驗(yàn)證;另一個(gè)是可以根據(jù)曲線擬合原理得出α和的關(guān)系式接近指數(shù)函數(shù)形式

    式中:x、y和z為待定常數(shù).

    3 擬板的應(yīng)用

    在應(yīng)用擬板理論求出等效彈性模量的過(guò)程中要求擬實(shí)心板和密排孔板的外形尺寸一致,以保證兩者在承受相同的外載荷時(shí)把密排孔對(duì)剛度的削弱表現(xiàn)在彈性模量上.如果把原密排孔板和擬板都看由2種無(wú)關(guān)的材料制成的實(shí)心板,當(dāng)兩板的表面積一樣時(shí)就把密排孔對(duì)剛度的削弱表現(xiàn)為厚度的折減,求出擬板的厚度值稱為等效板厚.

    根據(jù)克?;舴?泰勒假設(shè),選取中面撓度為基本未知量,另外2個(gè)位移分量是中面撓度的函數(shù),在直角坐標(biāo)系中彎曲微分方程為

    式中:q為垂直于中面的外載.

    根據(jù)中面的定義,計(jì)算出中面位置h1

    此時(shí)的抗彎剛度Dh1計(jì)算公式變?yōu)?/p>

    下面推導(dǎo)計(jì)算等效板厚的公式,目的是在已知密排孔板的材料的情況下,快速判斷出相當(dāng)于多厚的實(shí)心板.要說(shuō)明的是等效板厚公式的前提:孔為圓形,正方形,正六邊形等規(guī)則幾何輪廓;每個(gè)孔的中心距一樣,顯正三角形或正方形排列,目的是保證能利用各向同性板的理論;厚度方向孔的尺寸不變化.如果存在階梯孔可以根據(jù)厚度方向的不同孔徑分段處理,由于以上所有推導(dǎo)都是在小撓度線彈性的前提下進(jìn)行,可以把分段處理的結(jié)果直接相加得到整個(gè)板的等效厚度.

    建立方程的依據(jù)是保證兩者在相同外載荷下的相同變形,把密排孔對(duì)剛度的削弱表現(xiàn)為厚度的折減.厚度H由彎曲剛度表征,所以直接令兩者的彎曲剛度相等即可求出等效厚度HH,H為密排孔板的實(shí)際厚度.

    借助ANSYS得到原厚度為10 mm的穿刺模板相當(dāng)于厚度為5.71 mm的實(shí)心板.

    穿刺時(shí)模板在隨工作臺(tái)上下往復(fù)運(yùn)動(dòng),這種工況下通常需求出穿固有頻率.借助ANSYS的模態(tài)組塊,可以求出任意階的自振頻率,用擬板建模從后處理器可提取的各階頻率的數(shù)值和圖形,得到的一階自振頻率即密排孔板的固有頻率.

    4 剩余彈性模量的驗(yàn)證試驗(yàn)

    在理論分析中宏觀和細(xì)觀2種方法的操作中必須用有限元方法離散,解的精確性和節(jié)點(diǎn)的數(shù)量、網(wǎng)格的形狀及單元的選取密切相關(guān),故在理論分析后設(shè)計(jì)了試驗(yàn)測(cè)試圖4得出剩余彈性模量和孔密度的關(guān)系.實(shí)驗(yàn)設(shè)備為天津工業(yè)大學(xué)機(jī)械學(xué)院材料力學(xué)實(shí)驗(yàn)室的WDW3100微機(jī)控制電子萬(wàn)能試驗(yàn)機(jī).

    至今為止沒(méi)有針對(duì)具有密排孔結(jié)構(gòu)的板的拉伸試驗(yàn)標(biāo)準(zhǔn),本次試驗(yàn)的試樣參照國(guó)家標(biāo)準(zhǔn)中短比例薄板試樣的規(guī)范,取試樣的總長(zhǎng)度為200mm,標(biāo)距為50 mm.根據(jù)理論分析的結(jié)果可知對(duì)于薄板可忽略厚度對(duì)剩余彈性模量影響,考慮到試樣的拉斷時(shí)間和加工方便,直接從2 mm厚的鋼板取樣.由于密排孔區(qū)域強(qiáng)度削弱較大,故試樣在拉伸過(guò)程中一定是在密孔區(qū)域而不會(huì)在夾持部分屈服,故沒(méi)有設(shè)計(jì)夾持過(guò)度段,斷裂后的試樣也驗(yàn)證了這點(diǎn).

    表1列出了本次試驗(yàn)的4種試樣.

    表1中,0號(hào)試樣是未開(kāi)孔板,目的是測(cè)出剩余彈性模量中分母的數(shù)值.1、2、3號(hào)試樣分別測(cè)出開(kāi)孔密度為0.2、0.3、0.6時(shí)的等效彈性模量的值.其中1號(hào)和3號(hào)試樣的孔為等距排布,2號(hào)試樣的孔為橫縱向中心距不相等的情況.受加工條件的制約2和3號(hào)試樣沒(méi)能做到孔密度嚴(yán)格一致,但是相近的孔密度也能反映出孔的排布方式對(duì)剩余彈性模量的影響程度,而且以橫縱向距離不相等的試樣做反向?qū)φ?,為進(jìn)一步分析不等距的密排孔板的提供實(shí)驗(yàn)數(shù)據(jù).同樣受加工條件的限制本次試驗(yàn)的試樣最大孔密度取到60%,表2為實(shí)驗(yàn)數(shù)據(jù).

    表1 試樣規(guī)格明細(xì)Tab.1  Specification list of sample

    表2 實(shí)驗(yàn)結(jié)果Tab.2 Experimental data

    每種試樣都取3件做試驗(yàn),表2的結(jié)果是3件試樣的算術(shù)平均值.屈服力是按照慣例取0.2%的殘余變形時(shí)的值.0號(hào)試樣測(cè)出了剩余彈性模量中分母的值,其余試樣的彈性模量都和其相比.1號(hào)試樣的剩余彈性模量的實(shí)驗(yàn)值為0.37,理論值為0.297,2號(hào)試樣的實(shí)驗(yàn)值為0.612,理論值為0.65,3號(hào)試樣的實(shí)驗(yàn)值為0.718,理論值為0.74.分析實(shí)驗(yàn)數(shù)據(jù)可以直接看出2點(diǎn):一是當(dāng)開(kāi)孔密度低于0.6時(shí),用宏觀薄板理論和細(xì)觀均勻化方法求出的理論值和實(shí)驗(yàn)值吻合的較好,但隨開(kāi)孔密度孔的增加數(shù)值的離散度加大,本次試驗(yàn)的數(shù)據(jù)更接近薄板理論的分析結(jié)果,可見(jiàn)均勻化方法的結(jié)果更趨于保守;二是像2號(hào)試樣這樣每排孔的中心距統(tǒng)一變化時(shí),孔間距對(duì)剩余彈性模量的影響不大.因?yàn)殡m然2號(hào)試樣把縱向的孔距增大,按照開(kāi)孔率的計(jì)算公式板上有孔處的總面積也隨之加大,降低了開(kāi)孔率,并且整排孔的變動(dòng)數(shù)值一樣時(shí)結(jié)果類似于等距的情況.

    觀察試樣的斷裂位置發(fā)現(xiàn)同一規(guī)格的三件試樣的斷裂位置不一樣,即并非都發(fā)生在中間一排,并且2號(hào)試樣在相鄰兩排孔都有斷裂,1號(hào)和3號(hào)試樣斷裂雖然不一定發(fā)生在哪一排,但都是同一排孔統(tǒng)一短.由此可以看出改變孔間距雖然對(duì)剩余彈性模量的影響不大,但對(duì)板的破壞形式有影響.同時(shí)可以推測(cè)出斷裂發(fā)生在孔間距的加工誤差較大的位置,那么隨孔密度加大,加工誤差影響成為導(dǎo)致實(shí)際值的離散度加大的重要因素.

    圖5為未開(kāi)孔試樣的實(shí)驗(yàn)力-位移曲線,圖6為孔密度為0.6的曲線,結(jié)合其他試樣的曲線圖可以發(fā)現(xiàn),隨著孔密度的增大密排孔結(jié)構(gòu)試樣的曲線沒(méi)有明顯的彈性階段、屈服階段、強(qiáng)化階段的分界區(qū)間.

    圖5 0號(hào)試樣的實(shí)驗(yàn)力-位移曲線Fig.5 Experimental force displacement curve of sample 0

    圖6 1號(hào)試樣的實(shí)驗(yàn)力-位移曲線Fig.6 Experimental force displacement curve of sample 1

    5 結(jié)語(yǔ)

    具有等距密排微小孔的穿刺模板的力學(xué)分析不能直接用有限元方法仿真,本文用擬板理論,用模擬出的實(shí)心板分析密排孔板的一些力學(xué)指標(biāo).為了保證結(jié)果的可靠性,用細(xì)觀均勻化方法和薄板理論求出了等效彈性模量并得出了幾種不同規(guī)格的開(kāi)孔數(shù)量對(duì)彈性模量弱度,即剩余彈性模量-孔密度曲線;隨后列舉了擬板的兩個(gè)用途;最后設(shè)計(jì)了試驗(yàn)測(cè)試?yán)碚摲治龅慕Y(jié)果,發(fā)現(xiàn)由細(xì)觀均勻化方法求出的數(shù)值均稍微低于實(shí)驗(yàn)所得,說(shuō)明此方法所得數(shù)據(jù)更可靠.

    需要指出的是,這種模擬的實(shí)心板不是在任何情況都能代替密排孔板,比如在強(qiáng)度分析中由于孔的存在會(huì)造成局部的應(yīng)力集中,在疲勞破壞分析中沖擊載荷作用時(shí)某些邊緣區(qū)域的孔會(huì)先出現(xiàn)破壞.同時(shí)實(shí)驗(yàn)表明加工精度對(duì)密排孔板力學(xué)性能影響很大,例如板的某些邊緣區(qū)域孔間距如果有很小的誤差,受彎曲外載荷時(shí)破壞失效很可能就發(fā)生于此,這些用模擬的實(shí)心板都難以得出.

    [1]朱建勛,何建敏,王海燕.正交疊層機(jī)布整體穿刺工藝的纖維彎曲伸長(zhǎng)機(jī)理[J].中國(guó)工程科學(xué),2003,5(5):59-62.

    ZHU J X,HE J M,WANG H Y.Fiber elongation mechanism orthogonal bending laminated cloth overall puncture process[J]. Engineering Sciences,2003,5(5):59-62(in Chinese).

    [2]謝桂蘭,趙錦梟,田杰,等.均勻化有限元法預(yù)測(cè)復(fù)合材料層合板宏觀有效彈性性能 [J].玻璃鋼/復(fù)合材料,2014(7):23-27(in Chinese).

    XIE G L,ZHAO J X,TIAN J,et al.Predicting macroscopic effective elastic properties of composite laminated plate by homogenization fem[J].Fiber Reinforced Plastics/Composites,2014(7):23-27(in Chinese).

    [3]杜修力,金瀏.細(xì)觀均勻化方法預(yù)測(cè)非飽和混凝土宏觀力學(xué)性質(zhì)[J].水利學(xué)報(bào),2013,47(3):355-371.

    DU X L,JIN L.Micro-scale homogenization for prediction of the mechanical properties of unsaturated concrete[J].Journal of Hydraulic Engineering,2013,47(3):355-371(in Chinese).

    [4]TERADA K,KATO J,HIRAYAMA H,et al.A method of two-scale analysis with micro-macro decoupling scheme:Application to hyperelastic composite materials[J].Comput Mech,2013,52:1199-1219.

    [5]HASSANI B,HINTON E.A Review of homogenization and topology optimization I-homogenization theory for media with periodic structures[J].Computers and Structures,1998(69):707-717.

    [6]DAVALOS-SOTELO Raymundo.Determination of elastic properties of clear wood by the homogenization method in two dimensions[J].Wood Sci Technol,2005,39:385-417.

    [7]董紀(jì)偉,孫良新,洪平.基于均勻化方法的三維編織復(fù)合材料等效彈性性能預(yù)測(cè)[J].宇航學(xué)報(bào),2003,26(4):482-486.

    DONG J W,SUN L G,HONG P.Homogenization-based method for predicting effective elastic properties of three-dimensional braided composites[J].Journal of Asronautics,2003,26(4):482-486(in Chinese).

    [8] HASSANI B,HINTON E.A review of homogenization and topology optimization II-analytical and numerical solution of homogenization equations[J].Computers and Structures,1998 (69):719-738.

    [9]伍莉,劉均,程遠(yuǎn)勝.方形蜂窩夾芯夾層板彎曲問(wèn)題的新解法[J].船舶力學(xué),2012(8):926-934.

    WU L,LIU J,CHENG Y S.A new algorithm for solving being response of square-honeycomb sandwich plates[J].CSSRC Reports,2012(8):926-934(in Chinese).

    [10]劉洋,何斌,傅潔,等.周期變截面梁剛度分析的均勻化方法[J].機(jī)械設(shè)計(jì)與制造,2014(1):21-24

    LIU Y,HE B,F(xiàn)U J,et al.Stiffness analysis of the periodic variablecross-sectionbeambasedonhomogenizationmethod[J]. Machinery Design&Manufacture,2014(1):21-24(in Chinese).

    [11]劉遠(yuǎn),金國(guó)光,劉又午,等.基于有限段方法的高速凸輪機(jī)構(gòu)動(dòng)力學(xué)分析[J].天津工業(yè)大學(xué)學(xué)報(bào),2007,26(6):66-69.

    LlU Y,JIN G G,LIU Y W,et al.Dynamic analysis of highspeed cam mechanisms based on finite-segment method[J]. Journal of Tianjin Polytechnic University,2007,26(6):66-69 (in Chinese).

    [12]王珂,劉彤,董其伍.復(fù)雜多孔板有效彈性常數(shù)的結(jié)構(gòu)分析方法[J].壓力容器,2010,27(5):10-14.

    WANG K,LIU T,DONG Q W.Strucual analysis methods to determine the effective elastic constants for perforated plates[J]. Pressure Vessel Technology,2010,27(5):10-14(in Chinese).

    [13]周亮,唐予遠(yuǎn),聶建斌.多胞織物結(jié)構(gòu)復(fù)合材料靜態(tài)彎曲性能與有限元分析[J].天津工業(yè)大學(xué)學(xué)報(bào),2008,27(5):4-7.

    ZHOU L,TANG Y Y,NIE J B.Experimental and FEM calculation of 3-D glass/polyester resin cellular woven composite under quasistatic three-point bending[J].Journal of Tianjin Polytechnic University,2008,27(5):4-7(in Chinese).

    Mechanical analysis of piercing template with equal center distance and closely aligned tiny hole

    ZHAO Shi-hai1,2,WANG Ting-hui1,2,DONG Jiu-zhi1,2,JIANG Xiu-ming1,2
    (1.School of Mechanical Engineering,Tianjin Polytechnic University,Tianjin 300387,China;2.Tianjin City Key Laboratory of Modern Mechatronics Equipment Technology,Tianjin Polytechnic University,Tianjin 300387,China)

    Piercing template is used in the integrated piercing process of orthogonal laminated woven fabrics.There are nearly a thousand closely aligned small holes,which make it difficult to analyze the mechanical behavior.The effective modulus was obtained from the pseudo-plates theory from two aspects,that is,homogenization theory and elastic theory,respectively.And the characteristics of each method are discussed.According to the engineering requirement,the equivalent plate thickness and the natural frequency were calculated by using the isoeffective modulus.

    piercing template;equal center distance and closely aligned tiny hole;homogenization theory;isoeffective

    TH112

    A

    1671-024X(2016)03-0073-05

    10.3969/j.issn.1671-024x.2016.03.014

    2015-12-15

    天津市科技支撐重點(diǎn)計(jì)劃項(xiàng)目(15ZCZDGX00840)

    趙世海(1970—),男,博士,副教授,主要研究方向?yàn)樾滦图徔棛C(jī)械設(shè)計(jì).E-mai:tjshzhao@163.com

    猜你喜歡
    單胞等距細(xì)觀
    基于NURBS的點(diǎn)陣材料參數(shù)化建模方法
    復(fù)合材料周期結(jié)構(gòu)數(shù)學(xué)均勻化方法的一種新型單胞邊界條件
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動(dòng)彈性模量和阻尼比試驗(yàn)研究
    地震研究(2021年1期)2021-04-13 01:05:24
    擬凸Hartogs域到復(fù)空間形式的全純等距嵌入映射的存在性
    考慮界面層影響的三維機(jī)織復(fù)合材料單胞模型研究
    保持算子束部分等距的映射
    等距延拓以及相關(guān)問(wèn)題
    基于四叉樹(shù)網(wǎng)格加密技術(shù)的混凝土細(xì)觀模型
    PBX炸藥的抗壓強(qiáng)度及抗拉強(qiáng)度細(xì)觀尺度的數(shù)值計(jì)算
    開(kāi)裂混凝土中水分傳輸過(guò)程的細(xì)觀模型
    老司机午夜福利在线观看视频| 亚洲五月天丁香| 波多野结衣高清无吗| 一本综合久久免费| 美女高潮的动态| 国产视频一区二区在线看| 亚洲成av人片在线播放无| 精品一区二区三区av网在线观看| 最近最新免费中文字幕在线| 亚洲av免费在线观看| 国内少妇人妻偷人精品xxx网站 | 中文字幕最新亚洲高清| 久久性视频一级片| 五月玫瑰六月丁香| 五月伊人婷婷丁香| 亚洲国产精品合色在线| 精品久久久久久久久久久久久| 国产淫片久久久久久久久 | 国产97色在线日韩免费| 久久久久久久精品吃奶| 日韩欧美三级三区| 老熟妇仑乱视频hdxx| 亚洲专区字幕在线| 麻豆一二三区av精品| 亚洲男人的天堂狠狠| 一边摸一边抽搐一进一小说| 久久久久久国产a免费观看| 午夜福利欧美成人| 日本免费一区二区三区高清不卡| 超碰成人久久| 成人亚洲精品av一区二区| 麻豆av在线久日| 国产欧美日韩精品亚洲av| 亚洲国产看品久久| 三级国产精品欧美在线观看 | 国产黄a三级三级三级人| 国产精品香港三级国产av潘金莲| 一区二区三区高清视频在线| 成人精品一区二区免费| 精品国产乱子伦一区二区三区| 国产成人啪精品午夜网站| 免费在线观看亚洲国产| 岛国在线免费视频观看| 久久精品影院6| 99热只有精品国产| 搡老岳熟女国产| 91字幕亚洲| 免费无遮挡裸体视频| 国产精品电影一区二区三区| 熟女少妇亚洲综合色aaa.| 国产欧美日韩一区二区精品| 两个人视频免费观看高清| 99精品欧美一区二区三区四区| 久久欧美精品欧美久久欧美| 午夜两性在线视频| 麻豆久久精品国产亚洲av| 免费观看精品视频网站| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利视频1000在线观看| 国产伦人伦偷精品视频| 国产精品一区二区三区四区久久| 午夜福利18| 夜夜看夜夜爽夜夜摸| 国产成人精品久久二区二区免费| 亚洲欧洲精品一区二区精品久久久| 男插女下体视频免费在线播放| 国产精品久久视频播放| 90打野战视频偷拍视频| 在线永久观看黄色视频| 91av网一区二区| av中文乱码字幕在线| 亚洲一区二区三区色噜噜| 欧美精品啪啪一区二区三区| www.熟女人妻精品国产| 99久久精品国产亚洲精品| 90打野战视频偷拍视频| 两个人的视频大全免费| 看片在线看免费视频| www国产在线视频色| 一进一出抽搐动态| 在线观看午夜福利视频| www.www免费av| 在线观看美女被高潮喷水网站 | 最近最新中文字幕大全免费视频| 搡老妇女老女人老熟妇| 亚洲18禁久久av| av欧美777| 欧美日本亚洲视频在线播放| 亚洲专区中文字幕在线| 精品不卡国产一区二区三区| 国产亚洲精品一区二区www| 两个人看的免费小视频| 老司机午夜福利在线观看视频| 久久久色成人| 国产欧美日韩一区二区精品| 亚洲精品中文字幕一二三四区| 99热只有精品国产| 女人高潮潮喷娇喘18禁视频| 一区二区三区激情视频| 午夜福利高清视频| 亚洲熟女毛片儿| www日本在线高清视频| 日韩欧美在线乱码| 99精品在免费线老司机午夜| 国产69精品久久久久777片 | 成人永久免费在线观看视频| 怎么达到女性高潮| 久久久精品大字幕| 法律面前人人平等表现在哪些方面| 日本黄色片子视频| 两人在一起打扑克的视频| 别揉我奶头~嗯~啊~动态视频| 国产亚洲精品久久久com| 久久久久久大精品| 深夜精品福利| 日韩欧美三级三区| 国产精品免费一区二区三区在线| 97人妻精品一区二区三区麻豆| 精品电影一区二区在线| 欧美三级亚洲精品| 色老头精品视频在线观看| 国产高清videossex| 亚洲第一电影网av| www.熟女人妻精品国产| 中出人妻视频一区二区| 麻豆成人午夜福利视频| www.www免费av| 最近最新免费中文字幕在线| 久9热在线精品视频| 成熟少妇高潮喷水视频| 亚洲男人的天堂狠狠| 99热这里只有是精品50| 欧美国产日韩亚洲一区| 丰满人妻熟妇乱又伦精品不卡| 亚洲avbb在线观看| 特大巨黑吊av在线直播| 少妇的丰满在线观看| 特大巨黑吊av在线直播| 99久久精品国产亚洲精品| 97超级碰碰碰精品色视频在线观看| 黄色片一级片一级黄色片| 中亚洲国语对白在线视频| 69av精品久久久久久| 午夜免费激情av| 操出白浆在线播放| 日韩有码中文字幕| av女优亚洲男人天堂 | 成人特级黄色片久久久久久久| 日韩精品中文字幕看吧| 午夜福利在线观看免费完整高清在 | 国产私拍福利视频在线观看| 精品国产三级普通话版| 精华霜和精华液先用哪个| 五月伊人婷婷丁香| 18禁国产床啪视频网站| 丰满的人妻完整版| 亚洲av熟女| 亚洲国产精品成人综合色| 少妇裸体淫交视频免费看高清| 日本免费a在线| 亚洲精品美女久久av网站| 少妇裸体淫交视频免费看高清| 夜夜爽天天搞| 99久久无色码亚洲精品果冻| 亚洲国产精品成人综合色| 免费人成视频x8x8入口观看| 亚洲五月天丁香| 国产精品99久久99久久久不卡| 亚洲av五月六月丁香网| 亚洲国产欧美网| 不卡av一区二区三区| 成人高潮视频无遮挡免费网站| 日韩精品中文字幕看吧| 婷婷精品国产亚洲av在线| 人人妻人人澡欧美一区二区| 午夜福利在线在线| 美女大奶头视频| 精品国产美女av久久久久小说| 美女黄网站色视频| 观看美女的网站| 国产精品久久久久久人妻精品电影| 日日干狠狠操夜夜爽| 色播亚洲综合网| 91久久精品国产一区二区成人 | 亚洲成人中文字幕在线播放| 美女免费视频网站| 97超视频在线观看视频| 99热只有精品国产| 亚洲九九香蕉| 国产精华一区二区三区| 日本黄色片子视频| 国产97色在线日韩免费| 久久午夜综合久久蜜桃| 国产精品野战在线观看| 亚洲专区字幕在线| 免费无遮挡裸体视频| av天堂中文字幕网| 99视频精品全部免费 在线 | 亚洲国产色片| 欧美一级毛片孕妇| 一本精品99久久精品77| 亚洲av成人一区二区三| 国内精品美女久久久久久| 亚洲av成人不卡在线观看播放网| 97人妻精品一区二区三区麻豆| 91麻豆精品激情在线观看国产| 国产伦人伦偷精品视频| 一边摸一边抽搐一进一小说| x7x7x7水蜜桃| 91麻豆av在线| 丰满人妻熟妇乱又伦精品不卡| 亚洲中文日韩欧美视频| 欧美激情久久久久久爽电影| 免费av不卡在线播放| 黄色 视频免费看| 欧美日韩黄片免| 亚洲美女黄片视频| 欧美日韩中文字幕国产精品一区二区三区| 噜噜噜噜噜久久久久久91| 成人无遮挡网站| 午夜精品久久久久久毛片777| 色综合欧美亚洲国产小说| 国产aⅴ精品一区二区三区波| 欧美最黄视频在线播放免费| 一区二区三区国产精品乱码| 9191精品国产免费久久| 国产亚洲精品久久久久久毛片| 久久精品人妻少妇| 国产精品久久久久久精品电影| 19禁男女啪啪无遮挡网站| 两个人视频免费观看高清| 精品熟女少妇八av免费久了| 在线免费观看不下载黄p国产 | 91老司机精品| 天天一区二区日本电影三级| 又黄又爽又免费观看的视频| 1000部很黄的大片| 国产主播在线观看一区二区| 免费高清视频大片| 国产高清有码在线观看视频| 久久九九热精品免费| 成熟少妇高潮喷水视频| 亚洲一区二区三区不卡视频| 亚洲成人免费电影在线观看| 亚洲真实伦在线观看| 精品乱码久久久久久99久播| 91久久精品国产一区二区成人 | 极品教师在线免费播放| 亚洲国产精品成人综合色| 欧美日韩福利视频一区二区| 特大巨黑吊av在线直播| 国产亚洲精品综合一区在线观看| 老汉色∧v一级毛片| 亚洲乱码一区二区免费版| 亚洲,欧美精品.| 亚洲精品久久国产高清桃花| 国产视频内射| 久久精品国产99精品国产亚洲性色| 久久久精品大字幕| 岛国在线观看网站| 此物有八面人人有两片| 久久婷婷人人爽人人干人人爱| 日本免费a在线| 国产精品亚洲av一区麻豆| 中文字幕久久专区| 三级男女做爰猛烈吃奶摸视频| 欧美日本亚洲视频在线播放| 精品乱码久久久久久99久播| 岛国在线免费视频观看| 国产美女午夜福利| 中文在线观看免费www的网站| 国产精品亚洲一级av第二区| 色噜噜av男人的天堂激情| 一进一出好大好爽视频| 性色avwww在线观看| 久久伊人香网站| 人妻夜夜爽99麻豆av| 视频区欧美日本亚洲| 两人在一起打扑克的视频| 国产精品永久免费网站| 1000部很黄的大片| 最近视频中文字幕2019在线8| 法律面前人人平等表现在哪些方面| 欧美日韩黄片免| 国内少妇人妻偷人精品xxx网站 | 国产精品亚洲av一区麻豆| 日日摸夜夜添夜夜添小说| 在线观看免费视频日本深夜| 久久久国产欧美日韩av| 国产一区二区激情短视频| 国产精品久久视频播放| 久久久国产成人精品二区| 天天一区二区日本电影三级| 日韩精品青青久久久久久| 亚洲精品美女久久av网站| 亚洲aⅴ乱码一区二区在线播放| 欧美色视频一区免费| 亚洲国产精品久久男人天堂| 成人精品一区二区免费| 国语自产精品视频在线第100页| 久久欧美精品欧美久久欧美| 国产三级黄色录像| 麻豆成人午夜福利视频| 国产成人av激情在线播放| 欧美不卡视频在线免费观看| 搡老岳熟女国产| 国产精品九九99| 熟女人妻精品中文字幕| 精品一区二区三区av网在线观看| 白带黄色成豆腐渣| 日韩欧美一区二区三区在线观看| 嫩草影院入口| 精品日产1卡2卡| 欧美乱色亚洲激情| h日本视频在线播放| 亚洲国产看品久久| 性色av乱码一区二区三区2| 精品久久蜜臀av无| 国产精品野战在线观看| 人人妻人人澡欧美一区二区| 综合色av麻豆| 夜夜夜夜夜久久久久| 免费看光身美女| 大型黄色视频在线免费观看| 国产黄a三级三级三级人| 成年人黄色毛片网站| 欧美最黄视频在线播放免费| 婷婷精品国产亚洲av在线| 亚洲性夜色夜夜综合| 午夜福利在线观看吧| 亚洲第一电影网av| 国产亚洲av高清不卡| 国产私拍福利视频在线观看| 毛片女人毛片| 国产欧美日韩精品一区二区| 免费看a级黄色片| 性色av乱码一区二区三区2| 国产单亲对白刺激| 搞女人的毛片| 最好的美女福利视频网| 欧美一区二区国产精品久久精品| 69av精品久久久久久| 欧美国产日韩亚洲一区| 色老头精品视频在线观看| 国产精品一及| 后天国语完整版免费观看| 国产探花在线观看一区二区| 一进一出抽搐gif免费好疼| 国产成年人精品一区二区| 日韩高清综合在线| 国产av一区在线观看免费| 亚洲色图av天堂| 国产精品av视频在线免费观看| 18美女黄网站色大片免费观看| 很黄的视频免费| 亚洲一区二区三区不卡视频| 97人妻精品一区二区三区麻豆| 免费大片18禁| 成人亚洲精品av一区二区| 天堂av国产一区二区熟女人妻| 日本五十路高清| 大型黄色视频在线免费观看| 久久久久久久久久黄片| 又大又爽又粗| 国产又黄又爽又无遮挡在线| 小说图片视频综合网站| 亚洲国产精品sss在线观看| av视频在线观看入口| 亚洲欧洲精品一区二区精品久久久| 老熟妇乱子伦视频在线观看| 美女午夜性视频免费| 欧美性猛交╳xxx乱大交人| 美女 人体艺术 gogo| 免费观看精品视频网站| ponron亚洲| 免费看美女性在线毛片视频| aaaaa片日本免费| h日本视频在线播放| 美女黄网站色视频| 美女 人体艺术 gogo| 欧美日韩国产亚洲二区| 国产av在哪里看| 狠狠狠狠99中文字幕| 级片在线观看| 夜夜看夜夜爽夜夜摸| 国产成人福利小说| 欧美丝袜亚洲另类 | 亚洲欧美精品综合久久99| 成人永久免费在线观看视频| 国产av在哪里看| 黄频高清免费视频| 精品日产1卡2卡| 黄色成人免费大全| 他把我摸到了高潮在线观看| 亚洲欧美日韩卡通动漫| 校园春色视频在线观看| 精品福利观看| 在线永久观看黄色视频| 成人性生交大片免费视频hd| 亚洲在线观看片| 成人国产一区最新在线观看| 国产精品一区二区免费欧美| 老司机在亚洲福利影院| 亚洲av免费在线观看| 在线观看免费午夜福利视频| 欧美绝顶高潮抽搐喷水| 美女午夜性视频免费| 色播亚洲综合网| 日韩欧美国产在线观看| 日本免费一区二区三区高清不卡| 国产精品女同一区二区软件 | 国产激情久久老熟女| 成年女人看的毛片在线观看| 亚洲精品粉嫩美女一区| 亚洲成a人片在线一区二区| 99国产精品一区二区三区| 亚洲av成人av| 桃色一区二区三区在线观看| 免费看光身美女| 哪里可以看免费的av片| 亚洲一区高清亚洲精品| 一级黄色大片毛片| АⅤ资源中文在线天堂| 99久久综合精品五月天人人| 亚洲中文字幕一区二区三区有码在线看 | 淫秽高清视频在线观看| 少妇的丰满在线观看| 国产精品98久久久久久宅男小说| 免费高清视频大片| 狂野欧美激情性xxxx| 夜夜看夜夜爽夜夜摸| 亚洲av电影在线进入| www.熟女人妻精品国产| 亚洲人成网站在线播放欧美日韩| 国产av一区在线观看免费| 婷婷精品国产亚洲av| 99精品欧美一区二区三区四区| 在线观看日韩欧美| 少妇熟女aⅴ在线视频| 99在线人妻在线中文字幕| 亚洲国产精品久久男人天堂| 好男人电影高清在线观看| 精品久久久久久久末码| 国产精品电影一区二区三区| 精品国产三级普通话版| 精品熟女少妇八av免费久了| 亚洲欧美精品综合一区二区三区| 国产精品日韩av在线免费观看| 国产高清videossex| 人人妻人人澡欧美一区二区| 99riav亚洲国产免费| 亚洲国产欧美网| 亚洲av免费在线观看| 免费看日本二区| 老熟妇仑乱视频hdxx| 18美女黄网站色大片免费观看| 9191精品国产免费久久| 午夜免费观看网址| 女人被狂操c到高潮| 老熟妇仑乱视频hdxx| 日本 欧美在线| 男女那种视频在线观看| 无遮挡黄片免费观看| 欧美乱妇无乱码| 亚洲色图av天堂| 淫秽高清视频在线观看| 欧美一级a爱片免费观看看| 丁香欧美五月| 久久这里只有精品19| 精品久久蜜臀av无| 国产精品久久久久久亚洲av鲁大| 亚洲五月天丁香| 国产亚洲精品综合一区在线观看| 国产伦精品一区二区三区视频9 | 亚洲精品在线美女| 亚洲熟妇熟女久久| 欧美日韩亚洲国产一区二区在线观看| 一区二区三区高清视频在线| 国产成人影院久久av| 免费看光身美女| 亚洲av成人av| 国产精品久久电影中文字幕| 亚洲国产精品sss在线观看| 麻豆国产av国片精品| 久久久久国产精品人妻aⅴ院| 露出奶头的视频| 看免费av毛片| 亚洲激情在线av| 特大巨黑吊av在线直播| 日韩精品青青久久久久久| 成人性生交大片免费视频hd| 全区人妻精品视频| 最近最新免费中文字幕在线| 手机成人av网站| 久久中文字幕人妻熟女| 真实男女啪啪啪动态图| 黄色 视频免费看| 久久精品91无色码中文字幕| 成人永久免费在线观看视频| 亚洲精品美女久久av网站| 久久人妻av系列| 国产主播在线观看一区二区| 无人区码免费观看不卡| 精品国产超薄肉色丝袜足j| 中国美女看黄片| 老司机在亚洲福利影院| 一夜夜www| 亚洲国产欧美网| 亚洲18禁久久av| 国产精品 欧美亚洲| 亚洲国产精品sss在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲一区二区三区色噜噜| 精品久久久久久久末码| 国产成+人综合+亚洲专区| 亚洲国产欧美人成| 高潮久久久久久久久久久不卡| 国产不卡一卡二| 国产成人啪精品午夜网站| 免费看光身美女| 少妇人妻一区二区三区视频| 久久这里只有精品19| 成人av一区二区三区在线看| 99久久无色码亚洲精品果冻| 欧美日韩国产亚洲二区| 女生性感内裤真人,穿戴方法视频| 级片在线观看| 久久99热这里只有精品18| 精品熟女少妇八av免费久了| 精品国产超薄肉色丝袜足j| 中国美女看黄片| 亚洲色图av天堂| 亚洲aⅴ乱码一区二区在线播放| 搡老熟女国产l中国老女人| 琪琪午夜伦伦电影理论片6080| 成人亚洲精品av一区二区| 欧美又色又爽又黄视频| 两性夫妻黄色片| 小说图片视频综合网站| 亚洲aⅴ乱码一区二区在线播放| 亚洲专区字幕在线| 国产精品自产拍在线观看55亚洲| 成人亚洲精品av一区二区| 欧美日本亚洲视频在线播放| 操出白浆在线播放| 每晚都被弄得嗷嗷叫到高潮| 成年人黄色毛片网站| 麻豆成人av在线观看| 夜夜爽天天搞| h日本视频在线播放| 亚洲欧美日韩卡通动漫| 欧美乱色亚洲激情| 久久久精品欧美日韩精品| 日日夜夜操网爽| 美女高潮喷水抽搐中文字幕| 99热这里只有是精品50| 国产精品免费一区二区三区在线| 午夜成年电影在线免费观看| 男女下面进入的视频免费午夜| 成人三级做爰电影| 亚洲成人中文字幕在线播放| 欧美3d第一页| 18禁黄网站禁片免费观看直播| 亚洲第一欧美日韩一区二区三区| 日韩 欧美 亚洲 中文字幕| e午夜精品久久久久久久| 国产精品影院久久| 真人一进一出gif抽搐免费| 免费av不卡在线播放| 老鸭窝网址在线观看| 婷婷亚洲欧美| 久久久久久久精品吃奶| 99久久99久久久精品蜜桃| 男女下面进入的视频免费午夜| 99久国产av精品| 久久国产乱子伦精品免费另类| 国内精品久久久久久久电影| 亚洲 国产 在线| 日韩有码中文字幕| 一二三四在线观看免费中文在| 国产成人影院久久av| 好男人电影高清在线观看| 久久久久性生活片| 女人高潮潮喷娇喘18禁视频| 成年免费大片在线观看| 久久精品人妻少妇| 在线观看免费午夜福利视频| 日本免费一区二区三区高清不卡| 亚洲一区高清亚洲精品| 亚洲成人中文字幕在线播放| 制服人妻中文乱码| a级毛片a级免费在线| 曰老女人黄片| 禁无遮挡网站| 亚洲av成人不卡在线观看播放网| 欧美3d第一页| 久久99热这里只有精品18| 成人一区二区视频在线观看| 日韩欧美免费精品| 在线观看66精品国产| 一级作爱视频免费观看| 村上凉子中文字幕在线| 色综合欧美亚洲国产小说| 欧美黄色片欧美黄色片| 看黄色毛片网站| 日日干狠狠操夜夜爽| 美女被艹到高潮喷水动态| 叶爱在线成人免费视频播放| 无限看片的www在线观看| www.www免费av| 日本五十路高清| 国产久久久一区二区三区| 变态另类丝袜制服| 中出人妻视频一区二区| 一二三四社区在线视频社区8|