(57)【要約】  繊維強化複合材料である試験体の引張強度を評価するための強度検査装置は、ＡＥセンサ１と対象波特定部２と演算部３と評価データ生成部４を備える。ＡＥセンサ１は、引張荷重を増加させながら試験体に与えている試験期間において、引張荷重により試験体において生じるＡＥ波を検出してＡＥ波の波形データを生成する。対象波特定部２は、生成された波形データに基づいて、持続時間が時間閾値よりも長いＡＥ波を対象波として特定する。演算部３は、各対象波の重心周波数を求める。評価データ生成部４は、各対象波について、対象波の重心周波数と、対象波の検出時点で試験体に与え引張荷重の大きさとを対応付けた強度評価データを生成する。
2020170360.tif 000002
【特許請求の範囲】
【請求項１】
繊維強化複合材料である試験体の引張強度を評価するための強度検査装置であって、
引張荷重を増加させながら前記試験体に与えている試験期間において、該引張荷重により前記試験体において生じるＡＥ波を検出してＡＥ波の波形データを生成するＡＥセンサと、
前記波形データに基づいて、持続時間が時間閾値よりも長いＡＥ波を対象波として特定する対象波特定部と、
各対象波の重心周波数を求める演算部と、
各対象波について、該対象波の前記重心周波数と、該対象波の検出時点で前記試験体に与えた前記引張荷重の大きさとを対応付けた強度評価データを生成する評価データ生成部と、を備える強度検査装置。
【請求項２】
波形データにおいて、設定時間が経過する度に、当該設定時間中に大きさが少なくとも1回は変位閾値以上になっている一連の変位を１つのＡＥ波として、対象波特定部は、前記時間閾値よりも長く持続しているＡＥ波を対象波として特定する、請求項１に記載の強度検査装置。
【請求項３】
前記対象波特定部は、前記波形データに基づいて、持続時間が前記時間閾値よりも長いＡＥ波を対象波として特定し、各対象波と、当該対象波の検出時点とを対応付けた対象波データを生成し、
前記演算部は、前記対象波データに基づいて、各対象波の重心周波数を求め、各対象波の重心周波数と、該対象波の検出時点とを対応付けた重心周波数データを生成し、
前記評価データ生成部は、前記重心周波数データと、経過時間に対する前記引張荷重の大きさを表わす荷重データとに基づいて、前記強度評価データを生成する、請求項１に記載の強度検査装置。
【請求項４】
前記演算部は、
各対象波の波形に基づいて該対象波のスペクトルデータを生成するスペクトル生成部と、
各対象波の前記スペクトルデータに基づいて、該ＡＥ波の重心周波数を求める重心周波数算出部と、を有する、請求項１に記載の強度検査装置。
【請求項５】
前記強度評価データを表示するディスプレイ装置を備える、請求項１に記載の強度検査装置。
【請求項６】
繊維強化複合材料である試験体の引張強度を評価するための強度検査方法であって、
引張荷重を増加させながら前記試験体に与え、これにより前記試験体において生じるＡＥ波の波形データを、ＡＥセンサにより生成し、
前記波形データに基づいて、持続時間が時間閾値よりも長いＡＥ波を、対象波特定部により対象波として特定し、
各対象波の重心周波数を演算部により求め、
各対象波について、該対象波の前記重心周波数と、該対象波の検出時点で前記試験体に与えた前記引張荷重の大きさとを対応付けた強度評価データを、評価データ生成部により生成する、強度検査方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本開示は、繊維強化複合材料（ＦＲＰ：Ｆｉｂｅｒ  Ｒｅｉｎｆｏｒｃｅｄ  Ｐｌａｓｔｉｃ）である試験体の引張強度を検査する技術に関する。より詳しくは、本開示は、試験体に引張荷重を与えることにより、試験体に生じたＡＥ波（ａｃｏｕｓｔｉｃ  ｅｍｉｓｓｉｏｎ）に基づいて、試験体の引張強度を検査する技術に関する。
【背景技術】
【０００２】
ＦＲＰは、ロケットや航空機などに用いられている。特に、炭素繊維強化複合材料（ＣＦＲＰ：Ｃａｒｂｏｎ  Ｆｉｂｅｒ  Ｒｅｉｎｆｏｒｃｅｄ  Ｐｌａｓｔｉｃ）は、強度と剛性に優れている。ＦＲＰでは、積層の剥離や繊維の断線が生じた後に破壊に至る。【０００３】
このようなＦＲＰの引張強度を検査する技術が、特許文献１に開示されている。特許文献１の強度検査方法では、次のようにＦＲＰの引張強度を評価している。試験体に引張荷重を与え、この引張荷重を、時間の経過に従って増やす。この引張荷重により試験体に生じるＡＥ波を検出する。試験期間に含まれる複数の荷重印加区間の各々におけるＡＥ波の複数の周波数成分を求める。各荷重印加区間について、複数の周波数成分に基づいてＡＥ波の重心周波数を求め、複数の荷重印加区間のうち、先の荷重印加区間よりも重心周波数が下がっている荷重印加区間を特定し、特定した荷重印加区間において試験体に与えた引張荷重の大きさを、試験体の引張強度と判定する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第５８４１０８１号
【発明の概要】
【発明が解決しようとする課題】
【０００５】
試験体に引張荷重を与えている時に、多数のＡＥ波（一連の変位）が断続的に生じるが、これらのＡＥ波の中には、試験体の破断の有無を表わさない傾向にあるＡＥ波も含まれる。そのため、試験体の破断の有無を表わす傾向にあるＡＥ波を特定して、特定したＡＥ波に基づいて試験体の引張強度を評価することが望まれる。【０００６】
すなわち、本開示は、試験体の破断の有無を表わす傾向にあるＡＥ波を特定し、特定したＡＥ波に基づいて、試験体の引張強度を評価する技術を提供することにある。
【課題を解決するための手段】
【０００７】
本開示による強度検査装置は、繊維強化複合材料である試験体の引張強度を評価するための装置であって、
引張荷重を増加させながら前記試験体に与えている試験期間において、該引張荷重により前記試験体において生じるＡＥ波を検出してＡＥ波の波形データを生成するＡＥセンサと、
前記波形データに基づいて、持続時間が時間閾値よりも長いＡＥ波を対象波として特定する対象波特定部と、
各対象波の重心周波数を求める演算部と、
各対象波について、該対象波の前記重心周波数と、該対象波の検出時点で前記試験体に与えた前記引張荷重の大きさとを対応付けた強度評価データを生成する評価データ生成部と、を備える。【０００８】
本開示による強度検査方法は、繊維強化複合材料である試験体の引張強度を評価するための方法であって、
引張荷重を増加させながら前記試験体に与え、これにより前記試験体において生じるＡＥ波の波形データを、ＡＥセンサにより生成し、
前記波形データに基づいて、持続時間が時間閾値よりも長いＡＥ波を、対象波特定部により対象波として特定し、
各対象波の重心周波数を演算部により求め、
各対象波について、該対象波の前記重心周波数と、該対象波の検出時点で前記試験体に与えた前記引張荷重の大きさとを対応付けた強度評価データを、評価データ生成部により生成する、強度検査方法。
【発明の効果】
【０００９】
本開示によると、引張荷重を増加させながら試験体に与えている試験期間において、試験体において生じるＡＥ波の波形データを生成し、この波形データに基づいて、持続時間が時間閾値よりも長いＡＥ波を対象波として特定し、各対象波について、その重心周波数と引張荷重の大きさとを対応付けた強度評価データを生成する。【００１０】
これについて、持続時間が時間閾値よりも長いＡＥ波は、試験体の破断の有無を表わす周波数情報を含んでいる傾向にある。したがって、試験体の破断の有無を表わす傾向にあるＡＥ波を用いて生成した強度評価データに基づいて、試験体の引張強度を評価することができる。
【図面の簡単な説明】
【００１１】
【図１】本開示の実施形態による強度検査装置を示すブロック図である。
【図２】ＡＥセンサが生成する波形データの一例を示す模式図である。
【図３Ａ】重心周波数データの一例を示す概念図である。
【図３Ｂ】荷重データの一例を示す概念図である。
【図３Ｃ】強度評価データの一例を示す概念図である。
【図４】本開示の実施形態による強度検査方法を示すフローチャートである。
【図５Ａ】比較例の強度検査方法により得られた強度評価データを示す。
【図５Ｂ】図５Ａにおける重心周波数のプロットの密度をグレースケールで表わした図である。
【図６Ａ】本開示の実施形態による強度検査方法により得られた強度評価データを示し、時間閾値が３０マイクロ秒である場合を示す。
【図６Ｂ】本開示の実施形態による強度検査方法により得られた強度評価データを示し、時間閾値が２００マイクロ秒である場合を示す。
【発明を実施するための形態】
【００１２】
本開示の実施形態を図面に基づいて説明する。なお、各図において共通する部分には同一の符号を付し、重複した説明を省略する。
【００１３】
（強度検査装置の構成）
図１は、本開示の実施形態による強度検査装置１０を示すブロック図である。強度検査装置１０は、繊維強化複合材料である試験体に引張荷重を与えた時に試験体に生じるＡＥ波に基づいて、試験体の引張強度を評価するための装置である。ここで、繊維強化複合材料は、例えばＣＦＲＰである。試験体は、例えば、ロケット又は航空機を構成するＣＦＲＰ、あるいは、水素を燃料として走行する自動車に積載され水素を蓄える水素タンクを構成するＣＦＲＰであってよいが、これらに限定されない。
【００１４】
強度検査装置１０は、ＡＥセンサ１と対象波特定部２と演算部３と評価データ生成部４とを備える。
【００１５】
ＡＥセンサ１は、試験体に取り付けられ、試験体に生じたＡＥ波を検出する。より詳しくは、ＡＥセンサ１は、引張荷重を増加させながら試験体に与えている試験期間において、該引張荷重により試験体において生じるＡＥ波の波形データを生成する。波形データは、各時点におけるＡＥ波の変位を表わす。
【００１６】
図２は、波形データの一例を示す模式図である。図２において、横軸は経過時間を示し、縦軸は、ＡＥ波の変位を示す。図２では、波形データにおける一連の変位を１つのＡＥ波として２つのＡＥ波Ｗ１，Ｗ２を示している。なお、試験期間においては、多数のＡＥ波が断続的に生じる。ＡＥセンサ１は、後述する時間計測部６が計測した経過時間に基づいてＡＥ波の波形データを生成してよい。この時間計測部６は、ＡＥセンサ１に組み込まれていてもよい。
【００１７】
対象波特定部２は、ＡＥセンサ１が生成した波形データに基づいて、持続時間が時間閾値よりも長いＡＥ波を対象波として特定する。より詳しくは、波形データにおいて、設定時間（図２ではＴｓ）が経過する度に、当該設定時間Ｔｓ中に大きさが少なくとも1回は変位閾値（図２ではＤｔ）以上になっている一連の変位を１つのＡＥ波として、対象波特定部２は、時間閾値よりも長く持続しているＡＥ波を対象波として特定する。すなわち、対象波特定部２は、各ＡＥ波を特定し、当該各ＡＥ波の持続時間を求め、これらのＡＥ波から、持続時間が時間閾値よりも長い各ＡＥ波を対象波として特定する。例えば、図２において、ＡＥ波Ｗ１は、持続時間Ｔ１が時間閾値以下であるので対象波として特定されず、ＡＥ波Ｗ１は、持続時間Ｔ２が時間閾値よりも長いので対象波として特定される。
【００１８】
一例では、対象波特定部２は、次のように処理を行う。対象波特定部２は、上述の波形データから各ＡＥ波の持続時間の開始時点を特定する。この開始時点は、図２のように、ＡＥ波の変位の大きさが閾値Ｄｔ未満のまま、設定時間Ｔｓを経過した後に、変位の大きさが変位閾値Ｄｔに最初に達した時点である。また、対象波特定部２は、図２のように、持続時間の開始時点の後、ＡＥ波の変位の大きさが変位閾値Ｄｔ未満のまま最初に設定時間Ｔｓを経過した時点を特定し、当該設定時間Ｔｓの開始時点を持続時間の終了時点として特定する。対象波特定部２は、このように各ＡＥ波の持続時間を求め、各ＡＥ波について、その持続時間が時間閾値よりも長いかを判断する。これにより、対象波特定部２は、多数のＡＥ波のうち、持続時間が時間閾値よりも長い各ＡＥ波を対象波として特定する。
【００１９】
なお、適切な時間閾値は、試験体としての繊維強化複合材料の種類により変わるので、試験体の種類毎に実験的に予め求められてよい。また、一例では、試験期間において、持続時間が時間閾値以下のＡＥ波が、持続時間が時間閾値よりも長いＡＥ波（対象波）よりも多く生じるように、時間閾値が設定されてよい。この設定は、試験体と同じ構造の繊維強化複合材料に対する実験に基づいてなされてよい。
【００２０】
対象波特定部２は、例えば、ＡＥセンサ１により生成された波形データに基づいて、上述のように対象波を特定し、各対象波（すなわち対象波の波形）と、当該対象波の検出時点とを対応付けた対象波データを生成する。なお、対象波の波形は、ＡＥセンサ１により生成された波形データに含まれており、経過時間に対する対象波の変位の大きさを表わす波形である（以下同様）。また、対象波の検出時点は、例えば、対象波の持続時間の開始時点、終了時点、又は中間時点であってよい。
【００２１】
演算部３は、特定された各対象波の重心周波数を求める。演算部３は、例えば、上述の対象波データに基づいて、各対象波の重心周波数を求め、各対象波の重心周波数と、該対象波の検出時点とを対応付けた重心周波数データを生成する。図３Ａは、重心周波数データの一例を示す概念図である。図３Ａにおいて、横軸は、試験期間における経過時間を示し、縦軸は、各対象波について算出した重心周波数を示し、各白丸は、それぞれ、対応する対象波について算出された重心周波数を示す。演算部３は、スペクトル生成部３ａと重心周波数算出部３ｂとを有する。
【００２２】
スペクトル生成部３ａは、対象波データに基づいて、特定された各対象波の波形に基づいて、当該対象波のスペクトルデータを生成する。より詳しくは、スペクトル生成部３ａは、当該ＡＥ波の持続時間にわたって、経過時間に対して当該ＡＥ波の変位（強度）を表わした波形を、周波数に対してＡＥ波の強度を表わしたスペクトルデータに変換する。この変換は、ＦＦＴ（Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）により行われてよい。
【００２３】
重心周波数算出部３ｂは、各対象波について、当該対象波のスペクトルデータに基づいて、当該対象波の重心周波数を求める。重心周波数Ｆｇは、次の式により表される。
Ｆｇ＝Σ（Ｆi×Ｐｉ）／ΣＰi
ここで、Ｆiは、スペクトルデータにおける各周波数を示し、Ｐｉは、スペクトルデータにおける対象波の周波数成分（すなわち、周波数ＦiでのＡＥ波の強度）を示し、ＦｉとＰｉの添え字ｉは、複数の周波数を互いに区別するための指標値であって、１〜ｎ（ｎは、２以上の整数であり、好ましくは、十分に大きい値）までの値をとり、Σは、ｉのすべての値についての総和を示す。
【００２４】
なお、重心周波数算出部３ｂは、各対象波について、当該対象波のスペクトルデータの全ての周波数成分に基づいて当該対象波の重心周波数を求めてもよいし、当該対象波のスペクトルデータにおける所定の周波数域に含まれる各周波数成分のみに基づいて当該対象波の重心周波数を求めてもよい。ここで、所定の周波数域は、ＡＥセンサ１の共振周波数を除く（例えば当該共振周波数よりも低い）周波数域であってよい。
【００２５】
評価データ生成部４は、各対象波について、当該対象波の重心周波数と、当該対象波の検出時点で試験体に与えた引張荷重の大きさとを対応付けた強度評価データを生成する。評価データ生成部４は、例えば、重心周波数データと、試験期間における経過時間に対する引張荷重の大きさを表わす荷重データとに基づいて、強度評価データを生成する。そのために、評価データ生成部４には、重心周波数データ以外に荷重データが入力される。
【００２６】
図３Ｂは、荷重データの一例を示す概念図である。図３Ｂにおいて、横軸は、荷重増加時間を示し、縦軸は、試験期間において試験体に与えた引張荷重を示す。荷重データは、荷重データ生成部５により生成されてよい。荷重データ生成部５は、試験期間において試験体に与えた引張荷重の計測値と、試験期間において時間計測部６が計測した経過時間とに基づいて、荷重データを生成し、この荷重データを評価データ生成部４に入力する。ここで、荷重の計測値は、適宜の荷重センサ７（例えば試験体に取り付けたひずみゲージ）により計測された値であってよい。なお、荷重データ生成部５は、強度検査装置１０の構成要素であってもよい。また、時間計測部６と荷重センサ７も、強度検査装置１０の構成要素であってもよい。
【００２７】
図３Ｃは、強度評価データの一例を示す概念図である。図３Ｃの横軸は、試験体に与えた引張荷重を示し、経過時間に対応する。図３Ｃの縦軸は、重心周波数を示す。また、図３Ｃにおいて、多数の白丸は、それぞれ、対応する引張荷重について算出された重心周波数を示す。例えば、図３Ｂの荷重データによる引張荷重が試験体に与えられた結果、図３Ａの重心周波数データと図３Ｃの強度評価データが得られる。
【００２８】
評価データ生成部４は、例えば、強度評価データをディスプレイ装置８に出力する。この場合には、ディスプレイ装置８は、その画面に強度評価データを表示する。強度評価データは、図３Ｃのように、引張荷重の増加に伴う、重心周波数の変化を表わす。このような強度評価データにおいて、引張荷重の増加により重心周波数が（例えば顕著に）低下する直前の引張荷重の大きさを、試験体の引張強度として判定することができる。すなわち、繊維強化複合材料は、引張荷重により、その繊維が破断して破壊されると、この材料には、高い周波数のＡＥ波が伝わり難くなると考えられる。したがって、引張荷重の増加過程で試験体が破壊されると、以降において引張荷重により発生するＡＥ波の重心周波数が下がる。したがって、重心周波数が低下する直前の引張荷重の大きさを試験体の引張強度として判定できる。例えば、図３Ｃにおいて破線で示す引張荷重を引張強度と判定できる。
【００２９】
なお、評価データ生成部４は、強度評価データを他の装置（例えば、プリンタ装置または記憶装置）に出力してもよい。この場合、プリンタ装置は、強度評価データを用紙に印刷し、当該記憶装置は、強度評価データを記憶する。
【００３０】
（強度検査方法）
図４は、本開示の実施形態による強度検査方法を示すフローチャートである。強度検査方法は、繊維強化複合材料である試験体の引張強度を検査するために、以下のステップＳ１〜Ｓ５を有する。強度検査方法は、上述した強度検査装置１０を用いて行われる。
【００３１】
ステップＳ１では、試験期間において引張荷重を増加させながら試験体に与え、これにより試験体において生じるＡＥ波の波形データを、ＡＥセンサ１により生成する。
【００３２】
ステップＳ２では、ステップＳ１で生成された波形データに基づいて、持続時間が時間閾値よりも長いＡＥ波を、対象波特定部２により対象波として特定する。例えば、対象波特定部２は、波形データに基づいて上述の対象波データを生成する。
【００３３】
ステップＳ３では、ステップＳ２で特定された各対象波の重心周波数を演算部３により求める。例えば、演算部３は、ステップＳ２で生成された対象波データに基づいて、上述の重心周波数データを生成する。
【００３４】
ステップＳ３は、ステップＳ３１，Ｓ３２を有する。ステップＳ３１では、ステップＳ２で生成された対象波データにおける各対象波の波形に基づいて、当該対象波のスペクトルデータを、スペクトル生成部３ａにより生成する。ステップＳ３２では、ステップＳ３１で生成した各対象波のスペクトルデータに基づいて、当該ＡＥ波の重心周波数を、重心周波数算出部３ｂにより求める。
【００３５】
ステップＳ４では、各対象波について、当該対象波の重心周波数と、当該対象波の検出時点で試験体に与えた引張荷重の大きさとを対応付けた強度評価データを、評価データ生成部４により生成する。例えば、評価データ生成部４は、ステップＳ３で生成された重心周波数データと、ステップＳ１における経過時間に対する引張荷重の大きさを表わす荷重データとに基づいて、強度評価データを生成する。
【００３６】
ステップＳ５では、ステップＳ４で生成された強度評価データが評価データ生成部４により出力される。例えば、評価データ生成部４は、強度評価データをディスプレイ装置８に出力する。この場合、ディスプレイ装置８は、その画面に、図３Ｃのような強度評価データを表示する。人が、表示された強度評価データを見て、引張荷重の増加過程において重心周波数の低下（例えば顕著に）が始まる直前の引張荷重の大きさを、試験体の引張強度として判定することができる。
【００３７】
（実施例）
図５Ａは、比較例の強度検査方法により実際に得られた強度評価データを示す。図５Ａの場合では、上述のステップＳ２において、ステップＳ１で生成した波形データにおいて、大きさが変位閾値以上となる一連の変位を１つのＡＥ波として、ＡＥ波の持続時間に係わらず、全てのＡＥ波の各々について重心周波数を求め、他の点は、上述した本開示の実施形態による強度検査方法と同じである。図５Ａにおいて、横軸は、ステップＳ１で与えた引張荷重を示し、縦軸は、重心周波数を示す。図５Ａにおいて、小さい各白抜きの丸印は、１つのＡＥ波の重心周波数のプロットである。
【００３８】
図５Ａの強度評価データでは、試験体の破断を表わしていない各白抜きの丸印が多いので、引張荷重がどの値を超えたら重心周波数が低下していると言えるのかを判定し難い。図５Ｂは、図５Ａにおいて、重心周波数のプロット（多数の白抜きの丸印）の密度をグレースケールで表わした図である。図５Ｂにおいて、濃い色の領域が、白抜きの丸印の密度が高い領域である。この領域から、図５Ｂの矢印が示すように、引張荷重が概ね４７５ＭＰａを超えたら重心周波数が低下していると言える。したがって、試験体の引張強度は約４７５ＭＰａであることが分かる。
【００３９】
図６Ａと図６Ｂは、本開示の実施形態による強度検査方法により実際に得られた強度評価データを示す。図６Ａと図６Ｂは、図５Ａの場合と同じ、ステップＳ１で生成した波形データから得られたデータである。図６Ａは、持続時間が３０マイクロ秒を超えるＡＥ波を対象波として特定した場合（すなわち、上記時間閾値を３０マイクロ秒とした場合）を示す。図６Ｂは、持続時間が２００マイクロ秒を超えるＡＥ波を対象波として特定した場合を示す。
【００４０】
図６Ａと図６Ｂの強度評価データにおいて、これらの図の矢印が示すように、引張荷重が概ね４７５ＭＰａを超えたら重心周波数が低下していることが分かる。持続時間がより長い各対象波から生成した図６Ｂの強度評価データでは、引張荷重がどの値を超えたら重心周波数が低下しているのかを、図６Ａの場合よりも精度よく判定することができる。
【００４１】
（実施形態による効果）
持続時間が時間閾値以下であるＡＥ波には、試験体の破断とは関係しない比較的高い周波数成分が含まれている傾向がある。したがって、このような持続時間の短いＡＥ波の重心周波数は、試験体の破断の有無を表わさない傾向にある。
これに対し、本開示の実施形態では、上述のように、このような持続時間が短いＡＥ波を用いずに、重心周波数が試験体の破断の有無を表わす傾向にある持続時間の長いＡＥ波（対象波）を用いて強度評価データを生成している。したがって、このような強度評価データに基づいて、試験体の引張強度を精度よく判定することが可能となる。
【００４２】
また、持続時間が時間閾値以下のＡＥ波に対する処理（例えば重心周波数を求める処理）が不要となるので、強度評価データを生成するための処理量が低減される。例えば、図５Ａのように膨大な数のＡＥ波の重心周波数を生成することが不要になり、図５Ａのデータから図５Ｂのデータを生成する処理が不要になる。
【００４３】
本開示は上述した実施の形態に限定されず、本開示の要旨を逸脱しない範囲で種々変更を加え得ることは勿論である。
【符号の説明】
【００４４】
１  ＡＥセンサ、２  対象波特定部、３  演算部、３ａ  スペクトル生成部、３ｂ  重心周波数算出部、４  評価データ生成部、５  荷重データ生成部、６  時間計測部、７  荷重センサ、８  ディスプレイ装置、１０  強度検査装置、Ｔｓ  設定時間、Ｄｔ  変位閾値
【図１】
2020170360.tif 000003
【図２】
2020170360.tif 000004
【図３Ａ】
2020170360.tif 000005
【図３Ｂ】
2020170360.tif 000006
【図３Ｃ】
2020170360.tif 000007
【図４】
2020170360.tif 000008
【図５Ａ】
2020170360.tif 000009
【図５Ｂ】
2020170360.tif 000010
【図６Ａ】
2020170360.tif 000011
【図６Ｂ】
2020170360.tif 000012
【国際調査報告】
2020170360.tif 000013
2020170360.tif 000014
2020170360.tif 000015
2020170360.tif 000016