( A . I O )
(A.ll)
一方、半径方向および接線方向の変位UrおよびU()はそれぞれ変位とひずみの関係を積分すること より求められる。すなわち
ここで
Ur
=
F ](N rN(S‑2)+lur(8)旬
。 =
F ](N rN(S‑2)+lug(8)Ur( 8) =
u
g
(8) =ã~
‑1 { S( s ‑ D ~ ‑ ~ } ‑
Ur (B)さらに (A.2)式のひず‑みんJは次式で与えられる。
また (A.2)式中のIは次式で与えられる。
{7r ( ( N ¥. I¥Tけ
= J / ‑ π l
~( ¥N 一 一 ‑ + : ‑)
1) a‑a :
eq +1 cos 8 ‑sin 8 δ[r (U()一弘)‑
Tr( )(仏+
Ug)](A.12)
(A.13)
(A.14)
一(古)
cos B (a r ur+ 同 ) }
dB (A.15)ただし、上記の IとJ値、 Kの聞には以下に示す関係がある。
J = F ](N+lriN+1)(S‑2)+1] (A.16)
ここで rl :き裂先端から J積分の積分経路までの距離
以上の結果よりゆが与えられればσtJ、εりが求められる。ここで実際に解くことになる微分方程T を示しておく。これは (A.9)式を(A.6)式に代入することによって求められる。
(S ‑2)ーこす
I
lã~-l ~ S(S ‑1)ふ ‑
2~" ~[J2l a
d()2 I L ‑eq l ‑¥ ‑ ‑/ T ‑T J
+ {N(S ‑2)
+ 川
N(S‑2)}ã~ぺ S(2S
‑3) ~
‑~.. }
+6{N(S‑2)+
山
Sー 1)(〈リ)'
= 0 (A.17)ただし
(2)平面ひずみ状態
ψ ( 0 )
=ゆ( 0 ) = 0
ゆ(π)=ゆ(π)
=
0 (A.18)平面ひずみ条件についても全く平面ひずみ条件と同様な解法を用いて求めることができる。 平 面応力条件との違いは、平面ひずみ条件における相当応力が
σ ? q = ; ぽ ‑ a t ) + 3 r : e
となり、特異項の支配方程式が以下のようになることである。
(品一号‑;守 } { σ 2 ‑ 1 ( ; 〆 + シ ‑ < t " ) }
+4 占 ( 4 ‑ 1 ( ; け ) f = o
(A.19)
( A . 2 0 )
上記の境界条件は平面応力状態の場合と同様に (A.7)式であり、また同様にσT、旬、アreは応力関 数制こより (A.8)式のように表される。同様な手法で(A.9)式を(A.8)式、 (A.19)式に代入しδり
を求めると以下のようになる。
また、変位Uijは
δr(())=Sゆ(())+ゆ(())
δ 。 ( ( ) )
=5(5‑1)
ゆ(())子re(())= (1 ‑
S )
ゆ(())δeq(()) =
J~ ( ム ( ( ) )‑
ae(()))2+
3子f o
S =‑
竺土 1
N
+
13ã~-1
Ur(q V J ‑()) = 4(2‑S)
L1;~ e~
.Ci {S(2 ‑S)れふ)
3(S ‑3)δN‑l
Ua(()) 二 eq(S(2‑S)
れる)
e¥V) ‑ 4(2 ‑s)
( A . 2 1 )
( A . 2 2 )
ひずみ九は
3 ã-~-1
σEL ハ {5(2‑5)
れふ)
(A.23) ie =
~ã-~-l
{5(5 ‑2)~ +ふ)
のようになる。また、 (A.2)式中のIは平面応力状態と同様に(A.15)で求められ、最終的にσリ、
つ ハ リ
よ2こA ル
入 す る ﹂
代 日 川
に
1 φ 式 +
︑1
j 7 A V
FE︐ ︑
︑︑ ︐ ︐
J
2 S
A一
を や 式 わ い
Lー
j 1 0
M M J
=
f
ぃ
δγi /
は ー ー ー
JAV
式 斗
↓
程+Nrq
方 り ム
¥
分 d
一 け
微
S +
T ι / 1
¥
︑1
︐ ノ
ル 砂 川 い 斗
f 幻
S︐ ︐ ︑ ︑
〆' s
t︑︑
I r
一N角
S J t
ψ
ふN ] る れ 一
︑ ト
し
︑ り
2
一 戸 わ し い
用一ぽ何同﹂什
め に た う る よ め の 求 下 を 句 以 て
また、上式の境界条件も平面応力状態と同様に(A.18)で与えられる。
HutchinsonはIについてはN 二 3
,
5ヲ9,
13の場合についてのみ解を与えている。またんJに ついては N = 3,
13の場合のみ解を示している。しかし九などについては解を示しておらず¥(A.2)式を用いて種々の強度レベルの応力 ひずみ分布を求めることが簡単にはできない。そこで
は
.18)式の境界条件下で (A.17)、(A.24)式をルンゲ・クッタ法により解き、上述した平面応力状 態、平面ひずみ状態の種々のa ‑
ij、iij、Uijの値を求めた。破壊挙動を検討するとき、き裂線上のひずみや応力が問題となるので、。 =0におけるムハムj、Uij
などにつき、
N '
こ対する多項式近似を求めた。平面応力状態の結果、平面ひずみ状態の結果をそれぞ れまとめ、前者を(A.25)式に後者を(A.26)式に示す。またこれらの結果をFig.A‑1rv Fig.A‑5 に示す。図中の・印は Hutchinsonが与えている結果である。1
= 2.5034+
5.1180/N ‑3.3272/N2a ‑ r ( O )
0.57230+
0.62131/ N ‑0.15221/ N2a
‑ e ( O )
一 1.1516 ‑0.0098350/N ‑0.37225/N2子
r e ( O ) 。
ι(0) ‑0.0071570十0.63674/N ‑0.020181/ N2
ゐ(0) 0.82020 ‑0.21718/ N ‑0.41565/ N2 (A.25)
U r ( O )
一 0.47560+
1.3502/N ‑0.92073/N2弘(0)
。
む(0)
。
U o ( O )
一 0.22550 ‑4.4307/ N2+
9.0064/ N4九
q ( O )
一 0.99755 ‑0.013207/ N ‑0.10327/ N23.8830
+
6.8988/N ‑ 5.6685/N21. 7865 ‑0.60084/ N
+
0.39441/ N22.8730 ‑3.9361/ N
+
3.1427/ N2 一一
(A.26)
。
‑0.0057542 ‑0.10816/N
+
0.21970/N20.0057542
+
0.10816/N ‑0.21970/N2‑0.27209
+
1.0545/N ‑1.1219/N2。 。
一 一
一
。
r(OI ) δ。
(0)子rO(O) εr(O) [e(O) Ur(O) ι(0) ω(0)
む
(0) 0.27785 ‑0.94631/N +
0.90215/N
20.94099 ‑2.8884/ N十2.3801/N2
。
εq(O)Plane strain condition Plane stress condition
4
/
• After Hutchinson o Present analysis 5ωコ
一町
﹀一
• After Hutchinson Present analysis
。
3
ω
コ一回
﹀
4 8 12 16 20 24 28 Strain hardening exponent : N
。
1=2.5034+5.1180/N‑3.3272/N2
4 8 12 16 Strain hardening exponent : N
2 Q
Fig.A‑l Calculation results of 1 value .
g
l b
民 1.1
l b
0.8 Jゴ
lb 0.6 1.2
0.9a
O
.4a
Plane stress condition
1ト吉岡=0.94099‑2.8884/N+2.3801/N2
\~円~,....,円台 eeee
σo=1.1516‑0.0098350/N‑0.37225/N2
¥
gl b
~-\
σo=2.8730‑3.9361/N+3.1427/N2
│
: : │ ;
l b ""'‑J10 σ。q1
""'‑J
16σ。│
27
伺=0.99755‑0.013207/N‑0.10327/N2
。
Strain hardening exponent : N Strain hardening exponent : N
Fig.A‑2 Calculation results ofσ() and σeq along a crack line .
Plane stress condition
ur=O.4 7560+ 1.3502/N‑0.92073/N
同
2日
iゴ σr=1.7865‑0.60084/N+0.39441/N2l b
""'‑J
│ σ
i
lへ J
1 0 ur 1
。
ur=/
ー0.27209+1.0545/N‑1.1219/N2σr=0.57230‑0.62131/N‑0.15221/N2
。
Strain hardening exponent : N Strain hardening exponent : N
Fig.A‑3 Calculation results of ur and σr along a crack line .
いご 0.4~
0.2ト
。
0
.40.2 2:)'
jゴ O
Jコ
‑0.2
‑O.4
a
Plane stress condition