When we have 2 massive bodies coming close together say 2 black holes or 2 massive stars, how do their respective metrics/spacetime curvature combine in the space in between them?
Do we write
$$G_{\mu\nu}^1 + G_{\alpha\beta}^2 = \kappa (T_{\mu\nu}^1 + T_{\alpha\beta}^2)$$
With $G_{\mu\nu}^1$ being star 1 and $G_{\alpha\beta}^2$ being star 2 and same for the stress-energy tensors.
Or do we write it as something else?
I have not seen an example of this so please excuse if it sounds obvious.
A core feature of GR is that it is a nonlinear theory, meaning you cannot naively combine two solutions to get a new one. This is one of the aspects that make it particularly difficult to handle.
To give a first example, consider a Schwarzschild black hole with some mass. We can try to consider the scenario with two black holes now. However, the resulting gravitational field is not the sum of the two previous ones: in the Schwarzschild solution for a single black hole, the spacetime is time-translation invariant, but in this new setup the two black holes are going to move toward each other and eventually merge. Curiously, in both of these solutions (one black hole or two black holes) the stress tensor is the same ($T_{\mu\nu} = 0$ everywhere).
Even in what regards the Einstein tensor your expression may not be correct. Although it initially may seem like simply summing the Einstein equation for each case, notice the right-hand side of the Einstein equations has an implicit dependence on the actual metric of spacetime. When you bring two sorts of matter together, the metric changes, and hence the stress tensor itself ends up changing as well.
In short, there is no clear cut rule to combining solutions in GR. You have to start from scratch and solve the equations again whenever you change the matter content.
As an addition to the other great answers: The equation you wrote doesn't make sense, because you can't add tensors with different indices and expect the result to be Lorentz-invariant. That's like adding the $x$-component of one vector to the $y$-component of another vector. Sure, you can do that but the result you get will depend on the coordinate system you chose.
The Einstein Field Equations are not linear, in the sense that if you have two metric tensors that are solutions for a given stress-energy tensor, their linear combination is not (generally) also a solution. This is contrast to e.g. the Schrödinger equation, which is linear.
So beyond a very limited set of known analytic solutions you will have to resort to numerical methods. For the specific case you have mentioned (the 2-body problem in General Relativity) there are a few known solutions for special cases. See e.g. General Relativity 2-Body Closed Form
Only in the weak limit(e.g Earth+Moon system),generally no because EFE are not linear.
