Not Applicable.
1. Field of the Disclosure
This disclosure relates generally to ultrasonic transducers, and more particularly, to a system and method for driving ultrasonic transducers and improvements to a phase track controller for reducing loss of lock occurrence.
2. Background of the Disclosure
Ultrasonic transducers have been in use for many years. Current technology typically depends on resonant circuits to drive ultrasonic transducers. Resonant circuits are, by definition, designed to operate in a very narrow range of frequencies. Because of that, the transducer tolerances are held very tightly to be able to operate with the driving circuitry which results in higher costing circuits. In addition, there is no possibility of using the same driving circuit for transducers with different frequencies because of the static nature of typical drivers, and the circuit must be changed for every transducer frequency. There is also a need for a system and method for driving any transducer regardless of the resonant frequency of the transducer.
U.S. Pat. No. 8,115,366 B2 describes a linear phase track controller which can accurately generate a wide range of frequencies and is capable of individually driving multiple transducers with different resonant frequencies.
A disadvantage of a linear phase track controller occurs when it is no longer able to track the output phase, due to noise or some other perturbation in the system. This condition is called loss of lock. This disclosure describes an enhancement to a linear phase track controller which greatly reduces loss of lock occurrences.
The features and advantages of the disclosure will be set forth in the description which follows, and in part will be apparent from the description, or may be learned by the practice of the disclosure without undue experimentation. The features and advantages of the disclosure may be realized and obtained by means of the instruments and combinations particularly pointed out herein.
Briefly and in general terms, the disclosure is directed to improvements to a linear phase track controller used to drive an ultrasonic transducer for reducing loss of lock occurrence.
The disclosure will describe an implementation of methods and systems to drive ultrasonic transducers. It will be appreciated that a method is often required to generate a wide range of frequencies with high accuracy and very high frequency shifting speed. Such a system and method may individually drive multiple transducers each having a different frequency, thereby allowing device manufacturers to take advantage of economies of scale by implementing the same driver with various transducers having different frequencies.
In aspects of the disclosure, a system may comprise a controller adapted to provide a voltage and a frequency, the controller may be configured to vary the voltage based on a current error signal derived from a drive current through a transducer and from a current command, the controller may be configured to vary the frequency based on at least one parameter indicative of whether the transducer is at or near a resonance state. The system may also comprise a drive adapted to receive the voltage and the frequency from the controller, and may be adapted to provide a drive voltage at a drive frequency to the transducer based on the voltage and the frequency received from the controller, the drive voltage being at a level that maintains the drive current at substantially the current command, the drive frequency being at substantially a resonant frequency of the transducer. In further aspects of the disclosure, the at least one parameter includes a phase angle between the drive current and the drive voltage.
In aspects of the disclosure, a method may comprise providing a drive voltage at a drive frequency to a transducer, the drive voltage causing a drive current through the transducer. The method may further comprise sensing the drive current and determining a current error from the sensed drive current and from a current command. The method may further comprise adjusting the drive voltage based on the current error, and determining at least one parameter from the sensed drive current and from the voltage level, the at least one parameter indicative of whether the transducer is at or near a resonance state, the at least one parameter including a phase angle between the drive current and the drive voltage. The method may further comprise adjusting the drive frequency based on the at least one parameter, including maintaining the drive frequency at or substantially at a resonant frequency of the transducer.
The “measured phase” of a system may be defined as the phase angle between the transducer voltage and transducer current. The “commanded phase” of the system is the phase angle at which a user wants the transducer to be driven. A linear phase track controller operates by adjusting the drive frequency to drive the difference between the measured phase and a commanded phase (this difference is called the “phase error”) to zero. In implementations of the disclosure, a commanded phase, selected to be well below the anticipated minimum peak phase and to be near the resonant frequency of the transducer, is specified. The measured phase of the system, defined as the phase angle between the transducer voltage and transducer current, is fed back and instantaneously subtracted from the commanded phase to form the phase error. In response to a positive phase error (the actual phase is less than the commanded phase), the phase track controller generates a positive frequency step (ΔFrequency), increasing the drive frequency and thereby decreasing the phase error. Similarly, in response to a negative phase error (the actual phase is greater than the commanded phase), the phase track controller generates a negative ΔFrequency, decreasing the drive frequency and thereby again decreasing the phase error.
When the phase error is at or near zero, the linear phase track controller is correctly driving the transducer. If the phase error becomes too great, phase tracking may be lost. Loss of phase tracking (also called loss of lock) can occur when the controller is used in systems with significant phase noise. In addition, as the physical load on the transducer changes, the phase response (and the admittance response also) can shift in frequency, become wider or narrower, and/or can have a greater or lesser peak value. If this physical load change occurs rapidly with respect to the system sampling interval, measured phase can transition to the high frequency (or anti-resonance) side of the transducer phase curve, also ultimately resulting in loss of phase lock.
In aspects of the disclosure, the phase response curve with phase operating regions may be defined. As long as the sensed phase remains within the phase band (i.e., near the operating point), the linear phase track controller executes. However, if phase is sensed to be outside the phase band, conditionally, linear phase track control may be suspended and, one of two nonlinear mechanisms may be executed to drive the phase back into the phase band.
In aspects of the disclosure, the frequency source select module determines the source of the frequency step, ΔFrequency, applied to the frequency generator, switching among the original phase track controller and the “up search” and “down search” generators. When the sensed phase is within the specified phase band, the phase track controller generates ΔFrequency. If the sensed phase is greater than the upper edge of the phase band and the frequency is increasing, the down search generator is activated. The down search generator generates a pre-defined number of negative ΔFrequency steps, while also looking for a maximum admittance. After completion of the pre-defined number of steps, the frequency is set to the frequency of the sensed maximum admittance and the phase track controller re-enabled to as the source of the frequency stepping.
Similarly, if the sensed phase is less than the upper edge of the phase band and the frequency is decreasing, the up search generator is activated. The up search generator generates a pre-defined number of positive ΔFrequency steps, while also looking for a maximum admittance. After completion of the pre-defined number of steps, the frequency is set to the frequency of the sensed maximum admittance and the phase track controller re-enabled to as the source of the frequency stepping.
In aspects of the disclosure, various implementations of the up search and down search generators and the frequency source select components may be utilized. This logic is executed periodically, nominally at 1 to 100 millisecond intervals.
The features and advantages of the disclosure will be more readily understood from the following detailed description, which should be read in conjunction with the accompanying drawings.
Non-limiting and non-exhaustive implementations of the disclosure are described with reference to the following figures, wherein like reference numerals refer to like parts throughout the various views unless otherwise specified. Advantages of the disclosure will become better understood with regard to the following description and accompanying drawings where:
a is a flow diagram showing waveforms into and out of an exclusive OR gate of the circuit of
In the following description of the disclosure, reference is made to the accompanying drawings, which form a part hereof, and in which is shown by way of illustration specific implementations in which the disclosure may be practiced. It is understood that other implementations may be utilized and structural changes may be made without departing from the scope of the disclosure.
Some embodiments of the disclosure may involve hardware and software. The hardware may include a switching amplifier to create a sine wave output to an ultrasonic transducer. The ultrasonic transducer can or may be a piezoelectric transducer. The switching amplifier can be run with high efficiency over a broad range of frequencies and can, therefore, be used to drive transducers of many frequencies. The switching amplifier can also drive transducers that do not have tightly held frequency tolerances thereby reducing transducer production cost. This allows for reduction of production cost due to economies of scale and allows for customers that use different frequency transducers to always be able to use the same driver.
Previous ultrasonic generators have relied on resonant power sources or analog amplifiers to drive the transducer. In some embodiments of the disclosure, a class D or class E amplifier may be used to amplify the output of a digitally controlled AC source. This technique frees the manufacturer and user from the requirement of designing a resonant system around a specific transducer. Instead, this system is usable for any transducer over a broad range of frequencies.
Previous class D and class E amplifiers have used traditional LC or cascaded LC filters to significantly reduce the effects of the class D or E carrier frequency on the signal frequency. In some embodiments of the disclosure, a two phase output signal may be used in conjunction with a coupled transformer to reduce the effect of the carrier frequency to several times lower than could be done with similar size and cost components with the traditional LC type filters.
In some embodiments of the disclosure, software could run entirely on low cost, 16-bit, integer-only microcontrollers. The more powerful DSP (digital signal processor) modules typically required in the field are not required in the disclosure, although DSP modules could be used in some embodiments.
A method may be required to generate a wide range of frequencies with high accuracy and very high frequency shifting speed. A digital synthesizer could be used in an ultrasonic system to allow rapid and flexible frequency control for output of a frequency generator.
In some embodiments, dead time may be minimized in switching circuits in order to minimize the output impedance to the transducer. As used herein, the phrase “dead time” is the time in power switching circuits when all switching elements are off to prevent cross conduction. A minimum or maximum admittance may be used when determining the resonant frequency. The admittance measured will vary much less between in resonance and out of resonance in a low Q system than in a high Q system. The dimensionless parameter “Q” refers to what is commonly referred to in engineering as the “Q factor” or “quality factor.” Because Q is directly affected by the impedance of the driving circuit, this impedance must be kept very low. In addition to the commonly considered impedances of the output transformer, driving semiconductors, PCB (printed circuit board) and other directly measureable impedances, Applicants have found that the dead time has a very strong effect on the output impedance of the driver. As such, the switching circuit may be configured to have a very small (approximately 50 nanoseconds) dead time. In some embodiments, the switching circuit may have a dead time that may be greater than or less than 50 nanoseconds.
For optimum operation, the transducer may be run at or near its resonant frequency point. The resonant frequency point of the transducer is defined as the frequency at which maximum real power is transferred from the drive amplifier to the transducer.
Applicants have found that the admittance of the transducer provides a reliable indication of the proximity of the transducer to its resonant frequency point. Admittance is defined as the RMS (root-mean-square) amplitude of the transducer drive current divided by the RSM amplitude of the transducer drive voltage.
The circuit 10 shown in
The circuit in
Applicants have also found that the phase of the transducer also provides a reliable indication of the proximity of the transducer to its resonant frequency point. Phase is defined as the phase angle between the transducer drive voltage and transducer drive current.
The circuit shown in
As shown in
The drive 208 may provide a drive signal of controlled voltage and controlled frequency to the transducer 210. An output parameter sense circuit 212 senses transducer drive voltage and transducer drive current and generates a measure of current 218, admittance 220, and a frequency control parameter 222. The frequency control parameter may be different in different embodiments.
Current 218 is applied as an input to the current controller 202, which may generate a voltage 214 applied to the drive 208. The current controller 202 may set the voltage 214 to maintain the current required for correct operation of the transducer 210 in its given application.
The frequency controller 206 may perform two functions: frequency scanning and frequency tracking. The frequency scanning function searches for a frequency that is at or near the resonant frequency of the transducer. The frequency tracking function maintains the operating frequency at or near the resonant frequency of the transducer.
When the frequency controller 206 is frequency scanning, admittance 220 is applied to it as an input. The frequency controller sweeps the drive frequency over a range of frequencies appropriate for the transducer and application, searching for the resonant frequency.
When the frequency controller 206 is frequency tracking, a frequency control parameter 222 is applied to it as an input. The frequency controller sets the frequency required for correct operation of the transducer in its given applications.
When the frequency controller 206 performs either frequency scanning or frequency tracking, it applies the calculated frequency 216 to the drive 208.
The drive 208 may include the switching amplifier and switching circuits described above. The frequency controller 206 may include the digital synthesizer described above.
As previously mentioned, the Linear Phase Track Controller frequency controller 206 may perform two functions: frequency scanning and performs phase frequency tracking.
In an implementation, initial application of drive to the transducer at its resonant frequency may be highly desirable. Due to variations in transducer characteristics, applied power levels, and the mechanical load the transducer connects to, the resonant frequency may not be known and the frequency controller may perform a frequency scan to establish the drive frequency at or near the resonant frequency for the transducer.
When performing a frequency scan, the frequency controller searches a predefined range of frequencies for the frequency at which the transducer admittance is maximum. As shown in
In some embodiments, admittance may be detected after each narrow scan settling time and, at completion of the narrow scan, the drive frequency is set to the frequency of maximum detected admittance.
In some embodiments, phase is detected after each narrow scan settling time and, at completion of the narrow scan, the drive frequency is set to the frequency with detected phase closest to the phase required for correct operation of the transducer in its given application.
An ultrasonic transducer may often have multiple frequencies at which the commanded phase is measured. The frequency of maximum admittance will always be at or close to the resonant frequency, the frequency of maximum real power transfer. For this reason, maximum admittance may be used for wide and medium scans for the operating point, regardless of the method used in the narrow scan.
In an implementation, the frequency scanner 300 can be executed at either full power (as defined by the user) or at a predefined low power of less than 5 watts, measured at transducer resonance.
The frequency controller 206 may optionally perform a fast scan 308 as part of its operation, immediately prior to initiation of a frequency track algorithm. The fast scan includes a .+−.10 Hz sweep about the current frequency, in 2 Hz steps, with a 10 msec settling time after each step.
In some embodiments, admittance is detected after each fast scan settling time and, at completion of the fast scan, the drive frequency is set to the frequency of maximum detected admittance.
In some embodiments, phase may be detected after each fast scan settling time and, at completion of the fast scan, the drive frequency is set to the frequency with detected phase closest to the phase required for correct operation of the transducer in its given application. The fast scan 308 can be executed at either full power or at less than 5 watts power.
The transducer resonant frequency may fluctuate during normal operation. This fluctuation may occur due to changes in operating conditions of the transducer, such as changes in temperature of the transducer and mechanical load on the transducer. Phase tracking may be performed to compensate for this fluctuation in resonant frequency.
If the detected admittance has increased by greater than a predefined amount, the next step 418 is taken in the same direction as the previous step, with step size based on the magnitude of the increase in admittance. For example, the magnitude of the step can be proportional to the detected increase in admittance. If the detected admittance has decreased by greater than a predefined amount, the next step 418 is taken in the opposite direction, with the magnitude of the step being based on the magnitude of the increase in admittance. If the detected admittance has neither increased by greater than a predefined amount nor decreased by greater than a predefined amount, the admittance is assumed to be at its peak and a zero magnitude “step” is taken. The frequency tracker delays a short time period to allow the transducer to settle and the peak detection and step sequence is repeated.
The maximum admittance of a transducer may increase, remain unchanged, or decrease, depending on changes in operating conditions of the transducer. Frequency tracking for increasing and unchanging maximum admittance values is performed by the above-described frequency tracking method. Tracking the resonant frequency associated with a decreasing admittance maximum is performed by stepping quickly in equal magnitude steps in both directions about the current frequency until the decrease in admittance stops and increased admittance values are again detected. The Frequency Controller then changes the frequency to again lock on the point of maximum admittance.
The frequency tracking method described above can be implemented with an algorithm within software being run by the hardware of the system 200.
Another embodiment of the frequency tracker, shown in
The frequency tracker 500 performs frequency tracking by applying a phase angle error term 520 to a Proportional-Derivative (PD) controller 502 at regular sampling intervals of between 5 and 20 msecs. The phase angle error term is calculated to be the difference between the phase track command 518 and the measured transducer phase 516. The PD controller 502 includes a differentiator, A 502a, a proportional gain, KFP 502b, a differential gain, KFD 502c, and an output gain, KFO 502d. The output from the PD controller 502 in response to a phase error 520 is a step in frequency, A frequency 512, of magnitude and sign necessary to drive the phase error 520 toward zero. The step in frequency 512 is applied to the frequency generator 504 which calculates the new frequency 514. The driver drives the transducer 508 at the frequency 514 from the frequency generator 504.
The current controller 600 varies the current through the transducer by varying the drive voltage applied across the transducer. Increasing the drive voltage increases the transducer current and decreasing the drive voltage decreases the transducer current. In some embodiments, the current controller 600 provides a voltage 610 to the drive 604, and this voltage is provided by the drive 604 to the transducer 606.
At a regular sampling intervals, ranging between 5 and 20 msecs, the current controller 600 samples the transducer current and converts it to an RMS current value 612 by an RMS converter 608. At each sampling interval the current controller 600 calculates a current error term 616 by subtracting the sample of the output RMS current 612 from the commanded current 614.
The current controller 600 applies a current error term 616 to a Proportional-Integral-Derivative (PID) controller 602, which generates a response 610 to the error 616. The error 616 is integrated by an integrator 602a and differentiated by a differentiator 602b. The error 616 and its integral and differential are multiplied respectively by the P, I, and D gains, 602c, 602d, 602e internal to the PID controller, summed, and their sum multiplied by the controller output impedance factor KCO 602f to form the controller output voltage 610. Controller gains, 602c, 602d, 602e, 602f are set to achieve maximum rise time with an approximately 10% overshoot in the output response to a step in the input. The output impedance factor 602f provides both scaling and translation from current to voltage. The controller output voltage 610 is applied to driver 604 to be amplified to become the transducer drive voltage.
In some embodiments, the current controller 600 employs two output impedance factors 602f. A larger output impedance factor may be used for the first period of time (nominally 500 msecs) to assure the transducer reaches its steady-state behavior at the given drive power, physical load, and temperature as rapidly as possible. A smaller output impedance factor may be used once the transducer has reached its steady-state behavior. When the switch from the first to the second output impedance factor occurs, the integral of the current error maintained by the PID controller is modified to prohibit an undesired transient in the transducer drive voltage.
In
To achieve balanced operation, the controller scheduler 204 interleaves the operation of the frequency controller 206 and the current controller 202.
When the frequency controller is performing a scan or search operation, the controller scheduler disables the current controller.
When the frequency controller is tracking frequency, in some embodiments the controller scheduler alternates the operation of the two controllers. That is, a controller will execute every 5N msecs, with the current controller executing for odd N and the frequency controller executing for even N.
In some embodiments, both controllers are allowed to operate simultaneously, except immediately after a frequency step. When the frequency controller is tracking frequency, the controller scheduler disables the current controller for the first M 5-msec periods after a frequency step. The number of periods, M, is typically 2, but can be more or less than 2. At the end of the M periods, the frequency control parameter is now only a result of the step in frequency and not of control exerted by the current controller. The frequency control parameter is sampled at this time and stored for the next frequency controller calculation, and the controller scheduler re-enables the current controller.
The output of the processor running the code discussed previously is a small signal with all the characteristics of necessary to drive and ultrasonic transducer except for the amplitude. The drive circuit 208, 408, 506 can be broken down into two sections as shown in
Prior art has used linear amplifiers for this drive section. These have the disadvantages of being large, inefficient and costly. The illustrated embodiment of
In some embodiments, the drive 208, 408, 506 includes filter circuitry. In some embodiments with a transducer operational range of 20 kHz to 60 kHz, the filter circuitry may be configured to have a corner frequency higher than 60 kHz to avoid excessive resonant peaking depending on the type of transducer and its intended use, it will be appreciated that the transducer operational range can be lower than 20 kHz and/or higher than 60 kHz, and the filter circuitry can be configured to have a corner frequency higher than the transducer operational range. The carrier frequency used may be about 10 times that of the transducer resonance frequency in some implementations.
In some embodiments the filter circuitry is configured to reduce transmission of the carrier frequency (Fs) from a switching amplifier of the drive 208, 408, 506. Non-limiting examples of filter circuitry are described below.
In previous art, the output filter of a switching amplifier is typically implemented with an LC or cascaded LC filter. An example of a cascaded LC filter is shown in
Part of this disclosure is a new form of output filter that includes a coupled inductor as part of the output filter. An example schematic of this new coupled LCLC filter is shown in
To take advantage of the coupled inductor, a second change is made to the system. The class D or E amplifier from
The described phase shift between two or more channels can be found in prior art, for example in multiphase buck converter applications, or in U.S. Pat. No. 6,362,986 to Shultz et al., entitled “Voltage converter with coupled inductive windings, and associated methods.” U.S. Pat. No. 6,362,986 represents closer prior art, as it has phase shift together with magnetic coupling between inductors, as illustrated in
Notice that the output voltage of circuit in
Magnetic coupling between windings in
The coupled inductor from
Waveforms for the circuit in
The decreased current ripple offers several benefits to the circuit and its performance. Decreased current ripple makes it easier for the output filter to achieve low noise levels and low output voltage ripple at the output, in other words—either smaller attenuation could be used as compared to the case without magnetic coupling, or lower noise level can be achieved. Decreased amplitude of the current ripple also means that the RMS value of the current waveform is lower, which relates to lower conduction losses. Lower current ripple also implies lower peaks of the current, which relates to the lower stress in switching devices of the power circuits. As the DC component of the load current is the same in both coupled inductors (the outputs are connected to each other through the load so the load current is equal), and since these currents create opposite magnetic flux for arrangement shown in FIG. 14—cancellation of the DC component of the magnetic flux in the core is beneficial for the small core size and low core losses. The decrease of the current ripple is generally good for EMI decrease, and makes it easier to pass regulatory requirements. While the performance of the filter in terms of the amplifier signals is dependent on the leakage inductance values, the noise signals of the Common Mode (same in both output nets) will be attenuated by much larger magnetizing inductance. In this regard, Common Mode noise, often being present in circuits and representing a need for additional high frequency filtering for the output connections, will be attenuated at much higher degree in magnetically coupled inductor arrangement in
The phase shifted PWM2 signal for the second differential amplifier circuit in
The magnetic components from
The above described transducer can be a part of or contained in any type of apparatus, including without limitation a surgical device, a cutting tool, a fragmentation tool, an ablation tool, and an ultrasound imaging device.
Referring now to
The method described assumes the operating point is defined to be at or near the resonant frequency, where admittance is maximum. However, the same approach may be used, with appropriate sign reversals, for operation at the anti-resonant frequency, where admittance is minimum.
Similarly, in response to a negative phase error (the measured phase is greater than the commanded phase), the linear phase track controller generates a negative ΔFrequency, decreasing the drive frequency and thereby again decreasing the phase error.
If the phase error becomes too great, phase tracking can be lost. Loss of phase tracking (also called loss of lock) can occur when the controller is used in systems with significant phase noise. In addition, as the physical load on the transducer changes, the phase response (and the admittance response also) can shift in frequency, become wider or narrower, and/or can have a greater or lesser peak value. If this physical load change occurs rapidly with respect to the system sampling interval, measured phase can transition to the high frequency (or anti-resonance) side of the transducer phase curve, also ultimately resulting in loss of lock.
The following describes a method that anticipates loss of lock and corrects for it prior to its occurrence.
When the sensed and measured phase is within the specified phase band, the frequency source select selects the linear phase track controller for generating frequency steps (ΔFrequency) either increasing or decreasing.
If the sensed phase is greater than the upper edge of the phase band and the frequency is increasing, the frequency source select selects the down search generator. The down search generator generates a pre-defined number of negative ΔFrequency steps, while also looking for a maximum admittance. After completion of the pre-defined number of steps, the frequency is set to the frequency of the sensed maximum admittance and frequency source select re-selects the linear phase track controller as the source of the frequency stepping.
If the measured phase is greater than the upper edge of the phase band and the frequency is decreasing, the frequency source select continues to select the linear phase track controller as the source of the frequency stepping.
In an implementation, if the sensed measured phase is less than the upper lower edge of the phase band and the frequency is decreasing, the frequency source select selects the up search generator, which is activated. The up search generator generates a pre-defined number of positive ΔFrequency steps, while also looking for a maximum admittance. After completion of the pre-defined number of steps, the frequency is set to the frequency of the sensed maximum admittance and the frequency source select re-selects the linear phase track controller re-enabled to as the source of the frequency stepping.
In an implementation, if the measured phase is less than the lower edge of the phase band and the frequency is increasing, the frequency source select continues to select the linear phase track controller as the source of the frequency stepping.
In an implementation, if the measured phase falls below the loss of lock threshold, loss of lock has occurred and methods outside the scope of this description must be implemented to regain phase lock.
In an implementation, frequency steps may change linearly, exponentially, and logarithmically.
In an implementation, it may be determined that phase is corrected when it deviates more than 30 degrees out of phase and that ten (10) steps of up and down searching will be performed, wherein each step is 5 Hz. However, any number of steps may be used as desired.
The foregoing description has been presented for the purposes of illustration and description. It is not intended to be exhaustive or to limit the disclosure to the precise form disclosed. Many modifications and variations are possible in light of the above teaching. Further, it should be noted that any or all of the aforementioned alternate implementations may be used in any combination desired to form additional hybrid implementations of the disclosure.
Further, although specific implementations of the disclosure have been described and illustrated, the disclosure is not to be limited to the specific forms or arrangements of parts so described and illustrated. The scope of the disclosure is to be defined by the claims, if any, any future claims submitted here and in different applications, and their equivalents.
This application claims the benefit of U.S. Provisional Patent Application No. 61/843,306, filed on Jul. 5, 2013, which is hereby incorporated by reference herein in its entirety, including but not limited to those portions that specifically appear hereinafter, the incorporation by reference being made with the following exception: In the event that any portion of the above-referenced application is inconsistent with this application, this application supersedes said above-referenced application.
Number | Date | Country | |
---|---|---|---|
61843306 | Jul 2013 | US |