贡献者: ACertainUser
1
图 1:斜切面 $ \boldsymbol{\mathbf{n}} $ 上的应力 $ \boldsymbol{\mathbf{\sigma}} _\nu$。$ \boldsymbol{\mathbf{n}} $ 是面的单位法向量
在应力中,我们已经处理了微元体上的应力问题。那么新的问题来了,如图 1 所示,假设我们已经在某个坐标系中算出了微元体的各个应力,如果我们在微元体中取一个切面,那么要如何计算这个切面上的应力?
跨界大牛柯西给出了一个漂亮的回答,因此该公式也称柯西应力公式。柯西应力公式写成矩阵乘法的形式:
\begin{equation}
\begin{pmatrix}
\sigma_{\nu,1}\\
\sigma_{\nu,2}\\
\sigma_{\nu,3}\\
\end{pmatrix}
=
\begin{pmatrix}
\sigma_{11} & \tau_{12} & \tau_{13} \\
\tau_{21} & \sigma_{22} & \tau_{23} \\
\tau_{31} & \tau_{32} & \sigma_{33} \\
\end{pmatrix}
\begin{pmatrix}
n_1\\
n_2\\
n_3\\
\end{pmatrix}~.
\end{equation}
还可以写成更紧凑的张量乘法形式,
\begin{equation}
\boldsymbol{\mathbf{\sigma}} _\nu = \sigma \cdot n ~,
\end{equation}
或指标形式
\begin{equation}
\boldsymbol{\mathbf{(}} \boldsymbol{\mathbf{\sigma}} _\nu)_i = \sum_j \sigma_{ij} n_j~.
\end{equation}
因为已经被投影,此处的 $ \boldsymbol{\mathbf{\sigma}} _\nu$ 已经是一个矢量,而不是原先的张量。此外,斜面上的应力 $ \boldsymbol{\mathbf{\sigma}} _\nu$ 不一定垂直于该平面,而可以朝任意方向。该公式的论证方法大致是在三角形台中运用力的平衡方程,此处按下不表。
图 2:$\sigma_\nu$ 中各个分量的含义
回过头来,我们思考一下,我们依据上文方法写出的 $\sigma_\nu=
\begin{pmatrix}
\sigma_{\nu,1}\\
\sigma_{\nu,2}\\
\sigma_{\nu,3}\\
\end{pmatrix}$ 中各个分量的含义是什么?显然(但又有点让人困惑),各个分量仍基于原先的坐标系,如图 2 所示。
图 3:$\sigma_\nu$ 的斜面表示
我们写出斜面所受应力,总希望他和斜面有点关系(而不是原先的坐标系);因此,我们得想办法把它改写为和斜面有关系的、更 “直观” 的形式,例如表达为沿斜面法向与切向的应力。
只需再做一步简单的向量投影就能实现这个目的:法向应力 $ \boldsymbol{\mathbf{\sigma}} _n$ 即是 $ \boldsymbol{\mathbf{\sigma}} _\nu$ 在 $ \boldsymbol{\mathbf{n}} $ 上的投影:
\begin{equation}
\boldsymbol{\mathbf{\sigma}} _n = ( \boldsymbol{\mathbf{\sigma}} _\nu \cdot \boldsymbol{\mathbf{n}} ) \hat{ \boldsymbol{\mathbf{n}} }~.
\end{equation}
他的数值大小写成张量乘法的形式
\begin{equation}
\sigma_n = n \cdot \sigma \cdot n
\qquad \sigma_n =\sum_i \sum_j \sigma_{ij} n_i n_j~,
\end{equation}
而切向应力就是
\begin{equation}
\boldsymbol{\mathbf{\tau}} = \boldsymbol{\mathbf{\sigma}} _\nu - \boldsymbol{\mathbf{\sigma}} _n~.
\end{equation}
至此我们已经解决了文章开头的问题,剩下的只是繁琐的数学计算。不过,我们还可以额外思考另一个问题:如式 1 所示,柯西斜面应力公式相当于一个矩阵乘向量的矩阵乘法。
$$\sigma_\nu = \sigma \cdot n~.$$
如果你的数理直觉足够敏锐,就会意识到这里存在一个矩阵的本征值与本征向量问题:存不存在一种 $n_0$,使
$$
\sigma \cdot n_0 = \lambda n_0~.
$$?
更进一步的思考后,你会发现如果 $n_0$ 存在,那么 $ \boldsymbol{\mathbf{\sigma}} _\nu = \lambda \boldsymbol{\mathbf{n}} _0$ 将平行于 $ \boldsymbol{\mathbf{n}} _0$,也即在这个面上将只有正应力。那么,$n_0$ 将是一个很特殊的平面。到这一步,你其实就已经想明白了主应力问题。
1. ^ 本文参考了陆明万的《弹性力学》。
致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者
热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。