  We present a model for the squeezed dark matter bispectrum, where the short
modes are deep in the non-linear regime. We exploit the consistency relations
for large-scale structures combined with a response function approach to write
the squeezed bispectrum in terms of a few unknown functions of the short modes.
We provide an ansatz for a fitting function for these response functions,
checking that the resulting model is reliable when compared to the one-loop
squeezed bispectrum. We then test the model against measured bispectra from
numerical simulations for short modes ranging between $k \sim 0.1 \, h/$Mpc,
and $k \sim 0.7 \, h/$Mpc at redshift $z=0$. To evaluate the goodness of the
fit of our model we implement a non-Gaussian covariance and find agreement
within $1$-$\sigma$ standard deviation of the simulated data.
