Structural Acoustic Coupling Processing Method, Device and Terminal Equipment Based on Model Order Reduction

Abstract

The present application provides a structural acoustic coupling processing method, apparatus, and terminal device based on model order reduction, applicable to the field of computer-aided structural acoustics technology. The method comprises: obtaining a response frequency of an engineering structure system in response to an excitation signal in a preset frequency band; obtaining a structural acoustic coupling system equation for the engineering structure system based on the response frequency; and analyzing and processing the response frequency based on a reduced-order model corresponding to the structural acoustic coupling system equation to obtain a target structural displacement and target sound field pressure for the engineering structure system. This structural acoustic coupling processing method based on model order reduction can timely adjust the target structural displacement and target sound field pressure of the engineering structure system by adjusting the parameters of the engineering structure system, thereby controlling the vibration of the structure in the engineering structure system and controlling the noise of the sound field. It is highly user-friendly and practical.

Description

Technical Field

[0001] The present application belongs to the field of computer-aided structural acoustics technology, and in particular relates to a structural acoustic coupling processing method, apparatus, and terminal device based on model order reduction. Background Art

[0002] With increasingly stringent requirements for noise emissions, the dynamics and acoustic qualities of a system have become important performance evaluation indicators in fields such as military equipment and aerospace. For problems other than structural-acoustic coupling, the finite element–boundary element (FE-BE) coupling method is widely used to evaluate the vibroacoustic characteristics of a system, leveraging their respective advantages. Structural vibrations in a fluid medium generate a radiated acoustic field. At the same time, the acoustic pressure of the acoustic field reacts on the structure, causing additional vibrations. Therefore, when addressing the coupling problem between structural vibration and the fluid acoustic field, it is necessary to combine the structural and fluid dynamics equations and process the structural-acoustic coupling system equations to obtain the acoustic pressure and structural displacement.

[0003] When processing the equations for coupled structural-acoustic systems, frequency sweep calculation and analysis are essential for performance estimation and optimal design. However, existing methods for processing coupled structural finite element (FEM) and acoustic boundary element (BEM) systems often employ iterative methods. While these methods reduce memory requirements, they suffer from long computation times and low efficiency. Furthermore, due to the frequency-dependent nature of BEMs, the frequency response function (FRF) requires continuous integration and calculation of the system's coefficient matrix for each frequency of interest. Consequently, in engineering structural systems, adjusting the parameters in the coupled structural-acoustic system makes it difficult to timely control structural displacement and acoustic pressure. Therefore, it is necessary to develop a solution to these problems. Summary of the Invention

[0004] In view of this, the embodiments of the present application provide a structural acoustic coupling processing method, apparatus and terminal device based on model order reduction, which can solve the problem of difficulty in timely regulating structural displacement and sound field sound pressure when dealing with external problems of structural acoustic coupling in engineering structural systems.

[0005] A first aspect of an embodiment of the present application provides a structural acoustic coupling processing method based on model order reduction, comprising:

[0006] In response to an excitation signal in a preset frequency band, a response frequency of the engineering structure system is obtained;

[0007] Obtaining a structural acoustic coupling system equation of the engineering structure system according to the response frequency;

[0008] Based on the reduced-order model corresponding to the structural-acoustic coupling system equation, the response frequency is analyzed and processed to obtain the target structural displacement and target sound field pressure of the engineering structure system.

[0009] In a possible implementation of the first aspect, obtaining a structural-acoustic coupling system equation of the engineering structure system according to the response frequency specifically includes:

[0010] Acquiring first basic parameters of the engineering structure system, wherein the first basic parameters include structural stiffness, structural damping, structural mass, structural load, acoustic load, homogeneous fluid density, a first coupling term, and a second coupling term;

[0011] Establishing a geometric model of the engineering structure system based on the first basic parameters;

[0012] Performing triangular meshing on the geometric model to obtain second basic parameters of the engineering structure system, wherein the second basic parameters include a source point, a field point, an interface between the structure and the acoustic field, a normal direction of the interface between the structure and the acoustic field, and an average area of an acoustic boundary element triangular mesh;

[0013] Determining a truncation order according to the size of the geometric model and the range of the response frequency;

[0014] Determining a first boundary element kernel function and a second boundary element kernel function based on the second basic parameter and the response frequency;

[0015] Based on the truncation order, the first boundary element kernel function and the second boundary element kernel function are respectively expanded using Taylor's theorem to determine third basic parameters of the engineering structure system, wherein the third basic parameters include a first boundary element coefficient matrix and a second boundary element coefficient matrix;

[0016] determining the structural acoustic coupling system equation according to the first basic parameter, the second basic parameter, and the third basic parameter;

[0017] Wherein, the first boundary element kernel function is expressed as: The second boundary element kernel function is expressed as: Wherein, k is the wave number, k is expressed as k=ω / c, ω is the response frequency, c is the speed of sound in air, and j is an imaginary unit; x is the source point, y is the field point, r is the distance between the source point x and the field point y, and r is expressed as r=|xy|; n(y) is the normal direction of the interface between the structure and the sound field;

[0018] Among them, the first boundary element coefficient matrix is expressed as The second boundary element coefficient matrix is expressed as Where, ξ=jk, L is the truncation order, Γ a is the boundary surface between the structure and the sound field.

[0019] In a possible implementation of the first aspect, analyzing and processing the response frequency based on a reduced-order model corresponding to the structural-acoustic coupling system equation to obtain a target structural displacement and a target sound field pressure of the engineering structure system includes:

[0020] According to the structural-acoustic coupling system equation, a global orthogonal basis of a structural finite element and a global orthogonal basis of an acoustic boundary element are obtained;

[0021] Determining a second structural displacement and a 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;

[0022] Determining the target structural displacement and the target sound field pressure 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;

[0023] The reduced-order model is expressed as follows:

[0024]

[0025] Among them, K n is the second structural stiffness, K n Expressed as: K n =V H KV;

[0026] D n is the second structural damping, D n Expressed as: D n =V H DV;

[0027] M n is the second structural mass, M n Expressed as: M n =V H MV;

[0028] f s,n is the second structural load, f s,n Expressed as: f s,n =V H f s ;

[0029] f a,m is the second acoustic load, f a,m Expressed as: f a,m =W H f a ;

[0030] C sa,nm is the third coupling term, C sa,nm Expressed as: C sa,nm =V H C sa W;

[0031] C as,n is the fourth coupling term, C as,n Expressed as: C as,n =C as V;

[0032] G m is the third boundary element coefficient matrix, G m Expressed as: in,

[0033] H m is the fourth boundary element coefficient matrix, H m Expressed as: in, where the superscript H denotes the conjugate transpose; and

[0034] ρ a is the density of the homogeneous fluid; u n is the second structural displacement; p m is the sound pressure of the second sound field;

[0035] The target structure displacement is expressed as: u = Vu n , the target sound field sound pressure is expressed as: p = Wp m ;

[0036] Wherein, V is the global orthogonal basis of the structural finite element, and W is the global orthogonal basis of the acoustic boundary element.

[0037] In a possible implementation of the first aspect, obtaining the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis according to the structural-acoustic coupling system equation specifically includes:

[0038] When the structural load is not zero and the acoustic load is zero, the first coupling term of the structural-acoustic coupling system equation is ignored, and a first structural finite element global orthogonal basis is obtained by calculation using a second-order Arnoldi algorithm;

[0039] Based on the second basic parameters, truncate the first boundary element coefficient matrix and the second boundary element coefficient matrix of the structural acoustic coupling system equation to obtain a first sparse coefficient matrix and a second sparse coefficient matrix;

[0040] Based on the first sparse coefficient matrix and the second sparse coefficient matrix, a first acoustic boundary element global orthogonal basis is calculated by an Arnoldi algorithm or an eigenorthogonal decomposition algorithm;

[0041] The first structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the first acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

[0042] In a possible implementation of the first aspect, obtaining the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis according to the structural-acoustic coupling system equation specifically includes:

[0043] When the structural load is zero and the acoustic load is not zero, truncating the first boundary element coefficient matrix and the second boundary element coefficient matrix of the structural acoustic coupling system equation based on the second basic parameter to obtain a third sparse coefficient matrix and a fourth sparse coefficient matrix;

[0044] Based on the third sparse coefficient matrix and the fourth sparse coefficient matrix, the second coupling term of the structural acoustic coupling system equation is ignored, and a second acoustic boundary element global orthogonal basis is calculated by an Arnoldi algorithm or an eigenorthogonal decomposition algorithm;

[0045] The second structure finite element global orthogonal basis is obtained by calculating the second-order Arnoldi algorithm;

[0046] The second structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the second acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

[0047] In a possible implementation of the first aspect, obtaining the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis according to the structural-acoustic coupling system equation specifically includes:

[0048] When both the structural load and the acoustic load are non-zero, the first coupling term of the structural-acoustic coupling system equation is ignored, and a third structural finite element global orthogonal basis is obtained by calculation using a second-order Arnoldi algorithm;

[0049] Based on the second basic parameters, truncate the first boundary element coefficient matrix and the second boundary element coefficient matrix of the structural acoustic coupling system equation to obtain a fifth sparse coefficient matrix and a sixth sparse coefficient matrix;

[0050] Based on the fifth sparse coefficient matrix and the sixth sparse coefficient matrix, a third acoustic boundary element global orthogonal basis is calculated by an Arnoldi algorithm or an eigenorthogonal decomposition algorithm;

[0051] The third structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the third acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

[0052] In a possible implementation of the first aspect, obtaining the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis according to the structural-acoustic coupling system equation specifically includes:

[0053] When both the structural load and the acoustic load are non-zero, truncating the first boundary element coefficient matrix and the second boundary element coefficient matrix of the structural-acoustic coupling system equation based on the second basic parameter to obtain a seventh sparse coefficient matrix and an eighth sparse coefficient matrix;

[0054] Based on the seventh sparse coefficient matrix and the eighth sparse coefficient matrix, the second coupling term of the structural acoustic coupling system equation is ignored, and a fourth acoustic boundary element global orthogonal basis is calculated by an Arnoldi algorithm or an eigenorthogonal decomposition algorithm;

[0055] The fourth structure finite element global orthogonal basis is obtained by calculating the second-order Arnoldi algorithm;

[0056] The fourth structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the fourth acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

[0057] A second aspect of an embodiment of the present application provides a structural acoustic coupling processing device based on model order reduction, comprising:

[0058] A first acquisition module is used to obtain a response frequency of the engineering structure system in response to an excitation signal of a preset frequency band;

[0059] A second acquisition module is used to acquire a structural acoustic coupling system equation of the engineering structure system according to the response frequency;

[0060] A processing module is used to analyze and process the response frequency based on a reduced-order model corresponding to the structural-acoustic coupling system equation to obtain a target structural displacement and a target sound field pressure of the engineering structure system.

[0061] A third aspect of an embodiment of the present application provides a terminal device, comprising a memory and a processor, wherein the memory stores a computer program that can be run on the processor, and when the processor executes the computer program, the steps of the structural acoustic coupling processing method based on model order reduction as described in any one of the first aspects above are implemented.

[0062] A fourth aspect of an embodiment of the present application provides a computer-readable storage medium, comprising: storing a computer program, characterized in that when the computer program is executed by a processor, the steps of the structural acoustic coupling processing method based on model order reduction as described in any one of the first aspects above are implemented.

[0063] Compared with the prior art, the embodiments of the present application have the following advantages: by constructing a corresponding reduced-order model for the structural-acoustic coupling system equation, accelerated calculation is performed when analyzing and processing the response frequency in the engineering structure system, thereby improving the efficiency of solving and calculating the target structure displacement and the target sound field sound pressure in the engineering structure system. Therefore, in the engineering structure system, the target structure displacement and the target sound field sound pressure can be timely regulated by adjusting the parameters in the structural-acoustic coupling system equation; and the embodiment of the present application has strong ease of use and practicality. BRIEF DESCRIPTION OF THE DRAWINGS

[0064] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following briefly introduces the drawings required for use in the embodiments or descriptions of the prior art. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without any creative work.

[0065] Figure 1 Schematic diagram of an application scenario of the structural acoustic coupling processing method based on model order reduction provided in an embodiment of the present application;

[0066] Figure 2 1 is a schematic diagram of the implementation process of the structural acoustic coupling processing method based on model order reduction provided in an embodiment of the present application;

[0067] Figure 3 1 is a schematic diagram of an implementation flow of a structural acoustic coupling processing method based on model order reduction provided in an embodiment of the present application;

[0068] Figure 4 1 is a schematic diagram of the implementation process of the structural acoustic coupling processing method based on model order reduction provided in an embodiment of the present application;

[0069] Figure 5 Schematic diagram of the structure of a structural acoustic coupling processing device based on model order reduction provided in an embodiment of the present application;

[0070] Figure 6 It is a schematic diagram of the terminal device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0071] In the following description, specific details such as specific system structures and techniques are provided for purposes of illustration rather than limitation to facilitate a thorough understanding of the embodiments of the present application. However, it will be apparent to those skilled in the art that the present application may be implemented in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted to avoid obscuring the description of the present application with unnecessary detail.

[0072] like Figure 1 The figure shows a schematic diagram of an application scenario of the structural acoustic coupling processing method based on model order reduction provided by an embodiment of the present application. A steady-state structure is immersed in a homogeneous fluid. The structural vibration in the homogeneous fluid generates a radiated sound field. At the same time, the sound pressure of the sound field reacts on the structure, causing additional vibration of the structure and generating structural displacement. The structural and fluid dynamics equations must be combined to solve the structural acoustic coupling system equations to obtain the sound field pressure and structural displacement of the engineering structure system.

[0073] When addressing the coupling of structural vibration and fluid acoustic fields, sweep frequency analysis of the structural-acoustic coupling equations of the engineering structure system is required to determine the target acoustic pressure and displacement of the engineering structure system. Existing methods, most of which use iterative or direct methods to solve these targets, require extensive computation, take long calculation times, and suffer from low processing efficiency. Consequently, when parameters in the engineering structure system are adjusted, it is difficult to timely adjust the target displacement and pressure in the system.

[0074] This application proposes a structural acoustic coupling processing method based on model order reduction. It accelerates the analysis and processing of the response frequency through the reduced-order model corresponding to the structural acoustic coupling system equation, thereby improving the efficiency of solving the target structural displacement and target sound field sound pressure in the engineering structure system. The response frequency in the engineering structure system is analyzed and processed to obtain the target structural displacement and target sound field sound pressure. Images can be constructed based on the response frequency and the target structural displacement, and the response frequency and the target sound field sound pressure, respectively. It can estimate the dynamic performance of the engineering structure system, perform topological optimization design on the structure, identify the parameters of the engineering structure system, and perform acoustic sensitivity analysis on the engineering structure system. In addition, the structural acoustic coupling processing method based on model order reduction can adjust the parameters of the structural acoustic coupling system equation in the engineering structure system to adjust the target structural 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 controlling the noise of the sound field.

[0075] In order to illustrate the technical solution described in this application, specific embodiments are provided below.

[0076] Figure 2The following is a flowchart illustrating an implementation of a structural acoustic coupling processing method based on model order reduction provided by an embodiment of the present application, which is described in detail as follows: A structural acoustic coupling processing method based on model order reduction includes:

[0077] Step 101, obtaining a response frequency of the engineering structure system in response to an excitation signal of a preset frequency band;

[0078] Step 102: obtaining a structural acoustic coupling system equation of the engineering structure system according to the response frequency;

[0079] Step 103 : Analyze and process the response frequency based on the reduced-order model corresponding to the structural-acoustic coupling system equation to obtain a target structural displacement and a target sound field pressure of the engineering structure system.

[0080] In the above steps, a corresponding reduced-order model is constructed based on the structural-acoustic coupling system equation, and the analysis and processing process of the response frequency is accelerated to improve the computational efficiency of solving the engineering structure system to obtain the target structural displacement and target sound field sound pressure. After the structural-acoustic coupling system equation is adjusted in the engineering structure system, the target structural displacement and target sound field sound pressure can be timely regulated.

[0081] The engineering structure system includes a homogeneous fluid and a structure immersed in the homogeneous fluid. In the engineering structure system, the structure generates a target sound field pressure due to vibration. Simultaneously, the target sound field pressure reacts on the structure, causing additional structural vibration and generating target structural displacement.

[0082] In one embodiment, the excitation signal may be an excitation signal of a vibrator, which may be a mechanical vibrator, an electromagnetic vibrator, or an electro-hydraulic vibrator. In this case, the vibrator is used to excite the structure in the engineering structure system to generate vibrations, thereby generating an excitation signal.

[0083] In one embodiment, the excitation signal is a signal collected by a signal collector. In the engineering structure system, the signal collector is used to collect signals from the vibrating structure to obtain a collected signal, namely, the excitation signal.

[0084] In one embodiment, in step 101, specifically, in response to an excitation signal of a preset frequency band, a frequency domain analysis is performed on the excitation signal to obtain a response frequency of the engineering structure system.

[0085] In one embodiment, in step 102, obtaining a structural acoustic coupling system equation of the engineering structure system according to the response frequency specifically includes:

[0086] Step 201: Acquire first basic parameters of the engineering structure system, wherein the first basic parameters include structural stiffness, structural damping, structural mass, structural load, acoustic load, homogeneous fluid density, a first coupling term, and a second coupling term.

[0087] Step 202: Establish a geometric model of the engineering structure system based on the first basic parameters.

[0088] Step 203: Perform triangular meshing on the geometric model to obtain second basic parameters of the engineering structure system, wherein the second basic parameters include source points, field points, the interface between the structure and the acoustic field, the normal direction of the interface between the structure and the acoustic field, and the average area of the acoustic boundary element triangular mesh.

[0089] In the engineering structure system, in step 203, the geometric model is triangularly meshed to determine a first geometric model corresponding to the acoustic field in the engineering structure system. Based on the first geometric model of the acoustic field, each acoustic boundary element is determined. Based on all acoustic boundary elements, the center point of any one acoustic boundary element is determined as an acoustic boundary element collocation point, and the acoustic boundary element collocation point is used as a source point. The center points of the remaining acoustic boundary elements are determined as field points.

[0090] 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, wherein the third basic parameters include a first boundary element coefficient matrix and a second boundary element coefficient matrix.

[0091] In one embodiment, step 204 specifically includes determining a truncation order based on the size of the geometric model and the range of the response frequency, and determining a third basic parameter of the engineering structure system based on the truncation order, the response frequency, and the second basic parameter.

[0092] Specifically, based on the second basic parameter and the response frequency, a first boundary element kernel function and a second boundary element kernel function are determined, wherein the first boundary element kernel function is expressed as: The second boundary element kernel function is expressed as: Wherein, k is the wave number, k is expressed as k = ω / c, ω is the response frequency, c is the speed of sound in air, and j is an imaginary unit; x is the source point, y is the field point, r is the distance between the source point and the field point, and r is expressed as r = |xy|; n(y) is the normal direction of the interface between the structure and the sound field.

[0093] And, based on the truncation order, the first boundary element kernel function and the second boundary element kernel function are respectively expanded by Taylor's theorem to determine the first boundary element coefficient matrix and the second boundary element coefficient matrix, wherein the first boundary element coefficient matrix is expressed as The second boundary element coefficient matrix is expressed as Where, ξ=jk, L is the truncation order, Γ a is the boundary surface between the structure and the sound field. Specifically, based on the truncation order, the exponential coefficients of the first boundary element kernel function and the exponential coefficients of the second boundary element kernel function are expanded using Taylor's theorem, and the expanded 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.

[0094] In the above steps, based on the truncation order, the first boundary element kernel function and the second boundary element kernel function are respectively expanded using Taylor's theorem to obtain the first boundary element coefficient matrix and the second boundary element coefficient matrix. In the prior art, when analyzing and processing the response frequency (i.e., frequency sweep calculation), it is necessary to re-integrate and calculate each frequency of interest to form the first boundary element coefficient matrix and the second boundary element coefficient matrix, which is cumbersome and computationally intensive. In this method, only one integration process is required 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 shortcomings of the traditional coupling processing scheme that requires multiple integration calculations using the first boundary element kernel function and the second boundary element kernel function, resulting in low processing efficiency.

[0095] Step 205: Determine the structural acoustic coupling system equation based on the first basic parameter, the second basic parameter, and the third basic parameter. The structural acoustic coupling system equation is expressed as follows:

[0096]

[0097] Among them, K is the structural stiffness, D is the structural damping, M is the structural mass, and f s is the structural load, f a is the acoustic load, ρ a is the density of the homogeneous fluid, C sa is the first coupling term, C as is the second coupling term, u is the displacement of the target structure, and p is the sound pressure of the target sound field.

[0098] In the above steps, the first fundamental parameters of the engineering structure system are first obtained. Based on the first fundamental parameters, a geometric model of the engineering structure system is established. The geometric model is then triangularly meshed to obtain the second fundamental parameters, which are then used to construct a reduced-order model corresponding to the structural-acoustic coupling system equation. Based on the geometric model, the response frequency, and the second fundamental parameters, the first and second boundary element coefficient matrices are determined. Thus, the structural-acoustic coupling system equation of the engineering structure system is determined using the first, second, and third fundamental parameters, allowing the subsequent construction of a global orthogonal basis for structural finite elements and an acoustic boundary element. The reduced-order model is then used to improve the processing efficiency of the engineering structure system.

[0099] In one embodiment, the reduced-order model corresponding to formula (1) is used to analyze and process the response frequency to obtain the target structural displacement and target sound field pressure of the engineering structure system, including:

[0100] Step 301: According to formula (1), a global orthogonal basis of the structural finite element and a global orthogonal basis of the acoustic boundary element are obtained.

[0101] 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:

[0102]

[0103] Among them, K n is the second structural stiffness, K n Expressed as: K n =V H KV;

[0104] D n is the second structural damping, D n Expressed as: D n =V H DV;

[0105] M n is the second structural mass, M n Expressed as: M n =V H MV;

[0106] f s,n is the second structural load, f s,n Expressed as: f s,n =V H f s ;

[0107] f a,m is the second acoustic load, f a,mExpressed as: f a,m =W H f a ;

[0108] C sa,nm is the third coupling term, C sa,nm Expressed as: C sa,nm =V H C sa W;

[0109] C as,n is the fourth coupling term, C as,n Expressed as: C as,n =C as V;

[0110] u n is the second structural displacement; p m is the sound pressure of the second sound field;

[0111] G m is the third boundary element coefficient matrix, G m Expressed as: in,

[0112] H m is the fourth boundary element coefficient matrix, H m Expressed as: in, Wherein, 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 H represents the conjugate transpose.

[0113] Step 303: Determine the target structural displacement and the target sound field pressure based on the structural finite element global orthogonal basis, the acoustic boundary element global orthogonal basis, the second structural displacement, and the second sound field pressure. The target structural displacement is expressed as: u = Vu n , the target sound field sound pressure is expressed as: p = Wp m .

[0114] In the above steps, global orthogonal bases for the structural finite element and the acoustic boundary element are obtained based on the structural-acoustic coupling system equations. Based on these global orthogonal bases, the acoustic boundary element, and the reduced-order model, the reduced-order second structural displacement and second acoustic field sound pressure are determined. This avoids directly processing the original structural-acoustic coupling system equations and determines the reduced-order second structural displacement and second acoustic field sound pressure through the reduced-order model, thereby improving the processing efficiency of the engineering structure system. Furthermore, based on the second structural displacement, second acoustic field sound pressure, the global orthogonal bases for the structural finite element and the acoustic boundary element, the target structural displacement and target acoustic field sound pressure of the engineering structure system are obtained. This structural-acoustic coupling processing method based on model order reduction can adjust the parameters of the engineering structure system to timely adjust the target structural displacement and target acoustic field sound pressure of the engineering structure system, thereby controlling the vibration of the structure within the engineering structure system and controlling the noise of the acoustic field, thereby achieving subsequent structural optimization design of the engineering structure system under the constraints of structural and acoustic loads.

[0115] In this structural acoustic coupling processing method based on model order reduction, the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained according to formula (1). This application provides four embodiments for obtaining the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis under different constraints of structural loads and acoustic loads:

[0116] Example 1

[0117] In the first embodiment, according to formula (1), the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained, which specifically include:

[0118] When the structural load is not zero and the acoustic load is zero, according to formula (1), a first structural finite element global orthogonal basis and a first acoustic boundary element global orthogonal basis are obtained;

[0119] The first structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the first acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

[0120] In one embodiment, according to formula (1), the first structural finite element global orthogonal basis and the first acoustic boundary element global orthogonal basis are obtained, which specifically include:

[0121] Step 401: Ignore the first coupling term in formula (1), and calculate the first structural finite element global orthogonal basis using a second-order Arnoldi algorithm.

[0122] In step 401, specifically, the first coupling term of formula (1) is ignored to obtain a first structural acoustic coupling system equation, wherein the first structural acoustic coupling system equation is expressed as:

[0123]

[0124] Among them, K is the structural stiffness, D is the structural damping, M is the structural mass, and f s is the structural load, u is the target structural displacement, and ω is the response frequency.

[0125] According to formula (3), the first structure finite element global orthogonal basis is calculated by the second-order Arnoldi algorithm, specifically including:

[0126] According to formula (3), the first matrix, the second matrix, and the first starting vector are obtained, wherein the first matrix is expressed as: The second matrix is expressed as: The first starting vector is expressed as: in, s * =jω * is the selected expansion point, ω is the response frequency, and j is the imaginary unit;

[0127] Recursively and modified Gram-Schmidt orthogonalization processing is performed on the first matrix, the second matrix, and the first starting vector to obtain a set of first structure finite element local orthogonal bases; based on the set of first structure finite element local orthogonal bases, the first structure finite element local orthogonal bases are generated for different extension points; the first structure finite element local orthogonal bases are re-orthogonalized to obtain the first structure finite element global orthogonal base.

[0128] Step 402: Based on the second basic parameter, the first boundary element coefficient matrix and the second boundary element coefficient matrix of formula (1) are truncated to obtain a first sparse coefficient matrix and a second sparse coefficient matrix. Specifically, based on the second basic parameter and according to the characteristics of the first boundary element kernel function and the second boundary element kernel function attenuating with distance, the first boundary element coefficient matrix and the second boundary element coefficient matrix of formula (1) are truncated to obtain a first sparse coefficient matrix and a second sparse coefficient matrix.

[0129] Specifically, the cutoff radius is determined based on the average area of the acoustic boundary element triangle mesh, wherein the cutoff radius can be expressed as mean(A i ) is the average area of the acoustic boundary element triangular mesh.

[0130] In the engineering structure system, in the aforementioned step 203, the geometric model is subjected to triangular meshing processing to determine a first geometric model corresponding to the sound field in the engineering structure system. Based on the first geometric model of the sound field, a plurality of triangular units are obtained. For all triangular units, the triangular units are respectively used as the aforementioned acoustic boundary units. For all acoustic boundary units, the centers of the acoustic boundary units are respectively used as boundary element points. Among all acoustic boundary units, one of the acoustic boundary units is used as a source unit, the boundary element point of the acoustic boundary unit is used as a source point, and the source point is used as the center of the sphere; the boundary element points of the remaining acoustic boundary units are used as field points.

[0131] Based on the truncation radius and the center of the sphere, a truncated sphere geometric model is established; based on the truncated sphere geometric model, the acoustic boundary units located within the truncated sphere are determined to be strong interaction units; and the kd tree search algorithm is used to determine the number of adjacent strong interaction units; based on the first boundary element coefficient matrix, the second boundary element coefficient matrix, and the number of adjacent strong interaction units, the first sparse coefficient matrix and the second sparse coefficient matrix are determined.

[0132] Step 403: Based on the first sparse coefficient matrix and the second sparse coefficient matrix, the first acoustic boundary element global orthogonal basis is calculated by using an Arnoldi algorithm or an eigenorthogonal decomposition algorithm.

[0133] Specifically, based on the first sparse coefficient matrix and the second sparse coefficient matrix, a second structural acoustic coupling system equation is obtained, wherein the second structural acoustic coupling system equation is expressed as:

[0134]

[0135] Among them, ρ a is the density of the homogeneous fluid, C as is the second coupling term, ω is the response frequency, p is the target sound field pressure, is the first sparse coefficient matrix, is the second sparse coefficient matrix, with the superscript ◇ The representation is a sparse approximation of the original dense matrix.

[0136] In one embodiment, according to formula (4), the first acoustic boundary element global orthogonal basis is calculated using the Arnoldi algorithm.

[0137] Specifically, according to formula (4), a third matrix and a second starting vector are obtained, wherein the third matrix is expressed as: The second starting vector is expressed as: According to the third matrix and the second starting vector, 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 acoustic boundary element local orthogonal bases are generated for different extension points; the first acoustic boundary element local orthogonal bases are re-orthogonalized to obtain the first acoustic boundary element global orthogonal base.

[0138] Example 2

[0139] In the second embodiment, according to formula (1), the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained, which specifically include:

[0140] When the structural load is zero and the acoustic load is not zero, according to formula (1), a second structural finite element global orthogonal basis and a second acoustic boundary element global orthogonal basis are obtained;

[0141] The second structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the second acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

[0142] In one embodiment, according to formula (1), the second structural finite element global orthogonal basis and the second acoustic boundary element global orthogonal basis are obtained, which specifically include:

[0143] Step 501: Based on the second basic parameters, the first boundary element coefficient matrix and the second boundary element coefficient matrix of formula (1) are truncated to obtain a third sparse coefficient matrix and a fourth sparse coefficient matrix.

[0144] Step 502: Based on the third sparse coefficient matrix and the fourth sparse coefficient matrix, the second coupling term of formula (1) is ignored, and the second acoustic boundary element global orthogonal basis is calculated by Arnoldi algorithm or eigenorthogonal decomposition algorithm.

[0145] In steps 501 to 502, specifically, based on the third sparse coefficient matrix and the fourth sparse coefficient matrix, the second coupling term of formula (1) is ignored to obtain a third structural acoustic coupling system equation, wherein the third structural acoustic coupling system equation is expressed as:

[0146]

[0147] Where p is the sound pressure of the target sound field, is the third sparse coefficient matrix, f a is the acoustic load, ◇ The representation is a sparse approximation of the original dense matrix.

[0148] In one embodiment, according to formula (5), the second acoustic boundary element global orthogonal basis is calculated using the Arnoldi algorithm.

[0149] Specifically, according to formula (5), the fourth matrix and the third starting vector are obtained, wherein the fourth matrix is expressed as: The third starting vector is expressed as: According to the fourth matrix and the third starting vector, 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 second acoustic boundary element local orthogonal bases are generated for different extension points; the second acoustic boundary element local orthogonal bases are re-orthogonalized to obtain the second acoustic boundary element global orthogonal base.

[0150] Step 503: Obtain a global orthogonal basis of the second structure finite element by using a second-order Arnoldi algorithm.

[0151] In step 503, specifically, a fourth structural acoustic coupling system equation is obtained according to the third starting vector, wherein the fourth structural acoustic coupling system equation is expressed as:

[0152]

[0153] Where K is the structural stiffness, D is the structural damping, M is the structural mass, u is the target structural displacement, ω is the response frequency, C sa is the first coupling term.

[0154] According to formula (6), the second structure finite element global orthogonal basis is calculated by the second-order Arnoldi algorithm.

[0155] Example 3

[0156] In the third embodiment, according to formula (1), the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained, which specifically include:

[0157] When the structural load and the acoustic load are both non-zero, according to formula (1), the third structural finite element global orthogonal basis and the third acoustic boundary element global orthogonal basis are obtained;

[0158] The third structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the third acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

[0159] In one embodiment, according to formula (1), a third structural finite element global orthogonal basis and a third acoustic boundary element global orthogonal basis are obtained, specifically including:

[0160] Step 601: Ignore the first coupling term in formula (1), and calculate the third structure finite element global orthogonal basis using the second-order Arnoldi algorithm.

[0161] In step 601, specifically, the first coupling term of formula (1) is ignored to obtain the fifth structural acoustic coupling system equation, wherein the fifth structural acoustic coupling system equation is expressed as:

[0162]

[0163] Among them, K is the structural stiffness, D is the structural damping, M is the structural mass, and f s is the structural load, u is the target structural displacement, ω is the response frequency, and j is the imaginary unit;

[0164] According to formula (7), the third structure finite element global orthogonal basis is calculated by the second-order Arnoldi algorithm.

[0165] Specifically, according to formula (7), the fifth matrix, the sixth matrix, and the fourth starting vector are obtained, wherein the fifth matrix is expressed as: The sixth matrix is expressed as: The fourth starting vector is expressed as: in, s * =jω * is the selected expansion point, and ω is the response frequency. A set of local orthogonal bases of the third structure finite element is obtained based on the fifth matrix, the sixth matrix, and the fourth starting vector. Based on the set of local orthogonal bases of the third structure finite element, the local orthogonal bases of the third structure finite element are generated for different expansion points. The local orthogonal bases of the third structure finite element are re-orthogonalized to obtain the global orthogonal bases of the third structure finite element.

[0166] Step 602: Based on the second basic parameters, the first boundary element coefficient matrix and the second boundary element coefficient matrix of formula (1) are truncated to obtain a fifth sparse coefficient matrix and a sixth sparse coefficient matrix.

[0167] Step 603: Based on the fifth sparse coefficient matrix and the sixth sparse coefficient matrix, obtain the third acoustic boundary element global orthogonal basis by calculating using an Arnoldi algorithm or an eigenorthogonal decomposition algorithm.

[0168] In one embodiment, in 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 truncated to obtain a fifth sparse coefficient matrix and a 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, wherein the sixth structural acoustic coupling system equation is expressed as:

[0169]

[0170] Among them, ρ a is the density of the homogeneous fluid, C as is the second coupling term, ω is the response frequency, p is the target sound field pressure, j is the imaginary unit, is the fifth sparse coefficient matrix, is the sixth sparse coefficient matrix, f a is the acoustic load, ◇ The representation is a sparse approximation of the original dense matrix.

[0171] In one embodiment, according to formula (8), the third acoustic boundary element global orthogonal basis is calculated by using the Arnoldi algorithm.

[0172] Example 4

[0173] In the fourth embodiment, according to formula (1), the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis are obtained, which specifically include:

[0174] When the structural load and the acoustic load are both non-zero, according to formula (1), the fourth structural finite element global orthogonal basis and the fourth acoustic boundary element global orthogonal basis are obtained;

[0175] The fourth structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the fourth acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

[0176] In one embodiment, according to formula (1), the fourth structural finite element global orthogonal basis and the fourth acoustic boundary element global orthogonal basis are obtained, which specifically include:

[0177] Step 701: 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 seventh sparse coefficient matrix and an eighth sparse coefficient matrix;

[0178] 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 acoustic boundary element global orthogonal basis is calculated by the Arnoldi algorithm or the eigenorthogonal decomposition algorithm.

[0179] In one embodiment, in 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 truncated to obtain a seventh sparse coefficient matrix and an eighth sparse coefficient matrix;

[0180] Based on the seventh sparse coefficient matrix and the eighth sparse coefficient matrix, the second coupling term of formula (1) is ignored to obtain the seventh structural acoustic coupling system equation, wherein the seventh structural acoustic coupling system equation is expressed as:

[0181]

[0182] Where p is the sound pressure of the target sound field, is the seventh sparse coefficient matrix, f a is the acoustic load, ◇ The representation is a sparse approximation of the original dense matrix.

[0183] In one embodiment, according to formula (9), the fourth acoustic boundary element global orthogonal basis is calculated by using the Arnoldi algorithm.

[0184] Specifically, according to formula (9), the seventh matrix and the fifth starting vector are obtained, wherein the seventh matrix is expressed as: The fifth starting vector is expressed as: According to the seventh matrix and the fifth starting vector, a set of fourth acoustic boundary element local orthogonal bases is obtained; based on the set of fourth acoustic boundary element local orthogonal bases, the fourth acoustic boundary element local orthogonal bases are generated for different extension points; the fourth acoustic boundary element local orthogonal bases are re-orthogonalized to obtain the fourth acoustic boundary element global orthogonal base.

[0185] Step 703: Obtain the fourth structure finite element global orthogonal basis by calculating using a second-order Arnoldi algorithm.

[0186] In step 703, specifically, an eighth structural acoustic coupling system equation is obtained according to the fifth starting vector, wherein the eighth structural acoustic coupling system equation is expressed as:

[0187]

[0188] Where K is the structural stiffness, D is the structural damping, M is the structural mass, u is the target structural displacement, ω is the response frequency, j is the imaginary unit, C sa is the first coupling term.

[0189] According to formula (10), the fourth structure finite element global orthogonal basis is calculated by the second-order Arnoldi algorithm.

[0190] In the above four embodiments, the construction of the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis under different structural load and acoustic load conditions is respectively shown. In the present application, embodiment one, embodiment two, embodiment three and embodiment four are exemplarily shown, wherein embodiment one shows the specific load condition with only structural load, embodiment two shows the specific load condition with only acoustic load, and embodiment three and embodiment four show the specific load condition with both structural load and acoustic load. It can be seen that the structural acoustic coupling processing method based on model order reduction can meet the processing requirements under different load conditions, by truncating the first boundary element coefficient matrix and the second boundary element coefficient matrix, and by choosing to ignore the first coupling term or the second coupling term of the structural acoustic coupling system equation, so as to reduce the amount of data required to process the structural acoustic coupling system equation, and improve the efficiency of constructing the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis.

[0191] It should be understood that the size of the serial numbers of the steps in the above embodiments 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 embodiments of this application.

[0192] Corresponding to the method of the above embodiment, Figure 5 A structural block diagram of a structural acoustic coupling processing device based on model order reduction provided in an embodiment of the present application is shown. For ease of explanation, only the parts related to the embodiment of the present application are shown. Figure 5 The exemplary structural acoustic coupling processing apparatus based on model order reduction may be an execution body of the structural acoustic coupling processing method based on model order reduction provided in the aforementioned embodiment.

[0193] Reference Figure 5 The structural acoustic coupling processing device 50 based on model order reduction includes:

[0194] A first acquisition module 51 is configured to acquire a response frequency of the engineering structure system in response to an excitation signal in a preset frequency band;

[0195] A second acquisition module 52 is configured to acquire a structural acoustic coupling system equation of the engineering structure system according to the response frequency;

[0196] The processing module 53 is configured to analyze and process the response frequency based on a reduced-order model corresponding to the structural-acoustic coupling system equation to obtain a target structural displacement and a target sound field pressure of the engineering structure system.

[0197] The process of each module in the structural acoustic coupling processing device based on model order reduction provided in the embodiment of the present application realizing its own function can be specifically referred to the aforementioned Figure 1 The description of the embodiment of the method shown will not be repeated here.

[0198] It should be understood that when used in the present specification and the appended claims, the term "comprising" indicates the presence of described features, integers, steps, operations, elements and / or components, but does not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or collections thereof.

[0199] It will also be understood that the term "and / or" used in this specification and the appended claims refers to and includes any and all possible combinations of one or more of the associated listed items.

[0200] As used in this specification and the appended claims, the term "if" can be interpreted as "when" or "upon" or "in response to determining" or "in response to detecting," depending on the context. Similarly, the phrase "if it is determined" or "if [described condition or event] is detected" can be interpreted as meaning "upon determination" or "in response to determining" or "upon detection of [described condition or event]" or "in response to detecting [described condition or event]," depending on the context.

[0201] In addition, in the description of the present specification and the appended claims, the terms "first", "second", "third", etc. are only used to distinguish descriptions and should not be understood as indicating or implying relative importance. It should also be understood that although the terms "first", "second", etc. are used in the text to describe various elements in some embodiments of the present application, these elements should not be limited by these terms. These terms are only used to distinguish one element from another. For example, a first table can be named a second table, and similarly, a second table can be named a first table without departing from the scope of the various described embodiments. Both the first table and the second table are tables, but they are not the same table.

[0202] References to "one embodiment" or "some embodiments" in this specification mean that a particular feature, structure, or characteristic described in conjunction with that embodiment is included in one or more embodiments of the present application. Thus, phrases such as "in one embodiment," "in some embodiments," "in other embodiments," and "in other embodiments" appearing in various places in this specification do not necessarily refer to the same embodiment, but rather mean "one or more but not all embodiments," unless otherwise specifically emphasized. The terms "including," "comprising," "having," and variations thereof all mean "including but not limited to," unless otherwise specifically emphasized.

[0203] Figure 6 This is a schematic diagram of the structure of a terminal device provided by an embodiment of the present application. Figure 6 As shown, the terminal device 60 of this embodiment includes: at least one processor 61 ( Figure 6 Only one is shown), a memory 62, wherein the memory 62 stores a computer program 63 that can be run on the processor 61. When the processor 61 executes the computer program 63, the steps in the above-mentioned various embodiments of the structural acoustic coupling processing method based on model order reduction are implemented, such as Figure 2 Alternatively, when the processor 61 executes the computer program 63, the functions of the modules / units in the above-mentioned device embodiments are realized, for example, Figure 5 The functions of modules 51 to 53 are shown.

[0204] The terminal device 60 can be a computing device such as a desktop computer, a notebook, a PDA, or a cloud server. The terminal device can include, but is not limited to, a processor 61 and a memory 62. Those skilled in the art will understand that Figure 6 It is only an example of the terminal device 60 and does not constitute a limitation of the terminal device 60. It may include more or fewer components than shown in the figure, or a combination of certain components, or different components. For example, the terminal device may also include an input and sending device, a network access device, a bus, etc.

[0205] The processor 61 may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any conventional processor.

[0206] In some embodiments, 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 media card (SMC), a secure digital (SD) card, a flash card, etc. equipped on the terminal device 60. Furthermore, the memory 62 may include both an internal storage unit of the terminal device 60 and an external storage device. The memory 62 is used to store an operating system, application programs, a boot loader, data, and other programs, such as the program code of the computer program. The memory 62 may also be used to temporarily store data that has been sent or is about to be sent.

[0207] In addition, the functional units in the various embodiments of the present application may be integrated into a single processing unit, or each unit may exist physically separately, or two or more units may be integrated into a single unit. The aforementioned integrated units may be implemented in the form of hardware or software functional units.

[0208] 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. When the processor executes the computer program, the terminal device implements the steps of any of the above-mentioned method embodiments.

[0209] An embodiment of the present application further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps in the above-mentioned various method embodiments can be implemented.

[0210] If the integrated module / unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the present application implements all or part of the process in the above-mentioned embodiment method, and can also be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium, and the computer program can implement the steps of the above-mentioned various method embodiments when executed by the processor. Wherein, the computer program includes computer program code, and the computer program code can be in source code form, object code form, executable file or some intermediate form, etc. The computer-readable medium may include: any entity or device that can carry the computer program code, recording medium, U disk, mobile hard disk, magnetic disk, optical disk, computer memory, read-only memory (ROM), random access memory (RAM), electric carrier signal, telecommunication signal and software distribution medium, etc.

[0211] In the above embodiments, the description of each embodiment has its own focus. For parts that are not described or recorded in detail in a certain embodiment, reference can be made to the relevant description of other embodiments.

[0212] Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0213] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of these units may be selected to achieve the purpose of this embodiment according to actual needs.

[0214] The above-described embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. These modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present application, and should all be included in the scope of protection of the present application.

Claims

1. A structural acoustic coupling processing method based on model order reduction, characterized in that: include: In response to an excitation signal in a preset frequency band, a response frequency of the engineering structure system is obtained; Obtaining a structural acoustic coupling system equation of the engineering structure system according to the response frequency; Analyzing and processing the response frequency based on a reduced-order model corresponding to the structural-acoustic coupling system equation to obtain a target structural displacement and a target sound field pressure of the engineering structure system; The obtaining of the structural acoustic coupling system equation of the engineering structure system according to the response frequency specifically includes: Acquire first basic parameters of the engineering structure system, wherein the first basic parameters include structural stiffness, structural damping, structural mass, structural load, acoustic load, homogeneous fluid density, first coupling term, and second coupling term; establish a geometric model of the engineering structure system based on the first basic parameters; perform triangular meshing on the geometric model to obtain second basic parameters of the engineering structure system, wherein the second basic parameters include source points, field points, the interface between the structure and the acoustic field, the normal direction of the interface between the structure and the acoustic field, and the average area of the acoustic boundary element triangular mesh; determine a truncation order according to the size of the geometric model and the range of the response frequency; determine a first boundary element kernel function and a second boundary element kernel function based on the second basic parameters and the response frequency; based on the truncation order, expand the first boundary element kernel function and the second boundary element kernel function respectively by Taylor's theorem to determine third basic parameters of the engineering structure system, wherein the third basic parameters include a first boundary element coefficient matrix and a second boundary element coefficient matrix; determine the structural acoustic coupling system equation according to the first basic parameters, the second basic parameters, and the third basic parameters; According to the structural-acoustic coupling system equation, the global orthogonal basis of the structural finite element and the global orthogonal basis of the acoustic boundary element are obtained, which specifically include: When the structural load is not zero and the acoustic load is zero, the first coupling term of the structural-acoustic coupling system equation is ignored, and a first structural finite element global orthogonal basis is calculated by a second-order Arnoldi algorithm; based on the second basic parameter, the first boundary element coefficient matrix and the second boundary element coefficient matrix of the structural-acoustic coupling system equation are truncated to obtain a first sparse coefficient matrix and a second sparse coefficient matrix; based on the first sparse coefficient matrix and the second sparse coefficient matrix, a first acoustic boundary element global orthogonal basis is calculated by an Arnoldi algorithm or an eigenorthogonal decomposition algorithm; the first structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the first acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

2. The structural acoustic coupling processing method based on model order reduction according to claim 1, characterized in that: The first boundary element kernel function is expressed as: , the second boundary element kernel function is expressed as: ,in, is the wave number, Expressed as , is the response frequency, is the speed of sound in air, is an imaginary unit; is the source point, For the field point, The source point With the site The distance between Expressed as ; is the normal direction of the interface between the structure and the sound field; Among them, the first boundary element coefficient matrix is expressed as , ; The second boundary element coefficient matrix is expressed as , ,in, , is the truncation order, is the boundary surface between the structure and the sound field.

3. The structural acoustic coupling processing method based on model order reduction according to claim 2, characterized in that: The analyzing and processing 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 pressure of the engineering structure system includes: According to the structural-acoustic coupling system equation, a global orthogonal basis of a structural finite element and a global orthogonal basis of an acoustic boundary element are obtained; Determining a second structural displacement and a 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; Determining the target structural displacement and the target sound field pressure 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; The reduced-order model is expressed as follows: ， in, is the second structural stiffness, Expressed as: ; is the second structural damping, Expressed as: ; is the second structural mass, Expressed as: ; is the second structural load, Expressed as: ; is the second acoustic load, Expressed as: ; is the third coupling term, Expressed as: ; is the fourth coupling term, Expressed as: ; is the third boundary element coefficient matrix, Expressed as: ,in, ; is the fourth boundary element coefficient matrix, Expressed as: ,in, , where the superscript represents the conjugate transpose; and, is the density of the homogeneous fluid; is the second structural displacement; is the sound pressure of the second sound field; The target structure displacement is expressed as: , the target sound field sound pressure is expressed as: ; in, is the global orthogonal basis of the structural finite element, is the acoustic boundary element global orthogonal basis; is the structural damping, For the structural quality, is the structural load, is the acoustic load.

4. The structural acoustic coupling processing method based on model order reduction according to claim 3, characterized in that: The step of obtaining the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis according to the structural-acoustic coupling system equation specifically includes: When the structural load is zero and the acoustic load is not zero, truncating the first boundary element coefficient matrix and the second boundary element coefficient matrix of the structural acoustic coupling system equation based on the second basic parameter to obtain a third sparse coefficient matrix and a fourth sparse coefficient matrix; Based on the third sparse coefficient matrix and the fourth sparse coefficient matrix, the second coupling term of the structural acoustic coupling system equation is ignored, and a second acoustic boundary element global orthogonal basis is calculated by an Arnoldi algorithm or an eigenorthogonal decomposition algorithm; The second structure finite element global orthogonal basis is obtained by calculating the second-order Arnoldi algorithm; The second structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the second acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

5. The structural acoustic coupling processing method based on model order reduction according to claim 3, characterized in that: The step of obtaining the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis according to the structural-acoustic coupling system equation specifically includes: When both the structural load and the acoustic load are non-zero, the first coupling term of the structural-acoustic coupling system equation is ignored, and a third structural finite element global orthogonal basis is obtained by calculation using a second-order Arnoldi algorithm; Based on the second basic parameters, truncate the first boundary element coefficient matrix and the second boundary element coefficient matrix of the structural acoustic coupling system equation to obtain a fifth sparse coefficient matrix and a sixth sparse coefficient matrix; Based on the fifth sparse coefficient matrix and the sixth sparse coefficient matrix, a third acoustic boundary element global orthogonal basis is calculated by an Arnoldi algorithm or an eigenorthogonal decomposition algorithm; The third structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the third acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

6. The structural acoustic coupling processing method based on model order reduction according to claim 3, characterized in that: The step of obtaining the structural finite element global orthogonal basis and the acoustic boundary element global orthogonal basis according to the structural-acoustic coupling system equation specifically includes: When both the structural load and the acoustic load are non-zero, truncating the first boundary element coefficient matrix and the second boundary element coefficient matrix of the structural-acoustic coupling system equation based on the second basic parameter to obtain a seventh sparse coefficient matrix and an eighth sparse coefficient matrix; Based on the seventh sparse coefficient matrix and the eighth sparse coefficient matrix, the second coupling term of the structural acoustic coupling system equation is ignored, and a fourth acoustic boundary element global orthogonal basis is calculated by an Arnoldi algorithm or an eigenorthogonal decomposition algorithm; The fourth structure finite element global orthogonal basis is obtained by calculating the second-order Arnoldi algorithm; The fourth structural finite element global orthogonal basis is used as the structural finite element global orthogonal basis, and the fourth acoustic boundary element global orthogonal basis is used as the acoustic boundary element global orthogonal basis.

7. A structural acoustic coupling processing device based on model order reduction, characterized in that: include: A first acquisition module is used to obtain a response frequency of the engineering structure system in response to an excitation signal of a preset frequency band; A second acquisition module is used to acquire a structural acoustic coupling system equation of the engineering structure system according to the response frequency; a processing module, configured to analyze and process the response frequency based on a reduced-order model corresponding to the structural-acoustic coupling system equation to obtain a target structural displacement and a target sound field pressure of the engineering structure system; The second acquisition module is also used to obtain the first basic parameters of the engineering structure system, wherein the first basic parameters include structural stiffness, structural damping, structural mass, structural load, acoustic load, homogeneous fluid density, first coupling term, and second coupling term; based on the first basic parameters, a geometric model of the engineering structure system is established; the geometric model is subjected to triangular mesh division processing to obtain the second basic parameters of the engineering structure system, wherein the second basic parameters include source points, field points, the interface between the structure and the sound field, the normal direction of the interface between the structure and the sound field, and the average area of the acoustic boundary element triangular mesh; according to the size of the geometric model and the range of the response frequency, the truncation order is determined; based on the second basic parameters and the response frequency, the first boundary element kernel function and the second boundary element kernel function are determined; based on the truncation order, the first boundary element kernel function and the second boundary element kernel function are respectively expanded by Taylor's theorem to determine the third basic parameters of the engineering structure system, wherein the third basic parameters include the first boundary element coefficient matrix, the second ... third boundary element coefficient matrix, the third boundary element coefficient matrix, the third boundary element coefficient matrix, the third boundary element coefficient matrix, the third boundary element coefficient matrix, the third boundary element coefficient matrix, the third boundary element coefficient matrix, the third boundary Boundary element coefficient matrix; determining the structural acoustic coupling system equation according to the first basic parameters, the second basic parameters and the third basic parameters; obtaining a structural finite element global orthogonal basis and an acoustic boundary element global orthogonal basis according to the structural acoustic coupling system equation, specifically including: when the structural load is not zero and the acoustic load is zero, ignoring the first coupling term of the structural acoustic coupling system equation, and calculating the first structural finite element global orthogonal basis by a second-order Arnoldi algorithm; truncating the first boundary element coefficient matrix and the second boundary element coefficient matrix of the structural acoustic coupling system equation based on the second basic parameters to obtain a first sparse coefficient matrix and a second sparse coefficient matrix; calculating the first acoustic boundary element global orthogonal basis based on the first sparse coefficient matrix and the second sparse coefficient matrix by an Arnoldi algorithm or an eigenorthogonal decomposition algorithm; using the first structural finite element global orthogonal basis as the structural finite element global orthogonal basis, and using the first acoustic boundary element global orthogonal basis as the acoustic boundary element global orthogonal basis.

8. A terminal device, characterized in that: The terminal device includes a memory and a processor, wherein the memory stores a computer program that can be run on the processor, and when the processor executes the computer program, the steps of the method according to any one of claims 1 to 6 are implemented.

9. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

Images