ansys疲劳分析基本方法

疲劳是指结构在低于静态极限强度载荷的重复载荷\猀断裂破坏的一根能够承受300KN拉力，在200KN循环载荷作用下，经历1,000,000次循环后亦会破坏。导致疲劳破坏的主要因素如下：

载荷的循环次数；

每一个循环的应力幅；

每一个循环的平均应力；

存在局部应力集中现象。

真正的疲劳计算要考虑所有这些因素，因为在预测其生命周期时，它计算“消耗”的某个部件是如何形成的。

3.1.1ANSYS程序处理疲劳问题的过程

ANSYS疲劳计算以ASME锅炉和压力容器规范(ASMEBoilerandPressureVesselCode)第三节(和第八节第二部分)作为计算的依据，采用简化了的弹塑性假设和Mimer累积疲劳准则。

除了根据ASME规范所建立的规则进行疲劳计算外，用户也可编写自己的宏指令，或选用合适的，利用ANSYS计算的结果进行疲劳计算。《ANSYSAPDLProgrammer‘sGuide》讨论了上述二种功能。

ANSYS程序的疲劳计算能力如下：

对现有的应力结果进行后处理，以确定体单元或壳单元模型的疲劳寿命耗用系数(fatigueusagefactors)(用于疲劳计算的线单元模型的应力必须人工输入)；

可以在一系列预先选定的位置上，确定一定数目的事件及组成这些事件„吀把这些位置上的应力储存起来；

可以在每一个位置上定义应力集中系数和给每一个事件定义比例系数。

3.1.2基本术语

位置tion)：在模型上储存疲劳应力的节点。这些节点是结构上某些容易产生疲劳破坏的位置。

事件(Event)：是在特定的应力循环过程中，在不同时刻的一系列应力状态，见本章§3.2.3.4。

载荷(Loading)：是事件的一部分，是其中一个应力状态。

应力幅：两个载荷之间应力状态之差的度量。应力平均值对结果的影响。

3.2疲劳计算

完成了应力计算后，就可以在通用后处理器POST1中进行疲劳计算。一般有五个主要步骤：

1、进入后处理POST1，恢复数据库；

2、建立位置、事件和载荷的数目，定义材料疲劳性质，确定应力位置和定义应力集中系数；

3、存储不同事件和不同载荷下关心位置的应力，并指定事件的重复次数和比例系数；

4、激活疲劳计算；

5、查看结果。

3.2.1进入POST1和恢复数据库

依照下列步骤进行疲劳计算：

POST1

命令：POST1

GUI：MainMenu>GeneralPostproc

2、把数据库文件(Jobname.DB)读入到在内存中(如果所要做的疲劳计算是正在进行的ANSYS计算过程的继续,则Jobname.DB文件已在内存中)。结果文件(Jobname.RST)必须已经存在并将其读入内存。

命令：RESUME

GUI：Utility>File>Resumefrom

3.2.2建立疲劳计算的规模、材料疲劳性质和疲劳计算的位置

定义下列p

位置、事件和载荷的最大数目；

材料的疲劳性质；

应力位置与应力集中系数(SCFs)。

1、定义位置、事件和载荷的最大数目

缺省情况下，疲劳计算最多包括5个节点位置，10个事件，每个事件中3个载荷。如果需要，可以通过下面的命令来设置较大的规的位置荷)。

命令：FTSIZE

GUI：MainMenu>GeneralPostproc>Fatigue>SizeSettings

2、定义材料的疲劳性质

为了计算各种耗用系数，以及为了包含简化弹塑性效应，必须定义材料的疲劳性质。在疲劳计算中，感兴趣的材料性质有：

S-N曲线：应力幅[(Smax-Smin)/2]-疲劳循环次数的关系曲线。ASMES-N曲线考虑了最大平均应力的影响。如果需要，应把输入的S-N曲线进行调节以便考虑平均应力强度效应。如果不输入S-N曲线，那么对于各种可能的应力状态的组合，应力幅将降序排列，但不计算耗用系数。

命令：FP

GUI：MainMenu>GeneralPostproc>Fatigue>S-NTable

Sm-T曲线：设计应力强度值-温度曲线。如要考虑检查应力范围是否进入塑性，就必须定义该曲线。

命令：FP

GUI：MainMenu>GeneralPostproc>Fatigue>Sm_TTable

弹塑性材料参数M和N(应变强化

圀需要使用简化的弹塑性准则时，才输入M、N。这些参数可以从ASME规范中获得。

命令：FP

GUI：MainMenu>GeneralPostproc>Fatigue>Elas-plasPar

下述例子说明了用于输入疲劳材料性质的FP命令的使用方法：

!DefinetheS-Ntable:

FP,1,10,30,100,300,1000,10000!AllowableCycles,N

FP,7,100000,1000000!"

FP,21,650,390,240,161,109,59!AlternatingStress-

FP,27,37,26!IntensityRange,S,ksi

!DefinetheSm-Ttable:

FP,41,100,200,300,400,500,600!Temperature,degF

FP,47,650,700,750,800!"

FP,51,20,20,20,18.7,17.4,16.4!"DesignStress-Intensity

FP,57,16.1,15.9.15.5,15.1!Value",Sm(=2/3*Syor

!1/3*Su),ksi

!Definetheelastic-plasticmaterialparameters:

FP,61,1.7,.3!MandN

3、定义应力位置和应力集中系数

下面的选项允许用户显式地定义疲劳计算中关心的节点位置、该位置的应力集中系数，以及在该位置的一个短的标题(可用20个字母)。

命令：FL

GUI：MainMenu>GlPostproc>Fatigue>StressLocations

注意--不是所有的疲劳分析都需要使用FL命令。如果使用FS、FSNODE或FSSECT等命令(见下)，则疲劳节点位置是自动定义的。假使在建模时包含有足够细的网格，则所计算的应力是准确的，因此不必指定应力集中系数SCFs(但是如果考虑表面影响、尺寸影响和腐蚀影响，则仍然需要指定SCFs)。在计算疲劳时如果只需要考察一个位置，则可以省略标题。假使定位明确，或是不需要应力集中系数和标题，则可完全不使用FL命令。

这里给出了在一个圆柱筒分析中应用FL命令的例子。整体坐标Y轴为旋转轴。在不同璧厚交界处的外壁给出了应力集中系数SCFs(针对轴向线性化应力)。

FL,1,281,,,,Line1atinside

FL,2,285,,1.85,,Line1atoutside

FL,3,311,,,,Line2atinside

FL,4,315,,2.11,,Line2atoutside图3-1圆柱筒应力集中系数SCFs3.2.3储存应力、指定事件循环次数和比例因子

3.2.3.1储存应力

为了进行疲劳计算，程序必须知道每一个位置上不同事件和载荷时的应力，以及每一个事件的循环次数。可采用下列选项来存储每一个位置、事件和载荷组合情况下的应›

人工储存应力；

从(Jobname.RBT)文件中取得节点应力；

横截面应力。

[警告]：程序从不假定存在0应力条件。如果一定要考虑零应力条件，就必须在每一个事件中明确地输入何处产生零应力。

下面的命令序列，说明如何存储应力。在某些情况下，用户也可以用LCASE命令代替SET命令。

人工储存应力：FS

从Jobname.RST中取出节点应力：SET，FSNODE

横截bATH，PPATH，SET，FSSECT

(横截面计算也需从Jobname.RST的数据中取得)

可以用不同的方法在一个事件中储存应力。下面说明各种不同的方法。

3.2.3.1.1人工储存应力

可以人工存储应力和温度(不是直接从Jobname.RST的结果文件取得)。在这种情况下，实际并没有将POST1的疲劳模块作为后处理器，而是仅仅作为疲劳计算器使用。线单元(如梁单元)的应力必须人工输入，因为疲劳模块不能如体元或壳元那样从结果文件中读取数据。

命令：FS

GUI：MainMenu>GeneralPostproc>Fatigue>-StoreStresses-

SpecifiedVal

下面的例子说明用上述命令输入的方式：

FS,201,1,2,1,-2.0,21.6,15.2,4.5,0.0,0.0

FS,201,1,2,7,450.3

在这个例子中，只输入了全应力(1～6项)和温度。如果还要输入线性化应力，它们将紧随在温度后面3项。

注意--对只有轴向应力的梁单元，只需输入一个应力分量(SX)，其余各项空白。

3.2.3.1.2从结果文件中提取应力

该方法把包含有6个分量的节点应力向量直接储存在结果的数据库内。随后可以用FS命令修正存入的应力分量。

注意--在执行FSNODE命令之前，必须使用SET命令，可能还有SHELL命令。SET库的Jobname.RST文件中读取某一特殊载荷子步下的结果，SHELL命令可选择从壳单元的顶面、中面或底面读取结果(缺省是从顶面读取结果)。

命令：FSNODE

GUI：MainMenu>GeneralPostproc>Fatigue>-StoreStresses-

FromrstFile

下面给出用FSNODE命令在一个事件的一个节点位置输入应力结果的例子：

SET,1!Definedatasetforloadstep

FSNODE,123,1,1!Stressvectoratnode123assignedtoevent1,

!loading1.

SET,2!Definedatasetforloadstep2

FSNODE,123,1,2!...event1,loading2

SET,3!...loep3

FSNODE,123,1,3!..t1,loading3图3-2一个事件中的三个荷载

3.2.3.1.3横截面应力

本选项计算和存储截面路径(它是由以前的PATH和PPATH命令定义的)端点的线性化应力。因为通常线性化应力计算是在能代表两个表面的最短距离的线段上进行的，因此，只需在两个表面上各取一个点来描述PPATH命令中的路径。这一步骤将从计算结果的数据库中获得应力在SET命令之前使用FSSECT命令。用FSSECT命令储量可用FS命令修正。

命令：FSSECT

GUI：MainMenu>GeneralPostproc>Fatigue>-StoreStresses-

AtCrossSect

使用FSSECT命令的方法可用下述例子说明。如果节点位置没有在FL命令中指定，则在该例子中命令FSSECT将自动给两个路径节点391和395设定位置编号。见图3-3。图3-3在执行FSSECT命令之前由PPATH命令识别的表面节点3.2.3.2列表、显示或删除储存的应力

用下列选项对储存的应力进行列表、绘图或删除。

1、列出每一个位置、每一个事件、每一种载荷或每一种应力状态下的储存应力：

命令：FSLIST

GUI：MainMenu>GeneralPoc>Fatigue>-StoreStresses-

ListStresses

2、对某一位置和事件，以载荷号的函数来>＀

命令：FSPLOT

GUI：MainMenu>GeneralPostproc>Fatigue>-StoreStresses-

PlotStresses

3、删除储存、事件和载荷下的应力状态：

命令：FSDELE

GUI：MainMenu>GeneralPostproc>Fatigue>-StoreStresses-

DeleStresses

4、删除某一位置上的所有应力：

命令：FL

GUI：MainMenu>GeneralPostproc>Fatigue>StressLocations

5、删除在某一事件中各种载荷下的所有应力：

命令：FE

GUI：MainMenu>GeneralPostproc>Fatigue>EraseEventData

3.2.3.3设定事件重复次数和比例系数

本选项确定事件重复次数(对该事件的所有荷载和所有位置)。也可对该事件中构成载荷的所有应力施加

命令：FE

GUI：MainMenu>GeneralPostproc>Fatigue>AssignEvents

下面给出该步骤的一个例子：

FE,1,-1!Eraseallparametersandfatiguestressesformerly

!usedforevent1.

FE,2,100,1.2！Assign100occurrencestoevents2,3and4,

*REPEAT,3,1!andscaleby1.2.

FE,5,500!Assi0occurrencestoevent5.3.2.3.4获得准确的耗用系数

结构常常承受各种最大和最小应力，它们发生的顺序是未知的(甚至是随机的)。因此就必须小心地考虑如何在各种可能的应力范围内，得到正确的重复循环次数，以获得有效的疲劳寿命耗用系数。

ANSYS程序自动计—瘀应力范围，同时采用我们通常所知道的“雨流”法记数，跟踪这些应力发生的次数。在选定的节点有事件进行搜索，以寻找产生最大应力幅的载荷对(应力矢量)。记录这些应力幅的重复次数，同时包含这些载荷的事件的剩余重复次数随之减少。最终至少有一个事件在某一位置被“用光”，而属于这一事件的其它应力状态，在随后的过程中将被忽略。这一过程一直进行着，直到所有的应力幅及重复次数都被计及后结束。

[警告]通常容易误用疲劳模块的记数功能。如果需要使疲劳计算得到正确的耗用系数，就必须小心地设计好事件。

设计事件时遵循下列规则：

1、了解ANSYS的内部逻辑记数算法，参见《ANSYSTheoryReference》§19.5，详细地讲述ANSYS程序如何记数。

2、由于在三维应力状态，很难预测哪一个载荷步具有极值应力，因此可以对每一个事件采用多个载荷群，以便成功获得极值应力(见图3-2)。

3、如果在某一给定的事件中，只包含一个极值应力，则将获得一个较保守的结果。如果在一个事件中，引入不止一个极值应力，则有时将产生不保守的结果。如下例所示：

考察一个加载历程，由两个稍有不同的循环构成：

载荷循环1：做500次循环，Sx=+50.0～-50.1ksi。

载荷循环2：做1000次循环，Sx=+50.1～-50.0ksi。

显然这是应力幅为50Ksi的1500次循环。然而，如果不小心把这些载荷划分成两个事件，则将导致不正确的记数结果。请看该结果如何：

事件1载荷1：Sx=50.0500次循环

茀2：Sx=-50.1

事件2载荷1：Sx=50.11000次循环

载荷2：Sx=-50.0

可能的应力幅是：

a、从E1,L1到E1,L2:50.05ksi

b、从E1,L1到E2,L1:0.05ksi

c、从E1,L1到E2,L2:50.00ksid、从E1,L2到E2,L1:50.10ksi

e、从E1,L2到E2,L2:0.05ksi

f、从E2,L1到E2,L2:50.05ksi

把这些应力幅从大到小排列：

d、从E1,L2到E2,L1:50.10ksi

a、从E1,L1到E1,L2:50.0

f、从E2,L1到E2,L2:50.05ksi

c、从E1,L1到E2,L2:50.00ksi

b、从E1,L1到E2,L1:0.05ksi

e、从E1,L2到E2,L2:0.05ksi

则记数将如此：

d、从E1,L2到E2,L1500次循环-E1和E2用去500次循环

a、从E1,L1到E1,L20次循环-E1已用光

f、从E2,L1到E2,L2500次循环-E2又用去500次循环

c、从E1,L1到E2,L20次循环-两个事件都用光

b、从E1,L1到E2,L10次循环-两个事件都用光

e、从E1,L2到E2,L20次循环-两个事件都用光

上面记录了50ksi应力幅的1000次循环，而不是已知的1500次循环。这个错误的结果是由于不恰当的事件组合导致„

假使载荷被描述为各自分离的事件(如E1,L1=>E1;E1,L2=>E2;E2,L1=>E3;E2,L2=>E4)，则将得到下列的记数范围：

d、从E2到E3的500次循环，用去500次循环

a、从E1到E2的0次循环，因为E2“用光”了

f、从E3到E4的500次循环，E3又用去500次循环，E4用去500次循环

c、从E1到E4的500次循环，E4又用去500次循环

d、从E1到E3的0次循环，因为E3“用光”了

e、从E2到E4的0次循环，因为E2和E4都“用光”了。

上述情况恰当地记录了累积疲劳损伤被，即50ksi应力幅的1500次循环。

4、相反，对每一个最大和最小应力条件应用分离的事件，则将使记录变得太保守。在这种情况下，应小心地选择应被一起记数的载荷，将它们划分为同一事件。下面的示例说明一些事件如何才能够包括多个极值应力条件。

考虑由两个载荷循环组成的载荷历程：

载荷循环1：从Sx=+100.1到+100.0ksi，做500次循环

载荷循环2：从Sx=+50.50.0ksi，做1000次循环

显然这些循环的最不利组合为以大约25ksi的应力幅做500次循环。

在本例中，如果将载荷划分为两个事件，则产生了25ksi范围内进行500次循环的正确记录。如每一个载荷作为一个独立的事件，将产生一种过于保守的记录，形成25ksi应力幅的1000次循环。

3.2.4激活疲劳计算

现在有了位置、应力、事件及所有指定的材料参数，可在指定位置执行疲劳计算。位置的确定可用节点本身，也可用位置编号。

命令：FTCALC

GUI：MainMenu>Generalpostproc>Fe>CalculateFatig3.2.5查看计算结果

疲劳计算结果被输出在输出窗口。如果你将输出转向[/OUTPUT]到文件中(如Jobname.OUT)，就可打开该文件查看结果。

命令：*LIST

GUI：UtilityMenu>List>Files>Other>Jobname.OUT

如果用户已输入S-N曲线，则所有部位的应力幅(从大到小排列)与相应的事件/载荷、使用循环次数、允许循环次数、温度和阶段耗用系数一起以表格形式输出。随后输出累积寿命耗用系数。

正如前面提到的，FTCALC输出时，将对任一给定的应力幅显示相应的事件和载荷，这些信息能帮助分离出那些引起最大的疲劳损伤的事件与载荷。

修正事件的一个方便的方法是把所有的疲劳数据写入Jobname.FATG文件中(它可以在执行FTCALC命令前或后做)。写入Jobname.FATG文件的数据是疲劳模块的命令流。可通过编name.FATG文件(文本文件)来修改事件，然后用/INPUT命令重新读入被修改的疲劳命令。

命令：FTWRITE

GUI：MainMenu>Generalpostproc>Fatigue>WriteFatigData.

3.2.6其它记数方法

前面曾介绍过“雨流”记数法。在时间-载荷历程未知的情况下，这一技术才有用处。如果时间-载荷历程已知，对每一个接连的事件分别做疲劳分析[FTCALC]，然后人工相加这些耗用系数，则可避免得到不恰当的保守结果。

3.2.7疲劳分析示例(}

下面是一个疲劳计算输入命令流：

!EnterPOST1andResumetheDatabase:

/POST1

RESUME,...

!NumberofLocations,Events,andLoadings

FTSIZE,...

!MaterialFatigueProperties:

FP,1,....!Nvalues

FP,21,...!Svalues

FP,41,...!Tvalues

FP,51,...!Smvalues

FP,61,...!Elastic-plasticmaterialparameters

!Locations,StressConcentrationFactors,andLocationTitlesFL,...

!StoreStresses(3DifferentMethods)

!StoreStressesManually:

FS,...

!RetrieveStressesfromtheResultsFile:

SET,...

FSNODE,...

!StoreStressesataCross-Section:

PPATH,...

SET,...

FSSECT,...

!EventRepetitionsandScaleFactors

FE,...

!ActivatetheFatigueCalculations

FTCALC,...

!RetheResults(Listtheoutputfile)

FINISH

其中的FTSIZE，FP，FL，FS，FSNODE，PPATH，FSSECT，FE和FTCALC考《ANSYSCommandsReference》。ANSYS疲劳计算以ASME锅炉和压力容器规范(ASMEBoilerandPressureVesselCode)第三节(和第八节第二部分)作为计算的依据，采用简化了的弹塑性假设和Mimer累积疲劳准则。

除了根据ASME规范所建立的规则进行疲劳计算外，用户也可编写自己的宏指令，或选用合适的第三方程序，利用ANSYS计算的结果进行疲劳计算。《ANSYSAPDLProgrammer‘sGuide》讨论了上述二种功能。

ANS瘀疲劳计算能力如下：

对现有的应力结果进行后处理，以确定体单元或壳单元模型的疲劳寿命耗用系数(fatigueusagefactors)(用于疲劳计算的线单元模型的应力必须人工输入)；

可以在一系列预先选定的位置上，确定一定数目的事件及组成这些事件的载荷，然后把这些位置上的应力储存起来；

可以在每一个位置上定义应力集中系数事件定义比例系数。

3.1.2

基本术语

位置(Location)：在模型上储存疲劳应力的节点。这些节点是结构上某些容易产生疲劳破坏的位置。

事件(Event)：是在特定的应力循环过程中，在不同时刻的一系列应力状态，见本章§3.2.3.4。

载荷(Loading)：是事件的一部分，是其中一个应力状态。

应力幅：两个载荷之间应力状态之差的度量。程序不考虑应力平均值对结果的影响。3.2

疲劳计算

完成了应力计算后，就可以在通用后处理器POST1中进行疲劳计算。一般有五个主要步骤：

1、进入后处理POST1，恢复数据库；

2、建立位置、事件和载荷的数目，定义材料疲劳性质，确定应力位置和定义应力集中系数；

3、存储不同事件和不同载荷下关心位置的应力，并指定事件的重复次数和比例系数；

4、激活疲劳计算；

5、查看结果。

3.2.1

进入POST1和恢复数据库

依照下列步骤进行疲劳计算：

1、进入POST1

命令：POST1

GUIMenu>GeneralPostproc

2、把数据库文件(Jobname.DB)读入到在内存中(如果所要做的疲劳计算是正在进行的ANSYS计算过程的继续,则Jobname.DB文件已在内存中)。结果文件(Jobname.RST)必须已经存在并将其读入内存。

命令：RESUME

GUI：UtilityMenu>File>Resumefrom3.2.2

建立疲劳计算的规模、材料疲劳性质和疲劳计算的位置

定义下列数据：

位置、事件和载荷的最大数目；

材料的疲劳性质；

应力位置与应力集中系数(SCFs)。

1、定义位置、事件和载荷的最大数目

缺省情况下，疲劳计算最多包括5个节点位置，10个事件，每个事件中3个载荷。如果需要，可以通过下面的命令来设置较大的规模(即较多的位置、事件和载荷)。

命令：FTSIZE

GUI：MainMenu>GeneralPostproc>Fatigue>SizeSettings

2、定义材料的疲劳性质

吀种耗用系数，以及为了包含简化弹塑性效应，必须定义材料的疲劳性质。在疲劳计算中，感兴趣的材料性质有：

S-N曲线：应力幅-疲劳循环次数的关系曲线。ASMES-N曲线考虑了最大平均应力的影响。如果需要，应把输入的S-N曲线进行调节以便考虑平均应力强度效应。如果不输入S-N曲线，那么对于各种