【論 文
1
UPO ;624.
131.
55 :624.
15 :624.
042.
7 日本 建築 学 会 構 造 系 論 文報告集 第 367 号・
昭和 61 年9 月基
本 物
理 量
の
相
互
関係
に
基
づ く
地 盤
と
基礎
の
動 的
相
互
作
用
解
析
正 会員 正 会 員 正 会 員吉
藤
藤
長
義
大
田
幽
谷
井
行
*信
**地
* * * 1.
序 地 盤と基 礎の動的相互作 用 問 題に は,
インピー
ダン スの実部と虚部
,
及び地 震 強 制 力 (Seismic
excitingfor−
ce
,
Driving force).
なる 3個の基本物理量が存在す る。 これ らの 量か ら,
コ ンプライアン ス はイン ピー
ダン スの 逆数と して,
基 礎入力 動は地 震 強 制 力 を インピー
ダンス で割っ た 量 と して求める ことが できる。
基 本物 理 量を解 析 する手 法の 1つ に, 動 的
Green
関 数 (点 加 振解,
リング状線
加 振 解 )を用い る境 界 積 分 方 程 式 法が あ る。
こ の手 法では,
特に半 無 限 弾 性 地 盤 を対 象とする場 合, 逸散 波動の効 果 を厳 密に評 価す ることが で き る。1965年に
Lysmer
瑠 〕は , 円 板 を 同 心 リングに分 割し, Reissnerの解3) を応用 し た リン グ内で接 触 圧 を一
定とす る垂 直 方 向加 振解を重ね合わ せ ること によ り,.
剛 板 仮 定 を満足 する円板 基礎の垂直 方 向コ ンプライア ン ス を求め ている。
ま た, 1976年Wong とLuco4 〕は,
基 礎 領 域 を 接 触 圧一
定の矩 形小要素に分 割し,
.Thomson
とKobori
の並 進 運 動の解5 〕・
6〕を重ね合わ せる こ とに よっ て接 触圧 を離散 的に求め,
矩 形 剛 基 礎の コ ンプライアン ス を評価 し てい る。
さ ら に,
実体 波に よ る基 礎 入 力 動の解 析7)も 行っ てい る。
’
これに対し 1980 年頃よ り, 地表面点加 振 解 をGreen 関 数とする積 分 方 程 式 を立て,
これ を有限 要素法 的に離 散化し て解く境 界 積分方 程 式 法の研 究 が 盛 んに行 わ れて い る。 川 瀬,
吉 田ら は, 成 層 地 盤 上の剛 基 礎 8} や複数基 礎の相 互 連 成 系9〕の インピー
ダン スお よび基 礎入力 動の 解 析を対象に こ の方 法の適 用 を 試み てい る。
ま た,
小 林】°)は,
地 表 面 点 加 振 解の計 算法 を改 良す るこ と に よ り 計 算 量 を短 縮し,
基 礎の柔 性 を考 慮したインピー
ダン ス 解 析11}を行っ てい る。Iguchi
とLucOlz
} も矩形 基 礎の 剛 性をパ ラ メー
タ とするよ り詳細な解 析 例を示して い る。 これらの研 究は,
境 界 積 分 方 程 式 法に よ る半解 析 的な離 散 化 手 法が動 的 相互作用 問 題 にお け る基本 物理 量 を 混合 本 論文 は参 考文献.
15)に 加筆し、
ま と め た もの であ る。
ホ 広 島大 学 助手・
工博 *1 広 島大 学 助教 授・
工 博*
*
*
広 島 大 学 大 学院 生・
工 修 {昭和6]年4月17日原 稿受理} 境 界 値 問 題の解 として総合 的に解析し得ること を示し た もの と言える。’
一
方,
地 中 点 加 振 解を 用い た境界
積分 方程 式法の埋 込 み基 礎へ の適用も, 1982 年 Maedai3 )に よ っ て試み られ て い る。 しか し, 埋込 み基礎に対 する境 界 積 分 方 程 式 法 の適 用 例は,Green
関 数の代わ り に基 本 解で ある動 的 Kelvin解 を用い たDominguezz5
)の研 究な どを 除き,
い ま だ極 めて少 ない。
こ の 要因 は地中に載 荷 点 を もつGreen
関 数が地 表 面 載 荷に よ る もの に比 して複雑にな る こ と,
また埋込み基 礎の解析で はGreen
関 数の変 位 解 ば か りでな く応 力 解の計 算 を も要 求さ れ, 計算量 が大幅 に増大す ること な どの点にある。
そこで, 本 論 文で は求め るべ き基本 物理 量の間に存在 す る相互関 係14)を境 界 積 分 方 程 式 法に利用 す ること に よっ て地 盤の応 力に関 する情 報 を.
直 接 用いずにすべて の 基本 物理 量 を求め る解 析 法 を提案する15 )。 境 界 積 分 方 程 式 法の定 式 化に は直 接 法と間接法が あ る。
変位と応 力の関 係 を 直接 定 式 化す る ものが直接法で あ る。 これに対し,
ソー
ス分 布 (加 振強度分布 )を介し て変 位と応 力 を別 個に表 現 するものが間接法である。 本 論で はGreen
関 数の応 力 解の 計 算 を 避け る た めに 間 接 法 を採 用す る。 最 初に,Green
関 数の変 位解とソー
ス分 布に よっ て表現さ れ る変位 場 を 導 く。 イン ピー
ダン ス解析にお け る放 射 場の変位は,
基 礎と地 盤の接 触 面に 密 着 条 件を導入し,
ソー
ス分 布 を 決 定 す れ ば求まる。 通 常の 間 接 法で は, さ らにGreen
関 数の応 力 解を計 算し,
これを決 定 されたソー
ス分 布に従っ て重ね合わ せ ること に よっ てインピー
ダンス を 求 め る。
本 論で はこの よ う な 計 算 量の新たな負 荷 を避け, 基本物 理 量の相 互 関係に基 づき,
ソー
ス分 布か ら直接ま た は間接 的な方 法に よりす べ て の基 本 物理 量 を求め る。 イン ピー
ダン ス の定 義 式と ソー
ス分布の関係を論じ,
イン ピー
ダン ス虚 部 と静 的 イ ンピー
ダン スがソー
ス分布か ら直接 求め られ ることを示 す。
次に,
インピー
ダン ス虚部か らHilbert
変換を用い て イン ピー
ダン ス実 部を求める計 算 法 を 示す。Hilbert
変 換の 計算を実行する に は上 述の イン ピー
ダン ス虚 部と静 的 イン ピー
ダン ス が 必要であ る が,
Hilbert変 換が無 限一 103 一
積 分で表さ れ てい る た め
,
イン ピー
ダン ス虚 部は全 振 動 数 領 域にわ たっ て計算さ れ ていな けれ ば な らな い。
しか し,
これ は 数 値 解 析 上 不 可能で あ る。
ところで,
基礎表面か ら放 出さ れ る進 行 波 成 分の基 礎 表面近傍に お け る挙 動は, 高 振 動 数 域で平 面 波 状の実 体 波 と な る’
6レ。
こ の こと は地 盤の質 量と実 体 波の伝 播 速度 の積を 粘 性 係 数 とす る Winkler 型の 三次 元 粘 性 ダン パー
に よっ て基 礎 が 支 持 されて いること を意 昧す る。 こ の モデル よ り計 算さ れ る基 礎の粘 性抵抗は高振動数域に お け るインピー
ダン ス虚 部の漸 近 解であ る。 そこで解 析 の対 象 とし た振 動 数 領 域におけるインピー
ダン ス虚部の 数値解とこの漸 近 解 を 中 間 振 動 数 領 域で滑ら か に接 続 し,Hilber
し変 換の積 分を実行す る。
これにより通 常の 数値 解析法では不 可 能であっ た全 振 動 数 領 域に わ た る イ ン ピー
ダンス の実 部と虚 部が算 定で き る。
ソー
ス分布は基 礎の内部領域 と外 部領域の各 領 域にお け る応 力の基 礎 表 面に おける不 連 続 性 を表 す量で あ る。 地 震 強 制 力は,LucoM
の示 し た イン ピー
ダン ス解析で 得ら れ る地 中応 力 と地 震 強 制 力との 関 係を表す式に,
こ の不 連 続 量 を導 入し,
地 中応 力の情 報を消 去す る と,ソー
ス分 布と 入反射場変位の 積と して表現 さ れ る。
最 後に以上の定 式 化を 地中リング状ね じれ線加 振解IS} を 用いた軸 対 称 埋 込み剛基礎の ね じ れ動 解析に適用す る。
まず半 球 剛 基 礎の解析を行い,
厳密解と比較す るこ とにより定 式 化の妥 当 性 を確 認 する。 次に円 筒 剛 基礎の 解析 を行い,
既 往の数 値 解と比 較 することに よ り提案す る解析 法の有 効 性 を 明らか にする。
な お本解 析 法は, 次の仮 定に立 脚 する。 1) 地 盤は半 無 限三次 元 線 形 弾 性 体で,
かつ 均 質等 方 性を有す る もの とす る。 2) 地 盤に材 料 減 衰は無い もの と す る。3
>地 盤のボアソ ン比は0.5
未満と す る。
4)解 析に用い る各 量の 時 間 因 子 をettdt (i:虚 数 単 位, ω :円振 動 数,t
:時 間 )とし, 以 下の記 述で はこれを省略す る。
2.
境 界 積 分 方 程 式 法による変 位 場 表 現 図一
1に示す よ う なP
点に おける i方 向 点 加 振によ るQ
点で の ノ方 向の 変 位 と 表 面 力 をそれぞれ Ut丿(P
; 図一
1 半 無限場の領 域 と境 界一
104
一
Q
),p“(P
;Q
)と表記す る。表 面 力の方 向 を 定める法 線は,
境 界か ら地 盤 領 域 V 内に向か う もの とする。 ま た,
方 向 を 表す添 字 i,
ノは,
任 意の三次 元 直交座 標系におい て定 義さ れ, 例え ばj
=
1,
2,
3とする。S
ノ上の 自 由 地 表 面 条 件と遠 方S.
にお け る波 動の進 行 波 条 件を考 慮し て,
外 部 領域 V の P 点に お け るi
方 向変 位Ut を 直 接 法によっ て表 現 する と次 式と な る19)。u・(・)一
厶
{恥 (・・Q
)蝋Q
)−
u ・(・・Q
)・・(Q
)}dS
{Q
)・
………・
………・
……・
……・
・
(1 ) 上式で,
添 字は総 和規 約に従う。
また,
uKQ )とρJ(Q
) は そ れ ぞ れ領 域y
におい て定義さ れ た地 盤と基 礎の接 触 面S
.上の変 位と表 面 力で あ る。一
方,
内 部 領 域 V’
で は次 式が成 立する。o−
f
,.
IP
・(P ・Q
)u;(Q
)−
u・」(P
・Q
)P;(Q
)ldS
(Q
}・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
〈2
) こ こ に,
u:・
,
p;は領域y ’
内で定義さ れ る量である。S
.上で Ui=
u2 と し,
(1 >式か ら (2)式 を 差し引 く とソー
ス分布 η∫を 用い た間 接 法によ る変 位 場の表 現 が え ら れ る。
蝋・)痂
・’(・・Q
}nJ(Q
}・S
〈・
…・
………・
(・) こ こ に.
、
ηJ(
Q
)一・
ρ1
(Q
)−
PJ(Q
}・
…・
……・
…・
………・
・
…
(4) (3)式 をP
点の座標で微分す ればひずみ や応 力の表 現が得られ るが,
本 論では そ れ らの情 報を用いずに解 析 を進め る。
3.
静 的 イン ピー
ダンス と イ ン ピー
ダン ス虚部 イン ピー
ダン スを求め る た めに,
(3
)式で点P
をSu
上に移 動 し,
Su
上にお け る基礎と地 盤の密着 条 件ttt(
P
)=uf(P
)=
鉱島(P
)ノt丿・
………・
……・
…・
…
(5 ) を 導入 す る。
こ こ に,
罐 は基 礎のi
方向加振変位,
uS は基礎が加 振モー
ドん で運 動する と きのi
方向 変位成 分。
ま た, モー
ドを表 す 添 字 ノは,
基 礎が剛 体で あれ ばj
; 1−
6と な る。
上 式に対 応する放 射 場の ソー
ス分 布を (4 )式 か ら次 の よ うに置 く。
ηt= ηuAJ ; (P 訂
一
ρb
)ん…・
一 …・
・
……・
・
…・
・
(6
) こ のと き, (3
)式は次 式 とな る。
ufJ(・)
一
五
… (・・Q
)・・(Q
)・・(Q
)………一
(・) 上 式はS.
上の ソー
ス分 布 η届Q
)を未知 関数と す る積 分 方 程 式である。
これを解く と放 射 場の ソー
ス分布が得 られ る。 インピー
ダンス Kw(ω)は次 式で定 義され る。
K
、(・)==k
・’(・)・i
・・“(・)一一
/
.
・
k
・kdS
…………・
…一 ・
………
(8 )F
k
t
’
→ m’
1n
I
図一
2 内 部 抵 抗と外 部 抵 抗 Ok MoC (8 )式に (6 )式の関係を用いる と,
K
・(・)一
惹
儡 ・一 嘱 )dS
・
・
………・
一
(・) 右 辺 第 2項は領 域 y’
の 内 部 インピー
ダン ス で ある。
(9>式 を 概 念 的に示 し た図一
2によれば, (9
)式は内 部 抵 抗[h’
一
ω27π’
]と外 部 抵 抗 [h
(=
h
。一
ω 2m 。)+ ‘ωc]の 総 和である全 インピー
ダン ス か ら内 部 インピー
ダンス を 差 し引い た ものが求 めるインピー
ダン スであること を表 して い る。
第2項は,
仮 定2)よ り実 数であrl
t さ らに ω=
・
Oの とき領 域 V’
が剛 体 移 動するこ とか ら零とな る。
したがっ て,
次 式が成 立す る。
・w(・)一
厶
・竃… [・ ・(・一
・)]・・…・
・
……・
…・
(1・)・・J(・〉
一
(1/・)惹
撒 [・・]・S ・
…・
・
…・
…・
…
(ll ) すなわ ち,
静 的 イン ピL ダ ン ス とイン ピー
ダン ス の虚 部は地 盤 応 力 を 求める ま で も な く, (7 )式で決定 さ れ る ソー
ス分 布よ り直ちに計算す ること が で き る。4.
イン ピー
ダン ス’
実 部 複 素 振Ut
数 z=
ω十 iyの関 数H
(z)=HKz
)十iHl
(z) が次の条 件1
)ー の下半 平面におい て正 則。
il
) 実軸上で有界。
iii
)z の下半 平面でlzl
→ 。 。 の と き H(z)→
0 を満足 する とき,
関 数の実 部と虚部の間に Hilbert変 換 が成 立する20〕。
HK
・・一一
1
・∬
鷙
}・…一 ・
一
・・2−
1・碗
一
争
∬
駕
d
飾…一 ・
……
(・2−
2) こ こ ・.
・f
はC
・… h・ の主 値 積 分 を表す.
インピー
ダン スか ら構 成 され る次の関 数hw
(ω)−
h
、,(0
)+it・lc
、J(ω)−
c、,(0
)} Hu(ω);
ω2………・
…・
…………・
「
………
(13
) は仮定1
)− 4
)の下で,Hilbert
変換が成立す る た めの 条件i
)〜
i
の を満足 する。
そ の理由は次の ようである。
仮 定 3 )よ り縦波 速度が 無 限 大 と な ら ない た め,
ω → ・ 。 で 勧 (ω)は有 限 値に収 束する。 ま た,
こ の とき後 述 するi よ うに Cw(ω)も有 限な漸 近解に収 束す る 。 仮 定2)に よ り 砺 (ω)
,
Cw (a))はと もに ω の偶関 数とな る の で, (13 ) 式の 分 子は ω=
0の近傍でω2以 上のオー
ダー
である。
し た がっ て,
関 数H
‘J(ω1
はω≡O
で有
界で ある。
また,
仮 定1)より,
h
“(ω), Ce」〔ω)は実 軸上の そ の他の点にお い て も有 界である。
以上か ら条件
iD
とi
の が満 足さ れ る。 さ らに, 仮 定4)よ り条件1
)が満 足さ れ る。・
(13)式に (12>式の関係を適用 す る と次式を得る
。
k
・J・・ト箏
・∬
、篶
、d
… +・・(・・・
−9…
9−・
−7・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(ユ4−
1>c“… 一
箏
・∬
騾
≡
筈
甼
)d
・・… u(O)・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(14−
2 ) (14−
1)式に よれば, 静 的 イン ピー
ダン スとインピー
ダ ン ス の虚 部か ら, インピー
ダンス の実 部を算 定する こと がで き る。
さらに,
(14−1
)式に恒等式・
・
c‘」(。 。)・
f
。da
,o=
・
O
(ω1
一
ω2} を付 加し, ω → 。 。 と す る と次の関 係が得ら れ る。
h
・〈O)一
・・j(・・)一
(・/・)∬
颱 )一
・ii(・)}d
・・
・
・
・
・
・
・
・
・
・
・
・
・
…
t−一
一
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(15) 上 式は, 振 動数
を横軸に とっ て描い たインピー
ダン スの 図に お ける実 部 と虚 部の幾 何 学 的関係を示唆してい る (図一
3)。
ま た,
上式は高振 動 数 域に お け る イン ピー
ダ ン ス の実 部 を推 定す る手 掛か りを 与え ている が,
この場 kil(0) ) ) oooo ( ( kc Cij (0> Klj(ω )_
}
”
,
一
一
呷
一
一
c.
、
(ω) 0 ω → oo 図一
3 インピー
ダン ス の実 部と虚 部の幾 何 学 的関係 ω + ooも
図一
4 インピー
ダン ス虚 部の漸 近 解モデル一 105 一
合
,
インピー
ダン ス の虚部に関し 全振動数 領 域にわ た る 解析 値が必 要と な る。 高 振 動 数 域で基 礎 表 面か ら放 出さ れ る進 行 波 成 分は基 礎 表 面 近 傍で平 面 波 状の実 体 波とな る16}。 す な わち,
基 礎 表 面の任 意 点に おける外 向き法 線 n とこれに直 交 す る2
本の接 線 s,
tに よっ て構 成 さ れ る局 所 座 標系 (s,
t,n)を考える と, 基 礎はそ の法 線 方 向に は P波, 接 線 方 向に はS
波の粘 性 ダンパー
に よっ て支 持さ れ る (図一4
)。
これ を 式で表す と,
Ps 垢 P,=
‘ω [C
〕u ;・
………・
………・
(16) Pn 嬬 こ こ に,咐
i
÷
1
ゑ
]
p,Vs,
Vp
は,
そ れ ぞ れ 基礎近傍地 盤の質量,
横 波速度,
縦 波 速 度であ る。 局 所 座標系の 変位は全体座標系の変位と座標 変換マ ト リク ス [T
]に よっ て, 次の よ うに書け る もの と す る。
e e 秘ε Ul uf=
[T
] 鑪・
…・
一 …・
・
一 ……・
……
(17) を e u” u3 この とき,
基 礎の運 動に よ る粘性抵 抗係数,
す な わち,
インピー
ダン ス の虚 部は (16)式 と (17)式 を (8>式 に代入 す ることによ り次の よ うに表 現さ れ る。
c’(・・)
・
=一
五
1
嚠 丁] T [c
][T
]1
・g
’ldS
・
……
(・8
) イン ピー
ダンス の虚 部か ら実 部へ のHiibert
変換は, 積 分 領 域 を インピー
ダン ス の虚 部の数 値解が得ら れ てい る第 1区間0
≦ω。≦ a)N と そ れ以 外の第2
区間 WN ≦ω。≦ 。。 の 2つ に分けて行われ る。 第 1区 間で は数 値 解を最 小二乗 法によっ て有 理 代 数関数に近 似す る か,
ま た は区 分的に連 続な代 数 多項 式 (3次ス プラ イン関数)で補間し, Gauss
−
Legendre 積 分 法に よっ て Hilbert変換の積分を行 う。 その際
,
Cauchy の主 値 積 分 周 りで は,
座 標 変換に よっ て被積分関 数の挙 動を滑ら か に す る操 作を施〔
ε 評一
106
一
図一
5 数 値 解と 漸 近解の接 続 漸 近 解 す。
第 2区 間で は,
高振動数域で の イン ピー
ダン ス虚 部 の漸 近 解 c‘丿(・ 。)と 数値解を接 続す る近 似 関 数 を 求め,
積 分 を解 析 的に実 行す る (図一5
)。
近 似 関 数は CCJが ω の偶関数であ るこ とを 考 慮し て,
次のよ う な漸近 展 開形 式で表 現する。
,、、(。)一
耄咢
・
……・
…・
・
…・
……・
…・
・
…・
一
(・9 ) n=
0 ω an は,
第 1 区間にお け る後 半 部分の数 値 解 と漸 近 解 Cw(。 。)を 用い て,
最 小二乗 法によ り決定する。
M は第 1区 間の数値解が漸近解に滑ら かに接 続す る最 小 値 を 選 択す る。
す な わち,
tON で の数 値 解 と漸 近 解 を結ぶ 直 線 の傾きの正負が WN に お ける数 値 解の接 線こ う配の正負 と一
致す る場 合にM =
1, そ う で ない場 合は M=
2とす る。
(19>式 を式 (14−
1)式に代 入 し,
第 2区 間に お け る 積 分を解析的に行うと インピー
ダン ス の実 部は次式 とな る。kw
…一
尉争
艨
粤
・・・
齢 {
。嘉
、1
・gl器
1
罕
器
1
・ ・
翫
。一
,。一
}
凪 弾 ・}
・
…・
………・
一 …一 …一
〔20) (20
)式で ω → 。 。 と す る とイン ピー
ダン ス実 部の漸 近 解に対する近似式が得られ る。
h
・,(・・)・
h,・(・}一
髪
[
∫
一
・“(ah)ld
…・
嵩
( an2n−
D
ω解一
1]
・
………・
…一
(21) 上 式は, (15)式に示 した関係の近 似 表 現で あ る。5.
地震 強 制 力 地震 強制力は次式で定 義され るIT)。
昨
ゑ
(瀬一
pk ・t
)dS …・
t…………・
…・
・
(22
) こ こ に, uf ,ρ1
は そ れ ぞ れ入反 射 場の変 位と表 面 力。 領 域V
’
の内 部 問 題と領 域V
の外 部 問 題のS
。上 にお け る変 位の連 続 性 u:t= uk’
= u 急………・
一 ………・
・
…・
(23) よ り (22 )式は次の よ うに書く こ ともで き る。
・・
−
fSu
(・鰭一
醐 ・・・
・
…………・
…・
…・
(・4 ) Su上における応 力の不 連 続 性を表す (6>式を 上式 に用い ると ,鳳
(pi・・k
’
一
鋼 …ひ
・{・・…・
・
(25 ) 上 式の右 辺 第1
項は,
領 域 γ’
における相 反 定 理によ り消え る ので次式を得る。
・・
−
X
.
… {dS
………・
一………一 …t・
…・
・
(・・〉す な わ ち
,
地 震 強制力はソー
ス分布と 入 反 射 場変 位の積 を基 礎表面 全 体にわ たっ て積 分し た もの と な る。 (26 ) 式は,
計算 上,
地中応ガに関す る情報が 必要と さ れず,
簡 便であ る。
.
.
6.
軸 対 称 剛 基 礎のね じれ 動 解 析 前節までに論じて き た本 解 析 法の有 効 性 を検 討す る た め,
半無限一
様 弾 性地盤の地 中リング状ね じれ線加振 解 を利 用し た軸
対 称 埋 込み剛 基 礎の ね じれ動 解析 を行う。
解析には,
円 筒 座 標 系 (r,g,
z) を用いる (図一
6)。
6.
1 解 析の準 備 軸 対 称 剛 基 礎のね じれ 動 問 題の変位 場は, ep成分の み と なり,
し かもg方 向に沿っ て変 化し ない。
そこで,
(7 ) 式で epに関す る積 分を実 行し て お く と,u$(・)
一
・・f
、. ・。 。(PIQ
)・,(Q
)・C
(Q
)…一
(27) こ こで,P
= (re,
z。},
Q
’=
(r,
z),
Cu は軸 対 称 基 礎の 母線, Upp は次式に示す地 中リング状ね じれ線 加 振 解で ある18)。
財 (re,
2e ;r,z)= 屹9〔γθ,
9ε;r,
z) 十t4ちP(re,
−
Ze ;r。
z)’
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(28
) こ こ に,
鳩P(r。,
・
z。;r,
z)= u Smp (re,
Ze ;r,
z) 十U3e(re,
Ze;r,
z)・
・
・
…
(29)晦 嘔 ・…
一
、毒
纛
[−
hK
・h
・’
十(2/h
)IK
(h
)− E
(k
)}]・
…・
…・
…・
一 …・
…・
・
…・
……
(30) π呂ワ(re,
ze ;7・
,
z)=
11
/(4
πμ)}∫ {
詣
exp ・一研
ア
1・−
Ze・)
−
exp (− qIz−
ze))IJt
〔qr}Ji
(qre }dq ・
…・
(31
)こ こ に
,
μ :地 盤の剪 断剛性, ノ= ωf
Vs,
J
, :第 1種 1 次 Bessel関 数, K :第 1種 完 全楕円積分,
E
:第2
種 完 全 楕 円 積 分,k2=4rr
。/[(r+r。) 1 +(z−
z。} t ]。蜘
SH WAVE 図一
6 解 析モデルと 座標系 (27 )式は次に説 明する離 散化手法に よっ て近似 的に 解かれる。 ま ず 母 線Cu
をN
個の直 線 要素Cn (n=1
,2,…,N
)に分 割 する。
各 要 素 上の ソー
ス分 布 を一
定 と し, そ の未 知 分 布 強 度を ηen とする。
そこ で, リン グ状 ね じれ線 加 振の分 布 位 置 を要 素の中 央 点Pn
に と ると, (27 >式は次のよ うな連 立 方程 式に置 換さ れ る。
u$胙 ・・
名∫
・轟Q
)・Pt;(Q
)姻 );m
=1
,2
,…
,
ハ1…
99・
・
・
・
・
…
一・
一・
・
・
…
(32 ) 上式 を解く と
,
ソー
ス分 布は区 分 的に連 続な柱 状 関 数と し て求め ら れ る。
これより,
静 的ね じ れ イン ピー
ダン スとね じ れイン ピー
ダン ス の虚 部は,
(27)式と同様 な操 作 を (10)式 と (11 )式に施し て次式の よ うに表さ れ る。
h
・・C
・)一
・・ゑ
垢・・[・・(・一
・)]・C …………
(・・)・TT(・)
一
(・・/・).
L
..
・$・m [・・]・C・
……・
・
……
(34) (33),
(34
) 式から (20)式に より ねじれインピー
ダ ン ス の実 部k
。r(ω)を計 算 することに な る が,
そ の と き (18)式か ら計 算 さ れる CTT(∞ )が 必要とな る。
半 径 R。 の半 球 剛 基 礎 と半 径R
。,
埋 込み深 さD
の円 筒 剛 基 礎の ね じ れインピー
ダン スK
. を無 次 元 化 振 動 数 a。(=
ωR
。/ V。)を用い て,Krr(aa}
;2
’
πμR
言lhT7
(α o)十iα oCrT (αo)卜……・
(35 ) と表 現 し た と き,
c。r (a。→ 。。)は次の よ うにな る。
た だし, 回 転 中 心は基 礎の中 心 軸と地 表 面の交点と する。
c・…e・一
牌
。、β1
二
搬
、。 筒、・
・
・
・
・… ま た,
ね じ れ 地震 強 制 力は (26)式より次 式と なる。F
・一
・・厶
…5
・C …・
一 ……・
…・
…・
…・
…
(・・) こ こ に, μ5
は入反 射 場の ね じれ動 成 分である。
入射 波 を平面SH
波とする と,Neumann
展 開 21 } を 用い て, u6=2uticos
(j2
cos θ’ )J
,(ゴrsin θ ’ )・
…
(38
) こ こ に,
u、は入射波 振 幅,
θ ’ は鉛直下方よ り計っ た 入 射角で あ る。6.
2
解析例 最初に,
本解析法の精度 を調べ る た め,
半球剛基 礎の ね じれ イン ピー
ダン ス と斜めSH
波入射を受け る場合の 基 礎入力 動の解 析 を 行う。
離 散 化には半 球の母 線 を16
等 分し た直 線 要 素 を 用い,
要 素 上の 積 分はGauss−
Legendre
積 分 法の 2点 積 分に よっ て行う。
た だし, 加 振 点 を含む要 素で は,
要 素 を2分し,
それ ぞ れにおける 積 分を リング状ね じ れ線 加 振 解の静 的Kelvin成 分に相 当する (30>式に対し て 26点,
そ の他の成 分に対 し て 2点 積 分で評 価す る。 (32 )式に よ りソー
ス分 布の 解 析 値を無 次 元 化 振 動 数 α。に関し て6.
0
まで求め て お き,
(33)〜
(36) 式お よ び (19)一
(21) 式の手 順に従っ て,一
107
一
ね じれインピ
ー
ダン ス を算定する。
これ を表一
1に示す。
表に よ れば, 解析 値と厳 密 解 22) に は, 極め て良 好な一
致 が 見られる。
また, (37 ), (38 )式に よっ て SH 波に よ る ね じ れ地 震 強 制 力 を計 算し, こ れ を インピー
ダン ス の解 析 値で割 るとね じれ 基礎入 力 動 鵡 が得 られ る。
その絶 対 値 を表一
2に示す。SH
波の入射 角は, 45°
で ある。
表 中 無 次 元化に用い たu。は 入射 波 振 幅の 2倍, す な わ ち, u。=
2Ui
であ る。
表よ り解 析 値は ほぽ完 全に厳 密 解zz }を再 現 し てい る。
次に, 円筒 剛 基 礎の解 析 結 果 を 示 す。 埋 込み パラメー
タは β。 (;D
/R
。) =1である。
要 素 数は,
円 筒の側 面,
お よ び底 盤の各 部 分で 8,
合 計16 とし,
隅 角 部の応 力 表一
1 半 球 剛 基 礎のね じれインピー
ダン ス の 解 析 値と厳密解の比較 脚 听 籔 厳 密 鯛 5ロ
kτ
匸
c匸
匸
k匸
rCrr
0.
OL9940,
0002、
0000.
OCO o、
5L 聞2o.
1331.
8670.
133 1.
o1.
6640.
3321、
6670333 2,
0L4590、
531 且.
467o.
533 3.
o1.
397o.
5胛 ■.
400o,
600 タ・
o1.
3臼3o.
魄 弓 ■.
珊 o,
527 5.
o1、
3α90、
639L3590.
6 朝 6.
01.
3唱0o,
5腸 1.
35ユ 0.
5四 go ■.
33DO、
6651.
3認 o,
667 kH=
Re [KTTt (IIUR.
り ],
tT.
=
1m【KTTノ(2n幽R.
ヨ
a.
)コ 表一
2 斜め SH 波入射 (θ’
=
4S”
)を受け る 半球剛基礎の ね じ れ基礎入力 動の解 析 値と厳密解の比較一
108
一
集 中を考 慮 し て不 均 等 分 割す る。 解 析の手 順は前 述の半 球 剛 基 礎の場 合 と まっ た く同じ である。
解 析 値とGreen
関 数 分 布 法による Apse1 の結 果2s)を表一
3で比 較し て い る。 表に記 載 されて いるLuco の解 24)は,
静 的ねじれイ ンピー
ダン ス の厳密 解であ る。 Apsel の解は,
地 盤に材 料 減 衰 を導入 して得ら れ た もの である の で,
低 振 動 数 域 (α・=0.
5
)に おけるイン ピー
ダン ス虚 部に関し, 本 解 析 値とApsel
の解に は大き な差が見 られ るが,
その外 表一一
3 円 筒 剛 基 礎のね じ れインピー
ダンス の 解 析 値と既 往の解の比較 4.
o 3.
o 2.
0 1.
0 瞬 析 懾 byA 卩6d 5●
に丁
粐
CrrkH ご TT 0.
03.
■〒o0.
0ひ03.
166 (b野
Luc9 ) o.
52.
90L0.
3■32、
0.
ら27 1.
o2.
556o.
7■幽 2.
5L90、
750 2.
02.
246Lo502.
2田 LO50 3.
o2 」79L.
1552,
128L.
1嗣 1.
02 」23 匸.
1聞 2,
lo9L.
18 且 5.
02.
n21.
21臨 2,
LO61,
197 8.
o2JH 」.
2222.
u51.
206 092.
0941.
248.
,
kTr冨
R鱈
[KH ノ(2nμRo自>1,
c叮=
【皿【Krr!(2nμR.
,
a日
}1 krTτ ]
一
/ /
’
\
“
}
,
!
Ct1〆
〆
〆
!
!
゜’
亀
DL 。 2
.
。 S.
。 、.
囗 5、
。 、.
。 B.
=
uR.
!v■
図一
7 円筒 剛 基 礎のね じ れインピー
ダン ス 0,
50 e,
z5 0,
00一
〇,
Z5 R[R
囓
△
Tワ
u.
】d
叮
細
丶 〆、
ノ
Kヒ
!
ノ
丶
\ _
/ぺ
. 【m【R,
AT■
’u,
ユ’
°’
‘宅
.
。1
.
。 、.
。 3、
。 、,
。 5.
。 、.
。 A・
=
け
R,
’
Vg 図一8
斜 めSH 波 入 射による 円 筒 剛基礎の ねじれ基 礎 入 力 動で は, 比 較 的 良 好な
一
致 を示して い る。 表 中の本 解 析 値 を図 示す る と図一
7と な る。 ま た,
参 考の た め斜め 45Q 入 射 するSH
波によるね じれ基礎入力 動を図一
8に示し て おく。 7.
結 論 基本 物理 量の相互関係を境界積 分方程 式法に応用す る と,
放 射 場に関 する内 部 問 題 と外 部 問 題の境 界 面におけ る応 力の食い違い量で あるソー
ス分 布と基 本 物 理 量の間 に次の関 係が成 立する。
1 > 静 的 インピー
ダン ス は ソー
ス分 布 と加 振 変 位モー
ドの 積 を 基 礎 と 地 盤の接 触 面 上で積 分し た量で あ る。
2
) イン ピー
ダン ス の虚 部は ソー
ス分 布の虚 部と加 振 変位モー
ドの積 を基 礎と地 盤の接 触 面 上で積 分し た 量で ある。
3) 地 震 強 制 力は ソー
ス分 布 と入 反 射 場 変 位の積 を基 礎と地 盤の接 触 面上で積 分し た量で あ る。 インピー
ダン ス の実部は,
静的イン ピー
ダン ス とイン ピー
ダン ス の虚部および解 析的に得ら れ るイン ピー
ダン ス虚 部の高 振 動 数 領 域に お け る漸 近 解を利 用す れば,Hi
旦bert
変 換に よっ て算 定す るこ と がで き る。
し た がっ て,
放 射 場の ソー
ス分 布を決 定す れば,
すべ て の基本物 理 量 を求 めること ができ る。 このよ うに,
基本物理量の 相 互 関 係 を境 界 積 分 方 程 式 法に利 用す る と,
地 盤の応 力 その もの に関する情 報 を一
切 必 要 とし ない解 析 法が成 立 する。
解 析 法の特 徴を列 挙 する と以 下の 2点であ る。 1) 動 的Green
関 数の応 力,
すな わ ち変 位の微 分値 を計 算 する必 要が なく,
計 算 量は通 常の境 界 積 分 方 程 式 法の約1/4で あ る。 地中リング上線 加 振 解を利 用する軸 対称問 題で は約 1/3と な る。
2) 離散化手法では一
般に不 可能な全振 動数領域に お け るイン ピー
ダン ス 関 数の 推 定がHilbert
変 換に よっ て可 能である。
以 上に要約し た解析法の妥 当性を検 証するた め半 無 限一
様地 盤に埋 込ま れ た軸 対称剛 基 礎のね じ れ動 解 析 を 行っ た。 基本物理 量の 解析値を既往の厳密解お よび 他の手
法に よ る数 値 解と比較し良 好な一
致を確認 し た。
な お,
本 研 究の一
部は,
昭 和 60年 度文部省科 学研 究 費 補 助 金 (奨 励 研 究 (A))に よっ たこ と を付 記する。 参考 文 献1)Lysmer
,
J
.
:Vertical Motion of Rigid Fo6tings,
thesispresented to Unlversity of Michigan
,
at Ann Arbor,
in1965
,
in partlal fulfillrnent of requirements for the degreeof Doctor of Philosophy
2)
Lysmer
,
」.
and Richart,
F.
E.
:Dynamic Response ofFootings tQ Vertical Loadings
,
Journal
of 出e SoilMechanics and Foundations Division
,
ASCE,
Vol.
92,
) 3 > 4 ) 5 ) 6 ) 7 ) 8 > 9 10) 11) 12) 13) 14) 15} 16) 17) No
.
SMI,
pp.
65−
91,
Jan.
,
1966Reissner
,
E.
:Stationare,
axialsymmetrischedurch
eine Schuttelnde Masse erregteSchwingungen
eineshomogenen elastischen Halbraumes
,
Ingenieur−
Archiv,WBand ,
pp,
381−
396,
1936Wong
,
H,
L.
and Luco,
J.
E.
:Dynamic Response ofRigid Foundations of Arbitrary Shape
,
InternationalJournal
ofEarthquake
EngineeTing and StructuralDynamics
,
Vol.
4,
pp.
579−
587,
1976 Thomson,
W.
T.
and Kobori, T
.
;Dynamical Com−
pliance of Rectangular Foundations on an Elastic Half
−
Space
,
Journal
Qf Applied Mechanics,
TTans.
ASME,
Vol
.
30,
Series E,
NQ.
4,
pp.
579−
584,
1963小堀 鐸二
,
南 井 良一
郎,
鈴 木 有,
日下 部 馨 :長 方 形 基 礎の Dynamical Ground Compliance(その 1)一
半 無 限 弾 性 地 盤の場 合一,
京 都 大 学 防 災 研 究 所 年 報,
第10号A,
PP
.
283−
314,
昭 和42,
3Weng
, H
.
L.
and Luco,
J.
E.
;Dynarnic Response ofRectangular Foundations to ObLiquely Incident Seismic
Waves
,
International Journal efEarthquake Engineeringand
Structural
Dynamics,
Vol.
6,
No,
1,
pp.
3−
16,
1978吉田
一
博,
川 瀬 博,
佐 藤 俊 明 ;境 界 要 素 法 を用いた剛 基 礎一
地 盤相互作用解析 (その 2).
日本建築 学 会 大 会 学 術 講 演 梗 概 集 (北 陸〉, pp.
805−
806, 昭 和58年9月 川瀬 博,
吉田一
博 :境界 要素法 を用い た剛 基礎一
地 盤 相互作用解析 (その4),
日本 建築 学会大会学術 講演梗概 集 (関 東 ),
pp.
723−
724,
昭和 59年10月 小 林 俊 夫 :半 無 限 弾 性 体 地 表 面 点 加 振 解の無 限 積 分 を有 限 積 分で表す方法,
日本建 築学 会論 文報告 集,
第302,
pp.
29−
35,
1981 小 林 俊 夫 :Green関 数の離 散 化 手 法 を用い た建 屋と地 盤 の 動的 相互作 用の 研究,
日 本 建 築学会 論 文 報 告 集,
第 302,
pp.
2]−
27,
19811guchi
,
M.
and Luco,
J.
E.
:Dynamic Respense of Flexible Rectangular Foundation on a職Elastic Half・
Space
,
lnternationalJournal
of Earthquake Eagineeringand Stfuctural Dyna皿ics
,
Vol.
9,
pp.
239−
249,
1981Maeda
,
T.
:AStudy on the Dynamic Interaction Qf aRigid Embedded Foundation and Soil by the Boundary
Element Method in 3
−
Dimensions,
P【Qceedings Qf the 6thJapan
Earthquake Engineering Symposium,
、
pp.
1761−
1768,
1982 吉 田長 行, 藤 谷 義 信 :埋 込み構 造 物 と 地 盤の動 的 相 互 作 用 問 題に関す る基 本式,
構造 工学論文 集,
Vol.
31 B,
pp.
ll−
22,
1985.
3 吉田長 行,
藤 谷 義 信,
藤 井 大 地 :境 界 積 分 方 程 式 法によ る地 盤と基 礎の動 的 相互作 用 解析 (そ の 1》,
(その 2),
日本建 築 学 会 大 会 学 術 講 演 梗 概 集 (東 海 ),
pp、
215−
218,
昭和60年1b月Day
,
S,
M.
:Finite Element Analysis of Seismic Scat−
tering Prob【ems,
thesis presented to the University ofCalifornia
,
atSan
Diego,
Calif
.
,
in I977,
in partialfulfillment of the requirements for thedegree of Dectorof
Philosophy
Luco
,
].
E.
:Seismic Safety Margines Research Prog−
ram (PhaseI
),
Linear Soil・
Structure Interaction,
Lawrence LivermoreLaboratory, Califoinia,UCRL-
'oftheSeisrnologicalSociety
ofAmerica, Vol.66, No,1,15272,
July,
1980 pp.109-IZ3,197618)
:Nftfi,
paArsts
:NttheasV:nct6V
)xtitlten"Nfi 23) Apsel,R.J. :Dynamic Green'sFunctionsforLayered
ee,
Mliti
I\atRN,
Vol,31B, pp,23-34, 1985.3 Media ancl Applicationsto Boundary-Va!ueProblem,
19) Eringen, A.C. and Suhubi, E.S. :Elastodynamics, PhD dissertation,Univ.of Califernia,SanDiego,La Vol.[:Linear Theory,Academic Press, New Yerk,
,
Jolla,
California,1979pp.430-437, 1975 24} Luco, J.E.:Tersion efaRigid CylinderEmbedded in
20)
IWE-sc
:MmsP\gOktocoX\eeut,
ptrava,
#
an ElasticHalfSpace,Journal
ef AppliedMechanics,
utSE,
pp.542-543, 1979Transactiens
of the ASME, Vol.43,pp.419-423,197621) Watson, G.N.:Theery ef
Bessel
Functions,
Carn-
25)Dominguez.
J.:ResponseofErnbeddedFoundationsto
bridgeUniversityPrgss,London, pp.22-23, 1922
Trayelling
Wayes,M.I.T.
Research
Report,
22) Luco,
J.E.:Torsional
Responseef StructuresforSH No.R78-24, CivilEngineeringDepartment,1978Waves :the case of bemispherical foundations, Bulletin
SYNOPSIS
UDC:624.131.55:624.15:624.042.7
DYNAMIC
SOIL-FOUNDATION
INTERACTION
ANALYSIS
BASED
ON
INTERRELATIONS
OF
FUNDAMENTAL
PHYSICAL
QUANTITIES
byDr.NAGAYUKI YOSHIDA, ResearchAssociateof shima
Univ.,
Dr.YOSHINOBU FUJITANI,Associate
Prof.of HiroshimaUniy,, and DAUI FUJ", Graduate
Studentof HiroshimaUniv., Members ofA.I.
J.
There
exist threefundamental
phy$icalquantitiesin
thedynamic
soil-foundationinteraction
problem.These
are therealandimaginary
partsofimpedance,
and theseismie excitingforce,
An integralequation formulationbased on the dynamic Green'sfunctionsforan elastic medium is
developed
for
use as a tool inanalyzing thefundamental
physicalquantities.A limitednumber of theevaluationfor
themare available at present,most of which are restricted to thecase of
flat
ioundations
because
incase of theembed-ded
foundation
the enormous computational time isneededfor
calculation of thedynamic
Green'sfunctions.
Itis the objective of thiswork to rerneve thelimitation
by
applyinginterrelations
of the physical quantitiesto theboundary
integral
equation method.
The
displacement
field
is
representedby
thedisplacement
component of theGreen's
functions
and the sourcedistribution
at theinterface
of soil andfoundations.
Bbth
of theimaginary
partof theimpedance
and the staticimpedance
aredirectly
obtainedby
the sourcedistribution.
The
real partof theimpedance
is
indirectly
calcu-lated
from
themby
the Hilberttransform.The
seismic excitingforce
isrepresented as the productof the sourcedistribution
and thefree-field
motion through use of the reciprocal theorem, Therefere, theevaluation of the physica!quantitiesreduces totheproblem ofdetermining
theintensity
of sourcedistribution.
The pioposed method resolves a prevailingshortcoming inthe analysis of the embedded foundation;the ne-cessity of calculation