CN113673092B - Calculation method for window period proportion of pile driving ship capable of operating - Google Patents
Calculation method for window period proportion of pile driving ship capable of operating Download PDFInfo
- Publication number
- CN113673092B CN113673092B CN202110881061.3A CN202110881061A CN113673092B CN 113673092 B CN113673092 B CN 113673092B CN 202110881061 A CN202110881061 A CN 202110881061A CN 113673092 B CN113673092 B CN 113673092B
- Authority
- CN
- China
- Prior art keywords
- ship
- piling
- piling ship
- freedom
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Placing Or Removing Of Piles Or Sheet Piles, Or Accessories Thereof (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
The invention discloses a calculation method of window period proportion of a pile driving ship capable of working, which comprises the following steps: step one, collecting wave data of a piling ship every 3 hours during construction of an open sea area; calculating amplitude response operators of each degree of freedom of the piling ship by using hydrodynamic analysis software; step three, calculating the motion response energy density spectrum of each degree of freedom of the piling ship by adopting a frequency domain calculation method; step four, calculating the amplitude corresponding to the harmonic wave of the motion response of each degree of freedom of the piling ship according to the motion response energy spectrum of each degree of freedom of the piling ship; calculating motion responses of all degrees of freedom of the piling ship by adopting a time domain calculation method; and step six, calculating the operable window period of the piling ship according to the operable threshold balance of the piling ship. The invention adopts a time-frequency domain joint calculation method, reasonably utilizes the advantages of the two methods, calculates the floating body dynamics of the piling ship, and judges the operable window of the piling ship.
Description
Technical Field
The invention belongs to the technical field of theoretical calculation of operable conditions of engineering ships, and particularly relates to a calculation method of operable window period proportion of a piling ship.
Background
In the offshore wind power construction engineering of China, a pile driving ship occupies a relatively important position. Since the beginning of construction of the wind farm of the sea in China, namely the east China sea bridge wind farm, a large number of offshore fans in the form of high pile cap foundations with piling ships as construction equipment appear. The high pile cap foundation is a foundation type widely applied to the fields of port and channel engineering, a foundation pile group is generally formed by 6-8 steel pipe inclined piles with the diameter of 2-3 m, and a high-performance concrete poured bearing cap is arranged at the top of the foundation pile group to connect the pile groups, so that the construction of the high pile cap generally needs to adopt a professional pile driving ship and a concrete stirring ship to respectively finish pile sinking and concrete pouring construction of the pile group. At present, pile driving equipment of professional pile driving vessels in China is arranged on a fixed pile frame of a ship bow, so that if construction of 6-8 inclined piles uniformly distributed in 360 degrees is completed, the ship bow of the pile driving vessel needs to be laid out along a plurality of angles, and the bow direction of the ship cannot be selected freely.
The current offshore wind farms in China are concentrated in a water depth range of 30m to 50m, the offshore distance is up to 75km, the environmental conditions in deep water and rough sea areas are worse than those in offshore areas, the operable period of the piling ship is difficult to judge based on experience, and the operable window period of the piling ship needs to be calculated by adopting a quantitative analysis method.
In the world, because the professional piling ship is mainly used for offshore area engineering, the stormy waves are small, and the operable window period is mainly judged based on experience, the study on the motion response analysis of the piling ship at home and abroad is relatively less. But the pile driving ship is used as a floating engineering ship, and the working research under the complex sea condition can be used for referencing some research methods of other floating engineering ships. Based on the construction process of piling operation, whether the piling ship can operate or not can be judged, wherein the operation is mainly limited by the motion amplitude of three vertical degrees of freedom of ship heave, roll and pitch. The overall idea of the analysis method in this respect is: firstly, judging a risk event which possibly causes accidents or construction failure, judging critical limiting parameters of the risk event, then adopting a floating body dynamics simulation method to carry out numerical simulation on offshore construction under various sea conditions, and judging the sea conditions which can be operated and not operated by combining simulation results and the critical limiting parameters.
Numerical simulation of offshore construction is mainly performed by calculation of floating body dynamics. In the calculation, the construction ship is assumed to be a rigid body, and the six-degree-of-freedom motion amplitude of any position of the construction ship can be obtained through six-degree-of-freedom motion of the ship. Through decades of development, the current floating body dynamics simulation analysis is mature, especially the simulation of construction period motions with smaller relative motion amplitude, and the calculation results based on commercial software such as AQWA, SESAM and the like have larger reliability and reference value.
The key steps in the calculation of the operable window period of the piling ship are floating body hydrodynamics calculation, and the main method of the current floating body hydrodynamics simulation analysis comprises time domain calculation and frequency domain calculation, which have advantages and disadvantages respectively: the frequency domain calculation is more demanding, while the time domain calculation is extremely time consuming.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provide a calculation method for the operable window period proportion of a piling ship, which adopts a time-frequency domain combined calculation method, reasonably utilizes the advantages of a time domain calculation method and a frequency domain calculation method, calculates the floating body dynamics of the piling ship, and judges the operable window period of the piling ship.
The purpose of the invention is realized in the following way: a calculation method of a window period proportion of a pile driving ship capable of working comprises the following steps:
step one, collecting wave data S of a piling ship for every 3 hours in construction of an open sea area W (ω) the wave data is expressed as a function of the wave frequency ω;
calculating an amplitude response operator RAO (omega) of each degree of freedom of the piling ship by using hydrodynamic analysis software, wherein the amplitude response operator is expressed as a function of wave frequency omega;
step three, calculating the motion response energy density spectrum of each degree of freedom of the piling ship by adopting a frequency domain calculation method, namely calculating by adopting the following formula (1):
S X (ω)＝RAO 2 (ω)·S W (ω) (1)
in the formula (1), S X (omega) is the motion response energy density spectrum of a degree of freedom of the piling ship, S W (ω) is the data obtained in step one, and RAO (ω) is the data obtained in step two;
step four, calculating an amplitude A corresponding to a harmonic wave of the motion response of each degree of freedom of the piling ship according to the motion response energy spectrum of each degree of freedom of the piling ship obtained in the step three k I.e. S is firstly X (omega) is divided into n equal parts corresponding to the wave frequency range delta omega in the frequency domain, and then calculated by the following formula (2)Amplitude A of wave k ：
In the formula (2), S X (ω k ) For the kth part frequency omega k Corresponding motion response energy density spectrum, A k The amplitude corresponding to the harmonic of the kth motion response;
calculating motion responses of all degrees of freedom of the piling ship by adopting a time domain calculation method, namely calculating by adopting the following formula (3):
in the formula (3), X (t) is the motion response of the piling ship in a certain degree of freedom, t is time, and the motion response X (t) of the piling ship in a certain degree of freedom is a function of time t, epsilon k To obey [0,2 pi ]]The kth element in the uniformly distributed array; substituting different time t values to obtain the maximum value of the freedom degree motion response;
step six, calculating the operable window period of the piling ship according to the operable threshold balance of the piling ship, and performing the following steps:
step six-1, repeating the step three to the step five, and calculating the motion response extremum of each degree of freedom of the piling ship under the operation threshold balance;
step six-2, carrying out normalization processing on the motion response maximum value of each degree of freedom of the piling ship obtained by calculation in the step five, namely dividing the motion response maximum value of each degree of freedom of the piling ship calculated in the step five by the motion response extremum of each degree of freedom of the piling ship calculated in the step six-1 in a one-to-one correspondence manner;
step six-3, calculating the operable window period of the piling ship, namely judging whether the piling ship is operable or not by using the added numerical values of the results calculated in the step six-2, and if the sum is less than or equal to 1, indicating that the condition is the operable construction window period of the piling ship; if the sum is greater than 1, the condition is not the window period of the workable construction of the piling ship;
and step six-4, taking the ratio of the number of 3 hours for which the piling ship can operate to the total number of 3 hours as the ratio of the operable window period of the piling ship.
According to the method for calculating the proportion of the window period of the pile driving ship capable of working, the step one is carried out, the wave data are historical wave post-report data of 3 hours by 3 hours, the wave data are obtained through a local weather bureau or a marine bureau, and the wave data consist of stormy waves, main swells and secondary swells.
According to the calculation method of the window period proportion of the pile driving ship, in the second step, the hydrodynamic force analysis software is marine engineering hydrodynamic force analysis software based on the winding radiation principle.
The method for calculating the proportion of the operable window period of the piling ship comprises the steps of at least rolling the piling ship and swaying the pile gripper of the bow of the piling ship.
The calculation method of the window period proportion of the pile driving ship has the following effects:
1. according to the practical operation condition and the environmental condition characteristics of the piling ship, the invention relatively accurately analyzes the limit environmental condition of the piling ship, and on the basis of the limit environmental condition, the condition of whether the piling ship can operate in a certain period in the future can be predicted, so that the waiting time of the piling ship is effectively reduced, and the construction efficiency of the piling ship is improved;
2. unlike other environmental condition window calculations, the calculation model constructed in the invention is simple and convenient, and because waves appear in the form of a multi-peak spectrum mixed wave under most conditions for open sea, the energy density spectrum S of the motion response is not the same X (omega) wave spectrum S of irregular wave W (omega) it is in most cases difficult to meet the assumption requirements of Gaussian signals, so that the motion response value of the piling ship cannot be calculated directly by using the full frequency domain method, and at this time, the calculated energy density spectrum S of the motion response can be obtained X (omega) into a time domain sequence X (t) and counting the motion responses corresponding to different override probabilities in the time domainThe method is a time-frequency domain joint calculation method provided by the invention. The method has the main advantages that the floating body dynamics equation is prevented from being directly solved in the time domain, the calculation time is greatly reduced, and meanwhile, the harsher assumption condition relied on the full frequency domain calculation is avoided.
Drawings
FIG. 1 is a flow chart of a method of calculating the proportion of window periods during which a pile driving vessel of the present invention can operate;
FIG. 2 is a graph of wave data collected during the performance of the steps of the present invention as a function of wave frequency;
FIG. 3a is a graph of roll amplitude response operator of a pile driving vessel as a function of wave frequency, obtained when step two of the present invention is performed;
FIG. 3b is a graph of heave amplitude response operator of a pile driving vessel as a function of wave frequency, obtained when step two of the present invention is performed;
FIG. 4a is a plot of roll motion response energy density spectrum as a function of wave frequency for a pile driving vessel when step three of the present invention is performed;
FIG. 4b is a plot of heave motion response energy density spectrum of a pile driving vessel as a function of wave frequency, obtained when step three of the present invention is performed;
FIG. 5 is a graph comparing the maximum 3 hour roll motion response of the pile driving vessel as per step six-1 of the present invention with the extreme roll motion response at the operational threshold level;
FIG. 6 is a graph comparing the maximum value of heave motion response of the pile driving vessel per 3 hours and the extreme value of heave motion response under the operable threshold balance, obtained when carrying out step six-1 of the present invention;
FIG. 7 is a graph showing normalized results of maximum values of 3-hour roll motion response of the piling ship obtained when step six-2 of the present invention is performed;
FIG. 8 is a graph showing normalized results of maximum value of heave motion response of a piling ship for 3 hours, which is obtained when step six-2 of the invention is performed;
FIG. 9 is a schematic representation of the operational window period of a pile driving vessel according to step six-3 of the present invention;
fig. 10 is a schematic view showing the proportion of the operable window period of the pile driving vessel obtained by carrying out the step six-4 of the present invention.
Detailed Description
The invention will be further described with reference to the accompanying drawings.
Referring to fig. 1, the method for calculating the operable window period ratio of the piling ship of the present invention comprises the following steps:
step one, collecting wave data S of a piling ship for every 3 hours in construction of an open sea area W (ω) representing wave data as a function of wave frequency ω, the wave data being historical wave post-report data from 3 hours by 3 hours, the wave data being obtained by a local weather or ocean, the wave data being comprised of storms, major swells and minor swells;
calculating an amplitude response operator RAO (omega) of each degree of freedom of the piling ship by using hydrodynamic analysis software, wherein the amplitude response operator is expressed as a function of wave frequency omega; the hydrodynamic force analysis software is marine engineering hydrodynamic force analysis software based on the principle of winding radiation, such as AQWA and WAMIT (analysis software for calculating interaction between a zero-navigational floating structure and waves by Wave Analysis MIT) and the like; the degree of freedom of the piling ship is at least the roll of the piling ship and the heave of the pile gripper position of the ship bow of the piling ship;
step three, calculating the motion response energy density spectrum of each degree of freedom of the piling ship by adopting a frequency domain calculation method, namely calculating by adopting the following formula (1):
S X (ω)＝RAO 2 (ω)·S W (ω) (1)
in the formula (1), S X (omega) is the motion response energy density spectrum of a degree of freedom of the piling ship, S W (ω) is the data obtained in step one, and RAO (ω) is the data obtained in step two;
step four, calculating an amplitude A corresponding to a harmonic wave of the motion response of each degree of freedom of the piling ship according to the motion response energy spectrum of each degree of freedom of the piling ship obtained in the step three k I.e. S is firstly X (omega) dividing the wave frequency into n equal parts corresponding to the wave frequency range Deltaomega in the frequency domain, and calculating the amplitude A by the following formula (2) k ：
In the formula (2), S X (ω k ) For the kth part frequency omega k Corresponding motion response energy density spectrum, A k The amplitude corresponding to the harmonic of the kth motion response;
calculating motion responses of all degrees of freedom of the piling ship by adopting a time domain calculation method, namely calculating by adopting the following formula (3):
in the formula (3), X (t) is the motion response of the piling ship in a certain degree of freedom, t is time, and the motion response X (t) of the piling ship in a certain degree of freedom is a function of time t, epsilon k To obey [0,2 pi ]]The k-th element in the uniformly distributed array has the important function of guaranteeing the randomness of X (t), thus epsilon k Also known as random seeds; substituting different time t values to obtain the maximum value of the freedom degree motion response;
step six, calculating the operable window period of the piling ship according to the operable threshold balance of the piling ship, and performing the following steps:
step six-1, repeating the step three to the step five, and calculating the motion response extremum of each degree of freedom of the piling ship under the operation threshold balance;
step six-2, carrying out normalization processing on the motion response maximum value of each degree of freedom of the piling ship obtained by calculation in the step five, namely dividing the motion response maximum value of each degree of freedom of the piling ship calculated in the step five by the motion response extremum of each degree of freedom of the piling ship calculated in the step six-1 in a one-to-one correspondence manner;
step six-3, calculating the operable window period of the piling ship, namely judging whether the piling ship is operable or not by using the added numerical values of the results calculated in the step six-2, and if the sum is less than or equal to 1, indicating that the condition is the operable construction window period of the piling ship; if the sum is greater than 1, the condition is not the window period of the workable construction of the piling ship;
and step six-4, taking the ratio of the number of 3 hours for which the piling ship can operate to the total number of 3 hours as the ratio of the operable window period of the piling ship.
The present invention will now be described in detail with reference to one embodiment:
referring to fig. 2 to 10, the method for calculating the ratio of the window period of operation of the piling ship according to the present invention comprises the following steps:
step one, collecting wave data of a piling ship during construction in an open sea area, wherein the sea area is located at 27.25 degrees north latitude, 121 degrees east longitude and the water depth is about 35m; the third-generation wave model MFWAM based on the global wave system of the French weather agency is adopted to obtain the wave data S of the sea area, which is 3 hours by 3 hours between 1 month and 1 day in 2016 and 12 months and 31 days in 2016 W (ω) the range of collected wave frequencies ω is (0,3.3); will collect wave data S W (ω) are listed in table 1 below:
TABLE 1
| Frequency ω (rad/s) | Wave data S W (ω)(m) |
| 0.1 | 1.560 |
| 0.2 | 1.786 |
| 0.3 | 2.102 |
| … | … |
And collect wave data S W (ω) as a function of wave frequency ω is plotted as FIG. 2;
analyzing the operation environment of the piling ship, wherein the restriction factors of whether the piling ship can operate are the rolling amplitude of the ship and the heave amplitude of the ship bow pile gripper; in the linear wave theory, the response frequency of the floating structure is the same as the excitation frequency of waves, and when the floating structure is excited at a certain frequency, the motion amplitude of the floating body and the wave amplitude form a fixed ratio, which is called an amplitude response operator RAO, and for a large floating structure, the amplitude response operator RAO is calculated by adopting commercial hydrodynamic analysis software AQWA or WAMIT based on a boundary element method; therefore, firstly, a model of the piling ship is built in hydrodynamic analysis software AQWA, and then a roll amplitude response operator RAO of the piling ship is calculated i (ω) and heave amplitude response operator RAO j (ω) the calculation results are listed in table 2 below:
TABLE 2
And roll amplitude response operator RAO of piling ship i (ω) as a function of wave frequency ω as shown in FIG. 3a, while heave amplitude response operator RAO of the piling vessel j (ω) as a function of wave frequency ω is plotted in FIG. 3b;
step three, calculating the rolling motion response energy density spectrum of the piling ship by adopting a frequency domain calculation method, namely using a rolling amplitude response operator RAO corresponding to the wave frequency omega i Square of (ω) multiplied by wave data S at the same wave frequency ω W (omega) obtaining a roll response energy density spectrumIs represented by the following formula (1) -1:
calculating heave motion response energy density spectrum of pile gripper position of pile driving ship by using frequency domain calculation method, namely using heave amplitude response operator RAO correspondent to wave frequency omega j Square of (ω) multiplied by wave data S at the same wave frequency ω W (omega) obtaining heave motion response energy density spectrumIs represented by the following formula (1) -2:
S Xj (ω)＝RAO j 2 (ω)·S W (ω) (1)-2
the calculation results are listed in table 3 below:
TABLE 3 Table 3
And responds to the roll motion of the piling ship by an energy density spectrumFIG. 4a is plotted as a function of wave frequency ω while the heave motion response of the piling vessel is plotted as an energy density spectrum +.>FIG. 4b is plotted as a function of wave frequency ω;
step four, since the motion response of the piling ship is regarded as the sum of the superposition of a plurality of harmonics, the amplitude A corresponding to the harmonics is calculated k ，A k The value of (2) is obtained by frequency domain analysis; according to the energy density spectrum of the motion response of each degree of freedom of the piling ship obtained in the third step, calculating an amplitude A corresponding to the harmonic wave of the motion response of each degree of freedom of the piling ship k The method comprises the steps of carrying out a first treatment on the surface of the The amplitude corresponding to the harmonic of the roll motion response is calculated using the following equation (2) -1:
the amplitude corresponding to the harmonic of the heave motion response is calculated using the following formula (2) -1:
Δω is 0.1, different wave frequencies ω k Obtaining different wave amplitudes A k The calculation results are listed in table 4 below:
TABLE 4 Table 4
Calculating motion responses of all degrees of freedom of the piling ship by adopting a time domain calculation method, namely calculating roll motion responses of the piling ship by adopting the following formula (3) -1:
the rolling motion response of the piling ship is formed by superposing n harmonic waves, delta omega is obtained to be 0.1 in the fourth step, and the range of wave frequency omega is obtained to be (0,3.3) in the first step, so n=3.3/0.1=33; when the wave frequency ω takes different values, then the roll motion response of the pile driving vessel is calculated by equation (3) -1:
in the formula (3) -1, ε k To obey [0,2 pi ]]The kth element in the uniformly distributed array is randomly selected by MATLAB programming during calculation to obtain an array with total n=33, and the total n=33 is listed in the following table 5:
TABLE 5
| ε 1 | 6.0172 |
| ε 2 | 1.1819 |
| ε 3 | 3.3166 |
| … | … |
| ε 33 | 5.4390 |
Substituting into (3) -1 to obtain roll motion response X of piling ship i (t)，X i (t) is a function of time t, i.e. X i (t)＝0.194sin(0.1t+6.0172)+0.290sin(0.2t+1.1819)+0.344sin(0.3t+3.3166)+...
Substituting different times t to obtain maximum values of roll motion response corresponding to the different times t, and marking the maximum values as delta max ；
Similarly, calculating heave motion response of pile gripper position of pile driving ship by adopting the following formula (3) -2:
the heave motion response of the pile gripper position of the piling ship is formed by superposition of n harmonic waves, the delta omega is obtained to be 0.1 in the fourth step, the range of the wave frequency omega is obtained to be (0,3.3) in the first step, so that n=3.3/0.1=33, and when the wave frequency omega takes different values, the heave motion response of the piling ship is calculated by the following formula (3) -2:
in the formula (3) -2, epsilon k To obey [0,2 pi ]]The kth element in the uniformly distributed array is randomly selected by MATLAB programming during calculation, so as to obtain an array with total number of n=33, and the array is listed in the following table 6:
TABLE 6
| ε 1 | 1.8746 |
| ε 2 | 4.1785 |
| ε 3 | 2.3885 |
| … | … |
| ε 33 | 5.5899 |
Substituting into (3) -2 to obtain heave motion response X of piling ship j (t)，X j (t) is a function of time t, i.e. X j (t)＝0.077sin(0.1t+1.8746)+0.116sin(0.2t+4.1785)+0.138sin(0.3t+2.3885)+...
Substituting different time t values to obtain heaveThe maximum value of motion response, denoted as eta max ；
The calculated maximum value theta of the roll motion response max And heave motion response maximum η max Listed in table 7 below:
TABLE 7
| Time t (3 hours) | Maximum value of roll motion response θ max (°) | Heave motion response maximum η max (m) |
| 10 | 6.10 | 1.75 |
| 20 | 6.89 | 2.10 |
| 30 | 8.69 | 2.25 |
| … | … | … |
Step six, calculating the operable window period of the piling ship according to the operable threshold balance of the piling ship; wherein the empirically operable threshold balance of heave given by the piling ship in the design stage is the spectrum peak period of 7s, and the sense wave height is 1.5m of random sea conditions; the empirically operable threshold balance of the roll given by the piling ship in the design stage is 7s of spectrum peak period, and the sense wave height is 0.8m of random sea conditions;
step six-1, repeating the step three to the step five, calculating the heave motion response of the pile driving ship at the pile gripper position under the wave-facing condition and the random sea condition with the sense wave height of 1.5m, and marking the heave motion response as a heave motion response extremum as H; calculating the roll motion response of the piling ship under the condition of the roll and the spectrum peak period 7s and the random sea condition with the sense wave height of 0.8m, and recording as a roll motion response extremum value as theta; the calculated heave motion response extremum H and roll motion response extremum θ are listed in table 8 below:
TABLE 8
| Time t (3 hours) | Extreme value θ (°) of roll motion response | Heave motion response extremum H (m) |
| 10 | 2.20 | 1.02 |
| 20 | 2.20 | 1.02 |
| 30 | 2.20 | 1.02 |
| … | … | … |
Drawing the calculated rolling motion response extremum and the piling ship rolling motion response maximum value calculated in the step five for 3 hours into a graph 5; drawing the calculated heave motion response extremum and the piling ship obtained by the calculation in the step five into a graph 6 every 3 hours; FIG. 5 illustrates the maximum roll motion response of the piling vessel under the back-reported wave data and the extreme roll motion response of the piling vessel under an empirically operable threshold balance based on a time-frequency domain joint calculation; FIG. 6 illustrates the heave motion response maxima of the piling vessel under the post-report wave data and the heave motion response extremum of the piling vessel under the empirically operable threshold balance based on the time-frequency domain joint calculation;
step six-2, normalizing the maximum value of the rolling motion response and the maximum value of the heave motion response, namely using the maximum value theta of the rolling motion response calculated in the step five max And heave motion response maximum η max Dividing the rolling motion response extremum theta and the heave motion response extremum H calculated in the step six-1 in a one-to-one correspondence manner, listing the calculation results in the following table 9, and drawing the calculation results into fig. 7 and 8;
TABLE 9
| Time t (3 hours) | θ max /θ | η max /H |
| 10 | 2.773 | 1.716 |
| 20 | 3.132 | 2.059 |
| 30 | 3.950 | 2.206 |
| … | … | … |
Step six-3, the rolling movement of the piling ship and the heave movement at the pile gripper have the same weight on whether the piling ship can work or not; and (3) judging whether the piling ship is operable or not by using the added values of the results calculated in the step (six-2), namely judging whether the piling ship is operable or not by adopting the calculation result of the following formula (4):
if the calculated result is less than or equal to 1, the condition is that the piling ship can work for a construction window period; if the calculation result is greater than 1, the condition is not the window period of the pile driving ship which can work; the calculation results are listed in the following table 9 and plotted as fig. 9;
TABLE 9
And step six-4, taking the ratio of the number of 3 hours for which the piling ship can operate to the total number of 3 hours as the operable window period of the piling ship, calculating the average operable window period ratio of each month in the period from 1 month in 2016 to 31 months in 12 months in 2016 according to the month, and calculating the calculation result as shown in figure 10.
The above embodiments are provided for illustrating the present invention and not for limiting the present invention, and various changes and modifications may be made by one skilled in the relevant art without departing from the spirit and scope of the present invention, and thus all equivalent technical solutions should be defined by the claims.
Claims (4)
1. A method for calculating the proportion of a window period of a pile driving ship capable of operating, which is characterized by comprising the following steps:
step one, collecting wave data S of a piling ship for every 3 hours in construction of an open sea area W (ω) as a function of wave frequency ω:
calculating an amplitude response operator RAO (omega) of each degree of freedom of the piling ship by using hydrodynamic analysis software, wherein the amplitude response operator is expressed as a function of wave frequency omega;
step three, calculating the motion response energy density spectrum of each degree of freedom of the piling ship by adopting a frequency domain calculation method, namely calculating by adopting the following formula (1):
S X (ω)＝RAO 2 (ω)·S W (ω) (1)
in the formula (1), S X (omega) is the motion response energy density spectrum of a degree of freedom of the piling ship, S W (ω) is the data obtained in step one, and RAO (ω) is the data obtained in step two;
step four, calculating an amplitude A corresponding to a harmonic wave of the motion response of each degree of freedom of the piling ship according to the motion response energy spectrum of each degree of freedom of the piling ship obtained in the step three k I.e. S is firstly X (omega) dividing the wave frequency into n equal parts corresponding to the wave frequency range Deltaomega in the frequency domain, and calculating the amplitude A by the following formula (2) k ：
In the formula (2), S X (ω k ) For the kth part frequency omega k Corresponding motion response energy density spectrum, A k The amplitude corresponding to the harmonic of the kth motion response;
calculating motion responses of all degrees of freedom of the piling ship by adopting a time domain calculation method, namely calculating by adopting the following formula (3):
in the formula (3), X (t) is the motion response of the piling ship in a certain degree of freedom, t is time, and the motion response X (t) of the piling ship in a certain degree of freedom is a function of time t, epsilon k To obey [0,2 pi ]]The kth element in the uniformly distributed array; substituting different time t values to obtain the maximum value of the freedom degree motion response;
step six, calculating the operable window period of the piling ship according to the operable threshold balance of the piling ship, and performing the following steps:
step six-1, repeating the step three to the step five, and calculating the motion response extremum of each degree of freedom of the piling ship under the operation threshold balance;
step six-2, carrying out normalization processing on the motion response maximum value of each degree of freedom of the piling ship obtained by calculation in the step five, namely dividing the motion response maximum value of each degree of freedom of the piling ship calculated in the step five by the motion response extremum of each degree of freedom of the piling ship calculated in the step six-1 in a one-to-one correspondence manner;
step six-3, calculating the operable window period of the piling ship, namely judging whether the piling ship is operable or not by using the added numerical values of the results calculated in the step six-2, and if the sum is less than or equal to 1, indicating that the condition is the operable construction window period of the piling ship; if the sum is greater than 1, the condition is not the window period of the workable construction of the piling ship;
and step six-4, taking the ratio of the number of 3 hours for which the piling ship can operate to the total number of 3 hours as the ratio of the operable window period of the piling ship.
2. A method of calculating the proportion of the operable window period of a piling vessel according to claim 1, wherein the wave data is historical wave post-report data from 3 hours to 3 hours, the wave data being obtained by a local weather station or a marine station, the wave data being constituted by storms, major swells and minor swells.
3. The method for calculating the operable window period ratio of the piling ship according to claim 1, wherein in the second step, the hydrodynamic analysis software is ocean engineering hydrodynamic analysis software based on the principle of radiation.
4. A method of calculating the operable window period ratio of a pile driving vessel according to claim 1, wherein the degrees of freedom of the pile driving vessel are at least roll of the pile driving vessel and heave of the pile carrier position of the pile carrier bow.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110881061.3A CN113673092B (en) | 2021-08-02 | 2021-08-02 | Calculation method for window period proportion of pile driving ship capable of operating |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110881061.3A CN113673092B (en) | 2021-08-02 | 2021-08-02 | Calculation method for window period proportion of pile driving ship capable of operating |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113673092A CN113673092A (en) | 2021-11-19 |
| CN113673092B true CN113673092B (en) | 2023-08-22 |
Family
ID=78541090
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110881061.3A Active CN113673092B (en) | 2021-08-02 | 2021-08-02 | Calculation method for window period proportion of pile driving ship capable of operating |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113673092B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115456232A (en) * | 2022-06-16 | 2022-12-09 | 中交第三航务工程局有限公司 | Method for predicting ratio of operable window in offshore wind power single-pile foundation construction of crane ship |
| CN116500704A (en) * | 2023-06-05 | 2023-07-28 | 中交第三航务工程局有限公司 | Operability forecasting method for offshore wind power construction ship |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106383940A (en) * | 2016-09-08 | 2017-02-08 | 大连理工大学 | Method for calculating dangerous rocking condition of LNG independent C-type cabin for ship |
| CN108399279A (en) * | 2018-01-25 | 2018-08-14 | 交通运输部天津水运工程科学研究所 | Wave effect lower terminal mooring ship loses working days analysis and control method |
| CN111709625A (en) * | 2020-06-04 | 2020-09-25 | 上海佳豪船海工程研究设计有限公司 | System and method for evaluating offshore operation rate of pile driving barge |
| AU2020102354A4 (en) * | 2020-09-21 | 2020-10-29 | Tianjin Research Institute For Water Transport Engineering.M.O.T. | Morning and early warning method for coastal port ship operation conditions |
-
2021
- 2021-08-02 CN CN202110881061.3A patent/CN113673092B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106383940A (en) * | 2016-09-08 | 2017-02-08 | 大连理工大学 | Method for calculating dangerous rocking condition of LNG independent C-type cabin for ship |
| CN108399279A (en) * | 2018-01-25 | 2018-08-14 | 交通运输部天津水运工程科学研究所 | Wave effect lower terminal mooring ship loses working days analysis and control method |
| CN111709625A (en) * | 2020-06-04 | 2020-09-25 | 上海佳豪船海工程研究设计有限公司 | System and method for evaluating offshore operation rate of pile driving barge |
| AU2020102354A4 (en) * | 2020-09-21 | 2020-10-29 | Tianjin Research Institute For Water Transport Engineering.M.O.T. | Morning and early warning method for coastal port ship operation conditions |
Non-Patent Citations (1)
| Title |
|---|
| 涌浪海域起重船运动响应特性与作业窗口分析;刘修成;华晓涛;杨秀礼;;中国港湾建设(第02期);77-81 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113673092A (en) | 2021-11-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113673092B (en) | Calculation method for window period proportion of pile driving ship capable of operating | |
| CN111950211B (en) | Seabed foundation local scouring depth determination method and system based on ocean current design flow velocity | |
| CN107729656A (en) | Dike method for prewarning risk is overflow based on the storm tide under Typhoon Process Wind | |
| CN201953134U (en) | Self supporting maritime wind measurement tower structure | |
| CN105116165A (en) | System for observing and simulating wind-wave-flow coupling field for sea-crossing bridge | |
| CN102214262A (en) | Tide predicting method | |
| CN102486633A (en) | Motion parameter optimization method for optimal energy consumption of underwater glider | |
| Ortega-Achury et al. | Hydrokinetic power review | |
| Abdraman et al. | Wind Resource Assessment in the City of N'djamena in Chad | |
| CN103243732A (en) | Installing method of marine anemometer tower base structure | |
| CN101844606B (en) | Shallow sea wind power construction system, shallow sea wind power platform and shallow sea wind power construction method | |
| CN112075369A (en) | Integrated development device of single-pile type offshore wind turbine and multi-layer vertically-displaceable aquaculture net cage | |
| CN115456232A (en) | Method for predicting ratio of operable window in offshore wind power single-pile foundation construction of crane ship | |
| Guo et al. | Integrated floating method based on four-bucket jacket foundation for offshore substations and converter stations | |
| Loughney et al. | Application of a multiple-attribute decision-analysis methodology for site selection of floating offshore wind farms off the West coast of Ireland | |
| Liščić et al. | Off shore Wind Power Plant in the Adriatic Sea: An Opportunity for the Croatian Economy | |
| CN202450543U (en) | Foundation structure of novel guide frame-type marine anemometer tower | |
| CN210315637U (en) | Offshore wind turbine single pile foundation with water purification function in pile | |
| CN112163970A (en) | Site selection method for wave energy power generation device | |
| Barthelmie | A brief review of offshore wind energy activity in the 1990's | |
| CN107503879B (en) | Method for developing ocean tidal current energy resources based on' rock drilling and tunneling | |
| Orme et al. | Analysis and comparison of support structure concepts for tidal stream turbines | |
| Hsing-Yu et al. | Analysis of Wind Wave Characteristics of Offshore Wind Farm-Case Study in Miaoli Sea Area | |
| Rizal | Measurement of Depth of Seabed and Sea Flow Rate Under SuraMadu Bridge Madura Side for Determination of Mooring Turbine Construction Sea Power Plants | |
| Puyang et al. | Preliminary design and hydrodynamic analysis on the offshore integrated anemometer mast during wet-towing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |