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

    陽離子捕收劑對高嶺石的捕收性能及動力學(xué)模擬

    2017-05-16 11:39:52郭麗娜李志紅朱張磊劉彥麗樊民強
    中國礦業(yè) 2017年5期
    關(guān)鍵詞:高嶺石收劑電位

    郭麗娜,李志紅,朱張磊,劉彥麗,樊民強

    (太原理工大學(xué)礦業(yè)工程學(xué)院,山西 太原 030024)

    陽離子捕收劑對高嶺石的捕收性能及動力學(xué)模擬

    郭麗娜,李志紅,朱張磊,劉彥麗,樊民強

    (太原理工大學(xué)礦業(yè)工程學(xué)院,山西 太原 030024)

    通過單礦物浮選試驗、Zeta動電位測定、紅外光譜測試和MS 6.0分子動力學(xué)模擬,分析了十二胺鹽酸鹽(DAH)和十二胺聚氧乙烯醚(AC1201)對煤系高嶺石的浮選行為和作用機理。結(jié)果表明,藥劑濃度為750 g/t時,AC1201可較好地浮選高嶺石。借助MS軟件模擬,構(gòu)建DAH和AC1201單分子在高嶺石(001)面的最優(yōu)吸附構(gòu)型,計算結(jié)果表明:AC1201分子距高嶺石(001)面的最小距離小于DAH分子距高嶺石表面的最小距離,而且AC1201分子與礦物表面的接觸面積大于DAH,吸附能更低。理論計算和浮選實驗結(jié)果相互一致,可為探索高嶺石浮選藥劑提供理論基礎(chǔ)。

    煤系高嶺石;浮選;十二胺聚氧乙烯醚;動力學(xué)模擬

    煤系高嶺巖是一種寶貴的自然資源,儲量豐富,開發(fā)后可用于建筑、印染、塑料、化工電子、農(nóng)業(yè)及耐火材料等行業(yè)[1]。然而我國大部分煤系高嶺巖被視為尾礦直接廢棄,數(shù)量巨大,既污染環(huán)境又浪費資源[2]。山西省大同塔山煤礦、孝義、忻州等地的煤矸石品位好,屬于煤系優(yōu)質(zhì)高嶺巖,開發(fā)前景極好[3],對其進行回收利用對工業(yè)發(fā)展具有很強的指導(dǎo)意義。

    工業(yè)中常用浮選法生產(chǎn)高嶺石。陰離子、陽離子捕收劑均可用于浮選高嶺石,因陰離子捕收劑必須預(yù)先活化[4-6],目前多用陽離子捕收劑浮選高嶺石。胡岳華等[7]以一種新型季銨鹽(DTAL)為捕收劑,反浮選脫除鋁土礦中的高嶺石、伊利石、葉蠟石等硅酸鹽脈石礦物。研究表明,相比十二胺而言,DTAL對三種硅酸鹽礦物具有更好的選擇性。李海普等[8]、劉長淼等[9]、徐華龍等[10]、陳攀等[11]、岳彤等[12]、郭靜楠等[13]等借助分子動力學(xué)構(gòu)建陽離子捕收劑在高嶺石表面的吸附模型,在分子水平上解釋實驗中難以觀察和檢測到的吸附結(jié)構(gòu)和動力學(xué)過程[14-18]。高嶺石端面的荷電機理與石英的相同[19],其Si-O、Al-O鍵的斷裂顯示了單純氧化物的行為,可以當(dāng)作是一個石英與一個剛玉的表面。所以,二者的可浮性既有相似之處,又有不同之點[20],用于浮選石英的捕收劑一般均可用于浮選高嶺石。國外浮選石英使用的陽離子捕收劑主要為醚胺、酰胺、多胺、縮合胺及其鹽等[21],對石英單礦物的捕收能力AC1201優(yōu)于十二烷基丙基醚胺和十二胺[22]。余新陽等[23]利用醚胺對高嶺石進行了浮選試驗研究,結(jié)果表明烷氧丙基醚胺對其捕收性能比十二烷胺好。本文嘗試使用AC1201對高嶺石進行浮選,并結(jié)合分子動力學(xué)模比較AC1201與DAH對高嶺石單礦物的捕收效果的強弱。

    1 試驗

    1.1 試驗原料

    試驗中所用高嶺石單礦物選自山西大同塔山煤礦。高嶺石經(jīng)人工手揀塊礦,破碎挑選后,用瓷球磨至粒度小于0.074 mm,篩分備用。礦樣經(jīng)化學(xué)分析(表1)和XRD圖譜分析(圖1),純度在95%以上,達到純礦物研究要求。

    表1 單礦物化學(xué)成分分析結(jié)果(單位:%)

    圖1 高嶺石XRD圖譜

    表1中高嶺石單礦物的燒失量為16.20%,高于其理論值11%~15%,為分析高嶺石燒失量高的原因,對其作了熱重分析,結(jié)果如圖2所示。礦物生料受熱在30 ℃附近有約2%左右的質(zhì)量損失,是礦物表面水分的揮發(fā);250 ℃以下?lián)]發(fā)了0.2%左右的吸附水;250~450 ℃之間有約0.8%左右的損失,是有機組分的燒失。高嶺石的脫水發(fā)生在450 ℃ 以上,有大量結(jié)構(gòu)水的損失。

    圖2 高嶺石熱重圖

    試驗中的捕收劑為十二胺鹽酸鹽(DAH)和十二胺聚氧乙烯醚(AC1201),質(zhì)量濃度為1%的鹽酸(HCl)和氫氧化鈉(NaOH)作為pH調(diào)整劑。其中DAH由十二胺(分析純)與鹽酸混合配置而成(摩爾比為1∶1),化學(xué)純AC1201由江蘇海安石化提供。

    1.2 試驗方法

    1.2.1 浮選試驗

    單礦物浮選試驗在XFGⅡ3-35型掛槽浮選機上進行,主軸轉(zhuǎn)速為1 602 r/min。每次稱取10.0 g單礦物放入60 mL浮選槽中,加入50 mL水,攪拌2 min,用HCl或NaOH調(diào)整pH值2 min,再加入捕收劑攪拌3 min,然后刮泡浮選4 min。泡沫和槽內(nèi)產(chǎn)物于75 ℃烘箱中烘干,稱重,計算回收率。

    1.2.2 動電位測定及紅外光譜測試

    將單礦物樣磨至<5 μm,每次取50 mg置于100 ml燒杯中,加入50 ml蒸餾水,用磁力攪拌2 min,再用HCl或NaOH調(diào)pH值,測定pH值。最后加入浮選藥劑,用磁力攪拌器攪拌6 min,使礦漿充分分散,沉降10 min后取上層稀釋的礦漿注入Zeta電位分析儀的矩形電池內(nèi)進行電位測定。每個樣品測量3次,取其平均值。

    采用德國布魯克公司制造的TENSOR27傅里葉變換紅外光譜儀對單礦物及吸附表面活性劑后礦樣進行光譜掃描。測量時,加入KBr粉料(KBr與待測樣品質(zhì)量比為100∶1),在瑪瑙研缽中研磨至一定粒度并混合均勻;加壓制成片狀樣品,測量波數(shù)范圍為500~4 000 cm-1。測量完后根據(jù)測試數(shù)據(jù),繪制紅外圖譜。

    1.2.3 動力學(xué)模擬

    基于COMPASS力場,模擬主要在Material Studio 6.0軟件中的Modules-Discover模塊下進行。通常,分子動力學(xué)(MD)方法使用由實驗或計算獲得的力場參數(shù),并以經(jīng)典牛頓運動方程為基本原理,所以它們較量子力學(xué)(QM)方法有更小的計算量而適用于更大的體系[24-25]。結(jié)合高嶺石與捕收藥劑作用的實際情況,與高嶺石作用模擬過程中采用周期性邊界條件,構(gòu)建了二層結(jié)構(gòu)的模型:首先,構(gòu)建礦物表面模型,由晶胞模型沿解理面切面得到,本論文選擇高嶺石的(001)面為研究對象[8],對礦物表面模型進行能量優(yōu)化,優(yōu)化過程中固定底層原子,只對表面兩層原子進行馳豫,并將優(yōu)化后的表面晶胞擴展到2.73 nm×1.57 nm×3.64 nm的超晶胞,模擬過程中凍結(jié)礦物表面體系的所有原子;其次,利用Amorphous Cell模塊構(gòu)建包含300個水分子的水分子層,用于研究捕收劑在固液界面的吸附行為[21];再次,將優(yōu)化好的藥劑分子放置到礦物表面,分布COMPASS力場,其中硅氧四面體中的O、Si分別分布為O---Si力場和Si4---O力場,鋁氧八面體中的O、Al分別分布為Al---O力場和Al3---O力場,與Si和Al同時連接的O則均勻分為Si---O力場和Al---O力場[9]。

    通過Discover模塊下的Smart Minimizer方法對體系進行優(yōu)化尋找吸附平衡位置,用 vdw & Coulomb 計算非鍵合作用(靜電力和范德華作用力),smart minimizer最小化法優(yōu)化幾何構(gòu)型,能量收斂精度為 medium。能量收斂后,在298 K,正則系統(tǒng)(NVT)下進行分子動力學(xué)模擬,模擬時間均為 100 ps,步長為1.0 fs,計算完成后,動力學(xué)模擬輸出 200個軌跡結(jié)構(gòu),其中能量最低的構(gòu)型視為吸附最優(yōu)模型,對此模型進行相互作用能[26]計算,見式(1)。

    (1)

    式中:Ecomplex是優(yōu)化后吸附絡(luò)合模型的能量;Esurface和Ecollector分別是被吸附表面模型和吸附劑模型的能量。ΔE的負(fù)值表示吸附體系的穩(wěn)定性,ΔE越負(fù),吸附能越大,吸附體系越穩(wěn)定,越容易吸附于表面模型,ΔE為0或正值,吸附難于發(fā)生。

    2 結(jié)果與討論

    2.1 陽離子捕收劑對高嶺石的浮選行為

    圖3為DAH和AC1201作為捕收劑,藥劑用量為750 g/t時,礦漿pH值對高嶺石單礦物浮選回收率的影響。由圖3可知,pH值為2~6的范圍內(nèi),DAH對高嶺石回收率由54%降到52%左右,AC1201的回收率由59%降到57%左右,對高嶺石的浮選效果均最好,且回收率基本不變。之后隨著pH值的升高,高嶺石的回收率呈現(xiàn)降低的趨勢。當(dāng)pH≤7時,AC1201對高嶺石的捕收能力略強于DAH;而當(dāng)pH>7時,DAH對高嶺石的捕收能力明顯強于AC1201。

    在強酸性條件下,高嶺石具有較高的浮選回收率,但強酸性環(huán)境設(shè)備腐蝕嚴(yán)重,實際生產(chǎn)條件一般控制在pH值為4、5左右的低酸性和中性條件。

    圖4所示為pH=7.0時,捕收劑DAH和AC1201用量對高嶺石單礦物浮選回收率的影響。從圖4可以看出,兩種藥劑對高嶺石的回收率都是隨著藥劑用量的增加而升高。當(dāng)藥劑用量為100 g/t時,兩種藥劑對高嶺石的捕收能力都很低,AC1201捕收高嶺石時有少量泡沫產(chǎn)生,而DAH幾乎沒有;當(dāng)藥劑用量≤750g/t時,AC1201對高嶺石的捕收能力略強于DAH,當(dāng)藥劑用量>1 000 g/t時,隨著藥劑用量的增加,DAH對高嶺石的回收率略強于AC1201,高嶺石回收率達到74%,可能是由于AC1201吸附量已達極限所致。

    圖3 DAH和AC1201為捕收劑時,礦漿pH值對高嶺石浮選回收率的影響

    圖4 pH=7.0時,DAH和AC1201用量對高嶺石回收率的影響

    2.2 陽離子捕收劑對高嶺石Zeta電位的影響

    圖5所示為高嶺石與兩種捕收劑作用前后的Zeta電位的變化關(guān)系曲線,溶液中藥劑用量都控制在750 g/t。從圖5可以看出,高嶺石在pH=4.0至pH=10.0的區(qū)間內(nèi)均荷負(fù)電,其等電點接近pH=3。添加DAH和AC1201后,高嶺石等電點由3.2分別移到3.7和5.2左右,均使高嶺石Zeta電位向正方向偏移;與DAH相比,pH≤7時,AC1201作用下的高嶺石電位較高,對高嶺石電位的影響較為明顯;而pH>8時,AC1201對高嶺石Zeta電位的影響不明顯,電位低于與DAH作用后的高嶺石電位,AC1201對高嶺石的吸附不明顯,與礦漿pH值對浮選結(jié)果的影響相吻合。

    圖5 DAH和AC1201浮選高嶺石Zeta電位與pH值的關(guān)系

    2.3 陽離子捕收劑與礦物作用的紅外光譜分析

    圖6所示為高嶺石及吸附藥劑后礦樣的紅外光譜圖。高嶺石的紅外光譜圖中,937.4 cm-1是內(nèi)羥基的擺動吸收峰,914.2 cm-1處是外羥基的擺動吸收峰,1 112.1 cm-1和1 033.8 cm-1是Si---O鍵兩個伸縮振動收峰,788.8 cm-1和754.3 cm-1是Si---O---Si的對稱伸縮振動吸收峰,694.3 cm-1是Si---O---Al的不對稱伸縮振動吸收峰,540.0 cm-1則是Si---O---Al的伸縮振動吸收峰。

    對比與藥劑作用前后的紅外光譜可知,高嶺石與藥劑作用后沒有出現(xiàn)新的吸收峰,無明顯峰的遷移,但650~950 cm-1之間的吸收峰峰強度均變?nèi)?,且AC1201的降低程度相對較大。表明兩種藥劑在高嶺石表面的吸附均為物理吸附,且AC1201的吸附較強。

    圖6 高嶺石、DAH吸附高嶺石、AC1201吸附高嶺石的紅外光譜圖

    圖7 pH=7.0時,DAH和AC1201在高嶺石表面吸附的初始模型

    圖8 pH=7.0時,DAH和AC1201在高嶺石表面吸附的最終模型

    2.4 單分子動力學(xué)模擬

    分子動力學(xué)模擬可以得到藥劑分子在高嶺石(001)表面的吸附能以及吸附模型。吸附能反映分子與礦物表面的作用能力,能量越大說明結(jié)合越緊密,越有可能是化學(xué)吸附作用,相反就越可能是物理吸附作用。

    經(jīng)能量優(yōu)化和動力學(xué)模擬優(yōu)化,系統(tǒng)中溫度和能量逐漸趨于穩(wěn)定,最終得到pH=7.0時,DAH和AC1201單分子于高嶺石(001)面吸附的最低能量構(gòu)象,其初始結(jié)構(gòu)模型見圖7,最終模型見圖8。

    圖7(a)和圖7(b)分別為pH=7.0時,單個DAH和AC1201分子在高嶺石表面的初始吸附模型,圖8(a)和(b)分別為pH=7.0時,單個DAH和AC1201分子在高嶺石(001)面的最佳吸附模型。對比初始構(gòu)象發(fā)現(xiàn),當(dāng)藥劑分子在高嶺石(001)面發(fā)生吸附作用后,DAH和AC1201分子的疏水碳鏈發(fā)生彎曲、扭轉(zhuǎn),DAH分子距高嶺石表面的最小距離為2.660 ?,AC1201 分子距高嶺石表面的最小距離為2.398 ?。從圖8所示的最佳吸附構(gòu)象中看,DAH分子基本垂直于礦物表面,與高嶺石(001)面的接觸面積十分有限;AC1201分子在高嶺石(001)面上最終吸附時,疏水基與高嶺石(001)面形成一定的夾角,增加了AC1201分子與高嶺石(001)面的接觸面積。

    從表2中的計算結(jié)果可知,水相環(huán)境中,pH=7.0的環(huán)境下,DAH和AC1201對高嶺石(001)面的吸附作用能分別為-90.22 kCal/mol和-105.39 kCal/mol,說明該兩種藥劑均能夠穩(wěn)定地吸附在高嶺石表面。同時也可看出,中性環(huán)境下,AC1201在高嶺石(001)表面的吸附作用能比 DAH的更低,表明此環(huán)境下AC1201比DAH更牢固的吸附于高嶺石(001)面。因此,中性條件下,AC1201對高嶺石的捕收能力強于DAH,與單礦物浮選實驗中,pH=7.0環(huán)境下AC1201對高嶺石捕收能力強于DAH相符。

    表2 DAH和AC1201對高嶺石(001)面的吸附作用能(kCal/mol)

    3 結(jié) 論

    1)單礦物浮選中,AC1201在酸性和中性條件下對高嶺石的浮選效果優(yōu)于DAH,堿性條件下低于DAH。

    2)單礦物浮選中,中性條件下,藥劑用量較低時,AC1201可較有效地浮選高嶺石,而藥劑用量較高時,DAH可比AC1201有效地浮選高嶺石。

    3)紅外光譜分析結(jié)果表明,兩種藥劑在高嶺石礦物表面均為物理吸附,且AC1201的吸附能力較強。

    4) 在動力學(xué)模擬中,中性環(huán)境下,AC1201在高嶺石(001)表面的吸附能比DAH低,且其距離(001)表面距離小,吸附能低,更容易吸附高嶺石表面。

    [1] 李家毓,周興龍,雷力.我國煤系高嶺土的開發(fā)利用現(xiàn)狀及發(fā)展趨勢[J].云南冶金,2009,38(1):23-26.

    [2] 杜振寶,路邁西,丁靖洋.煤系高嶺土資源開發(fā)利用現(xiàn)狀[J].煤炭加工與綜合利用,2010(2):47-49.

    [3] 李曉輝.山西省應(yīng)大力發(fā)展煤矸石煅燒高嶺土精品[J].科技情報開發(fā)與經(jīng)濟,2014,14(11):101-102.

    [4] Zhao S M,Wang D Z,Hu Y H,et al.A series of aminoamides used for flotation of kaolinite[J].Journal of University and Technology Beijing,2005,12(03):208-212.

    [5] Cao X F,Liu C M,Hu Y H.Flotation of kaolinite with dodecyl tertiary amines[J].Journal of Central South University of Technology,2009(5):749-752.

    [6] 張國范,馮其明,盧毅屏,等.油酸鈉對一水硬鋁石和高嶺石的捕收機理[J].中國有色金屬學(xué)報,2001,11(2):298-301.

    [7] Wang Y H,Hu Y H,Chen X Q.Aluminum-silicates flotation with quaternary ammonium salts[J].Transactions of Nonferrous Metals Society of China,2003(3):715-719.

    [8] 李海普,胡岳華,王淀佐,等.陽離子表面活性劑與高嶺石的相互作用機理[J].中南大學(xué)學(xué)報:自然科學(xué)版,2004(2):228-233.

    [9] Liu C.M.,F(xiàn)eng A.S.,Guo Z.X.,et al.Dynamics simulation of tertiary amines adsorbing on kaolinite (001) plane[J].Transactions of Nonferrous Metals Society of China,2011(8):1874-1879.

    [10] 徐華龍,蔣昊,董發(fā)勤,等.DTAC對不同粒級高嶺石浮選行為的影響[J].中國礦業(yè)大學(xué)學(xué)報,2013,42(5):832-837.

    [11] 陳攀,孫偉,岳彤.季鏻鹽在高嶺石(001)面上的吸附動力學(xué)模擬[J].中國礦業(yè)大學(xué)學(xué)報,2014,43(2):294-299.

    [12] 岳彤,孫偉,陳攀.季銨鹽類捕收劑對鋁土礦反浮選的作用機理[J].中國有色金屬學(xué)報,2014,24(11):2873-2878.

    [13] 郭靜楠,王方平.羧基羥肟酸浮選一水硬鋁石高嶺石的理論研究[J].廣州華工,2012,40(11):6-8.

    [14] 羅思崗,王福良.分子力學(xué)在研究浮選藥劑與礦物表面作用中的應(yīng)用[J].礦冶,2009,18(1):1-4.

    [15] Reimer U.,Wahab M.,Schiller P.,M?gel H.-J.Monte Carlo study of surfactant adsorption on heterogeneous solid surfaces[J].Langmuir,2005,21(4):1640-1646.

    [16] Tarek M.,Bandyopadhyay S.,Klein M.L.Molecular dynamics studies of aqueous surfactants systems[J].Journal of Molecular Liquids,1998,78:1-6.

    [17] Leermakers F.A.M.,Koopal L.K.,Goloub T.P.,et al.Confinement induced symmetry breaking of interfacial surfactant layers[J].The Journal of Physical Chemistry Part B:Condensed Matter,Materials,Surfaces,Interfaces & Biophysical,2006,110(17):8756-8763.

    [18] 宋其圣,郭新利.十二烷基苯磺酸鈉在SiO2表面聚集的分子動力學(xué)模擬[J].物理化學(xué)學(xué)報,2009,25(6):1053-1058.

    [19] Huang Z.Q.,Zhong H.,Wang S.,et al.Gemini trisiloxane surfactant:Synthesis and flotation of aluminosilicate minerals[J].Minerals Engineering,2014(56):145-154.

    [20] 魏克武,牟杏妹.高嶺石和石英的分離[J].非金屬礦,1992(2):53-57.

    [21] 劉安.非極性油輔助十二胺混溶捕收劑提效機理研究[D].太原:太原理工大學(xué),2015:67-68

    [22] 劉安,樊民強.十二胺聚氧乙烯醚對石英的捕收性能[J].金屬礦山,2014(10):51-55.

    [23] 余新陽,鐘宏,劉廣義.陽離子反浮選脫硅捕收劑研究現(xiàn)狀[J].輕金屬,2008(6):6-10.

    [24] Fox J R,Andersen H C.Molecular dynamics simulations of a supercolled monatomic liquid and glass[J].J.Phys.Chem,1984,88(18):4019-4027.

    [25] Brown D,Clarke J H R.Molecular dynamics simulation of an amorphous polymer under tension.1.Phenomenology[J].Macromolecules,1991,24(8):2075-2082.

    [26] Pradip,Rai B.Design of tailor-made surfactants for industrial applications using a molecular modeling approach[J].Colloids Surfaces A:Physicochemical and Engineering Aspects,2002,205:139-148.

    Study and dynamics simulation of cationic collecting performance on kaolinite

    GUO Lina,LI Zhihong,ZHU Zhanglei,LIU Yanli,F(xiàn)AN Minqiang

    (College of Mining Engineering,Taiyuan University of Technology,Taiyuan 030024,China)

    Pure kaolinite flotation test,Zeta potential measurements,IR spectra and the MS 6.0 molecular dynamics simulation are conducted to study the flotation behavior and interaction mechanism of dodecyl-amine ethoxylate (AC1201) and dodecy-lamine hydrochloride (DAH) on coal kaolinite.The results show that AC1201 can separate kaolinite effectively at lower concentration of 750 compared with DAH.Based on the MS simulation,the adsorption models of DAH and AC1201 on kaolinite (001) plane were constructed and optimized.It shows that the contact area of AC1201 on kaolinite (001) plane is larger than DAH,the minimum distance between AC1201 and the kaolinite (001) plane is less than DAH.In addition,higher adsorption energy will be released when AC1201 is adsorbed on kaolinite surface.The theoretical computations match remarkably well with the experimental results.Therefore,it can provide a theoretical basis for exploring new reagent in kaolinite flotation.

    coal kaolinite;flotation;dodecyl-amine ethoxylate;dynamics simulation

    2017-01-22

    郭麗娜(1990-),女,碩士研究生,E-mail:511899516@qq.com。

    李志紅(1970-),女,漢族,山西太原人,博士,副教授,目前主要從事礦物浮選方面的教學(xué)科研工作,E-mail:lizhihong9123@163.com。

    TD923

    A

    1004-4051(2017)05-0112-05

    猜你喜歡
    高嶺石收劑電位
    新型捕收劑BK610和BK612在利民選煤廠的應(yīng)用
    電位滴定法在食品安全檢測中的應(yīng)用
    MG-9#捕收劑在極細(xì)粒煤泥浮選中的應(yīng)用
    二氧化碳在高嶺石孔隙中吸附的分子模擬
    內(nèi)蒙高嶺石燒結(jié)性能與應(yīng)用試驗
    山東冶金(2018年6期)2019-01-28 08:14:42
    電鍍廢水處理中的氧化還原電位控制
    高嶺石電子結(jié)構(gòu)的密度泛函理論研究
    淺談等電位聯(lián)結(jié)
    新型捕收劑DTX-1常溫分步浮選東鞍山鐵礦混磁精
    金屬礦山(2014年7期)2014-03-20 14:19:52
    基于二維自適應(yīng)hp有限元的自然電位高精度計算
    久久亚洲国产成人精品v| 91精品伊人久久大香线蕉| 简卡轻食公司| 国产色婷婷99| 大话2 男鬼变身卡| 又大又黄又爽视频免费| 毛片一级片免费看久久久久| 18禁动态无遮挡网站| 成人午夜精彩视频在线观看| 国产乱人视频| 亚洲欧洲日产国产| 亚洲国产精品成人久久小说| 久久精品熟女亚洲av麻豆精品| 精品国产露脸久久av麻豆| 欧美bdsm另类| 91在线精品国自产拍蜜月| 亚洲精品成人av观看孕妇| 国产探花在线观看一区二区| 午夜福利高清视频| 一级毛片电影观看| 久久99热6这里只有精品| 人妻 亚洲 视频| 欧美日韩综合久久久久久| 久久精品国产自在天天线| 超碰av人人做人人爽久久| 最近的中文字幕免费完整| 亚洲天堂国产精品一区在线| 99精国产麻豆久久婷婷| 久久久久久伊人网av| 久久97久久精品| 热re99久久精品国产66热6| 80岁老熟妇乱子伦牲交| 免费高清在线观看视频在线观看| 偷拍熟女少妇极品色| 在线免费观看不下载黄p国产| 亚洲三级黄色毛片| 22中文网久久字幕| 国产视频首页在线观看| 国产黄色视频一区二区在线观看| 国产成人午夜福利电影在线观看| 2021天堂中文幕一二区在线观| 一级毛片电影观看| 国产黄色免费在线视频| 国产黄色免费在线视频| 久久99热这里只频精品6学生| 亚洲精品成人av观看孕妇| 久久这里有精品视频免费| av在线蜜桃| 亚洲国产欧美在线一区| 水蜜桃什么品种好| 在线播放无遮挡| 午夜福利网站1000一区二区三区| 精品一区在线观看国产| 美女高潮的动态| 日韩在线高清观看一区二区三区| av免费观看日本| 色视频在线一区二区三区| 天堂俺去俺来也www色官网| 国产精品不卡视频一区二区| 亚洲欧美日韩无卡精品| 深夜a级毛片| 一级二级三级毛片免费看| 一级黄片播放器| 亚洲三级黄色毛片| 欧美亚洲 丝袜 人妻 在线| 有码 亚洲区| 国产淫片久久久久久久久| 精品视频人人做人人爽| 亚洲真实伦在线观看| 我要看日韩黄色一级片| 免费看日本二区| 精品少妇黑人巨大在线播放| 精华霜和精华液先用哪个| 国精品久久久久久国模美| 国产v大片淫在线免费观看| 国产一区二区三区综合在线观看 | 在线观看一区二区三区| 一区二区三区乱码不卡18| 国产亚洲av片在线观看秒播厂| 亚洲色图av天堂| 最新中文字幕久久久久| 一本久久精品| 久久这里有精品视频免费| 男女边摸边吃奶| 国国产精品蜜臀av免费| 亚洲精品中文字幕在线视频 | 嘟嘟电影网在线观看| 天天一区二区日本电影三级| 亚洲怡红院男人天堂| 欧美成人一区二区免费高清观看| 联通29元200g的流量卡| 神马国产精品三级电影在线观看| 国产精品爽爽va在线观看网站| 麻豆精品久久久久久蜜桃| 中文资源天堂在线| 亚洲欧美成人综合另类久久久| 亚洲高清免费不卡视频| 极品少妇高潮喷水抽搐| 日韩av不卡免费在线播放| 欧美另类一区| 女的被弄到高潮叫床怎么办| 人人妻人人澡人人爽人人夜夜| 精品国产三级普通话版| 在线 av 中文字幕| 80岁老熟妇乱子伦牲交| 少妇丰满av| 欧美xxxx黑人xx丫x性爽| 久久精品国产亚洲av天美| 精品国产一区二区三区久久久樱花 | 亚洲欧美精品自产自拍| 欧美少妇被猛烈插入视频| 亚洲精品成人av观看孕妇| 人妻一区二区av| 亚洲精品日韩av片在线观看| 王馨瑶露胸无遮挡在线观看| 成年女人在线观看亚洲视频 | 欧美 日韩 精品 国产| 成人免费观看视频高清| 日韩av免费高清视频| 久久精品夜色国产| 国产精品.久久久| 久久久成人免费电影| 性色avwww在线观看| 天堂网av新在线| 亚洲精品亚洲一区二区| 美女内射精品一级片tv| 男女无遮挡免费网站观看| 亚洲精品自拍成人| 欧美激情在线99| 97在线人人人人妻| 22中文网久久字幕| 国产免费视频播放在线视频| 天堂网av新在线| 神马国产精品三级电影在线观看| 王馨瑶露胸无遮挡在线观看| 久久久欧美国产精品| 欧美日韩在线观看h| 在线观看免费高清a一片| 极品少妇高潮喷水抽搐| 国产爱豆传媒在线观看| 69人妻影院| 亚洲图色成人| av在线蜜桃| 久久影院123| 香蕉精品网在线| .国产精品久久| 国产男人的电影天堂91| 亚洲成人一二三区av| 亚洲色图av天堂| 欧美精品国产亚洲| 又爽又黄无遮挡网站| 成人亚洲精品一区在线观看 | 久久久精品94久久精品| 欧美极品一区二区三区四区| 欧美xxⅹ黑人| 视频中文字幕在线观看| 真实男女啪啪啪动态图| freevideosex欧美| 日韩欧美精品免费久久| 国产精品不卡视频一区二区| 亚洲一级一片aⅴ在线观看| 人妻少妇偷人精品九色| 欧美人与善性xxx| 国产精品秋霞免费鲁丝片| 日韩亚洲欧美综合| 在线精品无人区一区二区三 | 性插视频无遮挡在线免费观看| 男女那种视频在线观看| 亚洲欧美日韩无卡精品| 精品一区二区三区视频在线| 国产伦理片在线播放av一区| 我的女老师完整版在线观看| 免费看av在线观看网站| 肉色欧美久久久久久久蜜桃 | 国产毛片a区久久久久| 久久久久久伊人网av| 国产 精品1| 国产精品99久久久久久久久| 亚洲av福利一区| 韩国高清视频一区二区三区| 夜夜看夜夜爽夜夜摸| 亚洲精品国产色婷婷电影| 午夜精品国产一区二区电影 | 一个人观看的视频www高清免费观看| 99热这里只有是精品50| 91aial.com中文字幕在线观看| 在线观看免费高清a一片| 国产色婷婷99| 久久精品国产自在天天线| 色婷婷久久久亚洲欧美| 国产成人精品福利久久| xxx大片免费视频| 久久精品久久久久久久性| 成人亚洲精品一区在线观看 | 精品少妇久久久久久888优播| 天堂俺去俺来也www色官网| 亚洲国产日韩一区二区| 国产精品三级大全| 日日啪夜夜爽| 亚洲av国产av综合av卡| 最近手机中文字幕大全| 免费大片黄手机在线观看| 精品久久久久久电影网| 亚洲欧美日韩另类电影网站 | 身体一侧抽搐| 国产成人免费无遮挡视频| 欧美一级a爱片免费观看看| 日本与韩国留学比较| 国产免费一区二区三区四区乱码| 色5月婷婷丁香| 美女国产视频在线观看| 国产乱来视频区| 三级男女做爰猛烈吃奶摸视频| 国产 一区精品| 国产高清国产精品国产三级 | 99久久九九国产精品国产免费| 777米奇影视久久| 香蕉精品网在线| 亚洲欧美精品自产自拍| 久久6这里有精品| 精品视频人人做人人爽| 九九久久精品国产亚洲av麻豆| 亚洲精品第二区| 久久女婷五月综合色啪小说 | 天堂网av新在线| 日本色播在线视频| 亚洲国产色片| 人体艺术视频欧美日本| 女人十人毛片免费观看3o分钟| av专区在线播放| 少妇人妻 视频| 国产黄色视频一区二区在线观看| 亚洲国产高清在线一区二区三| 亚洲av电影在线观看一区二区三区 | 只有这里有精品99| 99久久精品国产国产毛片| 一个人看视频在线观看www免费| 欧美精品国产亚洲| 老师上课跳d突然被开到最大视频| 亚洲一区二区三区欧美精品 | 国产黄色免费在线视频| 国产老妇伦熟女老妇高清| 亚洲最大成人中文| 嫩草影院精品99| 丝袜喷水一区| 国产精品av视频在线免费观看| 国产成人aa在线观看| 国产黄片美女视频| 一级毛片久久久久久久久女| av专区在线播放| 亚洲国产精品专区欧美| 久热久热在线精品观看| 亚洲精品亚洲一区二区| 久久久久精品性色| 精品国产一区二区三区久久久樱花 | 亚洲精品色激情综合| 免费看日本二区| 男的添女的下面高潮视频| 嫩草影院入口| 十八禁网站网址无遮挡 | 菩萨蛮人人尽说江南好唐韦庄| 狂野欧美白嫩少妇大欣赏| 精品国产露脸久久av麻豆| 久久久久久久精品精品| 亚洲精华国产精华液的使用体验| 亚洲av福利一区| 亚洲不卡免费看| 免费高清在线观看视频在线观看| 国产精品久久久久久精品电影小说 | 在线a可以看的网站| 我要看日韩黄色一级片| 免费人成在线观看视频色| 国内揄拍国产精品人妻在线| 香蕉精品网在线| 亚洲精品视频女| 99久久精品国产国产毛片| 久久久久国产网址| 欧美另类一区| 美女被艹到高潮喷水动态| 久久久久久伊人网av| 如何舔出高潮| av播播在线观看一区| av女优亚洲男人天堂| 国产精品一区二区在线观看99| 久久人人爽av亚洲精品天堂 | 久久久久久九九精品二区国产| 精品久久久精品久久久| 在线精品无人区一区二区三 | 亚洲欧美日韩另类电影网站 | 亚洲自拍偷在线| 午夜爱爱视频在线播放| 身体一侧抽搐| 久久久色成人| 欧美日韩综合久久久久久| av在线播放精品| 另类亚洲欧美激情| 亚洲自偷自拍三级| 97在线人人人人妻| 99热网站在线观看| 久久久久久久大尺度免费视频| 天天躁夜夜躁狠狠久久av| 大香蕉97超碰在线| 国产有黄有色有爽视频| videossex国产| 大香蕉久久网| 亚洲av一区综合| 亚洲av中文字字幕乱码综合| 亚洲av福利一区| 又爽又黄a免费视频| 26uuu在线亚洲综合色| 韩国av在线不卡| 精品熟女少妇av免费看| 欧美日韩国产mv在线观看视频 | 美女主播在线视频| 成人二区视频| 97人妻精品一区二区三区麻豆| 亚洲人与动物交配视频| 亚洲激情五月婷婷啪啪| 午夜激情福利司机影院| 日韩av在线免费看完整版不卡| 热99国产精品久久久久久7| 国产精品一区二区在线观看99| 91在线精品国自产拍蜜月| 亚洲国产精品999| 91久久精品国产一区二区三区| av在线app专区| 26uuu在线亚洲综合色| 成人二区视频| 亚洲成人中文字幕在线播放| 国产精品久久久久久精品古装| 男女啪啪激烈高潮av片| 亚洲欧美一区二区三区国产| 亚洲精品乱码久久久久久按摩| 日本欧美国产在线视频| av.在线天堂| 王馨瑶露胸无遮挡在线观看| 国产老妇女一区| 国产免费视频播放在线视频| 亚洲av二区三区四区| 亚洲,一卡二卡三卡| 人妻一区二区av| 在线观看国产h片| 99久久精品一区二区三区| 亚洲自偷自拍三级| 男人添女人高潮全过程视频| 一级a做视频免费观看| 国产色爽女视频免费观看| 国产成人免费观看mmmm| 国产大屁股一区二区在线视频| 日本爱情动作片www.在线观看| 欧美变态另类bdsm刘玥| 搡老乐熟女国产| 18+在线观看网站| 麻豆国产97在线/欧美| 一级爰片在线观看| 亚洲一级一片aⅴ在线观看| av黄色大香蕉| 亚洲精品乱久久久久久| 免费少妇av软件| 亚洲自偷自拍三级| 精品一区二区三卡| 一级a做视频免费观看| 国产淫语在线视频| 精品久久久久久久久av| 亚洲三级黄色毛片| 国产国拍精品亚洲av在线观看| 亚洲人成网站高清观看| 777米奇影视久久| 成人综合一区亚洲| 国产黄a三级三级三级人| 精品久久久噜噜| xxx大片免费视频| 夜夜爽夜夜爽视频| 少妇的逼好多水| 欧美97在线视频| 中文精品一卡2卡3卡4更新| 又黄又爽又刺激的免费视频.| 久久久久久久久大av| av线在线观看网站| 国产精品一及| 久久久久久久亚洲中文字幕| 国产亚洲91精品色在线| 一级毛片aaaaaa免费看小| 精品久久久久久久久av| 在线播放无遮挡| 18禁在线播放成人免费| 日日摸夜夜添夜夜添av毛片| 欧美一区二区亚洲| 国产成人freesex在线| 久久精品久久精品一区二区三区| h日本视频在线播放| 久久午夜福利片| 熟女电影av网| 日韩亚洲欧美综合| 亚洲一区二区三区欧美精品 | 99久久精品热视频| 国产大屁股一区二区在线视频| 18禁裸乳无遮挡动漫免费视频 | 99九九线精品视频在线观看视频| 国产一区二区三区综合在线观看 | 亚洲精品亚洲一区二区| 热re99久久精品国产66热6| 亚洲天堂av无毛| 欧美极品一区二区三区四区| 亚洲av.av天堂| 老师上课跳d突然被开到最大视频| 一本色道久久久久久精品综合| 在线观看一区二区三区激情| 久久久久久国产a免费观看| a级一级毛片免费在线观看| 狂野欧美激情性xxxx在线观看| 国产毛片在线视频| 国产黄片视频在线免费观看| 成年版毛片免费区| 性色avwww在线观看| 永久网站在线| 秋霞在线观看毛片| 美女主播在线视频| 涩涩av久久男人的天堂| 国产成人福利小说| 亚洲国产精品999| 久久精品久久精品一区二区三区| 国产男女超爽视频在线观看| 一本一本综合久久| 神马国产精品三级电影在线观看| 如何舔出高潮| 亚洲av.av天堂| 国产色婷婷99| 亚洲av电影在线观看一区二区三区 | 2018国产大陆天天弄谢| 国产成人午夜福利电影在线观看| 久久精品久久精品一区二区三区| 亚洲成人av在线免费| 一本一本综合久久| 国产精品国产三级专区第一集| 99热6这里只有精品| 久久久亚洲精品成人影院| 亚洲精品第二区| 亚洲精品日韩av片在线观看| 高清在线视频一区二区三区| 日韩av免费高清视频| kizo精华| 久久韩国三级中文字幕| 国产精品久久久久久久电影| 国产成人aa在线观看| 亚洲真实伦在线观看| 久久久久久久大尺度免费视频| 少妇人妻久久综合中文| 欧美97在线视频| 噜噜噜噜噜久久久久久91| 伦精品一区二区三区| 少妇猛男粗大的猛烈进出视频 | 少妇 在线观看| 欧美日韩视频精品一区| 国产在线一区二区三区精| 自拍欧美九色日韩亚洲蝌蚪91 | 真实男女啪啪啪动态图| 女人久久www免费人成看片| 视频中文字幕在线观看| 网址你懂的国产日韩在线| 日韩免费高清中文字幕av| 午夜福利网站1000一区二区三区| 少妇熟女欧美另类| av播播在线观看一区| 九九爱精品视频在线观看| 亚洲精品色激情综合| 国产精品国产av在线观看| 啦啦啦中文免费视频观看日本| 亚洲成色77777| 大香蕉97超碰在线| 亚洲真实伦在线观看| 99热这里只有精品一区| 啦啦啦在线观看免费高清www| 日韩在线高清观看一区二区三区| 亚洲,一卡二卡三卡| 国产精品久久久久久久电影| 亚洲精品成人av观看孕妇| 精品国产乱码久久久久久小说| 亚洲精品一区蜜桃| 少妇的逼好多水| 久久久久久久久久成人| 日韩欧美精品免费久久| 午夜日本视频在线| 国产欧美日韩一区二区三区在线 | 男女那种视频在线观看| 搡老乐熟女国产| 看十八女毛片水多多多| 下体分泌物呈黄色| 五月伊人婷婷丁香| 免费少妇av软件| 欧美性感艳星| 国产男女超爽视频在线观看| 51国产日韩欧美| 亚洲精品日韩av片在线观看| 国产高清三级在线| 嫩草影院精品99| 黄色日韩在线| 简卡轻食公司| 26uuu在线亚洲综合色| 免费少妇av软件| 久久人人爽av亚洲精品天堂 | 七月丁香在线播放| 亚洲精华国产精华液的使用体验| 大片免费播放器 马上看| 国产免费一级a男人的天堂| 国产69精品久久久久777片| 亚洲av成人精品一二三区| 日韩视频在线欧美| 亚洲av一区综合| 在线a可以看的网站| 欧美xxxx性猛交bbbb| 美女主播在线视频| 中文字幕久久专区| 久久97久久精品| 国内精品美女久久久久久| 欧美丝袜亚洲另类| 视频区图区小说| 国产探花在线观看一区二区| 国产高清有码在线观看视频| 国产精品一区二区性色av| 男人狂女人下面高潮的视频| 国产精品人妻久久久影院| 人妻一区二区av| 亚洲综合精品二区| 欧美亚洲 丝袜 人妻 在线| 高清在线视频一区二区三区| 国产高清国产精品国产三级 | 精品久久久久久电影网| 国产又色又爽无遮挡免| 97热精品久久久久久| 亚洲av国产av综合av卡| 又爽又黄a免费视频| 久久久精品欧美日韩精品| 色网站视频免费| 日韩一区二区三区影片| 亚洲精品影视一区二区三区av| 1000部很黄的大片| 亚洲精品国产av蜜桃| 有码 亚洲区| 女人久久www免费人成看片| 国产精品一区www在线观看| 婷婷色综合大香蕉| 精品一区二区三卡| 69av精品久久久久久| 1000部很黄的大片| 狠狠精品人妻久久久久久综合| 亚洲精品乱久久久久久| 国产探花极品一区二区| 国产片特级美女逼逼视频| 黄片无遮挡物在线观看| 天堂网av新在线| 国产高清有码在线观看视频| 中文天堂在线官网| 日本三级黄在线观看| 国产精品爽爽va在线观看网站| 国产男女超爽视频在线观看| 国产乱来视频区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 最新中文字幕久久久久| 欧美变态另类bdsm刘玥| 精品久久久精品久久久| 日韩欧美精品v在线| 黄色欧美视频在线观看| 国产av国产精品国产| 久久久久久九九精品二区国产| 欧美成人午夜免费资源| 国产高清国产精品国产三级 | 爱豆传媒免费全集在线观看| 国产成人freesex在线| 国产精品久久久久久精品电影小说 | 免费电影在线观看免费观看| 国产精品99久久久久久久久| 国产片特级美女逼逼视频| 国内揄拍国产精品人妻在线| 国产精品人妻久久久久久| 一区二区三区免费毛片| 国产黄频视频在线观看| 日本一本二区三区精品| 一区二区三区四区激情视频| 国产高潮美女av| 国产爽快片一区二区三区| 91久久精品国产一区二区三区| 国产有黄有色有爽视频| 蜜桃亚洲精品一区二区三区| 精品久久久久久久久av| 久久久久久久久久人人人人人人| 国产精品国产三级国产专区5o| 亚洲精品色激情综合| 国产成人精品福利久久| 99热6这里只有精品| 亚州av有码| 色综合色国产| 少妇 在线观看| 免费不卡的大黄色大毛片视频在线观看| 国产探花极品一区二区| 日韩 亚洲 欧美在线| 99久久中文字幕三级久久日本| 日韩国内少妇激情av| 国产精品久久久久久精品电影小说 | 久久久久久伊人网av| 内射极品少妇av片p| 国产美女午夜福利| 丰满人妻一区二区三区视频av| 99re6热这里在线精品视频| 久久久久精品性色| 免费av毛片视频| 伦精品一区二区三区| 国产精品一区www在线观看| 国产成人免费观看mmmm| 成人高潮视频无遮挡免费网站| 午夜激情福利司机影院| 亚洲精品色激情综合| 国产白丝娇喘喷水9色精品| 色视频www国产| 国产av不卡久久|