非线性基础教分解.ppt

July 30, 1999 非线性概述 第二章 2. 非线性概述 什么是“非线性”行为? 如果载荷引起刚度的显著变化, 则结构是非线性的. 刚度变化的主要原因是: 应变超出弹性极限 (塑性) 大挠度, 如承受载荷的钓鱼杆 两物体间的接触 …非线性概述 本章通过以下主题介绍非线性有限元分析(FEA)的基础: A. 什么是非线性? B. 用线性求解器求解非线性 C. 三类非线性 D. 非线性 FEA 难点 E. 一般的非线性分析过程 目的是了解非线性 FEA 的基本性质. 如何模拟各种非线性现象的更多细节包括在该书的其它章节和下面的两本书中: 高级接触和螺栓预紧培训手册 高级结构非线性培训手册 非线性概述A. 什么是非线性? 17世纪, Robert Hooke 发现了力(F)和位移(u)之间的简单线性关系, 即 Hooke 定律: F = Ku 常数 K 代表结构刚度. 线性结构服从这种线性关系. 常见的例子是简单弹簧: 非线性概述... 什么是非线性? 然而, 相当多结构的力和位移不是线性关系 因为这类结构的 F-u 曲线不是直线, 所以这样的结构称为非线性. 刚度不再是常数 K; 而成为所施加载荷的函数, KT (切向刚度). 常见的例子是韧性金属的拉伸实验: 非线性概述B.用线性求解器求解非线性 非线性分析中, 响应不能用一组线性方程式直接预报. 然而, 非线性结构可以用具有修正的线性逼近迭代系列进行分析. ANSYS 采用 Newton-Raphson方法 迭代. 每一次迭代称为平衡迭代. 非线性概述… 用线性求解器求解非线性 Newton-Raphson 法用下列方程迭代求解: [KT]{?u} = {F} - {Fnr} 式中: [KT] = 切向刚度矩阵 {?u} = 位移增量 {F} = 外载荷矢量 {Fnr} = 内力矢量 (单元应力的和) 每一次迭代都通过方程求解器独立求解. 一次迭代等价于一个线性静态分析. 非线性概述... 用线性求解器求解非线性 外载荷与内载荷之差{F} - {Fnr}称为残差, 是结构中力不平衡的度量. 目的是残差迭代至足够小, 即直到解收敛. 若收敛, 则解在允许的容差范围内平衡. 非线性概述... 用线性求解器求解非线性 但 Newton-Raphson 法不能保证在所有情况下都收敛! 仅当初始构形在收敛半径 内时 Newton-Raphson 才收敛. 非线性概述... 用线性求解器求解非线性 非线性概述... 用线性求解器求解非线性 如果初始构形在收敛半径外部, 有两种技术可帮助获得收敛解: 非线性概述... 用线性求解器求解非线性 一般的规律是系统任何方面的突变会导致收敛困难. 刚度突变. 载荷突变. 最佳收敛行为是把突变分成一系列很多小的递增的变化. 采用渐变加载. 采用小的时间步. 非线性概述B. 求解非线性 … 练习 请参考附加练习: W1. 非线性求解 – 浅拱形梁 非线性概述 C. 三类非线性 结构非线性由很多原因引起, 可归为三大类: 几何非线性 (大应变, 大挠度, 应力刚化) 材料非线性 (塑性, 超弹性, 蠕变) 状态改变非线性 (接触, 单元死活) 非线性概述 … 三类非线性 几何非线性 如果结构经历大变形, 则几何构形的改变会导致非线性行为. 非线性概述  … 三类非线性 材料非线性 非线性的应力-应变关系是非线性结构行为的常见原因. 非线性概述  … 三类非线性 状态变化非线性 由状态改变导致的刚度突变是非线性的另一常见原因: 电缆由松弛到张紧的状态改变 装配中两部件可以进入接触状态 加工可去除预应力化的材料 非线性概述  … 三类非线性 当然, 通常遇到三类非线性的组合. ANSYS 能处理组合的非线性效应. 非线性概述  D. 非线性 FEA难点 非线性有限元分析带来三个主要难点: 获得收敛 权衡代价与精度 验证 非线性概述  … 非线性 FEA难点 克服收敛困难 通常是最大的挑战. 求解必须在收敛半径内开始. 但无法确定收敛半径! 若解收敛, 初始点在收敛半径内. 若不收敛, 初始点在收敛半径外. 需要试错法获得收敛. 经验和训练可减少试错工作. 复杂的问题获得收敛需要多个载荷增量, 每个载荷增量需多次迭代. 多次迭代会增加总的求解时间. 非线性概述  …非线性 FEA难点 权衡代价与精度 所有 FEA 涉及代价(占用的时间, 磁盘和内存) 与精度的权衡. 更多细节和网格细化一般能得到更精确的解, 但需要更多时间和系统资源. 非线性分析增加一个附加因素--载荷增量数, 它影响精度和代价. 更