A new enrichment strategy with high convergence rate and low condition number. Applying enrichment strategy to model stress singularity at multi-material junctions. Modeling stress singularity at the crack tip terminating by bimaterial interface. Applying the eigen-function expansion method to derive enrichment functions. Calculation of the generalized SIFs using direct and interaction integral methods.