scikit-FEM
Posted wangshixi12
tags:
篇首语：本文由小常识网(cha138.com)小编为大家整理，主要介绍了scikit-FEM相关的知识，希望对你有一定的参考价值。
from skfem import *
m = MeshTri()
m.refine(4)
e = ElementTriP1()
basis = InteriorBasis(m, e)
@bilinear_form
def laplace(u, du, v, dv, w):
return du[0]*dv[0] + du[1]*dv[1]
@linear_form
def load(v, dv, w):
return 1.0*v
A = asm(laplace, basis)
b = asm(load, basis)
I = m.interior_nodes()
x = 0*b
x[I] = solve(*condense(A, b, I=I))
m.plot3(x)
m.show()
Please see the directory examples for more instructions.
Contributors
以上是关于scikit-FEM的主要内容，如果未能解决你的问题，请参考以下文章
scikit-FEM-例2-用Morley元在方形区域上解板弯曲问题