第8章 空间实体单元.pptVIP

第8章 空间实体单元 ; 空间离散化后的单元模型主要有：四面体单元、长方体单元、直边六面体单元、曲边六面体单元，如图8-1所示。 ; (b)图w为长方体单元，可以类似平面四节点矩形单元进行分析。 (c)图为任意八节点六面体单元，可以类似平面四节点任意四边形等参元进行分析。 (d)图为20节点曲边六面体单元，可以类似平面八节点曲边四边形等参元进行分析 。;单元节点位移向量为;（8-4）;（8-7）;在式（8-8）中，V为四面体的体积。为使其计算值不为负，单元的节点（i,j,m,n）编号次序应遵循右手法则。(p4) 采用同样的方法，可得;式中[N]为单元形函数矩阵，其维数为3×12。进一步可写为与平面问题式（2-21）、（2-22）类似的子块形式;在空间问题中，每点有6个应变分量。几何方程为： ;上式（8-15）、（8-16）与平面问题式（2-24）～（2-26）类似。与平面问题三节点三角形单元相同，在四节点四面体单元中，[B]的元素都是常量，因此是常应变单元。 ;（8-19）; 将式（8-20）所表示的[D]和式（8-15）、（8-16）所表示的[B]代入式（8-22），并将[S]定成分块矩阵的形式，有;式中 ;（8-26）;式中子矩阵[krs]为3×3的矩阵 ;（8-32）;如图8-3所示，在整体坐标系的二十节点六面体的实际单元与中心在局部坐标的原点、边长为2的立方体基本单元相对应。;（8-33）;根据几何方程，单元中的应变为;（8-38）;根据复合函数求导规则，有;3、应力矩阵;5、等价节点力矩阵;设某边界面上作用表面力 ，则等价节点力为 ;则式（8-44）可写为单元坐标系下的积分公式;8.4 空间轴对称单元;图8-4 轴对称构件及三角形环形单元;当然，也可以划分为矩形，构成矩形截面的环形单元。总之，单元截面形状，可采用三角形、四边形等平面有限元法所用单元形状。本节只讨论三角形截面的环形单元，其它截面形状环形单元可参照本节方法进行分析。 ;1、位移模式 ;（8-47）;（8-50）;与平面问题中???只含有3个应变分量不同，这里???含有4个应变分量。将u, w的表达代入式（8-51），得 ;由式（8-52）～（8-54）可见，矩阵中含有变量r, z，因此它不是常数矩阵。即轴对称问题的三角形环形单???不是常应变单元。 ;（8-56）;4、单元刚度矩阵 ;可见，单元刚度矩阵计算比平面三角形单元麻烦。须对其进行数值计算或近似计算。现在，讨论面积分?A g(r, z)drdz的计算问题。;（8-59）;则式（8-54）成为常数： ;其中，子矩阵可近似为 ;5、等价节点力 ;由式（8-64），有 ;以及Ni与Li等价的关系，利用积分公式（7-27），有 ; 由式（8-69）计算等价节点 力时，其中r由（8-66）式确定。并注意到在ij边上有Lm=0。利用式（7-28）进行积分，最后得 ;（8-70）;i点的集度为 ;