CN112487679A - 一种音箱声学特性的数值仿真方法、系统、终端及介质 - Google Patents
一种音箱声学特性的数值仿真方法、系统、终端及介质 Download PDFInfo
- Publication number
- CN112487679A CN112487679A CN202011313927.2A CN202011313927A CN112487679A CN 112487679 A CN112487679 A CN 112487679A CN 202011313927 A CN202011313927 A CN 202011313927A CN 112487679 A CN112487679 A CN 112487679A
- Authority
- CN
- China
- Prior art keywords
- sound
- sound box
- acoustic
- field
- loudspeaker
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000004088 simulation Methods 0.000 claims abstract description 56
- 238000010586 diagram Methods 0.000 claims abstract description 42
- 238000010168 coupling process Methods 0.000 claims abstract description 28
- 238000004458 analytical method Methods 0.000 claims abstract description 27
- 238000005859 coupling reaction Methods 0.000 claims abstract description 27
- 230000008878 coupling Effects 0.000 claims abstract description 26
- 239000000463 material Substances 0.000 claims abstract description 14
- 238000011160 research Methods 0.000 claims abstract description 8
- 238000000926 separation method Methods 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 230000005611 electricity Effects 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 5
- 229920000742 Cotton Polymers 0.000 claims description 4
- 239000007787 solid Substances 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000004140 cleaning Methods 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 9
- 238000011161 development Methods 0.000 abstract description 7
- 230000008569 process Effects 0.000 abstract description 6
- 230000035945 sensitivity Effects 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 208000005647 Mumps Diseases 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 208000010805 mumps infectious disease Diseases 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Obtaining Desirable Characteristics In Audible-Bandwidth Transducers (AREA)
Abstract
本发明公开了一种音箱声学特性的数值仿真方法、系统、终端及介质，涉及音箱技术领域，其技术方案要点是：获取扬声器的集总参数；建立音箱的三维几何模型图；根据三维几何模型图、材料属性参数以及耦合参数建立音箱的数值仿真模型；将数值仿真模型导入求解器中，根据分析问题的不同选择不同的研究类型；并根据模型大小选择全耦合或分离，以及选择直接或迭代算法进行求解；通过求解器计算得到数据集，并根据分析需要选择不同维度的绘图组和解，得到音箱的声学特性数值。本发明对音箱声学特性如灵敏度、指向性、箱体结构力学特性等的分析，缩短了开发周期，降低开发成本，且能对音箱工作过程物理过程有更清晰的认识，有利于提高音箱设计的品质。
Description
技术领域
本发明涉及音箱技术领域，更具体地说，它涉及一种音箱声学特性的数值仿真方法、系统、终端及介质。
背景技术
目前，音箱设计的的方法大多仍采用传统的经验法，一般为样品试做、测试、再改善样品、再测试的反复循环过程，这种方式开发周期长、成本高，且对测试场所测试硬件等均有较高要求。随着计算机技术的发展，使用数值仿真分析方法进行辅助设计变得越来越广泛，将数值仿真分析方法应用于音箱设计可以在样品试做前期预估音箱的性能，从而加快音箱的设计进度，减少开发成本。
现已有企业利用数值仿真方法对单独的扬声器、单独的箱体或声场进行分析，但对于音箱，其工作过程涉及多个物理场和多个系统，包括扬声器系统、箱体固体力学、声场等，其重要的声学特性参数也是多物理场相互耦合得到的结果，故在对音箱进行数值仿真分析时，须综合考虑各系统的相互耦合才能使仿真精度更高。
然而，对于音箱，特别是专业音箱产品，往往尺寸都较大且产品不具有轴对称性，导致无法对仿真模型进行对称简化。在进行多物理场耦合数值计算仿真时，若方法选择不当或处理方式不妥，则会导致对仿真模型进行数值计算时因计算能力不足而无法求解,。因此，如何研究设计一种音箱声学特性的数值仿真方法、系统、终端及介质是我们目前急需解决的问题。
发明内容
为解决现有音箱设计中对音箱组成部件单独分析，不能真实反应音箱产品整体声学特性，以及对音箱产品多物理场耦合数值计算时，易出现的因计算量超大而无法求解的问题，本发明的目的是提供一种音箱声学特性的数值仿真方法、系统、终端及介质。
本发明的上述技术目的是通过以下技术方案得以实现的：
第一方面，提供了一种音箱声学特性的数值仿真方法，包括以下步骤：
S101：获取扬声器的集总参数；
S102：通过三维绘图软件中或在数值仿真分析软件建立音箱的三维几何模型图；
S103：根据三维几何模型图、材料属性参数以及耦合参数建立音箱的数值仿真模型；
S104：将数值仿真模型导入求解器中，根据分析问题的不同选择不同的研究类型；并根据模型大小选择全耦合或分离，以及选择直接或迭代算法进行求解；
S105：通过求解器计算得到数据集，并根据分析需要选择不同维度的绘图组和解，得到音箱的声学特性数值。
进一步的，所述集总参数包括直流电阻Re、LR2电感损失模型中理想电感Le、LR2电感损失模型中电感L2、LR2电感损失模型中电阻R2、振动质量Mms、振动质量(不含推动的空气)Mmd、机械阻抗Rms、柔度Cms、磁路因数BL，获取具体为：
通过文字识别直接读取扬声器供应商提供的集中参数值；
或，通过扬声器测试设备测量得到集总参数。
进一步的，所述数值仿真模型建立具体为：
将三维几何模型图传输至数值计算软件中，对三维几何模型图进行网格划分，并清理三维几何模型图中多余的点、线、面、体等几何构造，同时生成包裹住音箱的空气域；
根据材料属性参数对三维几何模型图中各组成部分的材料属性进行修改或添加；
根据音箱扬声器、音箱箱体、声场及其相互耦合在三维几何模型图的仿真环境中建立与设置物理场、物理场的边界、载荷。
进一步的，所述扬声器物理场通过导纳类比法整合磁路系统、振动系统；
振动系统的计算方程为：
Fcoil＝FR+FM+FK
FR＝Rms·v
FM＝Mms·a＝Mms·j2πfv
式中，Fcoil为音圈和磁路系统总推动力，FR为机械和声学元件的摩檫力，FM为振动质量惯性力，FK为定心支片和折环的恢复力；Rms为机械阻抗，v为速度，Mms为振动质量，a为加速度，f为频率，j为虚数单位，Kms为刚性，x为位移；
Fcoil是流量变量，类似于电学中的电流；速度v为强度变量，类似于电学中的电压；则电路元件的力、速度、电压和电流的关系式具体为：
FM＝jωMms.v，
FR＝Rms v
I＝jωCU，
由导纳类比，扬声器由等效电路来表示，式中ω为角频率，I为电流，U为电压，C为电容，L为电感，R为电阻；
通过电力声类比，在数值计算软件中建立对应电路即可设置为扬声器的物理场；为电路各元件赋值，工作电压V0根据扬声器阻抗值匹配，其余参数由集总参数得到。
进一步的，所述音箱箱体物理场为固体力学物理场，建立音箱的力学方程为：
则简谐一频域为：
具体到音箱物理场为：
式中，ρ为音箱箱体密度，ω为箱体振动角频率，u为音箱振动位移矢量，S为箱体应力张量，Fv为箱体所受体积力(声介质中的压力作为载荷)，
为相位，kZ为面外波数，λ为特征值；
音箱箱体设置为线弹性模型，并指定杨氏模量、泊松比和密度。
进一步的，所述声场的物理场具体为：
通过在近场加上空气域包裹音箱，并通过有限元方法分析近场声场的压力声学频域；近场压力声学频域方程由连续流体的控制方程考虑无粘性和热损耗且无流动的情况，可得：
pt＝p+pb
选取设定声场内部声源，此声源由扬声器产生；并选取设定空气密度、声速、内部硬声场边界；
对音箱中铺设的吸音棉采用多孔声学Delany-Bazley-Mike模型，设置流阻率即可；
通过边界元方法分析远场声场的压力声学，声学方程为：
pt＝p2+pb
选取设定空气密度、声速；
式中，ρc为流体密度(此处为空气密度)，Pt为总声压，p和p2表示环境压力的声学变化(声学扰动)，pb为背景压力，qd为偶极源，ω为声波角频率，Cc为介质中的声速(此处介质为空气)，Qm为单极源，Keq为声波波数包含普通波数k和平面外波数kz和azimuthal波数km。
进一步的，所述多物理场耦合具体为：
根据扬声器、音箱箱体及声场相互作用的物理原理，声场与箱体用声-结构耦合，方程如下：
外部：
内部：
式中，utt是结构加速度，n是表面法线方向，pt为总声压，qd是声偶极子域声源，FA是作用在结构上的载荷(每单位面积的力)，up、down标识是指内部边界的两侧；
近场声场与远场声场用声学边界元BEM-有限元FEM耦合，方程如下：
pt，bem＝pt，fem
式中，Pt，bem为边界元边界总声压，Pt，fem为有限元边界总声压；
指定耦合场之间的边界。
第二方面，提供了一种音箱声学特性的数值仿真系统，包括：
数据获取模块，用于获取扬声器的集总参数；
图型建立模块，用于通过三维绘图软件中或在数值仿真分析软件建立音箱的三维几何模型图；
模型建立模块，用于根据三维几何模型图、材料属性参数以及耦合参数建立音箱的数值仿真模型；
数值求解模块，用于将数值仿真模型导入求解器中，根据分析问题的不同选择不同的研究类型；并根据模型大小选择全耦合或分离，以及选择直接或迭代算法进行求解；
结果分析模块，用于通过求解器计算得到数据集，并根据分析需要选择不同维度的绘图组和解，得到音箱的声学特性数值。
第三方面，提供了一种计算机终端，包含存储器、处理器及存储在存储器并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如第一方面中任意一项所述的一种音箱声学特性的数值仿真方法。
第四方面，提供了一种计算机可读介质，其上存储有计算机程序，所述计算机程序被处理器执行可实现如第一方面中任意一项所述的一种音箱声学特性的数值仿真方法。
与现有技术相比，本发明具有以下有益效果：
1、本发明通过数值仿真技术实现了对音箱声学特性如灵敏度、指向性、箱体结构力学特性等的分析，缩短了音箱开发设计的周期，降低了开发成本，且能对音箱工作过程物理过程有更清晰的认识，有利于提高音箱设计的品质；
2、本发明在建立耦合仿真模型时，将扬声器通过电力声类比，简化为电路物理场分析，省去了求解计算时大量的计算量，通常扬声器的集中参数易于得到，而扬声器磁路系统、振动系统的材料参数却不易得到，故采用电类比电路更易于仿真的实现；
3、本发明将组成音箱的扬声器、音箱箱体、吸音棉等与声场相互耦合进行数值仿真，较全面真实地考虑了影响音箱声学特性的各个方面，能更精确地分析音箱声学特性；
4、本发明在分析音箱声场时，使用了有限元FEM和边界元BEM耦合的方式，这样避免了只使用有限元进行数值仿真时，音箱体积过大而无法计算的难题。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：
图1是本发明实施例中的流程图；
图2是本发明实施例中三维几何模型图的示意图；
图3是本发明实施例中数值仿真模型的示意图；
图4是本发明实施例中音箱的频响曲线图；
图5是本发明实施例中音箱的辐射方向图
图6是本发明实施例中音箱箱体的位移变形示意图
图7是本发明实施例中扬声器的类比电路原理图；
图8是本发明实施例中的系统架构图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图1-8，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。
实施例1：一种音箱声学特性的数值仿真方法，如图1所示。
(一)获取扬声器的集总参数；集总参数包括Re、Le、L2、R2、Mms、Mmd、Rms、Cms、BL；集总参数通过文字识别直接读取扬声器供应商提供的集中参数值，或通过扬声器测试设备测量得到集总参数。现以一款音箱为例，集总参数如表1所示。
表1集总参数
(二)如图2所示，根据输入的集总参数通过三维绘图软件中或在数值仿真分析软件建立音箱的三维几何模型图；三维绘图软件包括但不限于ProE、Solidworks、UG软件。
(三)根据三维几何模型图、材料属性参数以及耦合参数建立音箱的数值仿真模型。如图3所示，数值仿真模型建立具体为：将三维几何模型图传输至数值计算软件COMSOL中，对三维几何模型图进行网格划分，并清理三维几何模型图中多余的点、线、面、体等几何构造，同时生成包裹住音箱的空气域；根据材料属性参数对三维几何模型图中各组成部分的材料属性进行修改或添加；根据音箱扬声器、音箱箱体、声场及其相互耦合在三维几何模型图的仿真环境中建立与设置物理场、物理场的边界、载荷。其中，网格划分为：指定网格单元类型及网格大小生成有限元网格单元，指定单元大小时，需注意单元大小需适配分析声波的最短波长，最大单元不大于最小波长的五分之一。对数值仿真模型划分网格，网格最大单元为CO/fmax/5，其中CO为声速，fmax为分析的最高频率。
(1)扬声器物理场通过导纳类比法整合磁路系统、振动系统。
振动系统的计算方程为：
Fcoil＝FR+FM+FK
FR＝Rms·v
FM＝Mms·a＝Mms·j2πfv
式中，Fcoil为总推动力，FR为摩檫力，FM为惯性力，FK为恢复力；Rms为机械阻抗，v为速度，Mms为振动质量，a为加速度，f为频率，j为虚数单位，Kms为刚性，x为位移。
Fcoil是流量变量，类似于电学中的电流；速度v为强度变量，类似于电学中的电压；则电路元件的力、速度、电压和电流的关系式具体为：
FM＝jωMms.v，
FR＝Rms v
I＝jωCU，
如图7所示，通过电力声类比，在数值计算软件中建立对应电路即可设置为扬声器的物理场；为电路各元件赋值，V0根据扬声器阻抗值匹配，其余参数由集总参数得到。如表2所示。
表2各物理场初始值及载荷
(2)音箱箱体物理场为固体力学物理场，建立音箱的力学方程为：
则简谐-频域为：
具体到音箱物理场为：
式中，ρ为音箱箱体密度，ω为箱体振动角频率，u为音箱振动位移矢量，S为箱体应力张量，Fv为箱体所受体积力，
为相位，kZ为面外波数，λ为特征值。
音箱箱体设置为线弹性模型，并指定杨氏模量、泊松比、密度、结构化损耗因子，如表2所示。
(3)声场的物理场具体为：通过在近场加上空气域包裹音箱，并通过有限元方法分析近场声场的压力声学频域；近场压力声学频域方程由连续流体的控制方程考虑无粘性和热损耗且无流动的情况，可得：
pt＝p+pb
选取设定声场内部声源，此声源由扬声器产生，设定声源为内部法向速度，其值为cir.R3_i。并选取设定空气密度、声速、内部硬声场边界。
对音箱中铺设的吸音棉采用多孔声学Delany-Bazley-Mike模型，设定流阻率Rf为20000Pa*s/m2。
通过边界元方法分析远场声场的压力声学，声学方程为：
pt＝p2+pb
选取设定空气密度、声速；式中，ρc为流体密度，Pt为总声压，p和p2表示环境压力的声学变化，pb为背景压力，qd为偶极源，ω为声波角频率，Cc为介质中的声速，Qm为单极源，Keq为声波波数包含普通波数k和平面外波数kz和azimuthal波数km。
(4)多物理场耦合根据扬声器、音箱箱体及声场相互作用的物理原理，声场与箱体用声-结构耦合，方程如下：
外部：
内部：
式中，utt是结构加速度，n是表面法线方向，pt为总声压，qd是声偶极子域声源，FA是作用在结构上的载荷，up、down标识是指内部边界的两侧。
近场声场与远场声场用声学边界元BEM-有限元FEM耦合，方程如下：
pt，bem＝pt，fem
指定耦合场之间的边界。
(四)将数值仿真模型导入求解器中，根据分析问题的不同选择不同的研究类型；并根据模型大小选择全耦合或分离，以及选择直接或迭代算法进行求解；为避免出现计算不收敛问题，可优选全耦合和直接求解器，如消耗内存过大，可考虑分离加迭代的方式。以进行频率研究为例，设定计算的频率点，此列采用1/12Oct，选择全耦合方式，使用直接求解器MUMPS进行求解。
(五)通过求解器计算得到数据集，并根据分析需要选择不同维度的绘图组和解，得到音箱的声学特性数值；声学特性数值包括但不限于声压分布、灵敏度、声压等值面、指向性参数、箱体的应力、位移。因远场声场分析使用了边界元，故数据集须根据分析需求建立三维格栅。在音箱正前方1m处取一三维截点分析音箱的灵敏度，其频响曲线如图4所示，建立极坐标绘图组，绘制音箱前方垂直平面和水平面的辐射方向图，如图5所示。建立三维绘图组分析音箱箱体的位移变形图，如图6所示。根据分析需求的不同还可以建立声压分布图、声压级分布图、声压等值面图、箱体应力图等。
实施例2
一种音箱声学特性的数值仿真系统，如图8所示，包括：数据获取模块，用于获取扬声器的集总参数；图型建立模块，用于根据输入的集总参数通过三维绘图软件中或在数值仿真分析软件建立音箱的三维几何模型图；模型建立模块，用于根据三维几何模型图、材料属性参数以及耦合参数建立音箱的数值仿真模型；数值求解模块，用于将数值仿真模型导入求解器中，根据分析问题的不同选择不同的研究类型；并根据模型大小选择全耦合或分离，以及选择直接或迭代算法进行求解；结果分析模块，用于通过求解器计算得到数据集，并根据分析需要选择不同维度的绘图组和解，得到音箱的声学特性数值。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种音箱声学特性的数值仿真方法，其特征是，包括以下步骤：
S101：获取扬声器的集总参数；
S102：通过三维绘图软件中或在数值仿真分析软件建立音箱的三维几何模型图；
S103：根据三维几何模型图、材料属性参数以及耦合参数建立音箱的数值仿真模型；
S104：将数值仿真模型导入求解器中，根据分析问题的不同选择不同的研究类型；并根据模型大小选择全耦合或分离，以及选择直接或迭代算法进行求解；
S105：通过求解器计算得到数据集，并根据分析需要选择不同维度的绘图组和解，得到音箱的声学特性数值。
2.根据权利要求1所述的一种音箱声学特性的数值仿真方法，其特征是，所述集总参数包括直流电阻Re、LR2电感损失模型中理想电感Le、LR2电感损失模型中电感L2、LR2电感损失模型中电阻R2、振动质量Mms、振动质量Mmd、机械阻抗Rms、柔度Cms、磁路因数BL，获取具体为：
通过文字识别直接读取扬声器供应商提供的集中参数值；
或，通过扬声器测试设备测量得到集总参数。
3.根据权利要求1所述的一种音箱声学特性的数值仿真方法，其特征是，所述数值仿真模型建立具体为：
将三维几何模型图传输至数值计算软件中，对三维几何模型图进行网格划分，并清理三维几何模型图中多余的点、线、面、体等几何构造，同时生成包裹住音箱的空气域；
根据材料属性参数对三维几何模型图中各组成部分的材料属性进行修改或添加；
根据音箱扬声器、音箱箱体、声场及其相互耦合在三维几何模型图的仿真环境中建立与设置物理场、物理场的边界、载荷。
4.根据权利要求3所述的一种音箱声学特性的数值仿真方法，其特征是，所述扬声器物理场通过导纳类比法整合磁路系统、振动系统；
振动系统的计算方程为：
Fcoil＝FR+FM+FK
FR＝Rms·v
FM＝Mms·a＝Mms·j2πfv
式中，Fcoil为音圈和磁路系统总推动力，FR为机械和声学元件的摩檫力，FM为振动质量惯性力，FK为定心支片和折环的恢复力；Rms为机械阻抗，v为速度，Mms为振动质量，a为加速度，f为频率，j为虚数单位，Kms为刚性，x为位移；
Fcoil是流量变量，类似于电学中的电流；速度v为强度变量，类似于电学中的电压；则电路元件的力、速度、电压和电流的关系式具体为：
FM＝jωMms.v，
I＝JωCU，
由导纳类比，扬声器由等效电路来表示，式中ω为角频率，I为电流，U为电压，C为电容，L为电感，R为电阻；
通过电力声类比，在数值计算软件中建立对应电路即可设置为扬声器的物理场；为电路各元件赋值，工作电压V0根据扬声器阻抗值匹配，其余参数由集总参数得到。
5.根据权利要求4所述的一种音箱声学特性的数值仿真方法，其特征是，所述音箱箱体物理场为固体力学物理场，建立音箱的力学方程为：
则简谐—频域为：
具体到音箱物理场为：
式中，ρ为音箱箱体密度，ω为箱体振动角频率，u为音箱振动位移矢量，S为箱体应力张量，Fv为箱体所受体积力，
为相位，kZ为面外波数，λ为特征值；
音箱箱体设置为线弹性模型，并指定杨氏模量、泊松比和密度。
6.根据权利要求5所述的一种音箱声学特性的数值仿真方法，其特征是，所述声场的物理场具体为：
通过在近场加上空气域包裹音箱，并通过有限元方法分析近场声场的压力声学频域；近场压力声学频域方程由连续流体的控制方程考虑无粘性和热损耗且无流动的情况，可得：
pt＝p+pb
选取设定声场内部声源，此声源由扬声器产生；并选取设定空气密度、声速、内部硬声场边界；
对音箱中铺设的吸音棉采用多孔声学Delany-Bazley-Mike模型，设置流阻率即可；
通过边界元方法分析远场声场的压力声学，声学方程为：
pt＝p2+pb
选取设定空气密度、声速；
式中，ρc为流体密度，Pt为总声压，p和p2表示环境压力的声学变化，pb为背景压力，qd为偶极源，ω为声波角频率，Cc为介质中的声速，Qm为单极源，Keq为声波波数包含普通波数k和平面外波数kz和azimuthal波数km。
7.根据权利要求6所述的一种音箱声学特性的数值仿真方法，其特征是，所述多物理场耦合具体为：
根据扬声器、音箱箱体及声场相互作用的物理原理，声场与箱体用声-结构耦合，方程如下：
外部：
内部：
式中，utt是结构加速度，n是表面法线方向，pt为总声压，qd是声偶极子域声源，FA是作用在结构上的载荷，up、down标识是指内部边界的两侧；
近场声场与远场声场用声学边界元BEM—有限元FEM耦合，方程如下：
pt，bem＝pt，fem
式中，Pt,bem为边界元边界总声压，Pt,fem为有限元边界总声压；
指定耦合场之间的边界。
8.一种音箱声学特性的数值仿真系统，其特征是，包括：
数据获取模块，用于获取扬声器的集总参数；
图型建立模块，用于通过三维绘图软件中或在数值仿真分析软件建立音箱的三维几何模型图；
模型建立模块，用于根据三维几何模型图、材料属性参数以及耦合参数建立音箱的数值仿真模型；
数值求解模块，用于将数值仿真模型导入求解器中，根据分析问题的不同选择不同的研究类型；并根据模型大小选择全耦合或分离，以及选择直接或迭代算法进行求解；
结果分析模块，用于通过求解器计算得到数据集，并根据分析需要选择不同维度的绘图组和解，得到音箱的声学特性数值。
9.一种计算机终端，其特征是，包含存储器、处理器及存储在存储器并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如权利要求1-7中任意一项所述的一种音箱声学特性的数值仿真方法。
10.一种计算机可读介质，其特征是，其上存储有计算机程序，所述计算机程序被处理器执行可实现如权利要求1-7中任意一项所述的一种音箱声学特性的数值仿真方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011313927.2A CN112487679A (zh) | 2020-11-20 | 2020-11-20 | 一种音箱声学特性的数值仿真方法、系统、终端及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011313927.2A CN112487679A (zh) | 2020-11-20 | 2020-11-20 | 一种音箱声学特性的数值仿真方法、系统、终端及介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112487679A true CN112487679A (zh) | 2021-03-12 |
Family
ID=74932863
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011313927.2A Pending CN112487679A (zh) | 2020-11-20 | 2020-11-20 | 一种音箱声学特性的数值仿真方法、系统、终端及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112487679A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113286233A (zh) * | 2021-05-20 | 2021-08-20 | 维沃移动通信有限公司 | 扬声器仿真方法、仿真装置、扬声器以及电子设备 |
| CN113630686A (zh) * | 2021-08-13 | 2021-11-09 | 南京工程学院 | 一种基于模式识别的高强度Helmholtz声源设计方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102004823A (zh) * | 2010-11-11 | 2011-04-06 | 浙江中科电声研发中心 | 一种扬声器振动和声学特性的数值模拟方法 |
| CN102368281A (zh) * | 2011-11-14 | 2012-03-07 | 浙江中科电声研发中心 | 一种扬声器磁路系统的数值模拟方法 |
| CN104252559A (zh) * | 2014-08-29 | 2014-12-31 | 浙江中科电声研发中心 | 一种扬声器多场耦合的数值仿真分析方法 |
| CN104408227A (zh) * | 2014-10-28 | 2015-03-11 | 浙江中科电声研发中心 | 一种扬声器失真特性的数值仿真分析方法 |
| CN108460204A (zh) * | 2018-02-26 | 2018-08-28 | 浙江中科电声研发中心 | 一种通过扬声器振动部件的受力和位移反推其材料动态力学参数的方法 |
| CN109145514A (zh) * | 2018-09-30 | 2019-01-04 | 浙江中科电声研发中心 | 一种扬声器失真的数值仿真分析方法 |
| CN110188428A (zh) * | 2019-05-20 | 2019-08-30 | 武汉大学 | 一种基于有限元法的扬声器联合仿真方法 |
| CN110427650A (zh) * | 2019-07-02 | 2019-11-08 | 浙江中科电声研发中心 | 动铁式扬声器基本特性的数值仿真分析方法 |
| CN110442907A (zh) * | 2019-07-02 | 2019-11-12 | 浙江中科电声研发中心 | 压电式mems扬声器基本特性的数值仿真分析方法 |
-
2020
- 2020-11-20 CN CN202011313927.2A patent/CN112487679A/zh active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102004823A (zh) * | 2010-11-11 | 2011-04-06 | 浙江中科电声研发中心 | 一种扬声器振动和声学特性的数值模拟方法 |
| CN102368281A (zh) * | 2011-11-14 | 2012-03-07 | 浙江中科电声研发中心 | 一种扬声器磁路系统的数值模拟方法 |
| CN104252559A (zh) * | 2014-08-29 | 2014-12-31 | 浙江中科电声研发中心 | 一种扬声器多场耦合的数值仿真分析方法 |
| CN104408227A (zh) * | 2014-10-28 | 2015-03-11 | 浙江中科电声研发中心 | 一种扬声器失真特性的数值仿真分析方法 |
| CN108460204A (zh) * | 2018-02-26 | 2018-08-28 | 浙江中科电声研发中心 | 一种通过扬声器振动部件的受力和位移反推其材料动态力学参数的方法 |
| CN109145514A (zh) * | 2018-09-30 | 2019-01-04 | 浙江中科电声研发中心 | 一种扬声器失真的数值仿真分析方法 |
| CN110188428A (zh) * | 2019-05-20 | 2019-08-30 | 武汉大学 | 一种基于有限元法的扬声器联合仿真方法 |
| CN110427650A (zh) * | 2019-07-02 | 2019-11-08 | 浙江中科电声研发中心 | 动铁式扬声器基本特性的数值仿真分析方法 |
| CN110442907A (zh) * | 2019-07-02 | 2019-11-12 | 浙江中科电声研发中心 | 压电式mems扬声器基本特性的数值仿真分析方法 |
Non-Patent Citations (5)
| Title |
|---|
| JINLAN HUANG: "在 COMSOL® 中模拟声-结构相互作用", pages 1 - 11, Retrieved from the Internet <URL:https://cn.comsol.com/blogs/introduction-to-modeling-acoustic-structure-interactions-in-comsol/> * |
| SØREN VEDEL: "用于声学仿真的新阻抗边界条件", pages 1 - 8, Retrieved from the Internet <URL:https://cn.comsol.com/blogs/new-impedance-boundary-conditions-for-acoustics-simulations/> * |
| THOMAS FORRISTER: "使用贝塞尔面板基准模型分析扬声器阵列", pages 1 - 7, Retrieved from the Internet <URL:https://cn.comsol.com/blogs/analyzing-a-loudspeaker-array-with-a-bessel-panel-benchmark-model/> * |
| ZHANGZI0231: "电声基础-第一章 电力声类比", pages 1 - 38, Retrieved from the Internet <URL:https://book.duxiu.com/searchDoc?sw=%E7%94%B5%E5%8A%9B%E5%A3%B0%E7%B1%BB%E6%AF%94&allsw=&searchtype=&channel=searchDoc&bCon=&ecode=utf-8&Field=all> * |
| ZHANGZI0231: "电声基础-第一章 电力声类比", 独秀, 7 August 2011 (2011-08-07), pages 1 - 38 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113286233A (zh) * | 2021-05-20 | 2021-08-20 | 维沃移动通信有限公司 | 扬声器仿真方法、仿真装置、扬声器以及电子设备 |
| CN113630686A (zh) * | 2021-08-13 | 2021-11-09 | 南京工程学院 | 一种基于模式识别的高强度Helmholtz声源设计方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110427650B (zh) | 动铁式扬声器基本特性的数值仿真分析方法 | |
| Liu et al. | A point interpolation mesh free method for static and frequency analysis of two-dimensional piezoelectric structures | |
| Jin et al. | A meshless method for some inverse problems associated with the Helmholtz equation | |
| CN102970647B (zh) | 一种扬声器振动中非线性特性的仿真计算方法 | |
| Amir et al. | Approximate reanalysis in topology optimization | |
| CN111106783B (zh) | 一种信号制作方法、信号制作装置、振动马达及触屏设备 | |
| CN112487679A (zh) | 一种音箱声学特性的数值仿真方法、系统、终端及介质 | |
| CN104408227A (zh) | 一种扬声器失真特性的数值仿真分析方法 | |
| CN111079278B (zh) | 三维时域杂交间断伽辽金方法外加电磁源项的处理方法 | |
| Zhang et al. | Meshless regular hybrid boundary node method | |
| Caplan et al. | Advancing parabolic operators in thermodynamic MHD models: Explicit super time-stepping versus implicit schemes with Krylov solvers | |
| Gao et al. | An immersed free-surface boundary treatment for seismic wave simulation | |
| Takeuchi et al. | Source directivity approximation for finite-difference time-domain simulation by estimating initial value | |
| TWI480522B (zh) | 電聲換能器之參數測量方法 | |
| Nielsen et al. | Estimation of optimal values for lumped elements in a finite element—Lumped parameter model of a loudspeaker | |
| CN105103155B (zh) | 分析装置以及分析方法 | |
| WO2020114088A1 (zh) | 一种扬声器振动系统的运动过程分析方法 | |
| Brunet et al. | State-space modeling and identification of loudspeaker with nonlinear distortion | |
| Dilgen et al. | Shape optimization of the time-harmonic response of vibroacoustic devices using cut elements | |
| CN108520110A (zh) | 电容式微加工超声传感器三类薄膜模型的有限元分析方法 | |
| JP7113168B2 (ja) | 音圧計算方法、および音圧計算プログラム | |
| CN105138746A (zh) | 通过扬声器振动部件的劲度系数反推其杨氏模量的方法 | |
| Rubinacci et al. | A fast algorithm for solving 3-D eddy current problems with integral formulations | |
| Wei et al. | Acceleration of free-vibrations analysis with the Dual Reciprocity BEM based on ℋ-matrices and CUDA | |
| Lee et al. | A three-dimensional meshless scheme with background grid for electrostatic-structural analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |