日本建築学会大会学術講演梗概集	 （関東）	 2015 年	 9	 月	 21290

摩擦機構を有する免震構造の実用的な解析法とモデル化に関する一考察

正会員
同
同
同－その 1 単体モデルによるパラメタリック解析－

免震構造
HHT-α法摩擦機構モデル化

地震応答解析

1. はじめに減衰機能を併せ持つ免震支承には様々な種類があり，摩擦や滑り機構（以下，摩擦機構という）で減衰性能を発揮するタイプもそれに含まれる。免震構造の応答評価には主として地震応答解析が用いられるが，その時間積分にはNewmark-法が多用されていると思われる。免震層に摩擦機構を含むようなモデルを解析すると，特に応答加速度において摩擦固着時に関わらず入力加速度以上の，また固着時と滑動時との移行の瞬間に高いピークの高振動モード（以下，不要な高振動モードという）が現れることがある。
これは積分時間間隔を極めて細かく設定すれば回避できることが多いが，逆に解析効率を低下させることになる。
一方，解析精度を保ちつつ数値減衰効果で不要な高振動モードが抑制できる解法として HHT-法 1)が知られている。
本報および次報で，摩擦機構を有する免震構造の実用的な解析をめざし，HHT-法や粘性モデルを用いることによる不要な高振動モードの抑制効果を確認する。本報は固着時（剛体）と摩擦機構単体モデルのパラメタリック解析を，次報で建物モデルによる比較検討を行う。
2. HHT-α法
HHT-法は次式で表される 1)。
M an+1+(1+) C vn+1-C vn+(1+) K dn+1-K dn=-M aG,n+1(1a)vn+1=vn+t(1-)an+tan+1(1b)dn+1=dn+tvn+t2(1/2-)an+t2an+1(1c)=1/2- ，=1/4(1-)2 ，=0～-1/3(1d,e,f)
ここに，M,C,K=質量,減衰係数,剛性，a,v,d=加速度,速度,変位，aG=入力加速度，t=時間間隔である。本検討では=-1/3 とする。なお，式(1)より HHT-法は=0 のとき Newmark-法（=1/4）に一致する。
3. 摩擦固着時（剛体モデル）の応答
まず摩擦固着時（剛体）の応答を確認する。固着時の応答
加速度は理想的には入力加速度に一致する。検討モデルは
質量 M=1 とし，固有周期 T=0.01s となる剛性 K=2，または粘性 C=を与える（=2/T，v=u として C=K/）。入力波はJMA KOBE NS，内部粘性減衰（C 以外）は非減衰とした。
図 2.1 は剛性 K=2 を与え，Newmark-法で解析した応答加速度の時刻歴で，積分時間間隔をt=0.002，0.001s とした結果である。t=0.002s では 10s 以後に不要な高振動モードが現れているが，t=0.001s と細かくすればそれが抑制できた。
図 2.2 は同じく K=2 を与え，HHT-法で解析した結果である。HHT-法ではt=0.01s と，粗い積分時間間隔でも不要な高振動モードが抑制できた。

山崎 久雄*1○ 大原 和之*2
中郡 一雄*3
中村 秀司*4

一方，図 2.3 は粘性 C=を与え，Newmark-法で解析した結果である。この場合もt=0.01s で抑制できた。
固着時の特性を剛性で与える場合，Newmark-法ではtを十分細かくする必要があるが，HHT-法では比較的粗いtでも理想的な応答として求まった。また粘性で与える場合は，Newmark-法でも比較的粗いt で良好な結果が得られた。
4. 摩擦機構単体モデルの応答
次に前述の摩擦固着（剛体）モデルに対して摩擦係数=0.2を与えた単純な Bi-Linear モデルを考える。ここでも固着時の特性に固有周期 T0=0.01s となる剛性 K0=2 または粘性 C0=を与える（添え字 0 は固着時の特性を表す）。これ以外の解析条件は前述と同じである。
K0=2 および C0=を与えたモデル模式図を図 3.1 に示す。
図 3.2 は固着時の特性を K0=02 で与え，Newmark-法で解析した応答加速度で，t=0.005，0.002，0.001s とした結果である。t=0.002s は 0.005s に対し改善されているが，所々に不要な高振動モードが現れている。t=0.001s では良好となる。
1000 (cm/s2)

t=0.002s

Input
Response

0010203040Time (s)-10001000 (cm/s2)

t=0.001s

Input
Response

0010203040Time (s)-1000

図 2.1K=2 による剛性モデル・Newmark-法

1000 (cm/s2)

t=0.01s

Input
Response

0010203040Time (s)-1000

K=2 による剛性モデル・HHT-法

図 2.21000 (cm/s2)

t=0.01s

Input
Response

0010203040Time (s)

-1000

図 2.3

A Sutdy on Practical Analysis and Modeling of Seismic Isolation with Friction System(Part 1) Parameter Analysis by Simple Model― 579 ―

C=による粘性モデル・Newmark-法
YAMAZAKI Hisao，OHARA Kazuyuki
NAKAGOORI Kazuo，NAKAMURA Hideji図 3.3 は同じく K0=02 を与え，HHT-法でt=0.01，0.005sとした結果である。t=0.01s では不要な高振動モードがやや残っているがt=0.005s では改善される。
図 3.4 は固着時の特性を C0=0 で与え，Newmark-法でt=0.005，0.002s とした結果である。t=0.005s では不要な高振動モードが完全に除去されないがt=0.002s で良好となる。
図 3.5 は同じく C0=0 で与え，HHT-法でt=0.01，0.005sとした結果である。t=0.01s でも良好であるが，t=0.005s ではほぼ完全に不要な高振動モードが除去できている。
Q/W
K0=02を与えた
剛性モデル

Q/Wu

K0=2u1000 (cm/s2)

01000

摩擦モデルの模式図（K0,C0 は M=1 の場合）

(cm/s2)t=0.005s

101000 (cm/s2)

Input
Response2030

t=0.002s

40Time (s)103010

図 3.43040
Time (s)t=0.01s

Input
Response

040Time (s)-1000

20C0=0 による粘性モデル・Newmark-法

1000 (cm/s2)

Input
Response
20Input
Response-10000040Time (s)0-10001000 (cm/s2)30

t=0.002s

01020-100000Input
Response0u

C0=0v

図 3.1

t=0.005s

Q/W

Q/W
C0=0を与えた
粘性モデル

一方，すべてのケースの荷重-変位関係についてはいずれもほぼ同等の結果を得た。図 3.6 に，不要な高振動モードが最も現れた図 3.2 のt=0.005s と良好な結果のひとつである図 3.5のt=0.005s での結果を示す。
5. まとめ摩擦機構の固着時の特性を大きな剛性または粘性で与えた場合の積分時間間隔t と積分法（Newmark-法と HHT-法）に着目した単体モデルのパラメタリック解析を行った。固着時を剛性で模擬した場合，Newmark-法に対し HHT-法による不要な高振動モードの抑制効果が確認できた。また粘性で模擬した場合でも同様の効果はあるものの Newmark-法でも比較的粗いt で良好な結果が得られることを確認した。

010203040Time (s)-10001000 (cm/s2)

t=0.001s

1000 (cm/s2)

Input
Response0t=0.005s

Input
Response

0010203040Time (s)-1000010203040Time (s)

-1000

図 3.2

K0=02 による剛性モデル・Newmark-法

1000 (cm/s2)

t=0.01s

図 3.5Input
Response

00102030C0=0 による粘性モデル・HHT-法

K0=02，Newmark-，t=0.005s

C0=0，HHT-，t=0.005s

0.30.3

Q/W-1000

01000 (cm/s2)

t=0.005s

-20

Input
Response

10203040Time (s)

-1000

図 3.3

K0=02 による剛性モデル・HHT-法*1 ユニオンシステム㈱*2 (株)建築構造研究所*3 (有) NCU 一級建築士事務所*4 新日鉄住金エンジニアリング(株)002040(cm)図 3.6

荷重-変位関係

-0.3

00Q/W

40Time (s)-2002040(cm)

-0.3

参考文献
1) Hans M.Hilber ， Thomas J.R.Hughes and Robert L.Taylor ： Improved
Numerical Dissipation for Time Integration Alogorithms in StructuralDynamics，Earthquake Engineering and Structural Dynamics，Vol.5，1977*1 UNION SYSTEM INC*2 Building Structure Institute*3 NCU Structural Design,Inc*4 Nippon Steel ＆ Sumikin Engineering Co.,Ltd― 580 ―