yield_slip_circle.h

The yield slip problem on a circle – exact solution.

The yield slip problem on a circle – exact solution

struct u {
  Float operator() (const point& x) const { return (1-norm2(x))/4 + us; }
  u (Float S, Float n, Float Cf) : us(pow(max(Float(0),(0.5-S)/Cf), 1/n)) {}
  protected: Float us;
};
struct grad_u {
  point operator() (const point& x) const { return -x/2; }
  grad_u (Float S, Float n, Float Cf) {}
};
struct lambda {
  Float operator() (const point& x) const { return 1./2; }
  lambda (Float S, Float n, Float Cf) {}
};