“师傅,上一讲您说有一个神仙公式,不用每次都切微元、列平衡方程,就能直接算出任意斜面上的应力向量。赶紧教教我吧!”
地点:深夜的力学研讨室。黑板上的斜截面推导图还未擦去,旁边堆着几本《弹性力学》教材。人物:
小生:师傅,上一讲您说有一个神仙公式,不用每次都切微元、列平衡方程,就能直接算出任意斜面上的应力向量。赶紧教教我吧!
老师傅:别急。在交给你这把终极武器之前,你得先彻底认清 Traction Stress(应力向量 )的局限性。 无论在单向拉伸还是纯扭转中,应力向量 都有一个致命缺陷——它太“短视”了。它只能忠实地描述你切开的那一个特定截面上的受力。但过物质内部的某一点,有着无穷多个方向的截面。如果我们想全面掌握这颗点被挤压、拉扯的“全局受力环境”,单靠一个面的 是远远不够的。我们需要一个能涵盖无穷多截面信息的终极物理量——这就是 Stress State(应力状态)。
老师傅:为了建立这个全局模型,数学家们引入了矩阵。但在写矩阵之前,咱们得先统一切应力的符号规矩。 材料力学教材里,常规定“使微元顺时针旋转的切应力为正”。但这在矩阵运算里会不方便。今天咱们要给出一个更常见的应力符号定义:
小生:这规矩好严谨!这样做有什么好处吗?
老师傅:最大的好处是,根据切应力互等定理,它能保证 。这样一来,我们可以用一个矩阵(二阶张量)来描述一点的应力状态了。在二维平面应力状态下,这个描述该点“受力基因”的应力张量 表达式如下,这是一个对称矩阵:
老师傅:好,现在见证奇迹的时刻到了。咱们切出一个最一般的二维斜面微元。左边面受 ,底面受 。 首先,我们要精确描述这个斜面的“指向”。引入单位法向量 ,它的分量由截面夹角 决定:
老师傅:斜面上的未知应力向量是 。我们老老实实列 向和 向的静力平衡方程:
小生:(拿起草稿纸化简)约掉 ,把 和 留在一边:
老师傅:敏锐!现在,把你算出的结果写成矩阵形式!
小生:
因为 ,所以中间那个矩阵是对称的。利用这个对称性,咱们可以写出最终最优雅的形式:
老师傅:这就是名震天下的柯西应力公式(Cauchy Stress Formula):
小生:师傅,这公式看着太简单了,它不仅能重算咱们之前的单向拉伸和纯剪切,简直像个“翻译机”!
老师傅:没错。咱们拿上一讲把你绕晕的“纯剪切状态”来验证一下。
在纯扭转(纯剪切)工况下,根据咱们新的张量符号规矩,右侧面(正 面)的切应力向下(负 向),所以 。顶面(正 面)的切应力向左(负 向),所以 。没有任何正应力。 所以,它的应力张量矩阵是:
现在,求法向量为 的斜面上的应力向量 。你用矩阵乘法算算看。
小生:(动笔飞快计算)
(笔尖一顿,猛地抬起头)天哪!!这和咱们上一讲通过画微元体、分解投影、列力学平衡方程辛辛苦苦算出来的应力向量分量( )完全一模一样!
第62夜 | 换个角度看世界:应力变换3-扭转工况下纯剪切的迷局
老师傅:哈哈!感受到数学的降维打击了吗?这个矩阵就像一个“应力生成器”。 它存储了该点处客观存在的、全方位的“应力基因”。当你给定一个截面的方向(输入法向量 ),这个矩阵就能瞬间完成线性变换,将其“变换”为该截面上的应力向量(输出 Traction Stress )。
老师傅:有了这个算子,你根本不需要再画微元体,也不需要去分析力的投影。哪怕在三维空间中,算子升级为 矩阵,魔法依然奏效:
小生:状态是矩阵本体,截面是输入向量,面力是输出结果。这套逻辑简直就像完美的计算机程序!
老师傅:没错。不过,如果换一个观察的坐标系x'-y',这个应力状态矩阵的每个分量会怎么变化呢?我们下一讲再讲。