材料数据的验证和不确定性
理论：拟合的试验与需要预测的变形
概述
材料常数是按试验数据拟合曲线确定的，这样还不够吗？
与拟合试验相同的变形时可能足够。问题在于零件中的变形与试验不同时。橡胶这类材料在拉伸、压缩、剪切、面内双向拉伸时响应不同，某些常数组合对一种变形吻合，对另一种却偏差很大。曲线拟合得好会让人放心，但那只是再现了拟合用的数据，不能证明能够预测。所以要用拟合中未使用的试验确认预测，这就是材料数据验证。
基准材料（Treloar橡胶，三项Ogden式）
P为名义应力，λ为伸长比。常数μ = 0.618、0.00118、−0.00981 MPa，α = 1.3、5.0、−2.0（Ogden，1972年），初始剪切模量0.414 MPa。用此式算出的值视为“试验数据”。
Mooney-Rivlin模型
C01 = 0时即为Neo-Hookean模型。用最小二乘法拟合，单轴伸长比1.05～3.0，等双轴1.05～2.5。
Treloar的数据
1944年，Treloar对硫化天然橡胶进行了单轴拉伸、等双轴拉伸和纯剪切三种试验，测量了直到大伸长的应力。这组数据作为同一材料不同变形模式的代表，至今仍是比较新超弹性模型性能的基准。Ogden在1972年提出的模型就是用这组数据证明了一组常数能很好地表示三种变形模式。
算例
例：拟合方式对预测的影响
| 模型与拟合方式 | C10 (MPa) | C01 (MPa) | 单轴最大误差 | 等双轴最大误差 | 等双轴λ=2.0应力 |
|---|---|---|---|---|---|
| Neo-Hookean（仅单轴） | 0.160 | 0 | 22.2% | 26.6% | 0.629 MPa |
| Mooney-Rivlin（仅单轴） | 0.112 | 0.112 | 6.6% | 309.5% | 2.209 MPa |
| Mooney-Rivlin（单轴＋等双轴） | 0.160 | 0.009 | 17.9% | 20.1% | 0.774 MPa |
基准等双轴应力在λ = 2.0时为0.806 MPa。单轴λ = 3.0时，仅按单轴拟合的Mooney-Rivlin为0.863 MPa，与基准完全一致。
单轴中吻合最好的Mooney-Rivlin，在双轴中偏差最大。
两个常数C10和C01在单轴中作用形式相似，仅凭单轴数据几乎确定不了两者的分配。最小二乘法选择最贴合单轴曲线的分配，结果使C01变大。在等双轴中C01要乘以λ²，作用方式大不相同，C01过大就会把应力估大好几倍。单轴吻合在6.6%以内，双轴却偏差2.7倍，原因就在这里。同时使用双轴数据，C01被压小，两者误差都在两成左右。橡胶密封件、轮胎、膜鼓胀这类以双轴变形为主的零件，使用仅凭单轴数据确定的常数很危险。
材料数据验证步骤
- 根据零件分析结果，确认主要变形模式以及应变、温度、应变率范围。
- 准备覆盖这些范围的试验（单轴、双轴、剪切、体积）。
- 同时拟合多项试验，确认常数物理上合理（稳定性等）。
- 用拟合中未使用的试验或零件试验确认预测。
- 根据试验离散性估算材料常数不确定性，评估其对结果的影响。
“密封面压不够”
某橡胶O形圈在分析中接触面压充足，样品却发生泄漏。材料常数只按材料厂商的单轴拉伸数据拟合。被压缩的O形圈截面中，表面附近的变形接近双轴。增加双轴和平面应变试验重新拟合后，接触面压降低约三成，解释了泄漏。改为增大压缩量的设计后解决。
常见错误
错误与对策
| 错误 | 后果 | 对策 |
|---|---|---|
| 只按单轴拟合 | 其他变形偏差大 | 多模式拟合 |
| 以曲线吻合当作已验证 | 预测能力不明 | 用未用试验确认 |
| 超出试验范围使用 | 外推不稳定 | 确认范围 |
| 忽略温度和速率依赖 | 条件变化即偏差 | 按条件试验 |
| 不考虑离散性 | 过信余量 | 评估不确定性 |
我想了解相关内容。
相关页面参考ASME V&V指南、边界条件检查、模型审查清单、NAFEMS质量体系和量纲分析。
价值
更详细
错误