WO1998049524A1 - Procede permettant d'etablir une distinction entre des erreurs de forme d'une surface courbe a forme libre - Google Patents

Procede permettant d'etablir une distinction entre des erreurs de forme d'une surface courbe a forme libre Download PDF

Info

Publication number
WO1998049524A1
WO1998049524A1 PCT/JP1998/001842 JP9801842W WO9849524A1 WO 1998049524 A1 WO1998049524 A1 WO 1998049524A1 JP 9801842 W JP9801842 W JP 9801842W WO 9849524 A1 WO9849524 A1 WO 9849524A1
Authority
WO
WIPO (PCT)
Prior art keywords
shape
principal curvatures
principal
free
regions
Prior art date
Application number
PCT/JP1998/001842
Other languages
English (en)
French (fr)
Inventor
Kiwamu Kase
Akitake Makinouchi
Tetsuya Kondo
Naomichi Mori
Original Assignee
Riken
Amada Metrecs Company Limited
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Riken, Amada Metrecs Company Limited filed Critical Riken
Priority to JP2000571733A priority Critical patent/JP3937212B2/ja
Priority to AU70791/98A priority patent/AU7079198A/en
Priority to DE69838581T priority patent/DE69838581T2/de
Priority to US09/403,589 priority patent/US6683985B1/en
Priority to EP98917621A priority patent/EP0977009B1/en
Priority to CA002288177A priority patent/CA2288177C/en
Publication of WO1998049524A1 publication Critical patent/WO1998049524A1/ja

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/30—Polynomial surface description
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/20—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring contours or curvatures, e.g. determining profile
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
    • Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y10S715/00—Data processing: presentation processing of document, operator interface processing, and screen saver display processing
    • Y10S715/961—Operator interface with visual structure or function dictated by intended use
    • Y10S715/964—CAD or CAM, e.g. interactive design tools

Definitions

  • the present invention relates to a shape error determination method for a free-form surface. Description of related technology
  • FIG. 1 shows a forming sample showing a defect in plate forming
  • FIG. 2 shows an example of a conventional method for evaluating a curved surface shape error. This example is described in "Simulation of 3-D sheet bending process" (Takizawa, et al., 1991, VDI BE ICHETE NR.
  • the white portion indicates the reference shape (for example, a mold shape)
  • the mesh portion indicates the shape obtained by the molding simulation.
  • the simulation shape is displayed at the same position, and only the one located on the front side is displayed.Therefore, from the displayed white portion and the mesh portion, the difference between the simulation shape and the reference shape is obtained. The shape error can be roughly determined. And, this man The law had the following problems.
  • Fig. 3 shows the CMM data (approximately 400 points) measured using the three-dimensional measuring machine "Mityos-Par-BHN506".
  • Fig. 4 shows the top view of Fig. 3. (Approximately 800 points).
  • the shape of a molded product actually molded using a mold can be displayed as an image as shown in Figs. 3 and 4 by measuring with a 3D measuring machine. From the image, shape errors such as unevenness and twist can be roughly discriminated visually.
  • (4) when the reference shape is not flat but has a complicated curve three-dimensional measurement is possible. There was a problem that the difference from the result could hardly be distinguished visually.
  • the main object of the present invention is to accurately grasp a difference between two three-dimensional shapes such as an actual molded shape, a simulation shape by computer simulation, or a reference shape by CAD.
  • the reference position can be applied without exact match, part batchwise bend like can cause grasp an error occurs in the It is possible to grasp by objective numerical values how much the shape matches as a whole, and the reference shape is complicated. It is an object of the present invention to provide a method of determining a shape error of a free-form surface which can be easily determined even in such a case.
  • the inventors of the present invention have devised a new evaluation model “extended Gaussian curvature” which does not depend on a coordinate system. This classifies the local shape error of a free-form surface into three types (peak, valley, and twist) by comparing it with a reference CAD surface, for example. We also devised a method to calculate the ratio of the same label using image processing technology. The present invention is based on such a novel idea.
  • the principal curvatures at the corresponding positions of the target surface S ′ and the reference surface S are respectively obtained, and from the difference in the principal curvature, each part is (a) when the two principal curvatures are increased, b) When the two principal curvatures are decreasing, (c) when the one is increasing, and the other is decreasing, the classification is made and displayed.
  • the present invention further extends the Gaussian curvature K.
  • the principal curvatures at the corresponding positions of the target surface S ′ including the error and the reference surface S are respectively obtained, and the respective portions are expressed as (a) (b) ( c)
  • the classification is displayed. It is possible to accurately grasp the difference between the two three-dimensional shapes.
  • the shape error can be determined by obtaining the principal curvatures of the corresponding positions, the method can be applied without accurately matching the reference positions of the two three-dimensional shapes, such as partial bending. It is possible to understand the cause of the error occurrence.
  • Figure 1 is a halftone image displayed on the display of a processed sample showing a failure in plate forming.
  • FIG. 2 is a diagram showing an example of a conventional method for evaluating a curved surface shape error using a simulation.
  • FIG. 3 is a measurement diagram by a three-dimensional measuring machine.
  • FIG. 4 is a top view of FIG.
  • FIG. 5 is a halftone image displayed on a display showing a perspective view of the impeller.
  • FIG. 6 is a diagram in which bending deformation is added to the originally twisted shape according to the method of the present invention.
  • FIG. 7 is a diagram in which a twisted deformation is applied to an originally twisted shape according to the method of the present invention.
  • FIG. 8 is an example in which the method of the present invention is applied based on measured data.
  • FIG. 9 is an example in which the method of the present invention is applied based on the result of numerical simulation.
  • FIG. 10 shows another example in which the method of the present invention is applied based on measured data.
  • FIG. 11 shows another example in which the method of the present invention is applied based on the result of numerical simulation. Description of the preferred embodiment
  • Equation 1 is a relational expression in differential geometry. This relationship is disclosed, for example, in "Curves and Surfaces for Computer Aided Geometric Desigm (Far in, G, 1988, A Practial Guide. Academic Press).
  • Equation 2 the normal surface / c at an arbitrary point S (u, v) is expressed as in Equation 2 according to the description in Equation 1.
  • the principal curvature is / c ,, ⁇ 2, which is obtained by solving the following (Equation 3), and the Gaussian curvature K is defined as / c, ⁇ 2.
  • Equation 3 the Gaussian curvature K is defined as / c, ⁇ 2.
  • the above is the conventional definition of the Gaussian curvature K.
  • FIG. 5 is a perspective view of the impeller
  • FIGS. 6 and 7 show the results of applying the method of the present invention.
  • 6 and 7 show the case where the force shown by the arrow in the figure is applied
  • FIG. 6 shows the case where the blade (blade) is subjected to bending deformation
  • FIG. 7 shows the case where torsional deformation is applied.
  • the labels “peak” according to the present invention are indicated by “10”, “valleys” are indicated by “ ⁇ ”, and the twist is indicated by a symbol of%.
  • the ratio of the same label is calculated from the labels of (a), (b), and (c) above to obtain the matching ratio. That is, the labels of “mountain”, “valley”, and “twist” are mapped onto the pixels of the u, V parameter overnight plane ([0, 1] X [0, 1]), and this parameter plane is By dividing the grid with an appropriate pitch d and applying the collation rate ( ⁇ ) in (Equation 5), it is possible to grasp objectively how much the shape matches as a whole by numerical values. can do.
  • This collation rate ( ⁇ ) indicates the precision of the type of label, and since there was no index for molding defects in the past, it can be used as a simple and clear index. It should be noted that this can be further developed, and further subdivided according to the magnitude of the extended Gaussian curvature ( ⁇ ).
  • the method can be applied without accurately matching the reference positions of the two three-dimensional shapes, such as partial bending. It is possible to understand the cause of the error occurrence.
  • FIGS. 10 and 11 are comparative examples for the “side view” of FIG. 1
  • FIG. 10 is an example of applying the method of the present invention based on measurement data
  • FIG. 11 is a numerical simulation. This is an example in which the method of the present invention is applied based on the results of the review.
  • the present invention proposes a simple and general-purpose method for defining a local shape error of a free-form surface.
  • This method calculates the principal curvature between the surface containing the error and the reference surface. Use the difference.
  • the reference surface is usually represented by CAD data.
  • Surfaces with errors can be obtained by surface approximation of discrete points such as measurement points and node points in numerical simulation.
  • the principal curvature is used in the evaluation of curved surfaces and the generation of curved surfaces from the aspect of design, but was not used numerically in the formulation and comparison of shape errors.
  • the inventors of the present invention extended the Gaussian curvature, formulated the local shape error, and compared the FEM simulation with the CMM data based on the CAD data.
  • the “extended Gaussian curvature” is newly defined, and the local shape error of the free-form surface can be classified into three types (peak, valley, and twist) by comparison with a reference CAD surface. Classify into.
  • a matching method using image processing technology was also proposed.
  • the method of the present invention is proved to be effective and powerful through comparison of deviation patterns between actual measurement data and numerical simulation data from reference CAD data. Proven.
  • the method of the present invention can be widely used in a wide range of industrial fields, such as forming and measuring / evaluating fields, in combination with a three-dimensional measuring device, a CAD device, a CAM device, or a shy-mechanical device. it can.
  • the shape error determination method for a free-form surface can accurately grasp a difference between two three-dimensional shapes, and can determine a shape error by determining a main curvature at a corresponding position.
  • the cause of error occurrence such as partial bending, etc.
  • it has excellent effects such as easy identification.

Description

明細書 自由曲面の形状誤差判別方法 発明の背景 発明の技術分野
本発明は、 自由曲面の形伏誤差判別方法に関する。 関連技術の説明
自動車のボディ等のようにプレス加工により成形する成形物では、 成形の不具 合は従来習慣的に経験によって評価されてきた。 すなわち、 自由曲面の評価は、 従来、 主として 「目視」 で行われていた。 しかし、 近年、 コンピュータを用いた 設計 （ C A D) が普及し、 かつ加工時の変形をシユミ レーシヨ ンできるようにな つてきており、 これに伴い、 自由曲面の評価、 すなわち成形不具合を客観的に定 義し、 表示する手段が望まれている。 図 1 は板成形の不具合を示す成形加工サンプルであり、 図 2は、 従来の曲面形 状誤差の評価方法の一例を示す図である。 この例は、 "Simulation of 3-D sheet bending process" (Takizawa, et al.， 1991, VDI BE ICHETE NR. 894), "Some advances in FEM simulation of sheet metal forming processes using shell elements (Kawka et al. , 1995, Simulation of Materials Processing, Shen & DawsonCeds. ), Balkema, Rotterdam, 735-740) 等で報告されたシユミ レーシ ョ ンソフ ト （ I TA S— 3 D ) を使用した数値シユ ミ レ一ショ ンの結果である。 図 2 において、 白色部分は基準となる形状 （例えば金型形状） を示し、 メ ッ シ ュ部分は、 成形シュミ レーシヨ ンで得られた形状を示している。 また基準形状と シユミ レーショ ン形状とは同一位置に表示され、 前側に位置するもののみを表示 している。 従って、 表示された白色部分とメ ッ シュ部分から、 シユミ レ一シ ヨ ン 形状と基準形状との形状誤差をおおよそ判別することができる。 しかし、 この方 法では、 以下の問題点があった。
①基準形状とシユミ レーショ ン形状とを基準位置を決めて正確に一致させる必要 があるが、 その基準位置の決め方で結果が大きく影響を受ける。
②部分的な曲がりによりその他の部分の位置が大きく変位するため、 誤差発生の 原因把握が難しい。
③全体としてどの程度の割合で形状が一致しているかを客観的な数値で把握する ことができない。
図 3は 3次元測定機 「ミッ トヨス一パ一 B H N 5 0 6」 を用いて測定された C M Mデ一夕 （約 4 0 3 0 0点） を示しており、 図 4は図 3の上面図を示している (約 8 0 0 0点） 。 これらの図に示すように、 金型を用いて実際に成形した成形 品は、 3次元測定機で計測することにより、 その形状を図 3、 4のように画像表 示することができ、 この画像から、 凹凸、 ねじれ等の形状誤差を目視によりおお よそ判別することができる。 しかし、 この方法でも、 数値シユ ミ レ一シヨ ンにお ける上述した②③の問題点があると共に、 ④基準形状が平坦ではなく複雑な曲が りを有している場合には、 3次元計測結果との差異を目視ではほとんど判別でき ない問題点があつた。
上述したように、 形状誤差の評価のための実験、 計測、 結果の表示方法などは、 いまだ体系化されておらず、 成形不具合の指標となる単純で明確な定義がなく、 更に、 繰り返し可能な評価手法が従来存在しなかった。 発明の要約 本発明は、 かかる問題点を解決するために創案されたものである。 すなわち、 本発明の主目的は、 実際の成形形状、 コンピュータシユ ミ一シヨ ンによるシユミ —ショ ン形状、 或いは C A Dによる基準形状等の 2つの 3次元形状の相違部分を 的確に把握することができる自由曲面の形状誤差判別方法を提供することにある c また、 本発明の別の目的は、 基準位置を正確に一致させることなく適用でき、 部 分的な曲がり等の誤差発生の原因把握が可能であり、 全体と してどの程度の割合 で形状が一致しているかを客観的な数値で把握することができ、 基準形状が複雑 な場合でも容易に判別できる自由曲面の形状誤差判別方法を提供することにある。 本発明の発明者等は、 新しく、 座標系に依存しない評価モデル 「拡張ガウス曲 率」 を創案した。 これは、 例えば基準となる C A D曲面との比較により、 自由曲 面の局所的な形状誤差を 3つのタイプ （山、 谷、 ねじれ） に分類するものである。 また、 画像処理技術を利用して同一ラベルの比率を算出する手法も創案した。 本 発明は、 かかる新規の創案に基づく ものである。
すなわち本発明によれば、 対象曲面 S ' と基準曲面 Sの対応する位置の主曲率 をそれぞれ求め、 主曲率の差から、 各部分を （ a ) 2つの主曲率が増加している 場合、 （b ) 2つの主曲率が減少している場合、 （ c ) 一方が増加、 他方が減少 している場合に分類して表示する、 ことを特徴とする自由曲面の形状誤差判別方 法が提供される。
本発明の好ましい方法によれば、 対象曲面 S ' の主曲率 （ /c ， /c 2' ) と、 対 応する位置の基準曲面 Sの主曲率 （ ι , κ 2 ) とから、 Δ /c！ = /c ,' - /c 1 、 Δ κ 2 = Λ 2' - c 2 を求め、 ① ≥ 0かつ Δ κ 2 ≥ 0の場合に、 （ a ) 2つ の主曲率が増加していると判別し、 ② Δ /c , ≤ 0かっ 0 2 ≤ 0の場合に、 （b ) 2つの主曲率が減少していると判別し、 ③ Δ /c , · Δ / 2 < 0の場合に、 （ c ) 一方が増加、 他方が減少していると判別する。 また、 （ a ) を山、 （b ) を谷、 ( c ) をねじれと判別し、 それぞれ異なる記号又は色で画像表示する、 ことが好 ましい。 更に、 上記 （ a ) ( b ) ( c ) のラベルから、 同一ラベルの比率を算出 して照合率とする、 ことが好ましい。
ガウス曲率 Kは、 3次元曲面の主曲率/ , κ 2 の積 κ 2 であり、 ®K > 0の場合に、 楕円的（elliptic)であり、 ② K= 0の場合に、 放物的（parabol ic) であり、 ®K < 0の場合に、 双曲的 （hyperbolic, 鞍型） であることがわかる。 本発明はこのガウス曲率 Kを更に拡張したものである。 すなわち、 本発明の方 法によれば、 誤差を含む対象曲面 S ' とその基準曲面 Sの対応する位置の主曲率 をそれぞれ求め、 主曲率の差から、 各部分を （ a ) ( b ) ( c ) に分類して表示 することにより、 2つの主曲率が増加している場合、 2つの主曲率が減少してい る場合、 一方が増加、 他方が減少している場合に分類して表示することができ、 2つの 3次元形状の相違部分を的確に把握することができる。 また、 この方法によれば、 対応する位置の主曲率をそれぞれ求めことにより形 状誤差を判別できるので、 2つの 3次元形状の基準位置を正確に一致させること なく適用でき、 部分的な曲がり等の誤差発生の原因把握が可能である。
更に、 （ a ) ( b ) ( c ) のラベルから、 同一ラベルの比率を算出して照合率 とすることにより、 全体としてどの程度の割合で形状が一致しているかを客観的 な数値で把握することができ、 かつ基準形状が複雑な場合でも容易に判別するこ とができる。
本発明のその他の目的及び有利な特徴は、 添付図面を参照した以下の説明から 明らかになろう。 図面の簡単な説明
図 1 は、 板成形の不具合を示す成形加工サンプルのディスプレー上に表示した 中間調画像である。
図 2は、 シユミ レ一シヨ ンを用いた従来の曲面形状誤差の評価方法の一例を示 す図である。
図 3は、 3次元測定機による測定図である。
図 4は、 図 3の上面図である。
図 5は、 インペラの斜視図を示すディ スプレー上に表示した中間調画像である 図 6は、 本発明の方法による元々ねじれている形状に曲げ変形が加えられた図 である。
図 7は、 本発明の方法による元々ねじれている形状にねじれ変形が加えられた 図である。
図 8は、 測定データを基に本発明の方法を適用した例である。
図 9は、 数値シユミ レ一ショ ン結果を基に本発明の方法を適用した例である。 図 1 0は、 測定データを基に本発明の方法を適用した別の例である。
図 1 1 は、 数値シユミ レーショ ン結果を基に本発明の方法を適用した別の例で める。 好ま しい実施例の説明
先ず、 本発明の方法の原理を説明する。
自由曲面 S = S ( u， V ) はパラメ一タ u , Vによって表現される。 （式 1 ) は、 微分幾何学における関係式である。 この関係は、 例えば、 "Curves and Surf aces for Computer Aided Geometric Desigm (Far in, G, 1988， A Practial Gu i de. Academic Press) に開示されている。
【数 1】
E = S„S„， = S,,S、,，G = S,,Sv,i = nS,
(式 1)
M = nSm,,N = nS„ S„ xS,
λ = d vノ d uとすると、 任意の点 S ( u， v ) における法曲面 / c は、 式 1 の 記述に従って式 2のように表現される。
【数 2】
L + 2Aa + NA2
(式 2)
E + 2FX + GXl
主曲率は / c , ， κ 2 であり、 以下の （式 3 ) を解く ことにより求められ、 ガウ ス曲率 Kは / c , κ 2 と定義される。 以上が、 従来のガウス曲率 Kの定義である。
【数 3】
NE— 1MF+LG , LN - M1 n
κ - κ + = 0 (式 3)
EG-F2 EG-F1 次に、 本発明の発明者等が創案した拡張ガウス曲率 （Λ) の定義は、 （式 4 ) で表現することができる。
【数 4】
V -κ )(κ2 -κ2)
-κ2 ))
(式 4) た
Figure imgf000008_0001
すなわち、 本発明の方法によれば、 対象曲面 S ' の主曲率 （ /C ， κ 2 ) と、 対応する位置の基準曲面 Sの主曲率 （ /c , ， κ 2 ) とから、 まず、 Δ , = κ Γ - κ 1 、 Δ κ 2 = c 2' - c 2 を求め、 ① ≥ 0かつ A /c 2 ≥ 0の場合に、 ( a ) 2つの主曲率が増加していると判別し、 ② Δ ≤ 0かっ八 2 Oの場 合に、 （b ) 2つの主曲率が減少していると判別し、 ③△ , · Δ /c 2 く 0の場 合に、 （ c ) 一方が増加、 他方が減少していると判別し、 （ a ) を山、 （b) を 谷、 （ c ) をねじれと判別し、 それぞれ異なる記号又は色で画像表示する。 言い換えれば、
① （ i ' - κ I ) ( κ 2 ' - κ 2 ) ≥ 0かつ ' - κ 1 ) ≥ 0のとき、 ラ ベル 「山」 が Λに付加される。
② （ l ' - K 1 ) ( K 2 ' - K 2 ) ≥ 0かつ ' - Κ 1 ) く 0のとき、 ラ ベル 「谷」 が Λに付加される。
③ （ / , ' - κ 1 ) ( κ 2 ' - κ 2 ) く 0のとき、 らベる 「ねじれ」 が Λに付加 される。 【実施例】
以下、 本発明の方法を適用した実施例を図面を参照して説明する。
(実施例 1 )
図 5は、 イ ンペラの斜視図であり、 図 6及び図 7は、 本発明の方法を適用した 結果を示している。 なお図 6及び図 7は図中に矢印で示す力を作用させて、 図 6 では翼 （ブレー ド） に曲げ変形、 図 7ではねじれ変形を与えた場合である。 図 6及び図 7において、 本発明によるラベル 「山」 を十、 「谷」 を—、 ねじれ を％の記号で表示している。 なお、 実際の画像表示では、 「山」 を茶色、 「谷」 を水色、 ねじれを赤色、 等に色分けすることが好ま しい。 この形状誤差の判別表 示により、 元の 3次元曲面との相違部分をラベルの相違又は色の違いにより簡単 かつ的確に把握することができる。
本発明の方法によれば、 上記 （ a ) (b ) ( c ) のラベルから、 同一ラベルの 比率を算出して照合率とするのがよい。 すなわち、 「山」 「谷」 「ねじれ」 のラ ベルを、 u, Vのパラメ一夕平面 （ [ 0 , 1 ] X [ 0 , 1 ] ) の画素上にマップ し、 このパラメ一タ平面を、 適当なピッチ dにより格子に分割して、 （式 5 ) の 照合率 （ψ) を適用することにより、 全体と してどの程度の割合で形状が一致し ているかを客観的な数値で把握することができる。 この照合率 （Ψ) は、 ラベル の種類の適合率をあらわしており、 従来成形不具合の指標となるものが全く なか つたことから、 単純で明確な指標と して用いることができる。 なお、 これを更に 発展させて、 拡張ガウス曲率 （Λ) の大きさにより、 更に細かく 区分すること も できる。
【数 5】
Γ - {n^ber_of_the_Same_labeJ) m{n ) (式 5)
Figure imgf000009_0001
' (実施例 2 )
図 8は図 4の測定データを基に本発明の方法を適用した例であり、 図 9は、 図 2の数値シュミ レ一ショ ン結果を基に本発明の方法を適用した例である。 すなわ ち図 8 と図 9はそれぞれ、 C M Mデータと F E Mデータを示している。 これらは ソリ ツ ドモデラ （ R i c 0 h D E S I G N B A S E ) によって、 元の点に対し て、 0 . 0 0 2 m mの許容誤差で曲面フィティ ングされた。 この場合、 照合率は Ψ = 5 0 . 2 3 %となった。
なお、 図 8、 9においても、 本発明によるラベル 「山」 を十、 「谷」 を一、 ね じれを％の記号で表示しているが、 実際の画像表示では、 「山」 を茶色、 「谷」 を水色、 ねじれを赤色、 等に色分けするのがよい。
図 8及び図 9から、 実際の成形形状 （図 8 ) 或いはコ ンピュータシユミ一ショ ンによるシュミーショ ン形状 （図 9 ) とそれらの基準となる基準形状との形状誤 差、 すなわち 2つの 3次元形状の相違部分を的確に把握することができることが わかる。 この点において、 図 2乃至図 4に示した従来の方法に比べて格段に優れ ている。
また、 この方法によれば、 対応する位置の主曲率をそれぞれ求めことにより形 状誤差を判別できるので、 2つの 3次元形状の基準位置を正確に一致させること なく適用でき、 部分的な曲がり等の誤差発生の原因把握が可能である。
更に、 （ a ) ( b ) ( c ) のラベルから、 同一ラベルの比率を算出して照合率 とすることにより、 全体としてどの程度の割合で形状が一致しているかを客観的 な数値で把握することができ、 かつ基準形状が複雑な場合でも容易に判別するこ とができる。
(実施例 3 )
図 1 0、 図 1 1 は、 図 1 の 「側面」 に対する比較例であり、 図 1 0は測定デ一 夕を基に本発明の方法を適用した例であり、 図 1 1 は、 数値シユミ レ一シヨ ン結 果を基に本発明の方法を適用した例である。 この場合の照合率は、 Ψ = 5 2 · 4 7 %となった。
上述したように、 本発明は、 単純で汎用的な、 自由曲面の局所的な形状誤差の 定義方法を提案している。 この方法は、 誤差を含む曲面と基準曲面との主曲率の 差を使用する。 基準曲面は、 通常 C A Dデータで表現される。 誤差を含む曲面は、 測定点や数値シュミ レーショ ンにおけるノー ド点などの離散的な点群を曲面近似 することによって得られる。 主曲率は、 意匠的な側面から曲面の評価や曲面の生 成などで使用されるが、 形状誤差の定式化や比較などで数値的に使用されること はなかった。 本発明の発明者等は、 ガウス曲率を拡張し、 局所的な形状誤差を定 式化し、 かつ C A Dデータを基準にして、 F E Mシュミ レーシヨ ンと C M Mデー 夕の比較を行った。
上述した本発明の方法によれば、 「拡張ガウス曲率」 を新たに定義し、 基準と なる C A D曲面との比較により、 自由曲面の局所的な形状誤差を 3つのタイプ （ 山、 谷、 ねじれ） に分類する。 また、 画像処理技術を使用する照合手法も提案さ れた。 成形不具合の評価例では、 参照される C A Dデータからの、 実際の測定デ —タと数値シュミ レーショ ンデ一夕の偏差パターンの比較を通して、 本発明の方 法が、 有効かつ強力であることが証明された。
本発明の方法を用いることにより、 自動車のボディ等のようにプレス加工によ り成形する成形物の不具合を、 例えば実際の測定データから客観的 ·数値的に評 価することができ、 或いは数値シユミ レ一シヨ ンの精度を客観的 ·数値的に評価 することができる。 従って、 本発明の方法は、 広範囲の成形加工分野や形状測定 •評価分野において、 3次元測定機、 C A D装置、 C A M装置、 或いはシユミ レ —シヨ ン装置との組み合わせにより、 広く産業上役立てることができる。
上述したように、 本発明の自由曲面の形状誤差判別方法は、 2つの 3次元形状 の相違部分を的確に把握することができ、 対応する位置の主曲率をそれぞれ求め ことにより形状誤差を判別でき、 部分的な曲がり等の誤差発生の原因把握が可能 であり、 全体としてどの程度の割合で形状が一致しているかを客観的な数値で把 握することができ、 かつ基準形状が複雑な場合でも容易に判別することができる 等の優れた効果を有する。
なお、 本発明をいくつかの好ま しい実施例により説明したが、 本発明に包含さ れる権利範囲は、 これらの実施例に限定されないことが理解できょう。 反対に、 本発明の権利範囲は、 添付の請求の範囲に含まれるすべての改良、 修正及び均等 物を含むものである。

Claims

請求の範囲
1. 対象曲面 S ' と基準曲面 Sの対応する位置の主曲率をそれぞれ求め、 主曲率の差から、 各部分を （ a ) 2つの主曲率が増加している場合、 （b ) 2つ の主曲率が減少している場合、 （ c ) 一方が増加、 他方が減少している場合に分 類して表示する、 ことを特徴とする自由曲面の形状誤差判別方法。
2. 対象曲面 S ' の主曲率 （ /c ， /c 2' ) と、 対応する位置の基準曲面 S の主曲率 （ V 1 ， ん 2 ) とから、
Δ Κ 1 = Κ 1 ― Κ 1 Δ ί 2 ― Κ 2 ― Κ 2 を求め、
① Δ Λ; , ≥ 0かっ厶 !；2 ≥ 0の場合に、 （ a ) 2つの主曲率が增加していると判 別し、
② ≤ 0かっ厶 2 Oの場合に、 （b ) 2つの主曲率が減少していると判 別し、
③ · Δ /c 2 < 0の場合に、 （ c ) 一方が増加、 他方が減少していると判別 する、 ことを特徴とする請求項 1 に記載の自由曲面の形状誤差判別方法。
3. ( a ) を山、 （b ) を谷、 （ c ) をねじれと判別し、 それぞれ異なる 記号又は色で画像表示する、 ことを特徴とする請求項 1又は 2に記載の自由曲面 の形状誤差判別方法。
4. 上記 （ a ) ( b ) ( c ) のラベルから、 同一ラベルの比率を算出して 照合率とする、 ことを特徴とする請求項 1乃至 3に記載の自由曲面の形状誤差判 別方法。
5. Δ κ 1 · Δ / 2 の絶対値を前記 （ a ) ( b ) ( c ) のラベル毎に計算 し定量化する、 ことを特徴とする請求項 2乃至 4に記載の自由曲面の形状誤差判 別方法。
PCT/JP1998/001842 1997-04-25 1998-04-22 Procede permettant d'etablir une distinction entre des erreurs de forme d'une surface courbe a forme libre WO1998049524A1 (fr)

Priority Applications (6)

Application Number Priority Date Filing Date Title
JP2000571733A JP3937212B2 (ja) 1997-04-25 1998-04-22 自由曲面の形状誤差判別方法
AU70791/98A AU7079198A (en) 1997-04-25 1998-04-22 Method of discriminating shape errors of free-form curved surface
DE69838581T DE69838581T2 (de) 1997-04-25 1998-04-22 Verfahren zur determinierung von fehlern im form einer willkürlichen krummen oberflache
US09/403,589 US6683985B1 (en) 1997-04-25 1998-04-22 Method of discriminating shape of free-form curved surface
EP98917621A EP0977009B1 (en) 1997-04-25 1998-04-22 Method of discriminating shape errors of free-form curved surface
CA002288177A CA2288177C (en) 1997-04-25 1998-04-22 Method of determining shape error of free-form surface

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP10913997 1997-04-25
JP9/109139 1997-04-25

Publications (1)

Publication Number Publication Date
WO1998049524A1 true WO1998049524A1 (fr) 1998-11-05

Family

ID=14502595

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP1998/001842 WO1998049524A1 (fr) 1997-04-25 1998-04-22 Procede permettant d'etablir une distinction entre des erreurs de forme d'une surface courbe a forme libre

Country Status (7)

Country Link
US (1) US6683985B1 (ja)
EP (1) EP0977009B1 (ja)
JP (1) JP3937212B2 (ja)
AU (1) AU7079198A (ja)
CA (1) CA2288177C (ja)
DE (1) DE69838581T2 (ja)
WO (1) WO1998049524A1 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000161946A (ja) * 1998-11-27 2000-06-16 Mitsutoyo Corp ３次元データ解析方法及び装置
JP2006201164A (ja) * 2005-01-18 2006-08-03 General Electric Co <Ge> 加工された構成要素の検査方法およびシステム
JP2010169395A (ja) * 2008-07-08 2010-08-05 Samsung Heavy Ind Co Ltd 曲形部材の加工完成度評価システム及びその方法
JP2011513727A (ja) * 2008-02-28 2011-04-28 コーニング インコーポレイテッド シート材の基準面に対する共形化可能度を予測する方法
JP2019090675A (ja) * 2017-11-14 2019-06-13 株式会社荏原製作所 検査装置、検査システム、プログラムを記憶した記憶媒体

Families Citing this family (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3344649B2 (ja) * 1998-08-19 2002-11-11 理化学研究所 自由曲面の形状誤差評価方法
CA2314794A1 (en) * 2000-08-01 2002-02-01 Dimitre Hristov Apparatus for lesion or organ localization
US20020065637A1 (en) * 2000-11-30 2002-05-30 Thompson John S. Method and apparatus for simulating the measurement of a part without using a physical measurement system
US7139432B2 (en) * 2002-04-10 2006-11-21 National Instruments Corporation Image pattern matching utilizing discrete curve matching with a mapping operator
JP4794160B2 (ja) * 2004-03-31 2011-10-19 日東電工株式会社 表面形状測定装置および表面形状測定方法
US7672705B2 (en) * 2004-07-19 2010-03-02 Resonant Medical, Inc. Weighted surface-to-surface mapping
US7729744B2 (en) * 2004-07-20 2010-06-01 Resonant Medical, Inc. Verifying lesion characteristics using beam shapes
WO2006007716A2 (en) * 2004-07-20 2006-01-26 Resonant Medical Inc. Calibrating imaging devices
CA2581127C (en) * 2004-09-20 2013-12-24 Resonant Medical Inc. Radiotherapy treatment monitoring using ultrasound
TWI267761B (en) * 2004-12-03 2006-12-01 Hon Hai Prec Ind Co Ltd A simulative measurement demo system and method
EP1940515A4 (en) * 2005-09-06 2010-05-26 Resonant Medical Inc SYSTEM AND METHOD FOR CONFIGURING PATIENT DATA FOR RADIOTHERAPY TREATMENT
US8929621B2 (en) * 2005-12-20 2015-01-06 Elekta, Ltd. Methods and systems for segmentation and surface matching
US8395626B2 (en) * 2006-04-08 2013-03-12 Alan Millman Method and system for interactive simulation of materials
US8786613B2 (en) 2006-04-08 2014-07-22 Alan Millman Method and system for interactive simulation of materials and models
US9451928B2 (en) * 2006-09-13 2016-09-27 Elekta Ltd. Incorporating internal anatomy in clinical radiotherapy setups
US10531858B2 (en) 2007-07-20 2020-01-14 Elekta, LTD Methods and systems for guiding the acquisition of ultrasound images
WO2009012577A1 (en) * 2007-07-20 2009-01-29 Resonant Medical Inc. Methods and systems for compensating for changes in anatomy of radiotherapy patients
US8135198B2 (en) * 2007-08-08 2012-03-13 Resonant Medical, Inc. Systems and methods for constructing images
US8189738B2 (en) * 2008-06-02 2012-05-29 Elekta Ltd. Methods and systems for guiding clinical radiotherapy setups
US10542962B2 (en) * 2009-07-10 2020-01-28 Elekta, LTD Adaptive radiotherapy treatment using ultrasound
US9248316B2 (en) 2010-01-12 2016-02-02 Elekta Ltd. Feature tracking using ultrasound
US20110172526A1 (en) 2010-01-12 2011-07-14 Martin Lachaine Feature Tracking Using Ultrasound
FR2966956B1 (fr) * 2010-10-27 2012-12-14 Michelin Soc Tech Methode de pre-traitement d'une image tri dimensionnelle de la surface d'un pneumatique a l'aide de deformations b-spline successives
US20130085734A1 (en) * 2011-10-03 2013-04-04 Toyota Motor Engineering & Manufacturing North America, Inc. Model and computer based coolant flow diagnostic system
JP6528308B2 (ja) * 2015-02-05 2019-06-12 国立大学法人神戸大学 形状評価方法および形状評価装置
CN110769731B (zh) * 2017-06-15 2022-02-25 奥林巴斯株式会社 内窥镜系统、内窥镜用处理系统、图像处理方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04276513A (ja) * 1991-03-04 1992-10-01 Matsushita Electric Ind Co Ltd 形状測定処理方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4679331A (en) * 1985-08-26 1987-07-14 Ppg Industries, Inc. Apparatus and method for determining contour characteristics of a contoured article
JPH03176703A (ja) * 1989-12-05 1991-07-31 Yoshiaki Kakino 数値制御装置
US5715166A (en) * 1992-03-02 1998-02-03 General Motors Corporation Apparatus for the registration of three-dimensional shapes
US5844801A (en) * 1994-12-08 1998-12-01 Honda Giken Kogyo Kabushiki Kaisha Method of inspecting and manufacturing vehicle body
JP3054108B2 (ja) * 1997-08-15 2000-06-19 理化学研究所 自由曲面の測定データ合成方法
US6285320B1 (en) * 1999-09-03 2001-09-04 Sikorsky Aircraft Corporation Apparatus and method for mapping surfaces of an object

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04276513A (ja) * 1991-03-04 1992-10-01 Matsushita Electric Ind Co Ltd 形状測定処理方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000161946A (ja) * 1998-11-27 2000-06-16 Mitsutoyo Corp ３次元データ解析方法及び装置
JP2006201164A (ja) * 2005-01-18 2006-08-03 General Electric Co <Ge> 加工された構成要素の検査方法およびシステム
JP2011513727A (ja) * 2008-02-28 2011-04-28 コーニング インコーポレイテッド シート材の基準面に対する共形化可能度を予測する方法
JP2010169395A (ja) * 2008-07-08 2010-08-05 Samsung Heavy Ind Co Ltd 曲形部材の加工完成度評価システム及びその方法
JP2019090675A (ja) * 2017-11-14 2019-06-13 株式会社荏原製作所 検査装置、検査システム、プログラムを記憶した記憶媒体

Also Published As

Publication number Publication date
EP0977009B1 (en) 2007-10-17
CA2288177C (en) 2006-07-18
DE69838581D1 (de) 2007-11-29
CA2288177A1 (en) 1998-11-05
EP0977009A4 (en) 2006-10-18
DE69838581T2 (de) 2008-02-07
US6683985B1 (en) 2004-01-27
JP3937212B2 (ja) 2007-06-27
EP0977009A1 (en) 2000-02-02
AU7079198A (en) 1998-11-24

Similar Documents

Publication Publication Date Title
WO1998049524A1 (fr) Procede permettant d&#39;etablir une distinction entre des erreurs de forme d&#39;une surface courbe a forme libre
US6353768B1 (en) Method and apparatus for designing a manufacturing process for sheet metal parts
Kase et al. Shape error evaluation method of free-form surfaces
JP4633625B2 (ja) メタルシート成形段階のジオメトリモデルの決定
US6785640B1 (en) Surface evaluation in a stamping manufacturing process utilizing true reflection line methodology and computer graphics technology
CN101253385A (zh) 应变评价装置及应变评价方法
CN102735204B (zh) 一种基于弦线的航空薄壁叶片加工扭曲度误差测量方法
US20100114350A1 (en) Method of determining mesh data and method of correcting model data
US6338000B1 (en) Method of generating shape data method of verifying shape data
CN100589126C (zh) 利用阻光度传递函数的ct图像体素成像方法
Hongzhi et al. Investigation of sheet metal forming by numerical simulation and experiment
JPH08339396A (ja) 金属板の変形過程の数値シミュレート結果の処理装置
JP4604762B2 (ja) 特徴図の比較評価方法およびシミュレーションモデル作成方法
US20210407064A1 (en) Method and device for geometric analysis of a part surface
JP2016203216A (ja) 鍛造型の補正方法
Eder et al. Validation of material models for sheet metals using new test equipment
CN111914350A (zh) 一种飞机框肋类钣金件弯边成形性判定方法
CN107609255A (zh) 车身特征线消失位置的检测方法
JP2021508290A (ja) 曲げるべきシート工作物の曲げ短縮の決定
JP4097164B2 (ja) 立体形状記述方法
JP2000122996A (ja) ワーク成形面の形状不良評価方法
JP2007299167A (ja) 設計検証装置
Carasusán et al. An automated procedure for non-contact strain analysis of sheet metal parts
CN117367516B (zh) 基于多维测试的铝塑复合板性能检测系统
JP3945311B2 (ja) プレス成形品寸法精度評価方法およびその装置

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AL AM AT AU AZ BB BG BR BY CA CH CN CZ DE DK EE ES FI GB GE HU IS JP KE KG KR KZ LK LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK TJ TM TR TT UA UG US UZ VN

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW SD SZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN ML MR NE SN TD TG

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 1998917621

Country of ref document: EP

ENP Entry into the national phase

Ref document number: 2288177

Country of ref document: CA

Ref country code: CA

Ref document number: 2288177

Kind code of ref document: A

Format of ref document f/p: F

WWE Wipo information: entry into national phase

Ref document number: 09403589

Country of ref document: US

WWP Wipo information: published in national office

Ref document number: 1998917621

Country of ref document: EP

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

ENP Entry into the national phase

Ref country code: JP

Ref document number: 2000 571733

Kind code of ref document: A

Format of ref document f/p: F

WWG Wipo information: grant in national office

Ref document number: 1998917621

Country of ref document: EP