Numerical Modeling of Fracture Height Propagation in Multilayer Formations Considering the Plastic Zone and Induced Stress