有限元与数值方法-讲稿8-等参元与数值积分.pptVIP

Number of integration points Numerical integration is simpler than analytical, but it is not exact. [k] is only approximately integrated regardless of the number of integration points Should we use fewer integration points for quick computation Or more integration points to improve the accuracy of calculations. Hmm…. 减缩积分 减缩积分通常指在计算单元内力和单元刚度阵时采用的数值积分方案减小一阶 有限元模型通常是过刚的。采用比较多的积分点，通常使得单元过刚更为严重，因为附加的积分点将反映单元刚度阵中的更高阶项，而这些项抵抗某些变形模式的发生 减缩积分可减小单元刚度阵的计算量 但是，采用过少的积分点可导致由于刚度阵缺秩引起的更为严重的问题，包括：不稳定性、多余的零能模式、沙漏(hourglass)模式等 如果一个或多个变形模式在所有积分点都产生零应变，则出现不稳定性； 对于特定的变形模式，如果积分点不能感受到应变的存在，则单元刚度阵对该变形模式无法抵抗。 显然，该单元对引发这两种变形模式的力不能抵抗，系统刚度阵将是奇异的。 四节点平面单元的减缩积分 例：如图4个四节点四边形Q4单元，采用单点积分（高斯点在单元中心） 假设右上角单元发生如下两种变形模式（坐标原点取为单元中心）： （1） u=cxy, v=0; （2） u=0, v=-cxy 两种变形模式下，单元在高斯点处应变均为： 零能模式 6.4.1 四节点四边形平面单元 (-1,1) (-1,-1) (1,-1) (-1,-1) 映射 坐标插值 位移插值 Ni 是相同的 x y 为了求出位移模式，以及局部坐标与整体坐标之间的变换式，在局部坐标系中定义一个四节点正方形单元（基本单元）。基本单元的位移模式为 局部坐标与整体坐标的变换式也取为 四边形平面单元 单元分析 对形函数求导，可得到应变的矩阵表达式 可得应变矩阵的分块矩阵， 由于等参单元的形函数是局部坐标 的函数，因此应变矩阵[B]也是局部坐标 的函数。形成等参单元的单元刚度矩阵需要在整体坐标系中对局部坐标的函数进行积分，包括以下三个基本步骤 （1）计算用局部坐标表示的形函数对整体坐标x,y的偏导数； （2）将整体坐标系中的面积积分转换为在局部坐标系中的面积积分； （3）用数值积分计算出单元刚度矩阵中的元素 （1）计算形函数对整体坐标x，y的偏导数 或 定义雅可比矩阵（Jacobian Matrix） 则有 （2）将整体坐标系中面积积分转换为在局部坐标系中的面积积分 （3）用数值积分计算出单元刚度矩阵中的元素 等参单元刚度矩阵的每个元素都是局部坐标的函数，通常采用高斯积分方法计算单元刚度矩阵中的元素。 将作用在单元上的外载荷同样表示为局部坐标的函数，就可以在局部坐标下完成单元的载荷移置。体力移置的公式为 在点 集中力移置的公式为， 6.4.2 八节点四边形等参单元 四结点四边形等参单元的局限性： 不能准确拟合曲线边界 位移模式的阶次仅为双线性 八节点四边形等参数单元具有更高阶的精度 基本单元 ? ? ?= -1 ?= 1 ?= -1 ?= 1 1 3 ? ? 5 6 7 8 2 4 ? ? 实际单元 x y 1 2 3 4 5 6 7 8 ? ? ? ? ? ? ? ? ?= -1 ?= 1 ?= -1 ?= 1 基本单元的位移模式为： 即 八节点四边形等参单元 坐标变换为 基本单元 ? ? ?= -1 ?= 1 ?= -1 ?= 1 1 3 ? ? 5 6 7 8 2 4 ? ? 实际单元 x y 1 2 3 4 5 6 7 8 ? ? ? ? ? ? ? ? ?= -1 ?= 1 ?= -1 ?= 1 以上形函数还可表示为如下形式： 八节点四边形等参单元 双二次形函数 显然，该坐标变换式将??平面上的正方形映射为xy平面上的曲边四边形。xy平面上每一条边都是一条二次曲线，它完全由对应边上3个结点的坐标唯一确定。因此，单元是协调的，同时也可证明，单元的位移函数反映刚体位移和常应变，具有完备性。 有关八节点四边形单元的特性分析和等价节点力计算过程与四节点四边形单元完全相同，具体公式形式也一致。区别仅在于两种单元有关矩阵的维数不同，具体是： （1）四节点单元的???e是一个8×1列阵；而