The purpose of this work is to identify simultaneously two irregular interfacial boundaries configurations between the components of three connected domains using a discrete number of displacement measurements obtained by a simple tension test. A unique combination of a global optimization method, i.e., the imperialist competitive algorithm and two local optimization methods, i.e., the conjugate gradient method and Simplex method along with the inverse application of the boundary elements method are employed in an inverse software package. A fitness function, which is the summation of squared differences between calculated displacements and measured displacements at identical locations on the outer boundary, is minimized. Due to the complexity and the ill-posed nature of this identification problem, imperialist competitive algorithm is used to find the best initial guesses of the unknown interface boundaries in order to be used by the local optimization techniques, i.e., conjugate gradient method and then Simplex method to accurately converge to the optimal shape of two irregular interfacial boundaries between the components of an inhomogeneous body. Several examples are selected, and the accuracy of the obtained results is discussed. The effect of experimental measurement errors and the influence of material properties of the sub regions on the identification process are also investigated.