非线性有限元基础.docVIP

§1.2 线性有限元的回顾 线性有限元的理论虽然不是本课程的重点，而线性有限元法是非线性有限元法的基础；非线性有限元法又是线性有限元法的发展。因为非线性问题的求解通常是采用分段线性化的思想，使其成为一系列的线性问题。因此有必要首先回顾线性有限元的一些基本内容。主要是线性有限元的基本理论和方法，以及当前应用最广的等参元理论。 固体力学的理论（材力、弹力、结构力学，无论是线性还是非线性）是建立在本构方程、几何（运动）方程以及平衡方程三方面方程的基础上的，由此导出的控制方程的解是满足上述充分、必要条件的唯一解，而且是反映了结构真实受载后的运动状态（变形）。在结构分析中许多情况下，本构和几何方程可以处理成线性，使控制方程线性化，求解也大大简化，同时可达到工程上的精度要求，这就是线性有限元的基本理论。 §1.2.1 线性本构关系（广义虎克定律） 影响结构材料性能有诸多因素（应力、应变、变形率、温度、湿度、时效等），而通常建立本构方程时仅考虑应力、应变两个物理参数，认为两者成线性关系的经典的弹性理论，即著名的虎克定律。 各向同性材料的Hooke定律 或 其中和分别为弹性阵和柔度阵。 由剪应力互等定理，弹性阵独立材料参数的个数由81个减少为21个。进而对于正交异性的材料参数为9个独立的参数，对于各向同性材料： 仅有两个独立的材料常数，即和，其中为弹性模量，为泊松比。 §1.2.2线性几何方程（小变形情况） 线性（小变形）关系： 位移边界条件：边界上 其中为位移向量，为边界上的指定位移，为微分算子。 实际结构可根据它们的几何特点，将三维问题简化为二维问题，这里主要有：平面应力、平面应变和轴对称状态。 1）平面应力（薄壁结构） 外力仅作用在平面内，两表面无外力作用（离心力作用下圆盘）两表面应力分量（且在整个厚度方向上），===0，而、、沿厚度均匀分布。 按Hooke定律， 得弹性阵 2）平面应变（沿轴线的几何形状和外载无明显变化的物体），如水坝，挡土墙等很长的物体，在离两端一定距离处，可以假定任意横截面上的位移、应力、应变等力学量仅是x、y的函数，沿物体纵向无变化，进一步假设沿z方向位移w为常数或为0，则有 3）轴对称问题（广义的平面应变） 用轴坐标代替直角坐标，（例如厚壁筒，高压罐，炮筒）。 ，，， §1.2.3 平衡方程（虚功原理） 经典力学给出了结构的平衡方程： or 边界条件 or 在上 为应力张量，为体力，为结构内力，为外部表面力，直接由平衡方程出发导出适合各类边界条件的有限元控制方程是困难的，但是从经典的虚功原理或加权余量原理出发，可得到广泛使用的有限元控制方程。 虚功原理认为：一个处于平衡的物体当发生满足位移边界条件的连续（微小）虚位移时，其外力虚功等于物体中应力在虚应变上产生的变形能。其表达式为： 其中为对应（虚位移）的虚应变，而虚位移应满足位移—应变的几何关系 和位移边界条件。以位移作为基本未知量的方程是当今有限元的主流称为位移法（位移模型），一个变形体总位能为：= ，式中是变形能，是外力功，按虚功原理 按照变分运算规则，仅对位移取变分，而力和应力不变，因此 ， 代入式即式。 按照上述能量原理，产生了多种求解结构位移的方法，如李兹法，伽辽金法等，有限元法与经典李兹法的一个重要区别是：有限元法不寻求在结构整个域上的连续位移函数，而是寻找各个子域（单元）上，满足域内及域边界上连续的位移函数。 设将物体分成若干个互不重迭的单元（包括三角单元，矩形单元等），对每个单元，设节点位移未知量，在单元内部和边界上满足位移协调条件的插值函数。于是体内位移场 代入几何关系可得 其中为式中微分算子。将其代入方程则有 进一步将本构关系代入，上式成为： 其中 单元刚度矩阵 作用于节点的当量体力 作用于节点的表面力 这里是指整个物体中所有单元求和，积分仅对每个单元域，这样在实际计算中对各单元，可用单元自身坐标