(57)【要約】【課題】各種鋼材の長期の腐食量を簡単に且つ高精度に予測することを可能にした鋼材の寿命予測方法及びその装置、並びにそらの方法又は装置をコンピュータを用いて実現するためのコンピュータプログラムを提供する。【解決手段】鋼材の腐食量予測式Ｙ＝ＡＸB（Ｙ：腐食量、Ｘ：年数、Ａ，Ｂ：材料と環境に依存する係数、べき数）を用いて鋼材の寿命を予測する方法において、次式に示すように、前記係数Ａを、鋼種に応じて予め設定された係数（α，β，γ）と、構造物の設置環境に依存した環境データ（Ｔ，ＰＷまたはＴＯＷ，Ｓａ）とに基づいて求め、前記べき数係数Ｂを前記係数Ａの関数として求める。    Ａ＝（α・Ｔ＋β）・ＰＷ（Ｔ，Ｈ）・（Ｓaγ）            Ｔ：温度（℃），Ｈ：相対湿度（％），Ｓａ：飛来塩分量（mdd）            ＰＷ（Ｔ，Ｈ）：濡れ確率、ＴＯＷ：年間濡れ時間（ｈ）            α、β、γ：鋼種に応じて設定された係数    Ｂ＝ｆ（Ａ）、ただし、ｆ（Ａ）はＡの関数であることを示す。【選択図】        図８
2006208346.tif 000002
【特許請求の範囲】
【請求項１】
鋼材の腐食量予測式Ｙ＝ＡＸB（Ｙ：腐食量、Ｘ：年数、Ａ，Ｂ：材料と環境に依存する係数、べき数）を用いて鋼材の寿命を予測する方法において、
前記係数Ａを、鋼種に応じて予め設定された係数（α，β，γ）と、前記鋼材を用いる構造物の設置環境に依存した環境データ（Ｔ，ＰＷまたはＴＯＷ，Ｓａ）とに基づいて求め、前記べき数係数Ｂを前記係数Ａの関数として求めることを特徴とする鋼材の寿命予測方法。ただし、Ｔ：温度（℃）、ＰＷ：濡れ確率、Ｓａ：飛来塩分量(mdd)、ＴＯＷ：年間漏れ時間（ｈ）である。
【請求項２】
前記Ａ及びＢを次式より求めることを特徴とする請求項１記載の鋼材の寿命予測方法。
Ａ＝（α・Ｔ＋β）・ＰＷ（Ｔ，Ｈ）・（Ｓaγ）  
Ｔ：温度（℃），Ｈ：相対湿度（％），Ｓａ：飛来塩分量（mdd）
ＰＷ（Ｔ，Ｈ）：濡れ確率、
α、β、γ：鋼種に応じて設定された係数
Ｂ＝ｆ（Ａ）、ただし、ｆ（Ａ）はＡの関数であることを示す。
【請求項３】
前記Ａ及びＢを次式より求めることを特徴とする請求項１記載の鋼材の寿命予測方法。
Ａ＝ｋ（α・Ｔ＋β）・ＴＯＷ・（Ｓaγ）
Ｔ：温度（℃），Ｓａ：飛来塩分量（mdd）
ＴＯＷ：年間濡れ時間（ｈ）、ｋ：係数
α、β、γ：鋼種に応じて設定された係数
Ｂ＝ｆ（Ａ）、ただし、ｆ（Ａ）はＡの関数であることを示す。
【請求項４】
前記Ｂを、０．５≦Ｂ≦１．１の範囲で設定することを特徴とする請求項１〜３の何れかに記載の鋼材の寿命予測方法。
【請求項５】
前記Ｂを、
Ａ＜0.083のとき、Ｂ＝−4611.3Ａ3＋769.19Ａ2−32.421Ａ＋1.0109
0.083≦Ａのとき、Ｂ＝0.9〜1.1
に設定することを特徴とする請求項４記載の鋼材の寿命予測方法。
【請求項６】
前記Ｂを、
Ａ≦0.03のとき、Ｂ＝0.5〜0.7、
0.03＜Ａ＜0.083のとき、Ｂ＝−4611.3Ａ3＋769.19Ａ2−32.421Ａ＋1.0109
0.083≦Ａのとき、Ｂ＝0.9〜1.1
に設定することを特徴とする請求項４記載の鋼材の寿命予測方法。
【請求項７】
鋼種を入力手段を介して入力する工程と、
前記鋼種に対応した前記係数（α、β、γ）を特定する工程と、
構造物を設置する場所についての情報を入力手段を介して入力する工程と、
前記構造物の設置場所又はその周辺の温度（Ｔ）及び相対湿度（Ｈ）を求める工程と、
前記の温度（Ｔ）及び相対湿度（Ｈ）に基づいて濡れ確率（ＰＷ）を求める工程と、
前記構造物の設置場所の離岸距離を入力手段を介して入力する工程と、
前記構造物の設置場所を含む領域又はその周辺の領域において予め計測されて設定されている飛来塩分量及び前記離岸距離に基づいて前記構造物の設置場所の飛来塩分量（Ｓａ）を求める工程と
を有することを特徴とする請求項２、４、５又は６記載の鋼材の寿命予測方法。
【請求項８】
鋼種を入力手段を介して入力する工程と、
前記鋼種に対応した前記係数（α、β、γ）を特定する工程と、
構造物を設置する場所についての情報を入力手段を介して入力する工程と、
前記構造物の設置場所又はその周辺の温度（Ｔ）及び相対湿度（Ｈ）を求める工程と、
前記の温度（Ｔ）及び相対湿度（Ｈ）に基づいて年間濡れ時間（ＴＯＷ）を求める工程と、
前記構造物の設置場所の離岸距離を入力手段を介して入力する工程と、
前記構造物の設置場所を含む領域又はその周辺の領域において予め計測されて設定されている飛来塩分量及び前記離岸距離に基づいて前記構造物の設置場所の飛来塩分量（Ｓａ）を求める工程と
を有することを特徴とする請求項３、４、５又は６記載の鋼材の寿命予測方法。
【請求項９】
前記構造物の設置場所を含む領域又はその周辺の領域を、前記構造物の設置場所の入力とは別に、入力手段を介して入力する工程を更に有することを特徴とする請求項７又は８記載の鋼材の寿命予測方法。
【請求項１０】
前記構造物の設置場所の飛来塩分（Ｓａ）に構造物の部位に対応した部位係数を乗算して、その乗算結果を前記飛来塩分（Ｓａ）とすることを特徴とする請求項７又は８記載の鋼材の寿命予測方法。
【請求項１１】
鋼材の腐食量予測式Ｙ＝ＡＸB（Ｙ：腐食量、Ｘ：年数、Ａ，Ｂ：材料と環境に依存する係数、べき数）を用いて鋼材の寿命を予測する装置において、
鋼種、前記鋼材を用いる構造物の設置場所及び構造物の設置場所の離岸距離をそれぞれ取り込む入力手段と、
鋼材の鋼種に対応した係数（α、β、γ）、各地域の温度（Ｔ）及び相対湿度（Ｈ）、及び特定地域の飛来塩分量がそれぞれ格納された記憶手段と、
前記入力手段を介して入力された情報と、前記記憶手段に格納されたデータとに基づいて前記係数Ａ及び前記べき係数Ｂを求め、更に前記係数Ａ及び前記べき係数Ｂに基づいて前記腐食量Ｙを求める演算手段と
を備えたことを特徴とする鋼材の寿命予測装置。
【請求項１２】
前記演算手段は、前記Ａ及びＢを次式により求めることを特徴とする請求項１１記載の鋼材の寿命予測装置。
Ａ＝（α・Ｔ＋β）・ＰＷ（Ｔ，Ｈ）・（Ｓaγ）
Ｔ：温度（℃），Ｈ：相対湿度（％），Ｓａ：飛来塩分量（mdd）
ＰＷ（Ｔ，Ｈ）：濡れ確率、
α、β、γ：鋼種に応じて設定された係数
Ｂ＝ｆ（Ａ）、ただし、ｆ（Ａ）はＡの関数であることを示す。
【請求項１３】
前記演算手段は、前記Ａ及びＢを次式により求めることを特徴とする請求項１０記載の鋼材の寿命予測装置。
Ａ＝（α・Ｔ＋β）・ＴＯＷ・（Ｓaγ）
Ｔ：温度（℃），Ｓａ：飛来塩分量（mdd）
ＴＯＷ：年間濡れ時間（ｈ）、
α、β、γ：鋼種に応じて設定された係数
Ｂ＝ｆ（Ａ）、ただし、ｆ（Ａ）はＡの関数であることを示す。
【請求項１４】
前記演算手段は、上記の腐食量の予測値の精度を、予測値に付帯して表示装置に表示させることを特徴とする請求項１１〜１３の何れかに記載の鋼材の寿命予測装置。
【請求項１５】
請求項１〜６記載の何れかに記載の鋼材の寿命予測方法における演算処理をコンピュータによる実行させることを特徴とするコンピュータプログラム。
【請求項１６】
請求項１１〜１４の何れかに記載の鋼材の寿命予測装置の演算手段をコンピュータによる実行させることを特徴とするコンピュータプログラム。
【請求項１７】
請求項１１〜１４の何れかに記載の鋼材の寿命予測装置をコンピュータにより実現し、演算結果を表示装置に表示させて顧客に提示し、顧客の購入を促すことを特徴とする鋼材の営業方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、構造物用の鋼材の寿命予測方法及びその装置並びにコンピュータプログラムに関する。
【背景技術】
【０００２】
従来、建築・土木分野で使用される鋼材として、ＳＳ、ＳＭと呼ばれる普通鋼とともに、無塗装で使用されるＳＭＡと呼ばれる耐候性鋼があり、最近ではＮｉ系高耐候性鋼等がある。これらの材料を用いた構造物の適用基準は、これまでには
(１）飛来海塩による地域区分          耐候性鋼
(２）鉄骨構造建築物の耐久性向上技術  普通鋼
(３）当該鋼種の暴露試験              任意の鋼種
(４）腐食量予測式                    耐候性鋼
等があり、これらは必要に応じて使い分けられている。
【０００３】
上記の（１）及び（２）は従来の暴露結果を纏め上げたもので、適用可否判断が瞬時に可能である。
上記の（３）の暴露試験による鋼材の腐食量は、
Ｙ＝ＡＸB、Ｙ：腐食量、Ｘ：年数
の式で表されることが知られている。腐食量のしきい値Ｙ1imを設定することにより、そのときのＸ1im年を寿命とする方法である。腐食寿命の判断の目安は、例えば１００年の推定片側板厚減少量が０．５mm以下である。この式のＡ、Ｂは環境や鋼種によって変化するため、Ａ、Ｂを決定するために、実環境又は実環境に近い環境に試験片を暴露し、試験片の腐食量の経年変化（Ｘ，Ｙ）から累乗近似する方法が用いられている。
【０００４】
ところで、腐食速度は、海塩量や亜硫酸ガス量によって影響を受けることが多数の文献に示されている。また、腐食現象は基本的には水溶液中の化学反応であるので、気温、湿度や、濡れ時間にも依存する。したがって、これらの環境因子をパラメータとする関数として記述することができる。上記の（４）の予測方法に関して、建設省（当時）土木研究所においては、ＳＭＡ（ＪＩＳ規格の耐候性鋼）に関し、暴露試験の結果から、上記のＡ，Ｂを次のように決定している（例えば非特許文献１）。
Ａ＝ＣＳaγ       Ｓa：飛来塩分量、Ｃ，γ：回帰係数
Ｂ＝０．７３
即ち、Ａは飛来塩分と相関があるとし、Ｂは鋼種・暴露方向、設置場所等によらず０．７３としている。
【０００５】
また、上記（４）の予測方法に関して、例えば次のように求める方法も提案されている（例えば特許文献１）。この方法では、腐食性指標Ｚの２次回帰式から初年度腐食量を推定している。 
Ｚ=α・ＴＯＷ・exp（-κ・W）・（C+δ・S）/（1+ε・C・S）・exp（-Ea/RT）
α＝106, κ＝-0.1, δ＝0.05, ε＝10.0, Ea＝50ｋJ/mol・K
Ｚ：腐食性指標，R：気体定数, C：飛来塩分, S：硫黄酸化物量, 
TOW：年間濡れ時間（ｈ），W：年平均風速（ｍ／sec），Ｔ：年平均気温（Ｋ）
【非特許文献１】「耐候性鋼材の橋梁への適用に関する共同研究報告書(XVIII)」（建設省土木研究所，（社）鋼材倶楽部，（社）日本橋梁建設協会、平成５年３月発行）
【特許文献１】国際公開03/006957号パンフレット
【発明の開示】
【発明が解決しようとする課題】
【０００６】
上記の従来の寿命予測方法の内、上記の（１）及び（２）の方法は鋼種が限定される上、或る時点でのその材料の適用可否を判断するだけであって、鋼種の選定は出来なかった。一方、上記の（３）の方法は、各鋼種の調査を行なうことによって鋼種の選定が可能である。しかしながら、新たな高耐候性を有する鋼が開発された場合には、開発後にただちに暴露試験を長期間種々の場所で行って、検量線あるいは式を作成しなければならない。また、上記の（３）の方法では、未知数Ａ，Ｂ２つを決定するために、最低２点以上の腐食量データが必要である。日本では四季があり、腐食は季節によって進行速度が異なるために、最低２年間の計測を必要とする。更に、近年、土木建築用鋼材の長寿命化が叫ばれており、例えば１００年といった長期寿命を精度よく予測するには、実質５〜２０年といった長期試験を行っており、実際に用いられるまでに時間がかかるという課題があった。また、上記の（４）の寿命予測方法においては、架設地の環境の影響は取り入れて演算式を構成しているものの、長期の暴露試験が必要であり、更に、非特許文献１の方法ではＢを一定としているので、さびの保護性に必要な飛来塩分量や鋼種の影響が反映されていない、という問題点があった。
【０００７】
また、上記の特許文献１において提案されている方法では、Ｓ（硫黄酸化物量）を予測式のパラメータとしているが、Ｓ（硫黄酸化物量）と板厚減少量との間には有意な相関が見られないことが知られている。
【０００８】
本発明は、このような問題点を解決するためになされたものであり、各種鋼材の長期の腐食量を簡単に且つ高精度に予測することを可能にした鋼材の寿命予測方法及びその装置並びにそれらの方法又は装置をコンピュータを用いて実現するためのコンピュータプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明に係る鋼材の寿命予測方法は、鋼材の腐食量予測式Ｙ＝ＡＸB（Ｙ：腐食量、Ｘ：年数、Ａ，Ｂ：材料と環境に依存する係数、べき数）を用いて鋼材の寿命を予測する方法において、前記係数Ａを、鋼種に応じて予め設定された係数（α，β，γ）と、前記鋼材を用いる構造物の設置環境に依存した環境データ（Ｔ，ＰＷまたはＴＯＷ，Ｓａ）とに基づいて求め、前記べき数係数Ｂを前記係数Ａの関数として求める。
但し、Ｔ：温度（℃）、ＰＷ：濡れ確率、Ｓａ：飛来塩分量(mdd)、ＴＯＷ：年間漏れ時間（ｈ）である。
【００１０】
本発明に係る鋼材の寿命予測方法は、前記Ａ及びＢを次式により求めるものである。
Ａ＝（α・Ｔ＋β）・ＰＷ（Ｔ，Ｈ）・（Ｓaγ）
Ｔ：温度（℃），Ｈ：相対湿度（％），Ｓａ：飛来塩分量（mdd）
ＰＷ（Ｔ，Ｈ）：濡れ確率、
α、β、γ：鋼種に応じて設定された係数
Ｂ＝ｆ（Ａ）、ただし、ｆ（Ａ）はＡの関数であることを示す。
【００１１】
本発明に係る鋼材の寿命予測方法は、前記Ａ及びＢを次式により求めるものである。
Ａ＝ｋ（α・Ｔ＋β）・ＴＯＷ・（Ｓaγ）
Ｘ：暴露時間（ｙ），Ｙ：板厚減少量（mm）
Ｔ：温度（℃），Ｓａ：飛来塩分量（mdd）
ＴＯＷ：年間濡れ時間（ｈ）、ｋ：係数
α、β、γ：鋼種に応じて設定された係数
Ｂ＝ｆ（Ａ）、ただし、ｆ（Ａ）はＡの関数であることを示す。
【００１２】
本発明に係る鋼材の寿命予測方法は、前記Ｂを、０．５≦Ｂ≦１．１の範囲で設定する。
【００１３】
本発明に係る鋼材の寿命予測方法は、前記Ｂを、
Ａ＜0.083のとき、Ｂ＝−4611.3Ａ3＋7,9.19Ａ2−32.421Ａ＋1.0109
0.083≦Ａのとき、Ｂ＝0.9〜1.1
に設定する。
【００１４】
本発明に係る鋼材の寿命予測方法は、前記Ｂを、
Ａ≦0.03のとき、Ｂ＝0.5〜0.7、
0.03＜Ａ＜0.083のとき、Ｂ＝−4611.3Ａ3＋769.19Ａ2−32.421Ａ＋1.0109
0.083≦Ａのとき、Ｂ＝0.9〜1.1
に設定する。
【００１５】
本発明に係る鋼材の寿命予測方法は、鋼種を入力手段を介して入力する工程と、前記鋼種に対応した前記係数（α、β、γ）を特定する工程と、構造物を設置する場所についての情報を入力手段を介して入力する工程と、前記構造物の設置場所又はその周辺の温度（Ｔ）及び相対湿度（Ｈ）を求める工程と、前記の温度（Ｔ）及び相対湿度（Ｈ）に基づいて濡れ確率（ＰＷ）を求める工程と、前記構造物の設置場所の離岸距離を入力手段を介して入力する工程と、前記構造物の設置場所を含む領域又はその周辺の領域において予め計測されて設定されている飛来塩分量及び前記離岸距離に基づいて前記構造物の設置場所の飛来塩分量（Ｓａ）を求める工程とを有する。
【００１６】
本発明に係る鋼材の寿命予測方法は、鋼種を入力手段を介して入力する工程と、前記鋼種に対応した前記係数（α、β、γ）を特定する工程と、構造物を設置する場所についての情報を入力手段を介して入力する工程と、前記構造物の設置場所又はその周辺の温度（Ｔ）及び相対湿度（Ｈ）を求める工程と、前記の温度（Ｔ）及び相対湿度（Ｈ）に基づいて年間濡れ時間（ＴＯＷ）を求める工程と、前記構造物の設置場所の離岸距離を入力手段を介して入力する工程と、前記構造物の設置場所を含む領域又はその周辺の領域において予め計測されて設定されている飛来塩分量及び前記離岸距離に基づいて前記構造物の設置場所の飛来塩分量（Ｓａ）を求める工程とを有する。
【００１７】
本発明に係る鋼材の寿命予測方法は、前記構造物の設置場所を含む領域又はその周辺の領域を、前記構造物の設置場所の入力とは別に、入力手段を介して入力する工程を、更に有する。
【００１８】
本発明に係る鋼材の寿命予測方法は、前記構造物の設置場所の飛来塩分（Ｓａ）に構造物の部位に対応した部位係数を乗算して、その乗算結果を前記飛来塩分（Ｓａ）とする。
【００１９】
本発明に係る鋼材の寿命予測装置は、構造物の鋼材の腐食量予測式Ｙ＝ＡＸB（Ｙ：腐食量、Ｘ：年数、Ａ，Ｂ：材料と環境に依存する係数、べき数）を用いて鋼材の寿命を予測する装置において、鋼種、前記鋼材を用いる構造物の設置場所及び構造物の設置場所の離岸距離をそれぞれ取り込む入力手段と、鋼材の鋼種に対応した係数（α、β、γ）、各地域の温度（Ｔ）及び相対湿度（Ｈ）、及び特定の地域の飛来塩分量がそれぞれ格納された記憶手段と、前記入力手段を介して入力された情報と、前記記憶手段に格納されたデータとに基づいて前記係数Ａ及び前記べき係数Ｂを求め、更に前記係数Ａ及び前記べき係数Ｂに基づいて前記腐食量Ｙを求める演算手段とを備えたものである。
【００２０】
本発明に係る鋼材の寿命予測装置において、前記演算手段は、前記Ａ及びＢを次式により求める。
Ａ＝（α・Ｔ＋β）・ＰＷ（Ｔ，Ｈ）・（Ｓaγ）
Ｔ：温度（℃），Ｈ：相対湿度（％），Ｓａ：飛来塩分量（mdd）
ＰＷ（Ｔ，Ｈ）：濡れ確率、
α、β、γ：鋼種に応じて設定され係数
Ｂ＝ｆ（Ａ）、ただし、ｆ（Ａ）はＡの関数であることを示す。
【００２１】
本発明に係る鋼材の寿命予測装置において、前記演算手段は、前記Ａ及びＢを次式により求める。
Ａ＝ｋ（α・Ｔ＋β）・ＴＯＷ・（Ｓaγ）
Ｔ：温度（℃），Ｓａ：飛来塩分量（mdd）
ＴＯＷ：年間濡れ時間（ｈ）、
α、β、γ：鋼種に応じて設定された係数
Ｂ＝ｆ（Ａ）、ただし、ｆ（Ａ）はＡの関数であることを示す。
【００２２】
本発明に係る鋼材の寿命予測装置において、前記演算手段は、上記の腐食量の予測値の精度を、予測値に付帯して表示装置に表示させる。
【００２３】
本発明に係るコンピュータプログラムは、上記の鋼材の寿命予測方法における演算処理をコンピュータによる実行させるものである。
【００２４】
本発明に係るコンピュータプログラムは、上記の鋼材の寿命予測装置の演算手段をコンピュータによる実行させるものである。
【００２５】
本発明に係る営業方法は、上記の鋼材の寿命予測装置をコンピュータにより実現し、演算結果を表示装置に表示させ顧客に提示し、顧客の購入を促す。
【発明の効果】
【００２６】
本発明によれば、鋼材の腐食量予測式Ｙ＝ＡＸB（Ｙ：腐食量、Ｘ：年数、Ａ，Ｂ：材料と環境に依存する係数、べき数）を用いて鋼材の寿命を予測する方法において、前記係数Ａを、鋼種に応じて予め設定された係数（α，β，γ）と、前記鋼材を用いる構造物の設置環境に依存した環境データ（Ｔ，ＰＷまたはＴＯＷ，Ｓａ）とに基づいて求め、前記べき数係数Ｂを前記係数Ａの関数として求めるようにしたので、各種鋼材の長期の腐食量を簡単に且つ高精度に予測することが可能になっている。
【発明を実施するための最良の形態】
【００２７】
実施形態１．
本発明の実施形態１として、鋼材の腐食量予測式Ｙ＝ＡＸB（Ｙ：腐食量、Ｘ：年数、Ａ，Ｂ：材料と環境に依存する係数、べき数）を用いた鋼材の寿命予測をする方法を説明するが、それに先だって、まず、本実施形態１に係る寿命予測方法の計測原理を説明する。
【００２８】
図１は期間Ｘと鋼材の腐食量Ｙとの関係を示した特性図であり、腐食寿命の判断の目安は、例えば１００年の推定片側板厚減少量が０．５mm以下である。この期間Ｘと鋼材の腐食量Ｙとの関係は、上述のように次式により表される。なお、Ａ値は構造物の架設地点の環境での鋼材自体の耐食性を示しており、Ｂ値はさびの保護性を表している。鋼材の耐食性が高ければ図１の特性の初期の傾きは小さく、さびの保護性が高ければ腐食量Ｙの長期の年数経過後の値は小さな値を示す。
Ｙ＝ＡＸB、Ｙ：腐食量、Ｘ：年数                        …（１）
【００２９】
特に、本実施形態１は、上記のＡ値は鋼種に応じた係数（α，β，γ）と環境データ（Ｔ，ＰＷまたはＴＯＷ，Ｓａ）とによって特定されること、また、Ｂ値はＡ値の関数として特定されることが見出された結果実現されたものであり、例えばＡ値及びＢ値はそれぞれ次式により求められる。
Ａ＝ＣＳaγ＝（α・Ｔ＋β）・ＰＷ（Ｔ，Ｈ）・（Ｓaγ）  …（２）
Ｂ＝ｆ（Ａ）
Ａ≦０．０３のとき、
Ｂ＝０．５〜０．７の範囲内で任意に設定、望ましくは０．６（理論的には放
物線則0.5乗と考えられるが、実際に形成れるさびは完全に緻密でないため、
０．５〜０．７で変動する。）
例えば、初期値としてＢ＝０．６と設定し、寿命予測精度の更なる向上等の必要に応じて、Ｂ＝０．５〜０．７の範囲内で調整すればよい。
０．０３＜Ａ＜０．０８３のとき、
Ｂ=−４６１１．３Ａ3＋７６９．１９Ａ2−３２．４２１Ａ＋１．０１０９
ただし、より簡便とするためにＡ≦０．０３のときにおいても、上記の式（Ｂ=−４６１１．３Ａ3＋７６９．１９Ａ2−３２．４２１Ａ＋１．０１０９）を適用してもよい。
０．０８３≦Ａのとき、
Ｂ＝０．９〜１．１の範囲内で任意に設定、望ましくは１      …（３）
例えば、初期値としてＢ＝１．０と設定し、寿命予測精度の更なる向上等の必要に応じて、Ｂ＝０．９〜１．１の範囲内で調整すればよい。
Ｔ：温度（℃），Ｈ：相対湿度（％），Ｓａ：飛来塩分量（mdd），
ＰＷ（Ｔ，Ｈ）：濡れ確率
α，β，γ：鋼種に応じた係数
なお、上記（２）式において、温度Ｔに関する項は、実際の対象となる温度範囲が比較的狭いので、直線近似にしてある。また、湿度Ｈに関する項は、腐食量が濡れ時間に比例すると考え、KuceraらがＩＳＯに提案した濡れ確率関数を導入した（Kucera, Tidblad, Mikhailov: ISO/TC156/WG4-N314, Annex A (1999)）。これは温度、湿度ごとにＰ（Ｔ），Ｐ（Ｈ）が数値表で与えられ、両者の積が濡れ確率（ＰＷ＝Ｐ（Ｔ）×Ｐ（Ｈ））となる。年間の濡れ時間は８７６６ｈ×ＰＷ（濡れ確率）である。なお、Ｐ（Ｔ）及びＰ（Ｈ）の数値表は次のとおりである。なお、この濡れ確率ＰＷ（Ｔ，Ｈ）は、kuceraの式によると、濡れ確率ＰＷ（Ｔ，Ｈ）＝N (T；0；9.96)*β(H/100；4.67；1.78)で表される。表１にＰ（Ｔ），Ｐ（Ｈ）の数値表の一部を抜枠して示す。
【００３０】
【表１】
2006208346.tif 000003
【００３１】
また、濡れ確率ＰＷ（Ｔ，Ｈ）と年間濡れ時間TOW（ｈ）とは、
TOW（ｈ）＝8766×ＰＷ（Ｔ，Ｈ）
で表されるから、上記（２）式は年間濡れ時間TOWの関数として次のように表される。
Ａ＝ｋ（α・Ｔ＋β）・ＴＯＷ・（Ｓaγ）                       …（２ａ）
ＴＯＷ：年間濡れ時間（ｈ）、ｋ：係数
【００３２】
なお、上記の温度Ｔ、相対湿度Ｈ及び飛来塩分量Ｓaは、例えば橋梁が設置される環境の値である。飛来塩分量Ｓaには、その環境での風の影響や方向（どの方角を向いているか）、橋桁の高さ等の影響が含まれた値であり、事前に測定されたものである。例えば日本のある地点（緯度・経度又は住所）を決定すれば、その位置情報に対応した温度Ｔ、相対湿度Ｈ、飛来塩分量Ｓaを求めることができる。例えば地点情報と温度Ｔ、相対湿度Ｈの気象庁データ平年値とが既にデータベース化されており、それらのデータベースを利用することにより必要な情報が得られる。また、飛来塩分量Ｓaについては、約３００地点での測定データを集約した関係式（各地方の飛来海塩量・離岸距離の関係。具体的にはＳa＝Ｓa0・Ｌ-0.6、Ｌ：離岸距離）があり、この関係式を用いることにより、該当地点における飛来塩分量Ｓaを求めることができる。なお、これらのデータをデータベースから入手できない場合には、該当する地点において実測により求めてもよい。このようにして、該当する地点における温度Ｔ、相対湿度Ｈ及び飛来塩分量Ｓaを求めることができるので、係数α、β、γを実験室での試験で決定することができれば、Ａ値が求まり、Ｂ値はＡ値の関数であり、Ｙ＝ＡＸBによりＹ（板厚減少量）を求めることができる。
【００３３】
ここで、鋼種に応じた係数α，β，γの求め方についてその概要を説明する。図２はＳＭＡ実暴露の試験データにより求めたＡ値の特性図である。上記の（１）式において、Ｘ＝１のときに、Ｙ＝Ａとなるため、Ａ値は１年分の板厚減少量（腐食量）に相当するものであり、図２の縦軸のＡ値（１年）は１年分の板厚減少量により求められたＡ値を示している。また、横軸のＡ値（〜９年）は暴露試験を９年間行ったときに、各経過年の板厚減少量に基づいて求められたＡ値を示している。図２の特性から明らかなように、暴露期間が９年の試験データによって求められたＡ値（Ａ値（〜９年））と、暴露期間が１年の試験データによって求められたＡ値（Ａ値（１年））とは相関があり、Ａ値は短期間（例えば１年間）の暴露試験により算出が可能であることが分かる。したがって、Ａ値を短期間の暴露試験により算出することができるので、同一の鋼種において複数の環境（Ｔ，ＰＷまたはＴＯＷ，Ｓａ）の下で、Ａ値を複数求めることにより、鋼種特有の係数α，β，γを求めることができる。勿論、このＡ値は実験室において求めてもよいことはいうまでもない。
【００３４】
図３はＳＭＡの実暴露の試験データ（ＳＭＡの９年間等）により求めたＡ値とＢ値との相関関係を示した特性図ある。この特性図からＢ値はＡ値で回帰することができ、それを表したのが上記の（３）式である。
【００３５】
上記のＡ値はその環境での鋼材自体の耐食性を示す係数であるが、腐食環境にも依存する性質をもっており、腐食環境の影響はＡ値にも含まれている。Ａ値が小さい範囲では、安定さびが形成されるためＢ値は一定値となるが、腐食環境の厳しさがある量を超えると（例えばＡ＝０．０３程度）、さびが剥がれやすくなり安定化しないためＢ値は上昇し、Ｂ＝１となって腐食曲線は直線のままとなる。但し、Ａ値が小さい範囲において、暴露年数が短い鋼材ではさびの生成量が少ないため保護効果が小さく、Ｂ値が高い値となる場合もある。このＳＭＡのデータに鋼種１（１．５Ｎｉ−０．３Ｍｏ鋼）及び鋼種２（２．５Ｎｉ−極低Ｃ鋼）の実暴露の試験データ（２年間）をプロットすると、この場合においても上記の関係式と合致していることが分かる。したがって、この関係式から鋼種に関係なく、Ｂ値をＡ値から求めることが可能であることが分かる。
前述のようにＡ値を短期間の暴露試験により算出することができるので、鋼種に応じた係数α，β，γに関しては、同一の鋼種において３条件以上の環境（Ｔ，ＰＷまたはＴＯＷ，Ｓａ）の下で、Ａ値を求めることにより各値を特定することができる。
さらに実験室での試験により係数α，β，γを決定することもできる。以下に実験室の試験に基づいてα，β，γを求める方法の一例について詳しく説明する。
【００３６】
まず、γの求め方について説明する。上記のγは、Ａ≦０．０３とＡ＞０．０３とではその値が異なる。Ａ≦０．０３のときには、例えばＳＭＡに関しては、図１４の暴露試験結果（「耐候性鋼材の橋梁への適用に関する共同研究報告書(XVIII)」（建設省土木研究所，（社）鋼材倶楽部，（社）日本橋梁建設協会、平成５年３月発行）によれば、γ＝０．４８７が得られている。また、Ｎｉ系耐候性鋼については発明者らが独自に暴露試験を行ったところ、図１５に示される特性が得られ、γ＝０．４８７が得られている。両データが一致しており、信頼性が高いものであることが分かる。よって、Ａ≦０．０３では鋼種によらずγ＝０．４９で一定値とした。
【００３７】
また、Ａ＞０．０３のときの上記のγは、恒温恒湿槽（ADVANTEC製 AGX-325）内での乾湿繰り返し試験により決定した。周期は２４ｈで、乾湿サイクルは次の６条件である。乾湿の移行時間は１ｈであり、これは１２ｈの中に含まれる。塩分付着はマイクロピペットを用いてＮａＣｌ水溶液を滴下した。滴下量は４０μＬ／ｃｍ2とし、水溶液濃度により付着塩分量を制御した。試験は最長５２週間行った。付着塩分量は、ここでは例えば０．１mdd、０．２mdd及び０．４mddの３種類について行う。（なお、この試験はα、βを求める際においても同様な条件で行われるものとする。但し、（１）〜（６）の条件である。）
（１）１３℃／９５％×１２ｈ−２０℃／６５％×１２ｈ
（２）２０℃／９５％×１２ｈ−２７℃／６５％×１２ｈ
（３）２５℃／９５％×１２ｈ−３２℃／６５％×１２ｈ
（４）２０℃／９５％×１２ｈ−３５℃／４０％×１２ｈ
（５）２５℃／９５％×１２ｈ−４０℃／４０％×１２ｈ
（６）１３℃／９５％×１２ｈ−２８℃／４０％×１２ｈ
【００３８】
なお、上記の実験の条件は、実環境において、温度が上がると相対湿度が下がる傾向があること、また、その温度範囲や湿度範囲についても上記の範囲内にあること等から設定されている。このため、この試験は所謂促進試験ではなく、実環境と同オーダーの腐食速度が得られるものである。塩分量、温度及び湿度に関して橋梁内桁の腐食をよく再現していることが既に確認されており、以下、再現腐食試験（又は再現試験）と称するものとする。
【００３９】
次に、上記の再現腐食試験のデータに基づいてγを求める方法を説明する。
図１６（Ａ）（Ｂ）（Ｃ）はＳＭＡ、鋼種１及び鋼種２の付着塩分量ごとの腐食量の時間変化を示した特性図である。同図の特性から各付着量に対応したＡ値が求められる。
【００４０】
図１７（Ａ）（Ｂ）（Ｃ）は、横軸に上記付着塩分量を、縦軸に付着塩分量に対応したＡ値をとり、その両対数をプロットした特性図であり、直線の傾きがγとなる。何れの鋼種においても、γ＝０．９の値が得られた。
【００４１】
次に、係数α及びβを求める方法について説明する。α及びβは、同一の付着塩分量でいくつかの腐食条件について連立方程式を当てることにより求められる。以下、具体的に説明する。上記の（２）式の塩分量は飛来塩分量であるが、再現腐食試験では付着塩分量であり、ともに単位はｍｄｄであるが、同じ数値でも影響度は異なる。そのため、再現試験で得られるα、βはα’、β’とおいて区別するが、α／βはα’／β’とは等しい。
【００４２】
同一付着塩分量で、かつ上記の（１）〜（６）の試験条件について、乾燥ステップ、湿潤ステップの腐食量の和が試験で得られる腐食量となる。再現試験では塩分は付着塩分で与えられるが、上記の（２）式の塩分の変数は飛来塩分であり、上述のように、ともに単位はmdd であるが、同じ数値でも影響度は異なる。上記の（２）式の飛来塩分に相当する量がわからないので、ここではSa’とおいて、計算する。なお、Ａ値が０．０３より大きいか小さいかでγが変わるので、Ａ値が０．０３超/以下でそれぞれα、βを設定する。
【００４３】
例えば、付着塩分が0.4mddで、13℃/95%×12h-20℃/65%×12hの試験の場合には、付着塩分は0.4mddでは、図１７からγ＝０．９であるので、これを用いる。
A' = (α'T+β')・Pw（T, H）・(Saγ)に塩分、温度、湿度を代入して、
乾燥ステップの腐食量  A'(13℃/95%、0.4mdd、乾燥) = (13α'+β')・Pw(13, 95)（Sa’0.9）
湿潤ステップの腐食量  A'(20℃/65%、0.4mdd、湿潤) = (20α'+β')・Pw(20, 65)（Sa’0.9）
である。
A'(13℃/95%×12h-20℃/65%×12h、0.4mdd) =（乾燥ステップの腐食量）+（湿潤ステップの腐食量）＝｛(13α'+β')・Pw(13, 95)+(20α'+β')・Pw(20, 65)}（Sa’0.9）
＝[{Pw(13, 95)*13+Pw(20, 65)*20}α'+{Pw(13, 95)+Pw(20, 65)}β']（Sa’0.9）
【００４４】
同じ0.4mddの塩分を載せた再現試験では、Saの値は不明であるが、載せた塩分量が同じであることから、少なくともSaを同じと見なすことができる。他の試験も同様に下記のように求めることができる。上記の式の下線部のA’の比から
{Pw(13,95)*13+Pw(20,65)*20}α'+{Pw(13,95)+Pw(20,65)}β'=1
{Pw(20,95)*20+Pw(27,65)*27}α'+{Pw(20,95)+Pw(27,65)}β'=1.403
{Pw(25,95)*25+Pw(32,65)*32}α'+{Pw(25,95)+Pw(32,65)}β'=1.661
{Pw(25,95)*25+Pw(40,40)*40}α'+{Pw(25,95)+Pw(40,40)}β'=0.858
のようになる。ここで腐食量として、Ａ値比（基準条件のＡ値に対する相対値）をとっている。同一付着塩分量でいくつかの試験条件について同様の方程式を立て、２元連立１次方程式を解く。数学的には２試験条件あればα'、β'が得られるが、精度を上げるために４試験条件を用いて、1 次回帰によりα'、β'を求める。得られたα'、β'を用いてＡ値比を求める。ＳＭＡについて得られたα'、β'を用い、全国４１橋試験の暴露地の飛来塩分量、年平均温度・湿度（最寄の気象庁観測所データ）を代入して、
A'(予測値)= (α'T1+β')・Pw（T1, H1）・（Sa10.9）
を求め、暴露試験から得られたＡ値との関係を図１８（Ｂ）に示す。両者の回帰式を求め、その1 次係数の比Ａ'/Ａ = 1.895 が得られた。Pw（T1, H1）・（Sa10.9）= K1
とおくと、
A'(予測値)= (α'T１+β')・Pw（T1, H1）・(Sa10.9)＝(α'T1+β') K1
A(実測値)= (αT1+β)・Pw（T1, H1）・(Sa10.9)＝(αT1+β)K1
任意のT1で成立するためには、
A'/A=α'/α=β'/β=1.895＝ｋ
【００４５】
一方、A≦0.03のときは、γ＝0.487として、同様に再現腐食試験の結果より方程式をたて、1 次回帰によりα'、β'を求める。さらに図１８（Ａ）より、A'/A = 4.658が得られた。A'/A=α'/α=β'/β=4.658
【００４６】
鋼種１及び鋼種２についても、同様に、これまでの再現腐食試験の結果から
A＞0.03のとき、γ＝0.9
A≦0.03のとき、γ＝0.487
で分類して、試験条件に毎の方程式をたて、1 次回帰によりα'、β'を求め、以下同様ににしてA'/A=α'/α=β'/βがそれぞれ求められる。
【００４７】
表２は、以上から得られた各耐候性鋼の腐食量予測式の係数、べき数を纏めたものである。
【００４８】
【表２】
2006208346.tif 000004
【００４９】
上記の（２）式に、表２に示された係数及びべき数と該当する地域の環境データを当てはめるとＡ値が求められ、そのＡ値を（３）式に当てはめてＢ値を求め、そのＡ値及びＢ値を鋼材の腐食量予測式に当てはめることにより寿命予測が可能になる。
【００５０】
本実施形態１に係る寿命予測方法の計測原理の概要が明らかになったところで、次に、本実施形態１の詳細を説明する。
【００５１】
図４は本実施形態１に係る寿命予測方法の具現化した装置を概念的に示したブロック図である。この寿命予測装置は、入力手段１０、ＰＣ本体２０及び表示装置３０から構成されており、ＰＣ本体２０は後述の各種の演算処理を行う演算装置（ＣＰＵ）２１、演算装置（ＣＰＵ）２１の演算処理を司るプログラムが格納された記憶装置２２及び後述の各種の気象データ等が格納された記憶装置２３から構成されている。また、表示装置３０は操作者が入力手段１０を操作してデータを入力する際の入力画面や、演算結果である腐食量の経時変化等を表示する。
【００５２】
図５は記憶装置２３に格納されている鋼種に対応して設定された係数（α，β，γ）の例であり、鋼種に対応した係数値が設定されている。この係数値は上述のように短期間の暴露試験や実験室での試験によって求めておくものとする。
【００５３】
図６は記憶装置２３に格納されている気象データの例であり、各観測地ごとに平均気温、最高気温、最低気温、平均湿度及び平均風力とそれらの標準偏差σが格納されている。
【００５４】
図７は記憶装置２３に格納されている各地域ごとの飛来塩分量のデータの例であり、ここでは、例えば日本海（北）、日本海（南）、太平洋側、瀬戸内海、及び沖縄における飛来塩分量が示されている。
【００５５】
図８は図４の寿命予測装置の演算過程を示したフローチャートである。以下、このフローチャートに従って説明する。
（Ｓ１１）操作者が入力装置１０を操作することにより鋼種を入力すると、演算装置２１はその鋼種を取り込む。この鋼種の入力に際しては、例えば鋼種の名称を直接入力させたり、或いは図５に示される鋼種を表示装置３０の画面にリストしてそれを選択させるようにしてもよい。
（Ｓ１２）演算装置２１は、入力された鋼種に基づいて、記憶装置２３に格納されている係数（α，β，γ）を検索する。
（Ｓ１３）演算装置２１は、操作者が入力装置１０を操作して例えば橋梁架設地点の入力を促すような入力画面を表示装置３０に表示させる。操作者が入力装置１０を操作して橋梁架設地点を入力するとそれを取り込む。なお、橋梁架設地点の入力に際しては、住所又は緯度・経度を入力するものとする。
（Ｓ１４）演算装置２１は、入力された橋梁架設地点に基づいて記憶装置２３に格納された気象データから、その地点又はその地点に近接した位置にある観測地の平均温度及び湿度を検索して求める。
【００５６】
図９は表示装置３０の入力画面の例を示したものであり（但し、この段階では飛来塩分量等は表示されていない）、「橋梁架設地点の近隣の気象観測値を使用する」という操作ボタンをクリックすることにより、該当する観測地の平均温度及び平均湿度が検索されて表示される。なお、橋梁架設地点の近隣の気象観測値を使用せずに、平均温度及び平均湿度を手入力することも可能である。
【００５７】
（Ｓ１５）演算装置２１は、例えば図９の「飛来塩分量の推定へ」という操作ボタンがクリックされると、表示装置３０に図１０の入力画面を表示させて必要な情報の入力を促す。この例では、日本海（北）、日本海（南）、太平洋側、瀬戸内海及び沖縄の各地域の名称がリストされ、何れかが選択される。
（Ｓ１６）また、「橋梁架設地点の海岸からの距離」（Ｌ）が入力され、「推定実行」の操作ボタンがクリックされることにより、選択された地点の飛来塩分量Ｓ0が求められ、更に、飛来塩分量Ｓ0と離岸距離Ｌとにより、Ｓａ＝Ｓ0・Ｌ-0.6という演算式によって、橋梁架設地点の飛来塩分量を求めて表示装置３０に表示させる。図１０の「飛来塩分量確定」の操作ボタンがクリックされると、図９の入力画面に戻って、飛来塩分量及び選択された地域がそれぞれ表示される。
（Ｓ１７）図９の入力画面の「腐食量推定計算実行」の操作ボタンをクリックすることにより、演算装置２１は上記の（２）式、（３）式及び（１）式を演算して腐食量Ｙを求める。
（Ｓ１８）演算装置２１は腐食量Ｙを求めると、それを表示装置３０に表させる。
【００５８】
図１１は表示装置３０に表示された腐食量Ｙの経年変化を示した図であり、例えば１００年で腐食量Ｙが０．５mm以下であれば、ここで入力された鋼種は、橋梁架設地点において使用に耐えられることを意味する。
【００５９】
図１２は上述の腐食量予測式に基づいて予測した100年後の腐食量予測値と、1,3,5,7,9年の値から直接Y=AXBにより回帰して求めた100年後の値とを比較した特性図である。図１２示されるように、正規確率紙で良い直線性を示し、正規分布に従うことがわかる。これによりもとめた標準偏差σはσ＝0.301（≒0.3）となった。したがって、直接回帰した値に対する本腐食量予測式の予測値の精度として、予測値に付帯して表すことができる。図１１においては、図１２の予測精度を上記の表示装置３０に予測量とともに表示しており、本実施形態１による計算方法が計算精度が高いことをその表示により知ることができる。
【００６０】
図１３は暴露データによる腐食量（点）と上記の計算による腐食量Ｙ（線）とを対比させたものであり、両者が良く一致しており、本実施形態１による計算方法が計算精度が高いことが分かる。
【００６１】
実施形態２．
なお、上記の実施形態１においては飛来塩分量を求めるのに構造物の部位に関係なく求める例について説明したが、構造物の部位によっては飛来塩分が付着しにくい部位もあるので、構造物の部位に応じた部位係数を予め設定しておいて、その部位係数を上記の飛来塩分量に乗算して「飛来塩分量」を求めるようにしてもよい。
さらに、構造物の設置場所又はその周辺の温度（Ｔ）、相対湿度（Ｈ）のうちの少なくとも一方に関しても、構造物の部位に応じた部位係数をそれぞれ予め設定しておいて、その部位係数を上記の温度、相対湿度にそれぞれ乗算して、「温度（Ｔ）」、「相対湿度（Ｈ）」を求めるようにしてもよい。
また、構造物の架設地点における地域的特徴、例えば海浜地域、内陸平地、山岳等の地域的特徴を飛来塩分量の算出に反映させて飛来塩分量を高精度に求めるようにすると、腐食量を更に高精度に求めることができる。
また、上記の実施形態１においては、橋梁架設地点の飛来塩分量を求めるのに、日本海（北）、日本海（南）、太平洋側、瀬戸内海及び沖縄の何れを選択する例について説明したが、橋梁架設地点の住所又は緯度・経度に基づいて該当する地域を自動的に選択するようにしてもよい。
【００６２】
実施形態３．
また、上述の実施形態１においては、Ａ値を求めるのに、何れも暴露試験又は実験室での実験を１年程度行う例について説明したが、それよりも期間が短い期間の試験結果であってもよく、例えば数ケ月程度の試験結果に外挿法を適用することにより１年間分の腐食量を求めることにより（予測することにより）Ａ値を求めることができる。
【００６３】
実施形態４．
上記のようにして求められた腐食量（寿命データ）を鋼材に添付することにより、実構造の設計の際に利用に供することが可能になっている。なお、このデータの添付とは紙などだけでなく、電子データとして管理する状態も含むものである。また、鉄骨構造物の設計に際して、このような寿命予測がなされた鋼材を適宜選択して設計することにより構造物の寿命についても適切に予測することが可能になっている。
【００６４】
実施形態５．
なお、図４に示された装置を携帯型パソコンにより実現したり、或いはネットワーク化して、顧客が使用しようとしている鋼材の寿命予測を行って、それを表示装置３０に表示させて顧客に提示し、顧客の購入を促すことも可能である。
【図面の簡単な説明】
【００６５】
【図１】期間Ｘと鋼材の腐食量Ｙとの関係を示した特性図。
【図２】ＳＭＡ実暴露の試験データにより求めたＡ値の特性図。
【図３】ＳＭＡ実暴露の試験データ（９年間）により求めたＡ値とＢ値との相関関係を示した特性図。
【図４】本発明の実施形態１に係る寿命予測方法の具現化した装置を概念的に示したブロック図。
【図５】図４の記憶装置に格納されている鋼種に対応して設定された係数（α，β，γ）の例を示した図。
【図６】図４の記憶装置に格納されている気象データの例を示した図である。
【図７】図４の記憶装置に格納されている各地域ごとの飛来塩分量のデータの例を示した図。
【図８】図４の寿命予測装置の演算過程を示したフローチャート。
【図９】図４の表示装置の入力画面の例を示した図。
【図１０】飛来塩分量を求めるための入力画面の例を示した図。
【図１１】図４の表示装置に表示された腐食量Ｙの経年変化を示した図。
【図１２】本実施形態１による１００年後の腐食量予測値と、Ｙ＝ＡＸBにより１００年後の回帰をしたときの値の差を分布を示した図。
【図１３】暴露データによる腐食量（点）と上記の計算による腐食量Ｙ（線）とを対比させた特性図。
【図１４】ＳＭＡの実暴露試験の試験データを示した特性図。
【図１５】Ｎｉ系高耐候性鋼の試験（実験室）の試験データを示した特性図。
【図１６】ＳＭＡ、鋼種１及び鋼種２の付着塩分量ごとの時間変化を示した特性図。
【図１７】横軸に上記付着塩分量を、縦軸に付着塩分量に対応したＡ値をとり、その両対数をプロットした特性図。
【図１８】Ａ’値と暴露試験から得られたＡ値との関係を示した特性図。
【図１】
2006208346.tif 000005
【図２】
2006208346.tif 000006
【図３】
2006208346.tif 000007
【図４】
2006208346.tif 000008
【図５】
2006208346.tif 000009
【図６】
2006208346.tif 000010
【図７】
2006208346.tif 000011
【図８】
2006208346.tif 000012
【図９】
2006208346.tif 000013
【図１０】
2006208346.tif 000014
【図１１】
2006208346.tif 000015
【図１２】
2006208346.tif 000016
【図１３】
2006208346.tif 000017
【図１４】
2006208346.tif 000018
【図１５】
2006208346.tif 000019
【図１６】
2006208346.tif 000020
【図１７】
2006208346.tif 000021
【図１８】
2006208346.tif 000022