CN113673092A - Method for calculating operable window period proportion of pile driving barge - Google Patents
Method for calculating operable window period proportion of pile driving barge Download PDFInfo
- Publication number
- CN113673092A CN113673092A CN202110881061.3A CN202110881061A CN113673092A CN 113673092 A CN113673092 A CN 113673092A CN 202110881061 A CN202110881061 A CN 202110881061A CN 113673092 A CN113673092 A CN 113673092A
- Authority
- CN
- China
- Prior art keywords
- pile driving
- calculating
- freedom
- ship
- piling
- 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.)
- Granted
Links
Images
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 method for calculating the operable window period proportion of a pile driving barge, which comprises the following steps: collecting wave data of a piling ship in 3 hours per hour during construction in an open sea area; step two, calculating amplitude response operators of all degrees of freedom of the piling ship by adopting hydrodynamic analysis software; calculating the motion response energy density spectrum of each degree of freedom of the pile driving ship by adopting a frequency domain calculation method; step four, calculating the wave 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 density spectrum of each degree of freedom of the piling ship; step five, calculating the motion response of each degree of freedom of the piling ship by adopting a time domain calculation method; and step six, calculating the operable window period of the pile driving ship according to the operable threshold calibration of the pile driving ship. The invention adopts a time-frequency domain combined 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 method for calculating an operable window period proportion of a pile driving ship.
Background
In offshore wind power construction projects in China, the piling ship occupies a more important position. Since the first offshore wind farm in China-the east bridge wind farm is constructed, a large number of high-pile bearing platform foundation type offshore wind turbines which take pile driving ships as construction equipment appear. The high pile cap foundation is a foundation type widely applied to the field of port and channel engineering, generally, a foundation pile group is formed by 6-8 steel pipe inclined piles with the diameter of 2-3 m, a cap poured by high-performance concrete is arranged at the top of the foundation pile group to connect the pile group, and therefore, the construction of the high pile cap generally needs to adopt a professional pile driving ship and a concrete mixing ship to respectively finish pile sinking of the pile group and concrete pouring construction. As the piling equipment of the professional piling barge in China is arranged on the fixed pile frame of the bow, if the construction of 6-8 inclined piles uniformly arranged in 360 degrees horizontally is to be finished, the bow of the piling barge needs to be laid at a plurality of angles, and the heading of the ship cannot be freely selected.
At present, offshore wind power plants in China are mostly concentrated in a water depth range of 30m to 50m, the offshore distance can reach 75km farthest, the environmental conditions of the offshore wind power plants in deep and heavy wave sea areas are worse than those in near sea and sea areas, the operational period of a pile driving ship is difficult to judge based on experience, and the operational window period of the pile driving ship needs to be calculated by adopting a quantitative analysis method.
In the world, the professional piling ship is mainly used for offshore area engineering, the stormy waves are small, and the operable window period is judged mainly based on experience, so that the research on the movement response analysis of the piling ship at home and abroad is relatively less. However, the pile driving ship as a floating engineering ship can work under complex sea conditions and research methods of other floating engineering ships can be used for reference. 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, namely ship heave, roll and pitch. The general idea of the analysis method in this respect is: firstly, determining a risk event which may cause an accident or construction failure, determining critical limiting parameters of the risk event, then performing numerical simulation on offshore construction under various different sea conditions by adopting a floating body dynamics simulation method, and determining the sea conditions of operable and inoperable by combining a simulation result and the critical limiting parameters.
Numerical simulation of offshore construction is mainly carried out by calculating 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 the six-degree-of-freedom motion of the ship. Through decades of development, the dynamic simulation analysis of the floating body is mature at present, particularly the simulation of the construction period movement with small relative movement amplitude is realized, and the calculation result based on commercial software such as AQWA, SESAM and the like has high reliability and reference value.
The key step in calculating the operable window period of the pile driving ship is floating body hydrodynamics calculation, the main methods of the existing floating body hydrodynamics simulation analysis comprise time domain calculation and frequency domain calculation, and the two methods respectively have the advantages and the disadvantages: the conditions for the frequency domain calculation are harsh, while the time domain calculation is very time consuming.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a method for calculating the operable window period proportion of a piling ship.
The purpose of the invention is realized as follows: a method for calculating the operable window period proportion of a pile driving barge comprises the following steps:
step one, collecting wave data S of every 3 hours of construction of a pile driving ship in an open sea areaW(ω), the wave data being expressed as a function of the wave frequency ω;
step two, calculating an amplitude response operator RAO (omega) of each degree of freedom of the pile driving ship by adopting hydrodynamic analysis software, wherein the amplitude response operator is expressed as a function of the wave frequency omega;
step three, calculating the motion response energy density spectrum of each degree of freedom of the pile driving ship by adopting a frequency domain calculation method, namely calculating by adopting the following formula (1):
SX(ω)＝RAO2(ω)·SW(ω) (1)
in the formula (1), SX(omega) is the energy density spectrum of the motion response of the piling vessel in a certain degree of freedom, SW(omega) is the data obtained in the first step, and RAO (omega) is the data obtained in the second step;
step four, calculating the amplitude A corresponding to the harmonic wave of the motion response of each degree of freedom of the piling ship according to the motion response energy density spectrum of each degree of freedom of the piling ship obtained in the step threekImmediately first of all, SX(omega) is divided into n equal parts corresponding to the wave frequency range delta omega in the frequency domain, and the wave amplitude A is calculated by the following formula (2)k：
In the formula (2), SX(ωk) Is the k-th frequency omegakCorresponding energy density spectrum of movement response, AkAmplitude corresponding to the harmonic of the kth motion response;
step five, calculating the motion response of each degree 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 a certain degree of freedom of the piling ship, t is time, and the motion response X (t) of the certain degree of freedom of the piling ship is shown as a function of the time t, and epsilonkTo obey [0,2 π]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 pile driving vessel according to the operable threshold balance of the pile driving vessel, and performing the following steps:
step six-1, repeating the step three to the step five, and calculating the motion response extreme value of each degree of freedom of the piling ship under the operable threshold balance;
step six-2, the motion response maximum value of each degree of freedom of the piling ship obtained by calculation in the step five is subjected to normalization processing, namely the motion response maximum value of each degree of freedom of the piling ship calculated in the step five is divided by the motion response extreme value of each degree of freedom of the piling ship calculated in the step six-1 in a one-to-one correspondence mode;
step six-3, calculating the operable window period of the pile driving barge, namely judging whether the pile driving barge can operate by using the numerical value obtained by adding 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 pile driving barge; if the sum is more than 1, the condition is not the construction window period of the pile driving barge in operation;
and step six-4, taking the ratio of the number of 3 hours which the piling barge can work to the total number of 3 hours as the proportion of the window period which the piling barge can work.
The method for calculating the operable window period proportion of the pile driving vessel comprises the step one, wherein the wave data is reported after 3 hours of historical waves, the wave data is obtained by a local meteorological office or a local oceanographic office, and the wave data comprises wind waves, main swell and secondary swell.
In the method for calculating the ratio of the operable window period of the pile driving vessel, when the step two is performed, the hydrodynamic analysis software is ocean engineering hydrodynamic analysis software based on a radiation principle.
The method for calculating the ratio of the operable window periods of the pile driving barge, wherein the degree of freedom of the pile driving barge is at least the rolling of the pile driving barge and the heaving of the position of a bow embracing pile device of the pile driving barge.
The method for calculating the operable window period proportion of the pile driving barge has the following effects:
1. according to the practical operation condition and the environmental condition characteristics of the pile driving barge, the operable limiting environmental condition of the pile driving barge is accurately analyzed, and on the basis of the limiting environmental condition, the condition that whether the pile driving barge can operate in a certain period in the future can be predicted, so that the waiting time of the pile driving barge is effectively reduced, and the construction efficiency of the pile driving barge is improved;
2. different from other environment condition window calculation, the calculation model constructed in the invention is simpler and more convenient, and because waves appear in the form of multi-peak spectrum mixed waves under most conditions for open sea, the energy density spectrum S of motion response is not onlyX(omega) or the wave spectrum S of an irregular waveW(omega) in most cases, the assumed requirement of Gaussian signals is difficult to meet, so that the motion response value of the pile driving ship cannot be directly calculated by using a full frequency domain method, and at the moment, the energy density spectrum S of the calculated motion response can be usedXAnd (omega) is converted into a time domain sequence X (t), and the motion response corresponding to different transcendental probabilities is calculated in the time domain, namely the time-frequency domain combined calculation method provided by the invention. The method has the main advantages that the floating body dynamic equation is prevented from being directly solved in the time domain, the calculation time is greatly reduced, and the harsh assumed conditions depending on full-frequency-domain calculation are avoided.
Drawings
FIG. 1 is a flow chart of a method of calculating the workable window period ratio of the pile driving vessel of the present invention;
FIG. 2 is a graph of wave data collected as a function of wave frequency during the performance of the steps of the present invention;
FIG. 3a is a plot of the roll amplitude response operator of the piling vessel obtained during step two of the present invention as a function of wave frequency;
FIG. 3b is a diagram showing the heave amplitude response operator of the piling vessel obtained in step two of the present invention as a function of the wave frequency;
FIG. 4a is a plot of the roll motion response energy density spectrum of the pile driving vessel obtained when step three of the present invention is performed as a function of wave frequency;
FIG. 4b is a plot of the heave motion response energy density spectrum of the pile driving vessel obtained when step three of the present invention is performed as a function of wave frequency;
FIG. 5 is a comparison graph of the maximum roll motion response per 3 hours of the piling vessel obtained when step six-1 of the present invention is performed and the extreme roll motion response under the workable threshold balance;
FIG. 6 is a graph comparing the maximum heave motion response and the extremum of the operational threshold constant heave motion response for 3 hours for a piling vessel, obtained when step six-1 of the present invention is performed;
FIG. 7 is a graph showing the normalized results of the response of the piling vessel to roll motion for 3 hours after performing step six-2 of the present invention;
FIG. 8 is a graph showing the normalized result of the heave motion response maximum for 3 hours for the piling vessel obtained when step six-2 of the present invention is performed;
FIG. 9 is a schematic view of a pile driving vessel operational window period resulting from performing step six-3 of the present invention;
figure 10 is a schematic diagram of the proportion of the window period of operability of the pile driving vessel resulting from carrying out step six-4 of the present invention.
Detailed Description
The invention will be further explained with reference to the drawings.
Referring to fig. 1, the method for calculating the operable window period ratio of the pile driving vessel of the present invention includes the following steps:
step one, collecting wave data S of every 3 hours of construction of a pile driving ship in an open sea areaW(ω) the wave data expressed as a function of the wave frequency ω, the wave data being 3 hour by hour historical post-wave report data, the wave data being obtained by the local meteorological office or the oceanic office, the wave data being comprised of storms, primary swells and secondary swells;
step two, calculating an amplitude response operator RAO (omega) of each degree of freedom of the pile driving ship by adopting hydrodynamic analysis software, wherein the amplitude response operator is expressed as a function of the wave frequency omega; the hydrodynamic Analysis software is ocean engineering hydrodynamic Analysis software based on a radiation principle, such as AQWA, WAMIT (Wave Analysis MIT Analysis software for calculating interaction between a zero-navigational-speed floating structure and waves) and the like; the freedom degree of the pile driving ship is at least the rolling of the pile driving ship and the heaving of the position of a bow pile embracing device of the pile driving ship;
step three, calculating the motion response energy density spectrum of each degree of freedom of the pile driving ship by adopting a frequency domain calculation method, namely calculating by adopting the following formula (1):
SX(ω)＝RAO2(ω)·SW(ω) (1)
in the formula (1), SX(omega) is the energy density spectrum of the motion response of the piling vessel in a certain degree of freedom, SW(omega) is the data obtained in the first step, and RAO (omega) is the data obtained in the second step;
step four, calculating the amplitude A corresponding to the harmonic wave of the motion response of each degree of freedom of the piling ship according to the motion response energy density spectrum of each degree of freedom of the piling ship obtained in the step threekImmediately first of all, SX(omega) is divided into n equal parts corresponding to the wave frequency range delta omega in the frequency domain, and the wave amplitude A is calculated by the following formula (2)k：
In the formula (2), SX(ωk) Is the k-th frequency omegakCorresponding energy density spectrum of movement response, AkAmplitude corresponding to the harmonic of the kth motion response;
step five, calculating the motion response of each degree 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 a certain degree of freedom of the piling ship, t is time, and the motion response X (t) of the certain degree of freedom of the piling ship is shown as a function of the time t, and epsilonkTo obey [0,2 π]The important role of the k-th element in the uniformly distributed array is to ensure the randomness of X (t), so that epsilonkAlso 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 pile driving vessel according to the operable threshold balance of the pile driving vessel, and performing the following steps:
step six-1, repeating the step three to the step five, and calculating the motion response extreme value of each degree of freedom of the piling ship under the operable threshold balance;
step six-2, the motion response maximum value of each degree of freedom of the piling ship obtained by calculation in the step five is subjected to normalization processing, namely the motion response maximum value of each degree of freedom of the piling ship calculated in the step five is divided by the motion response extreme value of each degree of freedom of the piling ship calculated in the step six-1 in a one-to-one correspondence mode;
step six-3, calculating the operable window period of the pile driving barge, namely judging whether the pile driving barge can operate by using the numerical value obtained by adding 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 pile driving barge; if the sum is more than 1, the condition is not the construction window period of the pile driving barge in operation;
and step six-4, taking the ratio of the number of 3 hours which the piling barge can work to the total number of 3 hours as the proportion of the window period which the piling barge can work.
The invention will now be described in detail with an embodiment:
referring to fig. 2 to 10, the method for calculating the operable window period ratio of the pile driving vessel of the present invention includes the following steps:
collecting wave data of a pile driving ship during construction in an open sea area, wherein the sea area is located at 27.25 degrees of north latitude and 121 degrees of east longitude, and the water depth is about 35 m; adopting a third-generation wave model MFWAM based on a French weather bureau global wave system to obtain wave data S of the sea area by 3 hours between 1/2016 and 12/31/2016W(ω), the collected wave frequency ω ranges from (0, 3.3); collecting wave data SW(ω) is shown in Table 1 below:
TABLE 1
| Frequency omega (rad/s) | Wave data SW(ω)(m) |
| 0.1 | 1.560 |
| 0.2 | 1.786 |
| 0.3 | 2.102 |
| … | … |
And collecting wave data SW(ω) is plotted as a function of the wave frequency ω in 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 heaving amplitude of the position of the bow embracing pile device; in the linear wave theory, the response frequency of the floating structure is the same frequency as the excitation frequency of waves, when the floating structure is excited corresponding to a certain frequency, the motion amplitude of the floating body and the wave amplitude are in a fixed ratio and are called as an amplitude response operator RAO, and for a large floating structure, the amplitude response operator RAO adopts commercial hydrodynamic analysis software AQWA or WAMIT based on a boundary element method to calculate; therefore, a model of the piling ship is established in the water power analysis software AQWA, and then the rolling amplitude response operator RAO of the piling ship is calculatedi(ω) and heave amplitude response operator RAOj(ω), the calculation results are shown in table 2 below:
TABLE 2
And will pile the shipRoll amplitude response operator RAOi(omega) is plotted as figure 3a as a function of the wave frequency omega, while the heave amplitude response operator RAO of the pile driving vessel is usedj(ω) is plotted as a function of the wave frequency ω in FIG. 3 b;
thirdly, calculating the rolling motion response energy density spectrum of the pile driving ship by adopting a frequency domain calculation method, namely using a rolling amplitude response operator RAO corresponding to the wave frequency omegaiThe square of (omega) is multiplied by the wave data S at the same wave frequency omegaW(omega) to obtain the energy density spectrum of roll motion response
Formula (1) -1:
calculating the heave motion response energy density spectrum of the position of the pile embracing device of the pile driving ship by adopting a frequency domain calculation method, namely using a heave amplitude response operator RAO corresponding to the wave frequency omegajThe square of (omega) is multiplied by the wave data S at the same wave frequency omegaW(omega) to obtain the energy density spectrum of the heave motion response
Formula (1) -2:
SXj(ω)＝RAOj 2(ω)·SW(ω) (1)-2
the results of the calculations are shown in table 3 below:
TABLE 3
And responding the rolling motion of the piling ship to an energy density spectrum
The relationship as a function of the wave frequency omega is plotted in FIG. 4a, while the heaving motion of the piling vessel is plottedResponse energy density spectrum
The functional relationship with the wave frequency ω is plotted in fig. 4 b;
step four, because the motion response of the pile driving ship is regarded as the sum of the superposition of a plurality of harmonics, the amplitude A corresponding to the plurality of harmonics is calculatedk，AkThe value of (a) is obtained by frequency domain analysis; calculating the amplitude A corresponding to the harmonic wave of the motion response of each degree of freedom of the pile driving barge according to the energy density spectrum of the motion response of each degree of freedom of the pile driving barge obtained in the step threek(ii) a The amplitude corresponding to the harmonic of the roll motion response is calculated using the following equation (2) -1:
calculating the amplitude corresponding to the harmonic of the heave motion response using the following equation (2) -1:
Δ ω is 0.1, different wave frequencies ωkObtain different amplitudes AkThe results of the calculations are shown in table 4 below:
TABLE 4
Step five, calculating the motion response of each degree of freedom of the pile driving ship by adopting a time domain calculation method, namely calculating the roll motion response of the pile driving ship by adopting the following formula (3) -1:
the rolling motion response of the pile driving ship is formed by overlapping n harmonics, the step four shows that delta omega is 0.1, the step one shows that the wave frequency omega is in the range of (0, 3.3), and therefore n is 3.3/0.1 is 33; when the wave frequency ω takes different values, then the rolling motion response of the pile driving vessel is calculated by equation (3) -1:
in the formula (3) -1,. epsilonkTo obey [0,2 π]The k-th element in the uniformly distributed array is randomly selected by MATLAB programming during calculation, and an array with a total number of n-33 is obtained and is listed in the following table 5:
TABLE 5
| ε1 | 6.0172 |
| ε2 | 1.1819 |
| ε3 | 3.3166 |
| … | … |
| ε33 | 5.4390 |
Calculation of the roll response X of a pile-driving vessel in the formula (3) -1i(t)，Xi(t) is a function of time t, i.e. Xi(t)＝0.194sin(0.1t+6.0172)+0.290sin(0.2t+1.1819)+0.344sin(0.3t+3.3166)+...
Substituting different time t to obtain the maximum roll motion response value corresponding to different time t, and recording the maximum roll motion response value as deltamax；
Similarly, the heave motion response of the position of the pile gripper of the piling ship is calculated by adopting the following formula (3) -2:
the heave motion response of the position of the pile driver of the pile driving vessel is formed by superposition of n harmonics, the step four shows that delta omega is 0.1, the step one shows that the range of the wave frequency omega is (0, 3.3), therefore, n is 3.3/0.1 is 33, and when the wave frequency omega is different values, the heave motion response of the pile driving vessel is calculated by the formula (3) -2:
in the formula (3) -2,. epsilonkTo obey [0,2 π]The k-th element in the uniformly distributed array is randomly selected by MATLAB programming during calculation, and an array with a total number of n-33 is obtained and is listed in table 6 below:
TABLE 6
| ε1 | 1.8746 |
| ε2 | 4.1785 |
| ε3 | 2.3885 |
| … | … |
| ε33 | 5.5899 |
Substituting formula (3) -2 to obtain heave motion response X of piling shipj(t)，Xj(t) is a function of time t, i.e. Xj(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 the maximum value of the response of the heave motion, and recording the maximum value as etamax；
The calculated maximum value theta of the rolling motion responsemaxAnd heave motion response maximum ηmaxAre shown in Table 7 below:
TABLE 7
| Time t (3 hours) | Roll motion response maximum θ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 pile driving barge according to the operable threshold balance of the pile driving barge; the experimental operational threshold balance of heave given by the pile driving ship in the design stage is the random sea state with the spectrum peak period of 7s and the sense wave height of 1.5 m; the empirical operational threshold balance of rolling given by the pile driving vessel in the design stage is the random sea state with the spectrum peak period of 7s and the sense wave height of 0.8 m;
step six-1, repeating the step three to the step five, and calculating the heaving motion response of the pile driving ship at the position of the pile embracing device under the random sea condition that the spectrum peak period is 7s under the condition of head wave and the sense wave height is 1.5m, wherein the heaving motion response is used as a heaving motion response extreme value and is marked as H; calculating the rolling motion response of the piling ship under random sea conditions of a spectrum peak period of 7s and a sense wave height of 0.8m under the condition of rolling waves, and taking the rolling motion response as a rolling motion response extreme value which is marked 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 roll response theta (°) | 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 extreme value and the piling ship rolling motion response maximum value calculated in the fifth step into a graph 5 every 3 hours; drawing the calculated heave motion response extreme value and the calculated heave motion response maximum value of the piling ship for 3 hours into a figure 6; FIG. 5 shows the maximum response value of the rolling motion of the piling ship under the post-reporting wave data and the extreme response value of the rolling motion of the piling ship under the empirical workable threshold balance, which are obtained based on the time-frequency domain joint calculation; FIG. 6 shows the heaving motion response maximum value of the piling ship under the back wave data and the heaving motion response extreme value of the piling ship under the empirical workable threshold balance, which are obtained based on the time-frequency domain combined calculation;
step six-2, normalizing the maximum value of the response of the rolling motion and the maximum value of the response of the heaving motion, namely the maximum value theta of the response of the rolling motion calculated in the step fivemaxAnd heave motion response maximum ηmaxDividing the roll motion response extreme value theta and the heave motion response extreme value H calculated in the step six-1 one by one, listing the calculation results in the following table 9, and drawing into a graph 7 and a graph 8;
TABLE 9
| Time t (3 hours) | θmax/θ | ηmax/ |
| 10 | 2.773 | 1.716 |
| 20 | 3.132 | 2.059 |
| 30 | 3.950 | 2.206 |
| … | … | … |
Step six-3, the influence of the rolling motion of the pile driving barge and the heaving motion at the pile embracing device on whether the pile driving barge can work has the same weight; and (3) judging whether the pile driving barge can work or not by using the value obtained by adding the results calculated in the step six-2, namely judging whether the pile driving barge can work or not by using the calculation result of the following formula (4):
if the calculation result is less than or equal to 1, the condition is that the piling ship can work and construct the window period; if the calculation result is larger than 1, the condition is not the construction window period of the pile driving barge in operation; the results of the calculations are set forth in Table 9 below and plotted in FIG. 9;
TABLE 9
Step six-4, taking the ratio of the number of 3 hours in which the pile driving vessel can operate to the total number of 3 hours as the operable window period of the pile driving vessel, and calculating the average operable window period ratio of each month from 1/2016 to 12/31/2016 by month in parts per month, with the calculation result shown in fig. 10.
The above embodiments are provided only for illustrating the present invention and not for limiting the present invention, and those skilled in the art can make various changes and modifications without departing from the spirit and scope of the present invention, and therefore all equivalent technical solutions should also fall within the scope of the present invention, and should be defined by the claims.
Claims (4)
1. A method for calculating the operable window period proportion of a pile driving barge is characterized by comprising the following steps:
step one, collecting wave data S of every 3 hours of construction of a pile driving ship in an open sea areaW(ω) the wave data being expressed as a function of the wave frequency ω:
step two, calculating an amplitude response operator RAO (omega) of each degree of freedom of the pile driving ship by adopting hydrodynamic analysis software, wherein the amplitude response operator is expressed as a function of the wave frequency omega;
step three, calculating the motion response energy density spectrum of each degree of freedom of the pile driving ship by adopting a frequency domain calculation method, namely calculating by adopting the following formula (1):
SX(ω)＝RAO2(ω)·SW(ω) (1)
in the formula (1), SX(omega) is the energy density spectrum of the motion response of the piling vessel in a certain degree of freedom, SW(ω)RAO (omega) is the data obtained in the step two;
step four, calculating the amplitude A corresponding to the harmonic wave of the motion response of each degree of freedom of the piling ship according to the motion response energy density spectrum of each degree of freedom of the piling ship obtained in the step threekImmediately first of all, SX(omega) is divided into n equal parts corresponding to the wave frequency range delta omega in the frequency domain, and the wave amplitude A is calculated by the following formula (2)k：
In the formula (2), SX(ωk) Is the k-th frequency omegakCorresponding energy density spectrum of movement response, AkAmplitude corresponding to the harmonic of the kth motion response;
step five, calculating the motion response of each degree 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 a certain degree of freedom of the piling ship, t is time, and the motion response X (t) of the certain degree of freedom of the piling ship is shown as a function of the time t, and epsilonkTo obey [0,2 π]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 pile driving vessel according to the operable threshold balance of the pile driving vessel, and performing the following steps:
step six-1, repeating the step three to the step five, and calculating the motion response extreme value of each degree of freedom of the piling ship under the operable threshold balance;
step six-2, the motion response maximum value of each degree of freedom of the piling ship obtained by calculation in the step five is subjected to normalization processing, namely the motion response maximum value of each degree of freedom of the piling ship calculated in the step five is divided by the motion response extreme value of each degree of freedom of the piling ship calculated in the step six-1 in a one-to-one correspondence mode;
step six-3, calculating the operable window period of the pile driving barge, namely judging whether the pile driving barge can operate by using the numerical value obtained by adding 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 pile driving barge; if the sum is more than 1, the condition is not the construction window period of the pile driving barge in operation;
and step six-4, taking the ratio of the number of 3 hours which the piling barge can work to the total number of 3 hours as the proportion of the window period which the piling barge can work.
2. The method for calculating the operable window period proportion of the pile driving vessel as recited in claim 1, wherein in the step one, the wave data is 3-hour-by-3-hour historical post-wave report data, the wave data is obtained by a local meteorological office or a local oceanographic office, and the wave data is composed of wind waves, main swells and secondary swells.
3. The method for calculating the ratio of the operable window period of the pile driving vessel as recited in claim 1, wherein in the second step, the hydrodynamic analysis software is ocean engineering hydrodynamic analysis software based on the radiation-induced principle.
4. The method for calculating the proportion of the workable window period of the pile driving barge according to claim 1, wherein the degree of freedom of the pile driving barge is at least a roll of the pile driving barge and a heave of the position of a stem embracing piler of the pile driving barge.
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 true CN113673092A (en) | 2021-11-19 |
| CN113673092B 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) |
Cited By (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 |
|---|
| 刘修成;华晓涛;杨秀礼;: "涌浪海域起重船运动响应特性与作业窗口分析", 中国港湾建设, no. 02, pages 77 - 81 * |
Cited By (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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113673092B (en) | 2023-08-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Guo et al. | Review of integrated installation technologies for offshore wind turbines: Current progress and future development trends | |
| CN107729656A (en) | Dike method for prewarning risk is overflow based on the storm tide under Typhoon Process Wind | |
| CN113673092A (en) | Method for calculating operable window period proportion of pile driving barge | |
| WO2022036752A1 (en) | Current design flow rate-based seabed foundation local scour depth determination method and system | |
| CN108413926B (en) | High-precision measurement method for underwater topography elevation of pile foundation of offshore wind farm | |
| CN102214262A (en) | Tide predicting method | |
| CN105116165A (en) | System for observing and simulating wind-wave-flow coupling field for sea-crossing bridge | |
| CN109472065A (en) | Forecast the three-dimensional numerical method of Low Speed Full Ships Added Resistance under severe sea condition | |
| CN115795970A (en) | Method for predicting lateral displacement of pile foundation of high-pile wharf | |
| Soares | Renewable energies offshore | |
| CN115456232A (en) | Method for predicting ratio of operable window in offshore wind power single-pile foundation construction of crane ship | |
| Nagai et al. | GPS buoy and seabed installed wave gauge application to offshore tsunami observation | |
| Guo et al. | Integrated floating method based on four-bucket jacket foundation for offshore substations and converter stations | |
| CN203472445U (en) | Walking type semi-submersible amphibious engineering ship | |
| Duan et al. | Wind characteristics and wind energy assessment in the Barents Sea based on ERA-Interim reanalysis | |
| 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 | |
| Lubis et al. | Characteristics of waters during transitional season, senimba waters | |
| Zheng et al. | Establishment of wave climate datasets: Case study for the Maritime Silk Road | |
| Meng et al. | A new integrated wet-towing method of offshore wind turbine with tri-bucket foundation | |
| KR102540462B1 (en) | System to establish a transportation plan for floating wind turbine based on forecasting of marine conditions | |
| CN116552731A (en) | Method for determining main scale of marine wind power construction ship | |
| Xiao et al. | New Approach for Penetration Resistance Calculation of Suction Bucket Foundation | |
| Govoni et al. | Caisson foundations for competitive offshore wind farms in Italy |
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 |