Unable to find the steady state

  1. Only do exp()-substitutions once the model runs. Even better, use an alternative. See Question about understanding irfs in dynare - #4 by jpfeifer
  2. It seems you are setting up stochastic simulations, but with inequality constraints. That won’t work.