[go: up one dir, main page]

CN115420407B - 基于fbg监测响应信息的cfrp板荷载识别方法 - Google Patents

基于fbg监测响应信息的cfrp板荷载识别方法 Download PDF

Info

Publication number
CN115420407B
CN115420407B CN202210754669.4A CN202210754669A CN115420407B CN 115420407 B CN115420407 B CN 115420407B CN 202210754669 A CN202210754669 A CN 202210754669A CN 115420407 B CN115420407 B CN 115420407B
Authority
CN
China
Prior art keywords
strain
cfrp plate
point
fbg
cfrp
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210754669.4A
Other languages
English (en)
Other versions
CN115420407A (zh
Inventor
王花平
陈聪
王省哲
张虎元
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Lanzhou University
Original Assignee
Lanzhou University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Lanzhou University filed Critical Lanzhou University
Priority to CN202210754669.4A priority Critical patent/CN115420407B/zh
Publication of CN115420407A publication Critical patent/CN115420407A/zh
Application granted granted Critical
Publication of CN115420407B publication Critical patent/CN115420407B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01LMEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
    • G01L1/00Measuring force or stress, in general
    • G01L1/24Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet
    • G01L1/242Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet the material being an optical fibre
    • G01L1/246Measuring force or stress, in general by measuring variations of optical properties of material when it is stressed, e.g. by photoelastic stress analysis using infrared, visible light, ultraviolet the material being an optical fibre using integrated gratings, e.g. Bragg gratings
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/16Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
    • G01B11/165Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by means of a grating deformed by the object

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开一种基于FBG监测响应信息的CFRP板荷载识别方法,首先在CFRP板上粘贴或内嵌FBG传感器;基于FBG传感器监测到的波长时程信号,通过波长与应变的转化公式计算CFRP板上各FBG测点的应变值;然后根据CFRP板上离散点的应变值,通过基于格林函数的CFRP板应变场重构算法得到CFRP板的应变场;再利用基于峰突强度的波峰识别算法得到CFRP板应变场的突变信息,对CFRP板应变场和突变信息进行二维笛卡尔平面表达并评估,实现了CFRP板的荷载识别。本发明的特点是具有实时性高、高精度,准确和直观在线实时测量CFRP板的应变特性,进行CFRP板的荷载识别和在线检测。

Description

基于FBG监测响应信息的CFRP板荷载识别方法
技术领域
本发明属于结构智能健康监测和检测领域,涉及到的是一种基于FBG监测响应信息的CFRP板荷载识别方法。
背景技术
碳纤维增强复合材料(CFRP)板由于具有比强度高、比刚度高、减振性能好、耐腐蚀、质量轻等诸多优点,被广泛应用于航天航空、汽车、轨道交通、能源储运、海洋管道等领域。然而,在工程应用中复合材料易遭受各种荷载作用,如飞鸟的撞击、冰雹的砸落等,这些不经意间发生的荷载作用极有可能在荷载作用点处造成结构的损伤;此外,CFRP板在遭受冲击、疲劳交变载荷、雷击等外界因素影响后发生损伤,不进行及时有效的检测和修复,损伤会进一步扩展,造成严重危害。
结构健康监测兴起于航天领域,20世纪80年代开始广泛应用于土木工程、机械等领域。结构健康监测最初是进行结构的荷载监测,逐渐向结构损伤诊断(损伤部位识别和损伤程度估计等)、结构剩余寿命的预估等方面发展,所以,结构荷载监测是结构健康监测的一个技术基础。另外,一些结构损伤识别问题可转化为结构荷载识别问题,因此,结构荷载识别成为结构损伤识别的一个间接手段。研究结构荷载监测对于促进结构健康监测的发展具有重要学术和现实意义。
现有的CFRP板荷载识别方法包括超声波检测、射线检测和声发射检测。超声检测和射线检测都是静态检测方法,不能实时动态检测结构运动过程中的荷载。在声发射测试方法中,通过压电传感器检测荷载,以捕获材料发生变形时释放的弹性波。然而,弹性波在复合材料中传播过程中,由于频散、模式转换、衰减和界面反射,容易引入干扰和噪声,而压电传感器和复合材料在连接表面容易发生脱胶失效,因此这种方法很难适应CFRP板的复杂工作条件。光纤布拉格光栅(FBG)传感器具有质量轻、体积小、易于嵌入复合材料内部、易于实现准分布式测量等优点,可以在线检测复合材料板。将FBG传感器技术引入复合材料板荷载识别后,可以有效克服超声检测和射线检测等传统荷载识别技术的缺点。
因此,可以利用FBG传感器的这些优点将系列FBG传感器粘贴在CFRP 板上或者嵌入CFRP板中构建FBG传感器网络,在线识别其荷载信息。
发明内容
本发明所要解决的技术问题是要提供一种基于FBG监测响应信息的CFRP 板荷载识别方法。
为了解决以上的技术问题,本发明提供了一种CFRP板荷载的实时识别方法,首先通过基于格林函数的应变场重构方法得到CFRP板的应变响应特征场,然后辨识荷载作用下CFRP板应变响应特征场的突变信息,进而得到荷载信息,包括以下的步骤:
步骤S1:根据CFRP板的尺寸大小,在CFRP板上粘贴若干个FBG传感器形成FBG传感器网络,或者制造内嵌FBG传感网络的CFRP板。将FBG传感网络依次与解调仪和计算机连接,组成CFRP板的荷载识别监测系统。
步骤S2:基于FBG传感器监测到的波长时程信号,通过波长与应变的转化公式计算CFRP板上各FBG测点的应变值εi
步骤S3:重构CFRP板上点P0的应变值。用于重构的最接近P0的已知点数为k,k应根据FBG传感器总数量和所需重构速度来取值。点P0到其他已知点的距离应首先通过以下公式计算:
r0i=|xi,x0|,i=1,2,...,n (1)
其中,r0i是点P0与第i个已知应变点之间的距离,x0是P0的位置向量,xi是第i个已知应变点的位置向量。所有距离都使用合并排序算法按升序排序。排序后,离P0最近的前k个应变已知点的坐标和应变(xj,yj,zj)j=1,2,...,k就获得了。
步骤S4:计算格林函数矩阵。设置坐标矩阵 X=[x1 x2 ... xk]T,Y=[y1 y2 ...yk]T,应变矩阵Z=[z1 z2 ... zk]T。k×k格林函数矩阵G为:
在上式中,元素dij表示换算变量,其计算公式为:
其中,表示xy平面内点i到点j的直线距离。使用双调和格林函数/>作为基函数,以获得最小曲率曲面。根据格林函数矩阵 G和应变矩阵Z,权重矩阵W可以计算为
W=G-1Z (4)
步骤S5:计算重构应变点P0的1×k格林函数矩阵GP,公式如下:
GP=[d01 d02 ... d0k] (5)
其中d01,d02,……,d0k的计算公式与步骤S4的dij的计算公式相同。然后,应变重构点P0的应变值为
zP=GPW (6)
步骤S6:重复步骤S3-S6,就可以得到CFRP板上所有点Pi的应变值zPi,即对CFRP板的应变场进行了重构。
步骤S7:荷载识别:
重构完成的CFRP板应变场,若荷载产生,CFRP板的应变场会出现突变,采用基于峰突强度的波峰识别算法得到CFRP板应变场的突变位置,该算法先计算出CFRP板应变场的所有极值点,然后计算各个极值点与该极值点周围极值点的应变差值(即峰突强度)是否大于预先设置的阈值,若大于则该极值点为突变点。突变位置既是荷载作用位置,CFRP板应变场突变位置的应变大小可以反映荷载作用的大小。
本发明基于CFRP板荷载作用时的应变场突变特性,采用了基于格林函数的应变场重构方法将离散FBG测点的监测信号进行处理转化为CFRP板的连续应变场,当荷载产生时,该重构的应变场会呈现突变信息,清晰、直观地识别荷载的发生。
本发明的优越功效在于:
1)本发明基于FBG监测数据的实时处理和应变场重构,可实时显示出CFRP 板的应变场,并且CFRP板有荷载作用和无荷载作用下的显示结果有明显差异,能直观、准确地感知荷载的发生,服务于CFRP板的荷载识别和运营维护工作,并且满足在线监测环境下的实时性、连续性的要求,易于实现,具有很高的工程应用价值和广阔的应用前景;
2)本发明方法过程简单,通过实时监测数据处理分析结果表明,极易识别和感知CFRP板荷载的产生,通过计算应变场的瞬时指标可实现荷载实时识别和在线测量,并且计算效率高,可持续稳定运行;
3)本发明方法具有实时性高(秒级)﹑高精度﹑准确和直观等特点;
4)本发明方法的应用场景广泛。
附图说明
附图1是本发明方法的流程图。
附图2是没有荷载作用时CFRP板的重构应变场图。
附图3是荷载作用时CFRP板的重构应变场图。
附图4是荷载作用时CFRP板的应变场出现突变信息。
图中:1荷载作用点;2FBG传感器;3CFRP板
具体实施方式
以下结合技术方案(和附图)详细叙述本发明的具体实施方式。
附图1示出了本发明实施的方法流程图。本发明提供了一种基于FBG监测响应信息的CFRP板荷载识别方法,包括以下步骤:
S1:根据CFRP板3的尺寸大小,在CFRP板上粘贴若干FBG传感器2形成FBG传感网络,如附图2,或者制造内嵌FBG传感网络的CFRP板3。将FBG 传感网络与解调仪和计算机连接,组成CFRP板3的荷载识别监测系统。
S2:基于FBG传感器2监测到的波长时程信号,通过波长与应变的转化公式计算CFRP板3上各FBG测点的应变值εi
S3:重构CFRP板3上点P0的应变值。用于重构的最接近P0的已知点数为k, k应根据FBG传感器2总数量和所需重构速度来取值。点P0到其他已知点的距离应首先通过以下公式计算:
r0i=|xi,x0|,i=1,2,...,n
其中,r0i是点P0与第i个已知应变点之间的距离,x0是P0的位置向量,xi是第i个已知应变点的位置向量。所有距离都使用合并排序算法按升序排序。排序后,离P0最近的前k个应变已知点的坐标和应变(xj,yj,zj)j=1,2,...,k就获得了。
S4:计算格林函数矩阵。设置坐标矩阵X=[x1 x2 ... xk]T和 Y=[y1 y2 ... yk]T,应变矩阵Z=[z1 z2 ... zk]T。k×k格林函数矩阵G为:
在上式中,元素dij表示换算变量,其计算公式为:
其中,表示xy平面内点i到点j的直线距离。使用双调和格林函数/>作为基函数,以获得最小曲率曲面。根据格林函数矩阵 G和应变矩阵Z,权重矩阵W可以表示为:W=G-1Z
步骤S5:计算重构应变点P0的1×k格林函数矩阵GP,公式如下:
GP=[d01 d02 ... d0k]
其中,d01,d02,……,d0k的计算公式与步骤S4的dij的计算公式相同。然后,应变重构点P0的应变值为zP=GPW。
S6:重复步骤S3-S6,就可以得到CFRP板3上所有点Pi的应变值zPi,即对 CFRP板的应变场进行了重构。
S7:荷载识别
通过计算重构出CFRP板3的应变场,并且在CFRP板3平面内建立二维笛卡尔坐标系,x轴和y轴表示CFRP板3上的横纵坐标,用颜色表示CFRP板3 的应变场大小;若没有荷载作用,CFRP板3的应变场是均匀的,如图3所示;在荷载作用点1位置,CFRP板3的应变场会出现突变信息,如图4所示。
本发明可用于航天航空、轨道交通等复合材料结构(如飞机机翼、高速铁路机车结构等)的CFRP板3类结构荷载识别,为CFRP板3的监控维护提供服务;也可用于CFRP板3实验室内的性能测试和疲劳检验过程监控;亦可用于具有荷载识别需求的其它板类结构,如建筑物中的楼盖、桥梁面板等。

Claims (1)

1.一种基于FBG监测响应信息的CFRP板荷载识别方法,其特征在于,包括:
步骤S1:根据CFRP板(3)的尺寸大小,在CFRP板(3)上粘贴若干个FBG传感器(2)形成FBG传感网络,或者制造内嵌FBG传感网络的CFRP板(3),将FBG传感网络依次与解调仪和计算机连接,组成CFRP板(3)的荷载识别监测系统;
步骤S2:基于FBG传感器(2)监测到的波长时程信号,通过波长与应变的转化公式计算CFRP板(3)上各FBG测点的应变值εi
步骤S3:重构CFRP板(3)上点P0的应变值,用于重构的最接近P0的已知点数为k,k应根据FBG传感器(2)总数量和所需重构速度来取值,点P0到其他已知点的距离应首先通过以下公式计算:
r0i=|xi,x0|,i=1,2,...,n
其中,r0i是点P0与第i个已知应变点之间的距离,x0是P0的位置向量,xi是第i个已知应变点的位置向量,所有距离都使用合并排序算法按升序排序,排序后,获得离P0最近的前k个应变已知点的坐标和应变(xj,yj,zj)j=1,2,...,k;
步骤S4:计算格林函数矩阵:设置坐标矩阵X=[x1 x2 ... xk]T和Y=[y1 y2 ... yk]T,应变矩阵Z=[z1 z2 … zk]T,k×k格林函数矩阵G为:
在上式中,元素dij表示换算变量,其计算公式为:
其中,表示xy平面内点i到点j的直线距离。使用双调和格林函数/>作为基函数,以获得最小曲率曲面。根据格林函数矩阵G和应变矩阵Z,权重矩阵W可以计算为:W=G-1Z
步骤S5:计算重构应变点P0的1×k格林函数矩阵GP,公式如下:
GP=[d01 d02 … d0k]
其中,d01,d02,……,d0k的计算公式与步骤S4的dij的计算公式相同,然后,应变重构点P0的应变值为:zP=GPW;
步骤S6:重复步骤S3-S6,就可以得到CFRP板(3)上所有点Pi的应变值zPi,即对CFRP板(3)的应变场进行了重构;
步骤S7:荷载识别:
重构完成的CFRP板(3)应变场,若荷载产生,CFRP板(3)的应变场会出现突变,采用基于峰突强度的波峰识别算法得到CFRP板(3)应变场的突变位置,该算法先计算出CFRP板(3)应变场的所有极值点,然后计算各个极值点与该极值点周围极值点的应变差值是否大于预先设置的阈值,若大于则该极值点为突变点,突变位置既是荷载作用点(1)位置,CFRP板(3)应变场突变位置的应变大小可以反映荷载作用的大小。
CN202210754669.4A 2022-06-30 2022-06-30 基于fbg监测响应信息的cfrp板荷载识别方法 Active CN115420407B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210754669.4A CN115420407B (zh) 2022-06-30 2022-06-30 基于fbg监测响应信息的cfrp板荷载识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210754669.4A CN115420407B (zh) 2022-06-30 2022-06-30 基于fbg监测响应信息的cfrp板荷载识别方法

Publications (2)

Publication Number Publication Date
CN115420407A CN115420407A (zh) 2022-12-02
CN115420407B true CN115420407B (zh) 2024-04-09

Family

ID=84195623

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210754669.4A Active CN115420407B (zh) 2022-06-30 2022-06-30 基于fbg监测响应信息的cfrp板荷载识别方法

Country Status (1)

Country Link
CN (1) CN115420407B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115629130B (zh) * 2022-12-21 2023-04-28 国网天津市电力公司电力科学研究院 盆式绝缘子残余应力的成像方法、系统及试块制备方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2334051A1 (en) * 2000-02-02 2001-08-02 Wolfgang Gleine Method and sensor arrangement for measuring temperature and strain using an optical fiber embedded in a cover layer on a substrate
CN105588673A (zh) * 2015-12-31 2016-05-18 武汉理工大学 一种光纤光栅传感器监测模具与构件作用力的方法
WO2017202139A1 (zh) * 2016-05-26 2017-11-30 东南大学 一种基于长标距应变影响包络线的桥梁损伤识别方法
CN110068406A (zh) * 2019-04-01 2019-07-30 南京航空航天大学 基于静载识别的四边简支薄板结构光纤应变场重构方法
CN111174934A (zh) * 2020-01-08 2020-05-19 兰州大学 复合结构界面损伤劣化过程的光纤监测方法
CN112613129A (zh) * 2020-12-30 2021-04-06 大连理工大学 一种基于逆有限元与微元动态响应法的复合材料板壳结构健康监测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2334051A1 (en) * 2000-02-02 2001-08-02 Wolfgang Gleine Method and sensor arrangement for measuring temperature and strain using an optical fiber embedded in a cover layer on a substrate
CN105588673A (zh) * 2015-12-31 2016-05-18 武汉理工大学 一种光纤光栅传感器监测模具与构件作用力的方法
WO2017202139A1 (zh) * 2016-05-26 2017-11-30 东南大学 一种基于长标距应变影响包络线的桥梁损伤识别方法
CN110068406A (zh) * 2019-04-01 2019-07-30 南京航空航天大学 基于静载识别的四边简支薄板结构光纤应变场重构方法
CN111174934A (zh) * 2020-01-08 2020-05-19 兰州大学 复合结构界面损伤劣化过程的光纤监测方法
CN112613129A (zh) * 2020-12-30 2021-04-06 大连理工大学 一种基于逆有限元与微元动态响应法的复合材料板壳结构健康监测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
复合层合板低速冲击应变测试及其状态监测;王利恒;;压电与声光;20120815(04);全文 *

Also Published As

Publication number Publication date
CN115420407A (zh) 2022-12-02

Similar Documents

Publication Publication Date Title
Du et al. Damage detection techniques for wind turbine blades: A review
Yang et al. Ultrasonic detection methods for mechanical characterization and damage diagnosis of advanced composite materials: A review
Ye et al. A multi-scale model for studying failure mechanisms of composite wind turbine blades
Qing et al. Machine learning based quantitative damage monitoring of composite structure
Alokita et al. Recent advances and trends in structural health monitoring
CN105956216B (zh) 大跨钢桥基于均匀温度响应监测值的有限元模型修正方法
CN105930571B (zh) 基于单位温度响应监测值的大跨钢桥有限元模型修正方法
CN103076394B (zh) 基于振动识别频率和振型综合的海洋平台安全评定的方法
CN106770648B (zh) 一种复合材料刚度系数的无损测量方法
CN107301282B (zh) 基于多源监测时序数据的混凝土坝力学参数反演方法
CN113720907B (zh) 一种轮廓与深度顺序辨识的复合材料分层损伤辨识方法
Hiche et al. A strain amplitude-based algorithm for impact localization on composite laminates
CN101706355A (zh) 基于NExT/ARMA的结构响应分析方法
Altabey et al. Tensile creep monitoring of basalt fiber-reinforced polymer plates via electrical potential change and artificial neural network
CN115420407B (zh) 基于fbg监测响应信息的cfrp板荷载识别方法
CN110108430B (zh) 一种纤维增强复合材料梁的高精度损伤定位方法
CN104155092B (zh) 一种风力机叶片静态分析方法
Zhou et al. Review on Structural Health Monitoring in Metal Aviation Based on Fiber Bragg Grating Sensing Technology
Kirikera et al. Monitoring multi-site damage growth during quasi-static testing of a wind turbine blade using a structural neural system
Malekimoghadam et al. A critical review on the structural health monitoring methods of the composite wind turbine blades
Kim et al. Impact localization for composite plate using the modified error-outlier algorithm with Pugh’s concept selection under various temperatures
Wang et al. Predicting delamination in composite laminates through semi-analytical dynamic analysis and vibration-based quantitative assessment
CN113688544B (zh) 一种主被动结合的复合材料损伤定量辨识方法
Derkevorkian et al. Computational studies of a strain-based deformation shape prediction algorithm for control and monitoring applications
Pillai et al. Monitoring Structural Health Using Digital Twin

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant