4. 材料卡片配置实战

在HyperMesh中，材料卡片是求解器输入文件里MAT1、MATS1、MATHE、MAT8等关键字段的图形化映射。卡片填错，求解器要么直接报错终止，要么给出看似收敛却完全失真的结果。本节以OptiStruct/Nastran格式为主线，把四类常用材料模型的字段含义、填写顺序和校验方法一次讲清。

一、先分清四类材料模型各自解决什么问题

线弹性材料只描述“应力与应变成正比”的阶段，适合刚度分析、模态分析和大多数金属在屈服前的工况；塑性材料在屈服后引入硬化规律，用于承载能力、残余变形和碰撞类分析；超弹性材料用应变能势函数描述大变形下的非线性弹性行为，典型对象是橡胶、硅胶、密封圈；复合材料则要同时表达各向异性刚度和铺层方向，单层用MAT8，叠层顺序交给PCOMP或PCOMPG。

选择卡片前先问三个问题：材料是否各向同性？是否允许永久变形？是否存在大变形或铺层方向？答案决定了卡片类型，也决定了后续必须填哪些字段。

二、线弹性MAT1：字段少，但单位与量纲最容易出错

MAT1是使用频率最高的卡片。在HyperMesh的Materials面板中新建材料，Card Image选择MAT1，至少填写三个字段：弹性模量E、泊松比NU、密度RHO。若材料为各向同性，只需E和NU；若已知剪切模量G，可单独填写，否则求解器按G=E/[2(1+NU)]自动计算。

以结构钢为例，采用mm-t-N-MPa一致单位制：E=210000 MPa，NU=0.3，RHO=7.85e-9 t/mm³。这里密度必须换算，因为1 t/mm³=1e9 kg/m³，钢材7850 kg/m³对应7.85e-12 t/mm³？实际换算为：7850 kg/m³=7.85e-9 t/mm³。填错量级会让重力载荷和模态结果完全失真。

校验方法：在HyperMesh中打开Card Edit，确认E、NU、RHO均非空；导出为.fem后搜索MAT1，检查字段位置是否与求解器手册一致。对于温度相关材料，还需通过TAB1定义E随温度变化，并在MAT1的TID字段引用。

三、塑性MATS1：屈服与硬化曲线的填写逻辑

塑性分析不能只填MAT1。正确做法是：先建MAT1给出弹性段参数，再建MATS1描述屈服后的行为。MATS1的关键字段包括：TID指向对应的MAT1编号；TYPE选择PLASTIC（理想塑性）或HARD（等向硬化）；YF指定屈服准则，常用VONM（von Mises）；H为硬化斜率，若使用表格曲线则置0；TAB1引用应力-塑性应变数据表。

示例：某低碳钢屈服强度235 MPa，弹性模量210000 MPa，硬化段取塑性应变0、0.02、0.05、0.10对应应力235、300、360、420 MPa。在HyperMesh中通过Curves面板创建TAB1，X列为塑性应变，Y列为真实应力。注意曲线必须单调递增，且第一个点的应力应等于屈服应力，否则求解器会提示屈服面定义不一致。

校验时重点看两处：一是MATS1的TID是否与MAT1的MID一致；二是TAB1的X值是否从0开始。若做的是理想塑性，H填0且不引用TAB1，此时材料一旦屈服即无限流动，仅适用于极限载荷估算。

四、超弹性MATHE：从试验数据到系数拟合

橡胶类材料用MATHE。HyperMesh支持Mooney-Rivlin、Ogden、Neo-Hookean等模型。以Mooney-Rivlin为例，需要填写C10、C01和D1。其中D1控制可压缩性，近似不可压时D1取0或极小值；C10与C01之和的两倍等于初始剪切模量，即2(C10+C01)=G。

假设某橡胶邵氏硬度50A，估算剪切模量G≈0.7 MPa，若取C01=0.25C10，则2(C10+0.25C10)=0.7，解得C10≈0.28 MPa，C01≈0.07 MPa。这组数值属于教学示例，实际工程中应通过单轴拉伸、等双轴拉伸和平面剪切试验数据拟合得到。

校验要点：检查C10+C01是否大于0；若填写D1，确认其值使体积模量远大于剪切模量；导出卡片后确认MATHE的MID与属性卡片引用一致。超弹性分析常伴随大变形，还需在分析步中打开大变形选项（如NLPARM或NLGEOM）。

五、复合材料MAT8与PCOMP：各向异性与铺层方向

复合材料单层用MAT8定义，字段包括E1、E2、NU12、G12、G1Z、G2Z，以及强度参数XT、XC、YT、YC、S。E1为纤维方向模量，E2为垂直纤维方向模量，通常E1远大于E2。NU12为主泊松比，需满足NU12·E2=NU21·E1，否则刚度矩阵不对称。

示例：T300/环氧单层，E1=135000 MPa，E2=9000 MPa，NU12=0.3，G12=5000 MPa，G1Z=G2Z=4500 MPa，单层厚度0.125 mm。强度参数XT=1500 MPa，XC=1200 MPa，YT=50 MPa，YC=200 MPa，S=70 MPa。这些数值为典型教学示例，实际应以材料供应商数据为准。

铺层顺序由PCOMP或PCOMPG定义，通过PLYID引用MAT8，并给出每层厚度和铺层角。校验时检查：所有PLYID是否已定义；铺层角是否与材料主方向一致；总厚度是否等于各层厚度之和；对称铺层是否满足均衡性要求。若使用PCOMPG，还需确认全局铺层编号与区域定义匹配。

六、统一校验流程与常见报错

无论哪类材料，导出求解器文件后都建议做四步检查：第一，搜索卡片关键字，确认字段数量与顺序正确；第二，检查所有引用关系（TID、PLYID、MID）是否闭合；第三，核对单位制是否统一；第四，用一个小模型做单步静力分析，观察位移和应力是否在合理量级。

常见报错中，“negative slope in TAB1”通常指塑性曲线下降；“invalid material id”多为引用缺失；“singularity”可能来自超弹性系数导致刚度为负。遇到这些提示，回到卡片逐字段核对，比反复调整网格更有效。