Nonzero initial conditions for a circuit solved using superposition theorem?
Since the system has initial conditions, the equations are differential equations. With your example, there are three 2nd order differential equations.

Since the equations are linear, the voltage sources (inputs) will always appear as the inhomogeneous terms (aka, input functions, forcing functions). When solved, each equation will produce a particular solution and a homogeneous solution with constant(s). With your example, 3 pairs of particular solutions and homogeneous solutions each with 2 integrating constants.

So it appears that there are two many constants to be determined by the initial conditions, but the homogeneous solutions will be linearly dependent, so the number of constants collapse back to just one set. In your example, the 6 constants collapse back to 2. Knowing this, you could just take the 3 particular solutions and any one of the homogeneous solutions, add them together and that is the general solution for the system. Then fix the 2 constants with the initial conditions.