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

State Key Laboratory of Structural Analysis for Industrial Equipment Dalian University of Technology, China 第二篇：弹性力学有限元的基本理论和格式 2.1 弹性力学有限元的一般格式 平面问题的有限元格式 弹性力学有限元的一般格式和求解步骤 有限元解的性质和收敛准则 2.2 单元与插值函数的构造 2.3 等参元与数值积分 等参元与数值积分 典型的等参单元 2.4 有限元法应用中的实际考虑 建立有限元模型： 计算结果的性质和处理 子结构法 对称性和周期性的利用 非协调元与分片实验 Isoparametric Element 数值积分 等参元中，单元面积积分转化为规则区域上的积分 积分的解析表达式有时难以给出；一般采用数值积分方法 Gauss积分法是数值积分的常用方法，适宜于有限元方法中计算积分 Gauss积分法将函数的积分近似为有限个点的函数值的加权和 Gauss积分 Gauss积分 Gauss积分 Gauss积分 积分点数目 两种数值积分策略：完全积分、减缩积分 完全积分：在单元内对所有的刚度系数进行精确积分（足够多点） 对于二次单元：高斯积分非常合适：多项式的乘积. Number of Integration Points The differences between the uniform strain formulation and the centroidal strain formulation can be shown as follows: Number of Integration Points 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)模式等 如果一个或多个变形模式在所有积分点都产生零应变，则出现不稳定性； 对于特定的变形模式，如果积分点不能感受到应变的存在，则单元刚度阵对该变形模式无法抵抗。 四节点平面单元的减缩积分 6.4.1 四节点四边形平面单元 四边形平面单元 单元分析 6.4.3 三维问题的有限元法 8 节点六面体等参单元 8节点六面体等参单元 8节点六面体等参单元 宝山壁画 宝山壁画是引人注目的昂贵文物。此壁画发现于阿鲁科尔沁旗东沙布乡境内。1994年列为“全国十大考古新发现”之一。宝山壁画中最引人注目的是《杨贵妃教鹦鹉图》。该画高0.7米、宽2.3米，用于笔重彩绘制，最突出的表现了 晚唐风格。唐代擅长绘贵妇仕女的大师周昉绘制了《杨贵妃教鹦鹉图》，不仅享誉中原，而且还影响全国各地。发现于阿旗宝山古墓里的这幅画，就是契丹人聘请中原画家按照周氏风格绘制的， 技法深得周氏画风的真传。在唐人真迹稀如星风的今天，能够从中完整了解唐代人物画的杰出成就，堪称美术史研究的辛事。这幅壁画现今保存在阿鲁科尔沁旗博物馆，历经千年，恍如新绘，是该馆的镇馆之宝。 欢迎大家观看！ 显然，该单元对引发这两种变形模式的力不能抵抗，系统刚度阵将是奇异的。 例：如图4个四节点四边形Q4单元，采用单点积分（高斯点在单元中心） 假设右上角单元发生如下两种变形模式（坐标原点取为单元中心）： （1） u=cxy, v=0; （2） u=0, v=-cxy 两种变形模式下，单元在高斯点处应变均为： 零能模式 (-1,1) (-1,-1) (1,-1) (-1,-1) 映射 坐标插值 位移插值 Ni 是相同的 x y 为了求出位移模式，以及局部坐标与整体坐标之间的变换式，在局部坐标系中定义一个四节点正方形单元（基本单元）。基