[go: up one dir, main page]

CN102004823B - Numerical value simulation method of vibration and acoustic characteristics of speaker - Google Patents

Numerical value simulation method of vibration and acoustic characteristics of speaker Download PDF

Info

Publication number
CN102004823B
CN102004823B CN201010538802A CN201010538802A CN102004823B CN 102004823 B CN102004823 B CN 102004823B CN 201010538802 A CN201010538802 A CN 201010538802A CN 201010538802 A CN201010538802 A CN 201010538802A CN 102004823 B CN102004823 B CN 102004823B
Authority
CN
China
Prior art keywords
loudspeaker
finite element
vibration
model
element model
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.)
Active
Application number
CN201010538802A
Other languages
Chinese (zh)
Other versions
CN102004823A (en
Inventor
温周斌
李军
金向锋
徐楚林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiashan Yiendi Electroacoustic Technology Service Co., Ltd.
Zhejiang Electro-Acoustic R&D Center, CAS
Original Assignee
ZHEJIANG ELECTRO-ACOUSTIC R&D CENTER CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by ZHEJIANG ELECTRO-ACOUSTIC R&D CENTER CAS filed Critical ZHEJIANG ELECTRO-ACOUSTIC R&D CENTER CAS
Priority to CN201010538802A priority Critical patent/CN102004823B/en
Publication of CN102004823A publication Critical patent/CN102004823A/en
Application granted granted Critical
Publication of CN102004823B publication Critical patent/CN102004823B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Stereophonic System (AREA)

Abstract

The invention provides a numerical value simulation method of the vibration and acoustic characteristics of a speaker. The method comprises the following steps of: firstly, drawing a 3D (Three-Dimensional) geometric model diagram of the speaker by using 3D drawing software; then, adding the 3D geometric model diagram of the speaker to meshing software to mesh the 3D geometric model diagram into body elements, defining element types, materials and boundary conditions and applying a load to acquire a finite element model; establishing a boundary element model matched with the finite element model; and finally, solving the finite element model with a finite element solver to acquire the vibration characteristics of the speaker, and solving the boundary element model with a boundary element solver to acquire the acoustic characteristics of the speaker, wherein the vibration characteristics include the natural frequency, the vibration mode (vibration type), the displacement, the strain and the stress, and the acoustic characteristics include a frequency response curve and a directivity curve.

Description

The method for numerical simulation of a kind of loudspeaker vibration and acoustic characteristic
Technical field
The invention belongs to the loudspeaker field, relate to the analogy method of loudspeaker vibration and acoustic characteristic, particularly relate to the method for numerical simulation of a kind of loudspeaker vibration and acoustic characteristic.Adopt the vibration and the acoustic characteristic of this methods analyst loudspeaker, can access natural frequency, mode of oscillation (vibration shape), displacement, strain, stress, frequency response curve, the directive property curve of loudspeaker.Can be widely used in the speaker design exploitation.
Background technology
Traditional loudspeaker analysis design method has two kinds: (1) empirical method.Be generally sample and have a fling at, test, improve sample, the cyclic process repeatedly of test again, this mode must design the later stage by the time could find the loudspeaker problem, and the construction cycle is long, cost is high; (2) equivalent circuit method.When speaker operation during at low frequency, the diaphragm of loudspeaker is done whole piston motion, and at this moment loudspeaker can be regarded lumped-parameter system as; According to acousto-electro-mechanical analogue, each parts equivalence is the circuit component of lumped parameter, can obtain analysis result rapidly and accurately; But when frequency is higher; Because the diaphragm of loudspeaker can produce and cut apart vibration, loudspeaker becomes distributed parameter system, can't analyze loudspeaker performance with equivalent electrical circuit.
Denmark LoudSoft company has developed loudspeaker design softwares such as FINEMotor, FINECone; Can analyze magnetic circuit of loudspeaker, voice coil loudspeaker voice coil, vibrating diaphragm; Obtain parameters such as T/S parameter, impedance curve, frequency response curve and directive property curve; But this two software institute analysis loudspeaker model is necessary for axisymmetric model, for racetrack loudspeaker and vibrating diaphragm surface the loudspeaker of reinforcement is arranged, and its limitation is very big.
Summary of the invention
The objective of the invention is to design the method for numerical simulation of a kind of loudspeaker vibration and acoustic characteristic.
The present invention will solve is the design later stage of must waiting until that empirical method exists could find loudspeaker problem and the construction cycle is long, cost the is high problem and the problem that when frequency is higher, can't analyze loudspeaker performance with equivalent electrical circuit of equivalent circuit method existence.
The vibration of a kind of loudspeaker of the present invention and the method for numerical simulation of acoustic characteristic, concrete steps are:
(1) use three-dimensional drawing software to draw three-dimensional (3D) the geometric model figure of loudspeaker.
(2) set up the loudspeaker finite element model.Concrete steps are following:
A, loudspeaker 3D geometric model figure is imported grid dividing software, it is carried out grid dividing, obtain the used body unit of finite element analysis.If need to consider the influence of air, also in model, add air and be body unit its grid dividing to loudspeaker vibration.
B, definition unit type.Define the cell type that is complementary for each several part unit in the A step.
C, definition material parameter.The material parameter of each parts of definition model comprises Young modulus, Poisson ratio, density and material damping.If comprise air in the model, also need define the material parameter of air, comprise the velocity of sound, density and acoustical absorption coefficient.
D, definition boundary condition.The boundary condition of loudspeaker finite element model is following: 1) fixed boundary condition.Because diaphragm of loudspeaker edge is fixed on the basin frame, need at its edge definition fixed boundary condition, if loudspeaker also includes locating branch slice, the edge that it and basin bridge joint touch also should define fixed boundary condition; 2) the solid coupled boundary condition of stream.If consider air influence when analyzing the vibration characteristics of loudspeaker, coupling can take place at the surface of contact of loudspeaker and air, need to flow coupled boundary condition admittedly in the interphase definition of loudspeaker and air.
E, give the finite element model imposed load, apply driving force at the voice coil loudspeaker voice coil position of model.
F, above-mentioned model is saved as the finite element model file.
(3) set up the loudspeaker boundary element model.The surface of extraction loudspeaker finite element model obtains the veil lattice with finite element model body unit coupling; Because baffle has very big influence to the acoustic characteristic of loudspeaker; Also to set up the veil lattice of baffle, define the cell type of above-mentioned veil lattice, above-mentioned model is saved as the boundary element model file.
(4) with the finite element solving device finite element model file is found the solution.The finite element model file is imported the finite element solving device and carries out model analysis, can obtain the natural frequency and the mode of oscillation (vibration shape) of loudspeaker, carry out frequency analysis again, can obtain displacement response, strain-responsive and the stress response of loudspeaker.
A, finite element model is carried out model analysis
The fundamental equation of dynamic analysis is following:
[ M ] { u · · } + [ C ] { u · } + [ K ] { u } = { F }
In the following formula; [M] is the mass of system matrix; [C] is the system damping matrix; [K] is the system stiffness matrix; is the node vector acceleration,
Figure BSA00000340545900023
be the node speed vector, { u} is the nodal displacement vector; { F} is the load vector that applies, with time correlation.
Suppose that loudspeaker does undamped-free vibration, equation is:
[ M ] { u · · } + [ K ] { u } = { 0 }
For linear system, suppose that system does simple harmonic motion
{u}={u} i?cosω it
The above-mentioned equation of motion becomes
([K]-ω i 2[M]){u} i={0}
Find the solution above-mentioned equation and can obtain ω i{ u} i, ω iBe the intrinsic circular frequency in i rank, { u} iBe proper vector, i.e. the corresponding mode of oscillation (vibration shape) of the intrinsic circular frequency in i rank.
Figure BSA00000340545900025
is i rank natural frequency.
B, finite element model is carried out frequency analysis
In frequency analysis, { F}, { u} is a simple harmonic quantity.Have:
{F}=({F 1}+i{F 2})e iΩt
{u}=({u 1}+i{u 2})e iΩt
([K]-Ω 2[M]+iΩ[C])({u 1}+i{u 2})={F 1}+i{F 2}
In the following formula, { u 1, { u 2Be respectively the real part and the imaginary part of displacement, { F 1, { F 2Being respectively the real part and the imaginary part of loading force, Ω is the angular frequency of loading force.
Find the solution above-mentioned equation and can obtain the real part { u of displacement 1And imaginary part { u 2, further can derive strain and stress by displacement.
(5) the finite element solving device is found the solution the displacement and the boundary element model file that obtain and be submitted to the boundary element solver, can find the solution the frequency response curve and the directive property curve that obtain loudspeaker.
In isotropic ideal fluid medium, for the sound wave that simple harmonic quantity changes, acoustic pressure satisfies Helmholtz equation:
2p(r)+k 2p(r)=0
According to divergence theorem, can above-mentioned equation be become borderline integral equation:
p ( r ) = ∫ S [ p ( r q ) ∂ G ( r , r q ) ∂ n q - ∂ p ( r q ) ∂ n q G ( r , r q ) ] dS
In the formula, k is the wave number of sound wave, and p (r) is the acoustic pressure of any arbitrarily in the outer sound field of loudspeaker, r qFor on the border a bit, G (r, r q) be Green function, n qUnit normal vector for the border.
Sound pressure level
Figure BSA00000340545900032
In the formula, p 0=2 * 10 -5Pa is a reference acoustic pressure.
Solve the sound pressure level of loudspeaker under ad-hoc location (referring generally to the position along acoustic axis line direction apart from loudspeaker acoustic centre of source 1m) different frequency and can obtain frequency response curve, the sound pressure level that solves loudspeaker space diverse location under a certain CF can obtain the directive property curve.
Used loudspeaker comprises various dynamic loudspeakers.
Described through finding the solution the undamped-free vibration secular equation analysis of loudspeaker finite element model, the multiple mode of oscillation of analysis loudspeaker comprises plug vibration, rocking vibration, flexural vibrations, twisting vibration and compound vibration thereof.
Described three-dimensional drawing software is Pro/E, draws the 3-D geometric model figure of loudspeaker with its, and Pro/E is the incorporate three-dimensional software of a CAD/CAM/CAE, and function comprises that Element Design, mould design, sheet metal component design, virtual assembling and engineering publish picture.
Described grid dividing software is HyperMesh, can divide grid to loudspeaker 3D geometric model with it.
Described finite element solving device is finite element software ANSYS, is used for finding the solution structure, fluid, heat, electromagnetic problem.The boundary element solver is LMS Virtual.Lab, sunykatuib analysis vibration noise characteristic.
Advantage of the present invention is: the present invention designs finite element method and Application of Boundary Element Methods in the loudspeaker analysis; Overcome the shortcoming that equivalent circuit method is confined to low frequency effectively; And can have a fling at the performance of estimating loudspeaker early stage at sample, thus accelerate the loudspeaker design progress, reduce cost of development.
Description of drawings
Fig. 1 is an implementation method process flow diagram of the present invention.
Fig. 2 is the 3D geometric model figure of a 6.5 inches auto loud hailers.
Fig. 3 is the diagrammatic cross-section of Fig. 2.
Fig. 4 is the finite element model of 6.5 inches auto loud hailers.
Fig. 5 is the boundary element model of 6.5 inches auto loud hailers.
Fig. 6 A is one of typical mode of oscillation of 6.5 inches auto loud hailers.
Fig. 6 B be 6.5 inches auto loud hailers typical mode of oscillation two.
Fig. 6 C be 6.5 inches auto loud hailers typical mode of oscillation three.
Fig. 7 is the displacement response curve of 6.5 inches auto loud hailers.
Fig. 8 is the sound pressure frequency response curve of 6.5 inches auto loud hailers.
Fig. 9 is 6.5 inches auto loud hailer directive property utmost point figure (frequency is 1kHz).
Embodiment
Below in conjunction with accompanying drawing and embodiment the present invention is further described.
The present invention is an example with 6.5 inches auto loud hailers, its vibration of sunykatuib analysis and acoustic characteristic.As shown in Figure 1, step of the present invention is following:
Step 1: use three-dimensional drawing software Pro/E to create loudspeaker 3D geometric model figure.
Fig. 2 is the 3D geometric model figure of this loudspeaker 1, and Fig. 3 is the diagrammatic cross-section of this loudspeaker 1, can be known by figure; This loudspeaker 1 is made up of three parts: (1) magnetic circuit system; Comprise T iron 11, magnet steel 10 and train wheel bridge 9, mainly produce magnetic field, can produce driving force energized coils; (2) vibrational system comprises paper awl 6, dust cap 2, ring 3, locating branch slice (bullet ripple) 7, skeleton 8 and voice coil loudspeaker voice coil 8-1, vibration and radiate sound under driving force; (3) backup system comprises basin frame 5 and flanging 4, supports to reach fixedly vibrational system.
Step 2: set up loudspeaker finite element model (see figure 4).Concrete steps are following:
1) the loudspeaker three-dimensional model is imported in the HyperMesh software,, driving force is provided, in the numerical simulation process, be equivalent to the effect of load to speaker vibration system because the magnetic circuit of loudspeaker mainly produces uniform magnetic field.Loudspeaker backup system significant feature is fixing vibrational system, in the numerical simulation process, may instead be fixed boundary condition.Like this, only need in HyperMesh, paper awl 6, ring 3, dust cap 2, voice coil loudspeaker voice coil 8-1 and skeleton 8 to be carried out the body unit that grid dividing obtains each several part.
2) definition unit type.Cell type is selected the Solid45 unit for use.
3) definition material parameter.For this routine loudspeaker 1; It is polypropylene (pp) that 2 of vibrating diaphragm and dust caps are selected material for use, and ring 3 materials are rubber, and locating branch slice (bullet ripple) 7 materials are cotton; Skeleton 8 materials are polyimide (kapton); Voice coil loudspeaker voice coil 8-1 material is a copper, and the material properties of definition above-mentioned material comprises Young modulus, Poisson ratio, density and material damping.
4) definition boundary condition.Because loudspeaker ring 3 is fixed on the basin frame 5 with locating branch slice 7 edges, at the edge of loudspeaker ring 3 and locating branch slice 7 definition fixed boundary condition (along coordinate system x, y, the translation of z axle and rotation displacement are 0).
5) imposed load.On voice coil loudspeaker voice coil 8-1, apply size and be the power of Bil, direction is parallel to voice coil loudspeaker voice coil 8-1 direction of motion, and wherein B is a magnetic density, and i is the electric current of flowing through voice coil 8-1, and l is the length of voice coil loudspeaker voice coil 8-1 lead.
6) preserve the finite element model file.
Step 3: set up the boundary element model (see figure 5).From paper awl 6, ring 3, the body unit surface extraction of dust cap 2 and the veil lattice that body unit is complementary, set up baffle veil lattice again.The veil lattice are selected the shell63 unit for use, preserve the boundary element model file.
Step 4: the finite element model file is imported ANSYS, carry out model analysis, obtain the natural frequency and the mode of oscillation of loudspeaker, several kinds of typical mode of oscillations are seen Fig. 6 A~Fig. 6 C; Carry out humorous response analysis, calculate loudspeaker under a series of frequencies response and obtain displacement, stress, strain curve with change of frequency.Fig. 7 is the displacement curve of analog computation.
Step 5: the displacement result of boundary element model file and step 4 is imported LMS Virtual.Lab, find the solution the frequency response curve and the directive property curve that obtain loudspeaker.Fig. 8 is the frequency response curve of analog computation.Fig. 9 is the directive property utmost point figure of the 1kHz of analog computation.
Above embodiment only in order to the explanation the present invention and and unrestricted technical scheme described in the invention.Therefore, although this instructions has carried out detailed explanation to the present invention with reference to each above-mentioned enforcement,, those of ordinary skill in the art should be appreciated that still and can make amendment or be equal to replacement the present invention; And all do not break away from the technical scheme and the improvement thereof of the spirit and scope of the present invention, and it all should be encompassed in the claim scope of the present invention.

Claims (7)

1. the method for numerical simulation of loudspeaker vibration and acoustic characteristic is characterized in that this method may further comprise the steps at least:
(1) use three-dimensional drawing software to draw the 3-D geometric model figure of loudspeaker;
(2) set up the loudspeaker finite element model, concrete steps are following:
A, loudspeaker 3-D geometric model figure is imported grid dividing software, it is carried out grid dividing, obtain the used body unit of finite element analysis; If need to consider the influence of air, also in model, add air and be body unit its grid dividing to loudspeaker vibration;
B, definition unit type: define the cell type that is complementary for each several part unit in the A step;
C, definition material parameter: the material parameter of each parts of definition model comprises Young modulus, Poisson ratio, density and material damping; If comprise air in the model, also need define the material parameter of air, comprise the velocity of sound, density and acoustical absorption coefficient;
D, definition boundary condition: the boundary condition of loudspeaker finite element model is following: 1) fixed boundary condition: because diaphragm of loudspeaker edge is fixed on the basin frame; Need at its edge definition fixed boundary condition; If loudspeaker also includes locating branch slice, the edge that it and basin bridge joint touch also should define fixed boundary condition; 2) the solid coupled boundary condition of stream: if consider air influence when analyzing the vibration characteristics of loudspeaker, coupling can take place, need to flow coupled boundary condition admittedly in the interphase definition of loudspeaker and air at the surface of contact of loudspeaker and air;
E, give the finite element model imposed load: the voice coil loudspeaker voice coil position at model applies driving force;
F, above-mentioned model is saved as the finite element model file;
(3) set up the loudspeaker boundary element model: the surface of extraction loudspeaker finite element model obtains the veil lattice with finite element model body unit coupling; Because baffle has very big influence to the acoustic characteristic of loudspeaker; Also to set up the veil lattice of baffle; Define the cell type of above-mentioned veil lattice, above-mentioned model is saved as the boundary element model file;
(4) with the finite element solving device finite element model file is found the solution: the finite element model file is imported the finite element solving device and carries out model analysis; Can obtain the natural frequency and the mode of oscillation of loudspeaker; Carry out frequency analysis again, obtain displacement response, strain-responsive and the stress response of loudspeaker;
A, finite element model is carried out model analysis
The fundamental equation of dynamic analysis is following:
[ M ] { u · · } + [ C ] { u · } + [ K ] { u } = { F }
In the following formula; [M] is the mass of system matrix; [C] is the system damping matrix; [K] is the system stiffness matrix;
Figure FSB00000798078100012
is the node vector acceleration,
Figure FSB00000798078100013
be the node speed vector, { u} is the nodal displacement vector; { F} is the load vector that applies, with time correlation;
Suppose that loudspeaker does undamped-free vibration, equation is:
[ M ] { u · · } + [ K ] { u } = { 0 }
For linear system, suppose that system does simple harmonic motion,
{u}={u} i?cosω it
The above-mentioned equation of motion becomes
([K]-ω i 2[M]){u} i={0}
Find the solution above-mentioned equation and can obtain ω i{ u} i, ω iBe the intrinsic circular frequency in i rank, { u} iBe proper vector, i.e. the corresponding mode of oscillation of the intrinsic circular frequency in i rank;
is i rank natural frequency;
B, finite element model is carried out frequency analysis
In frequency analysis, F}, u} is a simple harmonic quantity, has:
{F}=({F 1}+i{F 2})e iΩt
{u}=({u 1}+i{u 2})e iΩt
([K]-Ω 2[M]+iΩ[C])({u 1}+i{u 2})={F 1}+i{F 2}
In the following formula, { u 1, { u 2Be respectively the real part and the imaginary part of displacement, { F 1, { F 2Being respectively the real part and the imaginary part of loading force, Ω is the angular frequency of loading force;
Find the solution above-mentioned equation and can obtain the real part { u of displacement 1And imaginary part { u 2, further can derive stress by displacement;
(5) the finite element solving device is found the solution the displacement and the boundary element model file that obtain and be submitted to the boundary element solver, can find the solution the frequency response curve and the directive property curve that obtain loudspeaker;
In isotropic ideal fluid medium, for the sound wave that simple harmonic quantity changes, acoustic pressure satisfies Helmholtz equation:
▿ 2 p ( r ) + k 2 p ( r ) = 0
According to divergence theorem, can above-mentioned equation be become borderline integral equation:
p ( r ) = ∫ S [ p ( r q ) ∂ G ( r , r q ) ∂ n q - ∂ p ( r q ) ∂ n q G ( r , r q ) ] dS
In the formula, k is the wave number of sound wave, and p (r) is the acoustic pressure of any arbitrarily in the outer sound field of loudspeaker, r qFor on the border a bit, G (r, r q) be Green function, n qUnit normal vector for the border;
Sound pressure level L p = 20 Log ( p p 0 )
In the formula, p 0=2 * 10 -5Pa is a reference acoustic pressure;
Solve the sound pressure level of loudspeaker under the ad-hoc location different frequency and can obtain frequency response curve, the sound pressure level that solves loudspeaker space diverse location under a certain CF promptly obtains the directive property curve.
2. the method for numerical simulation of loudspeaker vibration according to claim 1 and acoustic characteristic is characterized in that loudspeaker comprises various dynamic loudspeakers.
3. the method for numerical simulation of loudspeaker vibration according to claim 1 and acoustic characteristic; It is characterized in that through finding the solution the undamped-free vibration secular equation analysis of loudspeaker finite element model; Analyze the multiple mode of oscillation of loudspeaker, comprise plug vibration, rocking vibration, flexural vibrations, twisting vibration and compound vibration thereof.
4. the method for numerical simulation of loudspeaker vibration according to claim 1 and acoustic characteristic; It is characterized in that three-dimensional drawing software is Pro/E; Draw the 3-D geometric model figure of loudspeaker with it; Pro/E is the incorporate three-dimensional software of a CAD/CAM/CAE, and function comprises that the design of Element Design, mould, sheet metal component design, virtual assembling and engineering publish picture.
5. the method for numerical simulation of loudspeaker vibration according to claim 1 and acoustic characteristic is characterized in that grid dividing software is HyperMesh, can divide grid to loudspeaker 3D geometric model with it.
6. the method for numerical simulation of loudspeaker vibration according to claim 1 and acoustic characteristic is characterized in that the finite element solving device is finite element software ANSYS, is used for finding the solution structure, fluid, heat, electromagnetic problem.
7. the method for numerical simulation of loudspeaker vibration according to claim 1 and acoustic characteristic is characterized in that the boundary element solver is LMS Virtual.Lab, sunykatuib analysis vibration noise characteristic.
CN201010538802A 2010-11-11 2010-11-11 Numerical value simulation method of vibration and acoustic characteristics of speaker Active CN102004823B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010538802A CN102004823B (en) 2010-11-11 2010-11-11 Numerical value simulation method of vibration and acoustic characteristics of speaker

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010538802A CN102004823B (en) 2010-11-11 2010-11-11 Numerical value simulation method of vibration and acoustic characteristics of speaker

Publications (2)

Publication Number Publication Date
CN102004823A CN102004823A (en) 2011-04-06
CN102004823B true CN102004823B (en) 2012-09-26

Family

ID=43812183

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010538802A Active CN102004823B (en) 2010-11-11 2010-11-11 Numerical value simulation method of vibration and acoustic characteristics of speaker

Country Status (1)

Country Link
CN (1) CN102004823B (en)

Families Citing this family (58)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221502B (en) * 2011-05-31 2012-10-31 哈尔滨工业大学 Multi-beam laser heterodyne second harmonic Young modulus measurement method
US9084058B2 (en) 2011-12-29 2015-07-14 Sonos, Inc. Sound field calibration using listener localization
US9219460B2 (en) 2014-03-17 2015-12-22 Sonos, Inc. Audio settings based on environment
US9106192B2 (en) 2012-06-28 2015-08-11 Sonos, Inc. System and method for device playback calibration
US20140019099A1 (en) * 2012-07-16 2014-01-16 Livermore Software Technology Corp Determination Of Failure In Sheet Metal Forming Simulation Using Isotropic Metal Failure Criteria
CN102968543B (en) * 2012-12-13 2015-09-02 嘉善恩益迪电声技术服务有限公司 A kind of Numerical simulation on temperature characteristics method of loudspeaker voice coil and magnetic circuit
CN103294871B (en) * 2013-06-09 2016-07-06 嘉善恩益迪电声技术服务有限公司 A kind of Numerical Simulation Analysis method of loudspeaker voice coil inductance nonlinear characteristic
EP2816332B1 (en) * 2013-06-17 2016-11-30 Siemens Industry Software NV Method and computer product for modeling the sound emission and propagation of systems over a wide frequency range
AT512755B1 (en) * 2013-06-18 2015-03-15 Avl List Gmbh Method for creating a three-dimensional acoustic network
CN103394561B (en) * 2013-08-01 2015-04-15 武汉科技大学 Method for determining bending angle of each pass of cold roll forming of plate and strip
FR3018025B1 (en) * 2014-02-26 2016-03-18 Devialet DEVICE FOR CONTROLLING A SPEAKER
US9264839B2 (en) 2014-03-17 2016-02-16 Sonos, Inc. Playback device configuration based on proximity detection
CN104252559B (en) * 2014-08-29 2018-04-17 浙江中科电声研发中心 A kind of Numerical Simulation Analysis method of loudspeaker multi- scenarios method
US9952825B2 (en) 2014-09-09 2018-04-24 Sonos, Inc. Audio processing algorithms
EP4243450B1 (en) * 2014-09-09 2024-08-07 Sonos, Inc. Method for transmitting a determined audio processing algorithm to a playback device, corresponding playback device, system and computer readable storage medium
CN104408227A (en) * 2014-10-28 2015-03-11 浙江中科电声研发中心 Loudspeaker distortion feature value simulation analysis method
CN105138745B (en) * 2015-08-07 2018-08-10 苏州上声电子股份有限公司 A kind of modification method of vibration component geometrical model in loud speaker simulation analysis
US9693165B2 (en) 2015-09-17 2017-06-27 Sonos, Inc. Validation of audio calibration using multi-dimensional motion check
EP3531714B1 (en) 2015-09-17 2022-02-23 Sonos Inc. Facilitating calibration of an audio playback device
CN105246018B (en) * 2015-10-29 2019-02-26 浙江中科电声研发中心 The approximate calculation method of limited big baffle axial symmetry harmonic distortion of loudspeaker
US9743207B1 (en) 2016-01-18 2017-08-22 Sonos, Inc. Calibration using multiple recording devices
US10003899B2 (en) 2016-01-25 2018-06-19 Sonos, Inc. Calibration with particular locations
US11106423B2 (en) 2016-01-25 2021-08-31 Sonos, Inc. Evaluating calibration of a playback device
US9860662B2 (en) 2016-04-01 2018-01-02 Sonos, Inc. Updating playback device configuration information based on calibration data
US9864574B2 (en) 2016-04-01 2018-01-09 Sonos, Inc. Playback device calibration based on representation spectral characteristics
US9763018B1 (en) 2016-04-12 2017-09-12 Sonos, Inc. Calibration of audio playback devices
CN107315895B (en) * 2016-04-20 2019-12-20 中国人民解放军国防科学技术大学 Semi-analytic method for vibration characteristic analysis of graphene nano-mass sensor
CN107341276A (en) * 2016-04-29 2017-11-10 宝沃汽车(中国)有限公司 A kind of method and system assessed engine performance
CN106055798B (en) * 2016-06-02 2017-09-19 东南大学 An Analysis Method of Acoustic Vibration Response Under Non-stationary Random Dynamic Load
US9794710B1 (en) 2016-07-15 2017-10-17 Sonos, Inc. Spatial audio correction
US10372406B2 (en) 2016-07-22 2019-08-06 Sonos, Inc. Calibration interface
US10459684B2 (en) 2016-08-05 2019-10-29 Sonos, Inc. Calibration of a playback device based on an estimated frequency response
CN106777630A (en) * 2016-12-07 2017-05-31 同济大学 For the determination method of the mode finite element simulation precision of machine components
CN106951628A (en) * 2017-03-16 2017-07-14 吉林航盛电子有限公司 Loudspeaker frequency tracing analysis method and device based on COMSOL softwares
CN107515972B (en) * 2017-08-10 2020-06-23 苏州上声电子股份有限公司 Numerical simulation analysis method for loudspeaker falling process
CN107563063B (en) * 2017-09-05 2021-02-02 无锡吉兴汽车声学部件科技有限公司 Method for establishing acoustic material CAE simulation parameter database
CN108460204B (en) * 2018-02-26 2022-02-18 浙江中科电声研发中心 Method for reversely deducing dynamic mechanical parameters of material of loudspeaker through stress and displacement of vibrating part of loudspeaker
CN108520110A (en) * 2018-03-21 2018-09-11 华南理工大学 Finite Element Analysis Method for Three Types of Thin Film Models of Capacitive Micromachined Ultrasonic Sensors
CN108614921B (en) * 2018-03-30 2022-04-12 北京空间飞行器总体设计部 Low-frequency sound and vibration response prediction method in spacecraft
CN109063343B (en) * 2018-08-08 2023-08-22 苏州上声电子股份有限公司 Finite element method-based loudspeaker screen design method
US10299061B1 (en) 2018-08-28 2019-05-21 Sonos, Inc. Playback device calibration
US11206484B2 (en) 2018-08-28 2021-12-21 Sonos, Inc. Passive speaker authentication
CN109344498B (en) * 2018-09-05 2020-04-21 成都数模码科技有限公司 Design implementation method of wedge module
CN108959835B (en) * 2018-09-28 2023-10-13 苏州上声电子股份有限公司 Strength simulation analysis method of speaker basket during screw installation
CN109145514A (en) * 2018-09-30 2019-01-04 浙江中科电声研发中心 A kind of Numerical Simulation Analysis method of loudspeaker distortions
CN109362020B (en) * 2018-09-30 2020-09-22 浙江中科电声研发中心 Numerical simulation analysis method for dynamic stiffness of loudspeaker frame
CN109684683B (en) * 2018-12-07 2024-03-01 苏州上声电子股份有限公司 A method for analyzing the motion process of loudspeaker vibration system
CN109766613B (en) * 2018-12-29 2020-10-27 西安交通大学 Probability theory-based method for predicting bullet jump of thin-wall structure under thermal noise load
CN110162738A (en) * 2019-03-18 2019-08-23 浙江大学 A kind of calculation method for the radiated assumed based on modality portions coupling
CN110442907B (en) * 2019-07-02 2023-04-28 浙江中科电声研发中心 Numerical simulation analysis method for basic characteristics of piezoelectric MEMS loudspeaker
US10734965B1 (en) 2019-08-12 2020-08-04 Sonos, Inc. Audio calibration of a portable playback device
CN111237348B (en) * 2020-01-14 2021-07-06 上海交通大学 Small waterplane area catamaran type upside-down thrust bearing base and its design method
CN112487679B (en) * 2020-11-20 2025-01-24 四川湖山电器股份有限公司 A numerical simulation method, system, terminal and medium for acoustic characteristics of a speaker
CN113378291B (en) * 2021-05-13 2023-07-14 北京宇航系统工程研究所 Liquid path natural frequency simulation method based on acoustic unit
WO2023056258A1 (en) 2021-09-30 2023-04-06 Sonos, Inc. Conflict management for wake-word detection processes
CN114048649B (en) * 2021-11-09 2024-04-12 华能澜沧江水电股份有限公司 Method for weakening stator low-frequency vibration of hydraulic generator
CN115659759B (en) * 2022-11-11 2023-06-06 西南交通大学 Prediction method of sound transmission loss in high-speed train profile structure based on 2.5-dimensional finite element-boundary element method
CN116132882B (en) * 2022-12-22 2024-03-19 苏州上声电子股份有限公司 Method for determining installation position of loudspeaker

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
毛崎波等.对声辐射模态的讨论.《振动工程学报》.2000,第13卷(第4期),633-637. *
黎胜等.用有限元/边界元方法进行结构声辐射的模态分析.《升学学报》.2001,第26卷(第2期),174-179. *

Also Published As

Publication number Publication date
CN102004823A (en) 2011-04-06

Similar Documents

Publication Publication Date Title
CN102004823B (en) Numerical value simulation method of vibration and acoustic characteristics of speaker
CN106202653B (en) A kind of vehicle body damping distribution optimization method and system
CN104408227A (en) Loudspeaker distortion feature value simulation analysis method
CN109145514A (en) A kind of Numerical Simulation Analysis method of loudspeaker distortions
CN102368281B (en) Numerical value simulating method for loudspeaker magnetic circuit system
CN108959835A (en) Strength Simulation analysis method of the speaker frame in screw installation procedure
CN104504215A (en) Automobile interior acoustic field prediction method based on partition-of-unity finite element-meshless cell
CN104252559A (en) Numerical simulation analysis method for multi-field coupling of loudspeaker
CN110427650A (en) The Numerical Simulation Analysis method of moving iron type loud-speaker fundamental characteristics
Zheng et al. Topology optimization of PCLD on plates for minimizing sound radiation at low frequency resonance
CN105184018A (en) Method for calculating mid-frequency transmission loss of subsystems and optimizing acoustic packages
CN108460204A (en) A method of pushing away its dynamic mechanics parameter of material by the way that the stress of loudspeaker vibration component and displacement are counter
Zhang et al. Three-dimensional coupled vibration theory for the longitudinally polarized piezoelectric ceramic tube
CN113868929B (en) An analytical calculation method for vibration and noise of an outer rotor axial flux hub motor
Kim et al. Modeling tire treadband vibration
CN116011124B (en) A Simulation Modeling Method for Vibration Test System
Fang et al. Design of vibration isolation device for ship power equipment and investigation on its dynamic properties
Pile Méthodes numériques appliquées au calcul des vibrations d’origine électromagnétique: schémas de projection des efforts magnétiques et développement de modeles réduits
Pile Numerical methods applied to the computation of electromagnetic vibrations: magnetic force projection and model reduction
JENSEN et al. Finite element model coupled with lumped parameter elements
WO2020248460A1 (en) Damping analysis method for pedestrian warning loudspeaker
Yao et al. Numerical Study of the Compartment Cavity Problem Using a Novel Cell-Based Smoothed Radial Point Interpolation Method
Beigmoradi et al. Enhancement of Vibration Characteristics of an Air Filter Box Utilizing Numerical Analysis
Shi et al. Dynamic analysis of annular sector plate with general boundary supports
Hansen Investigation of passive control devices for potential application to a launch vehicle structure to reduce the interior noise levels during launch

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: JIASHAN ENYIDI ELECTROACOUSTICAL TECHNOLOGY SERVIC

Effective date: 20121127

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20121127

Address after: Jiashan Jinyang road Jiaxing city Zhejiang province 314100 No. 568 Branch Center No. 1 Building 1 layer

Patentee after: Zhejiang Electro-Acoustic R&D Center, CAS

Patentee after: Jiashan Yiendi Electroacoustic Technology Service Co., Ltd.

Address before: Jiashan Jinyang road Jiaxing city Zhejiang province 314100 No. 568 Branch Center No. 1 Building 1 layer

Patentee before: Zhejiang Electro-Acoustic R&D Center, CAS