(57)【要約】【課題】停止又は惰性回転している電動機の速度センサレストルク制御を実施する場合、一次抵抗誤差により、トルク制御開始時の初期速度及び初期二次磁束推定値に誤差を生じる。その結果、トルク制御手段１の初期値に誤差があるために、電動機６のトルク制御を高精度に行うことができず、電動機６のトルク制御が不能となることもある。【解決手段】同定用タイマ、磁束演算器、回転座標変換器、回転磁束メモリ、極値時間推定器、極値磁束抽出器、R1B同定器で構成される一連のブロックから得られる一次抵抗同定量を用いて、トルク制御開始時の初期速度及び初期二次磁束推定値を演算する。【選択図】図１
2006014488.tif 000002
【特許請求の範囲】
【請求項１】
電動機に流れる電流ベクトルを検出する電流検出器と、前記電流ベクトルを入力し前記電動機に直流電流を流す拾い上げ電圧指令を出力する拾い上げ制御手段と、前記電流ベクトルと前記拾い上げ電圧指令と一次抵抗と拾い上げ制御開始指令とを入力して初期速度と初期二次磁束を出力する初期値推定器と、前記拾い上げ制御開始指令で同定時間ｔｘをカウントアップする同定用タイマと、前記電流ベクトルと前記拾い上げ電圧指令と一次抵抗ノミナル値を入力し演算磁束ベクトルを出力する磁束演算器と、時間ｔｘにおける前記演算磁束ベクトルを記憶する演算磁束メモリと、前記演算磁束メモリから任意の複数時点の演算磁束ベクトルを抽出する演算磁束抽出器と、前記任意の複数時点の演算磁束ベクトルと前記電流ベクトルを入力し前記一次抵抗を出力するＲ１Ａ同定器を持つ速度センサレス電動機制御装置において、
前記電流検出器と、前記拾い上げ制御手段と、前記初期値推定器と、前記同定用タイマと、前記磁束演算器と、前記演算磁束ベクトルを電流軸方向へ回転座標変換し回転磁束ベクトルを出力する回転座標変換器と、前記回転磁束ベクトルが極値となったときの極値時間txb1及びtxb2を出力する極値時間推定器と、前記同定時間ｔｘにおける前記回転磁束ベクトルを記憶する回転磁束メモリと、前記回転磁束メモリから前記極値時間を用いて複数時点の回転磁束ベクトルを抽出する極値磁束抽出器と、前記極値磁束抽出器の出力と前記同定時間ｔｘと前記電流ベクトルを入力し前記一次抵抗を出力するＲ１Ｂ同定器とから構成することを特徴とする速度センサレス電動機制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電動機のトルク制御に関するものであり、特に、電動機の初期速度と初期二次磁束の推定を高精度にし、電動機トルク制御を高精度にするものである。
【背景技術】
【０００２】
図４は、従来の一例を示すブロック図である。同図において、電動機６の拾い上げ制御の一例を説明する。１はトルク制御手段、２は拾い上げ制御手段、３は切替器、４は電力変換器、５は電流検出器、６は電動機、７は初期値推定器である。
電流検出器５は、電動機６に流れる電流ベクトルｉを検出する。拾い上げ制御手段２は、電流検出器５で検出した電流ベクトルｉを入力し、電動機６に流れる電流ｉを直流にする拾い上げ電圧指令ｖ０を出力する。
【０００３】
初期値推定器７は、電流ベクトルｉと拾い上げ電圧指令ｖ０と拾い上げ制御開始指令ＳＴと一次抵抗Ｒ１を入力し、初期速度ωｍ０と初期二次磁束φ２０を出力する。
【０００４】
トルク制御手段１は、トルク制御開始指令ＳＷが立つと、初期値推定器７の出力である初期速度ωｍ０と初期二次磁束φ２０を初期値として、電流ベクトルｉを基に、電動機６のトルクを制御するトルク制御電圧指令Ｖ１を出力する。
【０００５】
切替器３は、トルク制御開始指令ＳＷにより、拾い上げ電圧指令ｖ０とトルク制御電圧指令Ｖ１を切替える。すなわち、トルク制御開始指令ＳＷが立つまでは拾い上げ電圧指令ｖ０を電圧指令Ｖ＊とし、トルク制御開始指令ＳＷが立てばトルク制御電圧指令Ｖ１を電圧指令Ｖ＊として出力する。
電力変換器４は、電圧指令Ｖ＊を増幅して電動機６に電力を供給する。
【０００６】
以上の構成とすることにより、トルク制御開始指令ＳＷが立つまでは、拾い上げ制御手段２と初期値推定器７で電動機６の初期速度ωｍ０と初期二次磁束φ２０を推定する。トルク制御開始指令ＳＷが立てば、ＳＷが立った時点の初期速度ωｍ０と初期二次磁束φ２０を初期値として、電動機６のトルク制御が行われる。
なお、トルク制御開始指令ＳＷを立てるタイミングは、拾い上げ制御時間、初期速度ωｍ０と初期二次磁束φ２０の状態から決める。
【０００７】
以下、トルク制御開始指令ＳＷが立つまでの拾い上げ制御に関して、詳細に説明する。
【０００８】
拾い上げ制御手段２の一構成としては、図５に示すように、直流電流指令Ｉと電流位相角θと電流ベクトルｉから電動機６に流れる電流を直流に制御する拾い上げ電圧指令ｖ０を出力する電流制御手段８で構成するものがある。
初期値推定器７では、拾い上げ電圧指令ｖ０の動きから電動機６の初期速度ωｍ０と初期二次磁束φ２０を推定する。図６は、初期値推定器７の一構成例である。９は実磁束推定器、１０は実磁束メモリ、１１は実磁束抽出器、１２は初期速度推定器、１３は初期磁束推定器、１４は演算用タイマである。
【０００９】
演算用タイマ１４は、拾い上げ制御開始指令ＳＴのエッジで０クリアされるタイマカウンタであり、拾い上げ時間ｔ０を出力する。実磁束推定器９は、電流ベクトルｉ、拾い上げ電圧指令ｖ０、一次抵抗Ｒ１から式（Ａ）で実磁束推定ベクトルφ２ｒを演算する。
【数１】
2006014488.tif 000003
ここで、Ｌ１は一次自己インダクタンス、Ｌ２は二次自己インダクタンス、Ｍは相互インダクタンスである。
【００１０】
実磁束メモリ１０は、時々刻々変化する実磁束推定ベクトルφ２ｒの０〜ｔ０区間を記憶する。実磁束抽出器１１は、０〜ｔ０区間から、任意の３時点の実磁束推定ベクトルφ２ｒの値φ（ｔ００）、φ（ｔ０１）、φ（ｔ０２）を抽出する。
【００１１】
初期速度推定器１２は、一例として、式（Ｂ）〜式（Ｇ）で初期速度ωｍ０を演算する。最初に、3点を通る円の中心Ｒを式（Ｂ）〜式（Ｅ）で求める。
【数２】
2006014488.tif 000004
ここで、（Ｆ１Ａ、Ｆ１Ｂ）、（Ｆ２Ａ、Ｆ２Ｂ）、（ＲＡ、ＲＢ）は、各々、Ｆ１、Ｆ２、Ｒの各成分である。
次に、円の中心Ｒから見たφ（ｔ００）とφ（ｔ０２）の角度θＣを式（Ｆ）で求める。
【数３】
2006014488.tif 000005
最後に、式（Ｇ）で初期速度ωｍ０を求める。
【数４】
2006014488.tif 000006
【００１２】
初期磁束推定器１３は、式（Ｈ）にて初期二次磁束φ２０を求める。
【数５】
2006014488.tif 000007
ここで、ｊは虚数単位、Ｔ２は二次時定数で、Ｌ２／Ｒ２で得られる。Ｒ２は二次抵抗である。
【００１３】
以上の拾い上げ技術においては、以下に示す問題点がある。
電動機温度の変動等により一次抵抗Ｒ１に誤差があると、実磁束演算器９の出力値である演算磁束ベクトルφ２ｒの誤差が式（Ａ）により積算されていく。結果、実磁束抽出器１１の出力であるφ（ｔｘ０）、φ（ｔｘ１）、φ（ｔｘ２）に誤差が存在する。φ（ｔｘ０）、φ（ｔｘ１）、φ（ｔｘ２）を基に式（Ｂ）〜式（Ｈ）で推定する初期速度ωｍ０および初期二次磁束φ２０に誤差が生じる。
【００１４】
以上より、一次抵抗に誤差があると初期速度ωｍ０および初期二次磁束φ２０に誤差が生じ、トルク制御手段１の初期値に誤差があるために、電動機６のトルク制御が高精度に行うことができない。
【００１５】
そこで、図４の２１〜２５を追加し、式（Ａ）の一次抵抗Ｒ１を推定できるようにする。以下、図４の２１〜２５について説明する。
図４において、２１は同定用タイマ、２２は磁束演算器、２３は演算磁束メモリ、２４は演算磁束抽出器、２５はＲ１Ａ同定器である。
【００１６】
同定用タイマ２１は、拾い上げ制御開始指令ＳＴのエッジで０クリアされるタイマカウンタであり、同定時間ｔｘを出力する。磁束演算器２２は、電流ベクトルｉ、拾い上げ電圧指令v0、一次抵抗ノミナル値Ｒ１Ｃから式（Ｉ）で演算磁束ベクトルφ２ｓを演算する。
【数６】
2006014488.tif 000008
【００１７】
演算磁束メモリ２３は、時々刻々変化する演算磁束ベクトルφ２ｓの０〜ｔｘ区間を記憶する。演算磁束抽出器２４は、０〜ｔｘ区間から、等間隔の３時点の演算磁束ベクトルφ２ｓの値φｓ（ｔｘ０）、φｓ（ｔｘ１）、φｓ（ｔｘ２）を抽出する。特に、ｔｘ２を同定時間ｔｘとして、
【数７】
2006014488.tif 000009
と置く。
【００１８】
Ｒ１Ａ同定器２５は、φｓ（ｔｘ０）、φｓ（ｔｘ１）、φｓ（ｔｘ２）と同定時間ｔｘと電流ベクトルｉと一次抵抗ノミナル値Ｒ１Ｃを入力し、一次抵抗Ｒ１を出力する。以下、Ｒ１Ａ同定器２５について述べる。
【００１９】
式（Ｉ）において、一次抵抗ノミナル値Ｒ１Ｃが誤差ΔＲだけ真値Ｒ１からずれていたとする。
【数８】
2006014488.tif 000010
式（Ｉ）のＲ１Ｃに式（Ｋ）を代入すると、電流ベクトルｉは一定なので、式（Ｉ）は、式（Ａ）を用いて、式（Ｌ）で表される。
【数９】
2006014488.tif 000011
ここで、ｔはφ２ｓの演算時間である。また、Ｌ２≒Ｍとした。
【００２０】
そこで、式（Ｍ）〜式（Ｏ）でＥ０を計算する。
【数１０】
2006014488.tif 000012
ここで、Ｅ０ＡとＥ０ＢはＥ０をΔＲについてまとめた場合の係数である。
【００２１】
φ２ｒは円を描くので、原点とＦＳ１とＦＳ２は同一円上にある。さらに、原点とＦＳ１とＦＳ２は等間隔でサンプルしたものであるので、Ｅ０は０になる。ここで、一次抵抗誤差ΔＲはＥ０＝０を満たす式（Ｐ）で求めることができる。
ΔＲ＝−E0A/E0B                            式（Ｐ）
【００２２】
Ｒ１同定機能２１〜２５を新たに追加したことにより、一次抵抗Ｒ１の変動を同定することができ、実磁束推定器９で推定する実磁束推定ベクトルφ２ｒを高精度に求めることができる。その結果、電動機６の初期速度ωｍ０及び初期二次磁束φ２０を高精度に推定することができ、電動機トルクを高精度に制御することができる。（例えば、特許文献１参照。）。
【００２３】
【特許文献１】特開２００１−２１１６９７
【発明の開示】
【発明が解決しようとする課題】
【００２４】
従来技術においては、以下に示す問題点がある。
式(P)において、一次抵抗誤差ΔＲを求めることが発明の主点であったが、
・式(P)の分母であるE0Bが０となる場合がある、
・起動時に残留磁束が存在したり電動機速度が高いとき、あるいは、φｓ(tx0)とφｓ(tx2)間の相対ベクトルが電流ベクトルｉにほぼ平行の場合、E0A及びE0Bの大きさが小さくなる、
ことから、演算精度が著しく悪化することがある。このとき、初期速度ωｍ０および初期二次磁束φ２０に誤差が生じる。
【００２５】
初期速度ωｍ０および初期二次磁束φ２０に誤差が生じると、トルク制御手段１の初期値に誤差があるために、電動機６のトルク制御を高精度に行うことができない。初期速度ωｍ０および初期二次磁束φ２０の誤差が大きくなると、電動機６のトルク制御が不能となることもある。
本発明は、以上の問題点を解決するためになされたものである。
【課題を解決するための手段】
【００２６】
前述の問題点を解決するために以下の構成とする。
１、電動機６に流れる電流ベクトルｉを検知する電流検出器５
２、電流ベクトルｉを入力し電動機６に直流電流を流す拾い上げ電圧指令ｖ０を出力する拾い上げ制御手段２
３、電流ベクトルｉと拾い上げ電圧指令ｖ０と一次抵抗Ｒ１と拾い上げ制御開始指令ＳＴとを入力して初期速度ωｍ０と初期二次磁束φ２０を出力する初期値推定器７
４、拾い上げ制御開始指令ＳＴで同定時間ｔｘをカウントアップする同定用タイマ21
５、電流ベクトルｉと拾い上げ電圧指令ｖ０と一次抵抗ノミナル値Ｒ１Ｃを入力し演算磁束ベクトルφ２ｓを出力する磁束演算器22
６、演算磁束ベクトルφ２ｓを電流軸方向θへ回転座標変換し回転磁束ベクトルφ２ｔを出力する回転座標変換器31
７、回転磁束ベクトルφ２ｔが極値となったときの極値時間txb1及びtxb2を出力する極値時間推定器32
８、同定時間ｔｘ時点の回転磁束ベクトルφ２ｔを記憶する回転磁束メモリ35
９、回転磁束メモリ35から極値時間txb1及びtxb2を用いて複数時点の回転磁束ベクトルを抽出する極値磁束抽出器33
１０、極値磁束抽出器33の出力と同定時間ｔｘと電流ベクトルｉを入力し一次抵抗R1を出力するＲ１Ｂ同定器34
【００２７】
以上の構成により、一次抵抗変動ΔＲを高精度に同定することができ、その結果、電動機６の初期速度ωｍ０及び初期二次磁束φ２０を高精度に推定することができ、電動機トルクを高精度に制御することができる。
【発明の効果】
【００２８】
一次抵抗変動ΔＲを高精度に同定することができ、その結果、電動機６の初期速度ωｍ０及び初期二次磁束φ２０を高精度に推定することができ、電動機トルクを高精度に制御することができる。
【発明を実施するための最良の形態】
【００２９】
一次抵抗変動ΔＲを高精度に同定することにより、電動機トルクを高精度に制御する。
【実施例１】
【００３０】
図1は、本発明の一実施例を示すブロック図であり、31は回転座標変換器、35は回転磁束メモリ、32は極値時間推定器、33は極値磁束抽出器、34はＲ１Ｂ同定器である。
【００３１】
回転座標変換器31は、演算磁束ベクトルφ２ｓと電流位相角θを入力し、式(Q)と式(R)にて回転磁束ベクトルφ２ｔを演算する。
Ｆ２ｔＡ＝Ｆ２ｓＡ・cosθ＋Ｆ２ｓＢ・sinθ                     式(Q)
Ｆ２ｔＢ＝−Ｆ２ｓＡ・sinθ＋Ｆ２ｓＢ・cosθ                    式(R)
ここで、（Ｆ２ｓＡ、Ｆ２ｓＢ）、（Ｆ２ｔＡ、Ｆ２ｔＢ）は、各々、φ２ｓ及びφ２ｔの各成分である。φ２ｓに式(L)を適用すると、
Ｆ２ｔＡ＝Ｆ２ｒｔＡ＋ΔＲ・Ｉ・ｔ                              式(S)
Ｆ２ｔＢ＝Ｆ２ｒｔＢ                                            式(T)
となる。ここで、（Ｆ２ｒｔＡ、Ｆ２ｒｔＢ）は、実磁束推定ベクトルφ２ｒをθ方向へ回転座標変換したときの各成分である。回転座標変換器31を通すことにより、一次抵抗誤差ΔＲの影響がＦ２ｔＡのみとなり、Ｆ２ｔＢには一次抵抗誤差ΔＲの影響が無くなる。
【００３２】
回転磁束メモリ35は、同定時間ｔｘ時点の回転座標変換器31の出力である回転磁束ベクトルφ２ｔを記憶する。
【００３３】
極値時間推定器32は、図２に示すように、Ｆ２ｔＢの隣り合った極値地点の時間txb1及びtxb2を抽出する。なお、F2tB極値が1回しか抽出できていないときは0及び極値地点時間、F2tB極値を1回も抽出していないときはtxb1=txb2=0とする。
以下、txb1≦txb2として説明する。
【００３４】
極値磁束抽出器33は、回転磁束メモリ35から極値時間txb1及びtxb2を用いて3地点の回転磁束ベクトルφt(tx00)、φt(tx11)、φt(tx2)を抽出する。その詳細を以下に説明する。
同定時間ｔｘ=tx2がtxb2と(txb1+3*txb2)/2の間にある場合、極値時間txb1及びtxb2から式(U)、式(V)を用いて、tx00とtx11を計算する。
tx00=txb1                                                       式(U)
tx11=(txb1+tx2)/2                                               式(V)
同定時間ｔｘ=tx2が (txb1+3*txb2)/2以降にある場合、極値時間txb1及びtxb2から式(W)、式(X)を用いて、tx00とtx11を計算する。
tx00=txb2                                                       式(W)
tx11=(txb2+tx2)/2                                               式(X)
そして、回転磁束メモリ35から時間tx00、tx11、tx2の回転磁束ベクトル値φt(tx00)、φt(tx11)、φt(tx2)を抽出する。
【００３５】
Ｒ１Ｂ同定器31では、極値磁束抽出器33で抽出したφt(tx00)、φt(tx11)、φt(tx2)を式(M)〜式(P)に適用する。すなわち、φｓ（ｔｘ０）、φｓ（ｔｘ１）、φｓ（ｔｘ２）の代わりに、φt(tx00)、φt(tx11)、φt(tx2)を用いる。すると、ΔＲは式(Y)〜式(AA)で求められ、式(AB) にてＲ１を同定する。
F11A=Fx11A−Fx00A、F11B=Fx11B−Fx00B                            式(Y)
F22A=Fx22A−Fx00A、F22B=Fx22B−Fx00B                            式(Z)
ΔＲ＝{F22A*(2*F11A−F22A)+F22B*(2*F11B−F22B)}/(2*F11A−F22A)/tx2/I  式(AA)
Ｒ１＝Ｒ１Ｃ＋ΔＲ                                                  式(AB)
ここで、Fx00AとFx00Bはφt(tx00)の成分、Fx11AとFx110Bはφt(tx11)の成分、Fx22AとFx22Bはφt(tx2)の成分である。
【００３６】
以上の過程で一次抵抗変動ΔＲを高精度に推定することができる。その理由を以下に述べる。
図３は、従来と本発明における一次抵抗変動ΔＲを推定する際に用いる式(AA)に相当する演算式の分母値状態を比較した図である。R1A同定器演算式分母は(2*FmA−FeA)の値をプロットしたもの、R1B同定器演算式分母は(2*F11A−F22A)の値をプロットしたものである。ここで、
FmA=F2tA(tx1)−F2tA(tx0)                                        式(AC)
FeA=F2tA(tx2)−F2tA(tx0)                                        式(AD)
である。F2tA(tx0)は時間tx0でのF2tA値、F2tA(tx1)は時間tx1でのF2tA値、F2tA(tx2)は時間tx2でのF2tA値である。
【００３７】
tx1とtx0は、従来のＲ１同定式で用いた時間であり、(2*FmA−FeA)が従来のR1A同定器の演算式分母値に相当する。式(AA)で示したように、(2*F11A−F22A)は本発明のR1B同定器の演算式分母値となる。
【００３８】
R1A同定器演算式分母とR1B同定器演算式分母を比較すると、R1A同定器演算式分母値がよく０近くの値になるのに対し、R1B同定器演算式分母値は始点を除いて０となることはない。始点部分でＲ１推定をしなければ、R1B同定器34で一次抵抗変動ΔＲを高精度に推定することができる。
【００３９】
また、始点部分でＲ１推定が必要となるのは、低速域である。しかし、低速では、演算磁束ベクトルφ２ｓの軌跡がほぼ直線であるので、Ｒ１推定を実施しなくとも初期速度ωｍ０に大きな誤差は発生しない。よって、始点部分でＲ１推定をしなくとも問題は無い。
【００４０】
以上の構成とすることにより、一次抵抗変動ΔＲを高精度に推定することができ、その結果、電動機６の初期速度ωｍ０及び初期二次磁束φ２０を高精度に推定することができ、電動機トルクを高精度に制御することができる。
【産業上の利用可能性】
【００４１】
速度センサレスにおいても、電動機の運転開始時の惰行速度を速やかに、かつ高精度に推定することができる。
例えば、速度センサレス車両制御において、惰行走行状態からの再加速、惰行走行状態からのブレーキが可能となる。
【図面の簡単な説明】
【００４２】
【図１】図1は、本発明の一実施例を示すブロック図である。
【図２】図２は、本発明の極値時間推定器の一演算例を示す図である。
【図３】図３は、本発明と従来のＲ１同定演算分母値の比較をした図である。
【図４】図４は、一従来例を示すブロック図である。
【図５】図５は、一従来例の拾い上げ制御手段の一構成例を示すブロック図である。
【図６】図６は、一従来例の初期値推定器の一構成例を示すブロック図である。
【符号の説明】
【００４３】
１    トルク制御手段
２    拾い上げ制御手段
３    切替器
４    電力変換器
５    電流検出器
６    電動機
７    初期値推定器
８      電流制御手段
９      実磁束推定器
１０    実磁束メモリ
１１    実磁束抽出器
１２    初期速度推定器
１３    初期磁束推定器
１４    演算用タイマ
２１  同定用タイマ
２２  磁束演算器
２３  演算磁束メモリ
２４  演算磁束抽出器
２５  Ｒ１Ａ同定器
３１  回転座標変換器
３２  極値時間推定器
３３  極値磁束抽出器
３４  Ｒ１Ｂ同定器
３５  回転磁束メモリ
v0・・・・拾い上げ電圧指令
V1・・・・トルク制御電圧指令
Ｖ＊・・・・電圧指令
ｉ・・・・電流ベクトル
ωｍｒ・・・・実速度
ＳＷ・・・・トルク制御開始指令
ＳＴ・・・・拾い上げ制御開始指令
Ｒ１・・・・一次抵抗
ωｍ０・・・・初期速度
φ２０・・・・初期二次磁束
Ｉ・・・・直流電流指令
θ・・・・電流指令位相角
ｔ０・・・・拾い上げ時間
φ（ｔ０）・・・・ｔ０時点の実磁束推定ベクトル
ｔ００、ｔ０１、ｔ０２・・・・０以上ｔ０以下の任意の値
φ（ｔ００）・・・・ｔ００時点の実磁束推定ベクトル
φ（ｔ０１）・・・・ｔ０１時点の実磁束推定ベクトル
φ（ｔ０２）・・・・ｔ０２時点の実磁束推定ベクトル
φ２ｒ・・・・実磁束推定ベクトル
ｔｘ・・・・同定時間
φ２s・・・・演算磁束ベクトル
Ｒ１Ｃ・・・・一次抵抗ノミナル値
ΔＲ・・・・一次抵抗誤差
φs(tx)・・・・tx時点の演算磁束ベクトル
tx0、tx1、tx2・・・・0以上tx以下の任意の値
φs(tx0)・・・・tx0時点の演算磁束ベクトル
φs(tx1)・・・・tx1時点の演算磁束ベクトル
φs(tx2)・・・・tx2時点の演算磁束ベクトル
φ２ｔ・・・・回転磁束ベクトル
φｔ(tx)・・・・tx時点のφ２ｔ
F2tA・・・・回転磁束ベクトルのθ方向成分
F2tB・・・・回転磁束ベクトルのθ方向直交成分
txb1、txb2・・・・F2tBの極値時間
tx00、tx11・・・・Ｒ１同定用磁束抽出時間
φｔ(tx00)・・・・tx00時点の回転磁束ベクトル
φｔ(tx11)・・・・tx11時点の回転磁束ベクトル
φｔ(tx2)・・・・tx2時点の回転磁束ベクトル
【図１】
2006014488.tif 000013
【図２】
2006014488.tif 000014
【図３】
2006014488.tif 000015
【図４】
2006014488.tif 000016
【図５】
2006014488.tif 000017
【図６】
2006014488.tif 000018