6248402 データの表示方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6248402

(24)【登録日】2017年12月1日

(45)【発行日】2017年12月20日

(54)【発明の名称】データの表示方法

(51)【国際特許分類】

   G06F 17/50 20060101AFI20171211BHJP

【ＦＩ】

   G06F17/50 612C

   G06F17/50 612A

   G06F17/50 680Z

【請求項の数】6

【全頁数】17

(21)【出願番号】特願2013-56804(P2013-56804)

(22)【出願日】2013年3月19日

(65)【公開番号】特開2014-182605(P2014-182605A)

(43)【公開日】2014年9月29日

【審査請求日】2016年3月10日

(73)【特許権者】

【識別番号】000006714

【氏名又は名称】横浜ゴム株式会社

(74)【代理人】

【識別番号】100080159

【弁理士】

【氏名又は名称】渡辺 望稔

(74)【代理人】

【識別番号】100090217

【弁理士】

【氏名又は名称】三和 晴子

(74)【代理人】

【識別番号】100152984

【弁理士】

【氏名又は名称】伊東 秀明

(74)【代理人】

【識別番号】100148080

【弁理士】

【氏名又は名称】三橋 史生

(72)【発明者】

【氏名】小石 正隆

【審査官】 松浦 功

(56)【参考文献】

【文献】 特開２００５−０７０８４９（ＪＰ，Ａ）

【文献】 特開２０１３−０１２１８５（ＪＰ，Ａ）

【文献】 特開２０１２−０４５９７５（ＪＰ，Ａ）

【文献】 特開２００４−１１８７１９（ＪＰ，Ａ）

【文献】 特開２０１０−００９３４１（ＪＰ，Ａ）

【文献】 片田喜章，進化ロボティクスにおける多目的ＧＡに関する一考察，第５４回システム制御情報学会研究発表講演会講演論文集［ＣＤ−ＲＯＭ］，システム制御情報学会，２０１０年 ５月１９日，ｐｐ．６０９−６１０

【文献】 山地秀美 外２名，２目的ＧＡにおけるパレート解進化の視覚化システム，（社）日本経営工学会平成１７年度春季大会予稿集，社団法人日本経営工学会，２００５年 ５月１２日，ｐｐ．５６−５７

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０６Ｆ １７／５０

(57)【特許請求の範囲】

【請求項1】

構造体および構造体を構成する材料を規定するパラメータのうち設計変数として定められた少なくとも１つのパラメータと、構造体および構造体を構成する材料を規定するパラメータのうち特性値として定められた少なくとも１つのパラメータとの２種類のデータを対象としたデータの表示方法であって、
コンピュータが、前記設計変数と前記特性値との間の非線形応答関係を定める第１の工程と、
前記設計変数の定義域を定め、前記第１の工程で定めた非線形応答関係を用いて、特性値を目的関数とする最適化を実施しパレート解を算出する第２の工程と、
前記パレート解を目的関数空間で散布図として表示する際、前記設計変数の値に応じ、前記散布図で前記設計変数の値を表すシンボルを、その色、種類および大きさのうち、少なくとも１つを変えて表示する第３の工程を実行し、
前記コンピュータが、更に特性値として定められたパラメータを追加し、前記設計変数と前記パラメータが追加された特性値との間の非線形応答関係を定め、前記設計変数の定義域を定め、その非線形応答関係を用いて、特性値を目的関数とする最適化を実施して、拡張パレート解を算出する工程を実行し、
前記第３の工程において、前記コンピュータが、前記拡張パレート解を前記パレート解とともに前記目的関数空間で散布図として表示することを特徴とするデータの表示方法。

【請求項2】

前記コンピュータが、前記第２の工程と前記第３の工程の間に、
前記設計変数の前記定義域内で少なくとも１つの離散値を設定し、前記離散値に基づき算出された前記パレート解を縮約する工程を実行し、
前記第３の工程は、前記コンピュータが、前記縮約されたパレート解の値を表すシンボルを、シンボルの色、種類および大きさのうち、少なくとも１つを変えて前記散布図として表示する請求項１に記載のデータの表示方法。

【請求項3】

前記第２の工程は、前記コンピュータが、前記設計変数の前記定義域内で少なくとも１つの離散値を設定し、前記離散値に基づき、前記第１の工程で定めた非線形応答関係を用いて前記パレート解を算出する請求項１に記載のデータの表示方法。

【請求項4】

前記第２の工程で、前記コンピュータが、前記パレート解を算出する際に、前記パレート解の探索過程で得られる特性値を目的関数データとして保持しておき、
前記第３の工程で、前記コンピュータが、前記目的関数データを前記パレート解とともに前記目的関数空間で散布図として表示する請求項１〜３のいずれか１項に記載のデータの表示方法。

【請求項5】

前記第３の工程は、前記コンピュータが、更に、前記設計変数の値毎にパレートフロントにあるパレート解を通る線を表示する工程を実行する請求項１〜４のいずれか１項に記載のデータの表示方法。

【請求項6】

前記設計変数は、タイヤの材料挙動、タイヤの形状およびタイヤの構造のうち、少なくとも１つのパラメータであり、前記特性値は、タイヤの特性値である請求項１〜５のいずれか１項に記載のデータの表示方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、所定の設計因子（設計変数）とその特性値からなるデータの表示方法に関し、特に、パレート解を目的関数空間で散布図表示し、設計変数の値に応じて散布図上で設計変数の値として表わされるシンボルを、その色、種類および大きさのうち、少なくとも１つを変化させて表示させるデータの表示方法に関する。

【背景技術】

【0002】

シミュレーションでは、所定の設計因子、すなわち、設計変数に対する応答、すなわち、特性値を算出することができる。また、最適化では所望の特性値を得るための設計変数の値を算出することができる。
一方、設計探査と呼ばれる手法は、最適化とデータマイニングを用いて設計変数と特性値との因果関係を見出そうとする手法である。設計探査では、複数の特性値、すなわち、目的関数を対象とし、それらにトレードオフ関係が存在することが少なくない。したがって、パレート解と呼ばれる最適解は集合を形成する。そのため、パレート解のデータを、特性値の側面と設計変数の側面の両面から分析するため、自己組織化マップと呼ばれる手法が用いられている（特許文献１、非特許文献１参照）。

【0003】

特許文献１には、タイヤ等の構造体の数値シミュレーションを用いて、目的関数に対するパレート解の情報を付与した自己組織化マップを生成するとともに、設計変数の情報を付与した自己組織化マップを生成していることが記載されている。このため、設計者は、これらのマップの面上のパレート解の全体像を見ながら、トレードオフの関係にある性能バランスを考慮しつつ、自己組織化マップ上の位置を定めることにより最適設計案を決定することができることが、特許文献１には記載されている。

【0004】

また、非特許文献１では、自己組織化マップを用い、この自己組織化マップ上に目的関数、設計変数を表示することができ、それらを並べて表示することで、視覚的に目的関数間の相関関係を把握できるだけでなく、目的関数と設計変数との因果関係も理解できることが記載されている。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】特許第４３３９８０８号公報

【非特許文献】

【0006】

【非特許文献1】小石正隆、日本ゴム協会誌、Vol.85、2012、289-295.

【発明の概要】

【発明が解決しようとする課題】

【0007】

しかしながら、特許文献１および非特許文献１に開示されている自己組織化マップは、目的関数間のパレートフロントに沿った分析、およびパレートフロント近傍に着目した分析に適していないという問題点がある。一方、目的関数空間でのパレート解の散布図表示は、パレートフロント、およびその近傍に着目し易い反面、設計変数との因果関係を見出し難く、更にデータの重なりが多く視認性が悪いという問題点がある。

【0008】

本発明の目的は、前記従来技術に基づく問題点を解消し、設計変数と特性値（目的関数）との因果関係を容易に理解しやすい散布図を表示することができるデータの表示方法を提供することにある。

【課題を解決するための手段】

【0009】

上記目的を達成するために、構造体および構造体を構成する材料を規定するパラメータのうち設計変数として定められた少なくとも１つのパラメータと、構造体および構造体を構成する材料を規定するパラメータのうち特性値として定められた少なくとも１つのパラメータとの２種類のデータを対象としたデータの表示方法であって、設計変数と特性値との間の非線形応答関係を定める第１の工程と、設計変数の定義域を定め、第１の工程で定めた非線形応答関係を用いて、特性値を目的関数とする最適化を実施しパレート解を算出する第２の工程と、パレート解を目的関数空間で散布図として表示する際、設計変数の値に応じ、散布図で設計変数の値を表すシンボルを、その色、種類および大きさのうち、少なくとも１つを変えて表示する第３の工程を有することを特徴とするデータの表示方法を提供するものである。

【0010】

第２の工程と第３の工程の間に、設計変数の定義域内で少なくとも１つの離散値を設定し、離散値に基づき算出されたパレート解を縮約する工程を有し、第３の工程は、縮約されたパレート解の値を表すシンボルを、シンボルの色、種類および大きさのうち、少なくとも１つを変えて散布図として表示することが好ましい。
第２の工程は、設計変数の定義域内で少なくとも１つの離散値を設定し、離散値に基づき、第１の工程で定めた非線形応答関係を用いてパレート解を算出することが好ましい。
第２の工程で、パレート解を算出する際に、パレート解の探索過程で得られる特性値を目的関数データとして保持しておき、第３の工程で、目的関数データをパレート解とともに目的関数空間で散布図として表示することが好ましい。

【0011】

更に特性値として定められたパラメータを追加し、設計変数とパラメータが追加された特性値との間の非線形応答関係を定め、設計変数の定義域を定め、その非線形応答関係を用いて、特性値を目的関数とする最適化を実施して、拡張パレート解を算出する工程を備え、第３の工程において、拡張パレート解をパレート解とともに目的関数空間で散布図として表示することが好ましい。
第３の工程は、更に、設計変数の値毎にパレートフロントにあるパレート解を通る線を表示する工程を備えることが好ましい。
例えば、設計変数は、タイヤの材料挙動、タイヤの形状およびタイヤの構造のうち、少なくとも１つのパラメータであり、特性値は、タイヤの特性値である。

【発明の効果】

【0012】

本発明によれば、設計変数と特性値（目的関数）との因果関係を容易に理解しやすい散布図を表示することができる。

【図面の簡単な説明】

【0013】

【図1】本発明の実施形態の表示方法に利用される表示処理装置を示す模式図である。

【図2】本発明の実施形態の表示方法の第１の例を工程順に示すフローチャートである。

【図3】（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解の一例を示す散布図であり、（ｂ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解の他の例を示す散布図である。

【図4】（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解について設計変数の一部だけ表示を変えた一例を示す散布図であり、（ｂ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解について設計変数の一部だけ表示を変えた他の例を示す散布図である。

【図5】本発明の実施形態の表示方法の第２の例を工程順に示すフローチャートである。

【図6】（ａ）は、設計変数の定義域を説明する模式図であり、（ｂ）は、設計変数の定義域の離散値の一例を説明する模式図であり、（ｃ）は、設計変数の定義域の離散値の他の例を説明する模式図である。

【図7】（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解を示す散布図であり、（ｂ）は、図７（ａ）に示すパレート解を設計変数の定義域の離散値に応じて縮約した例を示す散布図である。

【図8】（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解を示す散布図であり、（ｂ）は、図８（ａ）に示すパレート解を設計変数の定義域の離散値に応じて縮約した例を示す散布図である。

【図9】（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解を示す散布図であり、（ｂ）は、図９（ａ）に示すパレート解にパレートフロントを明確にした例を示す散布図である。

【図10】本発明の実施形態の表示方法の第３の例を工程順に示すフローチャートである。

【図11】（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解を示す散布図であり、（ｂ）は、図１１（ａ）に示すパレート解にパレート解探索過程の特性値データを追加した例を示す散布図である。

【図12】本発明の実施形態の表示方法の第４の例を工程順に示すフローチャートである。

【図13】（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解を示す散布図であり、（ｂ）は、図１３（ａ）に示すパレート解に、他の特性値を考慮したパレート解を追加した例を示す散布図である。

【発明を実施するための形態】

【0014】

以下に、添付の図面に示す好適実施形態に基づいて、本発明のデータの表示方法を詳細に説明する。
図１は本発明の実施形態の表示方法に利用される表示処理装置を示す模式図である。図２は本発明の実施形態の表示方法の第１の例を工程順に示すフローチャートである。

【0015】

図１に示す表示処理装置１０は、本発明の実施形態の表示方法を実施に用いられる装置の一例である。表示処理装置１０は、コンピュータ等のハードウェアを用いて構成されるものである。
表示処理装置１０は、処理部１２と、入力部１４と、表示部１６とを有する。処理部１２は、条件設定部２０、モデル生成部２２、演算部２４、パレート解探索部２６、メモリ２８、表示制御部３０および制御部３２を有する。この他に図示はしないがＲＯＭ等を有する。
処理部１２は、制御部３２により制御される。また、処理部１２において条件設定部２０、モデル生成部２２、演算部２４およびパレート解探索部２６はメモリ２８に接続されており、条件設定部２０、モデル生成部２２、演算部２４およびパレート解探索部２６のデータがメモリ２８に記憶される。

【0016】

入力部１４は、マウスおよびキーボード等の各種情報をオペレータの指示により入力するための各種の入力デバイスである。表示部１６は、例えば、本発明の表示方法で得られた図を表示するものであり、公知の各種のディスプレイが用いられる。また、表示部１６には各種情報を出力媒体に表示するためのプリンタ等のデバイスも含まれる。

【0017】

表示処理装置１０は、ＲＯＭ等に記憶されたプログラム（コンピュータソフトウェア）を、制御部２６で実行することにより、条件設定部２０、モデル生成部２２、演算部２４、パレート解探索部２６の各部を機能的に形成する。表示処理装置１０は、上述のように、プログラムが実行されることで各部位が機能するコンピュータによって構成されてもよいし、各部位が専用回路で構成された専用装置であってもよい。

【0018】

条件設定部２０は、本実施形態の表示方法により、パレート解を目的関数空間で散布図として表示する際に必要な各種の条件、情報が入力され、設定する。各種の条件、情報は、入力部１２を介して入力される。条件設定部２０で設定する各種の条件、情報はメモリ２８に記憶される。

【0019】

条件設定部２０には、構造体および構造体を構成する材料を規定するパラメータのうち設計変数として定めた少なくとも１つのパラメータが設定される。なお、設計変数のパラメータには、荷重、境界条件などのばらつき因子を設定してもよい。
また、構造体および構造体を構成する材料を規定するパラメータのうち特性値（目的関数）として定めた少なくとも１つのパラメータが設定される。特性値には、コスト等の物理的および化学的な特性値以外の、構造体および構造体を構成する材料を評価する指標を用いてもよい。
構造体および構造体を構成する材料は、構造体単体ではなく、構造体を構成するパーツ、構造体のアッセンブリ形態等の構造体を含むシステム全体、またはその一部を対象としてもよい。

【0020】

条件設定部２０に設定される特性値は、評価しようとする物理量である。目的関数は、評価しようとする物理量を求めるための関数である。
構造体がタイヤである場合、特性値はタイヤの特性値である。この場合、特性値としては、タイヤ性能として評価しようとする物理量であり、例えば、操縦安定性の指標となるスリップ角１度における横力であるＣＰ（コーナリングパワー）、乗心地性の指標となるタイヤの１次固有振動数、転動抵抗の指標となる転がり抵抗、操縦安定性の指標となる横ばね定数、耐摩耗性の指標となるタイヤトレッド部材の摩耗エネルギ等が挙げられる。目的関数は、それらを求めるための関数である。目的関数は、性能として好ましい方向があり、値が大きくなる、小さくなる、または所定の値に近づく等がある。

【0021】

設計変数は、構造体の形状、構造体の内部構造および材料特性等を規定するものである。タイヤの場合、設計変数は、タイヤの材料挙動、タイヤの形状およびタイヤの構造のうち、少なくとも１つのパラメータである。設計変数としては、例えば、タイヤのトレッド部におけるクラウン形状を規定する曲率半径、タイヤ内部構造を規定するタイヤのベルト幅寸法等が挙げられる。これ以外にも、例えば、トレッド部における材料特性を規定するフィラー分散形状、フィラー体積率等が挙げられる。
制約条件は、目的関数の値を所定の範囲に制約したり、設計変数の値を所定の範囲に制約するための条件である。
また、構造体がタイヤである場合、タイヤの負荷荷重、タイヤの転動速度を初めとする走行条件、タイヤが走行する路面条件、例えば、凹凸形状、摩擦係数等、車両の走行シミュレーションに用いるための車両諸元の情報等が設定される。

【0022】

また、条件設定部２０に、設計変数のパラメータと特性値のパラメータとの間の非線形応答関係を定めるための情報が設定される。この非線形応答関係には、例えば、ＦＥＭ等の数値シミュレーション、理論式および近似式等が含まれる。
条件設定部２０では、非線形応答関係により生成するモデル、そのモデルの境界条件、ＦＥＭ等の数値シミュレーションする場合には、そのシミュレーション条件、シミュレーションにおける制約条件を設定する。更には、パレート解を得るための最適化条件、例えば、パレート解探索のための条件等を設定する。

【0023】

パレート解探索のための条件は、パレート解を探索するための手法、パレート解探索における各種条件である。本実施形態では、例えば、パレート解を探索するための手法として、遺伝的アルゴリズムを用いることができる。一般に、目的関数の増大と共に、遺伝的アルゴリズムの探査能力が低下することが知られている。それを解決する方法の一つが、個体数を増加させる方法である。一方、個体数を増加させ、パレート解を探査すると、多くのパレート解が算出される。したがって、多くの特性値データと設計パラメータとの因果関係を視認性良く表示する方法が設計探査の一つの課題となっているが、本発明ではこれを解決することができる。
これ以外に、条件設定部２０で、設計変数の定義域を設定する。また、条件設定部２０では後述するようにパレート解を縮約する際に用いられる離散値を設定する。

【0024】

モデル生成部２２は、設定された非線形応答関係に基づいて、各種の計算モデルを作成するものである。非線形応答関係は、上述のようにＦＥＭ等の数値シミュレーションが含まれており、この場合、モデル生成部２２で、設計変数を表わす設計パラメータ、特性値を表わす特性値パラメータに応じたメッシュモデルが生成される。また、理論式および近似式等の場合にも、設計パラメータ、特性値パラメータに応じた理論式および近似式等が作成される。なお、構造体がタイヤの場合には、タイヤモデルが作成される。演算部２４でタイヤモデルを用いてシミュレーション演算がなされる。

【0025】

なお、モデル生成部２２で作成されるタイヤモデルは、条件設定部２０で設定された各種類の設計パラメータを用いて作成されるが、タイヤモデルの作成には公知の作成方法を用いることができる。なお、タイヤモデルは、少なくとも、このタイヤモデルを転動させる対象である路面モデルも併せて生成する。また、タイヤが装着されるリム、ホイール、およびタイヤ回転軸を再現するものをタイヤモデルとしてもよい。また、必要に応じて、タイヤが装着される車両を再現するモデルをタイヤモデルに組み込んでもよい。この際、タイヤモデル、リムモデル、ホイールモデル、およびタイヤ回転軸モデルを、予め設定された境界条件に基づいて一体化したモデルを作成することもできる。

【0026】

これら各モデルは数値計算可能な離散化モデルであればよく、例えば、公知の有限要素法（ＦＥＭ）に用いるための有限要素モデル等であればよい。なお、タイヤモデルを用いて、例えば、タイヤウエット性能を初めとするタイヤ性能を最適化するタイヤ設計案を求める場合など、路面モデルとタイヤモデルの他に、路面上に存在する介在物を再現するモデルを生成しておけばよい。例えば、介在物モデルとして、路面上の水、雪、泥、砂、砂利、氷等を再現する各種モデルを、数値計算可能な離散化モデルで生成しておけばよい。なお、路面モデルも、表面が平坦な路面を再現するモデルに限らず、必要に応じて、表面に凹凸を有する路面形状を再現するモデルであってもよい。

【0027】

演算部２５は、モデル生成部２２で作成された各種のモデルを用いて特性値を算出するものである。これにより、設定変数に対する特性値が得られる。この特性値の中に、パレート解が存在する。得られた特性値は、メモリ２８に記憶される。
演算部２５では、例えば、路面上を転動するタイヤの転動を再現するシミュレーション条件を、モデル生成部２２で生成したタイヤモデル、または路面モデル等に与えたときの、タイヤモデルの挙動、またはタイヤモデルに作用する力などの物理量を時系列に求める。演算部２５は、例えば、公知の有限要素ソルバーによるサブルーチンを実行することで機能するものである。
また、演算部２５では、モデル生成部２２で理論式および近似式等を作成した場合には、理論式および近似式等を解き、特性値を算出する。

【0028】

パレート解探索部２６は、条件設定部２０で設定されたパレート解探索の条件に応じて、演算部２５で得られた特性値の中から、パレート解を探索し、パレート解を算出するものである。得られたパレート解は、メモリ２８に記憶される。

【0029】

ここで、パレート解は、トレードオフの関係にある複数の目的関数において、他の任意の解よりも優位にあるとはいえないが、より優れた解が他に存在しない解をいう。一般にパレート解は集合として複数個存在する。
パレート解探索部２６は、例えば、遺伝的アルゴリズムを用いてパレート解を探索する。
遺伝的アルゴリズムとしては、例えば、解集合を目的関数に沿って複数の領域に分割し、この分割した解集合毎に多目的ＧＡを行うＤＲＭＯＧＡ（Divided Range Multi-Objective GA）、ＮＣＧＡ（Neighborhood Cultivation GA），ＤＣＭＯＧＡ（Distributed Cooperation model of MOGA and SOGA）、ＮＳＧＡ（Non-dominated Sorting GA）、ＮＳＧＡ２（Non-dominated Sorting GA-II）、ＳＰＥＡII（Strength Pareto Evolutionary Algorithm-II）法等の公知の方法を用いることができる。その際、解集合が解空間に幅広く分布し、精度の高いパレート解の集合を求める必要がある。このため、パレート解探索部２６では、例えば、ベクトル評価遺伝的アルゴリズム（Vector Evaluated Generic Algorithms：ＶＥＧＡ）、パレートランキング法、またはトーナメント法を用いた選択が行われる。

【0030】

本発明では、設計変数と特性値との間で定める非線形応答関係、すなわち、設計変数を用いて特性値を求める場合に利用されるものは、ＦＥＭ等のシミュレーションに限定されるものではなく、上述のように理論式および近似式等を用いることもできる。例えば、シミュレーションモデルを用いた演算ではなく、シミュレーション近似式を用いて目的関数の値を算出してもよい。この場合、実験計画法に基づいて得られる実験結果から設計変数と目的関数との間の近似式、例えば、シミュレーション近似式を用いてパレート解を得ることができる。このシミュレーション近似式としては、多項式またはニューラルネットワーク等により得られる公知の非線形関数を用いることができる。

【0031】

表示制御部３０は、得られたパレート解について、設計変数のパラメータの値に応じて、設計変数のパラメータの値を表すシンボルの色、種類および大きさのうち、少なくとも１つを変えるものである。表示形態を変更したパレート解の情報はメモリ２８に記憶される。得られたパレート解は、表示制御部３０で表示形態が変えられて表示部１６で表示される。また、表示制御部３０では、設計変数のパラメータの値毎に、そのパレート解を結んだ線を表示させる機能も有する。

【0032】

次に、本実施形態の表示方法について説明する。
本実施形態の表示方法は、図１に示す表示処理装置１０を用いられるが、表示方法をコンピュータ等のハードウェアおよびソフトウェアを用いて実行することができれば、表示処理装置１０に限定されるものではない。
まず、対象となる構造体について設計変数および特性値を設定する。本実施形態では、構造体を、例えば、タイヤとした。なお、タイヤのサイズは１９５／６５Ｒ１５である。
タイヤに対して、設計変数として、タイヤの形状パラメータを設定する。そして、特性値として、ころがり抵抗、横ばね定数の２つを設定する。本実施形態では、入力がタイヤの形状パラメータであり、出力がころがり抵抗と横ばね定数となる。タイヤの形状パラメータの値により、ころがり抵抗と横ばね定数がどのように変化するかを表示する。タイヤの形状パラメータ、ころがり抵抗と横ばね定数が条件設定部２０に設定される。

【0033】

次に、図２に示すように、設計変数から特性値を求める際に用いる非線形応答を定める（ステップＳ１０）。すなわち、設計変数と特性値との関係を定める。この非線形応答の種類は、例えば、メモリ２８に記憶される。具体的には、タイヤの形状パラメータと、ころがり抵抗と横ばね定数との関係を設定する。タイヤの形状パラメータを入力とし、ころがり抵抗または横ばね定数を出力とした場合、設定する関係は、例えば、ころがり抵抗がタイヤの形状パラメータを変数とする二次多項式等の非線形関数を用いて表わされるものである。また、横ばね定数がタイヤの形状パラメータを変数とする二次多項式等の非線形関数を用いて表わされるものである。

【0034】

次に、設計変数の定義域を設定する（ステップＳ１２）。この場合、設計変数のパラメータに対して、上限値と下限値を設定し、下限値〜上限値の間が連続であるとする。例えば、タイヤの形状パラメータであれば、サイズの上限と下限を、下限値〜上限値の間が連続であるとして、設計変数の定義域として設定する。また、タイヤのゴム組成であれば、弾性率の上限と下限を設計変数の定義域として設定する。この設計変数の定義域の設定は、条件設定部２０でなされ、例えば、メモリ２８に記憶される。本実施形態では、タイヤの形状パラメータについて上限値と下限値を設定する。

【0035】

次に、非線形応答関係に基づいてモデル作成部２２でモデル作成を実施し、演算部２４にてステップＳ１０で設定した非線形応答関係に基づいて特性値を算出する（ステップＳ１４）。このとき、設定した設計変数の定義域がメモリ２８から読み出されて特性値が算出される。特性値の算出結果は、例えば、メモリ２８に記憶される。ＦＥＭ等のシミュレーションであれば、メッシュモデルがモデル作成部２２で作成され、演算部２４にて、ＦＥＭ等により入力に対する応答をシミュレーションが実施される。具体的には、タイヤの形状パラメータに対するころがり抵抗と横ばね定数が算出される。

【0036】

次に、パレード解探索部２６にて特性値の演算結果に対して、特性値を目的関数とする最適化を実施し、パレート解を得る（ステップＳ１６）。このパレート解の算出には、例えば、遺伝的アルゴリズムが用いられる。得られたパレート解は、メモリ２８に記憶される。

【0037】

次に、表示制御部３０にて、設計変数の値に応じて、設計変数の値を表すシンボルの色、種類および大きさのうち、少なくとも１つを変える。具体的には、設計変数Ｘ１、Ｘ２について、例えば、シンボルの種類を◇、＋、△、×、□、○とする。なお、シンボルの色、種類および大きさについては、特に限定されるものではない。
そして、図３（ａ）に示すように、縦軸に目的関数（特性値）Ｙ１をとり、横軸に目的関数（特性値）Ｙ２をとって設計変数Ｘ１の値を表すシンボルの種類（◇、＋、△、×、□、○）と色を変えてパレート解を、目的関数空間で散布図として表示部１６に示す（ステップＳ１８）。例えば、図３（ａ）では、領域Ｄにあるパレート解の方が好ましい。
なお、図３（ａ）において、目的関数（特性値）Ｙ１がタイヤの横ばね定数、目的関数（特性値）Ｙ２がタイヤの転がり抵抗、設計変数Ｘ１がタイヤの形状パラメータである。

【0038】

本実施形態では、図３（ａ）に示すように、パレート解を設計変数に対して層別化することができる。これにより、設計変数と特性値との関係を認識しやすくなる。設計変数Ｘ１の値毎に、パレート解は散布図上での位置が変わることが分かる。このように、タイヤの形状パラメータの変化による特性値の変化を容易に認識することができ、ひいては設計方針の指標を得ることができる。

【0039】

なお、図３（ａ）以外にも、設計変数Ｘ１のタイヤの形状パラメータに代えて、他のタイヤの形状パラメータである設計変数Ｘ２と、特性値としてころがり抵抗（目的関数（特性値）Ｙ２）と横ばね定数（目的関数（特性値）Ｙ１）を用いて、上述のステップＳ１０〜ステップＳ１８を実施することにより、図３（ｂ）に示すようにパレート解を表示することができる。
図３（ｂ）に示すパレート解も、設計変数Ｘ２の値を表すシンボルの種類（◇、＋、△、×、□、○）と色を変えて目的関数空間で散布図として表示部１６に示すこともできる。

【0040】

本発明において、表示方法は、図３（ａ）、（ｂ）に限定されるものに限定されるものではない。例えば、設計変数の値のうち、ある特定の値のものだけを、シンボルの色、種類および大きさのうち、少なくとも１つを変えて表示してもよい。この場合、図４（ａ）、（ｂ）に示す散布図が表示部１６に表示される。表示方法としては、得られたパレート解に対して、設計変数の値のうち、特定の値のものの色と大きさを変える。そして、表示部１６に表示する。これにより、設計変数の値のうち、特定の値のものが、散布図上のどこにあるかを明確に認識させることができる。
なお、図４（ａ）では、図３（ａ）のシンボル○を特定の値のものとして表示形態を●に変え、図４（ｂ）では、図３（ｂ）のシンボル◇を特定の値のものとして表示形態を●に変えて表示したものである。

【0041】

次に、本発明の表示方法の第２の例について説明する。
図５は、本発明の実施形態の表示方法の第２の例を工程順に示すフローチャートである。図６は、（ａ）は、設計変数の定義域を説明する模式図であり、（ｂ）は、設計変数の定義域の離散値の一例を説明する模式図であり、（ｃ）は、設計変数の定義域の離散値の他の例を説明する模式図である。
図７は、（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解を示す散布図であり、（ｂ）は、図７（ａ）に示すパレート解を設計変数の定義域の離散値に応じて縮約した例を示す散布図である。図７（ａ）は、図３（ａ）と同じである。
図８は、（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解を示す散布図であり、（ｂ）は、図８（ａ）に示すパレート解を設計変数の定義域の離散値に応じて縮約した例を示す散布図である。図８（ａ）は、図３（ｂ）と同じである。

【0042】

表示方法の第２の例は、図２に示す表示方法の第１の例に比して、パレート解を縮約した後に、パレート解を表示する点以外は、第１の例と同様の工程であるため、その詳細な説明は省略する。

【0043】

第２の例の表示方法では、表示処理装置１０において、パレート解探索部２６でパレート解を算出した後、このパレート解を演算部２４にて、後述する設計変数の離散値に基づいてパレート解を離散的に縮約する。結果としてパレート解が間引かれる。間引かれたパレート解が表示部１６に上述の第１の例と同様にして表示される。なお、縮約されたパレート解について表示制御部３０で表示形態が変更され、メモリ２８に記憶されるとともに、表示部１６に表示される。
図５に示すステップＳ２０〜２６は、図２に示すステップＳ１０〜Ｓ１６と同様の工程であるため、その詳細な説明は省略する。第２の例では、ステップＳ２８のパレート解の縮約工程から説明する。

【0044】

ここで、設計変数の定義域を設定する場合、第１の例では、図６（ａ）に示すように設計変数の定義域Ｘｉに対して下限値α〜上限値βまでの間を連続として設定する。すなわち、設計変数の定義域は、α≦Ｘｉ≦βで連続である。
これに対して、図６（ｂ）に示すように、定義域α〜βに対して、例えば、等間隔に５つの離散値を設定する。パレート解のうち、この離散値に該当するものを演算部２４にて抽出して、パレート解を縮約し、縮約されたパレート解はメモリ２８に記憶される。

【0045】

次に、縮約されたパレート解は、表示制御部３０で表示形態が変更されてパレート解を表示部１６に表示する。なお、パレート解の表示方法は、上述のように図３、図４に示すように表示することができる。図７（ａ）は、ステップＳ２６で得られたパレート解を示す散布図であり、パレート解を縮約することにより、図７（ｂ）に示すようなパレート解の散布図が得られる。
これ以外にも、図８（ａ）は、ステップＳ２６で得られたパレート解を示す散布図であり、パレート解を縮約することにより、図８（ｂ）に示すパレート解の分布を示す散布図が得られる。

【0046】

このように、パレート解を離散的に縮約することにより効果的にパレート解のデータ数を軽減することができるため、散布図において目的関数と設計変数との因果関係を見出しやすくなる。なお、離散値は、等間隔にすると変化率を見積もりやすくなるため更に好ましい。
なお、離散値については、図６（ｂ）に示す離散値に限定されるものではなく、図６（ｃ）に示すように、各離散値が１つの値ではなく、各離散値が所定の範囲ｔを有するものとしてもよい。

【0047】

表示方法の第２の例は、パレート解を算出した後に縮約したが、これに限定されるものではなく、縮約したパレート解を算出するようにしてもよい。例えば、表示方法の第１の例において、設計変数の定義域を設定する際（ステップＳ１２）、図６（ｂ）、（ｃ）に示す離散値として、特性値の算出（ステップＳ１４）、パレート解の算出（ステップＳ１６）、およびパレート解の表示（ステップＳ１８）をすることができる。この場合、第２の例の同様の図７（ｂ）に示すようなパレート解の散布図を得ることができる。

【0048】

さらには、表示方法として、設計変数毎に、特性がよいとされる領域にあるパレート解のうち、最も特性がよいパレート解のシンボルをそれぞれ線で結んで表示することもできる。
図９（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解を示す散布図であり、（ｂ）は、図９（ａ）に示すパレート解にパレートフロントを明確にした例を示す散布図である。図９（ａ）は、図８（ｂ）と同じであり、特性が良好な領域は領域Ｄである。
図９（ａ）において、設計変数毎のパレート解のうち、領域Ｄに最も近いシンボルをそれぞれ結ぶことにより、図９（ｂ）に示すように、ラインＥ_１〜Ｅ_４が得られる。このラインＥ_１〜Ｅ_４は、設計変数の値毎のパレートフロントを示すものである。これにより、設計変数の値の変化によるパレートフロントの変化を明確にすることができる。

【0049】

なお、パレート解探索部２６で得られたパレート解について、パレートフロントを求める。パレートフロントを求めるための公知の方法を用いることができる。パレートフロントの情報をメモリ２８に記憶させる。
次に、表示制御部３０において、パレートフロントの情報をメモリ２８から読み出し、この情報に基づいて、設計変数の値毎に、パレートフロントにあるパレート解を表わすシンボルを結ぶ線の情報を作成する。その線の情報をメモリ２８に記憶させる。
表示制御部３０により表示部１６に、パレート解とともに線の情報に基づいて、ラインＥ_１〜Ｅ_４を表示させる。

【0050】

また、パレート解探索部２６で得られたパレート解について、設計変数の値毎に、図９（ａ）に示すように目的関数（特性値）Ｙ２の各値において目的関数（特性値）Ｙ１が最も大きな値のパレート解を抽出する。抽出したパレート解の情報をメモリ２８に記憶させる。
そして、表示制御部３０において、抽出したパレート解を表わすシンボルを結ぶ線の情報を作成する。その線の情報をメモリ２８に記憶させる。
そして、表示制御部３０により、表示部１６に、パレート解とともに線の情報に基づいて、ラインＥ_１〜Ｅ_４を表示させる。このようにしても図９（ｂ）に示すように、パレートフロントを示すラインＥ_１〜Ｅ_４を得ることができる。

【0051】

次に、表示方法の第３の例について説明する。
図１０は、本発明の実施形態の表示方法の第３の例を工程順に示すフローチャートである。図１１（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解を示す散布図であり、（ｂ）は、図１１（ａ）に示すパレート解にパレート解探索過程の特性値データを追加した例を示す散布図である。

【0052】

表示方法の第３の例は、図２に示す表示方法の第１の例に比して、表示データ数が所定の数未満であるとき、パレート解探索過程で得られる特性値を目的関数データとして保持おき、この目的関数データをパレート解とともに表示する点以外は、第１の例と同様の工程であるため、その詳細な説明は省略する。

【0053】

第３の例の表示方法において、図１０に示すステップＳ３０〜３６は、図２に示すステップＳ１０〜Ｓ１６と同様の工程であるため、その詳細な説明は省略する。第３の例では、ステップＳ３８のパレート解の数の判定工程から説明する。
演算部２４で特性値を算出した後、メモリ２８に記憶させる。その後、パレート解探索部２６でパレート解を算出した後、パレート解探索部２６からパレート解が移動されてメモリ２８に記憶されるとともに、パレート解探索部２６でパレート解の数をカウントする。このパレート解の数が制御部３２に出力され、制御部３２において表示データ数が所定の数を超えていれば（ステップＳ３８）、パレート解が表示部１６表示される（ステップＳ４０）。この場合、例えば、図１１（ａ）に示すように表示される。なお、表示方法としては、上述の表示方法を用いることができる。

【0054】

一方、ステップＳ３８において、表示データ数が所定の数未満であるとき、メモリ２８からステップＳ３４で算出された特性値のうち、パレート解ではない特性値が目的関数データとして、パレート解とともに呼び出され、パレート解に目的関数データが追加される（ステップＳ４２）。
そして、表示部１６にパレート解とともに特性値が目的関数データとして表示される。この場合、例えば、図１１（ｂ）に示すように表示される。なお、表示方法としては、上述の表示方法を用いることができる。
このように、パレート解とパレート解以外のもの、例えば、パレート解検索過程でパレート解ではなかった特性値を目的関数データとして同時に表示することにより、パレート解とパレート解近傍のデータとの関係を可視化することができる。図１１（ｂ）に示す表示方法は、特に、特性値が２つの場合に有効である。

【0055】

次に、表示方法の第４の例について説明する。
図１２は、本発明の実施形態の表示方法の第４の例を工程順に示すフローチャートである。図１３（ａ）は、タイヤの設計変数とタイヤの特性値とについて得られたパレート解を示す散布図であり、（ｂ）は、図１３（ａ）に示すパレート解に、他の特性値を考慮したパレート解を追加した例を示す散布図である。
表示方法の第４の例は、図２に示す表示方法の第１の例に比して、表示データ数が所定の数未満であるとき、特性値の種類を追加したパレート解を別途算出し、この別途算出したパレート解（以下、拡張パレート解という）と、既に算出されたパレート解をまとめて表示する点以外は、第１の例と同様の工程であるため、その詳細な説明は省略する。

【0056】

第４の例の表示方法において、図１２に示すステップＳ３０〜３６は、図２に示すステップＳ１０〜Ｓ１６と同様の工程であるため、その詳細な説明は省略する。第３の例では、ステップＳ３８のパレート解の数の判定工程から説明する。
演算部２４で特性値を算出した後、メモリ２８に記憶させる。その後、パレート解探索部２６でパレート解を算出した後、パレート解探索部２６からパレート解が移動されてメモリ２８に記憶されるとともに、パレート解探索部２６でパレート解の数をカウントする。このパレート解の数が制御部３２に出力され、制御部３２において表示データ数が所定の数を超えていれば（ステップＳ３８）、パレート解が表示部１６表示される（ステップＳ４０）。この場合、例えば、図１３（ａ）に示すように表示される。なお、表示方法としては、上述の表示方法を用いることができる。

【0057】

一方、ステップＳ３８において、表示データ数が所定の数未満であるとき、条件設定部２０に設定されている特性値として定められたパラメータを追加し、設計変数と追加された特性値との非線形応答関係、おおび設計変数の定義域を定める。そして、モデル生成部２２でモデル生成を行い、演算部２４でパラメータ数の異なる特性値を算出する。そして、パレート解探索部２６で拡張パレート解を算出し、メモリ２８に記憶させる。具体的には、現時点で設定された特性値、例えば、ころがり抵抗と横バネ定数に、更に別の特性値、例えば、トレッド部の摩耗、ロードノイズを加えて、タイヤの形状パラメータ（設計変数）とこれらの４つの特性値との非線形応答関係およびタイヤの形状パラメータの定義域を定める。そして、モデル生成部２２にてモデル生成がされ、演算部２４にて４つの特性値が算出される。そして、パレート解探索部２６で４つの特性値を目的関数とする最適化を実施し、拡張パレート解を算出する（ステップＳ４４）。拡張パレート解をメモリ２８に記憶させる。そして、ステップＳ３８において、表示データ数が所定の数を超えるまで、繰り返し拡張パレート解の算出がなされる（ステップＳ４４）。なお、追加する特性値として定められたパラメータ値は、特に限定されるものではない。

【0058】

ステップＳ３８において、表示データ数が所定の数を超えると、パレート解とともに拡張パレート解がメモリ２８から呼び出され、表示部１６にパレート解とともに拡張パレート解を表示させる。この場合、例えば、図１３（ｂ）に示すように表示される。なお、表示方法としては、上述の表示方法を用いることができる。
このように、パレート解と、特性値の数が異なる拡張パレート解を同時に表示することにより、パレート解とパレート解近傍のデータとの関係を可視化することができる。特に、特性値が２つの場合に有効である。

【0059】

以上のように、タイヤの特性値とタイヤの形状パラメータ（設計変数）との関係を可視化することで、その因果関係を理解でき、得られる設計情報を設計へ生かすことができる。特に設計の初期段階の方向性を定める際にこの情報を用いることで、設計終盤の詳細設計においても大きく修正することのない商品開発を進めることができる。それにより開発コストを削減できるとともに商品のリードタイムを短縮することができる。
本実施形態では、タイヤを例にして説明したが、本発明の表示方法は、これに限定されるものではない。例えば、ゴム製品、家電製品、自動車、および飛行機等の構造設計にも適用することができる。この場合でも、設計変数と特性値（目的関数）との因果関係を理解でき、得られる設計情報を設計へ生かすことができる。更には上述のように、商品開発を進めることができ、開発コストを削減できるとともに開発のリードタイムを短縮することができる。

【0060】

本発明は、基本的に以上のように構成されるものである。以上、本発明のデータの表示方法について詳細に説明したが、本発明は上記実施形態に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良または変更をしてもよいのはもちろんである。

【符号の説明】

【0061】

１０ 表示処理装置
１２ 処理部
１４ 入力部
１６ 表示部
２０ 条件設定部
２２ モデル生成部
２４ 演算部
２６ パレート解探索部
２８ メモリ
３０ 表示制御部
３２ 制御部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】