A multidomain SGBEM approach for shakedown analysis as a conic quadratic optimization problem