Combining trajectory surface hopping (TSH) method with constraint molecular dynamics, we have extended TSH method from full to flexible dimensional potential energy surfaces. Classical trajectories are carried out in Cartesian coordinates with constraints in internal coordinates, while nonadiabatic switching probabilities are calculated separately in free internal coordinates by Landau-Zener and Zhu-Nakamura formulas along the seam. Two-dimensional potential energy surfaces of groundS0and excitedS1states are constructed analytically in terms of torsion angle and one dihedral angle around the central ethylenic C=C bond, and the other internal coordinates are all fixed at configuration of the conical intersection. At this conical intersection, the branching ratio from the present simulation is 48 : 52 (33 : 67) initially starting from trans(cis)-Stilbene in comparison with experimental value 50 : 50. Quantum yield for trans-to-cis isomerization is estimated as 49% in very good agreement with experimental value of 55%, while quantum yield for cis-to-trans isomerization is estimated as 47% in comparison with experimental value of 35%.