(57)【特許請求の範囲】
【請求項１】
長尺のレール体の上をこれに沿った進行方向へ向けて転動する車輪体の接触による応力解析のための有限要素解析システムであって、
前記レール体をその長手方向に沿って単位ブロックの複数に仮想分割し、少なくとも前記車輪体及び前記レール体の接触部を含む前記単位ブロックを挟んで前後所定数ずつの前記単位ブロックからなる解析区間を定めた上で、
前記解析区間内を前記接触部が前記進行方向に所定進行距離だけ移動するまで有限要素解析を行い、前記解析区間を前記進行方向に前記単位ブロックの１つ分だけずらして新たな解析区間を定め、更に前記有限要素解析を行うことを繰り返す解析手段を含むことを特徴とする有限要素解析システム。
【請求項２】
前記単位ブロック毎に独立したメッシュモデルを付与することを特徴とする請求項１記載の有限要素解析システム。
【請求項３】
前記有限要素解析を行うに際し、新たな解析区間を定めたときに、新たに解析に組み込まれた前記単位ブロックに対してメッシュモデルを追加付与していくことを特徴とする請求項２記載の有限要素解析システム。
【請求項４】
仮想分割に続いて、前記レール体の全体に亘ってメッシュモデルを付与することを特徴とする請求項２記載の有限要素解析システム。
【請求項５】
前記単位ブロック毎の前記メッシュモデルは同一であることを特徴とする請求項２乃至４のうちの１つに記載の有限要素解析システム。
【請求項６】
前記解析区間の前記メッシュモデルにはその要素節点を特定するための特定符号を与え、新たな解析区間を定めたときに、新たに解析に組み込まれた前記単位ブロックの対応する要素節点に解析から外れた前記単位ブロックの前記特定符号を付与することを特徴とする請求項５記載の有限要素解析システム。
【請求項７】
前記所定進行距離は前記レール体の長手方向への前記単位ブロックの長さであることを特徴とする請求項１乃至６のうちの１つに記載の有限要素解析システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、有限要素法を用いた有限要素解析システムに関し、特に、長尺のレール体の上をこれに沿った進行方向へ向けて転動する車輪体による接触応力解析のための有限要素解析システムに関する。
【背景技術】
【０００２】
有限要素法では、一般的に、解析モデルの要素を細かく設定するほど解析精度が向上し得る。その一方で要素の数が増大し、計算処理に時間が必要となってしまう。そこで大規模な解析モデルとなる長尺物の解析、例えば、路盤を転動する車輪体の移動方向に沿った路盤への応力などの解析の如きでは、許容される解析精度と解析時間との関係で各種の工夫が必要とされる。
【０００３】
例えば、特許文献１では、ローラに当接して移動せしめられる長尺のコンベヤベルトについて、ローラとの接触によりコンベヤベルトに生じる応力及び／又は歪みについて有限要素法を用いて解析する方法が開示されている。ここでは、コンベヤベルトの要素をローラの当接位置において細かく、当接位置から離間するほど粗く設定して、解析精度の向上と処理時間の短縮の両立を図っている。なお、コンベヤベルトは中心間隔Ｌで並ぶローラによって支持されるため、有限要素法による計算処理は長さＬのコンベヤベルトに対してのみ行うことで十分である。
【０００４】
更に、特許文献２では、玉軸受けやころ軸受けといった転がり軸受けの転動疲労による剥離等の損傷状態の推定を目的とした有限要素法による解析方法が開示されている。ここで、剥離等は転がり軸受けの金属製基体に内在する非金属介在物に起因するとされる。また、転動体が金属製基体の表面上を一回だけ通過するのではなく、繰り返し転がる場合には、歪（内部損傷）が金属製基体の内部に蓄積される。そこで、金属製基体の表面近傍に内在する歪の蓄積の原因となる非金属介在物の周囲でマトリクスの応力計算を有限要素法で行っている。かかる有限要素法において、金属製基体に与えられる要素は、非金属介在物の近傍及び転動体の移動範囲内での接触位置において細かく、これらの位置から離れるにつれ粗く設定している。非金属介在物に基づく転動疲労のシミュレーションであるから、転動体の移動範囲を非金属介在物の周囲に限定して計算処理を行ったとしても解析精度を損なうことがなく、また解析時間も短縮出来るのである。
【０００５】
一方で、例えば、特許文献３では、転動体であるタイヤが水膜を与えられた路面上を移動する挙動をシミュレーションする方法が開示されている。タイヤが高速で転動した場合、タイヤとともに、路面を長距離に亘ってモデル化してシミュレートする必要がある。そこで、流体解析領域をタイヤの進行方向に対して所定の長さの範囲に抑えた上で、タイヤの移動に追従するように流体解析領域を移動させる解析方法を採用している。詳細には、タイヤモデルをその移動速度と解析時間の刻み幅に応じた距離だけ移動する毎に流体解析領域を同じ距離だけ移動させ、流体解析領域中に設定された各要素に対して境界条件を再設定するのである。かかる方法によれば、タイヤモデルの移動方向に長い流体解析領域を設定する必要がなく、解析時間を短縮できる。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００５−７５６１７号公報
【特許文献２】特開２０１１−６４６４４号公報
【特許文献３】特開２００４−３３８６６０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
ところで、長尺のレール体の上をこれに沿った進行方向へ向けて転動する車輪体による接触応力解析においても、許容される解析精度と解析時間との関係で各種の工夫が必要とされる。特に、鉄道車両の高速化が進む中で、より長尺のレール体における解析を必要とされ、かかる要求は顕著となっている。
【０００８】
本発明は、かかる事情に鑑みてなされたものであって、その目的とするところは、長尺のレール体の上をこれに沿った進行方向へ向けて転動する車輪体による接触応力解析のための有限要素解析システムにおいて、高い解析精度を保ちつつ解析時間を短縮させたシステムを提供することにある。
【課題を解決するための手段】
【０００９】
本発明による有限要素解析システムは、長尺のレール体の上をこれに沿った進行方向へ向けて転動する車輪体の接触による応力解析のための有限要素解析システムであって、前記レール体をその長手方向に沿って単位ブロックの複数に仮想分割し、少なくとも前記車輪体及び前記レール体の接触部を含む前記単位ブロックを挟んで前後所定数ずつの前記単位ブロックからなる解析区間を定めた上で、前記解析区間内を前記接触部が前記進行方向に所定進行距離だけ移動するまで有限要素解析を行い、前記解析区間を前記進行方向に前記単位ブロックの１つ分だけずらして新たな解析区間を定め、更に前記有限要素解析を行うことを繰り返す解析手段を含むことを特徴とする。
【００１０】
かかる発明によれば、長尺のレール体に亘って与えられた要素間の所定の計算を繰り返し行う有限要素解析システムであっても、解析精度を維持したまま計算回数を減じ得る。すなわち、解析時間を短縮できるのである。また、解析区間を複数の単位ブロックからなるようにすることで、解析区間を新たに定めるステップを簡素化出来て、解析時間を短縮できる。
【００１１】
上記した発明において、前記単位ブロック毎に独立したメッシュモデルを付与することを特徴としてもよい。かかる発明によれば、解析区間を新たに定めるステップの処理を簡素化でき、解析時間を短縮できる。
【００１２】
上記した発明において、前記有限要素解析を行うに際し、新たな解析区間を定めたときに、新たに解析に組み込まれた前記単位ブロックに対してメッシュモデルを追加付与していくことを特徴としてもよい。かかる発明によれば、解析区間を新たに定めるステップにおいて、メッシュモデルの付与の処理を簡素化し、解析時間を短縮できる。
【００１３】
上記した発明において、仮想分割に続いて、前記レール体の全体に亘ってメッシュモデルを付与することを特徴としてもよい。かかる発明によれば、解析区間を新たに定めるステップにおいて、メッシュモデルの付与の処理を簡素化出来て、解析時間を短縮できる。
【００１４】
上記した発明において、前記単位ブロック毎の前記メッシュモデルは同一であることを特徴としてもよい。かかる発明によれば、解析区間を新たに定めるステップにおいて、メッシュモデルの付与の処理を簡素化出来て、解析時間を短縮できる。
【００１５】
上記した発明において、前記解析区間の前記メッシュモデルにはその要素節点を特定するための特定符号を与え、新たな解析区間を定めたときに、新たに解析に組み込まれた前記単位ブロックの対応する要素節点に解析から外れた前記単位ブロックの前記特定符号を付与することを特徴としてもよい。かかる発明によれば、単位ブロックに与えられる特定符号の数を減じて、解析区間を新たに定めるステップにおいて、メッシュモデルへの特定符号の付与の処理を簡素化出来て、解析時間を短縮できる。
【００１６】
上記した発明において、前記所定進行距離は前記レール体の長手方向への前記単位ブロックの長さであることを特徴としてもよい。かかる発明によれば、解析区間を新たに定めるステップを簡素化出来て、解析時間を短縮できる。
【図面の簡単な説明】
【００１７】
【図１】本発明による有限要素解析システムを示すブロック図である。
【図２】本発明による有限要素解析システムにおける有限要素解析の手順を示すフロー図である。
【図３】有限要素解析システムに使用される解析モデルの図である。
【図４】有限要素解析システムに使用される解析モデルの図である。
【図５】有限要素解析システムによる解析モデルの荷重条件の図である。
【図６】有限要素解析システムによる解析結果のうちの車輪の並進速度履歴を示す図である。
【図７】有限要素解析システムによる解析結果を示す図である。
【発明を実施するための形態】
【００１８】
本発明の１つの実施例である有限要素解析システムについて、図１を用いてその詳細を説明する。
【００１９】
図１に示すように、有限要素解析システム３０は、レール上を転動しながら進行する車輪をシミュレートし、レールと車輪との接触に伴う応力解析（若しくは、歪み解析）を行うコンピュータによるシステムである。有限要素解析システム３０は、制御部３１及び記憶部３２を含み、有限要素解析を行うための解析モデルの作成を支援するための解析モデル作成部３３、後述する解析区間を定めるための解析区間設定部３４、及び、作成された解析モデルの解析区間における解析を行うための解析部３５を含む。
【００２０】
制御部３１は所定のプログラムを実行する図示しないＣＰＵとそれに付随するＲＯＭ及びＲＡＭなどを含む。制御部３１において解析モデル作成部３３、解析区間設定部３４及び解析部３５の各プログラムを実行することで計算処理を行う。すなわち、制御部３１は所定のプログラムを実行することで有限要素解析システム３０の解析手段として機能し得る。なお、解析モデル作成部３３、解析区間設定部３４、及び解析部３５の詳細については後述する。
【００２１】
記憶部３２は、解析モデル作成部３３によって作成された解析モデルを記憶し、解析部３５によって算出された解析結果を記憶する。つまり、制御部３１の動作に必要となる入力データと、制御部３１の算出した出力データとを必要に応じて記憶することができる。
【００２２】
次に、有限要素解析システム３０における解析処理の手順について、図２及び図３を用いて説明する。
【００２３】
図２に示すように、有限要素解析システム３０では、解析モデル作成部３３を起動させ、解析モデルの作成を行う（Ｓ１）。図３に示すように、解析モデル４０の対象としてレール４２の上を転動する車輪４１が選択されると、まず、レール４２の所定の区間を長手方向に沿って複数の単位ブロックに仮想分割する。ここでは、単位ブロックの形状を全て同一とする。この解析モデル４０の対象は記憶部３２に記憶される。なお、後述する解析区間を定めるための指定値など、解析に必要な所定値はこの時点で入力されて、併せて記憶部３２に記憶させてもよい。
【００２４】
次に、解析区間設定部３４は、レール４２についての解析区間の設定を行う（Ｓ２）。解析区間は車輪４１のレール４２に接触する接触部Ｐの位置に基づいて定められる。すなわち、図３（ａ）に示すように、接触部Ｐが単位ブロック４２ｂ上にあるときに、単位ブロック４２ｂと、これを挟んで前後所定数ずつの単位ブロックとからなる単位ブロックを解析区間として定める。このとき、所定数は、解析精度を高く維持できるよう、解析結果に与える影響の比較的大きな単位ブロックを解析区間に少なくとも含むように定めることが好ましい。ここでは、車輪４１の進行方向（Ｘ方向）の前方に１つ、後方に１つの、計３つの連続した単位ブロック４２ａ、４２ｂ及び４２ｃをレール４２の解析区間として定めた。
【００２５】
続いて、解析モデル作成部３３は、車輪４１及びレール４２にメッシュモデルを付与する（Ｓ３）。なお、車輪４１には、あらかじめ入力された条件に従ってメッシュモデルが与えられるが、図面ではこのメッシュモデルを省略した。レール４２には、解析区間として定められた単位ブロック４２ａ〜４２ｃと、ダミーブロックである単位ブロック４２ｄとの４つの単位ブロックに対して、単位ブロック毎に独立するようなメッシュモデルを与える。つまり、メッシュモデルを構成する１つの要素が複数の単位ブロックに亘って跨ることなく、単位ブロック内に納められるようにメッシュモデルが与えられる。また、単位ブロック４２ａ〜４２ｄに付与されるメッシュモデルは同一である。
【００２６】
続いて、解析部３５は、接触部Ｐを所定距離だけ移動させるまで有限要素法による計算処理を行う（Ｓ４）。詳細には、接触部Ｐの位置が、単位ブロック４２ｂの後方端にあるときから、車輪４１の転動により単位ブロック４２ｂの前方端に到達するまで計算処理を行うのである。
【００２７】
更に、図３（ｂ）に示すように、接触部Ｐが単位ブロック４２ｂの前方端に到達するまでの計算処理をしたら、終了条件を満たしたかを判定する終了判定を行う（Ｓ５）。終了判定が否定された場合には（Ｓ５；Ｎ）、解析区間の設定（Ｓ２）に戻る。なお、終了条件については、例えば、接触部Ｐが解析の終了地点として定めた位置まで移動した場合や、解析により算出される値が所定値へ到達した場合など、種々の場合の設定が可能である。
【００２８】
終了判定が否定された場合には、続いて、解析区間設定部３４は、解析区間の設定（Ｓ２）を再び行う。このとき、接触部Ｐは単位ブロック４２ｃの後方端に到達し、単位ブロック４２ｃの上にある。すなわち、レール４２の解析区間は、単位ブロック４２ｃと、これを挟んで前後に１つずつの単位ブロック４２ｂ、４２ｃ及び４２ｄに定められる。つまり、解析区間の後方端と前方端をそれぞれ進行方向（Ｘ方向）に単位ブロック１つ分だけずらして新たに設定する。このとき、解析区間から外れた単位ブロック４２ａについての解析結果は記憶部３２に保存される。
【００２９】
続いて、解析モデル作成部３３はメッシュモデルの付与（Ｓ３）を再び行う。新たに定められた解析区間には単位ブロック４２ｂ〜４２ｃが含まれるが、これには先の計算処理において既にメッシュモデルが付与されている。つまり、ダミーブロックとして新たに解析に組み込まれる単位ブロック４２ｅにだけ、単位ブロック４２ａと同一のメッシュモデルが追加付与されるのである。なお、車輪４１のメッシュモデルについても既に付与されたものと同一である。
【００３０】
続いて、上記と同様に、解析部３５により新たに定められた解析区間において接触部Ｐが単位ブロック４２ｃの前方端に到達するまで計算処理を行い、以後、同様に計算処理を続ける。
【００３１】
一方、終了判定（Ｓ５）において、終了条件を満たした場合、終了判定は肯定され（Ｓ５；Ｙ）、有限要素解析システム３０における解析処理が終了する。
【００３２】
このように、レール４２の全長に亘る要素全てについて常に解析を行うのではなく、接触部Ｐの位置を基準に解析区間をレール４２の一部の範囲に限定し、この解析区間を接触部Ｐの位置に応じて移動させて、解析の計算規模を小さくできる。つまり、メッシュモデルの細かさを維持して高い解析精度を保ちつつ、解析時間を短縮できる。
【００３３】
なお、解析区間の設定（Ｓ２）において解析区間を移動させる処理は、接触部Ｐが所定距離だけ移動する毎に行えばよく、接触部Ｐが単位ブロックの端部に到達したとき以外に行うこともできる。このとき、所定距離を単位ブロックの長さと同一にすると、解析区間と車輪４１との位置関係を常に一定の範囲内とできて、解析精度を安定させ得る。
【００３４】
ここで、メッシュモデルにより定められた要素の節点のうち、解析区間に加えられている単位ブロックの各節点にはこれらを特定するための節点番号が付与される。更に、単位ブロックの前方端と後方端に位置する節点においては、隣り合う単位ブロック同士で節点番号を共有している。すなわち、図３（ａ）を再度参照すると、解析区間の単位ブロック４２ａ〜４２ｃの各節点にはそれぞれ節点番号１〜２０が付与されているが、単位ブロック４２ａと単位ブロック４２ｂとで節点番号７及び８を共有しており、単位ブロック４２ｂと単位ブロック４２ｃとで節点番号１３及び１４を共有している。
【００３５】
また、解析区間の移動に際して節点番号を新たに付してこれを増やさないため、解析区間の単位ブロックにおいて同一の節点番号が重複しないよう、解析からはずれた節点番号を新たに解析に組み込まれた単位ブロックの節点番号として使用する。なお、解析に組み込まれる単位ブロックには、後述するダミーブロックも含む。また、節点番号の各節点への付与を効率よく行うため、節点番号の単位ブロック内での組み合わせを定めておくとよい。
【００３６】
詳細には、図３（ａ）において、解析区間にない単位ブロック４２ｄは、ダミーブロックとして、解析区間にある単位ブロック４２ｃと節点番号１９及び２０を共有し、これに加え節点番号１９から２４までが付与されている。つまり、単位ブロック４２ｄは、解析区間にはないが、その一部の節点番号１９及び２０は単位ブロック４２ｃと共有するダミーブロックとして解析に組み込まれている。また、図３（ｂ）では、単位ブロック４２ｄが解析区間に新たに加えられ、これと同時にダミーブロックとして新たに解析に組み込まれた単位ブロック４２ｅに単位ブロック４２ａに対応する同一のメッシュモデルが追加付与される。つまり、節点番号１から６までを設定するが、節点番号１及び２は単位ブロック４２ｄと共有するのである。これにより、節点番号の数を抑制できる。続いて、図３（ｃ）では、解析区間に新たに単位ブロック４２ｅを加え、同様の処理が行われる。
【００３７】
次に、上記した有限要素解析システム３０により解析を行って得た計算処理結果を、レールの全長の要素について常に計算を行う従来式での解析による計算処理結果と比較した例について、図４乃至図７を用いて説明する。なお、有限要素解析システム３０による解析を「移動式」と表示する。
【００３８】
図４に示すように、移動式として、メッシュモデルを付与した車輪５１と、レール５２の全長に亘って単位ブロックに分割した上でメッシュモデルを所定の単位ブロックに付与したレール５２とを含む解析モデル５０において、レール５２は同一形状の複数の単位ブロック５２ａ〜５２ｆに分割されている。他の単位ブロックについてはその表示を省略した。この移動式での解析において、車輪５１の接触する単位ブロック５２ｂと、この後方の単位ブロック５２ａと、その前方の単位ブロック５２ｃ〜５２ｅとの合計５つの単位ブロックが解析区間である。なお、最前方に表示されている単位ブロック５２ｆは、解析区間に加えられている単位ブロック５２ｅと一部の節点番号を共有するダミーブロックとして解析に組み込むべく設定されている。また、各単位ブロックの進行方向の長さは任意に設定できるが、ここでは０．１２５ｍとした。なお、車輪５１には、側面に突出する要素５４が設定されており、解析に影響のない程度に剛性及び密度を小さく設定することで、車輪５１の転動を明示するためのマーカーとして用いた。
【００３９】
なお、移動式での解析では、解析に組み込まれることのない単位ブロックを設定する必要はない。つまり、解析の開始から終了までに解析区間に加えられる単位ブロックと最後にダミーブロックとして設定される単位ブロックとを合わせた長さがレールの全長となる。
【００４０】
図５に示すように、解析モデル５０にはいくつかの解析条件を与える。ここでは、車軸５３にトルクを与えるため、その外周において転動する向き（周方向）に０．７５ｋＮの集中荷重を４カ所均等に与え、その中心において下向きに車体重量に相当する５０．０ｋＮの荷重を与えた。また、レール５２の底面を完全固定する拘束条件を与えるとともに、車軸５３がＹ軸方向に移動しないよう、車軸５３の端面５３ａをＹ軸方向に固定する拘束条件を与えた。
【００４１】
また、従来式での解析では、レールの全長を５ｍと１０ｍの２通りとし、上記した解析モデル５０と同じメッシュモデルをレールの全長に亘って与えた。更に、従来式での解析では、解析区間を定めずに常にレールの全長について解析を行った。それ以外の解析条件は移動式での解析条件と同じである。
【００４２】
図６は、上記した解析条件で、車輪の転動をシミュレートしたときの車輪の中心軸の進行方向への移動速度の履歴、すなわち車輪の並進速度の履歴を示したものである。車軸には同一のトルクが常にかかるため、車輪の並進速度は直線的に増加している。かかる並進速度において、従来式のレールの全長を５ｍとした解析結果（従来式（５ｍ））及び１０ｍとした解析結果（従来式（１０ｍ））と移動式の解析結果とを比較すると、図６に示すような一致を得られることが判る。
【００４３】
図７に示すように、移動式での解析を従来式での解析と比較すると、解析に使用する要素数及び接点数が大幅に削減できる。なお、移動式の要素数及び接点数はダミーブロックも含めて算出している。また、車輪の並進速度が１０ｋｍ／ｈとなるまでをシミュレートした解析において、移動式での計算時間は５４０２３秒を要し、従来式（５ｍ）での１８８８８０秒及び従来式（１０ｍ）での５８３４４８秒と比較して、それぞれ０．２９倍及び０．０９倍と非常に小さい計算時間比である。つまり、移動式では、解析時間を大幅に短縮できる。また、レールの全長を長くすると計算時間比は更に小さくなるのである。
【００４４】
一方で、車輪の並進速度が１０ｋｍ／ｈとなったときにレールにかかるＺ方向の最大応力は、移動式で−７６８．０２５３ＭＰａ、従来式（５ｍ及び１０ｍ）で−７６８０．２５２ＭＰａとなった。つまり、高い解析精度を維持していた。
【００４５】
本実施例によれば、レール４２（５２）について接触部Ｐの位置を基準に所定の範囲に解析区間を限定してこれを移動させることで、解析時間を短縮できる。また、解析区間の移動に際しては、接触部Ｐの所定距離の移動毎に単位ブロックの１つ分だけずらすので、解析区間の移動の回数を低減できる。更に、単位ブロックには独立したメッシュモデルを与えて、例えば複数の単位ブロックに跨る要素についての処理のような煩雑な処理を回避し、解析区間を新たに定めるステップの処理を簡素化できる。また、単位ブロック毎のメッシュモデルを同一にすることで、解析区間を新たに定めるステップにおいてメッシュモデルを変化させるなどの処理を不要とし、メッシュモデルの付与の処理を更に簡素化できる。よって、レール４２（５２）のメッシュモデルを細かくして、解析の精度を高く保つことができ、その一方で、解析時間を短縮できる。
【００４６】
なお、解析区間の移動においては、解析区間を接触部Ｐの位置に応じて単位ブロックの１つ分だけ進行方向にずらせばよい。つまり、単位ブロックを解析区間に加えて解析区間から外すパターンを上記以外のパターンとすることもできる。例えば、前方に新たに１つの単位ブロックを加えて解析を行った後、後方から単位ブロックの１つを外して解析を行い、更にこれを繰り返すこともできる。予め定めたパターンに従い接触部Ｐの位置に応じて解析区間に加えられる単位ブロックの数を変更することも可能である。
【００４７】
また、レールを複数の単位ブロックに仮想分割した際に、レールの全長について予めメッシュモデルを作成しておいてもよい。メッシュモデルを与えられた単位ブロックを解析区間に選択するのである。上記した実施例との比較において、メッシュモデルの付与の処理をより簡素化して解析時間を短縮できる方法を選択することが好ましい。
【００４８】
ここまで本発明による代表的実施例及びこれに基づく変形例を説明したが、本発明は必ずしもこれらに限定されるものではなく、当業者であれば、本発明の主旨又は添付した請求項の範囲を逸脱することなく種々の代替実施例及び改変例を見出すことができるであろう。
【符号の説明】
【００４９】
３０        有限要素解析システム
３４        解析区間設定部
３５        解析部
４１、５１  車輪
４２、５２  レール
Ｐ          接触部
【図１】
0005902986.tif 000002
【図２】
0005902986.tif 000003
【図３】
0005902986.tif 000004
【図５】
0005902986.tif 000005
【図７】
0005902986.tif 000006
【図４】
0005902986.tif 000007
【図６】
0005902986.tif 000008