什么是频率响应?
一般而言,系统的频率响应显示系统的某些属性对以频率形式输入的激励产生响应的函数。COMSOL Multiphysics中的频率响应,通常是指对谐波激励的线性(或线性化)响应。为了生成频率响应曲线,需要进行频率扫描; 也就是对很多不同频率进行求解。通常,频率响应曲线将表现出许多与系统固有频率对应的不同峰值。![62eaa284-271e-eb11-8da9-e4434bdf6706.png](http://p01.5ceimg.com/content/62eaa284-271e-eb11-8da9-e4434bdf6706.png)
重新审视单自由度系统
阻尼的固有频率为 该频率是指在没有其他外部激励的情况下,将系统从变形状态释放之后,产生振动(具有衰减幅度)的频率。此时,一个有趣的问题出现了:“哪个激励频率会给出最大幅度的响应?”一般会认为恰好是阻尼的固有频率,但如下图所示,情况并非如此。![65eaa284-271e-eb11-8da9-e4434bdf6706.jpeg](http://p03.5ceimg.com/content/65eaa284-271e-eb11-8da9-e4434bdf6706.jpeg)
由于系统是谐波运动,因此采用复数表示法。分解出常见的谐波乘数 ,可得运动方程为
载荷f的相位角可以作为参考,因此f是实值。除以刚度k可以得到归一化形式:
等式的右边是静态位移。因此,动态解和静态解之比是
H 函数通常称作传递函数。其中 β 表示激励频率与无阻尼固有频率的比值。传递函数的大小为。因此,给出最大响应的激励频率为
最大响应的激励频率低于阻尼固有频率, 实际上,频移接近两倍。 最大响应的激励频率与自由频率不同,这似乎是矛盾的,但可以解释为由阻尼引起的力与位移之间的相移所导致。 在无阻尼的情况下,载荷和位移在低于固有频率时完全同相,高于 固有频率时翻转到180°异相。 而有了阻尼之后,相移是平滑过渡,如下图所示。 无论阻尼水平如何,无阻尼固有频率的相移始终为90°位移的相移作为频率的函数关系。
当存在阻尼时,由于力可能向系统提供能量,因此造成力和位移产生相移。
损耗因子的阻尼
对具有阻尼损耗因子的单自由度系统进行分析,此时系统的运动方程为阻尼的固有频率为
令人惊讶的是,增加阻尼并不会使固有频率减小,而是增加。这是因为这种形式的损耗因子阻尼实际上会增加刚度。复数刚度的绝对值是
由于这种损耗因子阻尼,传递函数为
其大小为
可以看出,当β= 1时(即,无阻尼固有频率),振幅最大。同样,最大响应频率低于有阻尼时的固有频率。
前一篇博客文章中提到,损耗因子阻尼的另一种定义具有如下属性:复数刚度的绝对值与阻尼水平无关具有以下属性:复杂刚度的绝对值与阻尼水平无关。这可以使用对复数刚度进行归一化的定义方法来获得如下复平面的纯旋转。
¥$
由上述公式可知,固有频率随着阻尼而减小:
本文不做详细分析,结论是相应的激励频率下降,将给出最大振幅,因此仍低于有阻尼的固有频率。当考虑损耗因子阻尼时,激励和响应之间的相移非常有趣:即使在非常低的激励频率下,相移仍然存在,且其渐近于arctan(η)
当使用损耗因子阻尼时,位移的相移与频率的函数关系。低频渐近线用虚线表示。
当两个表面之间的摩擦提供阻尼机制时,由于系统的非线性,对谐波输入的响应不再是谐波,可能仍然会存在周期性,但并不是谐波响应。频域法无法解决这些问题,因为其假设是输入输出为线性关系。
在 COMSOL 软件中进行频率响应建模
建立研究
在模型向导中添加结构力学物理接口后,您将看到许多研究类型,其中四种可用于计算频率响应:- 频域
- 频域,预应力
- 频域,模态
- 频域,预应力,模态
![6deaa284-271e-eb11-8da9-e4434bdf6706.png](http://p02.5ceimg.com/content/6deaa284-271e-eb11-8da9-e4434bdf6706.png)
![6eeaa284-271e-eb11-8da9-e4434bdf6706.png](http://p02.5ceimg.com/content/6eeaa284-271e-eb11-8da9-e4434bdf6706.png)
是否存在扰动?
对于频域研究,求解器序列中的稳态节点中有一个非常重要的设置,即线性。![6feaa284-271e-eb11-8da9-e4434bdf6706.png](http://p03.5ceimg.com/content/6feaa284-271e-eb11-8da9-e4434bdf6706.png)
![72eaa284-271e-eb11-8da9-e4434bdf6706.png](http://p03.5ceimg.com/content/72eaa284-271e-eb11-8da9-e4434bdf6706.png)
解释结果
频域分析的结果是复值,其中隐含了谐波变化。复数的相位角描述了与参考相位之间的相移(可以任意选择,但通常作为主载荷的相位)。它还提供有关结构中不同点之间相移的信息。需要注意的是,由于单个有限元内的位移分量可以具有不同的相位角,因此应力张量的分量也很可能彼此不同相。这在疲劳分析中会比较重要。在许多情况下,如在彩色图中,只能显示实数。所有展示的结果约定如下:如果期望查看一个复值变量v的实数值时,则使用实数部分。 其中,数据集的属性(如相角Φ)是可以修改的。
![73eaa284-271e-eb11-8da9-e4434bdf6706.png](http://p05.5ceimg.com/content/73eaa284-271e-eb11-8da9-e4434bdf6706.png)
![75eaa284-271e-eb11-8da9-e4434bdf6706.png](http://p02.5ceimg.com/content/75eaa284-271e-eb11-8da9-e4434bdf6706.png)
为了更详细地了解发生的情况,我们可以将结果的虚部和辐角添加到图中:
![76eaa284-271e-eb11-8da9-e4434bdf6706.png](http://p02.5ceimg.com/content/76eaa284-271e-eb11-8da9-e4434bdf6706.png)
在低频时,实部接近绝对值。在固有频率附近,则主要是虚部。这意味着响应与激励几乎不同相。考虑一下,如果将数据集里面的相角改为45°,频率响应会怎样变化。
![78eaa284-271e-eb11-8da9-e4434bdf6706.png](http://p01.5ceimg.com/content/78eaa284-271e-eb11-8da9-e4434bdf6706.png)
由上图可知,正如预期,振幅不会改变。但实部和虚部的各个值均会改变,其中相角曲线向上移动π/ 4。实际上,如果向载荷添加45°相位角,我们将获得完全相同的结果图。
![7beaa284-271e-eb11-8da9-e4434bdf6706.png](http://p01.5ceimg.com/content/7beaa284-271e-eb11-8da9-e4434bdf6706.png)
在模型中,除了上述相角输入,您还可以使用复数表示法直接输入载荷:
![7eeaa284-271e-eb11-8da9-e4434bdf6706.png](http://p03.5ceimg.com/content/7eeaa284-271e-eb11-8da9-e4434bdf6706.png)
当所有载荷并非彼此同相时,能够规定相角非常重要。例如,可以通过使y方向上的载荷相对于x方向上的载荷具有90°相移,来方便地描述一个旋转的非平衡受力物体。
扰动研究的结果
如果研究是扰动类型,则会有两组结果:预应力解和扰动解。在这种情况下,可以在各种结果表示功能中访问额外的选择::待求表达式。![80eaa284-271e-eb11-8da9-e4434bdf6706.png](http://p04.5ceimg.com/content/80eaa284-271e-eb11-8da9-e4434bdf6706.png)
在这里,您可以选择研究扰动解,预应力解或其组合。对于扰动解,您还可以获得一个选项:计算微分复选框。
![82eaa284-271e-eb11-8da9-e4434bdf6706.png](http://p04.5ceimg.com/content/82eaa284-271e-eb11-8da9-e4434bdf6706.png)
此设置会影响非线性表达式的处理方式。如果未选择计算微分,则会在面值处获取非线性数量。例如,u^2将会很方便地从扰动解中获取变量u的平方。由于u通常是复值的,因此通常该操作是无意义的。
当选择计算微分时,则该非线性量将预应力状态附近线性化。表达式u^2将计算2*u0*u,其中u0是线性化点处的值。将频率响应结果转换为时域
在某些情况下,您可能希望在时域中进行频域分析的谐波响应可视化,尤其希望有多个激励频率。![83eaa284-271e-eb11-8da9-e4434bdf6706.png](http://p01.5ceimg.com/content/83eaa284-271e-eb11-8da9-e4434bdf6706.png)
您可以使用频域到时域 FFT研究步骤将频率响应结果转换到时域。
![85eaa284-271e-eb11-8da9-e4434bdf6706.png](http://p02.5ceimg.com/content/85eaa284-271e-eb11-8da9-e4434bdf6706.png)
此技术用于以下教程模型:
- 厚梁的振动分析
- 黏弹性结构阻尼器 – 瞬态分析
结束语
频域分析是分析受谐波激励的线性系统的有力工具。实际上,通过对负载进行初始傅里叶变换,可以使用频率响应分析来研究任何类型的周期性激励。 COMSOL案例库中还有更多机械频率响应分析示例,例如:- 弯头支架(完整版)
- 支架(bracket_frequency.mph模型)
- 复合材料层合板的强迫振动分析
- 复合压电传感器
以上模拟案例,可以通过以下方式进行查找:首先进入COMSOL 中国的官网,然后在首页的右上角的查找框里输入您想要查找的关键字,如:弯头支架。就可以查到相关的所有资料了。
推荐阅读
COMSOL 免费开放 15 个专题培训视频资料评估壳厚度对消声器性能的影响CFD 仿真讨论——有限元法(FEM)VS. 有限体积法(FVM)焊接中的结构疲劳分析地震中建筑物的稳定性分析本文内容来自 COMSOL 博客,点击“”,阅读更多延伸文章。