为了开始检查非商业有限元程序的结果，我试图从程序计算的应力和应变状态，例如，计算给定方向上的杨氏模量E。将E视为一些角坐标 θ 的函数，它可以简单地计算为
E (θ) = 1 / ( v (θ)⊗<strong>v(θ): S : v (θ)⊗<strong>v(θ))
在张量符号中，其中v (θ) 是根据 θ 定义指向的酉向量。我在这里用张量表示法写，但这些张量显然在数字上被处理为矩阵和向量。
问题是我广泛尝试从应力和应变张量（张量 ← σ，ε ≡ {σ}，{ε} → 向量）计算 [S] 但没有成功。
有没有办法做到这一点？我找到了计算 Lamé 参数的各向同性介质的方法，但不适用于一般的各向异性介质。
请注意，由于处理的介质主要是各向异性的，我不是在寻找弹性参数（Lamé 参数）来（重新）构造 [S]，而是从给定的应力-应变状态计算 [S] 并得出在模拟开始时预先构造相同的[S]。
此外，数值实验不应被限制为单轴左右。我想知道是否有一个通用的程序。