機能
- ファイバー要素を利用した非弾性時刻歴解析を行うために鉄筋とコンクリートの応力-ひずみ関係を定義します。 各モデルは提案者の示方規定によって異なります。
経路
メインメニュー:[材料/断面] タブ > [非線形特性] > [ファイバー要素] > [ファイバー要素の材料定義]
入力
梁要素の断面をファイバーに分割し、各ファイバーセルが特定の応力-ひずみ関係を持つようにします..ファーバー材料の特性を定義します。
ファイバー要素の材料定義 ダイアログボックス
追加
ファイバー材料の属性を新規に入力または追加します。
修正
入力済みファイバー材料の属性を修正します。
削除
入力済みファイバー材料の属性を削除します。
名称
定義するファイバ要素モデルの名称
材料タイプ
履歴モデルを定義する材料を選択します。
履歴ループタイプ
履歴モデルを定義するタイプを選択します。
コンクリート
Modified Kent & Park Concrete モデルであり横拘束鉄筋の拘束効果(Confinement Effect)などを考慮できます。
Kent and Park(1973) の基本モデルを Scott et al.(1982) などが修正したモデルであり構成の明瞭さと解析の正確性を十分に備えたモデルと評価されています。
: コンクリートの圧縮強度
: 拘束効果により圧縮強度の増加効果を表現する係数
: スターラップで囲まれたコンクリートの体積に対する帯鉄筋の体積比、 :スターラップ鉄筋の降伏強度
: 圧縮Crushing発生時の変形率
: 最大圧縮強度発現時の変形率
Z : 圧縮降伏後のコンクリートの軟化区間の剛性を示すための係数
( : スターラップされたコンクリートの短辺長さ)
(S h : スターラップまたはフープセットの中心間隔)
以下に示す Scott et al が提案した式がある。しかし、この式は保守的に判断した値なので、使用者が適切に判断して任意の値を入力しても差し支えない。
本モデルでは部材全般にわたるその影響が少ないため、コンクリートの引張強度は無視しています。
本モデルの履歴特性を詳しく整理すると、次のようになります:
a. 圧縮状態で除荷(unloading)が発生した場合、除荷(unloading)が発生し始めたひずみとひずみ軸上の一点をつなぐ直線上で挙動をすると考えます。この直線の傾きは、圧縮強度が高くなるほど小さくなり、その実装式は以下の通りです。
b. 本モデルでは引張強度を無視しているため、完全除荷(Complete Unloading or Open Crack)状態では応力を「0」と考えます。(上図参照)
c. 再載荷(reloading)状態では、εpに到達するまでは応力を「0」と考えます。まず、εpより大きい圧縮ひずみが発生すると、前の除荷(unloading)時の直線状に再載荷が進行すると考えます。実際には、このように直線状に除荷(Unloading)と再載荷(Reloading)が発生せず、非線形的な軌跡に沿って挙動します。しかし、そのような非線形的挙動の影響が小さく、解析上の効率性と正確性を考慮すると、本モデルのように直線起動に沿って、除荷(Unloading)と再載荷(Reloading)が発生すると考えるのも十分であると言われています。
圧縮強度を超えた後、軟化区間での剛性で実際の履歴の傾きである E=K fc Z を使用する方式もありますが、この値は負の数であるため、極限状態での収束効率性を低下させることがあるので E=1.d-15 に設定されています。
日本コンクリート標準示方書「耐震性能比較調査編 p.23」のモデルとして最大応力点を超えた軟化領域と残留塑性変形、除荷後の再載荷時の剛性低下効果を反映しています。
fc' : コンクリートの最大圧縮強度
εpeak' : 最大圧縮強度発現時のひずみ
本履歴モデルの 応力-ひずみ の関係は以下のとおりです。
ここで、
本モデルにおいて、Kは弾性剛性の残存率を意味します。圧縮ひずみが大きくなるほど、除荷(unloading)時の傾きが減少する現象を模写するための値です。
一般的な梁部材の場合、引張強度は無視するように提案しています。上記履歴規則による挙動特性を表すと、下図のようになります。
参考として、日本のコンクリート標準示方書によるコンクリートモデルを用いて、構造物の極限挙動を評価する際に、上記の剛性低減率値を基準とした状態を調査する式があります。
圧縮ひずみが大きくなるほど、載荷時の剛性低減率が小さくなることを確認し、部材の破壊程度を判断する基準です。現在、MIDASプログラムでは、この基準に基づいた調査方式も実装しています。
(Section Damage Check with Fiber Element 機能参照)
日本道路橋示方書同解説·V耐震設計編【鉄筋拘束コンクリート·p.161】のモデルです。地震動の種類によって、応力-ひずみの関係が異なります。本モデルは地震動と拘束鉄筋による影響、断面形状などの影響を反映しています。
地震タイプI:極限ひすみと最大圧縮強度時のひずみが同じになり、下降勾配(Edes)区間を持ちません。つまり、圧縮ひずみが最大圧縮強度に達すると、即時にクラッシュが発生すると考えます。
地震タイプⅡ:極限ひずみを示方規定式により算定し、下降勾配(Edes)区間を持つことになります。
:コンクリートのヤング係数
:コンクリートの設計基準強度
:横拘束鉄筋の降伏点
:断面補正係数
円形断面では α=1.0 , β=1.0、台形断面,中空円形断面及び中空台形断面では α=0.2 , β=0.4
:横拘束鉄筋1本あたりの断面積
:横拘束の間隔
:横拘束の拘束長で、帯筋または中間帯筋により分割拘束された内部コンクリートの辺の長さのうち最も長い値とします。
:コンクリートの引張強度
:横拘束鉄筋で拘束されたコンクリートの強度
一方、横拘束鉄筋の体積比を意味するΡsは、0.018以下になるよう基準で示しています。
本モデルの 応力-ひずみの関係は以下の通りです。
ここで、
除荷(Unloading)と再載荷(Reloading)時の場合は、弾性剛性で挙動すると考えます。(下図参照)
圧縮強度を超えた後、下降区間での剛性で実際の履歴の傾きである E=-Edes を用いる方式もあります。しかし、この値は負の数であり、極限状態での収束効率を低下させる可能性があるため、E=1.d-15に設定されています。
名古屋高速道路公社(名高社)のモデルで、「コンクリートを部分的に充填した鋼製橋脚の耐震性能調査(案) p.7」のモデルです。
:コンクリートの圧縮強度
:コンクリートの圧縮強度到達時のひずみ
:圧縮強度の増加を反映するための係数
:コンクリートの極限圧縮ひずみ
:コンクリートの最大引張強度発現時のひずみ
:コンクリートの引張破壊発生時のひずみ
:コンクリートの極限引張ひずみ
本モデルでは図に示すように、圧縮ひずみが極限圧縮ひずみに到達するまでは一定の最大応力を示し、変形率がεcuを超えるとクラッシュが発生すると考えます。
一方、除荷(Unloading) , 再載荷(Reloading)時には弾性剛性で挙動すると考えます。実際の仕様規定では引張強度に対する提案はありませんが、汎用性を備えるために使用者の判断によって適切な引張強度を入力できるように実装されています。(下図参照)
引張強度を超えた後、下降区間での剛性に実際の履歴の傾きを使用する方式があります。しかし、この値は負の数であり、極限状態での収束効率を低下させる可能性があるため、E=1.d-15に設定されています。
引張部と圧縮部の両方を定義できるモデルで、圧縮部は3線形履歴を持ちます。3 線形履歴を定義するために 応力-ひずみ で入力する方式と、応力-剛性低減率で入力する方式の2 つの形式があります。
応力-ひずみ 定義方式 応力-剛性低減率 定義方式
: コンクリートの1次圧縮降伏強度
: コンクリートの2次圧縮降伏強度
: コンクリートの2次圧縮降伏後の強度(K3算定時に必要)
: コンクリートの最大引張強度発現時のひずみ
: コンクリートの引張破壊発生時のひずみ
: コンクリートの極限引張ひずみ
: コンクリートの1次圧縮降伏ひずみ
: コンクリートの2次圧縮降伏ひずみ
: コンクリートの2次圧縮降伏後のひずみ(K3算定時に必要)
: コンクリートの3次圧縮降伏ひずみ
: コンクリートの初期剛性
: コンクリートの1次降伏後の剛性と初期剛性の比
: コンコリートの2次降伏後の剛性と初期剛性の比
引張強度を超えた後、下降区間での剛性で実際の履歴の傾きを使用する方式があります。しかし、この値は負の数であり、極限状態での収束効率を低下させる可能性があるため、E=1.d-15に設定しました。
「σ-ε」入力方式で、ε_c1~ε_c3を入力した状態で、「σ-α」入力方式を選択すると、自動的にそのひずみに該当する剛性 K1、K2/K1、K3/K1 を計算します。その逆も自動的に計算されます。
特徴
Mander (1988) によって提案された、横方向鉄筋で拘束されたコンクリートモデルです。横拘束鉄筋はコンクリートの拘束効果だけでなく、主筋の座屈を防ぎ、せん断破壊を防止する役割をします。横拘束鉄筋によって拘束された圧縮コンクリートは、横拘束効果によって強度と延性を大きく増加させる効果があります。
Mander モデルは拘束された圧縮コンクリートの応力-ひずみ関係を提案したモデルで、コンクリートの断面形状に関係なく適用できる長所があります。また、縦方向鉄筋および拘束鉄筋の量、拘束鉄筋の降伏強度と配筋形態などを考慮し、コンクリートの拘束力である有効拘束応力を算定することができます。
Mander モデルは、Popovic (1973) によって提案された 1軸応力 - ひずみ曲線に基づいています。ただし、Mander モデルを使用すると、多軸の有効拘束応力を 1軸の有効拘束応力に変換できます。
制限事項
Mander モデルはコンクリートの断面形状に関係なく適用可能です。ただし、材料データと断面データの Import による強度自動算定は、円形 , 四角形の断面のみをサポートし、円形 , 四角形の断面以外の断面については、強度自動計算に必要な情報をユーザーが入力すれば自動計算による強度が得られます。
コンクリートタイプ
非拘束:かぶり部のコンクリートを定義するときに使用します。断面の横拘束効果を考慮せず、横拘束に関する事項は無効になります。
拘束:横拘束鉄筋によって拘束されたコアコンクリートを定義する時に使用します。断面の横拘束効果を考慮します。
非拘束コンクリートデータ
:非拘束コンクリートの圧縮強度を入力します。直接入力する場合は、チェックボタンを選択して入力します。
:f'co に対応するコンクリートのひずみを入力します。直接入力する場合は、チェックボタンを選択して入力します。
:コンクリートの弾性係数を入力します。Mander モデルの場合は内部的に計算され、ユーザー入力の場合は直接考慮することができます。
:コンクリートの引張強度を入力します。引張強度を直接入力したり(ユーザー入力)、無視(Neglect Tensile strength)することができます。
:引張強度に応じて自動計算されます。
断面 & 拘束鉄筋タイプ:コンクリートタイプが "拘束" と定義されると、このオプションが有効になります。横方向鉄筋の断面タイプおよび鉄筋タイプを定義します。
断面データ:拘束断面の大きさと横拘束鉄筋の間隔などを入力します。
鉄筋データ:断面の主鉄筋と横方向拘束鉄筋の鉄筋面積に必要な情報を入力します。
拘束有効係数, Ke:断面形状 , 縦方向鉄筋および拘束鉄筋の量 , 配筋形態等を考慮して算定したコンクリートの有効拘束係数
コンクリートの有効横拘束応力
flx:x方向の有効拘束応力
fly:y方向の有効拘束応力
fl:多軸の有効拘束応力を1軸に還元した有効拘束応力
拘束コンクリートの強度 & ひずみ:横拘束効果を考慮した圧縮コンクリート強度とひずみ
鉄骨
Menegotto and Pinto の 鉄骨モデルを Filippou などが修正したモデルで、解析効率性にも優れ、実験的挙動とも高い正確性を持つモデルとして知られています。
:鉄筋の降伏強度
:鉄筋の初期剛性
:降伏後の鉄筋の剛性と初期剛性の比
:降伏後の鉄筋の応力-ひずみ曲線挙動形態を定義する定数
本モデルの詳細な 応力-ひずみ 履歴モデルの説明は以下の通りです。
上記の関係式を通じて、降伏後の履歴挙動から弾性と降伏区間を転移する際に曲線形態を具現しています。
ここで、σ' と ε' は正規化(Normalized)された値であり、それぞれ次の数式によって計算されます。
,
本モデルは、下図のように2つの漸近線を有することになります。1つは弾性剛性を傾斜に持つ直線で、もう1つは降伏後剛性を持つ直線です。2つの直線の交点を意味するのが(ε0 , σ0) です。(εr , σr) は最後に除荷(Unloading)が発生したときの位置を意味します。これらの値は、除荷,再載荷時の度に更新され、転移曲線形成に影響を与えます。
Rは、転移曲線の形状に影響を与える値であり、Bauschinger effect を反映する役割をします。R値を決定するために、次の式が提案されています。
係数 a1 , a2 , R0 は実験から得られた履歴結果から求められ、MIDASでは、原文献(Menegotto and Pinto , 1973)で使用された値である18.5 , 0.15 , 20 をそれぞれの基本値として使用しています。
ξ は正規化されたひずみ軸上で2つの漸近線と除荷発生位置との間の距離を意味する値で、除荷(Unloading)発生時ごとに新しく更新される値です。
一般的な対称バイリニア鉄筋モデルです。
:鉄筋の降伏強度
:鉄筋の初期剛性
:降伏後の鉄筋の剛性と初期剛性の比
下図のように降伏後の除荷(Unloading)、再載荷(Reloading)時には弾性剛性で挙動します。
一般的なバイリニア鉄筋モデルを発展させて、降伏後の鉄筋の剛性を任意に定義でき、鉄筋の座屈や破断などを考慮できます。
:引張側降伏強度
:圧縮側降伏強度
:鉄筋の圧縮座屈発生時のひずみ
:引張降伏後鉄筋破断発生時のひずみ
:鉄筋の初期剛性
:引張降伏後の鉄筋の剛性
:載荷時に鉄筋の降伏後の剛性 (圧縮領域で引張力が無負荷になり再負荷されると、E3 ラインが圧縮荷重下での剛性を制限します。)
:圧縮降伏後の鉄筋の剛性(負値を指定すると負勾配を考慮できます)
:圧縮降伏が生じた鉄筋の座屈後の剛性
本モデルは鉄筋の圧縮-引張降伏、引張破断、圧縮座屈など様々な状態まで再現できるモデルです。発生可能な履歴状態を下図に示します。
各状態に対する詳しい説明は次の通りです。
引張挙動が主要な鉄筋履歴なので、荷重載荷方向は引張を基本として考慮します。 したがって、制荷(Unloading)は引張状態から圧縮状態に荷重方向が変わったことを意味し、再載荷(Reloading)は圧縮状態から引張状態に荷重方向が進むことを意味します。
State 1:弾性挙動状態で、E1の傾きで弾性挙動します。
1 → 2 = 降伏状態に転移することを意味し、引張降伏時にはE2の傾きで、圧縮降伏時にはE4の傾きで挙動します。
State 2:降伏後の状態でE2を傾きで動きます。
2 → 4 = 降伏以降に除荷または再載荷される状態を意味します。
2 → 8 = 降伏後に引張が持続して引張破断が起きた状態を意味します。破断後は応力が無条件で「0」と計算されます。
State3:除荷が続き、圧縮部でも降伏が発生した状態を意味します。傾きE2で挙動し、E3直線とひずみ軸との交点が ε2 より大きいようにE3を入力します。
3 → 4 = 再載荷が起きる状態を意味し、E1の傾きで挙動することになります。
3 → 5 = 除荷が続き、圧縮ひずみが ε1 を超過して圧縮座屈が始まることを意味します。E5の傾きで挙動します。
State 4:除荷、再載荷時に弾性剛性E1を傾きで挙動する状態を意味します。この状態での除荷 , 再載荷は E1の傾きで進行されます。
4 → 2 = 再載荷が進行され、引張降伏状態に移る過程を意味します。 または、除荷が進行して圧縮降伏状態に移る過程を意味します。
4 → 3 = 除荷が進み圧縮降伏が発生した状態を意味します。E3 が E2 と異なる場合でも、上記図に示すように、E3 直線と一致すると圧縮降伏が発生したと判断します。
4 → 5 = 除荷が進み圧縮座屈状態に入る過程を意味します。
State 5:圧縮ひずみが座屈発生ひずみを超えて鉄筋の座屈が発生した状態を意味します。傾きE5で挙動します。
5 → 4 = 圧縮座屈の途中で再載荷が進む状態を意味します。
5 → 7 = 圧縮が続き、鉄筋が完全に圧縮座屈した状態を意味します。State7 以降は、圧縮荷重下では応力が「0」となります。
State 6:圧縮座屈以降に再載荷される状態を意味します。この時、引張側降伏が発生する以前であれば σy に向けて再載荷が進行されます。
引張側降伏後であれば、引張側最大点に向かって再載荷が進行されます。
6 → 2 = 再載荷が進み引張降伏が発生した状態を意味します。
6 →-1 = 再載荷が発生した状態を意味します。
State 7:完全な圧縮座屈が発生した状態を意味し、これ以上圧縮応力を発生させられない状態です。圧縮応力は「0」となりますが、引張力には抵抗できます。
7 → 6 = 再載荷が発生し、引張最大点に向かって進行する状態を意味します。
State 8:引張破断が発生し、これ以上引張応力を発生させられない状態を意味します。 引張破断が発生すると、これ以上引張 , 圧縮のいかなる応力も発生させられません。
State-1:圧縮座屈後の再載荷過程(state6)で除荷が発生した状態を意味します。弾性剛性E1で挙動します。
-1 → 6:再載荷が発生して再び引張最大点に向かって進行する状態を意味します。
-1 → 7:除荷が行われ、完全に圧縮座屈状態に転移する状態を意味します。
以上のように多様な項の状態と転移を考慮したモデルなので、使用者の十分な理解に基づいて実験的数値を適切に使用すれば、効率的なツールになります。しかし、引張破断や圧縮座屈後のような極限状態では抵抗応力が全くなく「0」になるので注意が必要です。
3つの傾きを持つトリリニアモデルです。
履歴入力方式の利便性のために、履歴入力を 応力-ひずみ の座標で入力する方式と、剛性の低減率で入力する方式の2つを実装しています。
除荷、再載荷時の弾性剛性を基本に挙動することになります。
応力-ひずみ 定義方式 応力-剛性低減率 定義方式
:引張側1次降伏強度
:引張側2次降伏強度
引張側の2次降伏後の強度 (K3算定時に必要)
:圧縮側1次降伏強度
:圧縮側2次降伏強度
:圧縮側の2次降伏後の強度 (K5算定時に必要)
:引張側1次降伏ひずみ率
:引張側2次降伏ひずみ率
:引張側2次降伏後変形率 (K3算定時に必要)
:圧縮側1次降伏ひずみ率
:圧縮側2次降伏ひずみ率
:圧縮側の2次降伏後の変形率 (K5算定時に必要)
:鉄筋の初期剛性
:1次引張降伏後の鉄筋の剛性と初期剛性の比
:2次引張降伏後の鉄筋の剛性と初期剛性の比
:1次圧縮降伏後の鉄筋の剛性と初期剛性の比
:2次圧縮降伏後の鉄筋の剛性と初期剛性の比
「σ-ε」の入力方法において、ε1y~ε'3y の値が入力され、「σ-α」の入力方法が選択されると、そのひずみに該当する剛性 K1 , K2/K1 , K3/K1 などが自動的に計算されます。また、その逆計算も自動的に行われます。
Kent&Park(1973)によって行われた繰り返し荷重を受ける軟鋼の実験を通じて提案されたモデルです。本モデルは軟鋼の弾性区間と塑性区間とひずみ-硬化(Strain Hardening)区間の表現が可能で、Ramberg-Osgood式によってBauschinger Effectを精密に示し、実験的結果と高い一致性を示すモデルです。
:鉄骨の降伏強度
:鉄骨の終局強度
:鉄骨の初期剛性(弾性係数)
:鉄骨の降伏ひずみ
:ひずみ-硬化開始点のひずみ
:鋼材の終局ひずみ(破断時)
Case1:載荷時の挙動
載荷時の挙動は次のように区分されます。応力-ひずみの関係は、Thompson & Park(1980)が提案した式を適用します。
同じ符号の荷重を載荷した鋼材の応力-ひずみ曲線
:鉄骨のひずみ
:鉄骨の応力度
:鉄骨の降伏強度
:鉄骨の終局強度
:鉄骨の初期剛性(弾性係数)
:鉄骨の降伏ひずみ
:ひずみ-硬化開始点のひずみ
:鋼材の終局ひずみ(破断時)
Case 2:除荷および再載荷時の挙動
除荷時および載荷時の挙動は、Ramberg-Osgood 関係によって定義され、Newton's Method による反復計算によって応力を求めます。
逆荷重を受ける鋼材の応力-ひずみ曲線
:Ramberg-Osgood関数の特性応力度
:以前の載荷時の塑性ひずみ (0< εip <0.7097)
:Ramberg-Osgood Parameter
:Loading Run Number (ただし、圧縮側の場合は1、引張側の場合は2の固定値を使用)
:載荷時点で応力0に対するひずみ
履歴モデル
ファイバー要素に使用できる履歴モデルは鉄筋用4本、コンクリート用5本の計9本です。
スケルトン曲線
各モデルの履歴ループを参照して、材料の挙動履歴に対するバックボーン曲線の特性値を入力できます。