特許 6502803 断層横断埋設管路の挙動推定方法及び断層横断埋設管路の挙動推定装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6502803

(24)【登録日】2019年3月29日

(45)【発行日】2019年4月17日

(54)【発明の名称】断層横断埋設管路の挙動推定方法及び断層横断埋設管路の挙動推定装置

(51)【国際特許分類】

   G01L 1/00 20060101AFI20190408BHJP

   F16L 1/028 20060101ALI20190408BHJP

   F16L 1/024 20060101ALI20190408BHJP

【ＦＩ】

   G01L1/00 M

   F16L1/028 Z

   F16L1/028 F

   F16L1/024

【請求項の数】6

【全頁数】20

(21)【出願番号】特願2015-177923(P2015-177923)

(22)【出願日】2015年9月9日

(65)【公開番号】特開2017-53725(P2017-53725A)

(43)【公開日】2017年3月16日

【審査請求日】2018年6月15日

(73)【特許権者】

【識別番号】000001052

【氏名又は名称】株式会社クボタ

(74)【代理人】

【識別番号】100107478

【弁理士】

【氏名又は名称】橋本 薫

(72)【発明者】

【氏名】岸 正蔵

(72)【発明者】

【氏名】小田 圭太

【審査官】 公文代 康祐

(56)【参考文献】

【文献】 特開平０７−２４８１００（ＪＰ，Ａ）

【文献】 特開２００３−１７７０６７（ＪＰ，Ａ）

【文献】 特開２００５−３５１４２９（ＪＰ，Ａ）

【文献】 米国特許出願公開第２０１１／０３０３０１２（ＵＳ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｌ １／００

Ｇ０１Ｌ ５／００

Ｆ１６Ｌ １／０２４

Ｆ１６Ｌ １／０２８

Ｇ０６Ｆ １７／５０

(57)【特許請求の範囲】

【請求項1】

断層横断埋設管路の挙動推定方法であって、
所定の継手回転ばねモデルで規定される許容屈曲角度θ_ａ、管の有効長Ｌに対して、断層変位量のうち管軸直角方向の断層変位量Ｈを吸収するのに必要な断層面から片側の最少の耐震継手数Ｎを数式〔数１〕に基づいて算出するとともに、そのときの管軸方向の屈曲範囲Ｌ_０を数式〔数２〕に基づいて算出する第１ステップと、

【数1】

【数2】

所定の相対変位δ_ｇｙを境にばね定数ｋ_１ｙ，ｋ_２ｙ（ｋ_１ｙ＞ｋ_２ｙ）で定義される管軸直角方向の地盤ばねモデルに従う数式〔数３〕に基づいて、管と地盤の相対変位ｙに対して管の受ける荷重ｐ（ｙ）を台形分布荷重として算出する第２ステップと、

【数3】

数式〔数４〕で求まる管軸方向の位置ｘでの台形分布荷重の曲げモーメントＭ（ｘ）から各継手位置の曲げモーメント分布を求め、求めた曲げモーメント分布から前記継手回転ばねモデルに従って各継手位置の管の屈曲角度θを求める第３ステップと、

【数4】

前記第３ステップで求めた屈曲角度θが許容屈曲角度θ_ａ以内に収まるか否かに基づいて屈曲性能を評価する屈曲性能評価ステップと、
を含む断層横断埋設管路の挙動推定方法。

【請求項2】

断層変位量のうち管軸方向の断層変位量の半値となる断層面での管と地盤の相対変位量Ｘ_ｇ、継手伸縮量δ、数式〔数５〕で規定され断層面を基準とする管軸方向の位置ｘでの管と地盤の相対変位量ｙ（ｘ）、数式〔数６〕で規定される管軸方向の影響範囲Ｘに対して、所定の相対変位δ_ｇを境にばね定数ｋ_１，ｋ_２（ｋ_１＞ｋ_２）で定義される管軸方向の地盤ばねモデルに対応した数式〔数７〕に基づいて軸力ｆ（ｙ）を算出する第４ステップと、

【数5】

【数6】

【数7】

数式〔数８〕に基づいて断層位置での軸力ｆ_ｍａｘを算出する第５ステップと、

【数8】

前記第５ステップで求めた軸力ｆ_ｍａｘが所定の基準値に収まるか否かに基づいて軸力を評価する軸力評価ステップと、
を備えている請求項１記載の断層横断埋設管路の挙動推定方法。

【請求項3】

前記軸力評価ステップで軸力ｆ_ｍａｘが所定の基準値を超えると評価されると、前記第３ステップで求めた屈曲角度θが所定の角度閾値θ_ｔ以下となる位置を継手伸縮量Δの大変位吸収ユニットの設置位置として設定する第６ステップと、
大変位吸収ユニットの設置間隔ｓ、管軸方向影響範囲内のユニット数Ｎ_ｇ、及び断層面から一つ目の大変位吸収ユニットまでの継手数ｎ_１として、数式〔数９〕で規定される管軸方向の影響範囲Ｘに対して、数式〔数１０〕に基づいて軸力ｆ_ｍａｘを算出する第７ステップと、

【数9】

【数10】

前記第７ステップで求めた軸力ｆ_ｍａｘが所定の基準値に収まるか否かに基づいて軸力を評価する大変位吸収ユニット対応軸力評価ステップと、
を備えている請求項２記載の断層横断埋設管路の挙動推定方法。

【請求項4】

前記軸力ｆ_ｍａｘと管の断面積Ａから軸応力σ_ａ＝ｆ_ｍａｘ／Ａを算出する軸応力算出ステップと、
前記曲げモーメントＭと管の断面係数Ｚから曲げ応力σ_ｂ＝Ｍ／Ｚを算出する曲げ応力算出ステップと、
前記軸応力算出で求めた軸応力σ_ａと前記曲げ応力算出で求めた曲げ応力σ_ｂを加算して応力σ＝σ_ａ＋σ_ｂを算出する応力算出ステップと、
前記応力算出ステップで求めた応力σが所定の耐力に収まるか否かに基づいて応力を評価する応力評価ステップと、
を備えている請求項２または３記載の断層横断埋設管路の挙動推定方法。

【請求項5】

請求項１から４の何れかに記載の断層横断埋設管路の挙動推定方法を実行する簡易解析実行ステップと、
前記簡易解析実行ステップで所定の評価が得られた後に、有限要素法を用いた構造解析方法を実行する詳細解析実行ステップと、
を備えている断層横断埋設管路の挙動推定方法。

【請求項6】

断層横断埋設管路の挙動推定装置であって、
請求項１から４の何れかに記載の断層横断埋設管路の挙動推定方法を実行する挙動推定演算部と、
前記挙動推定演算部による演算条件を設定する条件入力部と、
前記挙動推定演算部による演算結果を記憶する記憶部と、
前記記憶部に記憶された演算結果の何れかを表示する表示部と、
を備えている断層横断埋設管路の挙動推定装置。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、断層横断埋設管路の挙動推定方法及び断層横断埋設管路の挙動推定装置に関する。

【背景技術】

【0002】

図１（ａ），（ｂ）に示すように、耐震継手ダクタイル鉄管等、耐震継手を介して複数の管１が接合された管路が埋設された地盤に地盤沈下や地割れが発生した場合に、一つの継手の伸縮量や屈曲角が限界に達しても、隣の継手が挙動することで大きな地盤変位が吸収される。

【0003】

図１（ｃ）には、このような耐震継手ダクタイル鉄管の一例であるＮＳ形ダクタイル鉄管の継手部の伸縮挙動が示されている。上段は標準状態、中段は収縮状態、下段は伸長状態がそれぞれ示されている。図中、符号２は受口、符号３は挿口、符号４はゴム輪、符号５はロックリングを示している。当該耐震継手の伸縮量は管長の±１％であり、引抜耐力は３Ｄ（ｋＮ）、Ｄは呼び径ｍｍ、配管施工時の許容屈曲角度αは４°、最大許容屈曲角度βは８°である。

【0004】

特許文献１には、メカニカル継手を含む配管系の不等沈下に対する健全性の評価が可能な維持管理方法として、メカニカル継手を有する地中埋設配管路に沿って、地盤の沈下分布を地表で計測し、沈下分布から局所的な相対沈下量δｒと、その発生範囲の長さＬとを求め、発生範囲に含まれるメカニカル継手の最大曲げ角度θｍａｘをθｍａｘ≦２ａｒｃｔａｎ（２δｒ／Ｌ）として、最大曲げ角度θｍａｘとメカニカル継手の許容曲げ角度とを比較することにより、配管系の健全性を評価する維持管理方法が提案されている。

【先行技術文献】

【特許文献】

【0005】

【特許文献1】特開平７−２４８１００号公報

【発明の概要】

【発明が解決しようとする課題】

【0006】

上述した配管系の健全性を評価する維持管理方法によれば、地表で計測された地盤の沈下分布に基づいてメカニカル継手の曲げ角度を算出し、算出した曲げ角度と許容曲げ角度とを比較することにより既設の配管系の健全性を評価することができる。

【0007】

しかし、地震の活動期に入ったといわれている我が国では、断層横断埋設管路の挙動を推定して健全性を評価することの必要性が認識されながら、大掛かりなシミュレーション装置を用いて評価する方法以外に、未だそのための実用的な挙動推定方法が無い。また、既設の管路のみならずこれから敷設する計画管路に対しても、断層横断埋設管路の挙動を推定して十分な安全性を見込んで設計するために断層横断埋設管路の挙動推定方法が求められている。

【0008】

本発明の目的は、上述した問題点に鑑み、大掛かりなシミュレーション装置を用いず簡易に精度が得られる断層横断埋設管路の挙動推定方法及び断層横断埋設管路の挙動推定装置を提供する点にある。

【課題を解決するための手段】

【0009】

上述の目的を達成するため、本発明による断層横断埋設管路の挙動推定方法の第一の特徴構成は、特許請求の範囲の書類の請求項１に記載した通り、断層横断埋設管路の挙動推定方法であって、所定の継手回転ばねモデルで規定される許容屈曲角度θ_ａ、管の有効長Ｌに対して、断層変位量のうち管軸直角方向の断層変位量Ｈを吸収するのに必要な断層面から片側の最少の耐震継手数Ｎを数式〔数１〕に基づいて算出するとともに、そのときの管軸方向の屈曲範囲Ｌ_０を数式〔数２〕に基づいて算出する第１ステップと、

【数1】

【数2】

所定の相対変位δ_ｇｙを境にばね定数ｋ_１ｙ，ｋ_２ｙ（ｋ_１ｙ＞ｋ_２ｙ）で定義される管軸直角方向の地盤ばねモデルに従う数式〔数３〕に基づいて、管と地盤の相対変位ｙに対して管の受ける荷重ｐ（ｙ）を台形分布荷重として算出する第２ステップと、

【数3】

数式〔数４〕で求まる管軸方向の位置ｘでの台形分布荷重の曲げモーメントＭ（ｘ）から各継手位置の曲げモーメント分布を求め、求めた曲げモーメント分布から前記継手回転ばねモデルに従って各継手位置の管の屈曲角度θを求める第３ステップと、

【数4】

前記第３ステップで求めた屈曲角度θが許容屈曲角度θ_ａ以内に収まるか否かに基づいて屈曲性能を評価する屈曲性能評価ステップと、を含む点にある。

【0010】

第１ステップでは、管軸方向の断層変位量とは無関係に、管軸直角方向の断層変位量Ｈを吸収するのに必要な断層面から片側の最少の耐震継手数Ｎの全てが許容屈曲角度θ_ａ屈曲するとの仮定の下で、数式〔数１〕に基づいて最少の耐震継手数Ｎが求められ、そのときの管軸方向の屈曲範囲Ｌ_０が数式〔数２〕により求められる。
第２ステップでは、管と地盤の相対変位ｙが断層面を中心に線形に分布するものと仮定し、所定の管軸直角方向の地盤ばねモデルに従う数式〔数３〕により、管と地盤の相対変位ｙに対して管の受ける荷重ｐ（ｙ）が台形分布荷重として算出される。
第３ステップでは、数式〔数４〕により台形分布荷重に対する各継手位置の曲げモーメント分布が求められ、継手回転ばねモデルに従って各継手位置の管の屈曲角度θが求められる。
屈曲性能評価ステップでは、算出された屈曲角度θが許容屈曲角度θ_ａ以内に収まるか否により耐屈曲性能が評価される。即ち、管軸方向の断層変位量が異なる様々な断層変位に対して管軸直角方向の断層変位量Ｈのみで耐屈曲性能が評価される。

【0011】

同第二の特徴構成は、同請求項２に記載した通り、上述の第一の特徴構成に加えて、断層変位量のうち管軸方向の断層変位量の半値となる断層面での管と地盤の相対変位量Ｘ_ｇ、継手伸縮量δ、数式〔数５〕で規定され断層面を基準とする管軸方向の位置ｘでの管と地盤の相対変位量ｙ（ｘ）、数式〔数６〕で規定される管軸方向の影響範囲Ｘに対して、所定の相対変位δ_ｇを境にばね定数ｋ_１，ｋ_２（ｋ_１＞ｋ_２）で定義される管軸方向の地盤ばねモデルに対応した数式〔数７〕に基づいて軸力ｆ（ｙ）を算出する第４ステップと、

【数5】

【数6】

【数7】

数式〔数８〕に基づいて断層位置での軸力ｆ_ｍａｘを算出する第５ステップと、

【数8】

前記第５ステップで求めた軸力ｆ_ｍａｘが所定の基準値に収まるか否かに基づいて軸力を評価する軸力評価ステップと、を備えている点にある。

【0012】

第４ステップでは、管軸方向の地盤ばねモデルに対応した数式〔数７〕に基づいて軸力ｆ（ｙ）が算出され、第５ステップでは、断層面から管軸方向の影響範囲Ｘまで軸力ｆ（ｙ）を積分した軸力ｆ_ｍａｘが算出される。
そして、軸力評価ステップでは、軸力ｆ_ｍａｘが所定の基準値に収まるか否かに基づいて軸力が評価される。即ち、管軸直角方向の断層変位量が異なる様々な断層変位に対して管軸方向の断層変位量Ｘ_ｇのみで耐軸力性能が評価される。

【0013】

同第三の特徴構成は、同請求項３に記載した通り、上述の第二の特徴構成に加えて、前記軸力評価ステップで軸力ｆ_ｍａｘが所定の基準値を超えると評価されると、前記第３ステップで求めた屈曲角度θが所定の角度閾値θ_ｔ以下となる位置を継手伸縮量Δの大変位吸収ユニットの設置位置として設定する第６ステップと、大変位吸収ユニットの設置間隔ｓ、管軸方向影響範囲内のユニット数Ｎ_ｇ、及び断層面から一つ目の大変位吸収ユニットまでの継手数ｎ_１として、数式〔数９〕で規定される管軸方向の影響範囲Ｘに対して、数式〔数１０〕に基づいて軸力ｆ_ｍａｘを算出する第７ステップと、

【数9】

【数10】

前記第７ステップで求めた軸力ｆ_ｍａｘが所定の基準値に収まるか否かに基づいて軸力を評価する大変位吸収ユニット対応軸力評価ステップと、を備えている点にある。

【0014】

軸力評価ステップで軸力ｆ_ｍａｘが所定の基準値を超えると評価されると、第６ステップでは、第３ステップで求めた屈曲角度θが所定の角度閾値θ_ｔ以下となる位置が継手伸縮量Δの大変位吸収ユニットの設置位置として設定され、第７ステップでは、数式〔数９〕で規定される管軸方向の影響範囲Ｘに対して、数式〔数１０〕に基づいて軸力ｆ_ｍａｘが算出される。そして、大変位吸収ユニット対応軸力評価ステップでは、第７ステップで求められた軸力ｆ_ｍａｘが所定の基準値に収まるか否かに基づいて耐軸力性能が評価される。

【0015】

同第四の特徴構成は、同請求項４に記載した通り、上述の第二または第三の特徴構成に加えて、前記軸力ｆ_ｍａｘと管の断面積Ａから軸応力σ_ａ＝ｆ_ｍａｘ／Ａを算出する軸応力算出ステップと、前記曲げモーメントＭと管の断面係数Ｚから曲げ応力σ_ｂ＝Ｍ／Ｚを算出する曲げ応力算出ステップと、前記軸応力算出で求めた軸応力σ_ａと前記曲げ応力算出で求めた曲げ応力σ_ｂを加算して応力σ＝σ_ａ＋σ_ｂを算出する応力算出ステップと、前記応力算出ステップで求めた応力σが所定の耐力に収まるか否かに基づいて応力を評価する応力評価ステップと、を備えている点にある。

【0016】

応力評価ステップでは、軸力ｆ_ｍａｘと管の断面積Ａから算出された軸応力σ_ａと、曲げモーメントＭと管の断面係数Ｚから算出された曲げ応力σ_ｂとの和で求まる応力が許容値に収まるか否かが評価される。

【0017】

同第五の特徴構成は、同請求項５に記載した通り、上述の第一から第四の何れかの特徴構成を備えた断層横断埋設管路の挙動推定方法を実行する簡易解析実行ステップと、前記簡易解析実行ステップで所定の評価が得られた後に、有限要素法を用いた構造解析方法を実行する詳細解析実行ステップと、を備えている点にある。

【0018】

短時間で繰り返して実行可能な第一から第四の何れかの特徴構成を備えた断層横断埋設管路の挙動推定方法を実行して、満足の得られた解析結果に基づいて有限要素法を用いた構造解析方法を用いた解析を行なうことで、信頼性を十分に確保した解析が実現でき、時間の掛かる有限要素法を用いた構造解析の繰返し回数を大きく低減できるので、より迅速な評価作業環境が整えられるようになった。

【0019】

本発明による断層横断埋設管路の挙動推定装置の特徴構成は、同請求項６に記載した通り、断層横断埋設管路の挙動推定装置であって、上述した第一から第四の何れかの特徴構成を備えた断層横断埋設管路の挙動推定方法を実行する挙動推定演算部と、前記挙動推定演算部による演算条件を設定する条件入力部と、前記挙動推定演算部による演算結果を記憶する記憶部と、前記記憶部に記憶された演算結果の何れかを表示する表示部と、を備えている点にある。

【0020】

このような断層横断埋設管路の挙動推定装置を用いることにより、演算時間のかかる有限要素法を用いたシミュレーション演算を行なうことなく、短時間である程度の信頼性のある評価が行なえるようになり、既設管路の耐震性評価、耐震性の高い未設管路の設計が行なえるようになる。

【発明の効果】

【0021】

以上説明した通り、本発明によれば、大掛かりなシミュレーション装置を用いず簡易に精度が得られる断層横断埋設管路の挙動推定方法及び断層横断埋設管路の挙動推定装置を提供することができる。

【図面の簡単な説明】

【0022】

【図1】（ａ）は耐震継手で接合された管路の地盤沈下時の挙動説明図、（ｂ）は同地割れ時の挙動説明図、（ｃ）は耐震継手の伸縮動作説明図

【図2】断層が管路に作用する軸力及び継手屈曲角の説明図

【図3】（ａ）は継手ばね及び地盤ばねの説明図、（ｂ）〜（ｄ）は継手ばねモデルの説明図、（ｅ），（ｆ）は地盤ばねモデルの特性図

【図4】断層近傍の継手から屈曲する様子を示したシミュレーション結果の説明図

【図5】（ａ）は継手屈曲角度の挙動推定方法に用いる継手回転ばねモデルの説明図、（ｂ）は断層変位Ｈに対する管の屈曲範囲Ｌ_０の説明図管路

【図6】（ａ）は継手屈曲角度の挙動推定方法に用いる管軸直角方向の地盤ばねモデルの説明図、（ｂ）は管路と地盤の相対変位ｙの分布の説明図、（ｃ）は管路に掛かる荷重p（ｙ）の分布の説明図

【図7】（ａ）は本発明の挙動推定方法による推定結果の継手回転ばねの説明図、（ｂ）は本発明の挙動推定方法による推定結果とＦＥＭ解析結果を対比した継手屈曲角度の分布の説明図

【図8】断層交差角度が異なっても軸力−管軸方向地盤変位の曲線が一致することを示すＦＥＭ解析結果の説明図

【図9】（ａ）は断層変位により継手が圧縮される状態の説明図、（ｂ）は断層面を境に左右対称に地盤と管が相対変位することを示す説明図、（ｃ）は単位長さ当たりの軸力分布の説明図

【図10】（ａ）は管と地盤の相対変位量ｙ（ｘ）の説明図、（ｂ）は軸力を推定するために用いる管軸方向地盤ばねモデルの説明図、（ｃ）継手の屈曲角度分布の説明図

【図11】本発明の挙動推定方法による推定結果とＦＥＭ解析結果を対比した軸力特性の説明図

【図12】（ａ）は大変位対応ユニットを用いた場合の管と地盤の相対変位量ｙ（ｘ）の説明図、（ｂ）は軸力の影響範囲Ｘの説明図

【図13】（ａ），（ｂ）は軸力の影響範囲Ｘの算定方法の説明図

【図14】大変位対応ユニットを用いた場合の本発明の挙動推定方法による推定結果とＦＥＭ解析結果を対比した軸力特性の説明図

【図15】断層横断埋設管路の挙動推定装置の機能ブロック構成図

【図16】断層横断埋設管路の挙動推定装置の解析条件入力画面の説明図

【図17】断層横断埋設管路の挙動推定装置の解析結果表示画面の説明図

【図18】断層横断埋設管路の挙動推定装置の解析結果表示画面の説明図

【図19】断層横断埋設管路の挙動推定方法の手順を示すフローチャート

【発明を実施するための形態】

【0023】

以下に本発明による断層横断埋設管路の挙動推定方法及び断層横断埋設管路の挙動推定装置を、耐震継手ダクタイル鉄管の一例であるＮＳ形ダクタイル鉄管の直管を例に説明する。尚、本発明は、ＮＳ形ダクタイル鉄管に限らず、耐震継手を介して複数の管が接合される断層横断埋設管路全般に適用できる。

【0024】

図２に示すように、本発明は、断層角度φの断層変位を管軸方向と管軸直角方向に二成分に分けて「継手屈曲角」、「軸力」及び「応力」を簡易的に算定することにより管路挙動を推定し評価する方法及び装置であり、継手部を含む管路をモデル化してＦＥＭ解析を行なった結果に基づいて構成されている。ＦＥＭ解析とは、有限要素法を用いた構造解析のことである。

【0025】

以下、順に説明する。
図３（ａ）には、継手ばね及び地盤ばねモデルが示されている。解析対象となる耐震継手で接合される管路を弾性床上の梁(剛体)と見なし、実験に基づいて得られたばね定数に基づいて継手および地盤をモデル化した。

【0026】

図３（ｂ）から（ｄ）には継手ばねの特性が示され、図３（ｅ），（ｆ）には地盤ばねの特性が示されている。継手ばねのうちの回転ばねは、受口側の管の管軸と挿口側の管の管軸が屈曲した状態となって受口側の管の継手部内面に挿口側の管の継手部外面が接触して曲がりにくくなる角度(図３（ｃ）のθ_ａ)が設定されている。当該角度θ_ａを許容屈曲角度といい、本実施形態では、許容屈曲角度を一律に４．０°に設定している。

【0027】

呼び径７５〜２５０のＮＳ形ダクタイル鉄管は地震時に曲がりうる最大屈曲角度が８°と規定され、８°以内で継手性能が保持される。尚、許容屈曲角度及び最大許容屈曲角度の値は、本実施形態に記載した値に限るものではなく、呼び径により異なる値であってもよいし、耐震継手鉄管の種類により異なる値に設定されていてもよい。

【0028】

軸方向ばねは、管とゴム輪が滑り継手部が伸縮する領域、継手の抜け出し防止機構が働き伸縮が止まる領域でそれぞれ異なるばね定数が設定されている。

【0029】

地盤ばねのうちの管軸方向地盤ばねは、管と地盤の滑りを考慮して管と地盤の相対変位が限界値を超えるとばね定数が小さくなるバイリニアモデルで設定され、管軸直角方向地盤ばねは、管が地盤に対し相対的に下方へ移動する場合は地盤反力を考慮し、管が地盤に対し相対的に上方へ移動する場合は地盤の崩壊を考慮して設定されている。何れも管と地盤の相対変位が限界値を超えるとばね定数が小さくなるバイリニアモデルで設定されている。

【0030】

３次元骨組構造物非線形動的解析システム「ＤＹＮＡ２Ｅ」（伊藤忠テクノソリューションズ株式会社）を用いて、このような継手ばねおよび地盤ばねモデルに対して断層位置や管長を変えてＦＥＭ解析を行なった。

【0031】

先ず、管軸直角方向地盤変位による継手屈曲角度の評価方法について説明する。
図４に示すように、管路に生じる軸力が小さい場合、断層近傍の継手から順に継手が屈曲し、断層変位が吸収されることが判明している。

【0032】

図５（ａ），（ｂ）に示すように、ＦＥＭ解析結果から、断層の影響範囲にある継手は全て許容屈曲角度θ_ａだけ屈曲するとともに断層面にある継手は屈曲しないとの仮定の下で、上述した所定の継手回転ばねモデルで規定される許容屈曲角度θ_ａ、管の有効長Ｌに対して、断層変位量のうち管軸直角方向の断層変位量Ｈを吸収するのに必要な断層面から片側の最少の耐震継手数Ｎを数式〔数１１〕に基づいて算出するとともに、そのときの管軸方向の屈曲範囲Ｌ_０を数式〔数１２〕に基づいて算出する第１ステップを実行する。

【数11】

【数12】

【0033】

次に、図６（ｂ）に示すように、管路と地盤の相対変位ｙは断層面を中心に安全側となる線形に分布すると仮定するとともに、断層面位置の相対変位はＨ／２、断層面からＬ_０離れた位置の相対変位は０になるとの仮定の下で、図６（ａ）に示すような所定の相対変位δ_ｇｙを境にばね定数ｋ_１ｙ，ｋ_２ｙ（ｋ_１ｙ＞ｋ_２ｙ）で定義される管軸直角方向の地盤ばねモデルに従う数式〔数１３〕に基づいて、図６（ｃ）に示すような管と地盤の相対変位ｙに対して管の受ける荷重ｐ（ｙ）を台形分布荷重として算出する第２ステップを実行する。

【数13】

【0034】

次に、数式〔数１４〕で求まる管軸方向の位置ｘでの台形分布荷重の曲げモーメントＭ（ｘ）から各継手位置の曲げモーメント分布を求め、求めた曲げモーメント分布から図７（ａ）に示すような継手回転ばねモデルに従って各継手位置の管の屈曲角度θを求める第３ステップを実行する。

【数14】

【0035】

さらに、第３ステップで求めた屈曲角度θが許容屈曲角度θ_ａ以内に収まるか否かに基づいて屈曲性能を評価する屈曲性能評価ステップを実行する。その結果、管軸方向の断層変位量が異なる様々な断層変位に対して管軸直角方向の断層変位量Ｈのみで耐屈曲性能が評価できるようになる。

【0036】

図７（ｂ）には、ＦＥＭ解析の結果と、本発明による解析結果が示されており、断層近傍の継手屈曲角度は概ね一致することが判明し、簡易的に継手屈曲角度を評価する本発明による方法の有効性が証明された。

【0037】

次に、簡易的な軸力の評価方法について説明する。
図８に示すように、ＦＥＭ解析の結果、断層交差角度が異なっても、軸力−管軸方向地盤変位の曲線はほぼ一致することが判明している。
図９（ａ），（ｂ）に示すように、断層変位後に破線で示される管路の継手が圧縮され、各継手部で相対変位を吸収するとの仮定の下で、図９（ｃ）に示すように、地盤ばねモデルを用いて単位長さ当たりの軸力ｆ（ｙ）分布を求め、軸力ｆ（ｙ）を図中ハッチングで示される範囲で積分することにより、軸力ｆ_ｍａｘを求める。

【0038】

具体的に、図１０（ａ），（ｂ）に示すように、断層変位量のうち管軸方向の断層変位量の半値となる断層面での管と地盤の相対変位量Ｘ_ｇ、継手伸縮量δ、数式〔数１５〕で規定され断層面を基準とする管軸方向の位置ｘでの管と地盤の相対変位量ｙ（ｘ）、数式〔数１６〕で規定される管軸方向の影響範囲Ｘに対して、所定の相対変位δ_gを境にばね定数ｋ_１，ｋ_２（ｋ_１＞ｋ_２）で定義される管軸方向の地盤ばねモデルに対応した数式〔数１７〕に基づいて軸力ｆ（ｙ）を算出する第４ステップを実行する。

【数15】

【数16】

【数17】

【0039】

次に、数式〔数１８〕に基づいて断層位置での軸力ｆ_ｍａｘを算出する第５ステップを実行する。

【数18】

【0040】

さらに、第５ステップで求めた軸力ｆ_ｍａｘが所定の基準値に収まるか否かに基づいて軸力を評価する軸力評価ステップを実行する。所定の基準値は、３ＤｋＮ（ここに、Ｄは呼び径である。）により設定される値が好適に用いられる。即ち、管軸直角方向の断層変位量が異なる様々な断層変位に対して管軸方向の断層変位量Ｘ_ｇのみで耐軸力性能が評価される。

【0041】

図１１に示すように、断層角４５°，６０°に対するＦＥＭ解析の結果と、本発明による方法を対比すると、低軸力時には曲線はほぼ一致し軸力３ＤｋＮ付近で安全側に評価できることが証明された。

【0042】

上述した軸力評価ステップで、軸力ｆ_ｍａｘが所定の基準値を超えると評価されると、図１０（ｃ）に示すように、第３ステップで求めた屈曲角度θが所定の角度閾値θ_ｔ以下となる位置を継手伸縮量（長尺継ぎ輪の伸縮量）Δの大変位吸収ユニットの設置位置として設定する第６ステップを実行する。図１０（ｃ）では、角度閾値θ_ｔ以下となるＤ−Ｅ間に大変位吸収ユニットが設置される。

【0043】

管軸直角方向の変位の影響を受けた範囲、つまり継手が屈曲した範囲を挟むように、大変位対応ユニットを設ける必要がある。そこで、具体的に角度閾値θ_ｔは安全を考慮して、図７（ａ）に示す継手回転ばねのモーメントが非常に大きくなる角度θ_ａよりも十分に小さな角度に設定され、本実施形態では、角度θ_ａの３．２°に対して角度閾値θ_ｔを１°に設定されている。しかし、角度閾値θ_ｔは１°に制限される値ではなく、安全を考慮して適宜設定される値であればよい。

【0044】

図１２（ａ）に示すように、大変位吸収ユニットを用いることにより、軸方向の影響範囲を狭めることができ、大変位吸収ユニットの外側の管路の相対変位量を低減することができ、これにより軸力が低減するようになる。

【0045】

さらに、数式〔数１９〕で規定される管軸方向の影響範囲Ｘに対して、数式〔数２０〕に基づいて軸力ｆ_ｍａｘを算出する第７ステップを実行する。数式〔数１９〕で、ｓは大変位吸収ユニットの設置間隔、Ｎ_ｇは管軸方向影響範囲内の大変位吸収ユニットの数である。

【数19】

【0046】

数式〔数１９〕で、ｎ_１は断層面から一つ目の大変位吸収ユニットまでの継手数である。Ｎ_ｇ，ｎ_１の値は、任意の実数に対して最大の整数値に変換する天井関数で定まる値である。また、影響範囲Ｘを示す式の前半部はユニットが縮み切っており、大変位吸収ユニットＮ_ｇとＮ_ｇ＋１の間の継手が縮んでいる場合（図１３（ａ）参照）を示し、後半部はユニットＮ_ｇが縮んでいる場合（図１３（ｂ）参照）を示す。

【数20】

【0047】

ｆ_ａｂは図１２（ｂ）に示す領域ａ＋ｂの相対変位によって生じる軸力である。また、ｆ_ｂは領域ｂの相対変位によって生じる軸力である。軸力ｆ_ａｂから軸力ｆ_ｂを減算することによりｆ_ｍａｘが求まる。

【0048】

上述した第７ステップで求めた軸力ｆ_ｍａｘが所定の基準値に収まるか否かに基づいて軸力を評価する大変位吸収ユニット対応軸力評価ステップが実行される。

【0049】

図１４に示すように、ＦＥＭ解析と上述の軸力の解析結果とはほぼ一致することが判明している。

【0050】

さらに、軸力ｆ_ｍａｘと管の断面積Ａから軸応力σ_ａ＝ｆ_ｍａｘ／Ａを算出する軸応力算出ステップと、曲げモーメントＭと管の断面係数Ｚから曲げ応力σ_ｂ＝Ｍ／Ｚを算出する曲げ応力算出ステップと、軸応力算出で求めた軸応力σ_ａと曲げ応力算出で求めた曲げ応力σ_ｂを加算して応力σ＝σ_ａ＋σ_ｂを算出する応力算出ステップと、応力算出ステップで求めた応力σが所定の耐力に収まるか否かに基づいて応力を評価する応力評価ステップとが実行される。

【0051】

応力評価ステップでは、軸力ｆ_ｍａｘと管の断面積Ａから算出された軸応力σａと、曲げモーメントＭと管の断面係数Ｚから算出された曲げ応力σ_ｂとの和で求まる応力が許容値に収まるか否かが評価される。

【0052】

このような簡易的な解析を行なって管路を見直した後に、最終的にＦＥＭ解析を実行することにより、簡易解析の信頼性が担保される。大変位吸収ユニットを付加する場合には、予め上述した簡易解析で大変位吸収ユニットの間隔を定めて最終的にＦＥＭ解析を実行すればよい。

【0053】

図１５に示すように、本発明による断層横断埋設管路の挙動推定装置１００は、汎用のコンピュータで構成され、上述した断層横断埋設管路の挙動推定方法を実行する挙動推定演算部２０と、挙動推定演算部２０による演算条件を設定する条件入力部１０と、挙動推定演算部２０による演算結果を記憶する記憶部３０と、記憶部３０に記憶された演算結果の何れかを表示する表示部４０と、を備えて構成されている。条件入力部１０と表示部４０とがタッチパネル式液晶表示装置１０１を用いて実現されている。

【0054】

このような断層横断埋設管路の挙動推定装置を用いることにより、演算時間のかかる有限要素法を用いたシミュレーション演算を行なうことなく、短時間である程度の信頼性のある評価が行なえるようになり、既設管路の耐震性評価、耐震性の高い未設管路の設計が行なえるようになる。

【0055】

即ち、汎用のコンピュータを断層横断埋設管路の挙動推定装置として機能させるために、所定の継手回転ばねモデルで規定される許容屈曲角度θ_ａ、管の有効長Ｌに対して、断層変位量のうち管軸直角方向の断層変位量Ｈを吸収するのに必要な断層面から片側の最少の耐震継手数Ｎを数式〔数１１〕に基づいて算出するとともに、そのときの管軸方向の屈曲範囲Ｌ_０を数式〔数１２〕に基づいて算出する第１ステップと、所定の相対変位δ_ｇｙを境にばね定数ｋ_１ｙ，ｋ_２ｙ（ｋ_１ｙ＞ｋ_２ｙ）で定義される管軸直角方向の地盤ばねモデルに従う数式〔数１３〕に基づいて、管と地盤の相対変位ｙに対して管の受ける荷重ｐ（ｙ）を台形分布荷重として算出する第２ステップと、数式〔数１４〕で求まる管軸方向の位置ｘでの台形分布荷重の曲げモーメントＭ（ｘ）から各継手位置の曲げモーメント分布を求め、求めた曲げモーメント分布から前記継手回転ばねモデルに従って各継手位置の管の屈曲角度θを求める第３ステップと、前記第３ステップで求めた屈曲角度θが許容屈曲角度θ_ａ以内に収まるか否かに基づいて屈曲性能を評価する屈曲性能評価ステップと、をコンピュータに実行させるためのプログラムがインストールされている。

【0056】

プログラムは、また、断層変位量のうち管軸方向の断層変位量の半値となる断層面での管と地盤の相対変位量Ｘ_ｇ、継手伸縮量δ、数式〔数１５〕で規定され断層面を基準とする管軸方向の位置ｘでの管と地盤の相対変位量ｙ（ｘ）、数式〔数１６〕で規定される管軸方向の影響範囲Ｘに対して、所定の相対変位δ_gを境にばね定数ｋ_１，ｋ_２（ｋ_１＞ｋ_２）で定義される管軸方向の地盤ばねモデルに対応した数式〔数１７〕に基づいて軸力ｆ（ｙ）を算出する第４ステップと、数式〔数１８〕に基づいて断層位置での軸力ｆ_ｍａｘを算出する第５ステップと、前記第５ステップで求めた軸力ｆ_ｍａｘが所定の基準値に収まるか否かに基づいて軸力を評価する軸力評価ステップと、を備えている。

【0057】

さらに、前記軸力評価ステップで軸力ｆ_ｍａｘが所定の基準値を超えると評価されると、前記第３ステップで求めた屈曲角度θが所定の角度閾値θ_ｔ以下となる位置を継手伸縮量Δの大変位吸収ユニットの設置位置として設定する第６ステップと、数式〔数１９〕で規定される管軸方向の影響範囲Ｘに対して、数式〔数２０〕に基づいて軸力ｆ_ｍａｘを算出する第７ステップと、前記第７ステップで求めた軸力ｆ_ｍａｘが所定の基準値に収まるか否かに基づいて軸力を評価する大変位吸収ユニット対応軸力評価ステップと、を備えている。

【0058】

さらにまた、前記軸力ｆ_ｍａｘと管の断面積Ａから軸応力σ_ａ＝ｆ_ｍａｘ／Ａを算出する軸応力算出ステップと、前記曲げモーメントＭと管の断面係数Ｚから曲げ応力σ_ｂ＝Ｍ／Ｚを算出する曲げ応力算出ステップと、前記軸応力算出で求めた軸応力σ_ａと前記曲げ応力算出で求めた曲げ応力σ_ｂを加算して応力σ＝σ_ａ＋σ_ｂを算出する応力算出ステップと、前記応力算出ステップで求めた応力σが所定の耐力に収まるか否かに基づいて応力を評価する応力評価ステップと、を備えている。

【0059】

このようなプログラムは汎用の表計算ソフトに、例えばマクロ命令を用いて組み込むことができる。

【0060】

図１６には、断層横断埋設管路の挙動推定装置の入力画面が示されている。断層交差角、管の呼び径、伸縮量、管長等の管路条件や各種のばねモデルの定数が入力されると、簡易解析が実行され、図１７，図１８に示すような演算結果が表示される。この出力画面によって解析結果が良好であるか否かが一目瞭然で判断できるように構成されている。

【0061】

図１９には、断層横断埋設管路の挙動推定装置及びＦＥＭ解析装置を用いた挙動推定手順が示されている。先ず、挙動推定装置１００により定尺管に対する屈曲角度の簡易解析が行なわれ（Ｓ１）、結果が評価され（Ｓ２）、ＮＧであるなら管路を見直して（Ｓ３）、再度屈曲角度の簡易解析が行なわれる。

【0062】

屈曲角度の簡易解析結果がＯＫなら（Ｓ２）、軸力の簡易解析が行なわれ（Ｓ４）、結果が評価され（Ｓ５）、ＮＧであるなら管路を見直して（Ｓ６）、再度軸力の簡易解析が行なわれる。

【0063】

軸力の簡易解析結果がＯＫなら（Ｓ５）、応力の簡易解析が行なわれ（Ｓ７）、結果が評価され（Ｓ８）、ＮＧであるなら管路を見直して（Ｓ９）、再度ステップＳ１からの簡易解析が行なわれる。

【0064】

応力の簡易解析結果がＯＫなら（Ｓ８）、ＦＥＭ解析が行なわれ（Ｓ１０）、結果が評価され、ＮＧなら（Ｓ１１）、管路を見直して再度ＦＥＭ解析が行なわれる（Ｓ１０）。ＯＫなら解析を終了する。このように、挙動推定装置１００を用いて、短時間で繰り返して実行されるステップＳ１からステップＳ９の解析を行なった結果、満足の得られた解析結果に基づいてＦＥＭ解析を行なうことで信頼性を十分に確保した解析が実現できる。そして、最終の時間の掛かるＦＥＭ解析の繰返し回数を大きく低減できるので、より迅速な評価作業環境が整えられるようになった。

【0065】

図１９に示したフローチャートでは、ステップＳ１，Ｓ４，Ｓ７の各解析の何れでも満足できる評価結果が得られた後にステップＳ１０のＦＥＭ解析を実行する例が示されているが、ステップＳ１，Ｓ４，Ｓ７の何れか一つの簡易解析で満足できる評価結果が得られた後にステップＳ１０のＦＥＭ解析を実行するように構成されていてもよい。つまり、本発明は、屈曲角度、軸力の何れか一つを評価するものであってもよいし、屈曲角度、軸力の何れか一つを評価した後にＦＥＭ解析を実行するように構成されていてもよい。

【0066】

上述した説明は、断層横断埋設管路の挙動推定方法及び断層横断埋設管路の挙動推定装置の一実施形態の説明であり、該記載により本発明の範囲が限定されるものではなく、目的変量及び目的変量のカテゴリの係数が上述した種類及び数値に限るものではなく、本発明の作用効果が奏される範囲で適宜変更設計可能であることはいうまでもない。

【符号の説明】

【0067】

１００：断層横断埋設管路の挙動推定装置
１０：条件入力部
２０：挙動推定演算部
３０：記憶部
４０：表示部

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】

【図19】