下载文档原格式
inconel600材料norton蠕变本构系数Inconel 600是一种镍铬合金，具有很好的耐高温、耐腐蚀性能，被广泛应用于航空航天、化学工业、核工业等领域。
在工程设计和分析中，对材料的蠕变特性进行准确描述是十分重要的，因为蠕变现象对结构的稳定性和寿命有很大影响。
蠕变是指材料在高温和常载荷下逐渐发生变形的现象。
在材料中，蠕变是由于晶体的滑移和再结晶而导致的。
蠕变行为受到一些因素的影响，其中一个重要的因素是材料的蠕变本构系数。
本构系数是描述材料变形规律的参数，通过建立合适的蠕变本构模型，可以预测和分析材料在不同工作条件下的蠕变行为。
Norton模型是描述金属材料蠕变行为的一种经验模型，以James Norton的名字命名。
该模型是最早的一种用于描述金属蠕变行为的本构模型，基于实验数据曲线进行拟合得到。
Norton模型的表达式为：ε = Aσ^nexp(Q/RT)其中，ε是蠕变应变，σ是应力，A是材料的常数，n是蠕变指数，Q是材料的激活能，R是普适气体常数，T是绝对温度。
对于Inconel 600材料而言，其蠕变本构系数可以通过实验数据进行拟合得到。
在实际应用中，通常需要通过试验来确定Inconel 600材料的蠕变指数、激活能等参数。
通过测量不同温度和应力下的蠕变应变，可以得到一组应变-应力数据点，然后利用最小二乘法等拟合方法来确定Norton模型中的参数。
Inconel 600材料的蠕变行为对于不同应力和温度下的应变和时间的关系进行了广泛的研究。
通过实验，可以得到Inconel 600材料在不同温度和应力下的蠕变曲线。
这些曲线可以用来拟合Norton模型，从而确定Inconel 600材料的蠕变本构系数。
值得注意的是，由于Inconel 600材料在高温下有良好的耐蠕变性能，因此在实际应用中通常需要进行较长时间的蠕变实验，以获得准确的蠕变本构系数。
除了Norton模型外，还有其他一些用于描述金属材料蠕变行为的本构模型，如NP模型、DD Visco Model等。
岩石材料的蠕变实验及本构模型研究蠕变是指材料在一定温度和应力条件下，随着时间的推移发生的持续变形。
在地质和工程领域，岩石是一种典型的蠕变材料。
岩石的蠕变行为对工程结构的长期稳定性和可靠性具有重要影响。
因此，对岩石材料的蠕变实验及本构模型研究具有重要的理论和实际意义。
岩石材料的蠕变实验主要分为应力松弛实验和恒定应力蠕变实验两种。
应力松弛实验是通过对材料施加一定的应力后，观察材料的应力随时间的变化，以及应变随时间的变化。
这种实验常常用来研究岩石材料的蠕变速率和蠕变变形的领导指数。
恒定应力蠕变实验则是在一定的应力水平下，观察材料的应变随时间的变化，并且通过实验数据拟合来得到本构模型。
岩石材料的蠕变行为可以通过多种本构模型来描述，其中最常用的是Norton、Burgers、Power-law以及Generalized Kelvin-Voigt模型。
这些模型可以通过实验数据进行参数拟合，从而得到对应的本构关系。
这些本构关系可以用来预测岩石材料在不同应力和温度下的蠕变行为。
此外，还可以通过拟合这些本构模型的参数，来研究岩石材料的蠕变机制。
研究表明，岩石材料的蠕变行为是由多种因素共同影响的，包括温度、应力水平、孔隙水压力、孔隙率等。
因此，在进行蠕变实验时，需要对这些因素进行控制和监测，以保证实验数据的可靠性。
同时，还需要考虑到实际工程环境中的应力和温度条件，从而得到更准确的本构关系。
总之，岩石材料的蠕变实验及本构模型研究对于预测岩石在地下工程中的蠕变行为具有重要的理论和实际意义。
通过研究岩石材料的蠕变行为及其本构关系，可以为地质和工程领域提供重要的科学依据，从而保证工程结构的长期稳定性和可靠性。
岩石材料的蠕变实验及本构模型研究引言：岩石是地球上最基础的构造材料之一，其性质的研究对于地质科学以及岩土工程领域具有重要意义。
岩石在地壳中扮演着起支撑与保护作用，因此了解岩石的变形行为以及蠕变性质对于地质灾害的预测与评估具有重要的指导意义。
本文将就岩石材料的蠕变实验及本构模型研究进行详细阐述。
一、岩石材料的蠕变实验蠕变是指物质在长时间内受到持续应力下的变形现象。
岩石材料由于具有多种类型的孔隙和裂隙，因此其蠕变行为比一般材料更为复杂。
蠕变实验是研究岩石材料蠕变性质的主要手段之一，其目的是了解岩石在不同应力、不同温度和不同时间下的蠕变特性。
1.实验设备蠕变实验一般需要使用蠕变试验机，该仪器能够提供连续加载并测量样品的应力和应变，同时控制温度。
实验所需的试样通常需要根据具体需要制备。
此外，还需要一些测量设备，如蠕变计和应变测量仪等。
2.实验过程蠕变实验的过程包括准备试样、加载试样、施加应力、保持应力和测量应变等步骤。
首先，需要根据实验要求制备符合标准的试样。
然后，将试样放置在蠕变试验机上，施加适当的负载并开始加载。
在加载过程中，需要保持恒定的应力并测量试样的应变，常用的应变测量方法有外部应变计和内部传感器等。
最后，根据实验结果绘制蠕变曲线，分析蠕变行为。
本构模型是描述材料力学性质的数学模型，通过建立岩石材料的本构模型，可以预测岩石的变形行为并进行力学仿真研究。
目前常用的岩石本构模型有线性弹性模型、弹塑性模型和粘弹性模型等。
1.线性弹性模型线性弹性模型是最简单的本构模型，它假设岩石材料的应力应变关系是线性的，即满足胡克定律。
这种模型适用于小应变范围内的岩石变形，但无法描述岩石的时间依赖性和非线性特性。
2.弹塑性模型弹塑性模型考虑了岩石在加载时的弹性变形和塑性变形，常用的模型有Mohr-Coulomb模型、Drucker-Prager模型等。
这些模型能够更准确地描述岩石的变形行为，但在蠕变时间很长的情况下，塑性本构模型可能会失效。
岩石材料的蠕变实验及本构模型研究小伙伴们！今天咱们来捣鼓一下岩石材料的蠕变实验以及本构模型研究这个超有趣的事儿。
首先呢，咱们得知道啥是岩石材料的蠕变。
简单来说，就好比你长时间弯腰干活，腰会慢慢变得更酸更累，岩石在受到持续的力的时候，它也会慢慢变形，这个过程就是蠕变。
那我们为啥要研究这个呢？我觉得吧，这就像是医生要研究病人的慢性病一样。
在很多工程里，像修隧道啊，建大坝如果岩石偷偷地在蠕变，那可就像一颗定时炸弹，说不定啥时候就出大问题了。
那咱们开始做这个实验吧。
你得先找一些合适的岩石样本。
这可不像在路边随便捡块石头就行，得是那种有代表性的，可能从矿山里采出来的就比较靠谱。
我曾经就跟着老师去矿山找样本，那地方啊，又脏又乱，我还差点被一块小石头绊倒，真是“出师未捷身先死”啊，不过好在最后找到了合适的岩石。
把岩石样本带回来后，就把它放到专门的实验设备里。
这个设备就像一个超级大力士，能够持续地给岩石施加压力。
然后呢，我们就在旁边眼巴巴地看着，就像等着蛋糕出炉一样，盼着岩石一点点变形。
这个过程可能得持续好长好长时间，也许几天，也许几周，甚至几个月。
这时候我就想啊，这岩石可真有耐心，要是我被这么压着，早就受不了啦。
在这个过程中，我们得不停地记录数据。
就像侦探在记录犯罪现场的每一个细节一样，一点都不能马虎。
比如说岩石在某个时刻变形了多少毫米啊，这个数据可金贵着呢。
然后就是本构模型的研究啦。
这玩意儿就像是给岩石的蠕变行为画一幅画像。
我们要根据收集到的数据，建立一个数学模型，来描述岩石在蠕变过程中的各种规律。
这就有点像猜谜语，我们知道了一些线索，然后要猜出整个谜底。
不过这个谜底可不好猜，也许我们建立的模型一开始漏洞百出，就像一件破了好多洞的衣服。
但是没关系不断地修改，不断地完善，就像补衣服一样，总有一天能让这个模型完美地描述岩石的蠕变。
那这个本构模型有啥用呢？我觉得它就像是一个导航仪。
在工程建设中，我们可以根据这个模型来预测岩石的蠕变，提前做好防范措施。
Hans Journal of Civil Engineering 土木工程, 2023, 12(9), 1197-1202 Published Online September 2023 in Hans. https:///journal/hjce https:///10.12677/hjce.2023.129138水环境下岩石蠕变特性及本构模型研究综述石棋军重庆科技学院建筑工程学院，重庆收稿日期：2023年8月6日；录用日期：2023年8月27日；发布日期：2023年9月7日摘要 岩石的重要力学特性之一是蠕变，而水又是影响岩石的蠕变力学性能的一个关键因素。
本文主要阐述了近年来国内外学者对岩石在水环境下进行的蠕变试验和本构方程方面的研究成果，并对未来的岩石蠕变及本构方程研究方向进行了展望。
关键词岩石，蠕变试验，本构方程，水环境Research Review on Creep Characteristics and Constitutive Model of Rock under Water EnvironmentQijun ShiSchool of Civil Engineering and Architecture, Chongqing University of Science and Technology, Chongqing Received: Aug. 6th , 2023; accepted: Aug. 27th , 2023; published: Sep. 7th , 2023AbstractOne of the important mechanical properties of rocks is creep, and water is a key factor affecting the creep mechanical properties of rocks. In this paper, the research results of creep test and con-stitutive equation of rock in water environment carried out by domestic and foreign scholars in recent years are mainly described, and the research direction of rock creep and constitutive equa-tion in the future is prospected.石棋军KeywordsRock, Creep Test, Constitutive Equation, Water EnvironmentCopyright © 2023 by author(s) and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License (CC BY 4.0)./licenses/by/4.0/1. 引言随着我国经济实力的快速发展，如地下洞室、大型水电站、道路桥梁等基础工程开始大规模建设，在建设这些大工程中，水环境中岩土工程类问题急需解决，水对岩石的作用对岩土类工程的时效变形和稳定性有着举足轻重的地位，因此这类问题成为国内外学者的关注重点。
高模量沥青混凝土蠕变本构方程的有限元分析第一章：引言随着路面交通的不断发展，人们对于道路的要求也越来越高。
此时，高模量沥青混凝土作为一种新型的路面材料已经出现并受到了广泛的应用。
高模量沥青混凝土由于其优秀的性能，被广泛应用于公路、机场跑道、码头等工程建设中。
但是高模量沥青混凝土的蠕变现象对于该材料的使用和设计是具有一定的影响的。
因此，本文提出了高模量沥青混凝土蠕变本构方程的有限元分析。
第二章：高模量沥青混凝土及其蠕变行为的研究本章介绍了高模量沥青混凝土的性质和特点。
其中，重点分析了高模量沥青混凝土在蠕变方面的表现和特性，探讨了蠕变行为和影响其蠕变行为的因素。
第三章：蠕变本构方程的基本原理本章主要介绍了蠕变本构方程的发展和基本原理。
在介绍相关理论的同时，本章还特别关注了蠕变曲线和蠕变本构方程的关系以及常用的蠕变本构方程型式。
第四章：高模量沥青混凝土蠕变本构方程的有限元分析本章主要介绍了基于有限元的高模量沥青混凝土蠕变本构方程分析方法。
该方法采用有限元模拟软件ANSYS APDL进行模拟实验，对高模量沥青混凝土蠕变本构方程进行分析与计算。
本章详细介绍了高模量沥青混凝土的有限元分析流程，从材料属性的定义到计算结果的解读都进行了一定程度的阐述。
第五章：结论与展望基于第四章的模拟实验和分析结果，本章对高模量沥青混凝土的蠕变规律和本构方程进行了总结和归纳。
并通过对其优缺点的分析和评估，着重探讨了高模量沥青混凝土在今后工程建设中的应用前景及其未来的研究方向。
总结本文基于高模量沥青混凝土的特性和性质，深入探讨了高模量沥青混凝土在蠕变方面的行为和表现。
在此基础上，提出了高模量沥青混凝土蠕变本构方程的有限元分析方法，并通过实验验证了模拟结果的准确性和可行性。
本文的研究结果为今后的道路建设和工程设计提供了一定的理论依据和技术支持。
ansys蠕变本构模型蠕变是指物体在一定温度和应力条件下逐渐发生形变的现象。
蠕变本构模型是用来描述材料的蠕变性能的数学模型。
在工程领域中，蠕变的研究对于可靠性和寿命预测等方面具有重要意义。
在ANSYS软件中，有多种蠕变本构模型可供选择，如Norton本构模型、Manson-Haftka本构模型、Power Law本构模型等。
Norton本构模型是最常用、最简单的蠕变本构模型之一。
它基于实验事实，即材料蠕变速率与应力的指数幂关系。
Norton本构模型可以用下面的公式表示：εc = Aσ^n exp(-Q/RT)其中，εc是蠕变应变，A是一个与材料的本构参数有关的常数，σ是应力，n是指数，Q是激活能，R是气体常数，T是材料的温度。
这个模型适用于大多数金属和合金的蠕变行为。
Manson-Haftka本构模型是一种更复杂的蠕变本构模型。
它考虑了应力和温度的交互作用，并使用一个修正系数来描述蠕变应变的非线性行为。
Manson-Haftka本构模型可以用下面的公式表示：εc = Aσ^n exp(-Q/RT) [1 + β(εP)^m]其中，εc是蠕变应变，A、n、Q、R、T的含义与Norton本构模型相同，β是一个非线性修正系数，εP是塑性应变，m是一个与材料有关的常数。
这个模型适用于不同应力水平下的蠕变行为。
Power Law本构模型是基于实验事实，即材料蠕变速率与应力的幂函数关系。
Power Law本构模型可以用下面的公式表示：εc = Aσ^n其中，εc是蠕变应变，A是一个与材料的本构参数有关的常数，σ是应力，n是指数。
这个模型适用于高负荷条件下的蠕变行为。
以上是ANSYS中常用的几种蠕变本构模型。
根据实际情况和材料性质的不同，可以选择适合的蠕变本构模型来模拟蠕变行为。
这些模型可以帮助工程师更好地理解和预测材料的蠕变性能，从而优化设计和提高产品的可靠性。
同时，在ANSYS软件中，还可以根据实验数据对蠕变本构模型参数进行拟合和优化，从而更准确地描述材料的蠕变行为。
蠕变是指材料在恒定应力或荷载下，随时间的推移而发生形变。
在有限元模拟计算中，利用蠕变本构方程可以计算材料在热处理过程中的应力松弛行为。
蠕变本构方程描述了材料的应力-应变关系，常用的蠕变本构方程有Norton、Dorn和Sherby-Dorn等形式。
在有限元模拟计算中，针对不同材料的蠕变行为，需要选择适当的蠕变本构方程。
这些方程以实验数据为基础，通过曲线拟合等方法得到。
下面是一些常用的参考内容，用于模拟材料在热处理过程中的应力松弛行为。
1.Norton本构方程： Norton本构方程是最早应用于蠕变研究的方程之一，它假设材料的蠕变变形由渐进微裂纹的扩展引起。
Norton方程可以写成如下形式：ε ̇=ε ̇s(σ) 其中，ε ̇为应变速率，ε ̇s为材料的引导应变速率，σ为应力。
Norton方程常用于描述金属和陶瓷等材料的蠕变行为。
2.Dorn本构方程： Dorn本构方程常用于描述塑料材料的蠕变行为，特别是在中高温下。
Dorn方程可以写成如下形式：σ=σ0(1+ε ̇/˙ε0)^n 其中，σ为应力，σ0为参考应力，ε ̇为应变速率，˙ε0为参考应变速率，n为指数。
Dorn方程通过参数σ0和n来描述材料的蠕变行为。
3.Sherby-Dorn本构方程： Sherby-Dorn本构方程适用于金属和合金的高温蠕变行为。
Sherby-Dorn方程可以写成如下形式：ε ̇=Aσ^mexp(-Q/RT) 其中，ε ̇为应变速率，σ为应力，A为常数，m为指数，Q为激活能，R为气体常数，T为温度。
Sherby-Dorn方程通过参数A、m、Q来描述材料的蠕变行为。
除了选择合适的蠕变本构方程外，有限元模拟计算中还需要考虑材料的初始条件、边界条件以及模拟的时间步长等因素。
此外，要进行有限元模拟计算还需要确定材料的杨氏模量、泊松比等弹性性质，这些参数可以通过实验测量或者经验公式来估计。
需要注意的是，蠕变行为是一个复杂的物理过程，受多种因素影响，如温度、应力、时间等。
plaxis 蠕变本构模型
Plaxis是一个常用的岩土工程有限元分析软件，用于分析地下结构的行为。
在Plaxis中，蠕变是指土体在长期荷载作用下发生的变形现象，通常需要使用本构模型来描述土体的蠕变行为。
在Plaxis中，常用的蠕变本构模型包括线性本构模型、多项式本构模型、本构模型和黏塑性本构模型等。
这些本构模型可以用来描述土体在长期荷载作用下的变形特性，包括蠕变变形、应力松弛等现象。
线性本构模型假设土体的蠕变行为是线性的，可以通过弹性模量和蠕变模量来描述土体的蠕变特性。
多项式本构模型则通过多项式函数来描述土体的蠕变行为，可以更准确地反映土体的实际蠕变特性。
本构模型和黏塑性本构模型则更加复杂，可以考虑土体的非线性蠕变行为和黏塑性特性，更适用于描述复杂的土体行为。
这些本构模型在Plaxis中都有相应的参数，用户可以根据实际情况选择合适的本构模型来模拟土体的蠕变行为。
总的来说，在Plaxis中，蠕变本构模型是用来描述土体在长期荷载作用下的变形特性的重要模型，不同的本构模型适用于不同类型的土体和工程条件，工程师需要根据实际情况选择合适的本构模型来进行分析和设计。
一种非定常参数的岩石蠕变本构模型罗润林;阮怀宁;孙运强;朱昌星【摘要】在前人的研究基础上提出了一种新的非定常参数的岩石蠕变模型,该模型由一个牛顿体和一个类似于开关的SO元件并联后与一个广义Kelvin体串联,通过SO元件的作用,新模型可以在Burgers模型与广义Kelvin模型之间切换,并将Burgers模型中牛顿体的粘滞系数η1看成与时间有关的非定常参数,使得模型不仅可以反映稳态蠕变过程也可以反映加速蠕变过程.【期刊名称】《桂林理工大学学报》【年(卷),期】2007(027)002【总页数】4页(P200-203)【关键词】非定常参数;蠕变;本构模型;Burgers模型;Kelvin 模型【作者】罗润林;阮怀宁;孙运强;朱昌星【作者单位】河海大学,岩土工程研究所,南京,210098;河海大学,岩土工程研究所,南京,210098;新疆石油管理测井公司研究所,新疆,克拉玛依,834000;河海大学,岩土工程研究所,南京,210098【正文语种】中文【中图分类】TU452大量的试验数据和工程实践表明，岩石的流变破坏普遍存在，许多岩体工程的破坏都因岩石流变破坏而引起，并在破坏前都有加速流变过程.根据目前的研究结果，描述岩石蠕变的本构模型主要有经验模型、元件组合模型、内时理论模型和屈服面模型等［1］. 元件模型是用模型元件线性组合来模拟岩土的流变行为，可以描述岩石的定常蠕变和衰减蠕变［2］. 但是，其力学性质单一，通过参数调整仍然无法定量模拟实测的应力－应变－时间曲线，特别是无法表现岩石的破坏特性，也即不能描述岩石的加速蠕变. 于是，研究人员将目光投向非线性理论，采用非线性元件( 如非线性弹性元件、非线性粘性元件等) 代替传统模型中部分或者全部线性元件. 金丰年等［3］基于试验结果，结合传统线性粘弹性模型的分析，提出了非线性粘弹性模型. 邓荣贵等［4］引进一种非线性粘滞阻尼器，该阻尼器所受应力与其蠕变加速度大小成正比，将该元件与其他元件组合形成一种新的复合元件，并在工程中得到了应用. 韦立德［5］等根据粘聚力在流变中的作用，建议采用一个类似于开关的SO 元件，并由此建立了一维粘弹塑性本构模型. G. N. Boukharov［6］于1995 年提出了一种具有一定质量的延迟阻尼器元件，该元件有一应变门槛值，当应变大于该值时，模型发生加速运动. 徐卫亚［7］提出了一个考虑粘聚力与内摩擦系数的岩石非线性粘塑性体，并与五元件粘弹性模型串联建立非线性粘弹塑性流变模型. 陈沅江［8］等引入一种变截面积的阻尼器来模拟岩石的加速蠕变. 曹树刚等［9］采用非牛顿体粘性元件构成五元件的改进西原正夫模型，探讨了与时间有关的软岩一维和三维本构方程和蠕变方程. 张向东［10］等基于泥岩的三轴蠕变试验结果，建立了泥岩的非线性蠕变方程，并以此分析了围岩的应力场和位移场. 王来贵等［11］以改进的西原正夫模型为基础，利用岩石全程应力-应变曲线与蠕变方程中参数的对应关系，建立了参数非线性蠕变模型.许宏发［12］通过对等时应力－应变曲线分析认为，弹模与时间符合指数的关系. 受此启发，丁志坤［13］提出了一种非定常参数的H－K 模型，与定常参数的模型对比发现，非定常模型能够更为准确地反映岩石的粘弹性变形性能. 本文为了能够同时描述岩石的稳态蠕变和加速蠕变过程，提出一种新的非定常参数岩石蠕变模型.1 岩石的定常参数蠕变模型大量的试验显示，岩石在不同应力下的蠕变特性可以归纳为图1 的3 种曲线. 曲线I 和曲线II为应力较低时的岩石的稳态蠕变曲线，曲线I 在初始蠕变后进入速率为零的稳态蠕变，曲线II 在初始蠕变后进入速率为一常数的稳态蠕变. 曲线III具有初始蠕变、稳态蠕变和加速蠕变3 个阶段.对于曲线I 和曲线II，已有传统的元件组合模型可以描述，如Burgers 模型可以描述曲线II 的蠕变过程.Burgers 模型是一种粘弹性蠕变模型，由一个广义Kelvin 体和一个牛顿体串联而成 ( 图2) .Burgers 模型的蠕变方程可以表示为［14］：从图3 的Burgers 模型的蠕变曲线可以看出，该模型描述图1 中的曲线Ⅱ时，稳态蠕变速率趋于σ0/η1，也就是说，稳态蠕变速率的大小取决于牛顿体的粘滞系数η1. 由于该模型有串联牛顿体元件，所以不能反映I 型曲线. 同样，由于牛顿体的粘滞系数η1是一常数，从而限制了Burgers 模型对加速蠕变过程描述.图1 3 种典型的岩石蠕变曲线Fig.1 3 kinds of typical creep curves图2 Burgers 模型Fig.2 Burgers model图3 Bugers 模型的蠕变曲线Fig.3 Creep curve of Burgers model2 非定常参数的Burgers 蠕变模型前人的研究结果表明，岩石在不同荷载下都存在图1 中的3 种曲线，当荷载较小时蠕变过程如曲线I 所示，当荷载大于某一门槛值σs 时，岩石的蠕变过程可以由图1 中的曲线II 来描述. 如上所分析，由于串联有牛顿体元件，传统的Burgers模型不能描述图1 中稳态速率为0 的蠕变曲线I.然而广义Kelvin 可以很好地描述稳态蠕变过程［14］.从图2 及式( 1) 可以看出Burgers 模型是由广义Kelvin 模型和一个牛顿体串联而成. 因此，可以采用一个开关元件与牛顿体并联：当荷载σ ＞σs时，开关开，牛顿体起作用，模型变成Burgers 模型，可以描述II 型曲线;σ≤σs 时开关关闭，牛顿体不起作用，模型变成广义Kelvin，可以描述I型曲线. SO 元件就是这样一个开关元件［5］，元件图形如图4 所示.图4 SO 元件Fig.4 SO component其应力应变关系：当σ≤σs 时，元件的应变为0; 当σ ＞σs 时，元件的应变为任意值. 该元件与圣维南体有点类似，所不同的是，圣维南体的σ's 是达到塑性应变时的应力，并且分担大小为σ's 的荷载. 而SO 元件的σs 则是一个门槛值，反映蠕变进入II 型曲线时的应力值，可能会小于σ's，且不分担荷载.将SO 模型与Burgers 模型组合形成的新模型如图5 所示. 该模型的数学表示为图5 本文提出的非定常参数的蠕变模型Fig.5 Non-stationary parameters creep model当σ ≤σs 时，新模型为广义Kelvin 体当σ ＞σs 时，模型为式(1) 的Burgers 体.这样，新模型可以反映岩石蠕变过程的I 型曲线和II 型曲线.丁志坤等［13］为了更好反映岩石的粘弹性变形性能，将广义Kelvin 体的E2参数采用与时间有关的非定常参数E2( t) 来表示，并取得了很好的效果.本文中E2( t) 仍然采用此表述，即从上一节分析知道，由于Burgers 模型中的η1为一常数而不能描述岩石蠕变的加速过程( 图1中曲线III) ，如果将η1 也采用与时间有关的非定常参数η1( t) 即可以解决这一问题.由于粘滞系数随时间不断减小，初始减小不明显，而当应变大于某一值时，粘滞系数突然加速减小，最终导致粘滞力为0 使得岩石破坏［6，8，13］，楼志文［15］、金丰年［16］指出岩石的损伤变量D 与时间t呈指数型关系. 由此，可以假设η1( t) 符合指数衰减的规律，即由式(2) 、(5) 和(6) 可以得到图5 模型的表达式H( σ－σs) 的含义与式(3) 一致.3 模型应用研究图6 中的虚线是某泥质板岩的单轴压缩蠕变试验结果［12］. 在不同荷载下，岩石出现不同的蠕变过程，在该试验中出现了图1 的3 种类型曲线.采用粒子群算法［17］，利用本文提出的模型( 式(7) ) 对该试验数据进行参数识别. 参数识别时，根据先验知识给出模型的范围，在此范围随机生成20 个粒子形成粒子群，当全局最优粒子的拟合误差小于1%或者迭代40 次，辨识完成. 如果没有先验知识，可以通过试算的方法进行辨识，即先给出模型的范围，如果某个模型参数与其边界接近，那么可以放宽模型的边界( 如果与下界接近，减小下界; 如果跟上界接近，则增大上界)再进行识别. 从拟合结果( 图6) 看出，本文提出的非定常参数蠕变模型与试验结果非常吻合，表明该蠕变模型不仅可以反映岩石的稳态蠕变过程，也能很好地反映岩石的加速蠕变过程. 从拟合误差( 表1) 来看，不同荷载下的拟合误差都很小，最大的拟合误差只有1.22%，最小拟合误差为0.56%.图6 蠕变试验曲线和理论曲线Fig.6 Curves of experiment vs. theory表1 参数识别结果Table 1 Identified results of parameters荷载/MPa E1/Gpa p1/GPa p2/GPa p3/1·d－1 η2/( GPa·d－1)η10/( GPa·d－1) a/10－5 b/( d－1)拟合误差/%5.08 4.46 2.24 6.55 0.627 58.60 0.66 7.06 4.82 2.14 8.29 0.617 61.49 1.20 10.10 4.82 2.11 8.24 0.604 61.19 1.07 13.16 4.19 2.53 7.56 0.17158.54 2 007.52 1.79 1.268 0.56 14.29 4.05 2.36 7.27 0.158 57.57 2 302.840.015 9 1.66 1.16 15.82 4.39 1.85 6.42 0.118 54.74 2 189.94 0.010 3 1.7551.22从拟合所得的参数( 表1) 看出，E1、p1、p2 和η10 与荷载大小基本无关，而与岩石的特性及试件的状态有关.在σ ≤σs 时，η2 也与荷载大小无关，这与文献［7］的结果一致，p3 随着荷载的增加而减小，也即随时间增大，E2( t) 的减小速度减慢，说明在荷载较小时，岩石蠕变初始时间较短，容易过渡到稳态蠕变( 图6) ，荷载越大，由初始蠕变过渡到稳态蠕变的时间越长，即越不容易达到稳态流变，说明该模型也能很好地反映岩石由初始蠕变进入稳态蠕变的过程.当σ ＞σs 时，p3和η2 都随着荷载的增大而减小，p3 的减小尤为明显.b 随荷载的增大而增大，由式(3) 和式(6) 看出，b 是主要反映加速蠕变的系数，b 越大，加速蠕变过程的蠕变加速度就越大.而a 则随荷载的增大而减小.4 结论(1) 通过类似于开关的SO 元件与牛顿体并联后与广义Kelvin 体串联形成一个新的岩石蠕变模型，该模型通过SO 元件的作用，可以实现Burgers模型与广义Kelvin 模型之间的切换，当σ ＞σs 时，新模型成为Burgers 模型，以描述速率不为0 的稳态蠕变. σ≤σs 时新模型转化为广义Kelvin 模型，以反映速率为0 的稳态蠕变.(2) 将牛顿体的粘滞系数假设为与时间有关并满足式( 7) 的关系式，与SO 并联后与文献［13］的非定常参数的广义Kelvin 体串联形成新的非定常参数的蠕变模型，根据以上分析，该模型可以很好地反映岩石的3 种蠕变过程.(3) 经过分析，参数p3、a、b 对荷载反映较为敏感.η2 只有在σ ＞σs 时随荷载的增加而减小.而其余参数受荷载的影响很小. 分析结果还表明，本文提出的新模型能很好地反映不同荷载下岩石由初始蠕变进入稳态蠕变的过程.【相关文献】［1］袁静，龚晓南，益德清. 岩土流变模型的比较研究［J］. 岩石力学与工程学报，2001，20(6) ：772－ 779.［2］夏才初，孙钧.蠕变试验中流变模型辨识及参数确定［J］. 同济大学学报，1996，24(5) ：498－503.［3］金丰年，浦奎英. 关于粘弹性模型的讨论［J］. 岩石力学与工程学报，1995，14(4) ：335－361.［4］邓荣贵，周德培，张倬元，等.一种新的岩石流变模型［J］.岩石力学与工程学报，2001，20(6) ： 780－784.［5］韦立德.岩石力学损伤和流变本构模型研究［D］.南京：河海大学，2003.［6］Boukharov G N，Chanda M W. The three processes of brittle crystalline rock creep ［J］. Int. J. Rock Mech. Min. Sci. ＆Geomech.Abstr.，1995，32(4) ：325－335.［7］徐卫亚，杨圣奇，褚卫江.岩石非线性粘弹塑性流变模型( 河海模型) 及其应用［J］. 岩石力学与工程学报，2006，25(3) ：433-447.［8］陈沅江，潘长良，曹平，等. 软岩流变的一种新力学模型［J］.岩土力学，2003，24(2) ：209－214.［9］曹树刚，边金，李鹏. 岩石蠕变本构关系及改进的西原正夫模型［J］. 岩石力学与工程学报，2002，21(5) ：632-634.［10］张向东，李永靖，张树光，等. 软岩蠕变理论及其工程应用［J］. 岩石力学与工程学报，2004，23(10) ：1635-1639.［11］王来贵，何峰，刘向峰，等. 岩石试件非线性蠕变模型及其稳定性分析［J］. 岩石力学与工程学报，2004，23(10) ：1640-1642.［12］许宏发.软岩的强度和弹模时间效应研究［J］.岩石力学与工程学报，1997，16(3) ：246－251.［13］丁志坤，吕爱钟. 岩石粘弹性非定常蠕变方程的参数识别［J］.岩土力学，2004，25( S) ：37－40.［14］蔡美峰，何满潮，刘东燕. 岩石力学与工程［M］. 北京：科学出版社，2002.［15］楼志文. 损伤力学基础［M］. 西安：西安交通大学出版社，1990.［16］金丰年，范华林，浦奎源. 岩石蠕变损伤模型研究［J］. 工程力学，2001( S) ：227－231. ［17］Kennedy J，Eberhart R C. Particle swarm optimization［C］//IEEE International Conference on Neural Networks，Perth，1995：1942－1948.。
第15卷第7期李俞韦，等：2024铝合金孔洞缺陷搅拌摩擦点焊修复数值模拟与实验研究95Finite-Element Studies on the Effect of Tool Shape inFriction Stir Welding[J]. Journal of Engineering Manu-facture, 2010, 224(8): 1161-1173.[21] BENSON D J, OKAZAWA S. Contact in a Multi-Mater-ial Eulerian Finite Element Formulation[J]. ComputerMethods in Applied Mechanics and Engineering, 2004,193(39/40/41): 4277-4298.[22] LIU Qi-peng, LI Wen, ZHU Lei, et al. Tempera-ture-Dependent Friction Coefficient and Its Effect onModeling Friction Stir Welding for Aluminum Alloys[J].Journal of Manufacturing Processes, 2022, 84: 1054-1063.[23] ZAHMATKESH B, ENAYATI M H, KARIMZADEH F.Tribological and Microstructural Evaluation of FrictionStir Processed Al2024 Alloy[J]. Materials & Design,2010, 31(10): 4891-4896.[24] AL-BADOUR F, MERAH N, SHUAIB A, et al.Thermo-Mechanical Finite Element Model of FrictionStir Welding of Dissimilar Alloys[J]. The InternationalJournal of Advanced Manufacturing Technology, 2014,72(5): 607-617.[25] MISHRA R S, MA Z Y. Friction Stir Welding and Proc-essing[J]. Materials Science and Engineering: R: Re-ports, 2005, 50(1/2): 1-78. [26] 赵华夏, 董继红, 孟强. 7050高强铝合金搅拌摩擦焊典型宏观缺陷试验研究[J]. 精密成形工程, 2019, 11(6): 141-148.ZHAO Hua-xia, DONG Ji-hong, MENG Qiang. Ex-perimental Study on Typical Macroscopic Defects of 7050 High Strength Aluminum Alloy Friction Stir Welding[J]. Journal of Netshape Forming Engineering, 2019, 11(6): 141-148.[27] 肖旋, 秦鼎强, 倪昱, 等. 铝合金薄板搭接高速FSW缺陷及断裂行为[J]. 精密成形工程, 2019, 11(6): 135- 140.XIAO Xuan, QIN Ding-qiang, NI Yu, et al. Defects ofHigh Speed Friction Stir Welding and Fracture Behaviorof Aluminum Alloy Thin Plate Lap Joints[J]. Journal ofNetshape Forming Engineering, 2019, 11(6): 135-140. [28] 邓运来, 邓舒浩, 叶凌英, 等. 焊后热处理对AA7204-T4铝合金搅拌摩擦焊接头组织与力学性能的影响[J].材料工程, 2020, 48(4): 131-138.DENG Yun-lai, DENG Shu-hao, YE Ling-ying, et al.Effects of Post-Weld Heat Treatment on Microstructuresand Mechanical Properties of AA7204-T4 Aluminum Alloy FSW Joint[J]. Journal of Materials Engineering, 2020, 48(4): 131-138.责任编辑：蒋红晨精 密 成 形 工 程第15卷 第7期96 JOURNAL OF NETSHAPE FORMING ENGINEERING2023年7月收稿日期：2023–03–09 Received ：2023-03-09基金项目：航空科学基金（2020Z047056003）；江西省重点研发计划（20202BBEL53012） Fund ：Aeronautical Science Foundation of China(2020Z047056003); Key Research and Development Project of Jiangxi Province (20202BBEL53012)作者简介：徐显强（1996—），男，硕士生。
所以在有限元模拟计算中，利用蠕变本构方程计算热处理过程中的应力松弛行为，蠕变是材料在恒温条件下经受恒定应力作用而产生的时间依赖性变形现象。
蠕变本构方程描述了材料在蠕变过程中应力与应变之间的关系。
在有限元模拟计算中，利用蠕变本构方程可以计算热处理过程中的应力松弛行为。
蠕变本构方程通常包含一个弹性部分和一个蠕变部分。
弹性部分描述了材料在蠕变过程中的弹性回复，而蠕变部分描述了材料在蠕变过程中的变形行为。
常用的蠕变本构方程有Norton、Bailey-Norton、Power-Law等形式。
Norton本构关系是最常用的一种蠕变本构方程。
它描述了蠕变速率和应力之间的关系。
Norton本构方程的表达式为：ε˙(t) = Aσ^n exp(-Q/RT(t))其中，ε˙(t)是时间为t时刻的应变速率，A是材料的蠕变速率常数，σ是应力，n是蠕变指数，Q是活化能，R是气体常数，T(t)是时间为t时刻的温度。
Bailey-Norton本构关系是一种修正的Norton本构关系，它考虑了蠕变对材料的纹理演化的影响。
Bailey-Norton本构方程的表达式为：ε˙(t) = Aσ^n(εp, γ, T) exp(-Q/RT(t))其中，εp是塑性应变，γ是蠕变变形量，T(t)是时间为t时刻的温度。
Power-Law本构关系是一种广义的蠕变本构方程，它能够描述材料在不同温度和应力条件下的蠕变行为。
Power-Law本构方程的表达式为：ε˙(t) = Aσ^m(t)其中，ε˙(t)是时间为t时刻的应变速率，A是材料的蠕变常数，σ是应力，m(t)是时间为t时刻的蠕变指数。
在有限元模拟中，利用蠕变本构方程可以将材料的蠕变行为纳入计算，从而预测热处理过程中的应力松弛行为。
计算过程通常需要确定材料的材料参数，如蠕变速率常数、蠕变指数等，这需要通过实验测试或者文献数据获取。
总之，在有限元模拟计算中，利用蠕变本构方程可以计算热处理过程中的应力松弛行为。
