A new type of steel moment resisting frame with bottom ? ange friction devices (BFFDs) has been developed to provide self-centering capacity and energy dissipation, and to reduce permanent deformations under earthquakes. This paper presents a numerical simulation of self-centering beam-column connections with BFFDs, in which the gap opening /closing at the beam-column interfaces is simulated by using pairs of zero-length elements with compression-only material properties, and the energy dissipation due to friction is simulated by using truss elements with specified hysteretic behavior. In particular, the effect of the friction bolt bearing against the slotted plate in the BFFDs was modeled, so that the increase in lateral force and the loss of friction force due to the bolt bearing could be taken into account. Parallel elastic-perfectly plastic gap (ElasticPPGap) materials in the Open System for Earthquake Engineering Simulation (OpenSees) were used with predef ined gaps to specify the sequence that each bolt went into the bearing and the corresponding increase in bending stiffness. The MinMax material in OpenSees is used to specify the minimum and maximum values of strains of the ElasticPPGap materials. To consider the loss of friction force due to bolt bearing, a number of parallel hysteretic materials were used, and the failure of these materials in sequence simulated the gradual loss of friction force. Analysis results obtained by using the proposed numerical model are discussed and compared with the test results under cyclic loadings and the seismic loading, respectively. This article presents a numerical simulation of self-centering beam- column connections with BFFDs, where the gap opening / closing at the beam-column interfaces is simulated by using pairs of zero-length elements with compression-only material properties, and the energy dissipation due to friction is simulated by using truss elements with specified In particular, the effect of the friction bolt bearing against the slotted plate in the BFFDs was modeled, so that the increase in lateral force and the loss of friction force due to the bolt bearing could be taken into account. Parallel elastic- perfectly plastic gap (Elastic PPGap) materials in the Open System for Earthquake Engineering Simulation (OpenSees) were used with predefined gaps to specify the sequence that each bolt went into the bearing and the corresponding increase in bending stiffness. The MinMax material in OpenSees is used to specify the minimum and maximum values ​​of strains of the Elastic PPGap materials. To consider the loss of friction force due to bolt bearing, a number of parallel hysteretic materials were used and the failure of these materials in used in the gradual loss of friction force. .