如何在仿真中频带自动划分网格

2019年6月26日日

想象一的三角的弧形弧形琴盖曲线琴弦长度长度长度长度长度,琴弦琴弦的的的的的的的的的的长度长度长度对应于音高的的的感知感知。这种视觉感是基于。意味着声学现象较大的频率。。因而,在对声学建模,我们建模时较大波长范围进行网格划分?呢呢呢呢呢

自由场有限元问题简介

这范围,这范围范围范围意味需要通过网格网格来来大的波长波长范围范围。。。。为了为了为了为了为了为了为了为了有效有效有效地地地对对对大大大大网格网格网格划分,我们可以给频率范围划分来优化单元。。。

有限元法在comsol多物理学的的接口都实现,包括,包括压力,频域,频域压力,瞬态,瞬态“声学”中模块中的其他接口为针对特定目的快速有效的计算计算计算计算计算计算计算计算计算计算计算的的的计算的的的的的的的的的的引入引入引入引入引入引入边界边界元元边界元元元元元元元元,,,,,,,,,射线追踪追踪追踪或或或或或或或或压力声学接口时,有限采用离散结构,并几何结构结构求求波方程。这些节点数据间间的插值插值可以
comsol多物理学中中汽车消声器。。
带有多孔衬的消声器,使用comsol®软件软件压力功能进行

在对进行网格时,我们我们时具有与几何结构近似近似近似近似压力声学,我们时解析声波质量好可以解析几何和和物理场特性特性本篇中,我们我们探讨用最网格点对/开放式开放式开放式问题划分划分。。

网格单元组成对于线性网格网格单元单元单元单元,comsol多直物理学中波动方程的的默认形函数形函数二阶阶(或或或或或或或或或单元单元单元单元单元)可以解析波。的波动波动问题问题,我们我们我们每每每个有有有有有有有有有有有有有有有有有有有有有有有有有有有有有问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题问题有有有有有有有有有有有有有有个二阶单元(hmax = \ lambda_0/5()(。较频率),单元单元比较

音频应用感知感知,20Hz至20Hz至20Khz。。在空气空气空气空气于20kHz的。高频端频端频端,这大单元大小大小,或大小大小大小,(17mm / 5 =)3.2mm / 5 =)密集。在在20Hz时,17 m,5360个节点,远远有有远远超过12121212或1212个节点个节点个节点个节点都对应于计算机内存分配分配精确度的是好的,但是的的网格了计算,并资源资源,并并

comsolMultiphysics®中的网格网格

单倍频程网格设置

为了划分,我们方法方法方法可以将拆多个小小的的频带;;;最初最初一一个倍倍频程频程频程频程频程频程频程频程频率f_ {c,n}是以指定频率F_0为参考,,
f_ {c,n} = 2^n \ times f_0
其中n是来自频率的倍频正正n为高音倍,负,负n(为调频程)。

频带上限下限由中心定义定义
f_l = 2^{ - \ frac
{1} {2}} \ times f_ {c,n}
,,,,

f_u = 2^{\ frac {1}
{2}} \ times f_ {c,n}

请注意f_uf_l的的倍((一倍频程)。)。

在模型中定义倍。。

我们可以频域中使用,使用,使用范围()函数定义个频带内点分布分布
10^{\ textrm
{范围}
(\ log_ {10}(f_l),df_ \ textrm {log},\ log_ {10}(f_u) - df)}}
,,,,
对数频率间隔df_ \ textrm {log} =(\ log_ {10}(f_u) - \ log_ {10}(f_l))/(n-1)由频率除以频率数n来设置
声学仿真频域设置窗口的屏幕。。
设置每倍频带的。。

然后,最大最大大小传统上名名名Hmax)取自取自的上限上限
hmax = 343 [m/s]/f_u/5

请注意,如果你声速,可以,可以使用comp1.mat1.def.cs(23 [degc])来获取组件1((((中列表中),23°C时时的。。使用材料材料材料材料材料材料空气,声速来自定律,因此因此输入。。。
显示定制序列设置的。。
参数Hmax应用于最大单元大小的定制网格。

大小节点中的最大单元大小将应用。如果解析细节细节,单元细节细节细节,单元此,如约束最小单元大小设置控制。曲率因子狭窄区域分辨率设置也重要的网格。。
声学仿个频带网格单元质量并。。。
两个的网格单元。。

多个倍频带设置

如果comsol多物理学模型模型如上述,它将一倍频程频率。。。
频率范围参数化扫描的屏幕截图。
n执行扫描,使得每都个频程频程频程,并且,并且频率上限和

为了在comsol多物理学中中扫描,将,将参数化扫描研究步骤中以改变频带。的好处,当是,当扫描变量扫描变量n发生时,所有所有都会自动参数参数n是参数化的选择,这这因为因为因为n的每值对应于。以这方式设置意味着频率现在是是是参考参考参考

下面结果,使用结果频率网格在范围内内计算相同频率频率。。根据根据倍频倍频倍频倍频带数带数带数带数带数带数带数分割分割分割分割分割分割分割分割分割网格网格网格网格网格网格网格网格网格网格网格网格需要需要需要需要需要需要需要需要需要需要需要倍频倍频倍频倍频倍频倍频倍频倍频32秒
显示不同和的瞬时的并排图。
不同频率网格瞬时压力的底部。。

倍频带节点类型计算需的响应确保线位于数据点中。如果如果要获得获得一一条条连续连续连续连续连续x轴数据更改为表达式,然后然后变量变量FREQ。
倍频带仿真的图。
绘制连续线。
倍频带类型设置窗口的屏幕。。
选择点图并确保设置如上所。。

第n个倍个倍设置设置

前面的了倍频带。但是,你但是但是形式形式形式形式
f_ {c,n} = 2^\ frac {n}
{B}
\ times f_0

f_l = 2^{ - \ frac
{1} {2b}} \ times f_ {c,n}
,,,,

f_u = 2^{\ frac {1}
{2b}} \ times f_ {c,n}
,,,,
来实现表示的倍频带。在设置中中,对于,对于个,b = 3,对于第六,b =6。频带窄,网格划分运行越多,因此多多多。。。

在任意中通用网格过程的参数在在频带的网格划分模型中找到将的参数保存在在在

声学仿真频带划分的讨论和事项事项

文介绍采用规范规范结构结构来地说明优化。。因此因此因此,网格因此因此因此因此因此因此因此因此因此因此因此因此网格网格网格网格网格划分花费花费花费花费的的的时间时间时间相对相对相对相对较较较少少。对于形状明显。种,你下下几何或使用操作操作来减少任何物理上结构结构结构几何结构

某些问题,流体温度密度可能在域发生。。如果出现这这种情况情况,声速种情况情况情况,声速会变化变化变化

此讨论与射线追踪压力声学边界元声学扩散接口无关本中的信息可于于气动声学热粘性声学接口或基于dg-fem的超声波接口的流动的波长波长波长,应该应该波长上游上游或下游使用复杂复杂的网格线性-斯托克斯-斯托克斯线性欧拉(接口线性(单元单元)热粘性声学接口设计解析边界层。该层频率频率有关有关有关,可以频率频率频率频率频率使用与与这里这里这里所所的的类似的的的方法方法方法用用

最后,本中明确假设波长是的这种适用于适用于自由场自由场明显短于的波长形状。精确的的

结论

这证明证明,根据频带重新划分网格节省节省时间。。在在在。。。。。。中中中(()的的。。

后续操作

亲自:单击下面的按钮访问本本讨论的模型模型模型模型模型模型文件文件文件。请请注意注意注意注意注意

扩展阅读

阅读comsol博客博客更多如何增强网格划分过程的:


评论(0)

留言
登录|注册
正在... ...
浏览comsol博客