工艺与检测Techno1I0 gy and CO2激光直写PMMA微流道凸起仿真研究 相恒富 (中国石油大学机电工程学院,山东东营257061) 摘要:建立了激光烧蚀流道凸起模型,获得了流道凸起与激光功率和作用时间的关系。模拟结果表明:增 大激光扫描速度,减小激光照射时间,选择较低的激光功率可以减小流道凸起高度。模拟结果对激 光直写PMMA微流道表面质量控制有指导意义。 关键词:激光加工微流道有限元仿真 凸起 中图分类号:TP391.9,TK124 文献标识码:A Simulation research on bulge formation by CO2 laser micromachining PMMA substrate microchannel XIANG Hengfu (College ot。Eleetromechanical Engineering,China University of Petroleum,Dongying 25706 1,CHN) Abstract:In this paper a bulge model was gained and then the relation between bulge and laser power and process— ing time.The simulation results show that increasing laser scanning speed,reducing the irradiation time and preferring lower laser power can reduce the bulge height.The simulation results can direct the im— provement of surface quality for laser machining PMMA microchanne1. Keywords:Laser Machining;Microchannel;Finite Element Simulation;Bulge 聚合物PMMA(聚甲基 烯酸甲酯,也称有机玻 璃)微流道的CO 激光直写加工技术,是近年来国际上 的研究热点。一些专家运用实验方法列‘CO 激光直写 PMMA微流道进行了研究。丹麦工业大学微系统研究 1 CO:激光直写PMMA微流道凸起有限元模 型 CO 激光直写PMMA微流控芯片的方法是利用红 所通过大量实验结果,获得了CO!激光加工PMMA微 流控芯片的流道深度模型,并分析了流道宽度、流道长 外激光的高能束在PMMA基片上直接加工微流道,加 工原理如图1所示。 度对流道深度的影响 nf 3 里 0 。Martin F.Jensen等人应用 当激光束聚焦在PMMA 光栅扫描技术优化激光光斑来提高微流道的加工质 但是,CO,激光直写出的流道表面粗糙度值较大, 在CO,激光产生的高温条件下,流道表面的材料发生 沸腾和起泡,冷凝后形成不规则的凸起和凹陷,在芯片 流道边缘会形成凸起的冷凝区带,直接影响芯片在该 区域的键合质量,并且可能导致流道的堵塞。针对流 基片的表面,升温开始,到达 大约115 玻璃态温度之前 一… {{{ 、 i MMA 直是固态。继续升温,PM. MA变成呵翅和橡皮状并开 ~ { 一 始热降解,PMMA的聚合链 } { _// i 断裂发展为单体MMA并挥 l ; 发出来,热降解主要发生在熔鬯 —型Lj 温度为370℃。CO2激光器图l激光加工PMMA原理图 以一定的速度在PMMA基片 道边缘凸起问题,Chung等人对流道凸起形成机理进 行了研究,认为流道边缘凸起是加工过程中熔化的材 料在气化物压力及表面张力作用于流道的结果,他们 通过在PMMA材料表瓶覆盖一层保护层来抑制流道 凸起 。为了预测流道凸起的成形规律,本文运用有 上根据给定的轨迹就可制作成微流控芯片的微流道。 为了用数学模型来描述激光烧蚀材料的过程,不 影响最终结果,我们做了如下假设: 限元方法研究_r激光功率、作用时问与流道凸起高度 的关系,对改善微流道表面质量提供依据。 (1)密度、比热和热传导系数等热物理参数不随 温度变化; 睁.,一 j £I U I耋耳 0删 Techn010gy and Tesf工艺与检测 (2)对比入射光强,材料表面的对流和气化作用 叮以忽略不计; 环境温度取20℃。降解潜热L:1.8X10。J/kg,弹性 模量E=3 300 MPa,泊松比 :0.4,材料的热膨胀系 (3)激光材料相互作用导致在气化温度时去除材 料,任何熔融材料都能被气流完全驱除,不会影响激光 与材料的作用 数 =7×10。、。K~。CO 激光 聚焦光斑半径为I20 m,功率在0~30 w变动。有限元方法用来研究激光 直写PMMA的流道形状和温度分布,通过调整激光功 率和扫描速度等JJU. ̄参数来研究微流道的影响因素。 激光功率遵循高斯分布,激光功率密发为 .,: C-毫 (1) 图3是激光功率为20 w,照射时间分别为0.2 qTW0 式中: 。为激光束的聚焦半径;P为激光功率;“为 PMMA吸收系数;r为点距离激光中心的距离 流道边缘凸起与激光烧蚀流道的温度场有关。这 里先用有限元分析软件建立流道截面的温度场模型, 在此基础上建立热应力模型,最终得到流道边缘^起 与激光工艺参数之问的 系。由于激光功率遵循高斯 分布,流道形状也是轴埘称的,大大简化r模型,节约 了模型运算时间,建 的有限元模 如 2所示。模 型尺寸为180 m×120 m,共有单元21 600个,节点 21 901个 图2激光直写PMMA有限元模型 根据激光直写PMMA原理,当温度升高到降解温 度370℃时,材料吸热但温度不再升高,等到吸收完降 解升华所需的热量后,材料气化掉,我们用“歹匕”单元 技术来模拟材料的气化去除。有多少材料被汽化掉, 对应多少“死”单元,生死单元的分界面即为流道形 状。对于连续激光作用,采用小时『Hj步长逐步施加负 载近似模拟,利用每一一微小时问步计算得到的单元温 度作为单元生死判据,不断循环迭代,最终得到流道形 状及其温度分布,、 2计算结果及分析 聚合物PMMA的参数j“J下:密度p=1 l90 kg/m’, 等压比热c 1 420 J/(kg・K),热传导系数 =0.19 W/(m・K),材料的吸收系数n:0.92,PMMA在受热 时的分解温度为350—380℃,具体依赖于加热速率, Ills、0.4 ms、0.6 ms和0.8 IIIS时的温度场分布。我们 运朋生死单元技术获得了流道外形。图中模 上部单 元缺少,缺少的部分 元为死单元,被激光热源所气 化,剩余的部分为加 后的形状,生死单元的交界部分 为流道形状j采用生死单元不仪能够得到激光烧蚀的 流道外形,还可以得到流道温度场分布,这为分析激光 烧蚀流道热应力热变形打下r基础。从 【11 以看 出,随着照射时问的增加,激光烧蚀的流道外形和深度 都在变化,可见烧蚀流道深度 照射时 有关。 。r 一盯 TEMP(AVG) RSYS=0 SMN=293 SMX=714.12 - { 寸寸 ∞ ∞甘∽ 卜 口 ∞一寸卜 ∞卜 甘∞ { 寸o。 卜_^ 0寸oo ■ 寸ooo0 ∞c。 oo 0 c。 oo cN 0 n 0 0@0卜一 n O0∞。。 卜 卜 n寸 0卜 n寸寸 0 fal 0 2 ms 『bl 0 4 mS TEMP(AVG) TEMP(AVG) RSYS=O RSYS=0 SMN=293 SMN=293 s ∞Ln 一_ 卜 卜 ∞一 寸 卜。o 0 H ∽ O0 O0 寸 0。。 ∽∞ 可 ∽。。 卜 卜 甘 寸 寸。。 。。∞ 0卜 0卜 寸 卜∞0 寸 寸∞∞o。 b- 卜 0 寸 0 0 ∞∽ 甘 卜 ∞ 寸寸 卜【、 (c)0 6 ms 【d】0 8 m 图3 20 W时照射时间分别为0.2 ms、0.4 ms、0.6 ms 和O.8 ms时的温度场 图4是当功率为2O w时特征时问分别为0.2 mS、0.4 HIS、0.6 mS和0.8 IriS时的凸起高度值。从图 中可以看出,随着照射时间的增加,凸起高度也逐渐增 加,fjIJ受热材料的运动变形更为剧烈,与实际的加工结 果是符合的。当照射时间为0.2 Ills时,凸起最大高度 为0.207 m,0.4 Ills时为0.275 m,0.6 HIS H寸为 0.290 m。0.8 Ills时为0.316 n1。这为实际加工中 0 0 控制流道凸起高度提供了依据,对CO 激光烧蚀PM— 卜0. H 甘.0 卜0.∞_【 .0 MA微流道并控制流道表面质量有指导意义。 卜0. 0 0. _[ 0 (AVG) 卜0. 0 0 控制流道凸起高度提供了依据,只需改变激光功率参 卜0. 寸 ,0 数和照射H,卜。 /l ̄,就可以获得大量的流道凸起高度值,卜 .。 忙 对 卜0.∞ ,0 流道凸起高度结果进行分析即可得到流道凸起的变化 规律。可见,合理控制激光的输出功率和照射时间对 减小凸起高度有着重要的作用。 USUM(AVG) RSYS;0 DMX=O 207E一0 6 SMX=O.207E一0.6 USUMRSYS:0 0 ∞∞∞_【.0 DMX=O.275E一0.6 0. _【 _【.0 SMX=O.275E—O.6 0. 寸∞'【.0 一 。越 ,∞ 。 _- 0 出卜0 .0 0 卜0 ∞ o .0 卜0.∞∞_【 .0 3 结语 本文主要研究了CO,激光直写PMMA微流道引 卜0. ∞_【 0 0. _[.0 0 ∞∞H.0 0。 寸∞ .0 (a1 0.2mS 0.b1 0 4 _[ 0 ms USUM(AVG) USUM(AVG)0, .0 RSYS=0 RSYS=O DMX=0.290E一0.6 DMX=O.316E—0。∞ 卜 .O.6 0 SMX=O.290E一0.6 SMX=O 3l6E一0.6 _ 0 0 0 fr.,2出 口 口 0 0 0 0 6 ms (d)0.8ms 图4功率为20 W时照射时间分别为0.2 ms、0.4 ms、 0.6 ms和0.8 ms时流道凸起高度 O 3 O 量o 2 \ 莲o 2 型o 2 魁 耀O・2 O O.1 8 t/ms 图5不同功率不同照射时间对比图 为了研究激光功率与流道凸起高度的关系,通过 改变激光功率,得到了一系列流道凸起高度,当照射时 间分别为0.2 Ills、0.4 mS、0.6 Ills和0.8 ms时,流道凸 起高度的变化如图5所示。流道凸起高度随照射时间 越长而越高,随着激光功率增大而增高。仿真结果为 起表面凸起的现象,并进行了有限元仿真。仿真结果 眶 表明,在流道形成的过程中,热应力作用会在流道边缘 产生凸起,而且流道凸起高度随激光照射时问越长而 越高,随着激光功率增大而增高。所以通过增大激光 。. 寸 。 扫描速度,减小激光照射时问,选择较低的激光功率可 以减小流道凸起高度。~ 参考文献 [1]KLANK H,KUTTER J P,GESCHKE O.CO2-laser mieromachining and back——end processing for rapid production of PMMA—・based microflu— idic systems[J].Lab on a Chip,2002(2):242-246. [2]SNAKENBORG D,KLANK H,KUTrER J P.Microstructure fabrication with a CO2 laser system[J].Journal of Micromechanics and Microengi— neefing,2004(14):182—189. [3]Martin F Jensen,Mikkel Noerholm,et a1.Microstructure fabricalion with a CO2 laser system:characterization and fabrication of cavities pro— duced by raster scanning of the laser beam[J].Lab on a Chip,2003 (3):302—307. [4]C K Cbung,Y C Lin,G R Huang.Bulge formation and improvement of the polymer in CO2 laser micromachining[J].J.Micromeeh.Micro— eng.,2005,55(15):1878一】884. [5]j M Li,C Liu,L Y Zhu.The formation and elimination of polymers bul- ges in CO2 laser microfabrication[J].Journal of Materials Processing Technology,2009,209(10):4814. [6]Chryssolouris G.1aser machining:theor.yand practicef M j.New York: Springer,1991. [7]D Hanselman,B Littlefield.Mastering MATLAB 6[M].北京:清华大 学出版社,2002:310—343. [8]MEUNG JUNG KIM{.Transient evaporative laser-cutting with bound ary element method[J].Applied Mathematical Modelling,2000(25): 25—39. [9]Barbara A Fogarty,Kathleen E lteppert,Theodore J Cory,et a1.Rapid fabrication of poly(din1ethy1si10xane)~based microchip capillary electro— phoresis devices using CO2 laser ablation[J].The Analyst,2005,1 30 (6):924—930. (编辑余捷) (收稿¨期:2o】1—02—13) 文章编号:1 10639 如果您想发表对本文的看法。请将文章编号填入读者意见调查表中的相应位置。 剧童箍. 和稼 20tt年第6期