WELCOME TO PARABOLIC ELEMENT BEM PROGRAM
I N P U T D A T A
BOUNDARY ELEMENT METHOD APPLIED TO POTENTIAL PROBLEMS
BOUNDARY CONDITIONS ARE ASSIGNED TO ELEMENTS.
IF ELM TYPE(I) = 1, THEN DIRICHLIT.
IF ELM TYPE(I) = 2, THEN NEUMANN.
====> NUMBER OF ELEMENTS= 13
ELM # ---> ELEMENT NUMBER
ELM # I J K ELM TYPE BOUNDARY VALUES
1 1 2 3 2 0.00000 0.00000 0.00000
2 3 4 5 2 0.00000 0.00000 0.00000
3 5 6 7 2 0.00000 0.00000 0.00000
4 7 8 9 2 0.00000 0.00000 0.00000
5 9 10 11 2 0.00000 0.00000 0.00000
6 11 12 13 1 10.00000 10.00000 10.00000
7 13 14 15 1 10.00000 10.00000 10.00000
8 15 16 17 2 0.00000 0.00000 0.00000
9 17 18 19 2 0.00000 0.00000 0.00000
10 19 20 21 1 3.00000 3.00000 3.00000
11 21 22 23 1 3.00000 4.16660 5.33330
12 23 24 25 1 5.33300 6.50000 7.66660
13 25 26 1 1 7.66660 8.83333 10.00000
====> NUMBER OF BOUNDARY NODES= 26
NODE X-COORDINATE Y-COORDINATE
1 10.00000000 10.00000000
2 9.00000000 10.00000000
3 8.00000000 10.00000000
4 7.00000000 10.00000000
5 6.00000000 10.00000000
6 5.00000000 10.00000000
7 4.00000000 10.00000000
8 3.00000000 10.00000000
9 2.00000000 10.00000000
10 1.00000000 10.00000000
11 0.00000000 10.00000000
12 0.00000000 7.50000000
13 0.00000000 5.00000000
14 0.00000000 2.50000000
15 0.00000000 0.00000000
16 2.50000000 0.00000000
17 5.00000000 0.00000000
18 7.50000000 0.00000000
19 10.00000000 0.00000000
20 10.00000000 1.50000000
21 10.00000000 3.00000000
22 10.00000000 4.16666000
23 10.00000000 5.33333000
24 10.00000000 6.50000000
25 10.00000000 7.66666000
26 10.00000000 8.83333000
====> NUMBER OF NODES ON FREE SURFACE = 11
NODAL NUMBER
1 2 3 4 5 6 7 8 9 10 11
====> NUMBER OF NODES ON SEEPAGE FACE = 7
NODAL NUMBER
21 22 23 24 25 26 1
END OF INPUT-DATA ECHO PRINT
*** N U M E R I C A L S O L U T I O N ***
-- FLUX ON BOUNDARY -- (QN = HEAT FLUX IN NORMAL DIRECTION)
ELEMENT NODE-I NODE-J NODE-K QN(I) QN(J) QN(K)
1 1 2 3 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
2 3 4 5 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
3 5 6 7 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
4 7 8 9 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
5 9 10 11 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
6 11 12 13 -0.2050216200 -0.3901906321 -0.4784522229
7 13 14 15 -0.4784522229 -0.5443844408 -0.5688938182
8 15 16 17 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
9 17 18 19 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
10 19 20 21 1.006694331 0.9372603091 1.896682334
11 21 22 23 1.896682334 1.283239982 1.097676341
12 23 24 25 1.097676341 0.9205527644 0.7462903332
13 25 26 1 0.7462903332 0.6380661255 -0.8189196023E-02
-- FREE TERM AND POTENTIAL VALUES AT NODAL POINTS --
NODE FREE TERM POTENTIAL
1 0.416188 4.3027969
2 0.500000 4.7286685
3 0.489084 5.3685438
4 0.500000 6.0005474
5 0.499132 6.6484857
6 0.500000 7.2879760
7 0.500009 7.9138979
8 0.500000 8.5154806
9 0.500864 9.0829588
10 0.499999 9.5915969
11 0.212718 10.000000
12 0.500000 10.000000
13 0.500000 10.000000
14 0.500000 10.000000
15 0.250000 10.000000
16 0.500000 8.5630088
17 0.500000 7.0057208
18 0.500000 5.1822220
19 0.250000 3.0000000
20 0.500000 3.0000000
21 0.500000 3.0000000
22 0.500000 3.2171328
23 0.500000 3.4342656
24 0.500000 3.6513984
25 0.500000 3.8685312
26 0.500000 4.0856641
COORDINATES OF SEEPAGE FACE
X= 10.0000000000000 Y= 3.00000000000000
X= 10.0000000000000 Y= 3.21713281206609
X= 10.0000000000000 Y= 3.43426562413219
X= 10.0000000000000 Y= 3.65139843619828
X= 10.0000000000000 Y= 3.86853124826437
X= 10.0000000000000 Y= 4.08566406033047
X= 10.0000000000000 Y= 4.30279687239656
COORDINATES OF FREESURFACE
X= 10.0000000000000 Y= 4.30279687239656
X= 9.68377223398316 Y= 4.72866847047851
X= 9.10557280900008 Y= 5.36854377025939
X= 8.35683232748450 Y= 6.00054741891475
X= 7.47017787186530 Y= 6.64848574594796
X= 6.46446609406726 Y= 7.28797599162892
X= 5.35241998455110 Y= 7.91389792879939
X= 4.14337981426147 Y= 8.51548055777747
X= 2.84458247200067 Y= 9.08295875353491
X= 1.46185031754538 Y= 9.59159685391281
X= 0.000000000000000E+000 Y= 10.0000000000000