We many times analysed large welded structures or automotive panel spot welds. One of the methods we used was include beam elements with appropriate properties then extract forces. you can have coincident nodes connecting beams also on contact surfaces.