Method of ice distribution and data processing for ship ice resistance experiment in broken ice field

Abstract

In a method of ice distribution and data processing for ship ice resistance experiment in broken ice field, the ice floe is uniformly distributed in the broken ice field before each ship ice resistance experiment with eliminating the mutually overlapped broken ice to make it close to the experimental design working condition, and independent repeated experiments are carried out on each ship ice resistance experiment to reduce the influence of accidental factors; after finishing the experiment, the experimental data are divided into intervals and the experimental data corresponding to unreasonable subdomains are eliminated by mathematical statistics, and the uncertainty analysis is carried out; image correction is carried out for pictures of in broken ice field in the pre-experimental data, and the pictures are divided to obtain subdomain pictures to calculate the actual coverage rate of each subdomain; the ice resistance in the experimental data corresponding to each remaining subdomain is modified.

Claims

1. A method of ice distribution and data processing for ship ice resistance experiment in broken ice field, comprising the following steps: S1. uniformly distributing ice floes in a broken ice field, eliminating overlapped ice floes, and recording pre-experiment data; carrying out a ship ice resistance experiment by passing a ship through the broken ice field, and recording experimental data; S2. repeating step S1 at least three times; S3. dividing the broken ice field into a plurality of subdomains; S4. conducting an outlier data point judgment for the experimental data corresponding to each subdomain, discarding experimental data corresponding to abnormal subdomains among the plurality of subdomains, and conducting uncertainty analysis for the experimental data corresponding to remaining subdomains among the plurality of subdomains; S5. performing a geometric correction on pictures of the broken ice field in the pre-experimental data to eliminate perspective deformation; assigning at least one subdomain picture to each of the plurality of subdomains; performing image binarization processing to the at least one subdomain picture for each subdomain; calculating an actual coverage rate for each subdomain picture after binarization processing, wherein the actual coverage rate equals a number of white pixels in the picture after binarization processing divided by a total number of pixels in the picture after binarization processing; and S6. according to the actual coverage rate and a design coverage rate of each subdomain, modifying ice resistance in the experimental data corresponding to each remaining subdomain to obtain a final ice resistance corresponding to each remaining subdomain.

2. The method according to claim 1, wherein the ship passes through the broken ice field at a constant speed during the ship ice resistance experiment in the broken ice field, and the ice floes are polyethylene plates with a uniform thickness.

3. The method according to claim 1, wherein, in step S3: a length of each subdomain is 1.5-2.5 times a length of the ship; and the ship has at least 10 collisions with the ice floes in each subdomain.

4. The method according to claim 3, wherein a number of the plurality of subdomains in step S3 is 10-15.

5. The method according to claim 4, wherein step S4 comprises: selecting one kind of data from the experimental data corresponding to each subdomain, and judging the data according to equation (1), $\begin{array}{cc}{\left(\mathrm{Chauv}\#\right)}_{\mathrm{Mean}}=.Math.\frac{F_{T\_\mathrm{mean}}-\left(\mathrm{Mean\_}{F}_{T\_\mathrm{mean}}\right)}{\left(\mathrm{STD}\_F_{T\_\mathrm{mean}}\right)}.Math.,& \left(1\right)\end{array}$ wherein, (Chauv #).sub.Mean represents a Chauvenet value of the kind of data, F.sub.T_mean represents a mean value of the selected data of each subdomain, Mean_F.sub.T_mean represents a mean value of the mean values of the selected data of all subdomains, STD_F.sub.T_mean represents a variance of the mean values of the selected data of all subdomains; when the value of (Chauv #).sub.Mean of the subdomain is not within the range of 1.96-2.13, determining the subdomain as abnormal and discarding the corresponding experimental data, and conducting uncertainty analysis of the experimental data corresponding to the remaining subdomains according to equation (2): $\begin{array}{cc}U\left(F_{T\_\mathrm{mean}}\right)=\frac{t*\left(\mathrm{STD}\_F_{T\_\mathrm{mean}}\right)}{\sqrt{N}},& \left(2\right)\end{array}$ wherein, U(F.sub.T_mean) represents an uncertainty of one kind of data in the experimental data corresponding to the remaining subdomains, t represents an empirical coefficient, (STD_F.sub.T_mean) represents a variance of the mean values of one kind of data in the experimental data corresponding to all the remaining subdomains, and N represents a number of the remaining subdomains; and expressing U(F.sub.T_mean) as a percentage value according to equation (3): $\begin{array}{cc}{U}_p\left(F_{T\_\mathrm{mean}}\right)=\frac{U\left(F_{T\_\mathrm{mean}}\right)}{\mathrm{Mean\_}{F}_{T\_\mathrm{mean}}}*100,& \left(3\right)\end{array}$ wherein Mean_F.sub.T_mean represents the mean value of one kind of data in the experimental data corresponding to all the remaining subdomains.

6. A method of ice distribution and data processing for ship ice resistance experiment in broken ice field, comprising: S1. uniformly distributing ice floes in a broken ice field, eliminating overlapped ice floes, and recording pre-experiment data; carrying out a ship ice resistance experiment by passing a ship through the broken ice field, and recording experimental data; S2. repeating step S1 at least three times; S3. dividing the broken ice field into a plurality of subdomains; S4. conducting an outlier data point judgment for the experimental data corresponding to each subdomain, discarding experimental data corresponding to abnormal subdomains among the plurality of subdomains, and conducting uncertainty analysis for the experimental data corresponding to remaining subdomains among the plurality of subdomains; S5. conducting an image correction for pictures of the broken ice field in the pre-experiment data and assigned one or more pictures to each of the plurality of subdomains, and calculating an actual coverage rate of each subdomain corresponding to the one or more pictures assigned thereto; and S6. according to the actual coverage rate and a design coverage rate of each subdomain, modifying ice resistance in the experimental data corresponding to each remaining subdomain to obtain a final ice resistance corresponding to each remaining subdomain, wherein the final ice resistance corresponding to each remaining subdomain is obtained according to equation (4): $\begin{array}{cc}{F}_{\mathrm{ice}\_\mathrm{final}}=F_{\mathrm{ice}\_\mathrm{test}}\frac{\mathrm{coverage\_design}}{\mathrm{coverage\_modified}},& \left(4\right)\end{array}$ wherein, F.sub.ice_final represents the final ice resistance corresponding to each remaining subdomain, F.sub.ice_test represents an ice resistance in the experimental data corresponding to each remaining subdomain, coverage_design represents the design coverage rate of the subdomain, coverage_modificed represents the actual coverage rate of the subdomain.

7. The method according to claim 6, wherein the ship passes through the broken ice field at a constant speed during the ship ice resistance experiment in the broken ice field, and the ice floes are polyethylene plates with a uniform thickness.

8. The method according to claim 6, wherein, in step S3, a length of each subdomain is 1.5-2.5 times a length of the ship; and the ship has at least 10 collisions with the ice floes in each subdomain.

9. The method according to claim 8, wherein a number of the plurality of subdomains in step S3 is 10-15.

10. The method according to claim 6, wherein step S4 comprises: selecting one kind of data from the experimental data corresponding to each subdomain, and judging the data according to equation (1), $\begin{array}{cc}{\left(\mathrm{Chauv}\#\right)}_{\mathrm{Mean}}=.Math.\frac{F_{T\_\mathrm{mean}}-\left(\mathrm{Mean\_}{F}_{T\_\mathrm{mean}}\right)}{\left(\mathrm{STD}\_F_{T\_\mathrm{mean}}\right)}.Math.,& \left(1\right)\end{array}$ wherein, (Chauv #).sub.Mean represents a Chauvenet value of the kind of data, F.sub.T_mean represents a mean value of the selected data of each subdomain, Mean_F.sub.T_mean represents a mean value of the mean values of the selected data of all subdomains, STD_F.sub.T_mean represents a variance of the mean values of the selected data of all subdomains; when the value of (Chauv #).sub.Mean of the subdomain is not within the range of 1.96-2.13, determining the subdomain as abnormal and discarding the corresponding experimental data, and conducting uncertainty analysis of the experimental data corresponding to the remaining subdomains according to equation (2): $\begin{array}{cc}U\left(F_{T\_\mathrm{mean}}\right)=\frac{t*\left(\mathrm{STD}\_F_{T\_\mathrm{mean}}\right)}{\sqrt{N}},& \left(2\right)\end{array}$ wherein, U(F.sub.T_mean) represents an uncertainty of one kind of data in the experimental data corresponding to the remaining subdomains, t represents an empirical coefficient, (STD_F.sub.T_mean) represents a variance of the mean values of one kind of data in the experimental data corresponding to all the remaining subdomains, and N represents a number of the remaining subdomains; and expressing U(F.sub.T_mean) as a percentage value according to equation (3): $\begin{array}{cc}{U}_p\left(F_{T\_\mathrm{mean}}\right)=\frac{U\left(F_{T\_\mathrm{mean}}\right)}{\mathrm{Mean\_}{F}_{T\_\mathrm{mean}}}*100,& \left(3\right)\end{array}$ wherein Mean_F.sub.T_mean represents the mean value of one kind of data in the experimental data corresponding to all the remaining subdomains.

11. The method according to claim 6, wherein step S5 comprises: performing a geometric correction on pictures of the broken ice field in the pre-experimental data to eliminate perspective deformation; assigning at least one subdomain picture to each of the plurality of subdomains; performing image binarization processing to the at least one subdomain picture for each subdomain; and calculating an actual coverage rate for each subdomain picture after binarization processing, wherein the actual coverage rate equals a number of white pixels in the picture after binarization processing divided by a total number of pixels in the picture after binarization processing.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) In order to illustrate more clearly the technical solutions in the embodiments of the present disclosure or in the prior art, a brief description is given hereinafter to the accompany drawings that may be used in the embodiments or the prior art, the drawings described below are some embodiments of the present disclosure, for the ordinary person skilled in the art, other drawings can be obtained according to these drawings without paying any creative labor.

(2) FIG. 1 is a flow diagram of the method of ice distribution and data processing for ship ice resistance experiment in broken ice field in the embodiments of the present disclosure.

(3) FIG. 2 is a flow diagram of image correction, subdomain division, binarization processing and actual coverage rate calculation of broken ice field pictures in the pre-experiment data in the embodiment of the present disclosure.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

(4) In order to make the purpose, technical solutions and advantages of the present disclosure clearer, a clear and complete description of the embodiments in the present disclosure is given herein after in combination with the accompany drawings of the embodiment, the embodiments described below are part embodiments of the present disclosure and not all of them. Based on the embodiments in the present disclosure, all other embodiments obtained by ordinary person skilled in the art without paying any creative labor belong to the protection scope of the present disclosure.

(5) As shown in FIG. 1 and FIG. 2, a method of ice distribution and data processing method for ship ice resistance experiment in broken ice field, includes the following steps:

(6) S1. the ice floes were uniformly distribute in the broken ice field and the overlapped ice floes was eliminate, then recorded the pre-experiment data;

(7) the ship ice resistance experiment was carried out in the broken ice field, and the experimental data was recorded;

(8) S2. step S1 was repeated at least three times;

(9) S3. based on the segmentation hypothesis and stable state requirement in ITTC ice testing guidance document (ICE Testing Experimental Uncertainty Analysis for Ship Resistance in Ice Tank Testing), the broken ice field was divided into multiple subdomains;

(10) S4. an outlier data point judgment was conducted for the experimental data corresponding to each subdomain to discard the experimental data corresponding to the abnormal subdomain, and uncertainty analysis was conducted for the experimental data corresponding to the remaining subdomains;

(11) S5. an image correction was conducted for the pictures of the broken ice field in the pre-experiment data and the pictures were divided by subdomains to obtain subdomain pictures, then the actual coverage rate of the subdomain corresponding to each subdomain picture was calculated;

(12) S6. according to the actual coverage rate and the design coverage rate of the subdomain, the ice resistance in the experimental data corresponding to each remaining subdomain was modified to obtain the final ice resistance in broken ice field corresponding to each remaining subdomain.

(13) During the ship ice resistance experiment in broken ice field, the ship passed through the broken ice field at a constant speed, and the ice floes were polyethylene plates with uniform thickness.

(14) The principle of dividing the broken ice field into multiple subdomains in step S3 is as follows:

(15) a length of the subdomain is 1.5-2.5 times a length of the ship;

(16) the ship has at least 10 collisions with the ice floes in the subdomain;

(17) the experimental data corresponding to the subdomain should be able to indicate a general variation trend.

(18) The number of the subdomains in step S3 was 10-15.

(19) The specific steps of step S4 was as follows:

(20) one kind of data was selected from the experimental data corresponding to each subdomain and judged according to formula (1); the data selected in this embodiment was the towing force;

(21) $\begin{array}{cc}{\left(\mathrm{Chauv}\#\right)}_{\mathrm{Mean}}=.Math.\frac{F_{T\_\mathrm{mean}}-\left(\mathrm{Mean\_}{F}_{T\_\mathrm{mean}}\right)}{\left(\mathrm{STD}\_F_{T\_\mathrm{mean}}\right)}.Math.& \left(1\right)\end{array}$

(22) wherein, (Chauv #).sub.Mean was the Chauvenet value of the mean towing force, F.sub.T_mean was the mean value of the towing force of each subdomain, Mean_F.sub.T_mean was the mean value of the mean values (i.e., F.sub.T_mean) of the towing forces of all subdomains, STD_F.sub.T_mean was the variance of the mean values (i.e., F.sub.T_mean) of the towing forces of all subdomains;

(23) if the value of (Chauv #).sub.Mean of the subdomain is not within the range of 1.96-2.13, the subdomain is abnormal and the corresponding experimental data is discarded, then an uncertainy analysis of the experimental data corresponding to the remaining subdomain is conducted:

(24) $\begin{array}{cc}U\left(F_{T\_\mathrm{mean}}\right)=\frac{t*\left(\mathrm{STD}\_F_{T\_\mathrm{mean}}\right)}{\sqrt{N}}& \left(2\right)\end{array}$

(25) wherein, U(F.sub.T_mean) was the uncertainty of the towing force in the experimental data corresponding to the remaining subdomains, t was the empirical coefficient, in this embodiment, t=2, (STD_F.sub.T_mean) was the variance of the mean values of the towing forces in the experimental data corresponding to all the remaining subdomains, N was the number of the remaining subdomains;

(26) U.sub.p(F.sub.T_mean) was expressed as the following form of percentage;

(27) $\begin{array}{cc}{U}_p\left(F_{T\_\mathrm{mean}}\right)=\frac{U\left(F_{T\_\mathrm{mean}}\right)}{\mathrm{Mean\_}{F}_{T\_\mathrm{mean}}}*100& \left(3\right)\end{array}$

(28) wherein, Mean_F.sub.T_mean was the mean value of the mean values of the towing forces in the experimental data corresponding to all the remaining subdomains.

(29) The detailed steps of step S5 are as follows:

(30) geometric correction was carried out by using MATLAB on the broken ice field pictures in the pre-experimental data to eliminate the perspective deformation, and the pictures were divided by subdomain to obtain subdomain pictures;

(31) the actual coverage rate of each subdomain picture was calculated:

(32) picture binarization processing was carried out by using MATLAB for the subdomain pictures, and the threshold value was adjusted as high as possible;

(33) the actual coverage rate of the picture after binarization processing was calculated: the actual coverage rate=the number of white pixels in the picture after binarization processing divided by the total number of pixels in the picture after binarization processing; the actual coverage rate obtained was equal to the actual coverage rate of the subdomain corresponding to the picture of the subdomain.

(34) In step S6, the final ice resistance corresponding to each remaining subdomain is obtained by the following formula:

(35) $\begin{array}{cc}{F}_{\mathrm{ice}\_\mathrm{final}}=F_{\mathrm{ice}\_\mathrm{test}}\frac{\mathrm{coverage\_design}}{\mathrm{coverage\_modified}}& \left(4\right)\end{array}$

(36) wherein, F.sub.ice_final was the final ice resistance corresponding to each remaining subdomain, F.sub.ice_test was the ice resistance in the experimental data corresponding to each remaining subdomain, coverage_design was the design coverage rate of the subdomain, coverage_modificed was the actual coverage rate of the subdomain.

(37) Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present disclosure without limiting; although the present disclosure is described in detail with reference to the foregoing embodiments, the ordinary skilled in the art shall understand that they may still make amendments to the technical solutions disclosed in the foregoing embodiments, or make equal replacements for some or all of their technical characteristics; these amendments or replacements do not remove the essence of the corresponding technical solutions from the scope of the technical solutions of each embodiment of the present disclosure.