イ
JJfFlり い )以下の式を用いて節点変位が次のように計算されるO
( 8
.12 )
‑ 9 1 ‑
U ; l ) ( t + 1 ' 1 t ) = u ? ) ( t ) + d ω ? ) ( t + 1 ' 1 t ) ( 8 . 1 3 )
(h)幾何学的更新の解析
大変形モードで、節点位置は以下の式を用いて更新される。その後に更新位置に基づい て、グリッドの変形が起こる。
x ? ) ( t +
l'1t )二x i
(l)( t ) + d t V ? ) ( t + d t ) ( 8 . 1 4 )
ただし、微小変形モードでは、幾何学的修正が無視される。
8 . 2
.4 モール・クーロンモデル2)土の供試体に外力を加えると変形が生じる。外力が等方的であれば体積が収縮するだけ であるが、せん断力が加わると形状変化を生じる。一般に、ひずみの小さい範囲ではせん 断応力とひずみの関係は直線的であるが、せん断応力が増大するにつれてひずみの増分が 大きくなり、ついにはせん断応力は増大しないでひずみのみが増大するようになる。この ような状態を一般に土の破壊という。モールの破壊規準線は一般には曲線であるが、土質 力学において、土の破壊包絡線が直線的になると数学的扱いが簡単になることから、この モデ、ルの破壊包絡線は直線としている。つまりモール・クーロン規準は、土のせん断強さは 応力状態に関係しない粘着力成分とせん断面に作用する垂直応力に比例する摩擦成分の和 として考えられる。モール・クーロン規準において引張力を考えるため、引張強度と、粘着 力
/ t a n φ
の小さい方を引張強度とする。( a )
応力せん断要素FLAC3D
でのモール・クーロン規準は、応力σ],a 2,σ3によって表される。それはこのモ デルの、3
つの応7Jベクトルの構成要素である。またそれに対応するせん断ベクトルの構成 要素は、ひずみ E ,]E 2, E 3である。( b )
増分弾性則応力または増分応力の式はフックの法則を適用し、次式で表される。
d
O"]= α ] d & t + α 2 ( d & ; + d & ; )
Aσ2
二α ] d & ;+ α
2( d & ] e + d & ; ) ( 8 . 1 5 ) Aσα]d&; + α 2 ( d & t + d & ; )
ここにα l、α 2はせん断弾性係数
G
と体積弾性係数K
により定義される物質の定数でありα
,=K+fG
3 ( 8 . 1 6 )
‑ 92 ‑
α
刊 二K‑2Gゐ
3
とする。
また、次式により
AEz=St(AE;)
i=
l,n ( 8 . 1 7 )
次のようになる。S I (~&; ,~&; ,~&;)=α1~&; + α 2 (~&; + ~&;)
S 2 (
企&;,~&;
,~&; ) = α I
企ε ; + α 2 (
企Ef+U;) ( 8 . 1 8 ) S 3 (~&; ,~&; ,~&;)=α1~&; + α 2 (~&; + ~&;)
( c ) 2
つの破壊規準FLAC
モデ、ルで、用いられる破壊規準は、引張り部分を切り離したモール・クーロン規準で ある。その場合、3
つの主要応力を次のように決める。σ1::;σ2 : : ; σ ( 8.1
9 )
図・8. 2
が示すように、この規準は( σ 1
,σ 3 )
平面で表される。(圧縮応力を負とする。)破壊 包絡線K σ 1
,σ 3 ) = 0
はモール・クーロン規準により、点Aから点 Bまでf"= 0
となる。f S =
σ1 ‑0'"3N O + 2c . j i i ; ( 8
.20 )
引張り破壊規準により点
B
から点C
においてft= 0となる。
λ=σ3 一 σr ( 8 . 2 1 )
ここに、
φ
は摩擦角、c
は粘着力、σt
は引張り応力であり、N
ーと皇位}
世 一
1‑ s i n ( O ) ( 8 . 2 2 )
材料の引張り強度が直線fS =
0
と( σ 1
,σ 3 )
平面の交点の縦軸σ 3
の値を越えることができ ないことに注目する。この応力の最大値は次式で与えられる。σ ; … 一 一 = 一 一 t 一 a c 一 n o ( 8 . 2 3 )
位置の関数は、
2
つの関数により記述される。それはg sと g t
であるが、これらはそれぞ れせん断塑性流動と引張り塑性流動を定義したものである。関数g s
は、次式で与えられる。g S
ニσ σ3N
'I'( 8
.24 )
守U
Qd
. o
L/ /
/
/ / /
尚一
駒 場 毒 毒 殺 E 8
・ 図 /
o e
/ /
, /¥
喧 /司
、
/
i 鞠島
家 百平 醤 組/
‑ f i‑ て ロ
・ ぺ/
ー モb 困 Z 8‑
l D
/ '
~ ~ d "ふ う' 晦傭 繊暢
/
プ ず/
0 '
お "
崎9 /' 予 j / J q
、;J/
/
/ /
D S
/ /
/事¥」
/司
、
〆'
れソ
あで角
シン
タ ー
レ
¥ V
一 一
⁝
v
イ削一刷
︑ダ .臼 一. 幻
土+一一
pb '' za
一'EA︐ψ
に
N
v 一 一( 8 . 2 5 )
である。関数
g t
は、次式となる。gt= σ3 ( 8 . 2 6 )
流動則は、以下の技術適用で固有の定義を与えられる。関数
h(σ]
,σ3)=0
は、図・8.3が示 すように、(σ]
,a3 )
平面上で、fS=
0とfl=Oに挟まれた斜線部分により表現される。図が示すように、関数は正か負の領域で選ばれ、次式でとなる。
h= σ 3
一σ
t+aP ヤ
1ーσ P ) ( 8 . 2 7 )
。
P二戸可 +N o
( 8 . 2 8 )
σPニσt凡 な 町
合成降伏関数の中での弾性則は、
(σ]
,a3 )
平面において、それぞれ正か負に対応する領域1
または領域2
のどちらでも表現できる。もし領域1
ならばせん断破壊と断言され、ポテン シヤル関数どを用いて導かれる流動則を用いる曲線fS=
0にその応力点はある。もし領域2 ならば引張り破壊となり、g l
を用いて導かれる流動則を使ったfl= 0
と応力点は一致する。(d)塑性調整
はじめに、せん断破壊を考えると、式(
8
.23 )
の偏微分より、。 gs1
0σl
金二。
θσ2 ( 8
.29 )
。
IgS一 一=‑N... Bσ '1'
。
IgS/θ σ
l' 8gS /θ σ
2' 8gS/8σ3
の変わりに企ε . ;
, 1'1&;, 1'1ε ;
を代入し、それぞれ(8 . 1 8 )
より、
s,(笠笠竿)~一一一 。 σ ] '
8σ2 ラー ~2" f//‑ 95 ‑
任まま)ヰ
21 ( 一九)
s j 8g
S生:生~1 一一_.
'A. T 盆α
3 18σθσ2θσ3 ) ‑ " " 1 " '
'11 I" " 2
( 8
.20 )
式から、σ
NσJ‑f(αl 一 α 凡)
σf=σ;‑λsα2 ( 1 ‑ N )
σ
Nσ;λ S ( ̲ α 1N
'II+ αJ
λ
fS (σLσ:)
一
( α l
一α 凡 ) ‑ ( ‑ α 1N
'II+ α
2) N
O次に、引張り破壊を考える。