CN116011297B - Method and system for calculating dynamic stiffness of bow based on energy - Google Patents
Method and system for calculating dynamic stiffness of bow based on energy Download PDFInfo
- Publication number
- CN116011297B CN116011297B CN202310084840.XA CN202310084840A CN116011297B CN 116011297 B CN116011297 B CN 116011297B CN 202310084840 A CN202310084840 A CN 202310084840A CN 116011297 B CN116011297 B CN 116011297B
- Authority
- CN
- China
- Prior art keywords
- bow
- static
- curve
- internal energy
- deformation
- 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
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000003068 static effect Effects 0.000 claims abstract description 70
- 238000006243 chemical reaction Methods 0.000 claims description 18
- 230000000694 effects Effects 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000010276 construction Methods 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 2
- 230000006835 compression Effects 0.000 description 3
- 238000007906 compression Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011960 computer-aided design Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003116 impacting effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method and a system for calculating the dynamic stiffness of a ship bow based on energy, which establish the relationship between static stiffness and dynamic stiffness of the ship bow under different structures, different tonnages and impact speeds, and further convert the static stiffness into the dynamic stiffness under different impact energies through direct description of the static stiffness, so that various parameters for determining the dynamic stiffness are omitted, and the method is simpler and more definite in form and easy to master for engineers.
Description
Technical Field
The invention relates to the technical field of computer aided design using a finite element method, in particular to a method and a system for calculating dynamic stiffness of a bow based on energy.
Background
In the process of the bridge collision, the rigidity characteristic of the bow directly influences the collision result. The stiffness of the bow changes dynamically during the impact, even for the same vessel, once the vessel mass, the impact speed change, the stiffness changes again. Therefore, how to characterize the dynamic stiffness of the bow is of great significance for researching the crashworthiness of the bridge, and further for designing a reasonable and efficient crashproof structure in a targeted manner.
At present, related documents are used for researching the dynamic stiffness curves and providing corresponding simplified formulas, specifically, a ship collision stiffness wall analysis is carried out by establishing refined finite element models of different ships to obtain a plurality of dynamic stiffness curves with different tonnages and different speeds, mathematical statistics is carried out on the curves to obtain a plurality of parameters with different tonnages and speeds, and the parameters are further fitted to obtain the dynamic stiffness curve. However, up to 36 mathematical statistical parameters per vessel are too complex to be mastered by engineers.
Disclosure of Invention
Aiming at the defects existing in the prior art, the invention provides a method and a system for calculating the dynamic stiffness of a bow based on energy. Various parameters for determining the dynamic stiffness can be omitted, and the dynamic stiffness is simpler and more definite in form and is easy to master by engineers. The specific technical scheme is as follows:
in a first aspect, there is provided an energy-based method for calculating the dynamic stiffness of a bow, comprising:
respectively acquiring various basic time course data of the foreship with different tonnages under the pseudo static force effect and different impact effects;
determining a bow static stiffness curve F (D), a bow quasi static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D), a bow dynamic stiffness curve unloading section and a bow impact internal energy-deformation curve E (D) of the ship with different tonnages through the acquired various basic time course data;
determining a quasi-static-dynamic conversion coefficient according to a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D) and a bow impact internal energy-deformation curve E (D) of ships with different tonnagesDetermining an unloading section slope k according to the bow dynamic stiffness curve unloading sections of the ships with different tonnages;
constructing a bow quasi-static internal energy-deformation function e (d) based on the bow quasi-static internal energy-deformation curve e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple ;
Combined with the static energy-deformation function e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple Unloading segment slope k and pseudo static-dynamic conversion coefficientDetermination of the impact internal energy of the foreship deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple ;
By said bow impact internal energy-deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple And calculating the dynamic stiffness of the bow.
With reference to the first aspect, in a first implementation manner of the first aspect, a finite element analysis method is adopted to obtain various basic time course data of foreship with different tonnages under the pseudo static force effect and different impact effects.
With reference to the first aspect, in a second implementation manner of the first aspect, determining the unloading segment slope k includes:
counting average slopes of bow dynamic stiffness curve unloading sections of bows with different tonnages at different impact speeds;
and determining the slope k of the unloading section according to the average slope of the unloading section of the bow dynamic stiffness curve.
With reference to the first aspect, in a third implementation manner of the first aspect, a static stiffness reduction analytical function f (d) of the foreship is constructed simple Comprising:
based on the bow quasi-static internal energy-deformation curve e (d), establishing a bow quasi-static internal energy-deformation function e (d) simple ;
e(d) simple =ax b ;
Static internal energy-deformation function e (d) of bow based on law of conservation of energy simple Solving to obtain a and b;
quasi-static internal energy-deformation function e (d) for bow simple Deriving and combining the obtained a and b to establish a static stiffness simplified analytical function f (d) of the bow simple ;
f(d) simple =a·bx b-1 。
In a second aspect, there is provided an energy-based bow dynamic stiffness calculation system comprising:
the data acquisition module is configured to acquire various basic time course data of the foreship with different tonnages under the quasi-static force effect and under different impact effects respectively;
the data processing module is configured to determine a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D), a bow dynamic stiffness curve unloading section and a bow impact internal energy-deformation curve E (D) of the ship with different tonnages through the acquired various basic time course data;
the parameter determining module is configured to determine a quasi-static-dynamic conversion coefficient according to a static stiffness curve F (D), a quasi-static internal energy-deformation curve E (D), a dynamic stiffness curve F (D) and an impact internal energy-deformation curve E (D) of the bowDetermining an unloading section slope k according to the bow dynamic stiffness curve unloading sections of the ships with different tonnages;
a function construction module configured to construct a bow quasi-static internal energy-deformation function e (d) based on the bow quasi-static internal energy-deformation curve e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple ;
A curve determination module configured to combine the bow pseudo-static internal energy-deformation function e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple Unloading segment slope k and pseudo static-dynamic conversion coefficientDetermination of the impact internal energy of the foreship deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple ;
A stiffness calculation module configured to pass through the bow impact internal energy-deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple And calculating the dynamic stiffness of the bow.
With reference to the second aspect, in a first implementation manner of the second aspect, the data acquisition module acquires various basic time course data of foreship with different tonnages under the pseudo static force effect and different impact effects by adopting a finite element analysis method.
With reference to the second aspect, in a second implementation manner of the second aspect, the parameter determining module includes:
the slope statistics unit is configured to count average slopes of unloading sections of the bow dynamic stiffness curves of the bows with different tonnages at different impact speeds;
and the slope determining unit is configured to determine the slope k of the unloading section according to the average slope of the unloading section of the bow dynamic stiffness curve.
With reference to the second aspect, in a third implementation manner of the second aspect, the function building module includes:
a function construction unit configured to establish a bow quasi-static internal energy-deformation function e (d) based on the bow quasi-static internal energy-deformation curve e (d) simple ;
e(d) simple =ax b ;
A function solving unit configured to simulate the static internal energy-deformation function e (d) of the ship bow based on the law of conservation of energy simple Solving to obtain a and b;
a function building unit configured to simulate a static internal energy-deformation function e (d) for the stem simple Deriving and combining the obtained a and b to establish a static stiffness simplified analytical function f (d) of the bow simple ;
f(d) simple =a·bx b-1 。
The beneficial effects are that: by adopting the method and the system for calculating the dynamic stiffness of the bow based on energy, the relationship between the static stiffness and the dynamic stiffness of the bow under different structures, different tonnages and impact speeds is established, and the static stiffness is directly described so as to be further converted into the dynamic stiffness under different impact energies, so that various parameters for determining the dynamic stiffness are omitted, and the method and the system are simpler and more definite in form and are easy to master for engineers.
Drawings
In order to more clearly illustrate the embodiments of the present invention, the drawings that are required to be used in the embodiments will be briefly described. Throughout the drawings, the elements or portions are not necessarily drawn to actual scale.
FIG. 1 is a flow chart of a method for calculating dynamic stiffness of a bow based on energy according to an embodiment of the present invention;
FIG. 2 is a graph showing the results of comparison of F (D), E (D), F (D), E (D) for a 3000DWT bow at an impact velocity of 1 m/s;
FIG. 3 is a graph showing the results of comparison of F (D), E (D), F (D), E (D) for a 3000DWT bow at an impact velocity of 1.5 m/s;
FIG. 4 is a graph showing the results of comparison of F (D), E (D), F (D), E (D) for a 3000DWT bow at an impact velocity of 2 m/s;
FIG. 5 is a graph showing the results of comparison of F (D), E (D), F (D), E (D) for a 3000DWT bow at an impact velocity of 2.5 m/s;
FIG. 6 is a graph showing the results of comparison of F (D), E (D), F (D), E (D) for a 3000DWT bow at an impact velocity of 3 m/s;
FIG. 7 is a graph showing the results of comparison of F (D), E (D), F (D), E (D) for a 3000DWT bow at an impact velocity of 3.5 m/s;
FIG. 8 is a graph showing the results of comparison of F (D), E (D), F (D), E (D) for a 3000DWT bow at an impact velocity of 4 m/s;
FIG. 9 is a graph showing the results of comparison of F (D), E (D), F (D), E (D) for a 3000DWT bow at an impact velocity of 4.5 m/s;
FIG. 10 is a graph showing the results of comparison of F (D), E (D), F (D), E (D) for a 3000DWT bow at an impact velocity of 5 m/s;
FIG. 11 is a graphical illustration of bow stiffness unloading of a 3000DWT bow at different impact speeds;
FIG. 12 is a schematic view of a bow quasi-static internal energy-deformation curve of a 3000DWT bow;
FIG. 13 is a simplified analytical curve schematic of static stiffness of a 3000DWT bow;
FIG. 14 is a graph of the dynamic stiffness of the bow of a 3000DWT bow at an impact velocity of 3 m/s;
FIG. 15 is a graph of the dynamic stiffness of the bow of a 3000DWT bow at an impact velocity of 4 m/s;
FIG. 16 is a graph of the dynamic stiffness of the bow of a 3000DWT bow at an impact velocity of 5 m/s;
FIG. 17 is a system block diagram of an energy-based bow dynamic stiffness calculation system.
Detailed Description
Embodiments of the technical scheme of the present invention will be described in detail below with reference to the accompanying drawings. The following examples are only for more clearly illustrating the technical aspects of the present invention, and thus are merely examples, and are not intended to limit the scope of the present invention.
A flow chart of a method of energy-based calculation of the dynamic stiffness of the bow of a ship, as shown in fig. 1, the method comprising:
step 1, respectively acquiring various basic time course data of a bow with different tonnages under the pseudo static force effect and different impact effects;
step 2, determining a bow static force stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D), a bow dynamic stiffness curve unloading section and a bow impact internal energy-deformation curve E (D) of the ship with different tonnages through the acquired various basic time course data;
step 3, determining a quasi-static-dynamic conversion coefficient according to a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D) and a bow impact internal energy-deformation curve E (D) of the ship with different tonnagesDetermining an unloading section slope k according to the bow dynamic stiffness curve unloading sections of the ships with different tonnages;
step 4, constructing a bow quasi-static internal energy-deformation function e (d) based on the bow quasi-static internal energy-deformation curve e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple ;
Step 5, combining the bow pseudo-static internal energy-deformation function e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple Unloading segment slope k and pseudo static-dynamic conversion coefficientDetermination of the impact internal energy of the foreship deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple ;
Step 6, impacting the internal energy-deformation curve E (D) through the bow simple And a bow dynamic stiffness curve F (D) simple And calculating the dynamic stiffness of the bow.
Specifically, first, various basic time course data of the foreship under the action of a constant quasi-static force with different tonnages, such as a foreship quasi-static force deformation time course curve, a quasi-static force compression time course curve and a foreship quasi-static force internal energy time course curve, can be acquired. And the foreship with different tonnages impacts various time-course data under the rigid wall at different impact speeds, such as a foreship impact deformation time-course curve, an impact force time-course curve and a foreship impact internal energy time-course curve.
In this embodiment, various basic time-course data may be obtained by using a finite element analysis method, specifically:
firstly, a high-precision finite element model and a rigid wall model of ships with different tonnages can be established.
And then, carrying out finite element analysis based on the constructed high-precision finite element model and the rigid wall model, and compressing the high-precision finite element model of each bow at a constant speed of 0.1m through the rigid wall model, thereby obtaining bow quasi-static deformation time course curves, quasi-static compression time course curves and bow quasi-static internal energy time course curves corresponding to ships with different tonnages.
And then the rigid wall model is impacted positively by the high-precision finite element model of the ship with different tonnages at different impact speeds of 1 m/s-5 m/s, so that a bow impact deformation time course curve, an impact force time course curve and a bow impact internal energy time course curve corresponding to the ship with different tonnages are obtained.
And then, eliminating the intermediate quantity t of the acquired various basic time course data to obtain a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D), a bow dynamic stiffness curve unloading section and a bow impact internal energy-deformation curve E (D) corresponding to the ships with different tonnages.
Then, a static stiffness curve F (D), a quasi-static internal energy-deformation curve e (D), a dynamic stiffness curve F (D) and an impact internal energy-deformation curve are carried out on the bow corresponding to the ships with different tonnagesCurve E (D) is compared to determine the pseudo-static-dynamic conversion coefficient
In this example, a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D), and a bow impact internal energy-deformation curve E (D) corresponding to a bow with tonnage in the range of 800DWT to 50000DWT are obtained in total. By comparison, the static force is equal to the impact result when the impact speed is lower than 2m/s, which shows that the threshold strain rate speed of the bridge collision is 2m/s. When the speed exceeds 2m/s, the result of the impact is substantially equal to 1.2 times of the static result, and the pseudo-static-dynamic conversion coefficient can be determinedThe method comprises the following steps:
because of the large data volume, the comparison results are not conveniently displayed in full, and the embodiment provides a comparison of the results of a 3000DWT bow, and the comparison results are shown in figures 2-10.
The stiffness curve unloading section is positioned at the end of energy exchange, and the ship moves away from the bridge pier in the opposite direction, so that the unloading section has less influence on the collision process. Meanwhile, the comparison of the curves of the unloading sections of the ships can be described by using straight lines, and the slopes are not greatly different, as shown in figure 11, so that the average slope of the unloading section of the dynamic stiffness curve of the ships with different tonnages under different impact speeds can be counted as the slope k of the unloading section,
then, based on the bow quasi-static internal energy-deformation curve e (d), a bow quasi-static internal energy-deformation function e (d) can be constructed simple And a static stiffness reduced analytical function f (d) for the bow simple 。
Specifically, based on the static internal energy-deformation curve e (d), a bow pseudo static internal energy-deformation function e (d) is established simple =ax b Based on energy conservationAnd (5) solving a and b constantly. The result of solving the 3000DWT bow is given here, as shown in fig. 12, with the corresponding values: a=6.488e6, b=1.25. Then pair e (d) simple =ax b Deriving, bringing a and b into, thus obtaining the static rigidity simplified analysis type f (d) of the bow simple See fig. 13.
Thereafter, the pseudo-static internal energy-deformation function e (d) of the stem can be combined simple And a static stiffness reduced analytical function f (d) for the bow simple Unloading segment slope k and pseudo static-dynamic conversion coefficientDetermination of the impact internal energy of the foreship deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple 。
Established bow impact internal energy-deformation curve E (D) simple The following is shown:
established bow dynamic stiffness curve F (D) simple The following is shown:
wherein x is the speed, D m The maximum collision depth of the bow is given, and c is a constant.
Finally, the internal energy of impact-deformation curve E (D) can be passed through the stem simple And a bow dynamic stiffness curve F (D) simple And calculating the dynamic stiffness of the bow.
Specifically, the impact internal energy-deformation curve E (D) of the ship bow is based on the initial impact kinetic energy of the ship simple Determining the maximum depth D of the bow during the collision m From the bow dynamic stiffness curve F (D) simple Determining the impact force F (D) m ). The unloading segment curve is obtained from the statistically determined average slope k. Taking a 3000DWT bow as an example, the simplified power steel at different speeds is givenThe degree curve is shown in fig. 14-16.
A system block diagram of an energy-based bow dynamic stiffness calculation system as shown in fig. 17, the system comprising:
the data acquisition module is configured to acquire various basic time course data of the foreship with different tonnages under the quasi-static force effect and under different impact effects respectively;
the data processing module is configured to determine a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D), a bow dynamic stiffness curve unloading section and a bow impact internal energy-deformation curve E (D) of the ship with different tonnages through the acquired various basic time course data;
a parameter determination module configured to determine a quasi-static-dynamic conversion coefficient according to a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D) and a bow impact internal energy-deformation curve E (D) of the ship with different tonnagesDetermining an unloading section slope k according to the bow dynamic stiffness curve unloading sections of the ships with different tonnages;
a function construction module configured to construct a bow quasi-static internal energy-deformation function e (d) based on the bow quasi-static internal energy-deformation curve e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple ;
A curve determination module configured to combine the bow pseudo-static internal energy-deformation function e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple Unloading segment slope k and pseudo static-dynamic conversion coefficientDetermination of the impact internal energy of the foreship deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple ;
A stiffness calculation module configured to pass through the bow impact internal energy-deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple Calculating the bowDynamic stiffness.
The computing system comprises a data acquisition module, a data processing module, a parameter determination module, a function construction module, a curve determination module and a rigidity calculation module. The data acquisition module can acquire various basic time course data of foreship with different tonnages under the quasi-static force effect and different impact effects by adopting a finite element analysis method.
Specifically, the data acquisition module can establish a high-precision finite element model and a rigid wall model of ships with different tonnages. And performing finite element analysis based on the constructed high-precision finite element model and the rigid wall model, and compressing the high-precision finite element model of each bow at a constant speed of 0.1m through the rigid wall model, thereby obtaining bow quasi-static deformation time course curves, quasi-static compression time course curves and bow quasi-static internal energy time course curves corresponding to ships with different tonnages.
The data acquisition module positively impacts the rigid wall model at different impact speeds of 1 m/s-5 m/s through high-precision finite element models of the foreship with different tonnages, so that a foreship impact deformation time course curve, an impact force time course curve and a foreship impact internal energy time course curve corresponding to the ships with different tonnages are obtained.
The data processing module can remove intermediate quantities of various basic time course data acquired by the data acquisition module, so that a bow static stiffness curve F (D), a bow quasi static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D), a bow dynamic stiffness curve unloading section slope and a bow impact internal energy-deformation curve E (D) corresponding to ships with different tonnages are obtained.
The parameter determination module can determine the quasi-static-dynamic conversion coefficient according to a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D) and a bow impact internal energy-deformation curve E (D) of ships with different tonnagesAnd determining the slope k of the unloading section according to the unloading sections of the bow dynamic stiffness curves of the ships with different tonnages.
Specifically, the parameter determination module includes a curveThe device comprises a comparison unit, a slope statistics unit and a slope determination unit. The curve comparison unit can compare the static force stiffness curve F (D), the quasi-static internal energy-deformation curve E (D), the dynamic stiffness curve F (D) and the impact internal energy-deformation curve E (D) of the bow corresponding to ships with different tonnages so as to determine the quasi-static-dynamic conversion coefficient
The slope statistics unit can count average slopes of bow dynamic stiffness curve unloading sections of ships with different tonnages at different impact speeds. The slope determination unit may determine the unloading segment slope k from an average slope of the unloading segment of the bow dynamic stiffness curve.
The function construction module can construct the bow quasi-static internal energy-deformation function e (d) based on the bow quasi-static internal energy-deformation curve e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple 。
Specifically, the function construction module comprises a function construction unit, a function solving unit and a function establishing unit. Wherein the function construction unit can establish a bow pseudo-static internal energy-deformation function e (d) based on the static internal energy-deformation curve e (d) simple =ax b . The function solving unit may solve a, b in the bow pseudo-static internal energy-deformation function based on conservation of energy. The result of solving the 3000DWT bow is given here, with the corresponding values: a=6.488e6, b=1.25. The function building unit can be applied to e (d) simple =ax b Deriving and bringing a, b into, to obtain a static stiffness reduction function f (d) of the bow simple 。
The curve determination module can be combined with a bow pseudo-static internal energy-deformation function e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple Unloading segment slope k and pseudo static-dynamic conversion coefficientDetermination of the impact internal energy of the foreship deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple 。
Determined bow impact internal energy-deformation curve E (D) simple The following is shown:
determined bow dynamic stiffness curve F (D) simple The following is shown:
wherein x is the speed, D m The maximum collision depth of the bow is given, and c is a constant.
The rigidity calculation module can pass through the bow impact internal energy-deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple And calculating the dynamic stiffness of the bow.
Specifically, the rigidity calculation module can pass through a bow impact internal energy-deformation curve E (D) according to the initial impact kinetic energy of the ship simple Determining the maximum depth D of the bow during the collision m Then pass through the bow dynamic stiffness curve F (D) simple Determining the impact force F (D) m ). The unloading segment curve is obtained from the statistically determined average slope k.
The above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention, and are intended to be included within the scope of the appended claims and description.
Claims (6)
1. The method for calculating the dynamic stiffness of the foreship based on energy is characterized by comprising the following steps of:
respectively acquiring various basic time course data of the foreship with different tonnages under the pseudo static force effect and different impact effects;
determining a bow static stiffness curve F (D), a bow quasi static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D), a bow dynamic stiffness curve unloading section and a bow impact internal energy-deformation curve E (D) of the ship with different tonnages through the acquired various basic time course data;
determining a quasi-static-dynamic conversion coefficient according to a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D) and a bow impact internal energy-deformation curve E (D) of ships with different tonnagesDetermining an unloading section slope k according to the bow dynamic stiffness curve unloading sections of the ships with different tonnages;
constructing a bow quasi-static internal energy-deformation function e (d) based on the bow quasi-static internal energy-deformation curve e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple Constructing a static stiffness simplified analytical function f (d) of the bow simple Comprising:
based on the bow quasi-static internal energy-deformation curve e (d), establishing a bow quasi-static internal energy-deformation function e (d) simple ;
e(d) simple =ax b ;
Static internal energy-deformation function e (d) of bow based on law of conservation of energy simple Solving to obtain a and b;
quasi-static internal energy-deformation function e (d) for bow simple Deriving and combining the obtained a and b to establish a static stiffness simplified analytical function f (d) of the bow simple ;
f(d) simple =a·bx b-1 ;
Combined with the static energy-deformation function e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple Unloading segment slope k and pseudo static-dynamic conversion coefficientDetermination of the impact internal energy of the foreship deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple The built-up bow impact internal energy-deformation curve E (D) simple The following is shown:
established bow dynamic stiffness curve F (D) simple The following is shown:
wherein x is the speed, D m The maximum collision depth of the bow is set, and c is a constant;
by said bow impact internal energy-deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple And calculating the dynamic stiffness of the bow under different tonnages and different speeds.
2. The method for calculating the dynamic stiffness of the foreship based on the energy according to claim 1, wherein various basic time course data of foreship with different tonnages under the pseudo static force and different impact actions are obtained by adopting a finite element analysis method.
3. The energy-based bow dynamic stiffness calculation method of claim 1, wherein determining the unloading segment slope k comprises:
counting average slopes of bow dynamic stiffness curve unloading sections of bows with different tonnages at different impact speeds;
and determining the slope k of the unloading section according to the average slope of the unloading section of the bow dynamic stiffness curve.
4. An energy-based bow dynamic stiffness computing system, comprising:
the data acquisition module is configured to acquire various basic time course data of the foreship with different tonnages under the quasi-static force effect and under different impact effects respectively;
the data processing module is configured to determine a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D), a bow dynamic stiffness curve unloading section and a bow impact internal energy-deformation curve E (D) of the ship with different tonnages through the acquired various basic time course data;
a parameter determination module configured to determine a quasi-static-dynamic conversion coefficient according to a bow static stiffness curve F (D), a bow quasi-static internal energy-deformation curve E (D), a bow dynamic stiffness curve F (D) and a bow impact internal energy-deformation curve E (D) of the ship with different tonnagesDetermining an unloading section slope k according to the bow dynamic stiffness curve unloading sections of the ships with different tonnages;
a function construction module configured to construct a bow quasi-static internal energy-deformation function e (d) based on the bow quasi-static internal energy-deformation curve e (d) simple And a static stiffness reduced analytical function f (d) for the bow simple The function construction module includes:
a function construction unit configured to establish a bow quasi-static internal energy-deformation function e (d) based on the bow quasi-static internal energy-deformation curve e (d) simple ;
e(d) simple =ax b ;
A function solving unit configured to simulate the static internal energy-deformation function e (d) of the ship bow based on the law of conservation of energy simple Solving to obtain a and b;
a function building unit configured to simulate a static internal energy-deformation function e (d) for the stem simple Deriving and combining the obtained a and b to establish a static stiffness simplified analytical function f (d) of the bow simple ;
f(d) simple =a·bx b-1 ;
A curve determination module configured to combineStatic energy-deformation function e (d) of bow simple And a static stiffness reduced analytical function f (d) for the bow simple Unloading segment slope k and pseudo static-dynamic conversion coefficientDetermination of the impact internal energy of the foreship deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple The built-up bow impact internal energy-deformation curve E (D) simple The following is shown:
established bow dynamic stiffness curve F (D) simple The following is shown:
wherein x is the speed, D m The maximum collision depth of the bow is set, and c is a constant;
a stiffness calculation module configured to pass through the bow impact internal energy-deformation curve E (D) simple And a bow dynamic stiffness curve F (D) simple And calculating the dynamic stiffness of the bow.
5. The energy-based bow dynamic stiffness calculation system according to claim 4, wherein the data acquisition module acquires various basic time course data of the bow under the quasi-static force and different impact forces by adopting a finite element analysis method.
6. The energy-based bow dynamic stiffness calculation system of claim 4, wherein the parameter determination module comprises:
the slope statistics unit is configured to count average slopes of unloading sections of the bow dynamic stiffness curves of the bows with different tonnages at different impact speeds;
and the slope determining unit is configured to determine the slope k of the unloading section according to the average slope of the unloading section of the bow dynamic stiffness curve.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310084840.XA CN116011297B (en) | 2023-02-01 | 2023-02-01 | Method and system for calculating dynamic stiffness of bow based on energy |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310084840.XA CN116011297B (en) | 2023-02-01 | 2023-02-01 | Method and system for calculating dynamic stiffness of bow based on energy |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116011297A CN116011297A (en) | 2023-04-25 |
CN116011297B true CN116011297B (en) | 2024-01-30 |
Family
ID=86029942
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310084840.XA Active CN116011297B (en) | 2023-02-01 | 2023-02-01 | Method and system for calculating dynamic stiffness of bow based on energy |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116011297B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106372345A (en) * | 2016-09-06 | 2017-02-01 | 广船国际有限公司 | Method and device for analyzing stress of ship docking blocks |
CN107256289A (en) * | 2017-05-18 | 2017-10-17 | 吉林大学 | The method for building up of car crass reduced parameter FEM model |
CN109670252A (en) * | 2018-12-25 | 2019-04-23 | 中南大学 | A kind of head vehicle contracting mould construction method and head vehicle contracting mould based on power and stiffnes s equivalent |
-
2023
- 2023-02-01 CN CN202310084840.XA patent/CN116011297B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106372345A (en) * | 2016-09-06 | 2017-02-01 | 广船国际有限公司 | Method and device for analyzing stress of ship docking blocks |
CN107256289A (en) * | 2017-05-18 | 2017-10-17 | 吉林大学 | The method for building up of car crass reduced parameter FEM model |
CN109670252A (en) * | 2018-12-25 | 2019-04-23 | 中南大学 | A kind of head vehicle contracting mould construction method and head vehicle contracting mould based on power and stiffnes s equivalent |
WO2020135425A1 (en) * | 2018-12-25 | 2020-07-02 | 中南大学 | Head vehicle miniature construction method based on force and rigidity equivalence and head vehicle miniature |
Non-Patent Citations (2)
Title |
---|
船桥碰撞时船艏的刚度简化模型;冒一锋;刘伟庆;周叮;方海;庄勇;;南京工业大学学报(自然科学版)(第05期);全文 * |
转筒式防撞设施防船撞性能研究;陈巍;耿波;沈锐利;郑植;;船舶力学(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116011297A (en) | 2023-04-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109670252B (en) | Force and rigidity equivalence-based head turning shrinkage model construction method and head turning shrinkage model | |
CN104298869B (en) | A kind of fluid structurecoupling Numerical prediction method of elastic hydrofoil | |
Ma et al. | Shock-capturing non-hydrostatic model for fully dispersive surface wave processes | |
Hu et al. | Numerical simulation and experiment on dam break problem | |
WO2020187113A1 (en) | Method for predicting emission amount of single vehicle | |
Fife et al. | Multiscaling in the presence of indeterminacy: wall-induced turbulence | |
Huxley et al. | TUFLOW GPU-best practice advice for hydrologic and hydraulic model simulations | |
Takami et al. | Efficient FORM-based extreme value prediction of nonlinear ship loads with an application of reduced-order model for coupled CFD and FEA | |
Ma et al. | Nonlinear waves in a floating thin elastic plate, predicted by a coupled SPH and FEM simulation and by an analytical solution | |
Alberello et al. | Three dimensional velocity field underneath a breaking rogue wave | |
CN116011297B (en) | Method and system for calculating dynamic stiffness of bow based on energy | |
CN106768763A (en) | A kind of dash-board injury method for quick | |
Yang et al. | A robust solver for incompressible high-Reynolds-number two-fluid flows with high density contrast | |
CN110530589B (en) | Beam bridge impact coefficient detection method considering continuous traffic flow coupling effect | |
Bu et al. | Unified viscous and potential prediction method for the coupled motion of damaged ship and floodwater in calm water | |
CN109063401B (en) | Method for realizing structural contact force optimization by adopting equivalent static displacement | |
Zhang et al. | Investigations on the hydroelastic slamming of deformable wedges by using the smoothed particle element method | |
CN106777850A (en) | A kind of automobile component design method based on simplified assessment | |
Attili et al. | Scaling approaches and scale effects in wave–flexible structure interaction | |
CN111625987B (en) | Turbulence model generation method based on deep learning | |
Sun et al. | A new rational-based optimal design strategy of ship structure based on multi-level analysis and super-element modeling method | |
Hao | A song of atmosphere and ocean: modeling and simulation-based study of nonlinear waves | |
Peng et al. | Study on the Influence of Ship Speed and Ship Weight on Ship-Bridge Collision Force | |
CN116882327B (en) | Array OWC hydrodynamic analysis method on complex terrain | |
CN113392450B (en) | Method, device and system for evaluating health state of railway pier |
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 |