A discrete interfacial bonding element(BS1) was developed using the User Element (UEL) interface in ABAQUS, ant it’s capable of automatically calculating element length based on a specified bond-slip relationship. Numerical simulations of beam-cast-in-place column subassemblages were carried out using the simplified Eligehausen bond-slip constitutive model. The results were consistent with experimental data. The BS1 element demonstrated a 25% improvement in computational efficiency due to its automatic length calculating feature. A total of 30 finite element models of PC beam-column subassemblages were established for multi-parameter analysis, considering variables such as the thickness of the cast-in-place layer, the slab width of T-composite beams, and the direction of load application. Results show that subassemblages with 90° hooked bars and headed bars exhibit nearly equivalent bearing capacities, though the ductility index differs by 12% , favoring the 90° hooked bar. The slab width specified in Chinese design codes proved more reliable for bearing capacity performance than those in foreign codes. Furthermore, changes in load directions result in a 12% variation of bearing capacity and a 55% difference in deformation capacity, highlighting the need for further investigation. This study provides valuable insights and reliable references for the numerical analysis of prefabricated concrete structures considering bond-slip effects.