The part I'm still hesitant about is that the manifold you call C is n-dimensional. Codimension arguments in the infinite dimensional setting are always a little sticky. So, I'm just going to treat it as an assumption.
Then the answer to your question is yes, and doesn't depend on what type of infinite dimensional vector spaces you use as your local model. Let's see why.
First you have a space E which is a fiber bundle over B, via the projection p. Then at any point x in E we have the sequence of vector spaces,
TxF = ker dp --> Tx E ----> Tp(x)B
You also have a map f: E --> R, and hence a map df: TxE --> R = Tf(x)R for all points x in E.
C is defined to be the points in E such that the restriction of df to TxF vanishes. By assumption, C is an n-dimensional manifold.
Now let's talk about second derivatives. There are two problems with second derivatives. The first arrises because we are working in the infinite dimensional setting. This means that we might not be able to think of the second derivative as a bilinear form due to possible convergence problems. However, as you pointed out, we can think of it as a map,
TxE --> TxE*
Now depending on the type of infinite dimensional manifold you are considering, this dual space may take on different meanings (banach dual, Frechet dual, etc). You at least get something in the algebraic dual (of possible non-continuous functionals).
The second problem with second derivatives is that they are usually not defined independent of a choice of coordinates. Again, as you pointed out, the pair (df, second der of f) transforms as a section of the second jet bundle, i.e. the second derivative doesn't change like a tensor but in an affine fashion depending on df.
This means that in general there is no intrinsic way to say the second derivative is "non-degenerate" at a random point of E. However, along C, a portion of the second derivative is still well defined, independently of coordinate choices. It is not the whole second derivative, but just the composite,
H: TcE --> TcE* --> TcF*
(Here "H" is for Hessian or some such thing). Under coordinate changes, H transforms like a tensor because, along C, the restriction of df to TcF vanishes.
Here is the key fact: The tangent space of C is the kernel of H. This is an easy calculation, which can be done in local coordinates.
Finally, we are supposed look at the points of C and we are supposed to consider the restriction of H to TcF. We want to consider a point c in C, where the kernel of this restriction is zero.
But this condition is equivalent to saying that the intersection of TcF with the kernel of H is zero, i.e.
(**) At the point c, TcF has zero intersection with TcC.
Now since TcF is the kernel of dp, this means the the restriction of dp to TcC also has zero kernel, and hence for dimension reasons is an isomorphism. Hence at the point c in C, under your assumptions, the map f: C --> B is a local diffeomorphism.