非线性有限元读书研究.docx

非线性有限元读书报告姓名：黄江益日期：2016年6月6日目录非线性有限元读书报告1目录21. 线性有限元回顾32. 非线性有限元概括32.1. 定义32.2. 使用分类32.3. 单元选择43. 非线性有限元求解方法63.1. 显式求解63.2. 隐式求解73.3. 求解步骤73.4. 收敛性84. 引用资料9线性有限元回顾核心点：影响结构材料性能有诸多因素（应力、应变、变形率、温度、湿度、时效等），而通常建立本构方程时仅考虑应力、应变两个物理参数，认为两者成线性关系的经典的弹性理论，即著名的虎克定律。在构建线性几何方程的时候，按照我们所学的课本所述，最为简单的是下面三种：平面应力（薄壁结构），平面应变（沿轴线的几何形状和外载无明显变化的物体），轴对称问题（广义的平面应变）。我们在构建平衡方程的时候，用的原则是虚功原理，即一个处于平衡的物体当发生满足位移边界条件的连续（微小）虚位移时，其外力虚功等于物体中应力在虚应变上产生的变形能。我们解线性有限元方程的时候，分为以下三个步骤：离散化单元分析（取位移插值函数，建立单元几何关系，这里假定的线性函数可以不同）总体分析（刚度叠加）我们还学了等参单元，尤其是基于坐标值变化获得的等参单元，最后利用大杀器高斯求积方法避免了复杂的积分函数的求解，极大降低了计算量。非线性有限元概括定义在有限元分析中的线性假设包含下列含义：结点位移为无限小量：应变——位移关系是线性的（几何方程）材料为线弹性：应力——应变关系是线性的（物理方程）加载时边界条件的性质保持不变：接触状态是已知且不随时间变化的（边界条件）而上述假设中的任何一条不能得到满足，那么就属于非线性有限元分析。结构力学问题，从本质上讲都是非线性的，线性假设只是实际工程问题的一种简化。使用分类在实际应用中，如图，主要遇到三类非线性问题：图 1 材料非线性图 2 几何非线性图 3 状态变化单元选择可以选择的有两种单元：线性单元和二次单元。这里依据的是ANSYS单元库和相关文档资料。图 4 （a）线性单元 （b）特殊形状的线性等参元 （c）二次单元选择线性单元（不带中间节点）的场合：在非线性分析中由于计算方程的复杂性，多采用线性单元。在对结构分析中，带有附加形函数的角点单元会在合理的计算时间内通常能得到准确的结果。当使用这些单元时，要注意防止在关键区域的退化形式。即避免在结果梯度很大或其它关注的区域使用二维三角形单元和楔形或四面体形的三维线单元。还应避免使用过于扭曲的线性单元，对于非线性结构分析，如果使用线性单元细致地而不是用二次单元相对粗糙的进行网格划分，那么将以很少的花费获得很好的精度。对弯曲壳体建模时，必须选用弯曲的（二次的）或平面（线性）的壳单元，每种选择都有其优缺点，对于多数的实际情况，主要问题利用平面单元以很少的计算时间，即可获得很高精度的结果。但是，必须保证使用足够多的平面单元来创建曲面。明显地，单元越小，准确性越好。推荐三维平面壳单元延伸不要超过 15 度的弧，圆锥壳（轴对称线）单元应限制在 10 度的弧以内（或 5 度如果离 Y 轴较近）。对多数非结构分析（热、电磁等），线性单元几乎与高次单元有同样好的结果，而且求解费用较低。退化单元（三角形和四面体）通常在非结构分析中产生准确结果。选择二次单元（带中间节点）的场合：在求解更高精度分析或者传统线性单元偏差较大时选用二次单元。对于用退化的单元形式进行的线性结构分析 （即二维三角形单元和楔形或三维四面体单元），二次单元通常会以比线性单元的求解费用更低且产生良好的结果。可是，为正确地使用这些单元，需要注意它们的特殊的性质。二次单元与常见的线性单元有诸多不同之处：对于分布载荷和边压力不象线性单元按“一般意义上”分配到单元节点上（见图5所示），单元的中间节点对反力也表现出相同的非直观的解释。由于中间节点的质量也大于角节点的质点， 所以通常将中间节点选为主自由度 （对于减缩自由度分析）。由于质量分配不均匀，在动力分析中，如果对波传波感兴趣，不推荐使用带中间节点的单元不要使用带中间节点的面单元定义以节点为基础的接触单元或连接间隙单元。当约束一个单元的边缘（或表面）自由度，面上所有的节点包括中间节点在内都要被约束。单元的角点只能与单元的角点相连，而不能与相邻单元的中间节点相连。相邻的单元应该有相连（或共同的）中间节点。图 5 载荷分配非线性有限元求解方法显式求解中心差分法是一种最为简单的求解方法，是用有限差分代替位移对时间的求导（对位移一阶求导得到速度，对位移二阶求导得加速度）。Runge-Kutta法，简称R-K法，是一种应用广泛的高精度单步法。 所谓单步法就是在计算yi时只用到前一步信息yi-1的方法。图 6 一种R-K法的求解结果比较优点：方法简单，容易实现，运算失效的可能性较小缺点：条件稳定