"コンピュータシミュレーションにより数値解析用の繊維材料のフルアトム分子動力学モデルを作成する繊維材料モデル作成方法であって、解析対象となる繊維材料のモノマーモデルを作成する工程と、前記モノマーモデルを重合させて単位モデルとなる重合体モデルを作成する工程と、前記重合体モデルの複数本を予め定めた空間にランダムに配置したランダム配置モデルを作成する工程と、前記ランダム配置モデルを分子動力学計算により構造緩和させて未配向のモデルを作成する工程と、前記未配向のモデルに対し、1軸方向にせん断エネルギーを与え、分子鎖を1軸方向に伸張する第1の配向処理を行う工程と、前記第1の配向処理により、伸びきり状態の1本の分子鎖を選択し、選択した前記分子鎖を複数本複写した分子鎖を作成し、前記各複写した分子鎖を繊維材料の結晶構造に対応した座標に配置したモデルを作成する第2の配向処理を行う工程と、前記第2の配向処理直後の前記モデルを分子動力学計算により構造緩和させて繊維材料モデルを作成する工程と、を備えることを特徴とする繊維材料モデル作成方法。コンピュータシミュレーションにより数値解析用の接着処理が施された繊維材料とゴムまたは樹脂、あるいはこれらの混合物からなる接着液の界面モデルを作成する繊維材料と接着液の界面モデル作成方法であって、分子動力学計算により解析対象となる接着液を構成するゴムまたは樹脂、あるいはこれらのモノマーモデルを作成する工程と、前記ゴムまたは樹脂、あるいはこれらのモノマーモデルをそれぞれ所定の重合度で重合させて単位モデルとなる前記ゴムまたは樹脂、あるいはこれらの重合体モデルを作成する工程と、前記ゴムまたは樹脂、あるいはこれらの重合体モデルの所定の混合比率で予め定めた空間に配置する工程と、分子動力学計算により構造緩和させて接着液モデルを作成する工程と、前記接着液モデルを、1軸配向した請求項1に記載の前記繊維材料モデルに隣接して配置する工程と、分子動力学計算により構造緩和させて繊維材料と接着液の界面モデルを作成する工程と、を備えることを特徴とする繊維材料と接着液の界面モデル作成方法。コンピュータシミュレーションにより繊維材料と接着液の界面の剥離特性を予測する繊維材料と接着液の界面の剥離特性予測方法であって、請求項3に記載の前記繊維材料と接着液の界面モデルを、界面の厚さ方向に1軸伸張させて応力-変位曲線を計算する工程と、前記応力-変位曲線に基づいて、剥離応力の最大値を算出する工程と、前記応力-変位曲線により与えられる変位に対して応力を積分することにより、剥離エネルギーを算出する工程と、を備えることを特徴とする繊維材料と接着液の界面の剥離特性予測方法。コンピュータ上に実装され、コンピュータシミュレーションにより数値解析用の繊維材料のフルアトム分子動力学モデルを作成する繊維材料モデル作成装置であって、解析対象となる繊維材料のモノマーモデルを作成するモノマーモデル作成手段と、前記モノマーモデルを重合させて単位モデルとなる重合体モデルを作成する重合体モデル作成手段と、前記重合体モデルの複数本を予め定めた空間にランダムに配置したランダム配置モデルを作成するランダム配置モデル作成手段と、前記ランダム配置モデルを分子動力学計算により構造緩和させて未配向のモデルを作成する未配向モデル作成手段と、前記未配向のモデルに対し、1軸方向にせん断エネルギーを与え、分子鎖を1軸方向に伸張する第1の配向処理を行う第1の配向処理手段と、前記第1の配向処理により、伸びきり状態の1本の分子鎖を選択し、選択した前記分子鎖を複数本複写した分子鎖を作成し、前記各複写した分子鎖を繊維材料の結晶構造に対応した座標に配置したモデルを作成する第2の配向処理を行う第2の配向処理手段と、前記第2の配向処理直後の前記モデルを分子動力学計算により構造緩和させて繊維材料モデルを作成する繊維材料モデル作成手段と、を備えることを特徴とする繊維材料モデル作成装置。コンピュータ上に実装され、コンピュータシミュレーションにより数値解析用の接着処理が施された繊維材料とゴムまたは樹脂、あるいはこれらの混合物からなる接着液の界面モデルを作成する繊維材料と接着液の界面モデル作成装置であって、分子動力学計算により解析対象となる接着液を構成するゴムまたは樹脂、あるいはこれらのモノマーモデルを作成するモノマーモデル作成手段と、前記ゴムまたは樹脂、あるいはこれらのモノマーモデルをそれぞれ所定の重合度で重合させて単位モデルとなる前記ゴムまたは樹脂、あるいはこれらの重合体モデルを作成する重合体モデル作成手段と、前記ゴムまたは樹脂、あるいはこれらの重合体モデルの所定の混合比率で予め定めた空間に配置するランダム配置モデル作成手段と、分子動力学計算により構造緩和させて接着液モデルを作成する接着液モデル作成手段と、前記接着液モデルを、1軸配向した請求項4に記載の前記繊維材料モデルに隣接して配置する接着液モデル及び繊維材料モデル配置手段と、分子動力学計算により構造緩和させて繊維材料と接着液の界面モデルを作成する界面モデル作成手段と、を備えることを特徴とする繊維材料と接着液の界面モデル作成装置。コンピュータ上に実装され、コンピュータシミュレーションにより繊維材料と接着液の界面の剥離特性を予測する繊維材料と接着液の界面の剥離特性予測装置であって、請求項5に記載の前記繊維材料と接着液の界面モデルを、界面の厚さ方向に1軸伸張させて応力-変位曲線を計算する1軸伸張手段と、前記応力-変位曲線に基づいて、剥離応力の最大値を算出する最大剥離応力算出手段と、前記応力-変位曲線により与えられる変位に対して応力を積分することにより、剥離エネルギーを算出する剥離エネルギー算出手段と、を備えることを特徴とする繊維材料と接着液の界面の剥離特性予測装置。そこで、コンピュータシミュレーションにより、繊維材料/接着液の接着メカニズムを解明して、繊維材料と接着液の界面の剥離特性を予測するには、各種材料繊維材料と接着液などが共存する状態においても、繊維材料が配向状態を維持できるような安定した配向構造及び界面構造を作成する必要がある。即ち、前提として、短時間で1軸方向に配向した繊維材料の分子モデルを作成する必要がある。本発明が解決しようとする課題は、コンピュータシミュレーションにより繊維材料/接着液の接着メカニズムの解析と繊維材料と接着液の界面の剥離特性の予測を行うことを目標とし、その前提となる繊維材料のフルアトム分子動力学モデルを作成する繊維材料モデル作成方法と繊維材料モデル作成装置、及び、作成された繊維材料モデルとゴムまたは樹脂、あるいはこれらの混合物が混在した接着液のフルアトム分子動力学モデルに基づいて繊維材料と接着液の界面モデルを作成する繊維材料と接着液の界面モデル作成方法と繊維材料と接着液の界面モデル作成装置、並びに、作成された繊維材料と接着液の界面モデルに基づいて、繊維材料と接着液の界面の剥離特性予測方法と繊維材料と接着液の界面の剥離特性予測装置を提供するものである。本発明に係る繊維材料モデル作成方法は、コンピュータシミュレーションにより数値解析用の繊維材料のフルアトム分子動力学モデルを作成する繊維材料モデル作成方法であって、解析対象となる繊維材料のモノマーモデルを作成する工程と、前記モノマーモデルを重合させて単位モデルとなる重合体モデルを作成する工程と、前記重合体モデルの複数本を予め定めた空間にランダムに配置したランダム配置モデルを作成する工程と、前記ランダム配置モデルを分子動力学計算により構造緩和させて未配向のモデルを作成する工程と、前記未配向のモデルに対し、1軸方向にせん断エネルギーを与え、分子鎖を1軸方向に伸張する第1の配向処理を行う工程と、前記第1の配向処理により、伸びきり状態の1本の分子鎖を選択し、選択した前記分子鎖を複数本複写した分子鎖を作成し、前記各複写した分子鎖を繊維材料の結晶構造に対応した座標に配置したモデルを作成する第2の配向処理を行う工程と、前記第2の配向処理直後の前記モデルを分子動力学計算により構造緩和させて繊維材料モデルを作成する工程と、を備えることを特徴とする。本発明に係る繊維材料モデル作成装置は、コンピュータ上に実装され、コンピュータシミュレーションにより数値解析用の繊維材料のフルアトム分子動力学モデルを作成する繊維材料モデル作成装置であって、解析対象となる繊維材料のモノマーモデルを作成するモノマーモデル作成手段と、前記モノマーモデルを重合させて単位モデルとなる重合体モデルを作成する重合体モデル作成手段と、前記重合体モデルの複数本を予め定めた空間にランダムに配置したランダム配置モデルを作成するランダム配置モデル作成手段と、前記ランダム配置モデルを分子動力学計算により構造緩和させて未配向のモデルを作成する未配向モデル作成手段と、前記未配向のモデルに対し、1軸方向にせん断エネルギーを与え、分子鎖を1軸方向に伸張する第1の配向処理を行う第1の配向処理手段と、前記第1の配向処理により、伸びきり状態の1本の分子鎖を選択し、選択した前記分子鎖を複数本複写した分子鎖を作成し、前記各複写した分子鎖を繊維材料の結晶構造に対応した座標に配置したモデルを作成する第2の配向処理を行う第2の配向処理手段と、前記第2の配向処理直後の前記モデルを分子動力学計算により構造緩和させて繊維材料モデルを作成する繊維材料モデル作成手段と、を備えることを特徴とする。これによると、分子動力学モデルを用いたコンピュータシミュレーションにより、比較的短時間で1軸方向に配向した繊維材料のフルアトム分子動力学モデルを作成することができる。そして、作成した繊維材料モデルを用いて、繊維材料/ゴムまたは樹脂、あるいはこれらの混合物からなる接着液の密着性を予測することが可能になる。また、本発明に係る繊維材料と接着液の界面モデル作成方法は、コンピュータシミュレーションにより数値解析用の接着処理が施された繊維材料とゴムまたは樹脂、あるいはこれらの混合物からなる接着液の界面モデルを作成する繊維材料と接着液の界面モデル作成方法であって、分子動力学計算により解析対象となる接着液を構成するゴムまたは樹脂、あるいはこれらのモノマーモデルを作成する工程と、前記ゴムまたは樹脂、あるいはこれらのモノマーモデルをそれぞれ所定の重合度で重合させて単位モデルとなる前記ゴムまたは樹脂、あるいはこれらの重合体モデルを作成する工程と、前記ゴムまたは樹脂、あるいはこれらの重合体モデルの所定の混合比率で予め定めた空間に配置する工程と、分子動力学計算により構造緩和させて接着液モデルを作成する工程と、前記接着液モデルを、1軸配向した請求項1に記載の前記繊維材料モデルに隣接して配置する工程と、分子動力学計算により構造緩和させて繊維材料と接着液の界面モデルを作成する工程と、を備えることを特徴とする。本発明に係る繊維材料と接着液の界面モデル作成装置は、コンピュータ上に実装され、コンピュータシミュレーションにより数値解析用の接着処理が施された繊維材料とゴムまたは樹脂、あるいはこれらの混合物からなる接着液の界面モデルを作成する繊維材料と接着液の界面モデル作成装置であって、分子動力学計算により解析対象となる接着液を構成するゴムまたは樹脂、あるいはこれらのモノマーモデルを作成するモノマーモデル作成手段と、前記ゴムまたは樹脂、あるいはこれらのモノマーモデルをそれぞれ所定の重合度で重合させて単位モデルとなる前記ゴムまたは樹脂、あるいはこれらの重合体モデルを作成する重合体モデル作成手段と、前記ゴムまたは樹脂、あるいはこれらの重合体モデルの所定の混合比率で予め定めた空間に配置するランダム配置モデル作成手段と、分子動力学計算により構造緩和させて接着液モデルを作成する接着液モデル作成手段と、前記接着液モデルを、1軸配向した請求項4に記載の前記繊維材料モデルに隣接して配置する接着液モデル及び繊維材料モデル配置手段と、分子動力学計算により構造緩和させて繊維材料と接着液の界面モデルを作成する界面モデル作成手段と、を備えることを特徴とする。これによると、作成した1軸方向に配向した繊維材料のフルアトム分子動力学モデルとゴムまたは樹脂、あるいはこれらの混合物が混在した接着液のフルアトム分子動力学モデルに基づいて繊維材料と接着液の界面モデルを作成することができる。そして、作成された繊維材料と接着液の界面モデルを使って剥離解析を行うことにより、界面における材料の配合と剥離特性の関係とを予測することができ、ひいては、実験では把握しがたい界面におけるナノミクロ領域の挙動と剥離特性の関係を予測することができる。更に、繊維材料と接着液の界面の剥離特性予測方法は、コンピュータシミュレーションにより繊維材料と接着液の界面の剥離特性を予測する繊維材料と接着液の界面の剥離特性予測方法であって、請求項3に記載の前記繊維材料と接着液の界面モデルを、界面の厚さ方向に1軸伸張させて応力-変位曲線を計算する工程と、前記応力-変位曲線に基づいて、剥離応力の最大値を算出する工程と、前記応力-変位曲線により与えられる変位に対して応力を積分することにより、剥離エネルギーを算出する工程と、を備えることを特徴とする。本発明に係る繊維材料と接着液の界面の剥離特性予測装置は、コンピュータ上に実装され、コンピュータシミュレーションにより繊維材料と接着液の界面の剥離特性を予測する繊維材料と接着液の界面の剥離特性予測装置であって、請求項5に記載の前記繊維材料と接着液の界面モデルを、界面の厚さ方向に1軸伸張させて応力-変位曲線を計算する1軸伸張手段と、前記応力-変位曲線に基づいて、剥離応力の最大値を算出する最大剥離応力算出手段と、前記応力-変位曲線により与えられる変位に対して応力を積分することにより、剥離エネルギーを算出する剥離エネルギー算出手段と、を備えることを特徴とする。これによると、作成した繊維材料と接着液の界面モデルを界面の厚さ方向に1軸伸張させて、剥離挙動を解析することにより、界面における材料の配合と剥離特性の関係とを予測することができ、実験では把握しがたい繊維材料と接着液の界面におけるナノミクロ領域の挙動と剥離特性の関係を予測することができる。本発明の繊維材料モデル作成方法と繊維材料モデル作成装置、及び、繊維材料と接着液の界面モデル作成方法と繊維材料と接着液の界面モデル作成装置、並びに、繊維材料と接着液の界面の剥離特性予測方法と繊維材料と接着液の界面の剥離特性予測装置は、分子動力学モデルを用いたコンピュータシミュレーションにより、比較的短時間で1軸方向に配向した繊維材料のフルアトム分子動力学モデルを作成して、繊維材料/ゴムまたは樹脂、あるいはこれらの混合物からなる接着液の密着性を予測することが可能になると共に、作成した1軸方向に配向した繊維材料のフルアトム分子動力学モデルとゴムまたは樹脂、あるいはこれらの混合物が混在した接着液のフルアトム分子動力学モデルに基づいて繊維材料と接着液の界面モデルを作成して、繊維材料と接着液の界面モデルを使って剥離解析を行うことにより、界面における材料の配合と剥離特性の関係とを予測することができ、ひいては、実験では把握しがたい繊維材料と接着液の界面におけるナノミクロ領域の挙動と剥離特性の関係を予測することができる。以下、図面を参照しつつ、本発明に係る繊維材料モデル作成方法と繊維材料モデル作成装置、及び、繊維材料と接着液の界面モデル作成方法と繊維材料と接着液の界面モデル作成装置、並びに、繊維材料と接着液の界面の剥離特性予測方法と繊維材料と接着液の界面の剥離特性予測装置を実施するための形態について、具体的な一例に即して説明する。尚、以下に説明するものは、例示したものにすぎず、本発明に係る繊維材料モデル作成方法と繊維材料モデル作成装置、及び、繊維材料と接着液の界面モデル作成方法と繊維材料と接着液の界面モデル作成装置、並びに、繊維材料と接着液の界面の剥離特性予測方法と繊維材料と接着液の界面の剥離特性予測装置の適用限界を示すものではない。すなわち、本発明に係る繊維材料モデル作成方法と繊維材料モデル作成装置、及び、繊維材料と接着液の界面モデル作成方法と繊維材料と接着液の界面モデル作成装置、並びに、繊維材料と接着液の界面の剥離特性予測方法と繊維材料と接着液の界面の剥離特性予測装置は、下記の実施形態に限定されるものではなく、特許請求の範囲に記載した限りにおいてさまざまな変更が可能なものである。[繊維材料モデル作成方法]まず、本実施形態に係る繊維材料モデル作成方法の一例について、図1及び図2に基づいて説明する。図1は、本実施形態に係る繊維材料モデル作成方法の手順を示す概略フローチャート図であり、図2は、本実施形態に係る繊維材料モデル作成方法の手順を示す詳細フローチャート図である。ここで、本実施形態において作成する繊維材料モデルは、繊維材料のフルアトム分子動力学モデルを意味している。図1に示すように、最初に、解析対象となる繊維材料について、原子間の相互作用を厳密に再現したモノマーモデルを、分子動力学計算により作成するステップS1。より詳細には、図2に示すように、原子の座標,結合の設定を行いステップS101、Dreiding,GAFF,AMBER等の力場を用いて、結合,結合角,二面角エネルギー,非結合エネルギーVanderWaalsを設定データベースを選択して各原子に割り当てるするステップS102。次に、モンテカルロ法又はミニマイズ法等を用いて構造最適化を行い、分子鎖のエネルギーが最小となる原子配置を計算するステップS103。そして、分子鎖の中の原子が保有する電荷を分子軌道法で計算し、分子鎖を形成する各電子が保有する電荷を、密度汎関数法により求めたモノマーの電荷を用いて計算するS104。以上により、繊維材料のモノマーモデルのモノマーファイルを作成するS105。次に、図1に示すように、ステップS1で作成したモノマーモデルを所定の重合度で重合させて単位モデルである重合体モデルを作成するステップS2。より詳細には、図2に示すように、モノマーモデルの重合度重合させる数と、重合体モデルの立体規則性の種類を設定するステップS201。ここで、立体規則性の種類として、立体規則性の影響を調査するために、立体規則性がアイソタクチックすべての不斉炭素が同じ絶対配置を持つ構造、シンジオタクチック絶対配置が交互に並ぶ構造、アタクチック全くランダムな構造の3種類を作成することが好ましい。そして、Dreiding,GAFF,AMBER等の力場を用いて、結合,結合角,二面角エネルギー,非結合エネルギーVanderWaalsを設定データベースを選択して各原子に割り当てるするステップS202。次に、モンテカルロ法又はミニマイズ法等を用いて構造最適化を行い、分子鎖のエネルギーが最小となる原子配置を計算するステップS203。そして、分子鎖の中の原子が保有する電荷を分子軌道法で計算し、密度汎関数法またはモノマーモデルの値を利用して、重合後の分子鎖の電荷の総和がゼロとなるように調整しながら任意の原子が保有する電荷を計算するステップS204。そして、図1に示すように、ステップS2で作成した重合体モデルの複数本を、予め定めたセル空間内のランダムな位置に配置したランダム配置モデルを作成するステップS3。より詳細には、図2に示すように、重合体モデルの数と、重合体モデルの立体規則性の種類を選択しステップS301、分子動力学計算用のファイルここで、ファイルとは原子の種類、座標、原子間の結合、相互作用に関する情報、温度、圧力等の計算条件等を含むものである。を作成するステップS302。次に、図1に示すように、ステップS3で作成したランダム配置モデルを、分子動力学計算により構造緩和させて、温度,圧力等が指定した動的な計算の中でエネルギーが最小となる原子配置を計算し、未配向のモデルを作成するステップS4。より詳細には、図2に示すように、分子鎖を熱平衡に近い状態にするために、高圧条件で、NPT粒子数,圧力,温度一定条件計算を実施する緩和計算1:ステップS401。次に、圧力を常圧に設定して、再度NPT計算を実施する緩和計算2:ステップS402。最後に、圧力が常圧で、NVT粒子数,体積,温度一定条件計算を実施する緩和計算3:ステップS403。なお、上記NPT計算及びNVT計算において、例えば、温度制御はNose-Hooverの方法を、圧力制御はAndersenの方法を用いる。以上により、未配向のモデルが作成されるステップS5。そして、ステップS5で作成された未配向のモデルに対し、セルの境界部分端部に歪を与えるLees-Edwards法等を用いて、モデルの1軸方向にせん断エネルギーを与え、分子鎖を1方向z軸方向=配向軸に引き伸ばす第1の配向処理を行うステップS6。そして、ステップS6で第1の配向処理を行った分子鎖の中から伸びきり状態の1本の分子鎖を選択し、この分子鎖を複写して同じ伸張状態の複数本の分子鎖を作成し、それぞれの分子鎖を繊維材料の結晶構造に対応した座標に配置したモデルを作成する第2の配向処理を行うステップS7。そして、ステップS7で行った第2の配向処理直後のモデルに対して、構造緩和計算を行うステップS8。より詳細には、図2に示すように、ステップS7で行った第2の配向処理後の分子鎖はz軸方向に引き伸ばされており、xyz軸の各方向に異方応力が作用しているため、第2の配向処理直後のモデルに対して、セル長さのxy方向xy断面積を固定して、高圧条件で、NPT計算を実施する緩和計算1:ステップS801。そして、z軸方向だけを固定して、再度、高圧条件で、NPT計算を実施する緩和計算2:ステップS802。次に、セルの固定を解除して、圧力を常圧に設定して、再度NPT計算を実施する緩和計算3:ステップS803。最後に、圧力が常圧で、NVT計算を実施する緩和計算4:ステップS804。尚、上記NPT計算及びNVT計算において、例えば、温度制御はNose-Hooverの方法を、圧力制御はParrinello-Rahmanの方法を用いる。以上により、z軸方向に配向した繊維材料モデルが作成されるステップS9。[繊維材料モデルの評価]以上により作成された繊維材料モデルの妥当性の評価として、まず、作成された繊維材料モデルの比重と配向度と現実の繊維材料の比重と配向度とを比較して、作成した繊維材料モデルの比重と配向度が、現実の繊維材料の比重と配向度の文献値に近い値が得られているのが妥当と評価することにより行う。また、作成された繊維材料モデルの妥当性を判断するために、力学物性の評価も必要である。そこで、作成された繊維材料モデルの1軸伸張シミュレーションから得られた応力-ひずみ曲線から弾性率を算出し、その結果を繊維材料の理論弾性率及びケブラーフィラメントの弾性率と比較することによりモデルの妥当性を評価する。即ち、作成された繊維材料モデルの弾性率を算出するために、作成された繊維材料モデルに対して、所定の伸張速度,ポアソン比,NVTアンサンブルの条件で伸張シミュレーションを実施する。次に、伸張シミュレーション結果から得られたモデルの応力-ひずみ曲線に基づいて、線形領域の傾きから弾性率を算出する。そして、算出された弾性率の値と、繊維材料の理論弾性率及びケブラーフィラメントの弾性率と比較する。そして、作成した繊維材料モデルは分子鎖が複数本なので、算出された弾性率の値は、1本鎖で計算する理論弾性率例えば、繊維材料がアラミド繊維の場合は、182GPa日本ポリイミド芳香族系高分子研究会編新訂最新ポリイミド-基礎と応用-1292010参照よりも低くなるのが妥当と評価される。また、算出された弾性率の値が、ケブラーフィラメントの弾性率の物性値に近い値が得られているのが妥当と評価される。このように、本実施形態の繊維材料モデル作成方法によれば、分子動力学モデルを用いたコンピュータシミュレーションにより、比較的短時間で1軸方向に配向した繊維材料のフルアトム分子動力学モデルを作成することができる。そして、作成した繊維材料モデルを用いて、繊維材料/ゴムまたは樹脂、あるいはこれらの混合物からなる接着液の密着性を予測することが可能になる。[繊維材料と接着液の界面モデル作成方法]次に、本実施形態に係る繊維材料と接着液の界面モデル作成方法の一例について、図3に基づいて説明する。図3は、本実施形態に係る繊維材料と接着液の界面モデル作成方法の手順を示すフローチャート図である。ここで、本実施形態では、接着液を用いて接着処理を施した繊維材料を解析対象とし、上述の実施形態により作成された繊維材料を用いて、繊維材料と接着液の界面モデルを作成する。図3に示すように、まず、接着液を構成するゴムまたは樹脂、あるいはこれらのモノマーモデルを、上述した繊維材料モデル作成方法において説明した図1及び図2におけるステップS1と同様の方法で、それぞれ作成するステップS11。尚、図3に示すステップS11~15において、接着液のフルアトムモデル以下、接着液モデルと略する。を作成するが、上述した繊維材料モデル作成方法において説明した図1及び図2におけるステップS1~5と同様の方法で作成するため、以下、詳細な説明は省略する。そして、ゴムまたは樹脂、あるいはこれらの重合体モデルを、上述した繊維材料モデル作成方法において説明した図1及び図2におけるステップS2と同様の方法で、それぞれ作成するステップS12。次に、上述した繊維材料モデル作成方法において説明した図1及び図2におけるステップS3と同様の方法で、それぞれ作成したゴムまたは樹脂、あるいはこれらの重合体モデルを、所定の混合比混合比として、例えば、重量比や分子数比となるように予め定めたセル空間にランダムに配置して、ランダム配置モデルを作成するステップS13。次に、上述した繊維材料モデル作成方法において説明した図1及び図2におけるステップS4と同様の方法で、ランダム配置モデルを、分子動力学計算により構造緩和させて、温度,圧力等が指定した動的な計算の中でエネルギーが最小となる原子配置を計算し、未配向の接着液モデルを作成するステップS14。ここで、分子鎖を熱平衡に近い状態にするために、高圧条件で、NPT粒子数,圧力,温度一定条件計算を実施する緩和計算1:ステップS401。次に、圧力を常圧に設定して、再度NPT計算を実施する緩和計算2:ステップS402。最後に、圧力が常圧で、NVT粒子数,体積,温度一定条件計算を実施する緩和計算3:ステップS403。以上により、上述した繊維材料モデル作成方法において説明した図1及び図2におけるステップS1~5と同様の方法で、接着液モデルが作成されるステップS15。尚、上述したゴムまたは樹脂、あるいはこれらの重合度、混合比を変化させることにより、種類の異なるモデルを作成することが可能である。そして、作成された接着液モデルを、z軸方向に配向した繊維材料モデル上述した本実施形態に係る繊維材料モデル作成方法において作成の下部に配置し、界面モデルを作成する。尚、配置において、繊維配向軸をz軸とし、厚さ方向をx軸とする。また、x軸の境界は反射境界条件を、その他の境界は周期境界条件を用いる。ここで、分子鎖の配置は、繊維材料モデルの配向軸であるz軸と垂直な軸であるy軸を基準として、これらのセルの長さを固定し、接着液モデルを再配置するステップS16。更に、作成された界面モデルに対して、分子鎖を熱平衡化させるために、以下のステップ緩和計算1~4で緩和計算を実施するステップS17。緩和計算1:接着液分子だけを熱平衡化させるために、繊維材料モデルの分子鎖のxyz座標を固定し、セルのyz軸断面積を固定して、高圧条件で、NPT計算を実施する。緩和計算2:繊維材料モデルの分子鎖の配向を維持し、分子鎖とx軸の境界面間の空隙を削除するために、分子全体をx軸方向に移動させ、隙間をつめる。高圧条件で、分子鎖のxyz座標の固定を解除し、セルのyz断面積を固定した状態で、NPT計算を実施する。緩和計算3:圧力を常圧に設定して、セルのyz断面積を固定した状態で、再度NPT計算を実施する。緩和計算4:セルの固定を解除して、圧力を常圧に設定して、NVT計算を実施する。以上により、繊維材料と接着液の界面モデルが作成されるステップS18。[繊維材料と接着液の界面モデルの評価]以上により作成された繊維材料と接着液の界面モデルの評価は、界面モデルの厚さ方向x軸方向。尚、繊維材料は、z軸方向に配向。のゴムまたは樹脂、あるいはこれらの混合物の分布状態を確認することにより行う。[繊維材料と接着液の界面モデルの剥離特性予測方法]次に、本実施形態に係るア繊維材料と接着液の界面の剥離特性予測方法の一例について、図4に基づいて説明する。図4は、本実施形態に係る繊維材料と接着液の界面の剥離特性予測方法の手順を示すフローチャート図である。実施形態に係る繊維材料と接着液の界面モデルの剥離特性予測方法では、以上により作成された繊維材料と接着液の界面モデルに対して、1軸伸張計算を実施し、界面の剥離挙動を解析することにより予測を行う。図4に示すように、界面の剥離挙動の解析するため、上述の本実施形態に係る繊維材料と接着液の界面モデル作成方法で作成された繊維材料と接着液の界面モデルに対して、所定の伸張速度,ポアソン比,NVTアンサンブルの条件で伸張させ、剥離シミュレーションを実施するS21。そして、剥離シミュレーション結果から得られた繊維材料と接着液の界面モデルの応力-変位曲線に基づいて、剥離応力の最大値を示す変位量を算出するS22。また、剥離シミュレーション結果から得られた繊維材料と接着液の界面モデルの応力-変位曲線により与えられる変位に対して応力を積分することにより、剥離エネルギーを算出するS23。更に、図4には記載していないが、S21における剥離シミュレーションの過程における伸張挙動を示すスナップショットにより、繊維材料と接着液の界面モデルの剥離解析を行っても良い。また、S22における最大剥離応力の算出とS23の剥離エネルギーの算出は、両方を行っても良いし、いずれかのみを行っても良い。尚、ゴムまたは樹脂、あるいはこれらの混合物の分布状態、ゴムまたは樹脂、あるいはこれらの混合物の分子量が異なる繊維材料と接着液の界面モデルを作成して剥離シミュレーションを実施して、繊維材料と接着液の界面の剥離挙動を比較することが可能である他、様々な材料を用いた繊維材料と接着液に基づいて、作成された繊維材料と接着液の界面モデルに基づいて、剥離シミュレーションを実施することで繊維材料と接着液の界面の剥離挙動を解析することができる。[繊維材料モデル作成装置または繊維材料と接着液の界面モデル作成装置または繊維材料と接着液の界面の剥離特性予測装置]次に、本実施形態に係る繊維材料モデル作成装置または繊維材料と接着液の界面モデル作成装置または繊維材料と接着液の界面の剥離特性予測装置の一例について、図5~図8に基づいて説明する。図5は、本実施形態に係る繊維材料モデル作成装置または繊維材料と接着液の界面モデル作成装置または繊維材料と接着液の界面モデルの剥離解析装置として使用されるコンピュータのハードウェア構成の一例を示すブロック図であり、図6は、本実施形態に係る繊維材料モデル作成装置の機能構成を示すブロック図である。図7は、本実施形態に係る繊維材料と接着液の界面モデル作成装置の機能構成を示すブロック図である。図8は、本実施形態に係る繊維材料と接着液の界面の剥離特性予測装置の機能構成を示すブロック図である。図5に示すように、本実施形態に係る繊維材料モデル作成装置100または繊維材料と接着液の界面モデル作成装置200または繊維材料と接着液の界面の剥離特性予測装置300は、演算部401と、記憶部402~404と、入力装置405と、出力装置406とがシステムバス407に接続されて構成される計算機例えば汎用のパーソナルコンピュータ等400上に実装される。かかる計算機400には、図5に示されているCPU401、ROM402、RAM403、ハードディスクなどの補助記憶装置404、キーボードやマウスなどの入力装置405、ディスプレイやプリンタなどの出力装置406などのハードウェアが収納されており、ハードディスクには、プログラムこのプログラムは、リムーバブルな記憶媒体に記録しておくことにより、様々なコンピュータにインストールすることが可能であるを含む各種のソフトウェアが記録されている。そして、これらのハードウェア及びソフトウェアが組み合わされることによって、上述の各部が構築されている。まず、本実施形態に係る繊維材料モデル作成装置100について、図6に基づいて説明する。図6に示すように、本実施形態に係る繊維材料モデル作成装置100は、モノマーモデル作成手段101と、重合体モデル作成手段102と、ランダム配置モデル作成手段103と、未配向モデル作成手段104と、第1の配向処理手段105と、第2の配向処理手段106と、繊維材料モデル作成手段107と、から構成されている。モノマーモデル作成手段101は、解析対象となる繊維材料について、原子間の相互作用を厳密に再現したモノマーモデルを、分子動力学計算により作成するためのものである。より詳細には、本実施形態に係る繊維材料モデル作成方法のステップS1の内容と同様であり、その説明を省略する。重合体モデル作成手段102は、モノマーモデル作成手段101で作成したモノマーモデルを所定の重合度で重合させて単位モデルである重合体モデルを作成するためのものである。より詳細には、本実施形態に係る繊維材料モデル作成方法のステップS2の内容と同様であり、その説明を省略する。ランダム配置モデル作成手段103は、重合体モデル作成手段102で作成した重合体モデルの複数本を、予め定めたセル空間内のランダムな位置に配置したランダム配置モデルを作成するためのものである。より詳細には、本実施形態に係る繊維材料モデル作成方法のステップS3の内容と同様であり、その説明を省略する。未配向モデル作成手段104は、ランダム配置モデル作成手段103で作成したランダム配置モデルを、分子動力学計算により構造緩和させて、温度,圧力等が指定した動的な計算の中でエネルギーが最小となる原子配置を計算し、未配向のモデルを作成するためのものである。より詳細には、本実施形態に係る繊維材料モデル作成方法のステップS4、5の内容と同様であり、その説明を省略する。第1の配向処理手段105は、未配向モデル作成手段104で作成された未配向のモデルに対し、セルの境界部分端部に歪を与えるLees-Edwards法等を用いて、モデルの1軸方向にせん断エネルギーを与え、分子鎖を1方向z軸方向=配向軸に引き伸ばす第1の配向処理を行うためのものである。より詳細には、本実施形態に係る繊維材料モデル作成方法のステップS6の内容と同様であり、その説明を省略する。第2の配向処理手段106は、第1の配向処理手段105で第1の配向処理を行った分子鎖の中から伸びきり状態の1本の分子鎖を選択し、この分子鎖を複写して同じ伸張状態の複数本の分子鎖を作成し、それぞれの分子鎖を繊維材料の結晶構造に対応した座標に配置したモデルを作成する第2の配向処理を行うためのものである。より詳細には、本実施形態に係る繊維材料モデル作成方法のステップS7の内容と同様であり、その説明を省略する。繊維材料モデル作成手段107は、第2の配向処理手段106で行った第2の配向処理直後のモデルに対して、構造緩和計算を行い、z軸方向に配向した繊維材料モデルを作成するためのものである。より詳細には、本実施形態に係る繊維材料モデル作成方法のステップS8、9の内容と同様であり、その説明を省略する。次に、本実施形態に係る繊維材料と接着液の界面モデル作成装置200について、図7に基づいて説明する。図7に示すように、本実施形態に係る繊維材料と接着液の界面モデル作成装置200は、モノマーモデル作成手段201と、重合体モデル作成手段202と、ランダム配置モデル作成手段203と、接着液モデル作成手段204と、接着液モデル及び繊維材料モデル配置手段205と、界面モデル作成手段206と、から構成されている。モノマーモデル作成手段201は、接着液を構成するゴムまたは樹脂、あるいはこれらのモノマーモデルを、それぞれ作成するためのものである。より詳細には、本実施形態に係る繊維材料と接着液の界面モデル作成方法のステップS11の内容と同様であり、その説明を省略する。重合体モデル作成手段202は、モノマーモデル作成手段201で作成された接着液を構成するゴムまたは樹脂、あるいはこれらのモノマーモデルに基づいて、ゴムまたは樹脂、あるいはこれらの重合体モデルを、それぞれ作成するためのものである。より詳細には、本実施形態に係る繊維材料と接着液の界面モデル作成方法のステップS12の内容と同様であり、その説明を省略する。ランダム配置モデル作成手段203は、重合体モデル作成手段202でそれぞれ作成したゴムまたは樹脂、あるいはこれらの重合体モデルを、所定の混合比混合比として、例えば、重量比や分子数比となるように予め定めたセル空間にランダムに配置して、ランダム配置モデルを作成するためのものである。より詳細には、本実施形態に係る繊維材料と接着液の界面モデル作成方法のステップS13の内容と同様であり、その説明を省略する。接着液モデル作成手段204は、ランダム配置モデル作成手段203で作成されたランダム配置モデルを、分子動力学計算により構造緩和させて、温度,圧力等が指定した動的な計算の中でエネルギーが最小となる原子配置を計算し、未配向の接着液モデルを作成するためのものである。より詳細には、本実施形態に係る繊維材料と接着液の界面モデル作成方法のステップS14、15の内容と同様であり、その説明を省略する。接着液モデル及び繊維材料モデル配置手段205は、接着液モデル作成手段204で作成された接着液モデルを、z軸方向に配向した繊維材料モデル上述した本実施形態に係る繊維材料モデル作成装置100において作成の下部に配置し、界面モデルを作成するためのものである。より詳細には、本実施形態に係る繊維材料と接着液の界面モデル作成方法のステップS16の内容と同様であり、その説明を省略する。界面モデル作成手段206は、接着液モデル及び繊維材料モデル配置手段205で作成された界面モデルに対して、分子鎖を熱平衡化させるために、緩和計算を実施し、繊維材料と接着液の界面モデルを作成するためのものである。より詳細には、本実施形態に係る繊維材料と接着液の界面モデル作成方法のステップS17、18の内容と同様であり、その説明を省略する。次に、本実施形態に係る繊維材料と接着液の界面の剥離特性予測装置300について、図8に基づいて説明する。図8に示すように、本実施形態に係る繊維材料と接着液の界面の剥離特性予測装置300は、1軸伸張手段301と、最大剥離応力算出手段302と、剥離エネルギー算出手段303と、から構成されている。1軸伸張手段301は、界面の剥離挙動の解析するため、上述の本実施形態に係る繊維材料と接着液の界面モデル作成装置200で作成された繊維材料と接着液の界面モデルに対して、所定の伸張速度,ポアソン比,NVTアンサンブルの条件で伸張させ、剥離シミュレーションを実施するためのものである。より詳細には、本実施形態に係る繊維材料と接着液の界面モデルの剥離解析方法のステップS21の内容と同様であり、その説明を省略する。ここで、1軸伸張手段301において、剥離シミュレーションの過程における伸張挙動を示すスナップショットにより、繊維材料と接着液の界面モデルの剥離解析を行っても良い。最大剥離応力算出手段302は、剥離シミュレーション結果から得られた繊維材料と接着液の界面モデルの応力-変位曲線に基づいて、剥離応力の最大値を示す変位量を算出するためのものである。より詳細には、本実施形態に係る繊維材料と接着液の界面モデルの剥離解析方法のステップS22の内容と同様であり、その説明を省略する。剥離エネルギー算出手段303は、また、剥離シミュレーション結果から得られた繊維材料と接着液の界面モデルの応力-変位曲線により与えられる変位に対して応力を積分することにより、剥離エネルギーを算出するためのものである。より詳細には、本実施形態に係る繊維材料と接着液の界面モデルの剥離解析方法のステップS23の内容と同様であり、その説明を省略する。また、最大剥離応力算出手段302及び剥離エネルギー算出手段303は、両方が繊維材料と接着液の界面の剥離特性予測装置300に構成されるものでも良いし、いずれかのみ繊維材料と接着液の界面の剥離特性予測装置300に構成されるものであっても良い。このように、本実施形態に係る繊維材料モデル作成方法と繊維材料モデル作成装置、及び、繊維材料と接着液の界面モデル作成方法と繊維材料と接着液の界面モデル作成装置、並びに、繊維材料と接着液の界面の剥離特性予測方法と繊維材料と接着液の界面の剥離特性予測装置によれば、作成した1軸方向に配向した繊維材料のフルアトム分子動力学モデルとゴムまたは樹脂、あるいはこれらの混合物が混在した接着液フルアトム分子動力学モデルに基づいて繊維材料と接着液の界面モデルを作成することができる。そして、作成された繊維材料と接着液の界面モデルを使って剥離解析を行うことにより、界面における材料の配合と剥離特性の関係とを予測することができ、ひいては、実験では把握しがたい繊維材料と接着液の界面におけるナノミクロ領域の挙動と剥離特性の関係を予測することができる。次に、本実施形態の繊維材料モデル作成方法、及び、繊維材料と接着液の界面モデル作成方法、並びに、繊維材料と接着液の界面モデルの剥離解析方法を実際に用いた場合の実施例について説明する。本実施例では、トワロン,ケブラーなどに代表される、パラ系アラミド繊維ポリパラフェニレンテレフタラミドを解析対象となる繊維材料とした。尚、本実施形態において、解析対象とする繊維材料は、パラ系アラミド繊維に限らず、メタ系アラミド繊維にも適用することができる他、有機繊維アラミド繊維、ナイロン繊維のようなポリアミド繊維、ポリエステル繊維等と無機繊維カラス繊維、炭素繊維等も含まれる。まず、上述した本実施形態に係る繊維材料モデル作成方法及び本実施形態に係る繊維材料モデル作成装置を用いて、アラミド繊維モデルを作成する。本実施例では、分子動力学解析ソフトJ-OCTAを用い、周期境界条件を採用して以下のパラ系アラミド繊維のフルアトムモデルを作成した。尚、上述した本実施形態の繊維材料モデル作成方法図1及び図2に基づいて、本実施例について以下で説明する。最初に、原子間の相互作用を厳密に再現したモノマーモデルを作成しステップS1、次に、10個のモノマーモデルを重合させアラミド分子としたステップS2。尚、アラミド繊維を構成する原子間の結合,角度,二面角ポテンシャル、およびファンデルワールスポテンシャルを設定するための力場はDREIDINGとGAFFを検討した。そして、分子鎖を形成する各原子が保有する電荷は、汎密度関数法B3LYP/631Gdで求めたモノマーの電荷を用い、重合後の分子鎖の電荷の総和がゼロとなるように任意の原子が保有する電荷を調整した。更に、立体規則性の影響を調査するために、立体規則性がアイソタクチック,シンジオタクチック,アタクチックとなる3種類の構造を作成した。そして、上記の条件で作成したアラミド分子6本をセル内のランダムな位置に設置したモデルを作成したステップS3。さらに、分子鎖を熱平衡に近い状態にするために、高圧100MPa条件でNPT-AndersenNoseHoover計算を実施した緩和計算1:ステップS401。その後、圧力を常圧0.1MPaに設定し、再度NPT計算を実施した緩和計算2:ステップS402。最後に、圧力が常圧で、NVT-AndersenNoseHoover計算を実施して緩和計算3:ステップS403、未配向のアラミドモデルが作成されるステップS5。尚、緩和計算における設定温度は300K、ステップ間隔は1fsec、ステップ数は10万ステップ100psecとし、周期的境界条件を用いて計算した。次に、分子鎖の第一の配向処理として、未配向のアラミドモデルに対しセルの境界部分端部に歪を与えるLees-Edwards法を使ってモデルのz軸配向軸方向にせん断エネルギーを与えずり又は1軸伸張、分子鎖を1方向へ引き延ばしたステップS6。なお、シェア0.05、初期歪1、設定温度は300K、ステップ間隔は1fsec、ステップ数は10万ステップ100psecとし、周期的境界条件を用いて計算した。そして、分子鎖の第二の配向処理として、配向させた分子の中から伸びきり状態の1本鎖を選択し、この分子を複写して同じ伸張状態の5本のアラミド分子を作成し、それぞれの分子鎖をアラミドの結晶構造に対応した座標に配置したステップS7。ここで、アラミド分子の結晶の格子定数はa=7.87Å,b=5.18Å,c=12.9Å,γ=90と考えられているため、配向方向z軸方向c軸に垂直な座標がa=7.87Å,b=5.18Åとなるようにセル内に5本の分子鎖を配置した。そして、緩和計算として、配向処理直後のモデルを、高圧100MPa条件で、NPT-ParrinelloRahmanNoseHoover計算により構造緩和させた。配向後の分子鎖はz方向に引き伸ばされており、xyz軸の各方向に異方応力が作用しているため、セル長さのxy方向xy断面積だけを固定して、高圧100MPa条件でNPT-AndersenNoseHoover計算を実施して緩和後緩和計算1:ステップS801、z軸方向だけを固定して、再度、高圧100MPa条件でNPT-AndersenNoseHoover計算を実施して緩和した緩和計算2:ステップS802。次に、セルの固定を解除して、圧力を常圧に設定して、再度NPT計算を実施して緩和した緩和計算3:ステップS803。最後に、圧力が常圧で、NVT計算を実施する緩和計算4:ステップS804。尚、緩和計算における設定温度は300K、ステップ間隔は1fsec、ステップ数は10万ステップ100psecとし、周期的境界条件を用いて計算した。これにより、本実施例のアラミド繊維モデルが作成されるステップS9表1に本実施例により作成したアラミド繊維モデルの比重と配向度を示す。ここで、現実のアラミド繊維の比重の文献値は、約1.4であり、配向度の文献値は、X線回折法で求めたものである。そして、アラミド繊維モデルの比重は現実のアラミド繊維比重約1.4よりも小さい値となった。モデルR4の比重は1.331となり、最も実材料に近い値となった。一方、配向度は立体配座がランダムなモデルでは0.9以上を示した。図9に、配向後のアラミド繊維モデルモデルR4を示す。図9左側に示すyz平面と図9右側に示すyx面から見た分子の配置が明らかに異なっており、z軸方向にアラミド分子が配向している様子が確認できる。ここで、ケブラーなどのパラ系アラミド繊維の水分率は約7である.現実のアラミド繊維の水分率sはアラミド分子の重量Aと水の重量Bを用いて以下の式で表される。s=B/AB1また、このときの繊維の比重ρexは、式1の関係と体積Vを用いて、以下の式で表される。ρex=AB/V=1s/1-sA/V2一方、吸水による繊維の体積変化が微小であると仮定すると、シミュレーションにおけるアラミド繊維の比重ρsmは以下の式となる。ρsm=A/V3吸水率が7s=0.07のときの比重ρexは、式2,3より、以下のようになる.ρex=AB/V=10.07/0.93×A/V=1.075×ρsmこの式に配向後のアラミド繊維のシミュレーションから得られた比重ρsm=1.331を代入すると、ρex=1.431となる。以上のことから、シミュレーションモデルには水の吸収が考慮されていないため、実材料に比べ比重が小さくなったものと推定される。以上から、本実施例では、パラ系アラミド繊維のフルアトム分子動力学モデルを作成する際に、力場や配向処理の方法を検討することで、現実の繊維材料に近い比重と配向度を有する配向後のアラミド繊維モデルR4を作成することができた。ここで、R4のセルサイズは2.46585nmZ軸方向×2.44005nmY軸方向×2.49245nmX軸方向である。本実施例で作成されたアラミド繊維モデルの弾性率を評価する。作成されたアラミド繊維モデルR4に対して、伸張速度0.5m/s0.005Å/ps、ポアソン比0.3、NVTアンサンブルの条件で伸張シミュレーションを実施した。図10に、伸張シミュレーション結果から得られたアラミド繊維モデルの応力-ひずみ曲線を示す。図10では、ひずみ量0.025付近に降伏点が確認できる。そして、図10の線形領域strain=0~0.017,ひずみ量0~1.7の傾きから、アラミド繊維モデルの弾性率は、47.409GPaと算出された。また、図10の線形領域の応力の最大値から降伏応力は1.44GPaと算出される。表2に、ケブラーフィラメントの物性値とモデルの物性値を示す。尚、表2において、伸度=伸び/元の長さ×100としている。本実施例において、降伏点において材料が破断するものと仮定し、降伏応力を引張破断とすると、引張強度は表2の約半分の値、破断時の伸度図10において、線形領域の応力が最大値となる伸度は、2.96と見積もられ、表2にけるアラミド繊維の破断伸度2.4~4.4とほぼ一致する結果となった。シミュレーション結果の弾性率はアラミド繊維の理論弾性率182GPa日本ポリイミド芳香族系高分子研究会編新訂最新ポリイミド-基礎と応用-1292010の約1/3の値となっている.作成したアラミド繊維モデルは、分子鎖が5本なので1本鎖で計算する理論弾性率よりも低くなるものと考えられる。従って、表2に示したフィラメントの物性値と比較するのが妥当と考えると、アラミド繊維モデルの弾性率は、ケブラー117フィラメントの弾性率と同程度、引張強度は1/2程度となる。次に、作成されたアラミド繊維モデルを用いて上述した本実施形態に係る繊維材料と接着液の界面モデル作成方法及び本実施形態に係る繊維材料と接着液の界面モデル作成装置図3を用いて、繊維材料と接着液の界面モデルを作成する。本実施例では、解析対象となる繊維材料を上述の通りアラミド繊維とし、接着液をRF樹脂とラテックスから構成されるRFLとして、分子動力学解析ソフトJ-OCTAを用い、周期境界条件を採用して以下のRFLのフルアトムモデル、及び、アラミド繊維とRFLの界面モデルのフルアトムモデルを作成した。尚、本実施形態において、解析対象とする接着液は、RF樹脂とラテックスから構成されるRFLに限らず、エポキシ樹脂、ブロックドイソシアネート等を含むRFLにも適用することができる他、RFL以外にエポキシ樹脂単体、またゴム組成物が含まれ、ゴム組成物としては、クロロプレンゴム、EPDMエチレンプロピレンジエンモノマー等のゴム組成物を溶剤メチルエチルケトン等に溶かしたゴム糊も含まれる。尚、上述した本実施形態に係る繊維材料と接着液の界面モデル作成方法図3に基づいて、本実施例について以下で説明する。まず、RF樹脂とラテックスのモノマーモデルを作成し、RF樹脂とラテックスのモノマーモデルを所定の重合度で重合したステップS11~12。尚、RF樹脂の重合体モデルは、RF樹脂のモノマーモデルを所定の重合度複数量体で重合するものとする。また、ラテックスの重合体モデルは、解析対象とするラテックスのモノマーの配合比に基づいて算出した重合比と、1分子の分子量に基づいて算出する。本実施例では、RF樹脂は7量体とし、また、ラテックスはVpラテックスであるNipole2518FS日本ゼオン製をモデル化し、各モノマーの配合比がスチレン/ブタジエン/ビニルピリジン=70/15/15wtであることから日本ゼオンゴムと繊維の接着剤組成物特願昭60-13869号公報参照、本実施例に係るRFLモデルの重合比はスチレン/ブタジエン/ビニルピリジン=3/24/3、1分子の分子量は1919455ビーズとした。尚、RFLモデルにおけるスチレンS,ブタジエンB,ビニルピリジンPの配列順序はランダムとした結果、BSPBBBBBBBSBBBBBPBBBBSBBBBBPBBの順序となった。そして、分子を構成する原子間の結合、角度、二面角ポテンシャル、およびファンデルワールスポテンシャルを設定するための力場は、水素結合を再現できるGAFFを使用した。そして、分子鎖を形成する各原子が保有する電荷は、汎密度関数法B3LYP/631Gdで求めたモノマーの電荷を用い、重合後の分子鎖の電荷の総和がゼロとなるように任意の原子が保有する電荷を調整した。上記の条件でそれぞれ作成したRF分子とVp分子を、所定の混合比本実施例では、重量比で1:1分子数比で4分子:10分子の割合となるように1つセルに配置した2860ビーズステップS13。さらに、分子鎖を熱平衡に近い状態にするために高圧100MPa条件でNPT-AndersenNoseHoover計算を実施した緩和計算1。その後、圧力を常圧0.1MPaに設定し、再度NPT計算を実施した緩和計算2。最後に、圧力が常圧で、NVT-AndersenNoseHoover計算を実施して緩和計算3、RFLモデルが作成されるステップS14~15。尚、緩和計算における設定温度は300K、ステップ間隔は1fsec、ステップ数は10万ステップ100psecとし、周期的境界条件を用いて計算した。次に、作成したRFLモデルをz軸方向に配向させた上記で作成したアラミド繊維モデルの下部に配置し、アラミド繊維とRFLの界面モデルを作成するステップS16~18。本実施例では、分子鎖の配置はアラミド繊維モデルの配向軸であるx軸と垂直な軸であるy軸を基準として、これらのセルの長さを固定し、RFLモデルを再配置したステップS16。このときのセルサイズは2.446585nmz軸方向×2.440055nmy軸方向×11.992445nmx軸方向である。尚、界面モデルの初期配置は、J-OCTAシステムモデラを使用した。さらに、分子鎖を熱平衡化させるために、以下のステップ緩和計算1~5で緩和計算を実施したステップS17。尚、緩和計算における設定温度は300K、ステップ間隔は1fsecとした。緩和計算1:接着液分子だけを熱平衡化させるために、繊維材料モデルの分子鎖のxyz座標を固定し、セルのyz軸断面積を固定して、高圧条件100MPaで、NPT計算を実施した。なお、セルのyz軸は固定し、10万ステップ100psec、x軸の境界は反射境界条件、電荷計算はcutoffCoulombを選択した。尚、本実施例において、アラミド繊維は、熱平衡化済みである。緩和計算2:繊維材料モデルの分子鎖の配向を維持し、分子鎖とx軸の境界面間の空隙を削除するために、分子全体をx軸方向に-0.5nm移動させ、隙間をつめる。高圧条件100MPaで、分子鎖のxyz座標の固定を解除し、セルのyz断面積を固定した状態で、NPT計算を10万ステップ100psec実施した。緩和計算3:圧力を常圧0.1MPaに設定して、セルのyz断面積を固定した状態で、再度NPT計算を10万ステップ100psec実施した。緩和計算4:セルの固定を解除して、圧力を常圧0.1MPaに設定して、NVT計算を100万ステップ1nsec実施した。以上により、アラミド繊維とRFLの界面モデルを作成したステップS18。このときのx軸方向のセルサイズは、6.86914nm、比重は1.1354g/cm3、計算時間は137時間であった。以上により作成された本実施例に係るアラミド繊維とRFLの界面モデルを、図11に示す。図11において、x軸方向左右方向が界面の厚さ方向を示している。図11から、アラミド繊維はz軸方向上下方向に配向しており、アラミド繊維表面付近にRF樹脂が多く分布していることが確認できる。図12に、本実施例に係るアラミド繊維とRFLの界面モデルにおける厚さ方向深さ方向のアラミド繊維、RF樹脂、Vpラッテクスの分布を示すものであり、左側はモデル1を、右側はモデル2を示している。図12の横軸はx軸座標厚さ方向、縦軸は任意の厚さに存在する各原子の数である。図12左側のモデル1はアラミド繊維表面Depth=25Å付近にRF樹脂が多く分布しており、RF樹脂とVpラテックスが混合状態であることが確認できる。一方、図12右側のモデル2ではアラミド繊維表面付近にVpラッテクスが多く分布しており、RF樹脂とVpラテックスが混合状態であることが確認できる。尚、モデル1のセルサイズは2.446585nmz軸方向×2.440055nmy軸方向×6.83731nmx軸方向、モデル2のセルサイズは2.446585nmz軸方向×2.440055nmy軸方向×6.89904nmz軸方向であり、アラミド繊維とRFLの厚さ方向の長さは、モデル1が2.5nm、モデル2が4.5nmと算出された。次に、作成されたアラミド繊維とRFLの界面モデルを用いて上述した本実施形態に係る繊維材料と接着液の界面の剥離特性予測方法及び本実施形態に係る繊維材料と接着液の界面の剥離特性予測装置を用いて、アラミド繊維とRFLの界面の剥離特性を予測する。尚、上述した本実施形態に係る繊維材料と接着液の界面の剥離特性予測方法図4に基づいて、本実施例について以下で説明する。本実施例に係るアラミド繊維とRFLの界面モデルモデル1,2に対して、1軸伸張計算を実施し界面の剥離挙動を確認したS21。尚、伸張計算はNVTアンサンブル、伸張速度25m/s、ポアソン比0の条件で伸張し、モデルの両端0.2nmの原子のy座標、z座標伸張方向に垂直な面の座標を固定しx軸伸張方向の座標を伸張速度と同じ速度で移動させた。本実施例に係るアラミド繊維とRFLの界面モデルの剥離挙動について、図13及び図14に基づいて、説明する。図13は、実施例に係るアラミド繊維とRFLの界面モデルの剥離シミュレーション結果から得られた応力-変位曲線を示すものであり、左側はモデル1を、右側はモデル2を示している。図13に示すように、剥離応力は変位量約1nm付近で最大値を示したS22。また、図14は、実施例に係るアラミド繊維とRFLの界面モデルの剥離シミュレーションの剥離過程における伸張挙動を示すものであり、左側はモデル1を、右側はモデル2を示している。図14に示すように、両モデルの剥離形態はRFLおそらくVpラテックスが破壊して剥離接着強度は、アラミド繊維とRFLの界面>RFLしている。尚、今回モデル化したVpラッテクスは未架橋状態である。また、モデル1はモデル2より剥離後のアラミド繊維表面へのRF樹脂の付着量が多いことが分かる。これは、RF樹脂とアラミドの相互作用の強さが大きいためと考えられる。さらに、図13の応力-変位曲線の積分により剥離エネルギーを算出すると、モデル1の剥離エネルギーは、12.55kJ/m2モデル2の剥離エネルギーは、12.38kJ/m2となり、アラミド繊維表面にRF樹脂が多く分布するモデル1の方が剥離するために、僅かに大きなエネルギーが必要なことが確認されるS23。本実施例においては、RF樹脂の分布状態の異なるアラミド繊維とRFLの界面モデルモデル1,2を作成し、剥離シミュレーションを実施した。そして、アラミド繊維と相互作用の強いRF樹脂がアラミド繊維表面に多く分布したモデル2の剥離エネルギーはRF樹脂の分布が小さいモデル1に比べ大きな剥離エネルギーを示すことが確認された。このように、本実施形態の繊維材料と接着液の界面モデル作成方法及び剥離シミュレーションを様々な材料に適用することで、繊維材料と接着液の界面の剥離挙動を解析できると考えられる。例えば、RF樹脂の分布、RF樹脂の分子量、Vpラテックス分子量5本鎖から1本鎖へ変更が異なるモデルを作成し、その剥離挙動を解析することが考えられる。"