Skin damage in mice induced by different power densities of 1064 nm laser
-
摘要: 通过实验和理论分析研究1064 nm激光不同输出功率对小鼠皮肤的热损伤规律。利用皮肤镜与光学相干断层成像（OCT）评估小鼠皮肤组织热损伤程度；参考Arrhenius热损伤方程进行理论分析，并与实验结果对比。结果显示，当激光连续辐照时间为400 ms时，激光输出功率密度小于958 W/cm2时，激光辐照处泛红；激光输出功率密度为958～1160 W/cm2时，损伤呈白色水疱状；激光输出功率密度为1160～1370 W/cm2时，损伤呈浅坑状焦黄斑，损伤斑周围伴一圈鼓起的白色皮肤水疱；激光的功率密度在1370～2190 W/cm2，损伤呈红色坑状斑，损伤斑周围伴黑黄色焦痂。
-
关键词:
- 1064 nm激光 /
- 皮肤热损伤 /
- Arrhenius热损伤方程 /
- 热损伤模型 /
- 热损伤分级
Abstract: The thermal damage of mouse skin under different output power of 1064 nm laser was studied by experiment and theoretical analysis. The degree of thermal injury of mouse skin tissue was evaluated by dermoscope and optical coherence tomography (OCT), and the theoretical analysis was made according to Arrhenius thermal damage equation, and the results were compared with the experimental results. The results show that when the laser’s continuous irradiation time is 400 ms and the laser output power density is less than 958 W/cm2, there is no obvious damage; when the laser output power density is 958−1160 W/cm2, the damage is white blister-like; when the laser output power density is 1160−1370 W/cm2, the damage is a shallow pit-like macula with a circle of bulging white skin blisters around the damage spot; when the laser power density is at 1370−2190 W/cm2, the damage is a red pit-like spot, and there is black-yellow eschar around the injury spot. -
现代光学系统对光学元件表面质量要求越来越高，射流抛光技术可以获得高质量、超光滑的光学元件。1988年，荷兰Delft大学的O.W.Fahnle等人首先提出把射流抛光技术应用到光学元件表面上，并通过实验证明射流抛光技术应用于加工精密的光学元件是可行的[1]。射流抛光技术是近几年发展起来的一种新型的光学加工技术，中国科学院光电技术研究所的施春燕等通过数值分析和模拟，发现速度分布曲线与材料去除曲线相对应，并对喷射距离和冲击角度对材料去除面型的影响做了研究[2-3]。四川大学李秀龙等利用单颗粒冲击去除模型，分析了在塑性加工条件下射流冲击去除效应，并对射流中粒子与元件发生塑性接触的临界速度进行了推导，引入了塑性转入脆性加工的临界速度，分析了粒子速度和半径对冲击去除量的影响[4]。由于掺有一定质量分数磨粒的抛光液的去除效果和效率会更好，但抛光时磨粒的不均匀性可能会导致去除量有所变化。在理想状态分析下认为射流抛光中颗粒尺度分布是均匀的，而在实际的工艺中，冲击到壁面上的颗粒大小不可能完全一致。目前鲜见相关文献考虑到抛光颗粒冲击到石英玻璃壁面上，颗粒分布的不均匀性对去除特性的影响。本文从含有SiC磨粒的抛光液出发，理论研究了在理想状态下(抛光颗粒均匀分布)改变颗粒直径和质量分数时的冲击去除分布的情况，并通过改变抛光粉颗粒的尺度和质量分数来分析颗粒分布的不均匀性对去除量的影响。
1. 射流去除理论
以射流抛光中最常用的单喷嘴结构为研究对象[5-6]，如图 1所示: 喷嘴直径1 mm，入口压强设置为0.8 MPa，网格大小为0.5。射流作用下，抛光粉颗粒作用到材料壁面上时的冲击去除分布为[7]
f(x)=a×K1mpu0exp[−0.693(xb)2]r×(ρρp)2H3/2E−1/2K−2c (1) 式中：a是单位时间内冲击到材料表面的抛光颗粒个数；K1是与材料特性、压力、速度等无关的比例常数；ρ为光学材料的密度；ρp为磨粒的密度；E是材料弹性模量；H为工件材料的硬度；Kc为工件材料的断裂韧度；b为壁面压强降为正对喷头位置压强一半时的位置；x为与壁面中心点的距离；u0为射流喷嘴出口处的速度；r为磨粒半径；mp为颗粒的质量。由此可以计算得到射流抛光颗粒直径改变时，去除量的冲击去除分布曲线。
由于单位时间内作用在材料上的抛光颗粒个数对去除量具有直接影响，因此，考虑抛光液质量分数变化时，材料的冲击去除分布可写为[8]
fR1(x)=K1u0ρscstρpπ2r4[5.76p00.86hd−2.5−pk](ρρp)1/2H3/2E−1/2Kcexp[−0.693(xb)2] (2) 式中：ρs为抛光液密度；p0为喷嘴出口压力；pk为能实现材料磨损的临界压力；c为抛光液质量分数；h为喷射距离；d为喷嘴口径；s是壁面面积；t是冲击时间。
2. 理想状态下材料的去除特性
依据射流抛光的去除机理，即磨粒通过在垂直压力作用下，对元件壁面进行反复冲击和剪切来实现材料的去除，可知有多种因素影响冲击去除量，包括射流速度、喷射距离、压强、冲击角度等[9]。因其他因素在文献[8]中已经加以分析，本文重点探讨抛光液质量分数和抛光颗粒直径对去除量的影响。抛光液采用碳化硅与水按一定质量比混合而成，以石英玻璃作为分析对象，其参数列表如表 1所示[4]，其中J代表比热容，ω代表热导率。抛光颗粒直径选定为20 μm，射流速度选定为20 m/s，抛光时间为100 s。
表 1 石英玻璃参数Table 1. Quartz glass parametersmaterial H/GPa E/GPa J/(kg·K) ρ/(g·cm-3) ω/(m·K) SiC - 450 1 266.93 3.22 16.7 fused silicon 7.1 73 670 2.2 1.4 由于抛光液质量分数对射流抛光中磨粒发挥的冲击作用和在壁面上流动时的剪切作用有着重要影响，因此，我们分析了在理想状态下(即冲击到壁面上的颗粒是均匀分布的)质量分数变化对冲击去除分布的影响。这里质量分数依次选取1%，4%，5%，6%。
由公式(2)模拟计算得到图 2，由此看出在抛光颗粒直径和冲击速度恒定时，不同质量分数的抛光粉颗粒的冲击去除分布具有相似性，服从高斯分布，且随着抛光液质量分数的增大，冲击去除相应增大。这是因为质量分数增大时，单位时间内作用到材料壁面上的颗粒个数增多，颗粒对壁面的撞击作用与撞击频率随之增强，冲击去除量就会变大。
图 2 冲击速度为20 m/s时, 不同质量分数的冲击去除分布Figure 2. Impinging removal distribution of different concentrations at impinging velocity of 20 m/s由抛光机理可知，抛光时颗粒的直径对冲击效果有着直接影响，因此我们进一步分析在一定质量分数下，不同的抛光颗粒直径对冲击去除的影响。这里颗粒直径分别选取10，20，30，40 μm，射流冲击速度为20 m/s，抛光液质量分数为5%，利用FLUENT软件仿真得到壁面冲蚀云图, 如图 3所示。
图 3 不同抛光颗粒直径下的冲蚀云图Figure 3. Erosion nephogram for different polishing particle diameter从图 3可以看出，纵坐标颜色变化代表颗粒冲击到壁面上，去除量的变化。随着颗粒直径的增大，冲蚀效果变得越好，越靠近中心位置对应的去除量越大，冲击去除量随着远离中心位置在减小，而图 3(a)是刚开始流场不稳定造成的结果。从冲蚀云图只能看出不同抛光颗粒直径下去除量分布的差异，而无法看出其差异大小，针对这一问题，我们画出了对应的二维分布图。在抛光液质量分数为5%，颗粒直径分别为10，20，30，40 μm时，每10 s作用到壁面上的颗粒数分别为13 950，3487，1550，872，对应公式(1)中的a值，由公式(1)模拟计算得到图 4。图 4表明：抛光颗粒冲击到石英壁面上时，不同粒径的粒子冲击去除分布具有相似性，冲击去除量随着颗粒直径的增大而减小。这是因为在质量分数和冲击速度为定值时，颗粒直径越小，恒定尺寸的喷嘴喷出的粒子数越多，作用到壁面上的颗粒数就会越多，其同工件壁面的碰撞切削作用随之增强，材料的去除量也就变大。
图 4 冲击速度为20 m/s时，不同抛光颗粒直径下的冲击去除分布Figure 4. Impinging removal distribution of multiple particle diameters at impinging velocity of 20 m/s通过图 4可以看出，抛光粉颗粒的大小会直接影响到材料的去除量，因此如果抛光过程中有不同尺度的抛光粉同时参与作用，材料的去除量有明显的增加趋势。
3. 抛光颗粒尺度均匀性对去除特性的影响
从理论上说，如果是单一尺度的抛光颗粒作用到材料表面，在确定时间内材料的去除量是确定的。而在实际的抛光中选择的抛光粉颗粒直径不可能完全相同，因此可将颗粒直径看成一个有一定波动大小的随机变量。由中心极限定理可知，在颗粒数足够多的情况下，其尺度变化满足随机高斯分布，其方差大小对应了颗粒分布的不均匀性，通过模拟计算不同高斯分布方差σ2下的去除量。
3.1 颗粒直径分布不均匀对去除量的影响
为分析不均匀的颗粒直径的冲击去除特性，把其他物理量设为定值，此时随机变化的不同抛光颗粒直径(即颗粒的不均匀性)满足的冲击去除分布公式应为
f(x)=N∑i=1K1mpu0exp[−0.693(xb)2]ri(ρρp)2H3/2E−1/2K−2c (3) 公式(3)中不均匀分布的半径ri是在公式(1)的基础上变换而成，公式(1)是恒定的抛光颗粒直径(即颗粒是均匀分布的)距离材料中心位置不同时对应的冲击去除分布。
当大量颗粒作用在抛光过程中时，式(3)中半径满足随机高斯分布，其概率密度函数为
f(r)=1√2πσe−(r−ˉr)22σ2 (4) 式中：σ2为方差，表示颗粒间尺度的差异大小；r为颗粒平均半径。
计算了不同颗粒直径情况下，抛光颗粒分布的不均匀对材料的去除特性。首先根据抛光颗粒尺度的随机分布规律生成一个满足公式(3)的随机数组，用其表征颗粒半径r的方差，其次利用公式(3)模拟得到图 5，其中横坐标表示方差，纵坐标表示冲击到材料位置x=0.001 5 m处对应的冲击去除量。
图 5 颗粒半径均值不同时，颗粒分布不均匀性与去除量的关系Figure 5. Relationship between particle distribution inhomogeneity and removal amoun for different average particle radius从图 5看出，当抛光颗粒的直径不再绝对均匀时，冲击后材料的去除量不是一个稳定的值，而是出现了明显的随机波动。由波动图形看出，当抛光颗粒分布较为均匀时，高斯分布的方差较小，对应得到的去除量波动小。随着抛光颗粒平均直径变化的增大，即方差增大，此时材料去除量的波动性也变得越来越大。由数值看出，随着抛光颗粒平均直径的增大，相对应的去除量在逐渐减小，通过计算每个图的去除量增值，可以得到其越来越小，增加的平均数值为0.46×10-4 m。以图 5(b)为例，选取方差为0，1×10-12，2×10-12和3×10-12时，计算材料去除量的增加量分别为0.090 3×10-4，0.251 2×10-4和0.114 8×10-4 m，由此可见波动变化明显。另外，图 5(a)和图 5(b)显示的去除分布随着颗粒分布方差的增大，去除量有增加趋势，而图 5(c)和图 5(d)去除量增加并不明显，这是因为在方差一定的情况下(即所代入的平均半径误差量n∑i=1(ri−ˉr)2相同)，随着r的增大，理论上半径误差量与r的相对误差量越小，因此对去除量的影响也越小，即去除量分布的波动变化也就越小。
3.2 不同质量分数下SiC颗粒不均匀性对去除量的影响
不仅颗粒直径会对去除量有影响，质量分数也是其影响因素，因此还要考虑到颗粒不均匀性在不同质量分数下的去除量变化情况。图 6是不同质量分数下的抛光粉颗粒分布的不均匀性与x=0.001 5 m位置时去除量的关系。抛光液质量分数依次选取1%，4%，5%，6%，颗粒直径设为20 μm，射流速度为20 m/s。
图 6 抛光液质量分数不同时，颗粒分布的不均匀性与去除量的关系Figure 6. Relationship between particle distribution inhomogeneity and removal amount for different concentration of polishing liquid从计算结果可知，当颗粒直径分布确定时，质量分数变化也会对材料的去除特性有明显的影响。随着方差的增大，去除量的波动和数值也在变大，且随着质量分数的增大，去除量的增加数值也随之增大，平均增加了1.322×10-4 m。
通过比较图 6和图 5的图形波动大小，得到在质量分数恒定颗粒直径改变时，相邻方差对应的颗粒去除量的波动变化较大，上下抖动较剧烈。而相较之质量分数改变时，由于抛光粉颗粒直径的分布是一样的，因此具有相近的抖动变化。但是由于质量分数变化直接导致了参与抛光作用的颗粒数量不同，因此去除量受同一粒径不同质量分数下抛光颗粒分布不均匀性的影响较大。
4. 结论
为了分析颗粒分布的不均匀性对冲击去除量的影响，本文从抛光颗粒冲击到石英玻璃壁面上出发，根据冲击去除分布曲线，研究了理想状态下抛光粉颗粒直径和质量分数对材料去除特性的影响。计算结果表明：冲击去除随着直径的增大而减小，随着质量分数的增大而增大。考虑到实际加工中掺杂有不同尺度的颗粒，颗粒直径分布的不均匀性和质量分数改变下的颗粒分布对去除量造成一定的影响，可以得出：随着方差的增大，去除量的波动也很大。在颗粒直径一定时，随着质量分数的增大，方差的增大，去除量的波动变化也增大，去除量的变化数值增大。颗粒直径和质量分数分布的不均匀性都会对去除量造成影响，从波动大小来看，粒径分布不均匀对去除量波动影响较大；从数值增量看出，质量分数改变时的颗粒不均匀性对去除量数值影响较大。
-
图 1 激光辐照小鼠皮肤实验的示意图
Figure 1. Experimental diagram of laser irradiation on the skin of mouse
图 2 1064 nm激光辐照400 ms后小鼠皮肤创面的皮肤镜影像
Figure 2. Dermoscope images of injury in mice skin at different duration of 1064 nm laser irradiation (400 ms after irradiation; scale bar of dermoscope images: 1 mm)
图 3 不同功率密度激光辐照对小鼠皮肤损伤面积的影响
Figure 3. Effect of different power density laser irradiation on the skin lesion area of mice
图 4 1064 nm激光辐照损伤后OCT结果观察（辐照400 ms，标尺：1 mm）
Figure 4. OCT images of injury in mice skin at different duration of 1064 nm laser irradiation（400 ms after irradiation; scale bar of OCT images：1 mm）
图 5 不同功率密度激光辐照对小鼠皮肤损伤深度的影响
Figure 5. Effect of laser irradiation with different power density on the depth of the skin damage in mice
图 6 小鼠经1064 nm激光照射后皮肤创面的病理形态学变化（HE染色，×100）
Figure 6. Pathomorphological changes of skin wounds in mice irradiated by 1064nm laser (hematoxylin-eosin staining ×100)
表 1 1064 nm激光致生物皮肤组织皮肤损伤分级
Table 1. Grading of skin lesions in biological skin tissue by 1064 nm laser
laser irradiation power
density/(W·cm−2)damage appearance ＜958 firstdegree injury Redness at laser irradiation, recovered after some time. 958～1160 seconddegree injury White blisters, deformed and bulging (crumpled) skin around the injury. 1160～1370 thirddegree injury Vaporization of small amounts (or few lesions) on the superficial layer of the skin, blanching skin blisters with a bulging circle around the injury spot, and a shallow pit like focal macule. 1370～2190 fourthdegree injury Through the skin up to the fat layer, superficial vaporization was severe, splashing oil droplets, and around the injury with dark yellow eschar as red pit like spots. 
下载: 导出CSV
表 2 皮肤热损伤程度与激光辐照功率对应关系
Table 2. The corresponding relationship between the degree of skin thermal injury and laser irradiation power
laser irradiation power density/(W·cm−2) Ω 374 1.80×10−5 561 3.00×10−4 762 1.9 958 14 1163 27 1369 111 1554 1200 1738 1.90×108 1934 5.50×1014
下载: 导出CSV
-
[1] 激光军事应用[J]. 中国光学与应用光学文摘, 2001, 15(5): 31-33Laser military application[J]. Chinese Optics, 2001, 15(5): 31-33 [2] 张延令. 激光在外科中的应用[J]. 国外医学. 外科学分册, 1984(3):143-146Zhang Yanling. Application of laser in surgery[J]. Foreign Medical Sciences. Surgery Section, 1984(3): 143-146 [3] 宋秋伦, 洪伟, 王倩. 逍遥丸联合调Q Nd: YAG 1064 nm激光治疗黄褐斑疗效分析[J]. 中国美容医学, 2022, 31(4):114-116Song Qiulun, Hong Wei, Wang Qian. Effect of Xiaoyao pill combined with Nd: YAG 1 064 nm laser on chloasma[J]. Chinese Journal of Aesthetic Medicine, 2022, 31(4): 114-116 [4] 左娜, 陶凯. 1064 nm Nd: YAG激光联合水光注射治疗面部皮肤老化疗效观察[J]. 中国美容医学, 2022, 31(5):23-25Zuo Na, Tao Kai. Efficacy of 1064 nm Nd: YAG laser combined with hyaluronic acid injection in the treatment of facial rejuvenation[J]. Chinese Journal of Aesthetic Medicine, 2022, 31(5): 23-25 [5] 施添霖. 积雪苷联合长脉宽1064 nm激光治疗瘢痕疙瘩的疗效观察[D]. 杭州: 浙江中医药大学, 2019Shi Tianlin. The therapeutic effect of asiaticoside ointment combined with long pulse 1064 nm laser on keloid[D]. Hangzhou: Zhejiang Chinese Medical University, 2019 [6] 冯微, 冯啸, 陈燕, 等. 长脉宽1064 nm激光在整形美容应用中的进展[J]. 中国美容整形外科杂志, 2021, 32(11):696-698 doi: 10.3969/j.issn.1673-7040.2021.11.016Feng Wei, Feng Xiao, Chen Yan, et al. Development of a 1064 nm long pulse width laser for cosmetic and cosmetic applications[J]. Chinese Journal of Aesthetic and Plastic Surgery, 2021, 32(11): 696-698 doi: 10.3969/j.issn.1673-7040.2021.11.016 [7] 李忠明, 骆清明. 激光对生物组织热和热致机械损伤的物理分析[J]. 激光生物学报, 2000, 9(2):81-84Li Zhongming, Luo Qingming. Physical analysis on tissues’ heat and heatinduced injury for laser irradiation[J]. Acta Laser Biology Sinica, 2000, 9(2): 81-84 [8] 王振华, 单健, 陈铀, 等. 激光对生物组织的热损伤机理研究及临床应用探讨[J]. 激光杂志, 2003, 24(1):60-62Wang Zhenhua, Shan Jian, Chen You, et al. A study on thermal damage of laser tissue and clinical application[J]. Laser Journal, 2003, 24(1): 60-62 [9] 吕晨阳, 战仁军. 激光-生物组织光热效应模型综述[J]. 激光杂志, 2021, 42(1):17-24Lv Chenyang, Zhan Renjun. Model review of laser-biological tissue photo-thermal effect[J]. Laser Journal, 2021, 42(1): 17-24 [10] 陈安宁, 吴家辉, 李希靖. 激光对生物组织热损伤的计算机模拟[J]. 光电子·激光, 2001, 12(2):212-214Chen Anning, Wu Jiahui, Li Xijing. Computer simulation of thermal damage to biological tissue during laser irradiation[J]. Journal of Optoelectronics·Laser, 2001, 12(2): 212-214 [11] 朱光明, 朱丹, 骆清铭, 等. 激光热疗中生物组织的光热特性响应以及动态热损伤研究[J]. 中国生物医学工程学报, 2004, 23(2):157-162Zhu Guangming, Zhu Dan, Luo Qingming, et al. Study on the optical-thermal response and dynamic thermal damage of bio-tissue during laser thermotherapy[J]. Chinese Journal of Biomedical Engineering, 2004, 23(2): 157-162 [12] 吴秀山, 侯宇, 李希靖. 强激光作用下生物组织热损伤问题的研究[J]. 中国计量学院学报, 2002, 13(2):135-138Wu Xiushan, Hou Yu, Li Xijing. The research on thermal damage in biological tissues induced by high-irradiance-laser[J]. Journal of China Jiliang University, 2002, 13(2): 135-138 [13] 周巡, 马琼, 刘智搏, 等. 1064 nm激光不同辐照时间对小鼠皮肤热损伤的实验与理论研究[J]. 强激光与粒子束, 2022, 34:011012 doi: 10.11884/HPLPB202234.210338Zhou Xun, Ma Qiong, Liu Zhibo, et al. 1064 nm laser induced thermal injure in mice skin with different laser duration[J]. High Power Laser and Particle Beams, 2022, 34: 011012 doi: 10.11884/HPLPB202234.210338 [14] Kagan R J, Peck M D, Ahrenholz D H, et al. Surgical management of the burn wound and use of skin substitutes: an expert panel white paper[J]. Journal of Burn Care & Research, 2013, 34(2): e60-e79. [15] Ganguly M, Mitra K. A modeling study to analyze thermal and mechanical effects of pulsed laser irradiation on tissues[J]. Computational Thermal Sciences: an International Journal, 2015, 7(5/6): 459-465. -
下载:
点击查看大图
计量
- 文章访问数: 1158
- HTML全文浏览量: 378
- PDF下载量: 130
- 被引次数: 0
