技术领域本发明涉及一种高温蠕变预测的方法，特别涉及一种基于短时高温蠕变数据预测材料高温蠕变变形的预测方法。
背景技术：
高温蠕变变形是指材料在高温和低于材料宏观屈服极限的应力下发生的缓慢的塑性变形，其描述的是材料抵抗高温变形的能力。耐热钢在高温高压的服役环境中，其稳定性是高温结构件安全运行的重要保障。高温蠕变变形是导致高温构件失效主要形式之一。因此需要对耐热钢的蠕变变形行为进行高效准确的分析和预测。传统的高温蠕变变形行为预测方法有theta投影法等，是通过高应力短时的蠕变曲线来外推低应力长时间的蠕变曲线。但是传统方法存在诸多弊端，如实验材料本身性能具有分散性，此材料性能分散性会体现在材料的蠕变曲线中，进而对拟合方程中的参数产生影响，最终影响预测的精度。
技术实现要素：
本发明的目的在于，解决现有材料高温蠕变变形预测结果受材料性能分散性的影响的问题，提高预测结果的准确性和稳定性，降低高温蠕变变形预测的波动性。为达到上述目的，本发明提供了一种高温蠕变变形预测方法，包括如下步骤：S1、目标样品短时测试：在设计温度Tj(j＝1、2、3…n＇，n＇≤4)及设计应力σk(k＝1、2、3…n＂，n＂≤4n＇)条件下，对金属测试样品进行蠕变测试，记录测试过程中，所述测试样品蠕变应变值为0.1～20％时，所对应的测试时间；所述设计温度的取值范围为400～1200℃，所述设计应力的取值范围为1MPa～1000MPa；S2、按由小到大的顺序选取步骤S1测得的蠕变应变值εi(i＝1,2,…x，x≤5)，记录步骤S1的各测试中，相同蠕变应变值εi所对应的测试时间tijk，根据公式(1)计算得到综合参数Qijk：Qijk＝10-3×Tj×(logtijk+c)(1)其中，tijk表示在设计温度为Tj、设计应力为σk的蠕变测试中，蠕变应变值为εi所对应的测试时间；c为材料常数；Qijk表示在设计温度为Tj、设计应力为σk的蠕变测试中，蠕变应变值为εi时的综合参数；将设计应力σk与综合参数Qijk形成的数据点(σk,Qijk)，按蠕变应变值εi由小到大的顺序，制成σ-Q关系图；S3、建立σ-Q关系模型：σ＝ai+b×Q+c×Q2+d×Q3：(2)其中，σ为蠕变应力，ai为与蠕变应变εi相关的常数；b,c,d为与蠕变应变εi不相关常数；S4、将步骤S3得到的σ-Q关系模型用于步骤S2得到的σ-Q关系图的数据进行线性或者非线性的拟合，得到一组相互平行的σ-Q关系拟合线。以蠕变应变值为ε1的拟合线作为基准线，其它蠕变应变值εi的拟合线相对于所述基准线的垂直偏移量即为Zci：Zci＝ai-a1(3)；S5、对不同蠕变应变值εi与步骤S4得到的相应的Zc进行拟合，得到Zc-ε拟合方程；S6、根据步骤S5得到的Zc-ε拟合方程，推导蠕变应变值为εm(m＝1、2、3…n)时偏移量Zcm值；根据公式(1)-(3)反推，得到目标温度、目标应力、蠕变应变值为εm所对应的蠕变时间，进而实现目标样品高温蠕变变形预测。优选方式下，步骤S1所述的蠕变测试过程中，测试温度波动小于1℃。优选方式下，为方便快速得到步骤S1中各测试、同一蠕变应变值所对应的测试时间，将步骤S1中设计温度Tj、设计应力σk条件下获得的测试结果绘制成蠕变应变值-测试时间曲线，获得包含n＇×n＂条蠕变实验曲线的蠕变应变值-测试时间图。优选方式下，步骤S2还可以按公式(4)计算得到综合参数Qijk：Qijk＝logtijk+c×Tj(4)。本发明的有益效果在于：1、本发明提出了一种全新的蠕变变形预测方法；据此方法，可以利用已知一定应力下的蠕变实验，来预测已知应力附近的，所求应力下蠕变实验的蠕变变形随时间的变化。2、本发明方法基于短时实验结果即可进行预测，缩短了实际实验过程中的实验时间，能够在相同的实验进度时间内，以更加接近真实环境的的温度，和应力条件操作蠕变实验，提高了预测的准确性，和说服力。3、本发明适用于所有金属材料，有效的解决了现有技术中测试结果波动性大的问题，本发明方法能够在材料性能分散性大的情况下，得出一个合理的预测的结果，提高了预测结果的准确性和稳定性。附图说明图1为实施例1典型蠕变曲线，黑线表示本预测方法中所选取的应变εi在蠕变曲线中的实际位置。图2是实施例1将数学模型σ＝1250.8+Zc-34.4×Q应用于各个应变下的σ-Q关系图。图3是实施例1中Zc参数与蠕变应变ε的关系图。图4是实施例1中蠕变应力为1％时预测时间t*m与真实蠕变实验时间texp的对比。图5是实施例1中蠕变应力为4％时预测时间t*m与真实蠕变实验时间texp的对比。图6为实施例2典型蠕变曲线，黑线表示本预测方法中所选取的应变εi在蠕变曲线中的实际位置。图7是实施例2将数学模型σ＝1318.9+Zc-31.4×Q应用于各个应变下的σ-Q关系图。图8是实施例2中Zc参数与蠕变应变ε的关系图。图9是实施例2中蠕变应力为1％时预测时间t*m与真实蠕变实验时间texp的对比。图10是实施例2中蠕变应力为4％时预测时间t*m与真实蠕变实验时间texp的对比。图11为实施例3典型蠕变曲线，黑线表示本预测方法中所选取的应变εi在蠕变曲线中的实际位置。图12是实施例3将数学模型σ＝285.7+Zc-10.14×Q应用于各个应变下的σ-Q关系图。图13是实施例3中Zc参数与蠕变应变ε的关系图。图14是实施例3中蠕变应力为1％时预测时间t*m与真实蠕变实验时间texp的对比。图15是实施例3中蠕变应力为5％时预测时间t*m与真实蠕变实验时间texp的对比。图16为实施例4典型蠕变曲线，黑线表示本预测方法中所选取的应变εi在蠕变曲线中的实际位置。图17是实施例4将数学模型σ＝320.7+Zc-11.05×Q应用于各个应变下的σ-Q关系图。图18是实施例4中Zc参数与蠕变应变ε的关系图。图19是实施例4中蠕变应力为1％时预测时间t*m与真实蠕变实验时间texp的对比。图20是实施例4中蠕变应力为5％时预测时间t*m与真实蠕变实验时间texp的对比。图21为实施例5典型蠕变曲线，黑线表示本预测方法中所选取的应变εi在蠕变曲线中的实际位置。图22是实施例5将数学模型σ＝2809.3+Zc-100.1×Q应用于各个应变下的σ-Q关系图。图23是实施例5中Zc参数与蠕变应变ε的关系图。图24是实施例5中蠕变应力为0.1％时预测时间t*m与真实蠕变实验时间texp的对比。图25是实施例5中蠕变应力为0.5％时预测时间t*m与真实蠕变实验时间texp的对比。具体实施方式以下实施例用于说明本发明，但不用来限制本发明的范围。在不背离本发明精神和实质的情况下，对本发明方法、步骤或条件所作的修改或替换，均属于本发明的范围。实施例1本发明中一种基于短时蠕变实验的高温蠕变变形预测方法具体方案由以下实施例详细给出。选取P92，开展实验和寿命预测，具体包括以下步骤：1)分别于575℃1下，235MPa8；600℃2下，215MPa6,225MPa7；650℃3下，125MPa1,135MPa2,150MPa3,165MPa4,175MPa5；进行蠕变实验，获的如图1所示得蠕变实验数据。实验过程中温度波动小于1℃；所述的蠕变实验数据的应变范围为0.1％-20％，所述的蠕变实验数据为试样变形值及试样变形值所对应的时间。2)将步骤1)中每个蠕变实验中的蠕变实验数据按蠕变应变从小到大选取不同的蠕变应变，即：0.5％1，1％2，4％3，7％4，10％5蠕变应变，如图1所示；由蠕变应变所对应的时间,获得每个蠕变实验下数据点(0.5％1，12.2118)、(1％2，43.4218)、(4％3，387.5318)、(7％4，495.8418)、(10％5，519518)；(0.5％1，1.6126)、(1％2，9.9226)…(εi，tijk)，其中tijk表示温度为Tj(j＝1、2、3)、应力为σk(k＝1、2、3…8)的蠕变实验中蠕变应变εi(i＝1、2、3…5)所对应的蠕变时间。由数据点(0.5％1，12.2118)，结合实际实验过程中蠕变实验温度575℃1，计算综合参数Q118＝10-3×T1×(c+log(t118))式中c为常数参数，此处所取的具体数值为34.34；得到(235MPa8,30118)数据点，求出步骤2)中其他数据点对应的(σk,Qijk)，按蠕变应变值εi由小到大的顺序制成σ-Q关系图。3)提出描述σ和Q关系的数学模型：σ＝ai+b×Q+c×Q2+d×Q3其中ai(i＝1、2、3…n)为与蠕变应变εi相关的常数，此处a1为1250.8；b,c,d为与蠕变应变εi不相关常数，此处b为-34.4；c、d为0。4)将步骤3)中提出的数学模型σ＝ai-34.4×Q拟合于σ和Q的关系图的各个蠕变应变εi下的数据点上，得到图2中的粗黑线，选取蠕变应变0.5％1的拟合线作为基准线，其它拟合线相对于该基准线的偏移量即为Zci即：Zc1＝1250.81-1250.81、Zc2＝1269.32-1250.81、Zc3＝1300.73-1250.81…Zci＝ai-a1。本发明各步骤根据数据的特点选择具体拟合方式，此处拟合选用的软件为origin9.0，也可以选其它相关软件拟合，5)由图2，得到的不同εi所对应的Zci值，以Zc为纵坐标蠕变应变ε为横坐标绘图，并且拟合Zc-ε数据点，得到图3，从图3中可以看出，Zc在蠕变过程中的具体变化，主要分成三个阶段，从图3中选取预测所需的应变1％*1，得到该蠕变应变所对应的18.4*1值；如果所预测的蠕变应变较大，可以延长Zc-ε曲线到所预测的蠕变应变值，再求出相应的Zcm值。6)将由步骤5)中所得到的18.4*1值代入建立的σ-Q关系模型σ＝1250.8+18.4-34.4×Q中，画出σ＝1250.8+18.4-34.4×Q关系线，由该关系线可以可得到设定温度575℃，设定蠕变应力235MPa下1％*1所对应的预测蠕变时间12.8*1,即数据点[1％*1,12.8*1]。7)循环步骤5)和步骤6)得到其它温度，其它蠕变应力下的数据点[ε*m,t*m]。8)由步骤7)中得到的数据点[1％*1,12.8*1]…[ε*m,t*m]，在1％*1和4％*2下，把由该Zc参数方法所得到的t*m与真实蠕变时间texp进行对比得到图4、图5。从图4、图5中可以看出该Zc参数法预测的时间与真实试验时间吻合良好。实施例2本发明中一种基于短时蠕变实验的高温蠕变变形预测方法具体方案由以下实施例详细给出。选取P92钢，开展实验和寿命预测，具体包括以下步骤：1)分别于575℃1下，235MPa8；600℃2下，215MPa6,225MPa7；650℃3下，125MPa1,135MPa2,150MPa3,165MPa4,175MPa5；进行蠕变实验，获的如图6所示得蠕变实验数据。实验过程中温度波动小于1℃；所述的蠕变实验数据的应变范围为0.1％-20％，所述的蠕变实验数据为试样变形值及试样变形值所对应的时间。2)将步骤1)中每个蠕变实验中的蠕变实验数据按蠕变应变从小到大选取不同的蠕变应变，即：0.5％1，1％2，4％3，7％4，10％5蠕变应变，如图6所示；由蠕变应变所对应的时间,获得每个蠕变实验下数据点(0.5％1，12.2118)、(1％2，43.4218)、(4％3，387.5318)、(7％4，495.8418)、(10％5，519518)；(0.5％1，1.6126)、(1％2，9.9226)…(εi，tijk)，其中tijk表示温度为Tj(j＝1、2、3)、应力为σk(k＝1、2、3…8)的蠕变实验中蠕变应变εi(i＝1、2、3、4、5)所对应的蠕变时间。由数据点(0.5％1，12.2118)，结合实际实验过程中蠕变实验温度575℃1，计算综合参数Q118＝logt118+c×T1式中c为常数参数，此处所取的具体数值为0.04；得到(235MPa8,35118)数据点，求出步骤2)中其他数据点对应的(σk,Qijk)，按蠕变应变值εi由小到大的顺序制成σ-Q关系图。3)提出描述σ和Q关系的数学模型：σ＝ai+b×Q+c×Q2+d×Q3其中ai(i＝1、2、3…n)为与蠕变应变εi相关的常数，此处a1为1318.9；b,c,d为与蠕变应变εi不相关常数，此处b为-31.4，c、d为0。4)将步骤3)中提出的数学模型σ＝ai-31.4×Q拟合于σ和Q的关系图的各个蠕变应变εi下的数据点上。得到图7中的粗黑线，选取蠕变应变0.5％1的拟合线作为基准线，其它拟合线相对于该基准线的偏移量即为Zci即：Zc1＝1318.91-1318.91、Zc2＝1337.62-1318.91、Zc3＝1369.43-1318.91…Zci＝ai-1318.91。5)由图7，得到的不同εi所对应的Zci值，以Zc为纵坐标蠕变应变ε为横坐标绘图，并且拟合Zc-ε数据点，得到图8，从图8中可以看出，Zc在蠕变过程中的具体变化，主要分成三个阶段，从图8中选取预测所需的应变1％*1，得到该蠕变应变所对应的18.6*1值；如果所预测的蠕变应变较大，可以延长Zc-ε曲线到所预测的蠕变应变值，再求出相应的Zcm值。6)将由步骤5)中所得到的18.6*1值代入建立的σ-Q关系模型σ＝1318.9+18.6-31.4×Q中，画出σ＝1318.9+18.6-31.4×Q关系线，由该关系线可以可得到设定温度650℃，设定蠕变应力175MPa下1％*1所对应的预测蠕变时间1.24*1,即数据点[1％*1,1.24*1]。7)循环步骤5)和步骤6)得到其它温度，其它蠕变应力下的数据点[ε*m,t*m]。8)由步骤7)中得到的数据点[1％*1,1.24*1]…[ε*m,t*m]，在1％*1和4％*2下，把由该Zc参数方法所得到的t*m与真实蠕变时间texp进行对比得到图9、图10。从图9、图10中可以看出该Zc参数法预测的时间与真实试验时间吻合良好。实施例3本发明中一种基于短时蠕变实验的高温蠕变变形预测方法具体方案由以下实施例详细给出。选取HP耐热钢，开展实验和寿命预测，具体包括以下步骤：1)分别于1000℃1下，32MPa1、35MPa2、37MPa3、38MPa4；980℃2下，42MPa5；950℃3下，50MPa6,900℃4,60MPa7；进行蠕变实验，获的如图11所示的蠕变实验数据。实验过程中温度波动小于1℃；所述的蠕变实验数据的应变范围为0.1％-20％，所述的蠕变实验数据为试样变形值及试样变形值所对应的时间。2)将步骤1)中每个蠕变实验中的蠕变实验数据按蠕变应变从小到大选取不同的蠕变应变，即：0.1％1，0.2％2，0.5％3，1％4，2％5，5％6蠕变应变，如图11所示；由蠕变应变所对应的时间,获得每个蠕变实验下数据点(0.1％1，1.2112)、(0.2％2，7.6212)、(0.5％3，56.2312)、(1％4，89.3412)、(2％5，136.2512)、(5％6，181612)；(0.1％1，1.4113)、(0.2％2，9.8213)、(0.5％3，36.3313)、(1％4，70.9413)、(2％5，104.1513)、(5％6，143.3613)；(0.1％1，1.9125)、(0.2％2，9.5225)、(0.5％3，39.7325)、(1％4，86.4425)、(2％5，132.9525)、(5％6，183.8625)…(εi，tijk)，其中tijk表示温度为Tj(j＝1、2、3…n＇，n＇≥1)、应力为σk(k＝1、2、3…n＂，n＂≥5)的蠕变实验中蠕变应变εi(i＝1、2、3…n)所对应的蠕变时间。由数据点(0.1％1，1.2112)，结合实际实验过程中蠕变实验温度1000℃1，计算综合参数Q112＝10-3×T1×(c+log(t112))式中c为常数参数，此处所取的具体数值为19；得到(35MPa2,24.3112)数据点，求出步骤2)中其他数据点对应的(σk,Qijk)，按蠕变应变εi由小到大的顺序，制成σ-Q关系图。3)提出描述σ和Q关系的数学模型：σ＝ai+b×Q+c×Q2+d×Q3其中ai(i＝1、2、3…n)为与蠕变应变εi相关的常数，此处a1为285.7；b,c,d为与蠕变应变εi不相关常数，此处b为-10.14，c、d为0。4)将步骤3)中提出的数学模型σ＝ai-10.14×Q拟合于σ和Q的关系图的各个蠕变应变εi下的数据点上。得到图12中的粗黑线，选取蠕变应变0.1％1的拟合线作为基准线，其它拟合线相对于该基准线的偏移量即为Zci即：Zc1＝285.71-285.71、Zc2＝295.12-285.71、Zc3＝303.43-285.71…Zci＝ai-285.71。5)由图12，得到的不同εi所对应的Zci值，以Zc为纵坐标蠕变应变ε为横坐标绘图，并且拟合Zc-ε数据点，得到图13，从图13中可以看出，Zc在蠕变过程中的具体变化，主要分成三个阶段，从图13中选取预测所需的应变1％*1，得到该蠕变应变所对应的21.1*1值；如果所预测的蠕变应变较大，可以延长Zc-ε曲线到所预测的蠕变应变值，再求出相应的Zcm值。6)将由步骤5)中所得到的21.1*1值代入建立的σ-Q关系模型σ＝285.7+21.1-10.14×Q中，画出σ＝285.7+21.1-10.14×Q关系线，由该关系线可以可得到设定温度1000℃，设定蠕变应力35MPa下1％*1所对应的预测蠕变时间114.8*1,即数据点[1％*1,114.8*1]。7)循环步骤5)和步骤6)得到其它温度，其它蠕变应力下的数据点[ε*m,t*m]。8)由步骤7)中得到的数据点[1％*1,114.8*1]…[ε*m,t*m]，在1％*1和5％*2下，把由该Zc参数方法所得到的t*m与真实蠕变时间texp进行对比得到图14、图15。从图14、图15中可以看出该Zc参数法预测的时间与真实试验时间吻合良好。实施例4本发明中一种基于短时蠕变实验的高温蠕变变形预测方法具体方案由以下实施例详细给出。选取HP耐热钢，开展实验和寿命预测，具体包括以下步骤：1)分别于1000℃1下，32MPa1、35MPa2、37MPa3、38MPa4；980℃2下，42MPa5；950℃3下，50MPa6,900℃4,60MPa7；进行蠕变实验，获的如图16所示的蠕变实验数据。实验过程中温度波动小于1℃；所述的蠕变实验数据的应变范围为0.1％-20％，所述的蠕变实验数据为试样变形值及试样变形值所对应的时间。2)将步骤1)中每个蠕变实验中的蠕变实验数据按蠕变应变从小到大选取不同的蠕变应变，即：0.1％1，0.2％2，0.5％3，1％4，2％5，5％6蠕变应变，如图16所示；由蠕变应变所对应的时间,获得每个蠕变实验下数据点(0.1％1，1.2112)、(0.2％2，7.6212)、(0.5％3，56.2312)、(1％4，89.3412)、(2％5，136.2512)、(5％6，181612)；(0.1％1，1.4113)、(0.2％2，9.8213)、(0.5％3，36.3313)、(1％4，70.9413)、(2％5，104.1513)、(5％6，143.3613)；(0.1％1，1.9125)、(0.2％2，9.5225)、(0.5％3，39.7325)、(1％4，86.4425)、(2％5，132.9525)、(5％6，183.8625)…(εi，tijk)，其中tijk表示温度为Tj(j＝1、2、3、4)、应力为σk(k＝1、2、3…7)的蠕变实验中蠕变应变εi(i＝1、2、3…6)所对应的蠕变时间。由数据点(0.1％1，1.2112)，结合实际实验过程中蠕变实验温度1000℃1，计算综合参数Q112＝log(t112)+c×T1式中c为常数参数，此处所取的具体数值为0.02；得到(35MPa2,25.5112)数据点，求出步骤2)中其他数据点对应的(σk,Qijk)，按蠕变应变εi由小到大的顺序，制成σ-Q关系图。3)提出描述σ和Q关系的数学模型：σ＝ai+b×Q+c×Q2+d×Q3其中ai(i＝1、2、3…n)为与蠕变应变εi相关的常数，此处a1为320.7；b,c,d为与蠕变应变εi不相关常数，此处b为-11.05，c、d为0。4)将步骤3)中提出的数学模型σ＝ai-11.05×Q拟合于σ和Q的关系图的各个蠕变应变εi下的数据点上。得到图17中的粗黑线，选取蠕变应变0.1％1的拟合线作为基准线，其它拟合线相对于该基准线的偏移量即为Zci即：Zc1＝320.71-320.71、Zc2＝328.92-320.71、Zc3＝336.23-320.71…Zci＝ai-320.71。5)由图17，得到的不同εi所对应的Zci值，以Zc为纵坐标蠕变应变ε为横坐标绘图，并且拟合Zc-ε数据点，得到图18，从图18中可以看出，Zc在蠕变过程中的具体变化，主要分成三个阶段，从图18中选取预测所需的应变1％*1，得到该蠕变应变所对应的18.4*1值；如果所预测的蠕变应变较大，可以延长Zc-ε曲线到所预测的蠕变应变值，再求出相应的Zcm值。6)将由步骤5)中所得到的18.4*1值代入建立的σ-Q关系模型σ＝320.7+18.4-11.05×Q中，画出σ＝320.7+18.4-11.05×Q关系线，由该关系线可以可得到设定温度1000℃，设定蠕变应力35MPa下1％*1所对应的预测蠕变时间115.8*1,即数据点[1％*1,115.8*1]。7)循环步骤5)和步骤6)得到其它温度，其它蠕变应力下的数据点[ε*m,t*m]。8)由步骤7)中得到的数据点[1％*1,115.8*1]…[ε*m,t*m]，在1％*1和5％*2下，把由该Zc参数方法所得到的t*m与真实蠕变时间texp进行对比得到图19、图20。从图19、图20中可以看出该Zc参数法预测的时间与真实试验时间吻合良好。实施例5本发明中一种基于短时蠕变实验的高温蠕变变形预测方法具体方案由以下实施例详细给出。选取DZ125实验材料，开展实验和寿命预测，具体包括以下步骤：1)分别于750℃1下，695MPa1、646MPa2；850℃2下，430MPa3、370MPa4下，348MPa5；950℃3,250MPa6、227MPa7、197MPa8、170MPa9、155MPa10；进行蠕变实验，获的如图21所示的蠕变实验数据。实验过程中温度波动小于1℃；所述的蠕变实验数据的应变范围为0.1％-20％，所述的蠕变实验数据为试样变形值及试样变形值所对应的时间。2)将步骤1)中每个蠕变实验中的蠕变实验数据按蠕变应变从小到大选取不同的蠕变应变，即：0.1％1，0.2％2，0.5％3蠕变应变，如图21所示；由蠕变应变所对应的时间,获得每个蠕变实验下数据点(0.1％1，0.14111)、(0.2％2，0.27211)、(0.5％3，0.71311)；(0.1％1，1.14112)、(0.2％2，2.39212)、(0.5％3，6.32312)…(εi，tijk)，其中tijk表示温度为Tj(j＝1、2、3)、应力为σk(k＝1、2、3…10)的蠕变实验中蠕变应变εi(i＝1、2、3)所对应的蠕变时间。由数据点(0.1％1，0.14111)，结合实际实验过程中蠕变实验温度750℃1，计算综合参数Q111＝log(t111)+c×T1式中c为常数参数，此处所取的具体数值为0.02；得到(695MPa1,19.6111)数据点，求出步骤2)中其他数据点对应的(σk,Qijk)，按蠕变应变εi由小到大的顺序制成σ-Q关系图。3)提出描述σ和Q关系的数学模型：σ＝ai+b×Q+c×Q2+d×Q3其中ai(i＝1、2、3…n)为与蠕变应变εi相关的常数，此处a1为2809.3；b,c,d为与蠕变应变εi不相关常数，此处b为-100.1，c、d为0。4)将步骤3)中提出的数学模型σ＝ai-100.1×Q拟合于σ和Q的关系图的各个蠕变应变εi下的数据点上。得到图22中的粗黑线，选取蠕变应变0.1％1的拟合线作为基准线，其它拟合线相对于该基准线的偏移量即为Zci即：Zc1＝2809.31-2809.31、Zc2＝2839.62-2809.31、Zc3＝2880.43-2809.31…Zci＝ai-2809.31。5)由图22，得到的不同εi所对应的Zci值，以Zc为纵坐标蠕变应变ε为横坐标绘图，并且拟合Zc-ε数据点，得到图23，从图23中可以看出，Zc在蠕变过程中的具体变化，从图23中选取预测所需的应变0.1％*1，得到该蠕变应变所对应的0*1值；如果所预测的蠕变应变较大，可以延长Zc-ε曲线到所预测的蠕变应变值，再求出相应的Zcm值。6)将由步骤5)中所得到的0*1值代入建立的σ-Q关系模型σ＝2809.3+0-100.1×Q中，画出σ＝2809.3+0-100.1×Q关系线，由该关系线可以可得到设定温度750℃，设定蠕变应力695MPa下0.1％*1所对应的预测蠕变时间*1,即数据点[0.1％*1,4.56*1]。7)循环步骤5)和步骤6)得到其它温度，其它蠕变应力下的数据点[ε*m,t*m]。8)由步骤7)中得到的数据点[0.1％*1,4.56*1]…[ε*m,t*m]，在0.1％*1和0.5％*2下，把由该Zc参数方法所得到的t*m与真实蠕变时间texp进行对比得到图24、图25。从图24、图25中可以看出该Zc参数法预测的时间与真实试验时间吻合良好。以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。