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

《有限元理论与方法-第5讲new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、讲授内容备注第5讲(第5周)附录1弹性力学的基本方程在有限单元法中经常要用到弹性力学的基本方程，关于它们的详细推导可从弹性力学的有关教材中查到。弹性体在载荷作用下，体内任意一点的应力状态可由6个应力分量，，，，，来表示。其中，，为正应力；，，为剪应力。应力分量的正负号规定如下：如果某一个面的外法线方向与坐标轴的正方向一致，这个面上的应力分量就以沿坐标轴正，方向为正，与坐标轴反向为负；相反，如果某一个面的外法线方向与坐标轴的负方向一致，这个面上的应力分量就以沿坐标轴负方向为正，与坐标轴同向为负。应力分量及其正方向见图1。图1应力分量应力分量的矩阵表示称为应力列阵或应力向量(1)弹性体在

2、载荷作用下，还将产生位移和变形，即弹性体位置的移动和形状的改变。弹性体内任一点的位移可由沿直角坐标轴方向的3个位移分量u，v，w来表示。它的矩阵形式是(2)称作位移列阵或位移向量。弹性体内任意一点的应变，可以由6个应变分量，，，，，来表示。其中，，为正应变；，，为剪应变。应变的正负号与应力的正负号相对应，即应变以伸长时为正，缩短为负；剪应变是以两个沿坐标轴正方向的线段组成的直角变小为正，反之为负。图2的(a)，(b)分别为和前正应变状态。7a正应变b剪应变图2应变的正方向应变的矩阵形式是(3)称作应变列阵或应变向量。对于三维问题，弹性力学基本方程可写成如下形式。1.平衡方程弹性体V域

3、内任一点沿坐标轴x，y，z方向的平衡方程为(4)其中，，为单位体积的体积力在x，y，z方向的分量。平衡方程的矩阵形式为（5）其中，A是微分算子，即是体积力矢量，。2.几何方程—应变-位移关系在微小位移和微小变形的情况下，略去位移导数的高次幂，则应变向量和位移向量间的几何关系有7(6)几何方程的矩阵形式为(7)其中L为微分算子，即(8)3．物理方程—应力-应变关系弹性力学中应力-应变之间的转换关系也称弹性关系（又称为广义虎克定律）。对于各向同性的线弹性材料，应力通过应变的表达式可用矩阵形式表示：(9)其中(10)称为弹性矩阵。它完全取决于弹性体材料的弹性模量E和泊桑比μ。表征弹性体的弹

4、性，也可以采用拉梅(Lam’e)常数G和λ：(11)G也称为剪切弹性模量。物理方程中的弹性矩阵D亦可表示为7(12)4.边界条件弹性体V的全部边界为S。一部分边界上已知单位面积外力，，称为力的边界条件，这部分边界用Sσ表示；另一部分边界上弹性体的位移，，已知，称为几何边界条件或位移边界条件，这部分边界用Su表示。这两部分边界构成弹性体的全部边界，即S=Sσ+Su(13)弹性体在边界上单位面积的内力为，，，根据平衡应有=，=，=(14)设边界外法线为N，其方向余弦为nx，ny，nz，则边界上弹性体的内力可由下式确定(15)式即为力的边界条件，矩阵形式为(16)在Su上有位移边界条件u=

5、，v=，w=(17)用矩阵形式表示为(18)以上是三维弹性力学问题中的一组基本方程和边界条件。同样，对于平面问题，轴对称问题等也可以得到类似的方程和边界条件。5.弹性体的应变能和应变余能单位体积的应变能（应变能密度）为(19)应变能是个正定函数，只有当弹性体内所有的点都没有应变时，应变能才为零。单位体积的余能（余能密度）为(20)余能也是个正定函数。在线性弹性力学中弹性体的应变能等于余能。第2章有限单元法基础理论73.1结构静力学问题的有限元法梁结构问题比桁架复杂一些，也可用矩阵分析法（线性代数方程组）得到问题的精确解。在ANSYS软件中上述两类问题的建模和求解较为简单。定义单元用E

6、T指令，考虑材料各种特性选不同的Link、Beam单元。有限单元法把杆系结构的矩阵分析方法推广应用于连续介质：把连续介质离散化，用有限个单元的组合体代替原来的连续介质，这样一组单元只在有限个节点上相互连接，因而包含有限个自由度，可用矩阵方法进行分析。3.1.1平面问题有限元法对一些特殊情况可把空间问题近似地简化为平面问题，只须考虑平行于某个平面的位移分量、应变分量与应力分量，且这些量只是两个坐标的函数。平面问题分平面应力问题和平面应变问题两类。设有很薄的均匀薄板，只在板边上受有平行于板面并且不沿厚度变化的面力，同时，体力也平行于板面并且不沿厚度变化，记薄板的厚度为t，以薄板的中面为x

7、y面，以垂直于中面的任一直线为z轴．由于板面上不受力，且板很薄，外力不沿厚度变化，可以认为恒有，，不为零的应力分量为、、，这种问题就称为平面应力问题。设有无限长的柱形体，在柱面上受有平行于横截面而且不沿长度变化的面力，同时，体力也平行于横截面且不沿长度变化。以任一横截面为xy面，任一纵线为z轴，由于对称性（任一横截面都可以看做对称面），此时，不为零的应变分量为、、，这种问题就称为平面应变问题。二维连续介质，用有限单元法分析的步骤如下：·用虚拟的直线把原介质