(57)【要約】      （修正有）【課題】水路トンネルの内側の覆工コンクリートの応力状態、変形状態を高い精度で推定し、覆工コンクリートの寿命予測を高い精度で行う外力推定方法を提供する。【解決手段】実測された実測線距離と、共役勾配法により算出される仮想測線距離との残差を評価関数により算出し、前記評価関数から算出される残差が最小となる仮想測線距離を満足する覆工コンクリートに作用する外力を算出し、寿命予測を行うことを特徴とする水路トンネルに作用する外力推定方法。【選択図】図１
2011043339.tif 000002
【特許請求の範囲】
【請求項１】
水路トンネルの内面に覆工されたコンクリート及び水路トンネルの地山の初期物性値を実測または推定により設定し、
前記水路トンネルの内空変位を計測し、
前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する解析モデルを作成し、
実測された実測線距離と、共役勾配法により算出される仮想測線距離との残差を評価関数とし、
前記評価関数から算出される残差が最小となる作用外力に基づき、前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出し、
算出された前記コンクリートの圧縮荷重、せん断荷重、曲げモーメントと、コンクリートの許容圧縮荷重、許容せん断荷重、許容曲げモーメントを比較し、前記コンクリートの寿命予測を行う
ことを特徴とする水路トンネルに作用する外力推定方法。
【請求項２】
前記評価関数が、前記仮想測線距離と前記実測線距離との二乗和であり、共役勾配法、準ニュートン法により前記評価関数を最小にする請求項１記載の水路トンネルに作用する外力推定方法。
【請求項３】
複数の前記実測線距離がそれぞれ連続する測線である請求項１または２記載の水路トンネルに作用する外力推定方法。
【請求項４】
前記解析モデルが、前記コンクリートを平面はり、前記水路トンネルの地山をばね支承モデルと近似した請求項１乃至３記載の水路トンネルに作用する外力推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、水路トンネルの覆工コンクリートに作用する外力を推定する方法に関する発明であり、特に、水路トンネルの抜水点検時における内空変位測定の測定値を利用し、水路トンネルの内面に覆工されたコンクリートに作用する外力を高い精度で推定可能な方法に関する発明である。
【背景技術】
【０００２】
従来、水路トンネルの覆工コンクリートに作用する外力を推定する方法としては、縮小模型を製作し、縮小模型の覆工コンクリートに作用する荷重及び覆工コンクリートに発生したひび割れパターンを実測することにより、実際の水路トンネルの覆工コンクリートの残存耐力を推定する方法があった。（非特許文献１）
【０００３】
しかし、上記推定方法にあっては、簡易に概略的な残留耐力を推定できるという利点を有するものの、実際に施行された水路トンネルの覆工コンクリートの形状、地山等の物性値とは異なることから、精度が高い残留耐力を推定することは困難であった。
【０００４】
また、水路トンネルの覆工コンクリートのひび割れが発生した箇所に、ひび割れを発生させる荷重分布を算出する既設トンネルの健全性診断方法が提案されている。（特許文献１）
【０００５】
しかし、上記健全性診断方法にあっては、ひび割れを発生する種々の荷重分布の中から実際の水路トンネルの覆工コンクリートに作用している荷重分布を選択するのは困難であった。また、上記健全性診断方法にあっては、ひび割れを発生している部位に作用する荷重分布を測定できるが、ひび割れを発生していない部位に作用している荷重分布については推定できず水路トンネルの覆工コンクリートの寿命予測を行うには限界があった。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００４−２９４２６１号
【非特許文献】
【０００７】
【非特許文献１】「こうえいフォーラム第９号」  日本工営出版、２００１年１月
【発明の概要】
【発明が解決しようとする課題】
【０００８】
本発明が解決しようとする課題は、水路トンネルの内側に覆工されたコンクリート内の応力状態、覆工コンクリートの変形状態を高い精度で推定し、覆工コンクリートの寿命予測を的確に行うことにある。
【課題を解決するための手段】
【０００９】
［請求項１記載の発明］
水路トンネルの内面に覆工されたコンクリート及び水路トンネルの地山の初期物性値を実測または推定により設定し、
前記水路トンネルの内空変位を計測し、
前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する解析モデルを作成し、
実測された実測線距離と、共役勾配法により算出される仮想測線距離との残差を評価関数とし、
前記評価関数から算出される残差が最小となる作用外力に基づき、前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出し、
算出された前記コンクリートの圧縮荷重、せん断荷重、曲げモーメントと、コンクリートの許容圧縮荷重、許容せん断荷重、許容曲げモーメントを比較し、前記コンクリートの寿命予測を行う
ことを特徴とする水路トンネルに作用する外力推定方法。
【００１０】
［作用効果］
発明者の知見によれば、水路トンネルの内空変位を計測すれば、覆工コンクリートに発生した変形状態、すなわち歪み、たわみを精度良く推定できる。
【００１１】
本発明は、内空変位により実測された実測線距離と、仮想測線距離との残差を評価関数とし、前記評価関数の残差を最小とする仮想測線距離を満足する覆工コンクリートに作用する外力を共役勾配法、準ニュートン法により求めることにより、容易に実測できる内空変位を用い実際に覆工コンクリートの各箇所に発生する歪み、たわみを高い精度で推定することができる。
すなわち、前述した内空変位により、覆工コンクリートの各箇所に発生する歪み、たわみを計測することは論理的には可能である。しかし、内空変位により、覆工コンクリートの測定には、時間的、経済的、空間的な制限があることから、覆工コンクリートの各箇所に発生する歪み、たわみを実測することは現実的には困難である。
また、本発明にあっては、覆工コンクリートからコアを採取し、前記覆工コンクリートの初期物性値、すなわち弾性係数、断面二次モーメント、断面積、反発係数を実測し、それらの初期物性値に基づき共役勾配法、準ニュートン法を用いて前記評価関数を最小にする覆工コンクリートに作用する外力を推定し、実際に覆工コンクリートの各箇所に発生する歪み、たわみを高い精度で推定することができる。
【００１２】
本発明は、前述した評価関数により算出される残差を最小とする作用外力を解析モデルに代入し、覆工コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントを算出することから、覆工コンクリートの各箇所に作用する圧縮荷重、せん断荷重、曲げモーメントを高い精度で求めることができ、覆工コンクリートに使用されたコンクリートの許容圧縮荷重、許容せん断荷重、許容曲げモーメントを比較・検討し、覆工コンクリートの寿命予測を的確に行うことができる。
【００１３】
［請求項２記載の発明］
前記評価関数が、前記仮想測線距離と前記実測線距離との残差二乗和であり、共役勾配法、準ニュートン法により前記評価関数を最小にする請求項１記載の水路トンネルに作用する外力推定方法。
【００１４】
［作用効果］
請求項２記載の発明にあっては、評価関数が、仮想測線距離と実測線距離との残差二乗和からなる関数であり、共役勾配法、準ニュートン法により評価関数が最小となる覆工コンクリートに作用する外力を算出することからより高い精度で覆工コンクリートの各箇所に発生する歪み、たわみを推定できる。
【００１５】
［請求項３記載の発明］
複数の前記実測線距離がそれぞれ連続する測線である請求項１または２記載の水路トンネルに作用する外力推定方法。
【００１６】
［作用効果］
請求項３記載の発明にあっては、内空変位により実測された複数の実測線距離がそれぞれ連続し相互に関係を有し、各実測線距離が長いことからより高い精度で覆工コンクリートの各箇所に発生する歪み、たわみを推定できる。
【００１７】
［請求項４記載の発明］
前記解析モデルが、前記コンクリートを平面はり、前記水路トンネルの地山をばね支承モデルと近似し、前記平面はりの中間部と両連結とにそれぞれ対応する個別の解析モデルを作成する請求項１乃至３記載の水路トンネルに作用する外力推定方法。
【００１８】
［作用効果］
請求項４記載の発明にあっては、覆工コンクリートを平面はりに近似し、平面はりには作用する圧縮応力、せん断応力、曲げモーメントを算出することから覆工コンクリートの寿命予測を的確に行うことができる。
また、水路トンネルの地山を２方向のばね支承モデルに近似し、垂直応力、水平応力を算出することから水路トンネルの地山から覆工コンクリートに作用する垂直応力、水平応力を的確に把握することができる。
【発明の効果】
【００１９】
本発明は、水路トンネルの内側の覆工コンクリート内の応力状態、変形状態を高い精度で推定し、覆工コンクリートの寿命予測を高い精度で行うことができる。
【図面の簡単な説明】
【００２０】
【図１】本発明の実施形態を説明するフローチャートである。
【図２】覆工コンクリートに作用する外力の推定方法を説明する模式図である。
【図３】覆工コンクリートを平面はりに近似した模式図である。
【図４】水路トンネルの地山を支承ばねに近似した模式図である。
【図５】測線距離の模式図である。
【図６】断面１の実施結果である。
【図７】断面２の実施結果である。
【発明を実施するための形態】
【００２１】
以下、図面に基づいて、本発明の水路トンネルに作用する外力推定方法の実施形態について説明する。
水路トンネルに作用する外力推定方法は、図１のフローチャートに示す、水路トンネルの内面に覆工されたコンクリートからコアを採取し、前記コンクリートに作用する初期物性値を実測により、または推定により定める工程１と、前記水路トンネルのコンクリートに生じた変形による内空変位を実測する工程２と、前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する解析モデルを作成する工程３と、実測された実測線距離と、前記評価関数との残差を最小にする覆工コンクリートに作用する外力を共役勾配法、準ニュートン法により算出する工程４と、前記評価関数から算出される残差が最小となる推定外力作用時の前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する工程５と、算出された前記コンクリートの圧縮荷重、せん断荷重、曲げモーメントと、コンクリートの許容圧縮荷重、許容せん断荷重、許容曲げモーメントを比較し、前記コンクリートの寿命予測を行う工程６を有する。
なお、図１のフローチャートには、前述した工程１、工程２、工程３の順番に記載しているが、工程１〜工程３の順番は任意の順番で行うことができる。
【００２２】
［工程１：初期物性値の実測］
工程１の初期物性値とは、覆工コンクリートの弾性係数（Ｅ）、断面二次モーメント（Ｉ）、断面積（Ａ）、水路トンネルの地山２の法線方向反発係数（ｋx）、せん断方向反発係数（ｋy）をいう。
覆工コンクリートの弾性係数（Ｅ）、断面二次モーメント（Ｉ）、断面積（Ａ）は、後述する覆工コンクリートの解析モデルに用いる係数であり、弾性係数は、ボーリングコアから求める他、超音波端子法等の非破壊調査結果から現地でも推定でき、断面２次モーメント、断面積は施工図面から算出できる。
また、水路トンネル１の地山２の法線方向反発係数（ｋx）、せん断方向反発係数（ｋy）は、後述する水路トンネルの地山２の解析モデルに用いる係数であり、法線方向反発係数、せん断方向反発係数は、水路トンネル施行時に行うボーリング調査等の調査結果から推定できる。
【００２３】
［工程２：内空変位の実測］
覆工コンクリート４の内空変位を実測する。
一般的な内空変位を実測にあっては、覆工コンクリート４のスプリングライン（ＳＬ）を測定するが、スプリングライン（ＳＬ）は比較的短く測定誤差の影響を受けやすいことから、本実施形態にあっては、図２に示すとおり覆工コンクリート４の左右対称位置で、且つ、それぞれの測定箇所が連続する実測線距離（ＭＬ１〜ＭＬ４）を測定した。      
実測線距離（ＭＬ１）と実測線距離（ＭＬ４）、実測線距離（ＭＬ４）と実測線距離（ＭＬ３）、実測線距離（ＭＬ３）と実測線距離（ＭＬ２）が、それぞれ連続していることから測定箇所相互の関係が評価関数（Ｊ（ｐ））に反映できる点で好適である。
【００２４】
［工程３：解析モデルの作成］
本実施形態の図２に示した例にあっては、覆工コンクリート４と水路トンネル１の地山２との間には隙間がなく一体物と近似し、また、水路トンネル１の地山２から覆工コンクリート４に作用する外力は、図２に示すとおり覆工コンクリート４の左下部に荷重Ｐ１、左上部に荷重Ｐ２、右上部に荷重Ｐ３、右下部に荷重Ｐ４と設定した。
解析モデルについては、荷重は任意の数に近似することが可能であり、覆工コンクリート４にひび割れが生じている場合、ひび割れ位置にジョイント要素を用いることで、適切な解析モデルを作成することができる。
また、解析モデルには特に制約条件はなく、外力推定を行うトンネル断面の状態に合わせた任意の形状、拘束条件で作成でき、特殊な種々の条件に合わせ任意に設定することができる。
【００２５】
本実施形態にあっては、図３に示すとおり、覆工コンクリート４を平面はりと近似し、平面はりに作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）と歪み（μa、μb）、たわみ（νa、νb）、たわみ角（θa、θb）との関係を数１の近似式で算出し、平面はりのａ端点を数２の近似式で算出し、平面はりのｂ端点を数３の近似式で算出する。
なお、数１〜数３の近似式中のＥは覆工コンクリート４の弾性係数［ｔｆ／ｍ2］、Ｉは覆工コンクリート４の断面二次モーメント［ｍ4］、Ａは覆工コンクリート４の縦断面積［ｍ2］、ｌは平面はりの長さ［ｍ］、ｗは平面はりに作用する荷重［ｔｆ］である。なお、ｗはそれぞれの部位における荷重Ｐ１〜Ｐ４と同値である。
【００２６】
【数１】
2011043339.tif 000003
【００２７】
【数２】
2011043339.tif 000004
【００２８】
【数３】
2011043339.tif 000005
【００２９】
本実施形態にあっては、図４に示すとおり、水路トンネル１の地山２を２方向のばね支承モデルに近似し、垂直荷重（Ｘa、Ｘb）、水平荷重（Ｙa、Ｙb）と歪み（μa、μb）、たわみ（νa、νb）との関係を数４の近似式で算出する。
なお、ｋxは法線方向反発係数、ｋyはせん断方向反発係数である。
【００３０】
【数４】
2011043339.tif 000006
【００３１】
［工程４：評価関数］
本実施形態にあっては、図５に内空変位の実測で得られた実測線距離（ｌ*）と、外力（ｐ）の関数から算出される仮想測線距離（ｌ）の残差の二乗和からなる数５、数６で示す評価関数（Ｊ（ｐ））が最小となる、すなわち、内空変位の実測で得られた実測線距離（ｌ*）と仮想測線距離（ｌ）の残差の二乗和が最小となるような外力（ｐ）を共役勾配法、準ニュートン法を用いて求める。
そして、算出された外力（ｐ）を解析モデルに代入し、歪み（μa、μb）、たわみ（νa、νb）を前述した数１〜数４に代入して覆工コンクリート４に作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）、水路トンネル１の地山２に作用する垂直荷重（Ｘa、Ｘb）、水平荷重（Ｙa、Ｙb）を算出する。
なお、数１〜数３のたわみ角（θa、θb）は、たわみ（νa、νb）を微分することにより得られる。
【００３２】
【数５】
2011043339.tif 000007
【００３３】
【数６】
2011043339.tif 000008
【００３４】
本実施形態において数５、数６で示す評価関数（Ｊ（ｐ））が最小となるような外力（ｐ）をFletcher-Reevees法により求める場合に必要となる感度方程式は数７に示したとおりになる。
なお、数７のＫは覆工コンクリート４の剛性マトリックス、Γは覆工コンクリート４に作用する合成力である。
【００３５】
【数７】
2011043339.tif 000009
【００３６】
本実施形態にあっては、前述した数６で示される評価関数（Ｊ（ｐ））の最小値を求めるアルゴリズムとして共役勾配法の一種であるFletcher-Reevees法とBFGS公式による準ニュートン法を使用した。
Fletcher-Reevees法を用いる際の評価関数の勾配は数８のようになる。また、繰り返し計算のステップ幅の初期値は数９のようになり、２回目以降のステップ（β）は数１０であり、評価関数（Ｊ（ｐ）が任意に定めた収束判定定数εJ以下になるまで繰り返し計算を行う。
なお、Fletcher-Reevees法を用いた場合にあっては、他の方法に比較し上述した勾配、感度方程式の計算が容易であることから、繰り返し計算のに係る時間が短くすることができる。
【００３７】
【数８】
2011043339.tif 000010
【００３８】
【数９】
2011043339.tif 000011
【００３９】
【数１０】
2011043339.tif 000012
【００４０】
［工程５：覆工コンクリートの外力算出］
前述した方法によって、評価関数（Ｊ（ｐ））が最小値を示す、すなわち、評価関数（Ｊ（ｐ））の計算値が収束判定定数εJ以下となる外力(ｐ)の値を、数１〜数４に代入し、覆工コンクリート４に作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）及び水路トンネル１の地山２に作用する垂直荷重（Ｘa、Ｘb）、水平荷重（Ｙa、Ｙb）を算出する。
【００４１】
［工程６：覆工コンクリートの寿命予測］
算出された覆工コンクリート４に作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）と覆工コンクリート４に用いられた材料の許容圧縮荷重、許容せん断荷重、許容曲げモーメントを対比し寿命予測、あるいは次回の点検時期を定める。
【実施例】
【００４２】
本発明の水路トンネルに作用する外力推定方法を用い、内空高さ２．２４ｍ、コンクリート厚さ５０ｃｍの無圧水路トンネルの寿命予測を行った。
表１にコンクリート、水路トンネルの地山の初期条件を示す、表１の数値は１年前にコンクリート、水路トンネルの地山を実測した値である。
【００４３】
【表１】
2011043339.tif 000013
【００４４】
水路トンネルの断面１、２について、図６，７に示すとおり内空変位で測線１〜測線４を計測した。
次に、コンクリートに作用する外力をP１〜P４の分布荷重と推定し、繰り返し計算によって求めた評価関数（Ｊ（ｐ））が最小となる作用外力による測線の計算値と、測線の実測による測定値を表２に、断面１のコンクリートに作用する分布荷重及びコンクリートの変位を図６に、断面２のコンクリートに作用する分布荷重及びコンクリートの変位を図７に示す。  
【００４５】
【表２】
2011043339.tif 000014
【００４６】
表２から明らかなように、本発明の評価関数（Ｊ（ｐ））により繰返し計算を行うことにより、内空変位の各測線の測定値と、評価関数（Ｊ（ｐ））を最小とする計算から得られる各測線の計算値とが最大差異９．６２×１０-7ｍｍとなり、各測線の測定値を極めて近似する作用外力を求められることが判った。
また、図６，７のコンクリートに作用する外力、コンクリートの変位を経過観察することにより水路トンネルの内側の覆工コンクリートの寿命予測が可能であることが判った。
【産業上の利用可能性】
【００４７】
本発明は、道路トンネル、鉄道トンネル、電力系の水路トンネルの内側の覆工コンクリートに作用する応力状態、変形状態を高い精度で推定し、覆工コンクリートの寿命予測に利用できる。
【符号の説明】
【００４８】
１…水路トンネル
２…地山
３…覆工コンクリート
４…測線
【図１】
2011043339.tif 000015
【図２】
2011043339.tif 000016
【図３】
2011043339.tif 000017
【図４】
2011043339.tif 000018
【図５】
2011043339.tif 000019
【図６】
2011043339.tif 000020
【図７】
2011043339.tif 000021