JP2011043339A - 水路トンネルに作用する外力推定方法 - Google Patents
水路トンネルに作用する外力推定方法 Download PDFInfo
- Publication number
- JP2011043339A JP2011043339A JP2009189852A JP2009189852A JP2011043339A JP 2011043339 A JP2011043339 A JP 2011043339A JP 2009189852 A JP2009189852 A JP 2009189852A JP 2009189852 A JP2009189852 A JP 2009189852A JP 2011043339 A JP2011043339 A JP 2011043339A
- Authority
- JP
- Japan
- Prior art keywords
- concrete
- external force
- tunnel
- acting
- load
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
【課題】水路トンネルの内側の覆工コンクリートの応力状態、変形状態を高い精度で推定し、覆工コンクリートの寿命予測を高い精度で行う外力推定方法を提供する。
【解決手段】実測された実測線距離と、共役勾配法により算出される仮想測線距離との残差を評価関数により算出し、前記評価関数から算出される残差が最小となる仮想測線距離を満足する覆工コンクリートに作用する外力を算出し、寿命予測を行うことを特徴とする水路トンネルに作用する外力推定方法。
【選択図】図１
【解決手段】実測された実測線距離と、共役勾配法により算出される仮想測線距離との残差を評価関数により算出し、前記評価関数から算出される残差が最小となる仮想測線距離を満足する覆工コンクリートに作用する外力を算出し、寿命予測を行うことを特徴とする水路トンネルに作用する外力推定方法。
【選択図】図１
Description
本発明は、水路トンネルの覆工コンクリートに作用する外力を推定する方法に関する発明であり、特に、水路トンネルの抜水点検時における内空変位測定の測定値を利用し、水路トンネルの内面に覆工されたコンクリートに作用する外力を高い精度で推定可能な方法に関する発明である。
従来、水路トンネルの覆工コンクリートに作用する外力を推定する方法としては、縮小模型を製作し、縮小模型の覆工コンクリートに作用する荷重及び覆工コンクリートに発生したひび割れパターンを実測することにより、実際の水路トンネルの覆工コンクリートの残存耐力を推定する方法があった。（非特許文献１）
しかし、上記推定方法にあっては、簡易に概略的な残留耐力を推定できるという利点を有するものの、実際に施行された水路トンネルの覆工コンクリートの形状、地山等の物性値とは異なることから、精度が高い残留耐力を推定することは困難であった。
また、水路トンネルの覆工コンクリートのひび割れが発生した箇所に、ひび割れを発生させる荷重分布を算出する既設トンネルの健全性診断方法が提案されている。（特許文献１）
しかし、上記健全性診断方法にあっては、ひび割れを発生する種々の荷重分布の中から実際の水路トンネルの覆工コンクリートに作用している荷重分布を選択するのは困難であった。また、上記健全性診断方法にあっては、ひび割れを発生している部位に作用する荷重分布を測定できるが、ひび割れを発生していない部位に作用している荷重分布については推定できず水路トンネルの覆工コンクリートの寿命予測を行うには限界があった。
「こうえいフォーラム第９号」 日本工営出版、２００１年１月
本発明が解決しようとする課題は、水路トンネルの内側に覆工されたコンクリート内の応力状態、覆工コンクリートの変形状態を高い精度で推定し、覆工コンクリートの寿命予測を的確に行うことにある。
［請求項１記載の発明］
水路トンネルの内面に覆工されたコンクリート及び水路トンネルの地山の初期物性値を実測または推定により設定し、
前記水路トンネルの内空変位を計測し、
前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する解析モデルを作成し、
実測された実測線距離と、共役勾配法により算出される仮想測線距離との残差を評価関数とし、
前記評価関数から算出される残差が最小となる作用外力に基づき、前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出し、
算出された前記コンクリートの圧縮荷重、せん断荷重、曲げモーメントと、コンクリートの許容圧縮荷重、許容せん断荷重、許容曲げモーメントを比較し、前記コンクリートの寿命予測を行う
ことを特徴とする水路トンネルに作用する外力推定方法。
水路トンネルの内面に覆工されたコンクリート及び水路トンネルの地山の初期物性値を実測または推定により設定し、
前記水路トンネルの内空変位を計測し、
前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する解析モデルを作成し、
実測された実測線距離と、共役勾配法により算出される仮想測線距離との残差を評価関数とし、
前記評価関数から算出される残差が最小となる作用外力に基づき、前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出し、
算出された前記コンクリートの圧縮荷重、せん断荷重、曲げモーメントと、コンクリートの許容圧縮荷重、許容せん断荷重、許容曲げモーメントを比較し、前記コンクリートの寿命予測を行う
ことを特徴とする水路トンネルに作用する外力推定方法。
［作用効果］
発明者の知見によれば、水路トンネルの内空変位を計測すれば、覆工コンクリートに発生した変形状態、すなわち歪み、たわみを精度良く推定できる。
発明者の知見によれば、水路トンネルの内空変位を計測すれば、覆工コンクリートに発生した変形状態、すなわち歪み、たわみを精度良く推定できる。
本発明は、内空変位により実測された実測線距離と、仮想測線距離との残差を評価関数とし、前記評価関数の残差を最小とする仮想測線距離を満足する覆工コンクリートに作用する外力を共役勾配法、準ニュートン法により求めることにより、容易に実測できる内空変位を用い実際に覆工コンクリートの各箇所に発生する歪み、たわみを高い精度で推定することができる。
すなわち、前述した内空変位により、覆工コンクリートの各箇所に発生する歪み、たわみを計測することは論理的には可能である。しかし、内空変位により、覆工コンクリートの測定には、時間的、経済的、空間的な制限があることから、覆工コンクリートの各箇所に発生する歪み、たわみを実測することは現実的には困難である。
また、本発明にあっては、覆工コンクリートからコアを採取し、前記覆工コンクリートの初期物性値、すなわち弾性係数、断面二次モーメント、断面積、反発係数を実測し、それらの初期物性値に基づき共役勾配法、準ニュートン法を用いて前記評価関数を最小にする覆工コンクリートに作用する外力を推定し、実際に覆工コンクリートの各箇所に発生する歪み、たわみを高い精度で推定することができる。
すなわち、前述した内空変位により、覆工コンクリートの各箇所に発生する歪み、たわみを計測することは論理的には可能である。しかし、内空変位により、覆工コンクリートの測定には、時間的、経済的、空間的な制限があることから、覆工コンクリートの各箇所に発生する歪み、たわみを実測することは現実的には困難である。
また、本発明にあっては、覆工コンクリートからコアを採取し、前記覆工コンクリートの初期物性値、すなわち弾性係数、断面二次モーメント、断面積、反発係数を実測し、それらの初期物性値に基づき共役勾配法、準ニュートン法を用いて前記評価関数を最小にする覆工コンクリートに作用する外力を推定し、実際に覆工コンクリートの各箇所に発生する歪み、たわみを高い精度で推定することができる。
本発明は、前述した評価関数により算出される残差を最小とする作用外力を解析モデルに代入し、覆工コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントを算出することから、覆工コンクリートの各箇所に作用する圧縮荷重、せん断荷重、曲げモーメントを高い精度で求めることができ、覆工コンクリートに使用されたコンクリートの許容圧縮荷重、許容せん断荷重、許容曲げモーメントを比較・検討し、覆工コンクリートの寿命予測を的確に行うことができる。
［請求項２記載の発明］
前記評価関数が、前記仮想測線距離と前記実測線距離との残差二乗和であり、共役勾配法、準ニュートン法により前記評価関数を最小にする請求項１記載の水路トンネルに作用する外力推定方法。
前記評価関数が、前記仮想測線距離と前記実測線距離との残差二乗和であり、共役勾配法、準ニュートン法により前記評価関数を最小にする請求項１記載の水路トンネルに作用する外力推定方法。
［作用効果］
請求項２記載の発明にあっては、評価関数が、仮想測線距離と実測線距離との残差二乗和からなる関数であり、共役勾配法、準ニュートン法により評価関数が最小となる覆工コンクリートに作用する外力を算出することからより高い精度で覆工コンクリートの各箇所に発生する歪み、たわみを推定できる。
請求項２記載の発明にあっては、評価関数が、仮想測線距離と実測線距離との残差二乗和からなる関数であり、共役勾配法、準ニュートン法により評価関数が最小となる覆工コンクリートに作用する外力を算出することからより高い精度で覆工コンクリートの各箇所に発生する歪み、たわみを推定できる。
［請求項３記載の発明］
複数の前記実測線距離がそれぞれ連続する測線である請求項１または２記載の水路トンネルに作用する外力推定方法。
複数の前記実測線距離がそれぞれ連続する測線である請求項１または２記載の水路トンネルに作用する外力推定方法。
［作用効果］
請求項３記載の発明にあっては、内空変位により実測された複数の実測線距離がそれぞれ連続し相互に関係を有し、各実測線距離が長いことからより高い精度で覆工コンクリートの各箇所に発生する歪み、たわみを推定できる。
請求項３記載の発明にあっては、内空変位により実測された複数の実測線距離がそれぞれ連続し相互に関係を有し、各実測線距離が長いことからより高い精度で覆工コンクリートの各箇所に発生する歪み、たわみを推定できる。
［請求項４記載の発明］
前記解析モデルが、前記コンクリートを平面はり、前記水路トンネルの地山をばね支承モデルと近似し、前記平面はりの中間部と両連結とにそれぞれ対応する個別の解析モデルを作成する請求項１乃至３記載の水路トンネルに作用する外力推定方法。
前記解析モデルが、前記コンクリートを平面はり、前記水路トンネルの地山をばね支承モデルと近似し、前記平面はりの中間部と両連結とにそれぞれ対応する個別の解析モデルを作成する請求項１乃至３記載の水路トンネルに作用する外力推定方法。
［作用効果］
請求項４記載の発明にあっては、覆工コンクリートを平面はりに近似し、平面はりには作用する圧縮応力、せん断応力、曲げモーメントを算出することから覆工コンクリートの寿命予測を的確に行うことができる。
また、水路トンネルの地山を２方向のばね支承モデルに近似し、垂直応力、水平応力を算出することから水路トンネルの地山から覆工コンクリートに作用する垂直応力、水平応力を的確に把握することができる。
請求項４記載の発明にあっては、覆工コンクリートを平面はりに近似し、平面はりには作用する圧縮応力、せん断応力、曲げモーメントを算出することから覆工コンクリートの寿命予測を的確に行うことができる。
また、水路トンネルの地山を２方向のばね支承モデルに近似し、垂直応力、水平応力を算出することから水路トンネルの地山から覆工コンクリートに作用する垂直応力、水平応力を的確に把握することができる。
本発明は、水路トンネルの内側の覆工コンクリート内の応力状態、変形状態を高い精度で推定し、覆工コンクリートの寿命予測を高い精度で行うことができる。
以下、図面に基づいて、本発明の水路トンネルに作用する外力推定方法の実施形態について説明する。
水路トンネルに作用する外力推定方法は、図１のフローチャートに示す、水路トンネルの内面に覆工されたコンクリートからコアを採取し、前記コンクリートに作用する初期物性値を実測により、または推定により定める工程１と、前記水路トンネルのコンクリートに生じた変形による内空変位を実測する工程２と、前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する解析モデルを作成する工程３と、実測された実測線距離と、前記評価関数との残差を最小にする覆工コンクリートに作用する外力を共役勾配法、準ニュートン法により算出する工程４と、前記評価関数から算出される残差が最小となる推定外力作用時の前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する工程５と、算出された前記コンクリートの圧縮荷重、せん断荷重、曲げモーメントと、コンクリートの許容圧縮荷重、許容せん断荷重、許容曲げモーメントを比較し、前記コンクリートの寿命予測を行う工程６を有する。
なお、図１のフローチャートには、前述した工程１、工程２、工程３の順番に記載しているが、工程１〜工程３の順番は任意の順番で行うことができる。
水路トンネルに作用する外力推定方法は、図１のフローチャートに示す、水路トンネルの内面に覆工されたコンクリートからコアを採取し、前記コンクリートに作用する初期物性値を実測により、または推定により定める工程１と、前記水路トンネルのコンクリートに生じた変形による内空変位を実測する工程２と、前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する解析モデルを作成する工程３と、実測された実測線距離と、前記評価関数との残差を最小にする覆工コンクリートに作用する外力を共役勾配法、準ニュートン法により算出する工程４と、前記評価関数から算出される残差が最小となる推定外力作用時の前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する工程５と、算出された前記コンクリートの圧縮荷重、せん断荷重、曲げモーメントと、コンクリートの許容圧縮荷重、許容せん断荷重、許容曲げモーメントを比較し、前記コンクリートの寿命予測を行う工程６を有する。
なお、図１のフローチャートには、前述した工程１、工程２、工程３の順番に記載しているが、工程１〜工程３の順番は任意の順番で行うことができる。
［工程１：初期物性値の実測］
工程１の初期物性値とは、覆工コンクリートの弾性係数（Ｅ）、断面二次モーメント（Ｉ）、断面積（Ａ）、水路トンネルの地山２の法線方向反発係数（ｋx）、せん断方向反発係数（ｋy）をいう。
覆工コンクリートの弾性係数（Ｅ）、断面二次モーメント（Ｉ）、断面積（Ａ）は、後述する覆工コンクリートの解析モデルに用いる係数であり、弾性係数は、ボーリングコアから求める他、超音波端子法等の非破壊調査結果から現地でも推定でき、断面２次モーメント、断面積は施工図面から算出できる。
また、水路トンネル１の地山２の法線方向反発係数（ｋx）、せん断方向反発係数（ｋy）は、後述する水路トンネルの地山２の解析モデルに用いる係数であり、法線方向反発係数、せん断方向反発係数は、水路トンネル施行時に行うボーリング調査等の調査結果から推定できる。
工程１の初期物性値とは、覆工コンクリートの弾性係数（Ｅ）、断面二次モーメント（Ｉ）、断面積（Ａ）、水路トンネルの地山２の法線方向反発係数（ｋx）、せん断方向反発係数（ｋy）をいう。
覆工コンクリートの弾性係数（Ｅ）、断面二次モーメント（Ｉ）、断面積（Ａ）は、後述する覆工コンクリートの解析モデルに用いる係数であり、弾性係数は、ボーリングコアから求める他、超音波端子法等の非破壊調査結果から現地でも推定でき、断面２次モーメント、断面積は施工図面から算出できる。
また、水路トンネル１の地山２の法線方向反発係数（ｋx）、せん断方向反発係数（ｋy）は、後述する水路トンネルの地山２の解析モデルに用いる係数であり、法線方向反発係数、せん断方向反発係数は、水路トンネル施行時に行うボーリング調査等の調査結果から推定できる。
［工程２：内空変位の実測］
覆工コンクリート４の内空変位を実測する。
一般的な内空変位を実測にあっては、覆工コンクリート４のスプリングライン（ＳＬ）を測定するが、スプリングライン（ＳＬ）は比較的短く測定誤差の影響を受けやすいことから、本実施形態にあっては、図２に示すとおり覆工コンクリート４の左右対称位置で、且つ、それぞれの測定箇所が連続する実測線距離（ＭＬ１〜ＭＬ４）を測定した。
実測線距離（ＭＬ１）と実測線距離（ＭＬ４）、実測線距離（ＭＬ４）と実測線距離（ＭＬ３）、実測線距離（ＭＬ３）と実測線距離（ＭＬ２）が、それぞれ連続していることから測定箇所相互の関係が評価関数（Ｊ（ｐ））に反映できる点で好適である。
覆工コンクリート４の内空変位を実測する。
一般的な内空変位を実測にあっては、覆工コンクリート４のスプリングライン（ＳＬ）を測定するが、スプリングライン（ＳＬ）は比較的短く測定誤差の影響を受けやすいことから、本実施形態にあっては、図２に示すとおり覆工コンクリート４の左右対称位置で、且つ、それぞれの測定箇所が連続する実測線距離（ＭＬ１〜ＭＬ４）を測定した。
実測線距離（ＭＬ１）と実測線距離（ＭＬ４）、実測線距離（ＭＬ４）と実測線距離（ＭＬ３）、実測線距離（ＭＬ３）と実測線距離（ＭＬ２）が、それぞれ連続していることから測定箇所相互の関係が評価関数（Ｊ（ｐ））に反映できる点で好適である。
［工程３：解析モデルの作成］
本実施形態の図２に示した例にあっては、覆工コンクリート４と水路トンネル１の地山２との間には隙間がなく一体物と近似し、また、水路トンネル１の地山２から覆工コンクリート４に作用する外力は、図２に示すとおり覆工コンクリート４の左下部に荷重Ｐ１、左上部に荷重Ｐ２、右上部に荷重Ｐ３、右下部に荷重Ｐ４と設定した。
解析モデルについては、荷重は任意の数に近似することが可能であり、覆工コンクリート４にひび割れが生じている場合、ひび割れ位置にジョイント要素を用いることで、適切な解析モデルを作成することができる。
また、解析モデルには特に制約条件はなく、外力推定を行うトンネル断面の状態に合わせた任意の形状、拘束条件で作成でき、特殊な種々の条件に合わせ任意に設定することができる。
本実施形態の図２に示した例にあっては、覆工コンクリート４と水路トンネル１の地山２との間には隙間がなく一体物と近似し、また、水路トンネル１の地山２から覆工コンクリート４に作用する外力は、図２に示すとおり覆工コンクリート４の左下部に荷重Ｐ１、左上部に荷重Ｐ２、右上部に荷重Ｐ３、右下部に荷重Ｐ４と設定した。
解析モデルについては、荷重は任意の数に近似することが可能であり、覆工コンクリート４にひび割れが生じている場合、ひび割れ位置にジョイント要素を用いることで、適切な解析モデルを作成することができる。
また、解析モデルには特に制約条件はなく、外力推定を行うトンネル断面の状態に合わせた任意の形状、拘束条件で作成でき、特殊な種々の条件に合わせ任意に設定することができる。
本実施形態にあっては、図３に示すとおり、覆工コンクリート４を平面はりと近似し、平面はりに作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）と歪み（μa、μb）、たわみ（νa、νb）、たわみ角（θa、θb）との関係を数１の近似式で算出し、平面はりのａ端点を数２の近似式で算出し、平面はりのｂ端点を数３の近似式で算出する。
なお、数１〜数３の近似式中のＥは覆工コンクリート４の弾性係数［ｔｆ／ｍ2］、Ｉは覆工コンクリート４の断面二次モーメント［ｍ4］、Ａは覆工コンクリート４の縦断面積［ｍ2］、ｌは平面はりの長さ［ｍ］、ｗは平面はりに作用する荷重［ｔｆ］である。なお、ｗはそれぞれの部位における荷重Ｐ１〜Ｐ４と同値である。
なお、数１〜数３の近似式中のＥは覆工コンクリート４の弾性係数［ｔｆ／ｍ2］、Ｉは覆工コンクリート４の断面二次モーメント［ｍ4］、Ａは覆工コンクリート４の縦断面積［ｍ2］、ｌは平面はりの長さ［ｍ］、ｗは平面はりに作用する荷重［ｔｆ］である。なお、ｗはそれぞれの部位における荷重Ｐ１〜Ｐ４と同値である。
本実施形態にあっては、図４に示すとおり、水路トンネル１の地山２を２方向のばね支承モデルに近似し、垂直荷重（Ｘa、Ｘb）、水平荷重（Ｙa、Ｙb）と歪み（μa、μb）、たわみ（νa、νb）との関係を数４の近似式で算出する。
なお、ｋxは法線方向反発係数、ｋyはせん断方向反発係数である。
なお、ｋxは法線方向反発係数、ｋyはせん断方向反発係数である。
［工程４：評価関数］
本実施形態にあっては、図５に内空変位の実測で得られた実測線距離（ｌ*）と、外力（ｐ）の関数から算出される仮想測線距離（ｌ）の残差の二乗和からなる数５、数６で示す評価関数（Ｊ（ｐ））が最小となる、すなわち、内空変位の実測で得られた実測線距離（ｌ*）と仮想測線距離（ｌ）の残差の二乗和が最小となるような外力（ｐ）を共役勾配法、準ニュートン法を用いて求める。
そして、算出された外力（ｐ）を解析モデルに代入し、歪み（μa、μb）、たわみ（νa、νb）を前述した数１〜数４に代入して覆工コンクリート４に作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）、水路トンネル１の地山２に作用する垂直荷重（Ｘa、Ｘb）、水平荷重（Ｙa、Ｙb）を算出する。
なお、数１〜数３のたわみ角（θa、θb）は、たわみ（νa、νb）を微分することにより得られる。
本実施形態にあっては、図５に内空変位の実測で得られた実測線距離（ｌ*）と、外力（ｐ）の関数から算出される仮想測線距離（ｌ）の残差の二乗和からなる数５、数６で示す評価関数（Ｊ（ｐ））が最小となる、すなわち、内空変位の実測で得られた実測線距離（ｌ*）と仮想測線距離（ｌ）の残差の二乗和が最小となるような外力（ｐ）を共役勾配法、準ニュートン法を用いて求める。
そして、算出された外力（ｐ）を解析モデルに代入し、歪み（μa、μb）、たわみ（νa、νb）を前述した数１〜数４に代入して覆工コンクリート４に作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）、水路トンネル１の地山２に作用する垂直荷重（Ｘa、Ｘb）、水平荷重（Ｙa、Ｙb）を算出する。
なお、数１〜数３のたわみ角（θa、θb）は、たわみ（νa、νb）を微分することにより得られる。
本実施形態において数５、数６で示す評価関数（Ｊ（ｐ））が最小となるような外力（ｐ）をFletcher-Reevees法により求める場合に必要となる感度方程式は数７に示したとおりになる。
なお、数７のＫは覆工コンクリート４の剛性マトリックス、Γは覆工コンクリート４に作用する合成力である。
なお、数７のＫは覆工コンクリート４の剛性マトリックス、Γは覆工コンクリート４に作用する合成力である。
本実施形態にあっては、前述した数６で示される評価関数（Ｊ（ｐ））の最小値を求めるアルゴリズムとして共役勾配法の一種であるFletcher-Reevees法とBFGS公式による準ニュートン法を使用した。
Fletcher-Reevees法を用いる際の評価関数の勾配は数８のようになる。また、繰り返し計算のステップ幅の初期値は数９のようになり、２回目以降のステップ（β）は数１０であり、評価関数（Ｊ（ｐ）が任意に定めた収束判定定数εJ以下になるまで繰り返し計算を行う。
なお、Fletcher-Reevees法を用いた場合にあっては、他の方法に比較し上述した勾配、感度方程式の計算が容易であることから、繰り返し計算のに係る時間が短くすることができる。
Fletcher-Reevees法を用いる際の評価関数の勾配は数８のようになる。また、繰り返し計算のステップ幅の初期値は数９のようになり、２回目以降のステップ（β）は数１０であり、評価関数（Ｊ（ｐ）が任意に定めた収束判定定数εJ以下になるまで繰り返し計算を行う。
なお、Fletcher-Reevees法を用いた場合にあっては、他の方法に比較し上述した勾配、感度方程式の計算が容易であることから、繰り返し計算のに係る時間が短くすることができる。
［工程５：覆工コンクリートの外力算出］
前述した方法によって、評価関数（Ｊ（ｐ））が最小値を示す、すなわち、評価関数（Ｊ（ｐ））の計算値が収束判定定数εJ以下となる外力(ｐ)の値を、数１〜数４に代入し、覆工コンクリート４に作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）及び水路トンネル１の地山２に作用する垂直荷重（Ｘa、Ｘb）、水平荷重（Ｙa、Ｙb）を算出する。
前述した方法によって、評価関数（Ｊ（ｐ））が最小値を示す、すなわち、評価関数（Ｊ（ｐ））の計算値が収束判定定数εJ以下となる外力(ｐ)の値を、数１〜数４に代入し、覆工コンクリート４に作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）及び水路トンネル１の地山２に作用する垂直荷重（Ｘa、Ｘb）、水平荷重（Ｙa、Ｙb）を算出する。
［工程６：覆工コンクリートの寿命予測］
算出された覆工コンクリート４に作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）と覆工コンクリート４に用いられた材料の許容圧縮荷重、許容せん断荷重、許容曲げモーメントを対比し寿命予測、あるいは次回の点検時期を定める。
算出された覆工コンクリート４に作用する圧縮荷重（Ｎa、Ｎb）、せん断荷重（Ｑa、Ｑb）、曲げモーメント（Ｍa、Ｍb）と覆工コンクリート４に用いられた材料の許容圧縮荷重、許容せん断荷重、許容曲げモーメントを対比し寿命予測、あるいは次回の点検時期を定める。
本発明の水路トンネルに作用する外力推定方法を用い、内空高さ２．２４ｍ、コンクリート厚さ５０ｃｍの無圧水路トンネルの寿命予測を行った。
表１にコンクリート、水路トンネルの地山の初期条件を示す、表１の数値は１年前にコンクリート、水路トンネルの地山を実測した値である。
表１にコンクリート、水路トンネルの地山の初期条件を示す、表１の数値は１年前にコンクリート、水路トンネルの地山を実測した値である。
水路トンネルの断面１、２について、図６，７に示すとおり内空変位で測線１〜測線４を計測した。
次に、コンクリートに作用する外力をP１〜P４の分布荷重と推定し、繰り返し計算によって求めた評価関数（Ｊ（ｐ））が最小となる作用外力による測線の計算値と、測線の実測による測定値を表２に、断面１のコンクリートに作用する分布荷重及びコンクリートの変位を図６に、断面２のコンクリートに作用する分布荷重及びコンクリートの変位を図７に示す。
次に、コンクリートに作用する外力をP１〜P４の分布荷重と推定し、繰り返し計算によって求めた評価関数（Ｊ（ｐ））が最小となる作用外力による測線の計算値と、測線の実測による測定値を表２に、断面１のコンクリートに作用する分布荷重及びコンクリートの変位を図６に、断面２のコンクリートに作用する分布荷重及びコンクリートの変位を図７に示す。
表２から明らかなように、本発明の評価関数（Ｊ（ｐ））により繰返し計算を行うことにより、内空変位の各測線の測定値と、評価関数（Ｊ（ｐ））を最小とする計算から得られる各測線の計算値とが最大差異９．６２×１０-7ｍｍとなり、各測線の測定値を極めて近似する作用外力を求められることが判った。
また、図６，７のコンクリートに作用する外力、コンクリートの変位を経過観察することにより水路トンネルの内側の覆工コンクリートの寿命予測が可能であることが判った。
また、図６，７のコンクリートに作用する外力、コンクリートの変位を経過観察することにより水路トンネルの内側の覆工コンクリートの寿命予測が可能であることが判った。
本発明は、道路トンネル、鉄道トンネル、電力系の水路トンネルの内側の覆工コンクリートに作用する応力状態、変形状態を高い精度で推定し、覆工コンクリートの寿命予測に利用できる。
１…水路トンネル
２…地山
３…覆工コンクリート
４…測線
２…地山
３…覆工コンクリート
４…測線
Claims (4)
- 水路トンネルの内面に覆工されたコンクリート及び水路トンネルの地山の初期物性値を実測または推定により設定し、
前記水路トンネルの内空変位を計測し、
前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出する解析モデルを作成し、
実測された実測線距離と、共役勾配法により算出される仮想測線距離との残差を評価関数とし、
前記評価関数から算出される残差が最小となる作用外力に基づき、前記コンクリートに作用する圧縮荷重、せん断荷重、曲げモーメントの値を算出し、
算出された前記コンクリートの圧縮荷重、せん断荷重、曲げモーメントと、コンクリートの許容圧縮荷重、許容せん断荷重、許容曲げモーメントを比較し、前記コンクリートの寿命予測を行う
ことを特徴とする水路トンネルに作用する外力推定方法。 - 前記評価関数が、前記仮想測線距離と前記実測線距離との二乗和であり、共役勾配法、準ニュートン法により前記評価関数を最小にする請求項１記載の水路トンネルに作用する外力推定方法。
- 複数の前記実測線距離がそれぞれ連続する測線である請求項１または２記載の水路トンネルに作用する外力推定方法。
- 前記解析モデルが、前記コンクリートを平面はり、前記水路トンネルの地山をばね支承モデルと近似した請求項１乃至３記載の水路トンネルに作用する外力推定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2009189852A JP2011043339A (ja) | 2009-08-19 | 2009-08-19 | 水路トンネルに作用する外力推定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2009189852A JP2011043339A (ja) | 2009-08-19 | 2009-08-19 | 水路トンネルに作用する外力推定方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2011043339A true JP2011043339A (ja) | 2011-03-03 |
Family
ID=43830886
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2009189852A Pending JP2011043339A (ja) | 2009-08-19 | 2009-08-19 | 水路トンネルに作用する外力推定方法 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2011043339A (ja) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016194460A (ja) * | 2015-04-01 | 2016-11-17 | 富士通株式会社 | 応力解析方法、応力解析装置及び応力解析プログラム |
| JP2017096684A (ja) * | 2015-11-19 | 2017-06-01 | 日本工営株式会社 | 粗石コンクリートの強度評価方法及び粗石コンクリートの強度評価プログラム |
| CN107577836A (zh) * | 2017-07-24 | 2018-01-12 | 同济大学 | 一种软土地层中隧道上方覆土压力的确定方法 |
| CN109989786A (zh) * | 2019-05-21 | 2019-07-09 | 东华理工大学 | 废石-尾砂高浓度充填料浆屈服应力预测模型的建立方法及应用 |
| CN112610277A (zh) * | 2020-12-08 | 2021-04-06 | 中铁隧道局集团有限公司 | 一种富水软弱围岩隧道地质灾害风险快速预测方法 |
| CN112818510A (zh) * | 2021-01-04 | 2021-05-18 | 浙江省大成建设集团有限公司 | 一种预测基坑开挖对邻近隧道竖向位移的精确方法 |
| CN115235366A (zh) * | 2022-07-20 | 2022-10-25 | 福州大学 | 基于内轮廓曲率变化的隧道结构安全性评估方法 |
-
2009
- 2009-08-19 JP JP2009189852A patent/JP2011043339A/ja active Pending
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2016194460A (ja) * | 2015-04-01 | 2016-11-17 | 富士通株式会社 | 応力解析方法、応力解析装置及び応力解析プログラム |
| JP2017096684A (ja) * | 2015-11-19 | 2017-06-01 | 日本工営株式会社 | 粗石コンクリートの強度評価方法及び粗石コンクリートの強度評価プログラム |
| CN107577836A (zh) * | 2017-07-24 | 2018-01-12 | 同济大学 | 一种软土地层中隧道上方覆土压力的确定方法 |
| CN109989786A (zh) * | 2019-05-21 | 2019-07-09 | 东华理工大学 | 废石-尾砂高浓度充填料浆屈服应力预测模型的建立方法及应用 |
| CN112610277A (zh) * | 2020-12-08 | 2021-04-06 | 中铁隧道局集团有限公司 | 一种富水软弱围岩隧道地质灾害风险快速预测方法 |
| CN112610277B (zh) * | 2020-12-08 | 2022-04-12 | 中铁隧道局集团有限公司 | 一种富水软弱围岩隧道地质灾害风险快速预测方法 |
| CN112818510A (zh) * | 2021-01-04 | 2021-05-18 | 浙江省大成建设集团有限公司 | 一种预测基坑开挖对邻近隧道竖向位移的精确方法 |
| CN112818510B (zh) * | 2021-01-04 | 2023-12-22 | 浙江省大成建设集团有限公司 | 一种预测基坑开挖对邻近隧道竖向位移的精确方法 |
| CN115235366A (zh) * | 2022-07-20 | 2022-10-25 | 福州大学 | 基于内轮廓曲率变化的隧道结构安全性评估方法 |
| CN115235366B (zh) * | 2022-07-20 | 2024-06-07 | 福州大学 | 基于内轮廓曲率变化的隧道结构安全性评估方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | Finite element modeling of concrete cover crack propagation due to non-uniform corrosion of reinforcement | |
| Ghorbani et al. | Full-field deformation measurement and crack mapping on confined masonry walls using digital image correlation | |
| Bagge et al. | In-situ methods to determine residual prestress forces in concrete bridges | |
| Orakcal et al. | Flexural modeling of reinforced concrete walls-experimental verification | |
| JP2011043339A (ja) | 水路トンネルに作用する外力推定方法 | |
| CN102928580B (zh) | 混凝土结构中钢筋锈蚀监测装置及方法 | |
| Alam et al. | Crack propagation and size effect in concrete using a non-local damage model | |
| Srivaranun et al. | Effect of the interaction of corrosion pits among multiple tensile rebars on the reliability of RC structures: Experimental and numerical investigation | |
| Kuntal et al. | Estimation of internal corrosion degree from observed surface cracking of concrete using mesoscale simulation with Model Predictive Control | |
| Rajaram et al. | Full‐field deformation measurements during seismic loading of masonry buildings | |
| Dias-da-Costa et al. | Curvature assessment of reinforced concrete beams using photogrammetric techniques | |
| KR101382381B1 (ko) | 표면 반발 경도계를 이용한 암석의 풍화도지수 산출방법 | |
| Ye et al. | Master S-N Curve-Based Fatigue Life Assessment of Steel Bridges Using Finite Element Model and Field Monitoring Data | |
| Wang et al. | Experimental study on the determination of welding residual stress in rib-deck weld by sharp indentation testing | |
| Tanaka et al. | Remaining fatigue life assessment of damaged RC decks-Data assimilation of multi-scale model and site inspection | |
| Kuntal et al. | Estimating corrosion levels along confined steel bars in concrete using surface crack measurements and mesoscale simulations guided by model predictive control | |
| He et al. | Performance assessment of partially corrosion-damaged RC segment incorporating the spatial variability of steel corrosion | |
| Wang et al. | Multiple laboratory characterization methods to identify the D-Load of reinforced concrete pipes based on three edge bearing tests | |
| Khodaie et al. | Parametric analyses on the initial stiffness of the SHS column base plate connections using FEM | |
| Pinto et al. | Use of ultrasound to estimate depth of surface opening cracks in concrete structures | |
| Yue et al. | Prediction of fatigue crack propagation in bulb stiffeners by experimental and numerical methods | |
| Xie et al. | Out-of-plane behavior of clay brick masonry infills contained within RC frames using 3D-Digital image correlation technique | |
| Sang et al. | Curvature-enabled analytical method for bending-induced crack location prediction | |
| Wang et al. | Experimental and numerical investigation on penetrating cracks growth of rib-to-deck welded connections in orthotropic steel bridge decks | |
| Mohammadipour | Interface fracture in masonry composites: a lattice approach |