dsh-engineering-ui 0.0.0-stage → 1.5.0

package/engineering/skills/physics-in-loop/SKILL.md

@@ -0,0 +1,238 @@
+---
+name: physics-in-loop
+description: >-
+  物理仿真驱动的结构设计闭环（Physics-in-the-Loop）。
+  从载荷工况 JSON 出发，自动进行参数化建模、有限元验证、迭代优化，
+  直到满足安全系数、位移、体积等约束。
+  触发关键词：力学验证、强度校核、安全系数、FEA、应力分析、拓扑优化、
+  载荷工况、静力分析、迭代设计、结构验证。
+whenToUse: >-
+  当用户要求对零件进行力学/结构验证、指定了力/压力/约束条件、
+  或询问"这个零件在载荷下会不会坏"时使用。
+  如果用户未提供材料、约束或载荷，必须先追问再开始。
+disable-model-invocation: false
+user-invocable: true
+source: engineering
+provider: filesystem
+---
+
+# Physics-in-the-Loop 物理验证工作流
+
+> 基于 IJCAI 2026 论文 "Physics-in-the-Loop: A Hybrid Agentic Architecture
+> for Validated CAD Engineering Design" 的工程化实现。
+
+## 核心思想
+
+LLM 擅长生成设计方案，但不擅长计算物理真实性。
+本系统把**确定性物理工具**（解析解/FEA）嵌入设计决策循环：
+`LLM 生成 → 工具验证 → 反馈修正 → 重新生成`，直到满足安全门槛。
+
+```
+用户工程需求
+    ↓
+结构化载荷工况 JSON（load_case）
+    ↓
+[Generate] LLM 提出设计参数（厚度、圆角、筋板…）
+    ↓
+[Simulate] 解析解 / Gmsh+CalculiX / SolidWorks Simulation
+    ↓
+[Refine]   若 SAFETY_FACTOR ∉ [min, max] → 生成修正建议
+    ↓
+循环直到 PASS 或达到 max_iterations
+    ↓
+生成 release_manifest（含局限性声明 + 人工复核要求）
+```
+
+## 支持的求解等级（自动降级）
+
+| 等级 | 方法 | 依赖 | 适用场景 |
+|------|------|------|---------|
+| L0 | 解析解（悬臂梁/简支梁/扭转轴） | 无 | 快速初判、回归测试 |
+| L1 | 杆系 truss 求解 | 无（预留接口） | 简单框架结构 |
+| L2 | Gmsh + CalculiX | gmsh + ccx | 一般 3D 零件 |
+| L3 | SolidWorks Simulation API | SW Simulation 授权 | 生产级验证 |
+
+第一版默认使用 L0，检测到 Gmsh/CalculiX 后自动升级。
+
+## 命令速查
+
+```powershell
+# 1. 先运行演示，理解工作流
+python "<工程模式根目录>\tools\sw_bridge.py" physics-demo
+
+# 2. 检查可用求解器
+python "<工程模式根目录>\tools\sw_bridge.py" physics-status
+
+# 3. 校验你的载荷工况文件
+python "<工程模式根目录>\tools\sw_bridge.py" physics-validate-case my_load_case.json
+
+# 4. 自动迭代优化（最多 5 轮）
+python "<工程模式根目录>\tools\sw_bridge.py" physics-optimize my_load_case.json --max-iter 5
+
+# 5. 查看某轮结果
+python "<工程模式根目录>\tools\sw_bridge.py" physics-report opt_iter_3
+
+# 6. 获取下轮修正建议
+python "<工程模式根目录>\tools\sw_bridge.py" physics-recommend opt_iter_3
+```
+
+## 载荷工况文件格式（load_case.json）
+
+```json
+{
+  "schema_version": "1.0",
+  "meta": {
+    "problem_id": "L_BRACKET_001",
+    "title": "L 型安装支架",
+    "revision": "A"
+  },
+  "units": { "length": "mm", "force": "N", "stress": "MPa", "mass": "kg" },
+
+  "design_domain": {
+    "bounds": {
+      "x_min": 0, "x_max": 120,
+      "y_min": 0, "y_max": 80,
+      "z_min": 0, "z_max": 60
+    }
+  },
+
+  "material": {
+    "id": "AL_6061_T6",
+    "youngs_modulus_mpa": 68900,
+    "poissons_ratio": 0.33,
+    "yield_strength_mpa": 276,
+    "density_kg_m3": 2700
+  },
+
+  "spatial_selectors": [
+    {
+      "id": "fixed_end",
+      "type": "box",
+      "bounds": { "x_min": 0, "x_max": 10, "y_min": 0, "y_max": 80, "z_min": 0, "z_max": 60 },
+      "selection_rule": "faces_intersecting_region"
+    },
+    {
+      "id": "load_face",
+      "type": "box",
+      "bounds": { "x_min": 110, "x_max": 120, "y_min": 20, "y_max": 60, "z_min": 10, "z_max": 50 },
+      "selection_rule": "faces_intersecting_region"
+    }
+  ],
+
+  "boundary_conditions": [
+    {
+      "id": "bc_fixed",
+      "spatial_selector_id": "fixed_end",
+      "type": "fixed_displacement",
+      "dof_lock": { "x": true, "y": true, "z": true, "rx": true, "ry": true, "rz": true }
+    }
+  ],
+
+  "loads": [
+    {
+      "id": "payload",
+      "spatial_selector_id": "load_face",
+      "type": "distributed_force",
+      "magnitude_n": 1500,
+      "direction": [0, 0, -1]
+    }
+  ],
+
+  "acceptance": {
+    "analysis_type": "linear_static",
+    "min_safety_factor": 2.0,
+    "target_safety_factor_max": 5.0,
+    "max_displacement_mm": 0.5,
+    "min_wall_thickness_mm": 3.0
+  },
+
+  "optimization": {
+    "primary": "minimize_mass",
+    "secondary": ["minimize_max_displacement"]
+  }
+}
+```
+
+## 设计参数格式（design_parameters.json）
+
+```json
+{
+  "thickness_mm": {
+    "id": "thickness_mm",
+    "value": 5.0,
+    "min": 2.0,
+    "max": 15.0,
+    "step": 0.5
+  },
+  "height_mm": {
+    "id": "height_mm",
+    "value": 60.0,
+    "min": 30.0,
+    "max": 150.0,
+    "step": 5.0
+  },
+  "fillet_mm": {
+    "id": "fillet_mm",
+    "value": 2.0,
+    "min": 0.5,
+    "max": 10.0,
+    "step": 0.5
+  },
+  "mount_hole_diameter_mm": {
+    "id": "mount_hole_diameter_mm",
+    "value": 10.0,
+    "locked": true,
+    "source": "interface_requirement"
+  }
+}
+```
+
+**注意**：`locked: true` 的参数不会被自动优化修改，用于保护接口尺寸。
+
+## Agent 必须遵守的流程规则
+
+1. **载荷工况缺失时不得开始仿真**
+   - 缺少材料屈服强度 → 必须询问用户或标 TBD
+   - 缺少固定约束位置 → 必须询问
+   - 缺少载荷大小/方向 → 必须询问
+
+2. **仿真通过后必须给出 release_manifest**
+   ```json
+   {
+     "status": "CANDIDATE_AFTER_HUMAN_REVIEW",
+     "passed_gates": ["CAD_BUILD", "DESIGN_SPACE", "MESH_QUALITY", "STATIC_FEA", "SAFETY_FACTOR"],
+     "not_evaluated": ["fatigue", "buckling", "contact_non_linearity"],
+     "human_review_required": true
+   }
+   ```
+
+3. **禁止的行为**
+   - 不得在无载荷数据时声称"结构安全"
+   - 不得修改 locked 参数
+   - 不得无限迭代（默认上限 5 次，最多 10 次）
+   - 不得将"解析解通过"等同于"FEA 通过"——两者结论应一致标注
+
+4. **何时停止迭代**
+   - `overall_status == "PASS"` → 通过
+   - 达到 `max_iterations` → REVIEW_REQUIRED
+   - 连续 2 轮无改善 → STALLED
+   - 网格化/求解连续失败 → UNSIMULATABLE
+   - 工况/材料/约束缺失 → BLOCKED_BY_INPUT
+
+## 可用材料速查
+
+```
+铝合金：AL_6061_T6 (σy=276MPa), AL_7075 (σy=503MPa), AL_5052 (σy=193MPa)
+结构钢：ST_API_S235 (σy=235MPa), ST_API_A36 (σy=250MPa), ST_API_4140 (σy=655MPa)
+不锈钢：ST_STAINLESS_304 (σy=205MPa), ST_STAINLESS_316 (σy=170MPa)
+钛合金：TI_GRADE5 (σy=880MPa)
+塑料：PC_POLYCARBONATE (σy=65MPa), ABS (σy=45MPa)
+```
+
+完整列表：`material_db.list_materials()`
+
+## 参考文档
+
+- `references/gbt-drafting.md` — GB/T 机械制图国家标准
+- `references/autocad-workflow.md` — AutoCAD 工作流与 CADX 验证
+- `../tools/physics/` — 完整源码（load_case / fea_solver / refine_rules 等）

package/engineering/skills/physics-in-loop/references/feadocs.md

@@ -0,0 +1,142 @@
+# Physics-in-the-Loop 参考文档
+
+## 1. 载荷工况 Schema（load_case_schema.md）
+
+### 字段说明
+
+| 字段 | 类型 | 必填 | 说明 |
+|------|------|------|------|
+| `schema_version` | string | 是 | 当前版本 `"1.0"` |
+| `meta.problem_id` | string | 是 | 问题唯一标识（大写下划线） |
+| `meta.title` | string | 否 | 人类可读标题 |
+| `units` | object | 否 | 单位约定（默认 mm/N/MPa/kg） |
+| `design_domain.bounds` | object | 是 | 6 个边界坐标（mm） |
+| `material.id` | string | 是 | 材料 ID（见 material_db.py） |
+| `material.youngs_modulus_mpa` | float | 是 | 弹性模量 |
+| `material.yield_strength_mpa` | float | 是 | 屈服强度 |
+| `material.density_kg_m3` | float | 是 | 密度 |
+| `spatial_selectors[]` | array | 是 | 区域选择器（固定面/载荷面） |
+| `spatial_selectors[].id` | string | 是 | 选择器唯一 ID |
+| `spatial_selectors[].type` | string | 是 | `box` / `cylinder` / `sphere` |
+| `spatial_selectors[].bounds` | object | box 时必填 | 6 个边界坐标 |
+| `boundary_conditions[]` | array | 是 | 约束列表 |
+| `loads[]` | array | 是 | 载荷列表 |
+| `acceptance.min_safety_factor` | float | 是 | 最小安全系数（通常 2.0） |
+| `acceptance.target_safety_factor_max` | float | 是 | 最大安全系数（通常 5.0） |
+| `acceptance.max_displacement_mm` | float | 否 | 最大允许位移 |
+
+### 单位约定
+
+所有坐标、尺寸单位为 **mm**，力为 **N**，应力为 **MPa**（= N/mm²），质量为 **kg**。
+求解器内部自动完成单位换算。
+
+---
+
+## 2. 材料策略（material_policy.md）
+
+### 材料来源优先级
+
+1. **用户指定**：用户提供材料牌号或屈服强度 → 直接使用
+2. **material_db 匹配**：ID 匹配内置材料库 → 使用库值
+3. **TBD 标记**：均不可得 → 标 TBD，禁止给出确定性安全结论
+
+### 安全系数与材料的关系
+
+| 材料类型 | 推荐 min SF | 推荐 max SF | 备注 |
+|---------|------------|------------|------|
+| 结构钢（静载） | 1.5 | 3.0 | 延性材料，失效前有明显变形 |
+| 铝合金（静载） | 2.0 | 4.0 | 需注意腐蚀和疲劳 |
+| 钛合金 | 2.0 | 4.0 |  aerospace 常用 |
+| 工程塑料 | 3.0 | 6.0 | 蠕变敏感，保守设计 |
+| 铸铁 | 3.0 | 5.0 | 脆性材料，SF 需更大 |
+| 疲劳载荷 | 按 S-N 曲线 | — | 本系统暂不支持疲劳分析 |
+
+### 不支持的材料行为
+
+以下情况必须声明 `NOT_EVALUATED`：
+- 塑性大变形（σ > σy 后行为）
+- 疲劳寿命预测
+- 屈曲（薄壁件失稳）
+- 冲击/动态载荷
+- 热-结构耦合
+- 接触非线性（螺栓预紧、摩擦）
+
+---
+
+## 3. 求解有效性（fea_validity.md）
+
+### 解析解适用条件
+
+解析解（Euler-Bernoulli 梁理论）仅在以下条件下可靠：
+1. 长细比 L/h > 10（梁近似成立）
+2. 材料线弹性（σ < 0.5×σy）
+3. 小变形（δ/L < 0.01）
+4. 均匀截面，无突变
+
+**不满足上述条件时，解析解结果仅作参考，必须升级求解器。**
+
+### 网格质量门控
+
+| 指标 | 合格阈值 | 处理方式 |
+|------|---------|---------|
+| 单元最小 Quality | ≥ 0.1 | 失败则细化网格 |
+| 单元数量 | 1,000 ~ 1,000,000 | 超出范围需人工审查 |
+| 雅可比行列式 | > 0 | 负值表示单元反转，网格无效 |
+| 自由边界数 | = 0（封闭体） | 开放体需确认设计意图 |
+
+### 收敛判定
+
+连续两轮（n, n+1）同时满足：
+```
+|SF_n - SF_{n+1}| / SF_n < 5%   AND   |V_n - V_{n+1}| / V_n < 3%
+```
+则视为收敛，即使 SF 仍在目标区间外也建议人工介入。
+
+---
+
+## 4. 修正策略（refinement_policy.md）
+
+### 修正动作优先级
+
+```
+1. 安全系数过低（SF < min）
+   → 优先增加壁厚（最有效，对 I 影响最大）
+   → 次选增大圆角（降低应力集中）
+   → 最后加筋（增加局部刚度）
+
+2. 安全系数过高（SF > max）
+   → 减小壁厚（在 min 之上）
+   → 减小截面高度（若位移允许）
+
+3. 位移超限
+   → 增大截面高度（I ∝ h³，效果最显著）
+   → 增加支撑点（改变边界条件）
+   → 增大壁厚（次选）
+
+4. 设计空间越界
+   → 拒绝本轮，缩小参数范围重新搜索
+   → 不允许越界设计通过
+```
+
+### 禁止的修正动作
+
+- 修改 locked 参数（接口尺寸、安装孔径）
+- 改变设计空间边界
+- 改变材料（除非用户明确允许）
+- 添加未声明的支撑或约束
+- 减少实体特征（如去掉加强筋）以降低质量但牺牲强度
+
+### 修正透明度
+
+每次修正必须记录：
+```json
+{
+  "change_id": "iter_03_change_01",
+  "param_id": "thickness_mm",
+  "from": 5.0,
+  "to": 7.0,
+  "reason": "safety factor 1.2 < 2.0; thickness scales ~sqrt(2.0/1.2)",
+  "expected_effect": "SF ≈ 2.1, displacement ≈ -30%",
+  "verified": true
+}
+```