DSolve 2nd Order Coupled Partial Differential Equations

I am trying to use Mathematica to solve 2 coupled differential equations.

My equations are of the form $$\ddot{x}_i + A_{il} \partial^l A^{jk} ( \dot{x}_b \dot{x}_c - y_b y_c ) =0$$

where $x_a=x_a(t)$ are functions of $t$ and $y_a$ are constant vectors (I'm happy to set these to zero to begin with if it will help get a solution!)

For the problem at hand, $a=1,2$ and so I have two coupled PDEs in $x_1$ and $x_2$.

The $2 \times 2$ matrix $A^{jk}$ is given by

$$A^{jk} = \frac{1}{(x_1^3-x_1 x_2^2)^2} \begin{pmatrix} \frac{x_1^4+x_2^4}{4} & x_1^3 x_2 \\ x_1^3 x_2 & \frac{x_1^4-2x_1^2x_2^2}{2} \end{pmatrix}$$

and lastly, the matrix $A_{ad}$ is just the inverse of the above matrix.

I have attempted to solve these using DSolve and have produced the following code (with $y_a=0$):

A = 1/((x1^3 - x1*x2^2)^2)*{{(x1^4 + x2^4)/4, x1^3*x2}, {x1^3*
x2, (x1^4 - 2*x1^2*x2^2)/2}}

xdot = {xdot4, xdot5}

y = {y4, y5}

S = {xdotdot1, xdotdot1} + 1/2*Inverse[A].{D[#, x1], D[#, x2]} &[
xdot.A.xdot - y.A.y];

St = S /. {xdotdot1 -> x1''[t], xdotdot2 -> x2''[t], xdot1 -> x1'[t],
xdot2 -> x2'[t], x1 -> x1[t], qx2 -> qx2[t], y1 -> 0, y2 -> 0};

system = {St[[1]] == 0, St[[2]] == 0};

DSolve[system, {x1, x2}, t]


I am fairly new to Mathematica and produced this code with a lot of help from a friend. The code works fine on simple examples where the PDEs are not coupled but in this case I just get the output

DSolve[.....]


i.e. it is not able to solve them. I understand that coupled PDEs are much harder for it to solve but given that it's only a 2x2 matrix, I hope it should be doable!

I have thought about trying to use NDSolve to get a numerical solution since apparently Mathematica is good at this but the problem is I'm not sure about the initial conditions - if there's no way to get an exact solution then perhaps we can give this approach a bit more thought.

I hope I've explained the problem well. Any advice/solutions greatly appreciated.

Thanks.