弹性力学空间问题有限元法.ppt

《弹性力学空间问题有限元法.ppt》由会员分享，可在线阅读，更多相关《弹性力学空间问题有限元法.ppt（19页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、弹性力学空间问题有限元法现在学习的是第1页，共19页 当分析结构的形状、尺寸和边界条件不具备某种特殊性时，当分析结构的形状、尺寸和边界条件不具备某种特殊性时，这种结构便属于空间问题。这时结构内的任一点这种结构便属于空间问题。这时结构内的任一点i都具有三个都具有三个位移分量，即位移分量，即和六个应变分量和六个应变分量以及六个应力分量以及六个应力分量以上以上15个分量都是空间坐标（个分量都是空间坐标（x，y，z）的函数。）的函数。现在学习的是第2页，共19页空间问题：规模大、网格划分困难，网格划分占有较长时间。空间问题：规模大、网格划分困难，网格划分占有较长时间。分析这类结构时要充分利用求解问题的

2、特点（对称性、相似性和分析这类结构时要充分利用求解问题的特点（对称性、相似性和重复性等），尽量减少有限元模型规模。重复性等），尽量减少有限元模型规模。有些结构可以先近似为轴对称问题或平面问题进行粗算，然后再有些结构可以先近似为轴对称问题或平面问题进行粗算，然后再按空间问题求解。按空间问题求解。现在学习的是第3页，共19页工程中的绝大部分结构都属于空间问题工程中的绝大部分结构都属于空间问题这类结构的几何模型为三维实体模型，网格划分时可采用这类结构的几何模型为三维实体模型，网格划分时可采用四面体、五面体和六面体形式的空间实体单元。四面体、五面体和六面体形式的空间实体单元。现在学习的是第4页，共19

3、页第二节第二节 空间问题有限元法空间问题有限元法现在学习的是第5页，共19页一、结构离散一、结构离散在空间结构的离散中，常用的单元是四节点四面体单元。在空间结构的离散中，常用的单元是四节点四面体单元。zyxOpupvpwpmumvmwmiuiviwijujvjwj现在学习的是第6页，共19页1、位移函数、位移函数四面体单元，其节点为四面体的四面体单元，其节点为四面体的4个顶点，节点坐标分别为个顶点，节点坐标分别为i（xi,yi,zi），），j（xj,yj,zj），m（xm,ym,zm），p（xp,yp,zp）。每个节点有。每个节点有3个位移分量，个位移分量，即三个自由度，一个单元共有即三个自由

4、度，一个单元共有12个自由度。个自由度。二、单元分析二、单元分析单元的位移列阵为单元的位移列阵为现在学习的是第7页，共19页假设单元的位移函数为线性函数假设单元的位移函数为线性函数将将4个节点的坐标值和位移分量分别代入上式，即可求得待定系数个节点的坐标值和位移分量分别代入上式，即可求得待定系数 ，然后用形函数表达，然后用形函数表达现在学习的是第8页，共19页单元的位移可写成单元的位移可写成式中式中称为形函数矩阵。其中称为形函数矩阵。其中I是三阶单位矩阵；是三阶单位矩阵；称为形函数，称为形函数，系数为已知，系数为已知，V的四面体的体积。的四面体的体积。现在学习的是第9页，共19页2.单元应变矩阵

5、单元应变矩阵知道单元内各点的位移后，就可确定单元内任一点的应变，位移表达式代入知道单元内各点的位移后，就可确定单元内任一点的应变，位移表达式代入几何方程式，得几何方程式，得式中各矩阵元素为已知，应变矩阵式中各矩阵元素为已知，应变矩阵B是常量矩阵，因此这类单元是常应变是常量矩阵，因此这类单元是常应变单元。单元。应变矩阵应变矩阵B现在学习的是第10页，共19页3.单元应力矩阵单元应力矩阵由物理方程，得由物理方程，得式中式中S为应力矩阵，也是常数矩阵。为应力矩阵，也是常数矩阵。现在学习的是第11页，共19页4.单元刚度矩阵单元刚度矩阵由虚位移原理推导出单元的刚度方程由虚位移原理推导出单元的刚度方程写

6、成分块矩阵形式为写成分块矩阵形式为式中，单元刚度是由节点坐标和材料弹性常数确定的，它是一个常数矩阵。式中，单元刚度是由节点坐标和材料弹性常数确定的，它是一个常数矩阵。现在学习的是第12页，共19页三、总刚集成三、总刚集成 将每个单元在总体坐标系中的刚度矩阵进行叠加，将每个单元在总体坐标系中的刚度矩阵进行叠加，便可得到结构的总刚度矩阵。便可得到结构的总刚度矩阵。现在学习的是第13页，共19页四、载荷移置四、载荷移置 作用在结构的各类载荷同样需要移置为等效作用在结构的各类载荷同样需要移置为等效的节点载荷，移置方法与平面问题相同。的节点载荷，移置方法与平面问题相同。1、集中力的移置、集中力的移置 设

7、作用在单元上的集中力为设作用在单元上的集中力为由由 移置后产生的等效节点载荷为移置后产生的等效节点载荷为根据第二章介绍的方法，上式各矩阵元素的表达式为根据第二章介绍的方法，上式各矩阵元素的表达式为式中式中 是形函数在载荷作用点处的值。是形函数在载荷作用点处的值。现在学习的是第14页，共19页2、表面力的移置、表面力的移置 设作用在单元上的表面力为设作用在单元上的表面力为由由 移置后产生的等效节点载荷为移置后产生的等效节点载荷为则有则有式中式中 是形函数在载荷作用点处的值。是形函数在载荷作用点处的值。现在学习的是第15页，共19页因此，由各类载荷移置产生的单元总的等效节点载荷列阵为因此，由各类载荷移置产生的单元总的等效节点载荷列阵为由此可求出整个结构的节点载荷列阵，由此可求出整个结构的节点载荷列阵，现在学习的是第16页，共19页3、体积力的移置、体积力的移置 设作用在单元上的表面力为设作用在单元上的表面力为由由 移置后产生的等效节点载荷为移置后产生的等效节点载荷为则有则有因此，由各类载荷移置产生的单元总的等效节点载荷列阵为因此，由各类载荷移置产生的单元总的等效节点载荷列阵为现在学习的是第17页，共19页五、约束处理和求解线性方程组五、约束处理和求解线性方程组 现在学习的是第18页，共19页ANSYS 算例算例 4-1现在学习的是第19页，共19页