Method for calculating three-dimensional edge diffraction sound pressure of rectangular sound barrier of transformer substation
By establishing a three-dimensional model and discretizing each edge of the substation rectangular acoustic barrier, the diffraction sound pressure formed by the sound source point through each diffraction point is calculated and superimposed, and it is solved by the difficulty in accurately estimating the rectangular acoustic barrier of the substation in the prior art. The problem of three-dimensional edge diffraction sound pressure achieves high-accuracy noise control.
Patent Information
- Application Number
- CN202411465995.9
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-10-21
- Publication Date
- 2025-02-11
AI Technical Summary
The prior art is difficult to accurately estimate the three-dimensional edge diffraction sound pressure of the rectangular acoustic barrier of the substation, resulting in sound insulation errors, and the traditional method is not suitable for rectangular acoustic barrier structures.
By establishing a three-dimensional model, including sound source point, sound barrier and target receiving point, each edge of the sound barrier is discretized, the diffraction sound pressure formed by the sound source point through each diffraction point is calculated, and superimposed on it Obtain the total diffraction sound pressure.
The accurate calculation of the three-dimensional edge diffraction sound pressure of the substation rectangular acoustic barrier is achieved, which reduces the sound insulation error, is suitable for the rectangular acoustic barrier structure, and improves the effect of noise pollution control.
Smart Images
Figure CN119416476A_ABST
Abstract
Description
Technical Field
[0001] The present application relates to the field of electric power environmental protection technology, and in particular to a method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier in a substation. Background Art
[0002] Substations are very important conversion sites in power systems. With the continuous acceleration of urbanization, more and more substations are being built in densely populated cities. During operation, substations cause serious noise pollution to urban areas where residents are relatively concentrated. In order to effectively control noise pollution from substations, sound barriers are built in substations to block the spread of noise. Therefore, the sound insulation of sound barriers has become a very important technical indicator. However, in actual use, the edges of sound barriers will produce diffraction, resulting in errors between their actual sound insulation and the sound insulation measured in the laboratory. Therefore, in the application of substation noise control, if the edge diffraction effect of the sound barrier can be effectively estimated, thereby reducing the sound insulation error caused by sound wave diffraction, it will be very beneficial to the more accurate design of substation sound barriers, thereby improving the effect of noise pollution control, which is of great significance for promoting the sustainable development of urban construction.
[0003] The noise emitted by substations has the characteristic of being low-frequency. Due to its long wavelength, it is easy to diffract at the edge of the sound barrier when it encounters it, that is, the sound wave bypasses the sound barrier and continues to propagate. There are currently three main methods for estimating the diffraction of obstacles: numerical solution, empirical formula method, and theoretical solution based on geometric diffraction formula.
[0004] The numerical solution is suitable for the calculation of low-frequency noise in three-dimensional scenes. However, since the space needs to be discretized into grids for processing, the amount of calculation is large. Especially when the scene size is large, the grid scale is large, which will cause a huge computational burden and is not conducive to on-site application in actual engineering.
[0005] Numerical solutions include finite element and boundary element methods, which refer to discretizing the problem domain using a grid, and then solving the problem based on finite element or boundary element theory. The empirical formula method refers to first conducting a typical sound barrier edge diffraction measurement experiment, and then using this experimental data as a benchmark to construct an empirical formula that has been fitted with data for sound diffraction estimation in other scenarios. This method has low accuracy and is not suitable for other scenarios with large environmental differences.
[0006] The theoretical solution is mainly a solution to estimate the diffraction amount based on the geometric diffraction formula. The geometric diffraction formula is a calculation method for the diffraction attenuation of radio waves on wedge-shaped obstacles proposed by Keller in 1951. Because of the similar principle, it is also often used in the diffraction estimation of sound waves. However, the classical geometric diffraction theory is a theory within a two-dimensional plane. In this case, the sound barrier needs to be processed in sections. The sound source point, diffraction point and receiving point are all within this plane, and the calculation results are also limited to this plane. In practical applications, the position of the target receiving point is random. There may be a situation where the sound source point, the sound barrier section and the receiving point are not within a two-dimensional plane. At this time, the traditional geometric diffraction formula is no longer applicable. The traditional geometric diffraction formula is a formula for wedge-shaped obstacles, but in practice, the substation sound barrier is usually a rectangular structure, which makes the traditional geometric diffraction formula also inapplicable.
[0007] Therefore, how to propose a method for calculating edge diffraction sound pressure in three-dimensional space and make it applicable to rectangular sound barrier structures is of great significance for improving its practicality. Summary of the invention
[0008] An embodiment of the present application provides a method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier in a substation to solve the technical problems existing in the prior art.
[0009] Other features and advantages of the present application will become apparent from the following detailed description, or may be learned in part by the practice of the present application.
[0010] According to a first aspect of an embodiment of the present application, a method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier in a substation is provided, comprising:
[0011] A three-dimensional model is obtained based on a real application scenario of a substation sound barrier, wherein the three-dimensional model includes a sound source point, a sound barrier and a target receiving point;
[0012] Discretizing each edge of the sound barrier in the three-dimensional model to obtain a plurality of diffraction points;
[0013] Calculate the diffracted sound pressure formed when the sound source passes through a certain diffraction point;
[0014] The diffracted sound pressure formed by the sound source point at each diffraction point is superimposed to obtain the diffracted sound pressure at the target receiving point after the sound source point is diffracted by the edge of the sound barrier.
[0015] In some embodiments of the present application, based on the above solution, each edge of the sound barrier in the three-dimensional model is discretized to obtain multiple diffraction points, including:
[0016] Each edge of the sound barrier in the three-dimensional model is evenly divided into a series of diffraction points, denoted as P: {P1, P2, ... P n}, where the subscript is the serial number of the point, and there are n points in total.
[0017] In some embodiments of the present application, based on the aforementioned solution, the distance between adjacent diffraction points is less than 10 cm.
[0018] In some embodiments of the present application, based on the above solution, the step of calculating the diffracted sound pressure formed when the sound source point passes through a certain diffraction point includes:
[0019] Determine a projection plane where the certain diffraction point is located, where the projection plane is perpendicular to the sound barrier;
[0020] Determine the projection points of the sound source point, the certain diffraction point and the target receiving point on the projection plane;
[0021] Forming an auxiliary plane based on the sound source point, the projection point of the sound source point, the target receiving point and the projection point of the target receiving point;
[0022] Determine the angular relationship between each projection point on the projection plane of the sound source point, the certain diffraction point and the target receiving point by using the projection plane and the auxiliary plane;
[0023] The diffracted sound pressure formed by the sound source point when passing through the certain diffraction point is calculated based on the distance relationship between the sound source point, the certain diffraction point and the target receiving point and the angle relationship between each projection point.
[0024] In some embodiments of the present application, based on the aforementioned solution, the calculating based on the distance relationship between the sound source point, the certain diffraction point and the target receiving point and the angle relationship between each projection point to obtain the diffracted sound pressure formed by the sound source point when passing through the certain diffraction point includes:
[0025] Obtaining a straight-line distance from the sound source point to the certain diffraction point, and calculating a direct sound pressure from the sound source point to the certain diffraction point;
[0026] Obtaining a straight-line distance from the certain diffraction point to the target receiving point, and calculating a direct sound pressure from the sound source point to the certain diffraction point;
[0027] Calculating a distance proportionality coefficient when the sound wave propagates through the certain diffraction point based on a straight-line distance from the sound source point to the certain diffraction point and a straight-line distance from the certain diffraction point to the target receiving point;
[0028] Obtaining the angular relationship between each projection point, and calculating the diffraction coefficient when the sound wave propagates through the diffraction point;
[0029] Based on the direct sound pressure from the sound source point to the certain diffraction point, the direct sound pressure from the sound source point to the certain diffraction point, the distance proportional coefficient when the sound wave propagates through the certain diffraction point, and the diffraction coefficient when the sound wave propagates through the certain diffraction point, the diffraction sound pressure p formed by the sound source point when passing through the certain diffraction point is calculated using formula (1): R ;
[0030]
[0031] In formula (1), Represents the distance from the sound source to the diffraction point P i The direct sound pressure, Indicates that the sound wave passes through the diffraction point P i The distance proportionality factor during propagation, Indicates that the sound wave passes through the diffraction point P i The diffraction coefficient during propagation is Represents the distance from the sound source to the diffraction point P i Direct sound pressure.
[0032] In some embodiments of the present application, based on the above solution, obtaining the straight-line distance from the sound source point to the certain diffraction point and calculating the direct sound pressure from the sound source point to the certain diffraction point includes:
[0033] Obtain the straight-line distance from the sound source point to the diffraction point, and use formula (2) to calculate the direct sound pressure from the sound source point to the diffraction point:
[0034]
[0035] In formula (1), j represents an imaginary unit, k = 2πf / c represents the wave number, f is the frequency, c is the speed of sound, Represents the distance from the sound source point S to a diffraction point P i The straight-line distance.
[0036] In some embodiments of the present application, based on the above solution, obtaining the straight-line distance from the certain diffraction point to the target receiving point and calculating the direct sound pressure from the sound source point to the certain diffraction point includes:
[0037] Obtain the straight-line distance from the diffraction point to the target receiving point, and use formula (3) to calculate the direct sound pressure from the sound source point to the diffraction point
[0038]
[0039] In formula (3), j represents an imaginary unit, k = 2πf / c represents the wave number, f is the frequency, and c is the speed of sound. Represents a diffraction point P i The straight-line distance to the receiving point R.
[0040] In some embodiments of the present application, based on the aforementioned solution, the distance proportional coefficient when the sound wave propagates through the certain diffraction point is calculated based on the straight-line distance from the sound source point to the certain diffraction point and the straight-line distance from the certain diffraction point to the target receiving point, including:
[0041] Based on the straight-line distance from the sound source point to the certain diffraction point and the straight-line distance from the certain diffraction point to the target receiving point, the distance proportional coefficient when the sound wave propagates through the certain diffraction point is calculated using formula (4):
[0042]
[0043] In formula (4), Indicates the distance from the sound source point S to a diffraction point P i The straight-line distance, Represents a diffraction point P i The straight-line distance to the receiving point R.
[0044] In some embodiments of the present application, based on the above solution, the step of obtaining the angle relationship between the projection points and calculating the diffraction coefficient when the sound wave propagates through the diffraction point includes:
[0045] Obtain the angular relationship between each projection point, and use formula (5) to calculate the diffraction coefficient when the sound wave propagates through the diffraction point
[0046]
[0047] In formula (5), j represents an imaginary unit, k represents a wave number, θ1 represents the angle between the projection point of the sound source point and the projection point of the certain diffraction point in the projection plane, θ2 represents the angle between the projection point of the certain diffraction point and the projection point of the target receiving point in the projection plane, and the expressions of M1, M2, M3, and M4 are as follows:
[0048]
[0049] In formula (6), k = 2πf / c represents the wave number, γ - With γ + represents the angle related to the incident angle, a + (γ + ), a - (γ+ ) is a function of the angle.
[0050] In some embodiments of the present application, based on the above solution, superimposing the diffracted sound pressure formed by the sound source point at each diffraction point includes:
[0051] Use formula (7) to superimpose the diffracted sound pressure formed by the sound source point at each diffraction point;
[0052]
[0053] In formula (7), spl represents the diffracted sound pressure at the target receiving point after the sound source point is diffracted by the edge of the sound barrier, and p R It represents the diffraction sound pressure formed when the sound source point passes through the diffraction point.
[0054] According to a second aspect of an embodiment of the present application, a computer-readable storage medium is provided, wherein the storage medium stores computer instructions, and when the computer instructions are executed on a computer, the computer executes the method according to the first aspect.
[0055] According to a third aspect of an embodiment of the present application, there is provided an electronic device, including: a memory and a processor;
[0056] The memory is used to store computer instructions;
[0057] The processor is used to call the computer instructions stored in the memory so that the electronic device executes the method as described in the first aspect.
[0058] The technical solution of the present application comprehensively considers the relationship between factors such as the sound source point, diffraction point, and receiving point, and can calculate the sound diffraction in three dimensions, solving the problem that traditional methods can only solve two-dimensional diffraction.
[0059] It should be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the present application. BRIEF DESCRIPTION OF THE DRAWINGS
[0060] The drawings herein are incorporated into the specification and constitute a part of the specification, showing embodiments consistent with the present application, and together with the specification, are used to explain the principles of the present application. Obviously, the drawings described below are only some embodiments of the present application, and for ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work. In the drawings:
[0061] Figure 1 A schematic diagram of a process for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier in a substation according to an embodiment of the present application is shown;
[0062] Figure 2 A schematic diagram of three-dimensional modeling of a substation sound barrier according to an embodiment of the present application is shown;
[0063] Figure 3 The acoustic wave according to one embodiment of the present application passes through the diffraction point P. i Schematic diagram of the diffraction of ;
[0064] Figure 4 A schematic diagram showing angles between projection points in a projection plane according to an embodiment of the present application is shown;
[0065] Figure 5 A schematic diagram of angles within an auxiliary plane according to an embodiment of the present application is shown;
[0066] Figure 6 A schematic diagram of a result verification model according to an embodiment of the present application is shown;
[0067] Figure 7 A schematic diagram showing comparison of verification results according to an embodiment of the present application is shown;
[0068] Figure 8 A block diagram of an electronic device according to an embodiment of the present application is shown;
[0069] Fig. 9 A schematic diagram of the structure of a computer system suitable for implementing an electronic device of an embodiment of the present application is shown. DETAILED DESCRIPTION
[0070] Example embodiments will now be described more fully with reference to the accompanying drawings. However, example embodiments can be implemented in a variety of forms and should not be construed as limited to the examples set forth herein; rather, these embodiments are provided so that this application will be more comprehensive and complete and fully convey the concept of the example embodiments to those skilled in the art.
[0071] In addition, described feature, structure or characteristic can be combined in one or more embodiments in any suitable manner. In the following description, many specific details are provided to provide a full understanding of the embodiments of the present application. However, those skilled in the art will appreciate that the technical scheme of the present application can be put into practice without one or more of the specific details, or other methods, components, devices, steps, etc. can be adopted. In other cases, known methods, devices, realizations or operations are not shown or described in detail to avoid blurring the various aspects of the application.
[0072] The block diagrams shown in the accompanying drawings are merely functional entities and do not necessarily correspond to physically independent entities. That is, these functional entities may be implemented in software form, or in one or more hardware modules or integrated circuits, or in different networks and / or processor devices and / or microcontroller devices.
[0073] The flowcharts shown in the accompanying drawings are only exemplary and do not necessarily include all the contents and operations / steps, nor must they be executed in the order described. For example, some operations / steps can be decomposed, and some operations / steps can be combined or partially combined, so the actual execution order may change according to actual conditions.
[0074] It should be noted that the “plurality” mentioned in this article refers to two or more.
[0075] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution in the embodiment of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiment of the present invention. Obviously, the described embodiment is only a part of the embodiment of the present invention, not all the embodiments. Based on the embodiment of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.
[0076] Some embodiments of the present application will be described in detail below in conjunction with the accompanying drawings. In the absence of conflict, the following embodiments and features in the embodiments may be combined with each other.
[0077] See also Figure 1 , shows a flow chart of a method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier in a substation according to an embodiment of the present application.
[0078] like Figure 1 As shown, a method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier in a substation is presented, which specifically includes steps S100 to S400.
[0079] refer to Figure 1 , step S100, modeling is performed based on the real application scenario of the substation sound barrier to obtain a three-dimensional model, wherein the three-dimensional model includes a sound source point, a sound barrier and a target receiving point.
[0080] For example, see Figure 2 , shows a schematic diagram of three-dimensional modeling of a substation sound barrier. Figure 2As shown, there is a sound barrier in this scene, and its four vertices are ABCD, A and B represent the two bottom vertices of the sound barrier, C and D represent the two top vertices of the sound barrier, S represents the sound source, R represents the target receiving point, and plane O1 is the plane passing through the bottom edge AB of the sound barrier and perpendicular to the sound barrier.
[0081] Continue to refer Figure 1 , step S200, discretize each edge of the sound barrier in the three-dimensional model to obtain multiple diffraction points.
[0082] It should be noted that in Figure 2 In the scenario shown, the sound wave will diffract along the edge of the sound barrier, and the diffracted wave will form a sound field through point R. Since the sound wave has the characteristic of continuous expansion of the wavefront, it will pass through all the edges of the sound barrier. In order to calculate the impact of the edge of the sound barrier on the sound field of the target point, the present invention discretizes the edge, that is, assumes the edge to be a series of points.
[0083] In some feasible embodiments, based on the above solution, each edge of the sound barrier in the three-dimensional model is discretized to obtain multiple diffraction points, including:
[0084] Each edge of the sound barrier in the three-dimensional model is evenly divided into a series of diffraction points, denoted as P: {P1, P2, ... P n}, where the subscript is the serial number of the point, and there are n points in total.
[0085] For example, Figure 2 As shown, the edges where edge diffraction will occur in this scene are AC, CD, and BD. Then these three edges are discretized into a series of points, represented by P: {P1, P2, …P n}, where the subscript is the serial number of the point, there are n points in total, and the distance between adjacent points is less than 10cm.
[0086] Continue to refer Figure 1 , step S300, calculating the diffracted sound pressure formed when the sound source point passes through a certain diffraction point.
[0087] In some feasible embodiments, based on the above solution, the step of calculating the diffracted sound pressure formed when the sound source point passes through a certain diffraction point includes:
[0088] Determine a projection plane where the certain diffraction point is located, where the projection plane is perpendicular to the sound barrier;
[0089] Determine the projection points of the sound source point, the certain diffraction point and the target receiving point on the projection plane;
[0090] Forming an auxiliary plane based on the sound source point, the projection point of the sound source point, the target receiving point and the projection point of the target receiving point;
[0091] Determine the angular relationship between each projection point on the projection plane of the sound source point, the certain diffraction point and the target receiving point by using the projection plane and the auxiliary plane;
[0092] The diffracted sound pressure formed by the sound source point when passing through the certain diffraction point is calculated based on the distance relationship between the sound source point, the certain diffraction point and the target receiving point and the angle relationship between each projection point.
[0093] It can be understood that, since the sound barrier has multiple different sides, when the diffraction point is on different sides of the sound barrier, the position of the projection surface is also different.
[0094] For example, see Figure 3 , showing the sound wave passing through the diffraction point P i Schematic diagram of the diffraction of . Figure 3 As shown, P i On the CD side, the projection surface is the surface O perpendicular to the sound barrier. If P i When it is on the AC side, its projection surface should be a plane passing through the BD side and perpendicular to the sound barrier. Figure 3 As shown, the projection of S on plane O1 is S', the projection of R on plane O1 is R', and P i The projection on plane O is P i '. Make an auxiliary plane through the four points S, S', R, and R'. The intersection point of this auxiliary plane with the side CD is J, and the intersection point with the side AB is J'.
[0095] In some feasible embodiments, based on the above scheme, the calculating based on the distance relationship between the sound source point, the certain diffraction point and the target receiving point and the angle relationship between each projection point to obtain the diffracted sound pressure formed by the sound source point when passing through the certain diffraction point includes:
[0096] Obtaining a straight-line distance from the sound source point to the certain diffraction point, and calculating a direct sound pressure from the sound source point to the certain diffraction point;
[0097] Obtaining a straight-line distance from the certain diffraction point to the target receiving point, and calculating a direct sound pressure from the sound source point to the certain diffraction point;
[0098] Calculating a distance proportionality coefficient when the sound wave propagates through the certain diffraction point based on a straight-line distance from the sound source point to the certain diffraction point and a straight-line distance from the certain diffraction point to the target receiving point;
[0099] Obtaining the angular relationship between each projection point, and calculating the diffraction coefficient when the sound wave propagates through the diffraction point;
[0100] Based on the direct sound pressure from the sound source point to the certain diffraction point, the direct sound pressure from the sound source point to the certain diffraction point, the distance proportional coefficient when the sound wave propagates through the certain diffraction point, and the diffraction coefficient when the sound wave propagates through the certain diffraction point, the diffraction sound pressure p formed by the sound source point when passing through the certain diffraction point is calculated using formula (1): R ;
[0101]
[0102] In formula (1), Represents the distance from the sound source to the diffraction point P i The direct sound pressure, Indicates that the sound wave passes through the diffraction point P i The distance proportionality factor during propagation, Indicates that the sound wave passes through the diffraction point P i The diffraction coefficient during propagation is Represents the distance from the sound source to the diffraction point P i Direct sound pressure.
[0103] In some feasible embodiments, based on the above solution, obtaining the straight-line distance from the sound source point to the certain diffraction point and calculating the direct sound pressure from the sound source point to the certain diffraction point includes:
[0104] Obtain the straight-line distance from the sound source point to the diffraction point, and use formula (2) to calculate the direct sound pressure from the sound source point to the diffraction point:
[0105]
[0106] In formula (1), j represents an imaginary unit, k = 2πf / c represents the wave number, f is the frequency, c is the speed of sound, Represents the distance from the sound source point S to a diffraction point P i The straight-line distance.
[0107] In some feasible embodiments, based on the above solution, obtaining the straight-line distance from the certain diffraction point to the target receiving point and calculating the direct sound pressure from the sound source point to the certain diffraction point includes:
[0108] Obtain the straight-line distance from the diffraction point to the target receiving point, and use formula (3) to calculate the direct sound pressure from the sound source point to the diffraction point
[0109]
[0110] In formula (3), j represents an imaginary unit, k = 2πf / c represents the wave number, f is the frequency, c is the speed of sound, Represents a diffraction point P i The straight-line distance to the receiving point R.
[0111] In some feasible embodiments, based on the above scheme, the distance proportional coefficient when the sound wave propagates through the certain diffraction point is calculated based on the straight-line distance from the sound source point to the certain diffraction point and the straight-line distance from the certain diffraction point to the target receiving point, including:
[0112] Based on the straight-line distance from the sound source point to the certain diffraction point and the straight-line distance from the certain diffraction point to the target receiving point, the distance proportional coefficient when the sound wave propagates through the certain diffraction point is calculated using formula (4):
[0113]
[0114] In formula (4), Indicates the distance from the sound source point S to a diffraction point P i The straight-line distance, Represents a diffraction point P i The straight-line distance to the receiving point R.
[0115] In some feasible embodiments, based on the above solution, obtaining the angle relationship between the projection points and calculating the diffraction coefficient when the sound wave propagates through the certain diffraction point includes:
[0116] Obtain the angular relationship between each projection point, and use formula (5) to calculate the diffraction coefficient when the sound wave propagates through the diffraction point
[0117]
[0118] In formula (5), j represents an imaginary unit, k represents a wave number, θ1 represents the angle between the projection point of the sound source point and the projection point of the certain diffraction point in the projection plane, θ2 represents the angle between the projection point of the certain diffraction point and the projection point of the target receiving point in the projection plane, and the expressions of M1, M2, M3, and M4 are as follows:
[0119]
[0120] In formula (6), k = 2πf / c represents the wave number, γ - With γ + represents the angle related to the incident angle, a + (γ + ), a - (γ + ) is a function of the angle.
[0121] For example, see Figure 4 , which shows a schematic diagram of the angles between the projection points in the projection plane. Figure 4 As shown, θ1 represents the difference between the projection point S' and the projection point P i ', where θ1 = ∠P i ′S′J′; θ2 represents the projection point P i ' and the angle between the projection point R', where θ2 = ∠P i 'R'J'.
[0122] In formula (6), cot(x) represents the cotangent function, and F(x) is a variation of the Fresnel integral, defined as:
[0123]
[0124] Where j represents the imaginary unit and τ is the independent variable of the integral.
[0125] γ - With γ + represents the angle related to the incident angle and is expressed as:
[0126]
[0127] Where θ3 and θ4 represent the angles related to incidence and diffraction in the auxiliary plane.
[0128] For example, see Figure 5 , shows a schematic diagram of the angles in the auxiliary plane. Figure 5 As shown, θ3 = 2π-∠RJJ′, θ4 = ∠SJJ′.
[0129] In formula (6), a + (γ + ), a - (γ + ) is as follows:
[0130]
[0131] In formula (9), N + and N - is a variable related to the x value, and its value is:
[0132]
[0133] Continue to refer Figure 1 In step S400, the diffracted sound pressure formed by the sound source point at each diffraction point is superimposed to obtain the diffracted sound pressure at the target receiving point after the sound source point is diffracted by the edge of the sound barrier.
[0134] In some feasible embodiments, based on the above solution, superimposing the diffracted sound pressure formed by the sound source point at each diffraction point includes:
[0135] The diffracted sound pressure formed by the sound source point at each diffraction point is superimposed using formula (11);
[0136]
[0137] In formula (11), spl represents the diffracted sound pressure at the target receiving point after the sound source point is diffracted by the edge of the sound barrier, and p R It represents the diffraction sound pressure formed when the sound source point passes through the diffraction point.
[0138] It should be noted that after completing the calculation of the sound pressure at all diffraction points, the diffracted sound pressure at the target receiving point after the sound source point is diffracted by the edge of the sound barrier can be calculated using formula (11).
[0139] The following is an example to verify this method.
[0140] by Figure 6 The three-dimensional sound barrier scene shown in the figure was simulated and verified. The sound barrier is 3m long, 3m high and 10cm thick. The projection of the sound source on the sound barrier is at its center, the sound source is 1m away from the sound barrier, and its sound power level is assumed to be 100dB. The projection of the receiving point on the sound barrier is also at its center, and the receiving point is 1m away from the sound barrier. In order to verify the accuracy of the present invention, the boundary element method was used to calculate this scene. The grid divided by the boundary element is shown in the figure. The calculation results are compared as shown in the figure. Figure 7 shown.
[0141] Figure 7 The calculation results shown show that the diffraction simulation algorithm proposed in the present invention can handle three-dimensional diffraction problems well. Compared with the boundary element method, the error at the main frequency point does not exceed 2dB, but the calculation time is much shorter than that of the boundary element method.
[0142] To summarize, this application aims at the three-dimensional diffraction estimation problem of substation sound barriers. According to the relative positions of the sound source point, a point on the edge of the sound barrier, and the target receiving point, a method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier in a substation is proposed. This method, based on the rectangular sound barrier, comprehensively considers the relationship between factors such as the sound source point, diffraction point, and receiving point. It can calculate the sound diffraction in three-dimensional cases and can solve the rectangular sound barriers that are more common in engineering, solving the problem that traditional methods can only solve wedge-shaped obstacles.
[0143] like Figure 8As shown, an embodiment of the present application also provides an electronic device 800, including a memory 810, a processor 820, and a computer program 811 stored in the memory 810 and executable on the processor. When the processor 820 executes the computer program 811, the steps of the above-mentioned method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier in a substation are implemented.
[0144] Since the electronic device introduced in this embodiment is a device used to implement a method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier in a substation in the embodiment of the present application, based on the method introduced in the embodiment of the present application, technical personnel in this field can understand the specific implementation mode of the electronic device of this embodiment and its various variations. Therefore, how the electronic device implements the method in the embodiment of the present application is not introduced in detail here. As long as the equipment used by technical personnel in this field to implement the method in the embodiment of the present application is within the scope of protection of this application.
[0145] During the specific implementation process, when the computer program 811 is executed by a processor, any implementation method in the embodiments corresponding to the first aspect can be implemented.
[0146] Fig. 9 A schematic diagram of the structure of a computer system suitable for implementing an electronic device of an embodiment of the present application is shown.
[0147] It should be noted that Fig. 9 The computer system 900 of the electronic device shown is only an example and should not bring any limitation to the functions and scope of use of the embodiments of the present application.
[0148] like Fig. 9 As shown, the computer system 900 includes a central processing unit (CPU) 901, which can perform various appropriate actions and processes according to the program stored in the read-only memory (ROM) 902 or the program loaded from the storage part 908 to the random access memory (RAM) 903, such as executing the method described in the above embodiment. In the RAM 903, various programs and data required for system operation are also stored. The CPU 901, ROM 902 and RAM 903 are connected to each other through a bus 904. An input / output (I / O) interface 905 is also connected to the bus 904.
[0149] The following components are connected to the I / O interface 905: an input section 906 including a keyboard, a mouse, etc.; an output section 907 including a cathode ray tube (CRT), a liquid crystal display (LCD), etc., and a speaker, etc.; a storage section 908 including a hard disk, etc.; and a communication section 909 including a network interface card such as a LAN (Local Area Network) card, a modem, etc. The communication section 909 performs communication processing via a network such as the Internet. A drive 910 is also connected to the I / O interface 905 as needed. A removable medium 911, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is installed on the drive 910 as needed so that a computer program read therefrom is installed into the storage section 908 as needed.
[0150] In particular, according to an embodiment of the present application, the process described above with reference to the flowchart can be implemented as a computer software program. For example, an embodiment of the present application includes a computer program product, which includes a computer program carried on a computer-readable medium, and the computer program includes a program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from a network through a communication section 909, and / or installed from a removable medium 911. When the computer program is executed by a central processing unit (CPU) 901, various functions defined in the system of the present application are executed.
[0151] It should be noted that the computer-readable medium shown in the embodiment of the present application may be a computer-readable signal medium or a computer-readable storage medium or any combination of the above two. The computer-readable storage medium may be, for example, - but not limited to - an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination of the above. More specific examples of computer-readable storage media may include, but are not limited to: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM), a flash memory, an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In the present application, a computer-readable storage medium may be any tangible medium containing or storing a program, which may be used by an instruction execution system, device or device or used in combination with it. In the present application, a computer-readable signal medium may include a data signal propagated in a baseband or as part of a carrier wave, wherein a computer-readable program code is carried. Such propagated data signals may take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. Computer-readable signal media may also be any computer-readable medium other than computer-readable storage media, which may send, propagate, or transmit programs for use by or in conjunction with an instruction execution system, apparatus, or device. The program code contained on the computer-readable medium may be transmitted using any appropriate medium, including but not limited to: wireless, wired, etc., or any suitable combination of the above.
[0152] The flowchart and block diagram in the accompanying drawings illustrate the possible architecture, functions and operations of the system, method and computer program product according to various embodiments of the present application. Wherein, each box in the flowchart or block diagram can represent a module, a program segment, or a part of the code, and the above-mentioned module, program segment, or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the box can also occur in a different order from the order marked in the accompanying drawings. For example, two boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram or flowchart, and the combination of boxes in the block diagram or flowchart can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.
[0153] The units involved in the embodiments described in this application may be implemented by software or hardware, and the units described may also be set in a processor. The names of these units do not, in some cases, constitute limitations on the units themselves.
[0154] As another aspect, the present application also provides a computer program product or a computer program, which includes a computer instruction stored in a computer-readable storage medium. A processor of a computer device reads the computer instruction from the computer-readable storage medium, and the processor executes the computer instruction, so that the computer device executes a method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier of a substation described in the above embodiment.
[0155] As another aspect, the present application also provides a computer-readable medium, which may be included in the electronic device described in the above embodiment; or may exist independently without being assembled into the electronic device. The above computer-readable medium carries one or more programs, and when the above one or more programs are executed by an electronic device, the electronic device implements a method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier of a substation described in the above embodiment.
[0156] It should be noted that, although several modules or units of the equipment for action execution are mentioned in the above detailed description, this division is not mandatory. In fact, according to the embodiments of the present application, the features and functions of two or more modules or units described above can be embodied in one module or unit. On the contrary, the features and functions of one module or unit described above can be further divided into being embodied by multiple modules or units.
[0157] Through the description of the above implementation methods, it is easy for those skilled in the art to understand that the example implementation methods described here can be implemented by software or by combining software with necessary hardware. Therefore, the technical solution according to the implementation methods of the present application can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (which can be a CD-ROM, a USB flash drive, a mobile hard disk, etc.) or on a network, and includes several instructions to enable a computing device (which can be a personal computer, a server, a touch terminal, or a network device, etc.) to execute the method according to the implementation methods of the present application.
[0158] Those skilled in the art will readily appreciate other embodiments of the present application after considering the specification and practicing the embodiments disclosed herein. The present application is intended to cover any variations, uses or adaptations of the present application, which follow the general principles of the present application and include common knowledge or customary techniques in the art that are not disclosed in the present application. It should be understood that the present application is not limited to the precise structures described above and shown in the accompanying drawings, and various modifications and changes can be made without departing from the scope thereof. The scope of the present application is limited only by the appended claims.
Claims
1. A method for calculating the three-dimensional edge diffraction sound pressure of a rectangular sound barrier in a substation, characterized in that: include: A three-dimensional model is obtained based on a real application scenario of a substation sound barrier, wherein the three-dimensional model includes a sound source point, a sound barrier and a target receiving point; Discretizing each edge of the sound barrier in the three-dimensional model to obtain a plurality of diffraction points; Calculate the diffracted sound pressure formed when the sound source passes through a certain diffraction point; The diffracted sound pressure formed by the sound source point at each diffraction point is superimposed to obtain the diffracted sound pressure at the target receiving point after the sound source point is diffracted by the edge of the sound barrier.
2. The method according to claim 1, characterized in that The discretization process is performed on each edge of the sound barrier in the three-dimensional model to obtain a plurality of diffraction points, including: Each edge of the sound barrier in the three-dimensional model is evenly divided into a series of diffraction points, denoted as P: {P1, P2, ... P n }, where the subscript is the serial number of the point, and there are n points in total.
3. The method according to claim 2, characterized in that The distance between adjacent diffraction points is less than 10 cm.
4. The method according to claim 1, characterized in that: The calculating of the diffracted sound pressure formed when the sound source point passes through a diffraction point includes: Determine a projection plane where the certain diffraction point is located, where the projection plane is perpendicular to the sound barrier; Determine the projection points of the sound source point, the certain diffraction point and the target receiving point on the projection plane; Forming an auxiliary plane based on the sound source point, the projection point of the sound source point, the target receiving point and the projection point of the target receiving point; Determine the angular relationship between each projection point on the projection plane of the sound source point, the certain diffraction point and the target receiving point by using the projection plane and the auxiliary plane; The diffracted sound pressure formed by the sound source point when passing through the certain diffraction point is calculated based on the distance relationship between the sound source point, the certain diffraction point and the target receiving point and the angle relationship between each projection point.
5. The method according to claim 4, characterized in that The calculating, based on the distance relationship between the sound source point, the certain diffraction point and the target receiving point and the angle relationship between each projection point, the diffraction sound pressure formed by the sound source point when passing through the certain diffraction point includes: Obtaining a straight-line distance from the sound source point to the certain diffraction point, and calculating a direct sound pressure from the sound source point to the certain diffraction point; Obtaining a straight-line distance from the certain diffraction point to the target receiving point, and calculating a direct sound pressure from the sound source point to the certain diffraction point; Calculating a distance proportionality coefficient when the sound wave propagates through the certain diffraction point based on a straight-line distance from the sound source point to the certain diffraction point and a straight-line distance from the certain diffraction point to the target receiving point; Obtaining the angular relationship between each projection point, and calculating the diffraction coefficient when the sound wave propagates through the diffraction point; Based on the direct sound pressure from the sound source point to the certain diffraction point, the direct sound pressure from the sound source point to the certain diffraction point, the distance proportional coefficient when the sound wave propagates through the certain diffraction point, and the diffraction coefficient when the sound wave propagates through the certain diffraction point, the diffraction sound pressure p formed by the sound source point when passing through the certain diffraction point is calculated using formula (1): R ; In formula (1), Represents the distance from the sound source point to the diffraction point P i The direct sound pressure, Indicates that the sound wave passes through the diffraction point P i The distance proportionality factor during propagation, Indicates that the sound wave passes through the diffraction point P i The diffraction coefficient during propagation, Represents the distance from the sound source point to the diffraction point P i Direct sound pressure.
6. The method according to claim 5, characterized in that The obtaining of the straight-line distance from the sound source point to the certain diffraction point and the calculation of the direct sound pressure from the sound source point to the certain diffraction point includes: Obtain the straight-line distance from the sound source point to the diffraction point, and use formula (2) to calculate the direct sound pressure from the sound source point to the diffraction point: In formula (1), j represents an imaginary unit, k = 2πf / c represents the wave number, f is the frequency, c is the speed of sound, Indicates the distance from the sound source point S to a diffraction point P i The straight-line distance.
7. The method according to claim 5, characterized in that The obtaining of the straight-line distance from the certain diffraction point to the target receiving point and the calculation of the direct sound pressure from the sound source point to the certain diffraction point includes: Obtain the straight-line distance from the diffraction point to the target receiving point, and use formula (3) to calculate the direct sound pressure from the sound source point to the diffraction point In formula (3), j represents an imaginary unit, k = 2πf / c represents the wave number, f is the frequency, and c is the speed of sound. Represents a diffraction point P i The straight-line distance to the receiving point R.
8. The method according to claim 5, characterized in that The calculating of the distance proportional coefficient when the sound wave propagates through the certain diffraction point based on the straight-line distance from the sound source point to the certain diffraction point and the straight-line distance from the certain diffraction point to the target receiving point includes: Based on the straight-line distance from the sound source point to the certain diffraction point and the straight-line distance from the certain diffraction point to the target receiving point, the distance proportional coefficient when the sound wave propagates through the certain diffraction point is calculated using formula (4): In formula (4), Indicates the distance from the sound source point S to a diffraction point P i The straight-line distance, Represents a diffraction point P i The straight-line distance to the receiving point R.
9. The method according to claim 5, characterized in that The step of obtaining the angle relationship between the projection points and calculating the diffraction coefficient when the sound wave propagates through the diffraction point includes: Obtain the angular relationship between each projection point, and use formula (5) to calculate the diffraction coefficient when the sound wave propagates through the diffraction point In formula (5), j represents an imaginary unit, k represents a wave number, θ1 represents the angle between the projection point of the sound source point and the projection point of the certain diffraction point in the projection plane, θ2 represents the angle between the projection point of the certain diffraction point and the projection point of the target receiving point in the projection plane, and the expressions of M1, M2, M3, and M4 are as follows: In formula (6), k = 2πf / c represents the wave number, γ - With γ + represents the angle related to the incident angle, a + (γ + ), a - (γ + ) is a function of the angle.
10. The method according to claim 1, characterized in that The step of superimposing the diffracted sound pressure formed by the sound source point at each diffraction point comprises: Use formula (7) to superimpose the diffracted sound pressure formed by the sound source point at each diffraction point; In formula (7), spl represents the diffracted sound pressure at the target receiving point after the sound source point is diffracted by the edge of the sound barrier, and p R It represents the diffraction sound pressure formed when the sound source point passes through the diffraction point.