There are TWO critical steps here!
1. The divergence of a general STRESS tensor (which must be symmetrical, by arguments given at first by Cauchy and others) appears in the PRIMITIVE equation of motion for the continuum.
2. We may "always" extract an isotropic pressure element out of this, by algebraic shuffling, which yields the the gradient of the pressure as appearing in our eqs of motion. Then, for the rest, we must MODEL a relation between our stresses, and our strains or strain rates. Simple, NEWTONIAN fluids are those saying that there exists a constant of proportionality between the strain rate tensor and the stress tensor, a more complicated model would be that of an anisotropic Newtonian fluid, in which the single constant of proportionality is replaced with a suitable constant tensor of "proportionality. Other models will be for "non-Newtonian" fluids, for example that the viscosity proportionality factor depends on temperature and pressure, thereby linking the thermodynamics of the fluid with its overall behaviour of motion. Or, we might have fluids where the viscosity is itself dependent on the strain rates, introducing an additional non-linearity. Or, we might have fluids in which not only strain rates, but also strains themselves generates stress; this makes typically, the fluid exhibit elastic properties.
--------------------------
The classical Navier-Stokes equations is gained as the simplest relationship, a single constant of proportionality between the stress tensor and the strain rate tensor is sufficient to describe the motion in the fluid.