CN111223173A - 基于光纤瑞利散射的柔性面板形状构造方法 - Google Patents
基于光纤瑞利散射的柔性面板形状构造方法 Download PDFInfo
- Publication number
- CN111223173A CN111223173A CN202010193429.2A CN202010193429A CN111223173A CN 111223173 A CN111223173 A CN 111223173A CN 202010193429 A CN202010193429 A CN 202010193429A CN 111223173 A CN111223173 A CN 111223173A
- Authority
- CN
- China
- Prior art keywords
- optical fiber
- axis
- optical
- flexible panel
- point
- 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.)
- Pending
Links
- 239000013307 optical fiber Substances 0.000 title claims abstract description 247
- 238000010276 construction Methods 0.000 title claims abstract description 10
- 230000003287 optical effect Effects 0.000 claims abstract description 77
- 230000005540 biological transmission Effects 0.000 claims abstract description 18
- 238000012545 processing Methods 0.000 claims description 25
- 238000000034 method Methods 0.000 claims description 24
- 239000000835 fiber Substances 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000002168 optical frequency-domain reflectometry Methods 0.000 claims 3
- 238000005259 measurement Methods 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 3
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 206010000369 Accident Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011284 combination treatment Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
-
- 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/002—Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Optical Transform (AREA)
Abstract
本发明提供一种基于光纤瑞利散射的柔性面板形状构造方法,本发明针对平铺的柔性面板建立三维坐标系,使X‑Y坐标系位于平铺柔性面板上,在平铺柔性面板的下表面铺设多条平行于Y轴且相连的第一光纤,上表面铺设多条平行于X轴且相连的第二光纤,使第一光纤与第二光纤垂直正交,下表面的第一光纤与上表面的第二光纤相连,形成分布式光纤传感网络,对光信号在分布式光纤传感网络的传输中产生的瑞利散射信号进行处理,获得各条第一光纤和第二光纤上各个位置处的曲率,根据各条第一光纤和第二光纤上各个位置处的曲率,对应构建出该柔性面板的下表面曲面和上表面曲面,进而构建出该柔性面板的三维曲面,由此可以提高柔性面板形变测量准确度。
Description
技术领域
本发明属于分布式光纤传感领域,具体涉及一种基于光纤瑞利散射的柔性面板形状构造方法。
背景技术
曲面板状结构广泛应用于空间站,人造卫星,航空航天等重要领域之中,如构成太阳能帆板,航空航天机翼。太空空间由于无阻尼,宇宙粒子流冲击,各种宇宙微粒扰动,长期如此,极易对曲面板状结构造成损坏,导致系统的损坏。除此之外,随着航天航空事业的迅猛发展,由此引出的安全问题也越来越多,飞机失事案例屡见不鲜。飞机在飞行过程中,机翼极易受到气流的影响,导致机翼结构不稳定,长时间的影响将导致其机翼损坏,导致飞行事故。因此,对曲面板状结构的实时结构监测具有重要的意义。
目前所使用的曲面重构传统方法多是使用光纤光栅、应变片等离散方式,无法在曲面上获得多点的传感信息,因此空间分辨率不足,重构出的形状误差大,无法满足现代航空航天的要求。
发明内容
本发明提供一种基于光纤瑞利散射的柔性面板形状构造方法,以解决目前采用分布式传感光纤进行形变测量的方法的空间分辨率和测量准确度较低的问题。
根据本发明实施例的第一方面,提供一种基于光纤瑞利散射的柔性面板形状构造方法,包括:
步骤S110、针对待测的柔性面板,选定该柔性面板的上表面的一点作为零点,将平铺的柔性面板上表面上通过该零点的一条线作为X轴,将平铺的柔性面板上表面上通过该零点并与该X轴垂直的线作为Y轴,将通过该零点并垂直于该平铺的柔性面板的线作为Z轴,从而建立三维坐标系;
步骤S120、在平铺的柔性面板的下表面和上表面对应铺设多条平行于Y轴的第一光纤和多条平行于X轴的第二光纤,各条相邻第一光纤之间以及各条相邻第二光纤之间相连,该下表面的一条第一光纤的自由端作为下表面光信号出口,与该上表面的一条第二光纤的自由端相连,从而形成分布式光纤传感网络,每条第一光纤和第二光纤上的各个位置都对应于唯一的X-Y轴坐标,且针对每条第一光纤,其都与对应第二光纤垂直正交,该下表面的另一条第一光纤的自由端作为下表面光信号入口,与光频域反射仪OFDR(OpticalFrequency Domain Reflectometer)系统相连;
步骤S130、在测量该柔性面板的形变时,所述OFDR系统向该下表面光信号入口输入光信号,以使所述光信号在该分布式光纤传感网络的传输过程中产生瑞利散射信号,该瑞利散射信号向后传输给所述OFDR系统,所述OFDR系统与处理装置相连,并将该瑞利散射信号传输给该处理装置,该处理装置对该瑞利散射信号进行处理,获得各条第一光纤和第二光纤上各个位置处的曲率;
步骤S140、该处理装置根据各条第一光纤上各个位置处的曲率,构建出对应条数的Y-Z轴二维曲线,根据各条Y-Z轴二维曲线构建出该柔性面板的下表面曲面,并且根据各条第二光纤上各个位置处的曲率,构建出对应条数的X-Z轴二维曲线,根据各条X-Z轴二维曲线构建出该柔性面板的上表面曲面;
步骤S150、所述处理装置对该柔性面板的下表面曲面和上表面曲面进行结合处理,获得该柔性面板的三维曲面。
在一种可选的实现方式中,所述步骤S140中根据各条Y-Z轴二维曲线构建出该柔性面板的下表面曲面包括:采用二维插值算法对该Y-Z轴二维曲线进行二维插值,获得该下表面曲面;
所述步骤S140中根据各条X-Z轴二维曲线构建出该柔性面板的上表面曲面包括:采用二维插值算法对该X-Z轴二维曲线进行二维插值,获得该上表面曲面。
在另一种可选的实现方式中,所述步骤S120中该上表面上与该下表面光信号出口相连的一条第二光纤的自由端为上表面光信号入口,选取所述零点为该下表面光信号入口,X轴上的一点(x0,0)为该下表面光信号出口,该零点为该上表面光信号入口,x0为非零值。
在另一种可选的实现方式中,所述步骤S130中该处理装置对该瑞利散射信号进行处理,获得各条第一光纤和第二光纤上各个位置处的曲率包括:
步骤S131、针对每条第一光纤,若光信号在该第一光纤中的传输方向与Y轴正方向相反,则将该第一光纤的光信号输入端作为该第一光纤的末端,光信号输出端作为该第一光纤的始端,若光信号在该第一光纤中的传输方向与Y轴正方向相同,则将该第一光纤的光信号输入端作为该第一光纤的始端,光信号输出端作为该第一光纤的末端;针对每条第二光纤,若光信号在该第二光纤中的传输方向与X轴正方向相反,则将该第二光纤的光信号输入端作为该第二光纤的末端,光信号输出端作为该第二光纤的始端,若光信号在该第二光纤中的传输方向与X轴正方向相同,则将该第二光纤的光信号输入端作为该第二光纤的始端,光信号输出端作为该第二光纤的末端;
步骤S132、该处理装置对该瑞利散射信号进行处理,获得各条第一光纤和第二光纤上各个位置处的波长漂移值;
步骤S133、从始端到末端,对各条第一光纤和第二光纤上各个位置处的波长漂移值进行重排;
步骤S134、根据各条第一光纤和第二光纤上重排后的波长漂移值,计算出各条第一光纤和第二光纤上各个位置处的曲率。
在另一种可选的实现方式中,所述步骤S132包括:该处理装置对该瑞利散射信号进行频域转换和互相关运算,获得各条第一光纤和第二光纤上各个位置处的波长漂移值;所述步骤S134包括:根据以下公式计算出各条第一光纤和第二光纤上各个位置处的曲率c:c=Δλ/κ=Δλ/(1-pe)λch
其中Δλ表示第一光纤或第二光纤对应位置处的波长漂移值,κ表示应变系数,pe表示光纤的弹光系数,λc表示中心波长,h表示柔性面板的厚度。
在另一种可选的实现方式中,所述步骤S140包括根据以下向量参数方程构建下表面各条第一光纤变形后的二维曲线:r(s)=y(s)j+z(s)k,其中y(s)表示下表面二维曲线上各点的Y轴坐标函数,z(s)表示下表面二维曲线上各点的Z轴坐标函数,s表示下表面各条第一光纤长度的参量,0≤s≤Ly,Ly是下表面各条第一光纤的总长度;
该下表面二维曲线上任意位置对应的点表示为:
式中,yn+1表示该点的Y轴坐标值,yn表示沿Y轴正方向,其上一点的Y轴坐标值,zn+1表示该点的Z轴坐标值,zn表示该上一点的Z轴坐标值,ds表示二维曲线的微小弧度;
式中,cy表示该点的曲率;
根据以下向量参数方程构建上表面各条第二光纤变形后的曲面:r(s)=x(s)i+z(s)k,其中x(s)表示上表面二维曲线上各点的X轴坐标函数,z(s)表示上表面二维曲线上各点的Z轴坐标函数,s表示上表面各条第二光纤长度的参量,0≤s≤Lx,Lx是上表面各条第二光纤的总长度;
该上表面二维曲线上任意位置对应的点表示为:
式中,xn+1表示该点的X轴坐标值,xn表示沿X轴正方向,其上一点的X轴坐标值,zn+1表示该点的Z轴坐标值,zn表示该上一点的Z轴坐标值,ds表示二维曲线的微小弧度;
式中,cx表示该点的曲率。
在另一种可选的实现方式中,所述ds等于两相邻点之间的连线长度。
在另一种可选的实现方式中,所述步骤S150包括:从下表面曲面和上表面曲面中查找到对应的点P(x,y,z1)和Q(x,y,z2),使z=(z1+z2),将点(x,y,z)作为该柔性面板的三维曲面上的对应点。
在另一种可选的实现方式中,所述柔性面板在正常使用时处于弯曲状态。
本发明的有益效果是:
1、相较于传统方式,光纤具有抗电磁干扰,抗高温等优点,具有极高的实用性,在重构柔性曲面时,本发明将光纤上的各个点都可以作为传感点,重构出来的图形误差小;本发明针对平铺的柔性面板建立三维坐标系,使X-Y坐标系位于平铺柔性面板上,在平铺柔性面板的下表面铺设多条平行于Y轴的第一光纤,上表面铺设多条平行于X轴的第二光纤,使第一光纤与第二光纤垂直正交,并且各个相邻第一光纤之间、各个相邻第二光纤之间以及具有自由端的对应第一光纤与第二光纤之间相连,构成分布式光纤传感网络,光信号在分布式光纤传感网络的传输过程中产生向后传输的瑞利散射信号,根据各条第一光纤对应的瑞利散射信号构建出对应条数的Y-Z轴二维曲线,从而构建出该柔性面板的下表面曲面,根据各条第二光纤对应的瑞利散射信号构建出对应条数的X-Z轴二维曲线,从而构建出该柔性面板的上表面曲面,对该柔性面板的下表面曲面和上表面曲面进行结合处理,获得该柔性面板的三维曲面,由此可以减小应力形变带来的微小误差,提高柔性面板形变测量准确度,本发明采用光纤对柔性曲面测量的方式,将柔性曲面离散成多条光纤曲线进行重构,空间分辨率极高,本发明重构出来的整个图形,空间分辨率极高,准确率极高,适用于航空航天等高精度领域;
2、本发明通过采用二维插值算法对该Y-Z轴二维曲线进行二维插值,获得该下表面曲面,采用二维插值算法对该X-Z轴二维曲线进行二维插值,获得该上表面曲面,采用二维插值的方式,在柔性曲面相邻光纤之间插入多条曲线,提高相邻两条光纤的空间分辨率,可以提高柔性面板的上下表面曲面构建准确度,提高形变测量的准确度和空间分辨率。
附图说明
图1是本发明基于光纤瑞利散射的柔性面板形状构造方法的一个实施例流程图;
图2是本发明分布式光纤传感网络的一个实施例结构示意图;
图3是本发明单条曲线重构算法的示意图;
图4是采用本发明方法构造出的三维曲面效果图。
具体实施方式
为了使本技术领域的人员更好地理解本发明实施例中的技术方案,并使本发明实施例的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明实施例中技术方案作进一步详细的说明。
在本发明的描述中,除非另有规定和限定,需要说明的是,术语“连接”应做广义理解,例如,可以是机械连接或电连接,也可以是两个元件内部的连通,可以是直接相连,也可以通过中间媒介间接相连,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
参见图1,为本发明基于光纤瑞利散射的柔性面板形状构造方法的一个实施例流程图。该方法可以包括以下步骤:
步骤S110、针对待测的柔性面板,选定该柔性面板的上表面上的一点作为零点,将平铺的柔性面板的上表面上通过该零点的一条线作为X轴,将平铺的柔性面板上通过该零点并与该X轴垂直的线作为Y轴,将通过该零点并垂直于该平铺的柔性面板的线作为Z轴,从而建立三维坐标系。
步骤S120、结合图2所示,在平铺的柔性面板4的下表面铺设多条平行于Y轴的第一光纤1(如图2中虚线所示),各条相邻第一光纤1之间通过过渡光纤3连接;在平铺的柔性面板4的上表面铺设多条平行于X轴的第二光纤2(如图2中实线所示),各条相邻第二光纤2之间通过过渡光纤3连接,该下表面的一条第一光纤1-2的自由端作为下表面光信号出口(x0,0),通过光纤连接段5与该上表面的一条第二光纤2-1的自由端相连,x0为非零值,从而形成分布式光纤传感网络,每条第一光纤1和第二光纤2上的各个位置都对应于唯一的X-Y轴坐标,且针对每条第一光纤1,其都与对应第二光纤2垂直正交,该下表面的另一条第一光纤1-1的自由端作为下表面光信号入口,与光频域反射仪OFDR(Optical Frequency DomainReflectometer)系统相连。
步骤S130、在测量该柔性面板的形变时,所述OFDR系统向该下表面光信号入口提供光信号,以使所述光信号在该分布式光纤传感网络的传输过程中产生瑞利散射信号,该瑞利散射信号向后传输给所述OFDR系统,所述OFDR系统与处理装置相连,并将该瑞利散射信号传输给该处理装置,该处理装置对该瑞利散射信号进行处理,获得各条第一光纤和第二光纤上各个位置处的曲率。
步骤S140、该处理装置根据各条第一光纤上各个位置处的曲率,构建出对应条数的Y-Z轴二维曲线,根据各条Y-Z轴二维曲线构建出该柔性面板的下表面曲面,并且根据各条第二光纤上各个位置处的曲率,构建出对应条数的X-Z轴二维曲线,根据各条X-Z轴二维曲线构建出该柔性面板的上表面曲面。
步骤S150、所述处理装置对该柔性面板的下表面曲面和上表面曲面进行结合处理,获得该柔性面板的三维曲面。
本实施例中,由于各条第一光纤和第二光纤上各个位置处的X-Y轴坐标是唯一的,而且各个位置距离该下表面光信号入口的距离也是固定不变的,因此处理装置根据接收到的瑞利散射信号的时间,可以确定该瑞利散射信号是在哪个位置产生的,从而根据该瑞利散射信号,获得各个位置处的曲率。
由上述实施例可见,相较于传统方式,光纤具有抗电磁干扰,抗高温等优点,具有极高的实用性,在重构柔性曲面时,本发明将光纤上的各个点都可以作为传感点,重构出来的图形误差小;本发明针对平铺的柔性面板建立三维坐标系,使X-Y坐标系位于平铺柔性面板上,在平铺柔性面板的下表面铺设多条平行于Y轴的第一光纤,上表面铺设多条平行于X轴的第二光纤,使第一光纤与第二光纤垂直正交,并且各个相邻第一光纤之间、各个相邻第二光纤之间以及具有自由端的对应第一光纤与第二光纤之间相连,构成分布式光纤传感网络,光信号在分布式光纤传感网络的传输过程中产生向后传输的瑞利散射信号,根据各条第一光纤对应的瑞利散射信号构建出对应条数的Y-Z轴二维曲线,从而构建出该柔性面板的下表面曲面,根据各条第二光纤对应的瑞利散射信号构建出对应条数的X-Z轴二维曲线,从而构建出该柔性面板的上表面曲面,对该柔性面板的下表面曲面和上表面曲面进行结合处理,获得该柔性面板的三维曲面,由此可以减小应力形变带来的微小误差,提高柔性面板形变测量准确度,本发明采用光纤对柔性曲面测量的方式,将柔性曲面离散成多条光纤曲线进行重构,空间分辨率极高,本发明重构出来的整个图形,空间分辨率极高,准确率极高,适用于航空航天等高精度领域。
另外,结合图2所示,所述步骤S120中该上表面上与该下表面光信号出口相连的一条第二光纤的自由端为上表面光信号入口,选取所述零点为该下表面光信号入口,X轴上的一点(x0,0)为该下表面光信号出口,该零点为该上表面光信号入口,x0为非零值。本发明通过对上表面光信号入口、下表面光信号入口和下表面光信号出口的坐标进行选择,可以使上表面和下表面的坐标相对应,避免后续重构形状时需要对波长漂移值做进一步重排,从而在一定程度上可以提高重构效率。
所述步骤S130中所述处理装置对该瑞利散射信号进行处理,获得各条第一光纤和第二光纤上各个位置处的曲率包括:
步骤S131、针对每条第一光纤,若光信号在该第一光纤中的传输方向与Y轴正方向相反,则将该第一光纤的光信号输入端作为该第一光纤的末端,光信号输出端作为该第一光纤的始端,若光信号在该第一光纤中的传输方向与Y轴正方向相同,则将该第一光纤的光信号输入端作为该第一光纤的始端,光信号输出端作为该第一光纤的末端;针对每条第二光纤,若光信号在该第二光纤中的传输方向与X轴正方向相反,则将该第二光纤的光信号输入端作为该第二光纤的末端,光信号输出端作为该第二光纤的始端,若光信号在该第二光纤中的传输方向与X轴正方向相同,则将该第二光纤的光信号输入端作为该第二光纤的始端,光信号输出端作为该第二光纤的末端;
步骤S132、该处理装置对该瑞利散射信号进行处理,获得各条第一光纤和第二光纤上各个位置处的波长漂移值;
步骤S133、从始端到末端,对各条第一光纤和第二光纤上各个位置处的波长漂移值进行重排;
步骤S134、根据各条第一光纤和第二光纤上重排后的波长漂移值,计算出各条第一光纤和第二光纤上各个位置处的曲率。
本发明通过对各条第一光纤和第二光纤的始端和末端进行确定,并从始端到末端,对各条第一光纤和第二光纤上各个位置处的波长漂移值进行重排,可以使各条第一光纤和第二光纤上各个位置处的波长漂移值对应沿着Y轴正方向和X轴正方向排序,由此可以便于后续对柔性面板的上下表面曲面进行结合处理。
其中,所述步骤S132包括:该处理装置对该瑞利散射信号进行频域转换和互相关运算,获得各条第一光纤和第二光纤上各个位置处的波长漂移值;所述步骤S134包括:根据以下公式计算出各条第一光纤和第二光纤上各个位置处的曲率c:c=Δλ/κ=Δλ/(1-pe)λch
其中Δλ表示第一光纤或第二光纤对应位置处的波长漂移值,κ表示应变系数,pe表示光纤的弹光系数,λc表示中心波长,h表示柔性面板的厚度。
具体地,波长漂移与应力的对应关系为:
Δλ=λc(1-pe)ε;
式中,Δλ表示为应变ε引起波长漂移,λc表示中心波长,pe表示光纤的弹光系数;其中ε=h/R,h表示待测柔性面板的厚度,R表示曲率半径;
弯曲曲率与曲率半径对应的关系为:c=1/R;
通过求解方程组,可以得到曲率与波长漂移的关系式为:
c=Δλ/κ=Δλ/(1-pe)λch;
其中,Δλ表示的是波长漂移量,κ=(1-pe)λch,定义为应变系数。
另外,所述步骤S140中根据各条Y-Z轴二维曲线构建出该柔性面板的下表面曲面包括:采用二维插值算法对该Y-Z轴二维曲线进行二维插值,获得该下表面曲面;所述步骤S140中根据各条X-Z轴二维曲线构建出该柔性面板的上表面曲面包括:采用二维插值算法对该X-Z轴二维曲线进行二维插值,获得该上表面曲面。本发明通过采用二维插值算法对该Y-Z轴二维曲线进行二维插值,获得该下表面曲面,采用二维插值算法对该X-Z轴二维曲线进行二维插值,获得该上表面曲面,采用二维插值的方式,在柔性曲面相邻光纤之间插入多条曲线,提高相邻两条光纤的空间分辨率,可以提高柔性面板的上下表面曲面构建准确度,提高形变测量的准确度和空间分辨率。
所述步骤S140包括根据以下向量参数方程构建下表面各条第一光纤变形后的二维曲线:r(s)=y(s)j+z(s)k,其中y(s)表示下表面二维曲线上各点的Y轴坐标函数,z(s)表示下表面二维曲线上各点的Z轴坐标函数,s表示下表面各条第一光纤长度的参量,0≤s≤Ly,Ly是下表面各条第一光纤的总长度;
该下表面二维曲线上任意位置对应的点表示为:
式中,yn+1表示该点的Y轴坐标值,yn表示沿Y轴正方向,对于yn+1对应点,其上一点的Y轴坐标值,zn+1表示该点的Z轴坐标值,zn表示该上一点的Z轴坐标值,ds表示二维曲线的微小弧度;
式中,cy表示该点(该点在对应Y-Z轴二维曲线上)的曲率;
根据以下向量参数方程构建上表面各条第二光纤变形后的二维曲线:r(s)=x(s)i+z(s)k,其中x(s)表示上表面二维曲线上各点的X轴坐标函数,z(s)表示上表面二维曲线上各点的Z轴坐标函数,s表示上表面各条第二光纤长度的参量,0≤s≤Lx,Lx是上表面各条第二光纤的总长度;
该上表面二维曲线上任意位置对应的点表示为:
式中,xn+1表示该点的X轴坐标值,xn表示沿X轴正方向,其上一点的X轴坐标值,zn+1表示该点的Z轴坐标值,zn表示该上一点的Z轴坐标值,ds表示二维曲线的微小弧度;
式中,cx表示该点(该点在对应X-Z轴二维曲线上)的曲率。结合图3所示,ds可以等于两相邻点之间的连线长度。
所述步骤S150包括:从下表面曲面和上表面曲面中查找到对应的点P(x,y,z1)和Q(x,y,z2),使z=z1+z2,将点(x,y,z)作为该柔性面板的三维曲面上的对应点。本发明尤其适用于正常使用时处于弯曲状态的柔性面板。本发明的空间分辨率为毫米至厘米量级。另外,举例说明,在柔性面板上取一个40cm*40cm的方形区域铺设,下表面x轴每间隔5cm铺设一条光纤,y轴空间分辨率为1cm,一共40个传感点,计算相应的z1;上表面y轴每隔5cm铺设一条光纤,x方向的空间分辨率为1cm,一共40个传感点,计算相应的z2;将上下表面均二维插值,分别求得插值之后的Z1和Z2轴坐标,然后进行相加,即为Z轴坐标,此时插值计算之后的X,Y,Z三个坐标绘制出的三维效果图如图4所示。
本领域技术人员在考虑说明书及实践这里公开的发明后,将容易想到本发明的其它实施方案。本申请旨在涵盖本发明的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本发明的一般性原理并包括本发明未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本发明的真正范围和精神由下面的权利要求指出。
应当理解的是,本发明并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。本发明的范围仅由所附的权利要求来管制。
Claims (9)
1.一种基于光纤瑞利散射的柔性面板形状构造方法,其特征在于,包括:
步骤S110、针对待测的柔性面板,选定该柔性面板的上表面的一点作为零点,将平铺的柔性面板上表面上通过该零点的一条线作为X轴,将平铺的柔性面板上表面上通过该零点并与该X轴垂直的线作为Y轴,将通过该零点并垂直于该平铺的柔性面板的线作为Z轴,从而建立三维坐标系;
步骤S120、在平铺的柔性面板的下表面和上表面对应铺设多条平行于Y轴的第一光纤和多条平行于X轴的第二光纤,各条相邻第一光纤之间以及各条相邻第二光纤之间相连,该下表面的一条第一光纤的自由端作为下表面光信号出口,与该上表面的一条第二光纤的自由端相连,从而形成分布式光纤传感网络,每条第一光纤和第二光纤上的各个位置都对应于唯一的X-Y轴坐标,且针对每条第一光纤,其都与对应第二光纤垂直正交,该下表面的另一条第一光纤的自由端作为下表面光信号入口,与光频域反射仪OFDR(Optical FrequencyDomain Reflectometer)系统相连;
步骤S130、在测量该柔性面板的形变时,所述OFDR系统向该下表面光信号入口输入光信号,以使所述光信号在该分布式光纤传感网络的传输过程中产生瑞利散射信号,该瑞利散射信号向后传输给所述OFDR系统,所述OFDR系统与处理装置相连,并将该瑞利散射信号传输给该处理装置,该处理装置对该瑞利散射信号进行处理,获得各条第一光纤和第二光纤上各个位置处的曲率;
步骤S140、该处理装置根据各条第一光纤上各个位置处的曲率,构建出对应条数的Y-Z轴二维曲线,根据各条Y-Z轴二维曲线构建出该柔性面板的下表面曲面,并且根据各条第二光纤上各个位置处的曲率,构建出对应条数的X-Z轴二维曲线,根据各条X-Z轴二维曲线构建出该柔性面板的上表面曲面;
步骤S150、所述处理装置对该柔性面板的下表面曲面和上表面曲面进行结合处理,获得该柔性面板的三维曲面。
2.根据权利要求1所述的基于光纤瑞利散射的柔性面板形状构造方法,其特征在于,所述步骤S140中根据各条Y-Z轴二维曲线构建出该柔性面板的下表面曲面包括:采用二维插值算法对该Y-Z轴二维曲线进行二维插值,获得该下表面曲面;
所述步骤S140中根据各条X-Z轴二维曲线构建出该柔性面板的上表面曲面包括:采用二维插值算法对该X-Z轴二维曲线进行二维插值,获得该上表面曲面。
3.根据权利要求1或2所述的基于光纤瑞利散射的柔性面板形状构造方法,其特征在于,所述步骤S120中该上表面上与该下表面光信号出口相连的一条第二光纤的自由端为上表面光信号入口,选取所述零点为该下表面光信号入口,X轴上的一点(x0,0)为该下表面光信号出口,该零点为该上表面光信号入口,x0为非零值。
4.根据权利要求3所述的基于光纤瑞利散射的柔性面板形状构造方法,其特征在于,所述步骤S130中该处理装置对该瑞利散射信号进行处理,获得各条第一光纤和第二光纤上各个位置处的曲率包括:
步骤S131、针对每条第一光纤,若光信号在该第一光纤中的传输方向与Y轴正方向相反,则将该第一光纤的光信号输入端作为该第一光纤的末端,光信号输出端作为该第一光纤的始端,若光信号在该第一光纤中的传输方向与Y轴正方向相同,则将该第一光纤的光信号输入端作为该第一光纤的始端,光信号输出端作为该第一光纤的末端;针对每条第二光纤,若光信号在该第二光纤中的传输方向与X轴正方向相反,则将该第二光纤的光信号输入端作为该第二光纤的末端,光信号输出端作为该第二光纤的始端,若光信号在该第二光纤中的传输方向与X轴正方向相同,则将该第二光纤的光信号输入端作为该第二光纤的始端,光信号输出端作为该第二光纤的末端;
步骤S132、该处理装置对该瑞利散射信号进行处理,获得各条第一光纤和第二光纤上各个位置处的波长漂移值;
步骤S133、从始端到末端,对各条第一光纤和第二光纤上各个位置处的波长漂移值进行重排;
步骤S134、根据各条第一光纤和第二光纤上重排后的波长漂移值,计算出各条第一光纤和第二光纤上各个位置处的曲率。
5.根据权利要求4所述的基于光纤瑞利散射的柔性面板形状构造方法,其特征在于,所述步骤S132包括:该处理装置对该瑞利散射信号进行频域转换和互相关运算,获得各条第一光纤和第二光纤上各个位置处的波长漂移值;所述步骤S134包括:根据以下公式计算出各条第一光纤和第二光纤上各个位置处的曲率c:c=Δλ/κ=Δλ/(1-pe)λch
其中Δλ表示第一光纤或第二光纤对应位置处的波长漂移值,κ表示应变系数,pe表示光纤的弹光系数,λc表示中心波长,h表示柔性面板的厚度。
6.根据权利要求1所述的基于光纤瑞利散射的柔性面板形状构造方法,其特征在于,所述步骤S140包括根据以下向量参数方程构建下表面各条第一光纤变形后的二维曲线:r(s)=y(s)j+z(s)k,其中y(s)表示下表面二维曲线上各点的Y轴坐标函数,z(s)表示下表面二维曲线上各点的Z轴坐标函数,s表示下表面各条第一光纤长度的参量,0≤s≤Ly,Ly是下表面各条第一光纤的总长度;
该下表面二维曲线上任意位置对应的点表示为:
式中,yn+1表示该点的Y轴坐标值,yn表示沿Y轴正方向,其上一点的Y轴坐标值,zn+1表示该点的Z轴坐标值,zn表示该上一点的Z轴坐标值,ds表示二维曲线的微小弧度;
式中,cy表示该点的曲率;
根据以下向量参数方程构建上表面各条第二光纤变形后的曲面:r(s)=x(s)i+z(s)k,其中x(s)表示上表面二维曲线上各点的X轴坐标函数,z(s)表示上表面二维曲线上各点的Z轴坐标函数,s表示上表面各条第二光纤长度的参量,0≤s≤Lx,Lx是上表面各条第二光纤的总长度;
该上表面二维曲线上任意位置对应的点表示为:
式中,xn+1表示该点的X轴坐标值,xn表示沿X轴正方向,其上一点的X轴坐标值,zn+1表示该点的Z轴坐标值,zn表示该上一点的Z轴坐标值,ds表示二维曲线的微小弧度;
式中,cx表示该点的曲率。
7.根据权利要求5所述的基于光纤瑞利散射的柔性面板形状构造方法,其特征在于,所述ds等于两相邻点之间的连线长度。
8.根据权利要求1所述的基于光纤瑞利散射的柔性面板形状构造方法,其特征在于,所述步骤S150包括:从下表面曲面和上表面曲面中查找到对应的点P(x,y,z1)和Q(x,y,z2),使z=(z1+z2),将点(x,y,z)作为该柔性面板的三维曲面上的对应点。
9.根据权利要求1所述的基于光纤瑞利散射的柔性面板形状构造方法,其特征在于,所述柔性面板在正常使用时处于弯曲状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010193429.2A CN111223173A (zh) | 2020-03-18 | 2020-03-18 | 基于光纤瑞利散射的柔性面板形状构造方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010193429.2A CN111223173A (zh) | 2020-03-18 | 2020-03-18 | 基于光纤瑞利散射的柔性面板形状构造方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111223173A true CN111223173A (zh) | 2020-06-02 |
Family
ID=70826928
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010193429.2A Pending CN111223173A (zh) | 2020-03-18 | 2020-03-18 | 基于光纤瑞利散射的柔性面板形状构造方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111223173A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115950372A (zh) * | 2023-03-13 | 2023-04-11 | 山东省科学院激光研究所 | 一种分布式多维传感光纤形状传感测量方法及系统 |
CN117490734A (zh) * | 2023-12-21 | 2024-02-02 | 之江实验室 | 一种同步监测温度和三轴应变的光纤传感网络系统 |
US11971310B2 (en) * | 2018-06-22 | 2024-04-30 | Commissariat A L'energie Atomique Et Energies Alternatives | Method for producing an optical fibre for a distributed measurement of temperature or deformation in a harsh environment using the Rayleigh backscattered signal |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5357813A (en) * | 1991-10-17 | 1994-10-25 | The State Of Israel, Ministry Of Defence | Pressure sensor |
CN101050948A (zh) * | 2007-05-11 | 2007-10-10 | 上海大学 | 太空柔性帆板结构形态感知与可视化方法和装置 |
CN103339467A (zh) * | 2011-01-28 | 2013-10-02 | 皇家飞利浦电子股份有限公司 | 用于确定3d形状的光纤传感器 |
CN104949628A (zh) * | 2015-04-14 | 2015-09-30 | 上海大学 | 基于二维正交曲率的柔性板状结构复杂形态重建方法 |
CN109186490A (zh) * | 2018-09-21 | 2019-01-11 | 北京航空航天大学 | 一种基于多芯光纤的三维形状传感测量方法及装置 |
CN109307482A (zh) * | 2018-12-05 | 2019-02-05 | 重庆大学 | 温度不敏感的光纤形状重构系统 |
CN109613005A (zh) * | 2018-12-20 | 2019-04-12 | 武汉隽龙科技股份有限公司 | 基于ofdr的损伤检测方法 |
-
2020
- 2020-03-18 CN CN202010193429.2A patent/CN111223173A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5357813A (en) * | 1991-10-17 | 1994-10-25 | The State Of Israel, Ministry Of Defence | Pressure sensor |
CN101050948A (zh) * | 2007-05-11 | 2007-10-10 | 上海大学 | 太空柔性帆板结构形态感知与可视化方法和装置 |
CN103339467A (zh) * | 2011-01-28 | 2013-10-02 | 皇家飞利浦电子股份有限公司 | 用于确定3d形状的光纤传感器 |
CN104949628A (zh) * | 2015-04-14 | 2015-09-30 | 上海大学 | 基于二维正交曲率的柔性板状结构复杂形态重建方法 |
CN109186490A (zh) * | 2018-09-21 | 2019-01-11 | 北京航空航天大学 | 一种基于多芯光纤的三维形状传感测量方法及装置 |
CN109307482A (zh) * | 2018-12-05 | 2019-02-05 | 重庆大学 | 温度不敏感的光纤形状重构系统 |
CN109613005A (zh) * | 2018-12-20 | 2019-04-12 | 武汉隽龙科技股份有限公司 | 基于ofdr的损伤检测方法 |
Non-Patent Citations (2)
Title |
---|
张俊康等: "柔性蒙皮形状监测的光纤光栅传感及重构方法", 《传感器与微系统》 * |
郭艳丽等: "高密度分布式光纤传感技术在FRP复合材料结构健康监测中的应用", 《复合材料学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11971310B2 (en) * | 2018-06-22 | 2024-04-30 | Commissariat A L'energie Atomique Et Energies Alternatives | Method for producing an optical fibre for a distributed measurement of temperature or deformation in a harsh environment using the Rayleigh backscattered signal |
CN115950372A (zh) * | 2023-03-13 | 2023-04-11 | 山东省科学院激光研究所 | 一种分布式多维传感光纤形状传感测量方法及系统 |
CN117490734A (zh) * | 2023-12-21 | 2024-02-02 | 之江实验室 | 一种同步监测温度和三轴应变的光纤传感网络系统 |
CN117490734B (zh) * | 2023-12-21 | 2024-04-05 | 之江实验室 | 一种同步监测温度和三轴应变的光纤传感网络系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111223173A (zh) | 基于光纤瑞利散射的柔性面板形状构造方法 | |
US7813599B2 (en) | Method and apparatus for shape and end position determination using an optical fiber | |
CN111023991A (zh) | 一种基于多芯光纤的三维形状感知监测方法和系统 | |
CN109186490A (zh) | 一种基于多芯光纤的三维形状传感测量方法及装置 | |
US8746076B2 (en) | Shape sensing using a multi-core optical fiber having an arbitrary initial shape in the presence of extrinsic forces | |
Floris et al. | Measurement uncertainty of multicore optical fiber sensors used to sense curvature and bending direction | |
CN104807414A (zh) | 基于分布式光纤传感技术的地铁隧道沉降变形监测方法 | |
CN110069832B (zh) | 基于几何坐标转换的非连续板结构变形反演与拼接方法 | |
CN107917717B (zh) | 具有光学变形传感器的飞行参数测量装置和相应测量方法 | |
CN106932337B (zh) | 多光纤光栅钢筋腐蚀传感器系统的数据处理方法 | |
CN111982000A (zh) | 一种基于Beta标架的光纤形状重构方法及装置 | |
Moore | Shape sensing using multi-core fiber | |
CN110006563A (zh) | 直升机桨叶挥舞与摆振弯矩分布式光纤解耦测量方法 | |
CN107389284A (zh) | 一种基于应变的框架结构弹性变形的测量方法 | |
CN110532708A (zh) | 一种残余应力预测方法、装置、服务器及存储介质 | |
KR101498381B1 (ko) | 광섬유 브래그 격자 센서를 이용한 파이프 구조의 3차원 형상 모니터링 시스템 | |
CN110044682B (zh) | 基于fbg传感器的单边缺口铝合金试件疲劳裂纹扩展监测方法 | |
CN110887448B (zh) | 一种基于光纤应变测量的梁结构形态重构方法 | |
CN208984056U (zh) | 一种实时自动测量隧道断面变形的装置 | |
Derkevorkian et al. | Computational studies of a strain-based deformation shape prediction algorithm for control and monitoring applications | |
KR20050062308A (ko) | 광섬유 센서를 이용한 터널의 2차원 내공변위 측정장치 및시스템 | |
JP7376052B2 (ja) | 形状測定システム及び形状測定方法 | |
Rogge et al. | Shape sensing using a multi-core optical fiber having an arbitrary initial shape in the presence of extrinsic forces | |
CN112433306B (zh) | 分布式光纤三维形状快速重构装置 | |
CN214893116U (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200602 |
|
RJ01 | Rejection of invention patent application after publication |