空间问题与板壳单元

第五章空间问题与板壳单元第五章空间问题与板壳单元§5.1空间轴对称问题§5.2常用的体单元§5.3板壳问题§5.4算例5.1空间轴对称问题在用有限元法求解轴对称问题时，采用的单元一般为整圆环，如图所示。它们和子午面rz面相交的截面，可以是直边三角形或矩形，也可以是任意四边形、曲边三角形、曲边四边形等。各个单元之间以圆环形的铰相互连接，而每一个铰与子午面rz面的交点就称为节点。如图上的i、j、m等。所有单元将在子午面rz面上形成有限元网格，与在平面问题中形成的网格一样。因为在轴对称问题中采用的单元是一个整圆环，所以在计算单元的体积时要注意到这一点。下面以三角形截面单元为例，说明如何求解轴对称问题。如图所示子午面rz面上的一个三角形单元，设单元上任一点的径向位移（沿r向位移）分量为，轴向位移（沿z向位移）为量为。由于这两个位移分量仅是r和z的函数，故令与平面问题一样，可将位移用形函数及节点位移表示为

5.1.1位移模式即式中：——二阶单位矩阵；

、、——形函数矩阵，如下

（i、j、m轮换）式中：5.1.1位移模式由弹性力学知，轴对称问题中除平面内的应变分量外，还有环向应变。因此其几何方程为 将位移代入上式，得

5.1.2单元应变引入则

5.1.2单元应变应变矩阵中的元素不全是常量，因此单元内的应变也不是常量，这是因为轴对称问题中采用的单元是圆环，径向的位移必引起环向应变，而此应变的大小与点的位置有关。另外，由于中含有项，使单元的应变、应力及单元刚度矩阵的计算比平面问题复杂得多。根据弹性力学可知，对于轴对称问题，有

5.1.3单元应力根据虚功原理或用最小位能原理，可以和平面问题一样推得其单元刚度矩阵的表达式为，在轴对称问题中，由于单元是一圆环，上述积分式中的微分体积可取为微分圆环的体积，即，故单元刚度矩阵为与平面问题一样，单元刚度矩阵是一个阶的方阵，矩阵可分成三块，故也可分成个子矩阵，每个子矩阵为阶的方阵，其表达式为

5.1.4单元刚度矩阵因为矩阵与坐标有关，且坐标r处于分母上，因此积分不像平面问题中那么简单，常采用三种办法进行计算：1、显式积分；2、数值积分；3、简单的近似积分。一般采用第3种简单的积分，它不仅在程序上简单，而且还回避了节点在极轴上时带来的奇异问题。实践证明，在精度方面它并不比精确的积分公式法差。具体做法是取节点坐标平均值，即单元中心坐标5.1.4单元刚度矩阵并取在式刚度矩阵中以代替，以代替可得这样就使得单元刚度矩阵中的被积函数化为常数，然后积分即可求得，具体表达式这里不再给出。5.1.4单元刚度矩阵与平面问题一样，无论使用虚功原理或最小位能原理可以得到相同的载荷移置公式，其形式与平面问题相似。1.集中载荷轴对称问题中的集中载荷实质上是沿着圆周线作用，均匀分布的一圈力。在子午面单元上任一点处作用集中力，其中、为单位弧长上分布力的合力，其对应的等效节点载荷为5.1.5等效节点载荷

2.面力设单元的分布面力为，则其对应的等效节点载荷为3.体积力设单元是分布体力为，其中、为单位体积的体力分量，其对应的等效节点载荷为例如，在体力为自重的情况下，有、其中为容重，于是有5.1.5等效节点载荷

1.位移模式如图所示的四面体单元，单元节点的编码为i，j，m，n。每个节点具有三个位移分量，单元节点的位移列阵可表示为单元的位移模式采用线性多项式的形式，如下 5.2.1

四面体单元将单元上四个节点坐标代入上式，求解相应参数，则单元内任一点的位移可用节点位移和形函数表示，如下

用矩阵表示为式中：I——三阶单位阵；、、、——四面体单元的形函数，如下5.2.1

四面体单元

式中：（i、j、m、n轮换）

式中：V——四面体的体积。5.2.1

四面体单元2.单元应变根据几何方程，得

（i、j、m、n轮换）由于单元中的应变是常量，所以四面体单元是常应变单元。5.2.1

四面体单元3.单元应力由物理方程得单元的应力为 式中：——空间问题的弹性矩阵，如下5.2.1

四面体单元式中：——四面体单元的应力矩阵，表示为

（i、j、m、n轮换）式中：，，。由于单元中的应力是常量，所以四面体单元是常应力单元。5.2.1

四面体单元4.单元刚度矩阵由虚功原理，按照平面问题的类似推导，可得空间问题的单元刚度矩阵为

单元刚度矩阵可表示为分块形式，如下

式中任一子块由下式计算5.2.1

四面体单元式中：

5.2.1

四面体单元5.等效节点载荷（1）集中载荷若集中载荷为，则等效节点载荷可表示为（3）体力若分布体力为，则等效节点载荷可表示为5.2.1

四面体单元（2）面力对于单元的某一边界上的分布面力，则等效节点载荷可表示为单元自重的等效节点载荷可由上式计算，因为，所以等效节点载荷表示为由上式可知，单元等效节点力可看作将单元自重W均匀移置到每个节点上。5.2.1

四面体单元1.四面体单元由上节可知，常应变四面体单元中的各点应力为常量，然而实际工程结构中的各点应力是随着坐标变化的，为了反映真实情况，可以假定高次位移模式，相应的单元称为高阶单元。例如，在四节点四面体单元的基础上增加节点数，形成四面体单元族，如图所示。图中的4节点、10节点和20节点四面体单元，与左侧的四面体比较，可以看出节点的数目与多项式的项数完全一致，因此相应的位移模式可取为一次、二次和三次完全多项式。不难验证10节点二次四面体单元的应力、应变随着坐标线性变化，20节点三次四面体的应力、应变是坐标的二次函数，对于弹性力学问题其计算精度大大高于4节点常应变四面体单元。5.2.2高阶体单元高次四面体单元5.2.2高阶体单元2.六面体单元在空间问题中经常使用六面体单元，常见的六面体单元有8节点、20节点和32节点的立方体单元，如图5.5所示。8节点六面体单元精度明显高于同样节点数的四面体单元，20节点的六面体单元精度更高一些，但六面体形状难以适应工程结构的复杂外形，实际中使用不多。

空间六面体单元5.2.2高阶体单元板壳结构是工程结构中常见的一种结构，如图所示，中面为平面的称为板，中面为曲面的称为壳。本节以四节点矩形弯曲薄板单元为例，说明板壳单元的基本理论。板壳结构5.3板壳单元1.应变在板的小挠度理论中，板单元的的位移表示为

板单元的应变表示为式中：5.3.1基本方程2.应力应力应变关系则应力为式中：——弹性矩阵。5.3.1基本方程3.内力分量由弹性理论得式中：——薄板弯曲问题的弹性矩阵，如下

则薄板应力可由内力表示，如下5.3.1基本方程四节点矩形薄板单元如图所示，节点位移为5.3.2单元位移模式

（i、j、m、n轮换）矩形薄板有4个节点，整个单元的位移列阵为矩形薄板单元的位移模式可取为将薄板单元的坐标和节点位移代入上式，求解得到，并代回上式，则上式表示为式中：式中：5.3.2单元位移模式式中：

（i、j、m、n轮换）式中：5.3.2单元位移模式将位移代入应变表达式，可得到单元应变式中：

（i、j、m、n轮换）5.3.3单元应变与应力化简后得到单元应变矩阵，如下则单元应力可由物理方程求得，如下5.3.3单元应变与应力单元刚度矩阵用分块形式表示为将应变矩阵、弹性矩阵代入单元刚度一般表达式中，得单元刚度矩阵为式中子矩阵为3×3的矩阵，如下5.3.4单元刚度矩阵式中：5.3.4单元刚度矩阵式中：至于单元节点等效载荷的求法以及总体刚阵的组装均与平面问题是类似的，这里不再给出。把四节点矩形平面应力单元和该节讨论的四节点矩形弯曲薄板单元结合起来，可以构造同时考虑中面的伸展和板弯曲的四节点矩形板单元，进而分析壳体单元。5.3.4单元刚度矩阵如图（a）表示一个边长为1m的正方形薄板，四边固支，中心受垂直板面的集中载荷P作用，计算薄板中心挠度。边长l=1m，厚度h=0.01m，弹性模量E=210GPa，泊松比，载荷kN。（a）（b）5.4算例

解：由于是对称结构，取1/4板进行分析，共划分4个单元，9个节点，网格划分如图（b）所示。由于对称结构，则