CN114943630B - Method for calculating emergency escape and evacuation time of personnel under ship water inflow condition - Google Patents
Method for calculating emergency escape and evacuation time of personnel under ship water inflow condition Download PDFInfo
- Publication number
- CN114943630B CN114943630B CN202111288697.3A CN202111288697A CN114943630B CN 114943630 B CN114943630 B CN 114943630B CN 202111288697 A CN202111288697 A CN 202111288697A CN 114943630 B CN114943630 B CN 114943630B
- Authority
- CN
- China
- Prior art keywords
- personnel
- evacuation
- ship
- calculating
- time
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000004364 calculation method Methods 0.000 claims abstract description 33
- 239000012530 fluid Substances 0.000 claims description 19
- 238000006073 displacement reaction Methods 0.000 claims description 17
- 230000001133 acceleration Effects 0.000 claims description 14
- 230000005484 gravity Effects 0.000 claims description 12
- 230000003068 static effect Effects 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 claims description 6
- 230000035484 reaction time Effects 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000005096 rolling process Methods 0.000 description 2
- 239000013535 sea water Substances 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 230000003213 activating effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
-
- 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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Mathematical Optimization (AREA)
- Development Economics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Fluid Mechanics (AREA)
- Computer Security & Cryptography (AREA)
- Computing Systems (AREA)
- Mathematical Analysis (AREA)
- Educational Administration (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for calculating personnel emergency escape and evacuation time under the condition of ship water inflow. Step 1: based on CFD, simulating water inflow of the ship, acquiring a ship body movement angle, and providing data for subsequent calculation, wherein the ship body movement angle is roll and pitch; step 2: and (3) calculating the emergency escape and evacuation time of the personnel under the condition of water inflow based on the data in the step (1). The invention is used for solving the problem that the damage water inflow of the ship has influence on the emergency escape and evacuation time of personnel.
Description
Technical Field
The invention relates to the field of escape and evacuation time, in particular to a method for calculating personnel emergency escape and evacuation time under the condition of ship water inflow.
Background
With the development of ocean energy development technology, the demand of various countries for ships is increasing, and the ship is becoming larger and has become a trend of the development of the ship industry. The increase of marine structures on the sea surface increases the probability of damage accidents of the ship, and has great potential safety hazards. Once the outer plate of the large ship is damaged, seawater flows in, so that the cargo is damaged, the personnel are panicked, and economic loss and casualties are caused.
The water intake process after the breakage of the ship is a very complicated problem, which includes exchange of liquid inside and outside the cabin, escape and compression of air inside the cabin, water flowing to other cabins through the internal opening, and the like. According to the statistical analysis of IMO, the process of breakage of water inflow over time can be divided into three main phases: and a transient water inlet stage, a continuous water inlet stage and a stable stage.
Instantaneous water ingress may be defined as the period from when a breach occurs to when the breach chamber fills with water, typically for a period of less than a few tens of seconds, and as seawater flows through the breach chamber, the water ingress within the chamber may cause a significant intermediate transient roll angle during extreme times due to the particular size and location of certain breach openings, causing the vessel to rapidly tip to the side being breached. The occurrence of large amounts of ballast and free water in the hull due to breakage may result in the vessel capsizing before static equilibrium is reached and the duration of the capsizing process is very short. There is no research significance for the case where capsizing occurs in a very short time, as opposed to the evacuation of people. In the case of a ship which is not tilted in a short time due to the broken water inflow after the occurrence of the breakage, the ship gradually becomes stable and reaches a new balance after the ship is severely rocked, at the moment, the ship is balanced at a certain inclination angle or slightly fluctuates within a certain inclination angle range, the ship is gradually stable and does not tilt, but people on the ship need to be evacuated rapidly due to the fact that a plurality of unstable factors exist on the sea, the people can evacuate according to a planned path, the people are panic due to the severe rocking of the water inflow early stage of the ship, even the influence of the people on the escape time of the people due to the fact that the people cannot walk and the influence of the ship inclination caused by the water inflow of the ship on the escape time of the people are caused, and therefore the calculation of the escape time of the people is also changed.
According to the invention, the CFD technology is used for simulating the ship damage water inlet process by utilizing FLUENT software, and the influence of ship damage water inlet on the emergency escape and evacuation time of personnel is considered, so that the calculation method of the emergency escape and evacuation time of personnel under the water inlet condition is constructed.
Disclosure of Invention
The invention provides a method for calculating the emergency escape and evacuation time of personnel under the condition of ship water inflow, which is used for solving the problem that the emergency escape and evacuation time of the personnel is influenced by the damaged water inflow of a ship.
The invention is realized by the following technical scheme:
the method for calculating the emergency escape and evacuation time of the personnel under the condition of water inflow of the ship comprises the following steps:
step 1: based on CFD, simulating water inflow of the ship, acquiring a ship body movement angle, and providing data for subsequent calculation, wherein the ship body movement angle is roll and pitch;
step 2: and (3) calculating the emergency escape and evacuation time of the personnel under the condition of water inflow based on the data in the step (1).
Further, the step 1 specifically includes the following steps:
step 1.1: establishing a CFD numerical model; the CFD numerical model only considers a mass conservation equation and a momentum conservation equation;
step 1.2: solving the CFD numerical model in the step 1.1, dividing a fluid solving domain into a limited number of mutually adjacent control bodies, and applying a conservation equation to each control body;
step 1.3: capturing the interface of water and gas by using a VOF method;
step 1.4: the six-degree-of-freedom solver and the dynamic grid technology are utilized to process the ship movement;
step 1.5: monitoring a ship movement angle and obtaining a ship movement curve; setting calculation time and step length, and performing iterative calculation
Further, the mass conservation equation and the momentum conservation equation in the step 1.1 are specifically,
the mass conservation equation is:
the conservation of momentum equation is:
wherein: u, v, w are velocity components of velocity vector v along x, y, z; x, Y, Z is the unit mass force applied to the microcell in the x, y, z directions; p fluid pressure; a coefficient of dynamic viscosity;
the step 1.2 specifically comprises the following conservation equations:
wherein: phi represents a general variable; v represents a control volume; Γ represents the generalized diffusion coefficient; s represents a generalized source term.
Further, the step 1.3 is specifically that each ratio F of the fluid volume in the grid to the grid volume is calculated; when F is 1, the fluid in the grid is the main phase fluid; when F is 0, the grid is filled with secondary phase fluid; when 0 < F < 1: representing a multiphase flow boundary within the mesh; the interface of the water and air is denoted as,
r w +r a =0
wherein: r is (r) w Representing the volume fraction of water; r is (r) a Representing the volume fraction of air.
Further, the step 1.4 is specifically to calculate the displacement and the rotation angle about the center of gravity, and use the stress and the moment of the object; the displacement control equation for the center of gravity is as follows:
wherein: m represents a rigid body mass;indicating the gravity center stress; />Representing the center of gravity displacement acceleration;
angular acceleration:
wherein: l represents an inertial tensor;representing a moment tensor; />Representing the rigid body angular velocity; />The method is the conversion from an inertial coordinate system to a satellite coordinate system; after the angular acceleration and the displacement acceleration are calculated, the dynamic grid calculation can update the rigid body position using the angular velocity and the displacement velocity.
Further, the step 2 specifically includes the following steps:
step 2.1: planning the shortest escape route of personnel, and obtaining the number N of the branch routes of each escape route and the inclination angle of the ship body;
step 2.2: calculating the walking speed v and the specific flow F of the personnel on the evacuation route S ;
Step 2.3: calculating personnel evacuation travel time T;
step 2.4: the time REST required for the evacuation of the person in the case of water intake is calculated, rest=r+t.
Further, the step 2.1 of planning the shortest escape route of the personnel to escape utilizes Dijstra algorithm to plan the shortest escape route and obtain the number N of the branch route.
Further, step 2.2 calculates the walking speed v and the specific flow F of the person on the evacuation route S The method comprises the following steps:
personnel walking speed v:
v=v 0 ·k 1 ·k 2 '·k″ 2
wherein: v 0 Indicating the initial speed, k of the person 1 Representing the influence of personnel density on walking speed, the MSC gives personnel density D asRelationship with the speed of travel of the person by a reduction factork 1 The form represents:
wherein: n is the total number of people in the path; l is the path length; w (W) C Is the static width of the path;
reduction factor k 'in the case of transverse tilting' 2 The calculation formula is, wherein, ψ is the transverse inclination:
Reduction factor k″ in the case of a longitudinal inclination 2 The calculation formula is shown as a formula in which,the pitch angle is:
Thus, the person specific flow rate F S Obtained by interpolation method, the calculation formula is:
in the middle of: n is the total number of people in the path; l is the path length; w (W) C Is the static width of the path.
Further, the specific method for calculating the personnel evacuation travel time T in the step 2.3 is as follows:
T=γ·t l
wherein: gamma is a correction coefficient; l (L) j In the path scheme, i is an evacuated people stream, and j branch paths are included in the evacuation route; l (L) y For the length of the y-th branch path in the evacuation route; v y For the speed of the person in the y-th branch path in the evacuation route; n (N) y For the total number of people in the y-th branch path in the evacuation route; f (F) sy A specific flow in the y-th branch path in the evacuation route; w (W) cy For the static width in the y-th branch path in the evacuation route.
Further, the specific formula for calculating the evacuation time in the step 2.4 is as follows:
wherein: r is the reaction time; t is the person evacuation travel time.
The beneficial effects of the invention are as follows:
according to the invention, the CFD technology is utilized to simulate the water inflow process of the ship, the influence of severe shaking in the early stage of water inflow of the ship on the escape time of personnel is considered, and the influence of ship inclination caused by water inflow on the walking speed of personnel is also considered on the basis of calculating the evacuation time by the MSC.
Has great value and social significance for ship design, emergency situation setting, evacuation plan making and even promotion of the improvement of the safety of shipping industry.
Drawings
FIG. 1 is a diagram of a broken water intake process of a ship.
Fig. 2 is a flow chart of the evacuation time under the water inflow calculation in step 2.
FIG. 3 is a flow chart of the FLUENT software simulation of water intake for a ship.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
According to the invention, the ship is not overturned in a short time due to broken water inflow, but gradually approaches to a stable state after violent shaking of the ship, and a new balanced water inflow scene is achieved. The method aims at providing the method, on the basis of calculating the escape and evacuation time of the personnel based on the MSC, the influence of panic or even incapability of walking of the personnel on the escape time caused by severe shaking at the early stage of ship damage and water inflow and the influence of ship inclination on the escape time caused by ship water inflow are also considered.
The method for calculating the emergency escape and evacuation time of the personnel under the condition of water inflow of the ship comprises the following steps:
step 1: based on CFD and FLUENT software, simulating water inflow of the ship, acquiring a ship body movement angle, and providing data for subsequent calculation, wherein the ship body movement angle is roll and pitch; step 1.1-1.5 is an illustration of the step 1 of using CFD to simulate the water intake of a ship to obtain the hull movement angle;
step 2: and (3) calculating the emergency escape and evacuation time of the personnel under the condition of water inflow based on the data in the step (1).
The method for calculating the emergency escape and evacuation time of personnel under the condition of water inflow of a ship specifically comprises the following steps:
step 1.1: establishing a CFD numerical model; the CFD numerical model only considers a mass conservation equation and a momentum conservation equation;
step 1.2: solving the CFD numerical model in the step 1.1, wherein discretization is needed for solving the partial differential equation, and CFD software utilizes a Finite Volume Method (FVM); dividing the fluid solution domain into a limited number of mutually adjacent control volumes, applying a conservation equation to each control volume;
step 1.3: capturing the interface of water and gas by using a VOF method;
step 1.4: processing hull motion using a Six degree of freedom (Six-DOF) solver and a dynamic grid technique;
step 1.5: monitoring a ship movement angle through a 'monitor' option in FLUENT, and activating 'print to control', 'plot' and 'file', so as to obtain a ship movement curve; setting calculation time and step length, and clicking the 'calculation' button to perform iterative calculation
A method for calculating the emergency escape and evacuation time of personnel under the condition of ship water inflow, wherein the mass conservation equation and the momentum conservation equation in the step 1.1 are specifically as follows,
the mass conservation equation is:
the conservation of momentum equation is:
wherein: u, v, w are velocity components of velocity vector v along x, y, z; x, Y, Z is the unit mass force applied to the microcell in the x, y, z directions; p fluid pressure; a coefficient of dynamic viscosity;
the step 1.2 specifically comprises the following conservation equations:
wherein: phi represents a general variable; v represents a control volume; Γ represents the generalized diffusion coefficient; s represents a generalized source term.
The invention adopts a separation type solution to correspond to a 'Pressure Based' solver in Fluent, is Based on a classical SIMPLE algorithm and is applicable to solving incompressible flow and micro-compressible flow.
The method for calculating the emergency escape and evacuation time of personnel under the condition of ship water inflow comprises the following steps of 1.3, determining a free surface by researching a fluid and grid volume ratio function F in a grid unit by a VOF model, wherein the principle of the VOF method is as follows: each calculating a ratio F of fluid volume within the grid to grid volume; when F is 1, the fluid in the grid is the main phase fluid; when F is 0, the grid is filled with secondary phase fluid; when 0 < F < 1: representing a multiphase flow boundary within the mesh; the interface of the water and air is denoted as,
r w +r a =0
wherein: r is (r) w Representing the volume fraction of water; r is (r) a Representing the volume fraction of air.
A calculation method of personnel emergency escape and evacuation time under the condition of ship water inflow, the step 1.4 is that the dynamic grid technology comprises an active grid and a passive grid; the active grid is used for defining the movement of the object boundary before calculation, and presetting the movement condition of the object and the grid by writing UDF or related settings in software. Calculating displacement and rotation angle about the center of gravity, and using the force and moment of the object; the displacement control equation for the center of gravity is as follows:
Wherein: m represents a rigid body mass;indicating the gravity center stress; />Representing the center of gravity displacement acceleration;
angular acceleration:
wherein: l represents an inertial tensor;representing a moment tensor; />Representing the rigid body angular velocity; />The method is the conversion from an inertial coordinate system to a satellite coordinate system; after the angular acceleration and the displacement acceleration are calculated, the dynamic grid calculation can update the rigid body position using the angular velocity and the displacement velocity.
B represents a transformation matrix:the method is the conversion from an inertial coordinate system to a satellite coordinate system;
here, C χ =cos(χ),S χ =sin (χ). Phi, theta, phi denote Euler angles that rotate along the x, y, z axes. After the angular acceleration and the displacement acceleration are calculated by the above calculation, the dynamic grid calculation can update the rigid body position by using the angular velocity and the displacement velocity.
The method for calculating the emergency escape and evacuation time of personnel under the condition of water inflow of a ship specifically comprises the following steps:
step 2.1: planning the shortest escape route of personnel, and obtaining the number N of the branch routes of each escape route and the inclination angle of the ship body;
step 2.2: calculating the walking speed v and the specific flow F of the personnel on the evacuation route S ;
Step 2.3: calculating personnel evacuation travel time T;
step 2.4: the time REST required for the evacuation of the person in the case of water inflow is calculated, rest=r+t, R representing the reaction time.
A calculation method of emergency escape and evacuation time of personnel under the condition of ship water inflow is disclosed, wherein in the step 2.1, the shortest escape path of personnel is planned, the Dijstra algorithm is utilized to plan the shortest escape path, and the number N of people on branch paths is obtained.
From the viewpoints of safety margin and simplified calculation, the inclination angle is taken as the inclination angle in the stable stage, the specific inclination angle is obtained in the step 1, and the water inflow process diagram is shown in the attached figure 1.
A method for calculating personnel emergency escape and evacuation time under the condition of ship water inflow, which comprises the following steps of 2.2, calculating personnel walking speed v and specific flow F on an evacuation route S The method comprises the following steps:
personnel walking speed v (m/s):
v=v 0 ·k 1 ·k 2 '·k″ 2
wherein: v 0 Representing the initial speed (m/s), k of the person 1 Representing the influence of personnel density on walking speed, the MSC gives the personnel densityRelationship with the speed of travel of the person to reduceCoefficient k 1 The form represents:
wherein: n is the total number of people (p) in the path; l is the path length (m); w (W) C Is the static width (m) of the path.
k 2 The invention summarizes relevant literature, statistical data of relevant literature are shown in attached tables 1-4, experimental data of longitudinal inclination angle ranging from-20 degrees to 20 degrees and transverse inclination angle ranging from 0 degrees to 20 degrees are summarized, each 10 degrees is provided with a node, the node value is a relevant experimental mean value, and for the non-obtained inclination angle, the reduction coefficient is determined by linear interpolation between two known points, and statistical values of the nodes are shown in tables 1 and 2:
table 1 corresponding roll reduction coefficient
Type/heel | 0° | 10° | 20° |
Stairs/walks upwards | 1 | 0.9174 | 1.1174 |
Stairs/walk downwards | 1 | 0.8816 | 0.8006 |
Corridor | 1 | 0.9639 | 0.8710 |
Table 2 trim corresponding reduction factor
Type/pitch | -20° | -10° | 0° | 10° | 20° |
Corridor | 0.9111 | 0.9807 | 1 | 0.8969 | 0.7895 |
Stairs/walks upwards | 1.0655 | 1.0801 | 1 | 0.8536 | 0.6598 |
Stairs/walk downwards | 0.7751 | 0.8564 | 1 | 0.9901 | 0.817 |
Reduction factor k 'in the case of transverse tilting' 2 The calculation formula is, wherein, ψ is the transverse inclination:
Reduction factor k in the case of longitudinal tilting 2 The "calculation formula" is defined as the formula,the pitch angle is:
Thus, the person specific flow rate F S Obtained by interpolation method, the calculation formula is:
wherein: n is the total number of people (p) in the path; l is the path length (m); w (W) C Is the static width (m) of the path.
The method for calculating the personnel evacuation travel time T in the step 2.3 specifically comprises the following steps:
T=γ·t l
wherein: gamma is a correction coefficient, and 1.3 is taken; l (L) j In the path scheme, i is an evacuated people stream, and j branch paths are included in the evacuation route; l (L) y For the length of the y-th branch path in the evacuation route; v y For the speed of the person in the y-th branch path in the evacuation route; n (N) y For the total number of people in the y-th branch path in the evacuation route; f (F) sy A specific flow in the y-th branch path in the evacuation route; w (W) cy For the static width in the y-th branch path in the evacuation route.
Static width (W) C ): the width measured from the handrail of the corridor and stairway, and the width (m) of the passageway with the door in an open condition.
Path length (L): when people escape in emergency, the distance length (m) of the personnel passing through the corridor, the stairway and other facilities.
A method for calculating the emergency escape and evacuation time of personnel under the condition of ship water inflow is provided, wherein the specific formula for calculating the evacuation time in the step 2.4 is as follows:
wherein: r is reaction time, the influence of ship movement caused by water inflow is considered, the ship body can shake violently in the initial stage of water inflow, the time from the beginning of water inflow to the ending of the violent shaking comprises the alarm time of a water inflow alarm and the reaction time after personnel receive alarm signals, so that the stage can be taken as evacuation preparation time under the condition of water inflow, the research on the water inflow process of ship breakage finds that the stage corresponds to the instantaneous water inflow stage in the water inflow process of ship breakage, the R value method can be judged through a ship rolling curve, the specific judging method refers to the figure 1, and the rolling curve is obtained by a part I; t is the person evacuation travel time.
Table 3 shows the reduction coefficient of the transverse inclined corridor
Table 4 shows the reduction coefficient of the laterally inclined stairs
Table 5 shows the reduction coefficient of the longitudinally inclined corridor
Table 6 shows the reduction coefficient of the longitudinally inclined stairs
Claims (8)
1. The method for calculating the emergency escape and evacuation time of the personnel under the condition of water inflow of the ship is characterized by comprising the following steps of:
step 1: based on CFD, simulating water inflow of the ship, acquiring a ship body movement angle, and providing data for subsequent calculation, wherein the ship body movement angle is roll and pitch;
step 2: calculating the emergency escape and evacuation time of the personnel under the condition of water inflow based on the data in the step 1;
the step 2 specifically comprises the following steps:
step 2.1: planning the shortest escape route of personnel, and obtaining the number N of the branch routes of each escape route and the inclination angle of the ship body;
step 2.2: calculating the walking speed v and the specific flow F of the personnel on the evacuation route S ;
Step 2.3: calculating personnel evacuation travel time T;
step 2.4: calculating the time REST required for personnel evacuation under the condition of water inflow, wherein rest=r+t, and R is the reaction time;
step 2.2 calculates the walking speed v and the specific flow F of the personnel on the evacuation route S The method comprises the following steps:
personnel walking speed v:
v=v 0 ·k 1 ·k 2 '·k 2 "
wherein: v 0 Indicating the initial speed, k of the person 1 Representing the influence of personnel density on walking speed, the MSC gives personnel density D asIn relation to the speed of travel of the person, by reducing the coefficient k 1 The form represents:
wherein: n is the total number of people in the path; l is the path length; w (W) C Is the static width of the path;
reduction factor k 'in the case of transverse tilting' 2 The calculation formula is, wherein, ψ is the transverse inclination:
Reduction factor k in the case of longitudinal tilting 2 The calculation formula is as follows, wherein phi is the pitch angle:
Thus, the person specific flow rate F S Obtained by interpolation method, the calculation formula is:
wherein: n is the total number of people in the path; l is the path length; w (W) C Is the static width of the path.
2. The method for calculating the emergency escape and evacuation time of personnel under the condition of water intake of a ship according to claim 1, wherein the step 1 specifically comprises the following steps:
step 1.1: establishing a CFD numerical model; the CFD numerical model only considers a mass conservation equation and a momentum conservation equation;
step 1.2: solving the CFD numerical model in the step 1.1, dividing a solution domain of the fluid into a limited number of mutually adjacent control bodies, and applying a mass conservation equation and a momentum conservation equation to each control body;
step 1.3: capturing the interface of water and gas by using a VOF method;
step 1.4: the six-degree-of-freedom solver and the dynamic grid technology are utilized to process the ship movement;
step 1.5: monitoring a ship movement angle and obtaining a ship movement curve; setting calculation time and step length, and performing iterative calculation.
3. The method for calculating the emergency escape and evacuation time of personnel under the condition of water inflow of a ship according to claim 2, wherein the mass conservation equation and the momentum conservation equation in the step 1.1 are specifically,
the mass conservation equation is:
the conservation of momentum equation is:
wherein: u, v, w are velocity components of velocity vector v along x, y, z; x, Y, Z is the unit mass force applied to the microcell in the x, y, z directions; p fluid pressure; a coefficient of dynamic viscosity;
the step 1.2 specifically comprises the following conservation equations:
4. The method for calculating the emergency escape and evacuation time of personnel in the case of water intake of a ship according to claim 2, wherein the step 1.3 is specifically that the ratio F of the fluid volume in each calculation grid to the grid volume; when F is 1, the fluid in the grid is the main phase fluid; when F is 0, the grid is filled with secondary phase fluid; when 0 < F < 1: representing a multiphase flow boundary within the mesh; the interface of the water and air is denoted as,
r w +r a =0
wherein: r is (r) w Representing the volume fraction of water; r is (r) a Representing the volume fraction of air; v is the control volume.
5. The method for calculating the emergency escape and evacuation time of personnel under the condition of water intake of a ship according to claim 2, wherein the step 1.4 is specifically to calculate the displacement and the rotation angle of the center of gravity, and use the stress and the moment of an object; the displacement control equation for the center of gravity is as follows:
wherein: m represents a rigid body mass;indicating the gravity center stress; />Representing the center of gravity displacement acceleration;
angular acceleration:
wherein: l represents an inertial tensor;representing a moment tensor; />Representing the rigid body angular velocity; />The method is the conversion from an inertial coordinate system to a satellite coordinate system; after the angular acceleration and the displacement acceleration are calculated, the dynamic grid calculation can update the rigid body position using the angular velocity and the displacement velocity.
6. The method for calculating the emergency escape and evacuation time of personnel under the condition of water intake of a ship according to claim 1, wherein the step 2.1 is to plan the shortest escape path of personnel, plan the shortest escape path by using Dijstra algorithm, and obtain the number N of people on the branch path.
7. The method for calculating the emergency escape and evacuation time of personnel under the condition of water intake of a ship according to claim 1, wherein the step 2.3 of calculating the personnel evacuation travel time tspeculiarly comprises the following steps:
T=γ·t l
wherein: gamma is a correction coefficient; l (L) j In the path scheme, i is an evacuated people stream, and j branch paths are included in the evacuation route; l (L) y For the length of the y-th branch path in the evacuation route; v y For the speed of the person in the y-th branch path in the evacuation route; n (N) y For the total number of people in the y-th branch path in the evacuation route; f (F) sy A specific flow in the y-th branch path in the evacuation route; w (W) cy For the static width in the y-th branch path in the evacuation route.
8. The method for calculating the emergency escape and evacuation time of personnel under the condition of water intake of a ship according to claim 1, wherein the specific formula for calculating the evacuation time in the step 2.4 is as follows:
wherein: r is the reaction time; t is the person evacuation travel time.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111288697.3A CN114943630B (en) | 2021-11-02 | 2021-11-02 | Method for calculating emergency escape and evacuation time of personnel under ship water inflow condition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111288697.3A CN114943630B (en) | 2021-11-02 | 2021-11-02 | Method for calculating emergency escape and evacuation time of personnel under ship water inflow condition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114943630A CN114943630A (en) | 2022-08-26 |
CN114943630B true CN114943630B (en) | 2023-06-13 |
Family
ID=82906173
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111288697.3A Active CN114943630B (en) | 2021-11-02 | 2021-11-02 | Method for calculating emergency escape and evacuation time of personnel under ship water inflow condition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114943630B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110059336A (en) * | 2019-02-19 | 2019-07-26 | 中国舰船研究设计中心 | The appraisal procedure of ship personnel evacuation |
CN110516331A (en) * | 2019-08-15 | 2019-11-29 | 哈尔滨工程大学 | Marine personnel's emergency escape departure time calculation method under a kind of fire condition |
CN110533232A (en) * | 2019-08-15 | 2019-12-03 | 哈尔滨工程大学 | A method of calculating the most short evacuation path of passenger boat personnel |
CN112836878A (en) * | 2021-02-04 | 2021-05-25 | 厦门市建筑科学研究院有限公司 | BIM technology-based emergency evacuation method and device |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10192411B2 (en) * | 2012-12-13 | 2019-01-29 | Oneevent Technologies, Inc. | Sensor-based monitoring system |
GB2554075B (en) * | 2016-09-15 | 2021-05-19 | Equinor Energy As | Optimising fire protection for an offshore platform |
CN110503257B (en) * | 2019-08-15 | 2021-09-03 | 哈尔滨工程大学 | MSC-based marine personnel emergency escape path planning method |
-
2021
- 2021-11-02 CN CN202111288697.3A patent/CN114943630B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110059336A (en) * | 2019-02-19 | 2019-07-26 | 中国舰船研究设计中心 | The appraisal procedure of ship personnel evacuation |
CN110516331A (en) * | 2019-08-15 | 2019-11-29 | 哈尔滨工程大学 | Marine personnel's emergency escape departure time calculation method under a kind of fire condition |
CN110533232A (en) * | 2019-08-15 | 2019-12-03 | 哈尔滨工程大学 | A method of calculating the most short evacuation path of passenger boat personnel |
CN112836878A (en) * | 2021-02-04 | 2021-05-25 | 厦门市建筑科学研究院有限公司 | BIM technology-based emergency evacuation method and device |
Also Published As
Publication number | Publication date |
---|---|
CN114943630A (en) | 2022-08-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tezduyar et al. | Modelling of fluid–structure interactions with the space–time finite elements: solution techniques | |
Muradoglu et al. | A front-tracking method for computation of interfacial flows with soluble surfactants | |
CN112100835A (en) | High-efficiency high-precision numerical simulation method suitable for complex flow | |
CN114943630B (en) | Method for calculating emergency escape and evacuation time of personnel under ship water inflow condition | |
Bagué et al. | Dynamic stability analysis of a hydrofoiling sailing boat using CFD | |
Snyder et al. | Comparison of experimental and computational ship air wakes for YP class patrol craft | |
CN116415463A (en) | Ship simulation state data analysis method | |
Dooley et al. | Effects of waves, motions and atmospheric turbulence on ship airwakes | |
Xu et al. | Experimental study on inertial hydrodynamic behaviors of a complex remotely operated vehicle | |
O’Sullivan et al. | Air-flow distortion and wave interactions on research vessels: An experimental and numerical comparison | |
Lu et al. | A preliminary investigation of airflow field in designated refuge floor | |
TABIB et al. | TOWARDS UNDERSTANDING WIND IMPACT FOR DRONE OPERATIONS: A COMPARISON OF WIND MODELS OPERATING ON DIFFERENT SCALES IN A NESTED MULTISCALE SET-UP. | |
Söding | Flow computations for ship safety problems | |
Wu et al. | A new approach to aircraft ditching analysis by coupling free surface lattice Boltzmann and immersed boundary method incorporating surface tension effects | |
Kelly | The development, validation, and integration of aircraft carrier airwakes for piloted flight simulation | |
CN114943191B (en) | Personnel emergency escape path planning method under ship water inflow condition and computer equipment | |
Dietz et al. | Application of vorticity confinement to the prediction of flows over complex bodies-A survey of recent results | |
Polsky et al. | A computational study of outwash for a helicopter operating near a vertical face with comparison to experimental data | |
Valenko et al. | Dynamic model of scuba diver buoyancy | |
JPH11124078A (en) | Method and device for simulating bubble covering ship body and storage medium containing bubble covering simulation program readable by computer | |
Chen et al. | Modelling and simulation for underwater hovering control based on ballast tank | |
Jansson et al. | Adaptive simulation of unsteady flow past the submerged part of a floating wind turbine platform | |
CN109063282A (en) | A kind of foundation pit layout method constrained based on largest passages and density of personnel | |
CN117610462B (en) | Sailing boat operation simulation method, sailing boat operation simulation system, sailing boat operation simulation application and electronic equipment | |
CN114740894B (en) | Aircraft guidance method and system based on attention mechanism and gated cycle unit |
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 |