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

荷重の確率過程の組み合わせにおける超過確率の理論解とその応用に関する考察

N/A
N/A
Protected

Academic year: 2021

シェア "荷重の確率過程の組み合わせにおける超過確率の理論解とその応用に関する考察"

Copied!
11
0
0

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

全文

(1)

【論  文

1

UDC :624

042

2 :519

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

1989 年 lt月

確率

過 程

超 過 確 率

      

論解

と そ

応 用

考察

正 会 員 正 会 員 正 会 員

* *

** *  

1.

ま えがき  建 築 物に作 用する積 載

風お よ び地 震荷重などは

そ れ らの生 起, 作用 時 間, 強さ に不 規 則 性

時 間動 を もつ こと か ら確率過 程 (stochastic  Plocess )と して モ デル化さ れ る。 マ クロ的 な 時 間変 動過程とし て特 徴づけ ら れる荷 重の確 率 過 程モデル につ いては 文 献 1)に概 述さ れて いる

これ らの荷重の確 率過程の組み合わ せの レベル横 断 率, お よ び ある時 間 区 間にお け る 超 過確 率あ るい は最大値の計 的 性 質 を求 める問 題は

時 間 依 存 信 頼 度 解 析 (time

dependent

 reliability  analysis )や統 計

確率 論的手法に よ る限界状態 設 計 法に おける荷 重 係 数の 設定に関 連す る問 題と し て理 論お よび応 用 的 研 究が重ね ら れて い る21

n} 。

Turkstra

に よ る

Turkstra

s rule 法2 ’

Wen に よ る 同 時 生 起 (率 ) 法 (

load

 coincidence method )4】

5)

  Larrabee

 Cornellに よ る点横断 (率 )法

(point

crossing  method i 〕な ど が こ の問題の代表 的な

解析 手 法と し てあ げら れ る

 

TurkstTa

s rule 法は

2個 あるい は それ以 ヒの時 間 変 動 過 程の組み合わ せ に お ける最大 値 を推 定す る経 験 則 的 な手 法であ り

実用

ヒの簡便さ を具え てい る が, 理 論 的 根 拠と 適用 性の不 明確さに不 満足 さ が あ る

。Wen

の手 法 も 多数の確 率 過 程の み合わ せに お け る 同時 生 起 を経 験 的に組み入れ た手 法で あ り

比 較 的に粗 〔sparse な 時 間 変 動 過 程の組み合わ せ に お け る超 過確 率の推 定に 適 用され る

 

Larrabee.

 

Comell

t

再 生 矩 形ル ス過 程の組み合 わ せ につ い て期 待 横 断 率の 理 論 解 を 導き点 横 断 公 式

(point

crossing  

formula

)を 示 し て い る

こ の公 式は

ほ かの確 率過程のみ合わ せに も拡 張 適 用され るが

点 横 断 率とレ ベ ル横 断を ボ アソ ン事 象と する仮 定に よる超 過確 率の定法は超過確 率の

つの上界値を 与え る

点  ’ 名 古屋 大 学   助 手

工 修 牌 名古屋 大 学 教 授

工 博

林零

株 式 会 社 竹 中 工 務 店

工 修   〔1999年5月8H 原 槁 受理

1989年8月Zl日採 用 決 定 } 横 断公式は 理論 的に明 快で あるが

こ の手 法では多 数の 時間変動 過 程の組み合わ せにおい て畳 込み多 重 積 分を必 要と し

応用 上の雑さ を伴う

 統 計

確 率 論 的 考え方に基づ く限 界状 態 設 計 法の構 成 で は

荷 重の組み合わ せ につ い て の合 理 的な取 扱いが重 要と なる

本 稿で は

マ クロ 的な時間変動 を伴う荷重の 確率過程のみ合わ せ に お け る超過確 率の解 析 手 法と理 論解につ い て述べ

本稿 2

方 形パ ル ス いはイン パ ル スに より構 成さ れ る確 率 過 程の組み合わ せ にお け る状 態 確 率を適 切に組み人れた超 過 確 率の理 論解 を示し て い る

こ れ ら の理 論 解は打ち 切 り定 常過 程に っ い て の近 似 理 論 解であるが 表 現の簡潔さ と応用 ヒの 柔 軟 性 を具えている

3節で は

理論解の限 界態設 計 法にお け る荷 重 係数の 設定問 題へ の 応用につ い て要 約し て述べ

4数 値 計 算 例て い る

本 槁内 容 は部分的に文 献

10

,11

>に報 告さ れて いる が

本 稿で は これら の報 告 資 料を整理し て ま とめ たもの で ある

な お

本 稿に示した超 過 確 率の理 論 解と そ れ らの限界状態 設計 法へ の応 用につ い て のよ り詳 細な考 察は

文 献1 )

9 ) の 期待 横断 率 法と

Turkstra

s ruie に基づ い た実 用 的 簡便さを配 慮し た荷重係数の設定手法との対 比 を含めて 続槁で述べ ら れ る

 

2.

超 過 確 率の理 論 解  

2.

1 確 率 過 程と その組み合わせ  本 稿で は

荷 重の確 率 過 程モ デ ルを対 象と し て

非 負 の値 域 を もつ 定 常 確 率 過 程 を 取 扱

ま た

確率過 程 おけるパル ス あるいは インパル ス の生 起

継 続 (作 用 ) 時間

強 さは統計的に 互いに独 立 と し

かつ パ ル ス は 零 再 帰 性 をもつ も の と す る

マ クロ な時 間 変 動 を 伴う 確 率過程モ デル b’ と して図

1に示す ボア ソ ン方 形 波 過

程 (

Poisson

 square  wave  process

以下で は

 PSW 程と 略記す る)

矩 形パ ル ス過 程 〔rectanguLar  pulse

process

 

RP

程 )

続矩 形パ ル ス過 程 intermittent reCtangular  process

 

IRP

過 程)

お よ び ボア ソ ン イン

(2)

w t PSW  process   A

μ A RP process T

 

TT λ lRP proo巳ss    A

fi

t w t 7 τ丁 λ PI process   A t w {t   λ 図

1 マ クロ的 時 間 変 動 過 程 t t り あげ, これ らの確率 過 程の スカ ラ

和の組み合わ せに お け る超 過確率の理 論 解につ いて述べ

本稿で の記 述 で は

各 確 率 過 程に お けるパ ルス あ るい はイン パ ル ス の 単 位 時 間 当た りの生 起回数の期待を λ

パ ル ス の継続 (作用 )時間の 期待値をμ とし て

前 者をパ ル ス ある いはイン パ ル ス の生 起 率, 後 者 を継 続 時 間と略 記する

な お

,RP

過 程お よび

IRP

過 程は

定 時 間 間隔τ ご と に規 則 的に生 起する パ ル ス によ り 構 成 さ れ る

こ れ ら の過 程では パ ル ス の生 起率お よ び継 続 時 間は確 定 変量 と し

λ お よ び

P

の よ う に 符 号を付して期 待 値の λ

お よ び μ と区 別す る。 また

各 確 率 過 程におけ る パ ル ス あ るいはインパ ル ス の さを 表 す 確 率 変 数 を W , そ の 累積 分 布 関 数をF試ω)と す る

 2

2PSW 過 程と

PI

過 程の組み合わ せ   PSW 過 程 を

W

,(

t

),  PI 過 程 を

W

,( 

t

)と し て, 両過程 の ス カラ

和の確 率 過 程を

W

(t)と す る

 

PSW

過程 W,(

t

)の

つの パ ル ス の継 続 時 間を表す確 率 変 数を

Tl

と す る と

 

PI

W

,(t)の イン パ ル ス が T,時 間 区 間に h回生 起す る確 率

qn

h

;T,)は 次 式で表 され る

 

 

 

・1掘

・xp

 

 

 

 

 

 

 

 

exp

・t)d・

、.

 λtXtll 十 λ2μ1

t                    ;

h

= 0

1 ,2,

… …・

…・

1  以 下で は

q、、〔

k

T

,)を確率過程 酬 孟)と W,(t>の組

32

み合わ せにおける状 態 確 率 とい う

 また

PI 過 程の イン パ ル ス の T,時 間 区 間 最 大 値 (極 値 〉を表す確 率 変 数 を w,

T

と し,  T,時 間 区間に イン パ ル ス が生 起 し な い と きに は

WaT、

を 零 値 (

fictitiOUS

impulse

)と して取 り扱う と

の累 積 分 布 関 数 は 次 式でさ れ る。

       の

     

Fw 、

rl(w )

q、2(0;

Tl

)ひ(ω)+Σ qn〔

h

;Tl)置照ω避       k

1

 

 

 

 

 

 

 

 

1

 

、恥 )

       

1

+ λ,μ1 喇 + 、+ λ,μ1醐 ω)                

 

9・

…99・

 (2 ) ま た

耽7

の確率密度 関 数は

2

)式の導 関 数と して 下 式となる

 

 

 

f

・・

{・)

q・・(・…

T

,)・麟

1

li

ω ・瓦 ω      

 

t・

 

tt・

 

tt・

 (3) こ こ に

U

(w)お よ び δ(w)は 

HeaviSlde

の ステッ プ関 数お よび

Dirac

の デル タ関数

 

F

賑ω)

1

(w)で あ る。  

PSW

過 程は ボア ソ ン過 程で あ る か ら

W (t>が (0

 T 時 間 区 間に レ ベ ル r を少な く とも1度 超 過 する確 率 (以 下で は

超 過 確 率と略 記す る)∫P(r ;T)は

下 記の よ うに得られ る

 

 

 

P

(・ ・η=

1一

γ・譎 ・臥 . 。

(・)

il

       

1

expi

λユT(1

Fw,+w,

T1(r))

1

     

 

一・

(4) こ cに

,Fw,

+隔 (ω)は確 率 変 数 W,と鵬

r

の和の累 積 分 布 関数で あり, (2)式を用い る と次 式で表され る。

   

F… 一 ・

h(zv)

x

1

+画

(。

副 鉚 砺 )

d

・ ,      

tt・

tt・

 

 

 (5> (5>式 を (4) 式に用い ると確 率 過 程 W (

t

)の超 過確 率の理 論 解は

次 式と なる

  

 

・P(・ ・T)

exp

・,・

1

f

,  r

  

    

  

1+ λ3μ,

(r

ndlJ

 

fWi

〔ω1)

dw

     

…・

……・

……・

……・

…・

(6 }   上 式の 超 過 確 率の理 論 解に基づい て レベ ル超過 率 〔exceedance  rate , failure rate >

hXr

)が, 次 式で定義 され る

1

..

一 ..

..

..

..

、,、

・,

1

nd

,。 ・ ;

,(。

Wl 細 〉伽

     

   

9−・

・・

  (8 )

(3)

こ こに

ん∫(r)

dt

確率 過 程 障(t>が (0

 

t

) 時間に レ ベ ル r を超 過 し ない条 件の 下で

t,t

dt

時 間に レ ベ ル r を超 過す る確 率 を表す13)

 (6)式の超 過 確 率の理 論解は

文 献 12 >に異な る解 析 手 順に よっ てか れ て い る ま た

(8) 式の レベ ル 超 過 率は, 次 式の よ う に表さ れ る。

 

 

 

h

・(・)一 r

f

, ’

1

      

fWi

(Wl

d

ω,

 

tt−・

 

一・

 

tt・

… 

(9 )

     

窟 (λ1

λ2)F 魯1(r)十 λ2F 移、

Wt(r)                ;λ2μ1《1

0

の とき

一 ……・

(10 ) (10 )式お よ び (9) 式の 超 過 率 の 近似理 論解 は

Larrabee,

 

Cornell

に よっ て 導か れ て い る点 横 断 率

(point

crossing  rate vr), お よ び reduced  point

crossing  rate VR(r〕7) )の理論解と同じ で ある

(10式 の近 似 を用いる と 超 過確 率は

次 式とな る

   

rP (r ;

T

)=rl

expl

h

ノ(r)

Ti

            

1

expl

(〔λ,

λ2)

F

(r)              十 λ2F 昏、

w,(r)}

Ti・

 

tt・

 

t…

 (ll) 〔6 )式と (11)式の超 過 確 率の理 論 解の精 度 と適 用 にっ い て は 文 献 12)に モ ンテカル ロ シ ミュ レ

ショ ン に よる計 算 結 果 との比 較 が 示さ れ てい る。 ま た

こ の

PSW

過 程と

PI

過 程の組み合わ せにお ける超 過 確 率の 厳密な 理論 解の 重積 分 表 現は

文 献8)に述べ ら れ て いる。  2

3 PSW 過 程と RP 過 程の組み 合わ せ   PSW 過 程

W

,(

t

)とRP 過 程

W

,(

t

)の組み合わ せに お け る状態 確 率は

近 似 的に次 式で表 さ れ る (

Appendix

参照)

     

O

 

 

 

q

’2(

h

 ‘ T ’)”

1

ha

、Xtl

 λ2μll 十 λ2ttl

k

1 な お

RP

過程で は β,

=1

/λ2 である。 ;k

0, ;h

1,2,

12

 

12

)式を用い る とRP 過 程の パ ル ス の

T、

時 間 区 間 最大値を表す確 率 変 数 鵬

,,の累 積分布閧 数は

近 似 的 に ド式で表 され る

   

凡、

(ω)駕 Σ σ、2(ん;

71

1

凡 (ω)ド       配

il

 

 

 

 

 

 

 

F

(・v>

1−

1

1 恥

ガ}

                  1十 λ2Pt,FC} ,(w)                   

 (13 )        

Fw、

{zv}

      

。)

”… 『

”齟

4 ただ し,

N

7

ソβ

(整 数〉と する

また, (14 )式は

N

が大きい と きの近似で あ る。  確 率 過 程

W

t

)(

 

W

t

既 (tの超 過 確 率の理 論 解 は

下 記の よ うに得られる

   

P

(r ;

T

)・・

1−

expF λ

T

1− F

.w

(r))}      

一・

 (ユ5)

       

・ ・

exp

・,・T

1

 

 

 

 

 

 

 

 

1

£

1

砺 )

fWL

( w

)・・,

     

……・

…t・

…t…t・

(16> (16)式の理 論 解は

(14)式の近 似 を用い て い る。 こ の とき の 恥 1

w、

T

(ω〉は

次 式であ る

 

 

 

Fw]

・w

(・)・

f

,  W

 

1

 

筋 )縞 )・筋      

 

一・

 (17)  こ の確率過 程み合わ せの と き に は

超過 率 源 の は下式と な る

   

h∫(r)

λ111

− F

+w2

T

(r)}

 

 

 

 

 

A

1

1+

f

(・v,)

dw

     

一 ……一 一 ・

…・

………

(18)  2

4PSW 過 程 とIRP 過 程の組み合 わせ  

PSW

過 程 Wi(

t

)と

IRP

過 程 W,(

t

)の組み合 わせ にお ける状 態 確 率は

,Appendix

を参照 して次 式の よ うに表 さ れ る

σ1、(

h

T

、)=

 

 

 

ユ+孀

       

1・ 角

气菰

 λ2μll 十 λ2μ1

’t       ;

h =

2

3

             

 (19 >

i

・μi

μ ’

1

4i

       

λ,μi

堕 +ρ2

 

2+λ2μ・ こ こ に

IRP

過程の パ ルス の継 続 時 間を表し

τ

島 と す る。  この確率過程の み合わせ に おける超 過の理 論 解 は

(19> 式の状 態確 率を用い て導くこと がで き る が

解 表 現が複 雑と な る

μ1>角 の と き に は, 2

3

に述べ た

PSW

過程と

RP

過 程の組み合わ せにお け る超 過 確 率の 理 論 解 を近似 的に適 用す ること がで き る

 2

5RP 過 程と

PI

過 程の組み合わ せ

 RP

過程

Wi

(t}とPI 過 程

W2

{ 

t

み合わ せに お け る状 態確率は

次 式で表さ れ る

 

 

 

・1師 1)

PL

 

ft1

A

(・・)

 

また

PI 過 程の インパ ル スの β,時 間 区 間 最 大 値を表 す確 率 変 数 を W,

n

と す る と

そ の累 積 分 布関 数 は 次 式

(4)

と なる。

   

F・ 、

fi

(w)

e

A

U

ω)+e

λ 2 ’

PF

Cus

e

λ

11

     

…・

………・

………一 ・

(21)  

RP

過程の パル ス が (

Oi

 

T

)時

N

=T

i

>1;整 数 〉 回 生 起 す る と す る と

こ の わせ 過 程 W (

t

) (

wrt

)+

W

,〈t))に お け る 超 過確 率と 超過 率の理 論解 は

次 式の よ うに得られる。

   

ノP(r ;T)

1

IFWi

Ψ2

a

(のド

………・

………

〔22)

       

= 1

exp 亅

λ , 

T

1− F

+w

fi

(r))}        

 

tt・

 

一・

 

(23)

   

h∫(r〕

λlll

FWI

w2

a1(r)

1

 

(24) こ こに

   

F

(Ul)

exp  

l

・鱗 (w

WI )}      

f

剛i(w,)

dw

 

一・

 (25)  2

61RP 過 程と PI過 程のみ合わせ  IRP 過 程 蝋 ‘)と PI 過 程

1

愀‘)の組み合わせ にお け る状 態確 率は

次 式とな る

 

 

 

・1軅

μ

 

…・

………一 ・

…・

(26)

 

 

 

q。2・鱒

e

一 ・

…・

一 …

27) こ こ に

q。t(

k

;ρわ は

時間 区 間 鮮 (

τ

鳥 )に PI 過 程の イン パ ル スが

h

回 生 起す る確 率 を表す

 PI 過 程の イン パ ル ス の

ft

、お よび 貯時間 区間 最 大 値 を表す確率 変数を そ れ ぞ れ 耽魚

お よび

Wl,

fi

と す る と

両 者の 累 積 分 布 関 数は次式で表さ れ る

   F

曜 、

fiL(ω)

e

λ2μ

σ(ω)十e

λ2

μ

殉閣

A’

                 

 

r・

 

r・

 

r・

7・

 

r…

 

7…

 〔28 >

   F

畷 ω

e

A

P「

u

z ・)+e

λ2 μ臨 佃Le

λ2 μ’                

一 一 ・

……・

…・

………・

29

)   確 率 過 程 W (t) (

W,(t)十隅(t))の 超 過確 率と 超 過 率の理論 解は

(28)

(29>式 を 用いて次式 と な る

   

P

(ア ;

T

1

IFw

+曜

a

(7}鬥

F

肱 卸〔r)} 削

……

(30 )            富 1

exp  

l

λ1 T [λ2(1

ft

,/τ)zrF 昏2(r)

        

十1

− F

膨 1 +w

fi

(7)コト

…・

− r…・

(31)

   h1

(r)= λ ,

1

λ,(

1−

Pi

/r>τ

F

魯2(r>十1

− Fw、

+ 剛

Db(r)}        

 

r・

 (32) こ こ に

   

・・

・)

・ ・p 卜為鳳 (・

wDl              

fWi

〔ω 1)dw1

………・

…………

(33>  

2.

7PSW 過 程

 

RP

過 程お よ び

PI

過 程の組み合わ       せ  

PSW

過 程 を

W

,(

t

), 

RP

過 程 を Wz(t)

  PI 過 程 を

W

,(

t

)と して, これ らの 3つの確 率 過 程の スカラ

確 率 過 程 を W (t)とす る

 こ の組み合わ せ で は

の状 態 確 率を用い る。

一一

34

qn (κ、飼

 

、h

。,1,,,

_

(34) q、z(

k

;T,)=

h

0

1十 λ2μ1

 λ2μ11 十 λ2μ1

s” ;

h =0,

k =

1

2,

35 } こ こ に q23(

k

ft2

は RP 過 程 隅(t)の 1つ の パ ル ス継 続 時 間 区 間 ρ2 に

PI

過 程 W,(

t

)の インパル スが

h

回 生 起する確 率

同 様に q12

h

T,

〉は

PSW

過程 隅(

t

}の 1 つ のパ ル ス継続時間区間

T

,に

RP

過程

W2

t

)の パ ル ス が

k

回 生 起す る確率を表す

 超過 確 率の 理論 解を導く 上で

RP 過 程とPI 過 程の ス カラ

和の確 率 過 程 を

W

,,(

t

)(

 

W

,(

t

)十

W

,(

t

)) とし て

W,,(t)の T,時 間 区 間 最 大 値 を 表す確 率 変 数を 隅、

, 、とする

確率 変数 嬲,

7

の 累積分布関 数

Fw ,

話ω) は

14

)式 と同様に次 式の よ うに表さ れ る

        

n

   Fw,

馬 .ω = Σ q、2(h;T,)

iFw

rv,

ll

{w)

IR

      s

1                

Fw、

 

 一

a

(ω〉

      

= ・+

a

,。、F

”『

… ”’

36  上 式 中の W31fi

 

PI

過程の インパ ル スの β2 時 間区 間最大値を表す確 率変数で あ る。 ま た

,F

.ve、

o

(ω)は 確 率 変数

W

, と 既

a,の和の累 積 分 布 関 数であり, 次 式 で表さ れ る

   

F・

・w、

il

(・V)

ω expi

為鱈 (・

w

)}              

fw

,(ω2)

d

ω2

 

7P− ・

rP・

 

 37 上記の諸 式を用い る と

確 率過程

W

(t)の超過 確 率お よ び超 過 率の理 論 解は次 式の よ うに得られ る。

 

 

P

(r ・

T

)・ 1

・1・

1

・。 、

w。

T

(r)

il

       

= ユ

expl

λ ,T (1

FWt+w

1(r))}        

 

一・

 

rP・

 

(38)

   hr

(r)

λ,

1

F障 L

w

r

(r)

1

tt・

 

t・

 

tt・

t・

(39) こ こ に

F… v・・

・・

1

fW

、+

1

。 、

f

・ ・(Wi ・

d

・,      

9・

一・

 

一・

9・

 く40)  その過 程の組み合わ せ に おける超過確 率の理 論 解につ い て も上 述と同 様の解 析 手 順によっ て導くこと ができる

こ こ で示した理 論 解の応 用 上の柔 軟 性 と精 度 に っ いて は次節以下に述べ る が そ れ ら の

文 献

O

,ll

,12

}に報 告さ れている

 

3.

超 過 確 率の理 論 解の応 用  前節に示 し た確 率過 程の組み合わ せ に おける超過確 率 の理論解は近 似解であるが, 応 用 上の柔軟性を 具え てい る

こ こ で は

荷 重

耐 力 係 数 方 式の限 界 状 態 設 計 法に お ける設 計規 範の構 成と荷 重

耐 力 係 数の決 定 問 題へ の

(5)

応 用につ いて述べ

 建築 物に作 用す る積 載

風 お よ び 地震 荷 重な どは 生 起

作 用 時 間

強 さに時 間 変 動

不 規 則 性 を もつ こと から

前 節で述べ たマ クロ 的な時 間 変 動 過 程 (確 率 過 程 ) とし て モデル化さ れ る 文 献 1), 9)参 照

 限 界 状 態 設 計 法で は, 構 造 物 あるい は構 造要素が設計 に おいて意 図さ れ た機 能 また は条件に適さ な く な る状態 を限 界状態と して 築 物で は安 全 性と使 用 性に っ い ての限 界態を考え るの が

搬 的で あ る。 これ らの 2つ の限界状態に は終局耐力 限 界状態と使 用 限 界 状 態が 設 定 され

設 計 条 件は 「構 造 物 または構 造 要 素 が, ある 定 め られた時 間 区 間 (以 下で は 基 準 期 間 とい う)中に 少な く と も

度限 界状 態に達する で あ ろ う確率 (以 下で は, 限 界状態確率とい う)を設計に おい て指定さ れ た確 率に等しいか

ま た は小さ く す る」設計 規範によっ て構 成 さ れ る

 基 準 期 間 を丁年

構 造 要 素の 限 界 耐 力 を

R ,

荷 重 効 果の組み 合わ せ を Σ Cx Wx

t

と する と

確率 論 的考え 方に莚つ く設 計 規 範は次式の よ うにさ れ る。 ・・

b

R

・・W・(t)… T

・ ・(

β}      

 

t・

 (41)      Ps

1

P/≧ φ(β)

………・

…・

…・

 

………・

42

) こ こ に

Prob

E

は事 象

E

の生 起確 率

φ(

)は標 準 正 規分布関 数

β は信 頼 性 指標を表す。 また

(41) 式 で は

R は確 率 変 数,  Wx (

t

)は荷重の確 率過程であ る。

P

∫および

Ps

は, 

T

年 期 間に お け る限 界状態確 率お よ び信 頼 度 を表 し

Φ(

β)お よ び φ(β)は設計におい て目 標と す る限 界 状 態 確 率および信 頼 度 を表 す。  上 記の設 計 規 範に基づい て荷 重

耐 力係 数 方 式の設 計 条件式が次式の よ うに導か れ

同 式 中の荷 重 係 数

7x

は 荷 重の組み合わせ と 設 計におい て目標と す る信 頼 度レベ ル に基づい て決 定さ れ る

     φ

R

。≧ΣCズ み

臥 。

…………・

………・

(43 )       濁 こ こに

添字

X

,D ,

 

L ,

 

S ,

 

W ,

お よび

E

と し

以 下のではこれ らの添字 記号 積 載 風, お よ び 地震 荷 重に関す る諸量 を表す もの と す る

  (43 )式 左 辺の φ お よ び Rn は 「耐力 係 数」お よび 「限 界 耐 力の公 称 値 」 を 表し

右 辺の γx は 「荷 重 係 数 」を

W.

n は 「荷 重の公 称 値 」 を 表 す。 ま た

Σ]は

「荷 重の 組み合わ せ」を表 し

Cx は係 数 倍 荷重 γx Wx

n と そ れに よっ て構 造 物に生 じ る力 (応 力 )と を関 連づ け る係 数で ある

  本 稿で は

2節で述べた荷 重の確 率 過 程モ デル とそれ らの組み合わ せにおけ る超 過 確 率の理 論 解に基づ く荷 重 係 数の算 定 式の誘 導と

そ の計 算ア ル ゴリ ズム につ い て 要 約して述べ る

以 下に は

荷 重の組み合わ せ の 1例に つ い て述べ

建 築 物

般 的な設計に想 定さ れ る荷 重の 組み合わせ と各 荷 重の具 体 的な統 計 資 料に よ る荷重 係 数 値の算定例につ いての考察は続稿で述べ る

 3

1 固定荷重と積載 荷 重の組み合わ せ  固 定 荷 重 W。と積 載 荷 重 肌 の組み合わ せ で は, (41) 式は次 式の よ うに記 述さ れ る

 

  

.Pω

一P

・ ・

b

叺(

t

… T

        ≦φ(

β)

…・

…・

…・

…・

……・

 

…・

…・

(44)     ΣコWx(

t

WLs(t>十WLe(t)= 既 (t

 (45) こ こ で は

固 定 荷 重 鴎 は時間変動が相対的に小さい も の と み な し て確 率 変 数とし て取 扱い

積 載 荷確率過 程 蹴(t)に つ い て は

持 続 的な積 載 荷重

WLS

(t)と

時 的な積 載荷重 蹴。  の確 率 過 程の ス カ ラ

確 率 程モ デル を考え る

文献

9

)参照

  (44 )式 は

下 式の よ うに表さ れ る

    

Prob

Z

− Z2

0

]≦Φ(

β)

……・

……・

……

(46 ) ま た は

    Prob [Zi

Z2≧0]≧ φ(β)

…・

…・

……・

……

(47> こ こ に      

Z

=R − W

…・

………・

…・

……・

………・

……

(48)      Z,

max Σ Wx(

t

……・

………・

……・

(49)       『     x 上 式におい て

,Z

、は時 間に依存し ない 確率 変数

 

Zi

は 時 間変動荷 重のみ 合 わ せの基準 期 間 最大値を 表 す 確率 変数であ り

以 下に両 者を分 離して取り扱 う信 頼 度 解 析 手法を述べ る

  (

46

)ま た は (47}式の信 頼 度 設 計 問 題における限 界 状 態 関 数を G(Zi

 Z)と す る と限 界 状 態 方 程 式は次 式と な る

    

G

ZI,

Z,)

Zl

Zt

o

 

r

(50)  信 頼度設計問題に お け る解析 手 法とし て は幾つ か の実 用 的な手法が あ る が, こ こ で は確 率 変 数の標 準 正 規 変 数 へ の変 換お よび標準正 規 変 数 空 間 (σ

空 間 )にお け る 限 界 状 態 方 程 式と線形化と線形 化点に基づ く設計 条 件 式 の構 成 手 法 を適 用す る

文 献

9

)参 照。  確 率 変 数 Zτお よび Z2 の累 積 分布関 数を FZi(zD およ び

Fz、

(zと す る と 標 準正規 変 数へ の変 換お よ び U

空 間に お け る線形 化 点 座 標 (u密

磁 )は下 記の よ うに 表され る

     Fz

(Zl)

φ(u1)← うZl

=Fi ,

LΦUl

r・

7…

 

P・

r・

 (

51

   

Fz,(z2>

Φ(u2) ∈→ x :

ニF

轟 1 (φ(u2))

 

(52 )

   

u苧

az

β

…・

………t・

……・

…………・

……

53)     u穿

az

β

……・

…・

…・

…………・

………

54

) こ こ に

α z

(ノ

2)は U

空間にお け る 座標 原点に 対す る方 向 余 弦であり

α、L〈0

α z、>0とする。   〔51)

(52) 式か ら線 形 化 点を設 計 点と す る設 計 条 件 式は

次 式で表さ れ る

(6)

     

Zf − Z

夛≧

0 ・

 

 (55 ) こ こ に,      2

1

Fi;(φ(az

fi

))

…・

…・

……・

…・

……・

…・

(56)      Z夛

Fを21(φ(α z2β))

……・

…………・

…・

……・

(57 ) ま た

(46)式 より      R* ≧ 回!十Zf

 

9・

 

一・

 

tt・

 

58

  (56 ),(

57

)式の逆 変 換 表現 は, そ れ ぞ れ (59)

(60) 式の よ うに記述さ れ る

   Prob

R −

w

.≦

z

裘]

Φ(aZlfl)

………・

{59)

  

 

b

9

剛 ・Zf

・(・・,

B

……・

(・・)   (

59

)式お よ び (

60

) 式か ら荷 重 係 数と耐 力 係 数の算 定 式を導くに は

(59)式につ い て は静 的 信頼度解析 手 法 を

(60)式につ い て は時間変動 荷重の組み 合わ せ に 対 する信 頼 度 解 析 手 法 を適 用す ればよい

 

(59)式に前 述と同様の正規 変数 変 換 と 線 形 化手法を 再び適 用すると, 線 形 化 点, 設 計 点は下 式のように表さ れる

   

u育

=一

αRαz

β

一・

 

tt…

 

r・

9・

 

 一・

 一

(61)

   

u

=−

aoaz

β

………・

………

(62)

   

R

FEL(φ(

aRα z

β))

………一

63 )

   

叫zま

=F

1

(φ(

αDaz

β))

………・

…………・

 

64

 

確 率 変数

R

お よ び 既 が対 数正規変数の と きに は (63}

(64)式は次 式と なる

 

 

 

・*

expl

a・a・

β… Rln

(・5 )

 

 

 

W

・一

exp トa・a・・

BaL

…・

… 6・ こ こ に

X

 Vxお よび σx は

確 率 変 数の X の平 均 値 変 動 係 数お よび標 準偏差を表す。

 

(65>

(66)式か ら

限 界 耐 力の公 称 値

Rn

と 固定 荷重 の公称 値

W 。

n に乗じ る耐 力係 数 φと荷重係数 γDは

次 式と な る。

 

 

 

exp ・娩 酬

………

(…           1     

万。

   

7

’D

7

expl

α・α漁 ・・

門’

(68) こ こ に

α R お よ び aD は

    αR

=−

a,。 ,/ σ

1

。、+晶 。ゲ

…・

………

69

)     ap

σ、n . 。/ σ m ,+σ]

w

 

…・

…・

…・

……t−tt

(70) 前 述の 60)式

同 式 中の時 間 変 動 荷 重の組みわ せ に (45 )式を代 入 する と

 

  

b

9

1

・W ・・(・t)・W・

・(・)

1

・zr

・崛      

・・

 

一・

 (71) … b

9

・ 脚 ・職 )

1

・・ }

     

9・

 

一一・

 (72) (72) 式 中の

w

,(t)お よ び vv,(t)は

そ れ ぞれ w、

( t)

36

および 蹴。(

t

)を表し,

2

節で述べ た確 率 過程の組み合 わ せ に おけ る付 加 番 号 表 示と し た もの であ る

 (72}式は 2節で述べ

PSW

過 程 とPI 過 程み 合わ せ に お け る超過確 率の 理論 解の (6)式を用い る と 次 式の よ うにさ れる

  

 

1−

T

1

隔 。

1

exp  

l

T

1

  

  

1+

(。穿

副 ム (w,・

d

・,

・(a・

β)      

………・

− tt− tt− ・

…・

…・

(73)  な お

こ こ では

持 続 的 積 載 荷 重 嘱(t)(

既湿

D

に は確 率1で生 起する条 件 を付し て (6)式の超 過 確 率 の理 論 解 を修 正し た次式 を用いている

P

〔・・

r

÷

1−

exp

・,・

T

1−

1+

(剛 畑

一 ・

4 ) (73 )式は 次 式の よ うにき換え られ る

 

 

 

+厨

z

Wl)

f

(・,,)

d

・,

Φ(

R

) こ こ に

 

7・

−7r・

 

7・

 

 

P・

 (75

β

一 一

i

Yt

 

l

1

(・

・ う婦 ・・

     

…一 ……一 …・

……・

……

(76) さ ら に

(75)式は

等 価 的に下記の よ うに取り扱うこ とができる

    

Prob

W

,十

W2

Z

φ(β)

…・

…………・

t

(77 ) こ こ に

W ,は

次 式で表され る累 積分布 関数 を もつ 確 率 変 数である

 

  

F・

(・V)

、+。。

i

。免(。)

一 ・

……・

一 一

8   (77 )式に先に述べ た 正規 変 数 変 換 と線 形 化 手 法を適 用すると

U

空 問に おける線 形 化 点は

    u許

a 、β

…・

…・

……・

…・

…………・

…・

………

(79)     u才; aiβ

………・

…『

…・

………・

……・

…・

…・

80 こ こ に

   

al 十α茎

=1,

 

 

α 1

α:>

0・

 

99…

 

一・

一・

 

81

) と な る

ま た

設 計 点で は次 式が成り立っ

 

 

 

F・ ,・・

rl

1n W

ln Wi     σln 曜 1

・・u・・

……・

(…

   

F盈,(vik;

1

 

      

1

      = φ(u才)

 

 

 

 

 

 

 

1+ λ・μ1φ

1nW 才

亘n 鵬     σIn Wt

     

…………・

……・

………

〔83)

t

だ し

82

83

)式で は

確率 変数

W

,お よ び 隅 は 対数

IE

規 分布に従う ものと してい る。  (82>

(83)式か らWf お よび 孵 は

次 式で表 さ れる。

(7)

 

 

 

wr

− −

i

exp

1

死 一・

一 …・

(・・)

 

 

 

exp

Φ

φ(

α 2β) λ2μ,Φ(α調 )

al

  ・ t

        

W2

 tt

 (85 ) したがっ て 持 続 的 積 載 荷 重 既お よび

時 的積載荷 重 WLeの公称 値に乗 じる荷 重 係 数は

下式と な る

 

 

 

7Ls一

 

π

exp {

……

…86・

 

 

 

XLe− 、

臨 exp

1

 Φ(

αL εβ) λL9μz3φ(α乙eβ)

 

 

 

 

 

1

ー ……・

…………・

87

・ ま た, 積載 荷重 呪 につ い て は,      

W

π

VLs.

n十

WLe.

n

7・

 (88)       WL

n        既

。        

 (

89

   

「・

 7L

+ γ・

匿 .

 上 述の諸 式か ら 固定 荷 重と積載 荷重のみ合わ せに お け る設計 条件式と耐 力 係 数お よ び荷重係数の算 定 式が

下 記の よ うにま とめ られ る。      φRπ≧ γbしVqn十 冫乞yレZL

n

−・

 

P−P9・

 

 (

90

 

 

 

exp ・

一 鰯

…一 ・

91

・ 7・D

exp ト 噛 ・恥

・・

・… 7・・

exp 襯 嚇

1

 

..

1

  

佩 ・ ・

p

φ

1

λLettLsdiφ(

aLe(α Leβ)β)

al・・w

e

1

……・

…・

…・

………一

(・・)

β一一

1

1

。呻 昭

λ 1 ・ 〉・(・・、β)・ e

・ ,・

1

 (

94

α彡 1 十 α多2 

 1 α託s十 ale = 1 α釜1+α}

lals

+ai。

1

=l

aR

=−

al

 n/

an

].R  in D α,; a、. 。 。/ σ

i

。,+σ隔

 

一…

 (95 た だ し

α z、

α R<O

α z、

α D

α LS

α Le>0である

 3

2 荷重 係 数お よび耐 力 係 数の計 算ア ル ゴ リ ズム   前 項で述 べ た荷 重 係 数お よび耐 力 係 数の 数 値 計 算に は

これ らの算 定 式に含 まれ る方 向 余 弦 a (以 ドで は

分 離 係 数とい う) を効 率よ く計 算するアル ゴ リ ズムが 必 要とな る。

 

こ の計 算ア ル ゴ リ ズム は

il

α

1

1の条 件の下で Rn(a) の 大値を見い だ す制 約 条 件 付 最 大 値 問 題と し て定 式 化

WD

WL

L

s

2 R。(a)曲面 D

 αZl さ れる

ある い は

  

 

/(・)

。)

………・

一 ……・

…一 ・

……

6

と し て最小 値問題 と して取 り扱うこと がで き る

 (90)

(95) 式の例で は

    

R

。(α、 、,α。。,α埆 α Lε}       7。(α 。

);v。

。+7,

s(α

α

w

n+rLα 、

α 、) 

WL

。       φ(az

)      

 

t・

一・

tS・

一・

 (97> と な り

上式 中に は4個の 分 離 係 数が含ま れるが

aLe E 鳳 お よ びα z

=驫

の条 件 を用い る と下 記 の

2

数 関 数最 小 値 問 題と して定 式 化さ れ る。

  

 

∫(・)

……一 ・

…一 ………・

…・

(・8)  こ の数 値 計 算に は例えば

Nelder−Mead

の シン プレ ッ クス法に よ る非 線 形 関 数の小 値 化 プログラム が使 用で き る

 シンプレッ クス法に より探 索さ れ る非 線 形 関 数の最 小 値は

,一

般 的に局 所 的な最 小 値を含み

正し い最 小 値を 求め る た めに は適 切な初期値のが 必要と され る

ま た

数値 計算のら れ た最小値につ い て は そ れ が大 域 的な性 質の最 小 値で あ る か否か確か め る 必要が あ る

  図

2は

上 述の荷 重み合わ せにつ いて

最 大 値 点を含む領 域にお け る R

(a)の値の性 状を図 示し たもの で ある

こ の例で は Rn(a)の曲 面は滑ら か に変 化し て お り, 図に描か れ て い る曲 面の最 大 値 点 以 外に は, 局 所 的にも最 大 値 を とる点は存 在し ない

 

4.

数値計 算例  本 稿で は 3節に述べ 関 連す る数 値 計 算 例 を 示 す

より全 般 的な数 値 計 算 結果 とその につ い て は 続 稿で報 告する

  4

1 超 過 確率の理論 解につ い ての計算 例  次 項に述べ る計 算 例と関 連して

2

2で示 し た PSW 過 程とPI 過 程の組み合わ せ に おける超過確率の理論 解

(8)

9

ε

§

8

§

9 こ cD

 .

0 頃

O も

白 も

  0

b

 

 

 

T

50

  o   Simuhati・n

; P【OP(馮ed o 2

4

6

6

ID

( &

ε

§

§

§

Φ

 .

0

0

b

も も

0 も

o

 

 

 

T

8

  o     S正mulation

: PI。P。6ed o 2

      4 

      6

      8

      10       Thre5ho〃  r      

Threshold

 T

一3PSW

過 程 とPI 過 程の組み合わ せ過 程の超 過 確 率, 

PSW

:(W

 Vw;λ

μ)

(1

0

0

4;1/8

8)

 PI W

     V

;λ

μ)

〔0

7

0

5;1

0

0) 表

一1

確 率 過 程モデル の パ ラ メ

タ お よび統 計量 荷 重     生 起 率  作 用 時 間   平 均 値   変 動 係 数        A〆年)  μ(年)  レVx/Wx

   Vx Wp 1

0 0

1D WL   WL

     1〆8 WL

1 8 D

360

25 0

40D

55 限 界 耐 力R 110 015 ((6 )式 )につ い て モ ンテカ ルロ シ ミュ レ

シ ョ ン手 法に よっ て算定さ れ た超 過 確 率の計 算 例 との比 較 を図

3に示す

 図

3は 横 軸に レ ベ ル r を と り

縦 軸に は超 過 確 率 を対 数 座 標で図 示し たもの で あ る

実 線は 理論 解によ る超過確 率を 表しtO 印は シ ミュ レ

シ ョ ン に よる超 過 確率を表し てい る

 各 確 率過程の パ ラメ

タは

1に示 す 持 続 的 積 載 荷重と

時的積 載荷重の確 率 過 程モ デル の パ ラメ

タ と し た。  シミュ レ

シ ョ ン手 法は

文 献 1)に述べ られて い る 手 法と同じ で あ る が

こ こ で は シミュ レ

シ ョ ン標 本 数 を

1Vs=loq

個と し ている

 図

3は基準 期間

T =

50

お よび T

8の解 析 例で あ る が

い ず れ も 理論 解は シミュ レ

シ ョ ン結 果と よ く対 応 しており

理 論 解は満 足な精 度 を具え てい る。 基 準 期 間の短い T

8の例で は

シ ミュ レ

ショ ン結 果と理 論 解とはやや差 異が認め られ る

これ は

PSW

過 程の継 続 時 間μ1 と基 準 期 間の長さ が同 程度であ る ため

打 ち 切 り定 常過 程にお け る 超 過確 率の近 似 的 取 扱い に よると考 え ら れ る

な お

こ の ほ かの確 率 過 程の組み合わ せ に お け る超過確率の理論 解とモ ンテカル ロ シ ミュ レ

ショ ン 結 果との比 較 例は

文 献10 )

1上に示さ れて いる

 4

2 荷 重 係 数および耐 力係 数の計算 例  3節に述べ た固 定 荷 重 と 積 載 荷 重のみ合わせ に お け O

N ← 導 0

0

O       O 5        ]

0       1

5       2

0       2

5 ;VL

1WD

b ぜ D

O 鄲1D +Pシ毛 ; β置 2

5 T

50

yr 創L TD

P一 尸 一

凾 F − P

網一一

O

0       0

5

0      】

5     2

0     2

5 xrL

IWFD

4 荷 重係数 お よび 耐 力 係 数 (終 局 限 界状態

固 定 荷 重+積 載 荷 重 )

一 38 一

参照

関連したドキュメント

状態を指しているが、本来の意味を知り、それを重ね合わせる事に依って痛さの質が具体的に実感として理解できるのである。また、他動詞との使い方の区別を一応明確にした上で、その意味「悪事や欠点などを

状態を指しているが、本来の意味を知り、それを重ね合わせる事に依って痛さの質が具体的に実感として理解できるのである。また、他動詞との使い方の区別を一応明確にした上で、その意味「悪事や欠点などを

An example of a database state in the lextensive category of finite sets, for the EA sketch of our school data specification is provided by any database which models the

Nonlinear systems of the form 1.1 arise in many applications such as the discrete models of steady-state equations of reaction–diffusion equations see 1–6, the discrete analogue of

This paper presents an investigation into the mechanics of this specific problem and develops an analytical approach that accounts for the effects of geometrical and material data on

Then the strongly mixed variational-hemivariational inequality SMVHVI is strongly (resp., weakly) well posed in the generalized sense if and only if the corresponding inclusion

 

Amount of Remuneration, etc. The Company does not pay to Directors who concurrently serve as Executive Officer the remuneration paid to Directors. Therefore, “Number of Persons”