DO 20 J=1,7 REFF(J)=1+((J‐1)*1.5/6.0) YDEL(J)=CLR*REFF(J) YVR(J)=‐YDEL(J) YRR(J)=XRR*YDEL(J) NDEL(J)=‐XRR*YDEL(J) NVR(J)=XRR*YDEL(J) NRR(J)=‐XRR**2*YDEL(J) DO 30 K=1,7 SEFF(K)=(K‐1)/2.0 YVS(K)=‐CLS*SEFF(K) YRS(K)=‐XSS*YVS(K) NVS(K)=‐XSS*YVS(K) NRS(K)=XSS**2*YVS(K) NV(I,J,K)=NVT(I)+NVR(J)+NVS(K) NR(I,J,K)=NRT(I)+NRR(J)+NRS(K) YV(I,J,K)=YVT(I)+YVR(J)+YVS(K) YR(I,J,K)=YRT(I)+YRR(J)+YRS(K) S(I,J,K)=(NR(I,J,K)/(YR(I,J,K)‐M))‐(NV(I,J,K)/YV (I,J,K)) RAD(I,J,K)=(L*((YV(I,J,K)*NR(I,J,K))‐(NV(I,J,K)* (YR(I,J,K) $ ‐M))))/(DEL*((NV(I,J,K)*YDEL(J))‐(YV(I,J,K) *NDEL(J)))) WRITE(1,5)TR(I),REFF(J),SEFF(K),YR(I,J,K),NR (I,J,K) $ ,YV(I,J,K),NV(I,J,K),S(I,J,K),RAD(I,J,K) 5 FORMAT(1X,F5.2,2X,F5.2,2X,F5.2,X, F6.3,2X,F6.3,2X,F6.3,2X $ ,F6.3,2X,F6.3,2X,F8.3 30 CONTINUE 20 CONTINUE 10 CONTINUE STOP END <strong>MIMET</strong> Technical Bulletin Volume 1 (2) 2010 Ship’s Data The above program was run based on the following ship’s input data as shown in Table 2; Table 2: Ship’s input data Distance of rudder center from Longitudinal Centre Gravity (LCG), a 50m aft of LCG Distance of skeg center from LCG, b 45m aft of LCG Length between Perpendiculars (LBP), L 115m Draught, T 3.92m Longitudinal position of the centre of buoayancy, LCB ‐5.0m Density, ρ 1.023 <strong>to</strong>nnes/m 3 Trim, t ‐0.5m < t < 1.0m Rudder Effectiveness, Reff 1.0 < (δCL/δα)r < 2.5 Skeg Effectiveness, Seff Non‐dimensionalised first derivative of sway force of the 0.0 < (δCL/δα)s < 3.0 bare hull with respect <strong>to</strong> sway Y � v0 velocity, Non‐dimensionalised first derivative of sway force of the ‐0.00495 bare hull with respect <strong>to</strong> turning 0.000973 rate, Y � r0 Non‐dimensionalised first derivative of yaw moment of the bare hull with respect <strong>to</strong> sway velocity, N � v0 Non‐dimensionalised first derivative of yaw moment of the bare hull with respect <strong>to</strong> N � r0 ‐0.00165 ‐0.000754 rate of turning, Rudder effectiveness fac<strong>to</strong>r, (δCL/δα)r 0.00045 Skeg effectiveness fac<strong>to</strong>r, (δCL/ δα)s (δCL/δα)r/2 Displacement 3708 <strong>to</strong>nnes | MARINE FRONTIER @ <strong>UniKL</strong> 9
Computation Output Extracts from the computation output based on the ship’s data input for t = ‐0.5, 1.0 < Reff < 2.5 and 0.0 < Seff < 3.0 are given below; Skeg Effectiveness Skeg Effectiveness 3.5 3 2.5 2 1.5 1 0.5 0 <strong>MIMET</strong> Technical Bulletin Volume 1 (2) 2010 Skeg Effectiveness 3.5 3 2.5 2 1.5 1 0.5 0 3.5 3 2.5 2 1.5 1 0.5 0 Directional Stability t=-0.5, Reff=1.5 -0.057 -0.034 -0.012 0.009 0.03 0.051 0.072 Directional Stability Criteria (a) Directional Stability t=-0.5, Reff=1.75 -0.032 -0.01 0.012 0.033 0.054 0.074 0.095 Directional Stability Criteria (b) Directional Stability t=-0.5, Reff=2 -0.008 0.014 0.035 0.056 0.077 0.097 0.118 Directional Stability Criteria (c) | MARINE FRONTIER @ <strong>UniKL</strong> 10