默认参数 (λ=45°, φ=45°, σ=200 MPa, τcrss=50 MPa) 时，m = 0.500，分解剪应力 τr = 100 MPa，降伏应力 σy = 100 MPa，安全系数 SF = σy/σ = 0.50。SF < 1 表示已发生滑移。当 λ 或 φ 接近 0°/90° 时，m 变小，σy 快速上升。
蓝箭头：拉伸轴 σ（竖直方向） / 半透明蓝色平面：滑移面 / 绿箭头：滑移方向 / 红箭头：滑移面法线 / 黄色弧：λ（拉伸轴与滑移方向的角度） / 紫色弧：φ（拉伸轴与法线的角度）。当滑移面倾角为 45° 时，λ=φ=45°，施密特因子达到最大值 0.5。反之当面垂直或平行于拉伸轴时，m=0。
λ ) 曲线
横轴：λ (deg, 0~90) / 纵轴：施密特因子 m / 蓝实线：φ=45° 固定时的 m(λ) = cos(λ) × cos(45°) 余弦曲线 / 黄色标记：当前 λ 的工作点 / 绿色虚线：理论最大值 m_max = 0.5。当 φ=45° 时，λ=0° 时 m=0.7071，λ=45° 时 m=0.5，λ=90° 时 m=0。
施密特因子（Schmid factor）：
$$m = \cos\lambda \, \cos\phi$$
分解剪应力与 Schmid 定律滑移启动条件：
$$\tau_r = \sigma\,m, \qquad \tau_r \geq \tau_{\mathrm{crss}}\;\Rightarrow\;\text{滑移启动}$$
单晶体降伏应力与安全系数：
$$\sigma_y = \frac{\tau_{\mathrm{crss}}}{m}, \qquad \mathrm{SF} = \frac{\sigma_y}{\sigma}$$
其中 $\lambda$ 是拉伸轴与滑移方向的夹角，$\phi$ 是拉伸轴与滑移面法线的夹角，$\sigma$ 是施加应力（MPa），$\tau_{\mathrm{crss}}$ 是临界分解剪应力（MPa）。物理约束条件为 $0 \leq m \leq 0.5$，最大值 $m=0.5$ 在 $\lambda=\phi=45°$ 时达到。
施密特因子模拟器简介
🙋
我是第一次听说"施密特因子"。这是用来做什么的？
🎓
这是一个几何系数，用来预测单晶金属沿不同方向拉伸时的降伏难度。德国物理学家 Erich Schmid 在 1924 年提出了这个概念。公式是 m = cosλ cosφ，其中 λ 是拉伸轴与滑移方向的夹角，φ 是拉伸轴与滑移面法线的夹角。本工具默认参数 (λ=45°, φ=45°, σ=200 MPa, τcrss=50 MPa) 时，m=0.500，分解剪应力 τr=100 MPa，降伏应力 σy=100 MPa，安全系数 SF=0.50。SF=0.5 表示"施加的应力已经是降伏应力的 2 倍"，即已发生滑移的状态。
🙋
为什么施密特因子的最大值是 0.5 而不是 1.0？
🎓
这是一个很好的问题。你应该知道，根据 Mohr 圆理论，与拉伸轴垂直的面上的"最大剪应力"是 σ/2。施密特因子遵循同样的原理：当滑移面相对轴倾斜 45° 时（λ=φ=45°），cos(45°)×cos(45°) = 0.7071×0.7071 = 0.500 达到最大。物理上，当面垂直或平行于拉伸方向时，滑移方向或法线之一会与轴直交，cos=0，所以 m=0 时不会发生滑移。在本工具右边的 m(λ) 曲线中，当你把 λ 从 0° 改到 90° 时，会看到一条以 0.5 为顶点的余弦型曲线。
🙋
单晶体示意图中的蓝箭头、绿箭头、红箭头和蓝色平面分别代表什么？
🎓
蓝色竖箭头是拉伸轴 σ，半透明蓝色平面是晶体内的滑移面，绿箭头是滑移面上的滑移方向，红箭头是滑移面的法线 n。λ 是蓝箭头（拉伸轴）与绿箭头（滑移方向）的夹角，φ 是蓝箭头与红箭头（法线）的夹角。当你用滑块把 φ 改到 0° 时，平面变得水平，法线与拉伸轴平行；改到 90° 时，平面变得竖直，法线变成水平。实际 FCC 金属（铜、铝）有 12 个滑移系（3 个滑移面 × 4 种方向），每个都有不同的 m 值。在多晶中，具有最大 m 值的滑移系会首先激活，这就是 Schmid 定律。
🙋
τcrss 是什么？当我拖动滑块改变 τcrss 时，σy 会怎样变化？
🎓
τcrss 是临界分解剪应力（Critical Resolved Shear Stress），是滑移系开始滑移所需的最小剪应力。这是材料本身的一个特性，取决于材料种类、温度和位错密度。纯铜大约是 0.5 MPa，加工硬化后升到 50 MPa，BCC 钢是 30~100 MPa，HCP 镁的底面滑移大约 0.5~2 MPa。Schmid 定律告诉我们 σ m = τcrss 时开始滑移，所以 σy = τcrss / m。本工具默认值中 τcrss=50 MPa, m=0.5，所以 σy=100 MPa。如果你把 τcrss 增加一倍到 100 MPa，σy 也会增加一倍到 200 MPa。这个线性关系反映了"加工硬化""析出强化""晶粒细化"等强化机制的本质——都是在提高 τcrss。
🙋
右边 m(λ) 曲线显示"黄色标记在 λ=45°，m=0.5"，安全系数 SF=0.5 应该怎样理解？
🎓
SF = σy / σ，是降伏应力相对施加应力的比值。SF > 1 表示"还没有降伏"，SF = 1 表示"刚好降伏"，SF < 1 表示"已经发生塑性变形"。本工具默认值中 σy=100 MPa，σ=200 MPa，所以 SF=0.50，意思是"施加的应力已经是降伏应力的 2 倍"，也就是已经大幅塑性变形了。如果你把 σ 改成 100 MPa，SF 就变成 1.00（刚好降伏点）；改成 50 MPa 时 SF=2.00（安全裕度是 2 倍）。或者你可以把 λ 改成 60°，那么 m 会变成 0.354（当 φ=45° 时）→ σy=141 MPa → SF=0.71。通过调节角度就能改变强度，这展现了结晶塑性的本质。
常见问题
什么是施密特因子？
施密特因子 m 是单晶体拉伸试验中，轴应力 σ 在滑移面上分解为剪应力的比例系数，定义为 m = cosλ cosφ。其中 λ 是拉伸轴与滑移方向的夹角，φ 是拉伸轴与滑移面法线的夹角。物理上 m 的范围为 0 到 0.5，最大值 0.5 在 λ = φ = 45° 时达到。本工具默认参数 (λ=45°, φ=45°, σ=200 MPa, τcrss=50 MPa) 下，m=0.500，分解剪应力 τr=100 MPa，降伏应力 σy=100 MPa，安全系数=0.50，表示刚好开始滑移的条件。
为什么施密特因子的最大值是 0.5？
m = cosλ cosφ 的最大值在 λ = φ = 45° 时达到，cos(45°) ≈ 0.7071，所以 m_max = 0.7071 × 0.7071 = 0.500。这与 Mohr 圆理论一致，即当滑移面相对于拉伸轴倾斜 45° 时会产生最大剪应力。实际的多晶金属中，来自最活跃滑移系（施密特因子接近 0.5）的晶粒会优先开始塑性变形。
什么是临界分解剪应力 τcrss？
临界分解剪应力（Critical Resolved Shear Stress, τcrss）是滑移系开始滑移所需的最小剪应力，取决于材料、温度和位错密度。FCC 纯金属典型值为 0.5 MPa（纯铜）到 50 MPa（加工硬化后），BCC 铁为 30~100 MPa，HCP 镁的基面滑移为 0.5~2 MPa。根据 Schmid 定律，当分解剪应力 τr = σ m 达到 τcrss 时，该滑移系开始滑移。本工具允许在 5~300 MPa 范围内调节 τcrss 来观察对降伏应力的影响。
单晶体的降伏应力 σy 如何计算？
根据 Schmid 定律 σ m = τcrss，单晶体的降伏应力为 σy = τcrss / m。这表明施密特因子越大（接近 45° 取向），单晶体在越低的应力下就会发生降伏。本工具默认参数 (λ=45°, φ=45°, τcrss=50 MPa) 给出 σy = 50 / 0.5 = 100 MPa。反之，当 λ 或 φ 接近 0° 或 90° 时，m 变小，σy 会趋向无穷大，这称为"困难取向（hard orientation）"，体现了晶体塑性变形对取向的强烈依赖性。
现实应用
镍基单晶超合金涡轮叶片： 喷气发动机高压涡轮叶片（如 CMSX-4 和 René N5）使用单晶叶片，其晶体生长方向（通常是 [001]）被精确控制为与主应力方向一致。FCC 结构 {111}<110> 滑移系在 [001] 拉伸方向上的最大施密特因子约为 0.408，比多晶平均值（m≈0.5）高约 20%，从而获得更高的降伏应力。通过本工具改变 λ 和 φ 的值，你可以直观地看到这种各向异性强化的效果，理解为什么实际叶片要如此精确地控制晶体取向。
硅单晶圆片的位错控制： 半导体制造中的 Czochralski 法和浮区法生长的硅单晶圆片，在引上过程中受到热应力，容易在 {111}<110> 滑移系上导入位错，损害器件性能。通过模拟圆片主应力方向 [100] 相对各滑移系的施密特因子分布，可以找出最容易活化的滑移系，进而优化引上条件（温度梯度、引上速度等）来最小化位错。本工具是这种分析的教学入门，让你理解取向与 m 的对应关系。
镁合金板材的各向异性预测： 轻量化材料 Mg-Al-Zn 系（AZ31、AZ91）的 HCP 结构表现出极端的各向异性：底面 (0001) 滑移的 τcrss 只有 0.5~2 MPa，而非底面滑移（锥面、柱面）的 τcrss 则达 30~80 MPa，相差几十倍。压延板的集合组织（织构）与拉伸方向的夹角决定了有多少晶粒的滑移系被激活。通过本工具的滑块，你可以体验当改变取向时降伏应力如何剧烈变化（这在铜等 FCC 金属中没有这么极端），从而理解 HCP 材料的特殊强化机理。
结晶塑性有限元法（CPFEM）的教学工具： 现代高级 CAE 如 ABAQUS UMAT 和 DAMASK 中的 CPFEM 在每个积分点都设置多个滑移系，逐个计算它们的施密特因子、分解剪应力和位错密度演化。本工具正是这种计算的最基本单元——一个滑移系的施密特因子和分解剪应力——以实时、可视化的形式呈现，可作为讲座和自学教材，帮助学生建立对 CPFEM 物理基础的直观理解。
常见误解与注意事项
最常见的误解是"施密特因子可以直接应用到多晶材料" 。实际上，Schmid 定律本来是单晶的单个滑移系的规则。在多晶中，每个晶粒方位不同，相邻晶粒还有变形协调的约束，所以需要 Taylor 的多晶模型（M=3.06 是典型值）或 von Mises 的 5 滑移系准则。多晶的降伏应力大约是 M × τcrss，相当于在单晶结果上乘以一个因子。本工具展示的是单晶模型，实际多晶分析需要更复杂的结晶塑性有限元。
另一个常见误解是"λ + φ = 90° 总是成立" 。实际上 λ 和 φ 是两个独立的量：λ 是滑移方向与拉伸轴的角度，φ 是滑移面法线与拉伸轴的角度。由于滑移方向在滑移面内，面法线垂直于面，所以两者在球面三角学上是独立参数。一般的约束是 λ + φ ≥ 90°（Schwarz 不等式）。本工具没有强制这种物理约束，所以你可以输入 (λ=10°, φ=10°) 这样物理上不可能的组合，但这只应该用于教学的敏感性分析，不要用于定量预测。
最后一个误解是"m 大 = 材料弱" 。正确的理解是："在某个取向下，m 大就意味着该取向的单晶在低应力下就会降伏"，这是因为几何关系，不是材料本身变弱。实际工程中，除了这种几何效应外，加工硬化（位错的纠缠）、析出强化（粒子钉扎位错）、晶粒细化（Hall-Petch）等都会提高 τcrss 本身。本工具假定 τcrss 恒定，只显示"在固定强度下，取向对降伏应力的影响"，要记住这一点。