Sorry for being imprecise. I am saying that in the context of a nonlinear model Dynare treats expressions that involve something like f( c(+1), r(+1) ) as E_t[ f( c(+1), r(+1) ) ]. Essentially, you can think of an conditional expectation around both sides of every equation. However, as is standard in perturbation techniques (see Schmitt-Grohé/Uribe (2004)
http://www.columbia.edu/~mu2166/2nd_order.htm) the function f (and the policy functions) are approximated with a Taylor polynomial inside of the expected value. Hence, up to first order E [ f( c(+1), r(+1) ) ] is the same as f( E[c(+1)], E[r(+1)] ) because the system is linear and there is no Jensen's Inequality effect. For higher order approximations, see the linked article.
If you search the forum, you will find several answers regarding this issue. Unfortunately, I am not aware of any official documentation regarding this issue. This convention is rather taken for granted.