Ogden — CAE术语解释
Ogden的理论基础
Ogden的基本概念
经常听到Ogden模型这个词,但它与胡克定律或Neo-Hookean模型的根本区别是什么?
根本区别在于应变能的表示形式。胡克定律是微小变形下的线性弹性,Neo-Hookean是用主伸长比平方和表示的最简单超弹性模型。而Ogden模型直接用主伸长比
当你说"主伸长比的幂次"时,我很难想象。
是的,这种情况确实存在。例如,天然橡胶单轴拉伸试验数据用3项Ogden模型(
材料常数
Ogden的数值计算方法
实施和稳定性
在有限元法中实现Ogden模型时,有什么特别需要注意的数值问题吗?
主要问题是"体积锁定"和"非凸性"。由于Ogden模型通常假设不可压缩,如果不使用混合形式(位移-压力)的单元(如Abaqus的C3D8H或C3D10H),就会发生过度刚性(体积锁定)。另外,当
软件是如何检测和防止因"非凸性"导致的不稳定性的?
主要求解器会检查每个积分点的切线刚度矩阵的正定性。例如,Abaqus/Standard在材料雅可比(切线刚度)不再正定时会发出"负特征值"警告,并自动缩小步长。在实施中,还对主伸长比
推导应变能函数的第二变分,即材料雅可比
你说对了,手工计算非常复杂。在实践中,通常采用自动微分(Automatic Differentiation, AD)或数值微分。例如,COMSOL Multiphysics的"非线性结构材料"模块将用户输入的应变能公式进行符号处理,自动生成必要的一阶和二阶偏导数。而Abaqus等商业代码对Ogden模型有预先优化的闭式雅可比公式,这些公式基于主伸长比和主方向的谱分解,计算效率很高。
Ogden的实际应用
材料常数的确定
实际上想确定橡胶材料的Ogden常数时,需要什么试验数据,如何进行拟合?
至少必须有单轴拉伸试验数据。理想情况下,应有单轴拉伸、等轴拉伸(或压缩)、平面拉伸(纯剪切)三种试验数据,这样能在多轴状态下很好地规定材料行为。拟合通常使用Abaqus/CAE附带的"材料标定"功能、MSC Marc的"材料拟合"工具,或专门的软件如MCalibration(PolymerFEM公司)。这些工具通过非线性优化算法(如Levenberg-Marquardt法),自动计算
项数
好问题。增加项数过多会导致"过度拟合"问题。由于实验数据含有噪声,如果用N=6,就会完美拟合噪声,表观精度提高但在实际多轴状态下的预测精度反而下降,数值不稳定性也增加。实际指导原则是从N=2或3开始。基于ISO 17200等橡胶材料试验规范的数据,通常N=3就足够了。拟合后必须进行"验证"过程:用未使用的变形模式(如仅用平面拉伸数据拟合时的等轴压缩行为)进行模拟预测,看是否合理。
如何判断拟合结果的"质量"?有定量指标吗?
主要有两个指标。一个是决定系数
Ogden的软件比较
各软件中的实施
在Ansys、Abaqus、COMSOL中定义Ogden模型时,输入方法或选项有差异吗?
有很大差异。在Abaqus中,选择"超弹性"材料模型的子选项"Ogden",直接输入项数N(最多6)和各项的
"三阶Ogden模型"这个词在Abaqus和COMSOL中的含义不同?
是的,很容易混淆。在Abaqus和学术文献中,"N项Ogden模型"的N指级数的项数。而COMSOL中的"阶数(Order)"为1、2、3等,它与应变能公式的展开阶数相关,COMSOL的"阶数N"实际对应"2N-1项Ogden模型"。也就是说,COMSOL选择"阶数3"时,内部实际使用5项Ogden模型。因此,从其他软件或论文中使用的常数值输入到COMSOL时需要进行转换。查阅手册中的公式是必须的。
免费或开源的CAE软件(CalculiX、Code_Aster等)也能使用Ogden模型吗?
可以。主要的开源求解器都支持。CalculiX(预处理器是PREPO7)用"*HYPERELASTIC"卡的"OGDEN"类型定义。Code_Aster中在`DEFI_MATERIAU`命令内的`ELAS_HYPER`关键字指定`OGDEN`。但与商业软件不同,这些通常不附带图形化材料拟合工具,需用外部工具(如MCalibration的免费版或Python脚本)先求出常数再手动输入。另外,用户需要深入理解求解器设置才能进行收敛性调优,难度相对较高。
Ogden故障排除
常见错误和对策
用Ogden材料的模拟出现"负雅可比"或"过度变形"错误导致不收敛。首先应该怀疑什么?
首先应怀疑"材料常数的合理性"。特别是在发生多轴变形(包括压缩或剪切)的区域,材料模型可能在物理上变得不稳定(非凸)。对策步骤是:1) 评估常数拟合的变形模式(通常是拉伸)之外的情况。进行简单的单元试验(如悬臂梁弯曲的小变形),检查是否出现异常柔软。2) 确认初始刚度
请解释"体积锁定"的具体症状以及在Ogden模型中如何确认。
症状是模型异常坚硬,不能变形,反力过大。特别是弯曲变形中很明显。由于Ogden模型严格处理不可压缩性,如用标准位移单元(如Abaqus的C3D8),必然发生锁定。确认方法:在Abaqus中,将网格改为"混合"单元(C3D8H、C3D10H)重新计算,看结果是否明显变软。另外,可视化单元体积应变(VOLUMETRIC STRAIN),在锁定情况下应该接近零。如果发生锁定,体积应变会出现不自然的波动或大值。
拟合工具的拟合精度很好,但实际部件模拟的应力与试验不符。除了材料模型还有什么原因?
材料模型以外的主要原因有三个:1) Mullins效应(应力软化):试验片是初次加载数据,但部件模拟是往复动作,存在历史依赖的软化被忽略。Abaqus可将"Mullins效应"与Ogden结合定义。2) 应变率依赖性(粘弹性):橡胶加载速度越快越硬。稳定超弹性模型无法考虑这点。3) 边界条件和摩擦:接触面摩擦系数(如0.5~0.7)是否现实。这些都是单纯试验片拟合无法捕捉的结构模拟特有问题。
想给Ogden常数加上温度依赖性时怎么办?
直接能"定义Ogden模型常数作为温度函数"的软件很少。通常的做法有两种:1) 在各温度(如-20°C、23°C、80°C)分别进行材料试验,为各温度独立拟合Ogden常数集合。然后用软件的"温度依赖表"功能(Abaqus中结合`*ELASTIC, DEPENDENCIES=1`和`*EXPANSION`等),让常数集合随温度插值。2) 物理基础方法:假定Ogden模型的
相关主题
详情
错误