dealii中unit_symmetric_tensor的用法

在deal.II中，unit_symmetric_tensor是一个函数，用于返回一个单位对称张量。它的用法如下：

#include <deal.II/base/symmetric_tensor.h>
using namespace dealii;

SymmetricTensor<2,dim> C = unit_symmetric_tensor<dim>();

其中，dim是维度，如dim=2表示二维空间，dim=3表示三维空间。unit_symmetric_tensor()返回一个dim x dim的单位对称张量C，即C[i][j] = 1.0 (i == j)，C[i][j] = 0.0 (i != j)。这个张量可以用于定义材料的应力/应变关系，比如Hooke定律中的弹性常数矩阵。