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

上下地震動を受ける矩形平面状の屋根型偏平球殻の応答性状 : 応答値の推定と位相差の影響の検討

N/A
N/A
Protected

Academic year: 2021

シェア "上下地震動を受ける矩形平面状の屋根型偏平球殻の応答性状 : 応答値の推定と位相差の影響の検討"

Copied!
11
0
0

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

全文

(1)

【論   文】 UDC ;624

074

43 :624

042

7 :620

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

昭 和 63 年1月

上 下

地 震 動

け る

応 答性 状

答値

推定

位相差

影響

      正 会 員 正 会 員 正 会 員

西

 

* 紳 ホ   ホ   ホ

 

1.

序   構 造 物の耐 震 設 計に関 する構 造 規 定の多く は

主に水 平地 震 動を受け るもの を対 象とし て い る

しか し な が ら

地 動は本 来 3次 元 的で あり

直 下 型 地 震を受け る可 能 性 が ある場 合には

水 平 成 分に対して上下成分 も大き く な ると想さ れ る

そ の ため 上下 動の影 響が大きい と考 え られる大ス パ ン構 造 物の耐 震 設 計に際しては

上下 地 震 動 が構 造 物にお よぼす影 響 を究 明 すること が重要と な る

 上下地 震 動 を受け る シエ ル構 造 物の挙 動 を分 析

検 討 し たもの とし て, 以 下の研 究があ げら れる

 真下

田中

原1 )は 片 持コ ノイドお よ び筒 形シエ ル 屋 根の静 的 応 力解 析 と動 的 応 答 解 析の比 較から

片 持シ エ ル の応 答 性 状の分 析 を行っ て い る

西 村

新 宮Z }は

上 下 動な ら びに水平動 地 震を受け る円 錐 形お よ び球 形シ エ ル につ いて, 動 的 応 答 解 析 を進 め

応 答にお よ ぼ す裾 は りの影 響 を 検 討 して い る

國 枝3L4 〕

5}は

球 形シエ ル を対象に 固有振 動 数

固 有モ

ドを 正確に求め る計 算 法 を示し, 上下動に対す る実 用 的な結 果 を 与えて い る。 加 藤

吉 川

横 尾fi) , 長 浜

加 藤

横 尾 7 )

球 形 シエ ル が

短周期成分のい上 下 地 震 動 を受けると

シエ ル 頂部に大き な加 速 度 応 答 が生じる こと を 明 らかに して い る

 

方, 地 震 動の入力 位 相 差が構 造 物に与え る影 響につ いて,近 年 盛ん に問題 として取り上 げ られて いる。松井

瀬 谷8} , 位 相 差の ある水 平 地 震 動を受け る リング状基 礎をもつ シエ ル構造の 応 答 性 状 を 検 討 して い る

近 藤

田 中9)

Paraboric 

Velaroidall

 

Shell

を対 象 位 相 差入力に対す る簡便で実用的な上 下 動 応 答 解 析 を提 案して いる

 そこ で論文で は, 大ス パ ン構 造 物の基 本 的なもの の

1

つ と考え ら れ る屋 根 型 偏 平 球 形シエ ル (矩 形 平 面)に つ いて その地 震 応答の基本性状を全般 的にす る た め, 以 下の よ う な検 討 を進め る

  (

i

) 位 相 差を有する

Vihite

 noise 上 下 動 加 速 度

   

入力を

け る場 合につ い て弾性 応 答解 析を行い

     入力位相差が応 答値に お よ ぼ す影 響につ い て検討      す る1ω

  (

ii

) TAFT

  EL

−CENTRO ,

 HACHINOHE の 3

     つ の 上 下地震動に よ る弾性 応 答 解 析 をパラ メ ト      リッ ク に行い

応 答 値の中で も設 計上 特に問 題に      な る と考え ら れ る直 応 力と曲 げモ

メ ン トの応 答      性 状を分 析 する

 上記の検 討

分 析に基づい て

上下 地 震 動に対する軸 力および曲 げモ

メ ン ト応 答の最 大 値 を推 定 する略 算 式 を求め る。  

2.

振 動 方 程 式の誘 導  本 論 文で は

解 析モ デル と し て 図

一1

に示す剛な は り 上の 4辺で単 純 支 持 されて いる屋 根 型 偏 平 球 形シエ ルを 取り扱 う

x  y z 軸 方 向の変位をそ れ ぞ れ u

 v

ω

伸びひずみ をex, ag せ ん断ひずみ を γ, 曲げひずみ を Kx

 Ky

 Kxy

シエ ル 面の曲 率 を 屍

捩 りを kxy で表す と, ひずみ

変 位 関 係 式次 式え ら れ る1°} 。 日本建 築 学会大 会 (昭 和61年 )に て

部 発 表  . 技 術 科 学 大 学 大 学 院 生 1* 橋 技 術 科 学 大学   大 学 院 * * * 橋 技 術 科 学大 学   教 授

  〔昭 和 62 年 6 月 25日原 稿受 理 ) 靴

一h

 Ey

9tl

h

  ∂u         ∂v γ

∂y + ∂x

2     ∂2w     ∂2ω       ∂!w

Kx=

∂x ・

・Ky=

K・・

xOy

1

偏 平 球 形シエ ル を対 象と して い る か ら

      妬

=1

R

hmy≡

O であ る

 偏平 シエ ルに生 ずる 面内 応力 を

Nx,

 

Ny,

 Nxy

げ 応 力 を

Mx

,脇, 

Mxy

, 面 外せ ん断 力 を

Qx

, 

Qy

で表し, 応 カ

ひずみ関 係式 を次式の よ うに仮 定す る。 微 小 部 分に 作 用する各 応 力 成 分と単 位 体 積 当りの 荷 重 ρx

Py

ρ。 の 正 負の方 向 を 図

2a }

b)に示す

     ハJm

D且(ε n

十リ ε Ny

=D1

(εy十レεエ),

58

(2)

・       図

1a 解 析モ デル

At

一一一一一

一一

一一

     ’

        ’       ’ b    /     ノ     ’   ノ  A

      ●

R

      !     !     ノ   ’ 

  ’ ’

      X

  ,

R 図

1 b) 解 析モデル の幾 何 形 状    

Nxy

D

置(

1−

v)γ/

2

   

M

』= 

D2

Kx

+vKy ), 

My ‘’

Dz

Ky

+ vKx )

   M = 。

=D2

(1

Kxy

 

 

D

Kx

・猷 鮓 ・・

(K・+Ky)      

E

濃        

Es・

ts        

D1

=        D2=

     

1−

v2 t        

12

1一

ノ)      

 (2> こ こ で

Es

 p tsは そ れ ぞ れ ヤング係 数, ボア ソ ン比, シエ ル厚を表 す。  上下 地 震 動 鴨 が作 用 する場 合

偏 平 性 を考 慮 して

上下 地 震 動に よ る シエ ル面の慣 性 力の うち シエ ル の接 線 方 向の慣 性 力は無 視で きる とす れ ば, 次 式の平衡 方 程 式を得る

     ∂醗 ∂

Nxy

   

ax

+ ∂y

− h・’

Qx

=o

     ∂酌  ∂

Nxv

     

hv

Qs =

0

   

y 十 ∂x       ∂

Qx

   

k・

Nx+2北・ ・

晦 + κ・

Ng+

SiL

 

 

 

 

・(・・+

吻 一

……

3 こ こ で

ρは単 位 面 積 当り の質 量を表し

単 位 面 積 当た り の重量 を

M

, 重 力 加速 度をg と す る と

p

;M

g

で あ る

ま た

剛 な4 本の は り が

1あ るい は図

3 に示す よ うに 4で ピン接合さ れてい る場 合で

シエ ル周 辺で単 純 支持さ れ てい る場 合 を 解 析の対 象と す る の で

上下 地 震 動 を 次 式で仮 定する

z   y    Nx    Ny dx  Px d

 

亀      Nxy       Nx       X         Z      X   図

2 微 小 部 分に作 用す る応 力 A     D 図

3 位 相 差を有す る 上下動外乱の仮 定 X

  

 

w。

・w・

(x

一b

) + w・

  

  

 

・Wc

(x

穿

……・

……

(・ こ こ で 鵬,

W

, , 

vac

, 

JV

,は

4

隅点か ら 入力さ れ る 上下 地 震 動である

4辺で の境 界 条件

    v

ω

0 (at t=

0,

α

 u= w = 0 (at y

0

 

b

を満 足す る よ う に, 変位を ダ ブル

リエ 表示 して次 式 を得る

   

U

・  … S

・血」

   

v

・一

i

・・s

一・

(・)

   

一・

・i・

こ こ で

πmn

 Vmn

ωmn は そ れ ぞ れ (肌

 n)次

ドに 対す るx

,y

お よ び2 方 向の変 位 を表す

 弾 性 応 答 解 析に用いる基本式 を次 式よ り導 出す る。

 

 

 

a

1−

fJ

∂Nx  ∂Nxy

ax

ay

κ

・’

Qx

・…

dy

 

 

 

 

 

∂Ny  ∂ハlxy ∂y + ∂x

κ

Qy

 

 

 

 

 

f

L

lh

Nx・・

h

・y

N

・e ・

     

・働+

勘 蜘                      

e

…・

…・

………・

 

一・

 (6 ) た だ し

    

m

 

 

 

 

m

   δU

Σ Σ 

aUmn・

COS

        M11n

1

    

co

 

 

 

o

   δt丿

Σ]ΣユδVn,n

sin         M

In

1

    

tO

 

 

 

m    δω= ΣユΣ]δtOmn

sin         銅

1n

1 mnX    

 nπy

 

α

Sln

5

… S

7

in

     

59

(3)

(6)式 より

偏 平シエル の 剛 性トリス を導 出し

x

y 方 向の変 位 u

 v をz 方 向の変位 w で表し た後

減 衰 項を付加 す る と

1質 点系の 振 動 方 程 式 と等価 な 次 式を得る。

1 定 常white noise 入 力 時の解 析に仮 定し た諸数値

Span−1ength

a

b冨

99

64(m

Radius of  curvatureR

  99

64m

ThickRes5 ts=  1

25m Young

s

m。dulus Es

3

23x103 kgfcm2 Poisson ,s ratio V

    1/3 Weight peτ unit surface M=

 

40,90 ,160 (kgf/m2 〕 Damping factor

hmn言

0

02 Power  spectral densities S1=S2=

 

9c皿2/sec3

Cross poweτ 5pectral

densities

S12冨 

9

 

0r

 

O〔cm2 /sec3 Nu皿ber of sample waves

 

  

  

 

  

  

 

  

1 300 N  ber of 皿ode5 m葛 1

3

n=1

3 Non

dimellsional circular  frequency Ω11

  0

94 (φ

0 独立

φ

1 た だ し 勿

n十

2

んmn

ω

ガ 勿

十ω翫

ZVmn

− −

4

+(

− 1

)・ ・

+(

− 1

)・ C       πしηπ   十(

− 1

+nw 厨

……一 ・

……・

……一 …・

(8)        

K3=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Lt

・+

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

       λ

b/α        (

8

)式に対 する固 有 値 解 析の結 果

τ (

開 角 半 径   表

2 定 常white  noise に よ る解析結果 同

一,

(  )内は変動係数 %

Mxyは端 部

その他は中 央 部での最 大値 〉 ft}mn 一

・一

K

,・跳 ) κ1

lh

弖十2レhxky十

hi

十2(1

レ)

hk

lab

 

 

t

3

・+

9LM

Ln2

1

K

[(

kx

+v

hy

励 π

i

hy

+レ

ηαπ]      

1

T11 (sec  M (

kgflm2

)  W 〔cm   宙 (

〕   赫 〔 、

1

、)  Mx 

(t

m)

 

My 〔t

m)  Mxγ (ti

m) Nx (

) Ny 〔

) T11

02

φ

00

174 (23

3) 5

49 (23

9)  236 〔

18.

2

〕 0

0399 (2L4 )

0 .

0399 (

21 .

4)

0 .

0224

〔15

7) 0

00352

〔23

4)

0.

00352

〔23

4) M 冨 40 φ

10

242 (20

6〕 7

65 〔20

0〕   331 〔16

0) 0

055ア 〔20

6) 0

0557 〔20

6) 0

0274 〔15

9〕 0

00490 〔20

6) 0

00490 〔20

6〕 TlrO

5 φ

00

291 (24

6) 6

14 (24

6)   201 〔16

8

) 0

0675 〔24

0

) 0

0675 (24

0) 0

0386 〔16r9) 0

00590 〔24

6

〕 0

00590 (24

6

) M

90 φ二10

406 〔25

4〕 8

58 〔25

5)

 280 〔17

0) 0

0950 〔24

0) 0

0950 (24

0) 0

0478 〔17

7〕 0

00822 (25

3) 0

00822 (25

3〕 T11噐0

4 φ =00

408

〔26

0) 6

4

ア (26

1〕  

180

〔15

5) 0

0955 〔26

0) 0

0955 (26

0) 0

0563 〔18

8〕 0

00827 (26

0〕 0

00827 〔26

0) M

160 φ鵠10

572 〔27

0) 9

09 (27

4)  255 〔15

5) 0

133

(26

6) 0

133 〔26

6) 0

0695 (18

7) O』 116 (27

2) 0

0116 〔27

2) (Tt

入力位 相差

( 表

3 定 常white noise による解 析 結 果 }内は変 動係 数 %

Mxy は端 部

その他は中 央部で の最 大 値) TllM (kgf/m2) 〔sec  w 〔cm 〕   歯 (

)   w   cm 〔    sec 2) 呶 t

cm 〔 cm 〕  My 〔t

m〕  崎(y (t

cm 〕   cm Nx (tcm) Nγ 〔

τt

1T110

(25

386

2) 8

10 〔

25.

8)  232 (19

2〕 0

0895 〔23

9

0.

0895

23.

9)

0

0462 (17

4) 0

00782 (24

9

0 .

00782

〔24

9) T11

0

3 τt

25T110

302 〔24

9) 6

25 〔25

9)   200 (17

5 0

0674 (23

9) 0

0674 〔23

9) 0

0394 (17

4) O

00611 〔24

9〕 0

00611 (24

9〕 M

90 τt

5T11

0528 〔14

6)

1.11

(13

6)   128 〔9

69) 0

0159 〔18

0) 0

0159 (18

0) 0

0292 (15

9〕 0

00107 (14

6〕 0

00107 〔14

6)

一 60 一

(4)

Ω11L4L2 1

0o

8 0

6 o

4O

20

0 λ

1

0 一 一       r

1

4 − −       r

1

2 − −        r

1

O o

一一

一一

〇                               r

O

8 Cl

−−

o

− −CNH

,        r

0

6 9111

41

2 1

00

8 0

6 0

4 0

2 λ

0

7 一 一        r

1

4 0

一・

cwo

−・

o       r

1

2 0

−・

一・

(>        r

1

O

8                          r

0

6                                    0

0 40  60  80  100  120  140 160  180  τ        aλ

1

0        図

4 無 次 元

次 固 有円振動 数 Ωll         〔r

a/R

τ

R/te

λ

b/a) 40  60  80  100  120  工40 160  180  τ          b)λ

0

7   

5

  

4

 

53

〕 自 

2

の   

1

  

0

  

0 .0

54321

〔 巨 o

q の

00

0

0 ・

2

  

0 。4

  

0 ●

6

 

0

8

   Period

 

T

sec

1

0

TAFT

 

UD

 motion

0

5432

1

00

    〔 霞 O )

8

0 ,2

   

0 .4

   

0

6

 

0

8

   

Period

 

T

 

〔sec

1

0

HACHINOHE

 

UD

 motion 図

5

0 .2 .

 

O .4

   

0 .6

 

0 .8

   

1 .

O

    

Period

 

T

(sec ) 使用し た地 震動の変位応答ス ペ ク トル (入 力 加 速 度100gal

減 衰 定 数h

G

02} 〔 。

) 〉 の

0482

7 」

211

6

304

   

8

   

2

   

 

6

2

   

1

     

1

8

5

〕 〉 の 〔 oo の 丶 目 り ) 〉 の

iL−

CENTRO

 

UD

 motion

0

   

0 .2

   

0

4

   

0 ●6

   

0 .

8

   

1

Period

 

T (

sec

1 。

0

TAFT

 

UD

 lnotion

0 .

0

048

21

12

6

0 .

2

   

0 .4

   

0 .6

   

0 .8

   Period

 

T 〔

sec

1 .0

HACHINOHE

 

UD

      

雷.

   motlon  

0

0

  

0 。2

 

0

4

  

0 。6

  

0 .

8

  

1 .O

      Period  T 

(sec

6 使 用し た地 震 動の速 度 応 答ス ペク トル     (入力 最 大 加 速 度 コOO gai

減 衰 定 数 h

0

02

一 61 一

(5)

R

/シエ ル厚

ts

)にす る無 次元 固有円振 動 数 Ω を 次 式の よ う に置いた

      aJmn                      

 

tt・

tt・

9

>      

9mn=

      Es

 ts           (1

レ 2 )

α

b・

ρ こ こ に tOinn は偏 平シエ ル の固有 円 振 動 数で あ る

無 次 元

次固有円振 動 数 Ω11を 図

4に示 し て い る。  

3.

位相 差 を有 す る上 下 地 震 動に よ る応 答 解 析

 

偏平シエル の 4隅 支 持 点

A ,B ,

 

C

 

D

(図

1)に 用す る上下地 震 動は, それぞれ 独 立の場 合 も考えられ る が

こ こ では地 震 波が x 方 向に進 行 する と仮 定 し

上 下地 震 動 を 次 式の よ うに仮 定する

    

FF

琳 = 既 , 

F

=Ws =

…・

…・

…一 一 ・

(10) ま た

上 下動の性 質を統 計 的に取 り扱う た め,

F

、,

F

, を次 式の 相 関 特 性 を有 す る 平均 値 零の 定常white noisei3 〕 とする。 E [FI (t)

 

F

,(8)コ

=S

δ(彦

8)

E

Fi

(の

 

F

s

s2・

δ

t−

s)

E

F

、(

t

.F

,(s)]

;E

F

、(8),

F

、{

t

)]      

=Sn ’

δ(

t−

s)

11 こ こ で

,E

]は期待値

δ(

)は

Dirac

の デル タ関 数

S

,S2

はパ ワ

ー ・

ス ペ ク トル密 度, 

Slt

は クロ ス

パ ワ

ー ・

ス ペ ク トル密度を表す

 外 乱 F、と F,が

数 値 的に同

な 場合と, 独 立な場 合につ い て

それ ぞ れ

8

)式に対す る弾性応 答 解 析 を 行った結 果 を表

2に示す。 F,

 

F

,の発生はモ ン テ カ ル ロ

シ ミュ レ

シ ョ ンを用い た。 ま た

変 位あ るいは 応 力の応 答量は

ま ず, 平 均 加 速 度 法に よ り各 時 刻毎に denn

π

 tVmnまた

対 応する Umn

η  を求め

次に (

5

)お よび (2 )式 を用い て シエ ル の平 面上の 座標毎 に時刻歴 応 答 量 を 求める。 そして

これ ら の内の最 大 値 を求めて応 答 性状を 分析した

こ の時 使 用し たデ

タ を 表

一1

に示す

 

な お

一1

の ヤング係 数および シエ ル厚は

重層 ト ラス シエ ル を仮 定して おり

,.

下 弦 材 間 隔

d =

t。/

VS

72cmコア サ ン ドイッチ板を連 続 体 置 換し て求め て い る12〕 。 ま た

パ ワ

ー ・

スペ ク トル 密 度は

white noise の加 速 度の平 均 最 大値がほ ぼ 100 gal に な るよ う に決めてい る。   表

一2

よ り

,3

つ のケ

ス とも外 乱 F1 と 凡 を同

と し た 場合の応 答 値が

独 立とし た場合の応 答 値の ほ ぼ 厘 倍に なっ て いる。 こ の 結 果は

文 献 (14)の は りに 対する同様な解 析か らも知ら れて いる

 次に

外 乱 F

1 とF2 を 同

と仮 定した場 合におい て

入 力 位 相 差が偏 平 シエ ル の応 答に お よ ぼ す影 響につ い て 検 討 する。 こ こ で

,F

,と

F

,の位 相差 を次の よ うに仮 定 す る。      

F2

t一

τt)

0

 if t〈τt(τt≧0>

 

 

−t…

 (12 )

一 62 一

式 中, rt は入力 時 間 差あ るいは 入 力位 相差を表す 12}

τt を 偏 平シエ ル の 1次 固 有 周 期

T

,,の 1/10, 1/4, 1/2と し た場 合の弾 性 応 答 解 析結果を表

一3

に示す

こ の時, モ デル は表

1に示す Tn

O

3(sec

 

M =

90 (

kgf

m2 の時の デ

タを 使 用し た。 表

3よ り

τ尸

TII

/2と し た場 合の応 答 値が最 も 小さ く

τt

=TII

/4と し た場合

外 乱

F

、 と

F2

を独 立と し た場 合の応 答 値に ほぼ

致 し てい る

外 乱

F

,と

F2

を 同

と し

入 力 位 相 差を考慮 し ない場 合 と較べ

入 力 位 相 差を考 慮し た時の応答 値のが すべ てさ くな っ て い る

 本 解 析よ り確認で き たことを以 下にPtす

 

i

) 外 乱

F

,と

F

, を数 値 的に 同

あるい は独 立と     して解 析 し た場合, 同

時の応 答 値が独 立 時の応     答 値の ほ ぼ   倍と な る

な お, 弾性は りの場 合

   

Gauss 型whit  noise

語 倍になる こ と       が 文 献 (14)に示さ れて いる。   (

ii

> 外乱

F

,と

Fz

を 同

と し 入力位 相差 を考慮 し     た 場合の応 答 値は

,F

 

F

,を 同

と し て 入力 位     相差 を考慮し ないの応答値よ り すべ て小さく       な る

  (

1

ii

 よ り

外 乱が数値的に同

で入力 位相差 を考 慮し ない場 合応 答が最 大と な り, こ の場合が上 下 地 震 動を受け る偏平シエルの安全性を判断す る目安と な る こ と が確認で き た

 

4

実 地 震 動 に よ る応 答 解 析  前 章で は

外 乱が数 値 的に同

入力 位 相 差を考慮 し な い場 合に

偏 平シエ ルOP

応 答値が最 大に な る こと が 確 認で き た

  し たがって

本 章の解 析にお ける偏 平シエ ル の 4隅 支 持 点か ら 入力さ れ る上 下 地 震 動は

すべ て入力 位 相 差を 考 慮しない数 値 的に同

なもの を用い る。 入 力 地 震 波と

して は

EL −CENTRO

 

UD

(1940)

 

TAFT

CALI ・

FORNIA

UD

(1952)

 

HACHINOHE

 

UD

1968の 3 つ のそれぞ れ特 性の異 なる地 震 動 を用いた。 それぞれの 地 震 波の変 位, 速 度 応 答ス ペ ク トル を 図

5

6 てい る

  図

一4

に示し た無 次 元

次 固 有 円 振 動 数 Ωll を用い て

単位 面 積当りの重 量

M

と曲 率半径と シエ ル 厚の比 τ (==  

R

ts

) と

  x 方 向 パ ン と曲 率 半径の比 r (

α/

R

)をパ ラメ

タとして,

Es =

3230 kgf/cm2

  ts

86

6cm と な るよ うな α

 

T

,, を選んで弾 性 応 答 解 析を 行っ た

こ こで

パ ラメ

タの値と して は

,M

を40

90

160kgf/m2 τ を40か ら 20間 隔で 180 まで, r を

0.

6,e.

 

S,1.

0,1.2,1.4

と し た

ま た シエ の平 面 形 状を表すパ ラ メ

タ λ(

 bα)は 1

O

0

7と設 定し た

その ほ か, 減 衰 定 数 妬

= 0

02, ボ アソ ン比 v= 1/3

地 震 波の入力 最 大 加 速 度 A

100gal とし た

採 用モ

ド数は

地 震 波の ピ

クとモデル の固 有 周 期を考 慮して,

(6)

b

2 b

5 00 b

2 b

5 0

64    0

59 00 図

7 b

2 b

5  

Of4

      ψ Wmax

 

〔cm  5

7       7

9        彑    2 437 448

OO

      a14      a!2

  

〔W’Wg m、x 〔ga1)       120 b

5       O    

Of4

     評2    0     a!4     》2 Nxm 。x (kg/・m)

   

M.皿。x 〔 kg

c「n/cm シエ ル平面上の最 大 応 答 値 分 布 (TAFT  UD , M

40kgf/m2

τ

120

 r

1

0

λ

LO

τ

R/ts

 r

α/R

 λ

b/α) 対 称モ

ドにつ い て m  n 共最大25 次迄と し た。 こ れ は曲げモ

メ ン ト応 答が きわめて高 次モ

ドに敏 感で あ ること に よ る

ま た

固 有 周 期が0

05秒 未 満の短 周 期部分につ い ては解 析 対 象 外とし た

 解 析 結 果の

として

M 竺

40 

kgf

m2 τ

=120

,  r

=1.0,

λ

1

0の 場 合の

TAFT

 

UD

に対す る変 位

絶 対 加 速 度

直 応 力お よ び曲 げモ

メ ン トの最 大 応 答 値 分 布が図

7に示 されて い る

これ ら は, 各時刻毎に 甑 η

inmn

 tVmn

 Umn

 Vmn を求めた後, シエ ル の平面 上の座標 毎に時 刻歴 応答を求め 各 座標毎に大値を抽 出し たも の であ る か ら

これら の生 起する時 刻は必ずし も

致し ない

 

本論 文では

これ らの応 答 値の中でも設計上特に問 題 に なる と考え ら れ る直 応 力と曲 げモ

メ ン トの大 値

Nr.

 

,1

%囮 x

max

,1

晦  、の 応 答 性 状につ い て の分 析

検 討を行う

な お

Nx.

  、

 Nymax

』晦

  , 

My.

ma

f は シエ ル の平 面 上の座標 毎の時 刻 歴 応 答か ら各 座 標 毎 大 応 答 値 を抽出し

更に その 中の最 大 値 をとっ たもの で あ る

 こ こで, 上 下 地 震 動に対 する直応 力 および 曲 げモ

メ ン トの

般 的 応 答 性 状 を分析す る た め

そ れ ぞ れ次 式の よ うに無 次 元 化を行う

     

     1

v2      

Nx .

max

     

Nr.

rm [        Es

ts        

………・

…t・

t

13

)      

      1

vt

   

N・

flkl

ts

°

N

max a 1

2 0

8 D

4 0

0 EL

CENTRO UD motlon        

 

 

 

 

 

「△

 

s ::

{ピ

 

7

  

M

40kgf皿2

7r

 M

90kgf /m2

−一

 M

160kgf /m2

λ

1

Or

1

O α 1

2 O

8 40  60  80  100 120 140 160 180 τ        a

4 0

0        40   60  80 a 1

2 0

8 O

4 0

0 100 120   140  160 180  τ b) 40  60  80  100 120  140  160 180 τ       c} 図

8 地震 波 毎の a

τ関 係 (λ

1

0

r

1

0     (M

単 位面積当り の 重 量

r

a/R, τ

R/t、)

 

 

 

一1

     

…・

…・

…・

14

 

 

 

砿 一 一

11

り弧 一 ・   無 次 元 化さ れ た直応 力および曲 げモ

メ ン トか ら更に 構造物の固有 周 期

開 角 半 径

ス パ ン長お よ び入 力地 震 波の影 響 を抽 出す る た め

直応 力お よ び曲げモ

メ ン ト に関す る係 数を次 式の よ うに定義する

      Aκ

      R

   

α・

s

。(

h,

 T,、

Ax)

i4

N ・max                         

 (15 )      

A

κ

      R

   

α y=

s

、(

h,T

、、

 

AK

i

 

 

N

max

一 63 一

(7)

β 0

6 0

4 0

2 0

0        40  60  80   100 120 140 160 180 τ β 0

6 0

4 0

2 0

040   60   80  100  120  140 160 180 τ        b) α L2 O

8 0

4 O

0 図

10 賑 1

2

O

8 40  60  80  100  120  140 160 180 τ a の全 解 析 結 果 及び平 均 値 i

標準 偏差 ∂〔λ

1

0) 0

4              O

0       40  

60

  

80

  100  120  140 160  180  τ a β 0

6 O

4 0

2                O

0       40  60  80  100 120 140 160 180  T       c

9 地 震 波 毎の β

τ関係 (λ

LO

 r= 1

0      (M

±

単 位 面 積 当りの重 量

r==a/R

τ

R/ts> Cty 1

2 0

8 0

4 0

0 λ

0

7       ay +σy

αy

0

641 σy

O

096 毳y

6y

1

OTAFT  UD motion △EL

CENTRO UD motion 口HACHINOHE UD motion 図

11 40  60  80  100  120  140 160  180 τ       b} a

ay の全 解 析 結 果 及び平 均 値 醗

標 準偏差δr

iy (λ

0

7)

flx

、、

M

一 角一

s

A,)

M

max

 

一・

 

16

 こ こ で

,S 。

h,

 

T、

Ax

)は 入力 地 震 波の最 大 加 速 度が Ax (gal)時の

次固有周期

T

,,の変 位 応 答スペ クトル の 値 (

5)で あり

,一

,SI

h,

 

TI1,Ak

)はス ペ ク トル強さを表し

速 度 応 答スペ ク トル

Sv

h,

 

T .

Ax)(図

6

よ う 定 義 す

   

s・(ん

嗚・

A

s

・仇

T ・

A

・)

dT …・

(17) 本論 文で は仮に

ax

 av を直応 力 推 定 係 数

βx

角 を 曲 げモ

メ ン ト推 定 係 数と呼ぶ こ とにす る

こ こ で β。

禽 の表 示に

Sl

h,

 

Tn ,

 

A

,)を用 いて い る が, これ は,

SI

の代りに

S

.(

h,

 T,,

Aκ)を用い た表 現で は

地 震 波 毎で呈す る βの性 状に比較 的 大き な ば らつ き が 生

曲げモ

メ ン ト応 答の 高 次モ

ドに対す る敏 感 性が十 分評価で き なか っ た ためである

SI は Sb と同 じ次 元 を も ち, (17)式の よ うに短 周 期 成 分の影 響が含 ま れて いる形で示さ れ る か ら

高 次モ

ドに対す る 地震 波の 影 響 が ある程 度 反 映 され る もの と 考え られ る

した が っ て βの表 示に は SI を採用す る こ と と し たe   図

8 に各地震 波に対す るα と τ の関 係が λ

LO

r

1

0の 場合に つ い て示さ れて い る

こ こ で

λ

1

0

一 64 − ・

(8)

の場合 a

ax

ay

β

βエ

βシと な る か ら

以 降 α

β の 表 現 を 用い る

8を み る と

各ケ

ス とも

M に よ る変 動は ほ と ん ど み ら れず, τ に関し てα が

定 値 を取る傾 向に ある。 ま た

8a )

c)を較べ れば

地 震 波に よる α の値の ば らつ き も 小 さい こ と が わ か る

これ らの傾 向は

r

O

6

0

8,1.

2,1.

4の場合につ い て も同 様であっ た

。一

9に は各地震 波 毎の βと τ の関 係を λ

1

0

r

1

0の 場 合に つ いて示 し た

こ れも

α 同 様

,M

に よ る変動 は余り み ら れず, そ れ ぞ れ τ に関し て右下 りの

定傾 向を示して いる

地 震 波 に よる ぱ らつ き も 小 さ く

16

式の表 示で曲 げモ

メ ン ト応 答を十 分 評 価でき た と もの と 思 わ れ る

r

0

6

0

8

L2

1

4の場合につ い ても同 様な傾 向 がみ られた

  ま た

λ

O

7の場合に つ いて も

ax

 ay

βτ

燭 それ ぞ れにつ い て同 様なこと が確 認で き た。  本 解 析結果の r ごとの平 均 値お よ び 分散 をま と めて 表

4に示し た

本 解 析か ら

以 トの こと を確認 し た

 (

i

) 直 応 力 推定係数 ax

 as は

地 震波お よび 重 量      M に ほ と ん ど関 係せず

τ に関 して

定 値 を取       る 傾 向に ある

 (

i

冂 曲げモ

メ ン ト推 定 係 数 βx

角 も同 様に, 地     震波

M に ほ と んど 関 係せず

τ に関 して は

    右 下 りの傾向を 示 す

  本解析で は, 短 周 期 範 囲に主 成 分 を有する EL

CEN

TRO

 

UD ,

比 較 的 長 周 期 範 囲に主 成 分 を 有する

HA −

CHINOHE

 

UD

お よびその両 者の 中間的な性 質 を もつ

TAFT

 

UD

性 質の 異な る 3つの 地 震 波 を外 乱 と して採用 してい る

この 3つ の地 震 波 を用い て応 答 解 析を行っ た結果

(15)

〔16)式に示し た係 数 α

βは 地 震波の性質に か か わ らず, お よ そ

傾 向 を 呈 す こ と が わか り

a

βに よる無 次 元 表 示が

直 応 力

曲 げモ

メ ン ト応 答の評 価に有 効で ある こと が確 認で き た

 

5.

応 答 値の推 定 法  本 章で は

章の解 析 結 果に基づ い て

ヒ下地 震動を 受ける偏 平シエ ル の応 答 値の う ち

直 応 力と曲げモ

メ ン トの最 大 応 答 値 を 推 定 することを試み る。   5

1

αx

ay を用い た最 大 直 応 力応 答の推 定 法   (13)

(15)式 より

上 ド地 震 動を受け る偏 平シエ ル の直 応 力の最 大 値 Nx

 

 Ny

  が次 式の よ うに表 示で き よ う

       

Es ・ts

  SD

h,

 T1 、

Ax)A

   

N・

1−

. ・

 

R

   

°

晦      

(18)            

E 。・ts

 

SD

h,T

 Ar)

A

   

乢一

1

. ・

’ .

R

−一一

ZR

αy 3つ の地 震 波 分の デ

タ を すべ て ま と め て , 図

10

11 に示し た。 こ の と き

λ= LO の場 合 α の 平均 値 ∂

0

809 を用い て

直応 力最 大 応 答 値

Nx _ ,

1覧皿 . の平 均 値

E

Nx、

max ]

E [N。

max ]は次 式の よ う に表せ る

     

E

[N

c

max ]

E [N鮎   ]              

Es 。

 

ts

 

S

。(ん

T,】

 

A

,)

A

       

ff

 

R

 

h

E

 XO

809      

 

tt・

 

tt・

 

 (

19

, λ

O

7の 場 合に つ い て は, d. 

 O

 989, δy

0

641 と な るの で

 

 

 

ll

e

. ・欄

 

 

 

E

匹 一 ]

1

・s

°(九

1・

A

×・

641        

 (20) とで き る

し たがっ て

(19}

(20}式に お い て

対 象 と する地 震 動の変位応 答スペ ク トル Sp(

h

, T,、,A。)を用 いれば

地 震 時 軸 力の最 大応 答 値が推 定で き よ う。  5

2

βx

角 を 用い た最 大 曲げモ

メ ン ト応 答の推 定 法   直応 力 応 答の場 合 と 同 様に

(14>

(ユ6}式 より

,.

ヒ 下 地 震 動 を 受 ける偏 平 シエ ル の曲げモ

メ ン トの最 大 値

Mx .

 

, 

My.

max は次 式の よ うに表せ る

      

E

Mx.

  =

      12(1

レ2)        

E 。

tk

M

μ

  【

      12(1

レ2)

SI

(ん, TII

As)

A

 

 

 

 

α

 

 

4

κ

β x s∬(九

TH

4κ)      

β。       α      14K

21 ) βx

角 につ い て

,3

つ の 地 震波分 の デ

タ を すべ て ま と め て

12

13に示 し た

前 章で述べ た よ う

βy は τ に つ い て右 下りのき を もっ て い る か ら, v に 関す る回 帰 直 線で β、,

防 を表現し た

図中の 直線お よ

破 線がそれに あた り

βx

fiv

につ い て τ ご とに均 値と標 準 偏 差を求め

これ らのか ら最小二法にっ て求め た回 帰 直線で あ る

ドシエ ル で は

般に軸 力に 対し て曲げモ

メ ン トの値 は そ れ 程 大 き く な ら ない こ と か ら

12

13 に し た 回帰直線の う ち, 各 τ に対 す るβx

β。の平均 値か ら求め た回 帰 直 線 式に よ り

曲 げモ

メ ン ト最 大 応 答の 平 均値

E

[Mx

 

 E [Ms

max ] が推 定で き よ う。 (21 )式に基づ いて, λ

1

0の場 合は,     

E

Mx.

max ]

=E

M

猷 mx ユ               

E

ε

雌 

SI

h,

 

T

[1

ん}

A

              l2(1

v2) α

A

,      

0

0017τ十 〇

598}

 (

22

) ま た

λ

O

7の場 合

・[

M

max ]

1

51

必亙・

£

     

r

− 0,

0017

τ十 〇

569) ・ [My

max ]

藷葺

・s

” ん )

      《

0

0015τ十 〇

470)

 

9・

 (23

65

iL

(9)

 表 4 実 地 震 動による解 析結果 (E (

):平 均 値

σ (

):標準 偏差) λ = 1

0 λ

0

7 earthquakefE (α σα Eβ σβ〕 E〔αx σαx)E〔αyσαyEβx)σ〔βx)E(βy)σβ

y

0

60

807 0

1020

3680

0811

0360

0940

585 0

0960

5710

1320

3190

144 0

80

7870

0790

3930

074LO570

1020

5820

1020

3460

0930

3020

083 TAFT1

00

8080

0870 .

4310 .

0830 .

9950 .

1090

6580 .

1040

3720

0780

2900

051 1

20

8150

0700

4330

0950

9770

0740

6690

0870

4000

0910

3180

069 1

40

8090

062

0

4450

1310

9360

0750

6850

0810

4180

0910

3410

067 0

60

8390

0940

3110

0641

0350

0650

634 0

0920

3ユ

20 .

0940 .

2540

093 0

80

8510

1080

3800

0771

0450

1000

6430

10S0

3320

0680

2660

056 EL

CENTRO1

00

8190

0650

4080 .

0720

9980

0880

6630

0830

3640

0610

3000

056 1

20

8210

0810

4230

0850

9720

1040

6840

0880

3950

077O

3000

061 1

40

8360

0700

4270

0900

9480

0730

7190

0740

4090

0730

3Z70

064

0

60

8170

1070

3430

0731

0350

0860

5950 .

1200

3510

0980

3080

111 0

80

7890

0730

3930

072LO18  0

0800

6050

1150 .

3680 。

0710

2910

072

HACH 工NOHE1

00

7820

0820

4430

0850

9690

0760

6180 .

0940

4000

0720

3060

056 1

20

7740

0480

4510

0740 .

9300

0520

6360

0890

4170

0520

3280

047 1

40

7820

0660 .

4620

1070

9130

0590

6590

0830

4350

0710

3470

067

β O

8 0

6 0

4 0

2 0

0

40608010012014 。 16。 180 τ 図

12 βの全 解 析 結 果 及 び平 均 値 β

標 準 偏 差∂ 〔λ

LO) βx 0

8 0

6 0

4 0

2 0

0 コ ℃ λ OOO ◎

0

OTAFT  UD motion △ EL

CENTRO  UD motiop 口 HACHINOHE  UD motion 菖x+δx

0021τ+0

684

へも

      口 菖x

δx

0

0013τ+0

454

Bx=

0

0017τ+0

S69 40   60   80001 } a120   140 160 180 τ とで き る

し た がっ て

(22}

(23}式に おい て

対 象 と する地 震 波の スペ ク トル強さ

SI

h

, 

Tl1

, 

A

κ)を 用い れ ば

地 震 時 曲 げモ

メ ン トの最 大応 答 値が推 定でき よ う

 

6,

ま と め  本 論 文では 矩 形 平 面 を有す る屋 根 型 偏平球 形シエ ル の上 下 地 震 動に対する基 本 的な地 震 応 答 性 状 を分 析する た め

まず

上下動 定 常 white  noise による弾 性 応 答 解 析 を 行い

入力 位 相 差 が 応 答 値にお よ ぼ す影 響 を調べ

次に 実 上 下 地 震 動に よる弾 性 応 答 解 析 をパ ラメ トリッ ク に行い

直 応 力 推 定 係 数 a

ay と曲げモ

メ ン ト推 定 係 数 β

,角 の分 析

検 討に よ り, 上 卜地 震 動 を受け た場 合の 直応 力と曲 げモ

メ ン トの応 答 値 を近 似 的に推 定する方 法を提 案し た

以 下に得ら れ た結 果 を要 約すれ ば,   〔

i

) 上

F

動 定 常white  noise に よ る弾 性 応 答 解 析の       結 果

偏 平シエ ルの 4隅 支 持 点か ら入力さ れる外       乱を数 値 的に同

.・

と し

入 力 位 相 差 を考 慮し ない       場 合の応 答 値が

外 乱が独 立な場 合や 入力 位 相 差 y β 0

8 O

6 0

4 0

2 0

0 O λ=O

7 06 白 OTAFT  UD motion △ EL

〔:ENTRO UD motion 口 HACHINOHE  UD motion

By

ay=−

O

OO20τ・0

591   口 。

   

O

丶 \

言广 。

eOl5T・0

470

By−

By

=−

0

0010τ+0

35° 図

13 40  60  80  100  120  140 160 180 τ        b〕 β躍

βy の全解析結果 及 び平均値 βx

β。

標 準 偏 差 ar

iy (λ

0

7}    を考 慮し た場 合の応 答 値よりも大き く な る

ii

 ) 直 応 力 推 定 係 数 αx

 ay は

地 震 波

偏 平シエ    ル の重 量

ス パ ン長に関 係せず

ほ ぼ

定 値 を取

一 66 一

(10)

     る

 

qii

) 曲 げモ

メ ン ト推 定係 数

fl

角 は 地震 動, 偏    平シエル の 重 量

ス パ ン長に関係せ ず

T (

開    角半径 / シエ ル 厚 )に対 して

右 下り の

定 傾 向      を示す

 (

iV

) 用 い た 入 力 地 震 波は

EL −CENTRO

 

UD ,

     

TAFT

 

UD ,

 

HACHINOHE

 

UD

3

波と 限 ら れ    たもので はあるが, 偏平シエ ルの地震 時の直応 力    およ び曲 げモ

メン トの応 答 値 を略 算 的にめ る       方 法 (19

20

 22

 23式 )提 案 し た

  謝   辞  本 論 文 を ま とめ るに あた り

豊橋技術科 学大学助手

田 坂誠

博 士, 石 川 工業 高 等 専門 学校 助手

石 川 浩

郎 氏に貴 重な助 言 を頂き ま し た

又, 数 値解析の

部は 元 豊 橋 技 術 科 学 大 学 学 部 生

佐 藤 健 氏 (現東北 大学大 学 院 生 )の御 協 力 を得ま し た

ここに 厚 く 謝意を表 しま す。又,数 値 計 算に は豊 橋 技 術 科 学 大 学計算機セン タ

ー ・

MELCOM

 

M

 800

名 古 屋 大 学 大 型 計 算 機セ ン タ

ーFACOM

 

M

 320 を 使用 し たこ と を付記 し ま す

参 考 文 献 1> 真下 和彦

田中 彌 壽 雄

原 道 也 :片 持コ ノイ ドお よび    筒形シエ ル の動 的特性に関す る研 究 (その1)

日本 建 築    学 会 論 文 報 告集

第247号

昭和51年9月

pp

93

   100

2) 西 村 敏雄

新宮 清 志 :上下 動お よ び水 平 動 地 震を受け る    裾 梁 付 回 転 体シエ ル の動 的 応 答に関す る 研究

日本建築    学 会 論 文 報 告 集

第326号

昭 和58年 4月, pp

47

   59

3) 國 枝 治 郎 :球 形 シエ ル の上 下 地 震 動 応 答 解 析

日本 建 築    学会 大会 学 術 講 演 梗 概 集

昭 和55年9月

pp

1005

    1006

4) 國 枝 治 郎 :上下 動を受け る球 形シエ ル の非 線 形 応 答 解 析

    日本建築 学 会 大会 学 術 講 演 梗 概集

昭 和56年 9月

    pp

1157

1158

5) HaruQ Kunieda :SQlutions of Free Vibrations of Spher

   ica且Shells

日 本 建 築 学 会 論 文報告集

第325号

昭 和     58年3月

PP

57

66

6) 加藤 史郎

吉川 健二

横 尾 義 貫 :上 下 地 震 動 を 受け る球    殻の振 動特性

日本建築学会大会学 術 講 演 梗 概 集

昭 和    55年9月

pp

1003

1004

7) 長 浜 哲 史

加 藤 史 郎

横 尾 義 貫 :上 下 地 震 動 を 受け る球    殻の振 動特 性

日本建築 学 会 大 会 学 術 講 演 梗 概 集, 昭 和    56年9月

PP

1159

−−

1160

8)松 井 源 吾

瀬 谷 均 ;シエ ル構 造 物の位相差 を考慮し た    地 震応答につ い て

日 本 建築 学 会 論 文 報 告 集

第266号

   昭和53年4月

pp

73

85

9) 近 藤

平, 田中弥称 雄 :Parabolic Velaroidall Shellの

    上下 動 地 震 応 答に関 する考 察

日本 建 築 学 会 大 会 学 術 講    演 梗 概集 (構造 B)

昭 和61年8月

pp

265

266

10) 加 藤 史 郎

高 島 英 幸

西 薗 博 美

石 川浩

郎:位 相差の     ある上 下 地 震 動 を受 け る偏 平 殻の振 動 特 性  癧 動 方 程 式    の誘 導と解 析 法につ い て

日本 建 築 学 会 大 会 学 術 講 演 梗    概 集 (構 造B}

昭和61年8月

pp

263

264

11) 坪 井 善 勝 ;曲面 構 造

丸 善

昭 和40年

pp

145

187

12) 日本建築 学 会 :建 築 構 造 物の応 力 解 析

丸 善

昭 和49年

    pp

125

145

13) 田 坂 誠

一,

西 薗 博 美

加 藤 史 郎 :多重ホ ワイト外乱を受     け る構 造 物の平 均 応 答スペ ク トル 日本 建 築学会 大会学    術 講 演 梗 概 集 (構造 B>

昭 和60年10月, pp

265

    266

14)

S.

Tasaka

 

A.

 Muto

 S

 Kato;Nonstationary Stochas

   tic Response of a Beam to support  Accelerations

日本    建 築 学 会 論 文 報 告 集

第345号

昭和59年 年11月

    pp

30

40

(11)

SYNOPSIS

UDC:624.074.43:624.042.7:620.1

DYNAMIC

ANALYSIS

OF

ROOF

SHELLS

ON

A

SQUARE

PLAN

SUBJECTED

.

TQ

VERTICAL

EARTHQUAKE

MOTIONS

-Estimation

of non-dimensionalized maximum responses of

in-plane

and

bending

byHIDEYUKI TAKASHIMA, GraduateStudent of Toyohashi

Univ. of Technology, HIROMI NISHIZONO, Graduate

Studentof Toyohashi Univ. of Tenhnology and Dr.

SHIRO KATO, P[of.of ToyehashiUniv.of Technology,

Members of A.I.

J.

The effect of vertica'1 earthquake motions playsan

important

role

for

the

design

ofa

long-span

structure,

This

paper considered a simply supported roof shell on a squaTe plan and investigatedthe effects of vertical earth-quake motions considering phase

differences

of appLied earthquake motions at the'bases.

As a result of dynamic analysis, the

following

results

have

been

obtained,

i>

The uniform support excitations

from

earthquake, which are

identical

and have no phase differences,

give rise to

larger

values of responses than those

induced

from

some

different

classes of support tlons.

ii

)

By the analysis under three earthquake motions :

EL:CENTRO

UD, TAFT UD and HACHINOHE UD,

the presumptive coeffic'ient a

for

non-dimensional maximum

in

-plane

stresses remains constant with respect

to T

(=Rlt.),

and the presumptive coefficient

fi

for

non-dimensional maximum

bending

stresses varies

linearly

with respect to r.

In accordance with above-mentioned results, thispaper

has

proposed the presumptive equati6ns

for

non-dimensional maximum responses of in-planeand

bending

stresses,

'

参照

関連したドキュメント

(実被害,構造物最大応答)との検討に用いられている。一般に地震動の破壊力を示す指標として,入

ところで、モノ、ヒト、カネの境界を越え た自由な往来は、地球上の各地域の関係性に

絡み目を平面に射影し,線が交差しているところに上下 の情報をつけたものを絡み目の 図式 という..

 原子炉建屋(R/B)及びタービン建屋(T/B)の汚染状況は、これら

9 時の館野の状態曲線によると、地上と 1000 mとの温度差は約 3 ℃で、下層大気の状態は安 定であった。上層風は、地上は西寄り、 700 m から 1000 m付近までは南東の風が

地下水採取等対象物 質と地下水採取を行う

(1) 建屋海側に位置するサブドレンのポンプ停止バックアップ位置(LL 値)は,建屋滞留 水水位の管理上限目標値 T.P.2,064mm ※1

認知症の周辺症状の状況に合わせた臨機応変な活動や個々のご利用者の「でき ること」