(57)【要約】【課題】刃振れ計測を行うことなく全ての切れ刃それぞれの振れ量を推定して送り軸の制御量を算出する。【解決手段】Ｓ１で加工中の切削負荷波形と工具回転位相の信号を工具回転２周期以上取得し、Ｓ２で、得られた切削負荷波形に対し、決定した解析区間ＴＦにおいてフーリエ変換を行い、フーリエ変換結果から解析点ＴＡにおける切れ刃通過周波数ｆＣの切削負荷振幅Ａを取得する（Ｓ３）。この処理を工具回転１周期分１／Ｓ（ｓ）より大きくなるまで繰り返し実行し（Ｓ４，Ｓ５）、Ｓ６で、工具回転１周期分の各解析点ＴＡにおける切れ刃通過周波数ｆＣの切削負荷振幅Ａをプロットし、Ｓ７で、送り軸重畳制御の制御量Ｒｉおよび制御位置（工具回転位相）Ｐｉを算出する。【選択図】図２
2019214092.tif 000002
【特許請求の範囲】
【請求項１】
同心円上に複数配置される切れ刃を備えた工具を回転させながら送り軸を移動させて加工を行う工作機械において、各前記切れ刃の振れ量に基づいて算出した制御量を前記送り軸の軸送り量に重畳する制御を行うために前記制御量を算出する方法であって、
一定時間前記加工を行って加工中の切削負荷を取得する切削負荷取得ステップと、
取得した切削負荷を周波数解析して前記切れ刃に起因する前記切削負荷の振幅を取得する切れ刃負荷取得ステップと、
取得した前記切削負荷の振幅に基づいて各前記切れ刃の振れ量をそれぞれ推定して前記制御量を算出する制御量算出ステップと、
を実行することを特徴とする工作機械の送り軸の制御量算出方法。
【請求項２】
前記切れ刃負荷取得ステップでは、前記切削負荷取得ステップで取得した加工中の切削負荷波形に対し、所定の解析区間でフーリエ変換を行い、その結果から、工具回転速度と刃数との積で表わされる切れ刃通過周波数における前記切削負荷の振幅の増減の推移を工具１回転分にわたり取得し、
前記制御量算出ステップでは、前記切削負荷の振幅の増減の推移に基づいて、各前記切れ刃のそれぞれの振れ量及び工具回転位相を推定して、前記制御量とその制御位置とを算出することを特徴とする請求項１に記載の工作機械の送り軸の制御量算出方法。
【請求項３】
同心円上に複数配置される切れ刃を備えた工具を回転させながら送り軸を移動させて加工を行う工作機械において、各前記切れ刃の振れ量に基づいて算出した制御量を前記送り軸の軸送り量に重畳する送り軸の制御方法であって、
前記制御量として、請求項１又は２に記載の工作機械の送り軸の制御量算出方法で算出された制御量を用いることを特徴とする工作機械における送り軸の制御方法。
【請求項４】
同心円上に複数配置される切れ刃を備えた工具を回転させながら送り軸を移動させて加工を行うと共に、各前記切れ刃の振れ量に基づいて算出した制御量を前記送り軸の軸送り量に重畳する制御を行う工作機械であって、
一定時間前記加工を行って加工中の切削負荷を取得する切削負荷取得手段と、
取得した切削負荷を周波数解析して前記切れ刃に起因する前記切削負荷の振幅を取得する切れ刃負荷取得手段と、
取得した前記切削負荷の振幅に基づいて各前記切れ刃の振れ量をそれぞれ推定して前記制御量を算出する制御量算出手段と、
を備えることを特徴とする工作機械。
【請求項５】
同心円上に複数配置される切れ刃を備えた工具を回転させながら送り軸を移動させて加工を行うと共に、各前記切れ刃の振れ量に基づいて算出した制御量を前記送り軸の軸送り量に重畳する制御を行う工作機械のコンピュータに、請求項１又は２に記載の送り軸の制御量算出方法を実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、工作機械において工具チッピングが生じるような切削加工、たとえばチタン合金といった難削材の重切削加工に対して行われる送り軸の制御に際し、当該送り軸の制御量を算出する制御量算出方法と、その制御量算出方法を用いた送り軸の制御方法と、その制御量算出方法を実施する工作機械と、その制御量算出方法を実行させるプログラムとに関する。
【背景技術】
【０００２】
例えばミーリングで難削材を加工する場合、加工コスト低減のためにスローアウェイやインサートと呼ばれる脱着式の切れ刃を装着するタイプの工具を使用するが、工具本体の切れ刃取付け座面や切れ刃自身の加工精度の影響で、装着した刃には振れ量（各刃間の相対取付け誤差）が生じる。この値は小さくないため、刃振れ量の大きな切れ刃から工具チッピングが生じて工具寿命が短くなるといった問題があった。
この対策として、特許文献１において、工具の刃振れ量を予め測定し、加工進行方向に刃振れ量だけ送り軸を制御することでその影響を抑制するといった加工方法が提案されている。これは、各刃振れ量に応じた送り軸制御を、加工進行方向に対し、各切れ刃を取り付けた主軸の回転位相に同期させて重畳し、加工を行うものである。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第５９０８３４２号公報
【特許文献２】特開２０１６−８３７２８号公報
【特許文献３】特開２０１７−１６６２３号公報
【特許文献４】特開平９−１７４３８３号公報
【特許文献５】特開２０１０−６９５４０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上記特許文献１の加工方法において、重畳する制御量は、各切れ刃の刃振れ量と取り付け位置の情報とに基づいて、例えば上記特許文献２に記載の手法により算出できる。しかし、現状、切れ刃振れ量および取り付け位置の取得には、センサによる事前計測を行っているため、刃振れ計測に時間を要する他、センサにコストを要するという問題があった。
そこで、センサによる刃振れ計測の代替手法として、上記特許文献３に記載のように、切削負荷波形のピーク値を検出し、各切れ刃の振れ量を推定する方法が知られている。しかし、切削負荷波形は刃振れ以外にも加工中の機械振動や電気ノイズの影響を含んでいるため、正確な刃振れ量を推定するのが困難である。また、複数切れ刃が同時に切削する場合、切削負荷は複数刃の重ね合わせとなるため、切削負荷波形のピーク値と切れ刃振れ量との相関がとれず、すべての刃それぞれの振れ量が取得できないという問題がある。上記特許文献４、特許文献５には、加工中の切削負荷波形を周波数解析することで刃振れの影響を推定する方法が開示されているが、どちらも、工具すべての切れ刃のうちの最大値の振れ量の取得のみであり、前記制御量の決定に必要な、切れ刃それぞれの振れ量および取り付け位置の特定までは行っていない。
【０００５】
そこで、本発明は、刃振れ計測を行うことなく全ての切れ刃それぞれの振れ量を推定して送り軸の制御量を算出できる制御量算出方法と、その制御量算出方法を用いた送り軸の制御方法と、その制御量算出方法を実施する工作機械と、その制御量算出方法を実行させるプログラムとを提供することを目的としたものである。
【課題を解決するための手段】
【０００６】
上記目的を達成するために、請求項１に記載の発明は、同心円上に複数配置される切れ刃を備えた工具を回転させながら送り軸を移動させて加工を行う工作機械において、各前記切れ刃の振れ量に基づいて算出した制御量を前記送り軸の軸送り量に重畳する制御を行うために前記制御量を算出する方法であって、
一定時間前記加工を行って加工中の切削負荷を取得する切削負荷取得ステップと、
取得した切削負荷を周波数解析して前記切れ刃に起因する前記切削負荷の振幅を取得する切れ刃負荷取得ステップと、
取得した前記切削負荷の振幅に基づいて各前記切れ刃の振れ量をそれぞれ推定して前記制御量を算出する制御量算出ステップと、を実行することを特徴とする。
請求項２に記載の発明は、請求項１の構成において、前記切れ刃負荷取得ステップでは、前記切削負荷取得ステップで取得した加工中の切削負荷波形に対し、所定の解析区間でフーリエ変換を行い、その結果から、工具回転速度と刃数との積で表わされる切れ刃通過周波数における前記切削負荷の振幅の増減の推移を工具１回転分にわたり取得し、
前記制御量算出ステップでは、前記切削負荷の振幅の増減の推移に基づいて、各前記切れ刃のそれぞれの振れ量及び工具回転位相を推定して、前記制御量とその制御位置とを算出することを特徴とする。
上記目的を達成するために、請求項３に記載の発明は、同心円上に複数配置される切れ刃を備えた工具を回転させながら送り軸を移動させて加工を行う工作機械において、各前記切れ刃の振れ量に基づいて算出した制御量を前記送り軸の軸送り量に重畳する送り軸の制御方法であって、
前記制御量として、請求項１又は２に記載の工作機械の送り軸の制御量算出方法で算出された制御量を用いることを特徴とする。
上記目的を達成するために、請求項４に記載の発明は、同心円上に複数配置される切れ刃を備えた工具を回転させながら送り軸を移動させて加工を行うと共に、各前記切れ刃の振れ量に基づいて算出した制御量を前記送り軸の軸送り量に重畳する制御を行う工作機械であって、
一定時間前記加工を行って加工中の切削負荷を取得する切削負荷取得手段と、
取得した切削負荷を周波数解析して前記切れ刃に起因する前記切削負荷の振幅を取得する切れ刃負荷取得手段と、
取得した前記切削負荷の振幅に基づいて各前記切れ刃の振れ量をそれぞれ推定して前記制御量を算出する制御量算出手段と、を備えることを特徴とする。
上記目的を達成するために、請求項５に記載の発明は、プログラムであって、同心円上に複数配置される切れ刃を備えた工具を回転させながら送り軸を移動させて加工を行うと共に、各前記切れ刃の振れ量に基づいて算出した制御量を前記送り軸の軸送り量に重畳する制御を行う工作機械のコンピュータに、請求項１又は２に記載の送り軸の制御量算出方法を実行させることを特徴とする。
【発明の効果】
【０００７】
本発明によれば、制御量を取得するための刃振れ計測作業が省略できるため、作業時間が短縮できる他、センサレス化によって送り軸重畳制御システムを低コストで提供することができる。また、工具の切れ刃に基づく切削負荷の振幅によって各切れ刃の振れ量の推定を行うため、周波数帯の異なる加工中の機械振動や電気ノイズの影響を排除して各切れ刃の振れ量を推定することができる。よって、より適切な制御量を算出することが可能となる。
【図面の簡単な説明】
【０００８】
【図１】工作機械の構成図である。
【図２】送り軸の制御量算出方法のフローチャート図である。
【図３】ミーリング工具における加工の説明図である。
【図４】刃列数４枚の工具の刃振れ量測定結果および実加工代増減分の一例である。
【図５】図４の刃振れを持つ工具によって、加工を行った際の主軸モータトルク波形である。
【図６】図５の波形について、所定の解析データ区間でフーリエ変換を行った結果である。
【図７】図６のフーリエ変換を所定の間隔で実行し、切れ刃通過周波数における振幅の推移をプロットしたものである。
【図８】図７の結果から、横軸を工具回転位相に、縦軸を実加工代増減分に換算した結果である。
【発明を実施するための形態】
【０００９】
以下、本発明の実施の形態を図面に基づいて説明する。
図１は、本発明に係る送り軸の制御量算出方法及び制御方法を実施する工作機械の一例を示す構成図である。同図において、１はベッド、２はコラムで、コラム２の前面には、主軸頭３が、Ｘ軸制御ユニット４及びＺ軸制御ユニット５によってＸ軸方向及びＺ軸方向へ移動制御可能に設けられて、主軸頭３の下部で下向きに設けられて回転制御される主軸６に、工具７が装着されている。一方、ベッド１上には、Ｙ軸制御ユニット８によってＹ軸方向へ移動制御可能なテーブル９が設けられて、テーブル９上に被加工物１０が固定可能となっている。
【００１０】
工作機械の制御系（コンピュータ）は、主軸６の回転速度を制御する主軸回転制御装置１１、送り軸（各制御ユニット４、５、８）を制御する数値制御装置１２の他、主軸回転制御装置１１による主軸６の制御データ等を記憶する記憶部１４と、記憶部１４に記憶された加工中の切削負荷を周波数解析するための解析区間を決定する解析区間演算部１５と、決定された解析区間に対して切削負荷を周波数解析して切れ刃に基づく切削負荷の振幅を演算する周波数解析部１６と、算出された切削負荷の振幅に基づいて送り軸に対する制御量を算出する制御量演算装置１３とを備えている。ここでは主軸回転制御装置１１及び数値制御装置１２が切削負荷取得手段として、解析区間演算部１５及び周波数解析部１６が切れ刃負荷取得手段として、制御量演算装置１３が制御量算出手段として、それぞれ記憶部１４に設定されたプログラムに従って機能する。
【００１１】
このように構成された工作機械において、各制御装置１１，１２によるＮＣプログラム指令の軸送り量に対して計算した制御量を重畳させる。例えばＸ−Ｙ平面における加工であれば、加工進行方向に対して計算した制御量をＸ軸、Ｙ軸方向に分配して送り軸を制御する。この制御を行うためには工具切れ刃の各位置とその刃振れ量に応じた、最適な制御量と制御位置（工具回転位相）を決定しておく必要がある。従来の制御量取得方法においては、工具７の切れ刃毎の振れ量をセンサにより測定し、制御量演算装置１３において、所定計算式を用いて工具回転位相における制御量を算出している。
しかし、刃振れ計測には時間を要する他、センサにコストを要するという問題があった。
【００１２】
そこで、本形態では、刃振れ計測センサを要せず、一定時間取得した加工中の切削負荷を周波数解析し、工具回転速度と刃列数から算出される切れ刃通過周波数の振幅増減分から切れ刃の振れ量を推定し、制御量を取得する。以下、この制御量決定の流れを図２のフローチャートに従って説明する。
まず、加工中の切削負荷波形と工具回転位相の信号を工具回転２周期以上取得し、記憶部１４に記憶する（Ｓ１：切削負荷取得ステップ）。このとき取得する切削負荷には、主軸回転制御装置１１から取得可能な主軸モータトルクを用いてもよい。また、工具回転位相の信号には主軸６に取り付けられたエンコーダから発生する電圧信号が利用できる。
続いて、得られた切削負荷波形に対して周波数解析を行う（Ｓ２）。解析区間演算部１５において、解析点ＴＡと、解析点ＴＡを中心にとる解析データ区間ＴＦを算出し、得られた解析区間に対して周波数解析部１６によりフーリエ変換を行う。このときの解析データ区間ＴＦは、工具１回転以内の区間が好適であり、例えば、外部入力装置１７により入力された径方向切込量ａｅ（ｍｍ）と工具半径Ｒ（ｍｍ）をもとに、以下の式１で得られる切れ刃一刃が切削にかかる時間を用いるとよい。
ＴＦ＝ｃｏｓ−１（１−ａｅ／Ｒ）／（３６０×Ｓ）    （ｓ） ・・式１
【００１３】
この解析区間に関してフーリエ変換を行い、切れ刃通過周波数の振幅を取得することで、切れ刃ごとの切削負荷が取得可能となる。解析データに関しては、フーリエ変換を行う前に、窓関数処理や、データ末尾へのゼロ付加を行っておき、フーリエ変換時の漏れ誤差低減や周波数分解能を向上するようにしておくのが望ましい。
解析点ＴＡには初期値ＴＳ＝ＴＦ／２を用い、以降は、後述の繰り返し処理時（Ｓ４）に時間ΔＴ（ｓ）ずつ大きくし（Ｓ５）、解析区間をオーバーラップさせて繰り返し解析を実行する。このとき、各切れ刃ごとの切削負荷を取得する必要があるため、ΔＴは後述する切れ刃通過周波数ｆＣの周期１／ｆＣの半分以下とする。
【００１４】
フーリエ変換結果から解析点ＴＡにおける切れ刃通過周波数ｆＣの切削負荷振幅Ａを取得する（Ｓ３）。切れ刃通過周波数ｆＣは、工具回転速度Ｓ（ｒｅｖ／ｓ）と刃列数Ｚとの積で表わされ、この周波数に注目することで、加工中の機械振動や切削負荷波形取得時の電気ノイズ等の影響を排除することができる。
以上の処理を、解析点ＴＡをΔＴだけ変化させ、初期値ＴＳから工具回転１周期分１／Ｓ（ｓ）より大きくなるまで繰り返し実行し（Ｓ４，Ｓ５）、Ｓ４で工具回転１周期分より大きくなったことを判別したら、工具回転１周期分の各解析点ＴＡにおける切れ刃通過周波数ｆＣの切削負荷振幅Ａをプロットする（Ｓ６）。これが、工具１回転中における刃振れによる切削負荷の増減である。ここまでのＳ２〜Ｓ６が切れ刃負荷取得ステップとなる。
【００１５】
以上の結果を元に、送り軸重畳制御の制御量Ｒｉおよび制御位置（工具回転位相）Ｐｉを制御量演算装置１３により算出する（Ｓ７：制御量算出ステップ）。
まず、切削負荷振幅Ａを実加工代増減分Ｄに換算する。実加工代増減分Ｄは、切れ刃通過周波数ｆＣにおける切削負荷振幅Ａと、その工具回転１周期平均値ＡＡＶＥ、外部入力装置１７から得られる一刃送りＦＺ（μｍ／刃）から、以下の式２によって算出する。
Ｄ＝Ａ／ＡＡＶＥ×ＦＺ   （μｍ）  ・・式２
次に、主軸６に接続されているエンコーダから各解析点における工具回転位相ＰＡを取得する。主軸モータが所定の位相ＰＥになった際、エンコーダから主軸回転制御装置１１に電圧信号が発生するため、位相ＰＥを基準として各解析点ＴＡにおける工具回転位相ＰＡを求めることができる。得られた工具回転位相ＰＡと実加工代増減分Ｄの結果から特徴点を、上記重畳制御における制御回数に応じて抽出し、制御位置Ｐｉと制御量Ｒｉを決定する。特徴点には、例えば、実加工代増減分Ｄが最大値となる１点をとり、そこから上記重畳制御の回数で等分した位相間隔ごとに点を抽出し、抽出した各点における工具回転位相を制御位置Ｐｉ、実加工代増減分Ｄを制御量Ｒｉとする。
【００１６】
以上より、上記重畳制御における制御位置Ｐｉおよび制御量Ｒｉが決定する。また、上記図２の本実施例フローチャートＳ１〜Ｓ７までをさらに繰り返し、実加工代増減分Ｄの変動幅が、より小さく収束するまで、すなわち、全ての切れ刃における切削負荷がより均一となるまで実行し、制御量をより好適な値へ収束させてもよい。決定した制御量は、数値制御装置１２に転送され、送り軸の重畳制御に用いることができる。
【００１７】
上記計算手法による手順の具体的な実施例を、図３の説明図に示すミーリング工具における加工例をもとに説明する。
まず、図３の加工は、４枚の切れ刃７(1)、７(2)、７(3)、７(4)が、９０°間隔で同心円上に配置された工具７により、被加工物１０を＋Ｘ方向に工具回転速度Ｓ＝５（ｒｅｖ／ｓ）、一刃送りＦＺ＝１００（μｍ／刃）で溝切削している。今回の実施例は溝切削であるため、同時に複数の切れ刃（図３における７(1)、７(2)）が切削に寄与することとなる。このとき使用している工具の各切れ刃取り付け位相及び刃振れ量、実加工代増減分は図４に示すとおりである。実加工代増減分は刃振れ量による加工代の差分であり、刃列の番号を添え字ｉ（１≦ｉ≦Ｚ）、測定した各切れ刃の振れ量をＣｉ（μｍ）とすると、各切れ刃における実加工代増減分Ｄｉの実測結果は、以下の式３により算出することができる。
Ｄｉ＝Ｃｉ−Ｃｉ−１    （但し、ｉ−１＝０の場合はＺに置き換える）  ・・式３
【００１８】
以降、上述の図２のフローチャートの手順に従い、実際に上記制御量を決定する。
加工中の切削負荷波形として、切削中の主軸モータトルク波形を工具回転２周期分取得した結果を図５に示す（Ｓ１）。上述の通り、今回の加工条件においては、同時に複数の切れ刃が切削に寄与するため、モータトルク波形は複数の切れ刃による負荷の重ね合わせとなり、各切れ刃の刃振れ量は推定できない。
続いて、得られた波形に対し周波数解析を行う（Ｓ２）。溝切削であるため、径方向切込量ａｅ＝２Ｒ（ｍｍ）と工具半径Ｒ（ｍｍ）から、上述の式２により得られる解析区間ＴＦ＝０．１（ｓ）、解析点ＴＡの初期値ＴＳ＝０．０５（ｓ）としてフーリエ変換を行った結果が図６である。このとき、解析データに対しては、フーリエ変換前にそれぞれ、漏れ誤差低減と周波数分解能向上のため、ハニング窓関数処理とデータ長２（ｓ）のゼロ付加を行っている。
【００１９】
得られたフーリエ変換結果から、加工条件である、工具回転速度Ｓ＝５（ｒｅｖ／ｓ）、刃列数Ｚ＝４の工具切れ刃通過周波数ｆＣ＝Ｓ×Ｚ＝２０（Ｈｚ）における、切削負荷振幅Ａを取得した（Ｓ３）。これを、工具回転１周期分０．２（ｓ）間のΔＴ＝０．０１（ｓ）で繰り返し実行し（Ｓ４，Ｓ５）、切削負荷振幅Ａの推移をプロットした結果が図７である（Ｓ６）。
図７の結果を、工具回転位相ＰＡと実加工代増減分Ｄに換算し、上記重畳制御における制御位置Ｐｉおよび制御量Ｒｉを決定する（Ｓ７）。実施例では、解析点ＴＳ＝０．０５（ｓ）において、上述の主軸エンコーダから電圧信号が発生し、このときの工具回転位相ＰＥ＝０（°）である。これより、解析点における工具回転位相ＰＡが分かる。
続いて、加工条件の一刃送りＦＺ＝１００（μｍ／刃）から、式３により、実加工代増減分Ｄが得られる。
【００２０】
以上の換算結果を図８に示す。図中(1)〜(4)で示される点はそれぞれ図４において示した切れ刃(1)〜(4)の取り付け位相となり、このとき、図８の実加工代増減分Ｄの換算値と図４の切れ刃の実加工代増減分の事前計測結果とがおおよそ一致しているのが分かる。すなわち、実加工代増減分Ｄを切れ刃それぞれの振れ量の推定値と見なすことができる。したがって、図８の結果が、送り軸重畳制御における制御位置および制御量となる。実加工代増減分Ｄが最大となる工具回転位相ＰＡ＝０（°）を基準とし、制御回数４回分９０（°）間隔で等分した点の位相を制御位置Ｐｉ、加工代増減分を制御量Ｒｉと決定する。
【００２１】
このように、上記形態の送り軸の制御量算出方法、送り軸の制御方法、工作機械、プログラムによれば、切削負荷を周波数解析して切れ刃による切削負荷の振幅を求め、その振幅に基づいて切れ刃の振れ量を推定して制御量を算出するので、刃振れ計測を行うことなく全ての切れ刃それぞれの振れ量を推定して送り軸の制御量を算出可能となる。よって、作業時間が短縮できる他、センサレス化によって送り軸重畳制御システムを低コストで提供することができる。また、工具の切れ刃に基づく切削負荷の振幅によって各切れ刃の振れ量の推定を行うため、周波数帯の異なる加工中の機械振動や電気ノイズの影響を排除して各切れ刃の振れ量を推定することができる。よって、より適切な制御量を算出することが可能となる。
【００２２】
なお、切れ刃の数や位置、形状は上記形態に限らず、適宜変更可能である。切れ刃が工具の軸方向に複数段あってもよい。工作機械も、複合加工機やマシニングセンタ等、特に機種を限定するものではない。推定する切れ刃の振れ量も、実加工代増減分以外の物理量とすることもできる。
また、本発明の制御量算出方法は、工作機械自身のコンピュータにより実行する必要はなく、工作機械と無線或いは有線接続された別のコンピュータに本発明のプログラムを格納して実行させたりしてもよい。
【符号の説明】
【００２３】
１・・ベッド、２・・コラム、３・・主軸ヘッド、４・・Ｘ軸制御ユニット、５・・Ｚ軸制御ユニット、６・・主軸、７・・工具、７(1)〜７(4)・・切れ刃、８・・Ｙ軸制御ユニット、９・・テーブル、１０・・被加工物、１１・・主軸回転制御装置、１２・・数値制御装置、１３・・制御量演算装置、１４・・記憶部、１５・・解析区間演算部、１６・・周波数解析部、１７・・外部入力装置。
【図１】
2019214092.tif 000003
【図２】
2019214092.tif 000004
【図３】
2019214092.tif 000005
【図４】
2019214092.tif 000006
【図５】
2019214092.tif 000007
【図６】
2019214092.tif 000008
【図７】
2019214092.tif 000009
【図８】
2019214092.tif 000010