CN111380476B - Beam type structure deformation measuring method and device based on strain measurement data - Google Patents

Beam type structure deformation measuring method and device based on strain measurement data Download PDF

Info

Publication number
CN111380476B
CN111380476B CN201811612819.8A CN201811612819A CN111380476B CN 111380476 B CN111380476 B CN 111380476B CN 201811612819 A CN201811612819 A CN 201811612819A CN 111380476 B CN111380476 B CN 111380476B
Authority
CN
China
Prior art keywords
optical fiber
beam structure
fiber
ith
structural
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811612819.8A
Other languages
Chinese (zh)
Other versions
CN111380476A (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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201811612819.8A priority Critical patent/CN111380476B/en
Publication of CN111380476A publication Critical patent/CN111380476A/en
Application granted granted Critical
Publication of CN111380476B publication Critical patent/CN111380476B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/16Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
    • G01B11/18Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge using photoelastic elements

Landscapes

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

Abstract

The invention relates to a beam type structure deformation measuring method based on strain measurement data, which is characterized by being applicable to an indirect measuring method for beam type structure deformation monitoring, wherein the strain data are acquired by adopting a fiber grating sensor. The measurement method provided by the invention can be used for measuring the large deformation of the beam structure, wherein the large deformation is defined to be more than 15% of the length of the beam structure. The method can be widely applied to the fields of real-time deformation measurement, structural health monitoring and the like of beam type structures such as high-aspect-ratio flexible wings, large bridge structures, large buildings and the like.

Description

Beam type structure deformation measuring method and device based on strain measurement data
Technical Field
The invention relates to a method and a device for measuring beam structure deformation based on strain measurement data.
Background
Aircraft structures require extensive testing during the design process, including ground and flight testing in terms of structural function, structural strength, load testing, and the like. In these tests, the amount of deformation of the wing structure is one of the important test data. The method is used for accurately and efficiently measuring the structural deformation, and is very important for analyzing the performance and the state of the structure and providing correct basis for designers.
In addition, composite materials have found widespread use in aerospace structures in recent years, and wing structures are generally of low mass. Meanwhile, for high-altitude long-endurance unmanned aerial vehicles, large civil aircraft, large transport aircraft and other aircrafts, the wing aspect ratio is large, and the flexibility is large. The aircraft can generate large deformation when being subjected to air action during flight. Therefore, the structural deformation measurement method needs to take into account the effect of large deformation of the structure.
The existing structural deformation measuring means commonly comprise displacement meter measurement, optical measurement and infrared instrument measurement. The displacement meter is used for measuring, the cost is low, the operation is simple, and the measuring process is visual. However, when a displacement meter is used for measurement, if the structural deformation is large (for example, the deformation of the flexible wing can reach more than 15% of half-span length), especially the direction change caused by large deformation can cause inaccurate measurement results; in addition, the displacement gauge measurement is hardly suitable for real-time measurement of structural deformation in an airborne test. The optical measurement system is a method developed in recent years for capturing the three-dimensional position of the structure space by a plurality of cameras and then calculating the deformation of the structure, and the method also has the problem of imaging distortion, and the measurement result can be inaccurate if the deformation of the structure is large. The infrared instrument can measure three-dimensional deformation, but the cost is higher, and the operation is complicated, and the measurement process is not directly perceived, can't in time discover when having the problem, does not utilize real-time deformation measurement of structure when airborne is tested.
Therefore, a three-dimensional configuration measuring method applicable to large structural deformation is needed when an aviation structural test is carried out.
Disclosure of Invention
According to an aspect of the present invention, there is provided a beam structure deformation measurement apparatus based on strain measurement data, characterized by comprising:
first to fourth optical fibers S1, S2, S3, S4, each fiber being connected in series with N fiber grating sensors for strain measurement, wherein the X-axis coordinates of the set point of the ith fiber grating sensor among the N fiber grating sensors on the first to fourth optical fibers are the same,
the first optical fiber S1 and the second optical fiber S2 are respectively disposed on the upper surface of the beam structure,
the third optical fiber S3 is disposed on the upper surface of the beam structure,
the fourth optical fiber S4 is disposed on the lower surface of the beam structure,
the X axis points to the length direction of the beam structure, the Y axis points to the normal direction of the upper surface of the beam structure, the Z axis is defined by the right-hand rule,
fiber bragg grating sensor(1)ε1(1)ε2(1)ε3(1)ε4(1)ε5Connect in series on the first fiber S1, the fiber grating sensor(2)ε1(2)ε2(2)ε3(2)ε4(2)ε5The fiber bragg grating sensors are connected in series on a second optical fiber S2, the arrangement directions of the fiber bragg grating sensors are all directions parallel to the X axis,
fiber bragg grating sensor(3)ε1(3)ε2(3)ε3(3)ε4(3)ε5The third optical fiber S3 is connected in series, the arrangement direction of the sensors alternately forms an angle of +45 degrees with the X axis,
fiber bragg grating sensor(4)ε1(4)ε2(4)ε3(4)ε4(4)ε5The fourth optical fiber S4 are connected in series, the arrangement directions of the sensors are all parallel to the Z axis,
wherein the content of the first and second substances,
with E1,E2,E3Three basis vectors, E, respectively, representing a coordinate system1Along the length direction of the beam;
in g1(s),g2(s),g3(s) three basis vectors representing the local coordinate system at length s before beam deformation, respectively;
with G1(s),G2(s),G3(s) represent three basis vectors of the local coordinate system at length s after the beam is deformed, and have:
G1(s)=[G11(s) G12(s) G13(s)]T
G2(s)=[G21(s) G22(s) G23(s)]T
G3(s)=[G31(s) G32(s) G33(s)]T
representing a space position vector of the central reference axis s in the reference coordinate system before the beam is deformed by r(s);
and R(s) represents the space position vector of the center reference axis s in the reference coordinate system after the beam is deformed, and the space position vector comprises:
R(si)=[Ri1 Ri2 Ri3]T. Column vector in inventive content step five
Figure BDA0001925172680000023
Can be defined according to the position vector and the base vector of the coordinate system:
Figure BDA0001925172680000021
the beam structure deformation measuring device is adapted to perform the following operations:
A) the arrangement position of the fiber grating sensor on each optical fiber is sequentially represented as:
Figure BDA0001925172680000022
wherein the content of the first and second substances,
Figure BDA0001925172680000031
the distance between the first optical fiber S3 and the fourth optical fiber S4 and the center line of the beam type structure at the ith measuring point is shown, wherein the ith measuring point is the position of the ith fiber grating sensor on each optical fiber;
Figure BDA0001925172680000032
a coordinate component in the Y direction relative to the structural beam structure centerline at the ith measurement point of the first optical fiber S1;
Figure BDA0001925172680000033
represents the coordinate component of the first optical fiber S1 in the Z direction relative to the structural beam structure centerline at the ith test point;
Figure BDA0001925172680000034
represents the coordinate component of the second optical fiber S2 in the Y direction relative to the structural beam structure centerline at the ith measuring point;
Figure BDA0001925172680000035
represents the Z-direction coordinate component of the second optical fiber S2 relative to the structural beam structure centerline at the ith test point,
B) adopting optical fiber measurement technology to measure structural strain at the position of each group of 4 sensors with the same X coordinate(1)εi}i=1,2,…N,{(2)εi}i=1,2,…N,{(3)εi}i=1,2,…N,{(4)εi}i=1,2,…NWherein: {(1)εi}i=1,2,…NThe structural strain measured at the ith test point on the first optical fiber S1; {(2)εi}i=1,2,…NThe structural strain measured at the ith test point on the second optical fiber S2; {(3)εi}i=1,2,…NThe structural strain measured at the ith test point on the third optical fiber S3; {(4)εi}i=1,2,…NThe structural strain measured at the ith test point on the fourth fiber S4,
C) calculating the three-dimensional curvature of the beam structure center line at each ith measurement point according to the measured structural strain and the following relation
Figure BDA00019251726800000311
Figure BDA0001925172680000036
Figure BDA0001925172680000037
D) The structural deformation-curvature relationship is calculated according to the following system of equations:
Figure BDA0001925172680000038
Figure BDA0001925172680000039
Figure BDA00019251726800000310
Figure BDA0001925172680000041
wherein
Figure BDA0001925172680000042
E) Calculating the structural deformation at the measuring point according to the relation in the step (C) and the step (D):
Figure BDA0001925172680000043
Figure BDA0001925172680000044
is a column vector of 12X1, the 1 st, 5 th and 9 th elements of which respectively correspond to the three-dimensional positions of the structure at the measuring point i in the space,
Figure BDA0001925172680000045
representing the position vector of the structure at the end point position, thereby representing the three-dimensional position of the ith measuring point after the beam structure is deformed as
Figure BDA0001925172680000046
According to another aspect of the present invention, there is provided a beam structure deformation measurement method based on strain measurement data, comprising:
the first optical fiber S1 and the second optical fiber S2 are respectively disposed on the upper surface of the beam structure,
the third optical fiber S3 is disposed on the upper surface of the beam structure,
the fourth optical fiber S4 is disposed on the lower surface of the beam structure,
each optical fiber is connected with N fiber grating sensors in series for strain measurement, wherein the X-axis coordinates of the set points of the ith fiber grating sensor in the N fiber grating sensors on the first to fourth optical fibers are the same, the X-axis points to the length direction of the beam structure, the Y-axis points to the normal direction of the upper surface of the beam structure, and the Z-axis is defined by the right-hand rule, which specifically includes:
opening and closing optical fiber grating sensor(1)ε1(1)ε2(1)ε3(1)ε4(1)ε5Connect in series with the first fiber S1 to make the fiber grating sensor(2)ε1(2)ε2(2)ε3(2)ε4(2)ε5The fiber bragg grating sensors are connected in series on a second optical fiber S2, wherein the arrangement directions of the fiber bragg grating sensors are all directions parallel to the X axis,
opening and closing optical fiber grating sensor(3)ε1(3)ε2(3)ε3(3)ε4(3)ε5Are connected in series to a third optical fibre S3, wherein the sensors are arranged alternately at +45 deg. to the X axis,
opening and closing optical fiber grating sensor(4)ε1(4)ε2(4)ε3(4)ε4(4)ε5Are connected in series to a fourth optical fibre S4, wherein the sensors are all arranged in a direction parallel to the Z axis,
with E1,E2,E3Three basis vectors, E, respectively, representing a coordinate system1Along the length direction of the beam;
in g1(s),g2(s),g3(s) three basis vectors representing the local coordinate system at length s before beam deformation, respectively;
with G1(s),G2(s),G3(s) represent three basis vectors of the local coordinate system at length s after the beam is deformed, and have:
G1(s)=[G11(s) G12(s) G13(s)]T
G2(s)=[G21(s) G22(s) G23(s)]T
G3(s)=[G31(s) G32(s) G33(s)]T
representing a space position vector of the central reference axis s in the reference coordinate system before the beam is deformed by r(s);
and R(s) represents the space position vector of the center reference axis s in the reference coordinate system after the beam is deformed, and the space position vector comprises:
R(si)=[Ri1 Ri2 Ri3]T
the beam structure deformation measurement method further includes:
A) the arrangement position of the fiber grating sensor on each optical fiber is sequentially represented as:
Figure BDA0001925172680000051
wherein the content of the first and second substances,
Figure BDA0001925172680000052
the distance between the first optical fiber S3 and the fourth optical fiber S4 and the center line of the beam type structure at the ith measuring point is shown, wherein the ith measuring point is the position of the ith fiber grating sensor on each optical fiber;
Figure BDA0001925172680000057
a coordinate component in the Y direction relative to the structural beam structure centerline at the ith measurement point of the first optical fiber S1;
Figure BDA0001925172680000058
represents the coordinate component of the first optical fiber S1 in the Z direction relative to the structural beam structure centerline at the ith test point;
Figure BDA0001925172680000059
represents the coordinate component of the second optical fiber S2 in the Y direction relative to the structural beam structure centerline at the ith measuring point;
Figure BDA00019251726800000510
represents the Z-direction coordinate component of the second optical fiber S2 relative to the structural beam structure centerline at the ith test point,
B) adopting optical fiber measurement technology to measure X-rayStructural strain retaining tone at the same position of each group of 4 sensors(1)εi}i=1,2,…N,{(2)εi}i=1,2,…N,{(3)εi}i=1,2,…N,{(4)εi}i=1,2,…NWherein: {(1)εi}i=1,2,…NThe structural strain measured at the ith test point on the first optical fiber S1; {(2)εi}i=1,2,…NThe structural strain measured at the ith test point on the second optical fiber S2; {(3)εi}i=1,2,…NThe structural strain measured at the ith test point on the third optical fiber S3; {(4)εi}i=1,2,…NThe structural strain measured at the ith test point on the fourth fiber S4,
C) calculating the three-dimensional curvature of the beam structure center line at each ith measurement point according to the measured structural strain and the following relation
Figure BDA0001925172680000053
Figure BDA0001925172680000054
Figure BDA0001925172680000055
D) The structural deformation-curvature relationship is calculated according to the following system of equations:
Figure BDA0001925172680000056
Figure BDA0001925172680000061
Figure BDA0001925172680000062
Figure BDA0001925172680000063
wherein
Figure BDA0001925172680000064
E) Calculating the structural deformation at the measuring point according to the relation in the step (C) and the step (D):
Figure BDA0001925172680000065
Figure BDA0001925172680000066
is a column vector of 12X1, the 1 st, 5 th and 9 th elements of which respectively correspond to the three-dimensional positions of the structure at the measuring point i in the space,
wherein the column vector
Figure BDA0001925172680000067
Is determined as:
Figure BDA0001925172680000068
Figure BDA0001925172680000069
representing the position vector of the structure at the end point position, thereby representing the three-dimensional position of the ith measuring point after the beam structure is deformed as
Figure BDA00019251726800000610
Drawings
Fig. 1 is a schematic diagram of a strain sensor arrangement of a beam structure deformation measurement device based on strain measurement data according to one embodiment of the present invention.
Fig. 2 is a schematic view of a deformation configuration of a beam structure suitable for applying the beam structure deformation measurement apparatus based on strain measurement data according to the present invention.
Detailed Description
In order to overcome the defects of the prior art and achieve the aim, the invention provides a structural deformation measuring method based on strain measurement data, which is suitable for a beam type structure such as a wing, a tail wing and the like. In the method, strain measurement data is acquired by a fiber grating sensor. The invention comprises the following steps:
step one, selecting N measuring points on a test piece along the length direction of a beam, respectively pasting 1 fiber grating sensor on each optical fiber at each measuring point, and determining the pasting position of the fiber grating sensors:
Figure BDA0001925172680000071
wherein the content of the first and second substances,
Figure BDA0001925172680000072
represents the distance from the beam structure centerline at the ith test point for fiber number S3 and fiber number S4;
Figure BDA0001925172680000073
a coordinate component in the Y direction of the optical fiber denoted by the number S1 with respect to the center line of the structural beam structure at the i-th measurement point;
Figure BDA0001925172680000074
a coordinate component in the Z direction from the center line of the structural beam structure at the i-th measuring point of the optical fiber with the number S1;
Figure BDA0001925172680000075
a coordinate component in the Y direction of the optical fiber denoted by the number S2 with respect to the center line of the structural beam structure at the i-th measurement point;
Figure BDA0001925172680000076
the optical fiber denoted by the number S2 is in the opposite structural beam structure at the ith test pointThe coordinate component of the center line in the Z direction.
Step two, adopting an optical fiber measurement technology to measure the structural strain of the positions of the 4 sensors at each measuring point: {(1)εi}i=1,2,…N,{(2)εi}i=1,2,…N,{(3)εi}i=1,2,…N,{(4)εi}i=1,2,…N{(1)εi}i=1,2,…NRepresents the structural strain measured at the ith test point in the optical fiber numbered S1; {(2)εi}i=1,2,…NRepresents the structural strain measured at the ith test point in the optical fiber numbered S2; {(3)εi}i=1,2,…NRepresents the structural strain measured at the ith test point in the optical fiber numbered S3; {(4)εi}i=1,2,…NThe structural strain measured at the ith station in the fiber numbered S4 is shown.
Step three, calculating the three-dimensional curvature of the central line of the beam structure at each measuring point i according to the measured strain and the following relational expression
Figure BDA00019251726800000710
Figure BDA0001925172680000077
Figure BDA0001925172680000078
Step four, calculating a structural deformation-curvature relation according to the following equation:
Figure BDA0001925172680000079
Figure BDA0001925172680000081
Figure BDA0001925172680000082
Figure BDA0001925172680000083
wherein
Figure BDA0001925172680000084
Step five, calculating the structural deformation at the measuring point according to the relational expression in the step three and the step four:
Figure BDA0001925172680000085
among them, in the above-mentioned case,
Figure BDA0001925172680000086
is a 12X1 column vector defined as follows:
Figure BDA0001925172680000087
its 1 st, 5 th and 9 th elements respectively correspond to the three-dimensional deformation of the structure at the measuring point i in the space. Wherein the definition of each element is shown in figure 2.
Figure BDA0001925172680000088
Shown is the position vector of the structure at the end position (the left-most root shown in fig. 1).
The method is not only suitable for the condition of beam root solid support, but also can directly measure the beam structure deformation under different constraint conditions only by determining boundary conditions.
The invention is further described below with reference to the accompanying drawings.
Fig. 1 is a schematic diagram of a strain sensor arrangement. The strain sensor is a fiber grating sensor, and N strain sensors (N is 5 in figure 1)) The fiber grating sensor is connected in series in an optical fiber, and specifically, the sensor(1)ε1(1)ε2(1)ε3(1)ε4(1)ε5Connected in series in the optical fiber S1; transducer opening(2)ε1(2)ε2(2)ε3(2)ε4(2)ε5Connected in series in the optical fiber S2; transducer opening(3)ε1(3)ε2(3)ε3(3)ε4(3)ε5Connected in series in the optical fiber S3; transducer opening(4)ε1(4)ε2(4)ε3(4)ε4(4)ε5Is connected in series in the optical fiber S3. The sensor has the advantages of light weight, electromagnetic interference resistance and the like, and is convenient for airborne arrangement. The coordinate system shown in fig. 1 is defined as follows: the origin is located at the center of the cross section of the root of the beam, the X axis points to the length direction of the beam, the Y axis points to the out-of-plane bending direction, and the Z axis is defined by the right-hand rule. And the coordinate positions of the sensors at the ith measuring point in the four optical fibers in the X-axis direction are the same. The optical fibers S1, S2 are arranged on the upper surface of the beam, respectively, in a direction parallel to the X axis. The optical fibers S3 are arranged on the upper surface of the beam, and the directions of the sensors at all measuring points on the S3 and the X axis are alternately arranged at an angle of +/-45 degrees; the optical fiber S4 is disposed on the lower surface of the beam with the sensor direction parallel to the Z-axis.
FIG. 2 is a schematic view of a beam structure deformation configuration, wherein E1,E2,E3Three base vectors of a reference coordinate system can be selected at will, and the origin of the coordinate system is the center of the cross section of the root part of the beam, E1Along the length direction of the beam; g1(s),g2(s),g3(s) three basis vectors of a local coordinate system at the length s before the beam is deformed respectively; g1(s),G2(s),G3(s) three basis vectors of the local coordinate system at the length s after the beam is deformed, and having G1(s)=[G11(s) G12(s) G13(s)]T,G2(s)=[G21(s) G22(s) G23(s)]T,G3(s)=[G31(s) G32(s) G33(s)]T(ii) a r(s) is a spatial position vector of the central reference axis s in the reference coordinate system before the beam is deformed; r(s) is a space position vector of the center reference axis s in the reference coordinate system after the beam is deformed, and R(s)i)=[Ri1 Ri2 Ri3]T. Column vector in inventive content step five
Figure BDA0001925172680000091
Can be defined according to the position vector and the base vector of the coordinate system:
Figure BDA0001925172680000092
Figure BDA0001925172680000094
shown is the position vector of the structure at the end position (the left-most root shown in fig. 1).
Step one, selecting N measuring points on a test piece along the length direction of a beam, respectively pasting 1 fiber grating sensor on each optical fiber at each measuring point, and determining the pasting position of the fiber grating sensors:
Figure BDA00019251726800000911
wherein the content of the first and second substances,
Figure BDA0001925172680000095
represents the distance from the beam structure centerline at the ith test point for fiber number S3 and fiber number S4;
Figure BDA0001925172680000096
a coordinate component in the Y direction of the optical fiber denoted by the number S1 with respect to the center line of the structural beam structure at the i-th measurement point;
Figure BDA0001925172680000097
a coordinate component in the Z direction from the center line of the structural beam structure at the i-th measuring point of the optical fiber with the number S1;
Figure BDA0001925172680000098
a coordinate component in the Y direction of the optical fiber denoted by the number S2 with respect to the center line of the structural beam structure at the i-th measurement point;
Figure BDA0001925172680000099
the coordinate component in the Z direction of the optical fiber denoted by number S2 with respect to the center line of the structural beam structure at the i-th measurement point.
Step two, adopting an optical fiber measurement technology to measure the structural strain of the positions of the 4 sensors at each measuring point: {(1)εi}i=1,2,…N,{(2)εi}i=1,2,…N,{(3)εi}i=1,2,…N,{(4)εi}i=1,2,…N{(1)εi}i=1,2,…NRepresents the structural strain measured at the ith test point in the optical fiber numbered S1; {(2)εi}i=1,2,…NRepresents the structural strain measured at the ith test point in the optical fiber numbered S2; {(3)εi}i=1,2,…NRepresents the structural strain measured at the ith test point in the optical fiber numbered S3; {(4)εi}i=1,2,…NThe structural strain measured at the ith station in the fiber numbered S4 is shown.
Step three, calculating the three-dimensional curvature of the central line of the beam structure at each measuring point i according to the measured strain and the following relational expression
Figure BDA0001925172680000101
Figure BDA0001925172680000102
Figure BDA0001925172680000103
Step four, calculating a structural deformation-curvature relation according to the following equation:
Figure BDA0001925172680000104
Figure BDA0001925172680000105
Figure BDA0001925172680000106
Figure BDA0001925172680000107
wherein
Figure BDA0001925172680000108
Step five, calculating the structural deformation at the measuring point according to the relational expression in the step three and the step four:
Figure BDA0001925172680000109
among them, in the above-mentioned case,
Figure BDA00019251726800001010
is a column vector of 12X1, whose elements 1,5 and 9 correspond to the three-dimensional positions of the structure at the measurement point i in space, respectively.
Figure BDA00019251726800001011
Shown is the position vector of the structure at the end position (the left-most root shown in fig. 1).
Thus, the three-dimensional position of the beam structure at the measuring point i after deformation is
Figure BDA0001925172680000111
After the three-dimensional positions of the measuring points are known, a continuous beam structure deformation curve can be obtained through a spline interpolation method.

Claims (4)

1. A beam structure deformation measuring device based on strain measurement data is characterized by comprising:
first to fourth optical fibers S1, S2, S3, S4, each fiber being connected in series with N fiber grating sensors for strain measurement, wherein the X-axis coordinates of the set point of the ith fiber grating sensor among the N fiber grating sensors on the first to fourth optical fibers are the same,
the first optical fiber S1 and the second optical fiber S2 are respectively disposed on the upper surface of the beam structure,
the third optical fiber S3 is disposed on the upper surface of the beam structure,
the fourth optical fiber S4 is disposed on the lower surface of the beam structure,
the X axis points to the length direction of the beam structure, the Y axis points to the normal direction of the upper surface of the beam structure, the Z axis is defined by the right-hand rule,
fiber bragg grating sensor(1)ε1(1)ε2(1)ε3(1)ε4(1)ε5Connect in series on the first fiber S1, the fiber grating sensor(2)ε1(2)ε2(2)ε3(2)ε4(2)ε5The fiber bragg grating sensors are connected in series on a second optical fiber S2, the arrangement directions of the fiber bragg grating sensors are all directions parallel to the X axis,
fiber bragg grating sensor(3)ε1(3)ε2(3)ε3(3)ε4(3)ε5The third optical fiber S3 is connected in series, the arrangement direction of the sensors alternately forms an angle of +45 degrees with the X axis,
fiber bragg grating sensor(4)ε1(4)ε2(4)ε3(4)ε4(4)ε5The fourth optical fiber S4 are connected in series, the arrangement directions of the sensors are all parallel to the Z axis,
wherein the content of the first and second substances,
with E1,E2,E3Three basis vectors, E, respectively, representing a coordinate system1Along the length direction of the beam;
in g1(s),g2(s),g3(s) three basis vectors representing the local coordinate system at length s before beam deformation, respectively;
with G1(s),G2(s),G3(s) represent three basis vectors of the local coordinate system at length s after the beam is deformed, and have:
G1(s)=[G11(s) G12(s) G13(s)]T
G2(s)=[G21(s) G22(s) G23(s)]T
G3(s)=[G31(s) G32(s) G33(s)]T
representing a space position vector of the central reference axis s in the reference coordinate system before the beam is deformed by r(s);
and R(s) represents the space position vector of the center reference axis s in the reference coordinate system after the beam is deformed, and the space position vector comprises:
R(si)=[Ri1 Ri2 Ri3]T
column vector
Figure FDA0002824012280000021
According to the spatial position vector and the base vector of the coordinate system, defining:
Figure FDA00028240122800000212
the beam structure deformation measuring device is adapted to perform the following operations:
A) the arrangement position of the fiber grating sensor on each optical fiber is sequentially represented as:
Figure FDA0002824012280000022
wherein the content of the first and second substances,
Figure FDA0002824012280000023
the distance between the first optical fiber S3 and the fourth optical fiber S4 and the center line of the beam type structure at the ith measuring point is shown, wherein the ith measuring point is the position of the ith fiber grating sensor on each optical fiber;
Figure FDA0002824012280000024
a coordinate component in the Y direction relative to the structural beam structure centerline at the ith measurement point of the first optical fiber S1;
Figure FDA0002824012280000025
represents the coordinate component of the first optical fiber S1 in the Z direction relative to the structural beam structure centerline at the ith test point;
Figure FDA0002824012280000026
represents the coordinate component of the second optical fiber S2 in the Y direction relative to the structural beam structure centerline at the ith measuring point;
Figure FDA0002824012280000027
represents the Z-direction coordinate component of the second optical fiber S2 relative to the structural beam structure centerline at the ith test point,
B) adopting optical fiber measurement technology to measure structural strain at the position of each group of 4 sensors with the same X coordinate(1)εi}i=1,2,…N,{(2)εi}i=1,2,…N,{(3)εi}i=1,2,…N,{(4)εi}i=1,2,…NWherein: {(1)εi}i=1,2,…NThe measured structure at the i-th measuring point on the first optical fiber S1Changing; {(2)εi}i=1,2,…NThe structural strain measured at the ith test point on the second optical fiber S2; {(3)εi}i=1,2,…NThe structural strain measured at the ith test point on the third optical fiber S3; {(4)εi}i=1,2,…NThe structural strain measured at the ith test point on the fourth fiber S4,
C) calculating the three-dimensional curvature of the beam structure center line at each ith measurement point according to the measured structural strain and the following relation
Figure FDA0002824012280000028
Figure FDA0002824012280000029
Figure FDA00028240122800000210
D) The structural deformation-curvature relationship is calculated according to the following system of equations:
Figure FDA00028240122800000211
Figure FDA0002824012280000031
Figure FDA0002824012280000032
Figure FDA0002824012280000033
wherein
Figure FDA0002824012280000034
E) Calculating the structural deformation at the measuring point according to the relation in the step (C) and the step (D):
Figure FDA0002824012280000035
Figure FDA0002824012280000036
is a column vector of 12X1, the 1 st, 5 th and 9 th elements of which respectively correspond to the three-dimensional positions of the structure at the measuring point i in the space,
Figure FDA0002824012280000037
representing a position vector of the structure at the location of the end point,
thus, the three-dimensional position of the i-th measurement point after deformation of the beam structure is represented as
Figure FDA0002824012280000038
2. The beam structure deformation measuring device according to claim 1, wherein:
after determining the three-dimensional position of each measuring point
Figure FDA0002824012280000039
And then, obtaining a continuous beam structure deformation curve by a spline interpolation method.
3. A beam structure deformation measurement method based on strain measurement data is characterized by comprising the following steps:
the first optical fiber S1 and the second optical fiber S2 are respectively disposed on the upper surface of the beam structure,
the third optical fiber S3 is disposed on the upper surface of the beam structure,
the fourth optical fiber S4 is disposed on the lower surface of the beam structure,
each optical fiber is connected with N fiber grating sensors in series for strain measurement, wherein the X-axis coordinates of the set points of the ith fiber grating sensor in the N fiber grating sensors on the first to fourth optical fibers are the same, the X-axis points to the length direction of the beam structure, the Y-axis points to the normal direction of the upper surface of the beam structure, and the Z-axis is defined by the right-hand rule, which specifically includes:
opening and closing optical fiber grating sensor(1)ε1(1)ε2(1)ε3(1)ε4(1)ε5Connect in series with the first fiber S1 to make the fiber grating sensor(2)ε1(2)ε2(2)ε3(2)ε4(2)ε5The fiber bragg grating sensors are connected in series on a second optical fiber S2, wherein the arrangement directions of the fiber bragg grating sensors are all directions parallel to the X axis,
opening and closing optical fiber grating sensor(3)ε1(3)ε2(3)ε3(3)ε4(3)ε5Are connected in series to a third optical fibre S3, wherein the sensors are arranged alternately at +45 deg. to the X axis,
opening and closing optical fiber grating sensor(4)ε1(4)ε2(4)ε3(4)ε4(4)ε5Are connected in series to a fourth optical fibre S4, wherein the sensors are all arranged in a direction parallel to the Z axis,
with E1,E2,E3Three basis vectors, E, respectively, representing a coordinate system1Along the length direction of the beam;
in g1(s),g2(s),g3(s) three basis vectors representing the local coordinate system at length s before beam deformation, respectively;
with G1(s),G2(s),G3(s) represent three basis vectors of the local coordinate system at length s after the beam is deformed, and have:
G1(s)=[G11(s) G12(s) G13(s)]T
G2(s)=[G21(s) G22(s) G23(s)]T
G3(s)=[G31(s) G32(s) G33(s)]T
representing a space position vector of the central reference axis s in the reference coordinate system before the beam is deformed by r(s);
and R(s) represents the space position vector of the center reference axis s in the reference coordinate system after the beam is deformed, and the space position vector comprises:
R(si)=[Ri1 Ri2 Ri3]T
the beam structure deformation measurement method further includes:
A) the arrangement position of the fiber grating sensor on each optical fiber is sequentially represented as:
Figure FDA0002824012280000041
wherein the content of the first and second substances,
Figure FDA0002824012280000042
the distance between the first optical fiber S3 and the fourth optical fiber S4 and the center line of the beam type structure at the ith measuring point is shown, wherein the ith measuring point is the position of the ith fiber grating sensor on each optical fiber;
Figure FDA0002824012280000043
a coordinate component in the Y direction relative to the structural beam structure centerline at the ith measurement point of the first optical fiber S1;
Figure FDA0002824012280000044
represents the coordinate component of the first optical fiber S1 in the Z direction relative to the structural beam structure centerline at the ith test point;
Figure FDA0002824012280000045
represents the coordinate component of the second optical fiber S2 in the Y direction relative to the structural beam structure centerline at the ith measuring point;
Figure FDA0002824012280000046
represents the Z-direction coordinate component of the second optical fiber S2 relative to the structural beam structure centerline at the ith test point,
B) adopting optical fiber measurement technology to measure structural strain at the position of each group of 4 sensors with the same X coordinate(1)εi}i=1,2,…N,{(2)εi}i=1,2,…N,{(3)εi}i=1,2,…N,{(4)εi}i=1,2,…NWherein: {(1)εi}i=1,2,…NThe structural strain measured at the ith test point on the first optical fiber S1; {(2)εi}i=1,2,…NThe structural strain measured at the ith test point on the second optical fiber S2; {(3)εi}i=1,2,…NThe structural strain measured at the ith test point on the third optical fiber S3; {(4)εi}i=1,2,…NThe structural strain measured at the ith test point on the fourth fiber S4,
C) calculating the three-dimensional curvature of the beam structure center line at each ith measurement point according to the measured structural strain and the following relation
Figure FDA0002824012280000051
Figure FDA0002824012280000052
Figure FDA0002824012280000053
D) The structural deformation-curvature relationship is calculated according to the following system of equations:
Figure FDA0002824012280000054
Figure FDA0002824012280000055
Figure FDA0002824012280000056
Figure FDA0002824012280000057
wherein
Figure FDA0002824012280000058
E) Calculating the structural deformation at the measuring point according to the relation in the step (C) and the step (D):
Figure FDA0002824012280000059
Figure FDA00028240122800000510
is a column vector of 12X1, the 1 st, 5 th and 9 th elements of which respectively correspond to the three-dimensional positions of the structure at the measuring point i in the space,
wherein the column vector
Figure FDA00028240122800000511
Is determined as:
Figure FDA00028240122800000512
Figure FDA0002824012280000061
representing a position vector of the structure at the location of the end point,
thus, the three-dimensional position of the i-th measurement point after deformation of the beam structure is represented as
Figure FDA0002824012280000062
4. A beam structure deformation measurement method according to claim 3, characterized by further comprising:
after determining the three-dimensional position of each measuring point
Figure FDA0002824012280000063
And then, obtaining a continuous beam structure deformation curve by a spline interpolation method.
CN201811612819.8A 2018-12-27 2018-12-27 Beam type structure deformation measuring method and device based on strain measurement data Active CN111380476B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811612819.8A CN111380476B (en) 2018-12-27 2018-12-27 Beam type structure deformation measuring method and device based on strain measurement data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811612819.8A CN111380476B (en) 2018-12-27 2018-12-27 Beam type structure deformation measuring method and device based on strain measurement data

Publications (2)

Publication Number Publication Date
CN111380476A CN111380476A (en) 2020-07-07
CN111380476B true CN111380476B (en) 2021-03-09

Family

ID=71216540

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811612819.8A Active CN111380476B (en) 2018-12-27 2018-12-27 Beam type structure deformation measuring method and device based on strain measurement data

Country Status (1)

Country Link
CN (1) CN111380476B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113916123B (en) * 2021-11-22 2022-05-27 北京航空航天大学 Strain-based large flexible beam structure three-dimensional space configuration determination method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61104237A (en) * 1984-10-26 1986-05-22 Nippon Telegr & Teleph Corp <Ntt> Apparatus and method for measuring strain of optical fiber due to bending of optical cable
JP2001050781A (en) * 1999-08-11 2001-02-23 Nippon Telegr & Teleph Corp <Ntt> Disturbed state measuring device for ground
CN105783758A (en) * 2016-04-07 2016-07-20 华北电力大学(保定) Self-heterodyne single-ended vector BOTDA-based dynamic strain measurement method and apparatus

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103123254B (en) * 2013-02-07 2015-05-27 南开大学 Two-dimensional bending vector sensor based on fiber grating interference structure
CN104794284B (en) * 2015-04-22 2017-08-22 西安电子科技大学 A kind of smart skins antenna electric compensation method based on embedded fiber grating

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS61104237A (en) * 1984-10-26 1986-05-22 Nippon Telegr & Teleph Corp <Ntt> Apparatus and method for measuring strain of optical fiber due to bending of optical cable
JP2001050781A (en) * 1999-08-11 2001-02-23 Nippon Telegr & Teleph Corp <Ntt> Disturbed state measuring device for ground
CN105783758A (en) * 2016-04-07 2016-07-20 华北电力大学(保定) Self-heterodyne single-ended vector BOTDA-based dynamic strain measurement method and apparatus

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Strain-based Shape Prediction for Flexible;Yang Meng1;《AIAA SciTech Forum》;20190131;全文 *

Also Published As

Publication number Publication date
CN111380476A (en) 2020-07-07

Similar Documents

Publication Publication Date Title
JP2913005B2 (en) Flight velocity vector detection system using a truncated polygonal pitot tube probe and a truncated polygonal pitot tube probe
CN108413887B (en) Wing-shaped deformation measuring method, device and platform of fiber bragg grating assisted distributed POS
CN109977448A (en) A kind of rotor aerodynamics analysis method based on actual measurement structural loads
BR102017020272B1 (en) AIRCRAFT MONITORING SYSTEM AND METHOD FOR MONITORING THE MOVEMENT OF AN AIRCRAFT STRUCTURE
CN110006563B (en) Distributed optical fiber decoupling measurement method for blade flapping and shimmy bending moment of helicopter
CN107869960B (en) A kind of aircraft wing deformation measuring device, installation method and test method
Damljanović et al. T-38 wind-tunnel data quality assurance based on testing of a standard model
CN111043979A (en) Plate-shaped structure deformation detection device and method based on fiber bragg grating sensor array
CN104597907A (en) Method for accurately evaluating flight of UAV (unmanned aerial vehicle) inspection system of overhead transmission line
CN107917717A (en) Flight parameter measurement device and corresponding measuring method with optical skew sensor
Bakunowicz et al. Measuring structure deformations of a composite glider by optical means with on-ground and in-flight testing
CN111380476B (en) Beam type structure deformation measuring method and device based on strain measurement data
Abrego et al. Blade displacement measurement technique applied to a full-scale rotor test
Raab et al. Dynamic flight load measurements with MEMS pressure sensors
Raab et al. In-flight testing of MEMS pressure sensors for flight loads determination
CN106323587B (en) The monocular video high-precision measuring method of wing model in wind tunnel flexible deformation
Ocokoljić et al. Testing of a standard model in the VTI's large-subsonic wind-tunnel facility to establish users' confidence
CN110849358B (en) Measuring device, measuring method and mounting method for array antenna phase center
Preisighe Viana Time-domain system identification of rigid-body multipoint loads model
Lau et al. Ground Vibration Testing Master Class: modern testing and analysis concepts applied to an F-16 aircraft
Sturm et al. Morphing shell design of a sailplane with a morphing forward wing section
Luczak et al. Static and dynamic testing of the full scale helicopter rotor blades
CN115169162B (en) Method and device for predicting airplane vibration environment and computer readable storage medium
CN112484751B (en) Method for measuring position and attitude of spacecraft verifier in relatively large space test field coordinate system
Heaney et al. Distributed sensing of a cantilever beam and plate using a fiber optic sensing system

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