excel-codegen 0.9.0__py3-none-any.whl

excel_codegen/examples/abs_fpi/generated/EXT-T15-mu90-kf+1.js

@@ -0,0 +1,91 @@
+// EXT-T15-mu90-kf+1
+// The rules live in the RULE frame: x from the A.P., y from the centre line
+// positive starboard, z above the baseline.  The model usually is not, so the
+// three mapping cells (ax_* / neg_*) convert it.  All three coordinates go
+// through temporaries first: assigning them in place would read a coordinate
+// that has already been overwritten.
+var t1 = x;
+var t2 = y;
+var t3 = z;
+x = t1;
+y = t2;
+z = t3;
+// --- global parameters (same for every case)
+var L = 340 m;
+var B = 64 m;
+var D = 32 m;
+var r_b = 0.703 m;
+var freeboard = 0.5 m;
+var rho_sea = 1025 kg/m^3;
+var g = 9.81 m/s^2;
+var Cb = 1;
+// --- this case
+var draft = 15 m;
+var mu_deg = 90 deg;
+var k_c = 1;
+var k_f0 = 1;
+var k_u = 1.1;
+var beta_EPS = 0.733;
+var beta_EPP = 0.733;
+var k_esf = 1.1;
+var x_o = 51 m;
+// --- the calculation
+var u = mu_deg;   // the heading, as an Angle -- cos/sin take it directly
+var C_1 = 0;   // wave coefficient, 0 until a length bracket matches
+// C_1 brackets of 3-2-1/3.5, closed at both ends; outside them C_1 is zero
+if      (L >= 61 m  && L <= 90 m ) C_1 = 0.044 * (L / 1 m) + 3.75;
+else if (L >  90 m  && L <  300 m) C_1 = 10.75 - Math.pow((300 m - L) / 100 m, 1.5);
+else if (L >= 300 m && L <  350 m) C_1 = 10.75;
+else if (L >= 350 m && L <= 500 m) C_1 = 10.75 - Math.pow((L - 350 m) / 150 m, 1.5);
+else                               C_1 = 0;
+var h_do = 1.36 * C_1 * 1 m;   // reference hydrodynamic head, k = 1 (5.5.1)
+var z_bilge = r_b;   // rule z of the bilge (the bilge radius)
+var z_freeboard = D - freeboard;   // rule z of the highest deck at side
+var xL = x / L;   // x / L, dimensionless
+var k_lo = 1;   // k_lo = 1.0 amidships
+// Fig. 9: k_lo = 1 amidships, rising towards both ends
+if      (xL < 0.2) k_lo = 1.5 - xL * 2.5;
+else if (xL > 0.7) k_lo = 1 + (xL - 0.7) * 5;
+else               k_lo = 1;
+var k_l = 1 + (k_lo - 1) * Math.cos(u);   // k_l = 1 + (k_lo - 1) cos mu; GeniE converts the angle inside cos
+var k_f = k_f0 * (1 - (1 - Math.cos((x - x_o) / L * 360 deg)) * Math.cos(u));   // k_f, the reference station against the wave
+var alpha_1 = (1 - 0.25 * Math.cos(u)) * k_f;   // i = 1, at the waterline, starboard
+var alpha_2 = (0.40 - 0.10 * Math.cos(u)) * k_f;   // i = 2, at the bilge, starboard
+var alpha_3 = (0.30 - 0.20 * Math.sin(u)) * k_f;   // i = 3, at the bottom centre line
+var alpha_5 = (0.75 - 1.25 * Math.sin(u)) * k_f;   // i = 5, at the waterline, port
+var alpha_4 = 2 * alpha_3 - alpha_2;   // i = 4 = 2 alpha_3 - alpha_2, at the bilge, port
+var h_d1 = k_l * alpha_1 * h_do;   // hydrodynamic head at i = 1
+var h_d5 = k_l * alpha_5 * h_do;   // hydrodynamic head at i = 5
+var h_star_s = Math.min(z_freeboard - draft, k_u * h_d1);   // Fig. 10: the head above the waterline, capped at the freeboard, starboard
+var h_star_p = Math.min(z_freeboard - draft, k_u * h_d5);   // Fig. 10: the head above the waterline, capped at the freeboard, port
+var alpha_i = 0;   // alpha at the point evaluated, assigned by the girth branches
+// Fig. 8: alpha at the point evaluated
+if (z >= z_bilge) {
+    if (z <= draft) {
+        // bilge to waterline, on the side of the girth the point is on
+        if (y >= 0 m) alpha_i = alpha_2 + (alpha_1 - alpha_2) * (z - z_bilge) / (draft - z_bilge);
+        else          alpha_i = alpha_4 + (alpha_5 - alpha_4) * (z - z_bilge) / (draft - z_bilge);
+    } else {
+        // above the waterline, faded out over that side's h_star
+        if (y >= 0 m) {
+            if (h_star_s <= 0 m) alpha_i = 0;
+            else                 alpha_i = alpha_1 - (z - draft) / h_star_s * alpha_1;
+        } else {
+            if (h_star_p <= 0 m) alpha_i = 0;
+            else                 alpha_i = alpha_5 - (z - draft) / h_star_p * alpha_5;
+        }
+    }
+} else {
+    // below the bilge: linear across the flat of bottom
+    alpha_i = alpha_4 + (alpha_2 - alpha_4) * (z + z_bilge) / (2 * z_bilge);
+}
+// 5.5.1: which side of the girth, and so which ESF
+// the ESF is a property of the SIDE the point is on: EPS to starboard, EPP to port
+var ESF_side = 0;
+if (y >= 0 m) ESF_side = beta_EPS * k_esf;
+else          ESF_side = beta_EPP * k_esf;
+var h_di = k_l * alpha_i * h_do;   // hydrodynamic head at the point evaluated (5.5.1)
+var h_de = k_c * h_di;   // h_de = k_c h_di (5.5.1 nominal form)
+var p_s = Math.max(rho_sea * g * (draft - z), 0 Pa);   // still-water pressure, zero above the waterline
+var p_d = rho_sea * g * ESF_side * k_u * h_de;   // hydrodynamic pressure, with the ESF of the side the point is on
+return Math.max(p_s + p_d, 0 Pa);   // the pressure the function returns

excel_codegen/examples/abs_fpi/generated/EXT-T15-mu90-kf+1_summary.md

@@ -0,0 +1,23 @@
+| 项目 | 值 | 出处 |
+| --- | --- | --- |
+| case_name | EXT-T15-mu90-kf+1 | 工况名（Excel 表头） |
+| kind | EXT | 规则集（多规则集共用一本工作簿时用 case_filter 分流） |
+| draft | 15 m | 装载工况 |
+| mu_deg | 90 deg | Table 1A D |
+| k_c | 1 | Table 1A B |
+| k_f0 | 1 | Table 1A B |
+| k_u | 1.1 | 5A-3-2/7, /9 |
+| beta_EPS | 0.733 | 5A-3-A1/3 |
+| beta_EPP | 0.733 | 5A-3-A1/3 |
+| k_esf | 1.1 | 项目 |
+| x_o | 51 m | Fig. 9 脚注 |
+| L | 340 m | 3-1-1/3 |
+| B | 64 m | 3-1-1/3 |
+| D | 32 m | 3-1-1/3 |
+| r_b | 0.703 m | Fig. 11 |
+| freeboard | 0.5 m | Fig. 10 |
+| rho_sea | 1025 kg/m^3 | 5.5.1 |
+| g | 9.81 m/s^2 | 5.5.1 |
+| Cb | 1 | 3-1-1/13.3 |
+| axis ax | (1, 2, 3) | 模型 |
+| axis neg | ( /  / ) | 模型 |

excel_codegen/examples/abs_fpi/generated/EXT-T20.559-mu90-kf+1.js

@@ -0,0 +1,91 @@
+// EXT-T20.559-mu90-kf+1
+// The rules live in the RULE frame: x from the A.P., y from the centre line
+// positive starboard, z above the baseline.  The model usually is not, so the
+// three mapping cells (ax_* / neg_*) convert it.  All three coordinates go
+// through temporaries first: assigning them in place would read a coordinate
+// that has already been overwritten.
+var t1 = x;
+var t2 = y;
+var t3 = z;
+x = t1;
+y = t2;
+z = t3;
+// --- global parameters (same for every case)
+var L = 340 m;
+var B = 64 m;
+var D = 32 m;
+var r_b = 0.703 m;
+var freeboard = 0.5 m;
+var rho_sea = 1025 kg/m^3;
+var g = 9.81 m/s^2;
+var Cb = 1;
+// --- this case
+var draft = 20.559 m;
+var mu_deg = 90 deg;
+var k_c = 1;
+var k_f0 = 1;
+var k_u = 1.1;
+var beta_EPS = 0.733;
+var beta_EPP = 0.733;
+var k_esf = 1.1;
+var x_o = 51 m;
+// --- the calculation
+var u = mu_deg;   // the heading, as an Angle -- cos/sin take it directly
+var C_1 = 0;   // wave coefficient, 0 until a length bracket matches
+// C_1 brackets of 3-2-1/3.5, closed at both ends; outside them C_1 is zero
+if      (L >= 61 m  && L <= 90 m ) C_1 = 0.044 * (L / 1 m) + 3.75;
+else if (L >  90 m  && L <  300 m) C_1 = 10.75 - Math.pow((300 m - L) / 100 m, 1.5);
+else if (L >= 300 m && L <  350 m) C_1 = 10.75;
+else if (L >= 350 m && L <= 500 m) C_1 = 10.75 - Math.pow((L - 350 m) / 150 m, 1.5);
+else                               C_1 = 0;
+var h_do = 1.36 * C_1 * 1 m;   // reference hydrodynamic head, k = 1 (5.5.1)
+var z_bilge = r_b;   // rule z of the bilge (the bilge radius)
+var z_freeboard = D - freeboard;   // rule z of the highest deck at side
+var xL = x / L;   // x / L, dimensionless
+var k_lo = 1;   // k_lo = 1.0 amidships
+// Fig. 9: k_lo = 1 amidships, rising towards both ends
+if      (xL < 0.2) k_lo = 1.5 - xL * 2.5;
+else if (xL > 0.7) k_lo = 1 + (xL - 0.7) * 5;
+else               k_lo = 1;
+var k_l = 1 + (k_lo - 1) * Math.cos(u);   // k_l = 1 + (k_lo - 1) cos mu; GeniE converts the angle inside cos
+var k_f = k_f0 * (1 - (1 - Math.cos((x - x_o) / L * 360 deg)) * Math.cos(u));   // k_f, the reference station against the wave
+var alpha_1 = (1 - 0.25 * Math.cos(u)) * k_f;   // i = 1, at the waterline, starboard
+var alpha_2 = (0.40 - 0.10 * Math.cos(u)) * k_f;   // i = 2, at the bilge, starboard
+var alpha_3 = (0.30 - 0.20 * Math.sin(u)) * k_f;   // i = 3, at the bottom centre line
+var alpha_5 = (0.75 - 1.25 * Math.sin(u)) * k_f;   // i = 5, at the waterline, port
+var alpha_4 = 2 * alpha_3 - alpha_2;   // i = 4 = 2 alpha_3 - alpha_2, at the bilge, port
+var h_d1 = k_l * alpha_1 * h_do;   // hydrodynamic head at i = 1
+var h_d5 = k_l * alpha_5 * h_do;   // hydrodynamic head at i = 5
+var h_star_s = Math.min(z_freeboard - draft, k_u * h_d1);   // Fig. 10: the head above the waterline, capped at the freeboard, starboard
+var h_star_p = Math.min(z_freeboard - draft, k_u * h_d5);   // Fig. 10: the head above the waterline, capped at the freeboard, port
+var alpha_i = 0;   // alpha at the point evaluated, assigned by the girth branches
+// Fig. 8: alpha at the point evaluated
+if (z >= z_bilge) {
+    if (z <= draft) {
+        // bilge to waterline, on the side of the girth the point is on
+        if (y >= 0 m) alpha_i = alpha_2 + (alpha_1 - alpha_2) * (z - z_bilge) / (draft - z_bilge);
+        else          alpha_i = alpha_4 + (alpha_5 - alpha_4) * (z - z_bilge) / (draft - z_bilge);
+    } else {
+        // above the waterline, faded out over that side's h_star
+        if (y >= 0 m) {
+            if (h_star_s <= 0 m) alpha_i = 0;
+            else                 alpha_i = alpha_1 - (z - draft) / h_star_s * alpha_1;
+        } else {
+            if (h_star_p <= 0 m) alpha_i = 0;
+            else                 alpha_i = alpha_5 - (z - draft) / h_star_p * alpha_5;
+        }
+    }
+} else {
+    // below the bilge: linear across the flat of bottom
+    alpha_i = alpha_4 + (alpha_2 - alpha_4) * (z + z_bilge) / (2 * z_bilge);
+}
+// 5.5.1: which side of the girth, and so which ESF
+// the ESF is a property of the SIDE the point is on: EPS to starboard, EPP to port
+var ESF_side = 0;
+if (y >= 0 m) ESF_side = beta_EPS * k_esf;
+else          ESF_side = beta_EPP * k_esf;
+var h_di = k_l * alpha_i * h_do;   // hydrodynamic head at the point evaluated (5.5.1)
+var h_de = k_c * h_di;   // h_de = k_c h_di (5.5.1 nominal form)
+var p_s = Math.max(rho_sea * g * (draft - z), 0 Pa);   // still-water pressure, zero above the waterline
+var p_d = rho_sea * g * ESF_side * k_u * h_de;   // hydrodynamic pressure, with the ESF of the side the point is on
+return Math.max(p_s + p_d, 0 Pa);   // the pressure the function returns

excel_codegen/examples/abs_fpi/generated/EXT-T20.559-mu90-kf+1_summary.md

@@ -0,0 +1,23 @@
+| 项目 | 值 | 出处 |
+| --- | --- | --- |
+| case_name | EXT-T20.559-mu90-kf+1 | 工况名（Excel 表头） |
+| kind | EXT | 规则集（多规则集共用一本工作簿时用 case_filter 分流） |
+| draft | 20.559 m | 装载工况 |
+| mu_deg | 90 deg | Table 1A D |
+| k_c | 1 | Table 1A B |
+| k_f0 | 1 | Table 1A B |
+| k_u | 1.1 | 5A-3-2/7, /9 |
+| beta_EPS | 0.733 | 5A-3-A1/3 |
+| beta_EPP | 0.733 | 5A-3-A1/3 |
+| k_esf | 1.1 | 项目 |
+| x_o | 51 m | Fig. 9 脚注 |
+| L | 340 m | 3-1-1/3 |
+| B | 64 m | 3-1-1/3 |
+| D | 32 m | 3-1-1/3 |
+| r_b | 0.703 m | Fig. 11 |
+| freeboard | 0.5 m | Fig. 10 |
+| rho_sea | 1025 kg/m^3 | 5.5.1 |
+| g | 9.81 m/s^2 | 5.5.1 |
+| Cb | 1 | 3-1-1/13.3 |
+| axis ax | (1, 2, 3) | 模型 |
+| axis neg | ( /  / ) | 模型 |

excel_codegen/examples/abs_fpi/generated/WBT6-d15.059-mu0.js

@@ -0,0 +1,124 @@
+// WBT6-d15.059-mu0
+// See ext_body.js.j2 for why the coordinates go through temporaries.
+var t1 = x;
+var t2 = y;
+var t3 = z;
+x = t1;
+y = t2;
+z = t3;
+// --- global parameters (same for every case)
+var L = 340 m;
+var B = 64 m;
+var D = 32 m;
+var rho_sea = 1025 kg/m^3;
+var g = 9.81 m/s^2;
+var Cb = 1;
+var d_f = 20.559 m;
+var rho_g_lo = 10050 N/m^3;
+var rho_g_hi = 11180 N/m^3;
+// --- this case
+var draft = 15.059 m;
+var mu_deg = 0 deg;
+var k_c = 0.4;
+var k_u = 1.1;
+var beta_VAC = 0.66;
+var beta_LAC = 0.758;
+var beta_TAC = 0.557;
+var beta_PMO = 0.686;
+var beta_RMO = 0.492;
+var k_esf = 1.1;
+var w_v = 0.75;
+var w_l = 0.25;
+var w_t = 0.75;
+var xi = 21 m;
+var delta_b = 0 m;
+var delta_h = 0 m;
+// --- the tank this case names: WBT6
+var rho_tank = 1025 kg/m^3;
+var l_tank = 42 m;
+var b_tank = 32 m;
+var h_tank = 32 m;
+var eta_deck = 0 m;
+var eta_overflow = 0 m;
+var C_dp = 1;
+var C_ru = 1;
+var p_vp = 0 N/cm^2;
+var GM_full_in = 0 m;
+var k_r_in = 0 m;
+var tank_is_ballast = 1;
+var member_11_17 = 0;
+// --- the calculation
+var C_theta = 1;   // C_theta, the roll angle factor of 5.7.2(b)
+var C_phi = -0.45;   // C_phi, the pitch angle factor of 5.7.2(b)
+var k_theta = 0.005;   // k_theta of the theta expression (5.7.1(b))
+var C_R = 1.05;   // C_R of the theta expression (5.7.1(b))
+var k_o = 1.34 - 0.47 * Cb;   // k_o (5.7.1(c))
+var a_o = (k_o * (2.4 / Math.pow(L / 1 m, 0.5) + 34 / (L / 1 m) - 600 / Math.pow(L / 1 m, 2))) * 1 m;   // a_o, the vertical acceleration parameter
+var D_1 = (10 / Cb);   // the (10 / C_b) term of phi
+var phi = beta_PMO * k_esf * 1030 * Math.pow(D_1, 0.25) / (L / 1 m) * 1 deg;   // phi, the pitch amplitude
+if (phi > 10 deg) phi = 10 deg;
+else              phi = phi;
+var GM_full = GM_full_in > 0 m ? GM_full_in : 0.12 * B;   // 5.7.1(b): 0.12 B when GM(full) is not available
+var k_r = k_r_in > 0 m ? k_r_in : (draft <= d_f ? 0.35 * B : 0.45 * B);   // roll radius of gyration
+var Delta = 10.05 kN/m^3 * L * B * d_f * Cb;   // displacement, k_d = 10.05 kN/m^3
+var C_di = 1.06 * (draft / d_f) - 0.06;   // C_di = 1.06 (d_i/d_f) - 0.06
+var GM = GM_full;   // GM at the draft evaluated, filled in below
+if      (draft / d_f >= 1)   GM = GM_full;
+else if (draft / d_f >= 0.9) GM = 1.1 * GM_full;
+else if (draft / d_f >= 2/3) GM = 1.5 * GM_full;
+else if (draft / d_f >= 0.5) GM = 2.0 * GM_full;
+else                         GM = GM_full;
+var T_r = 2 * k_r * Math.pow(GM / 1 m, -0.5) / 1 m * 1 s;   // T_r, the roll natural period, in seconds
+var theta = C_R * beta_RMO * k_esf * (35 - k_theta * C_di * Delta / 1000 kN) * 1 deg;   // theta, the roll amplitude
+if      (T_r > 20 s)              theta = theta;
+else if (T_r >= 12.5 s)           theta = theta * (1.5375 - 0.027 * T_r / 1 s);
+else                              theta = theta * (0.8625 + 0.027 * T_r / 1 s);
+if (theta > 30 deg) theta = 30 deg;
+else                theta = theta;
+var k_v = Math.pow(1 + 0.65 * Math.pow(5.3 - 45 / (L / 1 m), 2) * Math.pow(x / L - 0.45, 2), 0.5);   // k_v (5.7.1(c))
+var C_v = Math.cos(mu_deg) + (1 + 2.4 * z / B) * Math.sin(mu_deg) / k_v;   // C_v (5.7.1(c))
+var C_l = 0.35 - 0.0005 * (L / 1 m - 200);   // C_l (5.7.1(c))
+var k_l = 0.5 + 8 * y / L;   // k_l (5.7.1(c)), the VERTICAL distribution
+var C_t = 1.27 * Math.pow(1 + 1.52 * Math.pow(x / L - 0.45, 2), 0.5);   // C_t (5.7.1(c))
+var k_t = 0.35 + y / B;   // k_t (5.7.1(c)), the VERTICAL distribution
+var a_v = C_v * beta_VAC * k_esf * k_v * a_o / 1 m * g;   // a_v, positive downward
+var a_l = C_l * beta_LAC * k_esf * k_l * a_o / 1 m * g;   // a_l, positive forward
+var a_t = C_t * beta_TAC * k_esf * k_t * a_o / 1 m * g;   // a_t, positive starboard
+var rho_g = rho_tank * g;   // the liquid specific weight (5.7.2(a))
+var eta_overflow_add = Math.max(0.6666667 * eta_overflow, 760 mm);   // Fig. 12: 2/3 of the overflow distance, at least 760 mm
+var eta = h_tank + eta_deck + (eta_overflow > 0 m ? eta_overflow_add : 0 m) - z;   // eta, measured down from the tank top
+var k_s = 1;   // k_s, filled in below
+if (tank_is_ballast != 0 || member_11_17 == 0) {
+    k_s = 1;
+} else {
+    k_s = 0.878 + (rho_g - rho_g_lo) / (rho_g_hi - rho_g_lo) * 0.122;
+    if (k_s < 0.878) k_s = 0.878;
+    else             k_s = k_s;
+    if (k_s > 1)     k_s = 1;
+    else             k_s = k_s;
+}
+var p_o = p_vp > 0 N/cm^2 ? Math.max(p_vp - 2.06 N/cm^2, 0 N/cm^2) : 0 N/cm^2;   // p_o = (p_vp - p_n), p_n = 2.06 N/cm^2
+var a_i = 0.71 * C_dp * (w_v * a_v + w_l * (l_tank / h_tank) * a_l + w_t * (b_tank / h_tank) * a_t);   // a_i (5.7.2(a))
+var theta_e = 0.71 * C_theta * theta;   // theta_e (5.7.2(b))
+var phi_e = 0.71 * C_phi * phi;   // phi_e (5.7.2(b))
+var zeta = Math.abs(b_tank / 2 - z);   // zeta, from the tank centre line
+var zeta_e = zeta;   // zeta_e, filled in below
+var eta_e = eta;   // eta_e, filled in below
+var dh_i = 0 m;   // delta h_i, filled in below
+if (phi_e < 0 deg && theta_e > 0 deg) {
+    // branch i of 5.7.2(b)
+    zeta_e = b_tank - zeta;
+    eta_e = eta;
+    dh_i = xi * Math.sin(-phi_e) + C_ru * (zeta_e * Math.sin(theta_e) * Math.cos(phi_e)
+         + eta_e * Math.cos(theta_e) * Math.cos(phi_e) - eta);
+} else {
+    // branch ii of 5.7.2(b)
+    zeta_e = zeta - delta_b;
+    eta_e = eta - delta_h;
+    dh_i = (l_tank - xi) * Math.sin(-phi_e) + C_ru * (zeta_e * Math.sin(theta_e) * Math.cos(phi_e)
+         + eta_e * Math.cos(theta_e) * Math.cos(phi_e) - eta);
+}
+var h_d = k_c * (eta * a_i / g + dh_i);   // h_d = k_c (eta a_i / g + delta h_i) (5.7.2(a))
+var p_s = k_s * rho_g * eta;   // static pressure (5.7.2(a))
+var p_d = k_s * rho_g * k_u * h_d;   // dynamic pressure (5.7.2(a))
+return p_s + p_d + p_o;   // the pressure the function returns

excel_codegen/examples/abs_fpi/generated/WBT6-d15.059-mu0_summary.md

@@ -0,0 +1,30 @@
+| 项目 | 值 | 出处 |
+| --- | --- | --- |
+| case_name | WBT6-d15.059-mu0 | 工况名（Excel 表头） |
+| kind | INT | 规则集（多规则集共用一本工作簿时用 case_filter 分流） |
+| tank_ref | WBT6 | 指向成员表 `Tank Data` 的哪一行（指南 §18） |
+| draft | 15.059 m | 装载工况 |
+| mu_deg | 0 deg | Table 1A D |
+| k_c | 0.4 | Table 1A C |
+| k_u | 1.1 | 5A-3-2/7, /9 |
+| beta_VAC | 0.66 | 5A-3-A1/3 |
+| beta_LAC | 0.758 | 5A-3-A1/3 |
+| beta_TAC | 0.557 | 5A-3-A1/3 |
+| beta_PMO | 0.686 | 5A-3-A1/3 |
+| beta_RMO | 0.492 | 5A-3-A1/3 |
+| k_esf | 1.1 | 项目 |
+| w_v / w_l / w_t | 0.75 / 0.25 / 0.75 | Table 1A B |
+| xi | 21 m | 5.7.2(b) |
+| delta_b / delta_h | 0 m / 0 m | 5.7.2(b) |
+| rho_tank | 1025 kg/m^3 | 成员表 Tank Data |
+| l_tank / b_tank / h_tank | 42 m / 32 m / 32 m | Fig. 11 |
+| eta_deck / eta_overflow | 0 m / 0 m | Fig. 12 |
+| C_dp / C_ru | 1 / 1 | 5.7.2(d) |
+| p_vp | 0 N/cm^2 | 成员表 Tank Data |
+| GM_full_in / k_r_in | 0 m / 0 m | 稳性计算书 |
+| tank_is_ballast / member_11_17 | 1 / 0 | Table 3 |
+| L / B / D / d_f | 340 m / 64 m / 32 m / 20.559 m | 3-1-1/3 |
+| rho_sea / g / Cb | 1025 kg/m^3 / 9.81 m/s^2 / 1 | 5.7.1 |
+| rho_g_lo / rho_g_hi | 10050 N/m^3 / 11180 N/m^3 | 5.7.2(a) |
+| axis ax | (1, 2, 3) | 模型 |
+| axis neg | ( /  / ) | 模型 |

excel_codegen/examples/abs_fpi/generated/WBT6-d8-mu90.js

@@ -0,0 +1,124 @@
+// WBT6-d8-mu90
+// See ext_body.js.j2 for why the coordinates go through temporaries.
+var t1 = x;
+var t2 = y;
+var t3 = z;
+x = t1;
+y = t2;
+z = t3;
+// --- global parameters (same for every case)
+var L = 340 m;
+var B = 64 m;
+var D = 32 m;
+var rho_sea = 1025 kg/m^3;
+var g = 9.81 m/s^2;
+var Cb = 1;
+var d_f = 20.559 m;
+var rho_g_lo = 10050 N/m^3;
+var rho_g_hi = 11180 N/m^3;
+// --- this case
+var draft = 8 m;
+var mu_deg = 90 deg;
+var k_c = 1;
+var k_u = 1.1;
+var beta_VAC = 0.66;
+var beta_LAC = 0.758;
+var beta_TAC = 0.557;
+var beta_PMO = 0.686;
+var beta_RMO = 0.492;
+var k_esf = 1.1;
+var w_v = 0.75;
+var w_l = 0.25;
+var w_t = 0.75;
+var xi = 21 m;
+var delta_b = 0 m;
+var delta_h = 0 m;
+// --- the tank this case names: WBT6
+var rho_tank = 1025 kg/m^3;
+var l_tank = 42 m;
+var b_tank = 32 m;
+var h_tank = 32 m;
+var eta_deck = 0 m;
+var eta_overflow = 0 m;
+var C_dp = 1;
+var C_ru = 1;
+var p_vp = 0 N/cm^2;
+var GM_full_in = 0 m;
+var k_r_in = 0 m;
+var tank_is_ballast = 1;
+var member_11_17 = 0;
+// --- the calculation
+var C_theta = 1;   // C_theta, the roll angle factor of 5.7.2(b)
+var C_phi = -0.45;   // C_phi, the pitch angle factor of 5.7.2(b)
+var k_theta = 0.005;   // k_theta of the theta expression (5.7.1(b))
+var C_R = 1.05;   // C_R of the theta expression (5.7.1(b))
+var k_o = 1.34 - 0.47 * Cb;   // k_o (5.7.1(c))
+var a_o = (k_o * (2.4 / Math.pow(L / 1 m, 0.5) + 34 / (L / 1 m) - 600 / Math.pow(L / 1 m, 2))) * 1 m;   // a_o, the vertical acceleration parameter
+var D_1 = (10 / Cb);   // the (10 / C_b) term of phi
+var phi = beta_PMO * k_esf * 1030 * Math.pow(D_1, 0.25) / (L / 1 m) * 1 deg;   // phi, the pitch amplitude
+if (phi > 10 deg) phi = 10 deg;
+else              phi = phi;
+var GM_full = GM_full_in > 0 m ? GM_full_in : 0.12 * B;   // 5.7.1(b): 0.12 B when GM(full) is not available
+var k_r = k_r_in > 0 m ? k_r_in : (draft <= d_f ? 0.35 * B : 0.45 * B);   // roll radius of gyration
+var Delta = 10.05 kN/m^3 * L * B * d_f * Cb;   // displacement, k_d = 10.05 kN/m^3
+var C_di = 1.06 * (draft / d_f) - 0.06;   // C_di = 1.06 (d_i/d_f) - 0.06
+var GM = GM_full;   // GM at the draft evaluated, filled in below
+if      (draft / d_f >= 1)   GM = GM_full;
+else if (draft / d_f >= 0.9) GM = 1.1 * GM_full;
+else if (draft / d_f >= 2/3) GM = 1.5 * GM_full;
+else if (draft / d_f >= 0.5) GM = 2.0 * GM_full;
+else                         GM = GM_full;
+var T_r = 2 * k_r * Math.pow(GM / 1 m, -0.5) / 1 m * 1 s;   // T_r, the roll natural period, in seconds
+var theta = C_R * beta_RMO * k_esf * (35 - k_theta * C_di * Delta / 1000 kN) * 1 deg;   // theta, the roll amplitude
+if      (T_r > 20 s)              theta = theta;
+else if (T_r >= 12.5 s)           theta = theta * (1.5375 - 0.027 * T_r / 1 s);
+else                              theta = theta * (0.8625 + 0.027 * T_r / 1 s);
+if (theta > 30 deg) theta = 30 deg;
+else                theta = theta;
+var k_v = Math.pow(1 + 0.65 * Math.pow(5.3 - 45 / (L / 1 m), 2) * Math.pow(x / L - 0.45, 2), 0.5);   // k_v (5.7.1(c))
+var C_v = Math.cos(mu_deg) + (1 + 2.4 * z / B) * Math.sin(mu_deg) / k_v;   // C_v (5.7.1(c))
+var C_l = 0.35 - 0.0005 * (L / 1 m - 200);   // C_l (5.7.1(c))
+var k_l = 0.5 + 8 * y / L;   // k_l (5.7.1(c)), the VERTICAL distribution
+var C_t = 1.27 * Math.pow(1 + 1.52 * Math.pow(x / L - 0.45, 2), 0.5);   // C_t (5.7.1(c))
+var k_t = 0.35 + y / B;   // k_t (5.7.1(c)), the VERTICAL distribution
+var a_v = C_v * beta_VAC * k_esf * k_v * a_o / 1 m * g;   // a_v, positive downward
+var a_l = C_l * beta_LAC * k_esf * k_l * a_o / 1 m * g;   // a_l, positive forward
+var a_t = C_t * beta_TAC * k_esf * k_t * a_o / 1 m * g;   // a_t, positive starboard
+var rho_g = rho_tank * g;   // the liquid specific weight (5.7.2(a))
+var eta_overflow_add = Math.max(0.6666667 * eta_overflow, 760 mm);   // Fig. 12: 2/3 of the overflow distance, at least 760 mm
+var eta = h_tank + eta_deck + (eta_overflow > 0 m ? eta_overflow_add : 0 m) - z;   // eta, measured down from the tank top
+var k_s = 1;   // k_s, filled in below
+if (tank_is_ballast != 0 || member_11_17 == 0) {
+    k_s = 1;
+} else {
+    k_s = 0.878 + (rho_g - rho_g_lo) / (rho_g_hi - rho_g_lo) * 0.122;
+    if (k_s < 0.878) k_s = 0.878;
+    else             k_s = k_s;
+    if (k_s > 1)     k_s = 1;
+    else             k_s = k_s;
+}
+var p_o = p_vp > 0 N/cm^2 ? Math.max(p_vp - 2.06 N/cm^2, 0 N/cm^2) : 0 N/cm^2;   // p_o = (p_vp - p_n), p_n = 2.06 N/cm^2
+var a_i = 0.71 * C_dp * (w_v * a_v + w_l * (l_tank / h_tank) * a_l + w_t * (b_tank / h_tank) * a_t);   // a_i (5.7.2(a))
+var theta_e = 0.71 * C_theta * theta;   // theta_e (5.7.2(b))
+var phi_e = 0.71 * C_phi * phi;   // phi_e (5.7.2(b))
+var zeta = Math.abs(b_tank / 2 - z);   // zeta, from the tank centre line
+var zeta_e = zeta;   // zeta_e, filled in below
+var eta_e = eta;   // eta_e, filled in below
+var dh_i = 0 m;   // delta h_i, filled in below
+if (phi_e < 0 deg && theta_e > 0 deg) {
+    // branch i of 5.7.2(b)
+    zeta_e = b_tank - zeta;
+    eta_e = eta;
+    dh_i = xi * Math.sin(-phi_e) + C_ru * (zeta_e * Math.sin(theta_e) * Math.cos(phi_e)
+         + eta_e * Math.cos(theta_e) * Math.cos(phi_e) - eta);
+} else {
+    // branch ii of 5.7.2(b)
+    zeta_e = zeta - delta_b;
+    eta_e = eta - delta_h;
+    dh_i = (l_tank - xi) * Math.sin(-phi_e) + C_ru * (zeta_e * Math.sin(theta_e) * Math.cos(phi_e)
+         + eta_e * Math.cos(theta_e) * Math.cos(phi_e) - eta);
+}
+var h_d = k_c * (eta * a_i / g + dh_i);   // h_d = k_c (eta a_i / g + delta h_i) (5.7.2(a))
+var p_s = k_s * rho_g * eta;   // static pressure (5.7.2(a))
+var p_d = k_s * rho_g * k_u * h_d;   // dynamic pressure (5.7.2(a))
+return p_s + p_d + p_o;   // the pressure the function returns

excel_codegen/examples/abs_fpi/generated/WBT6-d8-mu90_summary.md

@@ -0,0 +1,30 @@
+| 项目 | 值 | 出处 |
+| --- | --- | --- |
+| case_name | WBT6-d8-mu90 | 工况名（Excel 表头） |
+| kind | INT | 规则集（多规则集共用一本工作簿时用 case_filter 分流） |
+| tank_ref | WBT6 | 指向成员表 `Tank Data` 的哪一行（指南 §18） |
+| draft | 8 m | 装载工况 |
+| mu_deg | 90 deg | Table 1A D |
+| k_c | 1 | Table 1A C |
+| k_u | 1.1 | 5A-3-2/7, /9 |
+| beta_VAC | 0.66 | 5A-3-A1/3 |
+| beta_LAC | 0.758 | 5A-3-A1/3 |
+| beta_TAC | 0.557 | 5A-3-A1/3 |
+| beta_PMO | 0.686 | 5A-3-A1/3 |
+| beta_RMO | 0.492 | 5A-3-A1/3 |
+| k_esf | 1.1 | 项目 |
+| w_v / w_l / w_t | 0.75 / 0.25 / 0.75 | Table 1A B |
+| xi | 21 m | 5.7.2(b) |
+| delta_b / delta_h | 0 m / 0 m | 5.7.2(b) |
+| rho_tank | 1025 kg/m^3 | 成员表 Tank Data |
+| l_tank / b_tank / h_tank | 42 m / 32 m / 32 m | Fig. 11 |
+| eta_deck / eta_overflow | 0 m / 0 m | Fig. 12 |
+| C_dp / C_ru | 1 / 1 | 5.7.2(d) |
+| p_vp | 0 N/cm^2 | 成员表 Tank Data |
+| GM_full_in / k_r_in | 0 m / 0 m | 稳性计算书 |
+| tank_is_ballast / member_11_17 | 1 / 0 | Table 3 |
+| L / B / D / d_f | 340 m / 64 m / 32 m / 20.559 m | 3-1-1/3 |
+| rho_sea / g / Cb | 1025 kg/m^3 / 9.81 m/s^2 / 1 | 5.7.1 |
+| rho_g_lo / rho_g_hi | 10050 N/m^3 / 11180 N/m^3 | 5.7.2(a) |
+| axis ax | (1, 2, 3) | 模型 |
+| axis neg | ( /  / ) | 模型 |

excel_codegen/examples/abs_fpi/generated/WBT7-d20.559-mu90.js

@@ -0,0 +1,124 @@
+// WBT7-d20.559-mu90
+// See ext_body.js.j2 for why the coordinates go through temporaries.
+var t1 = x;
+var t2 = y;
+var t3 = z;
+x = t1;
+y = t2;
+z = t3;
+// --- global parameters (same for every case)
+var L = 340 m;
+var B = 64 m;
+var D = 32 m;
+var rho_sea = 1025 kg/m^3;
+var g = 9.81 m/s^2;
+var Cb = 1;
+var d_f = 20.559 m;
+var rho_g_lo = 10050 N/m^3;
+var rho_g_hi = 11180 N/m^3;
+// --- this case
+var draft = 20.559 m;
+var mu_deg = 90 deg;
+var k_c = 1;
+var k_u = 1.1;
+var beta_VAC = 0.66;
+var beta_LAC = 0.758;
+var beta_TAC = 0.557;
+var beta_PMO = 0.686;
+var beta_RMO = 0.492;
+var k_esf = 1.1;
+var w_v = 0.75;
+var w_l = 0.25;
+var w_t = 0.75;
+var xi = 7.25 m;
+var delta_b = 0 m;
+var delta_h = 0 m;
+// --- the tank this case names: WBT7
+var rho_tank = 1025 kg/m^3;
+var l_tank = 14.5 m;
+var b_tank = 10.15 m;
+var h_tank = 32 m;
+var eta_deck = 0 m;
+var eta_overflow = 0 m;
+var C_dp = 1;
+var C_ru = 1;
+var p_vp = 0 N/cm^2;
+var GM_full_in = 0 m;
+var k_r_in = 0 m;
+var tank_is_ballast = 1;
+var member_11_17 = 0;
+// --- the calculation
+var C_theta = 1;   // C_theta, the roll angle factor of 5.7.2(b)
+var C_phi = -0.45;   // C_phi, the pitch angle factor of 5.7.2(b)
+var k_theta = 0.005;   // k_theta of the theta expression (5.7.1(b))
+var C_R = 1.05;   // C_R of the theta expression (5.7.1(b))
+var k_o = 1.34 - 0.47 * Cb;   // k_o (5.7.1(c))
+var a_o = (k_o * (2.4 / Math.pow(L / 1 m, 0.5) + 34 / (L / 1 m) - 600 / Math.pow(L / 1 m, 2))) * 1 m;   // a_o, the vertical acceleration parameter
+var D_1 = (10 / Cb);   // the (10 / C_b) term of phi
+var phi = beta_PMO * k_esf * 1030 * Math.pow(D_1, 0.25) / (L / 1 m) * 1 deg;   // phi, the pitch amplitude
+if (phi > 10 deg) phi = 10 deg;
+else              phi = phi;
+var GM_full = GM_full_in > 0 m ? GM_full_in : 0.12 * B;   // 5.7.1(b): 0.12 B when GM(full) is not available
+var k_r = k_r_in > 0 m ? k_r_in : (draft <= d_f ? 0.35 * B : 0.45 * B);   // roll radius of gyration
+var Delta = 10.05 kN/m^3 * L * B * d_f * Cb;   // displacement, k_d = 10.05 kN/m^3
+var C_di = 1.06 * (draft / d_f) - 0.06;   // C_di = 1.06 (d_i/d_f) - 0.06
+var GM = GM_full;   // GM at the draft evaluated, filled in below
+if      (draft / d_f >= 1)   GM = GM_full;
+else if (draft / d_f >= 0.9) GM = 1.1 * GM_full;
+else if (draft / d_f >= 2/3) GM = 1.5 * GM_full;
+else if (draft / d_f >= 0.5) GM = 2.0 * GM_full;
+else                         GM = GM_full;
+var T_r = 2 * k_r * Math.pow(GM / 1 m, -0.5) / 1 m * 1 s;   // T_r, the roll natural period, in seconds
+var theta = C_R * beta_RMO * k_esf * (35 - k_theta * C_di * Delta / 1000 kN) * 1 deg;   // theta, the roll amplitude
+if      (T_r > 20 s)              theta = theta;
+else if (T_r >= 12.5 s)           theta = theta * (1.5375 - 0.027 * T_r / 1 s);
+else                              theta = theta * (0.8625 + 0.027 * T_r / 1 s);
+if (theta > 30 deg) theta = 30 deg;
+else                theta = theta;
+var k_v = Math.pow(1 + 0.65 * Math.pow(5.3 - 45 / (L / 1 m), 2) * Math.pow(x / L - 0.45, 2), 0.5);   // k_v (5.7.1(c))
+var C_v = Math.cos(mu_deg) + (1 + 2.4 * z / B) * Math.sin(mu_deg) / k_v;   // C_v (5.7.1(c))
+var C_l = 0.35 - 0.0005 * (L / 1 m - 200);   // C_l (5.7.1(c))
+var k_l = 0.5 + 8 * y / L;   // k_l (5.7.1(c)), the VERTICAL distribution
+var C_t = 1.27 * Math.pow(1 + 1.52 * Math.pow(x / L - 0.45, 2), 0.5);   // C_t (5.7.1(c))
+var k_t = 0.35 + y / B;   // k_t (5.7.1(c)), the VERTICAL distribution
+var a_v = C_v * beta_VAC * k_esf * k_v * a_o / 1 m * g;   // a_v, positive downward
+var a_l = C_l * beta_LAC * k_esf * k_l * a_o / 1 m * g;   // a_l, positive forward
+var a_t = C_t * beta_TAC * k_esf * k_t * a_o / 1 m * g;   // a_t, positive starboard
+var rho_g = rho_tank * g;   // the liquid specific weight (5.7.2(a))
+var eta_overflow_add = Math.max(0.6666667 * eta_overflow, 760 mm);   // Fig. 12: 2/3 of the overflow distance, at least 760 mm
+var eta = h_tank + eta_deck + (eta_overflow > 0 m ? eta_overflow_add : 0 m) - z;   // eta, measured down from the tank top
+var k_s = 1;   // k_s, filled in below
+if (tank_is_ballast != 0 || member_11_17 == 0) {
+    k_s = 1;
+} else {
+    k_s = 0.878 + (rho_g - rho_g_lo) / (rho_g_hi - rho_g_lo) * 0.122;
+    if (k_s < 0.878) k_s = 0.878;
+    else             k_s = k_s;
+    if (k_s > 1)     k_s = 1;
+    else             k_s = k_s;
+}
+var p_o = p_vp > 0 N/cm^2 ? Math.max(p_vp - 2.06 N/cm^2, 0 N/cm^2) : 0 N/cm^2;   // p_o = (p_vp - p_n), p_n = 2.06 N/cm^2
+var a_i = 0.71 * C_dp * (w_v * a_v + w_l * (l_tank / h_tank) * a_l + w_t * (b_tank / h_tank) * a_t);   // a_i (5.7.2(a))
+var theta_e = 0.71 * C_theta * theta;   // theta_e (5.7.2(b))
+var phi_e = 0.71 * C_phi * phi;   // phi_e (5.7.2(b))
+var zeta = Math.abs(b_tank / 2 - z);   // zeta, from the tank centre line
+var zeta_e = zeta;   // zeta_e, filled in below
+var eta_e = eta;   // eta_e, filled in below
+var dh_i = 0 m;   // delta h_i, filled in below
+if (phi_e < 0 deg && theta_e > 0 deg) {
+    // branch i of 5.7.2(b)
+    zeta_e = b_tank - zeta;
+    eta_e = eta;
+    dh_i = xi * Math.sin(-phi_e) + C_ru * (zeta_e * Math.sin(theta_e) * Math.cos(phi_e)
+         + eta_e * Math.cos(theta_e) * Math.cos(phi_e) - eta);
+} else {
+    // branch ii of 5.7.2(b)
+    zeta_e = zeta - delta_b;
+    eta_e = eta - delta_h;
+    dh_i = (l_tank - xi) * Math.sin(-phi_e) + C_ru * (zeta_e * Math.sin(theta_e) * Math.cos(phi_e)
+         + eta_e * Math.cos(theta_e) * Math.cos(phi_e) - eta);
+}
+var h_d = k_c * (eta * a_i / g + dh_i);   // h_d = k_c (eta a_i / g + delta h_i) (5.7.2(a))
+var p_s = k_s * rho_g * eta;   // static pressure (5.7.2(a))
+var p_d = k_s * rho_g * k_u * h_d;   // dynamic pressure (5.7.2(a))
+return p_s + p_d + p_o;   // the pressure the function returns