A new computational approach that combines the extended finite element method associated with variable-node elements and cohesive zone model is developed. By using a new enriched technique based on sign function, the proposed model using 4-node quadrilateral elements can eliminate the blending element problem. It also allows modeling the equal stresses at both sides of the crack in the crack-tip as assumed in the cohesive model, and is able to simulate the arbitrary crack-tip location. The multiscale mesh technique associated with variable-node elements and the arc-length method further improve the efficiency of the developed approach. The performance and accuracy of the present approach are illustrated through numerical experiments considering both mode-I and mixed-mode fracture in concrete.