The present disclosure belongs to the technical field of an ultrasonic sound field measuring system, and relates to a manual acoustic axis aligning method, in particular to a manual acoustic axis aligning method for an ultrasonic sound field measuring system.
A hydrophone method, is one of the standard methods of measuring sound fields. The method uses a hydrophone to perform three-dimensional scanning on a target sound field in a multi-axis water tank system to obtain the sound field distribution, and then calculates a plurality of acoustic parameters based thereon. The method can measure many parameters and is most widely used. In order to scan the three-dimensional sound, field, the multi-axis water tank system comprises at least three moving axes perpendicular to each other, as shown in
In order to realize the alignment of acoustic axes, the multi-axis water tank system adds two angle adjustment axes perpendicular to each other on the basis of three linear motion axes to adjust the angle of acoustic axes in space. At present, there are two categories of multi-axis water tank systems: a fully automatic five-axis system and a hybrid five-axis (automatic three-axis and manual two-axis) (angle adjustment) system. The fully automatic five-axis water tank system, such as AIMS system of ONDA, uses an angle adjusting shaft driven by a stepping motor via a belt to align the acoustic axis. The whole system operates under the control of a computer. The computer controls the rotation angle of each rotating shaft. After each adjustment, the acoustic axis angle is measured, and the automatic alignment of the acoustic axis is completed by a feedback iterative algorithm. It is convenient to use AIMS, but its application is limited because of its high price. The hybrid five-axis system uses two manual axes to adjust the acoustic axis angle and realize the alignment of the acoustic axis, and is high in cost performance, which is widely used in the field of sound field measurement. Such system needs to adjust the angle manually in the process of acoustic axis alignment. Compared with automatic adjustment, the biggest problem is that the angle of each adjustment cannot be accurately obtained, so that a manual adjusting method must be used to gradually align the acoustic axis.
As shown in
The purpose of the present disclosure is to overcome the shortcomings of the prior art and provide a manual acoustic axis aligning method for an ultrasonic sound field measuring system, so as to improve the reliability of the acoustic axis alignment and reduce the number of iterative adjustments, and avoid the problems of long measurement time, large measurement error and poor measurement repeatability.
To achieve the above purpose, the present disclosure provides the following technical scheme.
A manual acoustic axis aligning method for an ultrasonic sound field measuring system is provided, specifically comprising the following cycles:
1) determining three planes according to the focus position, wherein the three planes are a focusing plane Zfocus, an acoustic axis angle detection plane Z1 and an acoustic axis angle detection plane Z2; defining the distances from the probe surface to the focusing plane Zfocus, the acoustic axis angle detection plane Z1 and the acoustic, axis angle detection plane Z2 as Dfocus, a and b, respectively;
2) determining a peak point of an underwater acoustic signal of the focusing plane Zfocus and taking the peak point as an initial center point PC1, and denoting the coordinate of the initial center point PC1 as [PC_X1, PC_Y1]:
3) in the acoustic axis angle detection plane Z1, scanning with PCi as the center, changing the coordinates of X and Y, and, acquiring sound field signals to obtain the intersection C1=[C1X, C1Y] between the acoustic axis and the acoustic axis angle detection plane Z1, where i is a positive integer;
4) in the acoustic axis angle detection plane Z2, scanning with PCi as the center, changing the coordinates of X and Y, and acquiring sound field signals to obtain the intersection C2=[C2X, C2Y] between the acoustic axis and the acoustic axis angle detection plane Z2, where i is a positive integer;
5) calculating the inclination angle θi=[θxi, θyi] of the current acoustic axis according to the intersection points C1 and C2; θxi is the included angle between the acoustic axis and Z axis after the acoustic axis is projected on XOZ, plane, and θyi is the included angle between the acoustic axis and Z axis after the acoustic axis is projected on YOZ plane;
6) judging whether the inclination angle θi is less than the adjustment threshold, if so, aligning the acoustic axis, ending the adjustment, otherwise, proceeding to step 7);
7) according to the inclination angle θi, calculating the next coordinate PCi=[PC_Xi, PC_Yi] by using an iterative algorithm, where i is an integer not less than 2;
8) keeping PCi of the X axis and Y axis unchanged, moving the Z axis to the focusing plane Zfocus, manually adjusting two corner adjustment knobs of a probe clamping device until the amplitude of a hydrophone measurement signal at the current coordinate point is the maximum, and then proceeding to step 3).
Further, in step 1), the acoustic axis angle detection plane Z1 and the acoustic axis angle detection plane Z2 are located on the same side or different sides of the focusing plane Zfocus.
Further, in step 2), the peak point is the maximum point of the underwater acoustic signal of the focusing plane.
Further, in step 3) and step 4), according to the acquired sound field signals, the intersection C1=[C1X, C1Y] between the acoustic axis and the acoustic axis angle detection plane Z1, and the intersection C2=[C2X, C1Y] between the acoustic axis and the acoustic axis angle detection plane Z2 are obtained by a simple amplitude criterion or the criterion for calculating the symmetry center based on symmetry, respectively.
Further, the simple amplitude criterion is the amplitude maximum criterion, and the criterion for calculating the symmetry center based on symmetry is the centroid criterion.
Further, in step 5), the specific calculation formula of the coordinate value of the inclination angle θi is as follows:
θxi=arcctan((C2X−C1X)/(b−a)) (1)
θyi=arcctan((C2Y−C1Y)/(b−a)) (2)
Further, in step 7), calculating the next coordinate PCi=[PC_Xi, PC_Yi] by using an iterative algorithm has the specific calculation process as follows:
if the inclination angle of the current acoustic axis is obtained after the first adjustment, calculating the adjustment amount according to the triangular relationship, the coordinate PC2=[PC_X2, PC_Y2] of a new center point is calculated according to the depth Dfocus of the focusing plane Zfocus with the following formula:
PC_X2=PC_X1+tan(θx1)×(Dfocus−a) (3)
PC_Y2=PC_Y1+tan(θy1)×(Dfocus−a) (4)
in the above formulas (3) and (4), the inclination angle θ1=[θx1, θy1] of the current acoustic axis is obtained after the first adjustment;
otherwise, the calculation formula of the coordinate PCi+1=[PC_Xi+1, PC_Yi+1] of the new center point is as follows:
PC_Xi+1=(PC_Xi−PC_Xi−1)/(θxi−θxi−1)×θxi (5)
PC_Yi+1=(PC_Yi−PC_Yi−1)/(θyi−θyi−1)×θyi (6)
in the above formulas (5) and (6), i is an integer not less than 2, the inclination angle of the current acoustic axis is θi=[θxi, θyi], and the inclination angle before the inclination angle of the current acoustic axis is adjusted is θi−1=[θxi−1, θyi−1].
Further, the order of step 3) and the order of step 4) are interchangeable.
Further, in step 3) and step 4), a computer controls a water tank system to automatically scan the acoustic axis angle detection plane Z1 or the acoustic axis angle detection plane Z2.
Further, the larger the interval between the acoustic axis angle detection plane Z1 and the acoustic axis angle detection plane Z2, the higher the alignment accuracy.
Compared with the prior art, the technical scheme according to the present disclosure has the following beneficial effects: first, the depths of the acoustic axis angle detection plane Z1 and the acoustic axis angle detection plane Z2 are randomly chosen, which can be both in front of the focusing plane, on both sides of the focusing plane or both behind the focusing plane. Moreover, the choice of the acoustic axis angle detection plane Z1 and the acoustic axis angle detection plane Z2 does not affect the convergence and the convergence speed of the acoustic axis aligning method, which can be applied to various applications. Then, the angle of the acoustic axis is calculated by the intersection between the acoustic axis and two XOY planes with different depths. On these two planes, the computer controls the water tank system to automatically scan while acquiring the data of the sound field, and then the intersection position of the acoustic axis is judged and calculated according to the data of the sound field. Finally, the iterative algorithm combined with geometric information can converge quickly, and the performance of the aligning method does not depend on the specific mechanical implementation of adjusting an acoustic axis angle.
In addition, the criterion of intersection between the acoustic, axis and the plane can be a simple amplitude criterion or other criteria for calculating the symmetry center based on symmetry. Considering that the criterion based on amplitude (maximum amplitude) can only be used around the focusing plane, the criterion is invalid when being far away from the focusing plane (the signal amplitude on the acoustic axis is no longer the maximum value of the XOY plane). The criterion of intersection between the acoustic axis and the plane used in the present disclosure is not limited to the amplitude criterion, and other criteria can be introduced by increasing the distance between two planes to increase the robustness of measuring the acoustic axis angle.
The present disclosure will be further described in detail with reference to the drawings hereinafter.
As shown in
S1, determining three planes according to the focus position, wherein the three planes are a focusing plane Zfocus, an acoustic axis angle detection plane Z1 and an acoustic axis angle detection plane Z2; defining the distances from the probe surface to the focusing plane Zfocus, the acoustic axis angle detection plane Z1 and the acoustic axis angle detection plane Z2 as Dfocus, a and b, respectively; for convenience of description,
S2, determining a peak point of an underwater acoustic signal of the focusing plane Zfocus> and taking the peak point as an initial center point PC1, and denoting the coordinate of the initial center point PC1 as [PC_X1, PC_Y1]; wherein the peak point is the maximum point, of the underwater acoustic signal of the focusing plane;
S3, scanning the Z1 plane: in the Z1 plane, scanning with PCi as the center, changing the coordinates of X and Y, and acquiring sound field signals to obtain the intersection C1=[C1X, C1Y] between the acoustic axis and the Z1 plane through the criteria, such as centroid or maximum value;
S4, scanning the Z2 plane: in the Z2 plane, scanning with PCi as the center, changing the coordinates of X and Y, and acquiring sound field signals to obtain the intersection C2=[C2X, C1Y] between the acoustic axis and the Z2 plane through the criteria such as centroid or maximum value;
S5, calculating the inclination angle θi=[θxi, θyi] of the current acoustic axis according to the intersection points C1 and C2; wherein θxi is the included angle between the acoustic axis and Z axis after the acoustic axis is projected on XOZ plane, and θyi is the included angle between the acoustic axis and Z axis after the acoustic axis is projected on YOZ plane;
θxi=arcctan((C2X−C1X)/(b−a)) (1)
θyi=arcctan((C2Y−C1Y)/(b−a)) (2)
S6, determining an adjustment threshold according to the alignment accuracy requirement of sound field measurement, and judging whether the current inclination angle θi is smaller than the adjustment threshold; if so, aligning the acoustic axis, ending the adjustment, otherwise, proceeding to step S7;
S7, according to the inclination angle θi, calculating the next coordinate PCi=[PC_Xi, PC_Yi] by using an iterative algorithm, and the specific calculation process is as follows:
setting the inclination angle θi=[θxi, θyi] of the current acoustic axis, and the result after iteration is as follows:
if the inclination angle of the current acoustic axis is obtained after the first adjustment, calculating the adjustment amount according to the triangular relationship, the coordinate PC2=[PC_X2, PC_Y2] of a new center point is calculated according to the depth. Dfocus of the focusing plane Zfocus with the following formula:
PC_X2=PC_X1+tan(θx1)×(Dfocus−a) (3)
PC_Y2=PC_Y1+tan(θy1)×(Dfocus−a) (4)
if the inclination angle of the current acoustic axis is not obtained after the first adjustment, the calculation formula of the coordinate PC2=[PC_X2, PC_Y2] of the new center point is as follows:
PC_Xi+1=(PC_Xi−PC_Xi−1)/(θxi−θxi−1)×θxi (5)
PC_Yi+1=(PC_Yi−PC_Yi−1)/(θyi−θyi−1)×θyi (6)
in the above formulas (5) and (6), i is an integer not less than 2, the inclination angle of the current acoustic axis is θi=[θxi, θyi], and the inclination angle before the inclination angle of the current acoustic axis is adjusted is θi−1=[θxi−1, θyi−1];
S8, keeping PCi of the X axis and Y axis unchanged, moving the Z axis to the focusing plane Zfocus, manually adjusting two corner adjustment knobs of a probe clamping device until the amplitude of a hydrophone measurement signal at the current coordinate point is the maximum, and then proceeding to step S3.
Further, the position of the intersection of the acoustic axis can be judged and calculated according to the acquired sound field signal. The specific criterion can be a simple amplitude criterion or other criteria for calculating the symmetry center based on symmetry. Considering that the criterion based on amplitude (maximum amplitude) can only be used around the focusing plane, the criterion is invalid when being far away from the focusing plane (the signal amplitude on the acoustic axis is no longer the maximum value of the XOY plane). The criterion of intersection between the acoustic axis and the plane used in the present disclosure is not limited to the amplitude criterion, and other criteria can be introduced by increasing the distance between two planes to increase the robustness of measuring the acoustic axis angle. Since the acoustic axis angle measurement is a control parameter that affects the next iteration, the increase of the measurement robustness will enhance the robustness of the whole, iteration process.
Further, in S3 and S4, the computer controls the water tank system to automatically scan the current plane.
Further, the order of step S3 and the order of step S4 are interchangeable. When the order of S3 and the order of S4 are exchanged, correspondingly, S8 is performed, keeping PCi of the X axis and Y axis unchanged, moving the Z axis to the focusing plane Zfocus, manually adjusting two corner adjustment knobs of a probe clamping device until the amplitude of a hydrophone measurement signal at the current coordinate point is the maximum, and then proceeding to step S3, that is, the exchanged step S4.
To sum up, the method can quickly and effectively realize the manual acoustic axis accurate alignment.
A specific implementation process of aligning the acoustic axis using this method is given hereinafter, specifically comprising the following cycles:
S1, determining three planes Zfocus, Z1 and Z2 according to the focus position, wherein Dfocus is 13.3 cm, Z1 and Z2 are taken as the two planes at the upper side of the focus, and the corresponding distances a=2.8 cm and b=10.8 cm;
S2, finding a peak point of an underwater acoustic signal on the focusing plane Zfocus, and taking the peak point as an initial center point Pei, in which the value of the coordinate [PC_X1, PC_Y1] is [−8.18, −7.48] cm;
S3, in the Z1 plane, taking PC1 as the center, carrying out two-dimensional scanning in the range of [−5, 5] cm from the center, in X direction and in the range of [−5, 5] cm from the center in Y direction;
the specific process is as follows: setting the moving range of Y direction to move from −5 cm to 5 cm away from the center, for each value of Y, moving from −5 cm to 5 cm away from the center in X direction, recording the waveform point by point, then adding Y, and then moving in X direction of the next line; after moving in the range of Y direction, counting the maximum amplitude value of the received signal at each position point, and calculating the centroid position of the whole scanning plane according to the maximum amplitude value, in which the centroid position is C1=[C1X, C1Y]=[−9.02, −8.75] cm;
S4, in the Z2 plane, taking PC2 as the center, carrying out two-dimensional scanning in the range of [−5, 5] cm from the center in X direction and in the range of [−5, 5] cm from the center in Y direction; the specific process is as follows: setting the moving range of Y direction to move from −5 cm to 5 cm away from the center, for each value of Y, moving from −5 cm to 5 cm away from the center in X direction, recording the waveform point by point, then adding Y, and then moving in X direction of the next fine; after moving in the range of Y direction, counting the maximum amplitude value of the received signal at each position point, and calculating the centroid position of the whole scanning plane according to the maximum amplitude value, in which the centroid position is C2=[C2X, C1Y]=[−8.38, −7.78] cm;
S5, calculating the inclination angle θi=[θxi, θyi]=[4.57°, 6.91° ] of the current acoustic axis according to the intersection points C1 and C2; where θxi is the included angle between the acoustic axis and Z axis after the acoustic axis is projected on XOZ plane, sand θyi is the included angle between the acoustic axis and Z axis after the acoustic axis is projected on YOZ plane.
θxi=arcctan((C2X−C1X)/(b−a)) (1)
θyi=arcctan((C2Y−C1Y)/(b−a)) (2)
S6, selecting an adjustment threshold of 0.001 degrees according to the accuracy requirement of acoustic axis adjustment, judging whether the inclination angle θi is less than the adjustment threshold, if so, aligning the acoustic axis, ending the adjustment, otherwise, proceeding to step S7;
S7, referring to the iterative algorithm given in the above specific implementation, calculating the inclination angle θi=[θxi, θyi] of the acoustic axis and the coordinate PCi=[PC_Xi, PC_Yi];
S8, keeping PCi of the X axis and Y axis unchanged, moving the Z axis to the focusing plane Zfocus, manually adjusting two corner adjustment knobs of a probe clamping device until the amplitude of a hydrophone measurement signal at the current coordinate point is the maximum, and then proceeding to step S3.
The coordinates of θi, PCi and C1, C2 under different numbers of iterations are shown in Table 1:
As can be seen from Table 1, after three iterations, the inclination angle θi of the acoustic axis is smaller than the current adjustment threshold, 0.001 degrees, that is, the acoustic axis has been aligned.
In this embodiment, according to the measured data in the actual acoustic axis alignment process, the acoustic axis deviation angle projection on XOZ plane and YOZ plane and the results of a single adjustment amount are recorded, and
The above contents are only specific embodiments of the present disclosure to enable those skilled in the art to understand or realize the present disclosure. Many modifications to the above embodiments will be obvious 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 disclosure.
It should be understood that the present disclosure is not limited to what has been described above, and various modifications and changes can be made without departing from the scope. The scope of the present disclosure is limited only by the appended claims.
Number | Date | Country | Kind |
---|---|---|---|
202110214437.5 | Feb 2021 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2022/071128 | 1/10/2022 | WO |