从循环载荷试验中探索背应力与有效应力:程控电子万能试验机的奇妙之旅
先把话说在前面:如果你做材料力学性能测试,尤其是接触过金属材料低周疲劳或者循环变形行为研究,那你大概率听过“背应力”和“有效应力”这两个词。这俩概念听起来很理论,但它实际上是解释包辛格效应、棘轮效应、循环软化/硬化等一系列宏观现象的钥匙。
我最近在一批车用高强钢的循环载荷试验里,认认真真地用程控电子万能试验机跑了一整套对称应变循环和单轴棘轮试验,把背应力(back stress)和有效应力(effective stress)从滞回曲线里分离出来,整个过程踩了不少坑,也把概念从书本上的公式落到了实实在在的数据曲线上。
这篇文章就把整个过程整理出来,包括测试方案怎么定、设备怎么调、数据怎么处理、背应力和有效应力怎么拆解,以及我在实际操作中遇到的那些说明书里绝对不会写的坑。无论你是刚接触循环加载试验的研究生,还是需要评估材料疲劳性能的工程师,这篇内容应该都能给你一些参考。
1. 为什么循环载荷下要把应力拆成“背应力”和“有效应力”?
在单次单调拉伸试验里,我们习惯用屈服强度、抗拉强度、延伸率这些指标评价材料。只要材料没有明显的屈服平台,工程上通常取0.2%残余应变对应的应力作为屈服强度。这个做法简洁高效,但进入循环加载领域之后,问题就变得复杂了。
1.1 包辛格效应:循环加载里第一个绕不开的现象
如果你把同一根试样先拉伸到塑性应变区,然后卸载再反向压缩,你会发现一个有意思的结果:反向压缩的屈服应力明显低于正向拉伸时的流变应力。也就是说,材料经过一次正向塑性变形后,反向变形的阻力变小了。这就是包辛格效应(Bauschinger effect)。
这个现象在单调拉伸里完全看不出来,因为单调拉伸只向一个方向加载,材料内部的位错结构演化是“单向推进”的。而循环加载是反复正反向变形,位错在滑移面上来回运动,材料内部的应力状态分布就变得很不均匀。
我在这批高强钢上跑循环试验时,第一个明显的现象就是正向拉伸到应变幅1.0%再反向压缩,反向屈服应力比正向流变应力低了不少。刚开始以为是设备或者引伸计出了问题,后来确认这是材料真实的力学响应。
1.2 背应力是什么:材料内部的“长程应力场”
要解释包辛格效应,传统塑性力学引入了背应力这个概念。
从物理图像上说,材料发生塑性变形后,位错在晶粒内部运动,但位错运动并不是均匀的。晶界、析出相、第二相颗粒这些障碍物会阻挡位错,导致位错在局部堆积。这些堆积的位错会产生一个长程的内应力场,这个内应力场的方向与外加应力的方向部分相反。
我们可以把这个内应力想象成一根弹簧:当你拉伸材料时,位错堆积产生的内应力会“拉”着原子往回复方向走。当你卸载甚至反向加载时,这个弹簧会帮助你继续产生塑性变形,所以反向屈服应力看起来就变低了。
这个内应力场的方向性特征,决定了它在外加应力不断反转时表现出“方向记忆”。在连续介质力学里,我们用背应力张量X来表征它,用随动硬化(kinematic hardening)来描述背应力的演化。
1.3 有效应力是什么:真正驱动位错运动的那部分力
与背应力相对,有效应力是外加应力扣除掉背应力之后,真正驱动位错越过局部障碍(如晶格阻力、固溶原子、短程障碍)的那部分应力。在率相关塑性理论中,有效应力被认为是决定塑性应变速率大小的关键参数。
可以打个比方:背应力像是一个宏观层面的“趋势领导者”,它代表材料内部的各向异性内应力状态;有效应力则像是“微观的推进器”,直接负责驱动位错在局部范围内运动。两者之和(在单轴应力状态下,代数和)等于外加应力。
所以,把流动应力拆分成:
σ = X + σ*
其中X为背应力,σ*为有效应力。这个分解不是为了数学上的方便,而是为了区分材料循环变形中两种截然不同的物理机制——背应力的演化反映位错堆积与湮灭的净效应,有效应力的变化则反映短程障碍的阻碍状态。
理解了这一点之后,我意识到一个问题:想要系统评估材料的循环变形行为,光看滞回环的宽度和高度是不够的,必须把流动应力拆开,分别跟踪X和σ*随循环周次的演化。
2. 试验平台搭建:程控电子万能试验机的角色与边界
确定了要测背应力和有效应力,接下来就是选设备和搭测试方案。
2.1 为什么选程控电子万能试验机而不是液压伺服疲劳机
这里要说明白,传统的低周疲劳试验通常使用液压伺服疲劳试验机,因为液压机刚度大、频率响应快,适合做应变控制的快速循环。但我这次做的是准静态循环加载(strain rate在10⁻³/s量级),频率极低,重点在于精确控制应变和采集完整的滞回环数据,这种情况下高质量的程控电子万能试验机完全够用,而且成本低很多。
电子万能试验机的特点在于:
- 横梁位移控制精度高,速度范围宽,从每分钟几微米到几百毫米都能精确控制
- 载荷传感器线性度好,适合做低循环次数下的精确载荷测量
- 配合数字控制器可以做三角波、梯形波等复杂波形加载
- 数据采样率对低频循环完全够用
- 试验过程可控性强,随时可以暂停、修改参数
我用的是一台100kN级电子万能试验机,搭配了动态引伸计和全套循环控制软件。在低应变速率条件下,这套设备做50-100周次的循环加载完全没问题。
2.2 试样制备:几何尺寸与表面质量对结果的影响
循环加载试验的试样制备比单调拉伸要严格得多。我这次用的是板状试样,标距段宽度、厚度、过渡圆弧都按照ASTM E606标准来设计。标距段的表面粗糙度直接影响裂纹萌生寿命和应力响应,一般要求Ra不超过0.2μm。
特别要提醒的是打磨方向。循环加载过程中,试样表面的打磨痕迹如果垂直于加载方向,会形成微小的应力集中源,加速表面裂纹萌生。所以打磨时一定要让划痕方向平行于加载方向,最后一道砂纸建议用600目以上的细砂纸。
试样尺寸的均匀性也很关键。如果标距段宽度偏差超过0.02mm,在应变控制模式下会导致局部应变分布不均,得到的应力响应数据误差会明显偏大。我一般在试验前用千分尺在标距段三个位置量宽度和厚度,取平均值作为计算截面积的输入。
2.3 夹持方案与对中调整
板状试样的夹持看起来很简单,但实际陷阱不少。循环加载要求每一周次拉压对称,如果试样装夹不正,压缩时候试样会微弯,导致压缩载荷偏小,数据全废。
我通常的做法是:
- 先用手拧紧上夹头夹具,让试样依靠自重自然下垂
- 用水平仪校准试样侧面,使试样轴线与加载轴线重合
- 用记号笔在试样表面画一条竖直线,观察夹紧过程中试样是否发生扭转
- 将设备载荷清零,用很小的力(比如50N)预拉一次,确认载荷-位移曲线无异常
- 完全夹紧后,再确认一次引伸计的读数稳定
对中问题在拉伸试验里不太明显,但在循环加载的压缩段会暴露得很彻底。如果试样的初始弯曲导致压缩时失稳,得到的压缩应力数据会系统性偏低,背应力计算结果自然也不可靠。
2.4 引伸计安装与标距选择
引伸计是循环应变控制的关键。很多初学者有个误区:以为用横梁位移换算应变就够了。这在弹性段没问题,一旦进入塑性段,试样标距段以外的部分也会发生弹性变形,更关键的是夹持端的滑移会让横梁位移与试样真实应变完全对不上。
所以塑性循环加载必须用引伸计直接测量标距段应变。我这台试验机的引伸计标距是25mm,分辨率为1μm,满足试验需求。安装引伸计的时候要注意两个细节:
- 引伸计的刀刃要对准试样标距段的对称中心,先用橡皮筋轻轻固定,确保刀刃与试样表面垂直接触
- 在循环加载过程中,引伸计刀刃与试样表面会产生微小的摩擦滑动,如果试样表面有油污,刀刃打滑会导致应变信号跳变,所以试样表面要清洁干净
3. 循环加载程序设计与控制参数选择
设备准备好了,接下来就是设计加载程序。这一步的核心问题是:采用什么样的加载路径、应变幅值、加载速率和循环周次,才能有效提取背应力与有效应力。
3.1 应变控制模式与对称三角波路径
我这次的核心试验是应变控制的对称循环试验,加载波形为三角波,应变幅值设定在±0.6%、±1.0%、±1.5%三档。为什么用三角波不用正弦波?
关键原因在于数据处理的方便性。三角波加载时,应变速率在每一段内是恒定值,应力响应不会因为应变速率的变化产生额外的粘性效应。更重要的是,从滞回环的加载段和卸载段提取特征应力时,三角波加载下的匀速变形状态更容易定义“屈服点”。
当然,三角波在每一段的起始和结束位置有加速度突变,会对应变速率产生瞬时扰动。对绝大多数金属材料来说,这种扰动的影响可以忽略。如果材料对应变速率极敏感,可以考虑梯形波(中间平台段保持恒应变速率),但需要增加保载时间,试验周期会变长。
3.2 应变速率的选择依据
应变速率直接影响到材料的粘塑性响应。我这次选的应变速率是2×10⁻³/s,相当于每秒钟应变变化0.2%。在±1.0%应变幅下,一个完整循环(拉至峰值再回到负峰值再回到正峰值)大约需要20秒,50周次循环约17分钟,试验时间完全可控。
为什么不选更快的速率?因为电子万能试验机毕竟不是动态试验机,横梁高速反转时会存在一定的反向间隙和惯性过冲,这些机械因素会叠加到真实的应力信号上,影响数据质量。而准静态速率下,横梁的加减速过程很平缓,应力响应几乎不受机械惯性影响。
3.3 循环周次:多少周次才算够
从背应力和有效应力的演化规律来看,材料通常在前10-20周次内发生明显的循环软化或硬化,之后逐渐进入稳定状态或准稳定状态。我通常做到50周次,已经足够看到背应力从初始状态快速增长到接近饱和的完整过程。
如果条件允许,做100周次自然更好,能更清楚地看到稳定阶段的背应力饱和值。但也要平衡时间成本,特别是应变幅较大的工况,试样可能在几十周次内就萌生裂纹,数据后段的可靠性反而下降。
3.4 棘轮试验作为补充
除了对称循环,我还加做了一组非对称应力控制的棘轮试验:应力均值固定为正,应力幅值低于屈服强度,但平均应力+应力幅值超过屈服极限。这种情况下,每个循环都会产生一次正向塑性应变累积,形成“棘轮效应”。
棘轮试验的好处是,它能在非对称循环条件下进一步验证背应力演化的方向特征:当平均应力不为零时,背应力会不断调整以适应该方向上的塑性应变累积,有效应力的变化则更直接地反映材料对持续塑性流动的阻力变化。
4. 从滞回环中分离背应力与有效应力:数据处理全流程
分离背应力与有效应力的实际操作,核心是找到滞回环上的特征点。这里我把整个数据处理流程完整写出来,包括每一步的物理定义和操作细节。
4.1 滞回环的基本形态与特征应力定义
一次应变控制的对称循环中,应力-应变关系会形成一条闭合滞回环。典型的滞回环有四个关键特征点:
- 峰值拉伸点A(ε = εmax,σ = σmax)
- 峰值压缩点C(ε = εmin,σ = σmin)
- 正向加载屈服点B(拉伸支路上塑性应变增量开始明显增大)
- 反向加载屈服点D(压缩支路上塑性应变增量开始明显增大)
更常用的做法是用0.2%残余应变偏移法确定正向屈服应力和反向屈服应力。以正向加载支路为例:
- 取该支路的线性弹性段斜率作为弹性模量E
- 沿弹性段偏置0.2%应变,画一条平行线
- 平行线与该支路应力-应变曲线的交点即为正向屈服应力σy
反向屈服应力以同样的方式从卸载支路上取。这里要特别注意,正向屈服应力是σy+,反向屈服应力是σy-。
4.2 背应力与有效应力的经典计算式
有了正向和反向屈服应力,背应力和有效应力可以通过以下关系计算:
X = (σy+ + σy-) / 2
σ* = (σy+ - σy-) / 2
关于这两个公式的物理意义,我简单解释一下:
- 如果材料中不存在背应力(X=0),那么正向和反向屈服应力应该对称相等,σy+和σy-大小一样。但实际上由于背应力的“方向记忆”,正向加载时背应力协助位错运动,导致σy+偏低;反向加载时背应力阻碍位错运动,导致σy-偏高(绝对值上)。
- 把两者的平均值看作背应力,差值的一半看作有效应力,这是一个很巧妙的线性分解。
- 这个计算方法的隐含假设是:短程障碍对位错运动的阻力在正向和反向加载时大小相等。对大多数金属材料,这个假设在工程误差范围内成立。
实际上,这个公式也常写作X = (σy+ + σy-) / 2,当σy-取负值时相加再除以2,所得X可以为正也可以为负。在对称循环稳定阶段,X会在一个正值附近摆动(因为拉伸塑性变形历史略多于压缩),所以我没有把X视为常量,而是逐周次计算后取峰值,得到背应力随循环周次的变化曲线。
4.3 数据处理中的取舍:为什么不用峰值应力而用屈服应力
有些人计算背应力时直接用峰值拉伸应力和峰值压缩应力来代替屈服应力。这样做虽然省事,但物理上并不严谨。
峰值应力包含了大量的弹性分量和复杂的塑性演化信息,它并不能反映位错刚开始发生反向运动时的临界状态。而屈服应力(用0.2%偏移法定义)更接近位错重新开始大量运动的临界点,用它来计算背应力和有效应力更有机理意义。
不过,0.2%偏移法在滞回环处理中有个麻烦:弹性模量要从整个加载支路的斜率来拟合。如果弹性模量取值偏大或偏小,屈服应力的位置会明显漂移。我的经验是至少取加载支路的前30%的应力-应变数据点做线性拟合,并且剔除掉首尾的数据点,以减少引伸计位移起步阶段的非线性影响。
4.4 从峰值应力到全曲线分解
仅仅计算屈服点处的背应力和有效应力还不够,我更想看到整个滞回环上每一点的背应力-应变关系。这需要用到经典的非线性随动硬化框架下的数据处理方法。
一种实用的做法是:在滞回环正向加载段,对每一个数据点,将当前应力减去相邻反向加载段对应应变处的应力再除以2,近似得到该应变下的背应力。但这个方法要求两个加载段的应变路径完全对齐,实际操作时要用线性插值将不同循环的应变坐标统一。
另一种更简单但更工程化的方法:
- 将滞回环按应变坐标等间隔离散为N个点(我一般取200个点)
- 取加载支路和卸载支路上相同应变位置的应力值之差的一半作为该应变水平的背应力增量贡献
这样做的好处是能得到背应力-应变曲线,观察背应力在塑性变形过程中如何随应变增加而演化。对于建立材料本构模型,这条曲线的饱和趋势和变化速率直接决定了随动硬化模型的参数选择。
5. 试验结果分析:背应力与有效应力的演化规律
经过50周次对称循环和几十组棘轮试验数据提取,我把这些曲线画在一起对照分析。下面把几组典型的、有代表性的结果规律写出来。
5.1 初期循环软化与背应力的快速调整
这批高强钢在初始单调拉伸时表现得很“硬”,屈服强度大约在680MPa。但在第一周次循环加载后,应力峰值开始明显下降,呈现典型的循环软化特征。这个软化过程主要对应背应力的重新分布:初次循环时,位错结构从单调塑性变形模式(方向单一)转向循环模式(正反向反复滑移),原本单向堆积的位错墙和胞状结构被部分拆解,背应力迅速下降到新的平衡水平。
从数据上看,第一周次到第五周次之间背应力从初始值约150MPa下降到120MPa左右,降幅约20%。有效应力则相对稳定,从260MPa微降到255MPa。这说明循环软化主要是背应力贡献的下降,而非短程障碍阻力的变化。这一结果让我对“循环软化来源”有了更直接的认识——很多文献里说循环软化与位错重排有关,落在背应力和有效应力的分解上,就是指背应力下降了。
5.2 背应力随循环周次的饱和行为
在10周次之后,背应力基本进入缓慢变化阶段。以±1.0%应变幅的试验为例,从第10周到第50周,背应力从118MPa缓慢增加到124MPa,增量和噪声已经接近。这说明材料内部长程内应力场已经达到某种动态平衡:正向加载时位错堆积产生的内应力增量,与反向加载时位错湮灭产生的内应力松弛基本持平。
不同应变幅下背应力的饱和值有明显差异。我对比了±0.6%、±1.0%、±1.5%三组试验数据:
| 应变幅 | 第10周背应力 | 第50周背应力 | 饱和趋势 |
|---|---|---|---|
| ±0.6% | 82MPa | 85MPa | 早饱和 |
| ±1.0% | 118MPa | 124MPa | 缓慢上升 |
| ±1.5% | 152MPa | 163MPa | 仍上升 |
这个趋势和随动硬化的饱和机制是对应的:背应力的饱和值取决于位错结构的稳定间距,而这个间距在小应变幅时细密且稳定,在大应变幅时粗大且演化更活跃,所以背应力饱和值更高且饱和过程更慢。
5.3 有效应力的“黏滞”特征
有效应力σ*在整个循环过程中的变化幅度远小于背应力。三组试验中,有效应力都在245-265MPa范围内窄幅波动,且随循环周次的演化趋势不明显。
这说明有效应力主要取决于材料本征的晶格摩擦力和短程障碍强度,而这些微观特征在循环过程中虽然也会发生变化(比如位错密度升高会增加局部阻力),但净效应远不如长程内应力的调整那么显著。
有效应力的另一个重要特征是速率敏感性。我尝试在同一个试样上改变加载速率(从2×10⁻³/s提高到5×10⁻³/s),有效应力立即增大了约8MPa。这个现象符合热激活位错运动理论:应变速率增加意味着位错运动需要更大的有效应力来克服短程障碍。把背应力和有效应力区分开以后,这种速率效应就能干净地归因到有效应力头上,不会和背应力的演化耦合在一起。
5.4 滞回环形状变化:从矩形到梭形的信息量
随着循环进行,滞回环的形状也在变化。初始循环的滞回环面积较大,几乎呈“矩形”;稳定循环后的滞回环面积缩小,更接近“梭形”。
滞回环面积的减小对应的是循环软化:同样应变幅下,材料在稳定循环中耗散掉的塑性功更少。我用图像积分的方法计算了每个循环的滞回环面积,发现第一周的滞回环积分能量密度约为3.2MPa(应变单位下),第10周下降到约2.1MPa,之后基本稳定。
这个滞回环形状变化也影响了背应力的判定。第一周循环时,由于包辛格效应明显,正向和反向屈服应力的差值较大,有效应力计算值相对高。稳定循环后,包辛格效应减弱,有效应力计算值回落。这提醒我一点:如果做的是低周疲劳损伤评估,滞回环面积随循环周次的变化本身就是材料损伤程度的重要判据。
6. 实操中那些最容易被忽视的“坑”与应对经验
这一节写我在整个过程中遇到的实际问题,尤其是那些说明书上不会写、但直接影响数据可靠性的细节。
6.1 引伸计的“零点漂移”与循环累积
引伸计零点漂移在循环加载中是个比较隐蔽的问题。因为循环试验时间长,温度变化、引伸计刀刃活动部分的摩擦、电子信号温漂都会造成零点漂移,表现为滞回环不闭合或者整体向应力正方向偏移。
应对措施:
- 试验前给引伸计充分的热机时间(至少30分钟),让电子元件温度稳定
- 试验中间若发现滞回环明显不闭合,要紧固刀刃固定螺丝,必要时重新标定
- 处理数据时以每个循环的峰值压应力点作为参考,将零漂扣除
这些小操作对单调拉伸没什么影响,但在需要精确计算正向反向屈服应力的循环试验中,是决定数据是否可信的底线。
6.2 夹持滑移的检测与处理
电子万能试验机的楔形夹具在夹紧试样后,理论上不会滑动。但循环加载的反复拉压会造成夹持面的微松动,特别是在高应力段。
我测试时的一个有效检测方法是:在每个指定循环(比如第5、10、20、50周)暂停试验,做一次小载荷下的加载-卸载循环(载荷不超过屈服载荷的30%),观察横梁位移与引伸计应变的关系。如果横梁位移增加但引伸计应变不变,说明夹持段存在滑移。这时候要松开夹具重新装夹,并且检查试样夹持段是否有压痕。
6.3 应变幅超调与三角波拐点过冲
在应变控制模式下,三角波在方向反转时会有一个轻微的“过冲-回调”过程。这个现象在控制器PID参数不当的时候特别明显。应变超调会导致实际应变幅大于设定值,直接影响背应力计算。
我的处理方式是:
- 在正式试验前用同批次废样做一次调试,观察控制器的应变跟踪效果
- 适当降低PID的P值,增大I值,可以有效抑制过冲(具体参数因设备而异)
- 数据整理时排除每个加载支路前后各1%的数据点,减少拐点非线性区对屈服点拟合的干扰
6.4 试样表面温度变化带来的“假强化”
循环塑性变形会产生热量。在快速循环时,试样温度上升会改变材料的力学响应,表现为明显的循环软化或应力下降。虽然我们的准静态加载产热不大,但在高应变幅(±1.5%)下,50周次循环后试样表面温度仍然上升了约5-8℃。
对于高应变速率试验,试样温升会极其严重。如果必须在较高应变速率下做试验,建议在试样表面贴热电偶监测温度,确保温度变化不改变材料的本征响应。如果观察到明显的温升,应该降低加载速率或者增加循环间的间歇时间。
6.5 数据采样率与滤波截止频率的选择
数据采样率太低,滞回环的尖角处会丢失细节;采样率太高,数据文件巨大且噪声明显。我的经验是每循环至少采200个点,同时记录应力和应变。对于±1.0%应变幅、三角波加载,采样率设为20Hz即可满足要求。
后处理时如果应力信号噪声较大,可以用滑动平均或Butterworth低通滤波。但滤波截止频率必须设置合理,太低会滤掉真实的材料响应(尤其是屈服点附近应力的快速变化)。我一般是先画原始数据曲线,确认噪声水平,再做滤波。
7. 背应力与有效应力分解结果对工程评估的启发
试验做到最后,数据摆在那里,其实给人的启发已经超出了单一试验本身的范畴。这里聊聊我的一些体会。
从材料本构建模的角度看,背应力和有效应力的分离结果直接对应随动硬化与各向同性硬化参数的标定。以Chaboche模型为例:
- 背应力的演化用非线性随动硬化律描述,其中的饱和值C/γ可以由背应力随循环周次的饱和行为直接标定
- 有效应力代表的各向同性硬化部分则可以用滞回环宽度、峰值应力的缓慢演化来校准
对于工程上的疲劳寿命估算,背应力演化数据可以用来解释平均应力松弛效应。非对称循环下,背应力在主应力方向上累积,导致平均应力逐渐松弛到某个稳定值,这对疲劳寿命的预测非常重要。传统基于峰值应力的疲劳模型往往预测偏保守,而把背应力纳入考虑后,可以更准确地评估循环软化材料的真实寿命。
对材料研发人员来说,背应力占比(X/(X+σ*))是一个很有价值的评价指标。背应力占比高的材料,循环变形时更容易表现出包辛格效应,在成形-服役一体化工艺中需要特别关注回弹和疲劳性能的关联。
这次试验还让我对设备选型有了更明确的判断:手头没有液压疲劳机的时候,程控电子万能试验机配合低应变速率循环控制,完全可以把背应力与有效应力的演化规律吃透。即使是带动态功能的电子万能试验机,做低频率、低周次的循环载荷试验,与纯粹的疲劳试验机相比,在数据精度与成本控制之间能拿到一个很好的折中。
试验做完后我把方法和结果整理成了内部报告,也给实验室的同事分享了一些数据处理上的经验。如果以后再有人问我背应力到底是“什么东西”,我会直接说:去跑一次循环载荷试验,从滞回环里把它算出来。比任何理论描述都直观。