Body Forces - The Boundary Element Method with Programming - page 374

would have to cover the entire domain, including the part where no swelling occurs.
Also note that the number of unknowns at the nodes of the boundary elements was not
increased by the fact that cells were required to compute the right hand side of the
system of equations.
The implementation of body forces requires additional fundamental solutions which
have been added to the library. It is obvious that volumetric loading effects also occur in
potential problems but the discussion of these in more detail is beyond the scope of this
topic.
We have only shown one example of implementation: the treatment of initial strains
as they may occur in problems where part of the domain, is subjected to swelling.
However with the knowledge of programming gained in previous chapters and the
explanation of the theory in this chapter the reader should be able to perform the
implementation of the other types of body forces into program General_purpose_BEM.

13.11 EXERCISES

Exercise 13.1
Write a Subroutine similar to SUBROUTINE Body_force that computes the right hand
side for gravity as discussed in section 13.2. Implement this into program
General_Purpose_BEM and test on an example of a cube subjected to self weight.

1, 0 N

1, 0 m

1, 0 N

1, 0 m

Note: Apply this
load as traction

Figure 13.10 Test example for exercise 13.1