Solving a system of two coupled ordinary differential equations up to fourth-order

Question

I have the following (and terrible) system of ordinary differential equations up to fourth order (one of the equations is like a restriction for the other) for the one variable real functions $a(t)$ and $b(t)$:

$0=\frac{-12 \gamma (3 \alpha -\beta ) a(t)^2 b(t)^2 a''(t)^2-24 \gamma (3 \alpha -\beta ) a(t)^2 b(t) a'(t)^2 b''(t)+60 \gamma (3 \alpha -\beta ) a(t)^2 a'(t)^2 b'(t)^2-36 \gamma (3 \alpha -\beta ) b(t)^2 a'(t)^4+6 \gamma a(t)^2 b(t)^4 a'(t)^2+24 \gamma (3 \alpha -\beta ) a^{(3)}(t) a(t)^2 b(t)^2 a'(t)-24 b'(t) \left(\gamma (3 \alpha -\beta ) a(t) b(t) a'(t)^3+2 \gamma (3 \alpha -\beta ) a(t)^2 b(t) a'(t) a''(t)\right)+24 \gamma (3 \alpha -\beta ) a(t) b(t)^2 a'(t)^2 a''(t)-\Lambda a(t)^4 b(t)^2}{2 \gamma a(t)^4 b(t)^4}$

$0=-\frac{8 \gamma (3 \alpha -\beta ) a(t)^3 a^{(4)}(t) b(t)^3+12 \gamma (3 \alpha -\beta ) a(t)^2 b(t)^3 a''(t)^2-8 \gamma (3 \alpha -\beta ) a(t)^3 b^{(3)}(t) b(t)^2 a'(t)-120 \gamma (3 \alpha -\beta ) a(t)^3 a'(t) b'(t)^3+12 \gamma (3 \alpha -\beta ) b(t)^3 a'(t)^4+2 \gamma a(t)^2 b(t)^5 a'(t)^2+16 \gamma (3 \alpha -\beta ) a(t)^2 a^{(3)}(t) b(t)^3 a'(t)+60 b'(t)^2 \left(2 \gamma (3 \alpha -\beta ) a(t)^3 b(t) a''(t)+\gamma (3 \alpha -\beta ) a(t)^2 b(t) a'(t)^2\right)-16 b''(t) \left(2 \gamma (3 \alpha -\beta ) a(t)^3 b(t)^2 a''(t)-5 \gamma (3 \alpha -\beta ) a(t)^3 b(t) a'(t) b'(t)+\gamma (3 \alpha -\beta ) a(t)^2 b(t)^2 a'(t)^2\right)+4 a''(t) \left(\gamma a(t)^3 b(t)^5-12 \gamma (3 \alpha -\beta ) a(t) b(t)^3 a'(t)^2\right)-4 b'(t) \left(12 \gamma (3 \alpha -\beta ) a(t)^3 a^{(3)}(t) b(t)^2-12 \gamma (3 \alpha -\beta ) a(t) b(t)^2 a'(t)^3+\gamma a(t)^3 b(t)^4 a'(t)+18 \gamma (3 \alpha -\beta ) a(t)^2 b(t)^2 a'(t) a''(t)\right)-\Lambda b(t)^7}{2 \gamma a(t)^2 b(t)^7}$

Where I am using prime to denote derivative with respect to time $t$. In Mathematica software they look like this

Ec1 = 1/(
  2 \[Gamma] a[t]^4 b[t]^4) (6 \[Gamma]*a[t]^2 b[t]^4 a'[t]^2 - 
    36 (3 \[Alpha] - \[Beta]) \[Gamma]*b[t]^2 a'[t]^4 + 
    24 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t] b[t]^2 a'[t]^2 a''[t] - 
    12 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^2 b[t]^2 a''[t]^2 + 
    24 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^2 b[t]^2 a'[t] a'''[t] + 
    60 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^2 a'[t]^2 b'[t]^2 - 
    24 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^2 b[t] a'[t]^2 b''[t] - 
    a[t]^4 b[t]^2 \[CapitalLambda][t] - 
    24 ((3 \[Alpha] - \[Beta]) \[Gamma]*a[t] b[t] a'[t]^3 + 
       2 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^2 b[t] a'[t] a''[t]) b'[
      t])

Ec2 = Ecii = -(1/(
   2 \[Gamma] a[t]^2 b[t]^7)) (2 \[Gamma]*a[t]^2 b[t]^5 a'[t]^2 + 
    12 (3 \[Alpha] - \[Beta]) \[Gamma]*b[t]^3 a'[t]^4 + 
    12 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^2 b[t]^3 a''[t]^2 + 
    16 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^2*b[t]^3 a'[t]*a'''[t] + 
    8 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^3*b[t]^3 a''''[t] - 
    120 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^3 a'[t]*b'[t]^3 - 
    8 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^3*b[t]^2*a'[t]*b'''[t] + 
    b[t]^7 (-\[CapitalLambda]) + 
    60 ((3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^2 b[t]*a'[t]^2 + 
       2 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^3 b[t]*a''[t]) b'[
       t]^2 + 4 (\[Gamma]*a[t]^3 b[t]^5 - 
       12 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]*b[t]^3*a'[t]^2) a''[
      t] - 4 (\[Gamma]*a[t]^3 b[t]^4 a'[t] - 
       12 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]*b[t]^2*a'[t]^3 + 
       18 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^2*
        b[t]^2 a'[t] a''[t] + 
       12 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^3*b[t]^2 a'''[t]) b'[
      t] - 16 ((3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^2*b[t]^2*
        a'[t]^2 + 
       2 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^3*b[t]^2*a''[t] - 
       5 (3 \[Alpha] - \[Beta]) \[Gamma]*a[t]^3*b[t]*a'[t]*b'[t]) b''[
      t])

Here, $\alpha$,$\beta$,$\gamma$ and $\Lambda$ are constants of $O(1)$. Of course, a lot of initial conditions need to be given. Some reasonable initial conditions from the context of the problem are $a(t=1)=b(t=1)=1$, and their first derivatives should be positive $a'(t=1),b'(t=1)>0$ but not so big (0.5, for example). I don't have more information in order to give solid initial conditions for the higher derivatives, but it is reasonable again to choose them positive and the initial condition for the second derivative $b''(t=1)$ can be obtained from the first equation by consistency. One could choose, for example:

{a[1] == 1, a'[1] == 1, a''[1] == 1/2, a'''[1]==1/2 , b[1] == 1, 
 b'[1] == 1, b''[1] == -(103/24)}

The problem is that Mathematica is cannot solve this system, giving errors like "NDSolve::ndcf: Repeated convergence test failure at t == 1.`; unable to continue.".

I have tried to use AsymptoticDSolveValue too (https://reference.wolfram.com/language/ref/AsymptoticDSolveValue.html), but it doesn't work. I would like to know if it's possible to know something about the solution of this system, if one can check its singularities, behaviour near $t=0$, etc.