ex5_7.out

an introduction to boundary element methods一书源码

Example 5.7 (Linear Discontinuous Elements)
 
Input 
 
Number of Boundary Elements = 16 
 
Number of InteriorPoints = 5 
 
Coordinates of the Extreme Points  of the  Boundary Elements  
 
Point      X          Y
 1       0.0010     0.0000
 2       2.0000     0.0000
 3       4.0000     0.0000
 4       5.9990     0.0000
 5       6.0000     0.0010
 6       6.0000     2.0000
 7       6.0000     4.0000
 8       6.0000     5.9990
 9       5.9990     6.0000
10       4.0000     6.0000
11       2.0000     6.0000
12       0.0010     6.0000
13       0.0000     0.5990
14       0.0000     4.0000
15       0.0000     2.0000
16       0.0000     0.0010
 
Boundary  Conditions 

               First Node             Second Node
               Prescribed             Prescribed
Element     CODE   Value           CODE    Value
 1			 0     0.0000			 0     0.0000
 2			 0     0.0000			 0     0.0000
 3			 0     0.0000			 0     0.0000
 4			 0     0.0000			 1    46.1900
 5			 1    46.1900			 1    46.1900
 6			 1    46.1900			 1    46.1900
 7			 1    46.1900			 1    46.1900
 8			 1    46.1900			 0     4.2000
 9			 0     4.2000			 0    84.0000
10			 0    84.0000			 0   168.0000
11			 0   168.0000			 0   247.8000
12			 0   247.8000			 1   -46.1900
13			 1   -46.1900			 1   -46.1900
14			 1   -46.1900			 1   -46.1900
15			 1   -46.1900			 1   -46.1900
16			 1   -46.1900			 0     0.0000

Interior Point Coordinates

       Xi         Yi
     2.0000     2.0000
     2.0000     4.0000
     3.0000     3.0000
     4.0000     2.0000
     4.0000     4.0000
 


Solution: 

Boundary Nodes 
          (X,Y)               u       Prenodal q  Postnodal q
(    0.0010,     0.0000)     0.0000   156.5300   156.5300
(    2.0000,     0.0000)     0.0000   -16.8707   -16.8707
(    4.0000,     0.0000)     0.0000   -17.7104   -17.7104
(    5.9990,     0.0000)     0.0000  -103.2299  -103.2299
(    6.0000,     0.0010)     6.1780    46.1900    46.1900
(    6.0000,     2.0000)   114.0960    46.1900    46.1900
(    6.0000,     4.0000)   131.4528    46.1900    46.1900
(    6.0000,     5.9990)    12.4986    46.1900    46.1900
(    5.9990,     6.0000)     4.2000  -151.6351  -151.6351
(    4.0000,     6.0000)    84.0000    11.2954    11.2954
(    2.0000,     6.0000)   168.0000    44.0099    44.0099
(    0.0010,     6.0000)   247.8000    71.9776    71.9776
(    0.0000,     0.5990)   -89.3071   -46.1900   -46.1900
(    0.0000,     4.0000)    52.9451   -46.1900   -46.1900
(    0.0000,     2.0000)  -116.3532   -46.1900   -46.1900
(    0.0000,     0.0010)    -8.5014   -46.1900   -46.1900


Interior Points 
          (Xi,Yi)            u          q_x          q_y
(    2.0000,     2.0000)    16.6571    25.3385     20.5396
(    2.0000,     4.0000)    84.4683    -1.0487     41.9870
(    3.0000,     3.0000)    59.1773    12.0620     23.8646
(    4.0000,     2.0000)    52.5296    17.4905     22.7641
(    4.0000,     4.0000)    83.9398     4.0897      7.2784