利用ANSYS随机振动分析功能实现随机疲劳分析.

合集下载

Ansys nCode疲劳分析

Ansys nCode疲劳分析
• 拉伸平均应力减少疲劳寿命 • 压缩平均应力增加疲劳寿命
Smax Sa Sm Sa Smin
© 2015 ANSYS, Inc. 36
DS
3.3 平均应力对疲劳寿命的影响
对于应力-疲劳寿命,平均应力的影响通常使用Goodman 或者 Gerber平均应力 修正理论 Goodman 平均应力修正理论假设一条直线
© 2015 ANSYS, Inc.
34
3.2 缺口应力集中对疲劳寿命的影响
应力梯度影响可以通过Advanced Edit进行定义
© 2015 ANSYS, Inc.
35
3 平均应力对疲劳寿命的影响
循环峰值应力范围(总应力变化)是影响疲劳寿命的主要因素
– 在应力循环中平均应力是第二个影响因素
在nCode中采用FKM的应力梯度法
– 确定校正系数
在雨流计数前得到有效应力分量:
© 2015 ANSYS, Inc. 33
3.2 缺口应力集中对疲劳寿命的影响
− User:从查找表中用户自定义应力梯度校正 自定义格式:
Stress gradient correction factor file v1.0 Dimension=mm #Normalized stress gradient G, correction factor n StartCorrectionData 0,1 1,1.1 10,1.2 100,1.3 EndCorrectionData
37
3.4 平均应力对疲劳寿命的影响
Sm Sa 1 Se Su
• Goodman 更保守 • 实际的数据往往介于Goodman 和Gerber之间
© 2015 ANSYS, Inc.

ansys疲劳分析基本方法

ansys疲劳分析基本方法

疲劳是指结构在低于静态极限强度载荷的重复载荷作用下,出现断裂破坏的现象;例如一根能够承受 300 KN 拉力作用的钢杆,在 200 KN 循环载荷作用下,经历 1,000,000 次循环后亦会破坏;导致疲劳破坏的主要因素如下:载荷的循环次数;每一个循环的应力幅;每一个循环的平均应力;存在局部应力集中现象;真正的疲劳计算要考虑所有这些因素,因为在预测其生命周期时,它计算“消耗”的某个部件是如何形成的;ANSYS程序处理疲劳问题的过程ANSYS 疲劳计算以ASME锅炉和压力容器规范ASME Boiler and Pressure Vessel Code第三节和第八节第二部分作为计算的依据,采用简化了的弹塑性假设和Mimer累积疲劳准则;除了根据 ASME 规范所建立的规则进行疲劳计算外,用户也可编写自己的宏指令,或选用合适的第三方程序,利用 ANSYS 计算的结果进行疲劳计算;ANSYS APDL Programmer‘s Guide讨论了上述二种功能;ANSYS程序的疲劳计算能力如下:对现有的应力结果进行后处理,以确定体单元或壳单元模型的疲劳寿命耗用系数fatigue usage factors用于疲劳计算的线单元模型的应力必须人工输入;可以在一系列预先选定的位置上,确定一定数目的事件及组成这些事件的载荷,然后把这些位置上的应力储存起来;可以在每一个位置上定义应力集中系数和给每一个事件定义比例系数;基本术语位置Location:在模型上储存疲劳应力的节点;这些节点是结构上某些容易产生疲劳破坏的位置;事件Event:是在特定的应力循环过程中,在不同时刻的一系列应力状态,见本章§;载荷Loading:是事件的一部分,是其中一个应力状态;应力幅:两个载荷之间应力状态之差的度量;程序不考虑应力平均值对结果的影响;疲劳计算完成了应力计算后,就可以在通用后处理器 POST1 中进行疲劳计算;一般有五个主要步骤:1、进入后处理 POST1,恢复数据库;2、建立位置、事件和载荷的数目,定义材料疲劳性质,确定应力位置和定义应力集中系数;3、存储不同事件和不同载荷下关心位置的应力,并指定事件的重复次数和比例系数;4、激活疲劳计算;5、查看结果;进入POST1和恢复数据库依照下列步骤进行疲劳计算:1、进入POST1命令:POST1GUI:Main Menu>General Postproc2、把数据库文件读入到在内存中如果所要做的疲劳计算是正在进行的 ANSYS 计算过程的继续,则文件已在内存中;结果文件必须已经存在并将其读入内存; 命令:RESUMEGUI:Utility Menu>File>Resume from建立疲劳计算的规模、材料疲劳性质和疲劳计算的位置定义下列数据:位置、事件和载荷的最大数目;材料的疲劳性质;应力位置与应力集中系数SCFs;1、定义位置、事件和载荷的最大数目缺省情况下,疲劳计算最多包括5个节点位置,10个事件,每个事件中3个载荷;如果需要,可以通过下面的命令来设置较大的规模即较多的位置、事件和载荷; 命令:FTSIZEGUI:Main Menu>General Postproc>Fatigue>Size Settings2、定义材料的疲劳性质为了计算各种耗用系数,以及为了包含简化弹塑性效应,必须定义材料的疲劳性质;在疲劳计算中,感兴趣的材料性质有:S-N 曲线:应力幅Smax -Smin/2-疲劳循环次数的关系曲线;ASME S-N 曲线考虑了最大平均应力的影响;如果需要,应把输入的 S-N 曲线进行调节以便考虑平均应力强度效应;如果不输入S-N曲线,那么对于各种可能的应力状态的组合,应力幅将降序排列,但不计算耗用系数;命令:FPGUI:Main Menu>General Postproc>Fatigue>S-N TableSm-T曲线:设计应力强度值-温度曲线;如要考虑检查应力范围是否进入塑性,就必须定义该曲线;命令:FPGUI:Main Menu>General Postproc>Fatigue>Sm_T Table弹塑性材料参数 M 和 N应变强化指数;只在需要使用简化的弹塑性准则时,才输入 M、N;这些参数可以从 ASME 规范中获得;命令:FPGUI:Main Menu>General Postproc>Fatigue>Elas-plas Par下述例子说明了用于输入疲劳材料性质的FP命令的使用方法:Define the S-N table:FP,1,10,30,100,300,1000,10000 Allowable Cycles, NFP,7,100000,1000000 "FP,21,650,390,240,161,109,59 Alternating Stress-FP,27,37,26 Intensity Range, S, ksiDefine the Sm-T table:FP,41,100,200,300,400,500,600 Temperature, deg FFP,47,650,700,750,800 "FP,51,20,20,20,,, "Design Stress-IntensityFP,57,,, Value", Sm =2/3Sy or1/3 Su, ksiDefine the elastic-plastic material parameters:FP,61,,.3 M and N3、定义应力位置和应力集中系数下面的选项允许用户显式地定义疲劳计算中关心的节点位置、该位置的应力集中系数,以及在该位置的一个短的标题可用20个字母;命令:FLGUI:Main Menu>General Postproc>Fatigue>Stress Locations注意--不是所有的疲劳分析都需要使用FL命令;如果使用FS、FSNODE或FSSECT等命令见下,则疲劳节点位置是自动定义的;假使在建模时包含有足够细的网格,则所计算的应力是准确的,因此不必指定应力集中系数 SCFs 但是如果考虑表面影响、尺寸影响和腐蚀影响,则仍然需要指定 SCFs ;在计算疲劳时如果只需要考察一个位置,则可以省略标题;假使定位明确,或是不需要应力集中系数和标题,则可完全不使用FL命令;这里给出了在一个圆柱筒分析中应用FL命令的例子;整体坐标 Y 轴为旋转轴;在不同璧厚交界处的外壁给出了应力集中系数 SCFs 针对轴向线性化应力; FL,1,281,,,,Line 1 at insideFL,2,285,,,,Line 1 at outsideFL,3,311,,,,Line 2 at insideFL,4,315,,,,Line 2 at outside图3-1 圆柱筒应力集中系数 SCFs储存应力、指定事件循环次数和比例因子储存应力为了进行疲劳计算,程序必须知道每一个位置上不同事件和载荷时的应力,以及每一个事件的循环次数;可采用下列选项来存储每一个位置、事件和载荷组合情况下的应力:人工储存应力;从文件中取得节点应力;横截面应力;警告:程序从不假定存在0应力条件;如果一定要考虑零应力条件,就必须在每一个事件中明确地输入何处产生零应力;下面的命令序列,说明如何存储应力;在某些情况下,用户也可以用LCASE命令代替SET命令;人工储存应力: FS从中取出节点应力: SET,FSNODE横截面应力: PATH,PPATH,SET,FSSECT横截面计算也需从的数据中取得可以用不同的方法在一个事件中储存应力;下面说明各种不同的方法;.1 人工储存应力可以人工存储应力和温度不是直接从的结果文件取得;在这种情况下,实际并没有将 POST1 的疲劳模块作为后处理器,而是仅仅作为疲劳计算器使用;线单元如梁单元的应力必须人工输入,因为疲劳模块不能如体元或壳元那样从结果文件中读取数据;命令:FSGUI:Main Menu>General Postproc>Fatigue>-Store Stresses-Specified Val下面的例子说明用上述命令输入的方式:FS,201,1,2,1,,,,,,FS,201,1,2,7,在这个例子中,只输入了全应力1~6项和温度;如果还要输入线性化应力,它们将紧随在温度后面,即8~13项;注意--对只有轴向应力的梁单元,只需输入一个应力分量SX,其余各项空白;.2 从结果文件中提取应力该方法把包含有6个分量的节点应力向量直接储存在结果的数据库内;随后可以用FS命令修正存入的应力分量;注意--在执行FSNODE命令之前,必须使用SET命令,可能还有SHELL命令;SET命令从数据库的文件中读取某一特殊载荷子步下的结果,SHELL命令可选择从壳单元的顶面、中面或底面读取结果缺省是从顶面读取结果;命令:FSNODEGUI:Main Menu>General Postproc>Fatigue>-Store Stresses-From rst File下面给出用FSNODE命令在一个事件的一个节点位置输入应力结果的例子: SET,1 Define data set for load step 1FSNODE,123,1,1 Stress vector at node 123 assigned to event 1,loading 1.SET,2 Define data set for load step 2FSNODE,123,1,2 ...event 1, loading 2SET,3 ...load step 3FSNODE,123,1,3 ...event 1, loading 3图3-2 一个事件中的三个荷载.3 横截面应力本选项计算和存储截面路径它是由以前的PATH和PPATH命令定义的端点的线性化应力;因为通常线性化应力计算是在能代表两个表面的最短距离的线段上进行的,因此,只需在两个表面上各取一个点来描述PPATH命令中的路径;这一步骤将从计算结果的数据库中获得应力;因此必须在SET命令之前使用FSSECT 命令;用FSSECT命令储存的应力分量可用FS 命令修正;命令:FSSECTGUI:Main Menu>General Postproc>Fatigue> -Store Stresses -At Cross Sect使用FSSECT命令的方法可用下述例子说明;如果节点位置没有在FL命令中指定,则在该例子中命令FSSECT将自动给两个路径节点 391 和 395 设定位置编号;见图3-3;图3-3 在执行 FSSECT 命令之前由 PPATH 命令识别的表面节点列表、显示或删除储存的应力用下列选项对储存的应力进行列表、绘图或删除;1、列出每一个位置、每一个事件、每一种载荷或每一种应力状态下的储存应力:命令:FSLISTGUI:Main Menu>General Postproc>Fatigue>-Store Stresses- List Stresses2、对某一位置和事件,以载荷号的函数来显示应力项:命令:FSPLOTGUI:Main Menu>General Postproc>Fatigue>-Store Stresses- Plot Stresses3、删除储存在某一位置、事件和载荷下的应力状态:命令:FSDELEGUI:Main Menu>General Postproc>Fatigue>-Store Stresses- Dele Stresses4、删除某一位置上的所有应力:命令:FLGUI:Main Menu>General Postproc>Fatigue>Stress Locations 5、删除在某一事件中各种载荷下的所有应力:命令:FEGUI:Main Menu>General Postproc>Fatigue>Erase Event Data设定事件重复次数和比例系数本选项确定事件重复次数对该事件的所有荷载和所有位置;也可对该事件中构成载荷的所有应力施加比例系数;命令:FEGUI:Main Menu>General Postproc>Fatigue>Assign Events下面给出该步骤的一个例子:FE,1,-1 Erase all parameters and fatigue stresses formerlyused for event 1.FE,2,100, Assign 100 occurrences to events 2, 3 and 4,REPEAT,3,1 and scale by .FE,5,500 Assign 500 occurrences to event 5.获得准确的耗用系数结构常常承受各种最大和最小应力,它们发生的顺序是未知的甚至是随机的;因此就必须小心地考虑如何在各种可能的应力范围内,得到正确的重复循环次数,以获得有效的疲劳寿命耗用系数;ANSYS 程序自动计算所有可能的应力范围,同时采用我们通常所知道的“雨流”法记数,跟踪这些应力发生的次数;在选定的节点位置,对所有事件进行搜索,以寻找产生最大应力幅的载荷对应力矢量;记录这些应力幅的重复次数,同时包含这些载荷的事件的剩余重复次数随之减少;最终至少有一个事件在某一位置被“用光”,而属于这一事件的其它应力状态,在随后的过程中将被忽略;这一过程一直进行着,直到所有的应力幅及重复次数都被计及后结束;警告通常容易误用疲劳模块的记数功能;如果需要使疲劳计算得到正确的耗用系数,就必须小心地设计好事件;设计事件时遵循下列规则:1、了解 ANSYS 的内部逻辑记数算法,参见ANSYS Theory Reference§,详细地讲述 ANSYS 程序如何记数;2、由于在三维应力状态,很难预测哪一个载荷步具有极值应力,因此可以对每一个事件采用多个载荷群,以便成功获得极值应力见图3-2;3、如果在某一给定的事件中,只包含一个极值应力,则将获得一个较保守的结果;如果在一个事件中,引入不止一个极值应力,则有时将产生不保守的结果;如下例所示:考察一个加载历程,由两个稍有不同的循环构成:载荷循环1:做 500次循环,Sx = + ~ ksi;载荷循环2:做1000次循环,Sx = + ~ ksi;显然这是应力幅为 50 Ksi的 1500 次循环;然而,如果不小心把这些载荷划分成两个事件,则将导致不正确的记数结果;请看该结果如何:事件1 载荷1: Sx = 500次循环载荷2: Sx =事件2 载荷1: Sx = 1000次循环载荷2: Sx =可能的应力幅是:a、从 E1,L1 到 E1,L2: ksib、从 E1,L1 到 E2,L1: ksic、从 E1,L1 到 E2,L2: ksid、从 E1,L2 到 E2,L1: ksie、从 E1,L2 到 E2,L2: ksif、从 E2,L1 到 E2,L2: ksi把这些应力幅从大到小排列:d、从 E1,L2 到 E2,L1: ksia、从 E1,L1 到 E1,L2: ksif、从 E2,L1 到 E2,L2: ksic、从 E1,L1 到 E2,L2: ksib、从 E1,L1 到 E2,L1: ksie、从 E1,L2 到 E2,L2: ksi则记数将如此:d、从 E1,L2 到 E2,L1 500 次循环 - E1 和 E2 用去500 次循环a、从 E1,L1 到 E1,L2 0 次循环 - E1 已用光f、从 E2,L1 到 E2,L2 500 次循环 - E2 又用去 500 次循环c、从 E1,L1 到 E2,L2 0 次循环 - 两个事件都用光b、从 E1,L1 到 E2,L1 0 次循环 - 两个事件都用光e、从 E1,L2 到 E2,L2 0 次循环 - 两个事件都用光上面记录了 50 ksi 应力幅的 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 都“用光”了;上述情况恰当地记录了累积疲劳损伤被,即 50 ksi 应力幅的 1500 次循环;4、相反,对每一个最大和最小应力条件应用分离的事件,则将使记录变得太保守;在这种情况下,应小心地选择应被一起记数的载荷,将它们划分为同一事件;下面的示例说明一些事件如何才能够包括多个极值应力条件;考虑由两个载荷循环组成的载荷历程:载荷循环1:从 Sx = + 到 + ksi,做 500 次循环载荷循环2:从 Sx = + 到 + ksi,做 1000 次循环显然这些循环的最不利组合为以大约 25 ksi 的应力幅做 500 次循环;在本例中,如果将载荷划分为两个事件,则产生了 25 ksi 范围内进行 500 次循环的正确记录;如每一个载荷作为一个独立的事件,将产生一种过于保守的记录,形成 25 ksi 应力幅的 1000 次循环;激活疲劳计算现在有了位置、应力、事件及所有指定的材料参数,可在指定位置执行疲劳计算;位置的确定可用节点本身,也可用位置编号;命令:FTCALCGUI:Main Menu>General postproc>Fatigue>Calculate Fatig查看计算结果疲劳计算结果被输出在输出窗口;如果你将输出转向/OUTPUT到文件中如,就可打开该文件查看结果;命令:LISTGUI:Utility Menu>List>Files>Other>如果用户已输入 S-N 曲线,则所有部位的应力幅从大到小排列与相应的事件/载荷、使用循环次数、允许循环次数、温度和阶段耗用系数一起以表格形式输出;随后输出累积寿命耗用系数;正如前面提到的,FTCALC输出时,将对任一给定的应力幅显示相应的事件和载荷,这些信息能帮助分离出那些引起最大的疲劳损伤的事件与载荷;修正事件的一个方便的方法是把所有的疲劳数据写入文件中它可以在执行FTCALC命令前或后做;写入文件的数据是疲劳模块的命令流;可通过编辑文件文本文件来修改事件,然后用/INPUT命令重新读入被修改的疲劳命令;命令:FTWRITEGUI:Main Menu>General postproc>Fatigue>Write Fatig Data.其它记数方法前面曾介绍过“雨流”记数法;在时间-载荷历程未知的情况下,这一技术才有用处;如果时间-载荷历程已知,对每一个接连的事件分别做疲劳分析FTCALC,然后人工相加这些耗用系数,则可避免得到不恰当的保守结果;疲劳分析示例命令流方法下面是一个疲劳计算输入命令流:Enter POST1 and Resume the Database:/POST1RESUME,...Number of Locations, Events, and LoadingsFTSIZE,...Material Fatigue Properties:FP,1,.... N valuesFP,21,... S valuesFP,41,... T valuesFP,51,... Sm valuesFP,61,... Elastic-plastic material parametersLocations, Stress Concentration Factors, and Location TitlesFL,...Store Stresses 3 Different MethodsStore Stresses Manually:FS,...Retrieve Stresses from the Results File: SET,...FSNODE,...Store Stresses at a Cross-Section:PPATH,...SET,...FSSECT,...Event Repetitions and Scale FactorsFE,...Activate the Fatigue CalculationsFTCALC,...Review the Results List the output fileFINISH其中的 FTSIZE,FP,FL,FS,FSNODE,PPATH,FSSECT,FE 和FTCALC 命令,请参考ANSYS Commands Reference;ANSYS 疲劳计算以ASME锅炉和压力容器规范ASME Boiler and Pressure Vesse l Code第三节和第八节第二部分作为计算的依据,采用简化了的弹塑性假设和M imer累积疲劳准则; 除了根据 ASME 规范所建立的规则进行疲劳计算外,用户也可编写自己的宏指令,或选用合适的第三方程序,利用 ANSYS 计算的结果进行疲劳计算;ANSYS APDL Programmer‘s Guide讨论了上述二种功能; ANSYS 程序的疲劳计算能力如下:对现有的应力结果进行后处理,以确定体单元或壳单元模型的疲劳寿命耗用系数fatigue usage factors用于疲劳计算的线单元模型的应力必须人工输入;可以在一系列预先选定的位置上,确定一定数目的事件及组成这些事件的载荷,然后把这些位置上的应力储存起来;可以在每一个位置上定义应力集中系数和给每一个事件定义比例系数; 基本术语位置Location:在模型上储存疲劳应力的节点;这些节点是结构上某些容易产生疲劳破坏的位置; 事件Event:是在特定的应力循环过程中,在不同时刻的一系列应力状态,见本章§; 载荷Loading:是事件的一部分,是其中一个应力状态; 应力幅:两个载荷之间应力状态之差的度量;程序不考虑应力平均值对结果的影响;疲劳计算完成了应力计算后,就可以在通用后处理器 POST1 中进行疲劳计算;一般有五个主要步骤: 1、进入后处理 POST1,恢复数据库; 2、建立位置、事件和载荷的数目,定义材料疲劳性质,确定应力位置和定义应力集中系数;3、存储不同事件和不同载荷下关心位置的应力,并指定事件的重复次数和比例系数;4、激活疲劳计算;5、查看结果; 进入POST1和恢复数据库依照下列步骤进行疲劳计算: 1、进入POST1 命令:POST1 GUI:Main Menu>General Postproc 2、把数据库文件读入到在内存中如果所要做的疲劳计算是正在进行的 ANSYS 计算过程的继续,则文件已在内存中;结果文件必须已经存在并将其读入内存; 命令:RESUME GUI:Utility Menu>File>Resume from 建立疲劳计算的规模、材料疲劳性质和疲劳计算的位置定义下列数据:位置、事件和载荷的最大数目;材料的疲劳性质;应力位置与应力集中系数SCFs; 1、定义位置、事件和载荷的最大数目缺省情况下,疲劳计算最多包括5个节点位置,10个事件,每个事件中3个载荷;如果需要,可以通过下面的命令来设置较大的规模即较多的位置、事件和载荷; 命令:FTSIZE GU I:Main Menu>General Postproc>Fatigue>Size Settings 2、定义材料的疲劳性质为了计算各种耗用系数,以及为了包含简化弹塑性效应,必须定义材料的疲劳性质;在疲劳计算中,感兴趣的材料性质有: S-N 曲线:应力幅-疲劳循环次数的关系曲线;ASME S-N 曲线考虑了最大平均应力的影响;如果需要,应把输入的 S-N 曲线进行调节以便考虑平均应力强度效应;如果不输入S-N曲线,那么对于各种可能的应力状态的组合,应力幅将降序排列,但不计算耗用系数; 命令:FP GUI:Main Menu>General Postproc>Fatigue>S-N Table Sm-T曲线:设计应力强度值-温度曲线;如要考虑检查应力范围是否进入塑性,就必须定义该曲线; 命令:FP GUI:Main Menu>General Postproc>Fatigue>Sm_T Ta ble 弹塑性材料参数 M 和 N应变强化指数;只在需要使用简化的弹塑性准则时,才输入 M、N;这些参数可以从 ASME 规范中获得; 命令:FP GUI:Main Menu>General Postproc>Fatigue>Elas-plas Par 下述例子说明了用于输入疲劳材料性质的FP命令的使用方法: Define the S-N table: FP,1,10,30,1 00,300,1000,10000 Allowable Cycles, N FP,7,100000,1000000 " FP,2 1,650,390,240,161,109,59 Alternating Stress- FP,27,37,26 Intensity Range, S, ksi Define the Sm-T table: FP,41,100,200,300,400,500,600 Temperature, deg F FP,47,650,700,750,800 " FP,51,20,20,20,,, " Design Stress-Intensity FP,57,,, Value", Sm =2/3Sy or 1/3 Su, ksi Define the elastic-plastic material parameters: FP,61,,.3 M and N 3、定义应力位置和应力集中系数下面的选项允许用户显式地定义疲劳计算中关心的节点位置、该位置的应力集中系数,以及在该位置的一个短的标题可用20个字母; 命令:FL GUI:Main Menu>General Postproc>Fatigue>Str ess Locations 注意--不是所有的疲劳分析都需要使用 FL 命令;如果使用 F S、FSNODE 或 FSSECT 等命令见下,则疲劳节点位置是自动定义的;假使在建模时包含有足够细的网格,则所计算的应力是准确的,因此不必指定应力集中系数 SC Fs 但是如果考虑表面影响、尺寸影响和腐蚀影响,则仍然需要指定 SCFs ;在计算疲劳时如果只需要考察一个位置,则可以省略标题;假使定位明确,或是不需要应力集中系数和标题,则可完全不使用 FL 命令; 这里给出了在一个圆柱筒分析中应用 FL 命令的例子;整体坐标 Y 轴为旋转轴;在不同璧厚交界处的外壁给出了应力集中系数 SCFs 针对轴向线性化应力;。

(完整版)疲劳分析的数值计算方法及ANSYS疲劳分析实例

(完整版)疲劳分析的数值计算方法及ANSYS疲劳分析实例

第十四章疲劳分析的数值计算方法及实例第一节引言零件或构件由于交变载荷的反复作用,在它所承受的交变应力尚未达到静强度设计的许用应力情况下就会在零件或构件的局部位置产生疲劳裂纹并扩展、最后突然断裂。

这种现象称为疲劳破坏。

疲劳裂纹的形成和扩展具有很大的隐蔽性而在疲劳断裂时又具有瞬发性,因此疲劳破坏往往会造成极大的经济损失和灾难性后果。

金属的疲劳破坏形式和机理不同与静载破坏,所以零件疲劳强度的设计计算不能为经典的静强度设计计算所替代,属于动强度设计。

随着机车车辆向高速、大功率和轻量化方向的迅速发展,其疲劳强度及其可靠性的要求也越来越高。

近几年随着我国铁路的不断提速,机车、车辆和道轨等铁路设施的疲劳断裂事故不断发生,越来越引起人们的重视。

疲劳强度设计及其研究正在成为我国高速机车车辆设计制造中的一项不可缺少的和重要的工作。

金属疲劳的研究已有近150年的历史,有相当多的学者和工程技术人员进行了大量的研究,得到了许多关于金属疲劳损伤和断裂的理论及有关经验技术。

但是由于疲劳破坏的影响因素多而复杂并且这些因素互相影响又与构件的实际情况密切相关,使得其应用性成果尚远远不能满足工程设计和生产应用的需要。

据统计,至今有约90%的机械零部件的断裂破坏仍然是由直接于疲劳或者间接疲劳而引起的。

因此,在21世纪的今天,尤其是在高速和大功率化的新产品的开发制造中,其疲劳强度或疲劳寿命的设计十分重要,并且往往需要同时进行相应的试验研究和试验验证。

疲劳断裂是因为在零件或构件表层上的高应力或强度比较低弱的部位区域产生疲劳裂纹,并进一步扩展而造成的。

这些危险部位小到几个毫米甚至几十个微米的范围,零件或构件的几何缺口根部、表面缺陷、切削刀痕、碰磕伤痕及材料的内部缺陷等往往是这种危险部位。

因此,提高构件疲劳强度的基本途径主要有两种。

一种是机械设计的方法,主要有优化或改善缺口形状,改进加工工艺工程和质量等手段将危险点的峰值应力降下来;另一种是材料冶金的方法,即用热处理手段将危险点局部区域的疲劳强度提高,或者是提高冶金质量来减少金属基体中的非金属夹杂等材料缺陷等局部薄弱区域。

ansys疲劳可靠性及随机振动分析

ansys疲劳可靠性及随机振动分析

ANSYS疲劳分析
5. 调入载荷步计算并储存节点应力分量
Main Menu>General Postproc>Read Results>First Set(Next Set) Main Menu>General Postproc>Fatigue>Store Stresses>From rst File
ANSYS疲劳分析
3. 定义材料应变强化指数M和N
Main Menu>General Postproc>Fatigue>Property Table>Elas-plas Par
ANSYS疲劳分析
4. 定义关键点及应力集中系数
Main Menu>General Postproc>Fatigue>Stress Locations
!定义载荷变量 !设置实常数 !定义材料参数
!建立平面 !建立中间圆孔面 !面的布尔操作 !设置单元划分的分段数 !定义单元为四边形 !网格自由划分 !划分网格 !结束前处理
ANSYS可靠性分析
二.进入概率分析模块,选定分析文件
Main Menu>Prob Design>Analysis File>Assign
ANSYS可靠性分析
② 显示输入随机变量的采样曲线
(Main Menu>Prob Results>Sampl History)
ANSYS可靠性分析
③ 显示输入随机变量的累计概率分布曲线
(Main Menu>Prob Results>CumulativeDF)
ANSYS可靠性分析
④ 还可以显示输出随机变量的概率敏感图以及输出随机变 量和输入随机变量的相关性等。

ANSYS疲劳分析的应用

ANSYS疲劳分析的应用

ANSYS疲劳分析的应用在传统的设计过程中,设计人员在概念或详细设计阶段通常使用简单而不真实的计算来估计产品的寿命,而对这些估计寿命的验证通常是通过一定量物理样机的耐久试验得到。

不但试验周期长、耗资大,而且许多参数与失效的定量关系也不能在试验中得出,试验结论还可能受许多偶然因素的影响。

因此对于产品疲劳寿命的仿真分析方法越来越受到产品设计人员的关注。

在塑料机械中,模板是注塑机最重要的零件之一,它的成本是注塑机成本的主要组成部分,模板断裂,注塑机就不能正常工作。

从强度出发,当然是选用高质量的材料,而且尽量将模板做得厚一些,但这两点均提高模板造价,影响整机成本。

目前模板大部分采用球墨铸铁铸造。

这主要考虑:(1) 在模板上铸出加强筋或将模板掏空,可有效减少质量;(2) 由于球铁较易于精铸(树脂砂铸造),使加工余量大大减少,可有效减少加工成本;(3) 球铁刚性较好,也具有一定强度。

虽然设计者充分考虑了模板的强度、刚度,但仍然有许多模板断裂的事故发生,其原因在于模板断裂不是因为静力破坏而是因为疲劳破坏。

一、元原理及模型建立当材料或结构受到多次重复变化的载荷作用后,在应力值虽然始终没有超过材料的强度极限,甚至比弹性极限还低的情况下就可能发生破坏。

这种在交变载荷作用下材料或结构的破坏现象称为疲劳破坏[1 ] 。

结构的疲劳破坏,首先在局部区域产生裂缝,一般是在零件和构件的表面,也可能在零件内部有缺陷处,即应力最高的区域。

由于该区域代表了整个结构的疲劳强度,所以该区域称为危险区,危险区的应力、应变变化情形为结构疲劳分析中所需的应力或应变2时间历程。

因此,结构疲劳应力分析的目的,就是要求得结构在承受各种负荷时,对其危险区的应力或应变响应,作为结构疲劳设计的依据之一。

在进行工程结构疲劳分析时,常应用ANSYS 软件为分析工具来确定结构的高应力危险区,并进行负荷谱转化为应力谱或应变谱的工作。

本工作将引用基本理论[2]:其中,式中,[ B ] 为应变矩阵;[ D] 为弹性矩阵;{ f e} 及[ Ke ] 为单元节点力及单元刚度矩阵。

ANSYSworkbench教程—疲劳分析

ANSYSworkbench教程—疲劳分析

第一章简介1.1 疲劳概述结构失效的一个常见原因是疲劳,其造成破坏与重复加载有关。

疲劳通常分为两类:高周疲劳是当载荷的循环(重复)次数高(如1e4 -1e9)的情况下产生的。

因此,应力通常比材料的极限强度低,应力疲劳(Stress-based)用于高周疲劳;低周疲劳是在循环次数相对较低时发生的。

塑性变形常常伴随低周疲劳,其阐明了短疲劳寿命。

一般认为应变疲劳(strain-based)应该用于低周疲劳计算。

在设计仿真中,疲劳模块拓展程序(Fatigue Module add-on)采用的是基于应力疲劳(stress-based)理论,它适用于高周疲劳。

接下来,我们将对基于应力疲劳理论的处理方法进行讨论。

1.2 恒定振幅载荷在前面曾提到,疲劳是由于重复加载引起:当最大和最小的应力水平恒定时,称为恒定振幅载荷,我们将针对这种最简单的形式,首先进行讨论。

否则,则称为变化振幅或非恒定振幅载荷。

1.3 成比例载荷载荷可以是比例载荷,也可以非比例载荷:比例载荷,是指主应力的比例是恒定的,并且主应力的削减不随时间变化,这实质意味着由于载荷的增加或反作用的造成的响应很容易得到计算。

相反,非比例载荷没有隐含各应力之间相互的关系,典型情况包括:σ1/σ2=constant在两个不同载荷工况间的交替变化;交变载荷叠加在静载荷上;非线性边界条件。

1.4 应力定义考虑在最大最小应力值σmin和σmax作用下的比例载荷、恒定振幅的情况:应力范围Δσ定义为(σmax-σmin)平均应力σm定义为(σmax+σmin)/2应力幅或交变应力σa是Δσ/2应力比R是σmin/σmax当施加的是大小相等且方向相反的载荷时,发生的是对称循环载荷。

这就是σm=0,R=-1的情况。

当施加载荷后又撤除该载荷,将发生脉动循环载荷。

这就是σm=σmax/2,R=0的情况。

1.5 应力-寿命曲线载荷与疲劳失效的关系,采用的是应力-寿命曲线或S-N曲线来表示:(1)若某一部件在承受循环载荷, 经过一定的循环次数后,该部件裂纹或破坏将会发展,而且有可能导致失效;(2)如果同个部件作用在更高的载荷下,导致失效的载荷循环次数将减少;(3)应力-寿命曲线或S-N曲线,展示出应力幅与失效循环次数的关系。

利用ANSYS随机振动分析功能实现随机疲劳分析

利用ANSYS随机振动分析功能实现随机疲劳分析

利用ANSYS随机振动分析功能实现随机疲劳分析ANSYS是一款常用的工程仿真软件,具有强大的分析功能。

在进行随机疲劳分析时,可以利用ANSYS的随机振动分析功能来模拟随机加载下的疲劳损伤。

随机疲劳分析是一种考虑工作载荷随机性对结构疲劳寿命影响的方法。

通过采用随机振动分析,可以考虑到工作载荷的随机特性,进一步分析结构的疲劳损伤。

ANSYS中的随机振动分析功能可以通过以下步骤来进行:1.几何建模:首先,需要进行结构的几何建模。

使用ANSYS的几何建模工具可以创建出要进行疲劳分析的结构。

2.材料属性定义:在进行材料属性的定义时,需要确定材料的弹性模量、泊松比、密度和疲劳参数等。

可以根据材料的材料数据手册来获取这些参数。

3.边界条件设置:在进行随机振动分析时,需要设置结构的边界条件。

这些边界条件可以是结构受到的随机外载荷或者是结构与其他部件的接触情况。

4.加载设置:在进行随机振动分析时,需要设置结构受到的随机载荷。

这些载荷可以是来自于实际工况的随机载荷,也可以通过振动台试验数据等手段获取。

5.随机振动分析:利用ANSYS的随机振动分析功能,可以进行频域分析或时域分析。

频域分析可以用于计算结构的响应功率谱密度,时域分析可以用于计算结构的随机响应。

6.疲劳寿命计算:在获得结构的随机响应后,可以进行疲劳寿命计算。

根据结构的随机响应和材料的疲劳性能参数,可以使用ANSYS的疲劳分析功能来计算结构的疲劳寿命。

通过以上步骤,可以利用ANSYS的随机振动分析功能实现随机疲劳分析。

这种方法能够更全面地考虑结构在实际工作环境下的疲劳寿命,为结构的设计和改进提供准确的参考。

需要注意的是,在进行随机疲劳分析时,需要对随机载荷进行合理的统计分析,获取载荷的概率密度函数。

如果没有足够的载荷数据,也可以使用统计模型进行估计。

此外,还需要对材料的疲劳性能参数进行准确的测定,以保证疲劳寿命计算的准确性。

总之,利用ANSYS的随机振动分析功能进行随机疲劳分析是一种有效的方法,可以更准确地评估结构在随机工作载荷下的疲劳性能,为结构的设计和改进提供有力的支持。

ANSYS 高清晰 精品资料:第15章 随机振动和随机疲劳分析实例

ANSYS 高清晰 精品资料:第15章 随机振动和随机疲劳分析实例

第15章 随机振动和随机疲劳分析实例 谱分析是一种将模态分析的结果与一个已知的谱联系起来计算模型的位移和应力的分析技术,主要用于确定结构对随机载荷或随时间变化载荷(如地震、风载、海洋波浪、喷气发动机推力、火箭发动机振动等等)的动力响应情况。

谱是谱值与频率的关系曲线,它反映了时间—历程载荷的强度和频率信息。

谱分析分为:响应谱分析、动力设计分析方法(Dynamic Design Analysis Method ,DDAM)和功率谱密度(Power Spectral Density —PSD,也称为随机振动分析)。

其中,一个响应谱代表单自由度系统对一个时间—历程载荷函数的响应,它是一个响应与频率的关系曲线,其中响应可以是位移、速度、加速度、力等。

下面以一个梁—板壳结构在地震位移激励作用下的随机振动分析为例,来将讲解用ANSYS6.1进行随机振动分析的具体过程,对于其它类型的谱分析的基本分析步骤可参阅本书第10章。

15.1 问题描述某板—梁结构如图15.1所示,计算在Y方向的地震位移激励谱作用下整个结构的响应情况。

板—梁结构的基本尺寸如图13.1所示,地震谱如表15.1所示,其它材料属性和几何特性数据如下:图15.1 梁-板壳结构模型A3钢的材料特性:杨氏模量EX =2.1 E11 N/m 2,泊松比PRXY =0.3,密度DENS =7.8 E 3 Kg/m 3。

板壳:厚度=2E-3 m梁几何特性:截面面积=1.6E-5 m 2, 惯性矩=21.333E-12 m 4 , 宽度=4E-3 m ,高度=4E-3 m表15.1梁—板结构所受的谱表位移激励谱频率(Hz) 0.5 1.0 2.4 3.8 17 18 20 32 位移(×10-3m) 0.01 0.016 0.03 0.02 0.005 0.01 0.015 0.0115.2 建立模型跟所有在ANSYS6.1中进行的结构分析一样,建立本实例的有限元模型同样需要完成如下工作:指定分析标题,定义材料性能,定义单元类型,定义单元实常数,建立几何模型并进行有限元网格划分等。

Ansys nCode疲劳分析

Ansys nCode疲劳分析
ANSYS nCode Designlife 高级疲劳寿命分析
© 2015 ANSYS, Inc.
1
目录
1 疲劳理论介绍
2 ANSYS nCode Designlife功能特色 3 S-N疲劳分析概述及案例介绍 4 随机振动疲劳分析概述
© 2015 ANSYS, Inc.
2
疲劳分析理论
只通过静力来确定 材料的机械性能时 ,机械性能没有充 分反映材料在交变 载荷作用下的特性 。使用过程中往往 会发生突如其来的 破坏。
nCode分析-五框图
© 2015 ANSYS, Inc. 19
ANSYS nCode 的功能特色
时间序列 恒幅载荷
时间步载荷
温度载荷 Hybird载荷 振动载荷 Duty Cycle
nCode分析-五框图
© 2015 ANSYS, Inc. 20
ANSYS nCode 的功能特色
ANSYS nCode DesignLife 具有全面 的疲劳分析能力
临界平面法计算平面应力和在几个径向平面进行重新求解
− 默认是每10度 (18 个平面)
− 雨流计数法在每个平面上计算疲劳损伤 − 准则平面是最大损伤的平面
© 2015 ANSYS, Inc.
Sm Sa 1 Su Se
Sa =应力幅 Sm = 平均应力 Su = 极限抗拉强度 Se = 疲劳极限
Sm 2 Se Sa 1 Su
Gerber平均应力修正理论假定一个抛物线
© 2015 ANSYS, Inc.
理对疲劳的影响
− 表面因子用来调整材料s-n曲线
Ksur是三个用户定义因子的乘积
− Ksur = KTreatment * Kuser * Kroughness

利用ANSYS随机振动分析功能实现随机疲劳分析

利用ANSYS随机振动分析功能实现随机疲劳分析

利用ANSYS随机振动分析功能实现随机疲劳分析
ANSYS随机振动分析功能可以获得结构随机振动响应过程的各种统计参数(如:
均值、均方根和平均频率等),根据各种随机疲劳寿命预测理论就可以成功地预测结
构的随机疲劳寿命。

本文介绍了ANSYS随机振动分析功能,以及利用该功能,按照Steinberg提出的基于高斯分布和Miner线性累计损伤定律的三区间法进行ANSYS随机疲劳计算的具体过程。

1.随机疲劳现象普遍存在
在工程应用中,汽车、飞行器、船舶以及其它各种机械或零部件,大多是在随机
载荷作用下工作,当它们承受的应力水平较高,工作达到一定时间后,经常会突然发
生随机疲劳破坏,往往造成灾难性的后果。

因此,预测结构或零部件的随机疲劳寿命
是非常有必要的。

2.ANSYS随机振动分析功能介绍
ANSYS随机振动分析功能十分强大,主要表现在以下方面:
1.具有位移、速度、加速度、力和压力等PSD类型;
2.能够考虑a阻尼、阻尼、恒定阻尼比和频率相关阻尼比;
3.能够定义基础和节点PSD激励;
4.能够考虑多个PSD激励之间的相关程度:共谱值、二次谱值、空间关系和波传
播关系等;
5.能够得到位移、应力、应变和力的三种结果数据: 1位移解,1速度解和
1加速度解;
3.利用ANSYS随机振动分析功能进行疲劳分析的一般原理
在工程界,疲劳计算广泛采用名义应力法,即以S-N曲线为依据进行寿命估算的
方法,可以直接得到总寿命。

下面围绕该方法举例说明ANSYS随机疲劳分析的一般原理。

ANSYSWORKBENCH疲劳分析指南第三章

ANSYSWORKBENCH疲劳分析指南第三章

ANSYSWORKBENCH疲劳分析指南第三章发表时间:2009-2-21 作者: 安世亚太来源: e-works关键字: CAE ansys Workbench疲劳分析第三章不稳定振幅的疲劳在前面一章中,考察了恒定振幅和比例载荷的情况,并涉及到最大和最小振幅在保持恒定的情况下的循环或重复载荷。

在本章将针对不定振幅、比例载荷情况,尽管载荷仍是成比例的,但应力幅和平均应力却是随时间变化的。

3.1 不规律载荷的历程和循环(History and Cycles)对于不规律载荷历程,需要进行特殊处理:计算不规律载荷历程的循环所使用的是“雨流”rainflow循环计算,“雨流”循环计算(Rainflowcycle counting)是用于把不规律应力历程转化为用于疲劳计算的循环的一种技术(如右面例子),先计算不同的“平均”应力和应力幅(“range”)的循环,然后使用这组“雨流”循环完成疲劳计算。

损伤累加是通过Palmgren-Miner 法则完成的,Palmgren-Miner法则的基本思想是:在一个给定的平均应力和应力幅下,每次循环用到有效寿命占总和的百分之几。

对于在一个给定应力幅下的循环次数Ni,随着循环次数达到失效次数Nfi时,寿命用尽,达到失效。

“雨流”循环计算和Palmgren-Miner损伤累加都用于不定振幅情况。

因此,任何任意载荷历程都可以切分成一个不同的平均值和范围值的循环阵列(“多个竖条”),右图是“雨流”阵列,指出了在每个平均值和范围值下所计算的循环次数,较高值表示这些循环的将出现在载荷历程中。

在一个疲劳分析完成以后,每个“竖条”(即“循环”)造成的损伤量将被绘出,对于“雨流”阵列中的每个“竖条”(bin),显示的是对应的所用掉的寿命量的百分比。

在这个例子中,即使大多数循环发生在低范围/平均值,但高范围(range)循环仍会造成主要的损伤。

依据Per Miner法则,如果损伤累加到1(100%),那么将发生失效。

基于ANSYS的专用喷洒车主管路随机疲劳分析

基于ANSYS的专用喷洒车主管路随机疲劳分析

减小 , 但超过一百万 次后最大应力不再 发生变化 , 曲线维持 且 在 2 × 0P 处 。该应力称为完全交变 应力 下的材料疲劳持 . 18a 3 久极 限。对 于铝合金试件 , 能找到这样 的材料疲劳持久极 就不 限。有几种方程可 以近似描述 s N曲线 , 简单、 — 最 最可接受 的
PD S 激励谱 r. 1 系统的传递函数I I 而△j一 — 正

=二] 旺
l 频率响应分析 l l 矩 } 惯 强 直至 苫
l型 模
堙 劳析 I I命 疲分 寿l
结果分析
本文采用疲劳寿命数值仿真分析方法 , 通过动力 学仿 真计
算, 得到结构在频域上 的动应力分布和传递 函数 , 根据这些 动 应 力分布 , 应用 A S S中 F tu 疲劳 分析模块 , NY ag e i 计算管路 的
是幂 函数形式 :
性能参数和损伤模 型计算疲劳 寿命 。 疲劳分析的简单流程如图
1 所示。 文根据线性 累积损伤理论 , 本 采用功率谱法 , 在获得结
构 动态 应力响应 P D 功率 谱密度 ) S( 和 N曲线后 , 计算结构
收稿 日期 :0 8 0 - 5 20 - 2 0 作者简介 : 申祖武 (9 & ) , 南固始人 , 16 一 , 河 男 博士 , 副研究 员, 主要研究领域为机械 动力学 、 动控 制、 振 防护工程等 ; 许志贵 (9 4 )男 , 17 一 , 湖北
t ~ ~ 一
式 中 , 是疲劳强度系数 , 由 S N曲线外推到第一个半 , 是 -
循 环 2 : ) ( 1 的应力 幅值 ; 6是疲 劳强度 指数 , 双对数 坐标 是 中 S N曲线的斜率 。 -
疲劳寿 u 。

ANSYS疲劳分析

ANSYS疲劳分析

ANSYS疲劳分析ANSYS是一种流行的工程仿真软件,用于进行各种工程问题的有限元分析。

在工程实践中,疲劳分析是一个非常重要的领域。

疲劳是指材料在重复载荷作用下逐渐破坏的过程。

疲劳分析的目的是评估结构在实际使用条件下的寿命和性能。

ANSYS可以用来进行疲劳分析,通过确定应力和应变的分布,评估结构在长期使用中可能出现的问题。

在进行疲劳分析之前,首先要进行有限元模型的建立。

这包括将结构模型导入到ANSYS中,确定边界条件和加载条件等。

在进行疲劳分析时,首先要确定疲劳载荷的类型和大小。

这可以通过实验测量或数值模拟来获取。

然后,将载荷应用在结构模型上,并进行动态分析。

ANSYS可以模拟不同的载荷情况,例如正弦载荷、随机载荷和脉冲载荷等。

通过分析结果,可以获得结构在不同位置的应力和应变分布。

在完成动态分析后,可以对结果进行验证和修正。

如果分析的结果与实际测量不符,可能需要对模型进行修正。

修正的方法包括调整材料的本构模型、改变模型的几何形状或重新定义载荷条件等。

完成验证后,可以进行疲劳分析。

在ANSYS中,可以使用不同的疲劳分析模块进行分析。

其中最常用的是疲劳寿命评估模块。

该模块可以根据疲劳参数和材料的S-N曲线,预测结构在给定载荷下的疲劳寿命。

这可以帮助工程师评估结构的安全性和可靠性,并采取适当的措施来延长结构的使用寿命。

疲劳分析还可以进行应力寿命曲线分析。

该分析方法可以通过建立不同应力水平和循环数的组合,预测结构的疲劳寿命。

这对于识别结构中的关键部位和进行寿命预测非常有帮助。

此外,还可以使用应变寿命方法进行疲劳分析。

该方法通过应变历程和损伤累积,评估结构在疲劳载荷下的性能。

在完成疲劳分析后,可以对结果进行后处理。

这包括评估结构的疲劳寿命、疲劳裕度和故障位置等。

通过分析结果,可以确定哪些部位可能会在疲劳过程中发生破坏,并采取适当的措施来加强这些部位。

总之,ANSYS是进行疲劳分析的强大工具。

它可以用于建立结构模型、应用载荷、进行动态分析和预测结构的疲劳寿命。

ansys-workbench疲劳分析流程

ansys-workbench疲劳分析流程

ansys-workbench疲劳分析流程ansys workbench困乏分析流程基于S-N曲线的困乏分析的最终目的是将变化无规律的多轴应力转化为简单的单轴应力循环,以便查询S-N曲线,得到相应的困乏寿命。

ansysworkbench 的困乏分析模块采用如下流程,其中r=Smin/Smax,Sa为应力幅度,Sm应力循环中的应力均值,注意后一个m不是大写:):(1)无规律多轴应力-->无规律单轴应力这个转换其实就是采用何种应力(或分量)。

只能有以下选择:Von-Mises等效应力;最大剪应力;最大主应力;或某一应力分量(Sx,Syz 等等)。

有时也采用带符号的Mises应力(大小不变等于Mises应力,符号取最大主应力的符号,好处是可以考虑拉或压的影响(反映在平衡应力或r上))。

同强度理论类似,Von-Mises等效应力和最大剪应力转换适用于延展性较好的材料,最大主应力转换用于脆性材料。

(2)无规律单轴应力-->简单单轴应力循环其本质是从无规律的高高低低的等效单轴应力--时间曲线中提取出一系列的简单应力循环(用Sa,Sm表征)以及对应的次数。

有很多种方法可以完成此计数和统计工作,其中又分为路径相关方法和路径无关方法。

用途最广的雨流法(rainflowcountingmethod)就是一种路径相关方法。

其算法和原理可见“Downing, S., Socie, D. (1982) Simplified rain flow counting algorithms. Int J Fatigue,4, 31–40“。

经过雨流法的处理后,无规律的应力--时间曲线转化为一系列的简单循环(Sa,Sm和ni,ni为该循环的次数,Sm如果不等于0,即r!=-1,需要考虑r的影响)。

然后将r!=-1的循环再转化到r=-1对应的应力循环(见下),这样就可以根据损伤累计理论(Miner准则)计算分析了:Sum(ni/Ni) Ni为该应力循环对应的寿命(考虑Sa,Sm)。

ANSYS随机振动分析教程

ANSYS随机振动分析教程

ANSYS随机振动分析教程随机振动分析(PSD: Power Spectral Density)是一种分析结构在随机动力加载下的响应特性的方法。

它通常应用于评估结构在实际工作环境中的可靠性和耐久性。

在ANSYS中进行随机振动分析,可以帮助我们理解结构在随机加载下的响应特性,评估结构的可靠性,并优化结构以提高其性能。

下面是一个基于ANSYS的随机振动分析的教程,通过该教程,你可以学习如何进行随机振动分析并分析结果。

步骤1:设置工程环境首先,打开ANSYS软件,并创建一个新的工程。

选择适当的单位制和求解器(如Mechanical APDL)。

步骤2:定义结构模型在这个教程中,我们将使用一个简单的悬臂梁作为结构模型。

创建一个梁模型,定义边界条件和加载条件。

确保模型准确代表了你想要分析的实际结构。

步骤3:定义随机负载在随机振动分析中,我们需要定义随机负载。

常见的随机负载包括自然地震、风荷载、机械振动等。

在这个教程中,我们以自然地震为例进行分析。

在ANSYS中,我们可以通过定义Power Spectral Density (PSD)函数来表示随机负载。

PSD函数描述了随机振动的能量分布,并用频率域表征。

使用命令“PSDZONE”创建一个PSD区域,然后通过命令“PSDFCN”定义一个PSD函数,并将PSD函数分配给PSD区域。

例如,你可以使用如下命令定义一个具有特定频率和幅值的PSD函数: /PSDZONE,1,FREQUENCY,1,200,AMP,0.1/PSDFCN,1,PSD,1步骤4:随机分析设置在进行随机振动分析之前,我们需要进行一些设置。

首先,我们需要定义分析的频率范围和步长。

可以使用命令“FREQSEP”来定义频率范围和步长。

例如,你可以使用如下命令定义频率范围为1Hz到200Hz,步长为1Hz:/FREQSEP,1HZ,1接下来,我们需要定义求解器参数。

使用命令“MODAL”定义模态分析参数:/COMBINATION,PSD/PSD,1,UNDEF然后,定义DAREA区域,并通过命令“SDPOINT”为每个频率分配一个积木节点。

ANSYSworkbench教程—疲劳分析

ANSYSworkbench教程—疲劳分析

第一章简介1.1 疲劳概述结构失效的一个常见原因是疲劳,其造成破坏与重复加载有关。

疲劳通常分为两类:高周疲劳是当载荷的循环(重复)次数高(如1e4 -1e9)的情况下产生的。

因此,应力通常比材料的极限强度低,应力疲劳(Stress-based)用于高周疲劳;低周疲劳是在循环次数相对较低时发生的。

塑性变形常常伴随低周疲劳,其阐明了短疲劳寿命。

一般认为应变疲劳(strain-based)应该用于低周疲劳计算。

在设计仿真中,疲劳模块拓展程序(Fatigue Module add-on)采用的是基于应力疲劳(stress-based)理论,它适用于高周疲劳。

接下来,我们将对基于应力疲劳理论的处理方法进行讨论。

1.2 恒定振幅载荷在前面曾提到,疲劳是由于重复加载引起:当最大和最小的应力水平恒定时,称为恒定振幅载荷,我们将针对这种最简单的形式,首先进行讨论。

否则,则称为变化振幅或非恒定振幅载荷。

1.3 成比例载荷载荷可以是比例载荷,也可以非比例载荷:比例载荷,是指主应力的比例是恒定的,并且主应力的削减不随时间变化,这实质意味着由于载荷的增加或反作用的造成的响应很容易得到计算。

相反,非比例载荷没有隐含各应力之间相互的关系,典型情况包括:σ1/σ2=constant在两个不同载荷工况间的交替变化;交变载荷叠加在静载荷上;非线性边界条件。

1.4 应力定义考虑在最大最小应力值σmin和σmax作用下的比例载荷、恒定振幅的情况:应力范围Δσ定义为(σmax-σmin)平均应力σm定义为(σmax+σmin)/2应力幅或交变应力σa是Δσ/2应力比R是σmin/σmax当施加的是大小相等且方向相反的载荷时,发生的是对称循环载荷。

这就是σm=0,R=-1的情况。

当施加载荷后又撤除该载荷,将发生脉动循环载荷。

这就是σm=σmax/2,R=0的情况。

1.5 应力-寿命曲线载荷与疲劳失效的关系,采用的是应力-寿命曲线或S-N曲线来表示:(1)若某一部件在承受循环载荷, 经过一定的循环次数后,该部件裂纹或破坏将会发展,而且有可能导致失效;(2)如果同个部件作用在更高的载荷下,导致失效的载荷循环次数将减少;(3)应力-寿命曲线或S-N曲线,展示出应力幅与失效循环次数的关系。

振动疲劳—ansys

振动疲劳—ansys

振动疲劳—ansys随机振动疲劳分析随机振动疲劳分析流程图随机振动疲劳分析将第一步频率响应分析得到的结果文件作为输入,并在疲劳软件中输入振动过程中的PSD曲线,经计算得到零件的振动疲劳寿命。

故随机振动疲劳分析可分为如下步骤:1.频率响应分析结果输入2.功率谱密度PSD输入3.材料疲劳特性设置4.各工况与PSD关联设置5.振动疲劳求解器参数设置6.输出设置7.分析结果处理频率响应分析结果输入功率谱密度PSD输入振动疲劳求解器Ncode云图显示输出设置Ncode随机振动疲劳分析流程图1.频率响应分析结果输入频率响应分析应与PSD 的单位相对应,比如PSD 单位为g^2/Hz ,则进行频率响应分析时可输入1g 的加速度激励来分析。

(如采取单位制ton-mm-s-N ,此时1g 的加速度激励为9800mm/s^2,应在分析中输入9800大小的加速度激励)1.1单位问题1.2频率响应分析结果输出设置为了避免输出结果过大,可以在输出中设置需要进行疲劳分析的部件,以set 形式输出,同时可设置输出频次Frequency=n ,只输出频响分析应力结果即可。

*OUTPUT, FIELD, Frequency=5 *ELEMENT OUTPUT, ELSET = ele_setS, 以Abaqus 进行频响分析为例,输出设置如下:每5步输出一次只输出单元集合名为ele_set 的应力结果2.功率谱密度PSD输入PSD可以用以下2种方式输入:1.通过MultiColumnInput读入定义好的CSV文件输入2.通过VibrationGenerator生成PSD 曲线CSV文件格式如下:(可在帮助文档中找一个PSD的CSV文件作为模板)。

振动疲劳—ansys

振动疲劳—ansys

随机振动疲劳分析随机振动疲劳分析流程图随机振动疲劳分析将第一步频率响应分析得到的结果文件作为输入,并在疲劳软件中输入振动过程中的PSD曲线,经计算得到零件的振动疲劳寿命。

故随机振动疲劳分析可分为如下步骤:1.频率响应分析结果输入2.功率谱密度PSD输入3.材料疲劳特性设置4.各工况与PSD关联设置5.振动疲劳求解器参数设置6.输出设置7.分析结果处理频率响应分析结果输入功率谱密度PSD输入振动疲劳求解器Ncode云图显示输出设置Ncode随机振动疲劳分析流程图1.频率响应分析结果输入频率响应分析应与PSD 的单位相对应,比如PSD 单位为g^2/Hz ,则进行频率响应分析时可输入1g 的加速度激励来分析。

(如采取单位制ton-mm-s-N ,此时1g 的加速度激励为9800mm/s^2,应在分析中输入9800大小的加速度激励)1.1单位问题1.2频率响应分析结果输出设置为了避免输出结果过大,可以在输出中设置需要进行疲劳分析的部件,以set 形式输出,同时可设置输出频次Frequency=n ,只输出频响分析应力结果即可。

*OUTPUT, FIELD, Frequency=5*ELEMENT OUTPUT, ELSET = ele_setS, 以Abaqus 进行频响分析为例,输出设置如下:每5步输出一次只输出单元集合名为ele_set 的应力结果2.功率谱密度PSD输入PSD可以用以下2种方式输入:1.通过MultiColumnInput读入定义好的CSV文件输入2.通过VibrationGenerator生成PSD 曲线CSV文件格式如下:(可在帮助文档中找一个PSD的CSV文件作为模板)。

基于Ansys二次开发的随机振动疲劳寿命可视化程序开发

基于Ansys二次开发的随机振动疲劳寿命可视化程序开发

基于Ansys二次开发的随机振动疲劳寿命可视化程序开发
丁杰;唐玉兔
【期刊名称】《机械强度》
【年(卷),期】2014(36)5
【摘要】为解决现有随机振动分析中存在的结构各节点的振动频率计算、三个方
向振动结果的叠加、全域可视化显示等问题,进行了基于Ansys软件的二次程序开发。

通过程序在某变流装置中的应用,表明了二次开发的程序可以更加快速、全面、有效地进行虚拟样机的疲劳寿命分析,为产品结构的优化设计提供理论指导。

【总页数】6页(P813-818)
【关键词】Ansys;随机振动;加速度谱密度;疲劳损伤;疲劳寿命;可视化
【作者】丁杰;唐玉兔
【作者单位】南车株洲电力机车研究所有限公司南车电气技术与材料工程研究院【正文语种】中文
【中图分类】TB123;U264.37
【相关文献】
1.基于ANSYS的机床滑合随机振动疲劳寿命分析 [J], 孟行;孟广耀;黄居鑫;石英训;生开明
2.基于Ansys的车桥结构随机振动疲劳寿命分析 [J], 王鹏利
3.基于ANSYS二次开发的大跨桥梁疲劳寿命与失效概率的可视化技术 [J], 侯士通;王莹
4.基于ANSYS的连接件随机振动疲劳寿命分析及优化设计 [J], 白永明; 邱恩举; 王宏建
5.基于ANSYS的连接件随机振动疲劳寿命分析及优化设计 [J], 白永明; 邱恩举; 王宏建
因版权原因,仅展示原文概要,查看原文内容请购买。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

利用ANSYS随机振动分析功能实现随机疲劳分析
ANSYS随机振动分析功能可以获得结构随机振动响
应过程的各种统计参数(如:均值、均方根和平均频率等),根据各种随机疲劳寿命预测理论就可以成功地预测结构
的随机疲劳寿命。

本文介绍了ANSYS随机振动分析功能,以及利用该功能,按照Steinberg提出的基于高斯分布和Miner线性累计损伤定律的三区间法进行ANSYS随机疲劳计算的具体过程。

1.随机疲劳现象普遍存在
在工程应用中,汽车、飞行器、船舶以及其它各种机械或零部件,大多是在随机载荷作用下工作,当它们承受的应力水平较高,工作达到一定时间后,经常会突然发生随机疲劳破坏,往往造成灾难性的后果。

因此,预测结构或零部件的随机疲劳寿命是非常有必要的。

2.ANSYS随机振动分析功能介绍
ANSYS随机振动分析功能十分强大,主要表现在以下方面:
1.具有位移、速度、加速度、力和压力等PSD类型;
2.能够考虑a阻尼、 阻尼、恒定阻尼比和频率相关阻
尼比;
3.能够定义基础和节点PSD激励;
4.能够考虑多个PSD激励之间的相关程度:共谱值、二
次谱值、空间关系和波传播关系等;
5.能够得到位移、应力、应变和力的三种结果数据: 1σ
位移解,1σ速度解和1σ加速度解;
3.利用ANSYS随机振动分析功能进行疲劳分析的一般原

在工程界,疲劳计算广泛采用名义应力法,即以S-N 曲线为依据进行寿命估算的方法,可以直接得到总寿命。

下面围绕该方法举例说明ANSYS随机疲劳分析的一般原理。

当应力历程是随机过程时,疲劳计算相对比较复杂。

但已经有许多种分析方法,这里仅介绍一种比较简单的方法,即Steinberg提出的基于高斯分布和Miner线性累计损伤定律的三区间法(应力区间如图1所示):
应力区间
发生的时

-1σ ~+1σ68.3%的时间
-2σ ~+2σ27.1%的时间
-3σ ~+3σ4.33%的时间
99.73%
大于3σ的应力仅仅发生在0.27%的时间内,假定其不造成任何损伤。

在利用Miner定律进行疲劳计算时,将应力处理成上述3个水平,总体损伤的计算公式就可以写成:
其中:
:等于或低于1σ水平的实际循环数目(0.6831 );
:等于或低于2σ水平的实际循环数目(0.271 );
:等于或低于3σ水平的实际循环数目(0.0433 );
, , :根据疲劳曲线查得的1σ、2σ和3σ应力水平分别对应许可循环的次数。

综上所述,针对Steinberg提出的基于高斯分布和Miner 线性累计损伤定律的三区间法的ANSYS随机疲劳分析的一般过程是:
(1) 计算感兴趣的应力分量的统计平均频率(应力速度/应力);
(2) 基于期望(工作)寿命和统计平均频率,计算1 ,2 和3 水平下的循环次数、和;
(3) 基于S-N曲线查表得到、和;
(4) 计算疲劳寿命使用系数。

显然,根据其他随机疲劳分析方法和ANSYS随机振动分析结果,我们还可以进行许多类似的疲劳分析计算。

相关文档
最新文档