CN101692072A - Calculation method of sound field of circular arc-shaped linear ultrasonic phased array - Google Patents
Calculation method of sound field of circular arc-shaped linear ultrasonic phased array Download PDFInfo
- Publication number
- CN101692072A CN101692072A CN200910148336A CN200910148336A CN101692072A CN 101692072 A CN101692072 A CN 101692072A CN 200910148336 A CN200910148336 A CN 200910148336A CN 200910148336 A CN200910148336 A CN 200910148336A CN 101692072 A CN101692072 A CN 101692072A
- Authority
- CN
- China
- Prior art keywords
- array element
- array
- sound field
- circular arc
- transducer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
The invention provides a calculation method of the sound field of a circular arc-shaped linear ultrasonic phased array, comprising a beam control focusing algorithm, a sound field calculation method and the influence of structural parameters on sound field characteristics. The calculation method uses an angle method to show the position information of an energy converter array element so as to deduct the beam control focusing algorithm of a circular arc shape linear phased array energy converter; a circular array element is similar to a rectangular array element and is combined with coordinate conversion to calculate the sound field of the circular arc-shaped linear phased array; the method is performed with simulation verification respectively by Reily integration and non-paraxial approximation multi-Gaussian models, and the obtained result is very identical. In addition, the calculation method analyzes the influence of structural parameters on sound field characteristics and provides theoretical basis for designing and optimizing the energy converter.
Description
One, technical field
The present invention relates to is a kind of calculation method of sound field of circular arc-shaped linear ultrasonic phased array, and this method provides theoretical foundation for the array energy transducer of design and research special shape.
Two, technical background
In Ultrasonic Detection, generally adopt contact measurement method.When the testee surface configuration is on-plane surface, if adopt linear phase controlled array transducer to detect,, contact can not reduce detection sensitivity because matching, can't accurately locate defective.Under these circumstances, be necessary to develop the array energy transducer of special shape to satisfy complicated detection requirement.Because computer technology rapid development, modeling analysis and simulation calculation provide very effective method for developing and analyze transducer performance, can reflect the sound field characteristic of transducer more intuitively, and optimize the structural parameters of transducer.
The circular arc-shaped linear phase array transducer is a kind of curved surface linear phase controlled array transducer simple in structure, and each array element has identical radius-of-curvature, and the positional information of array element is represented easily.And the linear phase controlled array transducer of irregular surface can be approximated to be by a plurality of circular array unit and forms, and different is that circular arc array element has different radius-of-curvature.Therefore, the circular arc-shaped linear phase array transducer is the basis of research and design complicated shape phase array transducer, has important researching value.
Three, summary of the invention
The objective of the invention is in order to overcome when the object of detection of complex surface configuration, the deficiency that linear phase controlled array transducer exists, and a kind of calculation method of sound field of regular surfaces linear phase controlled array transducer is provided, for design with optimize transducer the theoretical reference foundation is provided.
The object of the present invention is achieved like this: it comprises the influence to the phased array sound field of focusing algorithm, calculation method of sound field and the array parameter of the control of circular arc-shaped linear ultrasonic phased array wave beam; The focusing algorithm of described circular arc-shaped linear ultrasonic phased array wave beam control adopts preset angle configuration to represent the positional information of array element, thereby derives the focusing law of wave beam control; Described calculation method of sound field adopts circular arc array element to be approximately the method that rectangle array element and coordinate transform combine, and at first calculates the sound field of single circular arc array element, and then according to Huygens' principle, the sound field of phase array transducer equals the stack of a plurality of array element sound fields; Described array parameter is mainly studied the influence to sound field of array number, array element width and array element distance to the influence of phased array sound field.
Four, description of drawings
Fig. 1 is a circular arc-shaped linear phase array transducer structural representation;
Fig. 2 is the result of calculation in different spatial linear array and circular arc-shaped linear array time delay;
Fig. 3 is a coordinate system transformation;
Fig. 4 is the geometric configuration of single circular arc array element;
Fig. 5 is two-dimensional acoustic field result of calculation (F=20, θ
s=0 °);
Fig. 6 is two-dimensional acoustic field result of calculation (F=20, θ
s=5 °);
Fig. 7 is the asynchronous focusing acoustic field figure of array number;
Fig. 8 is the asynchronous transducer acoustic field of array element width;
Fig. 9 is a transducer deflection focusing sound field;
Figure 10 is the influence of array element height to sound field.
Five, embodiment
For a more detailed description to the present invention for example below in conjunction with accompanying drawing:
In conjunction with Fig. 1, suppose that the circular arc-shaped linear phased array is formed by the linear phase controlled array distortion, according to isometric principle, the aperture D of linear array equals the arc length L of circular arc-shaped linear array
D=L=θ
c/180×π×R (1)
Wherein R is a radius-of-curvature, θ
cBe the pairing central angle of the arc length of circular arc-shaped linear array.
In conjunction with Fig. 1, the spacing of the width of circular arc array element and adjacent two array elements is used symbol a respectively
c, d
cExpression, their pairing central angles are defined as θ respectively
aAnd θ
d, the relational expression below satisfying
(N-1)θ
d+θ
a=θ
c (2)
Wherein N is an array number.I the pairing central angle in array element center is
Under rectangular coordinate system, the centre coordinate of i array element is: x
i=Rsin θ
i, z
i=R-Rcos θ
i
In conjunction with Fig. 1, according to top condition, in the xoz plane, the computing formula of each array element time delay is
τ wherein
iBe the time delay of i array element, F is a focal length, θ
sBe deflection angle, c is the speed of ultrasound wave in medium, θ
iIt is the angle between i array element center and the z axle.
When F>>during R, τ
i≈ F/c+t.
In conjunction with Fig. 1, the circular arc-shaped linear phase array transducer of forming with 8 circular array units is an example, according to top research method, works as θ
c=60 °, θ
a=4 °, θ
dIn the time of=8 °, the locus is calculated the time delay of linear phase controlled array and each array element of circular arc-shaped linear phased array respectively not simultaneously.
In conjunction with Fig. 2, when the shape of array element and position change, be different the time delay of each array element that calculates at the space point, and the focusing effect that is produced also is different; When the testee surface configuration changes, if still adopt traditional linear array transducer to detect, will influence actual focusing effect, therefore should study new focusing delay algorithm to adapt to different detection requirements according to the situation on testee surface.
Because each array element of circular arc-shaped linear phased array is in same plane, in order to study conveniently, when calculating the acoustic pressure that each array element produces at spatial point P place, the employing coordinate transform is transformed into a P under the new coordinate system and calculates.
In conjunction with Fig. 3, after translation and rotational transform, the new coordinate system (ξ of foundation
iO ' η
i).Wherein true origin o ' is i array element center, η
iAxle is the central axis of i array element.According to geometric relationship, some P conversion formula under new coordinate system is
(x wherein
i, z
i) be the centre coordinate of i array element at former coordinate system, (ξ
Pi, η
Pi) be that some P is at (ξ
iO ' η
i) coordinate under the coordinate system.
In conjunction with Fig. 4, single array element according to the chord length formula, obtains the width a of array element through after the coordinate transform
cWith array element distance d
cFor
a
c=2Rsin(θ
a/2),d
c=2Rsin(θ
d/2)(6)
Work as θ
aWhen very little, sin (θ
a/ 2)~θ
a/ 2, a
c~R θ
a, that is: single circular arc array element can be approximated to be single rectangle array element.
With a radius of curvature R=20mm, central angle θ
c=60 ° of circular arc-shaped linear phased arrays are example, as the centre frequency f of transducer
0=5MHz, array number N=8, θ
a=4 °, θ
d=7 °, focal length F=20mm, deflection angle is respectively θ
s=5 ° and θ
sThe two-dimensional acoustic field figure of the phase array transducer that=0 ° of employing rayleigh integral (RSI) and non-paraxial approximate (NMG) Model Calculation obtain.
In conjunction with Fig. 5 and Fig. 6, the two-dimensional acoustic field result of calculation that adopts RSI and two kinds of models of NMG to obtain is very identical.
In order to study circular arc-shaped linear phased array elements spacing d
cTo the influence of sound field characteristic, get array element width a
cBe fixed value, i.e. θ
aConstant, by changing array number N, observe the influence of array element distance to acoustic pressure and beam feature.
In conjunction with Fig. 7, with θ
a=3 ° is example, in the xoz plane, when wave beam in (0,20) when locating to focus on, along with the increase of array number N, the acoustic pressure of focusing place increases gradually, main lobe width narrows down, graing lobe is eliminated gradually.Therefore can eliminate graing lobe by increasing the method for array number, strengthen the energy of main lobe.
In conjunction with Fig. 8,, be θ by increasing array element width when array number fixedly the time
a, observe the influence of array element distance and array element width to beam feature and acoustic pressure.When N=10, θ
aWhen getting different values, satisfying θ
a<θ
dSituation under, increase θ
aNot only can increase acoustic pressure but also can also change beam feature.
In conjunction with Fig. 9, when wave beam during in (5,0,20) point focusing, the sound field energy has obtained enhancing after the wave beam deflection focusing; Increase array element width, the sound field energy also is enhanced.
By last surface analysis, work as θ
dIn the time of<6 °, can eliminate graing lobe substantially; Under the certain clear condition of array number, can increase the energy that array element width improves sound field.
According to top analysis result, the circular arc-shaped linear phase array transducer of forming with 10 array elements is an example, changes the value of array element height b, observes the change of sound field situation.
In conjunction with Figure 10, when b not simultaneously, adopt the two-dimensional acoustic field of this transducer of NMG Model Calculation; The array element height is also not obvious to the influence of sound field, can consider to select suitable value from the processing and fabricating aspect.
Claims (7)
1. circular arc-shaped linear ultrasonic phased array transducer acoustic field computing method, it comprises the influence to sound field of beam control focusing algorithm, calculation method of sound field and structural parameters.It is characterized in that: the circular arc-shaped linear ultrasonic phased array transducer is formed by a plurality of circular array are first, the position of each array element is not in same plane, in order to derive the focusing algorithm of wave beam control, must adopt new method to represent the position of transducer array element; According to Huygens' principle, when calculating the sound field of circular arc-shaped linear phase array transducer, at first calculate the sound field of single circular arc array element, superpose after postponing accordingly then; The analytical structure parameter comprises the influence to sound field of array element width, array element distance, array number and array element height.
2. method according to claim 1 is characterized in that, the circular arc-shaped linear ultrasonic phased array transducer is formed by the linear ultrasonic phased array transducer distortion.
3. method according to claim 1 is characterized in that: adopt preset angle configuration to represent the position of array element in the array energy transducer, promptly represent with the radius-of-curvature of circular arc-shaped linear phase array transducer and the angle of array element center and z axle.
4. method according to claim 1 is characterized in that: circular arc array element is approximately rectangle array element, and calculates the sound field of single array element in conjunction with coordinate transformation method in newly-established coordinate system.
5. method according to claim 1 is characterized in that: when calculating the sound field of single circular arc array element, each circular arc array element is approximately rectangle array element, and in conjunction with coordinate transformation method, calculates the distribution situation of the whole sound field of transducer.
6. method according to claim 4 is characterized in that: the center with each array element is an initial point, and the central axis of array element is set up new coordinate system for the z axle, with the space arbitrarily any coordinate conversion in new coordinate system, calculate.
7. method according to claim 1 is characterized in that: when analyzing array element width and array element distance to the influencing of circular arc-shaped linear phase array transducer sound field, mainly study the influence of these two pairing central angles of variable to sound field.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910148336A CN101692072A (en) | 2009-06-16 | 2009-06-16 | Calculation method of sound field of circular arc-shaped linear ultrasonic phased array |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910148336A CN101692072A (en) | 2009-06-16 | 2009-06-16 | Calculation method of sound field of circular arc-shaped linear ultrasonic phased array |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101692072A true CN101692072A (en) | 2010-04-07 |
Family
ID=42080780
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910148336A Pending CN101692072A (en) | 2009-06-16 | 2009-06-16 | Calculation method of sound field of circular arc-shaped linear ultrasonic phased array |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101692072A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102824190A (en) * | 2012-09-24 | 2012-12-19 | 深圳大学 | Two-dimensional annular phased array ultrasonic transducer structure |
CN102967659A (en) * | 2012-10-30 | 2013-03-13 | 广东电网公司电力科学研究院 | Calculation method of sound field distribution of phased array ultrasonic probe in flaw detection in multilayer medium |
CN106501367A (en) * | 2016-10-10 | 2017-03-15 | 清华大学 | Phased array supersonic echo-wave imaging method based on elliptic arc scan transformation |
CN107271549A (en) * | 2013-11-07 | 2017-10-20 | 三菱日立电力***株式会社 | The control device and ultrasonic fault detection system of ultrasonic fault detection system |
CN107657078A (en) * | 2017-08-28 | 2018-02-02 | 天津大学 | Ultrasonic phase array floating-point focus emission implementation method based on FPGA |
-
2009
- 2009-06-16 CN CN200910148336A patent/CN101692072A/en active Pending
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102824190A (en) * | 2012-09-24 | 2012-12-19 | 深圳大学 | Two-dimensional annular phased array ultrasonic transducer structure |
CN102824190B (en) * | 2012-09-24 | 2015-02-04 | 深圳大学 | Two-dimensional annular phased array ultrasonic transducer structure |
CN102967659A (en) * | 2012-10-30 | 2013-03-13 | 广东电网公司电力科学研究院 | Calculation method of sound field distribution of phased array ultrasonic probe in flaw detection in multilayer medium |
CN102967659B (en) * | 2012-10-30 | 2014-11-26 | 广东电网公司电力科学研究院 | Calculation method of sound field distribution of phased array ultrasonic probe in flaw detection in multilayer medium |
CN107271549A (en) * | 2013-11-07 | 2017-10-20 | 三菱日立电力***株式会社 | The control device and ultrasonic fault detection system of ultrasonic fault detection system |
CN107271549B (en) * | 2013-11-07 | 2019-11-22 | 三菱日立电力***株式会社 | The control device and ultrasonic fault detection system of ultrasonic fault detection system |
CN106501367A (en) * | 2016-10-10 | 2017-03-15 | 清华大学 | Phased array supersonic echo-wave imaging method based on elliptic arc scan transformation |
CN106501367B (en) * | 2016-10-10 | 2019-02-22 | 清华大学 | Implementation method is imaged in phased array supersonic based on elliptic arc scan transformation |
CN107657078A (en) * | 2017-08-28 | 2018-02-02 | 天津大学 | Ultrasonic phase array floating-point focus emission implementation method based on FPGA |
CN107657078B (en) * | 2017-08-28 | 2021-01-19 | 天津大学 | Ultrasonic phased array floating point focusing transmission implementation method based on FPGA |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104730491B (en) | A kind of virtual array DOA estimation method based on L-type battle array | |
CN101692072A (en) | Calculation method of sound field of circular arc-shaped linear ultrasonic phased array | |
EP3261361A1 (en) | Beamforming method based on arrays of microphones and corresponding apparatus | |
CN103529364B (en) | Many Partial Discharge Sources localization method in oil based on ultrasonic thinned array sensor | |
CN101995574B (en) | Near field focusing beam forming positioning method | |
CN103999151B (en) | In calculating, effective wideband filtered and addition array focus on | |
CN101644773B (en) | Real-time frequency domain super-resolution direction estimation method and device | |
CN103743469B (en) | Piezoelectric underwater sound transducer acoustic radiation mode measurement method and system | |
CN102043145B (en) | Rapid broadband frequency domain beamforming method based on acoustic vector sensor uniform linear array | |
CN104931919A (en) | Direction finding method for two dimensional plane annular array | |
CN102736063A (en) | Near-field sound source positioning method | |
CN104730146A (en) | Ultrasonic array composite imaging method for quantitative evaluation of plate structure cracks | |
CN104678359B (en) | A kind of porous sound holographic method of sound field identification | |
CN106855619A (en) | A kind of method of the resolution ratio of acquisition MIMO imaging radar system all directions | |
CN103344706B (en) | Design method for linear array phased array probe | |
CN110736976B (en) | Method for estimating performance of sonar beam former of any array | |
CN110376546B (en) | Far field and near field mixed information source positioning method based on covariance matrix difference | |
Liu et al. | Application and challenges of signal processing techniques for Lamb waves structural integrity evaluation: Part B—Defects imaging and recognition techniques | |
CN105572629A (en) | Two-dimensional direction measuring method of low operation complexity and applicable to any array structure | |
CN101515033A (en) | Multilayer stereoscopic grid array for recognizing noise source in a beam shaping method | |
Øygard et al. | Accurate prediction of transmission through a lensed row-column addressed array | |
CN106124044A (en) | Medicine ball identification of sound source low sidelobe ultrahigh resolution acoustic picture fast acquiring method | |
CN101672824A (en) | Method for calculating sound field of irregular surface linear phase controlled array transducer | |
CN104809282A (en) | Method for constructing shaping field source based on metal scattering body assisting | |
CN112255429B (en) | Three-dimensional wind parameter measuring method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20100407 |