ANSYSWORKBENCH疲劳分析指南

1、第一章 简介1.1 疲劳概述  结构失效的一个常见原因是疲劳，其造成破坏与重复加载有关。疲劳通常分为两类：高周疲劳是当载荷的循环（重复）次数高(如1e4 -1e9)的情况下产生的。因此，应力通常比材料的极限强度低，应力疲劳（Stress-based）用于高周疲劳；低周疲劳是在循环次数相对较低时发生的。塑性变形常常伴随低周疲劳，其阐明了短疲劳寿命。一般认为应变疲劳（strain-based）应该用于低周疲劳计算。  在设计仿真中，疲劳模块拓展程序（Fatigue Module add-on）采用的是基于应力疲劳（stress-based）理论，它适用于高周

2、疲劳。接下来，我们将对基于应力疲劳理论的处理方法进行讨论。1.2 恒定振幅载荷  在前面曾提到，疲劳是由于重复加载引起：  当最大和最小的应力水平恒定时，称为恒定振幅载荷，我们将针对这种最简单的形式，首先进行讨论。  否则，则称为变化振幅或非恒定振幅载荷。1.3 成比例载荷  载荷可以是比例载荷，也可以非比例载荷：  比例载荷，是指主应力的比例是恒定的，并且主应力的削减不随时间变化，这实质意味着由于载荷的增加或反作用的造成的响应很容易得到计算。  相反，非比例载荷没有隐含各应

3、力之间相互的关系，典型情况包括：  1/2=constant  在两个不同载荷工况间的交替变化；  交变载荷叠加在静载荷上；  非线性边界条件。1.4 应力定义    考虑在最大最小应力值min和max作用下的比例载荷、恒定振幅的情况：    应力范围定义为(max-min)    平均应力m定义为(max+min)/2    应力幅或交变应力a是/2    应力比R是min/

4、max    当施加的是大小相等且方向相反的载荷时，发生的是对称循环载荷。这就是m=0，R=-1的情况。    当施加载荷后又撤除该载荷，将发生脉动循环载荷。这就是m=max/2，R=0的情况。1.5 应力-寿命曲线    载荷与疲劳失效的关系，采用的是应力-寿命曲线或S-N曲线来表示：    （1）若某一部件在承受循环载荷, 经过一定的循环次数后,该部件裂纹或破坏将会发展，而且有可能导致失效；    （2）如果同个部件作用在更高的载荷下,导致失

5、效的载荷循环次数将减少；    （3）应力-寿命曲线或S-N曲线,展示出应力幅与失效循环次数的关系。    S-N曲线是通过对试件做疲劳测试得到的弯曲或轴向测试反映的是单轴的应力状态，影响S-N曲线的因素很多，其中的一些需要的注意，如下：    材料的延展性，材料的加工工艺，几何形状信息,包括表面光滑度、残余应力以及存在的应力集中，载荷环境，包括平均应力、温度和化学环境，例如，压缩平均应力比零平均应力的疲劳寿命长，相反，拉伸平均应力比零平均应力的疲劳寿命短，对压缩和拉伸平均应力，平均应力将分别提高和降低S

6、-N曲线。    因此，记住以下几点：一个部件通常经受多轴应力状态。如果疲劳数据(S-N 曲线)是从反映单轴应力状态的测试中得到的，那么在计算寿命时就要注意：（1）设计仿真为用户提供了如何把结果和S-N曲线相关联的选择，包括多轴应力的选择；（2）双轴应力结果有助于计算在给定位置的情况。    平均应力影响疲劳寿命，并且变换在S-N曲线的上方位置与下方位置(反映出在给定应力幅下的寿命长短)：（1）对于不同的平均应力或应力比值，设计仿真允许输入多重S-N曲线(实验数据)；（2）如果没有太多的多重S-N曲线(实验数据)，那么设计仿真也允许采

7、用多种不同的平均应力修正理论。    早先曾提到影响疲劳寿命的其他因素，也可以在设计仿真中可以用一个修正因子来解释。1.6 总结    疲劳模块允许用户采用基于应力理论的处理方法，来解决高周疲劳问题。    以下情况可以用疲劳模块来处理：    恒定振幅，比例载荷(参考第二章)；    变化振幅，比例载荷(参考第三章)；    恒定振幅，非比例载荷(参考第四章)。    需要输入的数据是

8、材料的S-N曲线：    S-N曲线是疲劳实验中获得，而且可能本质上是单轴的，但在实际的分析中，部件可能处于多轴应力状态。    S-N曲线的绘制取决于许多因素，包括平均应力，在不同平均应力值作用下的S-N曲线的应力值可以直接输入，或可以执行通过平均应力修正理论实现。第二章 疲劳程序2.1 基本情况指定材料特性，包括S-N曲线；    定义接触区域(若采用的话)；    定义网格控制(可选的)；    包括载荷和支撑；  &#

9、160; (设定)需要的结果，包括Fatigue tool；  由于有线性静力分析,所以需要用到杨氏模量和泊松比：如果有惯性载荷,则需要输入质量密度；如果有热载荷,则需要输入热膨胀系数和热传导率；如果使用应力工具结果(Stress Tool result),那么就需要输入应力极限数据,而且这个数据也是用于平均应力修正理论疲劳分析。 疲劳模块也需要使用到在工程数据分支下的材料特性当中S-N曲线数据：数据类型在“疲劳特性”(“Fatigue Properties”)下会说明；S-N曲线数据是在材料特性分支条下的“交变应力与循环”(“Alternating Stress vs.

10、Cycles”)选项中输入的。    添加和修改疲劳材料特性:    在材料特性的工作列表中,可以定义下列类型和输入的S-N曲线，插入的图表可以是线性的（“Linear”）、半对数的（“Semi-Log”即linear for stress, log for cycles）或双对数曲线（“Log-Log”）。  记得曾提到的，S-N曲线取决于平均应力。如果S-N曲线在不同的平均应力下都可适用的，那么也可以输入多重S-N曲线，每个S-N曲线可以在不同平均应力下直接输入，每个S-N曲线也可以在不同应力比下输入

11、。  可以通过在“Mean Value”上点击鼠标右键添加新的平均值来输入多条SN曲线。  材料特性信息可以保存XML文件或从XML文件提取，保存材料数据文件，在material条上按右键，然后用“Export ”保存成XML外部文件，疲劳材料特性将自动写到XML文件中，就像其他材料数据一样。   对于应力分析的任何类型结果，都可能需要用到：应力、应变和变形接触结果(如果版本支持)；应力工具（Stress Tool）。    在疲劳计算被详细地定义以后，疲劳结果可下在Fatigue T

12、ool下指定；等值线结果（Contour）包括Lifes(寿命)，Damage(损伤)，Safety Factor（安全系数），BiaxialityIndication（双轴指示），以及Equivalent Alternating Stress（等效交变应力）；曲线图结果(graph results)）仅包含对于恒定振幅分析的疲劳敏感性(fatigue sensitivity)；这些结果的详细分析将只做简短讨论。2.2 Fatigue Tool  当Fatigue Tool在求解子菜单下插入以后，就可以在细节栏中输入疲劳说明：载荷类型可以在“Zero-Based”、“Ful

13、ly Reversed”和给定的“Ratio”之间定义；也可以输入一个比例因子，来按比例缩放所有的应力结果。    Goodman理论适用于低韧性材料，对压缩平均应力没能做修正，Soderberg理论比Goodman理论更保守，并且在有些情况下可用于脆性材料，Gerber理论能够对韧性材料的拉伸平均应力提供很好的拟合，但它不能正确地预测出压缩平均应力的有害影响，如下图所示。    缺省的平均应力修正理论可以从“Tools > Control Panel：Fatigue>Analysis Type”中进行设置如果存在多重S-

14、N曲线，但用户想要使用平均应力修正理论，那么将会用到在m=0或R=-1的S-N曲线。尽管如此，这种做法并不推荐。 2.3 求解疲劳分析    疲劳计算将在应力分析实施完以后自动地进行，与应力分析计算相比，恒定振幅情况的疲劳计算通常会快得多。如果一个应力分析已经完成，那么仅选择Solution或Fatigue Tool 分支并点击Solve，便可开始疲劳计算。在求解菜单中（solution branch）的工作表将没有输出显示，疲劳计算在Workbench中进行，ANSYS的求解器不会执行分析中的疲劳部分，疲劳模块没有使用ANSYS /POST1的疲劳命令(FSxx

15、xx, FTxxxx)。2.4 查看疲劳结果    对于恒定振幅和比例载荷情况，有几种类型的疲劳结果供选择：    Life（寿命）：等值线显示由于疲劳作用直到失效的循环次数，如果交变应力比S-N曲线中定义的最低交变应力低，则使用该寿命（循环次数）(在本例中，S-N曲线失效的最大循环次数是1e6，于是那就是最大寿命。    Damage（损伤）：设计寿命与可用寿命的比值，设计寿命在细节栏（Details view）中定义，设计寿命的缺省值可通过下面进行定义“Tools > Control Panel

16、:Fatigue > Design Life。      Safety Factor（安全系数）：安全系数等值线是关于一个在给定设计寿命下的失效，设计寿命值在细节栏（Details view）输入，给定最大安全系数SF值是15。    BiaxialityIndication：应力双轴等值线有助于确定局部的应力状态，双轴指示（Biaxialityindication）是较小与较大主应力的比值(对于主应力接近0的被忽略)。因此，单轴应力局部区域为B值为0，纯剪切的为-1，双轴的为1。    

17、  等效交变应力(Equivalent Alternating Stress）：等值线在模型上绘出了部件的等效交变应力，它是基于所选择应力类型，在考虑了载荷类型和平均应力影响后，用于询问（query）S-N曲线的应力。疲劳敏感性( Fatigue Sensitivity )：一个疲劳敏感曲线图显示出部件的寿命、损伤或安全系数在临界区域随载荷的变化而变化，能够输入载荷变化的极限(包括负比率)，曲线图的缺省选项，“Tools menu > OptionsSimulation：Fatigue>Sensitivity”。       任何疲

18、劳选项的范围可以是选定的部件（parts）和/或部件的表面，收敛性可用于等值线结果。收敛和警告对疲劳敏感性图是无效的，因为这些图提供关于载荷的敏感性(例如，没有为了收敛目的而指定的标量选项)。    疲劳工具也可以与求解组合一起使用，在求解组合中，多重环境可能被组合。疲劳计算将基于不同环境的线性组合的结果。2.5 总结   a 建立一个应力分析(线性,比例载荷)  b 定义疲劳材料特性,包括S-N曲线  c 定义载荷类型和平均应力影响的处理     d 求解和后处理

19、疲劳结果Solve and postprocessfatigue results第三章 不稳定振幅的疲劳    在前面一章中,考察了恒定振幅和比例载荷的情况,并涉及到最大和最小振幅在保持恒定的情况下的循环或重复载荷。在本章将针对不定振幅、比例载荷情况,尽管载荷仍是成比例的,但应力幅和平均应力却是随时间变化的。3.1 不规律载荷的历程和循环(History and Cycles)    对于不规律载荷历程,需要进行特殊处理：    计算不规律载荷历程的循环所使用的是“雨流”rainflow循环计算，“雨流”

20、循环计算(Rainflowcycle counting)是用于把不规律应力历程转化为用于疲劳计算的循环的一种技术(如右面例子)，先计算不同的“平均”应力和应力幅(“range”)的循环，然后使用这组“雨流”循环完成疲劳计算。    损伤累加是通过Palmgren-Miner 法则完成的，Palmgren-Miner法则的基本思想是：在一个给定的平均应力和应力幅下，每次循环用到有效寿命占总和的百分之几。对于在一个给定应力幅下的循环次数Ni，随着循环次数达到失效次数Nfi时，寿命用尽，达到失效。    “雨流”循环计算和Palmgren-

21、Miner损伤累加都用于不定振幅情况。    因此，任何任意载荷历程都可以切分成一个不同的平均值和范围值的循环阵列(“多个竖条”)，右图是“雨流”阵列，指出了在每个平均值和范围值下所计算的循环次数，较高值表示这些循环的将出现在载荷历程中。    在一个疲劳分析完成以后，每个“竖条”(即“循环”)造成的损伤量将被绘出，对于“雨流”阵列中的每个“竖条”(bin)，显示的是对应的所用掉的寿命量的百分比。在这个例子中，即使大多数循环发生在低范围/平均值，但高范围(range)循环仍会造成主要的损伤。依据Per Miner法则，如果损伤累加到1

22、(100%)，那么将发生失效。3.2 不定振幅程序  a 建立引领分析(线性,比载荷)   b 定义疲劳材料特性(包括S-N曲线)   a 定义载荷历程数据,并以及平均应力的影响的处理 b 为“雨流”循环次数的计算定义bins的数量    e 求解并查看疲劳结果(例如，损伤matrix，损伤等值线图，寿命等值线等)    对于建立基于不定振幅、比例载荷情况下疲劳分析的过程，与前面讲过的第二章中介绍非常相似，但有两个例外：载荷类型的定义不同，查看的疲劳结果中包括变化的

23、“雨流”和损伤阵列。    前面的讨论非常清楚地指出“bins”的数目影响求解精度。这是因为交互和平均应力在计算部分损伤前先被输入到“bins”中。这就是“Quick Counting”技术。    默认方法（因为其效率高）“Quick RainflowCounting”可以在“Details view”中关闭，在这种情况下，部分损伤发现前数据不会被输入到“bins”，因此“bins”的数目不会影响结果。    虽然这种方法很准确，但它会耗费更多的内存和计算时间。3.4 查看疲劳结果  

24、;  定义了需要的结果以后,不定振幅情况就可以采用恒定振幅情况相似的方式，与应力分析一起或在应力分析以后进行求解。由于求解的时间取决于载荷历程和竖条尺寸，所在进行的求解可能要比恒定振幅情况的时间长，但它仍比常规FEM的求解快。结果与恒定振幅情况相似：    代替疲劳循环次数,寿命结果报告了直到失效的载荷块的数量。举个例子，如果载荷历程数据描述了一个给定的时间块(假设是一周的时间)，以及指定的最小寿命是50，那么该部件的寿命就是50块或50周。    损伤和安全系数(Damageand Safety Factor)基于在Det

25、ails栏中输入的设计寿命(Design Life)，但仍然是以块形式出现,而不是循环。    BiaxialityIndication(双轴指示)与恒定振幅情况一样，对于不定振幅载荷均可用。    对于不定振幅情况，Equivalent Alternating Stress(等效交变应力)，不能作为结果输出。这是因为单个值不能用于决定失效的循环次数，因而采用基于载荷历程的多个值。    Fatigue Sensitivity(疲劳敏感性)对于寿命块也是可用的。    在不定

26、振幅情况中也有一些自身独特的结果：    Rainflow阵列，虽然不是真实的结果，对于输出是有效的，在前面已经讨论了，它提供了如何把交变和平均应力从载荷历程划分成竖条的信息。    损伤阵列显示的是指定的实体(scoped entities)的评定位置的损伤。它反映了所生成的每个竖条损伤的大小。注意，结果是在指定的部件或表面的临界位置上的结果。   第四章 非比例载荷的疲劳分析    在第二章中,讨论了恒定振幅和比例载荷情况，本节将针对恒定振幅非比例载荷情况进行介绍。其基本

27、思想是用两个加载环境代替单一加载环境，进行疲劳计算，不采用应力比，而是采用两个载荷环境的应力值来决定最大最小值。由于同一组应力结果不并不成比例，这就是为什么这种方法称为非比例(non-proportional)的原因,但是两组结果都会使用到，由于需要两个解，所以可以采用求解组合来实现。    对于恒定振幅,非比例情况的处理过程与恒定振幅、比例载荷的求解基本相同,除了下面所提出的以外:    1.建立两个带不同载荷条件的环境(two Environment )分支条。    2.增加一个求解组合分支条( Solution Combination branch)，并