Abstract: A set of non-linear finite element formulae with incremental form was derived on the basis of the fluid-solid coupling equations of porous fluid diffusion and deformation. The pore pressure cohesive elements were used to simulate the initiation and propagation of hydraulic pre-defined fractures resulted from hydraulic loading. A finite element model was proposed for hydraulic fracturing of oil well.Both the viscosity of fracturing fluid and volume concentration of proppant varied with time in the model. All of the simulation parameters are from data of oil wells. Some parameters including pressure at the fracture mouth, perforation friction, string friction, hydrostatic pressure of fluid column in well bore and wellhead pressure have changed with fracture time. The presented "curves of stimulation treatment history" are in good coincidence with those obtained from the field measurements.