d)(线性与非线性区别

简介近十年以来，人们已不再将有限元分析(FEA)视为仅供分析师使用的工具，它已进入到实际的设计工作中。如今，CAD软件中都内置了FEA功能，设计工程师可使用FEA作为日常设计工具，协助完成产品设计过程。但是，直到最近，设计工程师所采用的大多数FEA应用程序还仅仅局限于线性分析。对于设计工程师所遇到的大多数问题，此类线性分析所得到的结果均与其实际特征大体接近。但是，有时也会出现需要采用非线性方法解决的更具挑战性的问题。过去，工程师们不愿意使用非线性分析，因为使用这种方法对问题进行公式表示非常复杂并且需要很长的求解时间。现在，随着非线性FEA软件与CAD结合，情况有所改观，软件的使用也更加简便。此外，改进的求解算法辅之以强大的台式计算机性能，使求解时间大大缩短。十年前，工程师将FEA视为极具价值的设计工具。现在，他们开始认识到非线性FEA的优点并更深刻地理解了它对设计过程所产生的影响。线性分析与非线性分析的区别术语“刚度”定义了线性分析与非线性分析间的根本区别。刚度是零件或装配体的特性，用于表征其对所施加载荷的反应。影响刚度的因素有很多：1.形状：工字形横梁与槽形横梁具有不同的刚度。2.材料：与相同尺寸的钢制横梁相比，铁制横梁的刚度较低。3.零件支撑方式：对于相同的横梁，在只有一端支撑时的刚度要比带有两端嵌入式支撑时的刚度低，弯曲程度更大，如图1所示。图1：悬臂横梁（上图）比具有两端支撑的相同横梁（下图）的刚度要低。术语“刚度”定义了线性分析与非线性分析间的根本区别。刚度是零件或装配体的特性，用于表征其对所施加载荷的反应。影响刚度的三个主要因素为：形状、材料和零件的支撑方式。如果刚度变化足够小，则可以假定在变形过程中形状或材料属性没有任何变化。此假设是线性分析的基本原理。了解非线性分析2当结构在载荷的作用下发生变形时，由于上述一个或多个因素的影响，其刚度会发生变化。如果变形很大，则其形状会发生变化。如果材料达到失效极限，则材料属性将会变化。另一方面，如果刚度变化足够小，则可以假定在变形过程中形状或材料属性没有任何变化。此假设是线性分析的基本原理。这意味着，在整个变形过程中，所分析的模型将保持施加载荷前尚未变形时所具有的刚度。无论模型的变形程度如何、载荷是一次施加还是逐渐施加、对载荷的反抗应力有多大，模型都将保持最初的刚度。此假设极大地简化了问题的公式表示和求解过程。回想一下FEA的基本方程式：[F]=[K]*[d]其中：[F]是已知的节点载荷向量[K]是已知的刚度矩阵[d]是未知的节点位移向量此矩阵方程式描述了FEA模型的状态。它包含大量的线性代数方程式，方程式的数量从数千到数百万不等，具体取决于模型大小。刚度矩阵[K]取决于几何形状、材料属性和约束。在线性分析假设下，模型的刚度永远不变，这些方程式只需组建并求解一次。当模型发生变形时，无需进行任何更新。因此，在线性分析中，从问题的公式表示到求解完成遵循的是一条直线道路。即使是非常庞大的模型，也只需要数秒钟或数分钟就能得到结果。在非线性分析的世界中，任何事物都会发生变化，因为非线性分析要求工程师放弃刚度不变的假设，而是认为刚度在变形过程中会发生变化，并且在迭代求解过程中，刚度矩阵[K]必须随非线性解算器的进展而更新。这些迭代运算增加了获得准确结果所需的时间。了解不同类型的非线性行为尽管刚度变化的过程对于各类非线性分析都相同，但非线性行为的来源可能不同。因此，按照造成非线性的主要来源对非线性分析进行分类是合理的。由于在很多问题中无法指出造成非线性行为的单一原因，因此某些分析可能必须解释多种非线性类型。几何非线性如前所述，当零件的刚度在操作条件下发生变化时，就需要采用非线性分析。如果刚度的变化只是由形状变化造成，则将该非线性行为定义为几何非线性。当零件产生肉眼可看到的较大变形时，就会发生此类由形状导致的刚度变化。一条广为接受的经验是，当变形大于零件最大尺寸的1/20时，应进行几何非线性分析。另一个需要注意的重要因素是，当发生较大的变形时，载荷方向会随着模型的变形而发生变化。大多数FEA程序会提供两种选择来解释此种方向变化：跟随载荷和非跟随载荷。跟随载荷保持与变形后模型的相对方向，如图2所示。而非跟随载荷保持其最初的方向。图2：跟随载荷（或称非保守载荷）在变形过程中会改变其方向并与变形的横梁保持垂直（左图）。非跟随载荷（或称保守载荷）保持其原始方向（右图）。一条广为接受的经验是，当变形大于零件最大尺寸的1/20时，应进行几何非线性分析。另一个需要注意的重要因素是，当发生较大的变形时，载荷方向会随着模型的变形而发生变化。了解非线性分析3压力容器在高压条件下发生巨大的变形是后一种情况的典型示例。该压力载荷总是垂直施加于压力容器的器壁。尽管对此种情况的线性分析假定容器的形状没有改变，但实际对压力容器进行分析时，需要分析在非保守载荷（或称跟随载荷）下的几何非线性。当变形很小时，也可能发生由于形状造成的刚度变化。一个典型的示例是，原本平坦的薄膜在压力的作用下发生弯曲（参见图3）图3：尽管变形的幅度非常小，但压力载荷下的平面薄膜仍需要采用非线性几何分析。最初，薄膜仅以折弯刚度对抗压力载荷。当压力载荷造成一定的曲率后，除原始折弯刚度外，变形的薄膜还会产生附加刚度（图4）。变形会改变薄膜的刚度，这使变形后的薄膜具有比平面薄膜具有更大的刚度。有些FEA程序使用令人混淆的术语，将所有几何非线性分析统称为“大变形分析”。此说法忽略了对较小变形执行非线性分析的必要性。图4：平面薄膜对载荷仅作出折弯刚度效应。由于发生变形，因此还产生了薄膜刚度。因此，薄膜的刚度比线性分析预测的要大。非线性材料如果刚度的变化仅仅是由于材料属性在操作条件下的变化造成，则该问题属于材料非线性。线性材料模型假定应力与应变成比例（下面的图5）。这意味着，其假定条件为施加的载荷越大，应力越大，并且变形将与载荷的变化成比例。它还假定不会产生永久的变形，一旦载荷消失，模型将总是恢复到其原始形状。图5当变形很小时，也可能发生由于形状造成的刚度变化。一个典型的示例是，原本平坦的薄膜在压力的作用下发生弯曲。最初，薄膜仅以折弯刚度对抗压力载荷。当压力载荷造成一定的曲率后，除原始折弯刚度外，变形的薄膜还会产生附加刚度。了解非线性分析4尽管这种简化是可以接受的，但是如果载荷高到足以导致某些永久变形（例如，大多数塑料所发生的变形），或者如果应变非常高（有时会50%，例如橡胶和弹性体所发生的应变），则必须使用非线性材料模型。由于不同种类材料在操作条件下的行为有很大的差异，FEA程序开发了专门的技术和材料模型来模拟这些行为。下表简要介绍各材料模型最适用的问题。弹性塑料超弹性材料粘弹性材料蠕变材料超弹性材料（形状记忆合金）这些模型适用于在达到最终应力前，其应变-应力曲线呈现“平稳”形状的材料。通过此材料模型可以很好地表征大多数工程用金属和一些塑料。此模型适用于土壤和颗粒状材料。适用于不能压缩的弹性体，例如橡胶。此模型适用于可压缩的聚氨酯泡沫橡胶。此模型适用于硬橡胶或玻璃。蠕变是指在恒定应力状态下产生的与时间相关的应变。大多数工程材料，尤其是高温下的金属、高分子塑料、混凝土和火箭发动机中的固体推进剂都会发生蠕变。形状记忆合金(SMA)（例如镍钛诺）表现超弹性效应。这种材料在装载-卸载周期内会发生较大的变形，但不会发生永久变形。VonMises或TrescaDrucker-PragerMooney-Rivlin和OgdenBlatz-Ko多种（可选择其他模型）多种（可选择其他模型）镊钛诺材料分类模型备注图6：理想弹塑性材料模型的应力-应变曲线。在此材料模型中，最大应力值不能超过塑性应力的极限（屈服应力）。在分析理想弹塑性材料模型（即在变形后完全失去恢复到原有形状的能力的材料）时，应力在超过特定应变值时将保持恒定。它描述了一个由8个螺栓很好地固定到位的铸铁材料挡盖（此模型是一个最简单的非线性材料模型，其应变-应力曲线如图6所示）。如果载荷高到足以导致某些永久变形（例如，大多数塑料所发生的变形），或者如果应变非常高（有时会50%，例如橡胶和弹性体所发生的应变），则必须使用非线性材料模型。了解非线性分析5线性分析表明，最大vonMises应力为614MPa(89,600psi)，而材料屈服力为206MPa(30,000psi)。线性分析的结果如图7所示。图7：线性应力求解表明，挡盖有非常高的局部应力集中。如果应力超过屈服力，挡盖将会破裂吗？要回答此问题，需要使用弹塑性材料模型来分析材料变为塑性的程度。图8显示当最大应力等于屈服应力时的非线性求解。塑性区域仍是局部的，这表明挡盖不会破裂。当然，需要通过严谨的工程评估来判断此设计是否可以接受。图8：通过理想弹性塑性材料模型得到的非线性应力求解。红色区域表示材料变为塑性。塑性区域的范围仍为局部。在分析理想弹塑性材料模型（即在变形后完全失去恢复到原有形状的能力的材料）时，应力在超过特定应变值时将保持恒定。了解非线性分析6图9显示了铝制支架的线性应力求解。最大应力为44MPa(6,400psi)，忽略材料在28MPa(4,100psi)下屈服这一事实。图9：带孔支架的线性应力求解表明应力远远超过了材料屈服应力的限制。非线性材料分析可以解释这些结果，即当最大应力保持在28MPa(4,100psi)时材料产生屈服（下面的图10）。非线性应力结果表明支架即将断裂。塑性区域几乎占据了悬臂的整个截面，稍微增加载荷量即可导致截面完全塑化，形成塑性铰，进而导致支架断裂。图10：非线性应力求解显示最大应力不超过屈服应力。塑性区域的范围表明支架即将形成塑性铰。它已处在载荷承受能力的极限。如果应力超过屈服力，挡盖将会破裂吗？要回答此问题，需要使用弹性塑性材料模型来分析材料变为塑性的程度。塑性区域仍是局部的，这表明挡盖不会破裂。了解非线性分析7在为普通钢制曲别针的简单动作（“拉直”然后“弯回”）建模时，需要同时考虑材料非线性分析和几何非线性分析。图11显示了使用理想弹塑性材料模型获得的曲别针变形后的形状。图12显示了在曲别针恢复原有形状后残余的应力。图11：曲别针弯曲分析需要采用材料非线性和几何非线性分析。“拉直”位置的曲别针显示塑性应力。图12：弯回到原有形状的曲别针显示残余应力。失去弹性稳定性（屈曲）零件的刚度还会因所施加的载荷而变化。有时，根据载荷施加方式的不同，载荷会增加刚度（拉伸载荷）或减小刚度（压缩载荷）。例如，绷紧的绳索可以承载杂技演员的体重，而松弛的绳索将会使其坠落。在压缩载荷的情况下，如果刚度的改变足以导致结构刚度降为零，即会发生屈曲并且结构会发生快速变形。然后，结构可能坍塌或在其屈曲后的状态中获得新的刚度。可通过线性屈曲分析计算导致结构发生屈曲的载荷（欧拉载荷）。但是，线性屈扭曲分析的结果是不保守的。此外，FEA模型中的理想化条件可能导致FEA模型预测的屈曲载荷远远高于实际失稳载荷。因此，应该谨慎使用线性屈曲分析的结果。屈曲并不一定意味着灾难性故障，在发生屈扭曲后结构仍有可能承载载荷。非线性分析将解释屈曲后的行为。了解非线性分析8图13和14显示了翘曲效应。零件甚至在发生扭曲后仍保持其承载能力。图13图14：屈曲效应分析需要采用非线性分析。接触应力和非线性支撑如果支撑条件（包括接触条件）在应用工作载荷时发生变化，则需要采用非线性分析。接触应力在两个接触面间形成。因此，接触面积及接触区域的刚度在求解前是未知的。图15显示了典型接触问题的应力求解。尽管与总体的模型尺寸相比，接触应力面积非常小，但是接触区域的刚度变化仍需要采用非线性分析。图15：为两个球面（只显示了两个接触零件中的一个）间产生的应力建模的接触应力分析属于非线性支撑分析类别。图16显示非线性支撑的示例。有效横梁长度及其后续的刚度取决于横梁的变形量。当横梁接触支撑物时，其刚度因有效支撑长度增加而增加。图16：此支撑（当启用时）会改变横梁的有效长度，从而导致横梁的刚度发生变化，因此该问题需要采用非线性分析。如果支撑条件（包括接触条件）在应用工作载荷时发生变化，则