Neo-Hookean: $\sigma = 2C_{10}(\lambda^2 - \lambda^{-1})$
Mooney-Rivlin: $\sigma = 2(\lambda^2 - \lambda^{-1})(C_{10}+ C_{01}/\lambda)$
Ogden: $\sigma = \mu_1(\lambda^{\alpha_1 - 1}- \lambda^{-\alpha_1/2 - 1})$
剪切模量: NH: μ=2C₁₀, MR: μ=2(C₁₀+C₀₁)
实时对比Neo-Hookean、Mooney-Rivlin、Ogden三种本构模型,计算单轴、等双轴、纯剪切三种变形模式下的应力-伸长率曲线和应变能密度。
Neo-Hookean: $\sigma = 2C_{10}(\lambda^2 - \lambda^{-1})$
Mooney-Rivlin: $\sigma = 2(\lambda^2 - \lambda^{-1})(C_{10}+ C_{01}/\lambda)$
Ogden: $\sigma = \mu_1(\lambda^{\alpha_1 - 1}- \lambda^{-\alpha_1/2 - 1})$
剪切模量: NH: μ=2C₁₀, MR: μ=2(C₁₀+C₀₁)
所有超弹性模型的核心都是定义一个应变能密度函数$W$,应力则由$W$对变形的导数得到。对于单轴拉伸这种最常见的变形,工程应力$\sigma$的计算公式如下:
$$ \sigma = \frac{\partial W}{\partial \lambda}$$其中,$\lambda$是伸长比(变形后长度/原始长度),$W$是应变能密度,不同的模型有不同的$W$表达式。
本模拟器中对比的三种模型,其单轴拉伸应力公式具体为:
$$ \begin{aligned}&\text{Neo-Hookean:}& \sigma &= 2C_{10}(\lambda^2 - \lambda^{-1}) \\ &\text{Mooney-Rivlin:}& \sigma &= 2(\lambda^2 - \lambda^{-1})(C_{10}+ C_{01}/\lambda) \\ &\text{Ogden (N=1):}& \sigma &= \mu_1(\lambda^{\alpha_1 - 1}- \lambda^{-\alpha_1/2 - 1}) \end{aligned}$$变量定义:
$C_{10}, C_{01}$:材料常数,与剪切模量相关。
$\mu_1, \alpha_1$:Ogden模型的材料参数,$\mu$的量纲是应力,$\alpha$是无量纲的指数。
这些参数的值决定了曲线的“陡峭”程度和形状,需要通过实验数据拟合确定。
汽车工业 - 轮胎与密封件:轮胎在滚动时,橡胶材料经历复杂的多轴大变形。使用Mooney-Rivlin或Ogden模型进行有限元分析,可以精确预测其接地压力分布、生热和磨损,从而优化轮胎结构和配方。
消费电子 - 防水密封圈:手机、手表的防水橡胶圈在装配和使用中会被压缩和拉伸。通过超弹性模型模拟其接触压力和回弹力,是确保长期有效密封的关键设计步骤。
生物医学工程 - 人工软组织:人工心脏瓣膜、血管支架的涂层等生物相容性材料,其力学行为与生物软组织类似。Ogden模型因其在高变形下的准确性,常被用于模拟这些植入体在体内的力学响应。
减震与隔振产品:机械设备下的橡胶减震垫、建筑隔震支座等,需要吸收能量并控制振动。准确的超弹性模型能帮助工程师预测产品在不同载荷下的刚度变化和疲劳寿命。
首先,要明确“材料常数并非材料的产品目录值”。例如,当听说某种橡胶的C10=0.5 MPa时,人们容易认为这是拉伸强度或硬度的直接指标。然而,这些常数仅仅是针对“应变能密度函数”这一特定数学模型进行拟合的结果。即使是同一种材料,采用Neo-Hookean模型拟合得到的C10与采用Mooney-Rivlin模型拟合得到的C10在物理意义上也是不同的值。因此,直接使用文献或他公司的数据常数是危险的,必须基于自身的测试数据重新进行拟合。
其次,要注意“仅凭单轴拉伸数据无法决定一切”这个陷阱。单轴试验确实简单,但仅凭此数据常常会低估双轴或剪切行为的响应。例如,仅从单轴数据确定Mooney-Rivlin模型的C01时,可能导致实际部件承受双轴变形时的应力被低估50%以上。在实际工作中,理想的做法是尽可能获取压缩或平面拉伸(接近双轴)的试验数据,并调整参数使曲线在多种变形模式下都能吻合。在本模拟器中切换变形模式,可以直观地看到仅凭单轴数据拟合的曲线在其他模式下会产生多大偏差。
最后,“Ogden模型的阶数N并非越大越好”。Ogden模型确实可以通过增加项数(N=1,2,3...)来拟合复杂曲线,但存在过拟合的风险。当N≥3时,虽然可以完美贴合试验数据的误差范围,但在数据点之间(插值)或之外(外推)可能出现脱离实际的行为。建议从N=1或2开始,根据拟合情况判断是否需要增加阶数。在CAE分析中,计算成本也是不可忽视的因素。
以天然橡胶为例,设Neo-Hookean参数C10=0.35 MPa,在单轴拉伸工况下λ=2.0时,第一不变量I1=λ²+2/λ=4+1=5,Cauchy应力σ=2C10(λ-1/λ²)=2×0.35×(2-0.25)=1.225 MPa;而Mooney-Rivlin模型(C10=0.25 MPa,C01=0.10 MPa)同条件下σ=2(C10+C01/λ)(λ-1/λ²)=1.575 MPa,低应变能密度约10%,更精准拟合实验曲线