JP6041151B2 - 粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体 - Google Patents

粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体 Download PDF

Info

Publication number
JP6041151B2
JP6041151B2 JP2013124669A JP2013124669A JP6041151B2 JP 6041151 B2 JP6041151 B2 JP 6041151B2 JP 2013124669 A JP2013124669 A JP 2013124669A JP 2013124669 A JP2013124669 A JP 2013124669A JP 6041151 B2 JP6041151 B2 JP 6041151B2
Authority
JP
Japan
Prior art keywords
subsidence
viscous fluid
settlement
velocity
calculation
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.)
Active
Application number
JP2013124669A
Other languages
English (en)
Other versions
JP2015001399A (ja
Inventor
石井　卓
卓 石井
英行 櫻井
英行 櫻井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimizu Corp
Original Assignee
Shimizu Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shimizu Corp filed Critical Shimizu Corp
Priority to JP2013124669A priority Critical patent/JP6041151B2/ja
Publication of JP2015001399A publication Critical patent/JP2015001399A/ja
Application granted granted Critical
Publication of JP6041151B2 publication Critical patent/JP6041151B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Description

本発明は、粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびこのプログラムを格納したコンピュータ読み取り可能な記録媒体に関し、特に、ベントナイト系粘土材料などの湿潤粘土系の緩衝材に取り囲まれた放射性廃棄体のような物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体に関するものである。
従来、高レベル放射性廃棄物の地層処分では、図１０に示すように、緩衝材と呼ばれる厚さ０．３ｍ〜１ｍのベントナイト系粘土材料１で廃棄体２周囲を取り囲み、周囲からの地圧を緩衝させるとともに、地下水の浸入を抑制し、廃棄体２からの放射性物質の漏えいを抑止することが考えられている（例えば、非特許文献１を参照）。
廃棄体の周囲を０．３ｍ〜１ｍの厚さのベントナイト系材料で取り囲むことのメリットは、この材料が有する下記の特長にある。
（１）実現可能なエネルギーで達成できる高密度なベントナイト系材料の透水係数は１Ｅ−１０〜１Ｅ−１３ｍ／ｓであり、一般的な地盤材料の透水係数（１Ｅ−５〜１Ｅ−８ｍ／ｓ）に比べて難透水性に優れている。ここで、符号Ｅは１０のべき乗を意味している。例えば「Ｅ−５」は「１０のマイナス５乗」を意味し、「１Ｅ−５」は「０．００００１」を意味している。
（２）ベントナイトはコンクリートに比べて柔軟性があるため、ひび割れが発生しにくく、かつ、存在していたひび割れは、周囲から地下水が浸入していた時点で、ベントナイトの有する吸水膨潤挙動によりシールされるので、卓越する水みちにならない。
（３）ベントナイトは天然に存在する粘土材料であることから、数万年以上の長期間を経ても材料劣化が生じにくい。
次に、本発明の予測対象となる現象について説明する。
図１１に示すように、放射性廃棄物を格納した廃棄体２は、その周囲を取り囲んでいる粘土材料１よりも大きな密度を有しているので、長期間の間に次第に沈下して位置が変化し、やがてはバリア材（粘土材料）の厚さが小さく変化するので、放射性物質の移流拡散を抑制する能力が減じることになる。したがって、将来の安全性を評価する上で廃棄体２の沈下を正しく予測することが求められていた。
こうした現象に関連し、粘土の長期沈下挙動については、粘土が荷重を受けることによって、粘土の間隙水が絞り出される現象を圧密現象として評価予測する方法が既存の方法として存在している（例えば、非特許文献２を参照）。
また、このような沈下挙動を抑制するための構造として、例えば下記の特許文献１に示される粘土系止水材による躯体支持構造が知られている。
「わが国における高レベル放射性廃棄物地層処分の技術的信頼性−地層処分研究開発第２次取りまとめ−総論レポート」、核燃料サイクル開発機構発行、１９９９年１１月 「土の力学」、河野伊一郎、八木則男、吉国洋著、技報堂出版発行、１９９０年５月
特許第４３２９０５４号公報
ところで、土の圧密現象として沈下速度と沈下量を予測する方法では、物体の下に存在していた粘土は密度増加して厚さを減じるものの、粘土粒子は物体の下から側方には移動しないので、沈下量を少なく見積もることが懸念される。
実際に乾燥密度０．２ｇ／ｃｍ３のベントナイトを水で飽和状態にした粘土材料は柔らかくて、容易に粘性変形する。図１２はそのような粘土に円柱状の物体（貫入用円柱５０ｍｍφ）を一定荷重（８００ｇ）で貫入させたときの挙動を写真で示したものであるが、明らかに粘性流体に近い挙動をしていることがわかる。すなわち、粘土を粘性流体とみなして、物体の沈下速度と沈下量をより保守的な値として推定する方法が求められていた。
また、図１１や図１２に示すような極端に形状変化を伴う挙動を正確に予測計算する簡便な計算方法が存在しなかった。
一方、粘土を純粘性流体として計算する際に、物体は剛体であるから、粘性流体と物体の境界は常に一定の速度で形状を変化させることなく、沈下変位のみを生じるはずであるが、そのような条件で純粘性流体の運動を計算するには、単純に粘性流体に一定の力を作用させる境界条件では計算できないという問題があった。
本発明は、上記に鑑みてなされたものであって、粘性材料で囲まれた物体の沈下挙動をより正確かつ簡便に予測計算することができる粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体を提供することを目的とする。
上記した課題を解決し、目的を達成するために、本発明に係る粘性材料で囲まれた物体の沈下予測計算方法は、粘性材料で囲まれた物体の沈下を予測計算する方法であって、粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定するステップ１と、物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算するステップ２と、計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ３と、計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ４と、前記ステップ３と前記ステップ４を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求めるステップ５とを有することを特徴とする。
また、本発明に係る他の粘性材料で囲まれた物体の沈下予測計算方法は、上述した発明において、あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、上述した粘性材料で囲まれた物体の沈下予測計算方法によって沈下速度値を求めるステップ６と、複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算するステップ７と、前記ステップ７で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を前記ステップ７と同様にして計算するステップ８と、前記ステップ８をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算するステップ９とをさらに有することを特徴とする。
また、本発明に係る粘性材料で囲まれた物体の沈下予測計算装置は、粘性材料で囲まれた物体の沈下を予測計算する装置であって、粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定する第１演算手段と、物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算する第２演算手段と、計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算する第３演算手段と、計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算する第４演算手段と、第３演算手段と第４演算手段による計算を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求める第５演算手段とを有することを特徴とする。
また、本発明に係る他の粘性材料で囲まれた物体の沈下予測計算装置は、上述した発明において、あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、上述した粘性材料で囲まれた物体の沈下予測計算装置によって沈下速度値を求める第６演算手段と、複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算する第７演算手段と、第７演算手段で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を第７演算手段と同様にして計算する第８演算手段と、第８演算手段をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算する第９演算手段とをさらに有することを特徴とする。
また、本発明に係る粘性材料で囲まれた物体の沈下予測計算プログラムは、粘性材料で囲まれた物体の沈下を予測計算する方法をコンピュータに実行させる計算プログラムであって、粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定するステップ１と、物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算するステップ２と、計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ３と、計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ４と、前記ステップ３と前記ステップ４を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求めるステップ５とをコンピュータに実行させることを特徴とする。
また、本発明に係る他の粘性材料で囲まれた物体の沈下予測計算プログラムは、上述した発明において、あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、上述した粘性材料で囲まれた物体の沈下予測計算方法によって沈下速度値を求めるステップ６と、複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算するステップ７と、前記ステップ７で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を前記ステップ７と同様にして計算するステップ８と、前記ステップ８をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算するステップ９とをコンピュータに実行させることを特徴とする。
また、本発明に係るコンピュータ読み取り可能な記録媒体は、上述した計算プログラムを記録したコンピュータ読み取り可能な記録媒体である。
本発明に係る粘性材料で囲まれた物体の沈下予測計算方法によれば、粘性材料で囲まれた物体の沈下を予測計算する方法であって、粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定するステップ１と、物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算するステップ２と、計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ３と、計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ４と、前記ステップ３と前記ステップ４を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求めるステップ５とを有するので、粘性材料で囲まれた物体の沈下挙動をより正確かつ簡便に予測計算することができるという効果を奏する。
また、本発明に係る他の粘性材料で囲まれた物体の沈下予測計算方法によれば、あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、上述した粘性材料で囲まれた物体の沈下予測計算方法によって沈下速度値を求めるステップ６と、複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算するステップ７と、前記ステップ７で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を前記ステップ７と同様にして計算するステップ８と、前記ステップ８をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算するステップ９とをさらに有するので、上記効果に加え、任意の時間経過後（例えば１００年、１０００年、１万年、１０万年後）における物体の沈下量の累計値をより正確かつ簡便に予測計算することができるという効果を奏する。
また、本発明に係る粘性材料で囲まれた物体の沈下予測計算装置によれば、粘性材料で囲まれた物体の沈下を予測計算する装置であって、粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定する第１演算手段と、物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算する第２演算手段と、計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算する第３演算手段と、計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算する第４演算手段と、第３演算手段と第４演算手段による計算を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求める第５演算手段とを有するので、粘性材料で囲まれた物体の沈下挙動をより正確かつ簡便に予測計算することができるという効果を奏する。
また、本発明に係る他の粘性材料で囲まれた物体の沈下予測計算装置によれば、あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、上述した粘性材料で囲まれた物体の沈下予測計算装置によって沈下速度値を求める第６演算手段と、複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算する第７演算手段と、第７演算手段で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を第７演算手段と同様にして計算する第８演算手段と、第８演算手段をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算する第９演算手段とをさらに有するので、上記効果に加え、任意の時間経過後（例えば１００年、１０００年、１万年、１０万年後）における物体の沈下量の累計値をより正確かつ簡便に予測計算することができるという効果を奏する。
また、本発明に係る粘性材料で囲まれた物体の沈下予測計算プログラムは、粘性材料で囲まれた物体の沈下を予測計算する方法をコンピュータに実行させる計算プログラムであって、粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定するステップ１と、物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算するステップ２と、計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ３と、計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ４と、前記ステップ３と前記ステップ４を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求めるステップ５とをコンピュータに実行させるので、粘性材料で囲まれた物体の沈下挙動をより正確かつ簡便に予測計算することができるという効果を奏する。
また、本発明に係る他の粘性材料で囲まれた物体の沈下予測計算プログラムは、上述した発明において、あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、上述した粘性材料で囲まれた物体の沈下予測計算方法によって沈下速度値を求めるステップ６と、複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算するステップ７と、前記ステップ７で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を前記ステップ７と同様にして計算するステップ８と、前記ステップ８をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算するステップ９とをコンピュータに実行させるので、上記効果に加え、任意の時間経過後（例えば１００年、１０００年、１万年、１０万年後）における物体の沈下量の累計値をより正確かつ簡便に予測計算することができるという効果を奏する。
図１は、粘性体内に埋設される物体に作用する荷重と抗力の説明図である。 図２は、本発明に係る粘性材料で囲まれた物体の沈下予測計算方法の実施例であり、沈下速度を評価するための反復計算のフローチャート図である。 図３は、ストークスの法則に基づく検証の解析モデルを示す図である。 図４は、本発明を適用した有限要素解析方法で求めた球体の沈下速度の一例を示す図である。 図５−１は、本発明による計算方法で沈下現象を予測計算するモデル（その１）を示す図であり、計算開始時における流動場の幾何学形状を示す図である。 図５−２は、本発明による計算方法で沈下現象を予測計算するモデル（その１）を示す図であり、計算途中段階における流動場の幾何学形状を示す図である。 図６は、図５−１のモデル（その１）に対する計算結果を時間の経過と沈下量の関係でグラフ表示した図である。 図７は、本発明による計算方法で沈下現象を予測計算するモデル（その２）を示す図である。 図８は、図７のモデル（その２）に対する計算結果を時間の経過と沈下量の関係でグラフ表示した図である。 図９は、沈下速度と粘性体の流速を鉛直断面上において図化した例を示す図であり、（１）は図５−１のモデル（その１）に関する図、（２）は図７のモデル（その２）に関する図である。 図１０は、従来の高レベル放射性廃棄物の地層処分におけるベントナイト系粘土材料の利用例を示す図であり、（１）は廃棄体を縦置き定置した場合の模式図、（２）は廃棄体を横置き定置した場合の模式図である。 図１１は、従来の高レベル放射性廃棄物を縦置き定置したときの長期沈下挙動のイメージ図である。 図１２は、水で飽和している粘土に円柱状の物体（貫入用円柱５０ｍｍφ）を一定荷重（８００ｇ）で貫入させたときの挙動写真図である。
以下に、本発明に係る粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体の実施の形態を図面に基づいて詳細に説明する。なお、この実施例によりこの発明が限定されるものではない。
本発明に係る粘性材料で囲まれた物体の沈下予測計算方法は、粘性材料で囲まれた物体の沈下を予測計算する方法であって、ステップ１〜５を有するものである。
ステップ１は、粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定するものである。
ステップ２は、物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算するものである。
ステップ３は、計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するものである。
ステップ４は、計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するものである。
ステップ５は、ステップ３とステップ４を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求めるものである。
上記構成によれば、粘性材料で囲まれた物体の沈下挙動をより正確かつ簡便に予測計算することができる。
また、本発明に係る粘性材料で囲まれた物体の沈下予測計算方法は、上記の構成に加え、次のステップ６〜９をさらに有してもよい。
ステップ６は、あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、上述した粘性材料で囲まれた物体の沈下予測計算方法によって沈下速度値を求めるものである。
ステップ７は、複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算するものである。
ステップ８は、ステップ７で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、次の時間区間の一定時間間隔における沈下量をステップ７と同様にして計算するものである。
ステップ９は、ステップ８をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算するものである。
上記構成によれば、任意の時間経過後（例えば１００年、１０００年、１万年、１０万年後）における粘性流体の流動場の違いを計算条件に反映することにより物体の沈下量の累計値をより正確かつ簡便に予測計算することができるという効果を奏する。
また、本発明に係る粘性材料で囲まれた物体の沈下予測計算装置としては、上記のステップ１〜５の処理内容をそれぞれコンピュータにより演算処理する第１〜第５演算手段を備えるように構成することができる。また、この構成に加えて、さらに上記のステップ６〜９の処理内容をそれぞれコンピュータにより演算処理する第６〜第９演算手段を備えるように構成してもよい。
また、本発明に係る粘性材料で囲まれた物体の沈下予測計算プログラムとしては、上記のステップ１〜５をそれぞれコンピュータに実行させるように構成することができる。また、この構成に加えて、さらに上記のステップ６〜９をそれぞれコンピュータに実行させるように構成してもよい。
また、本発明に係るコンピュータ読み取り可能な記録媒体は、上述した計算プログラムを記録したコンピュータ読み取り可能な記録媒体により構成することができる。
次に、上述した本発明の計算方法の基本原理等についてより具体的に説明する。
（１）粘土材料を純粘性流体としてモデル化
本発明では、以下のように、ベントナイト系粘土材料を純粘性流体とみなしてモデル化する。ここで、弾性、あるいは塑性的な性質を有していない粘性流体を純粘性流体という。純粘性流体の構成方程式は次式で表現される。
ここに、Ｓは応力テンソル、−ｐは応力の等方成分、Ｉは単位テンソルである。
また、Ｄはひずみ速度テンソル、Θは体積膨張速度であり、μを粘性率、λを第２粘性率と呼ぶ。
ここで、ｖは速度ベクトルである。
非圧縮流体では、Θ≡０であるため、上記の構成方程式は、次のようになる。
水、空気、アルコール、グリセリンなどの低分子量の流体は、ニュートン流体として分類される。ニュートン流体では、粘性率μはせん断ひずみ速度Ｄｉｊに依存しない。すなわち、
一方、溶融高分子やゴム、粒子懸濁液、石油などの流動は、ニュートン流体とは異なり、粘性率μはせん断ひずみ速度に依存するため、粘性率をせん断粘度、せん断依存粘度、あるいは、単に粘度、みかけの粘度とも呼ばれる。ここでは、以降、単に粘度ηと呼ぶ。せん断ひずみ速度Ｄｉｊの大きさをγドット（・）で表せば、
粘度ηのＳＩ単位は、［Ｐａ・ｓ］、ｃｇｓ単位は［Ｐｏｉｓｅ］である。
表１に代表的な流体の粘度のオーダーを示す。
（２）計算方法
（２−１）支配方程式
計算に用いる支配方程式としては、非圧縮性流体を仮定し、以下に示すナビエ・ストークスの運動方程式と連続式を扱う。
ここで、ρは流体の密度、ｆは荷重ベクトルである。
上記の純粘性流体の構成方程式を再記すると、
非圧縮性の仮定は、水、粘土ペーストや溶融樹脂などを取り扱う多くの流動問題では許容し得る近似とされている。また、一般に高粘度流動問題では、粘性力と圧力勾配が流動状況を支配的に決定するため、ナビエ・ストークス運動方程式の左辺第２項の非線形項を無視した近似（ストークス方程式）も多用されている。
（２−２）有限要素法
上述の問題を有限要素法により解析するため、粘性流体の流速と圧力を未知量とし、流速に対しては要素の頂点の節点による一次補間、圧力については要素内一定とした、いわゆる混合補間を採用するものとする。
（２−３）物体の沈下速度の計算方法
粘性流体（以下、粘性体ということがある。）内に埋設される物体の沈降速度は、図１に示すように物体に作用する荷重（自重−浮力）と流体抗力との釣り合い条件から、次式を用いて評価する。
ここで、左辺は、物体に作用する荷重であり、ｇは重力加速度、ρｗは物体の密度、ρｂは粘性体の湿潤密度、Ｖｗは物体の体積である。右辺はベントナイトの流動に伴って発生する抗力を表し、Ｓ・ｎは粘性体から物体の表面（表面積Ｓｗ）に作用する表面力である。
物体は変形しない剛体とみなし、解析対象は流動領域となる粘性体とする。粘性体から物体に作用する表面力分布は未知であり、その積分値が荷重と釣り合う条件を満足する。したがって、流動境界条件として、物体表面を一定沈降速度ｕｗの速度規定境界とみなし、図２のフローチャート図に示す反復計算を通じて釣り合い条件を満足するｕｗを求める。
すなわち、図２に示すように、（１）まず、物体の仮の沈下速度の初期値ｕ０を設定し（ステップＳ１１）、この初期値ｕ０を仮の沈下速度ｕｗとして設定するとともに、物体の沈下現象の原動力となる沈下力ｆｗ（物体に作用する重力から、粘性体から受ける浮力を差し引いた値）を設定する（ステップＳ１２）。（２）粘性体とその境界を設定し、粘性体を有限個の要素に分割して、要素を構成する節点の座標情報や物性情報等の計算情報を整理する。（３）設定した仮の沈下速度ｕｗで物体が沈下するという条件で粘性体の流動変形挙動を粘性流体として有限要素法により計算し（ステップＳ１３）、物体の表面が粘性体から受ける抗力ｆｒの積分値を求める（ステップＳ１４）。（４）計算された物体に作用する抗力ｆｒが仮設定した沈下力ｆｗに一致していない場合には（ステップＳ１５でＮｏ）、その差が縮まるような沈下速度ｕ＊ ｗに変更して（ステップＳ１６）、ステップＳ１３に戻り、（３）の計算を行う。（３）と（４）の計算を繰り返して、物体に作用する抗力ｆｒの計算値が、沈下力ｆｗに対して事前に設定してある所定の誤差範囲であった場合には計算を終了する（ステップＳ１５でＹｅｓ）。
なお、上記のステップＳ１６において、仮の沈下速度を変更する場合の補正方法としては、例えば、ｕ＊ ｗ＝ωｕ’ｗ＋（１−ω）ｕｗのような計算式を用いることができる。ここで、ｕ’ｗ＝ｕｗ（ｆｗ／ｆｒ）、０≦ω≦１である。
（２−４）沈下量の計算方法
前述の方法である時点における物体の沈下速度を精度良く推定計算することができるので、次は下記の方法で所定時間経過後の沈下量の累計値を計算する。
（１）沈下現象を予測したい期間（例えば１万年間）を複数の区間（例えば１００年×１００区間）に分割する。
（２）最初の時間区間における物体の粘性体の中における位置を境界条件に反映して、物体の沈下速度を前述（図２のフローチャート）の方法で計算する。
（３）最初の時間区間における物体の沈下速度は上記（２）で計算された一定の沈下速度で沈下するものとみなして、沈下量を時間積分して計算する。
（４）最初の時間区間における物体の沈下量を反映して、物体と粘性体の境界の座標位置をシフトした場合の各要素の幾何学形状に節点情報を修正する（有限要素解析では「リメッシュ」と称す。）。この修正によって、例えば図５−１と図５−２の違いのように、粘性体の流動する領域幅が狭い条件に変わる部位と広い条件に変わる部位が計算モデルに反映される。
（５）物体および粘性体の修正した位置情報に基づいて、次の第２の時間区間における上記（２）（３）の計算を実施し、再び上記（４）の計算条件の修正を実施して、さらにその次の第３の時間区間の沈下予測計算を実施する。
（６）上記（５）の計算を繰り返して、所定の期間までの沈下量計算を実施したら、計算を終了する。
（２−５）沈下速度の計算方法の検証
ここでは、前節（２−４）で示した本発明の抗力評価方法の妥当性をストークスの法則に基づき検証した結果を示す。
ストークスの法則によれば、粘度ηの流体内を一定速度ｕで移動する半径Ｒの球体に作用する抗力ｆｒは、次式で与えられる。
ストークスの法則は、運動量方程式において慣性項を無視し、無限遠で流体速度がゼロになる仮定を採用することによって、解析的に導かれるものであり、この式から理論的な計算をすることができる。例えば、沈下する物体は半径２ｍの球体とし、粘性体の粘度ηを１Ｅ＋１４Ｐａ・ｓとすると、沈下速度ｖ＝１Ｅ−３ｍ／ｓとなるときの球体が受ける抗力はストークスの法則により、１．８８５Ｅ＋１２Ｎと計算される。
本発明の計算方法では、計算する座標系として軸対称円柱座標系を採用し、図３に示すような解析モデルを採用した。左辺、上辺、右辺を流速ゼロの壁面境界、底辺を軸対称境界、球体表面を一定速の速度規定境界として設定した。
収束計算の目標となる抗力は１．８８５Ｅ＋１２Ｎとし、粘度ηを１Ｅ＋１４Ｐａ・ｓとした。仮の沈下速度を０．０２ｍ／ｓに設定して、図２に示す反復計算を実施した結果は図４に示すように、理論計算で設定した沈下速度ｖ＝１Ｅ−３ｍ／ｓに収束することを確認できた。
表２は、ストークスの法則より得られる抗力の厳密解と本発明による有限要素（ＦＥＭ）解析により計算された抗力の比較である。両者は良好な一致を示すことが確認できた。
次に、本発明による沈下量の予測計算例（その１、その２）について説明する。
（３）予測計算例その１
図５−１は、本発明による計算方法で沈下現象を予測計算する解析モデル（その１）の計算開始時における流動場の幾何学形状を示したものである。図５−１に示すように、解析対象は高レベル放射性廃棄物を粘性材料１内に縦置き定置した場合の廃棄体オーバーパック２の沈下である。計算条件を表３に示す。なお、図５−２は計算ステップの途中における流動場の違いの一例として、廃棄体オーバパックが０．３ｍ沈下したときの幾何学形状を示している。
計算結果を時間の経過と沈下量の関係でグラフ表示したものが図６である。図６に示すように、５０万年間で約０．６５ｍの沈下量と計算された。
（４）予測計算例その２
図７は、本発明による計算方法で沈下現象を予測計算する別の解析モデル（その２）を示したものである。図７に示すように、このモデルは上記の特許文献１に示された沈下抑制対策を、高レベル放射性廃棄物を縦置きした場合に適用した場合を模擬している。左右下部に配置してある沈下抑制材３は、縦置き処分孔の底部にリング状に配置した岩石である。
解析対象は図５−１、図６と同様に高レベル放射性廃棄物を粘性材料１内に縦置き定置した場合の廃棄体オーバーパック２の沈下であり、計算条件は上記のその１と同じ（表３を参照）にしてある。ただし、縦置き処分孔の底部にリング状に配置した沈下抑制材３と粘性体（緩衝材と称するベントナイト系粘土材料）の境界面を不動点になるように境界条件を設定して計算した。
計算結果を時間の経過と沈下量の関係でグラフ表示したものが図８である。図８に示すように、５０万年間で約０．２ｍの沈下量と計算された。上記の「予測計算例その１」の結果と比較して沈下量が低減していることから、沈下抑制材３による沈下抑制対策が効果を発揮していることを予測評価できた。
図９には、上記の２つの計算条件（その１、その２）で予測した計算結果である沈下速度と粘性体の流速を鉛直断面上において図化した例を示している。矢印は粘性体の流動方向を示している。沈下抑制材３による沈下抑制対策によって粘性体の流れが抑制されていることの違いが明確に示されている。
以上説明したように、本発明に係る粘性材料で囲まれた物体の沈下予測計算方法によれば、粘性材料で囲まれた物体の沈下を予測計算する方法であって、粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定するステップ１と、物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算するステップ２と、計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ３と、計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ４と、前記ステップ３と前記ステップ４を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求めるステップ５とを有するので、粘性材料で囲まれた物体の沈下挙動をより正確かつ簡便に予測計算することができる。
また、本発明に係る他の粘性材料で囲まれた物体の沈下予測計算方法によれば、あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、上述した粘性材料で囲まれた物体の沈下予測計算方法によって沈下速度値を求めるステップ６と、複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算するステップ７と、前記ステップ７で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を前記ステップ７と同様にして計算するステップ８と、前記ステップ８をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算するステップ９とをさらに有するので、上記効果に加え、任意の時間経過後（例えば１００年、１０００年、１万年、１０万年後）における物体の沈下量の累計値をより正確かつ簡便に予測計算することができる。
また、本発明に係る粘性材料で囲まれた物体の沈下予測計算装置によれば、粘性材料で囲まれた物体の沈下を予測計算する装置であって、粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定する第１演算手段と、物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算する第２演算手段と、計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算する第３演算手段と、計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算する第４演算手段と、第３演算手段と第４演算手段による計算を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求める第５演算手段とを有するので、粘性材料で囲まれた物体の沈下挙動をより正確かつ簡便に予測計算することができる。
また、本発明に係る他の粘性材料で囲まれた物体の沈下予測計算装置によれば、あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、上述した粘性材料で囲まれた物体の沈下予測計算装置によって沈下速度値を求める第６演算手段と、複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算する第７演算手段と、第７演算手段で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を第７演算手段と同様にして計算する第８演算手段と、第８演算手段をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算する第９演算手段とをさらに有するので、上記効果に加え、任意の時間経過後（例えば１００年、１０００年、１万年、１０万年後）における物体の沈下量の累計値をより正確かつ簡便に予測計算することができる。
また、本発明に係る粘性材料で囲まれた物体の沈下予測計算プログラムは、粘性材料で囲まれた物体の沈下を予測計算する方法をコンピュータに実行させる計算プログラムであって、粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定するステップ１と、物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算するステップ２と、計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ３と、計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ４と、前記ステップ３と前記ステップ４を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求めるステップ５とをコンピュータに実行させるので、粘性材料で囲まれた物体の沈下挙動をより正確かつ簡便に予測計算することができる。
また、本発明に係る他の粘性材料で囲まれた物体の沈下予測計算プログラムは、上述した発明において、あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、上述した粘性材料で囲まれた物体の沈下予測計算方法によって沈下速度値を求めるステップ６と、複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算するステップ７と、前記ステップ７で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を前記ステップ７と同様にして計算するステップ８と、前記ステップ８をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算するステップ９とをコンピュータに実行させるので、上記効果に加え、任意の時間経過後（例えば１００年、１０００年、１万年、１０万年後）における物体の沈下量の累計値をより正確かつ簡便に予測計算することができる。
以上のように、本発明に係る粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体は、ベントナイト系粘土材料などの湿潤粘土系の緩衝材に取り囲まれた放射性廃棄体のような物体の沈下予測計算に有用であり、特に、物体の沈下挙動をより正確かつ簡便に予測計算するのに適している。
１ ベントナイト系粘土材料（緩衝材）
２ 廃棄体
３ 沈下抑制材

Claims (7)

  1. 粘性材料で囲まれた物体の沈下を予測計算する方法であって、
    粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定するステップ１と、
    物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算するステップ２と、
    計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ３と、
    計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ４と、
    前記ステップ３と前記ステップ４を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求めるステップ５とを有することを特徴とする粘性材料で囲まれた物体の沈下予測計算方法。
  2. あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、請求項１に記載の粘性材料で囲まれた物体の沈下予測計算方法によって沈下速度値を求めるステップ６と、
    複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算するステップ７と、
    前記ステップ７で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を前記ステップ７と同様にして計算するステップ８と、
    前記ステップ８をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算するステップ９とをさらに有することを特徴とする請求項１に記載の粘性材料で囲まれた物体の沈下予測計算方法。
  3. 粘性材料で囲まれた物体の沈下を予測計算する装置であって、
    粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定する第１演算手段と、
    物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算する第２演算手段と、
    計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算する第３演算手段と、
    計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算する第４演算手段と、
    第３演算手段と第４演算手段による計算を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求める第５演算手段とを有することを特徴とする粘性材料で囲まれた物体の沈下予測計算装置。
  4. あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、請求項３に記載の粘性材料で囲まれた物体の沈下予測計算装置によって沈下速度値を求める第６演算手段と、
    複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算する第７演算手段と、
    第７演算手段で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を第７演算手段と同様にして計算する第８演算手段と、
    第８演算手段をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算する第９演算手段とをさらに有することを特徴とする請求項３に記載の粘性材料で囲まれた物体の沈下予測計算装置。
  5. 粘性材料で囲まれた物体の沈下を予測計算する方法をコンピュータに実行させる計算プログラムであって、
    粘性材料を粘性流体とみなして、この粘性流体で囲まれている物体の沈下力を、物体に作用する重力から、粘性流体から受ける浮力を差し引いた値として設定するステップ１と、
    物体の沈下速度を任意の速度である仮の沈下速度に設定して、設定した仮の沈下速度で物体が沈下するときの物体周囲の粘性流体の粘性流れ速度を、粘性流体を単純な粘性流体とみなして計算するとともに、粘性流体から物体の表面に作用する抗力を計算するステップ２と、
    計算された抗力の積分値が設定した沈下力よりも大きい場合には、仮の沈下速度を若干小さく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ３と、
    計算された抗力の積分値が設定した沈下力よりも小さい場合には、仮の沈下速度を若干大きく設定して、改めて粘性流体から物体の表面に作用する抗力を計算するステップ４と、
    前記ステップ３と前記ステップ４を繰り返して、計算された抗力の積分値が事前に設定した沈下力の誤差範囲に収束した時点で計算を停止し、そのときの仮の沈下速度を沈下力に相当する沈下速度値として求めるステップ５とをコンピュータに実行させることを特徴とする粘性材料で囲まれた物体の沈下予測計算プログラム。
  6. あらかじめ沈下量を予測したい時間を複数の一定時間間隔に分割して、請求項５に記載の粘性材料で囲まれた物体の沈下予測計算方法によって沈下速度値を求めるステップ６と、
    複数の一定時間間隔に分割したときの最初の時間区間の一定時間間隔において、求めた沈下速度値に等しい沈下速度で物体が移動するものとみなして、沈下変位量を時間積分することにより、この一定時間間隔における沈下量を計算するステップ７と、
    前記ステップ７で計算した沈下量に応じて、粘性流体の移動に見合った位置に当初設定した有限要素解析における計算要素の節点が移動したとみなして、節点の座標を補正し、粘性流体の流動場の違いを計算条件に反映して次の時間区間の一定時間間隔における沈下量を前記ステップ７と同様にして計算するステップ８と、
    前記ステップ８をあらかじめ沈下量を予測したい時間として設定した時間まで繰り返して最終的な沈下量を計算するステップ９とをコンピュータに実行させることを特徴とする請求項５に記載の粘性材料で囲まれた物体の沈下予測計算プログラム。
  7. 請求項５または６に記載の計算プログラムを記録したコンピュータ読み取り可能な記録媒体。
JP2013124669A 2013-06-13 2013-06-13 粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体 Active JP6041151B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013124669A JP6041151B2 (ja) 2013-06-13 2013-06-13 粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013124669A JP6041151B2 (ja) 2013-06-13 2013-06-13 粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体

Publications (2)

Publication Number Publication Date
JP2015001399A JP2015001399A (ja) 2015-01-05
JP6041151B2 true JP6041151B2 (ja) 2016-12-07

Family

ID=52296017

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013124669A Active JP6041151B2 (ja) 2013-06-13 2013-06-13 粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体

Country Status (1)

Country Link
JP (1) JP6041151B2 (ja)

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3998051B2 (ja) * 2001-12-19 2007-10-24 清水建設株式会社 埋設廃棄物の止水充填法と充填材
JP4640671B2 (ja) * 2005-04-11 2011-03-02 清水建設株式会社 地盤・構造物変形量予測方法およびプログラム
JP4822098B2 (ja) * 2005-06-30 2011-11-24 清水建設株式会社 地盤・構造物沈下量予測方法およびプログラム

Also Published As

Publication number Publication date
JP2015001399A (ja) 2015-01-05

Similar Documents

Publication Publication Date Title
Soga et al. Trends in large-deformation analysis of landslide mass movements with particular emphasis on the material point method
Larese et al. Numerical and experimental study of overtopping and failure of rockfill dams
Cremonesi et al. A basal slip model for Lagrangian finite element simulations of 3D landslides
Egholm et al. The mechanics of clay smearing along faults
Hariri-Ardebili et al. Integrative seismic safety evaluation of a high concrete arch dam
Mehrzad et al. Centrifuge and numerical models to investigate liquefaction-induced response of shallow foundations with different contact pressures
Adamidis et al. Post-liquefaction reconsolidation of sand
Hamad et al. Modelling of membranes in the material point method with applications
Zhao et al. Simplified effective stress modeling of shear wave propagation in saturated granular soils
Hu et al. A coupled two‐phase fluid flow and elastoplastic deformation model for unsaturated soils: theory, implementation, and application
Huang et al. Computational fluid dynamics modeling of post-liquefaction soil flow using the volume of fluid method
Wang et al. 3D stability of partially saturated soil slopes after rapid drawdown by a new layer-wise summation method
Zhao et al. A finite element approach for predicting the full resistance profile of a spudcan deeply penetrating in dense sand overlying clay
Di et al. An operator‐split ALE model for large deformation analysis of geomaterials
Shahbodagh Khan et al. Dynamic analysis of strain localization in water‐saturated clay using a cyclic elasto‐viscoplastic model
Kim et al. Modeling of soil deformation and water flow in a swelling soil
Tamayo et al. On the validation of a numerical model for the analysis of soil-structure interaction problems
Becker et al. A multiresolution strategy for solving landslides using the Particle Finite Element Method
Wollny et al. A hierarchical sequential ALE poromechanics model for tire‐soil‐water interaction on fluid‐infiltrated roads
Salgado et al. Advances in the solution of geotechnical boundary-value problems
Peng et al. Application of a coupled hydro‐mechanical interface model in simulating uplifting problems
Nguyen et al. Fluidization of soil under increasing seepage flow: an energy perspective through CFD-DEM coupling
JP6238043B2 (ja) 地盤−構造物系の液状化時の残留変形簡易算定方法
KR101813903B1 (ko) 대변형 해석기법을 이용한 사방댐 설계 시스템 및 방법
JP6041151B2 (ja) 粘性材料で囲まれた物体の沈下予測計算方法、計算装置、計算プログラムおよびコンピュータ読み取り可能な記録媒体

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20151118

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160923

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20161004

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20161026

R150 Certificate of patent or registration of utility model

Ref document number: 6041151

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150