(57)【要約】　　　（修正有）【課題】溶接ビードの測定区間の全体を通じて、溶接ビードの品質を正しく評価する。【解決手段】溶接ビードの品質評価方法は、溶接ビードの長さ方向に所定間隔で溶接ビードの高さを測定して高さデータを取得し、高さデータから高さの分散値を演算し、高さの分散値に基づいて、溶接ビードの品質を評価するものである。【選択図】　図３
2004053427.tif 000002
【特許請求の範囲】
【請求項１】
溶接ビードの長さ方向に所定間隔で溶接ビードの高さを測定して高さデータを取得し、高さデータから高さの分散値を演算し、高さの分散値に基づいて、溶接ビードの品質を評価する溶接ビードの品質評価方法。
【請求項２】
溶接ビードの長さ方向に所定間隔で溶接ビードの幅を測定して幅データを取得し、幅データから幅の分散値を演算し、幅の分散値に基づいて、溶接ビードの品質を評価する溶接ビードの品質評価方法。
【請求項３】
溶接ビードの長さ方向に所定間隔で溶接ビードの高さを測定し、各測定か所において高さが最大となるところをビードの中心位置であると認定することにより、ビードの中心位置に係わる位置データを取得し、位置データからビードの中心位置の分散値を演算し、ビードの中心位置の分散値に基づいて、溶接ビードの品質を評価する溶接ビードの品質評価方法。
【請求項４】
溶接ビードの長さ方向に所定間隔で溶接ビードの幅方向の大きさおよび位置を測定し、溶接ビードの幅方向における一方の側に最大に張出したか所と、その他方の側に最大に張出したか所との溶接ビード幅方向の距離を測定し、これをビード最大余盛幅とし、ビード最大余盛幅と溶接ビードの幅方向の大きさの比に基づいて、溶接ビードの品質を評価する溶接ビードの品質評価方法。
【請求項５】
溶接ビードの長さ方向に所定間隔で溶接ビードの高さを測定して高さデータを取得し、高さデータをフーリエ変換し、フーリエ変換後の高さデータの平均値を求め、フーリエ変換後の高さデータの平均値に基づいて、溶接ビードの品質を評価する溶接ビードの品質評価方法。
【請求項６】
溶接ビードの長さ方向に所定間隔で溶接ビードの高さを測定して高さデータを取得し、高さデータをフーリエ変換し、フーリエ変換後の高さデータの分散値を求め、フーリエ変換後の高さデータの分散値に基づいて、溶接ビードの品質を評価する溶接ビードの品質評価方法。
【請求項７】
請求項１〜６に記載の溶接ビードの品質評価方法を組合せる溶接ビードの品質評価方法。
【請求項８】
基準面上の溶接ビードにスリット光を投影し、投影したスリット光をイメージブレーン上に結像させる溶接ビードの形状測定方法において、基準面上の溶接ビードの有無により、イメージブレーン上においてスリット光を移動させ、移動するスリット光の移動量を測定し、測定した移動量を用い、三角測量の原理に基づいて、基準面上における溶接ビードの３次元的データを取得することを特徴とする溶接ビードの形状測定方法。
【請求項９】
溶接ビードにスリット光を投影し、投影したスリット光を撮影して溶接ビードの形状に係わるデータを取得する溶接ビードの形状測定方法において、取得したデータ上において溶接ビードの高さが最大となるところを認定するとともに、同データを微分して微分値を取得し、溶接ビードの高さが最大となるところ境界として、取得した微分値を２つの領域に区分し、各領域内において、微分値の平均を求め、微分値の平均となるところから標準偏差分だけ移動したところを幅の端であると認定することにより、溶接ビードの幅を取得する溶接ビードの形状測定方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、例えば、鉄板同士の突合溶接において、溶接ビードの接合が強くて高品質なのか、接合が弱くて低品質なのかの品質を評価する方法に関する。
【０００２】
【従来の技術】
この種の方法としては、例えば、特開平５−８７５３９号公報に開示されているように、溶接ビードにスリット光を投影し、投影したスリット光を撮影して溶接ビードの形状に係わるデータを取得する溶接ビードの形状測定方法において、溶接ビードの長さ方向の測定区間に複数の測定点を設定し、隣合う測定点におけるデータ同士の差分を求め、求めた差分が設定値以下であるか、否かを判定し、溶接ビードの測定区間にわたって、設定値以下であると判定された判定結果の回数に基づいて、溶接ビードの品質を評価するものが知られている。
【０００３】
【発明が解決しようとする課題】
上記方法では、品質評価の基準となるデータが隣合う測定点におけるデータ同士の差分に基づいている。そのため、溶接ビードの測定区間の全体を通じて測定データに大きい変化が無い場合は、問題が無い。ところが、例えば、測定区間の一端から他端にかけて、測定データが漸次変化し、測定区間の全体では、測定データが大きく変化してしまうような場合は、品質を正しく評価できないという問題点がある。
【０００４】
この発明の目的は、溶接ビードの測定区間の全体を通じて、溶接ビートの品質を正しく評価することのできる溶接ビードの品質評価方法およびそのための形状測定方法を提供することにある。
【０００５】
【課題を解決するための手段】
この発明による第１の溶接ビードの品質評価方法は、溶接ビードの長さ方向に所定間隔で溶接ビードの高さを測定して高さデータを取得し、高さデータから高さの分散値を演算し、高さの分散値に基づいて、溶接ビードの品質を評価するものである。
【０００６】
この発明による第１の溶接ビードの品質評価方法では、高さの分散値が小さければ、高さが揃っていると評価され、高さの分散値が大きければ、高さが揃っていないと評価される。したがって、溶接ビードの測定区間の全体を通じて、ビードの高さが揃っているか否かを評価することができる。
【０００７】
この発明による第２の溶接ビードの品質評価方法は、溶接ビードの長さ方向に所定間隔で溶接ビードの幅を測定して幅データを取得し、幅データから幅の分散値を演算し、幅の分散値に基づいて、溶接ビードの品質を評価するものである。
【０００８】
この発明による第２の溶接ビードの品質評価方法では、幅の分散値が小さければ、幅が揃っていると評価され、幅の分散値が大きければ、幅が揃っていないと評価される。したがって、溶接ビードの測定区間の全体を通じて、ビードの幅が揃っているか否かを評価することができる。
【０００９】
この発明による第３の溶接ビードの品質評価方法は、溶接ビードの長さ方向に所定間隔で溶接ビードの高さを測定し、各測定か所において高さが最大となるところをビードの中心位置であると認定することにより、ビードの中心位置に係わる位置データを取得し、位置データからビードの中心位置の分散値を演算し、ビードの中心位置の分散値に基づいて、溶接ビードの品質を評価するものである。
【００１０】
この発明による第３の溶接ビードの品質評価方法では、ビードの中心位置の分散値が小さければ、ビードの中心位置が標準位置からずれていないと評価され、ビードの中心位置の分散値が大きければ、ビードの中心位置が標準位置からずれていると評価される。したがって、溶接ビードの測定区間の全体を通じて、ビードの中心位置が標準位置からずれているか否かを評価することができる。
【００１１】
この発明による第４の溶接ビードの品質評価方法は、溶接ビードの長さ方向に所定間隔で溶接ビードの幅方向の大きさおよび位置を測定し、溶接ビードの幅方向における一方の側に最大に張出したか所と、その他方の側に最大に張出したか所との溶接ビード幅方向の距離を測定して、これをビード最大余盛幅とし、ビード最大余盛幅と溶接ビードの幅方向の大きさの比に基づいて、溶接ビードの品質を評価するものである。
【００１２】
この発明による第４の溶接ビードの品質評価方法では、ビード最大余盛幅と溶接ビードの幅方向の大きさの比が小さいと、ビード最大余盛幅が小さいと評価され、ビード最大余盛幅と溶接ビードの幅方向の大きさの比が大きいと、ビード最大余盛幅が大きいと評価される。したがって、溶接ビードの測定区間の全体を通じて、ビード最大余盛幅の大小を評価することができる。
【００１３】
この発明による第５の溶接ビードの品質評価方法は、溶接ビードの長さ方向に所定間隔で溶接ビードの高さを測定して高さデータを取得し、高さデータをフーリエ変換し、フーリエ変換後の高さデータの平均値を求め、フーリエ変換後の高さデータの平均値に基づいて、溶接ビードの品質を評価するものである。
【００１４】
この発明による第５の溶接ビードの品質評価方法では、フーリエ変換後の高さデータの平均値が小さければ、ビードの波が小さいと評価され、フーリエ変換後の高さデータの平均値が大きければ、ビードの波が大きいと評価される。したがって、溶接ビードの測定区間の全体を通じて、ビードの波の大小を評価することができる。
【００１５】
この発明による第６の溶接ビードの品質評価方法は、溶接ビードの長さ方向に所定間隔で溶接ビードの高さを測定して高さデータを取得し、高さデータをフーリエ変換し、フーリエ変換後の高さデータの分散値を求め、フーリエ変換後の高さデータの分散値に基づいて、溶接ビードの品質を評価するものである。
【００１６】
この発明による第６の溶接ビードの品質評価方法では、フーリエ変換後の高さデータの分散値が小さければ、ビードの波の大きさが揃っていると評価され、フーリエ変換後の高さデータの分散値が大きければ、ビードの波が大きさが揃っていないと評価される。したがって、溶接ビードの測定区間の全体を通じて、ビードの波の大小が揃っているか否かを評価することができる。
【００１７】
この発明による第７の溶接ビードの品質評価方法は、上記第１〜第６の発明の溶接ビードの品質評価方法を組合せるものである。
【００１８】
この発明による第１の溶接ビードの形状測定方法は、基準面上の溶接ビードにスリット光を投影し、投影したスリット光をイメージブレーン上に結像させる溶接ビードの形状測定方法において、基準面上の溶接ビードの有無により、イメージブレーン上においてスリット光を移動させ、移動するスリット光の移動量を測定し、測定した移動量を用い、三角測量の原理に基づいて、基準面上における溶接ビードの３次元的データを取得することを特徴とするものである。
【００１９】
この発明による第１の溶接ビードの形状測定方法では、取得した３次元的データが現実の数値と一致する。したがって、３次元的データを様々なデータに容易に加工することができる。
【００２０】
この発明による第２の溶接ビードの形状測定方法は、溶接ビードにスリット光を投影し、投影したスリット光を撮影して溶接ビードの形状に係わるデータを取得する溶接ビードの形状測定方法において、取得したデータ上において溶接ビードの高さが最大となるところを認定するとともに、同データを微分して微分値を取得し、溶接ビードの高さが最大となるところ境界として、取得した微分値を２つの領域に区分し、各領域内において、微分値の平均を求め、微分値の平均となるところから標準偏差分だけ移動したところを幅の端であると認定することにより、溶接ビードの幅を取得するものである。
【００２１】
この発明による第２の溶接ビードの形状測定方法では、ノイズ、雑音等の影響を受けることなく、取得したデータから溶接ビードの幅を確実に認定することができる。
【００２２】
【発明の実施の形態】
この発明の実施の形態を図面を参照してつぎに説明する。
【００２３】
図１を参照すると、評価の対象となるワーク１１が示されている。ワーク１１は、開先同士を突合せた２枚の鉄板２１と、開先を埋めて盛上るように形成された溶接ビード２２とよりなる。
【００２４】
ビード２２の形状の特徴を表すものとして、高さＨおよび幅Ｗがある。高さＨとは、鉄板２１表面からビード２２頂点までの距離である。幅Ｗとは、ビード２２の一方の側における盛上り開始点から、その他方の側における盛上り終了点までの距離である。
【００２５】
ビード２２の高さＨおよび幅Ｗは、スリット光投影法によって測定される。そのための測定装置３１が図２に示されている。測定装置３１は、レーザ４１およびＣＣＤカメラ４２と、これらをビード２２の長さ方向と平行に移動させる移動手段４３とよりなる。
【００２６】
図３は、測定装置３１の光学的配置を示すものである。垂直基準面Ｓ上にビード２２の長さ方向が水平になるようにワーク１１が置かれている。基準面Ｓとそれぞれ相対するようにレーザ４１およびカメラ４２が並べられている。レーザ４１からはスリット光が基準面Ｓに向かって投射される。このときのスリット光は、基準面Ｓに対して所定の角度をもった垂直面に含まれている。カメラ４２は、基準面Ｓと直交させられたレンズ５１および基準面Ｓと平行なイメージプレーン５２を有している。基準面Ｓ上に投影された光はレンズ５１を通ってイメージプレーン５２上で結像させられる。
【００２７】
基準面Ｓ上にワーク１１が無いときに、基準面Ｓ上に投影されたスリット光がＰ１で示されている。基準面Ｓ上にワーク１１が有るときに、ワーク１１表面に投影されたスリット光がＰ２で示されている。これら２つのスリット光Ｐ１、Ｐ２は、イメージプレーン５２上における２つの位置Ｐ３、Ｐ４で結像する。ワーク１１の有無により、イメージプレーン５２上において、スリット光が２つの位置Ｐ３、Ｐ４間で移動し、これを、移動量δとして表される。この移動量δを、三角測量の原理に適用すると、基準面Ｓ上におけるビード２２の３次元データが取得される。
【００２８】
基準面Ｓ上におけるビード２２の位置（Ｘ、Ｙ、Ｚ）は、つぎの（１）〜（３）式で表される。位置（Ｘ、Ｙ、Ｚ）は、現実の空間における数値と一致している。
【００２９】
【数式１】
2004053427.tif 000003
ここに、（ｕ、ｖ）：イメージプレーン５２上の座標
ｈ：基準面Ｓおよびレンズ５１間の距離
ｄ：レーザ４１およびレンズ５１間の距離
ｌ：レンズ５１およびイメージプレーン５２間の距離である。
【００３０】
レーザ４１およびカメラ４２をビード２２の長さ方向と平行に移動させながら、複数の所定の測定位置でスリット光の投影画像を取得する。取得した投影画像には、ノイズ、雑音等の影響を無くすための様々な画像処理を施し、正確なスリット光の投影位置を抽出する。各測定点において抽出したスリットの光投影位置に基づいて、移動量δを求め、これを（１）〜（３）式に適用することにより、ビード２２の３次元情報を取得する。取得したビード２２の３次元情報の例を、図４に示す。
【００３１】
つぎに、取得したビード２２の３次元情報から、ビード２２の形状の特徴を表す高さＨおよび幅Ｗの情報を取得する。
【００３２】
図５に、ビード２２の高さ情報の求め方が図解されている。各測定点おいて、測定データの高さＨｍａｘが最大となっている地点を求める。その地点を中心として、左右両側に一定間隔離れた地点の高さＨｌ、Ｈｒを求め、これらの高さＨｌ、Ｈｒの平均値を最大高さＨｍａｘから引いた値をその測定点におけるビード２２の高さＨとする。すなわち、
Ｈ＝Ｈｍａｘ−（Ｈｌ−Ｈｒ）／２　　　　　（４）
である。
【００３３】
ビード２２幅Ｗの求め方は、図６に示されている。ここでは、ビード２２の測定データが正規分布をなすものと見なしている。測定データを微分して微分値のデータを取得する。先に求めたビード２２の高さＨが最大となる地点を中心として、左右２つの領域Ａｌ、Ａｒに分ける。左側の領域Ａｌでは微分値は正成分となり、右側の領域Ａｒでは負成分となっている。左側の領域Ａｌにおいて、正成分の微分値の平均Ｍｌを求め、その平均Ｍｌから左側に標準偏差分εだけ移動した地点をビード２２の幅Ｗ方向左側の端であると認定する。右側の領域においても同様に、負成分の微分値の平均Ｍｒから右側に標準偏差分εだけ移動した地点をビード２２の幅Ｗ方向右側の端であると認定する。このようにして得られた２つの地点間の距離をビード２２の幅Ｗであると認定する。
【００３４】
以上により求められたビード２２の高さＨおよび幅Ｗのデータの一例が図７および図８にそれぞれ示されている。
【００３５】
ビード２２の品質を評価するためには、つぎの６つの特徴量を抽出する。
【００３６】
（１）　高さの分散値
（２）　幅の分散値
（３）　ビード中心位置の分散値
（４）　ビード最大余盛幅に対するビード幅の比の加算値
（５）　高さデータの離散フーリエ変換後の平均値
（６）　高さデータの離散フーリエ変換後の分散値
高さＨおよび幅Ｗの分散値は、実際の評価においてビード２２の高さＨや幅Ｗの全体的なバランス等が評価されることから用いられる特徴量である。すなわち、溶接ビード２２の高さＨや幅Ｗが整ったものとなっているかを評価する場合、良好なものは高さＨおよび幅Ｗがほぼ一定であり、そのため、良好なものでは分散値が低くなり、逆に粗悪なものは分散値が高くなると推測できる。
【００３７】
ビード中心位置Ｃとは、図９に示すように、先に求めた測定データの高さＨが最大となっている地点のことである。良好なビード程、ビード中心位置Ｃはほぼ一直線上に並ぶこととなって、分散値は低くなり、粗悪なビード程、ビード中心位置Ｃは直線からはずれて大きく蛇行し、分散値は高くなると推測できる。
【００３８】
以上、３種類の分散値は、先に取得したビード２２の高さＨおよび幅Ｗのデータから、つぎの（５）式によって求められる。
【００３９】
【数式２】
2004053427.tif 000004
ここに、σは分散値、Ｎはサンプル数、ｍは平均値、Ｐｉはサンプル値である。
【００４０】
ビード最大余盛幅Ｗｍａｘについて、つぎのように定義する。図１０（ａ）には蛇行するビード２２、図１０（ｂ）には部分的に幅広部をもつビード２２がそれぞれ示されている。これらのビード２２の形態を考慮して、ビード２２の幅方向における上側に最大に張出したか所と、その下側に最大に張出したか所との上下方向の距離を測定し、これをビード最大余盛幅Ｗｍａｘと定義する。ビード最大余盛幅Ｗｍａｘに対するビード幅Ｗの比を各測定点において算出し、全ての測定点で算出した比を全て加算したものを、ビード最大余盛幅Ｗｍａｘに対するビード幅Ｗの比の加算値とする。この加算値が小さい程、ビード最大余盛幅Ｗｍａｘが整っていると考えられる。
【００４１】
高さデータの離散フーリエ変換後の平均値および分散値は、ビードピッチに着目したものである。ビードピッチとはビード２２を形成する波の間隔を指す。溶接を行う際、溶接棒を溶接部位に対し一定の速度で動かすことができればビードピッチが一定に保たれると考えられる。ビードピッチが一定であれば、溶接部位に均等に熱や圧力が加えられ、溶接の強度は高くて、高品質の溶接であるといえる。逆に、ビードピッチがまちまちであると、熱や圧力が均等に加わらず、強度の弱い溶接と考えられる。要するに、ビードピッチに周期性が見られる程、溶接が高品位であると推察される。
【００４２】
図１１（ａ）に示す高さＨのデータから、これらの平均値を求め、高さＨのデータから平均値を引いた平均値０の離散的な波のデータを取得する。これを、図１１（ｂ）に示す。この離散データに対して、波の周期性を特徴としてとらえるために離散フーリエ変換を施した結果が図１１（ｃ）に示されている。図１１（ｃ）は横軸に周波数が、縦軸にスペクトルがそれぞれとられている。
【００４３】
離散フーリエ変換は、時間的な離散データ（ａ０、ａ１、…、ａＮ）を周波数データ（Ａ０、Ａ１、…、ＡＮ　）に変換するものである。その変換式をつぎに示す。
【００４４】
【数式３】
2004053427.tif 000005
図１１（ｃ）に示す高さＨのデータの離散フーリエ変換後の値に基づいて、そのスペクトルの平均値および分散値を求める。平均値が小さい程、波が無いことを意味し、分散値が小さい程、波のピッチが一定であることを意味する。
【００４５】
つぎに、本願発明による評価方法の有効性を検証するための実験を行ったので、これを説明する。
【００４６】
ワーク１１として、４８枚のサンプルを用いた。これらのサンプルは、実際に行われた溶接技能試験に供せられかつ試験官によって評点されたものである。評点は、満点２０点で行われている。
【００４７】
４８枚のサンプルに対し、測定装置３１を１ｍｍずつ移動させながら、１サンプルあたり所定枚数ずつの画像を取得した。取得した画像に対し、上記データ処理を施し、それぞれのサンプルに対し、６つの特徴量をそれぞれ抽出した。
【００４８】
溶接技能試験における評点と、６つの特徴量との対比を図１２（ａ）〜（ｆ）に示す。図１２（ａ）〜（ｆ）の横軸には溶接技能試験における評点が、図１２（ａ）〜（ｆ）の縦軸には抽出した６つの特徴量の値がサンプル毎に示されている。図１２（ａ）〜（ｆ）中、点線は各評価点での特徴量の平均値を示している。図１２（ａ）〜（ｆ）の中でも、図１２（ａ）、（ｅ）および（ｆ）に示すビード高さＨの分散値、フーリエ変換後の平均値および分散値については、評点の大きいもの程、それらの値が低くなっていることが顕著に見てとれる。
【００４９】
さらに、６つの特徴量を総合して評価するケースについて説明する。
【００５０】
４８枚のサンプルからそれぞれ抽出した６つずつの特徴量は、その平均値や分散値が低いもの程、評価が良好である。そこで、６つの特徴量のそれぞれにおいて、４８枚のサンプルの値のうち、最小の値を良好な基準値とし、最大の値を粗悪な基準値として設定する。そして、これらの基準値から、評価対象の個々のサンプルの特徴量がどれだけ離れているか、その度合を定量評価する。すなわち、良好または粗悪な基準値に評価対象の特徴量を加減算したものと、良好または粗悪な基準値との比を求める。さらに、良好な基準値からの離れ度合いを減点、粗悪な基準値からの離れ度合いを加点とし、これらの点数を合計したものを最終的な定量評価の値とする。加点要素をＡｉ、減点要素をＤｉとすると、評価値Ｓｃｏｒｅは、つぎの式によって表される。ここに、ｉはサンプル数４８であり、Ｅは特徴量数６である。
【００５１】
【数式４】
2004053427.tif 000006
評価値は、比の値であるため、各特徴量の最高点は、１．０である。一方、試験官による評点は最高が２０であるため、評価値と評点を比較し易くするために、評価値Ａｏｌｄをつぎの式を用いて補正し、その補正値Ａｎｅｗをグラフにして、図１３に示す。
【００５２】
Ａｎｅｗ＝２０Ａｏｌｄ／Ｅ　　　　　　　　　　（１１）
図１３（ａ）では加点を黒○印とし、減点を＋印として別々に示し、これらを合計したものが図１３（ｂ）に示されている。若干のばらつきは見られるが、評価値と評点の間に関連性があることが見られるであろう。
【００５３】
以上、６つの特徴量に基づいて説明したが、これ以外の特徴量を抽出することも可能である。何故ならば、上記したように、ビード２２の３次元情報を現実の数値として取得しており、これは、どのようにも加工することができるため、必要に応じて、任意の情報の取得は容易に行える。
【００５４】
ビード２２の品質評価に際し、６つの特徴量を個別に用いても良いし、複数または全ての特徴量を組み合わせても良い。中でも、上記したように、良い傾向を示したビード高さＨの分散値、フーリエ変換後の平均値および分散値の３つの特徴量を組み合わせることにより、良い結果がもたらされることは推察される。複数の特徴量を組み合わせる場合、個々の特徴量の重みを一定にするか、特徴量の個々に異なる重みを持たせる選択肢も考えられる。
【００５５】
【発明の効果】
この発明によれば、溶接ビードの測定区間の全体を通じて、溶接ビードの品質を正しく評価することができる。
【図面の簡単な説明】
【図１】この発明による評価方法に用いられるワークの斜視図である。
【図２】この発明による評価方法に用いられる測定装置の斜視図である。
【図３】同測定装置の光学的配置説明図である。
【図４】この発明による評価方法によって取得されたビードの３次元データを示す説明図である。
【図５】同３次元データからビードの高さデータを取得するための解説図である。
【図６】同ヒー度の幅を取得するための解説図である。
【図７】取得したビードの高さデータの一例を示すグラフである。
【図８】取得したビードの幅データの一例を示すグラフである。
【図９】ビートの中心位置の説明図である。
【図１０】ビード最大余盛幅の説明図である。
【図１１】ビードの高さデータをフーリエ変換するための説明図である。
【図１２】６つの特徴量毎の評価結果を示すグラフである。
【図１３】６つの特徴量を加算した評価結果を示すグラフである。
【符号の説明】
１１　ワーク
２１　鉄板
２２　ビード
３１　測定装置
４１　レーザ
４２　カメラ
【図１】
2004053427.tif 000007
【図２】
2004053427.tif 000008
【図３】
2004053427.tif 000009
【図４】
2004053427.tif 000010
【図５】
2004053427.tif 000011
【図６】
2004053427.tif 000012
【図７】
2004053427.tif 000013
【図８】
2004053427.tif 000014
【図９】
2004053427.tif 000015
【図１０】
2004053427.tif 000016
【図１１】
2004053427.tif 000017
【図１２】
2004053427.tif 000018
【図１３】
2004053427.tif 000019