NII-Electronic Library Service
K
報 制 UDC :624.
042.
2 :519.
2 日本建 築 学 会 構 造 系論 文 報 告集 第 393 号・
昭 和 63 年II月荷
重
の
確 率 過 程
モ
デ
ル に
つい て
正 会員坂 本
順
*1
.
概 説単
一
の確 率 過 程 あるい は組 合せ確率過 程は,
多く の 工 学 分 野に おいて応用され てい る。
構造工学分 野で は, お もに不 規 則な動 的 負荷に対す る動 力学問 題の確率論的 解 析と その応 用につ い て多くの研究が なされて き た。 近 年,
社会施設や システ ム の確 率論的リス ク解析, さらに一
般 建 築 物や海 洋 構 造 物の確 率 論的構造 解析お よび設計 法に関連する問 題と して, さま ざまな荷 重 作 用 (loadactions ,
load
effects)と そ の組合せ につ いての統 計・
確率 論的 取り扱い に関心が注が れるように なっ た
。
例え ば, 地 震 荷重,
強 風荷重,
雪 荷重,
波浪荷重などは, そ れ らの生 起 特 性, 強さ, 作 用 時 間な どに不 規 則性・
時間 的変動 を伴う ものであ る こ と か ら, 確率過程 (stochastic process )として モ デル化さ れ る。
これ らの荷重過程は,
振 動 問 題に お け る時間連 続の確率過 程 (continuous tirne process)と異な りマ クロ 的 時間変動 を伴う確 率 過 程と し て特徴づ け られ る。Fig.
1−
1に示 す 再 生 矩 形パ ル ス過 程 (renewal rectangular pulse process), ボアソ ン矩 形 波 過 程 (Poisson rectangular wave Prosess)
,
ボア ソンパ ル ス過程 (Poissonpulse
process
)な ど は そのような確 率過程モ デル の例である
。
こ の ような荷 重の確 率過程モ デル に 基づ く構造信頼 度 解析は,
下記の ように記 述 され る。
統 計 的に独 立なm 個の荷 重の定常確 率過程をW
,(t
) (」=1
,2
,……
, m ) と して これらの組合せ確 率過程をW
(t)と す る。w
(t
)=
w
,(t
)十Wz
(t
)十・
・
…・
十Wm
(t
)・
・
……
(1−1
} こ の組 合せ確 率 過 程に対 する限 界 状 態 関 数をG
(R ,W
(t
);T。
)と する と,
構 造 信 頼 度 問 題は下 式の よ うに定 式 化さ れ る。
G
(R ,w
(t
);Ts
)=R −
1
臥(t
)+w,(t
}+肱 (t}};Ts…・
………・
……・
……
(1−
2)Ps=P
。。b[R −
lva
(t
)十Wa
{t
)……
+W .
{t
}1
≧0;Ts]・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(1−
3) こ こ に,
P
。=
構造信頼度R
黒
構 造 抵 抗Ts=
構 造物の使 用 予定期 間 P,
。b[E ]=
事象 E の生 起 確率 〔1−2
), (1−3
)式に おいて荷重の組 合せ過程W
(t
)を 確率変数に置き換え る手法が適用で きると信 頼 度 解析の 数理的 取扱いが容易とな る。 組 合せ確率過程W
(t )の時間 区 間Ts
における最大値 を表 す確 率 変 数は下 式で定義さ れ る。
W
,。=
maxlW
、(t)+vr
,(t)……
+ va。{t)1
・
一
(1−
4 > Ts こ の確率変数 貼。
の統計的 性 質が得られ れば, (1−
2>,
(1−
3)式は確 率過程 剛置)を陽に含ま ない で下式の よ う に記述で き る。
9(
R 、iVTs
)=R − WTs・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
tt・
・
・
・
・
…
(
1−5
)Ps;Preb
[R 一
陽s≧0
] x一
ル
肱 ωd
蜘・
…・
…………
(1−
6 ) 訓・
1≧ oCont}nUOUS τime proceSS
i
』
」
_
ll
−
m
⊥
一 .tRenewal rectangylar PUtse ProceSS
L
[
L
[
L
,
[
k
.ll
=
[
u
]
.[
b
.
.
.tPDiSsoll 「eCtangufa1
}
waVe p「oCeSSL
』
_砠
一
」
L
,掌 名 古屋 大 学 教 授
・
工博(昭和 63 年 3 月 23 日原 稿 受 理1
POiSSOn l1UISe P「OCeSS
Fig
.
1−
1 Schemahc representation of stochastic Processes一 62 一
こ こ に
,
fx
(x)は確 率 変 数 X の確 率 密 度 関 数 (以 下, p.
d,
f.
と略記 〉である。
(1−
4 )式に よっ て表さ れ る確 率 変 数 Wr。
は,
確 率 過 程W
(t)の 「Ts期 間 最 大 値 」 と 呼ばれ る。
こ の確 率 変 数 貼 。の確 率 的 性 質 (P.
d.
f.
もしくは c.
d.
f.
(確 率 分 布 関数 } を求め る 問題は,
確 率過 程の閾値横 断 (barrieI
crossing )および初通 過 時 刻 (
first
excursion time)の問 題 と 関 連づ け られ る
。
定 常確 率過程 (時間 連 続 過 程)X
(t
)が,
時 間 (t
、,
竕 に おい て閾 値 ξを 上側から,
あるい は下 側か らの い ずれ か か ら横 断する横 断回数を求め る問 題は,Fig.
1−
2 に示さ れ る よ う に不規則過 程Z
(t),
閾 値 横 断 位 置にお か れ た単位 イン パ ル ス過程 Z (t)を用い て時 間領 域(tl
, の に おけ る計 数 過程 娯ξ;tL,
t
,)の統計 的性 質を取り 扱う問題と な る (Lin,
Y .
K .
)]}。 確 率 過 程X (t)の単 位 時 間 当た り の閾 値 横 断 率の期 待 値 (mean crossing rate )は,
下 式で表さ れ る。
v(ξ)
一
∫
:
1thlf
・}(ξ,・tl
,・・}dth−
……
(1−
7)こ こ に, ル熟コじ.
ic
)は,X
とX
との結 合確 率密 度 関数である
。
と くに
,
下 側か らの 閾値 横 断のみを取り扱う と き に は,期 待 横 断 率 (mean upcrossing rate )v+
(ξ)は下 式 とな る。
・・(ξ)
イ
繍
(ξ,
±,t
)dtl …………
(1−
8) 定 常 確 率 過 程X
(t
)の (O
,T
)時間に お け る最 大 値の c.
d.
f.
Cま,
F
x、
Me(ξ)=P
γ。b[X
(t}一
ξ≦0;T ]……・
…
(1−
9 ) T ま た,
X (t)が閾 値 ξ を下 側か ら最初に横断す る時 刻を 表 す 確 率 変 数 をTf
とす る とFmax.
mt)(ξ}=
1−
Frr(T ;ξ)・
・
一・
・
・
…
tt・
一…
9・
(1−10
) T こ こ に, F.ノ(・
)は, 初通過時刻の確 率分布関 数である。 上 述の よ う にX
(t
)の最大値の確 率 分 布 関 数 を 求め る 問題は,
閾値横 断問 題と関 連づ け られ る。
しか し,一
般 的な確 率 過 程に おいて (1−9
), (1−10
)式の厳 密な解 析 x【t)9
ξ z1
同 lb, 2【己〕 り 蝿c ,Fig
.
1−
2 Barrier clessingt は数 理 的に難 しい問 題 とな る
。
その数理的取り扱いの因 難さの理由か ら,
しばしば近似 的手法が用い られ て い る。
閾 値 レ ベ ル ξが比 較 的 高く, それ を横 断す る事 象が ボ ア ソ ン生 起 で ある とみ な さ れ る仮 定を 用 い る と,
max.
X
(t
)の c.
d,
f.
は,
近似 的に下 式に よっ て推 定さ T れ る。F
x.
M跡(ξ)= expl−
1/1
(ξ)TI………
……
…
(1−11
) T マ クロ 時間変 動過 程 と して モ デル化される荷 重 過 程に つ い て も上 述の近 似 的 取り扱いが適 用さ れ る。
本報 告で は
,
荷 重の確 率過程モ デル と 上 述の閾 値 横 断 率お よ びTs
期間 最大値の c.
d.
f.
の 近 似 推 定 法の適 用 性につ い て検討を加え,
あ わ せて若 干の応 用 的考察につ い て述べ る。 本 報 告の内容は, 確率論 的構造 解 析あ るい は設 計 法に お け る荷重の確率過 程モデル につ い て の基 礎 的資料を提示す る もの で ある。
2.
荷 重 の 確 率 過 程モデル マ クロ 時 間 変 動 過程の代表 的な確率 過程モデル と閾 値 横 断につ いて 述べ る 。 こ こ では, 荷 重の確 率 過程モ デル を対象と して非 負の値 域 を もつ 定常 確 率過程を取 り扱 い, また, 1つ の確 率過 程 にお け るパ ル ス などの生 起, 強さ,
持続 時間は 互い に統 計 的に独 立 とす る。
単
一一
の 荷 重 過 程W
(t
)の 閾 値 ω* につ い て の期 待横 断率 vS (ωつ は,
下式で定義さ れ る。 ・あ(げ )一撫
毒
堀 [腕 )< w ・,
w (t+At )>w・ ]一 ・
………・
・
……
(2−
1) 上式の 喝(w ’ )を 用い る と 確 率 過程W
(t)の Ts期 間最 大値の c.
d.
f.
は,
近 似 的に F.
曜(ω)乞 expl
一
レ轟(w)Ts
}・
……一 …・
・
(2−
2) Ts と して推 定され る。
A
) 再生 矩 形パ ル ス過 程再生矩 形パ ル ス 過 程 (renewal rectangttlar pulse process
,
以下で は,
R−
R−
P 過 程と略 記 )は, 下 記の よ う に構成さ れ る (Larrabee
,
R .
0 .
,Cornel1,
C .
A .
)m。
Renewat poinl proceSS
劇
P
・
’
f〔w)Fig
.
2−
A Renewal rectangular pulse pToeess1) パ ル ス の 始 点 と終 点 は
,
再 生 点 過 程 (renewalpoint
process)と して規 定さ れ る。 再 生 点の期 待 生起 率 をλp と す る。
NII-Electronic Library Service
2
) 各パ ル ス の状態確 率をp, q と し て, 確 率 q でパ ル スが生起 (‘ on ’ 状 態 )し,
確率 p で ‘off’
状態を と る。
off状態で は, パ ル ス強さを零 値 と する。
3) パ ル ス (実パ ル ス〉の期 待生 起率は λα=
q・
λρ,
期 待 持 続 時 間は μd=
ユ/λp であ る。4
) 任意 時刻に お け るW
(t
)の p.
d.
f.
は,
fw
〔w}=P
δ(w)十qfw
(w)・
・
・
・
………・
…・
(2−A −1
)こ こ に,
fw
(w)は パ ル ス 強 さのp.
d.
f,
で あ る。(
2−A −1
)式のfw
(ω)は, 周 辺 分 布 (marginalp.
d
,
f
,
〉 とも呼ば れ る。R −R −P
過 程の 期 待 横 断 率は,
(2−
1 >式 よ りパ ル スが 零再帰す る と き と非 零 再 帰のときにつ い て下 式で表さ れ る。 yV (げ)一鵯
士
・
ガ
ゐ
ω・
A
。・
At ・
i1
−
ii
・(・V’}ldw
=
Xa“−
FAtL“)1
;零再帰のと き…・
(2−A −2
}=
ip
十q・
F
賦w’)トλall
− F
貳げ )i
;非零再帰のとき
…
(2−A −3
)B
) ボア ソン矩形波 過 程ボ ア ソ ン矩 形 波過 程 (
Poisson
rectangular waveprocess;
P −R −
W 過 程 }は,
指 数 分 布に従 うパ ルス生 起間隔丁。お よび持 続時 間丁,に よっ て構 成 され る
。
E
[Tti
]=
=
1/λα ,E
[T
』]= Ptd=1
/λa・
……
(2−B −
1)こ こに
,E
[・
]は期 待 値 を表 す。
瀏
M
lreturn te te【ol
Fig
.
2−
BPoisson
rectangular wave pfocess矩 形パ ル ス の強さの c
.
d.
f.
をFw
(ω)とす ると 期 待 横 断率は,
下式で表され る。
レ壱(w“ )
=
A
αFw
(tv * )11
− Fw
(w * )}………
(2−B −2
) 上式は, (2−A −
3)式におい て,p
= 0,
q= 1に対 応する。
各 矩 形パ ル スが零 再 帰の場 合に は,
(2−A −
2>式と同 様に,
レ訊”つ=
λα
{1− Fv
(盟*)1
………・
…・
・
・
…
(2−B −
3)C
> 矩形パ ル ス過 程矩形パ ル ス過 程 (rectangular pulse process;
R −P
過程)は,
一
定の単 位 時 間 τ ご とに生 起す る規 則 的なパ ル ス列で構成され る。
単 位 時 閻を τ;
1と す る と, 単 位 時間当た りの期 待 横 断 率は下 式と な る。 y高(ω*)=
1一
凡 (w*)・
一 ・
…・
・
……・
…一
(2−C −
1) なお, パ ル ス の持 続 時 間 ρ己は,0
〈P
.≦ 1と する。
^
記 号は確 定 変 数 を表 す。 D > ボア ソ ン パ ル ス過程お よび ボア ソ ンイン パ ル ス過 程一
64
一
血 曲
[
1
,,[
i
,[
L
図
Fig
.
2−
C Rectangu蓋a【pulse processボア ソ ン パ ル ス過 程 (
Poisson
pulse
process
;P −P
過 程 )は,
一
定の単 位 期 閾 をτ と して持 続 時 間魚,
期 待 生 起 率λのボ アソ ン生起の パ ル ス過程と し て構 成さ れ る。 ボアソ ン イン パ ル ス過 程 (P
∠工過 程)は,P −P
過程の パ ル スを イン パ ル ス に置き換え た もの であ る。
P−P
過 程お よびP −1
過程の期 待 横 断 率は, ン鼠ω)冨
λ{ユーFw
(Of
)1
・
…・
………・
・
……
(2−D −
1)劉
’L
一
L
一劃
Fig
.
2−
D Poisson pulse process and Poisson impulse pfocessE
)泝 波 され た矩 形パ ル ス過程
〜戸波 さ れ た矩形パ ル ス 過 程 (
filtered
Tectangttlarpulse
process;F −R −P
過程 )は, 各パ ル ス の期 待生起 率 λα
, 期 待 持 続 時 間μd に よっ て構 成さ れ る。
パル ス の 生起お よび持続時 間は, 指 数 分 布でな くて も よい。 こ の 確率過程モデルの期待横断 率は, 近似 的に下 式で推定さ れ る (Wen , Y.
K.
)。 レ存(が )=hll − Fw
(w ’ }i
’
・
……・
・
………
(2−E −1
)劃
Fig
.
2−
E F藍lte[ed recta皿gutar pulse process上 記の各 確 率 過 程モデル の
Ts
期 間 最 大 値の分 布 関 数は, 上述の期 待 横 断 率y占(w )を (2
−
2)式に用い て近 似 的に推 定さ れ る。な お
,P −
R−
W 過 程な ど で初期時 刻状態 を考 慮する場合に は,
Fmax
.
w(w )= oFv (w)exp1一
ソ缶(w)Ts卜……・
{2.
−
3}7s
こ こ に
,
oFw (ω)は初 期 時 刻にお ける非 超 過 確 率を表す。 また,R−P
過程に つ い て は,
Fmax.
v(w )zIFw
(w)iN
;ハ1=
Ts/τ,
integer,
Ts
・
…・
・
………・
(2−
4 >:expi(1
−
F
試ω})Ts
};τ=
1の と き・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
一・
・
・
・
・
・
・
・
・
・
・
…
(2−
5) つ ぎに,
組 合せ確 率 過 程の期 待 横 断 率につい て述べ る。
統 計 的に独 立な再 生 矩 形パ ル ス過 程W ,
(t)と 既(t)の ス カラー
和の確率過程をW
(t
)と す る。W
(t
)=
蝋 孟)十W
,(t
)……・
……・
…tt・
一 …・
(2−
6) 確率過 程 W (t
)の 期 待 横 断 率は, 下 式で表さ れ る (Larrabee,
R .
O .,CQrnell,
C .
A .
}zl。 v’ ”・’)一
蛎
置∫
°
°
∫
ガー
w’f
・ ,(・Vl)ノ
.,(・v,)。
1
λP,
・
At 。
(1− Fw ,
(w*−
Wz ))一
トλρ2・
△t・
(1− Fw
、(w *−
Wl))}d
.(∫w、・
・
一・
・
・
・
・
・
・
…
(2
−
7>∫
Utノ
・ ,(・)・
・朝(・f−
w)・dw
+f
, ”ノ
.,
(・)・
〃畆(げ一
ω)d
ω………・
・
……・
…
(2−8
) こ こ に,
fWi
(ωD
は 隅 (t)の 周 辺 分 布 (皿arginalp.
d.
f.
), 殉 (w )はW
,〈t
)の期待 横断 率で あ る。 (2−8
>式は
,
点 横 断公式 (point−
crossingformula
) と も 呼ば れ る。
(2−8
)式 右辺の積分項は, 密度関 数と期待 横断 率 関 数 の 「畳 込み積分」で ある。 隅ω け僻
1,
2)が零 再 帰 過 程の と き 同 式の 積分 を実 行す る と下式が導か れ る (Fig.
2−
1 参照〉。の
星
{L一
鳳 げ 吻 }}鳳 吻) F臨 〔ガ〕−
Fwけ、
【Ψ
L} ノ一
齢
一
一
一
一
『一
吻 1 旨ぴ
}■
丗
21吻 +d吻ψ
罪
∫耐呵Fig
.
2−
1 1ntegral domainレ訊ω*)
=
=
λai(Pe一
λ殉μ匸}F
篇1(ωつ十λ鈎(Pi一
λaiPt2)F
魯2(w * )→
一
λaiλα2(μ1十 μt>F
昏!や鐸t(ω ホ )}・
・
・
・
・
・
・
・
・
・
…
(2
−
9) こ こに,
F
叡ω)=
=
1− F
騨(W):補 分 布 関 数,
F
賑 駆ω)は,
和の確率分布関 数を表す。3
個 以 上の組合せ 過程の場 合 につ い て も期 待 横 断 率は 「畳 込み積 分 」に よって算 定 さ れ る。 し た がっ て,
確率 変数 隅 の特 性 関 数が存 在す る 条 件の下でfWj
(w )の フー
リェ 変 換が適 用で き る な ら ば,
期 待 横 断 率の算 定にお ける 「畳 込み積 分 」は容 易 化 さ れ る。
また, 前 述のF −R −P
過 程の和の確率過程の期待 横断 率につ い て は,
W,(t)(ノ亟
1, 2……,
m )が統 計 的に独 立の 場 合に は下 式の 近 似 評 価 法 が 提 案さ れ てい る (Wen , Y.
K.
> 3 )。
W (t)=
W,(t)十 隅(t
)+……
+ Wm(t
)・
……
(2−
10 ) m れ m レ藁(げ ):Σλ,F魯,(ガ1
十Σ]Σλ‘,F 昏1田 ,(げ) ’屆
王 ‘呂
且∫=
l t申 丿 m m nt 十Σ Σ Σ λ雛 F 昏‘+WJ+Wit(w*)十……
t!
ti=
kロ t幸丿幸k
・
…
一…
一・
一・
・
・
・
・
・
…
9・
・
・
…
(2−
11) こ こに,
λi=
λ画 λt」ニ
λ伍 λω (μdt 十 μdi),……
前述の (
2−9
)式に おい てp
=1,q
(=
λ。μ)= Oの と き に は,
(2−
9)式と (2−
11)式は同 じ評 価 式と な る。
組合せ確 率過 程のTs
期 間 最 大 値 分 布につ い て も上 記 の期待横断 率 を用い て (2−
2)式あるい は (2−
3)式によ り近似 的に推定さ れ る。 以上,
代表的な確率過 程モデル につ いて述べ た が,一
般 建築物の設計に用い られる荷重で は,
持続 的性 質をも つ 積 載 荷重な ど はP −R −W
過 程と し てモ デル化 さ れ,
気象現象とし て年 周 期 的に生 起す る雪 荷 重や 風荷 重は R−
P 過 程 とし て , また, ある強さ以 上の地震動を対象 とし た地 震荷重や強風荷重は,P −1
過程ある い はP −P
過 程と して モデル化さ れ る。
3
.
モンテ カル ロシ ミュ レー
ション による考 察前 節 に述べ たマ クロ 時 間 変 動 過 程の 期待 横断率茜 (w )お よ び
Ts
期 間最大 値 分布F
w(w )の近似評 価 法 TS の適 用性につ い て モ ンテ カル ロ シミュ レー
シ ョ ン手 法に よる検 討 資 料 を以 下に示 す。
モ ンテカル ロシ ミュ レー
ショ ン に よ る数値 解析は,
代 表 的な単一
確 率 過 程 と組 合せ過 程につ い て次の手 法に よっ て行 われた。
1) 本 稿で は,
確率過程W
(t)の荷重強さ を表す確 率変数
W
が対 数正規および極 値1
型 分 布に従う場 合の シ ミュ レー
ショ ン例につ い て示され てい る が,
対 数 正 規お よび極値1
型 乱数はBox ・
Muller
法に よ り生 成さ れ た 正規乱 数 を変換 して生 成 した。 2) 各 確 率 過 程の シ ミュ レー
ショ ン は無次元時間に よっ て行い, 各標 本関 数の無 次 元 時 間 長さをT
。=100
と して い る。
3) 標 本 関 数は, 各確率過程ごとにIO3
個を生成し,
期 待横 断率は点横 断 (point・
crossing )に基づ く計 数 法を,
Ts
期 間 最 大 値 分 布は平均ラン ク法 (Hazen
法)を適 用し て算 定。 最 大 値 分 布は,T
。;10
お よびT
。=50
の 2例につ い て計 算し てい る。
4 > シ ミュ レー
ショ ンに お け る基本パ ラメー
タ は,
確率 変数W
の平均値W =1.
0,2.
0,
変動係数V
− nO.
2,
O.
4,
0.
8,
パ ル ス生 起 率 λ=
O.
1,
02,
0.
3,
1.
0お よびこ れ らの組 合せを設 定してい る。 期 待 横 断 率 v+ (w )お よ びT
。期 間 最 大 値 分 布FWrs(ω ) の シ ミュ レー
ショ ン数値計算 結果と前節に述べ た評 価 式 に基づ く計 算 結 果とを 対 比 して Fig.
3−
1〜
Fig.
3−
5に示 す。 な お,
以 下で は,
前 者を シ ミュ レー
ショ ン解 (図 中NII-Electronic Library Service 冒α レ T (1 ) satspl■tunction
凾
ε ≧ …〕
B W 三s ; ヲw
●
e δ囗
UPCrOSS■
叩 「a ヒew c d 「 o「 Ts
.
囮ali靄u日P
三(
三,
wFig
.
3−1
Poisson irnpulse precessw
(
εこ
言即
、 w ■eanけpcrOSSI叩r 乱 o 冨 苓 3}
2 s繍 卩Ief恤net
匹
o 隠 零卩
O
【
, 冨w 鰓3n
u匹ros5i臆9 「凾te 脚 c
.
d l or Is・
ma罵1踟 墜 HFprOClsso
.
o「
耄 ま 9 9 tO■
0 5層
,
四
Cld 「 or IS」
蘭 民脚酣 臼 w.
図
Fig
.
3−
2 Rectangu且ar pulse p【ocessw OQ に破 線で示 す)
,
後 者を 理論 解 (図中に実 線で示 す } と 略 記する。
Fig.3−1〜Fig.
3−5
に示さ れて い る シミュ レー
ショ ン 資 料の性 状は,
下 記の よ う に要約さ れ る。
1
)Fig.
3−1
お よ びFig.
3−2
に示すP −1
過 程お よ びR −P
過程の期 待横断 率 とT
。期 間 最 大植分布の理論解 は近 似 を含 まない こと か ら,
両図は le3個の標本関 数 に よ る シ ミュ レー
ショ ン の 精度と有 効性を示し ている。
Ts
期 間最大 値分布の上 側 裾 野 域に お い て標 本数 の不 足に よる シ ミュ レー
ショ ン解の不安 定さが 認め られ る が
,
ほ ぼ満 足に シ ミュ レー
シ ョ ンが行わ れて い る と 言え る。Ts
期 間 最 大 値の 累積分布 値0.
999,0.
99
およ び0.9
は,
標 準正規 関数 φ(β)の β 指 標 値で表す とβ=3.
05
,2,
32
お よび1.
28に相 当 する。 シ ミュ レー
ション にお け る標 本数の制 約お よび
Fig.
3−
1〜Fig.3−2
の シミュ レー
ショ ン結果か ら, 以下の考察はβ :1.
00−
2.
50の範 囲を前提に し て述べ る。
2) P−1
過 程お よびR −P
過 程 と も荷 重 強さ W が対 数 正 規 分 布に従い,
かつ 変動係数がVw=
0.
2, 0,
4の場 合に は,T8
期間最 大値 分布は2
重 指 数 確 率 紙 上にお い て ほ ぼ直線分布と な り,
近 似 的に 2重 指 数 分 布と し て取り扱っ て よい こと を示し てい る。
しか し, 変 動 係 数が大きい 砺 =0.
8
の場 合に は , 2重 指 数 分 布の近 似的取り扱い は適 用で き な く なる。
冖
塁【
ε,
」
蝋 t 三k 5a薗
Ple ru旺
匸ion 辜゜
w 齠an
UlrcraSSlng
riLe
冖
三 、P
ら まW
W
w cdl Dlis
・
田∂篤
1明
●Fig
.
3−
3 Poisson【ecta夐gular wave Process荷重強 さ
W
が極 値1
型 分 布に 従 う 場 合に は,
Fig
.
3−2
に例示 さ れてい るよ うに W の変 動 係 数の値 にか か わ らずTs
期間最大 埴の 累 積 分 布は,
2重 指 数一
66
一
分布に ほぼ適合す る と言え る。 3)
Fig
,
3−3
に示 すP
−
R
−
w
過 程につ い て も理 論 解 とシ ミュ レー
ショ ン解は良い対 応を示 して お り,
上述のP −1
お よびR −P
過程の性状と同様な指摘が な さ れ る。 v[t) P刪←
円’
n腱
ヒ
“
賢
;”λ
’
°.
21ぬ
人【
’ ヨ 繍 黠 髭:,
サ
:1
.
1
O
一
「
ε 、.
O
向
露
; 8、
胛
膾
”
.
・
檜
.
麗
闇
F
鬮
、
幽
.
冨
1囲
【
5a踵
ple 「un匸
tio隠 壱90
一
.
6四
鹽
、
〇
一
「
ε 、 0囓
。
… 12.
, SY Iv』ean up匸rossi叩 ratc
.
, 言 } [dr oIls
.
tuaXtEU勵 loヤ
O一
▼
O「
畠
ε 、 ε−
ま P−
a一
鬨 λ・
02 曜o▼
FO
‘■
閃
n・
脂o P・
l▼
・
A・
01co 05盻
●
n甌
10 PR賢!
pnVrPI 門 ].
Pt●
.
願
■
’
tt・
o w Fig.
3−4
PRW 十PI process襁
.
■
鹽
・
1醍
L,
遍
・
”
”
闢
圏
卩
轄
5盈■P52 「u眠監嚠
Oh w一
こ、
円
と ど w・
罰
「
三 ま鬮
te
W
国
oanvp
[
ro5s十
no「
nle“
ε 』「
ε 塁 w w cd.
r oIlS一
日ax旧U輌 Fig.
3−
5 PRW 十RP process4
)Fig.3−
4およびFig.
3−
5 は, 組 合せ 確 率 過 程の 例 と してP −R −W
過 程 とP −1
過 程 あるいはR −P
過 程と の組 合せにつ い て の計 算 例 を図 示 した もの で ある。
両 図中の理 論 解は,
P
−
R
−W
過 程が常にbn ’
状 態にある もの と し て 〔2−
8)式の 考え方に基づい て近似 的に下 式に よ り算定し て い る。P −R −W
過 程 とP −1
過程の組合せの場 合; レ存(w*)= (λ蘆一
λ2)F
竇1(w * )十 λ2F 昏、+w3(tV * )・
・
・
・
…
(3−
1)P −R −W
過 程 とR −P
過 程の組 合せの場 合;レ存(wi)=
一
λ2F
昏1(ωf
}一
λ,F
昏2(ω1
) 十〔λ1→一
λ2)FI
レ匸+ w2(ω零)・
・
・
・
…
噛
・
・
噛
・
・
…
幽
・
…
(3−
2 ) こ こに,
添字1
お よび2
は,P −R −W
過 程お よびP −1
過 程 あるい はR−
P 過 程の を表 す。 な お, (3−
2 )式に おい て λ2=LO
C
記 号は確定 値を表す)とす る。Fig.3−4
お よ び Fig.
3−
5 で は,
変 動 係 数が大きい場 合が示さ れ ていること も あるが,
単一
の確 率 過 程の場 合 と比べ て理 論 解と シ ミュ レー
ション解との差 異が大き く なっ てい る。2
つ の組 合せ過 程に おい て各 過 程の パ ル ス の生起 率, 継 続 時 間が小 さ く, かつ 零 再 帰の場 合に は (2−
9 )式 あ るい は (2−
ll )式は,
期待横断率の正 しい評 価2)あ るい は良い近 似 評 価7)を 与え る 。 これらの評 価 式に準 じ た (3−
1), (3−
2 )式の近 似 評価は,
別稿で述べ る限界状態.
設計 法に おける荷重 係数の実用的な設定 法と して適用 す る こ と を 意 図し た もの である が, 非 零 再 帰のP −R −W
過 程を含む組 合せ過 程で は,
シ ミュ レー
ション解との差 異 を 生 じる傾 向がある。
また, こ こで示し た シ ミュ レー
ショ ン例の ように変動係数が大きい組 合せ過 程では,
シ ミュ レー
ショ ンにおける標 本 数 を さ らに増 して考 察 を加 え る必要が あ り, これ らの点につ いて は,
さ ら に詳 細な 考 察を加え稿を あ ら ためて報告する。4.
応 用 的 考 察 前節ま でに述べ た荷重の確 率過程モ デル の期待横 断率 と そ れ に基づ くTs
期間 最大値分布の応用的考 察と し て 確 率 論 的 限 界 状 態 設 計 (probability
−based
limit
statesdesign
;文 献8
),9
))へ の応 用につ い て概述 し て お く。
単一
の荷重の 確率過程W
(t)につ い て指定さ れ た期 待 横 断 率 φ(一
βw)に対 する閾値 勿は, 下式によっ て算 定さ れ る。
レ壱(勿)=
λ11
−
Fw (勿}}=
φ(一
βw)………・
…
(4−
1) こ こ に,
φ (・
)は標 準正規 関 数で あ る。 確 率 変 数 W が対 数 正 規 変 数で ある場 合に は,
荷 重 過 程 W (t)の 閾 値 {b
は, (4−
1 )式よ り下式の よ うに表さ れ るQ ・・一
論
・
exp{
・一
見(
i−
1
・(−
R
・))
al・ ・}
・
W
・
・
・
…
+
・
・
・
・
・
・
・
・
…
一
・
・
・
・
・
・
・
・
・
・
・
・
…
(4−
2 > こ こ に,
W および 砺 は,
W の期 待 値および変 動 係 数,
σt. 一・
=
viii
(i
’
Ft71JT
,
Φ一
1← )は標 準 正 規 関 数の逆 関 数でNII-Electronic Library Service ある
。
(4−2
>式に おいて単位 時間に年 (year
>を 用い, かつ,
Φ(一
βの己1
/50,
(X
?w=2.
05
)と す る と,
この時のdb
・
=
加5。は荷 重 過 程 W (t)の 50年 再 現 期 待 値を表す。単
一
の荷重 過 程W
(t
)に よ る 荷 重 効 果S
(t
)が変 換 係 数 c (こ こ で は非負の確 定変数 と す る)に よっ て下 式 で表さ れ る とする。S
(t
)=
c疊
W
(t
}・
・
r・
rr・
・
・
・
・
…
r…
−rr9
− ・
−r9・
r・
…
(4−3
> あ る1
つ の設定さ れ た限 界 状 態に対 する構 造 抵 抗 を rs,
荷 重 効 果S
(t
>の閾 値 rsにつ い ての期待横 断 率 を レ叡rs)と し て , レ喜(プs〉≦ φ(一
βs)・
…・
…………・
……・
・
…・
…
(4−
4) の設 計 条 件 を設 定す る と 下式が得ら れ る。
rs ・論
… p{
ll
−
・(
1一
去
・←Bs
))
a・…1
…Vi
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
¶
・
・
一・
・
・
・
…
(4−
5) あ るい は,
rs ≧ C・
γ(β.
)。
Wn・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(4−
6) こ こに,
φ(一
βs)は,
限 界状態に対する指 定 され た閾 値 横 断 率,Wn
は荷 重の公 称 値である。
(4−
6)式中の γ(β。) は, 指定さ れ た期待 横断 率によっ て決定 され る荷 重 係 数 (load
factor
>を表 す 。 ・ 一論
… p同
1一
去
・・一
・。1
)
剃
・
毘
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
t・
・
・
・
…
(4−
7 ) 構 造 抵 抗が確 率 変 数Rs
の場 合に は, (4−
4)式 中の期 待横 断率略(r。)を条 件付期待 横断 率略。
。
(rs)に置き換 え る と,
∫
鱒
軌 )漏(r・)・d
・rs・・〈−
fi
・1
・
・
一 ……
(4−
・)f
−
Fs
(砿 (r,)drs
≦1一
去
・(−
Bs
)・
・
……
(4−
9) こ こ に,fRs
(・
)お よびFs
(・
)は,Rs
の p.
d.
f.
お よ びS
の c.
d.f
を表す。 例え ば, 確 率 変 数 W お よ びRs
がと もに対 数 正 規 分 布 に従うときに は, (4−
9 )式よ り下記の荷重・
抵抗係 数設 計 方式の設計条 件式が導か れ る。
φ(β3)・
Rsn
≧ c・
γ(β8)・
Wn・
・
・
・
・
・
・
・
…
一・
・
・
・
・
・
…
《4−
10 ) こ こ に,・
一
詳
π
・
exp{
−
aRs・
・ “(
iTl ・(−
Bs
))
・
耐
驫
…・
一 ・
…・
・
………・
…
緯 ・・低)
一
纛
・
exp{
a・・
Φ一
’(
1−
i
・(一
・s))
・
嗣
・
幕
……・
・
……・
・
………・
…・
2・ αk ; の.
、 、/ σ葦。
tS +σ言nW , aw = a、。
W/ ai。
n。
+σi
。
w・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
一・
・
…
一・
・
・
・
…
(4−
13>一
68
一
φ(β3)および γ(β。
)は, 抵 抗 係 数 お よび荷 重 係 数, aR、
およびaw は, 分離係数である。T
。期間最 大値のc,
d.
f.
を用い る構 成で は,1
つ の設 定さ れた限 界状態に対する構造抵抗を r.
,
Ts期 間 中に おいて荷重 効 果S
(t
)が 構 造 抵 抗 ru を超 過し ない確率 を φ(βu>に指 定す る と,
設計 条件は下式の よ うに記述 さ れ る。Fman
【.
s(rti)= expl一
λ(1− F
』(ru)〉・
Ts}
≧ φ(βu} Ts…
一
・
一・
・
・
・
・
…
一
・
・
t・
・
・
・
…
一
一
・
・
…
(4−
14 > 確 率 変 数 W が対 数 正 規 変 数の と きに は,
(4−
14)式よ り下式 が導か れ る。
ru ・論
rexpl
φ一
匸(
1+纛
1
・φ輪1
厨
・・
iv
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(4−i5
) あ るい は, ru≧c・
γ(βの・
レ鷺ゴ・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
9幽
9・
・
・
・
…
(4−16
> こ こ に, ・…1
−
th
・
exp{
φ一
1(
1+ 、振
・・φ嚼
嗣
・
藷
…・
………・
……・
t…・
……
(4−
17 )な お,(4
−
17)式は,T
。=1
,ln
Φ(β.)=ln
}1一
φ(一
βu)1
=一
φ(一
βu),
βu=
β 。と すると (4−
7)式 となる。 構 造 抵 抗が確 率 変 数Ru
の 場 合に も (4−
10 )〜
(4−13
)式 と 同 型の設 計 条 件 式が導か れ る。
単一
の荷重 過程の例につ い て上に述べ た よ う に,
限 界 状 態 設 計 法にお け る設 計 規 範の構 成には, 限 界 状 態に対 す る閾値横 断率と 恥 期間非超 過 確 率に基づく2つ の構 成 法が考え ら れ る が,
こ れ に関連 す る問題 を付 記し て お く。一
般 的な建 築 物の限 界 状 態 設 計 法では,使
用 限界状 態と終局 限界 状 態に対する設 計 規 範が設 定され る。
前 者 は,
主と して建築物の常時の正常な使用状態に お け る機 能 性・
使 用 性 (serviceability )と関連する が,
そ の設 計規範の構成に は, 単位時 間ある い は任意時刻に おける 荷重の統計調査資料 (気象現象な ど に よ る荷重の年, 月 ある い は日最 大 値,
任意の時点にお ける積載荷重の調 査 資料など)をその まま利用で き る簡便さと使用 限界 状 態 の性格と を あ わ せ考えて閾値横 断率を用いるのが実用 上 便 利であ ろ う。 終 局 限 界 状態 は,
主 と して構造安 全 牲 (safety)に関 連 する性 格のもので あ る か ら建 築 物の使 用 予 定 期 間 中に おける非 超 過確率に基づ く設計規範の 構 成が用い ら れ る で あ ろう。 欧 米の統 計・
確 率 論 的考え方 を組み入れ た限 界 状 態 設 計 法の構 成におい てもこ の よ う な取り扱い が な され て い る。
一
般 的な建 築 物の使 用 限 界状 態 設 計 規 範の構 成に は,
そ の限 界 状 態の 多様 姓と超 過・
非 超 過の 2値 規 範 (binary
criteria )の適 用 性の問 題の ほ か に,
上に述べ た点 横 断 に基づ く閾値 横 断 率によ る設 計 規 範の実 用上の問題 が あ る。
例え ば,Fig.
4−
1に示 す 生 起 率 (λ),
持 続 時 間 (μd> N工 工一
Eleotronio Libraryの異な る
3
つ の荷重 過 程モデル の例を 上げよ う。 (4−
11 ) お よ び (4−12
)式に おいて eqB ・,
x一
φ一
・[
1一
去
・
・←Bs
〕1
・
……・
…・
…・
(4姻 と お く と,R
。,
x と 臥G
.
4s,
E)の期 待値に乗 じ る抵抗 係 数 お よ び荷重係数は,
下式で表さ れる。
5
・一
詳
π
・
exp 卜a・・
’
…3
…’
… R・
1
・
…・
(4−19
)7x
一
詳
π
…pl
…’
…9・・
x・
al・
W・
1
’
・
…・
・
(4−
2e) φx°
Rs、
x ≧ γx’
Wx ’
’
”…’
’
’
’
’
”「
匸
’
’
’
’
’
”・
r・
・
・
・
…
(4−
21 ) 3つ の 荷 重 過 程 隣 (t
),Ws
(t
)お よび 縣 (t
>の荷重強 さの期 待値,
変動 係 数お よ び確 率 分 布を仮に同一一
と し, かつ , λ,< λ.〈 心 とす ると, λが 最 も 小さ いW
,(t
)過 程に対するeaβs,
L,
γL は,
ほ かの2
つ の荷重の そ れ らに 比べ て小 さい埴が設 定 され る (Fig.
4−1
下図参照〉。
し か し,
荷 重WL
(t
)の持続 的作 用,
す な わ ち閾値超 過に よっ て生 じる支 障の 持 続 性に対する配 慮が実 際の設 計で は なされ るべ きで あ ろ う。
閾 値 横 断 率に基づ く設 計 規 範 の構 成には,
こ の点につ い て の配 慮 が 必要と な る。
一
般的な建築物の構造 設計に用い る荷 重の具 体 的 な 統 計 資 料に基づ く考 察は, 組 合せ荷 重の取り扱い と あわ せ WL〔U ,
λL −r
’
”一
一
L
. Ns(O ,
λS一
[
L
[
Hill
、[
L
[
L
[
ti
]
.[
Lall
.一
NE{り,
λE一
LL
,L
32q
σ O o / /ノ λ=
1,
0 // 〆 / / /,
・
0.4
/ //!・
o,
2 〆 Fig.
4−
1 12
3 βCrossing
【ate reliabilityて別稿に述べ る。
5.
結 語 本稿で は,
マ クロ時 間変動 過程と して特徴づ け ら れ る 荷重の確率過 程モ デル に つ い て閾 値 横 断 率 と そ れに 基づ くTs
期 間 最 大 値 分 布の理 論お よ びモ ンテ カル ロシ ミュ レー
ショ ンに よ る考 察 資 料につ い て述べ た。
国 際 的に も 限界状態設 計方式へ の移 行が構造 設計 法の趨勢と なっ て お り,
限 界 状 態 設計法の体系化に は建 築 物の使 用 性や安 全性につ い て の 限界状態がで き るだ け明確に, かつ 適切 に設定さ れ ること が 必要で あ る と と も に構造設計に用い る荷重の適 切な モデル化が必 要とな る。
本 稿は,
こ れ に 関連す る資 料を提 供し,
さ ら に我 国に おいて も こ の分 野 の調査・
研究が より・
一
層進められる ことを 期待し て記し たもの である。 な お,
本 報告の数値 計算に は,
名 古屋 大 学 助 手 河 野 守氏な ら びに大学院 生 間 瀬直人, 青 木和雄 両君の協力 を 得た。 参 考 文 献 1) Lin,
Y.
K.
,
森 大 吉 郎,
富 田 文 治ほ か共 訳 :構 造 動 力 学 の確 率 論 的 方 法, 培 風 館, 1972.
2) Larrabee
,
R.
0,
,
and Cornell,
C.
A.
:Combination ofVarious Load Processes
,
Jou
皿 al of the StTuctual Di−
vision
,
ASCE,
Vo亘,
103,
No.
ST
5,
January
,
1981.
3) Wen、
Y.
K.
,
:Statistical combination of ExtremeLoads
,
Journa
且of Structural Division,
ASCE,
Voi.
103,
No
.
ST S,
May,
1977.
4) Wen, Y
.
K,
:Methods forReliabi[ity of Structures underMultiple Time Varing Loads
,
Nuc且ear Engineering andDesign
,
60,
1980.
5) Tang
,
W.
H.
:Probabilistic Evaluation of LQads,
亅our.
nal of the Geotechnical Engineering Division
.
ASCE,
Vol
、
107,
No.
GT 3,
March,
1981.
6) Wen
,
Y.
K.
; Stochastic Dependencies in Load Com.
binatien
,
Preceedings of lCOSSAR’
85,
the 4 th Interna・
tional Conference on Stπuctual Safety and Reliability
,
May
,
1985.
7> Pearce
,
H.
T.
and Wen,
Y、
K.
:Stochastic Combinationof Load Effects
,
Journal
of Structural Division,
ASCE,
Vol、
llO,
No.
ST 7,
July,
1984.
8) 日本 建 築 学 会 鋼 構 造分科 会 :鋼 構造 荷重
・
耐 力 係 数 設 計 法試 案,
日本 建 築学 会,
1986.
9) 坂 本 順 ;信 頼 性 解 析 と信頼性 設計
,
建築にお け る確 率 統 計処理 (第2回セ ミナー
),
日本 建 築 学 会搆 造 委 員 会・
環 境工学 委員会・
確 率 統 計処 理合 同 委 員 会,
1986.
NII-Electronic Library Service
SYNOPSIS
UDC:624.e42.2:519.2
STOCHASTICPROCESSMODELSFORLOADS
AND
LOAD
by Dr.JVN SAKAMOro,
ofA.I.J.
COMBINATIONS'
Prof.of Nagoya
Univ.,Mernber
Stochastic
process modelsfor
loads
andload
combinations are reviewed, an emphasisbeing
placedonlevel-crossing and itsapplication toprobability
based-limit
statesdesign.
Applicability
and validity of theoreticalapproximations temean upcrossing rateand curnttlativedistribution
of the maximum value in a reference period are examined through comparingWith
Monte
Carlo
simulation results on the typicalstochastic models;i.e.,
Poisson
impulse
process.
rectangularpulse
process,Poisson
square wave processand cembined processes.Applications
tolimit
statesdesign
for
buildings
arebriefly
statedin
thispaper,/