(57)【要約】【課題】産業機械の形状を最適化するにあたり、探索の精度が悪化することを抑制しつつ、短時間で最適解の探索を完了する。【解決手段】最適化対象部の形状のうち、予め設定された要求性能以上の性能を有する形状を含む複数の形状を主成分分析することによって主成分を抽出する主成分抽出工程と、  主成分に基づいて最適化対象部の新規形状を生成する新規形状生成工程と、新規形状生成工程にて生成された新規形状の最適化対象部の性能を評価する評価工程とを有する。【選択図】図５
2013003971.tif 000002
【特許請求の範囲】
【請求項１】
産業機械の一部あるいは全部である最適化対象部の形状を最適化する形状最適化方法であって、
予め設定された要求性能以上の性能を有する形状を含む複数の形状を主成分分析することによって主成分を抽出する主成分抽出工程と、
前記主成分に基づいて前記最適化対象部の新規形状を生成する新規形状生成工程と
を有することを特徴とする形状最適化方法。
【請求項２】
前記主成分抽出工程にて、前記主成分分析で用いられる前記複数の形状は、全てが予め設定された要求性能以上の性能を有する形状であることを特徴とする請求項１記載の形状最適化方法。
【請求項３】
前記主成分分析で用いられる前記複数の形状を同一の座標系におけるベクトルに変換するベクトル変換工程を有し、
前記主成分抽出工程にて、変換後の前記ベクトルに基づいて前記主成分分析を行う
ことを特徴とする請求項１または２記載の形状最適化方法。
【請求項４】
前記主成分分析は、前記主成分を基底ベクトルとして抽出する固有直交分解であることを特徴とする請求項１〜３いずれかに記載の形状最適化方法。
【請求項５】
予め定められた数だけ低次側から前記基底ベクトルを選択する基底ベクトル選択工程を有し、
前記新規形状生成工程にて、前記基底ベクトル選択工程にて選択された前記基底ベクトルに基づいて前記最適化対象部の新規形状を生成する
ことを特徴とする請求項４記載の形状最適化方法。
【請求項６】
前記新規形状生成工程にて生成された前記新規形状の最適化対象部の性能を評価する評価工程を有することを特徴とする請求項１〜５いずれかに記載の形状最適化方法。
【請求項７】
産業機械の一部あるいは全部である最適化対象部の形状を最適化する形状最適化装置であって、
予め設定された要求性能以上の性能を有する形状を含む複数の形状を主成分分析することによって主成分を抽出する主成分抽出手段と、
前記主成分に基づいて前記最適化対象部の新規形状を生成する新規形状生成手段と
を備えることを特徴とする形状最適化装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、形状最適化方法及び形状最適化装置に関するものである。
【背景技術】
【０００２】
産業機械においては、流体に晒される翼や流体を用いる圧縮機等の流体に関係する様々な部品や装置が用いられている。
これらの部品や装置の形状（すなわち産業機械の形状）は、産業機械の性能を大きく左右する。このため、より高性能な産業機械を設計するために、産業機械の形状を最適化する試みが広く行われている。
【０００３】
流体の支配方程式は非線形性が強いため、線形的な法則に倣って形状を最適化することが難しい。このため、従来から流体に関係する部品や装置の形状最適化は、予め指定される条件を満足する仮想的な探索空間の中で発見的な手法を用いて行われている。
【０００４】
発見的な手法としては、例えば、非特許文献１に記載される遺伝的アルゴリズムを用いた手法、非特許文献２に記載される応答曲面法を用いた手法、非特許文献３に記載される逐次二次計画法を用いた局所最適化手法が行われている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】Deb, K.: Multi-objective optimization using evolutionary algorithm, John Wiley & Sons, New York, 2001.
【非特許文献２】Myers, R. H. and Montgomery, D. C.: Responce Surface Methodology: process and product optimization using designed experiments, second edition, A Willey-Interscience publication, New York, 2002.
【非特許文献３】勝井辰博，冨田高嗣: ＣＡＤ／ＣＦＤ  統合型高効率最適化システムの構築とその検証  −２次元翼型の最適化による最適化システムの検証−, 日本船舶海洋工学会論文集, 4, 175-183 (2006).
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、一般的に最低限の条件で規定される探索空間は非常に広い。探索空間が広い場合には、発見的な最適化手法を用いて探索空間の全領域を探索して最適解を見つけることは非常に時間がかかり、また探索の精度も悪化する。
【０００７】
なお、探索空間を絞る方法としては、予め基準となる形状を複数設定し、この複数の形状を組み合わせて探索を行うことが考えられる。
しかしながら、産業機械あるいはこの部品等は、様々な表現形式によって表されている。例えば、翼は、非特許文献３に示すベジェ曲線を用いた表現手法や、スプライン曲線を用いた表現手法で表されている。
このような異なる表現手法で表された形状同士を組み合わせることはできない。このため、解析可能な形状が限定されてしまい、探索の精度が悪化してしまう。
【０００８】
本発明は、上述する問題点に鑑みてなされたもので、産業機械の形状を最適化するにあたり、探索の精度が悪化することを抑制しつつ、短時間で最適解の探索を完了することを目的とする。
【課題を解決するための手段】
【０００９】
本発明は、上記課題を解決するための手段として、以下の構成を採用する。
【００１０】
第１の発明は、産業機械の一部あるいは全部である最適化対象部の形状を最適化する形状最適化方法であって、予め設定された要求性能以上の性能を有する形状を含む複数の形状を主成分分析することによって主成分を抽出する主成分抽出工程と、上記主成分に基づいて上記最適化対象部の新規形状を生成する新規形状生成工程とを有するという構成を採用する。
【００１１】
第２の発明は、上記第１の発明において、上記主成分抽出工程にて、上記主成分分析で用いられる上記複数の形状が、全てが予め設定された要求性能以上の性能を有する形状であるという構成を採用する。
【００１２】
第３の発明は、上記第１または第２の発明において、上記主成分分析で用いられる上記複数の形状を同一の座標系におけるベクトルに変換するベクトル変換工程を有し、上記主成分抽出工程にて、変換後の上記ベクトルに基づいて上記主成分分析を行うという構成を採用する。
【００１３】
第４の発明は、上記第１〜第３いずれかの発明において、上記主成分分析が、上記主成分を基底ベクトルとして抽出する固有直交分解であるという構成を採用する。
【００１４】
第５の発明は、上記第４の発明において、予め定められた数だけ低次側から上記基底ベクトルを選択する基底ベクトル選択工程を有し、上記新規形状生成工程にて、上記基底ベクトル選択工程にて選択された上記基底ベクトルに基づいて上記最適化対象部の新規形状を生成するという構成を採用する。
【００１５】
第６の発明は、上記第１〜第５いずれかの発明において、上記新規形状生成工程にて生成された上記新規形状の最適化対象部の性能を評価する評価工程を有するという構成を採用する。
【００１６】
第７の発明は、産業機械の一部あるいは全部である最適化対象部の形状を最適化する形状最適化装置であって、最適化対象部の形状のうち、予め設定された要求性能以上の性能を有する形状を含む複数の形状を主成分分析することによって主成分を抽出する主成分抽出手段と、上記主成分に基づいて上記最適化対象部の新規形状を生成する新規形状生成手段と、上記新規形状生成工程にて生成された上記新規形状の最適化対象部の性能を評価する評価手段とを備えるという構成を採用する。
【発明の効果】
【００１７】
本発明によれば、予め設定された要求性能以上の性能を有する形状を含む複数形状から主成分を抽出する。
このように抽出された主成分には、高性能を引き出す特徴が多く含まれていると考えられる。つまり、この主成分に基づいて作成する新規形状には、より高性能なものが含まれる可能性が高い。
そして、本発明によれば、より高性能なものが含まれる可能性が高い新規形状を中心として評価を繰り返すことにより、最適解を探索することとなる。したがって、本発明によれば、産業機械の形状を最適化するにあたり、探索の精度が悪化することを抑制しつつ、短時間で最適解の探索を完了することが可能となる。
【図面の簡単な説明】
【００１８】
【図１】本発明の一実施形態における形状最適化方法及び形状最適化装置にて用いられるＰＯＤの概略を説明する説明図である。
【図２】本発明の一実施形態における形状最適化方法及び形状最適化装置にて用いられるＰＯＤにより得られる設計空間を示す概念図である。
【図３】本発明の一実施形態における形状最適化方法の概略を示すフローチャートである。
【図４】本発明の一実施形態における形状最適化装置のハードウェア構成を示すブロック図である。
【図５】本発明の第１実施形態における形状最適化方法及び形状最適化装置を用いた翼形状最適化を説明するためのフローチャートである。
【図６】本発明の第１実施形態における形状最適化方法及び形状最適化装置で用いる翼形状を示す模式図である。
【図７】本発明の第１実施形態における形状最適化方法及び形状最適化装置で得られる固有ベクトルと基底ベクトルとによる形状を示すグラフである。
【図８】本発明の第１実施形態における形状最適化方法及び形状最適化装置で得られる基底ベクトルのノルムを示すグラフである。
【図９】本発明の第１実施形態における形状最適化方法及び形状最適化装置で得られる揚力係数ＣＬが大きく抗力係数ＣＤとの関係を示すグラフである。
【図１０】本発明の第１実施形態における形状最適化方法及び形状最適化装置によって選択された最適解に基づく形状と、他の従来手法にて最適化された最適解に基づく形状とを示すグラフである。
【図１１】本発明の第２実施形態における形状最適化方法及び形状最適化装置で用いる流路モデルを示す模式図である。
【図１２】本発明の第２実施形態における形状最適化方法及び形状最適化装置で得られる基底ベクトルのノルムを示すグラフである。
【発明を実施するための形態】
【００１９】
以下、図面を参照して、本発明に係る産業機械の形状最適化方法及び産業機械の形状最適化装置の一実施形態について説明する。
まず具体的な実施形態の説明の前に本発明の原理について説明する。
【００２０】
一般的な産業機械の歴史は古く、現在にたどり着くまでの改良によって、性能が高いとされる形状（以下、高性能形状）が複数発見されている。そして、このような高性能形状には、高い性能を引き出す要因となる特徴的な形状（以下、特徴形状）を含んでいると考えられる。
このような特徴形状を複数の高性能形状から抽出して組み合わせ、別の形状を生成した場合、この形状には多くの特徴形状が含まれることとなり、より高性能な形状となる可能性が高い。
つまり、探索空間において、高性能形状から抽出した特徴形状を含む領域に最適解が存在する可能性が高いと考えられる。
したがって、探索空間における探索領域を特徴形状が含まれる領域に限定することによって、探索領域が大幅に限定され、短時間で最適解を導き出すことができると考えられる。
【００２１】
そして、本発明においては、高性能形状から特徴形状を抽出する手法として、主成分分析を用いる。
この主成分分析は、複数の高性能形状をベクトルとして表現して同一の座標空間に配置し、この座標空間に配置された複数の高性能形状を示す座標の配置傾向から特徴形状を抽出する手法である。
このような主成分分析の手法は、固有直交分解（ＰＯＤ：Proper Orthogonal Decomposition）、カルーネン・レーベ展開（Karhunen-Loeve decomposition）、特異値分解（Singular value decomposition）、スペクトル分解、独立成分分析（Independent component analysis）、因子分析として用いられている。
【００２２】
ここで、固有直交分解（以下、ＰＯＤと称する）を用いて最適解の探索手法（以下、本探索方法と称する）について説明する。
まずＰＯＤについて説明する。なお、以下の説明において、ベクトルは下式（１）などの太字で表し、転置ベクトルを下式（２）で表すものとする。また、ベクトルのノルムは全て下式（３）で示すユークリッドノルムとする。また、扱う値は全て実数である。
【００２３】
【数１】
2013003971.tif 000003
【００２４】
【数２】
2013003971.tif 000004
【００２５】
【数３】
2013003971.tif 000005
【００２６】
ＰＯＤは、主成分分析の一手法であり、相関があることが期待される複数のベクトルを、基底ベクトル（主成分）に分解する方法である。なお、基底ベクトルは、特徴形状を示す。
【００２７】
図１は、上述したＰＯＤの概要を纏めた説明図である。
例えば、予めｎ個の形状の幾何学情報が下式（４）の形で与えられているものとする。このとき下式（５）のｊ番目の値を下式（６）で表すと約束すると、式（６）は、各形状に共通するｘ座標ｘｊに対するｙ座標値である。相関を求めるために、与えられた式（５）を平均成分とそれ以外の成分に下式（７）と下式（８）に示すように分解する。
【００２８】
【数４】
2013003971.tif 000006
【００２９】
【数５】
2013003971.tif 000007
【００３０】
【数６】
2013003971.tif 000008
【００３１】
【数７】
2013003971.tif 000009
【００３２】
【数８】
2013003971.tif 000010
【００３３】
ＰＯＤでは、下式（９）に示すように、下式（１０）を下式（１１）の線形和に分解する。
【００３４】
【数９】
2013003971.tif 000011
【００３５】
【数１０】
2013003971.tif 000012
【００３６】
【数１１】
2013003971.tif 000013
【００３７】
ただし、式（１１）は直交基底である。また、下式（１２）とすると、下式（１３）を満たす。
【００３８】
【数１２】
2013003971.tif 000014
【００３９】
【数１３】
2013003971.tif 000015
【００４０】
そして、ＰＯＤの目的は、式（９）及び式（１３）を満たし、添え字ｊが小さい順に下式（１４）が最大となるような式（１１）と下式（１５）を求めることである。
【００４１】
【数１４】
2013003971.tif 000016
【００４２】
【数１５】
2013003971.tif 000017
【００４３】
ここで、行列Ｖを下式（１６）と定義すると、式（１５）及び式（１１）は、式（１７）の固有値分解を用いて求めることができる。式（９）を用いると、式（１７）は、下式（１８）のように記述できる。
【００４４】
【数１６】
2013003971.tif 000018
【００４５】
【数１７】
2013003971.tif 000019
【００４６】
【数１８】
2013003971.tif 000020
【００４７】
ここで両辺に左から式（１９）、右から式（２０）をかけると、下式（２１）を得る。
【００４８】
【数１９】
2013003971.tif 000021
【００４９】
【数２０】
2013003971.tif 000022
【００５０】
【数２１】
2013003971.tif 000023
【００５１】
したがって、下式（２２）を固定したとき、下式（２３）を最大化する下式（２４）は式（２２）と直交する。さらにこのとき式（２３）を最大化する式（２２）は、式（１７）の最大固有値に対向する正規固有ベクトルである。同様にして、式（２４）は全て式（１７）の最大固有値に対応する正規固有ベクトルである。同様にして式（２４）は全て式（１７）のｊ番目に大きい固有値に対応する正規固有ベクトルであることが従う。最後に式（２４）が直交しているから、下式（２５）が従う。
【００５２】
【数２２】
2013003971.tif 000024
【００５３】
【数２３】
2013003971.tif 000025
【００５４】
【数２４】
2013003971.tif 000026
【００５５】
【数２５】
2013003971.tif 000027
【００５６】
以上のように行列の固有値分解を利用してＰＯＤを行うことができる。なお、数値実験においては、行列のサイズが比較的小さいこともあり、固有値分解の方法として実装が簡単なべき乗法を用いた。
【００５７】
そして、このようなＰＯＤによって得られた直交分解から新しい形状を生成することができる。
ここで、基底ベクトルである式（１１）は、ｎ個得られるが、ｊが大きくなるにつれてそのノルムは小さくなる。特に最初に与えた形状の個数が多く、さらに相関が大きければ高次の基底ベクトルのノルムは無視できるほど小さくなることが期待される。
このため、例えば基底ベクトルとして、そのノルムが閾値以上であるノルムのみを形状最適化のために用いることが考えられる。
具体的には、Ｒ＝０．０１に対して、下式（２６）などとして、下式（２７）であるような基底ベクトルを用いて新たな形状を生成する。
このとき、新たな形状を表す下式（２８）は、下式（２９）で表される。
【００５８】
【数２６】
2013003971.tif 000028
【００５９】
【数２７】
2013003971.tif 000029
【００６０】
【数２８】
2013003971.tif 000030
【００６１】
【数２９】
2013003971.tif 000031
【００６２】
ただし、下式（３０）は任意に設定するパラメータであり、与え方の自由度は大きい。今回は、下式（３１）の範囲で適当な格子点を選択した。ただし、制約条件が課されている場合にはその制約条件を満たす範囲で式（３２）を選択する。
【００６３】
【数３０】
2013003971.tif 000032
【００６４】
【数３１】
2013003971.tif 000033
【００６５】
【数３２】
2013003971.tif 000034
【００６６】
式（２９）により超平面が与えられる。これに拘束条件を考慮に入れることで超平面上の部分空間が設計空間（探索空間）となる。この設計空間は、予め与えられた高性能形状等の形状を示す形状ベクトルで張られる空間の部分空間でもある。
この設計空間と予め与えられた形状との関係は、図２のようになる。ここでは簡略化するために２次元平面上に描いているが、実際には次元はもっと大きい。
【００６７】
式（２６）におけるＲを適切に選ぶことでＰＯＤにより得られる設計空間は図２のように最初に与える幾何情報の点集合の十分近くを通る平面にできる。仮に元の幾何情報の間の相関が小さい場合であってもＲ＝０とすることで式（２９）の超平面が元の幾何情報ベクトル（形状ベクトル）をすべて含むようにできる。
ここでＲを変化させるとＮが変化し、ＰＯＤにより得られる設計空間の次元が変化することになる。
【００６８】
上述のように、超平面上の目的関数値の集合（高性能形状が有する特徴的な形状を含む形状の集合）は比較的良い値を含んでいる可能性が高い。したがって、超平面上の目的関数値の集合を探索するのみで最適解が得られると考えられる。
つまり、ＰＯＤを用いない場合には、図２の「feasible region」のような広い設計空間全体を探索しなければならないのに対して、ＰＯＤを用いた場合には図２の「POD-design-space」を探索すれば十分であると考えられる。
すなわち、ＰＯＤにより高性能な解の集合を抽出しようとしていると考えればよい。
特に目的関数が連続でかつ、Ｒを正として打ち切った基底ベクトルの影響が無視できるならば、超平面上の目的関数値は、少なくとも元形状に相当する目的関数値の最大値と最小値に挟まれた区間を含む。このとき、この超平面上を探索すれば元形状と同じかそれより高性能な解が求まる。
また、設計空間が狭められたことで応答曲面を生成する場合に曲面の精度が向上して、最適解の精度が向上し探索時間が短縮することが期待できる。
【００６９】
上述のような主成分分析を用いて最適解を探索する場合には、図３に示すように、最初に事前にいくつかの性能の良い形状を取得する（ステップＳ１１）。これは経験や知見に従って得るか、あるいはランダムサンプリングにより得ることができる。
次に、それらの形状についてＰＯＤを実行する（ステップＳ１２）。続いて、ＰＯＤにより得られた基底ベクトルのうち、高次の基底ベクトルを削除する（ステップＳ１３）。そして、低次の基底ベクトルだけを用いてパラメトリックに新しい形状を生成する（ステップＳ１４）。そしてパラメータを変数としてＲＳＭ等を用いて最適解を探索する（ステップＳ１５）。
【００７０】
（第１実施形態）
次に、本発明の第１実施形態として、２次元翼の形状最適化に本発明を適用した例について説明する。つまり、本実施形態は、産業機械の一部である翼を最適化対象部とするものである。
なお、翼形状設計において重要な因子は揚力係数ＣＬと抗力係数ＣＤである。そこで、本実施形態においては、表１に示すような多目的最適化問題を考え、下式（３３）が最大となる形状を最適解とする。
【００７１】
【表１】
2013003971.tif 000035
【００７２】
【数３３】
2013003971.tif 000036
【００７３】
また、一般に翼厚比が大きいほど揚力係数ＣＬが大きく抗力係数ＣＤが小さくなる傾向がある。ＰＯＤを用いた場合、式（２９）において式（３２）を大きくすると翼厚比が大きくなり、それだけで揚力係数ＣＬが大きくなることがありえる。すなわち、新規に高性能な形状を生成できたとしても、それがＰＯＤの効果であるのか翼厚比を大きくしたことによる効果であるのか判断が困難である。そこで、本実施形態においては、すべての翼について翼厚比が一定の値になるよう縮尺を変える。
【００７４】
図４は、本実施形態において形状最適化方法を行う形状最適化装置１０のハードウェア構成を概略的に示すブロック図である。
この図に示すように、本実施形態の形状最適化装置１０は、パーソナルコンピュータやワークステーション等のコンピュータによって構成されるものであり、ＣＰＵ１、記憶装置２、リムーバブルメディアドライブ３、入力装置４、出力装置５、及び通信装置６を備えている。
【００７５】
ＣＰＵ１は、記憶装置２、リムーバブルメディアドライブ３、入力装置４、出力装置５、及び通信装置６と電気的に接続されており、これらの各種装置から入力される信号を処理すると共に、処理結果を出力するものである。
そして、ＣＰＵ１は、後述の形状最適化プログラムＰに基づいて、形状最適化装置１０を、高性能形状を含む複数の翼形状を主成分分析することによって主成分を抽出する主成分抽出手段として機能させる。
また、ＣＰＵ１は、形状最適化プログラムＰに基づいて、形状最適化装置１０を、主成分に基づいて翼の新規形状を生成する新規形状生成手段として機能させる。
また、ＣＰＵ１は、形状最適化プログラムＰに基づいて、形状最適化装置１０を、新規形状生成工程にて生成された新規形状の最適化対象部の性能を評価する評価手段として機能させる。
【００７６】
記憶装置２は、メモリ等の内部記憶装置及びハードディスクドライブ等の外部記憶装置によって構成されており、ＣＰＵ１から入力される情報を記憶すると共にＣＰＵ１から入力される指令に基づいて記憶した情報を出力するものである。
そして、本実施形態において記憶装置２は、プログラム記憶部２ａとデータ記憶部２ｂとを備えている。
【００７７】
プログラム記憶部２ａは、形状最適化プログラムＰを記憶している。この形状最適化プログラムＰは、所定のＯＳにおいて実行されるアプリケーションプログラムであり、コンピュータから構成される本実施形態の形状最適化装置１０に形状最適化方法を実行させるものである。
【００７８】
データ記憶部２ｂは、翼形状を示す形状データＤ１を記憶する。なお、データ記憶部２ｂには、高性能形状を示す形状データを含む複数の形状データが記憶される。
また、データ記憶部２ｂは、目的関数データＤ２、制約条件データＤ３、計算結果データＤ４等も記憶する。さらに、データ記憶部２ｂは、ＣＰＵ１の処理過程において生成される中間データを一時的に記憶する。
【００７９】
リムーバブルメディアドライブ３は、ＤＶＤ（Digital Versatile Disc）やＵＳＢ（Universal Serial Bus）メモリ等のリムーバブルメディアＸを取り込み可能あるいは接続可能に構成されており、ＣＰＵ１から入力される指令に基づいて、リムーバブルメディアメディアＸに記憶されるデータを出力するものである。そして、本実施形態においては、リムーバブルメディアメディアＸに形状最適化プログラムＰが記憶されており、リムーバブルメディアドライブ３は、ＣＰＵ１から入力される指令に基づいて、リムーバブルメディアメディアＸに記憶される数値解析プログラムＰを出力する。
なお、形状最適化プログラムＰは、必ずしもリムーバルディスクＸに格納されている必要はない。例えば、ネットワークＮを介して形状最適化プログラムＰを取得し、この取得した形状最適化プログラムＰを記憶装置２に記憶させるようにしても良い。
【００８０】
入力装置４は、本実施形態の形状最適化装置１０と作業者とのマンマシンインターフェイスであり、ポインティングデバイスであるキーボード４ａやマウス４ｂを備えている。
出力装置５は、ＣＰＵ１から入力される信号を可視化して出力するものであり、ディスプレイ５ａ及びプリンタ５ｂを備えている。
【００８１】
通信装置６は、本実施形態の形状最適化装置１０とネットワークＮとを電気的に接続し、ネットワークＮに接続された他の機器との間においてデータの受け渡しを行うものである。
なお、ネットワークＮとしては、社内ＬＡＮ（Local Area Network）やインターネット等が考えられる。
【００８２】
次に、このように構成された本実施形態の形状最適化装置１０を用いた形状最適化方法について、図５のフローチャートを参照しながら説明する。
なお、以下の説明においては、プログラム記憶部２ａに既に形状最適化プログラムＰが格納されているものとする。
【００８３】
本実施形態の形状最適化方法では、図５に示すように、まず目的関数の設定及び制約条件の設定が行われる（ステップＳ２１）。
具体的には、作業者が入力装置４を操作し、表１に示す目的関数及び制約条件を入力し、この入力データが目的関数データＤ２と制約条件データＤ３としてデータ記憶部２ｂに記憶されることで目的関数の設定及び制約条件の設定が完了する。
【００８４】
続いて、基礎形状の設定が行われる（ステップＳ２２）。
ここでは、まず作業者が、性能が高いとされる既知の形状（高性能形状）を複数選択する。なお、ここでは、高性能形状がどのような表現手法で表現されていても構わない。つまり、複数の高性能形状のうち、いくつかがベジェ曲線で表現され、いくつかがスプライン曲線で表現され、いくつかがフリーハンドで表現されていても良い。
そして、作業者が入力装置４を操作し、選択した高性能形状を入力し、この入力データが形状データＤ１としてデータ記憶部２ｂに記憶されることで基礎形状の設定が完了する。
【００８５】
なお、本実施形態においては、「Abbot, I. H. and von Doenhoff, A. E.: Theory of wing sections including a summary of airfoil data, Dover publications, INC., New York, 1959.」を参考に、揚力係数ＣＬが比較的大きく抗力係数ＣＤが比較的小さい翼形状として、「ＮＡＣＡ２４１５」、「ＮＡＣＡ４４１５」、「ＮＡＣＡ２３０１２」、「ＮＡＣＡ６３１−４１２」、「ＮＡＣＡ６４１−４１２」、「ＮＡＣＡ６５１−２１２」を基礎形状として選択する（図６参照）。ただし、翼厚比は、これらの翼の平均翼厚比である「0.129730」とした。
これらの翼形状は、従来から高性能であると認められている形状であり、予め設定された要求性能以上の性能を有する形状である。つまり、本実施形態においては、後に行うＰＯＤで用いられる複数の形状の全てが予め設定された要求性能以上の性能を有する形状とされている。
【００８６】
続いて、異なる表現手法で表現された基礎形状を同一の座標系で表現するための幾何学ベクトル変換が行われる（ステップＳ２３）。
ここでは、異なる表現手法で表現された基礎形状が同一の座標系における座標値で表現される幾何学ベクトル（ベクトル）に変換される。この幾何学ベクトルへの変換によって、全ての基礎形状が共通の表現手法によって表されることとなる。
具体的には、各基礎形状をｘｙ直交座標系に投影し、ｘ軸を１００分割し、翼正圧側と負圧側のそれぞれに対してｘ座標に対応するｙ座標を並べたベクトルを生成し、これらのベクトルをつなげて幾何学ベクトルを生成する。
なお、この幾何学ベクトル変換は、ＣＰＵ１がプログラム記憶部２ａに記憶された形状最適化プログラムＰ及びデータ記憶部２ｂに記憶された形状データＤ１等に基づいて生成する。
そして、この幾何学ベクトル変換（ステップＳ２３）は、本発明におけるベクトル変換工程に相当する。
【００８７】
続いて、上記幾何学ベクトルを用いてＰＯＤが行われる（ステップＳ２４）。
ＰＯＤの詳細については、上述したため、ここでの説明は省略する。なお、このＰＯＤは、ＣＰＵ１が、幾何学ベクトルやプログラム記憶部２ａに記憶された形状最適化プログラムＰ等に基づいて実行する。
そして、本実施形態においては、予め設定された要求性能以上の性能を有する形状を含む複数の形状を主成分分析することによって基底ベクトルを主成分として抽出する。
つまり、本実施形態においてＰＯＤを行う工程（ステップＳ２４）は、本発明の主成分抽出工程に相当する。
【００８８】
このように、本実施形態においては、ＰＯＤで用いられる複数の形状を同一の座標系におけるベクトルである幾何学ベクトルに変換し、変換後の幾何学ベクトルに基づいてＰＯＤが行われる。
そして、本実施形態においては、ＰＯＤを実行することによって、図７に示すような固有ベクトルと基底ベクトルを得た。
【００８９】
続いて、ＰＯＤの実行によって得られた基底ベクトルの選択が行われる（ステップＳ２５）。この基底ベクトルの選択数が少ないほど、新規形状の数を少なくして後の工程での計算量を削減することができる。一方で、パラメータ基底ベクトルは、次数の低いものほど、特徴の傾向を強く示している。したがって、基底ベクトルを選択する際には、低次のものから順に必要最低限の数を選択することが好ましい。
具体的には、選択する基底ベクトルの数は、例えば、基底ベクトルのノルムが予め定められた閾値に基づいて設定したり、ＣＰＵ１の処理能力（処理時間等）に基づいて設定したりすることが考えられる。
なお、本実施形態においては、図８に示すように、４次までの基底ベクトルのノルムが卓越している。このため、４次までの基底ベクトルを選択することが考えられるが、後の計算量を考慮して本実施形態においては３次までの基底ベクトルを選択する。
また、基底ベクトルの選択は、作業者が入力装置４を操作することで行う。また、予め記憶される閾値等に基づいてＣＰＵ１が自動で行うことも可能である。
そして、この基底ベクトルの選択（ステップＳ２５）は、本発明における基底ベクトル選択工程に相当する。
【００９０】
続いて、選択した基底ベクトルを用いて新規形状が生成される（ステップＳ２６）。
具体的には、下式（３４）で示される係数ベクトルに対して、式（２９）に従って新規形状を生成する。
なお、本実施形態においては、式（３２）で示すパラメータを、下式（３５）、下式（３６）、下式（３７）で示すように選択し、新規形状を１０８個生成した。
また、新規形状の生成は、ＣＰＵ１が、選択された基底ベクトルやプログラム記憶部２ａに記憶された形状最適化プログラムＰ等に基づいて実行する。
そして、この新規形状が生成される工程（ステップＳ２６）は、本発明における新規形状生成工程に相当する。
【００９１】
【数３４】
2013003971.tif 000037
【００９２】
【数３５】
2013003971.tif 000038
【００９３】
【数３６】
2013003971.tif 000039
【００９４】
【数３７】
2013003971.tif 000040
【００９５】
最後に、新規形状に対して、揚力係数ＣＬが大きく抗力係数ＣＤが算出されて評価が行われる（ステップＳ２７）。
なお、新規形状の評価は、ＣＰＵ１が、新規形状やプログラム記憶部２ａに記憶された形状最適化プログラムＰ等に基づいて実行する。
そして、この評価を行う工程（ステップＳ２７）は、本発明の評価工程に相当する。
【００９６】
具体的には、本実施形態においては、揚力係数ＣＬと抗力係数ＣＤの計算には、パネル法を実装したソルバを用いた。また、計算条件はレイノルズ数を下式（３８）とし、迎角５°とした。
【００９７】
【数３８】
2013003971.tif 000041
【００９８】
そして、新規形状を元に応答曲面を生成し、最適解を探索し、最終的に式（３３）が最大となる翼を最適解とする。
なお、本実施形態においては、応答曲面の生成に、「ModeFRONTIER ver.4.2.1」（「多目的ロバスト設計最適化支援ツール modeFRONTIER,http://www.cdaj.co.jp/product/020000modefrontier/index.html.」参照）を用いた。
【００９９】
図９は、元形状と新規形状の揚力係数ＣＬと抗力係数ＣＤの計算結果を示すグラフである。この図に示すように、揚力係数ＣＬと抗力係数ＣＤの分布は、元形状の値の周囲に広範囲に分布しており、元形状よりも高性能な形状が多く存在していることが分かる。
つまり、本実施形態の形状最適化方法及び形状最適化装置１０によれば、設計空間をうまく限定することができたことが分かる。
【０１００】
図１０は、本実施形態によって選択された最適解に基づく形状と、他の従来手法にて最適化された最適解に基づく形状とを示すグラフであり、（ａ）が本実施形態によって選択された最適解に基づく形状を示し、（ｂ）が他の手法にて最適化された最適解に基づく形状とを示す。
なお、ここでの他の従来手法とは、非特許文献３に示される手法であり、ベジェ曲線を用いて翼を表現し、「ＮＡＣＡ００１２ 」翼を元にして逐次二次計画法による形状最適化を行う手法である。
なお、本実施形態との比較を行うために、翼厚比を本実施形態の条件に合わせて変え、パネル法を用いて再計算した。そのため、引用元の論文とは形状と揚力係数ＣＬ、抗力係数ＣＤ、揚坑比が異なっている。なお、レイノルズ数と迎角の計算条件は引用元においても本実施形態の条件と同じである。
【０１０１】
以下の表２は、本実施形態によって選択された最適解に基づく形状と、他の手法にて最適化された最適解に基づく形状とにおける、揚力係数ＣＬ、抗力係数ＣＤ、揚坑比の比較である。
【０１０２】
【表２】
2013003971.tif 000042
【０１０３】
図１０及び表２に示すように、本実施形態によって選択された最適解に基づく形状と、他の手法にて最適化された最適解に基づく形状とにおいて、幾何形状は大きく異なっている。
そして、揚力係数ＣＬはどちらもほぼ同じであるのに対して、抗力係数ＣＤはＰＯＤを用いた本実施形態の方が小さい。このため、ＰＯＤを用いた本実施形態の方が、揚抗比が大きい。したがって、本実施形態の形状最適化方法及び形状最適化装置１０の方が、他の従来手法よりも、高性能な解を得ることができていると言える。
【０１０４】
また、本実施形態の最適解と基底ベクトルの関係について述べる。
本実施形態において最適解を生成する係数パラメータは、下式（３９）であった。
【０１０５】
【数３９】
2013003971.tif 000043
【０１０６】
そして、元形状として採用した基礎形状の翼はすべてキャンバーが比較的小さいが、選択解はキャンバーが大きいという特徴がある。また翼正圧側後縁が凹な形状をしている。
ここで、１次の基底ベクトルを見るとキャンバーが大きい。また、２次の基底ベクトルは正圧側後縁で凸になっている。ＰＯＤを用いることでキャンバーを持つ成分と正圧側後縁が凸な成分を独立に抜き出し、これにパラメータを乗じることでキャンバーが大きく後縁側で凹な形状をした選択解を生成することができたことがわかる。また、平板翼よりもキャンバーが大きく後縁側が凹な翼の方が、揚力係数ＣＬが大きくなることが知られておいる。このため、本実施形態の最適解のような翼を選択することは妥当である捉えることができる。
このように、本実施形態の形状最適化方法及び形状最適化装置１０によれば、作業者は基底ベクトルを確認することで選択解の妥当性を確認することができる。
【０１０７】
以上のように、本実施形態の形状最適化方法及び形状最適化装置１０によれば、予め設定された要求性能以上の性能を有する形状を含む複数形状から主成分を抽出し、主成分に基づいて最適化対象部（翼）の新規形状を生成する。
つまり、本実施形態は、コンピュータを用いて産業機械の一部あるいは全部である最適化対象部の形状を最適化する形状最適化方法であって、コンピュータがＣＰＵと形状最適化プログラム及び形状データを記憶する記憶装置とを備え、ＣＰＵが、形状最適化プログラム及び形状データに基づいて、予め設定された要求性能以上の性能を有する形状を含む複数の形状を主成分分析することによって主成分を抽出する主成分抽出工程と、前記主成分に基づいて前記最適化対象部の新規形状を生成する新規形状生成工程とを行うものである。
このような本実施形態の形状最適化方法及び形状最適化装置１０によれば、より高性能なものが含まれる可能性が高い新規形状を中心として評価を繰り返すことにより、最適解を探索することとなる。
したがって、本実施形態の形状最適化方法及び形状最適化装置１０によれば、産業機械の形状を最適化するにあたり、探索の精度が悪化することを抑制しつつ、短時間で最適解の探索を完了することが可能となる。
【０１０８】
また、本記実施形態においては、ＰＯＤ（ステップＳ２４）において基底ベクトルの抽出に用いられる形状が、全て高性能形状（予め設定された要求性能以上の性能を有する形状）であるため、生成される新規形状に、より高性能の形状が含まれる可能性が高くなる。
【０１０９】
また、本実施形態においては、複数の翼形状を同一の座標系におけるベクトルに変換してからＰＯＤを行っている。
このため、異なる表現手法で表現された翼形状同士から基底ベクトルを抽出することが可能となっている。
【０１１０】
また、本実施形態においては、ＰＯＤによって主成分分析を行っている。このため、主成分を基底ベクトルとして抽出することができる。
そして、抽出された基底ベクトルのうち、予め定められた数だけ低次側から選択し、性能への影響度合いが低い高次側を捨てることによって、最適解を探索するための計算量を削減することができる。
【０１１１】
（第２実施形態）
上記実施形態においては、形状が座標値で与えられている例について説明した。一方で、実際の設計においては、部材の肉厚などの決められた部分の寸法を最適化したい場合も多い。
そこで、本発明の第２実施形態として、吸吐弁における寸法を最適化することによって吸吐弁（本実施形態においてはレシプロ圧縮機吸吐弁）の形状の最適化を図る吸吐弁の形状最適化に本発明を適用した例について説明する。つまり、本実施形態は、産業機械の全部である吸吐弁を最適化対象部とするものである。
なお、本実施形態において用いる形状最適化装置のハードウェア構成は、上記第１実施形態と同様であるため、ここでの説明は省略する。
【０１１２】
図１１は、本実施形態において形状の最適化を図る吸吐弁の流路モデルである。なお、図１１においては、吸吐弁が開放されている状態の流路モデルである。そして、本実施形態においては、図１１に示す１０個の寸法を設計変数とした。なお、各設計変数には、加工限界等に基づく制約条件を課している。
図１１に示す吸吐弁では、流路が合流する部分（図１１において最も拡大されている部分）でバックステップ流れがあり、全圧損失が高い。そこで、本実施形態においては、流路形状に対して最適化を行い、上記全圧損失を小さくすることを目標とする。
【０１１３】
本実施形態においては、設計変数の値を入力ベクトル（ベクトル）としてＰＯＤを実行する。具体的には、タグチメソッドを用いて得られた圧力損失が比較的小さい形状（高性能形状）を５つを用いてＰＯＤを実行し、基底ベクトルを得る。
この結果、基底ベクトルのノルムは、図１２に示すように、３次までの基底ベクトルが卓越している。
そこで、本実施形態においては、３次までの基底ベクトルを用いて式（２９）に従って新たにパラメトリックに４８個の設計変数を生成する。
なお、式（３２）に示すパラメータは、上記第１実施形態で説明したように設定する。具体的には、パラメータを下式（４０）及び下式（４１）に示すように選択する。
【０１１４】
【数４０】
2013003971.tif 000044
【０１１５】
【数４１】
2013003971.tif 000045
【０１１６】
続いて、４８個の設計変数のうち、制約条件を満たさないものを取り除き、残ったものを元に新規形状を生成する。そして、「Fluent」を用いて新規形状における全圧損失の値を求め、得られた全圧損失の値と式（３２）の情報を元に応答曲面を生成し、応答曲面のみを用いて最適解を探索する。
なお、本実施形態において探索された最適解は応答曲面の生成点であったため、再計算は行う必要がなかった。応答曲面を生成する際は、上記実施形態と同様に、「ModeFRONTIER ver.4.2.1」を用いた。
【０１１７】
表３は、本実施形態による最適解と、他の従来手法による最適解とにおける全圧損失係数との結果を示す表である。なお、表３において、ζが全圧損失係数を示す。また、他の従来手法としては、遺伝的アルゴリズム（ＧＡ）を用いた手法を用いた。また、遺伝的アルゴリズムを実行したときの条件は表４に示す。
【０１１８】
【表３】
2013003971.tif 000046
【０１１９】
【表４】
2013003971.tif 000047
【０１２０】
従来手法を用いた場合は、全圧損失係数が２７．３である。一方で、本実施形態を用いた場合は全圧損失係数が２５．６となっている。つまり、ＰＯＤを用いることでより高性能な解が探索できたことが確認できる。
このように限られた設計変数を元に形状が決まる場合に関してもＰＯＤを用いた本実施形態の最適化方法を用いて高性能な解を求めることができる。
【０１２１】
上述の第１実施形態及び第２実施形態では、まず、流体性能の良い構造物の幾何情報の集合を元にしてＰＯＤによりそれらを基底ベクトルに分解する。次に低次の基底ベクトルからパラメトリックに新規形状を表現し、応答曲面法などの非線形最適化手法を用いて最適解を探索する。
この手法を用いることで、上述のように、性能の良い形状に共通する特徴を持った形状が得られる。またパラメータを減らせるため、最適化を行う際の探索領域が狭くなり、発見的手法を用いた場合に、より良い解が求められることができる。
さらに低次の主成分（すなわち低次の基底ベクトル）を分析することで、高性能な形状の特徴を知ることができる。このためブラックボックス化しがちな最適設計の過程が作業者にわかりやすく、ノウハウの蓄積にも繋がる。
また逆に、これまでに得られている高性能な形状を用いて新たな形状を生成することで知識の有効利用が可能である。
【０１２２】
ＰＯＤを用いることの最大の利点は探索領域が小さくなることである。例えば、上記第１実施形態の翼形状の最適化では設計変数を３個のパラメータに減少させた場合でも、従来以上の高性能な形状が得られることが確認できた。この利点は特に設計変数の多い大規模な問題において効果が大きい。
【０１２３】
以上、添付図面を参照しながら本発明の好適な実施形態について説明したが、本発明は、上記実施形態に限定されないことは言うまでもない。上述した実施形態において示した各構成部材の諸形状や組み合わせ等は一例であって、本発明の趣旨から逸脱しない範囲において設計要求等に基づき種々変更可能である。
【０１２４】
例えば、上記実施形態においては、本発明を翼あるいは吸吐弁の形状最適化に適用した構成について説明した。
しかしながら、本発明はこれに限定されるものではなく、流体に関連する産業機械の全体あるいは一部の形状であれば、どのようなものであっても最適化を図ることができる。
例えば、インピンジ冷却に用いるインピンジ孔の形状を本発明によって最適化することも可能である。
【０１２５】
また、上記実施形態においては、主成分分析がＰＯＤである構成について説明した。しかしながら、本発明はこれに限定されるものではなく、主成分分析として、カルーネン・レーベ展開、特異値分解、スペクトル分解、独立成分分析あるいは因子分析を用いることも可能である。
【０１２６】
また、上記実施形態においては、ＰＯＤ（ステップＳ２４）において基底ベクトルの抽出に用いられる形状が、全て高性能形状である構成を採用した。
しかしながら、本発明はこれに限定されるものではなく、複数の形状のうち、いくつかの形状を高性能形状と異なる形状としても良い。これによって、思わぬ高性能な新規形状が生成される可能性もある。
【符号の説明】
【０１２７】
１０……形状最適化装置、１……ＣＰＵ（主成分抽出手段、新規形状生成手段、評価手段）、２……記憶装置、３……リムーバブルメディアドライブ、４……入力装置、５……出力装置、６……通信装置
【図１】
2013003971.tif 000048
【図２】
2013003971.tif 000049
【図３】
2013003971.tif 000050
【図４】
2013003971.tif 000051
【図５】
2013003971.tif 000052
【図６】
2013003971.tif 000053
【図７】
2013003971.tif 000054
【図８】
2013003971.tif 000055
【図９】
2013003971.tif 000056
【図１０】
2013003971.tif 000057
【図１１】
2013003971.tif 000058
【図１２】
2013003971.tif 000059