ex5_5.out

an introduction to boundary element methods一书源码

Example 5.5 (Linear Elements)
 
Input 
 
Number of Boundary Elements = 12 
 
Number of InteriorPoints = 5 
 
Coordinates of the Extreme Points  of the  Boundary Elements  
 
Point      X          Y
 1       0.0000     0.0000
 2       2.0000     0.0000
 3       4.0000     0.0000
 4       6.0000     0.0000
 5       6.0000     2.0000
 6       6.0000     4.0000
 7       6.0000     6.0000
 8       4.0000     6.0000
 9       2.0000     6.0000
10       0.0000     6.0000
11       0.0000     4.0000
12       0.0000     2.0000
 
Boundary  Conditions 

               First Node             Second Node
               Prescribed             Prescribed
Element     CODE   Value           CODE    Value
 1			 1     0.0000			 1     0.0000
 2			 1     0.0000			 1     0.0000
 3			 1     0.0000			 1     0.0000
 4			 0   200.0000			 0   200.0000
 5			 0   200.0000			 0   200.0000
 6			 0   200.0000			 0   200.0000
 7			 1     0.0000			 1     0.0000
 8			 1     0.0000			 1     0.0000
 9			 1     0.0000			 1     0.0000
10			 0     0.0000			 0     0.0000
11			 0     0.0000			 0     0.0000
12			 0     0.0000			 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.0000,     0.0000)     0.0000   -33.3333     0.0000
(    2.0000,     0.0000)    66.6666     0.0000     0.0000
(    4.0000,     0.0000)   133.3335     0.0000     0.0000
(    6.0000,     0.0000)   200.0000     0.0000    33.3333
(    6.0000,     2.0000)   200.0000    33.3334    33.3334
(    6.0000,     4.0000)   200.0000    33.3334    33.3334
(    6.0000,     6.0000)   200.0000    33.3333     0.0000
(    4.0000,     6.0000)   133.3335     0.0000     0.0000
(    2.0000,     6.0000)    66.6665     0.0000     0.0000
(    0.0000,     6.0000)     0.0000     0.0000   -33.3333
(    0.0000,     4.0000)     0.0000   -33.3334   -33.3334
(    0.0000,     2.0000)     0.0000   -33.3334   -33.3334


Interior Points 
          (Xi,Yi)            u          q_x          q_y
(    2.0000,     2.0000)    66.6667    33.3332     -0.0002
(    2.0000,     4.0000)    66.6667    33.3332      0.0002
(    3.0000,     3.0000)   100.0000    33.3333      0.0000
(    4.0000,     2.0000)   133.3339    33.3339     -0.0005
(    4.0000,     4.0000)   133.3339    33.3339      0.0005