The short answer is negative feedback - a very useful concept.
If you assume an ideal opamp with infinite gain (\$A_{OL}\$), infinite input impedances on the input pins, and zero output impedance; the concept of virtual ground or virtual short between the opamp input pins \$V_P\$ & \$V_N\$ can be used to show what is going on. Beware of the term virtual. Virtual short doesn't mean there's a physical short circuit between the positive and input pins, but looks like a short because of high open loop gain and negative feedback.
The output voltage is given by:
$$ V_O = A_{OL}(V_P - V_N)$$
When \$A_{OL} \to \infty\$, then \$(V_P - V_N) \to 0\$.
Since \$V_P = 0\$, \$V_N\$ will be zero when \$A_{OL} = \infty\$, i.e., a virtual short between \$V_P\$ and \$V_N\$
This says, assuming the opamp is not is saturation, the voltage at the negative input pin is zero volts for the ideal case.
This says that the current flowing through \$R_I\$ is the same magnitude of current flowing through \$R_F\$ which will keep the negative input node at zero volts. Thus, the current flowing through \$R_1\$ is zero, essentially, not a contributing part of the circuit.
Graphs are a good way to visualize what's going on. The opamp used in the simulation has an open loop gain (\$A_{OL}\$) of around 32 million which makes this nearly an ideal opamp at DC. Your simulation takes in to consideration non-ideal characteristics of the opamp: like the input bias currents which will affect the results.
The equation used for the upper graph is:
$$ {V_O \over V1} = {-R_F \over {{{R_I+R_F} \over A_{OL}}+R_I}} \tag{1}$$
The base equations for the equation 1 are:
$$ V_O = A_{OL}(V_P-V_N) \tag{2}$$
$$ I_I = I_F = {{V1 - V_N}\over R_I} = {{V_N-V_O}\over R_F} \tag{3}$$