Project 1

Istanbul Technical University

EFFECTS OF CAVITATION & NOISE ON PROPELLER DESIGN (GEM 606E)

PROJECT-1 KIVANÇ ALİ ANIL LTJG-MSME

508052003

May 08, 2007

KIVANÇ ALİ ANIL 508052003

May 08, 2007

POWER PREDICTION AND PROPELLER OPTIMIZATION PROGRAM

A .POWER PREDICTION PROGRAM:

A MATLAB code “powerpredict.m” is written in order to estimate the Effective Power of a ship. The program uses the Holtrop&Mennen algorithm.

An example case is

demonstrated here and the results are compared to the results of the Power Prediction Program (PPP) created by M. PARSONS from the University of Michigan, Department of Naval Architecture and Marine Engineering. The figures plotted using the results of both programs showed that the “powerpredict.m” works accurately.

The “powerpredict.m” is a user friendly program. It uses dialogue windows for the entry of the input data. Moreover, the program uses menu and list windows for the entry of the multiple choice data. However, user can enter the input data by changing the default values of the program using MATLAB or a text editor, since those are the particular values for the illustrated case.

The “powerpredict.m” plots the resistance components and the total resistance with respect to the ship speed. It also plots the Effective Power in trial and service conditions and shows the Effective Power value for the corresponding trial speed of the ship. Furthermore, the “powerpredict.m” writes all input and output data into an Excel file.

-1-


May 08, 2007

-2-


May 08, 2007

-3-


May 08, 2007

-4-


May 08, 2007

-5-


May 08, 2007

-6-


May 08, 2007

-7-


May 08, 2007

-8-


May 08, 2007

V (knots) RF (N) k1*RF (N) RA (N) (RF*(1+k1)+ RA) (N) RAPP (N) RB (N) RTR (N) (RAPP + RB + RTR) (N) RW (N) 2 .2 .2 20 20 45 45 EE- 16 16 1 .5 .5 09 09 24 24 EE- 29 29 3 .6 .68 33 33 6E 6E- 30 30 9 .7 .77 24 24 5E 5E- 30 30 2 .8 .8 54 54 82 82E -2 -2 9 4 .1 .1 41 41 1E 1E -3 -3 1 3 .5 .5 69 69 43 43 EE- 46 46 1. 06 06 97 97 8E 8E -2 -2 9 1 .1 .11 12 12 EE-2 9 0 0.20 0.2020 2020 2020 202 2 49. 49.44 4466 6666 6642 42 12.0 12.067 6770 7039 393 3 8.0 8.089 8933 3328 2879 79 69.6 69.603 0370 7032 323 3 1.35 1.3567 6736 3688 884 4 0.2 0.268 6864 6498 9821 21 8.77 8.7715 1556 5655 553 3 10.3 10.396 9694 9432 326 6 6.89 6.8927 27EE-14 149 9 0.40 0.4040 4040 4040 404 4 175. 175.97 9765 6592 927 7 42.9 42.947 4795 9593 936 6 32.3 32.357 5733 3315 152 2 251. 251.28 2818 1883 836 6 4.82 4.8285 8514 1426 261 1 2.14 2.1451 5113 1335 358 8 34.7 34.751 5108 0815 151 1 41.7 41.724 2470 7091 913 3 2.32 2.3226 26EE-76 76 0.60 0.6060 6060 6060 606 6 370 370.9 .914 1444 4447 47 90.5 90.523 2350 5079 799 9 72. 72.80 8039 3995 9591 91 534. 534.24 2419 1948 486 6 10.1 10.177 7729 2949 494 4 7.2 7.216 1688 8879 7935 35 77.4 77.435 3585 8578 782 2 94.8 94.830 3004 0406 069 9 9.34 9.3400 005E 5E-5 -51 1 0.80 0.8080 8080 8080 808 8 630 630.3 .339 3901 0199 99 153. 153.83 8373 7360 609 9 129 129.4 .429 2932 3261 61 913. 913.60 6057 5706 069 9 17.2 17.295 9548 4874 743 3 17. 17.03 0313 1345 4508 08 136. 136.32 3231 3168 684 4 170. 170.65 6500 0000 009 9 1.70 1.7082 821E 1E-3 -37 7 1.01 1.0101 0101 0101 01 951. 951.71 7157 5741 419 9 232. 232.27 2709 0942 423 3 202. 202.23 2333 3322 22 1386 1386.2 .220 2000 006 6 26.1 26.113 1354 5470 705 5 33.0 33.076 7691 9135 358 8 210. 210.91 9102 0296 963 3 270. 270.10 1007 0756 569 9 2.4 2.499 9944 44EE-29 29 1.21 1.2121 2121 2121 212 2 1333 1333.1 .192 9233 335 5 325. 325.37 3721 2195 95 291. 291.21 2159 5983 836 6 1949 1949.7 .780 8051 514 4 36.5 36.580 8065 6503 033 3 56.7 56.765 6569 6990 907 7 300. 300.69 6945 4524 243 3 394. 394.04 0408 0873 737 7 8.9 8.986 8644 44EE-24 24 1.41 1.4141 4141 4141 414 4 177 1773. 3.32 3207 0781 81 432. 432.78 7877 7722 229 9 396 396.3 .377 7731 3111 11 2602 2602.4 .485 8581 815 5 48.6 48.657 5706 0659 598 8 89. 89.41 4179 7925 2502 02 405. 405.17 1731 3135 355 5 543. 543.24 2481 8126 265 5 9.76 9.7635 353E 3E-2 -20 0 1.61 1.6161 6161 6161 616 6 2270 2270.9 .919 1916 161 1 554. 554.22 2290 9072 729 9 517. 517.71 7173 7304 043 3 3342 3342.8 .865 6553 538 8 62.3 62.310 1036 3640 405 5 132. 132.24 2485 8507 074 4 523. 523.84 8434 3412 128 8 718. 718.40 4022 2284 843 3 1.14 1.1498 98EE-16 16 1.81 1.8181 8181 8181 818 8 282 2824. 4.99 9932 3229 29 689. 689.45 4535 3594 946 6 655 655.2 .235 3596 9632 32 4169 4169.6 .682 8278 787 7 77.5 77.513 1326 2640 405 5 186 186.3 .355 5599 9932 32 656. 656.20 2026 2639 391 1 920. 920.07 0718 1896 963 3 3.01 3.0162 625E 5E-1 -14 4 2.02 2.0202 0202 0202 02 3434 3434.6 .687 8775 754 4 838. 838.25 2525 2564 642 2 808. 808.93 9332 3287 879 9 5081 5081.8 .873 7360 606 6 94.2 94.242 4229 2997 979 9 252. 252.71 7140 4005 05 801. 801.74 7480 8097 974 4 1148 1148.7 .704 0440 402 2 2.73 2.7308 085E 5E-1 -12 2 2.22 2.2222 2222 2222 222 2 4099 4099.2 .254 5433 332 2 1000 1000.4 .443 4333 33 978. 978.80 8092 9278 784 4 6078 6078.5 .506 0694 94 112. 112.47 4769 6936 363 3 332. 332.16 1652 5249 492 2 959. 959.97 9770 7070 706 6 1404 1404.6 .619 1925 256 6 1.1 1.131 3113 13EE-10 10 2.42 2.4242 4242 4242 424 4 4818 4818.0 .029 2904 042 2 1175 1175.8 .863 6386 86 1164 1164.8 .863 6393 935 5 7158 7158.7 .756 5683 836 6 132. 132.19 1989 8956 566 6 425. 425.41 4180 8058 589 9 1130 1130.3 .386 8684 842 2 1688 1688.0 .003 0385 857 7 2.5 2.591 9148 48EE-09 09 2.62 2.6262 6262 6262 626 6 5590 5590.4 .416 1632 329 9 1364 1364.3 .368 6880 804 4 1367 1367.0 .097 9725 257 7 8321 8321.8 .882 8239 39 153. 153.39 3920 2019 197 7 533. 533.04 0453 5371 712 2 1312 1312.4 .474 7469 694 4 1998 1998.9 .912 1208 085 5 3.74 3.7498 982E 2E-0 -08 8 2.82 2.8282 8282 8282 828 8 641 6415. 5.87 8770 7012 12 1565 1565.8 .826 2657 573 3 158 1585. 5.50 5092 9244 44 9567 9567.2 .212 1282 829 9 176. 176.04 0413 1331 313 3 655 655.4 .485 8597 9725 25 1505 1505.7 .737 3790 909 9 2337 2337.2 .265 6521 213 3 3.77 3.7704 048E 8E-0 -07 7 3.03 3.0303 0303 0303 03 7293 7293.9 .919 1913 134 4 178 1780. 0.11 1171 7104 04 1820 1820.0 .099 9989 898 8 1089 10894. 4.13 1361 614 4 200. 200.13 1333 3392 929 9 793 793.0 .047 4779 7942 42 1709 1709.6 .673 7377 772 2 2702 2702.8 .854 5495 959 9 2.82 2.8273 731E 1E-0 -06 6 3.23 3.2323 2323 2323 232 2 822 8224. 4.09 0908 0833 33 2007 2007.1 .130 3011 116 6 207 2070. 0.86 8692 9217 17 1230 12302. 2.09 0901 017 7 225. 225.65 6558 5806 061 1 945 945.9 .913 1300 0034 34 1923 1923.7 .779 7956 564 4 3095 3095.3 .348 4837 373 3 1.66 1.6676 768E 8E-0 -05 5 3.43 3.4343 4343 4343 434 4 9205 9205.9 .974 7468 68 2246 2246.7 .763 6373 733 3 2337 2337.8 .817 1720 202 2 1379 13790. 0.55 5556 562 2 252. 252.59 5971 7117 173 3 1114 1114.1 .144 4460 608 8 2147 2147.5 .552 5256 569 9 3514 3514.2 .294 9429 294 4 8.0 8.062 6208 08EE-05 05 3.63 3.6363 6363 6363 636 6 1023 10239. 9.18 1831 311 1 2498 2498.9 .923 2337 37 2620 2620.9 .943 4385 853 3 1535 15359. 9.05 0503 033 3 280. 280.94 9466 6692 924 4 1297 1297.6 .694 9429 29 2380 2380.4 .490 9006 069 9 3959 3959.1 .131 3105 051 1 0.00 0.0003 0329 2986 867 7 3.83 3.8383 8383 8383 838 8 113 11323 23.3 .354 5469 69 2763 2763.5 .520 2081 815 5 292 2920. 0.24 2491 9169 69 1700 17007. 7.12 1246 467 7 310. 310.69 6946 4614 143 3 149 1496. 6.41 4111 1165 65 2622 2622.0 .089 8934 348 8 4429 4429.1 .195 9512 127 7 0.0 0.001 0117 1718 1863 63 4.04 4.0404 0404 0404 04 1245 12458. 8.15 1510 103 3 304 3040. 0.47 4734 3485 85 3235 3235.7 .733 3315 152 2 1873 18734. 4.35 3576 767 7 341. 341.83 8315 1598 981 1 171 1710. 0.05 0512 1207 07 2871 2871.8 .847 4768 688 8 4923 4923.7 .730 3049 493 3 0.0 0.003 0368 6895 9596 96 4.24 4.2424 2424 2424 242 2 1364 13643. 3.25 2542 421 1 3329 3329.7 .703 0378 787 7 3567 3567.3 .395 958 8 2054 20540. 0.35 3537 379 9 374. 374.34 3489 8920 205 5 1938 1938.2 .287 8705 053 3 3129 3129.2 .262 6237 373 3 5441 5441.8 .898 9834 346 6 0.01 0.0104 0468 6897 977 7 4.44 4.4444 4444 4444 444 4 148 14878 78.3 .364 6457 57 3631 3631.1 .138 3859 591 1 391 3915. 5.23 2371 7113 13 2242 22424. 4.74 7402 027 7 408. 408.23 2383 8359 595 5 218 2180. 0.71 7179 7978 78 3393 3393.8 .830 3068 685 5 5982 5982.7 .787 8702 022 2 0.0 0.027 2714 1407 0775 75 4.64 4.6464 6464 6464 646 6 161 16163 63.1 .198 9889 89 3944 3944.7 .708 0877 774 4 427 4279. 9.25 2570 7093 93 2438 24387. 7.16 1647 476 6 443. 443.49 4921 2143 439 9 243 2436. 6.87 8798 9801 01 3665 3665.0 .049 4990 907 7 6545 6545.4 .421 2185 852 2 0.0 0.065 6502 0260 6058 58 4.84 4.8484 8484 8484 848 8 1749 17497. 7.48 4887 878 8 4270 4270.3 .348 4883 832 2 4659 4659.4 .455 5573 738 8 2642 26427. 7.29 2933 335 5 480. 480.10 1029 2909 093 3 2706 2706.2 .254 5457 579 9 3942 3942.4 .417 1732 322 2 7128 7128.7 .774 7481 81 0.14 0.1453 5360 6089 893 3 5.05 5.0505 0505 0505 051 1 188 18880 80.9 .979 7928 28 4607 4607.9 .996 9654 542 2 505 5055. 5.83 8330 3049 49 2854 28544. 4.80 8088 887 7 518. 518.06 0636 3660 601 1 298 2988. 8.27 2799 9905 05 4225 4225.4 .430 3021 214 4 7731 7731.7 .773 7377 779 9 0.3 0.305 0563 6342 4294 94 5.25 5.2525 2525 2525 253 3 2031 20313. 3.42 4276 768 8 4957 4957.5 .592 9267 672 2 5468 5468.3 .389 8902 026 6 3073 30739. 9.40 4093 938 8 557. 557.36 3677 7737 37 3282 3282.3 .357 5773 734 4 4513 4513.5 .585 8586 865 5 8353 8353.3 .311 1133 336 6 0.60 0.6085 8589 8935 357 7 5.45 5.4545 4545 4545 455 5 217 21794 94.6 .602 0249 49 5319 5319.0 .080 8072 723 3 589 5897. 7.12 1236 3669 69 3301 33010. 0.80 8068 688 8 598. 598.00 0087 8788 881 1 358 3587. 7.86 8626 2616 16 4806 4806.3 .381 8155 558 8 8992 8992.2 .252 5296 962 2 1.1 1.154 5440 4031 3134 34 5.65 5.6565 6565 6565 657 7 233 23324 24.2 .282 8248 48 5692 5692.4 .406 0670 705 5 634 6342. 2.03 0369 6977 77 3535 35358. 8.72 7261 616 6 639. 639.98 9807 0743 433 3 390 3904. 4.14 1492 9296 96 5103 5103.3 .314 1457 576 6 9647 9647.4 .444 4461 616 6 2.0 2.096 9647 4753 5304 04 5.85 5.8585 8585 8585 859 9 2490 24902. 2.25 2559 59 6077 6077.5 .518 1893 936 6 6803 6803.1 .128 2895 951 1 3778 37782. 2.90 9037 379 9 683. 683.27 2777 7792 923 3 4230 4230.5 .559 5964 644 4 5403 5403.8 .882 8220 202 2 1031 10317. 7.71 7196 964 4 3.66 3.6611 1124 2489 894 4 6.06 6.0606 0606 0606 061 1 265 26528 28.3 .319 1972 72 6474 6474.3 .367 6786 867 7 728 7280. 0.39 3995 9591 91 4028 40283. 3.08 0871 718 8 727. 727.89 8943 4364 644 4 456 4566. 6.42 4288 8885 85 5707 5707.5 .581 8171 719 9 1100 11001. 1.90 9049 497 7 6.1 6.171 7132 3295 9571 71 6.26 6.2626 2626 2626 263 3 2820 28202. 2.27 2790 901 1 6882 6882.9 .905 0592 925 5 7773 7773.8 .848 4889 897 7 4285 42859. 9.03 0338 383 3 773. 773.82 8251 5111 111 1 4911 4911.0 .091 9114 149 9 6013 6013.9 .910 1041 41 1169 11698. 8.82 8266 667 7 10.0 10.074 7445 4575 75 6.46 6.4646 4646 4646 465 5 2992 29923. 3.94 9463 634 4 7303 7303.0 .087 8736 369 9 8283 8283.4 .476 7686 868 8 4551 45510. 0.51 5105 058 8 821. 821.06 0648 4889 897 7 5263 5263.8 .884 8434 34 6322 6322.3 .365 6555 557 7 1240 12407. 7.31 3147 479 9 15.9 15.973 7375 7529 295 5 6.66 6.6666 6666 6666 667 7 3169 31693. 3.14 1413 13 7734 7734.8 .868 6816 164 4 8809 8809.2 .283 8350 505 5 4823 48237. 7.29 2929 297 7 869. 869.60 6087 8749 497 7 5624 5624.1 .154 5436 365 5 6632 6632.4 .444 4444 445 5 1312 13126. 6.20 2075 756 6 24.6 24.663 6315 1578 789 9 6.86 6.8686 8686 8686 869 9 3350 33509. 9.68 6899 997 7 8178 8178.2 .205 0586 867 7 9351 9351.2 .268 6880 808 8 5103 51039. 9.16 1646 464 4 919. 919.45 4519 1919 193 3 5991 5991.2 .258 5873 739 9 6943 6943.6 .644 4435 355 5 1385 13854. 4.35 3550 501 1 37.1 37.164 6489 8947 47 7.07 7.0707 0707 0707 071 1 3537 35373. 3.42 4245 455 5 8633 8633.0 .059 5952 525 5 9909 9909.4 .432 3277 777 7 5391 53915. 5.91 9168 685 5 970. 970.58 5897 9794 942 2 6364 6364.5 .569 6963 63 7255 7255.4 .462 6257 571 1 1459 14590. 0.62 622 2 54.7 54.769 6910 1087 871 1 7.27 7.2727 2727 2727 273 3 372 37284 84.1 .182 8293 93 9099 9099.3 .389 8957 577 7 104 10483 83.7 .775 7541 41 5686 56867. 7.34 3479 792 2 1023 1023.0 .017 1792 927 7 674 6743. 3.47 4763 6364 64 7567 7567.3 .396 9637 375 5 1533 15333. 3.89 8906 067 7 79. 79.07 0747 4778 7848 48 7.47 7.4747 4747 4747 475 5 3924 39241. 1.80 8084 84 9577 9577.1 .157 5777 771 1 1107 11074. 4.29 2967 671 1 5989 59893. 3.26 2628 288 8 1076 1076.7 .732 3201 016 6 7127 7127.3 .387 8745 457 7 7878 7878.9 .943 4305 051 1 1608 16083. 3.06 0625 252 2 112. 112.03 0310 1049 491 1 7.67 7.6767 6767 6767 677 7 412 41246 46.1 .149 4925 25 1006 10066. 6.32 3270 709 9 116 11680 80.9 .996 9668 68 6299 62993. 3.47 4730 301 1 1131 1131.7 .727 2790 901 1 751 7515. 5.73 7322 2202 02 8189 8189.5 .599 9988 882 2 1683 16837. 7.05 0599 999 9 155 155.9 .978 7812 1278 78 7.87 7.8787 8787 8787 879 9 4329 43297. 7.05 0585 853 3 1056 10566. 6.86 8616 166 6 1230 12303. 3.87 8753 531 1 6616 66167. 7.79 7955 55 1188 1188.0 .001 0154 548 8 7907 7907.9 .961 6187 878 8 8498 8498.8 .864 6415 15 1759 17594. 4.82 8275 758 8 213. 213.68 6869 6900 001 1 8.08 8.0808 0808 0808 081 1 4539 45394. 4.39 3937 379 9 1107 11078. 8.72 7267 672 2 1294 12942. 2.93 9326 261 1 6941 69416. 6.05 0531 312 2 1245 1245.5 .549 4905 05 8303 8303.5 .550 506 6 8806 8806.2 .233 3313 138 8 1835 18355. 5.33 3327 279 9 288. 288.39 3964 6471 71 8.28 8.2828 2828 2828 283 3 4753 47538. 8.01 0167 678 8 1160 11601. 1.88 8885 853 3 1359 13598. 8.16 1685 857 7 7273 72738. 8.07 0738 389 9 1304 1304.3 .366 6661 612 2 8701 8701.9 .995 9588 88 9111 9111.2 .204 0413 13 1911 19117. 7.56 5666 662 2 383. 383.84 8489 8909 095 5 8.48 8.4848 4848 4848 485 5 4972 49727. 7.79 7932 328 8 1213 12136. 6.31 3143 434 4 1426 14269. 9.58 5832 32 7613 76133. 3.69 6908 082 2 1364 1364.4 .450 5055 552 2 9102 9102.8 .818 1891 918 8 9413 9413.2 .274 7440 408 8 1988 19880. 0.54 5438 388 8 504. 504.32 3205 0561 619 9 8.68 8.6868 6868 6868 687 7 5196 51963. 3.59 5928 285 5 1268 12681. 1.97 9723 23 1495 14957. 7.17 1764 649 9 7960 79602. 2.74 7416 164 4 1425 1425.7 .797 9729 291 1 9505 9505.5 .564 6466 667 7 9711 9711.9 .941 4125 255 5 2064 20643. 3.30 3032 321 1 654. 654.64 6494 9403 03 8.88 8.8888 8888 8888 889 9 542 54245 45.2 .288 8865 65 1323 13238. 8.83 8314 146 6 156 15660 60.9 .948 4845 45 8314 83145. 5.06 0685 857 7 1488 1488.4 .403 0334 349 9 990 9909. 9.80 8017 1712 12 1000 10006. 6.70 7019 195 5 2140 21404. 4.90 9070 702 2 840 840.2 .258 5800 0029 29 9.09 9.0909 0909 0909 091 1 565 56572 72.7 .757 5726 26 1380 13806. 6.86 8617 171 1 163 16380 80.8 .899 9908 08 8676 86760. 0.51 5180 805 5 1552 1552.2 .265 6533 339 9 103 10315 15.1 .121 2198 98 1029 10297. 7.05 0537 379 9 2216 22164. 4.44 4411 111 1 106 1067. 7.17 1717 1796 96 9.29 9.2929 2929 2929 293 3 5894 58945. 5.87 8785 852 2 1438 14386. 6.03 0337 371 1 1711 17117. 7.02 0283 837 7 9044 90448. 8.94 9406 06 1617 1617.3 .379 7996 964 4 1072 10721. 1.14 1403 033 3 1058 10582. 2.49 4940 404 4 2292 22921. 1.01 0143 434 4 1342 1342.0 .032 3245 45 9.49 9.4949 4949 4949 495 5 6136 61364. 4.53 5353 534 4 1497 14976. 6.31 3188 889 9 1786 17869. 9.33 3363 633 3 9421 94210. 0.19 1905 057 7 1683 1683.7 .744 4401 011 1 1112 11127. 7.49 4940 401 1 1086 10862. 2.52 52 2367 23673. 3.75 7580 802 2 1672 1672.1 .106 0622 227 7 9.69 9.6969 6969 6969 697 7 6382 63828. 8.61 6136 362 2 1557 15577. 7.68 6894 941 1 1863 18637. 7.82 8229 295 5 9804 98044. 4.12 1259 598 8 1751 1751.3 .354 5435 35 1153 11533. 3.84 8420 207 7 1113 11136. 6.62 6289 893 3 2442 24421. 1.82 8253 535 5 2065 2065.2 .287 8737 372 2 9.89 9.8989 8989 8989 899 9 663 66338 38.0 .002 0205 05 1619 16190. 0.11 1180 808 8 194 19422 22.4 .488 8824 24 1019 101950 50.6 .608 084 4 1820 1820.2 .207 0792 927 7 119 11939 39.8 .864 6464 64 1140 11404. 4.31 3181 814 4 2516 25164. 4.39 3907 071 1 253 2530. 0.09 0965 6569 69 10.1 10.101 0101 0101 01 6889 68892. 2.59 592 2 1681 16813. 3.57 5784 842 2 2022 20223. 3.33 3322 22 1059 105929 29.5 .502 026 6 1890 1890.3 .301 0176 761 1 1234 12345. 5.26 2622 225 5 1166 11665. 5.08 0848 489 9 2590 25900. 0.64 6489 89 3075 3075.6 .674 7423 234 4 10.3 10.303 0303 0303 03 7149 71492. 2.27 2774 744 4 174 17448 48.0 .044 4453 53 2104 21040. 0.35 3548 482 2 1099 109980 80.6 .676 768 8 1961 1961.6 .632 3294 942 2 127 12749 49.7 .755 5501 01 1191 11918. 8.42 4264 648 8 2662 26629. 9.81 8144 443 3 371 3711. 1.76 7697 9795 95 10.5 10.505 0505 0505 051 1 7413 74136. 6.95 9547 477 7 1809 18093. 3.49 4911 114 4 2187 21873. 3.55 5561 61 1141 114104 04.0 .002 02 2034 2034.1 .198 9862 628 8 1315 13153. 3.08 0818 185 5 1216 12163. 3.84 8401 018 8 2735 27351. 1.12 1206 066 6 4448 4448.7 .729 2964 641 1 10.7 10.707 0707 0707 071 1 768 76826 26.5 .522 2276 76 1874 18749. 9.89 8935 353 3 227 22722 22.9 .936 3606 06 1182 118299 99.3 .352 524 4 2107 2107.9 .996 9604 042 2 135 13554 54.9 .999 9972 72 1240 12400. 0.82 8232 329 9 2806 28063. 3.81 8190 905 5 529 5297. 7.47 4769 6912 12 10.9 10.909 0909 0909 091 1 7956 79560. 0.88 8824 245 5 1941 19417. 7.22 2275 756 6 2358 23588. 8.49 4946 468 8 1225 122566 66.6 .604 047 7 2183 2183.0 .022 2246 467 7 1395 13955. 5.28 2827 275 5 1262 12628. 8.87 8730 307 7 2876 28767. 7.17 1782 829 9 6269 6269.4 .470 7055 55 11.1 11.111 1111 1111 111 1 8233 82339. 9.93 9370 702 2 2009 20095. 5.46 4696 96 2447 24470. 0.23 2319 196 6 1269 126905 05.6 .638 386 6 2259 2259.2 .275 7524 248 8 1435 14353. 3.72 7214 147 7 1284 12847. 7.48 4868 682 2 2946 29460. 0.48 4835 353 3 7376 7376.6 .672 7247 478 8 11.3 11.313 1313 1313 131 1 851 85163 63.5 .591 9174 74 2078 20784. 4.59 5965 652 2 253 25368 68.1 .147 4791 91 1313 131316 16.3 .336 362 2 2336 2336.7 .751 5178 786 6 147 14750 50.1 .121 2198 98 1305 13056. 6.16 1618 181 1 3014 30143. 3.03 0355 558 8 863 8631. 1.58 5888 8883 83 11.5 11.515 1515 1515 152 2 8803 88031. 1.75 7538 389 9 2148 21484. 4.58 5857 571 1 2628 26282. 2.24 2425 252 2 1357 135798 98.5 .582 821 1 2415 2415.4 .449 4953 536 6 1514 15144. 4.30 3052 522 2 1325 13254. 4.39 3953 534 4 3081 30814. 4.15 1500 009 9 1004 10047. 7.33 3393 93 11.7 11.717 1717 1717 172 2 9094 90944. 4.33 3326 266 6 2219 22195. 5.41 4149 499 9 2721 27212. 2.51 5158 58 1403 140352 52.2 .263 635 5 2495 2495.3 .366 6600 006 6 1553 15536. 6.10 1061 614 4 1344 13441. 1.68 6846 468 8 3147 31473. 3.15 1568 683 3 1163 11637. 7.51 5189 893 3 11.9 11.919 1919 1919 192 2 939 93901 01.2 .239 3908 08 2291 22917. 7.06 0626 265 5 281 28158 58.9 .967 6775 75 1449 144977 77.2 .269 695 5 2576 2576.4 .498 9875 756 6 159 15925 25.3 .373 7301 01 1361 13617. 7.52 5271 711 1 3211 32119. 9.39 3988 888 8 134 13415 15.7 .733 3333 33 12.1 12.121 2121 2121 212 2 969 96902 02.3 .385 8598 98 2364 23649. 9.50 5074 742 2 291 29121 21.5 .598 9836 36 1496 149673 73.4 .491 918 8 2658 2658.8 .845 4539 393 3 163 16311 11.9 .966 6666 66 1378 13781. 1.41 4199 993 3 3275 32752. 2.23 2319 198 8 153 15395 95.2 .218 1845 45 12.3 12.323 2323 2323 232 2 9994 99947. 7.68 6878 789 9 2439 24392. 2.72 7284 844 4 3010 30100. 0.40 4076 764 4 1544 154440 40.8 .824 24 2742 2742.4 .403 0357 572 2 1669 16695. 5.75 7597 979 9 1393 13932. 2.86 8604 041 1 3337 33371. 1.02 0237 377 7 1758 17589. 9.31 3114 147 7 12.5 12.525 2525 2525 253 3 1030 103037 37.0 .061 61 2514 25146. 6.70 7052 525 5 3109 31095. 5.39 3955 559 9 1592 159279 79.1 .161 618 8 2827 2827.1 .170 7099 994 4 1707 17076. 6.63 6363 631 1 1407 14071. 1.34 3458 583 3 3397 33975. 5.15 1531 314 4 2001 20012. 2.96 9647 473 3 12.7 12.727 2727 2727 273 3 1061 106170 70.4 .423 231 1 2591 25911. 1.41 4177 779 9 3210 32106. 6.56 5622 22 1641 164188 88.4 .403 031 1 2913 2913.1 .145 4540 402 2 1745 17454. 4.49 4907 071 1 1419 14196. 6.37 3734 348 8 3456 34564. 4.00 0095 959 9 2268 22684. 4.18 1888 88 12.9 12.929 2929 2929 293 3 109 10934 347. 7.69 6934 34 2668 26686. 6.84 8463 637 7 331 33133 33.9 .907 0747 47 1691 169168 68.4 .447 473 3 3000 3000.3 .324 2458 583 3 178 17829 29.2 .227 2739 39 1430 14307. 7.44 4406 064 4 3513 35136. 6.99 9926 262 2 256 25623 23.5 .535 3535 35 13.1 13.131 3131 3131 313 3 1125 112568 68.7 .792 929 9 2747 27472. 2.97 9716 165 5 3417 34177. 7.43 4314 141 1 1742 174219 19.1 .195 959 9 3088 3088.7 .706 0636 364 4 1820 18200. 0.76 7601 014 4 1440 14404. 4.04 0446 46 3569 35693. 3.51 5111 111 1 2885 28850. 0.61 6131 314 4 13.3 13.333 3333 3333 333 3 115 11583 833. 3.64 6437 37 2826 28269. 9.77 7746 466 6 352 35237 37.1 .134 3402 02 1793 179340 40.5 .552 524 4 3178 3178.2 .288 8861 612 2 185 18569 69.0 .011 1153 53 1448 14485. 5.68 6826 263 3 3623 36232. 2.98 9827 278 8 323 32378 78.9 .901 0157 57 13.5 13.535 3535 3535 354 4 119 11914 142. 2.16 1694 94 2907 29077. 7.23 2367 677 7 363 36313 13.0 .015 1529 29 1845 184532 32.4 .421 215 5 3269 3269.0 .069 6922 229 9 189 18933 33.9 .912 1239 39 1455 14551. 1.85 8520 203 3 3675 36754. 4.83 8336 365 5 362 36212 12.1 .198 9859 59 13.7 13.737 3737 3737 374 4 1224 122494 94.2 .295 951 1 2989 29895. 5.33 3396 965 5 3740 37405. 5.07 0752 523 3 1897 189794 94.7 .71 1 3361 3361.0 .046 4615 159 9 1929 19295. 5.40 4013 131 1 1460 14602. 2.05 0500 006 6 3725 37258. 8.49 4975 753 3 4034 40346. 6.23 2358 587 7 13.9 13.939 3939 3939 394 4 1258 125889 89.9 .947 47 3072 30724. 4.06 0653 532 2 3851 38513. 3.31 3138 384 4 1951 195127 27.3 .326 262 2 3454 3454.2 .217 1737 378 8 1965 19653. 3.42 4241 415 5 1463 14635. 5.77 7740 403 3 3774 37743. 3.41 4155 555 5 4477 44776. 6.74 7453 532 2 14.1 14.141 4141 4141 414 4 1293 129329 29.0 .052 527 7 3156 31563. 3.39 3960 609 9 3963 39637. 7.73 7311 111 1 2005 200530 30.1 .179 798 8 3548 3548.5 .580 8089 897 7 2000 20007. 7.93 9335 358 8 1465 14652. 2.52 5212 12 3820 38209. 9.03 0356 567 7 4951 49511. 1.58 5830 301 1 14.3 14.343 4343 4343 434 4 132 13281 811. 1.54 5408 08 3241 32413. 3.31 3145 456 6 407 40778 78.3 .327 2704 04 2060 206003 03.1 .182 824 4 3644 3644.1 .134 3476 761 1 203 20358 58.8 .888 8866 66 1465 14651. 1.78 7888 886 6 3865 38654. 4.81 8122 228 8 545 54581 81.5 .553 5381 81 14.5 14.545 4545 4545 455 5 1363 136337 37.3 .341 413 3 3327 33273. 3.80 8036 365 5 4193 41935. 5.10 1016 164 4 2115 211546 46.2 .246 466 6 3740 3740.8 .877 7704 048 8 2070 20706. 6.25 2544 441 1 1463 14633. 3.07 0743 43 3908 39080. 0.20 2057 576 6 6004 60044. 4.23 2368 687 7

-9-

KIVANÇ ALİ ANIL 508052003 V (knots) 0.00 0.0000 0000 0000 00 0.20 0.2020 2020 2020 20 0.40 0.4040 4040 4040 40 0.60 0.6060 6060 6061 61 0.80 0.8080 8080 8081 81 1.01 1.0101 0101 0101 01 1.21 1.2121 2121 2121 21 1.41 1.4141 4141 4141 41 1.61 1.6161 6161 6162 62 1.81 1.8181 8181 8182 82 2.02 2.0202 0202 0202 02 2.2222 2.2222222 2222 2 2.4242 2.4242424 4242 2 2.6262 2.6262626 6263 3 2.8282 2.8282828 8283 3 3.0303 3.0303030 0303 3 3.2323 3.2323232 2323 3 3.4343 3.4343434 4343 3 3.6363 3.6363636 6364 4 3.8383 3.8383838 8384 4 4.0404 4.0404040 0404 4 4.2424 4.2424242 2424 4 4.4444 4.4444444 4444 4 4.6464 4.6464646 6465 5 4.8484 4.8484848 8485 5 5.0505 5.0505050 0505 5 5.2525 5.2525252 2525 5 5.4545 5.4545454 4545 5 5.6565 5.6565656 6566 6 5.8585 5.8585858 8586 6 6.0606 6.0606060 0606 6 6.2626 6.2626262 2626 6 6.4646 6.4646464 4646 6 6.6666 6.6666666 6667 7 6.8686 6.8686868 8687 7 7.0707 7.0707070 0707 7 7.2727 7.2727272 2727 7 7.4747 7.4747474 4747 7 7.6767 7.6767676 6768 8 7.8787 7.8787878 8788 8 8.0808 8.0808080 0808 8 8.2828 8.2828282 2828 8 8.4848 8.4848484 4848 8 8.6868 8.6868686 6869 9 8.8888 8.8888888 8889 9 9.0909 9.0909090 0909 9 9.2929 9.2929292 2929 9 9.4949 9.4949494 4949 9 9.6969697 9.69696970 0 9.8989899 9.89898990 0

RT (N) 0.00 0.0000 0000 0000 00 80.0 80.000 0064 6464 649 9 293. 293.00 0065 6592 9268 68 629. 629.07 0719 1989 8926 26 1084 1084.2 .255 5570 7078 781 1 1656 1656.3 .320 2076 7630 300 0 2343 2343.8 .821 2138 3877 771 1 3145 3145.7 .733 3394 9414 145 5 4061 4061.2 .267 6782 8223 231 1 5089 5089.7 .754 5468 6834 343 3 6230 6230.5 .578 7800 0083 838 8 7483.1 7483.1261 261965 9650 0 8846.7 8846.7606 606934 9346 6 10320. 10320.794 794474 47415 15 11904. 11904.478 478042 04281 81 13596. 13596.991 991097 09704 04 15397. 15397.438 438556 55616 16 17304. 17304.849 849990 99020 20 19318. 19318.181 181711 71197 97 21436. 21436.320 320970 97094 94 23658. 23658.091 091847 84768 68 25982. 25982.262 262607 60777 77 28407. 28407.554 554437 43720 20 30932. 30932.651 651640 64031 31 33556. 33556.213 213521 52167 67 36276. 36276.888 888279 27996 96 39093. 39093.329 329301 30143 43 42004. 42004.214 214244 24433 33 45008. 45008.267 267254 25444 44 48104. 48104.284 284550 55000 00 51291. 51291.163 163474 47493 93 54567. 54567.934 934954 95485 85 57933. 57933.799 799118 11876 76 61388. 61388.163 163683 68358 58 64930. 64930.684 684551 55144 44 68561. 68561.307 307950 95077 77 72280. 72280.313 313366 36651 51 76088. 76088.356 356454 45476 76 79986. 79986.511 511123 12325 25 83976. 83976.309 309976 97697 97 88059. 88059.782 782375 37543 43 92239. 92239.489 489417 41758 58 96518. 96518.555 555258 25846 46 100900 100900.69 .69426 426017 017 105390 105390.23 .23358 358635 635 109992 109992.13 .13095 095695 695 114711 114711.98 .98738 738517 517 119556 119556.05 .05481 481468 468 124531.23 124531.238707 870799 99 129645.09 129645.095651 565176 76

May 08, 2007

PE (kW) trial PE (kW) service 0.00 0.0000 0000 0000 00 0.00 0.0000 0000 0000 00 0.00 0.0083 8313 1360 60 0.00 0.0099 9976 7632 32 0.06 0.0608 0898 9802 02 0.07 0.0730 3077 7762 62 0.19 0.1961 6117 1796 96 0.23 0.2353 5341 4155 55 0.45 0.4506 0699 9991 91 0.54 0.5408 0839 3989 89 0.86 0.8606 0617 1758 58 1.03 1.0327 2741 4109 09 1.46 1.4614 1408 0815 15 1.75 1.7536 3689 8978 78 2.28 2.2883 8314 1490 90 2.74 2.7459 5977 7789 89 3.37 3.3763 6349 4936 36 4.05 4.0516 1619 1923 23 4.76 4.7603 0308 0874 74 5.71 5.7123 2370 7049 49 6.47 6.4747 4766 6632 32 7.76 7.7697 9719 1958 58 8.5540 8.5540447 4470 0 10.264 10.264853 85364 64 11.032 11.032178 17867 67 13.238 13.238614 61440 40 13.942 13.942872 87208 08 16.731 16.731446 44650 50 17.319 17.319452 45234 34 20.783 20.783342 34281 81 21.194 21.194824 82491 91 25.433 25.433789 78989 89 25.601 25.601429 42996 96 30.721 30.721715 71595 95 30.571 30.571202 20246 46 36.685 36.685442 44296 96 36.135 36.135536 53699 99 43.362 43.362644 64439 39 42.325 42.325257 25791 91 50.790 50.790309 30949 49 49.170 49.170595 59574 74 59.004 59.004714 71489 89 56.701 56.701170 17042 42 68.041 68.041404 40451 51 64.945 64.945982 98223 23 77.935 77.935178 17868 68 73.933 73.933411 41173 73 88.720 88.720094 09408 08 83.691 83.691230 23023 23 100.42 100.42947 947628 628 94.246 94.246622 62288 88 113.09 113.09594 594746 746 105.62 105.62622 622695 695 126.75 126.75147 147234 234 117.85 117.85618 618804 804 141.42 141.42742 742565 565 130.96 130.96223 223736 736 157.15 157.15468 468483 483 144.96 144.96979 979297 297 173.96 173.96375 375156 156 159.90 159.90408 408783 783 191.88 191.88490 490539 539 175.79 175.79032 032686 686 210.94 210.94839 839223 223 192.65 192.65387 387486 486 231.18 231.18464 464983 983 210.52 210.52047 047599 599 252.62 252.62457 457119 119 229.41 229.41650 650516 516 275.29 275.29980 980619 619 249.36 249.36925 925017 017 299.24 299.24310 310021 021 270.40 270.40722 722324 324 324.48 324.48866 866789 789 292.56 292.56049 049914 914 351.07 351.07259 259897 897 315.86 315.86107 107681 681 379.03 379.03329 329218 218 340.34 340.34326 326065 065 408.41 408.41191 191278 278 366.04 366.04405 405700 700 439.25 439.25286 286840 840 393.00 393.00358 358134 134 471.60 471.60429 429760 760 421.26 421.26547 547124 124 505.51 505.51856 856549 549 450.87 450.87730 730030 030 541.05 541.05276 276036 036 481.89 481.89098 098806 806 578.26 578.26918 918567 567 514.36 514.36320 320149 149 617.23 617.23584 584179 179 548.35 548.35574 574350 350 658.02 658.02689 689219 219 583.93 583.93592 592445 445 700.72 700.72310 310934 934 621.17 621.176913 691337 37 745.41229 745.41229605 605 660.15 660.158065 806524 24 792.18967 792.18967829 829

V (knots) 10.1 10.101 0101 0101 010 0 10.3 10.303 0303 0303 030 0 10.5 10.505 0505 0505 051 1 10.7 10.707 0707 0707 071 1 10.9 10.909 0909 0909 091 1 11.1 11.111 1111 1111 111 1 11.3 11.313 1313 1313 131 1 11.5 11.515 1515 1515 152 2 11.7 11.717 1717 1717 172 2 11.9 11.919 1919 1919 192 2 12.1 12.121 2121 2121 212 2 12.323 12.323232 23232 32 12.525 12.525252 25253 53 12.727 12.727272 27273 73 12.929 12.929292 29293 93 13.131 13.131313 31313 13 13.333 13.333333 33333 33 13.535 13.535353 35354 54 13.737 13.737373 37374 74 13.939 13.939393 39394 94 14.141 14.141414 41414 14 14.343 14.343434 43434 34 14.545 14.545454 45455 55 14.747 14.747474 47475 75 14.949 14.949494 49495 95 15.151 15.151515 51515 15 15.353 15.353535 53535 35 15.555 15.555555 55556 56 15.757 15.757575 57576 76 15.959 15.959595 59596 96 16.161 16.161616 61616 16 16.363 16.363636 63636 36 16.565 16.565656 65657 57 16.767 16.767676 67677 77 16.969 16.969696 69697 97 17.171 17.171717 71717 17 17.373 17.373737 73737 37 17.575 17.575757 75758 58 17.777 17.777777 77778 78 17.979 17.979797 79798 98 18.181 18.181818 81818 18 18.383 18.383838 83838 38 18.585 18.585858 85859 59 18.787 18.787878 87879 79 18.989 18.989898 89899 99 19.191 19.191919 91919 19 19.393 19.393939 93939 39 19.595 19.595959 95960 60 19.7979798 19.79797980 0 20.0000000 20.00000000 0

- 10 -

RT (N) 1349 134905 05.8 .825 2575 7544 445 5 1403 140322 22.2 .261 6101 0105 050 0 1459 145903 03.8 .852 5231 3139 398 8 1516 151660 60.6 .648 4831 3123 231 1 1576 157603 03.2 .253 5352 5220 204 4 1637 163742 42.7 .794 9458 5833 332 2 1700 170090 90.9 .960 6063 6376 765 5 1766 176660 60.0 .071 7152 5219 192 2 1834 183462 62.9 .939 3921 2103 033 3 1905 190512 12.4 .401 0169 6982 822 2 1978 197820 20.9 .942 4220 2030 307 7 205401 205401.15 .15920 920858 858 213267 213267.27 .27969 969445 445 221436 221436.60 .60144 144835 835 229928 229928.97 .97523 523875 875 238763 238763.32 .32019 019890 890 247952 247952.43 .43671 671436 436 257499 257499.45 .45373 373760 760 267399 267399.44 .44341 341235 235 277647 277647.48 .48704 704739 739 288250 288250.79 .79853 853081 081 299239 299239.54 .54845 845397 397 310670 310670.68 .68918 918174 174 322621 322621.58 .58623 623898 898 335174 335174.06 .06620 620962 962 348392 348392.68 .68946 946301 301 362302 362302.77 .77618 618543 543 376874 376874.14 .14206 206000 000 392015 392015.92 .92255 255249 249 407586 407586.06 .06599 599273 273 423415 423415.73 .73913 913979 979 439344 439344.49 .49062 062506 506 455258 455258.16 .16557 557759 759 471120 471120.21 .21859 859629 629 486989 486989.15 .15024 024416 416 503019 503019.37 .37150 150444 444 519447 519447.66 .66401 401397 397 536570 536570.47 .47106 106153 153 554717 554717.72 .72129 129878 878 574227 574227.33 .33540 540239 239 595422 595422.25 .25836 836358 358 618589 618589.93 .93226 226113 113 643963 643963.16 .16850 850906 906 671701 671701.45 .45312 312663 663 701872 701872.55 .55653 653574 574 734435 734435.52 .52560 560055 055 769227 769227.29 .29830 830339 339 805955 805955.92 .92500 500724 724 844203.41 844203.4118757 187570 0 883440.37 883440.3765322 653229 9

PE (kW) trial PE (kW) service 700 700.9 .965 6521 2198 988 8 841. 841.15 1582 8263 6386 86 743. 743.69 6909 0974 7460 60 892. 892.42 4291 9169 6952 52 788. 788.43 4349 4942 4238 38 946. 946.12 1219 1930 3085 85 835 835.3 .303 0395 9569 698 8 1002 1002.3 .364 6474 7483 838 8 884.4 884.412 1214 1484 849 9 1061 1061.2 .294 9457 5781 819 9 935.8 935.881 8103 0392 926 6 1123 1123.0 .057 5724 2471 712 2 989.8 989.840 4005 0502 020 0 1187 1187.8 .808 0806 0602 025 5 1046 1046.4 .427 2719 1969 699 9 1255 1255.7 .712 1263 6363 638 8 1105 1105.7 .788 8858 5826 261 1 1326 1326.9 .946 4629 2991 913 3 1168 1168.0 .075 7579 7952 527 7 1401 1401.6 .690 9095 9543 432 2 1233 1233.4 .443 4354 5475 751 1 1480 1480.1 .132 3225 2570 701 1 1302.0 1302.0524 524715 7154 4 1562.4 1562.4629 629658 6584 4 1374.0 1374.0789 789288 2886 6 1648.8 1648.8947 947146 1463 3 1449.7 1449.7252 252990 9908 8 1739.6 1739.6703 703589 5890 0 1529.2 1529.2181 181315 3156 6 1835.0 1835.0617 617578 5787 7 1612.7 1612.7859 859341 3418 8 1935.3 1935.3431 431210 2101 1 1700.6 1700.6231 231126 1261 1 2040.7 2040.7477 477351 3513 3 1792.8 1792.8620 620551 5519 9 2151.4 2151.4344 344662 6622 2 1889.5 1889.5795 795173 1738 8 2267.4 2267.4954 954208 2085 5 1990.8 1990.8502 502719 7197 7 2389.0 2389.0203 203263 2637 7 2096.8 2096.8353 353037 0374 4 2516.2 2516.2023 023644 6448 8 2207.8 2207.8679 679766 7666 6 2649.4 2649.4415 415719 7199 9 2324.4 2324.4945 945820 8204 4 2789.3 2789.3934 934984 9845 5 2447.4 2447.4399 399412 4125 5 2936.9 2936.9279 279295 2950 0 2577.4 2577.4953 953403 4035 5 3092.9 3092.9944 944084 0841 1 2715.3 2715.3515 515069 0697 7 3258.4 3258.4218 218083 0836 6 2861.4 2861.4160 160915 9158 8 3433.6 3433.6993 993098 0989 9 3015.6 3015.6631 631349 3495 5 3618.7 3618.7957 957619 6195 5 3177.5 3177.5622 622755 7551 1 3813.0 3813.0747 747306 3061 1 3346.1 3346.1251 251546 5462 2 4015.3 4015.3501 501855 8555 5 3520.0 3520.0817 817165 1658 8 4224.0 4224.0980 980599 5990 0 3698.1 3698.1622 622796 7963 3 4437.7 4437.7947 947355 3556 6 3879.4 3879.4249 249758 7588 8 4655.3 4655.3099 099710 7105 5 4063.5 4063.5498 498903 9031 1 4876.2 4876.2598 598683 6837 7 4251.0 4251.0315 315932 9321 1 5101.2 5101.2379 379118 1185 5 4443.2 4443.2361 361615 6155 5 5331.8 5331.8833 833938 9386 6 4642.3 4642.3300 300080 0802 2 5570.7 5570.7960 960096 0963 3 4851.1 4851.1173 173691 6916 6 5821.3 5821.3408 408429 4299 9 5072.8 5072.8319 319259 2597 7 6087.3 6087.3983 983111 1117 7 5310.9 5310.9184 184198 1989 9 6373.1 6373.1021 021038 0387 7 5568.8 5568.8219 219945 9459 9 6682.5 6682.5863 863935 9350 0 5849.7 5849.7862 862959 9598 8 7019.7 7019.7435 435551 5518 8 6156.6 6156.6521 521529 5294 4 7387.9 7387.9825 825835 8353 3 6491.6 6491.6485 485164 1645 5 7789.9 7789.9782 782197 1974 4 6856.1 6856.1747 747171 1711 1 8227.4 8227.4096 096605 6053 3 7250.5 7250.5849 849020 0203 3 8700.7 8700.7018 018824 8244 4 7673.9 7673.9980 980072 0722 2 9208.7 9208.7976 976086 0866 6 8124.1 8124.1659 659795 7958 8 9748.9 9748.9991 991754 7549 9 8597.4357 8597.43576500 6500 10316.9229 10316.92291800 1800 9088.8345 9088.83459376 9376 10906.6015 10906.60151252 1252


May 08, 2007

- 11 -


May 08, 2007

- 12 -


May 08, 2007

- 13 -


May 08, 2007

Comparison of “powerpredict.m” and “PPP” 600000

500000

400000

) N ( F 300000 R

powerpredict.m PPP Result

200000

100000

0 0

5

10

15

20

25

30

35

V (knots)

140000

120000

100000

) N ( F R * 1 k

80000 powerpredict.m PPP Result 60000

40000

20000

0 0

5

10

15

20

V (knots)

- 14 -

25

30

35


May 08, 2007

200000

180000

160000

140000

120000

) N ( 100000 A R


80000

60000

40000

20000

0 0

5

10

15

20

25

30

35

V (knots)

16000

14000

12000

10000

) N ( P 8000 P A R


6000

4000

2000

0 0

5

10

15

20

V (knots)

- 15 -

25

30

35


May 08, 2007

40000

35000

30000

25000

) N ( 20000 B R


15000

10000

5000

0 0

5

10

15

20

25

30

35

V (knots)

16000

14000

12000

10000

) N ( R 8000 T R


6000

4000

2000

0 0

5

10

15

20

V (knots)

- 16 -

25

30

35


May 08, 2007

3500000

3000000

2500000

2000000 ) N ( W R


1500000

1000000

500000

0 0

5

10

15

20

25

30

35

V (knots)

4500000

4000000

3500000

3000000

2500000 ) N ( T R 2000000


1500000

1000000

500000

0 0

5

10

15

20

V (knots)

- 17 -

25

30

35


May 08, 2007

70000

60000

50000

40000 ) W k ( E P 30000


20000

10000

0 0

5

10

15

20

V (knots)

- 18 -

25

30

35


May 08, 2007

B .PROPELLER OPTIMIZATION PROGRAM:

A MATLAB code “propeller.m” is written for the basic design of a propeller using K T, K Q and open water efficiency curves for the Wageningen B-Series Propellers. The program uses the same algorithm of propopt.m and fcn.m from http://web.nps.navy.mil/ ~me/tsse/TS4001/docs_software.htm. User should run “ powerpredict.m” before running this program. There are four design options: 2

1. KT/J case – Known thrust (T), (T), propeller diameter (D) and advanced Nopt is required. 3 2. K Q/J case – Known power (PD), propeller diameter (D) and advanced Nopt is required. 4 3. K T/J case – Known thrust thrust (T), prop. rotation rotation rate (N) and advanced Dopt is required. 5 4. K Q/J case – Known power (PD), prop. rotation rate (N) and advanced Dopt is required.

velocity (V (VA), velocity (VA), velocity (VA), velocity (VA),

The program requires that function fcn.m is located in the same directory as this main program. “It uses a Lagrange multiplier method and utilizes the standard Matlab function fzero for the optimization problem”. The “Prediction of Performance in Service” part is added to the program for the prediction of the ship speed and propeller rate of rotation in service condition with the engine developing 85% of MCR. Blade Area Ratio (BAR) (BAR) is estimated using the Keller’s Method. Method. Besides, using Burrill’s method, method, cavitation control is performed

The “propeller.m” is a user friendly program similar to the “ powerpredict.m”. It uses dialogue windows for the entry of the input data. Moreover, the program uses menu and list windows for the entry of the multiple choice data. However, user can enter the input data by changing the default values of the program using MATLAB or a text editor, since those are the particular values for the illustrated case.

- 19 -


May 08, 2007

The “propeller.m” plots the Wageningen B-Series K T, K Q and open water efficiency curves for the selected number of blades (Z) and BAR. User can see the optimization steps on these curves by pressing the any key. It also plots the Effective Power in trial and service conditions and shows the predicted service speed of the ship for the corresponding Effective Power in service condition. The program plots the Burrill’s Cavitation Diagram and shows the thrust loading coefficient for the calculated local cavitation number at 0.7R. Furthermore, the “propeller.m” writes all input and output data into an Excel file. An example case is demonstrated here:

- 20 -


May 08, 2007

- 21 -


May 08, 2007

- 22 -


May 08, 2007

- 23 -


May 08, 2007

- 24 -


May 08, 2007

- 25 -


May 08, 2007

- 26 -


May 08, 2007

- 27 -


May 08, 2007

- 28 -


May 08, 2007

- 29 -


May 08, 2007

PROPELLER DESIGN Number of propellers - INPUT Number of blades - INPUT Wake Fraction - INPUT OR ESTIMATED Thrust Deduction Fraction - INPUT OR ESTIMATED Propeller Diameter in meters(initial) - INPUT OR ESTIMATED Open Water Diameter - CALCULATED AE/A0 ratio (BAR) - INPUT OR ESTIMATED USING KELLER Trial Speed (knots) from powerpredict.m Total Ship Resistance (N) for the Trial Speed from powerpredict.m Thrust for each propeller (Newton) - CALCULATED Advance Velocity (m/s) - CALCULATED Hull efficiency - CALCULATED TYPE=2 : KQ/J^3 case - Known power (PD), diameter (D) and ad vanced velocity (VA)- Nopt (required) Relative-rotative Efficiency - INPUT Maximum continous power at - INPUT Shaft transmission efficiency - INPUT Propulsive efficiency - CONVERGED Effective Power for the Trial Speed (kW) from powerpredict.m Delivered power (kW) - CALCULATED KQ/J^3 - CALCULATED Maximum Open Water Efficiency - OUTPUT Thrust Coefficient - OUTPUT Torque Coefficient - OUTPUT Pitch to Diameter - OUTPUT Advance Coefficient - OUTPUT Rotation Rate (rpm), N - OUTPUT Propeller Pitch (m) - OUTPUT Brake power (kW) - OUTPUT Installed maximum continous power (kW) - OUTPUT Propulsive efficiency for the service condition - CONVERGED Effective Power for the Service Condition (kW) - CALCULATED Ship speed (knots) in service condition - OUTPUT Propeller rate of rotation (rpm) in service condition - OUTPUT BURRILL CAVITATION CONTROL CONTROL :Suggested upper limit (1943) (for merchant ship propellers) Calculated BAR (0.51334) is LESS than the selected BAR (0.54454). THE DESIGN STAGE IS COMPLETED! COMPLETED!

- 30 -

1 5 0.3695 0.1869 4.355 4.584210526 0.54454 16 410736.7383 333970.0419 5.1892672 1.28961142 1 0.85 0.98 0.70405 3380.527651 4801.544849 0.34999 0.546569975 0.227196572 0.033771472 0.9 0.510474242 140.0539273 3.9195 4929.857795 5799.8327 0.7 3357.205235 15.26537031 137.4400086

powerpredict.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... ..

Page 1 May May 5, 2007 2007

% Kivanc Ali ANIL, 2007 % 508052003 %---------------------------------------------------------------------clear, clc, close all %---------------------------------------------------------------------% knt = 0.5144 m/sec % mile = knt.hr % metric ton = 1000 kg %---------------------------------------------------------------------% default values def0 = {'1.18831e-006' '1.18831e-006', ,... % Seawater Viscosity (m^2/sec) '1025', '1025' ,... % Seawater Density (kg/m^3) '9.8067'}; '9.8067' }; % g (m/sec^2) def1 = {'115.38' '115.38', ,... % Length (m), LWL '17.20', '17.20' ,... % Beam (m), B '6.70', '6.70' ,... % Draft (m),T '6.70', '6.70' ,... % Molded draft at the FP (m),TF '6.70', '6.70' ,... % Molded draft at the AP (m),TA '.783', '.783' ,... % Prismatic Coefficient, CP '.996', '.996' ,... % Midship Section Coefficient-CM '2.385', '2.385' ,... % LCB from midships as % of Length '15', '15' ,... % Service Speed (knt), Vs '1', '1' ,... % Speed Margin (knt), dV '1.2', '1.2' ,... % Sea Margin, (1+x) '10629'}; '10629' }; % Full Load Weight (metric ton), WFL def2 = {'15.24' '15.24', ,'3.6' '3.6'}; }; % Bulb Parameters def3 = {'.839' '.839'}; }; % Waterplane Area Coefficient, CW def4 = {'2900' '2900'}; }; % Wetted Surface Area (m^2), SS def5 = {'5.30e-004' '5.30e-004'}; }; % Model-Ship Correlation Allowance, CA def6 = {'8' '8'}; }; % Immersed area of the transom at zero speed def7 = {'30.00' '30.00'}; }; % The half angle entrance in degrees %---------------------------------------------------------------------% CONSTANTS prompt ={'Seawater ={ 'Seawater Viscosity (m^2/sec)' ,... 'Seawater Density (kg/m^3)' ,... 'g (m/sec^2) '}; ' }; data0 = inputdlg(prompt, ' CONSTANTS ',1,def0); ' ,1,def0); if isempty(data0)==1 clear return end cst1 = str2num(char(data0(1))); cst2 = str2num(char(data0(2))); cst3 = str2num(char(data0(3))); NuSW = cst1; % Seawater Viscosity (m^2/sec) RhoSW = cst2; % Seawater Density (kg/m^3) g = cst3; % (m/sec^2) %---------------------------------------------------------------------% BASIC DESIGN PARAMETERS prompt ={'Length ={ 'Length (m), LWL', LWL' ,... 'Beam (m), B', B' ,... 'Draft (m),T', (m),T' ,... 'Molded draft at the forward perpendicular (m),TF' ,... 'Molded draft at the aft perpendicular (m),TA' ,... 'Prismatic Coefficient, CP' ,... 'Midship Section Coefficient-CM (Maximum Section Coefficient, CX )' ,... 'LCB from midships as percent of Length, (+): forward, (-): abaft' ,...


Page 2 May May 5, 2007 2007

'Service Speed (knt), Vs' Vs', ,... 'Speed Margin (knt), dV' ,... 'Sea Margin, (1+x)', (1+x)' ,... 'Full Load Weight (metric ton), WFL' }; data1 = inputdlg(prompt, ' BASIC DESIGN PARAMETERS ' ,1,def1); if isempty(data1)==1 clear return end p1 = str2num(char(data1(1))); p2 = str2num(char(data1(2))); p3 = str2num(char(data1(3))); p4 = str2num(char(data1(4))); p5 = str2num(char(data1(5))); p6 = str2num(char(data1(6))); p7 = str2num(char(data1(7))); p8 = str2num(char(data1(8))); p9 = str2num(char(data1(9))); p10 = str2num(char(data1(10))); p11 = str2num(char(data1(11))); p12 = str2num(char(data1(12))); LWL = p1 ; % Length (m) B = p2 ; % Beam (m) T = p3 ; % Draft (m) TF = p4 ; % Molded draft at the forward perpendicular % Gemi basindaki su cekimi (m) TA = p5 ; % Molded draft at the aft perpendicular % Gemi kicindaki su cekimi (m) CP = p6 ; % Prismatic Coefficient: CX = p7 ; % Maximum Section Coefficient % (CM) Midship Section Coefficient LCB = p8 ; % LCB from midships % as percent of length, (-) = aft Vs = p9 ; % Service Speed (knt) dV = p10 ; % Speed Margin (knt) onepx = p11 ; % Sea Margin WFL = p12 ; % Full Load Weight (metric ton) %---------------------------------------------------------------------BULB = MENU('The MENU( 'The ship has a bulb ? ' ,... 'YES', 'YES' ,'NO' 'NO', ,'CANCEL' 'CANCEL'); ); if BULB == 1 % Bas dikeyden itibaren yumru basin yanal alani row1 = ['Cross-sectional area of bulb in the vertical' ,... ' plane intersecting the stem contour at the water' , ... ' surface (m^2), ABT ' ]; % Yumru bas, yanal alan (ABT) merkezinin omurgadan itibaren yuksekligi) row2 = ['Height of the centroid of the area ABT above' ,... ' the base line (m) ']; ' ]; prompt ={row1,row2}; data2 = inputdlg(prompt, ' BULB PARAMETERS ',1,def2); ' ,1,def2); if isempty(data2)==1 clear return end p12 = str2num(char(data2(1))); p13 = str2num(char(data2(2))); ABT = p12 ; % cross-sectional area of bulb in


hB

=

C2

=

% the vertical plane intersecting % the stem contour at the water % surface (Bas dikeyden itibaren % yumru basin yanal alani) p13 ; % height of the centroid of the % area ABT above the base line % (yumru bas, yanal alan (ABT) % merkezinin omurgadan itibaren % yuksekligi) exp(-1.89*sqrt(.56*ABT^1.5/(B*T*(.31*ABT^.5+TF-hB)))); exp(-1.89*sqrt(.56*ABT^1.5/ (B*T*(.31*ABT^.5+TF-hB)))); % C2 accounts for the effect of a % bulb. If no bulb is fitted C2 = 1 2 0 ; 0 ; 1; 3

Page 3 May May 5, 2007 2007

elseif BULB == ABT = hB = C2 = elseif BULB == clear return end %---------------------------------------------------------------------DELTAFL = WFL ; % Full Load Displacement % (equal to full load weight) VFL = DELTAFL*1000/RhoSW ; % Volume at LWL (m^3) CB = CP*CX ; % Block Coefficient LR = LWL*(1-CP+.06*CP*LCB/(4*CP LWL*(1-CP+.06*CP*LCB/(4*CP-1)); -1)); % Length of the run (Giris boyu) %---------------------------------------------------------------------CWSEL = MENU('Waterplane MENU( 'Waterplane Area Coefficient' ,... 'USER INPUT', INPUT' ,'ESTIMATED' 'ESTIMATED', ,'CANCEL' 'CANCEL'); ); if CWSEL == 1 prompt = { 'Waterplane Area Coefficient, CW' }; data3 = inputdlg(prompt, ' WATERPLANE AREA COEFFICIENT ' ,1,def3); if isempty(data3)==1 clear return end p14 = str2num(char(data3(1))); CW = p14 ; % Waterplane Area Coefficient elseif CWSEL == 2 CW = .236 + .836*CP; % Reference:"Hydrodynamics in Ship Design" % by Saunders, SNAME 1957 Vol II (pg 466) elseif CWSEL == 3 clear return end %---------------------------------------------------------------------SEL1 = {'- 25 for pram with gondola' , ... '- 10 for V-shaped sections' ,... ' 0 for normal section type' ,... '+10 for U-shaped sections with Hogner stern' }; [SELECTcstern,OK1]=listdlg( 'Name' 'Name', , ... 'Shape of the afterbody' ,... % KIC FORM 'SelectionMode' ,'single' 'single', ,... 'ListSize',[400 'ListSize' ,[400 200],... 200], ... 'ListString',SEL1); 'ListString' ,SEL1); if OK1 == 0


Page 4 May May 5, 2007 2007

clear return end %---------------------------------------------------------------------if SELECTcstern == 1 % - 25 for pram with gondola cstern = -25; namestern = [cellstr( '- 25 for pram with gondola' )]; elseif SELECTcstern == 2 % - 10 for V-shaped sections cstern = -10; namestern = [cellstr('[cellstr( '- 10 for V-shaped sections' )]; elseif SELECTcstern == 3 % 0 for normal section type cstern = 0; namestern = [cellstr('0 [cellstr( '0 for normal section type' )]; elseif SELECTcstern == 4 % +10 for U-shaped sections with Hogner stern cstern = 10; namestern = [cellstr('+10 [cellstr( '+10 for U-shaped sections with Hogner stern' )]; end %---------------------------------------------------------------------c = 1 + .011 * cstern ; % coefficient accounting for the specific % shape of the afterbody [2] k1 = -1 + ... .93+.487118*c*(B/LWL)^1.06806*(T/LWL)^.46106* ... (LWL/LR)^.121563*(LWL^3/VFL)^.36486*(1-CP)^-.604247; % k1: Form factor of bare hull [2] %---------------------------------------------------------------------WETS = MENU('Wetted MENU( 'Wetted Surface Area ' ,... 'USER INPUT', INPUT' ,'ESTIMATED' 'ESTIMATED', ,'CANCEL' 'CANCEL'); ); if WETS == 1 prompt = { 'Wetted Surface Area (m^2), SS' }; data4 = inputdlg(prompt, ' WETTED SURFACE ',1,def4); ' ,1,def4); if isempty(data4)==1 clear return end p15 = str2num(char(data4(1))); SS = p15 ; % Wetted Surface Area % (Tekne Islak Alani) elseif WETS == 2 SS = LWL*(2*T+B)*CX^.5*(.4530+.4 LWL*(2*T+B)*CX^.5*(.4530+.4425*CB425*CB- ... .2862*CX-.003467*B/T+.3696*CW)+2.38*ABT/CB; elseif WETS == 3 clear return end %---------------------------------------------------------------------Vt = dV+Vs; % Trial Speed V = [Vt,linspace(0+eps,Vs+5)] ; %---------------------------------------------------------------------% FRICTIONAL RESISTANCE ACCORDING TO THE ITTC-1957 FRICTION FORMULA % (ITTC 1957 YE GORE SURTUNME DIRENCI) RN = .5144*V*LWL/NuSW ; % Reynolds # CF = .075./(log10(RN)-2).^2; % Frictional Resistance Coefficient RF = .5*(RhoSW*SS*(.5144*V).^2. .5*(RhoSW*SS*(.5144*V).^2.*CF); *CF); %(Newton) %---------------------------------------------------------------------% MODEL-SHIP CORRELATION RESISTANCE (GEMI-MODEL KORELASYON DIRENCI) CASEL = MENU('Model-Ship MENU( 'Model-Ship Correlation Allowance' ,... 'USER INPUT', INPUT' ,'ESTIMATED' 'ESTIMATED', ,'CANCEL' 'CANCEL'); );

powerpredict.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. if

Page 5 May May 5, 2007 2007

CASEL == 1 prompt = { 'Model-Ship Correlation Allowance, CA' }; data5 = inputdlg(prompt, ' MODEL-SHIP CORRELATION RESISTANCE ' ,1,def5); if isempty(data5)==1 clear return end p16 = str2num(char(data5(1))); CA = p16 ; elseif CASEL == 2 if TF/LWL > .04 CA = .006*(LWL+100)^-.16-.00205 ; elseif TF/LWL <= .04 CA = .006*(LWL+100)^-.16-.00205 + .003*(LWL/7.5)^.5*CB^4*C2* .003*(LWL/7.5)^.5*CB^4*C2*(.04-TF/LWL); (.04-TF/LWL); end elseif CASEL == 3 clear return end RA = .5*(RhoSW*SS*(.5144*V).^2. .5*(RhoSW*SS*(.5144*V).^2.*CA); *CA); %(Newton) %---------------------------------------------------------------------% WAVE-MAKING AND WAVE-BREAKING RESISTANCE % (DALGA YAPMA VE DALGA KIRILMA DIRENCI) Fn = .5144*V/(g*LWL)^.5 ; % Froude # ATSEL = MENU('Immersed MENU( 'Immersed area of the transom at zero speed' ,... 'USER INPUT', INPUT' ,'ESTIMATED' 'ESTIMATED', ,'CANCEL' 'CANCEL'); ); if ATSEL == 1 prompt = {['Immersed {[ 'Immersed area of the transom at zero speed, (m^2) AT' ]}; data6 = inputdlg(prompt, ' WAVE-MAKING AND WAVE-BREAKING RESISTANCE ' ,1,def6); if isempty(data6)==1 clear return end p16 = str2num(char(data6(1))); AT = p16 ; % Sifir hizda ayna kicin suya batmis kisminin alani elseif ATSEL == 2 AT = B*T*CX*.3 ; % U/W (under water) transom area % approximation. elseif ATSEL == 3 clear return end C3 = 1-.8*AT/(B*T*CX) ; if B/LWL <= .11 C4 = .229577*(B/LWL)^0.33333; elseif B/LWL > .11 | B/LWL < .25 C4 = B/LWL; elseif B/LWL >= .25 C4 = .5-.0625*(LWL/B); end if CP < 0.8 C5 = 8.07981*CP-13.8673*CP^2+6. 8.07981*CP-13.8673*CP^2+6.984388*CP^3; 984388*CP^3; else C5 = 1.73014-.7067*CP; end if LWL^3/VFL <= 512 C6 = -1.69385;


Page 6 May May 5, 2007 2007

elseif LWL^3/VFL > 512 | LWL^3/VFL < 1727 C6 = -1.69385 + (LWL/VFL^(1/3)-8)/2.36; else LWL^3/VFL >= 1727 C6 = 0; end if LWL/B <12 lambda = 1.446*CP-.03*LWL/B ; else lambda = 1.446*CP-.36 ; end iESEL = MENU('The MENU( 'The half angle entrance in degrees' ,... 'USER INPUT', INPUT' ,'ESTIMATED' 'ESTIMATED', ,'CANCEL' 'CANCEL'); ); if iESEL == 1 prompt = { 'The half angle entrance in degrees, iE' }; data7 = inputdlg(prompt, ' WAVE-MAKING AND WAVE-BREAKING RESISTANCE ' ,1,def7); if isempty(data7)==1 clear return end p17 = str2num(char(data7(1))); iE = p17 ; elseif iESEL == 2 iE = 1+89*exp(-(LWL/B)^.80856*(1 1+89*exp(-(LWL/B)^.80856*(1-CW)^.30484* -CW)^.30484* ... (1-CP-.0225*LCB)^.6367*(LR/B)^.34574*(100*VFL/LWL^3)^0.16302); elseif iESEL == 3 clear return end d = -.9; C1Lo = 2223105*C4^3.78613*(T/B)^1 2223105*C4^3.78613*(T/B)^1.07961*(90-iE)^-1.37565 .07961*(90-iE)^-1.37565 ; % for Fn <0.4 C1Hi = 6919.3*CX^-1.3346*(VFL/LWL 6919.3*CX^-1.3346*(VFL/LWL^3)^2.00977*((LWL/B)-2)^1.40 ^3)^2.00977*((LWL/B)-2)^1.40692 692 ; % for Fn>0.55 m1Lo = .0140407*LWL/T-1.75254*(VF .0140407*LWL/T-1.75254*(VFL^(1/3))/LWL-4.79323*B/LWL-C L^(1/3))/LWL-4.79323*B/LWL-C5; 5; % for Fn <0.4 m1Hi = -7.2035*(B/LWL)^.326869*(T -7.2035*(B/LWL)^.326869*(T/B)^.605375 /B)^.605375 ; % for Fn>0.55 VFn4 = .4 *(g*LWL)^.5/.5144 ; % (knt) for Fn <0.4 VFn55 = .55*(g*LWL)^.5/.5144 ; % (knt) for Fn <0.55 m2 = C6*.4*exp(-.034*Fn.^-3.29); m2VFn4 = C6*.4*exp(-.034 *.4^-3.29); m2VFn55 = C6*.4*exp(-.034*.55^-3.29) C6*.4*exp(-.034*.55^-3.29); ; RWLoV = C1Lo*C2*C3*1000*DELTAFL*g* C1Lo*C2*C3*1000*DELTAFL*g*exp(m1Lo*Fn.^d+m2.*cos(lambd exp(m1Lo*Fn.^d+m2.*cos(lambda*Fn.^-2)); a*Fn.^-2)); % for Fn <0.4 RWHiV = C1Hi*C2*C3*1000*DELTAFL*g* C1Hi*C2*C3*1000*DELTAFL*g*exp(m1Hi*Fn.^d+m2.*cos(lambd exp(m1Hi*Fn.^d+m2.*cos(lambda*Fn.^-2)); a*Fn.^-2)); % for Fn>0.55 RWLoVFn4 = C1Lo*C2*C3*1000*DELTAFL*g*exp(m1Lo* .4^d+m2VFn4 *cos(lambda *.4^-2)); RWHiVFn55= C1Hi*C2*C3*1000*DELTAFL*g*exp(m1Hi*.55^d+m2VFn55*cos(lambda*.55^-2)); RWMidV = RWLoVFn4+(10*Fn-4)*(RWHiVF RWLoVFn4+(10*Fn-4)*(RWHiVFn55-RWLoVFn4)/1.5 n55-RWLoVFn4)/1.5 ; RW = []; % Wave Making Resistance *** for i = 1:length(V) if V(i) <= VFn4 RW(i)= RWLoV(i); elseif V(i) >= VFn55 RW(i)= RWHiV(i); else RW(i)= RWMidV(i);


Page 7 May May 5, 2007 2007

end end %---------------------------------------------------------------------% RESISTANCE OF APPENDAGES (TAKINTILARIN DIRENCI) SEL2 = {'rudder behind skeg (1+k2 = 1.5-2.0)' , ... 'rudder behind stern (1+k2 = 1.3-1.5)' ,... 'twin-screw balance rudders (1+k2 = 2.8)' ,... 'shaft brackets (1+k2 = 3.0)' ,... 'skeg (1+k2 = 1.5-2.0)' ,... 'strut bossing (1+k2 = 3.0)' ,... 'hull bossing (1+k2 = 2.0)' ,... 'shafts (1+k2 = 2.0-4.0)' ,... 'stabilizer fins (1+k2 = 2.8)' ,... 'dome (1+k2 = 2.7)', 2.7)' ,... 'bilge keels (1+k2 = 1.4)' }; [SELECTapp,OK2]=listdlg( 'Name' 'Name', ,... 'Appendages (Use CTRL key for multiselection !' ,... 'SelectionMode' ,'multiple' 'multiple', ,... 'ListSize',[400 'ListSize' ,[400 200],... 200], ... 'ListString',SEL2); 'ListString' ,SEL2); if OK2 == 0 clear return end defapponeplusk2 = {'1.75' { '1.75', ,'1.4' '1.4', ,'2.8' '2.8', ,'3.0' '3.0', ,'1.75' '1.75', ,'3.0' '3.0', ,... '2.0', '2.0' ,'3.0' '3.0', ,'2.8' '2.8', ,'2.7' '2.7', ,'1.4' '1.4'}; }; defarea = {'23.6' { '23.6', ,'23.6' '23.6', ,'10' '10', ,'4' '4', ,'5' '5', ,'6' '6', ,'2' '2', ,'10' '10', ,'40' '40', ,'12' '12', ,'32' '32'}; }; clear prompt for i = 1:length(SELECTapp) def9(i) = defapponeplusk2(SELECTapp(i) defapponeplusk2(SELECTapp(i)); ); promptapp(i)= SEL2(SELECTapp(i)); def8(i) = defarea(SELECTapp(i)); end data8 = inputdlg(promptapp, ' Wetted Area of Appendages (m^2), SAPP' ,1,def8); if isempty(data8)==1 clear return end SAPP = str2num(char(data8)); data9 = inputdlg(promptapp, ' 1 + k2',1,def9); k2' ,1,def9); if isempty(data9)==1 clear return end oneplusk2 = str2num(char(data9)); oneplusk2eq= sum(SAPP.*oneplusk2)/sum(SAPP); prompt

= { 'Tunnel diameter, (m) d' ,... 'Tunnel drag coefficient, CBTO 0.003-0.012' }; def10 = {'0' '0', , '0.004' '0.004'}; }; data10 = inputdlg(prompt, ... ' The Resistance of bow thruster tunnel openings' ,1,def10); if isempty(data10)==1 clear return end dtunnel = str2num(char(data10(1)));

powerpredict.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. CBTO Ctunnel RAPP

Page 8 May May 5, 2007 2007

= str2num(char(data10(2))); = RhoSW*(.5144*V).^2*pi*dtunne RhoSW*(.5144*V).^2*pi*dtunnel^2*CBTO; l^2*CBTO; = 0.5*RhoSW*(.5144*V).^2*sum(S 0.5*RhoSW*(.5144*V).^2*sum(SAPP)*oneplusk2eq.*CF+Ctunnel APP)*oneplusk2eq.*CF+Ctunnel; ; %(Newton) %---------------------------------------------------------------------% ADDITIONAL PRESSURE RESISTANCE OF BULBOUS BOW NEAR THE WATER SURFACE % (YUMRU BASIN ILAVE BASINC DIRENCI) PB = 0.56*(ABT^.5)/(TF-1.5*hB); % measure of immersion Fni = .5144*V./((g*(TF-hB-0.25*( .5144*V./((g*(TF-hB-0.25*(ABT^.5))+0.15*(.5144*V).^2). ABT^.5))+0.15*(.5144*V).^2).^.5); ^.5); % Froude Number based on the immersion RB = 0.11*exp(-3*(PB^(-2)))*(Fn 0.11*exp(-3*(PB^(-2)))*(Fni.^3)*(ABT^1.5)*RhoSW i.^3)*(ABT^1.5)*RhoSW *g./(1+Fni.^2); %(Newton) %---------------------------------------------------------------------% ADDITIONAL PRESSURE RESISTANCE OF IMMERSED TRANSOM STERN % (AYNA KICIN ILAVE BASINC DIRENCI) FnT = .5144*V./((2*g*AT/(B+B*CW) .5144*V./((2*g*AT/(B+B*CW))^.5); )^.5); for i = 1:length(FnT) if FnT(i) < 5 C7(i) = 0.2*(1-0.2*FnT(i)); else C7(i) = 0; end end RTR = 0.5*RhoSW*(0.5144*V).^2*AT 0.5*RhoSW*(0.5144*V).^2*AT.*C7; .*C7; %(Newton) %---------------------------------------------------------------------% TOTAL RESISTANCE (TOPLAM DIRENC) %(Newton) R1 = (RF*(1+k1)+ RA) ; R2 = (RAPP + RB + RTR) ; RT = R1 + R2 + RW ; %---------------------------------------------------------------------% OUTPUT Vplot = V(2:length(V)); R1plot = R1(2:length(V)); R2plot = R2(2:length(V)); RWplot = RW(2:length(V)); RTplot = RT(2:length(V)); % ***** R1Vt = R1(1); R2Vt = R2(1); RWVt = RW(1); RTVt = RT(1); % ***** PEtplot = .5144*Vplot.*RTplot; PEt = .5144*Vt.*RTVt; PEsplot = onepx * PEtplot ; PEs = onepx * PEt ; Pinterpdata = [PEsplot' Vplot']; OUTPUTDATA1 = [ V(2:length(V))' , RF(2:length(V))' , ... k1*RF(2:length(V))' , RA(2:length(V))' , ... R1(2:length(V))' , RAPP(2:length(V))' , ... RB(2:length(V))' , RTR(2:length(V))' , ... R2(2:length(V))', RW(2:length(V))' ]; OUTPUTDATA2 = [ V(2:length(V))' , RT(2:length(V))' , ... (PEtplot/1000)', (PEsplot/1000)' ]; OUTPUTTEXT1 = {'V { 'V (knots)', (knots)' ,'RF (N)', (N)','k1*RF (N)', (N)' ,'RA (N)', (N)' ,'(RF*(1+k1)+ RA) (N)', (N)' ,... 'RAPP (N)', (N)' ,'RB (N)', (N)','RTR (N)', (N)','(RAPP + RB + RTR) (N)' ,'RW (N)'}; (N)'}; OUTPUTTEXT2 = {'V { 'V (knots)', (knots)' ,'RT (N)', (N)','PE (kW) trial', trial' ,'PE (kW) service'}; service' }; %---------------------------------------------------------------------% FIGURE OUTPUT


Page 9 May May 5, 2007 2007

%---------------------------------------------------------------------h1 = figure(1); subplot(1,4,1),plot(Vplot,R1plot/1000, Vt,R1Vt/1000, 'or' 'or') ) axis([0 max(V) 0 max(RT/1000)]) title('\fontsize{10}\bf{(R_F*(1+k1)+ title('\fontsize{10}\bf{(R_F*(1+k1)+ R_A)}' ); xlabel('\fontsize{10}\bf{V xlabel( '\fontsize{10}\bf{V (knot)}' ); ylabel('\fontsize{10}\bf{R ylabel( '\fontsize{10}\bf{R (kN)}' ); text(Vt,R1Vt/1000,[ '\fontsize{10}\bf{ }',num2str(R1Vt/1000), }' ,num2str(R1Vt/1000), ' kN']); kN']); grid subplot(1,4,2),plot(Vplot,R2plot/1000, Vt,R2Vt/1000, 'or' 'or') ) axis([0 max(V) 0 max(RT/1000)]) title('\fontsize{10}\bf{+(R_{APP}+R_B+R_{TR})}' title('\fontsize{10}\bf{+(R_{APP}+R_B+R_{TR})}' ); xlabel('\fontsize{10}\bf{V xlabel( '\fontsize{10}\bf{V (knot)}' ); text(Vt,R2Vt/1000,[ '\fontsize{10}\bf{ }',num2str(R2Vt/1000), }' ,num2str(R2Vt/1000), ' kN']); kN']); grid subplot(1,4,3),plot(Vplot,RWplot/1000, Vt,RWVt/1000, 'or' 'or') ) axis([0 max(V) 0 max(RT/1000)]) title('\fontsize{10}\bf{+R_W}' title('\fontsize{10}\bf{+R_W}' ); xlabel('\fontsize{10}\bf{V xlabel( '\fontsize{10}\bf{V (knot)}' ); text(Vt,RWVt/1000,[ '\fontsize{10}\bf{ }',num2str(RWVt/1000), }' ,num2str(RWVt/1000), ' kN']); kN']); grid subplot(1,4,4), plot(Vplot,RTplot/1000, Vt,RTVt/1000, 'or' 'or') ) axis([0 max(V) 0 max(RT/1000)]) title('\fontsize{10}\bf{=R_T}' title('\fontsize{10}\bf{=R_T}' ); xlabel('\fontsize{10}\bf{V xlabel( '\fontsize{10}\bf{V (knot)}' ); text(Vt,RTVt/1000,[ '\fontsize{10}\bf{ }',num2str(RTVt/1000), }' ,num2str(RTVt/1000), ' kN']); kN']); grid set(h1,'Position' set(h1, 'Position',[1,1,1400,950]) ,[1,1,1400,950]) h2 = figure(2); plot(Vplot,PEtplot/1000, ... Vplot,PEsplot/1000, '--' '--', ,... 'linewidth',[2]) 'linewidth' ,[2]) hold on line([Vt Vt], [0 PEt/1000], 'color' 'color',[0 ,[0 0 1],'linewidth' 1],'linewidth',[2]) ,[2]) line([Vt 0], [PEt/1000 PEt/1000], 'color' 'color',[0 ,[0 0 1],'linewidth' 1],'linewidth',[2], ,[2],'marker' 'marker', ,'<' '<') ) text(0.2,PEt/1000+800,[ '\fontsize{10}\bf{ }',num2str(PEt/1000), }' ,num2str(PEt/1000), ' kW (Trial)']); (Trial)' ]); text(Vt+0.2,800,[ '\fontsize{10}\bf{ }',num2str(Vt), }' ,num2str(Vt), ' knots (Trial)']); (Trial)' ]); title('\fontsize{15}\bf{POWER}' title('\fontsize{15}\bf{POWER}' ); xlabel('\fontsize{12}\bf{V xlabel( '\fontsize{12}\bf{V (knot)}' ); ylabel('\fontsize{12}\bf{P_E ylabel( '\fontsize{12}\bf{P_E (kW)}' ); legend('\fontsize{15}\bf{Trial legend( '\fontsize{15}\bf{Trial Power}' , 'Service Power',2) Power' ,2) grid set(h2,'Position' set(h2, 'Position',[1,1,1400,950]) ,[1,1,1400,950]) %---------------------------------------------------------------------% EXCEL OUTPUT %---------------------------------------------------------------------question = questdlg('EXCEL questdlg( 'EXCEL OUTPUT ?', ?' ,... ' KIVANC', KIVANC' ,'YES' 'YES', ,'NO' 'NO', ,'NO' 'NO'); ); switch question case 'YES' exceloutput = 1; case 'NO' exceloutput = 2; end if exceloutput == 1 Excel = actxserver( 'Excel.Application' ); set(Excel, 'Visible' , 1);


Page 10 May May 5, 2007 2007

Workbook = invoke(Excel.Workbooks, 'Add' 'Add'); ); %-----------------------------------------------------values1 = [ ... NuSW % Seawater Viscosity (m^2/sec) RhoSW % Seawater Density (kg/m^3) g]; % (m/sec^2) %-----------------------------------------------------values2 = [ ... LWL % Length (m) B % Beam (m) T % Draft (m) TF % Molded draft at the forward perpendicular TA % Molded draft at the aft perpendicular CP % Prismatic Coefficient: CX % Maximum Section Coefficient,(CM) CW % Waterplane Area Coefficient CB % Block Coefficient (CP*CX) LCB % LCB from midships as percent of length Vs % Service Speed (knt) dV % Speed Margin (knt) onepx % Sea Margin WFL % Full Load Weight (metric ton) ABT % Bulb Parameter hB]; % Bulb Parameters %-----------------------------------------------------values3 = namestern; % Shape of the afterbody %-----------------------------------------------------values4 = [ ... SS % Wetted Surface Area (m^2), SS CA % Model-Ship Correlation Allowance, CA AT % Immersed area of the transom at zero speed iE % The half angle entrance in degrees dtunnel % Tunnel diameter, (m) d CBTO]; % Tunnel drag coefficient, CBTO 0.003-0.012 %-----------------------------------------------------values5 = [def8' def9']; % WETTED AREA and 1+k2 VALUES OF APPENDAGES (m^2) %-----------------------------------------------------text1 = [... cellstr('CONSTANTS' cellstr( 'CONSTANTS'); );... ... cellstr('Seawater cellstr( 'Seawater Viscosity (m^2/sec)' );... );... cellstr('Seawater cellstr( 'Seawater Density (kg/m^3)' );... );... cellstr('g cellstr( 'g (m/sec^2) '); ' );... ... cellstr('DESIGN cellstr( 'DESIGN PARAMETERS'); PARAMETERS' );... ... cellstr('Length cellstr( 'Length (m), LWL'); LWL' );... ... cellstr('Beam cellstr( 'Beam (m), B'); B' );... ... cellstr('Draft cellstr( 'Draft (m),T'); (m),T' );... ... cellstr('Molded cellstr( 'Molded draft at the forward perpendicular (m),TF' );... );... cellstr('Molded cellstr( 'Molded draft at the aft perpendicular (m),TA' );... );... cellstr('Prismatic cellstr( 'Prismatic Coefficient, CP' );... );... cellstr('Midship cellstr( 'Midship Section Coefficient-CM (Maximum Section Coefficient, CX)' );... );... cellstr('Waterplane cellstr( 'Waterplane Area Coefficient, CW' );... );... cellstr('Block cellstr( 'Block Coefficient (CP*CX)' );... );... cellstr('LCB cellstr( 'LCB from midships as percent of Length, (+): forward, (-): abaft' );... );... cellstr('Service cellstr( 'Service Speed (knt), Vs' );... );... cellstr('Speed cellstr( 'Speed Margin (knt), dV' );... );... cellstr('Sea cellstr( 'Sea Margin, (1+x)'); (1+x)' );... ... cellstr('Full cellstr( 'Full Load Weight (metric ton), WFL' );... );...


Page 11 May May 5, 2007 2007

cellstr(['Cross-sectional cellstr([ 'Cross-sectional area of bulb in the vertical plane intersecting' ,... ' the stem contour at the water surface (m^2)' ]);... ]);... cellstr('Height cellstr( 'Height of the centroid of the area ABT above the base line (m)' );... );... cellstr('Shape cellstr( 'Shape of the afterbody' );... );... cellstr('Wetted cellstr( 'Wetted Surface Area (m^2), SS' );... );... cellstr('Model-Ship cellstr( 'Model-Ship Correlation Allowance' );... );... cellstr('Immersed cellstr( 'Immersed area of the transom at zero speed, (m^2) AT' );... );... cellstr('The cellstr( 'The half angle entrance in degrees' );... );... cellstr('Tunnel cellstr( 'Tunnel diameter, (m) d' );... );... cellstr('Tunnel cellstr( 'Tunnel drag coefficient, CBTO 0.003-0.012' );... );... cellstr('WETTED cellstr( 'WETTED AREA and 1+k2 VALUES OF APPENDAGES (m^2)' )]; %-----------------------------------------------------text2 = promptapp'; %-----------------------------------------------------space = 2; Ltext1 = length(text1); firstcell = [ 'A' num2str(space)]; lastcell = [ 'A' num2str(space+Ltext1-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , text1, 'ColumnWidth' 'ColumnWidth',87); ,87); %-----------------------------------------------------space = space+Ltext1; Ltext2 = length(text2); firstcell = [ 'A' num2str(space)]; lastcell = [ 'A' num2str(space+Ltext2-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , text2, 'ColumnWidth' 'ColumnWidth',87); ,87); %-----------------------------------------------------space = 3; Lvalues1 = length(values1); firstcell = [ 'B' num2str(space)]; lastcell = [ 'B' num2str(space+Lvalues1-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , values1, 'ColumnWidth' 'ColumnWidth',21); ,21); %-----------------------------------------------------space = space+Lvalues1+1; Lvalues2 = length(values2); firstcell = [ 'B' num2str(space)]; lastcell = [ 'B' num2str(space+Lvalues2-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , values2, 'ColumnWidth' 'ColumnWidth',21); ,21); %-----------------------------------------------------space = space+Lvalues2; Lvalues3 = length(values3); firstcell = [ 'B' num2str(space)]; lastcell = [ 'B' num2str(space+Lvalues3-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , values3, 'ColumnWidth' 'ColumnWidth',21); ,21); %-----------------------------------------------------space = space+Lvalues3; Lvalues4 = length(values4); firstcell = [ 'B' num2str(space)];

powerpredict.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. lastcell = [ 'B' num2str(space+Lvalues4-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , values4, 'ColumnWidth' 'ColumnWidth',21); ,21); %-----------------------------------------------------space = space+Lvalues4+1; Lvalues5 = length(values5); firstcell = [ 'B' num2str(space)]; lastcell = [ 'C' num2str(space+Lvalues5-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , values5, 'ColumnWidth' 'ColumnWidth',21); ,21); %-----------------------------------------------------% Make the second sheet active. ActiveWorkbook = get(Excel, 'ActiveWorkBook' ); Sheets = get(ActiveWorkbook, 'Sheets' 'Sheets'); ); Sheet2 = Item(Sheets, 2); Activate(Sheet2); % Get a handle to the active sheet. ActiveSheet = get(Excel, 'ActiveSheet' 'ActiveSheet'); ); %-----------------------------------------------------space = 2; firstcell = [ 'A' num2str(space)]; lastcell = [ 'J' num2str(space)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , OUTPUTTEXT1, 'ColumnWidth' 'ColumnWidth',15); ,15); %-----------------------------------------------------space = space+1; [RowDat ColDat] = size(OUTPUTDATA1); firstcell = [ 'A' num2str(space)]; lastcell = [ 'J' num2str(space+RowDat-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , OUTPUTDATA1, 'ColumnWidth' 'ColumnWidth',15); ,15); %-----------------------------------------------------% Make the third sheet active. ActiveWorkbook = get(Excel, 'ActiveWorkBook' ); Sheets = get(ActiveWorkbook, 'Sheets' 'Sheets'); ); Sheet3 = Item(Sheets, 3); Activate(Sheet3); % Get a handle to the active sheet. ActiveSheet = get(Excel, 'ActiveSheet' 'ActiveSheet'); ); %-----------------------------------------------------space = 2; firstcell = [ 'A' num2str(space)]; lastcell = [ 'D' num2str(space)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , OUTPUTTEXT2, 'ColumnWidth' 'ColumnWidth',20); ,20); %-----------------------------------------------------space = space+1; [RowDat ColDat] = size(OUTPUTDATA2); firstcell = [ 'A' num2str(space)]; lastcell = [ 'D' num2str(space+RowDat-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet, 'Range' 'Range',firstcell,lastcell); ,firstcell,lastcell);

Page 12 May May 5, 2007 2007

powerpredict.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. set(ActivesheetRange, 'Value' 'Value', , OUTPUTDATA2, 'ColumnWidth' 'ColumnWidth',20); ,20); %-----------------------------------------------------prompt ={'ENTER ={ 'ENTER THE FILE NAME FOR THE RESULTS' }; defname = {['Project' {[ 'Project',num2str(1)]}; ,num2str(1)]}; fname = inputdlg(prompt, ' FILE NAME',1,defname); NAME' ,1,defname); filename = char(fname); if isempty(filename)==1 clc; return end w=pwd; invoke(Workbook, 'SaveAs' 'SaveAs', , [w,'\' [w,'\',filename, ,filename,'.xls' '.xls']); ]); delete(Excel); end %---------------------------------------------------------------------% REFERENCES: %---------------------------------------------------------------------% [1] HOLTROP,J., MENNEN,.G.J., AN APPROXIMATE % POWER PREDICTION METHOD % [2] HOLTROP,J., A STATISTICAL RE-ANALYSIS OF % RESISTANCE AND PROPULSION DATA % [3] PNA VOLUME II, RESISTANCE, PROPULSION AND % VIBRATION Chapter 5, Section 8.12, Pg 90-93 % [4] MIT FUNCTIONAL SHIP DESIGN SYNTHESIS MODEL % [5] DIKILI,C., ON DIZAYN ASAMASINDA GEMI % GUCUNUN BELIRLENMESI VE DEGISIK TIP GEMILER % ICIN MODEL DENEYLERI ILE KARSILASTIRILMASI % [6] PARSONS,M., Power Prediction Program (PPP), % University of Michigan, Department of Naval % Architecture and Marine Engineering, % Conceptual Ship Design Software Environment %----------------------------------------------------------------------

Page 13 May May 5, 2007 2007

propeller.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... ..

Page 1 May May 7, 7, 2007 2007

global Z A TYPE COEFF R ETA0 Kt Kq % Kivanc Ali ANIL, 2007 % 508052003 % % This program uses the same algorithm of propopt.m from % http://web.nps.navy.mil/~me/tsse/TS4001/docs_software.htm % Run powerpredict.m before running this program. % % Propeller B-Series Optimization % TYPE = 1 : KT/J^2 case % 2 : KQ/J^3 case % 3 : KT/J^4 case % 4 : KQ/J^5 case % % Requires that function fcn.m is located in the same directory % as this main program. % It uses a Lagrange multiplier method and utilizes the standard % Matlab function fzero for the optimization problem. % If the program does not converge, try a different initial guess. % %---------------------------------------------------------------------clear KTplot KQplot ETAOplot Jplot h3 = figure(3); h4 = figure(4); h5 = figure(5); h6 = figure(6); close(h3); close(h4); close(h5); close(h6); %---------------------------------------------------------------------% default values defp0 = {'1' '1', ,... % Number of propellers '5'}; '5' }; % Number of blades defp1 = {'0.3695' '0.3695'}; }; % Wake Fraction defp2 = {'0.1869' '0.1869'}; }; % Thrust Deduction Fraction defp3 = {'0.55' '0.55'}; }; % AE/A0 ratio defp4 = {'0.2' '0.2'}; }; % k constant for Keller's Criterion defp5 = {'0.8' '0.8'}; }; % initial guess for J % default values for the first step of iterations % (propeller optimization) defEth = {'0.002' '0.002'}; }; % Threshold value for the difference between % calculated and assumed Propulsive Efficiency defETAD = {'0.7' '0.7'}; }; % Propulsive Efficiency promptETAD = {'Propulsive Efficiency (Assumed)'}; (Assumed)'}; defETAS = {'1' '1', ,'0.85' '0.85', ,'0.98' '0.98'}; }; % Efficiencies defN = {'100' '100'}; }; % Rotation Rate (rpm), N % default values for the first step of iterations % (Prediction of Performance in Service) defETADs = {'0.7' '0.7'}; }; % Propulsive Efficiency promptETADs = {'Propulsive {'Propulsive Efficiency (Assume)'}; (Assume)'}; %---------------------------------------------------------------------prompt ={'Enter ={'Enter Number of propellers ', ',... 'Enter Number of blades '}; '}; datap0 = inputdlg(prompt,' inputdlg(prompt,' PROPELLER PARAMETERS ',1,defp0); ',1,defp0); if isempty(datap0)==1 %clear return end NoP = str2num(char(datap0(1))); Z = str2num(char(datap0(2))); %---------------------------------------------------------------------wSEL = MENU('Wake MENU('Wake Fraction', Fraction',...

propeller.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. 'USER INPUT', INPUT','ESTIMATED' 'ESTIMATED', ,'CANCEL' 'CANCEL'); ); if wSEL == 1 prompt = {'Wake Fraction, w'}; w'}; datap1 = inputdlg(prompt,' inputdlg(prompt,' Wake Fraction ',1,defp1); ',1,defp1); if isempty(datap1)==1 %clear return end w = str2num(char(datap1(1))); elseif wSEL == 2 w = 0.727-1.86*CB+1.749*CB^2; elseif wSEL == 3 %clear return end %---------------------------------------------------------------------tSEL = MENU('Thrust MENU('Thrust Deduction Fraction ', ',... 'USER INPUT', INPUT','ESTIMATED' 'ESTIMATED', ,'CANCEL' 'CANCEL'); ); if tSEL == 1 prompt = {'Thrust Deduction Fraction , t' t'}; }; datap2 = inputdlg(prompt,' inputdlg(prompt,' Thrust Deduction Fraction ',1,defp2); if isempty(datap2)==1 %clear return end t = str2num(char(datap2(1))); elseif tSEL == 2 t = 0.6*w; elseif tSEL == 3 %clear return end %---------------------------------------------------------------------Thrust = (RTVt*(1-t))/NoP; %(Newton) VAdvance = .5144*Vt*(1-w); %m/s etaHull = (1-t)/(1-w); %---------------------------------------------------------------------defD = {num2str(T*.65)}; prompt = {'Propeller Diameter (assumed to be %65 of the draft) '}; '}; dataD = inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller Diameter (m)',1,defD); (m)',1,defD); if isempty(dataD)==1 %clear return end D = str2num(char(dataD(1))); %---------------------------------------------------------------------% Open Water Diameter if NoP == 1 Dfactor = 0.05; else Dfactor = 0.03; end Do = D/(1-Dfactor); %---------------------------------------------------------------------ASEL = MENU('AE/A0 MENU('AE/A0 ratio (BAR)', (BAR)',... 'USER INPUT', INPUT','Keller"s Criterion', Criterion','CANCEL' 'CANCEL'); ); if ASEL == 1

Page 2 May May 7, 7, 2007 2007


Page 3 May May 7, 7, 2007 2007

prompt = {'AE/A0 ratio'}; ratio'}; datap3 = inputdlg(prompt,' inputdlg(prompt,' AE/A0 ratio',1,defp3); ratio',1,defp3); if isempty(datap3)==1 %clear return end A = str2num(char(datap3(1))); elseif ASEL == 2 % Keller"s Criterion h = (D/2)+0.2; % height of shaft centre-line above base Patm = 101300; % (N/m^2) Atmospheric pressure Pv = 1646; % (N/m^2) Vapour pressure of water at 15 °C Po = Patm+RhoSW*g*(T-h); % the static pressure at the shaft centerline prompt = {['Constant {['Constant varying from 0 for transom-stern naval vessels', vessels',... ' to 0.20 for high-powered single screw vessels ']}; ']}; datap4 = inputdlg(prompt,' inputdlg(prompt,' k',1,defp4); k',1,defp4); if isempty(datap4)==1 %clear return end k = str2num(char(datap4(1))); defA = ((1.3+0.3*Z)*Thrust/((Po-Pv ((1.3+0.3*Z)*Thrust/((Po-Pv)*D^2))+k; )*D^2))+k; %defA = round(defA*20)/20; %round to the nearest 0.05 defA = {num2str(defA)}; prompt = {'AE/A0 ratio (BAR) from Keller"s Criterion'}; Criterion'}; %dataA = inputdlg(prompt,' AE/A0 ratio (rounded to the nearest 0.05)',1,defA); dataA = inputdlg(prompt,' inputdlg(prompt,' AE/A0 ratio',1,defA); ratio',1,defA); if isempty(dataA)==1 %clear return end A = str2num(char(dataA(1))); elseif ASEL == 3 %clear return end %---------------------------------------------------------------------TYPE = MENU('Propeller MENU('Propeller B-Series Optimization', Optimization',... ['TYPE=1 : KT/J^2 case - Known thrust (T), diameter (D)', (D)',... ' and advanced velocity (VA)- Nopt (required)'], (required)'],... ... ['TYPE=2 : KQ/J^3 case - Known power (PD), diameter (D)', (D)',... ' and advanced velocity (VA)- Nopt (required)'], (required)'],... ... ['TYPE=3 : KT/J^4 case - Known thrust (T), rotation rate (N)', (N)',... ' and advanced velocity (VA)- Dopt (required)'], (required)'],... ... ['TYPE=4 : KQ/J^5 case - Known power (PD), rotation rate (N)', (N)',... ' and advanced velocity (VA)- Dopt (required)'], (required)'],... ... 'CANCEL'); 'CANCEL' ); %---------------------------------------------------------------------% PROPELLER OPTIMIZATION %---------------------------------------------------------------------check = 1; pauseN = 0; while check == 1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% pauseN = pauseN+1; %---------------------------------------------------------------------if TYPE == 1 % KT/J^2 CTN = Thrust/(RhoSW*VAdvance^2*Do^ Thrust/(RhoSW*VAdvance^2*Do^2); 2); defCOEFF = {num2str(CTN)};


Page 4 May May 7, 7, 2007 2007

prompt = {'KT/J^2 CONFIRMATION'}; CONFIRMATION'}; dataCOEFF= inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defCOEFF); Optimization',1,defCOEFF); if isempty(dataCOEFF)==1 %clear return end COEFF = str2num(char(dataCOEFF(1))) str2num(char(dataCOEFF(1))); ; elseif TYPE == 2 % KQ/J^3 % Threshold value for the difference between % calculated and assumed Propulsive Efficiency prompt = {['Threshold {['Threshold value for the difference between ', ',... 'calculated and assumed Propulsive Efficiency ']}; ']}; datapEth = inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defEth); Optimization',1,defEth); if isempty(datapEth)==1 %clear return end Ethreshold = str2num(char(datapEth(1))); prompt = {'Relative-rotative Efficiency', Efficiency',... 'Maximum continous power at = ', ', 'shaft transmission efficiency'}; efficiency'}; dataETAS= inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defETAS); Optimization',1,defETAS); if isempty(dataETAS)==1 %clear return end etaR = str2num(char(dataETAS(1))); MCPat = str2num(char(dataETAS(2))); etaS = str2num(char(dataETAS(3))); dataETAD = inputdlg(promptETAD,'Propeller inputdlg(promptETAD,'Propeller B-Series Optimization',1,defETAD); Optimization',1,defETAD); if isempty(dataETAD)==1 %clear return end etaD = str2num(char(dataETAD(1))); PD = PEt/etaD; % Delivered power CQN = PD/(2*pi*RhoSW*Do^2*VAdvanc PD/(2*pi*RhoSW*Do^2*VAdvance^3); e^3); defCOEFF = {num2str(CQN)}; prompt = {'KQ/J^3 CONFIRMATION'}; CONFIRMATION'}; dataCOEFF= inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defCOEFF); Optimization',1,defCOEFF); if isempty(dataCOEFF)==1 %clear return end COEFF = str2num(char(dataCOEFF(1))) str2num(char(dataCOEFF(1))); ; elseif TYPE == 3 % KT/J^4 prompt = {'Rotation Rate (rpm), N'}; N'}; dataN= inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defN); Optimization',1,defN); if isempty(dataN)==1 %clear return end N = str2num(char(dataN(1))); CTN = (Thrust*(N/60)^2)/(RhoSW*VA (Thrust*(N/60)^2)/(RhoSW*VAdvance^4); dvance^4); defCOEFF = {num2str(CTN)}; prompt = {'KT/J^4 CONFIRMATION'}; CONFIRMATION'}; dataCOEFF= inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defCOEFF); Optimization',1,defCOEFF); if isempty(dataCOEFF)==1


Page 5 May May 7, 7, 2007 2007

%clear return end COEFF = str2num(char(dataCOEFF(1))) str2num(char(dataCOEFF(1))); ; elseif TYPE == 4 % KQ/J^5 % Threshold value for the difference between % calculated and assumed Propulsive Efficiency prompt = {['Threshold {['Threshold value for the difference between ', ',... 'calculated and assumed Propulsive Efficiency ']}; ']}; datapEth = inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defEth); Optimization',1,defEth); if isempty(datapEth)==1 %clear return end Ethreshold = str2num(char(datapEth(1))); prompt = {'Rotation Rate (rpm), N'}; N'}; dataN = inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defN); Optimization',1,defN); if isempty(dataN)==1 %clear return end N = str2num(char(dataN(1))); prompt = {'Relative-rotative Efficiency', Efficiency',... 'Maximum continous power at = ', ', 'shaft transmission efficiency'}; efficiency'}; dataETAS= inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defETAS); Optimization',1,defETAS); if isempty(dataETAS)==1 %clear return end etaR = str2num(char(dataETAS(1))); MCPat = str2num(char(dataETAS(2))); etaS = str2num(char(dataETAS(3))); dataETAD = inputdlg(promptETAD,'Propeller inputdlg(promptETAD,'Propeller B-Series Optimization',1,defETAD); Optimization',1,defETAD); if isempty(dataETAD)==1 %clear return end etaD = str2num(char(dataETAD(1))); PD = PEt/etaD;% PEt/etaD;% Delivered power CQN = (PD*(N/60)^2)/(2*pi*RhoSW*V (PD*(N/60)^2)/(2*pi*RhoSW*VAdvance^5); Advance^5); defCOEFF = {num2str(CQN)}; prompt = {'KQ/J^5 CONFIRMATION'}; CONFIRMATION'}; dataCOEFF= inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defCOEFF); Optimization',1,defCOEFF); if isempty(dataCOEFF)==1 %clear return end COEFF = str2num(char(dataCOEFF(1))) str2num(char(dataCOEFF(1))); ; elseif TYPE == 5 % CANCEL %clear return end % COEFFICIENTS FOR KT POLYNOMIAL % C1 = .0088049600; C2 =-.2045540000; C3 = .1663510000; C4 = .1581140000; C5 =-.1475810000; C6 =-.4814970000; C7 = .4154370000; C8 = .0144043000; C9 =-.0530054000;


Page 6 May May 7, 7, 2007 2007

C10= .0143481000; C11= .0606826000; C12=-.0125894000; C13= .0109689000; C14=-.1336980000; C15= .0063840700; C16=-.0013271800; C17= .1684960000; C18=-.0507214000; C19= .0854559000; C20=-.0504475000; C21= .0104650000; C22=-.0064827200; C23=-.0084172800; C24= .0168424000; C25=-.0010229600; C26=-.0317791000; C27= .0186040000; C28=-.0041079800; C29=-.0006068480; C30=-.0049819000; C31= .0025983000; C32=-.0005605280; C33=-.0016365200; C34=-.0003287870; C35= .0001165020; C36= .0006909040; C37= .0042174900; C38= .0000565229; C39=-.0014656400; % % COEFFICIENTS FOR KQ POLYNOMIAL % D1 = .0037936800; D2 = .0088652300; D3 =-.0322410000; D4 = .0034477800; D5 =-.0408811000; D6 =-.1080090000; D7 =-.0885381000; D8 = .1885610000; D9 =-.0037087100; D10= .0051369600; D11= .0209449000; D12= .0047431900; D13=-.0072340800; D14= .0043838800; D15=-.0269403000; D16= .0558082000; D17= .0161886000; D18= .0031808600; D19= .0158960000; D20= .0471729000; D21= .0196283000; D22=-.0502782000; D23=-.0300550000; D24= .0417122000; D25=-.0397722000; D26=-.0035002400; D27=-.0106854000; D28= .0011090300; D29=-.0003139120; D30= .0035985000; D31=-.0014212100; D32=-.0038363700; D33= .0126803000; D34=-.0031827800; D35= .0033426800; D36=-.0018349100; D37= .0001124510; D38=-.0000297228; D39= .0002695510; D40= .0008326500; D41= .0015533400; D42= .0003026830; D43=-.0001843000; D44=-.0004253990; D45= .0000869243; D46=-.0004659000; D47= .0000554194; % i = 0; h3 = figure(3); set(h3,'Position' set(h3,'Position',[1,1,1400,950]) ,[1,1,1400,950]) title(['\fontsize{12}\bf{WAGENINGEN-B title(['\fontsize{12}\bf{WAGENINGEN-B SERIES PROPELLERS, FOR }', }', num2str(Z),... num2str(Z),... ' BLADES, A_E/A_O = ',num2str(A), ',num2str(A),', ', P/D = 0.50 TO 1.40']); 1.40']); xlabel('\fontsize{12}\bf{J}' xlabel('\fontsize{12}\bf{J}'); ); ylabel('\fontsize{12}\bf{K_T, ylabel('\fontsize{12}\bf{K_T, 5 x K_Q, \eta_0}'); \eta_0}'); grid axis([0 1.6 0 1.0]) set(gca,'XTick' set(gca,'XTick',[0.00;0.20;0.40;0.60;0.80;1.00;1.20;1.40;1.60]) ,[0.00;0.20;0.40;0.60;0.80;1.00;1.20;1.40;1.60]) set(gca,'XTickLabel' set(gca,'XTickLabel',[0.00;0.20;0.40;0.60;0.80;1.00;1.20;1.40;1.60]) ,[0.00;0.20;0.40;0.60;0.80;1.00;1.20;1.40;1.60]) %set(gca,'YTick',[0.00;0.50;1.25;2.50;4.00;6.00]) %set(gca,'YTickLabel',[0.00;0.50;1.25;2.50;4.00;6.00]) hold on if pauseN == 1 pause end for Ri = 0.5:0.1:1.4 i = i+1; j = 0; for Jj=0:0.001:1.6; j = j+1; % KT1=C1 + C2*Jj + C3*Ri + C4*Ri^2 + C5*Jj^2*A + C6*Jj*Ri*A; KT2=C7*Ri^2*A + C8*Z + C9*Jj^2*Z + C10*Ri*Z + C11*Jj*Ri*Z; KT3=C12*A*Z+C13*Jj*A*Z+C14*Ri^3+C15*Ri^6 + C16*Jj^2*Ri^6; KT4=C17*Jj^3*A+C18*A^2+C19*Jj^2*A^2 + C20*Jj^3*A^2;


Page 7 May May 7, 7, 2007 2007

KT5=C21*Jj*Ri^6*A^2 + C22*Jj^2*Ri^6*A^2 + C23*Ri^3*Z; KT6=C24*Jj*Ri^3*Z + C25*Jj^3*Ri^3*Z + C26*Ri^3*A*Z; KT7=C27*Jj*A^2*Z + C28*Ri^2*A^2*Z + C29*Z^2 + C30*Jj*Z^2; KT8=C31*Jj^2*Z^2 + C32*Jj^3*Z^2 + C33*Jj*Ri^2*Z^2; KT9=C34*Jj*Ri^6*Z^2 + C35*Jj^2*Ri^6*Z^2 + C36*A*Z^2; KT10=C37*Ri^3*A*Z^2+C38*Jj^3*Ri^6*A*Z^2+C39*Ri^3*A^2*Z^2; % % KT POLYNOMIAL % KT = KT1+KT2+KT3+KT4+KT5+KT6+KT7+KT8+KT9+KT10; % KQ1 =D1 + D2*Jj^2 + D3*Jj*Ri + D4*Ri^2 + D5*Ri*A + D6*Jj*Ri*A; KQ2 =D7*(Jj^2*Ri*A) + D8*(Ri^2*A) + D9*Jj*Z+D10*Ri*Z + D11*Jj*Ri*Z; KQ3 =D12*Jj^2*Ri*Z + D13*Jj^2*A*Z + D14*Jj*Ri*A*Z + D15*Ri^2*A*Z; KQ4 =D16*Jj^3*A + D17*Ri^3*A + D18*Jj*Ri^3*A + D19*A^2; KQ5 =D20*Jj*A^2 + D21*Jj^3*A^2 + D22*Ri*A^2 + D23*Jj^3*Ri*A^2; KQ6 =D24*Jj^2*Ri^2*A^2 + D25*Ri^3*A^2 + D26*Ri^6*A^2; KQ7 =D27*Jj^3*Z + D28*Jj^3*Ri^3*Z + D29*Ri^6*Z + D30*Jj^3*A*Z; KQ8 =D31*Ri^6*A*Z +D32*Jj*A^2*Z + D33*Ri^2*A^2*Z; KQ9 =D34*Jj^2*Ri^3*A^2*Z + D35*Ri^6*A^2*Z + D36*Jj*Ri*Z^2; KQ10=D37*Jj^3*Ri^2*Z^2 + D38*Jj^3*Ri^6*Z^2 + D39*Jj*A*Z^2; KQ11=D40*Jj^2*A*Z^2 + D41*Ri^2*A*Z^2 + D42*Ri^6*A*Z^2; KQ12=D43*A^2*Z^2 + D44*Ri^3*A^2*Z^2 + D45*Jj^3*Ri^3*A^2*Z^2; KQ13=D46*Ri^6*A^2*Z^2 + D47*Jj*Ri^6*A^2*Z^2; % % KQ POLYNOMIAL % KQ=KQ1+KQ2+KQ3+KQ4+KQ5+KQ6+KQ7+KQ8+KQ9+KQ10+KQ11+KQ12+KQ13; % ETA0=Jj*KT/(6.283185*KQ); if KT >=0 KTplot(j) = KT; KQplot(j) = KQ; ETAOplot(j) = ETA0; Jplot(j) = Jj; end end plot(Jplot,KTplot,'k' plot(Jplot,KTplot, 'k',Jplot,5*KQplot, ,Jplot,5*KQplot,'--' '--',Jplot,ETAOplot, ,Jplot,ETAOplot,'r' 'r'); ); axis([0 1.6 0 1.0]) legend('\fontsize{15}\bf{K_T}' legend('\fontsize{15}\bf{K_T}', ,'5 x K_Q', K_Q','\eta_0' '\eta_0',2) ,2) text(0.03,5*KQplot(4),['\fontsize{8}\bf{ text(0.03,5*KQplot(4),[ '\fontsize{8}\bf{ }',num2str(Ri)], }',num2str(Ri)],'color' 'color',[0 ,[0 0 1]); text(Jplot(length(KTplot)-3),0.01,... text(Jplot(length(KTplot)-3),0.01, ... ['\fontsize{8}\bf{ }',num2str(Ri)], }',num2str(Ri)],'color' 'color',[0 ,[0 0 0]); text(Jplot(length(KTplot)-30),ETAOplot(length(KTplot)-30),... text(Jplot(length(KTplot)-30),ETAOplot(length(KTplot)-30), ... ['\fontsize{8}\bf{ }',num2str(Ri)], }',num2str(Ri)],'color' 'color',[1 ,[1 0 0]); clear KTplot KQplot ETAOplot Jplot if pauseN == 1 pause end end Jplot = 0:0.005:1.6; if TYPE == 1 % KT/J^2 CTN=COEFF; plot(Jplot,CTN*Jplot.^2,'m' plot(Jplot,CTN*Jplot.^2, 'm', ,'linewidth' 'linewidth',[2]) ,[2]) text(0.8,CTN*0.8^2-0.02,['\fontsize{12}\bf{K_T text(0.8,CTN*0.8^2-0.02,[ '\fontsize{12}\bf{K_T = }', }',... num2str(CTN),'J^2' num2str(CTN),'J^2'], ],'color' 'color',[1 ,[1 0 1]); elseif TYPE == 2 % KQ/J^3


Page 8 May May 7, 7, 2007 2007

CQN=COEFF; plot(Jplot,5*CQN*Jplot.^3,'m' plot(Jplot,5*CQN*Jplot.^3, 'm', ,'linewidth' 'linewidth',[2]) ,[2]) text(0.8,5*CQN*0.8^3-0.02,['\fontsize{12}\bf{K_Q text(0.8,5*CQN*0.8^3-0.02,[ '\fontsize{12}\bf{K_Q = }', }',... num2str(CQN),'J^3' num2str(CQN),'J^3'], ],'color' 'color',[1 ,[1 0 1]); elseif TYPE == 3 % KT/J^4 CTN=COEFF; plot(Jplot,CTN*Jplot.^4,'m' plot(Jplot,CTN*Jplot.^4, 'm', ,'linewidth' 'linewidth',[2]) ,[2]) text(0.8,CTN*0.8^4-0.02,['\fontsize{12}\bf{K_T text(0.8,CTN*0.8^4-0.02,[ '\fontsize{12}\bf{K_T = }', }',... num2str(CTN),'J^4' num2str(CTN),'J^4'], ],'color' 'color',[1 ,[1 0 1]); elseif TYPE == 4 % KQ/J^5 CQN=COEFF; plot(Jplot,5*CQN*Jplot.^5,'m' plot(Jplot,5*CQN*Jplot.^5, 'm', ,'linewidth' 'linewidth',[2]) ,[2]) text(0.8,5*CQN*0.8^5-0.02,['\fontsize{12}\bf{K_Q text(0.8,5*CQN*0.8^5-0.02,[ '\fontsize{12}\bf{K_Q = }', }',... num2str(CQN),'J^5' num2str(CQN),'J^5'], ],'color' 'color',[1 ,[1 0 1]); end %---------------------------------------------------------------------prompt = {'Enter initial guess for J (recommended 0.8)'}; 0.8)'}; datap5 = inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defp5); Optimization',1,defp5); if isempty(datap5)==1 %clear return end if pauseN == 1 pause end GUESS = str2num(char(datap5(1))); iR = 0; for R=0.5:0.1:1.4, %P/D ratio iR = iR+1; PtoDR(iR)= R; Y = fzero('fcn' fzero('fcn',GUESS); ,GUESS); JR(iR) = Y; ETA0R(iR)= ETA0; KTR(iR) = Kt; KQR(iR) = Kq; if TYPE == 1 % KT/J^2 K = Kt; elseif TYPE == 2 % KQ/J^3 K = 5*Kq; elseif TYPE == 3 % KT/J^4 K = Kt; elseif TYPE == 4 % KQ/J^5 K = 5*Kq; end plot(Y,K,'^m' plot(Y,K,'^m', ,'linewidth' 'linewidth',[1.6]) ,[1.6]) if pauseN == 1 pause end line([Y Y], [K ETA0],'color' ETA0],'color',[1 ,[1 0 1]) plot(Y,ETA0,'^m' plot(Y,ETA0,'^m', ,'linewidth' 'linewidth',[1.6]) ,[1.6]) if pauseN == 1 pause end end [eta0,i] = max(ETA0R); PtoD = PtoDR(i); J = JR(i);


Page 9 May May 7, 7, 2007 2007

KT = KTR(i); KQ = KQR(i); plot(JR,ETA0R,'r' plot(JR,ETA0R,'r', ,'linewidth' 'linewidth',[2]) ,[2]) if pauseN == 1 pause end plot(J,eta0,'^' plot(J,eta0,'^', ,'linewidth' 'linewidth',[2]) ,[2]) text(J,eta0+0.02,['\fontsize{12}\bf{P/D text(J,eta0+0.02,[ '\fontsize{12}\bf{P/D = }' }',num2str(PtoD)], ,num2str(PtoD)],'color' 'color',[0 ,[0 0 1]); if pauseN == 1 pause end line([0 J], [eta0 eta0],'color' eta0],'color',[0 ,[0 0 1],'linewidth' 1],'linewidth',[2], ,[2], 'marker' 'marker', ,'<' '<') ) text(-0.2,eta0,['\fontsize{12}\bf{\eta_0 text(-0.2,eta0,['\fontsize{12}\bf{\eta_0 = }' }',num2str(eta0)], ,num2str(eta0)],'color' 'color',[0 ,[0 0 1]); if pauseN == 1 pause end line([J J], [eta0 0],'color' 0],'color',[0 ,[0 0 1],'linewidth' 1],'linewidth',[2], ,[2], 'marker' 'marker', ,'v' 'v') ) text(J,-0.03,['\fontsize{12}\bf{J text(J,-0.03,['\fontsize{12}\bf{J = }' }',num2str(J)], ,num2str(J)],'color' 'color',[0 ,[0 0 1]); if pauseN == 1 pause end line([0 J], [KT KT],'color' KT],'color',[0 ,[0 0 1],'linewidth' 1],'linewidth',[2], ,[2], 'marker' 'marker', ,'<' '<') ) text(-0.2,KT,['\fontsize{12}\bf{K_T text(-0.2,KT,['\fontsize{12}\bf{K_T = }',num2str(KT)], }',num2str(KT)],'color' 'color',[0 ,[0 0 1]); if pauseN == 1 pause end line([J 1.6], [5*KQ 5*KQ],'color' 5*KQ],'color',[0 ,[0 0 1],'linewidth' 1],'linewidth',[2], ,[2], 'marker' 'marker', ,'>' '>') ) text(1.61,5*KQ,['\fontsize{12}\bf{K_Q text(1.61,5*KQ,['\fontsize{12}\bf{K_Q = }',num2str(KQ)], }',num2str(KQ)],'color' 'color',[0 ,[0 0 1]); if TYPE == 1 % KT/J^2 N = 60* VAdvance/(J*D); %rpm P = PtoD*D; check = 0; elseif TYPE == 2 % KQ/J^3 N = 60* VAdvance/(J*D); %rpm P = PtoD*D; etaDnew = etaHull*etaR*eta0; EE = abs(etaDnew-etaD); if EE > Ethreshold wtxt = ['Propulsive Efficiency did not converged', converged',... ' --- etaDnew = ',num2str(etaDnew)]; ',num2str(etaDnew)]; warn = warndlg(wtxt ,' WARNING'); WARNING'); uiwait(warn) defETAD = {num2str(etaDnew)}; promptETAD = {'Propulsive {'Propulsive Efficiency (Iterated)'}; (Iterated)'}; defETAS = {num2str(etaR),num2str(MCPat),num2str(etaS)}; close(h3) check = 1; else PBt = PEt/(etaDnew*etaS); IMCP = PBt/MCPat; wtxt = ['Propulsive Efficiency CONVERGED', CONVERGED',... ' --- Installed maximum continous power (kW) = ', ',... num2str(IMCP/1000)]; warn = warndlg(wtxt ,' WARNING'); WARNING'); uiwait(warn) check = 0; end


Page 10 May May 7, 7, 2007 2007

elseif TYPE == 3 % KT/J^4 Donew = VAdvance/(J*(N/60)); Db = Donew*(1-Dfactor); P = PtoD*Db; check = 0; elseif TYPE == 4 % KQ/J^5 Donew = VAdvance/(J*(N/60)); Db = Donew*(1-Dfactor); P = PtoD*Db; etaDnew = etaHull*etaR*eta0; EE = abs(etaDnew-etaD); if EE > Ethreshold wtxt = ['Propulsive Efficiency did not converged', converged',... ' --- etaDnew = ',num2str(etaDnew)]; ',num2str(etaDnew)]; warn = warndlg(wtxt ,' WARNING'); WARNING'); uiwait(warn) defETAD = {num2str(etaDnew)}; promptETAD = {'Propulsive Efficiency (Iterated)'}; (Iterated)'}; defETAS = {num2str(etaR),num2str(MC {num2str(etaR),num2str(MCPat),num2str(etaS)}; Pat),num2str(etaS)}; defN = {num2str(N)}; close(h3) check = 1; else PBt = PEt/(etaDnew*etaS); IMCP = PBt/MCPat; wtxt = ['Propulsive Efficiency CONVERGED', CONVERGED',... ' --- Installed maximum continous power (kW) = ', ',... num2str(IMCP/1000)]; warn = warndlg(wtxt ,' WARNING'); WARNING'); uiwait(warn) check = 0; end end %---------------------------------------------------------------------end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %---------------------------------------------------------------------RESULTP1 = [... [... eta0 % Maximum Open Water Efficiency KT % Thrust Coefficient KQ % Torque Coefficient PtoD % Pitch to Diameter J]; % Advance Coefficient %---------------------------------------------------------------------% Prediction of Performance in Service % Prediction of the ship speed and propeller rate of rotation % in service condition with the engine developing 85% of MCR %---------------------------------------------------------------------if TYPE == 2 | TYPE == 4 wtxtBURRILL = 'NEXT STEP - PREDICTION OF PERFORMANCE IN SERVICE'; SERVICE'; warn = warndlg(wtxtBURRILL ,' ,' WARNING'); WARNING'); uiwait(warn) PDservice = PBt*etaS; wservice = 1.1*w; etaHulls = (1-t)/(1-wservice); %---------------------------------------------------------------------check = 1; while check == 1 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


Page 11 May May 7, 7, 2007 2007

dataETADs = inputdlg(promptETADs, inputdlg(promptETADs,... ... 'Prediction of Performance in Service',1,defETADs); Service',1,defETADs); if isempty(dataETADs)==1 %clear return end etaDs = str2num(char(dataETADs(1 str2num(char(dataETADs(1))); ))); PEservice = PDservice*etaDs; Vservice = interp1(Pinterpdata(:,1) interp1(Pinterpdata(:,1),Pinterpdata(:,2),PEservice) ,Pinterpdata(:,2),PEservice); ; VAdvances = .5144*Vservice*(1-wservi .5144*Vservice*(1-wservice); ce); %m/s CQN = PDservice/(2*pi*RhoSW*Do PDservice/(2*pi*RhoSW*Do^2*VAdvances^3); ^2*VAdvances^3); defCOEFF = {num2str(CQN)}; prompt = {'KQ/J^3 CONFIRMATION'}; CONFIRMATION'}; dataCOEFF= inputdlg(prompt,... inputdlg(prompt,... 'Prediction of Performance in Service',1,defCOEFF); Service',1,defCOEFF); if isempty(dataCOEFF)==1 %clear return end COEFF = str2num(char(dataCOEFF(1 str2num(char(dataCOEFF(1))); ))); i = 0; h4 = figure(4); set(h4,'Position' set(h4,'Position',[1,1,1400,950]) ,[1,1,1400,950]) title(['\fontsize{12}\bf{WAGENINGEN-B title(['\fontsize{12}\bf{WAGENINGEN-B SERIES PROPELLERS, FOR }', }',... num2str(Z),' num2str(Z),' BLADES, A_E/A_O = ', ',... num2str(A),', num2str(A),', P/D = 0.50 TO 1.40']); 1.40']); xlabel('\fontsize{12}\bf{J}' xlabel('\fontsize{12}\bf{J}'); ); ylabel('\fontsize{12}\bf{K_T, ylabel('\fontsize{12}\bf{K_T, 5 x K_Q, \eta_0}'); \eta_0}'); grid axis([0 1.6 0 1.0]) set(gca,'XTick' set(gca,'XTick',[0.00;0.20;0.40;0.60;0.80;1.00;1.20;1.40;1.60]) ,[0.00;0.20;0.40;0.60;0.80;1.00;1.20;1.40;1.60]) set(gca,'XTickLabel' set(gca,'XTickLabel',[0.00;0.20;0.40;0.60;0.80;1.00;1.20;1.40;1.60]) ,[0.00;0.20;0.40;0.60;0.80;1.00;1.20;1.40;1.60]) %set(gca,'YTick',[0.00;0.50;1.25;2.50;4.00;6.00]) %set(gca,'YTickLabel',[0.00;0.50;1.25;2.50;4.00;6.00]) hold on %pause for Ri = 0.5:0.1:1.4 i = i+1; j = 0; for Jj=0:0.001:1.6; j = j+1; % KT1=C1 + C2*Jj + C3*Ri + C4*Ri^2 + C5*Jj^2*A + C6*Jj*Ri*A; KT2=C7*Ri^2*A + C8*Z + C9*Jj^2*Z + C10*Ri*Z + C11*Jj*Ri*Z; KT3=C12*A*Z+C13*Jj*A*Z+C14*Ri^3+C15*Ri^6 + C16*Jj^2*Ri^6; KT4=C17*Jj^3*A+C18*A^2+C19*Jj^2*A^2 + C20*Jj^3*A^2; KT5=C21*Jj*Ri^6*A^2 + C22*Jj^2*Ri^6*A^2 + C23*Ri^3*Z; KT6=C24*Jj*Ri^3*Z + C25*Jj^3*Ri^3*Z + C26*Ri^3*A*Z; KT7=C27*Jj*A^2*Z + C28*Ri^2*A^2*Z + C29*Z^2 + C30*Jj*Z^2; KT8=C31*Jj^2*Z^2 + C32*Jj^3*Z^2 + C33*Jj*Ri^2*Z^2; KT9=C34*Jj*Ri^6*Z^2 + C35*Jj^2*Ri^6*Z^2 + C36*A*Z^2; KT10=C37*Ri^3*A*Z^2+C38*Jj^3*Ri^6*A*Z^2+C39*Ri^3*A^2*Z^2; % % KT POLYNOMIAL % KT = KT1+KT2+KT3+KT4+KT5+KT6+KT7+KT8+KT9+KT10; %


Page 12 May May 7, 7, 2007 2007

KQ1 =D1 + D2*Jj^2 + D3*Jj*Ri + D4*Ri^2 + D5*Ri*A + D6*Jj*Ri*A; KQ2 =D7*(Jj^2*Ri*A) + D8*(Ri^2*A) + D9*Jj*Z+D10*Ri*Z + D11*Jj*Ri*Z; KQ3 =D12*Jj^2*Ri*Z + D13*Jj^2*A*Z + D14*Jj*Ri*A*Z + D15*Ri^2*A*Z; KQ4 =D16*Jj^3*A + D17*Ri^3*A + D18*Jj*Ri^3*A + D19*A^2; KQ5 =D20*Jj*A^2 + D21*Jj^3*A^2 + D22*Ri*A^2 + D23*Jj^3*Ri*A^2; KQ6 =D24*Jj^2*Ri^2*A^2 + D25*Ri^3*A^2 + D26*Ri^6*A^2; KQ7 =D27*Jj^3*Z + D28*Jj^3*Ri^3*Z + D29*Ri^6*Z + D30*Jj^3*A*Z; KQ8 =D31*Ri^6*A*Z +D32*Jj*A^2*Z + D33*Ri^2*A^2*Z; KQ9 =D34*Jj^2*Ri^3*A^2*Z + D35*Ri^6*A^2*Z + D36*Jj*Ri*Z^2; KQ10=D37*Jj^3*Ri^2*Z^2 + D38*Jj^3*Ri^6*Z^2 + D39*Jj*A*Z^2; KQ11=D40*Jj^2*A*Z^2 + D41*Ri^2*A*Z^2 + D42*Ri^6*A*Z^2; KQ12=D43*A^2*Z^2 + D44*Ri^3*A^2*Z^2 + D45*Jj^3*Ri^3*A^2*Z^2; KQ13=D46*Ri^6*A^2*Z^2 + D47*Jj*Ri^6*A^2*Z^2; % % KQ POLYNOMIAL % KQ=KQ1+KQ2+KQ3+KQ4+KQ5+KQ6+KQ7+KQ8+KQ9+KQ10+KQ11+KQ12+KQ13; % ETA0=Jj*KT/(6.283185*KQ); if KT >=0 KTplot(j) = KT; KQplot(j) = KQ; ETAOplot(j) = ETA0; Jplot(j) = Jj; end end plot(Jplot,KTplot,'k' plot(Jplot,KTplot, 'k',Jplot,5*KQplot, ,Jplot,5*KQplot,'--' '--',Jplot,ETAOplot, ,Jplot,ETAOplot,'r' 'r'); ); axis([0 1.6 0 1.0]) legend('\fontsize{15}\bf{K_T}' legend('\fontsize{15}\bf{K_T}', ,'5 x K_Q', K_Q','\eta_0' '\eta_0',2) ,2) text(0.03,5*KQplot(4),['\fontsize{8}\bf{ text(0.03,5*KQplot(4),[ '\fontsize{8}\bf{ }',num2str(Ri)], }',num2str(Ri)],'color' 'color',[0 ,[0 0 1]); text(Jplot(length(KTplot)-3),0.01,['\fontsize{8}\bf{ text(Jplot(length(KTplot)-3),0.01,[ '\fontsize{8}\bf{ }', }',... num2str(Ri)],'color' num2str(Ri)],'color',[0 ,[0 0 0]); text(Jplot(length(KTplot)-30),ETAOplot(length(KTplot)-30),... text(Jplot(length(KTplot)-30),ETAOplot(length(KTplot)-30), ... ['\fontsize{8}\bf{ }',num2str(Ri)], }',num2str(Ri)],'color' 'color',[1 ,[1 0 0]); clear KTplot KQplot ETAOplot Jplot %pause end Jplot = 0:0.005:1.6; CQN =COEFF; plot(Jplot,5*CQN*Jplot.^3,'m' plot(Jplot,5*CQN*Jplot.^3, 'm', ,'linewidth' 'linewidth',[2]) ,[2]) text(0.8,5*CQN*0.8^3-0.02,['\fontsize{12}\bf{K_Q text(0.8,5*CQN*0.8^3-0.02,[ '\fontsize{12}\bf{K_Q = }', }',... num2str(CQN),'J^3' num2str(CQN),'J^3'], ],'color' 'color',[1 ,[1 0 1]); %---------------------------------------------------------------------prompt = {'Enter initial guess for J (recommended 0.8)'}; 0.8)'}; datap5 = inputdlg(prompt,'Propeller inputdlg(prompt,'Propeller B-Series Optimization',1,defp5); Optimization',1,defp5); if isempty(datap5)==1 %clear return end %pause GUESS = str2num(char(datap5(1))); iR=0; for R=0.5:0.1:1.4, %P/D ratio iR = iR+1; PtoDR(iR)= R; Y = fzero('fcn' fzero('fcn',GUESS); ,GUESS); JR(iR) = Y;


Page 13 May May 7, 7, 2007 2007

ETA0R(iR)= ETA0; KTR(iR) = Kt; KQR(iR) = Kq; K = 5*Kq; plot(Y,K,'^m' plot(Y,K,'^m', ,'linewidth' 'linewidth',[1.6]) ,[1.6]) %pause line([Y Y], [K ETA0],'color' ETA0],'color',[1 ,[1 0 1]) plot(Y,ETA0,'^m' plot(Y,ETA0,'^m', ,'linewidth' 'linewidth',[1.6]) ,[1.6]) %pause end [eta0,i] = max(ETA0R); PtoD = PtoDR(i); J = JR(i); KT = KTR(i); KQ = KQR(i); plot(JR,ETA0R,'r' plot(JR,ETA0R,'r', ,'linewidth' 'linewidth',[2]) ,[2]) %pause plot(J,eta0,'^' plot(J,eta0,'^', ,'linewidth' 'linewidth',[2]) ,[2]) text(J,eta0+0.02,['\fontsize{12}\bf{P/D text(J,eta0+0.02,[ '\fontsize{12}\bf{P/D = }' }',num2str(PtoD)], ,num2str(PtoD)],'color' 'color',[0 ,[0 0 1]); %pause line([0 J], [eta0 eta0],'color' eta0],'color',[0 ,[0 0 1],'linewidth' 1],'linewidth',[2], ,[2], 'marker' 'marker', ,'<' '<') ) text(-0.2,eta0,['\fontsize{12}\bf{\eta_0 text(-0.2,eta0,['\fontsize{12}\bf{\eta_0 = }' }', ,... num2str(eta0)],'color' num2str(eta0)],'color',[0 ,[0 0 1]); %pause line([J J], [eta0 0],'color' 0],'color',[0 ,[0 0 1],'linewidth' 1],'linewidth',[2], ,[2], 'marker' 'marker', ,'v' 'v') ) text(J,-0.03,['\fontsize{12}\bf{J text(J,-0.03,['\fontsize{12}\bf{J = }' }', ,... num2str(J)],'color' num2str(J)],'color',[0 ,[0 0 1]); %pause line([0 J], [KT KT],'color' KT],'color',[0 ,[0 0 1],'linewidth' 1],'linewidth',[2], ,[2], 'marker' 'marker', ,'<' '<') ) text(-0.2,KT,['\fontsize{12}\bf{K_T text(-0.2,KT,['\fontsize{12}\bf{K_T = }', }',... num2str(KT)],'color' num2str(KT)],'color',[0 ,[0 0 1]); %pause line([J 1.6], [5*KQ 5*KQ],'color' 5*KQ],'color',[0 ,[0 0 1],'linewidth' 1],'linewidth',[2], ,[2], 'marker' 'marker', ,'>' '>') ) text(1.61,5*KQ,['\fontsize{12}\bf{K_Q text(1.61,5*KQ,['\fontsize{12}\bf{K_Q = }', }',... num2str(KQ)],'color' num2str(KQ)],'color',[0 ,[0 0 1]); Nservice = 60* VAdvances/(J*D); %rpm etaDnews = etaHulls*etaR*eta0; EE = abs(etaDnews-etaDs); if EE > Ethreshold wtxt = ['Propulsive Efficiency did not converged', converged',... ' --- etaDnew = ',num2str(etaDnews)]; ',num2str(etaDnews)]; warn = warndlg(wtxt ,' WARNING'); WARNING'); uiwait(warn) defETADs = {num2str(etaDnews)}; promptETADs = {'Propulsive {'Propulsive Efficiency (Iterated)'}; (Iterated)'}; close(h4) check = 1; else PBt = PEt/(etaDnews*etaS); IMCP = PBt/MCPat; wtxt = ['Propulsive Efficiency CONVERGED', CONVERGED',... ' --- Propeller Rotation Rate in Service Condition (Nservice) rpm = '... ,num2str(Nservice)]; warn = warndlg(wtxt ,' WARNING'); WARNING'); uiwait(warn) check = 0; end


Page 14 May May 7, 7, 2007 2007

end h5 = figure(5); plot(Vplot,PEtplot/1000,... plot(Vplot,PEtplot/1000, ... Vplot,PEsplot/1000,'--' Vplot,PEsplot/1000, '--', ,... 'linewidth',[2]) 'linewidth' ,[2]) hold on line([Vservice Vservice], [0 PEservice/1000],... PEservice/1000],... 'color',[0 'color' ,[0 1 0],'linewidth' 0],'linewidth',[2], ,[2],'marker' 'marker', ,'v' 'v') ) line([Vservice 0], [PEservice/1000 PEservice/1000],... PEservice/1000],... 'color',[0 'color' ,[0 1 0],'linewidth' 0],'linewidth',[2]) ,[2]) text(0.2,PEservice/1000+800,['\fontsize{10}\bf{ text(0.2,PEservice/1000+800,[ '\fontsize{10}\bf{ }', }',... num2str(PEservice/1000),' num2str(PEservice/1000), ' kW (Service)']); (Service)']); text(Vservice+0.2,800,['\fontsize{10}\bf{ text(Vservice+0.2,800,[ '\fontsize{10}\bf{ }', }',... num2str(Vservice),' num2str(Vservice), ' knots (Service)']); (Service)']); title('\fontsize{15}\bf{POWER}' title('\fontsize{15}\bf{POWER}'); ); xlabel('\fontsize{12}\bf{V xlabel('\fontsize{12}\bf{V (knot)}'); (knot)}'); ylabel('\fontsize{12}\bf{P_E ylabel('\fontsize{12}\bf{P_E (kW)}'); (kW)}'); legend('\fontsize{15}\bf{Trial legend('\fontsize{15}\bf{Trial Power}', Power}', 'Service Power',2) Power',2) grid set(h5,'Position' set(h5,'Position',[1,1,1400,950]) ,[1,1,1400,950]) %---------------------------------------------------------------------end %---------------------------------------------------------------------% BURRILL CAVITATION CONTROL %---------------------------------------------------------------------wtxtBURRILL = 'NEXT STEP - BURRILL CAVITATION CONTROL'; CONTROL'; warn = warndlg(wtxtBURRILL ,' ,' WARNING'); WARNING'); uiwait(warn) % 50% back cavitation burrilldata1 = [... [... 0.0959401 0.074301914 0.1300000 0.113240036 0.1600000 0.146892628 0.2000000 0.185780446 0.2500000 0.227509743 0.3000000 0.260615355 0.4000000 0.323593657 0.5000000 0.371535229 0.5767665 0.408319386]; % 20% back cavitation burrilldata2 = [... [... 0.1056818 0.072110748 0.1300000 0.097498964 0.1600000 0.129121927 0.2000000 0.165958691 0.2500000 0.203235701 0.3000000 0.236047823 0.4000000 0.287078058 0.5000000 0.329609712 0.6000000 0.363915036 0.7000000 0.39536662]; % 10% back cavitation burrilldata3 = [... [... 0.1327394 0.072110748 0.1600000 0.092469817 0.2000000 0.12189896 0.2500000 0.152756606

propeller.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. 0.3000000 0.177827941 0.4000000 0.216271852 0.5000000 0.248313311 0.6000000 0.276057786 0.7000000 0.299226464 0.8000000 0.317687407 0.9000000 0.335737614 1.0000000 0.352370871 1.1000000 0.368977599 1.2331048 0.387257645]; % 5% back cavitation burrilldata4 = [... [... 0.1559553 0.072110748 0.1800000 0.086099375 0.2000000 0.095940063 0.2500000 0.118032064 0.3000000 0.138038426 0.4000000 0.167109061 0.5000000 0.191425593 0.6000000 0.212813905 0.7000000 0.23227368 0.8000000 0.247742206 0.9000000 0.262421854 1.0000000 0.27733201 1.1000000 0.289734359 1.2000000 0.302691343 1.3000000 0.314050869 1.4000000 0.324339617 1.5000000 0.335737614]; % 2.5% back cavitation burrilldata5 = [... [... 0.1800000 0.07014553 0.2000000 0.078523563 0.2500000 0.097498964 0.3000000 0.114287833 0.4000000 0.143548943 0.5000000 0.165576996 0.6000000 0.182810022 0.7000000 0.197696964 0.8000000 0.211348904 0.9000000 0.224905461 1.0000000 0.23659197 1.1000000 0.247742206 1.2000000 0.258226019 1.3000000 0.26915348 1.4000000 0.27733201 1.6943378 0.302691343]; % Upper limit for heavily loaded propellers (1943) % (Warship propellers with special sections) burrilldata6 = [... [... 0.1061696 0.072110748 0.1350000 0.088920112 0.1600000 0.103752842 0.2000000 0.125025903 0.2500000 0.151008015 0.3000000 0.172186857 0.3500000 0.1933657

Page 15 May May 7, 7, 2007 2007


Page 16 May May 7, 7, 2007 2007

0.4036454 0.210377844 0.4325138 0.218776162]; % Suggested upper limit (1943) % (for merchant ship propellers) burrilldata7 = [... [... 0.1241670 0.073790423 0.2000000 0.10280163 0.3000000 0.133352143 0.4000000 0.159587915 0.5000000 0.181970086 0.6000000 0.201372425 0.7000000 0.220800473 0.8000000 0.237137371 0.9000000 0.254683025 1.0000000 0.270395836 1.1000000 0.283791903 1.2000000 0.297851643 1.3000000 0.311171634 1.4000000 0.324339617 1.5000000 0.335737614]; % Suggested lower limit (1943) % (for tugs, trawlers, etc) burrilldata8 = [... [... 0.3111716 0.116412603 0.4000000 0.133659552 0.5000000 0.150314197 0.6000000 0.162554876 0.7000000 0.173780083 0.8000000 0.1840772 0.9000000 0.193196832 1.0000000 0.201372425 1.1000000 0.208929613 1.2000000 0.21677041 1.3000000 0.222843515 1.4000000 0.229086765 1.6032454 0.240990543]; %---------------------------------------------------------------------h6 = figure(6); loglog(burrilldata1(:,1),burrilldata1(:,2),'b' loglog(burrilldata1(:,1),burrilldata1(:,2), 'b', ,'linewidth' 'linewidth',[2]) ,[2]) hold on [Aburril map] = imread('burrill.JPG' imread('burrill.JPG'); ); image([0.05 2],[0.05 0.5],Aburril) loglog(burrilldata1(:,1),burrilldata1(:,2),'b' loglog(burrilldata1(:,1),burrilldata1(:,2), 'b', ,'linewidth' 'linewidth',[2]) ,[2]) loglog(burrilldata2(:,1),burrilldata2(:,2),'b' loglog(burrilldata2(:,1),burrilldata2(:,2), 'b', ,'linewidth' 'linewidth',[2]) ,[2]) loglog(burrilldata3(:,1),burrilldata3(:,2),'b' loglog(burrilldata3(:,1),burrilldata3(:,2), 'b', ,'linewidth' 'linewidth',[2]) ,[2]) loglog(burrilldata4(:,1),burrilldata4(:,2),'b' loglog(burrilldata4(:,1),burrilldata4(:,2), 'b', ,'linewidth' 'linewidth',[2]) ,[2]) loglog(burrilldata5(:,1),burrilldata5(:,2),'b' loglog(burrilldata5(:,1),burrilldata5(:,2), 'b', ,'linewidth' 'linewidth',[2]) ,[2]) loglog(burrilldata6(:,1),burrilldata6(:,2),'r--' loglog(burrilldata6(:,1),burrilldata6(:,2), 'r--', ,'linewidth' 'linewidth',[2]) ,[2]) loglog(burrilldata7(:,1),burrilldata7(:,2),'r--' loglog(burrilldata7(:,1),burrilldata7(:,2), 'r--', ,'linewidth' 'linewidth',[2]) ,[2]) loglog(burrilldata8(:,1),burrilldata8(:,2),'r--' loglog(burrilldata8(:,1),burrilldata8(:,2), 'r--', ,'linewidth' 'linewidth',[2]) ,[2]) axis([0.05 2 0.05 0.5]) axis tight set(gca,'XTick' set(gca,'XTick',[0.05;0.1;0.2;0.3;0.4;0.5;0.6;0.7;0.8;0.9;1.0;1.5;2.0]) ,[0.05;0.1;0.2;0.3;0.4;0.5;0.6;0.7;0.8;0.9;1.0;1.5;2.0]) set(gca,'XTickLabel' set(gca,'XTickLabel',[0.05;0.1;0.2;0.3;0.4;0.5;0.6;0.7;0.8;0.9;1.0;1.5;2.0]) ,[0.05;0.1;0.2;0.3;0.4;0.5;0.6;0.7;0.8;0.9;1.0;1.5;2.0]) set(gca,'YTick' set(gca,'YTick',[0.05;0.1;0.2;0.3;0.4;0.5]) ,[0.05;0.1;0.2;0.3;0.4;0.5]) set(gca,'YTickLabel' set(gca,'YTickLabel',[0.05;0.1;0.2;0.3;0.4;0.5]) ,[0.05;0.1;0.2;0.3;0.4;0.5]) xlabel('\fontsize{12}\bf{Local xlabel('\fontsize{12}\bf{Local cavitation number at 0.7R ( \sigma_{0.7} )}');


Page 17 May May 7, 7, 2007 2007

ylabel('\fontsize{12}\bf{The ylabel('\fontsize{12}\bf{The thrust loading coefficient ( \tau_C )}'); )}'); title('\fontsize{18}\bf{The title('\fontsize{18}\bf{The Burrill"s Cavitation Diagram}'); Diagram}'); set(h6,'Position' set(h6,'Position',[1,1,1400,950]) ,[1,1,1400,950]) %---------------------------------------------------------------------VR = sqrt(VAdvance^2+(2*pi*(N/6 sqrt(VAdvance^2+(2*pi*(N/60)*0.7*(D/2))^2); 0)*0.7*(D/2))^2); % Resultant velocity at r = 0.7*R qT = 0.5*RhoSW*VR^2; % Dynamic pressure at r = 0.7*R h = (D/2)+0.2; % height of shaft centre-line above base Patm = 101300; % (N/m^2) Atmospheric pressure Pv = 1646; % (N/m^2) Vapour pressure of water at 15 °C Po = Patm+RhoSW*g*(T-h); % the static pressure at the shaft centerline sigmaR = (Po-Pv)/qT; % local cavitation number at r = 0.7*R BURRILLTYPE = MENU('The MENU('The Burrill"s Cavitation Diagram', Diagram',... ' 50% back cavitation', cavitation',... ' 20% back cavitation', cavitation',... ' 10% back cavitation', cavitation',... ' 5% back cavitation' cavitation', ,... '2.5% back cavitation', cavitation',... ['Upper limit for heavily loaded propellers (1943)', (1943)',... ' (Warship propellers with special sections)'], sections)'],... ... 'Suggested upper limit (1943) (for merchant ship propellers)', propellers)',... 'Suggested lower limit (1943) (for tugs, trawlers, etc)', etc)',... 'CANCEL'); 'CANCEL' ); if BURRILLTYPE == 1 burrilldata = burrilldata1; elseif BURRILLTYPE == 2 burrilldata = burrilldata2; elseif BURRILLTYPE == 3 burrilldata = burrilldata3; elseif BURRILLTYPE == 4 burrilldata = burrilldata4; elseif BURRILLTYPE == 5 burrilldata = burrilldata5; elseif BURRILLTYPE == 6 burrilldata = burrilldata6; elseif BURRILLTYPE == 7 burrilldata = burrilldata7; elseif BURRILLTYPE == 8 burrilldata = burrilldata8; elseif BURRILLTYPE == 9 % CANCEL %clear return end tauC = interp1(burrilldata(:,1),b interp1(burrilldata(:,1),burrilldata(:,2),sigmaR, urrilldata(:,2),sigmaR,'spline' 'spline', ,'extrap' 'extrap'); ); Ap = Thrust/(tauC*qT); % projected area for the propeller AE = Ap/(1.067-0.229*P/D); Ap/(1.067-0.229*P/D);% % expanded area for the propeller Anew = AE/(pi*D^2/4); % This is the minimum BAR to avoid cavitation %---------------------------------------------------------------------line([sigmaR sigmaR], [0.05 tauC],'color' tauC],'color',[1 ,[1 0 1],'linewidth' 1],'linewidth',[4]) ,[4]) line([sigmaR 0.05], [tauC tauC],'color' tauC],'color',[1 ,[1 0 1],'linewidth' 1],'linewidth',[4], ,[4],'marker' 'marker', ,'<' '<') ) %---------------------------------------------------------------------if Anew < A wtxtBURRILL = ['Calculated BAR (', (', num2str(Anew),... num2str(Anew),... ') is LESS than the selected BAR (', (',... num2str(A), '). THE DESIGN STAGE IS COMPLETED!']; COMPLETED!']; else wtxtBURRILL = ['Calculated BAR (', (', num2str(Anew),... num2str(Anew),...


Page 18 May May 7, 7, 2007 2007

') is GREATER than the selected BAR (',num2str(A), (',num2str(A),... ... '). Choose a new BAR greater than the calculated BAR and ', ',... ' perform all calculations for the new bar']; bar']; end warn = warndlg(wtxtBURRILL ,' ,' WARNING'); WARNING'); uiwait(warn) %---------------------------------------------------------------------% EXCEL OUTPUT %---------------------------------------------------------------------question = questdlg('EXCEL questdlg('EXCEL OUTPUT ?', ?',... ' KIVANC', KIVANC','YES' 'YES', ,'NO' 'NO', ,'NO' 'NO'); ); switch question case 'YES' exceloutput = 1; case 'NO' exceloutput = 2; end if exceloutput == 1 Excel = actxserver('Excel.Application' actxserver('Excel.Application'); ); set(Excel, 'Visible' , 1); Workbook = invoke(Excel.Workbooks, 'Add' 'Add'); ); %-----------------------------------------------------values1 = [... NoP % Number of propellers Z % Number of blades w % Wake Fraction t % Thrust Deduction Fraction D % Propeller Diameter (initial) Do % Open Water Diameter A % AE/A0 ratio Vt % Trial Speed from powerpredict.m (knots) RTVt % Total Ship Resistance for the Trial Speed (N) Thrust % Thrust (Newton) VAdvance % Advance Velocity (m/s) etaHull]; % Hull efficiency %-----------------------------------------------------text1 = [... cellstr('PROPELLER cellstr('PROPELLER DESIGN'); DESIGN');... ... cellstr('Number cellstr('Number of propellers - INPUT'); INPUT');... ... cellstr('Number cellstr('Number of blades - INPUT'); INPUT');... ... cellstr('Wake cellstr('Wake Fraction - INPUT OR ESTIMATED'); ESTIMATED');... ... cellstr('Thrust cellstr('Thrust Deduction Fraction - INPUT OR ESTIMATED'); ESTIMATED');... ... cellstr('Propeller cellstr('Propeller Diameter in meters(initial) - INPUT OR ESTIMATED'); ESTIMATED');... ... cellstr('Open cellstr('Open Water Diameter - CALCULATED'); CALCULATED');... ... cellstr('AE/A0 cellstr('AE/A0 ratio (BAR) - INPUT OR ESTIMATED USING KELLER'); KELLER');... ... cellstr('Trial cellstr('Trial Speed (knots) from powerpredict.m'); powerpredict.m');... ... cellstr('Total cellstr('Total Ship Resistance (N) for the Trial Speed from powerpredict.m'); powerpredict.m');... ... cellstr('Thrust cellstr('Thrust for each propeller (Newton) - CALCULATED'); CALCULATED');... ... cellstr('Advance cellstr('Advance Velocity (m/s) - CALCULATED'); CALCULATED');... ... cellstr('Hull cellstr('Hull efficiency - CALCULATED')]; CALCULATED')]; burlist = {... {... ' 50% back cavitation', cavitation',... ' 20% back cavitation', cavitation',... ' 10% back cavitation', cavitation',... ' 5% back cavitation' cavitation', ,... '2.5% back cavitation', cavitation',... ['Upper limit for heavily loaded propellers (1943)', (1943)',...


Page 19 May May 7, 7, 2007 2007

' (Warship propellers with special sections)'], sections)'],... ... 'Suggested upper limit (1943) (for merchant ship propellers)', propellers)',... 'Suggested lower limit (1943) (for tugs, trawlers, etc)'}; etc)'}; %-----------------------------------------------------if TYPE == 1 % KT/J^2 %-----------------------------------------------------values2 = [... COEFF % KT/J^2 RESULTP1 N % Rotation Rate (rpm), N P]; % Propeller Pitch (m) %-----------------------------------------------------text2 = [... cellstr(['TYPE=1 cellstr(['TYPE=1 : KT/J^2 case - Known thrust (T), diameter (D)', (D)',... ' and advanced velocity (VA)- Nopt (required)']); (required)']);... ... cellstr('KT/J^2 cellstr('KT/J^2 - CALCULATED'); CALCULATED');... ... cellstr('Maximum cellstr('Maximum Open Water Efficiency - OUTPUT'); OUTPUT');... ... cellstr('Thrust cellstr('Thrust Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Torque cellstr('Torque Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Pitch cellstr('Pitch to Diameter - OUTPUT'); OUTPUT');... ... cellstr('Advance cellstr('Advance Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Rotation cellstr('Rotation Rate (rpm), N - OUTPUT'); OUTPUT');... ... cellstr('Propeller cellstr('Propeller Pitch (m) - OUTPUT'); OUTPUT');... ... cellstr(['BURRILL cellstr(['BURRILL CAVITATION CONTROL :', :',... char(burlist(BURRILLTYPE))]);... char(burlist(BURRILLTYPE))]); ... cellstr(wtxtBURRILL)]; %-----------------------------------------------------elseif TYPE == 2 % KQ/J^3 %-----------------------------------------------------values2 = [... etaR % Relative-rotative Efficiency MCPat % Maximum continous power at = etaS % shaft transmission efficiency etaD % Propulsive efficiency PEt/1000 % Effective Power for the Trial Speed (kW) PD/1000 % Delivered power (kW) COEFF % KQ/J^3 RESULTP1 N % Rotation Rate (rpm), N P % Propeller Pitch (m) PBt/1000 % Brake power (kW) IMCP/1000 % Installed maximum continous power (kW) etaDs % Propulsive efficiency for the service condition PEservice/1000 % Effective Power for the Service Condition (kW) Vservice % Ship speed (knots) in service condition Nservice]; % Propeller rate of rotation (rpm) in service condition %-----------------------------------------------------text2 = [... cellstr(['TYPE=2 cellstr(['TYPE=2 : KQ/J^3 case - Known power (PD), diameter (D)', (D)',... ' and advanced velocity (VA)- Nopt (required)']); (required)']);... ... cellstr('Relative-rotative cellstr('Relative-rotative Efficiency - INPUT'); INPUT');... ... cellstr('Maximum cellstr('Maximum continous power at - INPUT'); INPUT');... ... cellstr('Shaft cellstr('Shaft transmission efficiency - INPUT'); INPUT');... ... cellstr('Propulsive cellstr('Propulsive efficiency - CONVERGED'); CONVERGED');... ... cellstr('Effective cellstr('Effective Power for the Trial Speed (kW) from powerpredict.m'); powerpredict.m');... ... cellstr('Delivered cellstr('Delivered power (kW) - CALCULATED'); CALCULATED');... ... cellstr('KQ/J^3 cellstr('KQ/J^3 - CALCULATED'); CALCULATED');... ...


Page 20 May May 7, 7, 2007 2007

cellstr('Maximum cellstr('Maximum Open Water Efficiency - OUTPUT'); OUTPUT');... ... cellstr('Thrust cellstr('Thrust Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Torque cellstr('Torque Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Pitch cellstr('Pitch to Diameter - OUTPUT'); OUTPUT');... ... cellstr('Advance cellstr('Advance Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Rotation cellstr('Rotation Rate (rpm), N - OUTPUT'); OUTPUT');... ... cellstr('Propeller cellstr('Propeller Pitch (m) - OUTPUT'); OUTPUT');... ... cellstr('Brake cellstr('Brake power (kW) - OUTPUT'); OUTPUT');... ... cellstr('Installed cellstr('Installed maximum continous power (kW) - OUTPUT'); OUTPUT');... ... cellstr('Propulsive cellstr('Propulsive efficiency for the service condition - CONVERGED '); ');... ... cellstr('Effective cellstr('Effective Power for the Service Condition (kW) - CALCULATED'); CALCULATED');... ... cellstr('Ship cellstr('Ship speed (knots) in service condition - OUTPUT'); OUTPUT');... ... cellstr('Propeller cellstr('Propeller rate of rotation (rpm) in service condition - OUTPUT'); OUTPUT');... ... cellstr(['BURRILL cellstr(['BURRILL CAVITATION CONTROL :', :',... char(burlist(BURRILLTYPE))]);... char(burlist(BURRILLTYPE))]); ... cellstr(wtxtBURRILL)]; %-----------------------------------------------------elseif TYPE == 3 % KT/J^4 %-----------------------------------------------------values2 = [... N % Rotation Rate (rpm), N COEFF % KT/J^4 RESULTP1 Db % Propeller Diameter in meters P]; % Propeller Pitch (m) %-----------------------------------------------------text2 = [... cellstr(['TYPE=3 cellstr(['TYPE=3 : KT/J^4 case - Known thrust (T), rotation rate (N)', (N)',... ' and advanced velocity (VA)- Dopt (required)']); (required)']);... ... cellstr('Rotation cellstr('Rotation Rate (rpm), N - INPUT'); INPUT');... ... cellstr('KT/J^4 cellstr('KT/J^4 - CALCULATED'); CALCULATED');... ... cellstr('Maximum cellstr('Maximum Open Water Efficiency - OUTPUT'); OUTPUT');... ... cellstr('Thrust cellstr('Thrust Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Torque cellstr('Torque Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Pitch cellstr('Pitch to Diameter - OUTPUT'); OUTPUT');... ... cellstr('Advance cellstr('Advance Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Propeller cellstr('Propeller Diameter in meters - OUTPUT'); OUTPUT');... ... cellstr('Propeller cellstr('Propeller Pitch (m) - OUTPUT'); OUTPUT');... ... cellstr(['BURRILL cellstr(['BURRILL CAVITATION CONTROL :', :',... char(burlist(BURRILLTYPE))]);... char(burlist(BURRILLTYPE))]); ... cellstr(wtxtBURRILL)]; %-----------------------------------------------------elseif TYPE == 4 % KQ/J^5 %-----------------------------------------------------values2 = [... N % Rotation Rate (rpm), N etaR % Relative-rotative Efficiency MCPat % Maximum continous power at = etaS % shaft transmission efficiency etaD % Propulsive efficiency PEt/1000 % Effective Power for the Trial Speed (kW) PD/1000 % Delivered power (kW) COEFF % KQ/J^5 RESULTP1 Db % Propeller Diameter in meters P % Propeller Pitch (m) PBt/1000 % Brake power (kW)


Page 21 May May 7, 7, 2007 2007

IMCP/1000 % Installed maximum continous power (kW) etaDs % Propulsive efficiency for the service condition PEservice/1000 % Effective Power for the Service Condition (kW) Vservice % Ship speed (knots) in service condition Nservice]; % Propeller rate of rotation (rpm) in service condition %-----------------------------------------------------text2 = [... cellstr(['TYPE=4 cellstr(['TYPE=4 : KQ/J^5 case - Known power (PD), rotation rate (N)', (N)',... ' and advanced velocity (VA)- Dopt (required)']); (required)']);... ... cellstr('Rotation cellstr('Rotation Rate (rpm), N - INPUT'); INPUT');... ... cellstr('Relative-rotative cellstr('Relative-rotative Efficiency - INPUT'); INPUT');... ... cellstr('Maximum cellstr('Maximum continous power at - INPUT'); INPUT');... ... cellstr('Shaft cellstr('Shaft transmission efficiency - INPUT'); INPUT');... ... cellstr('Propulsive cellstr('Propulsive efficiency - CONVERGED'); CONVERGED');... ... cellstr('Effective cellstr('Effective Power for the Trial Speed (kW) from powerpredict.m'); powerpredict.m');... ... cellstr('Delivered cellstr('Delivered power (kW) - CALCULATED'); CALCULATED');... ... cellstr('KQ/J^5 cellstr('KQ/J^5 - CALCULATED'); CALCULATED');... ... cellstr('Maximum cellstr('Maximum Open Water Efficiency - OUTPUT'); OUTPUT');... ... cellstr('Thrust cellstr('Thrust Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Torque cellstr('Torque Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Pitch cellstr('Pitch to Diameter - OUTPUT'); OUTPUT');... ... cellstr('Advance cellstr('Advance Coefficient - OUTPUT'); OUTPUT');... ... cellstr('Propeller cellstr('Propeller Diameter in meters - OUTPUT'); OUTPUT');... ... cellstr('Propeller cellstr('Propeller Pitch (m) - OUTPUT'); OUTPUT');... ... cellstr('Brake cellstr('Brake power (kW) - OUTPUT'); OUTPUT');... ... cellstr('Installed cellstr('Installed maximum continous power (kW) - OUTPUT'); OUTPUT');... ... cellstr('Propulsive cellstr('Propulsive efficiency for the service condition - CONVERGED '); ');... ... cellstr('Effective cellstr('Effective Power for the Service Condition (kW) - CALCULATED'); CALCULATED');... ... cellstr('Ship cellstr('Ship speed (knots) in service condition - OUTPUT'); OUTPUT');... ... cellstr('Propeller cellstr('Propeller rate of rotation (rpm) in service condition - OUTPUT'); OUTPUT');... ... cellstr(['BURRILL cellstr(['BURRILL CAVITATION CONTROL :', :',... char(burlist(BURRILLTYPE))]);... char(burlist(BURRILLTYPE))]); ... cellstr(wtxtBURRILL)]; %-----------------------------------------------------end %-----------------------------------------------------space = 2; Ltext1 = length(text1); firstcell = ['A' num2str(space)]; lastcell = ['A' num2str(space+Ltext1-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet,'Range' get(Activesheet,'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , text1, 'ColumnWidth' 'ColumnWidth',87); ,87); %-----------------------------------------------------space = space+Ltext1; Ltext2 = length(text2); firstcell = ['A' num2str(space)]; lastcell = ['A' num2str(space+Ltext2-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet,'Range' get(Activesheet,'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , text2, 'ColumnWidth' 'ColumnWidth',87); ,87); %-----------------------------------------------------space = 3; Lvalues1 = length(values1); firstcell = ['B' num2str(space)]; lastcell = ['B' num2str(space+Lvalues1-1)]; Activesheet = Excel.Activesheet;

propeller.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. ActivesheetRange= get(Activesheet,'Range' get(Activesheet,'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , values1, 'ColumnWidth' 'ColumnWidth',21); ,21); %-----------------------------------------------------space = space+Lvalues1+1; Lvalues2 = length(values2); firstcell = ['B' num2str(space)]; lastcell = ['B' num2str(space+Lvalues2-1)]; Activesheet = Excel.Activesheet; ActivesheetRange= get(Activesheet,'Range' get(Activesheet,'Range',firstcell,lastcell); ,firstcell,lastcell); set(ActivesheetRange, 'Value' 'Value', , values2, 'ColumnWidth' 'ColumnWidth',21); ,21); %-----------------------------------------------------prompt ={'ENTER ={'ENTER THE FILE NAME FOR THE RESULTS'}; RESULTS'}; defname = {['Project' {['Project',num2str(1), ,num2str(1),'-propeller' '-propeller']}; ]}; fname = inputdlg(prompt,' inputdlg(prompt,' FILE NAME',1,defname); NAME',1,defname); filename = char(fname); if isempty(filename)==1 clc; return end w=pwd; invoke(Workbook, 'SaveAs' 'SaveAs', , [w,'\' [w,'\',filename, ,filename,'.xls' '.xls']); ]); delete(Excel); end %---------------------------------------------------------------------% REFERENCES: %---------------------------------------------------------------------% [1] propopt.m (Propeller Optimization Program in MATLAB) % http://web.nps.navy.mil/~me/tsse/TS4001/docs_software.htm % [2] fcn.m (Propeller Optimization Program in MATLAB) % http://web.nps.navy.mil/~me/tsse/TS4001/docs_software.htm % [3] PARSONS,M., Propeller Optimization Program (POP), % University of Michigan, Department of Naval % Architecture and Marine Engineering, % Conceptual Ship Design Software Environment %----------------------------------------------------------------------

Page 22 May May 7, 7, 2007 2007

fcn.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. function F = fcn(X) global Z A TYPE COEFF R ETA0 Kt Kq J = X; if J <=0 J=0; end % % This function fcn.m is taken directly from % http://web.nps.navy.mil/~me/tsse/TS4001/docs_software.htm % % COEFFICIENTS FOR KT POLYNOMIAL % C1 = .0088049600; C2 =-.2045540000; C3 = .1663510000; C4 = .1581140000; C5 =-.1475810000; C6 =-.4814970000; C7 = .4154370000; C8 = .0144043000; C9 =-.0530054000; C10= .0143481000; C11= .0606826000; C12=-.0125894000; C13= .0109689000; C14=-.1336980000; C15= .0063840700; C16=-.0013271800; C17= .1684960000; C18=-.0507214000; C19= .0854559000; C20=-.0504475000; C21= .0104650000; C22=-.0064827200; C23=-.0084172800; C24= .0168424000; C25=-.0010229600; C26=-.0317791000; C27= .0186040000; C28=-.0041079800; C29=-.0006068480; C30=-.0049819000; C31= .0025983000; C32=-.0005605280; C33=-.0016365200; C34=-.0003287870; C35= .0001165020; C36= .0006909040; C37= .0042174900; C38= .0000565229; C39=-.0014656400; % % COEFFICIENTS FOR KQ POLYNOMIAL % D1 = .0037936800; D2 = .0088652300; D3 =-.0322410000;

Page 1 May May 5, 2007 2007

fcn.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. D4 = .0034477800; D5 =-.0408811000; D6 =-.1080090000; D7 =-.0885381000; D8 = .1885610000; D9 =-.0037087100; D10= .0051369600; D11= .0209449000; D12= .0047431900; D13=-.0072340800; D14= .0043838800; D15=-.0269403000; D16= .0558082000; D17= .0161886000; D18= .0031808600; D19= .0158960000; D20= .0471729000; D21= .0196283000; D22=-.0502782000; D23=-.0300550000; D24= .0417122000; D25=-.0397722000; D26=-.0035002400; D27=-.0106854000; D28= .0011090300; D29=-.0003139120; D30= .0035985000; D31=-.0014212100; D32=-.0038363700; D33= .0126803000; D34=-.0031827800; D35= .0033426800; D36=-.0018349100; D37= .0001124510; D38=-.0000297228; D39= .0002695510; D40= .0008326500; D41= .0015533400; D42= .0003026830; D43=-.0001843000; D44=-.0004253990; D45= .0000869243; D46=-.0004659000; D47= .0000554194; % KQ1 =D1 + D2*J^2 + D3*J*R + D4*R^2 + D5*R*A + D6*J*R*A; KQ2 =D7*(J^2*R*A) + D8*(R^2*A) + D9*J*Z+D10*R*Z + D11*J*R*Z; KQ3 =D12*J^2*R*Z + D13*J^2*A*Z + D14*J*R*A*Z + D15*R^2*A*Z; KQ4 =D16*J^3*A + D17*R^3*A + D18*J*R^3*A + D19*A^2; KQ5 =D20*J*A^2 + D21*J^3*A^2 + D22*R*A^2 + D23*J^3*R*A^2; KQ6 =D24*J^2*R^2*A^2 + D25*R^3*A^2 + D26*R^6*A^2; KQ7 =D27*J^3*Z + D28*J^3*R^3*Z + D29*R^6*Z + D30*J^3*A*Z; KQ8 =D31*R^6*A*Z +D32*J*A^2*Z + D33*R^2*A^2*Z; KQ9 =D34*J^2*R^3*A^2*Z + D35*R^6*A^2*Z + D36*J*R*Z^2; KQ10=D37*J^3*R^2*Z^2 + D38*J^3*R^6*Z^2 + D39*J*A*Z^2; KQ11=D40*J^2*A*Z^2 + D41*R^2*A*Z^2 + D42*R^6*A*Z^2; KQ12=D43*A^2*Z^2 + D44*R^3*A^2*Z^2 + D45*J^3*R^3*A^2*Z^2;

Page 2 May May 5, 2007 2007

fcn.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. KQ13=D46*R^6*A^2*Z^2 + D47*J*R^6*A^2*Z^2; % % KQ POLYNOMIAL % KQ=KQ1+KQ2+KQ3+KQ4+KQ5+KQ6+KQ7+KQ8+KQ9+KQ10+KQ11+KQ12+KQ13; % PQR1 =D3*J + D4*2*R +D5*A +D6*J*A + D7*J^2.*A + D8*2.*R*A; PQR2 =D10*Z + D11*J*Z + D12*J^2*Z + D14*J*A*Z+D15*2.*R*A*Z; PQR3 =D17*3.*R^2*A+D18*3.*J*R^2*A+D22*A^2+D23*J^3*A^2; PQR4 =D24*2.*J^2*R*A^2 + D25*3.*R^2*A^2 +D26*6.*R^5*A^2; PQR5 =D28*3.*J^3*R^2*Z + D29*6.*R^5*Z + D31*6.*R^5*A*Z; PQR6 =D33*2.*R*A^2*Z+D34*3.*J^2*R^2*A^2*Z+D35*6.*R^5*A^2*Z; PQR7 =D36*J*Z^2+D37*2.*J^3*R*Z^2 + D38*6.*J^3*R^5*Z^2; PQR8 =D41*2.*R*A*Z^2+D42*6*R^5*A*Z^2+D44*3.*R^2*A^2*Z^2; PQR9 =D45*3.*J^3*R^2*A^2*Z^2 + D46*6.*R^5*A^2*Z^2; PQR10=D47*6.*J*R^5*A^2*Z^2; % % PARTIAL DERIVATIVE KQ/(P/D) % PKQR=PQR1+PQR2+PQR3+PQR4+PQR5+PQR6+PQR7+PQR8+PQR9+PQR10; % PKQJ1=D2*2*J+D3*R+D6*R*A+D7*2*J*R*A+D9*Z+D11*R*Z+D12*2*J*R*Z; PKQJ2= D13*2*J*A*Z + D14*R*A*Z + D16*3*J^2*A +D18*R^3*A; PKQJ3= D20*A^2 + D21*3*J^2*A^2 + D23*3*J^2*R*A^2; PKQJ4= D24*2*J*R^2*A^2 + D27*3*J^2*Z +D28*3*J^2*R^3*Z; PKQJ5=D30*3*J^2*A*Z + D32*A^2*Z + D34*2*J*R^3*A^2*Z; PKQJ6=D36*R*Z^2+D37*3*J^2*R^2*Z^2+D38*3*J^2*R^6*Z^2; PKQJ7=D39*A*Z^2+D40*2*J*A*Z^2+D45*3*J^2*R^3*A^2*Z^2; PKQJ8= D47*R^6*A^2*Z^2; % % PARTIAL DERIVATIVE KQ/J % PKQJ= PKQJ1+PKQJ2+PKQJ3+PKQJ4+PKQJ5+PKQJ6+PKQJ7+PKQJ8; % KT1=C1 + C2*J + C3*R + C4*R^2 + C5*J^2*A + C6*J*R*A; KT2=C7*R^2*A + C8*Z + C9*J^2*Z + C10*R*Z + C11*J*R*Z; KT3=C12*A*Z+C13*J*A*Z+C14*R^3+C15*R^6 + C16*J^2*R^6; KT4=C17*J^3*A+C18*A^2+C19*J^2*A^2 + C20*J^3*A^2; KT5=C21*J*R^6*A^2 + C22*J^2*R^6*A^2 + C23*R^3*Z; KT6=C24*J*R^3*Z + C25*J^3*R^3*Z + C26*R^3*A*Z; KT7=C27*J*A^2*Z + C28*R^2*A^2*Z + C29*Z^2 + C30*J*Z^2; KT8=C31*J^2*Z^2 + C32*J^3*Z^2 + C33*J*R^2*Z^2; KT9=C34*J*R^6*Z^2 + C35*J^2*R^6*Z^2 + C36*A*Z^2; KT10=C37*R^3*A*Z^2+C38*J^3*R^6*A*Z^2+C39*R^3*A^2*Z^2; % % KT POLYNOMIAL % KT = KT1+KT2+KT3+KT4+KT5+KT6+KT7+KT8+KT9+KT10; % PKTR1 = C3 + C4*2*R + C6*J*A + C7*2*R*A + C10*Z + C11*J*Z; PKTR2 = C14*3*R^2 + C15*6*R^5 + C16*6*J^2*R^5; PKTR3 = C21*6*J*R^5*A^2 + C22*6*J^2*R^5*A^2; PKTR4 = C23*3*R^2*Z + C24*3*J*R^2*Z + C25*3*J^3*R^2*Z; PKTR5 = C26*3*R^2*A*Z + C28*2*R*A^2*Z + C33*2*J*R*Z^2; PKTR6 = C34*6*J*R^5*Z^2 + C35*6*J^2*R^5*Z^2; PKTR7 = C37*3*R^2*A*Z^2 + C38*6*J^3*R^5*A*Z^2; PKTR8 = C39*3*R^2*A^2*Z^2;

Page 3 May May 5, 2007 2007

fcn.m C:\D C:\Doc ocum umen ents ts and and Sett Settin ings gs\x \xp\ p\My My Docu Docume ment nts\ s\ZO ZOVA VANC NC\T \Tur urki kiy. y... .. % % PARTIAL DERIVATIVE KT/(P/D) % PKTR = PKTR1+PKTR2+PKTR3+PKTR4+PKTR5+PKTR6+PKTR7+PKTR8; % PKTJ1 =C2+C5*2*J*A+C6*R*A+C9*2*J*Z+C11*R*Z + C13*A*Z; PKTJ2 = C16*2*J*R^6 + C17*3*J^2*A + C19*2*J*A^2; PKTJ3 = C20*3*J^2*A^2+C21*R^6*A^2+C22*2*J*R^6*A^2; PKTJ4 = C24*R^3*Z + C25*3*J^2*R^3*Z + C27*A^2*Z; PKTJ5 = C30*Z^2 + C31*2*J*Z^2 + C32*3*J^2*Z^2; PKTJ6 = C33*R^2*Z^2 + C34*R^6*Z^2; PKTJ7 = C35*2*J*R^6*Z^2 + C38*3*J^2*R^6*A*Z^2; % % PARTIAL DERIVATIVE KT/J % PKTJ = PKTJ1+PKTJ2+PKTJ3+PKTJ4+PKTJ5+PKTJ6+PKTJ7; % % SET-UP LAGRANGE MULTIPLIER EQUATIONS FOR % THE CORRESPONDING OPTIMIZATION PROBLEM % ETA0=J*KT/(6.283185*KQ); if TYPE==1 CTN=COEFF; F1= KT-CTN*J^2; F2= PKTR*(-3.*KQ+J*PKQJ)+PKQR*(-J*PKTJ+2.*KT); elseif TYPE==2 CQN=COEFF; F1= KQ-CQN*J^3; F2= PKTR*(-3.*KQ+J*PKQJ)+PKQR*(-J*PKTJ+2.*KT); elseif TYPE==3 CTD=COEFF; F1= KT-CTD*J^4; F2= PKTR*(J*PKQJ-5*KQ) + PKQR*(4*KT-J*PKTJ); elseif TYPE==4 CQD=COEFF; F1= KQ-CQD*J^5; F2= PKTR*(J*PKQJ-5*KQ) + PKQR*(4*KT-J*PKTJ); end F=F1; Kt=KT; Kq=KQ; %---------------------------------------------------------------------% REFERENCES: % [1] fcn.m Propeller Optimization Program % http://web.nps.navy.mil/~me/tsse/TS4001/docs_software.htm

Page 4 May May 5, 2007 2007

Project 1

Recommend Documents