def bilinf(u, v, w): def C(T): trT = T[0, 0] + T[1, 1] return E / (1. + nu) * \ np.array([[T[0, 0] + nu / (1. - nu) * trT, T[0, 1]], [T[1, 0], T[1, 1] + nu / (1. - nu) * trT]]) return t**3 / 12.0 * ddot(C(dd(u)), dd(v))
def bilinf(u, v, w): from skfem.helpers import dd, ddot, trace, eye d = 0.1 E = 200e9 nu = 0.3 def C(T): return E / (1 + nu) * (T + nu / (1 - nu) * eye(trace(T), 2)) return d**3 / 12.0 * ddot(C(dd(u)), dd(v))
def bilinf(u, v, w): from skfem.helpers import dd, ddot return ddot(dd(u), dd(v))
def bilinf(u, v, w): return ddot(dd(u), dd(v))
def biharmonic(u, v, w): from skfem.helpers import ddot, dd return ddot(dd(u), dd(v))
def penalty_2(u, v, w): return ddot(-dd(v), prod(w.n, w.n)) * dot(grad(u), w.n)
def a_load(u, v, w): ''' for $a_{h}$ ''' return ddot(dd(u), dd(v))
def biharmonic(u, v, w): return ddot(dd(u), dd(v))