WO2024020876A1 - Procédé et appareil de traitement de couplage structurel-acoustique basés sur une réduction d'ordre de modèle, et dispositif terminal - Google Patents

Procédé et appareil de traitement de couplage structurel-acoustique basés sur une réduction d'ordre de modèle, et dispositif terminal Download PDF

Info

Publication number
WO2024020876A1
WO2024020876A1 PCT/CN2022/108321 CN2022108321W WO2024020876A1 WO 2024020876 A1 WO2024020876 A1 WO 2024020876A1 CN 2022108321 W CN2022108321 W CN 2022108321W WO 2024020876 A1 WO2024020876 A1 WO 2024020876A1
Authority
WO
WIPO (PCT)
Prior art keywords
structural
acoustic
boundary element
orthogonal basis
coefficient matrix
Prior art date
Application number
PCT/CN2022/108321
Other languages
English (en)
Chinese (zh)
Inventor
谢祥
Original Assignee
中国科学院深圳先进技术研究院
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 中国科学院深圳先进技术研究院 filed Critical 中国科学院深圳先进技术研究院
Priority to PCT/CN2022/108321 priority Critical patent/WO2024020876A1/fr
Publication of WO2024020876A1 publication Critical patent/WO2024020876A1/fr

Links

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]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/10—Noise analysis or noise optimisation

Definitions

  • This application relates to the field of computer structural acoustics auxiliary technology, and specifically relates to a structural acoustics coupling processing method, device and terminal equipment based on model order reduction.
  • Finite element–boundary element coupling method Finite element–boundary element: FE-BE
  • FE-BE finite element–boundary element
  • Structural vibration in fluid media produces a radiated sound field.
  • the sound pressure of the sound field reacts on the structure, causing additional vibration of the structure. Therefore, when dealing with the coupling problem of structural vibration and fluid sound field, it is necessary to combine the structural and fluid dynamics equations to analyze the structural acoustics. Only by processing the coupled system equations can the sound field sound pressure and structural displacement be obtained.
  • frequency sweep calculation and analysis processing is an indispensable process for performance prediction and optimization design of the structural-acoustic coupling system.
  • most of the existing structural finite element-acoustic boundary element coupling processing methods use iterative methods to solve the problem. Although it reduces memory requirements, its calculation time is long and its efficiency is low. Moreover, due to the frequency-dependent characteristics of boundary elements, the solution is frequently solved.
  • it is necessary to continuously integrate and calculate the coefficient matrix of the system equation for each frequency of interest and solve it. In this way, in the engineering structural system, after adjusting the parameters in the structural acoustic coupling system equation, it is difficult to timely calculate the structural displacement. , sound field sound pressure is controlled. Therefore, it is necessary to study a solution to solve the above problems.
  • One of the purposes of the embodiments of the present application is to provide a structural acoustic coupling processing method, device and terminal equipment based on model reduction, aiming to solve the difficulty of timely analysis of structural displacement, structural acoustic coupling external problems in engineering structural systems.
  • the problem of regulating the sound pressure of the sound field is to provide a structural acoustic coupling processing method, device and terminal equipment based on model reduction, aiming to solve the difficulty of timely analysis of structural displacement, structural acoustic coupling external problems in engineering structural systems.
  • a structural-acoustic coupling processing method based on model order reduction including:
  • the response frequency is analyzed and processed to obtain the target structural displacement and target sound field sound pressure of the engineering structural system.
  • a structural acoustic coupling processing device based on model order reduction including:
  • the first acquisition module is used to acquire the response frequency of the engineering structural system in response to the excitation signal in the preset frequency band;
  • the second acquisition module is used to acquire the structural acoustic coupling system equation of the engineering structural system according to the response frequency;
  • a processing module configured to analyze and process the response frequency based on the reduced-order model corresponding to the structural acoustic coupling system equation, to obtain the target structural displacement and target sound field sound pressure of the engineering structural system.
  • a terminal device in a third aspect, includes a memory and a processor.
  • the memory stores a computer program that can run on the processor.
  • the processor executes the computer program, the computer program is implemented. The steps of the structural acoustic coupling processing method based on model order reduction as described in any one of the above first aspects.
  • the beneficial effect of the structural acoustic coupling processing method based on model reduction is that by constructing a corresponding reduced order model for the structural acoustic coupling system equation, the response frequency in the engineering structural system is analyzed and processed. Accelerate the calculation and improve the efficiency of solving and calculating the target structure displacement and target sound field sound pressure in the engineering structure system. Therefore, in the engineering structure system, the parameters in the structural acoustic coupling system equation are adjusted to timely calculate the target structure displacement and target sound pressure. Control the sound pressure of the sound field; it has strong ease of use and practicality.
  • Figure 1 is a schematic diagram of the application scenario of the structural acoustic coupling processing method based on model order reduction provided by the embodiment of the present application;
  • Figure 2 is a schematic flow chart of the implementation of the structural acoustic coupling processing method based on model order reduction provided by the embodiment of the present application;
  • Figure 3 is a schematic flow chart of the implementation of the structural acoustic coupling processing method based on model order reduction provided by the embodiment of the present application;
  • Figure 4 is a schematic flow chart of the implementation of the structural acoustic coupling processing method based on model order reduction provided by the embodiment of the present application;
  • Figure 5 is a schematic structural diagram of a structural acoustic coupling processing device based on model order reduction provided by an embodiment of the present application
  • Figure 6 is a schematic diagram of a terminal device provided by an embodiment of the present application.
  • first and second are only used for convenience of description and cannot be understood as indicating or implying relative importance or implicitly indicating the quantity of technical features.
  • “Plural” means two or more, unless otherwise expressly and specifically limited.
  • FIG. 1 a schematic diagram of the application scenario of the structural acoustic coupling processing method based on model order reduction provided by the embodiment of the present application.
  • a steady-state structure is immersed in a homogeneous fluid.
  • the vibration of the structure in the homogeneous fluid produces a radiated sound field.
  • the sound pressure of the sound field reacts on the structure, causing additional vibration of the structure and generating structural displacement.
  • the structure and fluid dynamics equations must be combined to solve the problem. Only by solving the structural acoustics coupling system equation can the sound field sound pressure and structural displacement of the engineering structural system be obtained.
  • the target sound field sound pressure and target structure displacement of the engineering structure system can be obtained by performing frequency sweep analysis on the structural acoustic coupling system equation of the engineering structure system.
  • Most of the existing processing methods use iterative methods or direct methods to solve the target sound field sound pressure and target structure displacement of the engineering structural system.
  • the calculation amount is large, the calculation time is long, and the processing efficiency is low. In this way, when the parameters in the engineering structural system are adjusted , it is difficult to timely control the target structure displacement and target sound field sound pressure in the engineering structural system.
  • This application proposes a structural acoustic coupling processing method based on model reduction. It accelerates the analysis and processing process of response frequency through the reduced order model corresponding to the structural acoustic coupling system equation, and improves the target structure displacement and target in engineering structural systems.
  • the efficiency of solving sound field sound pressure Analyze and process the response frequency in the engineering structural system to obtain the target structure displacement and target sound field sound pressure. Images can be constructed separately based on the response frequency and target structure displacement, response frequency and target sound field sound pressure, which can predict the dynamic performance of the engineering structure system. Conduct predictions, conduct topology optimization design of structures, identify parameters of engineering structural systems, and conduct acoustic and vibration sensitivity analysis of engineering structural systems.
  • this structural acoustic coupling processing method based on model order reduction can adjust the parameters of the structural acoustic coupling system equation in the engineering structural system to adjust the target structural displacement and target sound field sound pressure of the engineering structural system, thereby controlling the engineering structure. Vibration of the structure in the system, and noise controlling the sound field.
  • FIG. 2 shows the implementation flow chart of the structural acoustic coupling processing method based on model reduction provided by the embodiment of the present application. The details are as follows:
  • a structural acoustic coupling processing method based on model reduction includes:
  • Step 101 in response to the excitation signal in the preset frequency band, obtain the response frequency of the engineering structure system
  • Step 102 Obtain the structural acoustic coupling system equation of the engineering structural system according to the response frequency
  • Step 103 Based on the reduced-order model corresponding to the structural acoustic coupling system equation, analyze and process the response frequency to obtain the target structural displacement and target sound field sound pressure of the engineering structural system.
  • a corresponding reduced-order model is constructed based on the structural acoustics coupling system equation, and the analysis and processing process of the response frequency is accelerated and calculated, thereby improving the calculation efficiency of solving the engineering structural system to obtain the target structure displacement and target sound field sound pressure, making it possible to perform engineering work
  • the target structure displacement and target sound field sound pressure can be adjusted in a timely manner.
  • the engineering structural system includes a homogeneous fluid and a structure immersed in the homogeneous fluid.
  • the structure generates radiation target sound field sound pressure due to vibration.
  • the target sound field sound pressure reacts on the structure, causing additional vibration of the structure and causing displacement of the target structure.
  • the excitation signal may be an excitation signal of an exciter
  • the exciter may be a mechanical exciter, an electromagnetic exciter, or an electro-hydraulic exciter.
  • the exciter is used to excite the structure in the engineering structural system to generate vibration and generate an excitation signal.
  • the excitation signal is a collection signal of a signal collector.
  • a signal collector is used to collect signals from the vibrating structure to obtain the collected signal, that is, the excitation signal.
  • step 101 specifically, in response to an excitation signal in a preset frequency band, a frequency domain analysis is performed on the excitation signal to obtain the response frequency of the engineering structural system.
  • step 102 obtaining the structural acoustic coupling system equation of the engineering structural system according to the response frequency specifically includes:
  • Step 201 Obtain the first basic parameters of the engineering structural system, where the first basic parameters include structural stiffness, structural damping, structural mass, structural load, acoustic load, homogeneous fluid density, first coupling term, Two coupling terms.
  • Step 202 Establish a geometric model of the engineering structural system based on the first basic parameters.
  • Step 203 Perform triangular meshing processing on the geometric model to obtain second basic parameters of the engineering structural system, where the second basic parameters include source points, field points, interfaces between structures and sound fields, and structures. The normal direction of the interface with the sound field and the average area of the acoustic boundary element triangular mesh.
  • step 203 a triangle meshing process is performed on the geometric model to determine the first geometric model corresponding to the sound field in the engineering structure system. Based on the first geometric model of the sound field, each acoustic boundary element is determined. Based on all acoustic boundary units, the center point of any acoustic boundary unit is determined as the acoustic boundary element configuration point, and the acoustic boundary element configuration point is used as the source point; the center points of the remaining acoustic boundary units are determined as field points.
  • Step 204 Determine third basic parameters of the engineering structure system according to the size of the geometric model, the response frequency and the second basic parameters, where the third basic parameters include a first boundary element coefficient matrix , the second boundary element coefficient matrix.
  • the step 204 specifically includes: determining the truncation order according to the size of the geometric model and the range of the response frequency. According to the truncation order, the response frequency, and the second basic parameter, the third basic parameter of the engineering structure system is determined.
  • a first boundary element kernel function and a second boundary element kernel function are determined, where the first boundary element kernel function is expressed as:
  • the first boundary element kernel function and the second boundary element kernel function are respectively expanded through Taylor's theorem to determine the first boundary element coefficient matrix and the second boundary element coefficient matrix.
  • the exponential coefficient of the first boundary element kernel function and the exponential coefficient of the second boundary element kernel function are respectively expanded through Taylor's theorem, and the expanded coefficients are respectively expanded.
  • the first boundary element kernel function and the expanded second boundary element kernel function are integrated to determine the first boundary element coefficient matrix and the second boundary element coefficient matrix.
  • the first boundary element kernel function and the second boundary element kernel function are respectively expanded through Taylor's theorem to obtain the first boundary element coefficient matrix and the second boundary element coefficient matrix.
  • each frequency of interest needs to be re-integrated and calculated to form a first boundary element coefficient matrix and a second boundary element coefficient matrix.
  • the processing steps are cumbersome and It is computationally intensive.
  • only one integration process is performed for the first boundary element kernel function and the second boundary element kernel function, which facilitates improving the efficiency of subsequent analysis and processing of the response frequency and overcomes the need to use the first boundary element separately in the traditional coupling processing scheme.
  • the boundary element kernel function and the second boundary element kernel function perform multiple integral calculations and have low processing efficiency.
  • Step 205 Determine the structural acoustic coupling system equation according to the first basic parameter, the second basic parameter and the third basic parameter. Among them, the structural acoustic coupling system equation is expressed as follows:
  • K is the structural stiffness
  • D is the structural damping
  • M is the structural mass
  • f s is the structural load
  • f a is the acoustic load
  • ⁇ a is the homogeneous fluid density
  • C sa is the first coupling term
  • C as is the second Coupling term
  • u is the displacement of the target structure
  • p is the sound pressure of the target sound field.
  • the first basic parameters of the engineering structural system are first obtained; and based on the first basic parameters, a geometric model of the engineering structural system is established; the geometric model is triangulated and meshed to obtain the second basic parameters in order to construct structural acoustics
  • the reduced-order model corresponding to the coupled system equation determine the first boundary element coefficient matrix and the second boundary element coefficient matrix based on the geometric model, response frequency, and second basic parameters.
  • the structural-acoustic coupling system equation of the engineering structural system is determined through the first basic parameter, the second basic parameter, and the third basic parameter, so that the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis can be subsequently constructed, and then the reduced order model to improve the processing efficiency of engineering structural systems.
  • the response frequency is analyzed and processed based on the reduced-order model corresponding to formula (1) to obtain the target structural displacement and target sound field sound pressure of the engineering structural system, including:
  • Step 301 According to formula (1), the global orthogonal basis of the structural finite element and the global orthogonal basis of the acoustic boundary element are obtained.
  • Step 302 Determine the second structural displacement and the second sound field sound pressure based on the structural finite element global orthogonal basis, the acoustic boundary element global orthogonal basis and the reduced-order model.
  • the reduced-order model is expressed as follows:
  • K n V H KV;
  • D n V H DV
  • M n V H MV
  • f s,n V H f s ;
  • f a,m W H f a ;
  • C sa,nm V H C sa W;
  • u n is the displacement of the second structure
  • p m is the sound pressure of the second sound field
  • G m is the third boundary element coefficient matrix, G m is expressed as: in,
  • H m is the fourth boundary element coefficient matrix
  • H m is expressed as: in, Where, V is the global orthogonal basis of the structural finite element, W is the global orthogonal basis of the acoustic boundary element, and the superscript represents the conjugate transpose.
  • Step 303 Determine the target structure displacement, the target sound field according to the structural finite element global orthogonal basis, the acoustic boundary element global orthogonal basis, the second structural displacement and the second sound field sound pressure. Sound pressure.
  • the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained.
  • the reduced-order second structural displacement and second sound field sound pressure are determined; this avoids direct processing of the original structural-acoustic coupling system equations , determine the reduced-order second structural displacement and second sound field sound pressure through the reduced-order model, thereby improving the processing efficiency of the engineering structural system.
  • the target structure displacement and target sound field sound pressure of the engineering structural system are obtained.
  • This structural acoustics coupling processing method based on model order reduction can adjust the parameters of the engineering structure system to timely adjust the target structure displacement and target sound field sound pressure of the engineering structure system, thereby controlling the vibration of the structure in the engineering structure system, and Control the noise in the sound field to achieve subsequent structural optimization design of the engineering structural system under the constraints of structural load and acoustic load.
  • the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained.
  • This application provides four examples of obtaining the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis under different constraint conditions of structural load and acoustic load:
  • the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained, specifically including:
  • the first structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis
  • the first acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.
  • the first structural finite element global orthogonal basis and the first acoustic boundary element global orthogonal basis are obtained, specifically including:
  • Step 401 Ignore the first coupling term of formula (1), and calculate the first structural finite element global orthogonal basis through the second-order Arnoldi algorithm.
  • step 401 specifically, the first coupling term of formula (1) is ignored to obtain the first structural acoustic coupling system equation, where the first structural acoustic coupling system equation is expressed as:
  • K is the structural stiffness
  • D is the structural damping
  • M is the structural mass
  • f s is the structural load
  • u is the target structure displacement
  • ⁇ is the response frequency
  • the first structural finite element global orthogonal basis is calculated through the second-order Arnoldi algorithm, which specifically includes:
  • the first matrix, the second matrix, and the first starting vector are obtained, where the first matrix is expressed as:
  • the second matrix is expressed as:
  • Step 402 Based on the second basic parameters, truncate the first boundary element coefficient matrix and the second boundary element coefficient matrix of formula (1) to obtain a first sparse coefficient matrix and a second sparse coefficient matrix. . Specifically, based on the second basic parameters and according to the characteristics of the first boundary element kernel function and the second boundary element kernel function that decay with distance, the first boundary element coefficient matrix of formula (1), The second boundary element coefficient matrix is truncated to obtain a first sparse coefficient matrix and a second sparse coefficient matrix.
  • the cutoff radius is determined based on the average area of the acoustic boundary element triangular mesh, where the cutoff radius can be expressed as mean(A i ) is the average area of the acoustic boundary element triangular mesh.
  • a triangle meshing process is performed on the geometric model to determine the first geometric model corresponding to the sound field in the engineering structure system. Based on the first geometric model of the sound field, several triangular units are obtained. For all triangular units, the triangular units are used as the aforementioned acoustic boundary units. For all acoustic boundary units, the center of the acoustic boundary unit is used as the boundary element configuration point.
  • one of the acoustic boundary units is used as the source unit, the boundary element configuration point of the acoustic boundary unit is used as the source point, and the source point is used as the center of the ball; the boundary element configuration points of the remaining acoustic boundary units are used as field points. .
  • a truncated sphere geometric model is established; based on the truncated sphere geometric model, the acoustic boundary unit located in the truncated sphere is determined to be a strong interaction unit; and a k-d tree search algorithm is used to determine adjacent The number of strong interaction units; determining the first sparse coefficient matrix and the second sparse coefficient matrix based on the first boundary element coefficient matrix, the second boundary element coefficient matrix, and the number of adjacent strong interaction units .
  • Step 403 Based on the first sparse coefficient matrix and the second sparse coefficient matrix, calculate the first acoustic boundary element global orthogonal basis through the Arnoldi algorithm or the intrinsic orthogonal decomposition algorithm.
  • a second structural acoustic coupling system equation is obtained, wherein the second structural acoustic coupling system equation is expressed as:
  • ⁇ a is the homogeneous fluid density
  • C as is the second coupling term
  • ⁇ is the response frequency
  • p is the target sound field sound pressure
  • p is the target sound field sound pressure
  • ⁇ indicates that it is a sparse approximation of the original dense matrix
  • the first acoustic boundary element global orthogonal basis is calculated through the Arnoldi algorithm.
  • the third matrix and the second starting vector are obtained, where the third matrix is expressed as:
  • the second starting vector is expressed as:
  • a set of first acoustic boundary element local orthogonal bases is obtained; based on the set of first acoustic boundary element local orthogonal bases, the first set of acoustic boundary element local orthogonal bases are generated for different expansion points.
  • Acoustic boundary element local orthogonal basis perform re-orthogonalization processing on the first acoustic boundary element local orthogonal basis to obtain the first acoustic boundary element global orthogonal basis.
  • the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained, specifically including:
  • the second structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis
  • the second acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.
  • the second structural finite element global orthogonal basis and the second acoustic boundary element global orthogonal basis are obtained, specifically including:
  • Step 501 Based on the second basic parameters, truncate the first boundary element coefficient matrix and the second boundary element coefficient matrix of formula (1) to obtain a third sparse coefficient matrix and a fourth sparse coefficient matrix. .
  • Step 502 Based on the third sparse coefficient matrix and the fourth sparse coefficient matrix, ignore the second coupling term of formula (1), and calculate the second coupling term through the Arnoldi algorithm or the intrinsic orthogonal decomposition algorithm. Acoustic boundary element global orthonormal basis.
  • Equation where the third structure acoustic coupling system equation is expressed as:
  • p is the target sound field sound pressure
  • f a is the acoustic load
  • ⁇ indicates that it is a sparse approximation of the original dense matrix
  • the second acoustic boundary element global orthogonal basis is calculated through the Arnoldi algorithm.
  • the fourth matrix and the third starting vector are obtained, where the fourth matrix is expressed as:
  • the third starting vector is expressed as:
  • a set of second acoustic boundary element local orthogonal bases is obtained; based on the set of second acoustic boundary element local orthogonal bases, the first set of acoustic boundary element local orthogonal bases are generated for different expansion points.
  • Two acoustic boundary element local orthogonal bases perform re-orthogonalization processing on the second acoustic boundary element local orthogonal basis to obtain the second acoustic boundary element global orthogonal basis.
  • Step 503 Calculate the second-order finite element global orthogonal basis of the second structure through the second-order Arnoldi algorithm.
  • step 503 specifically, according to the third starting vector, a fourth structural-acoustic coupling system equation is obtained, wherein the fourth structural-acoustic coupling system equation is expressed as:
  • K is the structural stiffness
  • D is the structural damping
  • M is the structural mass
  • u is the target structure displacement
  • ⁇ is the response frequency
  • C sa is the first coupling term.
  • the second structural finite element global orthogonal basis is calculated through the second-order Arnoldi algorithm.
  • the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained, specifically including:
  • the third structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis
  • the third acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.
  • the third structural finite element global orthogonal basis and the third acoustic boundary element global orthogonal basis are obtained, specifically including:
  • Step 601 Ignore the first coupling term of formula (1), and calculate the third structural finite element global orthogonal basis through the second-order Arnoldi algorithm.
  • step 601 specifically, the first coupling term of formula (1) is ignored to obtain the fifth structural acoustic coupling system equation, where the fifth structural acoustic coupling system equation is expressed as:
  • K is the structural stiffness
  • D is the structural damping
  • M is the structural mass
  • f s is the structural load
  • u is the target structure displacement
  • ⁇ is the response frequency
  • j is the imaginary unit
  • the third structural finite element global orthogonal basis is calculated through the second-order Arnoldi algorithm.
  • the fifth matrix, the sixth matrix, and the fourth starting vector are obtained, where the fifth matrix is expressed as:
  • the sixth matrix is expressed as:
  • a set of third structural finite element local orthogonal bases is obtained; based on the set of third structural finite element local orthogonal bases, for different The expansion point generates the third structural finite element local orthogonal basis; the third structural finite element local orthogonal basis is re-orthogonalized to obtain the third structural finite element global orthogonal basis.
  • Step 602 Based on the second basic parameters, truncate the first boundary element coefficient matrix and the second boundary element coefficient matrix of formula (1) to obtain a fifth sparse coefficient matrix and a sixth sparse coefficient matrix. .
  • Step 603 Based on the fifth sparse coefficient matrix and the sixth sparse coefficient matrix, calculate the third acoustic boundary element global orthogonal basis through the Arnoldi algorithm or the intrinsic orthogonal decomposition algorithm.
  • steps 602 to 603 specifically, based on the second basic parameter, the first boundary element coefficient matrix and the second boundary element coefficient matrix of formula (1) are Perform truncation processing to obtain the fifth sparse coefficient matrix and the sixth sparse coefficient matrix. Based on the fifth sparse coefficient matrix and the sixth sparse coefficient matrix, a sixth structural acoustic coupling system equation is obtained, where the sixth structural acoustic coupling system equation is expressed as:
  • ⁇ a is the homogeneous fluid density
  • ⁇ is the response frequency
  • p is the target sound field sound pressure
  • j is the imaginary unit
  • f a is the sixth sparse coefficient matrix
  • f a is the acoustic load
  • the superscript ⁇ indicates that it is a sparse approximation of the original dense matrix.
  • the third acoustic boundary element global orthogonal basis is calculated through the Arnoldi algorithm.
  • the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained, specifically including:
  • the fourth structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis
  • the fourth acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.
  • the fourth structural finite element global orthogonal basis and the fourth acoustic boundary element global orthogonal basis are obtained, specifically including:
  • Step 701 Based on the second basic parameter, truncate the first boundary element coefficient matrix and the second boundary element coefficient matrix of formula (1) to obtain a seventh sparse coefficient matrix and an eighth sparse coefficient matrix. ;
  • Step 702 Based on the seventh sparse coefficient matrix and the eighth sparse coefficient matrix, the second coupling term of formula (1) is ignored, and the fourth coupling term is calculated through the Arnoldi algorithm or the intrinsic orthogonal decomposition algorithm. Acoustic boundary element global orthonormal basis.
  • steps 701 to 702 specifically, based on the second basic parameter, the first boundary element coefficient matrix and the second boundary element coefficient matrix of formula (1) are Perform truncation processing to obtain the seventh sparse coefficient matrix and the eighth sparse coefficient matrix;
  • the second coupling term of formula (1) is ignored to obtain a seventh structural acoustic coupling system equation, wherein the seventh structural acoustic coupling system equation
  • the equation is expressed as:
  • p is the target sound field sound pressure
  • f a is the acoustic load
  • ⁇ indicates that it is a sparse approximation of the original dense matrix
  • the fourth acoustic boundary element global orthogonal basis is calculated through the Arnoldi algorithm.
  • the seventh matrix and the fifth starting vector are obtained, where the seventh matrix is expressed as:
  • the fifth starting vector is expressed as:
  • a set of fourth acoustic boundary element local orthogonal basis is obtained; based on the set of fourth acoustic boundary element local orthogonal basis, the fourth acoustic boundary element is generated for different expansion points Boundary element local orthogonal basis; performing re-orthogonalization processing on the fourth acoustic boundary element local orthogonal basis to obtain the fourth acoustic boundary element global orthogonal basis.
  • Step 703 Calculate the fourth structural finite element global orthogonal basis through the second-order Arnoldi algorithm.
  • step 703 specifically, based on the fifth starting vector, an eighth structural acoustic coupling system equation is obtained, wherein the eighth structural acoustic coupling system equation is expressed as:
  • K is the structural stiffness
  • D is the structural damping
  • M is the structural mass
  • u is the target structure displacement
  • ⁇ is the response frequency
  • j is the imaginary unit
  • C sa is the first coupling term.
  • the fourth structural finite element global orthogonal basis is calculated through the second-order Arnoldi algorithm.
  • Example 1 shows specific load conditions where only structural loads exist
  • Embodiment 2 shows that only acoustic loads exist
  • Specific load conditions Example 3 and Example 4 show the specific load conditions that exist for both structural load and acoustic load.
  • the structural acoustic coupling processing method based on model order reduction can meet the processing needs under different load conditions, by truncating the first boundary element coefficient matrix and the second boundary element coefficient matrix, and by selecting the structure
  • the first coupling term or the second coupling term of the acoustic coupling system equation is ignored to reduce the amount of data that needs to be processed for the structural acoustic coupling system equation and improve the construction of structural finite element global orthogonal basis and acoustic boundary element global orthogonal basis. s efficiency.
  • sequence number of each step in the above embodiment does not mean the order of execution.
  • the execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiment of the present application.
  • Figure 5 shows a structural block diagram of the structural acoustic coupling processing device based on model reduction provided by the embodiment of the present application. For convenience of explanation, only the parts related to the embodiment of the present application are shown. .
  • the structure-acoustic coupling processing device based on model order reduction illustrated in Figure 5 may be the execution subject of the structure-acoustic coupling processing method based on model order reduction provided in the foregoing embodiments.
  • the structural acoustic coupling processing device 50 based on model order reduction includes:
  • the first acquisition module 51 is used to acquire the response frequency of the engineering structural system in response to the excitation signal in the preset frequency band;
  • the second acquisition module 52 is used to acquire the structural acoustic coupling system equation of the engineering structural system according to the response frequency;
  • the processing module 53 is configured to analyze and process the response frequency based on the reduced-order model corresponding to the structural acoustic coupling system equation, and obtain the target structural displacement and target sound field sound pressure of the engineering structural system.
  • Figure 6 is a schematic structural diagram of a terminal device provided by an embodiment of the present application.
  • the terminal device 60 of this embodiment includes: at least one processor 61 (only one is shown in Figure 6), a memory 62, and the memory 62 stores programs that can run on the processor 61.
  • Computer program 63 When the processor 61 executes the computer program 63 , the steps in each of the above embodiments of the structural acoustic coupling processing method based on model order reduction are implemented, such as steps 101 to 103 shown in FIG. 2 . Alternatively, when the processor 61 executes the computer program 63, it implements the functions of each module/unit in each of the above device embodiments, such as the functions of modules 51 to 53 shown in Figure 5 .
  • the terminal device 60 may be a computing device such as a desktop computer, a notebook, a PDA, a cloud server, etc.
  • the terminal device may include, but is not limited to, a processor 61 and a memory 62 .
  • FIG. 6 is only an example of the terminal device 60 and does not constitute a limitation on the terminal device 60. It may include more or fewer components than shown in the figure, or some components may be combined, or different components may be used. , for example, the terminal device may also include an input sending device, a network access device, a bus, etc.
  • the so-called processor 61 can be a central processing unit (Central Processing Unit, CPU), or other general-purpose processor, digital signal processor (Digital Signal Processor, DSP), application specific integrated circuit (Application Specific Integrated Circuit, ASIC), Ready-made field-programmable gate array (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc.
  • a general-purpose processor may be a microprocessor or the processor may be any conventional processor, etc.
  • the memory 62 may be an internal storage unit of the terminal device 60 , such as a hard disk or memory of the terminal device 60 .
  • the memory 62 may also be an external storage device of the terminal device 60, such as a plug-in hard disk, a smart memory card (Smart Media Card, SMC), or a secure digital (SD) equipped on the terminal device 60. Card, Flash Card, etc.
  • the memory 62 may also include both an internal storage unit of the terminal device 60 and an external storage device.
  • the memory 62 is used to store operating systems, application programs, boot loaders, data and other programs, such as program codes of the computer programs.
  • the memory 62 may also be used to temporarily store data that has been sent or is to be sent.
  • each functional unit in each embodiment of the present application can be integrated into one processing unit, each unit can exist physically alone, or two or more units can be integrated into one unit.
  • the above integrated units can be implemented in the form of hardware or software functional units.
  • An embodiment of the present application also provides a terminal device, which includes at least one memory, at least one processor, and a computer program stored in the at least one memory and executable on the at least one processor, so When the processor executes the computer program, the terminal device implements the steps in any of the above method embodiments.
  • Embodiments of the present application also provide a computer-readable storage medium.
  • the computer-readable storage medium stores a computer program.
  • the steps in each of the above method embodiments can be implemented.

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)
  • Soundproofing, Sound Blocking, And Sound Damping (AREA)

Abstract

La présente invention concerne un procédé et un appareil de traitement de couplage structurel-acoustique basés sur une réduction d'ordre de modèle, et un dispositif terminal, qui sont applicables au domaine technique de l'assistance acoustique structurelle informatique. Le procédé comprend : en réponse à un signal d'excitation d'une bande de fréquence prédéfinie, l'acquisition de la fréquence de réponse d'un système de structure d'ingénierie ; l'acquisition d'une équation de système de couplage structurel-acoustique du système de structure d'ingénierie en fonction de la fréquence de réponse ; et la conduite d'un traitement d'analyse sur la fréquence de réponse sur la base d'un modèle de réduction d'ordre correspondant à l'équation de système de couplage structurel-acoustique, de façon à obtenir un déplacement structurel cible et la pression sonore d'un champ sonore cible du système de structure d'ingénierie. Dans le procédé de traitement de couplage structural-acoustique basé sur une réduction d'ordre de modèle, des paramètres d'un système de structure d'ingénierie peuvent être ajustés, de sorte qu'un déplacement structurel cible et la pression sonore d'un champ sonore cible du système de structure d'ingénierie sont régulés et contrôlés en temps opportun, de façon à commander la vibration de la structure dans le système de structure d'ingénierie et contrôler le bruit du champ sonore. Par conséquent, le procédé de traitement de couplage structurel-acoustique présente une aptitude à l'emploi et une convivialité relativement élevées.
PCT/CN2022/108321 2022-07-27 2022-07-27 Procédé et appareil de traitement de couplage structurel-acoustique basés sur une réduction d'ordre de modèle, et dispositif terminal WO2024020876A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/CN2022/108321 WO2024020876A1 (fr) 2022-07-27 2022-07-27 Procédé et appareil de traitement de couplage structurel-acoustique basés sur une réduction d'ordre de modèle, et dispositif terminal

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2022/108321 WO2024020876A1 (fr) 2022-07-27 2022-07-27 Procédé et appareil de traitement de couplage structurel-acoustique basés sur une réduction d'ordre de modèle, et dispositif terminal

Publications (1)

Publication Number Publication Date
WO2024020876A1 true WO2024020876A1 (fr) 2024-02-01

Family

ID=89704822

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/108321 WO2024020876A1 (fr) 2022-07-27 2022-07-27 Procédé et appareil de traitement de couplage structurel-acoustique basés sur une réduction d'ordre de modèle, et dispositif terminal

Country Status (1)

Country Link
WO (1) WO2024020876A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117744455A (zh) * 2024-02-21 2024-03-22 大连理工大学 考虑加载偏差的试验数字孪生体动态校正方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102270252A (zh) * 2010-06-04 2011-12-07 利弗莫尔软件技术公司 对结构进行振动-声学分析的系统和方法
CN103559366A (zh) * 2013-11-18 2014-02-05 中国汽车工程研究院股份有限公司 基于光滑有限元边界元法的车身结构声场耦合预测方法
US20150127311A1 (en) * 2013-11-06 2015-05-07 Weidlinger Associates, Inc. Computer Implemented Apparatus and Method for Finite Element Modeling Using Hybrid Absorbing Element
CN104850721A (zh) * 2015-06-03 2015-08-19 湖南大学 一种基于概率和区间的混合进行外声场预测的方法和装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102270252A (zh) * 2010-06-04 2011-12-07 利弗莫尔软件技术公司 对结构进行振动-声学分析的系统和方法
US20150127311A1 (en) * 2013-11-06 2015-05-07 Weidlinger Associates, Inc. Computer Implemented Apparatus and Method for Finite Element Modeling Using Hybrid Absorbing Element
CN103559366A (zh) * 2013-11-18 2014-02-05 中国汽车工程研究院股份有限公司 基于光滑有限元边界元法的车身结构声场耦合预测方法
CN104850721A (zh) * 2015-06-03 2015-08-19 湖南大学 一种基于概率和区间的混合进行外声场预测的方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
XIE XIANG; LIU YIJUN: "Efficient multi-frequency solutions of FE–BE coupled structural–acoustic problems using Arnoldi-based dimension reduction approach", COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, NORTH-HOLLAND, AMSTERDAM, NL, vol. 386, 5 September 2021 (2021-09-05), AMSTERDAM, NL , XP086855056, ISSN: 0045-7825, DOI: 10.1016/j.cma.2021.114126 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117744455A (zh) * 2024-02-21 2024-03-22 大连理工大学 考虑加载偏差的试验数字孪生体动态校正方法

Similar Documents

Publication Publication Date Title
Abedi et al. An h-adaptive spacetime-discontinuous Galerkin method for linear elastodynamics
Mohammadi et al. A Galerkin-reproducing kernel method: Application to the 2D nonlinear coupled Burgers' equations
WO2024020876A1 (fr) Procédé et appareil de traitement de couplage structurel-acoustique basés sur une réduction d'ordre de modèle, et dispositif terminal
CN105469446A (zh) 点云网格简化系统及方法
Liew et al. Mesh‐free radial point interpolation method for the buckling analysis of Mindlin plates subjected to in‐plane point loads
Bodony et al. Provably stable overset grid methods for computational aeroacoustics
Titarev et al. Numerical modeling of the influence of the relative positions of a propeller and pylon on turboprop aircraft noise
Hu An efficient solution of time domain boundary integral equations for acoustic scattering and its acceleration by Graphics Processing Units
White et al. A dual mesh method with adaptivity for stress-constrained topology optimization
Nguyen et al. A novel hr-adaptive mesh refinement scheme for stress-constrained shape and topology optimization using level-set-based trimmed meshes
Chen et al. Recent advances on singular boundary method
TW202138999A (zh) 用於卷積運算的資料劃分方法及處理器
Cai et al. A high order semi-Lagrangian discontinuous Galerkin method for the two-dimensional incompressible Euler equations and the guiding center Vlasov model without operator splitting
CN111707262A (zh) 基于最近点向量投影的点云匹配方法、介质、终端和装置
Davis A simplified approach for predicting interaction between flexible structures and acoustic enclosures
CN114757347A (zh) 低比特量化神经网络加速器实现方法和系统
Zhang et al. An improved multiquadric collocation method for 3-D electromagnetic problems
CN111783351A (zh) 针对结构系统不确定性参数的非概率可信集合定量化方法
CN115345044A (zh) 基于模型降阶的结构声学耦合处理方法、装置及终端设备
Weinhardt et al. GFEM stabilization techniques applied to dynamic analysis of non-uniform section bars
CN112214842B (zh) 声衬设计方法、装置、设备及存储介质
Zhang et al. Fast realization of 3-D space-time correlation sea clutter of large-scale sea scene based on FPGA: From EM model to statistical model
Rumpler et al. An adaptive strategy for the bivariate solution of finite element problems using multivariate nested Padé approximants
Wei et al. RBF collocation method and stability analysis for phononic crystals
Chen et al. Implementation of parallel medical ultrasound imaging algorithm on CAPI-enabled FPGA

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22952339

Country of ref document: EP

Kind code of ref document: A1