📚 课程目录第3章 OptiStruct频响分析卡片体系

3. OptiStruct频响分析卡片体系

频响分析（Frequency Response Analysis）在 OptiStruct 中并不是靠一个“求解类型”开关就能跑起来的，它依赖一组分工明确的 Bulk Data 卡片协同工作：有的负责定义激励，有的负责把激励映射到有限元模型上，有的负责描述频率轴，有的负责指定求解方法与阻尼。理解这些卡片各自“管什么”、彼此“怎么连”，是搭建一个可运行频响 deck 的前提。本节围绕 RLOAD、DAREA、TABLED、FREQ、METHOD、SDAMPING 这几张核心卡片展开，并给出一个可直接改写使用的 deck 骨架。

一、频响分析的卡片分工

可以把频响 deck 的卡片体系理解为三条并行的“信息链”，它们在求解时被 OptiStruct 汇合：

激励链：RLOAD1/RLOAD2 定义载荷的幅值与相位（或实部/虚部），DAREA 定义载荷作用的自由度与缩放系数，TABLED1/TABLED2 提供随频率变化的载荷幅值曲线。三者组合后，才形成“在哪个节点、哪个自由度、以什么频率规律施加多大载荷”。
频率链：FREQ 或 FREQ1 定义求解的频率点集合，METHOD 通过引用 FREQ 卡片并指定求解方法（如模态法 MODAL 或直接法 DIRECT），把频率轴与求解策略绑定。
阻尼链：SDAMPING 定义模态阻尼（通常以 TABDMP1 表格给出），在模态法频响中决定各阶模态的阻尼比；直接法频响则更多依赖 PARAM、G 或 CDAMP 等定义。

这三条链缺一不可：没有激励链，模型没有输入；没有频率链，求解器不知道在哪些频率上计算；没有阻尼链，共振峰会被严重高估，结果失去工程意义。

二、关键卡片的作用与关联

RLOAD1 / RLOAD2：RLOAD1 用实部与虚部（或幅值与相位）描述载荷，RLOAD2 则通过 TLOAD 与相位表组合。两者都通过 EXCITEID 或 EXCITE 字段引用 DAREA 或直接引用节点。RLOAD2 的典型写法是：

RLOAD2  10      1       1.0     0.0     20      30
+       0.0     0.0

其中 10 是载荷集 ID，1 是 DAREA 集 ID，20 是幅值表 TABLED1 的 ID，30 是相位表 ID。

DAREA：把载荷“落到”具体自由度上。字段依次为节点 ID、方向（1~6 对应 X/Y/Z/RX/RY/RZ）、缩放系数。例如：

DAREA   1       1001    3       1.0

表示在节点 1001 的 Z 向施加系数为 1.0 的单位载荷。DAREA 本身不随频率变化，频率变化由 TABLED1 承担。

TABLED1 / TABLED2：以 (频率, 幅值) 对的形式给出载荷随频率的变化。TABLED1 是线性插值，TABLED2 是对数插值。例如：

TABLED1 20
+       0.0     1.0     50.0    1.0     100.0   2.0
+       200.0   2.0     ENDT

表示 0~50 Hz 幅值为 1.0，100~200 Hz 幅值为 2.0，中间线性过渡。

FREQ / FREQ1：FREQ 直接列出频率点，FREQ1 用起始频率、步长、步数生成。例如：

FREQ    5       0.0     10.0    20.0    50.0    100.0
+       200.0

表示在 0、10、20、50、100、200 Hz 六个频率点求解。

METHOD：把 FREQ 与求解方法绑定。例如：

METHOD  1       MODAL   1       1.0

其中 1 是 METHOD 集 ID，MODAL 表示模态法，后面的 1 引用 FREQ 集 ID。若用直接法，则写 DIRECT。

SDAMPING：定义模态阻尼。通常引用 TABDMP1 表格：

SDAMPING 1       1
TABDMP1 1
+       0.0     0.02    100.0   0.02    200.0   0.03
+       ENDT

表示 0~100 Hz 阻尼比 2%，200 Hz 阻尼比 3%。SDAMPING 只在模态法频响中生效，直接法需用 PARAM,G 或 CDAMP 定义。

三、卡片关联关系图

四、搭建可运行的频响 deck 骨架

下面给出一个模态法频响分析的最小可运行 deck 骨架。假设模型已包含网格、材料、属性与约束，节点 1001 的 Z 向施加单位载荷，频率范围 0~200 Hz，阻尼比 2%。

SUBCASE 1
  LABEL = Frequency Response
  METHOD = 1
  DLOAD  = 10
  SDAMPING = 1
  FREQ   = 1
  DISPLACEMENT = ALL

BEGIN BULK
DAREA   1       1001    3       1.0
RLOAD2  10      1       1.0     0.0     20      0
TABLED1 20
+       0.0     1.0     200.0   1.0     ENDT
FREQ    1       0.0     10.0    20.0    50.0
+       100.0   200.0
METHOD  1       MODAL   1
SDAMPING 1      1
TABDMP1 1
+       0.0     0.02    200.0   0.02    ENDT
ENDDATA

这个 deck 中，SUBCASE 通过 METHOD、DLOAD、SDAMPING、FREQ 四个字段把三条链“接”到一起。DLOAD 引用 RLOAD2 的载荷集 ID 10，RLOAD2 再引用 DAREA 集 1 和 TABLED1 集 20。METHOD 引用 FREQ 集 1，SDAMPING 引用 TABDMP1 集 1。任何一处 ID 写错，求解都会报“未找到引用卡片”或直接忽略激励。

五、常见错误与检查要点

DAREA 方向写错：方向 1~6 对应 X/Y/Z/RX/RY/RZ，写错会导致激励加在错误自由度上，频响曲线形状完全改变。
TABLED1 缺少 ENDT：表格未正常结束，求解器可能读入后续卡片作为表格数据，导致频率轴错乱。
METHOD 与 FREQ 的 ID 不一致：METHOD 第二个字段必须与 FREQ 集 ID 相同，否则频率轴为空。
SDAMPING 用于直接法：直接法频响不读 SDAMPING，需改用 PARAM,G 或 CDAMP 定义阻尼，否则共振峰无衰减。
RLOAD2 的相位表 ID 写 0：表示相位为零，若需要相位差激励，必须提供 TABLED1 相位表。

建议在提交计算前，先用 OptiStruct 的 CHECK 或 PARAM,CHECKOUT,YES 做一次卡片引用检查，确认所有 ID 都能被正确解析，再运行完整频响求解。

掌握这套卡片体系后，后续无论是做频响优化、还是把频响结果作为约束加入结构优化，都能在 deck 层面快速定位问题：先看激励链是否完整，再看频率链是否连续，最后确认阻尼链是否与求解方法匹配。这三步检查，基本能覆盖频响 deck 搭建中绝大多数“跑不起来”或“结果不对”的情况。