第6章 非线性有限元法(几何非线性)详解.pptVIP

1、变形体的运动描述 ２、变形梯度张量 ２、变形梯度张量 ３、应变与变形测度 ３、应变与变形测度 ３、应变与变形测度 ４、大变形的应力测度 3、大变形的应力测度 ４、大变形的应力测度 ４、大变形的应力测度 ５、几何非线性有限元方程的建立 ５、几何非线性有限元方程的建立 ５、几何非线性有限元方程的建立 ５、几何非线性有限元方程的建立 ５、几何非线性有限元方程的建立 6、非线性平衡增量方程的线性化 7、几何非线性问题的有限元基本方程 6、U.L列式下平面杆元的几何非线性切线刚度矩阵 第六章 非线性有限元法（几何非线性） x3 x1 x2 P0 t0=0 tn+1=tn+Δtn tn Pn Pn+1 A0 An+1 An 变形体上的质点的运动状态可以随不同的坐标选取以下几种描述方法： 1、全拉格朗日列式法（T.L列式法—Total Lagrangian Formulation)： 选取t0=0时刻未变形物体的构形A0作为参照构形进行分析。 2、修正拉格朗日列式法（U.L列式法—Updated Lagrangian Formulation）： 选取tn时刻的物体构形An作为参照构形。由于An随计算而变化，因此其构形和坐标值也是变化的，即与t有关。tn为非线性增量求解时增量步的开始时刻。 3、欧拉描述法(Eulerian Formulation)： 独立变量是质点当前时刻的位置xn+1与时间tn+1。 几何非线性的有限元方程一般采用T.L或U.L列式法建立！ x3 x1 x2 P’ P 初始/未变形 变形后 位移u x’ x 1、首先采用Lagrangian方法，将一个物体的加载过程划分为一系列平衡状态。 位移方程 初始状态与变形后状态之间坐标关系为： 2、然后，考虑材料方向矢量，这个矢量描述物体内一段无限小的单元。 x3 x1 x2 式中，Fij称为变形梯度张量。 初始状态与变形后状态之间材料方向矢量的关系： 由位移方程，得： 由二阶张量特性，变形梯度张量的三个不变量为： 由于Fij表示从初始状态到变形后状态的一个映射，其逆映射Fij-1一定存在，即： 或写为： 体积映射: 面积映射： 变形前面积dA’ Ni (初始面积法向矢量) 变形后面积dA ni(变形后面积法向矢量) 映射Fij 逆映射F-1ij Fij是一个二阶张量。 由于变形梯度张量Fij中包含了刚体运动，因此不能直接用于定义应变测度。而材料方向矢量则不包含刚体运动，因此它的平方值可以作为衡量从某一状态到变形后状态的一个测度，定义为： 初始状态: 一个应变测度应该能反映出材料一段长度发生的改变。因此，应变张量可以由下式定义： x3 x1 x2 变形后状态： 提醒：由于Green应变张量表达式中的变形梯度张量对应于初始状 态，因此该应变张量也应在初始状态下计算。 ２、Almanshi应变张量 1、Green 应变张量 Green应变张量采用Lagrangian运动描述方法，即按初始状态下的构形定义应变张量。 式中，eij称为Green应变张量或Green-Lagrangian应变张量。 Almanshi应变张量采用Eular运动描述方法，即按当前状态下的构形定义应变张量。 式中，Eij称为Almanshi应变张量或Almanshi –Eular应变张量。 　由于大变形问题有限元方程主要采用T.L列式法或U.L列式法建立，因此应在初始状态下定义应变张量，即采用Green应变张量。 可以证明Green应变张量和Almanshi应变张量都是二阶对称张量。 2、Green – Lagrangian应变张量eij与小应变张量εij的关系 将变形梯度张量表达式代入到Green应变张量公式中，得： 式中： 为小变形应变张量； 2、Green变形张量也可写为： 为非线性二次项 1、Green应变张量 为小应变张量与一个非线性二次项之和，这意味所有大变形分析都是非线性的。 式中，Cij是Cauchy变形张量 由于Cauchy变形张量是正定对称阵，因此该张量有三个实特征值；这些特征值的平方根记为材料的主轴拉伸。 1、柯西应力张量(Cauchy’s stress tensor) 取三维空间笛卡儿坐标系，在t时刻的现时构形中截取一个四面体素，斜面的法线为n，另外三个面元与所取坐标面平行。由四面体素的平衡条件得出其上的应力为： 这里σij=σji便是柯西应力张量，它是二阶对称张量。 ①、柯西(Cauchy)应力张量是一种采用欧拉描述法(是以质点的瞬