CN111380476A - 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 PDFInfo
- Publication number
- CN111380476A CN111380476A CN201811612819.8A CN201811612819A CN111380476A CN 111380476 A CN111380476 A CN 111380476A CN 201811612819 A CN201811612819 A CN 201811612819A CN 111380476 A CN111380476 A CN 111380476A
- 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.)
- Granted
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/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/18—Measuring 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
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)=[Ri1Ri2Ri3]T. Column vector in inventive content step fiveCan be defined according to the position vector and the base vector of the coordinate system:
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:
wherein the content of the first and second substances,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;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;indicating that the first optical fiber S1 opposes the structural beam junction at the ith test pointConstructing a coordinate component of the central line in the Z direction;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;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
D) The structural deformation-curvature relationship is calculated according to the following system of equations:
E) Calculating the structural deformation at the measuring point according to the relation in the step (C) and the step (D):
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,
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
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)=[Ri1Ri2Ri3]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:
wherein the content of the first and second substances,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;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;indicating that the first optical fiber S1 is at the i-th measurement point in the Z-direction relative to the structural beam structure centerlineThe coordinate component of (a);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;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
D) The structural deformation-curvature relationship is calculated according to the following system of equations:
E) Calculating the structural deformation at the measuring point according to the relation in the step (C) and the step (D):
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,
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:
wherein the content of the first and second substances,represents the distance from the beam structure centerline at the ith test point for fiber number S3 and fiber number S4;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;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;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;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
Step four, calculating a structural deformation-curvature relation according to the following equation:
Step five, calculating the structural deformation at the measuring point according to the relational expression in the step three and the step four:
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.
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. In the invention, the strain sensor is a fiber bragg grating sensor, N (N is 5 in fig. 1) fiber bragg grating sensors are connected in series in one optical fiber, and specifically, the sensor is a great distance(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)=[Ri1Ri2Ri3]T. Column vector in inventive content step fiveCan be defined according to the position vector and the base vector of the coordinate system:
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:
wherein the content of the first and second substances,represents the distance from the beam structure centerline at the ith test point for fiber number S3 and fiber number S4;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;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;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;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
Step four, calculating a structural deformation-curvature relation according to the following equation:
Step five, calculating the structural deformation at the measuring point according to the relational expression in the step three and the step four:
among them, in the above-mentioned case,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.
Showing the structure in the end position (left-most root shown in FIG. 1)) The position vector of (c).
Thus, the three-dimensional position of the beam structure at the measuring point i after deformation is
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)ε5Connected in series to a third fibre S3, sensorThe arrangement directions alternately form 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)=[Ri1Ri2Ri3]T. Column vector in inventive content step fiveCan be defined according to the position vector and the base vector of the coordinate system:
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:
wherein the content of the first and second substances,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;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;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;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;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
D) The structural deformation-curvature relationship is calculated according to the following system of equations:
E) Calculating the structural deformation at the measuring point according to the relation in the step (C) and the step (D):
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,
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)=[Ri1Ri2Ri3]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:
wherein the content of the first and second substances,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;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;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;indicating that the second fiber S2 is at the ith stationA coordinate component of the center line of the beam structure relative to the structure in the Y direction;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
D) The structural deformation-curvature relationship is calculated according to the following system of equations:
E) Calculating the structural deformation at the measuring point according to the relation in the step (C) and the step (D):
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,
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 true CN111380476A (en) | 2020-07-07 |
CN111380476B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113916123A (en) * | 2021-11-22 | 2022-01-11 | 北京航空航天大学 | Strain-based large flexible beam structure three-dimensional space configuration determination method |
Citations (5)
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 |
CN103123254A (en) * | 2013-02-07 | 2013-05-29 | 南开大学 | Two-dimensional bending vector sensor based on fiber grating interference structure |
CN104794284A (en) * | 2015-04-22 | 2015-07-22 | 西安电子科技大学 | Intelligent skin antenna electric compensation method based on embedded fiber bragg grating |
CN105783758A (en) * | 2016-04-07 | 2016-07-20 | 华北电力大学(保定) | Self-heterodyne single-ended vector BOTDA-based dynamic strain measurement method and apparatus |
-
2018
- 2018-12-27 CN CN201811612819.8A patent/CN111380476B/en active Active
Patent Citations (5)
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 |
CN103123254A (en) * | 2013-02-07 | 2013-05-29 | 南开大学 | Two-dimensional bending vector sensor based on fiber grating interference structure |
CN104794284A (en) * | 2015-04-22 | 2015-07-22 | 西安电子科技大学 | Intelligent skin antenna electric compensation method based on embedded fiber bragg grating |
CN105783758A (en) * | 2016-04-07 | 2016-07-20 | 华北电力大学(保定) | Self-heterodyne single-ended vector BOTDA-based dynamic strain measurement method and apparatus |
Non-Patent Citations (2)
Title |
---|
YANG MENG1: "Strain-based Shape Prediction for Flexible", 《AIAA SCITECH FORUM》 * |
安朝: "弹性欧拉梁大变形分析的椭圆积分统一形式", 《北京航空航天大学学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113916123A (en) * | 2021-11-22 | 2022-01-11 | 北京航空航天大学 | Strain-based large flexible beam structure three-dimensional space configuration determination method |
CN113916123B (en) * | 2021-11-22 | 2022-05-27 | 北京航空航天大学 | Strain-based large flexible beam structure three-dimensional space configuration determination method |
Also Published As
Publication number | Publication date |
---|---|
CN111380476B (en) | 2021-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7520176B1 (en) | Method for real-time structure shape-sensing | |
US8970845B1 (en) | In-situ three-dimensional shape rendering from strain values obtained through optical fiber sensors | |
KR102430843B1 (en) | Aircraft monitoring system | |
CN107869960B (en) | A kind of aircraft wing deformation measuring device, installation method and test method | |
CN109977448A (en) | A kind of rotor aerodynamics analysis method based on actual measurement structural loads | |
CN110006563B (en) | Distributed optical fiber decoupling measurement method for blade flapping and shimmy bending moment of helicopter | |
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 | |
Raab et al. | Dynamic flight load measurements with MEMS pressure sensors | |
US7999926B2 (en) | Method and device for determining anemometric parameters of an aircraft | |
Raab et al. | In-flight testing of MEMS pressure sensors for flight loads determination | |
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 | |
CN106323587B (en) | The monocular video high-precision measuring method of wing model in wind tunnel flexible deformation | |
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 | |
CN115169162B (en) | Method and device for predicting airplane vibration environment and computer readable storage medium | |
Luczak et al. | Static and dynamic testing of the full scale helicopter rotor blades | |
RU2396569C1 (en) | Method of determining aerial parametres in flight tests of aircraft flying at high angle of attack | |
RU2375690C1 (en) | Method for determination of pitot probe aerodynamic errors in flight tests of flying vehicle | |
CN112484751B (en) | Method for measuring position and attitude of spacecraft verifier in relatively large space test field coordinate 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 |