15-谐响应分析3-215-谐响应分析

谐响应分析,本章内容,本章中，DS的谐响应分析将会涉及到以下内容:假定用户已经了解了第四章的线性静态结构分析和第五章的模态分析所涉及的内容.在本章中，将会涵盖如下内容:建立谐响应分析谐响应求解方法阻尼查看结果本章所叙述的功能适用于ANSYS专业licenses及更高的licenses.一些特例将会在相应的位置注明,谐响应分析基础,谐分析常用于确定频率一定下的稳定正弦（简谐）载荷作用下的结构响应.一个谐分析或一个频率响应分析仅考虑单一频率的条件下。载荷之间可以具有不同的相位，但需已知激振频率。这项操作不能用于任意的瞬态载荷上.用户在进行谐分析前，可以先对模型进行一次模态分析，以了解模型的动力学特性.为了更好地理解谐分析，下列为运动学的通用方程:,谐响应分析基础,在一个谐分析中，结构的载荷与响应被假定为简谐的（循环）使用复数符号能有效地表示响应的状况。由于ejA可以简单地等于（(cos(A)+jsin(A),其中有虚部项(j=-1)，这就表示带有相位差的正弦运动。激振频率W是指加载时产生的频率。如果几个不同相位的载荷同时发生激振，将会产生一个力相位变换y；如果存在阻尼或力的相位变换，将会产生一个位移相位变换f。,谐响应分析基础,例如，考虑如右图所示的两力共同作用在同一结构上的工况两力都有受到同一频率W激励。但是.,”Force2”滞后于“Force1”45度的相位差，“Force2”的相位角y度。以上的叙述可通过复数标记的方法表示。因此，可写成:使用这种方法，F1表示复数的实部，F2为虚部.响应x类似于F,ModelshownisfromasampleSolidWorksassembly.,谐响应分析基础,对于谐分析，复数响应x1与x2可从矩阵方程的求解中获得：其结果有如下假设:M,C,与K恒定不变的:材料假设为线弹性的小变形，并不存在非线性包含有阻尼矩阵C.但若是激振频率W与结构的固有频率w相同，响应将变得无限大.虽然有相位的存在，但载荷F(与响应x)仍是按给定的频率W作正弦变化注意：要特别注意在DS中与谐分析相关的假设.,A.谐响应分析步骤,谐分析的操作很类似于线性静态的操作，因此我们不会涉及所有详细的操作步骤。在下列的文字叙述中将说明谐分析的步骤.获取几何模型设置材料属性定义接触域（若可用的话）定义网格控制（可选）施加载荷与约束的条件指定所要求谐分析选项设置谐分析选项求解模型查看结果,几何模型,在谐分析中，可适用任何类型的几何模型适用实体、面、线及其任意的组合对于线，将不能输出应力与应变的结果可以加入PointMass虽然只有加速度载荷对它起作用,材料属性,在谐分析中，要求输入杨氏弹性模量，泊松比和密度其它所有材料的属性可以指定，但它们不会参与谐分析后面将说明，阻尼不是作为材料的属性输入，而是作为全局属性被输入,接触区域,接触区域可用在模态分析中。但由于在线性分析中的接触行为不同于非线性接触问题，如下表所示:接触行为类似于模态分析（见第五章）。在模态分析中，由于简谐模拟是线性的，非线性接触相对于它的线性对应部件做出了简化.一般在谐分析中不要使用非线性接触,载荷和约束,除以下情况之外，载荷与约束皆可适用于谐分析中不支持热载荷不支持角速度不支持平移力载荷由于螺栓预紧力载荷是非线性的，因此也不适用“CompressionOnlySupport”是非线性的，也不适用。若施加了这种约束，它也变得类似于“FrictionlessSupport”记住：所有结构载荷都在同一激振频率作用下成正弦变化,载荷和约束,约束载荷表如下所示:在下一节中，我们将讨论“求解的方法”.要注意的是，在ANSYS专业licenses中并不支持完全法求解，因此它也并不支持在谐分析中给定位移约束.并不是所有可用的载荷皆可输入相位。螺栓载荷与力矩载荷的相位角为0.若有其它的载荷（如螺栓载荷、力矩载荷）存在，其相位角仍保持0。,载荷和约束,添加谐分析载荷:选择分析类型为“Harmonic”。输入载荷数据（矢量或者分量形式）。输入适当的相位角。若已知载荷的实部F1与虚部F2，模的大小与相位y便可根据如下的公式计算：,载荷和约束,选择载荷，可显示出载荷两个周期的图像，然后点击“Worksheet”菜单。模的大小与相位角便在载荷的视图中显示出来,B.谐响应分析求解,在求解前，需设定谐分析工具选项:选择求解分支条，并从相关的工具条中插入一个谐分析工具在谐分析工具的明细窗中，用户能通过输入最大值、最小值来确定激振频率域，并确定求解的步长。频率域fmax-fmin，间隔数n决定频率的步长DWDS将从W+DW.开始，求解n个频率,求解选项,在ANSYSStructural及其更高的licenses中，有两种求解方法。这两种方法各有优缺点，因此接下来针对它们进行讨论:模态叠加法是默认的求解选项，适用于ANSYS专业及其更高的licenses完全法适用于ANSYS专业及其更高的licenses在谐分析工具的明细窗中，“SolutionMethod”栏中只有两个的选项可以勾选（若可用的话）求解命令条的明细窗的选项不可用，因为它对分析并没有影响.,模态叠加法,模态叠加法是在模态的坐标中求解谐分析方程的谐分析方程如下对于线性系统，用户可以将x写成关于模态形状的fi的线性组合的表达式：在关系式中，yi指模态的坐标（系数）。例如，用户可以通过求解一个模态分析来确定固有频率wi和相应的模态形状因子fi。可以看到，包括的模态n越多，对x逼近越精确。,模态叠加法,前面的讨论主要介绍了关于模态叠加法的一些基本信息。在此，有三点需要注意的：1.由于采用了模态的坐标系，因此使用模态叠加法进行谐分析时，程序会首先自动地进行模态分析。这些过程能很清楚的记录在求解命令条的工作表菜单中，也能传给用户。DS能自动地确定模态数n，这对精确求解来说是必要的。虽然首先进行的是模态分析，但正如上个幻灯片所示，谐分析部分的求解还是很迅速且高效的，因此，总的来说，模态叠加法通常比完全法要快的多。,模态叠加法,由于进行了模态分析，DS将会获得结构的自然频率。在谐分析中，响应的峰值是与结构的固有频率相对应。由于自然频率已知，DS能够将结果聚敛到自然振动频率附近而不是使用evenlyspaced结果。,模态叠加法,3.由于模态叠加法的本质原因，所以并不支持给定位移的约束。由于求解在模态坐标系下完成的，所以不允许有非零的位移。在讨论载荷与约束期间，这早已被提及。,完全法,完全法也是求解谐分析的一种方法。谐分析的方程如下：在完全法中，直接在节点坐标系下求解矩阵方程，除了使用了复数基本类似于线性静态分析。,完全法,通过与模态叠加法几个不同的结果作比较：1.对每一个频率，完全法必须将Kc因式分解。在模态叠加法中是求解化简后的非耦合方程；在完全法中，必须将复杂的耦合矩阵KC因式分解。因此，完全法一般比模态叠加法更耗计算时间。2.支持给定位移约束由于对x直接求解，允许施加位移约束，并可以使用给定位移约束。,完全法,3.完全法没有使用模态的信息与模态叠加法不同的，完全法并不依赖模态形状与固有频率程序在内部并不执行模态分析对xC的求解是精确的。并没有模态形状的响应x近似的结果产生但是，由于在求解过程中，DS并没有产生模态信息，因此，不会产生频率聚集的结果，只有频率均匀分布的结果产生。,C.阻尼的输入,在谐方程中有一个阻尼矩阵C如前所述，阻尼是被指定为全局属性的。在ANSYS专业licenses中,只能输入阻尼比常量x在ANSYS结构licenses中可以输入阻尼比常量x或beta阻尼值注：两种阻尼选项同时被输入，影响是累积的。两种阻尼选项皆可在两种求解方法（完全法和模态叠加法）中使用。,阻尼的基础知识,动力系统中阻尼的结果是产生能量损耗对响应有影响的阻尼会改变结构固有频率并减弱响应的幅值。在任何的结构系统中，阻尼都是以多种形式存在的。由于多种因素的影响，阻尼是一个复杂的现象。然而，数学上对阻尼的表达却是很简单的，粘性阻尼被认为：粘性阻尼力Fdamp是与速度成正比的，其中c为阻尼常数。临界阻尼ccr是指没有振动时的阻尼值。阻尼比率x是指实际阻尼c与临界阻尼ccr的比值。,恒定阻尼比率,在DesignSimulation中，输入恒定的阻尼比率意味着x值在整个频率区间内是一个常量。x值在模态叠加法中会被直接使用。阻尼比率常数x是无量纲的在完全法中，阻尼率x不会被直接使用。在程序内部，它根据C转化为一个合适的值。,Beta阻尼,另一种模拟阻尼的方法是假设阻尼C是与刚度K成正比的，其中b为常量:那么相关的阻尼比率x为:由方程可知，在beta阻尼的作用下，阻尼对频率的影响成线性增长。与阻尼比率是常数不同，beta阻尼是随频率的增大而增大的。Beta阻尼往往衰减掉高频的影响。Beta阻尼的单位是时间。,Beta阻尼,有两种方法可输入Beta阻尼:直接输入Beta阻尼输入阻尼比率与频率，相应的Beta阻尼值，将会由DS中根据上一张幻灯片中的公式计算出来。,阻尼相互关系,一般来说，还有其它一些计算阻尼的方法。这些方法通常用于单自由度系统，因此，在将它们扩展到多自由度系统（例如有限元）中时，要慎重。品质因子Qi是1/(2xi)损耗因子hi是负的Q或2xi对数衰减量di近似于空阻尼情况，约为2pxi半量带宽Dwi近似于带有轻微阻尼的结构，约为2wixi这些计算阻尼的方程都进行了简化，仅适用于单自由度系统。如果用户能很好地理解在频率时间内的物理结构的响应与阻尼比率常数跟beta阻尼的差异，那么在DS中将能很好地模拟阻尼的情况。,D.谐响应结果获取工具,可以从谐响应工具条获取结果。可获取三种类型的结果:一定频率和相位角下，装配体、零件、表面的位移、应力、应变的结果云图。选定的点、线、面处的频率响应图，包括加速度、位移、应力、应变的最小值、最大值、平均值。特定频率下的零件的位移、应力、应变平均最大值、最小值相位响应。与线性静力分析不同，在获取结果之前，必须启动求解。否则，如果当一个求解完成时，想获取其它结果就必须重新进行求解。,谐响应结果获取工具,通过谐响应工具条，可以获取任意有用的结果。确认实体查询结果的范围。对于边和表面，要确定平均值、最小值、最大值是否显示。输入其它适合的参数。如果要查询的结果在求解的频率之间，可用线性插值计算响应。例如，如果在100和1000Hz之间以100Hz的间隔求解，若用户想获得333Hz处的结果,可以在300和400Hz之间进行线性插值。,RequestHarmonicToolResults,DesignSimulation假定响应是谐函数(正弦).如果分量不同相，获得的量，象等效应力/主应力不是简谐的。因此这些结果并不可用。收敛控制并不能用在谐分析结果中。执行模态分析，并在反映响应的模态形状上执行收敛。这有助于确认网格是否足够密以达到在后续的谐分析中捕捉动态响应。,求解模型,在谐分析中并不用到求解明细窗中的选项。只有将要被求解的分析类型的状态信息会显示出来.在设定谐响应分析选项与需求解的结果后，通常点击“solve”按钮，就可进行求解。,云图显示结果,在给定频率与相位角的条件下，可获得应力、应变或位移分量的结果云图,云图动画显示,这些结果可以制作成动画播放。这些动画使用到了实际的谐响应（实部与虚部的结果）,频率响应图,可以显示：应力、应变、位移或加速度的分量的图象。,频率响应图,在频率给定的条件下，能绘出带输入力载荷的应力、应变、或位移各分量的相位对比图。,RequestingResults,一个谐响应分析的求解通常需要进行多次求解：使用频率搜索器的模态分析总是首先执