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))
