杭応答変位法オプション 計算編のサブセクション

概要

機能一覧

機能一覧は以下の通りです。

項目仕様備考
解析モデル単杭 ・ 群杭フレーム(基礎梁剛または基礎梁剛性考慮、平面的なねじれ考慮)単杭の場合は弾性解析により各杭の慣性力分担を決定します。
水平地盤ばね基礎指針 ・ Francis基礎指針:塑性水平地盤反力度Pyを応力度の上限とする指定も可能です。
Francis:Bromsの耐力式により、双曲線モデルとします。
場所打ち鉄筋コンクリート杭,場所打ち鋼管コンクリート杭
断面切り替え杭頭,杭中間,杭先端など自由に切り替え可能
線形・非線形(Tri-Linear)
杭頭回転拘束条件:回転自由,回転固定,つなぎ梁の剛性を考慮した回転拘束
軸力短期・終局各ケースごとに設定可能です。
付加モーメント短期・終局各ケースごとに設定可能です。水平荷重の前に杭頭に作用させます。
地盤変位地盤No.ごとに地盤変位設定地層分割と異なる分割の場合は自動的に線形補間します。
応力解析結果の重ね合わせ単純累加/SRSSモーメント,せん断力で異なる重ね合わせ方法を設定することも可能です。各解析結果の重ね合わせではなく、同時載荷も可能です。同時載荷において、正位相・逆位相を考慮した検討も可能です。
M-N相関短期・終局終局耐力についてはコンクリート圧縮特性をe関数法により算出します。
図化機能モーメント、せん断力、変位、回転角分布図
M-N相関図
CSVファイル出力杭情報、地盤情報、解析結果
Wordファイル出力杭情報、地盤情報、解析・設計条件、応力・変位分布図、MN相関図、検定表

計算フロー

解析条件 - 応力解析モデル

解析モデルは単杭モデルと群杭フレームモデル(基礎梁剛)から選択できます。

以下は単杭モデルのイメージ図です。
後述の群杭フレームモデルを弾性解析することにより慣性力の分配率を計算し、分配された慣性力を各単杭モデルに作用させるモデルです。
慣性力の分配を決める計算の際には、杭体を弾塑性として指定した場合でも弾性として解析した結果を用います。

以下は群杭フレームモデルのイメージ図です。
全慣性力を剛床位置に作用させ、応力解析により慣性力が自動的に分配されます。

単杭モデル

単杭モデルの計算フローは以下となります。

    1. 共通
    • 1.1. 水平地盤特性の入力および深度の関係から、水平地盤ばね値およびばね耐力を算出します。
    • 1.2. 杭断面の入力から、杭の剛性、耐力および骨格曲線を算出します。
    • 1.3. 基礎梁の入力から、つなぎ梁の剛性を算出します。
    1. 単杭モデル
    • 2.1. 群杭フレームモデル(基礎梁剛)の弾性解析により、各杭のせん断力分担比を計算します。
    • 2.2. 計算された分担比に応じた慣性力により単杭モデルの慣性力解析を行います。
    • 2.3. 単杭モデルの地盤変位解析を行います。
    • 2.4. 慣性力解析の場合は指定により応力を補正します(非公開機能)。
    • 2.5. 慣性力、地盤変位結果を重ね合わせ検定応力を算出します。
    • 2.6. 検定応力を用いて断面検定を行います。

群杭フレームモデル

群杭フレームモデルの計算フローは以下となります。

応力重ね合わせの場合

    1. 共通
    • 1.1. 水平地盤特性の入力および深度の関係から、水平地盤ばね値およびばね耐力を算出します。
    • 1.2. 杭断面の入力から、杭の剛性、耐力および骨格曲線を算出します。
    • 1.3. 基礎梁の入力から、つなぎ梁の剛性を算出します。
    1. 群杭フレームモデル
    • 2.1. 群杭フレームモデル(基礎梁剛)により、慣性力解析を行います。
    • 2.2. 群杭フレームモデル(基礎梁剛)により、地盤変位解析を行います。
    • 2.3. 慣性力、地盤変位結果を重ね合わせ検定応力を算出します。
    • 2.4. 検定応力を用いて断面検定を行います。

同時載荷の場合

    1. 共通
    • 1.1. 水平地盤特性の入力および深度の関係から、水平地盤ばね値およびばね耐力を算出します。
    • 1.2. 杭断面の入力から、杭の剛性、耐力および骨格曲線を算出します。
    • 1.3. 基礎梁の入力から、つなぎ梁の剛性を算出します。
    1. 群杭フレームモデル
    • 2.1 群杭フレームモデル(基礎梁剛)により、慣性力+地盤変位同時載荷解析を行います。
    • 2.2 発生した応力により、検定を行います。

計算詳細

解析条件 - 杭

モデル化

以下の条件でモデル化します。

  • RC系の杭の場合、鉄筋剛性を考慮します。

  • 解析モデルにおいては、杭を1mピッチで分割します。

  • 地盤ばねは杭頭から0.1m下がった位置から配置されます(上部慣性力と杭頭せん断力が釣り合うようなモデル化とするため)。

  • 分割された要素において、断面は下図のように上部の断面を採用します。

弾塑性設定

弾塑性(Tri-Linear)を考慮する場合、以下に示す骨格曲線の各折点は以下のように計算されます。

曲げひび割れ耐力Mc : 耐力式により計算
曲げ降伏耐力Mu, 曲げ降伏曲率φu : 平面保持解析により計算

曲げひび割れ耐力Mc(第1折れ点)

曲げひび割れ耐力Mcは、長期軸力を用いて以下の計算により算出します。

$Mc=(0.56・√Fc+σ_L)・Ze$

Mc : 曲げひび割れ耐力(Nmm)
Fc : コンクリート基準強度(N/mm^2)
σL : 長期軸応力度(N/mm^2)
Ze : 断面係数(mm^3)

※ 2025年版 建築物の構造関係技術基準解説書参照。ただし、長期軸力の扱いが同資料では軸力N, 断面せいDに対してN・D/6となっているところを、本プログラムでは $σ_L・Ze$ とし、長期軸応力度は鉄筋を考慮した剛性による応力度としている。

曲げ終局耐力Mu(第2折れ点)

第2折れ点は、曲げ降伏時耐力Muとします。

曲げ降伏後剛性低下率 α3

曲げ降伏後の剛性低下率α3は初期剛性に対して0.001倍とします。

解析条件 - 地盤ばね

モデル化

杭の節点と地盤ばねの接続方法は、

杭の節点から上下に指定ピッチ(規定値は1m)の半分の範囲にある地層のばねが接続されます。地層が指定ピッチによる範囲をまたぐ場合、範囲の境界で分割された地層によるばねが当該節点に接続します。

また、範囲内に複数の地層がある場合、当該節点には複数のばねが接続されます。

各地層データの深度と杭頭位置の関係は、以下の4つの入力から計算されます。
構造芯よりも高い位置に設計GLがある場合が+となります。

A. 1FL構造心と設計GLとの距離
B. 最下FL構造心から杭頭までの距離デフォルト値
C. 最下FL構造心から杭頭までの距離個別指定値
D. 設計GLと地表面との距離

地盤変位の入力深度は地表面からの深度として解釈します。
解析モデル上の杭の分割位置(節点位置)に対して、入力した地盤変位を線形補間して解析上の強制変位を算出します。

地盤ばね値

基礎指針の方法による地盤ばね

日本建築学会「建築基礎構造設計指針」に基づき算出します。

A) 線形

水平地盤反力係数は以下の式で計算されます。
$k_{h0} = C_k・α・ξ・E_0・D ^ {-0.75}$

$k_{h0}$ : 基準水平地盤反力係数(kN/m^3)
$C_k$ : 剛性低減係数
$α$ : 評価法によって決まる定数(1/m) $ξ$ : 群杭の影響を考慮した係数
$E_0$ : 700N もしくは 直接入力 (kN/m^2)
$D$ : 杭径 (cm)

B) 非線形

水平地盤反力係数は以下の式で表され、塑性水平地盤反力Pyを考慮する場合はPyを応力度の上限とします。

(0.0 cm <= y <= 0.1 cm の場合)
$k_h = 3.16・k_{h0}$

(0.1 cm < y の場合)
$k_h = k_{h0}・y^{-1/2}$
$k_h・y ≦ p_y$ (※指定によりpyを考慮する場合のみ)

$k_h$ : 水平地盤反力係数(kN/m^2)
$y$ : 水平変位(cm)
$p_y$ : 単位面積当たりの塑性水平地盤反力(kN/m^2)

Francisの地盤ばね

日本建築学会「建物と地盤の動的相互作用を考慮した応答解析と耐震設計」を参考に以下のような計算で算出します。

A) 線形

水平地盤反力係数は以下の式で計算されます。
$kfs = Ck・\cfrac{1.3・Es}{ 1.0 - Vs^{2}}・(\cfrac{E_s・B^{4}}{E_p・I_p})^{1/12}$

$E_s = 2.0・(1.0+ν)・ρ・V_s^2$

$k_{fs}$ : 水平地盤反力係数(kN/m^2) $C_k$ : 剛性低減係数 $E_s$ : 地盤ヤング係数(kN/m^2)
$V_s$ : せん断波速度(m/s)
$B$ : 杭径(m)
$E_p$ : 杭ヤング係数(kN/m^2)
$I_p$ : 杭断面2次モーメント(m^4)

B) 非線形

水平地盤反力係数は双曲線モデルで表され、Bromsの耐力に漸近します。

$p(y) = \cfrac{k・y}{ 1.0 + \cfrac{k・y}{P_y} }$
$P_y = p_y・B$

$p(y)$ : 双曲線モデルの応力度
$k$ : 初期剛性(kN/m^2)
$y$ : 水平変位
$B$ : 杭径(m)
$P_y$ : Bromsによる塑性水平地盤反力(kN/m)

塑性水平地盤反力

砂質土、粘性土でそれぞれ以下に示す計算により単位面積当たりの塑性水平地盤反力を計算します。
計算値に対し、杭径Bと層厚l’ を乗じてばね耐力とします。

砂質土

$p_y = C_s・κ・Kp・σ_0$

$p_y$ : 塑性水平地盤反力(kN/m^2)
$C_s$ : 耐力低減係数
$κ$ : 群杭の影響を考慮した係数, 単杭の場合3.0(本プログラムでは常に3.0を採用)
$K_p$ : 受動土圧係数 = (1+sinφ)/(1-sinφ)
$φ$ : 内部摩擦角 = √(20・N)+15
$σ_0$ : 上載圧(kN/m^2), 地層の自重を上層から累積した値を用いる。地下水位を考慮した場合、地下水位以下の地層については単位体積重量から $9.8kN/m^{3}$ を減じた数値を用いて計算する。

粘性土

(z / B < 2.5の場合)
$p_y = C_s・2・(1.0 + \cfrac{μ・z}{B})・c_u$

(z / B >= 2.5の場合)
$p_y = C_s・λ・c_u$

$p_y$ : 塑性水平地盤反力(kN/m^2)
$C_s$ : 耐力低減係数
$μ,λ$ : 群杭の影響を考慮した係数, 単杭の場合 μ=1.4, λ=9.0(本プログラムでは常にμ=1.4,λ=9.0を採用)
$c_u$ : 非排水せん断強さ
$D$ : 杭径(m)
$z$ : 深さ(m)

解析条件 - つなぎ梁

つなぎ梁を考慮する場合、RESP-Dに入力された基礎梁について反曲点位置が梁の中央であると仮定し、材長の1/2までをモデル化します。

解析条件 - 外力

慣性力

慣性力は以下のように算出されます。

  • 短期は1次設計用外力、終局は2次設計用外力(終局検定用のステップ)における外力を用います。

  • 基礎震度は杭基礎解析条件の設定値とし、基礎重量は建物重量に加え、杭配置テーブルで設定した基礎重量を加算した重量を用います。

  • 直接入力の場合は入力した値をそのまま用います。

  • 単杭モデルの場合、各杭が負担する慣性力は杭頭変位を同一とした線形解析により算出します。単杭剛床モデルの場合、上記で計算した全慣性力を各杭の頂部節点を水平拘束した代表節点に対して作用させます。具体的には以下のような手順になります。

      1. 慣性力分配用の群杭フレームモデルによる解析を実行する
      1. 各杭頭位置におけるせん断力分担比により、各単杭モデルに対する慣性力を算出する

慣性力の分配計算

単杭モデルへの載荷

応答変位

応答変位は、ユーザーが設定した地盤変位の入力値に対し、地盤ばね位置において線形補完した値を強制変位として地盤ばね側節点に作用させます。

付加モーメント

付加モーメントは慣性力ケース、応答変位ケースの事前解析として載荷し、継続して各水平荷重ケースの解析を行います。 RESP-Dにおける付加モーメント自動計算機能は以下に従って算出されます。

考慮されるモーメントは以下の通りです。

1.支承材のP-δ効果によるモーメント
2.支承材のQ-hによるモーメント

以下、その他の事項を示します。

  • P-δ効果におけるPは、杭基礎解析条件において軸力の扱いについて「直接入力値を採用」としている場合には直接入力された軸力が考慮される。
  • δについては免震層設計条件に設定された「免震層層間変位」の入力値を採用する。
  • Q-hモーメントは応力解析におけるせん断力とは無関係に、「免震層層間変位」の入力値から支承材に生じる歪みを想定してQを算出することにより用いる。

断面検定

断面検定 - 応力の重ね合わせ

応力の重ね合わせは単純累加、SRSSから選択できます。

A) 単純累加

$Md = Mi±Mr$
$Qd = Qi±Qr$

$M_d$ : 重ね合わせモーメント
$Q_d$ : 重ね合わせせん断力
$M_i$ : 慣性力ケースによるモーメント
$Q_i$ : 慣性力ケースによるせん断力
$M_r$ : 応答変位ケースによるモーメント
$Q_r$ : 応答変位ケースによるせん断力

B) SRSS

$M_d = \sqrt{(M_i^{2}+M_r^{2})}$
$Q_d = \sqrt{(Q_i^{2}+Q_r^{2})}$

$M_d$ : 重ね合わせモーメント
$Q_d$ : 重ね合わせせん断力
$M_i$ : 慣性力ケースによるモーメント
$Q_i$ : 慣性力ケースによるせん断力
$M_r$ : 応答変位ケースによるモーメント
$Q_r$ : 応答変位ケースによるせん断力

断面検定 - 上下動係数・偏心の考慮

上下動係数を考慮した場合、検定用軸力に対して上下動係数×長期軸力を加算または減算した採用します。
偏心距離を設定した場合、偏心距離×軸力のモーメントを杭設計用応力に加算します。 考慮する軸力は、上下動係数を考慮します。

断面検定 - MN相関の計算方法

杭体コンクリートの許容応力度

杭体コンクリートの許容応力度は、2025年版建築物の構造関係技術基準解説書に準拠し、以下の表により算定されます。

種類長期短期
圧縮引張せん断圧縮引張せん断
掘削時に水または泥水を使用しない方法
または強度試験により
強度・寸法・形状を確認できる場合
1/4・Fc1/40・Fcかつ
3/4・(0.49+1/100・Fc)以下
長期に対する
値の2倍
長期に対する値の1.5倍
それ以外の場合1/4.5Fc
かつ6以下
1/45・Fcかつ
3/4・(0.49+1/100・Fc)以下
長期に対する
値の2倍
長期に対する値の1.5倍

短期許容曲げ耐力

短期許容曲げ耐力は、平面保持解析により、以下のいずれかが許容応力度に達した時点のモーメントとして算出します。

  • 鉄筋の引張
  • 鉄筋の圧縮
  • コンクリートの圧縮

ただし、ヤング係数比は実際の鉄筋のヤング係数とコンクリートのヤング係数の比率として考慮します。
以下にRC規準より抜粋した曲げ短期許容曲げモーメント算出の考え方を示します。ここで、Sn, Inは中立軸についての値であるので、中立軸位置xnによって変動します。実際の数値計算としてはxnを変数として2分法により(2)式において所定の軸力と釣り合う偏心距離eを算出し、(1)式により短期許容曲げモーメントを求めます。軸力については、(3)式により最小の軸力を採用します。

$M = N・e$ … (1)
$xn - \cfrac{D}{2} + e = \cfrac{In}{Sn}$ … (2)

(コンクリートで決まる場合)
$N1 = \cfrac{Sn}{xn}・fc$ … (3-a)

(圧縮鉄筋で決まる場合)
$N2 = \cfrac{Sn}{(n・xn-dc)}・rfc$ … (3-b)

(引張鉄筋で決まる場合)
$N3 = \cfrac{Sn}{n・(D - dt - xn)}・ft$ … (3-c)

$N = min(N1, N2, N3)$ … (3)

Sn  : 中立軸に関する有効等価断面の1次モーメント  
In  : 中立軸に関する有効等価断面の2次モーメント  
Ae  : 等価断面積  
g   : 等価断面の重心と断面の最大圧縮応力度を受ける位置との距離  
e   : 偏心距離  
xn  : 中立軸位置  

短期許容せん断耐力

※2026年4月時点ではSRC杭(内部鉄骨)のせん断検定は未対応となっており、SRC杭(内部鉄骨)のせん断耐力はRC杭のせん断耐力として算定されます。

鉄筋コンクリート杭

短期許容せん断耐力は杭断面を矩形置換して以下の式により計算します。

$Q_A = b・j・{fs + 0.5・ft・(pw-0.002) }$

  A 		: 杭断面積  
  dt 		: 鉄筋重心位置  
  b		    : 矩形置換した幅(=√A)  
  d		    : 矩形置換した有効せい(=b-dt)  
  j 		: 矩形置換した応力中心間距離(=7/8d)  
  fs 		: コンクリート短期せん断許容応力度  
  ft		: せん断補強筋短期許容応力度  
  pw 		: せん断補強筋比  
鋼管鉄筋コンクリート杭

鋼管部分のみの短期許容せん断耐力を用います。

$Q_A=_sA_s・F/\sqrt{3}$

  sAs     : 鋼管せん断断面積  
  F       : 鉄骨基準強度

曲げ終局耐力

曲げ終局耐力は平面保持解析により、コンクリート圧縮歪み0.003、もしくは鉄筋・鋼管の引張歪み0.01に達した時点として計算します。コンクリートの応力度-歪関係はe関数法を用いて以下の数式で定義します。ただし、耐力劣化は考慮せずε≧εyの場合はε=εyとして計算します。

$σ(ε) = 6.75・(e^{-0.812・(ε/εy)}- e^{-1.218(ε/εy)})・Fc$

εy : 基準化歪 (=0.002)  

せん断終局耐力

※2026年4月時点ではSRC杭(内部鉄骨)のせん断検定は未対応となっており、SRC杭(内部鉄骨)のせん断耐力はRC杭のせん断耐力として算定されます。

鉄筋コンクリート杭

場所打ちRC杭のせん断終局耐力は、計算条件より修正荒川式と塑性理論式が選択できる。 なお、せん断耐力算定に用いる軸力は上下動係数を考慮し、押し込みと引き抜きの小さい方の値を採用します。

(1) 修正荒川式

$Q_{su}=\lbrace \frac{k_{1}・pt^{0.23}(18+F_{c}’)}{M/(Q・d)+0.12}+k_{2}\sqrt{pw・σ_{wy}}+0.1σ_0 \rbrace・b・j$

$b$, $d$, $k_{2}$については、 i)ウルボン鉄筋以外のせん断補強筋 ii)高強度せん断補強筋 で異なる値となります。

  A 		: 杭断面積
  B 		: 杭径
  dt 		: 鉄筋重心位置(円形断面における引張縁から引張鉄筋重心までの距離)  
  b 		: i)の場合は正方形置換した幅(=√A), ii)の場合は長方形置換した短辺幅(=π/4*B)
  D 		: i)の場合は正方形置換したせい(=b)  , ii)の場合は長方形置換した長辺幅(=B)
  d 		: 有効せい(=D-dt)
  j 		: 矩形置換した応力中心間距離(=7/8d)
  k1 		: 0.053
  k2 		: i)の場合は0.85, ii)の場合は0.846
  M/(Q・d) 	: せん断スパン比(=3.0)  
  σ0 		: 軸方向応力度(設定された軸力から算出, 引き抜きの場合も負の値として考慮)  
  pw 		: せん断補強筋比(=aw/(D・x))で、ii)の場合は0.1%≦pw≦0.3%とします
  aw 		: せん断補強筋1組の断面積
  x 		: せん断補強筋ピッチ
  σwy 		: せん断補強筋のせん断強度
  Fc' 		: 杭のコンクリート強度(指定により、1.0xFc、長期圧縮許容応力度×3のいずれかが選択可能)  

(2) 塑性理論式

$Q_{su}=b・j_{t}・p_{w}・σ_{wy}+η・k_{1}・(1-k_{2})・b・D・ν・F_{c}'$

ただし、$Q_{su}≦ \frac{ν・F_{c}}{3}・b・j_{t}$     $k_{1}= \frac{\sqrt{(L/D)^{2}+1}-(L/D)}{2}$     $k_{2}= \frac{3・p_{w}・σ_{wy}}{ν・F_{c}}$     $ν=0.7- \frac{F_{c}}{200}$

$b$, $D$については、 i)ウルボン鉄筋以外のせん断補強筋 ii)高強度せん断補強筋 で異なる値となります。

  A 		: 杭断面積
  B 		: 杭径 
  b 		: i)の場合は正方形置換した幅(=√A), ii)の場合は長方形置換した短辺幅(=π/4*B)
  D 		: i)の場合は正方形置換したせい(=b)  , ii)の場合は長方形置換した長辺幅(=B)
  η 		: 寸法効果による補正係数ですが、現状は必ず1とします(寸法効果を考慮しません)
  jt 		: 引張側主筋および圧縮側主筋の重心間距離(=D-2dt)
  dt 		: 鉄筋重心位置(円形断面における引張縁から引張鉄筋重心までの距離) 
  pw 		: せん断補強筋比(=aw/(D・x))で、ii)の場合は0.1%≦pw≦0.3%とします
  aw 		: せん断補強筋1組の断面積
  x 		: せん断補強筋ピッチ
  σwy 		:せん断補強筋のせん断強度
  Fc' 		:杭のコンクリート強度(指定により、1.0xFc、長期圧縮許容応力度×3のいずれかが選択可能)
  L 		:アーチ機構の長さ
  ν 		:コンクリート圧縮強度の有効係数

本式は、日本建築学会:「鉄筋コンクリート造建物の終局強度型耐震設計指針・同解説」1990年版に記載されている柱および梁のせん断耐力式に対し、cotφ=1として導出したものです。ただし、k2についてはウルボン指針2024年版解説に記載の通り、靭性指針1999年版に示されているトラス機構の有効係数λを考慮し、一般的な杭の配筋がλ=0.7程度であることから、終局強度型設計指針に記載されているβの(1+(cotφ)^2)に相当する係数を3としています。

鋼管鉄筋コンクリート杭

鋼管部分のみの耐力を用います。

$Qsu=_sA_s・F/\sqrt{3}$

  sAs     : 鋼管せん断断面積  
  F       : 鉄骨基準強度