Usually Delunay triangulisation is sufficient for elements having one or more inclusions or either a crack or a tip element. But if bot the crack tip and inclusion are their triangulisation does not give correct stress distribution if we include both crack tip and inclusion data. I am coding in matlab.