The pressure crashes observed in shot No. 71326 in the Experimental Advanced Superconducting Tokamak are numerically investigated with the three-dimensional, toroidal, and full resistive-magnetohydrodynamics code (CLT). Based on the experimental observations, the pressure crash is caused by the nonlinear evolution of the m/n = 2/1 double tearing mode (DTM), where n and m are the toroidal and poloidal mode numbers, respectively. However, we find that the m/n = 2/1 DTM is stable based on the safety factor (q) profile from the equilibrium code EFIT, which indicates that the original q profile is somewhat inconsistent with the actual profile due to q measurement uncertainty. Since there is no motional Stark effect diagnostic for this shot, the local information of the magnetic field is missing, which leads to the largest contribution to the discrepancy. If other information is perfectly known and the q profile is the only uncertainty, then we could provide some information for the EFIT reconstruction by comparing our simulation results with electron cyclotron emission signals to constrain the uncertainty of the q profile to a much smaller region. The influence of plasma rotation and the two-fluids effect is also discussed.