• <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周期性測試要求開始實施
    亚洲国产中文字幕在线视频| 日韩欧美一区二区三区在线观看| 窝窝影院91人妻| 天堂√8在线中文| 不卡av一区二区三区| 又爽又黄无遮挡网站| 亚洲九九香蕉| 一边摸一边做爽爽视频免费| 99热6这里只有精品| 日韩精品免费视频一区二区三区| www日本在线高清视频| 可以免费在线观看a视频的电影网站| 国产欧美日韩一区二区三| 亚洲自偷自拍图片 自拍| 国产片内射在线| 人成视频在线观看免费观看| 12—13女人毛片做爰片一| 午夜久久久久精精品| 久久久久国产一级毛片高清牌| 男女床上黄色一级片免费看| 日韩免费av在线播放| 国产精品久久电影中文字幕| 精品欧美国产一区二区三| 18禁黄网站禁片免费观看直播| 99精品在免费线老司机午夜| 很黄的视频免费| 日本免费a在线| 国产一区二区三区视频了| 久久午夜综合久久蜜桃| 两人在一起打扑克的视频| 黄片小视频在线播放| 国产伦在线观看视频一区| 免费在线观看黄色视频的| 51午夜福利影视在线观看| 国产爱豆传媒在线观看 | 身体一侧抽搐| 国产激情久久老熟女| 我要搜黄色片| 日本精品一区二区三区蜜桃| 日本成人三级电影网站| 欧美中文日本在线观看视频| 人成视频在线观看免费观看| 亚洲熟妇熟女久久| 麻豆成人av在线观看| www.精华液| 桃色一区二区三区在线观看| 亚洲18禁久久av| 亚洲熟女毛片儿| 好看av亚洲va欧美ⅴa在| 三级毛片av免费| 男女床上黄色一级片免费看| 神马国产精品三级电影在线观看 | 美女黄网站色视频| 亚洲 国产 在线| 高潮久久久久久久久久久不卡| 99riav亚洲国产免费| 国产伦人伦偷精品视频| 色在线成人网| 91在线观看av| 国产单亲对白刺激| 99在线视频只有这里精品首页| 久久这里只有精品中国| 白带黄色成豆腐渣| 亚洲av成人一区二区三| 麻豆国产av国片精品| 久久久久久久久久黄片| 国产不卡一卡二| 女人爽到高潮嗷嗷叫在线视频| 亚洲人成电影免费在线| 久久久国产精品麻豆| 美女高潮喷水抽搐中文字幕| 又黄又爽又免费观看的视频| 亚洲av电影不卡..在线观看| 精品欧美一区二区三区在线| 亚洲男人的天堂狠狠| 大型av网站在线播放| 中文字幕熟女人妻在线| 午夜福利高清视频| 亚洲av中文字字幕乱码综合| 欧美久久黑人一区二区| 一边摸一边抽搐一进一小说| 国产精品久久久人人做人人爽| 此物有八面人人有两片| 亚洲专区字幕在线| 欧美色欧美亚洲另类二区| 一级片免费观看大全| 国产精品日韩av在线免费观看| 身体一侧抽搐| 人人妻人人看人人澡| 亚洲欧美日韩高清在线视频| 嫁个100分男人电影在线观看| 亚洲成人国产一区在线观看| 少妇人妻一区二区三区视频| 又黄又粗又硬又大视频| 国产单亲对白刺激| 国产97色在线日韩免费| 久久人妻福利社区极品人妻图片| 国产三级在线视频| 两个人的视频大全免费| 亚洲午夜理论影院| 国内精品久久久久久久电影| av视频在线观看入口| 亚洲片人在线观看| 国产精品国产高清国产av| 欧美一级毛片孕妇| 99久久精品国产亚洲精品| 亚洲av电影在线进入| 精品人妻1区二区| 国产精品综合久久久久久久免费| 免费在线观看日本一区| 国产成人一区二区三区免费视频网站| 精品熟女少妇八av免费久了| 日韩三级视频一区二区三区| 2021天堂中文幕一二区在线观| 99国产极品粉嫩在线观看| 国产伦在线观看视频一区| 伦理电影免费视频| 老司机福利观看| 久久精品成人免费网站| 91av网站免费观看| 国产精品一区二区三区四区久久| 亚洲精品中文字幕在线视频| 国产精品亚洲一级av第二区| 国产99久久九九免费精品| 久久这里只有精品19| 嫁个100分男人电影在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲专区中文字幕在线| 妹子高潮喷水视频| 午夜日韩欧美国产| 色av中文字幕| 婷婷丁香在线五月| 这个男人来自地球电影免费观看| 久久久国产精品麻豆| a级毛片a级免费在线| 午夜福利18| 老司机靠b影院| 亚洲va日本ⅴa欧美va伊人久久| 熟女少妇亚洲综合色aaa.| 成人三级做爰电影| 日韩高清综合在线| 久久久久久久午夜电影| 舔av片在线| 国产高清激情床上av| 午夜视频精品福利| 国产午夜福利久久久久久| 丁香六月欧美| 观看免费一级毛片| 成年女人毛片免费观看观看9| 亚洲精华国产精华精| 青草久久国产| 露出奶头的视频| 美女午夜性视频免费| 欧美黑人精品巨大| 两个人的视频大全免费| 日本黄色视频三级网站网址| 在线观看舔阴道视频| 校园春色视频在线观看| 国产激情久久老熟女| 精品欧美国产一区二区三| 午夜亚洲福利在线播放| 免费看十八禁软件| 亚洲欧美激情综合另类| 最近在线观看免费完整版| 欧美在线一区亚洲| 黄色视频不卡| av在线天堂中文字幕| 亚洲精品久久成人aⅴ小说| 亚洲全国av大片| 亚洲 国产 在线| 国产1区2区3区精品| 免费在线观看亚洲国产| a级毛片在线看网站| av国产免费在线观看| 久久人人精品亚洲av| 欧美国产日韩亚洲一区| 一卡2卡三卡四卡精品乱码亚洲| 日韩有码中文字幕| 国产成人精品无人区| 中国美女看黄片| 九九热线精品视视频播放| 国产蜜桃级精品一区二区三区| 色综合亚洲欧美另类图片| 男人舔女人的私密视频| 草草在线视频免费看| 12—13女人毛片做爰片一| 色哟哟哟哟哟哟| 婷婷亚洲欧美| 真人做人爱边吃奶动态| 在线十欧美十亚洲十日本专区| 国产在线精品亚洲第一网站| 制服人妻中文乱码| 国产精华一区二区三区| 免费观看人在逋| av超薄肉色丝袜交足视频| 九色成人免费人妻av| 99在线人妻在线中文字幕| av欧美777| 亚洲真实伦在线观看| 久久精品夜夜夜夜夜久久蜜豆 | 长腿黑丝高跟| 国产一区二区在线观看日韩 | 国产精品久久久久久精品电影| 18禁国产床啪视频网站| 国产精品久久久av美女十八| 俄罗斯特黄特色一大片| av福利片在线| 欧美大码av| 中文亚洲av片在线观看爽| 一级a爱片免费观看的视频| 成在线人永久免费视频| 两人在一起打扑克的视频| 人人妻人人看人人澡| 国产三级中文精品| 午夜福利免费观看在线| 亚洲成人久久爱视频| 伦理电影免费视频| 日韩中文字幕欧美一区二区| 亚洲精品粉嫩美女一区| 久久久久九九精品影院| 免费在线观看亚洲国产| 一级毛片精品| 欧美三级亚洲精品| 亚洲成人精品中文字幕电影| 精品国产美女av久久久久小说| 欧美性猛交黑人性爽| a级毛片在线看网站| 一级作爱视频免费观看| 欧美另类亚洲清纯唯美| 午夜日韩欧美国产| 亚洲熟女毛片儿| 国产精品av视频在线免费观看| 国产黄色小视频在线观看| 久久久久久亚洲精品国产蜜桃av| 一区二区三区国产精品乱码| 欧美日韩黄片免| 国产av一区二区精品久久| 听说在线观看完整版免费高清| 一进一出抽搐动态| 天堂av国产一区二区熟女人妻 | 一进一出好大好爽视频| 国产一区二区三区在线臀色熟女| 老司机午夜福利在线观看视频| 久热爱精品视频在线9| 在线免费观看的www视频| 男女视频在线观看网站免费 | 亚洲精品久久成人aⅴ小说| 亚洲精华国产精华精| 美女扒开内裤让男人捅视频| 女人爽到高潮嗷嗷叫在线视频| 午夜福利18| 亚洲无线在线观看| 久久久久国内视频| 熟女少妇亚洲综合色aaa.| 欧美日韩瑟瑟在线播放| 欧美黑人精品巨大| 俺也久久电影网| 久久久久精品国产欧美久久久| 成人特级黄色片久久久久久久| 精华霜和精华液先用哪个| 麻豆国产97在线/欧美 | 老司机午夜福利在线观看视频| 脱女人内裤的视频| 欧美乱色亚洲激情| 麻豆成人av在线观看| www.自偷自拍.com| 99在线人妻在线中文字幕| 男人的好看免费观看在线视频 | 又黄又粗又硬又大视频| 两性午夜刺激爽爽歪歪视频在线观看 | 非洲黑人性xxxx精品又粗又长| 成熟少妇高潮喷水视频| 夜夜躁狠狠躁天天躁| 特级一级黄色大片| 在线观看免费视频日本深夜| 久久久久久人人人人人| 99re在线观看精品视频| 国产高清视频在线观看网站| 香蕉丝袜av| 日韩大码丰满熟妇| 亚洲成a人片在线一区二区| 欧洲精品卡2卡3卡4卡5卡区| 国内精品久久久久精免费| 亚洲自偷自拍图片 自拍| 动漫黄色视频在线观看| 男人舔奶头视频| 国产蜜桃级精品一区二区三区| 亚洲国产欧美人成| 啪啪无遮挡十八禁网站| 一级作爱视频免费观看| 99热这里只有精品一区 | 日韩高清综合在线| 桃色一区二区三区在线观看| 国产又黄又爽又无遮挡在线| 国产亚洲精品综合一区在线观看 | 性色av乱码一区二区三区2| 给我免费播放毛片高清在线观看| 久久精品国产99精品国产亚洲性色| 成人18禁高潮啪啪吃奶动态图| 嫩草影视91久久| 2021天堂中文幕一二区在线观| 免费电影在线观看免费观看| 午夜亚洲福利在线播放| 午夜精品在线福利| 女人被狂操c到高潮| 亚洲人成网站高清观看| 日韩国内少妇激情av| 999久久久精品免费观看国产| 人人妻人人澡欧美一区二区| 美女高潮喷水抽搐中文字幕| 日韩欧美免费精品| 久久草成人影院| 亚洲一区中文字幕在线| 国产成人系列免费观看| 国产成人精品久久二区二区91| 性欧美人与动物交配| 国产精品乱码一区二三区的特点| 97人妻精品一区二区三区麻豆| 黄频高清免费视频| 国产精品亚洲av一区麻豆| 久久久久免费精品人妻一区二区| 亚洲成人中文字幕在线播放| 这个男人来自地球电影免费观看| 可以在线观看毛片的网站| 日韩精品青青久久久久久| 欧美性猛交黑人性爽| 夜夜夜夜夜久久久久| 亚洲最大成人中文| 亚洲精品久久成人aⅴ小说| 久9热在线精品视频| 成人精品一区二区免费| av在线播放免费不卡| 嫁个100分男人电影在线观看| 国产高清videossex| 三级男女做爰猛烈吃奶摸视频| 宅男免费午夜| 一本久久中文字幕| 男女那种视频在线观看| 午夜视频精品福利| 亚洲国产欧美人成| 国产成人av教育| 欧美乱码精品一区二区三区| av中文乱码字幕在线| 丰满的人妻完整版| 亚洲精品久久国产高清桃花| 国产av一区在线观看免费| 欧美精品亚洲一区二区| 日韩精品免费视频一区二区三区| 在线观看日韩欧美| 麻豆一二三区av精品| 成人av一区二区三区在线看| 亚洲av美国av| 很黄的视频免费| 欧美激情久久久久久爽电影| 久久精品国产清高在天天线| 在线看三级毛片| 久久精品91蜜桃| 狠狠狠狠99中文字幕| 亚洲欧洲精品一区二区精品久久久| 我的老师免费观看完整版| 亚洲av成人av| 亚洲人成电影免费在线| 又紧又爽又黄一区二区| 国产三级黄色录像| 国产成人啪精品午夜网站| 日韩精品青青久久久久久| 亚洲五月天丁香| 亚洲欧洲精品一区二区精品久久久| 手机成人av网站| 久久精品夜夜夜夜夜久久蜜豆 | 变态另类丝袜制服| 国产麻豆成人av免费视频| 久久精品国产综合久久久| 免费在线观看完整版高清| 中文字幕精品亚洲无线码一区| 两性午夜刺激爽爽歪歪视频在线观看 | 国产不卡一卡二| 岛国在线免费视频观看| 精品午夜福利视频在线观看一区| 日本熟妇午夜| 国产又色又爽无遮挡免费看| 男女视频在线观看网站免费 | av天堂在线播放| 亚洲欧美一区二区三区黑人| 日韩欧美在线二视频| 动漫黄色视频在线观看| 国产三级黄色录像| 成熟少妇高潮喷水视频| 日韩欧美 国产精品| 一夜夜www| 久久久久精品国产欧美久久久| 99久久久亚洲精品蜜臀av| 欧美一区二区国产精品久久精品 | 国产精品亚洲av一区麻豆| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美黑人欧美精品刺激| 成人av在线播放网站| 看片在线看免费视频| av天堂在线播放| 在线观看美女被高潮喷水网站 | 女生性感内裤真人,穿戴方法视频| 欧美色视频一区免费| 天堂av国产一区二区熟女人妻 | 精品国产乱子伦一区二区三区| 国产视频内射| 日日夜夜操网爽| 欧美三级亚洲精品| 久久天堂一区二区三区四区| 精品免费久久久久久久清纯| 久久久国产成人精品二区| 麻豆一二三区av精品| 国产精品久久久久久亚洲av鲁大| 中文字幕人妻丝袜一区二区| 欧美3d第一页| 欧美乱码精品一区二区三区| 日本五十路高清| 久99久视频精品免费| 午夜两性在线视频| 亚洲av第一区精品v没综合| 男男h啪啪无遮挡| 在线免费观看的www视频| √禁漫天堂资源中文www| 国产野战对白在线观看| 国产熟女午夜一区二区三区| 久久国产精品影院| 一个人免费在线观看电影 | 国产亚洲精品一区二区www| 老汉色av国产亚洲站长工具| svipshipincom国产片| 亚洲成人精品中文字幕电影| 亚洲人与动物交配视频| 中文亚洲av片在线观看爽| 美女黄网站色视频| 一进一出抽搐动态| 成人18禁高潮啪啪吃奶动态图| 男人舔女人的私密视频| 19禁男女啪啪无遮挡网站| 国产熟女午夜一区二区三区| 免费一级毛片在线播放高清视频| 麻豆久久精品国产亚洲av| 久久天堂一区二区三区四区| 久久精品成人免费网站| 人人妻人人看人人澡| 97碰自拍视频| 国产一区二区三区在线臀色熟女| 日本撒尿小便嘘嘘汇集6| 黑人操中国人逼视频| 中文在线观看免费www的网站 | 午夜免费观看网址| xxxwww97欧美| 国产亚洲精品一区二区www| 亚洲精品中文字幕在线视频| 亚洲熟妇熟女久久| 亚洲人与动物交配视频| 欧美绝顶高潮抽搐喷水| 18禁黄网站禁片免费观看直播| 日本一二三区视频观看| 日本一本二区三区精品| 在线观看www视频免费| 亚洲精品粉嫩美女一区| 午夜老司机福利片| 亚洲成人免费电影在线观看| 深夜精品福利| 麻豆国产av国片精品| 国产成人精品久久二区二区免费| 亚洲真实伦在线观看| 午夜福利视频1000在线观看| 黄色a级毛片大全视频| 老熟妇仑乱视频hdxx| av在线天堂中文字幕| 欧美丝袜亚洲另类 | 午夜福利欧美成人| 在线观看午夜福利视频| 久久国产精品人妻蜜桃| 国产男靠女视频免费网站| 亚洲av熟女| 成人18禁高潮啪啪吃奶动态图| 亚洲自拍偷在线| 亚洲七黄色美女视频| 久久伊人香网站| 国产一级毛片七仙女欲春2| 99在线人妻在线中文字幕| 在线播放国产精品三级| 美女午夜性视频免费| aaaaa片日本免费| 亚洲精品中文字幕在线视频| 一级片免费观看大全| 搡老熟女国产l中国老女人| 97人妻精品一区二区三区麻豆| 99久久国产精品久久久| 亚洲成人久久性| 午夜福利视频1000在线观看| 亚洲一区中文字幕在线| 一个人免费在线观看电影 | 国产成人精品久久二区二区91| 啦啦啦观看免费观看视频高清| 日韩 欧美 亚洲 中文字幕| 青草久久国产| 少妇熟女aⅴ在线视频| 亚洲熟妇熟女久久| 亚洲国产中文字幕在线视频| 午夜免费成人在线视频| 两性夫妻黄色片| 精品一区二区三区av网在线观看| 久久欧美精品欧美久久欧美| 嫁个100分男人电影在线观看| 成人高潮视频无遮挡免费网站| 久久久国产精品麻豆| 亚洲精品国产一区二区精华液| 亚洲熟女毛片儿| av视频在线观看入口| 久久人妻av系列| 老汉色∧v一级毛片| 亚洲成av人片免费观看| 在线观看免费视频日本深夜| 免费看十八禁软件| 不卡av一区二区三区| 少妇裸体淫交视频免费看高清 | 全区人妻精品视频| 亚洲五月天丁香| 午夜激情av网站| 日本 欧美在线| 亚洲午夜精品一区,二区,三区| 国产亚洲精品久久久久5区| 黑人操中国人逼视频| 欧美性猛交黑人性爽| 欧美高清成人免费视频www| 88av欧美| 亚洲男人的天堂狠狠| 欧美又色又爽又黄视频| 亚洲色图 男人天堂 中文字幕| 国产精品99久久99久久久不卡| 欧美性猛交黑人性爽| 国产又黄又爽又无遮挡在线| av国产免费在线观看| 狂野欧美白嫩少妇大欣赏| 在线a可以看的网站| 国产伦人伦偷精品视频| 黄色a级毛片大全视频| 亚洲精品久久成人aⅴ小说| 欧美色欧美亚洲另类二区| 三级国产精品欧美在线观看 | 亚洲人成伊人成综合网2020| 亚洲电影在线观看av| cao死你这个sao货| 欧美高清成人免费视频www| 中文字幕人妻丝袜一区二区| 蜜桃久久精品国产亚洲av| 99国产极品粉嫩在线观看| 亚洲欧美日韩高清专用| 久热爱精品视频在线9| 亚洲av第一区精品v没综合| 中文字幕av在线有码专区| 欧美黑人欧美精品刺激| 中文在线观看免费www的网站 | 桃色一区二区三区在线观看| 日本免费一区二区三区高清不卡| 桃红色精品国产亚洲av| 亚洲av五月六月丁香网| 亚洲精品国产精品久久久不卡| 99热这里只有是精品50| 国产91精品成人一区二区三区| svipshipincom国产片| 97人妻精品一区二区三区麻豆| 两人在一起打扑克的视频| 欧美又色又爽又黄视频| 热99re8久久精品国产| 好男人电影高清在线观看| 男插女下体视频免费在线播放| 亚洲色图av天堂| 欧美在线一区亚洲| 一本综合久久免费| 床上黄色一级片| 制服人妻中文乱码| 色在线成人网| 国产真实乱freesex| 麻豆久久精品国产亚洲av| 一本一本综合久久| 国产精品亚洲一级av第二区| 亚洲国产欧美一区二区综合| 亚洲欧美精品综合一区二区三区| 亚洲七黄色美女视频| 国产免费男女视频| 国产精品精品国产色婷婷| 亚洲国产欧美网| 国产精品精品国产色婷婷| 在线国产一区二区在线| 制服丝袜大香蕉在线| 欧美性长视频在线观看| а√天堂www在线а√下载| 老司机深夜福利视频在线观看| 最近视频中文字幕2019在线8| 久久久久免费精品人妻一区二区| 我要搜黄色片| 色老头精品视频在线观看| 这个男人来自地球电影免费观看| 日本成人三级电影网站| 51午夜福利影视在线观看| 成人高潮视频无遮挡免费网站| 国产亚洲精品久久久久5区| 熟妇人妻久久中文字幕3abv| 国产精品自产拍在线观看55亚洲| 看片在线看免费视频| 国产一区二区激情短视频| 黄色毛片三级朝国网站| a级毛片a级免费在线| 国产熟女xx| 国产精品精品国产色婷婷| 99久久99久久久精品蜜桃| 露出奶头的视频| 中文资源天堂在线| 亚洲人成77777在线视频| 亚洲国产日韩欧美精品在线观看 | 国产成人欧美在线观看| 亚洲国产精品999在线|