Method and system for simulating acupuncture twisting by using ultrasonic phased array frequency mixing vortex focusing method

Abstract

The invention belongs to the field of ultrasonic simulation, and provides a method for simulating acupuncture and moxibustion twisting by using an ultrasonic phased array frequency mixing vortex focusing method, which comprises the following steps: constructing a three-dimensional space, taking a plane where a probe of a two-dimensional ultrasonic phased area array is located as an x-y plane, and taking a direction in which the probe of the two-dimensional ultrasonic phased area array vertically points to the depth of an ultrasonic action object as a z axis; taking a part of array elements in the two-dimensional ultrasonic phased area array as vortex array elements, and taking the rest array elements as focusing array elements; according to a preset time delay, utilizing a first preset frequency to sequentially excite the focusing array elements, so that the ultrasonic signals emitted by the focusing array elements simultaneously arrive at a preset focus to obtain a focusing sound field; and on the focusing sound field, sequentially exciting the vortex array elements according to the preset phase difference by using a second preset frequency to generate sound vortexes, and realizing simulation when the strongest energy is concentrated at the preset focus in the vortex focusing process. Energy is concentrated at the preset focus in the vortex focusing process through frequency mixing vortex focusing, and acupuncture and moxibustion twisting simulation is achieved.

Description

Technical Field

[0001] The present application belongs to the field of ultrasonic simulation, and more specifically, to a method and system for simulating acupuncture twisting by using an ultrasonic phased array mixing vortex focusing method. Background Art

[0002] Ultrasonic waves with frequencies greater than 20 kHz have strong penetration and directivity. Their unique mechanical, thermal, and biological effects make them widely used in the medical field. In the field of traditional Chinese medicine, the combination of ultrasonic technology and acupuncture therapy has produced an innovative application of traditional acupuncture - ultrasonic "acupuncture", which uses the acoustic energy generated by ultrasound waves passing through the skin and converging under the skin to stimulate specific acupuncture points, and uses acoustic energy and thermal energy as stimulation methods to make the body produce corresponding ultrasonic biological effects. Current non-invasive ultrasonic "acupuncture" mostly uses self-focusing ultrasonic transducers, which can be ultrasonic transducers made of piezoelectric materials directly into concave spherical surfaces. It does not require other auxiliary means to achieve the focusing of the output sound field, but the focusing mode is single, and there are problems such as the inability to adjust the focus position and focal spot size in real time.

[0003] In the field of ultrasound simulation, researchers have developed phased array transducers such as rectangular arrays, fan vortex arrays, and concave spherical piston arrays that can be used for ultrasound therapy, solving the problems of traditional ultrasound focusing technology. Li Yuexin et al. optimized, simulated, and verified the axial sound field generated by the 6 equal-area rings divided by the concave spherical piezoelectric ceramic sheet, and finally obtained the "needle"-shaped sound field distribution that meets the requirements of ultrasonic "acupuncture". However, in this annular ultrasonic phased array acupuncture design method, the sound field can only be focused along the central axis of the ring. Although the sound field is in the shape of a "needle", it cannot achieve the twisting of the "needle" and the flexible change of the needle insertion angle. In 2013, Yang et al. proposed a phase encoding method to form an acoustic vortex, but due to the phase singularity, the central pressure is zero. During the vortex focusing process, the energy is not concentrated at the center, and the simulation of acupuncture twisting cannot be achieved. Summary of the invention

[0004] In view of the defects of the prior art, the purpose of this application is to provide a method and system for simulating acupuncture twisting using an ultrasonic phased array mixing vortex focusing method, aiming to solve the problem that the existing ultrasonic treatment methods cannot achieve functions such as acupuncture twisting and flexible change of needle insertion angle.

[0005] To achieve the above purpose, the present application provides a method for simulating acupuncture twisting using an ultrasonic phased array frequency mixing vortex focusing method, comprising: S1: the plane where the probe of the two-dimensional ultrasonic phased array is located is taken as the xy plane, and the direction in which the probe of the two-dimensional ultrasonic phased array is perpendicular to the depth of the ultrasonic action object is taken as the z axis; and a part of the array elements in the two-dimensional ultrasonic phased array are taken as vortex array elements, and the remaining array elements are taken as focusing array elements; S2, according to the preset delay corresponding to each focusing element, sequentially exciting the focusing elements using the first preset frequency, so that the ultrasonic signals emitted by each focusing element arrive at the preset focus at the same time, so as to obtain a focused sound field; S3, on the basis of focusing the sound field, utilizes the second preset frequency to sequentially excite the vortex array elements according to the preset phase difference to generate an acoustic vortex. When the strongest energy in the vortex focusing process is concentrated at the preset focus, the function of simulating acupuncture twisting by utilizing the ultrasonic phased array mixing vortex focusing method is realized.

[0006] The present application uses a frequency mixing vortex focusing method so that during the vortex focusing process, the energy is concentrated at a preset focal position, and the vortex phase changes more regularly, thereby achieving a simulation of the acupuncture twisting treatment method.

[0007] Furthermore, the two-dimensional ultrasonic phased array is a two-dimensional square array with an array element number of m×m, the number of focusing array elements is N, and the number of vortex array elements is M, wherein: when m is an even number, the number of vortex array elements M is 4, and the remaining array elements are used as focusing array elements; when m is an odd number, the number of vortex array elements M is 9, and the remaining array elements are focusing array elements.

[0008] Furthermore, the preset phase difference is .

[0009] Furthermore, the first preset frequency is 1 MHz; the second preset frequency is 0.15 MHz to 0.84 Hz.

[0010] Furthermore, the preset delay is set to:

[0011] in, To preset the delay, t max is the acoustic wave transmission time of the focusing element farthest from the preset focus, t i For the i The acoustic wave transmission time of a focusing array element, N is the number of focusing array elements, and i=1,2,3…,N .

[0012] Furthermore, the t i for:

[0013] Where i=1,2,3…,N; R i is the distance between the focusing element and the preset focus, c is the speed of sound waves.

[0014] Furthermore, theR i for: R i =

[0015] in, x i is the central abscissa of the i-th focusing element, is the central ordinate of the i-th focusing element, x is the horizontal coordinate of the preset focus, y is the ordinate of the preset focus, and z is the vertical coordinate of the preset focus.

[0016] Further, in step S2, the sound pressure generated by all focusing array elements at a preset focus in the focused sound field is:

[0017] in, p The sound pressure at the preset focus in the focused sound field, j is an imaginary unit, f is the emission frequency of the focusing element, is the medium density, is the angular frequency of the focusing element, To preset the delay, k is the wave number, R For the i The distance between the focusing element and the preset focus, e is an exponential function, u i For the i The normal velocity of the focusing element is, s i For the i The radiation surface of the focusing element, i =1,2,3…, N .

[0018] Furthermore, in step S3, the sound pressure generated by all vortex array elements at the preset focus is:

[0019] in, is the sound pressure generated by the vortex array element at the preset focus, x is the horizontal coordinate of the preset focus, y is the ordinate of the preset focus, z is the vertical coordinate of the preset focus, t is the time it takes for the sound wave generated by the vortex array element to reach the preset focus. A The transmission distance is R The sound pressure at 0, and A=(A 0 / R 0 ) , A 0 is the amplitude of the vortex array element, is an exponential function, k is the wave number, is the angular frequency of the focusing element, j is an imaginary unit, For the m The initial phase of the focusing elements, M is the number of focusing elements, and m=1,2,3…, M .

[0020] According to another aspect of the present application, a system for simulating acupuncture twisting by using an ultrasonic phased array frequency mixing vortex focusing method is also provided for implementing any of the methods described above, the system comprising: A three-dimensional space construction module is used to construct a three-dimensional space, taking the plane where the probe of the two-dimensional ultrasonic phased array is located as the xy plane, and the direction in which the probe of the two-dimensional ultrasonic phased array is perpendicular to the depth of the ultrasonic action object as the z axis; and taking a part of the array elements in the two-dimensional ultrasonic phased array as vortex array elements, and the remaining array elements as focusing array elements; A first excitation module is used to sequentially excite the focusing elements using a first preset frequency according to a preset delay corresponding to each focusing element, so that the ultrasonic signals emitted by each focusing element arrive at a preset focus at the same time to obtain a focused sound field; The second excitation module is used to excite the vortex array elements in sequence according to a preset phase difference to generate acoustic vortices based on a focused sound field and using a second preset frequency. It is also used to determine that the strongest energy in the vortex focusing process is concentrated at the preset focus, thereby realizing the function of simulating acupuncture twisting using the ultrasonic phased array mixing vortex focusing method.

[0021] It can be understood that the beneficial effects of the second aspect mentioned above can be found in the relevant description of the first aspect mentioned above, and will not be repeated here.

[0022] In general, the above technical solutions conceived by this application have the following beneficial effects compared with the prior art: (1) Under the premise of ensuring the focusing intensity, the present application adopts the delayed focusing principle to make the sound wave signals emitted by each array element reach the preset focus at the same time and superimpose in phase at the preset focus to achieve ultrasonic focusing; then, on the basis of ultrasonic focusing, the vortex of the sound field is realized by the phase encoding method, and the energy is concentrated at the preset focus (i.e., the needle tip) during the vortex focusing process by means of frequency mixing vortex focusing, thereby realizing the simulation of acupuncture twisting.

[0023] (2) The present application solves the problem of zero central pressure and excessive initial energy in the process of forming an acoustic vortex by a phase encoding method in the prior art through a method of mixing vortex focusing. This method makes the strongest energy appear at the focus during the mixing vortex focusing process, accurately simulating the acupuncture twisting treatment method.

[0024] (3) The present application adopts a frequency-mixing vortex focusing method and uses only one phased array to achieve focusing, vortexing, and vortex focusing, thereby saving the number of array elements and simplifying the operation method.

[0025] (4) The present application applies different phase delays to the M vortex array elements on the basis of the ultrasonic focusing delay, so that the phase of the sound wave can be distributed in a spiral shape along the circumferential direction. The focusing array element can achieve central focusing, thereby realizing the simulation of acupuncture twisting. BRIEF DESCRIPTION OF THE DRAWINGS

[0026] Figure 1 It is a schematic diagram of the process of simulating acupuncture twisting based on ultrasonic phased array mixing vortex focusing method provided in the embodiment of the present application; Figure 2 It is a schematic diagram of frequency mixing vortex focusing provided by the present application; Figure 3 It is a schematic diagram of the arrangement of the two-dimensional phased array provided in this application; Figure 4 is a phase distribution diagram of the focal plane cross section when the time is 19600ns in Example 1 of the present application; Figure 5 is a curve diagram of the sound intensity at the acoustic axis passing through the focus at a time of 22000ns in Example 1 of the present application; Figure 6 is a phase distribution diagram of the focal plane cross section at a time of 19600 ns in Example 2 of the present application; Figure 7 is a curve diagram of the sound intensity at the acoustic axis passing through the focus at a time of 22000ns in Example 2 of the present application; Figure 8 It is a phase distribution diagram of the focal plane cross section when the time is 19600ns in Comparative Example 1 of the present application; Fig. 9 is a curve diagram of the sound intensity at the acoustic axis passing through the focus when the time is 22000ns in Comparative Example 1 of the present application; Fig.10 It is a phase distribution diagram of the focal plane cross section when the time is 19600ns in Comparative Example 2 of the present application; Fig.11 It is a curve diagram of the sound intensity at the acoustic axis passing through the focus at a time of 22000ns in comparative example 2 of the present application. DETAILED DESCRIPTION

[0027] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0028] The term "and / or" in this article is a description of the association relationship of associated objects, indicating that there can be three relationships. For example, A and / or B can represent: A exists alone, A and B exist at the same time, and B exists alone. The symbol " / " in this article indicates that the associated objects are in an or relationship, for example, A / B means A or B.

[0029] The terms "first" and "second" in the specification and claims herein are used to distinguish different objects rather than to describe a specific order of the objects. For example, a first response message and a second response message are used to distinguish different response messages rather than to describe a specific order of the response messages.

[0030] In the embodiments of the present application, words such as "exemplary" or "for example" are used to indicate examples, illustrations or descriptions. Any embodiment or design described as "exemplary" or "for example" in the embodiments of the present application should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of words such as "exemplary" or "for example" is intended to present related concepts in a specific way.

[0031] In the description of the embodiments of the present application, unless otherwise specified, "multiple" means two or more than two. For example, multiple processing units refer to two or more processing units, etc.; multiple elements refer to two or more elements, etc.

[0032] The present application embodiment provides a method for simulating acupuncture twisting using an ultrasonic phased array mixing vortex focusing method, such as Figure 1 As shown, the method comprises the following steps: S1 constructs a three-dimensional space, wherein the plane where the probe of the two-dimensional ultrasonic phased array is located is taken as the xy plane, and the direction in which the probe of the two-dimensional ultrasonic phased array is perpendicular to the depth of the ultrasonic action object is taken as the z axis; and a part of the array elements in the two-dimensional ultrasonic phased array are taken as vortex array elements, and the remaining array elements are taken as focusing array elements; S2, according to the preset delay corresponding to each focusing element, sequentially excites the focusing elements using the first preset frequency, so that the ultrasonic signals emitted by each focusing element arrive at the same time. Figure 2 At the preset focus (0, F) shown in the figure, a focused sound field is obtained; On the basis of focusing the sound field, S3 uses the second preset frequency to sequentially excite the vortex array elements according to the preset phase difference to generate acoustic vortices. During the vortex focusing process, the strongest energy is concentrated at the preset focus, thereby realizing the function of simulating acupuncture twisting using the ultrasonic phased array mixing vortex focusing method.

[0033] In step S1, the aforementioned two-dimensional ultrasonic phased array is a two-dimensional square array with an array element number of m×m, the number of focusing array elements is N, and the number of vortex array elements is M, wherein: when m is an even number, the number of vortex array elements M is 4, and the remaining array elements are used as focusing array elements; when m is an odd number, the number of vortex array elements M is 9, and the remaining array elements are focusing array elements.

[0034] Specifically, Figure 3 As shown, if an 8*8 phased array is used, M is 4 and N is 60; if a 9*9 phased array is used, M is 9 and N is 72; if a 7*7 phased array is used, M is 9 and N is 40; if a 6*6 phased array is used, M is 4 and N is 32. In other preferred embodiments, a larger or smaller number of phased arrays may be used, and the focusing array elements and the vortex array elements may be divided according to the array element division method provided in the present application.

[0035] In step S2, the preset phase difference is In addition, the first preset frequency is 1 MHz, and the second preset frequency is 0.15 MHz to 0.84 Hz, such as 0.15 MHz, 0.25 MHz, 0.4 MHz, 0.55 MHz, 0.6 MHz, 0.84 MHz, or any frequency value between any two of the aforementioned frequency values.

[0036] The preset delay setting of each focusing element is as follows: (1) in, To preset the delay, t max is the acoustic wave transmission time of the focusing element farthest from the preset focus, t i For the i The acoustic wave transmission time of a focusing array element, N is the number of focusing array elements, and i=1,2,3…,N .

[0037] The transmission time of the acoustic wave generated by the aforementioned focusing element after being excited is: (2) in, t i For the i The acoustic wave transmission time of each focusing element, i=1,2,3…,N; R i For the iThe distance between the focusing element and the preset focal point, c is the speed of sound waves.

[0038] The distance between the aforementioned focusing element and the preset focus is: R i = (3) in, R i For the i The distance between the focusing element and the preset focal point, x i For the i The central abscissa of the focusing element is For the i The central ordinate of the focusing element is x is the horizontal coordinate of the preset focus, y is the vertical coordinate of the preset focus, z The vertical coordinate of the preset focus.

[0039] The sound pressure generated by all the aforementioned focusing array elements at the preset focus in the focused sound field is: (4) in, p The sound pressure at the preset focus in the focused sound field, j is an imaginary unit, f is the emission frequency of the focusing element, is the angular frequency of the focusing element, To preset the delay, is the medium density, k is the wave number, e is an exponential function, u i For the i The normal velocity of the focusing element is, s i For the i The radiation surface of the focusing element, i=1,2,3…, N .

[0040] In step S3, the sound pressure generated by all vortex array elements at the preset focus is:

[0041] in, is the sound pressure generated by the vortex array element at the preset focus, x is the horizontal coordinate of the preset focus, y is the ordinate of the preset focus, z is the vertical coordinate of the preset focus, t is the time it takes for the sound wave generated by the vortex array element to reach the preset focus. AThe transmission distance is R The sound pressure at 0, and A=(A 0 / R 0 ) , A 0 is the amplitude of the vortex array element, is an exponential function, is the angular frequency of the focusing element, j is an imaginary unit, For the m The initial phase of the focusing element, m=1,2,3…, M .

[0042] Next, the technical solution provided in this application is introduced through multiple embodiments and comparative examples.

[0043] Example 1 In this embodiment, the pressure acoustic transient physical field in COMSOL is used for simulation, which mainly includes the following steps: Step 1: Use a 64-element two-dimensional phased array as the transmitting element to excite the longitudinal wave signal, with the four central elements as vortex excitation elements (i.e. vortex elements) and the remaining elements as focusing excitation elements (i.e. focusing elements).

[0044] Step 2: Excite the focusing array element at the first preset frequency of 1 MHz, and set the sine wave sound source of the focusing array element according to the following formula: F 0 (t) = P 0 sin(ωt) (6) in, F 0 (t) is the initial sine wave sound source, P 0 is the initial set pressure, is the angular frequency of the focusing element.

[0045] Based on formula (3), the distance between the focusing array element and the preset focus (i.e., the depth of the needle tip in the ultrasonic action object during acupuncture twisting) is obtained, and then based on this distance, formula (1) is used to determine the delay (i.e., the preset delay) of each array element emitting the sound wave signal. By setting the preset delay, the focusing array element far from the preset focus emits the sound wave first, and the focusing array element close to the preset focus emits the sound wave later, so that the sound wave signals emitted by each array element reach the preset focus at the same time, and are superimposed in phase at the preset focus, further maximizing the amplitude of the sound field and forming a focused sound field.

[0046] The focusing element excitation signal is as follows: F i (t) = P 0 sin[ω(t-△t)] (7) in,△t For the i Focusing elements at any time t Sine wave sound source F i (t) Relative to F 0 (t) Trigger delay time, i is the serial number of the focusing element.

[0047] Step 3: Based on the ultrasound focusing, according to the preset phase difference , the vortex array element is excited at a second preset frequency of 0.84 MHz by means of spiral phase encoding to form an acoustic vortex. The vortex array element excitation signal is as follows: S n (t) = F n (t)exp(-iωt)exp(-iΦ n ) (8) in, F n (t) is the sinusoidal sound source of the nth vortex array element, Φ n is the phase difference of the nth vortex array element, Φ n = 2 π(n-1)m / N ; m is the topological charge number, which is 1 here; N is the number of vortex array elements, n is the serial number of the vortex array element.

[0048] The phase distribution of the focal plane section and the sound intensity curve at the focal position are collected. Figure 4 is the phase distribution diagram of the focal plane section at 19600ns, Figure 5 This is the sound intensity curve at the focal acoustic axis at 22000ns.

[0049] like Figure 4 As shown, the phase distribution of the mixed frequency vortex focusing under the excitation of the first preset frequency of 1MHz and the second preset frequency of 0.84MHz is more regular than that of the single frequency vortex focusing, but the change of the phase from -π to π is more chaotic, making the simulation effect of acupuncture twisting poor. Figure 5 As shown in Figure 2, the mixing vortex focusing under 1MHz and 0.84MHz excitation makes the energy concentrated at the tip of the needle, and the sound intensity at the focusing position is 9.565W / cm 2Therefore, it is concluded that the simulation of acupuncture twisting can be achieved by mixing vortex focusing under the excitation frequencies of 1 MHz and 0.84 MHz, but the effect is poor.

[0050] Example 2 In this embodiment, the pressure acoustic transient physical field in COMSOL is used for simulation, and the main steps are different from those in Embodiment 1: the vortex array element is excited at a second preset frequency of 0.15 MHz to form an acoustic vortex.

[0051] The phase distribution of the focal plane section and the sound intensity curve at the focal position are collected. Figure 6 is the phase distribution diagram of the focal plane section at 19600ns, Fig.10 is the sound intensity curve at the focal sound axis at 22000ns. Figure 6 As shown, the phase distribution of the mixing vortex focusing under the excitation of 1 MHz and 0.15 MHz is more regular than that in the above-mentioned embodiment 1, and the phase changes from -π to π significantly, which can realize the simulation of acupuncture twisting.

[0052] like Figure 7 As shown in Figure 2, the mixing vortex focusing under the excitation of 1 MHz and 0.15 MHz respectively makes the energy concentrated at the tip of the needle, and the sound intensity at the focusing position is 9.586 W / cm 2 It is proved that the frequency mixing vortex focusing method under 1MHz and 0.15MHz excitation can better realize the simulation of acupuncture twisting.

[0053] Comparative Example 1 The difference between this comparative example 1 and embodiment 1 is that the first preset frequency for exciting the focusing element is the same as the second preset frequency for exciting the vortex element, both of which are 1 MHz. After simulation according to the steps in embodiment 1, the phase distribution of the focal plane section and the sound intensity curve at the preset focus position are collected.

[0054] Figure 8 is the phase distribution diagram of the focal plane section at 19600ns, Fig. 9 is the sound intensity curve at the focal sound axis at 22000ns. Figure 8 As shown in the figure, the single-frequency vortex focusing under 1MHz excitation cannot form a regular phase spiral, and the phase change from -π to π cannot be clearly observed in the phase diagram. Therefore, it is actually impossible to simulate acupuncture twisting treatment. Fig. 9 As shown in Figure 2, the single-frequency vortex focusing under 1MHz excitation makes the energy concentrated at the central axis, and the sound intensity at the focusing position is 9.545W / cm 2From this, it can be seen that the mixing vortex focusing method makes the energy concentrated at the center during the vortex focusing process, and the change of the vortex phase cannot be clearly observed, and the simulation of acupuncture twisting cannot be achieved.

[0055] Comparative Example 2 The difference between Comparative Example 2 and Example 1 is that in Comparative Example 2, all array elements of the two-dimensional phased array are used as focusing array elements, and some focusing array elements located at the center of the two-dimensional phased array are excited by the second preset frequency of 0.15 MHz, and the remaining focusing array elements located at the edge of the two-dimensional phased array are excited by the first preset frequency of 1 MHz.

[0056] Collect the phase distribution of the focal plane section and the energy at the focal position, Fig.10 is the phase distribution diagram of the focal plane section at 19600ns, Fig.11 is the sound intensity curve at the focal sound axis at 22000ns. Fig.10 As shown in Figure 1, frequency mixing focusing under 1 and 0.15 MHz excitation cannot form phase spirals and cannot simulate acupuncture twisting. Fig.11 As shown in Figure 2, the frequency mixing focusing under 1 and 0.15 MHz excitation can concentrate the energy at the center. From the sound intensity curve, it can be seen that the sound intensity at the focus position is 9.240 W / cm 2 The frequency mixing focusing method under 1MHz and 0.15MHz excitation can only achieve the simulation of acupuncture lifting and insertion treatment, but cannot achieve the simulation of acupuncture twisting.

[0057] In summary, the method of simulating acupuncture twisting based on ultrasonic phased array mixing vortex focusing method provided in this application has an energy difference of only 0.346 W / cm2 compared with simple focusing. 2 Under the premise of this, the simulation of acupuncture twisting was achieved through the method of mixed vortex focusing. During the simulation process, the vortex phase changes were relatively regular, and the maximum energy was concentrated at the needle tip.

[0058] It should be understood that the above-mentioned device is used to execute the method in the above-mentioned embodiment. The implementation principle and technical effect of the corresponding program module in the device are similar to those described in the above-mentioned method. The working process of the device can refer to the corresponding process in the above-mentioned method, which will not be repeated here.

[0059] Based on the method in the above embodiment, an embodiment of the present application provides a computer-readable storage medium, which stores a computer program. When the computer program runs on a processor, the processor executes the method in the above embodiment.

[0060] Based on the method in the above embodiment, an embodiment of the present application provides a computer program product. When the computer program product runs on a processor, the processor executes the method in the above embodiment.

[0061] It is understandable that the processor in the embodiment of the present application may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application specific integrated circuits (ASIC), field programmable gate arrays (FPGA) or other programmable logic devices, transistor logic devices, hardware components or any combination thereof. The general-purpose processor may be a microprocessor or any conventional processor.

[0062] The method steps in the embodiments of the present application can be implemented by hardware or by a processor executing software instructions. The software instructions can be composed of corresponding software modules, and the software modules can be stored in random access memory (RAM), flash memory, read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), electrically erasable programmable read-only memory (EEPROM), registers, hard disks, mobile hard disks, CD-ROMs, or any other form of storage medium known in the art. An exemplary storage medium is coupled to the processor so that the processor can read information from the storage medium and write information to the storage medium. Of course, the storage medium can also be a component of the processor. The processor and the storage medium can be located in an ASIC.

[0063] In the above embodiments, it can be implemented in whole or in part by software, hardware, firmware or any combination thereof. When implemented by software, it can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions. When the computer program instructions are loaded and executed on a computer, the process or function described in the embodiment of the present application is generated in whole or in part. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions may be stored in a computer-readable storage medium or transmitted through the computer-readable storage medium. The computer instructions may be transmitted from a website site, computer, server or data center to another website site, computer, server or data center by wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.). The computer-readable storage medium may be any available medium that a computer can access or a data storage device such as a server or data center that includes one or more available media integrated. The available medium may be a magnetic medium (e.g., a floppy disk, a hard disk, a tape), an optical medium (e.g., a DVD), or a semiconductor medium (e.g., a solid state disk (SSD)), etc.

[0064] It should be understood that the various numerical numbers involved in the embodiments of the present application are only used for the convenience of description and are not used to limit the scope of the embodiments of the present application.

[0065] It will be easily understood by those skilled in the art that the above description is only a preferred embodiment of the present application and is not intended to limit the present application. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present application shall be included in the scope of protection of the present application.

Claims

1. A method for simulating acupuncture twisting by using an ultrasonic phased array mixing vortex focusing method, characterized in that: include: S1: the plane where the probe of the two-dimensional ultrasonic phased array is located is taken as the xy plane, and the direction in which the probe of the two-dimensional ultrasonic phased array is perpendicular to the depth of the ultrasonic action object is taken as the z axis; and a part of the array elements in the two-dimensional ultrasonic phased array are taken as vortex array elements, and the remaining array elements are taken as focusing array elements; S2, according to the preset delay corresponding to each focusing element, sequentially exciting the focusing elements using the first preset frequency, so that the ultrasonic signals emitted by each focusing element arrive at the preset focus at the same time, so as to obtain a focused sound field; S3, on the basis of focusing the sound field, utilizes the second preset frequency to sequentially excite the vortex array elements according to the preset phase difference to generate an acoustic vortex. When the strongest sound wave is concentrated at the preset focus during the vortex focusing process, the function of simulating acupuncture twisting by utilizing the ultrasonic phased array mixing vortex focusing method is realized.

2. The method for simulating acupuncture twisting by using ultrasonic phased array frequency mixing vortex focusing method as claimed in claim 1, characterized in that: The two-dimensional ultrasonic phased array is a two-dimensional square array with an array element number of m×m, the number of focusing array elements is N, and the number of vortex array elements is M, wherein: when m is an even number, the number of vortex array elements M is 4, and the remaining array elements are used as focusing array elements; when m is an odd number, the number of vortex array elements M is 9, and the remaining array elements are focusing array elements.

3. The method for simulating acupuncture twisting by using ultrasonic phased array frequency mixing vortex focusing method as claimed in claim 2, characterized in that: The preset phase difference is .

4. The method for simulating acupuncture twisting by using ultrasonic phased array frequency mixing vortex focusing method as claimed in claim 1, characterized in that: The first preset frequency is 1 MHz; the second preset frequency is 0.15 MHz to 0.84 Hz.

5. The method for simulating acupuncture twisting by using ultrasonic phased array frequency mixing vortex focusing method as claimed in claim 1, characterized in that: The preset delay is set to: in, To preset the delay, t max is the acoustic wave transmission time of the focusing element farthest from the preset focus, t i For the i The acoustic wave transmission time of a focusing array element, N is the number of focusing array elements, and i=1,2,3…,N .

6. The method for simulating acupuncture twisting by using ultrasonic phased array frequency mixing vortex focusing method as claimed in claim 5, characterized in that: Said t i for: in, R i For the i The distance between the focusing element and the preset focal point, c is the speed of sound waves.

7. The method for simulating acupuncture twisting by using ultrasonic phased array frequency mixing vortex focusing method as claimed in claim 6, characterized in that: Said R i for: R i = in, x i For the i The central abscissa of the focusing element is For the i The central ordinate of the focusing element, x is the horizontal coordinate of the preset focus, y is the vertical coordinate of the preset focus, z The vertical coordinate of the preset focus.

8. The method for simulating acupuncture twisting by using ultrasonic phased array frequency mixing vortex focusing method as claimed in claim 1, characterized in that: In step S2, the sound pressure generated by all focusing array elements at the preset focus in the focused sound field is: in, p The sound pressure at the preset focus in the focused sound field, j is an imaginary unit, f is the emission frequency of the focusing element, is the medium density, is the angular frequency of the focusing element, To preset the delay, k is the wave number, R i For the i The distance between the focusing element and the preset focal point, e is an exponential function, u i For the i The normal velocity of the focusing element is, s i For the i The radiation surface of the focusing element, i=1, 2,3…, N .

9. The method of simulating acupuncture twisting using ultrasonic phased array mixed frequency vortex focusing method as claimed in claim 1, characterized in that, In step S3, the sound pressure generated by all vortex array elements at the preset focus is: in, It is the sound pressure generated by the vortex array element at the preset focus. x is the horizontal coordinate of the preset focus, y It is the vertical coordinate of the preset focus, z is the vertical coordinate of the preset focus, t The time when the sound wave generated by the vortex array element is transmitted to the preset focus, A The transmission distance is R The sound pressure at 0, and A=(A 0 / R 0 ) , k is the wave number, A 0 is the amplitude of the vortex array element, is an exponential function, is the angular frequency of the focusing element, j is an imaginary unit, For the m The initial phase of the focusing array element, M is the number of the focusing array element, and m=1,2,3…, M .

10. A system for simulating acupuncture twisting by using an ultrasonic phased array frequency mixing vortex focusing method for implementing the method as claimed in any one of claims 1 to 9, characterized in that: The system comprises: A three-dimensional space construction module is used to construct a three-dimensional space, taking the plane where the probe of the two-dimensional ultrasonic phased array is located as the xy plane, and the direction in which the probe of the two-dimensional ultrasonic phased array is perpendicular to the depth of the ultrasonic action object as the z axis; and taking a part of the array elements in the two-dimensional ultrasonic phased array as vortex array elements, and the remaining array elements as focusing array elements; A first excitation module is used to sequentially excite the focusing elements using a first preset frequency according to a preset delay corresponding to each focusing element, so that the ultrasonic signals emitted by each focusing element arrive at a preset focus at the same time to obtain a focused sound field; The second excitation module is used to excite the vortex array elements in sequence according to a preset phase difference to generate an acoustic vortex based on a focused sound field and using a second preset frequency. It is also used to determine whether the strongest energy in the vortex focusing process is concentrated at the preset focus, thereby realizing the function of simulating acupuncture twisting using an ultrasonic phased array mixing vortex focusing method.

Images