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

    細編穿刺織物的數學建模理論與實例

    2017-04-19 08:09:04周鈺博李艷霞顧軼卓張佐光宋永忠余立瓊
    材料工程 2017年4期
    關鍵詞:貝塞爾單胞周期性

    周鈺博,李艷霞,李 敏,顧軼卓,張佐光,宋永忠,余立瓊,程 家

    (1 北京航空航天大學 材料科學與工程學院 空天材料與服役教育部重點實驗室,北京 100191;2 航天材料及工藝研究所,北京 100076)

    細編穿刺織物的數學建模理論與實例

    周鈺博1,李艷霞1,李 敏1,顧軼卓1,張佐光1,宋永忠2,余立瓊2,程 家2

    (1 北京航空航天大學 材料科學與工程學院 空天材料與服役教育部重點實驗室,北京 100191;2 航天材料及工藝研究所,北京 100076)

    針對細編穿刺織物,采用具有參數連續(xù)性的樣條線作為纖維絲束的軌跡特性函數,建立了三維立體織物結構的數學抽象與三維建模方法?;贛icroCT方法實際測量的細編穿刺織物結構特征參數,采用四種絲束軌跡特性函數和兩種絲束截面建立了八種穿刺織物結構的三維數字模型,通過絲束結構與孔隙率的對比,三維數字模型與織物真實結構符合較好。該織物結構數學抽象與三維建模方法,可用于小尺度的亞單胞模型與大尺度的宏觀模型分析,具有較高的適應性。

    細編穿刺織物;預成型體;數學抽象;計算機建模

    細編穿刺織物是一種三維立體織物,與其他織物相比[1],厚度方向纖維含量很高,因而在厚度方向具有獨特的性能,常作為C/C復合材料的增強體[2,3],用于飛機剎車裝置、火箭耐高溫部件等處。細編穿刺織物的加工工藝復雜,成本高,通過建模仿真方法預測其力學及物理特性[4](如滲透率等)具有重要的意義。目前,針對細編穿刺織物,僅在1~2束纖維的亞單胞尺度有數學模型的報道[5],由于其損失了織物的周期性結構,無法考慮到纖維滑動等因素,模擬仿真的準確性受到影響,而在單胞及若干的單胞的較大尺度上,未見有仿真模型的報道,僅有抽象結構模型[4]。本工作針對細編穿刺織物,建立了一種準確度較高的數學抽象與織物結構數學建模方法,可同時用于小尺度的亞單胞模型與大尺度的宏觀模型分析,具有較高的適應性。

    織物結構建模可以上溯到20世紀30年代,Peirce[6]提出經典的平紋織物的幾何模型,其將絲束的截面假定為橢圓形,并用絲束屈曲高度、絲束幾何密度、絲束紡織角等參數描述織物的模型。近年來,隨著預成型體編織技術和液體成型工藝的發(fā)展,預成型體織物結構的數字化、模型化受到關注。基于單胞模型的各項研究得到展開,如三維立體織物的結構優(yōu)化設計[7]與預成型體中纖維彎曲對復合材料性能的影響分析[8,9];基于含有若干單胞的三維縫合織物的模型,采用計算流體力學和有限元方法,預測三維縫合織物層內的滲透率[10,11]與預成型體的拉伸性能[10];采用有限單元方法,預測預成型體的導電性與導熱性[12],液體成型過程中預成型體受力分析[13],基于單胞模型的預成型體纖維體積分數和厚度預測[14-16]等。與傳統(tǒng)的跑道型建模方法[15]不同,本工作采用具有參數連續(xù)性的樣條線建立細編織物穿刺織物單胞模型,同時,對比了不同絲束截面與軌跡特性所建立織物數學模型,通過絲束結構與孔隙率的對比,織物數學模型與織物真實結構具有較好的符合性,為細編穿刺織物的滲透率預測及液體成型工藝模擬提供技術基礎。

    1 數學建模用織物樣例

    采用的細編穿刺織物由南京玻璃纖維研究院生產,纖維型號為T700碳纖維,通過Skyscan 1076 MicroCT獲取的斷面圖如圖1所示??椢飬禐椋篫向纖維束直徑0.603 mm±0.031mm,XY向纖維束高0.280 mm±0.023mm,XY向纖維束寬0.575 mm±0.046mm,孔隙率為46.67%。

    圖1 細編穿刺織物的MicroCT照片F(xiàn)ig.1 MicroCT photo of fine weave pierced fabric

    2 細編穿刺織物的數學抽象

    對于細編穿刺織物進行三維仿真建模,需將織物進行數學抽象,提取出織物的絲束軌跡特性與截面特性,再將截面特性與絲束軌跡特性結合,獲得絲束實體或表面的特性。針對細編穿刺織物,本文的織物仿真模型精確到絲束級別,總尺度為若干個絲束,為編織結構的若干個單胞組合。

    2.1 絲束軌跡特性

    為了定義絲束的軌跡,將三維絲束實體抽象為一維絲束的中心線,通過設定絲束中心線通過的“特殊點”(即絲束在織物中的“結點”),運用數學插值的方法,實現(xiàn)絲束中心線的連續(xù)性模擬。

    絲束的中心線是一條光滑、連續(xù)的曲線,其函數表達需要滿足至少一階參數連續(xù)性,因此,采用“樣條線”來對絲束的中心線進行擬合。樣條線是通過指定點的擬合曲線,對其進行描述的插值函數具有定義域上的連續(xù)性。樣條線方程通常用分段多項式方程來描述,對于一條定義域為[a,b]的樣條線,典型方程形式如下:

    a=t0

    (1)

    其中,給定k值,ti點為樣條線的結點,也是絲束在織物中的結點,Si(t)為樣條線的分段。

    本工作考慮了兩類四種樣條線:周期性與非周期性貝塞爾樣條線(具有一階參數連續(xù)性)、周期性與非周期性自然立方樣條線(具有二階參數連續(xù)性)。為了避免龍格現(xiàn)象,樣條線的分段參數方程采用三階:

    S(t)=P1(1-t)3+3P2(1-t)2t+

    3P3(1-t)t2+4P4t30≤t≤1

    (2)

    2.1.1 貝塞爾樣條線的參數定義

    貝塞爾樣條線的分段是按照零階連續(xù)性和一階連續(xù)性定義的,按照樣條線結點要求,有:

    P1i=Si(ti) 0≤i≤k-2

    P4i=Si(ti+1) 0≤i≤k-2

    (3)

    且滿足:

    (4)

    在樣條線兩端(t=a與t=b)滿足

    (1)若絲束不具有周期性,有:

    (5)

    結合式(3)和(4),就得到了(非周期性)貝塞爾樣條線。

    (2)若絲束是周期性的,滿足:

    (6)

    結合式(3)和(4),就得到了周期性貝塞爾樣條線。

    通常采用貝塞爾樣條線的模型與真實情況符合較好,同時在特殊情況下,通過提供特定的定義點使曲線更符合真實的絲束路徑。

    2.1.2 自然立方樣條線的參數定義

    與貝塞爾樣條線相似,通過樣條線零階、一階與二階參數連續(xù)性能得出以下方程:

    Si(ti+1)=Si+1(ti+1) 0≤i≤k-3

    S″i(ti+1)=S″i+1(ti+1) 0≤i≤k-3

    (7)

    對于自然立方樣條線,要求使得以下函數值最?。?/p>

    (8)

    函數J(S)的意義是樣條線S(t)的總曲率。因為自然立方樣條線的物理模型是彈性繩(樣條線)在若干個固定點(樣條線結點)的約束下處于總能量最低的狀態(tài),因此要使得樣條線具有最低的總曲率。

    (1)若絲束不具有周期性,則在樣條線端頭處(t=a與t=b),滿足:

    S″(a)=S″(b)=0

    (9)

    結合式(7)和(8),就得到了(非周期性)自然立方樣條線。

    (2)若針對織物的單胞建模,絲束具有周期性,則在樣條線端頭處(t=a與t=b),滿足:

    S′(a)=S′(b)S″(a)=S″(b)

    (10)

    結合式(7)和(8),就得到了周期性自然立方樣條線。由于S′(tk-1)=S′(t0),周期性自然立方樣條線的方程比(非周期性)自然立方樣條線少一個。

    由于在樣條線兩端具有相同的一階與二階導數,周期性自然立方樣條線能夠非常好的適應織物結構單元的建模,在結構單元重復堆積的地方,能夠同時滿足零階、一階與二階參數連續(xù)性,因此對于織物重復單元建模有著重要的意義。

    四種樣條線的對比如圖2所示。

    圖2 四種樣條線的對比Fig.2 Comparison of four kinds of spline

    2.2 絲束截面特性

    絲束的橫截面是二維的平面圖形,所在平面與絲束中心線垂直。絲束是含有若干單絲的三維實體,其橫截面為垂直于絲束中心線,能包含所有纖維的最小平面區(qū)域。由于真實絲束的橫截面外凸且邊緣不光滑,在建模中需將橫截面的最小區(qū)域的邊緣光滑化,在滿足精確性與便捷性的同時做出近似。

    通常,絲束的橫截面有兩種近似:冪次橢圓與透鏡形。

    2.2.1 冪次橢圓截面

    冪次橢圓截面是橢圓形截面的一種修正,通過一個冪次參數n對橢圓在高度h上進行修正,n越小,橫截面越接近于矩形,n越大,橫截面越接近于透鏡形,當冪次系數n=1時,冪次橢圓還原為普通橢圓,如圖3(a)所示。其數學表達式如下:

    (11)

    2.2.2 透鏡形截面

    透鏡形截面是半徑為r1與r2的兩個圓在垂直方向上各自偏移o1與o2所相交形成的,參數r1,r2,o1,o2可以通過截面的寬w、高h以及透鏡形截面的扭曲距離d計算而得,解析表達式為:

    (12)

    其中

    (13)

    圖3(b)所示是三種典型的透鏡形,其中,當透鏡形截面的扭曲距離d為0,兩個圓的半徑相同,偏移距離大小相同,方向相反。

    圖3 冪次橢圓(a)與透鏡形(b)示意圖Fig.3 Power ellipse(a) and lenticular(b)

    2.3 絲束實體與表面特性

    根據數學理論,由絲束中心線和絲束橫截面方程即可確定絲束的三維結構,但在實際建模中,需要繪制絲束實體或表面的三維圖形,而數學抽象的橫截面,是零厚度的二維平面,若簡單地將絲束中心線與絲束橫截面簡單合并起來,需要無數橫截面堆疊,此時數據量巨大。為簡化模型的建立,利用絲束的特性,將建模需要的橫截面減少為有限個。

    沿著絲束的中心線方向,絲束的橫截面是一個光滑變化的圖形,即絲束的橫截面是一個沿著絲束中心線變化的函數。固絲束實體可以定義為一個關于絲束路徑S與絲束橫截面C的函數P(u,v):

    (14)

    (15)

    在參數表面的函數P(u,v)的定義式中,橫截面C由兩個參數確定,分別為變量u與變量v,其中,u代表的是橫截面與絲束中心線的交點在中心線上的位置,v代表的是橫截面上的位置。

    圖4 建模坐標系Fig.4 Coordinate system of modeling

    由于絲束的橫截面是一個不斷變化的圖形,所以采用關鍵點與關鍵面的方法簡化建模:離散地在絲束的中心線上定義一些關鍵點,在關鍵點上定義該點的橫截面(關鍵面),在關鍵點之間采用插值運算的方法,補充橫截面,生成絲束的外表面。通常關鍵點為絲束中心線(樣條線)的結點,為了提高精確度,關鍵點也可以為任意點。當絲束的橫截面變化非常小時,可只在每條絲束中心線上定義一個關鍵面,簡化建模,提高運算性能。

    采用如下的插值方法:若A(t)與B(t)為定義在關鍵點上的關鍵面,則A(t)與B(t)之間的一個插值面C(t)滿足:

    C(t,μ)=A(t)+[B(t)-A(t)]μ

    0≤t≤1 0≤μ≤1

    (16)

    式中μ的值從0到1線性變化,代表插值面C(t)在關鍵點A(t)與B(t)之間線性變化。這種線性插值同時適用于相同與不同類型的截面,并生成滿足零階參數連續(xù)性的光滑過度表面。如圖5(a)所示,圖中顯示了線性插值的截面變化情況。模型從左到右,關鍵面分別為橫向放置的長軸長2單位、短軸長1單位的橢圓形截面,縱向放置的長軸長2單位、短軸長1單位的橢圓形截面以及橫向放置的長2單位、高1單位、扭曲距離為-0.25單位的透鏡型截面。

    線性插值通常與實際情況存在著較大的偏差,如在關鍵點處出現(xiàn)截面的突變,即截面不滿足一階參數連續(xù)性??梢酝ㄟ^適當的平滑修正,改變μ與絲束位置d的線性相關性,修正以上缺陷。可使用三次等式修正,使得在關鍵點處滿足:

    (17)

    即有

    (18)

    圖5(a)經過平滑修正之后的插值如圖5(b)所示。

    圖5 截面插值示意圖(a)線性插值;(b)平滑修正后的插值Fig.5 Sketch of section interpolation(a)linear interpolation;(b)interpolation after smooth

    3 細編穿刺織物的三維仿真模型建立

    基于細編穿刺織物的數學抽象模型,即可按照一定的方法與手段建立織物的三維數字模型。

    分別按照兩種截面,四種絲束軌跡特性曲線進行建模。其中,采用橢圓形截面、周期性貝塞爾樣條線的織物模型的透視圖如圖6所示。

    圖6 采用橢圓形截面、周期性貝塞爾樣條線的織物模型的透視圖Fig.6 Scenograph of models with ellipse section and periodic Bezier spline

    對比八種模型的幾何結構,由于建模所選取的區(qū)域為織物編織結構的單胞,已具有一定的周期重復性,因此對于同一種絲束軌跡的插值方法,樣條線的周期性與否,對幾何上的影響并不明顯,同時,由于絲束軌跡與絲束截面在數學處理上相對獨立,不同的截面形狀并不影響絲束的軌跡。以透鏡形截面為例,對比周期性貝塞爾樣條線與自然立方樣條線,如圖7所示,可以明顯看出,在緊鄰彎折區(qū)域的結點附近(圓圈處),自然立方樣條線出現(xiàn)了明顯的屈曲,而貝塞爾樣條線的屈曲幅度很小,這是由于自然立方樣條線所具有二階參數連續(xù)性而引起的。

    圖7 周期性貝塞爾樣條線(a)與自然立方樣條線(b)的區(qū)別Fig.7 Difference of periodic Bezier spline(a) with periodic natural cubic spline(b)

    將采用橢圓形截面、周期性貝塞爾樣條線的織物模型的斷面圖(如圖8)與真實織物獲取的斷面圖(圖1)進行比較,可以發(fā)現(xiàn),在排除因為三維編織造成的織物不規(guī)則變形之后,建立的模型與真實模型符合較好。

    圖8 采用橢圓形截面、周期性貝塞爾樣條線的織物模型的斷面圖Fig.8 Sectional view of models with ellipse section and periodic Bezier spline

    圖9 不同模型的孔隙率及其與實際值的比較Fig.9 Comparison of porosities from different models and with reality

    基于上述八種織物模型,進一步計算了預成型的孔隙率如圖9所示。由圖9計算數據對比可知:(1)與透鏡形截面的織物模型相比,橢圓形截面模型的孔隙率較低且與穿刺織物的孔隙率更接近,表明穿刺織物的絲束截面更接近于橢圓形;(2)同一種樣條線,周期性與非周期性的差別很小,這是由于織物的特性所決定的;(3) 對比真實織物與模型孔隙率,模型的孔隙率也均高于真實織物。

    孔隙率誤差主要由三個方面產生:(1)穿刺織物在制備過程中,由于外力作用,纖維出現(xiàn)大量的變形與損傷,織物整體呈現(xiàn)受擠壓狀態(tài)[5],破損的纖維填充了相當量的孔隙,導致纖維的有序性降低,絲束的規(guī)則性下降,造成孔隙率偏差;(2)模型將絲束作為一個整體處理,僅考慮了束間孔隙,從而造成孔隙率偏差;(3)數學模型的特征尺寸和結構特性來自于真實織物的測量,同時,模型中絲束的截面與軌跡特性也需要近似化處理,因此,實驗測量和數學近似都會造成數學抽象與真實情況之間的誤差,固需要采取多重測量,以及采取多種近似方法進行建模,從中選取最合理的方案,減小誤差。

    4 結論

    (1)與傳統(tǒng)的跑道型軌跡模型相比,本工作中采用的貝塞爾樣條線與自然立方樣條線軌跡模型的出發(fā)點是絲束的參數連續(xù)性。相對于只具有零階參數連續(xù)性的跑道型軌跡模型,貝塞爾樣條線與自然立方樣條線軌跡模型分別具有一階和二階參數連續(xù)性,因而其更接近于自然狀態(tài)下的織物,建立的模型更符合真實的絲束狀態(tài)。

    (2)以細編穿刺三維立體織物為樣例,本工作分別建立了基于兩種典型截面、四種軌跡樣條線的八種模型,并將它們的孔隙率與真實織物的實測孔隙率進行比較,發(fā)現(xiàn)模型與真實織物的孔隙率數據吻合較好,說明所建立的數學模型抽象與建模方法能夠較好地描述織物的真實狀況,為穿刺織物滲透率及力學性能等參數預測提供了重要基礎。

    (3)該織物微觀結構建模方法也可用于其他類型的2D,2.5D及3D織物等預成型體的抽象建模,并可用于預成型體的滲透率、導電性、導熱性等物理性能以及力學性能的預測與分析,可以提高針對預成型體相關性能研究的效率,降低研究的成本。

    [1] 于倩倩,陳剛,鄭志才. 縫合技術在復合材料上的應用及研究進展[J]. 工程塑料應用, 2009,37(5): 85-88.

    YU Q Q, CHEN G, ZHENG Z C. Application and review of research in stiching technology for composites[J]. Engineering Plastics Application, 2009,37(5): 85-88.

    [2] 朱建勛. 細編穿刺織物的結構特點及性能[J]. 宇航材料工藝, 1998,(1): 41-43.

    ZHU J X. The structural characteristics and properties of fine weave pierced fabric[J]. Aerospace Materials & Technology, 1998,(1): 41-43.

    [3] 朱建勛. 碳布整體穿刺織物編織工藝與結構參數優(yōu)化[D]. 南京: 東南大學, 2003.

    ZHU J X. The optimization of knitting technology and structure parameters for integrated piercing carbon fabric[D]. Nanjing: Southeast University, 2003.

    [4] 方國東,梁軍,王裕. 利用有限元方法對細編穿刺復合材料進行有效性能預報[C]// 第十五屆全國復合材料學術會議.哈爾濱: 中國力學協(xié)會, 2008:1397-1401.

    FANG G D,LIANG J,WANG Y. The finite element analysis of mechanical properties of fine woven pierced composites [C]//15th National Conference on Composite Materials.Harbin: Chinese Society of Theoretical and Applied Mechanics, 2008:1397-1401.

    [5] 陳盛洪. 細編穿刺C/C復合材料熱物理性能的模擬研究[D]. 哈爾濱: 哈爾濱工業(yè)大學, 2008.

    CHEN S H. Simulation on thermal physical properties of fine weave pierced C/C composite[D]. Harbin: Harbin Institute of Technology, 2008.

    [6] PEIRCE F T. The geometry of cloth structure[J]. Journal of the Textile Institute Transactions, 1937, 28(3): T45-T96.

    [7] OKUMURA T, YOKOYAMA A, NAGAI K, et al. Optimum design of weaving structure of 3-D woven fabric composites by using genetic algorithms[J]. Composite Structures, 1995, 32(1-4): 417-426.

    [8] STIG F, HALLSTR M S. Influence of crimp on 3D-woven fibre reinforced composites[J]. Composite Structures, 2013, 95: 114-122.

    [9] 燕瑛,成傳賢. 基于細觀結構的三維機織復合材料彈性性能的分析[J]. 航空學報, 1999,28(4): 2-6.

    YAN Y, CHENG C X. Analysis of elastic property for 3-D woven composites based on fabric microstructure[J]. Acta Aeronautica et Astronautica Sinica, 1999,28(4): 2-6.

    [10] ZENG X, BROWN L P, ENDRUWEIT A, et al. Geometrical modelling of 3D woven reinforcements for polymer composites: Prediction of fabric permeability and composite mechanical properties[J]. Composites Part A: Applied Science and Manufacturing, 2014, 56: 150-160.

    [11] 楊波,金天國,鄭龍. 微-細雙尺度單胞下織物預成型體滲透率預測[J]. 復合材料學報, 2013,30(5): 209-217.

    YANG B, JIN T G, ZHENG L. Permeability prediction for textile preform with micro-meso dual-scale unit cell[J]. Acta Materiae Compositae Sinica, 2013,30(5): 209-217.

    [12] SIDDIQUI M O R, SUN D. Finite element analysis of thermal conductivity and thermal resistance behaviour of woven fabric[J]. Computational Materials Science, 2013, 75: 45-51.

    [13] WENDLING A, HIVET G, VIDAL-SALLé E, et al. Consistent geometrical modelling of interlock fabrics[J]. Finite Elements in Analysis and Design, 2014, 90: 93-105.

    [14] 丁辛,易洪雷. 三維機織幾何結構的數值表征[J]. 東華大學學報(自然科學版), 2003,29(3): 15-19.

    DING X, YI H L. Representation of geometric architecture of three-dimensional woven structures by numerical method[J]. Journal of Donghua University(Natural Science),2003,29(3): 15-19.

    [15] 丁辛,易洪雷. 三維機織結構的幾何模型[J]. 復合材料學報, 2003,20(5): 108-113.

    DING X, YI H L. A geometric model of three dimensional woven structures[J]. Acta Materiae Compositae Sinica, 2003,20(5): 108-113.

    [16] 楊彩云,李嘉祿. 基于紗線真實形態(tài)的三維機織復合材料細觀結構及其厚度計算[J]. 復合材料學報, 2005,22(6): 178-182.

    YANG C Y, LI J L. Microstructure and thickness equation of 3D woven composites based on yarn’s true configuration [J]. Acta Materiae Compositae Sinica, 2005,22(6): 178-182.

    (本文責編:解 宏)

    Theory and Examples of Mathematical Modeling for Fine Weave Pierced Fabric

    ZHOU Yu-bo1,LI Yan-xia1,LI Min1,GU Yi-zhuo1,ZHANG Zuo-guang1,SONG Yong-zhong2,YU Li-qiong2,CHENG Jia2

    (1 Key Laboratory of Aerospace Materials and Performance (Ministry of Education), School of Materials Science and Engineering, Beihang University,Beijing 100191,China;2 Aerospace Research Institute of Materials & Processing Technology,Beijing 100076,China)

    A mathematical abstraction and three-dimensional modeling method of three-dimensional woven fabric structure was developed for the fine weave pierced fabric, taking parametric continuity splines as the track function of tow. Based on the significant parameters of fine weave pierced fabric measured by MicroCT, eight kinds of the three-dimensional digital models of the fabric structure were established with two kinds of tow sections and four kinds of tow trajectory characteristic functions. There is a good agreement between the three-dimensional digital models and real fabric by comparing their structures and porosities. This mathematical abstraction and three-dimensional modeling method can be applied in micro models for sub unit cell and macro models for macroscopic scale fabrics, with high adaptability.

    fine weave pierced fabric;preform;mathematical abstraction;computer modeling

    10.11868/j.issn.1001-4381.2015.000479

    TB33

    A

    1001-4381(2017)04-0102-06

    2015-04-23;

    2016-10-08

    李艷霞(1977-),女,講師,研究方向:樹脂基復合材料成型工藝及數值模擬研究,聯(lián)系地址:北京市海淀區(qū)學院路37號北京航空航天大學材料學院104教研室(100191),E-mail:liyanxia@buaa.edu.cn

    猜你喜歡
    貝塞爾單胞周期性
    基于NURBS的點陣材料參數化建模方法
    復合材料周期結構數學均勻化方法的一種新型單胞邊界條件
    看星星的人:貝塞爾
    少兒科技(2021年3期)2021-01-20 13:18:34
    數列中的周期性和模周期性
    中等數學(2019年5期)2019-08-30 03:52:20
    考慮界面層影響的三維機織復合材料單胞模型研究
    基于虛宗量貝塞爾函數的螺旋帶色散模型
    一類整數遞推數列的周期性
    中等數學(2018年12期)2018-02-16 07:48:40
    基于擴頻碼周期性的單通道直擴通信半盲分離抗干擾算法
    一種脈沖貝塞爾波的構造及其非線性聲場的仿真
    CPSIA周期性測試要求開始實施
    亚洲人与动物交配视频| 特大巨黑吊av在线直播| 日韩欧美在线乱码| 91午夜精品亚洲一区二区三区| 热99在线观看视频| 卡戴珊不雅视频在线播放| 美女高潮的动态| 日日摸夜夜添夜夜添小说| 你懂的网址亚洲精品在线观看 | 国产乱人视频| 露出奶头的视频| 国产淫片久久久久久久久| 久久久久精品国产欧美久久久| 国产亚洲精品av在线| 别揉我奶头~嗯~啊~动态视频| 亚洲熟妇中文字幕五十中出| 成人漫画全彩无遮挡| 国产一区二区在线观看日韩| 午夜福利在线观看免费完整高清在 | 久久午夜福利片| 一级毛片久久久久久久久女| 久久久精品94久久精品| 在线免费观看的www视频| 色综合亚洲欧美另类图片| www.色视频.com| 看非洲黑人一级黄片| 国产精品,欧美在线| 18禁在线播放成人免费| 中出人妻视频一区二区| 精品久久久久久久人妻蜜臀av| 亚洲中文字幕一区二区三区有码在线看| 久久99热这里只有精品18| 日韩欧美免费精品| 特大巨黑吊av在线直播| 亚洲精品一卡2卡三卡4卡5卡| 久久国产乱子免费精品| 少妇猛男粗大的猛烈进出视频 | 老熟妇乱子伦视频在线观看| 精华霜和精华液先用哪个| 欧美一区二区亚洲| 女的被弄到高潮叫床怎么办| 亚洲色图av天堂| 日韩中字成人| 少妇裸体淫交视频免费看高清| 精品无人区乱码1区二区| 亚洲人成网站在线播| 亚洲欧美清纯卡通| 久久久久久久久久黄片| 在线观看免费视频日本深夜| 国产精品三级大全| 亚洲欧美成人综合另类久久久 | 最近2019中文字幕mv第一页| 亚洲最大成人av| 波多野结衣高清作品| 哪里可以看免费的av片| 露出奶头的视频| 亚洲性夜色夜夜综合| 男女之事视频高清在线观看| 日韩欧美免费精品| 俺也久久电影网| 国产精华一区二区三区| 久久九九热精品免费| 观看免费一级毛片| 一区二区三区高清视频在线| 日韩强制内射视频| 极品教师在线视频| 精品不卡国产一区二区三区| 哪里可以看免费的av片| 亚洲一区高清亚洲精品| 国产白丝娇喘喷水9色精品| 色噜噜av男人的天堂激情| 亚洲av第一区精品v没综合| 一级av片app| 国产探花极品一区二区| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 亚洲av不卡在线观看| 成人高潮视频无遮挡免费网站| 国产一区二区在线观看日韩| 久久国内精品自在自线图片| 久久久久精品国产欧美久久久| 欧美人与善性xxx| 99精品在免费线老司机午夜| 欧美潮喷喷水| 午夜视频国产福利| 亚洲性久久影院| 成年版毛片免费区| 在线a可以看的网站| 网址你懂的国产日韩在线| 寂寞人妻少妇视频99o| 亚洲精品在线观看二区| 久久久成人免费电影| 国产 一区精品| 久久精品久久久久久噜噜老黄 | 国产在线男女| 特级一级黄色大片| 免费av不卡在线播放| 日韩制服骚丝袜av| 国产黄色视频一区二区在线观看 | 亚洲人成网站高清观看| 亚洲精品日韩在线中文字幕 | 性欧美人与动物交配| 又黄又爽又免费观看的视频| 久久中文看片网| 国内精品宾馆在线| 国产精品不卡视频一区二区| 国内揄拍国产精品人妻在线| 久99久视频精品免费| 国内精品宾馆在线| 长腿黑丝高跟| 欧美精品国产亚洲| 国产成人福利小说| 久久综合国产亚洲精品| 成人漫画全彩无遮挡| videossex国产| 热99re8久久精品国产| 成熟少妇高潮喷水视频| 精品日产1卡2卡| 真实男女啪啪啪动态图| 国产精品电影一区二区三区| 不卡视频在线观看欧美| 欧美日韩综合久久久久久| 又粗又爽又猛毛片免费看| 久久久成人免费电影| av视频在线观看入口| 中文字幕人妻熟人妻熟丝袜美| 亚洲经典国产精华液单| 国产av在哪里看| 免费人成在线观看视频色| 成年女人毛片免费观看观看9| 欧美潮喷喷水| 黄色一级大片看看| 国产av在哪里看| 久久久久久久久久黄片| 成人特级av手机在线观看| 91在线精品国自产拍蜜月| 精品欧美国产一区二区三| 国产成人福利小说| 亚洲性夜色夜夜综合| 国产精品嫩草影院av在线观看| 色视频www国产| 男女那种视频在线观看| 91狼人影院| 成年版毛片免费区| av天堂在线播放| 精品乱码久久久久久99久播| 最近视频中文字幕2019在线8| 男人和女人高潮做爰伦理| 亚洲第一电影网av| 亚洲av五月六月丁香网| 99国产精品一区二区蜜桃av| 嫩草影院入口| 99在线视频只有这里精品首页| 久久人人爽人人片av| 成年女人毛片免费观看观看9| 天堂网av新在线| 国产精品国产高清国产av| 波多野结衣高清无吗| www日本黄色视频网| 嫩草影院入口| 久久久久免费精品人妻一区二区| 色吧在线观看| 午夜激情欧美在线| 精品福利观看| 亚洲无线在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲激情五月婷婷啪啪| 卡戴珊不雅视频在线播放| 精品久久久噜噜| 蜜臀久久99精品久久宅男| 毛片一级片免费看久久久久| 可以在线观看毛片的网站| 精华霜和精华液先用哪个| 一本久久中文字幕| 搞女人的毛片| 精品午夜福利在线看| 一级黄片播放器| 国产精品国产高清国产av| 欧美精品国产亚洲| 别揉我奶头 嗯啊视频| 国产精品综合久久久久久久免费| 毛片一级片免费看久久久久| 日本五十路高清| 日本免费一区二区三区高清不卡| 18禁在线播放成人免费| 国产日本99.免费观看| 真实男女啪啪啪动态图| 日韩欧美在线乱码| 日本熟妇午夜| 免费在线观看影片大全网站| 最好的美女福利视频网| 哪里可以看免费的av片| 国产精品久久久久久av不卡| 美女高潮的动态| 国产色婷婷99| 小说图片视频综合网站| 精品一区二区三区人妻视频| 亚洲成人精品中文字幕电影| 久久人人爽人人爽人人片va| 亚洲五月天丁香| 成人二区视频| 国产精品不卡视频一区二区| 91麻豆精品激情在线观看国产| 能在线免费观看的黄片| 欧美区成人在线视频| 最近视频中文字幕2019在线8| 久久久久精品国产欧美久久久| 嫩草影视91久久| 女人十人毛片免费观看3o分钟| 亚洲三级黄色毛片| 欧美bdsm另类| 免费在线观看成人毛片| 亚洲一区二区三区色噜噜| 久久综合国产亚洲精品| 国产精品野战在线观看| 免费av毛片视频| 少妇熟女欧美另类| 亚洲成人久久性| 日日摸夜夜添夜夜爱| 联通29元200g的流量卡| 丰满的人妻完整版| 性欧美人与动物交配| 亚洲一区二区三区色噜噜| 两个人视频免费观看高清| 一区二区三区免费毛片| 床上黄色一级片| 国产精华一区二区三区| 久久久国产成人精品二区| 91在线观看av| 午夜福利在线在线| 我要看日韩黄色一级片| 国产美女午夜福利| 亚洲丝袜综合中文字幕| 久久久久国产精品人妻aⅴ院| 此物有八面人人有两片| 久久精品影院6| 嫩草影院入口| 内地一区二区视频在线| 色5月婷婷丁香| 国产三级在线视频| 亚洲精品色激情综合| videossex国产| 麻豆av噜噜一区二区三区| 久久久久久久久中文| 嫩草影视91久久| 淫秽高清视频在线观看| 精品一区二区三区视频在线| 男女那种视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产aⅴ精品一区二区三区波| 神马国产精品三级电影在线观看| 一本一本综合久久| 国产视频一区二区在线看| 成人漫画全彩无遮挡| 久久久精品大字幕| 色综合亚洲欧美另类图片| 成人鲁丝片一二三区免费| 亚洲中文字幕一区二区三区有码在线看| 69人妻影院| 国产在线男女| 麻豆国产97在线/欧美| 国产91av在线免费观看| 日产精品乱码卡一卡2卡三| 性插视频无遮挡在线免费观看| 国产高清三级在线| 内射极品少妇av片p| 99热只有精品国产| 久久精品影院6| 99在线视频只有这里精品首页| 看黄色毛片网站| 成人三级黄色视频| 搡老妇女老女人老熟妇| 91精品国产九色| 日本三级黄在线观看| 国产av不卡久久| 男插女下体视频免费在线播放| 成人永久免费在线观看视频| 九九热线精品视视频播放| 国产在视频线在精品| 无遮挡黄片免费观看| 国产美女午夜福利| av在线蜜桃| 在线观看午夜福利视频| 伦理电影大哥的女人| 亚洲最大成人av| 成人综合一区亚洲| 亚洲国产欧洲综合997久久,| 99热这里只有是精品50| 久久中文看片网| 亚洲久久久久久中文字幕| 亚洲精品一区av在线观看| 亚洲国产精品成人久久小说 | 国内精品美女久久久久久| 一夜夜www| 成人精品一区二区免费| 99在线人妻在线中文字幕| 尤物成人国产欧美一区二区三区| 男人舔女人下体高潮全视频| 亚洲精品456在线播放app| 精品人妻偷拍中文字幕| 成人性生交大片免费视频hd| 国内精品美女久久久久久| 一边摸一边抽搐一进一小说| 小说图片视频综合网站| 18禁在线无遮挡免费观看视频 | 欧美一区二区亚洲| 亚洲无线观看免费| 激情 狠狠 欧美| 国内精品久久久久精免费| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久久久免| 国产三级在线视频| 最新在线观看一区二区三区| 男女啪啪激烈高潮av片| 麻豆成人午夜福利视频| 国产精品一区二区性色av| 国产成人影院久久av| 亚洲经典国产精华液单| 国产精品久久视频播放| 国语自产精品视频在线第100页| 欧美性感艳星| 亚洲,欧美,日韩| 丰满乱子伦码专区| 国产在线男女| 久久精品91蜜桃| av在线蜜桃| 最后的刺客免费高清国语| 熟妇人妻久久中文字幕3abv| 俄罗斯特黄特色一大片| 国产亚洲精品综合一区在线观看| 中文字幕熟女人妻在线| 18+在线观看网站| 国产成人福利小说| 亚洲av中文字字幕乱码综合| 变态另类丝袜制服| 99久久精品一区二区三区| 亚洲综合色惰| 一本精品99久久精品77| 亚洲av免费在线观看| 中文在线观看免费www的网站| 久久精品国产自在天天线| 欧美色欧美亚洲另类二区| 又粗又爽又猛毛片免费看| 人妻少妇偷人精品九色| 97碰自拍视频| 桃色一区二区三区在线观看| 国产三级在线视频| 国产精品伦人一区二区| eeuss影院久久| 干丝袜人妻中文字幕| 亚洲精品在线观看二区| 美女高潮的动态| 最好的美女福利视频网| 亚洲精品456在线播放app| 欧美日韩在线观看h| 最新在线观看一区二区三区| 精华霜和精华液先用哪个| 久久草成人影院| 久久人人精品亚洲av| 日本黄色片子视频| 日韩欧美精品v在线| 久久久久性生活片| 午夜福利18| 在线看三级毛片| 日日啪夜夜撸| 色视频www国产| 亚洲美女黄片视频| 国产精品伦人一区二区| 午夜a级毛片| 老司机福利观看| 91av网一区二区| 久久草成人影院| 日本色播在线视频| 99热这里只有精品一区| 免费av毛片视频| 午夜久久久久精精品| 精品一区二区免费观看| 男女那种视频在线观看| 九色成人免费人妻av| 人妻丰满熟妇av一区二区三区| 亚洲熟妇中文字幕五十中出| 亚洲av中文字字幕乱码综合| 欧美性感艳星| 黄色一级大片看看| 国模一区二区三区四区视频| 不卡一级毛片| 国产人妻一区二区三区在| 日本-黄色视频高清免费观看| 99久久精品国产国产毛片| 亚洲一区二区三区色噜噜| 麻豆国产av国片精品| 人妻夜夜爽99麻豆av| 最近在线观看免费完整版| 身体一侧抽搐| 午夜激情福利司机影院| 少妇人妻一区二区三区视频| 亚洲美女搞黄在线观看 | 亚洲av成人精品一区久久| 日本精品一区二区三区蜜桃| 国产激情偷乱视频一区二区| 免费大片18禁| 色综合色国产| 日韩欧美精品免费久久| 亚洲成人久久性| 国产麻豆成人av免费视频| 国产乱人视频| 欧美中文日本在线观看视频| 久久这里只有精品中国| 国产在视频线在精品| 在线免费观看的www视频| 人妻丰满熟妇av一区二区三区| 国产成人福利小说| 欧美日韩国产亚洲二区| 69人妻影院| 国语自产精品视频在线第100页| 在线观看午夜福利视频| 国产高潮美女av| 非洲黑人性xxxx精品又粗又长| 国产精品福利在线免费观看| 精品日产1卡2卡| 色综合亚洲欧美另类图片| 亚洲人成网站高清观看| 色综合站精品国产| 伦精品一区二区三区| 搡老妇女老女人老熟妇| 一级黄色大片毛片| 亚洲四区av| 久久精品人妻少妇| 亚洲一区高清亚洲精品| 日韩制服骚丝袜av| 免费人成在线观看视频色| 噜噜噜噜噜久久久久久91| 91狼人影院| 亚洲国产精品成人综合色| 欧美最新免费一区二区三区| 99在线视频只有这里精品首页| 日韩精品青青久久久久久| 国产视频内射| 精品国产三级普通话版| 黄色欧美视频在线观看| av在线天堂中文字幕| 一个人看视频在线观看www免费| 国产成年人精品一区二区| 99热全是精品| 亚洲成人中文字幕在线播放| 中文在线观看免费www的网站| 国产欧美日韩精品亚洲av| 在线国产一区二区在线| а√天堂www在线а√下载| av福利片在线观看| 简卡轻食公司| 亚洲精品日韩av片在线观看| 尤物成人国产欧美一区二区三区| 床上黄色一级片| 国产成人一区二区在线| 欧美xxxx黑人xx丫x性爽| av.在线天堂| 国产精品嫩草影院av在线观看| 此物有八面人人有两片| 一卡2卡三卡四卡精品乱码亚洲| 久久中文看片网| 日日摸夜夜添夜夜添小说| 少妇裸体淫交视频免费看高清| 精品福利观看| 别揉我奶头~嗯~啊~动态视频| 国产爱豆传媒在线观看| 六月丁香七月| 国产精品国产高清国产av| 成年女人永久免费观看视频| 淫妇啪啪啪对白视频| 亚洲一级一片aⅴ在线观看| 午夜精品国产一区二区电影 | 国产精品久久久久久久久免| 欧美又色又爽又黄视频| 亚洲精品久久国产高清桃花| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久久av| 在线观看av片永久免费下载| 一夜夜www| 天堂av国产一区二区熟女人妻| 草草在线视频免费看| 欧美成人免费av一区二区三区| 少妇人妻一区二区三区视频| 亚洲天堂国产精品一区在线| 淫秽高清视频在线观看| 三级经典国产精品| 能在线免费观看的黄片| 久久精品国产自在天天线| 亚洲人成网站高清观看| 国产精品女同一区二区软件| 欧美区成人在线视频| 热99re8久久精品国产| 久久精品久久久久久噜噜老黄 | 直男gayav资源| 国产男人的电影天堂91| 内射极品少妇av片p| 成年女人看的毛片在线观看| 免费电影在线观看免费观看| 亚洲中文日韩欧美视频| 毛片一级片免费看久久久久| 国产高清视频在线播放一区| 一级av片app| 精品乱码久久久久久99久播| 国产精品人妻久久久久久| 国产精品久久久久久av不卡| 精品人妻偷拍中文字幕| 国产精品久久电影中文字幕| 久久人人爽人人爽人人片va| 日本熟妇午夜| 久久久国产成人免费| 久久人人精品亚洲av| 中文在线观看免费www的网站| 亚洲精品国产成人久久av| 亚洲色图av天堂| 简卡轻食公司| 国内揄拍国产精品人妻在线| 国产av麻豆久久久久久久| 男女之事视频高清在线观看| 国产成人精品久久久久久| 欧美中文日本在线观看视频| 国产一区亚洲一区在线观看| 搞女人的毛片| 国产欧美日韩一区二区精品| eeuss影院久久| 免费人成在线观看视频色| 成人午夜高清在线视频| 亚洲中文字幕一区二区三区有码在线看| 麻豆av噜噜一区二区三区| 久久中文看片网| 久久久久免费精品人妻一区二区| 俺也久久电影网| 最好的美女福利视频网| 国产色爽女视频免费观看| 免费大片18禁| 国内精品美女久久久久久| 蜜桃亚洲精品一区二区三区| 黑人高潮一二区| 久久久久免费精品人妻一区二区| 国产精华一区二区三区| 美女大奶头视频| 欧美一级a爱片免费观看看| 日韩欧美免费精品| 亚洲一区二区三区色噜噜| 极品教师在线视频| 精品久久久久久久久av| 亚洲av二区三区四区| 午夜激情欧美在线| 香蕉av资源在线| 欧美高清成人免费视频www| 一个人观看的视频www高清免费观看| 午夜福利高清视频| 欧洲精品卡2卡3卡4卡5卡区| 露出奶头的视频| 亚洲人成网站在线播放欧美日韩| 亚洲一区高清亚洲精品| 亚洲国产欧洲综合997久久,| 国产精品国产高清国产av| 国产美女午夜福利| 免费av毛片视频| 久久久久免费精品人妻一区二区| 国产亚洲精品久久久com| 亚洲成人中文字幕在线播放| 国产欧美日韩精品一区二区| 18禁在线无遮挡免费观看视频 | 色av中文字幕| 九九热线精品视视频播放| 女同久久另类99精品国产91| 欧美丝袜亚洲另类| 久久久久久久久久黄片| 少妇熟女aⅴ在线视频| 成人午夜高清在线视频| 国产精品野战在线观看| 一区二区三区免费毛片| 免费看a级黄色片| 精品久久久久久久久久免费视频| av在线播放精品| 大型黄色视频在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美免费精品| 99久久久亚洲精品蜜臀av| 成人亚洲精品av一区二区| 少妇的逼好多水| av福利片在线观看| 成人鲁丝片一二三区免费| 久久人人爽人人爽人人片va| 一个人看的www免费观看视频| 少妇人妻精品综合一区二区 | 国产一区二区亚洲精品在线观看| 国产高清不卡午夜福利| 成人二区视频| 国产一区二区在线av高清观看| 变态另类成人亚洲欧美熟女| 久久久国产成人精品二区| 九色成人免费人妻av| 亚洲成人av在线免费| 久久99热6这里只有精品| 亚洲熟妇中文字幕五十中出| 三级经典国产精品| 色视频www国产| 亚洲内射少妇av| 在线观看免费视频日本深夜| 国产精品久久久久久久久免| 久久久久久久午夜电影| 在线观看免费视频日本深夜| 亚洲欧美精品综合久久99| 综合色av麻豆| 不卡一级毛片| 日本在线视频免费播放| 欧美日韩乱码在线| 男女啪啪激烈高潮av片| 日本一本二区三区精品| 久久韩国三级中文字幕| 村上凉子中文字幕在线| 99久久精品一区二区三区| 春色校园在线视频观看|