(57)【要約】【課題】離散点で表された曲線の各点の曲率及び各点間の弧長を入力とし、それに対応する座標を簡単な処理で出力すること。【解決手段】二次元直交座標系の点列Ｐ（１）、Ｐ（２）、・・・Ｐ（ｎ）のそれぞれに対して弧長（点列の間隔）と曲率が既知であり、かつ、ｘｙ座標が未知である場合、隣接する３点が同一の円上に位置すると仮定し、幾何学的な手法（三平方の定理及び余弦定理）を用いることによって、それぞれの点列のｘｙ座標を算出する。【選択図】図２
2018151774.tif 000002
【特許請求の範囲】
【請求項１】
離散点で表された曲線の各点の曲率、前記各点間の弧長を入力する入力部と、
前記曲率に基づき円を構成し、前記弧長に基づき前記円上で座標を特定する座標算出部と、
前記座標算出部で算出された座標を出力する出力部と
を具備する曲線座標生成装置。
【請求項２】
請求項１に記載の曲線座標生成装置であって、
前記入力部は、前記離散点の第１の座標及び前記第１の座標の次の第２の座標を入力し、
前記座標算出部は、前記第１の座標、前記第２の座標及び前記第２の座標における前記曲率に基づき円を構成し、前記第２の座標から、前記第２の座標と前記第２の座標の次の第３の座標となる点との間の前記弧長だけ離れた前記円上の座標を第３の座標として算出し、
前記出力部は、前記座標算出部で算出された第３の座標を出力する
曲線座標生成装置。
【請求項３】
請求項２に記載の曲線座標生成装置であって、
前記第１の座標は、前記離散点で表された曲線における始点である
曲線座標生成装置。
【請求項４】
請求項３に記載の曲線座標生成装置であって、
前記入力部は、前記離散点の終点の座標及び前記曲率の修正のための重みを入力し、
前記座標算出部は、前記入力した重みに基づき仮想の各点の座標を算出し、仮想の終点の座標が前記入力した終点の座標と一致するように、前記入力した重みにパラメータを付与し、前記パラメータを付与した重みに基づき前記各点の曲率及び座標を算出する
曲線座標生成装置。
【請求項５】
請求項１から４のうちいずれか１項に記載の曲線座標生成装置であって、
前記各点で曲率を指定でき、かつ、前記各点間の曲率が単調となる曲線生成法を用いて、前記出力部により出力された座標の点列を通過する曲線を生成する曲線生成部
を更に具備する曲線座標生成装置。
【請求項６】
請求項１から５のうちいずれか１項に記載の曲線座標生成装置であって、
前記出力部により出力された座標の点列を表示する表示部と、
前記表示部に表示された座標の点列における各点の曲率及び弧長を前記表示部上で可変操作する操作部と
を更に具備する曲線座標生成装置。
【請求項７】
曲線を有する物体を生産する方法であって、
離散点で表された前記曲線の各点の曲率、前記各点間の弧長を入力し、
前記曲率に基づき円を構成し、前記弧長に基づき前記円上で座標を特定し、
前記特定された座標に基づき前記物体の曲線を設計し、前記物体を生産する
曲線を有する物体の生産方法。
【請求項８】
離散点で表された曲線の各点の曲率、前記各点間の弧長を入力するステップと、
前記曲率に基づき円を構成し、前記弧長に基づき前記円上で座標を特定するステップと、
前記特定された座標を出力するステップと
をコンピュータシステムに実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、例えば航空機の翼の設計に用いられる曲線座標生成装置、翼などの曲線を有する物体の生産方法及びプログラムに関する。
【背景技術】
【０００２】
航空機の翼の設計など、曲線を取り扱う技術分野においては、曲率に基づいた設計手法が必要とされる。その場合に曲率制御を可能とする曲線や点列を生成する技術が必要である。
【０００３】
非特許文献１には、離散曲率を入力として離散点列を出力する技術が記載されている。この技術は、座標Ｐｉ−2及び座標Ｐｉ−１の曲率から円を描き、この円上で次の座標Ｐｉのｘ座標に対応するｙ座標を求めることで、次の座標Ｐｉを求め、以下同様の演算を繰り返すことで、離散点で表された曲線の座標を求める。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Barger, R. L. and Brooks, C. W.: A Streamline Curvature Method for Design of Supercritical Airfoils, NASA TN D-7770, 1974.
【発明の概要】
【発明が解決しようとする課題】
【０００５】
典型的には、流体力学の分野において、空力性能向上のため弧長をパラメータとする形状曲率制御を行う必要があるが、非特許文献１では弧長による制御ができない。
【０００６】
また、非特許文献１では、ｘ座標に対してｙ座標が多価関数の場合には使えない。
【０００７】
以上のような事情に鑑み、本発明の目的は、離散点で表された曲線の各点の曲率及び各点間の弧長を入力とし、それに対応する座標を簡単な処理で出力することができる曲線座標生成装置及びプログラムを提供することにある。
【０００８】
本発明は、このような手法を用いて、翼などの曲線を有する物体の生産方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記目的を達成するため、本発明の一形態に係る曲線座標生成装置は、離散点で表された曲線の各点の曲率、前記各点間の弧長を入力する入力部と、前記曲率に基づき円を構成し、前記弧長に基づき前記円上の座標を算出する座標算出部と、前記座標算出部で算出された座標を出力する出力部とを具備する。
【００１０】
本発明の一形態に係る曲線座標生成装置は、離散点で表された曲線の曲率、前記離散点間の弧長、前記離散点の第１の座標及び前記第１の座標の次の第２の座標を入力する入力部と、前記第１の座標、前記第２の座標及び前記第２の座標における前記曲率に基づき円を構成し、前記第２の座標から、前記第２の座標と前記第２の座標の次の第３の座標となる点との間の前記弧長だけ離れた前記円上の座標を第３の座標として算出する座標算出部と、前記座標算出部で算出された第３の座標を出力する出力部とを具備する。
【００１１】
前記第１の座標は、前記離散点で表された曲線における始点である。
本発明では、第１の座標、第２の座標及び第２の座標における曲率に基づき円を構成し、第２の座標から、第２の座標と第３の座標となる点との間の弧長だけ離れた円上の座標を第３の座標として算出しているので、離散点で表された曲線の曲率及び離散点間の弧長を入力とし、それに対応する座標を簡単な処理で出力することができる。
【００１２】
前記入力部は、前記離散点の終点の座標及び前記曲率の修正のための重みを入力し、前記座標算出部は、前記入力した重みに基づき仮想の各点の座標を算出し、仮想の終点の座標が前記入力した終点の座標と一致するように、前記入力した重みにパラメータを付与し、前記パラメータを付与した重みに基づき前記各点の曲率及び座標を算出する。これにより、曲率に基づいた点列座標の修正が可能であると共に、曲率修正量の重みの入力によって修正後の各点列の曲率分布の制御が可能である。
【００１３】
本発明の一形態に係る曲線座標生成装置は、前記各点で曲率を指定でき、かつ、前記各点間の曲率が単調となる曲線生成法を用いて、前記出力部により出力された座標の点列を通過する曲線を生成する曲線生成部を更に具備する。
【００１４】
本発明の一形態に係る曲線座標生成装置は、前記出力部により出力された座標の点列を表示する表示部と、前記表示部に表示された座標の点列における各点の曲率及び弧長を前記表示部上で可変操作する操作部とを更に具備する。
【００１５】
本発明の一形態に係る曲線を有する物体を生産する方法は、離散点で表された前記曲線の各点の曲率、前記各点間の弧長を入力し、前記曲率に基づき円を構成し、前記弧長に基づき前記円上で座標を特定し、前記特定された座標に基づき前記物体の曲線を設計し、前記物体を生産する。
【００１６】
本発明の一形態に係るプログラムは、離散点で表された曲線の各点の曲率、前記各点間の弧長を入力するステップと、前記曲率に基づき円を構成し、前記弧長に基づき前記円上で座標を特定するステップと、前記特定された座標を出力するステップとをコンピュータシステムに実行させる。
【発明の効果】
【００１７】
本発明によれば、離散点で表された曲線の各点の曲率及び各点間の弧長を入力とし、それに対応する座標を簡単な処理で出力することができる。
【図面の簡単な説明】
【００１８】
【図１】本発明の一実施形態に係る曲線座標生成装置の構成を示すブロック図である。
【図２】本発明に係る原理を説明するための図である。
【図３】本発明の一実施形態に係る曲線座標生成装置による具体的な処理の流れを説明するための図である。
【図４】座標算出部による詳細なプロセスを数式により説明するための図である。
【図５】曲率の正負についての説明図である。
【図６】円曲線とその点列を示す図である。
【図７】クロソイド曲線とその点列を示す図である。
【図８】曲率に基づいた点列修正機能に係る実施形態を説明するための図である。
【図９】一様流に置かれた翼と、その周辺の流線と圧力分布を示す図である。
【図１０】流線上を移動する流体粒子と遠心力と圧力差との関係を示す図である。
【図１１】翼の曲率から揚力・抗力を求めるフロー図である。
【図１２】本発明に係る曲線座標生成装置をインタラクティブに座標を求めるツールとして用いる場合の処理フロー図である。
【図１３】図１２に示した実施形態に係る画面の表示例を示す図である。
【図１４】ＣＦＤ用の格子点列の一例を示す図である。
【発明を実施するための形態】
【００１９】
以下、図面を参照しながら、本発明の実施形態を説明する。
【００２０】
図１は本発明の実施形態に係る曲線座標生成装置の構成を示すブロック図である。
【００２１】
図１に示すように、本発明の一形態に係る曲線座標生成装置１０は、入力部１１と、座標算出部１２と、出力部１３と、曲線生成部１４とを具備する。この曲線座標生成装置１０は、典型的には、コンピュータシステムに本発明に係るプログラムを実行させることで実現することができる。
【００２２】
曲線座標生成装置１０は、図２に示すように、二次元直交座標系の点列Ｐ（１）、Ｐ（２）、・・・Ｐ（ｎ）のそれぞれに対して弧長（点列の間隔）と曲率が既知であり、かつ、ｘｙ座標が未知である場合、隣接する３点が同一の円上に位置すると仮定し、幾何学的な手法（三平方の定理及び余弦定理）を用いることによって、それぞれの点列のｘｙ座標を算出する。ただし、Ｐ（１）とＰ（２）のｘｙ座標は初期値として与える必要がある。
【００２３】
図２の例では、Ｐ（２）、Ｐ（３）、Ｐ（４）が同一円Ｃ上に位置するという仮定を用い、加えてＰ（２）とＰ（３）のｘｙ座標、Ｐ（３）Ｐ（４）の弧長Ｓ３、Ｐ（３）の曲率Ｒ３の情報を使って、Ｐ（４）のｘｙ座標を求めている。同様の手順を踏むことで、点列Ｐ（５）、Ｐ（６）、・・・Ｐ（ｎ）のｘｙ座標を求めることができる。
【００２４】
ここで、入力部１１は、離散点で表された曲線の曲率Ｒｎ、離散点間の弧長Ｓｎ、離散点の第１の座標（ｘｎ−１，ｙｎ−１）及び第１の座標の次の第２の座標（ｘｎ，ｙｎ）を入力する。
【００２５】
座標算出部１２は、第１の座標（ｘｎ−１，ｙｎ−１）、第２の座標（ｘｎ，ｙｎ）及び第２の座標（ｘｎ，ｙｎ）における曲率Ｒｎに基づき円Ｃを構成し、第２の座標（ｘｎ，ｙｎ）から、第２の座標（ｘｎ，ｙｎ）と第２の座標の次の第３の座標（ｘｎ＋１，ｙｎ＋１）となる点との間の弧長Ｓｎだけ離れた円Ｃ上の座標を第３の座標（ｘｎ＋１，ｙｎ＋１）として算出する。
【００２６】
出力部１３は、座標算出部１２で算出された第３の座標（ｘｎ＋１，ｙｎ＋１）を出力する。
【００２７】
図３はこの曲線座標生成装置１０による具体的な処理の流れを説明するための図である。  
入力部１１は、以下の値を入力する
（１）離散点における各点の曲率Ｒ
（２）離散点における各点間の弧長Ｓ
（３）初期条件として始点の座標（ｘ１，ｙ１）及び次の点の座標（ｘ２，ｙ２）
【００２８】
座標算出部１２は、以下の処理を実行する
（１）初期条件としての始点の座標（ｘ１，ｙ１）及び次の点の座標（ｘ２，ｙ２）をｘｙ座標上の既知の点丸１及び点丸２として用意する（図３（ａ））
（２）点丸２の曲率から点丸１及び点丸２のを通過する円Ｃを三平方の定理を用いて描く（図３（ｂ））。この円Ｃ上では第丸３の点は一意に求められていない
（３）点丸２と点丸３との間の弧長から余弦定理を用いて点丸２から弧長だけ離れた円Ｃ上の点丸３の座標（ｘ３，ｙ３）を決定する（図３（ｃ））。
【００２９】
座標算出部１２は、座標（ｘ２，ｙ２）及び次の点の座標（ｘ３，ｙ３）をｘｙ座標上の既知の点として用意し、上記の（２）と（３）の処理を実行し、以下この処理を繰り返し、離散点で表された曲線の各点におけるｘｙ座標を算出する。
【００３０】
出力部１３は、算出された各点のｘｙ座標を出力する。
曲線生成部１４は、各座標で曲率を指定でき、かつ、隣接する座標間の曲率が単調となる曲線生成法を用いて、出力部１３により出力されたｘｙ座標の点列を通過する曲線を生成して出力する。
【００３１】
次に、上記座標算出部１２による詳細なプロセスを数式により説明する。図４はそのプロセスを説明するための図である。
・定義
ｘ：ｘ座標
ｙ：ｙ座標
ｓ：弧長
ｋ：曲率
ｒ：曲率半径（＝１／ｋ）
【００３２】
・入力
ｘｉ−１，ｙｉ−１：Ｐi-1の座標
ｘｉ，ｙｉ：Ｐiの座標
ｓｉ：ＰiとＰi＋1との間の弧長
ｋｉ：Ｐiの曲率（図５に示すように、曲率ｋｉが正は下に凸、負は上に凸となる。）
ｒｉ：Ｐiの曲率半径（ｋi＝１／ｒi）
【００３３】
・出力
ｘｉ＋１，ｙｉ＋１：Ｐi＋1の座標
・変数
ｘｍ，ｙｍ：Ｐi−１とＰiとの間の中点Ｐｍの座標
ｘｃ，ｙｃ：Ｐi−１とＰiを通過し半径１／ｋiである円の中心Ｐｃの座標
ｘｍ２ｃ，ｙｍ２ｃ：ＰｍとＰｃに向かう単位ベクトルＶｍ２ｃの終点の座標
θ：三角形ＰiＰi＋1Ｐｃの中心Ｐｉにおける内角
【００３４】
・関数
ａｂｓ：絶対値
ｓｉｇｎ（ａ，ｂ）：ｂの符号にａを乗じた数値
ａｃｏｓ：ａｒｃｃｏｓ
座標算出部１２は以下のアルゴリズムを実行する
【００３５】
【数１】
2018151774.tif 000003
【００３６】
ここで、上記のアルゴリズムにおいて、
「ｉｆ  ａｂｓ（ｋi）＞0  ｔｈａｎ」は、曲率が0でない場合であり、「ｅｌｓｅ」は曲率が0の場合である。
【００３７】
・円曲線の点列
円曲線
ｘ＝ｃｏｓθ，ｙ＝ｓｉｎθ
について上記アルゴリズムにより点列を求めると以下に示す座標値が得られた。
【００３８】
【表１】
2018151774.tif 000004
【００３９】
図６は円曲線とその点列を示す図であり、点列は厳密解と一致することが確認された。
・クロソイド曲線の点列
クロソイド曲線
【００４０】
【数２】
2018151774.tif 000005
【００４１】
について上記アルゴリズムにより点列を求めると以下の座標値が得られた。
【００４２】
【表２】
2018151774.tif 000006
【００４３】
図７はクロソイド曲線とその点列を示す図であり、点列は厳密解と一致することが確認された。
【００４４】
次に、本発明の別の実施形態を説明する。
これまでに説明した実施形態では、始点（ｘ１,ｙ１）は指定できたが終点座標（ｘｎ,ｙｎ）を指定できないという課題がある。
これに対して、非特許文献１では、終点座標を所望の座標とするために、その「Ｆｉｇｕｒｅ  １．（ｂ）」に示されるように、点列全体に所望の終点座標に基づいた修正量を加えているが、修正後の曲率分布はその意図の有無にかかわらず入力した曲率分布と異なることになる。また、文献（Lin, W. F.: 3D TRANSONIC NACELLE AND WINGLET DESIGN, AIAA 90-3064, 1990）では、その「Ｆｉｇｕｒｅ  ６．」及び「Ｆｉｇｕｒｅ  ７．」に示されるように、局所的に角度に基づいた修正量を加えているが、上記の非特許文献１と同様に、修正後の曲率分布はその意図の有無にかかわらず入力した曲率分布と異なることになる。
【００４５】
そこで、この実施形態に係る曲線座標生成装置１０は、終点座標を所望の座標に一致させるために、以下に示す曲率に基づいた点列の修正及び曲率修正量の重みの入力による修正後の曲率分布の制御を実行する。以下、図８を参照しながら説明する。
【００４６】
（入力部１１）
・弧長ｓｉと曲率ｋｉに加えて、曲率修正量に関する重みｗｉを入力する。
・始点座標（ｘ１，ｙ１）及び次の座標（ｘ２，ｙ２）に加えて、終点座標（ｘｎ，ｙｎ）を入力する。
【００４７】
（座標算出部１２）
・次の関数ｆ（ａ）を作成する。
比例定数ａを入力とする。
弧長ｓｉと曲率ｋｉ（１＋ａｗｉ）とする。
［数１］に示したアルゴリズムを用いて点列座標を計算する。
計算によって得られた終点座標（ｘｔｎ，ｙｔｎ）と入力した座標（ｘｎ，ｙｎ）の距離を関数ｆ（ａ）の出力とする。
【００４８】
・ｆ（ａ）＝０となる比例定数ａ０を求める。例えば繰り返し計算で比例定数ａ０の収束値を求める。
・弧長ｓｉと曲率ｋｉ（１＋ａ０ｗｉ）とする。
・［数１］に示したアルゴリズムを用いて点列座標（ｘｉ，ｙｉ）を計算する。
・曲率ｋｉ（１＋ａ０ｗｉ）及び座標（ｘｉ，ｙｉ）を出力する。
【００４９】
従って、この実施形態に係る曲線座標生成装置１０では、曲率ｋｉに基づいた点列座標（ｘｉ，ｙｉ）の修正が可能であると共に、曲率修正量の重みｗｉの入力によって修正後の各点列（ｘｉ，ｙｉ）の曲率分布ｋｉ（１＋ａ０ｗｉ）の制御が可能である。
【００５０】
以上説明した本発明は様々な分野において有効である。典型的には、流体力学の分野（例えば空力性能向上のため形状曲率制御）、数値流体力学（例えばＣＦＤ用計算格子の点制御）、意匠設計（例えば意匠曲線(点列)を描くための曲率制御）、経路設計（例えば経路のコーナ部分曲率制御）などの分野において本発明は有効である。
【００５１】
・空力性能向上のため形状曲率制御に本発明を適用した場合
翼の性能を表す揚力・抗力は圧力分布を翼形状周りで積分することで得られる。図９に示すように、一様流に置かれた翼８１は、その存在によって周辺の流線８２と圧力分布８３を変化させる。
【００５２】
ここで、ある一本の流線８２に着目したとき、翼８１が存在することによって流線８２は曲げられるが、図１０に示すように、流線８２上を移動する流体粒子８４においては遠心力８５が発生すると同時にそれに釣り合う圧力差８６が発生することが空気力学の流線曲率定理から知られている。流線曲率定理は、例えば「流体力学（前編）」（今井功著  裳華房発行）に詳しく記載されている。
【００５３】
また、流れが翼８１に付着している限り、流線８２の曲率は翼８１の曲率と対応することが空気力学の微少擾乱理論から知られている。微少擾乱理論は、例えば「気体力学」（Liepmann, H.W., Roshko,著  吉岡書店発行）に詳しく記載されている。
【００５４】
図１１は以上の説明をまとめたものである。
【００５５】
ここで、翼の性能を向上させるための設計手法を考える。
これまでの設計手法は殆どが"翼の座標"に基づいた手法であったように思われる。しかしながら空気力学理論を考えたとき、圧力分布と"翼の座標"の関係性よりも、圧力分布と"翼の曲率"の関係性の方が密接と考えられる。つまり"翼の曲率"に基づいた設計手法の方が空気力学的には自然であると考えられる。
【００５６】
このように曲率の重要性を示す定理が存在しているにもかかわらず、"翼の曲率"に基づいた設計手法がこれまでに少なかったか理由を考えると、その理由は、曲率制御を可能とする曲線や点列の生成技術（非特許文献１)が周知されてなかったためと考えられる。
ここで、本発明と非特許文献１に記載された技術とを比較すると、両者ともに曲率制御を可能とする点列の生成技術であり、"翼の曲率"に基づいた設計手法を可能とするための基礎技術といえる。しかしながら、非特許文献１では弧長による制御ができず、またｘ座標に対してｙ座標が多価関数の場合には使えず、普遍的ではない。これに対して、本発明は、離散点で表された曲線の各点の曲率及び各点間の弧長を入力とし、それに対応する座標を簡単な処理で出力することができ、またｘ座標に対してｙ座標が多価関数の場合には使えるので、非特許文献１に記載された技術と比べてより普遍的な手法である。
【００５７】
・意匠設計に本発明を適用した場合
意匠面とは、製品の表面形状を構成する曲面であり、意匠的な意図を持って形状設計された曲面を指す用語である。意匠面形状の設計においては「光線の映り込み」、「曲がり具合の味わい」など、その意匠的な意図に曲面の曲率が深く関わる。この点については、例えば「清水保弘，意匠面創生のための曲率単調曲線生成技術，UNISYS TECHNOLOGY REVIEW vol.71, 2001.」に詳しく記載されている。
【００５８】
曲率を制御することは「光線の映り込み」、「曲がり具合の味わい」を制御することに繋がり、意匠曲線を描く際の重要な課題の一つといえる。そのため、曲率制御を目的とした曲線生成技術に関する研究は散見されるが、いまだ普遍的な手法が確立されていないものと考える。一方、本発明は曲率制御を可能とする"点列"生成技術であるが、具体例で示したように普遍性の高い技術であると考える。
【００５９】
本発明を用いて作成された"点列"を基にして"曲線"を描くための手法は、典型的には、以下のものがある。なお、本発明を用いることでそれぞれの点における座標、曲率が既知となるが、この時点で正接（ｔａｎｇｅｎｔ）の情報は欠けている。ただし、必要であれば両隣の点を用いて離散的に正接を求めることは可能である。
【００６０】
米国特許公開公報２００４−１０９００１号は、それぞれの点における座標、正接、曲率が既知の場合に曲率単調な曲線を描くための補間技術を開示する。そして、補間技術と本発明を組み合わせることで、所望の座標を通過し、かつ所望の正接、曲率を有する曲線を描くことができる。もちろん、米国特許公開公報２００４−１０９００１号は、本明細書の開示の範囲に含まれる。
【００６１】
もう一つの手法は、本発明から出力される点列の密度を十分に高くし、その点列の座標を通過するＮＵＲＢＳ（Ｎｏｎ−Ｕｎｉｆｏｒｍ−Ｒａｔｉｏｎａｌ−Ｂ−Ｓｐｌｉｎｅ)曲線を描くことである。実際にＣＡＤソフトウェアの「ＣＡＴＩＡ」（Ｄａｓｓｕｌｔ社,  ｈｔｔｐ：／／ｗｗｗ．３ｄｓ．ｃｏｍ／ＪＡ／ｐｒｏｄｕｃｔｓ−ｓｅｒｖｉｃｅｓ／ｃａｔｉａ／参照）上でＮＵＲＢＳ曲線の曲率と点列の離散曲率がほぼ一致することを確認している。後者の手法は前者の手法に比べて簡易的ではあるが、実用的であると考える。
【００６２】
・ＣＡＤに本発明を適用した場合
本発明は、ＣＡＤ（Ｃｏｍｐｕｔｅｒ  Ａｉｄｅｄ  Ｄｅｓｉｇｎ）上で点列の離散曲率と弧長を入力して、インタラクティブに座標を求めるツールとして用いることが可能である。
【００６３】
図１２は本発明に係る曲線座標生成装置をインタラクティブに座標を求めるツールとして用いる場合の処理フローの一例を示している。
コンピュータシステム上で構成された曲線座標生成装置は、例えば本発明に係る処理により生成された点列を取り込む（ステップ１１１）と、点列の条件判定を行い（ステップ１１２）、点列に基づき曲線分布を画面上に表示する（ステップ１１３）。
【００６４】
図１３はその画面表示の一例を示している。図１３に示す画面上の第１のウィンドウ１２１上には、点列に基づき曲線分布１２２が表示され、第２のウィンドウ１２３には、曲率とｘ座標との関係を示すグラフ１２４と、曲率を変更するためのスライダー１２５とが表示されている。なお、弧長についても同様にウィンドウに表示し、弧長を変更するためにスライダーを使うことができる。
【００６５】
曲線座標生成装置は、スライダー１２５などを使って曲率や弧長が変更されると（ステップ１１４）、取り込んだ点列の更新を行い（ステップ１１５）、更新された点列に基づき変更された曲線分布１２２を画面上に表示する（ステップ１１３）。この曲線座標生成装置に係るアプリケーションが終了されるまで、以上の処理を継続する（ステップ１１６）。
【００６６】
空力・流体設計で使用される数値流体力学（ＣＦＤ（Ｃｏｍｐｕｔｅｒ  Ｆｌｕｉｄ  Ｄｙｎａｍｉｃｓ））では、物体表面の格子点列間隔（弧長）（図１３参照）が計算精度に影響するが、上記のＣＡＤに本発明を適用した場合と同様の手法によってこの格子点列間隔（弧長）を制御することが可能となる。つまり、本発明に係る曲線座標生成装置を、ＣＦＤツールの入力ファイルである計算格子点列を出力する装置、更にＣＦＤツールの出力ファイルを読み込んで、その情報を基に弧長と曲率を変更する機能を有する装置として用いることができる。
【００６７】
以上説明したとおり、流体力学では、ベルヌーイの定理と並んで流線曲線の定理という基本的な定理が存在し、曲率と流れ場は密接な関係にある。加えて、数値流体力学では物体表面を離散点列で扱うことが多く、格子点列間隔（弧長）の設定は基本的な検討課題のひとつである。本発明は、流体力学と数値流体力学に関連する基本技術であり、両分野を扱う産業（航空宇宙、船舶、自動車等）で利用される可能性は高い。
【００６８】
意匠の分野では、Ａｅｓｔｈｅｔｉｃな曲線を描くために曲率制御は重要な検討課題のひとつである。本発明は、曲線ではなくて点列を対象としているが、その代わり、比較的容易な幾何学を用いた手法で曲率制御を可能とする。そのため、意匠を扱う産業（建築、自動車等）で利用される可能性が高い。
【００６９】
経路設計においても、経路のコーナーにおける曲率制御は重要な検討課題のひとつである。そのため、経路を扱う産業（インフラ、航空輸送等）で利用される可能性が高い。
【符号の説明】
【００７０】
１０  曲線座標生成装置
１１  入力部
１２  座標算出部
１３  出力部
１４  曲線生成部
１２１  第１のウィンドウ
１２３  第２のウィンドウ
１２５  スライダー
【図１】
2018151774.tif 000007
【図２】
2018151774.tif 000008
【図３】
2018151774.tif 000009
【図４】
2018151774.tif 000010
【図５】
2018151774.tif 000011
【図６】
2018151774.tif 000012
【図７】
2018151774.tif 000013
【図８】
2018151774.tif 000014
【図９】
2018151774.tif 000015
【図１０】
2018151774.tif 000016
【図１１】
2018151774.tif 000017
【図１２】
2018151774.tif 000018
【図１３】
2018151774.tif 000019
【図１４】
2018151774.tif 000020