Applicability of Calculation Formulae of Impact Force by Tsunami Driftage

<jats:p>The aftermath of the Indian Ocean tsunami on 26 December 2004 triggered by the off Sumatra earthquake (magnitude “M” = 9.1), and the Great East Japan earthquake of 11 March 2011 off the Pacific coast of Tohoku (M = 9.0), evidence the secondary damage from driftage collision due to large tsunami waves. To prevent this type of damage, the establishment of methods for predicting driftage movement and calculating the impact force by driftage is necessary. Several numerical models have been developed to predict the driftage movement of objects. Every year, these improve in accuracy and usability. In contrast, there are many calculation formulae for calculating the impact force. However, since there are considerable differences between values calculated using these formulae, the reliability of each formula is unknown. Therefore, in this research, one team of the committee on tsunami research of the Japan Society of Civil Engineers summarizes the main calculation formulae of impact forces that have been proposed until 2019. In addition, for each type of driftage (driftwood, containers, cars, ships), we compare calculation values of these formulae with measured data of large-scale experiments. Finally, we check the range of calculation values for each formula up to 15 m/s in collision velocity and clarify then the following facts: (1) In the case of driftwood, the formulae of Matsutomi, Federal Emergency Management Agency (FEMA) and National Oceanic and Atmospheric Administration (NOAA), and American Society of Civil Engineers (ASCE) are most reliable; (2) In the case of containers, the formulae of Matsutomi, Arikawa et al., FEMA and NOAA, Ikeno et al., and ASCE are most reliable; (3) In the case of cars, the formulae of FEMA and NOAA, and ASCE are most reliable; (4) In the case of ships, the formulae of Mizutani, FEMA and NOAA, and ASCE are most reliable.</jats:p>

section_content: When a large tsunami caused by an earthquake reaches coastal communities, secondary damage by driftage and collision of objects is observed. For instance, when the damage surveys in Thailand (January 2005) and in Sri Lanka (May 2005) took place, following the off Sumatra earthquake (magnitude "M" = 9. 1) and the Indian Ocean tsunami, this type of damage was confirmed. Later, damage surveys conducted in Japan (March and April 2011) in the aftermath of the Great East Japan earthquake (M = 9. 0) evidenced, once again, the significant damage caused by driftage and collision of objects. Some examples are shown in the following. 

(1) Driftwood Figure 1a shows fallen trees that collided against a hotel at Khao Lak in Thailand, while Figure 1b portraits some logs that pierced the windows of a building at Ofunato in Japan. 

Similarly, in Figure 2b, a container broke the wall of a large warehouse at the Soma port in Japan. 

(3) Cars Figure 3a shows a dump truck that collided against a warehouse at the Kirinda fishing port in Sri Lanka. In Figure 3b, a passenger car collided against a shop at Miyako in Japan. Figure 3c shows how cars were pushed against the corner of a school building at Sendai in Japan. 

(4) Ships In Figure 4a, a patrol boat stopped at the foot of a mountain 1 km inland of Khao Lak in Thailand. In addition, Figure 4b shows a building with a big hole made by the patrol boat in Khao Lak. On the other hand, Figure 4c shows a large fishing boat left behind in the ruins near a fishing port at Kesennuma in Japan. Figure 4d shows a broken steel frame building on the shortest line between the large fishing boat and the fishing port. (2) Containers In Figure 2a, a container collided against a house near the coast at Hitachi in Japan. Similarly, in Figure 2b, a container broke the wall of a large warehouse at the Soma port in Japan. 

(1) Driftwood Figure 1a shows fallen trees that collided against a hotel at Khao Lak in Thailand, while Figure 1b portraits some logs that pierced the windows of a building at Ofunato in Japan. 

(2) Containers In Figure 2a, a container collided against a house near the coast at Hitachi in Japan. Similarly, in Figure 2b, a container broke the wall of a large warehouse at the Soma port in Japan. 

(3) Cars Figure 3a shows a dump truck that collided against a warehouse at the Kirinda fishing port in Sri Lanka. In Figure 3b, a passenger car collided against a shop at Miyako in Japan. Figure 3c shows how cars were pushed against the corner of a school building at Sendai in Japan. 

(4) Ships In Figure 4a, a patrol boat stopped at the foot of a mountain 1 km inland of Khao Lak in Thailand. In addition, Figure 4b shows a building with a big hole made by the patrol boat in Khao Lak. On the other hand, Figure 4c shows a large fishing boat left behind in the ruins near a fishing port at Kesennuma in Japan. Figure 4d shows a broken steel frame building on the shortest line between the large fishing boat and the fishing port. (3) Cars Figure 3a shows a dump truck that collided against a warehouse at the Kirinda fishing port in Sri Lanka. In Figure 3b, a passenger car collided against a shop at Miyako in Japan. Figure 3c shows how cars were pushed against the corner of a school building at Sendai in Japan. (4) Ships In Figure 4a, a patrol boat stopped at the foot of a mountain 1 km inland of Khao Lak in Thailand. In addition, Figure 4b shows a building with a big hole made by the patrol boat in Khao Lak. On the other hand, Figure 4c shows a large fishing boat left behind in the ruins near a fishing port at Kesennuma in Japan. Figure 4d shows a broken steel frame building on the shortest line between the large fishing boat and the fishing port. To prevent damage by driftage, it is important to establish a method for predicting driftage movement and calculating the impact force by driftage. Nistor et al. [1] reviewed existing papers published until 2016 about the transportation and the impact force of debris caused by a tsunami. 

Experiential methods to predict driftage movement can be easily handled as the range at which driftage can reach is estimated based on simple assumptions, as demonstrated by Naito et al. [2]. In contrast, driftage movement can be evaluated through numerical simulations (Nistor et al. [1] ), though the further improvement of the simulation accuracy and the convenience of the handling for engineers is desired. Regarding evaluation methods by numerical models, Yoneyama et al. [3] conducted a study by which the 3-dimensional VOF method was used for fluid motion calculation, and Murase et al. [4] used the storm surge and tsunami simulator in oceans and coastal areas (STOC) model. 

For methods to calculate the impact force, although there are numerical models, such as the elastoplastic finite deformation analysis of Magoshi et al. [5], they have not reached the level at which they can be used simply and reliable for engineers. Thus, traditional methods generally used are semi-theoretical formulae and empirical formulae, as shown in Tables 1 2 3. However, since there are considerable differences between values calculated by these formulae, the reliability of each formula is unknown. In this regard, Kaida and Kihara [6] investigated calculation formulae of impact forces proposed until 2012, and they checked the performance range of calculated values by each formula up to a maximum value of collision velocity of 10 m/s. Finally, they summarized the calculation formulae that are effective for each type of driftage. Moreover, they compiled the appropriate values of shaft rigidity, which is an important parameter in some formulas presented in Tables 1 2 3. Furthermore, Stolle et al. [7, 8] indicated that the calculation accuracy of the impact force of low rigid debris (like driftwood, containers) could be improved by estimating the rigidity of the debris exactly.
section_name: Introduction

section_content: Driftwood

1+ εo r i (1+µ ro εo ) F i is the impact force (N); u max is the maximum velocity of driftwood; k is the effective contact stiffness; m d is the mass of driftwood; c M is the added mass coefficient; m f is the mass of displaced fluid;

β is the decreasing coefficient by the orientation angle ϕ; e is the decreasing coefficient by eccentricity; ε o is the length from the barycenter to the collided point; r i is the turning radius; µ is the friction coefficient between the driftwood and the collided body; r o is the radius of the sectional area of the driftwood. Haehnel and Daly (2004) [9] Driftwood and Containers

u is the collision velocity; c I is the important factor; c o is the orientation coefficient; c D is the coefficient of water depth; c B is the barrier coefficient; R max is the maximum response ratio to the collision load; ∆t is the continuing time of the collision; k is the effective stiffness.
section_name: Driftage Formula Reference

section_content: Driftwood (the Mass is about 450 kg) and Containers (the Empty Mass is 2200-3800 kg)

F i is the impact force (N); 1. 3 is the importance coefficient of risk category IV; u max is the maximum fluid velocity; m d is the mass of driftage; k is the effective stiffness (2. 4 × 10 6 N/m to driftwood; 6. 5 × 10 8 N/m to 12. 2 m empty containers; 1. 5 × 10 9 N/m to 6. 1 m empty containers); c is the hydrodynamic mass coefficient (0. 0 to the long axis direction for the collision of driftwood; 0. 2 to the long axis direction for the collision of 12. 2 m containers; 1. 0 to the transverse direction for the collision of 12. 2 m containers; 0. 3 to the long axis direction of 6. 1 m containers; 1. 0 to the transverse direction of 6. 1 m containers); g is the gravitational acceleration; R is the maximum run-up height × 1. 3; z is the ground height on the water surface. 

FEMA & NOAA (2012) [11] Table 2. Impact force formula (no. 2).
section_name: ASCE (2010) [10]

section_content: Driftwood, Containers, etc. 

F i is the impact force (N), u is the collision velocity; k d and m d are the stiffness and the mass of driftage; E is the elastic modulus of driftage; A d and L d are the sectional area and the length of driftage. 

The impulse duration of the steel tube can be obtained using

Aghl, Naito, and Riggs (2014) [12] Driftwood, Containers (6. 1-12. 2 m), Cars, Ships, etc. 

R max is the dynamic response ratio to the collision load (it is in the range of 0. 0-1. 8 according to the collision action time); I TSU is the importance factor (1. 0-1. 25); C o is the orientation coefficient (=0. 65); u max is the maximum fluid velocity; k is the smaller one of (i) the stiffness of the driftage; (ii) the transverse direction stiffness of the collided body; m d is the mass of the driftage. The stiffness of the driftwood can be calculated by using k = EA/L. E is the elasticity modulus of the long axis direction of driftage; A and L are the sectional area and the length. 

The impulse duration for empty containers t d can be calculated by using

The impulse duration for loaded containers t d can be calculated by using

Here, m contents is 50% of the maximum capacity of the container.
section_name: Driftage Formula Reference

section_content: Ships

F is the impact force (kN); W is the weight of the ship (kN); V is the collision velocity; l d is the stopping distance.
section_name: ASCE (2015) [13]

section_content: Fi is the impact force (N); Rmax is the dynamic response ratio to the collision load (it is in the range of 0. 0-1. 8 according to the collision action time); ITSU is the importance factor (1. 0-1. 25); Co is the orientation coefficient (=0. 65); umax is the maximum fluid velocity;

k is the smaller one of (i) the stiffness of the driftage;

(ii) the transverse direction stiffness of the collided body; md is the mass of the driftage. The stiffness of the driftwood can be calculated by using = /. 

E is the elasticity modulus of the long axis direction of driftage; A and L are the sectional area and the length. The impulse duration for empty containers td can be calculated by using = 2

The impulse duration for loaded containers td can be calculated by using =

Here, mcontents is 50% of the maximum capacity of the container. ASCE (2015) [13] Ships = 4 is the impact force (kN); is the weight of the ship (kN);

is the collision velocity; is the stopping distance. F is the impact force (kN); is the apparent mass coefficient (1. 7 for a bore or a surge, 1. 9 for a steady flow);

V is the collision velocity (m/s); D and L are the diameter and the length of driftwood (m);

is the yield stress of driftwood (kN/m 2 ]; is the unit weight of driftwood (kN/m 3 ); is the gravitational acceleration. 

Matsutomi (1999) [15] Bridge and Engineering Association (1978) [14] Driftwood (Long Axis Direction Collision)

F is the impact force (kN); C MA is the apparent mass coefficient (1. 7 for a bore or a surge, 1. 9 for a steady flow); V is the collision velocity (m/s); D and L are the diameter and the length of driftwood (m); σ f is the yield stress of driftwood (kN/m 2 ]; γ is the unit weight of driftwood (kN/m 3 ); g is the gravitational acceleration. 

Matsutomi (1999) [15] Table 3. Impact force formula (no. 3).
section_name: 𝐹 = 𝑅 𝐼 𝐶 𝑢 𝑘𝑚

section_content: Driftwood, Cars, Ships, etc. 

v is the surface fluid velocity (m/s). 

Japan Road Association (2002) [16] Driftwood

F is the impact force (kN); S is a coefficient (5. 0); C MA is the added mass coefficient [pillar oriented transversely: 2. 0 (two dimension), 1. 5 (three dimension); square pillar oriented transversely: 2. 0-4. 0 (two dimension); 1. 5 (three dimension); pillar oriented longitudinally: about 2. 0, sphere: about 0. 8]; V is the bore fluid velocity (m/s); D is the collision section size (m); L is the length of driftwood (m); M is the mass of driftwood (kg); g is the gravitational acceleration (m/s 2 ). 

Ikeno and Tanaka (2003) [17] Table 3. Cont.
section_name: Driftage Formula Reference

section_content: Containers

F is the impact force (kN); ρ w is the water density (t/m 3 ); η m is the maximum run-up height; B c is the width of the container; V X is the drift velocity (m/s); W is the container weight; g is the gravitational acceleration;

dt is the collision time.
section_name: Driftage Formula Reference

section_content: Small Ships F = 2 WVx gdt F is the impact force (kN); W is the weight of the small ship (kN); V X is the drifting velocity of the small ship (m/s); g is the gravitational acceleration; dt is the collision time. 

Mizutaniet al. ( 2007 ) [19] Driftwood, Containers, etc. 

F is the impact force (kN); a is 1/2 of the radius of the collision sectional area (m); E is the elastic modulus (kN/m 2 ]; ν i is Poisson's ratio; m i is the mass (kg); V is the collision velocity (m/s); γ p is the energy attenuation effect by plasticity (0. 25); the suffix i of m and k indicates the driftage and the collided body.
section_name: Mizutani et al. (2005) [18]

section_content: (2010) [21] Driftwood and Containers

4 F is the impact force at the time of the elastic collision (kN); k c is a constant 0. 243; C MA is the apparent mass coefficient (2. 0 for wooden pillar oriented longitudinally, 1. 0 for containers); M is the mass of driftage (kg); V is the collision velocity (m/s); D is the collision sectional area size (m); E is the representative stiffness =

Ikeno et al. ( 2013 ) [22] In this research, we summarize the main calculation formulae of impact forces proposed by 2019, and the calculation values by these formulae are compared with measured data of large-scale experiments for each type of driftage (driftwood, shipping containers, cars, ships). Moreover, the range of calculated values by each formula up to a collision velocity of 15 m/s is investigated. Finally, the reliability of each formula is discussed and clarified.
section_name: Arikawa et al. (2007) [20],

section_name: Accuracy of Formulae for Calculating the Impact Force

section_content: The main calculation formulae of impact forces by tsunami driftage proposed by 2019 are summarized in Tables 1 2 3. Although these calculation formulae were designed to reflect the proposer's experimental data, there are considerable differences among the proposed formulae. A possible reason for these differences might be that several formulae cannot make the difference in rigidity between the driftage and the collided body correctly into consideration. When the driftage collides with the collided body, the impact force is reduced by carrying out an elastic modification. Therefore, the impact force is different if the rigidity of the driftage or the collided body is not identical, even when the velocity at the time of collision is the same. In other words, if the differences of bodies with different rigidity from the experiment used to verify the proposed formula are not considered, the accuracy of the results is reduced. 

In the case of the tsunami force in which seawater collides against a rigid body made of concrete or steel, the seawater deforms easily compared with the collided body. Moreover, the inverse of whole rigidity is proportional to the inverse of each rigidity of the collision body and the collided body. Therefore, we can say that the rigidity of the seawater is very low, and the decrease of the tsunami force is governed not by the collided body but by seawater. Therefore, we do not need to make the difference in rigidity of the collided body into consideration.
section_name: List of Main Calculation Formulae

section_content: We searched for large-scale experimental data with the size, weight, rigidity of the driftage and the collided body, the collision velocity, and the impact force for each driftage type (driftage, containers, cars, ships) in the published literature. We compare the calculated values of the impact force and measurements. 

Here, as the improved formula of Haehnel and Daly [9] is the formula of FEMA and NOAA [11] and the improved formula of ASCE [10] and Aghl, Naito and Riggs [12] is the formula of ASCE [13], the formulae of Haehnel and Daly [9], ASCE [10], and Aghl, Naito and Riggs [12] are not examined. 

Because we could not find materials, which can determine the stop distance of driftage except for ships, the formula of the Bridge and Engineering Association [14] was examined only for the case of ships. Moreover, since we could not find materials, which can determine the yield stress of ships, the formula of Matsutomi [15] was not examined for the case of ships. The formula of Mizutani et al. [18] was examined for containers and cars, which can estimate a maximum run-up height easily because they do not float and move immediately. In addition, the formula of [19] was examined in the case of ships. Moreover, for the formula of ASCE [13], the values for R max, I TSU, C o were set as 1. 8, 1. 25, and 0. 65, respectively. Although ASCE [13] recommends some values for the effective stiffness k, when this was not available, the recommended values by FEMA and NOAA [11] were used.
section_name: Correlation Examination of the Main Calculation Formulae

section_content: Although there are many small-scale experiments, we selected the large-scale experiments of Matsutomi [15] and Ikeno et al. (Central Research Institute of Electric Power Industry) [23] because they specify necessary information, and we actually used the data for collision velocities larger than 2. 0 m/s in their experiments, which were measured in unsteady flow using a water flume. 

The calculation conditions are as follows:

(1) Experiment 1 of Matsutomi [15] ( (c) Coefficients in the formula Each apparent mass coefficient C MA was set up according to each proposer's recommended value (1. 7 for Matsutomi [15], 2. 0 for Ikeno and Tanaka [17], 1. 7 for Arikawa et al. [20], 2. 0 for Ikeno et al. [22] ). The constant 0. 0 was set up to the coefficient C of FEMA and NOAA [11]. The axial stiffness k 1 of the driftwood was 269 MN/m (since the axial stiffness changes according to the dimensions of the object, the recommended value of FEMA and NOAA [11] is usable only to specified dimensions. Therefore, this value was calculated using the elasticity modulus of driftwood), and the bending stiffness k 2 of the steel was 5 MN/m.
section_name: Driftwood

section_content: The measured impact force and the collision velocity are 0. 022 MN and 2. 0 m/s, respectively. 

(2) Experiment 2 of Matsutomi [15] (a) Parameters of the driftwood (larch) 0. 3 m in diameter, 2. 7 m long, 101. 15 kg, 25 MN/m 2 yield stress, 0. 4 Poisson's ratio, 8000 MN/m 2 elasticity modulus. 

(b) Parameters of the collided body (steel) and the apparent mass coefficients used the same values as (1) Matsutomi's experiment 1 [15]. (c) Coefficients in the formula;

The coefficient C of FEMA and NOAA [11] was chosen as 0. 0. The axial stiffness k 1 of the driftwood was 209 MN/m, and the bending stiffness k 2 of the steel was 5 MN/m.
section_name: (d) Measured impact force and the collision velocity

section_content: The measured impact force and the collision velocity are 0. 026 MN and 2. 0 m/s, respectively. 

(3) Matsutomi's experiment 3 [15] The coefficient C of FEMA and NOAA [11] was 0. 0. The axial stiffness k 1 of the driftwood was 141 MN/m, and the bending stiffness k 2 of the steel was 5 MN/m.
section_name: (d) Measured impact force and the collision velocity

section_content: The measured impact force and the collision velocity are 0. 036 MN and 2. 0 m/s, respectively. The coefficient C of FEMA and NOAA [11] was 0. 0. The axial stiffness k 1 of the driftwood was 776 MN/m, and the bending stiffness k 2 of the steel was 12. 8 MN/m. (d) Measured impact force and the collision velocity;

The measured impact force and the collision velocities are 0. 080 MN and 2. 0 m/s, 0. 082 MN and 4. 0 m/s, respectively. The coefficient C of FEMA and NOAA [11] was 0. 0. The axial stiffness k 1 of the driftwood was 1040 MN/m, and the bending stiffness k 2 of the steel was 12. 8 MN/m.
section_name: (d) Measured impact force and the collision velocity

section_content: The measured impact force and the collision velocity are 0. 054 MN and 4. 1 m/s, respectively. The coefficient C of FEMA and NOAA [11] was 0. 0. The axial stiffness k 1 of the driftwood was 521 MN/m, and the bending stiffness k 2 of the steel was 12. 8 MN/m.
section_name: (d) Measured impact force and the collision velocity

section_content: The measured impact force and the collision velocity are 0. 077 MN and 4. 1 m/s, respectively. The verification results are shown in Figure 5. 

The Japan Road Association's formula [16] overestimates, and Ikeno and Tanaka's formula [17] underestimates the measured data. The formulae of FEMA and NOAA [11], ASCE [13], and Matsutomi [15] result in values reasonably close to the measured values. Since the formulae of Arikawa et al. [20] and Ikeno et al. [22] correlate well with measured values, provided the formulae are improved so that parameters concerning the rigidity of the driftage and the collided body can be decided adequately, the calculation accuracy of the improved formulae will improve.
section_name: (d) Measured impact force and the collision velocity

section_content: Large-scale experiments with the necessary information were conducted by Aghl et al. [12] and Arikawa et al. [20] (Port and Airport Research Institute). Only data with collision velocities larger than 1. 4 m/s were used. The experiment of Aghl et al. [12] was implemented in the air by using a full-scale container, and the experiment of Arikawa et al. [20] was conducted in unsteady flow using a water flume and a container model of 1:5 in scale. Moreover, when using the formula of Mizutani et al. [18], 1000 kg/m 3 was used for ρw since freshwater was used in their experiment, and half of the height of their containers was used for ηm so that 2 × × becomes a collision sectional area. Furthermore, the impulse duration of ASCE [13] was used for dt. 

The calculation conditions are as follows: (1) Experiment of Arikawa et al. [20] ( [22] ). The coefficient C of FEMA and NOAA [11] was 0. 2. The axial stiffness k1 of the small container was 1080 MN/m (since the axial stiffness changes according to the dimensions of the object, the recommended value of FEMA and NOAA [11] cannot be used. Therefore, this value was calculated using the elasticity modulus of the small container), and the bending stiffness k2 of the concrete was 2220 MN/m. The Japan Road Association's formula [16] overestimates, and Ikeno and Tanaka's formula [17] underestimates the measured data. The formulae of FEMA and NOAA [11], ASCE [13], and Matsutomi [15] result in values reasonably close to the measured values. Since the formulae of Arikawa et al. [20] and Ikeno et al. [22] correlate well with measured values, provided the formulae are improved so that parameters concerning the rigidity of the driftage and the collided body can be decided adequately, the calculation accuracy of the improved formulae will improve.
section_name: Containers

section_content: Large-scale experiments with the necessary information were conducted by Aghl et al. [12] and Arikawa et al. [20] (Port and Airport Research Institute). Only data with collision velocities larger than 1. 4 m/s were used. The experiment of Aghl et al. [12] was implemented in the air by using a full-scale container, and the experiment of Arikawa et al. [20] was conducted in unsteady flow using a water flume and a container model of 1:5 in scale. Moreover, when using the formula of Mizutani et al. [18], 1000 kg/m 3 was used for ρ w since freshwater was used in their experiment, and half of the height of their containers was used for η m so that 2 × η m × B c becomes a collision sectional area. Furthermore, the impulse duration of ASCE [13] was used for dt. 

The calculation conditions are as follows:

( (c) Coefficients in the formula Each apparent mass coefficient C MA was set up according to each proposer's recommendation value (1. 7 for Matsutomi [15], 4. 0 for Ikeno and Tanaka [17], 1. 0 for Ikeno et al. [22] ). The coefficient C of FEMA and NOAA [11] was 0. 2. The axial stiffness k 1 of the small container was 1080 MN/m (since the axial stiffness changes according to the dimensions of the object, the recommended value of FEMA and NOAA [11] cannot be used. Therefore, this value was calculated using the elasticity modulus of the small container), and the bending stiffness k 2 of the concrete was 2220 MN/m.
section_name: Containers

section_content: The measured impact force and the collision velocities are 0. 033 MN and 1. 70 m/s, 0. 049 MN and 2. 15 m/s, 0. 057 MN and 2. 40 m/s, respectively. 

(2) Experiment of Aghl et al. [12] (a) Parameters of the 6. 1 m container (steel) (c) Coefficients in the formula Each apparent mass coefficient C MA was set up as the ones in the air (1. 0 for Matsutomi [15], 2. 0 for Ikeno and Tanaka [17], 1. 0 for Ikeno et al. [22] ). The constant 0. 0 and 85 MN/m, which FEMA and NOAA [11] recommends, were set up for the coefficient C and the effective stiffness k of FEMA and NOAA [11], respectively. Moreover, 42. 9 MN/m recommended by ASCE [13] was set up for the effective stiffness k of ASCE [13].
section_name: (d) Measured impact force and the collision velocity

section_content: The measured impact force and the collision velocities are 0. 795 MN and 1. 43 m/s, 0. 995 MN and 1. 80 m/s, 1. 195 MN and 2. 13 m/s, 1. 315 MN and 2. 42 m/s, respectively. 

The verification results are shown in Figure 6. (c) Coefficients in the formula Each apparent mass coefficient CMA was set up as the ones in the air (1. 0 for Matsutomi [15], 2. 0 for Ikeno and Tanaka [17], 1. 0 for Ikeno et al. [22] ). The constant 0. 0 and 85 MN/m, which FEMA and NOAA [11] recommends, were set up for the coefficient C and the effective stiffness k of FEMA and NOAA [11], respectively. Moreover, 42. 9 MN/m recommended by ASCE [13] was set up for the effective stiffness k of ASCE [13]. 

(d) Measured impact force and the collision velocity The measured impact force and the collision velocities are 0. 795 MN and 1. 43 m/s, 0. 995 MN and 1. 80 m/s, 1. 195 MN and 2. 13 m/s, 1. 315 MN and 2. 42 m/s, respectively. 

The verification results are shown in Figure 6. 

The formula of the Japan Road Association [16] overestimates, and the formula of Ikeno and Tanaka [17] underestimates the measured data. The formulae of FEMA and NOAA [11], ASCE [13], and Matsutomi [15] provide values near the measured values. Since the formulae of Mizutani et al. [18], Arikawa et al. [20], and Ikeno et al. [22] show good correlations with the measured values, provided the formulae are improved so that parameters concerning rigidity of the driftage and the collided body can be decided adequately, the calculation accuracy of the improved formulae will improve.
section_name: (d) Measured impact force and the collision velocity

section_content: Unique large-scale experiments with all necessary information were conducted by Takabatake et al. [24] (Central Research Institute of Electric Power Industry). On the other hand, when using the formula of Mizutani et al. [18], 1000 kg/m 3 was used for ρw, and half The formula of the Japan Road Association [16] overestimates, and the formula of Ikeno and Tanaka [17] underestimates the measured data. The formulae of FEMA and NOAA [11], ASCE [13], and Matsutomi [15] provide values near the measured values. Since the formulae of Mizutani et al. [18], Arikawa et al. [20], and Ikeno et al. [22] show good correlations with the measured values, provided the formulae are improved so that parameters concerning rigidity of the driftage and the collided body can be decided adequately, the calculation accuracy of the improved formulae will improve.
section_name: Car

section_content: Unique large-scale experiments with all necessary information were conducted by Takabatake et al. [24] (Central Research Institute of Electric Power Industry). On the other hand, when using the formula of Mizutani et al. [18], 1000 kg/m 3 was used for ρ w, and half of the height of their car was used for η m. Furthermore, the impulse duration of ASCE [13] was used for dt. 

The calculation conditions are as follows:

(a) Parameters of the actual car 1. 48 m high, 1. 40 m wide, 3. 30 m long, 316. 0 kg, 83,700 N/m 2 yield stress (from data when the buckling occurred in Takabatake et al. [24] ), 0. 3 Poisson's ratio. The elasticity modulus was calculated by using the empirical equation of the axial stiffness (kN/m) (=457. 41 × collision velocity-170. 79) obtained from measured data of Takabatake et al. [24]. 

(b) Parameters of the collided body (steel) 0. 3 Poisson's ratio, 200,000 MN/m 2 elasticity modulus.
section_name: Car

section_content: Each apparent mass coefficient C MA was set up according to each proposer's recommended value (1. 7 for Matsutomi [15], 2. 0 for Ikeno and Tanaka [17], 2. 0 for Ikeno et al. [22] ). The coefficient C of FEMA and NOAA [11] was 0. 4. The effective stiffness k was set up using the empirical equation obtained from measured data of Takabatake et al. [24].
section_name: (c) Coefficients in the formula

section_content: The measured impact force and the collision velocities are 0. 003 MN and 0. 39 m/s, 0. 024 MN and 0. 98 m/s, 0. 040 MN and 1. 78 m/s, 0. 042 MN and 1. 86 m/s, 0. 053 MN and 2. 07 m/s, 0. 039 MN and 2. 12 m/s, 0. 087 MN and 3. 90 m/s, respectively. 

The verification results are shown in Figure 7. 

of the height of their car was used for ηm. Furthermore, the impulse duration of ASCE [13] was used for dt. The calculation conditions are as follows:

(a) Parameters of the actual car 1. 48 m high, 1. 40 m wide, 3. 30 m long, 316. 0 kg, 83,700 N/m 2 yield stress (from data when the buckling occurred in Takabatake et al. [24] ), 0. 3 Poisson's ratio. The elasticity modulus was calculated by using the empirical equation of the axial stiffness (kN/m) (= 457. 41× collision velocity-170. 79) obtained from measured data of Takabatake et al. [24]. 

(b) Parameters of the collided body (steel) 0. 3 Poisson's ratio, 200,000 MN/m 2 elasticity modulus. 

(c) Coefficients in the formula Each apparent mass coefficient CMA was set up according to each proposer's recommended value (1. 7 for Matsutomi [15], 2. 0 for Ikeno and Tanaka [17], 2. 0 for Ikeno et al. [22] ). The coefficient C of FEMA and NOAA [11] was 0. 4. The effective stiffness k was set up using the empirical equation obtained from measured data of Takabatake et al. [24]. 

(d) Measured impact force and the collision velocity The measured impact force and the collision velocities are 0. 003 MN and 0. 39 m/s, 0. 024 MN and 0. 98 m/s, 0. 040 MN and 1. 78 m/s, 0. 042 MN and 1. 86 m/s, 0. 053 MN and 2. 07 m/s, 0. 039 MN and 2. 12 m/s, 0. 087 MN and 3. 90 m/s, respectively. 

The verification results are shown in Figure 7. 

The formula of the Japan Road Association [16] overestimates the experimental values. The formulae of FEMA and NOAA [11] and ASCE [13] result in values near the measurements. Although the formulae of Matsutomi [15], Ikeno and Tanaka [17], Mizutani et al. [18], Arikawa et al. [20], and Ikeno et al. [22] underestimate the measurements, their predictions show a good correlation with measured values. Therefore, provided their formulae are improved so that parameters concerning the rigidity of the driftage and the collided body can be decided adequately, the calculation accuracy of the improved formulae will improve.
section_name: (d) Measured impact force and the collision velocity

section_content: The large-scale experiments of Arita [25] (National Maritime Research Institute) were the only ones with all the necessary information. Even though in his experiment, a full- The formula of the Japan Road Association [16] overestimates the experimental values. The formulae of FEMA and NOAA [11] and ASCE [13] result in values near the measurements. Although the formulae of Matsutomi [15], Ikeno and Tanaka [17], Mizutani et al. [18], Arikawa et al. [20], and Ikeno et al. [22] underestimate the measurements, their predictions show a good correlation with measured values. Therefore, provided their formulae are improved so that parameters concerning the rigidity of the driftage and the collided body can be decided adequately, the calculation accuracy of the improved formulae will improve.
section_name: Ship

section_content: The large-scale experiments of Arita [25] (National Maritime Research Institute) were the only ones with all the necessary information. Even though in his experiment, a full-scale steel bow model collides against a steel pier model, it was also implemented in air, and the information, which we can use, is only a semi-theoretical formula based on data of this experiment. Moreover, since this formula consists only of the gross tonnage and the collision speed, this formula can only be applied when the conditions of the driftage and the collided body are in agreement with those of the experiment (other similar formulae cannot be used as general-purpose formulae owing to the same reason). When using the formula of Mizutani et al. [19], the impulse duration of ASCE [13] was used for dt. 

The calculation conditions are as follows:

(a) Parameters of the ship (500 GTs)

9. 6 m wide, 50. 0 m long, 1,500,000 kg, 0. 51 m in stop distance (from the recommended values of the Bridge and Engineering Association [14]), 0. 3 Poisson's ratio, 5,308,660 N/m in axial stiffness (from the data table for the axial stiffness of Takabatake et al. [24] ), 5,760,000 N/m 2 elasticity modulus calculated from the axial stiffness by inverse analysis.
section_name: Ship

section_content: Each apparent mass coefficient C MA was set up according to each proposer's recommended value (1. 5 for Ikeno and Tanaka [17], 1. 0 for Ikeno et al. [22] ). The coefficient C of FEMA and NOAA [11] was 0. 4. The effective stiffness k was set up using the data table for the axial stiffness of Takabatake et al. [24].
section_name: (c) Coefficients in the formula

section_content: The measured impact force and the collision velocities are 1. 96 MN and 0. 8 m/s, 2. 94 MN and 1. 2 m/s, 3. 92 MN and 1. 6 m/s, 4. 90 MN and 2. 0 m/s, 5. 88 MN and 2. 4 m/s, respectively. 

The formula of Ikeno and Tanaka [17] underestimates the experimental values. !The formula of Mizutani et al. [19] results in values very close to the experimental ones. 

Although the formulae of FEMA and NOAA [11] and ASCE [13] overestimate the experimental values, their formulae have given values close to the experimental ones. 

Although the formulae of the Bridge and Engineering Association [14], Japan Road Association [16], Arikawa et al. [20], and Ikeno et al. [22] underestimate the measured values, their formulae show good correlations with the experimental values. Therefore, provided their formulae are improved so that parameters concerning the rigidity of the driftage and the collided body can be decided adequately, the calculation accuracy of the improved formulae will improve. 

The reason why the calculated values of each formula in Figure 5 are not aligned along the straight line as in Figures 6 7 8 is because the accuracies of the measured values are low since wood is a natural material and the collision angle of wood relative to the collided body is a sensitive parameter for the result. Moreover, the elasticity modulus and the stiffness of wood also vary significantly. In other cases, the materials of the driftage were fixed, and especially in the case of the full-scale air experiments, since collision angles were right-angled relatively to the collided body, there is no difference due to the collision angles in the measured impact forces. 

The verification results are shown in Figure 8. 

Each apparent mass coefficient CMA was set up according to each proposer's recommended value (1. 5 for Ikeno and Tanaka [17], 1. 0 for Ikeno et al. [22] ). The coefficient C of FEMA and NOAA [11] was 0. 4. The effective stiffness k was set up using the data table for the axial stiffness of Takabatake et al. [24]. 

(d) Measured impact force and the collision velocity The measured impact force and the collision velocities are 1. 96 MN and 0. 8 m/s, 2. 94 MN and 1. 2 m/s, 3. 92 MN and 1. 6 m/s, 4. 90 MN and 2. 0 m/s, 5. 88 MN and 2. 4 m/s, respectively. 

The verification results are shown in Figure 8. The formula of Ikeno and Tanaka [17] underestimates the experimental values. The formula of Mizutani et al. [19] results in values very close to the experimental ones.
section_name: (d) Measured impact force and the collision velocity

section_content: The collision velocity of the experimental data used for comparison with the previous section was not fast (0. 39 m/s-4. 1 m/s), and the range was small. Therefore, to check the validity of the calculation range of each formula, in the collision velocity range of 1 m/s-15 m/s, a general size was set up for each driftage type, and the impact force was calculated. 

In addition, when the specific gravity of the driftage is larger than 1, the collision velocity must be lower than the tsunami velocity. Moreover, when a tsunami flow collides with a body with a large vertical plane, a backward flow occurs near the body. Therefore, as the driftage approaches the body, the collision velocity must become small by the backward flow. However, these influences are not considered. Each apparent mass coefficient C MA was set up according to each proposer's recommended value. (1. 7 for Matsutomi [15], 2. 0 for Ikeno and Tanaka [17], 1. 7 for Arikawa et al. [20], 2. 0 for Ikeno et al. [22] ). The coefficient C of FEMA and NOAA [11] was 0. 0. The effective stiffness k of the slender wood was 20. 1 MN/m (this value was calculated using the elasticity modulus of the slender wood), and the effective stiffness k of the thick wood was 27. 5 MN/m. 

The results are shown in Figure 9. The formula of Ikeno and Tanaka [17] gives too small values compared with the mean values of all formulae. Each apparent mass coefficient CMA was set up according to each proposer's recommended value (1. 7 for Matsutomi [15], 4. 0 for Ikeno and Tanaka [17], 1. 0 for Ikeno et al. [22] ). The constant 0. 2 and 60 MN/m recommended by FEMA and NOAA [11] were set up for the coefficient C and the effective stiffness k of FEMA and NOAA [11], respectively. Moreover, 29. 8 MN/m recommended by ASCE [13] was set up for the effective stiffness k of ASCE [13]. 

When using the formula of Mizutani et al. [18], the density ρw was 1000 kg/m 3 and half of the height of their containers was used for ηm. Furthermore, the impulse duration of ASCE [13] was used for dt. 

The results are shown in Figure 10. The formula of the Japan Road Association [16] gives too big values as compared with the mean values of all formulae. 2. 3. 3. Car (a) Parameters of an actual car 1. 80 m high, 2. 00 m wide, 5. 00 m long, 2500 kg, 47,800 N/m 2 yield stress (from data when the buckling occurred in Takabatake et al. [24] ), 0. 3 Poisson's ratio. The elastic modulus was calculated by using the empirical equation of the axial stiffness (kN/m) (=457. 41 Each apparent mass coefficient C MA was set up according to each proposer's recommended value (1. 7 for Matsutomi [15], 4. 0 for Ikeno and Tanaka [17], 1. 0 for Ikeno et al. [22] ). The constant 0. 2 and 60 MN/m recommended by FEMA and NOAA [11] were set up for the coefficient C and the effective stiffness k of FEMA and NOAA [11], respectively. Moreover, 29. 8 MN/m recommended by ASCE [13] was set up for the effective stiffness k of ASCE [13]. 

When using the formula of Mizutani et al. [18], the density ρ w was 1000 kg/m 3 and half of the height of their containers was used for η m. Furthermore, the impulse duration of ASCE [13] was used for dt. 

The results are shown in Figure 10. The formula of the Japan Road Association [16] gives too big values as compared with the mean values of all formulae. (c) Coefficients of the formula Each apparent mass coefficient CMA was set up according to each proposer's recommendation value (1. 7 for Matsutomi [15], 2. 0 for Ikeno and Tanaka [17], 2. 0 for Ikeno et al. [22] ). The coefficient C of FEMA and NOAA [11] was 0. 4. The effective stiffness k was set up using the empirical equation obtained from measured data of Takabatake et al. [24]. 

Moreover, when using the formula of Mizutani et al. [18], the density ρw was 1000 kg/m 3, and half of the height of the car was used for ηm. Furthermore, the impulse duration of ASCE [13] was used for dt. 

The calculation results are shown in Figure 11. The formula of the Japan Road Association [16] gives too big values as compared with the mean values of all formulae. 2. 3. 4. Ships (a) Parameters of a 5 GT ship 2. 8 m wide, 11. 0 m long, 15,000 kg, 0. 04 m in stop distance (from recommended values of the Bridge and Engineering Association [14]), 0. 3 Poisson's ratio, 1144,950 N/m in axial stiffness (from the data table for the axial stiffness of Takabatake et al. [24] ), 3210,000 N/m 2 elasticity modulus calculated from the axial stiffness by inverse analysis. 

(b) Parameters of a 500 GT ship 9. 4 m wide, 55. 0 m long, 1500,000 kg, 0. 51 m in stop distance (from recommended values of the Bridge and Engineering Association [14]), 0. 3 Poisson's ratio, 5308,660 N/m in axial stiffness (from the data table for the axial stiffness of Takabatake et al. [24] ), 6610,000 N/m 2 elasticity modulus calculated from the axial stiffness by inverse analysis. (c) Coefficients of the formula Each apparent mass coefficient C MA was set up according to each proposer's recommendation value (1. 7 for Matsutomi [15], 2. 0 for Ikeno and Tanaka [17], 2. 0 for Ikeno et al. [22] ). The coefficient C of FEMA and NOAA [11] was 0. 4. The effective stiffness k was set up using the empirical equation obtained from measured data of Takabatake et al. [24]. 

Moreover, when using the formula of Mizutani et al. [18], the density ρ w was 1000 kg/m 3, and half of the height of the car was used for η m. Furthermore, the impulse duration of ASCE [13] was used for dt. 

The calculation results are shown in Figure 11. The formula of the Japan Road Association [16] gives too big values as compared with the mean values of all formulae. [22] ). The coefficient C of FEMA and NOAA [11] was 0. 4. The effective stiffness k was set up using the data table for the axial stiffness of Takabatake et al. [24]. 

When using the formula of Mizutani et al. [19], the impulse duration of ASCE [13] was used for dt. 

The calculation results are shown in Figure 12. The formula of the Bridge and Engineering Association [14] gives too big values as compared with the mean values. The formula of the Japan Road Association [16] gives small values as compared with the mean values of all formulae.
section_name: Range of Impact Force of Each Formula

section_content: Concerning the calculated and measured values of the impact force in Section 2. 2, we classified formulae, which considerably overestimate or underestimate as "×", formulae, which overestimate moderately or give values near the measured values as "〇" and formulae, which give intermediate values of these formulae as "△". Moreover, in examining the range of the calculated values in Section 2. 3, since true values are not known, we classified formulae, which give too large or too small values compared with the mean values of all formulae as "×". Then, we judged that only formulae classified as "〇" in Section 2. 2 and were not classified as "×" in Section 2. 3 should be used because their reliability is high, and formulae, which were classified as "×" in both sections should not be used. These results are shown in Table 4. 

The following can be concluded based on Table 4 :

(1) In the case of driftwood, the calculation formulae of Matsutomi [15], FEMA and NOAA [11] and ASCE [13] are most reliable. 

(2) In the case of containers, the calculation formulae of Matsutomi [15], Arikawa et al. [20], FEMA and NOAA [11], Ikeno et al. [22], and ASCE [13] are most reliable. 

(3) In the case of cars, the calculation formulae of FEMA and NOAA [11] and ASCE [13] are the most reliable.
section_name: Reliability Examination of Main Formulae

section_content: Each apparent mass coefficient C MA was set up according to each proposer's recommended value (1. 5 for Ikeno and Tanaka [17], 1. 0 for Ikeno et al. [22] ). The coefficient C of FEMA and NOAA [11] was 0. 4. The effective stiffness k was set up using the data table for the axial stiffness of Takabatake et al. [24]. 

When using the formula of Mizutani et al. [19], the impulse duration of ASCE [13] was used for dt. 

The calculation results are shown in Figure 12. The formula of the Bridge and Engineering Association [14] gives too big values as compared with the mean values. The formula of the Japan Road Association [16] gives small values as compared with the mean values of all formulae. 

(〇 indicates "good", △ indicates "ordinary", and × indicates "bad").
section_name: (d) Coefficients in the formula

section_content: Concerning the calculated and measured values of the impact force in Section 2. 2, we classified formulae, which considerably overestimate or underestimate as "×", formulae, which overestimate moderately or give values near the measured values as " " and formulae, which give intermediate values of these formulae as " ". Moreover, in examining the range of the calculated values in Section 2. 3, since true values are not known, we classified formulae, which give too large or too small values compared with the mean values of all formulae as "×". Then, we judged that only formulae classified as " " in Section 2. 2 and were not classified as "×" in Section 2. 3 should be used because their reliability is high, and formulae, which were classified as "×" in both sections should not be used. These results are shown in Table 4. The following can be concluded based on Table 4 :

(1) In the case of driftwood, the calculation formulae of Matsutomi [15], FEMA and NOAA [11] and ASCE [13] are most reliable. (2) In the case of containers, the calculation formulae of Matsutomi [15], Arikawa et al. [20],

FEMA and NOAA [11], Ikeno et al. [22], and ASCE [13] are most reliable. (3) In the case of cars, the calculation formulae of FEMA and NOAA [11] and ASCE [13] are the most reliable. (4) In the case of ships, the calculation formulae of Mizutani [19], FEMA and NOAA [11],

and ASCE [13] are the most reliable. 

However, when FEMA and NOAA [11] and ASCE [13] are used, although the appropriate stiffness is very important, it is very difficult to estimate the appropriate stiffness of complicated structures. Moreover, the stiffness changes according to the change of collision velocity. Therefore, we used the stiffness recommended by FEMA and NOAA [11] and ASCE [13] in the case of a 6. 1 m container and a 12. 2 m container, and the recommended stiffness of Takabatake et al. [24] or Kaida and Kihara [6] for cars and ships. 

Moreover, since the formulae of the Bridge and Engineering Association [14], Matsutomi [15], Mizutani et al. [18, 19], Arikawa et al. [20], and Ikeno et al. [22] show good correlations with the measured values, if the formulae were improved so that parameters concerning the rigidity of the driftage and the collided body could be considered, the reliability of the formulae becomes better, and the applicability range also increases. 

On the other hand, since the formula of the Japan Road Association [16] cannot make the difference in rigidity of the driftage and the collided body into consideration and tends to give an excessive impact force, it is better to avoid its use. In addition, since the formula of Ikeno and Tanaka [17] can also not make the difference in rigidity of the driftage and the collided body into consideration and tends to give a too small impact force, it is better to avoid its use.
section_name: Reliability Examination of Main Formulae

section_content: There are many formulae for calculating the impact forces from driftage. Since there are considerable differences between the calculated values of these formulae, the reliability of each formula was assessed here. We summarized the main calculation formulae of the impact force proposed by 2019 and compared the calculation values of these formulae with the measured data from large-scale experiments. Furthermore, we checked the range of the calculation values of each formula up to 15 m/s in collision velocity and then clarified the following points:

(1) In the case of driftwood, the formulae of Matsutomi [15], FEMA and NOAA [11], and ASCE [13] are the most reliable; (2) In the case of containers, the formulae of Matsutomi [15], Arikawa et al. [20], FEMA and NOAA [11], Ikeno et al. [22], and ASCE [13] are most reliable; (3) In the case of cars, the formulae of FEMA and NOAA [11] and ASCE [13] are the most reliable; (4) In the case of ships, the formulae of Mizutani [19], FEMA and NOAA [11], and ASCE [13] are the most reliable. 

Here, except for the formulae of the Japan Road Association [16] and Ikeno and Tanaka [17], the stiffness or parameters governed by rigidity (stop distance of ships, yield stress, collision time, elasticity modulus) must be set appropriately:

(a) When FEMA and NOAA [11] and ASCE [13] are used, the stiffness recommended by FEMA and NOAA [11] and ASCE [13] can be used for the 6. 1 m container and the 12. 2 m container, and the stiffness recommended by Takabatake et al. [24] or Kaida and Kihara [6] can be used for cars and ships; (b) Since the formulae of the Bridge and Engineering Association [14], Matsutomi [15],

Mizutani et al. [18, 19], Arikawa et al. [20], and Ikeno et al. [22] show good correlations with the measured values, provided the formulae are improved so that parameters concerning the rigidity of the driftage and the collided body can be considered, the reliability of the formulae is expected to increase, and the applicable scope to become larger; (c) Since the formulae of the Japan Road Association [16] and Ikeno and Tanaka [17] cannot make the difference in rigidity of the driftage and the collided body into consideration, the applicable range of their formulae is limited. Moreover, the former tends to give an excessive impact force, and the latter tends to give a too small impact force. Therefore, it is desirable to avoid using their equations.
section_name: Conclusions

section_name: Data Availability Statement:

section_content: In this research, all photos were taken by authors and all data used for calculations are open data from the references.
section_name: Data Availability Statement:

section_name: Conflicts of Interest:

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest: