Wave load automatic division method based on flow field distribution and related equipment thereof

Abstract

The embodiment of the invention provides a wave load automatic division method based on flow field distribution and related equipment thereof, and the method comprises the steps: obtaining two-dimensional flow field data of a target flow field in a preset horizontal direction, and obtaining coordinate information and flow field speed information of two ends of an initial rod segment based on the two-dimensional flow field data, coordinate information, linear velocity information and flow field velocity information of the midpoint are calculated; under the condition that the midpoint speed information does not meet a preset condition, the initial rod section is divided into two sub-rod sections according to the coordinate information of the midpoint, and midpoint division is conducted on the sub-rod sections according to the coordinate information of the two ends of the sub-rod sections and the flow field speed information; and calculating equivalent point loads at the two ends of the initial rod section according to the coordinate information and the flow field speed information of the multiple division points until the difference value of the flow field speed information and the linear speed information of the midpoint corresponding to each sub rod section meets a preset condition, so that the load error can be reduced, and the design accuracy of the offshore wind power platform is improved.

Description

Technical Field

[0001] The embodiments of the present application relate to the field of electronic information technology, and in particular to a method for automatically dividing wave loads based on flow field distribution and related equipment. Background Art

[0002] In the design stage of offshore platforms, calculation methods based on traditional theories are currently widely used to quickly obtain flow field data. The main structure of an offshore platform is usually a steel frame duct, and in the simulation design process, a rod unit model is usually used for simulation. By applying the Morrison equation load calculation formula, a data scheme for wave loads can be obtained. However, in the above scheme, the grid calculation points of the platform model are discretely distributed in space, which leads to differences between the load sampling of discrete calculation points and the continuously distributed wave load field when wave loads are applied. This difference makes the wave load conditions considered in actual calculations too conservative, which may cause large errors. Summary of the invention

[0003] The embodiments of the present application provide a method for automatically dividing wave loads based on flow field distribution and related equipment, which can automatically divide according to the flow field velocity information and linear velocity information of the midpoint of the rod segment, so as to more accurately and effectively reflect the distribution of wave loads on the structure and reduce errors caused by overly conservative loads.

[0004] To achieve the above-mentioned purpose, a first aspect of an embodiment of the present application provides a method for automatically dividing wave loads based on flow field distribution, comprising: obtaining two-dimensional flow field data of a target flow field in a preset horizontal direction, the two-dimensional flow field data comprising velocity field data and acceleration field data; obtaining coordinate information and flow field velocity information at both ends of an initial rod segment based on the two-dimensional flow field data; calculating coordinate information, linear velocity information and flow field velocity information of a midpoint according to the coordinate information and flow field velocity information at both ends of the initial rod segment; if the difference between the flow field velocity information and the linear velocity information at the midpoint does not meet a preset condition, dividing the initial rod segment into two sub-rod segments according to the coordinate information of the midpoint, and dividing the sub-rod segments by midpoints according to the coordinate information and flow field velocity information at both ends of the sub-rod segments, until the difference between the flow field velocity information and the linear velocity information at the midpoint corresponding to each sub-rod segment meets the preset condition, so as to obtain a plurality of division points corresponding to the initial rod segment; calculating the equivalent point loads at both ends of the initial rod segment according to the coordinate information and flow field velocity information of the plurality of division points.

[0005] In some embodiments, the method of obtaining two-dimensional flow field data of the target flow field in a preset horizontal direction includes: obtaining velocity data of the target flow field, and obtaining a two-dimensional vector of the velocity field and a two-dimensional vector of the acceleration field based on the velocity data; obtaining angle information between a wave horizontal direction of the velocity data and a preset horizontal direction, and obtaining a horizontal rotation matrix based on the angle information; performing vector rotation calculation on the two-dimensional vector of the velocity field and the two-dimensional vector of the acceleration field according to the horizontal rotation matrix, to obtain the two-dimensional flow field data of the target flow field in the preset horizontal direction.

[0006] In some embodiments, the flow field velocity information includes a two-dimensional velocity component and a two-dimensional acceleration component at both ends of the initial rod segment, and obtaining the coordinate information and flow field velocity information at both ends of the initial rod segment based on the two-dimensional flow field data includes: obtaining the coordinate information of the two ends of the initial rod segment in the plane corresponding to the preset horizontal direction based on the two-dimensional flow field data; obtaining the two-dimensional velocity component at both ends of the initial rod segment based on the velocity field data, and obtaining the two-dimensional acceleration component at both ends of the initial rod segment based on the acceleration field data.

[0007] In some embodiments, the coordinate information, linear velocity information and flow field velocity information of the midpoint are calculated based on the coordinate information and flow field velocity information at both ends of the initial rod segment, including: calculating the coordinate information of the corresponding midpoint based on the coordinate information at both ends of the initial rod segment; calculating the linear velocity and linear acceleration of the midpoint based on the two-dimensional velocity components and two-dimensional acceleration components at both ends of the initial rod segment; and acquiring the flow field velocity information of the midpoint based on the two-dimensional flow field data and the coordinate information of the midpoint.

[0008] In some embodiments, the preset conditions include: the difference between the flow field velocity information and the linear velocity information of the midpoint is less than the absolute value of the linear velocity information multiplied by a preset ratio; or, the difference between the flow field velocity information and the linear velocity information of the midpoint is less than a preset unit velocity threshold.

[0009] In some embodiments, the initial rod segment is divided into two sub-rod segments according to the coordinate information of the midpoint, and the midpoint division of the sub-rod segments is performed according to the coordinate information of the two ends of the sub-rod segments and the flow field velocity information, including: dividing the initial rod segment into a first sub-rod segment corresponding to the starting end of the initial rod segment to the midpoint, and a second sub-rod segment corresponding to the midpoint to the end of the initial rod segment according to the coordinate information of the midpoint; calculating the linear velocity information and flow field velocity information of the midpoint of the first sub-rod segment according to the coordinate information of the starting end of the initial rod segment and the midpoint, the coordinate information and the flow field velocity information. degree information, and when the difference between the linear velocity information and the flow field velocity information of the midpoint of the first sub-rod segment does not meet the preset condition, repeatedly divide the first sub-rod segment at the midpoint to obtain a plurality of sub-rod segments; calculate the linear velocity information and the flow field velocity information of the midpoint of the second sub-rod segment according to the coordinate information, coordinate information and flow field velocity information of the midpoint of the initial segment, and when the difference between the linear velocity information and the flow field velocity information of the midpoint of the second sub-rod segment does not meet the preset condition, repeatedly divide the second sub-rod segment at the midpoint to obtain a plurality of sub-rod segments.

[0010] In some embodiments, obtaining the multiple division points corresponding to the initial rod segment includes: in the midpoint division process, determining the multiple midpoints whose differences do not meet the preset conditions as the multiple division points corresponding to the initial rod segment.

[0011] In some embodiments, the equivalent point loads at both ends of the initial rod segment are calculated based on the coordinate information and flow field velocity information of the multiple division points, including: calculating the wave loads of the multiple division points based on the coordinate information and flow field velocity information of the multiple division points through the Morrison equation; performing integral conversion on the wave loads of the multiple division points based on the principle of virtual work to obtain the equivalent point loads at both ends of the initial rod segment.

[0012] To achieve the above-mentioned purpose, the second aspect of the embodiments of the present application proposes a computer-readable storage medium storing computer-executable instructions, and the computer-executable instructions are used to execute the method for automatic division of wave loads based on flow field distribution as described in any embodiment of the first aspect.

[0013] To achieve the above-mentioned purpose, a third aspect of the embodiments of the present application provides a computer program product, which includes a computer program stored in a computer-readable storage medium. A processor of a computer device reads the computer program from the computer-readable storage medium, and the processor executes the computer program, so that the computer device executes the method for automatically dividing wave loads based on flow field distribution as described in any one of the embodiments of the first aspect.

[0014] The automatic wave load division method based on flow field distribution and related equipment provided by the embodiment scheme of the present application can obtain two-dimensional flow field data including velocity field and acceleration field according to wave theory and direction, and these data can be used to simulate flow field conditions in different horizontal directions; further, based on the two-dimensional flow field data, the position coordinates and flow field velocity information of the beginning and end of the rod segment are obtained; further, according to the coordinates and velocity information of the beginning and end, the coordinates, linear velocity and actual flow field velocity of the midpoint of the rod segment are calculated; further, the linear velocity and actual flow field velocity of the midpoint are compared, and if the difference between the linear velocity and the actual flow field velocity of the midpoint does not meet the preset condition, the rod segment is further divided into two sub-rod segments from the midpoint; further, the above process is repeated for each newly formed sub-rod segment until the difference between the midpoint flow field velocity and the linear velocity of all sub-rod segments meets the preset condition. This ensures that the flow field velocity at each division point is close enough to the velocity distribution of the actual flow field; further, once the division points of all sub-rod segments are determined, the wave loads can be calculated using the Morrison equation based on the coordinate information and flow field velocity information of these points, and these distributed loads can be converted into equivalent point loads acting on both ends of the rod segment through the principle of virtual work; it can be understood that the present application is able to dynamically adjust the division of the rod segments to ensure that the flow field distribution at any given position can be accurately represented, and through recursive division until the flow field velocity distribution of all sub-segments is consistent with the linear approximation, it can avoid conservative estimates caused by the difference between discrete point load sampling and the continuously distributed wave load field, and thus can more accurately reflect the distribution of wave loads on the structure, reduce unnecessary errors, and improve the accuracy of offshore wind power platform design. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] Figure 1 A flow chart of a method for automatically dividing wave loads based on flow field distribution provided in one embodiment of the present application;

[0016] Figure 2 A flow chart of a method for automatically dividing wave loads based on flow field distribution provided in one embodiment of the present application;

[0017] Figure 3 A flow chart of a method for automatically dividing wave loads based on flow field distribution provided in one embodiment of the present application;

[0018] Figure 4 A flow chart of a method for automatically dividing wave loads based on flow field distribution provided in one embodiment of the present application;

[0019] Figure 5 A flow chart of a method for automatically dividing wave loads based on flow field distribution provided in one embodiment of the present application;

[0020] Figure 6A schematic diagram of a method for automatically dividing wave loads based on flow field distribution provided in one embodiment of the present application;

[0021] Figure 7 A flow chart of a method for automatically dividing wave loads based on flow field distribution provided in one embodiment of the present application;

[0022] Figure 8 A flow chart of a method for automatically dividing wave loads based on flow field distribution provided in one embodiment of the present application;

[0023] Fig. 9 An example flow chart of a method for automatically dividing wave loads based on flow field distribution provided in one embodiment of the present application;

[0024] Fig.10 A schematic diagram of the structure of an electronic device provided in one embodiment of the present application. DETAILED DESCRIPTION

[0025] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0026] In some embodiments, although the functional units are divided in the system schematic diagram and the logical order is shown in the flowchart, in some cases, the steps shown or described may be performed in a different order than the unit division in the system or the order in the flowchart. The terms first, second, etc. in the specification, claims, and the above drawings are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence.

[0027] In addition, unless otherwise clearly specified and limited, the term "connected" should be understood in a broad sense. For example, it can be a fixed connection or a movable connection, a detachable connection or a non-detachable connection, or an integral connection; it can be a mechanical connection, an electrical connection, or the two can communicate with each other; it can be a direct connection or an indirect connection through an intermediate medium.

[0028] In the description of the embodiments of the present application, the description of reference terms "one embodiment / implementation", "another embodiment / implementation" or "certain embodiments / implementations", "in the above embodiments / implementations", etc. means that the specific features, structures, materials or characteristics described in combination with the implementation or example are included in at least two embodiments or implementations disclosed in the present application. In the disclosure of this application, the schematic representation of the above terms does not necessarily refer to the same embodiment or implementation. It should be noted that although the logical order is shown in the flowchart, in some cases, the steps shown or described may be performed in an order different from that in the flowchart.

[0029] In the design stage of offshore platforms, calculation methods based on traditional theories are currently widely used to quickly obtain flow field data. The main structure of an offshore platform is usually a steel frame duct, and in the simulation design process, a rod unit model is usually used for simulation. By applying the Morrison equation load calculation formula, a data scheme for wave loads can be obtained. However, in the above scheme, the grid calculation points of the platform model are discretely distributed in space, which leads to differences between the load sampling of discrete calculation points and the continuously distributed wave load field when wave loads are applied. This difference makes the wave load conditions considered in actual calculations too conservative, which may cause large errors.

[0030] Based on this, an embodiment of the present application provides a method for automatically dividing wave loads based on flow field distribution and related equipment. The method includes obtaining two-dimensional flow field data of the target flow field in a preset horizontal direction, and obtaining coordinate information and flow field velocity information at both ends of the initial rod segment based on the two-dimensional flow field data to calculate the coordinate information, linear velocity information and flow field velocity information of the midpoint; when the midpoint velocity information does not meet the preset conditions, the initial rod segment is divided into two sub-rod segments according to the coordinate information of the midpoint, and the sub-rod segments are midpoint-divided according to the coordinate information and flow field velocity information at both ends of the sub-rod segments, until the difference between the flow field velocity information and the linear velocity information of the midpoint corresponding to each sub-rod segment meets the preset conditions, so as to calculate the equivalent point loads at both ends of the initial rod segment according to the coordinate information and flow field velocity information of multiple division points, which can reduce load errors and improve the accuracy of offshore wind power platform design.

[0031] The embodiments of the present application are further described below in conjunction with the accompanying drawings.

[0032] refer to Figure 1 , Figure 1 A flowchart of a method for automatically dividing wave loads based on flow field distribution provided in an embodiment of the present application; in some embodiments, the present application discloses a method for automatically dividing wave loads based on flow field distribution, including but not limited to the following steps:

[0033] Step S110, acquiring two-dimensional flow field data of a target flow field in a preset horizontal direction, where the two-dimensional flow field data includes velocity field data and acceleration field data;

[0034] Among them, this step refers to the fluid motion state data obtained through theoretical calculations or other means based on the propagation direction and properties of the waves (such as wave speed, wave height, etc.). These data describe the changes in velocity and acceleration with position in two-dimensional space (usually in a plane perpendicular to the wave propagation direction).

[0035] Step S120, obtaining coordinate information and flow field velocity information of both ends of the initial rod segment based on the two-dimensional flow field data;

[0036] Among them, the coordinate information and flow field velocity information at both ends of the initial rod segment are the positions of the beginning and end of a rod segment in the offshore platform model in space, as well as the flow field velocity at these positions, in preparation for subsequent flow field distribution calculations.

[0037] Step S130, calculating the coordinate information, linear velocity information and flow field velocity information of the midpoint according to the coordinate information and flow field velocity information of both ends of the initial rod segment;

[0038] Among them, the position of the midpoint between the beginning and the end of the rod segment is calculated, and the linear velocity of the midpoint is deduced based on the velocity information at both ends of the rod segment. In addition, since the two-dimensional flow field data includes the flow field velocity information of each position of any rod segment in the offshore platform model, the flow field velocity information of the two points and the midpoint can be directly obtained from the two-dimensional flow field data.

[0039] Step S140, when the difference between the flow field velocity information and the linear velocity information of the midpoint does not meet the preset condition, the initial rod segment is divided into two sub-rod segments according to the coordinate information of the midpoint, and the sub-rod segments are divided at the midpoint according to the coordinate information and the flow field velocity information at both ends of the sub-rod segment, until the difference between the flow field velocity information and the linear velocity information of the midpoint corresponding to each sub-rod segment meets the preset condition, so as to obtain a plurality of division points corresponding to the initial rod segment;

[0040] This step describes a recursive process to ensure the accuracy of the calculation. If the difference between the flow field velocity and the linear velocity at the midpoint is too large, that is, it does not meet the pre-set conditions, the rod segment will be further subdivided until the difference between the flow field velocity and the linear velocity of all subdivided segments meets the set conditions. This is done to ensure the accuracy of the calculation results and reduce errors.

[0041] Step S150, calculating the equivalent point loads at both ends of the initial rod segment according to the coordinate information of the plurality of division points and the flow field velocity information.

[0042] Among them, this step is to convert the flow field velocity information at each division point into equivalent point loads acting on both ends of the rod segment. This involves the application of the Morrison equation, which is used to calculate the force exerted on the structure in the fluid. In this way, the distributed wave load can be simplified to a concentrated force acting on both ends of the rod segment, thereby simplifying the structural analysis.

[0043] In some embodiments, the present application can obtain two-dimensional flow field data including velocity field and acceleration field according to wave theory and direction, and these data can be used to simulate flow field conditions in different horizontal directions; further, based on the two-dimensional flow field data, the position coordinates and flow field velocity information of the beginning and end of the rod segment are obtained; further, according to the coordinates and velocity information of the beginning and end, the coordinates, linear velocity and actual flow field velocity of the midpoint of the rod segment are calculated; further, the linear velocity and actual flow field velocity of the midpoint are compared, and if the difference between the two does not meet the preset conditions, the rod segment is further divided into two sub-rod segments from the midpoint; further, the above process is repeated for each newly formed sub-rod segment until the difference between the midpoint flow field velocity and the linear velocity of all sub-rod segments meets the preset conditions. This ensures that the flow field velocity at each division point is close enough to the velocity distribution of the actual flow field; further, once the division points of all sub-rod segments are determined, the wave loads can be calculated using the Morrison equation based on the coordinate information and flow field velocity information of these points, and these distributed loads can be converted into equivalent point loads acting on both ends of the rod segment through the principle of virtual work; it can be understood that the present application is able to dynamically adjust the division of the rod segments to ensure that the flow field distribution at any given position can be accurately represented, and through recursive division until the flow field velocity distribution of all sub-segments is consistent with the linear approximation, it can avoid conservative estimates caused by the difference between discrete point load sampling and the continuously distributed wave load field, and thus can more accurately reflect the distribution of wave loads on the structure, reduce unnecessary errors, and improve the accuracy of offshore wind power platform design.

[0044] refer to Figure 2 , Figure 2 A flowchart of a method for automatically dividing wave loads based on flow field distribution provided in an embodiment of the present application; in some embodiments, obtaining two-dimensional flow field data of a target flow field in a preset horizontal direction includes but is not limited to the following steps:

[0045] Step S210, acquiring velocity data of the target flow field, and obtaining a velocity field two-dimensional vector and an acceleration field two-dimensional vector according to the velocity data;

[0046] Step S220, obtaining angle information between the wave horizontal direction of the velocity data and the preset horizontal direction, and obtaining a horizontal rotation matrix according to the angle information;

[0047] Step S230, performing vector rotation calculation on the velocity field two-dimensional vector and the acceleration field two-dimensional vector according to the horizontal rotation matrix to obtain two-dimensional flow field data of the target flow field in a preset horizontal direction.

[0048] In some embodiments, the velocity data of the target flow field refers to the velocity of the fluid at different positions under a certain wave condition, the two-dimensional vector of the velocity field represents the magnitude and direction of the fluid velocity on the two-dimensional plane, the two-dimensional vector of the acceleration field represents the magnitude and direction of the fluid acceleration on the two-dimensional plane, and the angle information is the angle difference between the horizontal direction of the wave and the preset horizontal direction. The above-mentioned horizontal rotation matrix is ​​used to multiply the two-dimensional vector of the velocity field and the two-dimensional vector of the acceleration field to adjust the directions of these vectors to conform to the velocity field and acceleration field data after adjustment in the preset horizontal direction. These data represent the velocity and acceleration of the fluid in the specified direction, thereby providing accurate data support for subsequent wave load calculations.

[0049] Specifically, flow field data is obtained according to traditional wave theory and wave direction. The flow field data includes velocity field and acceleration field data. In steps S210 to S230, further, the velocity field data and the acceleration field data can be expressed in the form of a two-dimensional plane vector (assuming that the plane is located in the XOZ plane), which are velocity field and the acceleration field where u x0 、u z0 and a x0 、a z0 are the horizontal and vertical components of the velocity field and acceleration field, respectively. By applying the horizontal rotation matrix to the above two-dimensional wave field, the flow fields in different horizontal directions can be considered. Assuming that the angle between the horizontal direction vector of the wave and the global horizontal x-axis is α (counterclockwise rotation is a positive angle), the rotation matrix R is:

[0050]

[0051] New velocity field

[0052] New acceleration field

[0053] refer to Figure 3 , Figure 3 A flow chart of a method for automatically dividing wave loads based on flow field distribution provided in an embodiment of the present application; in some embodiments, the flow field velocity information includes a two-dimensional velocity component and a two-dimensional acceleration component at both ends of the initial rod segment, and the coordinate information and flow field velocity information at both ends of the initial rod segment are obtained based on the two-dimensional flow field data, including but not limited to the following steps:

[0054] Step S310, obtaining coordinate information of both ends of the initial rod segment in a plane corresponding to a preset horizontal direction based on the two-dimensional flow field data;

[0055] Step S320, obtaining two-dimensional velocity components at both ends of the initial rod segment based on the velocity field data, and obtaining two-dimensional acceleration components at both ends of the initial rod segment based on the acceleration field data.

[0056] In some embodiments, the position coordinates of the two ends of the initial rod segment in a specified horizontal direction (i.e., the wave propagation direction) are determined based on the information in the two-dimensional flow field data. These coordinates are obtained after the rotation operation to ensure that they are consistent with the direction of the waves, and then the velocity information of the two ends of the initial rod segment can be extracted from the velocity field data, including the velocity components in the x-direction and the y-direction.

[0057] Specifically, in steps S310 to S320, the velocity and acceleration of the start end A and the end end B of each rod segment are further considered and recorded as and Where (x A ,z A ) and (x B ,z B ) are the local coordinates of the two end points of the rod segment in the horizontal wave plane obtained by the rotation R in S1, u cX and a cX (X=A or B, c=x or z) are the magnitudes of the velocity component and the acceleration component respectively.

[0058] refer to Figure 4 , Figure 4 A flowchart of a method for automatically dividing wave loads based on flow field distribution provided in an embodiment of the present application; in some embodiments, the coordinate information, linear velocity information and flow field velocity information of the midpoint are calculated according to the coordinate information of both ends of the initial rod segment and the flow field velocity information, including but not limited to the following steps:

[0059] Step S410, calculating the coordinate information of the corresponding midpoint according to the coordinate information of both ends of the initial rod segment;

[0060] Step S420, calculating the linear velocity and linear acceleration of the midpoint according to the two-dimensional velocity components and the two-dimensional acceleration components at both ends of the initial rod segment;

[0061] Step S430, acquiring flow field velocity information of the midpoint based on the two-dimensional flow field data and the coordinate information of the midpoint.

[0062] In some embodiments, the coordinate information of the two ends of the initial rod segment is used to calculate the coordinates of the midpoint between the two points. This can usually be achieved through a simple geometric average, based on the velocity and acceleration components at the two ends of the initial rod segment, using a linear interpolation method to estimate the velocity and acceleration of the midpoint. Similarly, the acceleration is calculated in a similar way; then, based on the coordinate information of the midpoint, the actual flow field velocity information at the midpoint can be found or calculated in the two-dimensional flow field data.

[0063] In some embodiments, the preset conditions include: the difference between the flow field velocity information and the linear velocity information of the midpoint is less than the absolute value of the linear velocity information multiplied by the preset ratio; or, the difference between the flow field velocity information and the linear velocity information of the midpoint is less than the preset unit velocity threshold, wherein the preset conditions are used to determine the consistency between the linear velocity information and the flow field velocity information of the midpoint: the difference between the flow field velocity information and the linear velocity information of the midpoint is less than the absolute value of the linear velocity information multiplied by the preset ratio: this means that the difference between the flow field velocity and the linear velocity at the midpoint cannot exceed a certain ratio of the linear velocity value. For example, if the preset ratio is 5%, then the difference must be less than 5% of the linear velocity value; the difference between the flow field velocity information and the linear velocity information of the midpoint is less than the preset unit velocity threshold: this means that the difference between the flow field velocity and the linear velocity at the midpoint must be less than a specific numerical threshold, such as 0.1 m / s. Through such preset conditions, it can be ensured that the velocity information of the midpoint in the calculation process is close enough to the real flow field velocity information, thereby improving the accuracy of the calculation.

[0064] Specifically, calculate the midpoint M between the starting point A and the end point B 0,0 Linear speed and linear acceleration At the same time, according to the midpoint M 0,0 The actual speed is obtained from the spatial coordinates and acceleration The midpoint M 0,0 The spatial coordinates of

[0065]

[0066] Midpoint M 0,0 The linear velocity is:

[0067]

[0068] Midpoint M 0,0 The linear acceleration is:

[0069]

[0070] Midpoint M 0,0 The flow field velocity is:

[0071]

[0072] Midpoint M 0,0 The flow field acceleration is:

[0073]

[0074] refer to Figure 5 , Figure 5 A flow chart of a method for automatically dividing wave loads based on flow field distribution provided in an embodiment of the present application; in some embodiments, the initial rod segment is divided into two sub-rod segments according to the coordinate information of the midpoint, and the sub-rod segments are divided at the midpoint according to the coordinate information at both ends of the sub-rod segment and the flow field velocity information, including but not limited to the following steps:

[0075] Step S510, dividing the initial rod segment into a first sub-rod segment corresponding to the initial rod segment start end to the midpoint and a second sub-rod segment corresponding to the midpoint to the initial rod segment end according to the coordinate information of the midpoint;

[0076] Step S520, calculating the linear velocity information and the flow field velocity information of the midpoint of the first sub-rod segment according to the coordinate information, coordinate information and flow field velocity information of the starting end and the midpoint of the initial rod segment, and when the difference between the linear velocity information and the flow field velocity information of the midpoint of the first sub-rod segment does not meet the preset condition, repeatedly dividing the midpoint of the first sub-rod segment to obtain a plurality of sub-rod segments;

[0077] Step S530, calculate the linear velocity information and flow field velocity information of the midpoint of the second sub-rod segment based on the coordinate information, coordinate information and flow field velocity information of the midpoint and the end of the initial rod segment, and when the difference between the linear velocity information and the flow field velocity information of the midpoint of the second sub-rod segment does not meet the preset conditions, repeat the midpoint division of the second sub-rod segment to obtain multiple sub-rod segments.

[0078] In some embodiments, based on the calculated midpoint coordinate information, the initial rod segment is divided into two to form two new sub-segments, the first sub-segment extends from the starting end of the initial rod segment to the midpoint, and the second sub-segment extends from the midpoint to the end of the initial rod segment. Based on the coordinate information and flow field velocity information at both ends of the first sub-segment, the linear velocity and flow field velocity of the midpoint of the sub-segment are calculated. If the difference between the linear velocity and the flow field velocity of the midpoint of the first sub-segment does not meet the preset conditions, the first sub-segment is further divided into two smaller sub-segments, and the same calculation and division process is repeated for these new sub-segments; based on the coordinate information and flow field velocity information at both ends of the second sub-segment, the linear velocity and flow field velocity of the midpoint of the sub-segment are calculated. If the difference between the linear velocity and the flow field velocity of the midpoint of the second sub-segment does not meet the preset conditions, the second sub-segment is further divided into two smaller sub-segments, and the same calculation and division process is repeated for these new sub-segments.

[0079] Specifically, the midpoint M 0,0 The linear flow field value and the actual flow field value are compared and judged according to two classification criteria. The two criteria are: criterion 1: the absolute value of the difference in flow field data is less than 5% of the absolute value of the linear flow field value. Criterion 2: the absolute value of the difference in flow field data is less than 0.1 unit (for velocity units, m / s; for acceleration units, m / s 2 ).

[0080] Criterion 1 requires:

[0081]

[0082] and

[0083]

[0084] Criterion 2 requires:

[0085]

[0086] and

[0087]

[0088] Wherein i, j are integers greater than or equal to 0. Through the above recursive process, the present application ensures that the difference between the linear velocity and the flow field velocity at the midpoint of each sub-rod segment is controlled within an acceptable range, thereby reducing the error caused by the difference between discrete point load sampling and the continuously distributed wave load field.

[0089] Through the above recursive process, the present application ensures that the difference between the linear velocity and the flow field velocity of each sub-rod segment midpoint is controlled within an acceptable range, thereby reducing the error caused by the difference between discrete point load sampling and the continuously distributed wave load field.

[0090] refer to Figure 6 , Figure 6 Schematic diagram of a wave load automatic division method based on flow field distribution provided in an embodiment of the present application; it can be understood that when i=0, j=0, it corresponds to the first automatic division. If the criterion 1 or criterion 2 is met, the midpoint M is returned. i,j Otherwise, save the midpoint M i,j The actual flow field value and Calculate the midpoint M i,j The midpoint of the left sub-bar segment M i+1,j and the midpoint M of the right sub-bar segment i,j+1, follow the above steps to calculate the linear flow field and actual flow field of each sub-segment midpoint to determine whether they meet criterion one or criterion two. If the sub-segment meets any criterion, it will not be divided further; if not, it is necessary to further subdivide the segment into two sub-segments and record the flow field value and coordinates of the midpoint. Repeat this process until all sub-segments meet the corresponding criterion, such as Figure 6 As shown, the flow field values ​​and coordinate information of all recorded points are finally returned.

[0091] refer to Figure 7 , Figure 7 A flowchart of a method for automatically dividing wave loads based on flow field distribution provided in an embodiment of the present application; in some embodiments, obtaining multiple division points corresponding to the initial rod segment includes but is not limited to the following steps:

[0092] Step S710, in the midpoint division process, a plurality of midpoints whose difference values ​​do not meet a preset condition are determined as a plurality of division points corresponding to the initial pole segment.

[0093] In some embodiments, when performing midpoint division, if the difference between the flow field velocity (or acceleration) at a midpoint and the flow field velocity (or acceleration) obtained by linear interpolation does not meet the preset conditions (i.e., criterion one: the absolute value of the difference is less than 5% of the absolute value of the linear flow field value, or criterion two: the absolute value of the difference is less than 0.1 units), then this midpoint will be marked as a division point. In other words, if there is a significant difference between the flow field value at the midpoint and the estimated value obtained by linear interpolation, it indicates that the current division is not fine enough and cannot accurately reflect the actual distribution of the flow field. At this time, the midpoint is determined as a position that needs to be further divided, and the algorithm will continue to divide these sub-bar segments that do not meet the conditions into smaller parts again, and recalculate the difference between the flow field value at the new midpoint and the flow field value of the linear interpolation until the midpoints of all sub-bar segments meet the preset conditions. This process ensures that each division point can represent the actual situation of the flow field distribution, thereby improving the accuracy of the overall calculation, and ensuring that the equivalent point load of the division point finally obtained can more accurately represent the effect of the actual wave load.

[0094] It can be understood that in the midpoint division process, the midpoints whose corresponding differences meet the preset conditions will not be used for rod segment division, that is, the midpoints whose differences meet the preset conditions will not be determined as the multiple division points corresponding to the initial rod segment. The initial rod segment or sub-rod segment whose midpoint difference meets the preset conditions can effectively approximate the actual flow field distribution.

[0095] refer to Figure 8 , Figure 8A flowchart of a method for automatically dividing wave loads based on flow field distribution provided in an embodiment of the present application; in some embodiments, the equivalent point loads at both ends of the initial rod segment are calculated based on the coordinate information of multiple division points and the flow field velocity information, including but not limited to the following steps:

[0096] Step S810, calculating the wave loads of the plurality of divided points according to the coordinate information and flow field velocity information of the plurality of divided points by using the Morrison equation;

[0097] Step S820, performing integral conversion on the wave loads at multiple division points based on the principle of virtual work to obtain equivalent point loads at both ends of the initial rod segment.

[0098] Here, the Morrison equation is written as:

[0099]

[0100] Where D is the outer diameter of the rod segment, ρ is the density of seawater, and Substitute the Morrison equation into the flow field of all saved points in sequence to obtain a one-dimensional rod distribution force vector array. The vector array constitutes a multilinear piecewise function as follows:

[0101]

[0102] Where N is the total number of points saved in the segment, i = 1,…,N,l i is the distance from the starting point of the segment to the i-th point, and the length of the segment is L, then l N =L, j = 2, ..., N, δ ij is the discrete Dirac function, that is:

[0103]

[0104] According to the principle of virtual work, the point loads at both ends of the rod segment and It can be calculated by the following formula:

[0105]

[0106] refer to Fig. 9 , Fig. 9 An example flow chart of a method for automatically dividing wave loads based on flow field distribution provided in an embodiment of the present application; in some embodiments, the present application may include the following steps:

[0107] Input velocity field and acceleration field: First, two-dimensional velocity field and acceleration field data for simulating the wave environment need to be provided. These data can come from physical experiments, numerical simulations or other reliable sources.

[0108] Traversing the segments in the simulation model that are subject to wave loads: In this stage, the segments to be studied are selected and placed in a given flow environment. This step usually involves placing the segments at specific positions and angles in order to simulate their behavior in actual sea conditions.

[0109] Calculate and record the velocity and acceleration at the beginning and end of the rod segment: Use the input flow field data to calculate the velocity and acceleration at the beginning and end of the rod segment. These values ​​reflect the state of the rod segment under the action of waves.

[0110] Calculate linear velocity, linear acceleration, and flow velocity and flow acceleration at the midpoint of the segment: Next, calculate the velocity and acceleration at the midpoint of the segment, taking into account the effects of the flow field. This is done by linear interpolation, which applies the velocity and acceleration information at the start and end points to the midpoint.

[0111] Determine whether the midpoint velocity meets the criteria: Check whether the difference between the linear velocity of the midpoint and the flow field velocity meets certain criteria. If the difference between the two is small, it means that the linear model of the midpoint is accurate enough and you can go directly to the next step. Otherwise, you need to continue to refine the division of the rod segment.

[0112] Output the flow field values ​​of all divided points: If the difference between the linear velocity and the flow field velocity of the midpoint is large, the rod segment needs to be further subdivided into more sub-segments, and the same operation is performed on the midpoint of each sub-segment until all midpoints meet the predetermined criteria. Then, output the flow field values ​​of all divided points.

[0113] Calculate wave loads according to Morrison's equation: Once all midpoints meet the requirements, you can start calculating wave loads. Morrison's equation is a commonly used mathematical model for calculating fluid dynamics loads, which takes into account the inertial force and resistance of an object in a fluid.

[0114] Calculate the equivalent point loads at both ends of the rod segment according to the strain energy equivalence principle: Finally, integrate all the calculated wave loads, and convert the loads distributed on the rod segment into equivalent point loads at both ends of the rod segment according to the strain energy equivalence principle. This method helps to simplify the analysis, making it easier for engineers to understand and handle complex load situations.

[0115] When the entire segment is divided into enough sub-segments and the midpoints of all sub-segments meet the requirements, the equivalent point loads at both ends of the segment can be obtained, which can be used to evaluate the strength and stability of the structure.

[0116] Some embodiments of the present application provide an electronic device, Fig.10 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application, referring to Fig.10The electronic device includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the method for automatically dividing wave loads based on flow field distribution according to any one of the above embodiments is implemented. For example, the method described above is executed. Figure 1 Steps S110 to S150 of the method, Figure 2 Steps S210 to S230 of the method, Figure 3 Steps S310 to S320 of the method, Figure 4 Steps S410 to S430 of the method, Figure 5 Steps S510 to S530 of the method, Figure 7 In the method step S710, Figure 8 Method steps S810 to S820.

[0117] The electronic device 1000 of the embodiment of the present application includes one or more processors 1010 and a memory 1020. Fig.10 In the figure, a processor 1010 and a memory 1020 are taken as an example.

[0118] The processor 1010 and the memory 1020 may be connected via a bus or other means. Fig.10 The example of connecting through bus is taken in the following.

[0119] The memory 1020, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs and non-transitory computer executable programs. In addition, the memory 1020 may include a high-speed random access memory, and may also include a non-transitory memory, such as at least one disk storage device, a flash memory device, or other non-transitory solid-state storage device. In some embodiments, the memory 1020 may optionally include a memory 1020 remotely disposed relative to the processor 1010, and these remote memories may be connected to the electronic device 1000 via a network. At the same time, examples of the above-mentioned network include but are not limited to the Internet, an intranet, a local area network, a mobile communication network, and a combination thereof.

[0120] In some embodiments, when the processor executes the computer program, it executes the method for automatic division of wave loads based on flow field distribution of any of the above embodiments at preset intervals.

[0121] Those skilled in the art will understand that Fig.10 The device structure shown in the figure does not constitute a limitation on the electronic device 1000, and may include more or less components than shown in the figure, or combine certain components, or arrange the components differently.

[0122] exist Fig.10In the electronic device 1000 shown, the processor 1010 can be used to call the automatic wave load division method based on flow field distribution stored in the memory 1020, so as to implement the automatic wave load division method based on flow field distribution.

[0123] Based on the hardware structure of the above-mentioned electronic device 1000, various embodiments of the multi-agent interaction system of the present application are proposed. At the same time, the non-transient software programs and instructions required to implement the automatic wave load division method based on flow field distribution of the above-mentioned embodiment are stored in the memory. When executed by the processor, the automatic wave load division method based on flow field distribution of the above-mentioned embodiment is executed.

[0124] The embodiment of the present application also provides a computer program product, which includes a computer program, which is stored in a computer-readable storage medium. The processor of the computer device reads the computer program from the computer-readable storage medium, and the processor executes the computer program, so that the computer device executes the method for automatically dividing wave loads based on flow field distribution according to any of the above embodiments, for example, executing the above described Figure 1 Steps S110 to S150 of the method, Figure 2 Steps S210 to S230 of the method, Figure 3 Steps S310 to S320 of the method, Figure 4 Steps S410 to S430 of the method, Figure 5 Steps S510 to S530 of the method, Figure 7 In the method step S710, Figure 8 Method steps S810 to S820.

[0125] The device embodiments described above are merely illustrative, and the units described as separate components may or may not be physically separated, that is, they may be located in one place, or they may be distributed on multiple network nodes. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0126] It will be appreciated by those skilled in the art that all or some of the steps and systems in the disclosed method above can be implemented as software, firmware, hardware and appropriate combinations thereof. Some physical components or all physical components can be implemented as software executed by a processor, such as a central processing unit, a digital signal processor or a microprocessor, or implemented as hardware, or implemented as an integrated circuit, such as an application-specific integrated circuit. Such software can be distributed on a computer-readable medium, and the computer-readable medium can include a computer-readable storage medium (or a non-transitory medium) and a communication medium (or a temporary medium). As known to those skilled in the art, the term computer-readable storage medium includes volatile and non-volatile, removable and non-removable media implemented in any method or technology for storing information (such as computer-readable instructions, data structures, program units or other data). Computer-readable storage media include, but are not limited to, RAM, ROM, EEPROM, flash memory or other memory technologies, CD to ROM, digital versatile disk (DVD) or other optical disk storage, magnetic cassettes, magnetic tapes, disk storage or other magnetic storage devices, or any other medium that can be used to store desired information and can be accessed by a computer. Furthermore, it is well known to those skilled in the art that communication media generally contain computer readable instructions, data structures, program elements, or other data in a modulated data signal such as a carrier wave or other transport mechanism, and may include any information delivery media.

[0127] The above is a specific description of the preferred implementation of the present application, but the present application is not limited to the above implementation mode. Technical personnel familiar with the field can also make various equivalent modifications or substitutions without violating the spirit of the present application. These equivalent modifications or substitutions are all included in the scope defined by the claims of the present application.

Claims

1. A method for automatic division of wave loads based on flow field distribution, characterized in that: include: Acquire two-dimensional flow field data of a target flow field in a preset horizontal direction, wherein the two-dimensional flow field data includes velocity field data and acceleration field data; Based on the two-dimensional flow field data, coordinate information and flow field velocity information of both ends of the initial rod segment are obtained; Calculate the coordinate information, linear velocity information and flow field velocity information of the midpoint according to the coordinate information and flow field velocity information of both ends of the initial rod segment; In the case where the difference between the flow field velocity information and the linear velocity information of the midpoint does not meet the preset condition, the initial rod segment is divided into two sub-rod segments according to the coordinate information of the midpoint, and the sub-rod segments are divided at the midpoint according to the coordinate information and the flow field velocity information at both ends of the sub-rod segment, until the difference between the flow field velocity information and the linear velocity information of the midpoint corresponding to each sub-rod segment meets the preset condition, so as to obtain a plurality of division points corresponding to the initial rod segment; The equivalent point loads at both ends of the initial rod segment are calculated according to the coordinate information of the multiple division points and the flow field velocity information.

2. The method for automatic division of wave loads based on flow field distribution according to claim 1 is characterized in that: The step of obtaining two-dimensional flow field data of a target flow field in a preset horizontal direction includes: Acquire velocity data of the target flow field, and obtain a velocity field two-dimensional vector and an acceleration field two-dimensional vector according to the velocity data; Acquire the angle information between the wave horizontal direction of the velocity data and the preset horizontal direction, and obtain the horizontal rotation matrix according to the angle information; Vector rotation calculation is performed on the two-dimensional vector of the velocity field and the two-dimensional vector of the acceleration field according to the horizontal rotation matrix to obtain two-dimensional flow field data of the target flow field in the preset horizontal direction.

3. The method for automatic division of wave loads based on flow field distribution according to claim 1 is characterized in that: The flow field velocity information includes a two-dimensional velocity component and a two-dimensional acceleration component at both ends of the initial rod segment. The coordinate information and flow field velocity information at both ends of the initial rod segment are obtained based on the two-dimensional flow field data, including: Based on the two-dimensional flow field data, coordinate information of both ends of the initial rod segment in the plane corresponding to the preset horizontal direction is obtained; The two-dimensional velocity components at both ends of the initial rod segment are obtained based on the velocity field data, and the two-dimensional acceleration components at both ends of the initial rod segment are obtained based on the acceleration field data.

4. The method for automatic division of wave loads based on flow field distribution according to claim 3 is characterized in that: The calculating of the coordinate information, linear velocity information and flow field velocity information of the midpoint according to the coordinate information and flow field velocity information of the two ends of the initial rod segment includes: Calculate the coordinate information of the corresponding midpoint according to the coordinate information of both ends of the initial rod segment; Calculating the linear velocity and linear acceleration of the midpoint according to the two-dimensional velocity components and the two-dimensional acceleration components at both ends of the initial rod segment; The flow field velocity information of the midpoint is acquired based on the two-dimensional flow field data and the coordinate information of the midpoint.

5. The method for automatic division of wave loads based on flow field distribution according to claim 1 is characterized in that: The preset conditions include: The difference between the flow field velocity information and the linear velocity information at the midpoint is less than the absolute value of the linear velocity information multiplied by a preset ratio; Alternatively, the difference between the flow field velocity information and the linear velocity information of the midpoint is smaller than a preset unit velocity threshold.

6. The method for automatic division of wave loads based on flow field distribution according to claim 1 is characterized in that: The step of dividing the initial rod segment into two sub-rod segments according to the coordinate information of the midpoint, and dividing the sub-rod segments by midpoints according to the coordinate information of both ends of the sub-rod segments and the flow field velocity information, comprises: Dividing the initial rod segment into a first sub-rod segment from the starting end of the initial rod segment to the midpoint and a second sub-rod segment from the midpoint to the end of the initial rod segment according to the coordinate information of the midpoint; Calculate the linear velocity information and the flow field velocity information of the midpoint of the first sub-rod segment according to the coordinate information, coordinate information and flow field velocity information of the starting end of the initial rod segment and the midpoint, and when the difference between the linear velocity information and the flow field velocity information of the midpoint of the first sub-rod segment does not meet the preset condition, repeatedly divide the first sub-rod segment at the midpoint to obtain a plurality of sub-rod segments; The linear velocity information and flow field velocity information of the midpoint of the second sub-rod segment are calculated based on the coordinate information, coordinate information and flow field velocity information of the midpoint and the end of the initial rod segment, and when the difference between the linear velocity information and the flow field velocity information of the midpoint of the second sub-rod segment does not meet the preset conditions, the midpoint division of the second sub-rod segment is repeated to obtain multiple sub-rod segments.

7. The method for automatic division of wave loads based on flow field distribution according to claim 6 is characterized in that: The obtaining of a plurality of division points corresponding to the initial rod segment comprises: In the midpoint division process, a plurality of midpoints whose difference does not meet the preset condition are determined as a plurality of division points corresponding to the initial rod segment.

8. The method for automatic division of wave loads based on flow field distribution according to claim 1 is characterized in that: The step of calculating the equivalent point loads at both ends of the initial rod segment according to the coordinate information of the plurality of division points and the flow field velocity information includes: Calculate the wave loads of the plurality of divided points according to the coordinate information and flow field velocity information of the plurality of divided points by using the Morrison equation; The wave loads at the plurality of divided points are converted by integration based on the principle of virtual work to obtain equivalent point loads at both ends of the initial rod segment.

9. A computer-readable storage medium storing computer-executable instructions, characterized in that: The computer executable instructions are used to execute the method for automatic division of wave loads based on flow field distribution as described in any one of claims 1 to 8.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the method for automatically dividing wave loads based on flow field distribution as described in any one of claims 1 to 8 is implemented.

Images