Method of Ice Distribution and Data Processing for Ship Ice Resistance Experiment in Broken Ice Field

Abstract

The present disclosure discloses 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, so as to make the experimental data more reliable.

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 with eliminating overlapped ice floes, and recording pre-experiment data; carrying out a ship ice resistance experiment in the broken ice field, and recording experimental data; S2. repeating step S1 at least three times; S3. based on the segmentation hypothesis and stable state requirement in ITTC ice testing guidance document, dividing the broken ice field into multiple subdomains; S4. conducting an outlier data point judgment for the experimental data corresponding to each subdomain to discard the experimental data corresponding to abnormal subdomains, and conducting uncertainty analysis for the experimental data corresponding to remaining subdomains; S5. conducting an image correction for pictures of the broken ice field in the pre-experiment data and dividing the pictures by subdomains to obtain subdomain pictures, and calculating an actual coverage rate of subdomain corresponding to each subdomain picture; and S6. according to the actual coverage rate and design coverage rate of the subdomains, 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 broken ice field, and the ice floes are polyethylene plates with uniform thickness.

3. The method according to claim 1, wherein the principle of dividing the broken ice field into multiple subdomains in step S3 is as follows: a length of the subdomain is 1.5-2.5 times a length of the ship; the ship has at least 10 collisions with the ice floes in each subdomain; and the experimental data corresponding to the subdomain should be able to indicate a general variation trend.

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

5. The method according to claim 4, wherein a detailed steps of step S4 are as follows: selecting one kind of data from the experimental data corresponding to each subdomain, and judging the data according to formula (1); $\begin{array}{cc}{\left(\mathrm{Chauv}.Math.\#\right)}_{\mathrm{Mean}}=.Math.\frac{F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}-\left(\mathrm{Mean\_}.Math..Math.F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}\right)}{\left(\mathrm{STD}.Math..Math.\_.Math..Math.F_{T.Math..Math.\_.Math..Math.\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; if the value of (Chauv #).sub.Mean of the subdomain is not within the range of 1.96-2.13, the sub domain is abnormal and discard the corresponding experimental data, then conduct uncertainty analysis of the experimental data corresponding to the remaining subdomains: $\begin{array}{cc}U\left(F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}\right)=\frac{t*\left(\mathrm{STD}.Math..Math.\_.Math..Math.F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}\right)}{\sqrt{N}}& \left(2\right)\end{array}$ wherein, U(F.sub.T_mean) represents a uncertainty of one kind of data in the experimental data corresponding to the remaining subdomains, t represents a empirical coefficient, (STD_F.sub.T_mean) represents the 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 Mean_F.sub.T_mean as the following form of percentage: $\begin{array}{cc}{U}_p\left(F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}\right)=\frac{U\left(F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}\right)}{\mathrm{Mean\_}.Math..Math.F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}}*1.Math.0.Math.0& \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. The method according to claim 5, wherein the detailed steps of step S5 are as follows: using MATLAB to carry out geometric correction on the broken ice field pictures in the pre-experimental data to eliminate perspective deformation and dividing the pictures by subdomains to obtain subdomain pictures; calculating the actual coverage rate of each subdomain picture: using MATLAB to carry out image binarization processing for the subdomain pictures, and adjusting the threshold value as high as possible; calculating the actual coverage rate of each subdomain picture after binarization processing: actual coverage rate=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 the obtained actual coverage rate equals the actual coverage rate of the subdomain corresponding to the picture of the subdomain.

7. The method according to claim 6, wherein in step S6, a final ice resistance corresponding to each remaining subdomain is obtained by the following formula: $\begin{array}{cc}{F}_{\mathrm{ice}.Math..Math.\_.Math..Math.\mathrm{final}}=F_{\mathrm{ice}.Math..Math.\_.Math..Math.\mathrm{test}}.Math.\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 a ice resistance in the experimental data corresponding to each remaining subdomain, coverage_design represents a design coverage rate of the subdomain, coverage_modified represents the actual coverage rate of the subdomain.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0036] 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.

[0037] 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.

[0038] 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

[0039] 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.

[0040] 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:

[0041] 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;

[0042] the ship ice resistance experiment was carried out in the broken ice field, and the experimental data was recorded;

[0043] S2. step S1 was repeated at least three times;

[0044] 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;

[0045] 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;

[0046] 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;

[0047] 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.

[0048] 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.

[0049] The principle of dividing the broken ice field into multiple subdomains in step S3 is as follows:

[0050] a length of the subdomain is 1.5-2.5 times a length of the ship;

[0051] the ship has at least 10 collisions with the ice floes in the subdomain;

[0052] the experimental data corresponding to the subdomain should be able to indicate a general variation trend.

[0053] The number of the subdomains in step S3 was 10-15.

[0054] The specific steps of step S4 was as follows:

[0055] 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;

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

[0056] 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;

[0057] 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 uncertainty analysis of the experimental data corresponding to the remaining subdomain is conducted:

[00006]$\begin{array}{cc}U\left(F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}\right)=\frac{t*\left(\mathrm{STD}.Math..Math.\_.Math..Math.F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}\right)}{\sqrt{N}}& \left(2\right)\end{array}$

[0058] 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;

[0059] Mean_F.sub.T_mean was expressed as the following form of percentage:

[00007]$\begin{array}{cc}{U}_p\left(F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}\right)=\frac{U\left(F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}\right)}{\mathrm{Mean\_}.Math..Math.F_{T.Math..Math.\_.Math..Math.\mathrm{mean}}}*1.Math.0.Math.0& \left(3\right)\end{array}$

[0060] 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.

[0061] The detailed steps of step S5 are as follows:

[0062] 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;

[0063] the actual coverage rate of each subdomain picture was calculated:

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

[0064] 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.

[0065] In step S6, the final ice resistance corresponding to each remaining subdomain is obtained by the following formula:

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

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_modified was the actual coverage rate of the subdomain.

[0066] 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.