Method and device for identifying indirect contact boundaries between parts of thermal equipment

Abstract

This invention provides a method and apparatus for identifying contact boundaries between parts of thermal equipment. The method includes: identifying multiple protrusions of the first part and the second part from multiple points on the contour boundaries of the first part and the second part; identifying multiple first-type breakpoints on the contour boundary of the first part based on the first type of shortest distance from each of the multiple protrusions of the second part to the contour boundary of the first part, and identifying multiple first-type breakpoints on the contour boundary of the second part based on the first type of shortest distance from each of the multiple protrusions of the first part to the contour boundary of the second part; identifying multiple second-type breakpoints on the contour boundary of the first part based on the second type of shortest distance from each of the multiple first-type breakpoints of the second part to the contour boundary of the first part, and identifying multiple second-type breakpoints on the contour boundary of the second part based on the second type of shortest distance from each of the multiple first-type breakpoints of the first part to the contour boundary of the second part; breaking the contour boundaries of the first part and the second part at the positions of the first-type breakpoints and the second-type breakpoints to obtain multiple broken boundaries included in the first part and the second part respectively; and determining multiple target boundaries that coincide with each other in position among the multiple broken boundaries of the first part and the second part as the contact boundary between the first part and the second part.

Description

Technical Field

[0001] This invention relates to the field of thermal analysis and thermal management technology for gas turbines, and specifically to a method and apparatus for identifying contact boundaries between mechanical parts. Background Technology

[0002] Temperature field analysis of the entire gas turbine is a key technology in the safety design of modern aero-engines and gas turbines. It achieves accurate simulation of the entire heat transfer process of the entire machine through multi-physics field coupling simulation.

[0003] Identifying the contact positions of multiple components in a gas turbine is a prerequisite for overall temperature field analysis. The definition of contact pairs directly impacts the accuracy of simulating thermal resistance and heat conduction paths within the overall temperature field. Currently, contact pairs are defined manually. Engineers pre-identify key contact areas (such as rotor disk contact surfaces, blade tenons and slots, bolt connection surfaces, and stop mating surfaces) based on assembly relationships. Then, geometric preprocessing is performed, rounding off sharp corners and edges that could lead to convergence failure to ensure the stability of the contact algorithm. This process consumes the majority of the time in the entire cycle of overall temperature field analysis, is not only time-consuming but also inherently prone to human error. Therefore, an automated and accurate rapid contact identification method is urgently needed to achieve efficient, fast, and accurate simulation of the overall temperature field. Summary of the Invention

[0004] In view of the above problems, the present invention provides a method for identifying the contact boundary between parts of thermal equipment.

[0005] One aspect of the present invention provides a method for identifying contact boundaries between parts of a thermal device, comprising identifying contact boundaries between a plurality of parts included in the thermal device, wherein, for any first part and second part that are positionally adjacent to each other among the plurality of parts, the contact boundary between the first part and the second part is identified by the following method:

[0006] Identify multiple protrusions of each of the first and second parts from multiple points on the contour boundaries of the first and second parts;

[0007] Based on the first type of shortest distance from each of the multiple protrusions of the second part to the contour boundary of the first part, multiple first type of breakpoints on the contour boundary of the first part are identified; and based on the first type of shortest distance from each of the multiple protrusions of the first part to the contour boundary of the second part, multiple first type of breakpoints on the contour boundary of the second part are identified.

[0008] Based on the second type of shortest distance from each of the first type of breakpoints of the second part to the contour boundary of the first part, multiple second type breakpoints on the contour boundary of the first part are identified; and based on the second type of shortest distance from each of the first type of breakpoints of the first part to the contour boundary of the second part, multiple second type breakpoints on the contour boundary of the second part are identified.

[0009] At the locations of the first and second type breakpoints, the contour boundaries of the first and second parts are broken to obtain multiple broken boundaries included in each of the first and second parts.

[0010] Among the multiple breaking boundaries of the first part and the second part, the multiple target boundaries that coincide with each other in position are determined as the contact boundaries between the first part and the second part.

[0011] Another aspect of the present invention provides an apparatus for identifying contact boundaries between parts of a thermal device. The apparatus is used to identify contact boundaries between a plurality of parts included in the thermal device. Specifically, for any first part and a second part that are positionally adjacent to each other among the plurality of parts, the apparatus is used to identify the contact boundary between the first part and the second part. The apparatus includes:

[0012] The first identification module is used to identify multiple protrusions of the first part and the second part from multiple points on the contour boundary of the first part and the second part respectively.

[0013] The second identification module is used to identify multiple first-type breakpoints on the contour boundary of the first part based on the first type of shortest distance from each of the multiple protrusions of the second part to the contour boundary of the first part, and to identify multiple first-type breakpoints on the contour boundary of the second part based on the first type of shortest distance from each of the multiple protrusions of the first part to the contour boundary of the second part.

[0014] The third identification module is used to identify multiple second-type breakpoints on the contour boundary of the first part based on the second-type shortest distance from each of the multiple first-type breakpoints of the second part to the contour boundary of the first part, and to identify multiple second-type breakpoints on the contour boundary of the second part based on the second-type shortest distance from each of the multiple first-type breakpoints of the first part to the contour boundary of the second part.

[0015] The break module is used to break the contour boundaries of the first part and the second part at the locations of the first type of breakpoint and the second type of breakpoint, so as to obtain multiple break boundaries included in the first part and the second part respectively.

[0016] The determination module is used to determine the multiple target boundaries that coincide with each other in position among the multiple breaking boundaries of the first part and the second part as the contact boundary between the first part and the second part. Attached Figure Description

[0017] The above-described features, other objects, and advantages of the present invention will become clearer from the following description of embodiments of the invention with reference to the accompanying drawings, in which:

[0018] Figure 1 A flowchart illustrating a method for identifying contact boundaries between parts of a thermal equipment according to an embodiment of the present invention is shown.

[0019] Figure 2 A schematic diagram illustrating the principle of identifying protrusions on a component according to an embodiment of the present invention is shown.

[0020] Figure 3 This schematic diagram illustrates the principle of determining the landing point on a part according to an embodiment of the present invention;

[0021] Figure 4 This schematic diagram illustrates the principle of determining a first type of breakpoint and a second type of breakpoint on a part according to an embodiment of the present invention.

[0022] Figure 5 A schematic diagram illustrating the principle of identifying contact boundaries between parts according to an embodiment of the present invention is shown.

[0023] Figure 6 The illustration schematically shows the identification results of the contact boundary between engine parts according to an embodiment of the present invention;

[0024] Figure 7 The diagram schematically illustrates a structural block diagram of a device for identifying contact boundaries between parts of a thermal equipment according to an embodiment of the present invention. Detailed Implementation

[0025] Hereinafter, embodiments of the present invention will be described with reference to the accompanying drawings. However, it should be understood that these descriptions are exemplary only and are not intended to limit the scope of the invention. In the following detailed description, numerous specific details are set forth to provide a thorough understanding of the embodiments of the invention for ease of explanation. However, it will be apparent that one or more embodiments may be practiced without these specific details. Furthermore, descriptions of well-known structures and techniques are omitted in the following description to avoid unnecessarily obscuring the concept of the invention.

[0026] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the invention. The terms “comprising,” “including,” etc., as used herein indicate the presence of the stated features, steps, operations, and / or components, but do not exclude the presence or addition of one or more other features, steps, operations, or components.

[0027] All terms used herein (including technical and scientific terms) have the meanings commonly understood by those skilled in the art, unless otherwise defined. It should be noted that the terms used herein are to be interpreted in a manner consistent with the context of this specification, and not in an idealized or overly rigid way.

[0028] When using expressions such as "at least one of A, B and C", they should generally be interpreted in accordance with the meaning that is commonly understood by those skilled in the art (e.g., "a system having at least one of A, B and C" should include, but is not limited to, a system having A alone, a system having B alone, a system having C alone, a system having A and B, a system having A and C, a system having B and C, and / or a system having A, B and C, etc.).

[0029] One aspect of the present invention provides a method for identifying contact boundaries between parts of a thermal equipment.

[0030] Figure 1 A flowchart illustrating a method for identifying contact boundaries between parts of a thermal equipment according to an embodiment of the present invention is shown.

[0031] This method includes identifying the contact boundaries between multiple parts comprising a thermal device, such as Figure 1 As shown, for any first and second parts that are geographically adjacent to each other among a plurality of parts, the contact boundary between the first and second parts is identified using the following method:

[0032] Step S101: Identify multiple protrusions of each of the first and second parts from multiple points on the contour boundaries of the first and second parts;

[0033] Step S102: Based on the first type of shortest distance from each of the multiple protrusions of the second part to the contour boundary of the first part, identify multiple first type of breakpoints on the contour boundary of the first part; and based on the first type of shortest distance from each of the multiple protrusions of the first part to the contour boundary of the second part, identify multiple first type of breakpoints on the contour boundary of the second part.

[0034] Step S103: Based on the second type of shortest distance from each of the multiple first type breakpoints of the second part to the contour boundary of the first part, identify multiple second type breakpoints on the contour boundary of the first part; and based on the second type of shortest distance from each of the multiple first type breakpoints of the first part to the contour boundary of the second part, identify multiple second type breakpoints on the contour boundary of the second part.

[0035] Step S104: At the locations of the first type of breakpoint and the second type of breakpoint, break the contour boundaries of the first part and the second part to obtain multiple break boundaries included in each of the first part and the second part.

[0036] Step S105: Among the multiple break boundaries of the first part and the second part, the multiple target boundaries that coincide with each other in position are determined as the contact boundary between the first part and the second part.

[0037] According to embodiments of the present invention, the thermal equipment may be a gas turbine or an aircraft engine, etc.

[0038] Temperature field analysis of the entire gas turbine is a key technology in the safety design of modern aero-engines and gas turbines. It achieves accurate simulation of the entire heat transfer process through multiphysics coupling simulation. This technology can precisely track the heat conduction path from high to low temperature, from internal to external, and from start-up to shutdown, which is of great significance for engine hot clearance control, cooling system efficiency evaluation, and strength analysis verification. In engineering practice, it is necessary to comprehensively consider the nonlinear thermophysical parameters of materials, transient boundary conditions, and contact thermal resistance models to achieve accurate temperature field simulation.

[0039] Identifying the contact locations of multiple components in a gas turbine is a prerequisite for overall temperature field analysis. This is because the contact locations determine the subsequent heat transfer calculation methods; heat transfer calculations are performed at these locations using methods based on heat conduction and heat transfer. If the contact locations are not accurately identified, errors will occur in the heat transfer calculation process. For example, if a contact location that should be used for heat conduction calculations is not identified, the temperature calculation result at that location will deviate significantly from the correct temperature distribution. Therefore, the definition of contact locations directly affects the simulation accuracy of thermal resistance and heat conduction paths in the overall temperature field of the gas turbine.

[0040] In the above method of the embodiments of the present invention, taking any first part and second part that are adjacent to each other in position among a plurality of parts as an example, the method for identifying the contact boundary between the two parts is described. Other parts are identified using the same method, and will not be repeated here.

[0041] In the above method, the protrusions of the first and second parts are first identified, and the breakpoints of the first and second parts are identified based on the distance from the protrusion (breakpoint) to the boundary of the adjacent part. At the breakpoint location, the boundary of the part is broken. After all parts of the whole machine have completed the breaking process, if the boundary of one part coincides with the boundary of another part, the location of the coincident boundary is the contact position of the temperature field of the whole machine.

[0042] The above method can be executed by a processor in an electronic device, wherein the input data are the position coordinates of multiple points on the contour boundaries of the first and second parts.

[0043] Here, a convex point refers to a point on the contour boundary of the first and second parts where the turning angle in the tangential direction is greater than or equal to a predetermined angle threshold (e.g., a turning angle greater than or equal to 15°). For example, the four corner points of a square part are convex points.

[0044] After identifying the protrusions, the breakpoints of the first and second parts are then identified. The breakpoints correspond to the endpoints of the contact boundaries between the first and second parts; therefore, only by accurately identifying the breakpoints can the contact positions be determined.

[0045] Furthermore, the breakpoints may include a first type of breakpoint and a second type of breakpoint. The first type of breakpoint is determined based on the distance from the protrusion to the boundary of the adjacent part, and the second type of breakpoint is determined based on the distance from the first type of breakpoint to the boundary of the adjacent part.

[0046] After identifying the breakpoints of the first and second parts, the contour boundaries of the first and second parts are broken at the breakpoint locations, resulting in multiple broken boundaries for each part. The breaking process refers to splitting one boundary of a part into two or more boundaries at the breakpoint location.

[0047] Finally, by executing step S105, multiple target boundaries that coincide with each other in position among the multiple break boundaries of the first and second parts are determined as the contact boundaries between the first and second parts. Here, a coincident boundary refers to a boundary where the shortest distance between all points on the boundary of the first part and the boundary of the second part is less than or equal to a predetermined distance threshold. The method for determining boundary coincidence is to collect sample points on both boundaries according to their length proportions. When the distance between the sample points corresponding to each length proportion of the two boundaries is less than the predetermined distance threshold, the two boundaries are considered to coincide.

[0048] It should be noted that the protrusions and breaks in the embodiments of the present invention do not refer to a single point, but rather to a set of points on the boundary of the part. That is, a part boundary is allowed to have multiple protrusions and breaks.

[0049] According to embodiments of the present invention, the above method can automatically identify the contact boundaries of parts, and can quickly identify parts regardless of their structural complexity. It possesses the capability to identify contact positions of a large number of parts, solving the technical problem of high workload in manual identification; furthermore, it offers high completeness, high efficiency, and low error, effectively avoiding the tediousness and repetitiveness of manual identification.

[0050] The ability to identify large-scale part contact positions is demonstrated by the fact that the contact position identification process uses significant feature points such as convex points and breakpoints for distance calculation. Compared to traditional methods that require pairwise calculations between all points on the boundary, this significantly reduces the computational load and saves computational costs. Furthermore, when a point breaks the boundary, it is determined by the shortest distance from the point to the part boundary, which has the characteristic of small error. Moreover, the shortest distance calculation has already been completed before the break, eliminating the need for repeated calculations. This achieves both low identification error and saves computational resources.

[0051] Furthermore, since the contact points of the gas turbine components only occur at the boundaries of the protrusions and breaks, contact identification based on these protrusions and breaks ensures the integrity of the identification and features high identification completeness. By calculating the shortest distance from the protrusions and breaks to the adjacent boundaries, the contact identification error is minimized. By requiring the shortest distance to be lower than a predetermined distance threshold, the range of contact identification is strictly guaranteed to be less than the defined geometric error, preventing the expansion of the contact position range and thus exhibiting low identification error.

[0052] Furthermore, by statistically analyzing all broken boundaries, it is ensured that the identification is not duplicated, and the contact boundaries of the two parts can be matched one-to-one. After the parts are broken, the contact position is determined by statistically analyzing all broken boundaries, so there will be no cases where a boundary should be broken but is not, thus the completeness of the identification is high.

[0053] As can be seen, the method of this invention solves the technical problems involved in the temperature field analysis of the whole gas turbine, such as cumbersome manual identification of contact positions, incomplete identification, low identification efficiency, large identification error, and repeated identification. At the same time, it improves the stability of the contact algorithm and provides a solution for accelerating the construction of intelligent and rapid design level of high-end gas turbine equipment.

[0054] According to an embodiment of the present invention, a convex point represents a point among a plurality of points on the contour boundary of the first part and the second part, where the turning angle in the tangential direction is greater than or equal to a predetermined angle threshold (e.g., the turning angle is greater than or equal to 15°).

[0055] According to an embodiment of the present invention, identifying a plurality of protrusions of the first part and the second part from a plurality of points on the contour boundaries of the first part and the second part includes: determining a plurality of reference points from the plurality of points on the contour boundaries of the first part and the second part, wherein the plurality of reference points include at least the inflection points of the contour boundaries; and then identifying a plurality of protrusions from the plurality of reference points.

[0056] The process of identifying multiple convex points from multiple reference points includes performing convex point identification processing on each reference point. The convex point identification processing includes:

[0057] First, determine the first and second adjacent reference points that are adjacent to the reference point;

[0058] Next, the angle between the first vector and the second vector is calculated, where the first vector represents the vector formed from the first adjacent reference point to the reference point, and the second vector represents the vector formed from the reference point to the second adjacent reference point.

[0059] Then, based on the angle between the first vector and the second vector, it is determined whether the reference point is a convex point. For example, if the angle between the first vector and the second vector is greater than or equal to a predetermined angle threshold (e.g., the predetermined angle threshold is 15 degrees), the reference point is determined to be a convex point; if the angle between the first vector and the second vector is less than the predetermined angle threshold, the reference point is determined not to be a convex point.

[0060] Figure 2 A schematic diagram illustrating the principle of identifying protrusions on a component according to an embodiment of the present invention is shown below. Figure 2 The method for identifying convex points will be explained.

[0061] According to an embodiment of the present invention, the multiple points on the contour boundaries of the first part and the second part can be all the points on the contour boundaries. In order to save computation, some representative points can be selected from all the points as reference points for calculation. The reference points need to include at least the inflection points of the contour boundaries to avoid missing the protrusions.

[0062] like Figure 2 Taking the first part as an example, this paper explains the method for identifying protrusions on the first part.

[0063] like Figure 2 Among all the points on the outline boundary of the first part, select some representative reference points, such as points 1, 2, 3, 4, 5, 6, and 7 in the figure.

[0064] For any given point, calculate the angle between the first vector and the second vector.

[0065] For example, for point 2, the first adjacent reference point and the second adjacent reference point are point 3 and point 1 respectively. Calculate the angle between the first vector and the second vector, that is, calculate the angle between the vector from point 3 to point 2 and the vector from point 2 to point 1. Points 1, 2 and 3 are located on the same straight line, so the angle between the vectors is zero degrees, which is less than 15 degrees. Therefore, point 2 is not a convex point.

[0066] For example, for point 1, the first adjacent reference point and the second adjacent reference point are point 2 and point 6 respectively. Calculate the angle between the first vector and the second vector, that is, calculate the angle between the vector from point 2 to point 1 and the vector from point 1 to point 6, as shown in angle α in the figure. Assuming that angle α is 30 degrees, which is greater than 15 degrees, then point 1 is a convex point.

[0067] According to an embodiment of the present invention, identifying multiple protrusions of the first part and the second part from multiple points on the contour boundaries of the first part and the second part may also include: directly identifying protrusions from multiple points on the contour boundaries of the first part and the second part.

[0068] The convexity identification process includes: First, for each point, determining the first and second adjacent points. Then, calculating the angle between the vector formed from the first adjacent point to the point and the vector formed from the point to the second adjacent point. If the angle is greater than or equal to a predetermined angle threshold (e.g., 15 degrees), the point is determined to be a convexity; if the angle is less than the predetermined angle threshold, the point is determined not to be a convexity.

[0069] For specific methods on convexity identification for each point, please refer to the aforementioned relevant information. Figure 2 The description will not be repeated here.

[0070] According to an embodiment of the present invention, in step S102 above, the method for identifying multiple first-type breakpoints on the contour boundary of the first part based on the first type of shortest distance from each of the multiple protrusions of the second part to the contour boundary of the first part includes the following method:

[0071] First, the point corresponding to the first type of shortest distance among the multiple points included in the contour boundary of the first part is determined as the first type of landing point; then, based on the first type of shortest distance and a predetermined distance threshold, multiple first type of breakpoints are determined from the multiple first type of landing points on the contour boundary of the first part.

[0072] Furthermore, the method for identifying multiple first-type breakpoints on the contour boundary of the second part based on the first type of shortest distance from each of the multiple protrusions of the first part to the contour boundary of the second part is the same as the aforementioned method for identifying multiple first-type breakpoints on the contour boundary of the first part, and will not be repeated here.

[0073] According to an embodiment of the present invention, in step S103 above, identifying multiple second-type breakpoints on the contour boundary of the first part based on the second-type shortest distance from each of the multiple first-type breakpoints of the second part to the contour boundary of the first part includes:

[0074] First, the point corresponding to the second type of shortest distance among the multiple points included in the contour boundary of the first part is determined as the second type of landing point; then, based on the second type of shortest distance and a predetermined distance threshold, multiple second type of breakpoints are determined from the multiple second type of landing points on the contour boundary of the first part.

[0075] Furthermore, the method for identifying multiple second-type breakpoints on the contour boundary of the second part based on the second-type shortest distance from each of the multiple first-type breakpoints of the first part to the contour boundary of the second part is the same as the aforementioned method for identifying multiple second-type breakpoints on the contour boundary of the first part, and will not be repeated here.

[0076] Figure 3 This schematic diagram illustrates the principle of determining the landing point on a part according to an embodiment of the present invention; Figure 4 A schematic diagram illustrating the principle of determining a first type of breakpoint and a second type of breakpoint on a part according to an embodiment of the present invention is shown below. Figure 3 , Figure 4 The methods for identifying the landing point and breakpoint described above are explained.

[0077] According to an embodiment of the present invention, the shortest distance from the protrusion of a part to the contour boundary of an adjacent part is defined as the first type of shortest distance, and the shortest distance from the breakpoint of a part to the contour boundary of an adjacent part is defined as the second type of shortest distance.

[0078] The outline boundary of adjacent parts refers to the complete loop formed by the outlines of adjacent parts (the outline of each part is a closed loop). The distance from the convex point of a part to the outline boundary of an adjacent part can be multiple sets of distance values ​​from the convex point of the part to multiple points that make up the outline boundary of the adjacent parts. The shortest distance is selected from these values ​​as the first type of shortest distance from the convex point of the part to the outline boundary of the adjacent part. Furthermore, among the multiple points on the outline boundary of the adjacent parts, the point corresponding to the first type of shortest distance is the first type of landing point.

[0079] The distance from the breakpoint of a part to the contour boundary of an adjacent part can be a set of distance values ​​from the breakpoint of the part to multiple points that make up the contour boundary of the adjacent part. The shortest distance is selected from these values ​​as the second type of shortest distance from the breakpoint of the part to the contour boundary of the adjacent part. Furthermore, among the multiple points on the contour boundary of the adjacent part, the point corresponding to the second type of shortest distance is the second type of landing point.

[0080] like Figure 3 The method for determining the first type of landing point of the contour boundary of the first part based on the protrusions on the second part, given two adjacent first and second parts, is as follows:

[0081] For example, the distance from protrusion 1 on the second part to the contour boundary of the first part can be multiple sets of distance values ​​from protrusion 1 to multiple points included in the contour boundary of the first part. Among these distances, the distance from point 1 to point 5 on the contour boundary of the first part is the shortest. Therefore, the distance from point 1 to point 5 is the first type of shortest distance from protrusion 1 on the second part to the contour boundary of the first part. Point 5, which corresponds to the first type of shortest distance among the multiple points included in the contour boundary of the first part, is determined as the first type of landing point. This process is repeated to determine the first type of landing points from other protrusions on the second part to the contour boundary of the first part, ultimately resulting in multiple first type of landing points on the contour boundary of the first part.

[0082] Subsequently, based on the first type of shortest distance and a predetermined distance threshold, multiple first type breakpoints are determined from multiple first type landing points on the contour boundary of the first part.

[0083] Specifically, if the shortest distance of the first type is less than a predetermined distance threshold (e.g., less than 1 micrometer, less than 2 micrometers, less than 10 micrometers, etc., where the predetermined distance threshold can be a value customized based on design error), the first type of landing point is converted into a first type of breakpoint. If the shortest distance of the first type is greater than or equal to the predetermined distance threshold, the first type of landing point is not converted into a first type of breakpoint.

[0084] like Figure 3 If the distance between point 1 and point 5 is greater than a predetermined distance threshold (e.g., greater than 1 micrometer), then point 5 is not a first-type breakpoint. Conversely, if the distance between point 1 and point 5 is less than the predetermined distance threshold (e.g., less than 1 micrometer), then point 5 is converted into a first-type breakpoint.

[0085] After identifying multiple first-type breakpoints on the contour boundary of the first part, multiple first-type breakpoints on the contour boundary of the second part can be identified in the same way based on the first-type shortest distance from each of the multiple protrusions of the first part to the contour boundary of the second part. This will not be elaborated further here.

[0086] After identifying the first type of breakpoints on the first and second parts, identify the second type of breakpoints on the first and second parts.

[0087] Specifically, for the first part, based on the second type of shortest distance from the first type of breakpoint on the second part to the contour boundary of the first part, the second type of breakpoint on the contour boundary of the first part is found; for the second part, based on the second type of shortest distance from the first type of breakpoint on the first part to the contour boundary of the second part, the second type of breakpoint on the contour boundary of the second part is found. The method for identifying the second type of breakpoint is similar to the aforementioned method for identifying the first type of breakpoint based on the shortest distance from the convex point of a part to the contour boundary of an adjacent part; the only difference is that when determining the second type of landing point on an adjacent part, it needs to be changed to determine the second type of landing point of the adjacent part based on the shortest distance from the first type of breakpoint of the part to the contour boundary of the adjacent part. Further details will not be elaborated here.

[0088] The following, combined with Figure 4 This example illustrates a method for determining first-type and second-type breakpoints on a part.

[0089] like Figure 4 For the first part, based on the first type of shortest distance from the protrusions on the second part to the contour boundary of the first part, the first type of landing points on the contour boundary of the first part are found, and the first type of breakpoints are determined. For the second part, based on the first type of landing points from the protrusions on the first part to the contour boundary of the second part, and after determining the first type of shortest distance, the first type of breakpoints on the contour boundary of the second part are found.

[0090] For example, based on the first type of shortest distance from the protrusion 4 on the second part to the contour boundary of the first part, the first type of landing point on the contour boundary of the first part is found, and the first type of breakpoint 3 is determined; based on the first type of shortest distance from the protrusion 1 on the first part to the contour boundary of the second part, the first type of landing point on the contour boundary of the second part is found, and the first type of breakpoint 2 is determined (other landing points on the two parts are not breakpoints, and the discrimination method is the same as the aforementioned embodiment, and will not be repeated).

[0091] Then, for the first part, based on the second type of shortest distance from the first type of breakpoint on the second part to the outline boundary of the first part, the second type of landing point on the outline boundary of the first part is found and the second type of breakpoint is determined; for the second part, based on the second type of shortest distance from the first type of breakpoint on the first part to the outline boundary of the second part, the second type of landing point on the outline boundary of the second part is found and the second type of breakpoint is determined.

[0092] For example, based on the second type of shortest distance from the first type of breakpoint 2 on the second part to the outline boundary of the first part, the second type of landing point on the outline boundary of the first part is found, and the second type of breakpoint 1 is determined; for the second part, based on the second type of landing point from the first type of breakpoint 3 on the first part to the outline boundary of the second part, and the second type of shortest distance is determined, the second type of breakpoint 4 on the outline boundary of the second part is found (other landing points on the two parts are not breakpoints, and the discrimination method is the same as the aforementioned embodiment, and will not be repeated).

[0093] Figure 5 A schematic diagram illustrating the principle of identifying contact boundaries between parts according to an embodiment of the present invention is shown.

[0094] Referring to the description of the foregoing embodiments, a method for identifying contact boundaries is exemplarily described. It includes the following steps:

[0095] (a) Calculate the protrusion positions of all parts.

[0096] (b) Calculate the shortest distance from all protrusions of the surrounding parts to the boundary of the parts, and determine the first type of landing point on the boundary of the parts based on the shortest distance.

[0097] (c) When the shortest distance is less than the predetermined distance threshold, the first type of landing point is the first type of breakpoint on the boundary of the part.

[0098] (d) Based on the locations of all breakpoints of the surrounding parts, calculate the new shortest distance from all breakpoints of the surrounding parts to the boundary of the part, and determine the second type of landing point on the boundary of the part based on the new shortest distance.

[0099] (e) When the new shortest distance is less than the predetermined distance threshold, the second type of landing point is the second type of breakpoint on the boundary of the part.

[0100] For example, such as Figure 5 Using the method described above, we determined that there are two Type II breakpoints on the first part: point 1 and point 2, and that there are two Type I breakpoints on the second part: point 3 and point 4.

[0101] (f) At the locations of the first and second type breakpoints, the contour boundaries of the first and second parts are broken to obtain the broken boundaries. For example... Figure 5 We obtain the broken boundary a and the broken boundary b.

[0102] It should be noted that boundary breaking can be performed either by identifying the first type of breakpoint and breaking at that breakpoint first, then identifying the second type of breakpoint and breaking at that breakpoint, or by identifying both types of breakpoints and breaking them together.

[0103] Following steps (a)-(f), identify all breakpoints and break the boundaries of all parts at the breakpoint locations.

[0104] Finally, count all the broken boundaries. When a broken boundary of one part coincides with a broken boundary of another part, then the broken boundary is the contact position of the two parts. Finally, output all the contact positions.

[0105] like Figure 5 The boundary breaks a and b coincide, which is a contact position.

[0106] Figure 6 The illustration schematically shows the identification results of the contact boundary between engine parts according to an embodiment of the present invention.

[0107] like Figure 6 As shown, in the process of analyzing the temperature field of a real gas turbine, given the two-dimensional geometry of the parts, contact identification was performed on all parts of the turbine according to the steps described in this invention, and the contact positions of all parts were finally obtained. The identification results are as follows. Figure 6 As shown.

[0108] According to an embodiment of the present invention, after identifying the contact boundaries between the multiple parts included in the thermal device, the method further includes: performing heat transfer calculations on the thermal device based on a predetermined algorithm, wherein the heat transfer calculations are performed at the contact boundaries of the parts of the thermal device using a predetermined heat conduction algorithm.

[0109] The method described in this invention can be applied to the analysis of the overall temperature field of a gas turbine. After identifying the contact positions between engine parts, different heat transfer algorithms can be applied to calculate the temperature field based on the contact conditions between the parts. For example, a heat conduction algorithm can be used for heat transfer calculation at contact positions, while a convection heat transfer algorithm can be used at non-contact positions. It is evident that accurate identification of the contact positions of multiple gas turbine parts is a prerequisite for overall temperature field analysis. The definition of the contact positions directly affects the simulation accuracy of thermal resistance and heat conduction paths in the overall temperature field.

[0110] Furthermore, under complex operating conditions such as high temperature and start-stop cycles in gas turbines, the nonlinear thermal expansion caused by the overall temperature field will significantly change the distribution of contact positions and contact stresses. It is necessary to update the contact in real time to ensure that the model accurately reflects the contact state under temperature change conditions.

[0111] It's important to clarify that contact in the overall temperature field and contact during overall thermal deformation are entirely different concepts. Contact in the overall temperature field aims to achieve heat transfer; solid-to-solid heat transfer occurs only at the contact points, with no heat transfer at non-contact locations. Contact during overall thermal deformation, on the other hand, ensures that the solid surfaces of different solids do not penetrate each other in the contact area. Therefore, when identifying contact in the overall temperature field, it's necessary to ensure that the contact surfaces of the two parts at the contact points are completely consistent. This requirement means that only some boundaries of the parts may participate in contact heat transfer, exacerbating the nonlinear transmission of the temperature field.

[0112] According to an embodiment of the present invention, since the engine may undergo thermal deformation during heat transfer calculations, the above method further includes:

[0113] During the operation of the thermal equipment, it is determined whether the parts of the thermal equipment have undergone thermal deformation; if the parts of the thermal equipment have undergone thermal deformation, the contour boundary of the parts of the thermal equipment is updated, that is, the left-hand data of the points of the contour boundary is re-determined; then, based on the updated contour boundary, the contact boundary between the multiple parts included in the thermal equipment is re-identified according to the method for identifying contact boundaries in the aforementioned embodiment.

[0114] According to embodiments of the present invention, by updating the contact position in real time to adapt to the influence of engine thermal deformation, the accuracy of heat transfer analysis is ensured. At the same time, the above-mentioned intelligent identification method of contact position can replace manual definition, reduce the cost of manual trial and error, and realize efficient, fast and accurate simulation of the temperature field of the whole machine through automated contact identification. It is especially suitable for scenarios where real-time thermal deformation of the engine requires a lot of manual calculation work.

[0115] Another aspect of the present invention provides a device for identifying contact boundaries between parts of a thermal equipment. Figure 7 The diagram schematically illustrates a structural block diagram of a device 700 for identifying contact boundaries between parts of a thermal equipment according to an embodiment of the present invention. The device 700 is used to identify contact boundaries between a plurality of parts included in a thermal equipment. Specifically, for any first part and a second part that are positionally adjacent to each other among the plurality of parts, the device is used to identify the contact boundary between the first part and the second part. The device 700 includes a first identification module 701, a second identification module 702, a third identification module 703, an interruption module 704, and a determination module 705.

[0116] The first identification module 701 is used to identify multiple protrusions of the first part and the second part from multiple points on the contour boundary of the first part and the second part respectively.

[0117] The second identification module 702 is used to identify multiple first-type breakpoints on the contour boundary of the first part based on the first type of shortest distance from each of the multiple protrusions of the second part to the contour boundary of the first part, and to identify multiple first-type breakpoints on the contour boundary of the second part based on the first type of shortest distance from each of the multiple protrusions of the first part to the contour boundary of the second part.

[0118] The third identification module 703 is used to identify multiple second-type breakpoints on the contour boundary of the first part based on the second-type shortest distance from each of the multiple first-type breakpoints of the second part to the contour boundary of the first part, and to identify multiple second-type breakpoints on the contour boundary of the second part based on the second-type shortest distance from each of the multiple first-type breakpoints of the first part to the contour boundary of the second part.

[0119] The interruption module 704 is used to interrupt the contour boundaries of the first part and the second part at the positions of the first type of interruption point and the second type of interruption point, so as to obtain multiple interruption boundaries included by the first part and the second part respectively.

[0120] The determination module 705 is used to determine the multiple target boundaries that coincide with each other in position among the multiple breaking boundaries of the first part and the second part as the contact boundary between the first part and the second part.

[0121] According to embodiments of the present invention, any plurality of modules among the first identification module 701, the second identification module 702, the third identification module 703, the interruption module 704, and the determination module 705 can be combined into one module, or any one of these modules can be split into multiple modules. Alternatively, at least part of the functionality of one or more of these modules can be combined with at least part of the functionality of other modules and implemented in one module. According to embodiments of the present disclosure, at least one of the first identification module 701, the second identification module 702, the third identification module 703, the interruption module 704, and the determination module 705 can be at least partially implemented as hardware circuitry, such as a field-programmable gate array (FPGA), a programmable logic array (PLA), a system-on-a-chip, a system-on-a-substrate, a system-on-package, an application-specific integrated circuit (ASIC), or implemented in hardware or firmware by any other reasonable means of integrating or packaging the circuitry, or implemented in any one of the three implementation methods of software, hardware, and firmware, or in a suitable combination of any of these. Alternatively, at least one of the first identification module 701, the second identification module 702, the third identification module 703, the interruption module 704, and the determination module 705 can be implemented at least partially as a computer program module, which can perform corresponding functions when the computer program module is run.

[0122] Those skilled in the art will understand that the features described in the various embodiments of the present invention can be combined and / or combined in various ways, even if such combinations or combinations are not explicitly described in the present invention. In particular, the features described in the various embodiments of the present invention can be combined and / or combined in various ways without departing from the spirit and teachings of the present invention. All such combinations and / or combinations fall within the scope of the present invention.

[0123] The embodiments of the present invention have been described above. However, these embodiments are merely illustrative and not intended to limit the scope of the invention. Although various embodiments have been described above, this does not mean that the measures in the various embodiments cannot be used advantageously in combination. Various substitutions and modifications can be made by those skilled in the art without departing from the scope of the invention, and all such substitutions and modifications should fall within the scope of the invention.

Claims

1. A method of identifying contact boundaries between parts of a thermal plant, comprising identifying contact boundaries between a plurality of parts comprising the thermal plant, wherein, For any first and second parts that are geographically adjacent to each other among a plurality of parts, the contact boundary between the first and second parts is identified using the following method: Identify multiple protrusions of each of the first and second parts from multiple points on the contour boundaries of the first and second parts; Based on the first type of shortest distance from each of the multiple protrusions of the second part to the contour boundary of the first part, multiple first type of breakpoints on the contour boundary of the first part are identified; and based on the first type of shortest distance from each of the multiple protrusions of the first part to the contour boundary of the second part, multiple first type of breakpoints on the contour boundary of the second part are identified. Based on the second type of shortest distance from each of the multiple first type breakpoints of the second part to the contour boundary of the first part, multiple second type breakpoints on the contour boundary of the first part are identified; and based on the second type of shortest distance from each of the multiple first type breakpoints of the first part to the contour boundary of the second part, multiple second type breakpoints on the contour boundary of the second part are identified. At the locations of the first type of breakpoint and the second type of breakpoint, the contour boundaries of the first part and the second part are broken to obtain multiple broken boundaries included in the first part and the second part respectively. Among the multiple breaking boundaries of the first part and the second part, the multiple target boundaries that coincide with each other in position are determined as the contact boundaries between the first part and the second part.

2. The method of claim 1, wherein, Identifying multiple convex points of each of the first and second parts from multiple points on the contour boundaries of the first and second parts includes: Multiple reference points are determined from multiple points on the contour boundaries of the first part and the second part, wherein the multiple reference points include at least the inflection points of the contour boundaries; Identify the plurality of convex points from the plurality of reference points.

3. The method of claim 2, wherein, Identifying the plurality of convex points from the plurality of reference points includes performing convex point identification processing on each reference point, the convex point identification processing including: Determine the first adjacent reference point and the second adjacent reference point that are adjacent to the reference point; Calculate the angle between the first vector and the second vector, wherein the first vector represents the vector formed from the first adjacent reference point to the reference point, and the second vector represents the vector formed from the reference point to the second adjacent reference point; Based on the angle between the first vector and the second vector, determine whether the reference point is a convex point.

4. The method of claim 3, wherein, Determining whether the reference point is a convex point based on the angle between the first vector and the second vector includes: If the angle between the first vector and the second vector is greater than or equal to a predetermined angle threshold, the reference point is determined to be a convex point. If the angle between the first vector and the second vector is less than a predetermined angle threshold, the reference point is determined not to be a convex point.

5. The method according to claim 4, wherein, The predetermined angle threshold is 15 degrees.

6. The method of claim 1, wherein, Based on the first type of shortest distance from each of the multiple protrusions of the second part to the contour boundary of the first part, multiple first type breakpoints on the contour boundary of the first part are identified, including: Among the multiple points included in the contour boundary of the first part, the point corresponding to the shortest distance of the first type is determined as the first type of landing point; Based on the first type of shortest distance and a predetermined distance threshold, multiple first type breakpoints are determined from multiple first type landing points on the contour boundary of the first part.

7. The method of claim 1, wherein, Based on the second type of shortest distance from each of the multiple first type breakpoints of the second part to the contour boundary of the first part, the multiple second type breakpoints on the contour boundary of the first part are identified as follows: Among the multiple points included in the contour boundary of the first part, the point corresponding to the shortest distance of the second type is determined as the second type of landing point; Based on the second type of shortest distance and a predetermined distance threshold, multiple second type breakpoints are determined from multiple second type landing points on the contour boundary of the first part.

8. The method according to claim 1, further comprising, after identifying the contact boundaries between the plurality of parts included in the thermal equipment: The heat transfer calculation of the thermal equipment is performed based on a predetermined algorithm, wherein the heat transfer calculation is performed at the contact boundary of the parts of the thermal equipment using a predetermined heat conduction algorithm.

9. The method of claim 8, wherein, The heat transfer calculation is performed during the operation of the thermal equipment, and the method further includes: During the operation of the thermal equipment, it is determined whether the parts of the thermal equipment undergo thermal deformation; When the parts of the thermal equipment undergo thermal deformation, the contour boundaries of the parts of the thermal equipment are updated. Based on the updated contour boundaries, the contact boundaries between multiple parts included in the thermal equipment are re-identified.

10. An apparatus for identifying indirect contact boundaries between parts of a thermal plant, characterized in that, The device is used to identify contact boundaries between multiple parts included in a thermal device, wherein, for any first part and second part that are positionally adjacent to each other among the multiple parts, the device is used to identify the contact boundary between the first part and the second part, the device comprising: The first identification module is used to identify multiple protrusions of the first part and the second part from multiple points on the contour boundary of the first part and the second part respectively. The second identification module is used to identify multiple first-type breakpoints on the contour boundary of the first part based on the first type of shortest distance from each of the multiple protrusions of the second part to the contour boundary of the first part, and to identify multiple first-type breakpoints on the contour boundary of the second part based on the first type of shortest distance from each of the multiple protrusions of the first part to the contour boundary of the second part. The third identification module is used to identify multiple second-type breakpoints on the contour boundary of the first part based on the second-type shortest distance from each of the multiple first-type breakpoints of the second part to the contour boundary of the first part, and to identify multiple second-type breakpoints on the contour boundary of the second part based on the second-type shortest distance from each of the multiple first-type breakpoints of the first part to the contour boundary of the second part. The interruption module is used to interrupt the contour boundaries of the first part and the second part at the positions of the first type of interruption point and the second type of interruption point, so as to obtain multiple interruption boundaries included in the first part and the second part respectively. The determination module is used to determine the multiple target boundaries that coincide with each other in position among the multiple breaking boundaries of the first part and the second part as the contact boundary between the first part and the second part.

Images