In this paper, we present a thorough treatment of a harpsichord plectrum-string interaction which allows for large end deflections and both transverse motions of the string. We start from the general equations of motion of a bent beam, and an accurate shape of the plectrum is calculated, agreeing with existing known cantilever beam models when end deflections are assumed small. All the governing forces on the string are considered, and the complete motion of the string up to its release is simulated, allowing for future implementation on physical model sound synthesis of strings. Simulation results agree with what is experienced playing a real harpsichord string.