<p>它确实被替换了,但是<code>subs</code>在矩阵上不可变地工作。<code>applyfunc</code>也不是可变的,很不幸。我明白了</p>
<pre><code>In [10]: pprint(stiff.subs({-nuxy-2*nuxz**2+1:m}))
⎡ ⎛ 2 ⎞ ⎛ 2⎞ ⎤
⎢ Eₓ⋅⎝nu_xz - 1⎠ -Eₓ⋅⎝-nu_xy - nu_xz ⎠ Eₓ⋅nu_xz ⎥
⎢ ────────────────────────────────────── ──────────────────────────────────────── ──────── ⎥
⎢ 2 2 2 2 2 2 m ⎥
⎢ nu_xy + 2⋅nu_xy⋅nu_xz + 2⋅nu_xz - 1 - nu_xy - 2⋅nu_xy⋅nu_xz - 2⋅nu_xz + 1 ⎥
⎢ ⎥
⎢ ⎛ 2⎞ ⎛ 2 ⎞ ⎥
⎢ -Eₓ⋅⎝-nu_xy - nu_xz ⎠ Eₓ⋅⎝nu_xz - 1⎠ Eₓ⋅nu_xz ⎥
⎢──────────────────────────────────────── ────────────────────────────────────── ──────── ⎥
⎢ 2 2 2 2 2 2 m ⎥
⎢- nu_xy - 2⋅nu_xy⋅nu_xz - 2⋅nu_xz + 1 nu_xy + 2⋅nu_xy⋅nu_xz + 2⋅nu_xz - 1 ⎥
⎢ ⎥
⎢ E_z⋅nu_xz E_z⋅nu_xz -E_z⋅(nu_xy - 1)⎥
⎢ ───────── ───────── ────────────────⎥
⎣ m m m ⎦
</code></pre>
<p>有计划使Matrix在默认情况下是不可变的,然后使MutableMatrix在所有操作中完全正常工作。见<a href="https://code.google.com/p/sympy/issues/detail?id=3410" rel="nofollow">https://code.google.com/p/sympy/issues/detail?id=3410</a>。但这还没有发生。</p>