[go: up one dir, main page]

CN110006563A - Distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment - Google Patents

Distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment Download PDF

Info

Publication number
CN110006563A
CN110006563A CN201910254397.XA CN201910254397A CN110006563A CN 110006563 A CN110006563 A CN 110006563A CN 201910254397 A CN201910254397 A CN 201910254397A CN 110006563 A CN110006563 A CN 110006563A
Authority
CN
China
Prior art keywords
blade
section
fiber grating
strain
moment
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.)
Granted
Application number
CN201910254397.XA
Other languages
Chinese (zh)
Other versions
CN110006563B (en
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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201910254397.XA priority Critical patent/CN110006563B/en
Publication of CN110006563A publication Critical patent/CN110006563A/en
Application granted granted Critical
Publication of CN110006563B publication Critical patent/CN110006563B/en
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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开了一种直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法,包括以下步骤:计算直升机桨叶截面静矩和计算桨叶剖面惯性矩;根据基尔霍夫中性层假设原理,初步选取挥舞/摆振弯矩对应的大致解耦节点位置。在该解耦节点位置沿桨叶展向粘贴分布式光纤光栅传感器;通过对直升机桨叶施加挥舞/摆振载荷,计算得到桨叶弦向任意位置对应的应变值;将应变值进行拟合,选择该拟合曲线应变值为零时所对应的坐标位置。将该坐标位置选作为直升机桨叶挥舞/摆振弯矩对应的解耦节点。本方法能够快速准确确定直升机桨叶挥舞/摆振弯矩解耦节点关键位置,能够实现对不同形式挥舞/摆振弯矩作用下直升机桨叶不同方向应变的准确测量。

The invention discloses a distributed optical fiber decoupling measurement method for helicopter blade flapping and swaying bending moment, comprising the following steps: calculating the static moment of the helicopter blade section and calculating the blade section inertia moment; according to the Kirchhoff neutral layer Assuming the principle, preliminarily select the approximate decoupling node positions corresponding to the swaying/swing moment. The distributed fiber grating sensor is pasted along the spanwise direction of the blade at the position of the decoupling node; the strain value corresponding to any position in the chordwise direction of the blade is calculated by applying a flapping/swing vibration load to the helicopter blade; the strain value is fitted, Select the coordinate position corresponding to the zero strain value of the fitting curve. This coordinate position is selected as the decoupling node corresponding to the flapping/swing moment of the helicopter blade. The method can quickly and accurately determine the key position of the decoupling node of the helicopter blade flapping/swinging moment, and can accurately measure the strain of the helicopter blade in different directions under the action of different forms of flapping/swinging moment.

Description

直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法Distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment

技术领域technical field

本发明属于结构健康监测的弯矩监测技术领域,具体提出了一种基于基尔霍夫中性层原理的直升机桨叶挥舞与摆振弯矩解耦以及光纤传感器布局方法。The invention belongs to the technical field of bending moment monitoring for structural health monitoring, and specifically proposes a method for decoupling a helicopter blade flapping and swaying bending moment and an optical fiber sensor layout method based on the Kirchhoff neutral layer principle.

背景技术Background technique

桨叶类机械是航空航天行业的关键设备,近年来,随着研究的深入,桨叶类机械性能不断提高,不断向高速、重载、髙效、高可靠性方向发展,因此对机械系统的要求也越来越高。桨叶作为直升机的核心部件,是主要的承力对象,桨叶的损坏事故(裂纹、折断等),绝大部分是由于振动引起的,如直升机桨叶的挥舞运动、摆振运动等。因此对直升机桨叶状态进行监测是评估设计优劣、机械损伤动态监测、机械故障诊断、预测的重要方法。Blade machinery is the key equipment in the aerospace industry. In recent years, with the deepening of research, the performance of blade machinery has been continuously improved, and it has continued to develop in the direction of high speed, heavy load, high efficiency and high reliability. The requirements are also getting higher and higher. As the core component of the helicopter, the blade is the main load-bearing object. Most of the damage accidents (cracks, breaks, etc.) of the blade are caused by vibration, such as the waving motion and the swaying motion of the helicopter blade. Therefore, monitoring the state of helicopter blades is an important method to evaluate the design quality, dynamic monitoring of mechanical damage, mechanical fault diagnosis and prediction.

目前,国内外针对结构变形状态监测方法的研究主要包括光电式测量系统和非光电式测量系统,其中光电式包括摄影测量法、摄像测量法、激光跟踪仪等;非光电式包括应变传感器、加速度传感器、位移传感器、光纤传感器等。光电式测量系统一般以影像捕捉和激光扫描跟踪为主,在目标结构给定位置上布置目标点或光源,测量系统以成像的方式感知目标点或光源,根据距离修正原理、反射原理和摄影测量原理,通过坐标转换得到目标点的空间位置,响应速度快,精度高。但该类方法只适用于外形比较平整、形式比较简单的结构,抗外界干扰能力比较差,对于大变形的柔性结构,光线传输易受影响,测量系统的可靠性很难保证。非光电式测量系统主要以接触式的测量方式为主,以应变信息或曲率信息为基础,适应较复杂的结构,可以直接将传感器嵌入到结构材料中,测量性能稳定。接触式的测量方式适应较复杂结构的形变重构,测量性能比较稳定。At present, the research on structural deformation state monitoring methods at home and abroad mainly includes photoelectric measurement systems and non-photoelectric measurement systems. Among them, photoelectric methods include photogrammetry, photogrammetry, laser trackers, etc.; Sensors, displacement sensors, fiber optic sensors, etc. The photoelectric measurement system is generally based on image capture and laser scanning tracking. The target point or light source is arranged at a given position of the target structure. The measurement system perceives the target point or light source in an imaging manner. According to the principle of distance correction, reflection principle and photogrammetry The principle is to obtain the spatial position of the target point through coordinate transformation, with fast response speed and high precision. However, this type of method is only suitable for structures with relatively flat shape and simple form, and has poor anti-interference ability. For flexible structures with large deformation, the light transmission is easily affected, and the reliability of the measurement system is difficult to guarantee. Non-photoelectric measurement systems are mainly based on contact measurement methods, based on strain information or curvature information, adapt to more complex structures, and can directly embed sensors into structural materials, with stable measurement performance. The contact measurement method is suitable for the deformation and reconstruction of more complex structures, and the measurement performance is relatively stable.

因此,针对直升机桨叶工作环境恶劣,外界干扰较多的特点,选择使用非光电的接触式测量方法,并且选择相对先进且技术日趋成熟的光纤布拉格光栅传感器为测试元件,与传统传感器相比,光纤光栅传感器具有极高的灵敏度和精度,并且本征安全、抗电磁干扰强、绝缘强度高、体积小、频带宽。由于桨叶剖面的不对称性,在实际弯矩状态监测过程中挥舞弯矩和摆振弯矩会存在耦合问题,针对目前桨叶弯矩状态监测方法的不足,需要研究无需大量先验知识,能够适用于采样频率较低的常规光纤光栅解调仪,简单快速、方便可靠、实用性强的新方法。Therefore, in view of the harsh working environment of helicopter blades and the characteristics of more external interference, a non-photoelectric contact measurement method is selected, and a relatively advanced and increasingly mature fiber Bragg grating sensor is selected as the test element. Compared with traditional sensors, Fiber Bragg grating sensors have extremely high sensitivity and precision, and are intrinsically safe, strong in anti-electromagnetic interference, high in dielectric strength, small in size, and wide in frequency. Due to the asymmetry of the blade profile, there will be a coupling problem in the actual monitoring process of the bending moment and the swaying moment. In view of the shortcomings of the current monitoring method of the blade bending moment, it is necessary to study without a lot of prior knowledge. It can be applied to conventional fiber grating demodulators with lower sampling frequency, and is a new method that is simple, fast, convenient, reliable and practical.

发明内容SUMMARY OF THE INVENTION

发明目的:为了克服现有技术中存在的不足,本发明提供一种基于基尔霍夫中性层原理的直升机桨叶挥舞与摆振弯矩解耦以及光纤传感器布局方法。Purpose of the invention: In order to overcome the deficiencies in the prior art, the present invention provides a method for decoupling the flapping and swaying bending moment of a helicopter blade and an optical fiber sensor based on Kirchhoff's neutral layer principle.

技术方案:为实现上述目的,本发明采用的技术方案为:Technical scheme: In order to realize the above-mentioned purpose, the technical scheme adopted in the present invention is:

一种直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法,包括以下步骤:A distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment, comprising the following steps:

步骤一、根据直升机桨叶形状计算直升机桨叶剖面椭圆部分的静矩;Step 1. Calculate the static moment of the elliptical part of the helicopter blade section according to the shape of the helicopter blade;

步骤二、根据静力学力矩定理计算桨叶剖面形心位置,根据形心位置确定桨叶剖面主坐标系,并根据所确定坐标系计算剖面惯性矩;Step 2: Calculate the position of the centroid of the blade section according to the static moment theorem, determine the main coordinate system of the blade section according to the centroid position, and calculate the moment of inertia of the section according to the determined coordinate system;

步骤三、根据基尔霍夫中性层假设原理,初步选取挥舞/摆振弯矩对应的大致解耦节点位置,在该解耦节点位置沿桨叶展向粘贴光纤光栅传感器,光纤光栅传感器轴向方向应与桨叶展向方向一致;Step 3. According to Kirchhoff's neutral layer hypothesis, preliminarily select the approximate decoupling node position corresponding to the waving/swing bending moment, and paste the fiber grating sensor along the span of the blade at the decoupling node position, and the fiber grating sensor axis The direction should be consistent with the spanwise direction of the blade;

在桨叶结构上布置分布式光纤光栅传感器,桨叶根部固定,选取一个剖面,在经过剖面主坐标系的桨叶表面依次布置三个光纤光栅传感器,分别为光纤光栅传感器一FBG1、光纤光栅传感器二FBG2、光纤光栅传感器三FBG3,根据基尔霍夫中性层假设原理,结构中性层应力为零,而剖面主坐标系的横向轴和纵向轴近似看作是该剖面初始中性层,坐标轴与桨叶剖面的交点即为光纤光栅传感器布置点,采用光纤跳线将光纤光栅传感器一FBG1、光纤光栅传感器二FBG2、光纤光栅传感器三FBG进行串行连接,形成波分复用型传感网络;Distributed fiber grating sensors are arranged on the blade structure, the root of the blade is fixed, a section is selected, and three fiber grating sensors are sequentially arranged on the blade surface passing through the main coordinate system of the section, which are fiber grating sensor-FBG1 and fiber grating sensor respectively. 2 FBG2, fiber grating sensor 3 FBG3, according to Kirchhoff’s neutral layer hypothesis, the structural neutral layer stress is zero, and the transverse axis and longitudinal axis of the main coordinate system of the section are approximately regarded as the initial neutral layer of the section, The intersection of the coordinate axis and the blade profile is the fiber grating sensor arrangement point, and the fiber grating sensor one FBG1, the fiber grating sensor two FBG2, and the fiber grating sensor three FBGs are serially connected by using fiber jumpers to form a wavelength division multiplexing type transmission. sense network;

步骤四:对直升机桨叶进行挥舞摆振加载,测得其剖面有限点的应变值,对采集到的桨叶离散点的展向应变数据进行连续化,求得桨叶弦向方向上不同位置对应的展向应变分布曲线方程,进而反演出桨叶弦向任意位置对应的应变值;Step 4: Load the helicopter blade with swaying vibration, measure the strain value of the finite point of its profile, and continuously collect the spanwise strain data of the discrete points of the blade to obtain different positions in the chordwise direction of the blade. The corresponding spanwise strain distribution curve equation, and then invert the strain value corresponding to any position in the chordwise direction of the blade;

对桨叶施加挥舞力和摆振力,通过采集三个节点上的光纤反射波长信号,根据光纤光栅应变传感原理将其转化为应变值,再将应变值进行应变连续化,求得应变变化方程,得到该剖面任意点应变值,应变与光栅波长关系如式(6)所示Swing force and sway force are applied to the blade, and the reflected wavelength signal of the optical fiber on the three nodes is collected, which is converted into a strain value according to the principle of fiber grating strain sensing, and then the strain value is subjected to strain continuity to obtain the strain change. Equation, the strain value at any point of the section is obtained, and the relationship between strain and grating wavelength is shown in formula (6)

ΔλBB=(1-Pe)Δε (6)Δλ BB =(1-P e )Δε (6)

其中,△λB为光纤光栅中心波长变化量,λB为光纤光栅中心波长,Pe为光纤光栅有效弹光系数,△ε为应变值变化量,应变值与弦向位置关系如式(7)所示Among them, Δλ B is the change of the center wavelength of the fiber grating, λ B is the center wavelength of the fiber grating, P e is the effective elastic-optical coefficient of the fiber grating, Δε is the change of the strain value, and the relationship between the strain value and the chord position is shown in Eq. (7 ) shown

其中,△λB为光纤光栅中心波长变化量,λB为光纤光栅中心波长,Pe为光纤光栅有效弹光系数,△ε为应变值变化量,应变值与弦向位置关系如式(7)所示Among them, Δλ B is the change of the center wavelength of the fiber grating, λ B is the center wavelength of the fiber grating, P e is the effective elastic-optical coefficient of the fiber grating, Δε is the change of the strain value, and the relationship between the strain value and the chord position is shown in Eq. (7 ) shown

其中,ε(γ)为应变值,εa为弦向起始应变值,εb为弦向末端应变值,γ为弦向位置,γa为弦向起始位置,Δγ为弦向长度差;Among them, ε(γ) is the strain value, εa is the chordwise initial strain value, εb is the chordwise end strain value, γ is the chordwise position, γa is the chordwise starting position, and Δγ is the chordwise length difference ;

步骤五:将步骤四所测桨叶弦向离散位置对应的应变值进行拟合,选择该拟合曲线应变值为零时所对应的坐标位置,得到挥舞摆振解耦节点,再对挥舞摆振解耦节点位置进行修正,进而确定桨叶剖面光纤光栅传感器具体布局位置,即挥舞测量节点和摆振测量节点;Step 5: Fit the strain value corresponding to the discrete position of the blade chord direction measured in step 4, select the coordinate position corresponding to the strain value of the fitting curve when it is zero, and obtain the decoupling node of the swaying and swaying vibration. The position of the vibration decoupling node is corrected, and then the specific layout position of the blade profile fiber grating sensor is determined, that is, the swing measurement node and the swing vibration measurement node;

直升机桨叶在旋转状态下,主要受到离心力Fc以及剖面弯矩Me作用,根据力平衡可得:When the helicopter blade is rotating, it is mainly affected by the centrifugal force F c and the profile bending moment Me. According to the force balance, it can be obtained:

Fc=∫AσdA=0 (8)F c = ∫ A σdA = 0 (8)

Me=∫AzσdA (9)Me = ∫ A zσdA (9)

其中,σ为应力,A为剖面面积,z为微元面到中性轴的距离,建立桨叶剖面坐标系存在坐标关系Zσ=Z0+ΔZ;Among them, σ is the stress, A is the section area, z is the distance from the micro-element surface to the neutral axis, and there is a coordinate relationship Z σ =Z 0 +ΔZ in the blade section coordinate system;

当桨叶受到弯矩载荷作用时,一部分剖面达到极限应力σm,根据几何关系可得应力与Zσ关系:When the blade is subjected to the bending moment load, a part of the section reaches the ultimate stress σ m , and the relationship between the stress and Z σ can be obtained according to the geometric relationship:

结合坐标关系,将式(10)代入式(8)式可得Combined with the coordinate relationship, substituting Equation (10) into Equation (8) can get

替换变量可得Substitute variables to get

再将式(10)、(12)代入式(9)可得Substitute equations (10) and (12) into equation (9) to get

其中,l为剖面屈服部分到中性轴的距离,△z为中性轴偏移量,Me为剖面弯矩,Sk为桨叶剖面屈服部分对剖面中性轴的静矩,St为桨叶剖面弹性部分对剖面中性轴的静矩,It为桨叶剖面弹性区域对剖面中性轴的惯性矩,Ak为桨叶剖面屈服部分面积,At为剖面弹性部分面积,计算出中性轴偏移量△z后,对桨叶中性轴位置进行修正,最终确定FBG传感器节点位置。Among them, l is the distance from the yielding part of the profile to the neutral axis, △z is the offset of the neutral axis, Me is the bending moment of the profile, Sk is the static moment of the yielding part of the blade to the neutral axis of the profile, S t is the static moment of the elastic part of the blade section to the neutral axis of the section, I t is the moment of inertia of the elastic section of the blade section to the neutral axis of the section, A k is the area of the yield part of the blade section, A t is the area of the elastic part of the section, After the neutral axis offset △z is calculated, the neutral axis position of the blade is corrected, and the node position of the FBG sensor is finally determined.

优选的:直升机桨叶剖面椭圆部分的静矩为:Preferably: the static moment of the elliptical part of the helicopter blade section is:

其中,Sy为横向(y方向)静矩,Sz为纵向(z方向)静矩,bi,hi分别为桨叶剖面坐标系中横向、纵向椭圆公式系数,θi为桨叶剖面坐标系中椭圆相对坐标原点所占跨度,θ为角度,n为剖面所划分图形个数。Among them, S y is the transverse (y-direction) static moment, S z is the longitudinal (z-direction) static moment, bi and hi are the transverse and longitudinal ellipse formula coefficients in the blade section coordinate system, respectively , and θ i is the blade section The span occupied by the ellipse relative to the origin of the coordinate in the coordinate system, θ is the angle, and n is the number of figures divided by the section.

优选的:步骤二中剖面形心坐标y和z分别是Preferably: in step 2, the coordinates y and z of the centroid of the section are respectively

其中,y为剖面形心横向坐标位置,z为剖面形心纵向坐标位置,A为剖面面积,Ai为剖面各部分的面积,yi为剖面各部分形状的形心的横向坐标,zi为剖面各部分形状的形心的纵向坐标,Among them, y is the transverse coordinate position of the centroid of the section, z is the longitudinal coordinate position of the centroid of the section, A is the area of the section, A i is the area of each part of the section, yi is the transverse coordinate of the centroid of the shape of each section of the section, zi i is the longitudinal coordinate of the centroid of the shape of each part of the section,

确定形心位置后,以形心为原点建立桨叶剖面主坐标系,再根据主坐标系计算桨叶剖面惯性矩:After determining the position of the centroid, the main coordinate system of the blade section is established with the centroid as the origin, and then the moment of inertia of the blade section is calculated according to the main coordinate system:

I=∫Aρ2dA=∫A(y2+z2)dA=Iy+Iz (5)I=∫ A ρ 2 dA=∫ A (y 2 +z 2 )dA=I y +I z (5)

其中,I为截面总惯性矩,A为剖面面积,ρ为微面积到坐标原点的距离,Iy为横向惯性矩,Iz为纵向惯性矩,y和z分别为桨叶微元剖面的坐标值。Among them, I is the total inertia moment of the section, A is the section area, ρ is the distance from the micro-area to the coordinate origin, I y is the transverse inertia moment, I z is the longitudinal inertia moment, y and z are the coordinates of the blade micro-element section, respectively value.

优选的:当桨叶发生弯曲屈服变形时,剖面中性轴会发生偏移,因此需要确定偏移量△z,修正FBG传感器节点位置,从而提高解耦精度。Optimal: When the blade is deformed by bending and yielding, the neutral axis of the profile will be offset, so it is necessary to determine the offset Δz to correct the node position of the FBG sensor, thereby improving the decoupling accuracy.

优选的:当挥舞应变与摆振应变比值小于直升机桨叶弦展比,认为该节点处能够实现挥舞与摆振解耦;Preferably: when the ratio of flapping strain to pendulum strain is less than the chord aspect ratio of the helicopter blade, it is considered that the decoupling of flapping and pendulum vibration can be achieved at this node;

第一阶段:给直升机桨叶施加一个挥舞载荷,测得光纤光栅传感器一FBG1、光纤光栅传感器二FBG2、光纤光栅传感器三FBG3各节点的应变值分别为ε1、ε2、ε3,再给桨叶施加一个摆振载荷,测得光纤光栅传感器一FBG1、光纤光栅传感器二FBG2、光纤光栅传感器三FBG3各节点的应变值分别为ε11、ε22、ε33The first stage: apply a waving load to the helicopter blade, and measure the strain values of each node of fiber grating sensor one FBG1, fiber grating sensor two FBG2, and fiber grating sensor three FBG3 are ε 1 , ε 2 , ε 3 , and then give A swing load is applied to the blade, and the measured strain values of each node of fiber grating sensor one FBG1, fiber grating sensor two FBG2, and fiber grating sensor three FBG3 are ε 11 , ε 22 , and ε 33 respectively;

第二阶段:设定光纤光栅传感器二FBG2节点处为预估挥舞解耦点,当挥舞载荷不大于摆振载荷一半时,将摆振应变与挥舞应变之比与直升机桨叶桨叶的弦长与展长之比相比较,若则认为该点满足挥舞解耦要求;The second stage: Set the FBG2 node of the fiber grating sensor as the estimated swing decoupling point. When the swing load is not more than half of the swing load, the ratio of the swing strain to the swing strain is compared with the chord length of the helicopter blade. Compared with the length ratio, if Then it is considered that this point meets the requirements of waving decoupling;

第三阶段:设定光纤光栅传感器一FBG1和光纤光栅传感器三FBG3节点处为预估的摆振解耦点,此种情况需要比较光纤光栅传感器一FBG1和光纤光栅传感器三FBG3所在节点处分别对应的解耦精度,最终确认最佳摆振解耦点。The third stage: set the node of fiber grating sensor 1 FBG1 and fiber grating sensor 3 FBG3 as the estimated swing decoupling point. In this case, it is necessary to compare the corresponding nodes of fiber grating sensor 1 FBG1 and fiber optic sensor 3 FBG3 respectively. The best decoupling accuracy is finally confirmed.

优选的:当挥舞载荷不大于摆振载荷一半时,比较两节点处挥舞应变与摆振应变之比即的大小,将其较小者与直升机桨叶的弦长与展长之比相比较,若小于则认为该点满足摆振解耦要求。Preferably: when the swing load is not more than half of the swing load, compare the ratio of the swing strain to the swing strain at the two nodes, i.e. , compare the smaller of it with the ratio of the chord to the span of the helicopter blade, if it is less than It is considered that this point satisfies the requirement of decoupling of pendulum vibration.

本发明相比现有技术,具有以下有益效果:Compared with the prior art, the present invention has the following beneficial effects:

本方法基于基尔霍夫中性层假设原理,能够快速准确确定直升机桨叶挥舞/摆振弯矩解耦节点关键位置,无需通过在桨叶表面大量布置应变片,借助大量试验,逐步迭代确定解耦节点位置。巧妙的利用光纤光栅传感器独特的应变-方向敏感特性,能够实现对不同形式挥舞/摆振弯矩作用下直升机桨叶不同方向应变的准确测量,克服了传统电阻应变式传感器应变响应方向敏感性差的限制。此外,采用分布式光纤应变传感网络,无需像电阻应变式传感系统那样需要大量信号传输线缆,不仅有利于旋转工况下的信息传输,也避免了对被测桨叶结构固有频率等力学性能造成负面影响。Based on Kirchhoff's neutral layer assumption principle, this method can quickly and accurately determine the key positions of the helicopter blade flapping/swing moment decoupling nodes, without arranging a large number of strain gauges on the blade surface, and with the help of a large number of experiments, iteratively determined step by step Decouple node locations. Clever use of the unique strain-direction sensitivity characteristics of fiber grating sensors can accurately measure the strain of helicopter blades in different directions under the action of different forms of waving/swing bending moments, overcoming the poor sensitivity of traditional resistance strain sensors to the direction of strain response. limit. In addition, the distributed optical fiber strain sensing network does not require a large number of signal transmission cables like the resistance strain sensing system, which is not only conducive to information transmission under rotating conditions, but also avoids the need for the natural frequency of the blade structure under test. Mechanical properties are negatively affected.

附图说明Description of drawings

图1分布式光纤布拉格光栅传感器布局图;Fig. 1 Layout of distributed fiber Bragg grating sensor;

图2桨叶剖面坐标系;Figure 2 Blade section coordinate system;

具体实施方式Detailed ways

下面结合附图和具体实施例,进一步阐明本发明,应理解这些实例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。Below in conjunction with the accompanying drawings and specific embodiments, the present invention will be further clarified. It should be understood that these examples are only used to illustrate the present invention and are not used to limit the scope of the present invention. Modifications in the form of valence all fall within the scope defined by the appended claims of the present application.

一种直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法,本方法是一种基于基尔霍夫中性层原理的直升机桨叶挥舞与摆振弯矩解耦以及光纤传感器布局的方法,包括以下步骤:步骤一:推导直升机桨叶截面静矩;步骤二:计算直升机桨叶剖面的形心位置,根据形心位置确定桨叶剖面主坐标系,并根据所确定的坐标系计算桨叶剖面惯性矩;步骤三:根据基尔霍夫中性层假设原理,初步选取挥舞/摆振弯矩对应的大致解耦节点位置。在该解耦节点位置沿桨叶展向粘贴光纤光栅传感器,光纤光栅传感器轴向方向应与桨叶展向方向一致;步骤四:通过对桨叶施加挥舞/摆振载荷,测得剖面若干离散位置的应变值,对采集到的桨叶离散点的展向应变数据进行连续化,可以求得桨叶弦向方向上不同位置对应的展向应变分布曲线方程,进而可反演计算得到桨叶弦向任意位置对应的应变值;步骤五:将步骤四所测桨叶弦向若干离散位置对应的应变值进行拟合,选择该拟合曲线应变值为零时所对应的坐标位置。将该坐标位置选作为直升机桨叶挥舞/摆振弯矩对应的解耦节点。本方法基于基尔霍夫中性层假设原理,能够快速准确确定直升机桨叶挥舞/摆振弯矩解耦节点关键位置,无需通过在桨叶表面大量布置应变片,借助大量试验,逐步迭代确定解耦节点位置。巧妙的利用光纤光栅传感器独特的应变-方向敏感特性,能够实现对不同形式挥舞/摆振弯矩作用下直升机桨叶不同方向应变的准确测量,克服了传统电阻应变式传感器应变响应方向敏感性差的限制。此外,采用分布式光纤应变传感网络,无需像电阻应变式传感系统那样需要大量信号传输线缆,不仅有利于旋转工况下的信息传输,也避免了对被测桨叶结构固有频率等力学性能造成负面影响。A distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment. The method includes the following steps: step 1: deriving the static moment of the helicopter blade section; step 2: calculating the centroid position of the helicopter blade section, determining the main coordinate system of the blade section according to the centroid position, and calculating according to the determined coordinate system Blade profile moment of inertia; Step 3: According to Kirchhoff's neutral layer assumption, preliminarily select the approximate decoupling node positions corresponding to the flapping/swing moment. The fiber grating sensor is pasted along the spanwise direction of the blade at the position of the decoupling node, and the axial direction of the fiber grating sensor should be consistent with the spanwise direction of the blade; Step 4: By applying a waving/swinging load to the blade, a number of discrete cross sections are measured. The strain value of the position, and the collected spanwise strain data of discrete points of the blade are continuous, the spanwise strain distribution curve equation corresponding to different positions in the chordwise direction of the blade can be obtained, and then the blade can be obtained by inversion calculation. The strain value corresponding to any position in the chord direction; Step 5: Fit the strain values corresponding to several discrete positions in the chord direction of the blade measured in step 4, and select the coordinate position corresponding to the fitting curve when the strain value is zero. This coordinate position is selected as the decoupling node corresponding to the flapping/swing moment of the helicopter blade. Based on Kirchhoff's neutral layer assumption principle, this method can quickly and accurately determine the key positions of the helicopter blade flapping/swing moment decoupling nodes, without arranging a large number of strain gauges on the blade surface, and with the help of a large number of experiments, iteratively determined step by step Decouple node locations. Clever use of the unique strain-direction sensitivity characteristics of fiber grating sensors can accurately measure the strain of helicopter blades in different directions under the action of different forms of waving/swing bending moments, overcoming the poor sensitivity of traditional resistance strain sensors to the direction of strain response. limit. In addition, the distributed optical fiber strain sensing network does not require a large number of signal transmission cables like the resistance strain sensing system, which is not only conducive to information transmission under rotating conditions, but also avoids the need for the natural frequency of the blade structure under test. Mechanical properties are negatively affected.

具体包括以下步骤:Specifically include the following steps:

步骤一、推导直升机桨叶截面静矩,采用碳纤维材料的直升机旋翼桨叶缩比模型;直升机桨叶的典型剖面一般可以看成是椭圆,三角形等二维图形的结合。其中,直升机桨叶剖面椭圆部分的静矩为:Step 1. Derive the static moment of the helicopter blade section, and use the carbon fiber material helicopter rotor blade scale model; the typical section of the helicopter blade can generally be regarded as a combination of two-dimensional figures such as ellipses and triangles. Among them, the static moment of the elliptical part of the helicopter blade section is:

其中,Sy为横向(y轴方向)静矩,Sz为纵向(z轴方向)静矩,bi,hi分别为桨叶剖面坐标系中椭圆公式系数,θi为桨叶剖面坐标系中椭圆相对坐标原点所占跨度。Among them, S y is the transverse (y-axis direction) static moment, S z is the longitudinal (z-axis direction) static moment, b i and hi are the ellipse formula coefficients in the blade section coordinate system, respectively, and θ i is the blade section coordinate The span occupied by the ellipse relative to the origin of the coordinates in the system.

步骤二、计算桨叶剖面形心位置,根据形心位置确定桨叶剖面主坐标系,并根据所确定坐标系计算剖面惯性矩;Step 2: Calculate the centroid position of the blade section, determine the main coordinate system of the blade section according to the centroid position, and calculate the moment of inertia of the section according to the determined coordinate system;

根据静力学力矩定理,剖面形心坐标y和z分别是According to the static moment theorem, the centroid coordinates y and z of the section are respectively

确定形心位置后,即可以形心为原点建立桨叶剖面主坐标系,再根据主坐标系计算桨叶剖面惯性矩After determining the position of the centroid, the main coordinate system of the blade section can be established with the centroid as the origin, and then the moment of inertia of the blade section is calculated according to the main coordinate system.

I=∫Aρ2dA=∫A(y2+z2)dA=Iy+Iz (5)I=∫ A ρ 2 dA=∫ A (y 2 +z 2 )dA=I y +I z (5)

即剖面惯性矩为剖面各部分惯性矩之和。其中,I为截面总惯性矩,Iy为横向(y轴方向)惯性矩,Iz为纵向(z轴方向)惯性矩,ρ为微面积到坐标原点的距离,A为剖面面积,y和z分别为桨叶微元剖面的坐标值。That is, the moment of inertia of the section is the sum of the moments of inertia of each part of the section. Among them, I is the total moment of inertia of the section, I y is the transverse (y-axis direction) moment of inertia, I z is the longitudinal (z-axis direction) moment of inertia, ρ is the distance from the micro area to the coordinate origin, A is the section area, y and z are the coordinate values of the blade micro-element profile, respectively.

步骤三、根据基尔霍夫中性层假设原理,初步选取挥舞/摆振弯矩对应的大致解耦节点位置。在该解耦节点位置沿桨叶展向粘贴光纤光栅传感器,光纤光栅传感器轴向方向应与桨叶展向方向一致;Step 3: According to Kirchhoff's neutral layer hypothesis, preliminarily select the approximate decoupling node positions corresponding to the waving/swing bending moment. Paste the fiber grating sensor along the spanwise direction of the blade at the position of the decoupling node, and the axial direction of the fiber grating sensor should be consistent with the spanwise direction of the blade;

如图1所示,在桨叶结构上布置分布式光纤光栅传感器,桨叶根部固定,选取一个剖面,在经过剖面主坐标系的桨叶表面依次布置三个光纤光栅传感器。根据基尔霍夫中性层假设原理,结构中性层应力为零,而剖面主坐标系的x轴、y轴可以近似看作是该剖面初始中性层,坐标轴与桨叶剖面的交点即为光纤光栅传感器布置点,采用光纤跳线将FBG1、FBG2、FBG3进行串行连接,形成波分复用型传感网络,传感网络所使用的光纤光栅解调仪采样频率为1KHz。As shown in Figure 1, a distributed fiber grating sensor is arranged on the blade structure, the root of the blade is fixed, a section is selected, and three fiber grating sensors are sequentially arranged on the blade surface passing through the main coordinate system of the section. According to Kirchhoff's neutral layer hypothesis, the structural neutral layer stress is zero, and the x-axis and y-axis of the main coordinate system of the profile can be approximately regarded as the initial neutral layer of the profile, and the intersection of the coordinate axis and the blade profile It is the arrangement point of the fiber grating sensor. The FBG1, FBG2, and FBG3 are serially connected by fiber jumpers to form a wavelength division multiplexing sensor network. The sampling frequency of the fiber grating demodulator used in the sensor network is 1KHz.

步骤四:对直升机桨叶进行挥舞摆振加载,测得其剖面有限点的应变值,对采集到的桨叶离散点的展向应变数据进行连续化,可以求得桨叶弦向方向上不同位置对应的展向应变分布曲线方程,进而可反演出桨叶弦向任意位置对应的应变值;Step 4: Load the helicopter blade with swaying vibration, measure the strain value at the finite point of its profile, and serialize the collected spanwise strain data at discrete points of the blade, so that the difference in the chordwise direction of the blade can be obtained. The spanwise strain distribution curve equation corresponding to the position, and then the strain value corresponding to any position in the chordwise direction of the blade can be inverted;

对桨叶施加挥舞力和摆振力,通过采集三个节点上的光纤反射波长信号,根据光纤光栅应变传感原理将其转化为应变值,再将应变值进行应变连续化,求得应变变化方程,得到该剖面任意点应变值,应变与光栅波长关系如式(6)所示Swing force and sway force are applied to the blade, and the reflected wavelength signal of the optical fiber on the three nodes is collected, which is converted into a strain value according to the principle of fiber grating strain sensing, and then the strain value is subjected to strain continuity to obtain the strain change. Equation, the strain value at any point of the section is obtained, and the relationship between strain and grating wavelength is shown in formula (6)

ΔλBB=(1-Pe)Δε (6)Δλ BB =(1-P e )Δε (6)

其中,△λB为光纤光栅中心波长变化量,λB为光纤光栅中心波长,Pe为光纤光栅有效弹光系数,△ε为应变值。应变值与弦向位置关系如式(7)所示Among them, Δλ B is the variation of the center wavelength of the fiber grating, λ B is the center wavelength of the fiber grating, P e is the effective elastic-optical coefficient of the fiber grating, and Δε is the strain value. The relationship between strain value and chord position is shown in formula (7)

步骤五:将步骤四所测桨叶弦向若干离散位置对应的应变值进行拟合,选择该拟合曲线应变值为零时所对应的坐标位置,得到挥舞摆振解耦节点。再对节点位置进行修正,进而可以确定桨叶剖面光纤光栅传感器具体布局位置,即挥舞测量节点和摆振测量节点。Step 5: Fit the strain values corresponding to several discrete positions in the blade chord direction measured in step 4, and select the coordinate position corresponding to the fitting curve when the strain value is zero to obtain the swinging vibration decoupling node. Then the node position is corrected, and then the specific layout position of the blade profile fiber grating sensor can be determined, that is, the swing measurement node and the sway measurement node.

其中,当桨叶发生弯曲屈服变形时,剖面中性轴会发生偏移,因此需要确定偏移量△z,修正FBG传感器节点位置,从而提高解耦精度。Among them, when the blade is deformed by bending and yielding, the neutral axis of the profile will be offset, so it is necessary to determine the offset △z and correct the node position of the FBG sensor, thereby improving the decoupling accuracy.

直升机桨叶在旋转状态下,主要受到离心力Fc以及弯矩Me作用,根据力平衡可得In the rotating state of the helicopter blade, it is mainly affected by the centrifugal force F c and the bending moment Me. According to the force balance, it can be obtained:

Fc=∫AσdA=0 (8)F c = ∫ A σdA = 0 (8)

Me=∫AzσdA (9)Me = ∫ A zσdA (9)

其中,σ为应力,A为剖面面积,z为微元面到中性轴的距离。建立桨叶剖面坐标系如图2所示,存在坐标关系Zσ=Z0+ΔZ。Among them, σ is the stress, A is the cross-sectional area, and z is the distance from the element surface to the neutral axis. The blade section coordinate system is established as shown in Figure 2, and there is a coordinate relationship Z σ =Z 0 +ΔZ.

当桨叶受到弯矩载荷作用时,一部分剖面达到极限应力σm,根据几何关系可得应力与Zσ关系When the blade is subjected to the bending moment load, a part of the section reaches the ultimate stress σ m , and the relationship between the stress and Z σ can be obtained according to the geometric relationship

结合坐标关系,将式(10)代入式(8)式可得Combined with the coordinate relationship, substituting Equation (10) into Equation (8) can get

替换变量可得Substitute variables to get

再将式(10)、(12)代入式(9)可得Substitute equations (10) and (12) into equation (9) to get

其中,l为剖面屈服部分到中性轴的距离,△z为中性轴偏移量,Me为剖面弯矩,Sk为桨叶剖面屈服部分对剖面中性轴的静矩,St为桨叶剖面弹性部分对剖面中性轴的静矩,It为桨叶剖面弹性区域对剖面中性轴的惯性矩,Ak为桨叶剖面屈服部分面积,At为剖面弹性部分面积。计算出中性轴偏移量△z后,对桨叶中性轴位置进行修正,最终确定FBG传感器节点位置。Among them, l is the distance from the yielding part of the profile to the neutral axis, △z is the offset of the neutral axis, Me is the bending moment of the profile, Sk is the static moment of the yielding part of the blade to the neutral axis of the profile, S t is the static moment of the elastic part of the blade section to the neutral axis of the section, I t is the moment of inertia of the elastic section of the blade section to the neutral axis of the section, Ak is the area of the yield part of the blade section, and At is the area of the elastic part of the section. After the neutral axis offset △z is calculated, the neutral axis position of the blade is corrected, and the node position of the FBG sensor is finally determined.

当桨叶受摆振力时,若节点处挥舞应变与摆振应变的比值为0,则认为摆振力解耦。但这一要求在实际工程中往往较难达到,因此,一般当挥舞应变与摆振应变比值小于直升机桨叶弦展比,即可认为该节点处能够实现挥舞与摆振解耦。When the blade is subjected to the sway force, if the ratio of the flapping strain to the sway strain at the node is 0, the sway force is considered to be decoupled. However, this requirement is often difficult to achieve in practical engineering. Therefore, when the ratio of the flapping strain to the swaying strain is less than the chord aspect ratio of the helicopter blade, it can be considered that the decoupling of flapping and swaying can be achieved at this node.

第一阶段:给桨叶施加一个挥舞载荷,测得FBG1、FBG2、FBG3各节点的应变值分别为ε1、ε2、ε3。再给桨叶施加一个摆振载荷,测得FBG1、FBG2、FBG3各节点的应变值分别为ε11、ε22、ε33The first stage: a swing load is applied to the blade, and the measured strain values of each node of FBG1, FBG2, and FBG3 are ε 1 , ε 2 , and ε 3 , respectively. A swing load is then applied to the blade, and the measured strain values of each node of FBG1, FBG2, and FBG3 are ε 11 , ε 22 , and ε 33 , respectively;

第二阶段:设定FBG2节点处为预估挥舞解耦点,当挥舞载荷不大于摆振载荷一半时,将摆振应变与挥舞应变之比与直升机桨叶的弦长与展长之比相比较,若则认为该点满足挥舞解耦要求;The second stage: Set the FBG2 node as the estimated swing decoupling point. When the swing load is not more than half of the swing load, compare the ratio of the swing strain to the swing strain with the ratio of the chord to the extension of the helicopter blade compare, if Then it is considered that this point meets the requirements of waving decoupling;

第三阶段:设定FBG1和FBG3节点处为预估的摆振解耦点,此种情况需要比较FBG1和FBG3所在节点处分别对应的解耦精度,最终确认最佳摆振解耦点。当挥舞载荷不大于摆振载荷一半时,比较两节点处挥舞应变与摆振应变之比即的大小,将其较小者与直升机桨叶的弦长与展长之比相比较,若小于则认为该点满足摆振解耦要求。The third stage: Set the nodes FBG1 and FBG3 as the estimated swing decoupling points. In this case, it is necessary to compare the decoupling accuracy corresponding to the nodes where FBG1 and FBG3 are located, and finally confirm the optimal swing decoupling point. When the swing load is not more than half of the swing load, compare the ratio of swing strain to swing strain at the two nodes, i.e. , compare the smaller of it with the ratio of the chord to the span of the helicopter blade, if it is less than It is considered that this point satisfies the requirement of decoupling of pendulum vibration.

以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。The above is only the preferred embodiment of the present invention, it should be pointed out that: for those skilled in the art, without departing from the principle of the present invention, several improvements and modifications can also be made, and these improvements and modifications are also It should be regarded as the protection scope of the present invention.

Claims (7)

1.一种直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法,其特征在于,包括以下步骤:1. a method for measuring the distributed optical fiber decoupling of helicopter blade waving and sway bending moment, is characterized in that, comprises the following steps: 步骤一、根据直升机桨叶形状计算直升机桨叶剖面椭圆部分的静矩;Step 1. Calculate the static moment of the elliptical part of the helicopter blade section according to the shape of the helicopter blade; 步骤二、根据静力学力矩定理计算桨叶剖面形心位置,根据形心位置确定桨叶剖面主坐标系,并根据所确定坐标系计算剖面惯性矩;Step 2: Calculate the position of the centroid of the blade section according to the static moment theorem, determine the main coordinate system of the blade section according to the centroid position, and calculate the moment of inertia of the section according to the determined coordinate system; 步骤三、根据基尔霍夫中性层假设原理,初步选取挥舞/摆振弯矩对应的大致解耦节点位置,在该解耦节点位置沿桨叶展向粘贴光纤光栅传感器,光纤光栅传感器轴向方向应与桨叶展向方向一致;Step 3. According to Kirchhoff's neutral layer hypothesis, preliminarily select the approximate decoupling node position corresponding to the waving/swing bending moment, and paste the fiber grating sensor along the span of the blade at the decoupling node position, and the fiber grating sensor axis The direction should be consistent with the spanwise direction of the blade; 在桨叶结构上布置分布式光纤光栅传感器,桨叶根部固定,选取一个剖面,在经过剖面主坐标系的桨叶表面依次布置三个光纤光栅传感器,分别为光纤光栅传感器一FBG1、光纤光栅传感器二FBG2、光纤光栅传感器三FBG3,根据基尔霍夫中性层假设原理,结构中性层应力为零,而剖面主坐标系的横向轴和纵向轴近似看作是该剖面初始中性层,坐标轴与桨叶剖面的交点即为光纤光栅传感器布置点,采用光纤跳线将光纤光栅传感器一FBG1、光纤光栅传感器二FBG2、光纤光栅传感器三FBG进行串行连接,形成波分复用型传感网络;Distributed fiber grating sensors are arranged on the blade structure, the root of the blade is fixed, a section is selected, and three fiber grating sensors are sequentially arranged on the blade surface passing through the main coordinate system of the section, which are fiber grating sensor-FBG1 and fiber grating sensor respectively. 2 FBG2, fiber grating sensor 3 FBG3, according to Kirchhoff’s neutral layer hypothesis, the structural neutral layer stress is zero, and the transverse axis and longitudinal axis of the main coordinate system of the section are approximately regarded as the initial neutral layer of the section, The intersection of the coordinate axis and the blade profile is the fiber grating sensor arrangement point, and the fiber grating sensor one FBG1, the fiber grating sensor two FBG2, and the fiber grating sensor three FBGs are serially connected by using fiber jumpers to form a wavelength division multiplexing type transmission. sense network; 步骤四:对直升机桨叶进行挥舞摆振加载,测得其剖面有限点的应变值,对采集到的桨叶离散点的展向应变数据进行连续化,求得桨叶弦向方向上不同位置对应的展向应变分布曲线方程,进而反演出桨叶弦向任意位置对应的应变值;Step 4: Load the helicopter blade with swaying vibration, measure the strain value of the finite point of its profile, and continuously collect the spanwise strain data of the discrete points of the blade to obtain different positions in the chordwise direction of the blade. The corresponding spanwise strain distribution curve equation, and then invert the strain value corresponding to any position in the chordwise direction of the blade; 对桨叶施加挥舞力和摆振力,通过采集三个节点上的光纤反射波长信号,根据光纤光栅应变传感原理将其转化为应变值,再将应变值进行应变连续化,求得应变变化方程,得到该剖面任意点应变值,应变与光栅波长关系如式(6)所示Swing force and sway force are applied to the blade, and the reflected wavelength signal of the optical fiber on the three nodes is collected, which is converted into a strain value according to the principle of fiber grating strain sensing, and then the strain value is subjected to strain continuity to obtain the strain change. Equation, the strain value at any point of the section is obtained, and the relationship between strain and grating wavelength is shown in formula (6) ΔλBB=(1-Pe)Δε (6)Δλ BB =(1-P e )Δε (6) 其中,△λB为光纤光栅中心波长变化量,λB为光纤光栅中心波长,Pe为光纤光栅有效弹光系数,△ε为应变值变化量,应变值与弦向位置关系如式(7)所示Among them, Δλ B is the change of the center wavelength of the fiber grating, λ B is the center wavelength of the fiber grating, P e is the effective elastic-optical coefficient of the fiber grating, Δε is the change of the strain value, and the relationship between the strain value and the chord position is shown in Eq. (7 ) shown 其中,ε(γ)为应变值,εa为弦向起始应变值,εb为弦向末端应变值,γ为弦向位置,γa为弦向起始位置,Δγ为弦向长度差;Among them, ε(γ) is the strain value, εa is the chordwise initial strain value, εb is the chordwise end strain value, γ is the chordwise position, γa is the chordwise starting position, and Δγ is the chordwise length difference ; 步骤五:将步骤四所测桨叶弦向离散位置对应的应变值进行拟合,选择该拟合曲线应变值为零时所对应的坐标位置,得到挥舞摆振解耦节点,再对挥舞摆振解耦节点位置进行修正,进而确定桨叶剖面光纤光栅传感器具体布局位置,即挥舞测量节点和摆振测量节点;Step 5: Fit the strain value corresponding to the discrete position of the blade chord direction measured in step 4, select the coordinate position corresponding to the strain value of the fitting curve when it is zero, and obtain the decoupling node of the swaying and swaying vibration. The position of the vibration decoupling node is corrected, and then the specific layout position of the blade profile fiber grating sensor is determined, that is, the swing measurement node and the swing vibration measurement node; 直升机桨叶在旋转状态下,主要受到离心力Fc以及剖面弯矩Me作用,根据力平衡可得:When the helicopter blade is rotating, it is mainly affected by the centrifugal force F c and the profile bending moment Me. According to the force balance, it can be obtained: Fc=∫AσdA=0 (8)F c = ∫ A σdA = 0 (8) Me=∫AzσdA (9)Me = ∫ A zσdA (9) 其中,σ为应力,A为剖面面积,z为微元面到中性轴的距离,建立桨叶剖面坐标系存在坐标关系Zσ=Z0+ΔZ;Among them, σ is the stress, A is the section area, z is the distance from the micro-element surface to the neutral axis, and there is a coordinate relationship Z σ =Z 0 +ΔZ in the blade section coordinate system; 当桨叶受到弯矩载荷作用时,一部分剖面达到极限应力σm,根据几何关系可得应力与Zσ关系:When the blade is subjected to the bending moment load, a part of the section reaches the ultimate stress σ m , and the relationship between the stress and Z σ can be obtained according to the geometric relationship: 结合坐标关系,将式(10)代入式(8)式可得Combined with the coordinate relationship, substituting Equation (10) into Equation (8) can get 替换变量可得Substitute variables to get 再将式(10)、(12)代入式(9)可得Substitute equations (10) and (12) into equation (9) to get 其中,l为剖面屈服部分到中性轴的距离,△z为中性轴偏移量,Me为剖面弯矩,Sk为桨叶剖面屈服部分对剖面中性轴的静矩,St为桨叶剖面弹性部分对剖面中性轴的静矩,It为桨叶剖面弹性区域对剖面中性轴的惯性矩,Ak为桨叶剖面屈服部分面积,At为剖面弹性部分面积,计算出中性轴偏移量△z后,对桨叶中性轴位置进行修正,最终确定FBG传感器节点位置。Among them, l is the distance from the yielding part of the profile to the neutral axis, △z is the offset of the neutral axis, Me is the bending moment of the profile, Sk is the static moment of the yielding part of the blade to the neutral axis of the profile, S t is the static moment of the elastic part of the blade section to the neutral axis of the section, I t is the moment of inertia of the elastic section of the blade section to the neutral axis of the section, A k is the area of the yield part of the blade section, A t is the area of the elastic part of the section, After the neutral axis offset △z is calculated, the neutral axis position of the blade is corrected, and the node position of the FBG sensor is finally determined. 2.根据权利要求1所述直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法,其特征在于:直升机桨叶剖面椭圆部分的静矩为:2. according to the described helicopter blade of claim 1, it is characterized in that: the static moment of the helicopter blade section ellipse part is: 其中,Sy为横向静矩,Sz为纵向静矩,bi,hi分别为桨叶剖面坐标系中横向、纵向椭圆公式系数,θi为桨叶剖面坐标系中椭圆相对坐标原点所占跨度,θ为角度,n为剖面所划分图形个数。Among them, S y is the transverse static moment, S z is the longitudinal static moment, b i , hi are the transverse and longitudinal ellipse formula coefficients in the blade section coordinate system, respectively, θ i is the relative coordinate origin of the ellipse in the blade section coordinate system Occupy the span, θ is the angle, and n is the number of figures divided by the section. 3.根据权利要求2所述直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法,其特征在于:步骤二中剖面形心坐标y和z分别是3. according to the described helicopter blade waving and swaying bending moment distributed optical fiber decoupling measurement method of claim 2, it is characterized in that: in step 2, the profile centroid coordinates y and z are respectively 其中,y为剖面形心横向坐标位置,z为剖面形心纵向坐标位置,A为剖面面积,Ai为剖面各部分的面积,yi为剖面各部分形状的形心的横向坐标,zi为剖面各部分形状的形心的纵向坐标,Among them, y is the transverse coordinate position of the centroid of the section, z is the longitudinal coordinate position of the centroid of the section, A is the area of the section, A i is the area of each part of the section, yi is the transverse coordinate of the centroid of the shape of each section of the section, zi i is the longitudinal coordinate of the centroid of the shape of each part of the section, 确定形心位置后,以形心为原点建立桨叶剖面主坐标系,再根据主坐标系计算桨叶剖面惯性矩:After determining the position of the centroid, the main coordinate system of the blade section is established with the centroid as the origin, and then the moment of inertia of the blade section is calculated according to the main coordinate system: I=∫Aρ2dA=∫A(y2+z2)dA=Iy+Iz (5)I=∫ A ρ 2 dA=∫ A (y 2 +z 2 )dA=I y +I z (5) 其中,I为截面总惯性矩,A为剖面面积,ρ为微面积到坐标原点的距离,Iy为横向惯性矩,Iz为纵向惯性矩。Among them, I is the total moment of inertia of the section, A is the area of the section, ρ is the distance from the micro-area to the origin of the coordinates, I y is the transverse moment of inertia, and I z is the longitudinal moment of inertia. 4.根据权利要求3所述直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法,其特征在于:当桨叶发生弯曲屈服变形时,剖面中性轴会发生偏移,因此需要确定偏移量△z,修正FBG传感器节点位置,从而提高解耦精度。4. The distributed optical fiber decoupling measurement method of helicopter blade waving and sway bending moment according to claim 3, characterized in that: when the blade undergoes bending and yielding deformation, the neutral axis of the profile will be offset, so it is necessary to determine The offset △z corrects the node position of the FBG sensor, thereby improving the decoupling accuracy. 5.根据权利要求4所述直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法,其特征在于:当挥舞应变与摆振应变比值小于直升机桨叶弦展比,认为该节点处能够实现挥舞与摆振解耦;5. The method for measuring the distributed optical fiber decoupling between the swinging and swaying bending moment of a helicopter blade according to claim 4, is characterized in that: when the ratio of swinging strain and swinging strain is less than the chord aspect ratio of the helicopter blade, it is considered that the node can be Realize the decoupling of swing and swing; 第一阶段:给桨叶施加一个挥舞载荷,测得光纤光栅传感器一FBG1、光纤光栅传感器二FBG2、光纤光栅传感器三FBG3各节点的应变值分别为ε1、ε2、ε3,再给桨叶施加一个摆振载荷,测得光纤光栅传感器一FBG1、光纤光栅传感器二FBG2、光纤光栅传感器三FBG3各节点的应变值分别为ε11、ε22、ε33The first stage: apply a waving load to the propeller, and measure the strain values of each node of fiber grating sensor 1 FBG1, fiber grating sensor 2 FBG2, and fiber grating sensor 3 FBG3, respectively ε 1 , ε 2 , ε 3 , and then apply the propeller A swing load is applied to the blade, and the measured strain values of each node of fiber grating sensor 1 FBG1, fiber grating sensor 2 FBG2, and fiber grating sensor 3 FBG3 are ε 11 , ε 22 , ε 33 , respectively; 第二阶段:设定光纤光栅传感器二FBG2节点处为预估挥舞解耦点,当挥舞载荷不大于摆振载荷一半时,将摆振应变与挥舞应变之比与直升机桨叶的弦长与展长之比相比较,若则认为该点满足挥舞解耦要求;The second stage: Set the FBG2 node of the fiber grating sensor as the estimated swing decoupling point. When the swing load is not more than half of the swing load, the ratio of the swing strain to the swing strain is compared with the chord length and spread of the helicopter blade. Compared with the length ratio, if Then it is considered that this point meets the requirements of waving decoupling; 第三阶段:设定光纤光栅传感器一FBG1和光纤光栅传感器三FBG3节点处为预估的摆振解耦点,此种情况需要比较光纤光栅传感器一FBG1和光纤光栅传感器三FBG3所在节点处分别对应的解耦精度,最终确认最佳摆振解耦点。The third stage: set the node of fiber grating sensor 1 FBG1 and fiber grating sensor 3 FBG3 as the estimated swing decoupling point. In this case, it is necessary to compare the corresponding nodes of fiber grating sensor 1 FBG1 and fiber optic sensor 3 FBG3 respectively. The best decoupling accuracy is finally confirmed. 6.根据权利要求5所述直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法,其特征在于:当挥舞载荷不大于摆振载荷一半时,比较两节点处挥舞应变与摆振应变之比即的大小,将其较小者与直升机桨叶的弦长与展长之比相比较,若小于则认为该点满足摆振解耦要求。6. The method for measuring the distributed optical fiber decoupling of helicopter blade flapping and swaying bending moment according to claim 5, characterized in that: when the flapping load is not more than half of the swaying load, compare the flapping strain and the swaying strain at the two nodes The ratio is , compare the smaller of it with the ratio of the chord to the span of the helicopter blade, if it is less than It is considered that this point satisfies the requirement of decoupling of pendulum vibration. 7.根据权利要求6所述直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法,其特征在于:传感网络所使用的光纤光栅解调仪采样频率为1KHz。7 . The distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment according to claim 6 , wherein the sampling frequency of the fiber grating demodulator used in the sensing network is 1KHz. 8 .
CN201910254397.XA 2019-03-31 2019-03-31 Distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment Active CN110006563B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910254397.XA CN110006563B (en) 2019-03-31 2019-03-31 Distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910254397.XA CN110006563B (en) 2019-03-31 2019-03-31 Distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment

Publications (2)

Publication Number Publication Date
CN110006563A true CN110006563A (en) 2019-07-12
CN110006563B CN110006563B (en) 2020-09-18

Family

ID=67169090

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910254397.XA Active CN110006563B (en) 2019-03-31 2019-03-31 Distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment

Country Status (1)

Country Link
CN (1) CN110006563B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110595373A (en) * 2019-09-06 2019-12-20 沈阳天眼智云信息科技有限公司 Power transformer winding structure strain detection method based on optical fiber sensing
CN111289020A (en) * 2020-03-05 2020-06-16 山东大学 A Perception Method for Large Deformation Spatial Pose of Flexible Body
CN111677635A (en) * 2020-04-30 2020-09-18 中国电力科学研究院有限公司 A method and device for measuring the bending moment of wind turbine blades considering orthogonal effects
CN112432746A (en) * 2020-10-21 2021-03-02 中国人民解放军总参谋部第六十研究所 Method for determining swing load coupling degree of helicopter blade
CN112525066A (en) * 2020-11-20 2021-03-19 中国直升机设计研究所 Helicopter blade surface strain gauge pasting and wire wiring method
CN113029405A (en) * 2021-02-04 2021-06-25 南京航空航天大学 Blade bending moment decoupling and calibrating method based on optical fiber strain combined bridge circuit
CN114166141A (en) * 2021-11-23 2022-03-11 中国直升机设计研究所 A multi-circle marker detection method for helicopter high-speed rotating blade images
CN115031647A (en) * 2022-04-12 2022-09-09 上海隧道工程有限公司 Deformation detection method for end plates
CN115165175A (en) * 2022-04-29 2022-10-11 南京航空航天大学 Bending moment decoupling method in rotor load identification based on fiber grating sensor
CN119223168A (en) * 2024-12-05 2024-12-31 南京信息职业技术学院 Determination method of helicopter blade flapping and shimmy decoupling position using fiber Bragg grating sensor
CN119468956A (en) * 2025-01-14 2025-02-18 中国航空工业集团公司哈尔滨空气动力研究所 A method for measuring blade spacing in wind tunnel tests of coaxial rigid rotor models

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008131303A3 (en) * 2007-04-20 2009-02-05 Hansen Medical Inc Optical fiber shape sensing systems
EP2351169A1 (en) * 2008-10-24 2011-08-03 Cymer, Inc. System method and apparatus for selecting and controlling light source bandwidth
CN102426457A (en) * 2011-11-20 2012-04-25 西北工业大学 Flight control navigation system of miniature flapping-wing aircraft
CN103010485A (en) * 2012-12-21 2013-04-03 南京航空航天大学 Simulation modeling method for tilt-rotor unmanned plane and system thereof
CN104648688A (en) * 2013-11-25 2015-05-27 中国直升机设计研究所 Paddle strain gauge arranging and decoupling method
EP2915752A1 (en) * 2012-11-05 2015-09-09 Kabushiki Kaisha Towani Aeroplane crushing device, and aeroplane demolition method
CN108120592A (en) * 2017-11-29 2018-06-05 中国直升机设计研究所 A kind of test method of helicopter blade static strength
CN108387183A (en) * 2018-01-31 2018-08-10 天津大学 Servo-actuated real-time measurement apparatus and its method for rotating vane whole audience dynamic deformation
US20180372778A1 (en) * 2017-06-26 2018-12-27 Corad Technology Inc. Detachable probe card interface

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008131303A3 (en) * 2007-04-20 2009-02-05 Hansen Medical Inc Optical fiber shape sensing systems
EP2351169A1 (en) * 2008-10-24 2011-08-03 Cymer, Inc. System method and apparatus for selecting and controlling light source bandwidth
CN102426457A (en) * 2011-11-20 2012-04-25 西北工业大学 Flight control navigation system of miniature flapping-wing aircraft
EP2915752A1 (en) * 2012-11-05 2015-09-09 Kabushiki Kaisha Towani Aeroplane crushing device, and aeroplane demolition method
CN103010485A (en) * 2012-12-21 2013-04-03 南京航空航天大学 Simulation modeling method for tilt-rotor unmanned plane and system thereof
CN104648688A (en) * 2013-11-25 2015-05-27 中国直升机设计研究所 Paddle strain gauge arranging and decoupling method
US20180372778A1 (en) * 2017-06-26 2018-12-27 Corad Technology Inc. Detachable probe card interface
CN108120592A (en) * 2017-11-29 2018-06-05 中国直升机设计研究所 A kind of test method of helicopter blade static strength
CN108387183A (en) * 2018-01-31 2018-08-10 天津大学 Servo-actuated real-time measurement apparatus and its method for rotating vane whole audience dynamic deformation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JUNSONG YU: "Short-time Energy Characterization of Low Velocity Impact Localization Method on CFRP Using FBG Sensors", 《THE 6TH ASIA PACIFIC OPTICAL SENSORS CONFERENCE 2016 ©》 *
余询: "一种高效率的桨叶载荷标定方法", 《直升机技术》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110595373A (en) * 2019-09-06 2019-12-20 沈阳天眼智云信息科技有限公司 Power transformer winding structure strain detection method based on optical fiber sensing
CN111289020A (en) * 2020-03-05 2020-06-16 山东大学 A Perception Method for Large Deformation Spatial Pose of Flexible Body
CN111289020B (en) * 2020-03-05 2021-07-13 山东大学 A Perception Method for Large Deformation Spatial Pose of Flexible Body
CN111677635A (en) * 2020-04-30 2020-09-18 中国电力科学研究院有限公司 A method and device for measuring the bending moment of wind turbine blades considering orthogonal effects
CN112432746A (en) * 2020-10-21 2021-03-02 中国人民解放军总参谋部第六十研究所 Method for determining swing load coupling degree of helicopter blade
CN112432746B (en) * 2020-10-21 2022-11-11 中国人民解放军总参谋部第六十研究所 Method for determining swing load coupling degree of helicopter blade
CN112525066A (en) * 2020-11-20 2021-03-19 中国直升机设计研究所 Helicopter blade surface strain gauge pasting and wire wiring method
CN113029405B (en) * 2021-02-04 2022-04-22 南京航空航天大学 Decoupling and calibration method of blade bending moment based on fiber-optic strain composite bridge
CN113029405A (en) * 2021-02-04 2021-06-25 南京航空航天大学 Blade bending moment decoupling and calibrating method based on optical fiber strain combined bridge circuit
CN114166141A (en) * 2021-11-23 2022-03-11 中国直升机设计研究所 A multi-circle marker detection method for helicopter high-speed rotating blade images
CN114166141B (en) * 2021-11-23 2023-05-26 中国直升机设计研究所 Multi-circular mark detection method for helicopter high-speed rotating blade image
CN115031647A (en) * 2022-04-12 2022-09-09 上海隧道工程有限公司 Deformation detection method for end plates
CN115165175A (en) * 2022-04-29 2022-10-11 南京航空航天大学 Bending moment decoupling method in rotor load identification based on fiber grating sensor
CN119223168A (en) * 2024-12-05 2024-12-31 南京信息职业技术学院 Determination method of helicopter blade flapping and shimmy decoupling position using fiber Bragg grating sensor
CN119223168B (en) * 2024-12-05 2025-07-11 南京信息职业技术学院 Fiber bragg grating sensor measuring method for decoupling position of flapping and shimmy of helicopter blade
CN119468956A (en) * 2025-01-14 2025-02-18 中国航空工业集团公司哈尔滨空气动力研究所 A method for measuring blade spacing in wind tunnel tests of coaxial rigid rotor models

Also Published As

Publication number Publication date
CN110006563B (en) 2020-09-18

Similar Documents

Publication Publication Date Title
CN110006563B (en) Distributed optical fiber decoupling measurement method for helicopter blade flapping and sway bending moment
CN108195554B (en) Six-component optical fiber aerodynamic force measurement balance and output signal combination method
CN103616157B (en) The quiet calibration system of wind-tunnel balance body axle system and method
CN110127078B (en) Optical fiber monitoring method for strain-deflection-bending moment state of helicopter blade structure
CN107271090B (en) A kind of aircraft wing moment of flexure method of real-time based on fiber grating
CN104613923B (en) A kind of deformation monitoring safety estimation system and appraisal procedure
CN105115438B (en) A kind of optical fiber sensing system temperature compensation
CN113029405B (en) Decoupling and calibration method of blade bending moment based on fiber-optic strain composite bridge
CN114754691B (en) Distributed optical fiber monitoring and inversion method for helicopter blade bending form
CN109556821B (en) Device for measuring relative position of balance center and wind tunnel test model
CN104111032A (en) Fiber bragg grating sensor network based large-scale structure body deformation measurement method
CN203116893U (en) Hidden weld joint stress monitoring device for orthotropic steel bridge decks
CN107505477B (en) A three-dimensional fiber Bragg grating wind speed and direction sensor and system
CN108317965A (en) A kind of measurement deformation structure and method with fiber grating
CN205373657U (en) Enamel amber strain sensor device based on sensitization structure temperature compensating
CN102252792B (en) Absolute axial force test method of bar
CN107367237A (en) A kind of deformation application of distributed optical fiber sensing system and Calibration Method
CN107044828A (en) A kind of optical fiber grating structure device for detecting deformation, network and method
CN114199156B (en) A Method of Airfoil Surface Shape Reconstruction Based on Dual FBG Shape Sensors
CN111649907A (en) Rib strip type micro-axial force annular two-component strain balance
CN106338384A (en) Wind generating set blade full-spanwise load measuring method
CN112763124B (en) Gear shaft torque monitoring device of lifting system and gear shaft
CN113970298A (en) FBG-based large-scale aerospace flexible structure overall deformation detection method
CN105277130A (en) Folded and unfolded structure fold area axial deformation optical fiber fast calculation method
CN210862556U (en) Bridge dynamic deflection monitoring system based on distributed optical fibers

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