空间有限元法.pptxVIP

Definition空间问题的有限单元法2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-Definition 空间问题的有限单元法用有限单元法求解弹性力学空间问题，首先也要将连续的空间物体用一系列的单元离散化。空间问题中，最简单的是四面体单元。离散的空间结构是这些单元只在节点处以空间铰相互连接的集合体。2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-Exercise空间问题的有限单元法（续）2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-Definition位移模式空间问题中，每一个结点有3个位移分量，单元结点位移向量由12个分量组成，分别表示为：2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-位移模式（续）假定单元内的位移分量为坐标的线性函数:2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-位移模式（续）将上式中的第一式应用于4个结点，则有：2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-位移模式（续） 由上式可解出a1,a2,a3和a4再代回位移分量的表达式，可得： 式中： 为形函数，其中：2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-位移模式（续）2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-位移模式（续）用同样的方法，可以得到：合并υ,ν,ω的表达式，可以将单元内任一点的位移写为：2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-Definition单元中的应变和应力在空间问题中，每点有6个应变分量，由几何关系：2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-单元中的应变和应力(续)将υ,ν,ω的表达式代入上式，得到:式中：2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-单元中的应变和应力(续)可以看出，应变矩阵[B]中的元素都是常量，从而单元中的应变都是常量，故线性位移模式的四面体单元是常应变单元。由应力-应变关系，得到单元中的应力为：式中[D]为一般空间问题的弹性矩阵从下面[D]的表达式可以看出，单元中的应力都是常数。2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-单元中的应变和应力(续)2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-Definition 单元刚度矩阵和结点载荷向量仿照平面问题的推导，可以得到四面体单元的刚度矩阵：分块形式：2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-单元刚度矩阵和结点载荷向量（续）式中子矩阵可以表达为：其中：2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-单元刚度矩阵和结点载荷向量（续）经过与平面问题中同样的推导，单元的体积力向量和表面力向量可以用下列公式计算：经叠加，组合，得有限元支配方程：代入约束条件，可解出结点位移向量，从而就可以求出各单元的应变和应力。2001年10月1日ANSYS培训教程 – 版本 5.5 – XJTU MSSV (001128)Lc-