第七章 阻力近似估算 船舶阻力 与推进.docx

1、第七章 阻力近似估算 船舶阻力 与推进第七章 阻力近似估算方法 在船舶设计过程中，特别是在方案设计的初期，当主尺度和船型系数被确定以后，必须要知道主机功率以预报船舶能达到设计航速；如果主机功率已知，则需要估计阻力，以确定船的航速，便于初步分析、比较各种方案的优劣。在此阶段，由于船舶线型尚未确定，因而还不能应用船模试验方法来确定阻力，所以只能用近似方法进行估算。此外在某些不准备作船模试验的小型船舶或航速不重要的船舶的设计过程中，只能用近似方法来确定其阻力值。 近似估算阻力的方法很多，但所有这些方法不外乎是根据船模系列试验结果或者是在总结、分析大量的船模试验和实船试验的基础上得出的。因此可以想象应

2、用近似估算法所得结果的准确程度取决于设计船与母型船或设计船与各图谱所依据的船模系列之间的相似程度。所以为了尽可能提高近似估算的准确性，应该有针对性地选择适当的估算方法。 阻力近似估算方法按计算内容可分为两类：一类是直接近似估算总阻力或有效功率；另一类是估算剩余阻力，而用相当平板公式计算摩擦阻力；如果依阻力近似估算方法的表达形式可分为图谱法和回归公式法两种；若根据估算方法的资料来源进行分类，则可分为船模系列资料估算法、归纳实船和船模资料估算法、母型船数据估算法等三类估算方法。 7-1 船模系列试验资料估算法 这类方法都是根据船模系列试验资料，直接给出阻力图表等，供实际估算应用。 一、泰洛(Tay

3、lor)法 泰洛估算法是根据泰洛标准系列船模试验结果整理得到的。其所用母型船虽为军舰(参见6-1)，但也可用于民用船，特别是双螺旋桨客船的阻力估算。最初的泰洛法其阻力数据绘制成单位排水量剩余阻力的等值线，并均采用英制单位。1954年盖脱勒(Gertler)将泰洛标准组阻力数据重新进行分析整理，并对水温、层流和限制航道的影响分别加以修正，最后整理出一套无量纲剩余阻力系数图表，其中摩擦阻力系数按桑海公式计算。计算所用的船体湿面积可以由无量纲湿面积系数图谱求得。因此该估算法又称泰洛-盖脱勒法。 1泰洛法给出的图谱形式和参数范围 (1) 剩余阻力系数Cr图谱：其给出的函数关系为： Cr = f 1 (

4、Cp，Fr) (7-1) 显然当B/d，Cp一定时，(7-1)式的函数关系可表达为： Cr = f 2 (，Fr) (7-1a)所以给出的图谱形式是：对每一组B/d，以每一棱形系数Cp给出一张图谱，图中以不同排水体积长度系数 /L3为参数，作出剩余阻力系数Cr对Fr(或)的曲线，如图7-1所示。 显然，根据设计船的B/d、Cp、 /L3和Fr值由泰洛图谱可以很方便地求得Cr值。表7-1(a)、表7-1(b)、表7-1(c)为重新分析后的泰洛标准组剩余阻力系数数据摘要。因此亦可由该表查得Cr值。 (2) 无量纲湿面积系数Cs图谱：其给出的函数关系为： Cs = f 3 (Cp，) (7-2)表7

5、-1(a) 泰洛标准系列剩余阻力系数( B/d 2.25)Crl03值Cp1000 /L3Fr0.180.200.220.240.260.280.300.320.501.02.03.04.05.06.07.00.200.270.310.360.420.47.0.530.200.280.320.380.440.490.560.200.280.360.420.490.550.640.200.320.430.520.600.680.800.200.370.500.590.700.820.990.200.400.530.640.730.881.040.240.440.610.770.930.400.7

6、21.001.261.540.601.02.03.04.05.06.07.00.230.280.330.380.450.500.560.230.290.330.400.470.520.590.230.300.380.460.530.600.680.250.370.470.560.640.710.810.280.500.660.800.901.011.100.380.731.021.321.491.671.820.451.011.451.862.180.461.031.542.002.410.701.02.03.04.05.06.07.00.230.300.370.420.490.560.640

7、.230.320.420.520.590.680.770.280.470.600.710.820.900.970.400.710.871.011.101.201.280.550.971.221.431.561.681.820.811.461.972.482.883.203.501.132.303.264.281.232.473.775.020.801.02.03.04.05.06.07.00.450.530.600.670.750.830.900.620.800.961.131.221.311.410.911.291.581.801.962.102.231.301.942.302.542.75

8、1.452.323.023.431.702.632.213.775.382.184.366.64 表7-1(b) 泰洛标准系列剩余阻力系数( B/d 3.00) Crl03值Cp1000 /L3Fr0.180.200.220.240.260.280.300.320.501.02.03.04.05.06.07.00.250.300.360.430.480.550.620.250.300.360.430.500.600.680.270.380.440.500.600.720.880.330.520.630.740.820.961.240.380.610.740.860.961.111.310.38

9、0.600.770.870.971.111.280.430.630.851.021.191.390.590.901.181.501.862.200.601.02.03.04.05.06.07.00.300.380.420.470.550.610.660.300.380.420.490.570.640.740.350.440.500.570.630.740.880.460.580.650.710.760.861.000.560.750.870.930.991.041.200.640.971.201.301.421.631.790.701.171.571.892.202.450.741.281.7

10、42.192.630.701.02.03.04.05.06.07.00.420.490.560.630.720.820.480.580.680.770.870.970.630.820.951.041.111.170.821.101.241.331.401.521.081.401.601.721.851.991.582.152.572.893.133.522.403.404.222.693.985.160.801.02.03.04.05.06.07.00.650.860.930.991.031.070.931.141.251.301.331.381.391.651.771.851.921.992

11、.052.602.822.933.053.142.613.483.034.154.97表7-1(c) 泰洛标准系列剩余阻力系数( B/d 3.75)Crl03值Cp1000 /L3Fr0.180.200.220.240.260.280.300.320.501.02.03.04.05.06.07.00.330.380.470.530.580.670.740.330.400.480.560.610.710.780.370.460.550.680.740.840.940.430.550.670.830.961.121.270.480.630.780.961.141.331.530.470.640.8

12、31.021.191.401.600.460.660.901.161.431.660.580.851.211.612.102.55（续上页）Crl03值Cp1000 /L3Fr0.180.200.220.240.260.280.300.320.601.02.03.04.05.06.07.00.420.500.550.610.700.770.840.420.500.560.630.720.810.900.430.520.590.700.800.901.000.450.560.660.770.870.971.100.490.710.870.961.071.22 1.340.560.921.161.

13、381.581.791.880.621.131.501.832.202.500.701.301.802.242.700.701.02.03.04.05.06.07.00.580.650.710.770.830.920.620.700.790.870.961.050.740.870.981.081.201.300.911.121.261.361.481.581.151.421.641.822.002.081.672.202.622.983.223.392.403.364.052.703.984.980.801.02.03.04.05.06.07.00.760.840.910.971.061.13

14、1.071.141.221.291.381.461.551.681.781.881.982.112.172.522.672.782.913.112.613.233.083.974.68当B/d一定时，(7-2)式成为： Cs = f 4 (Cp，) (7-2a)所以Cs图谱的形式是：以每一组B/d值给出一张图谱，表示为在不同Cp和/L3时的Cs等值曲线，如图7-2所示。 显然，根据设计船的B/d，Cp和/L3值，由该图谱可查得Cs值，进而可按下式求得船体湿面积： S = Cs (7-3) (3) 适用的参数范围 泰洛标准系列的船型参数范围分别为： 第组； B/d2.25， L/B5.014.5

15、， Cb0.440.80， Cp0.480.86 第二组： B/d3.00， L/B4.615.1， Cb0.440.74， Cp0.480.80 第三组： B/d3.75， L/B4.013.4， Cb0.440.80， Cp0.480.86显见，泰洛系列的适用范围很广。2泰洛法阻力估算的具体步骤 应用泰洛法估算阻力的步骤及其注意点可归纳如下： (1) 计算设计船的船型参数值B/d，Cp，/L3和Fr； (2) 求湿面积系数Cs值 一般可由船体型线图求得湿面积。但如果在估算设计船的实船有效功率时，尚无型线图，则可近似认为其湿面积与标准船型的湿面积相等。这样可按给定的B/d值选用与其最接近的两

16、张Cs图谱，并由给定的Cp和 /L3值查得Cs值，然后对B/d进行线性内插求得给定参数值时标准船型的Cs值。 (3) 计算摩擦阻力系数Cf值 在泰洛法中，Cf值是按桑海公式进行计算的，可先计算雷诺数，然后由查表或直接按公式计算得到。应该注意计算雷诺数的船长应为水线长度。粗糙度补贴Cf一般取0.410-3。 (4) 求剩余阻力系数Cr值 根据设计船的参数B/d，Cp，/L3和Fr值，可选对应的Cr图谱，经各参数内插求得Cr值。 应该注意，如果实船的湿面积S 已知时，则可求得实船之湿面积系数Cs为： Cs = S / (7-4)其值如与标准船型湿面积系数Cs不同时，必须对所得到的剩余阻力系数值Cr

17、进行修正，修正后的实船剩余阻力系数Cr为：Cr = Cr Cs / Cs (7-5) 若在计算时不知船的湿面积系数，则可假定CsCs，则有： Cr = Cr (7-5a) (5) 计算总阻力Rts，有效功率Pe值 总阻力系数：Cts CrCfCf 总阻力： Rts Cts 2 S (N) 有效功率： Pe Rts/1000 (kw) 对不同航速，重复上述计算步骤可以得到有效功率曲线。具体计算可列表进行。 最后，应该指出，泰洛标准船型的母型为一巡洋舰，其阻力性能较好，因此对航速较高、船型较瘦的双螺旋桨船用此法计算比较适当；而用于估算一艘商船，所得结果往往偏低。 二、陶德(Todd)法(又称系列6

18、0) 陶德用编号为60的一组不同方形系数母型船模为基础，首先系列变化其浮心纵向位置，然后对不同方形系数的母型船保持其浮心在最佳纵向位置不变，系统地变化L/B、B/d发展成系列60。系列60在民用船中应用十分广泛，在船型研究中常被用作标准船型。陶德法比较适用于估算单螺旋桨运输船的有效功率。由于系列60船型的阻力性能甚好，因而用此法估算的有效功率对一般船来讲往往偏低一些。 系列60的船型系数范围为： Lbp/B 6.58.5；B/d 2.53.5； Cb 0.600.80； 航速范围： Vs / 0.41.0(相应Fr0.120.30)。其中Lbp为垂线间长，Lwl为水线长，计算方形系数用垂线间长

19、，Vs为航速，其单位为节(kn)。所用单位均为英制，即长度以英尺(ft)计，排水量为英吨，阻力用磅力(1bf)计算。 系列60阻力试验结果已绘制成每英吨排水量剩余阻力(磅力)Rr/图谱。并以B/d2.50，3.00和3.75分成三组。每组对一定的Vs /给出一张以L/B为纵坐标，方形系数Cb为横坐标的Rr/等值图谱。图7-3是该系列图谱的一张。 其次摩擦阻力按桑海公式计算，粗糙度补贴系数Cf取0.410-3。为方便起见，已绘成诺模图，根据航速和船长可直接得到每平方英尺湿面积的摩擦阻力Rf /S。为了确定摩擦阻力，该方法还给出了三种B/d情况下的该系列湿面积系数图谱。 陶德法估算阻力的具体方法是

20、；在给定航速时，分别由相应图谱求得S、Rr/和Rf /S值之后，则总阻力为：Rt(Rr / )(Rf/S ) S如设计船的浮心纵向位置不在系列60母型的最佳位置时，特别是两者相差较大时，则需另加修正。 由于在应用电子计算机进行阻力分析、计算及船型参数的优选中，使用回归公式具有很多优越性，所以我们重点介绍系列60的阻力回归公式。系列60的阻力回归公式是用线性回归分析方法处理系列60的试验图谱得到的，其总阻力系数Ct表示为下列函数形式：Ct = f (7-6) 对于标准船长为122m的船，其总阻力系数Ct122的回归多项式可表示为：Ct122 = A1 + A2(L/B) + A3(B/d) +

21、A4Cb + A5xb + A6(L/B)2 + A7(B/d)2+ A8(Cb)2 + A9(xb)2 + A10(L/B)(B/d) + A11(L/B)Cb + A12(L/B)xb + A13(B/d)Cb + A14(B/d)xb + A15Cbxb + A16xb(Cb)2 (7-7) 式中， Ct = 0.9492 (7-8)其中 Rt 总阻力（N）； 排水量（t）；Vs 航速（kn）；L 垂线间长（m）；Lwl 水线长（m）；Ct122 垂线间长为122m船的总阻力系数；xb 浮心纵向位置，以距船中%L计，舯前为正。为了便于确定回归系数，将（7-7）式化为下列线性形式： Y12

22、2 = a1x1 + a2x2 + a3x3+ a4x4 + a5x5 + a6x6+ a7x7 + a8x8 +a9x9 + a10x10+ a11x11 + a12x12 + a13x13+ a14x14 + a15x15 + a16x16 (7-9)式中，x1=1；x2 = 2(L/B-7)/3；x3 = 2(B /d -3)；x4 = 10(Cb-0.7)；x5 =( xb -0.515)/ 2.995；x6 = x22；x7 = x32；x8 = x42；x9 = x52；x10 = x2 x3；x11 = x2 x4；x12 = x2 x5；x13 = x3 x4；x14 = x3

23、 x5；x15 = x4 x5； x16 = x5 x42；Y122 =(Ct 122-17.3505)/ 8.3375。 对于不同的速长比Vs/，(7-9)式中的系数列于表7-2。船长为L的设计船的总阻力系数CtL可通过对标准船长的总阻力系数加以修正得到：CtL = Ct122 + SFC (7-10)表7-2 陶德法的回归方程系数式（7-9）中系数Vs /0.90550.99611.08661.17721.26781.35831.44891.53941.63a1a2a3a4a5a6a7a8a9a10a11a12a13a14a15a16- 0.8244+ 0.1906+ 0.1164- 0.

24、0519+ 0.0057+ 0.0072- 0.0052+ 0.1134+ 0.0670+ 0.0483- 0.1276+ 0.1125- 0.0481+ 0.0372- 0.0954- 0.0629- 0.8249+ 0.1865+ 0.1133+ 0.0060- 0.0109+ 0.0198- 0.0036+ 0.1109+ 0.0917+ 0.0510- 0.0745+ 0.0971- 0.0213+ 0.0206- 0.1924+ 0.0108- 0.8278+ 0.2050+ 0.1042+ 0.0832- 0.0451+ 0.0211+ 0.0067+ 0.0933+ 0.0708

25、+ 0.0400- 0.0729+ 0.1269+ 0.0232- 0.0105- 0.0855+ 0.0036- 0.7970+ 0.2332+ 0.1116+ 0.1075- 0.0165+ 0.0172+ 0.0068+ 0.1041+ 0.0826+ 0.0409- 0.0879+ 0.1882+ 0.0265- 0.0049- 0.0189+ 0.0581- 0.7562+ 0.2496+ 0.1221+ 0.1494+ 0.0472+ 0.0216+ 0.0064+ 0.1585+ 0.1428+ 0.0414- 0.0744+ 0.2115+ 0.0188+ 0.0135+ 0.0018+ 0.0884- 0.6619+ 0.2607+ 0.1298+ 0.2603+ 0.1491+ 0.0361