The initial stages of creep in nickel base superalloy deserve a special attention since they prevail in most of the service live of critical jet engine components. In nickel base superalloys, a very sharp minimum in creep rate is observed at strains less than 0.5% in a wide range of temperatures: the lower the applied stress the steeper the minimum, often referred to as incubation period. In this paper, the primary creep stage of Ni superalloy is investigated using 3D discrete dislocation dynamics modeling. The simulated volume corresponds to a 3D periodic cell of the cubic arrangement of γ−γ′ NiAl loaded along the (001) direction. Misfit stresses induced by the lattice mismatch (δ = -3.10−3) between the precipitate and the matrix are taken into account through a coupling of the dislocation dynamics code with a finite element code. Frank-Read sources are positioned in the horizontal channel. It is shown that most of the dislocations propagate in the horizontal channels by repeated cross-slip on activated {111} glide planes and a few invade the vertical channels. This behavior is rationalized by the model. In the γ-γ′ interfaces, the zigzagging lines of dislocations straighten out into [100] and [010] directions leading to a dislocation microstructure similar to the ideal misfit dislocation network needed to compensate the lattice mismatch. Calculations of internal stresses in the channels show that this dislocation microstructure induces biaxial planar stresses in the horizontal γ channels leading to an elastic contraction of the sample along the loading direction. Three-dimensional DDD simulations show that the plastic strain induced, in the loading direction, by the movement and multiplication of dislocations is less important than the elastic relaxation it produces along that same direction. Hence, for very low applied tensile stresses, Ni superalloys display an abnormal behavior often referred to as negative creep that can be evidenced experimentally, even in polycrystalline materials. This last observation clearly implies that the first active batches of dislocations, during creep, or static high temperature ageing are mobilized to cancel out the numerous and pervading internal elastic strain fields generated within the crystals by the γ’ precipitates due to their misfit with the γ matrix.