4. 材料与截面参数化

在Abaqus参数化建模中，几何参数化只解决了“形状”问题，而真正决定仿真物理行为的是材料本构、截面属性以及它们随温度、场变量的变化规律。本节围绕四个相互关联的对象展开：材料本构参数、截面属性、复合材料铺层、温度与场变量依赖。核心思路是把这些原本写在CAE界面里的数值，转化为可由Python脚本统一管理的参数字典，使材料与截面成为可批量替换、可扫描、可优化的输入变量。

一、材料本构模型的参数化定义

材料本构的参数化，本质是把mdb.models['Model-1'].Material(...)及其下属的Elastic、Plastic、Hyperelastic等子对象，用变量替代硬编码数值。以最常见的金属弹塑性模型为例，线弹性需要杨氏模量E和泊松比ν，塑性需要屈服应力与塑性应变的硬化曲线。参数化时建议按“材料名—本构类型—参数集”三层组织：

mat_params = {
    "Steel_A": {
        "elastic": {"E": 210000.0, "nu": 0.30},
        "plastic": [(320.0, 0.0), (400.0, 0.05), (480.0, 0.15)]
    },
    "Alu_6061": {
        "elastic": {"E": 69000.0, "nu": 0.33},
        "plastic": [(276.0, 0.0), (310.0, 0.04), (350.0, 0.10)]
    }
}

其中plastic列表的每一项为(屈服应力, 塑性应变)，第一点塑性应变必须为0，这是Abaqus对Plastic表数据的硬性要求。脚本中通过循环写入：

for name, p in mat_params.items():
    m = mdb.models['Model-1'].Material(name=name)
    m.Elastic(table=((p["elastic"]["E"], p["elastic"]["nu"]),))
    m.Plastic(table=p["plastic"])

这样只需替换字典内容，即可在不改动建模逻辑的前提下切换材料。对于超弹性橡胶，参数化为C10、C01或Mooney-Rivlin系数；对于损伤模型，则把损伤起始应变、演化位移等一并纳入字典。关键原则是：参数字典的键名与Abaqus关键字保持可追溯的对应关系，避免出现“脚本里叫E1、界面里叫Young”的混乱。

二、截面属性赋值

截面是材料与几何之间的桥梁。参数化截面需要同时管理三件事：截面类型（Solid/Homogeneous、Shell、Beam、Membrane等）、材料指派、以及截面自身的几何参数（如壳厚、梁截面尺寸）。典型写法为：

sec_params = {
    "Sec_Shell": {"type": "HomogeneousShell",
                  "material": "Steel_A", "thickness": 2.0},
    "Sec_Solid": {"type": "HomogeneousSolid",
                  "material": "Alu_6061"}
}
for sname, s in sec_params.items():
    if s["type"] == "HomogeneousShell":
        mdb.models['Model-1'].HomogeneousShellSection(
            name=sname, material=s["material"],
            thickness=s["thickness"])
    else:
        mdb.models['Model-1'].HomogeneousSolidSection(
            name=sname, material=s["material"])

截面赋值到几何区域时，需要区分Set与Surface。实体单元用sectionAssignment(region=cellSet, sectionName=...)，壳单元则赋给面集合。参数化时建议把“区域集合名”也作为变量，使截面与区域的绑定关系可配置，而不是写死在脚本里。若同一部件不同区域需要不同厚度，可建立“区域—截面”映射表，循环完成赋值。

三、复合材料铺层参数化

复合材料是参数化收益最明显的场景。铺层角度、厚度、铺层顺序往往需要大量扫描。Abaqus中通过CompositeLayup定义铺层，核心参数包括：单层材料、单层厚度、铺层方向角、积分点数。参数化结构可设计为：

layup = {
    "material": "CFRP_Ply",
    "plies": [
        {"angle": 0.0,  "thickness": 0.125},
        {"angle": 45.0, "thickness": 0.125},
        {"angle": -45.0,"thickness": 0.125},
        {"angle": 90.0, "thickness": 0.125}
    ]
}

脚本中先创建单层材料（通常为Engineering Constants或Lamina类型，需给出E1、E2、ν12、G12、G13、G23），再按plies列表逐层写入CompositeLayup的plies参数。铺层角度以参考方向为基准，正负号遵循右手定则，全文与图示必须统一。若要做铺层顺序优化，可把plies列表作为设计变量，用循环生成不同排列组合，批量提交计算。

四、温度与场变量依赖

材料参数随温度或场变量变化时，Abaqus允许在Elastic、Plastic等表中增加温度列或场变量列。参数化时把“参数—温度”关系写成嵌套结构：

temp_dep = {
    "Steel_A": {
        "elastic": [
            (210000.0, 0.30, 20.0),
            (200000.0, 0.30, 200.0),
            (185000.0, 0.30, 400.0)
        ]
    }
}

表中每行末尾为温度值，Abaqus按线性插值获取中间温度下的参数。若同时依赖场变量，则温度列之后再加场变量列，顺序不能颠倒。参数化脚本只需把该列表直接传给table参数即可。需要注意的是，温度依赖与场变量依赖的列顺序、单位必须与模型整体单位制一致，否则插值结果会严重偏离。

五、参数化实施中的关键判断

把上述四类参数整合进一个脚本时，需要关注几个容易出错的点。第一，材料名与截面名的唯一性：批量生成时若名称重复，Abaqus会报错或覆盖，建议用“材料_温度_序号”的命名规则。第二，单位制一致性：E用MPa、厚度用mm、温度用℃时，热膨胀系数必须与之匹配，否则温度载荷结果失真。第三，铺层与截面的绑定：复合材料铺层通常通过CompositeLayup直接生成壳截面，再赋给壳区域，不能再叠加普通壳截面。第四，场变量依赖的插值：Abaqus默认线性插值，若数据点过少，中间温度下参数可能偏离实际，建议在关键温度区间加密数据点。

完成本节参数化后，材料与截面不再是CAE界面里手工填写的固定值，而是脚本中可替换、可循环、可优化的数据对象。下一阶段即可把这些字典作为设计变量，接入参数扫描或优化算法，实现从“能建模”到“能自动寻优”的跨越。