Method of controlling operating speed of a loom

Information

  • Patent Grant
  • 5060161
  • Patent Number
    5,060,161
  • Date Filed
    Wednesday, November 29, 1989
    35 years ago
  • Date Issued
    Tuesday, October 22, 1991
    33 years ago
Abstract
A method of controlling an operating speed of the loom for increasing the production rate by varying the present operating speed so as to be increased includes a step of estimating a stop rate at the final point of time during a predetermined period at a certain time during the predetermined period, and comparing the estimated stop rate with a predetermined limit stop rate, and setting the present operating speed to increase it when the estimated stop rate is less than the limit stop rate and to decrease it when the estimated stop rate exceeds the limit stop rate.
Description

BACKGROUND OF THE INVENTION
1. Field of the Invention:
The present invention relates to a method of controlling an operating speed of a loom in an optimum condition.
2. Prior Art:
The applicant proposed an idea which was described in U.S. Ser. No. 324,536 and European Patent Application No. 89 104 577.5 to increase the operating speed of the loom taking into account a weaving condition during operation of the loom to thereby increase the production rate as much as possible.
The method comprises the steps of judging the quality of a fabric with reference to past data during operation of the loom and setting the operating speed of the loom to increase it or to decrease it on the basis of the result of the judgment of the quality of the fabric. During the controlling of the operating speed of the loom, when the fabric having a deteriorated quality is woven or the operator can not cope with or share a plurality of stopped looms during the controlling of the operating speed of the loom, there occurs a case that the operating condition of the loom is not placed in the optimum weaving operation during the collecting of the past data.
Hence, it is an important factor in the case of the controlling the operating speed of the loom to estimate the operating condition of the loom while the loom is kept in the optimum operating condition.
Accordingly, the object of the present invention is to provide a method of controlling the operating speed of the loom comprising the steps of estimating a stop rate of the loom taking into account a predetermined quality of the fabric and an operator's sufficient time for sharing the loom (hereinafter referred to as an operator's sufficient time) on the basis of the present data during the operation of the loom and setting the operating speed of the loom to the optimum operating condition on the basis of the estimated result.
In view of the object of the present invention, the method comprises the steps of: estimating a stop rate of the loom at the final point of time, the estimation being carried out at a certain point of time during a predetermined period while the operating speed of the loom is controlled to be increased so that the production rate is increased; comparing the estimated stop rate with a predetermined limit stop rate so that the operating speed of the loom is increased when the estimated stop rate is less than the predetermined limit stop rate and the operating speed of the loom is decreased when the estimated stop rate exceeds the predetermined limit stop rate. The limit stop rate is previously determined at least from the quality of the fabric or the operator's sufficient time. The quality of the fabric is generally determined on the basis of stop rate of the loom per length of the fabric (hereinafter referred to as the downtime rate of the loom). That is, as the downtime rate of the loom is increase, the stepped portion is increased in the fabric to thereby deteriorate the quality of the fabric. Accordingly, it is judged that the less the downtime rate of the loom, the better the quality of the fabric. The operator's sufficient time is determined by an operating rate of the loom. That is, the downtime of the loom comprises a waiting time, namely, the arrival time of the operator to the loom and a time needed by the operator for sharing the repair of the loom after the arrival of the operator. If the operator's time is sufficient, the waiting time is increased and the downtime of the loom is increased, whereby the operating rate of the loom is reduced. Accordingly, it is judged that the operator's sufficient time is increased as the operating rate is increased.
The method of controlling the operating speed of the loom comprises the steps of: controlling the stop rate on the basis of the estimated stop rate during the change of operating speed of the loom so that the fabric is prevented from being woven inferiorly during the step of controlling the operating speed of the loom or the operating speed of the loom is prevented from being increased at the state where there is insufficient operator's time.
The above and other objects, features and advantages of the present invention will become more apparent from the following description taken in conjunction with the accompanying drawings.





BRIEF DESCRIPTION OF THE DRAWINGS
FIGS. 1 and 2 are graphs of assistance in explaining a relationship between a shifting time and operation and stoppage of the loom;
FIG. 3 is a graph of assistance in explaining a increment of an operating speed of the loom and an increment of the stop rate relative to the increase of a production rate;
FIG. 4 is a graph of assistance in explaining the increment of the operating speed of the loom and the increment of the stop rate relative to an operating rate;
FIG. 5 is a graph of assistance in explaining the increment of operating speed of the loom and the increment of the stop rate relative to a downtime rate of the loom;
FIG. 6 is a graph of assistance in explaining the increment of the operating speed of the loom and the increment of the stop rate relative to the operating rate, the downtime rate of the loom and the production rate;
FIG. 7 is a graph of assistance in explaining a weighting relative to an elapsed time;
FIG. 8 is a view of assistance in explaining an example of result made by a computer simulation;
FIG. 9 is a block diagram showing a control system; and
FIG. 10 is a flowchart showing a program for controlling the operating speed of the loom.





DESCRIPTION OF A PREFERRED EMBODIMENT OF THE INVENTION
Assume that the result of of i times of shift during operation of the loom has the following data.
Shifting time: T.sub.0 (min)
Operating speed: N.sub.i (RPM)
Stop rate: n.sub.i (time/shift)
Average downtime: .tau. (min)
Production rate: P.sub.i (cmpx)
Operating rate: E.sub.i (%)
Downtime rate of the loom: S.sub.i (stop rate/cmpx)
where cmpx is a unit representative of 0.1 million picks.
The production rate P.sub.i, the operating rate E.sub.i and the downtime rate of the loom S.sub.i are respectively determined as follows with reference to FIG. 1.
FIG. 1 shows the operation and stoppage of the loom during i times of shift in simplicity and the hatched portion corresponds to the production rate P.sub.i. ##EQU1## where
A=100000
T.sub.0 >n.sub.i .tau.
The downtime rate of the loom S.sub.i becomes a factor for determining the quality of the fabric.
Assuming that there is the operator's sufficient time and the average stoppage time .tau. is not varied provided that the stop rate n.sub.i increased for the increment y when the operating speed of the loom is increased for the increment of x, the production rate P.sub.i+1, the operting rate E.sub.i+1 and the downtime of the loom during i+1 times of shift S.sub.i+1 are expressed as follows. ##EQU2## where
N.sub.i +x>0 .thrfore.x>-N.sub.i
n.sub.i +y>0.thrfore.y>-n.sub.i
T.sub.0 -(n.sub.i +y).tau.>0 .thrfore.y<(T.sub.0 -.tau.n.sub.i).tau.
FIG. 2 shows the operation and stoppage of the loom during i+1 times of shift in simplicity. It is evident from FIG. 2 that with increase of the operating speed of the loom production rate is likely to increase but the production rate decreases with increase of the stop rate. A condition where the production rate increases during (i+1) times of shift where the operating speed of the loom is increased.
A difference Z of the production rate for the period between the i times of shift and (i+1) times of shift is expressed as follows. ##EQU3##
To meet Z>0, the following expression is to be established.
.tau.(N.sub.i +x)y<(T.sub.0 -n.sub.i .tau.)x
where, .tau.>0, N.sub.i +x>0, hence the above expression is expressed as follows. ##EQU4##
FIG. 3 shows this portion in the hatched portion. A fourth quadrant in the same figure is omitted as out of scope since the stop rate is generally increased with increase of the operating speed of the loom.
In fact since the stop rate ni, namely, y is a positive integer supposed that the increment y=k (k is a positive integer) of the stop rate is established if the operating speed of the loom N.sub.i is increased for the increment of x, the following expression is established. ##EQU5##
k.tau.x+k.tau.N.sub.i <(T.sub.0 -n.sub.i .tau.)x
[T.sub.o -(n.sub.i +k).tau.]x>k.tau.N.sub.i
Since T.sub.o -(n.sub.i +k)>0 ##EQU6##
That is, when the stop rate n.sub.i is increased for the increment of k, the increment of the operating speed of the loom is necessary to be greater than k.tau.N.sub.i /[T.sub.o -(n.sub.i +k).tau.] to increase the production rate P.sub.i. In other word, the production rate P is increased if the increment of the stop rate n.sub.i is less than k when the increment of the operating speed of the loom is increased to be greater than k.tau.N.sub.i /(T.sub.o -(n.sub.i +k).eta.).
For example, assuming that k=2 and the operating speed N.sub.i+1 of the loom is increased to become N.sub.i+1 =2.tau.N.sub.i /[T.sub.o -(n.sub.i+2).tau.]+N.sub.i, the stop rate N.sub.i+1 <n.sub.i +2, namely, inasmuch as the n.sub.i and n.sub.i+ 1 are positive integers the production rate P will be increased if the expression n.sub.i+1 .ltoreq.n.sub.i+1 is established.
Studying the operator's sufficient time with reference to the operating rate, if the operating rate is greater than E.sub.o, the operator can work with sufficient time and the following expressions are to be established provided that the average downtime is not varied.
E.sub.i .gtoreq.E.sub.o
E.sub.i+1 .gtoreq.E.sub.o
Accordingly, the following expression is established. ##EQU7## from the expression (1) ##EQU8##
The n.sub.i is defined as the limit stop rate Q.sub.E during one shift in view of the operator's sufficient time. From the expression (2)
100T.sub.o -100n.sub.i .tau.-100.tau.y.gtoreq.E.sub.o T.sub.o
100.tau.y.ltoreq.T.sub.o (100-E.sub.o)-100n.sub.i .tau.
Hence, the increment y of the stop rate for satisfying the expression E.sub.i+1 .gtoreq.E.sub.o is to satisfy the following expression. ##EQU9##
FIG. 4 shows the area satisfying the expression as the hatched portions.
Studying the quality of the fabric with reference to the downtime rate of the loom, provided that the standard that the fabric stands the test of the quality is determined if the downtime rate of the loom (stop rate/cmpx) is less than S.sub.o, the following expressions are to be established. ##EQU10## from the expression (3) ##EQU11## The n.sub.i is to be defined as the limit stop rate Q.sub.S during one shift of the loom relative to the quality of the fabric.
From the expression (4)
An.sub.i +Ay.ltoreq.S.sub.o (N.sub.i+x)(T.sub.o -.tau.n.sub.i)-S.sub.o (N.sub.i+x).tau.y
(S.sub.o .tau.(N.sub.i+x)+A)y.ltoreq.S.sub.o (T.sub.o -.tau.n.sub.i)(x+N.sub.i)-An.sub.i
Accordingly, the increment y of the stop rate for satisfying the expression S.sub.i+1 .ltoreq.S.sub.o is to satisfy the following expression. ##EQU12##
FIG. 5 shows an area satisfying the expression just above. The following expression will be established provided that y=k (k is a positive integer). ##EQU13##
kS.sub.o .tau.x+k(S.sub.o .tau.N.sub.i +A).ltoreq.S.sub.o (T.sub.o -.tau.n.sub.i)x+S.sub.o (T.sub.o -.tau.n.sub.i)N.sub.i A.sub.i
S.sub.o (k.tau.-T.sub.o +.tau.n.sub.i)x.ltoreq.S.sub.o (T.sub.o -.tau.n.sub.i)N.sub.i -k(S.sub.o .tau.N.sub.i +A)-An.sub.i
From (k+n.sub.i).tau.-T.sub.o <O ##EQU14##
That is, when the stop rate is increased for the increment of k times, no quality problem will occur if the increment x of the operating rate satisfies the expression set forth just above. In other words, no quality problem will occur if the operating speed is increased more than x establised just above provided that the increment of the stop rate is less than k time.
FIG. 6 is a single graph representing the combination of graphs of FIGS. 3, 4 and 5. Evident from FIG. 6 is the the area (hatched portion) where the conditions of both the operator's sufficient time (operating rate E.sub.o) and the quality of the fabric (stoppage level S.sub.o) are satisfied and the production rate P.sub.i is increased.
Generally the stop rate is increased as the operating speed is increased. The stop rate can be represented by, for example, straight lines L.sub.1 and L.sub.2 as illustrated in FIG. 6, provided that the stop rate is proportional to the increment of the operating speed of the loom although it is a positive integer.
The straight line L.sub.1 is deviated from the area defined by a curved line y.sub.p in case of x>0 and within the same area in case of x<0. When the stop rate is varied greatly, accompanyied by the variation of the operating speed, the production rate is not, increased even if the operating speed is increased; rather, the production rate is expected to be increased when the operating speed is decreased.
The straight line L.sub.2 is within the area defined by the curved line y.sub.p in case of x>0 and is out of the same area in case of x<0. In such case, the production rate is increased when the operating speed is increased but the production rate is not increased even if the operating speed is decreased.
In the case represented by the straight line contacting the curved line y.sub.p at the origin, the production rate will be decreased when the operating speed is increased or decreased and the same straight line becomes the optimum point of the operating speed in order to maximize the production rate.
Although the theory can be applicable to the case set forth above, the ratio of increase and decrease of the operating speed to those of the stop rate can not be determined since the straight lines L.sub.1 and L.sub.2 are practically not known.
Hence, the above tendency is estimated with reference to the algorithm set forth hereunder to control the operating speed.
The followings are definitions of estimating controls, judgment based on the estimated control and prosecution made by the judgment assuming that the production rate is P at the elapsed time t=T, the stop rate is n, the operating time is r, and the operating speed is N.
(1) Estimation of Parameter
Provided that the result of each parameter at present, namely, at the elapsed time t=T is kept advanced as it is until reaching to a cetain elapsed time t=T.sub.o which is a completion time of one shift, the production rate P, the stop rate n, the operating rate E and the downtime rate of the loom S are estimated as follows. Estimated values are marked at .
P=T.sub.o P/T
n=T.sub.o n/T
E=100AP/T.sub.o N
S=n/P
A limited stop rate Q.sub.S satisfying a limited downtime rate of the loom S.sub.o and during one shift and a limited stop rate Q.sub.E satisfying an operating rate during one shift are expressed as follows from the explanation set forth above.
Q.sub.S =S.sub.o T.sub.o N/(A+S.sub.o .tau.N)
Q.sub.E =T.sub.o (100-E.sub.o)/100.tau.(where .tau..noteq.0)
.tau.=(T-r)/n
(2) Setting of Weighting Coefficient (Refer to FIG. 7)
Inasmuch as the certainty of the estimated value will be increased as the elapsed time t=T draws to close to T.sub.o, the weighting coefficient value can be set, for example, to as follows corresponding to the elapsed time.
w=0 in case of T<T.sub.o /4
w=(2T/T.sub.o)-0.5 in case of T.sub.o /4.ltoreq.T<3T.sub.o /4
w=1 in case of T.ltoreq.3T.sub.o /4
(3) Judgement
(a) Operating speed down in case of n>Min(Q.sub.S, Q.sub.E)
When the estimated stop rate n exceeds one of less values of Q.sub.S or Q.sub.E, it will be changed in the direction to decrease the operating speed. At this time, a variation ratio .DELTA.N is, for example, determined in the following manner.
The increment Y of a desired stop rate after variation of the operating speed is expressed as follows.
Y=(Min(Q.sub.S, Q.sub.E)-n)/2Y<0
The result shows Y<0 which represents the decrement of the stop rate. The above expression shows, to prevent excessive control, that the increment Y of the stop rate after variation of the operating speed is set to be half of the difference between the limit stop rate and the estimated stop rate as an example. However, the increment Y for a provisional stop rate can be varied to be decreased or increased by allowing the denominator values to be greater or less than 2.
The increment X of the operating speed satisfying the increment Y can be determined from the following expression using the relationship between the increment x of the operating speed and the increment k of the stop rate.
x=Y.multidot..tau..multidot.N/(T.sub.o -(n+Y).tau.)
A discriminant expression W is defined as follows using the weighting coefficient w corresponding to the elapsed time and the weighting coefficient W.sub.N of the increment X of the operating speed determined by the expression just set forth above can be determined as follows.
Suppose that W=w[n-Min(Q.sub.S, Q.sub.E)].
W.sub.N =0 in case of W.ltoreq.0
W.sub.N =W in case of 0<W<1
W.sub.N =1 in case of W.gtoreq.1
The variation ratio .DELTA. N of the operating speed can be expressed as follows from the thus determined increment X and the weighting coefficient W.sub.N.
.DELTA.N=W.sub.N .multidot.X
(b) Operating speed up in case of n.ltoreq.Min(Q.sub.S, Q.sub.E)
When the estimated stop rate n is less than the limit stop rate Q.sub.S or the limit stop rate Q.sub.E, it will be varied in the direction so as to increase the operating speed. At this time, the variation ratio .DELTA. N can be determined, for example, as follows in the same way as (a).
Y=(Min(Q.sub.S, Q.sub.E)-n)/2 Y.gtoreq.0
x=Y.multidot..tau..multidot.N/(T.sub.O -(n+Y).tau.)
Suppose that W=w(n-Min(Q.sub.S,Q.sub.E)).
W.sub.N =0 in case of 0<W
W.sub.N =.vertline.W.vertline. in case of -1<W.ltoreq.0
W.sub.N =1 in case of W.ltoreq.-1
N=W.sub.N .multidot.X
(c) Suppose that the maximum ratio of the variation ratio .DELTA. N of the operating speed is .DELTA. N.sub.o and .DELTA.N=.DELTA.N.sub.o in case of .vertline..DELTA.N.vertline.>.DELTA.N.sub.o.
(d) Suppose that the upper limit and the lower limit of the variation range of the operating speed are respectively N.sub.max, N.sub.min.
N=N.sub.min in case of N.sub.min >N+.DELTA.N
N=N.sub.max in case of N.sub.max <N+.DELTA.N
N=N+.DELTA.N in case of N.sub.min .ltoreq.N+.DELTA.N.ltoreq.N.sub.max
The operating speed is to be varied in the course of shifting loom as set forth above so that the stop rate at the time of completion of one shift does not exceed the limit stop rates Q.sub.S and Q.sub.E. The control method is confirmed by the computer simulation which results in the following. The computer simulation was carried out under the following steps.
(1) Generation of downtime
The situation of generation of downtime is given by a Poisson distribution. Provided that the number of .lambda. of the downtime per unit time is generated, the probability of stoppage of the loom for the number of s times during the interval of T can be determined by the following expression. ##EQU15## where .lambda.T is an average downtime of the loom for the interval of T. The exponential distribution can be expressed as follows if the time interval per unit time is represented by the distribution type.
P(t)=.lambda.e.sup.-.lambda.T
where 1/.lambda. is an average value of the distribution.
(2) Service time
A service time (downtime of the loom) is given by the index distribution. The service time is expressed as follows, assuming that an average service ratio (average number of services per unit time) is .mu..
P(t)=.mu.e.sup.-.mu.T
where 1/.mu. is an average value of the distribution which accords to the average service time.
(3) In the case of a weaving mill
T.sub.o : shifting time
n: stop rate in total
N: operating speed
R: service time/stop
Suppose that the data set forth above is given, the following definitions are expressed. ##EQU16## where each distribution for each item is expressed as follows. ##EQU17##
FIG. 8 shows an example of the result of the computer simulation supposing that the limit operating rate E.sub.o =80.0(%) and the limit downtime rate S.sub.o =5.0 (stop rate/cmpx).
The simulation operation is effected while the control according to the present invention is turned off during the first 1 to 9 shifts. The simulation operation is effected while the control according to the present invention is turned on during next 10 to 18 shifts.
Inasmuch as each data is varied for each shift it is difficult to find out the variations of the operating rate, the downtime rate of the loom, the limit operating rate and the production rate. However, the stop rate n is increased from 5 to 11, the operating rate E is decreased from 9.6 to 9.0, the downtime rate of the loom S is increased from 1.45 to 3.18 and the production rate P is increased from 3.23 cmpx to 3.42 cmpx. This means that the production rate P is inceased 1.74 cmpx during the total of 9 shifts.
Although the operating rate and the downtime rate are both deteriorated, it is evident that the limit operating rate (STD Eff) does not exceed 8.0% and the limit downtime rate of the loom (STD stop rate/cmpx) does not exceed 5.0.
According to the present invention, the stop rate at the final point of time during the predetermined period is previously estimated and the operating speed is increased or decreased not to exceed the limit stop rate determined from the estimated value and the quality of the fabric, namely, from the downtime rate and the operator's sufficient time, namely from the limit stop rate determined from the operating rate. Hence, the production rate can be increased so as to be as high as possible since the predetermined set quality of the fabric or the operator's sufficient time are satisfied at the time of completion of the predetermined period and the operating speed can be controlled on the basis of the prospective estimation.
In case the operating speed is to be varied every time the fabric is stopped, the inferior influence on the quality of the fabric caused by the variation of operating speed during the weaving operation is eliminated. Furthermore, if the estimated value is weighted, the greater variation at the early shifting period having less data does not occur, thereby stabilizing the operating speed.
FIG. 9 shows a system when the method of controlling the operating speed according to the present invention is carried out.
A host computer 1 having stored therein a program of the method of controlling the operating speed is connected to computers 3 for controlling a plurality of looms 2 (hereinafter refered to as control computers) via a data line 4. The program stored in the host computer 1 specifies the loom 2 to be controlled within a predetermined time (set more than one time) during the predetermined period to thereby execute the program of the method of controlling the operating speed. The control can be carried out every time that the loom 2 is stopped or at a predetermined period.
FIG. 10 is a flowchart carrying out the method of contolling the operating speed.
In a first step, after starting the program, a determination is made as to whether or not a shift change has occurred. If the shift change has occurred the program goes to a second stop where the elapsed time is set to 0 and goes to a third step where the execution time for controlling the loop 2 is judged. If the shift change has not occurred at the first step, the program jumps the second step and goes to the third step. The execution time in the step 3 is set after a lapse of the predetermined time or the stop of the loom as mentioned above. The host computer 1 executes an estimation of parameters relative to the loom 2 supposing that the shift change has been completed in a fourth step when the predetermined time has elapsed after the completion of execution of the control program or when the loom 2 is stopped, thereby executing an arithmetic operation of data necessary for the parameter. The host computer 1 calculates the weighting coefficient w corresponding to the lapsed time from the shift change in a fifth step and then determines if determined estimated stop rate n is within the limit stop rate Q.sub.S or Q.sub.E in a sixth step. If the estimated stop rate n is within the limit stop rate Q.sub.S or Q.sub.E, the host computer 1 calculates the amount that the operating speed is to be incremented in a seventh step. However, if the estimated stop rate n is not within the limit stop rate Q.sub.S or Q.sub.E, the host computer 1 calculates the operating speed to be decremented in an eighth step. In a ninth step, the operating speed after the shift change is calculated and thereafter, the host computer 1 feeds the control computer 3 a new operating speed in a tenth step to thereby change the operating speed of the corresponding loom 2. After a series of programs are executed, the same control in one shift is repeated for the next execution time. In such a manner, the operating speed of the loom 2 to be controlled is set to increase the production rate as much as possible within the limited stop rate determined from the downtime rate of the loom or the opeating rate and within the operator's sufficient time. The variation rate of the operating speed is limited so as not to exceed the predetermined value. The operating speed is changed within a predetermined maximum or the minimum operating speed.
In case the host computer is used only for collecting the data of the operation of a plurality of looms 2 and for the storing thereof, the method of controlling the operating speed according to the present invention is respectively executed by the control computer 3 of the loom 2.
Although the invention has been described in its preferred form with a certain degree of particularlity, it is to be understood that many variations and changes are possible in the invention without departing from the scope thereof.
Claims
  • 1. A method of controlling an operating speed of a loom comprising the steps of:
  • estimating a stop rate of the loom at a final point of time during a predetermined period of a weaving operation, the estimation being carried out at a certain point of time during the predetermined period; and
  • comparing the estimated stop rate with a predetermined limit stop rate and setting the present operating speed of the loom to increase it when the estimated stop rate is less than the limit stop rate and setting the present operating speed of the loom to decrease it when the estimated stop rate exceeds the limit stop rate.
  • 2. A method of controlling an operating speed of a loom according to claim 1, wherein the limit stop rate is determined by one of a predetermined quality of the fabric and an operator's sufficient time for sharing the loom.
  • 3. A method of controlling an operating speed of a loom according to claim 1, wherein a weighting coefficient calculation operation is carried out corresponding to the lapsed time from a first point of time of the predetermined period unit the final point of time of the same predetermined period during the execution of an arithmetic, the weighting coefficient being used to calculate the change the present operating speed of the loom on the basis of the estimated stop rate.
  • 4. A method of controlling an operating speed of a loom according to claim 2, wherein the quality of the fabric is judged from a downtime rate of the loom.
  • 5. A method of controlling an operating speed of a loom according to claim 2, wherein the operator's sufficient time for sharing the loom is determined from an operating rate of the loom.
Priority Claims (1)
Number Date Country Kind
63-303183 Nov 1988 JPX
US Referenced Citations (10)
Number Name Date Kind
4494203 Suzuki et al. Jan 1985
4525787 Tzvetkov Jun 1985
4529012 Sainen et al. May 1985
4553569 Kimbara et al. Nov 1985
4556945 Fry Dec 1985
4605044 Sakano Aug 1986
4715411 Van Bogaert et al. Dec 1987
4736324 Sainen et al. Dec 1990
4884597 Tamura et al. Dec 1989
4893250 Sainen Jan 1990