ADINA学习交流之——有限元非线性静力平衡方程组迭代算法.pdfVIP

ADINA 学习交流之 有限元非线性静力平 衡方程组迭代算法 （讲稿） 主讲人：胡明祎 博士 于2009-06-06 整理 ADINA 学习交流群学习资料 作者简介： 姓名：胡明祎 网名：淮水木鱼 胡明祎，男，中国地震局工程力学研究所防灾减 灾专业2006 级博士。 研究领域：贮液结构液固耦合抗震分析 结构抗震分析 土结相互作用 结构动力参数优化设计 有限元数值模拟等 ADINA 学习交流群学习资料 有限元非线性静力平衡方程组迭代算法 胡明祎 说明：首先感谢苦苦创造这么一次学习交流机会。因为最近一直很忙，讲义没有 详细做，只是拍了几张图片，又copy 了自己文章中的几幅图片进行说明的，后 面粘贴的命令流属于本人原创，但是仅仅是实验而已，所以整体讲义比较粗糙， 敬请谅解。希望通过这次交流能够和大家多多互换心得，旨在互相提高。 1 问题的提出 对于材料非线性问题，最终的整体结构的平衡方程组可以写成如 下的形式： K δ R [ ]{ } { } K B D B dV [ ] ∑∫[ ][ ][ ] 当材料处于非线性阶段时候，上面的静平衡方程组是如何进行求 解的呢？材料特性的随变形而变化，刚度矩阵 K 是如何变化的呢？ [ ] 在有限元数值计算过程中刚度矩阵 K 又是如何计算出来的呢？现在 [ ] 有哪些方法能有效地求解材料非线性影响下的有限元结构静力方程 呢？ ADINA 学习交流群学习资料 2 有限元常用的的静力平衡方程组求解方 法 2.1 增量法 2.1.1 纯增量法 2.1.2 初始刚度增量法 基于纯增量法，刚度矩阵取上一步增量分析结束时的切线刚 度矩阵，图略。 2.1.3 平均刚度增量法 1 ⎡ ⎤ ( =+ ) K K K [ ] [ ] i i+1 i ⎣ ⎦ 2 ADINA 学习交流群学习资料 2.1.4 中点刚度增量法 1 K {Δδ } (=ΔR ) [ ] i+1 i−1/2 2 δ δ =+ Δδ