The magnetization of type-II superconductors has been usually described using the critical state models (CSMs). However, the CSMs does not consider the time-dependent magnetic relaxation. In this work, to study the influence of magnetic relaxation on the magnetic hysteresis loops (MHLs), critical current density Jc, and flux pinning force F
p, a numerical method proposed by Qin et al. [10.1103/PhysRevB.54.7536] was adopted and extended to the quasi-two-dimensional case. MHLs at different temperatures measured by using applied field with different sweeping rates can be well reproduced. The critical current density after relaxation Js, and Fp can also fit the experimental results well. Based on the numerical simulation, the non-scaling behavior of field-dependent normalized pinning force has been proved to be from the magnetic relaxation. By comparing the normalized pinning force with D. Dew-Hughes model, the dominant pinning type in Fe1+y
Te0.6Se0.4 has been confirmed to be the volume ∆κ pinning. In addition, the second peak effect is found to be related to the flux lattice transition from elastic lattice to plastic lattice (E-P transition).