excel-codegen 0.9.0__py3-none-any.whl

excel_codegen/examples/abs_fpi/abs_fpi_internal.yaml

@@ -0,0 +1,266 @@
+version: 1
+
+# =========================================================================== #
+# 规则集：ABS FPI 5A-3-2/5.7   内压（internal tank pressures）
+#
+# 与 abs_fpi_external.yaml 同构，可以单独生成一本工作簿，也可以被 compose.py
+# 合进一本项目工作簿（与内压共用一张 Global 表）。
+#
+# 舱数据用**第三层作用域**（variables.group 成员表，指南 §18）：
+# 一个舱的参数只写一遍，工况用 local 变量 tank_ref 指向自己用哪个舱。
+# 成员表在独立工作表 "Tank Data" 里：一行一个舱、B 列起一个变量一列。
+# =========================================================================== #
+excel:
+  output: "abs_fpi_internal.xlsx"
+  sheets:
+    global: "Global Parameter"
+    local: "Local Parameter"
+    outputs: ["Code", "Summary"]
+
+variables:
+  # ------------------------------------------------------------------ 船 / 海 --
+  global:
+    - name: L
+      description: "规范船长 [m]（3-1-1/3）"
+      default: 340
+      type: float
+      suffix: " m"
+    - name: B
+      description: "规范型宽 [m]"
+      default: 64
+      type: float
+      suffix: " m"
+    - name: D
+      description: "规范型深 [m]"
+      default: 32
+      type: float
+      suffix: " m"
+    - name: d_f
+      description: "结构吃水 [m]；规范表按它给出，GM 四点表也按它插值"
+      default: 20.559
+      type: float
+      suffix: " m"
+    - name: rho_sea
+      description: "海水密度 [kg/m^3]"
+      default: 1025
+      type: float
+      suffix: " kg/m^3"
+    - name: g
+      description: "重力加速度 [m/s^2]"
+      default: 9.81
+      type: float
+      suffix: " m/s^2"
+    - name: Cb
+      description: "方形系数 [-]"
+      default: 1
+      type: float
+    - name: rho_g_lo
+      description: "k_s = 0.878 对应的 rho g [N/m^3]（1.005 N/cm^2-m = 10050）"
+      default: 10050
+      type: float
+      suffix: " N/m^3"
+    - name: rho_g_hi
+      description: "k_s 达到 1.0 的 rho g [N/m^3]（1.118 N/cm^2-m = 11180）"
+      default: 11180
+      type: float
+      suffix: " N/m^3"
+
+    # ------------------------------------------------------- 坐标映射（纯替换）--
+    - name: ax_x
+      description: "规范 x（自 AP）取自模型哪个轴：1=x, 2=y, 3=z"
+      default: "1"
+      type: string
+    - name: ax_y
+      description: "规范 y（自中线，向右舷为正）取自模型哪个轴"
+      default: "2"
+      type: string
+    - name: ax_z
+      description: "规范 z（基线以上，向上为正）取自模型哪个轴"
+      default: "3"
+      type: string
+    - name: neg_x
+      description: "规范 x 是否取反：填 - 表示取反，留空表示不取反"
+      default: ""
+      type: string
+    - name: neg_y
+      description: "规范 y 是否取反：填 - 表示取反，留空表示不取反"
+      default: ""
+      type: string
+    - name: neg_z
+      description: "规范 z 是否取反：填 - 表示取反，留空表示不取反"
+      default: ""
+      type: string
+
+  # ---------------------------------------------------------------- 每个工况 --
+  local:
+    - name: kind
+      description: "这个 Case 属于哪一套规则：EXT / INT（多规则集合用时用 case_filter 分流）"
+      default: "INT"
+      type: string
+    - name: draft
+      description: "工况吃水 [m]"
+      default: 8
+      type: float
+      suffix: " m"
+    - name: mu_deg
+      description: "浪向角 [deg]（Table 1A D）"
+      default: 90
+      type: float
+      suffix: " deg"
+    - name: k_c
+      description: "组合工况相关系数 [-]（Table 1A C）"
+      default: 1
+      type: float
+    - name: k_u
+      description: "载荷系数 [-]"
+      default: 1.1
+      type: float
+    - name: beta_VAC
+      description: "垂向加速度 ESF [-]（5A-3-A1/3）"
+      default: 0.660
+      type: float
+    - name: beta_LAC
+      description: "纵向加速度 ESF [-]"
+      default: 0.758
+      type: float
+    - name: beta_TAC
+      description: "横向加速度 ESF [-]"
+      default: 0.557
+      type: float
+    - name: beta_PMO
+      description: "纵摇 ESF [-]"
+      default: 0.686
+      type: float
+    - name: beta_RMO
+      description: "横摇 ESF [-]"
+      default: 0.492
+      type: float
+    - name: k_esf
+      description: "项目对 ESF 的分项系数 [-]"
+      default: 1.1
+      type: float
+    - name: w_v
+      description: "垂向加速度权重 [-]（Table 1A B）"
+      default: 0.75
+      type: float
+    - name: w_l
+      description: "纵向加速度权重 [-]"
+      default: 0.25
+      type: float
+    - name: w_t
+      description: "横向加速度权重 [-]"
+      default: 0.75
+      type: float
+    - name: xi
+      description: "舱长方向参考位置 [m]（5.7.2(b)）"
+      default: 21
+      type: float
+      suffix: " m"
+    - name: delta_b
+      description: "舱宽方向修正 [m]（5.7.2(b)）；无则填 0"
+      default: 0
+      type: float
+      suffix: " m"
+    - name: delta_h
+      description: "舱高方向修正 [m]（5.7.2(b)）；无则填 0"
+      default: 0
+      type: float
+      suffix: " m"
+
+    # ---- 指向成员表的"指针"：这个 Case 用哪个舱 ----
+    # 舱本身的 13 个参数不在这里，在下面的 variables.group（工作表 "Tank Data"）
+    - name: tank_ref
+      description: "本 Case 用哪个舱（成员表的行名）"
+      default: "WBT6"
+      type: string
+      choices: ["WBT6", "WBT7", "COT1"]
+
+  # ------------------------------------------------------------------ 每个舱 --
+  # 一行一个舱，只写一遍：多个工况用同一个舱时不再重复
+  group:
+    sheet: "Tank Data"
+    key: tank_ref
+    members: ["WBT6", "WBT7", "COT1"]
+    variables:
+    - name: rho_tank
+      description: "舱内液体密度 [kg/m^3]"
+      default: 1025
+      type: float
+      suffix: " kg/m^3"
+    - name: l_tank
+      description: "舱长 [m]（Fig. 11）"
+      default: 42
+      type: float
+      suffix: " m"
+    - name: b_tank
+      description: "舱宽 [m]（Fig. 11）"
+      default: 32
+      type: float
+      suffix: " m"
+    - name: h_tank
+      description: "舱高 [m]（Fig. 11）"
+      default: 32
+      type: float
+      suffix: " m"
+    - name: eta_deck
+      description: "甲板以上高度 [m]"
+      default: 0
+      type: float
+      suffix: " m"
+    - name: eta_overflow
+      description: "溢流管高度 [m]；Fig. 12 要求取 2/3 且不低于 760 mm"
+      default: 0
+      type: float
+      suffix: " m"
+    - name: C_dp
+      description: "舱形系数 [-]（5.7.2(d)）：矩形 1.0，J 形 0.7"
+      default: 1
+      type: float
+    - name: C_ru
+      description: "舱形系数 [-]（5.7.2(d)）：矩形 1.0，U 形 0.7"
+      default: 1
+      type: float
+    - name: p_vp
+      description: "真空/压力阀整定压力 [N/cm^2]"
+      default: 0
+      type: float
+      suffix: " N/cm^2"
+    - name: GM_full_in
+      description: "完整稳性 GM [m]；填 0 时用规范估算 0.12 B"
+      default: 0
+      type: float
+      suffix: " m"
+    - name: k_r_in
+      description: "横摇惯性半径 [m]；填 0 时用规范估算（0.35 B / 0.45 B）"
+      default: 0
+      type: float
+      suffix: " m"
+    - name: tank_is_ballast
+      description: "压载舱 1 / 货舱 0；决定 k_s 是否取 1"
+      default: 1
+      type: int
+    - name: member_11_17
+      description: "构件是否属于 5.7 表 3 的 11..17 项（1 = 是）"
+      default: 0
+      type: int
+
+templates:
+  - name: genie_int
+    description: "ABS FPI 5A-3-2/5.7 内压 -- GeniE 强度函数体"
+    output_sheet: "Code"
+    start_cell: "B2"
+    direction: "horizontal"
+    write_case_headers: true
+    engine: "excel"
+    filename: "{{ case_name }}.js"
+    template_file: "templates/int_body.js.j2"
+
+  - name: int_summary
+    description: "内压工况输入汇总（Markdown）"
+    output_sheet: "Summary"
+    start_cell: "B2"
+    direction: "vertical"
+    write_case_headers: true
+    engine: "excel"
+    filename: "{{ case_name }}_summary.md"
+    template_file: "templates/int_summary.md.j2"

excel_codegen/examples/abs_fpi/abs_fpi_internal_render.bat

@@ -0,0 +1,22 @@
+@echo off
+REM ===========================================================================
+REM  由 excel_codegen 生成 —— 改完参数双击本文件即可把结果写回 Output 表。
+REM  重新生成工作簿（init）时会一并覆盖本文件。
+REM ===========================================================================
+cd /d "%~dp0"
+
+where uv >nul 2>nul
+if %errorlevel%==0 (
+  uv run excel-codegen render -c "abs_fpi_internal.yaml" -x "abs_fpi_internal.xlsx" --write-excel
+) else (
+  excel-codegen render -c "abs_fpi_internal.yaml" -x "abs_fpi_internal.xlsx" --write-excel
+)
+
+echo.
+if errorlevel 1 (
+  echo [失败] 上面有报错信息。常见原因：依赖没装（跑一次 setup.sh / uv sync）、
+  echo        或者 Excel 正开着这个文件（先关掉再试）。
+) else (
+  echo [完成] 回到 Excel 打开「Output」表看结果。
+)
+pause

excel_codegen/examples/abs_fpi/abs_fpi_internal_render.sh

@@ -0,0 +1,23 @@
+#!/usr/bin/env bash
+# ===========================================================================
+#  由 excel_codegen 生成 —— 改完参数跑一次本文件即可把结果写回 Output 表。
+#  重新生成工作簿（init）时会一并覆盖本文件。
+# ===========================================================================
+set -uo pipefail
+cd "$(dirname "$0")"
+
+if command -v uv >/dev/null 2>&1; then
+  uv run excel-codegen render -c "abs_fpi_internal.yaml" -x "abs_fpi_internal.xlsx" --write-excel
+else
+  excel-codegen render -c "abs_fpi_internal.yaml" -x "abs_fpi_internal.xlsx" --write-excel
+fi
+status=$?
+
+echo
+if [ "$status" -ne 0 ]; then
+  echo "[失败] 上面有报错信息。常见原因：依赖没装（跑一次 ./setup.sh 或 uv sync）、"
+  echo "       或者 Excel / WPS 正开着这个文件（先关掉再试）。"
+else
+  echo "[完成] 回到 Excel 打开「Output」表看结果。"
+fi
+exit "$status"

excel_codegen/examples/abs_fpi/compose.yaml

@@ -0,0 +1,20 @@
+# =========================================================================== #
+# 项目工作簿的配方：把哪些规则集合进一本工作簿、各自用什么 case_filter 分流。
+#
+#   python compose.py   →   abs_fpi.yaml   →   init / render
+#
+# 变量与模板的真源仍在各规则集 YAML 里；compose 会把它们合并、检查冲突，
+# 再给每套模板加上 case_filter（只在多于一个规则集时才加）。
+# =========================================================================== #
+project: "abs_fpi.yaml"
+output: "ABS_FPI_load_cases.xlsx"
+
+rulesets:
+  # 外压放前面：case_filter 模式下"横向相对列"要在每个规则集内部保持连续，
+  # 同一规则集的工况必须相邻（见 README「工况顺序」）。
+  - file: "abs_fpi_external.yaml"
+    set: "EXT"
+    filter: "kind == 'EXT'"
+  - file: "abs_fpi_internal.yaml"
+    set: "INT"
+    filter: "kind == 'INT'"

excel_codegen/examples/abs_fpi/generated/COT1-d20.559-mu60.js

@@ -0,0 +1,124 @@
+// COT1-d20.559-mu60
+// 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 = 60 deg;
+var k_c = 0.5;
+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: COT1
+var rho_tank = 900 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.21 N/cm^2;
+var GM_full_in = 0 m;
+var k_r_in = 0 m;
+var tank_is_ballast = 0;
+var member_11_17 = 1;
+// --- 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/COT1-d20.559-mu60_summary.md

@@ -0,0 +1,30 @@
+| 项目 | 值 | 出处 |
+| --- | --- | --- |
+| case_name | COT1-d20.559-mu60 | 工况名（Excel 表头） |
+| kind | INT | 规则集（多规则集共用一本工作簿时用 case_filter 分流） |
+| tank_ref | COT1 | 指向成员表 `Tank Data` 的哪一行（指南 §18） |
+| draft | 20.559 m | 装载工况 |
+| mu_deg | 60 deg | Table 1A D |
+| k_c | 0.5 | 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 | 900 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.21 N/cm^2 | 成员表 Tank Data |
+| GM_full_in / k_r_in | 0 m / 0 m | 稳性计算书 |
+| tank_is_ballast / member_11_17 | 0 / 1 | 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/EXT-T15-mu0-kf-1.js

@@ -0,0 +1,91 @@
+// EXT-T15-mu0-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 = 0 deg;
+var k_c = 0.5;
+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-mu0-kf-1_summary.md

@@ -0,0 +1,23 @@
+| 项目 | 值 | 出处 |
+| --- | --- | --- |
+| case_name | EXT-T15-mu0-kf-1 | 工况名（Excel 表头） |
+| kind | EXT | 规则集（多规则集共用一本工作簿时用 case_filter 分流） |
+| draft | 15 m | 装载工况 |
+| mu_deg | 0 deg | Table 1A D |
+| k_c | 0.5 | 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 | ( /  / ) | 模型 |