CN110127078A - 直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法 - Google Patents
直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法 Download PDFInfo
- Publication number
- CN110127078A CN110127078A CN201910254398.4A CN201910254398A CN110127078A CN 110127078 A CN110127078 A CN 110127078A CN 201910254398 A CN201910254398 A CN 201910254398A CN 110127078 A CN110127078 A CN 110127078A
- Authority
- CN
- China
- Prior art keywords
- deflection
- strain
- segment
- equation
- amount
- 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
Links
- 238000005452 bending Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000012544 monitoring process Methods 0.000 title claims abstract description 22
- 239000013307 optical fiber Substances 0.000 claims abstract description 5
- 230000008569 process Effects 0.000 claims abstract description 5
- 230000008859 change Effects 0.000 claims description 14
- 239000000835 fiber Substances 0.000 claims description 9
- 238000006467 substitution reaction Methods 0.000 claims description 9
- 102100038032 F-box only protein 17 Human genes 0.000 claims description 6
- 102100022116 F-box only protein 2 Human genes 0.000 claims description 6
- 102100022115 F-box only protein 27 Human genes 0.000 claims description 6
- 102100026080 F-box only protein 44 Human genes 0.000 claims description 6
- 102100024513 F-box only protein 6 Human genes 0.000 claims description 6
- 101000878584 Homo sapiens F-box only protein 17 Proteins 0.000 claims description 6
- 101000824158 Homo sapiens F-box only protein 2 Proteins 0.000 claims description 6
- 101000824171 Homo sapiens F-box only protein 27 Proteins 0.000 claims description 6
- 101000913298 Homo sapiens F-box only protein 44 Proteins 0.000 claims description 6
- 101001052796 Homo sapiens F-box only protein 6 Proteins 0.000 claims description 6
- 230000003287 optical effect Effects 0.000 claims description 4
- 230000007935 neutral effect Effects 0.000 claims description 3
- 229920000049 Carbon (fiber) Polymers 0.000 claims 1
- 239000004917 carbon fiber Substances 0.000 claims 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims 1
- 238000000691 measurement method Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000036244 malformation Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 208000037656 Respiratory Sounds Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000005622 photoelectricity Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/165—Measuring 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
- G01D5/3537—Optical fibre sensor using a particular arrangement of the optical fibre itself
- G01D5/35374—Particular layout of the fiber
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
- G01M5/0016—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings of aircraft wings or blades
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
- G01M5/0041—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress
- G01M5/005—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress by means of external apparatus, e.g. test benches or portable test systems
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种直升机桨叶结构应变‑挠度‑弯矩状态的光纤监测方法,包括以下步骤:在直升机桨叶结构表面特征点位置布设光纤光栅传感器,求得应变分布曲线方程,进而可以推导该曲线任意点的应变值;将桨叶沿展向分成若干测量段。求得第一段上任意点的挠度;计算第一段的末尾点即第二段的起始点处的挠度值和转角大小;利用对第二段的应变分布曲线方程的二重积分,可求得第二段上的挠度方程;重复上述过程,可得到结构上每一段的挠度方程,即可得知结构上任意点的挠度值;根据挠曲线微分方程,得到弯矩与挠度关系,即可得知桨叶结构表面任意点的弯矩值。本发明的算法无需大量先验知识、简单快速、方便可靠、精确性强。
Description
技术领域
本发明属于结构健康监测的变形监测技术领域,具体提出了一种直升机桨叶结构应变-弯矩状态监测方法。
背景技术
叶片类机械是航空航天行业的关键设备,近年来,随着研究的深入,叶片类机械性能不断提高,不断向高速、重载、髙效、高可靠性方向发展,因此对机械系统的要求也越来越高。桨叶作为直升机的核心部件,是主要的承力对象,桨叶的损坏事故(裂纹、折断等),绝大部分是由于振动引起的,如直升机桨叶的挥舞运动、摆振运动等。因此对旋转桨叶状态进行监测是评估设计优劣、机械损伤动态监测、机械故障诊断、预测的重要方法。
目前,国内外针对结构变形状态监测方法的研究主要包括光电式测量系统和非光电式测量系统,其中光电式包括摄影测量法、摄像测量法、激光跟踪仪等;非光电式包括应变传感器、加速度传感器、位移传感器、光纤传感器等。光电式测量系统一般以影像捕捉和激光扫描跟踪为主,在目标结构给定位置上布置目标点或光源,测量系统以成像的方式感知目标点或光源,根据距离修正原理、反射原理和摄影测量原理,通过坐标转换得到目标点的空间位置,响应速度快,精度高。但该类方法只适用于外形比较平整、形式比较简单的结构,抗外界干扰能力比较差,对于大变形的柔性结构,光线传输易受影响,测量系统的可靠性很难保证。非光电式测量系统主要以接触式的测量方式为主,以应变信息或曲率信息为基础,适应较复杂的结构,可以直接将传感器嵌入到结构材料中,测量性能稳定。接触式的测量方式适应较复杂结构的形变重构,测量性能比较稳定。
因此,针对直升机桨叶工作环境恶劣,外界干扰较多的特点,选择使用非光电的接触式测量方法,并且选择相对先进且技术日趋成熟的光纤布拉格光栅传感器为测试元件,与传统传感器相比,光纤光栅传感器具有极高的灵敏度和精度,并且本征安全、抗电磁干扰强、绝缘强度高、体积小、频带宽。针对目前桨叶弯矩状态监测方法的不足,需要研究无需大量先验知识,能够适用于采样频率较低的常规光纤光栅解调仪,简单快速、方便可靠、实用性强的新方法。
发明内容
发明目的:为了克服现有技术中存在的不足,本发明提供一种直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法,该方法在结构待测点上布置光纤光栅传感器,测得有限点的应变值,对采集的应变值进行应变连续化,求得应变变化方程,根据结构的边界条件,利用对应变方程的二重积分求得第一段上任意点的挠度,根据第一段求得的挠度方程和转角变化方程,可得知第一段的末尾点即第二段的起始点出的挠度值和转角大小,由此作为已知条件,利用对第二段的应变方程的二重积分可求的第二段上的挠度方程,重复上述过程,即可得到结构上每一段的挠度方程,即可得知结构上任意点的挠度值,再根据挠曲线微分方程得到弯矩与挠度关系,即可得知结构上任意点的弯矩值。本发明的算法无需大量先验知识、简单快速、方便可靠、精确性强。
技术方案:为实现上述目的,本发明采用的技术方案为:
一种直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法,包括以下步骤:
步骤一、在直升机桨叶结构表面特征点位置布设光纤光栅传感器,获得桨叶表面相关离散测点的实时应变数据,对离散分布的光纤光栅传感器采集的应变值进行应变连续化,求得应变分布曲线方程,进而可以推导该曲线任意点的应变值;
在纯弯曲变形或以弯曲变形为主要变形方式下,等厚度结构表面应变ε(x)与弯曲挠度ω(x)之间的关系表示为:
式中,d表示微分,c为结构表面到结构中性面的垂直距离,其定义如下:
式中,εt为结构i节点上表面测得的应变值,εb为结构i节点下表面测得的应变值,h为结构厚度,在弯曲变形下,结构上下表面的应变值数值上近似相等,且有εb≈-εt,则c的表达式简化为:
c=h/2 (1.3)
通过光栅传感器测得桨叶结构上第i段起始点与末尾点处的应变值分别为εi、εi+1;结构在受弯情况时表面应变呈线性分布变化,通过线性插值得到第i段起始点与末尾点之间的应变变化方程,表示为:
其中,,为应变值,Δli为第i段长度,xi为第i段起始位置,x为长度值,
步骤二、根据直升机桨叶结构的边界条件,将桨叶沿展向分成n个测量段l1、l2、l3、…ln,利用对应变分布曲线方程的二重积分,求得第一段上任意点的挠度;
在已知桨叶结构上表面的第i段的起始点的转角和挠度条件下,通过一次积分可得知第i段转角θ(x)变化方程,如下所示:
对截面转角进行一次积分可得弯曲挠度,所以有:
在单边固支板结构弯曲条件下,起始点处有边界条件为:所以,对于桨叶结构的第一段挠度计算可由下式得到:
由于所以上式表示为:
由上式可求得第1段上任意点的挠度大小;
步骤三、根据第一段求得的挠度方程和转角变化方程,推导第一段的末尾点即第二段的起始点处的挠度值和转角大小:
当时,代入式(2.4)得第1段末尾点即第二段起始点挠度值
代入式(2.1)得第二段起始点转角
步骤四、根据步骤三计算所得第二段的起始点处的挠度值和转角数据,利用对第二段的应变分布曲线方程的二重积分,求得第二段上的挠度方程;
根据步骤三所得为边界条件,可知第二段转角变化方程为:
第二段挠度变化方程为:
步骤五、重复上述过程,可得到结构上每一段的挠度方程,即可得知桨叶结构上任意点的挠度值;
因为第i段的末尾点即为第i+1段的起始点,则第i+1段的边界条件可由第i段求得,重复步骤二和步骤三的一重积分可得转角方程,二重积分可得挠度方程。
步骤六、根据挠曲线微分方程得到弯矩与挠度关系,即可得知结构上任意点的弯矩值,纯弯曲情况下,弯矩和曲率间的关系式为:
其中,ρ为曲率半径,M是弯矩,E是弹性模量,I是截面惯性矩;
又有ds=ρdθ (6.2)
于是(5.1)式化为
挠度与截面转角关系式:
根据上式:
注意到代入式(5.3)得
式(6.6)即为旋转桨叶挠曲线的非线性微分方程;
由此可知挠度与弯矩的关系,再结合根据步骤三、步骤四、步骤五推导所得的应变与挠度关系,即可得出应变-弯矩关系,将式(2.4)代入式(6.6)即可得第1段弯矩:
再根据步骤三、步骤四、步骤五依次类推,即可得第2段至第i段的弯矩公式。
优选的:步骤一中采用双通道空分复用方法,在桨叶结构监测区域内布置平行于展长方向的四个光纤光栅传感器,分别记为FBG1、FBG2、FBG3、FBG4,同时在靠近桨叶迎风边的位置分别布置平行展长方向的四个光纤布拉格光栅传感器,分别记为FBG5、FBG6、FBG7、FBG8;将这些光纤光栅传感器粘贴于桨叶结构表面,采用光纤跳线将FBG1、FBG2、FBG3、FBG4和FBG5、FBG6、FBG7、FBG8分别串行连接,以此构成分布式传感器网络。
优选的:步骤一中在直升机桨叶结构表面布置分布式光纤光栅传感器,桨叶一端固定,一端处于自由状态,光纤光栅传感器在靠近固定端区域的布局数量应相对较多,且间距较小,而在桨叶自由端方向,光纤光栅传感器数量适当减小,间距相应变大。
本发明相比现有技术,具有以下有益效果:
本发明采用分布式光纤布拉格光栅传感器网络获取结构中不同位置的应变信息,先使用基于应变—挠度转换计算得出结构变形挠度,再使用挠曲线微分方程将挠度转化为弯矩,该算法避免了应变弯矩直接转换关系中的近似方程导致的误差,适用于大展弦比叶片类结构,无需大量先验知识,简单快速、方便可靠、实用性强。
附图说明
图1分布式光纤布拉格光栅传感器布局图;
图2桨叶结构分段示意图;
图3弯矩状态监测算法流程图;
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
一种直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法,如图3所示,包括以下步骤:
步骤一、如图1所示,在直升机桨叶结构表面特征点位置布设光纤光栅传感器,获得桨叶表面相关离散测点的实时应变数据。对离散分布的光纤光栅传感器采集的应变值进行应变连续化,求得应变分布曲线方程,进而可以推导该曲线任意点的应变值。
在直升机桨叶结构表面布置分布式光纤光栅传感器,桨叶一端固定,一端处于自由状态。桨叶运动具有多种低阶振型,振型复杂。考虑四阶以下的振型,振阶越高,其变形越小。综合各阶振型,旋转桨叶越远离固定一端,其内部应力越小,应变也越小。因此,光纤光栅传感器在靠近固定端区域的布局数量应相对较多,且间距较小。而在桨叶自由端方向,光纤光栅传感器数量可以适当减小,间距相应变大。
采用双通道空分复用方法,在桨叶结构监测区域内布置平行于展长方向的四个光纤光栅传感器FBG1、FBG2、FBG3、FBG4,同时在靠近桨叶迎风边的位置分别布置平行展长方向的四个光纤布拉格光栅传感器FBG5、FBG6、FBG7、FBG8。将这些光纤光栅传感器粘贴于桨叶结构表面,采用光纤跳线将FBG1、FBG2、FBG3、FBG4和FBG5、FBG6、FBG7、FBG8分别串行连接,以此构成分布式传感器网络。
在纯弯曲变形或以弯曲变形为主要变形方式下,等厚度结构表面应变ε(x)与弯曲挠度ω(x)之间的关系可以表示为:
式中,c为结构表面到结构中性面的垂直距离,其定义如下:
式中,εt为结构i节点上表面测得的应变值,εb为结构i节点下表面测得的应变值,h为结构厚度。通常,在弯曲变形下,结构上下表面的应变值数值上近似相等,且有εb≈-εt。则c的表达式可简化为:
c=h/2 (1.3)
通过光栅传感器测得桨叶结构上第i段起始点与末尾点处的应变值分别为εi、εi+1。假设结构在受弯情况时表面应变呈线性分布变化,所以可通过线性插值得到第i段起始点与末尾点之间的应变变化方程,可以表示为:
步骤二、如图2所示,根据直升机桨叶结构的边界条件,将桨叶沿展向分成若干测量段l1、l2、l3、…ln。利用对应变分布曲线方程的二重积分,可以求得第一段上任意点的挠度;
在已知桨叶结构上表面的第i段的起始点的转角和挠度条件下,通过一次积分可得知第i段转角θ(x)变化方程,如下所示:
对截面转角进行一次积分可得弯曲挠度,所以有:
显然,为完成上式积分,必须要知道起始点的转角和挠度才能再次利用分段拟合方法完成上述计算。在单边固支板结构弯曲条件下,起始点处有边界条件为:所以,对于桨叶结构的第一段挠度计算可由下式得到:
由于所以上式可以表示为:
由上式可求得第1段上任意点的挠度大小。
步骤三、根据第一段求得的挠度方程和转角变化方程,可以推导第一段的末尾点即第二段的起始点处的挠度值和转角大小。
当时,代入式(2.4)得第1段末尾点即第二段起始点挠度值
代入式(2.1)得第二段起始点转角
步骤四、根据步骤三计算所得第二段的起始点处的挠度值和转角数据,利用对第二段的应变分布曲线方程的二重积分,可求得第二段上的挠度方程;
根据步骤三所得为边界条件,可知第二段转角变化方程为:
第二段挠度变化方程为:
步骤五、重复上述过程,可得到结构上每一段的挠度方程,即可得知桨叶结构上任意点的挠度值。
因为第i段的末尾点即为第i+1段的起始点,则第i+1段的边界条件可由第i段求得,重复步骤二和步骤三的一重积分可得转角方程,二重积分可得挠度方程。
步骤六、根据挠曲线微分方程以,得到弯矩与挠度关系,即可得知结构上任意点的弯矩值。纯弯曲情况下,弯矩和曲率间的关系式为:
其中,ρ为曲率半径,M是弯矩,E是弹性模量,I是截面惯性矩。横力弯曲时,梁截面上有弯矩也有剪力,此式只代表弯矩对弯矩变形的影响。对跨度远大于截面高度的梁,剪力对弯曲变形的影响可以忽略。
又有ds=ρdθ (6.2)
于是(5.1)式化为
挠度与截面转角关系式:
根据上式:
注意到代入式(5.3)得
式(6.6)即为旋转桨叶挠曲线的非线性微分方程,此方程适用于桨叶弯曲变形的绝大多数情况。
由此可知挠度与弯矩的关系,再结合根据步骤三、步骤四、步骤五推导所得的应变与挠度关系,即可得出应变-弯矩关系。将式(2.4)代入式(6.6)即可得第1段弯矩:
再根据步骤三、步骤四、步骤五依次类推,即可得第2段至第i段的弯矩公式。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法,其特征在于,包括以下步骤:
步骤一、在直升机桨叶结构表面特征点位置布设光纤光栅传感器,获得桨叶表面相关离散测点的实时应变数据,对离散分布的光纤光栅传感器采集的应变值进行应变连续化,求得应变分布曲线方程,进而可以推导该曲线任意点的应变值;
在纯弯曲变形或以弯曲变形为主要变形方式下,等厚度结构表面应变ε(x)与弯曲挠度ω(x)之间的关系表示为:
式中,c为结构表面到结构中性面的垂直距离,其定义如下:
式中,εt为结构i节点上表面测得的应变值,εb为结构i节点下表面测得的应变值,h为结构厚度,在弯曲变形下,结构上下表面的应变值数值上近似相等,且有εb≈-εt,则c的表达式简化为:
c=h/2 (1.3)
通过光栅传感器测得桨叶结构上第i段起始点与末尾点处的应变值分别为εi、εi+1;结构在受弯情况时表面应变呈线性分布变化,通过线性插值得到第i段起始点与末尾点之间的应变变化方程,表示为:
其中,为应变值,Δli为第i段长度,xi为第i段起始位置,x为长度值,
步骤二、根据直升机桨叶结构的边界条件,将桨叶沿展向分成n个测量段l1、l2、l3、…ln,利用对应变分布曲线方程的二重积分,求得第一段上任意点的挠度;
在已知桨叶结构上表面的第i段的起始点的转角和挠度条件下,通过一次积分可得知第i段转角θ(x)变化方程,如下所示:
对截面转角进行一次积分可得弯曲挠度,所以有:
在单边固支板结构弯曲条件下,起始点处有边界条件为:所以,对于桨叶结构的第一段挠度计算可由下式得到:
由于所以上式表示为:
由上式可求得第1段上任意点的挠度大小;
步骤三、根据第一段求得的挠度方程和转角变化方程,推导第一段的末尾点即第二段的起始点处的挠度值和转角大小:
当时,代入式(2.4)得第1段末尾点即第二段起始点挠度值
代入式(2.1)得第二段起始点转角
步骤四、根据步骤三计算所得第二段的起始点处的挠度值和转角数据,利用对第二段的应变分布曲线方程的二重积分,求得第二段上的挠度方程;
根据步骤三所得为边界条件,可知第二段转角变化方程为:
第二段挠度变化方程为:
步骤五、重复上述过程,可得到结构上每一段的挠度方程,即可得知桨叶结构上任意点的挠度值;
因为第i段的末尾点即为第i+1段的起始点,则第i+1段的边界条件可由第i段求得,重复步骤二和步骤三的一重积分可得转角方程,二重积分可得挠度方程。
步骤六、根据挠曲线微分方程得到弯矩与挠度关系,即可得知结构上任意点的弯矩值,纯弯曲情况下,弯矩和曲率间的关系式为:
其中,ρ为曲率半径,M是弯矩,E是弹性模量,I是截面惯性矩;
又有:ds=ρdθ (6.2)
于是(5.1)式化为:
挠度与截面转角关系式:
根据上式:
注意到代入式(5.3)得
式(6.6)即为旋转桨叶挠曲线的非线性微分方程;
由此可知挠度与弯矩的关系,再结合根据步骤三、步骤四、步骤五推导所得的应变与挠度关系,即可得出应变-弯矩关系,将式(2.4)代入式(6.6)即可得第1段弯矩:
再根据步骤三、步骤四、步骤五依次类推,即可得第2段至第i段的弯矩公式。
2.根据权利要求1所述直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法,其特征在于:步骤一中采用双通道空分复用方法,在桨叶结构监测区域内布置平行于展长方向的四个光纤光栅传感器,分别记为FBG1、FBG2、FBG3、FBG4,同时在靠近桨叶迎风边的位置分别布置平行展长方向的四个光纤布拉格光栅传感器,分别记为FBG5、FBG6、FBG7、FBG8;将这些光纤光栅传感器粘贴于桨叶结构表面,采用光纤跳线将FBG1、FBG2、FBG3、FBG4和FBG5、FBG6、FBG7、FBG8分别串行连接,以此构成分布式传感器网络。
3.根据权利要求3所述直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法,其特征在于:步骤一中在直升机桨叶结构表面布置分布式光纤光栅传感器,桨叶一端固定,一端处于自由状态,光纤光栅传感器在靠近固定端区域的布局数量应相对较多,且间距较小,而在桨叶自由端方向,光纤光栅传感器数量适当减小,间距相应变大。
4.根据权利要求3所述直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法,其特征在于:传感网络所使用的光纤光栅解调仪采样频率为2KHz。
5.根据权利要求4所述直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法,其特征在于:桨叶采用单边固支碳纤维结构。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910254398.4A CN110127078B (zh) | 2019-03-31 | 2019-03-31 | 直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910254398.4A CN110127078B (zh) | 2019-03-31 | 2019-03-31 | 直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110127078A true CN110127078A (zh) | 2019-08-16 |
CN110127078B CN110127078B (zh) | 2022-04-22 |
Family
ID=67568799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910254398.4A Active CN110127078B (zh) | 2019-03-31 | 2019-03-31 | 直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110127078B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111750790A (zh) * | 2020-04-22 | 2020-10-09 | 江苏法尔胜光电科技有限公司 | 一种基于分布式光纤的管道二维变形监测方法 |
CN112001098A (zh) * | 2020-06-30 | 2020-11-27 | 同恩(上海)工程技术有限公司 | 一种基于插值反演方式的构件截面内力监测方法及系统 |
CN112287454A (zh) * | 2020-10-28 | 2021-01-29 | 北京天骥空间科技有限公司 | 一种基于应力应变的箭体外形重构技术 |
CN112304756A (zh) * | 2020-10-16 | 2021-02-02 | 中国航发四川燃气涡轮研究院 | 一种纤维增强复合材料环状结构周向拉伸性能表征方法 |
CN112729147A (zh) * | 2020-12-28 | 2021-04-30 | 南京航空航天大学 | 一种复杂环境下桨叶应变测量试验系统及试验方法 |
CN112964459A (zh) * | 2021-02-05 | 2021-06-15 | 中国航发沈阳发动机研究所 | 一种考虑叶片扭转的光纤限制值计算方法 |
CN113221354A (zh) * | 2021-05-11 | 2021-08-06 | 西安电子科技大学 | 一种管道弯曲变形拟合算法 |
CN113945164A (zh) * | 2021-10-14 | 2022-01-18 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种直升机桨叶自动测量装置 |
CN116388859A (zh) * | 2022-10-31 | 2023-07-04 | 广东电网有限责任公司中山供电局 | 一种光纤状态监测数据采集装置、方法、设备和介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020025744A1 (en) * | 1998-05-14 | 2002-02-28 | Mccarthy Peter T. | Methods for creating large scale focused blade deflections |
WO2013091638A1 (en) * | 2011-12-20 | 2013-06-27 | Vestas Wind Systems A/S | Control method for a wind turbine, and wind turbine |
CN105043285A (zh) * | 2015-09-10 | 2015-11-11 | 北京航空航天大学 | 一种基于光纤传感器的测量变截面梁变形的方法 |
CN107271090A (zh) * | 2017-06-26 | 2017-10-20 | 北京航空航天大学 | 一种基于光纤光栅的飞机机翼弯矩实时监测方法 |
-
2019
- 2019-03-31 CN CN201910254398.4A patent/CN110127078B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020025744A1 (en) * | 1998-05-14 | 2002-02-28 | Mccarthy Peter T. | Methods for creating large scale focused blade deflections |
WO2013091638A1 (en) * | 2011-12-20 | 2013-06-27 | Vestas Wind Systems A/S | Control method for a wind turbine, and wind turbine |
CN105043285A (zh) * | 2015-09-10 | 2015-11-11 | 北京航空航天大学 | 一种基于光纤传感器的测量变截面梁变形的方法 |
CN107271090A (zh) * | 2017-06-26 | 2017-10-20 | 北京航空航天大学 | 一种基于光纤光栅的飞机机翼弯矩实时监测方法 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111750790A (zh) * | 2020-04-22 | 2020-10-09 | 江苏法尔胜光电科技有限公司 | 一种基于分布式光纤的管道二维变形监测方法 |
CN112001098A (zh) * | 2020-06-30 | 2020-11-27 | 同恩(上海)工程技术有限公司 | 一种基于插值反演方式的构件截面内力监测方法及系统 |
CN112001098B (zh) * | 2020-06-30 | 2023-11-17 | 同恩(上海)工程技术有限公司 | 一种基于插值反演方式的构件截面内力监测方法及系统 |
CN112304756A (zh) * | 2020-10-16 | 2021-02-02 | 中国航发四川燃气涡轮研究院 | 一种纤维增强复合材料环状结构周向拉伸性能表征方法 |
CN112287454A (zh) * | 2020-10-28 | 2021-01-29 | 北京天骥空间科技有限公司 | 一种基于应力应变的箭体外形重构技术 |
CN112729147A (zh) * | 2020-12-28 | 2021-04-30 | 南京航空航天大学 | 一种复杂环境下桨叶应变测量试验系统及试验方法 |
CN112964459A (zh) * | 2021-02-05 | 2021-06-15 | 中国航发沈阳发动机研究所 | 一种考虑叶片扭转的光纤限制值计算方法 |
CN113221354A (zh) * | 2021-05-11 | 2021-08-06 | 西安电子科技大学 | 一种管道弯曲变形拟合算法 |
CN113945164A (zh) * | 2021-10-14 | 2022-01-18 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种直升机桨叶自动测量装置 |
CN116388859A (zh) * | 2022-10-31 | 2023-07-04 | 广东电网有限责任公司中山供电局 | 一种光纤状态监测数据采集装置、方法、设备和介质 |
CN116388859B (zh) * | 2022-10-31 | 2023-08-18 | 广东电网有限责任公司中山供电局 | 一种光纤状态监测数据采集装置、方法、设备和介质 |
Also Published As
Publication number | Publication date |
---|---|
CN110127078B (zh) | 2022-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110127078A (zh) | 直升机桨叶结构应变-挠度-弯矩状态的光纤监测方法 | |
CN110006563A (zh) | 直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法 | |
CN107271090B (zh) | 一种基于光纤光栅的飞机机翼弯矩实时监测方法 | |
CN101620000B (zh) | 输电线路覆冰量的计算方法 | |
CN109766617B (zh) | 一种基于应变传感器的位移场重构方法 | |
CN201731984U (zh) | 基于无线倾角的桥梁结构挠度测量系统 | |
CN101852597B (zh) | 温度不敏感型光纤光栅倾斜传感器 | |
CN103776473A (zh) | 一种基于光纤光栅传感器的飞机燃油油量测量系统 | |
CN203688576U (zh) | 一种光纤Bragg光栅动态风速风向仪 | |
CN107765030A (zh) | 一种感知流速变化的压电纤维阵列传感器及其测量方法 | |
CN110133324A (zh) | 一种差动式光纤光栅加速度传感装置 | |
CN105841857A (zh) | 一种并联式五维力传感器 | |
CN105547554B (zh) | 一种用于测量小型螺旋翼升力和扭矩的装置及方法 | |
CN106525301A (zh) | 一种基于分布式光纤感测的力和位移测量方法及传感器 | |
CN105277130A (zh) | 折叠展开结构褶皱区域轴向变形光纤快速计算方法 | |
CN112649052A (zh) | 基于全光纤光栅传感网络的船体多参量监测系统及方法 | |
CN111664783A (zh) | 一种大变形位移传感器和测量方法 | |
CN209745726U (zh) | 一种微型室内静力触探测试系统 | |
CN201688838U (zh) | 一种温度不敏感型光纤光栅倾斜传感器 | |
CN211717557U (zh) | 一种基于光纤光栅压力传感阵列的液位检测系统 | |
CN107796421A (zh) | 快速确定fbg反射光谱异常模式的方法 | |
CN110487529B (zh) | 一种利用角速度传感器测量大展弦比翼面弯矩的方法 | |
CN105352445A (zh) | 单边固支板结构变形光纤模式辨识系统及标定方法与应用 | |
CN101995226B (zh) | 新型光纤光栅的多点裂缝监测器 | |
ITFI20080223A1 (it) | Dispositivo e procedimento per la rilevazione o la generazione di una curvatura |
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 |