A method for fast estimation of acoustic radiation from a fluid domain structure

By using a neural network model to map pulsating pressure information in the estimation of acoustic radiation in the fluid domain structure, the problems of low computational efficiency and poor convergence in the existing technology are solved, achieving more efficient acoustic radiation estimation and reducing the requirements for the quality of surface element division.

CN121031459BActive Publication Date: 2026-02-03WUHAN UNIV OF TECH
View PDF 2 Cites 0 Cited by

Patent Information

Application Number
CN202511564245.1
Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
Filing Date
2025-10-30
Publication Date
2026-02-03
Estimated Expiration
2045-10-30

AI Technical Summary

Technical Problem

Existing techniques suffer from low computational efficiency and poor convergence in the estimation of acoustic radiation in structures containing fluid domains, and the need to establish accurate boundary layer meshes increases computational complexity.

Method used

By establishing a coupling surface between the structure and the fluid domain, and using a neural network model to map the pulsating pressure information on the fluid domain model to the structural nodes, the process bypasses the step of importing the flow field information into the structural elements as boundary conditions to solve the vibration response in existing technologies. Only the consistency of the surface elements needs to be ensured, which reduces the requirements for the quality of the surface element division.

Benefits of technology

It improves computational efficiency, reduces computational time costs, and enhances computational convergence, thus achieving more efficient acoustic radiation estimation.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN121031459B_ABST
    Figure CN121031459B_ABST
Patent Text Reader

Abstract

The present application relates to a kind of acoustic radiation fast estimation method of fluid domain structure, belong to acoustic radiation estimation technical field, wherein, the method includes: establishing structure model;Structure model includes the coupling surface of structure and fluid domain, coupling surface includes multiple first surface elements, and first surface element includes structure node;Establish fluid domain model;Fluid domain model also includes above-mentioned coupling surface;The coupling surface of fluid domain model includes the second surface element corresponding to first surface element;First pulsating pressure information on second surface element is obtained by fluid dynamics simulation;Second pulsating pressure information of each structure node is obtained by mapping first pulsating pressure information on second surface element to the structure node on the first surface element corresponding to second surface element by neural network model, by neural network model to the structure node on the first surface element corresponding to second surface element;Second pulsating pressure information is used to determine the acoustic radiation of structure.The present application improves the estimation efficiency and the convergence of calculation of the acoustic radiation of structure with fluid domain.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention relates to the field of acoustic radiation estimation technology, and in particular to a fast method for estimating acoustic radiation in structures containing fluid domains. Background Technology

[0002] In engineering scenarios such as ships, aircraft, and vehicles, structures containing fluid domains are often included. Accurate estimation of the sound radiation of such structures can enable noise control and acoustic optimization in engineering scenarios.

[0003] Existing techniques typically employ numerical analysis methods to estimate the acoustic radiation of structures containing fluid domains. Specifically, this involves first establishing finite element fluid meshes and structural meshes, then calculating fluid flow field information (pulsating pressure, velocity, etc.) and importing it into the structural mesh as boundary conditions. The vibration response of the structure is then solved, and finally, the acoustic radiation information is calculated based on the vibration response. However, this method usually requires accurate boundary layer meshes to ensure computational convergence. Establishing accurate boundary layer meshes necessitates specialized preprocessing software such as ICEM and Hypermesh, significantly increasing computational complexity. Therefore, in summary, existing techniques suffer from high computational time costs and poor convergence. Summary of the Invention

[0004] In view of this, it is necessary to provide a fast method for estimating acoustic radiation of structures containing fluid domains, so as to solve the problems of low computational efficiency and poor convergence of existing acoustic radiation calculation methods for structures containing fluid domains.

[0005] To address the above problems, this invention provides a fast method for estimating acoustic radiation in structures containing fluid domains, the method comprising:

[0006] Establish a structural model; the structural model includes a coupling surface between the structure and the fluid domain, the coupling surface includes multiple first surface elements, and the first surface elements include structural nodes;

[0007] A fluid domain model is established; the fluid domain model includes the coupling surface; the coupling surface of the fluid domain model includes multiple second surface elements that correspond one-to-one with the multiple first surface elements;

[0008] Fluid dynamics simulation is performed based on the fluid domain model to obtain the first pulsating pressure information on the plurality of second surface units;

[0009] The first pulsating pressure information on the second surface unit is mapped to the structural node on the first surface unit corresponding to the second surface unit through a neural network model, thereby obtaining the second pulsating pressure information of each structural node;

[0010] Based on the second pulsating pressure information, the vibration information of each structural node is determined, and the acoustic radiation of the structure is determined based on the vibration information.

[0011] In one possible implementation, the first pulsating pressure information includes a target pulsating pressure frequency domain signal; the step of performing computational fluid dynamics simulation based on the fluid domain model to obtain the first pulsating pressure information on the plurality of second surface elements includes:

[0012] Computational fluid dynamics simulation is performed based on the fluid domain model to obtain the time-domain signal of the pulsating pressure on the plurality of second surface units;

[0013] The pulsating pressure time-domain signal is converted into a pulsating pressure frequency-domain signal, and a signal within a preset frequency range is extracted from the pulsating pressure frequency-domain signal as the target pulsating pressure frequency-domain signal.

[0014] In one possible implementation, the preset frequency range is 2-300Hz.

[0015] In one possible implementation, the structural model is constructed using finite element analysis software, and the fluid domain model is constructed using fluid dynamics simulation software.

[0016] In one possible implementation, the neural network model is a neural network model based on radial basis functions.

[0017] In one possible implementation, determining the vibration information of each structural node based on the second pulsating pressure information includes:

[0018] The second pulsating pressure information of each structural node is input into the finite element analysis software to obtain the vibration information of each structural node.

[0019] In one possible implementation, determining the acoustic radiation of the structure based on the vibration information includes:

[0020] Based on the vibration information, determine the equivalent sound source inside the structure and the source strength of the equivalent sound source;

[0021] Based on the equivalent sound source and the source strength, the superimposed sound radiation of the equivalent sound source at the target field point is determined by combining the wave superposition method.

[0022] In one possible implementation, the vibration information includes vibration velocity; determining the equivalent sound source inside the structure and the source intensity of the equivalent sound source based on the vibration information includes:

[0023] Based on the vibration velocity of each structural node, determine the volume velocity corresponding to each first surface element;

[0024] Based on the volume velocity corresponding to each first surface unit, the equivalent sound source inside the structure and the source intensity of the equivalent sound source are determined.

[0025] In one possible implementation, determining the volume velocity corresponding to each first surface element based on the vibration velocity of each structural node includes:

[0026] The volume velocity of each first-face unit is determined by the vibration velocity of the structural nodes at the center of each first-face unit.

[0027] The beneficial effects of this invention are:

[0028] This invention, after obtaining the first pulsating pressure information on multiple second-face elements on the coupling surface of the fluid domain model through computational fluid dynamics simulation using a fluid domain model, maps the first pulsating pressure information on the second-face elements to structural nodes on the first-face elements of the structural model according to the correspondence of the surface elements, thus obtaining the second pulsating pressure information of each structural node. In this mapping process, only the consistency between the first and second-face elements needs to be ensured; the quality of the surface element partitioning is not critical. Then, the vibration information of each structural node is determined based on the second pulsating pressure information, and finally, the sound radiation is determined based on the vibration information. The sound radiation estimation process of this invention bypasses the step in the prior art of importing flow field information into structural elements as boundary conditions to obtain the vibration response, and since the quality of the surface element partitioning is not critical, it has higher computational efficiency, lower computational time cost, and better computational convergence. Attached Figure Description

[0029] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0030] Figure 1 A flowchart illustrating an embodiment of the fast acoustic radiation estimation method for fluid-domain structures provided by the present invention;

[0031] Figure 2 This invention provides a schematic diagram of a complex structure containing a fluid domain;

[0032] Figure 3 A schematic diagram of coupling surface division provided by the present invention;

[0033] Figure 4 A schematic diagram of RBF mapping provided by the present invention;

[0034] Figure 5 This invention provides a schematic diagram of the mean square vibration velocity of a structural node.

[0035] Figure 6 Provided by the present invention Figure 1 A flowchart of one embodiment of S104;

[0036] Figure 7 A schematic diagram of a wave superposition method provided by the present invention;

[0037] Figure 8 This invention provides a schematic diagram of a discrete wave superposition method.

[0038] Figure 9 A comparison diagram of boundary element method and wave superposition method provided by the present invention;

[0039] Figure 10 This is a flowchart illustrating another embodiment of the fast acoustic radiation estimation method for a fluid domain structure provided by the present invention. Detailed Implementation

[0040] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of them. All other embodiments obtained by those skilled in the art based on the embodiments of the present invention without creative effort are within the scope of protection of the present invention.

[0041] In the description of the embodiments of this invention, unless otherwise stated, "a plurality of" means two or more. The terms "first," "second," etc., used in the embodiments of this invention are used to distinguish similar objects, and are not used to describe a specific order or sequence, nor to indicate or imply their relative importance or implicitly specify the number of indicated technical features. It should be understood that such data can be interchanged where appropriate so that embodiments of this application can be implemented in orders other than those illustrated or described herein, and the objects distinguished by "first," "second," etc., are generally of the same class, and the number of objects is not limited; for example, a first object can be one or more.

[0042] In this document, the term "embodiment" means that a particular feature, structure, or characteristic described in connection with an embodiment may be included in at least one embodiment of the invention. The appearance of this phrase in various places throughout the specification does not necessarily refer to the same embodiment, nor is it a separate or alternative embodiment mutually exclusive with other embodiments. It will be explicitly and implicitly understood by those skilled in the art that the embodiments described herein can be combined with other embodiments.

[0043] Reference Figure 1 The diagram illustrates a flowchart of an embodiment of the fast acoustic radiation estimation method for fluid-domain structures provided by the present invention. The method includes:

[0044] S101, Establish the structural model; the structural model includes the coupling surface between the structure and the fluid domain, the coupling surface includes multiple first surface elements, and the first surface elements include structural nodes.

[0045] This embodiment is applicable to acoustic radiation estimation of simple structures containing fluid domains, as well as complex structures containing fluid domains. Complex structures can refer to structures with large coupling surface areas or structures containing multiple coupling surfaces. It is also applicable to free-field acoustic radiation prediction when complex structures contain fluid and are simultaneously subjected to excitation from multiple devices.

[0046] The structural model of the structure can be created using Finite Element Method (FEM) software. The structural model includes a coupling surface, which contains multiple first-face elements, and each first-face element contains at least one structural node.

[0047] The coupling surface refers to the interface between the fluid domain and the structure. It is a key boundary region for the transfer of physical quantities between the fluid domain and the structural domain. Its function is to enable bidirectional or unidirectional data interaction between fluid loads (such as pressure and shear force) and structural responses (such as displacement and vibration).

[0048] S102, Establish a fluid domain model, which includes a coupling surface; the coupling surface of the fluid domain model includes multiple second surface elements that correspond one-to-one with multiple first surface elements.

[0049] A fluid domain model can be established using computational fluid dynamics (CFD) simulation software. ANSYS Fluent can be used for CFD simulation. The fluid domain model can be the SAS turbulence model, an improved version of the URANS model. This model not only accurately obtains sound source information in the flow field but also has low dependence on the computational grid, improving the efficiency of fluid simulation. Its governing equations are as follows:

[0050] (1)

[0051] (2)

[0052] in:

[0053]

[0054] The SAS correction term is used to adjust the turbulence model to accommodate different flow characteristics, where: , , .

[0055] S103, based on the fluid domain model, performs fluid dynamics simulation to obtain the first pulsating pressure information on multiple second surface elements.

[0056] After the fluid domain model is established, boundary conditions can be set, CFD simulation can be performed, and the first pulsating pressure information on the second surface element can be obtained.

[0057] In one example, the boundary condition could be: the inlet flow rate is set to 50 rpm, which is the rated flow rate of the centrifugal pump. m 3 / h The outlet reference pressure was set to 1 MPa to simulate the pressure environment of a spacecraft compartment at a diving depth of 100 m, and the centrifugal pump speed was set to 1450 rpm. The calculation step size and number of steps were set. The unsteady flow field calculation settings were modified from the steady flow field settings, changing the solution type from steady-state to transient. The total calculation time was set to 0.331 s, corresponding to 8 impeller rotations, and the time step was set to... s, meaning one solution is performed every 3° rotation of the impeller, resulting in a total of 960 time-step files. A high-resolution solution scheme is employed, with a convergence accuracy set to [value missing]. .

[0058] S104, the first pulsating pressure information on the second surface unit is mapped to the structural nodes on the corresponding first surface unit through a neural network model, thereby obtaining the second pulsating pressure information of each structural node.

[0059] The first pulsating pressure information can be mapped onto the structural nodes using a radial basis function (RBF) neural network.

[0060] S105, based on the second pulsating pressure information, determine the vibration information of each structural node, and determine the acoustic radiation of the structure based on the vibration information.

[0061] The second pulsating pressure information of each structural node can be input into the ANSYS software to obtain the vibration information of each structural node, which may include vibration velocity. Then, based on the vibration information of each structural node, the acoustic radiation at any field point near the structure can be determined.

[0062] The acoustic radiation estimation method for structures provided in this embodiment can be applied to acoustic radiation estimation systems for structures. These systems can be software systems running on terminal devices. Terminal devices can be tablets, augmented reality (AR) / virtual reality (VR) devices, laptops, ultra-mobile personal computers (UMPCs), netbooks, personal digital assistants (PDAs), mobile phones, etc. This embodiment does not impose any restrictions on the specific type of terminal device.

[0063] In summary, this embodiment, after obtaining the first pulsating pressure information on multiple second-face elements on the coupling surface of the fluid domain model through computational fluid dynamics simulation using a fluid domain model, maps the first pulsating pressure information on the second-face elements to structural nodes on the first-face elements of the structural model according to the correspondence of the surface elements, thereby obtaining the second pulsating pressure information of each structural node. In this mapping process, only the consistency between the first and second-face elements needs to be ensured; the quality of the surface element partitioning is not critical. Then, the vibration information of each structural node is determined based on the second pulsating pressure information, and finally, the sound radiation is determined based on the vibration information. The sound radiation estimation process in this embodiment bypasses the step in the prior art of importing flow field information into structural elements as boundary conditions to obtain the vibration response, and the quality of the surface element partitioning is not critical. Therefore, it has higher computational efficiency, lower computational time cost, and better computational convergence.

[0064] In some embodiments of the present invention, the coupling surface can be customized during preprocessing, as shown in the following reference. Figure 2 The diagram shows a structural schematic of a complex structure containing a fluid domain provided by the present invention. (Refer to...) Figure 3 This illustrates a schematic diagram of the coupling surface division provided by the present invention. Figure 2 On the complex structure shown, the coupling surface is custom-divided into multiple segments, such as coupling surface 1, coupling surface 2, and coupling surface 3. Therefore, one of the coupling surfaces, for example coupling surface 3, can be selected using the ACT plugin in ANSYS software. Then, the information of the first element on that coupling surface can be determined. This first element information includes the coordinates of the center node of the first element, the element number, and the time-domain signal of the pulsating pressure on the first element. Based on the extracted first element information, the acoustic radiation at the field points near the coupling surface is estimated.

[0065] In some embodiments of the present invention, the first pulsating pressure information includes a target pulsating pressure frequency domain signal. S103 includes: performing computational fluid dynamics simulation based on a fluid domain model to obtain pulsating pressure time domain signals on multiple second surface units; converting the pulsating pressure time domain signals into pulsating pressure frequency domain signals, and extracting signals within a preset frequency range from the pulsating pressure frequency domain signals as the target pulsating pressure frequency domain signal.

[0066] The time-domain signal of pulsating pressure can be converted into the frequency-domain signal of pulsating pressure using FFT transformation in MATLAB software. Then, the main frequency band signal, that is, the signal in the frequency range of 2-300 Hz, can be extracted.

[0067] In some embodiments of the present invention, when the first pulsating pressure information is a target pulsating pressure frequency domain signal, the mapping can be achieved through an RBF neural network. The RBF neural network mapping process is as follows: Figure 4 As shown.

[0068] In this embodiment, the first pulsating pressure information may include concentrated force. The specific process of the RBF neural network to implement the mapping is as follows: the target pulsating pressure frequency domain signal on the second facet unit is mapped to the corresponding first facet unit; the total concentrated force of the first facet unit is determined according to the target pulsating pressure frequency domain signal on the first facet unit; the total concentrated force of the first facet unit is evenly distributed to the structural nodes on the first facet unit to obtain the concentrated force of each structural node.

[0069] RBF neural networks obtain the concentrated force of structural nodes in the frequency domain through interpolation. The RBF interpolation function is:

[0070] (3)

[0071] in: The radial basis function represents the distance between the interpolation point and the target point. The interpolation coefficients are represented by the Wendland function, which characterizes the basis functions.

[0072] (4)

[0073] In some embodiments of the present invention, the above mapping process can be implemented based on the number and coordinates of the first face unit, the number and coordinates of the second face unit, and the number and coordinates of the structural nodes.

[0074] In some embodiments of the present invention, after mapping is completed, a macro file suitable for ANSYS software can be generated to facilitate subsequent calculations.

[0075] In some embodiments of the present invention, macro file commands can be entered in the ANSYS software interface to apply the mapped concentrated force to the structural nodes, set the fixed-support constraints on both sides, and calculate the harmonic response based on the complete method. ANSYS command flow code is written to directly extract the second-face element number, structural node number, structural node coordinates, and vibration velocity of the structural nodes from the structural model, and output as a MATLAB-compatible .mat file.

[0076] Vibration velocity can be the mean square vibration velocity, as referenced. Figure 5 The diagram shows a schematic of the mean square vibration velocity of a structural node provided by the present invention.

[0077] The wave superposition method, because it does not require boundary integral calculations, has advantages over the boundary element method in terms of higher computational efficiency and stronger stability. Its physical meaning is: it is assumed that the sound radiation of any complex radiator can be superimposed from the sound radiation of a series of equivalent sound sources located inside the radiator, and the source intensity of the equivalent sound sources can be obtained by matching the velocity boundary conditions of the radiator surface. Therefore, in some embodiments of the present invention, referring to... Figure 6 The steps for determining the acoustic radiation of a structure based on vibration information may include:

[0078] S601, based on vibration information, determine the equivalent sound source inside the structure and the source strength of the equivalent sound source.

[0079] S602, based on the equivalent sound source and source strength, the superimposed sound radiation of the equivalent sound source at the target field point is determined by combining the wave superposition method.

[0080] In some embodiments of the present invention, S601 may include: determining the volume velocity corresponding to each first surface unit based on the vibration velocity of each structural node; and determining the equivalent sound source inside the structure and the source intensity of the equivalent sound source based on the volume velocity corresponding to each first surface unit.

[0081] Volumetric velocity refers to the alternating flux per unit area per time produced by a wave on a specified surface, which is equal to the product of the particle velocity and the area of ​​that surface.

[0082] In one example, the calculation process for volumetric velocity can be as follows: determine the normal vector of the first facet element, then determine the projection of the vibration velocity of each structural node on the first facet element onto the normal vector, and finally integrate the projections of the vibration velocities of all structural nodes on the first facet element onto the entire first facet element to obtain the "volume throughput" of the first facet element to the fluid medium, i.e., the volumetric velocity. Alternatively, the volumetric velocity can be obtained by multiplying the projection of the vibration velocity of the structural node located at the center of the first facet element onto the normal vector with the area of ​​the first facet element.

[0083] Reference Figure 7This diagram illustrates a wave superposition method provided by the present invention. A series of simple sound sources located inside the radiator S (structure) are positioned at a certain point in space. r The radiated sound pressure at a certain point can be expressed as:

[0084] (5)

[0085] In the formula: i represents the imaginary part, The density of the fluid medium; It is the angular frequency of simple harmonic motion; This represents the source strength value of a simple source inside the radiator; The Green's function is defined as follows:

[0086] (6)

[0087] In the formula: Let be the distance between the field point and the simple source, and satisfy:

[0088] (7)

[0089] In the formula: Wave number; The function. Equation (5) is the integral equation of the wave superposition method theory used to solve the sound radiation problem. The linear Euler equation describing the relationship between particle velocity and sound pressure is: The field point can be obtained from equation (5). r The velocity of the particle at that point is:

[0090] (8)

[0091] Therefore, the normal velocity of the particle on the boundary surface S of the radiating body can be obtained as:

[0092] (9)

[0093] In the formula: v(r S ) This is the normal velocity of the particle, which is also the volume velocity of the first surface element. It is the coordinate vector of the particles on the boundary surface of the radiator; To calculate the normal gradient of particles on the boundary surface. Theoretically, a simple source can be placed anywhere inside the radiator, but for ease of calculation, it is generally placed inside the radiator with a thickness of... On the spherical shell. Therefore, equation (9) becomes

[0094] (10)

[0095] In the formula: The surface of the sphere where the simple source is located; for Surface particle coordinate vector. Because... Located inside the radiator, Always less than Therefore, equation (9) avoids the singularity problem of boundary integrals in the boundary element method. (Refer to...) Figure 8 The diagram illustrates a discrete wave superposition method provided by the present invention. Now, The surface is discretized into N segments, each segment is... If expressed in this way, equation (10) can be discretized as follows:

[0096] (11)

[0097] If each segment in equation (9) Small enough to display each segment The source of strength If we approximate it to be a constant, then equation (9) can be expressed as:

[0098] (12)

[0099] In the formula: For each segment The source strength. When the velocity boundary conditions of N points uniformly distributed on the surface of the radiator are known, the source strength of each simple source inside the radiator can be obtained by performing a matrix transformation on equation (11). Similarly, equation (5) is discretized, as shown. Figure 9 As shown, the sound pressure solution at any spatial point can be obtained as follows:

[0100] (13)

[0101] In structural acoustic radiation calculation, after obtaining the vibration velocity of the surface of the radiating body structure through the finite element method, the source intensity of the simple source inside the radiating body can be obtained according to equation (12). Then, the radiated sound field solution of the radiator can be obtained according to equation (13). Therefore, after obtaining the vibration velocity of the surface nodes of the structure by the finite element method, the position of the center node, the volume velocity and its normal vector of each element are calculated according to the extracted element and node information. In order to overcome the influence of the characteristic frequency, the wave superposition method changes the Green's function to a combination of monopole and dipole as shown in equation (14), so as to more accurately describe the propagation characteristics of the sound field and realize the prediction of sound radiation.

[0102] (14)

[0103] Reference Figure 9The diagram shows a comparison between the boundary element method and wave superposition method provided by this invention. By comparing this invention with traditional numerical methods, it can be seen that within the main characteristic frequency range, the method proposed in this invention exhibits good agreement and consistency with traditional numerical methods.

[0104] Reference Figure 10 This diagram illustrates a flowchart of another embodiment of the rapid acoustic radiation estimation method for a fluid-domain structure provided by the present invention. A fluid domain model is established and boundary conditions are set using ANSYS Fluent software. CFD simulation is then performed to obtain the time-domain signal of the pulsating pressure of each second-face element on the fluid domain coupling surface. Simultaneously, a finite element structural model is established using ANSYS software, and the structural domain coupling surface on the structural model is determined. The pulsating pressure time-domain signal is then frequency-domain transformed and truncated to obtain the target pulsating pressure frequency-domain signal of each second-face element. An RBF neural network is used to map the target pulsating pressure frequency-domain signal of the second-face element to the structural nodes on the structural domain coupling surface, obtaining the concentrated force corresponding to each structural node. Finally, the vibration velocity of each structural node is calculated based on the concentrated force corresponding to each structural node. The vibration source and source strength are determined based on the vibration velocity, and the acoustic radiation caused by each vibration source at the field point is determined using the wave superposition method.

[0105] Compared to traditional numerical methods, this invention can more efficiently capture the coupling effect between fluid and structure, thereby significantly reducing computation time and resource consumption. Traditional numerical methods are limited by computational resources and have high mesh requirements, while this invention has low requirements for both mesh and computational resources. Table 1 below provides a comparison of computation time and synthesized sound power level between this invention and traditional numerical methods.

[0106] Table 1. Comparison of computation time and synthesized sound power level between the present invention and traditional numerical methods.

[0107]

[0108] As shown in Table 1, this invention offers faster calculations while maintaining high accuracy. This invention has significant advantages in practical engineering applications, particularly for analyzing the acoustic radiation characteristics of structures containing fluid domains.

[0109] Those skilled in the art will understand that all or part of the processes of the methods described in the above embodiments can be implemented by a computer program instructing related hardware, and the program can be stored in a computer-readable storage medium. The computer-readable storage medium may be a disk, optical disk, read-only memory, or random access memory, etc.

[0110] The above description is only a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any changes or substitutions that can be easily conceived by those skilled in the art within the scope of the technology disclosed in the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for rapid estimation of acoustic radiation in structures containing fluid domains, characterized in that, The method includes: Establish a structural model; the structural model includes a coupling surface between the structure and the fluid domain, the coupling surface includes multiple first surface elements, and the first surface elements include structural nodes; A fluid domain model is established; the fluid domain model includes the coupling surface; the coupling surface of the fluid domain model includes multiple second surface elements that correspond one-to-one with the multiple first surface elements; Fluid dynamics simulation is performed based on the fluid domain model to obtain the first pulsating pressure information on the plurality of second surface units; The first pulsating pressure information on the second surface unit is mapped to the structural node on the first surface unit corresponding to the second surface unit through a neural network model, thereby obtaining the second pulsating pressure information of each structural node; Based on the second pulsating pressure information, the vibration information of each structural node is determined, and the acoustic radiation of the structure is determined based on the vibration information; The first pulsating pressure information includes the target pulsating pressure frequency domain signal, and the neural network model mapping process includes: The target pulsating pressure frequency domain signal on the second face unit is mapped to the corresponding first face unit; the total concentrated force of the first face unit is determined based on the target pulsating pressure frequency domain signal on the first face unit; the total concentrated force of the first face unit is evenly distributed to the structural nodes on the first face unit to obtain the concentrated force of each structural node.

2. The method for rapid estimation of acoustic radiation in a fluid-domain structure according to claim 1, characterized in that, The computational fluid dynamics simulation based on the fluid domain model, to obtain the first pulsating pressure information on the plurality of second surface elements, includes: Computational fluid dynamics simulation is performed based on the fluid domain model to obtain the time-domain signal of the pulsating pressure on the plurality of second surface units; The pulsating pressure time-domain signal is converted into a pulsating pressure frequency-domain signal, and a signal within a preset frequency range is extracted from the pulsating pressure frequency-domain signal as the target pulsating pressure frequency-domain signal.

3. The method for rapid estimation of acoustic radiation in structures containing fluid domains according to claim 2, characterized in that, The preset frequency range is 2-300Hz.

4. The method for rapid estimation of acoustic radiation in structures containing fluid domains according to claim 1, characterized in that, The structural model was constructed using finite element analysis software, and the fluid domain model was constructed using fluid dynamics simulation software.

5. The method for rapid estimation of acoustic radiation in a fluid-domain structure according to claim 1, characterized in that, The neural network model is a neural network model based on radial basis functions.

6. The method for rapid estimation of acoustic radiation in a fluid-domain structure according to claim 1, characterized in that, The step of determining the vibration information of each structural node based on the second pulsating pressure information includes: The second pulsating pressure information of each structural node is input into the finite element analysis software to obtain the vibration information of each structural node.

7. The method for rapid estimation of acoustic radiation in a fluid-domain structure according to claim 1, characterized in that, Determining the acoustic radiation of the structure based on the vibration information includes: Based on the vibration information, determine the equivalent sound source inside the structure and the source strength of the equivalent sound source; Based on the equivalent sound source and the source strength, the superimposed sound radiation of the equivalent sound source at the target field point is determined by combining the wave superposition method.

8. The method for rapid estimation of acoustic radiation in a fluid-domain structure according to claim 7, characterized in that, The vibration information includes vibration velocity; determining the equivalent sound source inside the structure and the source intensity of the equivalent sound source based on the vibration information includes: Based on the vibration velocity of each structural node, determine the volume velocity corresponding to each first surface element; Based on the volume velocity corresponding to each first surface unit, the equivalent sound source inside the structure and the source intensity of the equivalent sound source are determined.

9. The method for rapid estimation of acoustic radiation in a fluid-domain structure according to claim 8, characterized in that, The step of determining the volume velocity corresponding to each first surface element based on the vibration velocity of each structural node includes: The volume velocity of each first-face unit is determined by the vibration velocity of the structural nodes at the center of each first-face unit.

Citation Information

Patent Citations

  • Pipeline structure vibration sound radiation analogue simulation method

    CN116341154A

  • Wind power boosting rotor external field vibration sound radiation simulation method

    CN119670619A