Modified bond-based peridynamic approach for modeling the thermoviscoelastic response of bimaterials with viscoelastic–elastic interface

This study investigates the constitutive relationships for a modified bond-based peridynamic (MBB-PD) model to analyze bimaterial structures with a viscoelastic component subjected to mechanical and thermal loads. The hereditary integral viscoelastic constitutive relationship is used to model the viscoelastic behavior of a generalized Maxwell model, with viscoelastic constitutive coefficients obtained from the Prony series. To model the bimaterial interface between viscoelastic and elastic segments, various functions are proposed. Additionally, normal and shear bonds are defined as failure criteria based on critical stretches and angles. The fidelity of the proposed model is evaluated by comparing its thermoviscoelastic creep and recovery responses with finite element (FE) solutions for single and bimaterial plates. The sensitivity of the interface function to the horizon size is also examined. The results demonstrate that small horizon sizes still yield satisfactory results near the interface while reducing computational costs. Finally, the proposed MBB-PD model is applied to simulate a double cantilever beam (DCB) delamination test with a rubber interface, with predictions that are in good agreement with those from the extended finite element method (XFEM). The MBB-PD model proves to be a useful and efficient tool for analyzing bimaterial structures with viscoelastic components.