应用有限元工程实例-4解答.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 对于Newmark法，在积分参数满足无条件稳定要求的前提下，时间步长的选择 主要考虑解的精度要求。 原则上说，可以根据对结构响应有主要贡献的若干个振型（比如 p 个）所对应 的周期来决定，即从这 p个周期中选出最小的周期 ，则 可以取 的若 干分之一，这是目前依据的主要原则。这里建议取 因为作为数值计算来讲，一个响应周期用20个点描述应该是已经比较精确的，当然，对于一个实际工程问题这样可能会导致计算量很大，所以，有的商用有限元程序建议取 。 同时，还要考虑外载荷的周期，即如果外载荷的周期小于 ，那么式 8.56 中 的 就应该是外载荷的周期。 另外，这里对于稳定性和计算精度的讨论都是针对线性动力问题而言的，而对 非线性问题，还有附加的因素需要考虑。 7 非线性有限元分析 7.1 非线性静力问题分类 非线性静力问题主要分为两类，即材料非线性和几何非线性。 材料非线性 材料非线性 材料非线性是由结构材料的非线性性能（如非线性弹性、弹塑性或粘弹性）引起的非线性力学行为，它反映在非线性的应力-应变关系（物理方程）中，在此情况下，材料的线性应力-应变关系不再适用。几种常见的应力-应变模型： 几何非线性 几何非线性是指在外力作用下，结构的位移比较大时引起的非线性力学行为。 产生的位移不再远小于结构自身的几何尺度，产生的应变不再远小于1，即应变应该包括位移的二次项。 较大的变形会导致结构的几何形状发生显著的改变，以至平衡方程必须按照变形后的位置来建立。 这种非线性反映在非线性的应变－位移关系（几何方程）中。 几何非线性弹性问题可以分为四类： 大应变问题 小应变、大位移问题 小应变、小位移，但位移的平方和应变大小同量级 小应变、小位移，但位移和应变大小同量级 7.2 非线性静力问题的求解 （1）非线性问题的总体方程 在线性情况下，总体刚度矩阵与应力或应变无关，即总体刚度矩阵是常数矩阵： 对于非线性问题（物理非线性和几何非线性），总体刚度矩阵与应力或应变相关，因为位移有限元法中，应力和应变都是可由位移表示的，所以总体刚度矩阵不再是常数矩阵，而是与位移相关，即： 可以看出，非线性的求解主要归属非线性方程组的求解。 在有些情况下，荷载项也可能与位移相关，这里主要介绍总体刚度矩阵与位移相关的情况。 （2）非线性方程组的求解方法 非线性方程组的求解方法有很多，常见的有：直接迭代法、牛顿法、拟牛顿法（修正牛顿法）、荷载增量法、弧长法等。 直接迭代法（也称割线刚度法）形式简单 牛顿法（也称切线刚度法）具有收敛速度快的优点。 修正牛顿法则提高了计算效率。 几种非线性有限元方程组的常用解法见下页图。 7.3 非线性静力问题的求解策略 在求解非线性问题时，需要确定荷载步长、收敛准则、收敛精度等。求解策略指如何选择、调整这些参数。 （1）荷载步和子步 下图为ANSYS程序中的荷载步和子步示意图 * 图中表示将载荷分成两个载荷步，即进行两次加载，同时又分别 为这两个载荷步定义了一个和两个子步。 * 子步还可以由程序根据计算情况自动判断的自动步长加载方法。 * 一般说，载荷步分得越多，即每次增加的载荷越少越容易收敛， 但是载荷步增加就会增加求解时间 * 载荷需要逐步增加，需要确定每步增加多少。 * 第一个载荷步的选定需要估算，一般可以做一次线性分析来确定 第一个载荷步。 （2）平衡迭代控制 对于线性分析，直接求解总体平衡方程就可以得到节点位移，对于非线性分析，则需要通过平衡迭代才能得到节点位移。 在那个载荷步进行迭代、迭代多少次都是可以控制的。比如：第一荷载步不进行平衡迭代，第二荷载步进行平衡迭代，限定次数为2次。 （3）收敛精度 收敛精度一般定义为 第i次迭代的增量的模 ? 第i-1次迭代的增量的模?? ?为小于1的实数。 收敛精度直接关系到求解精度，但无严格的控制标准，可以 选用程序给出的推荐收敛精度。 如能顺利算出结果，并使精度再提高一些，则可再提高收敛精度； 如不能算出结果，可适当降低收敛精度 （4）收敛准则（收敛判断标准） 几种收敛准则：位移收敛准则、力收敛准则、能量收敛准则等。 收敛准则的选取：位移、力相比较，哪一种量变化大，作为收敛的判断标准比较容易收敛。 在不易判断哪一种量变化大时，可用能量收敛准则。 （5）迭代结束 迭代结束指程序不再继续迭代，结束迭代一般会有三种情况： （a 收敛-结束：已经满足收敛精度。 （b 未收敛-结束：尚未收敛，但达到了允许的迭代次数。 如果增大允许的最大迭代次数，可能会收