第2章-结构几何非线性.ppt

第2章 第2章 结构几何非线性 1.概述 在线弹性力学分析中，假定位移与应变关系是线性的，且应变为小量，由此而得到线性几何方程。当考虑位移与应变的非线性关系或采用大应变理论（有限变形理论）则都属于几何非线性问题，亦即非线性问题包括了大位移、小应变以及大位移、大应变等问题，此时均导致几何运动方程成为非线性。但材料的本构关系还是符合虎克定律，结构的弹性稳定问题是结构非线性分析的内容之一。 几何非线性问题至今尚未完全成熟，仍然是一门正在迅速发展与完善的课题。这主要表现在建立非线性有关基本方程方面存在不同学派的争论，他们各有优缺点，并未得到一个权威性的结论；由于大变形引起载荷的变动（非保守系统）对方程与解的影响问题研究很少，也无明确、相对肯定的结论，大多数问题的局限在保守系统内；大应变下，应变不可叠加性的讨论与研究也不成熟；几何非线性解法发展尚处在活跃阶段，尚未找到一种十分满意的适应性广、收敛快的解法。 二.变形体的运动描述 图示一变形体在to=0时有构形Ao，物体中一质点Po的坐标为（x10，x20，x30)，在t=tn时，物体有运动构形An，质点Po运动至Pn，在时间tn+1=tn+△tn时，物体运动有构形An+1，质点运动至Pn+1，对于变形体及其上的质点运动状态，可以随不同的坐标选取有以下几种描述方法： 全拉格朗日列式法（ T.L列式法－ Total Lagrangian Formulation）。选取to=0时刻未变形物体的构形Ao作为参照构形进行分析。 修正的拉格朗日列式法（ U.L列式法－ Updated Lagrangian Formulation）。选取tn时刻的物体构形An为参照构形。由于An随计算而变化，因此其构形和坐标值也是变化的，即与t有关。tn为非线性增量求解时增量步的开始时刻。 尤拉描述法（Eulerian Formulatlon） 独立变量是质点P当前时刻的位置xn+1与时间tn+1。 现在用得最为广泛的是Lagriangian列式，因此下面主要讲述T.L和U.L下几何非线性有限元方程的建立。 三. 线弹性杆单元刚度矩阵建立的步骤 1.选择插值函数描述体内任一点的位移 2. 根据几何方程及位移函数确定应变矩阵[B] 3. 求单元的总势能∏ 4. 根据势能驻值原理求单元刚度矩阵[k]  求和约定 Kroneker符号 连续体变形的描述方法 Lagrange法－用各质点在初始位置的坐标作为独立变量进行描述。即以变形前的初始构形为基准，然后确定它与变形构形间的相对变形，导出的应变张量称为Green应变张量。 Euler法－用各质点位置的即时坐标或我们需要的时刻的坐标作为独立变量进行描述。即以变形构形为基准，然后确定其余初始构形间的相对变形，由此导出的应变张量成为Almansi应变张量。 在连续介质力学中，变形体的运动描述仅分为Langrange和Euler两种方法，UL列式法归结到Euler列式法中。 Lagrange 描述－Green应变张量 Euler 描述－Almansi应变张量 2.杆元的几何运动方程 基于工程应变定义，其内隐含着l与l0的方向一致，因此其适用于转动可以忽略不计的小转动情形。 对于小位移和位移导数，两种应变张量及工程应变均趋于同一。 工程应变不具可加性，而对数应变具备可加性   对数应变－hanchi应变具有可加性 不同杆元运动方程的比较 图示杆元绕原点发生一刚性转动，转角为 ，因杆元仅发生刚性转动，故应变应为零，下面来考察不同的应变表达式对这一状态的描述结果。 1.7 ×10-2 7.7 ×10-3 1.9× 10-3 1.2 ×10-3 6.9 ×10-4 3.0 ×10-4 7.6 ×10-5 -5.8 ×10-4 -1.2 ×10-4 -7.2 ×10-6 -3 ×10-6 -9.4 ×10-7 -1.9 ×10-7 -1.2 ×10-8 -3.4 ×10-2 -1.5 ×10-2 -3.8 ×10-3 -2.4 ×10-3 -1.4 ×10-3 -6.1 ×10-4 -1.5 ×10-4 15° 10° 5° 4° 3° 2° 1° 由此可以看出： 五.T.L列式下平面杆元的 几何非线性切线刚度矩阵  1.杆元位移函数 2.几何方程 3.应变方程与应变矩阵 4.物理方程 5.杆元总势能 6.杆元切线刚度矩阵[K]T 六.U.L列式下平面杆元的 几何非线性切线刚度矩阵 1.杆元位移函数 注意下列推导过程中，杆元的节点位移和节点力实际为相应杆元的位移和节点力的增量。 * o x1 x3 x2 t0=0 tn tn+1=tn+△tn A0 An An+1 ·P0 ·Pn ·Pn+1 四. 杆元非线性几何运动方程1