ansys高级非线性分析七粘弹性.ppt

1、粘弹性,七,September 30, 2001 Inventory #001491 7-2,粘弹性本章综述,本章讨论 ANSYS 用于模拟玻璃(非晶态固体)和非晶态聚合物等材料的粘弹性能力。 因为粘弹性行为非常复杂, 本章将主要按如下步骤编排: 首先定义一些常用的术语 通过使用一维流变模型来解释线性粘弹性行为，这将有助于说明广义 Maxwell 模型的基本特征。 使用所有粘弹性材料通用的一般输入要求。 焦点将转至 WLF 偏移函数和聚合物。 将他讨论 TN偏移函数及其对玻璃材料的适用性。,September 30, 2001 Inventory #001491 7-3,粘弹性本章综述,本章包

2、括以下主题: A. 粘弹性理论背景 B. 流变模型 (Maxwell, Kelvin-Voigt, 标准线性) C. ANSYS 粘弹性模型输入 D. WLF 偏移函数 E. TN 偏移函数 F. 求解粘弹性模型 G. 实验数据的曲线拟合,September 30, 2001 Inventory #001491 7-4,粘弹性A. 粘弹性理论背景,有些非晶态聚合物的行为随温度而改变。 在玻璃转变温度以下, 材料行为与弹性固体类似。 在玻璃转变温度以上, 材料响应与一个橡胶固体类似。 在高温时, 材料行为与粘性液体类似。 在玻璃转变温度以上, 响应是弹性固体和粘性液体的结合(指上面的橡胶固体)，

3、这种行为是粘弹性 的特性。,September 30, 2001 Inventory #001491 7-5,粘弹性 . 粘弹性理论背景,粘弹性是率相关行为, 材料特性可能与时间和温度都有关，粘弹性响应可看作由弹性和粘性部分组成。 弹性部分是可恢复的, 且是瞬时的。 粘性部分是不可恢复的, 且在整个时间范围内发生。 ANSYS 中能模拟线性粘弹性，这导致如下假设: 应变率与瞬态应力成比例 瞬态应变与瞬态应力也成比例 限于小应变、小变形行为(NLGEOM,OFF),September 30, 2001 Inventory #001491 7-6,粘弹性 . 术语的定义,速率效应 对于标准的线性模

4、型(后面讨论), 极限行为是: 非常慢和非常快的应变率的弹性。,September 30, 2001 Inventory #001491 7-7,粘弹性 . 术语的定义,蠕变 在恒定的外加应力作用下, 应变单调增加。 右图所示, 线性和指数蠕变情况 应力松弛 在恒定外加应变作用下, 应力渐近降低。,September 30, 2001 Inventory #001491 7-8,粘弹性 . 术语的定义,简单热流变行为 如前所述, 粘弹性材料与时间和温度有关，两个相关性都需要考虑。 简单热流变 (TRS)行为是指时间和温度是同一现象，这意味着粘弹性响应-对数时间的函数关系随着温度变化而平移。 上

5、面说法的另一种解释是, 材料对高温载荷的短时间作用的响应与较低温度长时间作用的响应是相同的。,September 30, 2001 Inventory #001491 7-9,粘弹性 . 术语的定义,简单热流变行为(续) ANSYS的粘弹性能力采用 TRS 行为，该假设考虑了与时间和温度相关的关系，这足以描述很多非晶态聚合物。 这一假设的结果是, 在后面将看到, 这意味着短期G(t=0) 和长期G(t=) 模量保持一致, 与温度无关(即对于任何平移, 前面曲线图的上下限保持一致)。 这样可以在某一温度下定义粘弹性行为而捕捉其它温度下的响应。采用缩减的时间 和偏移函数(下面讨论)的概念, 粘弹性

6、响应曲线被移动来说明另一个温度下的行为。根据材料的不同而采用不同的偏移函数。,September 30, 2001 Inventory #001491 7-10,粘弹性 . 术语的定义,缩减的或伪时间(虚拟温度) 应该注意前述幻灯片的响应曲线的移动是由一个被称作为缩减(或伪)时间 的变量(x) 来完成的，用缩减时间后, 等温方程可用于描述非等温过程。一个虽然独立但却相关的概念是虚拟温度，后面将会讨论。 偏移函数 上面方程中, a(T) 代表偏移函数, 用来描述响应曲线的移动。正如后面所述, ANSYS 中可使用两个不同的偏移函数(以及一个用户自定义偏移函数)。,September 30, 20

7、01 Inventory #001491 7-11,粘弹性B. 流变模型,在详细讨论 ANSYS 中可用的粘弹性选项之前, 对一些常用的流变模型(1-D)进行总结有助于理解粘弹性行为。 下面的讨论依赖于两个基本模型, 弹簧和缓冲器 弹簧采用 Hooke 定律作为应力和应变的关系，这代表弹性固体。 缓冲器(或阻尼器)定义应力和应变率间的行为, 用于表示粘性流体。,September 30, 2001 Inventory #001491 7-12,粘弹性 . Maxwell 模型,Maxwell模型是一个串联的弹簧和阻尼器。 蠕变行为 松弛行为,September 30, 2001 Invento

8、ry #001491 7-13,粘弹性 . Maxwell 模型,Maxwell 模型与前面幻灯片所示的流体具有类似的特征。 呈现线性蠕变(在外加应力作用下连续变形)，线性蠕变不能代表大多数材料的第一阶段蠕变。 松弛时为指数应力响应。 因为线性蠕变, Maxwell 模型一般不能代表真实的行为。 速率效应如下: 快速加载极限: 弹性 缓慢加载极限: 没有抗力,September 30, 2001 Inventory #001491 7-14,粘弹性 . Kelvin-Voigt 模型,Kelvin-Voigt 模型是并联的弹簧和缓冲器。 蠕变行为 松弛行为,September 30, 2001

9、 Inventory #001491 7-15,粘弹性 . Kelvin-Voigt 模型,Kelvin-Voigt 模型有以下特征: 对蠕变, 在外载作用下, 应变是指数的并渐近至 so/E 在恒应变条件下, 没有可观察到的松弛行为，只有弹性部分在恒应变条件下对响应有贡献。 因为没有松弛行为, Kelvin-Voigt 模型一般不能捕捉材料真实的响应。 速率效应如下: 快速加载极限: 锁定 缓慢加载极限: 弹性,September 30, 2001 Inventory #001491 7-16,粘弹性 . 标准线性模型,标准线性模型(SLM) 是一个弹簧 与串联的弹簧/缓冲器的并联组合。,E

10、1,E0,h,这种表示法是 SLM 更常见的例子之一，然而, 也有其它模型。,September 30, 2001 Inventory #001491 7-17,粘弹性 . 标准线性模型,多数材料不能由 Maxwell 或 Kelvin-Voigt 模型本身充分地描述。 Maxwell 模型充分地描述了应力松弛而不是蠕变。 Kelvin-Voigt 模型充分地描述了蠕变而不是应力松弛。 SLM 提供了一个简单的蠕变和松弛行为的表示法。 速率效应如下: 快速加载极限: s=(E0+E1)e 缓慢加载极限: s=E0e,September 30, 2001 Inventory #001491 7-

11、18,粘弹性 . 广义 Maxwell 模型,广义 Maxwell 模型是由k 个并联的弹簧 和缓冲器组成.,ANSYS 中采用广义 Maxwell 模型 表示粘弹性行为，它是通用模型, Maxwell, Kelvin-Voigt 和 SLM 是其中的特殊情况。,September 30, 2001 Inventory #001491 7-19,粘弹性C. ANSYS 粘弹性模型,ANSYS 采用应力松弛函数的 Prony 级数表示法来模拟粘弹性。 这是广义 Maxwell 模型的数值实现过程。 仅能使用 VISCO88 (2D) 和 VISCO89 (3D) 单元类型。 这是一些高阶单元(能

12、使用退化形式)，虽然可能, 但不推荐用它们作为低阶单元(减去中间节点)。 目前 2D 平面应力、壳和梁单元对粘弹性不可用。 VISCO88/89 单元有应力刚化能力, 但假设为小应变、小位移行为。,September 30, 2001 Inventory #001491 7-20,粘弹性 . ANSYS 粘弹性模型,由TB,EVISC定义粘弹性材料模型，由材料GUI 或 TBDATA 命令输入数据。 通过常数 1-95 输入粘弹性特性数据。下面的讨论在适当的时候将用编号来表示这些常数。 用户应在适当的位置输入常数-常数1应在“C1”下输入, 常数5应在“C5”下输入, 等等。 材料常数与温度无

13、关, 虽然由于 TRS 假设, 这是不要求的。,September 30, 2001 Inventory #001491 7-21,粘弹性 . ANSYS 粘弹性模型,由材料 GUI 定义粘弹性模型的屏幕快照 Main Menu Preprocessor Material Props Material Models Structural Nonlinear Elastic Viscoelastic,September 30, 2001 Inventory #001491 7-22,粘弹性 . ANSYS 粘弹性模型,ANSYS 中考虑了与时间和温度有关的材料行为的三个方面: 剪切模量 (偏差刚

14、度) 体积模量 (体积刚度) 体积变化 (热膨胀) 如前所述, 偏移函数 a(T) 描述松弛响应曲线怎样关于温度移动。 Williams-Landau-Ferry 偏移函数用于聚合体。 Tool-Narayanaswamy 偏移函数应用于玻璃。 也可定义用户自定义虚拟温度/偏移函数。,September 30, 2001 Inventory #001491 7-23,粘弹性 . 剪切和体积模量,在粘弹性单元中, 体积和偏差行为是解耦的， 体积和剪切模量被独立地描述。 体积和/或剪切模量由下面的方程表示: 上面的方程中, “G” 指剪切或体积模量。 速率效应描述如下: 对短时间(x0), 材料特

15、性趋向于G(0)，这也可看作是固相特性。 对长时间(x ), 材料特性趋向于G()，这也可看作是液相特性。,September 30, 2001 Inventory #001491 7-24,粘弹性 . 剪切和体积模量,前面方程的常数输入如下: N 为 Maxwell 单元数, =10 G(0) 是初始剪切/体积模量(即固相) G(inf)或G()是最终剪切/体积模量(即液相) Ci 是剪切/体积模量松弛系数，使用时, 确保SCi = 1.0 常数li 是松弛时间,September 30, 2001 Inventory #001491 7-25,粘弹性 . 剪切和体积模量,为了提供前面讨论的

16、广义 Maxwell 单元的更直观的类比, 回顾关系 松弛时间是 li Ci(G(0)-G() = Gi G(0) = SGi G() = G0,September 30, 2001 Inventory #001491 7-26,粘弹性 . 剪切和体积模量,尽管剪切和体积模量松弛参数都能输入, 但用户不必为两个模量都指定松弛参数。 例如, 根据材料的不同, 可能很少或没有观察到的体积松弛，因此, 体积模量G(0) = G()，在这种情况下, 常数76-95不必输入。 必须 输入初始和最终体积和剪切模量，即使初始值和最终值是相同的。不存在缺省值，这些是常数C46-C49。 为了获得一个适当的或精

17、确的拟合, 用于拟合Prony 级数常数的常数数目N 由用户确定。 能够在每个角节点输出每个子步的有效(计算出的)体积和剪切模量。 参阅单元手册中的VISCO88 & 89, 输出数据一节。EFF BULK MOD 和 EFF SHEAR MOD 是在每个角节点输出的, 角节点可通过单元表获得。,September 30, 2001 Inventory #001491 7-27,粘弹性D. WLF 偏移函数 (C5=1),Williams-Landau-Ferry (WLF) 偏移函数通常用于聚合体，本节主要讨论这些材料。 对于很多聚合体, 若载荷在一固定值之下, 则变形仍是弹性的, 尽管整个

18、时间内发生能量耗散(由粘度表示)。 从分子角度看, 聚合体链节彼此相对移动，链间化学键在热激活过程中重新排列。 采用 WLF 偏移函数, 在ANSYS 中用粘弹性单元可以模拟该行为。,September 30, 2001 Inventory #001491 7-28,粘弹性 . WLF 偏移函数 (C5=1),WLF 偏移函数a(T) 描述如下: 设定常数C5=1, 选择 WLF 偏移函数 由常数 C1 和 C2 分别确定 C1 和 C2，它们是材料特性, 必须通过实验由不同温度下松弛时间来确定。 Tbase 为基础温度(在该温度下定义松弛特性)，这由常数 C4 输入。,September 3

19、0, 2001 Inventory #001491 7-29,粘弹性 . WLF 偏移函数 (C5=1),总之, 当使用 WLF 偏移函数时: 在给定温度 Tbase下定义剪切和/或体积模量的松弛参数, 正如在 C 节中讨论的那样。 基于不同温度下的实验松弛数据, 能够推导 WLF 常数 C1 和 C2，于是该偏移函数把Tbase 的松弛曲线和温度的变化相联系。,September 30, 2001 Inventory #001491 7-30,粘弹性E. TN 偏移函数 (C5=0),Tool-Narayanaswamy (TN) 偏移函数通常用于玻璃工业。下面主要讨论该偏移函数的用法。 在

20、加热/冷却条件下玻璃的粘弹性行为假设遵循TRS行为，然而与聚合体情况不同, 除结构松弛外, 还有明显的体积变化。 与液相有关的体积变化, 或生长应变, 通常远大于与固相有关的体积变化。 在玻璃转变区, 也有极强的时间相关性。原因是材料可能处于称作虚拟温度 Tf 的非平衡温度下, 该温度滞后于实际应用温度，这是与时间有极强的相关性的原因。,September 30, 2001 Inventory #001491 7-31,粘弹性 . TN 偏移函数 (C5=0),虚拟温度 Tf 用于说明在加热/冷却循环过程中生长应变与时间和温度有极强的相关性的函数。 N是体积衰减函数的 Maxwell 单元数,

21、 由常数3 输入。 由常数 6-15 输入 Cfi 由常数 16-25 输入松弛参数 tfi 对初始虚拟温度 Tf(0), 由常数 36-45 输入Tfi 虚拟温度 Tf 以 FICT TEMP 输出, 可由单元表获取。,September 30, 2001 Inventory #001491 7-32,粘弹性 . TN 偏移函数 (C5=0),增量生长应变(体积膨胀)可表示为: Degr 是生长应变 a(Tf)l 是基于虚拟温度 Tf 的液态热膨胀系数 a(T)g 是基于实际温度 T 的玻璃态热膨胀系数 总的(累积)生长应变 是对于子步总数Nt 而言的，当前(总)生长应变, 标记为 GR S

22、TRAIN, 可由单元表在每个子步输出，详情查阅单元手册。,September 30, 2001 Inventory #001491 7-33,粘弹性 . TN 偏移函数 (C5=0),前面提到的液相和玻璃相的热膨胀系数定义如下: 由常数 26-30 输入 Cli 由常数 31-35 输入 Cgi 注意 al 是虚拟温度 Tf 的函数, 而ag 是实际温度 T 的函数。,September 30, 2001 Inventory #001491 7-34,粘弹性 . TN 偏移函数 (C5=0),TN 偏移函数描述如下: H/R 为激活能除以理想气体常数R, 由常数C1输入。 X 代表参与定义

23、a(T) 的实际和虚拟温度值的数量, 由常数 C2 输入，根据定义, 0 x 1。 由TREF 命令输入参考温度 Tref。 记住指定 TOFFST 为绝对温度。,September 30, 2001 Inventory #001491 7-35,粘弹性 . TN 偏移函数 (C5=0),TN 偏移函数的特征总结如下: 虚拟温度 Tf 用来说明生长应变 (体积变化)与时间和温度的极强的相关性。 TN 偏移函数 a(T) 描述与时间和温度相关的剪切和/或体积模量的松弛参数, 如C节所述。这是对给定温度 Tref 定义的。同样, 定义 a(T) 时, 虚拟温度的一部分可以通过参数 x 与实际温度一

24、起使用。 TN 偏移函数涉及常数 H/R(材料的一种特性)，该偏移函数定义了温度变化时松弛模量的平移. 注意与 Arrhenius 方程的类似之处, 见第 4 章。,September 30, 2001 Inventory #001491 7-36,粘弹性F. 求解粘弹性模型,在 ANSYS 中进行粘弹性材料的分析时, 必须记住以下几点: 选择适当的单元类型 仅 VISCO88 和 VISCO89 可以用于粘弹性材料。 这些单元支持应力-刚化, 但仅适用于小位移、小应变分析。 定义粘弹性常数 除密度外, 仅需通过命令TB,EVISC 或 GUI 输入粘弹性常数。 指定温度时还需要规定TREF

25、和/或 TOFFSET。,September 30, 2001 Inventory #001491 7-37,粘弹性 . 求解粘弹性模型,使用粘弹性时还需记住: 指定适当的求解选项 需要指定一个足够小的时间步，这意味着在一个时间步内温度变化应非常小。同样, 在转变区, 时间步也应足够小以捕捉材料响应。 目前, CUTCONTROL 不支持粘弹性后退选项, 所以用户应保证时间步足够小以充分捕捉材料响应。 在通用后处理器中验证结果 后处理时, 粘弹性应变以“弹性应变” (EPEL) 给出。 当通过单元表对每个(存储的)子步的每个角节点进行输出时, 有效体积和剪切模量、虚拟温度和生长应变都可得到。,September 30, 2001 Inventory #001491 7-38,粘弹性 . 求解粘弹性模型,使用时粘弹性还需记住: 在时间-历程后处理器中验证结果: 下述量可以绘制为与时间的关系: 温度与时间 有效体积模量与时间 有效剪切模量与时间 虚拟温度与时间 生长应变与时间 上述时间-历程图, 若可用,