循环硬化材料本构模型的隐式应力积分和有限元实现

针对新发展的、能够描述循环硬化行为应变幅值依赖性的粘塑性本构模型,讨论了它的数值实现方法。首先,为了能够对材料的循环棘轮行为（Ratcheting）和循环应力松弛现象进行描述,对已有的本构模型进行了改进;然后,在改进模型的基础上,建立了一个新的、全隐式应力积分算法,进而推导了相应的一致切线刚度（Consistent Tangent Modulus）矩阵的表达式;最后,通过ABAQUS用户材料子程序UMAT将上述本构模型进行了有限元实现,并通过一些算例对一些构件的循环变形行为进行了有限元数值模拟,讨论了该类本构模型有限元实现的必要性和合理性。     

金属材料在复杂应力状态下的塑性流动特性及本构模型

工程应用中,金属材料和结构往往处于复杂应力状态。材料的塑性行为会受到应力状态的影响,要精确描述材料在复杂应力状态下的塑性流动行为,必须在本构模型中考虑应力状态效应的影响。然而,由于在动态加载下材料的应变率效应和应力状态效应相互耦合、难以分离,给应力状态效应的研究和模型的建立造成很大困难。通过对Ti-6Al-4V钛合金材料开展不同加载条件下的力学性能测试,提出了一个包含应力三轴度和罗德角参数影响的新型本构模型,并通过VUMAT用户子程序嵌入ABAQUS/Explicit软件。分别采用新提出的塑性模型和Johnson-Cook模型对压剪复合试样的动态实验进行了数值模拟。结果表明,新模型不仅在对材料本构曲线的拟合方面具有较强的优势,而且由该模型所得到的透射脉冲和载荷-位移曲线均更加准确。因此,该模型能够更精确地描述和预测金属材料在复杂应力状态下的塑性流变行为。     

非定常微分型黏弹性本构模型

在应力作用下, 材料的力学参数随着微观结构的变化而变化, 需要考虑参数的时间效应. 利用黏滞系数随时间变化的黏性元件, 构造出非定常Maxwell模型、非定常Kelvin模型和非定常Zener模型. 求解非定常模型的微分型本构方程得到它们的松弛模量、蠕变柔量和卸载方程. 结果表明, 可以把常见的经验松弛函数和经验蠕变函数视为非定常微分型本构模型.      

一个反映土体复杂加卸荷路径的弹黏塑性模型

基于Perzyna 过应力理论,构造了一个土体各向异性弹黏塑性本构模型. 模型引入Wheeler 旋转硬化法则,能够较好的描述土体的初始各向异性和应力诱发各向异性. 借助ABAQUS 软件中UMAT 子程序接口,模型采用隐式积分算法——图形返回算法实现. 通过对一组复杂加卸荷路径的三轴不排水剪切试验(HKMD) 及一组分级加载三轴不排水蠕变试验(Sachville clay) 的模拟,表明本模型能够合理反映土体的率效应、应力松弛、蠕变及土体各向异性现象,验证了模型的有效性.     

单轴循环载荷下镁合金宏观塑性本构关系研究

基于微观织构演化,本文对镁合金的宏观塑性本构模型进行了改进,得到了室温下镁合金AZ31B和AZ61A单轴循环拉压载荷作用下的应力应变曲线。模型中区分了滑移、孪晶与去孪三种不同微观塑性变形机制,并且研究了它们对材料宏观循环塑性行为的影响。通过引入初始非零背应力来模拟循环应力应变曲线中拉压不对称性,给出了材料在单轴循环拉压载荷下背应力及屈服面的演化规律,并运用线性预测塑性回拉技术与隐式积分法求解相应的非线性方程组,编写了与.ABAQUS接口的UMAT用户子程序。预测与实验结果的比较显示该模型能较好地模拟镁合金材料在单轴循环拉压载荷作用下的力学行为,具有较高的准确性与良好的实用性。     

考虑运动硬化的DD3镍基合金力学行为研究

详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,详细介绍了镍基合金的晶体塑性本构模型,在Asaro大变形晶体塑性框架下,引入了运动硬化规律,考虑了温度和应变率对晶体塑性变形的影响,通过针对每个滑移系考虑屈服准则和流动规律建立了晶体塑性模型. 对积分过程进行了推导,通过编写ABAQUS材料用户子程序(UMAT), 实现本构模型的有限元积分算法. 在此基础上模拟了DD3镍基单晶合金在单轴拉伸和循环载荷下的响应,并与实验数据进行了对比. 利用该模型可以很好地模拟镍基单晶所具有的各向异性特性,体现了镍基单晶在循环载荷作用下的拉-压不对称性.      

率相关本构方程积分新算法

提出一种积分率相关本构方程的隐式积分新算法,引入0～1范围内的缩放因子λ对本构方程进行间接求解,可以避免直接求解等效塑性应变或等效塑性应变率时,由于其数值过大或过小而造成的收敛困难或收敛失败,实现对率相关本构方程的快速准确求解.以B-P统一本构方程及双曲正弦本构方程为例,验证了算法的可行性.结果表明,新算法对于准静态变形条件下的无硬化本构方程也可以得出准确的解.     

高应变率下航空透明聚氨酯的动态本构模型

采用低阻抗分离式霍普金森压杆对航空透明聚氨酯进行了高应变率下的动态力学性能测试,得到的应力应变曲线表现出了显著的非线性黏弹性特征。基于本构理论和实验数据,构建了航空透明聚氨酯的松弛时间应变相关的超黏弹性本构形式。该本构模型由2部分组成:一部分表征准静态下的超弹性行为,另一部分描述非线性应变率的相关特性。利用超黏弹性本构模型对不同应变率下航空透明聚氨酯的动态应力应变曲线进行拟合,拟合曲线与实验曲线一致性良好。     

黏弹-Perzyna黏塑性有限元法应力更新隐式算法

黏弹-黏塑性耦合模型的黏弹性部分由弹簧、黏壶和Kelvin链串联而成,黏塑性部分为双曲线型DruckerPrager屈服函数、各向同性硬化和Perzyna黏塑性流动模型。基于黏弹性蠕变柔度,通过定义与弹性问题相对应的与时间增量相关的黏弹性剪切模量和体积模量,导出增量递推形式的本构方程。为保证算法的收敛和稳定性,把Perzyna黏塑性流动方程转化为与弹塑性相似的一致性条件,建立黏塑性增量因子单侧逼近其收敛值的N-R迭代算法。最后,给出应力更新完全隐式算法和最终计算公式。分别采用黏弹性、黏弹-塑性和黏弹-黏塑性本构关系对一地基蠕变模型进行三维有限元分析和比较,结果表明,本文算法具有较高的计算效率和稳定性。     

一种修正的形状记忆合金本构模型数值算法及应用

采用Hughes-Winget算法修正了已有文献所发展的考虑塑性和相变耦合效应的形状记忆合金(SMA)本构模型的积分算法,使其能适用于发生较大转动变形的问题。据此编制了ABAQUS用户子程序UMAT,对SMA弹簧拉伸和"三点弯"等发生较大转动的问题进行了模拟。结果表明,修正算法可减小结构在发生较大转动变形时采用小变形本构模型计算带来的误差,提高计算过程的收敛速度与稳定性。采用本文的修正算法模拟了Ni Ti SMA的大变形拉伸伪弹性和塑性、SMA板的大挠度弯曲伪弹性和SMA弹簧的大变形拉伸伪弹性行为,结果与试验和其他研究者的计算结果吻合较好,证明了该修正算法的有效性。     

A finite deformation thermomechanical constitutive model for triple shape polymeric composites based on dual thermal transitions

Shape memory polymers (SMPs) have gained strong research interests recently due to their mechanical action that exploits their capability to fix temporary shapes and recover their permanent shape in response to an environmental stimulus such as heat, electricity, irradiation, moisture or magnetic field, among others. Along with interests in conventional “dual-shape” SMPs that can recover from one temporary shape to the permanent shape, multi-shape SMPs that can fix more than one temporary shapes and recover sequentially from one temporary shape to another and eventually to the permanent shape, have started to attract increasing attention. Two approaches have been used to achieve multi-shape shape memory effects (m-SMEs). The first approach uses polymers with a wide thermal transition temperature whilst the second method employs multiple thermal transition temperatures, most notably, uses two distinct thermal transition temperatures to obtain triple-shape memory effects (t-SMEs). Recently, one of the authors’ group reported a triple-shape polymeric composite (TSPC), which is composed of an amorphous SMP matrix (epoxy), providing the system the rubber-glass transition to fix one temporary shape, and an interpenetrating crystallizable fiber network (PCL) providing the system the melt-crystal transition to fix the other temporary shape. A one-dimensional (1D) material model developed by the authors revealed the underlying shape memory mechanism of shape memory behaviors due to dual thermal transitions. In this paper, a three-dimension (3D) finite deformation thermomechanical constitutive model is presented to enable the simulations of t-SME under more complicated deformation conditions. Simple experiments, such as uniaxial tensions, thermal expansions and stress relaxation tests were carried out to identify parameters used in the model. Using an implemented user material subroutine (UMAT), the constitutive model successfully reproduced different types of shape memory behaviors exhibited in experiments designed for shape memory behaviors. Stress distribution analyses were performed to analyze the stress distribution during those different shape memory behaviors. The model was also able to simulate complicated applications, such as a twisted sheet and a folded stick, to demonstrate t-SME.     

A coupled temperature and strain rate dependent yield function for dynamic deformations of bcc metals

A coupled temperature and strain rate microstructure physically based yield function is proposed in this work. It is incorporated along with the Clausius–Duhem inequality and an appropriate free energy definition in a general thermodynamic framework for deriving a three-dimensional kinematical model for thermo-viscoplastic deformations of body centered cubic (bcc) metals. The evolution equations are expressed in terms of the material time derivatives of the elastic strain, accumulated plastic strain (isotropic hardening), and the back stress conjugate tensor (kinematic hardening). The viscoplastic multipliers are obtained using both the Consistency and Perzyna viscoplasticity models. The athermal yield function is employed instead of the static yield function in the case of the Perzyna viscoplasticity model. It is found that the static strain rate value, at which the material shows rate-independent behavior, varies with the material deformation temperature. Computational aspects of the proposed model are addressed through the finite element implementation with an implicit stress integration algorithm. Finite element simulations are performed by implementing the proposed viscoplasticity constitutive models in the commercial finite element program ABAQUS/Explicit [ABAQUS, 2003. User Manual, Version 6.3. Habbitt, Karlsson and Sorensen Inc., Providence, RI] via the user material subroutine coded as VUMAT. Numerical implementation for a simple compression problem meshed with one element is used to validate the proposed model implementation with applications to tantalum, niobium, and vanadium at low and high strain rates and temperatures. The analysis of a tensile shear banding is also investigated to show the effectiveness and the performance of the proposed framework in describing the strain localizations at high velocity impact. Results show mesh independency as a result of the viscoplastic regularization used in the proposed formulation.     

Thickness dependence of viscoelastic stress relaxation of laminated microbeams due to mismatch strain

Time-dependent behaviors due to various mismatch strains are very important to the reliability of micro-/nano-devices. This paper aims at presenting an analytical model to study the viscoelastic stress relaxation of the laminated microbeam caused by mismatch strain. Firstly, Zhang’s two-variable method is used to establish a mechanical model for predicting the quasi-static stress relaxation of the laminated microbeam. Secondly, the related analytical solutions are obtained by combining the diffe...     

Finite deformation of a polymer: experiments and modeling

The response of a polymer (polytetrafluoroethylene) to quasi-static and dynamic loading is determined and modeled. The polytetrafluoroethylene is extremely ductile and highly nonlinear in elastic as well as plastic behaviors including elastic unloading. Constitutive model developed earlier by Khan, Huang and Liang (KHL) is extended to include the responses of polymeric materials. The strain rate hardening, creep, and relaxation behaviors of polytetrafluoroethylene were determined through extensive experimental study. Based on the observation that both viscoelastic and viscoplastic deformation of polytetrafluoroethylene are time dependent and nonlinear, a phenomenalogical viscoelasto–plastic constitutive model is presented by a series connection of a viscoelastic deformation module (represented by three elements standard solid spring dashpot model), and a viscoplastic deformation module represented by KHL model. The KHL module is affected only when the stress exceeds the initial yield stress. The comparison between the predictions from the extended model and experimental data for uniaxial static and dynamic compression, creep and relaxation demonstrate that the proposed constitutive model is able to represent the observed time dependent mechanical behavior of polytetrafluoroethylene polytetrafluoroethylene qualitatively and quantitatively.     

Cohesive zone model for high-cycle fatigue of adhesively bonded joints under mode I loading

A cohesive zone model adequate for simulating the behaviour of adhesively bonded joints subjected to high-cycle fatigue and pure mode I loading is presented. The bilinear cohesive zone law with linear softening relationship was considered. The main advantage of the proposed formulation is the use of a unique damage parameter accounting for cumulative damage resulting from static and fatigue loading. The method was implemented in a user subroutine of the commercial finite element software Abaqus®. Two-dimensional numerical simulations of the double cantilever beam test using different representative combinations of the modified Paris law coefficients were performed. It was verified that the results of the model simulate with excellent agreement the several Paris laws used as input, thus demonstrating the good performance of the method as a predictive tool.     

双向加载条件下尼龙6-橡胶复合材料的应力松弛研究

在双向测试系统上进行了不同纵向应变与不同横向应变的双向松弛实验,研究了在双向拉伸载荷作用下单向尼龙6-橡胶复合材料的应力松弛特性．为了预测尼龙6-橡胶复合材料的应力松弛规律,提出了一个松弛型本构模型．当试件承受双向拉伸载荷作用时,将松弛型本构模型获得的理论曲线和实验数据进行了对比,二者取得了较好的一致性．     

60Si2Mn螺旋弹簧的压缩应力松弛行为与贮存寿命预测

为评价60Si2Mn螺旋压缩弹簧的室温松弛特性,利用InstronE3000K8953型小吨位电子动静态疲劳试验机,对其在不同温度和初始应力水平条件下进行了高温压缩加速应力松弛试验,研究了环境温度、初始应力水平对松弛行为的影响.基于粘弹性体模型,揭示了应力松弛过程中弹性应变向塑性应变的转化特性与塑性应变随松弛时间的变化规律.在对应力松弛前后弹簧丝材金相和TEM微结构进行对比分析的基础上,探讨了应力松弛的微观机制.结果表明,环境温度与初始应力水平对松弛速率具有显著影响.基于应力松弛过程的热激活特性,建立了60Si2Mn螺旋压缩弹簧的贮存寿命预测方程,并对不同应力水平下弹簧的室温和高温贮存寿命进行了合理预测.     

一种考虑界面相物性影响的界面单元法

界面是由复杂的界面相简化而成的,界面破坏实际是界面相材料的破坏。数值计算为了方便,如经典模型和内聚力模型等,都把很薄的界面相作无厚度化处理。导致只能考虑界面的面力,而无法考虑界面相内的应力(平行于界面方向的应力)。使界面失效准则先天性地排除了界面相内部应力的影响,从界面相材料失效机理的角度来看这是不够严谨的。本文将界面相材料等效为一种弹性连续体,由界面本构关系推导得到了一种新的界面单元。该单元具有界面参数易确定、对界面相物性可以进行等效描述等优点。通过商用有限元软件ABAQUS和用户子程序UEL实现了数值分析,并与直接物理模型的数值模拟结果进行对比,证明了本方法的简便及准确性。通过对不同界面相厚度结构的进一步分析,探讨了本文方法的可行范围。     

基于微结构动态演化机制的单晶镍基高温合金晶体塑性本构及其有限元模拟

因其优异的高温力学性能,镍基单晶高温合金在航空航天和能源等领域得到了广泛的应用.镍基单晶高温合金优异的高温性能来源于其特有的两相微结构.基于代表体胞模型及分块均匀化方法,以位错密度为主要内变量,发展了一个包含两相微结构和位错演化信息的单晶镍基高温合金塑性行为的本构模型.该本构模型充分考虑了镍基单晶合金中位错在基体相和沉淀增强相中的多种演化机制,例如,基体位错八面体滑移、立方滑移、位错攀移、交滑移、位错弓出、位错切过沉淀增强相以及位错Kear-Wilsdolf(K-W)锁形成与解锁等.在商用有限元软件ABAQUS的框架下,编制了UMAT用户材料子程序.利用该用户子程序,对单晶和多晶镍基高温合金在不同温度、不同加载方向下的单调塑性、循环塑性、蠕变等典型行为进行了计算模拟.结果表明:该晶体塑性本构模型能"统一地"刻画镍基高温合金在不同温度、不同方向下的多种变形行为,并与实验结果具有良好的一致性.