留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

玻璃钢激光烧蚀碳化致微波传输衰减的数值模拟

刘卫平 吕玉伟 吴丽雄 韦成华 王家伟 韩永超 张爽

刘卫平, 吕玉伟, 吴丽雄, 韦成华, 王家伟, 韩永超, 张爽. 玻璃钢激光烧蚀碳化致微波传输衰减的数值模拟[J]. 红外与激光工程, 2021, 50(12): 20210137. doi: 10.3788/IRLA20210137
引用本文: 刘卫平, 吕玉伟, 吴丽雄, 韦成华, 王家伟, 韩永超, 张爽. 玻璃钢激光烧蚀碳化致微波传输衰减的数值模拟[J]. 红外与激光工程, 2021, 50(12): 20210137. doi: 10.3788/IRLA20210137
Liu Weiping, Lv Yuwei, Wu Lixiong, Wei Chenghua, Wang Jiawei, Han Yongchao, Zhang Shuang. Numerical simulation of ablation and carbonization of GFRP irradiated by laser induced microwave transmission decay[J]. Infrared and Laser Engineering, 2021, 50(12): 20210137. doi: 10.3788/IRLA20210137
Citation: Liu Weiping, Lv Yuwei, Wu Lixiong, Wei Chenghua, Wang Jiawei, Han Yongchao, Zhang Shuang. Numerical simulation of ablation and carbonization of GFRP irradiated by laser induced microwave transmission decay[J]. Infrared and Laser Engineering, 2021, 50(12): 20210137. doi: 10.3788/IRLA20210137

玻璃钢激光烧蚀碳化致微波传输衰减的数值模拟

doi: 10.3788/IRLA20210137
基金项目: 激光与物质相互作用国家重点实验室基金(SKLLIM1701)
详细信息
    作者简介:

    刘卫平,男,副研究员,硕士,主要从事激光辐照效应与机理方面的研究

  • 中图分类号: TB303; O242.1

Numerical simulation of ablation and carbonization of GFRP irradiated by laser induced microwave transmission decay

  • 摘要: 当激光辐照玻璃钢烧蚀碳化至一定程度时,产生的树脂碳产物对微波传输产生衰减作用。针对该现象,开展了数值建模研究,将激光辐照-微波传输衰减效应分解为激光辐照、材料热响应、提取模型表征量、微波传输衰减分析等过程。通过玻璃钢材料的激光耦合特性和表面温度测试,对建立的玻璃钢层合板激光辐照温度场计算模型进行了验证;通过材料体温度分布的时间演进分析,提取了网格单元温度超过阈值温度的持续时间加权和St, Tc、网格单元温度超过阈值温度的持续时间与温度乘积的加权和STt, Tc两个模型表征量,采用单个实验数据标定系数、整体数据点匹配分析方法,对微波传输衰减实验与数值模拟结果进行了分析。分析结果表明,St, TcSTt, Tc更适于表征玻璃钢烧蚀碳化致微波传输衰减效应,温度阈值为873 ℃时计算结果与实验结果最匹配,其R2为0.9956。这说明,通过计算激光辐照玻璃钢温度响应并提取微波传输衰减效应表征量St, Tc,可实现对玻璃钢激光烧蚀碳化致微波传输衰减效应的模拟和预测。
  • 图  1  数值建模分析总体思路框图

    Figure  1.  Idea scheme of numerical modeling analysis

    图  2  玻璃钢的激光耦合特性测试实验布局

    Figure  2.  Experimental setup of laser coupling characteristic measurement of GFRP

    图  3  玻璃钢激光耦合特性测试实验结果

    Figure  3.  Experiment results of laser coupling characteristics test of GFRP

    图  4  激光辐照玻璃钢层合板温度场计算模型计算流程

    Figure  4.  Flow diagram of temperature field calculation model for laser irradiating GFRP laminated plate

    图  5  激光辐照玻璃钢层合板温度场计算模型计算结果与实验结果对比

    Figure  5.  Comparison between calculation results and experiment results of temperature field calculation model for laser irradiating GFRP laminated plate

    图  6  典型的玻璃钢烧蚀碳化形貌

    Figure  6.  Typical ablation and carbonization morphology of GFRP

    图  7  玻璃钢试样的介电常数ε和损耗角正切tanδ测试结果

    Figure  7.  Experiment results of permittivity ε and tangent of loss angle tanδ of GFRP samples

    图  8  微波传输衰减量随激光功率密度变化趋势

    Figure  8.  Tendency of microwave transmission decay varying with laser power density

    图  9  典型的模型计算结果与实验数据匹配结果

    Figure  9.  Typical matching results between model calculating results and experiment data

    图  10  计算结果与实验结果匹配R2随温度阈值的变化曲线

    Figure  10.  Results matching R2 curves of simulation and experiment varying with temperature threshold

  • [1] Penn L, Larsen F. Physicochemial properties of kevlar 49 fiber [J]. Journal of Applied Polymer Science, 1979, 23(1): 59-73. doi:  10.1002/app.1979.070230106
    [2] Yi Shoujun, Fan Zhen, Wu Chao, et al. Catalytic graphitization of furan resin carbon by yttrium [J]. Carbon, 2008, 46(2): 378-380. doi:  10.1016/j.carbon.2007.11.012
    [3] 余维新. 碳材料催化石墨化研究[D]. 长沙: 湖南大学, 2008: 30-32.

    Xu Weixin. Studies on catalytic graphitization of carbon materials[D]. Changsha: Hunan University, 2008: 30-32. (in Chinese)
    [4] 李威龙. 低温各向同性热解炭涂层结构和性能研究[D]. 杭州: 杭州电子科技大学, 2016: 17-21.

    Li Weilong. Research of microstructure and properties of low temperature isotropic pyrocarbon coating[D]. Hangzhou: Hangzhou Dianzi University, 2016: 17-21. (in Chinese)
    [5] Russell N V, Williamson J, Gibbins J R. Structural ordering in high temperature coal chars and the effect on reactivity [J]. Fuel, 1999, 78(7): 803-807. doi:  10.1016/S0016-2361(98)00210-5
    [6] Gao Rui, Wang Chengguo, Li Binpeng, et al. Preparation and microwave absorbing property of polyacrylonitrile based porous disorder carbon [J]. Fiber Reinforced Plastics/Composites, 2012(s1): 111-114. (in Chinese)
    [7] 徐樑, 煤焦微波介电性能的研究[D]. 太原: 太原理工大学, 2015: 25-28.

    Xu Liang. Research on microwave dielectric properties of coal chars[D]. Taiyuan: Taiyuan University of Technology, 2015: 25-28. (in Chinese)
    [8] Wang Jingjing, Liu Haiyu, Xu Liang, et al. Research on dielectric property of biomass and coal during carbonization [J]. Science Technology and Engineering, 2016, 16(14): 181-185. (in Chinese) doi:  10.3969/j.issn.1671-1815.2016.14.033
    [9] Xu Liang, Liu Haiyu, Jin Yan, et al. Variation of microstructure at elevated temperatures and its impact on dielectric properties of coal chars [J]. Journal of Engineering Thermophysics, 2015, 36(5): 1148-1152. (in Chinese)
    [10] Zubkova V, Prezhdo V. Change in electric and dielectric properties of some Australian coals during the processes of pyrolysis [J]. Journal of Analytical and Applied Pyrolysis, 2006, 75(2): 140-149. doi:  10.1016/j.jaap.2005.05.002
    [11] Monsef M P, Ravindran M, McWhinnie W R, et al. Rapid microwave pyrolysis of coal methodology and examination of the residual and volatile phases [J]. Fuel and Energy Abstracts, 1995, 36(2): 20-27.
    [12] Ge Lichao, Zhang Yanwei. Effects of microwave irradiation treatment on physico- chemical characteristics of Chinese low-rank coals [J]. Energy Conversion and Management, 2013, 71(7): 74-91.
    [13] Peng Zhiwei, Hwang Jiann-Yang, Kim Byoung-Gon, et al. Microwave absorption capability of high volatile bituminous coal during pyrolysis [J]. Energy and Fuels, 2012, 26(7): 5146-5151.
    [14] Semak V V, Miller T F. Modeling of laser charring and material removal in fiberglass materials [J]. Journal of Directed Energy, 2006, 2: 5-21.
    [15] Liu Zunfeng, Bai Gang, Huang Yi, et al. Reflection and absorption contributions to the electromagnetic interference shielding of single-walled carbon nanotube/polyurethane composites [J]. Carbon, 2007, 45(4): 821-827. doi:  10.1016/j.carbon.2006.11.020
    [16] Cao Maosheng, Song Weili, Hou Zhiling, et al. The effects of temperature and frequency on the dielectric properties, electromagnetic interference shielding and microwave-absorption of short carbon fiber/silica composites [J]. Carbon, 2010, 8(3): 788-796.
  • [1] 张子浩, 王旭, 黄怡晨, 李福泉, 李俐群, 蔺晓超, 杨诗瑞, 郭鹏.  供油孔流量数值模拟及其飞秒激光加工工艺优化 . 红外与激光工程, 2023, 52(4): 20220454-1-20220454-10. doi: 10.3788/IRLA20220454
    [2] 沈玲君, 宋英雄.  地星大气湍流路径三层传输模型数值模拟 . 红外与激光工程, 2023, 52(11): 20230125-1-20230125-13. doi: 10.3788/IRLA20230125
    [3] 傅杨奥骁, 江涛, 刘庆宗, 丁明松, 李鹏, 董维中, 高铁锁, 许勇.  高超声速飞行器热喷高温燃气红外辐射特性数值模拟 . 红外与激光工程, 2022, 51(10): 20220023-1-20220023-10. doi: 10.3788/IRLA20220023
    [4] 李丹梦, 金伟其, 李力, 鲁啸天, 裘溯.  水下运动目标的水面波纹数值模拟及分析 . 红外与激光工程, 2018, 47(11): 1126004-1126004(8). doi: 10.3788/IRLA201847.1126004
    [5] 李紫慧, 王续跃.  金属复合板激光弯曲过程中翘曲变形数值模拟 . 红外与激光工程, 2018, 47(5): 506004-0506004(7). doi: 10.3788/IRLA201847.0506004
    [6] 李紫慧, 王续跃.  层合板激光往复扫描弯曲中翘曲变形数值模拟 . 红外与激光工程, 2018, 47(12): 1206008-1206008(8). doi: 10.3788/IRLA201847.1206008
    [7] 彭进, 张文洁, 王星星, 郭国全, 张芙蓉.  焊丝熔化填充方式对激光焊接熔池影响的数值模拟 . 红外与激光工程, 2018, 47(3): 306005-0306005(7). doi: 10.3788/IRLA201847.0306005
    [8] 高铁锁, 江涛, 丁明松, 董维中, 刘庆宗.  高超声速拦截弹绕流红外辐射特性数值模拟 . 红外与激光工程, 2017, 46(12): 1204001-1204001(8). doi: 10.3788/IRLA201746.1204001
    [9] 李盾, 宁禹, 吴武明, 孙全, 杜少军.  旋转相位屏的动态大气湍流数值模拟和验证方法 . 红外与激光工程, 2017, 46(12): 1211003-1211003(7). doi: 10.3788/IRLA201746.1211003
    [10] 于学丽, 丁双红, 贾海旭, 辛磊.  主动调Q腔内和频拉曼激光器的数值模拟 . 红外与激光工程, 2017, 46(9): 906001-0906001(7). doi: 10.3788/IRLA201746.0906001
    [11] 刘家琛, 唐鑫, 巨永林.  微型红外探测器组件快速冷却过程数值模拟分析 . 红外与激光工程, 2015, 44(3): 816-820.
    [12] 戴毅斌, 樊玉杰, 李明尧, 刘晓宇, 蒋彭胜, 殷开婷.  多点激光微冲击成形的数值模拟研究 . 红外与激光工程, 2015, 44(S1): 50-56.
    [13] 张金玉, 孟祥兵, 杨正伟, 王冬冬, 陶胜杰.  红外锁相法涂层测厚数值模拟与分析 . 红外与激光工程, 2015, 44(1): 6-11.
    [14] 常金勇, 强希文, 胡月宏, 宗飞, 李志朝, 封双连.  激光传输光束抖动效应的数值模拟 . 红外与激光工程, 2015, 44(S1): 46-49.
    [15] 许兆美, 刘永志, 周建忠, 蒋素琴, 王庆安, 汪通, 洪宗海.  基于生死单元的激光铣削温度场数值模拟与验证 . 红外与激光工程, 2014, 43(6): 1755-1760.
    [16] 陈林, 杨立, 范春利, 吕事桂, 石宏臣.  线性调频激励的红外无损检测及其数值模拟 . 红外与激光工程, 2014, 43(5): 1385-1389.
    [17] 孙杜娟, 胡以华, 李乐.  气溶胶“沉降—扩散”联合动态数值模拟 . 红外与激光工程, 2014, 43(2): 449-453.
    [18] 崔文达, 杜少军.  角谱传输数值模拟中参数优化分析 . 红外与激光工程, 2014, 43(9): 2935-2940.
    [19] 吴慧云, 陈金宝, 李鑫, 黄志松, 黄书海, 赵四清, 孙振海, 徐卸古.  激光中继镜传输的等效菲涅耳数分析 . 红外与激光工程, 2013, 42(12): 3218-3222.
    [20] 唐力铁, 于志闯, 赵乐至, 尹飞, 郭士波, 谈斌.  DF/HF化学激光器喷管总压损失的数值模拟 . 红外与激光工程, 2013, 42(5): 1194-1197.
  • 加载中
图(10)
计量
  • 文章访问数:  350
  • HTML全文浏览量:  79
  • PDF下载量:  28
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-29
  • 修回日期:  2021-04-20
  • 网络出版日期:  2022-01-06
  • 刊出日期:  2021-12-31

玻璃钢激光烧蚀碳化致微波传输衰减的数值模拟

doi: 10.3788/IRLA20210137
    作者简介:

    刘卫平,男,副研究员,硕士,主要从事激光辐照效应与机理方面的研究

基金项目:  激光与物质相互作用国家重点实验室基金(SKLLIM1701)
  • 中图分类号: TB303; O242.1

摘要: 当激光辐照玻璃钢烧蚀碳化至一定程度时,产生的树脂碳产物对微波传输产生衰减作用。针对该现象,开展了数值建模研究,将激光辐照-微波传输衰减效应分解为激光辐照、材料热响应、提取模型表征量、微波传输衰减分析等过程。通过玻璃钢材料的激光耦合特性和表面温度测试,对建立的玻璃钢层合板激光辐照温度场计算模型进行了验证;通过材料体温度分布的时间演进分析,提取了网格单元温度超过阈值温度的持续时间加权和St, Tc、网格单元温度超过阈值温度的持续时间与温度乘积的加权和STt, Tc两个模型表征量,采用单个实验数据标定系数、整体数据点匹配分析方法,对微波传输衰减实验与数值模拟结果进行了分析。分析结果表明,St, TcSTt, Tc更适于表征玻璃钢烧蚀碳化致微波传输衰减效应,温度阈值为873 ℃时计算结果与实验结果最匹配,其R2为0.9956。这说明,通过计算激光辐照玻璃钢温度响应并提取微波传输衰减效应表征量St, Tc,可实现对玻璃钢激光烧蚀碳化致微波传输衰减效应的模拟和预测。

English Abstract

    • 玻璃纤维增强树脂基复合材料(简称玻璃钢,英文缩写GFRP)是一种使用广泛的功能复合材料,具有比质量轻、比强度高,易于成形制造等优点。玻璃钢中的增强纤维为玻璃纤维,树脂基体是一种有机高分子聚合物,均具有优良的介电特性,对微波传输的衰减作用很小,因而,玻璃钢常用于制造微波器件的容器或外形蒙皮,如雷达罩等。

      当玻璃钢受到激光辐照、气动加热等外界热流载荷作用时,树脂基体受热发生热降解,聚合物分子链被打断、缩聚,形成高度无序的热解碳产物[1-3]。热解碳是无定形碳,当其温度达到一定程度后,其微观结构产生变化,局部向石墨层状结构转变,形成石墨微晶[4-5]。当石墨微晶层数和尺寸增加时,石墨层状结构形成的π键电子增加,π键电子在微波作用下形成局域电流,从而将部分微波能量转换为电路的焦耳热,使微波传输产生损耗。在这个过程中,由于无定形碳结构向石墨结构转变存在阈值温度[6-9],热解碳的温度演变过程对石墨微晶产生与转化有重要影响。高瑞等[6]研究了碳化温度对PAN基多孔无定形碳介电性能和吸波性能的影响,发现热处理温度达到800 ℃及以上时,损耗角正切出现了明显增大。V.Zubkova[10]、M.P.Monsef[11]、Ge[12]、Peng[13]等人分别对碳材料热解、碳化等过程中介电性能变化进行研究,发现当温度小于500 °C时,介电性能没有明显变化,当温度更高时,介电常数和损耗角正切都出现明显变化。V.V.Semak等[14]根据玻璃钢热解过程,利用三级热解反应方程对激光烧蚀碳产物产生和移除过程进行了建模,并将热解碳体积分数关联到碳颗粒半径进而关联到射频波衰减表达式,但其模型复杂,参数众多,模型中也未考虑碳产物微观结构的影响。

      文中针对玻璃钢激光烧蚀碳化后树脂碳微观结构变化导致的微波传输衰减效应开展了数值模拟研究,基于玻璃钢温度场时间演变过程,对影响树脂碳微观结构变化过程的相关量进行分析,提取出微波传输衰减的模型表征量,对数值计算模型进行了验证。

    • 针对玻璃纤维/环氧树脂(E-51环氧树脂)复合材料层合板开展了玻璃钢激光烧蚀碳化致微波传输衰减效应的数值建模研究。

      为了模拟玻璃钢激光烧蚀碳化致微波传输衰减效应,将激光辐照-微波传输衰减效应分解为激光辐照耦合、材料热响应、提取表征量、微波传输衰减分析等过程,将激光辐照过程解耦,基于材料温度响应开展分析。

      数值建模分析的总体思路如图1所示,首先,设计实验获取材料对激光的耦合特性,作为温度响应计算模型的输入,构建激光辐照玻璃钢温度场的计算模型;其次,通过实验获得特定功率密度激光辐照下材料的表面温度响应,对计算模型进行校验;随后,采用模型计算得到不同激光功率密度参数下试样的温度响应,获得材料的体温度场,实现从实验测量可获得的点/面温度等有限信息向体温度信息的升维处理;根据材料体温度场的时间演变过程,分析提取出模型计算结果中可用于表征微波传输衰减效应的表征量,对该表征量的合理性进行分析。具体方法是,从玻璃钢烧蚀碳化微波传输衰减的实验结果中选取一个数据点(通常选择绝对值较大的数据点),对表征量进行标定,然后根据标定系数对计算获得的表征量变化曲线进行整体校准,分析校准后的表征量变化曲线与微波传输衰减测量结果间的匹配一致性,判断表征量是否具有合理性。如果校核后的表征量与实验结果一致,则说明该表征量可用于模拟微波传输衰减效应。

      图  1  数值建模分析总体思路框图

      Figure 1.  Idea scheme of numerical modeling analysis

    • 采用积分球对玻璃钢样品激光反射率、透过率进行测试,得到了材料耦合特性变化规律。实验布局如图2所示,激光经过扩束和光束整形,通过光阑、反射率测量积分球后辐照在试样上,试样后放置积分球测量透过的激光。

      图  2  玻璃钢的激光耦合特性测试实验布局

      Figure 2.  Experimental setup of laser coupling characteristic measurement of GFRP

      测量得到试样反射率、透过率随时间变化的曲线后,移除积分球,采用相同实验条件通过热像仪测量试样表面温度,从而得到材料在不同温度下的激光耦合特性曲线,如图3所示。从图3可知,材料对激光的吸收率在459 ℃之前有先下降后上升再下降的变化过程,这和玻璃钢在激光加热作用下的表面状态、树脂整体均匀度变化等有关;在459 ℃后,材料表面热解产生的热解碳产物使材料对激光的吸收由体吸收模式转变为面吸收模式,激光能量主要沉积在材料表面。

      图  3  玻璃钢激光耦合特性测试实验结果

      Figure 3.  Experiment results of laser coupling characteristics test of GFRP

      根据实验现象和材料耦合特性的测试结果,基于商业有限元软件建立了玻璃钢层合板激光辐照温度场有限元计算模型,模型基于傅里叶传热方程,涉及热传导、辐射散热、对流散热过程,激光辐照等效为材料体热源或表面热流。其控制方程组为:

      $$\left\{ {\begin{array}{*{20}{c}} {\begin{array}{*{20}{c}} {\rho c\dfrac{{\partial T}}{{\partial t}} = \nabla \cdot \left( {k\nabla T} \right) + \dddot q} \\ { - k\dfrac{{\partial T}}{{\partial n}} = {q_{{\rm{conv}}}} + {q_{{\rm{rad}}}} - {q_{{\rm{flux}}}}} \end{array}} \\ {{q_{{\rm{conv}}}} = {h_f}\left( {T - {T_0}} \right)} \\ {{q_{{\rm{rad}}}} = \varepsilon \sigma ({T^4} - T_0^4)} \\ {{q_{{\rm{flux}}}} = \alpha {I_0}} \end{array}} \right.$$ (1)

      式中:ρ为密度;c为比热容;T为材料温度;T0为环境温度;k为热导率;$\dddot q$为体生热率;n为表面法向单位矢量;qconv为对流热流;qrad为热辐射热流;qflux为激光等效面热流;hf为对流换热系数;ε为材料发射率;σ为Stefan-Boltzmann常数,其值为5.67×10−8 W/(m2·K4);α为材料对激光的吸收率;I0为入射激光功率密度。

      层合板厚度为2 mm,单层玻璃钢厚度为0.1 mm,层合层数为10层。实验试样长宽尺寸为22.8 mm×10.2 mm。激光光斑为均匀分布的矩形光斑,光斑尺寸为16 mm×10.2 mm,即在材料宽度方向为全覆盖。由于玻璃钢层合板热导率较低,激光辐照时间也较短,热量在层合板面内方向的扩散长度较小,故选取材料宽度方向的中截面为特征面,构建了二维计算模型,材料宽度方向任一截面上的温度分布可认为与中截面上的温度分布一致。计算模型根据玻璃钢的层压层数分别构建各层的几何形态,各层间通过接触换热系数(thermal contact conductance coefficient, TCC)定义层间热接触,TCC随材料热解状态变化。模型计算流程如图4所示。在模型中,采用热解温度判据(459 ℃)判断材料对激光的吸收是否出现体吸收模式到面吸收模式的转变。若未出现转变,模型中材料对激光的吸收以体热源形式进行模拟;否则,材料对激光的吸收以面热流形式进行模拟。

      图  4  激光辐照玻璃钢层合板温度场计算模型计算流程

      Figure 4.  Flow diagram of temperature field calculation model for laser irradiating GFRP laminated plate

      模型中,材料的比热容取1044 J/(kg·K),密度为2 000 kg/m3,热导率1.18 W/(m·K)。

      设计实验对计算模型进行了验证。验证实验采用图2所示的实验布局,用功率密度为75 W/cm2的近红外激光辐照在玻璃钢试样上,采用热像仪测量试样表面的最高温度。模型计算结果与实验结果对比如图5所示。从图中可得,模型计算结果与实验结果基本一致。另外,图中实验曲线在300 ℃时截止,是由于热像仪在选择高温档位时,其温度测量区间下限为300 ℃。

      图  5  激光辐照玻璃钢层合板温度场计算模型计算结果与实验结果对比

      Figure 5.  Comparison between calculation results and experiment results of temperature field calculation model for laser irradiating GFRP laminated plate

    • 采用不同功率密度的激光对玻璃钢试样进行烧蚀,获得不同烧蚀碳化程度的测试试样。典型的玻璃钢烧蚀碳化形貌如图6所示。

      图  6  典型的玻璃钢烧蚀碳化形貌

      Figure 6.  Typical ablation and carbonization morphology of GFRP

      采用矢量网络分析仪和波导管对试样烧蚀碳化后在10 GHz微波下的介电特性进行测试,得到试样介电常数ε和损耗角正切tanδ,结果如图7所示。

      图  7  玻璃钢试样的介电常数ε和损耗角正切tanδ测试结果

      Figure 7.  Experiment results of permittivity ε and tangent of loss angle tanδ of GFRP samples

      根据公式(1)可得到试样的微波传输衰减量[15-16]

      $$\left\{ \begin{gathered} A = 10\log \left(\frac{{{P_0}}}{P}\right) \\ P = {P_0}\exp ( - 2\gamma l) = {P_0}\exp \left( { - 2l(\alpha + \beta i)} \right) \\ \alpha = \frac{{2{\rm{\pi }}}}{{{\lambda _0}}}\sqrt {\frac{{{\varepsilon _r}}}{2}[\sqrt {1 + {{\tan }^2}\delta } - 1]} \\ \beta = \frac{{2{\rm{\pi }}}}{{{\lambda _0}}}\sqrt {\frac{{{\varepsilon _r}}}{2}[\sqrt {1 + {{\tan }^2}\delta } + 1]} \\ \end{gathered} \right.$$ (1)

      式中:P0为进入样品前的微波功率;P为经样品传输后的功率;γ为微波的传播系数,其实部与虚部分别为衰减系数α(单位N·p/m,其中1 N·p/m=8.686 dB/m)与相位系数β(单位rad/m);l为微波在材料中传输距离;εr为介电常数实部;λ0为微波波长。计算得到的试样微波传输衰减量变化趋势如图8所示。可见,微波传输衰减量变化趋势与损耗角正切tanδ的变化相似。这是由于,在公式(1)中,微波传输衰减量与介电常数实部的方根成正比,与损耗角正切近似成正比,而图7中介电常数实部方根的相对变化率较小(以75 W/cm2处的数据为比较基准,最大相对变化率为15%),损耗角正切tanδ相对变化率较大(最大相对变化率为60%),故微波传输衰减量曲线形态主要与损耗角正切曲线形态相关。

      图  8  微波传输衰减量随激光功率密度变化趋势

      Figure 8.  Tendency of microwave transmission decay varying with laser power density

      在环氧树脂热解碳化过程中,温度是影响树脂碳微观结构转变的重要因素。在这个过程中,存在温度阈值Tc,只有材料温度超过该温度阈值时,树脂碳微观结构才发生转变。树脂碳温度越高,微观结构转变速率越快。另外,树脂碳微观结构转变过程是一个与持续时间相关的过程,当树脂碳温度超过Tc后,其微观结构以一定速率从对微波无吸收的无定形碳结构转变为对微波存在吸收的石墨微晶结构。这个过程持续的时间越长,最终转变得到的石墨微晶结构含量越多,对微波传输衰减作用越强。

      根据上述分析,为模拟树脂烧蚀碳化所致微波传输衰减变化规律,预设了2个模型表征量St,TcSTt,Tc来提取整合玻璃钢层合板激光辐照温度场计算模型得到的材料体温度场信息,研究这2个表征量是否可以用于描述微波传输衰减变化趋势。其中,St,Tc表示网格单元温度超过Tc的持续时间加权和:

      $${S_{t,{T_{\rm{c}}}}} = \frac{1}{{N{t_{\rm{0}}}}}\sum\limits_{i = 1}^n {{t_i}({T_i} > {T_c})} $$ (2)

      式中:N为网格单元总数;t0为计算总时间;Ti为网格单元i的温度;ti为第i个网格单元温度超过Tc的总持续时间;n为温度超过Tc的网格单元数。STt,Tc表示网格单元温度超过Tc的温度与持续时间乘积加权和:

      $${S_{Tt,{T_{\rm{c}}}}} = \frac{1}{{{T_0}N{t_{\rm{0}}}}}\sum\limits_{i = 1}^n {\sum\limits_{j = 1}^m {{T_{ij}}} \Delta {t_j}({T_{ij}} > {T_c})} $$ (3)

      式中:T0为所有计算结果中的最高温度;Tij为分析时间步j(其时间步长为$\Delta {t_j}$)中网格单元i的温度。Nt0T0可视为表征量中的权重因子。

      假定St,TcSTt,Tc与微波传输衰减量A间存在比例关系,即

      $$\left\{ \begin{gathered} {A_t} = {c_t}{S_{t,{T_{\rm{c}}}}} \\ {A_{Tt}} = {c_{Tt}}{S_{Tt,{T_{\rm{c}}}}} \\ \end{gathered} \right.$$ (4)

      式中:ctcTt为比例系数。取一个实验数据点对系数ctcTt定标,观察AtATt与其他实验数据点的匹配程度。若AtATt与实验数据A匹配程度高,则说明该表征量具有合理性。

      以玻璃钢烧蚀碳化试样介电特性测试实验中的激光参数作为参考,采用玻璃钢层合板激光辐照温度场计算模型进行温度场计算,处理温度场计算结果时在数据处理程序中取不同的温度阈值Tc对温度场的演进过程数据进行处理,得到St,TcSTt,Tc。其中,温度阈值Tc取值时,依据文献研究结果取值,并覆盖较大的区间,以进行进一步的分析。以微波传输衰减实验数据中的最大值为标定点,分别对St,TcSTt,Tc进行标定,可得到标定系数,从而得到模型计算的材料微波传输衰减量变化曲线。典型的模型计算结果与实验结果对比如图9所示。由图可见,温度阈值为873 ℃时,计算结果与实验结果的匹配程度良好,而温度阈值为800、1000 ℃时,计算结果与实验结果在激光功率密度较低时匹配程度较差,产生微波传输衰减作用的功率密度阈值存在减小/增大的现象。分析出现这种现象的原因认为,若模型中温度阈值Tc取值高于树脂碳实际温度阈值Tc0,则当材料温度$T \in [{T_{{\rm{c}}0}},{T_{\rm{c}}})$时,实验中树脂碳已经产生石墨微晶,对微波传输产生衰减作用,但计算模型中提取的St,TcSTt,Tc却因未满足条件未能反映这种作用,表现为计算结果出现微波传输衰减的功率密度阈值增大,如图9(b)所示。反之,若模型中温度阈值Tc取值低于树脂碳实际温度阈值Tc0,则提取的St,TcSTt,Tc提前反映出微波传输衰减作用,表现为计算结果出现微波传输衰减的功率密度阈值减小,如图9(c)所示。另外,图9(b)(c)中计算结果曲线出现折线的原因尚需进一步分析研究。

      图  9  典型的模型计算结果与实验数据匹配结果

      Figure 9.  Typical matching results between model calculating results and experiment data

      为了量化计算结果与实验结果的匹配一致性,选择决定系数R2为表征量。决定系数R2常用于曲线拟合或回归分析,用于判断拟合或回归模型的优劣,其值越大代表拟合/回归模型越接近实际,其计算式为:

      $${R^2} = 1 - \frac{{\displaystyle\sum\limits_{i = 1}^n {{{\left( {{y_i} - {{\hat y}_i}} \right)}^2}} }}{{\displaystyle\sum\limits_{i = 1}^n {{{\left( {{y_i} - \bar y} \right)}^2}} }}$$ (5)

      式中:yi为实测值;${\hat y_i}$为拟合值;$\bar y$为实测值的平均值。对取不同温度阈值时的R2进行比对,得到如图10所示的结果。

      图  10  计算结果与实验结果匹配R2随温度阈值的变化曲线

      Figure 10.  Results matching R2 curves of simulation and experiment varying with temperature threshold

      可见,计算模型中采用表征量St,TcR2最大为0.995 6,采用STt,TcR2最大为0.9933,故采用表征量St,Tc更匹配实验结果。此时,对应的温度阈值为873 ℃。另外,采用表征量St,Tc得到的R2最大为0.9956,表明此时计算结果与实验结果间的匹配偏差很小。这说明,玻璃钢激光烧蚀碳化致微波传输衰减效应的数值计算模型是合理有效的。

    • (1)构建了玻璃钢激光烧蚀碳化致微波传输衰减数值模型分析总体方案,建立了激光辐照玻璃钢温度场计算模型,设计实验获得了材料的温度相关性耦合特性,并对模型进行了校验,结果表明温度场计算模型与实验结果基本一致。

      (2)采用不同功率密度激光辐照玻璃钢样品,获得了不同烧蚀碳化程度玻璃钢的介电常数实部、损耗角正切的变化曲线,从而计算得到了微波传输衰减量随激光功率密度的变化规律。

      (3)基于玻璃钢温度场演变过程,从数值模型计算结果中提取了网格单元温度超过阈值温度的持续时间加权和、网格单元温度超过阈值温度的持续时间与温度乘积加权和两个表征量进行计算结果与实验数据匹配性分析。结果表明,前者更适合作为烧蚀碳化致微波传输衰减效应的模型表征量,且碳化温度阈值为873 ℃时,可以很好地计算预测玻璃钢烧蚀碳化后的微波传输衰减量。

参考文献 (16)

目录

    /

    返回文章
    返回