【研 究 論 文】 UDC :624
.
042.
7;62〔}.
1 日本 建 築 学 会 構造 系 論 文 報 告 集 第 347.
号・
昭和 60 年 1月non
−white
入
力
を
受
け る
履 歴 型
多質点系
の
ラ
ン
ダ
ム
応
答解析
正 会 員 正 会 員和
木
泉
村
正
正
哲
*彦
**L
序 建 築 構造物の耐震 構 造安全性 評 価に おい て, 入力地震 動の不確 定性を考慮し た構 造 物モデル の応 答解析が 必要 であ る との視点か ら,
筆 者ら は文 献1
>,2
)に おい て マ ル コ フ連 鎖確率モ デル を用いた履歴系の ランダム応答を 求め る解析手法を展 開し,
そ れ ら を定常お よび振幅非定常な Gaussian white noise を受け る完全弾塑性型 1質
点 系の非定常ランダム応答を求め る場合に適用して き た。 こ の手 法は
,
荷 重一
変 位 (P一
δ)平 面を有 限 個の領 域に分 割する こ とに より P一
δ平 面上の応答 位 置を離 散 化 し,
これ らの離散化さ れ た応答 位 置を応答の状態 量 と み な し,
これら の状態上に応答が存在す る確率分布 (確 率 関数)を非定常に求め ること を 基本と す るe そ して,
こ の確率分布の推 移に対してマ ル コフ運鎖を仮 定し, 応 答の 状態 間の期 待推移回数に よ り規 定 さ れ る推移マ ト リックス を求め るこ とによ り, 応 答の状 態 上の確 率分布 を非定常に求め る の で あ る。 この手法に よ れば, 同時に 最大応答お よ び極値応答の確率分布,
初通過 破壊確率を も容 易に求め ること がで き る。
ところで,
履 歴系の ランダム応 答に関す る研究に おい て は,一
部の研 究 (文献3 ),
4 ))を除い て入力をwhite noise と して いるの が ほ と ん どであ る。
し か し な が ら,
一
般に建築構 造 物の受け る 地震 動が,
地 盤の存在等に よ り non−
white 性を示す こ と を考え る と,
入力と して non−
white noise を 対象と する必 要 が あ る。
ま た, 非 線形性の強い履 歴 型多質点系の ラン ダム応答を, 質点 数が
増え て も効 率よ く求め られ る こと が 望 ましい
。
本 報 告で は こ れ らの点 を鑑み,
マ ルコ フ連 鎖 確 率モデル による履 歴系の ラン ダム応 答解析 手 法 を, non−
white nQise の一
種で ある振 幅 非 定 常なfiltered
Gaussian
white noise を受け る,
Bilinear
型 多 質 点せん断 系に適用 する場 合の理 論的 拡 張お よびそ れ ら理論 解とMonte
Carlo
法に ょ る シ ミュ レー
シ ョ ン解との比 較 検 討 を行 う。 2.
解 析 手 順 * 東 北大 学 教 授・
工 博 ** 東北 大学 大学院 生 〔昭 和59年 4月 10日 原 稿 受 理 日,
昭 和59年8月6日 改 訂 原 稿 受 理 日,
討 論 期 限 昭 和60年4月末 日〉 マ ル コフ連鎖 確 率モ デル に よ り 履 歴 型質 点 系の ラ ンダ ム応 答 解 析を行う場 合の手順は, 以下の よ うにまとめる こと がで き る (文献1 ))。
} 1 > 2 ) ) 34 )5
) 6 ) 7 ) 8P 一
δ平 面を,
履 歴特 性 を考 慮 して有 限な n 個の 領域8
=iSi
,
S2,
…,
Si.
…,
S
。卜…・
…・
…・
…・
・
…
(1> に分 割 し, これ ら を応答の状態量 と み な す。
時 刻 tに応答が状態S
‘ (ま た は単に状 態i
と呼 ぶ )に存 在して いる確 率をPs
(t)と して,
ベ ク ト ルP
〔t)を次の ご と く定 義す る。 P〈t)=
IPI
(t},
P2(t),
…
,
P‘(t),
一
・
・
Pn (t}ト噛
幽
・
・
凾
(2 ) P(t}は応 答の状態S
上の確 率 分 布 (確 率 関数) その もの であ り,
応 答ベ ク トル と呼ぶこと に す る。 初 通 過 破 壊の 基 準と な る応 答の状態を定め る。
初期 状態にお け る 応 答の 状 態S
上の応 答ベ ク ト ルP
(t。)を設定す る。
応 答の状態i
か らj
へ の 単位 時間当た りの期 待 推 移 回 数μw(tm)を設 定す る。
時 点tm
か ら tm.
1 へ の 推 移マ トリッ クス πtM tm一 をμ‘,(の に基づき作 成す る。μ‘,(t”)(i,j
=
1−
n ) を む 要 素 とするマ トリッ ク スをA(tm
)と お く と πtM tm+
1 は 冶 …一
・xp(
∬
+
’
A
(・)d
・)
≒explA (tm)(tPt+1−
tm}}…
…・
…・
…
(3) と表 すこと がで き る。
時 点tm
に お ける応 答の状 態S
上の応 答ベ ク トル P(tm)(n 次の行ベ ク トル )に推 移マ ト リッ クス rrt”
’ tm+
1 (n 次の正方マ トリッ ク ス )を掛け ること に より,
P (tm
.1) を 求 める。 P(置m+’
1)=
P (tm)πtmtm+
ltt・
tt・
・
tt・
・
・
・
・
・
・
・
・
・
・
・
…
(4) 5 )〜
7 >の手 順 を入力の継 続 時 間 tm。
に し たがい,
m=
O−
(m 。−
1)に おい て繰り返す。
従っ て,
基 本 的に は履 歴 系が多 質 点 系であっ ても,
各 質 点ご とに応 答の状態間で の単 位 時 間 当た りの期 待推 移 回 数μw(t)を 適切に設 定す るこ と がで き れば,
各質点ご と に独 立に確 率 応 答P
(t
、)(i=
1−
m 。) を算 出する こと が で きる わ けで あ る。
こ こ で問 題とな る の は,
こ の手法一
39
一
を対 象 と している履 歴 系に適用す る場 合に
,
応答の状 態 間で の単 位 時 間 当た りの期 待 推 移回数 μv(t)(i,
j
;
1〜
n)す な わ ちマ ト リ ック ス A (t
)をい かに設 定す る か で ある。 こ のA
(t
)に は, 入 力 と履 歴系の特性が非定常 に取 込ま れ な け ればな ら ない。
本 論で対 象と す る,
振幅非 定 常なfiltered
Gaussian
white noise を受け るBilinear型 多 質 点せ ん断系におい て は, 次に述べ る ご と く μ”(
t
)を設 定す る も の と す る。ま ず
h
質 点の荷 重一
層 間 変 位 (P
,一
δ,)平 面をFig.
1 の よ うに n 個の領 域に分 割 し, 状 態1
お よ び状ma
n を 初 通 過 破 壊の基 準と す る。 こ の質 点の応 答の状態i
か らj
へ の単 位 時 間 当た り の期 待 推 移回数 Vtu(t)は,
既 報 (文献 1),2
))の場 合 と 同 様に, 状態i
に対 応づけら れ た線形系kLt (Fig.
2参 照 )に関す る閾値横断 問 題か ら近 似す る もの とす る。
こ の閾値は,
Bilinear系のエ ネ ル ギー一
定則に基づ き状 態i
と状me
j
の 幾 何 学 的 関 係 か ら決定さ れ る。 この ように し て与え ら れ る Pt‘Kt
) Li,
次の よ う に整 理さ れる (文 献1>)。
a.
(i=
2〜
n−
L
ノ≠i,
1,
n)kSttJ(
t
)=
髭レゴヒ〔t
)一
斥レ施(t
)・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(5 ) b.
(i=
2〜
n−
1,j
= 1,
n)hSttJ(
t
)=
なレ査(t
)・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(
6
) c.
(‘=1,
n,
」≠ の尾μ‘,〔
t
)=0・・
・
・
・
・
・
・
・
・
・
・
・
…
一・
・
・
・
…
一・
・
7r・
…
P呷
・
・
・
… r
(7
)d .
(i=
1〜
n ) kPt“〈t)=一
Σ二塩μ‘丿ω∵………・
…・
…・
…・
…・
…
(8) ’キt(
5
)〜
(8
)式に出て くる 産崛 (t
), 別氛の は,
平 均 値 0の定 常ガ ウス過 程が Fig.
2に示さ れ る よ う な線形系 kLs の変位上の閾 値 両 1,
k τ t、を,
単 位 時間当た りに正 勾 配で横 断す る期待回数の形で表せ る もの とする。
こ の 閾値hXJ、, i、XJ!は,
弾 塑 性 応 答 と線 形 弾 性 応 答のポテン シ ャ ルエ ネルギー
を等
置すう
ことによっ て,Fig.
2に示 さ れ る ご とく状 態i
と状態 ノの 幾 何 学 的 関 係か ら定めP
。
s1 認 舳 SL嵐
d6
、一
40
一
Fig
.
l Discretized Response.
State.
る。 閾値を
一
般に sXA とすれ ば, ,ul (t)は 次 式の ご と く表 現さ れ る。
… ω
一
、
黠騾
〕)・
exp{
、。蟲
ω,}
……・
(・)こ こ に
T
。(t
)は,・
入 力で あ るfiltered
Gaussian whitenoise Z(
t
)の 時 刻t
に お け る レベ ル ∬(t)に より規 定 され る量であり
,
対応す る線形系CLS
(t)(correspondinglinear
system ) (後 述)が零 初 期 条 件の も とで レベ ル1
(t
)な る定常なfiltered
Gaussian
white noise zo(τ)を受け た場 合の .
L
、の 非 定 常出 力応 答x(τ)が状態 iに平 均 的に滞 在してい る 時 間 を 意 味 す る。
(9)式に お ける 歹」匚(To
(t
)), 歹鼠τoω )は, そ れ ぞ れこ の x(τ),
叙τ)の標 準 偏 差 応 答 kσx(τ),
kα±(τ)を[0,
T。(t)]に お い て時 間 平 均し た もの であ る。T
。(t)は ・∫
(・{・)・・)・・
− 1 ……・
・
……・
…一 ・
・
(1
・) とし て,
状 態.
i
が他の状態に移 行す る場 合の基 準と なる 閾 値 Xv (Fig.
3
参照)をx(r)が超 過 する期 待 回 数 が1
P
, ’藷
’ 1∠]触冨
△ 細c’
’.
’
∠コハa 瓠=
ム ハED ’ ’ 1一
kXjl=
AC ’一
’.
k茂」2i 俎 ’ ’’
●
b e a ’ ’.
1 ’ ’
.
.
ll
撫 i/
1
ノ,
卩
卜
,
1
’i
’1
.
o ’4 CD ’ ’ c 、L
6
、 ’ ’ Sj ’ ’’
Si’
’ ’ んL
Fig
.
2 Thresho互d
Va且ues sc,匚
and fpxJ2 associa しed with ut‘i.
kP ’ ’ ’ ’ ’ ’
覧
’ B ’ C ’ ’●
’ ’ ’’
1
∠コ 3 △胚 E ’■
一
’ ’ Xy=
詔’
’ 驂 Si ’ Si+1.
1
’・
,A ED6
’と な る時 間を もっ て 定め る
。
こ こ に N ;(1(t);τ)は, x (τ)が閾値Xy を単位時 間 当た りに正 勾 配で横 断す る期 待 回 数であ る。
前 述の‘
対 応 する線 形系’
CLS
(t
)とは, 時 刻 tにお け るBihnear
型 多 質点せ ん断系NLS
(t
>の 各 質 点の 応 答 が状 態 i(h,
t)〈h
= 1〜K
,K
:質 点 数 )で あ る と き,
これ らの状 態に対 応す る線形 系inL、か ら構 成さ れ る線形 多 質点せ ん断 系を指す。
す な わ ち,CLS
(t)と は NLS (t)の 各質点の 応 答の 状 態i
(h
,
t)(h =
1〜K
)に 関して定ま る層 間 塑性 変形 σκω (Fig.
2参 照 )を介して,
,L‘をせ ん断 型に結合し た線形 系である。
ところで
,
このCLS
(t
}は,
各 、tLt 上 の原 点 (Fig.
2
のA
点 )か らの変位, 速 度を状 態 変 数とし て運動方程式 を たてる と, こ の運 動 方 程 式は,
ひ粛 )の ない場 合の線形 多質 点せ ん断 系CLSO
(t
)の層 間 変 位,
層 間速度に関す る運 動 方 程式 と等価 に な るの で,
h σx〔r),
hax(τ)はCLSO
(t
)のh
質 点の層 間 変 位,
層 間 速 度の 標 準偏差 応 答と して求め るこ と がで き る。
定 常
filtered
Gaussian
white noise 2。(τ}は,
以τ。)をパ ワ
ー
スペ ク トル レベ ルS
。(t
)の 定常Gaussian
white noise と す る と,
以 下の よ う に 以 諭 を 受ける線 形1自 由 度 系の絶 対 加 速 度 応 答の定常 部分と し て表さ れる (文 献5))。 U+2hgt・。む+ω加 =−
tLK r。)…・
・
………
(ll)・、(Te
)
一 ・(・}・ 臨 )……・
………・
…・
………
(12)x。(τ)
=
z、(τ)………・
…・
…・
・
…・
…・
…・
…・
……
(13 )r
=
To一
τOeo(τ≧
0
)・
・
・
・
・
・
・
・
・
・
・
…
曾
・
・
・
・
・
・
…
一・
・
…
(14 ) rOco≧
1
/(hg
ωg)・
・
…………
…・
………・
…・
(15> 従っ て,
z。(τ)の レ ベ ル ∬(t)を 示す指 標 として ω(τ。)の パ ワー
ス ペ ク トル レ ベ ルS
。(t)を用いるこ と が で き, 以 後S
。(t)をもっ て1
(t
)を表す も の と す る。
次 に レ ベ ル
1
(t
)=S
。(t)な る定 常filtered
Gaussian
white noise z。(τ)を受け る
CLSO
(零 初期 条 件 )のk
質 点の 層 間変 位, 層 間 速 度の標 準 偏 差 応 答κax(τ}
,
kOth(τ)を示す。一
般にkσx(τ),
磁τ)は, 次 式の よ うに 表 現す るこ と が で きる。
hσ
k
(τ)==κσis
(1− 9−
(τ))・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(
16
)、σ茎(T)三 、σ認
1− 9t
(r))・
・
一・
…・
…・
………・
一 ・
(17
) こ こ に kaxs,
、adrs は, それ ぞれk
質 点の層 間 変 位,
層 間 速 度の標 準偏 差応答の定常値である。g
.(τ),
g
£(τ)は 過 渡 応 答 性状を表す部分で,
τ の増 加 と と もに 1か ら0
へ 変 調し な が ら減 少し て い く関数で あ る。 icars,
kσdesは, 線形システム の周 波 数 領 域にお け る 入出 力 間の関係お よ びパ ワー
ス ペ ク トル密 度 関 数と二 乗平均 値の関 係か ら解 析 的に求め るこ と がで き る (文 献6})。一
方,S.
Krenk,
H.
0 ,Madsen
a皿d
P .H ,
Madsen
(文 献7))は,
定 常 過 程 を受け る零 初 期 条 件 下の線 形1自 由 度 系において, 応 答変 位の基 準 化さ れ た標準偏差応答 (ax(τ)/ a。
。)の過 渡応答性状が 入力の スペ ク トル内 容にほ と ん ど独 立で あ り,
主に系の減 衰定 数に依 存 するこ と を,
white noise 入力とfiltered
white noise 入力の場 合を例と して確認し て い る。 以上 の 点か ら本 法 で は
,
線 形 多質 点 系CLSO
の 応 答 に最 も 支 配 的 なモー
ドの 基準 座標が white noise 入力 を 受けた場 合に示す次 式 (文献8
)〉に より,
(16 ), (17 )式のg
.(τ),9t
(τ)を 近似す る もの と する。8
τ(τ)=
exp (− 2h
ωoτ)11
十2rZsint
(のdτ) + rsin (2ω。τ)卜…………・
………・
…・
(18)9s(τ)= exp (
− 2ha
,e τ)11
十2r ’sin1 〔ω ατ)−
rsin(2aJdT
)1
………・
・
…一……
(19)o
闇
二
ん/ /i
::7i
『・
・
・
・
・
・
・
…
一・
・
・
・
・
・
…
鹽
鹽
・
・
・
・
・
・
・
・
・
・
・
…
(20) こ こ に
,
h,
ω。,ωd(=
ω。Vi
=’
iii
)の モー
ダル パ ラ メー
タ は, そ れ ぞ れ減 衰 定 数, 非 減 衰 固 有 円 振 動数,
減衰固有円振 動 数であ る。
これ らの パ ラメー
タ値と して は,
Bilinear 型 多質点せ ん断 系NLS
の粘 性 減 衰 係 数を,
NLS
の質 量 お よ び第 1剛 性 分布か ら剛 性 比 例 型で与え る こと を前 提に,
線 形系CLSO
の1
次の モー
ダル パ ラメー
タ値 を 使 用す る。
こ の よ うに設 定す る理 由は以 下の ご と くであ る。 ま ず,
Bilinear型の履 歴系は剛性 低 下を考え て い な い の で,
線 形系CLSO
に おい てもその減 衰 係 数は剛 性 比例型の形 を保存して いる。
そ し て,
剛 性 比 例 型の減衰 係 数を有す る線形 多質 点 系に おいて は,
(18),
(19 )式 の &(τ},g
±(τ}の減 衰 性 状を最も規 定するh
ωn の項の 値が,
1 次モー
ドの 場 合に最 小 値 をとる の で,
こ の線 形 系の過 渡応答 性は1次モー
ドに よ りほと ん ど支 配さ れ る。 従っ て,
(18
}, (19
)式右辺のh
, ω。と して線 形 系 CLSO の 1次の モー
ダル パ ラメー
タ値 を用い る こ と が 許 容され るので あ る。
(IO) 式 に お け るN
。 ’ (1
(t
) ;τ)は,
対 応す る線 形 系CLS
(t)が零 初 期 条 件の も とで レ ベ ル1
(t
)なる定 常なfiltered Gaussian white noise xe(τ)を受ける場 合の 非 定 常 出 力応 答 頭τ)が
,
閾値Xy を単位 時 間 当た りに正 勾 配 で横 断す る 期待回 数で あ り, 次 式の ご と く表さ れ る。
N・(… )… 一 、
懇
,・xp(
x螽2h
σ1
(τ))
{
nVexp
(
一
、、i
−
・−
lll
:itlt
x ,iak
、.、)
・纛
、.、・坤
+ er ・(
ρXy2 〔1−
P2)kax (τ))
}
]
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
(21) ・一 ,覊
,.,・
一 ・
……・
……・
一 ………・
・
…2
・ K。de(τ)ニE
[x(τ)離 }]=
Dexp
(− 2
んωoτ)sin2(ωd τ)……・
……
(23 )一
一
D
= hax(τ)lath(τ)●
2ん/(1一
九2)・
・
∵・
・
・
・
・
・
・
・
・
・
・
・
…
』
・
(24
)erf(・)一
去
ガ
・・p(−
t・)・……・
…………
(・5) ここに,h,
an は (18) , (19) 式で用い ら れ るパ ラ メー
タ値を使 う。 kax (τ), tC ±(τ)は, (16
>, (17
)式で定義さ れて いる。3.
解析 結果ま ず解析モ デル につ い て説 明 する
。
入力は振幅非定 常な
filtered
Gaussian
white noise であ る が,
こ れ を便宜的に確定強 度関 数と定 常fih『red
Gaussian
white noiseの積と し て表す も の と す る。 し た がっ て
,
入力 z〈t)は 次式の よ うに表さ れ る。2(
t
);f
(t
)ζ(t)・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
…
曾
・
・
・
・
・
・
・
・
・
・
…
(
26
) こ こ に,f
(t
)は確 定強 度 関 数で次 式の よ う に与え る も のとし, その形状をFig.
4に示す。
∫(t
)=A =
0.
1609T4
= 30ζ〔t)は定
aS
filteredGaussian
white noise で あ り,
ζ(t
}のパ ワー
ス ペ ク トル形 状はFig.
5の ご と く設定し た。
Z(t}の継 続 時 間は 30秒と す る。
構造 物モデルは,
基礎固 定の Bilinear型 履 歴特性を有す る せ ん断 型の 5 質点系で,
クー
ロ ン ス ライ ダー
を用い たW .D.
Iwan
(文 献9 ))に よる表現 を 用い る とFig.
6
の よ うに表すこと ができる。 こ の履 歴 系の シ ス テ ムパ ラメー
タをTable
1 に示し,
線 形 時モー
ダル パ ラ メー
タ をTable
2に示す。
線形時の 1次 固 有 周 期が 0.
433秒,1
次の減 衰 定 数がO.
05
(減 衰 係 数は剛 性 比 例 型で与え た 。)1
で あ る。 従っ て 入力z(t)は,Fig.
5に示さ れ るごと く,
対 応する線 形系 (線形 時の履 歴 系 )に対 して相 対 的に 短周 期 成 分が 卓越す る 入力と な る。 ζ(t
)の レベ ル は, こ の 対応す る 線 形系が ζ(t
)を受け た場 合の 1質 点 目 (最下層 }の定 常r.
m.
s.
層 間 変 位 応 答が 1 と な る よ うに,
∬=S
。=
637.
35gal’
・
sec /【ad と し’
た。 〔t/7
)2 (0≦ 診≦ TI》1
(r1
く t≦ Tr)・
・
…
一・
・
一・
(27) exp (−
A (t−
T2))(Ti< t≦ Ts)B
(t−
Tnt十〇.
05
(Ts
<t
≦T
,)B =
0.
006 τ1=
10T2=
15 コ「li
=
25Table l Paraineters of a 5
−
DOF BilinearSystem
(shear type)
.
1m1 区1iK211KuC 且 d1 工st1.
03000.
00.
120.
6651.
0 2開d1.
o2700,
00.
1
P18 .
5991 』 3rdLO2400,
00.
116.
5321.
0 4Lh1,
0210D.
00.
114.
4661.
0 5th1.
01800,
00.
1 ユ2.
3991.
0(
ご』
1.
o.
0 5 10 15 20 25 TIME (sec )Fig
.
4 Deterministic IntensityFunction
.
10 8 コ
)
り
哨 o6 4 。 の、
(
3)
の 日 5 勘 90 Φ 幽 の 2 臼 Φ き O門
0 024681012 (tU!r) (1!sec )Fig
.
5 Power Spectral Density Function Qf ζ(t).
3G マ ルコフ連鎖確率モ デ ル によ り 導 か れ る解は
,
応答ベ ク トルP
(tt}(i=1
−
m 。)で あ る。
最大 応 答お よび極値 応 答の確 率 分布は,P
(t
‘)が求ま れ ば容易に算 出さ れ,
その 精 度 はMonte
Carlo
法 に よ る シ ミュ レー
シ ョ ン 解 と の 比 較 に お い て P (tnの精 度と同程 度で ある こと が確 認され て い る (文 献1
),2
))。 し た がっ て, こ こ では応 答ベ ク トルP
(tE
)につ い て の了able
2
Modal Parametersof a Linear S¥stem Correspond茸ng to a 5
−
DOF Bilinear‘
System,
mi f iaS・ ・f ⊥th Gt・rey.
(t 貿・
・ec’
lcm) K11;
ft「
st stlfrness of lth Btorey (tftcm ) K21:
seCond stlffness of jth stOrey (tf/cm):
:
:
;
1
:
!
:
;
:
、 d器
蠶
:
1
:
∵
flllh。
器
y 、謡
sec ’cm ’:
1
iTi 〔sec ) 五 1s ℃ o.
4330.
050 2ndO.
158O.
137 3rdO、
1010.
214 4th0.
0790,
275 5th0,
066o,
3z6 ith nat 口ral perSed圸 m◎dal dt叩 ⊥ng factor
K2i
’
冒
’
K
}
}
5ー
血 「.
ー
一
曽
2
1
.
F一
}
.
畿
Ci
r_
一
_
Kli− K2i
Fig
.
6 She6r Type 5−
DOFBilinear System
.
一 42 一
,
一
創.
一
.
ト一
.
吩口
.
0 苛 囚11NM.
攣
凵一
冨 QO 」 」韓
.
塗
O.
a 皐一
掣
.
霪
d 諸 コ 凵塁
瞿,
一
創,
g一
卩,
卜冖
,
の一
.
円一
.
一
一
.
の.
卜.
鴫.
O,
一
.
ONN11N 乂.
匿 卜乂 肛 00 一 」 Q 酌 創
,
竃
工 り.
繋尸
03 “ の マ,
口 巳 0 ←.
雨.
田゜
」二
巴d
略 」 漣讃
凵 肛・
卜・
ゆ・
° 二一
〇 〇 鴎 昨 望,
§
期一
乂 配 QQJL ← oor 偲,
靉
工 ゆ.
国 搴一
〇⊃掣
.
竃
卜 R 『 9 的 「 R 笛 ロ ロ.
国 N 国 从1「1日刪 d ゆロ
コ
寸.
ロ コ
oecロ
一
ロ 口 N 国 国 笆 e 《⊥nl臼ザ田 畄 頃
’ リ
マ の e占ヒ の
一
り
口 嘘コ 8 国 翰 NN 从 工「1臓 d唱
」 凵;
1
.
卜 ゆ m.
(
ho2
隅 唱 h6 り)
冠 > q」
一 口 o ロ o 国 } o口 o
墹
ぢ 口 5 隅 あ u = 三 幡 』 o」
氏,
9 」望
凵 」 凵羃
毘.
「
U.
の尸
.
ト
一
.
ゆ「
.
O一
,
「
尸
.
叨.
卜.
ゆ.
の.
「
誇
・晶
エ の・
国 n 从 マ.
ロ O,
国.
O τ N 田 夐.
蕁
一
¥ ¢ 00 」 恥 Q 馬 尉,
塞
工 の.
孚}
oo 創 の 寸,
国凸
.
(
ho2 臨 ℃ 自8
{
〉 ω 口 Φ g ワ ロ o ロ ω 国 } o切
‘ o嘱
ぢ 偶 5 」 か唱
霈 』 幡 阻 9 畠 o ◎.
9 」、
口旁
凵羃
瞿.
一
N,
の一
.
ト
一
.
の一
,
〇一
,
一
一
−
の.
卜.
り.
O.
一
謡
轟
エ O.
qn 杜.
寸.
国 O.
国.
ONN11N ¥.
圏 卜 窄一
図 さ D 」 匸 Ω 面 N.
塗
噂.
口皐
一
〇⊃掣
.
竃
ト.
〔
国
o 』 隅 冨 こ て 〉 口 。コ ロ o 島 切 Φ 国}
o。
。 口 o コ 9 口 5 餓 』 三刷
』 鳴 』 oh 山 卜.
O 匚皀
」蹇
昌 匡.
N.
田.
俺,
雪,
2.
二.
の.
卜.
り.
o.
ロ・
ロ鈩
.
回 O.
国 n 瓧 寸・
国゜
8n 閏 望゜
§
⊆§
ビ ト 野 鰻埀
の’
§ 萬窺
塞
ト ・,
国『
一
cu.
の一
,
」一
.
〇一
−
□ 己 」 凵豢
瞿PO
一
『
一
”
,
切,
卜,
O,
09.
一
.
0 茸 N 既 夐.
麈 ▽ 畔一
乂.
一
W.
oo冖
『
卜【
.
の一
巳 D」
匚 ( 匣 の 四.
竃
亠 O,
蜃 峯一
〇孥
,
国 皿 0 卜 昌 己 」護
毘.
『⊃「
,
一
冖
.
切.
卜.
頃.
O辱
一
一
〇 」 四 脾 N ヱ.
瑠 卜 個 陰一
乂 配 QO 」 」 自 7 潤 個.
竃
亠 切.
§一
〇
9 芻 寸,
塞
卜 」 凵 己 」 凵 oo罍
り
凵 匡.
一
周,
oウ一
.
ト
一
−
の一
,
〇一
、
冖
一
.
口り.
卜.
恥.
O、
一
の マ oov−
N N 匸 ≧ N 国 国 国 人1Σ「19∀eOUc] o’ ロ ロ ’
寸、
e 創一
国 N 口 e 国 国 N 从lr工日VEDecs ゆロ サ ロ ロ
寸 o く丶3−
e.
NeN 巴 N 口 A⊥1「1日V臼Dごd,
〇 一 N1 − N乂
.
園冖
α 聞一
y き ○ 」 匡 工ト
マ 創.
墨
噂.
睾一
u } eJ 等.
0 ←.
」 畫[
の 0 ヨ匡
く り]
bZOE「
]
」 O 乙 く冖
匡 卜 》匡
O ] = ト”
] 一 〕匡
同
Qり
凵
〇 四 O く ゜ ぐ Oo ( ロ
.
国一
.
O 的、
。 脳 o,
o 田 川 寸.
国 杜 の.
ロ O.
騒一
釦 雪 卜冖
面 己 」 〕窰
毘 の一
.
O7
一
冖
鹽
卜 ゆ,
卩
つ一
隅 占一
寸)
一
〉 口 O 翰 口 O ユ ひ》 国}
Oの
口 O = リ ロ , 叫 あβ
呂吶
』 鳴 』 Oh 山d
> 凵 」 凵詈
巳
の
凵 α冖
N.
自D一
,
ト
一
,
里FO
一
.
;、
−
ト めー
eっ.
一
壕
ー
線
ロ・
凸 ト9
.
O 匚 T一
1 」コ
≡の
oo ぐ ○ 0」
伍
可
U 国 ト ZO 冨一
〕
」 OZ く冖
匡
卜 》 匡 0 凵工
ト”
]
一 〇 配[
りり
凵 の 旨 Φ ぐQ.
OHN11 翌,
国 ロ一
窄}
ど營
ピ エ ト マ e9.
茎
頃,
国 灘冖
αり N 零.
.
ロ.
国,
O.
国 螢、
ロ.
O 槲 瓧 曷(
Ho 。 O 円 N11 鬯.
箇 国一
創 11一
乂 匡 ∩ σ 」 L 工 ← マ N.
国 咀 0 工 り.
国 案一
q ⊃ 芻 Ψ.
竃
ト 」 畫一
しっ O 」配
く り]
匝 ZO = 差 O ]エ
ト 凵」
OZ く[
配 ト 凵 d 匡 冒 団の
9
o 〈 Oq O国
,
〇一
.
O.
騒 跏 n.
邑 刪 n マ.
国 杜 ゆ.
国 O.
国 里 = 巴 1 口 〉 国 」 凵 の ZO 圧 し石 凵 配.
eり一
.
;.
,
卜 の n.
O°
〇
一
11 望.
§
了一
乂髫
ピ ≡ b 尉,
ロ亞
り.
国 躍一
め 〜 孚.
孚 0 ト.
一
碧圏
し
っ 〇 一 匡 く U 凵 ト = Oε
・
凵 一 〇 Z く[
匡 ト》
匡 O ] = ト.
]一
〕
匡[
〕 冨 の ON o oo O,
国 餌,
O 『畆
エ O,
011 杜 g2.
q 国.
「
濁 叨一
,
ト一
,
切卩
ー
口 己 」 凵豢
凵 匡,
O卩
.
一
一
『
.
卜.
ゆ め.
卜尸
.
の一
.
O一
.
一
一
,
O.
卜.
の,
O.
「
,
−
」 コ 冨[
の
〇一
匡 て り ] ト ZO =゜
〉 配 0凵
弄
ト
り
芻、
Ooり
一
1− NM,
匿 四一
幽「
Z = 00 」 匸 = 巳 創.
墨
O,
国 畠一
謬
,
孚 卜詩
・,
晶
エ o.
ロ ー− 紅 O.
9 国.
一
(
】
8
.
匡 砦 ε,
一
N.
2 届 さ 日 Φ 昌a
の 。 国 で の8
−
口 己 」護
凵 配 臣一
西 里 = 。 塁 」 ゐ噌
「
三 霊 o 」 =.
O 匚.
Ooo蔓
,
§
朝
斐
髫
ピ百
α,
守 望 O,
愛一
u }掣
.
甲 凸 ←り
”の
OH.
一 ∋ … し。〇
一 匡 くり
国 ト ZO 竃 差 O 凵 譯 凵 」 OZq[
匡 ト 国」
星 ご o 」一
.
口 己 」 凵篝
凵 匡藍
「
駐
国.
O 創.
国轟
エ o.
国 n 瓧 O.
国 ロ.
一
,
Oe み本 手法 に よ る理 論 解 とMonte
Carlo
法に よ る シミュ レー
ショ ン解 との比較 検討を行 う。 理 論 解とシ ミュ レr ショ ン解と の比 較は,
各応 答の状 態に対 応づ け ら れ た ナン バー
(Si
の の を 用い て,
応答 レ ベ ルX
‘を 次 式に よ り定義 し,Xt
に対 す る確 率分 布の比較と して行う。X‘
=
li
− iml
十1・
−s・
・
…
(28)
i
:応 答の状 態の ナンバー
(1≦i
≦n ) 畆 :P 一
δ平 面の 原 点 を 有 する応 答の状 態のナン バー
Xs
は,
応 答の 状 態S
‘のP 一
δ平 面の原 点か らの離れ具 合を示す確 率 変数であ り,
応 答 レ ベ ルの大 き さ を 表す指標と な る。 各 質点と も応 答の状態の 総数 n を 41と し,
応 答 の 状 態 の 幅d
(Fig.
1
参照)を各質 点の降伏 変位 値の1
/5
と し た。 こ の値の n,d
に 対して,
初 通過破 壊の 基準と な る 応 答の状 態Si,
S
πは, 原点か ら塑性 率で5,
44
の位置に あ る。
ま た,
(3
) 式右辺の (tm+「 tm)は,
すべ ての m (m=
O〜
29 >に わ たっ て 1 秒と し た。
Figs.
7−
U
に,
各 質 点につ い て 応 答レ ベ ルX
‘上の確 率 分 布 (確 率 関 数)の 理論 解 をMonteCarlo
法 に よ る シミュ レー
シ ョ ン解と ともに 示す。 各 質 点とも10,
15,
20秒の3
時点につ い て示 して あ り,
丸が理 論 解,
三角が シミュ レー
シ ョ ン解で あ る。 質 点の ナン バー
は,
下から 1,
2,……,5
であ る (Fig.6
参照 )。
Monte
Carlo
法に基づ く シミ ュ レー
シ ョ ン解 は, 模 擬 人 力 波 100波に よ る応 答 波 群の統 計 処 理 を各 時 点に お け る集 合 平 均に よ り算 出し た もの で あ り,
応 答 計算は線型 加 速度法に 拠っ た。Figs.
7〜
11の結 果 をみ ると, 1〜
4質点 につ い て は理論 解は シ ミュ レー
ショ ン解に対し て全 体 的にや や 安 全 側の評 価になっ てい る が, 工学 的 精 度に お い て はよく対 応し て い るとい え る。 特に入 力の減衰 部分にある
20
秒 時の確 率 分 布が, 主 要 動の最 終 時 刻で あ る15
秒時の確 率 分 布に対 して, P一
δ平面の原 点 方 向に指 向して いる 傾向を シ ミュ レー
ショ ン解は示して い るが,
これ を 理論 解が再 現し て い る点は,
理 論 解の妥 当性を考え る上で重要で ある。
最 上にある 5質 点 目につ い ては,
理論 解は シ ミュ レー
ショ ン解に対 し過 大に安 全 側の評 価を し て お り, 両者は よい 対 応を示 し て い な い。
本手法におO.
ては,
対 応す る線 形 系の確 率 応 答の情 報とエ ネルギー一
定 則か ら, 各 層の推 移マ トリックス を近 似し てい る が, 最 上層 に お け る理 論 解の過 大 評 価は,
最 上 層へ のエ ネル ギー
分 配 が 過 大で あっ た こと を意 味して いる。 この点につ い ては, 履 歴 型 多 質 点 系 を規 定 する システ ムパ ラ メー
タ との関連 も含め て, 今 後の検 討 課 題と し たい。 し か しな が ら,
この点 を 除け ば, 理 論 解は全 体と して シミュ レー
ショ ン解 を よ く 再 現 し て いる といえ る。 4.
結 語 本 報 告で は,
従 来 white noise 入 力 , 履 歴 型 1質 点 系 に対し適 用さ れて いたマル コ フ連鎖確率モ デル によ る履 歴 系の ラ ンダム応答解析手法 を,
non・
white noise 入 力,
履 歴 型多質 点 系へ 適 用する場 合の理 論 的 拡 張を行い
,
基礎 固 定せん 断 型の
Bilinear
型 5質 点 系が,
振 幅 非 定 常な
filtered
Gaussian
white noise (Bilinear系に対応す る線形系の 1次 固 有 周 期に対し, 相 対 的に短 周 期 成 分が 卓 越 ) を 受 け る場 合 につ い て,
その理 論 解 をMonteCarlo
法に基づ くシ ミュ レー
シ ョ ン解と比 較 検 討した。 こ の手法に よる ラ ンダム応 答 理 論 解は,
各 質 点の P一
δ 平 面に設 定され た応 答の状 態S
上の非 定 常 確 率 関 数(応 答ベ ク トル )によ り表され るが,
シ ミュ レー
ショ ン解と の 比 較に お い て は, これ と同 等な応 答レ ベ ルX
上の非 定 常 確 率 関 数により行っ た。 その結 果,
最上にあ る質点 の応 答が安 全 側に過 大 評 価す る点を除い て,
各 質点と も 理 論 解は シ ミュ レー
シ ョ ン解と お お む ね よい対 応を示し た。
特に, 入 力の減 衰 部 分に あた る時 点の確 率 分 布が, 主 要 動の最 終 時 点の確 率 分 布に対し,
P一
δ平 面の原 点 方 向に指 向する傾 向 を理 論 解が再 現 し て い る点は興 味 深 い。
今回扱っ た入力は, 履 歴 型 多 質 点 系に対応する線形 系 の1
次固有周期に対して相 対的に短 周期成分が卓越す るnon
・
white noise で あっ た が,
長 周期成分が卓越す る場合につ い て も
,
目下 検 討 中で あり,
別 途 報 告 する予 定で ある。 ま た,
周 波 数 非 定 常 入 力に対して も今 後 本 法を適 用する こと を考えて い き た い。 参 考 文献 1) 和泉正哲,
渡 辺孝英, 勝 倉 裕, 木 村 正 彦 :マ ル コ フ連 鎖 確 率モ デル によ る履 歴 構 造物のランダム 応答 解 析,
日 本建築 学会 論 文 報 告 集,
第 324号,
昭和58年 2月,
pp.
18−
27.
2)和 泉 正 哲,
木村 正 彦 :振 幅 非 定 常入力を受け る履歴系の ランダム応 答 解 析,
日本建築 学会論文報告集,
第335号, 昭和 59年1月,
pp.
15−
21.
3) 浅野幸一
郎:弾塑 性 形 履 歴 特 性 を有 する質 点 構 造 物の確 率 統 計 的 地 震 応 答 そ の 3.
擬 定 常Nen・
white Noise 入 力を受け る 1自 由 度 系,
日本建 築学会 論 文報 告集,
第257号
,
昭和52年7月,
pp.
67−
73.
4)M
.
Ozaki, S.
Yasuda and M.
Ishii:APractical Method
for the Earthquake Response Analysis of Multi
−
story Structures with Degrading Stiffness Systems by RandomVibration Theory
,
Trans.
of A.
1.
J.
,
No,
336,
Feb.
,
1984
,
pp.
21−
33.
5}H
.
Tajimi;AStandard Method of Determining theMaximum Response of a Building Structure Duling an
Earthquake
,
PToc.
ef 2 nd WCEE,
Tokyo,
Vol.
2,
July,
1960
,
pp.
781−
798.
6) A
.
D.
Kiureghian ;Structural Response to StationaTyExcitation
,
Joumal
of the Engineering Mechanics Di・
vision
,
ASCE,
Vol.
106,
No.
EM
6,
Dec.
,
1980,
pp
.
1195−
1213.
7) S
.
Krenk,
H.
0.
Madsen and P.
H,
Madsen :Stationaryand Transient Response Envelopes
,
Journa亘of Engineer−
ing Mechanics
,
Vol.
109,
No.
1,
Feb.
.
1983,
pp.
263−
278
.
8) T
.
K.
Caughey alldH.
J
.
Stumpf:Transient Response ofaDynamic Systern Under Random Excitation
,
Journal
of Applled Mechanics
,
Dec.
,
1961,
pp.
563−
566.
g) W
.
D.
Iwan : A Distributed Element Model forHysteresis and Its Steady
.
State Dアna皿ic Response,
Journal
of Applied Mechanics,
Dec.
,
1966,
Vol.
33,
SYNOPSIS
'
UDC:624.042.7:620,1
RANDOM
RESPONSE
ANALYSIS
OF
HYSTERETIC
MDOF
'
SYSTEM
TO
NON-WHITE
EXCITATION
'
byDr. MASANbRI IZUMI,Pr6fessor of Toheku UniveTsity
and MASAHIKO KIMURA, GraduateStudentof Tohoku
,,
'
University, Members of A.I.J.
'
Itisnecessary to estimate a
probabilistic
'seismic
safety ofbuilding
structures to a non-white excitation,because
most of expected earthquake motions will have non-white spectral properties forthe existence of thegrotind.
In
thispapera random response anqlysis methodby
Markov
chain stochastic modelhas
been
extended tonon-white excitation inputand hystereticMDOF system case and thisextended stochastic model has been applied to
a shear type s-DOF
bilinear
system tohmplitude-nenstationary
non-white noise whose predominant period is shorter thanthatof alinear
system egrresponding tothehysteretic
systerp.This theoreticalmethod givesi a nonstationary probability function
P(S)
of discretizedresponse stateS
inforce-deflection
plane.The
comparisonbetween
theoreticalresults'and' simulation onesbased
onMonte
Carlo
method
has
been
madeby
a nonstationaryprobability
functionP(X) of response levelX, which isequivalent toP(S).・The
theoreticalresults have showed relatively good agreements with,simulation ones exceptfpr
'
'
overestimation of thetop
floor
(s
thfloor).
・
'
t
t
'