1、M3-1谐响应分析谐响应分析M3-2第一节：谐响应分析的定义和目的第一节：谐响应分析的定义和目的第二节：关于谐响应分析的基本术语和概念第二节：关于谐响应分析的基本术语和概念第三节：谐响应分析在第三节：谐响应分析在ANSYS中的应用中的应用第四节：谐响应分析的实例练习第四节：谐响应分析的实例练习M3-3什么是谐响应分析？什么是谐响应分析？确定一个结构在已知频率的正弦（简谐）载荷作用下结构响应的技术。确定一个结构在已知频率的正弦（简谐）载荷作用下结构响应的技术。输入：输入：已知大小和频率的谐波载荷（力、压力和强迫位移）；已知大小和频率的谐波载荷（力、压力和强迫位移）；同一频率的多种载荷，可以是同相
2、或不同相的。同一频率的多种载荷，可以是同相或不同相的。输出：输出：每一个自由度上的谐位移，通常和施加的载荷不同相；每一个自由度上的谐位移，通常和施加的载荷不同相；其它多种导出量，例如应力和应变等。其它多种导出量，例如应力和应变等。M3-4谐响应分析用于设计：谐响应分析用于设计：旋转设备（如压缩机、发动机、泵、涡轮机械等）的支座、固定旋转设备（如压缩机、发动机、泵、涡轮机械等）的支座、固定装置和部件；装置和部件；受涡流（流体的漩涡运动）影响的结构，例如涡轮叶片、飞机机受涡流（流体的漩涡运动）影响的结构，例如涡轮叶片、飞机机翼、桥和塔等。翼、桥和塔等。M3-5为什么要作谐响应分析？为什么要作谐响应
3、分析？确保一个给定的结构能经受住不同频率的各种正弦载荷（例如：以确保一个给定的结构能经受住不同频率的各种正弦载荷（例如：以不同速度运行的发动机）；不同速度运行的发动机）；探测共振响应，并在必要时避免其发生（例如：借助于阻尼器来避探测共振响应，并在必要时避免其发生（例如：借助于阻尼器来避免共振）。免共振）。M3-6包含的主题：包含的主题：运动方程运动方程谐波载荷的本性谐波载荷的本性复位移复位移求解方法求解方法M3-7通用运动方程：通用运动方程：F矩阵和矩阵和 u矩阵是简谐的，频率为矩阵是简谐的，频率为 w w:谐响应分析的运动方程：谐响应分析的运动方程：FuKuCuM titiititiieui
4、ueeuueFiFeeFFwwww)()(21max21max )()(21212FiFuiuKCiMwwM3-8Fmax=载荷幅值载荷幅值 =载荷函数的相位角载荷函数的相位角F1 =实部实部,Fmaxcos F2 =虚部虚部,Fmaxsin umax=位移幅值位移幅值u1=实部实部,umaxcosf fu2=虚部虚部,umaxsinf fM3-9 在已知频率下正弦变化；在已知频率下正弦变化；相角相角 允许不同相的多个载荷允许不同相的多个载荷同时作用，同时作用，缺省值为零；缺省值为零；施加的全部载荷施加的全部载荷都假设是简都假设是简谐的，包括温度和重力。谐的，包括温度和重力。实部虚部M3-10
5、在下列情况下计算出的位移将是复数在下列情况下计算出的位移将是复数 具有阻尼具有阻尼 施加载荷是复数载荷（例如：虚部为非零的载荷）施加载荷是复数载荷（例如：虚部为非零的载荷）复位移滞后一个相位角复位移滞后一个相位角（相对于某一个基准而言）（相对于某一个基准而言）可以用实部和虚部或振幅和相角的形式来查看可以用实部和虚部或振幅和相角的形式来查看M3-11求解简谐运动方程的三种方法：求解简谐运动方程的三种方法：完整法完整法为缺省方法，是最容易的方法；为缺省方法，是最容易的方法；使用完整的结构矩阵，且允许非对称矩阵（例如：声学矩阵）。使用完整的结构矩阵，且允许非对称矩阵（例如：声学矩阵）。缩减法缩减法*
6、使用缩减矩阵，比完整法更快；使用缩减矩阵，比完整法更快；需要选择主自由度，据主自由度得到近似的需要选择主自由度，据主自由度得到近似的 M矩阵和矩阵和C矩阵。矩阵。模态叠加法模态叠加法*从前面的模态分析中得到各模态；再求乘以系数的各模态之和；从前面的模态分析中得到各模态；再求乘以系数的各模态之和；所有求解方法中最快的。所有求解方法中最快的。*未包括在本手册中；请参看未包括在本手册中；请参看ANSYS结构分析指南。结构分析指南。*将在第六章中讨论。将在第六章中讨论。M3-12The procedure for a full harmonic response analysis consists o
7、f three main steps:Build the model.Apply loads and obtain the solution.Review the results.M3-13The procedure for a reduced harmonic analysis consists of five main steps:Build the model.Apply the loads and obtain the reduced solution.Review the results of the reduced solution.Expand the solution(expa
8、nsion pass).Review the results of the expanded solution.M3-14The procedure to use the method consists of five main steps:Build the model.Obtain the modal solution.Obtain the mode superposition harmonic solution.Expand the mode superposition solution.Review the results.M3-15M3-16四个主要步骤：四个主要步骤：建模建模选择分
9、析类型和选项选择分析类型和选项施加谐波载荷并求解施加谐波载荷并求解观看结果观看结果M3-17模型模型只能用于线性单元和材料，忽略各种非线性；只能用于线性单元和材料，忽略各种非线性；记住要输入密度；记住要输入密度；注意：注意：如果如果ALPX（热膨胀系数）和（热膨胀系数）和 T均不为零，就有可能不经意均不为零，就有可能不经意地包含了简谐热载荷。为了避免这种事情发生，请将地包含了简谐热载荷。为了避免这种事情发生，请将ALPX设置为设置为零零.如果参考温度如果参考温度 TREF与均匀节点温度与均匀节点温度 TUNIF不一致不一致,那么那么 T为非零值为非零值；请参阅第一章中的建模需要考虑的问题。请参
10、阅第一章中的建模需要考虑的问题。M3-18/PREP7ET,.MP,EX,.MP,DENS,!建立几何模型建立几何模型!划分网格划分网格.M3-193 建模建模选择分析类型和选项选择分析类型和选项输入求解器，选择谐响应分析；输入求解器，选择谐响应分析；主要分析选项是求解方法主要分析选项是求解方法-在后面讨在后面讨论；论；规定阻尼规定阻尼-在后面讨论。在后面讨论。典型命令典型命令:/SOLUANTYPE,HARMIC,NEWM3-20分析选项分析选项求解方法求解方法-完整法、缩减法和模态叠加法。完整法、缩减法和模态叠加法。缺省为完整法；缺省为完整法；自由度输出格式自由度输出格式-主要用于批处理方
11、式中；主要用于批处理方式中；集中质量矩阵。集中质量矩阵。推荐用于如果结构的一个方向的尺寸远推荐用于如果结构的一个方向的尺寸远小于另两个方向的尺寸的情况中。例如小于另两个方向的尺寸的情况中。例如：细长梁与薄壳。：细长梁与薄壳。典型命令典型命令:HROPT,.HROUT,LUMPM,.M3-21阻尼阻尼从从-阻尼、阻尼、-阻尼和阻尼比中选取阻尼和阻尼比中选取阻尼比最常用阻尼比最常用典型命令典型命令:ALPHAD,BETAD,DMPRAT,.M3-223建模建模3选择分析类型和选项选择分析类型和选项施加谐波载荷并求解施加谐波载荷并求解所有施加的载荷以规定的频率（或频率所有施加的载荷以规定的频率（或频
12、率范围）简谐地变化范围）简谐地变化“载荷载荷”包括：包括：位移约束位移约束-零或非零的零或非零的作用力作用力压强压强注意：注意：如果要施加重力和热载荷，它如果要施加重力和热载荷，它们也被当作简谐变化的载荷来考虑！们也被当作简谐变化的载荷来考虑！典型命令典型命令:DK,!或或 D或或DSYMDA,.DL,M3-23规定谐波载荷时要包括：规定谐波载荷时要包括：振幅和相角振幅和相角频率频率阶梯载荷或线性变化载荷的说明阶梯载荷或线性变化载荷的说明振幅和相角振幅和相角载荷值（大小）代表振幅载荷值（大小）代表振幅 Fmax相角相角 f f 是在两个或两个以上谐波载荷间的是在两个或两个以上谐波载荷间的相位差
13、，单一载荷不需要相角相位差，单一载荷不需要相角 f f 实部虚部fF1maxF2maxM3-24振幅和相角（接上页）振幅和相角（接上页）ANSYS 不能直接输入振幅和相角，而是规定实部不能直接输入振幅和相角，而是规定实部和虚部分量；和虚部分量；例如，施加两个简谐力例如，施加两个简谐力 F1和和 F2，其相角相差其相角相差 f f:F1real=F1max(F1的振幅)F1imag=0F2real=F2maxcosfF2imag=F2maxsinf可以使用可以使用APDL语言计算，但要确保角度单位为度语言计算，但要确保角度单位为度（缺省为弧度）（缺省为弧度）。实部虚部fF1maxF2maxM3-
14、25*AFUN,DEGFK,F,SFA,SFL,SFE,SF,M3-26M3-27谐波载荷的频率：谐波载荷的频率：通过频率范围和在频率范围内的子步数量来规通过频率范围和在频率范围内的子步数量来规定每秒的循环次数（赫兹）；定每秒的循环次数（赫兹）；例如，在例如，在0-50频率范围内有频率范围内有10个子步时将给出个子步时将给出在在5，10，15.45和和50Hz等频率上的解；而同等频率上的解；而同一频率范围只有一个子步时，则只给出一频率范围只有一个子步时，则只给出50Hz频频率上的解率上的解。典型命令典型命令:HARFRQ,NSUBST,KBC,1M3-28阶梯载荷阶梯载荷/线性变化载荷：线性变
15、化载荷：采用若干子步，可以逐渐地施加载荷（线性变化载荷），或者在第采用若干子步，可以逐渐地施加载荷（线性变化载荷），或者在第一个子步立刻施加载荷（阶梯载荷）；一个子步立刻施加载荷（阶梯载荷）；谐波载荷通常是阶梯加载，因为载荷值代表的是最大振幅。谐波载荷通常是阶梯加载，因为载荷值代表的是最大振幅。M3-29在施加谐波载荷后，下一步就是开始求解在施加谐波载荷后，下一步就是开始求解通常采用一个载荷步，但是可以采用若干子通常采用一个载荷步，但是可以采用若干子步，且每个子步具有不同的频率范围步，且每个子步具有不同的频率范围典型命令典型命令:HARFRQ,NSUBST,KBC,1M3-303健摸健摸3选择
16、分析类型和选项选择分析类型和选项3施加谐波载荷并求解施加谐波载荷并求解观看结果观看结果分三步分三步 绘制结构上的特殊点处的位移绘制结构上的特殊点处的位移-频率曲线频率曲线 确定各临界频率和相应的相角确定各临界频率和相应的相角 观看整个结构在各临界频率和相角时的位移观看整个结构在各临界频率和相角时的位移和应力和应力采用采用POST26，时程后处理器，时程后处理器采用采用POST1，通用后处理器，通用后处理器M3-31位移位移-频率关系曲线频率关系曲线首先定义首先定义 POST26 变量变量 节点和单元数据表节点和单元数据表 用大于等于二的数据识别用大于等于二的数据识别 变量变量1包含各频率，并是
17、预先定义了的包含各频率，并是预先定义了的M3-32定义变量（接上页）定义变量（接上页）挑选可能发生最大变形的节点，然后选择自由度的方向；挑选可能发生最大变形的节点，然后选择自由度的方向；定义变量的列表被更新。定义变量的列表被更新。M3-33定义变量定义变量画变量关系曲线画变量关系曲线典型命令典型命令:/POST26NSOL,PLVAR,.M3-34确定各临界频率和相角确定各临界频率和相角用图形显示最高振幅发生时的频率；用图形显示最高振幅发生时的频率；由于位移与施加的载荷不同步（如果存在阻尼的话），需要确定出现最大由于位移与施加的载荷不同步（如果存在阻尼的话），需要确定出现最大振幅时的相角；振幅
18、时的相角；要进行上述工作，首先要选择振幅要进行上述工作，首先要选择振幅+相位选项。相位选项。M3-35 然后用表列出变量。然后用表列出变量。注意：最大振幅注意：最大振幅=3.7出现在出现在48Hz,85.7时时下一步就是观看整个模型在该频率和相角下一步就是观看整个模型在该频率和相角下的位移和应力（使用下的位移和应力（使用POST1)典型命令典型命令:PRCPLX,1PRVAR,FINISHM3-36观看整个结构的结果观看整个结构的结果进入进入POST1，且列出结果综述表，确定临界频率的载荷步和子步序号，且列出结果综述表，确定临界频率的载荷步和子步序号；典型命令典型命令:/POST1SET,LI
19、STM3-37使用使用 HRCPLX 命令读入在期望频率和相角时的命令读入在期望频率和相角时的结果：结果：HRCPLX,LOADSTEP,SUBSTEP,PHASE,.例如：例如：HRCPLX,2,4,85.7绘制变形图，应力等值线图和其它期望的结果。绘制变形图，应力等值线图和其它期望的结果。典型命令典型命令:HRCPLX,PLDISP,2PLNSOL,FINISHM3-383 建立模型建立模型3 选择分析类型和选项选择分析类型和选项3 施加谐波载荷和求解施加谐波载荷和求解3 观看结果观看结果M3-39吉他弦的谐响应分析。弦收到吉他弦的谐响应分析。弦收到F1=84N的拉力，在弦的的拉力，在弦的1/4长度处受长度处受到到F2=1N的激励力。计算结构基频，并表明在奇数频率激励时结构的激励力。计算结构基频，并表明在奇数频率激励时结构中点才有位移响应。中点才有位移响应。
