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= 11
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 1 20.00000 20.00000 20.00000
4 7 8 9 1 20.00000 20.00000 20.00000
5 9 10 11 2 0.00000 0.00000 0.00000
6 11 12 13 2 0.00000 0.00000 0.00000
7 13 14 15 1 10.00000 10.00000 10.00000
8 15 16 17 1 10.00000 11.25000 12.50000
9 17 18 19 1 12.50000 13.75000 15.00000
10 19 20 21 1 15.00000 16.25000 17.50000
11 21 22 1 1 17.50000 18.75000 20.00000
====> NUMBER OF BOUNDARY NODES= 22
NODE X-COORDINATE Y-COORDINATE
1 40.00000000 20.00000000
2 35.00000000 20.00000000
3 30.00000000 20.00000000
4 25.00000000 20.00000000
5 20.00000000 20.00000000
6 15.00000000 15.00000000
7 10.00000000 10.00000000
8 5.00000000 5.00000000
9 0.00000000 0.00000000
10 15.00000000 0.00000000
11 30.00000000 0.00000000
12 45.00000000 0.00000000
13 60.00000000 0.00000000
14 55.00000000 5.00000000
15 50.00000000 10.00000000
16 48.75000000 11.25000000
17 47.50000000 12.50000000
18 46.25000000 13.75000000
19 45.00000000 15.00000000
20 43.75000000 16.25000000
21 42.50000000 17.50000000
22 41.25000000 18.75000000
====> NUMBER OF NODES ON FREE SURFACE = 5
NODAL NUMBER
1 2 3 4 5
====> NUMBER OF NODES ON SEEPAGE FACE = 9
NODAL NUMBER
15 16 17 18 19 20 21 22 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.4629375508 -0.1812527288 -0.1357510541
4 7 8 9 -0.1357510541 -0.6192324992E-01 -0.2004327095E-02
5 9 10 11 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
6 11 12 13 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
7 13 14 15 0.1613558439E-01 0.1653745196 0.8097661705
8 15 16 17 0.8097661705 0.6407300384 0.5220634462
9 17 18 19 0.5220634462 0.4589920030 0.4002007453
10 19 20 21 0.4002007453 0.3457503392 0.2789858457
11 21 22 1 0.2789858457 0.2247377878 -0.6530291043E-03
-- FREE TERM AND POTENTIAL VALUES AT NODAL POINTS --
NODE FREE TERM POTENTIAL
1 0.470006 10.860264
2 0.500000 12.514603
3 0.496605 14.570977
4 0.500000 16.921124
5 0.326739 20.000000
6 0.500000 20.000000
7 0.500000 20.000000
8 0.500193 20.000000
9 0.125000 20.000000
10 0.500000 19.036930
11 0.500000 16.302199
12 0.500000 12.381187
13 0.125000 10.000000
14 0.500193 10.000000
15 0.500000 10.000000
16 0.500000 10.107533
17 0.500000 10.215066
18 0.500000 10.322599
19 0.499999 10.430132
20 0.499991 10.537665
21 0.499929 10.645198
22 0.499165 10.752731
COORDINATES OF SEEPAGE FACE
X= 50.0000000000000 Y= 10.0000000000000
X= 49.8924670128289 Y= 10.1075329871711
X= 49.7849340256577 Y= 10.2150659743423
X= 49.6774010384866 Y= 10.3225989615134
X= 49.5698680513155 Y= 10.4301319486845
X= 49.4623350641443 Y= 10.5376649358557
X= 49.3548020769732 Y= 10.6451979230268
X= 49.2472690898020 Y= 10.7527309101980
X= 49.1397361026309 Y= 10.8602638973691
COORDINATES OF FREESURFACE
X= 49.1397361026309 Y= 10.8602638973691
X= 45.4972690898020 Y= 12.5146027077983
X= 38.8372836025525 Y= 14.5709770646058
X= 30.2129223092912 Y= 16.9211242870010
X= 20.0000000000000 Y= 20.0000000000000