(57)【特許請求の範囲】
【請求項１】
初期形状及び構成粒子を含み、所定温度及び圧力を含む解析条件のもとで分子動力学計算により平衡化されたゴムモデルデータを取得するステップと、
第１軸の変形及び第２軸の変形を少なくとも含んだ変形後の形状になるように、前記平衡化したゴムモデルを構成する粒子の位置を移動させ、ゴムモデルを変形させるステップと、
変形後のゴムモデルを分子動力学計算により平衡化するステップと、
平衡状態にある変形後のモデルに基づき第１軸から第３軸の主応力（σ１，σ２，σ３）を算出するステップと、
第１軸から第３軸の主応力（σ１，σ２，σ３）および第１軸から第３軸の伸長比（λ１，λ２，λ３）、並びに式（１）〜（５）を用いて、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出するステップと、
を含む、ゴムの二軸引張をシミュレーションする方法。
【数１】
0006568782.tif 000012
【数２】
0006568782.tif 000013
Ｉ１＝λ１２＋λ２２＋λ３２      …（４）
Ｉ２＝λ１２λ２２＋λ２２λ３２＋λ３２λ１２      …（５）
【請求項２】
初期形状から目標となる最終形状までの前記ゴムモデルの変形を、変形１回あたりの変形量が所定量以下になるように複数段階に分け、前記ゴムモデルを変形させるステップと前記変形後のモデルを平衡化するステップとを繰り返し実行する、請求項１に記載の方法。
【請求項３】
初期形状及び構成粒子を含み、所定温度及び圧力を含む解析条件のもとで分子動力学計算により平衡化されたゴムモデルデータを取得する初期形状ゴムモデル取得部と、
第１軸の変形及び第２軸の変形を少なくとも含んだ変形後の形状になるように、前記平衡化したゴムモデルを構成する粒子の位置を移動させ、ゴムモデルを変形させる変形実行部と、
変形後のゴムモデルを分子動力学計算により平衡化する変形後モデル平衡化演算部と、
平衡状態にある変形後のモデルに基づき第１軸から第３軸の主応力（σ１，σ２，σ３）を算出する主応力算出部と、
第１軸から第３軸の主応力（σ１，σ２，σ３）および第１軸から第３軸の伸長比（λ１，λ２，λ３）、並びに式（１）〜（５）を用いて、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出する微分値算出部と、
を備える、ゴムの二軸引張をシミュレーションする装置。
【数３】
0006568782.tif 000014
【数４】
0006568782.tif 000015
Ｉ１＝λ１２＋λ２２＋λ３２      …（４）
Ｉ２＝λ１２λ２２＋λ２２λ３２＋λ３２λ１２      …（５）
【請求項４】
初期形状から目標となる最終形状までの前記ゴムモデルの変形を、変形１回あたりの変形量が所定量以下になるように複数段階に分け、前記変形実行部による変形処理と前記変形後モデル平衡化演算部による平衡化処理とを繰り返し実行する、請求項３に記載の装置。
【請求項５】
請求項１又は２に記載の方法をコンピュータに実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、分子動力学を用いた、ゴムの二軸引張シミュレーション方法、装置及びプログラムに関する。
【背景技術】
【０００２】
ゴムの力学特性を調べるために、一軸伸長試験が知られている。この試験では、一軸に沿ったゴムの伸長比λ及び引張応力σの関係を測定でき、ゴム弾性を知る上で標準的な測定方法として広く知られている。
【０００３】
ゴム変形の特徴として、変形による体積変化がないこと、及び、応力歪みの関係は方向によらないこと（等方性）が挙げられる。歪みエネルギー密度Ｗは、Ｗ＝Ｗ（Ｉ１，Ｉ２）で表される。歪み不変量Ｉ１，Ｉ２，Ｉ３は次の式（４）〜（６）で表される。
Ｉ１＝λ１２＋λ２２＋λ３２                  …（４）
Ｉ２＝λ１２λ２２＋λ２２λ３２＋λ３２λ１２    …（５）
Ｉ３＝λ１２λ２２λ３２＝１                  …（６）
それぞれ直交する第１軸、第２軸、第３軸を主軸とした場合に、λ１は第１軸の伸長比、λ２は第２軸の伸長比、λ３は第３軸の伸長比を示す。ゴムの特性である体積変化がないことによりＩ３＝１である。
【０００４】
しかしながら、上記に示すような一般的な変形に対する応力は、一軸伸長試験では分からないため、二軸引張試験が必要となる。二軸引張試験は、第１軸に沿ったゴムの伸長比λ１及び引張応力σ１と、第２軸に沿ったゴムの伸長比λ２及び引張応力σ２とを測定でき、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を得ることができる。二軸引張試験やその基礎理論については、非特許文献１〜３に開示されている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】川端季雄, “ゴムの力学 [III].”, 日本ゴム協会誌 53(4), p.208 (1980).
【非特許文献２】川端季雄, “ゴムの力学 [IV].”, 日本ゴム協会誌 53(7), p.408 (1980).
【非特許文献３】川端季雄, “ゴムの力学 [V].”, 日本ゴム協会誌 53(8), p.479 (1980).
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、二軸引張試験にて得られる結果を、分子動力学を用いたシミュレーションにて算出する方法については、知られていない。
【０００７】
本開示は、このような課題に着目してなされたものであって、その目的は、分子動力学を用いたゴムの二軸引張シミュレーション方法、装置及びプログラムを提供することである。
【課題を解決するための手段】
【０００８】
本開示は、上記目的を達成するために、次のような手段を講じている。
【０００９】
すなわち、本開示のゴムの二軸引張をシミュレーションする方法は、
初期形状及び構成粒子を含み、所定温度及び圧力を含む解析条件のもとで分子動力学計算により平衡化されたゴムモデルデータを取得するステップと、
第１軸の変形及び第２軸の変形を少なくとも含んだ変形後の形状になるように、前記平衡化したゴムモデルを構成する粒子の位置を移動させ、ゴムモデルを変形させるステップと、
変形後のゴムモデルを分子動力学計算により平衡化するステップと、
平衡状態にある変形後のモデルに基づき第１軸から第３軸の主応力（σ１，σ２，σ３）を算出するステップと、
第１軸から第３軸の主応力（σ１，σ２，σ３）および第１軸から第３軸の伸長比（λ１，λ２，λ３）、並びに式（１）〜（５）を用いて、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出するステップと、
を含む。
【００１０】
本開示のゴムの二軸引張をシミュレーションする装置は、
初期形状及び構成粒子を含み、所定温度及び圧力を含む解析条件のもとで分子動力学計算により平衡化されたゴムモデルデータを取得する初期形状ゴムモデル取得部と、
第１軸の変形及び第２軸の変形を少なくとも含んだ変形後の形状になるように、前記平衡化したゴムモデルを構成する粒子の位置を移動させ、ゴムモデルを変形させる変形実行部と、
変形後のゴムモデルを分子動力学計算により平衡化する変形後モデル平衡化演算部と、
平衡状態にある変形後のモデルに基づき第１軸から第３軸の主応力（σ１，σ２，σ３）を算出する主応力算出部と、
第１軸から第３軸の主応力（σ１，σ２，σ３）および第１軸から第３軸の伸長比（λ１，λ２，λ３）、並びに式（１）〜（５）を用いて、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出する微分値算出部と、
を備える。
【００１１】
このようにすれば、従来では実現方法が開発されていなかった、ゴムの二軸引張シミュレーションを実行でき、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出することが可能となる。また、実試験では試験できる領域が限定されているが、シミュレーションであれば、その制限を超えた範囲にて二軸引張におけるゴムの挙動を知ることができる。
【図面の簡単な説明】
【００１２】
【図１】本開示のゴムの二軸引張をシミュレーションする装置を模式的に示すブロック図。
【図２】ゴムの二軸引張変形に関する説明図。
【図３】装置が実行する処理ルーチンを示すフローチャート。
【図４】シミュレーション結果を示す図。
【発明を実施するための形態】
【００１３】
以下、本開示の一実施形態を、図面を参照して説明する。
【００１４】
［ゴムの二軸引張をシミュレーションする装置］
装置１は、ゴムの二軸引張試験を、分子動力学計算によりコンピュータでシミュレーションし、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出する装置である。
【００１５】
図１に示すように、装置１は、初期設定部１０と、初期形状ゴムモデル取得部１１と、変形実行部１２と、変形後モデル平衡化演算部１３と、主応力算出部１４と、微分値算出部１５と、を有する。これら各部１０〜１５は、ＣＰＵ、メモリ、各種インターフェイス等を備えたパソコン等の情報処理装置において予め記憶されている図示しない処理ルーチンをＣＰＵが実行することによりソフトウェア及びハードウェアが協働して実現される。
【００１６】
図１に示す初期設定部１０は、キーボードやマウス等の既知の操作部を介してユーザからの操作を受け付け、解析対象となるゴムモデルに関するデータの設定、分子動力学計算に必要な解析条件、二軸伸長変形条件などの各種設定を実行し、これら設定値をメモリに記憶する。ゴムモデルに関するデータとしては、分子動力学計算により平衡化されたゴムモデルデータそのものでもよいし、ゴムモデルデータを生成するために必要となる高分子、架橋ビーズ、それらのポテンシャルなどのデータでもよい。二軸伸長変形条件としては、少なくとも第１軸及び第２軸に対する変形方向、その変形率（伸長比）が挙げられる。図２に示す一例では、第１軸に対して伸長させ、第２軸に対して縮小させ、第３軸は維持する単純せん断線形であり、最終形状Ｓ_finalにおける第１軸の伸長比λ１は３．０、第２軸の伸長比λ２は０．３３３・・、第３軸の伸長比λ３は１．０に設定されている。
【００１７】
初期形状ゴムモデル取得部１１は、初期形状及び構成粒子を含み、所定温度及び圧力を含む解析条件のもとで分子動力学計算により平衡化されたゴムモデルデータＭ１を取得する。本実施形態では、初期形状ゴムモデル取得部１１は、分子動力学計算を実行して平衡化されたゴムモデルデータを生成するが、これに限定されず、装置１外で予め生成され、メモリに記憶されたゴムモデルデータを取得するようにしてもよい。
【００１８】
ゴムモデルデータＭ１は、図２に示すように初期形状Ｓ_initを有し、ゴムモデルを構成する粒子が配置されたデータである。本実施形態において初期形状Ｓ_initは立方体であるがこれに限定されず、種々の形状に設定できる。ゴムモデルデータＭ１を構成する粒子は、高分子に架橋ビーズが反応して結合されたものである。ゴムモデルデータＭ１は、未架橋高分子の分子鎖に反応可能なビーズを設定し、架橋ビーズを投入し、全ての反応可能なビーズが架橋ビーズと結合される平衡状態になるまで分子動力学計算を演算することで得られる。本実施形態では、一例として、２００ビーズからなるKremer-Grest分子鎖モデル１００本に対して、各分子鎖に等間隔に１０個の反応可能ビーズを設定しておき、２００個の架橋ビーズを添加して、架橋反応させながら平衡化して生成したゴムモデルデータＭ１を使用している。未架橋高分子モデルの結合ポテンシャルには、ＦＥＮＥ−ＬＪ（レナードジョーンズ）が設定され、非結合ポテンシャルには、ＷＣＡ（斥力のみのＬＪポテンシャル）が設定されている。勿論、これは一例であって、その他の設定が可能である。
【００１９】
変形実行部１２は、第１軸の変形及び第２軸の変形を少なくとも含んだ変形後の形状になるように、平衡化したゴムモデルを構成する粒子の位置を移動させ、ゴムモデルを変形させる。例えば、図２に示すように、初期形状Ｓ_initから変形後の形状Ｓ_iに変形させる場合には、変形後の形状Ｓ_iにおける各座標軸の伸長比に応じて、初期形状Ｓ_initにおける粒子ｇ１の位置を移動させることにより、ゴムモデルに歪みを与えて変形することができる。
【００２０】
変形後モデル平衡化演算部１３は、変形後のゴムモデルを分子動力学計算により平衡化する。すなわち、平衡状態になるまで分子動力学計算を実行する。変形後のゴムモデルは、粒子の位置を強制的に移動させているために応力変動が不安定な状態となっているため、熱平衡状態になるまで、変形後のゴムモデルの分子動力学計算を繰り返し実行する。熱平衡状態であるかの判定方法として、例えば、ゴムモデルの温度と圧力がほぼ一定になれば熱平衡状態であると判定することが挙げられる。
【００２１】
主応力算出部１４は、平衡状態にある変形後のモデルに基づき第１軸から第３軸の主応力（σ１，σ２，σ３）を算出する。系全体（ゴムモデル全体）の第１軸から第３軸の主応力（σ１，σ２，σ３）の算出は、全粒子の位置及び速度に基づき算出できるが、公知であるために説明を省略する。タイムステップ毎に応力変動があるために、複数のタイムステップに亘り主応力（σ１，σ２，σ３）を算出しておき、その平均値を採用することが好ましい。なお、現実世界における二軸引張で計測される応力は、引張でゴムに作用する応力のみであるが、分子動力学計算において算出される主応力（σ１，σ２，σ３）には、静水圧成分（大気圧）が含まれている。
【００２２】
微分値算出部１５は、第１軸から第３軸の主応力（σ１，σ２，σ３）および第１軸から第３軸の伸長比（λ１，λ２，λ３）、並びに式（１）〜（５）を用いて、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出する。
【数１】
0006568782.tif 000002
【数２】
0006568782.tif 000003
Ｉ１＝λ１２＋λ２２＋λ３２      …（４）
Ｉ２＝λ１２λ２２＋λ２２λ３２＋λ３２λ１２      …（５）
【００２３】
上記式（１）〜（３）について、以下で説明する。
【００２４】
歪みエネルギー密度関数と主応力の関係式は、次の通りである。
0006568782.tif 000004
【００２５】
上記で述べた通り、分子動力学計算において算出される主応力（σ１，σ２，σ３）には、静水圧成分（大気圧）が含まれているので、これを除く必要がある。静水圧をキャンセルするための式は、次の通りである。
0006568782.tif 000005
【００２６】
この式を３行２例の行列式に変形すれば、次の式（６）となる。
0006568782.tif 000006
【００２７】
ここで、
0006568782.tif 000007
とすれば、式（６）は、次の式で表すことができる。
0006568782.tif 000008
【００２８】
ここで、Λの転置行列を用いて、正規化すれば、次の通りとなる。
0006568782.tif 000009
この式を変形すると、次のような２行２列の式となる。  
0006568782.tif 000010
これを、歪みエネルギー密度関数の歪み不変量に関する微分値を求める式に変形すれば、次の式（１）、（２）が得られる。
【数３】
0006568782.tif 000011
【００２９】
ところで、図２に示すように、初期形状Ｓ_initから最終形状Ｓ_finalまでゴムモデルＭ１を変形させるにあたり、変形実行部１２による１回の変形処理で初期形状Ｓ_initから最終形状Ｓ_finalに変形させてしまうと、急激に多大な応力が作用することになる。分子動力学計算では運動方程式を差分法で解くために、差分法の能力を超えて計算が破綻してしまい、変形後モデル平衡化演算部１３による演算が不能になるおそれがある。
【００３０】
そこで、図２に例示するように、初期形状Ｓ_initから目標となる最終形状Ｓ_finalまでのゴムモデルＭ１の変形を、変形１回あたりの変形量が所定量以下になるように複数段階に分け、変形実行部１２による変形処理と変形後モデル平衡化演算部１３による平衡化処理とを繰り返し実行することが好ましい。なお、図２は図示しやすいように４段階にしているが、実際はより微細な変形に留めてある。
【００３１】
［ゴムの二軸引張をシミュレーションする方法］
図１に示す装置１を用いて、ゴムの二軸引張をシミュレーションする方法について、図３を用いて説明する。
【００３２】
まず、ステップＳＴ１１において、初期設定部１０は、解析対象となるゴムモデルに関するデータの設定、分子動力学計算に必要な解析条件、二軸伸長変形条件などの各種設定を実行し、これら設定値をメモリに記憶する。
【００３３】
次のステップＳＴ１２において、初期形状ゴムモデル取得部１１は、初期形状及び構成粒子を含み、所定温度及び圧力を含む解析条件のもとで分子動力学計算により平衡化されたゴムモデルデータＭ１を取得する。本実施形態では、分子動力学計算により生成する。
【００３４】
次のステップＳＴ１３において、初期形状Ｓ_initから目標となる最終形状Ｓ_finalまでのゴムモデルＭ１の変形を、変形１回あたりの変形量が所定量以下になるように複数段階に分ける。ステップＳＴ１４〜１７を繰り返す。
【００３５】
ステップＳＴ１４において、変形実行部１２は、第１軸の変形及び第２軸の変形を少なくとも含んだ変形後の形状になるように、平衡化したゴムモデルを構成する粒子の位置を移動させ、ゴムモデルＭ１を変形させる。
【００３６】
次のステップＳＴ１５において、変形後モデル平衡化演算部１３は、変形後のゴムモデルＭ１を分子動力学計算により平衡化する。すなわち、平衡化するまで分子動力学計算を実行する。
【００３７】
次のステップＳＴ１６において、主応力算出部１４は、平衡状態にある変形後のモデルに基づき第１軸から第３軸の主応力（σ１，σ２，σ３）を算出する。
【００３８】
次のステップＳＴ１７において、微分値算出部１５は、第１軸から第３軸の主応力（σ１，σ２，σ３）および第１軸から第３軸の伸長比（λ１，λ２，λ３）、並びに式（１）〜（５）を用いて、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出する。
【００３９】
なお、本実施形態では、初期形状から最終形状までのゴムモデルＭ１の変形を複数段階に分けているが、計算が破綻しなければ、複数段階に分けなくてもよい。
【００４０】
また、本実施形態では、ゴムの変形過程における微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出するために、変形過程においてステップＳＴ１６及びＳＴ１７を実行しているが、最終形状の微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）のみが必要であれば、変形過程におけるステップＳＴ１６及びＳＴ１７の処理を省略することも可能である。
【００４１】
図４は、図２に示す単純せん断変形のシミュレーション結果を示す図である。微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）が概ね一定値になり、架橋ゴムの微分値の形状に近いことから、シミュレーションにより微分値が適切に算出されていることが分かる。
【００４２】
以上のように、本実施形態のゴムの二軸引張をシミュレーションする方法は、
初期形状Ｓ_init及び構成粒子を含み、所定温度及び圧力を含む解析条件のもとで分子動力学計算により平衡化されたゴムモデルデータを取得するステップＳＴ１２と、
第１軸の変形及び第２軸の変形を少なくとも含んだ変形後の形状になるように、平衡化したゴムモデルＭ１を構成する粒子の位置を移動させ、ゴムモデルＭ１を変形させるステップＳＴ１４と、
変形後のゴムモデルＭ１を分子動力学計算により平衡化するステップＳＴ１５と、
平衡状態にある変形後のモデルに基づき第１軸から第３軸の主応力（σ１，σ２，σ３）を算出するステップＳＴ１６と、
第１軸から第３軸の主応力（σ１，σ２，σ３）および第１軸から第３軸の伸長比（λ１，λ２，λ３）、並びに式（１）〜（５）を用いて、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出するステップＳＴ１７と、
を含む。
【００４３】
本実施形態のゴムの二軸引張をシミュレーションする装置は、
初期形状Ｓ_init及び構成粒子を含み、所定温度及び圧力を含む解析条件のもとで分子動力学計算により平衡化されたゴムモデルデータＭ１を取得する初期形状ゴムモデル取得部１１と、
第１軸の変形及び第２軸の変形を少なくとも含んだ変形後の形状になるように、前記平衡化したゴムモデルＭ１を構成する粒子の位置を移動させ、ゴムモデルＭ１を変形させる変形実行部１２と、
変形後のゴムモデルＭ１を分子動力学計算により平衡化する変形後モデル平衡化演算部１３と、
平衡状態にある変形後のモデルに基づき第１軸から第３軸の主応力（σ１，σ２，σ３）を算出する主応力算出部１４と、
第１軸から第３軸の主応力（σ１，σ２，σ３）および第１軸から第３軸の伸長比（λ１，λ２，λ３）、並びに式（１）〜（５）を用いて、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出する微分値算出部１５と、
を備える。
【００４４】
このようにすれば、従来では実現方法が開発されていなかった、ゴムの二軸引張シミュレーションを実行でき、歪みエネルギー密度関数の歪み不変量に関する微分値（∂Ｗ／∂Ｉ１，∂Ｗ／∂Ｉ２）を算出することが可能となる。また、実試験では試験できる領域が限定されているが、シミュレーションであれば、その制限を超えた範囲にて二軸引張におけるゴムの挙動を知ることができる。
【００４５】
本実施形態の方法では、初期形状Ｓ_initから目標となる最終形状Ｓ_finalまでのゴムモデルＭ１の変形を、変形１回あたりの変形量が所定量以下になるように複数段階に分け、ゴムモデルＭ１を変形させるステップＳＴ１４と変形後のモデルを平衡化するステップＳＴ１５とを繰り返し実行する。
本実施形態の装置では、初期形状Ｓ_initから目標となる最終形状Ｓ_finalまでのゴムモデルＭ１の変形を、変形１回あたりの変形量が所定量以下になるように複数段階に分け、変形実行部１２による変形処理と変形後モデル平衡化演算部１３による平衡化処理とを繰り返し実行する。
【００４６】
このようにすれば、変形１回あたりの変形量が抑制されているので、分子動力学計算が不能になり、解が求まらない事態を回避することが可能となる。
【００４７】
本実施形態に係るプログラムは、上記方法をコンピュータに実行させるプログラムである。このプログラムを実行することによっても、上記方法の奏する作用効果を得ることが可能となる。
【００４８】
以上、本開示の実施形態について図面に基づいて説明したが、具体的な構成は、これらの実施形態に限定されるものでないと考えられるべきである。本開示の範囲は、上記した実施形態の説明だけではなく特許請求の範囲によって示され、さらに特許請求の範囲と均等の意味および範囲内でのすべての変更が含まれる。
【００４９】
例えば、図１に示す各部１０〜１５は、所定プログラムをコンピュータのＣＰＵで実行することで実現しているが、各部を専用回路で構成してもよい。
【００５０】
上記の各実施形態で採用している構造を他の任意の実施形態に採用することは可能である。各部の具体的な構成は、上述した実施形態のみに限定されるものではなく、本開示の趣旨を逸脱しない範囲で種々変形が可能である。
【符号の説明】
【００５１】
１１…初期形状ゴムモデル取得部
１２…変形実行部
１３…変形後モデル平衡化演算部
１４…主応力算出部
１５…微分値算出部
Ｍ１…ゴムモデル
Ｓ_init…初期形状
Ｓ_final…最終形状
【図１】
0006568782.tif 000016
【図２】
0006568782.tif 000017
【図３】
0006568782.tif 000018
【図４】
0006568782.tif 000019