Deliyannis-type narrowband filter is oscillating when I apply GND to the input
Regarding the shown circuit some comments are necessary. Without re-calculating the parts values (not possible because of uncomplete filter specification) I did some simulation runs based on your parts values.

Results: 

Opamp LF356: Midfrequency fo=67 kHz. 

Opamp LM741: Midfrequency fo=39 kHz.

What does this mean? The midfrequency of the bandpass is not determined by the external feedback elements only but also by the real opamp gain characteristic (in particular, unwanted additional phase shift). This is, of course, no surprise because the frequency characteristics of the 741-type opamp is not sufficient for realizing a 40kHz bandpass.

As a consequence, a redesign of the circuit is necessary based on a clear specification (bandwidth, Q-Value, gain requirements?) and a high-frequency opamp with a sufficient slew rate. In this context, one should realize that an ac analysis is a small-signal analysis which cannot reveal slew rate problems.    

EDIT/UPDATE

Your calculation of R2 seems to be false. The correct value is R2=210 ohms. In this case and using the LM741 model, the center frequency is app. at 30 kHz. Using another opamp model (TL072) with a much larger gain-bandwidth product, the desired value of 40kHz can be achieved with good accuracy.