AC-UNSAF3D

第四回地下水コロキュアム in OKAYAMA
操作の説明
実行時に必要なファイル
ACUNSF3D.EXE
WORK.LST
データファイル
結果の図化または取り出し
TRANS3D.EXE
GRAPH3D.EXE
すべて同じディレクトリに入れる！
WORK.LSTの内容
UNSAF使用するワークファイルのファイル名を指
定している。
第１ワークファイル名（A20）
第２ワークファイル名（A20）
第３ワークファイル名（A20）
必須
第４ワークファイル名（A20）
第５ワークファイル名（A20）
第６ワークファイル名（A20）（図化用出力）
実数配列数（I10）
整数配列数（I10）
入力データファイル名（A20）
出力データファイル名（A20）
プログラムとファイルの関係
work.lst
input.dat
file01.out file04.out
file02.out file05.out
file03.out file06.out
など
ACUNSF3D.exe
HTIME.TXT
print.prn
file06.out
GRAPH3D.exe
TRANS3D.exe
tecplot.txt
QTIME.TXT
（エクセル等で処理）
（印字イメージ・結果ファイル）
restart.dat
PTIME.TXT
femap.neu
avs.inp
output.txt
（結果ファイル）
エディター
TECPLOT
FEMAP
MicroAVS
表-2
MESH CHECK DIMENSION BD(
8002 )IB(
計算出力結果
必要メモリー(メッシュ情報読み取り時点)
11905)
要素構成情報
ELEMENT INFORMATION
要素番号
要素構成節点番号 1∼8
材質番号
M,
MTYP NODE1
NODE2
NODE3
NODE4
NODE5
NODE6
NODE7
NODE8
1
1
122
133
134
123
1
12
13
2
2
3
1
1
123
124
134
135
135
136
124
125
2
3
13
14
14
15
3
4
4
5
6
1
1
1
125
126
127
136
137
138
137
138
139
126
127
128
4
5
6
15
16
17
16
17
18
5
6
7
7
1
128
139
140
129
7
18
19
8
・
・
996
997
1
1
1315
1316
1326
1327
1327
1328
・
1316
1317
1194
1195
1205
1206
1206
1207
1195
1196
998
999
1000
1
1
1
1317
1318
1319
1328
1329
1330
1329
1330
1331
1318
1319
1320
1196
1197
1198
1207
1208
1209
1208
1209
1210
1197
1198
1199
UNITS ARE CONTROLED BY THE FOLLOWING CONSTANTS 基本単位
CONS1 =
1.000E+00
MATERIAL
K1
CONS2 =
1.000E+00
K2
K3
POPOSITY
SS
飽和材質情報
1
.10000E-05
.10000E-05
.10000E-05
.30000E+00
TABLE OF RELATIVE CONDUCTIVITY VERSUS MOISTURE CONTENT FOR MATERIAL
RELA-CON
MOIS-CON
RELA-CON
MOIS-CON
RELA-CON
MOIS-CON
.10000E-05
1
RELA-CON
MOIS-CON
RELA-CON
MOIS-CON
Kr-θ
.0000E+00
.0000E+00
.1000E+01
.3000E+00
TABLE OF PRESSURE HEAD VERSUS MOISTURE CONTENT FOR MATERIAL 1
PRES-HE
MOIS-CON
PRES-HE
MOIS-CON
PRES-HE
MOIS-CON
-.2000E+02
.0000E+00 -.1000E+02
.1500E+00
.0000E+00
PRES-HE
.1500E-01
MOIS-CON
.7500E-01
SPECI-CA
.1500E-01
MOIS-CON
.2250E+00
SPECI-CA
.0000E+00
PRES-HE
MOIS-CON
ψ-θ
.3000E+00
TABLE OF SPECIFIC MOISTURE CAPACITY VERSUS MOISTURE CONTENT FOR MATERIAL
SPECI-CA
MOIS-CON
MOIS-CON
.0000E+00
1
SPECI-CA
MOIS-CON
SPECI-CA
MOIS-CON
c-θ
節点情報
節点番号
KODE
節点流量
座標(x,y,z)
NODAL POINT INFORMATION
NODE NO.
圧力水頭
境界条件番号
-X-
-Y-
-Z-
.PSI.
Q
1
2
.000000E+00
.000000E+00
.000000E+00
.000000E+00
.000000E+00
2
2
.000000E+00
.100000E+01
.000000E+00
.000000E+00
.000000E+00
3
2
.000000E+00
.200000E+01
.000000E+00
.000000E+00
.000000E+00
4
5
2
2
.000000E+00
.000000E+00
.300000E+01
.400000E+01
.000000E+00
.000000E+00
.000000E+00
.000000E+00
.000000E+00
.000000E+00
6
2
.000000E+00
.500000E+01
.000000E+00
.000000E+00
.000000E+00
7
2
.000000E+00
.600000E+01
.000000E+00
.000000E+00
.000000E+00
8
2
.000000E+00
.700000E+01
.000000E+00
.000000E+00
.000000E+00
・
・
・
1324
1
.100000E+02
.300000E+01
.100000E+02
.000000E+00
.000000E+00
1325
1
.100000E+02
.400000E+01
.100000E+02
.000000E+00
.000000E+00
1326
1
.100000E+02
.500000E+01
.100000E+02
.000000E+00
.000000E+00
1327
1328
1
1
.100000E+02
.100000E+02
.600000E+01
.700000E+01
.100000E+02
.100000E+02
.000000E+00
.000000E+00
.000000E+00
.000000E+00
1329
1
.100000E+02
.800000E+01
.100000E+02
.000000E+00
.000000E+00
1330
1
.100000E+02
.900000E+01
.100000E+02
.000000E+00
.000000E+00
1331
1
.100000E+02
.100000E+02
.100000E+02
.000000E+00
.000000E+00
MESH CHECK DIMENSION BD(
10002 )IB(
11905)
必要メモリー
(一次元配列透水マトリックス必要容量検索時点)
MESH CHECK DIMENSION BD(
10002 )IB(
必要メモリー
41696)
(一次元配列透水マトリックス格納時点)
TOTAL DIMENSION
BD(
79675 )
総必要メモリー(倍精度実数配列)
TOTAL DIMENSION
CD(
44359 )
総必要メモリー(4 バイト整数配列)
非定常問題を選択
0***** UNSTEADY PROBLEM *******
HEAD タイトル
0---- DEMO DATA 01 ---NUMBER OF NODAL POINTS-------
1331
MAX NODES IN A ELEMENT-------
8
NUMBER OF ELEMENTS-----------
1000
NUMBER OF MATERIALS----------
1
MAX NUMBER OF ITERATIONS-----
5
MAX NUMBER OF ITERATIONS(T=0)
0
MAX NUMBER OF CG ITERATIONS--
2662
主な設定条件の出力
NUMBER OF MATRIX DIMENSIONS-- 29791
計算開始以下タイムステップ毎の結果が出力される
---- DEMO DATA 01 ---TIME =
1.00000E-01
HEAD タイトル
TIME STEP =
時間，タイムステップ数
1
(第一ステップは常に出力)
*** SUCCESS TO CONVERGENCE IN SOLVER -- IT =
*** SUCCESS TO CONVERGENCE IN SOLVER -- IT =
14 EPSX =
13 EPSX =
.1955855E-10
.8915536E-11
繰り返し計算毎の PCG 解法収束状況
*** SUCCESS TO CONVERGENCE IN SOLVER -- IT =
1 EPSX =
.2111744E-13
繰り返し回数 IT，収束値 EPSX
MAXIMUM CHANGE IN PRESSURE HEAD DURING ITERATION 3 WAS 7.10543E-15 AT NODE
非線形繰り返しステップと最大残差値
CUMULATIVE INFLOW INTO SYSTEM IS = -4.59962E-05
全水頭値
0 NODE
HEAD
圧力水頭値
PRESS HEAD
76
最大残差節点番号
このステップの時間増分間での変動水量収支
流量
DISCHARGE NODE
HEAD
PRESS HEAD
DISCHARGE NODE
HEAD
PRESS HEAD
DISCHARGE
節点番号
1 0.0000E+00 0.0000E+00 -2.3926E-06
2 0.0000E+00 0.0000E+00 -4.7852E-06
3 0.0000E+00 0.0000E+00 -4.7852E-06
4 0.0000E+00 0.0000E+00 -4.7852E-06
5 0.0000E+00 0.0000E+00 -4.7852E-06
6 0.0000E+00 0.0000E+00 -4.7852E-06
7 0.0000E+00 0.0000E+00 -4.7852E-06
8 0.0000E+00 0.0000E+00 -4.7852E-06
9 0.0000E+00 0.0000E+00 -4.7852E-06
10 0.0000E+00 0.0000E+00 -4.7852E-06
13 9.2186E+00 9.2186E+00 0.0000E+00
11 0.0000E+00 0.0000E+00 -2.3926E-06
14 9.2186E+00 9.2186E+00 0.0000E+00
12 9.2186E+00 9.2186E+00 0.0000E+00
15 9.2186E+00 9.2186E+00 0.0000E+00
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
1318 1.0000E+01 -4.1854E-10 0.0000E+00 1319 1.0000E+01 -4.1854E-10 0.0000E+00 1320 1.0000E+01 -4.1854E-10 0.0000E+00
1321 1.0000E+01 0.0000E+00 1.0984E-16 1322 1.0000E+01 0.0000E+00 2.1968E-16 1323 1.0000E+01 0.0000E+00 2.1968E-16
1324 1.0000E+01 0.0000E+00 2.1968E-16 1325 1.0000E+01 0.0000E+00 2.1968E-16 1326 1.0000E+01 0.0000E+00 2.1968E-16
1327 1.0000E+01 0.0000E+00 2.1968E-16 1328 1.0000E+01 0.0000E+00 2.1968E-16 1329 1.0000E+01 0.0000E+00 2.1968E-16
1330 1.0000E+01 0.0000E+00 2.1968E-16 1331 1.0000E+01 0.0000E+00 1.0984E-16
***
QP( + ) = 1.02561E-13
解析領域への総流入流量(>0)
***
QM( - ) = -4.59962E-04
解析領域から総排出流量(<0)
0 0.0000E+00 0.0000E+00 0.0000E+00
** * S U CC E SS TO CO N VE R GE NC E I N S O LV E R - - I T =
1 EP S X =
M A XI M UM CH A NG E I N P R ES SU R E H EA D D U RI N G I TE R AT I ON
3 WA S
. 6 62 4 52 3 E- 1 1
5 .4 0 18 4 E- 0 8 A T N OD E
1 21 1
C U MU L AT I VE IN F LO W I N TO S Y ST E M I S = - 2 .5 2 86 5 E- 0 4
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
・
* * ** * *G E T G OO D R E SU L T* ** * ** *
* ** G OO D , G OO D ** * *
計算正常終了
第四回地下水コロキュアム in OKAYAMA
演習問題の説明
AC-UNSAF3D演習
500ｍ
トンネル
500ｍ
掘削
20ｍ
揚水
30ｍ
15ｍ
20ｍ
斜面
50ｍ
50ｍ
5ｍ
（要素番号）＝（最上部の要素番号）＋４００×（（層番号）−１）
1層目
１
２層目
401
３層目
801
（節点番号）＝（最上部の節点番号）＋４４１×（（層番号）−１）
1層目
１
２層目
442
883
1324
３層目
４層目
トンネル
側面境界
（水位一定）
側面境界
（水位一定）
地表面
(降雨境界)
トンネル
（浸出点）
初期水位地表面に水面
降雨強度 1.6E-8m/sec≒505mm/year
V1
G1
10m
（10要素）
材質２
6m
（6要素）
材質１
４m
（４要素）
Z
Y
9要素
X
2要素
掘削
側面境界
（水位一定）
側面境界
（水位一定）
地表面
(降雨境界)
掘削面
（浸出点）
遮水壁
初期水位地表面に水面
降雨強度 1.6E-8m/sec≒505mm/year
V1
G1
材質２
10m （10要素）
材質１
５m （５要素）
５m （５要素）
Z
Y
X
材質３
掘削範囲４×９要素
斜面
地表面
(降雨境界)
水位一定境界
揚水
側面境界
（水位一定）
揚水井
（水位一定）
側面境界
（水位一定）
９要素
３要素
GL-14mの位置の節点を圧力水頭＝０で水位固定
トンネル全水頭コンター図
V1
G1
20.
19.
18.
17.
16.
15.
14.
13.
12.
11.
10.
9.
8.
7.
Z
Y
X
Output Set: 0 T= 0.0000E+00
Contour: Total Head
6.
5.
4.
掘削全水頭コンター図
V1
G1
20.
19.37
18.74
18.11
17.48
16.85
16.22
15.59
14.96
14.33
13.7
13.07
12.44
11.81
Z
Y
X
Output Set: 0 T= 0.0000E+00
Contour: Total Head
11.17
10.54
9.914
斜面全水頭コンター図
V1
30.
28.44
26.88
25.31
23.75
22.19
20.63
19.06
17.5
15.94
14.38
12.81
11.25
9.688
Z
Y
8.125
X
Output Set: 0 T= 0.0000E+00
Contour: Total Head
6.563
5.
揚水全水頭コンター図
V1
20.
19.13
18.25
17.38
16.5
15.63
14.75
13.88
13.
12.13
11.25
10.38
9.5
8.625
Z
Y
7.75
X
Output Set: 0 T= 0.0000E+00
Contour: Total Head
6.875
6.

Download Report