有限元理论与方法-第5讲.docVIP

讲 授 内 容 备 注 第5讲(第5周) 附录1 弹性力学的基本方程 在有限单元法中经常要用到弹性力学的基本方程关于它们的详细推导可从弹性力学的有关教材中查到。 弹性体在载荷作用下，体内任意一点的应力状态可由6个应力分量，，，，，，，， 图1 应力分量 应力分量的矩阵表示称为应力列阵或应力向量 (1) 弹性体在载荷作用下，还将产生位移和变形，即弹性体位置的移动和形状的改变。 弹性体内任一点的位移可由沿直角坐标轴方向的3个位移分量u，v，w (2) 称作位移列阵或位移向量。 弹性体内任意一点的应变，可以由6个应变分量，，，，，，，，，2的(a)，(b)分别为和前正应变状态。 a正应变 b 剪应变 图2 应变的正方向 应变的矩阵形式是 (3) 称作应变列阵或应变向量。 对于三维问题，弹性力学基本方程可写成如下形式。 1. 平衡方程 弹性体V域内任一点沿坐标轴x，y，z方向的平衡方程为 (4) 其中，，为单位体积的体积力在x，y，z方向的分量。 平衡方程的矩阵形式为 （5） 其中，A是微分算子，即 是体积力矢量，。 2. 几何方程—应变-位移关系 在微小位移和微小变形的情况下，略去位移导数的高次幂，则应变向量和位移向量间的几何关系有 (6) 几何方程的矩阵形式为 (7) 其中L为微分算子，即 (8) 3．物理方程—应力-应变关系 弹性力学中应力-应变之间的转换关系也称弹性关系（又称为广义虎克定律）。对于各向同性的线弹性材料，应力通过应变的表达式可用矩阵形式表示： (9) 其中 (10) 称为弹性矩阵。它完全取决于弹性体材料的弹性模量E和泊桑比μ。 表征弹性体的弹性，也可以采用拉梅(Lam’e)常数G和λ： (11) G也称为剪切弹性模量。物理方程中的弹性矩阵D亦可表示为 (12) 4. 边界条件 弹性体V的全部边界为S。一部分边界上已知单位面积外力，，称为力的边界条件，这部分边界用Sσ表示；另一部分边界上弹性体的位移，，已知，称为几何边界条件或位移边界条件，这部分边界用Su表示。这两部分边界构成弹性体的全部边界，即 S=Sσ+Su (13) 弹性体在边界上单位面积的内力为，，，根据平衡应有 =，=，= (14) 设边界外法线为N，其方向余弦为nx，ny，nz，则边界上弹性体的内力可由下式确定 (15) 式即为力的边界条件，矩阵形式为 (16) 在Su上有位移边界条件 u=，v=，w= (17) 用矩阵形式表示为 (18) 以上是三维弹性力学问题中的一组基本方程和边界条件。同样，对于平面问题，轴对称问题等也可以得到类似的方程和边界条件。 5. 弹性体的应变能和应变余能 单位体积的应变能（应变能密度）为 (19) 应变能是个正定函数，只有当弹性体内所有的点都没有应变时，应变能才为零。 单位体积的余能（余能密度）为 (20) 余能也是个正定函数。在线性弹性力学中弹性体的应变能等于余能。 第2章 有限单元法基础理论 结构静力学问题的有限元法 梁结构问题比桁架复杂一些，也可用矩阵分析法（线性代数方程组）得到问题的精确解。在ANSYS软件中上述两类问题的建模和求解较为简单。定义单元用ET指令，考虑材料各种特性选不同的Link、Beam单元。 有限单元法把杆系结构的矩阵分析方法推广应用于连续介质：把连续介质离散化，用有限个单元的组合体代替原来的连续介质，这样一组单元只在有限个节点上相互连接，因而包含有限个自由度，可用矩阵方法进行分析。 平面问题有限元法 t，以薄板的中面为xy面，以垂直于中面的任一直线为z轴．由于板面上不受力，且板很薄，外力不沿厚度变化，可以认为恒有，， 不为零的应力分量为、、，这种问题就称为平面应力问题。 设有无限长的柱形体，在柱面上受有平行于横截面而且不沿长度变化的面力，同