Multilayered metallic sheets fabricated by hot metal forming processes are attractive structural components in the aerospace and automobile industry. The strength of the mechanical bond between the different layers is critical to the applicability of such components in load bearing applications. While it is well understood that the process conditions such as temperature, deformation and rate of deformation influence the bond strength, a clear functional relationship between these conditions and the bond strength has yet to be established. The objective of the current study is to investigate the joining of Aluminum-6061 alloy specimens by hot compression bonding. Hot compression tests of beam specimen pairs were conducted at different temperatures (300-500℃), with different strains and strain rates (10-2-10-1 [1sec]). Coupled thermo-mechanical finite element models were utilized to compute the time dependent thermo-mechanical fields that develop along the specimen interface. The computed values were used to map the spatial distribution of a scalar parameter J which is hypothesized to govern bond formation. Mechanical tensile tests were used to determine the bond strength and expose the bonded surfaces. It is demonstrated that the computed J distribution can be correlated well with SEM observations of the debonded surface. It is also shown that the actual bonded area can be predicted from the computations for a certain range of critical J values.