Nonlinear structural sound intensity visualization method, device, equipment and readable storage medium

Abstract

The present application relates to a method, device, equipment and readable storage medium for visualizing sound intensity of nonlinear structures, and relates to the field of optical mechanical simulation technology, including inputting a nonlinear structure into a transient dynamics simulation model to obtain structural sound intensity simulation data of the nonlinear structure, the structural sound intensity simulation data including simulation result data of each entity unit in the nonlinear structure, the simulation result data including displacement response and stress response; performing vector processing on the simulation result data of each entity unit to obtain structural sound intensity vector field data of each entity unit; drawing a sound intensity vector image of the nonlinear structure based on the structural sound intensity vector field data of each entity unit; and importing the sound intensity vector image into a visualization interactive interface. The present application can realize the visualization of nonlinear structural sound intensity, and effectively improve the efficiency and accuracy of directional optimization of nonlinear structural dynamic characteristics.

Description

Technical Field

[0001] The present application relates to the field of optical mechanical simulation technology, and in particular to a method, device, equipment and readable storage medium for visualizing sound intensity of nonlinear structures. Background Art

[0002] With the development of modern precision manufacturing industry, vibration and noise prediction and control is a very necessary research topic, which involves the reliability of the function and performance of precision equipment during operation. On the one hand, vibration will cause the function and performance of precision equipment to fail to reach the expected design level; on the other hand, local vibration can easily cause fatigue damage to the structure of precision equipment, greatly reducing the service life of the equipment. As a typical precision equipment, the performance of nonlinear large-aperture optical machine is extremely sensitive to vibration, so the study of the vibration characteristics of nonlinear large-aperture optical machine structure is a very important task.

[0003] However, traditional methods lack precision and accuracy in evaluating the vibration performance of nonlinear structures, and are unable to accurately optimize the dynamic characteristics of nonlinear structures. Large-scale trial changes combined with iterative calculations are often required to gradually find the direction of structural optimization, resulting in low efficiency and poor accuracy. Summary of the invention

[0004] The present application provides a method, device, equipment and readable storage medium for visualizing sound intensity of a nonlinear structure, so as to solve the problems of low efficiency and poor precision in directional optimization of dynamic characteristics of a nonlinear structure in the related art.

[0005] In a first aspect, a nonlinear structural sound intensity visualization method is provided, comprising the following steps:

[0006] Inputting a nonlinear structure into a transient dynamics simulation model to obtain structural sound intensity simulation data of the nonlinear structure, wherein the structural sound intensity simulation data includes simulation result data of each entity unit in the nonlinear structure, and the simulation result data includes displacement response and stress response;

[0007] Perform vector processing on the simulation result data of each entity unit to obtain the structural sound intensity vector field data of each entity unit;

[0008] Drawing the sound intensity vector image of the nonlinear structure based on the structural sound intensity vector field data of each entity unit;

[0009] The sound intensity vector image is imported into a visual interactive interface.

[0010] In some embodiments, the transient dynamics simulation model includes Rayleigh damping.

[0011] In some embodiments, the vector processing of the simulation result data of each entity unit to obtain the structural sound intensity vector field data of each entity unit includes:

[0012] Performing time-frequency transformation on the stress response and displacement response of the first entity unit in the time domain respectively, to obtain the displacement response and stress response of the first entity unit corresponding to each frequency in the frequency domain;

[0013] A typical frequency is selected from multiple frequencies based on the frequency selection mechanism and the displacement response and stress response corresponding to each frequency;

[0014] Structural acoustic intensity vector field data of the first solid element is generated based on the displacement response and the stress response corresponding to the typical frequency.

[0015] In some embodiments, generating the structural sound intensity vector field data of the first entity unit based on the displacement response and the stress response corresponding to the typical frequency includes:

[0016] The displacement response and stress response corresponding to the typical frequency are calculated based on the structural sound intensity expression in the Cartesian coordinate system to obtain the structural sound intensity vector field data of the first entity unit.

[0017] In some embodiments, the structure-borne sound intensity expression in the Cartesian coordinate system is:

[0018]

[0019]

[0020]

[0021] Among them, I x ,I y and I z Respectively represent the structural sound intensity of the solid unit in the x-direction, y-direction and z-direction, σ x , σ y and σ z They represent the surface normal stress of the solid element in the x-direction, y-direction and z-direction respectively. and Respectively represent the conjugate displacement of the solid element in the x-direction, y-direction and z-direction, τ xy and τ yx represents the shear stress of the solid element on the xy plane, τ xz and τ zx represents the shear stress of the solid element on the xz plane, τ yz and τ zy It represents the shear stress of the solid element on the yz plane, and ω is the circular frequency of vibration.

[0022] In some embodiments, the step of drawing the sound intensity vector image of the nonlinear structure based on the structural sound intensity vector field data of each entity unit includes:

[0023] Based on the mapping function, the structural sound intensity vector field data of each entity unit is uniformly distributed to obtain the structural sound intensity vector field data after mapping changes corresponding to each entity unit;

[0024] The sound intensity vector image of the nonlinear structure is drawn according to the structural sound intensity vector field data after mapping changes corresponding to each entity unit.

[0025] In some embodiments, the expression of the mapping function is:

[0026]

[0027] Among them, I n represents the structural sound intensity vector field data after mapping changes, I represents the structural sound intensity vector field data before mapping changes, and I max Indicates the maximum value of the structure-borne sound intensity, I min Indicates the minimum value of the structure-borne sound intensity amplitude.

[0028] In a second aspect, a nonlinear structural sound intensity visualization device is provided, comprising:

[0029] A simulation unit, which is used to input a nonlinear structure into a transient dynamics simulation model to obtain structural sound intensity simulation data of the nonlinear structure, wherein the structural sound intensity simulation data includes simulation result data of each entity unit in the nonlinear structure, and the simulation result data includes displacement response and stress response;

[0030] A processing unit, which is used to perform vector processing on the simulation result data of each entity unit to obtain the structural sound intensity vector field data of each entity unit;

[0031] A drawing unit, which is used to draw an acoustic intensity vector image of the nonlinear structure based on the structural acoustic intensity vector field data of each entity unit;

[0032] A visualization unit is used to import the sound intensity vector image into a visualization interaction interface.

[0033] In a third aspect, a nonlinear structural sound intensity visualization device is provided, comprising: a memory and a processor, wherein the memory stores at least one instruction, and the at least one instruction is loaded and executed by the processor to implement the aforementioned nonlinear structural sound intensity visualization method.

[0034] In a fourth aspect, a computer-readable storage medium is provided, wherein the computer storage medium stores a computer program, and when the computer program is executed by a processor, the aforementioned nonlinear structural sound intensity visualization method is implemented.

[0035] The beneficial effects brought about by the technical solution provided in this application include: realizing the visualization of the sound intensity of nonlinear structures, and effectively improving the efficiency and accuracy of the directional optimization of the dynamic characteristics of nonlinear structures.

[0036] The present application provides a method, device, equipment and readable storage medium for visualizing sound intensity of nonlinear structures, including inputting a nonlinear structure into a transient dynamics simulation model to obtain structural sound intensity simulation data of the nonlinear structure, the structural sound intensity simulation data including simulation result data of each entity unit in the nonlinear structure, the simulation result data including displacement response and stress response; vector processing the simulation result data of each entity unit to obtain structural sound intensity vector field data of each entity unit; drawing the sound intensity vector image of the nonlinear structure based on the structural sound intensity vector field data of each entity unit; and importing the sound intensity vector image into a visualization interactive interface. The present application performs post-processing on the stress and displacement obtained by finite element simulation to synthesize the structural sound intensity vector field data of each entity unit in the nonlinear structure, and draws and visualizes the sound intensity vector image, so as to accurately identify the vibration noise source in the nonlinear structure according to the visualized sound intensity vector image, and analyze the main propagation path of the vibration noise, thereby realizing the precise directional optimization of the dynamic characteristics of the nonlinear structure, which does not require large-scale trial changes and iterative calculations, and effectively improves the efficiency and accuracy of the directional optimization of the dynamic characteristics of the nonlinear structure. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. 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 creative work.

[0038] Figure 1 A schematic diagram of a flow chart of a nonlinear structural sound intensity visualization method provided in an embodiment of the present application;

[0039] Figure 2 A schematic diagram of the structure of the optical-mechanical support structure of a large-aperture curved reflector provided in the embodiment of the present application, "Scheme 1";

[0040] Figure 3 A schematic diagram of the structure of the optical-mechanical support structure of the large-aperture curved reflector provided in the embodiment of the present application, "Scheme 2";

[0041] Figure 4 The sound intensity cloud diagram of the 7Hz vibration structure of the large-aperture curved reflective mirror in the first embodiment of the present application;

[0042] Figure 5 A schematic diagram of the structure of a nonlinear structural sound intensity visualization device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0043] In order to make the purpose, technical solution and advantages of the embodiments of the present application clearer, the technical solution in the embodiments of the present application will be clearly and completely described below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of this application.

[0044] The embodiments of the present application provide a method, device, equipment and readable storage medium for visualizing sound intensity of a nonlinear structure, which can solve the problems of low efficiency and poor precision in directional optimization of dynamic characteristics of nonlinear structures in related technologies.

[0045] Figure 1 A nonlinear structural sound intensity visualization method provided in an embodiment of the present application includes the following steps:

[0046] Step S10: inputting a nonlinear structure into a transient dynamics simulation model to obtain structural sound intensity simulation data of the nonlinear structure, wherein the structural sound intensity simulation data includes simulation result data of each entity unit in the nonlinear structure, and the simulation result data includes displacement response and stress response;

[0047] Furthermore, the transient dynamics simulation model includes Rayleigh damping.

[0048] Demonstratively, the structural sound intensity theory is based on the basic theory of continuous medium mechanics and is transplanted from the concept of sound intensity in acoustics. The structural sound intensity method focuses on the stress inside the structure and the movement of the medium, and takes into account the stress and displacement of the structure during the vibration process. Structural sound intensity is defined as the vibration power flow per unit width in a specified direction, also known as power flow density, which is a vector representing the magnitude and direction of energy flow through a certain point on the structure. It can be seen that the structural sound intensity method can be used to solve the problem of structural vibration noise. It is a power flow method that studies the transmission of structural vibration energy from the perspective of fluctuations. It can be used as an effective tool for diagnosing structural defects, identifying the source of vibration noise in the structure, and analyzing the main propagation path of structural vibration noise. Therefore, in this embodiment, the analysis of the vibration characteristics of nonlinear large-aperture optomechanical structures is realized based on the structural sound intensity method.

[0049] This embodiment first performs transient dynamic simulation analysis of the entity units in the nonlinear large-aperture optomechanical structure based on the transient dynamic simulation model in the finite element simulation analysis software, and obtains the structural sound intensity simulation data of the nonlinear large-aperture optomechanical structure. The structural sound intensity simulation data includes simulation result data of each entity unit in the nonlinear large-aperture optomechanical structure, and the simulation result data includes displacement response and stress response.

[0050] Since the energy represented by the structural sound intensity belongs to the energy transmitted by the real work part in the structure, that is, the energy diffused from the vibration (noise) source to the surrounding through the structure, it is not the energy converted without loss between the potential energy and kinetic energy of the solid unit structure itself. Therefore, in order to analyze the energy dissipated by the nonlinear large-aperture optomechanical structure during the vibration process, it is necessary to set the damping of the structure in the dynamic simulation analysis. Since nonlinear transient finite element analysis is considered, the calculated Rayleigh damping is added in the transient simulation analysis to calculate the real work part transmitted by the structure during the vibration process.

[0051] In the analytical calculation, Rayleigh damping is actually solved by dividing the damping matrix in the dynamic equation into the mass matrix and the stiffness matrix. The relationship is shown in equation (1).

[0052] [C]=α[M]+β[K] (1)

[0053] Where [C] is the damping matrix, [M] is the mass matrix, [K] is the stiffness matrix, α is the mass proportional coefficient, and β is the stiffness proportional coefficient.

[0054] It can be seen that when setting the Rayleigh damping, it is particularly important to select the mass proportional coefficient α and the stiffness proportional coefficient β, and the structural damping and these two proportional coefficients also satisfy the relationship of formula (2).

[0055]

[0056] Where c is the structural damping and ω is the circular frequency of vibration.

[0057] Considering that α is a coefficient proportional to the mass matrix, when the structure being analyzed is non-homogeneous or the material density does not change slowly with space, it is not suitable to use the mass proportional coefficient α. At this time, α should be set to zero, otherwise it will affect the solution accuracy; and β is a coefficient proportional to the stiffness matrix. Therefore, when studying the nonlinear dynamic response of a structure with a sudden change in stiffness, it is not suitable to use the stiffness proportional coefficient β. At this time, β should be set to zero, otherwise it will affect the solution accuracy. Since this application mainly involves the nonlinear dynamic response analysis of (homogeneous) structures, the coefficient β is set to zero and only the coefficient α is used.

[0058] Among them, this embodiment also performs secondary development of the command stream of the finite element simulation software, and can convert the transient simulation analysis results in binary form into structural sound intensity simulation data (including displacement, stress, position and other result data of each entity unit), and store it in the database, waiting for the call of the external program. Specifically, firstly, by writing the secondary development program of the command stream of the finite element simulation software, the total number of entity units of the nonlinear large-aperture optical-mechanical structure, the center coordinates of the entity units, and the stresses in six directions and displacements in three directions of each entity unit under each transient load step are derived from the transient simulation analysis result file; the entity unit in this embodiment is preferably a high-order entity unit with 20 nodes, which is a hexahedral unit with 8 corner nodes and 12 edge nodes (it can also be reduced to a tetrahedral unit), and then the displacement of each node of the derived entity unit is calculated by arithmetic averaging, as the displacement response of the entity unit (the displacement response is the displacement response in the time domain at this time) and stored in the database, as shown in formula (3).

[0059]

[0060] Where U i represents the displacement of the node, Represents the displacement of the solid element.

[0061] Step S20: performing vector processing on the simulation result data of each entity unit to obtain the structural sound intensity vector field data of each entity unit;

[0062] Furthermore, the simulation result data of each entity unit is vector-processed to obtain the structural sound intensity vector field data of each entity unit, including:

[0063] Performing time-frequency transformation on the stress response and displacement response of the first entity unit in the time domain respectively, to obtain the displacement response and stress response of the first entity unit corresponding to each frequency in the frequency domain;

[0064] A typical frequency is selected from multiple frequencies based on the frequency selection mechanism and the displacement response and stress response corresponding to each frequency;

[0065] Structural acoustic intensity vector field data of the first solid element is generated based on the displacement response and the stress response corresponding to the typical frequency.

[0066] Furthermore, the generating of the structural sound intensity vector field data of the first entity unit based on the displacement response and the stress response corresponding to the typical frequency includes:

[0067] The displacement response and stress response corresponding to the typical frequency are calculated based on the structural sound intensity expression in the Cartesian coordinate system to obtain the structural sound intensity vector field data of the first entity unit.

[0068] Furthermore, the structure-borne sound intensity expression in the Cartesian coordinate system is:

[0069]

[0070]

[0071]

[0072] Among them, I x ,I y and I z Respectively represent the structural sound intensity of the solid unit in the x-direction, y-direction and z-direction, σ x , σ y and σ z They represent the surface normal stress of the solid element in the x-direction, y-direction and z-direction respectively. and Respectively represent the conjugate displacement of the solid element in the x-direction, y-direction and z-direction, τ xy and τ yx represents the shear stress of the solid element on the xy plane, τ xz and τ zx represents the shear stress of the solid element on the xz plane, τ yz and τ zy It represents the shear stress of the solid element on the yz plane, and ω is the circular frequency of vibration.

[0073] Exemplarily, this embodiment reads the simulation result data of each entity unit in the database, and synthesizes the structural sound intensity vector field data according to the theoretical calculation formula of the structural sound intensity of the entity unit. Among them, the structural sound intensity SI (Structural Intensity) of the entity unit can be expressed as:

[0074]

[0075] In the formula, and are the surface normal, 1-direction and 2-direction displacement conjugates of the solid element respectively; σ n is the surface normal stress of the solid element, τ n1 and τ n2 are the shear stresses on the surface of the solid element in directions 1 and 2 respectively. It represents the structural sound intensity component contributed by normal stress and displacement in the corresponding direction to the structure. and It represents the structural sound intensity component contributed by shear stress and displacement in the corresponding direction to the structure. Therefore, in the Cartesian rectangular coordinate system, the structural sound intensity of the solid element in the x, y and z directions can be expressed as:

[0076]

[0077]

[0078]

[0079] Since the coordinate system direction of the entity element itself is parallel to the overall coordinate system by default, there is no process of transforming data between two different coordinate systems.

[0080] Since the structural sound intensity field is a vector field in the frequency domain, it is necessary to perform time-frequency transformation on the stress response and displacement response in the time domain to obtain the stress and displacement response data of the entity unit in the frequency domain. In this embodiment, the fast Fourier transform algorithm is used as the time-frequency transformation method, and the transformation process can be expressed by the following formula (8):

[0081] R f =FFT<R t > (8)

[0082] In the formula, R t Represents the stress response and displacement response in the time domain, R f represents the stress response and displacement response in the frequency domain and R f Is plural.

[0083] Thus, the stress response and displacement response of the entity unit at each frequency in the frequency domain can be obtained. However, what needs to be calculated in this embodiment is the structural sound intensity vector field data of the entity unit at a specific frequency, and the frequency must be representative, and the structural sound intensity vector field at this frequency can reflect the typical structural vibration mode and energy transfer information. Therefore, in order to obtain the representative characteristic frequency, this embodiment establishes a frequency screening algorithm mechanism.

[0084] The specific mechanism of the frequency screening algorithm is as follows: First, R f Perform complex modulo operations to obtain the amplitudes of the stress response and displacement response of the entity unit at each frequency; then, loop through the array sequence of the stress response and displacement response at each frequency corresponding to each entity unit, and mark the frequency bits that meet the following conditions at the same time:

[0085] ① The stress response amplitude corresponding to this frequency is greater than the stress response amplitude corresponding to the adjacent frequencies on both sides, or the displacement response amplitude corresponding to this frequency is greater than the displacement response amplitude corresponding to the adjacent frequencies on both sides;

[0086] ② The stress response amplitude corresponding to the frequency is not less than 10% of the maximum stress response amplitude of all entity units at the frequency, or the displacement response amplitude corresponding to the frequency is not less than 10% of the maximum displacement response amplitude of all entity units at the frequency.

[0087] Finally, the number of times each frequency is marked is counted, and the two to three frequencies with the most marks are recorded as the winning frequencies (i.e., typical frequencies). The displacement responses and stress responses at these winning frequencies are the stress responses and displacement responses required for the structural sound intensity of the solid unit. Therefore, the displacement responses and stress responses at these winning frequencies are substituted into equations (5) to (7), and the structural sound intensity vector field data of the solid unit at the selected typical frequency can be calculated and synthesized.

[0088] Step S30: drawing an acoustic intensity vector image of the nonlinear structure based on the structural acoustic intensity vector field data of each entity unit;

[0089] Furthermore, the step of drawing the sound intensity vector image of the nonlinear structure based on the structural sound intensity vector field data of each entity unit includes:

[0090] Based on the mapping function, the structural sound intensity vector field data of each entity unit is uniformly distributed to obtain the structural sound intensity vector field data after mapping changes corresponding to each entity unit;

[0091] The sound intensity vector image of the nonlinear structure is drawn according to the structural sound intensity vector field data after mapping changes corresponding to each entity unit.

[0092] Furthermore, the expression of the mapping function is:

[0093]

[0094] Among them, I n represents the structural sound intensity vector field data after mapping changes, I represents the structural sound intensity vector field data before mapping changes, and I max Indicates the maximum value of the structure-borne sound intensity, I min Indicates the minimum value of the structure-borne sound intensity amplitude.

[0095] For example, in order to achieve the purpose of making the area with smaller structure-borne sound intensity vector amplitude more obvious and the area with larger structure-borne sound intensity vector amplitude more distinguishable, this embodiment performs uniform distribution processing on the amplitude of the structure-borne sound intensity vector field data of each entity unit based on the mapping function to obtain the structure-borne sound intensity vector field data after the mapping change corresponding to each entity unit. Since it is necessary to ensure that the monotonicity of the mapping function remains unchanged, and the referenced data mapping function should satisfy the unchanged structure-borne sound intensity vector amplitude interval before and after the mapping transformation, this embodiment adopts the mapping function shown in formula (9):

[0096]

[0097] Then, based on the structural sound intensity vector field data after mapping transformation of each entity unit, the sound intensity vector spectrum of the nonlinear large-aperture optomechanical structure can be drawn.

[0098] Step S40: importing the sound intensity vector image into a visual interactive interface.

[0099] Exemplarily, the sound intensity vector image obtained in step S30 is imported into the visual interactive interface, and the user can intuitively see the vibration noise source of the nonlinear large-aperture optomechanical structure based on the sound intensity vector image presented in the visual interactive interface, and analyze the main propagation path of the vibration energy flow, thereby providing a design basis and evaluation criteria for the optimization of vibration and noise reduction of the nonlinear large-aperture optomechanical structure.

[0100] It can be seen that the embodiment of the present application studies and analyzes the calculation and visualization method of the structural sound intensity of a nonlinear structure during non-steady-state motion. Specifically, the stress and displacement obtained by finite element simulation are post-processed to synthesize the structural sound intensity vector field data of each entity unit in the nonlinear structure, and the sound intensity vector image of the nonlinear structure is drawn and visualized, that is, the structural sound intensity vector field data of a complex structure at any position and posture in three-dimensional space is calculated, generated and visualized, and then the vibration noise source in the nonlinear structure can be accurately identified according to the visualized sound intensity vector image, and the main propagation path of the vibration noise can be analyzed, thereby realizing accurate directional optimization of the dynamic characteristics of the nonlinear structure, without the need for large-scale trial changes and iterative calculations, effectively improving the efficiency and accuracy of the directional optimization of the dynamic characteristics of the nonlinear structure.

[0101] Combine the following Figure 2 and Figure 3 The embodiments of the present application are further explained.

[0102] Figure 2 The following is a scheme 1: an open optical-mechanical support structure of a large-aperture curved reflector with a flexible hinge support. Figure 3 The second scheme is shown: an open optical-mechanical support structure for a large-aperture curved reflector with a flexible hinge support. Compared with the first scheme, the position of the flexible hinge support in the second scheme has changed.

[0103] This embodiment performs time-domain transient fluid-solid coupling simulation on the optical-mechanical support structure of the open large-aperture curved reflector with flexible hinge support in Scheme 1 and Scheme 2, respectively, simulates the inclined open large-aperture curved reflector subjected to a constant wind speed load of 15m / s, and analyzes the dynamic response of the damping structure under the action of wind load; secondly, by writing a secondary development program of the finite element simulation software command flow, the total number of mirror entity units, the center coordinates of the entity units, and the stresses in six directions and displacements in three directions of each entity unit under each transient load step are derived from the transient simulation analysis result file; then the simulation result data in the database is read to perform A fast Fourier transform is performed, and according to the frequency screening algorithm mechanism, 3Hz, 7Hz, and 10Hz are screened out as the typical vibration frequencies of the optomechanical mirror of Scheme 1, and 3Hz, 7Hz, 10Hz, and 23Hz are screened out as the typical vibration frequencies of the optomechanical mirror of Scheme 2. Subsequently, the synthetic structural sound intensity vector field data are calculated using the theoretical calculation formula for the sound intensity of the solid unit structure. The mirror structural sound intensity vector field data are then uniformly distributed to achieve a visualization effect in which the areas with smaller structural sound intensity vector amplitudes are more obvious and the areas with larger structural sound intensity vector amplitudes have a greater degree of amplitude distinction, and the sound intensity vector image is visualized.

[0104] Among them, according to the analysis of the visualized sound intensity vector image, the maximum structural sound intensity amplitude of the 3Hz vibration form of Scheme 1 is 0.0210W / m 2 The maximum structural sound intensity amplitude of 7Hz vibration form is 0.00242W / m 2 (See Figure 4 The maximum structure-borne sound intensity in the form of 10Hz vibration is 0.000214W / m 2 ; The maximum structural sound intensity amplitude of the 3Hz vibration form of Scheme 2 is 0.0602W / m 2 The maximum structural sound intensity amplitude of 7Hz vibration form is 0.00465W / m 2 The maximum structural sound intensity amplitude of 10Hz vibration form is 0.000420W / m 2 The maximum structural sound intensity amplitude of 23Hz vibration form is 0.000150W / m 2 .

[0105] In addition, at 3 Hz, the structural sound intensity streamlines of the flat plate are relatively smooth, with few sound intensity vortices. The low-frequency vibration energy can almost all flow to the damping of the overall structure, without concentrating the energy at the center of the mirror surface due to the vortex, causing non-designed vibration of the structure. At 7 Hz (see Figure 4As shown in the figure), at frequencies such as 10 Hz and 10 Hz, due to the reinforcement of the mirror back structure, the energy path from the vibration source to the damping node of the overall structure is interrupted to a certain extent, and sound intensity vortices are easily formed at the intersection of each reinforcement, which hinders the energy flow to the damping area and dissipates it in the center of the mirror.

[0106] It can be seen that in terms of vibration magnitude, as the vibration frequency increases, the vibration amplitude is in a decreasing trend; however, the same-frequency vibration magnitude of Scheme 2 is significantly higher than that of Scheme 1, and Scheme 2 will excite a significant vibration form at a higher frequency (23Hz) than Scheme 1. Therefore, through comparative calculations in this embodiment, it can be seen that the optical-mechanical design of Scheme 1 is better than that of Scheme 2.

[0107] The embodiment of the present application also provides a nonlinear structure sound intensity visualization device, comprising:

[0108] A simulation unit, which is used to input a nonlinear structure into a transient dynamics simulation model to obtain structural sound intensity simulation data of the nonlinear structure, wherein the structural sound intensity simulation data includes simulation result data of each entity unit in the nonlinear structure, and the simulation result data includes displacement response and stress response;

[0109] A processing unit, which is used to perform vector processing on the simulation result data of each entity unit to obtain the structural sound intensity vector field data of each entity unit;

[0110] A drawing unit, which is used to draw an acoustic intensity vector image of the nonlinear structure based on the structural acoustic intensity vector field data of each entity unit;

[0111] A visualization unit is used to import the sound intensity vector image into a visualization interaction interface.

[0112] It can be seen that the embodiment of the present application performs post-processing on the stress and displacement obtained by the finite element simulation to synthesize the structural sound intensity vector field data of each entity unit in the nonlinear structure, and draws and visualizes the sound intensity vector image, and then can accurately identify the vibration noise source in the nonlinear structure according to the visualized sound intensity vector image, and analyze the main propagation path of the vibration noise, thereby realizing precise directional optimization of the dynamic characteristics of the nonlinear structure, which does not require large-scale trial changes and iterative calculations, and effectively improves the efficiency and accuracy of the directional optimization of the dynamic characteristics of the nonlinear structure.

[0113] Furthermore, in an embodiment of the present application, the transient dynamics simulation model includes Rayleigh damping.

[0114] Furthermore, in the embodiment of the present application, the processing unit is specifically used for:

[0115] Performing time-frequency transformation on the stress response and displacement response of the first entity unit in the time domain respectively, to obtain the displacement response and stress response of the first entity unit corresponding to each frequency in the frequency domain;

[0116] A typical frequency is selected from multiple frequencies based on the frequency selection mechanism and the displacement response and stress response corresponding to each frequency;

[0117] Structural acoustic intensity vector field data of the first solid element is generated based on the displacement response and the stress response corresponding to the typical frequency.

[0118] Furthermore, in the embodiment of the present application, the processing unit is further configured to:

[0119] The displacement response and stress response corresponding to the typical frequency are calculated based on the structural sound intensity expression in the Cartesian coordinate system to obtain the structural sound intensity vector field data of the first entity unit.

[0120] Furthermore, in the embodiment of the present application, the structure-borne sound intensity expression in the Cartesian coordinate system is:

[0121]

[0122]

[0123]

[0124] Among them, I x ,I y and I z Respectively represent the structural sound intensity of the solid unit in the x-direction, y-direction and z-direction, σ x , σ y and σ z They represent the surface normal stress of the solid element in the x-direction, y-direction and z-direction respectively. and Respectively represent the conjugate displacement of the solid element in the x-direction, y-direction and z-direction, τ xy and τ yx represents the shear stress of the solid element on the xy plane, τ xz and τ zx represents the shear stress of the solid element on the xz plane, τ yz and τ zy It represents the shear stress of the solid element on the yz plane, and ω is the circular frequency of vibration.

[0125] Furthermore, in the embodiment of the present application, the drawing unit is specifically used for:

[0126] Based on the mapping function, the structural sound intensity vector field data of each entity unit is uniformly distributed to obtain the structural sound intensity vector field data after mapping changes corresponding to each entity unit;

[0127] The sound intensity vector image of the nonlinear structure is drawn according to the structural sound intensity vector field data after mapping changes corresponding to each entity unit.

[0128] Furthermore, in the embodiment of the present application, the expression of the mapping function is:

[0129]

[0130] Among them, I n represents the structural sound intensity vector field data after mapping changes, I represents the structural sound intensity vector field data before mapping changes, and I max Indicates the maximum value of the structure-borne sound intensity, I min Indicates the minimum value of the structure-borne sound intensity amplitude.

[0131] It should be noted that those skilled in the art can clearly understand that, for the convenience and simplicity of description, the specific working process of the above-described device and each unit can refer to the corresponding process in the aforementioned embodiment of the nonlinear structural sound intensity visualization method, and will not be repeated here.

[0132] The apparatus provided in the above embodiment may be implemented in the form of a computer program. The computer program may be Figure 5 The nonlinear structure-borne acoustic intensity visualization device is shown running on it.

[0133] The embodiment of the present application also provides a nonlinear structural sound intensity visualization method device, comprising: a memory, a processor and a network interface connected via a system bus, wherein at least one instruction is stored in the memory, and the at least one instruction is loaded and executed by the processor to implement all or part of the steps of the aforementioned nonlinear structural sound intensity visualization method.

[0134] The network interface is used for network communication, such as sending assigned tasks, etc. Those skilled in the art will understand that Figure 5 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0135] The processor may be a CPU, or other general-purpose processors, a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor, or the processor may be any conventional processor, etc. The processor is the control center of the computer device, and uses various interfaces and lines to connect various parts of the entire computer device.

[0136] The memory can be used to store computer programs and / or modules. The processor realizes various functions of the computer device by running or executing the computer programs and / or modules stored in the memory and calling the data stored in the memory. The memory can mainly include a program storage area and a data storage area, wherein the program storage area can store an operating system, an application required for at least one function (such as a video playback function, an image playback function, etc.), etc.; the data storage area can store data created according to the use of the mobile phone (such as video data, image data, etc.), etc. In addition, the memory can include a high-speed random access memory, and can also include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart memory card (SmartMedia Card, SMC), a secure digital (Secure digital, SD) card, a flash card (Flash Card), at least one disk storage device, a flash memory device or other volatile solid-state storage device.

[0137] The embodiment of the present application further provides a computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, all or part of the steps of the aforementioned nonlinear structural sound intensity visualization method are implemented.

[0138] The embodiment of the present application implements all or part of the aforementioned process, 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 each of the above methods when executed by the processor. Among them, 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 capable of carrying computer program code, recording medium, USB flash drive, 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. It should be noted that the content contained in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, the computer-readable medium does not include electric carrier signals and telecommunication signals.

[0139] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, servers or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage and optical storage, etc.) that contain computer-usable program code.

[0140] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems) and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0141] It should be noted that, in this article, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or system including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or system. In the absence of further restrictions, an element defined by the sentence "comprises a ..." does not exclude the existence of other identical elements in the process, method, article or system including the element.

[0142] The above description is only a specific implementation of the present application, so that those skilled in the art can understand or implement the present application. Various modifications to these embodiments will be apparent to those skilled in the art, and the general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present application. Therefore, the present application will not be limited to the embodiments shown herein, but will conform to the widest range consistent with the principles and novel features applied for herein.

Claims

1. A method for visualizing nonlinear structural sound intensity, characterized in that: The following steps are involved: Inputting a nonlinear structure into a transient dynamics simulation model to obtain structural sound intensity simulation data of the nonlinear structure, wherein the structural sound intensity simulation data includes simulation result data of each entity unit in the nonlinear structure, and the simulation result data includes displacement response and stress response; Perform vector processing on the simulation result data of each entity unit to obtain the structural sound intensity vector field data of each entity unit; Drawing the sound intensity vector image of the nonlinear structure based on the structural sound intensity vector field data of each entity unit; Importing the sound intensity vector image into a visual interactive interface; Wherein, the drawing of the sound intensity vector image of the nonlinear structure based on the structural sound intensity vector field data of each entity unit includes: Based on the mapping function, the structural sound intensity vector field data of each entity unit is uniformly distributed to obtain the structural sound intensity vector field data after mapping changes corresponding to each entity unit; Draw the sound intensity vector image of the nonlinear structure according to the structural sound intensity vector field data after mapping changes corresponding to each entity unit; The expression of the mapping function is: in, Represents the structural sound intensity vector field data after mapping changes, represents the structural sound intensity vector field data before mapping change, Indicates the maximum value of the structure-borne sound intensity amplitude, Indicates the minimum value of the structure-borne sound intensity amplitude.

2. The nonlinear structural sound intensity visualization method according to claim 1, characterized in that: The transient dynamics simulation model includes Rayleigh damping.

3. The nonlinear structural sound intensity visualization method according to claim 1, characterized in that: The vector processing of the simulation result data of each entity unit to obtain the structural sound intensity vector field data of each entity unit includes: Performing time-frequency transformation on the stress response and displacement response of the first entity unit in the time domain respectively, to obtain the displacement response and stress response of the first entity unit corresponding to each frequency in the frequency domain; A typical frequency is selected from multiple frequencies based on the frequency selection mechanism and the displacement response and stress response corresponding to each frequency; Structural acoustic intensity vector field data of the first solid element is generated based on the displacement response and the stress response corresponding to the typical frequency.

4. The nonlinear structural sound intensity visualization method according to claim 3, characterized in that: The generating the structural sound intensity vector field data of the first entity unit based on the displacement response and the stress response corresponding to the typical frequency includes: The displacement response and stress response corresponding to the typical frequency are calculated based on the structural sound intensity expression in the Cartesian coordinate system to obtain the structural sound intensity vector field data of the first entity unit.

5. The nonlinear structural sound intensity visualization method according to claim 4, characterized in that: The structure-borne sound intensity expression in the Cartesian coordinate system is: in, , and They represent the structural sound intensity of the solid element in the x-direction, y-direction and z-direction respectively. , and They represent the surface normal stress of the solid element in the x-direction, y-direction and z-direction respectively. , and They represent the conjugate of the displacement of the solid element in the x-direction, y-direction, and z-direction, respectively. and represents the shear stress of the solid element on the xy plane, and represents the shear stress of the solid element on the xz plane, and represents the shear stress of the solid element on the yz plane, is the circular frequency of vibration.

6. A nonlinear structural sound intensity visualization device, used to implement the nonlinear structural sound intensity visualization method according to any one of claims 1 to 5, characterized in that: include: A simulation unit, which is used to input a nonlinear structure into a transient dynamics simulation model to obtain structural sound intensity simulation data of the nonlinear structure, wherein the structural sound intensity simulation data includes simulation result data of each entity unit in the nonlinear structure, and the simulation result data includes displacement response and stress response; A processing unit, which is used to perform vector processing on the simulation result data of each entity unit to obtain the structural sound intensity vector field data of each entity unit; A drawing unit, which is used to draw an acoustic intensity vector image of the nonlinear structure based on the structural acoustic intensity vector field data of each entity unit; A visualization unit, used for importing the sound intensity vector image into a visualization interaction interface; The drawing unit is specifically used for: Based on the mapping function, the structural sound intensity vector field data of each entity unit is uniformly distributed to obtain the structural sound intensity vector field data after mapping changes corresponding to each entity unit; Draw the sound intensity vector image of the nonlinear structure according to the structural sound intensity vector field data after mapping changes corresponding to each entity unit; The expression of the mapping function is: in, Represents the structural sound intensity vector field data after mapping changes, represents the structural sound intensity vector field data before mapping change, Indicates the maximum value of the structure-borne sound intensity amplitude, Indicates the minimum value of the structure-borne sound intensity amplitude.

7. A nonlinear structural sound intensity visualization device, characterized in that: include: A memory and a processor, wherein at least one instruction is stored in the memory, and the at least one instruction is loaded and executed by the processor to implement the nonlinear structure sound intensity visualization method according to any one of claims 1 to 5.

8. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the nonlinear structure-borne sound intensity visualization method according to any one of claims 1 to 5 is implemented.

Images