Accurate modelling of rigid and soft inclusions in 2D elastic solids by the boundary element method

In this work, a particular sub-region technique that does not require traction approximations along interfaces is used to model rigid and soft inclusions. The stability and the accuracy of the formulation have been studied for several applications, including the degenerated cases where one dimension is strongly reduced to represent an embedded bar. The rigid inclusions have been modelled by increasing their elastic modulus. On the other hand thin inclusions with elastic modulus reduced to nearly zero can accurately model crack problems. In all analysed cases, the obtained results have shown that this simple strategy is reliable and can be applied to complex stress analysis problems.