08 ANSYS13.0 Workbench 结构非线性培训 超弹性

6-17
Workbench Mechanical - Hyperelasticity
... 混合u-P 公式
Training Manual
• „Mixed u-P‟公式满足体积协调约束. Mixed u-P 代码在位移自由度(u)外自 动添加一内部静水压力自由度(P) 来强迫体积协调条件. 因此, 名为 Mixed u-P.
Training Manual
• 超弹性本构模型通过应变能密度函数来定义
– 与塑性不同, 超弹性不定义为速率公式 – 相反, 总应力与总应变的关系由应变势能 (W)来定义
D:ε σ
6-5
Workbench Mechanical - Hyperelasticity
... 延伸率定义
为理解应变能密度函数的形式，一些定义是必要的:
Unextended chain
Training Manual
Extended chain
– 橡胶强化通过橡胶硫化过程中分子链间形成交联 来实现(例如，可能的链结构数量减少).
• 这不象金属，例如，弹性行为是由于分子连接展开。 (见第5章)
Schematic of single molecular chain. In network, these chains are randomly oriented and often have crosslinks.
J total J th
6-9
Workbench Mechanical - Hyperelasticity
... 应变势能的定义
• 应变势能 (或应变能函数)通常表示为W
– 应变势能或者是主延伸率的直接函数，或者是应变不变量的函数
Training Manual
W W I1 , I 2 , I 3 or W W l1 , l2 , l3
N
o
a
i 1 i
i
2
o
2 d1
6-19
Workbench Mechanical - Hyperelasticity
... 体积容差
• 体积协调约束中的容差(vtol)可通过 Command Objects放松.
Training Manual
为接受后续的solc,,,vtol手动激活 Mixed u-P 是必要的
Training Manual
若材料是完全不可压缩的, 则I3 = 1.
• 由于假设材料为各向同性，应变势能的一些形式可以表示为这些标量不变量的函数。 换句话说, 应变不变量是应变的度量，与用于度量应变时使用的坐标系无关。
6-8
Workbench Mechanical - Hyperelasticity
Training Manual
– Vtol 默认值为1e-5.
• 如由于Mixed u-P 不能满足体积协调条件而导致收敛失败, 通过定义 ‘di‟ 参 数的一个小的非零值引入一不可压缩小量是有帮助的.
– di 反比于体积模量. 默认地, 如曲线拟合(下一部分)中没引入体积试验数据, 则材料 假定为完全不可压缩的 (di=0).
• 从前面幻灯片中知, 对完全不可压缩材料体积比 J 应为常数 (J=1), :
J 1 0
• 这导致下面的体积协调方程:

V
J 1 dV J Vtol V
6-18
Workbench Mechanical - Hyperelasticity
... 混合 u-P 公式
• 当不满足条件的时候，求解信息会记录.
第一应变不变量现象模型
Neo-Hookean
第一应变不变量现象模型
基于应变不变量的近似/完全不可压缩现象超弹性模型.
.
6-15
Workbench Mechanical - Hyperelasticity
... W的特殊形式
• 根据所选择的应变能函数W, 材料常数数目不同.
• 选择依赖于:
• 高弹体类型 • 加载条件 (% Strain)
Training Manual
• 可利用的材料试验数据.
• 总体而言，应变能密度函数是最接近应力应变实验数据曲线的函数。
– 在大多多挑战性的应用里，往往只有一个函数是合适的。 – 通常两个或更多函数产生同样的响应。
6-16
Workbench Mechanical - Hyperelasticity
第六章
超弹性
Workbench – Mechanical 结构非线性
6-1
Workbench Mechanical - Hyperelasticity
超弹性章节综述
• 本章包含以下内容:
– – – – 高弹体背景 超弹性理论 曲线拟合过程 查看结果
Training Manual
• 本章节描述的能力适用于ANSYS Structural及以上的 licenses.
F
Tension
u
Compression
6-4
Workbench Mechanical - Hyperelasticity
B.超弹性背景
• ANSYS中关于超弹性本构模型有一些关键假设
– 材料响应是各向同性、等温和弹性的 • 热膨胀是各向同性 • 变形完全可恢复(保守的) – 材料是完全或几乎不可压缩的 – 更复杂的真实橡胶行为理想化
– 例外情况会根据需要说明
6-2
Workbench Mechanical - Hyperelasticity
A. 高弹体背景
• 高弹体是一种聚合物, 具有如下性能
– 高弹体包括天然和合成橡胶, 它是非晶态的, 由 长的分子链组成
• 分子链高度扭转、卷曲, 且在未变形状态下取向任 意 • 在拉伸载荷作用下, 这些分子链部分变得平直、不 扭曲 • 去除载荷后， 这些分子链恢复最初的形态
1 2k W cij I1 3 I 2 3 J 1 i j 1 k 1 d k
i j
N



N
• 其中初始体积模量和初始剪切模量是
o 2c10 c01
2 o d1
• cij 和 di 通常定义为材料性质.
– 如果未知，这些值可以通过从试验数据用曲线拟合的方法得到（下节讨论）.
Training Manual
• 延伸率(或只是‘延伸’) 定义为
L Lo u l 1 e E Lo Lo
• 上述为一个定义橡胶试样单轴拉伸的延伸率的例子, 式中eE 为工程应 变。有三个主延伸率l1, l2, 和 l3, 它们用来度量变形，也用于定义应 变势能。
6-6
Workbench Mechanical - Hyperelasticity
... 混合u-P 公式
• 超弹性体的体积比 (J) 定义为:
Training Manual
V J Vo
– 其中 V和Vo 分别是单元的更新的和初始体积.
• 为保持不可压缩行为, 必须满足体积协调性约束.
– 完全-不可压缩超弹性材料, 不会发生体积改变. – 使用 J, 可确定体积变化数 – 完全-不可压缩状况, J 应等于 1. 换句话说, 最终和初始体积是相同的 (无体积变 化)
6-11
Workbench Mechanical - Hyperelasticity
... 应力和应变的计算
• 通过应变能密度函数，可计算应力应变.
– 需要采用应力-应变共轭测量 – 基于W,确定第二 Piola-Kirchoff 应力(和 Green-Lagrange 应变) :
Training Manual
– 式中初始体积和剪切模量被定义为
o
a
i 1 i
N
i
2
2 o d1
i, ai,和 di是用户定义的材料属性. – 如果未知，可以从试验数据经曲线拟合推导而来.
6-14
Workbench Mechanical - Hyperelasticity
... 模型分类
Training Manual
Example of Rubber boot, o-rings/seals
6-3
Workbench Mechanical - Hyperelasticity
...高弹体背景
• 宏观上，橡胶行为呈现下列特征
Training Manual
– 能承受大弹性(可恢复)变形，任何地方都可达100-700%。正如前面提及的，这 是由于交联分子链拆开的原因. – 由于分子链的拉直引起变形, 所以在外加应力作用下, 体积变化很小。 因此, 高 弹体几乎不可压缩. – 它们的应力-应变关系是高度非线性的 – 通常, 拉伸状态下, 材料先软化再硬化，而压缩时材料急剧硬化.
• 其他模型 (Mooney-Rivlin, Yeoh 和 Neo-Hookean) 都是多项式模 型的缩减形式。
Polynomial
第一和第二应变不变量现象模型
Mooney-Rivlin
第一和第二应变不变量线性模型
Yeoh
第一应变不变量现象模型
2-term M-R
第一第二应变不变量现象模型
3-term Yeoh
• 求解信息将记录这个改变
6-20
Workbench Mechanical - Hyperelasticity
... 混合 u-P考虑事项
Training Manual
• 对完全不可压缩问题, 如果所有的边界节点已指定位移，不存在唯一的求解.
– 这是由于静水压力 (内部DOF) 是独立于变形的. 静水压力需由力/压力边界条件 来确定. 否则不能计算静水压力，没有唯一求解. – 对这些情形, 至少一个节点没有应用边界条件将矫正这种情形.