(57)【要約】【課題】クリーニングプロセスにおける外添剤群の挙動を高速に計算する方法及び装置を提供する。【解決手段】電子写真方式の画像形成装置の感光ドラムに堆積したトナー及び外添剤をクリーニングブレードによりかき取るクリーニングプロセスにおける外添剤の挙動を計算する計算方法であって、クリーニングブレードを表す従動元壁、外添剤を表し前記従動元壁に従動して運動する従動壁、及びトナーを表す粒子に対し、従動元壁と従動壁とは感光ドラムのトナー堆積面上に存在し、かつ、従動元壁と従動壁との交点と感光ドラムのトナー堆積面との距離が一定とする拘束条件を課して、モデル化した系において、コンピュータが、粒子と従動壁との相互作用を粒子と従動元壁への相互作用に変換することにより従動元壁の運動を計算し、従動元壁の運動と前記拘束条件に基づいて従動壁の運動を計算することにより、外添剤の挙動を計算する計算方法。【選択図】図１
2014071127.tif 000002
【特許請求の範囲】
【請求項１】
電子写真方式の画像形成装置の感光ドラムに堆積したトナー及び外添剤をクリーニングブレードによりかき取るクリーニングプロセスにおける外添剤の挙動をコンピュータが計算する計算方法であって、
クリーニングブレードを表す従動元壁、外添剤を表し前記従動元壁に従動して運動する従動壁、及びトナーを表す粒子に対し、従動元壁と従動壁とは感光ドラムのトナー堆積面上に存在し、かつ、従動元壁と従動壁との交点と感光ドラムのトナー堆積面との距離が一定とする拘束条件を課して、モデル化した系において、コンピュータが、粒子と従動壁との相互作用を粒子と従動元壁への相互作用に変換することにより従動元壁の運動を計算し、従動元壁の運動と前記拘束条件に基づいて従動壁の運動を計算することにより、外添剤の挙動を計算する計算方法。
【請求項２】
前記系において、外添剤とトナーの粒子との接触面の位置に外添剤を表す従動壁が設定される請求項１に記載の計算方法。
【請求項３】
コンピュータは、
粒子が従動壁に接触して及ぼす力を従動元壁に作用する力に換算し、
粒子の従動壁へのめり込み量を従動元壁へのめり込み量に換算することにより、
粒子からの力及びめり込みは従動元壁にのみ作用するものとして従動元壁の運動を計算する請求項１または２に記載の計算方法。
【請求項４】
コンピュータに、請求項１〜３のいずれか１項に記載の計算方法を実行させるプログラム。
【請求項５】
請求項４に記載のプログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項６】
電子写真方式の画像形成装置の感光ドラムに堆積したトナー及び外添剤をクリーニングブレードによりかき取るクリーニングプロセスにおける外添剤の挙動を計算する計算装置であって、
クリーニングブレードを表す従動元壁、外添剤を表し前記従動元壁に従動して運動する従動壁、及びトナーを表す粒子に対し、従動元壁と従動壁とは感光ドラムのトナー堆積面上に存在し、かつ、従動元壁と従動壁との交点と感光ドラムのトナー堆積面との距離が一定とする拘束条件を課して、モデル化した系において、粒子と従動壁との相互作用を粒子と従動元壁への相互作用に変換することにより従動元壁の運動を計算し、従動元壁の運動と前記拘束条件に基づいて従動壁の運動を計算することにより、外添剤の挙動を計算する計算装置。
【請求項７】
前記系において、外添剤とトナーの粒子との接触面の位置に外添剤を表す従動壁が設定される請求項６に記載の計算装置。
【請求項８】
粒子が従動壁に接触して及ぼす力を従動元壁に作用する力に換算し、
粒子の従動壁へのめり込み量を従動元壁へのめり込み量に換算することにより、
粒子からの力及びめり込みは従動元壁にのみ作用するものとして従動元壁の運動を計算する請求項６または７に記載の計算装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、計算方法及び計算装置に関する。
【背景技術】
【０００２】
粉体や粒体などの粒子を取り扱う分野では、個々の粒子の挙動を知ることが重要である。従来これらの分野では、実験と観察から粒子の挙動を把握することが多かった。しかし近年、計算機内でその運動を再現させ、その挙動を把握する、いわゆるシミュレーション技術の適用が活発に行われている。正確な粒子挙動のシミュレーションを行うには、粒子のもつ形状、機械的特性等の情報をできるだけ詳細に計算モデルに反映させ、それらの因子を考慮して解を求めることが重要である。
【０００３】
粒子が密集した状態での各粒子の挙動をシミュレーションする代表的な方法として個別要素法（DEM、Discrete Element Method）がある。個別要素法の技術は非特許文献１に詳細に記載されているが、個々の粒子を球形で表現し、それらの接触に際してはそのめり込み量からヘルツの式とミンドリンの式を用いることで粒子の持つ機械特性を考慮した接触抗力を求める。そしてニュートンの運動方程式（並進方向）とオイラーの運動方程式（回転方向）を用いて、その接触抗力をもとに個々の粒子の時間進展に伴う運動を求めるものである。個別要素法の特徴は個々の粒子を要素と捉え、それらが接触する部分にバネ、ダンパー、スライダを仮定して、その接触力を求め、運動方程式を解くところにある。
【０００４】
このような粒子を扱う装置の１つに、プリンタ、複写機、ファクシミリ等の電子写真技術を用いた画像形成装置がある。これらの装置では、直径数ミクロンの粉体粒子からなるトナーを用いて、帯電、露光、現像、転写、定着、クリーニングという６つのプロセスを経て、紙上に画像を形成する。図１０はクリーニングプロセスを説明する図である。図１０の上図は、クリーニング部の感光体回転軸に垂直な断面である。回転している円筒形状の感光体に板状のゴムから成るクリーニングブレードがある力で押しあてられる。図１０の下図は、感光体とクリーニングブレードの当接部付近の拡大図である。当接部に到達したトナーはクリーニングブレードに接触することで機械的にかきとられる。クリーニングブレード当接部近傍では、感光体上のトナーは十分に密集している。
【０００５】
クリーニング部にやってきたトナーはクリーニングブレードによって全てかきとられることが望ましいが、実際にはすり抜けて画像に影響を及ぼす場合があり、問題となる。そこで、クリーニング部におけるトナーやブレードの挙動をシミュレーションで解析することが求められている。非特許文献２では、DEMを用いてトナーのすり抜けの様子を解析し
ている。トナー間の接触力だけでなく、トナーとブレード間の接触力も考慮し、ブレードの運動を合わせて解析している。本文献によると、トナーがブレードを押し上げ、ブレードと感光体間に隙間が生じるためにトナーがすり抜ける。また、球形トナーと楕円形トナーを比較し、球形トナーの方がすり抜け量が多いことを解析から求めている。
【０００６】
一方で、非特許文献３によれば、クリーニングブレードと感光体が接触しているニップ近傍には外添剤が堆積している。図１１は外添剤群を示す図である。トナーについた外添剤がクリーニングブレードと感光体の間に堆積している。感光体は図中左から右へ運動しており、外添剤群はその上を滑っている。また、外添剤群はトナーやクリーニングブレードと常に接触した状態を保っている。このような外添剤群の挙動を解析することでクリーニング部における現象をより正確にシミュレーションすることができる。しかし、外添剤はサブミクロンオーダーの大きさで、その数もトナーに比べて多く、これを考慮した解析には膨大な時間を要するという問題がある。
【先行技術文献】
【非特許文献】
【０００７】
【非特許文献１】粉体工学会編、粉体シミュレーション入門(1998)
【非特許文献２】中山信行、電子写真クリーニングプロセスのシミュレーション、第18回「電磁力関連のダイナミクス」シンポジウム講演論文集、pp.165-170(2006)
【非特許文献３】新井和彦、中山信行、小笠原正、織田康弘、勅使川原亨、ゴムブレードクリーニングにおける感光体摩耗ばらつきの抑制、Imaging Conference JAPAN 2010 論文集、pp.247-250(2010)
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明の目的は、クリーニングプロセスにおける外添剤群の挙動を高速に計算することにある。
【課題を解決するための手段】
【０００９】
本発明は、電子写真方式の画像形成装置の感光ドラムに堆積したトナー及び外添剤をクリーニングブレードによりかき取るクリーニングプロセスにおける外添剤の挙動をコンピュータが計算する計算方法であって、
クリーニングブレードを表す従動元壁、外添剤を表し前記従動元壁に従動して運動する従動壁、及びトナーを表す粒子に対し、従動元壁と従動壁とは感光ドラムのトナー堆積面上に存在し、かつ、従動元壁と従動壁との交点と感光ドラムのトナー堆積面との距離が一定とする拘束条件を課して、モデル化した系において、コンピュータが、粒子と従動壁との相互作用を粒子と従動元壁への相互作用に変換することにより従動元壁の運動を計算し、従動元壁の運動と前記拘束条件に基づいて従動壁の運動を計算することにより、外添剤の挙動を計算する計算方法である。
【００１０】
本発明は、電子写真方式の画像形成装置の感光ドラムに堆積したトナー及び外添剤をクリーニングブレードによりかき取るクリーニングプロセスにおける外添剤の挙動を計算する計算装置であって、
クリーニングブレードを表す従動元壁、外添剤を表し前記従動元壁に従動して運動する従動壁、及びトナーを表す粒子に対し、従動元壁と従動壁とは感光ドラムのトナー堆積面上に存在し、かつ、従動元壁と従動壁との交点と感光ドラムのトナー堆積面との距離が一定とする拘束条件を課して、モデル化した系において、粒子と従動壁との相互作用を粒子と従動元壁への相互作用に変換することにより従動元壁の運動を計算し、従動元壁の運動と前記拘束条件に基づいて従動壁の運動を計算することにより、外添剤の挙動を計算する計算装置である。
【発明の効果】
【００１１】
本発明によれば、クリーニングプロセスにおける外添剤群の挙動を高速に計算することができる。
【図面の簡単な説明】
【００１２】
【図１】外添剤群を表す壁を説明する図。
【図２】本装置の構成を示す図。
【図３】本装置の概略構成を示す図。
【図４】バネ支持壁モデルを示す図。
【図５】粒子バネと壁可動方向の関係を示す図。
【図６】従動壁に作用する粒子バネの効果の変換方法を示す図。
【図７】壁運動計算部の処理フローを示す図。
【図８】実施例の計算条件を示す図。
【図９】トナー形状によるブレード変位量の違いを示す図。
【図１０】クリーニングプロセスを示す図。
【図１１】外添剤群を示す図。
【発明を実施するための形態】
【００１３】
以下、本発明の実施例について図に基づき説明する。本実施例は、電子写真方式の画像形成装置の感光ドラムに堆積したトナー及び外添剤をクリーニングブレードによりかき取るクリーニングプロセスにおける外添剤の挙動を計算する計算方法に関する。
【００１４】
図１（a）は外添剤群の運動の様子を示したものである。左図のように外添剤群はトナ
ーやクリーニングブレードと接触した状態で、図の左から右へ動く感光体の上を滑っている。右図はクリーニングブレードがトナーによる力を受けて押し上げられた場合を示している。外添剤群もトナーによって押されるため、クリーニングブレードと接触した状態を保って運動する。一方で、外添剤群は引き続き感光体上を滑っている。このとき、外添剤群とトナーの接触面がクリーニングブレードと交わる点の運動に着目する。交点と感光体との距離が常に一定となるように運動することになる。
【００１５】
図１（b）は外添剤群の運動をモデル化したものである。外添剤群はトナーとの接触面
のみを考えれば十分である。よって、接触面のみとし、壁で表す。クリーニングブレードも壁で表す。クリーニングブレードを表す壁がトナー（粒子）による力を受けて運動した時、外添剤群を表す壁もその運動に従動して運動させる。以下では、外添剤群を表す壁を従動壁、クリーニングブレードを表す壁を従動元壁とする。従動壁と従動元壁は感光体のトナー堆積面上に存在し、その交点はトナー堆積面との距離を一定に保つという拘束条件を満たして運動する。また、粒子は従動壁に接触して力を及ぼす。この接触力は従動元壁に作用するものとする。さらに、接触力の基となる粒子の従動壁へのめり込み量は従動元壁へのめり込み量とする。結果、粒子からの力は従動元壁にのみ作用するとして運動を計算することになる。粒子からの力を受けて従動元壁が運動し、従動壁は従動元壁の運動に従動して運動する。このようにすることで、上記拘束条件を満たした従動壁の運動を容易に計算することが可能となる。また、外添剤群は単純な壁で表現するので、外添剤群に含まれる個々の粒子の運動を考えることなく解析が可能となり、計算の高速化が実現できる。
【００１６】
図２は、本実施例における解析装置が動作する計算装置を示すブロック図である。この計算装置の本体１０は、各種判断及び処理を行う中央処理装置（ＣＰＵ）１１と、プログラム及び固定データを格納するＲＯＭ１２と、データを記憶するＲＡＭ１３と、外部記憶装置との間でデータをやり取りする入出力回路（Ｉ／Ｏ）１４とから構成されている。この計算装置本体１０には、入力データ１５がＩ／Ｏ１４を介して入力される。またその計算結果を、Ｉ／Ｏ１４を介して出力データ１６として出力する。
【００１７】
図３は、粒子挙動解析を行うプログラム部分の構成を示すブロック図である。
【００１８】
B100は全体制御部であり、以下に説明する入力データ読込部、粒子・壁受力計算部、粒子運動計算部、壁運動計算部、計算結果出力部を制御する。
【００１９】
B110は入力データ読込部であり、粒子形状データ、壁形状データ、及び各種パラメータを読み込んでＲＡＭ１３に格納する。粒子形状データとは、粒子の表面を球形で表現したものである。壁形状データは、壁の表面を複数の多角形で表現したものである。粒子形状データは粒子の中心座標、半径と、各粒子の番号として入力されるものである。壁形状データは、多角形の頂点座標と、各多角形を構成する頂点の番号として入力されるものであ
る。なお、壁は粒子等から力を受けて運動する従動元壁と、従動元壁の運動に従動して運動する従動壁からなる。
【００２０】
各種パラメータには、粒子と壁の機械特性（ヤング率、ポアソン比、減衰係数、摩擦係数、比重）、粒子に働く力（重力、付着力、電磁力、空気の粘性特性等）、及び計算条件としての計算刻み時間（Δt）、計算終了時刻がある。
【００２１】
B120は、粒子・壁受力計算部であり、各粒子と壁が受ける力をそれぞれ計算するものである。粒子と壁は他の粒子または壁と接触している場合に、そのめり込み量に応じて接触抗力を受ける。
【００２２】
B130は粒子運動計算部であり、B120で求めた粒子に働く力による、各粒子の刻み時間後の速度と位置を求めるものである。
【００２３】
B140は壁運動計算部であり、B120で求めた粒子に働く力による、従動元壁の刻み時間後の速度と位置を求めるものである。また、従動元壁の運動に従って従動壁の刻み時間後の速度と位置を求めるものである。
【００２４】
B150は計算結果出力部であり、得られた各粒子及び壁の位置と速度の結果を出力データ１６として出力するものである。
【００２５】
以下に、B120〜140の詳しい内容と計算方法を示す。
【００２６】
まず、B120の粒子・壁受力計算部について説明する。粒子同士、または粒子と壁が接触状態にあるとき、接触面とめり込み量が定義される。めり込み量に対して、ヘルツの式を用いて法線方向の接触力（バネ力）を求める。また各接触点に対して、接触が始まってからの接触面（上記接触面法線方向に垂直な面）内方向の移動量を求め、ミンドリンの式を用いて粒子iに働く接触面内方向の接触力（バネ力）を求める。また接触面法線方向の運
動速度（めり込み速度）と面内方向の各運動速度（ずり速度）から、各方向に作用する粘性力（ダンパー力）を求める。そしてそれら接触力と粘性力の和に、粒子に働くその他の力（重力、付着力、電磁力、空気の粘性抵抗等）を加えて、各粒子が受ける並進力と回転力を求める。なお、本めり込み量から接触力を求める過程は、上記非特許文献１に記載されている従来技術と同じである。
【００２７】
次に、B130の粒子運動計算部における粒子の運動について説明する。B120で求めた粒子に働く力による、各粒子の刻み時間後の速度と位置を求めるものである。すなわち、並進力からニュートンの運動方程式を用いて並進方向の加速度を求め、更に速度と位置を求める。また回転力からオイラーの運動方程式を用いて角加速度を求め、更に回転速度と回転位置を求める。本処理は、非特許文献１で行われている従来技術と同じである。
【００２８】
次に、B140の壁運動計算部における従動元壁の運動について説明する。図４は壁板の運動を表した模式図である。壁板はバネ定数Kwalのバネによって支持される。一方、壁板は接触している粒子からもバネによる力を受ける。粒子iのバネ定数をkgrn(i)、壁板と接触している粒子n個のバネ定数の総和をKgrnとする。壁板は壁板支持バネと粒子バネからの
力を受けて運動し、釣り合いの位置で停止する。このとき、壁板と粒子が離れる瞬間があると、壁板は即座に初期位置に戻ることになり、粒子挙動計算が破綻する。したがって、壁板は粒子との接触を保ったまま各時間ステップで微小並進運動させる。非特許文献２によれば、時刻tからt+Δtの間の壁板の変位ベクトルΔUは数１で表される。
【数１】
2014071127.tif 000003
ここで、Fwalは壁板の支持バネからの受力ベクトル、Fgrnは壁板の粒子からの受力ベクトル、Fttlは壁板の総受力ベクトルである。また、上付き文字は時刻を表す。壁板の位置は支持バネと粒子バネの釣り合いの位置から決定されることから、
【数２】
2014071127.tif 000004
とした分だけ壁板を移動させれば、粒子との接触を保ったまま移動させることができる。しかし、数１で算出される変位量が大きすぎて解が振動することがある。この場合は、
【数３】
2014071127.tif 000005
として変位量を緩和した数４を用いるのが良い。
【数４】
2014071127.tif 000006
ここで、a（0≦a＜1）は壁運動の緩和係数に当たる。
【００２９】
数４の各値の算出方法について説明する。
壁板支持バネのバネ定数Kwalにはクリーニングブレードの機械物性に基づいた値を与える。
【００３０】
Fwalは数５で表される。
【数５】
2014071127.tif 000007
ここで、Lは壁板の初期配置（粒子と接触していない状態）からの変位ベクトルである。
【００３１】
Fgrnは壁板と接触している粒子iからの受力ベクトルfgrn(i)の総和である。壁板と接触している粒子の数をnとして数６で表される。
【数６】
2014071127.tif 000008
fgrn(i)は、粒子・壁受力計算部B120で算出される値である。
【００３２】
Kgrnは、壁板に接触する各粒子バネのバネ定数の総和である。各粒子バネは線形バネであるとして数７で表される。
【数７】
2014071127.tif 000009
【００３３】
数７について、図５を用いて説明する。fgrn(i)は、図５の左図のように粒子バネ力ベ
クトルfgrn(i)の壁可動方向成分である。Vmovは壁可動方向の単位ベクトルである。δ(i)は図５の右図のように粒子iの壁平面法線方向のめり込み量δ０(i)を壁可動方向のめり込み量に換算したものである。Vnorは壁平面法線方向の単位ベクトルである。
【００３４】
なお、ここで説明した従動元壁の変位量算出方法は非特許文献２に記載されている従来技術と同じである。
【００３５】
次に、従動壁の運動について説明する。従動壁は図１（a）のように従動元壁と交点を
持つように堆積面上に配置する。図は２次元であるが、従動壁と従動元壁は図の奥行き方向に一様に存在する。従動元壁が運動した時、交点と堆積面との距離が一定となるように従動壁が運動する。このようにすることで、従動壁は従動元壁との接触条件を保ったまま、移動する堆積面上を自在に並進することになる。
【００３６】
従動壁も粒子と接触するため、力を受ける。従動壁が粒子から受ける力は全て従動元壁が受ける力に変換する。変換について図６を用いて説明する。従動壁に接触している粒子iからの受力ベクトルfgrn(i)はそのまま従動元壁の受力ベクトルとする。数７で必要となる可動方向のめり込み量δ(i)は、図より数８で表される。
【数８】
2014071127.tif 000010
ここで、Umovは従動壁可動方向単位ベクトルであり堆積面に平行な向きを持つ。Unorは従動壁面法線方向単位ベクトル、δ’０(i)は粒子iの従動壁面法線方向のめり込み量、δ’(i)は粒子iの従動壁可動方向のめり込み量である。
【００３７】
以上により、従動壁に作用する粒子バネの効果が従動元壁に引き継がれる。直接、従動
元壁と接触している粒子バネの効果と合わせて、数４により従動元壁が運動する。そして、交点と堆積面との距離が一定という拘束条件のもとで従動壁が運動する。このとき、従動壁は堆積面に平行な方向に運動することになる。なお、ここで説明したのは従動壁に作用する粒子バネの力とめり込み量を、従動元壁に作用する力とめり込み量に変換する方法の一例であり、他の変換方法を用いてもよい。重要なのは、従動壁と粒子の相互作用を計算することなく従動壁の移動量を算出することである。
【００３８】
ここで、B140の壁運動計算部の処理フローを図７で説明する。
【００３９】
まず、従動元壁と接触している粒子iを設定する（S100）。次に従動元壁の可動方向単
位ベクトルVmovを設定する（S110）。数７を基に粒子iのバネ定数kgrn(i)を算出する（S120）。従動元壁と接触している粒子全てについてS100からS120の操作を繰り返す（S130）。次に従動壁と接触している粒子iを設定する（S140）。粒子iの従動壁への受力fgrn(i)
を従動元壁受力に変換する（S150）。また、数８を基に粒子iの従動壁へのめり込み量δ
’０(i)を従動元壁へのめり込み量に変換する（S160）。従動壁と接触している粒子全て
についてS140からS160の操作を繰り返す（S170）。ここまでで算出された各粒子の受力ベクトルfgrn(i)とバネ定数kgrn(i)を使って数４を基に従動元壁の変位ベクトルΔUを算出
する（S180）。従動元壁と従動壁の交点と粒子堆積面との距離が一定という拘束条件のもとで従動壁を運動させる（S190）。
【００４０】
なお、粒子の形状は球形としたが、本手法はめり込み量を指定すれば粒子の形状に依らず適用できるものである。
【００４１】
（実施例）
以上の壁運動計算方法を用いて粒子挙動計算を実施した例を、図８〜図９を用いて説明する。図８（a）は解析モデルを示す。解析モデルは電子写真技術を用いた画像形成装置
におけるクリーニングプロセスである。0.1m/secで矢印の方向に移動する感光体ドラムに対してクリーニングブレードをβ＝30°で当接させておき、感光体ドラム上に乗って運ばれてきたトナーの運動を見る。各トナーの形状を粒子形状データとして与え、感光体ドラム、クリーニングブレード及び外添剤群の形状を壁形状データとして与える。クリーニングブレードを支持するバネのバネ定数Kwalは20,000N/m、壁運動の緩和係数aは0.1である
。従動元壁であるクリーニングブレードと従動壁である外添剤群の交点は、感光体からの距離14μmに拘束される。クリーニングブレードの可動方向ベクトルVmovは各粒子iからの受力ベクトルfgrn(i)と同じ向きとしている。
【００４２】
図８（b）は、解析条件として与えた、トナー、感光体ドラム、クリーニングブレード
それぞれの機械特性（ヤング率、ポアソン比、反発係数、滑り摩擦係数、付着力）を示している。この条件で、刻み時間Δtを2.5nsecに設定し、約４千万ステップの計算を実行した。
【００４３】
図９は、（a）異形トナーと（b）球形トナーの結果を比較したものである。どちらも外添剤群がクリーニングブレードの変形に従動して、感光体上を自由に滑りながら運動している様子を計算できている。クリーニングブレードの変位量を比較すると、球形トナーの方が異形トナーより変位量が大きい。球形トナーの方が回転しやすいことからクリーニングブレードを押し上げる力が強く、トナーがすり抜けやすい。図９の計算結果では、この様子が再現できており、本手法で外添剤群の効果を精度良くシミュレーションできていることが分かる。
【００４４】
上記の実施例は本発明を説明するための一例である。上記の計算方法を実行するコンピュータ（計算装置）、コンピュータに上記の計算方法を実行させるプログラム、当該プロ
グラムを記録したコンピュータ読み取り可能な記録媒体は本発明に含まれる。
【符号の説明】
【００４５】
１０：コンピュータ、１１：ＣＰＵ
【図１】
2014071127.tif 000011
【図２】
2014071127.tif 000012
【図３】
2014071127.tif 000013
【図４】
2014071127.tif 000014
【図５】
2014071127.tif 000015
【図６】
2014071127.tif 000016
【図７】
2014071127.tif 000017
【図８】
2014071127.tif 000018
【図９】
2014071127.tif 000019
【図１０】
2014071127.tif 000020
【図１１】
2014071127.tif 000021