2014年有限元法及应用总结.pptVIP

对称性利用中的特殊问题 利用结构的对称性取某一部分建立有限元模型时，往往会产生约束不足现象。 例如，若取上例中图c的四分之一建立有限元时，根据上述分析，在两对称面上应加水平放置的滚动铰支座，因此模型在垂直方向存在刚体位移。对这种约束不足问题，利用有限元分析时，必须增加附加约束，以消除模型的刚体位移。在本例中，垂直方向可以用刚度很小的杆单元或边界弹簧单元连接到模型某节点上，使得既消除了模型的刚体位移，又不致于因附加的杆单元或边界弹簧单元刚度太大而影响结构原有的变形状态。 * 20.单元形态的选择原则 单元形态包括单元形状、边中节点的位置、细长比等，在结构离散化过程中必须合理选择。一般来说，为了保证有限元分析的精度，必须是单元的形态尽可能的规则。 对于三角形单元，三条边长尽量接近，不应出现大的钝角、大的边长。这是因为根据误差分析，应力和位移的误差都和单元的最小内角的正弦成反比。因而，等边三角形单元的形态最好，它与等腰直角三角形单元的误差之比为sin45°:sin60°=1:1.23。但是为了适应弹性体边界，以及单元由小到大逐渐过渡，不可能是所有的三角形单元都接近等边三角形。实际上，常常使用等腰直角三角形。 * 对于矩形单元来说，细长比不宜过大。细长比是指单元最大尺寸和最小尺寸之比。最优细长比在很大程度上取决于不同方向上位移梯度的差别。梯度较大的方向，单元尺寸要小些，梯度小的方向，单元尺寸可以大一些；如果各方向上位移梯度大致相同，则细长比越接近1，精度越高。有文献推荐，一般情况下，为了得到较好的位移结果，细长比不应超过7；为了获得较好的应力结果，细长比不应超过3。一般情况下，正方形单元的形态最好。 对于一般的四边形单元应避免过大的边长比，过大的边长比会导致病态的方程组。 20.单元形态的选择原则（续） * *21.边界条件的确定 确定边界条件是建立有限元模型的重要一环，合理确定有限元模型的边界条件是成功地进行结构有限元分析的基本要求。 一般情况下，建模对象的边界条件是明确的。根据力学模型的边界条件可以很容易确定其有限元模型的边界条件。例如电线杆插入地基的一端为固定端，桥梁一端为固定铰支座，另一端为滚动较支座。 但是，在机械工程中，建模对象往往是整个结构中的一部分，在建立有限元模型，确定其边界条件时，必须考虑其余部分的影响。这方面主要考虑如下两类问题。 * 1.边界位置的确定 在建立连续弹性体局部区域的有限元模型时，往往取该局部区域为隔离体，取其隔离边界条件为零位移约束，并通过试探校正确定零位移边界条件的位置。例如，进行齿轮齿有限元分析时，取一个轮齿的局部区域为隔离体，如图所示，设定PQRS的边界条件为零位移约束，通过改变边界深度PQ和边界宽度PS研究边界位置对齿根最大拉应力的影响，最后确定合理的边界条件。 * 2.边界条件的确定 有些分析对象的边界位置是零部件的连接部位。在建立有限元模型时，必须研究如何给定边界位置上的边界条件，以反映相连接结构的影响。确定这种问题的边界条件是用简单支撑连杆替代相连接结构的作用，使替代后结构的系统刚度等价于原结构的系统刚度。如分析机床主轴和传动轴时，可以利用等刚度的杆单元替代轴承和支座的作用，使轴的分析中包含有轴承和支座的影响。 * **22.单元和节点编号规则 当利用整体刚度矩阵的带状特征进行存贮和求解方程组时，单元节点编号直接影响系统整体刚度矩阵的半带宽，也就是影响在计算机中存贮信息的多少、计算时间和计算费用。因而，要求合理的节点编号使带宽极小化。半带宽的计算公式： 半带宽d=（单元节点号的最大差值+1）×节点自由度 由此，进行网格节点编号时应使网格中单元节点号的最大差值最小，这样才能保证半带宽最小。试比较下图。 * 图所示网格的四种编号方案中，单元节点标号的最大差值分别为5，3，5，9。显然，图2方案要合理。由此得出结论：沿着短边方向按列-列-列-列地顺序编号比沿着长度方向按行-行-行-行地顺序要合理（半带宽小） **22.单元和节点编号规则（续） * 然而，对于具有中间节点的单元或空间问题，须借助于带宽极小化的优化程序来对节点重新编号，先进的有限元程序包一般都配备有这样的程序。 对单元的编号只影响整体刚度矩阵的装配时间。由于这一时间在有限元运算时间中只占很小的比例，因而对于单元的编号并无特殊的要求。 **22.单元和节点编号规则（续） * **23.掌握分析三角形单元的位移模式求解方法 如图所示，在局部坐标系中，三角形平面单元的三个节点分别为1、2、3，其编号按逆时针方向进行，节点坐标分别为 * 24.求解平面问题中局部坐标系中的单元刚度矩阵 将几何方程和弹性方程代入虚功方程经整理后得：局部坐标系中 式中 —单元刚度矩阵 其中t——三角形单元平板的厚度， —三角形单元的