CN112985310A - 一种挖泥船耙管空间角度测量方法 - Google Patents

一种挖泥船耙管空间角度测量方法 Download PDF

Info

Publication number
CN112985310A
CN112985310A CN202110381675.5A CN202110381675A CN112985310A CN 112985310 A CN112985310 A CN 112985310A CN 202110381675 A CN202110381675 A CN 202110381675A CN 112985310 A CN112985310 A CN 112985310A
Authority
CN
China
Prior art keywords
hose
coordinate system
curvature
axis
fiber grating
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
Application number
CN202110381675.5A
Other languages
English (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.)
Guangzhou Maritime University
Original Assignee
Guangzhou Maritime University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangzhou Maritime University filed Critical Guangzhou Maritime University
Priority to CN202110381675.5A priority Critical patent/CN112985310A/zh
Publication of CN112985310A publication Critical patent/CN112985310A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • 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/26—Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes
    • 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/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
    • G01B11/255—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures for measuring radius of curvature

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明提供了一种挖泥船耙管空间角度测量方法，在耙管的每个软管上沿其周向均匀设有四个传感器装置，传感器装置包括设置在软管上记忆合金丝和若干个粘贴于记忆合金丝上光纤光栅传感器，光纤光栅传感器通过光纤串接在一起，其中方法包括：确定软管的软管测量节点，由第一预设计算规则得到软管测量节点的第一曲率标量和第二曲率标量，在软管上建立静坐标系{B}，将软管等效为静坐标系中的目标曲线，由第二预设计算规则计算得到每个微小圆弧等效的矢量，由预设公式计算得到目标曲线末端在静坐标系中的空间角度，以目标曲线末端在静坐标系中的空间角度作为耙管的空间角度。本发明利用光纤光栅传感器网络测量挖泥船耙管的空间角度，灵敏度高，寿命长。

Description

一种挖泥船耙管空间角度测量方法
技术领域
本发明属于挖泥船技术领域，尤其涉及一种挖泥船耙管空间角度测量方法。
背景技术
目前测量挖泥船耙管的孔径角度的方法主要有两种：一是测量绞车钢丝绳下放长度，二是基于计算机控制的角度传感器。而挖泥船施工过程中耙管一般都处于以下几种工况中：(1)挖泥船耙管的耙头工作在20-50米深的水下；(2)挖泥船航行时泥泵与高压水枪冲泥沙带来的振动与摩擦；(3)耙管实用钢缆悬空吊下容易遭受浪涌冲击摆动。在上述的工况下，由于角度传感器为电子传感器，容易受到损伤和失灵，且附在耙管上的电缆在耙管受到浪涌冲击、耙管振动而发生位移磨损。一旦水下电缆有破损口，在水压浸透下，海水顺着电线进入到传感器、分线盒等与电缆有链接的电气设备进而使电路短路，使敏感元件失灵，导致不能准确测量耙管的空间角度，使得疏浚的工作效率低，电子传感器在水下恶劣环境下的寿命并不长，常常需要维护或者更换。
发明内容
本发明的目的在于提供一种挖泥船耙管空间角度测量方法，利用光纤光栅传感器网络测量挖泥船耙管的空间角度，灵敏度高，寿命长。
本发明是通过以下技术方案实现的：
一种挖泥船耙管空间角度测量方法，包括耙管，耙管包括弯管、中段耙管和耙头，弯管通过一软管与中段耙管的一端连接，中段耙管的另一端通过另一软管与耙头的一端连接；每个软管上沿其周向均匀设有四个传感器装置，传感器装置包括记忆合金丝和若干个具有不同中心波长的光纤光栅传感器，记忆合金丝设置在软管上，且其长度方向与软管的长度方向一致，光纤光栅传感器粘贴于记忆合金丝上并通过光纤串接在一起；
方法包括以下步骤：
以软管上设有光纤光栅传感器的圆截面为软管测量节点，对于每个软管测量节点，根据第一预设计算规则计算得到软管测量节点的第一曲率标量和第二曲率标量，其中第一曲率标量为软管测量节点处的空间曲率在第一方向上的曲率分量的标量，第二曲率标量为软管测量节点处的空间曲率在第二方向上的曲率分量的标量，第一方向和第二方向互相垂直；
在软管上建立静坐标系{B}，静坐标系{B}中的原点O位于软管第一端的圆心；
将软管等效为在静坐标系{B}中的目标曲线，目标曲线的起始端为静坐标系{B}的原点，并将目标曲线分为多段依次连接的微小圆弧,当微小圆弧的长度取值趋向于0时，微小圆弧等效为该微小圆弧起始端指向其末端的矢量，基于每个软管测量节点的第一曲率标量和第二曲率标量，根据第二预设计算规则计算得到每个微小圆弧等效的矢量；
基于每个微小圆弧等效的矢量，根据预设公式计算得到目标曲线末端在静坐标系{B}中的空间角度，以目标曲线末端在静坐标系{B}中的空间角度作为耙管的空间角度。
进一步地，根据第一预设计算规则计算得到第一曲率标量和第二曲率标量的步骤包括：
在软管测量节点上建立平面坐标系O-YZ，其中平面坐标系的原点O为软管测量节点的圆心，Y轴和Z轴分别穿过位于软管测量节点上的两相对光纤光栅传感器，第一方向为Y轴的正方向，第二方向为Z轴的正方向；
通过公式(1)计算第一曲率标量，通过公式(2)计算第二曲率标量：
Figure BDA0003013216220000021
Figure BDA0003013216220000022
其中，k1为第一曲率标量，k2为第二曲率标量，Δλ1为位于Y轴正方向上的光纤光栅传感器的波长偏移量，Δλ2为位于Z轴负方向上的光纤光栅传感器的波长偏移量，Δλ3为位于Y轴负方向上的光纤光栅传感器的波长偏移量，Δλ4为位于Z轴正方向上的光纤光栅传感器的波长偏移量，a1和a2均为系数。
进一步地，公式(1)和公式(2)的具体推导过程为：
每个光纤光栅传感器的轴向应变引起中心波长漂移量与其轴向应变满足：
Figure BDA0003013216220000023
其中，Δλ为光纤光栅传感器的中心波长漂移量，λ为光纤光栅传感器的初始中心波长，Pe为光纤光栅传感器的有效弹光系数，εz为光纤光栅传感器的轴向应变；
根据材料力学，对于圆截面弹性梁发生纯弯曲时，其轴向应变εz与曲率之间的关系有：
Figure BDA0003013216220000024
其中，r为光纤光栅传感器中心到平面坐标系O-YZ原点的距离，R为测量点的曲率半径，k为测量点的曲率；
当软管弯曲时，每个光纤光栅传感器的轴向应变与曲率之间满足：
Figure BDA0003013216220000025
其中，εz,i为软管测量节点上的第i个光纤光栅传感器的轴向应变，ri为软管测量节点上的第i个光纤光栅传感器到平面坐标系O-YZ原点的距离，θb为软管弯曲方向与Y轴的夹角，θi为软管测量节点上的第i个光纤光栅传感器与Y轴之间的角度，i的取值为i＝1、2、3、4；
定义一个由平面坐标系O-YZ原点指向一个光纤光栅传感器的表观曲率向量
Figure BDA0003013216220000031
则
Figure BDA0003013216220000032
满足：
Figure BDA0003013216220000033
其中，
Figure BDA0003013216220000034
为沿Y轴的单位向量，
Figure BDA0003013216220000035
为沿Z轴的单位向量；
所有表观曲率向量之和为：
Figure BDA0003013216220000036
每个光纤光栅传感器中心到平面坐标系O-YZ原点的距离均为常量r，把θi带入公式(7)得到：
Figure BDA0003013216220000037
则软管曲率为：
Figure BDA0003013216220000038
其中，λB为位于Y轴上的两个光纤光栅传感器的初始中心波长，λA为位于Z轴上的两个光纤光栅传感器的初始中心波长；
软管曲率在Y轴上的分向量为：
Figure BDA0003013216220000039
由公式(10)得到第一曲率标量
Figure BDA00030132162200000310
软管曲率在Z轴上的分向量为：
Figure BDA00030132162200000311
由公式(11)得到第二曲率标量
Figure BDA00030132162200000312
进一步地，静坐标系{B}中X轴垂直于软管第一端的端面，Y轴和Z轴位于软管第一端所在平面上，且两者分别穿过两相对的传感器装置；
基于每个软管测量节点的第一曲率标量和第二曲率标量，根据第二预设计算规则计算得到每个微小圆弧的等效矢量的步骤包括：
以静坐标系{B}的原点为原点，并以单位向量
Figure BDA00030132162200000313
为A轴、单位向量
Figure BDA00030132162200000314
为B轴以及单位向量
Figure BDA00030132162200000315
为C轴，建立动坐标系{A}，动坐标系{A}在初始状态下，其A轴、B轴和C轴分别与静坐标系{B}的Y轴、Z轴和X轴重合，动坐标系{A}用齐次矩阵表示为：
Figure BDA0003013216220000041
以目标曲线起始端到末端的顺序，依次选取一段微小圆弧作为目标微小圆弧，目标微小圆弧等效的矢量的计算过程如下：
步骤1、基于目标微小圆弧两侧的软管测量节点的第一曲率标量，由线性插值法计算得到目标微小圆弧的第一圆弧曲率标量k′1，基于目标微小圆弧两侧的软管测量节点的第二曲率标量，由线性插值法得到目标微小圆弧的第二圆弧曲率标量k′2；
步骤2、判断目标微小圆弧是否为首段微小圆弧；
步骤3、若是首段微小圆弧，则执行步骤4至步骤9，若不是首段微小圆弧，则执行步骤10至步骤11；
步骤4、获取初始状态的动坐标系{A}，以初始状态的动坐标系{A}作为目标微小圆弧的目标动坐标系；
步骤5、基于获取的目标动坐标系，根据公式(13)合成目标微小圆弧的空间曲率：
Figure BDA0003013216220000042
其中，
Figure BDA0003013216220000043
为目标动坐标系中A轴的单位向量，其方向指向A轴的正方向，
Figure BDA0003013216220000044
为目标动坐标系中B轴，其方向指向B轴的正方向；
步骤6、根据公式(14)计算得到目标微小圆弧的空间曲率半径：
Figure BDA0003013216220000045
其中，D1为目标微小圆弧的起始端；
步骤7、在静坐标系{B}中，
Figure BDA0003013216220000046
由
Figure BDA0003013216220000047
绕平面D1RD2的单位法向量
Figure BDA00030132162200000414
旋转β角得到，则
Figure BDA0003013216220000048
满足：
Figure BDA0003013216220000049
其中，D2为目标微小圆弧的末端，
Figure BDA00030132162200000410
Δs为圆弧的长度；
步骤8、根据公式(16)计算得到目标微小圆弧等效的矢量：
Figure BDA00030132162200000411
步骤9、获取前一段微小圆弧的目标动坐标系，根据公式(17)计算获得的目标动坐标系绕平面D1RD2的单位法向量
Figure BDA00030132162200000412
旋转β角得到的动坐标系{A′m}，以动坐标系{A′m}作为目标微小圆弧的目标动坐标系：
Figure BDA00030132162200000413
步骤10、执行步骤5至步骤8。
进一步地，根据预设公式计算得到目标曲线末端在静坐标系{B}中的空间角度的步骤中，预设公式为：
Figure BDA0003013216220000051
进一步地，光纤光栅传感器采用环氧树脂粘贴于记忆合金丝上。
进一步地，记忆合金丝为镍钛记忆合金丝。
相比于现有技术，本发明的有益效果为：光纤光栅传感器无需电信号传输电路，无需现场供电，体积小，不受电磁干扰，可靠性高，灵敏度强，寿命长，具有多路传输能力，与光纤传输系统很好地兼容，大量光纤光栅传感器将数据通过多路传送构成一个光纤传感网络，利用光纤传感网络测量挖泥船耙管的空间角度，使得耙管的工作角度都处于要求的角度，让挖泥船保持在理想的工作状态，不仅提高疏浚工作的效率，同时可以避免施工人员的视觉盲区，保证设备和人员的安全。
附图说明
图1为本发明挖泥船耙管空间角度测量方法中耙管的结构示意图；
图2为本发明挖泥船耙管空间角度测量方法中耙管的侧视图；
图3为本发明挖泥船耙管空间角度测量方法中传感器装置安装在软管上的示意图；
图4为本发明挖泥船耙管空间角度测量方法中平面坐标系的示意图；
图5为本发明挖泥船耙管空间角度测量方法中曲率计算等效模型图；
图6为本发明挖泥船耙管空间角度测量方法中建立静坐标系{B}的示意图；
图7为本发明挖泥船耙管空间角度测量方法中软管在静坐标系{B}中的等效示意图；
图8为本发明挖泥船耙管空间角度测量方法中曲线拟合示意。
图中，1-弯管，2-中段耙管，3-耙头，4-软管，5-传感器装置，51-记忆合金丝，52-光纤光栅传感器，53-光纤。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。
因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
在本发明的描述中，需要说明的是，术语“上”、“下”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，或者是该发明产品使用时惯常摆放的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。
请参阅图1至图3，图1为本发明挖泥船耙管空间角度测量方法中耙管的结构示意图，图2为本发明挖泥船耙管空间角度测量方法中耙管的侧视图，图3为本发明挖泥船耙管空间角度测量方法中传感器装置安装在软管上的示意图。一种挖泥船耙管空间角度测量方法，包括耙管，耙管包括弯管1、中段耙管2和耙头3，弯管1通过一软管4与中段耙管2的一端连接，中段耙管2的另一端通过另一软管4与耙头3的一端连接；每个软管4上沿其周向均匀设有四个传感器装置5，传感器装置5包括记忆合金丝51和若干个具有不同中心波长的光纤光栅传感器52，记忆合金丝51设置在软管4上，且其长度方向与软管4的长度方向一致，光纤光栅传感器52粘贴于记忆合金丝51上并通过光纤53串接在一起；
耙管为现有设备，其具体结构可以参考现有耙管的结构，耙管的弯管1和中段耙管2之间设有一字关节，一字关节保证弯管1和中段耙管2之间连接的软管4仅能进行水平弯曲，中段耙管2与耙头3之间设有万向节，万向节保证中段耙管2与耙头3之间连接的软管4能在水平方向和垂直方向均可以弯曲。在每个软管4上安装有四个传感器装置5组成光纤传感网络，每个传感器装置5包括若干个具有不同中心波长的光纤光栅传感器52，将该若干个具有不同中心波长的光纤光栅传感器52间隔粘贴在一根记忆合金丝51上，整体形成一个传感器装置5。由于软管4表面不平整，其曲率变化在软管4表面并不能灵敏地反映，同时，软管4的体积大，不方便进行标定实验，也不方便对光纤光栅传感器52进行保护，因此将若干光纤光栅传感器52与记忆合金丝51整体形成一个传感器装置5，避免了将光纤光栅传感器52直接粘贴在软管4上时可能会因软管4是柔性材料导致应变分布不均产生啁啾效应，使光纤光栅解调仪无法解调出光纤光栅传感器52的中心波长，及避免光纤光栅传感器52可能长期处于应力状态下产生形变而缩短光纤光栅传感器52的工作寿命。同时光纤光栅传感器52无需电信号传输电路，无需现场供电，体积小，不受电磁干扰，可靠性高，灵敏度强，寿命长，具有多路传输能力，与光纤传输系统很好地兼容，可以灵敏且准确地检测到记忆合金丝51的曲率变化，由于记忆合金丝51弯曲情况跟软管4情况相同，从而利用光纤光栅传感器52和记忆合金丝51封装成传感器装置5能够实时获取耙管软管4的曲率，并且光纤光栅传感器52对轴向应变灵敏，提高其测量软管4空间曲率的精度和可靠性。优选地，各光纤光栅传感器52采用环氧树脂粘贴于记忆合金丝51上，进一步地，环氧树脂采用353ND环氧树脂胶水。在一实施例中，记忆合金丝51为镍钛记忆合金丝。在将四个传感器装置5均匀安装在软管4上之后，在软管4同一圆截面上的四个光纤光栅传感器52的波长相同。
基于上述装置，挖泥船耙管空间角度测量方法包括以下步骤：
S1、以软管4上设有光纤光栅传感器52的圆截面为软管4测量节点，对于每个软管4测量节点，根据第一预设计算规则计算得到软管4测量节点的第一曲率标量和第二曲率标量，其中第一曲率标量为软管4测量节点处的空间曲率在第一方向上的曲率分量的标量，第二曲率标量为软管4测量节点处的空间曲率在第二方向上的曲率分量的标量，第一方向和第二方向互相垂直；
S2、在软管4上建立静坐标系{B}，静坐标系{B}中的原点O位于软管4第一端的圆心；
S3、将软管4等效为在静坐标系{B}中的目标曲线，目标曲线的起始端为静坐标系{B}的原点，并将目标曲线分为多段依次连接的微小圆弧,当微小圆弧的长度取值趋向于0时，微小圆弧等效为该微小圆弧起始端指向其末端的矢量，基于每个软管4测量节点的第一曲率标量和第二曲率标量，根据第二预设计算规则计算得到每个微小圆弧等效的矢量；
S4、基于每个微小圆弧等效的矢量，根据预设公式计算得到目标曲线末端在静坐标系{B}中的空间角度，以目标曲线末端在静坐标系{B}中的空间角度作为耙管的空间角度。
请结合参阅图4，图4为本发明挖泥船耙管空间角度测量方法中平面坐标系的示意图。在上述步骤S1中，在每个软管4上均安装有四个传感器装置5，四个传感器装置5的安装方向均相同，软管4上的软管4测量节点的数量与传感器装置5包含的光纤光栅传感器52的数量相同，每个软管4测量节点上的四个光纤光栅传感器52的中心波长相同，且材料一致，如图4所示，在一实施例中，每个传感器装置5包括四个光纤光栅传感器52，则软管4上有四个软管4测量节点。获取软管4测量节点上的四个光纤光栅传感器52的中心波长，则可以根据第一预设计算规则计算得到软管4测量节点的第一曲率标量和第二曲率标量，其中第一曲率标量为软管4测量节点处的空间曲率在第一方向上的曲率分量的标量，第二曲率标量为软管4测量节点处的空间曲率在第二方向上的曲率分量的标量，第一方向和第二方向垂直，且第一方向和第二方向分别穿过位于软管4测量节点上的相邻两光纤光栅传感器52。
进一步地，在步骤S1中，根据第一预设计算规则计算得到第一曲率标量和第二曲率标量的步骤包括：
S11、在软管4测量节点上建立平面坐标系O-YZ，其中平面坐标系的原点O为软管4测量节点的圆心，Y轴和Z轴分别穿过位于软管4测量节点上的两相对光纤光栅传感器52，第一方向为Y轴的正方向，第二方向为Z轴的正方向；
S12、通过公式(1)计算第一曲率标量，通过公式(2)计算第二曲率标量：
Figure BDA0003013216220000081
Figure BDA0003013216220000082
其中，k1为第一曲率标量，k2为第二曲率标量，Δλ1为位于Y轴正方向上的光纤光栅传感器52的波长偏移量，Δλ2为位于Z轴负方向上的光纤光栅传感器52的波长偏移量，Δλ3为位于Y轴负方向上的光纤光栅传感器52的波长偏移量，Δλ4为位于Z轴正方向上的光纤光栅传感器52的波长偏移量，a1和a2均为系数。
请结合参阅图5，图5为本发明挖泥船耙管空间角度测量方法中曲率计算等效模型图。在上述步骤S11中，如图4所示，在软管4测量节点上建立平面坐标系O-YZ，由于软管4测量节点上的四个光纤光栅传感器52沿软管4的周向均匀设置，因此，四个光纤光栅传感器52在软管4测量节点上成90度等间分布，平面坐标系O-YZ的Y轴和Z轴可以分别穿过两个光纤光栅传感器52，得到软管4曲率计算模型，并且由于记忆合金丝51弯曲情况跟软管4情况相同，光纤光栅传感器52和记忆合金丝51可以看作整体，则得到如图5所示的软管4曲率计算等效模型，图中，WQ指示方向表示弯曲方向，ZX表示中性轴。
在上述步骤S12中，公式(1)和公式(2)可以预先设定，在获取光纤光栅传感器52的中心波长漂移量之后，直接计算出软管4测量节点处空间曲率在Y轴上的第一曲率标量，以及软管4测量节点处空间曲率在Z轴上的第二曲率标量。
进一步地，公式(1)和公式(2)的具体推导过程为：
每个光纤光栅传感器52的轴向应变引起中心波长漂移量与其轴向应变满足：
Figure BDA0003013216220000091
其中，Δλ为光纤光栅传感器52的中心波长漂移量，λ为光纤光栅传感器52的初始中心波长，Pe为光纤光栅传感器52的有效弹光系数，εz为光纤光栅传感器52的轴向应变；
根据材料力学，对于圆截面弹性梁发生纯弯曲时，其轴向应变εz与曲率之间的关系有：
Figure BDA0003013216220000092
其中，r为光纤光栅传感器52中心到平面坐标系O-YZ原点的距离，R为测量点的曲率半径，k为测量点的曲率，记忆合金丝可看作是圆截面弹性梁，其轴向应变与光纤光栅传感器轴向应变的比值等于光纤光栅传感器半径与记忆合金丝半径的比值；
当软管4如图5所示向WQ方向弯曲时，每个光纤光栅传感器52的轴向应变与曲率之间满足：
Figure BDA0003013216220000093
其中，εz,i为软管4测量节点上的第i个光纤光栅传感器52的轴向应变，ri为软管4测量节点上的第i个光纤光栅传感器52到平面坐标系O-YZ原点的距离，θb为软管4弯曲方向与Y轴的夹角，θi为软管4测量节点上的第i个光纤光栅传感器52与Y轴之间的角度，i的取值为i＝1、2、3、4，θ1＝0，θ2＝-pi/2，θ3＝-pi，θ4＝-3pi/2；具体地，在软管4测量节点上的四个光纤光栅传感器52中，以位于Y轴正方向的光纤光栅传感器52为第一光纤光栅传感器52，以位于Z轴负方向的光纤光栅传感器52为第二光纤光栅传感器52，以位于Y轴负方向的光纤光栅传感器52为第三光纤光栅传感器52，以位于Z轴正方向的光纤光栅传感器52为四个光纤光栅传感器52。
定义一个由平面坐标系O-YZ原点指向一个光纤光栅传感器52的表观曲率向量
Figure BDA0003013216220000094
则
Figure BDA0003013216220000095
满足：
Figure BDA0003013216220000096
其中，
Figure BDA0003013216220000097
为沿Y轴的单位向量，
Figure BDA0003013216220000098
为沿Z轴的单位向量；
所有表观曲率向量之和为：
Figure BDA0003013216220000099
每个光纤光栅传感器52中心到平面坐标系O-YZ原点的距离均为常量r，把θi带入公式(7)得到：
Figure BDA00030132162200000910
则软管4曲率为：
Figure BDA0003013216220000101
其中，λB为位于Y轴上的两个光纤光栅传感器52的初始中心波长，λA为位于Z轴上的两个光纤光栅传感器52的初始中心波长；
软管4曲率在Y轴上的分向量为：
Figure BDA0003013216220000102
传感器装置5的曲率和波长漂移量之间近似线性，因此由公式(10)得到第一曲率标量
Figure BDA0003013216220000103
软管4曲率在Z轴上的分向量为：
Figure BDA0003013216220000104
由公式(11)得到第二曲率标量
Figure BDA0003013216220000105
a1和a2均为计算出来的系数。
上述推导过程为理论关系，在计算之前，应先对传感器装置5进行标定，提高测量精度。
请结合参阅图6，图6为本发明挖泥船耙管空间角度测量方法中建立静坐标系{B}的示意图。在上述步骤S3中，在软管4上建立静坐标系{B}，用于搭建耙管三维空间运动模型。具体地，静坐标系{B}的原点O位于所述软管4第一端的圆心，X轴垂直于软管4第一端的端面，Y轴和Z轴位于软管4第一端所在平面上，且两者分别穿过两相对的传感器装置5，如图6所示。
请结合参阅图7，为本发明挖泥船耙管空间角度测量方法中。在上述步骤S4中，软管4弯曲后可以等效为在静坐标系{B}中的目标曲线，如图7所示，图中曲线c为目标曲线。目标曲线的起始端为静坐标系{B}的原点，基于微分原理，将目标曲线分为n段依次连接的微小圆弧，n为整数，每段微小圆弧的长度为Δs，以目标曲线起始端向末端方向，按顺序分别将微小圆弧记为
Figure BDA0003013216220000106
当微小圆弧的长度Δs取值趋向于0时，
Figure BDA0003013216220000107
等效为矢量
Figure BDA0003013216220000108
等效为矢量
Figure BDA0003013216220000109
等效为矢量
Figure BDA00030132162200001010
然后基于每个软管4测量节点的第一曲率标量和第二曲率标量，根据第二预设计算规则计算每个微小圆弧等效的矢量。
请结合参阅图8，图8为本发明挖泥船耙管空间角度测量方法中曲线拟合示意。在上述步骤S4中，基于每个软管4测量节点的第一曲率标量和第二曲率标量，根据第二预设计算规则计算得到每个微小圆弧的等效矢量的步骤包括：
S41、以静坐标系{B}的原点为原点，并以单位向量
Figure BDA0003013216220000111
为A轴、单位向量
Figure BDA0003013216220000112
为B轴以及单位向量
Figure BDA0003013216220000113
为C轴，建立动坐标系{A}，动坐标系{A}在初始状态下，其A轴、B轴和C轴分别与静坐标系{B}的Y轴、Z轴和X轴重合，动坐标系{A}用齐次矩阵表示为：
Figure BDA0003013216220000114
S42、以目标曲线起始端到末端的顺序，依次选取一段微小圆弧作为目标微小圆弧，目标微小圆弧等效的矢量的计算过程如下：
步骤1、基于目标微小圆弧两侧的软管4测量节点的第一曲率标量，由线性插值法计算得到目标微小圆弧的第一圆弧曲率标量k′1，基于目标微小圆弧两侧的软管4测量节点的第二曲率标量，由线性插值法得到目标微小圆弧的第二圆弧曲率标量k′2；
步骤2、判断目标微小圆弧是否为首段微小圆弧；
步骤3、若是首段微小圆弧，则执行步骤4至步骤9，若不是首段微小圆弧，则执行步骤10至步骤11；
步骤4、获取初始状态的动坐标系{A}，以初始状态的动坐标系{A}作为目标微小圆弧的目标动坐标系；
步骤5、基于获取的目标动坐标系，根据公式(13)合成目标微小圆弧的空间曲率：
Figure BDA0003013216220000115
其中，
Figure BDA0003013216220000116
为目标动坐标系中A轴的单位向量，其方向指向A轴的正方向，
Figure BDA0003013216220000117
为目标动坐标系中B轴，其方向指向B轴的正方向；
步骤6、根据公式(14)计算得到目标微小圆弧的空间曲率半径：
Figure BDA0003013216220000118
其中，D1为目标微小圆弧的起始端；
步骤7、在静坐标系{B}中，
Figure BDA0003013216220000119
由
Figure BDA00030132162200001110
绕平面D1RD2的单位法向量
Figure BDA00030132162200001115
旋转β角得到，则
Figure BDA00030132162200001111
满足：
Figure BDA00030132162200001112
其中，D2为目标微小圆弧的末端，
Figure BDA00030132162200001113
Δs为圆弧的长度；
步骤8、根据公式(16)计算得到目标微小圆弧等效的矢量：
Figure BDA00030132162200001114
步骤9、获取前一段微小圆弧的目标动坐标系，根据公式(17)计算获得的目标动坐标系绕平面D1RD2的单位法向量
Figure BDA0003013216220000121
旋转β角得到的动坐标系{A′m}，以动坐标系{A′m}作为目标微小圆弧的目标动坐标系：
Figure BDA0003013216220000122
步骤11、执行步骤5至步骤8。
在上述步骤S41中，以静坐标系{B}的原点为原点，以静坐标系{B}Y轴上的单位向量
Figure BDA0003013216220000123
为A轴、以Z轴上的单位向量
Figure BDA0003013216220000124
为B轴以及以X轴上的单位向量
Figure BDA0003013216220000125
为C轴，建立动坐标系{A}，因此动坐标系{A}在初始状态与静坐标系{B}重合，且
Figure BDA0003013216220000126
以静坐标系{B}和动坐标系{A}构建耙管的三位空间运动模型。由于
Figure BDA0003013216220000127
和
Figure BDA0003013216220000128
为单位向量，则动坐标系{A}用齐次矩阵表示为：
Figure BDA0003013216220000129
在上述步骤S42中，为计算
Figure BDA00030132162200001210
等效的矢量
Figure BDA00030132162200001211
等效的矢量
Figure BDA00030132162200001212
等效的矢量
Figure BDA00030132162200001213
需从首段矢量
Figure BDA00030132162200001214
开始计算，因此以目标曲线起始端到末端的顺序，依次选取一段微小圆弧作为目标微小圆弧，则第一次选取的目标微小圆弧为首段微小圆弧
Figure BDA00030132162200001215
首段微小圆弧
Figure BDA00030132162200001216
等效的首段矢量
Figure BDA00030132162200001217
的计算过程为：
首先，目标曲线的曲率变化不剧烈，因此基于首段微小圆弧
Figure BDA00030132162200001218
两侧的软管4测量节点的第一曲率标量，由线性插值法计算得到首段微小圆弧
Figure BDA00030132162200001219
的第一圆弧曲率标量k′11，基于首段微小圆弧
Figure BDA00030132162200001220
两侧的软管4测量节点的第二曲率标量，由线性插值法得到首段微小圆弧
Figure BDA00030132162200001221
的第二圆弧曲率标量k′21，线性插值法为现有计算方法，在此不赘述。
根据获取初始状态的动坐标系{A}，以动坐标系{A}作为首段微小圆弧
Figure BDA00030132162200001222
的目标动坐标系，即以初始状态的动坐标系{A}中A轴的单位向量
Figure BDA00030132162200001223
和以初始状态的动坐标系{A}中B轴的单位向量
Figure BDA00030132162200001224
分别确定第一圆弧曲率标量k′11和第二圆弧曲率标量k′21的方向，得到首段微小圆弧
Figure BDA00030132162200001225
的空间曲率
Figure BDA00030132162200001226
进而求得空间曲率半径
Figure BDA00030132162200001227
首段微小圆弧
Figure BDA00030132162200001228
的起始端为原点O。
在静坐标系{B}中，
Figure BDA00030132162200001229
由
Figure BDA00030132162200001230
绕平面ORQ1的单位法向量
Figure BDA00030132162200001231
旋转β角得到，单位法向量
Figure BDA00030132162200001232
当
Figure BDA00030132162200001233
的长度Δs趋向于0时，
Figure BDA00030132162200001234
根据空间坐标系的齐次变换知道点在空间直角坐标系中绕过原点的任意单位向量的一般变换的旋转算子为：
Figure BDA0003013216220000131
其中，versβ＝1-cosβ、cβ＝cosβ、sβ＝sinβ；
Figure BDA0003013216220000132
绕
Figure BDA0003013216220000133
旋转β角后得到
Figure BDA0003013216220000134
即
Figure BDA0003013216220000135
因此所求的首段矢量
Figure BDA0003013216220000136
为：
Figure BDA0003013216220000137
则首段微小圆弧
Figure BDA0003013216220000138
的末端Q1在静坐标系{B}上的位置为：
Figure BDA0003013216220000139
首段微小圆弧
Figure BDA00030132162200001310
等效的首段矢量
Figure BDA00030132162200001311
计算完之后，按目标曲线起始端到末端的顺序，取位于首段微小圆弧
Figure BDA00030132162200001312
后面的第二段微小圆弧，第二段微小圆弧
Figure BDA00030132162200001313
等效的第二段矢量
Figure BDA00030132162200001314
的计算过程为：
首先，基于第二段微小圆弧
Figure BDA00030132162200001315
两侧的软管4测量节点的第一曲率标量，由线性插值法计算得到第二段微小圆弧
Figure BDA00030132162200001316
的第一圆弧曲率标量k′12，基于第二段微小圆弧
Figure BDA00030132162200001317
两侧的软管4测量节点的第二曲率标量，由线性插值法得到第二段微小圆弧
Figure BDA00030132162200001318
的第二圆弧曲率标量k′22。
而在首段矢量
Figure BDA00030132162200001319
的计算过程中，
Figure BDA00030132162200001320
绕平面ORQ1的单位法向量
Figure BDA00030132162200001321
旋转β角后，第二段微小圆弧
Figure BDA00030132162200001322
的空间曲率分量的方向也应跟着旋转，所以需要知道首段微小圆弧
Figure BDA00030132162200001323
的目标动坐标系经过旋转后的位置，即初始状态的动坐标系{A}绕平面ORQ1的单位法向量
Figure BDA00030132162200001324
旋转β角得到的动坐标系{A′m}，以得到第二段微小圆弧
Figure BDA00030132162200001325
的目标动坐标系{A′m}：
Figure BDA00030132162200001326
根据获取的第二段微小圆弧
Figure BDA00030132162200001327
的目标动坐标系{A′m}，以目标动坐标系{A′m}中A轴上的单位向量
Figure BDA00030132162200001328
和以目标动坐标系{A′m}中B轴上的单位向量
Figure BDA00030132162200001329
分别确定第一圆弧曲率标量k′11和第二圆弧曲率标量k′21的方向，得到第二段微小圆弧
Figure BDA00030132162200001330
的空间曲率
Figure BDA00030132162200001331
Figure BDA00030132162200001332
进而求得空间曲率半径
Figure BDA00030132162200001333
第二段微小圆弧
Figure BDA00030132162200001334
的起始端为Q1。
在静坐标系{B}中，
Figure BDA00030132162200001335
由
Figure BDA00030132162200001336
绕平面Q1R′Q2的单位法向量
Figure BDA00030132162200001337
旋转β′角得到，单位法向量
Figure BDA00030132162200001338
当
Figure BDA00030132162200001339
的长度Δs趋向于0时，
Figure BDA00030132162200001340
即可得到
Figure BDA00030132162200001341
因此所求的第二段矢量
Figure BDA00030132162200001342
为：
Figure BDA00030132162200001343
则第二段微小圆弧
Figure BDA0003013216220000141
的末端Q2在静坐标系{B}上的位置为：
Figure BDA0003013216220000142
根据上述方式，第三段微小圆弧
Figure BDA0003013216220000143
等效的第三段矢量
Figure BDA0003013216220000144
的计算过程与第二段微小圆弧
Figure BDA0003013216220000145
等效的第二段矢量
Figure BDA0003013216220000146
的计算过程相同，在第二段矢量
Figure BDA0003013216220000147
的计算过程中，第三段微小圆弧
Figure BDA0003013216220000148
的空间曲率分量的方向也应跟着
Figure BDA0003013216220000149
绕平面Q1R′Q2的单位法向量
Figure BDA00030132162200001410
旋转β′旋转，即计算动坐标系{A′m}绕平面Q1R′Q2的单位法向量
Figure BDA00030132162200001411
旋转β′角得到的动坐标系{A″m}，得到第三段微小圆弧
Figure BDA00030132162200001412
的目标动坐标系
Figure BDA00030132162200001413
其余计算过程参照第二段微小圆弧
Figure BDA00030132162200001414
等效的第二段矢量
Figure BDA00030132162200001415
的计算过程，在此不再赘述。
第四段微小圆弧
Figure BDA00030132162200001416
等效的第四段矢量
Figure BDA00030132162200001417
至第n段微小圆弧
Figure BDA00030132162200001418
等效的第n段矢量
Figure BDA00030132162200001419
的计算过程参照第二段微小圆弧
Figure BDA00030132162200001420
等效的第二段矢量
Figure BDA00030132162200001421
的计算过程，在此不再赘述。
基于上述过程，得到
Figure BDA00030132162200001422
等效的矢量
Figure BDA00030132162200001423
等效的矢量
Figure BDA00030132162200001424
等效的矢量
Figure BDA00030132162200001425
以及得到Q1、Q2、……、Qn在静坐标系{B}上位置，将Q1、Q2、……、Qn连接起来便得到拟合曲线。
在上述步骤S5中，通过首段矢量
Figure BDA00030132162200001426
和尾段矢量
Figure BDA00030132162200001427
的矢量计算便可获得软管4弯曲角度θ：
Figure BDA00030132162200001428
由公式(12)计算得到目标曲线末端Qn在静坐标系{B}中的空间角度：
Figure BDA00030132162200001429
其中，
Figure BDA00030132162200001430
为静坐标系{B}中Y轴的单位向量，
Figure BDA00030132162200001431
为静坐标系{B}中Z轴的单位向量，
Figure BDA00030132162200001432
为静坐标系{B}中X轴的单位向量，
Figure BDA00030132162200001433
为第n段微小圆弧
Figure BDA00030132162200001434
等效的矢量。
相比于现有技术，本发明的有益效果为：光纤光栅传感器52无需电信号传输电路，无需现场供电，体积小，不受电磁干扰，可靠性高，灵敏度强，寿命长，具有多路传输能力，与光纤传输系统很好地兼容，大量光纤光栅传感器52将数据通过多路传送构成一个光纤传感网络，利用光纤传感网络测量挖泥船耙管的空间角度，使得耙管的工作角度都处于要求的角度，让挖泥船保持在理想的工作状态，不仅提高疏浚工作的效率，同时可以避免施工人员的视觉盲区，保证设备和人员的安全。
以上所述，仅是本发明的较佳实施例而已，并非对本发明做任何形式上的限制，故凡未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化与修饰，均仍属于本发明技术方案的范围内。

Claims (7)

1.一种挖泥船耙管空间角度测量方法，包括耙管，所述耙管包括弯管、中段耙管和耙头，所述弯管通过一软管与中段耙管的一端连接，所述中段耙管的另一端通过另一软管与耙头的一端连接；其特征在于，每个所述软管上沿其周向均匀设有四个传感器装置，所述传感器装置包括记忆合金丝和若干个具有不同中心波长的光纤光栅传感器，所述记忆合金丝设置在软管上，且其长度方向与软管的长度方向一致，所述光纤光栅传感器粘贴于记忆合金丝上并通过光纤串接在一起；
所述方法包括以下步骤：
以所述软管上设有光纤光栅传感器的圆截面为软管测量节点，对于每个软管测量节点，根据第一预设计算规则计算得到软管测量节点的第一曲率标量和第二曲率标量，其中第一曲率标量为软管测量节点处的空间曲率在第一方向上的曲率分量的标量，第二曲率标量为软管测量节点处的空间曲率在第二方向上的曲率分量的标量，第一方向和第二方向互相垂直；
在所述软管上建立静坐标系{B}，静坐标系{B}中的原点O位于所述软管第一端的圆心；
将所述软管等效为在静坐标系{B}中的目标曲线，目标曲线的起始端为静坐标系{B}的原点，并将目标曲线分为多段依次连接的微小圆弧,当微小圆弧的长度取值趋向于0时，微小圆弧等效为该微小圆弧起始端指向其末端的矢量，基于每个软管测量节点的第一曲率标量和第二曲率标量，根据第二预设计算规则计算得到每个微小圆弧等效的矢量；
基于每个微小圆弧等效的矢量，根据预设公式计算得到目标曲线末端在静坐标系{B}中的空间角度，以目标曲线末端在静坐标系{B}中的空间角度作为耙管的空间角度。
2.根据权利要求1所述的挖泥船耙管空间角度测量方法，其特征在于，所述根据第一预设计算规则计算得到第一曲率标量和第二曲率标量的步骤包括：
在软管测量节点上建立平面坐标系O-YZ，其中平面坐标系的原点O为软管测量节点的圆心，Y轴和Z轴分别穿过位于软管测量节点上的两相对光纤光栅传感器，第一方向为Y轴的正方向，第二方向为Z轴的正方向；
通过公式(1)计算第一曲率标量，通过公式(2)计算第二曲率标量：
Figure FDA0003013216210000021
Figure FDA0003013216210000022
其中，k1为第一曲率标量，k2为第二曲率标量，Δλ1为位于Y轴正方向上的光纤光栅传感器的波长偏移量，Δλ2为位于Z轴负方向上的光纤光栅传感器的波长偏移量，Δλ3为位于Y轴负方向上的光纤光栅传感器的波长偏移量，Δλ4为位于Z轴正方向上的光纤光栅传感器的波长偏移量，a1和a2均为系数。
3.根据权利要求2所述的挖泥船耙管空间角度测量方法，其特征在于，所述公式(1)和公式(2)的具体推导过程为：
每个光纤光栅传感器的轴向应变引起中心波长漂移量与其轴向应变满足：
Figure FDA0003013216210000023
其中，Δλ为光纤光栅传感器的中心波长漂移量，λ为光纤光栅传感器的初始中心波长，Pe为光纤光栅传感器的有效弹光系数，εz为光纤光栅传感器的轴向应变；
根据材料力学，对于圆截面弹性梁发生纯弯曲时，其轴向应变εz与曲率之间的关系有：
Figure FDA0003013216210000024
其中，r为光纤光栅传感器中心到平面坐标系O-YZ原点的距离，R为测量点的曲率半径，k为测量点的曲率；
当软管弯曲时，每个光纤光栅传感器的轴向应变与曲率之间满足：
Figure FDA0003013216210000025
其中，εz,i为软管测量节点上的第i个光纤光栅传感器的轴向应变，ri为软管测量节点上的第i个光纤光栅传感器到平面坐标系O-YZ原点的距离，θb为软管弯曲方向与Y轴的夹角，θi为软管测量节点上的第i个光纤光栅传感器与Y轴之间的角度，i的取值为i＝1、2、3、4；
定义一个由平面坐标系O-YZ原点指向一个光纤光栅传感器的表观曲率向量
Figure FDA0003013216210000026
则
Figure FDA0003013216210000027
满足：
Figure FDA0003013216210000028
其中，
Figure FDA0003013216210000029
为沿Y轴的单位向量，
Figure FDA00030132162100000210
为沿Z轴的单位向量；
所有表观曲率向量之和为：
Figure FDA0003013216210000031
每个光纤光栅传感器中心到平面坐标系O-YZ原点的距离均为常量r，把θi带入公式(7)得到：
Figure FDA0003013216210000032
则软管曲率为：
Figure FDA0003013216210000033
其中，λB为位于Y轴上的两个光纤光栅传感器的初始中心波长，λA为位于Z轴上的两个光纤光栅传感器的初始中心波长；
软管曲率在Y轴上的分向量为：
Figure FDA0003013216210000034
由公式(10)得到第一曲率标量
Figure FDA0003013216210000035
软管曲率在Z轴上的分向量为：
Figure FDA0003013216210000036
由公式(11)得到第二曲率标量
Figure FDA0003013216210000037
4.根据权利要求1所述的挖泥船耙管空间角度测量方法，其特征在于，静坐标系{B}中X轴垂直于软管第一端的端面，Y轴和Z轴位于软管第一端所在平面上，且两者分别穿过两相对的传感器装置；
所述基于每个软管测量节点的第一曲率标量和第二曲率标量，根据第二预设计算规则计算得到每个微小圆弧的等效矢量的步骤包括：
以静坐标系{B}的原点为原点，并以单位向量
Figure FDA0003013216210000038
为A轴、单位向量
Figure FDA0003013216210000039
为B轴以及单位向量
Figure FDA00030132162100000310
为C轴，建立动坐标系{A}，动坐标系{A}在初始状态下，其A轴、B轴和C轴分别与静坐标系{B}的Y轴、Z轴和X轴重合，动坐标系{A}用齐次矩阵表示为：
Figure FDA00030132162100000311
以目标曲线起始端到末端的顺序，依次选取一段微小圆弧作为目标微小圆弧，目标微小圆弧等效的矢量的计算过程如下：
步骤1、基于目标微小圆弧两侧的软管测量节点的第一曲率标量，由线性插值法计算得到目标微小圆弧的第一圆弧曲率标量k′1，基于目标微小圆弧两侧的软管测量节点的第二曲率标量，由线性插值法得到目标微小圆弧的第二圆弧曲率标量k′2；
步骤2、判断目标微小圆弧是否为首段微小圆弧；
步骤3、若是首段微小圆弧，则执行步骤4至步骤8，若不是首段微小圆弧，则执行步骤9至步骤10；
步骤4、获取初始状态的动坐标系{A}，以初始状态的动坐标系{A}作为目标微小圆弧的目标动坐标系；
步骤5、基于获取的目标动坐标系，根据公式(13)合成目标微小圆弧的空间曲率：
Figure FDA0003013216210000041
其中，
Figure FDA0003013216210000042
为目标动坐标系中A轴的单位向量，其方向指向A轴的正方向，
Figure FDA0003013216210000043
为目标动坐标系中B轴，其方向指向B轴的正方向；
步骤6、根据公式(14)计算得到目标微小圆弧的空间曲率半径：
Figure FDA0003013216210000044
其中，D1为目标微小圆弧的起始端；
步骤7、在静坐标系{B}中，
Figure FDA0003013216210000045
由
Figure FDA0003013216210000046
绕平面D1RD2的单位法向量
Figure FDA0003013216210000047
旋转β角得到，则
Figure FDA0003013216210000048
满足：
Figure FDA0003013216210000049
其中，D2为目标微小圆弧的末端，
Figure FDA00030132162100000410
Δs为圆弧的长度；
步骤8、根据公式(16)计算得到目标微小圆弧等效的矢量：
Figure FDA00030132162100000411
步骤9、获取前一段微小圆弧的目标动坐标系，根据公式(17)计算获得的目标动坐标系绕平面D1RD2的单位法向量
Figure FDA00030132162100000412
旋转β角得到的动坐标系{A′m}，以动坐标系{A′m}作为目标微小圆弧的目标动坐标系：
Figure FDA00030132162100000413
其中，{Am}为前一段微小圆弧的目标动坐标系；
步骤10、执行步骤5至步骤8。
5.根据权利要求4所述的挖泥船耙管空间角度测量方法，其特征在于，所述根据预设公式计算得到目标曲线末端在静坐标系{B}中的空间角度的步骤中，所述预设公式为：
Figure FDA0003013216210000051
6.根据权利要求1所述的挖泥船耙管空间角度测量方法，其特征在于，所述光纤光栅传感器采用环氧树脂粘贴于记忆合金丝上。
7.根据权利要求1所述的挖泥船耙管空间角度测量方法，其特征在于，所述记忆合金丝为镍钛记忆合金丝。
CN202110381675.5A 2021-04-09 2021-04-09 一种挖泥船耙管空间角度测量方法 Pending CN112985310A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110381675.5A CN112985310A (zh) 2021-04-09 2021-04-09 一种挖泥船耙管空间角度测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110381675.5A CN112985310A (zh) 2021-04-09 2021-04-09 一种挖泥船耙管空间角度测量方法

Publications (1)

Publication Number Publication Date
CN112985310A true CN112985310A (zh) 2021-06-18

Family

ID=76339581

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110381675.5A Pending CN112985310A (zh) 2021-04-09 2021-04-09 一种挖泥船耙管空间角度测量方法

Country Status (1)

Country Link
CN (1) CN112985310A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113959359A (zh) * 2021-10-28 2022-01-21 清华大学深圳国际研究生院 一种管道状连续机器人的形状测量装置和方法
CN114279365A (zh) * 2021-12-21 2022-04-05 浙江贝盛光伏股份有限公司 一种光伏组件弯曲测量系统及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102269582A (zh) * 2011-05-05 2011-12-07 华中科技大学 一种空间三维角度测量装置
CN102445172A (zh) * 2011-10-23 2012-05-09 江苏阳明船舶装备制造技术有限公司 一种用于空间物体位置关系的测量方法
CN109172990A (zh) * 2018-10-25 2019-01-11 河南大学 一种医疗软管路径监测及控制系统
CN109238220A (zh) * 2018-07-20 2019-01-18 山西大学 物体空间转动角度的测量方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102269582A (zh) * 2011-05-05 2011-12-07 华中科技大学 一种空间三维角度测量装置
CN102445172A (zh) * 2011-10-23 2012-05-09 江苏阳明船舶装备制造技术有限公司 一种用于空间物体位置关系的测量方法
CN109238220A (zh) * 2018-07-20 2019-01-18 山西大学 物体空间转动角度的测量方法
CN109172990A (zh) * 2018-10-25 2019-01-11 河南大学 一种医疗软管路径监测及控制系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘刚, 柯映林: "管道机器人全程定位理论和方法研究――基于光纤光栅空间曲率传感器", 浙江大学学报(工学版), no. 06, 30 June 2004 (2004-06-30) *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113959359A (zh) * 2021-10-28 2022-01-21 清华大学深圳国际研究生院 一种管道状连续机器人的形状测量装置和方法
CN113959359B (zh) * 2021-10-28 2024-01-26 清华大学深圳国际研究生院 一种管道状连续机器人的形状测量装置和方法
CN114279365A (zh) * 2021-12-21 2022-04-05 浙江贝盛光伏股份有限公司 一种光伏组件弯曲测量系统及方法
CN114279365B (zh) * 2021-12-21 2024-02-02 浙江贝盛光伏股份有限公司 一种光伏组件弯曲测量系统及方法

Similar Documents

Publication Publication Date Title
JP5851630B2 (ja) ３次元位置計測装置
CN101397903B (zh) 一种应用光纤光栅传感器监测套管周向应变的方法
EP1540272B1 (en) Monitoring of stress in steel catenary risers
CN112985310A (zh) 一种挖泥船耙管空间角度测量方法
US4654520A (en) Structural monitoring system using fiber optics
BRPI0720203A2 (pt) Monitoramento de um cabo de energia flexível
EP1869402A1 (en) Method of applying a strain sensor to a cylindrical structure
CA2716145A1 (en) Monitoring system for well casing
CN111879279B (zh) 高堆石坝的心墙变形监测方法
US10620018B2 (en) Method for measuring the displacement profile of buildings and sensor therefor
CN108020167A (zh) 一种基于光纤光栅的固定式测斜仪器
CN214470669U (zh) 一种基于光纤光栅的耙管
CN109764832A (zh) 一种船舶推进轴系校中方法
CN106525301B (zh) 一种基于分布式光纤感测的力和位移测量方法及传感器
CN109991593B (zh) 基于多芯光纤布里渊散射的有缆潜水器定位装置及方法
Zheng et al. Design, sensing principle and testing of a novel fiber optic displacement sensor based on linear macro-bending loss
Chen et al. Shape monitoring method of submarine cable based on fiber Bragg grating
EP0278143B1 (en) Structural monitoring system using fiber optics
JP2005534033A (ja) 光学的距離および角度の測定方法およびシステム
CN113960328A (zh) 感测装置及其感测二维流速、二维加速度的方法
Johny et al. Theoretical investigation of positional influence of FBG sensors for structural health monitoring of offshore structures
CN104374496A (zh) 导管架平台杆件应力测量的传感器系统及方法
Paloschi et al. Transformation matrices for 3D shape sensing with polyimide-coated multicore optical fiber
CN112986119A (zh) 一种管道腐蚀监测系统及方法
EP2583076B1 (en) System for measuring and monitoring the strain on objects stressed by external forces

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