図5.3圧縮二軸応力荷重の模式図(非定常変形の例)
1 .4
1 .2
1
0.8 o)
0.6 0.4 0.2
...d‑.... ¥
CFL=1 2 l' b¥ ' ' / f ' ' ¥¥ ¥ " I ¥
l:.‑..
.d‑.
f / p ¥ ¥ 't r
I' / . I = ' 5 I ¥
/ CFL O 8:¥t
t I ')
OFL I O i 1 IL ' ,1 , ¥
cl IL ,
l ¥ t,l ' I tlt ' I Il" i ' I alt t 'l t' l *
l CFL=0'5' 't l
(¥. tt
$ ; t:; L :
0.5
‑0.5
lgl
l 5.4 , ‑‑ ‑f o)i =tp ! 7 / )f l
CL1‑
¥
D: r
40 35 30 25 20 15 10 5
1
‑5
t I
I
, c as e 6
J‑llll"Illll'l'll‑1'1111‑Ill‑1'1'1'1111"I'lllall‑Ill‑ll
I a I
I I I
I :
I
vaSeu
a
l‑
lll‑'111'111'11111‑lll'lIIllIIllIII'll'lll'lI l
Qa )4
r
ll'll'lllIII l"I l' Il 'I 'I' Ill 'II 'I' l‑ 111 11 111 'IIrL O
vaseo
l l' 11 l' l' II' ll' I' I ll' l' l' Il ‑' Il I ll' ll' I a' I ll' l
l case2
['1111111111'11111111'11'11111111111lalllllllllllI
,
‑2
o 2 4 6 8 10
t [ms]
' 5.5 EE )B * :
Bi‑axial stress loading
l l I
P
IlIII III
<
Initial state
I
l
I I
I !I
l
l l l
・I
l l
I‑
=1
<
Length
l l l
1lt ll ll
‑ ‑ ‑ ‑ ‑ ‑
Oompression force
l 5.6 EE f I t; 4 ; f o) J! ) *.*',./' d l 'l
x
y
Bi‑axial stress loading End wall
LJ :
30 20
10 o‑1 O
‑20
‑30
r 11‑" 1'‑‑(r:b Sr ll"I:lll l': 1‑ l" "
rtT Ib l ll‑̲"b I '('d""'1 d , ' d db' l"I' 'I ,,... ll‑ It'
l I NI'll'
l """'J"' INl‑'̲Itl‑'1'
l l l (b ,(' al"N'I"
tbl I I l I l I I I I l l I
.
l.
l . l
. l
. l
. l
. I
. , I l l l , l l I , I l l I I l l l , l l l l l l l l I , I
I l
, l l : l l l
l l l
I I I l I
l l
I l l l I l
l l I l l l l l
l l l l ,
l・lll・1,II
lllll‑1・I
Ill‑1‑,
t=0.05 ms t=0.25 ms t=0.50 ms t=0.75 ms t=1 .OO ms t=1 .25 ms
,l
."tl"'IdL d'a"d"I dPIIPJ"d‑' II‑
, '1p
CIPa・1 ‑'1
l '‑1' l" ̲,1" I . I d' 'I d"'1 l l I I'dl I.
olp'l' ,.,,
J !':Zb. a=d' dlb :' t'H dla t: 1"'1 :‑':. I‑ ' " I L
o
20 40
60 80 1 OOx [mm]
t=0.05‑ I .25ms
l 5.7 LB 5 cUeC e J q) (t̲0.05‑1 25ms,P4 1 5MPa Case 4)
7
L‑1'
E
::
30 20
10 o‑1 o
.20
‑30
Bi‑axial stress loading End wall
o 20<
40 60
80 1 oox [mm]
t=1 .50‑2.75ms
I 5.8 A lJ}C b )E , ;; (t‑1.50‑2.75ms,P4=1.5MPa,Case 4)
‑ E E
Bi‑axial stress loading End wall
30 20
10 O‑1 o
‑20
‑30
displacement
r
:
l l l I
l l l l
l l l l I I l
1,
l I l
l , l
l I
, I I
t l l
I I l
I I l
I , l
I I I
I , l
t t I , l I I l l l I l , l l l , I I l I , l t I , l , l I , , l l I
I .
I .
. l I .
. I i
. l l .
l .
. l
l l I I I l l I l l l l I I I I l l
a.dll d'd' "" " 1' 1'‑ d'l
.1・‑ 'l‑ 1 1‑' l‑ I ‑IJ 't.dl‑11‑' ' '
Illlt・・Ill'l"It・・・1'1"'1'111""II'lll'll lb,‑'1‑ll:!‑・lll‑el.‑..llll'lll,.1'111'11'1lal'll'l l‑・.1・,lll‑It,ll'.III'l'l・1'‑‑,1'11'‑1'1'
l I I Il I ll・・・‑・It l・1・llll, IIllllll・・, lll‑,llll,
t=3.00 ms t=3.25 ms t=3.50 ms t=3.75 ms t=4.00 ms t=4.25 ms
Wd・lJ,III‑‑・ll‑,・‑‑' 'W・Nl,・・1,‑IIII‑
,1‑ ,d・11‑̲:‑1・ll・・・・1,11,・d‑,・lll・・・・‑1・1・‑・ll・llllll・lllIII l,Il・ll・II・1・1・・Jll,ll・la・a,Il・・1・1b・1・・・・lll
e・:II:It llb ,・・‑.Qlb・lNllp.ll‑lll‑1
"eb
"'・lb..,̲l,・̲.I
,‑W,‑1・・l
O 20 40 60 80 1 OO
x [mm]t=3.00‑4.25 ms
l 5 9 Ue 5by ) (t=3.00‑4.25ms,P4=1.5MPa,Case 4)
‑ E E
‑
Bi‑axial stress loading End wall
30 20
10 o‑1 o
‑20
‑30
displacement
l
: l l I l
: l I l
: l
I
1 l 'I
l I I I t
Il・IIII‑
l・1・llll, l,lllll・t
IIIllll
t=4.50 ms t=4.75 ms t=5.00 ms t=5.25 ms t=5.50 ms t=5.75 ms
O 20 40 60 1 OO
t=4.50‑5.75ms x [mm] 80J 5.lO l B / . cUeC b ; i q) ; (t 4.50‑5.75ms,P4 1 5MPa Case 4)
Bi‑axial stress loading End wall
E
IJ
>
30 20
10 o‑1 o
‑20
‑30
displacement
I l
l I
l l I I l I l l I l l l l I t l I l l I I I I l I , I , I I l I l I l I
l l
l l l : l I l l l l I , l I I l
, I l
I l l l l
l l l l l I l l I l l l l l l l l l l
RHII:lldrlllEllld 'Lll"'‑'11111'=11: :1' ' ::1'1 :llll :1'1'. .̲
l"'1111111'11"Ililll'lll'll"Il'llllllll'll"'1'1 "‑"I""‑1lt"' ̲1'1‑1'1‑'11'111:lll'l l"""'11'‑‑llil'llllll̲llll t=6.00 ms
‑ ‑ ‑ ‑ , ‑6 25 ms
t .
‑II・lll‑It Illlllllt
llllll,
I‑ I I I t
t=6.50 ms t=6.75 ms t=7.00 ms t=7.25 ms
l'Idlllll,1・1・1"I'l
.L .・・i:̲1'11'.1‑d・1' ‑1'1.1.̲,1 1 II‑ I 'I‑ I d‑1 ‑‑lIIIlIIll llll'l,.,lll"I'l'llll"'1'1111・IIIlIIIIlIIII'l ItJ llJl 1・llll・・1‑Ill‑・・ll・llll,ll・blll・lIIllal l‑ll・alllb・‑.1,llll,a.lb・It・It,‑・ll・ll・・・・l
O
20 40 60 80
1 oox [mm]
t=6.00‑7.25 ms
l 5 Il B 5 cUe b : O) d; (t 6 OO 7 25ms,P4‑1.5MPa,Case 4)
‑ E E
‑
)Bi‑axial stress loading End wall
30 20
10‑10
‑20
‑30
displacement
: '
l 11 It
・‑lllll,
lll・Wlllt
IlIII・・i
l‑ ,1, ‑・ l・ ,
If " 11 : "' a'l :・‑
t=7.50 ms t=7.75 ms t=8.00 ms t=8.25 ms t=8.50 ms t=8.75 ms
',・ aU at・ l・S ,la
O 20 40 60 80 1 OO
x [mm]t=7.50‑8.75ms
l 5.12 I !B 5i cUe e i o) d; (t 7 50 8 75ms,P4 1 5MPa Case 4)
Bi‑axial stress loading End wall of rubber
r‑l OO
ILI
COO 1 oO
80 60 40 20 O
P4 / P1 =2 ‑‑‑‑' P4 / P1 =1 5
(case I ) (Case 4)
P4 / P1 =5 P4 / P1 =20
(case 2) (case 5)
‑ ‑ ‑ ‑ ' ‑"‑"' P4 1 P1 =30
P4/ P1 =10: /¥. (case 3) (Case 6)
. ‑.,!;//// :;:1 ̲ s /' / ' ll / / ¥/ ¥
(f ! /"
!/ / / ¥ : l' ' f f f f
i/1:///f:':' f " ' " //j/ ' : i) :1bi . Ib / /
. f/ // . / "'
!. . I ‑d cc o""' '/ '$ ;'̲'̲̲̲ .̲1' l ,,, .
$O 2 6 8 4* t [ms]
l 5.13 o) l'y ec e , )p * ・.^・4
10
Bi‑axial stress loading
‑
caCL
er 12 10 8 6 4 2 o
t l
t
t
t
t
I'
I
I l
,
,
t
t l I I
I
I t
O
I RUbber
!1
l i
l """"'
i ‑‑‑‑'
.t,! i i "p 'l
l le
,i ILi
l li / . " i
. . ,‑ Il l l
. . , ,, l gll ̲ '̲ s" L ̲.̲̲a
:"I ' t
':1' ! . .L
:'il' l "t T :
fll "
. . .. . .
a"'/' '
l ,.,' ,
'l‑
Oase I ‑‑‑‑' Oase 4
Oase 2‑'‑" oase 5 Oase 3‑"‑"' case 6
It' ' L
i ' .t
ii '
lr'¥i
' i
! i / ¥P¥
/'X ¥ '
1 i
l. . l¥
/ : < ̲・; :: ":; ;t X:':¥
i /
l, . I L ' / /
/ f'/' /'/
'/i/1!/' ̲'/ ,/
" " : ' ?t [ms]
l 5.14 ) l ,y c e F : o)B * ・.^.4
10
Bl axlal stress loading End wall of rubber
)
1 OO
80
‑ 60
OOco 40 o
20 O
I
,
l
,
l
I I
i
,
t I
I
,
, I I I
, I ,
Rubber P / P =15
(Case 4)4 1
O I O 40 50 20
t [ms]30
l 5.15 q) ly f C e ; 10)B * .. 4 (P4=1.5MPa,Case 4)
Bi‑axial stress loading
‑
cr 12 10 8 6
4*
o
I
t
I
I
I
I
I
t
I
l
, ,
l
I
l.
,.
,
=1
t I
,
O
Rubber P /P =15
(Oase 4)
4 1
1 O 40 50 20
t [ms]30
l 5.16 q) l y+ j ec e F : )B* * '. '4 (P4 1.5MPa,Case 4)
L
−
o
4.5
4.0
3.5
3.0
2.5
2.0
1.5
1.0
◇
Rubber
◇
◇
◎
◇
◎
◇
051015202530
巳/R
図5.17ゴムの衝撃波管圧力比とDLF(動荷重係数)との関係
Nose(So)
iiiiiil難iiiill難ii難lii…iil難iiiiiiヨii蕪li…iiiiii難萎ii…ii
Poi nt・†P・int2P・int3 P・int4P・int5
・End Wa .
図5.18リーマン波の絶対速度のサンプリシグ点の概略図
Bi‑axial stress loading 1 20
1 OO
80
E̲ 60
O t 40
D
20 O
‑20
‑ Rubber
C ( Ur + c ) p4 / P1 =15 (case 4)
/ /1 "‑' !
' t';i ;} i¥L¥ :̲i:/ ,"/// " ;':
l ..:': L:
! ' ¥‑ ;'¥
l I "'//' " ll
" I ""// ' '$'b'. : :: '.¥
s i I I "I '1 1 " l)s . 'td::'/ t ' ; '
i: '" : ' : U POintl ‑‑‑‑' point4
r""""' ‑'‑" Point5
Point2‑"‑"' POint3 i/ . ff: L ¥¥i ::¥ I¥ '
'‑‑' "
'lpl/ .‑ ... ,
̲ItL :' '/‑ ‑ ‑ ‑ ‑
' ‑ ' ‑ !
¥"‑'.d:. i b ̲̲̲
" ¥; :b̲. ̲ .
‑,
O 2 4 8 10
t [ms]6
l 5.19 :,j o) EEl f; C l I e ,) 7/Y; q) L C,U.q) i * ..̲.4
5.3 気体衝撃波の衝突による柱状ゴムの一次元干渉解析,
及び実験との比較・
5.3.1 気体側の基礎方程式
気体衝撃波の衝突による柱状ゴムの一次元変動解析を行うに当たり,気体側の流 れは圧縮性の一次元非定常問題とし,流体の粘性は無視している.弾性体側の計算 条件,基礎式,初期条件などは柱状ゴムにステップ状の外力を印加した場合の計算に 使用したものと同じであり,気体側とゴム側は境界条件を介して同時に計算を行っ たものである.また,この解析はMazor et.a1.(9)のゴムに関する実験の条件に合わ せて行った.
気体側の基礎方程式の無次元化は弾性体側の無次元化と同じように気体の音速,
密度,長さを用いた.無次元化の各量は以下の通りである.
オα1 P
ρ9 ¢ 孟多二τ了,P二雇,ρ}ニk厚,ゴ二耳,砺1 σ.T ん9
¢=冴・7=rヌr,ん多二ρ91L1み9 但し,添字9は気体(gas)に関する量を示すもので,また
LTOα1 L1≡一.
CO
ここで,L1、は,衝撃波管の長さではないので注意しなければならない.
オイラー座標系で表された気体側の質量保存式は,
勉+∂(ρ9砺)一・ (5鋤
∂オ , ∂の
運動方程式は,
∂(ρ9砺)+∂(ρ9確+P)ニ。 〈5.5。)
∂オ ∂¢
エネルギー方程式は,
璽+∂(E+P)砺一。 (5.51)
∂オ ∂¢
P l E一胃+喜ρ9婿
である.ここで,Eは単位体積当たりの全エネルギー,7は気体の比熱比である.
オイラー座標系で表された気体側の質量保存式,運動方程式,エネルギー方程 式を弾性体側と同様な方法を用いて一般座標変換を行う.まず,ラグランジュ変数 4ん,・=ρ,(面一砺4オ)≧定義すると・質量と質量流束は以下のようになる・