从一块平板学会 HyperStudy：厚度、材料、DOE、Fit、优化与鲁棒性

42分钟前浏览0

一篇文章把hyperstudy的流程走完，由于笔者用word写的，直接去上传文档，图片被压缩了。

软件界面为 HyperStudy 2026；

结构求解由 OptiStruct 2025 完成。

文中所有数值均来自本案例的实际求解文件，单位体系为 N、mm、tonne，报告质量换算为 kg。

先看问题：怎样让板足够轻，同时不挠得太厉害？

我们建立一块 300 mm × 200 mm 的矩形钣金板，网格尺寸为 10 mm，采用 CQUAD4 壳单元。左边界固定，自由边中点的 GRID 341 施加 1000 N、沿 −Z 方向的集中力。线性静力求解给出各节点位移。设计要求是加载点合位移不大于 2 mm，目标是质量最轻。

只有两个真正需要改变的工程属性：壳属性厚度，以及材料。厚度是连续量；钢、6061-T6 铝合金和 AZ31B 镁合金是离散选项。因此我们为三种材料各建一个 HyperStudy Study，分别优化厚度，最后比较三条最优解。这样每个 Study 内的求解器输入、密度和弹性模量都明确，不会把“材料编号”错误地当作连续物理量插值。

图 1　平板工程在 HyperStudy 中的 DOE Integrity。max_node 一行的 Range = 0 以红色显示，表示这 11 次计算的最大位移节点编号没有变化；它不表示位移为零，也不表示求解失败，最大位移值始终是加载位置。

1　建立并检查有限元模型

先在 HyperMesh 中建立板的几何和壳网格，定义 PSHELL 厚度、MAT1 材料、固支 SPC 和集中力，再输出 FEM。在线性静力分析里，求解器处理的是。这里是厚度，是弹性模量，是泊松比；位移随厚度和材料刚度改变。读者复现时应先直接求解一次原始 FEM，确认约束没有造成刚体运动、加载节点正确、OUT 文件无致命错误、PCH/H3D 中有可提取的位移。OptiStruct-线性静力学分析理论基础

线性静力分析：手算有限元矩阵

三种材料参数如下。密度表中的 tonne/mm³ 与 FEM 单位体系一致，报告质量乘 1000 得 kg。

材料	弹性模量 E / MPa	泊松比	密度 / tonne·mm⁻³	在本教程中的作用
钢	210000	0.30	7.85×10⁻⁹	刚度高、密度也高
6061-T6 铝合金	69000	0.33	2.70×10⁻⁹	中等刚度、较轻
AZ31B 镁合金	45000	0.35	1.77×10⁻⁹	弹性模量较低、密度最低

壳板质量可用近似理解；面积固定，所以同一材料下质量随厚度线性增加。位移则随弯曲刚度快速变化，弯曲刚度与有关。这说明“减一点厚度”对位移的影响可能远大于对质量的影响；具体数值仍以有限元求解为准。

2　Define Models：让 HyperStudy 真正驱动 FEM

在 HyperStudy 的 Setup → Define Models 中建立求解模型，指定 OptiStruct 入口 FEM、运行命令和工作目录。输入模板中的 PSHELL 厚度位置要输出实际数值；仅声明一个 parameter(...) 并不会自动改写有限元卡片。每次 Evaluate 时，HyperStudy 创建独立运行目录，写入该设计的 FEM，调用求解器，再从结果文件提取响应。

图 2　左侧点 Setup → Definition → Define Models。右侧一行模型中，Type 是 Parameterized File，Resource 指向厚度模板，Solver Execution Script 为 Python。操作时先关联模板和运行脚本，再点下方 Import Variables。

本工程的模板 plate_input.tpl 实际只有三行：第一行 thickness 是输出 CSV 的列名；第二行 {parameter(thickness, "Plate Thickness", 5.0, 0.5, 15.0)} 声明可导入参数；第三行 {thickness, %16.8f} 将当前设计厚度写入plate_input.csv。Python 适配脚本读取该 CSV，重建 PSHELL 的厚度字段、运行 OptiStruct，并把响应写到 responses.csv。所以这里采用的是“模板 CSV → 适配脚本 → FEM”的链路；在 HyperStudy 里看到变量不等于 FEM 已被改动，必须检查单次运行目录中的实际 FEM。

在 Define Input Variables 中将板厚 thickness 设为 Numeric，填写下界、Nominal 和上界。这里的界限是可搜索设计域，并非求解器单位或制造公差。材料在三个 Study 中分别固定。先做 Test Models / Nominal Run：钢板 5 mm 的实算结果是质量 2.355 kg、加载点位移 19.90701 mm。这远超 2 mm 约束，但说明变量替换、求解和响应提取链路已经打通。Nominal 超限不等于流程失败。

图 3　Setup → Definition → Define Input Variables。钢板工程的 thickness 下界 0.5 mm、Nominal 5 mm、上界 15 mm。铝板和镁板采用各自更大的上界，以覆盖位移约束附近。Active 勾选表示该变量参与计算。

厚度输入行局部放大：Varname=thickness，Lower Bound=0.5，Nominal=5，Upper Bound=15，Data Type=Numeric。

在 Define Output Responses 中至少创建两个数值响应：mass（kg）和 load_disp（GRID 341 的合位移，mm）。max_disp 可以作为全模型校核量。加载点编号 GRID 341 本身固定，不应成为优化目标或约束：它是定位信息，不随设计改变。图 1 的红色 Range = 0 正说明这一点。若某次改造后最大位移位置可能迁移，则保留全局最大位移数值作为约束，比保存节点 ID 更有物理意义。

图 4　Setup → Definition → Define Output Responses。Value 列显示 Nominal 的提取结果：mass=2.355、load_disp=max_disp=19.90701、max_node=341。max_node 是早期保留的诊断列，正式优化仅用 mass 和 load_disp。

图 5　同一页面切换顶部 Data Sources Table。四行 ds_mass 等数据源都读取当前运行目录 m_1/responses.csv，Tool Settings 选对应 CSV 字段；其上层 Response Expression 使用 max(ds_mass) 等标量表达式。本例每个 CSV 只有一行，所以 max 返回该行数值。

读者手工定义响应时，顺序是：①在 Data Sources 中点 Add Data Source，选运行文件 run-file:///m_1/responses.csv；②在 Tool Settings 中指定字段，如 loading_point_displacement_mm；③返回 Define Output Responses，点 Add Output Response，取 Varname load_disp，Expression 指向 ds_load_disp；④点 Evaluate，核对 Value 与 CSV 原值一致。如果文件路径指向 Study 根目录的旧结果，Nominal 也可能显示一个“合理数字”，但每个 DOE 点读到的其实是同一份数据。

什么是“变量”，什么是“响应”？

输入变量是 HyperStudy 主动修改的量。这里的 thickness 由我们指定，并通过模板传给 FEM。响应是模型求解之后才知道的量：求解器计算出位移，适配脚本得到质量与位移，HyperStudy 从 CSV 读出 mass、load_disp。一个响应不是一个按钮，也不是节点集合，它必须有明确的数据来源、计算规则和单位。load_disp=19.90701 mm 表示名义钢板在指定加载点的合位移数值；max_node=341 只表示这个最大点的位置编号。优化目标和约束最终都作用在响应数值上。

把链路写成一行最容易检查：thickness → plate_input.csv → PSHELL 厚度 → OptiStruct → PCH 位移与质量计算 → responses.csv → ds_load_disp → load_disp → 2 mm 约束。如果响应不随厚度变化，沿这条链从右向左查：先看 Summary 是否恒定，再看每次的 responses.csv，最后比对每次 FEM 的 PSHELL 厚度。

3　DOE：先认识厚度与响应的关系

每种材料做 11 个全域 DOE 设计点。DOE 是 Design of Experiments，即按规定的采样方案批量改变厚度，调用真实求解器，得到厚度—质量—位移数据。进入 Post-Processing 的 Integrity，先看 Points、No Values 和 Bad Values：只有设计点完成且响应有效，后续拟合才有意义。再看 Scatter：横轴厚度、纵轴位移；厚度增大，质量上升而位移下降。材料不同会形成三条不同曲线。

DOE 的厚度上下界须覆盖 2 mm 约束附近。过窄的范围可能根本看不到可行点；过宽则可能把非常薄、强非线性的区间和最优附近混在同一个拟合模型中。建议根据 DOE 图先找到位移从大于 2 mm 过渡到小于 2 mm 的位置，再决定是否增加局部样本。

在界面中看结果：左侧选 DOE → Post-Processing；顶部点 Integrity，确认 mass、load_disp 的 Points 为 11，No Values 和 Bad Values 为 0；再点 Summary Table，看每一行厚度、质量与位移，而不是只看红绿标记。Scatter Table 可把 X Axis 设为 thickness、Y Axis 设为 load_disp，查看厚度增大时的位移趋势。

“正相关”和“负相关”应该怎样读？

若 DOE 中厚度大的点通常对应质量大，就称厚度与质量正相关；若厚度大的点通常对应位移小，就称厚度与位移负相关。钢板 11 点中，厚度与质量的样本 Pearson 相关系数约 +1.000，厚度与加载点位移约 −0.557。相关系数的范围是 −1 到 +1：符号表示同向或反向，绝对值描述这批样本的线性关系强弱。这里位移虽然随厚度单调下降，却不是直线，极薄区还有巨大位移，所以线性相关系数并不接近 −1；这不表示厚度对位移影响小。

读图时先检查横轴、纵轴和单位。Scatter 默认可能选中两个输入变量，图上有点并不代表已经画出了“厚度—位移”；需在右侧 X Axis 选 thickness，Y Axis 选 load_disp。Linear Effects 的斜线表示所选输入与响应的总体变化方向，Pareto Plot 的条纹/颜色表达效应方向和大小。只有一个输入变量时，Interactions 页不会给出真正的双变量交互结论，因为交互至少需要两个可变输入。相关性描述采样结果；在前罩那样的多变量案例中，不能直接把它当成“单独改某厚度就一定如此”的因果规律。

图 6　钢板 11 个 DOE 点逐行列出厚度和位移。厚度从约 1.16 mm 到 14.34 mm，位移从约 1597.93 mm 降到 0.8438 mm。这里极薄板的位移跨度巨大，直接给全域拟合增加了难度。

4　Fit：响应面是近似模型，必须验证

Fit 将有限个求解器样本拟合为和。它可快速预测新厚度，但必须检查模型质量。钢板全域 Fit 的质量训练，加载点位移的只有 0.4092176：质量随厚度线性变化，很好拟合；位移跨越从亚毫米到上千毫米，单个全域响应面表现差。随后又在最优附近做了 13 点局部 DOE，利用 25 个独立求解点检查 Fit。三种材料的加载点位移最大测试误差分别为：钢 8.83×10⁻⁷ mm、铝 1.33×10⁻⁶ mm、镁 1.50×10⁻⁶ mm。这些数字只证明该局部范围的预测精度，不能外推到未经测试的极薄区域。

图 7　左侧选 Fit from DOE → Post-Processing，顶部选 Diagnostics。mass 的 R² 为 1.000，而 load_disp 为 0.4092176；max_node 恒定，R² 显示 NaN，因此也不应拿它做有意义的响应面。可再点 Residuals 检查误差集中在哪些厚度。

复现时应在 Fit → Post-Processing 中同时看残差、训练拟合以及独立点误差。若响应面用于筛选候选厚度，最终仍应将候选设计回代到 OptiStruct。否则近似模型在约束边界处的一点偏差，就可能把超限设计误判为可行。

5　Optimization：在位移约束下找最轻厚度

Optimization 中定义目标为最小化 mass，约束为 load_disp ≤ 2 mm，设计变量为厚度。案例采用 GRSM 进行搜索，并对每种材料完成 18 次求解器评估。最轻的实算可行设计如下：

图 8　左侧选 GRSM minimum mass → Specifications。方法为 Global Response Search Method，Number of Evaluations 为 18。这里控制搜索预算；目标和约束在该 Approach 的 Definition → Objectives/Constraints 中定义。

目标行选择 Response=mass、Type=Minimize；约束行选择 Response=load_disp、Bound Type=≤、Bound Value=2.0。目标和约束都由前述 CSV 数据源得到，单位必须与上限一致。设置好后再 Evaluate。若只在 Setup 创建了响应，却没有在 Optimization 的目标/约束页引用它们，算法不会自动知道“要轻”和“位移不能超限”。

目标/约束表局部放大：两行 Active 均被勾选。Minimum mass 引用 mass，Loading point displacement limit 引用 load_disp，上界填 2.000000 mm。

图 9　左侧选 GRSM minimum mass → Post-Processing，再点顶部 Summary。每一行是一轮评估。红色 Violated 为位移超限，Feasible 为满足约束；绿色一行是当前显示的最优可行点。该钢板工程中绿色行厚度约 10.7636 mm、质量约 5.0697 kg、位移约 1.9955 mm。

顶部 Optima Table 可以只看最优点，并在右侧 X Axis/Y Axis 选择厚度、质量或位移；若要判断收敛过程，还应回到 Summary 看各轮设计是否反复越界。最后点击某一行的 Show in Study Directory，在对应 run__xxxxx/m_1 中打开 FEM、OUT、PCH 和 responses.csv，确认最优点经过求解器回代。只截一张最优值图不足以判断约束是否满足。

固定材料	最优厚度 / mm	质量 / kg	加载点位移 / mm	与 2 mm 上限的距离
钢	10.7636	5.0697	1.9955	0.0045 mm
6061-T6 铝合金	15.5390	2.5173	1.9991	0.0009 mm
AZ31B 镁合金	17.8796	1.8988	1.9981	0.0019 mm

镁合金在这组三种材料、这套载荷与位移约束下质量最低。它需要更厚的板，因为弹性模量更低；低密度仍使它取得最小质量。请勿把这个结果直接当作材料选型结论：屈服强度、局部屈曲、疲劳、连接、成形、成本和腐蚀均未纳入本次目标与约束。

6　Stochastic：名义最优不等于制造上稳健

在三个最优厚度附近，各做 25 个厚度 ±5% 的随机样本，并逐点调用求解器。三组结果都是 13/25，即 52% 满足 2 mm 位移限制。原因很直观：质量最小化会把设计推到位移边界附近；一旦实际板厚落在名义值以下，就容易超限。这里的 52% 是该 25 点样本的通过率，不是已经证明的总体可靠度。

若要考虑镁板厚度 ±5% 公差，可以把名义厚度提高到 18.83 mm。实算名义质量 1.999746 kg、位移 1.710575 mm；最薄边界 mm 的实算位移为 1.995089 mm。这只验证了当前模型中“厚度偏差”这一项的最薄边界；载荷、材料离散性和边界条件变化未纳入。选择 18.83 mm 是从“最低名义质量”走向“可承受所设厚度公差”的工程取舍。

复现时最容易漏掉的检查

按 Define Models → Nominal Run → DOE → Fit → Optimization → Stochastic 逐级走。每一级至少核对三件事：运行目录中存在求解器结果；HyperStudy 响应无缺失或 Bad Values；关键候选点由真实求解器回代。遇到 Integrity 中红色零 Range，先看它对应的是位移数值还是恒定编号，再判断是否异常。

工程入口：钢板 Study、铝板 Study、镁板 Study。三套工程各自保留运行目录与响应文件，便于逐点复核。

7　逐个打开结果 Table：平板工程读图手册

前文说明了求解链路和最终值。下面按左侧树节点整理其余 Table。配图均取自本工程 HyperStudy 界面；截图保留左侧树与顶部页签，移除了下方日志区。图中若保留软件默认的坐标轴选择，我会写明应如何切换，避免读者把默认图当成最终分析图。

DOE → Post-Processing：九个页签

Integrity 已见图 1：Points=11 表示 11 组样本，No Values=0 表示没有空值，Bad Values=0 表示没有坏值。max_node 的 Unique=1、Range=0，仅说明节点编号恒定。Summary 已见图 6：按行读取厚度、质量和位移，这是所有图形解释的原始依据。

Parallel Coordinate Table：一条折线对应一个 DOE 点。各竖轴刻度独立，线在相邻轴间向上或向下，只代表该样本在两列中的相对位置；不要把倾斜角度当成力学导数。可用右侧 Channel 选择显示 thickness、mass、load_disp。

Scatter Table：当前图的横轴是板厚、纵轴是质量，点落在上升直线上，直观显示正相关。右侧点 Y Axis 的 load_disp 后，才是“厚度—位移”散点；此时应看到整体下降且极薄区曲率很大。底部 Correlation 可查看相关性矩阵。

Pareto Plot Table：横向比较输入变量对所选响应的效应大小。当前只有一个厚度变量，因此只有一根柱；柱子向上表示对当前 mass 响应的正向效应。若切换响应为 load_disp，效应方向应反向。柱高的刻度是该图的统计敏感度指标，不能直接当作 kg/mm。

Linear Effects Table：蓝线随厚度上升，展示 thickness → mass 的主趋势。红点是实际 DOE 值，线是概括趋势。对位移应另选 load_disp，并结合 Scatter 判断它是否近似直线；本例显然不能用一条直线描述全厚度域。

Interactions Table：该图为空白是符合模型定义的，因为本 Study 只有一个可变厚度，无法选出两个不同输入变量 A/B。不要把空白解释为“两个物理效应不存在”；材料并未在同一个 Study 内作为输入变量。

Fit → Post-Processing：先辨认真实响应与预测响应

Fit Integrity：原始响应和带 Fit 图标的预测响应分别列出。它检查数据是否齐全，不评价代理模型是否准确；max_node 的零变化仍会以红色出现。

Fit Summary：同一设计行可比较 load_disp 与 load_disp@fit。两列接近才说明该点被拟合好；还需看样本之外的误差。对本例，不要因为 mass@fit 与 mass 一致就认为位移 Fit 也可靠。

Fit Parallel Coordinate：可把原始响应与 Fit 响应放进同一多轴图，快速发现明显偏离的样本。多轴折线只能定位异常行，定量误差仍看 Summary 与 Residuals。

Fit Scatter：可选 load_disp 作 X Axis、load_disp@fit 作 Y Axis。若预测理想，点应靠近 45° 线；画图前必须确认右侧两个轴选的是原始值和拟合值，而非两个相同变量。

Residuals Table：残差是“求解值－拟合值”。当前默认选中质量，点几乎在 0 线上，因为质量随厚度线性变化。右侧改选 load_disp 才能检查位移误差；观察误差是否在极薄点系统性增大，不能以默认质量残差替代位移验证。

Trade-off Table：拖动输入厚度滑块，下面的 Output Predictions 即时显示拟合预测。它用于交互探索候选值，Verify 才会触发进一步校核；单看预测数值不能替代真实 OptiStruct 求解。全域位移 Fit 的 R² 只有 0.4092176，因此尤其不能用这页极薄区预测值直接判定约束。