ANSYS 高清晰 精品资料：第15章 随机振动和随机疲劳分析实例
疲劳分析的数值计算方法及ANSYS疲劳分析实例
及实例第一节引言零件或构件由于交变载荷的反复作用，在它所承受的交变应力尚未达到静强度设计的许用应力情况下就会在零件或构件的局部位置产生疲劳裂纹并扩展、最后突然断裂。
这种现象称为疲劳破坏。
疲劳裂纹的形成和扩展具有很大的隐蔽性而在疲劳断裂时又具有瞬发性，因此疲劳破坏往往会造成极大的经济损失和灾难性后果。
金属的疲劳破坏形式和机理不同与静载破坏，所以零件疲劳强度的设计计算不能为经典的静强度设计计算所替代，属于动强度设计。
随着机车车辆向高速、大功率和轻量化方向的迅速发展，其疲劳强度及其可靠性的要求也越来越高。
近几年随着我国铁路的不断提速，机车、车辆和道轨等铁路设施的疲劳断裂事故不断发生，越来越引起人们的重视。
疲劳强度设计及其研究正在成为我国高速机车车辆设计制造中的一项不可缺少的和重要的工作。
金属疲劳的研究已有近150年的历史，有相当多的学者和工程技术人员进行了大量的研究，得到了许多关于金属疲劳损伤和断裂的理论及有关经验技术。
但是由于疲劳破坏的影响因素多而复杂并且这些因素互相影响又与构件的实际情况密切相关，使得其应用性成果尚远远不能满足工程设计和生产应用的需要。
据统计，至今有约90%的机械零部件的断裂破坏仍然是由直接于疲劳或者间接疲劳而引起的。
因此，在21世纪的今天，尤其是在高速和大功率化的新产品的开发制造中，其疲劳强度或疲劳寿命的设计十分重要，并且往往需要同时进行相应的试验研究和试验验证。
疲劳断裂是因为在零件或构件表层上的高应力或强度比较低弱的部位区域产生疲劳裂纹，并进一步扩展而造成的。
这些危险部位小到几个毫米甚至几十个微米的范围，零件或构件的几何缺口根部、表面缺陷、切削刀痕、碰磕伤痕及材料的内部缺陷等往往是这种危险部位。
因此，提高构件疲劳强度的基本途径主要有两种。
一种是机械设计的方法，主要有优化或改善缺口形状，改进加工工艺工程和质量等手段将危险点的峰值应力降下来；另一种是材料冶金的方法，即用热处理手段将危险点局部区域的疲劳强度提高，或者是提高冶金质量来减少金属基体中的非金属夹杂等材料缺陷等局部薄弱区域。
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疲劳分析
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 Mechanical 动力学_L05_随机振动
29 © 2016 ANSYS, Inc.
September 23, 2016
…PSD 曲线拟合
• 为了解决拟合优度问题： • 点击弹出的载荷数据选项，并选择改良的拟合
• 如果插值点在拟合优度近似中可用则显示
30 © 2016 ANSYS, Inc.
September 23, 2016
E. 求解设置
• 这些随机振动同时包含了所有频率 • 这些频率的振动幅值随时间随机变化
• 我们需要一些方法来描述和衡量这些激励
5
© 2016 ANSYS, Inc.
September 23, 2016
…定义和目的
• 如果幅值持续变化，如何评估一个随机激励
• 重点关注：在给定频率范围内，激励的幅值确实持续变化，但对很多 过程而言，其均值趋向于保持相对稳定
ANSYS Mechanical 动力学技术专题培训
随机振动分析
1
© 2016 ANSYS, Inc.
September 23, 2016
简介
在这一章中，将介绍随机振动分析 本章内容： A. 什么是随机振动 B. 功率谱密度 PSD C. 理论概述 D. PSD 曲线拟合 E. 求解设置 F. 作业4
• PSD 曲线定义为分段线性频率表，并以如下方式在对数空间下绘制
• 曲线拟合多项式用于曲线的积分解析表达式 • 为了得到较好的拟合，在连续点间的 PSD 值应在一个数量级内变化
28 © 2016 ANSYS, Inc.
September 23, 2016
…PSD 曲线拟合
• 一旦输入了载荷，图中提供一种 颜色标识以指示每一区段 • 绿色：认为值可信和准确 • 黄色：警告提示，认为产生 的结果不可信不准确 • 红色：产生的结果不可信。 建议你在求解前改变 PSD 输 入载荷
随机振动分析实例
ANSYS 动力分析(18) - 随机振动分析- 实例(1)2010-09-26 07:41:23| 分类：ANSYS 动力分析| 标签：随机振动实例模型飞机机翼psd|举报|字号订阅PSD 实例：模型飞机机翼的随机振动说明：确定由于施加在机翼根部的Y 向加速度PSD，在模型飞机机翼中造成的位移和应力。
假设机翼在Z=0 处固支。
操作指南1. 清除数据库并读入文件wing. inp 以创建几何模型和网格。
2. 定义材料属性：弹性模量= 38000 psi 泊松比= 0.3密度= 1.033E-3/12 lbf-sec2/in4 = 8.6083E-53. 施加边界条件。
提示：选择在areas 上施加位移约束，拾取Z=0 处所有的Areas，约束所有自由度。
4. 定义新分析为Model，使用Block Lanczos 方法，抽取和扩展前15 个自然模态。
然后求解Current LS。
5. 查看模态形状，如图为前 4 阶振型。
6. 使用所显示的 PSD 谱，执行 PSD Spectrum 分析。
首先定义分析类型为 Spectrum分析类型为 PSD，使用全部模态，计算单元应力：注意激活“Calculate elem stresses”选项。
7. 在基础上施加指定的 PSD 谱 (注意：确保 PSD 的单位是 G2/Hz)。
施加 Y 向激励 (方法是：在基础节点上施加单位 Y 向位移)。
设置常阻尼比 0.02：设置有关参数–重力加速度值注意：响应谱类型选择 Accel (g**2/Hz)，否则后面的 PSD 谱应该输入实际加速度值：定义 PSD 谱表格：绘制 PSD 谱表格曲线以检查输入值：激活和设置模态组合参数：选择所有模态。
设置计算内容：计算各模态的参与因子:在输出窗口中可以看到参与因子的计算结果：求解：10. 在一般后处理中，查看相对位移和应力 (载荷步 3)。
首先查看 Results Summary：载荷步 3 为 set 17查看载荷步 3 的相对位移和应力：在 Read Results 中读取载荷步 3 的结果：使用 By set Number，输入 Data Set Number 为 17。
利用ANSYS随机振动分析功能实现随机疲劳分析
利用ANSYS随机振动分析功能实现随机疲劳分析ANSYS是一款常用的工程仿真软件，具有强大的分析功能。
在进行随机疲劳分析时，可以利用ANSYS的随机振动分析功能来模拟随机加载下的疲劳损伤。
随机疲劳分析是一种考虑工作载荷随机性对结构疲劳寿命影响的方法。
通过采用随机振动分析，可以考虑到工作载荷的随机特性，进一步分析结构的疲劳损伤。
ANSYS中的随机振动分析功能可以通过以下步骤来进行：1.几何建模：首先，需要进行结构的几何建模。
使用ANSYS的几何建模工具可以创建出要进行疲劳分析的结构。
2.材料属性定义：在进行材料属性的定义时，需要确定材料的弹性模量、泊松比、密度和疲劳参数等。
可以根据材料的材料数据手册来获取这些参数。
3.边界条件设置：在进行随机振动分析时，需要设置结构的边界条件。
这些边界条件可以是结构受到的随机外载荷或者是结构与其他部件的接触情况。
4.加载设置：在进行随机振动分析时，需要设置结构受到的随机载荷。
这些载荷可以是来自于实际工况的随机载荷，也可以通过振动台试验数据等手段获取。
5.随机振动分析：利用ANSYS的随机振动分析功能，可以进行频域分析或时域分析。
频域分析可以用于计算结构的响应功率谱密度，时域分析可以用于计算结构的随机响应。
6.疲劳寿命计算：在获得结构的随机响应后，可以进行疲劳寿命计算。
根据结构的随机响应和材料的疲劳性能参数，可以使用ANSYS的疲劳分析功能来计算结构的疲劳寿命。
通过以上步骤，可以利用ANSYS的随机振动分析功能实现随机疲劳分析。
这种方法能够更全面地考虑结构在实际工作环境下的疲劳寿命，为结构的设计和改进提供准确的参考。
需要注意的是，在进行随机疲劳分析时，需要对随机载荷进行合理的统计分析，获取载荷的概率密度函数。
如果没有足够的载荷数据，也可以使用统计模型进行估计。
此外，还需要对材料的疲劳性能参数进行准确的测定，以保证疲劳寿命计算的准确性。
总之，利用ANSYS的随机振动分析功能进行随机疲劳分析是一种有效的方法，可以更准确地评估结构在随机工作载荷下的疲劳性能，为结构的设计和改进提供有力的支持。
Ansys培训-随机振动分析
15. In the Details of the PSD Load,
14
change “Direction” to “Y Axis” for
this particular XYZ orientation.
16. For >Load Data chose >New PSD Load
Acceleration
– The data points can be entered for each Freq & Amplitude, or a function can be entered.
Acceleration
A2
A3
A1 A4
F1 F2
F3
F4
Frequency
Workshop – 假定
• The Girder has fixed constraints along all lower edges.
2. Click OK, thus accepting the default number of modes
3. Choose the U.S. inch pound unit system.
– “Units > U.S. Customary (in, lbm, lbf, …)” 1
2
3
Workshop – 前处理-壳体厚度
PSD分析. • Steps: 进行模态和随
机振动分析，并显示 结果.
随机振动分析
随机振动分析流程 打开, Tower.dsdb.
• Browse to file if not in list
• 打开分析向导…
随机振动分析
随机振动分析流程
随机振动分析
ANSYS-随机振动分析功能实现随机疲劳分析
利用 ANSYS随机振动分析功能实现随机疲劳分析ANSYS 随机振动分析功能可以获得结构随机振动响应过程的各种统计参数（如：均值、均方根和平均频率等），根据各种随机疲劳寿命预测理论就可以成功地预测结构的随机疲劳寿命。
本文介绍了ANSYS随机振动分析功能，以及利用该功能，按照Steinberg 提出的基于高斯分布和Miner 线性累计损伤定律的三区间法进行 ANSYS随机疲劳计算的具体过程。
１．随机疲劳现象普遍存在在工程应用中，汽车、飞行器、船舶以及其它各种机械或零部件，大多是在随机载荷作用下工作，当它们承受的应力水平较高，工作达到一定时间后，经常会突然发生随机疲劳破坏，往往造成灾难性的后果。
因此，预测结构或零部件的随机疲劳寿命是非常有必要的。
2．ANSYS随机振动分析功能介绍ANSYS随机振动分析功能十分强大，主要表现在以下方面：1.具有位移、速度、加速度、力和压力等 PSD类型；2.能够考虑 a 阻尼、阻尼、恒定阻尼比和频率相关阻尼比；3.能够定义基础和节点 PSD激励；4.能够考虑多个 PSD激励之间的相关程度：共谱值、二次谱值、空间关系和波传播关系等；5.能够得到位移、应力、应变和力的三种结果数据： 1 位移解， 1 速度解和 1 加速度解；3．利用 ANSYS随机振动分析功能进行疲劳分析的一般原理在工程界，疲劳计算广泛采用名义应力法，即以S-N 曲线为依据进行寿命估算的方法，可以直接得到总寿命。
下面围绕该方法举例说明 ANSYS随机疲劳分析的一般原理。
当应力历程是随机过程时，疲劳计算相对比较复杂。
但已经有许多种分析方法，这里仅介绍一种比较简单的方法，即Steinberg 提出的基于高斯分布和 Miner 线性累计损伤定律的三区间法（应力区间如图 1 所示）：应力区间发生的时间68.3% 的时-1~+1间27.1% 的时-2~+2间4.33% 的时-3~+3间99.73%大于 3的应力仅仅发生在0.27%的时间内，假定其不造成任何损伤。
随机振动分析实例
ANSYS 动力分析(18) - 随机振动分析- 实例(1)2010-09-26 07:41:23| 分类：ANSYS 动力分析| 标签：随机振动实例模型飞机机翼psd|举报|字号订阅PSD 实例：模型飞机机翼的随机振动说明：确定由于施加在机翼根部的Y 向加速度PSD，在模型飞机机翼中造成的位移和应力。
假设机翼在Z=0 处固支。
操作指南1. 清除数据库并读入文件wing. inp 以创建几何模型和网格。
2. 定义材料属性：弹性模量= 38000 psi 泊松比= 0.3密度= 1.033E-3/12 lbf-sec2/in4 = 8.6083E-53. 施加边界条件。
提示：选择在areas 上施加位移约束，拾取Z=0 处所有的Areas，约束所有自由度。
4. 定义新分析为Model，使用Block Lanczos 方法，抽取和扩展前15 个自然模态。
然后求解Current LS。
5. 查看模态形状，如图为前 4 阶振型。
6. 使用所显示的 PSD 谱，执行 PSD Spectrum 分析。
首先定义分析类型为 Spectrum分析类型为 PSD，使用全部模态，计算单元应力：注意激活“Calculate elem stresses”选项。
7. 在基础上施加指定的 PSD 谱 (注意：确保 PSD 的单位是 G2/Hz)。
施加 Y 向激励 (方法是：在基础节点上施加单位 Y 向位移)。
设置常阻尼比 0.02：设置有关参数–重力加速度值注意：响应谱类型选择 Accel (g**2/Hz)，否则后面的 PSD 谱应该输入实际加速度值：定义 PSD 谱表格：绘制 PSD 谱表格曲线以检查输入值：激活和设置模态组合参数：选择所有模态。
设置计算内容：计算各模态的参与因子:在输出窗口中可以看到参与因子的计算结果：求解：10. 在一般后处理中，查看相对位移和应力 (载荷步 3)。
首先查看 Results Summary：载荷步 3 为 set 17查看载荷步 3 的相对位移和应力：在 Read Results 中读取载荷步 3 的结果：使用 By set Number，输入 Data Set Number 为 17。
Ansys培训_随机振动分析汇总
DYNAMICS 11.0
A2
A3
Acceleration
A1 A4
F1
F2
F3
F4
Frequency
Workshop – 假定
• 和确定性谱分析不同，随机振动不能用瞬态动 力学分析代替. • 应用基于概率的功率谱密度分析，分析载荷作 用过程中的统计规律
DYNAMICS 11.0
Image from “Random Vibrations Theory and Practice” by Wirsching, Paez and Ortiz.
定义和目的
输入:
– 结构的自然频率和阵型 – 功率谱密度曲线
Training Manual
DYNAMICS 11.0
输出:
– 1s位移和应力 （用于疲劳分析）.
• 载荷:
– 单点激励
Training Manual
• 得到结果:
– 相对或绝对的1s 输出 – 整体结构的结果，可以进行云图显示. – 1s位移, 速度或加速度
Training Manual
DYNAMICS 11.0
随机振动分析
随机振动分析流程
• • • • 插入一个PSD Base Excitation. 在弹出的PSD Base Excitation详情串口，选择新的PSD载荷. 选择带G的加速度PSD，单位G^2/Hz. 设置PSD曲线
Training Manual
打开, Tower.dsdb.
Training Manual
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中进行的结构分析一样，建立本实例的有限元模型同样需要完成如下工作：指定分析标题，定义材料性能，定义单元类型，定义单元实常数，建立几何模型并进行有限元网格划分等。
ansysworkbench疲劳分析实例.
AHSYS作业A12疲劳作业M2-目标Sha^cFMer AN、、、• H ■:d■気ro tanavM・律》枫fi人・・.aewwwRt *4) 口•■fMTEMir »|・•£ FulN *屮・少4W0WW-sASNV mmuMJiAua t胡尺・乙2彳卜曲…sASNViMlOUViin.imWV ・毘m却…^pwu**vssASNV I吹OIU¥Q5T2H 冷.・■ »sASNV■»*>« I«in»ws. Mfbrn W■哥MHV■・J» D «d tiw^im-iuun i•(• W ^VUI> WMM||^fK»diJW ・砂-"仇#Shape…沁Ad过程9. «命*朋广■rvmawwn,10. 止心疫”臭金K*一1化怪人一牛■址!l雯鼻mRi . tpw«1 -• J.DMMBOt CvMrieal Support *11!12. rtW RadW Eed .・U・・rwlrangemiBlz Fr<«'ANSYS••eSh^fx hnfef 作求解杖IF ANSY、w・、. Rxh4iT 皂1Z RMB * tn*cn > Stress » Equivalent (wn M MIB RMB > biMd > OUcHiMtion > Totol.IL•・Sha^FMer…•作业A12-过粹入r»«m.Sha^FMer作业M2-求解设需• MAWIguoioei (•.倉出•oluMon 苣ANSYS・.2AIMSYSZMCTANSY<A i2—求样设栄•从 Dewhof gRgr ■“：21.尼乂M ■良也 山3*54『巧尹》318)便""o ■停■补尤一*样♦««Sha^FMer作业M2-求解设胃•・ Fell0wTo (iaMI ； 0■八・3吋 Fictor (*±se>:RMS ・ kiMcl > Fatigue • 3My Ftclof• azamof “MtryF■:tor ■口： »-仗・缸论呼佢皿9“ UK .卯 * 环次■34 4 IV OR Mhoe W7)U«R1 ■为M 仇 ft I itW. Shape FMcr作WA12-求駢设国ANSYS••r>・ iQWZYgeWf 集); 27.}tA*Fm*gue Sensitivity ^«* Sh^pcFMr AN 、Y ，作业A12 -求解设SF• AFMI9M TocMMl<«): » «AmtSe^ttpo-mtB • kwel • FMgue ■ Baiwiv NGtoon・ MDeiwtoofl FmweSmsitivny Su ；28WR-亍钛两賢4<»嗖冲 S(—訪 tttSW ffiizg 九力祥殳化吩 伽加的文号内力.Sha^cFMef AIMS YY作业也-缺• Atf)HatnVft£力輪旻w*hf«iww・湍—EY' • ■力卄■山旬fi«itueAi«erj・H次啟micf S^pcFMer AN、YS 作业M2—箪杲• Art)M«inJtJ t • f 4kAtt« 9$t«« WWW 个A.K«Aa«V Rm Y V2«xrMm«O0tttiif FM QM SeoMW1• ftLiixa A«« «ita«n «；*•的・・■|1>M»N ・tHk・竹结限•A，»»0 MJiaWKtr^» ■•MityinactMo廿・4 •曲蚯fllQWftl ・力轶*精近巾■的(U) ■凶为H IWV-fl • f・an«9, -to btnMy, f*巨力#WU ・”LWfc>!妙为何■杯 dSh^fxFinkf 判如2-处ANSY、Sht^FMer 作业M2-箪杲ANSY1—Shape Fln^f徉WA12-求解设匮ANSYS SharpeFMer• nxn •个孙gtiwti■■细住dt«ft*(nSttf«MV.hH ・■亍朴“栄屮L2Q0 号■lr“SOOM・|・ «WH«：30 Hl0hli0hrt I9W solution U BAC DAN、Y\ANSY 、ge FMtfANS Y 、作心22-戒倩设實 •出 MMe 0t TgT Su （•<） 35. /上OoodiMf 堆U 沖・锂力・・.M.Von BRia^ltW （铁用*pn ・lie G>0dEW **®K 的HR 的$列虫力昭氏小H ）. 37 «tbU alwMa “ KRalMlow «««»> Oom»Qe mMnccs 是 32K 3?|. Sh^ptFMer作业A2-求解设拧•・ FMl0MtiTocU0MI«； 3a ・八•金•■”'•：f. P h«vi p F«0flvt f MHShape Mcrf 打M12 ■求解设国AN 、Y 、•均 Filigue Tool （M ）: ・18n fatigue 8en84iYity ；AMB > Insert > 层A H QIM > ・・Ugs •从 D<H*I» o4 *FabQue SeositivHy* tt ：l ： 42.程文■變化・咬为50M 〔交仇 ««为NK>“ 0?耳虫力为 WOOH ）.W.・ AFactor':RNB • IwY • FsnQiit • MHy FaGKr• U Deum of SafHY Factor f IJ ：SMpe^Antfer作业A2-求解设拧•创 Fallgut Tool:• mm • Fsoove ・•・・・iy M GXMMI44.・八■■如・ fUlMlow Matfix* :/ tan » fottQiM / AeMrw MitruANSYS—ANSYSAN 、、你如2-桝 • ■帕Ml 山加 OOO&I 掰的3i*FeW ・ShapeMtf作血M2-求解设海ANSY 、ySh^ptFMer 作业入、2-帘果• A) Fnigue Tort # Ml« • (Mt) --Ownafle Matrix-：RMB> toMd>FMKM > Oerao* Matrix•flre«e・ AA.tfOttl*#** -LMe^・点帕&QtW。
利用ANSYS随机振动分析功能实现随机疲劳分析
利用ANSYS随机振动分析功能实现随机疲劳分析
ANSYS随机振动分析功能可以获得结构随机振动响应过程的各种统计参数（如：
均值、均方根和平均频率等），根据各种随机疲劳寿命预测理论就可以成功地预测结
构的随机疲劳寿命。
本文介绍了ANSYS随机振动分析功能，以及利用该功能，按照Steinberg提出的基于高斯分布和Miner线性累计损伤定律的三区间法进行ANSYS随机疲劳计算的具体过程。
１．随机疲劳现象普遍存在
在工程应用中，汽车、飞行器、船舶以及其它各种机械或零部件，大多是在随机
载荷作用下工作，当它们承受的应力水平较高，工作达到一定时间后，经常会突然发
生随机疲劳破坏，往往造成灾难性的后果。
因此，预测结构或零部件的随机疲劳寿命
是非常有必要的。
2．ANSYS随机振动分析功能介绍
ANSYS随机振动分析功能十分强大，主要表现在以下方面：
1.具有位移、速度、加速度、力和压力等PSD类型；
2.能够考虑a阻尼、阻尼、恒定阻尼比和频率相关阻尼比；
3.能够定义基础和节点PSD激励；
4.能够考虑多个PSD激励之间的相关程度：共谱值、二次谱值、空间关系和波传
播关系等；
5.能够得到位移、应力、应变和力的三种结果数据： 1位移解，1速度解和
1加速度解；
3．利用ANSYS随机振动分析功能进行疲劳分析的一般原理
在工程界，疲劳计算广泛采用名义应力法，即以S-N曲线为依据进行寿命估算的
方法，可以直接得到总寿命。
下面围绕该方法举例说明ANSYS随机疲劳分析的一般原理。
随机振动分析实例
ANSYS 动力分析(18) - 随机振动分析- 实例(1)2010-09-26 07:41:23| 分类：ANSYS 动力分析| 标签：随机振动实例模型飞机机翼psd|举报|字号订阅PSD 实例：模型飞机机翼的随机振动说明：确定由于施加在机翼根部的Y 向加速度PSD，在模型飞机机翼中造成的位移和应力。
假设机翼在Z=0 处固支。
操作指南1. 清除数据库并读入文件wing. inp 以创建几何模型和网格。
2. 定义材料属性：弹性模量= 38000 psi 泊松比= 0.3密度= 1.033E-3/12 lbf-sec2/in4 = 8.6083E-53. 施加边界条件。
提示：选择在areas 上施加位移约束，拾取Z=0 处所有的Areas，约束所有自由度。
4. 定义新分析为Model，使用Block Lanczos 方法，抽取和扩展前15 个自然模态。
然后求解Current LS。
5. 查看模态形状，如图为前 4 阶振型。
6. 使用所显示的 PSD 谱，执行 PSD Spectrum 分析。
首先定义分析类型为 Spectrum分析类型为 PSD，使用全部模态，计算单元应力：注意激活“Calculate elem stresses”选项。
7. 在基础上施加指定的 PSD 谱 (注意：确保 PSD 的单位是 G2/Hz)。
施加 Y 向激励 (方法是：在基础节点上施加单位 Y 向位移)。
设置常阻尼比 0.02：设置有关参数–重力加速度值注意：响应谱类型选择 Accel (g**2/Hz)，否则后面的 PSD 谱应该输入实际加速度值：定义 PSD 谱表格：绘制 PSD 谱表格曲线以检查输入值：激活和设置模态组合参数：选择所有模态。
设置计算内容：计算各模态的参与因子:在输出窗口中可以看到参与因子的计算结果：求解：10. 在一般后处理中，查看相对位移和应力 (载荷步 3)。
首先查看 Results Summary：载荷步 3 为 set 17查看载荷步 3 的相对位移和应力：在 Read Results 中读取载荷步 3 的结果：使用 By set Number，输入 Data Set Number 为 17。
Ansys nCode疲劳分析及应用
© 2015 ANSYS, Inc.
30
3.2 缺口应力集中对疲劳寿命的影响
理论分析：
与材料相关的参数，可通过下式估算：
q：缺口敏感系数
© 2015 ANSYS, Inc.
缺口半径
Kf：疲劳缺口系数 光滑件疲劳极限与缺口件疲劳极限之比
Kt：弹性应力集中系数 缺口处最大实际应力与该处名义应力之比
一般Kf小于Kt ，两者的关系可以用缺口敏 感系数q表征。
40
3.4 多轴状态对疲劳寿命的影响
如果加载过程中主应力方向改变，计算应力范围使用一个派生的单轴应力可 能是不正确的
临界平面法可用于计算主应力方向改变的情况
临界平面法计算平面应力和在几个径向平面进行重新求解
− 默认是每10度 (18 个平面) − 雨流计数法在每个平面上计算疲劳损伤 − 准则平面是最大损伤的平面
• 拉伸平均应力减少疲劳寿命
• 压缩平均应力增加疲劳寿命
Smax
Sa
Sm
DS
Sa Smin
© 2015 ANSYS, Inc.
36
3.3 平均应力对疲劳寿命的影响
对于应力-疲劳寿命，平均应力的影响通常使用Goodman 或者 Gerber平均应力 修正理论
Goodman 平均应力修正理论假设一条直线
1. 避免设计缺陷引起的疲劳破坏
2. 避免过于保守的设计
产品出现不应当发生 的疲劳失效，会使企 业的信誉受损，经济 损失更大！
过于保守设计，使得产 品的成本增加，市场竞 争力下降！
© 2015 ANSYS, Inc.
17
ANSYS nCode 的功能特色
支持的有限元结果： 静态分析(线性/非线 性) 瞬态分析 模态分析 频谱响应
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”为每个频率分配一个积木节点。
振动疲劳—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 nCode疲劳分析
© 2015 ANSYS, Inc.
14
Ansys nCode疲劳分析优势（二）
ANSYS nCode DesignLife是ANSYS与HBM的nCode强强 联合的结晶。Ncode公司是国际著名的疲劳耐久性工程专 业公司和技术领导者。自80 年代公司成立至今，引导并推 动了疲劳理论在工业领域中的应用及其发展。
© 2015 ANSYS, Inc.
39
3.4 多轴状态对疲劳寿命的影响
von-Mises 应力 标记的 von-Mises 应力
− 以最大绝对主应力的符号标记
最大剪切应力 最大主应力
− 用于脆性材料
绝对最大主应力 任一分量的应力 临界平面法
标记的 Von-Mises 应力是最常用的
© 2015 ANSYS, Inc.
平滑区域
扰动载荷 局部性 发展过程
5
疲劳分析方法
应力疲劳分析（S-N）
第一种疲劳分析方法 循环应力水平低、寿命长
应变疲劳分析（E-N）
循环应力水平高、寿命短 屈服后应变变化大、应力变化小
高周疲劳：材料在低于其屈服强度 的循环应力作用下，经104-105以上 循环产生的失效；
低周疲劳：材料在接近或超过其屈服 强度的应力作用下，低于104-105次塑 性应变循环产生的失效；
© 2015 ANSYS, Inc.
30
3.2 缺口应力集中对疲劳寿命的影响
理论分析：
与材料相关的参数，可通过下式估算：
q：缺口敏感系数
© 2015 ANSYS, Inc.
缺口半径
Kf：疲劳缺口系数 光滑件疲劳极限与缺口件疲劳极限之比
Kt：弹性应力集中系数 缺口处最大实际应力与该处名义应力之比
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第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 E
11 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.01
15.2 建立模型
跟所有在ANSYS6.1中进行的结构分析一样，建立本实例的有限元模型同样需要完成
如下工作：指定分析标题，定义材料性能，定义单元类型，定义单元实常数，建立几何模
型并进行有限元网格划分等。
对于本实例由于其结构中有大量的相同部件，因此主要介绍
如何利用ANSYS6.1提供的实体拷贝功能来建立需要的有限元模型，并大量用运实体属性
进行选择实体选取，读者可对这两种技巧着重掌握。
下面将详细讲解分析过程。
15.2.1指定分析标题并设置分析范畴
根据本实例分析的问题，指定分析标题为“PSD analysis of the beam - shell structure”。
因为本实例进行的仍是结构分析，所以应指定为结构分析范畴：Structural，以便在进行分
析时程序能提供合适的菜单选项。
下面进行具体的操作。
1．选取菜单路径Utility Menu | File | Change Jobname，将弹出Change Jobname (修改文
件名)对话框，如图15.2所示。
图15.2 修改文件名对话框
2．在Enter new jobname (输入新文件名)文本框中输入文字“CH15”，为本分析实例
的数据库文件名。
单击对话框中的按钮，完成文件名的修改。
3．选取菜单路径Utility Menu | File | Change Title，将弹出Change Title (修改标题)对话
框，如图15.3所示。
图15.3 修改标题对话框
4．在Enter new title (输入新标题)文本框中输入文字“PSD analysis of the beam - shell structure”，为本分析实例的标题名。
单击对话框中的按钮，完成对标题名的指定。
5．选取菜单路径Utility Menu | Plot | Replot，指定的标题“PSD analysis of the beam - shell structure”将显示在图形窗口的左下角。
6．选取菜单路径Main Menu | Preference，将弹出Preference of GUI Filtering (菜单过滤参数选择)对话框。
单击Structural(结构)选项使之被选中，以将菜单设置为与结构分析相关的选项。
单击按钮，完成分析范畴的指定。
15.2.2 定义单元类型
对于本实例分析的问题需要定义两种单元类型，一种为模拟梁的弹性梁单元BEAM4，另一种为模拟板特性的壳单元SHELL63。
具体的定义过程如下：
1．选取菜单路径Main Menu | Preprocessor | Element Type | Add/Edit/Delete，将弹出
Element Types (单元类型定义)对话框。
单击对话框中的按钮，将弹出Library of Element Types (单元类型库)对话框，如图15.4所示。
图15.4 单元类型库对话框
2．在左边的滚动框中单击“Structural Shell”，选择结构壳单元类型。
在右边的滚动框中单击“Elastic 4node 63”，使其高亮度显示，选择4节点弹性壳单元。
在对话框中单击按钮，完成对这种单元的定义。
由于单击了按钮，单元类型库对话框将继续显示。
3．接着继续在Library of Element Types (单元类型库)对话框的左边滚动框中单击“Structural Beam”，在右边的滚动框中单击“3D elastic 4”，使其高亮度显示，选择3维弹性梁单元。
单击对话框中的按钮，完成单元定义并关闭Library of Element Types (单元类型库)对话框。
4．在单元类型定义对话框的定义的单元列表框中显示了刚才定义的两种单元类型：SHELL63,BEAM4 (见图15.5)。
单击Element Types (单元类型定义)对话框中的按钮，关闭对话框中，完成单元类型的定义。
图15.5 定义的单元类型
15.2.3 定义单元实常数
因为使用的单元是梁单元和壳单元，所以还需要定义相应的单元实常数才能完成对单元特性的描述。
具体的操作如下：
1．选取菜单途径Main Menu | Preprocessor | Real Constants，将弹出Real Constants (实常数定义)对话框。
单击对话框中的按钮，将弹出Element Type for Real Constants (选择定义实常数的单元类型)对话框，如图15.6所示。
图15.6选择定义实常数的单元类型对话框
3．在选择单元类型列表框中，单击“Type 1 SHELL63”使其高亮度显示，选择第一类单元SHELL63。
然后单击该对话框中的按钮，将弹出Real Constant Set Number 1,for SHELL63 (为SHELL63单元定义实常数) 对话框如图15.07所示。
图15.7为SHELL63单元定义实常数对话框
4．在对话框中的Shell thickness at node I TK(I) (壳的厚度)文本框中输入2E-3，定义板壳的厚度为2E-3 m。
5．其余参数保持缺省。
单击按钮，关闭Real Constants Set Number 1，for SHELL63 (单元SHELL63的实常数定义)对话框。
完成对单元SHELL63实常数的定义。
6．重复步骤2的过程，在弹出的Element Type for Real Constants (选择定义实常数的单元类型)对话框的列表框中单击“Type 2 BEAM4”，使其高亮度显示。
然后单击按钮，将弹出Real Constant Set Number 2,for BEAM4 (为BEAM4单元定义实常数) 对话框如图15.8所示。