Dynamic Strain Field Measuring Method and System for Rotor Blade Based on Blade Tip Timing

Abstract

The present invention discloses a dynamic strain field measuring method and system for a rotor blade based on blade tip timing. The method includes the following steps: establishing a three-dimensional finite element model of a to-be-measured rotor blade, and extracting modal parameters of the three-dimensional finite element model; determining the number and circumferential mounting positions of blade tip timing sensors; establishing a mapping relationship between single-point displacement and full-field dynamic strains of the blade; acquiring blade tip single-point displacement of the rotor blade based on the blade tip timing sensors; and realizing, by the single-point displacement, dynamic strain measurement in any position and direction of the rotor blade based on the mapping relationship.

Claims

1. A dynamic strain field measuring method for a rotor blade based on blade tip timing, comprising the following steps: in the first step (S1), establishing a three-dimensional finite element model of a to-be-measured rotor blade, and extracting modal parameters of the three-dimensional finite element model; in the second step (S2), determining the number and circumferential mounting positions of blade tip timing sensors; in the third step (S3), establishing a mapping relationship between single-point displacement and full-field dynamic strains of the blade; in the fourth step (S4), acquiring blade tip single-point displacement of the rotor blade based on the blade tip timing sensors; and in the fifth step (S5), realizing, by the single-point displacement, dynamic strain measurement in any position and direction of the rotor blade based on the mapping relationship.

2. The method according to claim 1, wherein preferably, in the first step (S1), first n.sub.m order modal parameters: a modal frequency f.sub.i, a displacement modal shape ϕ.sub.i with the size being n.sub.dof×1 and a strain modal shape ψ.sub.i with the size being 2 n.sub.dof×1, of the three-dimensional finite element model are extracted by modal analysis, and a full-field strain modal shape matrix, ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] with the size being 2 n.sub.dof×n.sub.m, of the rotor blade is constructed, wherein n.sub.m represents for a modal number, i represents for a modal order, n.sub.dof represents for the number of degrees of freedom of the finite element model of the rotor blade, n.sub.dof=3 n.sub.n, and n.sub.n represents for the number of nodes of the finite element model of the rotor blade.

3. The method according to claim 2, wherein in the first step (S1), strains of each node of the finite element model of the rotor blade comprise six strain components in total, comprising three positive strains ε.sub.z and three shearing strains γ.sub.xz.

4. The method according to claim 1, wherein in the second step (S2), a relationship between the number n.sub.btt of the blade tip timing sensors mounted in a circumferential direction of a rotor blade casing and the multimodal number n.sub.m of vibration is that: n.sub.btt≥2 n.sub.m+1.

5. The method according to claim 1, wherein in the second step (S2), a measuring point selection matrix S.sub.btt arranged in a circumferential direction of a blade tip timing sensor casing under multimodal excitation is constructed: $S_{\mathrm{btt}}=\left[\begin{array}{cccccc}1& \sin \left(E{O}_1{\theta }_1\right)& \cos \left(E{O}_1{\theta }_1\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_1\right)& \cos \left(E{O}_{n_m}{\theta }_1\right)\\ 1& \sin \left(E{O}_1{\theta }_2\right)& \cos \left(E{O}_1{\theta }_2\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_2\right)& \cos \left(E{O}_{n_m}{\theta }_2\right)\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ 1& \sin \left(E{O}_1{\theta }_{n_{\mathrm{btt}}}\right)& \cos \left(E{O}_1{\theta }_{n_{\mathrm{btt}}}\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_{n_{\mathrm{btt}}}\right)& \cos \left(E{O}_{n_m}{\theta }_{n_{\mathrm{btt}}}\right)\end{array}\right]$ wherein the size of the measuring point selection matrix S.sub.btt is n.sub.btt×(2 n.sub.m+1), θ.sub.j represents for layout angles of the blade tip timing sensors j(j=1, . . . n.sub.btt).sub.th in the casing, n.sub.btt represents for the total number of the j(j=1, . . . n.sub.btt).sub.th blade tip timing sensors arranged in the circumferential direction, and EO.sub.i represents for a concerned excitation order (i=1, . . . n.sub.m); a position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and a measuring point layout when the condition number κ of the matrix is minimum is selected therefrom.

6. The method according to claim 2, wherein in the third step (S3), a conversion matrix, $T=\left[\frac{{\psi }_1}{{\varphi }_{\mathrm{btt},1}},\frac{{\psi }_2}{{\varphi }_{\mathrm{btt},2}},\mathrm{.Math.},\frac{{\psi }_i}{{\varphi }_{\mathrm{btt},i}},\mathrm{.Math.},\frac{{\psi }_{n_m}}{{\varphi }_{\mathrm{btt},n_m}}\right]$ with the size being 2 n.sub.dof×n.sub.d, of the single-point displacement and the full-field dynamic strains of the blade is constructed, wherein ϕ.sub.btt,i represents for an i.sup.th order displacement modal shape of a blade tip timing measuring point of the blade.

7. The method according to claim 6, wherein in the fourth step (S4), blade tip multimodal vibration signals u.sub.btt(t) on the N.sup.th turn of the rotor blade are acquired by n.sub.btt blade tip timing sensors, and furthermore, an n.sub.m order vibration parameter β=(S.sub.btt).sup.†u.sub.btt(t) is acquired by utilizing a circumferential Fourier fitting algorithm, wherein β=[c,A.sub.1 cos(φ.sub.1),A.sub.1 sin(φ.sub.1), . . . ,A.sub.n.sub.m cos(φ.sub.n.sub.m),A.sub.n.sub.m sin (φ.sub.n.sub.m)].sup.T, A.sub.i represents for a vibration amplitude, f.sub.i represents for a modal frequency, and φ.sub.i represents for an initial phase; the superscript † represents for the inversion of the matrix; the superscript T represents for the transposition of a vector; and further, blade tip single-point multimodal vibration signals are acquired by reconstruction: $u_{\mathrm{btt}}\left(t\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{u}_{\mathrm{btt},i}\left(t\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{A}_i\sin \left({\omega }_{i}t+{\phi }_i\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{A}_i\cos \left({\phi }_i\right)\sin \left({\omega }_{i}t\right)+A_i\sin \left({\phi }_i\right)\cos \left({\omega }_{i}t\right)$ wherein u.sub.btt,i(t) represents for decoupled i.sup.th order vibration signals, c represents for static deformation of the blade, ω.sub.i represents for a circular frequency of multimodal vibration of the blade, and t represents for a vibration moment of the blade.

8. The method according to claim 7, wherein in the fifth step (S5), based on the conversion matrix T and the decoupled multimodal vibration signals u.sub.btt,i(t), strains of all the nodes on a surface of and inside the blade at each turn of the rotor blade or at the moment t are calculated according to a formula S(t)=T[u.sub.btt,1,u.sub.btt,2, . . . ,u.sub.btt,i, . . .,u.sub.btt,n.sub.m].sup.T to obtain S(t)=[ε.sub.1,x,ε.sub.1,y,ε.sub.1,z,γ.sub.1,xy,γ.sub.1,yz,γ.sub.1,xz, . . . ,ε.sub.i,x,ε.sub.i,y,ε.sub.i,z,γ.sub.i,xy,γ.sub.i,yz,γ.sub.i,xz, . . . ,ε.sub.n.sub.n.sub.,x,ε.sub.n.sub.n.sub.,y,ε.sub.n.sub.n.sub.,z,γ.sub.n.sub.n.sub.,xy,γ.sub.n.sub.n.sub.,yz,γ.sub.n.sub.n.sub.,xz].sup.T, and the strains S(t) comprise positive strains and shearing strains, wherein ε.sub.i,x represents for a positive strain of the i.sup.th node of the finite element model of the blade in an x direction, ε.sub.i,y represents for a positive strain of the i.sup.th node of the finite element model of the blade in a y direction, ε.sub.i,z represents for a positive strain of the i.sup.th node of the finite element model of the blade in a z direction, γ.sub.i,xy represents for a shearing strain of the i.sup.th node of the finite element model of the blade in an x-y direction, γ.sub.i,yz represents for a shearing strain of the i.sup.th node of the finite element model of the blade in a y-z direction, γ.sub.i,xz represents for a shearing strain of the i.sup.th node of the finite element model of the blade in an x-z direction, and the superscript T represents for the transposition of a vector.

9. A measuring system for implementing the method according to claim 1, comprising: a plurality of blade tip timing sensors arranged on a rotor blade casing; a blade tip timing vibration measuring module connected with the blade tip timing sensors so as to measure circumferential blade tip multimodal vibration signals of a rotor blade; and a calculation unit connected with the blade tip timing vibration measuring module, and the calculation unit comprising: a modal analysis module configured to perform modal analysis based on a three-dimensional finite element model of a to-be-measured rotor blade so as to acquire a displacement modal shape ϕ.sub.i of first n.sub.m order modal frequency f.sub.i of the rotor blade and construct a full-field strain modal matrix ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] of the rotor blade; a measuring point optimization module configured to optimize the number of measuring points of the blade tip timing sensors arranged on the rotor blade, wherein a measuring point selection matrix S.sub.btt arranged in the circumferential direction of the blade tip timing sensor casing under multimodal excitation is constructed, a position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and a measuring point layout when the condition number κ of the matrix is minimum is selected therefrom; a conversion matrix calculation module configured to construct a conversion matrix of single-point displacement and full-field dynamic strains of the blade; and a strain field reconstruction module configured to calculate positive strains and shearing strains of all nodes on the surface of and inside the blade at each turn of the rotor blade.

10. The measuring system according to claim 9, wherein the blade tip timing vibration measuring module comprises the rotating speed sensor, a signal conditioning module and a time-displacement conversion module, wherein the rotating speed sensor is configured to measure the rotating speed of the blade; the time-displacement conversion module is configured to acquire blade tip vibration displacement of the blade based on a speed signal; and the signal conditioning module is configured to extract a vibration frequency and amplitude value parameters of the blade from vibration displacement.

11. A measuring system for implementing the method according to claim 2, comprising: a plurality of blade tip timing sensors arranged on a rotor blade casing; a blade tip timing vibration measuring module connected with the blade tip timing sensors so as to measure circumferential blade tip multimodal vibration signals of a rotor blade; and a calculation unit connected with the blade tip timing vibration measuring module, and the calculation unit comprising: a modal analysis module configured to perform modal analysis based on a three-dimensional finite element model of a to-be-measured rotor blade so as to acquire a displacement modal shape ϕ.sub.i of first n.sub.m order modal frequency f.sub.i of the rotor blade and construct a full-field strain modal matrix ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] of the rotor blade; a measuring point optimization module configured to optimize the number of measuring points of the blade tip timing sensors arranged on the rotor blade, wherein a measuring point selection matrix S.sub.btt arranged in the circumferential direction of the blade tip timing sensor casing under multimodal excitation is constructed, a position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and a measuring point layout when the condition number κ of the matrix is minimum is selected therefrom; a conversion matrix calculation module configured to construct a conversion matrix of single-point displacement and full-field dynamic strains of the blade; and a strain field reconstruction module configured to calculate positive strains and shearing strains of all nodes on the surface of and inside the blade at each turn of the rotor blade.

12. A measuring system for implementing the method according to claim 3, comprising: a plurality of blade tip timing sensors arranged on a rotor blade casing; a blade tip timing vibration measuring module connected with the blade tip timing sensors so as to measure circumferential blade tip multimodal vibration signals of a rotor blade; and a calculation unit connected with the blade tip timing vibration measuring module, and the calculation unit comprising: a modal analysis module configured to perform modal analysis based on a three-dimensional finite element model of a to-be-measured rotor blade so as to acquire a displacement modal shape ϕ.sub.i of first n.sub.m order modal frequency f.sub.i of the rotor blade and construct a full-field strain modal matrix ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] of the rotor blade; a measuring point optimization module configured to optimize the number of measuring points of the blade tip timing sensors arranged on the rotor blade, wherein a measuring point selection matrix S.sub.btt arranged in the circumferential direction of the blade tip timing sensor casing under multimodal excitation is constructed, a position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and a measuring point layout when the condition number κ of the matrix is minimum is selected therefrom; a conversion matrix calculation module configured to construct a conversion matrix of single-point displacement and full-field dynamic strains of the blade; and a strain field reconstruction module configured to calculate positive strains and shearing strains of all nodes on the surface of and inside the blade at each turn of the rotor blade.

13. A measuring system for implementing the method according to claim 4, comprising: a plurality of blade tip timing sensors arranged on a rotor blade casing; a blade tip timing vibration measuring module connected with the blade tip timing sensors so as to measure circumferential blade tip multimodal vibration signals of a rotor blade; and a calculation unit connected with the blade tip timing vibration measuring module, and the calculation unit comprising: a modal analysis module configured to perform modal analysis based on a three-dimensional finite element model of a to-be-measured rotor blade so as to acquire a displacement modal shape ϕ.sub.i of first n.sub.m order modal frequency f.sub.i of the rotor blade and construct a full-field strain modal matrix ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] of the rotor blade; a measuring point optimization module configured to optimize the number of measuring points of the blade tip timing sensors arranged on the rotor blade, wherein a measuring point selection matrix S.sub.btt arranged in the circumferential direction of the blade tip timing sensor casing under multimodal excitation is constructed, a position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and a measuring point layout when the condition number κ of the matrix is minimum is selected therefrom; a conversion matrix calculation module configured to construct a conversion matrix of single-point displacement and full-field dynamic strains of the blade; and a strain field reconstruction module configured to calculate positive strains and shearing strains of all nodes on the surface of and inside the blade at each turn of the rotor blade.

14. A measuring system for implementing the method according to claim 5, comprising: a plurality of blade tip timing sensors arranged on a rotor blade casing; a blade tip timing vibration measuring module connected with the blade tip timing sensors so as to measure circumferential blade tip multimodal vibration signals of a rotor blade; and a calculation unit connected with the blade tip timing vibration measuring module, and the calculation unit comprising: a modal analysis module configured to perform modal analysis based on a three-dimensional finite element model of a to-be-measured rotor blade so as to acquire a displacement modal shape ϕ.sub.i of first n.sub.m order modal frequency f.sub.i of the rotor blade and construct a full-field strain modal matrix ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] of the rotor blade; a measuring point optimization module configured to optimize the number of measuring points of the blade tip timing sensors arranged on the rotor blade, wherein a measuring point selection matrix S.sub.btt arranged in the circumferential direction of the blade tip timing sensor casing under multimodal excitation is constructed, a position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and a measuring point layout when the condition number κ of the matrix is minimum is selected therefrom; a conversion matrix calculation module configured to construct a conversion matrix of single-point displacement and full-field dynamic strains of the blade; and a strain field reconstruction module configured to calculate positive strains and shearing strains of all nodes on the surface of and inside the blade at each turn of the rotor blade.

15. A measuring system for implementing the method according to claim 6, comprising: a plurality of blade tip timing sensors arranged on a rotor blade casing; a blade tip timing vibration measuring module connected with the blade tip timing sensors so as to measure circumferential blade tip multimodal vibration signals of a rotor blade; and a calculation unit connected with the blade tip timing vibration measuring module, and the calculation unit comprising: a modal analysis module configured to perform modal analysis based on a three-dimensional finite element model of a to-be-measured rotor blade so as to acquire a displacement modal shape ϕ.sub.i of first n.sub.m order modal frequency f.sub.i of the rotor blade and construct a full-field strain modal matrix ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] of the rotor blade; a measuring point optimization module configured to optimize the number of measuring points of the blade tip timing sensors arranged on the rotor blade, wherein a measuring point selection matrix S.sub.btt arranged in the circumferential direction of the blade tip timing sensor casing under multimodal excitation is constructed, a position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and a measuring point layout when the condition number κ of the matrix is minimum is selected therefrom; a conversion matrix calculation module configured to construct a conversion matrix of single-point displacement and full-field dynamic strains of the blade; and a strain field reconstruction module configured to calculate positive strains and shearing strains of all nodes on the surface of and inside the blade at each turn of the rotor blade.

16. A measuring system for implementing the method according to claim 7, comprising: a plurality of blade tip timing sensors arranged on a rotor blade casing; a blade tip timing vibration measuring module connected with the blade tip timing sensors so as to measure circumferential blade tip multimodal vibration signals of a rotor blade; and a calculation unit connected with the blade tip timing vibration measuring module, and the calculation unit comprising: a modal analysis module configured to perform modal analysis based on a three-dimensional finite element model of a to-be-measured rotor blade so as to acquire a displacement modal shape ϕ.sub.i of first n.sub.m order modal frequency f.sub.i of the rotor blade and construct a full-field strain modal matrix ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] of the rotor blade; a measuring point optimization module configured to optimize the number of measuring points of the blade tip timing sensors arranged on the rotor blade, wherein a measuring point selection matrix S.sub.btt arranged in the circumferential direction of the blade tip timing sensor casing under multimodal excitation is constructed, a position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and a measuring point layout when the condition number κ of the matrix is minimum is selected therefrom; a conversion matrix calculation module configured to construct a conversion matrix of single-point displacement and full-field dynamic strains of the blade; and a strain field reconstruction module configured to calculate positive strains and shearing strains of all nodes on the surface of and inside the blade at each turn of the rotor blade.

17. A measuring system for implementing the method according to claim 8, comprising: a plurality of blade tip timing sensors arranged on a rotor blade casing; a blade tip timing vibration measuring module connected with the blade tip timing sensors so as to measure circumferential blade tip multimodal vibration signals of a rotor blade; and a calculation unit connected with the blade tip timing vibration measuring module, and the calculation unit comprising: a modal analysis module configured to perform modal analysis based on a three-dimensional finite element model of a to-be-measured rotor blade so as to acquire a displacement modal shape ϕ.sub.i of first n.sub.m order modal frequency f.sub.i of the rotor blade and construct a full-field strain modal matrix ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] of the rotor blade; a measuring point optimization module configured to optimize the number of measuring points of the blade tip timing sensors arranged on the rotor blade, wherein a measuring point selection matrix S.sub.btt arranged in the circumferential direction of the blade tip timing sensor casing under multimodal excitation is constructed, a position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and a measuring point layout when the condition number κ of the matrix is minimum is selected therefrom; a conversion matrix calculation module configured to construct a conversion matrix of single-point displacement and full-field dynamic strains of the blade; and a strain field reconstruction module configured to calculate positive strains and shearing strains of all nodes on the surface of and inside the blade at each turn of the rotor blade.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0029] By reading the detailed description in the preferred specific implementation manners hereunder, various other advantages and benefits of the present invention will become clear for the ordinary skilled in the art. Accompanying drawings of the specification are merely intended to show the preferred implementation manners, rather than to limit the present invention. Obviously, the accompanying drawings described below are only some embodiments of the present invention, and other accompanying drawings can also be acquired by the ordinary skilled in the art based on the accompanying drawings without creative labor. Moreover, the same components are marked by the same symbols in the overall accompanying drawings.

[0030] In the accompanying drawings:

[0031] FIG. 1 is a schematic process diagram of a preferred example of a dynamic strain field measuring method for a rotor blade based on blade tip timing, provided by the present invention;

[0032] FIG. 2(a) to FIG. 2(b) are schematic structural diagrams of a dynamic strain field measuring system for the rotor blade based on blade tip timing, provided by the present invention, wherein FIG. 2(a) shows a composition of a dynamic strain field reconstruction system for the rotor blade; and FIG. 2(b) shows a circumferential mounting schematic diagram of a blade tip timing vibration measuring module and blade tip timing sensors;

[0033] FIG. 3 is a schematic diagram of simulating a dynamic load excitation position of the rotor blade and measuring point positions of blade tip timing sensors (BTT) in one embodiment;

[0034] FIG. 4(a) to FIG. 4(f) are a displacement modal shape and a strain modal shape of the rotor blade in one embodiment, wherein FIG. 4(a) is a first flexural displacement shape; FIG. 4(b) is a first flexural strain shape; FIG. 4(c) is a first torsional displacement shape; FIG. 4(d) is a first torsional strain shape; FIG. 4(e) is a second flexural displacement shape; and FIG. 4(f) is a second flexural strain shape;

[0035] FIG. 5 is blade tip displacement vibration signals actually measured by seven blade tip timing sensors of the rotor blade in one embodiment;

[0036] FIG. 6 is a multimodal decoupling result of a displacement signal of the rotor blade in one embodiment;

[0037] FIG. 7(a) to FIG. 7(b) are comparison results of the dynamic strain of a No. 748 node on a blade body in a reconstructed strain field of the rotor blade in one embodiment and a true dynamic strain, wherein FIG. 7(a) is a comparison result of three positive strain components and the true dynamic strain; and FIG. 7(b) is a comparison result of three shearing strain components and the true dynamic strain;

[0038] FIG. 8(a) to FIG. 8(b) are comparison results of the dynamic strain of a No. 1100 node on a blade root in a reconstructed strain field of the rotor blade in one embodiment and the true dynamic strain, wherein FIG. 8(a) is a comparison result of three positive strain components and the true dynamic strain; and FIG. 8(b) is a comparison result of three shearing strain components and the true dynamic strain.

[0039] The present invention is further explained below in combination with the accompanying drawings and embodiments.

DETAILED DESCRIPTION

[0040] The specific embodiments of the present invention will be described in more detail below by reference to accompanying drawings 1-8(b). Although the specific embodiments of the present invention are shown in the accompanying drawings, it should be understood that the present invention can be implemented in various forms, but should not be limited by the embodiments described herein. On the contrary, these embodiments are provided in order to more thoroughly understand the present invention and integrally deliver the scope of the present invention to the skilled in the art.

[0041] It should be noted that some words are used in the specification and claims to denote specific components. It should be understood by the skilled in the art that a technician may call the same component with different nouns. The functional difference of the component, rather than the difference of the nouns, is used as a distinction criterion in the specification and claims. For example, “including” or “includes” mentioned in the overall specification and claims is generally an open word so as to be explained as “including, but not being limited to”. The subsequent description of the specification is a preferred implementation manner of the present invention, however, the description is intended to conform to a general principle of the specification, rather than to limit the scope of the present invention. The protective scope of the present invention should be based on the scope defined by the appended claims.

[0042] In order to facilitate understanding the embodiments of the present invention, the specific embodiments serving as examples are further explained below in combination with the accompanying drawings, and each accompanying drawing is not regarded as a limit to the embodiments of the present invention.

[0043] In order to realize better understanding, FIG. 1 is a working process diagram of a dynamic strain field measuring method for a rotor blade based on blade tip timing, as shown in FIG. 1, the dynamic strain field measuring method for the rotor blade based on blade tip timing includes the following steps:

[0044] in the first step S1, a three-dimensional finite element model of the to-be-measured rotor blade is established, and modal parameters of the three-dimensional finite element model are extracted;

[0045] in the second step S2, the number and circumferential mounting positions of blade tip timing sensors are determined;

[0046] in the third step S3, a mapping relationship between single-point displacement and full-field dynamic strains of the blade is established;

[0047] in the fourth step S4, blade tip single-point displacement of the rotor blade is acquired based on the blade tip timing sensors; and

[0048] in the fifth step S5, dynamic strain measurement in any position and direction of the rotor blade is realized by the single-point displacement based on the mapping relationship.

[0049] In one implementation manner of the method, in the first step S1, first n.sub.m order modal parameters: the modal frequency f.sub.i, a displacement modal shape ϕ.sub.i with the size being n.sub.dof×1 and a strain modal shape ψ.sub.i with the size being 2 n.sub.dof×1, of the three-dimensional finite element model are extracted by modal analysis, and a full-field strain modal shape matrix, ψ=[ψ.sub.1,ψ.sub.2, . . . .ψ.sub.n.sub.m] with the size being 2 n.sub.dofn.sub.m, of the rotor blade is constructed, wherein n.sub.m represents for the modal number, i represents for the modal order, n.sub.dof represents for the number of degrees of freedom of the finite element model of the rotor blade, n.sub.dof=3 n.sub.n, and n.sub.n represents for the number of nodes of the finite element model of the rotor blade.

[0050] In one implementation manner of the method, in the first step S1, strains of each nodes of the finite element model of the rotor blade include six strain components in total, including three positive strains ε.sub.z and three shearing strains γ.sub.xz.

[0051] In one implementation manner of the method, in the second step S2, a relationship between the number n.sub.btt of the blade tip timing sensors mounted in the circumferential direction of a rotor blade casing and the multimodal number n.sub.m of vibration is that: n.sub.btt≥2 n.sub.m+1.

[0052] In one implementation manner of the method, in the second step S2, a measuring point selection matrix S.sub.btt arranged in the circumferential direction of a blade tip timing sensor casing under multimodal excitation is constructed:

[00004]$S_{\mathrm{btt}}=\left[\begin{array}{cccccc}1& \sin \left(E{O}_1{\theta }_1\right)& \cos \left(E{O}_1{\theta }_1\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_1\right)& \cos \left(E{O}_{n_m}{\theta }_1\right)\\ 1& \sin \left(E{O}_1{\theta }_2\right)& \cos \left(E{O}_1{\theta }_2\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_2\right)& \cos \left(E{O}_{n_m}{\theta }_2\right)\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ 1& \sin \left(E{O}_1{\theta }_{n_{\mathrm{btt}}}\right)& \cos \left(E{O}_1{\theta }_{n_{\mathrm{btt}}}\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_{n_{\mathrm{btt}}}\right)& \cos \left(E{O}_{n_m}{\theta }_{n_{\mathrm{btt}}}\right)\end{array}\right]$

wherein the size of the measuring point selection matrix S.sub.btt is n.sub.btt×(2 n.sub.m+1) , EO.sub.i represents for the concerned excitation order; the position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and the measuring point layout when the condition number κ of the matrix is minimum is selected therefrom.

[0053] In one implementation manner of the method, in the third step S3, a conversion matrix,

[00005]$T=\left[\frac{{\psi }_1}{{\varphi }_{\mathrm{btt},1}},\frac{{\psi }_2}{{\varphi }_{\mathrm{btt},2}},\mathrm{.Math.},\frac{{\psi }_i}{{\varphi }_{\mathrm{btt},i}},\mathrm{.Math.},\frac{{\psi }_{n_m}}{{\varphi }_{\mathrm{btt},n_m}}\right]$

with the size being 2 n.sub.dof×n.sub.d, of the single-point displacement and the full-field dynamic strains of the blade is constructed, wherein ϕ.sub.btt,i represents for an i.sup.th order displacement modal shape of a blade tip timing measuring point.

[0054] In one implementation manner of the method, in the fourth step S4, blade tip multimodal vibration signals u.sub.btt(t) on the N.sup.th turn of the rotor blade is acquired by the n.sub.btt blade tip timing sensors, and furthermore, an n.sub.m order vibration parameter is acquired by utilizing a circumferential Fourier fitting algorithm,

[0055] wherein β=[c,A.sub.1 cos(φ.sub.1),A.sub.1 sin(φ.sub.1), . . . ,A.sub.n.sub.m cos(φ.sub.n.sub.m),A.sub.n.sub.m sin (φ.sub.n.sub.m)].sup.T, A.sub.i represents for a vibration amplitude, f.sub.i represents for a modal frequency, and φ.sub.i represents for an initial phase; the superscript † represents for the inversion of the matrix; the superscript T represents for the transposition of a vector; and further, blade tip single-point multimodal vibration signals is acquired by reconstruction:

[00006]$u_{\mathrm{btt}}\left(t\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{u}_{\mathrm{btt},i}\left(t\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{A}_i\sin \left({\omega }_{i}t+{\phi }_i\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{A}_i\cos \left({\phi }_i\right)\sin \left({\omega }_{i}t\right)+A_i\sin \left({\phi }_i\right)\cos \left({\omega }_{i}t\right)$

[0056] wherein u.sub.btt,i(t) represents for the decoupled i.sup.th order vibration signals, c represents for static deformation of the blade, ω.sub.i represents for a circular frequency of the multimodal vibration of the blade, and t represents for a vibration moment of the blade.

[0057] In one implementation manner of the method, in the fifth step S5, the strains of all the nodes on the surface of and inside the blade at each turn of the rotor blade are calculated according to a formula S(t)=T[u.sub.btt,1,u.sub.btt,2, . . . ,u.sub.btt,i, . . .,u.sub.btt,n.sub.m].sup.T to obtain S(t)=[ε.sub.1,x,ε.sub.1,y,ε.sub.1,z,γ.sub.1,xy,γ.sub.1,yz,γ.sub.1,xz, . . . ,ε.sub.i,x,ε.sub.i,y,ε.sub.i,z,γ.sub.i,xy,γ.sub.i,yz,γ.sub.i,xz, . . . ,ε.sub.n.sub.n.sub.,x,ε.sub.n.sub.n.sub.,y,ε.sub.n.sub.n.sub.,z,γ.sub.n.sub.n.sub.,xy,γ.sub.n.sub.n.sub.,yz,γ.sub.n.sub.n.sub.,xz].sup.T, and the strains S(t) include the positive strains and the shearing strains, wherein ε.sub.i,x represents for the positive strain of the i.sup.th node of the finite element model of the blade in an x direction, ε.sub.i,y represents for the positive strain of the i.sup.th node of the finite element model of the blade in a y direction, ε.sub.i,z represents for the positive strain of the i.sup.th node of the finite element model of the blade in a z direction, γ.sub.i,xy represents for the shearing strain of the i.sup.th node of the finite element model of the blade in an x-y direction, γ.sub.i,yz represents for the shearing strain of the i.sup.th node of the finite element model of the blade in a y-z direction, represents for the shearing strain of the i.sup.th node of the finite element model of the blade in an x-z direction, and the superscript T represents for the transposition of a vector.

[0058] In order to further understand the present invention, the present invention is further described below in combination with the accompanying drawings 1-8 and the specific embodiments, and it should be emphasized that the following description is merely exemplary, while application objects of the present invention are not limited to the following examples.

[0059] FIG. 1 is a schematic process diagram of the dynamic strain field measuring method for the rotor blade based on blade tip timing, completed by the present invention, and according to the method, blade tip single-point vibration is measured by using blade tip timing sensors mounted in the circumferential direction of the casing, multimodal vibration decoupling is realized by utilizing the circumferential Fourier fitting method, and a conversion relationship between blade tip displacement measuring points of the rotor blade and the strains of all the nodes in a full field is established, so that reconstruction of the dynamic strain field of the blade is realized; and FIG. 2(a) to FIG. 2(b) are schematic structural diagrams of a dynamic strain field measuring system for the rotor blade based on blade tip timing, provided by the present invention, wherein 1-blade tip timing sensor, 2-rotor casing, 3-strain gauge, 4-rotor blade, 5-wheel disc, 6-rotor, and 7-rotating speed sensor, and the method includes the specific steps as follows: [0060] 1) the modal parameters of the three-dimensional finite element model of the blade are extracted: referring to FIG. 3, the three-dimensional finite element model simulating the rotor straight-panel blade is established by utilizing ANSYS finite element analysis software, wherein the three-dimensional finite element model is made of aluminum with the density being 2700 kg/m<3>, the Poisson ratio being 0.33 and the elasticity modulus being 72000 Mpa; the blade has the length of 48 mm, the thickness of 1 mm and the width of 20 mm; a finite element is an entity unit SOLID185 and has 3153 nodes in total; and the two side surfaces of a blade root are fixedly restricted to simulate the actual working state of the rotor blade; [0061] the first three order modal parameters, namely n.sub.m=3: modal frequencies f.sub.i, a displacement modal shape ϕ.sub.i with the size being n.sub.dof×1 and a strain modal shape ψ.sub.i with the size being 2 n.sub.dof×1, are extracted by ANSYS modal analysis, wherein the first three order modal frequencies are respectively f.sub.1=333.08 Hz, f.sub.2=1806.03 Hz and f.sub.3=2076.52 Hz; and a full-field strain modal shape matrix, ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] with the size being 2 n.sub.dof×n.sub.m, of the rotor blade is constructed, and the shape refers to FIG. 4(a) to FIG. 4(f), wherein i represents for the modal order, n.sub.dof=9459 represents for the number of degrees of freedom of the finite element model of the blade, n.sub.dof=3 n.sub.n, and n=3153 represents for the number of nodes of the finite element model of the blade; and the strains include the six strain components in total, including the three positive strains ε.sub.z and the three shearing strains γ.sub.xz, namely each node has the six strain modal shapes. [0062] 2) The number and circumferential mounting positions of the blade tip timing sensors are determined: the relationship between the number n.sub.btt of the blade tip timing sensors mounted in the circumferential direction of the rotor blade casing and the multimodal number n.sub.m of vibration is that: n.sub.btt≥n.sub.m+1; in the solution, first three order vibration modes of the simulated rotor blade are concerned, and n.sub.m=3; the number of the blade tip timing sensors in the circumferential direction of the casing is at least n.sub.btt=7;
the measuring point selection matrix S.sub.btt arranged in the circumferential direction of the blade tip timing sensor casing under multimodal excitation is constructed:

[00007]$S_{\mathrm{btt}}=\left[\begin{array}{cccccc}1& \sin \left(E{O}_1{\theta }_1\right)& \cos \left(E{O}_1{\theta }_1\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_1\right)& \cos \left(E{O}_{n_m}{\theta }_1\right)\\ 1& \sin \left(E{O}_1{\theta }_2\right)& \cos \left(E{O}_1{\theta }_2\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_2\right)& \cos \left(E{O}_{n_m}{\theta }_2\right)\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ 1& \sin \left(E{O}_1{\theta }_{n_{\mathrm{btt}}}\right)& \cos \left(E{O}_1{\theta }_{n_{\mathrm{btt}}}\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_{n_{\mathrm{btt}}}\right)& \cos \left(E{O}_{n_m}{\theta }_{n_{\mathrm{btt}}}\right)\end{array}\right]$

wherein the size of the measuring point selection matrix S.sub.btt is n.sub.btt×(2 n.sub.m+1)=7×7, θ.sub.j represents for the layout angles of the blade tip timing sensors j(j=1, . . . n.sub.btt) in the casing, n.sub.btt represents for the total number of the j(j=1, . . . n.sub.btt).sub.th blade tip timing sensors arranged in the circumferential direction, and EO.sub.i represents for the concerned excitation order (i=1, . . . n.sub.m); θ.sub.j represents for the mounting angle of the j.sup.th blade tip timing sensor in the casing; the three excitation orders concerned in the solution are respectively 4, 18 and 23, and the first three order vibration modes of the rotor blade are simultaneously excited at the same rotating speed; the position where the rotating speed sensor is located serves as the reference 0°, the limited mounting angle range of the casing is removed, and the n.sub.btt=7 angles are randomly selected in the circumferential direction of the casing to serve as the mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated; the random process is repeated for R=500 times, and a measuring point layout solution when the condition number κ of the matrix is minimum is selected therefrom. The circumferential mounting angles of the selected seven blade tip timing sensors are 3.12°, 117.33°, 183.42°, 189.58°, 303.71°, 315.14° and 351.08°, and the condition number of the corresponding measuring point selection matrix S.sub.btt is 3.4694. [0063] 3) the mapping relationship between the single-point displacement and the full-field dynamic strains of the blade is established: a conversion matrix,

[00008]$T=\left[\frac{{\psi }_1}{{\varphi }_{\mathrm{btt},1}},\frac{{\psi }_2}{{\varphi }_{\mathrm{btt},2}},\mathrm{.Math.},\frac{{\psi }_i}{{\varphi }_{\mathrm{btt},i}},\mathrm{.Math.},\frac{{\psi }_{n_m}}{{\varphi }_{\mathrm{btt},n_m}}\right]$

with the size being 2 n.sub.dof×n.sub.d=18918×3, of the single-point displacement and the full-field dynamic strains of the blade is constructed, wherein θ.sub.btt,i represents for the i.sup.th order displacement modal shape of the blade tip timing measuring point. [0064] 4) The blade tip single-point displacement of the rotor blade is acquired based on the blade tip timing sensors: the simulated rotor blade is subjected to transient analysis in the ANSYS finite element software, a mass damping coefficient is set as α=12.1380 , a stiffness damping coefficient is set as β=8.1986×10.sup.−8 a rotating speed is set to be 15000 RPM, multimodal vibration of an aerodynamic load to the rotor blade is simulated, multi-frequency harmonic excitation f(t)=cos(2πf.sub.1t)+10 cos(2πf.sub.2t)+20 cos(2πf.sub.3t) is applied to an X direction of a No. 1117 node of a blade tip, and a strain field which is subjected to transient analysis to be stabilized is taken as a reference of a reconstruction result; the blade tip multimodal vibration signals u.sub.btt(t) on the N.sup.th turn of the rotor blade is acquired by the n.sub.btt=7 blade tip timing sensors, referring to FIG. 5, wherein the blade tip timing sensors sample seven data of the vibration signal in the X-direction of the No. 1122 node of the blade tip every time when the rotor blade rotates for a cycle, 25 cycles of 175 data are acquired in total, and therefore, blade tip timing signals are seriously undersampled; and meanwhile, FIG. 5 shows that the data length of the vibration signal in the X-direction of the No. 1122 node of the blade tip under the sampling frequency f.sub.s=25000 Hz is N=2500 , and the sampling time is t=N/f.sub.s0.1 s.

[0065] The n.sub.m in order vibration parameter β=(S.sub.btt).sup.†u.sub.btt(t) is acquired by utilizing the circumferential Fourier fitting algorithm,

[0066] wherein β=[c,A.sub.1 cos(φ.sub.1),A.sub.1 sin(φ.sub.1), . . . ,A.sub.n.sub.m cos(φ.sub.n.sub.m),A.sub.n.sub.m sin (φ.sub.n.sub.m)].sup.T, A.sub.i represents for the vibration amplitude, J.sub.i represents for the modal frequency, and represents for the initial phase; the superscript represents for the inversion of the matrix; the superscript T represents for the transposition of the vector; and further, the blade tip single-point multimodal vibration signals are acquired by reconstruction:

[00009]$u_{\mathrm{btt}}\left(t\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{u}_{\mathrm{btt},i}\left(t\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{A}_i\sin \left({\omega }_{i}t+{\phi }_i\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{A}_i\cos \left({\phi }_i\right)\sin \left({\omega }_{i}t\right)+A_i\sin \left({\phi }_i\right)\cos \left({\omega }_{i}t\right)$

wherein u.sub.btt,i(t) represents for the decoupled i.sup.th order vibration signals, c represents for static deformation of the blade, ω.sub.i represents for the circular frequency of multimodal vibration of the blade, and t represents for the vibration moment of the blade. FIG. 6 gives first three order vibration modal decoupling results of a displacement signal of the rotor blade in the embodiment. [0067] 5) Strains of all the nodes on the surface of and inside the blade at each turn of the rotor blade are calculated according to the formula S(t)=T[u.sub.btt,1,u.sub.btt,2,u.sub.btt,i,u.sub.btt,n.sub.m].sup.T to obtain S(t)=[ε.sub.1,x,ε.sub.1,y,ε.sub.1,z,γ.sub.1,xy,γ.sub.1,yz,γ.sub.1,xz, . . . ,ε.sub.i,x,ε.sub.i,y,ε.sub.i,z,γ.sub.i,xy,γ.sub.i,yz,γ.sub.i,xz, . . . ,ε.sub.n.sub.n.sub.,x,ε.sub.n.sub.n.sub.,y,ε.sub.n.sub.n.sub.,z,γ.sub.n.sub.n.sub.,xy,γ.sub.n.sub.n.sub.,yz,γ.sub.n.sub.n.sub.,xz].sup.T,
wherein the strains S(t) include the positive strains and the shearing strains.

[0068] With the No. 748 node on a blade body and the No. 1100 node on a blade root of the rotor blade as typical examples (referring to FIG. 3) of high-precision reconstruction of the dynamic strain field, the conclusion is also suitable for the other nodes. FIG. 7(a) to FIG. 7(b) are comparison results of the dynamic strains of the No. 748 node on the blade body in a reconstructed strain field of the rotor blade in one embodiment and true dynamic strains; FIG. 8(a) to FIG. 8(b) are comparison results of the dynamic strains of the No. 1100 node on the blade root in the reconstructed strain field of the rotor blade in one embodiment and the true dynamic strains. Known from FIG. 7(a) to FIG. 7(b) and FIG. 8(a) to FIG. 8(b), reconstructed dynamic strain signals are highly consistent with the true dynamic strains; in order to quantitatively evaluate the performance of the dynamic strain field measuring method for the rotor blade in the present invention, a relative error of the reconstructed signals and the true strains is calculated within the range t∈[0,0.1] s, relative errors of the three positive strains of the No. 748 node on the blade body in FIG. 7(a) are respectively 4.41%, 4.69% and 5.77%, and relative errors of the three shearing strains γ.sub.xz of the No. 748 node on the blade body in FIG. 7(b) are respectively 6.24%, 6.21% and 9.63%; and relative errors of the three positive strains ε.sub.z of the No. 1100 node on the blade root in FIG. 8(a) are respectively 4.41%, 4.41% and 4.41%, and relative errors of the three shearing strains γ.sub.xz of the No. 1100 node on the blade root in FIG. 8(b) are respectively 4.40%, 7.81% and 8.09%. Therefore, the dynamic strain field measuring method for the rotor blade based on blade tip timing, provided by the present invention, may be used for high-precision reconstruction of the dynamic strain field of the blade.

[0069] According to the method provided by the present invention, the reconstruction of the overall dynamic strain field of the rotor blade is realized by only utilizing the blade tip finite measuring points, and the measurement of the positive strains and the shearing strains of all the nodes on the surface of and inside the rotor blade under multimodal vibration can be realized, the method is simple in calculation process, and online measurement is easy. The above descriptions are only preferred embodiments of the present invention, and the measuring method can be applied to vibration tests of fans/gas compressors/turbine blades of rotating machinery such as an aircraft engine, a gas turbine and a steam turbine, but is not intended to limit the present invention.

[0070] In another embodiment, the method includes the following steps: [0071] 1) modal parameters of a three-dimensional finite element model of a blade are extracted; [0072] 2) the number and circumferential mounting positions of blade tip timing sensors are determined; [0073] 3) a mapping relationship between single-point displacement and full-field dynamic strains of the blade is established; [0074] 4) blade tip single-point displacement of the rotor blade is acquired by utilizing the blade tip timing sensors; and [0075] 5) the dynamic strain measurement in any position and direction of the rotor

[00010]$S_{\mathrm{btt}}=\left[\begin{array}{cccccc}1& \sin \left(E{O}_1{\theta }_1\right)& \cos \left(E{O}_1{\theta }_1\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_1\right)& \cos \left(E{O}_{n_m}{\theta }_1\right)\\ 1& \sin \left(E{O}_1{\theta }_2\right)& \cos \left(E{O}_1{\theta }_2\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_2\right)& \cos \left(E{O}_{n_m}{\theta }_2\right)\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ 1& \sin \left(E{O}_1{\theta }_{n_{\mathrm{btt}}}\right)& \cos \left(E{O}_1{\theta }_{n_{\mathrm{btt}}}\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_{n_{\mathrm{btt}}}\right)& \cos \left(E{O}_{n_m}{\theta }_{n_{\mathrm{btt}}}\right)\end{array}\right]$

blade is realized.

[0076] Further, in step 1), the three-dimensional finite element model of the rotor blade is established, first n.sub.m order modal parameters: a modal frequency f.sub.i, a displacement modal shape ϕ.sub.i with the size being n.sub.dof×1 and a strain modal shape ψ.sub.i with the size being 2 n.sub.dof×1 are extracted by modal analysis; a full-field strain modal shape matrix, ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] with the size being 2 n.sub.dof×n.sub.m, of the rotor blade is constructed; i represents for a modal order, n.sub.dof represents for the number of degrees of freedom of the finite element model of the blade, n.sub.dof=3 n.sub.n, and n.sub.n represents for the number of nodes of the finite element model of the blade; and the strains include six strain components in total, including three positive strains ε.sub.z and three shearing strains γ.sub.xz, namely each node has the six strain modal shapes.

[0077] Further, in step 2), a relationship between the number n.sub.btt of the blade tip timing sensors mounted in a circumferential direction of a rotor blade casing and the multimodal number n.sub.m of vibration is that: n.sub.btt≥2 n.sub.m+1; and a measuring point selection matrix S.sub.btt arranged in a circumferential direction of a blade tip timing sensor casing under multimodal excitation is constructed:

wherein the size of the measuring point selection matrix S.sub.btt is n.sub.btt×(2 n.sub.m+1), EO.sub.i represents for a concerned excitation order, and θ.sub.j represents for a mounting angle of the j.sup.th blade tip timing sensor in the casing; the position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated; and the random process is repeated for R times, and a measuring point layout solution when the condition number κ of the matrix is minimum is selected therefrom.

[0078] Further, in step 3), a conversion matrix,

[00011]$T=\left[\frac{{\psi }_1}{{\varphi }_{\mathrm{btt},1}},\frac{{\psi }_2}{{\varphi }_{\mathrm{btt},2}},\mathrm{.Math.},\frac{{\psi }_i}{{\varphi }_{\mathrm{btt},i}},\mathrm{.Math.},\frac{{\psi }_{n_m}}{{\varphi }_{\mathrm{btt},n_m}}\right]$

with the size being 2 n.sub.dof×n.sub.d, of the single-point displacement and the full-field dynamic strains of the blade is constructed, wherein ϕ.sub.btt,i represents for an i.sup.th order displacement modal shape of a blade tip timing measuring point.

[0079] Further, in step 4), blade tip multimodal vibration signals u.sub.btt(t) on the N.sup.th turn of the rotor blade is acquired by the n.sub.btt blade tip timing sensors, and furthermore, an n.sub.m order vibration parameter β=(S.sub.btt).sup.†u.sub.btt(t) is acquired by utilizing a circumferential Fourier fitting algorithm,

[0080] wherein β=[c,A.sub.1 cos(φ.sub.1),A.sub.1 sin(φ.sub.1), . . . ,A.sub.n.sub.m cos(φ.sub.n.sub.m),A.sub.n.sub.m sin (φ.sub.n.sub.m)].sup.T, A.sub.i represents for a vibration amplitude, f.sub.i represents for a modal frequency, and φ.sub.i represents for an initial phase; the superscript † represents for the inversion of the matrix; the superscript T represents for the transposition of a vector; and further, a blade tip single-point multimodal vibration signal is acquired by reconstruction:

[00012]$u_{\mathrm{btt}}\left(t\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{u}_{\mathrm{btt},i}\left(t\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{A}_i\sin \left({\omega }_{i}t+{\phi }_i\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{A}_i\cos \left({\phi }_i\right)\sin \left({\omega }_{i}t\right)+A_i\sin \left({\phi }_i\right)\cos \left({\omega }_{i}t\right)$

wherein u.sub.btt,i(t) represents for decoupled i.sup.th order vibration signals, and c represents for static deformation of the blade.

[0081] Further, in step 5), the strains of all the nodes on the surface of and inside the blade at each turn of the rotor blade are calculated according to a formula S(t)=T[u.sub.btt,1,u.sub.btt,2, . . . ,u.sub.btt,i, . . .,u.sub.btt,n.sub.m].sup.T to obtain S(t)=[ε.sub.1,x,ε.sub.1,y,ε.sub.1,z,γ.sub.1,xy,γ.sub.1,yz,γ.sub.1,xz, . . . ,ε.sub.i,x,ε.sub.i,y,ε.sub.i,z,γ.sub.i,xy,γ.sub.i,yz,γ.sub.i,xz, . . . ,ε.sub.n.sub.n.sub.,x,ε.sub.n.sub.n.sub.,y,ε.sub.n.sub.n.sub.,z,γ.sub.n.sub.n.sub.,xy,γ.sub.n.sub.n.sub.,yz,γ.sub.n.sub.n.sub.,xz].sup.T, wherein the strains S(t) include the positive strains and the shearing strains.

[0082] On the other hand, a measuring system for implementing the method includes:

[0083] a plurality of blade tip timing sensors 1 arranged on a rotor blade casing 2;

[0084] a blade tip timing vibration measuring module connected with the blade tip timing sensors so as to measure circumferential blade tip multimodal vibration signals of a rotor blade 4; and

[0085] a calculation unit connected with the blade tip timing vibration measuring module, and the calculation unit including:

[0086] a modal analysis module configured to perform modal analysis based on a three-dimensional finite element model of a to-be-measured rotor blade so as to acquire a displacement modal shape φ.sub.i of first n.sub.m order modal frequency f.sub.i, of the rotor blade and construct a full-field strain modal shape matrix ψ=[ψ.sub.1ψ.sub.2, . . . ,ψ.sub.n.sub.m] of the rotor blade;

[0087] a measuring point optimization module configured to optimize the number of measuring points of the blade tip timing sensors arranged on the rotor blade, wherein a measuring point selection matrix S.sub.btt arranged in the circumferential direction of a blade tip timing sensor casing under multimodal excitation is constructed, the position where a rotating speed sensor is located serves as a reference 0°, a limited mounting angle range of the casing is removed, and n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as mounting positions of the blade tip timing sensors; and the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated, a random process is repeated for R times, and a measuring point layout when the condition number κ of the matrix is minimum is selected therefrom;

[0088] a conversion matrix calculation module configured to construct a conversion matrix of single-point displacement and full-field dynamic strains of the blade; and

[0089] a strain field reconstruction module configured to calculate the positive strains and shearing strains of all nodes on the surface of and inside the blade at each turn of the rotor blade according according to a formula S(t)=T[u.sub.btt,1,u.sub.btt,2, . . . ,u.sub.btt,i, . . .,u.sub.btt,n.sub.m].sup.T to obtain S(t)=[ε.sub.1,x,ε.sub.1,y,ε.sub.1,z,γ.sub.1,xy,γ.sub.1,yz,γ.sub.1,xz, . . . ,ε.sub.i,x,ε.sub.i,y,ε.sub.i,z,γ.sub.i,xy,γ.sub.i,yz,γ.sub.i,xz, . . . ,ε.sub.n.sub.n.sub.,x,ε.sub.n.sub.n.sub.,y,ε.sub.n.sub.n.sub.,z,γ.sub.n.sub.n.sub.,xy,γ.sub.n.sub.n.sub.,yz,γ.sub.n.sub.n.sub.,xz].sup.T, wherein the strains S(t) include the positive strains and the shearing strains.

[0090] In one embodiment, the blade tip timing vibration measuring module includes the rotating speed sensor, a signal conditioning module and a time-displacement conversion module.

[0091] In one embodiment, the measuring system further includes a display unit and a wireless communication device, and the wireless communication device includes a 4G/GPRS or an internet communication module.

[0092] In one embodiment, the modal analysis module, the measuring point optimization module, the conversion matrix calculation module or a dynamic strain field reconstruction module is a general-purpose processor, a digital signal processor, an application-specific integrated circuit (ASIC) or a field-programmable gate array (FPGA).

[0093] In one embodiment, the modal analysis module, the measuring point optimization module, the conversion matrix calculation module or the dynamic strain field reconstruction module includes a memory, and the memory includes one or more read-only memories (ROMs), random access memories (RAMs), flash memories or electronic erasable programmable read only memories (EEPROMs).

[0094] In one embodiment, a system for the dynamic strain field measuring method for the rotor blade based on blade tip timing includes:

[0095] on the other hand, the present invention further provides the system for the above-mentioned dynamic strain field measuring method for the rotor blade based on blade tip timing, including:

[0096] the modal analysis module configured to perform modal analysis on the three-dimensional finite element model of the blade by utilizing finite element analysis software to extract the first n.sub.m order modal parameters: the modal frequency f.sub.i, the displacement modal shape ϕ.sub.i with the size being n.sub.dof×1 and a strain modal shape ψ.sub.i with the size being 2 n.sub.dof×1; and to construct the full-field strain modal matrix, ψ=[ψ.sub.1,ψ.sub.2, . . . ,ψ.sub.n.sub.m] with the size being 2 n.sub.dof×n.sub.m, of the rotor blade, wherein i represents for the modal order, n.sub.dof represents for the number of degrees of freedom of the finite element model of the blade, n.sub.dof=3 n.sub.n, n.sub.n represents for the number of the nodes of the finite element model of the blade, and the strains include the six strain components in total, including the three positive strains ε.sub.z and the three shearing strains γ.sub.xz, namely each node has the six strain modal shapes;

[0097] The measuring point optimization module configured to optimize the relationship, n.sub.btt≥2 n.sub.m+1, between the number n.sub.btt of the blade tip timing sensors mounted in the circumferential direction of the rotor blade casing and the multimodal number n.sub.m of vibration; and to construct the measuring point selection matrix S.sub.btt arranged in the circumferential direction of the blade tip timing sensor casing under multimodal excitation:

[00013]$S_{\mathrm{btt}}=\left[\begin{array}{cccccc}1& \sin \left(E{O}_1{\theta }_1\right)& \cos \left(E{O}_1{\theta }_1\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_1\right)& \cos \left(E{O}_{n_m}{\theta }_1\right)\\ 1& \sin \left(E{O}_1{\theta }_2\right)& \cos \left(E{O}_1{\theta }_2\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_2\right)& \cos \left(E{O}_{n_m}{\theta }_2\right)\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ 1& \sin \left(E{O}_1{\theta }_{n_{\mathrm{btt}}}\right)& \cos \left(E{O}_1{\theta }_{n_{\mathrm{btt}}}\right)& \mathrm{.Math.}& \sin \left(E{O}_{n_m}{\theta }_{n_{\mathrm{btt}}}\right)& \cos \left(E{O}_{n_m}{\theta }_{n_{\mathrm{btt}}}\right)\end{array}\right]$

wherein the size of the measuring point selection matrix S.sub.btt is n.sub.btt×(2 n.sub.m+1), θ.sub.j represents for the layout angles of the blade tip timing sensors j(j=1, . . . n.sub.btt).sub.th the casing, n.sub.btt represents for the total number of the j(j=1, . . . n.sub.btt).sub.th blade tip timing sensors arranged in the circumferential direction, EO.sub.i represents for the concerned excitation order (i=1, . . . n.sub.m), and θ.sub.j represents for the mounting angle of the i.sup.th blade tip timing sensor in the casing; the position where the rotating speed sensor is located serves as the reference 0°, the limited mounting angle range of the casing is removed, and the n.sub.btt angles are randomly selected in the circumferential direction of the casing to serve as the mounting positions of the blade tip timing sensors; the condition number κ of the blade tip timing measuring point selection matrix S.sub.btt is calculated; and the random process is repeated for R times, and the measuring point layout solution when the condition number κ of the matrix is minimum is selected therefrom;

[0098] the conversion matrix calculation module configured to construct the conversion matrix,

[00014]$T=\left[\frac{{\psi }_1}{{\varphi }_{\mathrm{btt},1}},\frac{{\psi }_2}{{\varphi }_{\mathrm{btt},2}},\mathrm{.Math.},\frac{{\psi }_i}{{\varphi }_{\mathrm{btt},i}},\mathrm{.Math.},\frac{{\psi }_{n_m}}{{\varphi }_{\mathrm{btt},n_m}}\right]$

with the size being 2 n.sub.dof×n.sub.d, of the single-point displacement and the full-field dynamic strains of the blade, wherein ϕ.sub.btt,i represents for the i.sup.th order displacement modal shape of the blade tip timing measuring point;

[0099] the blade tip timing vibration measuring module including the plurality of blade tip timing sensors, the at least one rotating speed sensor, the signal conditioning module and the time-displacement conversion module; and being configured to acquire the blade tip multimodal vibration signals u.sub.btt(t) on the N.sup.th turn of the rotor blade by n.sub.btt blade tip timing sensors; to further acquire the n.sub.m order vibration parameter by β=(S.sub.btt).sup.†u.sub.btt(t) by utilizing the circumferential Fourier fitting algorithm,

[0100] wherein β=[c,A.sub.1 cos(φ.sub.1),A.sub.1 sin(φ.sub.1), . . . ,A.sub.n.sub.m cos(φ.sub.n.sub.m),A.sub.n.sub.m sin (φ.sub.n.sub.m)].sup.T, A.sub.i represents for the vibration amplitude, f.sub.i represents for the modal frequency, and φ.sub.i represents for the initial phase; and to further acquire the blade tip single-point multimodal vibration signals by reconstruction:

[00015]$u_{\mathrm{btt}}\left(t\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{u}_{\mathrm{btt},i}\left(t\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{A}_i\sin \left({\omega }_{i}t+{\phi }_i\right)=c+\underset{i=1}{\overset{n_m}{.Math.}}{A}_i\cos \left({\phi }_i\right)\sin \left({\omega }_{i}t\right)+A_i\sin \left({\phi }_i\right)\cos \left({\omega }_{i}t\right)$

wherein u.sub.btt,i(t) represents for the decoupled i.sup.th order vibration signals, and c represents for static deformation of the blade; and

[0101] the strain field reconstruction module configured to calculate the strains of all the nodes on the surface of and inside the blade at each turn of the rotor blade according to the formula