フジタ技術研究報告第 47 号 2011 年

地盤の非線形履歴特性のモデル化

概中 川 太 郎

小 林 勝 已佐 々 木 聡佐 々 木 仁要

広範囲のひずみレベルで土質試験結果を精度良く近似し、かつ容易に逐次非線形地震応答解析に取り込むことを目的として、骨格曲線には双曲線モデルを修正したモデル、履歴曲線には Masing 則を満足しつつ、パラメータ α にひずみ依存性をもたせた R-O モデルを利用した非線形モデル（以降 H-R モデルと称す：Hyperbolic & R-O）を定式化した。
このことにより、提案する非線形モデルが、双曲線モデルと R-O モデルの簡便性を保ちながらも広いひずみ領域において土質試験結果の G/G0-γ 関係と h-γ 関係を共に精度良く近似することができることを示した。
また、地震応答解析に適用した例を示した。

Modeling of non-linear hysteretic soilAbstract
To represent the shear modulus and damping ratio of soils under dynamic loading conditions
covering a wide range of shear strains, we proposed a hybrid method of modified
Ramberg-Osgood model and modified hyperbolic model.
A mathematical theory of skeleton curve was formulated by adding two parameters to thehyperbolic model, and the mathematical theory of hysteresis curve was formulated based on theparameter of Ramberg-Osgood model being dependent on strain.
By doing so, this nonlinear model combines accuracy and convenience. Examples of seismic
response analysis employing these models are also preseneted.

キーワード： 地盤 非線形 履歴特性 地震応答解析
双曲線モデル R-O モデル－59－

フジタ技術研究報告 第 47 号

§１．はじめに

3) 履歴曲線のパラメータ α を過去の経験最大ひずみの関数とすることにより、上記 1)、2)の条件を満足した上でMasing 則を適用させる。

1.1 背景土質試験結果のせん断剛性比 G/G0 と減衰定数 h のひずみ依存性を近似する非線形モデルとして、

2.2 H-R モデルの骨格曲線

R-O(Ramberg-Osgood)モデル 1)および双曲線モデル 2)が双曲線モデルは(1)式で示される 4)。

よく使われている。しかし R-O モデルでは、せん断ひずみ

τ=の増大にともない応力が上昇するため、大ひずみ領域での適用に注意が必要であり 3)、双曲線モデルは、履歴曲線がγ ⋅ G01+ γ / γ f(1)

ここに、τ：せん断応力、 γ：せん断ひずみ、G0：初期せ定式化されていないため理論的な減衰特性が定められな
い 4)など、小ひずみ～大ひずみの広い領域での剛性特性（G/G0-γ 関係）と減衰特性（h-γ 関係）の二つを共に精度良く近似することが難しい。

ん断剛性、γf：基準ひずみである。
なお後述の G は G=τ/γ で定義される割線剛性である。
本論で提案する修正HP モデルは、G/G0-γ 関係の近似精度を上げるために(1)式に 2 つのパラメータを追加したものであり、(2)式で表される。

1.2 目的本論では、広範囲のひずみレベルで土質試験結果を精度良く近似し、かつ容易に逐次非線形地震応答解析に取りτ=γ ⋅ G0(2)

1 + {(γ − γ m ) /(γ f − γ m )}λ込むことを目的として、骨格曲線には双曲線モデルを修正

ただし、γ<γm の場合は τ =γ･G0 とする。

したモデル、履歴曲線には Masing 則を満足しつつ、パラ

ここで、λ と γm は今回導入したパラメータであり、骨格曲線の形状の自由度を高めるものである。

メータ α にひずみ依存性をもたせた R-O モデルを利用しλ<1.0 としたときは、基準ひずみより小ひずみ領域のた非線形モデル（以降 H-R モデルと称す：Hyperbolic &G/G0 を双曲線モデルより小さくするとともに基準ひずみよりR-O）を定式化した。

大ひずみ領域での G/G0 を大きくする効果がある。また
λ>1.0 の場合は逆の効果を有する。

§２．定式化γm は γ＝γm で G/G0=1.0 とする効果を有する。すなわ
2.1 H-R モデルの特徴

ち剛性低下が始まるひずみを γm により規定でき、またそのH-R モデルの概念を図 1 に示し、H-R モデルの特徴を影響は剛性低下が始まる微小ひずみ近傍のみである。このことは、λ を導入することにより生じる小ひずみ領域の誤差以下に示す。

を補正する効果をもたらすことを意味する。
以上の λ および γm の導入により、小ひずみから非線形骨格曲線は双曲線モデル
を修正した修正 HP モデル
τ

解析で対象となると考えられる大きさのひずみ領域まで、土質試験結果の G/G0-γ、h-γ 関係を精度良く近似できることを後述の例で示す。

τa2.3 H-R モデルの履歴曲線-γaγa-τa

本論では G/G0-γ 関係の近似精度を確保するために双γ曲線モデルを修正して用いているが、これに直接 Masing則を適用すると、修正 H-D モデルと同様に減衰が過大に履歴曲線はひずみ依存性の α を用いた
R-O モデルなる。そこで、これを防ぐために R-O モデルを利用する。
R-O モデルの骨格曲線は(3)式で表される。

図１ H-R モデルの概念

γ=1) 骨格曲線に双曲線モデルを修正したモデル（以降、修正 HP(双曲線：Hyperbolic)モデルと称す）を用いる。

τ ⎧⎛⎞⎨1 + α ⎜ τ τ ⎟⎝ 0⎠
G0 ⎩r −1⎫⎬⎭(3)

ここに、τ0= G0･γf、α、r：パラメータである。
R-O モデルに Masing 則を適用すると、履歴曲線は骨格

2) 履歴曲線に R-O モデルを用いる。ただし、R-O モデル

曲線を τ 軸と γ 軸の両方に対して 2 倍したものとすることのパラメータ α にはひずみ依存性を考慮する。

から、(4)式で表せる 5) （τa、γa は図 1 を参照）。

－60－

地盤の非線形履歴特性のモデル化

τ −τ a2⎛γ −γa ⎞= f⎜⎟⎝ 2 ⎠

すなわち、パラメータ α は骨格曲線上ではひずみに依(4)存して刻々変化し、履歴曲線上では経験最大ひずみ γa で規定される定数となる。

パラメータ r は(3)式に Masing 則を適用した場合の、(5)式で示される減衰定数から求めることができる。
2(r − 1) ⎛ G ⎞⋅ ⎜1 − ⎟h=π (r + 1) ⎜⎝ G0 ⎟⎠2.4 パラメータαのひずみ依存性と近似精度(5)

ここで、後述する地震応答解析に用いた地盤モデルにおいて土質試験結果によるパラメータ α がひずみ依存性を次に、(3)式を α について解くと(6)式が得られる。
⎛⎝

α = ⎜ G0γ ⎞− 1⎟τ ⎠⎛τ ⎞⎜⎜ ⎟⎟⎝τ0 ⎠r −1

有することを示すとともに、H-R モデルで土質試験結果を良好に近似できることを示す。

(6)

図 2 に、(6)式で求めた土質試験結果(Test)の α と、(7)(6)式より、α は定数ではなく、せん断ひずみ γ が大きく式を用いて計算した H-R モデルの α の比較を示す。なお、なるに従い大きくなるというひずみ依存性があることがわか

基準ひずみ γf は粘性土で 0.13%、砂質土 1 で 0.09%、砂る。しかしながら、従来α は定数として扱われており、このこ質土 2 で 0.04%である。通常、R-O モデルでは基準ひずみとが R-O モデルの骨格曲線が土質試験結果を必ずしも精での値を定数として用いることが多いため、その値（R-O と度良く模擬できないことの一因と考えられる。

表記、α=2r-1）も併記する。
図の土質試験結果の α は基準ひずみより小さいひずみ

ここで、α をひずみ依存性として扱う際の α の求め方について、以下に示す 2 つの方法が考えられる。

領域では基準ひずみ時の α より小さく、基準ひずみを超え

1) (6)式の右辺の G0、γ、τ、τ0、r に土質試験から得らると急激に大きくなっており、H-R モデルはこれを精度良く近似できている。

れた値を代入して求める方法
2) (2)式の修正 HP モデルの骨格曲線を利用する方法
第一の方法は、地震応答のように任意のひずみに対応ひずみが土質試験結果の範囲を超える場合は外挿する必αするためには、内挿もしくは外挿する必要がある。特に応答要があるが、適切な外挿方法を見つけることは困難である。
本論で提案するのは第二の方法であり、α をひずみの関数として扱うものである。この場合にも第一の方法と同様、1412粘性土
10864200.0001 0.001

外挿方法が問題となるが、本論では骨格曲線として修正HP モデルを採用することにより対処する。
修正 HP モデルは非常に精度良く G/G0-γ 関係を近似できるため、α のひずみ依存性も精度良く近似できる。また、
α

骨格曲線が定義されているため、ひずみ依存性を有するα を定式化できる。すなわち、(2)式の τ を(6)式に代入すれば(7)式となる。
λλ⎛ γ − γ m ⎞ ⎡ γ f ⎧⎪ ⎛ γ − γ m ⎞ ⎫⎪⎤⎟ ⋅ ⎢ ⎨1 + ⎜⎟ ⎬⎥α = ⎜⎜⎟⎜⎟⎝ γ f − γ m ⎠ ⎢⎣ γ ⎪⎩ ⎝ γ f − γ m ⎠ ⎪⎭⎥⎦r −1(7)ここで地震応答解析においては、(7)式の α を刻々変化するひずみに依存する関数とすると、載荷から除荷あるいαは逆に除荷から載荷に変化して骨格曲線から外れる場合にも α は応答ひずみに応じて常に刻々と変化するため、履歴曲線を Masing 則で表すことができない。そこで、過去の正側あるいは負側の最大ひずみに達した時点の α の値を定数として、最大ひずみが更新されるまでその値を用いるとする。このことより、Masing 則の適用が可能となる。

14砂質土11210864200.0001 0.001

Test
H-R (7)式
R-O0.010.1γ(%)

110110110Test
H-R (7)式
R-O0.010.1γ(%)

14
砂質土21210Test
8H-R (7)式
R-O64200.0001 0.001 0.010.1γ(%)

図 2 パラメータαのひずみ依存性－61－フジタ技術研究報告 第 47 号

2.5 H-R モデルのパラメータの計算手順

衰定数を満足する Masing 則に基づくものとしているが、骨H-R モデルのパラメータの計算手順を以下に示す。

格曲線 G/G0 がそれぞれ異なる。パラメータ α に関しては、

1) (2)式において、基準ひずみ γf は G/G0=0.5 となるひず

H-R、HP モデルでは(7)式によるひずみ依存の値としておみ、最小ひずみ γm は G/G0=1.0 から剛性低下が始まり、R-O モデルは既往のモデルと同じく定数としている。

るひずみとすると、実質的に求めるパラメータは λ のみである。そこで、骨格曲線の実験値と(2)式による計算値の残差平方和が最小となるように最小二乗法で λ を求める。

表 2 非線形モデルの違い
モデル
曲線
骨格曲線

2) 減衰定数の実験値と(5)式の計算値の残差平方和が最小となるように、最小二乗法でパラメータ r を求める。

履歴曲線3) パラメータ α は、上記パラメータ値、および載荷履歴あるいは地震応答解析中の過去の最大ひずみ γa を(7)備 考式に代入することにより求められる。

H-R

R-OHP(2)式

(3)式
(1)式
Masing 則に基づく(4)式
パラメータ α はパラメータ α はパラメータ α はひずみ依存性を
ひずみ依存性をr-1定数：α=2示す(7)式
示す(7)式
骨格曲線を双曲
本論の提案モデ 既往の R-O モデ 線モデルとし、履ル
ルと同じ
歴曲線は H-R と同じモデル§３．G/G0-γ、h-γ 関係の近似精度と地震応答解析
3.3 G/G0-γ、h-γ 関係の近似精度
本章では、H-R モデルの土質試験結果の良好な近似精

H-R、R-O、HP 各モデルによる G/G0-γ、h-γ 関係の近度を R-O、HP(双曲線)の両モデルとの比較により示すとと似精度を比較検討する。図 3 に土質試験結果の G/G0-γ、もに、地震応答解析結果のτ-γ関係に及ぼす影響についh-γ 関係とそれぞれのモデルでの近似結果を示す。なお、て考察する。

H-R モデルのパラメータの計算手順は 2.5 節に示したが、他のモデルのパラメータは次のように扱った。
表 2 に示すように、HP モデルと H-R モデルの相違は骨3.1 地盤モデル
表 1 に地盤モデルを示す。地盤モデルは文献 3)を参考格曲線のモデル化にあり、HP モデルでは(1)式においてに表層地盤の固有周期 T=0.31(s)の Case1(地盤-1)と、固基準ひずみ時の G を用いた。また R-O モデルのパラメータ有周期 T=0.80(s)の Case2 の 2 種類設定した。

r は H-R モデルと同じであり、パラメータ α は定数として基

ただし、非線形特性を示す G/G0-γ、h-γ 関係は土質試準ひずみ時の値を用いた。

験結果を元に設定し、Case1 と Case2 のそれぞれの層に粘性土についてみると、H-R モデルでは全ひずみ領域対して、土質試験を行った地盤に似た非線形タイプをあてにおいて G/G0-γ 関係を良く近似できている。これは(2)はめた。非線形タイプは Case1 と Case2 共通で、粘性土 1式に示す 2 つのパラメータ λ と γm の効果と考えられ種類と砂質土 2 種類とした。

る。h-γ 関係は 0.01%より大ひずみ領域では試験結果を良く近似できており、パラメータ α のひずみ依存性の効果と考えられるが、0.01%より小ひずみ域では試験結

表１ 地盤モデル
層番号
12345
層番号
1234

深度 D(m)
3.2
5.7
10.0
17.6
基盤
深度 D(m)17.7
28.536.1
基盤

Case1(地盤-1)
ρ(t/m3)Vs(m/s)1.71301.8340
1.72801.93802.1510
Case2(地盤-3)
ρ(t/m3)Vs(m/s)1.8
1301.92501.93601.9430

果より小さい値である。これは Masing 則を用いている

非線形 Type
砂質土 1
砂質土 2
粘性土
砂質土 2―

ことによるもので、(5)式より G/G0 が 1 に漸近するに従い h が 0 に近づくためである。
R-O モデルの G/G0 は基準ひずみより小ひずみ領域では試験結果より小さく、基準ひずみより大きいひずみ領域では試験結果より大きい。この傾向は後述する砂質土 1、2 に非線形 Type
粘性土
砂質土 1
砂質土 2―比較して顕著である。
一方 h は基準ひずみより小さいひずみ領域では試験結果より大きく、基準ひずみより大きい領域では試験結果より小さい結果である。
hmax を変えることで R-O モデルの近似度が変わってくる

3.2 各モデルの骨格曲線と履歴曲線が、対象とした粘性土の R-O モデルの近似精度が良くない各モデルの骨格曲線と履歴曲線の概要を表 2 にまとめて理由の一つとして、α を基準ひずみ時の定数としたことが示す。
履歴曲線は H-R、R-O、HP 各モデルとも(5)式に示す減－62－挙げられる。

地盤の非線形履歴特性のモデル化
1.01.0粘性土

0.00.0001 0.001

0.6Test
H-R
R-O
HP

0.4
0.2

0.01
0.1
γ(%)

0.140.12
Test
0.10
H-R
R-O
0.08
HP
0.06
0.04
0.02
0.00
0.0001 0.001 0.01
0.1
γ(%)

10.00.0001 0.001

10G/G0Test
H-R
R-O
HP1Test
H-R
R-O
HP

0.10

0.150.050.010.1
γ(%)

1100.010.1
γ(%)

10.10

Test
H-RR-O
HP0.05

砂質土10.00
0.0001 0.001

10Test
H-R
R-O
HP

0.20

0.15

10.4

0.00.0001 0.001

100.20

粘性土

0.6

0.2

0.01
0.1
γ(%)

砂質土2

0.8

h0.4

G/G00.6

0.2

h0.8

h

G/G00.81.0
砂質土1

砂質土2

0.00
0.0001 0.001

100.010.1
γ(%)

110110

図３ G/G0-γ 関係，h-γ 関係

粘性土0.30.2

砂質土10.3

0.010.1

1100.2

Test
H-R
R-OHP

砂質土2

0.10.10.10.00.0001 0.001

0.2

Test
H-R
R-O
HP

τ/G00.3

Test
H-R
R-O
HP

τ/G0

τ/G00.40.40.40.50.00.0001 0.001

0.010.11100.00.0001 0.001

0.010.1

γ(%)

γ(%)

γ(%)

図４ 基準化骨格曲線 τ/G0HP モデルの G/G0-γ 関係は試験結果を概ね良好に近の近似結果の比較により、 H-R モデルは粘性土、砂質土似しているが、H-R モデルと比べると基準ひずみを境に大に関わらず、他のモデルに比べ全ひずみ領域で良い近似小関係が逆転しており、H-R モデルの方が近似精度は良く、精度であることを示した。

パラメータ λ の効果が見られる。

ここでは、土質試験の近似誤差が地震応答解析結果に砂質土 1 に関しても同様の傾向であるが、各モデル間の及ぼす影響について検討する。
各ケースとも用いた地震波は告示波レベル 2 のランダム差は他の土質に比べそれほど大きくはない。
砂質土 2 に関しては、H-R モデルと HP モデルの差は他位相模擬地震波である。の土質に比べ大きい。G/G0 は基準ひずみを境に両者の大図 5 にそれぞれ Case1、Case2 の場合の H-R、R-O、小関係が逆転しており、H-R モデルの方が試験データの近HP 各モデルの解析結果の最大ひずみ発生層での τ-γ似精度が良く、λ の効果が見られる。またこの差は h の近関係を示す。図中の太線は骨格曲線、細線は履歴曲線を似精度にも影響を及ぼしている。

示す。Case1、Case2 ともそれぞれ非線形モデルの違いに

ここで、改めて各モデルの違いを、τ を G0 で除した基準関わらず最大ひずみは同一層で発生しており、Case1 は化骨格曲線で比較したものが図 4 である。γ が 0.1%より大GL-2～3.2m の砂質土 1、Case2 は GL-16.7～17.7m のきいひずみ領域で各モデルの差が先の G/G0 に比べ、より粘性土である。
Case1 では、各モデル間の相違はそれほど大きくないが、明瞭に出ている。

R-O モデルは若干他より等価剛性が大きい。これは基準ひこのひずみ領域での差が地震応答結果に大きく影響すずみである γ=0.09%より大きいひずみ領域で R-O モデルることについては後述する。

が剛性を過大評価していることによる。
Case2 では各モデル間の相違が大きく、R-O は等価剛

3.4 地震応答解析結果
前節で G/G0-γ、h-γ 関係の土質試験結果と各モデル性が大きく急峻な形状、HP はひずみが 0.5％を超えると急－63－

フジタ技術研究報告 第 47 号

激に進み横長の形状、H-R は両者の中間的な形状であG/G0-γ 関係については全ひずみ領域において良く近る。

似でき、h-γ 関係については 0.01%より大きなひずみ最大ひずみは H-R は 1%強のひずみ、HP は 2%強のひ領域で良く近似できることを示した。
4) 逐次非線形地震応答解析で H-R、R-O、HP 各モデルずみとなっているが、この差異はこの層の大ひずみ領域での剛性評価の違いに因るものと考えられる。の比較を行った。土質試験結果に対する大ひずみ領域

すなわち、図 4 の粘性土の基準化骨格曲線をみると、での G/G0-γ、h-γ の近似精度は H-R が良く、R-O は0.3%以上のひずみでは H-R と HP の基準化骨格曲線の剛性を硬めに、HP は剛性を柔らかめ評価する傾向で差が徐々に大きくなり、さらに 1%を超えるとその差はかなりあったが、地盤状況によってはその結果が応答結果に大きくなり、割線剛性は H-R に比べ HP はかなり小さくなる。

顕著に表れることを示した。
5) H-R モデルの実用的な利点は次の点である。

したがってその差が応答結果に大きく影響したものと考えら

・ H-R モデルのパラメータは、R-O モデルのパラメータれる。と同程度の簡便さで求めることができる。
Case220-0.3 -0.2 -0.1
-10 0 0.1 0.2 0.3-20-1

H-RH-R

析センター平澤光春氏に深く感謝いたします。
参考文献

2τ (kN/m )0-20-10121) Jennings, P. C. : Periodic Response of General-40

R-O

R-O-30γ(%)-80γ(%)

3080

Yielding Structure , ASCE, Vol. 90, No. EM2,
pp.131-167, 19642) Hardin,B.O. and Drnevich,V.P. : Shear Modulus andDamping in Soils, Design Equations and Curves,

Case22040

Proc. ASCE, SM7, pp.667-692, 1972210τ (kN/m )2τ (kN/m )2

本報をまとめるにあたり、終始指導いただいた時系列解

40-0.3 -0.2 -0.1
-10 0 0.1 0.2 0.3

τ (kN/m )謝辞

80

200-0.3 -0.2 -0.1
-10 0 0.1 0.2 0.3-20だけで、H-R モデルを組み込むことができる。

2Case210Case11-80γ(%)

30-200-40-30γ(%)

Case1

直すように既存の R-O モデルのプログラムを改良する

22

τ (kN/m )0-20

テップの中で、最大ひずみが更新されるごとに規定し4010

τ (kN/m )Case1

・ 地震応答解析に関しては、パラメータ α を各時間ス

8030HP-23) 日本建築学会：建物と地盤の動的相互作用を考慮した応0-1012答解析と耐震設計，2006-404) 国生剛治，桜井彰雄：Modified Hardin-Drnevich モデル

HP-30
γ(%)

-80
γ(%)

(a) Case1(b) Case2について，土木学会第 33 回年次学術講演会梗概集 第Ⅲ部 pp.116-117，1978
5) 地盤工学会：地盤・基礎構造物の耐震設計，2001

図５ 最大ひずみ発生層でのτ-γ関係§４．おわりに

ひ と こ と以下にまとめを示す。

東北地方太平洋沖地震による建物

1) G/G0-γ 関係を精度良く近似するために、双曲線モデ

被害は記憶に新しいところですが、今ルの骨格曲線に 2 つのパラメータを追加した修正 HP

後発生する地震に対する被害低減にモデルを提案した。

関して、地盤解析技術という側面から

2) h-γ 関係を精度良く模擬するために、履歴曲線の R-O寄与できればと思います。

モデルのパラメータ α にひずみ依存性を考慮し，かつR-O モデルと修正 HP モデルを結合した H-R モデルを中川 太郎

提案した。
3) H-R モ デ ル は 土質試験結果と の 比較に お い て 、－64－