ex5_8.out

an introduction to boundary element methods一书源码

Example 5.8 (Discontinuous linear elements)
 
Input 
 
Number of Boundary Elements = 26 
 
Number of InteriorPoints = 6 
 
Coordinates of the Extreme Points  of the  Boundary Elements  
 
Point      X          Y
 1       0.0200     0.0000
 2       0.5000     0.0000
 3       1.0000     0.0000
 4       2.5000     0.0000
 5       4.0000     0.0000
 6       5.0000     0.0000
 7       7.5000     0.0000
 8       9.0000     0.0000
 9      10.0000     0.0000
10       9.2388     3.8268
11       7.0711     7.0711
12       3.8268     9.2388
13       0.0000    10.0000
14      -3.8268     9.2388
15      -7.0711     7.0711
16      -9.2388     3.8268
17     -10.0000     0.0000
18      -9.0000     0.0000
19      -7.5000     0.0000
20      -5.0000     0.0000
21      -4.0000     0.0000
22      -2.5000     0.0000
23      -1.0000     0.0000
24      -0.5000     0.0000
25      -0.0200     0.0000
26       0.0000     0.0200
 
Boundary  Conditions 

               First Node             Second Node
               Prescribed             Prescribed
Element     CODE   Value           CODE    Value
 1			 0   100.0000			 0   100.0000
 2			 0   100.0000			 0   100.0000
 3			 0   100.0000			 0   100.0000
 4			 0   100.0000			 0   100.0000
 5			 0   100.0000			 0   100.0000
 6			 0   100.0000			 0   100.0000
 7			 0   100.0000			 0   100.0000
 8			 0   100.0000			 0   100.0000
 9			 0   100.0000			 1     0.0000
10			 1     0.0000			 1     0.0000
11			 1     0.0000			 1     0.0000
12			 1     0.0000			 1     0.0000
13			 1     0.0000			 1     0.0000
14			 1     0.0000			 1     0.0000
15			 1     0.0000			 1     0.0000
16			 1     0.0000			 1     0.0000
17			 1     0.0000			 0     0.0000
18			 0     0.0000			 0     0.0000
19			 0     0.0000			 0     0.0000
20			 0     0.0000			 0     0.0000
21			 0     0.0000			 0     0.0000
22			 0     0.0000			 0     0.0000
23			 0     0.0000			 0     0.0000
24			 0     0.0000			 0     0.0000
25			 0     0.0000			 0     0.0000
26			 0     0.0000			 0     0.0000

Interior Point Coordinates

       Xi         Yi
    -4.9497     4.9497
    -2.2121     2.2121
     0.0000     3.0000
     0.0000     7.0000
     4.9497     4.9497
     2.2121     2.2121
 


Solution: 

Boundary Nodes 
          (X,Y)               u       Prenodal q  Postnodal q
(    0.0200,     0.0000)   100.0000  1317.7312  1317.7312
(    0.5000,     0.0000)   100.0000  -181.4681  -181.4681
(    1.0000,     0.0000)   100.0000    41.6826    41.6826
(    2.5000,     0.0000)   100.0000     6.1766     6.1766
(    4.0000,     0.0000)   100.0000     7.6256     7.6256
(    5.0000,     0.0000)   100.0000     5.6410     5.6410
(    7.5000,     0.0000)   100.0000     3.7712     3.7712
(    9.0000,     0.0000)   100.0000     3.3226     3.3226
(   10.0000,     0.0000)   100.0000     0.4944     0.4944
(    9.2388,     3.8268)    88.3650     0.0000     0.0000
(    7.0711,     7.0711)    76.0514     0.0000     0.0000
(    3.8268,     9.2388)    63.7051     0.0000     0.0000
(    0.0000,    10.0000)    51.2230     0.0000     0.0000
(   -3.8268,     9.2388)    38.6031     0.0000     0.0000
(   -7.0711,     7.0711)    25.8647     0.0000     0.0000
(   -9.2388,     3.8268)    13.0274     0.0000     0.0000
(  -10.0000,     0.0000)     1.3478     0.0000     0.0000
(   -9.0000,     0.0000)     0.0000    -5.0252    -5.0252
(   -7.5000,     0.0000)     0.0000    -4.3527    -4.3527
(   -5.0000,     0.0000)     0.0000    -6.7003    -6.7003
(   -4.0000,     0.0000)     0.0000    -8.4801    -8.4801
(   -2.5000,     0.0000)     0.0000   -12.6731   -12.6731
(   -1.0000,     0.0000)     0.0000   -33.2684   -33.2684
(   -0.5000,     0.0000)     0.0000   -87.2871   -87.2871
(   -0.0200,     0.0000)     0.0000  -226.3256  -226.3256
(    0.0000,     0.0200)     0.0000 -4947.9292 -4947.9292


Interior Points 
          (Xi,Yi)            u          q_x          q_y
(   -4.9497,     4.9497)    25.8225     3.2881      3.2935
(   -2.2121,     2.2121)    26.0081     7.5758      7.3495
(    0.0000,     3.0000)    51.6793    10.7238     -0.3675
(    0.0000,     7.0000)    51.1956     4.5655     -0.0131
(    4.9497,     4.9497)    75.9860     3.1497     -3.1411
(    2.2121,     2.2121)    76.4183     6.7149     -7.2595