高温合金材料循环相关和时间相关热机械疲劳循环性能模拟

从平衡热力学不可逆系统出发，用非线性粘弹塑性运动强化莱模拟高温合金材料的应变控制热机械疲劳循环特性。讨论了温度变化和应变循环的相位关系，循环相关和时间相关热机械疲劳损伤机制，蠕变和疲劳间的相互作用。在建立本构关系和状态方程时，均考虑了温度变化所产生的影响。     

冲击载荷下炸药装药动态响应的有限元分析及热点形成机理的数值模拟

对稳态温度场中受冲击载荷作用的炸药药柱进行了弹粘塑性分析。在Ｐｅｒｚｙｎａ本 构模型的基础上，作了适当的补充和修正，将流动参数、弹性模量Ｅ均视为温度的函数，动态 有限元计算结果表明，计算曲线和实验曲线有很好的近似。为模拟材料中不均匀性的影响，在 药柱中心引入一孔洞，有限元计算结果给出含孔洞药柱的粘塑性动态响应、药柱网格变形图以 及药柱等温线，可以清楚看出在孔洞附近区域有局部高温产生。本文的本构模型和计算方法对 于研究冲击载荷下炸药装药的力学响应以及炸药装药中热点形成机理的数值模拟提供了良好 的基础。     

含损伤材料的热粘塑性本构关系和柱壳破裂研究

以含内变量的本构关系理论为基础 ,结合材料损伤演化方程 ,并考虑了温度和损伤对材料参数的影响 ,得到了增量形式的热粘塑性本构关系的普适显式表达式。然后使用Bodner幂函数型粘塑性模型 ,具体推导了其增量形式的热粘塑性本构方程。接着结合在实践中有重要意义的内部爆炸载荷作用下的柱壳破裂问题 ,建立了含损伤热粘塑性柱壳破裂问题的完备方程组 ,使用有限差分方法 ,完成了对问题的数值模拟 ,并对结果进行了分析。计算结果与实验结果符合良好。     

冲击载荷下炸药装药动态响应的有限元分析及热点形成机理的数值模拟

对稳态温度场中受冲击载荷作用的炸药药柱进行了弹粘塑性分析。在Ｐｅｒｚｙｎａ本 构模型的基础上，作了适当的补充和修正，将流动参数γ、弹性模量Ｅ均视为温度的函数，动态 有限元计算结果表明，计算曲线和实验曲线有很好的近似。为模拟材料中不均匀性的影响，在 药柱中心引入一孔洞，有限元计算结果给出含孔洞药柱的粘塑性动态响应、药柱网格变形图以 及药柱等温线，可以清楚看出在孔洞附近区域有局部高温产生。本文的本构模型和计算方法对 于研究冲击载荷下炸药装药的力学响应以及炸药装药中热点形成机理的数值模拟提供了良好 的基础。     

高温合金材料循环相关热机械疲劳寿命预测

在变温非线性运动强化规律所描述的高温合金材料热机械寿命应力－应变循环特性的基础上，讨论了应变控制的循环相关热机械疲劳寿命预测技术，所建模型采用了由应变以密度表示的损伤参数，并且引入了温度损伤系数，考虑了温度变化范围以及温度循环和应变循环相位关系对疲劳寿命的影响，在确定模型的一些参数，采用等温力学试验和疲劳试验的数据，为了把等温疲劳研究成果推广到变温疲劳分析领域，开辟了新的途径。     

温度循环应力剖面对QFP焊点热疲劳寿命的影响

采用统一型粘塑性Anand本构方程描述了QFP(四方扁平封装)焊点的粘塑性力学行为,利用有限元分析软件建立组装在印制电路板上QFP的有限元模型,通过研究焊点内部总应变范围的变化进而研究温度循环应力剖面各个参数对焊点热疲劳寿命的影响,为设计合理的温度循环应力剖面提供了理论依据。     

一种考虑拉伸保持效应的多轴疲劳寿命模型

通过定义考虑拉伸保载效应的CFI因子(creep-fatigue interaction factor)，将拉伸蠕变损伤和疲劳损伤进行非线性耦合. 根据断裂实验的观察，针对拉伸主 导的裂纹萌生、扩展及破坏的多轴疲劳问题，给出了一个基于临界面方法的能量型高温多轴 疲劳寿命预测模型. 所给出的模型可对不同温度、不同载荷特点、不同保载时间的多轴疲劳 寿命进行预测，模型的材料参数不依赖于温度和载荷. 并且此方法可以很方便地推广到其它 因素主导破坏的高温多轴疲劳寿命预测. 通过拟合高温合金Udimet720Li单轴带保持时间的 低循环疲劳(low cycle fatigue, LCF)寿命试验数据，得到了材料常数. 结合黏 塑性有限元分析方法，对高温双轴带保载循环载荷下Cruciform试件的寿命进行了 预测，预测结果基本落在2倍分散带内，达到工程的要求，证明了该模型的有效性.     

材料的弹粘塑性本构关系和应力波的演化

试验表明，大多数工程材料在冲击载荷作用之下的变形一般都同时包含有可恢复的瞬态性弹性变形和不可恢复的粘滞性塑性变形，即其本构关系可以用弹粘塑性模型来描述。本文从内变量理论出发，探讨了时率相关材料的弹粘塑性本构关系的一般特性，建立了增量型的弹粘塑性本构关系的一般理论框架和普适的表达式，并且对两种最常用的本构模型——Bodner-Partom模型和Johnson-Cook模型给出了在一维应变条件下的具体形式。通过计算和讨论一维应变粘塑性靶板中冲击波的衰减机制和应力波的演化规律，特别是考察各种粘塑性本构模型中的材料参数对冲击波的衰减和应力波的演化的影响，得出了一些可以直接应用或具有一定借鉴价值的结果，为研究应力波的其他衰减机制以及在人防工程中智能防护层设计时新材料的选取奠定了基础。     

有限变形下单晶变温本构模型

本文构造了单晶热弹粘塑性的本构模型,模拟材料在不同温度下的力学行为。该模型以晶体热运动学作为分析变形的基础,即考虑温度变化情况下总体变形梯度的乘式分解,建立温度影响下的以弹性变形梯度为基本变量的控制方程来描述单晶材料的变形,算法采用隐式积分方法来求解控制方程以保证计算的稳定性。模型能反映单晶材料变形过程中温度对应力-应变响应的影响。     

连铸过程铸坯弹／黏塑性热应力的研究

本文推导并建立了高温铸坯黏塑性蠕变综合效应的非稳态热弹／黏塑性本构方程，在数值分析中考虑了对铸坯应力有重要影响的相变引起的体积变化，并编制了有限元程序。作通过对某厂板坯连铸机进行分析计算的结果表明，对板坯连铸过程铸坯热应力的模拟本模型和程序是有效的。     

高熵合金药型罩射流成型与稳定性

近年迅速兴起的多主元高熵合金因其具有很宽的成分-性能调控范围及一系列优异力学性能, 有望替代紫铜成为新一代药型罩材料. 本文基于CrMnFeCoNi五元高熵合金动态力学性能实验和数值模拟, 探索该合金用作药型罩的适用性. 基于分离式霍普金森拉杆和材料试验机研究了高熵合金不同应变率及温度下的力学行为, 获得了高熵合金Johnson-Cook热黏塑性动态本构. 利用流动速度与临界压垮角关系对凝聚性高熵合金射流形成边界进行界定. 结合数值模拟验证了高熵合金射流形成边界的合理性, 并进一步揭示了射流高速拉伸断裂演化规律. 研究表明: 射流断裂时间与材料强度成负相关, 材料动态强度增大, 将会引起射流断裂时间下降. 本工作可为新型高熵合金药型罩结构设计提供参考.      

SnPb钎料合金的粘塑性Anand本构方程

采用统一型粘塑性本构 Anand方程描述了电子封装焊点 Sn Pb钎料合金的非弹性变形行为 ,基于 Sn Pb 合金的弹塑性蠕变本构方程和实验数据 ,确定了6 2 Sn36 Pb2 Ag、6 0 Sn40 Pb、96 .5 Sn3.5 Ag和 97.5 Pb2 .5 Sn四种钎料合金 Anand方程的材料参数 ,验证了粘塑性 Anand本构方程对 Sn Pb合金在恒应变速率和稳态塑性流动条件下应力应变行为的预测能力。结果表明 ,Anand方程能有效描述 Sn Pb钎料的粘塑性本构行为 ,并可应用于电子封装 Sn Pb焊点的可靠性模拟和失效分析     

A thermodynamics based damage mechanics constitutive model for low cycle fatigue analysis of microelectronics solder joints incorporating size effects

Below certain length scales and in the presence of a non-uniform plastic strain field the mechanical behavior of many metals and its alloys is substantially different from that in bulk specimens. In particular, an increase in resistance with decreasing size has been observed in Pb/Sn eutectic solder alloys which are extensively used in microelectronics packaging interconnects. Due to the high homologous temperature, the Pb/Sn solder exhibits creep–fatigue interaction and significant time, temperature, stress and rate dependent material characteristics. The simultaneous consideration of all the above mentioned factors makes constitutive modeling an extremely difficult task. In this paper, a viscoplastic constitutive model unified with a thermodynamics based damage evolution model is embedded into a couple stress framework in order to simulate low cycle fatigue response coupled to size effects. The model is implemented into commercial finite element code ABAQUS. The microbending experiment on thin nickel foils is used to validate the model. Analyses are performed on a thin layer solder joint in bending under cyclic loading conditions.     

Modeling and validation of the large deformation inelastic response of amorphous polymers over a wide range of temperatures and strain rates

A robust physically consistent three-dimensional constitutive model is developed to describe the finite mechanical response of amorphous polymers over a wide range of temperatures and strain rates, including the rubbery region and for impact loading rates. This thermomechanical model is based on an elastic–viscoplastic rheological approach, wherein the effects of temperature, strain rate, and hydrostatic pressure are accounted for. Intramolecular, as well as intermolecular, interactions under large elastic–inelastic behavior are considered for the mechanisms of deformation and hardening. For a wide range of temperatures and strain rates, our simulated results for poly(methyl methacrylate) (PMMA) and polycarbonate (PC) are in good agreement with experimental observations.     

96.5Sn-3.5Ag钎料合金的时相关单轴应变循环变形行为及其本构关系研究

在室温下对96.5Sn-3.5Ag钎料合金进行了不同加载波形下的单轴应变循环实验。研究了在具有不同保持时间、不同应变率、不同应变幅值及其历史对材料的循环变形行为的影响。基于材料时相关变形行为,提出了统一粘塑性本构模型,并对该材料的变形行为进行本构模拟。实验结果表明:该钎料合金单轴变形行为具有应变率、保持时间以及应变幅值依赖性。本构关系的预言结果与实验结果吻合得一致性说明该种模型能够很好地描述材料的单轴应变循环变形行为。     

Numerical simulation and experimental investigation of the elastocaloric cooling effect in sputter-deposited TiNiCuCo thin films

The exploitation of the elastocaloric effect in superelastic shape memory alloys (SMA) for cooling applications shows a promising energy efficiency potential but requires a better understanding of the non-homogeneous martensitic phase transformation. Temperature profiles on sputter-deposited superelastic \({\mathrm {Ti_{55.2}Ni_{29.3}Cu_{12.7}Co_{2.8}}}\) shape memory alloy thin films show localized release and absorption of heat during phase transformation induced by tensile deformation with a strong rate dependence. In this paper, a model for the simulation of the thermo-mechanically coupled transformation behavior of superelastic SMA is proposed and its capability to reproduce the mechanical and thermal responses observed during experiments is shown. The procedure for experiment and simulation is designed such that a significant temperature change from the initial temperature is obtained to allow potential cooling applications. The simulation of non-local effects is enabled by the use of a model based on the one-dimensional Müller–Achenbach–Seelecke model, extended by 3D mechanisms such as lateral contraction and by non-local interaction, leading to localization effects. It is implemented into the finite element software COMSOL Multiphysics, and comparisons of numerical and experimental results show that the model is capable of reproducing the localized transformation behavior with the same strain rate dependency. Additionally to the thermal and the mechanical behavior, the quantitative prediction of cooling performance with the presented model is shown.     

An implicit consistent algorithm for the integration of thermoviscoplastic constitutive equations in adiabatic conditions and finite deformations

The so-called viscoplastic consistency model, proposed by Wang, Sluys and de Borst, is extended here to the integration of a thermoviscoplastic constitutive equation for J₂ plasticity and adiabatic conditions. The consistency condition in this case includes not only strain rate but also the effect of temperature on the yield function. Using the backward Euler integration scheme to integrate the constitutive equations, an implicit algorithm is proposed, leading to generalized expressions of the classical return mapping algorithm for J₂ plasticity, both for the iterative calculation of the plastic multiplier increment and for the consistent tangent operator when strain rate and temperature are considered also as state variables of the hardening equation. The model was implemented in a commercial finite element code and its performance is demonstrated with the numerical simulation of four Taylor impact tests.     

Temperature dependent bulge test for elastomers

The bulge test is a particularly convenient testing method for characterizing elastomers under biaxial loading. In addition, it is convenient to utilize this test for validating material models in simulation due to the heterogeneous strain field induced during inflation. During the bulge test the strain field for elastomers covers uniaxial tension at the border to pure shear and equibiaxial tension at the pole. Elastomeric materials exhibit a hyperelastic material behavior, with a dependency on temperature and loading rate. The temperature effect on the mechanical behavior during biaxial loading is considered in the present study. A bulge test setup combined with a temperature chamber is developed in order to characterize this effect, and an exemplary temperature dependent characterization of a poly(norbornene) elastomer is performed with this setup. The equibiaxial stress–strain curves measured at 60 °C, 20 °C and −20 °C are presented.     

Linear stability analysis for a thermoviscoplastic material under cyclic axial loading

Some mechanical properties exhibit a very strong dependence upon temperature; these evolutions can be properly analyzed by the steady state response in cyclic loading. To relate experimental conditions to thermomechanical characteristics, the existence and the stability of steady state solutions are studied for cylinders submitted to cyclic compression. The material, considered as rigid viscoplastic, is modeled by a non-Newtonian temperature dependent viscous law. Closed form solutions are obtained in the framework of a large deformation theory by neglecting thermal expansion and inertia effects. Steady state regime is analyzed. The stress versus strain rate response and the temperature distribution are established as functions of the geometry of the cylinder, the loading characteristics and the material parameters. The stability of steady state solutions is analyzed with use of a linear perturbation scheme.Received: 4 July 2002, Accepted: 5 August 2004, Published online: 24 February 2005PACS: 46.15.Ff, 83.60.St Correspondence to: F. Dinzart