潤滑基礎方程式を用いて潤滑計算を行うシミュレーション装置であって、エンジンのシリンダ内面に形成される微細形状を表す形状パラメータとして、少なくとも前記微細形状のプラトー率を入力する入力手段と、前記入力手段によって入力される前記形状パラメータに基づいて、前記潤滑基礎方程式に前記微細形状が油膜形成に及ぼす影響を導入するための係数である修正係数を算出する修正係数算出手段と、を具備することを特徴とするシミュレーション装置。前記修正係数算出手段によって算出される前記修正係数を導入した前記潤滑基礎方程式を、前記エンジンのピストンリング列と、前記微細形状を有する前記シリンダ内面との間の潤滑膜に適用し、潤滑計算を行う潤滑計算手段、を更に具備することを特徴とする請求項1に記載のシミュレーション装置。前記修正係数算出手段は、前記シリンダ内面を前記微細形状と仮定し、前記入力手段によって入力される前記形状パラメータに基づいて、前記ピストンリング列と前記シリンダ内面との間の油膜の流れの数値計算を行い、圧力流れ成分の流量、及び、せん断流れ成分の流量を算出する第1数値計算手段と、前記シリンダ内面を前記平滑面と仮定し、前記ピストンリング列と前記シリンダ内面との間の油膜の流れの数値計算を行い、圧力流れ成分の流量、及び、せん断流れ成分の流量を算出する第2数値計算手段と、前記第1数値計算手段によって算出される圧力流れ成分の流量と、前記第2数値計算手段によって算出される圧力流れ成分の流量との比を圧力流量係数として算出し、前記第1数値計算手段によって算出されるせん断流れ成分の流量と、前記第2数値計算手段によって算出されるせん断流れ成分の流量との比をせん断流量係数として算出する係数算出手段と、を含み、前記潤滑計算手段は、前記係数算出手段によって算出される前記圧力流量係数及び前記せん断流量係数を前記修正係数として導入した前記潤滑基礎方程式を用いて潤滑計算を行い、前記ピストンリング列と前記シリンダ内面との間の膜厚分布を算出することを特徴とする請求項2に記載のシミュレーション装置。前記第1数値計算手段は、更に、前記シリンダ内面を前記微細形状と仮定して、圧力流れ成分のせん断応力、及び、せん断流れ成分のせん断応力を算出し、前記第2数値計算手段は、更に、前記シリンダ内面を前記平滑面と仮定して、圧力流れ成分のせん断応力、及び、せん断流れ成分のせん断応力を算出し、前記係数算出手段は、前記第1数値計算手段によって算出される前記圧力流れ成分のせん断応力と、前記第2数値計算手段によって算出される前記圧力流れ成分のせん断応力との比を圧力流せん断応力係数、前記第1数値計算手段によって算出される前記せん断流れ成分のせん断応力と、前記第2数値計算手段によって算出される前記せん断流れ成分のせん断応力との比をせん断流せん断応力係数として算出し、前記圧力流せん断応力係数、前記せん断流せん断応力係数、及び前記膜厚分布に基づいて、前記ピストンリング列と前記シリンダ内面との間の摩擦力を算出する摩擦力算出手段、を更に具備することを特徴とする請求項3に記載のシミュレーション装置。前記微細形状は、ハッチ形状であり、前記入力手段は、更に、前記形状パラメータとして、前記微細形状のハッチ角度、及び溝深さを入力し、前記修正係数算出手段は、前記入力手段によって入力される前記プラトー率、前記ハッチ角度、及び前記溝深さに基づいて、前記修正係数を算出することを特徴とする請求項1乃至請求項4のいずれかに記載のシミュレーション装置。コンピュータが潤滑基礎方程式を用いて潤滑計算を行うシミュレーション方法であって、エンジンのシリンダ内面に形成される微細形状について、プラトー率を含む前記微細形状を表す形状パラメータとして、少なくとも前記微細形状のプラトー率を入力する入力ステップと、前記入力ステップによって入力される前記形状パラメータに基づいて、前記潤滑基礎方程式に前記微細形状が油膜形成に及ぼす影響を導入するための係数である修正係数を算出する修正係数算出ステップと、を含むことを特徴とするシミュレーション方法。コンピュータを、エンジンのシリンダ内面に形成される微細形状を表す形状パラメータとして、少なくとも前記微細形状のプラトー率を入力する入力手段と、前記入力ステップによって入力される前記形状パラメータに基づいて、潤滑計算に用いられる潤滑基礎方程式に前記微細形状が油膜形成に及ぼす影響を導入するための係数である修正係数を算出する修正係数算出手段と、して機能させるためのシミュレーションプログラム。従来技術1では、上記2パラメータが等しい2面は、全く等価なものとして扱われる。しかしながら、実際には、対象表面形状実測値から統計的手法により2つのパラメータを求める必要がある。このとき、例えばシリンダ内面に形成される微細形状がハッチ形状の場合、ハッチ角度は方向性パラメータγ、溝深さはRMSあらさσに反映されるが、プラトー率溝密度を反映させるパラメータがない。また、従来技術1では、ハッチ形状がクロスハッチ両ハッチの場合と片ハッチの場合を区別することもできない。以上のように、従来技術1では、ハッチ形状の情報を十分に計算に反映させることができず、シリンダ内面に形成される最適なハッチ形状を求めるには不十分である。従来技術2では、任意のハッチ形状の解を正確に得ることが出来るが、計算負荷が高く、実際のピストンリングとシリンダ内面との間の流れを解いて摩擦特性を求めるには膨大な計算時間を要する。従って、従来技術2は、シリンダ内面に形成される微細形状の設計ツールとしては不向きである。本発明は、前述した問題点に鑑みてなされたもので、その目的とすることは、シリンダ内面に形成される微細形状が油膜形成に及ぼす影響を短時間かつ精度良く算出することが可能なシミュレーション装置等を提供することである。前述した目的を達成するために第1の発明は、潤滑基礎方程式を用いて潤滑計算を行うシミュレーション装置であって、エンジンのシリンダ内面に形成される微細形状を表す形状パラメータとして、少なくとも前記微細形状のプラトー率を入力する入力手段と、前記入力手段によって入力される前記形状パラメータに基づいて、前記潤滑基礎方程式に前記微細形状が油膜形成に及ぼす影響を導入するための係数である修正係数を算出する修正係数算出手段と、を具備することを特徴とするシミュレーション装置である。第1の発明によって、シリンダ内面に形成される微細形状が油膜形成に及ぼす影響を短時間かつ精度良く算出することができる。第1の発明は、前記修正係数算出手段によって算出される前記修正係数を導入した前記潤滑基礎方程式を、前記エンジンのピストンリング列と、前記微細形状を有する前記シリンダ内面との間の潤滑膜に適用し、潤滑計算を行う潤滑計算手段、を更に具備することが望ましい。これによって、シリンダ内面に形成される微細形状が及ぼす油膜厚さへの影響を理論的に推定することが可能となる。例えば、第1の発明における前記修正係数算出手段は、前記シリンダ内面を前記微細形状と仮定し、前記入力手段によって入力される前記形状パラメータに基づいて、前記ピストンリング列と前記シリンダ内面との間の油膜の流れの数値計算を行い、圧力流れ成分の流量、及び、せん断流れ成分の流量を算出する第1数値計算手段と、前記シリンダ内面を前記平滑面と仮定し、前記ピストンリング列と前記シリンダ内面との間の油膜の流れの数値計算を行い、圧力流れ成分の流量、及び、せん断流れ成分の流量を算出する第2数値計算手段と、前記第1数値計算手段によって算出される圧力流れ成分の流量と、前記第2数値計算手段によって算出される圧力流れ成分の流量との比を圧力流量係数として算出し、前記第1数値計算手段によって算出されるせん断流れ成分の流量と、前記第2数値計算手段によって算出されるせん断流れ成分の流量との比をせん断流量係数として算出する係数算出手段と、を含み、前記潤滑計算手段は、前記係数算出手段によって算出される前記圧力流量係数及び前記せん断流量係数を前記修正係数として導入した前記潤滑基礎方程式を用いて潤滑計算を行い、前記ピストンリング列と前記シリンダ内面との間の膜厚分布を算出する。これによって、ピストンリング列と、所定の微細形状が形成されるシリンダ内面との間の膜厚分布を短時間かつ精度良く推定することができる。また、例えば、前記第1数値計算手段は、更に、前記シリンダ内面を前記微細形状と仮定して、圧力流れ成分のせん断応力、及び、せん断流れ成分のせん断応力を算出し、前記第2数値計算手段は、更に、前記シリンダ内面を前記平滑面と仮定して、圧力流れ成分のせん断応力、及び、せん断流れ成分のせん断応力を算出し、前記係数算出手段は、前記第1数値計算手段によって算出される前記圧力流れ成分のせん断応力と、前記第2数値計算手段によって算出される前記圧力流れ成分のせん断応力との比を圧力流せん断応力係数、前記第1数値計算手段によって算出される前記せん断流れ成分のせん断応力と、前記第2数値計算手段によって算出される前記せん断流れ成分のせん断応力との比をせん断流せん断応力係数として算出し、前記圧力流せん断応力係数、前記せん断流せん断応力係数、及び前記膜厚分布に基づいて、前記ピストンリング列と前記シリンダ内面との間の摩擦力を算出する摩擦力算出手段、を更に具備する。これによって、ピストンリング列と、所定の微細形状が形成されるシリンダ内面との間の摩擦力を短時間かつ精度良く推定することができる。また、例えば、第1の発明における前記微細形状は、ハッチ形状であり、前記入力手段は、更に、前記形状パラメータとして、前記微細形状のハッチ角度、及び溝深さを入力し、前記修正係数算出手段は、前記入力手段によって入力される前記プラトー率、前記ハッチ角度、及び前記溝深さに基づいて、前記修正係数を算出する。これによって、シリンダ内面に形成されるハッチ形状が油膜形成に及ぼす影響を短時間かつ精度良く算出することができる。第2の発明は、コンピュータが潤滑基礎方程式を用いて潤滑計算を行うシミュレーション方法であって、エンジンのシリンダ内面に形成される微細形状について、プラトー率を含む前記微細形状を表す形状パラメータとして、少なくとも前記微細形状のプラトー率を入力する入力ステップと、前記入力ステップによって入力される前記形状パラメータに基づいて、前記潤滑基礎方程式に前記微細形状が油膜形成に及ぼす影響を導入するための係数である修正係数を算出する修正係数算出ステップと、を含むことを特徴とするシミュレーション方法である。第2の発明によって、シリンダ内面に形成される微細形状が油膜形成に及ぼす影響を短時間かつ精度良く算出することができる。第3の発明は、コンピュータを、エンジンのシリンダ内面に形成される微細形状を表す形状パラメータとして、少なくとも前記微細形状のプラトー率を入力する入力手段と、前記入力ステップによって入力される前記形状パラメータに基づいて、潤滑計算に用いられる潤滑基礎方程式に前記微細形状が油膜形成に及ぼす影響を導入するための係数である修正係数を算出する修正係数算出手段と、して機能させるためのシミュレーションプログラムである。第3の発明を汎用のコンピュータにインストールすることによって、第1の発明のシミュレーション装置を得ることができる。本発明により、シリンダ内面に形成される微細形状が油膜形成に及ぼす影響を短時間かつ精度良く算出することが可能なシミュレーション装置等を提供することができる。以下図面に基づいて、本発明の実施形態を詳細に説明する。図1は、シミュレーション装置1を実現するコンピュータのハードウエア構成図である。尚、図1のハードウエア構成は一例であり、用途、目的に応じて様々な構成を採ることが可能である。シミュレーション装置1は、制御部11、記憶部12、メディア入出力部13、通信制御部14、入力部15、表示部16、周辺機器I/F部17等が、バス18を介して接続される。制御部11は、CPUCentralProcessingUnit、ROMReadOnlyMemory、RAMRandomAccessMemory等で構成される。CPUは、記憶部12、ROM、記録媒体等に格納されるプログラムをRAM上のワークメモリ領域に呼び出して実行し、バス18を介して接続された各装置を駆動制御し、シミュレーション装置1が行う後述する処理を実現する。ROMは、不揮発性メモリであり、コンピュータのブートプログラムやBIOS等のプログラム、データ等を恒久的に保持している。RAMは、揮発性メモリであり、記憶部12、ROM、記録媒体等からロードしたプログラム、データ等を一時的に保持するとともに、制御部11が各種処理を行う為に使用するワークエリアを備える。記憶部12は、HDDハードディスクドライブであり、制御部11が実行するプログラム、プログラム実行に必要なデータ、OSオペレーティングシステム等が格納される。プログラムに関しては、OSオペレーティングシステムに相当する制御プログラムや、後述する処理をコンピュータに実行させるためのアプリケーションプログラムが格納されている。これらの各プログラムコードは、制御部11により必要に応じて読み出されてRAMに移され、CPUに読み出されて各種の手段として実行される。メディア入出力部13ドライブ装置は、データの入出力を行い、例えば、CDドライブ-ROM、-R、-RW等、DVDドライブ-ROM、-R、-RW等等のメディア入出力装置を有する。通信制御部14は、通信制御装置、通信ポート等を有し、コンピュータとネットワーク間の通信を媒介する通信インタフェースであり、ネットワークを介して、他のコンピュータ間との通信制御を行う。ネットワークは、有線、無線を問わない。入力部15は、データの入力を行い、例えば、キーボード、マウス等のポインティングデバイス、テンキー等の入力装置を有する。入力部15を介して、コンピュータに対して、操作指示、動作指示、データ入力等を行うことができる。表示部16は、液晶パネル等のディスプレイ装置、ディスプレイ装置と連携してコンピュータのビデオ機能を実現するための論理回路等ビデオアダプタ等を有する。周辺機器I/Fインタフェース部17は、コンピュータに周辺機器を接続させるためのポートであり、周辺機器I/F部17を介してコンピュータは周辺機器とのデータの送受信を行う。周辺機器I/F部17は、USBやIEEE1394やRS-232C等で構成されており、通常複数の周辺機器I/Fを有する。周辺機器との接続形態は有線、無線を問わない。バス18は、各装置間の制御信号、データ信号等の授受を媒介する経路である。シミュレーション装置1は、潤滑基礎方程式レイノルズ方程式を用いて、エンジンのピストンとシリンダ内面との間の潤滑計算を行う装置である。エンジンのシリンダ内面には、摩擦及びオイル消費の減少などを目的として、細かな凹凸の繰返しである微細形状を形成させることがある。以下では、図2から図12を参照しながら、微細形状として、ハッチ形状を例にして説明する。ハッチ形状は、クロスハッチ両ハッチと片ハッチの2種類がある。図2は、シミュレーション装置1の処理の流れを示すフローチャートである。図2に示すように、シミュレーション装置1の制御部11は、入力手段メディア入出力部13、通信制御部14、入力部15、周辺機器I/F部17等を介して、形状パラメータを入力するS1。また、制御部11は、記憶部12にファイルとして記憶されている形状パラメータを入力しても良い。ハッチ形状の場合、形状パラメータは、ハッチ角度、溝深さ、プラトー率及びプラトー面あらさである。図3は、クロスハッチ面と形状パラメータを説明する図である。図3aは、シリンダ内面に形成されるクロスハッチ面と油膜の模式図である。クロスハッチ面が形成されることによって、シリンダ内面は無数の溝を有する。シリンダ内面に対するピストン不図示のすべり方向、すなわちピストンの移動方向は、紙面に対して左右方向である。図3bは、クロスハッチ面を正面方向から見たときの模式図と、この模式図の一部の拡大断面図である。図3に示すように、油膜の厚さは、平均膜厚hとして表現される。ハッチ形状は、ハッチ角度、溝深さ、プラトー率及びプラトー面あらさの形状パラメータによって定義される。ここで、プラトー台地領域とは、溝ではない領域を意味する。ハッチ角度とは、例えば、すべり方向と直交する方向に伸びる直線を基準線とし、基準線と溝とのなす角度である。尚、基準線は、すべり方向と平行する方向に伸びる直線としても良い。溝深さaは、プラトーから最深部までの長さである。プラトー率は、隣接する溝の最深部同士の距離をL、プラトーの幅をWとしたとき、プラトー率=W/Lとして定義される。プラトー面あらさは、プラトー台地領域の面あらさである。面あらさの定義は、最大高さ、算術平均あらさ、自乗平均平方根あらさRMSなどがあり、本発明では、例えば、自乗平均平方根あらさRMSを用いる。図3は、両ハッチの例であるが、片ハッチの場合も同様に、ハッチ角度、溝深さ、プラトー率及びプラトー面あらさの形状パラメータが定義される。制御部11が実行するプログラムには、微細形状の関数fとして、ハッチ角度、溝深さ、及びプラトー率の形状パラメータを変数とする関数が定義されている。関数fは、ピストンに嵌め込まれるピストンリング列不図示とシリンダ内面との間の膜厚を計算するための関数である。例えば、片ハッチにおける膜厚hshsのsはsingleの頭文字は、ハッチ角度をθ、溝深さをa、プラトー率をqとすると、hs=fθ、a、qとして計算される。また、例えば、両ハッチにおける膜厚hdhdのdはdoubleの頭文字は、第1の方向の片ハッチについて、ハッチ角度をθ1、溝深さをa1、プラトー率をq1とし、第1の方向の片ハッチと交差する第2の方向の片ハッチについて、ハッチ角度をθ2、溝深さをa2、プラトー率をq2とすると、hd=minfθ1、a1、q1、fθ2、a2、q2として計算される。図2の説明に戻る。潤滑膜は流れすべり方向に対して厚み方向が極めて小さい薄膜である。この流れは、通常圧力流れPoiseuilleFlowとせん断流れCouetteFlowに分離し、これらの和とみなすことができる。この圧力流れ、せん断流れそれぞれについて、微細形状が及ぼす影響を数値計算によって定量化する。そこで、制御部11は、シリンダ内面を微細形状と仮定し、S11において入力される形状パラメータ微細形状のハッチ角度、溝深さ、プラトー率及びプラトー面あらさに基づいて、ピストンリング列とシリンダ内面との間の油膜の流れについて数値計算を行い、圧力流れ成分の流量及びせん断応力、並びに、せん断流れ成分の流量及びせん断応力を算出するS2:第1数値計算。また、制御部11は、シリンダ内面を平滑面と仮定し、ピストンリング列とシリンダ内面との間の油膜の流れについて数値計算を行い、圧力流れ成分の流量及びせん断応力、並びに、せん断流れ成分の流量及びせん断応力を算出するS3:第2数値計算。そして、制御部11は、数値計算の結果に基づいて、修正係数を算出するS4。より詳しくは、制御部11は、以下の4つの修正係数を算出する。1圧力流量係数:S2において算出される圧力流れ成分の流量と、S3において算出される圧力流れ成分の流量との比2せん断流量係数:S2において算出されるせん断流れ成分の流量と、S3において算出されるせん断流れ成分の流量との比3圧力流せん断応力係数:S2において算出される前記圧力流れ成分のせん断応力と、S3において算出される圧力流れ成分のせん断応力との比4せん断流せん断応力係数:S2において算出されるせん断流れ成分のせん断応力と、S3において算出されるせん断流れ成分のせん断応力との比ここで、修正係数とは、潤滑基礎方程式レイノルズ方程式に、微細形状が油膜形成に及ぼす影響を導入するための係数である。制御部11が実行するプログラムには、各項に修正係数が含まれる潤滑基礎方程式レイノルズ方程式が定義されている。図4は、圧力流れ成分の計算パラメータについて説明する図である。図4に示すように、油膜長さLx、油膜幅Lyの仮想油膜を考える。油膜の厚さは、油膜長さLx、油膜幅Lyに対して十分小さく、薄膜近似ができるものとする。また、油膜長さLx×油膜幅Ly内の領域には、ハッチ形状の凹凸が十分な数だけ存在するとする。そして、油膜の左側面x=0における境界条件を圧力p=p0、油膜の右側面x=Lxにおける境界条件を圧力p=0とすると、圧力差p0に起因する流れ圧力流れ、Poiseuille流れが油膜内に生じる。Poiseuille流れの流量Qxは、次式の潤滑基礎方程式レイノルズ方程式を解くことによって求まる。尚、hは流体膜厚さ、ηは粘度、pは圧力である。図5は、圧力流量係数の計算結果の一例を示す図である。式1を解いて得られる流量は、下面シリンダ内面が微細形状の場合と、平滑面の場合とで異なる。制御部11は、平滑面の場合の流量に対する微細形状の場合の流量との比を、圧力流量係数φxと定義し、数値計算の結果から算出する。すなわち、制御部11は、φx=Qx{下面が微細形状の場合}/Qx{下面が平滑面の場合}を算出する。図5aには、圧力流量係数φxの計算結果の一例が示されている。また、図5bには、図5aに示す計算結果の物理的意味が示されている。図5aの横軸は油膜厚さ/溝深さであり、縦軸は圧力流量係数φxである。図5aには、3種の異なるハッチ形状A、B、Cの例が示されている。3つの例のそれぞれについて、油膜厚さ/溝深さが異なる数回の数値計算の結果に基づき、カーブフィットによって圧力流量係数φxの近似式を求める。図5aに示されている3つの曲線は、カーブフィットによって求めた近似式によって描画されている。図5aにおける3つの曲線のうち、実線数値計算結果を正方形にて図示がハッチ形状A、点線数値計算結果を菱形にて図示がハッチ形状B、一点鎖線数値計算結果を三角形にて図示がハッチ形状Cを示している。図5bに示すように、ハッチ形状によって下面シリンダ内面近傍の流動抵抗が増大し、φx分だけ流量が減少しようとする。これに逆らって同一流量を流そうとすると油膜圧力が発生し、平均膜厚を増大させる方向に作用する。図6は、圧力流せん断応力係数の計算結果の一例を示す図である。式1を解いて得られる下面シリンダ内面上の速度勾配から、下面シリンダ内面に作用するせん断応力が決まる。制御部11は、平滑面の場合のせん断応力に対する微細形状の場合のせん断応力との比を、圧力流せん断応力係数φfと定義し、数値計算の結果から算出する。すなわち、制御部11は、φf=dudy{下面が微細形状の場合}/dudy{下面が平滑面の場合}を算出する。図6aには、圧力流せん断応力係数φfの計算結果の一例が示されている。また、図6bには、図6aに示す計算結果の物理的意味が示されている。図6aの横軸は油膜厚さ/溝深さであり、縦軸は圧力流せん断応力係数φfである。図6aには、3種の異なるハッチ形状A、B、Cの例が示されている。3つの例のそれぞれについて、油膜厚さ/溝深さが異なる数回の数値計算の結果に基づき、カーブフィットによって圧力流せん断応力係数φfの近似式を求める。図6aに示されている3つの曲線は、カーブフィットによって求めた近似式によって描画されている。図6aにおける3つの曲線のうち、実線数値計算結果を正方形にて図示がハッチ形状A、点線数値計算結果を菱形にて図示がハッチ形状B、一点鎖線数値計算結果を三角形にて図示がハッチ形状Cを示している。図6bに示すように、ハッチ形状によって下面シリンダ内面近傍の速度勾配が減少し、φf分だけせん断応力粘性摩擦が減少する。図7は、せん断流れ成分の計算パラメータについて説明する図である。圧力流れ成分の計算と同様に、図7に示すように、油膜長さLx、油膜幅Lyの仮想油膜を考える。油膜の厚さhは、油膜長さLx、油膜幅Lyに対して十分小さく、薄膜近似ができるものとする。また、油膜長さLx×油膜幅Ly内の領域には、ハッチ形状の凹凸が十分な数だけ存在するとする。また、上面をxの正方向の速度Uによって平行に滑らせるとする。そして、油膜の左側面x=0における境界条件を圧力p=0、油膜の右側面x=Lxにおける境界条件を圧力p=0とすると圧力差を0とすると、上面のすべりUに起因する流れせん断流れ、Couette流れが油膜内に生じる。Couette流れの流量Qsは、次式の潤滑基礎方程式レイノルズ方程式を解くことによって求まる。尚、hは流体膜厚さ、ηは粘度、pは圧力、Uは速度である。図8は、せん断流量係数の計算結果の一例を示す図である。式2を解いて得られる流量は、下面シリンダ内面が微細形状の場合と、平滑面の場合とで異なる。制御部11は、平滑面の場合の流量に対する微細形状の場合の流量との比を、せん断流量係数φsと定義し、数値計算の結果から算出する。すなわち、制御部11は、φs=Qs{下面が微細形状の場合}/Qs{下面が平滑面の場合}を算出する。図8aには、せん断流量係数φsの計算結果の一例が示されている。また、図8bには、図8aに示す計算結果の物理的意味が示されている。図8aの横軸は油膜厚さ/溝深さであり、縦軸はせん断流量係数φsである。図8aには、3種の異なるハッチ形状A、B、Cの例が示されている。3つの例のそれぞれについて、油膜厚さ/溝深さが異なる数回の数値計算の結果に基づき、カーブフィットによってせん断流量係数φsの近似式を求める。図8aに示されている3つの曲線は、カーブフィットによって求めた近似式によって描画されている。図8aにおける3つの曲線のうち、実線数値計算結果を正方形にて図示がハッチ形状A、点線数値計算結果を菱形にて図示がハッチ形状B、一点鎖線数値計算結果を三角形にて図示がハッチ形状Cを示している。図8bに示すように、ハッチ形状によって下面シリンダ内面近傍の流動抵抗が増大し、φs分だけ流量が減少しようとする。これに逆らって同一流量を流そうとすると油膜圧力が発生し、平均膜厚を増大させる方向に作用する。図9は、せん断流せん断応力係数の計算結果の一例を示す図である。式2を解いて得られる下面シリンダ内面上の速度勾配から、下面シリンダ内面に作用するせん断応力が決まる。制御部11は、平滑面の場合のせん断応力に対する微細形状の場合のせん断応力との比を、せん断流せん断応力係数φfsと定義し、数値計算の結果から算出する。すなわち、制御部11は、φfs=dudy{下面が微細形状の場合}/dudy{下面が平滑面の場合}を算出する。図9aにはせん断流せん断応力係数φfsの計算結果の一例が示されている。また、図9bには、図9aに示す計算結果の物理的意味が示されている。図9aの横軸は油膜厚さ/溝深さであり、縦軸はせん断流せん断応力係数φfsである。図9aには、3種の異なるハッチ形状A、B、Cの例が示されている。3つの例のそれぞれについて、油膜厚さ/溝深さが異なる数回の数値計算の結果に基づき、カーブフィットによってせん断流せん断応力係数φfsの近似式を求める。図9aに示されている3つの曲線は、カーブフィットによって求めた近似式によって描画されている。図9aにおける3つの曲線のうち、実線数値計算結果を正方形にて図示がハッチ形状A、点線数値計算結果を菱形にて図示がハッチ形状B、一点鎖線数値計算結果を三角形にて図示がハッチ形状Cを示している。図9bに示すように、ハッチ形状によって下面シリンダ内面近傍の速度勾配が増加し、φfs分だけせん断応力粘性摩擦が増加する。図2の説明に戻る。次に、制御部11は、S4において算出される修正係数を導入した潤滑基礎方程式レイノルズ方程式を、ピストンリング列と、微細形状を有するシリンダ内面との間の潤滑膜に適用し、潤滑計算を行い、ピストンリング列とシリンダ内面との間の圧力分布や膜厚分布を算出するS5。更に、制御部11は、数値計算及び潤滑計算の結果に基づいて、摩擦力を算出するS6。より詳しくは、制御部11は、圧力流せん断応力係数、せん断流せん断応力係数、及び膜厚分布等に基づいて、ピストンリング列とシリンダ内面との間の摩擦力を算出する。そして、制御部11は、出力手段メディア入出力部13、通信制御部14、表示部16、周辺機器I/F部17等を介して、S5及びS6の計算結果を出力するS7。また、制御部11は、記憶部12に、ファイルとしてS5及びS6の計算結果を出力しても良い。図10は、潤滑計算及び摩擦計算を説明する図である。図10aは、ピストンとシリンダ内面の一部の模式図である。図10bは、図10aの一部を拡大した図であり、ピストンに嵌め込まれている3つのピストンリング、シリンダ内面、及び油膜の模式図である。図10cは、図10bの一部を拡大した図であり、1つのピストンリング、シリンダ内面、及び油膜の模式図である。図10bに示すTopは、ピストンリングの1つであるトップリング、2ndは、ピストンリングの1つであるセカンドリング、Oilは、ピストンリングの1つであるオイルリングを示している。トップリングやセカンドリングは、燃焼室で発生した燃焼ガスを外に出さないように、燃焼室内の機密性を保つ役割を果たしている。オイルリングは、シリンダ内面の潤滑油をかき落とし、オイルの使い過ぎを防ぐ役割を果たしている。図10bに示すように、シリンダ内面と対向する面の形状は、ピストンリングごとに異なる。本発明は、図10に示すピストンリングの例に限定されるわけではなく、例えば、2ピース構造のリングなどにも同様に適用可能である。図10cに示すように、制御部11は、ピストンリングの面をすべり方向x軸方向に速度Uによって平行に滑らせ、すべり方向と直交する方向y軸方向に荷重Wをかけると仮定する。そして、制御部11は、図5に示す圧力流量係数φxの近似式、図8に示すせん断流量係数φsの近似式を、次式に示すように、ピストンリング列の潤滑基礎方程式レイノルズ方程式に導入し、ピストンリングとシリンダ内面との間の油膜の潤滑計算を行い、膜厚分布を算出する。尚、hは流体膜厚さ、ηは粘度、pは圧力、Uは速度、tは時間である。更に、制御部11は、潤滑計算によって算出される膜厚分布に基づいて、ピストンリングとシリンダ内面が固体接触している箇所固体接触部と、油膜を介して接触している箇所とを判別し、固体接触部については公知の境界潤滑モデルによって境界摩擦の計算を行い、油膜を介して接触している箇所については粘性摩擦の計算を行う。境界摩擦の計算においては、公知の境界潤滑モデルにおいて、予め与えられるプラトー面あらさと、クロスハッチ溝によって形成される微細形状とを合成したあらさを用いる。粘性摩擦の計算については、制御部11は、潤滑計算によって算出される膜厚分布、図5に示す圧力流量係数φxの近似式、図8に示すせん断流量係数φsの近似式に基づいて、次式によって、粘性せん断応力τ摩擦力を算出する。尚、hは流体膜厚さ、pは圧力、Uは速度である。図11は、膜厚分布と摩擦力の計算結果の一例を示す図である。図11では、中速中負荷での計算結果を示している。図11aの横軸はCACrankAngle:クランク角であり、縦軸はリング列摩擦力である。図11bの横軸はCACrankAngle:クランク角であり、縦軸はTopリング油膜厚さである。1CA=-360°~-180°が吸気行程、2CA=-180°~0°が圧縮行程、3CA=0°~180°が燃焼膨張行程、4CA=180°~360°が排気行程である。図11aにおける3つの曲線のうち、実線がハッチ形状A、点線がハッチ形状B、一点鎖線がハッチ形状Cを示している。リング列摩擦力は、ハッチ形状Aが最も大きく、続いてハッチ形状B、ハッチ形状Cという結果となった。また、図11bにおいても、実線がハッチ形状A、点線がハッチ形状B、一点鎖線がハッチ形状Cを示している。Topリング油膜厚さは、ハッチ形状Aが最も薄く、ハッチ形状B及びハッチ形状Cはほとんど変わらず、ハッチ形状Aよりも厚いという結果となった。図12は、ハッチ角度による摩擦平均有効圧の比較結果の一例を示す図である。図12では、中速中負荷での計算結果を示している。図12の縦軸はFMEPFrictionMeanEffectivePressure:摩擦平均有効圧である。図12では、ハッチ形状A、ハッチ形状B、ハッチ形状Cについて、それぞれTopトップリング、2ndセカンドリング、Upperオイルリングの上部、及びLowerオイルリングの下部のFMEPの合算値を棒グラフとして表現している。図12に示す計算結果からは、摩擦を減らす為には、シリンダ内面の微細形状をハッチ形状Cとすることが望ましいと推察される。以上の通り、本発明のシミュレーション装置1によって、ハッチ形状の具体的な形状パラメータハッチ角度、溝深さ、プラトー率、及びプラトー面あらさに基づいて、油膜厚さ及び摩擦力への影響を理論的に推定することが可能となる。本発明は、例えば、摩擦及びオイル消費の最小化に最も有利なシリンダ内面に形成されるハッチ形状の形状パラメータの指針提示に極めて有用である。<変形例>以上では、シリンダ内面に形成される微細形状として、特に、ハッチ形状両ハッチ及び片ハッチを含む。を例にして説明したが、本発明は、この例に限定されない。ハッチ形状の例のように、シリンダ内面に形成される微細形状について、いくつかの形状パラメータを変数とする形状関数を定義することができれば、本発明を適用することができる。例えば、微細形状として、ディンプル形状小さな無数のくぼみ領域を有する形状。ゴルフボールなどに形成されることが多い。などについても、本発明を適用することができる。くぼみ領域が縦横に規則正しく配置されるディンプル形状の場合、くぼみ領域の半径をr、隣接するくぼみ領域の最深部同士の距離をdとすると、例えば、プラトー率=プラトー領域の面積/全体の面積=全体の面積-くぼみ領域の面積/全体の面積=d^2-π×r^2/d^2のように定義することができる。但し、x^yは、xのy乗を表す。そして、その他の形状パラメータとして、くぼみ領域の最深部の深さハッチ形状の溝深さに対応やプラトー面あらさなどを定義し、これらの形状パラメータを変数とする形状関数を定義することができる。その他、レーザ加工などによって形成される様々な加工溝の形状についても、本発明を適用することができる。以上、添付図面を参照しながら、本発明に係るシミュレーション装置等の好適な実施形態について説明したが、本発明はかかる例に限定されない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本発明の技術的範囲に属するものと了解される。