一种获取金属材料在静态压缩状态下应力应变曲线的方法与流程

本发明涉及一种金属材料性能检测领域，涉及一种获取金属材料在极限慢压缩即静态压缩状态下的不同温度环境下的真实应力应变曲线的方法。

背景技术

金属材料在长期承受压缩载荷服役状态下，尤其在高温环境下，因为材料强度的降低，往往会引发结构失效，造成巨大的经济损失，该方法可以有效的评价材料在该状态下性能的衰减。对材料的安全服役具有重要的指导意义。

在本发明之前，金属材料进行压缩实验，实验加载速率一般保持在10^-5s^-1以上，无法有效模拟长期受力的使用情况，而且压缩加载速率越快检测到的材料性能越高，与材料在长期承受压缩在状态下的真实材料性能差距越大。松弛实验可以测得单点的真实应力情况，但是其无法得到完整材料拉伸曲线，而且松弛实验耗费时间很长，实验时间成本较高；同时尤其对于各项异性材料，获取材料在拉伸状态下的应力应变曲线，并不能代表材料的压缩性能，因此获取极限慢压缩即静态压缩状态下的金属材料真实应力应变曲线具有重要意义。

技术实现要素：

本发明的目是提供一种获取金属材料在极限慢压缩即静态压缩状态下的不同温度环境下的真实应力应变曲线的方法，该方法可以得到金属材料在恒压缩服役状态下的真实应力应变曲线，同时可以比松弛实验节约大量时间成本。对长期处于压缩服役状态下的金属材料的使用具有重要的指导意义。

为实现上述目的，本发明提供的一种获取金属材料在极限慢压缩即静态压缩状态下的不同温度环境下的真实应力应变曲线的方法包括以下步骤：

1、改进型压缩实验方法，该方法可以获取后续处理分析所用的数据，具体包含一下步骤：

1)实验条件：首先，压缩实验机必须具备配套实验环境温度箱和可以应用在对应环境温度下的引伸计，如不具备此功能，该方法只可以获取常温下的实验数据；其次压缩实验机必须具备应变控制或者位移控制功能，如果不具备应变控制或者位移控制，该实验无法进行，推荐采用gleeble3500等热模拟试验机进行该试验。

2)制备金属材料压缩实验用样，并将实验用样加热到预定实验温度，并在整个实验过程中保持实验温度恒定不变。

3)在一定载荷速率情况下进行压缩，当压缩到材料的应变数值在-0.5％、-1.0％、-1.5％、-2.0％、-2.5％、-3.0％、-3.5％、-4.0％处进行应变控制或者位移控制，保持在该应变30分钟到2小时不等，采集该过程中所有应力、应变、时间数据。

2、实验数据处理，该方法可以获取金属材料在极限慢压缩即静态压缩状态下的不同温度环境下的真实应力应变曲线，具体实验步骤如下：

1)将从上一步实验数据导出应力应变曲线、应力时间曲线。

2)将工程应变转化为常规压缩状态下的真实应变。

3)将工程应力转化为常规压缩状态下的真实应力。

4)分别对-0.5％、-1.0％、-1.5％、-2.0％、-2.5％、-3.0％、-3.5％、-4.0％应变处时间应力曲线进行数据拟合求解。

5)根据拟合求解得出公式，求解时间对应无穷大情况下对应的应力值σs。

6)将不同应变点对应的σs连接绘制金属材料在极限慢压缩即静态压缩状态下的不同温度环境下的真实应力应变曲线。

本发明的效果是，通过一种获取金属材料在极限慢压缩即静态压缩状态下的不同温度环境下的真实应力应变曲线的方法，能够精确的获得金属材料在不同温度环境下承受恒定压缩时金属材料的实际应力应变曲线。本发明得到的材料性能曲线，对长期处于压缩服役状态下的金属材料的使用具有重要的指导意义。

附图说明

图1是由压缩实验机直接获取的工程压缩应力应变数据曲线；

图2是通过计算得出的真实应力应变曲线；

图3、图6、图9分别是-0.5％、-2.0％、-3.5％应变处时间应力曲线以及数据拟曲线；

图4、图7、图10分别是-0.5％、-2.0％、-3.5％应变处时间应变曲线；

图5、图8、图11分别是-0.5％、-2.0％、-3.5％应变处时间温度曲线；

图12是极限慢压缩情况下对应的真实应力应变曲线；

图13是工程压缩应力应变数据曲线、真实应力应变曲线和极限慢压缩下的应力应变曲线对比图。

具体实施方式

结合附图对本发明的一种获取金属材料在极限慢压缩(静态压缩)状态下的不同温度环境下的真实应力应变曲线的方法加以说明

本发明的一种获取金属材料在极限慢压缩即静态压缩状态下的不同温度环境下的真实应力应变曲线的方法包括一下步骤：

1、改进型压缩实验方法，该方法可以获取后续处理分析所用的数据，具体包含一下步骤：

1)实验条件：首先，压缩实验机必须具备配套实验环境温度箱和可以应用在对应环境温度下的引伸计，如不具备此功能，该方法只可以获取常温下的实验数据；其次压缩实验机必须具备应变控制或者位移控制功能，如果不具备应变控制或者位移控制，该实验无法进行。

2)制备金属材料压缩实验用样，优选采用圆棒试样，试样长度在保证能够正常试验的前提下尽量缩短，可以保证试样在压缩过程中的稳定性，并将实验用样加热到预定实验温度，并在整个实验过程中保持实验温度恒定不变，如图5、图8、图11时间温度曲线所示，在整个实验保温过程中应保持试验温度温差不超过±3℃，如果温差较大试验结果不建议采用。压缩实验应变加载速率保持在1.0×10^-5s^-1左右，试验加载速率越慢试验越稳定，得到的试验结果越接近真实值。

3)在一定载荷速率情况下进行压缩，当压缩到材料的应变数值在-0.5％、-1.0％、-1.5％、-2.0％、-2.5％、-3.0％、-3.5％、-4.0％处进行应变控制或者位移控制，如图4、图7、图10时间应变曲线所示，在每个保载时间内控制应变波动在0-0.01％范围内，如果应变波动超出上限，试验结果不建议采用，保持在该应变30分钟到2小时不等，优选时间为1小时，时间越长最终应力数值越稳定，但由于考虑到时间成本和实验设备承受能力，最终推荐时间1小时，采集该过程中所有应力、应变、时间数据。

2、实验数据处理，该方法可以获取金属材料在极限慢压缩即静态压缩状态下的不同温度环境下的真实应力应变曲线，具体实验步骤如下：

1)将从上一步实验数据导出应力应变曲线，如图1所示。

2)将工程应变转化为常规压缩状态下的真实应变。压缩试样在加载过程中存在由于压缩变形的情况(此情况是完全无法避免的)，同时实验设备采集到应力应变数据是没考虑这种变形存在的前提下的一种数据，因此称之为工程应力应变数据，所以要通过计算得到材料的真实应力应变数据，真实应变s真实与工程应变s工程的计算公式如下：

s真实＝s工程*(1+σ工程)

3)将工程应力转化为真实应力，真实应力σ真实与工程应力σ工程的计算公式如下：

σ真实＝ln(1+σ工程)*100

通过数据计算的到真实应力与真实应变数据，将真实应力与真实应变进行数据绘制得到真实应力应变曲线，如图2所示。

4)分别对-0.5％、-1.0％、-1.5％、-2.0％、-2.5％、-3.0％、-3.5％、-4.0％等应变处时间应力曲线进行数据拟合求解。

由于材料在保持位移或者应变恒定的情况下，随着时间的增长存在载荷降低的情况，而且随着时间增长载荷降低的趋势会逐渐降低，载荷值最终会趋近于一个恒定不变的较小稳定数值。图3蓝色曲线记录了在-0.5％应变处的时间真实应力曲线。

由于试验时间、成本、设备等的因素的考虑，试验时间不能过长，因此需要通过计算获取对应时间下的载荷值，因此根据曲线特点采用如下公式进行数据计算：

σ真实＝σ最终+a*exp(-t/b)+c*exp(-t/d)

式中σ真实是通过上步计算得到的各个点对应的真实应力值；σ最终是曲线最终趋近的稳定最小真实应力值，也是需要通过数据计算得出的常数值；a、b、c、d是需要通过数据计算得到的常数值；t是对应曲线里面的时间值。

利用以上公式进行数据计算计算得到如下数值与公式：

σ最终＝435.92，a＝141.85，b＝105.30，c＝15.42，d＝993.2。

σ真实＝435.92+141.85exp(-t/105.30)+15.42exp(-t/993.2)

该公式对应的曲线如图3所示，同时在该应变情况下材料长期服役的材料的真实性能σ最终＝435.92mpa。

5)根据拟合求解得出公式，求解时间对应无穷大情况下对应的应力值σs。

根据上步拟合求解得到的公式，时间t带入无限大，求得

σs＝σ最终＝435.92mpa

6)将不同应变点对应的σs连接绘制金属材料在极限慢压缩(静态压缩)状态下的不同温度环境下的真实应力应变曲线，如图12所示。

图13是图1、图2和图12分别所示的工程压缩应力应变数据曲线、真实应力应变曲线和极限慢压缩下的应力应变曲线对比图。

下面给出四个具体计算实例。

实施例1

如图6所示是压缩试样在应变-2.0％处的真实应力应变曲线，利用该公式σ真实＝σ最终+a*exp(-t/b)+c*exp(-t/d)对试验获取的真实应力应变数据进行数据拟合，得到如下结果：

σ最终＝474.34，a＝6.39e²⁹，b＝72.33，c＝12365.3，d＝657.71

σ真实＝474.34+6.39e²⁹*exp(-t/72.33)+12365.3*exp(-t/657.71)

拟合曲线如图6所示，根据上步拟合求解得到的公式，时间t带入无限大，求得

σs＝σ最终＝474.34mpa

实施例2

如图9所示是压缩试样在应变-3.5％处的真实应力应变曲线，利用该公式σ真实＝σ最终+a*exp(-t/b)+c*exp(-t/d)对试验获取的真实应力应变数据进行数据拟合，得到如下结果：

σ最终＝527.54，a＝2.80e³²，b＝128.92，c＝6120.9，d＝1402.10

σ真实＝527.54+2.80e³²*exp(-t/128.92)+6120.9*exp(-t/1402.10)

拟合曲线如图9所示，根据上步拟合求解得到的公式，时间t带入无限大，求得

σs＝σ最终＝527.54mpa。