The present invention is comprised in the field of material analysis, particularly in analysis methods using signal processing.
Specifically, the invention relates to methods which allow obtaining data relating to the elasticity of materials.
Torsional waves are a spatial distribution of transverse waves propagating along an axis in which particle movement occurs along a circumference the center of which is said axis, such that the amplitude of the movement in the generating plane is proportional to the distance from the axis within the diameter of the transducer.
These waves propagate through solid and semi-solid media, but not through perfect liquids, so measuring the speed of sound in media of this type can be very useful for studying their structural characteristics.
A transducer is a device capable of transforming or converting a specific type of input energy into another different type of output energy. These devices include, among others, electromechanical transducers which transform electrical energy into mechanical energy in the form of bidirectional displacements elastically coupled with stresses.
Ultrasonic transducers emit and receive sonic and/or ultrasonic waves, which allows, based on solid mechanics, identifying changes in tissue consistency which may indicate the presence of tumors, and quantifying mechanical or physical changes in the tissue can anticipate certain pathologies sooner than other diagnostic techniques can. Furthermore, the operating mode of the torsion sensor for obtaining non-linear parameters provides description and itemization in physical terms, rendering a value that is linked to the mode in which the fibers and the matrix in which they are embedded undergo deformation.
The problem with quasi-compressible materials (soft tissues and gels), the Poisson's ratio of which is about 0.5, lies in their bulk modulus and shear modulus being different. P- and S-type waves with different magnitudes propagate in these materials, and spurious P-waves that predominate and mask the S-waves are generated, not allowing commercial devices to read the S-waves which provide information about the shear modulus. If the operating mode is used as a diagnostic mode differentiating a healthy tissue from a diseased tissue at the level of tissue fibers and their support matrices, this mode has no ionizing effects which are indeed found in other diagnostic means such as X-rays.
The propagation of torsional waves is correlated with shear modulus by means of the elastic wave propagation equations, whereas the propagation of longitudinal waves is correlated with bulk modulus. In soft tissues, the non-linearity parameters vary by several orders of magnitude, so by using ultrasonic transducers based on non-linear torsional waves, a sensitivity that greater than that obtained with ultrasonic transducers based on P- and S-waves can be achieved.
Non-linear torsional wave generators were unknown up until the invention of the device, but the origin of devices for obtaining non-linear parameters based on P- and S-waves can be traced to a thesis published in 2009 by Muir entitled One-Sided Ultrasonic Determination of Third Order Elastic Constants using Angle-Beam Acoustoelasticity Measurements. The main limitation of this technique is that the P-wave transducers propagate at an angle so that, by means of mode conversion, an S-wave, the non-linear parameters of which are analyzed, is generated. It is extremely difficult to transfer this method to the field of tissues and quasi-fluids and almost impossible to extrapolate it to in-vivo assays. This has only been studied up until now in metals or homogenous materials.
Techniques such as those described in [Cristian Pantea, Curtis F Osterhoudt, and Dipen N Sinha. Determination of acoustical nonlinear parameter β of water using the finite amplitude method. Ultrasonics, 53(5):1012-1019, 2013] or [Pham Chi Vinh and Jos{tilde over (e)} Merodio. On acoustoelasticity and the elastic constants of soft biological tissues. Journal of Mechanics of Materials and Structures, 8(5): 359-367, 2013] are also known to be useful for obtaining parameters for measuring acoustic non-linearity in water by means of a mixture of waves and measuring non-linearity in tissue by means of dynamic acoustoelastic testing, DAET, in tissue, respectively.
These techniques cannot be readily extrapolated as a diagnostic method due to experimental difficulty and due to the fact that the non-linear parameters have never been separated depending on their physical and biological origin before.
There is therefore a need for a method for obtaining data relating to the elasticity or consistency of materials which allows differentiating the volumetric part (linked to deformations due to tensile and compression forces) from the offset part (linked to deformations due to shear forces) in the studied sample.
The present invention relates to a method or operating mode, hereinafter “method of the invention” which, using a device capable of emitting and receiving sonic and/or ultrasonic torsional waves, allows obtaining data relating to the consistency or elasticity of quasi-incompressible solid media (with Poisson's ratio close to 0.5), preferably biological tissues or quasi-fluids, based on the separation of non-linear parameters.
In particular, the method of the invention allows identifying changes in the consistency or elasticity of the analyzed sample depending on the behavior of the shear modulus thereof, by means of analyzing the torsional waves propagated through the material and received by the device of the invention.
The method of the invention uses a method for generating and measuring ultrasounds by means of the non-conventional use of shear waves and/or surface instead of longitudinal waves as they are several orders of magnitude more sensitive to variations in the microstructure of the relevant material which are closely linked to the non-linear parameters of the material.
The method of the invention improves the current technique, making it possible to differentiate the volumetric part (linked to deformations due to tensile and compression forces) from the offset part (linked to deformations due to shear forces) in the studied sample so as to allow determining its elasticity or consistency characteristics.
The use of non-linear sonic or ultrasonic waves as a physical magnitude has other fundamental advantages. Firstly, it is a mechanical wave that can be controlled, and is therefore more sensitive to mechanical properties than any other indirect measurement. Secondly, the wave is generated in a low-energy regimen which is more sensitive to variations in non-linear tissue parameters than waves generated at a high-energy regimen are.
Another object of the invention is a system, hereinafter “system of the invention,” comprising the means required for carrying out the method of the invention.
Another object of the invention is a computer program comprising instructions for making a computer carry out the method of the invention. Likewise, other objects of the invention are a computer-readable storage medium comprising program instructions capable of making a computer carry out the method of the invention and a transmissible signal comprising program instructions capable of making a computer carry out the method of the invention.
By way of example, it can be obtained using the method of the invention parameters that are useful for identifying, based on solids mechanics, changes in non-linear parameters that define behavioral and state changes in materials manifested through changes in their elastic parameters which in turn control the propagation of waves going through them.
By means of analyzing these changes in elastic parameters, the presence and the type of tumors and any disorder manifesting as changes in said elastic parameters can be identified.
Throughout the present invention, “wave train” will be understood as a set of two or more identical sinusoidal wave cycles. In other words, two or more cycles having the same sinusoidal wave will be emitted in order to emit a wave train.
“Selecting a time window” of a wave determined by {(t, f(t)), t≥0} defines the selection of an interval [t0, t1] of the time domain of the wave whose wave image or function, ƒ([t0, t1]), comprises one or more complete cycles.
Throughout the present description, “specimen” will be understood as material or sample of material, preferably a tissue, a tissue culture or a cell culture, through which waves emitted by the transducer are made to pass in order to learn about its structural characteristics (elastic parameters, viscoelastic parameters, microstructural geometry, porous, or energy dissipation models, among others). For purposes of the present document, an animal, preferably a human being, through which waves are made to pass, will also be considered a “specimen.”
“Electromechanical actuator” will be understood as a device capable of transforming electrical energy into a movement, particularly a rotational movement. In a particular embodiment suitable for this invention, the electromechanical actuator is stimulated with an electrical signal generated by an electrical pulse generator and is capable of transforming said signal into a minimum fraction of a rotation which will be used to generate the wave that is subsequently analyzed. An example of actuators of this type may consist of an electromagnetic motor.
“Electrical signal” will be understood as an electrical magnitude the value of which depends on time. For purposes of the present invention, constant magnitudes will be considered as particular cases of electrical signals.
Decimal notation: This document uses the symbol “.” as decimal point.
The first object of the present invention is a method (“method of the invention”) for obtaining data relating to the elasticity of materials using torsional waves comprising the following steps of:
By way of illustration,
The non-linearity parameters, such as constitutive non-linearity β or δ, or conventional mechanical non-linearity parameters such as TOECs (third order elastic constants) or others, are parameters that are linked to the elasticity of materials and are used, for example, as a marker for the diagnosis of processes and pathologies in soft tissues.
The method of the invention uses non-linear torsional waves at several frequencies, the propagation speed of which depends directly on shear modulus, as the main indicator of soft tissue consistency. The use of torsional waves offers a higher sensitivity in the detection of irregularities in tissue consistency and has the advantage of eliminating almost all compression waves which contaminate the signal due to their complex propagation modes.
The main advantage of the method of the invention lies in the attainment of non-linearity parameters which allow characterizing the materials using constitutive and/or mechanical non-linearity parameters, the resolution of which is between three and six orders greater than the conventional linear parameters, offering a view at the microstructural level with respect to the onset of pathologies or changes of histological state in the case of tissues, or of microarchitecture in the case of inert materials, which translates into an early and accurate detection of potentially pathological relevant changes in the material.
Each of the features of the method of the invention and the different alternatives giving rise to particular embodiments of said method are described in more detail below:
In a particular embodiment, the emitted wave train consists of between 2 and 80 cycles, preferably between 3 and 10 cycles.
Preferably, the excitation energy used for generating the wave train, in terms of the maximum amplitude of the preceding sinusoidal wave, is comprised between 0.1 V and 20 V, preferably between 2 and 10 V, and the sinusoidal excitation frequency is in the range between 100 Hz and 100 kHz, preferably between 500 Hz and 5 kHz.
It is not easy to obtain this excitation energy when generating torsional waves, so a torsional wave emitter device comprising an electromechanical actuator as will be described below is preferably used.
In a particular embodiment, the specimen is grounded such that the electromagnetic noise is minimized and the signal quality improves.
A time window of the received wave, originating from reflection on the specimen, must be selected such that it occupies a whole number of cycles to prevent the Fourier transform thereof from containing artifacts, where “artifact” is understood as those significant energies at frequencies different from the fundamental excitation frequency and its multiples or harmonics. Preferably, if C is the total number of cycles having reflected wave, a time window consisting of the domain with a length associated with a natural number c, c≤C, of wave cycles is selected, commencing at the moment of the transient fraction of the first cycle, particularly, halfway through the first cycle.
The determination of the moment in which the time window commences is in accordance with the heuristics determined by the type of material to be analyzed. Nevertheless, in another more particular embodiment, the time window commences at the moment of the transient fraction of the first cycle and is associated with a number of cycles comprised between C−2 and C where C is the total number of cycles having reflected wave. In other words, the time window commences at the moment after the commencement of the received wave cycles and is associated with a number of cycles comprised between C−2 and C, where C is the total number of cycles having reflected wave, excluding significantly transient components.
Generally, the greater the selection of the time window is, i.e., the larger the number of cycles associated with selected domain is, the higher the resolution used for analyzing the recorded signal will be.
By using the method of the invention, the behavior and the state of the tissue can be studied in terms of conventional constitutive non-linearity parameters (by analyzing the relationships existing between fundamental harmonics and harmonics of the order higher than or equal to two) such as β or δ, or conventional mechanical non-linearity parameters such as TOECs (third order elastic constants) or others.
In a particular embodiment, once the Fourier transform over the selected time window has been calculated and the amplitudes of harmonics extracted, the constitutive non-linearity or elasticity parameters can be calculated by means of the formula:
where x is the shortest distance between the emitter and receiver, and n is the order of the harmonics that is analyzed.
In another particular embodiment, the elasticity parameters of the specimen are obtained by repeating the method of the invention r times, where r≥2, using identical wave trains with a time interval T>0 between the emission of each wave train and calculating the mean of the calculated non-linearity parameters.
The time interval, T, is preferably greater than or equal to 5 times the duration of the emitted wave train.
Noise is significantly reduced with this average, whereas the time interval prevents the overheating of the emitter.
The emitted torsional waves must preferably have a high signal magnitude, preferably greater than 2 mV, more preferably greater than or equal to 5 mV, so the wave train will be emitted with a sonic and ultrasonic torsional wave emitter device comprising an electrical signal generator connected to an electromechanical actuator which is in turn attached to the contact element, such that when the actuator receives electrical signals, it induces rotational movement of the contact element and upon contacting the specimen, said contact element induces a torsional wave that goes through said specimen.
Any electronic circuit digitalizing the electrical signals at the desired frequencies can be used as the electrical signal generator. Another example of an electrical signal generator may be an oscilloscope, as it allows emitting an electrical signal with a variable voltage over a specific time.
In a preferred embodiment, the electrical signal used for stimulating the electromechanical actuator is an oscillating signal, more preferably a sine-wave signal, and even more preferably a sinusoidal signal, in the cycles claimed in “duty cycle” or work cycle between 1% and 20%, preferably 5%, to prevent overheating the device.
The electromechanical actuator comprising the emitter is preferably covered by a Faraday cage which eliminates electronic noise. Specifically, the electromechanical actuator is wrapped with a conductive covering acting as a Faraday cage.
Another object of the invention is a system comprising the means required for carrying out the method of the invention.
Particularly, the system comprises means for emitting torsional waves, means for receiving torsional waves, and a processor suitable for executing instructions which allow carrying out the method of the invention.
More particularly, the system (
In another particular embodiment, the system further comprises means which allow grounding the specimen. By way of example, these means can consist of a ground connection connected to the support on which the specimen is placed.
In a particular embodiment, the system includes a direct connection between the specimen and the ground connection of the transducer receiving channel. This connection will act as a reference for the input of the received signal.
A fourth object of the invention is a computer program comprising instructions making a computer, connected to the means making up the system of the invention, carry out the method of the invention.
The invention covers computer programs which can be in the form of a source code, an object code, or an intermediate code which is between a source code and an object code, such as in a partially compiled form, or in any other form suitable for use in the implementation of the processes according to the invention. In particular, the computer programs also cover cloud applications implementing the method of the invention.
These programs can be included in or on a medium suitable to be read, hereinafter, a “carrier medium” or a “carrier.” The carrier medium can be any entity or device capable of holding the program. When the program is incorporated in a signal that can be transported directly through a cable or another device or means, the carrier medium may be formed by said cable or another device or means. As a variant, the carrier medium may be an integrated circuit in which the program is included, the integrated circuit being adapted to execute or to be used for executing the corresponding processes.
By way of example, the programs may be incorporated in a storage medium, such as a ROM memory, a CD ROM memory, or a semiconductor ROM memory, a USB memory, or a magnetic recording medium, for example, a floppy disk or a hard disk. Alternatively, the programs may be supported in a transmissible carrier signal. For example, said signal may be an electrical or optical signal that may be transported through an electric or optical cable, by radio, or through any other means.
In this sense, another object of the invention is a computer-readable storage medium comprising program instructions capable of making a computer carry out the method of the invention.
Finally, the last object of the invention relates to a transmissible signal comprising program instructions capable of making a computer carry out the method of the invention.
The embodiment of the method of the invention for obtaining data relating to the elasticity of different materials including tissues, among others, is proposed in a non-exclusive manner.
A system (
Specifically, the emitter is connected, through an amplifier, to a wave generator controlled by the computer by means of an analog/digital converter. The receiver in turn sends the received signal to the analog/digital converter and it is received by the computer which processes the received signal according to the method of the invention.
In turn, the device for emitting and receiving of torsional waves (
For hygienic purposes, the transducer is coated with a latex membrane adapted to the shape of the device. The use of latex assures the dissipation of the wave travelling through it with an adapted involution between the emitter and the receiver.
Using the preceding device, the elasticity of a silicone sample, a connective tissue sample, and a chicken liver tissue sample has been analyzed by executing, by means of instructions interpreted by the computer, the method of the invention with the following characteristics:
where x is the shortest distance between the emitter and receiver, in this case 2.3 mm, and n is the order of the harmonics that is analyzed, in this case 2.
The method of the invention was carried out by emitting torsional waves with different energies on a silicone sample (
The wave train used was a sinusoidal wave which was repeated for 6 cycles and emitted with an excitation energy of 5 and 10 V and a sinusoidal excitation frequency of 800 Hz.
The selection of the time window commenced halfway (50%) through the first cycle, and the domain associated with 5 cycles of the received wave (c=5) is selected.
The method was cyclically repeated for 50 times (r=50) with a time interval, T, of 80 milliseconds during which there was no excitation (0 V) in order to calculate the average of 50 measurements and to thereby reduce noise.
The results obtained after the test are shown in Table 1.
Number | Date | Country | Kind |
---|---|---|---|
P201630123 | Feb 2016 | ES | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/ES2017/070065 | 2/6/2017 | WO | 00 |