c1 = {k1 ->  1.0289716,  
          k2 ->  1.23924503, s1 ->  0.27885572, 
          s2 ->  0.29125954, Z ->  3};
  
c2 = {k1 ->  1.02480344,  k2 ->  1.232443643, 
          s1 ->  0.27877144, s2 ->  0.28887416, 
          Z ->  3};

Vne = -3*Z/P - 2*Z/(P + Q) - 
          2/(P + Q + R) - (3 - (R/Q)^2)/
            Q;
Vee = 4/(P + Q) + 3*s1/P + 3*s2/Q;
Vnn = Z/(P + Q + R);
T = 2.25*k1/P^2 + 2.25*k2/Q^2;
func = T + Vne + Vee + Vnn  /.  c1
N[t = FindMinimum[func,  {P, 0.6}, {Q, 2.0}, {R, 0.5}], 
      10]

u1 = t[[2]];
N[-(Vne + Vee + Vnn)/T  /. u1  /. c1, 
      10]
N[(P + Q + R)*0.529177  /. u1  /. c1, 10]

func = 
      T + Vne + Vee + Vnn  /.  c2
N[t = FindMinimum[func,  {P, 0.5}, {Q, 2.0}, {R, 0.5}], 
      10]
u2 = t[[2]];
N[-(Vne + Vee + Vnn)/T  /. u2  /. c2, 
      10]
N[(P + Q + R)*0.529177  /. u2 /. c2, 
      10]
p = P  /. u2;
q = Q  /.  u2;
r = R /.  u2;

plot1 = 
      Graphics[{
Circle[{0, 0}, q], Circle[{0, q + p}, p], 
          Disk[{0, p + q}, 
            0.05], 
Line[{{0, -r}, {0, 
                p + q}}]}]
Show[
      plot1, {AspectRatio ->  
          Automatic}]