搬送パラメータに従い搬送経路に沿って搬送される柔軟媒体の挙動を、モデルによりシミュレートするシミュレーション装置であって、前記搬送経路は、搬送ローラを含む複数の搬送要素により形成されており、前記モデルは、複数の質点と、各々が、前記複数の質点のうち隣接する2つの質点の間を連結する複数のバネと、を有し、前記複数の搬送要素のそれぞれを、対応する位置情報にしたがって配置することにより、前記搬送経路を作成する経路作成部と、モデル情報に基づいて前記モデルを作成するモデル作成部と、前記搬送パラメータのうち、前記モデルの各質点と前記搬送ローラの表面との間の摩擦係数を演算するパラメータ演算部と、前記モデルが前記搬送パラメータに従って前記複数の搬送要素により搬送される場合において、前記モデルの運動の時間的変化を演算する運動演算部と、を備え、前記パラメータ演算部は、前記摩擦係数を、浮上量に従って時刻毎に演算するとともに、前記浮上量は、前記モデルの各質点と前記搬送ローラの表面との間の距離として、時刻毎に演算されることを特徴とするシミュレーション装置。請求項1に記載のシミュレーション装置において、前記運動演算部により演算された前記モデルの運動の時間的変化を、前記柔軟媒体の挙動として表示させる表示部、をさらに備えることを特徴とするシミュレーション装置。搬送パラメータに従い搬送経路に沿って搬送される柔軟媒体の挙動を、モデルによりシミュレートするコンピュータ読み取り可能なプログラムであって、前記搬送経路は、搬送ローラを含む複数の搬送要素により形成されており、前記モデルは、複数の質点と、各々が、前記複数の質点のうち隣接する2つの質点の間を連結する複数のバネと、を有し、前記コンピュータによる前記プログラムの実行は、前記コンピュータに前記複数の搬送要素のそれぞれを、対応する位置情報にしたがって配置することにより、前記搬送経路を作成する機能と、モデル情報に基づいて前記モデルを作成するモデル作成部と、前記搬送パラメータのうち、前記モデルの各質点と前記搬送ローラの表面との間の摩擦係数を演算するパラメータ演算部と、前記モデルが、前記搬送パラメータに従って前記複数の搬送要素により搬送される場合において、前記モデルの運動の時間的変化を演算する機能と、を実現させ、前記機能は、前記摩擦係数を、浮上量に従って時刻毎に演算するとともに、前記浮上量は、前記モデルの各質点と前記搬送ローラの表面との間の距離として、時刻毎に演算されることを特徴とするプログラム。しかしながら、特許文献1および非特許文献1のいずれの技術によっても、ウェブ搬送装置により搬送される柔軟媒体の挙動を、この摩擦係数の変化に基づいてシミュレートすることができなかった。その結果、柔軟媒体の搬送状態がシミュレーションとウェブ装置とで相違するという問題が、生じていた。そこで、本発明では、柔軟媒体の挙動を良好にシミュレートすることができるシミュレーション装置およびプログラムを提供することを目的とする。上記の課題を解決するため、請求項1の発明は、搬送パラメータに従い搬送経路に沿って搬送される柔軟媒体の挙動を、モデルによりシミュレートするシミュレーション装置であって、前記搬送経路は、搬送ローラを含む複数の搬送要素により形成されており、前記モデルは、複数の質点と、各々が、前記複数の質点のうち隣接する2つの質点の間を連結する複数のバネとを有し、前記複数の搬送要素のそれぞれを、対応する位置情報にしたがって配置することにより、前記搬送経路を作成する経路作成部と、モデル情報に基づいて前記モデルを作成するモデル作成部と、前記搬送パラメータのうち、前記モデルの各質点と前記搬送ローラの表面との間の摩擦係数を演算するパラメータ演算部と、前記モデルが前記搬送パラメータに従って、前記複数の搬送要素により搬送される場合において、前記モデルの運動の時間的変化を演算する運動演算部とを備え、前記パラメータ演算部は、前記摩擦係数を、浮上量に従って時刻毎に演算するとともに、前記浮上量は、前記モデルの各質点と前記搬送ローラの表面との間の距離として、時刻毎に演算されることを特徴とする。また、請求項2の発明は、請求項1に記載のシミュレーション装置において、前記運動演算部により演算された前記モデルの運動の時間的変化を、前記柔軟媒体の挙動として表示させる表示部、をさらに備えることを特徴とする。また、請求項3の発明は、搬送パラメータに従い搬送経路に沿って搬送される柔軟媒体の挙動を、モデルによりシミュレートするコンピュータ読み取り可能なプログラムであって、前記搬送経路は、搬送ローラを含む複数の搬送要素により形成されており、前記モデルは、複数の質点と、各々が、前記複数の質点のうち隣接する2つの質点の間を連結する複数のバネとを有し、前記コンピュータによる前記プログラムの実行は、前記コンピュータに前記複数の搬送要素のそれぞれを、対応する位置情報にしたがって配置することにより、前記搬送経路を作成する機能と、モデル情報に基づいて前記モデルを作成するモデル作成部と、前記搬送パラメータのうち、前記モデルの各質点と前記搬送ローラの表面との間の摩擦係数を演算するパラメータ演算部と、前記モデルが、前記搬送パラメータに従って前記複数の搬送要素により搬送される場合において、前記モデルの運動の時間的変化を演算する機能とを実現させ、前記機能は、前記摩擦係数を、浮上量に従って時刻毎に演算するとともに、前記浮上量は、前記モデルの各質点と前記搬送ローラの表面との間の距離として、時刻毎に演算されることを特徴とする。請求項1から請求項3に記載の発明によれば、各時刻における摩擦係数は、モデルの各質点と搬送ローラの表面との間の浮上量に従って演算される。これにより、搬送ローラ付近における柔軟媒体の挙動を、浮上量および摩擦係数の変化に応じて把握することができる。そのため、柔軟媒体の挙動を、より正確にシミュレートすることが可能となる。特に、請求項2に記載の発明によれば、柔軟媒体の挙動を表示部に表示することができる。そのため、使用者は、柔軟媒体の挙動を視覚的に観察することができ、柔軟媒体の挙動を容易に把握することができる。以下、図面を参照しつつ本発明の実施の形態について詳細に説明する。<1.シミュレーション装置の構成>図1は、本発明の実施の形態におけるシミュレーション装置1の構成の一例を示すブロック図である。シミュレーション装置1は、ウェブ搬送装置で搬送される柔軟媒体の挙動を、モデルによりシミュレートする。図1に示すように、シミュレーション装置1は、主として、CPU10と、RAM20と、ROM30と、大容量記憶部40と、表示処理部50と、入出力部60と、を有している。また、CPU10と、RAM20と、ROM30と、大容量記憶部40と、表示処理部50と、入出力部60と、は、バス8を介して電気的に接続されている。RAMRandomAccessMemory20は、揮発性の記憶部であり、例えば、CPU10の演算で使用されるデータが格納されている。また、図1に示すように、RAM20は、位置情報格納テーブル21およびモデル情報格納テーブル22を記憶する。位置情報格納テーブル21は、柔軟媒体の搬送に用いられる各搬送要素例えば、搬送ローラ、直線ガイド、または円弧ガイド等と、対応する位置情報と、を関連付けて格納する。ここで、搬送要素は、例えば、上述の搬送ローラ、直線ガイド、または円弧ガイド等を含み、柔軟媒体が搬送される搬送経路を形成する。モデル情報格納テーブル22は、柔軟媒体を模倣するためのモデルに関する情報を格納している。ここで、柔軟媒体のモデルは、主として、複数の質点と、各々が複数の質点のうち隣接する2つの質点の間を連結する複数のバネと、を有している。ROMReadOnlyMemory30は、いわゆる不揮発性の記憶部であり、例えば、プログラム31が格納されている。なお、ROM30としては、読み書き自在の不揮発性メモリであるフラッシュメモリが使用されてもよい。大容量記憶部40は、シリコンディスクドライブやハードディスクドライブ等のようにRAM20と比較して記憶容量の大きな素子により構成された記憶部であり、必要に応じてRAM20との間でデータの授受を行う。CPUCentralProcessingUnit10は、ROM30のプログラム31に従った動作制御やデータ処理を実行する。また、図1中のCPU10内に記載されているブロックそれぞれ符号11~14が付与されているに対応する演算機能は、CPU10により実現される。なお、各ブロック11~14に対応する演算機能については後述する。表示処理部50は、いわゆるビデオコントローラにより構成されており、信号線50aを介して表示部51と電気的に接続されている。表示処理部50は、描画処理を実行することにより表示部51の画面に文字や図形等を表示させる。入出力部60は、いわゆるI/OポートInputOutputPortにより構成されており、信号線60aを介してシミュレーション装置1の入出力機器例えば、キーボードやタッチパッドを有する操作部61と電気的に接続されている。したがって、入出力部60は、CPU10と入出力機器との間でデータの授受を実行する。<2.CPUにより実現される演算機能>図2は、搬送経路70bの一例を説明するための正面図である。図3は、運動演算部14により演算された柔軟媒体80の搬送速度と、時刻tとの関係を示すグラフである。ここでは、図1から図3を参照しつつ、経路作成部11、モデル作成部12、パラメータ演算部13、および運動演算部14に対応する機能を説明する。なお、図2および以降の各図には、方向関係を明確にすべく必要に応じて適宜、Z軸方向を鉛直方向とし、XY平面を水平面とするXYZ直交座標系が、付されている。経路作成部11は、RAM20に格納された位置情報格納テーブル21に基づいて、配置空間70aに搬送経路70bを作成する。例えば、図2に示すように、複数の搬送要素従動ローラ71~73、テンションローラ74、および搬送ローラ75等のそれぞれは、対応する位置情報にしたがって配置空間70a内に配置される。これにより、従動ローラ71~73、テンションローラ74、および搬送ローラ75に沿った搬送経路70bが形成される。モデル作成部12は、RAM20に格納されたモデル情報格納テーブル22に基づいて、柔軟媒体80のモデルを作成する。例えば、図2の柔軟媒体80は、上述のように、複数の質点と、複数のバネと、によって模倣される。ここで、経路作成部11で用いられる位置情報格納テーブル21と、モデル作成部12で用いられるモデル情報格納テーブル22は、シミュレーション装置1において、操作部61を介したシミュレーション装置1の使用者の指示に従い作成されても良い。また、位置情報格納テーブル21およびモデル情報格納テーブル22は、シミュレーション装置1とは異なった情報処理装置で、その情報処理装置の使用者により作成されても良い。さらに、位置情報格納テーブル21およびモデル情報格納テーブル22は、シミュレーションにより作成されてもよい。パラメータ演算部13は、搬送装置で用いられる搬送パラメータを演算する。例えば、パラメータ演算部13は、搬送パラメータのうち、柔軟媒体80のモデルの各質点と、搬送ローラ75の表面と、の間の有効摩擦係数μeffを演算する。なお、パラメータ演算部13によって実行される演算の詳細は、後述する。運動演算部14は、柔軟媒体80のモデルが搬送パラメータに従って複数の搬送要素により搬送される場合において、この柔軟媒体80のモデルの運動の時間的変化を演算する。例えば、図3には、有効摩擦係数μeffの時間的変化を考慮した場合における柔軟媒体80のモデルの搬送速度の時間的変化と、有効摩擦係数μeffの時間的変化を考慮しない場合における柔軟媒体80のモデルの搬送速度の時間的変化と、が、示されている。そして、表示部51は、運動演算部14により演算された柔軟媒体モデルの運動の時間的変化に基づいて、柔軟媒体の挙動を表示する。例えば、図3に示すように、運動演算部14は、各時刻tにおける柔軟媒体の速度を、表示部51に表示させることができる。すなわち、表示部51は、運動演算部14により演算されたモデルの運動の時間的変を、柔軟媒体の挙動として表示することができる。そのため、シミュレーション装置1の使用者は、柔軟媒体の挙動を視覚的に観察することができ、容易に柔軟媒体の挙動を把握することができる。<3.パラメータ演算部による演算>図4は、搬送ローラ75付近の構成の一例を示す正面図である。ここでは、図4を参照しつつ、有効摩擦係数μefftの演算手法について説明する。ここで、図4に示すように、搬送ローラ75が矢印R1方向に回転することによって、柔軟媒体80が矢印AR1方向以下、搬送方向とも称するに搬送される場合、搬送ローラ75の周速度Urによっては、搬送ローラ75と柔軟媒体80の間に空気が巻き込まれる。その結果、柔軟媒体80が搬送ローラ75に対して浮上する。この場合において、搬送ローラ75に対する柔軟媒体80の浮上量hが大きくなるに従って、柔軟媒体80と搬送ローラ75の表面との間の摩擦係数は小さくなる。そして、浮上量hの値が一定値以上となると、この摩擦係数が0になり、柔軟媒体80に対して搬送ローラ75が空転する。また、この浮上量hは、柔軟媒体80のモデルの各質点と搬送ローラ75の表面との間の距離であり、式1のように表すことができる。ただし、式1中において、aRは、搬送ローラ75の半径を、bΘは、巻角を、cσは、柔軟媒体80および搬送ローラ75の間の合成表面粗さを、dηは、空気粘度を、eUrは、搬送ローラ75の周速度を、fUwは、搬送経路70bに沿った柔軟媒体80の搬送速度を、gT1は、柔軟媒体80に負荷される張力を、hkは、柔軟媒体80の透気率を、itwは、柔軟媒体80の厚さを、ixは、搬送経路70b上における位置を、それぞれ示すものである。なお、式1中の位置xの値は、原点P0より搬送方向下流側巻付開始位置P1側が負数、原点P0より搬送方向上流側巻付終了位置P2側が正数となるように、設定されている。また、搬送ローラ75の周速度Urおよび柔軟媒体80の搬送速度Uwは、時刻tの関数と見ることができる。したがって、式1は、式2のように変形することができ、浮上量hは、各位置xにおいて時刻t毎に演算できる。ただし、式2のa、b、c、pt、およびqtは、それぞれ式3、式4、式5、式6、および式7のように表される。したがって、巻付開始位置P1から巻付終了位置P2までの浮上量hの平均値havetは、式8のように表される。ただし、式8中において、jμsは、初期摩擦係数を、示すものである。また、搬送ローラ75に対する柔軟媒体80の摩擦係数μ1は、式9のように表すことができる。よって、式9に式8を代入すると、有効摩擦係数μefftは、式10のように表すことができる。式10から分かるように、パラメータ演算部13は、予め定められた初期摩擦係数μsと、浮上量hの平均値havetと、から有効摩擦係数μefftを演算することができる。すなわち、パラメータ演算部13は、有効摩擦係数μefftを浮上量hの平均値havetに従って時刻t毎に演算することができる。これにより、空気浮上により有効摩擦係数μefftが変化することを考慮した柔軟媒体80の搬送速度Uwを、容易に演算することができる。そのため、柔軟媒体80の挙動を、計算コストを増大させることなく、正確にシミュレートすることができる。なお、空気浮上とは、搬送ローラ75および柔軟媒体80の間に空気が巻き込まれ、搬送ローラ75に対して柔軟媒体80が浮上することを言うものとする。<4.シミュレーション結果>図5は、パラメータ演算部13により演算された浮上量hの平均値havetおよび有効摩擦係数μefftと、時刻tとの関係を示すグラフである。図6は、運動演算部14により演算された柔軟媒体80の搬送速度Uwと、時刻tとの関係を示すグラフである。図5において、横軸は時刻tを示す。また、図5において、左側の縦軸は、有効摩擦係数μefftを示し、点線が対応する。一方、右側の縦軸は、浮上量hの平均値havetを示し、実線が対応する。図3および図6において、横軸は時刻tを示し、縦軸は柔軟媒体80の搬送速度Uwを示す。また、図3および図6中の点線は、空気浮上により有効摩擦係数μefftが変化することを考慮した場合おける柔軟媒体80の搬送速度Uwを示す。一方、図3および図6中の実線は、空気浮上を考慮しない場合における柔軟媒体80の搬送速度Uwを示す。なお、空気浮上とは、搬送ローラ75および柔軟媒体80の間に空気が巻き込まれ、搬送ローラ75に対して柔軟媒体80が浮上することを言うものとする。また、本シミュレーションにおいて、合成表面粗さσの値は、0.583μmに設定されている。図5および図6に示すように、時刻t1において、搬送ローラ75に対して柔軟媒体80が浮上し始める。ただし、時刻t1において、浮上量hの平均値havet1は合成表面粗さσより小さいため、有効摩擦係数μefft1は、初期摩擦係数μsとなる。この場合において、搬送ローラ75の周速度Urと、柔軟媒体80の周速度Urは、同一値である。次に、時刻t2において、浮上量hの平均値havet2は合成表面粗さσ以上となり、有効摩擦係数μeffが減少し始める。また、図6に示すように、時刻t2から時刻t3までの期間において、搬送ローラ75の周速度Urと、柔軟媒体80の周速度Urと、は、同一値となる。続いて、時刻t3において、搬送ローラ75に対して柔軟媒体80が滑り始める。これにより、図6に示すように、柔軟媒体80の搬送速度Uwが、搬送ローラ75の周速度Ur以下となる。続いて、時刻t4において、浮上量hの平均値havet4は、合成表面粗さσの3倍となり、有効摩擦係数μeffが0となる。すなわち、時刻t5において、柔軟媒体80の搬送速度Uwが一定値となる。すなわち、少なくとも時刻t5から時刻t7において、搬送ローラ75の周速度Urが増大しても、柔軟媒体80の搬送速度Uwは変化しない。そして、時刻t6において搬送ローラ75の周速度Urが一定になると図3参照、時刻t7以降において浮上量hの平均値havet7は一定値となる。<5.本実施の形態のシミュレーション装置の利点>以上のように、本実施の形態のシミュレーション装置1は、各時刻における有効摩擦係数μeffは、柔軟媒体80のモデルの各質点と搬送ローラ75の表面との間の浮上量hの平均値havetに従って演算される。これにより、搬送ローラ75付近における柔軟媒体80の挙動を、浮上量hの平均値havetおよび有効摩擦係数μeffの変化に応じて把握することができる。例えば、有効摩擦係数μeffを考慮に入れた柔軟媒体80の搬送速度を求めることが可能となる。そのため、柔軟媒体80の挙動が、より正確にシミュレートできる。<6.変形例>以上、本発明の実施の形態について説明してきたが、本発明は上記実施の形態に限定されるものではなく様々な変形が可能である。1本実施の形態のシミュレーションでは、一方向に延びる連続用紙図2参照が柔軟媒体80として用いられているが、これに限定されるものでない。例えば、縦横サイズが規格化された固定サイズの用紙が、柔軟媒体80として用いられても良い。2また、本実施の形態において、経路作成部11、モデル作成部12、パラメータ演算部13、および運動演算部14により実現される機能は、いずれもCPU10によりソフトウェア的に実現されるものとして説明したが、これに限定されるものでない。例えば、これらの機能は、電子回路等のハードウェアによって実現されてもよい。以下、本発明を実施例によりさらに詳細に説明するが、本発明はこれらの例に何ら限定されるものではない。<実施例1>図2は、搬送経路70bの一例を示すとともに、本実施の形態の有効性の確認を行うための実機試験を実行した実験装置70の構成の一例を示す正面図である。図7は、解析例および実機試験の試験条件を説明するための図である。図8は、解析例および実機試験の試験結果を示すグラフである。ここで、図8において、解析例と実機試験の駆動ローラの速度Urは同じ破線で表す。実線および点線は、それぞれ解析例および実機試験のウェブ速度Uwを示す。また、図7に示すように、解析例および実機試験では、PTMEで形成された連続用紙が柔軟媒体80として使用されている。また、解析例は、空気浮上により有効摩擦係数μefftが変化することを考慮して、柔軟媒体80の搬送速度Uwをシミュレートし、その演算結果をグラフにプロットしたものである。さらに、実機試験は、実験装置70により計測された柔軟媒体80の搬送速度Uwをプロットしたものである。図8に示すように、実機試験の試験結果は、解析例の試験結果に符合した。このように、空気浮上により有効摩擦係数μefftが変化することを考慮することによって、柔軟媒体80の挙動が、より正確にシミュレートできた。