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

    聚氨酯反應擠出過程的數值模擬

    2016-08-11 06:07:18豪,原,
    關鍵詞:雙螺桿聚合度擠出機

    唐 豪, 宗 原, 趙 玲

    (華東理工大學化學工程聯合國家重點實驗室,上海 200237)

    ?

    聚氨酯反應擠出過程的數值模擬

    唐豪,宗原,趙玲

    (華東理工大學化學工程聯合國家重點實驗室,上海 200237)

    采用三維有限元數值模擬結合標量混合分析方法,研究了反應物進料狀態(tài)、轉速及進料單體物質的量之比對聚氨酯在雙螺桿擠出機中反應擠出過程的影響。反應進程不僅受停留時間控制,還與組分的均化程度有關。受組分濃度分布的影響,相比于預混進料,采用非預混進料方式導致反應程度較低;提高螺桿轉速雖可加速均化過程,但卻導致物料在擠出機中的停留時間縮短,因此需結合反應混合特性選取最優(yōu)轉速。通過數值模擬發(fā)現,采用90 r/min轉速能達到較高反應程度。另外,模擬計算了不同預聚程度進料對反應擠出過程的影響。結果表明,采用二次投料法可階段性控制聚合反應程度,提高反應擠出過程的穩(wěn)定性,同時也便于為不同反應階段的黏度范圍匹配對應的反應器。

    聚氨酯; 反應擠出; 混合; 數值模擬

    反應擠出過程通常以雙螺桿擠出機作為連續(xù)化反應器,可在擠出加工的同時實現聚合物的合成或改性。Ganzeveld等[1]將發(fā)生在擠出機中的反應分為兩類:單組分反應和多組分反應。前者的進料組分單一,進料流可認為處于微觀完全混合狀態(tài),例如對鏈式聚合反應擠出過程的三維模擬將引發(fā)劑和反應物的混合物視為均相流體,產品質量決定于物料的宏觀混合和溫度分布。對多組分聚合反應而言,組分相界面是聚合反應發(fā)生的主要場所,因而物料間的混合程度將影響最終聚合產品的質量[2]。當混合時間與反應時間相當或更長時,混合的影響將占主導地位[3-4]。因此,深入了解多組分聚合過程中的反應流特征對最終聚合產品質量的控制具有重要的指導意義。

    聚氨酯廣泛應用于發(fā)泡、彈性體、涂層等領域,其聚合過程是典型的多組分逐步聚合反應。Vargas等[5]將雙螺桿擠出機處理為多釜串聯模型,考察了聚氨酯反應擠出過程中轉速、產量和配料比對反應結果的影響。Puaux等[6]采用一維模型對聚氨酯的反應擠出過程進行建模,比較了完全分隔模型、微觀完全混合模型和平推流模型在不同轉速、單體物質的量之比下的反應程度。實際上,由于雙螺桿擠出機結構的復雜性,反應物料在螺槽中的流動具有復雜的三維特征,徑向及軸向流動之間存在復雜的相互影響;另一方面,在擠出過程中,物料的黏度、相對分子質量等也隨反應發(fā)生變化。采用簡單數學模型對反應擠出過程進行研究與實際過程存在偏差,無法真實描述其傳遞特征。因此,需建立更為準確的三維模型,綜合考慮流動和物性對反應擠出過程的影響。

    三維模型中評價擠出機混合能力的參數如混合指數、應力場、面積伸展率、混合效率等可有效描述擠出機的流場特征,但對組分濃度均化過程的表征作用有限。本質上說,混合操作的目的是減少各組分的不均勻性。由于反應擠出過程中不存在湍流和渦擴散,其決定性機理是本體擴散或對流混合[7]。Wu[8]于1994年提出了基于混合物層狀結構模型的三維模擬方法,其本質是進行濃度標量的混合分析,但限于當時的計算條件,該方法并未有效實施。隨著計算能力和理論的發(fā)展,基于計算流體力學的工作方法已被廣泛應用于各類反應器及混合-反應過程的研究中,如微通道反應器、反應注射成型等[9-11]。近年來,Fard等[12-13]開始采用標量混合分析方法對雙螺桿擠出機中的混合過程進行可視化分析,并利用分隔強度定量表征了擠出機的混合能力,但其研究中并未探討混合與反應的關系。

    本文利用三維有限元數值模擬結合標量混合分析方法考察聚氨酯的反應擠出過程,討論了進料條件、轉速和進料單體物質的量之比對組分均化的影響,并從反應工程的角度分析了這些因素對聚氨酯聚合反應擠出的作用。

    1 計算模型

    1.1反應體系模型

    根據Cassagnau[14]提出的反應體系,熱塑性聚氨酯以相對分子質量為1 000的聚環(huán)氧乙烷(PEO)為軟段,氫化苯基甲烷二異氰酸酯(H12MDI)為硬段,1,4-丁二醇(BDO)為擴鏈劑合成,BDO與PEO的物質的量之比為2∶1,合成過程[15]如下:

    (1)

    該反應符合二級反應動力學,反應速率(rA)由式(2)計算得到:

    rA=kcAcB

    (2)

    其中:cA為異氰酸基(A)濃度;cB為羥基(B)濃度。

    反應速率常數(k)

    k=10.9×109ccatexp(-48 500/RT)

    (3)

    其中:ccat為催化劑濃度;R為理想氣體常數;T為熱力學溫度。當催化劑濃度為0.000 475 mol/L,溫度為433 K時,k=0.146 L/(mol·s)。

    Cassagnau等[16]同時給出了該聚氨酯反應體系的黏度(η)隨剪切速率的變化規(guī)律,根據實驗參數,采用Carreau模型回歸了其本構方程

    (4)

    (5)

    聚合度采用相對分子質量矩的方法計算[6]

    (6)

    其中:Xn與XW分別是數均及重均聚合度;μi為i階相對分子質量矩。

    相對分子質量矩可由式(7)獲得

    (7)

    其中:cA0為異氰酸基初始濃度;cB0為羥基初始濃度;r=cA0/cB0是進料單體異氰酸基與羥基的物質的量之比;PA是組分A(異氰酸基)的轉化率。

    1.2數值建模

    流體設為不可壓縮流體,全充滿操作,計算為擬穩(wěn)態(tài)流動且滿足等溫條件,控制方程包括:

    連續(xù)性方程

    (8)

    動量方程

    (9)

    (10)

    其中:p為壓力;U為速度矢量;τ為應力張量;ρ為聚氨酸體系密度。

    組分濃度ci的輸送方程本質為層流擴散方程

    (11)

    Di為組分i的擴散系數,模擬中設Di=10-9m2/s。υi為組分i的化學計量數,對反應物,υi<0;對產物,υi>0。

    1.3數值條件

    采用的螺桿元件螺距為40 mm,長300 mm,前后均添加5 mm的發(fā)展段。雙螺桿擠出機的具體幾何約束條件見圖1,db為螺桿根部直徑,ds為螺桿頂部直徑,d為機筒直徑,CL為中心距。

    圖1 雙螺桿擠出機的幾何約束條件Fig.1 Constrains of the twin screw extruder

    螺桿和機筒壁面均采用無滑移邊界條件。忽略組分間界面張力作用,體系黏度隨反應進行而升高,可由式(4)~式(5)確定。不同組分反應物料分別從左、右入口進入擠出機,其流率均為3 300 mm3/s,因此,可認為兩種反應物在入口處完全分隔。隨著擠出過程的進行,反應物料在流場的作用下通過對流擴散進行均化。擠出機出口處設為大氣壓。采用基于有限元的數值方法對反應擠出過程進行數值模擬,螺桿元件為非結構網格填充,網格數約130 000;流體域為結構網格,同時在壁面處設置邊界層網格,網格數約175 000,見圖2。計算的收斂判據為殘差小于0.000 1,在八核工作站上每次計算大約耗時9 h。

    圖2 計算域的離散Fig.2 Discretization of simulated domain

    此外,分別于擠出機的非嚙合區(qū)和嚙合區(qū)作兩條取樣線P1-P1′及P2-P2′,以便考察物理量沿擠出方向的變化情況。

    1.4混合的表征

    停留時間分布及其方差在均相進料的反應器混合研究中得到廣泛應用,但對多組分體系混合的表征有限[17]。分隔強度(Iseg)反映了不同區(qū)域中組分濃度與平均值的偏差,更適合于描述多組分黏性流體的層流微觀混合特征[18]

    (12)

    1.5驗證

    停留時間分布曲線在一定程度上可反映擠出機內的流動情況,因此被用于驗證模擬的正確性。采用文獻[19]中的KD2元件進行模型驗證,其機筒內徑為36 mm,螺桿根徑和外徑分別為23.6 mm及35.6 mm,中心距30 mm,捏合盤厚度為7.5 mm,錯列角60°,見圖3(a)。操作條件為進料流率3.06×10-6m3/s,轉速為120 r/min,出口壓力3 MPa。圖3(b)比較了該邊界條件下文獻[19]與模擬所得的停留時間分布曲線,可發(fā)現兩者十分接近,相對誤差在5%以內,說明獲得的數值結果真實可靠。

    圖3 數值方法的驗證 (a) 元件的幾何結構 (b) 停留時間分布曲線Fig.3 Validation of simulation (a) elements geometry (b) RTD curves

    2 結果與討論

    2.1初始濃度條件的影響

    根據反應動力學方程,反應組分的初始濃度及局部濃度分布可影響聚合反應的局部反應速率,這是影響多組分反應擠出過程中不可忽視的因素。本節(jié)比較了相同轉速(30 r/min)下預混進料及非預混進料對聚氨酯聚合反應的影響。表1列出了兩種進料情況下的組分初始濃度分布。

    表1 預混與非預混條件下的組分進料初始濃度分布Table 1 Initial concentration distribution of premixed and non-premixed feedings

    圖4給出了不同進料條件下y=13 mm平面處聚合度(Xn)的分布情況。由圖4可見,兩種進料情況下反應聚合度均沿擠出方向不斷提高,但預混進料下的聚合程度明顯高于非預混進料。這是由于當采用預混進料時,兩種反應物已混合均勻,各區(qū)域組分濃度分布相同;相反,在非預混進料條件下,由于初始物料完全分隔,組分間接觸不充分導致局部反應速率下降。在此情況下,隨著擠出過程的進行組分間逐漸混合,濃度分布趨于均勻,見圖5。由此說明,當采用非預混進料時,需要消耗一定的反應空間及時間完成組分均化??傮w來說,反應物間不充分的混合抑制了反應速率,進而導致非預混條件下較低的聚合度。

    圖6示出了不同進料方式下組分A濃度沿直線P1-P1′和P2-P2′的變化過程,定義量綱為一的軸向長度為|z|/L。其中|z|為取樣點或取樣面的軸向位置,L為擠出機的總長。預混進料中,由于反應物的不斷消耗,曲線逐漸平滑衰減,且兩條曲線近乎重合,表明擠出機中的徑向濃度分布較為均勻。非預混進料中,兩條曲線在進口段均呈現強烈波動,經過一段擠出長度后,嚙合區(qū)及非嚙合區(qū)的濃度趨于接近,且濃度曲線呈現出平滑衰減的趨勢。這些波動本質上是均化過程的體現。

    圖4 不同進料下聚合度分布情況Fig.4 Distribution of Xn with various feeding conditions

    圖5 不同進料下cA分布情況Fig.5 Distribution of cA with various feeding conditions

    以上分析揭示了流體的微觀混合狀態(tài)對反應擠出結果的影響:在預混進料條件下,由于微觀混合程度已達到最高,反應只受停留時間的影響;相反,當反應物為非預混進料時,反應不僅與停留時間有關,也與物料的均化程度有關。兩種反應物之間混合發(fā)生的遲早對聚合反應的產品質量有較大影響,需盡可能縮短混合所需時間,提高均化效率。

    圖6 不同進料條件下cA沿P1-P1′和P2-P2′的變化Fig.6 Evolution of cA along P1-P1′ and P2-P2′ with various feeding conditions

    2.2轉速的影響

    圖7示出了不同轉速下分隔強度(Iseg)沿擠出機的變化,即采用分隔強度定量表征擠出機的微觀混合效果。從圖7可見,在相同進料流率下,高轉速導致分隔強度快速下降,說明高轉速加速了組分間的均化過程。在一定的軸向距離之后(即經歷一定混合后),分隔強度減小至接近于零,可認為組分間已接近微觀完全混合。

    圖7 不同轉速下分隔強度沿擠出機的變化Fig.7 Evolution of segregation intensity along the extruder at various rotating speeds

    圖8示出了非嚙合區(qū)內組分A濃度沿軸向的變化情況。由前所述,組分濃度的衰減與反應消耗有關,其波動反映了擠出機內的混合過程。由圖8可發(fā)現高轉速抑制了波動的幅度,同時也縮小了波動的軸向范圍,證明提高轉速有利于組分均化,對濃度效應具有正向作用。

    加快組分均化的最終目的是提高反應程度,為此,考察了在y=13 mm平面處反應聚合度的分布(圖9)。轉速在90 r/min下的數均聚合度的值總是高于30 r/min下的值,說明加速均化作用有助于反應程度的提高。然而當轉速提高到150 r/min時,在進口段聚合度與90 r/min下的值十分接近,但沿著擠出方向其聚合度與90 r/min下的值相比出現明顯下降。這可能與高轉速造成停留時間的減少有關[1]。

    圖8 不同轉速下cA沿P1-P1′的變化Fig.8 Evolution of cA along line P1-P1′ at various rotating speeds

    圖9 不同轉速下聚合度的分布情況Fig.9 Distribution of Xn at various rotating speeds

    圖10給出了不同轉速下沿擠出方向的局部平均停留時間,證明在給定的流率下,提高轉速可縮短停留時間,進而導致反應時間不足。同時也可以發(fā)現,由于停留時間的累積特性,沿著擠出方向其差異不斷增大。反應擠出過程與其他化學反應過程類似,多組分反應受組分分布和反應時間的共同影響,因而單純提高轉速對反應程度存在兩面性:一方面能加快組分間的均化,促進反應進行;另一方面則因減少停留時間而抑制反應的進一步發(fā)生。模擬結果表明,非預混進料時,反應初期由于組分的均化程度較低,混合在反應控制中占主導作用;在擠出機的混合作用下,均化程度沿軸向不斷提高,反應時間逐漸成為主要控制因素。為獲得較高的反應程度,需在確?;旌闲Ч疤嵯卤M可能延長停留時間。綜合組分分布和停留時間兩方面因素,本節(jié)中獲得的較優(yōu)轉速為90 r/min。

    圖10 不同轉速下沿擠出方向的平均停留時間Fig.10 Mean residence time along the extruder at various rotating speeds

    2.3單體物質的量之比的影響

    聚氨酯的線性聚合反應中,單體物質的量之比是一個至關重要的因素??刂茊误w物質的量之比進行多步加料有利于控制產物的相對分子質量[20]。為此,模擬計算了非預混進料條件下單體物質的量之比(r)分別為1、0.90和0.85情況下,雙螺桿擠出機內聚合度(XW)沿擠出方向的變化情況,如圖11所示。從中可見,在聚合初期,3條曲線在進口段比較接近,即單體物質的量之比對聚合度影響不大,其影響主要反映在高反應程度區(qū)??拷隹谔?由于單體物質的量之比不為1,聚合度受到抑制,單體物質的量之比為0.90和0.85時,出口處聚合度相比于r=1進料時的聚合度分別降低了9%和16%。

    圖11 不同單體物質的量之比下聚合度沿擠出方向的變化Fig.11 Evolution of XW along the extruder with various feed molar ratios

    根據式(5)可知,當聚合度為5時,反應體系零切黏度僅為8×10-3Pa·s;而當聚合度提升至100時,體系黏度約為120 Pa·s??梢?聚氨酯聚合過程中黏度變化范圍較大,需進行反應器選型以提高設備使用率。預聚階段體系黏度較低,聚合可在釜式反應器中進行;聚合反應中后期黏度顯著上升,則可利用雙螺桿擠出機對高黏流體的混煉輸送特性進一步提升產品相對分子質量。其中預聚階段的反應程度可通過單體物質的量之比控制。本文進一步計算了采用二次投料方法時,預聚時單體物質的量之比對雙螺桿擠出機內重均聚合度的影響,見圖12。模擬中擠出機長度為900 mm,總流率為3 300 mm3/s,假設反應物在擠出機投料前已經歷了預聚,預聚時單體物質的量之比分別為0.80,0.85及0.87,預聚結束后體系黏度分別達到0.28,0.72,1.4 Pa·s。二次投料時,為控制最終的單體物質的量之比為0.90,需在擠出機投料時補充相應配比的異氰酸基。從圖中可見,在入口段,聚合度增長速率最快,且預聚度越低,聚合度增長斜率越大。單體物質的量之比為(0.80+0.10)、(0.85+0.05)及(0.87+0.03)時對應斜率分別為68.1,53.7和41.5。沿擠出方向,聚合度的增長幅度逐漸放緩,在達到出口處時,聚合度基本達到穩(wěn)定狀態(tài)??梢娡ㄟ^調節(jié)預聚度可控制擠出機中相對分子質量的增長速率,實現較高相對分子質量產品的穩(wěn)定擠出。

    圖12 不同單體物質的量之比下聚合度 沿擠出方向的變化(二次投料)Fig.12 Evolution of XW along the extruder with various feed molar ratios (two-step feeding)

    在雙螺桿擠出機中,物料傳送主要依靠螺槽內的壓力流及拖曳流,而這兩種流動均與物料的黏度成正比[21]。較低的黏度會減弱雙螺桿擠出機的拖曳性能;低預聚度也會引起物系在擠出過程中黏度存在較大的波動,不利于對反應擠出過程的控制。計算結果表明,采用二次投料法能實現相對分子質量的階段性增長,達到控制反應體系物性的目的。這對某些熱敏性或高反應活性單體的聚合產物意義尤為突出,且這種半連續(xù)化的生產方法將有效提高雙螺桿擠出機的使用效率和聚合產品質量。

    3 結 論

    聚氨酯的合成是典型的多組分聚合反應,其微觀混合情況是決定產品質量的重要因素。本文利用三維CFD模擬,并結合標量混合分析方法,考察了多組分在雙螺桿擠出機內的輸運-反應情況。

    與預混進料不同,非預混進料條件下擠出機中組分的均化行為表現為其濃度沿擠出方向波動衰減。從統(tǒng)計學的角度出發(fā),均化過程的實質為分隔強度的不斷減小。同時,受濃度效應的影響,非預混進料的反應程度低于預混進料的反應程度。另外,多組分聚合反應中為獲得較高的反應聚合度,需同時考慮混合程度和反應時間。較高轉速可以實現組分的快速均化,但也會縮短停留時間,從而影響反應程度。通過對不同轉速的考察發(fā)現,螺桿轉速為90 r/min時能獲得較高的聚合程度。

    單體物質的量之比對聚合度有顯著影響。非等量進料會抑制高反應程度區(qū)相對分子質量的進一步上升。利用該特性控制預聚程度可實現反應體系黏度的階段性控制。模擬結果表明,二次投料法中,體系相對分子質量增長速率對預聚度較為敏感,單體物質的量之比為(0.80+0.10)時,預聚度為15,相對分子質量相對增長速率為68.1;而單體物質的量之比為(0.87+0.03)時,預聚度為25,相對增長速率為41.5,說明調節(jié)預聚程度可有效控制擠出機中相對分子質量的增長速率,提高反應擠出過程的穩(wěn)定性。同時,二次進料方式也避免了為輸送低黏物料而進行的螺桿設計,有效提升了設備的使用效率。

    [1]GANZEVELD K J,JANSSEN L P B M.A mixing model for multicomponent reactions in twin screw extruders applied to the polymerization of urethanes[J].Polymer Engineering & Science,1992,32(7):457-466.

    [2]GANZEVELD K J,JANSSEN L P B M.Role of mixing and rheology in reactive extrusion[J].Industrial & Engineering Chemistry Research,1994,33(10):2398-2403.

    [3]SURESH A,CHAKRABORTY S,KARGUPTA K,etal.Low-dimensional models for describing mixing effects in reactive extrusion of polypropylene[J].Chemical Engineering Science,2008,63(14):3788-3801.

    [5]VARGAS R O,VIVALDO-LIMA E,MANERO O.Simulation of nonlinear polyurethane production in a twin-screw extruder[J].Polymer-Plastics Technology and Engineering,2006,45(1):9-21.

    [6]PUAUX J P,CASSAGNAU P,BOZGA G,etal.Modeling of polyurethane synthesis by reactive extrusion[J].Chemical Engineering and Processing,2006,45(6):481-487.

    [7]LAWAL A,KALYON D M.Mechanisms of mixing in single and co-rotating twin screw extruders[J].Polymer Engineering & Science,1995,35(17):1325-1338.

    [8]WU Juichuang.Micromixing of a polymer melt in twin screw extruders with kneading disks[D].Newark:University of Delaware,1994.

    [9]WANG Wentan,ZHAO Shufang,SHAO Ting,etal.Numerical study of mixing behavior with chemical reactions in micro-channels by a lattice boltzmann method[J].Chemical Engineering Science,2012,84(52):148-154.

    [10]SANTOS R J,TEIXEIRA A M,LOPES J C B.Study of mixing and chemical reaction in RIM[J].Chemical Engineering Science,2005,60(8-9):2381-2398.

    [11]AOKI N,HASEBE S,MAE K.Mixing in microreactors:Effectiveness of lamination segments as a form of feed on product distribution for multiple reactions[J].Chemical Engineering Journal,2004,101(1-3):323-331.

    [12]FARD A S,ANDERSON P D.Simulation of distributive mixing inside mixing elements of co-rotating twin-screw extruders[J].Computers & Fluids,2013,87:79-91.

    [13]FARD A S,HULSEN M A,MEIJER H E H,etal.Tools to simulate distributive mixing in twin-screw extruders[J].Macromolecular Theory and Simulations,2012,21(4):217-240.

    [14]CASSAGNAU P,NIETSCH T,BERT M,etal.Reactive blending by in situ polymerization of the dispersed phase[J].Polymer,1999,40(1):131-138.

    [15]SEMSARZADEH M A,NAVARCHIAN A H,MORSHEDIAN J.Reactive extrusion of poly(urethane-isocyanurate)[J].Advances in Polymer Technology,2004,23(3):239-255.

    [16]CASSAGNAU P,NIETSCH T,MICHEL A.Bulk and dispersed phase polymerization of urethane in twin screw extruders[J].International Polymer Processing,1999,14(2):144-151.

    [17]OTTINO J M,CHELLA R.Laminar mixing of polymeric liquids;a brief review and recent theoretical developments[J].Polymer Engineering & Science,1983,23(7):357-379.

    [18]TADMOR Z,GOGOS C G.Principles of polymer processing[M].2nd Ed.Hoboken:John Wiley & Sons,2006:385-389.

    [19]ZHANG Xianming,FENG Lianfang,CHEN Wenxing,etal.Numerical simulation and experimental validation of mixing performance of kneading discs in a twin screw extruder[J].Polymer Engineering & Science,2009,49(9):1772-1783.

    [20]姜濤,陸興軍,趙明,等.分步加料法合成聚酯多元醇[J].聚氨酯工業(yè),1997,12(3):11-14.

    [21]ZHU Linjie,NARH K A,HYUN K S.Investigation of mixing mechanisms and energy balance in reactive extrusion using three-dimensional numerical simulation method[J].International Journal of Heat and Mass Transfer,2005,48(16):3411-3422.

    Numerical Simulation of Reactive Extrusion of Polyurethane

    TANG Hao,ZONG Yuan,ZHAO Ling

    (State Key Laboratory of Chemical Engineering,East China University of Science and Technology,Shanghai 200237,China)

    Three dimensional finite element numerical simulation coupled with scalar mixing analysis was employed to investigate the effects of feeding condition,rotating speed and feed monomer ratio on reactive extrusion of polyurethane.The reaction course is controlled by the residence time and the degree of homogenization.The polymerization degree with non-premixed feeding is lower than that of permixed feeding.Higher rotating speed is beneficial to homogenization,while it also leads to the reduction of residence time.Therefore,the rotating speed should be optimized according to the mixing-reaction characteristics.Results show that higher polymerization degree in current study is obtained under 90 r/min.Moreover,the effect of prepolymer feeding on reactive extrusion is also simulated.Results indicate that the polymerization degree can be controlled via two-step feeding in order to improve the stability of reactive extrusion and to match the appropriate reactors for various reaction periods.

    polyurethane; reactive extrusion; mixing; numerical simulation

    A

    1006-3080(2016)03-0301-07

    10.14135/j.cnki.1006-3080.2016.03.002

    2015-11-06

    國家重點基礎研究發(fā)展(973)計劃(2011CB606100);國家自然科學基金(21406059)

    唐豪(1987-),男,浙江嘉興人,博士生,研究方向為聚合物反應擠出的數值模擬。E-mail:rain520737420@126.com

    通信聯系人:趙玲,E-mail:zhaoling@ecust.edu.cn

    TQ021.1

    猜你喜歡
    雙螺桿聚合度擠出機
    膜分離工藝提高產品中高聚合度ε-聚賴氨酸含量
    雙螺桿擠壓對沙米復合粉理化及糊化特性的影響
    耦合雙槽單螺桿擠出機停留時間研究
    中國塑料(2016年7期)2016-04-16 05:25:54
    聚醋酸乙烯聚合度的在線監(jiān)測
    安徽化工(2016年5期)2016-02-27 08:25:04
    三螺桿擠出機熔體輸送能力的評價
    中國塑料(2015年2期)2015-10-14 05:34:32
    聚乳酸用單螺桿擠出機熔體輸送段數值模擬研究
    中國塑料(2015年11期)2015-10-14 01:14:22
    自潔型非對稱同向雙螺桿擠出機混合分析
    中國塑料(2014年8期)2014-10-17 03:05:10
    聚合度7-15的殼寡糖抑制脂多糖刺激的單核細胞產生TNF-α和IL-8的作用研究
    雙螺桿擠出炸藥模擬物的流場模擬分析及實驗研究
    火炸藥學報(2014年5期)2014-03-20 13:17:54
    雙螺桿食品擠壓流動過程數值模擬
    亚洲国产欧洲综合997久久,| 精品久久久久久久久久久久久| 免费看日本二区| 亚洲av电影在线观看一区二区三区 | 男人狂女人下面高潮的视频| 欧美一区二区亚洲| 久久热精品热| 国产探花极品一区二区| av在线老鸭窝| 亚洲欧美一区二区三区国产| 国产成人精品久久久久久| 高清毛片免费看| 国产黄片视频在线免费观看| 免费观看a级毛片全部| 久久久久久久久久黄片| 国产伦理片在线播放av一区| 日韩在线高清观看一区二区三区| 日韩在线高清观看一区二区三区| 日日摸夜夜添夜夜爱| 少妇裸体淫交视频免费看高清| www.色视频.com| 国产真实伦视频高清在线观看| 免费av毛片视频| 久久久久久国产a免费观看| 国产黄色小视频在线观看| 大香蕉97超碰在线| 国产精品久久久久久久久免| 国产v大片淫在线免费观看| 精品人妻一区二区三区麻豆| 在线播放无遮挡| 99视频精品全部免费 在线| 久久99精品国语久久久| 亚洲一级一片aⅴ在线观看| 97在线视频观看| 嫩草影院入口| 久99久视频精品免费| 女的被弄到高潮叫床怎么办| 国产欧美日韩精品一区二区| h日本视频在线播放| 亚洲精品aⅴ在线观看| 婷婷色综合大香蕉| 亚洲欧洲日产国产| 亚洲国产最新在线播放| 色视频www国产| 十八禁国产超污无遮挡网站| 久久精品国产亚洲av涩爱| 国产免费一级a男人的天堂| 在线观看一区二区三区| 久久久久久大精品| 插阴视频在线观看视频| 欧美3d第一页| 亚洲图色成人| 韩国高清视频一区二区三区| 精品欧美国产一区二区三| 男人的好看免费观看在线视频| 亚洲国产精品成人久久小说| 成人性生交大片免费视频hd| 欧美成人一区二区免费高清观看| av视频在线观看入口| 亚洲欧美日韩卡通动漫| 亚洲精品影视一区二区三区av| av卡一久久| 欧美色视频一区免费| 精品无人区乱码1区二区| 国产69精品久久久久777片| 亚洲国产成人一精品久久久| 一级黄片播放器| 性插视频无遮挡在线免费观看| 最近中文字幕高清免费大全6| 狂野欧美激情性xxxx在线观看| 国产男人的电影天堂91| 色噜噜av男人的天堂激情| 久久久午夜欧美精品| 中文欧美无线码| 亚洲国产精品专区欧美| 欧美又色又爽又黄视频| 黄色一级大片看看| 精品久久久噜噜| 亚洲美女搞黄在线观看| 国产精品电影一区二区三区| videos熟女内射| 国产精品麻豆人妻色哟哟久久 | 亚洲一级一片aⅴ在线观看| 97超碰精品成人国产| 亚洲一区高清亚洲精品| 成人性生交大片免费视频hd| av国产久精品久网站免费入址| 免费看美女性在线毛片视频| 少妇人妻一区二区三区视频| ponron亚洲| 女的被弄到高潮叫床怎么办| 毛片一级片免费看久久久久| 色吧在线观看| 精华霜和精华液先用哪个| 你懂的网址亚洲精品在线观看 | 久久久久久伊人网av| 少妇熟女aⅴ在线视频| 可以在线观看毛片的网站| 五月伊人婷婷丁香| 国产白丝娇喘喷水9色精品| 99久久成人亚洲精品观看| 女的被弄到高潮叫床怎么办| 毛片女人毛片| 三级国产精品欧美在线观看| 老女人水多毛片| 国模一区二区三区四区视频| 亚洲人成网站高清观看| 欧美日本亚洲视频在线播放| 在线观看66精品国产| 一区二区三区高清视频在线| 蜜臀久久99精品久久宅男| 国产精品一区二区性色av| 亚洲国产高清在线一区二区三| 寂寞人妻少妇视频99o| 两个人的视频大全免费| 国产精品综合久久久久久久免费| 在线观看66精品国产| 夜夜看夜夜爽夜夜摸| 爱豆传媒免费全集在线观看| 久久精品国产亚洲网站| 免费无遮挡裸体视频| 久久久久精品久久久久真实原创| 1024手机看黄色片| 欧美日韩国产亚洲二区| 精华霜和精华液先用哪个| 免费搜索国产男女视频| 成人美女网站在线观看视频| 黄片无遮挡物在线观看| 亚洲婷婷狠狠爱综合网| 一级爰片在线观看| 成人毛片60女人毛片免费| 有码 亚洲区| 麻豆乱淫一区二区| 国产白丝娇喘喷水9色精品| 我的女老师完整版在线观看| 国产黄色视频一区二区在线观看 | 九色成人免费人妻av| 麻豆乱淫一区二区| 亚洲色图av天堂| 欧美bdsm另类| 日本-黄色视频高清免费观看| 国产精品熟女久久久久浪| 欧美成人免费av一区二区三区| 亚洲人成网站在线观看播放| 国产成人精品久久久久久| 亚洲欧洲国产日韩| 国产男人的电影天堂91| 亚洲国产最新在线播放| 亚洲欧美精品自产自拍| 国产黄a三级三级三级人| 国产精品一区二区在线观看99 | 日日摸夜夜添夜夜爱| 麻豆av噜噜一区二区三区| 欧美潮喷喷水| 亚洲电影在线观看av| av播播在线观看一区| 日韩强制内射视频| 丝袜美腿在线中文| 大香蕉97超碰在线| 欧美xxxx性猛交bbbb| 国产激情偷乱视频一区二区| 亚洲欧洲国产日韩| 亚洲四区av| 人人妻人人看人人澡| 五月伊人婷婷丁香| 女人十人毛片免费观看3o分钟| 午夜精品国产一区二区电影 | 男人和女人高潮做爰伦理| 欧美成人精品欧美一级黄| 美女高潮的动态| 伦理电影大哥的女人| 国产 一区精品| 18禁裸乳无遮挡免费网站照片| 亚洲丝袜综合中文字幕| 成人国产麻豆网| 麻豆成人午夜福利视频| 看免费成人av毛片| 亚洲欧洲日产国产| 亚洲国产精品久久男人天堂| 精品一区二区三区人妻视频| 一夜夜www| videos熟女内射| 波多野结衣高清无吗| 91精品伊人久久大香线蕉| 久久久精品大字幕| 亚洲欧美精品专区久久| 亚洲av免费高清在线观看| 久久亚洲国产成人精品v| 精品久久久久久电影网 | 国产免费一级a男人的天堂| 淫秽高清视频在线观看| 国产 一区 欧美 日韩| 国产精品国产三级专区第一集| 18禁在线无遮挡免费观看视频| 国产精华一区二区三区| 亚洲国产欧美人成| 亚洲av一区综合| 国产av不卡久久| 丝袜美腿在线中文| 国产色爽女视频免费观看| 国内精品美女久久久久久| 精品久久久久久电影网 | 最近视频中文字幕2019在线8| 97热精品久久久久久| 亚洲美女搞黄在线观看| www.av在线官网国产| 久久久久精品久久久久真实原创| 99久久无色码亚洲精品果冻| 色尼玛亚洲综合影院| 黄色一级大片看看| 最新中文字幕久久久久| 日日干狠狠操夜夜爽| 特大巨黑吊av在线直播| 91狼人影院| 国产av在哪里看| 日本av手机在线免费观看| 韩国高清视频一区二区三区| 91午夜精品亚洲一区二区三区| 99国产精品一区二区蜜桃av| 日本wwww免费看| 毛片一级片免费看久久久久| 韩国av在线不卡| 免费观看精品视频网站| 白带黄色成豆腐渣| 欧美bdsm另类| 国产av不卡久久| 一区二区三区四区激情视频| 高清在线视频一区二区三区 | 欧美另类亚洲清纯唯美| 一级毛片aaaaaa免费看小| 婷婷色麻豆天堂久久 | 尾随美女入室| av在线观看视频网站免费| 观看美女的网站| 成人二区视频| av福利片在线观看| 国产亚洲91精品色在线| 韩国高清视频一区二区三区| 国产精品乱码一区二三区的特点| 亚洲国产精品专区欧美| 日韩高清综合在线| av专区在线播放| 97人妻精品一区二区三区麻豆| 国产高清国产精品国产三级 | 国产精品野战在线观看| 亚洲精品久久久久久婷婷小说 | av线在线观看网站| 小说图片视频综合网站| 免费播放大片免费观看视频在线观看 | 丰满人妻一区二区三区视频av| av免费在线看不卡| 综合色av麻豆| 国产精品一区www在线观看| 欧美+日韩+精品| 国产三级中文精品| 中文精品一卡2卡3卡4更新| 国产精品久久视频播放| 中文字幕久久专区| 午夜福利网站1000一区二区三区| 97在线视频观看| 国产精品美女特级片免费视频播放器| 国产成人午夜福利电影在线观看| 久久热精品热| 亚洲国产精品专区欧美| 免费黄色在线免费观看| 亚洲在久久综合| 国产精品人妻久久久影院| 最近2019中文字幕mv第一页| 国产精品人妻久久久久久| 亚洲欧洲国产日韩| 狂野欧美激情性xxxx在线观看| 午夜福利视频1000在线观看| 久久人妻av系列| 小说图片视频综合网站| 国产成人精品婷婷| 非洲黑人性xxxx精品又粗又长| 一个人看视频在线观看www免费| 精品国产三级普通话版| 日韩欧美精品免费久久| 久久精品熟女亚洲av麻豆精品 | 免费看日本二区| 麻豆成人av视频| 亚洲av成人av| 亚洲人成网站高清观看| 色吧在线观看| 亚洲美女搞黄在线观看| 色噜噜av男人的天堂激情| 韩国av在线不卡| 插逼视频在线观看| 高清日韩中文字幕在线| 国产精品熟女久久久久浪| 日韩,欧美,国产一区二区三区 | 高清视频免费观看一区二区 | 国产美女午夜福利| 一区二区三区高清视频在线| 97热精品久久久久久| 少妇熟女aⅴ在线视频| 亚洲最大成人av| 欧美一级a爱片免费观看看| 精品国产三级普通话版| 欧美成人精品欧美一级黄| 身体一侧抽搐| 日韩亚洲欧美综合| 寂寞人妻少妇视频99o| 国产精品久久久久久av不卡| 黄色配什么色好看| 嫩草影院新地址| 国产又色又爽无遮挡免| 欧美一区二区国产精品久久精品| 久久这里有精品视频免费| 亚州av有码| 久久精品国产鲁丝片午夜精品| 精品久久久久久久久亚洲| 秋霞在线观看毛片| 日韩av在线免费看完整版不卡| 精品欧美国产一区二区三| 亚洲在久久综合| 青春草亚洲视频在线观看| 99热精品在线国产| 中文精品一卡2卡3卡4更新| 国产在视频线在精品| 亚洲欧洲国产日韩| 色尼玛亚洲综合影院| 少妇丰满av| 嘟嘟电影网在线观看| 蜜臀久久99精品久久宅男| 水蜜桃什么品种好| 日韩av不卡免费在线播放| 国产片特级美女逼逼视频| 久久精品91蜜桃| 亚洲高清免费不卡视频| 中文字幕av在线有码专区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲丝袜综合中文字幕| 村上凉子中文字幕在线| 日韩高清综合在线| 亚洲人成网站在线观看播放| 热99在线观看视频| 简卡轻食公司| 伦理电影大哥的女人| 国产男人的电影天堂91| 如何舔出高潮| 亚洲精品成人久久久久久| 国产亚洲一区二区精品| 国产真实伦视频高清在线观看| a级毛色黄片| 亚洲精品亚洲一区二区| 青春草国产在线视频| 日韩av在线免费看完整版不卡| av女优亚洲男人天堂| 精品一区二区三区人妻视频| 亚洲欧美中文字幕日韩二区| ponron亚洲| 日日摸夜夜添夜夜添av毛片| 久久久精品欧美日韩精品| 久久久久久久久中文| 久久久久久国产a免费观看| 可以在线观看毛片的网站| 99在线人妻在线中文字幕| 啦啦啦观看免费观看视频高清| 国产高清视频在线观看网站| 国产精品电影一区二区三区| 日本猛色少妇xxxxx猛交久久| 国产高清三级在线| 99在线人妻在线中文字幕| 大又大粗又爽又黄少妇毛片口| 中文精品一卡2卡3卡4更新| 久久久精品大字幕| 在线免费观看的www视频| 久久精品国产鲁丝片午夜精品| 日本wwww免费看| 成年av动漫网址| 日韩中字成人| 久久久久久久久久久丰满| 久久精品夜夜夜夜夜久久蜜豆| 免费黄网站久久成人精品| 3wmmmm亚洲av在线观看| 极品教师在线视频| 男女下面进入的视频免费午夜| 中文字幕制服av| 婷婷色综合大香蕉| 欧美另类亚洲清纯唯美| 亚洲伊人久久精品综合 | 精品一区二区三区人妻视频| 国产人妻一区二区三区在| 两个人的视频大全免费| 午夜激情福利司机影院| 日产精品乱码卡一卡2卡三| 日韩一区二区三区影片| 久久久午夜欧美精品| 非洲黑人性xxxx精品又粗又长| 日产精品乱码卡一卡2卡三| 国产麻豆成人av免费视频| 亚洲精品aⅴ在线观看| 精品人妻熟女av久视频| 成人毛片60女人毛片免费| 国产真实乱freesex| 日韩人妻高清精品专区| 久久精品久久久久久久性| 国产精品福利在线免费观看| 如何舔出高潮| 国产成人91sexporn| 老师上课跳d突然被开到最大视频| 亚洲av.av天堂| 国产一区有黄有色的免费视频 | 国产视频首页在线观看| 国产探花极品一区二区| 丝袜美腿在线中文| 老司机影院毛片| 男人的好看免费观看在线视频| 亚洲色图av天堂| av在线播放精品| 国产精品国产高清国产av| 嫩草影院入口| 天堂中文最新版在线下载 | 草草在线视频免费看| 99久久人妻综合| 亚洲av成人精品一二三区| 精品人妻一区二区三区麻豆| 在线免费观看的www视频| 精品国内亚洲2022精品成人| 日韩强制内射视频| 看免费成人av毛片| 国产在视频线精品| 国产爱豆传媒在线观看| a级毛色黄片| av免费观看日本| 日本三级黄在线观看| 日韩大片免费观看网站 | 日本午夜av视频| 最后的刺客免费高清国语| 免费观看在线日韩| 亚洲四区av| 久久精品人妻少妇| 国产一级毛片在线| 亚洲av.av天堂| 联通29元200g的流量卡| 国产精品女同一区二区软件| 欧美激情久久久久久爽电影| 性色avwww在线观看| 亚洲中文字幕一区二区三区有码在线看| 午夜免费激情av| 国产免费福利视频在线观看| 亚洲欧美日韩高清专用| 午夜福利在线在线| 一边摸一边抽搐一进一小说| 伦精品一区二区三区| 黑人高潮一二区| 99久久中文字幕三级久久日本| 最近手机中文字幕大全| 网址你懂的国产日韩在线| 男女视频在线观看网站免费| 亚洲在线观看片| 国产爱豆传媒在线观看| 51国产日韩欧美| 国产精品一区二区性色av| 欧美高清性xxxxhd video| 精品国产一区二区三区久久久樱花 | 国产精品乱码一区二三区的特点| 国产黄色小视频在线观看| 日韩,欧美,国产一区二区三区 | 久久精品国产亚洲av涩爱| 国产av不卡久久| 久久久a久久爽久久v久久| www日本黄色视频网| 色视频www国产| 久久久亚洲精品成人影院| 国产精品嫩草影院av在线观看| 日本黄色视频三级网站网址| 精品久久久久久久末码| 少妇熟女aⅴ在线视频| 校园人妻丝袜中文字幕| 色5月婷婷丁香| 国国产精品蜜臀av免费| 成人美女网站在线观看视频| 最近最新中文字幕免费大全7| 国产久久久一区二区三区| 亚洲精品日韩av片在线观看| 亚洲一区高清亚洲精品| 国产精品国产三级专区第一集| 国产成人aa在线观看| 国产高潮美女av| 午夜精品在线福利| 在线天堂最新版资源| 亚洲精品影视一区二区三区av| 伦理电影大哥的女人| 欧美激情在线99| 久久久久国产网址| 国产成人福利小说| 国产色婷婷99| 国产av不卡久久| 成年女人永久免费观看视频| 精品无人区乱码1区二区| 纵有疾风起免费观看全集完整版 | 日韩精品青青久久久久久| 精品人妻熟女av久视频| 国内少妇人妻偷人精品xxx网站| av女优亚洲男人天堂| 1000部很黄的大片| 好男人在线观看高清免费视频| 免费观看在线日韩| 2021少妇久久久久久久久久久| 亚洲成色77777| 男女国产视频网站| 十八禁国产超污无遮挡网站| 美女黄网站色视频| 超碰97精品在线观看| 亚洲欧美精品自产自拍| 男人的好看免费观看在线视频| 全区人妻精品视频| 99久久精品热视频| 国产精品永久免费网站| 国产一区二区在线av高清观看| 不卡视频在线观看欧美| 大香蕉97超碰在线| 欧美日韩在线观看h| 日韩欧美国产在线观看| 亚洲av电影在线观看一区二区三区 | 建设人人有责人人尽责人人享有的 | 一个人观看的视频www高清免费观看| 午夜福利在线观看免费完整高清在| 乱系列少妇在线播放| 在线观看一区二区三区| 99热精品在线国产| 亚洲欧美日韩卡通动漫| 欧美一区二区精品小视频在线| 亚洲av电影不卡..在线观看| 日本猛色少妇xxxxx猛交久久| 久久久久国产网址| 午夜爱爱视频在线播放| 亚洲欧洲国产日韩| 国产午夜精品一二区理论片| videossex国产| 午夜福利在线在线| 国产大屁股一区二区在线视频| 国产精品精品国产色婷婷| 老司机影院毛片| 久久人人爽人人片av| 亚洲综合色惰| 中文资源天堂在线| 亚洲欧美清纯卡通| 少妇人妻一区二区三区视频| 啦啦啦观看免费观看视频高清| av在线播放精品| 日韩欧美 国产精品| 精品久久久久久成人av| 国产亚洲午夜精品一区二区久久 | 国产精品99久久久久久久久| 美女cb高潮喷水在线观看| 亚洲成人精品中文字幕电影| 国产高清三级在线| 日本爱情动作片www.在线观看| 成人毛片60女人毛片免费| 成人亚洲精品av一区二区| 亚洲国产欧洲综合997久久,| 色5月婷婷丁香| av视频在线观看入口| 国产精品一区二区三区四区免费观看| 99国产精品一区二区蜜桃av| 国产精品1区2区在线观看.| 午夜激情福利司机影院| av免费在线看不卡| 亚洲在久久综合| 亚洲国产色片| 男女下面进入的视频免费午夜| 亚洲av中文字字幕乱码综合| 精品久久久久久久末码| 国产三级在线视频| 午夜精品一区二区三区免费看| 亚洲欧美日韩卡通动漫| 男女啪啪激烈高潮av片| 级片在线观看| 观看美女的网站| 永久免费av网站大全| 高清毛片免费看| 乱人视频在线观看| av天堂中文字幕网| 春色校园在线视频观看| 亚洲激情五月婷婷啪啪| 嘟嘟电影网在线观看| 亚洲伊人久久精品综合 | 亚洲国产精品专区欧美| 伦精品一区二区三区| 国产淫片久久久久久久久| 免费播放大片免费观看视频在线观看 | 男的添女的下面高潮视频| 国产91av在线免费观看| 国产视频首页在线观看| 成年女人看的毛片在线观看| 老司机福利观看| 淫秽高清视频在线观看| 亚洲第一区二区三区不卡| 久久久久久久国产电影| 激情 狠狠 欧美| 亚洲中文字幕一区二区三区有码在线看| 中文字幕亚洲精品专区| 久久精品综合一区二区三区| 国产熟女欧美一区二区| 91精品国产九色| 欧美性猛交黑人性爽| 亚洲国产高清在线一区二区三| 中文字幕亚洲精品专区| 日韩强制内射视频| 寂寞人妻少妇视频99o| 内地一区二区视频在线| 一级二级三级毛片免费看| 大话2 男鬼变身卡| 亚洲不卡免费看| 久久久精品欧美日韩精品| 国产精品一区二区性色av| 国产在线一区二区三区精 | 99久久精品一区二区三区| 一区二区三区免费毛片| 精品人妻视频免费看| 赤兔流量卡办理|