四面体单元和六面体单元积分器选择问题

[MegrezLi]

请问一下四面体单元和六面体单元的积分器和组装器应该选择哪种，也是用LinearElasticityOperatorIntegrator，和BilinearForm.assembly吗？

[weihuayi]

你的问题还不够明确，你要求解的方程是哪一类方程，想用什么方法？

[MegrezLi]

求解三维固体力学问题，想用六面体单元或者四面体单元，目前卡在了HexahedronMesh中的shape_function
变量中的bc，输入的bc的结构在函数内调用爱因斯坦求和时就会报错。我的bc输入的是，bc = (((1 - xi) / 2, (1 + xi)/2),
((1 - eta)/2, (1 + eta)/2), ((1 - zeta)/2, (1 + zeta)/2))。

[MegrezLi]

return np.einsum('im, jn, ko->ijkmno', phi[0], phi[1], phi[2]).reshape(-1, ldof) 六面体中的shape_function，为什么形函数是六阶张量ijkmno？不应该是三阶吗？

[MegrezLi]

bc似乎应该输入的是拉格朗日多项式？