Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
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) {}
};
see the Float page for the full documentation
see the point page for the full documentation
point operator()(const point &x) const
Float operator()(const point &x) const
Definition leveque.h:25
Float us
point operator()(const point &x) const
Definition leveque.h:26
u()
Definition zalesak.h:33