Characteristic features of extracellular biological fluids with complicated inner structure require a special approach to their modeling. This approach is demonstrated on the example of blood flow in large blood vessels. The non-Newtonian liquid and the semi-slip boundary conditions were proposed for modelling the blood pulse flow. As results of numerical calculations the pressure distribution in the blood flow and the stress distribution in the deformable blood vessel were obtained.