(57)【特許請求の範囲】
【請求項１】
コンピュータを用いた有限要素法に基づくプレス成形シミュレーション方法を用いて前記有限要素法の各要素節点における温度変化を求め、求めた各要素節点の温度変化に基づき、所定温度以上となっている高温領域部分を抽出し、該高温領域部分からブランク材が圧型内に入り込む方向線上で、ポンチプロファイルラインの外側にディスタンスブロックを配置することを特徴とするディスタンスブロックの配置方法。
【請求項２】
前記ポンチプロファイルラインの外側に配置したディスタンスブロックによって熱を逃がすことのできない前記ポンチプロファイルラインの内側にある前記高温領域部分には、その近傍に圧型リブを設けることを特徴とする請求項１に記載のディスタンスブロックの配置方法。
【請求項３】
前記プレス成形シミュレーション方法における前記各要素節点における温度変化は、ブランク材とプレス圧型の間の摺動熱、およびブランク材の塑性変形に起因する材料変形熱を算出することにより求めることを特徴とする請求項１または２に記載のディスタンスブロックの配置方法。
【請求項４】
前記摺動熱は、ブランク材とプレス圧型の接触部分における面圧、摩擦係数、および単位時間あたりの摺動距離から算出することを特徴とする請求項３に記載のディスタンスブロックの配置方法。
【請求項５】
前記摺動熱の算出時には、ブランク材とプレス圧型との接触部分における熱伝導による熱の移動を算出することを特徴とする請求項３または４に記載のディスタンスブロックの配置方法。
【請求項６】
前記材料変形熱は、相当塑性仕事と比例した量として近似することを特徴とする請求項３に記載のディスタンスブロックの配置方法。
【請求項７】
前記摺動熱および材料変形熱の算出時には、ブランク材からの大気放熱量を算出することを特徴とする請求項３〜６のいずれか１つに記載のディスタンスブロックの配置方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、ディスタンスブロックの配置方法に関する。
【０００２】
【従来の技術】
自動車用鋼板やその他の金属板材などをプレスして成形するプレス成形における形成不具合として、主なものに、寸法精度不良、われ不具合、しわ不具合などがある。近年数値解析技術の進歩により、これらの不具合を、コンピュータを用いたシミュレーションにより推定すること可能になってきている。例えば特開平１１−３１９９７１号公報では、有限要素法を用いたシミュレーションによりプレス成形時のわれやしわの発生をシミュレートする技術が開示されている。
【０００３】
【発明が解決しようとする課題】
しかしながら、このようなシミュレーションを実行して、その解析結果からプレス圧型を製造し、実際にプレス成形を行った場合、単発形成ではわれ、しわなどの不具合の発生しない圧型でも、量産時には形成変動により、われ、しわが発生することがあり、従来のシミュレーションだけでは、現場でさらに型の調整が必要であるといった問題があった。
【０００４】
本発明者による研究の結果、このような量産時における形成変動の主要因は、量産に伴う型の温度変化であることがわかってきた。
【０００５】
ところが、プレス成形シミュレーションでは、これまでプレス圧型の温度変化に着目して数値解析したものはなく、全くといっていいほど研究されていないのが現状である。
【０００６】
そこで、本発明の目的は、プレス成形シミュレーション方法によりプレス形成時におけるブランク材や圧型の温度上昇の原因となる熱発生を評価して、発熱分布の偏りのないプレス圧型を得ることができるディスタンスブロックの配置方法を提供することである。
【０００７】
【課題を解決するための手段】
本発明の目的は、下記する手段により達成される。
【０００８】
（１）コンピュータを用いた有限要素法に基づくプレス成形シミュレーション方法を用いて、前記有限要素法の各要素節点における温度変化を求め、求めた各要素節点の温度変化に基づき、所定温度以上となっている高温領域部分を抽出し、該高温領域部分からブランク材が圧型内に入り込む方向線上で、ポンチプロファイルラインの外側にディスタンスブロックを配置することを特徴とするディスタンスブロックの配置方法。
【０００９】
（２）前記ポンチプロファイルラインの外側に配置したディスタンスブロックによって熱を逃がすことのできない前記ポンチプロファイルラインの内側にある前記高温領域部分には、その近傍に圧型リブを設けることを特徴とする（１）に記載のディスタンスブロックの配置方法。
【００１０】
（３）前記プレス成形シミュレーション方法における前記各要素節点における温度変化は、ブランク材とプレス圧型の間の摺動熱、およびブランク材の塑性変形に起因する材料変形熱を算出することにより求めることを特徴とする（１）または（２）に記載のディスタンスブロックの配置方法。
【００１１】
（４）前記摺動熱は、ブランク材とプレス圧型の接触部分における面圧、摩擦係数、および単位時間あたりの摺動距離から算出することを特徴とする（３）に記載のディスタンスブロックの配置方法。
【００１２】
（５）前記摺動熱の算出時には、ブランク材とプレス圧型との接触部分における熱伝導による熱の移動を算出することを特徴とする（３）または（４）に記載のディスタンスブロックの配置方法。
【００１３】
（６）前記材料変形熱は、相当塑性仕事と比例した量として近似することを特徴とする（３）に記載のディスタンスブロックの配置方法。
【００１４】
（７）前記摺動熱および材料変形熱の算出時には、ブランク材からの大気放熱量を算出することを特徴とする（３）〜（６）のいずれか１つに記載のディスタンスブロックの配置方法。
【００１６】
【発明の効果】
本発明は請求項ごとに以下の効果を奏する。
【００１７】
請求項１記載の本発明によれば、プレス成形シミュレーションにより有限要素の各要素節点における温度変化を求め、求めた各要素節点の温度変化に基づき、所定温度以上となっている高温領域部分を抽出し、該高温領域部分からブランク材が圧型内に入り込む方向線上で、ポンチプロファイルラインの外側にディスタンスブロックを配置することとしたので、プレス成形シミュレーションにより解析された温度分布のうち、高温領域の熱をポンチプロファイルラインの外側に配置したディスタンスブロックによって逃がすことができるようになり、プレス成形時に、発熱分布の偏りがなく、熱による圧型変形や形成変動を抑えたプレス圧型を得ることが可能となる。
【００１８】
請求項２記載の本発明によれば、ポンチプロファイルラインの外側に配置したディスタンスブロックによって熱を逃がすことのできない前記ポンチプロファイルラインの内側にある前記高温領域部分の近傍に圧型リブを設けることとしたので、ポンチプロファイルラインの内側で、ディスタンスブロックの配置が難しい部分においてもリブを設けることで、熱を逃がし、プレス成形時の熱により圧型変形や形成変動を抑えることが可能となる。
【００１９】
請求項３記載の本発明によれば、温度変化を求めるために、ブランク材とプレス圧型の間の摺動熱、およびブランク材の塑性変形に起因する材料変形熱を算出することとしたので、プレス成形時の主な発熱原因による温度変化を適確に求めることが可能となる。
【００２０】
請求項４記載の本発明によれば、摺動熱をブランク材とプレス圧型の接触部分における面圧、摩擦係数、および単位時間あたりの摺動距離から算出することとしたので、増分型で解析を行う有限要素法によるシミュレーションに容易に導入することができる。
【００２１】
請求項５記載の本発明によれば、摺動熱の算出時にブランク材とプレス圧型の接触部分における熱伝導による熱の移動を算出することとしたので、接触部分からの放熱が考慮された正確な温度変化を求めることが可能となる。
【００２２】
請求項６記載の本発明によれば、材料変形熱を、相当塑性仕事と比例した量として近似することとしたので、増分型で解析を行う有限要素法によるシミュレーションに容易に導入することができる。
【００２３】
請求項７記載の本発明によれば、摺動熱および材料変形熱の算出時には、ブランク材からの大気放熱量を算出することとしたので、大気放熱が考慮された正確な温度変化を求めることが可能となる。
【００２５】
【発明の実施の形態】
以下、添付した図面を参照して、本発明の一実施の形態を説明する。
【００２６】
本実施の形態では、従来からある有限変形弾塑性有限要素法、特に動的陽解法を元にしたコンピュータによるシミュレーションにより、プレス圧型やブランク材（プレス材料）に加わる変位量、ひずみ量、および応力などを解析すると共に、本発明を適用してプレス圧型の温度変化を解析し、熱の発生領域を求めるものである。なお、シミュレーションは、後述する処理手順にしたがって作成されたプログラムをワークステーション、あるいはパソコンなどと称されているコンピュータが実行することにより行なわれる。また、このようなシミュレーションを行なうコンピュータをシミュレータと称する。
【００２７】
図１および図２は、本発明を適用したシミュレーションの処理手順を示すフローチャートである。
【００２８】
まず、シミュレーションに必要なデータの入力を行う。入力するデータは、例えばＣＡＤ装置などに記憶されている設計データをもとに、有限要素法解析に必要な有限要素メッシュデータ、プレス圧型の形状および質量、ブランク材の形状および質量、プレス圧型とブランク材との接触条件および支持条件などの幾何学境界条件、荷重、熱源などの力学境界条件などである（Ｓ１）。
【００２９】
続いて、変数の初期条件を設定する（Ｓ２）。ここでは、シミュレートするプレス圧型およびブランク材の座標、変位、ひずみ、応力、および温度などの初期値を設定する。
【００３０】
続いて、入力されたプレス圧型およびブランク材のデータから形状関数の設定し（Ｓ３）、質量マトリックスおよび減衰マトリックスを設定する（Ｓ４）。
【００３１】
そして、各種設定が済んだところで、増分計算の繰り返しとなる（Ｓ５）。増分量は、通常の動的陽解法においては、時間増分を用いるのであるが、プレス成形のシミュレーションでは、プレス圧型の移動によりブランク材が変形されて成形されるため、ここではプレス圧型の移動量を増分とし、あらかじめ決めた一定量を一増分量としている。また、増分計算の繰り返しの終了（後述するステップＳ２１）は、プレス圧型が下死点に到達した時点である。
【００３２】
増分計算では、まず、周知の方法（例えば、「非線形有限要素法」、１９９４年１２月２０日、コロナ社発行）により、材料構成式マトリックス（弾塑性構成式）の算出（Ｓ６）、ガウス求積法による節点内力の算出とペナルティ法による節点外力の算出（Ｓ７）、中央差分による変位算出（Ｓ８）、算出されたこの増分内での変位から、変位、ひずみ、応力の増分量の計算（Ｓ９）を順次実行して、有限要素節点ごとの変位量、ひずみ、および応力を求め、これらの値を更新する（Ｓ１０）。これにより、プレス圧型の変形量や、ひずみ、応力がシミュレートされる。
【００３３】
続いて、同じ増分計算の中で、本発明を適用した要素節点ごとの温度変化を求める計算に入る。
【００３４】
温度変化は、発熱量と放熱量を算出することにより求めることができる。発熱量としては、ブランク材とプレス圧型が接触する部分で発生する摺動熱と、ブランク材の塑性変形による発熱である。一方、放熱は、大気中への放熱と、プレス圧型と接触した部分からの熱伝導による放熱である。本実施の形態では、まず発熱量を算出した後、放熱量を算出して、ブランク材および圧型の温度を更新する。
【００３５】
まず、摺動熱を算出する。プレス圧型とブランク材との摺動熱算出のために、要素節点面圧を算出し（Ｓ１１）、続いて、算出された節点面圧を用いて節点摺動熱増分を算出する（Ｓ１２）。
【００３６】
摺動熱Ｔｓは、摩擦抵抗の仕事量と比例していると考えることができ、ブランク材の要素節点ごとに、摺動熱増分ΔＴｓは下記（１）式より求める。
ΔＴｓ＝ｆ（Ｐ，μ，ｌ）＝αＰｉμΔＬ    ・・・・・（１）
ただし、式中、Ｐｉは増分ステップｉにおける面圧（前記ステップＳ１２にて算出）、μはブランク材の摩擦係数、ΔＬは１増分における摺動距離、αは摩擦抵抗の仕事量が熱に変換される割合で、ブランク材を自動車車体用の鋼板とした場合、α＝０．７５〜０．９５を用いる。
【００３７】
なお、ステップＳ１２における面圧Ｐｉの算出は、例えばペナルティ法を用いて算出し、各増分ステップにおいて、ブランク材の有限要素が圧型モデルの方へ食い込む量を計算し、この値にペナルティ数を乗じることで、法線方向の接触力として面圧を算出する。
【００３８】
続いて塑性変形に伴う発生熱を算出する。このためには、相当応力および相当塑性ひずみ速度を算出し（Ｓ１３）、これを用いて有限要素の素材（ブランク材）の変形熱増分を算出する（Ｓ１４）。
【００３９】
ここで要素素材変形熱増分ΔＴｄは、相当塑性仕事と比例していると考えることができる。そこで、各材料要素に対して下記（２）式により、変形熱増分ΔＴｄを近似して求める。
【００４０】
【数１】
0004696345.tif 000002
【００４１】
ここで、式中、βは塑性仕事の熱に変換される割合で、β＝０．８５〜０．９５とした。また、式中、Ｗpは塑性仕事量であり、σは相当応力、
【００４２】
【外１】
0004696345.tif 000003
【００４３】
である。
【００４４】
続いて、求めた摺動熱増分ΔＴｓと変形熱増分ΔＴｄを各要素節点に分配する（Ｓ１５）。これにより、ブランク材の発熱による温度上昇が要素節点ごとに求まる。
【００４５】
次に、プレス形成中における放熱量を算出する。
ブランク材は、成形中に型との接触離脱を繰り返す。本実施の形態では、要素節点ごとに熱量計算を行うため、各要素節点の有効面積を定義しておく必要がある。
【００４６】
要素節点ごとの有効面積Ｓは、下記（３）式により、節点を含む各要素について、要素節点周りにおける形状関数の数値積分を行い、各要素の総和をとることで求める。
【００４７】
【数２】
0004696345.tif 000004
【００４８】
ここで，式中、ｓは各要素の接触面積、ｉは各要素番号、ξ，ηは局所座標系、Ｆ（ξ，η）は形状関数（ステップＳ３で定義）である。
【００４９】
また、ブランク材と型との接触した部分における熱伝導によるブランク材から型への熱移動速度は、下記（４）式により求める。
【００５０】
【数３】
0004696345.tif 000005
【００５１】
ここで、ｋはブランク材の熱伝導率、Δｘは熱源間距離、ΔＴは２つの熱源間の温度差である。
【００５２】
本実施の形態では、熱源、すなわち熱の発生部分をブランク材の内部と仮定して、Δｘを板厚の１／２とし、接触時間を熱移動速度にかけることで、下記（５）式により、各増分ステップでの熱移動量ΔＱを計算する。
【００５３】
【数４】
0004696345.tif 000006
【００５４】
ここで、式中、Δｔは１増分あたりの時間、ｗはブランク材の厚さである。
【００５５】
さらに、熱移動量から温度変化を求めるためには下記（６）式による計算を行う。
ΔＴ＝ΔＱ／（ｃｍ）    ・・・・・（６）
ここで、式中、ｃは比熱、ｍは質量である。
【００５６】
本実施形態において、比熱ｃは、ブランク材とプレス圧型をいずれも炭素鋼とした場合４８０Ｊ／ｋｇＫを使用する。また、質量ｍは、ブランク材においては、その質量であるが、プレス圧型は、ブランク材と同板厚の質量を仮定する。そのため、プレス圧型の質量ｍは、下記（７）式より求める。
【００５７】
ｍ＝Ｖｔ・ρ＝Ｓｔ・Ｗ・ρ    ・・・・・（７）
ここで、式中、Ｖｔは接触要素体積、ρはプレス圧型の密度、Ｓｔは接触要素面積である。なお、雰囲気温度は一定として取り扱う。
【００５８】
放熱量の算出には、これら各式を用いて、まず、プレス圧型と接触していない非接触材料節点の放熱量を算出する（Ｓ１６）。そしてこの非接触材料節点、すなわち、空気中への放熱量を先に算出した発熱量に加算する（実質的には大気中への放熱であるので、（６）式の算出結果はマイナスとなる）。すなわち、材料節点温度増分＝摺動熱増分＋材料変形熱増分＋放熱量増分を算出する（Ｓ１７）。
【００５９】
続いて、プレス圧型と接触しているブランク材の要素節点における吸放熱量による温度変化を算出し（Ｓ１８）、続いて、プレス圧型の大気放熱量による温度変化を算出する（Ｓ１９）。
【００６０】
そして、この増分ステップにける各要素節点の温度変数、すなわち、ブランク材の温度と、プレス圧型の温度を更新して（Ｓ２０）、プレス圧型が下死点へ到達するまでステップＳ５からＳ２１を繰り返し実行する（Ｓ２１）。これにより、プレス成形時のプレス圧型とブランク材の温度変化を有限要素法によるシミュレーションにより求めることができる。
【００６１】
次に、上述したプレス成形時の温度変化のシミュレーション結果を用いて、ディタンスブロック（ディスタンスピースなどとも称されている）をレイアウトする方法について説明する。
【００６２】
ディスタンスブロックは、プレス圧型の温度増分が大きな部分に配置することで、板厚クリアランス変動を抑え、成形変動によるわれやしわの発生を抑制するためのものである。したがって、ディスタンスブロックとしては、放熱の良い形状、例えばフィン形状のものを配置したり、あるいは熱容量の大きな部材を設けることが好ましい。
【００６３】
図３は、上記のシミュレーションによる熱解析の結果を用いたディスタンスブロックの配置方法の手順を示すフローチャートである。
【００６４】
まず、シミュレーションにプレス成形時の熱解析の結果と、初期ブランク形状のデータを入力する（Ｓ５１）。このとき入力されるデータを図４に示す。図４は、解析結果の温度分布図１とブランク材の初期外形線２（プレス成形前の初期形状である）。温度分布図１はシミュレーションの結果として、プレス形成後のブランク材の熱分布を示すもので、図４では、温度分布を示す等高線として示し、特に、温度の高くなっている部分をハッチングの密度の違い（ハッチング密度の密な方がより温度が高い）により示したが、実際のシミュレーションにおいては、例えば温度分布を色で表すなど、ブランク材の温度分布がよりわかりやすいように表示される。
【００６５】
続いて、図５に示すように、入力されたデータから一定温度以上となっているホットポイント３（高温領域部分）を抽出する（Ｓ５２）。このホットポイントの抽出は、あらかじめ閾値温度を設定しておき、その閾値温度以上となっている発熱点をホットポイント３としてシミュレータに抽出させる。なお、自動抽出させる代わりに、設計者などが温度分布図１を見て問題となそうな発熱点をホットポイント３として設定してもよい。
【００６６】
続いて、図６に示すように、ホットポイント３の流入線４上、ポンチプロファイルライン５の外側にディスタンスブロック６を配置するようにディスタンスブロック６のレイアウトを決定する（Ｓ５３）。発熱点、すなわち解析結果におけるホットポイント３の近傍では、発熱によって型変形が発生し、上下型間のクリアランスが変化する。そこで、発熱点から流入方向にディスタンスブロック６を配置することで、この部分の熱を下げ、型変形の発生を防止するものである。なお、流入線４とは、ブランク材が圧型内に入り込む方向線である。
【００６７】
続いて、配置されたディスタンスブロック６の効果を確かめるために熱発生点の再評価を行う（Ｓ５４）。この再評価は、図７に示すように、配置したディスタンスブロック６によりクリアランスが安定する効果範囲７（半径任意設定）内にあるホットポイント３を評価対象から外す。しかし、ダイフェース上の発熱範囲が広い場合、単独のディスタンスブロック６では発熱による影響を全て解決できないので、ここで、温度分布図１内における発熱部分を再評価して、新たに問題となるような高温域があればホットポイント１３として設定する。この新たなホットポイント１３の設定は、設計者などが温度分布図１を見て問題となりそうな発熱点を見出して、その部分をホットポイント１３として設定してもよいし、先のホットポイントを抽出する閾値を下げて、シミュレータに自動抽出させてもよい。なお、１つのディスタンスブロック６によりクリアランスが安定する効果範囲７は、型構造によって変化するため、ここでは、その半径を設計者などが任意（例えばこれまでの経験値など）に設定するものとする。
【００６８】
続いて、新たに設定したホットポイント１３に対して、図８に示すように、前記ステップＳ５３と同様に、新たにディスタンスブロック１６を配置する（Ｓ５５）。このディスタンスブロック１６によってクリアランスが安定する効果範囲１７により、新たに設定したホットポイント１３を消去することができる。
【００６９】
そして、ダイフェース上のポンチプロファイルライン５の外側に配置するディスタンスブロック６および１６によって消去できるホットポイント３および１３がなくなるまで、前記ステップＳ５４〜Ｓ５５を繰り返し実行して、ディスタンスブロック１６を配置する（Ｓ５６）。
【００７０】
続いて、ディスタンスブロック６および１６の配置が終了したなら、図９に示すように、圧型リブ必要線９および必要点１０を入力する（Ｓ５７）。これには、まず、ディスタンスブロックおよびネストピンなどのレイアウト要件から圧型外形線８を設定し、成形荷重などの要件から圧型リブ必要線（必ずリブを設けなくてはならない設計線）９を定義し、さらに、先にレイアウトしたディスタンスブロック６および１６の配置個所にも流入方向に、圧型外形線からポンチプロファイルラインまで、圧型リブ必要線９を設定する。
【００７１】
また、ポンチプロファイルライン５内のホットポイント３（ステップＳ５４〜Ｓ５６においてディスタンスブロック６の配置により消去できないホットポイント３）のうち、リブ必要線９が通らないホットポイントを圧型リブ必要点１０と定義し、必ずリブが通らなくてはならない設計線とする。これにより、ポンチプロファイルライン５内のホットポイント３に起因する型変形を抑制し、圧型の熱容量の増大を図る。
【００７２】
続いて、リブ必要点１０を再評価し、リブ必要点１０近傍には、圧型リブ必要線１９を設け、これにより、この部分での型変形などが抑制されるものとして、リブ必要点１０を評価対象から外す（Ｓ５８）。結果として、図１０に示すように、圧型リブ必要箇所を線分のみで表現する。
【００７３】
続いて、図１１に示すように、圧型リブ必要線を単に線分として設定した段階において、ところどころ他のリブと接続されていないような部分を、それぞれ他のリブと交差するまで追加延長して（Ｓ５９）、すべての圧型リブ必要線９および１９を接続し（Ｓ６０）、さらに、圧型リブ必要線９および１９が設計基準以上に離れている個所に新たに圧型補強リブ１１を追加する（Ｓ６１）。
【００７４】
このようにして出来上がったディスタンスブロックおよびリブ必要線のレイアウト図面データを圧型設計ＣＡＤなどへ出力して（Ｓ６２）、処理を終了する。
【００７５】
以上のように、本実施の形態では、有限要素法によるシミュレーションにより、プレス成形時の温度変化を求め、その結果に基づいてディスタンスブロックおよびリブを配置することで、ディスタンスブロックおよびリブにより、圧型の表面積および熱容量の増大を図り、発熱点から熱を効果的に逃がし、また、リブにより圧型の強度を増大させて、熱による型変形を抑制することができる。
【図面の簡単な説明】
【図１】  本発明を適用した実施の形態におけるシミュレーションの処理手順を示すフローチャートである。
【図２】  図１に続くシミュレーションの処理手順を示すフローチャートである。
【図３】  本発明を適用した実施の形態におけるディスタンスブロックの配置方法の手順を示すフローチャートである。
【図４】  前記シミュレーションの結果である温度分布図とブランク材の初期外形を示す図面である。
【図５】  前記シミュレーションの結果である温度分布図、ブランク材の初期外形、およびホットポイントを示す図面である。
【図６】  前記シミュレーションの結果である温度分布図、ブランク材の初期外形、ホットポイント、およびディスタンスブロックの配置を示す図面である。
【図７】  前記ディスタンスブロックの効果範囲を示す図面である。
【図８】  前記ディスタンスブロックの効果範囲から外れたホットポイントの設定を示す図面である。
【図９】  前記ディスタンスブロックの効果範囲から外れたホットポイントに新たに配置されたディスタンスブロックを示す図面である。
【図１０】  前記ディスタンスブロックと圧型リブの配置を示す図面である。
【図１１】  最終的に出力されるディスタンスブロックと圧型リブの配置を示す図面である。
【符号の説明】
１  温度分布図
２  ブランク材初期外形線
３、１３  ホットポイント
４  流入線
５  ポンチプロファイルライン
６、１６  ディスタンスブロック
７  ディスタンスブロックの効果範囲
８  圧型外形線
９、１９  圧型リブ必要線
１０  圧型リブ必要点
１１  圧型補強リブ
【図１】
0004696345.tif 000007
【図２】
0004696345.tif 000008
【図３】
0004696345.tif 000009
【図４】
0004696345.tif 000010
【図５】
0004696345.tif 000011
【図６】
0004696345.tif 000012
【図７】
0004696345.tif 000013
【図８】
0004696345.tif 000014
【図９】
0004696345.tif 000015
【図１０】
0004696345.tif 000016
【図１１】
0004696345.tif 000017