I am trying to obtain the state space model of the DC microgrid shown in the figure. It is made of three converters connected in parallel. Each one is modeled with a voltage source that outputs the nominal voltage \$V_{ref}\$, a CCVS that models a droop control, and an RLC network that models the dynamics of the output voltage caused by load changes. The differential equations of one converter (for example, nº 1) are the following:
$$C_1\frac{dv_{c1}}{dt} + I_{o1} = i_{l1} +\frac{(V_{ref} - K_1I_{o1} - v_{c1})}{R_1}$$ $$L_1\frac{di_{l1}}{dt} + v_{c} = V_{ref} - K_1I_{o1} $$ From there, we can write the state space representation of one converter as a subsystem:
$$\begin{pmatrix} \frac{dv_{c1}}{dt} \\ \frac{di_{l1}}{dt} \end{pmatrix} = \begin{pmatrix} -\frac{1}{C_1R_1} & \frac{1}{C_1} \\ -\frac{1}{L_1} & 0 \end{pmatrix} \begin{pmatrix} v_{c1}\\ i_{l1} \end{pmatrix} + \begin{pmatrix} \frac{1}{R_1C_1} & -\frac{K_1}{R_1C_1} - \frac{1}{C_1} \\ \frac{1}{L_1} & -\frac{K_1}{L_1} \end{pmatrix} \begin{pmatrix} V_{ref} \\ I_{o1} \end{pmatrix} $$
Now the issue is that we don't know the output currents of each converter \$I_{oi}\$, but we know the total load \$I_L\$, and we want to write the state space representation of the entire system, where the input variables are just \$I_L\$ and \$ V_{ref}\$. Is there any way to do that by combining the state space models of every converter? The connection in parallel implies the following restrictions:
$$I_{o1} + I_{o2} + I_{o3} = I_L$$ $$v_{c1} = v_{c2} = v_{c3}$$
The total load is equal to the sum of the output currents, and the capacitor voltages are the same for every converter.