Unsteady and two-dimensional numerical simulation is applied to study the transition process from steady convection to turbulence via subharmonic bifurcation in thermocapillary convection of a liquid bridge in the half-floating zone. The results of numerical tests show clearly the fractal structure of period-doubling bifurcations, and frequency-locking at f/4, f/8, f/16 with basic frequency f is observed with increasing temperature difference. The Feigenbaum universal constant is given by the present paper as delta(4) = 4.853, which can be compared with the theoretical value 4.6642016.