4. 线弹性各向同性材料参数标定

线弹性各向同性材料在 Ansys Mechanical 中只需要两个独立弹性常数:弹性模量 E 与泊松比 ν。它们共同决定刚度矩阵,进而决定结构在载荷下的位移、应变与应力响应。所谓“标定”,并不是把试验曲线硬套进软件,而是从单轴拉伸试验的弹性段中提取 E 和 ν,再通过参数化仿真与试验曲线对比,确认这两个参数在目标工况下能复现真实响应。本节以标准哑铃形试样为例,说明从原始数据到参数确认的完整过程。

一、从单轴拉伸数据中提取 E 与 ν

单轴拉伸试验中,试样标距段受轴向拉力 F,同时测量标距 L0 的伸长量 ΔL 和横向宽度 b0 的收缩量 Δb。工程应力 σ = F/A0,工程应变 ε轴向 = ΔL/L0,横向应变 ε横向 = Δb/b0。在弹性段内,弹性模量由轴向应力—应变曲线的斜率给出:

E = Δσ / Δε_轴向
ν = - ε_横向 / ε_轴向

实际操作中,不建议直接取“原点—断裂点”的割线斜率,而应在弹性段内选取一段线性良好的区间做最小二乘拟合。例如某次示例试验中,标距 L0=50 mm、宽度 b0=10 mm、厚度 t0=4 mm,在 2 kN 到 6 kN 之间取 8 个数据点,拟合得到斜率约 68.5 GPa;同时横向引伸计给出 ε横向/ε轴向 ≈ -0.31,则 ν 初值取 0.31。这里的数据属于教学示例,用于说明处理流程,并非实测结果。

需要注意两个细节:第一,若试验机柔度较大,位移测量会包含夹持段和机架变形,应优先使用引伸计或非接触式光学测量;第二,ν 对横向测量精度非常敏感,横向应变片或光学测量的分辨率不足时,ν 的波动可能达到 ±0.03 以上,因此应把 ν 视为一个需要标定的参数,而不是一次读数就固定下来。

二、在 Ansys Mechanical 中定义弹性参数

在 Engineering Data 中新建材料,展开 Linear Elastic → Isotropic Elasticity,输入 Young's Modulus 和 Poisson's Ratio。若后续要做参数化标定,应把这两个值设为参数:在数值输入框右侧点击参数化图标,使其变为 “D” 或 “P” 状态,这样它们会出现在 Parameter Set 中,可被 Design of Experiments 或直接优化驱动。

几何方面,建议按试验实际尺寸建立标距段模型,而不是只建一个等截面长方体。对于哑铃形试样,可只取标距段并施加对称边界,以减少计算量;但若研究夹持效应,则必须把夹持段和过渡圆弧一并建入。材料指派后,网格在标距段采用六面体主导,厚度方向至少两层单元,过渡圆角处适当加密。

三、参数化建模与曲线对比

参数化标定的核心是让仿真输出与试验输出在同一坐标下可比。常用做法是:在试样一端施加固定约束,另一端施加与试验一致的位移或力;提取标距段的轴向应变与横向应变,或者直接提取力—位移曲线。然后在 Parameter Set 中把 E、ν 设为输入参数,把“仿真曲线与试验曲线的偏差”设为输出参数。

曲线对比时,不要只看一条曲线是否“像”。建议同时观察三个量:弹性段斜率、横向收缩比例、卸载后的残余位移。若仿真斜率偏大,说明 E 偏高;若横向收缩明显小于试验,说明 ν 偏小。一个实用的迭代策略是:先固定 ν=0.3,用斜率比调整 E;再固定 E,用横向应变比调整 ν;通常两到三轮即可收敛。若使用 Ansys 的 Direct Optimization,可把目标设为最小化斜率误差与横向应变误差的加权和。

单轴拉伸参数标定:试验数据 → 参数化仿真 → 曲线对比 试验数据 L₀ = 50 mm,b₀ = 10 mm t₀ = 4 mm,标距段 F–ΔL 曲线(弹性段) 横向应变 ε_横向 拟合斜率 ≈ 68.5 GPa ν 初值 ≈ 0.31 (示例数据,非实测) Ansys Mechanical 参数化 Engineering Data: E = P1,ν = P2 几何:标距段 + 夹持段 边界:一端固定,一端位移 输出:斜率、横向应变比 Parameter Set 驱动 目标:复现试验弹性响应 曲线对比与判断 斜率偏大 → E 偏高 横向收缩偏小 → ν 偏小 残余位移偏大 → 约束/夹持 迭代 2~3 轮收敛 确认 E、ν 标定值 判断依据:斜率误差 + 横向比 参数更新后重新计算 试样尺寸与夹持效应示意 夹持段 标距段 L₀ = 50 mm 夹持段 夹持约束 夹持约束 过渡圆角处存在应力集中,夹持段刚度会“吃掉”部分位移

四、试样尺寸与夹持效应的影响

试样尺寸对 E 的提取影响相对小,因为 E 是应力—应变斜率,与绝对尺寸无关;但它对 ν 和局部应力状态影响明显。标距段越短,夹持端约束对应变分布的影响越容易侵入测量区,导致轴向应变被低估、横向收缩被约束,从而把 ν 标定得偏小。一般建议标距段长度不小于宽度或厚度的 5 倍,过渡圆角半径不小于厚度的一半。

夹持效应在仿真中表现为:若把夹持段简化为刚性固定,试样端部无法自由收缩,横向应变会偏小;若把夹持段建为弹性体并施加摩擦接触,则更接近真实,但接触刚度、摩擦系数又会引入新的不确定性。对于以标定 E、ν 为目标的模型,推荐只取标距段,并在两端施加“轴向位移 + 横向自由”的边界条件,避免夹持段干扰。若必须保留夹持段,应把夹持长度设为参数,检查夹持长度变化 10% 时 E、ν 的标定结果是否稳定;若变化超过 2%,说明夹持效应已经不可忽略。

标定顺序建议:先由弹性段斜率确定 E,再由横向应变比确定 ν,最后用完整加载—卸载曲线检查残余位移。不要同时调整两个参数,否则很难判断误差来源。
若试验曲线弹性段本身不直,或试样在夹持处过早屈服,应先检查试验有效性,而不是强行用仿真去“拟合”一条不可靠的曲线。参数标定的前提是试验数据本身可信。

五、可直接使用的标定记录表

下面给出一份简化的标定记录表,可用于记录每轮参数、曲线误差和判断结论。表中“斜率误差”定义为(仿真斜率 − 试验斜率)/试验斜率,“横向比误差”定义为(仿真横向比 − 试验横向比)/试验横向比。

轮次 E / GPa ν 斜率误差 横向比误差 判断
初值 68.5 0.31 +4.2% −6.8% E 偏高,ν 偏小
第 1 轮 65.8 0.33 +0.6% −1.1% 接近,继续微调
第 2 轮 65.4 0.335 −0.2% +0.3% 收敛,确认