CN113029405B - 基于光纤应变组合桥路的桨叶弯矩解耦与标定方法 - Google Patents
基于光纤应变组合桥路的桨叶弯矩解耦与标定方法 Download PDFInfo
- Publication number
- CN113029405B CN113029405B CN202110158338.XA CN202110158338A CN113029405B CN 113029405 B CN113029405 B CN 113029405B CN 202110158338 A CN202110158338 A CN 202110158338A CN 113029405 B CN113029405 B CN 113029405B
- Authority
- CN
- China
- Prior art keywords
- wave
- bending moment
- blade
- shimmy
- section
- 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
Links
- 238000005452 bending Methods 0.000 title claims abstract description 116
- 238000000034 method Methods 0.000 title claims abstract description 40
- 239000013307 optical fiber Substances 0.000 title claims abstract description 28
- 230000008878 coupling Effects 0.000 claims abstract description 20
- 238000010168 coupling process Methods 0.000 claims abstract description 20
- 238000005859 coupling reaction Methods 0.000 claims abstract description 20
- 239000000835 fiber Substances 0.000 claims description 43
- 102100022116 F-box only protein 2 Human genes 0.000 claims description 27
- 101000824158 Homo sapiens F-box only protein 2 Proteins 0.000 claims description 27
- 102100024513 F-box only protein 6 Human genes 0.000 claims description 26
- 101001052796 Homo sapiens F-box only protein 6 Proteins 0.000 claims description 26
- 238000005259 measurement Methods 0.000 claims description 12
- 230000009471 action Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 6
- 238000010998 test method Methods 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims 2
- 230000005484 gravity Effects 0.000 claims 2
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 230000036541 health Effects 0.000 abstract description 2
- 230000005540 biological transmission Effects 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 4
- 230000008054 signal transmission Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/24—Measuring 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/242—Measuring 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/246—Measuring 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L25/00—Testing or calibrating of apparatus for measuring force, torque, work, mechanical power, or mechanical efficiency
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于光纤应变组合桥路的桨叶弯矩解耦与标定方法,属于结构健康监测技术领域。包括以下步骤,步骤一:传感器网络布置;步骤二:测量桨叶剖面预扭角并实现挥舞弯矩解耦;步骤三:计算挥舞弯矩标定参量△λ 挥对桨叶挥舞弯矩进行标定;步骤四:计算耦合系数C,选择桨叶结构摆振弯矩传感器最佳解耦组合桥路;步骤五:计算摆振弯矩标定参量△λ 摆对桨叶摆振弯矩进行标定。本发明所述方法通过获取光纤Bragg光栅传感器中心波长偏移量特征参量组合△λ 挥和△λ 摆,能够实现挥舞与摆振弯矩解耦并成功运用于挥舞与摆振弯矩标定。
Description
技术领域
本发明属于结构健康监测的弯矩监测技术领域,具体提出了一种基于光纤Bragg光栅传感器的直升机桨叶结构挥舞弯矩与摆振弯矩解耦与标定方法。
背景技术
随着现代工业设备逐渐朝着自动化方向发展,为满足机械设备性能指标要求,叶片类旋转结构不断向高转速、高承载、高强度、柔性化及轻量化的方向发展,现代新型直升机对于桨叶结构设计要求也越来越高,其安全性及可靠性问题引起高度重视。直升机桨叶作为直升机起飞悬停的主要承载与提供升力的结构,其服役状态的稳定性直接关系着直升机飞行安全,需研究桨叶结构载荷测量技术对直升机桨叶结构在服役状态下的受载情况进行实时监测,为直升机飞行安全性评估及服役寿命预测提供依据。
以往试验研究中,直升机桨叶类旋转叶片一般是采用传统应变片来获取结构应变参数,基本流程是将应变片粘贴于桨叶结构表面,采集直升机模拟飞行状态下单个节点的应变响应信号。起初信号传输是基于集流环系统,可实现旋转桨叶电信号传输,之后将无线电Wi-Fi无线传输技术与应变测量结合,形成了一套适用于旋转桨叶结构的无线电遥测动态应变测试分析系统。采用传统电阻应变片对直升机桨叶动态应变测量存在一定缺陷,在直升机服役过程中,电阻应变片易受电磁干扰造成测量数据错误或缺失。为获得桨叶结构整体应变分布状态,需要配置较多应变片及传输线缆组成应变测量桥路,桥路接线工作量大且布线成本较高。同时,桨叶结构布设过多传输电缆可能会改变直升机桨叶自身结构刚度,影响桨叶表面气动特性,无法满足直升机桨叶服役状态下应变与弯矩载荷长期监测需求。
因此,针对直升机桨叶工作环境恶劣,外界干扰较多的特点,选择使用非光电的接触式测量方法,并且选择相对先进且技术日趋成熟的光纤Bragg光栅传感器为测试元件。与传统传感器相比,光纤Bragg光栅传感器具有多参量测量、传输纤线纤芯细、质量轻、体积小、抗拉刚度强、耐高温、耐腐蚀、集信号传感与传输一体化、抗电磁干扰且易于构成空分/波分复用型监测阵列等独特优点。由于桨叶剖面的不对称性,在实际弯矩状态监测过程中挥舞弯矩和摆振弯矩会存在耦合问题,针对目前桨叶弯矩状态监测方法的不足,需要研究无需大量先验知识、简单快速、方便可靠、实用性强的新方法。为此,本发明提出了一种基于光纤Bragg光栅传感器的桨叶结构挥舞/摆振弯矩解耦与标定方法。
发明内容
发明目的:针对传统应变片对直升机桨叶结构弯矩载荷测量易受外界干扰问题,本发明提供一种基于光纤Bragg光栅传感器的直升机桨叶结构挥舞/摆振弯矩解耦与标定方法。该方法采用光纤Bragg光栅传感器感知桨叶典型剖面结构表面弦向不同位置在挥舞与摆振弯矩载荷下的中心波长响应信号,通过不同传感器组合运算,实现挥舞与摆振弯矩解耦并提供一种挥舞与摆振弯矩标定方程计算方法。根据此方法可以分别可以实现桨叶挥舞和摆振弯矩解耦,在旋转飞行状态下分别对挥舞和摆振弯矩载荷进行测量。
技术方案:为实现上述目的,本发明的所采用的技术方案包括以下步骤:
一种基于光纤应变组合桥路的桨叶弯矩解耦与标定方法,其特征在于,包括以下步骤:
步骤一:桨叶结构光纤Bragg光栅传感器网络布置;具体为:
桨叶传感网络布局分为挥舞弯矩解耦布局与摆振弯矩解耦布局,其中用于挥舞弯矩解耦标定的光纤Bragg光栅传感器FBG1和FBG2分别粘贴于桨叶结构弯矩待测剖面的上、下表面距翼型后缘1/4弦线对称位置;用于摆振弯矩解耦标定的光纤Bragg光栅传感器除FBG1和FBG2外还包括分别粘贴于桨叶结构弯矩待测剖面上、下表面距翼型前缘1~10mm弦线对称位置的FBG3和FBG5,以及分别粘贴于桨叶结构弯矩待测剖面上、下表面距翼型前缘11~20mm弦线对称位置的FBG4和FBG6。上述光纤Bragg光栅传感器粘贴方向均平行于桨叶展向以获取桨叶结构展向应变,以此构成桨叶传感器网络,光纤Bragg光栅传感器粘贴位置所对应剖面即为桨叶结构弯矩待测剖面;
步骤二:采集摆振方向载荷作用下光纤光栅Bragg传感器FBG1与FBG2中心波长偏移量,测量桨叶剖面预扭角以实现挥舞弯矩解耦;具体为:
将桨叶根部固定于标定旋转台夹具,保持桨叶翼型前缘向下,调整桨叶根部翼型弦线与竖直面夹角为θ,采用翼型夹具下方悬挂标准砝码方式对桨叶结构施加摆振方向载荷m摆,依次改变夹角θ,分别记录同加载不同夹角θ条件下光纤Bragg光栅传感器FBG1与FBG2中心波长偏移量△λ1与△λ2;
上、下表面FBG1与FBG2传感器中心波长偏移量之差与夹角θ为线性关系,表达式为:
θ=kΔλ12+b (1)
式中△λ12=△λ1-△λ2,为FBG1与FBG2中心波长偏移量之差,k为拟合系数,b为拟合曲线截距。通过上述试验测量结果,当两个光纤Bragg光栅传感器中心波长偏移量之差△λ12为0时,得到对应角度θ即拟合曲线截距b为待测桨叶预扭角θ0;
步骤三:挥舞方向分级加载作用下采集光纤光栅Bragg传感器FBG1与FBG2中心波长偏移量,计算挥舞弯矩标定参量△λ挥对桨叶挥舞弯矩进行标定;具体为:
将标定旋转台旋转90°使得桨叶根部剖面与水平面夹角大小为预扭角θ0,此时待标定剖面翼型弦线保持水平,对桨叶结构施加挥舞方向载荷m挥,桨叶挥舞标定加载产生的挥舞弯矩为:
M挥=m挥gL (2)
式中L为挥舞载荷桨叶结构加载剖面与待标定剖面之间距离,FBG1与FBG2传感器中心波长偏移量之差△λ挥=△λ1-△λ2,该偏移量之差数值仅由剖面挥舞弯矩M挥产生。桨叶结构剖面挥舞弯矩与桨叶表面FBG传感器中心波长偏移量△λ关系式如下:
其中λ为光纤光栅Bragg传感器中心波长初始值,k为常量,h为距离,I为桨叶剖面惯性距,Pe为光纤光栅有效弹光常数。令K=I/λ(1-Pe)kh,桨叶结构剖面挥舞弯矩与FBG1传感器中心波长偏移量△λ1及FBG2传感器中心波长偏移量△λ2关系分别为:
M挥=K1Δλ1 M挥=K2Δλ2 (4)
用于标定挥舞弯矩的FBG传感器波长参量△λ挥=△λ1-△λ2,其中:
与上述试验方法相同,分别在施加挥舞方向载荷2m挥、3m挥、4m挥…N·m挥加载状态下记录波长参量△λ挥,计算N次分级加载下标定系数K挥的平均值即得到桨叶叶结构挥舞弯矩与FBG1与FBG2传感器中心波长偏移量之差的关系方程:
步骤四:分别采集挥舞方向载荷与摆振方向载荷作用下光纤光栅Bragg传感器中心波长偏移量,计算耦合系数C,选择桨叶结构摆振弯矩传感器最佳解耦组合桥路;具体为:
根据步骤二所测得桨叶剖面预扭角,将桨叶根部固定于标定旋转台夹具,保持桨叶翼型前缘向下,桨叶根部翼型弦线与竖直面夹角大小为剖面预扭角θ0,对桨叶结构施加摆振方向载荷m摆,分别记录光纤Bragg光栅传感器FBG1~FBG6在未加载与加载状态下的中心波长信号,计算FBG1~FBG6传感器中心波长偏移量△λ1~△λ6;利用标定旋转台将桨叶旋转90°,使得待标定剖面翼型弦线保持水平,对桨叶结构施加挥舞方向载荷m挥,分别记录光纤Bragg光栅传感器FBG1~FBG6在未加载与加载状态下的中心波长信号,计算FBG1~FBG6传感器中心波长偏移量△λ’1~△λ’6;
评估摆振解耦效果的耦合系数C计算公式为:
式中i取值为3或5,j取值为4或6,四组传感器组合方式计算得到耦合系数分别为:
比较耦合系数C1、C2、C3及C4数值大小,选择耦合系数最小的组合方式并将该传感器组合桥路偏移量作为用于标定摆振弯矩的FBG传感器波长参量,其计算公式为;
Δλ摆=Δλ1+Δλ2-(Δλi+Δλj) (10)
步骤五:采集摆振方向分级载荷作用下光纤光栅Bragg传感器FBG1~FBG6中心波长偏移量,计算摆振弯矩标定参量△λ摆对桨叶摆振弯矩进行标定;具体为:
将桨叶根部固定于标定旋转台夹具,保持桨叶翼型前缘向下,桨叶根部翼型弦线与竖直面夹角大小为剖面预扭角θ0,对桨叶结构施加摆振方向载荷m摆,桨叶摆振标定加载产生的摆振弯矩为:
M摆=m摆gL (11)
式中L为摆振载荷桨叶结构加载剖面与待标定剖面之间距离,桨叶结构剖面摆振弯矩与桨叶表面FBG1~FBG6传感器中心波长偏移量之差△λ关系式如下:
其中λ为光纤光栅Bragg传感器中心波长初始值,k为常量,h为距离,I为桨叶剖面惯性距,Pe为光纤光栅有效弹光常数。令K=I/λ(1-Pe)kh,桨叶结构剖面摆振弯矩与FBGn传感器中心波长偏移量△λn关系为:
用于标定摆振弯矩的FBG传感器波长参量△λ摆=△λ1+△λ2-(△λi+△λj),为步骤四所选择耦合系数最小的FBG传感器组合桥路波长偏移量,将式(10)代入得:
与上述试验方法相同,分别在施加摆振方向载荷2m摆、3m摆、4m摆…N·m挥加载状态下记录波长参量△λ摆,计算N次分级加载下标定系数K摆的平均值即得到桨叶叶结构摆振弯矩与FBG传感器组合桥路中心波长偏移量之差的关系方程:
所述的一种基于光纤Bragg光栅传感器的直升机桨叶结构挥舞与摆振弯矩解耦与标定方法,其特征在于:摆振方向载荷m摆与环境温度变化△T不影响挥舞弯矩标定参量△λ挥。
所述的一种基于光纤Bragg光栅传感器的直升机桨叶结构挥舞与摆振弯矩解耦与标定方法,其特征在于:挥舞方向载荷m挥与环境温度变化△T不影响摆振弯矩标定参量△λ摆。
有益效果:本发明适用于直升机桨叶弯矩载荷测量工程应用领域,所具备的优点是:仅需采用6个光纤Bragg光栅传感器构成传感网络,相较于传统传感方式,具有线路布置简单、抗电磁干扰能力强等优点。其次,本发明通过对挥舞与摆振载荷作用下FBG传感器中心波长信号进行数值运算,选取了挥舞弯矩标定参量△λ挥和摆振弯矩标定参量△λ摆,实现挥舞/摆振解耦的同时消除了温度对弯矩载荷测量的影响,简化了挥舞/摆振弯矩解耦过程。最后,采用光纤应变传感网络获取结构表面应变,相比于传统应变片测量方法无需布置大量信号传输线缆,不仅利于旋转工况下信号传输,还避免了对被测桨叶结构参数性能造成不利影响。
附图说明
图1是桨叶结构光纤Bragg光栅传感器布置图;
图2是桨叶待测剖面翼型示意图;
图3是挥舞方向加载示意图;
图4是摆振方向加载示意图;
图中标号名称:1.夹具2.光纤测量试验桨叶3.挥舞标定FBG1 4.挥舞标定FBG2 5.翼型夹具6.砝码7.摆振标定FBG。
具体实施方式
步骤一:桨叶结构光纤Bragg光栅传感器网络布置;具体为:
桨叶传感网络布局分为挥舞弯矩解耦布局与摆振弯矩解耦布局,其中用于挥舞弯矩解耦标定的光纤Bragg光栅传感器FBG1和FBG2分别粘贴于桨叶结构弯矩待测剖面的上、下表面距翼型后缘1/4弦线对称位置;用于摆振弯矩解耦标定的光纤Bragg光栅传感器除FBG1和FBG2外还包括分别粘贴于桨叶结构弯矩待测剖面上、下表面距翼型前缘1~10mm弦线对称位置的FBG3和FBG5,以及分别粘贴于桨叶结构弯矩待测剖面上、下表面距翼型前缘11~20mm弦线对称位置的FBG4和FBG6。上述光纤Bragg光栅传感器粘贴方向均平行于桨叶展向以获取桨叶结构展向应变,以此构成桨叶传感器网络,光纤Bragg光栅传感器粘贴位置所对应剖面即为桨叶结构弯矩待测剖面;
步骤二:采集摆振方向载荷作用下光纤光栅Bragg传感器FBG1与FBG2中心波长偏移量,测量桨叶剖面预扭角以实现挥舞弯矩解耦;具体为:
将桨叶根部固定于标定旋转台夹具,保持桨叶翼型前缘向下,调整桨叶根部翼型弦线与竖直面夹角为θ,采用翼型夹具下方悬挂标准砝码方式对桨叶结构施加摆振方向载荷m摆,依次改变夹角θ,分别记录同加载不同夹角θ条件下光纤Bragg光栅传感器FBG1与FBG2中心波长偏移量△λ1与△λ2;
上、下表面FBG1与FBG2传感器中心波长偏移量之差与夹角θ为线性关系,表达式为:
θ=kΔλ12+b (1)
式中△λ12=△λ1-△λ2,为FBG1与FBG2中心波长偏移量之差,k为拟合系数,b为拟合曲线截距。通过上述试验测量结果,当两个光纤Bragg光栅传感器中心波长偏移量之差△λ12为0时,得到对应角度θ即拟合曲线截距b为待测桨叶预扭角θ0;
步骤三:挥舞方向分级加载作用下采集光纤光栅Bragg传感器FBG1与FBG2中心波长偏移量,计算挥舞弯矩标定参量△λ挥对桨叶挥舞弯矩进行标定;具体为:
将标定旋转台旋转90°使得桨叶根部剖面与水平面夹角大小为预扭角θ0,此时待标定剖面翼型弦线保持水平,对桨叶结构施加挥舞方向载荷m挥,桨叶挥舞标定加载产生的挥舞弯矩为:
M挥=m挥gL (2)
式中L为挥舞载荷桨叶结构加载剖面与待标定剖面之间距离,FBG1与FBG2传感器中心波长偏移量之差△λ挥=△λ1-△λ2,该偏移量之差数值仅由剖面挥舞弯矩M挥产生。桨叶结构剖面挥舞弯矩与桨叶表面FBG传感器中心波长偏移量△λ关系式如下:
其中λ为光纤光栅Bragg传感器中心波长初始值,k为常量,h为距离,I为桨叶剖面惯性距,Pe为光纤光栅有效弹光常数。令K=I/λ(1-Pe)kh,桨叶结构剖面挥舞弯矩与FBG1传感器中心波长偏移量△λ1及FBG2传感器中心波长偏移量△λ2关系分别为:
M挥=K1Δλ1 M挥=K2Δλ2 (4)
用于标定挥舞弯矩的FBG传感器波长参量△λ挥=△λ1-△λ2,其中:
与上述试验方法相同,分别在施加挥舞方向载荷2m挥、3m挥、4m挥…N·m挥加载状态下记录波长参量△λ挥,计算N次分级加载下标定系数K挥的平均值即得到桨叶叶结构挥舞弯矩与FBG1与FBG2传感器中心波长偏移量之差的关系方程:
步骤四:分别采集挥舞方向载荷与摆振方向载荷作用下光纤光栅Bragg传感器中心波长偏移量,计算耦合系数C,选择桨叶结构摆振弯矩传感器最佳解耦组合桥路;具体为:
根据步骤二所测得桨叶剖面预扭角,将桨叶根部固定于标定旋转台夹具,保持桨叶翼型前缘向下,桨叶根部翼型弦线与竖直面夹角大小为剖面预扭角θ0,对桨叶结构施加摆振方向载荷m摆,分别记录光纤Bragg光栅传感器FBG1~FBG6在未加载与加载状态下的中心波长信号,计算FBG1~FBG6传感器中心波长偏移量△λ1~△λ6;利用标定旋转台将桨叶旋转90°,使得待标定剖面翼型弦线保持水平,对桨叶结构施加挥舞方向载荷m挥,分别记录光纤Bragg光栅传感器FBG1~FBG6在未加载与加载状态下的中心波长信号,计算FBG1~FBG6传感器中心波长偏移量△λ’1~△λ’6;
评估摆振解耦效果的耦合系数C计算公式为:
式中i取值为3或5,j取值为4或6,四组传感器组合方式计算得到耦合系数分别为:
比较耦合系数C1、C2、C3及C4数值大小,选择耦合系数最小的组合方式并将该传感器组合桥路偏移量作为用于标定摆振弯矩的FBG传感器波长参量,其计算公式为;
Δλ摆=Δλ1+Δλ2-(Δλi+Δλj) (10)
步骤五:采集摆振方向分级载荷作用下光纤光栅Bragg传感器FBG1~FBG6中心波长偏移量,计算摆振弯矩标定参量△λ摆对桨叶摆振弯矩进行标定;具体为:
将桨叶根部固定于标定旋转台夹具,保持桨叶翼型前缘向下,桨叶根部翼型弦线与竖直面夹角大小为剖面预扭角θ0,对桨叶结构施加摆振方向载荷m摆,桨叶摆振标定加载产生的摆振弯矩为:
M摆=m摆gL (11)
式中L为摆振载荷桨叶结构加载剖面与待标定剖面之间距离,桨叶结构剖面摆振弯矩与桨叶表面FBG1~FBG6传感器中心波长偏移量之差△λ关系式如下:
其中λ为光纤光栅Bragg传感器中心波长初始值,k为常量,h为距离,I为桨叶剖面惯性距,Pe为光纤光栅有效弹光常数。令K=I/λ(1-Pe)kh,桨叶结构剖面摆振弯矩与FBGn传感器中心波长偏移量△λn关系为:
用于标定摆振弯矩的FBG传感器波长参量△λ摆=△λ1+△λ2-(△λi+△λj),为步骤四所选择耦合系数最小的FBG传感器组合桥路波长偏移量,将式(10)代入得:
与上述试验方法相同,分别在施加摆振方向载荷2m摆、3m摆、4m摆…N·m挥加载状态下记录波长参量△λ摆,计算N次分级加载下标定系数K摆的平均值即得到桨叶叶结构摆振弯矩与FBG传感器组合桥路中心波长偏移量之差的关系方程:
Claims (1)
1.一种基于光纤应变组合桥路的桨叶弯矩解耦与标定方法,其特征在于,包括以下步骤:
步骤一:桨叶结构光纤Bragg光栅传感器网络布置;具体为:
桨叶传感网络布局分为挥舞弯矩解耦布局与摆振弯矩解耦布局,其中用于挥舞弯矩解耦标定的光纤Bragg光栅传感器FBG1和FBG2分别粘贴于桨叶结构弯矩待测剖面的上、下表面距翼型后缘1/4弦线对称位置;用于摆振弯矩解耦标定的光纤Bragg光栅传感器除FBG1和FBG2外还包括分别粘贴于桨叶结构弯矩待测剖面上、下表面距翼型前缘1~10mm弦线对称位置的FBG3和FBG5,以及分别粘贴于桨叶结构弯矩待测剖面上、下表面距翼型前缘11~20mm弦线对称位置的FBG4和FBG6;上述光纤Bragg光栅传感器粘贴方向均平行于桨叶展向以获取桨叶结构展向应变,以此构成桨叶传感器网络,光纤Bragg光栅传感器粘贴位置所对应剖面即为桨叶结构弯矩待测剖面;
步骤二:采集摆振方向载荷作用下光纤光栅Bragg传感器FBG1与FBG2中心波长偏移量,测量桨叶剖面预扭角以实现挥舞弯矩解耦;具体为:
将桨叶根部固定于标定旋转台夹具,保持桨叶翼型前缘向下,调整桨叶根部翼型弦线与竖直面夹角为θ,采用翼型夹具下方悬挂标准砝码方式对桨叶结构施加摆振方向载荷m摆,依次改变夹角θ,分别记录同加载不同夹角θ条件下光纤Bragg光栅传感器FBG1与FBG2中心波长偏移量△λ1与△λ2;
上、下表面FBG1与FBG2传感器中心波长偏移量之差与夹角θ为线性关系,表达式为:
θ=kΔλ12+b (1)
式中△λ12=△λ1-△λ2,为FBG1与FBG2中心波长偏移量之差,k为拟合系数,b为拟合曲线截距;通过上述试验测量结果,当两个光纤Bragg光栅传感器中心波长偏移量之差△λ12为0时,得到对应角度θ即拟合曲线截距b为待测桨叶预扭角θ0;
步骤三:挥舞方向分级加载作用下采集光纤光栅Bragg传感器FBG1与FBG2中心波长偏移量,计算挥舞弯矩标定参量△λ挥对桨叶挥舞弯矩进行标定;具体为:
将标定旋转台旋转90°使得桨叶根部剖面与水平面夹角大小为预扭角θ0,此时待标定剖面翼型弦线保持水平,对桨叶结构施加挥舞方向载荷m挥,桨叶挥舞标定加载产生的挥舞弯矩为:
M挥=m挥gL (2)
式中g为重力加速度,L为挥舞载荷桨叶结构加载剖面与待标定剖面之间距离,FBG1与FBG2传感器中心波长偏移量之差△λ挥=△λ1-△λ2,该偏移量之差数值仅由剖面挥舞弯矩M挥产生;桨叶结构剖面挥舞弯矩与桨叶表面FBG传感器中心波长偏移量△λ关系式如下:
其中λ为光纤光栅Bragg传感器中心波长初始值,k为常量,h为距离,I为桨叶剖面惯性距,Pe为光纤光栅有效弹光常数;令K=I/λ(1-Pe)kh,桨叶结构剖面挥舞弯矩与FBG1传感器中心波长偏移量△λ1及FBG2传感器中心波长偏移量△λ2关系分别为:
M挥=K1Δλ1,M挥=K2Δλ2 (4)
用于标定挥舞弯矩的FBG传感器波长参量△λ挥=△λ1-△λ2,其中:
与上述试验方法相同,分别在施加挥舞方向载荷2m挥、3m挥、4m挥…N·m挥加载状态下记录波长参量△λ挥,计算N次分级加载下标定系数K挥的平均值即得到桨叶叶结构挥舞弯矩与FBG1与FBG2传感器中心波长偏移量之差的关系方程:
步骤四:分别采集挥舞方向载荷与摆振方向载荷作用下光纤光栅Bragg传感器中心波长偏移量,计算耦合系数C,选择桨叶结构摆振弯矩传感器最佳解耦组合桥路;具体为:
根据步骤二所测得桨叶剖面预扭角,将桨叶根部固定于标定旋转台夹具,保持桨叶翼型前缘向下,桨叶根部翼型弦线与竖直面夹角大小为剖面预扭角θ0,对桨叶结构施加摆振方向载荷m摆,分别记录光纤Bragg光栅传感器FBG1~FBG6在未加载与加载状态下的中心波长信号,计算FBG1~FBG6传感器中心波长偏移量△λ1~△λ6;利用标定旋转台将桨叶旋转90°,使得待标定剖面翼型弦线保持水平,对桨叶结构施加挥舞方向载荷m挥,分别记录光纤Bragg光栅传感器FBG1~FBG6在未加载与加载状态下的中心波长信号,计算FBG1~FBG6传感器中心波长偏移量△λ’1~△λ’6;
评估摆振解耦效果的耦合系数C计算公式为:
式中i取值为3或5,j取值为4或6,四组传感器组合方式计算得到耦合系数分别为:
比较耦合系数C1、C2、C3及C4数值大小,选择耦合系数最小的组合方式并将该传感器组合桥路偏移量作为用于标定摆振弯矩的FBG传感器波长参量,其计算公式为;
Δλ摆=Δλ1+Δλ2-(Δλi+Δλj) (10)
步骤五:采集摆振方向分级载荷作用下光纤光栅Bragg传感器FBG1~FBG6中心波长偏移量,计算摆振弯矩标定参量△λ摆对桨叶摆振弯矩进行标定;具体为:
将桨叶根部固定于标定旋转台夹具,保持桨叶翼型前缘向下,桨叶根部翼型弦线与竖直面夹角大小为剖面预扭角θ0,对桨叶结构施加摆振方向载荷m摆,桨叶摆振标定加载产生的摆振弯矩为:
M摆=m摆gL (11)
式中g为重力加速度,L为摆振载荷桨叶结构加载剖面与待标定剖面之间距离,桨叶结构剖面摆振弯矩与桨叶表面FBG1~FBG6传感器中心波长偏移量之差△λ关系式如下:
其中λ为光纤光栅Bragg传感器中心波长初始值,k为常量,h为距离,I为桨叶剖面惯性距,Pe为光纤光栅有效弹光常数;令K=I/λ(1-Pe)kh,桨叶结构剖面摆振弯矩与FBGn传感器中心波长偏移量△λn关系为:
用于标定摆振弯矩的FBG传感器波长参量△λ摆=△λ1+△λ2-(△λi+△λj),为步骤四所选择耦合系数最小的FBG传感器组合桥路波长偏移量,将式(10)代入得:
与上述试验方法相同,分别在施加摆振方向载荷2m摆、3m摆、4m摆…N·m挥加载状态下记录波长参量△λ摆,计算N次分级加载下标定系数K摆的平均值即得到桨叶叶结构摆振弯矩与FBG传感器组合桥路中心波长偏移量之差的关系方程:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110158338.XA CN113029405B (zh) | 2021-02-04 | 2021-02-04 | 基于光纤应变组合桥路的桨叶弯矩解耦与标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110158338.XA CN113029405B (zh) | 2021-02-04 | 2021-02-04 | 基于光纤应变组合桥路的桨叶弯矩解耦与标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113029405A CN113029405A (zh) | 2021-06-25 |
CN113029405B true CN113029405B (zh) | 2022-04-22 |
Family
ID=76460078
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110158338.XA Active CN113029405B (zh) | 2021-02-04 | 2021-02-04 | 基于光纤应变组合桥路的桨叶弯矩解耦与标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113029405B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113719428B (zh) * | 2021-08-30 | 2022-11-04 | 中国华能集团清洁能源技术研究院有限公司 | 一种水平轴风力发电机组叶片寿命评估方法和系统 |
CN115014598B (zh) * | 2021-11-23 | 2023-05-23 | 中国直升机设计研究所 | 一种基于光纤光栅的旋翼主轴载荷测量方法 |
CN114754691B (zh) * | 2022-02-28 | 2023-03-31 | 南京航空航天大学 | 直升机桨叶弯扭形态分布式光纤监测与反演方法 |
CN115165175A (zh) * | 2022-04-29 | 2022-10-11 | 南京航空航天大学 | 基于光纤光栅传感器的旋翼载荷识别中的弯矩解耦方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012007005A2 (en) * | 2010-07-14 | 2012-01-19 | Vestas Wind Systems A/S | Ice detection method and system for wind turbine blades |
CN108120592A (zh) * | 2017-11-29 | 2018-06-05 | 中国直升机设计研究所 | 一种直升机桨叶静强度的试验方法 |
CN110006563A (zh) * | 2019-03-31 | 2019-07-12 | 南京航空航天大学 | 直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19842231C1 (de) * | 1998-09-15 | 2000-07-06 | Test Gmbh | Meßvorrichtung zur Bestimmung des Gesamtanzugsmoments, des Kopfreibungsmoments und der Vorspannkraft einer angezogenen Schraubverbindung |
US6915237B2 (en) * | 2002-05-14 | 2005-07-05 | Analysis And Measurement Services Corporation | Integrated system for verifying the performance and health of instruments and processes |
CN104648688B (zh) * | 2013-11-25 | 2016-11-09 | 中国直升机设计研究所 | 一种桨叶应变片布置及解耦方法 |
CN108082534B (zh) * | 2017-12-03 | 2021-03-30 | 中国直升机设计研究所 | 一种直升机旋翼复合载荷测试光纤布置和组桥方法 |
-
2021
- 2021-02-04 CN CN202110158338.XA patent/CN113029405B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012007005A2 (en) * | 2010-07-14 | 2012-01-19 | Vestas Wind Systems A/S | Ice detection method and system for wind turbine blades |
CN108120592A (zh) * | 2017-11-29 | 2018-06-05 | 中国直升机设计研究所 | 一种直升机桨叶静强度的试验方法 |
CN110006563A (zh) * | 2019-03-31 | 2019-07-12 | 南京航空航天大学 | 直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法 |
Non-Patent Citations (4)
Title |
---|
基于光纤光栅的直升机桨叶载荷测试技术研究;梁磊等;《光电子·激激光》;20191215(第12期);全文 * |
旋翼桨叶载荷测量中的解耦方法研究;余洵等;《直升机技术》;20101215(第04期);全文 * |
直升机旋翼桨叶弯矩与变形分布式光纤智能监测技术研究;陈铭杰;《中国优秀硕士学位论文全文数据库》;20200215;第32-36页 * |
直升机旋翼桨叶载荷标定新方法研究;韩玉旺等;《直升机技术》;20120915(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113029405A (zh) | 2021-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113029405B (zh) | 基于光纤应变组合桥路的桨叶弯矩解耦与标定方法 | |
CN107933957B (zh) | 一种无人直升机桨叶载荷飞行实测系统及其实测方法 | |
CN110006563B (zh) | 直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法 | |
US20090324409A1 (en) | Rotor blade monitoring | |
CN107820541B (zh) | 风力涡轮机叶片载荷传感系统 | |
CN109580057A (zh) | 基于埋入式光纤传感器的直升机旋翼载荷监测系统和方法 | |
CN108082534B (zh) | 一种直升机旋翼复合载荷测试光纤布置和组桥方法 | |
US9664506B2 (en) | High speed and high spatial density parameter measurement using fiber optic sensing technology | |
CN102937526B (zh) | 风机载荷的测量装置、系统和风机控制系统 | |
CN112796957B (zh) | 一种风机叶片的检测方法和装置以及设备 | |
CN107991060A (zh) | 基于自适应和迭代算法的载荷分布式光纤辨识方法 | |
CN110789733B (zh) | 一种直升机尾桨柔性梁挥舞变形段疲劳寿命评定方法 | |
Wu et al. | A high-frequency acceleration sensor based on fiber grating | |
James et al. | Fibre-optic measurement of strain and shape on a helicopter rotor blade during a ground run: 1. Measurement of strain | |
CN203177817U (zh) | 用于杆塔水平角度测量的带温度补偿光纤角度传感器 | |
CN201010859Y (zh) | 一种智能型桥梁拉索 | |
CN102175185A (zh) | 一种光纤光栅倾角传感器及标定方法 | |
CN112665766B (zh) | 一种基于fbg的风机叶片载荷测量方法及应用 | |
Barino et al. | Loading condition estimation using long-period fiber grating array | |
CN115752585A (zh) | 一种光纤复合传感器 | |
CN210862557U (zh) | 光纤光栅传感器装置 | |
CN100381845C (zh) | 光纤光栅门开关状态传感器探头 | |
CN113340350A (zh) | 光栅矢量传感器、架空线路覆冰弧垂状态监测装置及方法 | |
CN115126665A (zh) | 风力发电机的风机叶片监测方法及装置、存储介质、风力发电机 | |
CN205982113U (zh) | 风力发电机组的裂纹检测装置及风力发电机组 |
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 |