In this paper, a finite element model is formulated to simulate the crustal of deformation of faultbend folding with the large scale features using the large deformation theory in which the fault surface is represented with the master-slave method. The viscous-plastic material properties modeled by the Perzyna viscoplasticity theory are used to simulate the pressure solution creep and cataclasis respectively. Moreover, an interfacial viscous-plastic model is adopted to represent the nonlinear behavior of the fault surface. Effects on different structural models and parameters of fault-bend folds are simulated. Contours and development of the stress invariant, horizontal strain, equivalent viscous plastic strain by following an individual particle and traction on the fault surface are shown. The numerical results comparing with kinematic models are also discussed in order to infer a sequence of the deformation mechanisms.