光弹性应力分析 返回
实验应力分析

光弹性应力分析 — 等色条纹图案可视化

改变试件形状、荷载和条纹常数,实时可视化与主应力差对应的等色条纹(等色线)图案。体验FEM验证的基础原理。

试件与荷载条件
试件形状
荷载 P
N
试件径/宽 D
mm
厚度 t
mm
条纹常数 f
N/mm
计算结果
计算结果
σ₁ (MPa)
σ₂ (MPa)
τmax (MPa)
条纹级数 N
等色条纹图案(Canvas)
荷载轴上的应力分布
理论与主要公式

条纹级数:$N = \dfrac{(\sigma_1 - \sigma_2)\,t}{f}$

圆盘(赫兹):
$\sigma_y = -\dfrac{2P}{\pi D t},\quad \sigma_x = \dfrac{6P}{\pi D t}$

主应力差越大的区域条纹越密。

什么是光弹性应力分析

🙋
光弹性法是什么?听起来好神奇,用光就能看到应力?
🎓
简单来说,就是把一种特殊的透明塑料模型放在两个偏光片中间,然后给它加力。当光穿过这个受力模型时,模型内部因为应力变成了“光学各向异性”,导致光的传播速度不一样,最后就形成了彩色的条纹图案。这些条纹就像应力的“等高线图”,条纹越密的地方,应力变化越剧烈。你可以在模拟器里试着拖动“荷载P”的滑块,马上就能看到条纹图案变密,这就是应力增大的直观表现。
🙋
诶,真的吗?那这些彩色的条纹,我怎么知道哪条对应多大的应力呢?
🎓
问得好!这就需要用到“条纹级数”这个概念了。每一条颜色相同的线,代表一个固定的“主应力差”值。比如,从黑色的0级条纹开始,接下来出现的红色、黄色、绿色……分别对应1级、2级、3级条纹。在实际工程中,比如分析一个带圆孔的平板,我们最关心的就是孔边那条最密集的彩色环,那通常就是应力最大的地方。你可以在模拟器里把“试件形状”从“圆盘”切换到“带孔平板”,马上就能看到孔边出现密集的彩色环,这就是典型的应力集中现象。
🙋
原来如此!那这个方法和我们学的有限元模拟(FEM)有什么关系?它算出来的结果准吗?
🎓
这正是光弹性法的核心价值之一——它是验证FEM结果的“黄金标准”实验方法。工程现场常见的是,先用FEM软件算出一个零件的应力云图,然后做一个光弹性模型来实测对比。如果两者的条纹(或等值线)图案吻合得很好,就说明FEM的模型和边界条件设置是可靠的。你可以在模拟器里试试看:选择一个形状,然后改变“条纹常数f”,你会发现整个图案的疏密会整体变化,这其实就相当于在调整材料的灵敏度,帮助你理解实验参数如何影响最终的观测结果。

物理模型与关键公式

光弹性法的核心是应力-光学定律,它建立了材料内部的光程差(表现为条纹)与主应力差之间的定量关系。这是将光学现象转化为力学数据的桥梁。

$$N = \dfrac{(\sigma_1 - \sigma_2)\,t}{f}$$

$N$:条纹级数(无量纲),从0开始的整数或小数。
$\sigma_1 - \sigma_2$:最大与最小主应力之差(MPa)。
$t$:试件厚度(mm)。
$f$:材料条纹常数(N/mm/条),是材料的固有属性,表示产生一级条纹所需的主应力差与厚度的乘积。

对于经典的对径受压圆盘试件(赫兹接触问题),其中心点(应力最大处)的应力有理论解。这个解常被用来标定材料条纹常数 $f$,或者验证光弹性实验的准确性。

$$\sigma_y = -\dfrac{2P}{\pi D t},\quad \sigma_x = \dfrac{6P}{\pi D t}$$

$P$:施加的集中荷载(N)。
$D$:圆盘直径(mm)。
$t$:圆盘厚度(mm)。
中心点主应力差 $\sigma_x - \sigma_y = \dfrac{8P}{\pi D t}$,代入条纹公式即可求得该点的理论条纹级数。

现实世界中的应用

汽车与航空航天结构验证:在发动机连杆、飞机起落架等关键承力部件设计初期,会制作光弹性模型来观察应力分布,特别是圆角、螺栓孔等易产生应力集中的部位,确保FEM分析结果的可靠性,避免疲劳破坏。

土木工程与建筑节点分析:用于分析复杂的钢结构节点、混凝土预制件的连接处。通过光弹性条纹,可以直观地看到力流是如何在节点中传递的,从而优化焊缝位置或加强筋的布置。

微电子与封装应力检测:芯片封装过程中,不同材料的热膨胀系数不匹配会产生热应力。利用透明模型模拟封装结构,通过光弹性法可视化热应力分布,对防止芯片开裂或焊点失效至关重要。

生物力学与医疗器械设计:例如分析人工关节(如髋关节假体)与骨骼接触面的应力分布。通过光弹性模型可以优化假体形状,使应力传递更均匀,减少对骨骼的损害和植入物的松动。

常见误解与注意事项

首先,人们常认为“条纹密集=应力大”,但严格来说这是不准确的。等色条纹反映的是主应力“差”,例如均匀拉伸(σ1较大,σ2=0)与均匀剪切(σ1和σ2绝对值相同但符号相反)虽然应力差相同,但应力状态完全不同。在模拟器中选择“圆盘”并施加载荷,观察中心区域:出现的条纹是由拉伸和压缩组合产生的大应力差所致。而简单拉伸试件在加载初期几乎看不到条纹(因σ2接近0,应力差较小)。请始终牢记,我们观察的是应力的“差异”而非“绝对大小”。

其次,关于条纹灵敏度常数“f”的处理。该值为材料固有属性,例如环氧树脂约为3~5 kN/m·条纹,丙烯酸树脂约为10~15 kN/m·条纹,不同材料差异显著。在模拟器中更改此值,即使载荷相同,条纹数量也会大幅变化,对吗?实践中,不能将在已知材料(如环氧树脂)实验中获得的知识直接套用于其他材料(如铝合金)的设计,因为f值的差异会改变“表观危险程度”。

最后,必须理解该方法“仅能评估二维面内应力”的根本局限。该技术的前提是光透过方向的应力沿厚度恒定(平面应力状态)。实际厚壁构件中,表面与内部可能产生应力状态不同的“三维应力”。因此,该方法不适用于复杂三维形状的精确评估,其定位应限于趋势把握或有限元分析验证。

使用指南

  1. 在左侧面板设置试件几何参数:输入圆盘直径(vDNum单位mm)、厚度(vTNum单位mm)、材料条纹常数(sP单位MPa/条纹),选择荷载类型(径向压缩/偏心拉伸/弯曲)
  2. 调整荷载大小(vPNum单位kN)和施加位置参数(sPNum),实时观察等色条纹图案变化,监测σ₁、σ₂、τmax输出值的变化趋势
  3. 通过条纹级数N与实测应力差验证公式N=(σ₁-σ₂)/f_t(f_t为条纹常数),改变材料或荷载重复计算,建立条纹密度与应力集中的对应关系

具体计算示例

圆盘直径D=50mm、厚度t=10mm、条纹常数f_t=5.2MPa/条纹、径向压缩荷载P=8kN的情况:根据圆盘接触应力理论,接触应力最大处σ₁≈120MPa、σ₂≈-85MPa,则主应力差Δσ=205MPa,条纹级数N=205÷5.2≈39.4级。模拟器显示圆心附近出现密集的等色条纹,外围区域条纹稀疏,τmax≈102.5MPa。若增加荷载至P=12kN,σ₁升至165MPa,条纹级数增至56级,圆心区域条纹密度显著增加,清晰展现应力集中现象。

实务注意事项

  1. 偏心拉伸荷载(sPNum≠圆心)会产生非对称的等色条纹分布,试件一侧条纹密集(高拉应力)、另一侧稀疏,用于验证应力分布的不对称性
  2. 条纹常数f_t随材料选择变化:玻璃纤维复合材料f_t≈3.5MPa/条纹、环氧树脂f_t≈6.8MPa/条纹,相同荷载下条纹级数反比,直接影响实验分辨率
  3. 厚度参数(vTNum)变化会改变应力分布深度,薄片(t<5mm)易产生边界效应导致条纹畸变,应选择t≥8mm确保二维应力状态有效性
  4. 荷载施加位置精度控制在±0.5mm以内,否则引入附加弯矩导致实验结果与经典圆盘理论偏差超过15%