Ship berthing wind action measuring and calculating method based on self-adaptive Simpson algorithm

Abstract

The invention provides a ship berthing wind effect measuring and calculating method based on an adaptive Simpson algorithm, and the method comprises the steps: building a ship contour line function and a waterplane equation, and determining a waterplane segmentation position in real time; according to the ship contour coordinate point sequence, a self-adaptive Simpson algorithm is applied to calculate a contour wind area on a ship waterline, and a ship centroid coordinate and a centroid height are determined; and converting the wind speed data into the wind force of the windward action point corresponding to the centroid height by adopting a conversion model. According to the method, the centroid height and the centroid coordinate are calculated firstly, and then the predicted wind or the actually measured wind is converted to the wind speed at the centroid height of the wind area of the ship, so that the wind power and the wind level borne by the ship are judged, and the error of the influence of an anemograph or weather forecast on the actual wind power of the ship is reduced; and taking corresponding measures based on the conversion and analysis results, so as to ensure that the operation ship can timely and accurately make corresponding adjustment in a complex sea condition environment, and ensure that the ship can safely and efficiently carry out related operation.

Description

Technical Field

[0001] The present invention relates to the technical field of wind force conversion when a ship is berthing or leaving a berth, and in particular to a method for calculating the wind force when the ship is berthing or leaving a berth based on an adaptive Simpson algorithm. Background Art

[0002] When a ship berths or leaves a berth, the management department will specifically assess the wind force acting on the ship based on the predicted wind speed and the measured wind speed, which will be used to decide whether to stop port operations and whether the ship should leave the port. Under normal circumstances, the predicted wind speed from the weather forecast is used. The predicted wind speed usually refers to the wind at a height of 10m above the sea surface, and the sea surface is relatively smooth and unobstructed. This is the wind under relatively ideal conditions, but the predicted wind is not the wind at the center of force acting on the ship. Another method uses the wind speed measured on-site on the ship. Usually, the ship's anemometer is installed on an unobstructed mast, and the anemometer is often more than 20m away from the sea surface. Wind speed is highly correlated with altitude. The higher the altitude, the greater the wind speed. Therefore, the use of the above two types of wind speeds to measure the wind at the height of the center of the ship's wind action often has a certain speed deviation. The above-mentioned existing technologies have the following limitations: the existing technologies lack a dynamic adjustment mechanism for the real-time operating status of the ship and the external sea conditions, and it is difficult to consider the dynamic changes in the state of the ship during actual operation. Because during the operation of the ship, the dynamic characteristics of the hull and changes in the external environment will cause the waterplane equation of the ship to change at any time, which easily makes it difficult to accurately measure the position of the center of the ship's wind action point.

[0003] The above-mentioned record of background technology knowledge is intended to help ordinary technicians in this field understand the existing technology that is relatively close to the present invention, and at the same time facilitate the understanding of the inventive concept and technical solution of the present invention. It should be clear that in the absence of clear evidence that the above-mentioned content has been disclosed before the filing date of this patent application, the above-mentioned background technology should not be used to evaluate the novelty of the technical solution of this application. Summary of the Invention

[0004] Technical issues Based on the consideration of the defects of the existing technology in this field in the background technology, it is necessary to convert the predicted wind speed or the measured wind speed into the wind speed at the centroid of the ship's windward area, and then use it to judge the wind force and wind level the ship is experiencing, which is used for port loading and unloading operations and unberthing safety decisions. In addition, the wind force can also be used to calculate the numerical value of the wind force acting on the ship and the ship's turning moment caused by the wind.

[0005] Technical Solution In order to achieve the above-mentioned purpose, the inventors of this application have conducted in-depth research and found that by establishing a waterplane equation, the windward side area of the ship and the corresponding centroid height, that is, the vertical distance between the centroid coordinates and the waterplane, can be calculated in real time according to the draft change, so as to ensure that the ship can still be accurately measured in a harsh port environment, avoid misoperation of the ship, thereby enhancing the operational safety of the ship, and effectively preventing the safety hazards caused by not understanding the specific wind force of the ship. It is also possible to convert the predicted wind or measured wind into the wind speed at the centroid height of the ship's windward area, and then use it to judge the size of the wind force and the corresponding wind level received by the ship, thereby improving the accuracy of the wind force measurement at the ship's windward point, thereby reducing the error of the anemometer or weather forecast on the actual wind force of the ship, and reducing the risk of operational errors. Based on the above conversion and analysis results, countermeasures are taken to ensure that in a complex sea condition, the operating ship can make corresponding adjustments in a timely and accurate manner, ensuring that the ship can carry out related operations safely and efficiently.

[0006] That is, the present invention is: (1) A method for calculating the centroid coordinates or centroid height of a ship berthing based on an adaptive Simpson algorithm, including: Establish ship contour function and waterplane equation to determine waterplane segmentation position in real time; The adaptive Simpson algorithm is applied to calculate the windward area of the ship's waterline contour through the ship's contour coordinate point sequence, and the ship's centroid coordinates or centroid height are determined based on the windward area of the ship's waterline contour.

[0007] Furthermore, the wind-exposed area of the ship's waterline profile is calculated according to the following formula:

[0008] Where A represents the wind-exposed area of the ship's outline on the waterline; and Respectively represent the minimum and maximum values of the horizontal coordinates of the polygon vertices of the ship's contour on the waterline; Represents the contour function.

[0009] Furthermore, the centroid coordinates are calculated according to the following formula:

[0010] in, is the horizontal and vertical coordinates of the centroid coordinates; and is the coordinate point sequence of the ship's waterline contour The coordinates of the point; and is the coordinate point sequence of the ship's waterline contour The coordinates of the point; is the total number of contour coordinate points; It is the wind-exposed area of the ship's contour on the waterline.

[0011] Furthermore, the centroid height is calculated according to the following formula:

[0012] in, Indicates the centroid height, that is, the vertical distance from the centroid coordinate position to the waterline plane; It is the bow and stern water gauge position of the ship and the real-time water gauge reading; are the horizontal and vertical coordinates of the ship's centroid.

[0013] (2) A method for calculating the wind effect on a ship berthing based on an adaptive Simpson algorithm includes the following steps: Establish ship contour function and waterplane equation to determine waterplane segmentation position in real time; Through the ship's contour coordinate point sequence, the adaptive Simpson algorithm is applied to calculate the wind-exposed area of the ship's contour on the waterline, and determine the coordinates and height of the ship's centroid. A conversion model is used to convert wind speed data to the wind force at the wind-receiving point corresponding to the centroid height.

[0014] Furthermore, the ship profile function is established by the following steps: Import ship engineering drawing data into the software and establish a rectangular coordinate system to generate the ship contour function .

[0015] Furthermore, the waterplane equation is established by the following steps: According to the bow water gauge coordinates and the stern water gauge coordinates ,in, and The coordinate system of the ship's front and rear water gauge positions The projected abscissa of the axis; is the ship's first draft reading, is the ship's stern draft reading; Constructing the waterplane equation , where the coefficient are the coefficients given by the general equation of the line, is a constant, Not all zero at the same time; Substitute the bow and stern water gauge coordinates into the waterplane equation , determine the coefficients of the equation.

[0016] Furthermore, the adaptive Simpson algorithm is applied to calculate the wind-exposed area of the ship's waterline contour through the ship's contour coordinate point sequence, specifically including: The waterplane surface equation With the contour function Solve the intersection points together and split to obtain the area above the waterline; Get the ship outline coordinate point sequence The minimum and maximum values of the horizontal coordinates of the polygon vertices on the ship's waterline contour are and , using the adaptive Simpson algorithm to calculate the windward area above the waterline :

[0017] Where A represents the windward area of the ship's waterline profile to be calculated; Represents the contour function.

[0018] Furthermore, the calculation formula of the centroid coordinates is:

[0019] in, is the horizontal and vertical coordinates of the centroid coordinates; and is the coordinate point sequence of the ship's waterline contour The coordinates of the point; and is the coordinate point sequence of the ship's waterline contour The coordinates of the point; is the total number of contour coordinate points; It is the wind-exposed area of the ship's contour on the waterline.

[0020] Furthermore, using the waterplane equation and centroid coordinates , calculate the vertical distance from the centroid coordinate to the waterline plane, i.e. the centroid height, by the following formula: :

[0021] in, It is the ship's bow and stern water gauge positions and real-time water gauge readings; are the horizontal and vertical coordinates of the ship's centroid.

[0022] Furthermore, the obtaining of wind speed data includes obtaining predicted wind data and / or measured wind data, wherein the predicted wind data includes: predicted wind speed With predicted wind height ; The measured wind data includes: measured wind speed Measured wind height .

[0023] Furthermore, the conversion model is used to convert the predicted wind speed and / or measured wind speed at a certain height into an equivalent wind speed at the centroid height.

[0024] Furthermore, the equivalent wind speed at the centroid height is obtained according to the following first conversion model, and the formula is as follows:

[0025] in, is the equivalent wind speed at the centroid height obtained by converting the predicted wind speed based on the first conversion model; To predict wind speed; To predict wind height; is the equivalent wind speed at the centroid height obtained by converting the measured wind speed based on the first conversion model; is the measured wind speed; is the measured wind height; is the centroid height; The length of the rough ground.

[0026] Furthermore, the ground roughness length The value is usually 0.01m (relative to the sea surface) Furthermore, the equivalent wind speed at the centroid height is obtained according to the second conversion model described below, and the formula is as follows:

[0027] in, is the equivalent wind speed at the centroid height obtained by converting the predicted wind speed based on the second conversion model; To predict wind speed; To predict wind height; is the equivalent wind speed at the centroid height obtained by converting the measured wind speed based on the second conversion model; is the measured wind speed; is the measured wind height; is the centroid height; is the conversion index.

[0028] Furthermore, the value range of the conversion index α is shown in Table 1: Table 1. Value range of conversion index α

[0029] Furthermore, the conversion index α ranges from 0.1 to 0.25.

[0030] Furthermore, the equivalent wind speed at the centroid height is obtained according to the third conversion model described below, and the formula is as follows:

[0031] in, is the equivalent wind speed at the centroid height obtained by converting the predicted wind speed based on the third conversion model; To predict wind speed; To predict wind height; is the equivalent wind speed at the centroid height obtained by converting the measured wind speed based on the third conversion model; is the measured wind speed; is the measured wind height; is the centroid height; is the ground roughness length; α is the conversion index.

[0032] Furthermore, the ground roughness length The value is usually 0.01m (relative to the sea surface) Furthermore, the conversion index α ranges from 0.1 to 0.25.

[0033] Furthermore, the equivalent wind speed at the centroid height is obtained according to the fourth conversion model described below, and the formula is as follows:

[0034] in, To predict wind friction speed; is the coefficient; is the rough length of the ground; is the measured wind friction velocity; is the equivalent wind speed at the centroid height obtained by converting the predicted wind speed based on the fourth conversion model; To predict wind speed; To predict wind height; is the equivalent wind speed at the centroid height obtained by converting the measured wind speed based on the fourth conversion model; is the measured wind speed; is the measured wind height; is the centroid height. First, and Substitute the known information into the first equation to obtain the predicted wind friction speed , then and Substitute into the third formula to obtain , similarly we can find .

[0035] Furthermore, the coefficient The value is usually 0.4.

[0036] Furthermore, the ground roughness length The value is usually taken as 0.01m (relative to the sea surface).

[0037] In the method of the present invention, a ship contour line function and a water surface line equation are first constructed to determine the waterplane segmentation position in real time. Then, an adaptive Simpson algorithm is applied to obtain a sequence of ship contour coordinate points and the minimum and maximum values of the horizontal coordinates of the polygon vertices of the ship's contour on the waterline. The adaptive Simpson algorithm formula is used to calculate the wind-exposed area of the ship's contour on the waterline, and the coordinates and centroid height of the ship are determined. Finally, based on the above-mentioned multiple conversion models including the first conversion model, the second conversion model, the third conversion model and the fourth conversion model, corresponding different equivalent wind speeds can be obtained from the predicted wind data and / or the measured wind data. , that is, the above-mentioned predicted wind and measured wind are converted into the wind speed at the centroid height of the ship's windward area, which is then used to judge the size of the wind force and the corresponding wind level received by the ship, and used for port loading and unloading operations and unberthing safety decisions, to ensure that the ship can still be accurately measured in harsh port environments, thereby avoiding misoperation of the ship, thereby enhancing the operational safety of the ship, and effectively preventing safety hazards caused by not understanding the specific wind force received by the ship.

[0038] Furthermore, it also includes: the equivalent wind speed By comparing with the preset wind speed threshold, a signal of safe ship operation status is generated.

[0039] (3) A ship berthing wind effect measurement system based on an adaptive Simpson algorithm, wherein the system executes or does not execute the method described in (1) or (2); The system comprises: Ship contour function generation module, used to establish a rectangular coordinate system, import ship engineering drawing data and generate ship contour function ; Water gauge coordinate input module, used to receive bow water gauge coordinates and the stern water gauge coordinates ; Waterline equation generation module, which generates dynamic waterline equations based on bow and stern waterline coordinates , where the coefficient , , ; Contour segmentation module, the water plane equation With ship contour function Solve the intersection point by joint equations, and , segment to obtain the area above the waterline; The adaptive Simpson algorithm calculation module obtains the ship's contour coordinate point sequence and the minimum and maximum values of the horizontal coordinates of the polygon vertices of the ship's contour on the waterline, and uses the adaptive Simpson algorithm formula to calculate the wind-exposed area of the ship's contour on the waterline ; The centroid coordinate calculation module is based on the contour coordinate points on the waterline of the ship and uses the area weighted method to calculate the centroid coordinates of the wind-exposed area. ; Centroid height calculation module, using the waterplane equation and centroid coordinates The vertical distance from the centroid coordinate to the waterline plane, i.e. the centroid height, is calculated using the point-line distance formula. ; A wind data input module, configured to receive predicted wind data and / or measured wind data; Measured wind height conversion module, used to process the straight-line distance between the anemometer position and the waterline Height conversion; The conversion module uses a conversion model to convert the wind speed data to the wind force at the wind-receiving point corresponding to the centroid height.

[0040] Furthermore, a safety decision module is included for converting the equivalent wind speed By comparing with the preset wind speed threshold, a signal of safe ship operation status is generated.

[0041] (4) A computer device comprising a memory, a processor, a communication interface, and a communication bus; wherein the memory, the processor, and the communication interface communicate with each other via the communication bus; the memory is used to store a computer program; the processor is used to execute the computer program stored in the memory, and when the processor executes the computer program, it implements at least one step of the method described in (1) or (2) above.

[0042] (5) A computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements at least one step of the method described in (1) or (2) above.

[0043] On the basis of conforming to the common sense in this field, the above-mentioned preferred conditions can be combined with each other to obtain a specific implementation method.

[0044] Beneficial effects This method for measuring the wind effects on ships at berthing, based on the adaptive Simpson algorithm, establishes a waterplane equation and can calculate the windward side area of the ship and the corresponding centroid height in real time according to the draft changes, ensuring that the ship can still be accurately measured in severe sea conditions, thereby avoiding ship misoperation, thereby enhancing the operational safety of the ship and effectively preventing safety hazards caused by not understanding the specific wind force acting on the ship.

[0045] This method for calculating wind speed effects on ships during berthing, based on the adaptive Simpson algorithm, combines the adaptive Simpson algorithm with a conversion model to provide a comprehensive set of wind speed calculation methods. By comparing measured wind speeds with predicted wind data, the accuracy of wind speed calculations at the ship's wind-exposed point is improved, thereby reducing the error in the influence of anemometers or weather forecasts on the ship's actual wind speeds. This reduces the risk of operational errors and enables timely and accurate adjustments to ship berthing and unberthing operations in complex sea conditions, ensuring safe and efficient ship operations.

[0046] By converting measured and predicted winds to wind forces at the ship's center of gravity, this present invention addresses the ambiguity inherent in conventional methods, which ignore the variations in wind forces with altitude. Conventional methods typically assume that wind forces act uniformly on the ship, ignoring the variations in wind forces with altitude, resulting in significant deviations between calculated results and actual forces. The present invention employs the aforementioned technical solution to achieve this objective, addressing the shortcomings of existing technologies with a rational design and ease of operation. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to make the above-mentioned and / or other purposes, features, advantages and examples of the present invention more obvious and easy to understand, the following is a brief introduction to the drawings required for use in the specific embodiments of the present invention. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative work.

[0048] Figure 1 A simplified diagram showing the technical logic of this application; Figure 2 It represents the centroid coordinates and centroid high-speed calculation flow chart in the technical route of this application; Figure 3 A flowchart showing how to convert measured wind and predicted wind into equivalent wind speed at the wind action point in the technical route of this application; Figure 4 A simplified diagram showing the ship No. II; Figure 5 A schematic diagram showing the conversion of predicted wind into equivalent wind force by ship I based on the first conversion model; Figure 6 Schematic diagram showing how ship I converts measured wind into equivalent wind force based on the first conversion model; Figure 7 A schematic diagram showing the conversion of predicted wind into equivalent wind force based on the first conversion model for ship II; Figure 8 Schematic diagram showing how ship II converts measured wind into equivalent wind force based on the first conversion model; Figure 9 A schematic diagram showing the conversion of predicted wind into equivalent wind force by ship I based on the second conversion model; Figure 10 Schematic diagram showing how ship I converts measured wind into equivalent wind force based on the second conversion model; Figure 11 A schematic diagram showing how ship II converts the predicted wind into equivalent wind force based on the second conversion model; Figure 12 Schematic diagram showing how ship II converts measured wind into equivalent wind force based on the second conversion model; Figure 13 A schematic diagram showing the conversion of predicted wind into equivalent wind force by ship I based on the third conversion model; Figure 14 Schematic diagram showing how ship I converts measured wind into equivalent wind force based on the third conversion model; Figure 15 Schematic diagram showing the conversion of predicted wind into equivalent wind force based on the third conversion model for ship II; Figure 16 Schematic diagram showing how the measured wind is converted into equivalent wind force based on the third conversion model for ship II; Figure 17 A schematic diagram showing the conversion of predicted wind into equivalent wind force by ship I based on the fourth conversion model; Figure 18 Schematic diagram showing how ship I converts measured wind into equivalent wind force based on the fourth conversion model; Figure 19 Schematic diagram showing the conversion of predicted wind into equivalent wind force based on the fourth conversion model for ship II; Figure 20 Schematic diagram showing how Ship II converts measured wind into equivalent wind force based on the fourth conversion model. DETAILED DESCRIPTION

[0049] Those skilled in the art may refer to the contents herein and appropriately substitute and / or modify the process parameters to achieve the desired effect. However, it should be noted that all such substitutions and / or modifications are obvious to those skilled in the art and are considered to be included in the present invention. The products and preparation methods described herein have been described through preferred embodiments. It is obvious that those skilled in the art can modify or appropriately change and combine the products and preparation methods described herein to implement and apply the technology of the present invention without departing from the content, spirit, and scope of the present invention.

[0050] Unless otherwise defined, technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The present invention utilizes the methods and materials described herein; however, other suitable methods and materials known in the art may also be used. The materials, methods, and examples described herein are illustrative only and are not intended to be limiting. All publications, patent applications, patents, provisional applications, database entries, and other references mentioned herein are incorporated herein by reference in their entirety. In the event of a conflict, the present specification, including definitions, will control.

[0051] Unless otherwise specified, the materials, methods, and examples described herein are illustrative only and not limiting. Although methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present invention, suitable methods and materials are described herein.

[0052] The technical solutions in the embodiments of the present invention will be described clearly and completely below. Obviously, the embodiments described are only part of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention. At the same time, the embodiments in this application and the features in the embodiments can be combined with each other unless there is a conflict.

[0053] It should be understood that any technical solution claimed for protection in the present invention does not involve the diagnosis and treatment of diseases.

[0054] To facilitate understanding of the embodiments of the present invention, abbreviations and key terms that may be involved in the embodiments of the present invention are first explained or defined. Undefined abbreviations or key terms are generally understood by those skilled in the art.

[0055] In addition, the experimental methods used in the examples are conventional methods unless otherwise specified. The materials and reagents used are commercially available unless otherwise specified. Reagents or instruments used without manufacturer indication are commercially available conventional products. All publications and other references cited herein are incorporated by reference in their entirety.

[0056] The present invention is described in detail below Example 1: Provides a method for calculating the wind effect on ships berthing based on the adaptive Simpson algorithm. The technical route flow chart is as follows: Figure 1-Figure 3As shown in the figure, the ship's contour function and waterplane equation are first established to determine the waterplane segmentation location in real time. Then, the adaptive Simpson algorithm is applied to obtain the ship's contour coordinate point sequence and the minimum and maximum horizontal coordinates of the polygon vertices of the ship's waterline contour. The adaptive Simpson algorithm formula is used to calculate the wind-exposed area of the ship's contour on the waterline and determine the coordinates and height of the ship's centroid. Finally, a conversion model is used to convert the wind speed data to the wind force at the wind-exposed point corresponding to the centroid height. The specific steps are as follows.

[0057] Establish a rectangular coordinate system in CAD, import ship CAD drawing data, and draw the ship contour function ; According to the front and rear water gauge position data, establish the bow water gauge coordinates and the stern water gauge coordinates ,in, and The coordinate system of the ship's front and rear water gauge positions The projected abscissa of the axis; is the ship's first draft reading, It is the ship's stern draft reading.

[0058] Establishing the ship waterplane equation , record the front and rear draft through the water gauge function and , substitute the water gauge coordinates into , determine the equation coefficients, split the ship outline into upper and lower parts, The equation is shown as follows:

[0059] Where: are the coefficients given by the general equation of the line; is a constant, The above coefficients can be calculated using the coordinates of two points: , , .

[0060] The waterplane surface equation With the contour function The intersection point is obtained by the combined solution and , dividing the waterline into two areas above and below.

[0061] Using the adaptive Simpson algorithm and ship contour function Calculate the wind-exposed area of the ship's waterline contour and obtain the coordinate point sequence of the ship's contour within the regular shape area of the ship's engineering drawings. The minimum and maximum values of the horizontal coordinates of the polygon vertices on the ship's waterline contour are and , using the adaptive Simpson algorithm to calculate the windward area above the waterline :

[0062] Where A represents the windward area of the ship's waterline profile to be calculated; Represents the contour function.

[0063] It should be explained that the waterplane equation can calculate the windward side area of the ship and the corresponding centroid height and height from the waterplane in real time according to the draft changes.

[0064] Calculating the centroid coordinates of the ship's waterline using the weighted method , as shown in the formula:

[0065] in, and is the coordinate point sequence of the ship's waterline contour The coordinates of the point; and is the coordinate point sequence of the ship's waterline contour The coordinates of the point; is the total number of contour coordinate points; It is the wind-exposed area of the ship's contour on the waterline.

[0066] The centroid coordinates on the ship are the height of the center of the wind action point of the ship from the waterline surface. The waterline surface equation W and the centroid coordinates of the part on the waterline are used to calculate the centroid coordinates. , calculate the vertical distance from the centroid to the waterline, recorded as the centroid height , as shown in the following formula:

[0067] in, It is the ship's bow and stern water gauge positions and real-time water gauge readings; The obtained centroid height is an important parameter in the subsequent logarithmic model.

[0068] Wind speed data consists of the following two parts: predicted wind data (predicted wind speed With predicted wind height ), measured wind data (measured wind speed Measured wind height ). The predicted wind speed represents the wind force in the weather forecast; represents the actual wind speed measured by the ship's anemometer; the corresponding wind heights are and : is the wind reference height for weather forecast, in this embodiment, =10m; is the vertical distance from the anemometer height to the waterline; using the waterplane equation and anemometer location coordinates Calculate the vertical distance between the anemometer and the waterline as shown in the formula:

[0069] like Figure 1 As shown, The anemometer position coordinates (like Figure 1 ④) to the waterline surface equation W (such as Figure 1 The vertical distance of ③) (such as Figure 1 ⑥), where the waterline equation W can be recorded by the water gauge function through it. and The coordinate point is represented, so the vertical distance between the anemometer height and the waterline can be calculated by the above formula .

[0070] If the wind forecast under the weather forecast is used, then 、 Height to centroid , substituted into the first conversion model formula, and the equivalent wind speed at the centroid height obtained by converting the predicted wind speed based on the first conversion model is obtained. If the actual wind speed measured by the ship's own anemometer is used, then the 、 Height to centroid , substituted into the first conversion model formula, and the centroid height obtained by converting the measured wind speed based on the first conversion model is obtained. The first conversion model formula is as follows:

[0071] in, is the ground roughness length, which is 0.01m.

[0072] This example uses the actual ship drawings of a 208,947-ton bulk carrier (denoted as Ship I) to calculate wind effects on the ship. The ship has a length of 299.95 m, a length between perpendiculars of 294.75 m, a molded breadth of 50 m, and a molded depth of 25.20 m. The ship's profile dataset is shown in Table 2-4.

[0073] Table 2. The first part of the No. 1 ship profile dataset

[0074] Table 3. The second part of the No. 1 ship profile dataset

[0075] Table 4. Part 3 of the Ship Profile Dataset No. 1

[0076] Furthermore, this embodiment uses the actual ship drawings of a bulk carrier with a gross tonnage of 1314 tons and a displacement of 1768.4 tons (denoted as Ship II) as the basis for calculating the wind effects on the ship. The ship has an overall length of 63.50 meters, an overall length excluding bulwarks of 61.97 meters, a molded width of 12.60 meters, a molded depth of 4.60 meters, and a designed draft of 3.20 meters. The simplified ship drawings of Ship II are shown in the following figure: Figure 4 shown.

[0077] For ship I, the first conversion model is applied to convert the predicted wind into the equivalent wind speed at the centroid height (9.6m). Figure 5 shown.

[0078] For ship I, the first conversion model is applied to convert the measured wind into the equivalent wind speed at the centroid height (9.6m). Figure 6 shown.

[0079] For ship II, the first conversion model is applied to convert the predicted wind into the equivalent wind speed at the centroid height (2.74312m). Figure 7 shown.

[0080] For ship II, the first conversion model is applied to convert the measured wind into the equivalent wind speed at the centroid height (2.74312m). Figure 8 shown.

[0081] Example 2: Based on the above embodiment, the second conversion model is used to replace the first conversion model to obtain the equivalent wind speed at the centroid height. The formula is as follows:

[0082] in, is the equivalent wind speed at the centroid height obtained by converting the predicted wind speed based on the second conversion model; To predict wind speed; To predict wind height; is the equivalent wind speed at the centroid height obtained by converting the measured wind speed based on the second conversion model; is the measured wind speed; is the measured wind height; is the centroid height; is the conversion index, with a value of 0.1.

[0083] For ship I, the second conversion model is applied to convert the predicted wind into the equivalent wind speed at the centroid height. Figure 9 shown.

[0084] For ship I, the second conversion model is applied to convert the measured wind into the equivalent wind speed at the centroid height. Figure 10 shown.

[0085] For ship II, the second conversion model is applied to convert the predicted wind into the equivalent wind speed at the centroid height. Figure 11 shown.

[0086] For ship II, the second conversion model is applied to convert the measured wind into the equivalent wind speed at the centroid height. Figure 12 shown.

[0087] Example 3: Based on the above embodiment 1, the equivalent wind speed at the centroid height is obtained using the third conversion model, and the formula is as follows:

[0088] in, is the equivalent wind speed at the centroid height obtained by converting the predicted wind speed based on the third conversion model; To predict wind speed; To predict wind height; is the equivalent wind speed at the centroid height obtained by converting the measured wind speed based on the third conversion model; is the measured wind speed; is the measured wind height; is the centroid height; is the ground roughness length; α is the conversion index.

[0089] For ship I, the third conversion model is applied to convert the predicted wind into the equivalent wind speed at the centroid height. Figure 13 shown.

[0090] For ship I, the third conversion model is applied to convert the measured wind into the equivalent wind speed at the centroid height. Figure 14 shown.

[0091] For ship II, the third conversion model is applied to convert the predicted wind into the equivalent wind speed at the centroid height. Figure 15 shown.

[0092] For ship II, the third conversion model is applied to convert the measured wind into the equivalent wind speed at the centroid height. Figure 16 shown.

[0093] Example 4: Based on the above embodiment 1, the equivalent wind speed at the centroid height is obtained using the fourth conversion model, and the formula is as follows:

[0094] in, To predict wind friction speed; is the coefficient, take 0.4; is the ground roughness length, which is 0.01m; is the measured wind friction velocity; is the equivalent wind speed at the centroid height obtained by converting the predicted wind speed based on the fourth conversion model; To predict wind speed; To predict wind height; is the equivalent wind speed at the centroid height obtained by converting the measured wind speed based on the fourth conversion model; is the measured wind speed; is the measured wind height; is the centroid height. First, and Substitute the known information into the first equation to obtain the predicted wind friction speed , then and Substitute into the third formula to obtain , similarly we can find .

[0095] For ship I, the fourth conversion model is applied to convert the predicted wind into the equivalent wind speed at the centroid height. Figure 17 shown.

[0096] For ship I, the fourth conversion model is applied to convert the measured wind into the equivalent wind speed at the centroid height. Figure 18 shown.

[0097] For ship II, the fourth conversion model is applied to convert the predicted wind into the equivalent wind speed at the centroid height. Figure 19 shown.

[0098] For ship II, the fourth conversion model is applied to convert the measured wind into the equivalent wind speed at the centroid height. Figure 20 shown.

[0099] According to the technical solution of the present invention, the contour area, centroid coordinates and centroid height of the ship on the waterline are first obtained based on the adaptive Simpson algorithm. The centroid height of ship No. I is calculated to be 9.6m, which is relatively close to the predicted wind height of 10m and the measured wind height of 11.15121m; while the centroid height of ship No. II is calculated to be 2.74312m, which is far from the predicted wind height of 10m and the measured wind height of 11.15121m. It can be seen from ships No. I and No. II that for different ship types, the present invention can accurately obtain the centroid coordinates and centroid height by applying the adaptive Simpson algorithm, which is beneficial to determine the actual wind action point of the ship, thereby facilitating the accurate and real-time adjustment of ship operating parameters and improving ship operation safety.

[0100] from Figure 5-Figure 20 It can be seen that based on the first to fourth conversion models of the scheme of this embodiment, the predicted wind under the weather forecast and the wind speed measured by the ship's own anemometer are converted into the wind speed at the point where the ship is affected by the wind, that is, the centroid coordinates, and the equivalent wind speed at the centroid height is compared and judged with the wind force requirements of the terminal operating environment to achieve the expected effect.

[0101] This implementation uses normalization and comprehensive analysis of multiple key parameters of the vessel's berthing environment (such as the vessel's profile coordinates, draft, anemometer height, weather forecast, and anemometer data) to accurately calculate the wind force at the center of the ship's wind-exposed point. This avoids errors in different parameters due to environmental variations, thereby reducing the impact of each parameter on the actual wind conditions experienced by the ship during berthing and unberthing. Simultaneously, under the influence of multiple factors, the vessel's draft fully reflects the vessel's actual needs under different operating conditions, thereby avoiding errors caused by real-time changes in the vessel's waterplane equation. This allows the vessel to navigate stably under various navigation conditions, ensuring safety and efficiency. Furthermore, through the adaptive Simpson algorithm and conversion model, comprehensive consideration is taken into account for height conversion, thereby accurately optimizing the vessel's adjustment efficiency and stability in various complex environments during actual operational scenarios, thereby enhancing the vessel's operational safety.

[0102] The conventional techniques in the above embodiments are prior arts known to those skilled in the art, and thus will not be described in detail here.

[0103] The specific embodiments described herein are merely illustrative of the spirit of the present invention. Persons skilled in the art may make various modifications, additions, or substitutions to the described specific embodiments without departing from the spirit of the present invention or exceeding the scope defined by the appended claims.

[0104] Although the present invention has been described in detail and certain specific embodiments have been cited, it will be apparent to those skilled in the art that various changes and modifications can be made without departing from the spirit and scope of the invention.

[0105] Although the above-mentioned specific embodiments have shown, described and pointed out the novel features applied to various embodiments, it should be understood that various omissions, replacements and changes can be made to the form and details of the described devices or methods without departing from the spirit of the present disclosure. In addition, the various features and methods described above can be used independently of each other, or can be combined in various ways. All possible combinations and sub-combinations are intended to fall within the scope of the present disclosure. Many of the above-mentioned embodiments include similar components, and therefore, these similar components are interchangeable in different embodiments. Although the present invention has been disclosed in the context of certain embodiments and examples, it should be understood by those skilled in the art that the present invention can extend beyond the specifically disclosed embodiments to other alternative embodiments and / or applications and their obvious modifications and equivalents. Therefore, the present invention is not intended to be limited by the specific disclosure of the preferred embodiments herein.

[0106] Matters not covered in the present invention are all known technologies.

Claims

1. A method for calculating the centroid coordinates of a ship berthing based on an adaptive Simpson algorithm, characterized in that: include: Establish ship contour function and waterplane equation to determine waterplane segmentation position in real time; The adaptive Simpson algorithm is applied to the ship contour coordinate point sequence to calculate the windward area of the ship's contour on the waterline, and the coordinates of the ship's centroid are determined according to the windward area of the ship's contour on the waterline.

2. The method according to claim 1, wherein: The windward area of the ship's waterline contour is calculated using the following formula: Where A represents the wind-exposed area of the ship's outline on the waterline; and Respectively represent the minimum and maximum values of the horizontal coordinates of the polygon vertices of the ship's contour on the waterline; Represents the contour function.

3. The method according to claim 2, wherein: The centroid coordinates are calculated according to the following formula: in, is the horizontal and vertical coordinates of the centroid coordinates; and is the coordinate point sequence of the ship's waterline contour The coordinates of the point; and is the coordinate point sequence of the ship's waterline contour The coordinates of the point; is the total number of contour coordinate points; It is the wind-exposed area of the ship's contour on the waterline.

4. The method according to claim 2 or 3, characterized in that: The centroid height is calculated according to the following formula: in, Indicates the centroid height, that is, the vertical distance from the centroid coordinate position to the waterline plane; It is the bow and stern water gauge position of the ship and the real-time water gauge reading; are the horizontal and vertical coordinates of the ship's centroid.

5. A method for calculating wind effects on ships berthing based on an adaptive Simpson algorithm is characterized by: include: Establish ship contour function and waterplane equation to determine waterplane segmentation position in real time; Through the ship's contour coordinate point sequence, the adaptive Simpson algorithm is applied to calculate the wind-exposed area of the ship's contour on the waterline, and determine the coordinates and height of the ship's centroid. The conversion model is used to convert the wind speed data to the equivalent wind speed at the wind-affected point corresponding to the centroid height.

6. The method according to claim 5, characterized in that The method of calculating the wind-exposed area of the ship's waterline contour by applying the adaptive Simpson algorithm through the ship's contour coordinate point sequence specifically includes: The waterplane surface equation With the contour function Solve the intersection points together and split to obtain the area above the waterline; Get the ship outline coordinate point sequence The minimum and maximum values of the horizontal coordinates of the polygon vertices on the ship's waterline contour are and , using the adaptive Simpson algorithm to calculate the windward area above the waterline : Where A represents the windward area of the ship's waterline profile to be calculated; Represents the contour function.

7. The method according to claim 5, characterized in that The conversion model is used to convert the predicted wind speed or the measured wind speed into an equivalent wind speed at the centroid height.

8. The method according to any one of claims 5 to 7, characterized in that: The conversion model converts the equivalent wind speed at the centroid height according to the first conversion model shown in the following formula: in, is the equivalent wind speed at the centroid height obtained by converting the predicted wind speed based on the first conversion model; To predict wind speed; To predict wind height; is the equivalent wind speed at the centroid height obtained by converting the measured wind speed based on the first conversion model; is the measured wind speed; is the measured wind height; is the centroid height; The length of the rough ground.

9. The ship berthing wind effect measurement system based on the adaptive Simpson algorithm is characterized by: include: Software module for establishing a rectangular coordinate system, importing ship engineering drawing data and generating the ship contour function f(x); Water gauge coordinate input module, used to receive the water gauge coordinates of the bow of the ship and the stern water gauge coordinates ; Waterline equation generation module, which generates dynamic waterline equations based on bow and stern waterline coordinates ; Contour segmentation module, the water plane equation With ship contour function Solve the intersection points together to split the area above and below the waterline; The adaptive Simpson algorithm calculation module obtains the ship's contour coordinate point sequence and the minimum and maximum values of the horizontal coordinates of the polygon vertices of the ship's contour on the waterline, and uses the adaptive Simpson algorithm formula to calculate the wind-exposed area of the ship's contour on the waterline ; The centroid coordinate calculation module calculates the centroid coordinates of the wind-exposed area based on the partial contour coordinate points on the ship and using the area weighted method. ; Centroid height calculation module, using the waterplane equation and centroid coordinates The vertical distance from the centroid to the waterline plane, i.e. the centroid height, is calculated using the point-line distance formula. ; Measured wind height conversion module, used to process the straight-line distance between the anemometer position and the waterline Height conversion; A wind data acquisition module, configured to receive predicted wind data and / or measured wind data; The equivalent wind speed conversion module is used to convert the wind force data into the wind speed at the wind-receiving point corresponding to the centroid height, that is, the equivalent wind speed, in combination with the centroid height.

10. Use of the method according to any one of claims 5 to 8 or the system according to claim 9 in a ship operation process, characterized in that: The equivalent wind speed at the centroid height is compared with the preset wind threshold to generate a ship operation safety status signal, which is used for ship safety. If the wind speed exceeds the threshold, an adjustment command is generated to control the ship's power system or steering gear to dynamically correct the berthing posture.

Images