具有法向转动的膜单元.pptxVIP

第9章 具有法向转动的膜单元“梁、膜和壳单元之间的转动必须是协调的。”9.1 引言由于大多数建筑结构和土木工程结构的复杂特性，往往要求在同一计算模型中同时考虑框架单元、弯曲板单元和膜单元。为保证协调性，两个单元之间必须具有相同的自由度，即它们必须具有相同的个数和类型的平移自由度和旋转自由度。壳单元中，在每个节点上只有两个平面内的位移。对于结构系统中连接到边梁和许多其他普通类型壳单元，需要一种在每一个节点上具有一个法向转动作为基本自由度的膜单元。寻找具有法向转动特性膜单元的研究，在最近15年内，一种实用四边形单元得到了发展。本章将创建该单元的基本方程。9.2 基本假设 单元建立可以分为下列四个步骤：1.起点是图9.1中所示的8节点四边形单元（16个自由度）；2.将每边中间节点的相对位移旋转为垂直于和切向于相应边，并把相对切向位移设为零，从而将单元自由度减少到12个自由度；3.引入抛物线法线位移约束，以除去四个边中间节点法线位移，并在角部节点上引入四个相对法线转动；4.把相对法向转动转换成绝对值，并修改形函数以通过分片检验。9.3 位移近似法基本假设是平面内的x和y位移由下列方程所定义： 八个形函数可以从下列方程得出：9.4 节点转动介绍图9.2显示了一个典型的单元边ij.如果假设该边的相对法线位移是抛物线，则必须满足下列方程： 因为每边中间节点切向位移是零，所以整体相对每边中间节点位移可以从下列方程得出： 方程（9.4）可以应用于所有四条边和整体位移，方程（9.1）可以写为:9.5 应变—位移方程现在可以由下列基本方程构建应变—位移方程： 3*12个应变—位移方程可以使用下列子矩阵形式写成： 为使单元满足固定应力分片检验，必须对 矩阵做出以下修改：9.6 应力—应变关系正交各向异性平面应力材料的应力—应变关系可以写成： 应力—应变矩阵的唯一限制是它必须是对称且正定的。9.7 把相对转动转换为绝对转动具有法向转动的四边形单元的12*12刚度矩阵可以通过四点数字积分获得。即： 从方程（9.9）中计算出的膜单元刚度矩阵，在这些节点上有四个未知的相对转动。刚度矩阵的属性检验显示出它除了具有三个刚体模式外，还有一个零能量模式。图9.3中显示了相对于单元刚体旋转的伪变形模式。为消去此模式，只需把一个秩为1的矩阵加到含有与该模式相关的单元刚度矩阵中。根据转动的弹性定义，在单元中心的绝对转动，可由下列方程得到： 在单元中心的绝对转动与平均相对转动之间的差是： 将刚度k0指定给此变形，以便使用一点积分来创建下列秩为1的刚度矩阵： 很多问题求解的经验表明下列转动刚度数值是有效的： 当把该1秩矩阵添加到12*12刚度矩阵时，零能量模式可以被避免，节点转动可转换为绝对转动。9.8 三角形膜单元与四边形单元相同，把这一近似法应用于带有三个边中间节点的三角形单元，所得的刚度矩阵是9*9.四边形单元的计算机程序与三角形单元的计算机程序有大约90℅的相同之处，只是使用不同的形函数，并且省略与第四条边相关的约束。然而，三角形比四边形更具刚性。事实上，带有法向旋转自由度的三角形膜性能准确度与常规应变三角形几乎相同。9.9 计算范例图9.4中所示的梁是带有法向旋转自由度的两个膜单元模拟而成。表9.1中总结了位移和应力的结果。对于受到端点弯矩影响的矩阵单元，可以获得精确结果，而且不存在“剪力锁定”。随着单元扭曲的加重，位移和应力的情况更加恶化。9.10 小结本章介绍的膜平面应力单元，可以用于许多复杂结构系统的精确模型中，在这些结构系统中框架、膜和板单元间是相互连接的。分析中使用四边形单元可以产生理想的结果，但是相反，三角形膜单元的性能是很差的。