几何非线性有限元分析.ppt

可能必须将数据从一种度量转换为另一种. 输入数据要用正确的度量. 用一致的度量将 ANSYS 结果与已知的响应数据比较. 对于单轴应力-应变数据, 工程应力-工程应变可以通过以下公式转换为真实应力-对数应变: el = ln (1 + e) t = s (1 + e) 注意这种应力转换假设材料经历的大应变是不可压缩或几乎不可压缩的. 该假设对大塑性应变或超弹材料有效. 5.1 综述 用相反关系从真实转换为工程: 从工程应变转换为 Green-Lagrange: 相反关系是: 5.1 综述 5.2 三类几何非线性 一、大应变 : 当材料中的应变变“大”时(比如多于百分之几), 由于变形引起的几何形状改变不能再忽略了. “大” 是与问题相关的. 大应变分析不再假设应变是无穷小, 而是有限的或大的. 大应变理论考虑了形状改变 (如厚度, 面积等) 和任意大转动. 也固有地考虑了应力刚化效应. 大应变行为包括所有与几何非线性有关的三种现象. 大应变 大转动 应力刚化 二、大挠度 : 当一个单元的转动变“大”时 (比如大于1 到2 度), 单元的局部刚度转换为全局分量时将发生显著的改变. “大” 是与问题相关的. 在ANSYS中, 术语大挠度 和大转动 可以相互交换使用. 大挠度理论考虑了大转动, 但是它假设应变是小应变. 还固有地考虑了应力刚化效应. 大挠度理论是大应变理论的子集. 大应变 大转动 应力刚化 5.2 三类几何非线性 三、应力刚化: 一个零件中的应力状态会影响到该零件的刚度. 随着张力的增大, 电缆的横向刚度增加. 随着压缩量的增大, 柱体横向刚度下降 (最终导致完全丧失刚度 – 如, 屈曲). 当应力刚化被激活时, 程序计算应力刚度矩阵, 并将它添加到原始刚度矩阵去包含此效应. 应力刚度矩阵仅仅是应力和几何的函数. 应力刚度矩阵使切向刚度矩阵更加一致 (一般会改善收敛). 应力刚化理论是大挠度理论的子集. 5.2 三类几何非线性 大应变 大转动 应力刚化 因此, 在大位移分析中, 一个单元可用的特殊特征决定其行为: 大应变单元考虑形状改变和大转动, 并固有地包含应力刚化效应. 在这样的单元中, 不能将大应变效应和大转动效应分开. 大挠度单元考虑大转动, 并固有地包含应力刚化效应. 许多过去的杆、梁和壳单元有大挠度能力, 但是没有大应变( BEAM4, SHELL63 等). 将 应力刚化 作为特殊特征列出的单元, 缺省时在刚度矩阵中包括应力刚化项. 这一般有助于改善收敛速度. 5.2 三类几何非线性 程序几乎可以用任何形式的刚度矩阵, 收敛时仍然会获得精确解. 初始刚度 割线刚度 切向刚度 F u F u F u 精度由收敛容差决定, 而不是 [K] 的形式. 5.3 一致切向刚度矩阵 尽管精度不受刚度矩阵形式的影响, 但却严重影响收敛速度. 完全一致切向刚度矩阵 通常收敛速度最快. 完全一致切向刚度矩阵 [Kenl] 由四个分量组成: [Kenl] = [Keinc] + [Kes] + [Keu] + [Kea] [Keinc] 为主切向矩阵. [Kes] 为应力- 刚化矩阵. [Keu] 为初始位移- 转动矩阵, 它包括在刚度关系中改变几何形状的效应. [Kea] 为压力载荷刚度矩阵, 它包括在刚度关系中改变压力载荷取向的效应. 5.3 一致切向刚度矩阵 对于大多数单元的大位移分析, 缺省时自动包括前三个分量. 用求解控制去控制怎样包括第四个分量 (压力载荷刚度) : 缺省 (“程序选择”) 包括单元SURF153, SURF154, SHELL181, PLANE182, PLANE183, SOLID185, SOLID186, SOLID187, BEAM188和BEAM189的压力载荷刚度. 仅仅当遇到收敛困难问题时采用非缺省设置. 对不直接支持压力载荷刚度的单元, 可以通过在施加压力的表面上用 SURF154 包括该运算. 5.3 一致切向刚度矩阵 本章的剩余部分将描述在建立、运行和后处理大位移模型时经常有用的各种过程技巧. 首先研究建模的一些技巧. 单元选择. 网格划分. 耦合和约束方程. 5.4 建模 一、 用适当的单元类型： 并非所有的单元都支持几何非线性! 一些单元没有几何非线性能力. 例如 CONTAC52 和 PRETS179. 其他一些单元只有有限的几何非线性能力. 例如 VISCO88 不支持大应变或大转动, SHELL63 不支持大应变. 必须对计划采用的每一种单元类型检查单元描述中的特殊特征列表. 可以通过仔细地选择单元类型和单元选项来模拟特定的几何非线性行为。