编织复合材料低速冲击损伤超声热成像检测

李胤, 宋远佳, 江海军, 陈飞, 张凯

李胤, 宋远佳, 江海军, 陈飞, 张凯. 编织复合材料低速冲击损伤超声热成像检测[J]. 红外技术, 2023, 45(8): 876-883.
引用本文: 李胤, 宋远佳, 江海军, 陈飞, 张凯. 编织复合材料低速冲击损伤超声热成像检测[J]. 红外技术, 2023, 45(8): 876-883.
LI Yin, SONG Yuanjia, JIANG Haijun, CHEN Fei, ZHANG Kai. Low-velocity Impact Damage Detection of Woven Composites Based on Ultrasonic Infrared Thermography[J]. Infrared Technology , 2023, 45(8): 876-883.
Citation: LI Yin, SONG Yuanjia, JIANG Haijun, CHEN Fei, ZHANG Kai. Low-velocity Impact Damage Detection of Woven Composites Based on Ultrasonic Infrared Thermography[J]. Infrared Technology , 2023, 45(8): 876-883.

编织复合材料低速冲击损伤超声热成像检测

基金项目: 

国家自然科学基金 52005495

中国空气动力研究与发展中心基础和前沿技术研究基金 PJD20200223

详细信息
    作者简介:

    李胤(1990-),男,博士,工程师,主要从事结构检测与性能评估工作。E-mail:DIYLLLLY@163.com

    通讯作者:

    江海军(1988-),男,硕士,高级工程师,主要从事红外无损检测技术研发工作。E-mail:navy63@163.com

  • 中图分类号: V45

Low-velocity Impact Damage Detection of Woven Composites Based on Ultrasonic Infrared Thermography

  • 摘要: 编织复合材料低速冲击损伤主要为内部的分层损伤,采用目视检测无法有效检测损伤,损伤使得结构承载能力严重降低,威胁编织复合材料构件的安全使用。本文使用超声红外热成像技术对编织复合材料低速冲击损伤进行无损检测研究,使用10 J、20 J、30 J、40 J、50 J的冲击能量制作了5个试件。对超声激励过程的温升曲线、空间温度曲线进行了详细分析;对比不同冲击能量试件发现,低速冲击下损伤主要是内部损伤,冲击能量越大,损伤区域越大,且损伤具有延展性。采用曲线分类算法对损伤区域进行了定量识别,发现编织复合材料损伤面积和冲击能量成线性关系。
    Abstract: Low-velocity impact damage of woven composites is mainly caused by internal delamination damage, which cannot be effectively detected by visual inspection; this seriously reduces the structural load-bearing capacity and threatens the safe use of the compiled composite components. In this study, ultrasonic infrared thermography was used to perform nondestructive testing of the low-velocity impact damage of woven composites, and five specimens were produced using impact energies of 10, 20, 30, 40, and 50 J. The temperature increase and space temperature curves of the ultrasonic excitation process were analyzed. By comparing different impact energy specimens, it was found that the damage under low-speed impact was mainly internal, and the larger the impact energy, the larger the damage area. Moreover, the damage was ductile. The damage area was identified quantitatively using a curve classification algorithm. It was found that the damage area of woven composites and the impact energy were linearly correlated.
  • 星敏感器是当前广泛应用的天体敏感器[1-2],其工作环境必然会受到安装面温度变化等影响,光机结构的热变形会导致镜片面型变化,从而影响光学系统成像质量下降[3-4],会对弥散斑产生较大的影响,为了保证系统的成像质量,需要对结构完成光机热集成分析,分析环境温度对镜头的影响[5]

    随着航天事业的发展,星敏感器已经广泛应用于多种场合[6],由于我国星敏感器研究起步较晚,国内对长焦距大口径星敏感器的研究相对较少,孟祥月[7]等研制了焦距50 mm,入瞳直径40 mm的星敏感器。孙东起[8]等人研制了一种焦距200 mm,入瞳直径125 mm的双高斯光学系统的长焦距星敏感器。伍雁雄[9]等研制了焦距200 mm,入瞳直径100 mm的高精度星敏感器。

    本文设计了一种大口径热不敏星敏感器,光学系统焦距900 mm,入瞳直径200 mm,光谱范围450~750 nm,通过光机热集成分析方法对系统进行热分析,通过将Nastran计算的主次镜表面节点刚体位移代入Sigfit光机热耦合软件进行Zernike多项式拟合,再将主次镜表面Zernike系数导入Zemax光学设计软件中,分析了由于温度变化导致的光机结构刚体位移等变化。

    光学系统参数:焦距范围为900 mm,入瞳直径D≥200 mm,光谱范围为470~900 nm,在热不敏光学系统安装面温度为20℃±5℃时,其光轴偏角优于1″,0.8视场下80%能量集中在9.2~18.4 μm之间。光学系统结构如图 1所示,0.8视场下各波段圈入能量曲线如图 2所示,可以看出满足80%能量集中时,弥散斑直径满足指标要求。

    图  1  光学系统结构图
    Figure  1.  Optical system structure diagram
    图  2  0.8视场各波段弥散斑包围能量曲线
    Figure  2.  Surrounding energy curves of scattered spots in various bands of 0.8 field of view

    本系统采用改进型卡式系统,为保证主镜和后接透镜组的同轴度,选用中心固定形式,主镜材料选用微晶玻璃,为达到热不敏效果,减少温度变化对结构的影响,主镜轴材料应选用与主镜材料热膨胀系数相近的殷钢,主镜通过胶层与主镜轴固定连接,主镜轴作为整个系统的连接构件,具有一定的刚性,而胶层的柔性能够很好的减少重力、温度等对主镜产生的变形影响,主镜结构如图 3所示。

    图  3  主镜结构
    Figure  3.  Main mirror structure

    次镜是非常敏感的光学构件,微小的变化都会带来很大影响,并且支架的大小直接影响光学系统的中心遮挡大小,为保证结构稳定、中心遮挡小以及减小加工难度等原因采用三片殷钢片连接主次镜,能够有效减少温度等因素引起的主次镜间距的变化,支撑结构如图 4所示。

    图  4  主次镜支撑结构
    Figure  4.  Main and secondary mirror support structure

    透镜组通过压圈固定方式保证镜片间间距,镜筒材料采用A704能够减轻结构质量,并且在后端机械结构上留有两个接口方便后续探测器接入,系统整体结构如图 5所示。

    图  5  整体结构模型图
    Figure  5.  Overall structural model diagram

    在本系统中,主次镜结构的稳定性对成像质量的影响最大,本次分析只对主次镜结构进行仿真,分析目的是验证主次镜结构在20±5℃范围内是否满足光学系统设计指标要求。

    通过MSC.Patran建立模型如图 6所示,整个模型采用手工划分网格的方法,控制网格疏密,使得计算结果更加精确,模型主要六面体单元及少量的五面体建模,共有单元数12172个,节点数18707个,结构有限元建模计算中主次镜及支撑结构的材料及其属性参数如表 1所示。

    图  6  有限元模型图
    Figure  6.  Finite element model diagram
    表  1  选用材料属性参数
    Table  1.  Selected material property parameters
    Material Elasticity modulus Ea/MPa Poisson ratio μ Density ρ/(103 kg/m3) CTE α/
    (10-6mm/℃)
    Invar 141000 0.25 8.1 0.2
    TC4 114000 0.29 4.4 8.9
    Microcrystalline glass 90600 0.24 2.53 0.5
    D04 RTV 850 0.40 1.15 236
    下载: 导出CSV 
    | 显示表格

    按照指标要求的环境温度25℃,对主次镜模型施加温度载荷,利用Nastran软件计算得到刚体位移结果,主次镜刚体位移云图如图 7所示,可以看出主镜最大轴向位移为0.228 μm,次镜最大轴向位移为0.986 μm,目前来看热变形结果还在可控范围内。

    图  7  整体位移分布
    Figure  7.  Overall displacement distribution

    利用光机热耦合工具Sigfit输入系统主次镜的曲率半径、主次镜表面节点位置数据、热变形后主次镜表面节点变化数据等进行拟合。温度为25℃时,Sigfit拟合得到的Zernike多项式系数[10]表 2所示。

    表  2  Zernike系数
    Table  2.  Zernike coefficient
    Serial number Expression Value (The main mirror) Value(The secondary mirror)
    1 1 1.53E-05 7.30E-06
    2 ρcosθ 6.47E-10 1.15E-08
    3 ρsinθ 4.36E-10 1.09E-11
    4 2ρ2-1 -1.07E-04 1.26E-05
    5 ρ2cos2θ -8.13E-08 -1.41E-10
    6 ρ2sin2θ 4.33E-08 2.76E-10
    7 (3ρ2-2ρ)cos2θ 3.62E-09 1.74E-11
    8 (3ρ2-2ρ)sin2θ 5.39E-09 -1.85E-08
    9 6ρ4-6ρ2+1 3.63E-06 -3.7E-08
    下载: 导出CSV 
    | 显示表格

    将主次镜的Zernike多项式系数导入Zemax光学设计软件中,即可得到系统弥散斑直径以及光轴的变化,图 8给出了在环境温度25℃,0.8°视场下各波段的圈入能量曲线图,由图中信息可知,各波段80%能量弥散斑直径集中在9.2~18.4 μm之间,与图 2对比可知在温度的影响下,各波段的弥散斑直径也会增大。同时由图 9得到波前RMS(Root-Mean-Square)值为0.035λ<1/12λ,成像质量良好,调用评价函数RAID指令,在0°视场入射光线与像面法线夹角可以近似为光轴偏角约为0.033″优先于1″。

    图  8  0.8视场圈入能量曲线
    Figure  8.  0.8 field of view encirclement energy curves
    图  9  波前图
    Figure  9.  Wavefront diagram

    为检验光机热集成分析的准确性以及光机设计的合理性,设置实验室20±5℃的温度条件下,进行光学系统主镜、次镜以及透镜组系统装调,主镜及透镜组利用三坐标进行检测装调,保证其位置精度,然后利用干涉仪进行次镜的装调工作,系统整体装调结构如图 10所示。

    图  10  系统整体装调结构图
    Figure  10.  Overall system assembly and adjustment structure diagram

    在实验室室温25℃下,系统装调后的轴上视场波像差如图 11所示,RMS值为0.08λ,所测得RMS值与有限元分析结果相差很小,分析实例验证了本系统分析方法的有效性。

    图  11  0视场波像差
    Figure  11.  Zero field of view wave aberration

    本次测试温度环境分别设为15℃、20℃、25℃,采用平行光管照射,镜头放置在精密旋转的调整台上,通过对镜头的成像光斑与能量分布进行分析获得弥散斑,检测图如图 12所示,记录3组数据取平均值最终结果如图 13所示,由此可见各波段均符合在0.8视场下集中80%能量时,弥散斑直径在9.2~18.4 μm区间的指标要求。

    图  12  弥散斑测试现场图
    Figure  12.  Field diagram of diffuse spot testing
    图  13  弥散斑直径图
    Figure  13.  Diffuse spot diameter map

    在20±5℃温度范围内,通过对0°视场像点观测,由公式(1)可知:

    $$ \frac{a}{f} \times \frac{{180^\circ }}{{\text{π }}} \times 3600 < 1'' $$ (1)

    式中:像元大小a为4.6 μm,焦距f为900 mm,经过计算只要像点偏移小于一个像元即可认为光轴偏角优于1″。经过观察,像点最大位移小于一个像元,故可以判断光轴偏角优于1″,满足指标要求。通过对弥散斑直径以及光轴漂移量的检测结果与仿真分析结果对比发现光机热集成分析具有可靠性,所以有必要对系统进行光机热集成分析以快速检验设计的系统是否满足指标。

    本文通过对热不敏光学系统进行结构设计,并对结构进行有限元分析,结合光机热集成分析方法,通过sigfit计算出在20±5℃下主次镜RMS值为0.13λ,将拟合得到的Zernike系数代入光学设计软件Zemax中进行仿真模拟,设计结果表明光轴偏角为0.023″优于1″,波前RMS值为0.035λ,圈入能量80%集中度弥散斑直径在9.2~18.4 μm之间,最终进行装调检测,结果显示系统轴上视场波像差RMS值为0.08λ,实现弥散斑能量80%集中度的直径在9.2~18.4 μm内,像点最大位移小于一个像元,光轴偏角优先于1″,满足项目设计指标要求。该分析方法能够准确地验证系统是否满足指标要求,极大地缩短了研制周期,能够对系统性能进行有效的评估,同时可以将该方法运用到其他光学系统光机热集成分析中。

  • 图  1   超声红外热成像技术原理

    Figure  1.   Principle of ultrasonic infrared thermal imaging technology

    图  2   冲击能量-时间曲线

    Figure  2.   Impact energy-time curves

    图  3   冲击后试件光学图像

    Figure  3.   Optical image of the specimen after impact

    图  4   超声红外热成像系统

    Figure  4.   Ultrasonic infrared thermography system

    图  5   试件3#红外序列热图

    Figure  5.   Infrared sequence thermogram of specimen 3#

    图  6   试件3#不同点温升曲线

    Figure  6.   Temperature rise curves at different points of specimen 3#

    图  7   试件3#空间温度曲线

    Figure  7.   Space temperature curves of specimen 3#

    图  8   不同冲击能量试件检测图像

    Figure  8.   Images of test specimens with different impact energies

    图  9   试件3#双阈值分割效果

    Figure  9.   Effect of double threshold segmentation on specimen 3#

    图  10   不同冲击能量试件分割效果

    Figure  10.   Effect of splitting specimens with different impact energies

    图  11   编织复合材料损伤面积与冲击能量的关系

    Figure  11.   Damage area of braided composite versus impact energy

    表  1   不同冲击能量试件损伤面积

    Table  1   Damage area of specimens with different impact energy

    Specimen Red area/mm2 Total area/mm2
    1#-10J 630 926
    2#-20J 1475 2525
    3#-30J 2909 3929
    4#-40J 2896 3383
    5#-50J 4517 6291
    下载: 导出CSV
  • [1]

    SONG S J, Waas A M, Shahwan K W, et al. Compression response of 2D braided textile composites: single cell and multiple cell micromechanics based strength predictions[J]. Journal of Composite Materials, 2008, 42(3): 2461-2482.

    [2]

    SONG S J, Waas A M, Shahwan K W, et al. Effects of matrix microcracking on the response of 2D braided textile composites subjected to compression loads[J]. Journal of Composite Materials, 2010, 44(2): 221-240. DOI: 10.1177/0021998309341345

    [3]

    Tabiei A, YI Weitao. Comparative study of predictive methods for woven fabric composite elastic properties[J]. Composites Structures, 2002, 58: 149-164. DOI: 10.1016/S0263-8223(02)00028-4

    [4]

    MORIBIDINI M. The detectability of cracks using sonic IR[J]. Journal of Applied Physics, 2009, 105(9): 093530.

    [5] 陶永强, 矫桂琼, 王波, 等. 二维编织陶瓷基复合材料应力-应变行为[J]. 力学与实践, 2009, 31(3): 39-44. https://www.cnki.com.cn/Article/CJFDTOTAL-LXYS200903012.htm

    TAO Y Q, JIAO G Q, WANG B, et al. Stress-strain behavior in 2D weave ceramic matrix composite[J]. Mechanics in Engineering, 2009, 31(3): 39-44. https://www.cnki.com.cn/Article/CJFDTOTAL-LXYS200903012.htm

    [6] 董慧民, 安学锋, 益小苏, 等. 纤维增强聚合物基复合材料低速冲击研究进展[J]. 材料工程, 2015, 43(5): 89-100. https://www.cnki.com.cn/Article/CJFDTOTAL-CLGC201505016.htm

    DONG Huimin, AN Xuefeng, YI Xiaosu, et al. Progress in research on low velocity impact properties of fibre reinforced polymer matrix composites[J]. Journal of Materials Engineering, 2015, 43(5): 89-100. https://www.cnki.com.cn/Article/CJFDTOTAL-CLGC201505016.htm

    [7]

    Tsai K H, HWAN C L, CHEN W L, et al. A parallelogram spring model for predicting the effective elastic properties of 2D braided composites[J]. Composite Structures, 2008, 83: 273-283 DOI: 10.1016/j.compstruct.2007.04.021

    [8]

    MORIBIDINI M. A calibration procedure for sonic infrared nondestructive evaluation[J]. Journal of Applied Physics, 2009, 106(2): 1-9.

    [9]

    RENSHAW J, HOLLAND S D, THOMPSON R B, et al. Vibration-induced tribological damage to fracture surfaces via vibro thermography[J]. International Journal of Fatigue, 2011, 33(7): 849-857. DOI: 10.1016/j.ijfatigue.2011.01.005

    [10] 冯辅周, 张超省, 闵庆旭, 等. 超声红外热像技术中金属平板裂纹的生热特性[J]. 红外与激光工程, 2015, 44(5): 1456-1461. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201505012.htm

    FENG F Z, ZHANG C S, MIN Q X, et al. Heating characteristics of metal plate crack in sonic IR imaging[J]. Infrared and Laser Engineering, 2015, 44(5): 1456-1461. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201505012.htm

    [11] 秦雷, 刘俊岩, 龚金龙, 等. 超声红外锁相热像技术检测金属板材表面裂纹[J]. 红外与激光工程, 2013, 42(5): 1l23-1130. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201305002.htm

    QIN L, LIU J Y, GONG J L, et al. Testing surface crack defects of sheet metal with ultrasonic lock-in thermography[J]. Infrared and Laser Engineering, 2013, 42(5): 1123-1130. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ201305002.htm

    [12] 韩梦, 尹嘉雯, 黄军科, 等. 铝合金板疲劳微裂纹超声红外成像检测的数值及实验研究[J]. 应用声学, 2022, 41(5): 727-734.

    HAN M, YIN J W, HUANG J K, et al. A comprehensive numerical and experimental study on detecting fatigue microcrack in aluminum alloy plate by vibro-thermography[J]. Applied Acoustics, 2022, 41(5): 727-734.

    [13] 习小文, 苏清风, 袁雅妮, 等. 超声红外热成像技术在航空发动机叶片裂纹的对比研究[J]. 红外技术, 2021, 43(2): 186-191. http://hwjs.nvir.cn/article/id/1c71a583-324b-4a27-9a83-cc725603f156

    XI X W, SU Q F, YUAN Y N, et al. Comparative study of using ultrasonic infrared thermography for detecting aeroengine blade cracks[J]. Infrared Technology, 2021, 43(2): 186-191. http://hwjs.nvir.cn/article/id/1c71a583-324b-4a27-9a83-cc725603f156

    [14]

    Rantala J, WU D, Busse G. Amplitude-modulated lock-in vibro-thermography for NDE of polymers and composites[J]. Research in Nondestructive Evaluation, 1996, 21(7): 215-228.

    [15] 田干, 张炜, 金国锋, 等. 超声红外热波检测多模式激励的数值仿真[J]. 江苏大学学报, 2014, 35(2): 171-175. https://www.cnki.com.cn/Article/CJFDTOTAL-JSLG201402009.htm

    TIAN G, ZHANG W, JIN G F, et al. Numerical simulation of multi-mode excitation of ultrasonic infrared thermal wave detection[J]. Journal of Jiangsu University, 2014, 35(2): 171-175. https://www.cnki.com.cn/Article/CJFDTOTAL-JSLG201402009.htm

    [16] 金国锋, 张炜, 杨正伟, 等. 界面贴合型缺陷的超声红外热波检测与识别[J]. 四川大学学报: 工程科学版, 2013, 45(2): 167-175. https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201302028.htm

    JIN G F, ZHANG W, YANG Z W, et al. Application of ultrasonic infrared thermal wave technique in detection and recognition of interface kissing damage[J]. Journal of Sichuan University: Engineering Science Edition, 2013, 45(2): 167-175. https://www.cnki.com.cn/Article/CJFDTOTAL-SCLH201302028.htm

    [17] 张炜, 罗文源, 王冬冬, 等. 基于超声热波技术的复合材料冲击损伤检测研究[J]. 固体火箭技术, 2013, 36(6): 836-841. https://www.cnki.com.cn/Article/CJFDTOTAL-GTHJ201306027.htm

    ZHANG W, LUO W Y, WANG D D, et al. Impact damage detection of composites based on ultrasonic thermography technique[J]. Journal of Solid Rocket Technology, 2013, 36(6): 836-841. https://www.cnki.com.cn/Article/CJFDTOTAL-GTHJ201306027.htm

    [18] 李胤, 杨正伟, 张炜, 等. 复合材料低速冲击损伤超声红外热波检测能力评估[J]. 仪器仪表学报, 2016, 37(5): 1124-1130. https://www.cnki.com.cn/Article/CJFDTOTAL-YQXB201605021.htm

    LI Y, YANG Z W, ZHANG W, et al. Detection capability evaluation of low velocity impact damage in composites using ultrasonic infrared thermography[J]. Chinese Journal of Scientific Instrument, 2016, 37(5): 1124-1130. https://www.cnki.com.cn/Article/CJFDTOTAL-YQXB201605021.htm

    [19] 李晓霞, 伍耐明, 段玉霞, 等. 碳纤维层合板低速冲击后的红外热波检测分析[J]. 复合材料学报, 2010, 27(6): 88-93. https://www.cnki.com.cn/Article/CJFDTOTAL-FUHE201006016.htm

    LI X X, WU N M, DUAN Y X, et al. Infrared thermal wave imaging for carbon fiber laminated boards after low velocity impact[J]. Acta Materiae Compositae Sinica, 2010, 27(6): 88-93. https://www.cnki.com.cn/Article/CJFDTOTAL-FUHE201006016.htm

  • 期刊类型引用(0)

    其他类型引用(1)

图(11)  /  表(1)
计量
  • 文章访问数:  117
  • HTML全文浏览量:  40
  • PDF下载量:  38
  • 被引次数: 1
出版历程
  • 收稿日期:  2022-04-18
  • 修回日期:  2022-04-25
  • 刊出日期:  2023-08-19

目录

/

返回文章
返回