見る/開く

by user

on 28 марта 2017

Category: Documents

>> Downloads: 2

views

Report

Comments

Description

Download 見る/開く

Transcript

見る/開く

佐賀大学農学業報
第24号
1
1∼48
電子計算機の利用による土の物理試験
データの計算処理について
生島芳雄・細山田健三・中島明
（土地改良学研究室）
Ont
h
eC
a
l
c
u
l
a
t
i
v
eTreatmentso
ft
h
e DataofP
h
y
s
i
c
a
l
T
e
s
t
sf
o
rS
o
i
lbyu
s
i
n
ganE
l
e
c
t
r
o
n
i
c
D
i
g
i
t
a
lComputer
YoshioIKusmMA,KenzoHosoYAMADAandAkiraNAKAJIMA
(
L
a
b
o
r
a
t
o
r
yofLandM
e
l
i
o
r
a
t
i
o
n
)
Summary
Thea
u
t
h
o
r
sprogrammedt
h
ec
a
l
c
u
l
a
t
i
o
n
so
ft
h
ed
a
t
ao
fp
h
y
s
i
c
a
lt
e
s
t
sf
o
rs
o
i
lby
u
s
i
n
gt
h
ee
l
e
c
t
r
o
n
i
cd
i
g
i
t
a
lcomputer“
HITAC5020E
．
”
I
t
e
m
sp
r
o
v
i
d
e
df
o
rthema
r
e:
(
1
)
(
2
)
(
3
)
(
4
)
(
5
)
(
6
)
(
7
)
(
8
)
i
f
i
cg
r
a
v
i
t
ybyJ
I
SA 1202
s
p巴c
s
p
e
c
i
f
i
cg
r
a
v
i
t
yandm
o
i
s
t
u
r
eo
fc
o
a
r
s
ea
g
g
r
e
g
a
t
ebyJ
I
SA 1110
waterc
o
n
t
e
n
tbyJ
I
SA 1203
b
u
l
kd
e
n
s
i
t
ybymercurymethod
b
u
l
kd
e
n
s
i
t
ybyformingmethod
b
u
l
kd
e
n
s
i
t
ybyp
a
r
a
f
f
i
nwaxp
a
i
n
t
i
n
gmethod
m
e
c
h
a
n
i
c
a
la
n
a
l
y
s
i
sbyJ
I
SA 1204
l
i
q
u
i
dl
i
m
i
tbyJ
I
SA 1205
p
l
a
s
t
i
cl
i
m
i
tbyJ
I
SA 1206
c
e
n
t
r
i
f
u
g
em
o
i
s
t
u
r
ec
o
n
t
e
n
tbyJ
I
SA 1207
自e
l
dm
o
i
s
t
u
r
ee
q
u
i
v
a
l
e
n
tbyJ
I
SA 1208
s
h
r
i
n
k
a
g
ec
o
e
妊1
c
i
e
n
tbyJ
I
SA 1209
Treatmentso
fa
l
lc
a
l
c
u
l
a
t
i
o
n
sa
r
ea
r
i
t
h
m
e
t
i
c
a
l
,s
ot
h
a
tt
h
e
ydon
o
tneeds
p
e
c
i
a
lo
p
e
r
a
t
i
o
n
.
E
s
p
e
c
i
a
l
l
y
,t
h
eprogrammeso
f(
5
)m
e
c
h
a
n
i
c
a
la
n
a
l
y
s
i
sa
r
el
o
n
g
,b
u
tt
h
et
i
m
eo
ft
h
e
,
t
r
e
a
t
m
e
n
t
sf
o
ra
b
o
u
t10s
a
m
p
l
e
si
sw
i
t
h
i
n1m
i
n
u
t
e
.
Theprogrammeso
fc
a
l
c
u
l
a
t
i
o
n
so
fm
e
c
h
a
n
i
c
a
lt
e
s
t
sf
o
rs
o
i
la
r
enowunberi
n
v
e
s
t
i
g
a
t
i
o
n
.
1 まえがき
ほとんどの土木工事，耕地の改良等に土の物理試験あるいは力学試験を要することは今さら述
べるまでもない乙とであるが，ここで注意すべき乙とは土 l
とは均一性の保証はなく，大抵の場合
1
2
佐賀大学農学業報第 24
号
ある方向に層をなして堆積していることである．たとえば斜面の安定を計算するのに，予想さ
れるスベリ面 1カ所の土を採取し，単位体積重量，粘着カ，内部マサツ角等を試験して代表値
とする乙とは許されないのであって，多数カ所の土を採取して試験することが計算の精度を高め
る結果 l
となる．数多くの試験用土の採取と同時に，じん速な試験とそのデータの計算処理が必要
となるゆえんである．筆者等はこのような意味から電子計算機の利用を考えた．本報は土質試験
AC5020E を使って，
のうち，物理試験のみをとりあげ，東大大型計算機センターの HIT
FORTRANIVKよるプログラムを組んで計算した結果で−ある．
J
I
. 試験項目
プログラムは日本工業規格 (
J
I
S
)I
C
:
:基づいて作られた土質工学会の土質試験用データシートに
従って組んだ．以下プログラム化した物理試験項目を列挙するとつぎの 8項目になる．
(
I
) 土粒子の比重試験
J
I
S・
A
:1202
(
2
) 組骨材の比重および吸水量試験 J
I
SA1
1
1
0
I
SA1203
(
3）土の含水量試験 J
(
4
) 単位体積重量試験
(
5
) 粒度試験 J
I
SA 1204
;2000μ フノレイ残留部分
¥
2
0
0
0μ フルイ通過部分
(
6
)
波性限界 J
I
SA 1
2
0
5
ず·~in 限界
JIS
A1
2
0
6
(
7
) 遠心含水当 }
i
t試験
J
I
SA1207
現場合水当量試験 J
I
SA1208
(
8
) 収縮常数試験
J
I
SA 1
2
0
9
m.
温度に対する水の比重と補正係数＊ (F）の関係
(
1
）土粒子の比重試験と（ 5
）校度試験を行なう場合，前者では水の比重，
後おーでは水の比重と補
正係数＊ (
F）が関係するので，種々の温度 l
乙対する水の比重と補正係数（F）の関係式を求めて
,
2を得た．同
おく必要がある．電子計算機によって，それぞれに最小二乗法を通用して，表− 1
表の計算 I
,I
Iの差は既存の水の比重および補正係数（F）の温度関係 r
H
1総の読みとりの違いによ
るものであって，中等誤差，確率誤差の小さい計算 Iの方をとると
(
i
) 水の比重（
Gr
）と温度（X°C
）の関係は
Gr=-0.00000501X2- .
o000000006x+1
.0
0
0
1
7
9
• (
i
i
) 補正係数（ F）と温度（X°C
）の関係は
F=O.0
0
0
0
0
4
3
7X2+0
.0
0
0
0
0
0
0
0
2X - 0
.0
0
0
9
0
0
9
6
I
V
. フローチャー卜，プログラムおよび結果
(
1
) 土粒子の比重試験
＊水の比重ζ
l関する温度の補正係数（ K）ではない．
1
3
生島・細山・中島：電子計算機の利用による土の物理試験データ
表− I 温度と水の比重および F値
温
*
度
。
C
の
｜計算
4
5
6
7
8
9
1
0
1
1
1
2
1
3
1
4
1
5
1
6
1
7
1
8
1
9
2
0
2
1
2
2
2
3
2
4
2
5
2
6
2
7
比
I ｜計算
I
.0
0
0
0
9
9
I
.0
0
0
0
5
4
0
.
9
9
9
9
9
9
0
.
9
9
9
9
3
3
0
.
9
9
9
8
5
8
0
.9
9
9
7
7
3
0
.
9
9
9
6
7
8
0
.
9
9
9
5
7
3
0
.
9
9
9
4
5
7
0
.
9
9
9
3
3
2
0
.
9
9
9
1
9
7
0
.
9
9
9
0
5
1
0
.
9
9
8
8
9
6
0
.
9
9
8
7
3
1
0
.
9
9
8
5
5
5
0
.
9
9
8
3
7
0
0
.
9
9
8
1
7
4
0
.
9
9
7
9
6
9
0
.9
9
7
7
5
3
0
.
9
9
7
5
2
8
0
.
9
9
7
2
9
2
0
.
9
9
7
0
4
7
0
.
9
9
6
7
9
1
0
.
9
9
6
5
2
6
0
.
9
9
6
2
5
0
0
.
9
9
5
9
6
4
0
.
9
9
5
6
6
9
｜計算
［計算
I
.0
0
1
3
0
6
-0.0006
-0,0006
-0.0006
-0.0
0
0
6
-0.0006
-0.0
0
0
5
-0.0
0
0
5
0
.
0
0
0
4
-0.0
0
0
3
-0.0
0
0
2
-0.0
0
0
1
0
.
0
0
0
0
+o.0
0
0
1
+o.0
0
0
3
+o.0004
+o.0006
+0.0008
+o.0010
+o.0
0
1
2
十0
.
0
0
1
4
+0.0016
+o.0018
+o.0020
+o.0
0
2
3
+o.0
0
2
5
+o.0028
+o.0
0
3
1
I
.0
0
1
2
0
6
I
.0
0
1
0
8
3
I
.0
0
0
9
3
7
I
.0
0
0
7
7
0
I
.0
0
0
5
8
0
I
.0
0
0
3
6
8
I
.0
0
0
1
3
4
0
.
9
9
9
8
7
7
0
.
9
9
9
5
9
8
0
.
9
9
9
2
9
7
0
.
9
9
8
9
7
3
0
.
9
9
8
6
2
7
0
.
9
9
8
2
5
9
0
.
9
9
7
8
6
9
0
.
9
9
7
4
5
6
0
.
9
9
7
0
2
1
0
.
9
9
6
5
6
4
0
.
9
9
6
0
8
4
0
.
9
9
5
5
8
2
0
.
9
9
5
0
5
8
0
.
9
9
4
5
1
2
0
.
9
9
3
9
4
3
0
.
9
9
3
3
5
2
0
.
9
9
2
7
3
9
0
.
9
9
2
1
0
3
0
.
9
9
1
4
4
5
0
.
0
0
1
2
6
7
し／／｜
h
戸 J
0
.
0
0
0
0
6
8
QUQd
レ／｜
nO
0
.
0
0
1
8
9
2
nunu
nunU
0
.
0
0
0
1
0
1
-0.0
0
0
8
3
-0.0
0
0
7
9
-0.0
0
0
7
4
-0.0
0
0
6
9
-0.0
0
0
6
2
-0.0
0
0
5
5
-0.0
0
0
4
6
-0.0
0
0
3
7
-0.0
0
0
2
7
-0.0
0
0
1
6
-0.0
0
0
0
5
十0
.
0
0
0
0
8
+0.00022
+o.0
0
0
3
5
+o.o
o
o
s
1
十0
.
0
0
0
6
8
+0.00085
十0
.0
0
1
0
2
+o.0
0
1
2
1
+0.00141
+o.0
0
1
6
1
+o.0
0
1
3
3
+o.0
0
2
0
5
+o.0
0
2
2
8
+o.0
0
2
5
2
+o.0
0
2
7
7
+o.0
0
3
0
3
nunU
レ
ペ
レ／／｜
E
F 値
∞∞
等差一率差
中誤一確誤
1
.
0
0
0
0
0
0
0
.
9
9
9
9
9
2
0
.
9
9
9
9
6
8
0
.
9
9
9
9
3
0
0
.
9
9
9
8
7
7
0
.
9
9
9
8
0
9
0
.
9
9
9
7
2
8
0
.
9
9
9
6
3
4
0
.
9
9
9
5
2
6
0
.
9
9
9
4
0
6
0
.
9
9
9
2
7
3
0
.
9
9
9
1
2
9
0
.
9
9
8
9
7
2
0
.
9
9
8
8
0
4
0
.
9
9
8
6
2
5
0
.
9
9
8
4
3
5
0
.
9
9
8
2
3
4
0
.
9
9
8
0
2
2
0
.
9
9
7
8
0
0
0
.
9
9
7
5
6
8
0
.
9
9
7
3
2
7
0
.
9
9
7
0
7
5
0
.9
9
6
8
1
4
0
.
9
9
6
5
4
4
0
.
9
9
6
2
6
4
2
8
0
.
9
9
5
9
7
6
2
9
0
.
9
9
5
6
7
8
3
0
一一一
補正係数
重
H
-0.0
0
0
8
4
-0.0
0
0
8
0
-0.0
0
0
7
5
-0.0
0
0
7
0
-0.0
0
0
6
3
-0.0
0
0
5
6
-0.0
0
0
4
7
-0.0
0
0
3
8
-0.0
0
0
2
8
-0.00017
-0.0
0
0
0
6
+o.00001
+o.0
0
0
2
0
+o.0
0
0
3
5
+o.0
0
0
5
0
+o.0
0
0
6
6
+o.0
0
0
8
3
+0.00101
+0.00120
+o.0
0
1
3
9
+o.oo凶O
+o.0
0
1
8
1
十0
.
0
0
2
0
3
+o.0
0
2
2
6
+o.0
0
2
s
o
+o.0
0
2
7
5
+o.0
0
3
0
1
0
.
0
0
0
0
9
1
0
.
0
0
0
0
6
2
表 2 水の比重および F値の計算式の係数
水の比重の場合
計
算
計
算
補正係数 F の場合
I
I
計
算
計
算
I
I
a
-0.0
0
0
0
0
5
0
1
-0.00001114
0
.
0
0
0
0
0
4
3
7
0
.
0
0
0
0
0
4
3
5
b
-0
.0
0
0
0
0
0
0
0
6
-0.0
0
0
0
0
0
3
5
0
.
0
0
0
0
0
0
0
0
2
0
.
0
0
0
0
0
0
0
0
2
c
!
.0
0
0
1
7
9
I
.0
0
1
4
8
6
-0.0
0
0
9
0
0
9
6
-0.0
0
0
9
0
8
6
6
佐賀大学農学業報第24
号
1
4
(
a
) 記号の説明
N
：試料の数
M
ANl
羽T
Fl
：ピクノメーターの番号
WAI
Tl
WBl
T2
AN2
WD
WT
WO
Gl
：測定数
：ピクノメーターの重量
：（ピクノメーター＋蒸留水）の重量
:WAIを測った時の温度
：（ピクノメーター＋炉乾燥（または湿潤）土＋蒸留水）重量
:WBlを測った時の内容物の温度
：炉乾燥土秤量の容器番号
：（容器十炉乾燥土）重量
：容器重重
：ピクノメーターに入れた土の炉乾燥（または湿潤土）重量
羽TAF
：
〔T2における水の比重（G
γ）］／〔Tlにおける水の比重（G.
／
）
〕
:WAI一明T
Fl
G2
:Gl× WAF
:G2+WFl
羽TOAB :W O
十（羽TAl-WBl)
WA2
GSl
G3
:T21とおける比重＝ WO/WOAB
：補正係数 K
GS2
GS3
:15°Cにおけるよヒ lH=G3×GSl
:.
l
'
G
S
2
GR
:l5°Cにおける比重の平均航
(
h
)
フローチャート
11
~コ
↓
／
Nを読むつ
↓
N 恒l
繰り返す
／
M を必む
I
GS3=0.0
／
I
1
1
1川目繰り返す
A
:
'
¥I
,¥
¥
"
F
I
, W AI
, TI, ¥
¥
"
B1
, T2
,
A
:
'
¥2
,¥
¥
"
D
, WT を＂~＂む
↓
.
.
¥
'
.
'
¥
I
,¥
¥
"
F
I
, W A1
, TI,¥
¥
"
B1
,T2
,
A
:
'
¥2＼
＼
・D
, W Tを印刷する
↓
＼
＼
・0
,G 1
,¥
¥
"
A
F
,G2
,W
A2
,¥
¥
"
O
A
B
,
GSi
,G3
, GS2
, GS3を計算する
↓
＼
＼
・0
,G 1
,¥
¥
"
A
F
,G2
, W A2
,＼
＼
・O
AB,
GS1,G3,GS2，を印刷する
↓
1
5
生島・細山・中島：電子計算機の利用による土の物理試験データ
I GRを計算する
L GRを印刷する
｜
／
C豆二〉
(
c
) プログラムおよび結果
この場合のプログラムは Program1
1であって結果は表 3
1
<
:示す．
IProgram
scc
1
0
0
1
1
0
1
2
0
1
3
0
1
5
0
1-1 j
HARP
PHISICALTESTFORSOIL1
1
OPERATIONSONSPECIFICGRAVITYTESTBYJISA 1
2
0
2
WRITE(
6
,1
0
0
)
FORMAT(
1H 1
,5X, 1
3HSOILTEST 1-1/1HO,5X,3
4HOPERATION
ONSPECIFICGRAVITYTEST)
READ(
5
,1
1
0
)N
FORMAT(!4)
DOIOI=1
,N
READ(
5
,1
2
0
)M
FORMAT(
I4
)
GS3=0
.0
DO20]=l,M
,T 1
,W B1
,T2
,AN2
,WD,W T
READ(
5
,1
3
0
)AN1
,WF1
,W A1
FORMAT(
9
F8
.0
)
WRITE(6,1
5
0
)AN!,WFI,WAT,Tl,WBl,T2,AN2,WD,W T
FORMAT(!HO,3HN0.E
l
3
,5
,3X,3HWF=El5.7
,3X,4HWAl=El5.7
,
3X,3HTl=El5.7
,3X,4HWBl=El5.7/!Hb, 3X,3HT2=El5.7
,1
1X,
3HN0.E
l
3
.5
,3X,3HWD=El5.7
,3X,3HWT=E1
5
.7
)
WO=WD-WT
（
ー0
.5
0
1
1E-05）
キ T2**2+ （
ーo
.6E-OB)*T2+o.oool79E+o1)/
G 1=（
0
.5011Eー05)*TI**2十（ー 0.6Eー08)*Tl十0
.1000179E+Ol)
WAF=WAト WFI
G2=GlホWAF
WA2=G2+WFI
W0AB=WA2+W0-WBI
GSl=WO/WOAB
.6E-08
）
キ T2十1
.000179E+01)I
(
(
G3=(
( 0
.5
0
1
1Eー05hT2料 2＋（ー 0
0
.5
0
1
1E 0
5）
キ1
5
.0**2＋（ー 0
.6E-08)*1
5
.0+1
.000179E+Ol)
GS2=GSl*G3
GS3=GS3+GS2
WRITE(6,160)W0,Gl,WAF,G2,WA2,WOAB,GSI,G3,GS2
FORMAT(lHb, 3X, 3HW0=El5.7
, -3X, 7HGR/GRl=El5.7
, 3X,
7HWAl-WF=El5.7
,3X, 16HGR/GR1*(WAl-WF)=El5.7/lHb,3X,
3HWA=El5.7,3X, llHW0+(WA-WB)=El5.7,3X, 16HW0/(WO+
(WA-WB))=E1
5
.7
,
3
X
,2HK=E1
5
.7/1Hb,3X,3HGR=E1
5
.7
)
CONTINUE
GR=GS3/FL0AT(M)
WRITE(
6
,1
7
0
)GR
FORMAT(
IHb,2
5X,l4HHEIKINCHIGR=E1
5
.7
)
<
<
-
1
6
0
2
0
1
7
0
STI
STI
ST1
S TI
ST1
ST1
ST1
ST1
STI
STI
ST1
STI
STI
ST1
ST1
ST1
ST1
STI
ST1
ST1
ST1
ST1
S TI
S TI
S T1
S TI
ST1
S T1
S T1
S T1
ST1
S T1
S T1
ST1
7
1
72
7
3
7
4
7
5
7
6
7
8
7
9
80
8
1
8
2
8
3
84
8
5
8
8
8
9
9
0
9
1
9
2
9
3
94
9
5
9
6
9
8
1
0
0
1
0
3
1
0
6
1
0
7
1
0
8
1
1
3
1
1
8
1
1
9
1
2
0
1
2
1
ST1
S TI
S TI
ST1
1
2
2
1
2
3
1
2
4
1
2
5
1
6
4
号
佐賀大学農学蒙報第 2
STI
S TI
STI
STI
STI
T -I
1
0 CONTINUE
STOP
END
$
DATA
2
2
24.257
2
1
.0
7
91
6
4
29.0
2
5
.7
6
5
7
9
.
6
3
2
1
1
.8
1
1
.8
97.967
96.430
1
1
.7
1
1
.7
9.2
9.2
1
0
0
.
4
9
2
1
0
2
.
0
8
4
1
2
.
0
1
2
.
0
2
5.0
20.0
20.950
23.890
8
1
.284
8
2
.
5
6
1
表
3
.
3
9
1 24.257 T I
2
1
.0 5
29.0 5
1
.7
4
5 2
5
.7
6
5 T -I
y -I
5.0 53.212 20.950 y - I
20.0 5
6
.1
1
2 23.890 y -I
3 真比重の計算結果
容器番号
No. 2
1
N
o
. 29
No. 5
No. 20
GS2 の値
2.8106
2.8190
2.4749
2
.5410
平均値 GR
2.5079
2.8148
(
2
) 粗骨材の比重および吸水量試験
(
a
) 記号の説明
N
：試料の数
GRl
:~GR2
AMI
B
： ~AM2
：表面乾燥飽和状態の試料の空気中軍地
c
：試料の水中重量
A
：乾燥後の試料の重量
BC
:B-C
GR2
：粗骨材の比重
BA
:B-A
AM2
：粗骨材の吸水量
GR
AM
：平均比重
：平均吸水量
B
、
b
，
，
フローチャー卜
2-1
Nを読む
t
N同繰り返す
。
GRl=0.0
A M1=0.
↓
2回繰返す
一一了一−
B,C,Aを読む
1
2
6
1
2
7
1
2
8
1
2
9
1
3
0
2
3
2
3
生島・細山・中島：電子計算機の利用による土の物理試験データ
I B,C,Aを印刷する
1
7
／
B
C
,GR2
,B
A
,A恥12
を計算する
↓
B
C
,GR2
,B
A
,AM2
を印刷する
一
一
一-r
IGRl,AMlを計算する｜
T
IGR,AMを計算する｜
IGR,AMを印刷する／
L一
一
一
一
C二Eコ
二
(c) プログラムおよび結果
乙の場合のプログラムは P
rogram2-1であって，結果は表− 4!
C示す．
事
c
c
HARP
S T1 2
0
0
PHSICALTESTFORSOIL2
1
OPERATIONSONSPECIFICGRAVITYANDMOISTUREOFCOARSE S T1 2
0
1
AGGREGATEBYJISA 1
1
1
0
S T1 2
0
2
WRITE(
6
,2
0
0
)
2
0
0 FORMAT(IH!
,5X,13HS0ILTEST2
1
/
1HO,5X,65HOPERATIONON S TI 2
0
3
GRAVITYANDMOISTUREOFFINEAGGREGATEBYJISA 1
2
0
3
)
ST 1 2
0
4
S T1 2
0
5
READ(
5
,2
1
0
)N
S T1 2
0
6
210 FORMAT(
!
4
)
S TI 2
0
7
DO!OI=l,N
S TI 2
0
8
GRl=O.O
S T1 2
0
9
AMl=O.O
S TI 210
D020J=l,2
S T1 2
1
1
READ(
5
,2
2
0
)B
,C
,A
S T1 2
1
2
2
2
0 FORMAT(
3F1
2
.0
)
S T1 2
1
3
WRITE(
6
,2
3
0
)B
,C
,A
S T1 214
2
3
0 FORMAT(
1HO,2HB=E1
5
.7
,3X,2HC=E1
5
.7
,3X,2HA=E1
5
.7
)
S T1 2
1
5
BC=B-C
S T1 2
1
6
GR2=B/(B-C)
S T1 2
1
7
BA=B-A
S T1 2
1
8
AM2=(B-A)/A*100
S T1 2
1
9
WRITE(
6
,2
4
0
)BC,GR2
,BA,A M2
2
0
2
4
0 FORMAT(lH,3X,4HB-C=El5.7,3X,3HGR=E15.7,3X,4HB-A= S TI 2
El5.7,3X,3HMW=El5.7
)
S T1 2
2
1
GRl=GRI+B/(B-C)
S T1 2
2
2
A Ml=AMl+(B-A)/A*IOO.
S T1 2
2
3
2
0 CONTINUE
S T1 224
GR=GRl/2.
S T1 2
2
5
AM=AMl/2.
S T1 2
2
6
WRITE(
6
,2
5
0
)GR,A M
。
。
。
1
8
佐賀大学農学業報第 2
4
号
250 FORMAT(IH,5
2
7
:
X
,IOHHEIKINGR=El5.7
,5X,IOHHEIKINM W=El5.7
) S TI 2
2
8
ST I 2
1
0 CONTINUE
STOP
S TI 229
END
ST! 230
DATA
S TI 2
3
1
事
STI 232
3
3
2
7
2
6
.
0
4
4
1
8
.
5
4
3
2
9
.
5
S TI 2
2
7
3
0
.
4
3
2
3
.
4
4
1
2
.
3
STI 234
。
。
表− 4 真比重と吸水量の計算結果
料
試
比
GR
重
GR の平
均
吸
水
品1W
量
%
2
.
0
5
5
7
2
.6106
2
.
6
1
6
7
2
.
2
.
6
2
2
8
(3) 土の含水量試験
(a) 記号の説明
N
：試料数
M
：測定数
K
：容器番号
羽T
l
：湿潤重量
DW
：乾燥重量
TW
：容器重量
羽T
2
：水分重量＝ Wl-DW
ws
w
：土粒子重量＝ DW-TW
：含水比
AMW2 :2W
AM
羽f
：平均含水比
r
Lu
、
‘
，
フローチャート
31
の平
%
2
.
0
6
0
7
2
.
0
6
5
7
均
1
9
生島・細山・中島：電子計算機の利用による土の物理試験データ
IW2,ws,w を計算する｜
/w2,ws,w を印刷する／
IAMW2を計算する
｜
IAMWを計算する｜
I AMWを印刷する／
C二互コ
(c) プログラムおよび結果
乙の場合のプログラムは Program3 1であって，結果は表−5r
e示す．
$
c
c
1
0
0
1
1
0
1
2
0
1
3
0
1
4
0
1
5
0
2
0
1
6
0
1
0
事
HARP
PHISICALTESTFORSOIL3
1
OPERATIONONWATERCONTENTTESTBYJISA 1
2
0
3
WRITE(
6
,1
0
0
)
FORMAT(lHI,13HSOILTEST3-1/IHO, 3!HOPERATIONONWATER
CONTENTTEST)
READ(
5
,1
1
0
)N
FOR
！
＼
在AT(
1
4
)
DO101=1,N
AMW2=0.0
READ(
5
,1
2
0
)M
F0RMAT(l4)
0020J=l,M
READ(
5
,1
3
0
)K,W 1
,DW,T W
FORMAT(
1
4
,3F1
2
.
0
)
WRITE(
6
,1
4
0
)K,W 1
,DW,T W
F0RMAT(lH0,3HN0.14,5X,3HWW=El5.7,5X,3HDW=El5.7,5X,
3HTW=El5.7
)
W2=Wl-DW
WS=DW-TW
）
キ1
0
0
.
0
W=((W1-DW)/(DW T W）
WRITE(
6
,1
5
0
)W 2
,WS,W
FORMAT(lH,5X,3HWW=El5.7
,5X,3HWS=El5.7
,5X,2HW=El3.5,
BHPERCENT)
AMW2=AMW2+(((W1-DW)/(DW TW))*lOO.0
)
CONTINUE
AMW=AMW2/FL0AT(M)
WRITE(
6
,1
6
0
)AMW
FORMAT(
1H ,5X,IOHHEIKINCHI=E1
3
.5
,8HPERCENT)
CONTINUE
STOP
END
DATA
ST2
ST2
ST2
S T2
ST2
ST2
ST2
ST2
ST2
S T2
S T2
S T2
S T2
S T2
S T2
ST2
2
3
4
5
6
7
8
9
1
0
1
1
1
2
1
3
1
4
1
5
1
6
S T2
S T2
S T2
S T2
S T2
1
7
1
8
1
9
20
2
1
S T2
S T2
S T2
S T2
S T2
S T2
S T2
S T2
S T2
2
2
2
3
2
4
2
5
2
6
2
7
2
8
2
9
30
号
20
佐賀大学農学業報第2
4
6
5
1
1
23
59
78
43
5
47
229
1
6
3
2
3
1
33
5
44
1
3
95
50
62
5
1
9
1
3
1
5
6
1
1
0
47
5
1
3
47
1
1
1
43
1
9
5
1
3
1
1
1
1
8
8
1
1
0
220
26.574
27.347
24.674
27.270
29.832
2
1
.0
3
7
2
1
.608
1
9
.1
0
6
2
1
.229
22.845
1
0
.732
1
1
.820
9.848
1
1
.027
1
0
.
3
4
9
29.654
28.267
22.389
24.222
2
5
.7
6
2
25.275
24.572
1
9
.1
8
2
2
1
.1
5
0
2
1
.7
6
7
1
0
.9
5
1
1
2
.
9
8
2
1
0
.678
13.054
1
1
.278
22.520
22.534
2
7
.740
20.348
20.960
1
9
.7
4
5
1
9
.1
8
5
2
3
.1
7
2
1
7
.
5
7
3
1
8
.1
0
7
1
2
.
4
5
1
1
0
.878
1
1
.579
1
0
.1
2
2
1
0
.784
2
7
.080
28.067
26.279
26.114
2
6
.1
1
1
20.809
2
1
.622
1
9
.
9
1
5
1
9
.
5
4
9
2
0
.1
8
2
1
0
.1
1
4
1
0
.873
9.918
1
0
.1
4
4
1
0
.950
1
8
.
5
7
6
2
0
.410
20.282
2
1
.958
2
1
.324
1
5
.7
6
2
1
6
.
8
7
9
1
6
.
5
7
3
1
7
.7
8
9
1
7
.
1
2
8
1
0
.873
10.950
9.937
1
0
.
3
4
8
10.114
25.277
27.578
28.567
24.689
25.624
1
9
.
9
2
3
20.582
2
2
.7
4
5
1
9
.
0
7
1
1
9
.
4
6
9
1
0
.
8
7
3
9.937
1
3
.
3
2
9
10.144
9.981
S T2 3
1
T -I
T-I 2
T -I 3
T I 4
T I 5
T-1 6
T-I 7
T -I 8
T I 9
T -I I
O
T -I 1
1
T I 1
2
T 1 1
3
T -1 1
4
T 1 1
5
T -I 1
6
T -I 1
7
T -1 1
8
y -I
y -1 2
y -1 3
y 1 4
y -l 5
y -l 6
y -1 7
y -1 8
y -1 9
y -l 1
0
y -1 1
1
y -1 1
2
y -1 1
3
y l 1
4
y -l 1
5
y l 1
6
y l 1
7
y -1 1
8
表− 5 含水比の計算結果
業
｝
試
I
I
:
f
I
.
ι
器
番
号｜
一
含水
比
w
平均含水比 AMW
一一一一一
1
1
5
3
.731%
23
58.633
59
60.143
78
59.214
43
55.914
47
3
0
.
5
7
1
229
3
1
.8
8
1
1
6
3
3
7
.734
2
3
1
37.945
33
38.088
44
38.045
1
3
40.315
5
2
.527%
35.244
2
1
生島・細山・中島：電子計算機の利用による土の物理試験データ
9
5
i
l
l
I
V
v
V
I
3
8
.
9
2
7
50
62
3
9
.
4
0
3
3
7
.9
1
0
3
8
.
9
5
9
1
9
1
3
1
5
6
1
1
0
4
7
5
8
.
6
3
5
5
9
.
9
5
9
6
3
.
6
5
9
6
9
.
8
0
3
6
4
.
2
2
2
6
3
.
2
5
6
1
3
4
7
1
1
1
4
3
1
9
5
7
.
5
5
8
5
9
.
5
5
5
5
6
.7
4
2
5
6
.
0
2
7
5
9
.
8
2
3
5
7
.
9
4
1
1
3
1
1
1
1
8
8
1
1
0
220
5
9
.1
6
0
6
5
.7
2
1
6
1
.8
3
1
6
2
.
9
3
3
6
4
.
8
7
1
6
2
.
9
0
3
(4) 単位体積重量試験
水銀法の場合
(a) 記号の説明
N
,M,K
,Wl,TW,DW,W2,WS,A M Wは（3
）の説明ζ
l同じ．
WP
：含水比
AMWl :2WP
VO
：：ガラス容器の体積
Vl
W3
GS
：容器に残った水銀の体積
：供試体の重量
：供試体の比重
v
：供試体の体積
GMT
：湿潤単位体積重量
GMD
：乾燥単位体積重量
E
：間隙比
SR
：飽和度
(
b
) フローチャート 4-1
C二Eコ
／
I
I
Nを読む
／
N回繰り返す
I
M
ふむ／
｜
｜
22
佐賀大学農学業報第24号
L
AMWl=O.O
I
M 回繰り返す
IK,w 1
,DW,TWを読む／
IK,Wl,DW,TWを印刷する／
1
w
2
,WS,WPを計算する｜
/w2,WS,WPを印刷する／
IAMWl を計算する I
T
IAMW を計算する｜
IAMW‘を印刷する／
I
VO,v1
,w3
, GSを読む／
IVO,Vl,W 3，回を印刷する／
V
,GMT,GMD,E
, SR，を計
算する
V
,GMT, GMD, E
, SR を印
刷する
(c) プログラムおよび結果
C二Eコ
この場合のプログラムは Program4-1であって，結果は表− 61
1
:示す．
Program 4-1
$
c
c
1
0
0
1
1
0
1
2
0
1
3
0
1
4
0
HARP
PHISICALTESTFORSOIL4
1
OPERATIONSONBULKDENSITYBYMERCURYMETHOD
WRITE(
6
,1
0
0
)
FORMAT(lHl, 13HS0ILTEST4-1/!HO, 43H0PERATIONONBULK
DENSITYBYMERCURYMETHOD)
READ(
5
,1
1
0
)N
F0RMAT(14)
DO1
0l=l,N
READ(
5
,1
2
0
)M
FORMAT(
I4
)
AMWl=O.O
DO20]=1,M
READ(
5
,1
3
0
)K,W 1
,DW,T W
FORMAT(
I4
,3F1
2
.0
)
WRITE(6,1
4
0
)K,W 1
,mへ
んTW
FORMAT(
!H0
, 3HNO. I
4
, 5X,3HV¥
パIV=El5.7,5X, 3HDW=El5.7,
,5X,3HTW=El5.7
)
W2=Wl-DW
ST3
ST3 2
ST3 3
S T3 4
ST3 5
ST3 6
ST3 7
ST3 8
ST3 9
ST3 1
0
1
ST3 1
2
ST3 1
3
ST3 1
4
ST3 1
5
ST3 1
6
ST3 1
ST3
1
7
2
3
生島・細山・中島：電子計算機の利用による土の物理試験データ
WS=DW-TW
WP=((W1-DW)/(DW-TW))*lOO.0
WRITE(
6
,1
5
0
)W 2
,WS,W D
1
5
0 FORMAT(lH, 5X, 3HWW=El5.7, 5X, 3HWS=El5.7,5X,3HWP=
E1
3
.5
,8HPERCENT)
A M W1=AMW1 +(((W1-DW)/(DW-TW))*lOO.0
)
2
0 CONTINUE
AMW=AMWl/FLOAT(M)
WRITE(
6
,1
6
0
)AMW
1
6
0 FORMAT(lHb,5X,lOHHEIKNNCHI=E1
3
.5
,8HPERCENT)
READ(
5
,1
7
0
)VO,V 1
,W 3
,GS
1
7
0 FORMAT(4F1
2
.
0
)
WRITE(
6
,1
8
0
)VO,V1
,W 3
,GS
1
8
0 FORMAT(lHO, 3X, 3HV0=El3.5, 3X, 3HVl=El3.5, 3X, 2HW=
E1
5
.7
,3X,3HGS=E1
3
.5
)
V=VO-Vi
GMT=W3/(VO-VI
)
GMD=(W3/(VO V 1
)
)
/
(
1
.0+AMW/1
0
0
.0
)
キ
（GS/(W3/(VO VI）
）
）
ーI
.0
E=(l.O+AMW/
1
0
0
.0）
）
キ
（ GS/(W3/(VO-V 1
）
）
）
ーI
.0
)
SR=AMW*GS/((1.O+AMW/100.0
WRITE(
6
,1
9
0
)V,GMT,GMD,E
,SR
1
9
0 FORMAT(!HO,3X,2HV=El3.5
,3X,!OHGAMMERT=E1
3
.
5
,3X, 1
0
HGAMMERD=E1
3
.5
,3X,2HE=E1
3
.5
,3X,3HSR=E1
3
.5
)
1
0 CONTINUE
STOP
END
DATA
$
3
1
1
1
2
1
3
7
6
.
4
1
1
8
.
9
4
3
1
8
.
0
0
3
1
9
.
4
8
0
6
6
.
8
0
1
6
.
8
9
0
1
6
.
0
7
0
1
7
.
3
5
0
1
5
.
0
8
8
1
1
.6
0
1
1
1
.2
3
3
1
1
.9
2
2
2
.
6
4
ST3 1
8
ST3 1
9
ST3 2
0
ST3 2
1
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
2
2
2
3
24
2
5
26
2
7
2
8
2
9
3
0
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
A -1
A -I
A -I
A-1
A -1
3
1
3
2
3
3
3
4
3
5
3
6
3
7
3
8
3
9
40
4
1
4
2
4
3
2
3
4
5
表− 6 含水比，湿潤密度，乾燥密度，問ゲキ比，飽和度の計算結果（水銀法）
i平均含水比
試
AMW
l
湿潤密度｜乾燥密度！間ゲキ比
GMT
GMD
E
I
.5
7
0
0
!
.1
2
6
8
!
.3
4
3
飽
和
度
SR
7
7
.3
3
3
%
成形法の場合
(a) 記号の説明
N,M,K,Wl,W2,DW, TW, W2, WS,AMW は（3
)I
L
,WP, AMWl, W3,V,
,SRは水銀法の場合の説明に同じ
GS,GMD,GMT,E
AL
：供試体高さ
DT
：上部直径（平均）
DC
：中央部直径（平均）
DB
：下部直径（平均）
DAV
：全平均直径
A
：断面積
24
佐賀大学農学業報第 24
号
E
BEJ
、
、
b
，
，
フローチャート
4-2
こ
手
コ
り返す
読む
=0.0
A
L
,DT,DC,D
B
,W 3
, GS
を読む
↓
A
L
, DT, DC,D
B
,W3
, GS
を印刷する
DAV,A
,V
,GMT,GMD,E
,
SRを計算する
土
DAV,A
,V
,GMT,GMD,E
,
SRを印刷する
土
C三二〉
(c) プログラムおよび結果
と示す．
この場合のプログラムは Program4-2であって，結果は表一 71
Program 4-2
in403A
AFhd
仙
ST3
ST3
ST3
ST3
ST3
’
55555
scc
HARP
PHISICALTESTFORSOIL4
2
OPERATIONSO NBULKDENSITYBYFORMINGMETHOD
WRITE(
6
,2
0
0
)
200 FORMAT(
IH 1
,13HS0ILTEST4
2
/IHO, 43HOPERATION ON BULK
DENSITYBYMOLDINGMETHOD)
2
5
生島・細山・中島：電子計算機の利用による土の物理試験データ
、
1
9
.
8
3
3
2
6
.
9
9
0
2
5
.
8
8
4
3
.7
6
1
6
.
0
6
5
2
1
.3
5
5
2
0
.
4
3
0
3
.
8
6
9
.
9
3
9
1
1
.9
7
8
9
.
9
8
1
3
.
9
0
1
5
0
.
3
5
2
.
4
5
2
6
7
6
8
6
9
7
0
7
1
72
7
3
7
4
7
5
7
6
7
7
7
8
7
9
80
8
1
8
2
8
3
84
8
5
86
8
7
’
2
.4
5
2
3
3
4
3
.1
5
。
。
。
。
。
5
.
4
5
の
3
5
.
3
7
0
a
cdcaoDQdpD
5
.
2
8
TTTTT
1
0
.
6
7
3
1
1
1
2
3
7
2
2
0
8.84
5
6
5
7
5
8
5
9
60
6
1
62
6
3
6
4
6
5
6
6
oonynu A q L
oonmununwun3
READ(
5
,2
1
0
)N
S T3
210 FORMAT(
I
4
)
S T3
S T3
DO30I=l,N
S T3
READ(
5
,2
2
0
)M
S T3
2
2
0 FORMAT(
I
4
)
S T3
AMWI=O.O
S T3
0040J=l,M
S T3
READ(
5
,2
3
0
)K,W 1
,DW,T W
S T3
230 FORMAT(
I4
,3F1
2
.0
)
S T3
WRITE(6,2
4
0
)K,WI,DW,T W
240 FORMAT(IHO,3HNO.I
4
,5X,3HWW=El5.7
,5X,3HDW=El5.7
,5X, S T3
3HTW=El5.7
)
S T3
W2=Wlー D W
S T3
WS=DV T W
S T3
）
キ1
0
0
.
0
WP=((W1-DW)/(DW-TW）
S T3
WRITE(
6
,2
5
0
)W 2
,WS,WP
250 FORMAT(IH,5X, 3HWW=El5.7
, 5X, 3HWS=El5.7
,5X,3HWP=E S T3
1
3
.5
,8HPERCENT)
S T3
A M WI=AMWI+(
(
(
W1-DW)/(DW-TW))*IOO.0
)
S T3
4
0 CONTINUE
S T3
AMW=AMWI/FLOAT(M)
S T3
WRITE(
6
,2
6
0
)A M W
S T3
2
6
0 FORMAT(IH,5X,IOHHEIKINCHI=E13.5, 8HPERCENT)
S T3
READ(
5
,2
7
0
)AL,DT,DC,DB,W 3
,GS
S T3
270 FORI¥
在AT(
6F1
2
.0
)
S T3
WRITE(
6
,2
8
0
)AL,DT,DC,DB,W 3
,GS
280 FORMAT(!HO,3X,2HL=El3.5,3X,3HDT=E13.5,3X,3HDC=E13.5, S T3
S T3
3X,3HDB=El3.5,3X,2HW=E15.7
,3X,3HGS=El3.5)
S T3
DAV=(DT+2.O*DC+DB)*O・ 2
5
S T3
ホ0
.25**3*(DT+2
.O*DC+DB)**2
A=3.14
S T3
V=3
.1
4
*
0
.25**3*(DT+2
.O*DC+DB)≪*2*AL
S T3
）
料2*AL)
GMT=W3/(3.14*0.25林弘（DT+2.O*DC+DB
林2
*
A
L
)
/
(
l
.0十 AMW/100.0) S T3
GMD=W3
/
(
3
.1
4
*
0
・ 25**3*(DT+2.O*DC+DB）
料 3
*(DT+2.O*DC+DB)**2* S T3
E=(l.O+AMW/
1
0
0
.O)*GS/(W3/(3.1
4
*
0
.2
5
AL）
）
ーI
.0
（
キDT+2.0*DC
SR=AMW*GS/(l.O+AMW/
1
0
0
.O)*(GS/(W3
/
(
3
.1
4
*
0
・2
5林 3
+DB)**2*AL
）
）
ーI
.0
)
WRITE(
6
,2
9
0
)DAV,A
,V,GMT,GMD,E
,SR
290 FORMAT(IH,3X,4HDAV=E13.5,3X,2HA=El3.5,3X,2HV=E14.6,
3X, IOHGAMMERT=E13.5,3X, IOHGAMMERD=El3.5/IH, 5X,
2HE=El3.5,3X,3HSR=E13.5)
3
0 CONTINUE
S T3
S T3
STOP
S T3
END
DATA
$
S T3
S T3
1
0
3
Y -I
4
3
2
1
.4
3
1
1
7
.
6
2
9
1
0
.
3
4
8
Y -I
1
8
8
3
1
.910
2
4
.7
9
1
1
3
.
3
2
9
Y -I
2
5
.
6
0
7
2
0
.
6
8
9
1
8
2
1
2
.7
6
3
Y -I
9
3
94
9
5
96
9
7
I
2
3
4
'
{-I 5
Y -I 6
Y -I 7
Y I 8
Y -I 9
Y -I 1
0
26
佐賀大学農学業報第 2
勾
1
7
.1
0
9
20.970
2
1
.8
9
1
5
.
0
1
1
1
.978
1
3
.
3
2
9
1
2
.7
6
3
5
.
2
2
.
5
5
1
25.645
26.664
3
.70
1
8
.
0
7
4
2
0
.1
7
7
20.620
3
.7
0
9
.918
9
.
9
8
1
1
0
.1
4
4
3
.
6
9
1
5
1
.5
3
2
.
6
8
1
28.273
2
4
.
2
0
7
泡.
5
1
6
2
4
.
6
1
2
1
.248
1
8
.9
7
1
22.494
4.96
1
0
.
3
4
8
I
O
.
l
l
4
1
1
.978
5
.
0
1
3
6
3
.7
1
2
.
5
2
3
。
。
2
.
6
8
1
3
9
9
.
9
7
3
1
5
6
220
l
l
O
8.92
3
43
1
9
237
1
1
.82
y -I
y -I
y-I
y -I
y -I
Y -I
Y -I
Y -I
y -I
y -I
y -I
y I
Y-I
Y-I
y I
の
2
0
.1
8
4
25.030
26.405
4.97
1A’
2i’
3A4
5i’
6i’
7A 8
9A0
12345
qLn4n4nLqL 4
’
TA
YA
’
’
3
237
1
8
8
1
8
2
1
2
.
5
5
表ー 7 含水比，湿潤密度，乾燥密度，間ゲキ比，飽和度の計算結果（成形法）
料
試
平均含水比
AMW
乾燥密度
湿潤密度
間
GMT
GMD
I
.4
220g
/
c
c
0
.8
9
5
5g
/
c
c
除
、
比
E
I
5
8
.7
9
2 %
n
m
57.933
I
.4
655
0
.
9
2
7
9
I
.6
42
I
.6
2
5
6
I
.0
544
I
.5
43
N
5
4
.1
7
2
55.405
I
.5829
6
0
.2
7
7
I
.6
426
I
.0
186
!
.
0
2
4
9
I
.6
32
v
I
.7
3
8
I
.4
62
飽
和度
SR
。
ム
6
5
.756 川
60.545
6
1
.1
6
2
6
6
.3
l
l
50.854
パラフィン塗布法の場合
(a) 記号の説明
N
,M,K,Wl,DW,TW,W2,WS,WP,AMWl,AMW,GS,GMT,GMD,E
,SR
は（2
）の説明と同じ
W4
W5
：パラフィンを塗った供試体重量
：供試体の水中重量
Vl
WPW
VPW
V2
W6
V3
：パラフィンを塗った供試体体積
：パラフィンの重量
：パラフィンの体積
:Vlかまたは（W 4 羽T
5
)
：供試体の重量
：供試体の体積
J
LU
，
‘
、
フローチャート 4 3
2
7
生島・細山・中島：電子計算機の利用による土の物理試験データ
M 回繰り返す
IK,Wl,DW,TWを読む／
I
K Wl,DW TW を印刷／
する
I
I
W2,ws,WPを計算する｜
IW2,ws,WPを印刷する／
IAMWlを計算する｜
IAMW を計算する
IAMWを印刷する
｜
／
l
W4
,W 5
,V1
,WPW,VPW,
GSを読む
W4
,W 5
,V1
1WPW,VPW,
GSを印刷する
↓
I V2を計算する
｜
W6
,V3
,GMT,GMD,E
,
SRを計算する
l
W6
,V3
,GMT,GMD,E
,
SRを印刷する
C三二〉
(c) プログラムおよび結果
乙の場合のプログラムは Program4-3であって，結果は表− 8I
L示す．
事
c
c
HARP
pHISICALTESTFORSOIL4 3
OPERATIONSO N BULKDENSITY BY PARAFFIN WAX PAINTING
METHOD
WRITE(
6
,3
0
0
)
3
0
0 FORMAT(IHI,13HS0ILTESTι3/
IHO,5
7HOPERATIONONBULK
DENSITYBYPARAFFIN l
o
/
A
XPAINTINGMETHOD)
READ(
5
,3
1
0
)N
、
S T3 1
0
1
S T3 1
0
2
S T3
S T3
S T3
ST3
1
0
3
1
0
4
1
0
5
1
0
6
28
4
号
佐賀大学農学業報第 2
310 F0RMAT(I4)
D050I=l,N
READ(
5
,3
2
0
)M
320 FOR!¥
在A
T(I4)
AMWl=O.O
D060J=l,M
READ(
5
,3
3
0
)K,WI,DW,T W
,3F1
2
.0
)
3
3
0 FORMAT(I4
WRITE(
6
,3
4
0
)K,WI,DW,T W
3
4
0 FORMAT(IHO,3HN0.I
4
,5X,3HWW=El5.7
,5X,3
田）W=El5.7
,5X,
3HTW=El5.7
)
W 2＝もVl-DW
WS=DW-TW
WP=((W1-DW)/(DW-TW))*lOO.0
WRITE(
6
,3
5
0
)W 2
,WS,W P
3
5
0 FORMAT(IH,5X, 3HWW=El5.7
,5X,3HWS=El5.7
,5X, 3HWP=E
1
3
.5
,8HPERCENT)
A M W1=AMWl+(((W1 DW)/(DW-TW))*IOO.0
)
60 CONTINUE
A恥1W=A M W1
/FLOAT(M)
WRITE(
6
,3
6
0
)AMW
3
6
0 FORMAT(lH,5X,IOHHEIKINCHl=El3.5,8HPERCENT)
READ(
5
,3
7
0
)W 4
,W 5
,V 1
,WPW,VPW,GS
2
.0
)
3
7
0 FORMAT(6F1
WRITE(
6
,3
8
0
)W 4
,W 5
,VI,WPW,VPW,GS
3
8
0 F0RMAT(lH0,3X,3HWl=El3.5,3X,3HW2=El3.5,3X, 3HVI=E
1
3
.5
,3X,3HWP=E1
3
.5
,3X,3HVP=El3.5
,3X,3HGS=E1
3
.5
)
事
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
ST3
1
0
7
1
0
8
1
0
9
l
l
O
1
1
1
1
1
2
1
1
3
1
1
4
l
1
5
1
1
6
ST3
ST3
ST3
ST3
ST3
l
1
7
1
1
8
1
1
9
1
2
0
1
2
1
S T3
S TS
S T3
S T3
S T3
S T3
S T3
S T3
S T3
ST 3
S T3
V2=Vl
S T3
I
F(V2
.GT.0
.0
)GOTO1
S T3
V2＝
、
へT4 へ
、T5
S T3
W6=W4一
、
へTPW
S T3
V3=V2-VPW
S T3
GMT=(W4 WPW)/(V2-VPM)
S T3
GMD=((W4-WPW)/(V2 VPW))/(l.O+AMW/100.0
)
S T3
）
）
ー
！
.0
E=(l.O+AMW/
1
0
0
.O)*GS/((W4-WPW)/(V2-VPM
S T3
SR=AMW*GS/((l.O+AMW/
1
0
0
.O)*GS/((W4-WPW)/(V2-VPW)))
S
T3
WRITE(
6
,3
9
0
)W 6
,V3
,GMT,GMD,E
,SR
3
9
0 FORMAT(IH,3X,2HW=El3.5,3X,2HV=El3.5,3X,IOHGAMMER S T3
T=El3.5, 3X, lOHGAMMER D=El3.5/1H, 5X, 2HE=El3.5, 3X,3 S T3
HSR=El3.5
)
ST3
50 CONTINUE
ST3
STOP
ST3
END
ST3
DATA
ST3
1
2
2
1
2
3
1
2
4
1
2
5
1
2
6
1
2
7
1
2
8
1
2
9
1
3
0
1
3
1
1
3
2
1
3
3
1
3
4
1
3
5
1
3
6
1
8
7
1
3
8
1
3
9
1
4
0
1
4
1
1
4
2
1
4
3
1
6
.
8
9
0
1
6
.
0
7
0
1
7
.
3
5
0
0.0
1
1
.6
0
1
1
1
.2
3
3
1
1
.9
2
2
83.0
B-2
B -2
B -2
9
5
.
5
2
.
6
4
B-2
J
。
B-2
1
8
.943
1
8
.003
1
9
.
4
8
0
7
2
7
.
0
21qLoad
せに
3
l
l
1
2
1
3
2
0
4
8
.
1
4
4
1
4
5
1
4
6
1
4
7
1
4
8
2
9
生島・細山・中島：電子計算機の利用による土の物理試験データ
表− 8 合水比，湿潤密度，乾燥密度，間ゲキ比，飽和度の計算結果
（パラフィン塗布法）
｜平均含水比｜湿潤密度｜乾燥密度｜間
試
AMW
GMT
I
.6
03
GMD
I
.1
5
0
(5) 粒度試験
2
,0
0
0
μ フルイ残留部分
(a) 記号の説明
N
：試料数
ANZ
：含水比測定の容器番号
WWl ：含水比測定の湿潤重量
DWl
：含水比測定の乾燥重量
T羽T
l
：含水比測定の容器重量
W W 2 ：含水比測定の水分重量＝ WWl-DWl
WSl
：含水比測定の土粒子重量＝ DWl-TWl
Wl
：含水比
W3
:；
｝羽T
l
WPl
：平均含水比
ANl
：空気乾燥試料の容器番号
W W 3 ：（空気乾燥試料＋容器）重量
TW2
：空気乾燥試料の容器重量
BNl
：乾燥試料の容器番号
W W 4 ：（乾燥試料＋容器）重量
T羽T3 ：乾燥試料の容器重量
W2
：空気乾燥試料重量
WS2
：試料の乾燥重量
WO
BN2
BN3
：乾燥試料重量
w羽T5
：（残留土＋容器）重量
TW4
：容器重量
WS3
：残留土重量
羽T
4
：残留率
W5
WP2
：加積残留率
：フノレイ番号
：容器番号
：加積通過率
(
b
) フローチャート 5-1
際
E
比
I
.2
94
飽
和
度
SR
45.270
%
3
0
佐賀大学農学業報第 24
号
C平
コ
N回繰り返す
＋
W3=0.0
3回繰り返す
↓
2
,W W1
,D W1
,T W1を読む
↓
AN2
,W W1
,D W1
,TWl を印刷する
↓
W W2
,ws1
,w 1を計算する
7
7
L
A
N
ど
I
LW W2,ws1,w 1
を印刷する
I
I
｜
7
↓
W 3を計算する
よ
↓
WPl を計算する
AN1 W W3 T W2.BN1,W W4
TW3 を読む
I
I
ANl W W山川 1
,WW4,
TW3 を印刷する
IW2,WS2,W O
を計算する
↓
WP1
,W 2
,WS2
,W Oを印刷l
する
ど
7
W5=0.O
↓
7回繰り返す
LBN2,BN3;!!_W5,T W4を読む 7
↓
jws3,W4,w 刊
P2を計算する｜
[ws3,w リ刊P2を印刷する
7
!
亡
二Eコ
(c) プログラムおよび結果
この場合のプログラムは Program5-1であって，結果は表
sr
e示す．
回
官且
a
m
a
r
r
s
−
o
p
3
1
生島・細山・中島：電子計算機の利用による土の物理試験データ
事
c
c
c
HARP
PHISICALTESTFORSOIL5
1
OPERATIONSONMECHANICALANALYSISBYJISA 1
2
0
4
COARSEGRAINEDOFSIZEDISTRIBUTION
WRITE(
6
,1
1
0
)
1
1
0 FORMAT(lHl,5X,13HS0ILTEST5-1/lHO,5X,46H0PERATIONON
MECHANICALANALYSISBYJISA 1
2
0
4
)
READ(
5
,1
2
0
)N
1
2
0 FORMAT(
I4
)
DO1
0I=l,N
W3=0.0
D020J=l,3
READ(
5
,1
3
0
)AN2,W WI
,D Wl
,T Wl
1
3
0 FORMAT(
4F1
2
.0
)
,D Wl
,T Wl
WRITE(
6
,1
4
0
)AN2,W Wl
1
4
0 FORMAT(lH0,3HN0.EI3.5,3X,3HWW=El3.5,3X,3HDW=El3.5,
3X,3HTW=El3.5)
WW2=WWI-DWl
WSl=DWI-TWl
W !=(WW1-DWl)/(DW1-TW1
)
*
1
0
0
.
WRITE(
6
,1
5
0
)W W2
,WSI
,W l
.5
,3X,3HWS=E1
3
.5
,3X,2HW=El3.5
)
1
5
0 FORMAT(
lHb,3X,3HWW＝引 3
（ W W1-DWl)/(DW1-TW1
)
*
1
0
0
.
W3=W3十
20 CONTINUE
WPl=W3/3.
,T W2
,BNI
,W W4
,T W3
READ(
5
,1
6
0
)ANl
,W W3
1
6
0 FORMAT(6F1
2
.
0
)
WRITE(
6
,1
7
0
)ANI
,W W3
,T W2
,BNI
,W W4
,T W3
1
7
0 FORMAT(lHO,9HKIKANN0.EI3.5,3HWW=El5.7,3HTW=El5.7,
SX,9HKANSONO.El3.5,3HWW=El5.7
,3HTW=El5.7
)
W2=WW3 TW2
WS2=1
0
0
.O*(WW3-TW2
)
/
(
1
0
0
.O+WPI
)
W0=WW4-TW3
WRITE(
6
,1
8
0
)WPI
,W 2
,WS2,W O
1
8
0 FORMAT (!Hb, 5X, 19HHEIKIN GANSUIHI W=El3.5, 14HKIKAN
SIRYOW=El5.7
,14HSIRY0KANSOW=El5.7
,8HKANS0W=El5.7
)
W5=0,0
D030J=l, 7
,T W4
READ(
5
,1
9
0
)BN2
,BN3
,W W5
1
9
0 FORMAT(
4F1
2
.0
)
WRITE(
6
,2
0
0
)BN2
,BN3
,W W5
,T W4
乱1AT(IHO,6HFURUI.El3.5,8HMILLIM, 3X, 3HN0.El3.5,3X,
200 FOR
6HDW+TW=El5.7
,3X,3HTW=E15.7
)
WS3=WW5 TW4
W4=(WW5-TW4）
川V
S
2
*
1
0
0
.
川V
S
2
*
1
0
0
.
W 5=W5+(WW5-TW4）
WP2=1
0
0
.0-W5
WRITE(6,2
1
0
)WS3，もγ4,W5,WP2
210 FORMAT(!Hb,3X,3HWS=E15.7,3X,3HW4=El3.5,3X,3HW5=E
1
3
.5
,3X,3HWP=E1
3
.5
)
3
0 CONTINUE
1
0 CONTINUE
STOP
END
DATA
。
。
。
。。
$
ST4
ST4
ST4
ST4
S T4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
2
3
4
5
6
7
8
9
1
0
1
1
1
2
1
3
1
4
1
5
ST4
ST4
ST4
ST4
S T4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
1
6
1
7
1
8
1
9
2
0
2
1
2
2
2
3
24
2
5
2
6
2
7
2
8
2
9
3
0
3
1
3
2
3
3
3
4
3
5
3
6
3
6
3
8
3
9
40
4
1
4
2
4
3
44
4
5
4
6
4
7
4
8
4
9
5
0
5
1
5
2
5
3
5
4
3
2
侍
佐賀大学農学業報第2
1
5
.
4
8
4
1
5
.7
8
5
1
6
.1
1
6
5
1
1
0
.0
0.0
6.0
5.0
4.0
3.0
2.0
1
.0
1
4
.1
2
4
1
4
.
2
8
6
1
4
.7
4
2
9
5
5
.
0
0.0
4
1
0
.
0
5
7
7
.
0
3
8
4
.
0
4
3
4
.
0
4
7
9
.
0
3
8
2
.
0
8
.
3
6
2
7
.
9
3
1
8
.
9
2
2
4.0
0.0
2
5
3
.
0
2
4
8
.
0
2
7
3
.
0
2
3
2
.
0
2
8
1
.0
2
3
4
.
0
1
7
9
0
.
0
EXl
2
EXl
EXl
3
6
4
5
.0 E X1 4
EXl
5
EXl
6
EXl
7
EXl
8
EXl
9
EXl 1
0
EXl 1
1
表− 8 粒度分析の計算結果（2.0mm以
上）
Rd
ny
一−
2
,0
0
0
μ フルイ通過部分
(a) 記号の説明
N
：試料数
AN4
：含水比測定の容器番号
羽TW6
：合水比測定の湿潤亙−畳
DW2
T羽T5
：含水比測定の乾燥重量
w羽T7
：含水比測定の水分重量＝ WW6-DW2
：含水比測定の容器重量
WS4
W7
W9
WP3
AN3
WW8
TW6
GS
GT
CM
T2
：測定時の平均温度
W8
：：気乾試料重量
：含水比測定の土粒子重量＝ DW2-TW5
：合水比
： ~W7
：平均含水比
：気乾試料重量測定の容器番号
：（気乾試料＋容器）重量
：容器重量
：土粒子の比重
:T2における水の比重
：メニスカス補正
明T
S5
：試料の乾燥重量
vws
:1
0
0
/（
羽T
S
5
/
V
)=100／
（
羽T
S
5
/
l
0
0
0
)
:GS/(GS-GT)
:
v
w
s×GST2
GST2
AM
WP2
A
B
：粒径 2,0mm以上における加積通過率
：）有効深さと比重計の読みの関係式の係数
h
m一
4
3
一
m，
一$
;
J
o
o一
。
、
“
一
竺
9
i−
’・−
川一パ一
1
.0
2
.
0
3.0
3
.
0
5
0
.
8
3
8
.1
25.4
1
9
.1
9
.
5
2
4
.7
6
2
.
0
0
5
0
.
8m m
%
1
0
0
.
0
0
生島・細山・中島：電子計算機の利用による土の物理試験データ
Al
Bl
Cl
Tl
Rl
R2
AL
TL!
}
1
i
l
i
i
E
i
U
H… 式の係数
：測定時間
：比重計のよみ
:Rl+CM
：有効深さ
:AL/Tl
TL2
：仏L/Tl=.JTIτ
GSTl
／
〔9
8
0(GS-GT
〕
）
:,
,
/
3
0甲
D
：粒径
F
FR
：補正係数
p
:R2+F
:2
0
0メッシュ以下の加積通過率
:2
0
0メッシュ以下の補正加積通過率
WP4
BN4
WW9
BN5
TW7
：容器重量
羽T
S6
：残留土重量
WlO
Wll
P2
：加積残留率
WP5
：フルイ番号
：（残留土＋容器）重量
：容器番号
：残留率
:2
0
0メッシュ以上の加積通過率
:2
0
0メッシュ以上の補正加積通過率
(b) フローチャート 5-2
I
C
手
コ
Nを読む
I
N回繰り返す
I
W9=0,0
I
3回繰り返す
AN4,WW6
日DW2,TW5を読む
AN4,WW6,DW2,TW5を印刷
する
↓
IWW7,WS4,W7,W9を計算する｜
↓
IWW7, ws久 W 7 を印刷する l
↓
↓
3
3
佐賀大学農学会報第 24
号
34
l
IWP3を計算する｜
↓
IWP3を印刷する I
↓
AN3,WWB,TW6,G
S
,CM,T2
を読む
↓
AN3,WWB,TW6
日G
S
,CM,T2
を印刷する
↓
WB,WSS,VWS,GST~ ， A Mを
計算する
↓
WB
日WSS,VWS,GST2,A Mを
印刷する
↓
WP2 を読む
↓
A
,B を読む
I
I
／
／
↓
A 1= 0
.537E-05
B 1'=0.2E-O
唱
C1= -0.90096E-03
8回繰り返す
↓
T1
,R 1
，を説む
／
↓
R2
,AL, TL1
,TL2
, GST1
,D
,
F
,FR,P
,WP4 を計算する
↓
R2
,AL;TL1
, TL2
, GST1
,D
,
F
,FR,P
,WP4 を印刷する
I
止
I
↓
w u=0.o
I
↓
5回繰り返す
↓
J
"
T
ζ二
互コ
(c) プログラムおよび結果
この場合のプログラムは Program5-2であって，結果は表− 9I
乙示す．
3
5
生島・細山・中島：電子計算機の利用による土の物理試験データ
Program 5-2
$
c
c
c
HARP
PHISICALTESTFORSOIL5
2
OPERATIONSONMECHANICALANALYSISBYJISA 1
2
0
4
FINEGRAINEDOFSIZEDISTRIBUTION
WRITE(
6
,2
2
0
)
2
2
0 FORMAT(1H1,5X,13HSOIL_TEST52/1HO,5X,46HOPERATI0NON
h在ECHANICALANALYSISBYJ
I
SA 1
2
0
4
)
READ(
5
,2
3
0
)N
230 FORMAT(
I
4
)
004
0I=l,N
W9=0.
DO50J=l,3
READ(
5
,2
4
0
)AN4
,W W6
,D W2
,T W5
2
4
0 FORMAT(
4F1
2
.0
)
WRITE(
6
,2
5
0
)AN4
,W W6
,D W2
,T W5
250 FORMAT(lHO,3HNO.E1
3
.
5
,3X,3HWW=El5.7
,3X,3HDW=El5.7
,
3X,3HTW=El5.7
)
WW7=WW6-DW2
WS4=DW2-TW5
W 7=(WW6D W2)/(DW2-TW5
)
*
1
0
0
.0
W9=W9+(WW6-DW2)/(DW2-TW5
)
*
1
0
0
.
WRITE(
6
,2
6
0
)W W7
,WS4,W 7
2
6
0 FORMAT(
IHb,3X,3HWW=El5.7
,3X,3HWS=El5.7
,3X,2HW=El3.5
)
5
0 CONTINUE
WP3=W9/3.
WRITE(6,2
6
1
)WP3
2
6
1 FORMAT(IHb,IOX,17HHEIKINGANSUIHI=El5.7
)
READ(
5
,2
7
0
)AN3,WW8,TW6,GS,CM,T2
2
7
0 FORMAT(6F1
2
.0
)
WRITE(
6
,2
8
0
)AN3
,W W8
,T W6
,GS,CM,T2
280 FORMAT(!HO,3HN0.E
l
3
.
5
,3HWW=El5.7
,3HTW=El5.7
,3HGS=
El5.7
,3HCM=El5.7
,4HTC=El3.5)
W8=WW8-TW6
WS5=100.0*(WW8TW6)/(100.0+WP3)
VWS=1
0
0
.0
/
(
1
0
0
.O*(WW8-TW6)/(100.0+WP3)/IOOO.0
)
GST2=GS/(GS（
（
ー0
.5011E-05)*T2**2＋
（
ー0
.6E-08)*T2+0.1
0
0
0
1
7
9
E+
0
1
)
)
AM=1
0
0
.0
/
(
1
0
0
.O*(WW8-TW6)/(100.0+WP3）
パ0
0
0
.O)*GS/(GS−
（
（
ー0
.5
0
1
1
キ
）T2**2＋
（
ー0
.6E08)*T2+0.1000179E+Ol))
E-05
WRITE(
6
,2
9
0
)W8,WS5,VWS,GST2,A M
290 FORMAT(
IHb,3X,14HKIKANSIRYOW =E1
5
.7
,3HWS=E1
5
.7
,1
1H
1
0
0
.O/WS/V=E1
5
.7
,11HGS/(GS-GT)=E15.7
,2HM=E1
5
.7
)
READ(
5
,2
9
1
)WP2
2
9
1 FORMAT(
F1
2
.0
)
READ(
5
,3
0
0
)A
,B
3
0
0 FORMAT(2F1
2
.0
)
A 1=0.437E-05
Bl=O.2E-08
C1= 0
.9
0
0
9
6E-03
D060J=l,8
READ(
5
,3
1
0
)T1
,R 1
3
1
0 FORMAT(2F1
2
.0
)
R2=Rl+CM
AL=A
（
キRl十 CM)+B
。
。
。
ST4
ST4
ST4
ST4
S T4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
S T4
ST4
ST4
ST4
ST4
ST4
S T4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
S T4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST4
ST 4
S T4
S T4
S T4
6
1
6
2
6
3
6
4
6
5
6
6
6
7
6
8
6
9
7
0
7
1
7
2
7
3
7
4
7
5
7
6
7
7
7
8
7
9
8
0
8
1
8
2
8
3
8
3
1
832
8
4
8
5
86
8
7
8
8
8
9
9
0
9
1
9
2
9
3
94
9
5
9
6
9
7
9
7
1
972
9
8
9
9
9
9
1
992
9
9
3
1
0
0
1
0
1
1
0
2
1
0
3
1
0
4
3
6
事
佐賀大学農学業報第2
4
号
TLI=(A*(RI十 CM)+B)/TI
TL2=SQRT((A*(RI+CM)+B)/TI
)
(
I
.o+o.0337*T2+0.0
0
0
2
2
1＊~T2**2) /
(
9
8
0
.o
*
GST1=SQRT(
3
0
.o*o.0
1
7
8I
ー0
.5
1
1E 05)*T2**2+(0
.6E-08)*T2+0.1
0
0
0
1
7
9
E
+
O
l
)
)
)
)
(GS（
本o
.0
1
7
8
/
(
I
.o+o.0337*T2+
D=SQRT(
(
A
*
(
R1+CM)+B)/Tl)*SQRT(
3
0
.0
（
ー0
.6E-08)*T2+0.1
0.000221*T2**2)/(980.0*(GS(
(0.5011E05)*T2**2＋
000179E+ol))))
F=A1*T2**2+Bl*T2+C1
FR=R1+CM+A1*T2**2+B1*T2十C1
P=Al¥
在
＊ FR
WP4=P*WP2/100.0
WRITE(
6
,3
2
0
)R2,AL,TLl,TL2,GSTl,D
,F
,FR,P，
、
へfP4
320 FORMAT (lHb, 3X, 3HR2=El5.7
, 2HL=El5.7
, 4HL/T=El5.7
,1
0
HSQRT(L/T)=El5.7
,5HGST1=El5.7/lHb,3X,2HD=El5.7
,2HF=
El5.7
,5HR2+F=El5.7
,2HP=El5.7
,3HWP=El5.7
)
6
0 CONTINUE
Wll=O.O
D O70J=l,5
READ(
5
,3
3
0
)BN4
,W W9
,BN5
,T W7
330 FORMAT(4F1
2
.0
)
,BN5
,T W7
WRITE(
6
,3
4
0
)BN4
,W W9
IHO,6HFURUI.E1
3
.5
,3X,3HWW=El5.7
,3X,3HN0.E
l
3
.5
,
3
4
0 FORMAT(
3X,3HTW=El5.7
)
WS6=W W9-TW7
W 10=(WW9-TW7)/WS5*100.0
W 11=WI1+(WW9-TW7)/WS5*100.0
P2=1
0
0
.0−
、
へf1
1
WP5=WP2*P2/100.0
WRITE(
6
,3
5
0
)WS6,W 1
0
,WII
,P2
,WP5
IHb,3X,3HWS=El5.7
,3X,4HWIO=E1
5
.7
,3X,4HW1
1=E
3
5
0 FORMAT(
1
5
.7
,3X,2HP=E15.7
,3X,3HWP=El5.7
)
7
0 CONTINUE
40 CONTINUE
STOP
END
DATA
．
．
．
nUAUnUP3P3nUP3nuoa J
句 dTAUFD
戸コヮ“マ’つ＆ nUQUFhJ
QU 勺J F 3 q d Q u
qLqL’
l’
I’
inununun
日 nUAUAU ’
i
AUAUnunり nununU
“
ヮAり・
司
1
7
.3
5
0
1
6
.
8
9
0
1
6
.
0
7
0
3
2
.
3
5
1
1
.9
2
2
II
.6
1
0
II
.2
3
3
2
.8
5
0
4
0
1
.0
4
0
2
.
0
4
0
3
.
0
4
0
4
.
0
4
0
5
.
II
.8
6
0
1
1
.7
4
0
1
1
.010
1
1
.0
8
0
II
.6
1
0
l
必
仙
Aη4qL の
3η4
．0
．．
4 0．
0．
0．
0
0．
o o，
l
−−−−
，
勺
qL
CO
の
令d
AT
nu
5
1
0 0 0 0 91
0 0 0 0 0 0 0 0 0 0 0 5 4J
nLRJvnU
1 2 3 2 5 1 2 5 5 0 0 0 0 8 A4
210
6
1
3
6
4
4
・
4 d守口U n u n u n u n U
31
・
の
1
9
.
4
8
0
1
8
.
9
4
3
1
8
.
0
0
3
1
0
2
.
3
5
。
0
.0
0
0
5
2
1
.0
S T4
ST4
S T4
ST4
ST4
S T4
S T4
S T4
S T4
S T4
ST4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
S T4
EX
EX
EX
EX
EX
EX
EX
EX
EX
EX
EX
EX
EX
EX
EX
EX
EX
EX
EX
1
0
5
1
0
6
1
0
7
1
0
8
1
0
9
1
1
0
1
1
1
1
1
2
1
1
3
1
1
4
1
1
7
1
2
0
1
2
1
1
2
2
1
2
3
1
2
4
1
2
5
1
2
6
1
2
7
1
2
8
1
2
9
1
3
0
1
3
1
1
3
2
1
3
3
1
3
4
1
3
5
1
3
6
1
3
7
1
3
8
1
3
9
1
4
0
1
4
1
1
4
2
1
4
3
1
4
4
1
4
5
1
2
1
3
1
4
1
5
1
5
1
1
6
1
7
1
8
1
9
2
0
2
1
2
2
2
3
2
4
2
5
2
6
2
7
2
8
2
9
生島・細山・中島：電子計算機の利用による土の物理試験データ
3
7
表− 9 粒度分析の計算結果（2.0mm以下）
フィル分け試験
｜
比重計試験
I
μ420I
μぉoI
μlosI
.
μ 74I
μ44.t
i32.:
t
i2i.t
i 12.t
ia.9I
μ6.4I
μ 3.2i
μ i.3μ
a40
5.74165.30163.93161.37160.61153.56147.50137.40128.32124.28119.23I 1
4
.
1
9
1 7
.
1
2
(6) 液性限界，塑性限界試験
(a) 記号の説明
N
：試料数
M
：測定数（液性限界試験）
AN(J) ：落下回数（液性限界試験）
K
：容器番号（液性限界試験）
Wl
：湿潤重量（液性限界試験）
DWl
T羽T
l
：容器重量（液性限界試験）
：乾燥重量（液性限界試験）
W2
：水分重量（液性限界試験）
WSl ：土粒子重量（液性限界試験）
羽r
p(
J
) ：含水比（液性限界試験）
X2
Xl
XY
Yl
A
B
:~ (
l
o
g
1
0AN(
]
)
)
2 （液性限界試験）
:~ l
o
g
1
0AN(
J
) （液性限界試験）
:2 (
l
o
g
1
0AN(
J
)
)(WP(
J
)
) （液性限界試験）
:2WP(J) （液性限界試験）
：流動曲線の係数（液性限界試験）
：流動曲線の切片（液性限界試験）
vv
:2 {WP(J)一
（A l
o
g
1
0AN(])+B
)
2
} （液性限界試験）
ESM
EM
ALL
AIF
L
：中等誤差（液性限界試験）
：液性限界（液性限界試験）
：確率誤差（液性限界試験）
：流動指数（液性限界試験）
：測定数（塑性限界試験）
K
：容器番号（塑性限界試験）
W3
D羽T2
：湿潤重量（塑性限界試験）
TW2
W4
：容器重量（型性限界試験）
WS2
WPP
AMWl
PL
AIP
AIT
：土粒子重量（塑性限界試験）
：乾燥重量（塑性限界試験）
：水分重量（塑性限界試験）
：含水比（塑性限界試験）
:XWPP （塑性限界試験）
：塑性限界（塑性限界試験）
：塑性指数
：タフネス指数
，
、b
3
8
4
号
佐賀大学農学業報第 2
J
フローチャート
6-1
↓
AN(j),K,W 1
,D W1
,T W1を
印刷する
W2,WSl,WP（］）を計算する
W2
,WS1
,WP(
j〕を印刷する
X2
,X 1
,XY,Y 1を計算する
↓
A, B を
.
1
1
n
:する
A, B を印刷する
↓
VV=0.0
M 回繰り返す
↓
vvを計算する
↓
＋
ESM,EM を計算する
ESM,EM を印刷する
ALL,A
I
F
.を計算する
ALL,AIFを印刷する
Lを読む
↓
AMWl=O.0
↓
L回繰り返す
r
K
，
而3
,DW
瓦王山云品五7
と
.
.
.
！
！
＿3
,D W2
,TW2 を印刷する／
3
9
生島・細山・中島：電子計算機の利用による土の物理試験データ
W4
,W
S2,WPP,AMWlを計算
する
Iw4, ws2,WPP を印刷する／
IPL,AIP,AITを計算する
I PL,AIP,A汀を印刷する
｜
／
1
C二Eコ
(c) プログラムおよび結果
この場合のプログラムは P
rogram6 1であって，結果は表− 1
0K示す．
Program 6-1
$
c
c
c
1
0
0
1
1
0
1
2
0
1
3
0
1
4
0
1
5
0
HARP
PHISICALTESTFORSOIL6
1
OPERATIONSONLIQUIDLIMITBYJISA 1
2
0
5
OPERATIONSONPLASTICLIMITBYJISA 1
2
0
6
DIMENSIONAN(
2
0
)
,WP(
2
0
)
WRITE(
6
,1
0
0
)
FORMAT(IHI,13HSOILTEST6ー
！
／!
HO, 7
1HOPERATIONONLIQUID
LIMITBYJ
I
SA 1
2
0
5ANDPLASTICLIMITBYJ
I
SA 1
2
0
6
)
READ(
5
,l
l
O
)N
FORMAT(
I
4
)
DO1
0I
,N
READ(
5
,1
2
0
)M
FORMAT(
1
4
)
X2=0.0
Xl=O.O
XY=O.O
Yl=O.O
0020]=1,M
READ(
5
,1
3
0
)AN(
J
)
,K
,W 1
,O W1
,T W1
F0RMAT(Fl2.0,I4,3Fl2.0)
WRITE(
6
,1
4
0
)AN(
J
)
,K
,W 1
,O W1
,T W1
FORMAT(lHO,7HKAISUU=Ell.3,3X, 3HN0.I4,3X,3HWl=El5.7,
3X,3HDW=El5.7
,3X,3HTW=El5.7
)
W2=Wト DWl
WSl=DWト TWl
)
)
*
1
0
0
.0
WP(J)=((W1O Wl)/(DWト T W1
WRITE(
6
,1
5
0
)W2,WS1
,WP(J)
FORMAT(
IH ,3X,3HWW=El5.7
,3X,3HWS=El5.7
,3X,3HWP=El3.5
,
8HPERCENT)
］
）
）
）
料2
X2=X2+(ALOG1
0(AN（
X 1=X1+(ALOG10(AN(])))
XY=XY+(ALOG1
0(AN(])))*WP(])
J
)
Y1=Y1+WP(
ST5
ST5
ST5
ST5
ST5
S T5
ST5
ST5
ST5
ST5
S T5
S T5
S T5
S T5
S T5
S T5
S T5
S T5
S T5
S T5
ST5
S T5
S T5
S T5
S T5
S T5
S T5
2
3
3
1
4
5
6
7
8
9
1
0
l
l
1
2
1
3
1
4
1
5
1
6
1
7
1
8
1
9
2
0
2
1
2
2
2
3
2
4
2
5
2
6
S T5
S T5
S T5
S T5
2
7
2
8
2
9
3
0
40
S T5
S T5
S T5
S T5
S T5
3
1
3
2
33
34
3
5
57
58
59
60
6
1
62
S T5
S T5
S T5
S T5
S T5
S T5
S T5
S T5
S T5
S T5
S T5
K -I
K-I
K I
K-I
K -I
K I
K-I
K-I
K I
6
2
1
63
64
65
66
67
68
69
7
0
7
1
72
I
2
3
4
5
6
7
8
9
。。
崎
’
n
ヨ
GU 匂J n y n訪日UTAq403dTRJEUmJoon3nu lq4038TFhJcu
T d A a T d守 aTdyaA
3 3034TAT 凋
仙 ATP3FhJFhJFhJRJFhJF3
勺，
RJFaFakJFhJpapaF3F3F3FaFhJ5kJF3FaP3RJpapaFa
叫
S T5
S T5
S T5
S T5
S T5
S T5
cdpapδοδοδPDoδP303papacdoboδoacdoapDODO303
﹃
OOH30303
P昌広 U 氏U
nu
i n dnUFOτi
J’
。J
．．
．．
．
．
h u n，
inwdnu
’
，
同
1i1
1u’
1
1
20 CONTINUE
A=(FLOAT(M)*XY-XI
*Y !)/(FLOAT(M
）
材2-X1*XI
)
）
料2-Xl*XI
)
B=(X2*Yl X1*XY)/(FL0AT(M
WRITE(
6
,1
6
0
)A
:
,B
1
6
0 FORMAT(!HO,5X,IOHFLOWCURVE,5X,2HW=El5.7
,9HL0GIO(N)
+E1
5
.7
)
VV=O.O
D030J=l,M
（
ーA*ALOG1
0(AN(J))+B）
）
紳2
VV=VV+(WP(J)
30 CONTINUE
回 M=SQRT(VV/FLOAT(
M
3
)
)
E
l
¥
在＝ 0
.4
7
6
9
3
6
3
*
1
.4
1
4
2
1
4
本ESM
WRITE(
6
,1
7
0
)ESM,E M
1
7
0 FORMAT (IH, 7X, 20HMEAN SQUARE ERR0R=E15.7, 5X, 1
6
E
王PROBABLEERR0R=E1
5
.7
)
ALL=A*ALOG1
0(
2
5
.O)+B
AIF= A
WRITE(
6
,1
8
0
)ALL,AIF
1
8
0 FORMAT(
IHO,5X,3HLL=E1
3
.5
,5X,3HIF=E13.5
)
READ(
5
,1
9
0
)L
1
9
0 FORMAT(
!
4
)
AMWl=O.O
DO40J=l,L
,T W2
READ(
5
,2
0
0
)K,W 3
,D W2
200 F0RMAT(I4,3F12.0)
WRITE(
6
,2
1
0
)K
,W 3
,D W2
,T W2
210 F0RMAT(IH0,3HN0.I4,3X,3HWW=E15.7,3X,3HDW=E15.7,3X,
3HTW=E15.7
)
W4=W3DW2
WS2=DW2-T
へ
もr
2
WPP=((W3-DW2)/(DW2-TW2
)
)
*
1
0
0
.0
）
キ1
0
0
.
0
A M W1=AMW1+(((W3-DW2)/(DW2-TW2）
WRITE(
6
,2
2
0
)W 4
,WS2,WPP
220 FORMAT(lH,5X,3HWW=E15.7
,3X,3HWS=E15.7
,3X, 3HWP=E
1
3
.5
,8HPERCENT)
40 CONTINUE
PL=AMW1/FLOAT(L)
AIP=ALL-PL
AIT=(ALL-PL)/AIF
WRITE(
6
,2
3
0
)PL,AID,AIT
230 FORMAT(lH,5X,3HPL=E1
3
.5
,5X,3HIP=El3.5
,5X,3HIT=El3.5
)
1
0 CONTINUE
STOP
END
HARP
3
6
1
4
.
6
4
0
1
6
.1
4
4
39.0
4
7
36.0
1
3
.
4
3
4
1
4
.
9
6
4
34
29.0
1
4
.
9
0
2
1
6
.772
87
1
6
.
2
3
1
2
1
.0
1
6
3
1
8
.
8
0
1
1
6
.
0
1
4
.
0
8
1
1
2
.
5
8
9
6
5
1
5
.
5
2
7
1
3
.
0
1
3
.
6
5
1
1
0
3
1
3
.
3
7
1
1
2
.
5
2
8
1
0
.1
4
6
1
1
1
6
.
5
9
5
1
6
.
3
5
9
1
5
.6
8
1
5
TTTTTTTTTTTTTTTTTTTTT
$
佐賀大学農学業報第 24
号
K -I 1
0
4
1
生島・細山・中島：電子計算機の利用による土の物理試験データ
1
4
6
42.0
40.0
3
3
.
0
23.0
1
9
.
0
1
0
.
0
3
3
3
5
3
9
8
6
4
3
.
0
3
8
.
0
3
3
.
0
1
7
.
0
1
5
.
0
1
1
.0
3
7
2
1
0
6
4
4
表
1
6
.
3
8
4
1
6
.1
2
2
1
5
.
3
7
3
4
2
1
2
5
1
2
8
7
1
5
7
5
4
2
2
.
2
0
0
2
1
.294
2
0
.1
1
0
2
4
.
2
7
0
2
1
.7
2
5
2
4
.1
6
1
1
8
.
0
6
5
1
7
.
4
8
7
1
6
.
3
6
0
1
9
.
8
2
1
1
7
.
8
2
1
1
8
.
9
8
3
1
2
.
3
3
2
1
1
.9
7
1
1
0
.
9
1
9
1
2
.2
1
1
1
1
.6
5
2
1
0
.
6
3
9
1
1
.825
1
0
.6
1
5
9
.
6
6
3
9
8
5
9
2
1
3
1
2
3
9
1
2
5
1
6
.
6
0
9
1
9
.
3
7
4
1
6
.7
9
6
1
8
.
2
5
5
1
7
.
6
9
6
1
8
.9
1
1
1
4
.
3
8
7
1
6
.
6
1
5
1
4
.
5
5
8
1
5
.
3
0
2
1
4
.
9
4
4
1
5
.1
8
1
1
0
.
6
8
6
1
1
.7
7
5
1
1
.3
9
8
1
0
.
4
8
2
1
1
.5
3
3
1
1
.0
8
9
9
.
8
0
4
1
0
.7
1
9
1
0
.
0
4
8
1
0
.
0
9
6
1
0
.1
9
1
9
.
3
4
2
1
1
.7
2
0
1
0
.856
1
0
.
2
2
5
9
.7
4
6
1
0
.
9
3
8
1
0
.
0
2
7
9
.
8
5
8
9
.
9
7
1
8.654
K-1 1
1
K -1 1
2
K -1 1
3
K -1 1
4
K -1 1
5
K-1 1
6
K-1 1
7
K- 1 1
8
K -1 1
9
K - 1 20
K-1 2
1
K - 1 22
K -1 23
K - 1 24
K -1 25
K - 1 26
K-1 2
7
K - 1 28
K - 1 29
K -1 30
K-1 3
1
K-1 3
2
K-1 3
3
1
0 液性限界，流動指数，塑性限界，塑性指数，タフネス指数の計算結果
｜液性限界｜流動指数｜塑性限界 l塑性指数｜タフネス指数
試料｜
I
I
I
I
I
ALL
I
AIF
I
PL
I
AIP
I
AIT
% I
% I
% I
%
4
8
.7
1
6
1
4
.
0
2
9
3
5
.
0
5
9
1
3
.
6
5
7
0
.
9
7
3
5
I
I
5
4
.
5
9
3
1
1
.5
0
2
3
0
.
2
6
6
2
4
.
3
2
7
2
.1
1
5
1
i
l
l
5
1
.5
4
2
1
6
.1
0
1
2
9
.
8
3
4
2
1
.7
0
8
1
.3483
(7) （遠心含水当量試験
し現場含水当量試験
(a) 記号の説明
N
AN
WC
WE
WD
WA
WB
WAD
WBE
：試料数
：ノレツボ番号（遠心含水当量試験）
：ノレツボ重量（遠心含水当量試験）
：かわいたコシ紙の重量（遠心含水当量試験）
：湿ったコシ紙の重量（遠心含水当量試験）
：遠心分離後の（ノレツボ＋湿紙＋土）の重量（遠心含水当量試験）
：炉乾燥後の（ノレツボ十湿紙十土）の重量（遠心含水当量試験）
：
羽TA-WD （遠心含水当量試験）
:WB-WE （遠心含水当量試験）
羽TAB
；
（
羽TA一羽TD
）
一
（
羽TB-WE) （遠心含水当量試験）
WBC
:WB-(WC+WE) （遠心含水当量試験）
WC2(I
）：遠心含水当量（遠心含水当量試験）
WC4 ：平均遠心含水当量（遠心含水当量試験）
佐賀大学農学会報第 24
号
42
WC3
:IWC2(1)-WC2(
2
)
1（遠心含水当量試験）
AN
：容器番号（現場含水当量試験）
Wl
：湿潤重量（現場合水当量試験）
DW
：乾燥重量（現場合水当量試験）
TW
：容器重量（現場含水当量試験）
W2
：水分重量（現場合水当量試験）
ws
：土粒子重量（現場合水当量試験）
WF(I) ：現場含水当量（現場合水当量試験）
WF2
WF3
：平均現場合水当量（現場合水当量試験）
：川T
F(l)-WF(2)I
（現場合水当量試験）
(
b
) フローチャート 7-1
C二互コ
／
Nを読む
／
N回繰り返す
↓
2回繰り返す
j~N
we,w山
W A WB
読む
を
／
AN,WC,WE,WO,WA,WBを
印刷する
,WAD,WBE,WAB,WBC,WC2(1)
を計算する
一
一
一
一
ι
WAD,WBE,WAB,WBC,WC2(1)
を印刷する
f
IWC4, WC3 を計算する
｜
生島・細山・中島：電子計算機の利用による土の物理試験データ
43
,NO
IAN,w 1,DW,TW を印刷する／
I
W2, ws,WF(l）を計算する｜
/w2,w丸 WF(l）を印刷する／
C玉二〉
(c) プログラムおよび結果
との場合のプログラムは Program7-1であって，結果は表− 1
1
1
ζ 示す．
Ip呼
$
C
C
C
200
210
220
230
am 7-1
I
HARP
PHISICALTESTFORSOIL7
1
S T6
OPERATIONSO NCENTRIFUGEMOISTURECONTENTT回 .TBYJISAI207 S T6
ANDFIELDMOISTUREEQUIVALENTTESTBYJ
I
SA 1
2
0
8
S T6
2
)
,WF(
2
)
ST 6
DIMENSIONW C2(
WRITE(
6
,2
0
0
)
S T6
即・ ST 6
FORMAT(IHI, 13HSOILTEST7-1/IHO,82HOPERATIONONGET
FUGE MOISTURECONTENTτ'EST AND FIELD MOISTURE EQUIV- S T6
ALENTTFST)
READ(
5
,2
1
0
)N
ST 6
FORMAT(
I
4
)
S T6
DOIOJ=l,N
ST6
00201=1,2
ST6
READ(
5
,2
2
0
)AN,w
e,WE,WD,wA,W B
sT 6
FORMAT(6F1
2
.0
)
S T6
WRITE(
6
,2
3
0
)AN,WC,WE,WD,WA,W B
ST 6
F RMAT(lHO,3HN0. El3.5,2X,3HWC=El5.7
,2X,3HWE=El5.7
, S T6
2X,3HWD=El5.7
,2X,3HWA=El5.7
,2X,3HWB=El5.7
)
S T6
WAD=W AW D
S T6
WBE=WB-WE
S T6
。
I
2
3
4
4
1
5
6
7
8
9
1
0
1
1
1
2
1
3
1
4
1
5
1
6
1
7
44
佐賀大学農学柔報第2
4
号
WAB=(WA-WD
）
ー
（WB-WE)
S T6 1
8
WBC=WB-(WC+WE)
S T6 1
9
）
ー
（W BWE))/(WB-(WC+W
E))*IOO.0
S T6 2
0
WC2(I)=((WA-WD
WRITE(
6
,2
4
0
)WAD,WEE,WAB,WBC,W C2(
I
)
S T6 2
1
2
4
0 FORMAT(IH,6HWA-WD=El5.7
,2X,6HWB-WE=El5.7
,2X,16H(WA- S T6 2
2
WD
）
ー
（WB-WE)=E15.7
,2X,llHWB(WC+WE)=El5.7/IH,!OX,3HWC= S T6 2
3
E1
3
.5
,8HPERCENT)
2
0 CONTINUE
S T6 2
4
WC4=(WC2(!)+WC2(
2
)
)
/
2
.
S T6 2
5
WC3=ABS(WC2(
!
)-WC2(2))
S T6 2
6
I
F(WC4.LE. 1
5
.0
)GOTOI
S T6 2
7
I
F(WC3
.GT.2
.0
)GOTO2
S T6 2
8
WRITE(
6
,2
5
0
)WC4,WC3
S T6 2
9
2
5
0 FORMAT (IH, 5X, 14HHEIKINCHIWC=El3.5, 8HPER CENT, 3X, ST6 3
0
6HG0SA=E1
3
.5
,8HPERCENT)
S T6 3
1
GOT040
S T6 32
I
F(WC3.GT.I
.0
)GOTO2
S T6 3
3
WRITE(
6
,2
6
0
)W C4
,W C3
S T6 34
2
6
0 FORMAT (!H, 5X, 14HHEIKINCHIWC2=El3.5,8HPER CENT, 3X, S T6 3
5
6HG0SA=E1
3
.5
,8HPERCENT)
S T6 3
6
GOT040
S T6 3
7
2 WRITE(
6
,2
7
0
)WC4WC3
S T6 3
8
IH ,!OX, 14HHEIKINCHI WC2=El3.5,8HPERCENT,3X, S T6 39
2
7
0 FORMAT(
6HGOSA=E1
3
.5
,8HPERCENT,3X,18HSAISIKEN0 YOOSURU)
S T6 40
S T6 4
1
4
0 DO301=1,2
S T6 42
、1,TW
READ(
5
,2
8
0
)AN,WI,Dヘ
2
8
0 FORMAT(4F1
2
.0
)
S T6 4
3
S T6 44
I
F(AN.EQ.0
.0
)GOTO1
0
WRITE(
6
,2
9
0
)AN,WI,DW,TW,
S T6 4
5
2
9
0 FORMAT(!HO,3HN0.El3.5, 2X,3HWW=El5.7
,2X,3HDW=El5.7
, S T6 4
6
2X,3HTW=El5.7
)
S T6 4
7
W2=WI D W
S T6 4
8
WS=DW-TW
S T6 4
9
WF(I)=(WIDW)/(DW-TW)*IOO.0
S T6 5
0
WRITE(
6
,3
0
0
)W 2
,WS,W F(
I
)
3
0
0 FORMAT (!H,3HWW=El5.7
,2X,3HWS=El5.7
,2X,3HWF=El3.5, S T6 5
1
8HPERCENT)
3
0 CONTINUE
S T6 52
S T6 5
3
WF2=(WF(
I
)
+WF(
2
)
)
/
2
.
WF3=ABS(WF(l)-WF(
2
)
)
S T6 54
I
F(WF2
.LE. 1
5
.0
)GOTO3
S T6 5
5
I
F(WF3
.GT.2
.0
)GOTO4
S T6 56
WRITE(
6
,3
1
0
¥WF2,WF3
S T6 5
7
3
1
0 FORMAT (IH,5X, 14HHEIKINCHI WF2=El3.5
,8HPERCENT, 3X, S T6 5
8
6HG0SA=E1
3
.5
,8HPERCENT)
S T6 5
9
GOTO1
0
S T6 60
.0
)GOTO4
S T6 6
1
3 I
F(WF3.GT. I
WRITE(
6
,3
2
0
)WF2WF3
S T6 62
3
2
0 FORMAT (lH, 5X, 14HHEIKINCHI WF2=El3.5,8HPERCENT, 3X, ST6 6
3
6HG0SA=E1
3
.5
,8HPERCENT)
S T6 64
GOTO1
0
S T6 6
5
4 WRITE(
6
,3
3
0
)W F2
,WF3
S T6 66
3
3
0 FORMAT(lH,!OX,14HHEIKINCHIWF2=El3.5,5X,6HG0SA=El3.5, S T6 6
7
5X,8HPERCENT,3X,18HSAISIKEN0 YOOSURU)
S T6 6
8
1
0 CONTINUE
S T6 6
9
STOP
S T6 7
0
。
。
4
5
生島・細山田・中島：電子計算機の利用による土の物理試験データ
END
DATA
事
1
5
.
8
6
1
4
.
3
2
1
3
.
8
1
2
0
.1
5
4.0
5
.
0
1
.0
2
.
0
表
0
.
0
5
0
.
0
7
1
2
.7
3
1
7
.
9
1
0
.1
2
o
.11
定
調担
料
試
5
6
.735%
6
0
.335%
5
8
.985%
6
0
.054%
57.860%
6
0
.194%
2
.250%
0
.279%
均差
平誤
：容器番号
Wl
DW
：湿潤重量
TW
W2
：容器重量
：乾燥重量
：水分重量
ws
：土粒子重量
W3
W4
：含水比
WM
Vl
VO
DWI
TWl
WSl
：平均含水比
、
現場合水当量
NF
(a) 記号の説明
AN
1
0
.
9
4
1
4
.1
8
遠心含水当量
WC2
(8) 収縮常数試験
：試料数
2
3
.
2
6
2
1
.6
8
1
1 遠心含水当量，現場合水当量の計算結果
2
N
2
7
.
5
0
2
6
.
0
2
ST6
ST6
ST6
ST6
ST6
ST6
ST6
:2W3
：湿潤土体積
：乾燥土体積
：乾燥土重量
：容器重量
：乾燥土重量
R
：収縮比
V2
WS2
：収縮体積
:(Vl-VO)/WSl×1
0
0
WS3
：収縮限界
c
：体積変化
ALS
：線収縮
GSl
:l/R一羽TS3/100
GS2
WS4
：土粒子の比重
:1/R-l/GS2
WS5
:Rおよび GSより求めた収縮限界
7
1
7
2
7
3
7
4
7
5
7
6
7
7
46
佐賀大学農学業報第24
号
(
h
) フローチャート 8-1
~コ
I
N を読む
I
N回繰り返す
I
↓
w4
=
o
.o
I
3回繰り返す
IAN,w 1
,DW,TW を読む I
I
AN,WI,DW,TWを印刷する l
I
w2
,w
s
,w3
,w4 を計算する｜
IW2, ws,w 3を印刷する I
t
IW Mを計算する｜
IVI,VD,DWI,TWI を読む l
V
I
,VO,DWI,TWIを印刷する l
I
WSI,R
,V
2
,WS2,WS3,C
,A
L
S
,
G
S
I
,G
S
2
,WS4,WSSを2
十算する
l
VI
,VO,V2
,WS2
,WS3,WS4,
wssを印刷する
VO, WSI
,R
, WM, WS3
,C
,
A
L
S
,GSI
;GS2を印刷する
l
C二Eコ
(c) プログラムおよび結果
rogram8-1であって，結果は表− 1
2I
ζ 示す．
この場合のプログラムは P
Program 8 I
8CC
HARP
PHISICALTESTFORSOIL8
1
OPERATIONO NSHRINKAGECOEFFICIENTTESTBYJISA 1
2
0
9
WRITE(
6
,1
0
0
)
1
0
0 FORMAT(
IHI,13HSOILTEST8-1/IHO,39H0PERATIONO NSHRINKAGECOEFFICIENTTEST)
READ(
5
,1
1
0
)N
1
1
0 FORMAT(
1
4
)
D O1
0 I=l,N
,
y4=0.
、。
S T7
ST7
ST7
ST7
ST7
ST7
ST7
S T7
S T7
I
2
3
4
5
6
7
8
8
1
4
7
生島・細山田・中島：電子計算機の利用による土の物理試験データ
3nu
日
η，
a の
1’
1i1
13
ST7
DO20J=l,3
READ(
5
,1
2
0
)AN,W 1
,DW,T W
ST7
2
.0
)
1
2
0 FORMAT(4F1
ST7
WRITE(
6
,1
3
0
)AN,W 1
,DW,T W
ST7
1
3
0 FORMAT(lHO,3BN0.E
l
3
.
5
,2X,3HWW=El5.
.
;
'
,2X,3HDW=El5.7
, ST7
2X,3HTW=E1
5
.
7
)
W2=Wト D W
ST7
WS=DW-TW
ST7
W3=(Wl DW)/(DW-TW)*IOO.O
ST7
W 4=W4+((W1-DW)/(DW-TW)*lOO.0
)
ST7
WRITE(
6
,1
4
0
)W2,WS,W 3
ST7
1
4
0 FORMAT(lH,5X,3HWW=El5.7
,5X,3HWS=El5.7
,5X,2HW=El3.5, ST7
8HPERCENT)
2
0 CONTINUE
ST7
WM=W4/3.0
ST7
,T W1
READ(
5
,1
5
0
)V 1
,VO,D W1
ST7
1
5
0 FORMAT(4F1
2
.0
)
ST7
WRITE(
6
,1
6
0
)VI,VO,DWI,T WI
ST7
1
6
0 FORMAT(!HO,2HV=El5.7
,2X,3HV0=El5.7
,2X,3HDW=El5.7
,2X, ST7
3HTW=El5.7
)
WSl=DWl-TWI
S T7 2
6
R=(DW1-TW!)/VO
S T7 2
7
V2=Vl-VO
8
S T7 2
)
*
1
0
0
.0
WS2=(Vト VO)/(DW1-TW1
9
S T7 2
WS3=WM-((V1-VO)/(DW1-TW1
)
*
1
0
0
.0
)
0
S T7 3
C=((VIVO)/(DW1-TW1
)
*
1
0
0
.O)*((DW1T Wl)/VO)
1
S T7 3
ALS=(l.0
(
1
0
0
.0/(((VI-VO)/(DW1-TWI
）
キ1
0
0
.O)*((DWIT W1
)/VO)+ S T7 3
2
S T7 3
）
紳
（
！
.0
/
3
.0
)
)
*
1
0
0
.0
3
1
0
0
.0）
GSI=I
.0/((DW1-TW!)/VO)(WM-((V1VO)/(DWト T W1
)
*
1
0
0
.0
)
)
/1
0
0
.0 S T7 3
4
S T7 3
.0/GSI
5
GS2=I
WS4=V0/(DWIT WI
)(
!
.O/((DW1-TWl)/VO
）
ー
（W M
ー
（
（Vト VO)/(DWト S T7 3
6
T W1
)
*
1
0
0
.0
)
)
/1
0
0
.0
)
S T7 3
7
ー
（
！
.O/((DWIT Wl
)/VO) (WM(
(
V
I VO)/(DWI- S T7 3
8
WS5=(VO/(DWト TWI）
）
パ0
0
.0
)
)
*
1
0
0
.0
9
T W1
)
*
1
0
0
.0）
S T7 3
WRITE(
6
,1
7
0
)V I
,VO,V2
,WS2
,WS3
,WS4
,WS5
S T7 4
0
1
1
7
0 FORMAT(IHO, 5X, 15HSHRINKAGELIMIT//,3X,2HV=El5.7
,3X, S T7 4
2
3HV0=E15.7
,3X,5HV-V0=E1
5
.7
,3X,21H(V-VO)*GAMW/WS*IOO.O S T7 4
=E1
5
.7/
IH,3X,3HWS=El5.7
,3X, 13Hl.O
/
R
1
.O/GS=E15.7
,3X,9HWS S T7 4
3
(
R
,GS)=E1
5
.7
)
WRITE(
6
,1
8
0
)VO,WSI,R
,WM,WS3,C
,ALS,GS!,GS2
ST7
1
8
0 FORMAT(
IHO,5X,15HSHRINKAGERATIO／
ん 3X,3HV0=El5.7
,3X, ST7
3HWS=El5.7
,3X,2HR=El5.7/IHO,5X,20HV0LUMETRICSHRINK- ST7
AGE／
ん 3X,3HWl=El5.7
,3X,3HWS=El5.7
,3X,2HC=El5.7/IHO,5X, ST7
!6HLINEARSHRINKAGE//,5X,3HLS=El5.7
/
1HO,5X, 16HSPECIFIC ST7
.O/R=WS/100.0=El5.7
,3X,3HGS=E1
5
.7
)
GRAVITY//,3X,15HI
1
0 CONTINUE
ST7
STOP
ST7
END
ST7
DATA
ST7
dTRJCO 四
oonwd
1
1 1 1／ 1
1
−
−
。
nu
Ln34せに1d
，
a anLqL
qLq4n
，
。
A T R d p O 勾J o o
d T d T 4 せ S T A守
，h
nwdnu’
n
ARidpO
4
5 5i5
5のO5A 5
55
$
3
.
0
1
1
.860
1
1
.7
4
5
1
1
.010
2
6
.
2
6
1
1
.1
9
ST7
ST7
ST7
ST7
ST7
J
1
7
.
2
9
0
1
7
.
0
0
5
1
5
.
8
5
0
1
3
.
0
勾
2
2
.
8
1
0
2
2
.
3
7
0
2
0
.
8
4
5
2
0
.
8
地
I
.0
2
.
0
48
佐賀大学農学業報第2
4号
表ー 1
2 収縮限界，収縮比，容積変化，線収縮，土粒子の比重の計算結果
料
試
I
縮
WS3
5
0
.
5
2
7
比
R
容積変化｜線収縮｜土粒子の比重
C
ALS
GS2
l
.
1
5
9
v
.
2
.7
9
8
あとがき
現在，土の物理試験で規定されている 8項目の物理試験について，電子計算機 HITAC5020E
（言語 FORTRANIV
）による計算処理のプログラムを組んで， 2
,3のデータの処理を行なった．
すべて加減乗除の計算であるから，特別の工夫を要することではないが，手動は一切行なわずに，
生の測定値から答を求めるプログラムであるから，とくに（5
）粒度試験などは長いプログラムに
なった．それでも処理時聞は 1
0
試料程度でも 1分間以内ですむようである．次報では土の力学試
験の電子計算機による計算処理について報告する予定である．