橡胶疲劳≠金属疲劳 第2部分：线性叠加-1

疲劳分析中要考虑的载荷情况可能非常长，并且可能涉及多个载荷轴。通常，载荷工况比通过直接时域有限元分析（FEA）计算的载荷工况要长得多。

在金属疲劳分析中，线性叠加是一种广泛使用的技术，用于根据道路载荷生成应力应变历史[1]，[2]，[3]。当结构表现为线性时，这种方法是精确的和计算高效的，允许分析冗长的载荷信号。对于单轴问题，单个单元载荷情况的有限元（Fe）解根据输入载荷历史进行简单缩放。对于多轴问题，求解每个轴的单位载荷工况，然后根据输入载荷历史进行缩放和组合。

由于橡胶的1）非线性材料行为，2）非线性运动，3）非线性接触的可能性，线性叠加不能应用于橡胶疲劳分析。本文是研究橡胶疲劳分析程序与金属疲劳分析程序的不同之处的系列文章中的第二篇。在这里，我们介绍了Endurica EIETM（Endurica插值引擎）求解器，它是一种快速生成应力-应变历史的工具，用于在线性叠加失效的情况下进行疲劳分析。

1.金属线性叠加法的简要评述

对于线性结构，力[F]和位移[U]之间的关系可以写成矩阵乘法，其中[K]是刚度矩阵。

[F] = [k][u]

函数复合的结合性意味着将位移乘以标量a会产生成比例的较大力。

a[F] = [k](a[u])

加法的分布性质是指由组合位移[u]和[v]产生的力

[F] = [k][u] +[k][v]

也可以计算为

[F] = [k]([u] + [v])

类似地，可以通过线性叠加来缩放和组合应力和应变场。多年来，工程师们一直在金属疲劳分析中使用这一原理，特别是在处理由现场记录的载荷-位移历史引起的多轴情况时。

假设零件中的应力和应变场由单位载荷情况的线性组合产生，其中每个单位载荷情况的比例因子应用于给定输入通道对应的应力或应变场。

例如，对于图2所示的梁，如果通道1是大小为a（t）的单位位移u，通道2是结构中其他位置的大小为β（t）的另一个单位位移v，则可以通过线性叠加恢复梁中所有点的应力和应变的整个历史。

请注意，FE求解器只需为每个单位载荷工况生成一个与时间无关的解。完全通过比例因子a（t）和β（t）的时间变化来获得解的时间依赖性。这种极其有效的方法已在金属疲劳分析中使用多年。它允许快速分析由数百万个时间步长组成的完整道路载荷历史。

Endurica EIE：非线性情况下的载荷空间离散化和插值

解决非线性情况需要一种完全不同的方法。我们希望保留从预先计算的有限元解中有效构建应力-应变时间历程的优点。我们不是为每个输入通道预先计算单个单位载荷情况，而是从离散载荷空间中预先计算一组载荷情况。我们称这个集 合为地图。

图中载荷工况的数量必须足够，以便我们可以使用插值来获得图中任意点的非线性响应的合理近似值。图3显示了由x和z位移定义的两个通道的图。使用FE求解器如Ansys 或 LS-Dyna沿着蓝线所描绘的路径预先计算地图中的蓝点。一旦定义了映射，就可以从映射中预先计算的解中插入沿红线的应力-应变历史。

Endurica EIE 是一种通用工具，用于创建和使用非线性映射来生成疲劳分析的应力-应变历史[4]，[5]。EIE是Efficient Interpolation Engine（高效插值引擎）的缩写。EIE提供了一个简单的工作流程和强大的实用程序，用于创建和使用地图进行插值。它支持多达六个独立的输入通道。

整个EIE工作流程由三个主要步骤组成。第一步是创建地图。下一步是根据力或位移指定历史。请注意，任何可作为边界条件应用于有限元模型的量都可以设置为通道。最后一步是执行指定的插值。该过程为FE模型中的每个元素生成应变张量分量的时间历程。

地图创建过程包括四个步骤，如图4所示。首先，必须定义将用于指定历史的独立通道的数量。还必须指定映射类型。有几种类型可用，包括完全可自定义的地图。基于网格的地图通常适用于一维、二维和三维地图。对于高维映射，基于案例向量的映射通常是最方便的。

一旦定义了映射类型，EIE就会生成解决方案路径。这些由列举的载荷状态组成，这些载荷状态应作为边界条件应用于有限元模型，以生成映射。根据映射类型，可以生成一个或多个路径。每条路径称为一个分支。对于每个分支，EIE都会编写一个具有适当边界条件历史的文件，这是生成地图所必需的。接下来，使用EIE的边界条件建立并执行有限元模型。最后，FE结果的数据库链接到地图定义中的相应分支。

在这一点上，地图是完整的，并准备内插。请注意，当单位载荷情况解决方案被收集并定义为映射时，线性叠加可以作为EIE中的特殊情况实施。然而，一般来说，非线性映射将包含更多数量的求解步骤。

指定载荷历史与选择包含每个输入通道的时间历史的文件一样简单。在文件中，每行表示一个时间步长，每列表示一个输入通道。EIE 支持 .csv 和 .rsp 格式，这两种格式都是常见的数据格式。图5显示了具有和位移的示例历史。请注意，历史中的位移范围不应超过预先计算的标测图的范围。尽管内插解可能相当精确，但非线性问题的外推可能非常危险和不准确。

一旦指定了地图和历史，就可以开始插值。Endurica EIETM 支持多线程，这意味着插值计算可以在可用的CPU之间并行分布和执行。这使得插值速度非常快，并且可扩展到大型模型和冗长的历史。请注意，Endurica EIETM 生成的文件很大，因为它计算每个有限元的每个时间步长的应力和应变张量分量。因此，在运行Endurica EIETM 时，确保有足够的可用磁盘空间非常重要。

正文未完，请继续关注。

References

[1.] R. W. Landgraf, “Applications of fatigue a nalyses: transportation”, Fatigue ’87, vol. 3, pp. 1593–1610, 1987 

[2.] Moon, Seong-In et al, “Fatigue life evaluation of mechanical components using vibration fatigue a nalysis technique”, Journal of Mechanical Science and Technology, vol. 25, pp. 631–637, 2011. 

[3.] F. A. Conle and C. W. Mousseau, “Using vehicle dynamics simulations and finite-element results to generate fatigue life contours for chassis components”, International Journal of Fatigue, vol. 13(3), pp. 195–205, 1991. 

[4.] K. P. Barbash and W. V. Mars, “Critical plane a nalysis of rubber bushing durability under road loads”, SAE Technical Paper No. 2016-01-0393, 2016.

 [5.] W. V. Mars, “Interpolation engine for a nalysis of time-varying load data signals”. U.S. Patent 9, 645, 041, May 9, 2017. 

[6.] W. Mars,  K. Barbash et al, “Durability of Elastomeric Bushings Computed from Track-Recorded Multi-Channel Road Load Input”, SAE Technical Paper No. 2024-01-2253, 2024.