论文部分内容阅读
随着工业技术的发展,在现代常规电站、核电站及大型化工等行业中,工作介质的温度和压力大为提高,以超超临界汽轮机为例,其蒸汽的温度达到了600℃,压力达到了25MPa,在如此高的参数下长期运行,材料不可避免的要发生疲劳和蠕变破坏行为,当这两种行为同时发生时,产生一种具有更大破坏性的力学行为:棘轮变形。由于其平均应力不为零,导致在循环过程中累积塑性应变逐步增加,一直到材料破坏。在常规电站和核电站等行业中,压力容器的强度设计均要求进行棘轮效应分析。能否进一步提高高温材料的力学性能是限制超超临界汽轮机进一步提高参数的原因之一。材料力学性能的研究的主要内容之一是材料本构关系的研究,一个好的本构关系对于材料的寿命评估具有重要意义。目前大部分统一粘塑性本构模型(如Ohno-Wang、Chaboche等)在低应力区进行应力保持时,其粘性应力σ_v迅速衰减为零,导致塑性应变速率为零,因此不能对蠕变行为进行合理的描述。针对这一现象,提出在本构方程的塑性应变的演变项中加入Norton蠕变准则,蠕变项通过影响运动硬化和随动硬化实现与疲劳项的耦合。由于Norton方程只能描述蠕变的第一和第二阶段,体现不了第三阶段的蠕变加速阶段,根据勒梅特提出了应变等价原理,利用等效应力的概念引入蠕变损伤变量,并把等效应力应用到整个统一粘塑性本构方程中,在引入蠕变及蠕变损伤变量后,模型能够较为合理的描述棘轮与蠕变的交互作用。提出粒子群算法和有限元模块相结合的本构方程参数优化方法,使用面向对象的C++语言开发了有限元计算模块,与面向过程的Fortran语言相比,程序具有很好的扩展性能。提出的优化方法可以直接以实验数据作为输入得到模型参数,与传统的参数确定方法相比可以大大的减少工作量。针对有限元实现过程中的两个主要技术难点:隐式应力积分过程、一致性刚度矩阵推导,实现了耦合蠕变及损伤变量后的应力积分过程和刚度矩阵的表达式,并成功的嵌入到有限元软件ANSYS中。通过对内压管道受轴向循环载荷实验的模拟,结果表明,新的本构方程在描述棘轮蠕变交互作用方面的能力有所提高。