The simulations of the five historic pianos are based on a numerical model previously developed by Juliette Chabassier (See References [13] and [14]). The action of the hammer is given by a force pulse, localized in time and space. This force pulse is derived from measurements of the string velocity, following an original method especially developed for this project (see Section 5. Applications and Reference [1]). The nonlinear model of strings accounts for the dependence of piano tones upon string amplitude. In this model, the transverse polarization of the string is coupled to its longitudinal motion. The soundboard model includes the ribs, bridges and apron. These additional components are viewed as localized zones of the soundboard with varying mass and stiffness. The thickness effects of these components (such as the torsional motion, for example) are ignored. For simulating the pressure field, the soundboard is the unique considered radiating element. The other parts of the instrument (rim, lid,...) are supposed to be passive. The meshes take the complete shape of the piano into account. Examples of simulated tones are found in the subpage D. Sound examples, showing the influence of some parameters (string tension, soundboard thickness,...). Comparison with measurements are also given.