I have raw data with three independent variables: temperature, stress and creep strain and one dependent variable: creep rate (see plot below). A colleague is doing some FEM where he needs to be able to get the values between the points in the plot, so I used the neural platform in JMP to get these for him.

This worked very well, I got very good agreement between the predictions from the neural model and and the data. I tried several models with different numbers of nodes in the hidden layer. (I don't have JMP PRO, so I can only use 1 layer.)
To decide which model (e.g. number of nodes N) to use, so I generate a plot of R2 and RMSE as a function of the number of nodes in the model. I expected that for the training set the R2 curve would rise and the RMSE curve would drop monotonously (with some random noise). For the validation set, I expected the R2 to rise and the RSME curve to drop until a certain value of N at which point the trends for the validation curve would be reversed. This would be the point at which the model would start over-fitting the training data and hence the performance on the validation data would drop.
What I got was this:

Initially, the trends are pretty much as I expected, but I have some questions:
1) In a model with ca. 10 nodes or more, the RMSE for the training set is systemically lower than for the validation set. Is that already a sign for over-fitting?
2) When using more than ca. 43 nodes (for the RMSE) or ca. 63 nodes (for the R2) the performance of the model drops for the validation and the training set. Why is that? My only idea is that this might be a problem of the number of points in the data sets not being sufficient for training the neural model properly. The total data set consists of 1600 points of which I used 2/3 for the training and 1/3 for the validation set.