05-01 薄板弯曲问题的基本原理.docVIP

§4-2 空间等参数单元的数学分析 在进行空间等参数单元的力学分析时，需要用到 各个形函数对整体坐标的导数（求应变）； 局部坐标系中微分体的体积及微分面的面积（载荷移置、刚度矩阵）； 局部坐标面的法向余弦（载荷移置）。 现在来分别导出这些参数的表达式。 一、形函数对整体坐标的导数 由复合函数的求导法则，有 (x,y,z) 不过，由上节(4-17)和(4-18)式知，等参数单元的形函数Ni及x,y,z均只是局部坐标ξ，η，ζ的显函数，所以，利用上式无法求出形函数Ni对整体坐标x,y,z的导数。 我们如将Ni理解为局部坐标ξ，η，ζ的复合函数，则有 （ξ，η，ζ） 等等，所以有 (4-19) 从而有 (4-20) 这里 (4-21) 称为雅可比矩阵。 将(4-18)式代入上式。于是得 (4-22) 能够保证式(4-20)顺利运算的充要条件就是雅可比矩阵行列式｜J｜的值在整个单元上任意点均不为零，这也是等参单元可以进行坐标变换的必要条件。 具体地讲，这一条件就是前面所说的，在进行单元划分时，每个单元顶角的配置应限制在0°〈θ〈180°，一般应尽量控制在90°左右；节点的配置也应尽量和基本单元相似。这样，则既保证坐标变换的条件，又保证得到足够的计算精度。 二、局部坐标系中微分体的体积 在单元内的任一点P，沿局部坐标ξ、η、ζ的方向分别作微分矢量、、 (分别对应局部坐标的dξ，dη，dζ)，在实际空间形成一平行六面体(图4-9)。 由于在ξ方向，只是坐标ξ在变化，而η及ζ支持不变，所以微分矢量在整体坐标轴上的投影分别为 注：任一线段（dξ，dη，dζ）在XYZ上的投影为 ） 同样可得矢量及在整体坐标轴上的投影： 由矢量运算，可知此三矢量构成之平行六面体的体积为 此即为局部坐标系中微分体dξdηdζ映射到整体坐标系中相应微分体的体积。将上式化简，最后得。 (4-23) 即 (4-24) 式中｜J｜为前述雅可比矩阵的行列式，即 (4-25) 三、局部坐标系中微分面的面积 局部坐标系中平行于坐标 面的微分面dξdη，映射到整体坐标系中即为图4-9中矢量、所构成的平行四边形。 由矢量运算知，该四边形的面积为和矢量积的模，即 (4-26) 由前面导出的微分矢量在各坐标轴上的投影，极易求得： (a,b,c；ξ,η,ζ) (4-27) 其中 (ξ，η，ζ) 也可求得： (a,b,c;ξ,η,ζ) (4-28) 其中 (ξ，η，ζ) 将(4-27)式和(4-28)式代入(4-26)式，可得 (a,b,c;ξ,η,ζ) 此即为局部坐标系中微分面dξdη映射到实际坐标系中微分面的面积，即 (4-29) 四、局部坐标面的法向余弦 微分矢量与所构成平面的法线方向，即为矢量×的方向，也就是局部坐标面ξη的法线方向 (但一般并不是ζ坐标轴的方向，因为整体坐标下ξ，η，ζ一般并不是正交坐标)。 利用矢量运算公式，有 (4-30) 其模即为与所成平面的面积，也就是上面所求dAξη。因此，可由上列各式求得局部坐标面ξη的法向余弦表达式： (ξ，η，ζ) (4-31) 第五章 板壳弯曲问题的有限单元法 本章将讨论弹性板壳弯曲问题的有限单元法。为便于理解，首先介绍薄板弯曲问题理论。 弹性理论将板厚和平面尺寸相比很小的板称为薄板(工程上则将的板就认为是薄板)。当薄板受有垂直于板面(亦即垂直于中面)的外力时，横截面上将发生弯曲和扭转，引起弯曲应力和扭转应力。薄板的中面成为一个弹性曲面。这种问题称为薄板的弯扭问题，但通常简称为薄板的弯曲问题。 对于薄板小挠度(挠度)问题，经一定简化，可认为它的变形完全由横挠度所确定。因此，可以取和它的若干阶导数作为节点参数建立平板单元。目前已提出了非常多的板单元。下面将着重介绍比较常用的部分协调的矩形单元和一种三角形单元。虽然它们都不是完全协调的，但所得计算结果表明，它们都具有良好的收敛性和精确度。在这之后将介绍一种考虑了横向剪切变形影响的平板单元，它适用于薄板和厚板，由于是完全协调单元，它具有较好的性能。