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

    含水量對(duì)PAM/PVA互穿網(wǎng)絡(luò)水凝膠性能的影響

    2015-12-01 11:34:44魏慶華汪焰恩楊明明柴衛(wèi)紅張映鋒
    計(jì)算物理 2015年5期
    關(guān)鍵詞:結(jié)合能氫鍵含水量

    魏慶華,汪焰恩,楊明明,柴衛(wèi)紅,張映鋒

    (西北工業(yè)大學(xué)機(jī)電學(xué)院,西安 710072)

    文章編號(hào):1001?246X(2015)05?0572?07

    含水量對(duì)PAM/PVA互穿網(wǎng)絡(luò)水凝膠性能的影響

    魏慶華,汪焰恩?,楊明明,柴衛(wèi)紅,張映鋒

    (西北工業(yè)大學(xué)機(jī)電學(xué)院,西安 710072)

    應(yīng)用分子動(dòng)力學(xué)方法對(duì)PAM/PVA互穿網(wǎng)絡(luò)水凝膠體系進(jìn)行研究.發(fā)現(xiàn)該水凝膠體系的內(nèi)聚能密度、結(jié)合能隨著水含量的增加而增加,即體系的穩(wěn)定性更強(qiáng).由靜態(tài)力學(xué)性能分析發(fā)現(xiàn)隨體系含水量的增加,彈性系數(shù)、工程模量以及延展性均減小.此外,通過(guò)對(duì)相關(guān)函數(shù)的分析,發(fā)現(xiàn)H2O分子與其周圍原子、官能團(tuán)之間主要是通過(guò)氫鍵相互作用,形成氫鍵的強(qiáng)弱關(guān)系為Owater>OPVA>OPAM>NPAM,與形成氫鍵的可能性(難易)關(guān)系一致.

    聚丙烯酰胺;聚乙烯醇;分子動(dòng)力學(xué);內(nèi)聚能密度;力學(xué)性能

    0 引言

    聚丙烯酰胺(Polyacrylamide,PAM)是一種被廣泛應(yīng)用于制備生物水凝膠的原材料,具有典型三維網(wǎng)絡(luò)結(jié)構(gòu)、無(wú)毒副作用、性能穩(wěn)定等諸多優(yōu)點(diǎn),但是由于機(jī)械強(qiáng)度較低使得PAM的應(yīng)用受到了限制,特別是在人工肌肉、人工軟骨組織等臨床應(yīng)用方面.為了彌補(bǔ)材料的這一缺陷,目前最為有效的方法是通過(guò)引入另外一種材料制備互穿網(wǎng)絡(luò)水凝膠.

    互穿網(wǎng)絡(luò)水凝膠是由2種或2種以上高分子聚合通過(guò)網(wǎng)絡(luò)互穿而形成的具有獨(dú)特性能的聚合物共混物.在互穿網(wǎng)絡(luò)水凝膠中單組分之間有化學(xué)鍵合作用,并保持其各自的基本性能,從而使得其在性能或功能上產(chǎn)生特殊的協(xié)同作用[1-2],實(shí)現(xiàn)兩種材料的優(yōu)勢(shì)疊加,同時(shí)兩種網(wǎng)絡(luò)的互穿會(huì)使得共混體系的機(jī)械性能得到一定的提高[3-4].

    目前針對(duì)互穿網(wǎng)絡(luò)水凝膠的研究已有很多,并取得了大量的研究成果,但都是從實(shí)驗(yàn)及制備工藝的角度開(kāi)展的[5-7],而從理論角度對(duì)互穿網(wǎng)絡(luò)水凝膠的研究目前還鮮有報(bào)道.當(dāng)前,分子動(dòng)力學(xué)(Molecular Dynamic,MD)方法已成為在原子尺度上模擬研究材料微觀結(jié)構(gòu)-性能重要的理論方法,在各領(lǐng)域都得到了應(yīng)用,得到了很多通過(guò)傳統(tǒng)試驗(yàn)方法無(wú)法得到的結(jié)論[8-11].本文通過(guò)MD理論方法對(duì)PAM與聚乙烯醇(Polyvinyl alcohol,PVA)共混形成的PAM/PVA互穿網(wǎng)絡(luò)水凝膠體系進(jìn)行研究,從微觀分子相互作用的角度研究含水量對(duì)互穿網(wǎng)絡(luò)水凝膠性能的影響,為互穿網(wǎng)絡(luò)水凝膠的基礎(chǔ)研究提供理論依據(jù).

    1 模型及模擬方法

    1.1 力場(chǎng)的選擇

    COMPASS[12]力場(chǎng)是在早期第2代CFF9x和PCFF力場(chǎng)的基礎(chǔ)上開(kāi)發(fā)的一個(gè)ab initio力場(chǎng),能準(zhǔn)確計(jì)算物質(zhì)的熱物理性質(zhì),廣泛應(yīng)用于有機(jī)分子、無(wú)機(jī)分子、聚合物的分子模擬.COMPASS力場(chǎng)對(duì)高聚物體系的適用性已得到證實(shí)[13-15],本文采用Materials Studio(MS)程序包中的COMPASS力場(chǎng)對(duì)水凝膠的基本性能進(jìn)行MD模擬.

    1.2 模型及MD模擬

    根據(jù)PAM、PVA分子結(jié)構(gòu)式,通過(guò)MS軟件包中的Visualizer模塊分別構(gòu)建聚合度為25、40的PAM、PVA高分子鏈,端基碳原子加氫以達(dá)到飽和,采用Smart Minimizer方法對(duì)其進(jìn)行優(yōu)化,然后利用Amorphous Cell模塊構(gòu)建不同水含量的PAM/PVA互穿網(wǎng)絡(luò)水凝膠無(wú)定型分子模型.為了減少尺寸效應(yīng),而又不至于使計(jì)算量過(guò)大,所構(gòu)建的無(wú)定型晶胞所含高聚物鏈數(shù)各2條;同時(shí),為了探索水凝膠中水含量對(duì)對(duì)水凝膠性能的影響,分別構(gòu)建含水量質(zhì)量百分比為0%、20%、40%、60%、80%的無(wú)定型晶胞,所構(gòu)建的模型參數(shù)如表1所示.

    表1 無(wú)定型晶胞的參數(shù)Table 1 Parameters of am orphous unit cells

    初始模型通過(guò)MM能量最小化后,在COMPASS力場(chǎng)下,選取NVT正則系綜,采用Andeson恒溫器[16],用Velocity Verlet積分法求解牛頓運(yùn)動(dòng)方程,原子初始速度按Boltzmann隨機(jī)分布方法確定.計(jì)算采用的時(shí)間步長(zhǎng)為1 fs,模擬時(shí)間為150 ps,前100 ps用于平衡運(yùn)算,后50 ps用于計(jì)算分析,模擬過(guò)程中每500步記錄一次體系的軌跡.模擬過(guò)程中VanderWaals和Coulomb分別按Atom based[17]和Ewald[18]方法計(jì)算,截?cái)喟霃綖?.5?.而模擬溫度需要根據(jù)探索的目的進(jìn)行設(shè)定,即在計(jì)算含水量對(duì)水凝膠性能的影響時(shí),溫度設(shè)置為298 K.

    1.3 平衡判別和平衡結(jié)構(gòu)

    體系的平衡可由溫度和能量同時(shí)到達(dá)平衡予以確定.通常當(dāng)溫度和能量在5%~10%范圍內(nèi)波動(dòng)即可認(rèn)為體系已達(dá)到平衡.本研究的所有MD模擬都是在體系達(dá)到平衡的條件下進(jìn)行的.以含水量為40%的PAM/PVA體系在298 K下的MD模擬為例,從其溫度-時(shí)間(圖1)和能量-時(shí)間曲線(圖2)可見(jiàn),二者均已達(dá)到熱力學(xué)平衡.

    圖1 模擬溫度-時(shí)間關(guān)系Fig.1 Temperature vs simulation time

    經(jīng)MD模擬獲得了不同含水量水凝膠體系在298 K下及含水量為40%的水凝膠在不同溫度下的平衡結(jié)構(gòu).為簡(jiǎn)潔計(jì),這里只列出不同含水量PAM/PVA水凝膠在298 K下的平衡結(jié)構(gòu),參見(jiàn)圖3.

    圖2 模擬能量-時(shí)間關(guān)系Fig.2 Energy vs simulation time

    圖3 PAM/PVA互穿網(wǎng)絡(luò)水凝膠在298 K下的平衡結(jié)構(gòu)(PAM、PVA分別標(biāo)記為深色、淺色,其余為水分子.)Fig.3 Equilibrium structures of PAM/PVA hydrogels with differentwater contents at298 K(Deep and light color denotes PAM and PVA,respectively.The rest are watermolecules.)

    2 結(jié)果與討論

    2.1 內(nèi)聚能密度

    內(nèi)聚能(Cohesive Energy)是指1 mol凝聚態(tài)克服分子間相互作用力變?yōu)闅鈶B(tài)所需要的能量.內(nèi)聚能密度(Cohesive Energy Density,CED)是單位體積的內(nèi)聚能,是對(duì)物質(zhì)分子間相互作用力大小的度量,凡與物質(zhì)間相互作用有關(guān)的性質(zhì)和物質(zhì)屬性與內(nèi)聚能密度都具有內(nèi)在聯(lián)系,如物質(zhì)的溶解性、相容性、汽化熱等,

    式中Ecoh是指內(nèi)聚能,V為物質(zhì)的量混合體積.

    為了探索含水量對(duì)水凝膠體系內(nèi)聚能密度的影響,對(duì)所構(gòu)建的PAM/PVA水凝膠平衡體系進(jìn)行了內(nèi)聚能計(jì)算,計(jì)算結(jié)果如表2所示.由表2可以看出,在298 K溫度下水凝膠體系的CED隨著水含量的增加而增加,這是因?yàn)榘殡S著水含量的增加,單位體積內(nèi)所含水分子數(shù)目增加.發(fā)生相互作用的水分子數(shù)增加了,從而形成了更多的氫鍵及靜電力相互作用,致使體系內(nèi)聚能密度得到提高.所以PAM/PVA互穿網(wǎng)絡(luò)水凝膠體系的穩(wěn)定性也會(huì)隨著含水量的增加而加強(qiáng).

    2.2 結(jié)合能

    表2 不同含水量PAM/PVA水凝膠在298 K下的內(nèi)聚能密度Table 2 CED of PAM/PVA hydrogel system sw ith different water contents at 298 K

    在一定溫度下,PAM/PVA互穿網(wǎng)絡(luò)水凝膠的相互作用能,等于PAM/PVA水凝膠體系平衡結(jié)構(gòu)的總能量,減去除掉PAM和水后該結(jié)構(gòu)下的PVA能量,再減去除掉PVA后在該結(jié)構(gòu)下的PAM和水的能量,

    式中,Etotal、EPVA、EPAM+H2O分別表示PAM/PVA水凝膠體系、PVA、PAM+H2O的能量.定義結(jié)合能為相互作用能的負(fù)值,即Ebind=-Einter.結(jié)合能的大小可以反映水凝膠體系中PVA與PAM+H2O產(chǎn)生分子相互作用的強(qiáng)弱,結(jié)合能越大,產(chǎn)生的分子相互作用力越大.

    結(jié)合能由非鍵能和晶格限制能構(gòu)成,而非鍵能占主要部分.由表3可見(jiàn),體系的結(jié)合能隨含水量呈單調(diào)遞增,這是因?yàn)樗康脑龆啻偈贵w系內(nèi)產(chǎn)生了更多的分子間相互作用.水凝膠體系結(jié)合能隨含水量的變化趨勢(shì)與內(nèi)聚能密度隨水含量的變化趨勢(shì)相一致,這與事實(shí)相符,同時(shí)也說(shuō)明了該方法準(zhǔn)確性.

    表3 不同含水量PAM/PVA水凝膠在298 K下的結(jié)合能Table 3 Binding energy of PAM/PVA hydrogel systems w ith different water contents at 298 K

    2.3 力學(xué)性能

    力學(xué)性能可通過(guò)靜態(tài)模型分析獲得,在原子水平計(jì)算中,靜態(tài)模型計(jì)算內(nèi)應(yīng)力張量(σ)[19]采用維里方程

    式中,mi和Vi分別表示原子的質(zhì)量和速度;V0為無(wú)形變時(shí)的體積.

    彈性系數(shù)矩陣由數(shù)值法求解相應(yīng)應(yīng)力與應(yīng)變的一階導(dǎo)數(shù)而得到;采用最小二乘法擬合求得拉伸模量(E)和泊松比γ,進(jìn)而可求得其它有效的各向同性力學(xué)性能,如剪切模量(G)、體積模量(K)、柯西壓(C12-C44)等.對(duì)于各項(xiàng)同性材料,其模量可用拉梅參數(shù)表示[20],

    E、K和G均可用來(lái)衡量材料的剛性強(qiáng)弱,即抵抗彈性形變的能力.柯西壓通常用來(lái)衡量體系的延展性,若為負(fù)值則材料顯脆性,若是正值則表明材料韌性即延展性較好,其值越高表明體系的延展性越好[21-22].

    根據(jù)靜態(tài)力學(xué)原理,對(duì)平衡體系的結(jié)構(gòu)進(jìn)行軌跡分析,可得到PAM/PVA水凝膠在不同含水量體系下的力學(xué)性能,如表4所示.

    表4 不同含水量PAM/PVA水凝膠在298 K溫度下的力學(xué)性能Table 4 Mechanical properties of PAM/PVA hydrogel system sw ith different water contents at 298 K

    由表4可見(jiàn),隨著PAM/PVA水凝膠體系含水量的增加,體系的彈性系數(shù)和各項(xiàng)工程模量呈下降趨勢(shì)說(shuō)明其各向同性增長(zhǎng)和剛性下降.同時(shí)隨著體系含水量的增加,C11,C22和C33組、C12,C13和C23組、C44,C55和C66組內(nèi)的數(shù)據(jù)彼此也更為接近,表明水凝膠含水量的提高可以使材料更為接近各向同性彈性體.此外,柯西壓值隨著含水量的增加而減小,這說(shuō)明體系延展性隨著含水量的增加而有所降低.

    2.4 對(duì)相關(guān)函數(shù)分析

    結(jié)合能的大小體現(xiàn)了組分間相互作用的強(qiáng)弱,但未能對(duì)H2O的分布及其與水凝膠體系中其它組分間的作用方式進(jìn)行分析.為了揭示H2O分子在體系中與其它分子間的作用方式,本文對(duì)平衡結(jié)構(gòu)體系的對(duì)相關(guān)函數(shù)(Pair Correlation Function,PCF)進(jìn)行分析.對(duì)相關(guān)函數(shù)g(r)表示在距離指定參考原子A為r處,B原子出現(xiàn)的概率,

    式中nB表示在距離A原子半徑為r的范圍內(nèi)B原子的數(shù)目,NB表示整個(gè)體系中B原子的總數(shù),V表示整個(gè)體系的體積.

    H2O分子在水凝膠體系主要是通過(guò)氫鍵與PAM、PVA高分子形成相互作用,因此只需要對(duì)可能形成較強(qiáng)相互作用的原子對(duì)相關(guān)函數(shù)進(jìn)行分析.H2O分子作為一個(gè)整體,我們只需要對(duì)氧原子進(jìn)行標(biāo)記即可代表H2O分子,而PAM可能與H2O形成相互作用的主要是羰基(—CO)中的O原子,還有與羰基相連的N原子;PVA中主要是羥基(—OH)中的O原子;此外,還有H2O與周圍其它H2O形成的相互作用.依次將上述原子標(biāo)記為Owater、OPAM、NPAM、OPVA,如圖4所示(圖為H2O分子與高聚物的單體結(jié)構(gòu)).并對(duì)H2O與所標(biāo)記的原子進(jìn)行對(duì)相關(guān)聯(lián)函數(shù)分析,圖5為298K下的含水量為60%的PAM/PVA平衡結(jié)構(gòu)所得的對(duì)相關(guān)函數(shù).

    圖4 對(duì)相關(guān)聯(lián)函數(shù)的原子標(biāo)記示意圖Fig.4 PCF schematic ofmarked atoms

    圖5 水凝膠體系H2O中O原子與近鄰原子的對(duì)相關(guān)函數(shù)Fig.5 PCF between oxygen atoms in H2Omolecules and other atoms nearby

    在PAM/PVA水凝膠體系中分子間相互作用力主要分為氫鍵和范德華作用,徑向分布對(duì)相關(guān)函數(shù)g(r)圖中3.5?以內(nèi)的峰主要由化學(xué)鍵、氫鍵構(gòu)成,3.5?~5?則主要為范德華作用.由圖5可見(jiàn),H2O分子與PAM中羰基(—CO)O原子和N原子、PVA中羥基(—OH)O原子的對(duì)相關(guān)函數(shù)最高峰值范圍均在3.5?以內(nèi),即上述原子、官能團(tuán)主要是通過(guò)氫鍵形成相互作用;而對(duì)相關(guān)函數(shù)峰值在徑向(圖中橫軸)的分布范圍可以反映所形成氫鍵的強(qiáng)弱關(guān)系,H2O分子與上述原子、官能團(tuán)所形成氫鍵的強(qiáng)弱關(guān)系為OPVA>OPAM>NPAM;此外,對(duì)相關(guān)函數(shù)峰值的大小也可以反映出H2O分子與其周圍原子、官能團(tuán)形成氫鍵的可能性(難易),由圖可知,體系中H2O分子與周圍原子、官能團(tuán)形成氫鍵的可能性關(guān)系為Owater>OPVA>OPAM>NPAM,這與前面的強(qiáng)弱關(guān)系相一致.

    為了探索了水含量對(duì)PAM/PVA水凝膠中H2O分子分布的影響,節(jié)選了存在較強(qiáng)相互作用的原子進(jìn)行分析,即H2O分子與(a)PVA羥基(—OH)O原子、(b)PAM羰基(—CO)O原子以及(c)體系中其它H2O分子之間的對(duì)相關(guān)函數(shù).得到了H2O分子在不同含水量水凝膠體系下的對(duì)相關(guān)函數(shù),如圖6所示.

    通過(guò)分析可知,水凝膠體系的含水量對(duì)對(duì)相關(guān)函數(shù)的峰值有一定的影響,最高峰值隨著體系含水量的增加而增大,這主要是因?yàn)殡S著體系含水量的增加,高聚物單元體周圍的H2O分子數(shù)增加,從而形成較為充分的接觸.此外,圖6(c)H2O分子與其周圍其他H2O分子的對(duì)相關(guān)函數(shù)的分布最高峰出現(xiàn)在r=2.76?附近,這與Seung[23]、Seung[24-25]等人的研究結(jié)果一致.

    圖6 不同水含量體系的對(duì)相關(guān)函數(shù)Fig.6 PCFwith differentwater contents

    3 結(jié)論

    通過(guò)分子動(dòng)力學(xué)研究分析了含水量、溫度對(duì)PAM/PVA互穿網(wǎng)絡(luò)水凝膠的影響,得到以下結(jié)論:

    1)體系的內(nèi)聚能密度會(huì)隨著含水量的增加而單調(diào)遞增,物質(zhì)的穩(wěn)定性得到加強(qiáng).從內(nèi)聚能密度的層面闡述了水凝膠體系含水量對(duì)體系穩(wěn)定性的影響.

    2)體系含水量的增加使得高聚物鏈段與H2O分子形成了較為充分的相互作用,致使結(jié)合能隨含水量單調(diào)遞增,這一規(guī)律與體系的內(nèi)聚能密度相一致.

    3)基于靜態(tài)力學(xué)性能分析求得該水凝膠體系的力學(xué)性能,發(fā)現(xiàn)隨體系含水量的增加,彈性系數(shù)和工程模量均減小,表明體系的各向同性遞增,整體剛性下降.由柯西壓值的變化可知,體系的延展性也隨著含水量的增加而增加.

    4)通過(guò)分析H2O分子與PAM中羰基(—CO)O原子和N原子、PVA中羥基(—OH)O原子以及其它H2O分子之間的對(duì)相關(guān)函數(shù),發(fā)現(xiàn)它們之間主要是通過(guò)氫鍵形成相互作用,形成氫鍵的強(qiáng)弱關(guān)系為Owater>OPVA>OPAM>NPAM,同時(shí)形成氫鍵的可能性(難易)與強(qiáng)弱關(guān)系相一致.此外,對(duì)相關(guān)函數(shù)的峰值隨著水含量的增加而增大,這一變化規(guī)律與實(shí)際相符.

    [1] Tian W,F(xiàn)an X D,Chen W X,et al.Research progress in polymeric drug delivery carriers[J].Polymer Materials Science and Engineering,2006,22(4):19-22.

    [2] Rodrigu E IR,F(xiàn)orte M M C,Azambuja D S,etal.Synthesis and characterization ofhybrid polymeric networks(HPN)based on polyvinyl alcohol/chitosan[J].Reactive and Functional Polymers,2007,67(8):708-715.

    [3] Tang Y F,Du Y M,Hu X W,et al.Rheological characterisation of a novel thermo sensitive chitosan/poly(vinyl alcohol)blend hydrogel[J].Carbohydrate Polymers,2007,67(4):491-499.

    [4] Xiao C M,Yang M L.Controlled preparation of physical cross?linked starch?g?PVA hydrogel[J].Carbohydrate Polymers,2006,64(1):37-40.

    [5] Nakajima T,F(xiàn)urukawa H,Tanaka Y,Kurokawa T,et al.True chemical structure of double network hydrogel[J]. Macromolecules,2009,42(6):2184-2189.

    [6] Oldiges C,Tonsing T.Molecular dynamic simulation of structural,mobility effects between dilute aqueous CH3CN solution and crosslinked PAA?part1.Structure[J].Physical Chemistry Chemical Physics,2002,4:1628-1636.

    [7] Tamai Y,Tanaka H,Nakanishi K.Molecular dynamics study of polymer?water interaction in hydrogels.2.Hydrogen?bond dynamics[J].Macromolecules,1996,29:6761-6769.

    [8] Arenaza IM,Meaurio E,Coto B,et al.Molecular dynamics modelling for the analysis and prediction of miscibility inpolylactide/polyvinilphenol blends[J].Polymer,2010,51(19):4431-4438.

    [9] Luo Z L,Jiang JW.Molecular dynamics and dissipative particle dynamics simulations for the miscibility of poly(ethylene oxide)/poly(vinyl chloride)blends[J].Polymer,2010,51(1):291-299.

    [10] Motamedi P,Bagheri R.Investigation of the nanostructure and mechanical properties of polypropylene/polyamide 6/layered silicate ternary nanocomposites[J].Mater Des,2010,31(4):1776-1784.

    [11] WeiQ H,Wang Y E,Yang M M,Wei SM.Molecular dynamics simulation ofmechanical properties and surface interaction for HA/PLA[J].Journal of Functional Materials,2013,44(21):3089-3094.

    [12] Sun H.COMPASS:An ab initio force?field optimized for condensed?phase applications overview with details on alkane and benzene compounds[J].Journal of Physical Chemistry B.1998,102:7338-7364.

    [13] Li M M,Shen R Q,Li F S.Molecular dynamics simulation of binding energies,mechanical properties and energetic performance of the RDX/BAMO propellant[J].Acta Phys Chim Sin,2011,27(6):1379-1385.

    [14] Liao L Q,F(xiàn)u Y Z,Liang X Y,Liu Y Q.Molecular dynamics simulations on miscibility of polyetylene terephthalate/polylactide blends[J].Polymer Materials Science and Engineering,2012,28(5):170-174.

    [15] Jawalkar SS,Aminabhavi T M.Molecular modeling simulations and thermodynamic approaches to investigate compatibility/incompatibility of poly(1?lactide)and poly(vinyl alcohol)blends[J].Polymer,2006,47(23):8061-8071.

    [16] Andersen H C.Molecular dynamics simulations at constant pressure and/or temperature[J].Journal of Chemical Physics,1980,72:2374-2383.

    [17] TosiM P.Cohesion of ionic solids in the Born model[J].Solid State Physics,1964,16:1-120.

    [18] Ewald PP.Evaluation of optical and electrostatic lattice potentials[J].Annals of Physic,1921,64:253-287.

    [19] Swenson R J.Comments on virial theorems for bounded systems[J].Am JPhys,1983,51:940-942.

    [20] Watt JP,Davies G F,Connell R JO.The elastic properties of compositematerials[J].Rev Geophys Space Phys,1976,14:541-563.

    [21] 金日華,華幼卿.高分子物理[M].第二版.北京:化學(xué)工業(yè)出版社,2000.

    [22] Weiner JH.Statisticalmechanics of elasticity[M].Dover Publications Inc,New York,2002.

    [23] Seung G L,Giuseppe F B,Seung S J,David G B.Molecular dynamics simulation study of P(VP?co?HEMA)hydrogels:Effect of water content on equilibrium structures and mechanical properties[J].Biomaterials,2005,30:6130-6141.

    [24] Jang SS,Lin ST,Cagin T,Molinero V,etal.Nanophase segregation and water dynamics in the dendrion diblock copolymer formed from the Frechet polyaryl ethereal dendrimer and linear PTFE[J].Journal of Physical Chemistry B,2005,109:10154-10167.

    [25] Jang SS,Molinero V,Cagin T,Goddard W A.Nanophase?segregation and transport in Nafion 117 from molecular dynamics simulations:Effect ofmonomeric sequence[J].Journal of Physical Chemistry B,2004,108:3149-3157.

    Effects of W ater Content on PAM/PVA Interpenetrating Network Hydrogel Performance

    WEIQinghua,WANG Yanen,YANG Mingming,CHAIWeihong,ZHANG Yingfeng (Mechatronics School,Northwestern Polytechnical University,Xi’an 710072,China)

    To study interpenetrating network hydrogel system of PAM/PVA,molecular dynamics simulation is made to investigate molecular interaction inside a hydrogel system.Effects ofwater contenton PAM/PVA composite hydrogel performance are studied.It is found that cohesive energy density and binding energy of hydrogel system increase with water content increasing.Meanwhile elastic coefficients,engineeringmodulus and ductility decreased with increasing ofwater content.In addition,with analysis of pair correlation function,we found that there are mainly hydrogen bonding interactions between H2O molecules and surrounding atoms or functional groups.Strengths of hydrogen bonds formed are Owater>OPVA>OPAM>NPAM,which consists with possibility(difficulty)of forming hydrogen bond.

    polyacrylamide;polyvinyl alcohol;molecular dynamics;cohesive energy density;mechanical properties

    TB332

    A

    2014-10-18;

    2015-02-12

    國(guó)家自然科學(xué)基金(51175432);教育部高校博士點(diǎn)專項(xiàng)科研基金(20116102110046);中央高?;A(chǔ)研究基金(3102014JCS05007)及111工程(B13044)資助項(xiàng)目

    魏慶華(1988-),男,江西南昌,博士研究生,從事生物材料的設(shè)計(jì)及制備工藝研究.E?mail:weiqinghua@m(xù)ail.nwpu.edu.cn

    ?通訊作者:汪焰恩(1976-),男,江西婺源人,教授,從事人工骨支架的仿生設(shè)計(jì)及制造,E?mail:wangyanen@nwpu.edu.cn

    Received date: 2014-10-18;Revised date: 2015-02-12

    猜你喜歡
    結(jié)合能氫鍵含水量
    教材和高考中的氫鍵
    晶體結(jié)合能對(duì)晶格動(dòng)力學(xué)性質(zhì)的影響
    借鑒躍遷能級(jí)圖示助力比結(jié)合能理解*
    結(jié)合Sentinel- 1B和Landsat8數(shù)據(jù)的針葉林葉片含水量反演研究
    森林工程(2018年4期)2018-08-04 03:23:16
    一次冰雹過(guò)程成雹機(jī)理的數(shù)值模擬
    基于兩相混合流理論P(yáng)EMFC含水量特性分析
    ε-CL-20/F2311 PBXs力學(xué)性能和結(jié)合能的分子動(dòng)力學(xué)模擬
    對(duì)“結(jié)合能、比結(jié)合能、質(zhì)能方程、質(zhì)量虧損”的正確認(rèn)識(shí)
    二水合丙氨酸復(fù)合體內(nèi)的質(zhì)子遷移和氫鍵遷移
    銥(Ⅲ)卟啉β-羥乙與基醛的碳?xì)滏I活化
    欧美在线黄色| 日韩高清综合在线| 美女高潮的动态| 久久久久久久精品吃奶| 国产精品精品国产色婷婷| 国产成年人精品一区二区| 国产精品久久久久久精品电影| 三级男女做爰猛烈吃奶摸视频| 日韩精品中文字幕看吧| 日韩欧美在线乱码| 国产高清三级在线| 怎么达到女性高潮| 久久久精品大字幕| 尤物成人国产欧美一区二区三区| 一级作爱视频免费观看| 99精品久久久久人妻精品| 嫁个100分男人电影在线观看| 亚洲av熟女| 天美传媒精品一区二区| 少妇裸体淫交视频免费看高清| 90打野战视频偷拍视频| 久久精品综合一区二区三区| 国产欧美日韩精品亚洲av| 色播亚洲综合网| 欧美日本视频| 久久久久久久久中文| 国产激情偷乱视频一区二区| 久久伊人香网站| 色尼玛亚洲综合影院| www.熟女人妻精品国产| 中文字幕久久专区| 看黄色毛片网站| 欧美日本视频| 国产熟女xx| 免费av观看视频| 9191精品国产免费久久| 成人精品一区二区免费| 在线观看舔阴道视频| 婷婷精品国产亚洲av| 国产熟女xx| 欧美乱码精品一区二区三区| 少妇高潮的动态图| 一级黄色大片毛片| 国产亚洲精品久久久久久毛片| 麻豆久久精品国产亚洲av| 免费在线观看影片大全网站| 亚洲成人久久性| 每晚都被弄得嗷嗷叫到高潮| 伊人久久大香线蕉亚洲五| 无遮挡黄片免费观看| www日本黄色视频网| 国产精品av视频在线免费观看| 午夜免费成人在线视频| 久久国产乱子伦精品免费另类| 国产真实伦视频高清在线观看 | 欧美日韩乱码在线| 精品免费久久久久久久清纯| 久久久久久久精品吃奶| 免费看十八禁软件| 久久久国产成人免费| 日本撒尿小便嘘嘘汇集6| 在线看三级毛片| 成人鲁丝片一二三区免费| 99久久综合精品五月天人人| svipshipincom国产片| 日本三级黄在线观看| 久久久久精品国产欧美久久久| 欧美大码av| 两个人看的免费小视频| 日日摸夜夜添夜夜添小说| 色综合婷婷激情| 91九色精品人成在线观看| 特级一级黄色大片| 亚洲色图av天堂| 久久香蕉精品热| 午夜两性在线视频| 国产一级毛片七仙女欲春2| 精品一区二区三区人妻视频| 每晚都被弄得嗷嗷叫到高潮| 国产精品综合久久久久久久免费| 91久久精品电影网| 日韩精品青青久久久久久| 午夜福利18| 少妇人妻精品综合一区二区 | 久久国产乱子伦精品免费另类| 精品国产美女av久久久久小说| 免费高清视频大片| 黄色视频,在线免费观看| 毛片女人毛片| 国产亚洲欧美98| 亚洲av成人精品一区久久| 一本综合久久免费| 国产一区二区三区视频了| 国产一区二区在线观看日韩 | 亚洲久久久久久中文字幕| 欧美激情久久久久久爽电影| 91麻豆精品激情在线观看国产| 久久久久国内视频| 国产精品野战在线观看| 搡老岳熟女国产| 亚洲七黄色美女视频| 午夜福利成人在线免费观看| 欧美三级亚洲精品| 变态另类丝袜制服| 少妇人妻精品综合一区二区 | 又紧又爽又黄一区二区| 波多野结衣高清作品| 麻豆国产av国片精品| 国产伦精品一区二区三区四那| 波多野结衣高清无吗| 我要搜黄色片| 黄色丝袜av网址大全| 波野结衣二区三区在线 | 色综合欧美亚洲国产小说| 内射极品少妇av片p| 国产一区在线观看成人免费| 国产亚洲av嫩草精品影院| 在线播放国产精品三级| 国产精品综合久久久久久久免费| 精品久久久久久久毛片微露脸| 日本黄色片子视频| 在线天堂最新版资源| 深夜精品福利| 欧美日韩福利视频一区二区| 国产伦一二天堂av在线观看| 男女下面进入的视频免费午夜| 国产精品 国内视频| 真实男女啪啪啪动态图| 人人妻,人人澡人人爽秒播| 日韩欧美免费精品| 精品人妻一区二区三区麻豆 | 久久亚洲真实| 日日夜夜操网爽| 日韩欧美国产一区二区入口| 嫩草影院入口| 国产色爽女视频免费观看| 12—13女人毛片做爰片一| 噜噜噜噜噜久久久久久91| 首页视频小说图片口味搜索| 成人特级黄色片久久久久久久| 九色成人免费人妻av| 很黄的视频免费| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久久久大av| 给我免费播放毛片高清在线观看| 9191精品国产免费久久| 国产精品香港三级国产av潘金莲| 国产精品女同一区二区软件 | 免费在线观看成人毛片| 一个人看的www免费观看视频| 91麻豆精品激情在线观看国产| 麻豆成人午夜福利视频| 757午夜福利合集在线观看| 99久久精品热视频| 日韩有码中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 亚洲,欧美精品.| bbb黄色大片| 好看av亚洲va欧美ⅴa在| 色噜噜av男人的天堂激情| 国产精华一区二区三区| 欧美+日韩+精品| 男女做爰动态图高潮gif福利片| 中国美女看黄片| 精品国内亚洲2022精品成人| 美女大奶头视频| 中国美女看黄片| 免费一级毛片在线播放高清视频| 国内揄拍国产精品人妻在线| 国产精品美女特级片免费视频播放器| 99国产精品一区二区三区| 久久久成人免费电影| 1024手机看黄色片| 成人高潮视频无遮挡免费网站| 高清日韩中文字幕在线| 国产97色在线日韩免费| 久久香蕉精品热| 久久久国产精品麻豆| xxxwww97欧美| 国产熟女xx| 久久精品国产自在天天线| 给我免费播放毛片高清在线观看| 久久久久国产精品人妻aⅴ院| 国产私拍福利视频在线观看| 少妇的丰满在线观看| 欧美日韩综合久久久久久 | 亚洲精品在线观看二区| 嫁个100分男人电影在线观看| 看黄色毛片网站| 精品国产三级普通话版| 免费在线观看亚洲国产| 久久精品91蜜桃| 九九热线精品视视频播放| 亚洲av不卡在线观看| 中文字幕av在线有码专区| 色噜噜av男人的天堂激情| 亚洲欧美日韩高清专用| 老熟妇乱子伦视频在线观看| 综合色av麻豆| 亚洲欧美日韩高清在线视频| 国产不卡一卡二| 精品欧美国产一区二区三| 内射极品少妇av片p| 欧美黑人欧美精品刺激| 精品99又大又爽又粗少妇毛片 | 又紧又爽又黄一区二区| 观看美女的网站| 国产又黄又爽又无遮挡在线| 久久精品国产自在天天线| 午夜精品久久久久久毛片777| 天天一区二区日本电影三级| 久久久国产精品麻豆| 亚洲人成电影免费在线| 在线天堂最新版资源| 俄罗斯特黄特色一大片| 成人欧美大片| 高清日韩中文字幕在线| 亚洲av电影不卡..在线观看| 精品人妻1区二区| 久久精品91无色码中文字幕| 欧美+日韩+精品| 中文字幕人妻丝袜一区二区| 法律面前人人平等表现在哪些方面| 老熟妇仑乱视频hdxx| 欧美激情在线99| 窝窝影院91人妻| 91av网一区二区| 久久午夜亚洲精品久久| 国产精品98久久久久久宅男小说| 丰满乱子伦码专区| 久久久久久久亚洲中文字幕 | 国产黄色小视频在线观看| 日本一本二区三区精品| 草草在线视频免费看| 亚洲国产色片| 欧美不卡视频在线免费观看| 国产精品精品国产色婷婷| 两个人看的免费小视频| 久久精品国产亚洲av香蕉五月| x7x7x7水蜜桃| 99久久久亚洲精品蜜臀av| 免费搜索国产男女视频| 欧美乱色亚洲激情| 综合色av麻豆| 免费人成在线观看视频色| 国内揄拍国产精品人妻在线| 亚洲人成伊人成综合网2020| 两个人的视频大全免费| 俺也久久电影网| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久毛片微露脸| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 男女下面进入的视频免费午夜| 国产成人aa在线观看| 黄色日韩在线| 日日干狠狠操夜夜爽| 国产探花在线观看一区二区| 宅男免费午夜| 高清在线国产一区| 午夜福利欧美成人| 成年女人毛片免费观看观看9| 日韩欧美 国产精品| 69人妻影院| 国产伦精品一区二区三区四那| 日韩欧美一区二区三区在线观看| 最新在线观看一区二区三区| 此物有八面人人有两片| 在线观看免费午夜福利视频| 国内精品美女久久久久久| 久久亚洲真实| 色老头精品视频在线观看| 国产亚洲欧美在线一区二区| 脱女人内裤的视频| 一进一出好大好爽视频| 日本三级黄在线观看| 亚洲成人精品中文字幕电影| 很黄的视频免费| 不卡一级毛片| 久久性视频一级片| 国产亚洲精品久久久久久毛片| 日本a在线网址| 一区二区三区国产精品乱码| 九色成人免费人妻av| av黄色大香蕉| 99久久精品热视频| 欧洲精品卡2卡3卡4卡5卡区| 一级a爱片免费观看的视频| 三级毛片av免费| 动漫黄色视频在线观看| 桃红色精品国产亚洲av| 在线a可以看的网站| 香蕉丝袜av| 99精品在免费线老司机午夜| 美女高潮喷水抽搐中文字幕| 高清在线国产一区| 黄色女人牲交| 老司机午夜福利在线观看视频| 97碰自拍视频| 一进一出抽搐动态| 亚洲专区中文字幕在线| 国产精品美女特级片免费视频播放器| 91久久精品电影网| av女优亚洲男人天堂| 小说图片视频综合网站| 一个人观看的视频www高清免费观看| 91九色精品人成在线观看| 国产一区二区在线av高清观看| 最近最新中文字幕大全电影3| avwww免费| 桃色一区二区三区在线观看| 午夜福利在线在线| 老司机午夜十八禁免费视频| 国产精品久久久久久亚洲av鲁大| 国产成人av教育| 不卡一级毛片| 国产欧美日韩一区二区精品| 我的老师免费观看完整版| 亚洲精品久久国产高清桃花| 亚洲专区国产一区二区| 亚洲精华国产精华精| 老司机深夜福利视频在线观看| 国产亚洲欧美在线一区二区| 看片在线看免费视频| 日韩欧美三级三区| 真实男女啪啪啪动态图| 亚洲精品色激情综合| 日本在线视频免费播放| 国产免费av片在线观看野外av| 亚洲av不卡在线观看| 在线视频色国产色| www国产在线视频色| 欧美又色又爽又黄视频| 18禁黄网站禁片午夜丰满| 三级国产精品欧美在线观看| 亚洲av免费高清在线观看| 女同久久另类99精品国产91| 18禁裸乳无遮挡免费网站照片| 欧美不卡视频在线免费观看| 老鸭窝网址在线观看| 少妇熟女aⅴ在线视频| 日韩欧美免费精品| 日本熟妇午夜| 亚洲国产精品成人综合色| 一本综合久久免费| 国产国拍精品亚洲av在线观看 | 悠悠久久av| 免费看十八禁软件| 变态另类成人亚洲欧美熟女| www.www免费av| 99久久99久久久精品蜜桃| 中文字幕人成人乱码亚洲影| 又爽又黄无遮挡网站| 此物有八面人人有两片| a在线观看视频网站| 午夜日韩欧美国产| or卡值多少钱| 国产熟女xx| 天天添夜夜摸| 国产免费男女视频| 国产高清有码在线观看视频| 757午夜福利合集在线观看| 国产黄色小视频在线观看| 国产伦在线观看视频一区| a级毛片a级免费在线| 天堂av国产一区二区熟女人妻| 床上黄色一级片| 级片在线观看| 欧美黄色片欧美黄色片| 亚洲欧美日韩高清在线视频| 国产成人福利小说| 国产伦精品一区二区三区四那| 国产成人福利小说| 在线国产一区二区在线| 亚洲狠狠婷婷综合久久图片| 久久久久免费精品人妻一区二区| 国产伦人伦偷精品视频| 99国产精品一区二区三区| 男女之事视频高清在线观看| 日韩国内少妇激情av| 亚洲欧美一区二区三区黑人| 国产精品一区二区三区四区免费观看 | 97超视频在线观看视频| 怎么达到女性高潮| 欧美一区二区国产精品久久精品| 亚洲午夜理论影院| 精品免费久久久久久久清纯| 欧美日韩精品网址| 亚洲精品久久国产高清桃花| 欧美日韩黄片免| 久久久久九九精品影院| 在线观看免费视频日本深夜| 十八禁人妻一区二区| 小说图片视频综合网站| 一本一本综合久久| 精品人妻偷拍中文字幕| 免费在线观看成人毛片| 一级a爱片免费观看的视频| 精品日产1卡2卡| 在线观看一区二区三区| 国产精品永久免费网站| 免费电影在线观看免费观看| 人妻丰满熟妇av一区二区三区| 国产av不卡久久| 成人一区二区视频在线观看| 成人特级av手机在线观看| 看黄色毛片网站| 乱人视频在线观看| 757午夜福利合集在线观看| 欧美在线一区亚洲| 女同久久另类99精品国产91| 1024手机看黄色片| 国产精品久久久久久亚洲av鲁大| www国产在线视频色| 日韩 欧美 亚洲 中文字幕| 精品99又大又爽又粗少妇毛片 | 亚洲国产精品999在线| 久久草成人影院| 最新中文字幕久久久久| 最近最新中文字幕大全电影3| 日本精品一区二区三区蜜桃| 亚洲av中文字字幕乱码综合| 亚洲无线在线观看| 欧美日韩国产亚洲二区| 欧美黑人巨大hd| 三级男女做爰猛烈吃奶摸视频| 一本久久中文字幕| 日本免费一区二区三区高清不卡| 桃色一区二区三区在线观看| 最新中文字幕久久久久| 国产高清三级在线| 成人特级av手机在线观看| 搡老熟女国产l中国老女人| 日本免费a在线| 好看av亚洲va欧美ⅴa在| 91在线精品国自产拍蜜月 | av女优亚洲男人天堂| 97碰自拍视频| 精品日产1卡2卡| 听说在线观看完整版免费高清| 日本在线视频免费播放| 亚洲熟妇中文字幕五十中出| 9191精品国产免费久久| 日韩av在线大香蕉| 国产精品99久久99久久久不卡| 黄片大片在线免费观看| 日韩免费av在线播放| 99久久精品国产亚洲精品| 成人三级黄色视频| 久久伊人香网站| 88av欧美| 久久精品国产综合久久久| 18禁美女被吸乳视频| 亚洲精品一卡2卡三卡4卡5卡| 午夜福利视频1000在线观看| 亚洲专区中文字幕在线| 欧美一级毛片孕妇| 丰满人妻一区二区三区视频av | 国产欧美日韩一区二区三| 国产高清三级在线| 噜噜噜噜噜久久久久久91| 天天一区二区日本电影三级| 久久久久九九精品影院| 伊人久久大香线蕉亚洲五| 国产单亲对白刺激| 99久久精品热视频| 成人欧美大片| 51午夜福利影视在线观看| av天堂中文字幕网| 别揉我奶头~嗯~啊~动态视频| 午夜福利在线在线| 蜜桃久久精品国产亚洲av| 国产免费av片在线观看野外av| 国产色婷婷99| 深爱激情五月婷婷| av福利片在线观看| 99精品久久久久人妻精品| 国产三级在线视频| 国产精品永久免费网站| 麻豆国产97在线/欧美| 无限看片的www在线观看| 天天躁日日操中文字幕| 亚洲成人久久性| 国产亚洲欧美98| 男女做爰动态图高潮gif福利片| 高潮久久久久久久久久久不卡| 亚洲狠狠婷婷综合久久图片| 国产精品99久久99久久久不卡| 真人一进一出gif抽搐免费| 女人高潮潮喷娇喘18禁视频| 欧美国产日韩亚洲一区| 国产精品三级大全| 成人高潮视频无遮挡免费网站| 高清毛片免费观看视频网站| 久久久久亚洲av毛片大全| 亚洲av不卡在线观看| 美女高潮喷水抽搐中文字幕| 搡老熟女国产l中国老女人| 亚洲精华国产精华精| 日韩欧美 国产精品| www国产在线视频色| 国产精品一及| netflix在线观看网站| 免费搜索国产男女视频| 人妻夜夜爽99麻豆av| 日本黄色片子视频| 18禁国产床啪视频网站| 中文字幕av在线有码专区| 亚洲国产精品成人综合色| 人妻夜夜爽99麻豆av| 久久久国产成人精品二区| 91麻豆精品激情在线观看国产| 成人精品一区二区免费| 亚洲黑人精品在线| 国内少妇人妻偷人精品xxx网站| 久久草成人影院| 亚洲av免费在线观看| 国产亚洲精品综合一区在线观看| 91久久精品国产一区二区成人 | 中国美女看黄片| 中文字幕av成人在线电影| 一卡2卡三卡四卡精品乱码亚洲| 美女高潮喷水抽搐中文字幕| 国产伦一二天堂av在线观看| 黄片大片在线免费观看| 日韩欧美国产一区二区入口| 亚洲第一电影网av| 一个人免费在线观看的高清视频| 国产精品99久久99久久久不卡| 国产精品久久视频播放| 免费无遮挡裸体视频| 久久九九热精品免费| 日本五十路高清| 亚洲精品乱码久久久v下载方式 | 国产精品99久久久久久久久| 男插女下体视频免费在线播放| 欧美日本亚洲视频在线播放| 国产又黄又爽又无遮挡在线| 免费大片18禁| 亚洲中文字幕一区二区三区有码在线看| 成人鲁丝片一二三区免费| 亚洲av免费在线观看| 又黄又粗又硬又大视频| 久久久色成人| 欧美bdsm另类| 丁香六月欧美| 亚洲五月天丁香| 神马国产精品三级电影在线观看| 级片在线观看| 免费大片18禁| 身体一侧抽搐| 精品久久久久久成人av| 久久久国产精品麻豆| 三级男女做爰猛烈吃奶摸视频| 免费在线观看亚洲国产| 18禁国产床啪视频网站| 熟妇人妻久久中文字幕3abv| 最近在线观看免费完整版| 搡老岳熟女国产| 99精品在免费线老司机午夜| 欧美黑人巨大hd| 国产精品久久久久久亚洲av鲁大| 午夜久久久久精精品| 亚洲av一区综合| 午夜福利在线观看免费完整高清在 | 亚洲国产中文字幕在线视频| 欧美黄色片欧美黄色片| 免费看光身美女| av女优亚洲男人天堂| 精品国产三级普通话版| 欧美高清成人免费视频www| 久久久久九九精品影院| 亚洲精品日韩av片在线观看 | 亚洲,欧美精品.| 99国产精品一区二区三区| 国产真实乱freesex| 久久亚洲精品不卡| 在线国产一区二区在线| 日本一二三区视频观看| 91av网一区二区| 丁香六月欧美| 午夜老司机福利剧场| 亚洲成人精品中文字幕电影| 男女下面进入的视频免费午夜| 丁香六月欧美| 日韩高清综合在线| 97人妻精品一区二区三区麻豆| 国产欧美日韩精品一区二区| 亚洲精华国产精华精| 日本与韩国留学比较| 人妻丰满熟妇av一区二区三区| 91久久精品国产一区二区成人 | 久9热在线精品视频| 一a级毛片在线观看| 中文字幕av在线有码专区| 桃红色精品国产亚洲av| 啦啦啦免费观看视频1| 少妇的逼水好多| 午夜福利视频1000在线观看| 3wmmmm亚洲av在线观看| 小蜜桃在线观看免费完整版高清| a在线观看视频网站| 五月伊人婷婷丁香| 看黄色毛片网站| 免费观看的影片在线观看| av国产免费在线观看| 午夜免费男女啪啪视频观看 | 国产高清视频在线播放一区| 在线播放国产精品三级| 一夜夜www| 91字幕亚洲| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 久久亚洲真实| 村上凉子中文字幕在线| 窝窝影院91人妻| 嫁个100分男人电影在线观看|