(19)
(11) EP 0 977 009 B1

(12) EUROPEAN PATENT SPECIFICATION

(45) Mention of the grant of the patent:
17.10.2007 Bulletin 2007/42

(21) Application number: 98917621.9

(22) Date of filing: 22.04.1998
(51) International Patent Classification (IPC): 
G01B 21/20(2006.01)
G06T 17/00(2006.01)
(86) International application number:
PCT/JP1998/001842
(87) International publication number:
WO 1998/049524 (05.11.1998 Gazette 1998/44)

(54)

METHOD OF DISCRIMINATING SHAPE ERRORS OF FREE-FORM CURVED SURFACE

VERFAHREN ZUR DETERMINIERUNG VON FEHLERN IM FORM EINER WILLKÜRLICHEN KRUMMEN OBERFLACHE

PROCEDE PERMETTANT D'ETABLIR UNE DISTINCTION ENTRE DES ERREURS DE FORME D'UNE SURFACE COURBE A FORME LIBRE


(84) Designated Contracting States:
BE CH DE FI FR GB IT LI NL

(30) Priority: 25.04.1997 JP 10913997

(43) Date of publication of application:
02.02.2000 Bulletin 2000/05

(73) Proprietors:
  • RIKEN
    Wako-shi, Saitama 351-0198 (JP)
  • AMADA COMPANY, LIMITED
    Kanagawa-Ken (JP)

(72) Inventors:
  • KASE, Kiwamu-The Ins. of Physical and Chem. Res.
    Saitama 351-0198 (JP)
  • MAKINOUCHI, Akitake
    Wako-shi, Saitama 351-0198 (JP)
  • KONDO, Tetsuya
    Isehara-shi, Kanagawa 259-1196 (JP)
  • MORI, Naomichi
    Isehara-shi, Kanagawa 259-1196 (JP)

(74) Representative: Grünecker, Kinkeldey, Stockmair & Schwanhäusser Anwaltssozietät 
Maximilianstrasse 58
80538 München
80538 München (DE)


(56) References cited: : 
JP-A- 4 276 513
   
  • PANKANTI S ET AL: "Robust feature detection for 3D Object Recognition and Matching" PROCEEDINGS OF THE SPIE, SPIE, BELLINGHAM, VA, US, vol. 2031, 12 July 1993 (1993-07-12), pages 366-377, XP002211019 ISSN: 0277-786X
  • TANAKA H T ET AL: "Curvature-based face surface recognition using spherical correlation-principal directions for curved object recognition" PROCEEDINGS OF THE 13TH INTERNATIONAL CONFERENCE ON PATTERN RECOGNITION IEEE COMPUT. SOC. PRESS LOS ALAMITOS, CA, USA, vol. 3, 1996, pages 638-642 vol.3, XP002391855 ISBN: 0-8186-7282-X
  • HUANG X ET AL: "Localization and comparison of two free-form surfaces" COMPUTER AIDED DESIGN ELSEVIER UK, vol. 28, no. 12, December 1996 (1996-12), pages 1017-1022, XP002391856 ISSN: 0010-4485
   
Note: Within nine months from the publication of the mention of the grant of the European patent, any person may give notice to the European Patent Office of opposition to the European patent granted. Notice of opposition shall be filed in a written reasoned statement. It shall not be deemed to have been filed until the opposition fee has been paid. (Art. 99(1) European Patent Convention).


Description

BACKGROUND OF THE INVENTION


Technical Field



[0001] The present invention relates to a method of determining a shape error of a free-form surface.

Background Art



[0002] With respect to formings formed by press working such as the body of an automobile, discrepancy has been conventionally evaluated by experiences. That is, the free-form surface has been conventionally evaluated mainly by "visual observation". In recent years, however, designing using a computer (CAD) has been spread and a deformation upon processing can be simulated. Accordingly, means for evaluating a free-form surface, that is, for objectively defining discrepancy in a forming and displaying it is desired.

[0003] Fig. 1 shows a forming sample showing discrepancy in a plate forming. Fig. 2 is a diagram showing an example of a conventional method of evaluating a shape error in a free-form surface. The example relates to a result of numerical simulation using simulation software (ITAS-3D) reported in "Simulation of 3-D sheet bending process" (Takizawa et al., 1991, VDI BERICHTE NR. 894), "Some advances in FEM simulation of sheet metal forming processes using shell elements " (Kawka et al., 1995, Simulation of Materials Processing, Shen & Dawson (eds.), Balkema, Rotteerdam, pp. 735 - 740), and the like. The Article "Robust feature detection for 3D Object Recognition and Matching", Proceedings of the SPIE, Geometric Methods in Computer Vision, 1993, by Sharath Pankati et al relates to the matching of two free-form surfaces. In that disclosed matching process, the principle and mean curvatures of both free-form surfaces are computed which are then segmented into eight different types. The actual matching is based on a minimisation of the Hausdorff-distance between the two segmented curvatures.

[0004] In Fig. 2, a white part shows a shape as a reference (for example, the shape of a die) and a mesh part indicates a shape obtained by a forming simulation. The reference shape and the simulation shape are displayed at the same position and only a part positioning on the front side is displayed. Consequently, a shape error between the simulation shape and the reference shape can be roughly determined from the displayed white and mesh parts. The method has, however, the following problems.
  1. (1) It is necessary to determine a reference position and make the reference shape and the simulation shape accurately coincide with each other at the reference position. The result is largely influenced according to the way the reference position is determined.
  2. (2) Since the position of the other part is largely displaced due to a partial bending, it is difficult to find the cause of occurrence of an error.
  3. (3) How much the shapes coincide with each other as a whole cannot be shown by an objective numerical value.


[0005] Fig. 3 shows CMM data (about 40,300 points) measured by using a three coordinate measuring machine "Mitsutoyo Super BHN 506". Fig. 4 shows the top view of Fig. 3 (about 8,000 points). As shown in these figures, the shape of a forming item actually formed by using a die can be displayed as images as shown in Figs. 3 and 4 by measuring the forming item by the three coordinate measuring machine. From the images, shape errors such as projected and recessed parts and a twisted part can be roughly determined from the views. The method, however, also has the above-mentioned problems (2) and (3) in the numerical simulation and has a problem that (4) when the reference shape is not flat but has a complicated curve, the difference from the result of the three coordinate measurement can be hardly determined.

[0006] As mentioned above, methods of experiment, measurement, and display of result for evaluating the shape error have not been systematized yet. There has not been a simple and clear definition as an index of a forming discrepancy and, further, an evaluation method which can be repeatedly performed has not been existed conventionally.

SUMMMARY OF THE INVENTION



[0007] The present invention is made in order to solve the problems. That is, a principle object of the invention is to provide a method of determining a shape error of a free-form surface which can accurately grasp a different part between two three-dimensional shapes of an actual forming shape and a simulation shape by a computer simulation, a reference shape by CAD, or the like. It is another object of the invention to provide a method of determining a shape error of a free-form surface which can be applied without making the reference positions accurately coincide with each other, find the cause of occurrence of an error such as a partial bending, show how much the shapes coincide with each other as a whole by an objective numerical value, and easily determine the error even if the reference shape is complicated.

[0008] The inventors of the present invention have invented "extended Gaussian curvature" as an evaluation model which does not depend on the coordinate system. According to the invention, a local shape error of a free-form surface is classified into three types (mountain, valley, and twist) by comparing an actual curved surface with, for example, a CAD curved surface as a reference. A method of calculating the ratio of the same labels by using the image processing technique has also been invented. The invention is based on the novel ideas.

[0009] According to the invention, there is provided a method of determining a shape error of a free-form surface by obtaining a principle curvature of a target curved surface S' and a principle curvature of a corresponding position of a reference surface S; and displaying each part by classifying it from the difference between the principle curvatures into (a) a case where the two principle curvatures increase, (b) a case where the two principle curvatures decrease, and (c) a case where one of the principle curvatures increases and the other decreases.

[0010] That is, according to a preferred method of the invention, Δ κ1 = κ'1 - κ1, Δκ2 = κ2' - κ2 are obtained from the principle curvature (κ1', κ2') of the target curved face S' and the principle curvature (κ1, κ2) of the reference curved face S. (1) When Δκ1 ≧0 and Δ κ2 ≧0, (a) it is determined that the two curvatures increase. (2) When Δ κ1 ≦ 0 and Δ κ2≦0, (b) it is determined that two curvatures decrease. (3) When Δ κ1 · Δ κ2 < 0, (c) it is determined that one of the curvatures increases and the other decreases. Preferably, (a), (b), and (c) are determined as mountain, valley, and twist, respectively, and are displayed in different symbols or colors on an image. Further, it is preferable that the ratio of the same labels is calculated from the labels (a), (b), and (c) and is used as a coincidence ratio.

[0011] The Gaussian curvature K is a product κ1κ2 of the principle curvatures κ1 and κ2 of three-dimensional surfaces. (1) When K > 0, it is known that the shape is elliptic. (2) When K = 0, it is known that the shape is parabolic. (3) When K < 0, it is known that the shape is hyperbolic.

[0012] The invention relates to an extended Gaussian curvature. That is, according to the method of the invention, the principle curvature of a target curved surface S' including an error and the principle curvature of a corresponding position of a reference curved surface S are obtained and each part is classified into (a), (b), and (c) from the difference between the principle curvatures, namely, the parts can be displayed while being classified into the case where two principle curvatures increase, the case where two principle curvatures decrease, and the case where one of the principle curvatures increases and the other decreases. Thus, the different part between two three-dimensional shapes can be accurately grasped.

[0013] According to the method, the shape error can be determined by obtaining the principle curvatures of corresponding positions. Consequently, the invention can be applied without making the reference positions of two three-dimensional shapes accurately coincide with each other and the cause of occurrence of an error such as a partial bending can be found.

[0014] Further, by calculating the ratio of the same labels from the labels (a), (b), and (c) and using it as a coincidence ratio, how much the shapes coincide with each other as a whole can be grasped by an objective numerical value and an error can be easily determined even if the reference shape is complicated.

[0015] The other objects and advantageous features of the invention will be made clear from the following description with reference to the attached drawings.

BRIEF DESCRIPTION OF THE DRAWINGS



[0016] 

Fig. 1 is a halftone image displayed on a display of a forming sample showing a discrepancy of a plate forming;

Fig. 2 is a diagram showing an example of a conventional method of evaluating an error in a free-form surface using simulation;

Fig. 3 is a diagram of measurement by a three coordinate measuring machine;

Fig. 4 is the top view of Fig. 3;

Fig. 5 is a halftone image displayed on a display showing a perspective view of an impeller;

Fig. 6 is a diagram showing that a twist deformation according to the method of the invention is performed to an originally twisted shape;

Fig. 7 is a diagram showing that a twist deformation according to the method of the invention is performed to an originally twisted shape;

Fig. 8 shows an example of applying the method of the invention on the basis of measurement data;

Fig. 9 shows an example of applying the method of the invention on the basis of the result of numerical simulation;

Fig. 10 shows another example of applying the method of the invention on the basis of measurement data; and

Fig. 11 shows another example of applying the method of the invention on the basis of the result of numerical simulation.


DESCRIPTION OF THE PREFERRED EMBODIMENTS



[0017] The principle of the method of the invention will be first described.

[0018] Free-form surface S = S(u, v) is expressed by parameters u and v. Expression 1 is a relational expression in differential geometry. The relation is disclosed in, for example, "Curves and Surfaces for Computer Aided Geometric Design" (Farin. G. 1988, A Practical Guide, Academic Press).



[0019] When it is assumed that λ = dv/ du, a normal curved surface κ at an arbitrary point S (u, v) is expressed as Expression 2 in accordance with Expression 1.



[0020] The principle curvatures are κ1, κ2 which are obtained by solving the following Expression 3 and the Gaussian curvature K is defined as κ1 κ2. That is the definition of the conventional Gaussian curvature.



[0021] The definition of an extended Gaussian curvature (Λ) invented by the inventors of the present invention can be expressed by Expression 4.



[0022] That is, according to the method of the invention, first, Δκ1 = κ1' - κ1 and Δκ2 = κ2' - κ2 are obtained from the principle curvature (κ1', κ2') of the target curved surface S' and the principle curvature (κ1, κ2) of the reference curved surface S. (1) When Δ κ1≧0 and Δ κ2 ≧0, (a) it is determined that the two curvatures increase. (2) When Δ κ1≦0 and Δ κ2 ≦0, (b) it is determined that two curvatures decrease. (3) When Δκ1 · Δ κ2 < 0, (c) it is determined that one of the curvatures increases and the other decreases. (a), (b), and (c) are discriminated as mountain, valley, and twist, respectively, and are displayed by different symbols or colors on an image.

[0023] In other words,
  1. (1) when (κ1' - κ1)(κ'2 - κ2)≧0 and (κ'1 - κ1) ≧0, a label "mountain" is appended to Λ.
  2. (2) When (κ1' - κ1)( κ2' - κ2)≧0 and (κ1' - κ1) < 0, a label "valley" is appended to Λ.
  3. (3) When (κ1' - κ1))( κ2' - κ2)< 0, a label "twist" is appended to Λ.

[Embodiments]



[0024] Embodiments in which the method of the invention is applied will be described hereinbelow with reference to the drawings.

(Embodiment 1)



[0025] Fig. 5 is a perspective view of an impeller. Figs. 6 and 7 show the results of application of the method of the invention. Fig. 6 shows a case where a blade is bent and Fig. 7 shows a case where the blade is twisted by applying forces shown by the arrows.

[0026] In Figs. 6 and 7, the labels "mountain", "valley", and "twist" according to the invention are indicated by symbols "+", "-", and "%", respectively. In an actual image display, it is preferable to show the "mountain", "valley", and "twist" in colors such as brown, blue, and red, respectively. By displaying the determination of the shape error, a different part from the original three-dimensional curved face can be grasped easily and accurately by the different labels or different colors.

[0027] According to the method of the invention, the ratio of the same labels is calculated from the labels (a), (b), and (c) and is used as a coincidence ratio. That is, the labels of "mountain", "valley", and "twist" are mapped on pixels of a parameter plane ([0, 1] x [0, 1]) of u, v and the parameter plane is divided into lattice at a proper pitch (d). By applying a coincidence ratio (ψ) of Expression 5, how much the shapes coincide with each other can be shown by an objective numerical value. The coincidence ratio (ψ) indicates the relevance factor of the kind of the label. Since there is conventionally no index of discrepancy in a forming item, the coincidence ratio (ψ) can be used as a simple and clear index. It can be further developed and can be further finely divided in accordance with the magnitude of the extended Gaussian curvature (Λ).


(Embodiment 2)



[0028] Fig. 8 shows an application example of the method of invention on the basis of the measurement data of Fig. 4. Fig. 9 shows an application example of the method of the invention on the basis of the numerical simulation result of Fig. 2. That is, Figs. 8 and 9 show CMM data and FEM data, respectively. The data has been subjected to surface fitting of a solid modeler "DESIGNBASE" (manufactured by Richo) with an allowance of 0.002 mm from the original point. In this case, the coincidence ratio (ψ) is 50.23%.

[0029] Although the labels "mountain", "valley", and "twist" in the invention are shown by symbols of +, -, and %, respectively, also in Figs. 8 and 9, it is preferable that the mountain, valley, and twist are shown by colors such as brown, blue, and red, respectively, in an actual image display.

[0030] It will be understood from Figs. 8 and 9 that the difference between the actual forming shape (Fig. 8) or the simulation shape (Fig. 9) obtained by computer simulation and the reference shape, that is, a different part between the two three-dimensional shapes can be accurately grasped. In this regards, the method of the invention is more excellent than the conventional method shown in Figs. 2 to 4 by far.

[0031] Since the shape error can be determined by obtaining the principle curvatures of corresponding positions, the method can be applied without accurately coinciding the reference positions of two three-dimensional shapes with each other, and the cause of the shape error such as a partial curve can be found.

[0032] Further, by calculating the ratio of the same labels from the labels (a), (b), and (c) and using it as a coincidence ratio, how much the shapes coincide with each other as a whole can be objectively shown by a numerical value and an error can be easily discriminated even when the reference shape is complicated.

(Embodiment 3)



[0033] Figs. 10 and 11 show comparison examples with respect to the "side face" of Fig. 1. Fig. 10 shows an application example of the method of the invention on the basis of the measurement data. Fig. 11 shows an application example of the method of the invention on the basis of the result of the numerical simulation. The coincidence ratio (ψ) in this case is 52.47%.

POSSIBILITY OF INDUSTRIAL UTILIZATION



[0034] As mentioned above, the invention provides a simple and general method of defining a local shape error in a free-form surface. According to the method, the difference between the principle curvature of a curved surface including an error and that of a reference curved surface is used. The reference curved surface is usually expressed by CAD data. The curved surface including an error is obtained by approximating a group of discrete points such as measurement points and node points in numerical simulation. The principle curvature is used to evaluate a curved surface from the viewpoint of design, to form a curved surface, and the like, but is conventionally not used numerically to formulate or compare a shape error. The inventors of the present invention formulated a local shape error and performed FEM simulation and CMM data comparison by using the CAD data as a reference.

[0035] According to the method of the invention, "extended Gaussian curvature" is newly defined and a local shape error in a free-form surface is classified into three types (mountain, valley, and twist) by comparing the free-from surface with a CAD curved surface as a reference. The coincidence ratio calculating method using the image process technique is also proposed. With respect to the example of evaluating the forming discrepancy, the method of the invention was proved to be effective and mighty from the comparison of the deviation pattern of the actual measurement data and the numerical simulation data from the CAD data which is referred to.

[0036] By using the method of the invention, the discrepancy of the forming item which is formed by a press work such as the body of an automobile or the like can be evaluated objectively and numerically from, for example, the actual measurement data or the accuracy of the numerical simulation can be evaluated objectively and numerically. Consequently, the method of the invention can be widely industrially used with a three coordinate measuring machine, a CAD system, a CAM system, or a simulation system in the wide-ranged forming processing field or the shape measuring and evaluating field.

[0037] As mentioned above, according to the method of determining the shape error in the free-form surface of the invention, there are excellent effects such that the different part between two three-dimensional shapes can be accurately grasped, the shape error can be determined by obtaining the principle curvatures of the corresponding positions, the cause of the occurrence of an error such as a partial curve can be found, how much the shapes coincide with each other as a whole can be shown by an objective numerical value as a ratio, and the error can be easily determined even in a case where the reference shape is complicated.

[0038] Although the invention has been described in accordance with some preferred embodiments, it can be understood the range of the right included in the invention is not limited by the embodiments. On the contrary, the range of the right of the invention includes all of improvements, modifications, and equivalents included in the scope of the appended claims.


Claims

1. A method of determining a shape error in a free-form surface, comprising the steps of:

obtaining a principle curvature (K1' , K2') of a target curved surface S' and a principle curvature (K1. K2) of a corresponding position of a computer generated reference surface S; and

displaying each part on a display while classifying the difference between the principle curvatures into (a) a case where the two principle curvatures increase, (b) a case where the two principle curvatures decrease, and (c) a case where one of the principle curvatures increases and the other decreases.


 
2. The method of determining a shape error in a free-form surface according to claim 1, wherein Δκ1 = κ1' - κ1 and Δκ2 = κ2' - κ2 are obtained from the principle curvature (κ1', κ2') of the target curved surface S' and the principle curvature (κ1, κ2) at the corresponding position of the reference curved surface S,

(1) when Δκ1 ≥ 0 and Δκ2 ≥ 0, (a) it is determined that the two curvatures increase,

(2) when Δκ1 ≤ 0 and Δκ2 ≤ 0, (b) it is determined that the two curvatures decrease, and

(3) when Δ κ1 · Δ κ2 < 0, (c) it is determined that one of the principle curvatures increases and the other decreases.


 
3. The method of determining a shape error in a free-form surface according to claim 1 or 2, wherein (a), (b), and (c) are determined as labels mountain, valley, and twist, respectively, and are displayed in different symbols or colors on an image.
 
4. The method of determining a shape error in a free-form surface according to claim 3,wherein the ratio of the same labels is calculated from (a), (b), and (c) and is used as a coincidence ratio (ψ).
 
5. The method of determining a shape error in a free-form surface according to any one of claims 2 to 4, wherein the absolute value of Δκ1· Δκ2 is calculated for every label of (a), (b), and (c).
 


Ansprüche

1. Verfahren zur Bestimmung eines Formfehlers auf einer Freiformfläche, das die Schritte umfasst:

Erhalten einer Hauptkrümmung (k1', k2') einer gekrümmten Zielfläche S' und einer Hauptkrümmung (k1, k2) an einer entsprechenden Position einer computererzeugten Referenzfläche S; und

Darstellen eines jeden Bereiches auf einem Bildschirm, während der Unterschied zwischen den Hauptkrümmungen in (a) einen Fall, bei dem die beiden Hauptkrümmungen ansteigen, (b) einen Fall, bei dem die beiden Hauptkrümmungen abfallen, und (c) einen Fall, bei dem eine der Hauptkrümmungen ansteigt und die andere abfällt, unterteilt wird.


 
2. Verfahren zur Bestimmung eines Formfehlers auf einer Freiformfläche nach Anspruch 1, wobei Δ k1 = k1'- k1 und Δ k2 = k2'- k2 von der Hauptkrümmung (k1', k2') der gekrümmten Zielfläche S' und der Hauptkrümmung (k1, k2) an der entsprechenden Position der Referenzfläche S erhalten werden, und

(1) wenn Δ k1 ≥ 0 und Δ k2 ≥ 0, der Fall (a) eintritt, bei dem die beiden Krümmungen ansteigen,

(2) wenn Δ k1 ≤ 0 und Δ k2 ≤ 0, der Fall (b) eintritt, bei dem die beiden Krümmungen abfallen,

(3) wenn Δ k1 · Δ k2, der Fall (c) eintritt, bei dem eine der Hauptkrümmungen ansteigt und die andere abfällt.


 
3. Verfahren zur Bestimmung eines Formfehlers auf einer Freiformfläche nach Anspruch 1 oder 2, wobei (a), (b), und (c) jeweils mit den Bezeichnungen Berg, Tal, und Verwindung versehen und in unterschiedlichen Farben und mit unterschiedlichen Symbolen auf einem Bild dargestellt werden.
 
4. Verfahren zur Bestimmung eines Formfehlers auf einer Freiformfläche nach Anspruch 3, wobei das Verhältnis der gleichen Bezeichnungen aus (a), (b), und (c) berechnet und als ein Koinzidenzverhältnis (ψ) verwendet wird.
 
5. Verfahren zur Bestimmung eines Formfehlers auf einer Freiformfläche nach einem der Ansprüche 2 bis 4, wobei der Absolutwert von Δ k1 · Δ k2 für jede Bezeichnung von (a), (b), und (c) berechnet wird.
 


Revendications

1. Procédé de détermination d'une erreur de forme d'une surface à forme libre, comprenant les étapes suivantes :

obtenir une courbure de principe (K1', K2') d'une surface courbe cible S' et une courbure de principe (K1, K2) d'une position correspondante d'une surface de référence générée par ordinateur S ; et

afficher chaque partie sur un dispositif d'affichage tout en classifiant la différence entre les courbures de principe en (a) un cas dans lequel les deux courbures de principe augmentent, (b) un cas dans lequel les deux courbures de principe décroissent, et (c) un cas dans lequel une des courbures de principe augmente et l'autre décroît.


 
2. Procédé de détermination d'une erreur de forme d'une surface à forme libre selon la revendication 1, dans lequel ΔK1 = K1' -K1 et ΔK2 = K2' -K2 sont obtenus à partir de la courbure de principe (K1', K2') de la surface courbe cible S' et de la courbure de principe (K1, K2) à la position correspondante de la surface courbe de référence S,

(1) lorsque ΔK1 ≥ 0 et ΔK2 ≥ 0, (a) il est déterminé que les deux courbures augmentent,

(2) lorsque ΔK1 ≤ 0 et ΔK2 ≤ 0, (b) il est déterminé que les deux courbures décroissent,

(3) lorsque ΔK1 · ΔK2 < 0, (c) il est déterminé qu'une des courbures de principe augmente et l'autre décroît.


 
3. Procédé de détermination d'une erreur de forme d'une surface à forme libre selon la revendication 1 ou selon la revendication 2, dans lequel (a), (b) et (c) sont déterminés en tant qu'étiquettes montagne, vallée, et tournant, respectivement, et sont représentés par différents symboles ou couleurs sur une image.
 
4. Procédé de détermination d'une erreur de forme d'une surface à forme libre selon la revendication 3, dans lequel le rapport des étiquettes identiques est calculé à partir de (a), (b) et (c) et est utilisé comme rapport de coïncidence.
 
5. Procédé de détermination d'une erreur de forme d'une surface à forme libre selon l'une quelconque des revendications 2 à 4, dans lequel la valeur absolue de ΔK1 · ΔK2 est calculée pour chaque étiquette de (a), (b) et (c).
 




Drawing














Cited references

REFERENCES CITED IN THE DESCRIPTION



This list of references cited by the applicant is for the reader's convenience only. It does not form part of the European patent document. Even though great care has been taken in compiling the references, errors or omissions cannot be excluded and the EPO disclaims all liability in this regard.

Non-patent literature cited in the description