DOI： 10.11817/j.issn.1672-7207.2021.04.005

注塑成型微流控芯片通道脱模的变形机理

王宜磊，翁灿，杨锦，刘介珍，蒋炳炎

(中南大学 机电工程学院，高性能复杂制造国家重点实验室，湖南 长沙，410083)

摘要：为提高微流控芯片注塑成型的脱模质量，采用分子动力学(MD)对环烯烃共聚物(COC)、聚丙烯(PP)、聚甲基丙烯酸甲酯(PMMA)的脱模过程进行模拟，研究在7×10^-11 N脱模外力作用下，聚合物的平均速度、密度分布以及界面相互作用能的变化规律，分析通道脱模变形的分子演化机制。研究结果表明：在脱模过程中，通道底部最早与Ni模芯分离，随后通道肩部开始分离，COC和PP脱模速度最快；随着脱模过程的进行，界面相互作用能先增大后逐渐减小，PMMA的界面相互作用能最大，黏附能最大；在界面相互作用能作用下，聚合物层各部分的运动速度不一致，这是造成通道一定程度拉伸、凹陷和孔隙等变形行为的主要原因；脱模后聚合物整体密度减小和分子链回转半径增大，较好地解释了聚合物整体结构的变化情况。

关键词：注塑成型；微流控芯片；分子动力学；通道形貌；脱模变形

中图分类号：TQ320.66          文献标志码：A            开放科学(资源服务)标识码(OSID)

文章编号：1672-7207（2021）04-1083-07

Mechanism of channel deformation of microfluidic chips during demolding process of injection molding

WANG　Yilei, WENG　Can, YANG　Jin, LIU　Jiezhen, JIANG　Bingyan

(State Key Laboratory of High-Performance Complex Manufacturing, School of Mechanical and Electrical Enginnering, Central South University, Changsha 410083, China)

Abstract: In order to improve the demolding quality of injection molded microfluidic chips, molecular dynamics(MD) method was used to simulate the demolding process of cycloolefin copolymer(COC), polypropylene(PP) and polymethyl methacrylate(PMMA). Under the external force of 7×10^-11 N, the average velocity, density distribution and interfacial interaction energy of the polymer were analyzed to investigate the molecular evolution mechanism on the channel deformation. The results show that the bottom of channel is first separated from the Ni insert, and then the shoulders of the channel begin to separate. The demolding speeds of COC and PP are faster than those of PMMA. The interaction energy between PMMA and Ni is the largest with the highest adhesion energy and the longest interaction time, and the interface interaction energy first increases and then gradually decreases in the demolding process. The inconsistency of the velocity in each part of the polymer layer under the action of the interfacial interaction energy is the main cause of the deformation behaviors of channel, such as stretching, depression and pores. The change that the polymer overall density decreases and molecular chain radius of gyration increases after demolding explains the reason of the change in the overall structure of the polymer.

Key words: injection molding; microfluidic chips; molecular dynamics; channel morphology; demolding deformation

微流控芯片在生命科学、化学分析和医疗检测等研究领域中得到了越来越广泛的应用^[1]。微注射成型技术是目前微流控芯片制造的主要方法之一。微流控芯片微通道的作用主要是进样和分离样品，微通道的成型质量将直接影响芯片的功能和特性。注射成型过程主要包括充填、保压、冷却和脱模4个阶段，其中，脱模阶段对于微结构零件的成型和高精度复制至关重要。闵丽萍等^[2]发现微通道的复制度受模具温度和熔体温度影响显著。聚合物微结构在脱模阶段主要受到金属模芯壁面的摩擦力^[3]、黏附力^[4]和由热收缩系数差异所引起的热收缩力^[5]作用，导致微结构易出现断裂、拉伸、表面毛刺等缺陷^[6-7]。当零件的特征结构尺度缩小至纳米量级时，表面效应和界面效应将产生重要影响，分子间力或表面力的作用十分重要^[8]。经过数十年的发展，分子动力学(MD)模拟方法在纳米压印成型(NIL)^[9]、界面机理^[10]、复合材料^[11]等研究领域得到了广泛应用。TAKAI等^[12]建立了纳米压印(NIL)脱模过程的MD模型，并分析了MD理论和连续介质力学理论这2种模拟方法之间的差异。YANG等^[13]通过MD模拟研究了模芯与基底材料对纳米压印工艺中界面黏附特性的影响。KANG等^[14]使用MD模拟对比分析了NIL脱模过程中不同深宽比纳米结构受到的黏附力和摩擦力。对于微注射成型技术，ZHOU等^[15]利用MD方法研究了聚合物在纳米腔中的充填机理；吴辉等^[16]采用MD方法研究了不同模温下塑件表皮层形成的分子演化机制。近年来，MD模拟方法在聚合物充填纳米凹槽结构研究方面得到了成功运用，但对于微流控芯片中通道脱模变形行为机理的研究目前尚未见相关报道。为此，本文作者采用MD方法模拟了微流控芯片通道注塑成型的脱模过程。通过分析不同脱模时刻的通道形貌以及脱模过程中聚合物的平均速度、界面相互作用能和分子链回转半径等变化情况，研究通道的复制质量和脱模变形机理，旨在为微流控芯片的无损成型提供理论依据。

1  材料与方法

1.1　模型构建

选择热塑性较好的环烯烃共聚物(COC)、聚丙烯(PP)、聚甲基丙烯酸甲酯(PMMA)作为研究材料，构建上层聚合物、下层镍金属模芯的注射成型初始模型，如图1所示。从Materials Studio(MS)软件中导入Ni晶胞，并切割出晶胞的(100)表面，对(100)表面进行扩展得到长×宽×高为10.5 nm×4.5 nm×3.0 nm的超晶胞结构。去除超晶胞两侧区域部分Ni原子，建立了宽度为4.5 nm、深度为1.5 nm的纳米凸起结构。为了便于分子链在脱模模拟过程中发生滑移和解缠，以COC为例，建立聚合度为20的分子链，对单链能量最小化后建立COC聚合物的无定型态体系。为了减小聚合物体系的内应力，得到结构稳定的聚合物模型，先后对聚合物体系进行能量最小化、循环退火及高温驰豫模拟。最后，将聚合物的温度升高至550 K，得到熔融状态的聚合物体系。PP与PMMA聚合物模型的建立方法类似，具体的模型建立参数如表1所示。由于聚合物/Ni模芯界面体系采用周期性边界条件，为避免在z方向上受到周期性边界条件的影响，在聚合物层的上方添加厚度为2.0 nm的真空层。

图1　MD模拟中用于微注射成型过程的聚合物-镍模型

Fig. 1　Polymer-Ni model in MD simulation for micro-injection molding process

表1　3种聚合物结构参数

Table 1　Construction parameters of three polymers

1.2　分子间相互作用势

采用一致性力场PCFF分配聚合物和Ni模芯原子的势参数^[17]。聚合物与Ni模芯界面的相互作用能主要由范德华能和静电相互作用能组成^[18]，分别采用12-6 Lennard-Jones势和库仑相互作用势来描述，其计算式为

(1)

式中：E_non-bonded为界面相互作用能；E_vdw为范德华能；E_ele为静电相互作用能；为非键相互作用常数；为两原子间的平衡距离，nm；为两原子间任意时刻的距离，nm；q_i和q_j为第i个和第j个原子的电荷量；为截断半径，nm。截断半径设置为0.95 nm，当原子间距离超过截断半径时，非键相互作用能可忽略不计。

不同聚合物与金属Ni模芯界面的相互作用能可近似为黏附能，其计算式为

(2)

式中：E_adhesion为黏附能；E_interaction为聚合物和金属界面的相互作用能；E_total为聚合物-金属体系总能量；E_polymer为除去金属的聚合物总能量；E_Ni为除去聚合物的金属表面能。

1.3　分子模拟过程

使用开源软件LAMMPS(large-scale atomic/molecular massively parallel simulator)实现聚合物整个充填和脱模过程的分子动力学模拟。在仿真过程中，在x和y方向上选用周期性边界条件，在z方向上选用自由边界条件。采用等温等容(NVT)正则系综对聚合物体系进行控制^[19]，模拟时间步长设置为0.1 fs，模拟总步数25 000步。3种聚合物材料COC，PP和PMMA的注射成型充填温度分别设置为553，500和500 K，保压时间均设置为1.0 ps。选择充填模拟的最终结果作为脱模过程的初始阶段，在z方向上对聚合物的每个原子施加7×10^-11 N的脱模外力，体系的脱模温度均设置为350 K。

2  结果与讨论

2.1　脱模过程的分子模拟

3种聚合物通道的脱模过程图像如图2所示。从图2可见：0 ps时，聚合物处于完全充填状态；在0.6 ps前，聚合物和Ni模芯之间并没有明显分离，但在脱模外力作用下聚合物层已经沿z方向产生拉伸；0.6 ps时，COC-Ni界面处开始出现孔隙；随着脱模进行，聚合物通道的底部先与Ni模芯分离，随后通道肩部开始分离，而且整个聚合物层的运动速度逐渐变快；PP通道的肩部在1.1 ps时已经与Ni模芯分离，此时，COC和PMMA的通道肩部并没有与Ni模芯完全分离；1.8 ps时，COC和PP近乎同时完成脱模，由于分子运动的时间依赖性，聚合物体系与金属Ni模芯表面的非键相互作用能在界面分离之后继续作用一段时间；1.1 ps时，通道肩部与模芯分离，肩部的分子链结构继续发生变化；1.8 ps时，COC完成脱模，肩部的孔隙变形增大，通道破坏严重；PP通道表面出现凹陷和轻微的表面毛刺。这可能是界面相互作用能较大，PMMA在通道肩部与模芯分离之后，通道内侧壁面的拉伸变形逐渐增大，直到2.3 ps时也未能完全脱模，通道肩部与模芯表面仍有黏连，但这3种聚合物在脱模后都基本保持了通道的整体结构形貌。

图2　MD模拟COC，PP和PMMA脱模过程图像

Fig. 2　Snapshots of MD simulation of demolding process for COC, PP and PMMA

2.2　速度分布

为了研究微流控芯片通道的脱模变形机理，对COC，PP和PMMA整体的原子运动平均速度和通道两侧结构部分的原子运动平均速度分布进行计算，如图3所示。由于界面相互作用力的影响，脱模前期3种聚合物通道两侧区域的原子运动速度比较慢。金属模芯接触处原子的运动速度与内部原子速度不同是造成通道脱模变形的主要原因。由于脱模过程中施加在聚合物原子上的脱模外力大于聚合物本身的内聚力，又小于聚合物-金属界面的黏附力，导致COC通道肩部和底部的聚合物分子链与内部分子链发生严重分离。1.0 ps前，通道左右两侧区域的原子运动平均速度明显低于聚合物整体的原子运动平均速度，此时，PP通道肩部与Ni模芯分离之后出现了凹陷变形；0.2 ps后，PP整体的原子运动平均速度大于COC和PMMA聚合物整体的原子运动平均速度，略早于COC完成脱模，通道出现的拉伸变形量也最小。PMMA结构致密使得其整体的原子运动平均速度较低，因而，在0.2 ps后，PMMA整体的原子运动平均速度明显比COC和PP的平均速度低。当通道肩部和底部的分子链与金属模芯分离后，其两侧区域的原子运动平均速度开始增大，使得0.5 ps后聚合物整体的原子运动平均速度缓慢增大。当聚合物完全脱离模芯后，整个聚合物的原子运动速度基本一致。

图3　聚合物整体平均速度和通道两侧部分平均速度

Fig. 3　Average velocities of polymer layers and sides of channel for each polymer

2.3　界面相互作用能

为了更好地研究脱模过程中脱模变形行为的形成机理，对3种聚合物材料的界面相互作用能进行计算，结果如图4所示。图4中，界面相互作用能为负值表示聚合物-金属模芯表面之间存在的是黏附力，其绝对值越大，表明黏附力越大。由于脱模外力的作用使接触界面处的聚合物分子链受到拉扯，导致脱模前期接触界面的黏附力不断增大。随着脱模进行，通道的底部优先于肩部与Ni模芯分离，一些聚合物分子与金属界面分离而向上运动，通道底部聚合物分子与金属原子之间的距离增大，界面相互作用能开始减小，黏附力减小，此时，通道肩部和壁面处的界面相互作用能起主导作用。在0.2 ps左右时，COC，PP和PMMA接触界面的相互作用能达到峰值，分别为4.31×10^-17，3.53×10^-17和4.70×10^-17 J，且PMMA-Ni模芯界面的相互作用能最大；0.2 ps后，聚合物与Ni模芯界面的相互作用能开始减小，通道两侧结构部分的原子运动速度逐渐增大，使得整个聚合物的平均速度出现缓慢增大现象。结果表明，随着脱模的进行，聚合物分子链结构不断舒展，分子链被拉伸。3种聚合物通道脱模过程中的界面相互作用能变化情况基本一致，与聚合物材料的种类无关。

图4　3种材料界面相互作用能

Fig. 4　Non-bonded interaction energies of three interfaces

2.4　密度分布

沿z方向对脱模前和脱模后的聚合物结构进行水平分层，计算片层厚度为0.1 nm的聚合物平均密度，如图5所示。从图5可见：脱模后，3种材料片层密度都明显降低，且密度变化主要发生在脱模前期，由于聚合物整体结构被拉伸、体积增大，使通道部分聚合物的片层密度明显低于平均片层密度；充填之后聚合物分子在Ni模芯表面聚集，导致通道肩部和底部的聚合物片层密度发生突变，分别出现2个密度峰值。经计算得出，COC通道肩部和底部的片层密度峰值分别为1.72 g/cm³和1.85 g/cm³，PP片层密度峰值分别为1.62 g/cm³和1.90 g/cm³，PMMA片层密度峰值分别为1.48 g/cm³和1.86 g/cm³，且通道底部的片层密度峰值略比肩部的高。结果表明，脱模后越靠近通道的肩部，聚合物的片层密度就越低，结构就越松散。比较脱模前后通道部分的聚合物片层高度(如图5中虚线所示)，发现PMMA脱模后通道两侧结构部分的聚合物结构拉伸变形量最大，这是界面作用能和黏附力作用时间较长所致。

图5　脱模前和脱模后COC，PP和PMMA的片层密度

Fig. 5　Density in each slice of COC, PP and PMMA before and after demolding

为了更深入地研究聚合物脱模过程中的密度变化情况，引入回转半径(R_g)来表征分子链的延展程度，其计算公式为

(3)

式中：M为聚合物所有原子的总质量；m_i为第i个原子的质量；r_i为第i个原子的位置，nm；r_cm为聚合物的重心位置，nm。

脱模过程中聚合物分子链的回转半径如图6所示。从图6可见：0 ps时，COC，PP和PMMA的分子链回转半径分别为3.83，3.81和3.78 nm；在脱模过程中，由于界面相互作用能的影响使聚合物的分子链结构不断舒展；1.1 ps前，3种聚合物分子链的回转半径不断增大；随着脱模进行，聚合物-Ni模芯间的界面相互作用能不断减小，分子链回转半径的增加速率减小，最后趋于稳定值；脱模结束后，COC，PP和PMMA的分子链回转半径分别为4.00，3.99和3.93 nm。PMMA由于结构致密，分子间内聚能密度较大，脱模后其分子链的回转半径改变较小。相比于脱模前，这3种聚合物分子链的回转半径均增大，聚合物分子间距离变大，较好地解释了脱模后聚合物结构密度发生变化的原因。

图6　3种聚合物的回转半径

Fig. 6　Gyration radii of three polymers

3  结论

1) 在7×10^-11 N脱模外力下，通道底部先于肩部与Ni模芯发生分离，其肩部分离后脱模速度增大。COC和PP脱模速度最快，PMMA脱模用时最长。

2) 界面相互作用能的存在使得接触界面和边缘处原子的运动速度与内部原子的运动速度不一致，这是导致通道脱模变形的主要原因。在脱模过程中，COC通道肩部的分子链相互分离产生孔隙，变形最严重；PMMA的界面相互作用能最大，为4.70×10^-17 J，通道内侧壁面产生的拉伸变形最严重；PP纳通道脱模后出现肩部凹陷和轻微的表面毛刺。

3) 充填后，聚合物分子在Ni模芯接触表面聚集导致通道底部和肩部的片层密度发生突变；脱模后，聚合物结构均沿z方向拉伸，分子链结构舒展，回转半径增大，聚合物结构密度减小。

参考文献：

[1] 孙薇, 陆敏, 李立, 等. 微流控芯片技术应用进展[J]. 中国国境卫生检疫杂志, 2019, 42(3): 221-224.

SUN Wei, LU Min, LI Li, et al. Application progress on microfluidic chip technology[J]. Chinese Journal of Frontier Health and Quarantine, 2019, 42(3): 221-224.

[2] 闵丽萍, 蒋炳炎, 朱来余, 等. 微流控芯片微通道复制度的表征及在注塑成型中的应用[J]. 中南大学学报(自然科学版), 2019, 50(9): 2091-2099.

MIN Liping, JIANG Bingyan, ZHU Laiyu, et al. Characterization of micro-channel replication in microfluidic chip and its application in injection molding[J]. Journal of Central South University(Science and Technology), 2019, 50(9): 2091-2099.

[3] BAKHAREV A, ASTBURY D. Simulation of ejection of moldings using combination of mold filling and thermo-mechanical analyses[C]. AIP Conference Proceedings, 2013, 1532(1): 555-560.

[4] ZHOU Mingyong, XIONG Xiang, DRUMMER D, et al. Interfacial interaction and joining property of direct injection-molded polymer-metal hybrid structures: a molecular dynamics simulation study[J]. Applied Surface Science, 2019, 478: 680-689.

[5] KIM B, MIN J. Residual stress distributions and their influence on post-manufacturing deformation of injection-molded plastic parts[J]. Journal of Materials Processing Technology, 2017, 245: 215-226.

[6] STORMONTH-DARLING J M, PEDERSEN R H, HOW C et al. Injection moulding of ultra high aspect ratio nanostructures using coated polymertooling[J]. J Micromech Microeng, 2014, 24(7): 486-514.

[7] CHEN Jianyu, HWANG S J. Design and fabrication of an adhesion force tester for the injection moulding process[J]. Polymer Testing, 2013, 32(1): 22-31.

[8] TAKEDA M, IKEDA M, SATOH S, et al. Rab13 is involved in the entry step of hepatitis C virus infection[J]. Acta Medica Okayama, 2016, 70(2): 111-118.

[9] KIM S, LEE D E, LEE W I. Molecular dynamic simulation on the effect of polymer molecular size in thermal nanoimprint lithographic(T-NIL) process[J]. Tribology Letters, 2013, 49(2): 421-430.

[10] 熊孝东, 魏正英, 杜军, 等. 纳米压印中聚合物界面行为的分子动力学[J]. 高分子材料科学与工程, 2011, 27(6): 186-190.

XIONG Xiaodong, WEI Zhengying, DU Jun, et al. Molecular dynamics of polymer's interface behaviour in the nanoimprint lithography[J]. Polymer Materials Science & Engineering, 2011, 27(6): 186-190.

[11] 李健, 杨延清, 罗贤, 等. 分子动力学模拟在复合材料界面研究中的进展[J]. 稀有金属材料与工程, 2013, 42(3): 644-648.

LI Jian, YANG Yanqing, LUO Xian, et al. Developments of molecular dynamics simulation for composites' interface[J]. Rare Metal Materials and Engineering, 2013, 42(3): 644-648.

[12] TAKAI R, YASUDA M, TOCHINO T, et al. Computational study of the demolding process in nanoimprint lithography[J]. Journal of Vacuum Science & Technology B, 2014, 32(6): 06FG02.

[13] YANG S, YU S, CHO M. Influence of mold and substrate material combinations on nanoimprint lithography process: MD simulation approach[J]. Applied Surface Science, 2014, 301: 189-198.

[14] KANG J H, KIM K S, KIM K W. Molecular dynamics study of pattern transfer in nanoimprint lithography[J]. Tribology Letters, 2007, 25(2): 93-102.

[15] ZHOU Mingyong, JIANG Bingyan, WENG Can. Molecular dynamics study on polymer filling into nano-cavity by injection molding[J]. Computational Materials Science, 2016, 120: 36-42.

[16] 吴辉, 辛勇. 注塑成型聚合物表皮层充填行为的分子机制[J]. 高分子材料科学与工程, 2017, 33(12): 92-98.

WU Hui, XIN Yong. Filling behavior and molecular mechanisms of polymer epidermis in the injection molding process[J]. Polymer Materials Science & Engineering, 2017, 33(12): 92-98.

[17] SUN H. Force field for computation of conformational energies, structures, and vibrational frequencies of aromatic polyesters[J]. Journal of Computational Chemistry, 1994, 15(7): 752-768.

[18] LI Chunyu, STRACHAN A. Cohesive energy density and solubility parameter evolution during the curing of thermoset[J]. Polymer, 2018, 135: 162-170.

[19] ZHANG Meili, XIN Yong. Molecular mechanism research into the replication capability of nanostructures based on rapid heat cycle molding[J]. Applied Sciences, 2019, 9(8): 1683.

(编辑  陈灿华)

收稿日期： 2020 -06 -12; 修回日期： 2020 -09 -23

基金项目(Foundation item)：国家自然科学基金国际(地区)合作与交流项目(51920105008)；国家自然科学基金资助项目(51775562) (Project(51920105008) supported by the International Cooperation and Exchange Foundation of the National Natural Science Foundation of China; Project(51775562) supported by the National Natural Science Foundation of China)

通信作者：蒋炳炎，博士，教授，从事微型构件的精密成型与可控制造研究；E-mail: jby@csu.edu.cn

引用格式： 王宜磊, 翁灿, 杨锦, 等. 注塑成型微流控芯片通道脱模的变形机理[J]. 中南大学学报(自然科学版), 2021, 52(4): 1083-1089.

Citation: WANG Yilei, WENG Can, YANG Jin, et al. Mechanism of channel deformation of microfluidic chips during demolding process of injection molding[J]. Journal of Central South University(Science and Technology), 2021, 52(4): 1083-1089.