JP6358404B1 - 非線形多自由度系の動的応答解析法 - Google Patents
非線形多自由度系の動的応答解析法 Download PDFInfo
- Publication number
- JP6358404B1 JP6358404B1 JP2017551722A JP2017551722A JP6358404B1 JP 6358404 B1 JP6358404 B1 JP 6358404B1 JP 2017551722 A JP2017551722 A JP 2017551722A JP 2017551722 A JP2017551722 A JP 2017551722A JP 6358404 B1 JP6358404 B1 JP 6358404B1
- Authority
- JP
- Japan
- Prior art keywords
- degree
- response
- freedom
- mode
- linear
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
塑性領域におけるモード間の相互作用を無視し、モードの重ね合わせにモード重ね合わせ則を仮定した上で、弾性領域の固有振動モードにより非線形多自由度系をモード展開して、非線形一自由度系の時刻歴応答解析により作成した一定塑性率応答スペクトルから得たモード塑性率を用いることにより、各モードの最大変形を把握する。
Description
によって算出されるr次のモード降伏加速度が、r次の固有周期をTr、減衰比をhrとしたときの弾性応答スペクトルae(Tr,hr)以上であることを、多自由度系の線形応答解析を多自由度系の非線形応答解析の簡易法として用いることができる条件とすることが好ましい。
によって定義される想定荷重係数ε(ar y)が、閾値α以上であることを、多自由度系の線形応答解析を多自由度系の非線形応答解析の簡易法として用いることができる条件とすることが好ましい。
耐震設計において、機器・配管系は降伏後も剛性が低下しない仮想的な系、即ち、線形系として扱いその応答を線形応答解析により求めることが一般的である。線形多自由度系の応答解析法には、時刻歴応答解析法と応答スペクトル解析法がある。前者は系の最大応答の厳密解が得られるが、解析パラメータの変化が解析結果に及ぼす影響程度の予見性がない。後者は各モードが最大応答となる同時性について仮定し、モード重ね合わせ則を用いることにより系の最大応答が得られるため、厳密解は得られないが、解析に用いる加速度応答スペクトルから各モードの最大応答が把握でき、解析パラメータの変化が解析結果に及ぼす影響程度の予見性に優れている。また、加速度応答スペクトルへの保守性の付与も容易であり、機器・配管系の耐震設計においては、応答スペクトル解析法が広く用いられている。
(1)一自由度系の弾塑性応答は、剛性の低下や履歴減衰の影響を含め応答解析により厳密に計算されている。
(2)一自由度系の弾性応答の固有周期を横軸としており、弾塑性応答の程度に係わらず着目する系の弾性応答の固有周期により、応答スペクトル上で着目する周期を特定できる。
(3)一自由度系の降伏耐力(又はこれを系の質量で除した降伏加速度)を縦軸としており、これを与えることにより、弾性応答から弾塑性応答まで連続的に系の最大応答を塑性率という形で与える。
(1)弾性応答の固有振動モードでモード展開すると、弾塑性応答により剛性が低下した際にはモード間の相互作用を生じるが、これを無視する。
(2)一定塑性率応答スペクトルからは各モードの最大応答のみ得られるため、線形多自由度系の応答スペクトル解析と同様に、モードの重ね合わせにモード重ね合わせ則を仮定する。
2.1 一自由度系の運動方程式
荷重・変形関係が二直線骨格曲線であり移動硬化則に従う復元力特性をもつ一自由度系を考える。応答が弾性領域にある系の運動方程式は次式となる。
一定塑性率応答スペクトルは、外力F(t)が作用する減衰比h、剛性比γの一自由度系の降伏加速度ayを、特定の塑性率μに対して固有周期Tnの関数ay(Tn,h,γ,μ)として図示するものである。この図は、固有周期Tn毎に、外力F(t)が作用する減衰比h、剛性比γの一自由度系に対して降伏加速度を変化させ時刻歴応答解析を行い、その応答が特定の塑性率μとなる系の降伏加速度ayを収束計算により求め、作成する(非特許文献3)。島根原子力発電所2号機原子炉建物で観測した鳥取県西部地震の床応答を図3に、この一定塑性率応答スペクトルを図4に示す。なお、式(5)、(6)より、系の降伏加速度ayと入力加速度
弾性応答スペクトル解析法と弾塑性応答スペクトル解析法による一自由度系の最大変形の求め方を図5A及び図5Bを用いて示す。
xmax=μxy=(ae/ay)xy=ae/(ωn 2)
となる。ここで、式(11)よりμ=ae/ayである。
一定塑性率応答スペクトルは、一自由度系の弾性応答の固有周期と降伏加速度を与えることにより、弾性応答から弾塑性応答まで連続的に系の最大応答を塑性率という形で与える。この優れた特性を有する一定塑性率応答スペクトルを用いる非線形多自由度系の応答スペクトル法を考える。
N個の質点を線形ばねで、互いに又は床と結合した線形N自由度系を考える。質点iの質量をmi、質点i−j間のばねi−jの剛性をkijとすると、この系の運動方程式は次式となる(i,j=1,2・・・N)。ただし、質点i−i間のばねi−iは質点iと床の間のばねを表す。
前節の線形N自由度系のばねを非線形ばねとした非線形N自由度系を考える。非線形ばねは、荷重・変形関係が二直線骨格曲線であり移動硬化則に従う復元力特性を持つとする。一例として非線形N自由度系のモデル図を図6に示す。応答が弾性領域である系の運動方程式は、式(12)、(16)及び(17)に等しい。次に、非線形N自由度系の応答が弾性領域から塑性領域に入り、ばねi1−j1からばねin−jnでまで順次降伏する運動を考える。この運動は、ばねi−jの降伏変形をxij y、降伏したばねの二次剛性をk’ijとすると、以下のように記述することができる。
前節で得たモード座標{q}の塑性領域の運動方程式(22)は、左辺第3項の総和のうちs≠rであるN−1個の項が、r次モードと他のモードの相互作用を表す。また、降伏したばねi1−j1からばねin−jnは、r次モードの応答のみにより降伏したものではなく、各モードの応答を重ね合わせた全応答により、各ばねがその降伏変位となり降伏したものである。
地震による床応答に対して、一自由度の非線形系より対応する線形系の塑性率が大きくなるか判定する条件が確認されている(非特許文献10)。この一自由度系の判定条件を用い、地震による床応答に対して、多自由度の非線形系より対応する線形系の塑性率が大きくなるか判定する条件を検討する。
地震による床応答により弾塑性応答する、即ち、降伏加速度ayが
一自由度の非線形系と対応する線形系の塑性率の関係を、多自由度の非線形系と対応する線形系のモード塑性率の関係に当てはめる。式(37)から(42)より、モード降伏加速度ar y及び固有周期Trが次の判定条件
機器・配管系の設計に適用する許容塑性率として、重量機器に対し1.2〜2、配管系に対し1.5〜3、鋼材の引張りと曲げに対して2.5〜10等が提案されている(非特許文献8)。この程度の塑性率を生じる系の動的応答解析に、弾塑性応答スペクトル解析法を適用することを想定し、その有効性を確認するため、床応答に対する非線形二自由度系の塑性率を弾塑性応答スペクトル解析法及び時刻歴応答解析法により計算する。また、対応する線形二自由度系の塑性率を、拡幅率Δ=0.1の弾性床応答スペクトルを用いて弾性応答スペクトル解析法により計算する。
解析する非線形二自由度系は、2個の質点を非線形ばねで、互いに又は床と結合した系とし、非線形ばねは荷重・変形関係が二直線骨格曲線であり移動硬化則に従う復元力特性を持つとする。解析モデルと固有振動モード及びモデル図を図10、11に示す(m1=m2=1kg,k11=k12=k,k22=2k,f11 y=f12 y=f22 y=2N,γ=0.2,h=0.02)。外力は、図3に示す島根原子力発電所2号機原子炉建物(1次固有周期Ts:0.21s,2次固有周期:0.13s)で観測した鳥取県西部地震の床応答(NS−2R/Bch24を係数倍したものとした。床応答(NS−2R/Bch24)の一定塑性率床応答スペクトルに解析モデルの固有周期を追記した図を図12に示す。Model Aは建物の一次モードを超えた領域および建物の一次モードにより床応答の増幅が大きい領域(スペクトルの山領域)に固有周期を持つ系、Model Bは建物の一次及び二次モードにより床応答の増幅が大きい領域(スペクトルの山領域)に固有周期を持つ系、Model Cは建物の振動モードから外れ床応答の増幅が小さい領域(スペクトルの谷領域)に固有周期を持つ系である。これらのモデルで最初に降伏するばねは一次モード、二次モード共に、Spring2−2である。
ModelA,B及びCの塑性率の計算結果を図13A〜13Cに示す。EP−RSA(棒グラフとバー付き点線)、E−RSA(実線と一点鎖線)及びRHA(バー付き二点鎖線)は、各々、弾塑性応答スペクトル解析、弾性応答スペクトル解析及び時刻歴応答解析から得られた塑性率を示す。弾塑性応答スペクトル解析から得られた一次と二次のモード塑性率を各々棒グラフの黒と白で示しており、モード重ね合わせ則を絶対値和則とした塑性率は積み重ね棒グラフで、モード重ね合わせ則をSRSS則とした塑性率はバー付き点線で示す。弾性応答スペクトル解析から得られた塑性率は直線であり、モード重ね合わせ則を絶対値和則とした塑性率は実線で、モード重ね合わせ則をSRSS則とした塑性率は一点鎖線で示す。塑性率はモード重ね合わせ則を絶対値和則とした値よりSRSS則とした値が小さい。また、ELは系が弾性限界となる外力を示し、EL×1.5,2,3,6は系が弾性限界となる外力の各々1.5,2,3,6倍の外力を示す。あわせて、Model Aについては2次の想定荷重係数ε2=ε(a2 y)を、ModelB,Cについては1次と2次の想定荷重係数のうち値が小さい1次の想定荷重係数ε1=ε(a1 y)を、各ModelのSpring2−2の図に白丸付き点線で示した。なお、Model Aの1次固有周期はTs<T1であり想定荷重係数を定義していない周期帯にある。
Model A,Bは、時刻歴応答解析から得られた各ばねの塑性率が弾性応答スペクトル解析から得られた塑性率より小さく、地震力の増加に対する塑性率の増加が小さい傾向にある。これは、Model A,Bが床応答の増幅が大きい領域に固有周期を持つ系であり、履歴減衰の効果が大きいためであると理解できる。Model Cは、時刻歴応答解析から得られた各ばねの塑性率が弾性応答スペクトル解析から得られた塑性率より大きく、地震力の増加に対する塑性率の増加が大きい傾向にある。これは、Model Cが床応答の増幅が小さい領域に固有周期を持つ系であり、履歴減衰の効果が小さいためであると理解できる。
(1)系が弾性限界となる外力の6倍までの外力で加振していること
(2)Model A,B,Cは各々床応答スペクトルの特徴的な領域に固有周期を持つ系であること
(3) これらの解析結果にはばね1か所が降伏する応答、ばね2か所が降伏する応答、ばね3か所全てが降伏する応答を含んでいること
(4) 機器・配管系の耐震設計に適用する上限程度までの塑性率が系に生じていること
多自由度の非線形系より対応する線形系の塑性率が大きくなるか、判定する条件の有効性を解析結果から確認する。解析モデルと外力毎に(系が弾性限界であるEL×1は除く)、判定条件式(43)への適合性を確認した結果を表1に示す。ここでα=1.5とした。Model A,Bの全ての外力の解析ケース及びModel CのEL×1.5の外力の解析ケースが判定条件式(43)を満たす。図13A〜13Cより、これら判定条件を満たす全ての解析ケースで、弾塑性応答スペクトル解析より弾性応答スペクトル解析が大きな塑性率を与える(同じモード重ね合わせ則で比較)こと、さらに、時刻歴応答解析より弾性応答スペクトル解析が大きな塑性率を与えることがわかる。なお、Model A,EL×1.5,2,Model C,EL×1.5の弾塑性応答スペクトル解析においてモード重ね合わせ則がSRSS則の塑性率は、時刻歴応答解析の塑性率より少し小さいが、これはSRSS則の精度の影響と考える。
発明者らは、非線形多自由度系の動的応答解析法として弾塑性応答スペクトル解析法を開発した。この解析法は、線形多自由度系の動的応答解析法である弾性応答スペクトル解析法を、一定塑性率応答スペクトルを用いて非線形多自由度系の動的応答解析法に拡張したものであり、新たに定義した時間依存モード塑性率に対する非線形一自由度系の運動方程式を基礎とするものである。荷重・変形関係が二直線骨格曲線である非線形二自由度系について、弾塑性応答スペクトル解析法と時刻歴応答解析法により動的応答解析を行った結果は、一定の精度で合致した。弾塑性応答スペクトル解析法は、弾塑性応答する多自由度系の動的応答解析法として、また、このような系の応答を弾性領域のモード概念により理解する上で有効であると考える。
Claims (5)
- 塑性領域におけるモード間の相互作用を無視し、モードの重ね合わせにモード重ね合わせ則を仮定した上で、弾性領域の固有振動モードにより非線形多自由度系をモード展開して、非線形一自由度系の時刻歴応答解析により作製した一定塑性率応答スペクトルから得たモード塑性率を用いることにより、各モードの最大変形を把握する、非線形多自由度系の動的応答解析法。
- 更に、多自由度の非線形系に比べ、対応する線形系の塑性率が大きくなると判定された場合に、多自由度系の非線形応答解析の簡易法として、多自由度系の線形応答解析を用いる、請求項1に記載の動的応答解析法。
- 系を支持する構築物の一次固有周期をTs、非線形多自由度系のr次の固有周期をTrとしたとき、Ts≦Trであることを、多自由度系の線形応答解析を多自由度系の非線形応答解析の簡易法として用いることができる条件とする、請求項2に記載の動的応答解析法。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2017/015794 WO2018193564A1 (ja) | 2017-04-19 | 2017-04-19 | 非線形多自由度系の動的応答解析法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP6358404B1 true JP6358404B1 (ja) | 2018-07-18 |
JPWO2018193564A1 JPWO2018193564A1 (ja) | 2019-04-25 |
Family
ID=62904963
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017551722A Active JP6358404B1 (ja) | 2017-04-19 | 2017-04-19 | 非線形多自由度系の動的応答解析法 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP6358404B1 (ja) |
WO (1) | WO2018193564A1 (ja) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108460200B (zh) * | 2018-02-12 | 2019-07-16 | 青岛理工大学 | 基于性态双谱的多模态性能化抗震设计方法 |
JP7393968B2 (ja) | 2020-02-19 | 2023-12-07 | 三菱重工業株式会社 | 評価装置、評価方法及び評価プログラム |
CN113484911B (zh) * | 2021-04-30 | 2022-11-15 | 石家庄学院 | 一种基于地震动演化功率谱的地震位移反应计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001290848A (ja) * | 2000-04-05 | 2001-10-19 | Hitachi Ltd | 配管の地震応答解析方法 |
JP2015222232A (ja) * | 2014-05-23 | 2015-12-10 | 三菱重工業株式会社 | プラント機器の健全性評価装置及び方法 |
-
2017
- 2017-04-19 WO PCT/JP2017/015794 patent/WO2018193564A1/ja active Application Filing
- 2017-04-19 JP JP2017551722A patent/JP6358404B1/ja active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001290848A (ja) * | 2000-04-05 | 2001-10-19 | Hitachi Ltd | 配管の地震応答解析方法 |
JP2015222232A (ja) * | 2014-05-23 | 2015-12-10 | 三菱重工業株式会社 | プラント機器の健全性評価装置及び方法 |
Non-Patent Citations (1)
Title |
---|
田村伊知郎 外1名: "地震力に対する機器・配管系の弾塑性応答評価法", 日本機械学会論文集 [ONLINE], vol. 82, no. 835, JPN6017027684, 25 March 2016 (2016-03-25), ISSN: 0003714599 * |
Also Published As
Publication number | Publication date |
---|---|
JPWO2018193564A1 (ja) | 2019-04-25 |
WO2018193564A1 (ja) | 2018-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
De Angelis et al. | Optimal tuning and assessment of inertial dampers with grounded inerter for vibration control of seismically excited base-isolated systems | |
De Domenico et al. | Soil-dependent optimum design of a new passive vibration control system combining seismic base isolation with tuned inerter damper | |
Zhao et al. | Optimal design of an inerter isolation system considering the soil condition | |
Liu et al. | Effectiveness of fluid-viscous dampers for improved seismic performance of inter-storey isolated buildings | |
Hessabi et al. | Investigations of the application of gyro-mass dampers with various types of supplemental dampers for vibration control of building structures | |
Zordan et al. | Improved equivalent viscous damping model for base-isolated structures with lead rubber bearings | |
JP6358404B1 (ja) | 非線形多自由度系の動的応答解析法 | |
Gill et al. | Robustness of multi-mode control using tuned mass dampers for seismically excited structures | |
Akehashi et al. | Bounding of earthquake response via critical double impulse for efficient optimal design of viscous dampers for elastic‐plastic moment frames | |
Cimellaro et al. | Optimal softening and damping design for buildings | |
Losanno et al. | An optimal design procedure for a simple frame equipped with elastic-deformable dissipative braces | |
Xiang et al. | An extended modal pushover procedure for estimating the in-plane seismic responses of latticed arches | |
Li et al. | Effects of brace stiffness and nonlinearity of viscous dampers on seismic performance of structures | |
Talley et al. | Evaluation of steel building structures with inerter-based dampers under seismic loading | |
Li et al. | A novel tuned inerter eddy current damper: modeling, optimization, and evaluation | |
Kang et al. | Multi-modal seismic control design for multi-storey buildings using cross-layer installed cable-bracing inerter systems: Part 1 theoretical treatment | |
Rajana et al. | A novel nonlinear isolated rooftop tuned mass damper-inerter (IR-TMDI) system for seismic response mitigation of buildings | |
Behnamfar et al. | Enhanced nonlinear static analysis with the drift pushover procedure for tall buildings | |
Shi et al. | Second‐mode tuned mass dampers in base‐isolated structures for reduction of floor acceleration | |
Dai et al. | Inerter location effect on the generalized tuned mass damper inerter control | |
Ghorbanzadeh et al. | Lateral soil pile structure interaction assessment for semi active tuned mass damper buildings | |
Dezhdar et al. | Influence of ground motion scaling on seismic response of concrete shear wall buildings | |
Yao et al. | Real-time hybrid testing for active mass driver system based on an improved control strategy against control-structure-interaction effects | |
Tamura et al. | Seismic Evaluation Method of Piping Systems by Inelastic Response Spectrum Analysis: Part 1—Response Analysis | |
Bantilas et al. | Capacity spectrum method based on inelastic spectra for high viscous damped buildings |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20171002 |
|
A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20171002 |
|
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20171226 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20180109 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20180213 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20180522 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20180604 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6358404 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |