We present the first tests and results from a study of QCD with two fl
avours of dynamical Wilson fermions using the Hybrid Monte Carlo Algor
ithm (HMCA) on APE100 machines. The simulations have been performed on
6(4) lattice for the pure gauge HMCA and on 8(4), 12(3) X 32 lattices
for full QCD configurations. We discuss the inversion algorithm for t
he fermionic operator, the methods used to overcome the problems arisi
ng using a 32 bit machine and the implementation of a new random numbe
r generator for APE100 machines. We propose different scenarios for th
e simulation of physical observables, with respect to the memory capac
ity and speed of different APE100 configurations.