Information
-
Patent Grant
-
6435006
-
Patent Number
6,435,006
-
Date Filed
Friday, June 25, 199925 years ago
-
Date Issued
Tuesday, August 20, 200222 years ago
-
Inventors
-
Original Assignees
-
Examiners
Agents
- Townsend and Townsend and Crew LLP
-
CPC
-
US Classifications
Field of Search
US
- 073 37
- 073 375
- 033 837
- 033 655
- 033 517
-
International Classifications
-
Abstract
The invention relates to a method for measuring lengths pneumatically. According to the invention, the medium used for measuring is conducted through a front nozzle and a measuring nozzle, onto a deflection plate. The distance of the deflection plate from the measuring nozzle is determined from the change in pressure p2 occurring between the front nozzle and the measuring nozzle. The pressure p1 in front of the front nozzle is also measured and the length value resulting from the pressure measurement between the front nozzle and the measuring nozzle is corrected according to said pressure p1. This ensures that a high level of measuring accuracy can be achieved using a small quantity of compressed air.
Description
BACKGROUND OF THE INVENTION
1. Field of the Invention
The present invention relates to a method for pneumatic length measurement in which the medium used for the measurement is directed through a pre-nozzle and a measuring nozzle onto an impact plate and the spacing of the impact plate from the measurement nozzle is determined from a change of the pressure p
2
between the pre-nozzle and the measurement nozzle.
2. Description of the Prior Art
Methods for pneumatic length measurement have been known for a long time and are in particular used in the manufacture of workpieces because of the high measurement accuracy, the rapid determination of the measured value and the insensitivity of the measurement method with respect to external influences.
A disadvantage of this method lies in the fact that the supply pressure of the measurement apparatus must be constant in order to be able to guarantee an adequate measurement accuracy. So-called bridge circuits are thus also known, in which the pressure is measured between a measurement branch and a reference branch. Because air pressure fluctuations of the supply pressure act equally on the reference branch and on the measurement branch, these fluctuations are compensated. The disadvantage of the bridge circuit, however, lies in the fact that twice the quantity of air is consumed through the reference branch. Moreover, the manufacturing costs and the complexity of the adjustment are increased by the reference branch.
SUMMARY OF THE INVENTION
The invention is based on the object of setting forth a pneumatic method of length measurement of the initially named kind, which has a high measurement accuracy despite fluctuations of the supply pressure and which, at the same time, enables a low consumption of compressed air and low manufacturing costs.
This object is solved in that the pressure p
1
in front of the pre-nozzle is additionally measured and the length value resulting from the pressure measurement between the pre-nozzle and the measurement nozzle is corrected in dependence on the pressure p
1
.
In accordance with the invention a correction of the pressure measured between the pre-nozzle and the measurement nozzle of the measurement branch, i.e. of the length value derived therefrom, is effected instead of a bridge circuit. The correction takes place in this respect in dependence on the additionally measured pressure p
1
before the pre-nozzle. This additionally measured pressure p
1
is a measure of the fluctuations of the air supply pressure and can thus be used for the correction. However, the correct length value cannot be precisely calculated from the additionally measured pressure p
1
. An approximation function is, however, sufficient as a correction in order to increase the measurement accuracy in the desired manner.
In accordance with one embodiment of the invention, the following function is used to calculate the length instead of the pressure p
2
:
f(p
1, p
2
)=αp
1
-p
2
wherein a is a basically freely selectable constant. Through this approximation function a considerable improvement of the measured values is already achieved. It is preferred if a suitable value for a is determined by measurement of the pressures p
1
and p
2
at, in each case, two different air supply pressures A and B and by calculation of ac in accordance with the formula:
wherein P
1A
. . . p
1
at air supply pressure A
P
1B
. . . p
1
at air supply pressure B
P
2A
. . . p
2
at air supply pressure A
P
2B
. . . p
2
at air supply pressure B
is determined. It has been shown that a high measurement accuracy is achievable herewith.
In accordance with a further embodiment of the invention, the two air supply pressures A and B are selected symmetrically to a desired air pressure, in particular to the intended air pressure, for the determination of α. The measurement is executed at the desired spacing X
0
of the measurement nozzle. Through this design it is ensured that the highest measurement accuracy is achieved in the region of the intended air pressure and of the desired distance. A change of the desired distance can in this respect be realized simply by a change of α. That is to say, a zero point shift for the adjustment of the measuring apparatus can be effected via a: without mechanical changes at the measurement apparatus.
In accordance with a further embodiment of the invention, the correction is carried out in accordance with the formula:
x
e
=x−γΔpΔx
where Δp=p
1
−P
0
and Δx=x−X
0
, with
x
e
. . . being the corrected length value
x . . . being the length value calculated from p
2
γ . . . being a basically arbitrarily selectable constant.
Through this embodiment an even higher accuracy of the measured values can be achieved, in particular with larger air pressure fluctuations and larger deviations of the distance from the desired value.
In accordance with a further embodiment of the invention, a suitable value for γ is determined by determination of the indicated length value at, in each case, two different distances set by differing measurement and at two different pressures p
1
and by calculation of γ in accordance with the formula:
with
X
2Ph
. . . being the indicated length value at a distance X
2
and pressure P
h
X
2P0
. . . being the indicated length value at a distance X
2
and pressure P
0
X
0
P
0
. . . being the indicated length value at a distance X
0
and pressure P
0
.
It has been shown that particularly precise measurement values can be achieved with this γ value.
In accordance with a further design of the invention, X
0
is selected as the desired spacing, and P
0
as the desired pressure. Through this embodiment the highest accuracy is advantageously achieved at the desired spacing and at the desired pressure.
In accordance with a further embodiment of the invention, p
1
and p
2
can be calibrated values of the pressure sensors before the pre-nozzle and between the pre-nozzle and the measurement nozzle. In accordance with another design, p
1
and p
2
can be uncalibrated values of these pressure sensors, and in particular their voltage values v and u.
When using uncalibrated pressure values the calculation of the corrected length value takes place in accordance with a further embodiment of the invention, preferably in accordance with the formula:
wherein
f
e
(v, u)=f(v, u)−γΔvΔf(v, u)
Δv=v−V
0
Δf(v, u)=vα−u−F
0
F
1
=f
e
(v, u) at a distance x=X
1
and
F
2
=f
e
(v, u) at a distance x=X
2
with
x
e
. . . being the corrected length value
v . . . being the uncalibrated value of the pressure sensor before the pre-nozzle
u . . . being the uncalibrated value of the pressure sensor between the pre-nozzle and the measurement nozzle.
With this formula a very precise measurement value can be calculated without the calibrated pressure sensors being necessary. Accordingly, cost saving, uncalibrated pressure sensors can be used for the measurement apparatus.
BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1
illustrates a measurement arrangement for carrying out the method of the invention,
FIG. 2
illustrates a customary network without a bridge,
FIG. 3
a
illustrates the layout of a customary pneumatic length measurement system, which has to be fed with a stable air supply pressure,
FIG. 3
b
illustrates the named characteristic p
2
/p
1
,
FIG. 4
illustrates the pressure dividing circuit to
FIG. 3
,
FIG. 5
illustrates the dependence of the measured pressure p2 on the gap x, wherein the quasi-linear measurement range around the point of inflection X0 is utilized,
FIG. 6
illustrates the classical pressure measurement method with a pneumatic measurement bridge,
FIG. 7
illustrates an equivalent circuit to show that the network bridge circuit is independent of the air supply pressure p1,
FIG. 8
illustrates the characteristics of the measurement branches of
FIG. 2
, which are strongly influenced by the air supply pressure,
FIG. 9
illustrates the compensation of the measurement valve via an additional pressure sensor,
FIG. 10
illustrates the straight pressure lines p
1
and p
2
in relation to the reference pressure abscissa,
FIG. 11
illustrates the straight pressure lines p
1
and p
2
in relation to the air supply pressure p
1
,
FIG. 12
illustrates a representation of the data of Table 4, in accordance with which the measured distance x=X
0
is independent of the air supply pressure p
1
,
FIG. 13
illustrates the indication error, which increases with the increase in distance to the distance x=X
0
,
FIG. 14
illustrates the indicated distance in dependence on the air pressure,
FIG. 15
illustrates the deviations from the desired value, with an increase of the stability of the measured value of up to ten times being obtained by use of the equation 22,
FIG. 16
illustrates the transformation f(v, u), the crude value of which is independent of the air supply pressure only around the point x=X
0
,
FIG. 17
illustrates the differences f(v, u)−f(V
0
, u) which lead to the deviations of the air supply pressure p
1
to the reference pressure p
1
=P
0
and which are proportional to the distance from the reference point x=X
0
,
FIG. 18
illustrates the transformation for example(v, u), and
FIG. 19
illustrates the differences for example(v, u)−f(V
0
, u), in accordance with which the sensitivity achieved against air pressure fluctuations of p
1
is reduced by the factor 10.
DESCRIPTION OF SPECIFIC EXEMPLARY EMBODIMENTS
The measurement arrangement shown in
FIG. 1
includes a pressure source
1
, a filter
2
, a pressure regulator
3
, a pre-nozzle
4
a measurement nozzle
5
, a first pressure sensor
6
arranged before the pre-nozzle
4
and a second pressure sensor
7
arranged between the pre-nozzle
4
and the measurement nozzle
5
. An impact plate
8
representative for a workpiece to be measured is shown in front of the air outlet opening of the measurement nozzle
5
.
In the following the relationships in customary pressure measurement methods and in the pressure measurement method of the invention will be explained once again in detail, with special embodiments being used for the explanation of the invention.
1. Introduction
In the 70-year development history of the air pressure measurement method, a plurality of variants arose, of which the most important are nowadays summarized in DIN 2271[1]. The important advantages of the pneumatic length measurement technique are: simple mechanical design, no movable parts (measurement probe/measurement lever), self-cleaning of the measurement location, contact-free sensing, small dimensions of the measurement means (measurement nozzles), high accuracy with the resolution above 0.1 μm [3, 4, 8].
As a result of the advantages listed, a widespread use of the pneumatic measurement principle took place in past years. The high measurement accuracy, rapid detection of the measured value and insensitivity of the measurement method with respect to external influences (chips, emulsions, contamination) are greatly appreciated, particularly in manufacturing [2, 3, 5, 7].
A general disadvantage of modern pneumatic measurement lies in the fact that the supply pressure of the measurement apparatus must be constant. Embodiments without a pneumatic bridge (for example throughflow measurement methods, speed measuring methods and pressure pressuring methods [1, 7, 9]) are particularly affected. In order to guarantee a measurement accuracy in the μm range, air pressure fluctuations are permitted in the small mbar range between the supply unit and the measurement apparatus. In order to be able to guarantee constant air pressure conditions within pneumatic measurement instruments, a compensation bridge has to be added. This in turn caused higher manufacturing costs and increased the compressed air consumption of the measurement instrument.
In this manuscript a new embodiment of the pneumatic high pressure measurement method of DIN 2271 is described. The pneumatic length measurement method portrayed here achieves a measured value which is stable against the fluctuations of the air supply pressure, even without a pneumatic bridge. The pneumatic bridge is replaced by an additional air pressure sensor, with which the actual air supply pressure is measured. The actual measured value is thereby continually corrected. The derivation of the compensation calculation will be explained with reference to a computational example.
2. Theoretical Principles of the Pressure Measurement Method
Pneumatic measurement instruments are divided into two groups. This classification takes account of the supply pressure of the individual measurement instrument. For low pressure measurement instruments, the supply pressure amounts to ≦0.1 bar, and for high pressure measurement instruments the supply pressure lies at ≧0.5 bar. At this point the characteristics of the pneumatic networks for high pressure measurement methods will be described.
2.1 Network Without Bridge
The basic layout of a pneumatic length measurement system, which operates in accordance with the pressurized air measurement method, is shown in FIG.
2
. The supplied air is cleaned by a filter and set constant with the aid of an air pressure regulator. At a stable input pressure p
1
the area ratio A
2
/A
1
(measurement nozzle/measurement pre-nozzle) essentially influences the airflow and thereby the measurement pressure p
2
, which is indicated by a manometer.
A precise calculation of the measurement pressure p
2
in dependence on the gap x is made more difficult by the compressibility of the air, friction losses and heat exchange [2, 6, 9]. To a first approximation it can be assumed that the flow takes place with an adiabatic process, i.e. the heat exchange can be ignored (T
1
=T
2
=constant; see FIG.
3
). The areas or cross-sections are calculated in accordance with the following formulae:
A
1
=D
1
2
π4; A
2
=D
2
x (1)
wherein D
1
and D
2
are the smallest diameter of the pre-nozzle and the measurement nozzle. The theoretical measurement range is limited by the diameter of the measurement nozzle, i.e. x
max
=D
2
/4.
The geometry of the pre-nozzle and measurement nozzle, together with the selected feed pressure p
1
and the measurement pressure p
2
, play a decisive role for the numerical size of the outflow and the energy losses. Accordingly, the nozzles are frequently rounded at the inlet [7, 9], i.e. the flow passage of the pre-measurement nozzle has almost the shape of a Laval nozzle [9]. The normalized characteristic with respect to the supply pressure is shown in
FIG. 3
b
. If the measurement nozzle is pressed closed by an impact plate (x=0), the measurement pressure P2 is the same as p
1
.
An enlargement of the gap beyond a quarter of the measurement nozzle diameter results in practically no change of the measurement pressure p
2
(normally the measurement nozzle ends with the narrowest cross-section).
On the assumption that the cross-sections of the supply lines for the pre-nozzles and measurement nozzles are adequately large and that T
1
=T
2
is the same, the same quantity of air flows through the pre-nozzle and measurement nozzle in the stationary state. Thus, the following equation can be drawn up [9]:
wherein: κ=1.4; A
2max
=D
2
2
π/4; α
1
and α
2
are outflow figures from the pre-nozzle and measurement nozzle; A
1min
is the narrowest measurement nozzle cross-section (approximated A
1min
=A
1
); p
1min
and p
2min
are the pressures in the narrowest cross-sections of the pre-nozzle and measurement nozzle.
For a feed pressure p
1
between 2.5 and 3.2 bar the equation (2) can be simplified:
In
FIG. 4
a simple, analog model of a length measurement pressure converter (see
FIGS. 1 and 2
) is shown. The voltage is replaced with the air pressure and the resistance with the corresponding pneumatic resistance R
p
=k/A
2
, with k being a constant. From this the following equation follows:
wherein
The factor β is also dependent on the outflow figures α
1
and α
2
[2, 6]. In
FIG. 4
plots of the equation (4) are shown simultaneously with a recorded characteristic (D
1
=0.86 mm; D
2
=1 mm; β=0.5; p
1
=3.0 bar).
2.2 Network with Bridge
In the pressure difference measurement method, the difference δ of the measurement pressure p
2
relative to a comparison pressure p
3
is used for the formation of the measured value (see FIG.
6
). Since the comparison pressure p
3
is directly related to the feed pressure p
1
, the influence of the pressure fluctuations in the air supply is minimized.
If the reference pre-nozzle corresponds to the pre-measurement nozzle, and if the reference nozzle was so selected that its cross-section corresponds to the gap x at the point of inflection X
0
(see FIGS.
5
and
6
), then the measurement value around the “zero-point” becomes independent from the air supply pressure p
1
(for example in the range 3.0 bar±0.3 bar). This is illustrated with respect to the equivalent circuit in FIG.
7
. The pressure difference δ indicated is calculated in accordance with the equation (5):
The constant K
3
is determined by the cross-sections of the reference pre-nozzle and the reference nozzle. If the pressures p
2
and p
3
are the same, i.e. K
3
=1/(1+Bx
2
), then the measured value is independent of p
1
. The greater the difference δ, the more the stability of the measurement value depends on the air supply pressure (see Table 4).
3. The Calibration in the Pressure Measuring Method
In
FIG. 8
the characteristics found for a pneumatic length measuring system (network without bridge; see
FIGS. 2 and 3
) are shown for different air supply pressures (p
1
=2.9, p
1
=3.0, and p
1
=3.1 bar). The center of the linear measurement range (i.e. the point of inflection) wanders in this case around X
0
=160 μm. Assuming that the pressure p
1
can be kept at a constant value (for example p
1
=3.0 bar) and that the linear measurement range was selected for the measurement (for example in this case ±20 μm about the point of inflection X
0
) then the equation (6) can be used for the calibration. X
1
and X
2
are two specific distances, which can be measured with a micrometer screw (e.g. X
1
=X
0
−20 μm=140 μm and X
2
=X
0
+20 μm =180 pm). The value P
x1
is the value of the measurement pressure p
2
at a distance X
1
, and P
x2
is the value of p
2
at a distance of X
2
.
The bridge circuit can also be calibrated with the formula (6). The variable p
2
is now correspondingly replaced by the variable δ [8].
In practice, air pressure sensors are used, which deliver at the output a non-calibrated output signal with the linear dependency on the measured pressure (i.e. u=kp
2
+b (current or voltage)). The equation (6) also applies in this case; however, the unknowns p
2
, P
x1
and P
x2
are substituted by the set u, U
x1
and U
x2
.
4. Improved Embodiment of the Customary Air Pressure Measurement Method
It is evident from
FIG. 8
that the pressure fluctuations Δp
1
above ±100 mbar cause a falsification of the spacing Δx=25 μm. That is to say, for a pressure fluctuation of ca. ±3.3% of the air supply, this corresponds to 30% to 50% of the measurement range. A shift of the point of inflection X
0
by 10 to 15 μm can be clearly recognized. The calibration formula (6) is only usable for the case of a stable and constant air pressure supply p
1
(as a rule a stability of the supply pressure of ≦1% is required).
At this point a transformation of the crude measured value is achieved with the aid of an additional pressure measuring sensor (that is, the sensor with the output parameter v; see FIG.
9
). Thus, the calibration formula (6) can continue to be used. As a consequence, the stability of a bridge circuit is achieved. For the pressure change of the air supply of above ±100 mbar, the measured value may only be permitted to change by ±0.5 μm within the entire measurement range.
The direct path of solution of the object set out above is to form the measurement pressure p
2
with a transformation of p
2
and the actual value of the air pressure supply p
1
.
And thus the result can be inserted into the equation (6), i.e. p
2
→f(p
1
, p
2
). The compensation condition which must therefore satisfy the transformation is:
f(p
1
, p
2
)=const. at x=X
0
(7)
Thus a measured value can be formed in the middle of the measurement range independently of the air supply pressure.
TABLE 1
|
|
Changes of p
1
[bar] and p
2
[bar] at the distances X
1
, X
0
and X
2
|
(reference to FIG. 4 and FIG. 8)
|
p
1
p
2
p
1
− p
2
p
1
p
2
p
1
− p
2
p
1
p
2
p
1
− p
2
|
|
X
1
2.901
2.245
0.656
3.001
2.329
0.672
3.102
2.413
0.688
|
X
0
2.899
2.102
0.797
3.000
2.182
0.818
3.101
2.263
0.838
|
X
2
2.902
1.966
0.935
3.000
2.041
0.959
3.101
2.118
0.984
|
|
In Table 1 the results found in
FIG. 8
are shown. The intuitive nature of the compensation of the measurement pressure p
2
with the difference to p
1
, i.e. the transformation:
f(p
1
, p
2
)=p
1
−p
2
(8)
does not appear to be plausible, because it does not behave constant in accordance with Table 1 at a fixed spacing x (e.g. at x=X
0
the pressure difference p
1
−p
2
for Δp
1
=202 mbar varies by Δ(p
1
−p
2
)=41 mbar).
It can be derived from Table 1 that the pressure ratio p
1
to p
2
is not the same with reference to a rising reference pressure p. This is shown in the sketch in
FIG. 10. A
change of p with reference to P
A
to P
B
results at p
1
as a change from P
1A
to P
1B
, i.e. at P2 as a change from P
2A
to P
2B
.
The condition (7) and the transformation (8) can now be newly written in this system:
p
1
(p)−p
2
(p)=const. at any desired p (9)
Since one can assume linear relationships, the pressure conditions of
FIG. 10
can be further developed:
p
1
(p)−p
2
(p)=k
1p
+m
1
−k
2
p−m
2
=p(k
1
−k
2
)+m
1
−m
2
=const (10)
with k
1
and k
2
describing the straight line coefficients and m
1
and m
2
the points of intersection on the Y axis at X p=0 for the straight lines p
1
(p) and p
2
(p). Since m
1
and m
2
are constant, the equation (10) is satisfied if k
1
=k
2
. Since, however, k
1
is a fixed constant, a correction value is necessary in order to be able to satisfy the required condition, i.e.:
k
1
α=k
2
(11)
The coefficient cc can simply be calculated:
Since, in accordance with (9), no precise value but only a constant value is demanded, the coefficients of (10) m
1
and m
2
are omitted, i.e. the transformation (11) satisfies the condition (9):
p
1
(p)−p
2
(p)=k
1
αp−k
2
p (13)
The reference pressure p can be selected as desired. If we set p=p
1
, then the expressions (12) and (13) become more specific, as can be seen in
FIG. 11
, i.e. the required transformation reads:
f(p
1
, p
2
)=αp
1
−p
2
(14)
The calibration formula (6) is now:
with F
1
being the value of the transformation f(p
1
, p
2
) at the first (X
1
) and F
2
being the value of the transformation f(p
1
, p
2
) at the second adjustment master (X
2
). Instead of the precise air pressure values p
1
and p
2
, direct uncalibrated sensor outputs v and u can now be used in the equations (14) and (15): f(v, u)=αv−u.
5. Example and Results
The calibration of the P/E transducer in accordance with the newly described embodiment takes place in two steps. In the first step the compensation coefficient α is determined.
TABLE 2
|
|
Determination of α
|
p
1
[bar]
v [V]
u [V]
|
|
2.8
6.988
5.061
|
3.2
7.983
5.853
|
|
This step is carried out only once. First, the impact plate is fixed at the “0”-position (X
0
). For the nozzle combination D
1
=0.86 mm and D
2
=1 mm (see FIG.
9
), X
0
=160 μm. Then the air pressure supply is set to ca. 2.8 bar, and the output values of the sensors (v and u) are taken on (see FIG.
9
). The second value pair (v, u) is recorded at an increased pressure of ca. 3.2 bar. The constant a is calculated via the stored values in Table 2:
In the second step, the v and u values are adopted in the fixed positions X
1
=X
0
−20 μm and X
2
=X
0
+20 μm, with the air pressure now being set in the pressure air supply at ca. 3.0 bar. In this way Table 3 can be produced:
TABLE 3
|
|
Crude values at p
1
ca. 3.0 bar
|
v [V]
u [V]
f (v, u)
|
|
X
1
7.487
5.821
0.13850151
|
X
2
7.486
5.103
0.85570553
|
|
From this table F
1
=0.13850151 and F
2
=0.85570553 are calculated. The measurement formula (15) reads, after insertion of the calculated constants:
x(v, u)=v44.3934991−u55.7721358+132.28 [μm] (16)
In Table 4 the newly described pressure measuring method is compared with the classical variant (embodiment with bridge) (see FIG.
5
). The data found, which were determined via the derived transformation (14) and the measurement formula (15), correspond.
TABLE 4
|
|
Determined and desired measured values for the different air supply pressures
|
A: classic method
B: new method
|
p
1
desired:
desired:
desired:
desired:
desired:
desired:
|
[bar]
X
1
[μm]
X
0
[μm]
X
2
[μm]
X
1
[μm]
X
0
[μm]
X
2
[μm]
|
|
2.7
141.8
160.5
178.6
141.7
160.4
178.5
|
2.8
141.2
160.4
178.9
141.1
160.2
179.0
|
2.9
140.6
160.3
179.5
140.6
160.3
179.6
|
3
140
160.3
180
140
160.3
180
|
3.1
139.4
160.3
180.5
139.4
160.3
180.5
|
3.2
138.7
160.2
181
138.7
160.2
180.9
|
3.3
138.1
160.1
181.5
138.2
160.2
181.4
|
|
A: classic bridge circuit
|
B: proposed embodiment of the pressure measuring method without bridge shows the same results
|
6. Summary
The described new embodiment of the pressure measuring method has a synergy between the stability of a bridge circuit and the economy of air consumption of a classical pressure measurement method in accordance with DIN 2271. The main features are: lower mechanical effort in manufacturing the transducer, halved air consumption during the measurement operation, stable measurement signal, zero point shift without a mechanical change (i.e. range of insensitivity to the air supply pressure is determined by α).
7. Literature
[1] DIN 2271 part 1, part 2, part 3 and part 4
[2] Eler, R.: “Über den Stand der pneumatischen Feinmeβtechnik und zu Fragen des erforderlichen Gerätesortiments”, (Concerning the state of pneumatic fine measurement technology and questions of the required assortment of instruments), Feingerätetechnik 12/1961, pages 563 to 571
[3] Falkenberg, Y., Werner F.: “Hochleistungsschleifscheiben prozeβnah messen”, (Process orientated measurement of high performance grinding discs), Feinwerktechnik & Meβtechnik 102 (1994), 5-6, pages 256 to 262
[4] Gale, G. E.: “A thickness measuring device using pneumatic gauging to detect the sample”, Meas. Sci. Technol. 6/1995, pages 1566 to 1571
[5] Lehnert, E.: “Pneumatische In-Process-Meβtechnik”, Qualitätstechnik, (Pneumatic in-process measurement technology, Quality Technology), QZ 12/1987, pages 583 to 586
[6] Lotze, W.: “Neue Methoden zur Berechnung pneumatischer Feinzeiger”, (New method for calculating pneumatic fine indication apparatus), Feingerätetechnik 6/1966, pages 275 to 281
[7] Lotze/Kondaschewski: “Meβsteuergerate spanender Werkzeugmaschinen”, (Measurement control instruments of chip forming machine tools), Berlin, VEB Verlag Technik (1974)
[8] Stamenkovic, M.: “Genaues Kalibrieren und Messen ohne genaue Positionierung des Werkstüickes”, (Precise calibration and measurement without accurate positioning of the workpiece), Feinwerktechnik & Meβtechnik 9/1996, pages 651 to 654
[9] Wiemer, A.: “Pneumatische Längenmessung”, (Pneumatic length measurement), VEB Verlag 1970
Section: Compensation of the Supply Pressure Over the Whole Measurement Range
Variant 1: Compensation Related to the Pressures p
1
and p
2
The numerical data of Table 4 are shown graphically in FIG.
7
. The measurement values (which were found in accordance with the formula 14) around the point X
0
=160 mm are practically independent of the air supply pressure (Y axis in FIG.
12
). For the larger deviations of X
0
(e.g. for more than ±2-3 μm) the indication errors for the extreme air pressure fluctuations (±0.3 bar about P
0
=3.0 bar) are no longer negligible. In the illustrated example the read of the indicated measured values for X
1
and X
2
lies at ca. 3 μm (see Table 4 and FIGS.
12
and
13
).
In accordance with these data a general behavior of the pneumatic bridge is presented in FIG.
13
. For the pressure p
1
>P
0
=3.0 bar and for x>X
0
(e.g. x=X
2
) a falsified enlargement of the indication would result (see FIG.
13
). For the same distance (x>X
0
) a reduction of p
1
(p
1
<P
0
) causes a production of the indication of the measured distance. If, however, x is smaller than X
0
(e.g. if x=X
1
), the sign of the falsification will invert.
In general form, the value y, which is read off, can be considered as the sum of the true value x
e
at air pressure p
1
=P
0
and of the generated error ε, i.e.:
x=x
e
+ε (17)
As can be seen from
FIGS. 12 and 13
, it can be assumed that the error ε directly follows the deviation Δx=x−X
0
and the pressure difference Δp=p
1
P
0
, i.e.:
ε=γΔp Δx=γ(p
1
−P
0
)(x−X
0
) (18)
where γ is a constant of the proportionality. It is assumed that at the pressure p
1
=P
0
the error ε is equal to zero (ε=0), i.e. that in this case x=x
e
over the full measurement range. Thus, the equation for the corrected measured value results:
xe=x−γΔpΔx (19)
The constant γ is derived in relation to a specific pressure P
0
and to a specific mean value X
0
:
wherein X
2Ph
is the indicated measured value at the pressure p
1
=P
h
and x=X
2
; X
2Po
is the indicated measured value at the pressure p
1
=P
0
and x=X
2
; X
0Po
is the indicated measured value at the point x=X
0
measured at the pressure p
1
=P
0
. For the calculation of the equations (19) and (20) already 10 calibrated values are necessary, or the distance value y=x must first be formed in accordance with equation (15) and then Δx calculated. In the same way, the calibrated pressure p
1
is required for the calculation of Δp.
The developed expression (19) now reads:
x
e
=x(1+γP
0
−γp
1
)+X
0
γp
1
−P
0
X
0
γ (21)
with the value x being found in accordance with the equation (15). The equation (21) can be used for the compensation of the classic network circuit with the pneumatic bridge. In the summarized form of (15) and (21) x
e
is calculated in the following manner:
As an example, the values are taken over from the Table 4 (these data are shown graphically in FIG.
12
), with the constant γ being:
Since the already calibrated x values are available in the Table 4, the following expression is achieved after inserting the values of the Table 5 into the equation (15):
x
e
=x(1+γP
0
−γp
1
)+X
0
γp
1
−P
0
X
0
γ=x(1.82423−0.2747p
1
)+44.038p
1
−132.1 (23)
TABLE 5
|
|
Inserted values for the compensation of p
1
|
|
|
X
1
= 140 μm
X
2P0
= 180 μm
P
0
= 3 bar
|
X
0
= 160 μm
X
2Ph
= 178,5 μm
P
h
= 2, 7 bar
|
X
2
= 180 μm
X
0P0
= 160,3 μm
|
|
The value x is fed into (23) in μm and the pressure p
1
is fed into (23) in bar. The numerical data from the Table 4 (column B) were treated in accordance with (21) and shown as results in Table 6.
In the example the nozzle combination D
1
=0,86 mm and D
2
=1 mm was taken (see FIG.
9
). With this nozzle set a measurement range ±50 μm is measured in practice. The transducer was calibrated at the points X
1
=140 μm and X
2
=180 μm. Because of the linearity error, X
0
=160 μm is indicated falsified by 0.3 μm. This error does not influence the stability calculation. Accordingly, in Table 6 the compensation relates to the desired value X
0
=160,3 μm.
TABLE 6
|
|
Found values and desired measured values for the different air supply pressures
|
A: difference from desired values
B: measured values x
e
|
p
1
desired:
desired:
desired:
desired:
desired:
desired:
|
[bar]
0 [μm]
0 [μm]
0 [μm]
X
1
[μm]
X
0
1
[μm]
X
2
[μm]
|
|
2.7
0.17
0.11
0.00
140.17
160.41
180.00
|
2.8
0.04
−0.11
0.03
140.04
160.19
180.03
|
2.9
0.06
0.0
0.13
140.16
160.30
180.13
|
3
−0.11
0.0
0.00
140.00
160.30
180.00
|
3.1
0.03
0.0
−0.05
139.97
160.30
179.95
|
3.2
−0.11
−0.11
−0.23
139.89
160.21
179.77
|
3.3
0.02
−0.11
−0.34
140.02
160.21
179.66
|
|
A: difference from the desired values x (at p
1
= P
0
)
|
B: values found in accordance with the equation (23)
|
The calculated measured values xe have almost no changes for the customary air pressure fluctuations (above ±0,1 bar), i.e. the changes of x
e
that are found lie under ±0,1 μm (see Table 6).
The ±0,1 bar-band is also almost preserved for air pressure fluctuations of 2,8 bar to 3,2 bar, which in practice are exaggerated.
In comparison to a classic pneumatic pressure difference measurement method with a bridge circuit (or in comparison to the improved pressure measuring methods described in Sections 4 and 5), the compensation equation (22) shows an increase of the stability of the measured value by up to ten times. If the value x differs greatly from the X
0
, the relative gain in stability is even greater.
Variant 2: Compensation Related to the Pressure Sensor Values v and u:
For the calculation of x
e
in accordance with the formula (22) two calibrated pressure sensors are required for the precise measurement of p
1
and p
2
. If a classic PEW-transducer with a bridge circuit is to be compensated in this way, only one sensor is required in accordance with the equation (21). In this case the air consumption, however, doubles. Since these requirements make the use of the described compensation more expensive, a form of the equation (22) independent of p
1
and p
2
is desired. In this chapter the derivation of such a compensation formula will be explained.
As illustrated in
FIG. 9
, the pressure sensors deliver uncalibrated voltages v and u. The voltage v has a linear dependency on the pressure p
1
and the voltage u has a linear dependency on the pressure p
2
, i.e. the following can be written:
v=mp
1
+n; u=kp
2
+b (24)
where m, n k and b are constants. These constants are determined by the electronic amplifiers and characteristics of the pressure sensors. As was described in chapters 4 and 5, the transformation (14) can be calculated directly with these voltage values, i.e.:
f(v, u)=αv−u (25)
The coefficient a is calculated at the distance x=X
0
as was described in chapter 5:
wherein the indices h and 1 designate the overpressure and the underpressure. For example, v
h
and u
h
are output voltages at the supply pressure p
1
=P
h
; in the same way v
1
and u
1
are output voltages at the supply pressure p
1
=P
1
. The pressures P
h
And P
1
do not need to be accurately determined. They should, however, be located around the target pressure p
1
=P
0
(e.g. P
h
=2,8 bar and P
1
=3,2 bar or vice versa).
The expression (25) has an independence of v (i.e. independence from p
1
) only at a fixed point x=X
0
; i.e. in association with the formula (14) the measured value behaves as in a classic pneumatic bridge. This is illustrated in FIG.
16
. The values of f(v, u) at a constant value v=V
0
, i.e. the values of f(V
0
, u), are shown on the horizontal axis. This state corresponds to p
1
=P
0
. It can be seen from
FIG. 16
that all parametricized straight lines f(v, u) intersect at one point. This point is f(V
0
, U
0
), i.e. the value of f(v, u) determined at the distance x=X
0
and at the supply pressure p
1
=P
0
. In
FIG. 17
the difference between f(f, u) and f(V
0
, u) is shown on the vertical axis.
The differences f(v, u)−f(V
0
, u) are proportional to the deviations of the supply pressure p
1
relative to the reference pressure p
1
=P
0
and to the distance to the reference point x=X
0
.
The expression (25) with which f(v, u) is defined should now be so converted that the result of f(v, u) remains constant at a desired point x=X
i
, even with different values of p
1
. For this purpose a correction of f(v, u) is required. Since v and u are linear functions of p
1
and p
2
these corrections will have the form of equations (17, 18) and (19):
f(v, u)=f
e
(v, u)+ω (27)
ω=γΔvΔf(v, u)=γ(v−V
0
)(vα−u−F
0
)=γ(v−V
0
)(vα−u−F
0
)=γ(v−V
0
)(vα−u−V
0
α+U
0
) (28)
where F
0
=V
0
α−U
0
is the value of the transformation of f(v, u) at the point x=X
0
at the reference pressure p
1
=P
0
. The reference pressure P
0
corresponds to the output voltage V
0
. Under the same conditions the value U
0
is the voltage magnitude from the pressure measurement sensor.
The corrected transformation (25) now reads:
f
e
(v, u)=f(v, u)−γΔvΔf(v, u)=vα−u−γ(v−V
0
)(vα−u−V
0
α+U
0
) (29)
The constant γ is calculated in accordance with the following equation:
wherein: F
2Vh
is the value of f(v, u) at v=V
h
and x=X
2
; F
2Vo
is the value of f(v, u) at v=V
0
and x=X
2
; F
0Vo
is the value of f(v, u) at v=V
0
=x=X
0
.
In
FIGS. 18 and 19
results of (28) are entered on the vertical axis (the numerical data are shown in Table 9). The horizontal axis is the same as in
FIGS. 16 and 17
. A brief comparison of
FIGS. 16
with
18
and
17
with
19
documents a significant (ca. 10 times) reduction of the differences between the actual values f
e
(v, u) and the reference values f(V
0
, u).
Since the expression (29) provides a crude value independent of v (i.e. of p
1
), it can be directly substituted in the calibration formula (15):
with F
1
=f
e
(v, u) at the distance x=X
1
and F
2
=f
e
(v, u) at the distance x=X
2
.
This results in the calibration equation (32), which does not require any calibrated air pressure values as argument, but rather only uses the pressure sensor output values v and u, i.e.:
The required parameters (α, γ, V
0
, etc.) can be calculated from 6 or 5 calibration steps. If one sets X
1
=X
0
for the original calibration, then the step 6 is superfluous. This is in particular necessary if an analogue circuit solution of the equation (xx) is to occur. In order to be able to show a more accurate comparison between the compensation method and improved pressure measurement methods, X
1
is considered separately from X
0
in the described example, i.e. the step 6 is used. The required steps are shown in the Table 7 (in this Table the data of Table 2 are used). In this way the same α value is used and a direct comparison of the results is made easier.
TABLE 7
|
|
Carrying out the original calibration in 6 or 5 steps
|
Sensor value
Sensor value
Supply pressure
Distance
|
Step
v [V]
u [V]
[bar]
x [μm]
Remarks
|
|
1
V
1
= 6.988
U
1
= 5.061
2.8
x = X
0
calculation of α
|
2
V
2
= 7.983
U
2
= 5.853
3.2
x = X
0
calculation of α
|
3
V
3
= 7.484
U
3
= 5.454
3.0 (exact!)
x = X
0
|
4
V
4
= 7.486
U
4
= 5.103
3.0 (exact!)
x = X
2
|
5
V
5
= 6.738
U
5
= 4.534
2.7
x = X
2
|
6
V
6
= 7.487
U
6
= 5.821
3
x = X
1
Omitted if X
1
= X
0
|
|
Step 1: The impact plate is brought to a distance x = X
0
= 160 μm from the nozzle D
2
(see FIG. 9). The air supply pressure is set to ca. 2.8 bar. This values does not need to be accurately set. It is, however, important that the value is 0.1 to 0.3 bar smaller or greater than “zero”-pressure P
0
= 3.0 bar. The actual sensor voltage values v = V
1
and u = U
1
are stored.
|
Step 2: With the distance x = X
0
unchanged, the supply pressure is set in the opposition direction from the pressure set in step 1. In this case the value of ca. 3.2 was selected. The values V
2
and U
2
correspond to this state.
|
Step 3: With the distance x = X
0
unchanged, the supply pressure is accurately set to 3.0 bar. In this step the reference state is reached (x = X
0
; pressure = P
0
), i.e. the values V
0
and U
0
are recorded.
|
Step 4: With an unchanged “zero”-pressure P
0
= 3.0 bar, the impact plate is brought to a distance x = X
2
= 180 μm. The supply pressure should be precisely set up as long as the value V
4
corresponds with the value V
3
= V
0
.
|
Step 5: With the distance x = X
2
unchanged, the supply pressure is set to ca. 2.7 bar. Instead of 2.7 bar a larger pressure than P
0
= 3.0 bar can also be selected, for example 3.2 bar.
|
Step 6: With the distance x = X
1
unchanged, the supply pressure is set to ca. 3 bar (this step is superfluous if X
1
= X
0
).
|
If the calibration steps have been carried out, then the required parameters are calculated in accordance with the expressions and sequences which are shown in the Table 8:
TABLE 8
|
|
No.
Expression
|
|
|
1
|
|
2
V
h
= V
5
= 6.738
|
3
V
o
= V
3
= V
4
= 6.738
|
4
U
o
= U
3
= 5.454
|
5
F
2V
h
= V
5
α − U
5
= 0.92931257
|
6
F
2V
0
= V
4
α − U
4
= 0.85570553
2
|
7
F
0V
0
= V
3
α − U
3
= 0.50311357
|
|
8
|
|
9
F
1
− V
6
α − U
6
− γ(V
6
− V
3
)(V
6
α − U
6
− F
0V
0
) = 0.13862015
|
10
F
2
− V
4
α − U
4
− γ(V
4
− V
3
)(V
4
α − U
4
− F
0V
0
) = 0.85562905
|
|
11
|
|
The parameters of the equation (20) are found from the data of Table 8
|
In reality, the difference between V3 and V4 will already result. This is also the case in the illustrated example, where V3=7,484 [V] and V4=7.486 [mV]. This difference of ca. 2 mV can be ignored because the sensitivity of the sensors used lies at ca. 2.5 V/bar.
After the insertion of the calculated constants in the formula (12) a specific calibration formula can be written:
x
3
=140+55,78731 (f
3
(v, u)−0,13862) (33)
wherein the function f
3
(v, u) is:
f
3
(v, u)=v0,79598−u−0,108459(v−7,484)(v0,79598−u−0,50311) (34)
The uncompensated transformation f(v, u) is simpler:
f(v, u)=v0,79598−u (35)
TABLE 9
|
|
Measured pressure sensor values for the distances X
1
, X
0
and X
2
with air supply pressure
|
fluctuations from 2.7 to 3 bar.
|
x = X
1
= 140 μm
x = X
0
= 160 μm
x = X
2
= 180 μm
|
p
1
v
u
f (v, u)
f
e
(v, u)
v
u
f (v, u)
f
e
(v, u)
v
u
f (v, u)
f
e
(v, u)
|
|
2.7
6.739
5.196
0.16811
0.14104
6.739
4.86
0.50411
0.50419
6.38
4.534
0.82931
0.85571
|
2.8
6.987
5.404
0.15751
0.13888
6.988
5.061
0.50131
0.50121
6.985
4.723
0.83692
0.85499
|
2.9
7.239
5.613
0.14910
0.13969
7.234
5.255
0.50312
0.50312
7.24
4.915
0.84789
0.85702
|
3
7.487
5.821
0.13850
0.13862
7.484
5.454
0.50311
0.50311
7.486
5.103
0.85571
0.85563
|
3.1
7.739
6.033
0.12709
0.13749
7.738
5.657
0.50229
0.50232
7.738
5.294
0.86529
0.85531
|
3.2
7.985
6.241
0.11490
0.13599
7.983
5.853
0.50131
0.50141
7.986
5.485
0.87170
0.85163
|
3.3
8.234
6.448
0.10610
0.13839
8.237
6.055
0.50149
0.50162
8.238
5.677
0.88028
0.84944
|
|
In Table 9 the crude sensor values v and u are shown together with the calculated transformations f(v, u) (in accordance with 35) and f
3
(v, u) (in accordance with 34).
Around the point x=X
0
practically no differences between f(v, u) and f
3
(v, u) can be noted for pressure fluctuations of p
1
=2.7 to p
1
=3.3 bar. The differences are first evident when the distance x deviates from the center point X
0
(e.g. x=X
1
or x=X
2
). In this case the changes of f(v, u) are significant, with f
e
(v, u) still remaining constant (see Table 8).
The relatively complicated equation (30) can be written in the following general form after calculation of the constants:
x
e
=Av
2
+Bvu+Cv+Du+E (36)
with the constants A, B, C, D and E being:
Thus, the values A=−4.8162, B=6,05065, C=83.494167, D=−01.07039 and E=109.484 result, i.e. the calibration equation:
x
e
=−4.8162v
2
+6.05065vu+83.49417v−101.07039u+109.484 (42)
The Table 10 represents the calibrated values which were calculated in accordance with the expression (42).
TABLE 10
|
|
Derived values and desired measurement values for the different air supply pressures
|
A: difference to the desired values
B: measured values
|
p
1
desired:
desired:
desired:
desired:
desired:
desired:
|
[bar]
0 [μm]
0 [μm]
0 [μm]
X
1
[μm]
X
0
[μm]
X
2
[μm]
|
|
2.7
0.13
0.09
0.00
140.13
160.39
180.00
|
2.8
0.01
−0.07
−0.04
140.01
160.23
179.96
|
2.9
0.06
0.03
0.08
140.06
160.33
180.08
|
3
0.00
0.03
0.00
140.00
160.33
180.00
|
3.1
−0.06
−0.01
−0.02
139.94
160.29
179.98
|
3.2
−0.15
−0.06
−0.22
139.85
160.24
179.78
|
3.3
−0.01
−0.05
−0.35
139.99
160.25
179.65
|
|
A: difference from the desired values x
|
B: derived values in accordance with the equation (24)
|
The result i agrees with the results of Table 8.
Summary
Manuscript:
1. New pressure measuring method, which has the stability of the pressure difference measurement method.
Section:
1. Supply pressure compensation of the classic pneumatic bridge (pressure difference measurement method):
x
e
=x(1+γP
0
−γp
1
)+X
0
γp
1
−P
0
X
0
γ (21)
The above formula can generally be used for the compensation of a sensor bridge circuit.
2. Supply pressure compensation of the improved pressure measuring method, which was described in the manuscript:
Variant 1:
Input parameters are calibrated measured pressures p
1
and p
2
:
Variant 2:
Input parameters are uncalibrated outputs (current or voltage) of the pressure sensors:
in general form:
x
e
=Av
2
+Bvu+Cv+Du+E (36)
Claims
- 1. A method for pneumatic length measurement in which the medium used for the measurement is directed through a pre-nozzle and a measuring nozzle onto an impact plate and the distance of the impact plate from the measurement nozzle is determined from a change of a pressure p2 between the pre-nozzle and the measurement nozzle, wherein pressure p1 in front of the pre-nozzle is additionally measured and the length value resulting from the pressure measurement between the pre-nozzle and the measurement nozzle is corrected in dependence on the pressures p1, p2, and wherein as a correction, the following function is used for the calculation of the length:f(p1, p2)=αp1−p2, wherein α: is a determinable constant.
- 2. A method in accordance with claim 1, wherein a preferable value for α is determined by measurement of the pressures p1 and p2 for, in each case, two different supply pressures A and B and calculation of α is accordance with the formula: α=(P2B-P2AP1B-P1A)with P1A being p1 at the air supply pressure AP1B being p1 at the air supply pressure B P2A being p2 at the air supply pressure A P2B being p2 at the air supply pressure B.
- 3. A method in accordance with claim 2, wherein the two air supply pressures A and B for the determination of α are selected such thatA<Ps<B, (Ps−A)≈(B−Ps), and Ps is the supply pressure at which the corrected measurement is to later be carried out, and wherein the measurements to determine cc are carried out with the measuring nozzle a distance X0 from an impact plate, where X0 is the nominal distance between the measurement nozzle and an impact plate under test.
- 4. A method for pneumatic length measurement in which the medium used for the measurement is directed through a pre-nozzle and a measuring nozzle onto an impact plate and the distance of the impact plate from the measurement nozzle is determined from a change of a pressure p2 between the pre-nozzle and the measurement nozzle wherein pressure p1 in front of the pre-nozzle is additionally measured and the length value resulting from the pressure measurement between the pre-nozzle and the measurement nozzle is corrected in dependence on the pressures p1, p2, and wherein as a correction the following function is used for the calculation of the length:xe=x−γΔpΔx where Δp=p1−P0 and Δx=x−X0, withxe . . . being the corrected length value X0 . . . being the distance of the impact plate from the measurement nozzle x . . . being the length value calculated from p2 P0 . . . being the desired value of p1 at which the measurement is later to be carried out, and γ . . . being a determinable constant.
- 5. A method in accordance with claim 4, wherein a preferable value for γ is determined by calculating the expected spacing between the measurement nozzle and the impact plate while varying the actual spacing and pressure p1, and calculating γ in accordance with the formula γ=X2Ph-X2P0(Ph-P0)(X2Ph-X0P0)whereX2Ph . . . is the calculated length value at the known predetermined spacing X2 and the pressure Ph X2P0 . . . is the calculated length value at the known predetermined spacing X2 and pressure P0 X0P0 . . . is the calculated length value at the known predetermined spacing X0 and pressure P0 X0 . . . is the nominal expected distance between the measurement nozzle and the impact plate to be later measured and P0 . . . is the desired value of p1 at which the measurement is later to be carried out.
- 6. A method in accordance with claim 5, characterized in that X0 is selected as the nominal expected distance between the measurement nozzle and the impact plate to be later measured and P0 as the desired value of p1 at which the measurement is later to be carried out.
- 7. A method in accordance with claim 1, wherein the calculation of the corrected length value takes place in accordance with the formula xe=X1+(fe(v,u)-F1)X2-X1F2-F1,whereinxe . . . is the corrected length value X1, X2 . . . are predetermined values of x f(v, u)=αv−u γ . . . is a constant Vu . . . is v at X0 and P0 α=(uh-uivh-vi), wherein uh . . . is u at p1=Ph ui . . . is u at p1−P1 vh . . . is v at p1=Ph vi . . . is v at p1=PI Ph . . . is a predetermined value of p1>P0 PI . . . is a predetermined value of p1<P0 X0 . . . is a predetermined value of x at which the measurement is to take place (zero position) P0 . . . is a value of p1 at which the measurement is to take place (zero position) fe(v, u)−f(v, u)−γΔvΔf(v, u) Δv=v−Vu Δf(v, u)=vα−u−F0 F1=fe(v, u) at the distance x=X1 F2=fe(v, u) at the distance x=X2 xe . . . is a corrected length value v . . . is an uncalibrated value of the pressure sensor before the pre-nozzle u . . . is an uncalibrated value of the pressure sensor between the pre-nozzle and the measurement nozzle and F0 . . . is the function value at x=X0 and p=P0.
Priority Claims (1)
Number |
Date |
Country |
Kind |
197 33 984 |
Aug 1997 |
DE |
|
PCT Information
Filing Document |
Filing Date |
Country |
Kind |
PCT/EP98/04915 |
|
WO |
00 |
Publishing Document |
Publishing Date |
Country |
Kind |
WO99/08069 |
2/18/1999 |
WO |
A |
Foreign Referenced Citations (3)
Number |
Date |
Country |
2921989 |
Apr 1980 |
DE |
4200401 |
Jul 1993 |
DE |
0671599 |
Mar 1995 |
EP |