• 検索結果がありません。

基本物理量の相互関係に基づく地盤と基礎の動的相互作用解析

N/A
N/A
Protected

Academic year: 2021

シェア "基本物理量の相互関係に基づく地盤と基礎の動的相互作用解析"

Copied!
8
0
0

読み込み中.... (全文を見る)

全文

(1)

【論  文

1

UPO ;624

131

55 :624

15 :624

042

7 日本 建築 学 会 構 造 系 論 文報告集 第 367 号

昭和 61 年9 月

本 物

理 量

関係

づ く

地 盤

基礎

動 的

正 会員 正 会 員 正 会 員

**

* * *   1

序  地 盤と基 礎の動的相互作 用 問 題に は

インピ

ダン ス

実部と虚部

及び地 震 強 制 力 (

Seismic

 exciting  

for−

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 一

(2)

積 分で表さ れ てい る た め

イン ピ

ダン ス虚 部は全 振 動 数 領 域にわ たっ て計算さ れ ていな けれ ば な らな い

しか し

これ は 数 値 解 析 上 不 可能で あ る

 ところで

基礎表面か ら放 出さ れ る進 行 波 成 分の基 礎 表面近傍に お け る挙 動は 高 振 動 数 域で平 面 波 状の実 体 波 と な る

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 )

(3)

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 一

(4)

インピ

ダン ス の虚部に関し 全振動数 領 域にわ た る 解析 値が必 要と な る。   高 振 動 数 域で基 礎 表 面か ら放 出さ れ る進 行 波 成 分は基 礎 表 面 近 傍で平 面 波 状の実 体 波とな る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・

…・

(・・〉

(5)

す な わ ち

地 震 強制力はソ

ス分布と 入 反 射 場変 位の積 を基 礎表面 全 体にわ たっ て積 分し た もの と な る。 (26 ) 式は

計算 上

地中応ガに関す る情報が 必要と さ れず

簡 便であ る

     

       

 

6.

軸 対 称 剛 基 礎のね じれ 動 解 析   前節までに論じて き た本 解 析 法の有 効 性 を検 討す る た め

半無限

様 弾 性地盤の地 中リングね じれ加振 解 を利 用し た

対 称 埋 込み剛 基 礎の ね じれ動 解析 を行う

解析には

円 筒 座 標 系 (rg

 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

 Zer

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

(6)

ね じれインピ

ダン ス を算定する

これ を表

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

  〆

 〆

D

 

L 。 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)

で は 比 較 的 良 好な

致 を示して い る。 表 中の本 解 析 値 を図 示す る と図

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

 thesis

  presented to Unlversity of Michigan

 at Ann Arbor

 in

  1965

in partlal fulfillrnent of  requirements  for the degree

  of  Doctor of  Philosophy

2)

 

Lysmer

 and Richart

 F

E

:Dynamic Response of

  Footings tQ Vertical Loadings

 

Journal

 of 出e Soil

  Mechanics 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

1966

Reissner

  E

 :Stationare

  axialsymmetrische  

durch

eine Schuttelnde Masse erregte 

Schwingungen

 eines

homogenen elastischen  Halbraumes

 Ingenieur

Archiv

WBand ,

 pp

381

396

1936

Wong

 H

 L

 and Luco

 

J.

 E

:Dynamic Response  of

Rigid Foundations of Arbitrary Shape

 International

Journal

 of 

Earthquake

 EngineeTing and Structural

Dynamics

 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

3

Weng

, H

 L

 and  Luco

 

J.

E

;Dynarnic Response of

Rectangular Foundations  to ObLiquely Incident Seismic

Waves

 International Journal efEarthquake  Engineering

and 

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

 1981

1guchi

 M

 and Luco

 

J.

E

:Dynamic Respense  of Flexible Rectangular Foundation  on aElastic Half

Space

 lnternational 

Journal

 of Earthquake Eagineering

and Stfuctural Dyna皿ics

 Vol

9

 pp

239

249

1981

Maeda

 T

:AStudy on the Dynamic Interaction Qf a

Rigid Embedded Foundation and Soil by the Boundary

Element Method in 3

Dimensions

 P【Qceedings  Qf  the 6th 

Japan

 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 of

California

 at 

San

 Diego

 

Calif

 in I977

 in partial

fulfillment of the requirements  for thedegree of Dectorof

Philosophy

Luco

E

:Seismic Safety Margines Research Prog

ram  (Phase 

I

 Linear Soil

Structure Interaction

(8)

Lawrence LivermoreLaboratory, Califoinia,UCRL-

'oftheSeisrnologicalSociety

ofAmerica, Vol.66, No,1,

15272,

July,

1980 pp.109-IZ3,1976

18)

: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!ue

Problem,

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,1979

pp.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, 1979

Transactiens

of the ASME, Vol.43,pp.419-423,1976

21) 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,1978

Waves :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 three

fundamental

phy$icalquantities

in

the

dynamic

soil-foundation

interaction

problem.

These

are therealand

imaginary

partsof

impedance,

and theseismie exciting

force,

An integralequation formulationbased on the dynamic Green'sfunctionsforan elastic medium is

developed

for

use as a tool inanalyzing the

fundamental

physicalquantities.A limitednumber of theevaluation

for

them

are available at present,most of which are restricted to thecase of

flat

ioundations

because

incase of the

embed-ded

foundation

the enormous computational time isneeded

for

calculation of the

dynamic

Green's

functions.

Itis the objective of thiswork to rerneve the

limitation

by

applying

interrelations

of the physical quantitiesto the

boundary

integral

equation method.

The

displacement

field

is

represented

by

the

displacement

component of the

Green's

functions

and the source

distribution

at the

interface

of soil and

foundations.

Bbth

of the

imaginary

partof the

impedance

and the static

impedance

are

directly

obtained

by

the source

distribution.

The

real partof the

impedance

is

indirectly

calcu-lated

from

them

by

the Hilberttransform.

The

seismic exciting

force

isrepresented as the productof the source

distribution

and the

free-field

motion through use of the reciprocal theorem, Therefere, theevaluation of the physica!quantitiesreduces totheproblem of

determining

the

intensity

of source

distribution.

The pioposed method resolves a prevailingshortcoming inthe analysis of the embedded foundation;the ne-cessity of calculation

for

the stress components of the

dynamic

Green's

functions.

参照

関連したドキュメント

In particular, we consider a reverse Lee decomposition for the deformation gra- dient and we choose an appropriate state space in which one of the variables, characterizing the

This paper develops a recursion formula for the conditional moments of the area under the absolute value of Brownian bridge given the local time at 0.. The method of power series

Answering a question of de la Harpe and Bridson in the Kourovka Notebook, we build the explicit embeddings of the additive group of rational numbers Q in a finitely generated group

Then it follows immediately from a suitable version of “Hensel’s Lemma” [cf., e.g., the argument of [4], Lemma 2.1] that S may be obtained, as the notation suggests, as the m A

In our previous paper [Ban1], we explicitly calculated the p-adic polylogarithm sheaf on the projective line minus three points, and calculated its specializa- tions to the d-th

To derive a weak formulation of (1.1)–(1.8), we first assume that the functions v, p, θ and c are a classical solution of our problem. 33]) and substitute the Neumann boundary

In order to be able to apply the Cartan–K¨ ahler theorem to prove existence of solutions in the real-analytic category, one needs a stronger result than Proposition 2.3; one needs

Our method of proof can also be used to recover the rational homotopy of L K(2) S 0 as well as the chromatic splitting conjecture at primes p > 3 [16]; we only need to use the