CN117648779B - Method, device, equipment and computer storage medium for designing camber line of blade - Google Patents
Method, device, equipment and computer storage medium for designing camber line of blade Download PDFInfo
- Publication number
- CN117648779B CN117648779B CN202410123620.8A CN202410123620A CN117648779B CN 117648779 B CN117648779 B CN 117648779B CN 202410123620 A CN202410123620 A CN 202410123620A CN 117648779 B CN117648779 B CN 117648779B
- Authority
- CN
- China
- Prior art keywords
- coordinate
- camber line
- bezier curve
- points
- control points
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 86
- 238000003860 storage Methods 0.000 title claims abstract description 14
- 238000005070 sampling Methods 0.000 claims abstract description 101
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 41
- 238000013461 design Methods 0.000 claims abstract description 36
- 238000005452 bending Methods 0.000 claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims description 59
- 238000012545 processing Methods 0.000 claims description 11
- 238000003491 array Methods 0.000 claims description 10
- 238000004891 communication Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 230000001360 synchronised effect Effects 0.000 description 4
- 230000003068 static effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- KLDZYURQCUYZBL-UHFFFAOYSA-N 2-[3-[(2-hydroxyphenyl)methylideneamino]propyliminomethyl]phenol Chemical compound OC1=CC=CC=C1C=NCCCN=CC1=CC=CC=C1O KLDZYURQCUYZBL-UHFFFAOYSA-N 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 201000001098 delayed sleep phase syndrome Diseases 0.000 description 1
- 208000033921 delayed sleep phase type circadian rhythm sleep disease Diseases 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Landscapes
- Image Generation (AREA)
Abstract
The embodiment of the disclosure provides a design method, a device, equipment and a computer storage medium of a camber line of a blade; the design method comprises the following steps: collecting a plurality of sampling points on an initial camber line corresponding to the uniform section of the blade; acquiring second coordinates of a plurality of control points of the Bezier curve corresponding to the initial mean camber line by using a singular value decomposition method based on the acquired first coordinates corresponding to the plurality of sampling points; adjusting the bending degree of the Bezier curve corresponding to the initial mean camber line by adjusting the second coordinates of the control points; the Bezier curve after the bending is adjusted is used for representing a final camber line corresponding to the constant section of the blade.
Description
Technical Field
The embodiment of the disclosure relates to the technical field of design of axial flow impellers, in particular to a method, a device and equipment for designing camber lines of blades and a computer storage medium.
Background
The axial flow impeller machine is a power machine widely applied to the fields of aviation, ships, electric power, metallurgy, energy, chemical industry, medicine and the like, and is one of core equipment of many large-scale industrial production enterprises. However, because of the long period and high research and development cost of the structural design of the axial flow impeller machine, the development and market popularization of new products of the axial flow impeller machine are severely restricted.
In the blade design of an axial flow impeller machine, a plurality of equal sections are generally selected, parameters of the selected equal sections are respectively adjusted, and the equal sections are stacked in a specified manner, so as to obtain a geometric shape (hereinafter referred to as a "blade shape") of the blade. And then, simulating pneumatic performance and mechanical performance of the obtained blade to preliminarily judge whether the blade meets the design requirement. If the blade does not meet the design requirement, continuing to adjust the parameters of the plurality of equal sections, and then continuing to simulate, wherein the blade profile meeting the design requirement is expected to be obtained by adjusting the parameters of the plurality of equal sections for a plurality of times. In the process of adjusting the parameters of a plurality of uniform sections, the shape of the camber line is the basis for designing the blade profile, and plays a decisive role in the aerodynamic performance and the mechanical performance of the blade. By adjusting the shape of the camber line, parameters of the uniform section in the blade can be quickly adjusted, so that the aerodynamic performance and the mechanical performance of the blade are changed.
Currently, in the related art, a Bezier curve is mainly obtained by adopting least square fitting to represent a mean camber line. However, the least square method requires that as many sampling points as possible are acquired in the fitting process so as to ensure the accuracy of Bezier curve fitting, and the number of the sampling points increases the calculated amount to a certain extent. On the other hand, the least square method involves the inversion calculation of the matrix in the calculation process, but the inversion calculation of the matrix with high dimension takes a long time, resulting in low fitting efficiency of the bezier curve and low efficiency of adjusting the parameters of the above-mentioned constant cross section.
Disclosure of Invention
The embodiment of the disclosure provides a design method, a device, equipment and a computer storage medium of a camber line of a blade; the fitting precision of the camber line of the blade and the adjusting efficiency of the camber line shape can be improved, and the design efficiency of the blade is further improved.
The technical scheme of the embodiment of the disclosure is realized as follows:
In a first aspect, an embodiment of the present disclosure provides a method for designing a camber line of a blade, the method comprising:
collecting a plurality of sampling points on an initial camber line corresponding to the uniform section of the blade;
Acquiring second coordinates of a plurality of control points of the Bezier curve corresponding to the initial mean camber line by using a singular value decomposition method based on the acquired first coordinates corresponding to the plurality of sampling points;
Adjusting the bending degree of the Bezier curve corresponding to the initial mean camber line by adjusting the second coordinates of the control points; the Bezier curve after the bending is adjusted is used for representing a final camber line corresponding to the constant section of the blade.
In a second aspect, embodiments of the present disclosure provide a design apparatus for a camber line of a blade, the design apparatus comprising: a collection unit, an acquisition unit, and an adjustment unit; wherein,
The acquisition section is configured to: collecting a plurality of sampling points on an initial camber line corresponding to the uniform section of the blade;
the acquisition section is configured to: acquiring second coordinates of a plurality of control points of the Bezier curve corresponding to the initial mean camber line by using a singular value decomposition method based on the acquired first coordinates corresponding to the plurality of sampling points;
the adjustment section is configured to: adjusting the bending degree of the Bezier curve corresponding to the initial mean camber line by adjusting the second coordinates of the control points; the Bezier curve after the bending is adjusted is used for representing a final camber line corresponding to the constant section of the blade.
In a third aspect, the disclosed embodiments provide a computing device comprising: a communication interface, a memory and a processor; the components are coupled together by a bus system; wherein,
The communication interface is used for receiving and transmitting signals in the process of receiving and transmitting information with other external network elements;
The memory is used for storing a computer program capable of running on the processor;
The processor is configured to execute the steps of the method for designing a camber line of a blade according to the first aspect when the computer program is run.
In a fourth aspect, embodiments of the present disclosure provide a computer storage medium storing a design program of a camber line of a blade, which when executed by at least one processor, implements the steps of the method of designing a camber line of a blade according to the first aspect.
The embodiment of the disclosure provides a design method, a device, equipment and a computer storage medium of a camber line of a blade; and acquiring a plurality of sampling points on an initial camber line corresponding to the uniform section of the blade, and further obtaining second coordinates of a plurality of control points of the Bezier curve corresponding to the initial camber line by utilizing a singular value decomposition method based on first coordinates corresponding to the acquired plurality of sampling points, so as to design and obtain a final camber line by adjusting the bending degree of the Bezier curve through adjusting the second coordinates of the plurality of control points. The Bezier curve obtained by fitting through the design method is high in accuracy and high in calculation efficiency, the design efficiency of the equal section of the blade is greatly improved, and the design period of the blade is shortened.
Drawings
Fig. 1 is a flow chart of a design method of camber line of a blade according to an embodiment of the disclosure.
Fig. 2 is a schematic structural view of a blade with a uniform cross section according to an embodiment of the present disclosure.
Fig. 3 is a schematic diagram of the position of a control point of a cubic bezier curve.
Fig. 4 is a schematic diagram of a position of a control point corresponding to a bezier curve according to an embodiment of the disclosure.
Fig. 5 is a schematic diagram illustrating a design apparatus for camber line of a blade according to an embodiment of the disclosure.
Fig. 6 is a schematic diagram of a specific hardware structure of a computing device according to an embodiment of the disclosure.
Detailed Description
The technical solutions in the embodiments of the present disclosure will be clearly and completely described below with reference to the drawings in the embodiments of the present disclosure.
Referring to fig. 1, a method for designing a camber line of a blade according to an embodiment of the present disclosure is shown, and the method specifically includes the following steps.
In step S101, a plurality of sampling points are acquired on an initial camber line corresponding to a uniform section of the blade.
As shown in fig. 2, which shows a schematic structural view of a uniform cross section of the blade. As can be seen from fig. 2, the constant section of the blade is constituted by a pressure surface PS, a trailing edge line TE, a suction surface SS and a leading edge line LE. The camber line of the blade is a curve obtained by making a plurality of inscribed circles on the blade and then connecting the circle centers of the inscribed circles.
In general, during the design process of the blade, the blade profile corresponding to the constant cross section of the blade can be approximately obtained through the preliminary design, and then the camber line corresponding to the constant cross section of the current blade is determined, such as the curve camber in fig. 2. In the embodiment of the present disclosure, the camber line obtained from the blade profile corresponding to the substantially uniform cross section of the blade is referred to as an initial camber line. In order to design a blade with aerodynamic and mechanical properties that meet design requirements, the shape of the initial camber line also needs to be adjusted during implementation.
It should be noted that, the method for obtaining the initial camber line corresponding to the uniform section of the blade is not specifically limited in the embodiments of the present disclosure.
In some examples, the number of sampling points may be determined according to actual requirements. Optionally, the number of the sampling points is greater than 200, so as to improve the fitting precision of the Bezier curve. It should be noted that, the sampling points do not include two end points of the initial mean camber line.
In step S102, based on the acquired first coordinates corresponding to the plurality of sampling points, second coordinates of a plurality of control points of the bezier curve corresponding to the initial mean camber line are acquired by using a singular value decomposition method.
In the embodiment of the disclosure, based on the first coordinates corresponding to the collected multiple sampling points, a singular value decomposition (Singular Value Decomposition, SVD) method is used to perform eigenvalue decomposition, so as to obtain the second coordinates of multiple control points of the bezier curve corresponding to the initial mean camber line. The SVD method is adopted in the fitting method of the Bezier curve, and the SVD method can be used for carrying out dimension reduction treatment on the matrix in the fitting process, so that the problem of long time consumption caused by matrix inversion operation when the Bezier curve is fitted by using a least square method in the related art is avoided, and the calculation efficiency is improved.
It will be appreciated that the Bezier curve is defined by a plurality of control points P 0、P1、…、Pn, where n represents the order of the Bezier curve. Typically control point P 0 and control point P n are located at the two endpoints of the initial mean camber line described above. In the disclosed embodiment, control point P 0 is defined as an initial endpoint of the two endpoints of the initial mean camber line, and control point P n is defined as an end endpoint of the two endpoints of the initial mean camber line.
In general, a plurality of control points of the bezier curve are used to control the shape of the bezier curve, and the shape of the entire bezier curve can be adjusted by adjusting the positions of the control points. The bezier curve can be represented by the following parameterized equation:
the parameterized equations above are written in a matrix as follows:
wherein, Represents the basis of the Bezier curve at the parameter t, and/>; P j represents the j-th control point of the bezier curve, and the coordinates of P j are (x Pj,yPj); n represents the order of the above Bezier curve.
Illustratively, fig. 3 shows a cubic bezier curve. The cubic bezier curve corresponds to four control points (P 0,P1,P2,P3). In a specific implementation process, the four control points are sequentially connected into three line segments (not connected end to end), three new interpolation points (shown as black squares in fig. 3) are obtained by performing linear interpolation on the three line segments based on the parameter t, the three interpolation points can be sequentially connected into two line segments, then the two line segments are subjected to linear interpolation by using the parameter t, two new interpolation points (shown as black triangles in fig. 3) are obtained, and interpolation is performed on the line segments formed by connecting the two new interpolation points by using the parameter t, so that the control point of the three Bezier curves is finally obtained.
It will be appreciated that in the embodiments of the present disclosure, the second coordinates of the plurality of control points of the Bezier curve are solved, that is, the array in the parameterized equation corresponding to the Bezier curve is solved。
In some examples, the first coordinates corresponding to the plurality of sampling points are coordinates in a cartesian rectangular coordinate system.
In step S103, the degree of curvature of the bezier curve corresponding to the initial mean camber line is adjusted by adjusting the second coordinates of the plurality of control points; the Bezier curve after the bending is adjusted is used for representing a final camber line corresponding to the constant section of the blade.
It can be understood that the plurality of control points are used for controlling the shape of the bezier curve, so that in a specific implementation process, the curvature of the bezier curve can be adjusted by adjusting the second coordinates of the plurality of control points, and the bezier curve after the curvature is adjusted, namely, a final camber line corresponding to the constant cross section of the vane is represented. In addition, when the second coordinates of a plurality of control points of the Bezier curve are known, fewer adjustment times can be adopted to obtain a final mean camber line as much as possible in the implementation process, so that the negative influence caused by repeatedly adjusting the Bezier curve for many times is eliminated, and the design efficiency is improved.
For the technical scheme shown in fig. 1, a plurality of sampling points are collected on an initial camber line corresponding to a uniform section of a blade, and then based on first coordinates corresponding to the collected plurality of sampling points, second coordinates of a plurality of control points of a bezier curve corresponding to the initial camber line are obtained by utilizing an SVD method, so that the bending degree of the bezier curve is adjusted by adjusting the second coordinates of the plurality of control points, and a final camber line is obtained. The Bezier curve obtained by fitting through the design method is high in accuracy and high in calculation efficiency, the design efficiency of the equal section of the blade is greatly improved, and the design period of the blade is shortened.
For the technical solution shown in fig. 1, in some possible embodiments, the acquiring, based on the first coordinates corresponding to the acquired plurality of sampling points, the second coordinates of a plurality of control points of the bezier curve corresponding to the initial mean camber line by using a singular value decomposition method includes:
acquiring a first coordinate array corresponding to the plurality of sampling points based on the first coordinates;
Based on the first coordinate arrays corresponding to the plurality of sampling points, normalizing to obtain arc length parameter arrays corresponding to the plurality of sampling points;
Based on a third coordinate corresponding to an initial endpoint of the two endpoints of the initial mean camber line and a fourth coordinate corresponding to a terminal endpoint of the two endpoints, respectively acquiring slopes corresponding to an abscissa and an ordinate of the third coordinate and the fourth coordinate;
calculating a characteristic value decomposition matrix corresponding to the Bezier curve with a set order by utilizing a singular value decomposition method based on the arc length parameter array;
Calculating a second coordinate array corresponding to a plurality of control points of the Bezier curve based on the eigenvalue decomposition matrix; wherein the second coordinate array is used for representing the second coordinates.
It will be appreciated that after N sampling points are acquired on the initial mean camber line, the abscissa and the ordinate of the first coordinate corresponding to each of the N sampling points are stored in the first coordinate array [ [ X, Y ] ] N×1.
When the third coordinate of the initial endpoint and the fourth coordinate of the terminal endpoint in the two endpoints of the initial mean camber line are known, the slopes of the horizontal coordinate and the vertical coordinate in the third coordinate and the fourth coordinate can be calculated respectively. Specifically, when the third coordinates corresponding to the two end points of the initial mean camber line are A1 (x A1,yA1) and A2 (x A2,yA2), respectively, and the corresponding derivatives are obtained by deriving in the x direction and the y direction according to the curve equation corresponding to the initial mean camber line, respectively, the slopes of the horizontal coordinate and the vertical coordinate in the third coordinates and the fourth coordinates can be obtained by setting x=x A1,x=xA2 and y=y A1,y=yA2.
For the above embodiment, in some examples, the normalizing processing, based on the first coordinate arrays corresponding to the plurality of sampling points, obtains arc length parameter arrays corresponding to the plurality of sampling points, including:
Calculating a distance L i between the ith-1 sampling point and the ith sampling point based on a first coordinate corresponding to the ith-1 sampling point in the first coordinate array and a first coordinate corresponding to the ith sampling point in the first coordinate array; wherein i is more than or equal to 1 and less than or equal to N, N represents the number of sampling points;
Based on the distance L i between the i-1 th sampling point and the i-th sampling point, calculating according to t i'=ti-1+Li to obtain an initial arc length parameter t i' corresponding to the i-th sampling point; wherein t i-1 represents the arc length parameter after normalization processing corresponding to the i-1 th sampling point;
According to Obtaining the arc length parameter after normalization processing corresponding to the ith sampling point; wherein,Representing an initial arc length parameter corresponding to an Nth sampling point;
And obtaining an arc length parameter array corresponding to the plurality of sampling points according to the normalized arc length parameters corresponding to all the sampling points.
It will be appreciated that the distance L i between the i-1 th sample point and the i-th sample point may be based onAs a result, the embodiments of the present disclosure will not be described in detail.
It should be noted that, in the embodiment of the present disclosure, performing normalization processing on the arc length parameter of each sampling point can facilitate data processing in the subsequent technical solutions.
For the above embodiment, in some examples, the calculating, based on the arc length parameter array, a eigenvalue decomposition matrix corresponding to the bezier curve with a set order by using a singular value decomposition method includes:
Acquiring the number of control points of the Bezier curve, the column number and the line number of a coefficient matrix corresponding to the first coordinate based on the acquired number of sampling points and the order of the Bezier curve;
The basis corresponding to the different control points at each sampling point is obtained according to the following equation:
wherein, Representing the normalized arc length parameter/>, corresponding to the ith sampling pointThe base corresponding to the j-th control point below is more than or equal to 0 and less than or equal to n, wherein n represents the order of the Bessel curve;
obtaining the coefficient matrix based on the bases corresponding to different control points at each sampling point;
and calculating the eigenvalue decomposition matrix by using a singular value decomposition method based on the coefficient matrix.
Optionally, the acquiring, based on the number of collected sampling points and the order of the bezier curve, the number of control points corresponding to the bezier curve, the number of columns and the number of rows of the coefficient matrix corresponding to the first coordinate includes:
The number of control points of the Bezier curve is equal to the order +1 of the Bezier curve;
The number of columns of the coefficient matrix is equal to 2 times of the number of control points of the Bezier curve, and is reduced by 6;
the number of lines of the coefficient matrix is equal to 2 times the number of the sampling points.
For the above embodiment, in some examples, obtaining the coefficient matrix based on the base corresponding to the different control point at each sampling point includes:
based on the basis corresponding to different control points at each sampling point, the obtained coefficient matrix M is:
。
specifically, after the bases corresponding to different control points at each sampling point are obtained by solving one by one, the matrix form of the parameterized equation based on the bezier curve is as follows:
in the embodiment of the present disclosure, the coefficient matrix m=as described above 。
In addition, when the above coefficient matrix M is decomposed into eigenvalues by SVD, the coefficient matrix M is a matrix of nx (n+1), and is decomposed into products of three specific matrices, specifically as follows:
M=R×W×VT
Wherein, R is an N×N orthogonal matrix, and R is the eigenvector of MM T in the implementation process; w is a diagonal matrix of N x (n+1), and in the implementation is the non-negative square root of M T M or MM T, and is arranged in descending order; v is an orthogonal matrix of (n+1) × (n+1), and in the implementation is a eigenvector of M T M. In the embodiment of the present disclosure, the step of performing eigenvalue decomposition on the coefficient matrix M by using the SVD method is not described in detail.
The above-mentioned SVD method is used to decompose the characteristic value of the coefficient matrix, so that the characteristic of the sampling point can be largely maintained, and the dimension of the matrix operation can be effectively reduced, thereby improving the calculation efficiency.
Furthermore, in some examples, diagonal matrix W in the feature decomposition matrix obtained by embodiments of the present disclosure requires noise values to be removed during implementation to avoid affecting subsequent fitting results. The specific rejecting method comprises the following steps: the data minimum in the diagonal matrix W is set, for example, to a precision of 1e-6. In addition, since the diagonal matrix W is a matrix having a rank much smaller than that of the coefficient matrix M, the technical solution provided by the embodiment of the disclosure greatly reduces the calculation amount in the calculation process.
For the above embodiment, in some examples, the calculating, based on the eigenvalue decomposition matrix, a second coordinate array corresponding to a plurality of control points of the bezier curve includes:
Acquiring numerical contribution arrays of the two endpoints and control points adjacent to the two endpoints to all sampling points based on the base corresponding to different control points at each sampling point, the slopes corresponding to the third coordinate and the abscissa and the ordinate in the third coordinate, and the slopes corresponding to the fourth coordinate and the abscissa and the ordinate in the fourth coordinate; wherein the numerical contribution array is expressed as And (2) and，；/>An abscissa indicating the third coordinate; a1 represents a slope corresponding to an abscissa in the third coordinates; /(I)An abscissa in the fourth coordinates; a2 represents a slope corresponding to an abscissa in the fourth coordinate; /(I)Representing the ordinate of the third coordinates; a3 represents a slope corresponding to an ordinate in the third coordinate; /(I)Representing the ordinate in the fourth coordinate; a4 represents a slope corresponding to an ordinate in the fourth coordinate;
Calculating a third coordinate array corresponding to the plurality of sampling points based on the first coordinate array and the numerical contribution array;
Calculating to obtain a constant array based on the first coordinate array and the third coordinate array;
and obtaining a second coordinate array corresponding to the plurality of control points based on the eigenvalue decomposition matrix and the constant array.
Specifically, referring to fig. 3, the initial end point of the two end points of the initial mean camber line is P 0, the control point adjacent to the initial end point P 0 is P 1, the end point is P 3, and the control point adjacent to the end point P 3 is P 2, so that the numerical contributions of the two end points and the control points P 1 and P 2 adjacent to the two end points to each sampling point can be obtained, and the array formed by the numerical contributions corresponding to all the sampling points is called as a numerical contribution array [ [ X ', Y' ] ] ] N×1.
In some examples, subtracting the value contribution array [ [ X ', Y' ] ] N×1 from the first coordinate array [ [ X, Y ] ] N×1 described above yields a third coordinate array [ [ X ", Y" ] ] N×1.
Finally, subtracting the third coordinate array [ [ X ', Y' ] ] N×1 from the first coordinate array [ [ X, Y ] ] N×1 to obtain a constant array。
It will be appreciated that the constant arrayTo eliminate the array corresponding to the abscissa and the ordinate of the weight occupied by the two end points on the initial middle arc line and the control points adjacent to the two end points.
In some examples, after the eigenvalue matrix U, W, V and the constant array are obtained by using the SVD method, a second coordinate array of the plurality of control points can be calculated according to the following equation, specifically:
Wherein b represents the above-mentioned constant array.
The method further comprises the following steps:。
From this, from the solution A second coordinate of the plurality of control points can be determined.
Referring to fig. 4, a schematic diagram of the location of the control points (shown as black circles) in the bezier curve is shown. In the specific implementation process, the second coordinates of the control points are adjusted, so that the Bezier curve can be quickly adjusted, a final camber line corresponding to the uniform section of the blade is obtained, and the uniform section of the blade is adjusted. Through the technical scheme provided by the embodiment of the disclosure, not only the design precision of the camber line can be improved, but also the design efficiency of the blade in the embodiment of the disclosure is improved by 30% compared with that of the blade in the related art.
Based on the same inventive concept as the previous technical solution, referring to fig. 5, there is shown a design device 50 for camber line of a blade according to an embodiment of the present disclosure, where the design device 50 includes: an acquisition unit 501, an acquisition unit 502, and an adjustment unit 503; wherein,
The acquisition unit 501 is configured to: collecting a plurality of sampling points on an initial camber line corresponding to the uniform section of the blade;
The acquisition unit 502 is configured to: acquiring second coordinates of a plurality of control points of the Bezier curve corresponding to the initial mean camber line by using a singular value decomposition method based on the acquired first coordinates corresponding to the plurality of sampling points;
The adjusting unit 503 is configured to: adjusting the bending degree of the Bezier curve corresponding to the initial mean camber line by adjusting the second coordinates of the plurality of control points; the Bezier curve after the bending is adjusted is used for representing a final camber line corresponding to the constant section of the blade.
It should be noted that, in the design apparatus 50 for a camber line of a blade according to the foregoing embodiment, when implementing the function thereof, only the division of the foregoing functional modules is illustrated, and in practical application, the foregoing functional allocation may be implemented by different functional modules, that is, the internal structure of the terminal may be divided into different functional modules, so as to implement all or part of the functions described above. In addition, the device 50 for designing the camber line of the blade provided in the above embodiment belongs to the same concept as the embodiment of the method for designing the camber line of the blade, and the detailed implementation process of the device is shown in the method embodiment, which is not repeated here.
The components in this embodiment may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional modules.
The above-described integrated units, if implemented in the form of software functional modules, may be stored in a computer-readable storage medium, if not sold or used as separate products, and based on such understanding, the technical solution of the present embodiment may be embodied essentially or partly in the form of a software product, or all or part of the technical solution may be embodied in a storage medium, where the computer software product includes several instructions to cause a computer device (which may be a personal computer, a server, or a network device, etc.) or processor to perform all or part of the steps of the above-described method of the present embodiment. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read Only Memory (ROM), a random access memory (Random Access Memory, RAM), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
Accordingly, an embodiment of the present disclosure provides a computer storage medium storing a program for designing a camber line of a blade, which when executed by at least one processor, implements steps of a method for designing a camber line of a blade.
According to the design apparatus 50 of the camber line of the blade and the computer storage medium, referring to fig. 6, a specific hardware structure of a computing device 60 capable of implementing the design apparatus 50 of the camber line of the blade is shown in an embodiment of the disclosure, the computing device 60 may be a wireless device, a mobile or cellular phone (including a so-called smart phone), a Personal Digital Assistant (PDA), a video game console (including a video display, a mobile video game device, a mobile video conference unit), a laptop computer, a desktop computer, a television set-top box, a tablet computing device, an e-book reader, a fixed or mobile media player, etc. The computing device 60 includes: a communication interface 601, a memory 602 and a processor 603; the various components are coupled together by a bus system 604. It is understood that the bus system 604 is used to enable connected communications between these components. The bus system 604 includes a power bus, a control bus, and a status signal bus in addition to the data bus. But for clarity of illustration, the various buses are labeled as bus system 604 in fig. 6. Wherein,
The communication interface 601 is configured to receive and send signals during the process of receiving and sending information with other external network elements;
the memory 602 is configured to store a computer program that can be executed by the processor 603;
The processor 603 is configured to execute the following steps when executing the computer program:
collecting a plurality of sampling points on an initial camber line corresponding to the uniform section of the blade;
Acquiring second coordinates of a plurality of control points of the Bezier curve corresponding to the initial mean camber line by using a singular value decomposition method based on the acquired first coordinates corresponding to the plurality of sampling points;
Adjusting the bending degree of the Bezier curve corresponding to the initial mean camber line by adjusting the second coordinates of the plurality of control points; the Bezier curve after the bending is adjusted is used for representing a final camber line corresponding to the constant section of the blade.
It is to be appreciated that the memory 602 in embodiments of the disclosure may be either volatile memory or nonvolatile memory, or may include both volatile and nonvolatile memory. The nonvolatile Memory may be a Read-Only Memory (ROM), a Programmable ROM (PROM), an Erasable PROM (EPROM), an Electrically Erasable EPROM (EEPROM), or a flash Memory. The volatile memory may be random access memory (Random Access Memory, RAM) which acts as external cache memory. By way of example, and not limitation, many forms of RAM are available, such as static random access memory (STATIC RAM, SRAM), dynamic random access memory (DYNAMIC RAM, DRAM), synchronous Dynamic Random Access Memory (SDRAM), double data rate Synchronous dynamic random access memory (Double DATA RATE SDRAM, DDRSDRAM), enhanced Synchronous dynamic random access memory (ENHANCED SDRAM, ESDRAM), synchronous link dynamic random access memory (SYNCHLINK DRAM, SLDRAM), and Direct memory bus random access memory (DRRAM). The memory 602 of the systems and methods described herein is intended to comprise, without being limited to, these and any other suitable types of memory.
And the processor 603 may be an integrated circuit chip with signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuitry of hardware in the processor 603 or instructions in the form of software. The processor 603 may be a general purpose processor, a digital signal processor (DIGITAL SIGNAL processor, DSP), an Application SPECIFIC INTEGRATED Circuit (ASIC), a field programmable gate array (Field Programmable GATE ARRAY, FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware components. The various methods, steps and logic blocks of the disclosure in the embodiments of the disclosure may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of a method disclosed in connection with the embodiments of the present disclosure may be embodied directly in hardware, in a decoded processor, or in a combination of hardware and software modules in a decoded processor. The software modules may be located in a random access memory, flash memory, read only memory, programmable read only memory, or electrically erasable programmable memory, registers, etc. as well known in the art. The storage medium is located in the memory 602, and the processor 603 reads information in the memory 602, and in combination with its hardware, performs the steps of the method described above.
It is to be understood that the embodiments described herein may be implemented in hardware, software, firmware, middleware, microcode, or a combination thereof. For a hardware implementation, the processing units may be implemented within one or more application SPECIFIC INTEGRATED circuits (asics), digital signal processors (DIGITAL SIGNAL processing, dsps), digital signal processing devices (DSP DEVICE, DSPD), programmable logic devices (Programmable Logic Device, plds), field-programmable gate arrays (field-programmable GATE ARRAY, FPGA), general purpose processors, controllers, micro-controllers, microprocessors, other electronic units for performing the above-described functions of the application, or a combination thereof.
For a software implementation, the techniques described herein may be implemented with modules (e.g., procedures, functions, and so on) that perform the functions described herein. The software codes may be stored in a memory and executed by a processor. The memory may be implemented within the processor or external to the processor.
Specifically, the processor 603 is further configured to execute the steps of the method for designing a blade according to the foregoing technical solution when executing the computer program, which will not be described herein.
It should be noted that: the technical schemes described in the embodiments of the present disclosure may be arbitrarily combined without any conflict.
The foregoing is merely a specific implementation of the embodiments of the disclosure, but the protection scope of the embodiments of the disclosure is not limited thereto, and any person skilled in the art can easily think about changes or substitutions within the technical scope of the embodiments of the disclosure, and the changes or substitutions are covered by the protection scope of the embodiments of the disclosure. Therefore, the protection scope of the embodiments of the present disclosure shall be subject to the protection scope of the claims.
Claims (10)
1. A method of designing a camber line of a blade, the method comprising:
collecting a plurality of sampling points on an initial camber line corresponding to the uniform section of the blade;
Acquiring second coordinates of a plurality of control points of the Bezier curve corresponding to the initial mean camber line by using a singular value decomposition method based on the acquired first coordinates corresponding to the plurality of sampling points;
The acquiring, based on the first coordinates corresponding to the acquired plurality of sampling points, second coordinates of a plurality of control points of the bezier curve corresponding to the initial mean camber line by using a singular value decomposition method includes:
Obtaining an arc length parameter array corresponding to the plurality of sampling points based on the first coordinates;
Based on the arc length parameter array, obtaining bases corresponding to different control points at each sampling point;
Obtaining a coefficient matrix based on the bases corresponding to different control points at each sampling point;
Calculating to obtain the eigenvalue decomposition matrix by a singular value decomposition method based on the coefficient matrix;
calculating second coordinates corresponding to a plurality of control points of the Bezier curve based on the eigenvalue decomposition matrix;
Adjusting the bending degree of the Bezier curve corresponding to the initial mean camber line by adjusting the second coordinates of the control points; the Bezier curve after the bending is adjusted is used for representing a final camber line corresponding to the constant section of the blade.
2. The method according to claim 1, wherein the obtaining, based on the first coordinates corresponding to the collected plurality of sampling points, the second coordinates of the plurality of control points of the bezier curve corresponding to the initial mean camber line using a singular value decomposition method includes:
Acquiring a first coordinate array corresponding to the plurality of sampling points based on the first coordinates;
Based on a first coordinate array corresponding to the plurality of sampling points, normalizing to obtain an arc length parameter array corresponding to the plurality of sampling points;
Acquiring slopes corresponding to an abscissa and an ordinate in the third coordinate and the fourth coordinate respectively based on a third coordinate corresponding to an initial endpoint in two endpoints of the initial mean camber line and a fourth coordinate corresponding to an end endpoint in the two endpoints;
calculating to obtain a characteristic value decomposition matrix corresponding to the Bezier curve with a set order by utilizing a singular value decomposition method based on the arc length parameter array;
calculating a second coordinate array corresponding to a plurality of control points of the Bezier curve based on the eigenvalue decomposition matrix; wherein the second coordinate array is used for representing the second coordinates.
3. The design method according to claim 2, wherein the normalizing, based on the first coordinate arrays corresponding to the plurality of sampling points, obtains arc length parameter arrays corresponding to the plurality of sampling points, includes:
Calculating to obtain a distance L i between the ith sampling point and the ith sampling point based on a first coordinate corresponding to the ith-1 sampling point in the first coordinate array and a first coordinate corresponding to the ith sampling point in the first coordinate array; wherein i is more than or equal to 1 and less than or equal to N, N represents the number of sampling points;
Based on the i-1 th sampling point and the distance L i between the i-1 th sampling points, calculating according to t i'=ti-1+Li to obtain an initial arc length parameter t i' corresponding to the i-th sampling points; wherein t i-1 represents the normalized arc length parameter corresponding to the i-1 th sampling point;
According to Obtaining the arc length parameter after normalization processing corresponding to the ith sampling point; wherein/>Representing an initial arc length parameter corresponding to an Nth sampling point;
And obtaining arc length parameter arrays corresponding to the sampling points according to the normalized arc length parameters corresponding to all the sampling points.
4. The design method according to claim 2, wherein the calculating, based on the arc length parameter array, the eigenvalue decomposition matrix corresponding to the bezier curve of the set order by using a singular value decomposition method includes:
acquiring the number of control points of the Bezier curve, the column number and the line number of a coefficient matrix corresponding to the first coordinate based on the acquired number of sampling points and the order of the Bezier curve;
The basis corresponding to the different control points at each sampling point is obtained according to the following equation:
wherein, Representing the normalized arc length parameter/>, corresponding to the ith sampling pointThe base corresponding to the j-th control point below is more than or equal to 0, and is more than or equal to n, wherein n represents the order of the Bessel curve;
Obtaining the coefficient matrix based on the bases corresponding to different control points at each sampling point;
And calculating the eigenvalue decomposition matrix by using a singular value decomposition method based on the coefficient matrix.
5. The method according to claim 4, wherein the obtaining the number of control points corresponding to the bezier curve, the number of columns and the number of rows of the coefficient matrix corresponding to the first coordinate based on the number of collected sampling points and the order of the bezier curve includes:
The number of control points of the Bezier curve is equal to the order +1 of the Bezier curve;
The number of columns of the coefficient matrix is equal to 2 times the number of control points of the Bezier curve minus 6;
The number of rows of the coefficient matrix is equal to 2 times the number of the sampling points.
6. The design method according to claim 4, wherein the obtaining the coefficient matrix based on the basis corresponding to the different control point at each sampling point includes:
Based on the basis corresponding to different control points at each sampling point, the obtained coefficient matrix M is:
。
7. the method according to claim 4, wherein calculating a second coordinate array corresponding to a plurality of control points of the bezier curve based on the eigenvalue decomposition matrix comprises:
Acquiring numerical contribution arrays of the two endpoints and control points adjacent to the two endpoints to all sampling points based on the base corresponding to different control points at each sampling point, the slopes corresponding to the third coordinate and the abscissa and the ordinate in the third coordinate, and the slopes corresponding to the fourth coordinate and the abscissa and the ordinate in the fourth coordinate; wherein the numerical contribution array is expressed as And (2) and，；/>Representing an abscissa in the third coordinate; a1 represents a slope corresponding to an abscissa in the third coordinate; /(I)Representing an abscissa in the fourth coordinate; a2 represents a slope corresponding to an abscissa in the fourth coordinate; /(I)Representing an ordinate in the third coordinate; a3 represents a slope corresponding to an ordinate in the third coordinate; /(I)Representing an ordinate in the fourth coordinate; a4 represents a slope corresponding to an ordinate in the fourth coordinate;
Calculating a third coordinate array corresponding to the plurality of sampling points based on the first coordinate array and the numerical contribution array;
calculating to obtain a constant array based on the first coordinate array and the third coordinate array;
and obtaining a second coordinate array corresponding to the control points based on the eigenvalue decomposition matrix and the constant array.
8. A device for designing a camber line of a blade, the device comprising: a collection unit, an acquisition unit, and an adjustment unit; wherein,
The acquisition section is configured to: collecting a plurality of sampling points on an initial camber line corresponding to the uniform section of the blade;
the acquisition section is configured to: acquiring second coordinates of a plurality of control points of the Bezier curve corresponding to the initial mean camber line by using a singular value decomposition method based on the acquired first coordinates corresponding to the plurality of sampling points;
The adjustment section is configured to: adjusting the bending degree of the Bezier curve corresponding to the initial mean camber line by adjusting the second coordinates of the control points; the Bezier curve after the bending is adjusted is used for representing a final camber line corresponding to the constant section of the blade;
Wherein the acquisition section is configured to:
Obtaining an arc length parameter array corresponding to the plurality of sampling points based on the first coordinates;
Based on the arc length parameter array, obtaining bases corresponding to different control points at each sampling point;
Obtaining a coefficient matrix based on the bases corresponding to different control points at each sampling point;
Calculating to obtain the eigenvalue decomposition matrix by a singular value decomposition method based on the coefficient matrix;
and calculating second coordinates corresponding to a plurality of control points of the Bezier curve based on the eigenvalue decomposition matrix.
9. A computing device, the computing device comprising: a communication interface, a memory and a processor; the components are coupled together by a bus system; wherein,
The communication interface is used for receiving and transmitting signals in the process of receiving and transmitting information with other external network elements;
The memory is used for storing a computer program capable of running on the processor;
the processor is configured to perform the steps of the method for designing a camber line of a blade according to any one of claims 1 to 7 when the computer program is run.
10. A computer storage medium storing a program for designing a camber line of a blade, which when executed by at least one processor implements the steps of the method for designing a camber line of a blade according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410123620.8A CN117648779B (en) | 2024-01-30 | 2024-01-30 | Method, device, equipment and computer storage medium for designing camber line of blade |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410123620.8A CN117648779B (en) | 2024-01-30 | 2024-01-30 | Method, device, equipment and computer storage medium for designing camber line of blade |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117648779A CN117648779A (en) | 2024-03-05 |
| CN117648779B true CN117648779B (en) | 2024-04-19 |
Family
ID=90046413
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410123620.8A Active CN117648779B (en) | 2024-01-30 | 2024-01-30 | Method, device, equipment and computer storage medium for designing camber line of blade |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117648779B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20040043652A (en) * | 2002-11-19 | 2004-05-24 | 주식회사 에어로네트 | Method For Designig A Low Noise Impeller Using An Inver Design Method |
| CN112231901A (en) * | 2020-10-10 | 2021-01-15 | 北京航空航天大学 | Offset calculation method of T spline surface |
| CN115544878A (en) * | 2022-09-28 | 2022-12-30 | 珠海格力电器股份有限公司 | Design method of blade profile, blade, wind wheel, fan and oil fume suction device |
| CN117036422A (en) * | 2023-08-25 | 2023-11-10 | 新石器慧通（北京）科技有限公司 | Method, device, equipment and storage medium for tracking lane lines |
| CN117171924A (en) * | 2023-11-03 | 2023-12-05 | 陕西空天信息技术有限公司 | Blade design method, device, equipment and computer storage medium |
| CN117171923A (en) * | 2023-11-02 | 2023-12-05 | 陕西空天信息技术有限公司 | Blade design method, device, equipment and computer storage medium |
-
2024
- 2024-01-30 CN CN202410123620.8A patent/CN117648779B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR20040043652A (en) * | 2002-11-19 | 2004-05-24 | 주식회사 에어로네트 | Method For Designig A Low Noise Impeller Using An Inver Design Method |
| CN112231901A (en) * | 2020-10-10 | 2021-01-15 | 北京航空航天大学 | Offset calculation method of T spline surface |
| CN115544878A (en) * | 2022-09-28 | 2022-12-30 | 珠海格力电器股份有限公司 | Design method of blade profile, blade, wind wheel, fan and oil fume suction device |
| CN117036422A (en) * | 2023-08-25 | 2023-11-10 | 新石器慧通（北京）科技有限公司 | Method, device, equipment and storage medium for tracking lane lines |
| CN117171923A (en) * | 2023-11-02 | 2023-12-05 | 陕西空天信息技术有限公司 | Blade design method, device, equipment and computer storage medium |
| CN117171924A (en) * | 2023-11-03 | 2023-12-05 | 陕西空天信息技术有限公司 | Blade design method, device, equipment and computer storage medium |
Non-Patent Citations (3)
| Title |
|---|
| Renhui Zhang .etal.Knowledge Mining of Low Specific Speed Centrifugal Pump Impeller Based on Proper Orthogonal Decomposition Method.Journal of Thermal Science.2020,第30卷全文. * |
| 尚仁操 等.基于参数法和贝塞尔曲线的涡轮叶片造型及其优化.机械设计与制造.2007,(第08期),全文. * |
| 樊会元 等.扩压器叶片形线逆命题设计的一种新方法.航空动力学报.1999,(第01期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117648779A (en) | 2024-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110175641B (en) | Image recognition method, device, equipment and storage medium | |
| CN117648779B (en) | Method, device, equipment and computer storage medium for designing camber line of blade | |
| CN112149060B (en) | Steam turbine low-pressure final-stage temperature field prediction method and device | |
| CN117171924B (en) | Blade design method, device, equipment and computer storage medium | |
| EP4287061A1 (en) | Soil body state data calculation method and device based on boundary surface plasticity model | |
| CN115438571B (en) | Ground wave radar wave field calculation method and device based on machine learning | |
| CN114060102B (en) | Method and device for determining outlet metal angle of guide vane blade | |
| CN112149901B (en) | Steam turbine low-pressure final-stage pressure field prediction method and device | |
| CN113780523A (en) | Image processing method, image processing device, terminal equipment and storage medium | |
| CN113469506A (en) | User baseline load estimation method, terminal and computer readable storage medium | |
| CN116484771B (en) | Method and device for generating CFD grid of axial flow compressor | |
| CN116244865B (en) | Method and device for finite element modeling of axial flow impeller and computer storage medium | |
| CN115952834A (en) | Batch layer normalization method and device applied to Transformer | |
| CN116049658B (en) | Wind turbine generator abnormal data identification method, system, equipment and medium | |
| CN113051688B (en) | Blade profile shape geometry calculation method, device and equipment | |
| CN116259328A (en) | Post-training quantization method, apparatus and storage medium for audio noise reduction | |
| CN114065913A (en) | Model quantization method and device and terminal equipment | |
| CN112784840A (en) | License plate recognition method, device, equipment and medium | |
| CN114881913A (en) | Image defect detection method and device, electronic equipment and storage medium | |
| CN116520289B (en) | Intelligent control method and related device for hardware machining process | |
| CN116205605B (en) | Intelligent evaluation method, system and medium for quality of science and technology project file | |
| CN118094076A (en) | Apparatus, method, and associated product for floating point quantization operations in neural network models | |
| CN115984705B (en) | Large-scale remote sensing network image bridge segmentation method based on convolution network | |
| CN114419064B (en) | Mammary gland area image segmentation method based on RN-DoubleU-Net network | |
| CN111753391B (en) | Power transmission line design method and system based on maximum field intensity of surface of split conductor |
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 |