The present application claims priority under 35 U.S.C. §119 to Japanese Patent Application No. 2012-243949, filed Nov. 5, 2012. The contents of this application are incorporated herein by reference in their entirety.
1. Field of the Invention
The present invention relates to a controller of an AC motor.
2. Discussion of the Background
In controllers of AC (alternating current) motors, their driving control in a constant output region is generally the control of setting a current command on the d axis, which is parallel to the flux of the AC motor, into the negative direction, thereby weakening the flux. This control is also referred to as voltage limiting control (see, for example, Japanese Unexamined Patent Application Publication No. 2010-022165).
According to one aspect of the present invention, a controller of an AC motor includes a d-axis voltage command section, a d-axis non-interactive control section, a first current deviation arithmetic section, a q-axis integral control section, a q-axis voltage command section, a constant output control section, and a d-axis voltage command correction section. The d-axis voltage command section is configured to generate a d-axis voltage command on a d axis of a d-q coordinate system. The d axis is parallel to a flux of the AC motor and orthogonal to a q axis of the d-q coordinate system. The d-axis non-interactive control section is configured to remove, from the d-axis voltage command, an interference component resulting from a current on the q axis. The first current deviation arithmetic section is configured to perform an arithmetic operation to obtain a q-axis current deviation between a current command on the q axis and the current on the q axis flowing through the AC motor. The q-axis integral control section is configured to receive the q-axis current deviation and output an integral value of the q-axis current deviation. The q-axis voltage command section is configured to generate a q-axis voltage command on the q axis based on the q-axis current deviation and is configured to output the q-axis voltage command. The constant output control section is configured to output a correction voltage command relative to the d-axis voltage command based on an output of the q-axis integral control section. The d-axis voltage command correction section is configured to subtract the correction voltage command from the d-axis voltage command after the d-axis non-interactive control section has performed non-interactive control, so as to correct the d-axis voltage command.
A more complete appreciation of the present disclosure and many of the attendant advantages thereof will be readily obtained as the same becomes better understood by reference to the following detailed description when considered in connection with the accompanying drawings, wherein:
Embodiments of a controller of an AC motor (hereinafter referred to as “motor controller”) disclosed in the present application will be described in detail below by referring to the accompanying drawings. The following embodiments are provided for exemplary purposes only and are not intended in a limiting sense.
First, a motor controller according to the first embodiment will be described.
As shown in
The power conversion section 10 includes a three-phase inverter circuit 13 and a PWM signal generation section 14. The three-phase inverter circuit 13 is coupled between the DC power source 2 and the motor 3. The three-phase inverter circuit 13 is made up of, for example, six switching elements in three-phase bridge connection. Based on a control signal from the vector control section 12, the PWM signal generation section 14 generates a PWM signal to turn on and off the switching elements constituting the three-phase inverter circuit 13, and outputs the PWM signal to the three-phase inverter circuit 13. The configuration of the DC power source 2 may also be to convert AC power into DC power and output the DC power, examples including a combination of a rectifier circuit of diode and a smoothing capacitor that smoothes out DC output voltage. In this case, an AC power source is coupled to the input side of the rectifier circuit.
The current detection section 11 detects current flowing between the power conversion section 10 and the motor 3. Specifically, the current detection section 11 detects instantaneous values iu, iv, and iw of the currents flowing between the power conversion section 10 and a U phase, a V phase, and a W phase of the motor 3 (the instantaneous values being hereinafter referred to as output currents iu, iv, and iw). An example of the current detection section 11 is a current sensor that detects current using a Hall device, which is a magneto-electric converting device.
The vector control section 12 generates a control signal and outputs the control signal to the power conversion section 10. The control signal is based on the output currents iu, iv, and iw detected by the current detection section 11 and based on a rotor electrical angle phase θ of the motor 3 detected by a position detection section 4 (the electrical angle being defined as the mechanical angle of the rotor of the motor 3 multiplied by the number of magnetic pole pairs of the motor 3, which applies throughout the description that follows). In a d-q coordinate system in which an axis parallel to the flux of the motor 3 is a d axis and an axis having a direction orthogonal to the d axis is a q axis, the vector control section 12 divides the current components into a d-axis component and a q-axis component in performing vector control.
In the following description, the d-axis component and the q-axis component of the current command will be respectively referred to as a d-axis current command id* and a q-axis current command iq*, and the d-axis component and the q-axis component of the current flowing through the motor 3 will be respectively referred to as a d-axis current id
When the motor 3 is an IPM (Interior Permanent Magnet) motor, a voltage equation on the d-q coordinate system can be represented by Formula (1). While the following description will assume that the motor 3 is an IPM motor, the motor 3 will not be limited to the IPM motor. For example, when the motor 3 is an SPM (Surface Permanent Magnet) motor, then Ld=Lq.
In Formula (1), id and iq respectively denote the d-axis component and the q-axis component of the current flowing through the motor 3, and Vd and Vq respectively denote the d-axis component and the q-axis component of the voltage applied on the motor 3. Also, R denotes the coil resistance of the motor 3, Ld denotes the d-axis inductance of the motor 3, Lq denotes the q-axis inductance of the motor 3, ω denotes the electrical angular velocity of the motor 3, φ denotes an induced voltage constant, and p denotes a differential arithmetic section. It is noted that R, Ld, Lq, and φ are motor parameters.
The subtraction section 24 performs an arithmetic operation to obtain a d-axis current deviation, which is a deviation between the d-axis current command id* and the d-axis current id
The non-interactive control section 29 is provided to avoid interference between the d axis and the q axis. The non-interactive control section 29 generates a d-axis voltage compensation value vd
The constant output control section 34 subtracts, from an integral value ΣACRq of the q-axis current deviation, a value obtained by multiplying the q-axis current iq
The subtraction section 30 subtracts the d-axis voltage compensation value vd
In the high velocity region of the motor 3, a relationship represented by Formula (2) is established. Hence, a dominating factor of the voltage saturation in the constant output region is the q-axis voltage command vq**, which is a torque-axis voltage command. When the q-axis voltage command vq** becomes saturated, the current control of the q-axis current stops functioning, resulting in degraded torque responsivity.
|vd|<|vq| (2)
When the above-described motor parameters have no errors and the non-interactive control section 29 accurately performs its control with the q-axis voltage command vq** in non-saturation state, then a q-axis integral control section 52, described later, of the q-axis current control section 27b only outputs an amount equivalent to the voltage drop that is due to the coil resistance R. When the q-axis voltage command vq** becomes saturated, the q-axis integral control section 52 increases its output. Hence, the difference between the output of the q-axis integral control section 52 and the amount equivalent to the voltage drop that is due to the coil resistance R indicates the degree of saturation of the q-axis voltage command vq**.
Incidentally, when the d-axis current id
Specifically, the motor controller 1 includes the constant output control section 34 to perform the control of generating the correction voltage command Δvd* relative to the d-axis voltage command vd**′ based on the integral value ΣACRq of the q-axis current deviation, and subtracting the correction voltage command Δvd* from the d-axis voltage command vd**′. This ensures control of the q-axis current iq
In the block diagram shown in
It is also assumed that with a limiter of a d-axis integral control section 42, described later, of the d-axis current control section 27a being set to have a low upper limit, the output of the d-axis integral control section 42 is saturated as indicated by a value on the limiter in the constant output state. In this case, the d-axis integral control section 42 can be omitted in the d-axis current control section 27a. Also in the constant output state, the q-axis voltage command vq** is saturated, and the q-axis current control section 27b (excluding the q-axis integral control section 52) and the non-interactive control section 29 are stopping functioning on their q axis control. Hence, it is possible to omit this portion (which is the portion associated with generation of the q-axis voltage command vq* and the q-axis voltage command vq**). Additionally, when a voltage error compensation section 33, described later (see
Thus, the block diagram shown in
In Formula (3), Kp
When id*≈id
From Formulae (4) and (5), Formula (7) is derived, and further, from Formula (7), Formula (8) is derived.
In the constant output state, the q-axis voltage command vq** is saturated and thus is at a fixed value (hereinafter referred to as saturation limit value). Then, by making the saturation limit value of the q-axis voltage command vq** into ωφ, Formula (6) can be simplified into Formula (9).
From Formula (9), Formula (10) is derived.
Formula (10) is substituted into Formula (8), and thus Formula (11) is derived. Further, from Formula (11), Formula (12) is derived.
Thus, a transfer function GT(s) from the q-axis current command iq* to the q-axis current iq
Here, when a current control gain is set as represented by Formula (14) with a current control response ωACR [rad/s] as a parameter, the transfer function GT(s) represented by Formula (13) can be represented by Formula (15).
The transfer function GT(s) represented by Formula (15) contains the electrical angular velocity ω of the motor 3, and thus the torque response is velocity dependent. In view of this, in the motor controller 1 according to this embodiment, a constant output control gain G(s) in the constant output control section 34 is set to be inversely proportional to the electrical angular velocity ω. This eliminates the velocity dependency of the torque response. The electrical angular velocity ω is in proportional relationship with the electrical-angle rotational frequency of the motor 3 and with the output frequency of the voltage command (because when the motor 3 is a synchronous motor, the electrical-angle rotational frequency of the motor 3 matches the output frequency of the motor controller). Hence, setting the constant output control gain G(s) to be inversely proportional to the electrical-angle rotational frequency of the motor 3 and to the output frequency of the voltage command eliminates the velocity dependency of the torque response. As used herein, the output frequency of the voltage command refers to the frequency of the output voltage specified by the output voltage command.
When the control by the constant output control section 34 is P control (proportional control), as the voltage saturation develops, a steady-state deviation occurs in the integral value ΣACRq of the q-axis current control, resulting in degraded accuracy of the q-axis current control. In view of this, as the method of control by the constant output control section 34, such a control method is employed that I control (integral control) is added to P control.
When the control by the constant output control section 34 is PI control, the characteristic equation becomes one order higher into a fourth-order characteristic equation, which makes the designing complicated in PI control. In view of this, to facilitate the designing, the control by the constant output control section 34 may be PID control.
When the control by the constant output control section 34 is PID control, the constant output control gain G(s) of the constant output control section 34 can be represented by Formula (16). In Formula (16), Kp
From Formula (16), Formula (15) can be represented by Formula (17).
Thus, the motor controller 1 according to this embodiment includes the constant output control section 34, which outputs a correction voltage command relative to the d-axis voltage command vd**′ based on the output of the q-axis integral control section 52 of the q-axis current control section 27b. Then, the motor controller 1 subtracts the correction voltage command Δvd* from the d-axis voltage command vd**′, thereby obtaining the d-axis voltage command vd**′.
Specifically, the motor controller 1 controls the d-axis voltage command vd** based on the output of the q-axis integral control section 52, thereby performing constant output control of the q-axis current iq
An exemplary detailed configuration of the motor controller 1 according to this embodiment will be described in detail below by referring to
As shown in
The three-phase/two-phase conversion section 21 converts each of the output currents iu, iv, and iw into αβ components of two orthogonal axes on a fixed coordinate system, and obtains a fixed coordinate current vector on a αβ-axes coordinate system, which has, as vector components, an output current iα in the α axis direction and an output current iβ in the β axis direction.
Based on the rotor electrical angle phase θ, which is detected by the position detection section 4 and indicates the rotor position of the motor 3, the d-q coordinate conversion section 22 converts the components on the αβ-axes coordinate system output from the three-phase/two-phase conversion section 21 into components on a d-q axis coordinate system. In this manner, the d-q coordinate conversion section 22 obtains the d-axis current id
The subtraction section 24 subtracts the d-axis current id
The subtraction section 26 subtracts the q-axis current iq
The current control section 27 performs PI control of the d-axis current deviation, thereby generating the d-axis voltage command vd*, and outputs the d-axis voltage command vd* to the subtraction section 30. Also the current control section 27 performs PI control of the q-axis current deviation, thereby generating the q-axis voltage command vq*, and outputs the q-axis voltage command vq* to the addition section 31. Further, the current control section 27 outputs the integral value ΣACRq of the q-axis current deviation to the voltage error compensation section 33 and the constant output control section 34.
The d-axis proportional control section 41 performs proportional control at a proportional gain Kp
When the integral value of the integral section 45 exceeds the upper limit or falls below the lower limit, the limiter 46 limits the output of the integral value of the integral section 45 to the upper limit or the lower limit.
The d-axis current id
The d-axis voltage command section 43 adds the output of the d-axis integral control section 42 to the output of the d-axis proportional control section 41, thereby generating the d-axis voltage command vd*.
The q-axis current control section 27b includes a q-axis proportional control section 51, the q-axis integral control section 52, and a q-axis voltage command section 53. The q-axis proportional control section 51 performs proportional control at a proportional gain Kp
The q-axis voltage command section 53 adds the output of the q-axis integral control section 52 to the output of the q-axis proportional control section 51, thereby generating the q-axis voltage command vq*. Also the output of the q-axis integral control section 52, which is the integral value ΣACRq, is output to the voltage error compensation section 33 and the constant output control section 34.
Referring back to
The non-interactive control section 29 generates the d-axis voltage compensation value vd
The LPF 61 removes a high-frequency component of the d-axis current id
vq
The LPF 65 removes a high-frequency component of the q-axis current iq
vd
As described above, the d-axis current id
In view of this, the non-interactive control section 29 generates the q-axis voltage compensation value vq
Referring back to
The voltage error compensation section 33 identifies the components contained in the integral value ΣACRq of the q-axis current control as the voltage error Δv, excluding the component equivalent to the voltage drop due to the coil resistance R. A dominating factor of the voltage error Δv is an induced voltage constant error Δ, which is dependent on the rotational velocity of the motor 3. In view of this, the voltage error compensation section 33 obtains the induced voltage constant error Δφ based on the voltage error Δv.
The voltage error compensation section 33 determines whether a voltage saturation has occurred, and only in the state of no voltage saturation, executes voltage error compensation processing. Whether a voltage saturation has occurred is determined based on Formula (20). Specifically, the voltage error compensation section 33 executes the voltage error compensation processing when Kh is less than v1
Based on the q-axis current iq
The coefficient multiplication section 71 multiplies the coil resistance R by the q-axis current iq
The absolute value arithmetic section 73 performs an arithmetic operation to obtain the absolute value of the electrical angular velocity ω of the motor 3, and outputs the absolute value to the limiter 74. When the absolute value of the electrical angular velocity co reaches a limit value set in advance, the limiter 74 limits the absolute value of the electrical angular velocity ω to the limit value. For example, the limiter 74 limits the lower limit of the absolute value of the electrical angular velocity ω to 10 Hz×2π, and limits the upper limit of the absolute value of the electrical angular velocity ω to 100 Hz×2π. The sign function arithmetic section 75 performs an arithmetic operation to obtain the positivity or negativity of the electrical angular velocity ω of the motor 3, and outputs the arithmetic result to the multiplication section 76. The multiplication section 76 multiplies the output of the limiter 74 and the output of the sign function arithmetic section 75, and outputs the product to the division section 77. An example of the sign function processing is to output “+1” when the input is a positive, and to output “−1” when the input is a negative value.
The division section 77 divides the output of the subtraction section 72 by the output of the multiplication section 76, thereby obtaining the induced voltage constant error Δφ. A dominating factor of the voltage error Δv is the induced voltage constant error Δφ. Another voltage error factor is the voltage component, which is dependent on current differentiation. The voltage component is not easy to identify and is negligible in stationary state. In view of this, the voltage component is set to be removed by the LPF 78.
Specifically, the LPF 78 removes a high-harmonic component of the induced voltage constant error Δφ output from the division section 77, thereby generating the induced voltage constant error ΔφLPF, and outputs the induced voltage constant error ΔφLPF. An example of the LPF 78 is a primary-delay filter with a cutoff frequency adjustable as a setting parameter.
Referring back to
The absolute value arithmetic section 80 performs an arithmetic operation to obtain the absolute value of the electrical angular velocity ω of the motor 3, and outputs the absolute value to the limiter 81. When the absolute value of the electrical angular velocity co reaches a limit value set in advance, the limiter 81 limits the absolute value of the electrical angular velocity ω to the limit value. For example, the limiter 81 limits the lower limit of the absolute value of the electrical angular velocity ω to 10 Hz×2π, and limits the upper limit of the absolute value of the electrical angular velocity ω to 100 Hz×2π. This inhibits occurrences at the time when the motor 3 is at a super-low velocity, such as the division section 77 dividing by zero and the output of the division section 77 becoming excessive.
The sign function arithmetic section 83 performs sign function processing of the electrical angular velocity ω of the motor 3, and outputs the arithmetic result to the multiplication sections 82 and 85. The multiplication section 82 multiplies the output of the limiter 81 and the output of the sign function arithmetic section 83, and outputs the product to the division section 87.
The coefficient multiplication section 84 multiplies the q-axis current iq
The division section 87 divides the output of the subtraction section 86 by the output of the multiplication section 82, and outputs the division result to the subtraction section 88. The subtraction section 88 subtracts the induced voltage constant error ΔφLPF from the output of the division section 87, and outputs the subtraction result as an adjustment value φCPC
The integral value ΣACRq of the q-axis current deviation in the state of voltage saturation is considered to contain a voltage component represented by Formula (21). In Formula (21), ΔLd denotes a parameter error of d-axis inductance.
ΣACRq=Riq
In Formula (21), the first item on the right-hand side denotes the amount equivalent to the voltage drop due to resistance, and the second item on the right-hand side denotes the amount equivalent to the parameter error of induced voltage and inductance. Additionally, the third item on the right-hand side (ΔVst) denotes an item representing the voltage saturation.
The second item on the right-hand side is an item dependent on the rotational velocity of the motor 3, and as such, is compensated for by the induced voltage constant error ΔφLPF (which is a value obtained by an arithmetic operation in the state of no voltage saturation) output from the voltage error compensation section 33.
When the input adjustment value φCPC
Referring back to
The adjustment value φCPC
The output of the coefficient multiplication section 91 is integrated by the integral section 93 and input into the limiter 94. The limiter 94 limits the output of the integral section 93 within a predetermined range, and outputs the limited output to the addition section 96. Specifically, when the output of the integral section 93 reaches an upper limit or a lower limit set in advance, the limiter 94 limits the output of the integral section 93 to the upper limit or the lower limit, and outputs the limited output.
The output of the coefficient multiplication section 92 is differentiated by the differential section 95 and output to the addition section 96. The addition section 96 adds together the output of the coefficient multiplication section 90, the output of the limiter 94, and the output of the differential section 95, and outputs the sum to the limiter 97. The limiter 97 limits the correction voltage command Δvd* to keep the correction voltage command Δvd* from exceeding a predetermined range.
Referring back to
The amplitude command generation section 35 obtains an amplitude v1* of the output voltage command based on the q-axis voltage command vq** and the d-axis voltage command vd**. For example, the amplitude command generation section 35 obtains the amplitude v1* of the output voltage command from Formula (22), and outputs the amplitude v1* to the limiter 37. The limiter 37 limits the amplitude v1* of the output voltage command within a predetermined range, and outputs the limited amplitude v1*.
v1*=(vd**2+vq**2)1/2 (22)
The phase command generation section 36 obtains an output phase command θa* (phase difference as compared with the d axis) based on the q-axis voltage command vq** and the d-axis voltage command vd**. For example, the phase command generation section 36 obtains a phase command θa* of the output voltage from Formula (23), and outputs the phase command θa* to the addition section 38. The addition section 38 adds the rotor electrical angle phase θ detected by the position detection section 4 to the phase command θa of the output voltage, thereby generating an output phase command θ*, and outputs the output phase command θ* to the power conversion section 10.
θa*=tan−1(vq**/vd**) (23)
Based on the amplitude v1* of the output voltage command output from the vector control section 12 and based on the phase command θ* of the output voltage, the PWM signal generation section 14 of the power conversion section 10 generates a PWM signal by known PWM control so as to control the three-phase inverter circuit 13.
Thus, the motor controller 1 according to the first embodiment includes the constant output control section 34, which generates the correction voltage command Δvd* based on the integral value ΣACRq of the q-axis current deviation output from the q-axis integral control section 52 and which outputs the correction voltage command Δvd*. Then, the motor controller 1 subtracts the correction voltage command Δvd* from the d-axis voltage command vd**′, thereby obtaining the d-axis voltage command vd**. This inhibits both degradation of torque responsivity and degradation of the voltage utilization ratio in the constant output region of the motor 3.
Next, a motor controller according to the second embodiment will be described. The motor controller according to this embodiment is different from the motor controller 1 according to the first embodiment in that the vector control section includes a current limitation section. The elements with corresponding or identical functions to those of the elements of the motor controller 1 according to the first embodiment are assigned identical reference numerals, and these elements will not be elaborated here.
As described above, during the constant output control, the control of the q-axis current iq
im=√{square root over (id
In view of this, the vector control section 12A of the motor controller 1A includes the current limitation section 23 and the subtraction section 25. This configuration controls the output current im to turn the driving state toward alleviated voltage saturation. Specifically, when the output current im exceeds a current limitation command im*, which is a limit value, then the speed command is lowered (for example, the acceleration-deceleration rate of the speed command or the speed command value is lowered in accordance with a current deviation Δim* between the current limitation command im* and the output current im, or the acceleration-deceleration rate of the speed command or the speed command value is lowered in accordance with a q-axis current amend command Δiq*). At the same time, the q-axis current command iq* is also lowered. This inhibits overcurrent as early as possible while controlling the driving state toward alleviated voltage saturation. The speed command is a value that is proportional to the output frequency of the voltage command. When the current flowing through the motor 3 exceeds the limit value, the current limitation section 23 lowers, for example, the output frequency of the voltage command or the acceleration-deceleration rate of the output frequency.
In the current limitation section 23, with a transfer function Gc(s) defined between the output current deviation Δim* and the q-axis current amend command Δiq*, a block diagram ranging from the current limitation command im* to the output current im is established as shown in
The output current deviation Δim* is a deviation between the current limitation command im* and the output current im, and the current limitation command im* is generated dynamically within the current limitation section 23.
Here, the control block diagram shown in
Q=V×I=|V∥I| sin θφ=|V∥I|−ω×T (25)
T=Δiq*×Kt=ΔTlim* (26)
Also considering that the driving is in the constant output state, two parameters (V and ω) are set as fixed values as represented by Formula (27), and the transfer function from the torque T to the current |I| can be represented by Formula (28).
Hence, ensuring that the transfer function Gc(s) between the output current deviation Δim* and the q-axis current amend command Δiq* is an integral characteristic realizes stable control of the output current im. Thus, the arithmetic operator of the transfer function Gc(s) can be represented as an integral control section by Formula (29).
In Formula (29), ωAIC denotes control response [rad/s] of the current limitation section 23, and Vdc denotes DC bus line voltage [v (volts)]. It is noted that ωAIC is a setting parameter. From Formula (29), an integral gain of the transfer function Gc(s) can be represented by Formula (30).
Thus, the current limitation section 23 multiplies the output current deviation Δim*, which is a deviation between the current limitation command im* and the output current im, by an integral gain Ki, which is shown in Formula (30), and integrates the product, thereby generating the q-axis current amend command Δiq*. Then, as described above, the current limitation section 23 lowers the speed command value co in accordance with the output current deviation Δim*, while at the same time performing q-axis current control based on the q-axis current command iq*, which results from subtraction of the q-axis current amend command Δiq* by the subtraction section 25. This inhibits great changes in the output current im as early as possible.
The current limitation section 23 is thus configured, and the integral gain Ki of the current limitation section 23 changes in accordance with the power factor θφ. A specific operation is that when the power factor θφ is high, the gain becomes low, while when the power factor θφ is small, the gain becomes high. In
In both the first embodiment and the second embodiment, the position detection section 4 is used to detect the rotor electrical angle phase θ of the motor 3. This, however, should not be construed in a limiting sense. A configuration without the position detection section is also possible, in which case an arithmetic operation is performed for the rotor electrical angle phase θ based on the output current and the output voltage.
Obviously, numerous modifications and variations of the present disclosure are possible in light of the above teachings. It is therefore to be understood that within the scope of the appended claims, the present disclosure may be practiced otherwise than as specifically described herein.
Number | Date | Country | Kind |
---|---|---|---|
2012-243949 | Nov 2012 | JP | national |
Number | Name | Date | Kind |
---|---|---|---|
5689166 | Nagayama et al. | Nov 1997 | A |
5877607 | Masaki et al. | Mar 1999 | A |
20030090231 | Yoshimoto et al. | May 2003 | A1 |
20040183496 | Yoshimoto | Sep 2004 | A1 |
20050264255 | Maruyama et al. | Dec 2005 | A1 |
20060061310 | Takai et al. | Mar 2006 | A1 |
20070222409 | Kariatsumari et al. | Sep 2007 | A1 |
20090128081 | Yoneda et al. | May 2009 | A1 |
20110031907 | Takahashi | Feb 2011 | A1 |
20110080125 | Shimada et al. | Apr 2011 | A1 |
20120194113 | Lee et al. | Aug 2012 | A1 |
Number | Date | Country |
---|---|---|
2010-022165 | Jan 2010 | JP |
Number | Date | Country | |
---|---|---|---|
20140125261 A1 | May 2014 | US |