纳米单晶NiTi合金单程形状记忆效应的分子动力学模拟

采用基于第二近邻修正型嵌入原子势的分子动力学方法研究了纳米单晶NiTi合金的单程形状记忆效应,详细阐明了温度诱发马氏体相变和应力诱发马氏体重定向过程中纳米单晶的变形行为和微结构演化,进一步分析了加/卸载速率对NiTi合金单程形状记忆效应的影响。结果表明,NiTi纳米单晶在应力加载过程中发生马氏体重定向,卸载后存在残余应变;当加热到奥氏体转变结束温度以上时,马氏体逆相变为奥氏体相,残余应变逐渐减小,但未完全回复;随着应力加载速率的增加,重定向临界应力和模量逐渐增加;再次降温过程中不同加载速率下的原子结构演化各不相同。     

纳米单晶NiTi合金单程形状记忆效应的分子动力学模拟

采用基于第二近邻修正型嵌入原子势的分子动力学方法研究了纳米单晶NiTi合金的单程形状记忆效应，详细阐明了温度诱发马氏体相变和应力诱发马氏体重定向过程中纳米单晶的变形行为和微结构演化，进一步分析了加/卸载速率对NiTi合金单程形状记忆效应的影响。结果表明，NiTi纳米单晶在应力加载过程中发生马氏体重定向，卸载后存在残余应变；当加热到奥氏体转变结束温度以上时，马氏体逆相变为奥氏体相，残余应变逐渐减小，但未完全回复；随着应力加载速率的增加，重定向临界应力和模量逐渐增加；再次降温过程中不同加载速率下的原子结构演化各不相同。     

应力作用下NiTi形状记忆合金微结构演化的相场模拟及其本征应变率敏感性

基于Ginzburg-Landau动力学控制方程建立了NiTi形状记忆合金非等温相场模型,实现了对NiTi合金内应力诱导马氏体相变的数值模拟。同时将晶界能密度引入系统局部自由能密度,从而考虑多晶系统中晶界的重要作用。数值计算了单晶和多晶NiTi形状记忆合金在单轴机械载荷作用下微结构的动态演化过程和宏观力学行为,并重点研究了晶粒尺寸为60 nm的NiTi纳米多晶在低应变率下（0.0005～15 s?1）力学行为的本征应变率敏感性。研究结果表明,单晶NiTi合金系统高温拉伸-卸载过程中马氏体相变均匀发生,未形成奥氏体-马氏体界面。而纳米多晶系统在加载阶段出现了马氏体带的形成-扩展现象,在卸载阶段出现了马氏体带的收缩-消失现象。相同外载作用过程中,NiTi单晶系统的宏观应力-应变曲线具有更大的滞回环面积,拥有更优的超弹性变形能力。计算结果显示,在中低应变率下纳米晶NiTi形状记忆合金应力-应变关系表现出较明显的应变率相关性,应变率升高导致材料相变应力提升。这一应变率相关性主要源于相场模型中外加载荷速率与马氏体空间演化速度的相互竞争关系。     

梯度纳米晶NiTi形状记忆合金的超弹性和形状记忆效应相场模拟

通过建立考虑两个马氏体变体的二维相场模型,对梯度纳米晶镍钛(NiTi)合金系统的超弹性、单程和应力辅助双程形状记忆过程进行了模拟和预测.模拟结果显示: 在梯度纳米晶NiTi合金的超弹性过程中,较粗晶粒的区域保留了传统粗晶的马氏体相变和逆相变特征,即局部马氏体带的形核-扩展和缩减-消失, 而随着晶粒尺寸的减小,细晶粒区域表现为均匀相变的特点, 即无局部马氏体带产生; 此外,在超弹性和形状记忆过程中,马氏体相变和重取向都首先在较粗晶粒区域开始并逐步向细晶粒区域传播,而逆相变则相反.马氏体相变和重取向的逐步扩展使梯度纳米晶NiTi合金的应力-应变和应变-温度曲线呈现出“硬化状”,其可归因于纳米多晶NiTi合金中马氏体相变对晶粒尺寸的依赖性,即随着晶粒尺寸的减小, 相变或重取向壁垒逐渐增大,马氏体相变或重取向的形核、扩展越来越困难. 可见,梯度纳米晶结构具有比传统均匀晶粒尺寸NiTi合金更宽的相变应力区间、重取向应力区间和相变温度区间,可显著提高该合金非弹性变形的可控性.      

形状记忆合金伪弹性行为模拟

基于对伪弹性形状记忆合金(SMA)典型应力-应变曲线的特征分析,在原Graesser本构模型中增加简洁多项式来描述SMA应力诱发马氏体相变完成后在变形马氏体相下继续加载阶段的变形特征;并引入应变幅值与混相下SMA弹性模量的关系来改进不同应变幅值下卸载时SMA的应力-应变关系,从而提出了一种新的SMA一维本构关系模拟其伪弹性力学行为。该模型对直径为0.5mm的NiTi合金丝的拉伸加载、卸载试验曲线的模拟结果表明:改进本构模型与原Graesser模型相比,其能够准确地模拟SMA在不同应变幅值下加载和卸载应力-应变关系。此外,通过研究SMA本构模型的物理关系,推导出了控制SMA滞回曲线特征的关键参数fT与相变临界应力、弹性常数之间的明确关系,可利用该关系直接确定参数fT,摆脱了只靠试算获取该参数的传统做法,其准确性得到了试验验证。     

应力引起的镍钛单晶形状记忆合金相变的实验研究

在单晶形状记忆合金试样中，由于没有晶粒之间的约束，它的马氏体相界面移动比多晶容易，用实验方法研究其相变的特点，对建立新的理论模型有意义，因而对它的实验分析显得重要。本文利用高分辨率的CCD系统监测到NiTi单晶形状记忆合金在拉伸时的相变伪弹性的过程；利用X射线衍射法得到了NiTi单晶试样在拉伸方向的晶向；运用高分辨率的云纹干涉技术，获得了应力引起的NiTi单晶形状记忆合金相变时的变形场；利用高分辨率、高灵敏度的红外相机记录了NiTi单晶在拉伸状态下的温度变化规律；对低温下NiTi单晶的拉伸性能做了初步的研究，得到一些有意义的现象。     

NiTi合金力学和阻尼性质研究及伪弹性退化分析

试验考察了NiTi合金在不同应变幅值循环载荷作用下的力学特性和阻尼性能,并引入了残余内应力分析伪弹性退化机理。用马氏体相变开始应力、残余应变等参数表征NiTi合金伪弹性特征,用等效阻尼比表征NiTi合金的阻尼性能。试验结果分析表明:应变幅值增加会加快NiTi合金伪弹性随循环次数的退化;当应变幅值处于马氏体相变开始和结束应变之间时,不同应变幅值下NiTi合金的马氏体相变平台随循环次数增加同步降低,且当NiTi合金在相变中段卸载时其阻尼性能最好;结合残余内应力与残余应变正相关线性关系可分析NiTi合金伪弹性退化内在机理。该研究可为循环载荷下NiTi合金伪弹性行为的准确描述提供依据,并可为NiTi合金阻尼器的设计提供参考。     

NiTi形状记忆合金裂纹尖端热力耦合行为研究

本文对NiTi形状记忆合金I型裂纹尖端热力耦合行为进行了数值仿真分析和实验验证。建立了包含相变和热力耦合的本构模型,通过有限元计算得到了裂纹尖端附近的纵向应变、马氏体体积分数和温度场分布,依据马氏体相变情况对裂纹尖端有效应力强度因子进行了修正,揭示了加载速率对形状记忆合金裂纹尖端有效应力强度影子的影响规律。参数研究表明,随着加载频率的增加,裂纹尖端附近温度逐渐升高,马氏体相变区域逐渐缩小,有效应力强度因子呈下降趋势,形状记忆合金表现出增韧效应,有助于减缓裂纹扩展。本研究结果对于揭示热力耦合作用下超弹性形状记忆合金疲劳裂纹扩展规律具有重要参考意义。     

应变幅值对循环载荷下NiTi合金力学性质的影响

实验研究了应变幅值对循环载荷下NiTi合金伪弹性退化特征的影响规律,结果表明:当卸载发生在NiTi合金应力诱发马氏体相变阶段时,应变幅值对马氏体相变开始应力的退化规律影响较小,但此时应变幅值的增加会显著增大奥氏体弹性模量的退化程度,而其大变形可回复能力和阻尼特性在应变幅值大于6%时才有大幅度降低。对各参数退化程度进行定量分析,得到了NiTi合金具有较强可回复能力和阻尼性能的应变幅值范围。该研究可为NiTi合金阻尼器的设计提供参考。     

超弹性形状记忆合金管单向拉伸试验的数值模拟

NiTi形状记忆合金具有很强的超弹性行为，这种超弹性行为是由于材料在应力作用下发生可逆的马氏体相变所引起。最近Sun和Lee^[4]在NiTi形状记忆合金管单向拉伸试验中观测到，应力诱导马氏体相变具有螺旋带状的形貌特征，本文对此作了数值模拟研究。采用包含应变软化效应的三线性本构关系，建立了NiTi形状记忆合金管的三维有限元模型。通过迭代计算，成功地再现了试验中所观察到的螺旋状相变带从形成到长大的全过程。数值计算结果表明，产生这一独特现象的力学机制，在于NiTi形状记忆合金管在拉伸状态下出现的局部变形失稳极其传播。     

层状复合材料的弹塑性模型

通过将宏微观应力及应变分成面内及面外两部分,利用理想界面上的面外应力与面内应变连续条件,应用微观应力及应变均匀化方法得到了两层复合材料的弹塑性分析模型,并在此基础上构造了多层复合材料的全量型三维显式本构模型.各向同性层组成的复合材料算例表明了该模型的可靠性,同时该模型完全适用于各向异性层组成的复合材料弹塑性分析.     

SUPER-ELASTIC CONSTITUTIVE MODEL CONSIDERING PLASTICITY AND ITS FINITE ELEMENT IMPLEMENTATION

Based on the experimental results of super-elastic NiTi alloy, a three-dimensional super-elastic constitutive model including both of stress-induced martensite transformation and plasticity is constructed in a framework of general inelasticity. In the proposed model, transformation hardening, reverse transformation of stress-induced martensite, elastic mismatch between the austenite and martensite phases, and temperature-dependence of transformation stress and elastic modulus of each phase are considered. The plastic yielding of martensite occurred under high stress is addressed by a bilinear isotropic hardening rule. Drucker-Prager-typed transformation surfaces are employed to describe the asymmetric behavior of NiTi alloy in tension and compression. The prediction capability of the proposed model is verified by comparing the simulated results with the correspondent experimental ones. Based on backward Euler＇s integration, a new expression of consistent tangent modulus is derived. The proposed model is then implemented into a finite element package ABAQUS by user-subroutine UMAT. Finally, the validity of such implementation was verified by some numerical samples.     

Interfacial energy and dissipation in martensitic phase transformations. Part II: Size effects in pseudoelasticity

This paper is a continuation of the Part I (H. Petryk, S. Stupkiewicz, Interfacial energy and dissipation in martensitic phase transformations. Part I: Theory. J. Mech. Phys. Solids, 2010, doi:10.1016/j.jmps.2009.11.003). A fully three-dimensional model of an evolving martensitic microstructure is examined, taking into account size effects due to the interfacial energy and also dissipation related to annihilation of interfaces. The elastic micro-strain energy at microstructured interfaces is determined with the help of finite element computations and is approximated analytically. Three interface levels are examined: of grain boundaries attained by parallel martensite plates, of interfaces between austenite and twinned martensite, and of twin interfaces within the martensite phase. Minimization of the incremental energy supply, being the sum of the increments in the free energy and dissipation of the bulk and interfacial type at all levels, is used as the evolution rule, based on the theory presented in Part I. An example of the formation and evolution of a rank-three laminated microstructure of finite characteristic dimensions in a pseudoelastic CuAlNi shape memory alloy is examined quantitatively.     

考虑横法向热应变的C~0型Reddy板理论和三角形板单元

考虑横法向热变形,建议了C0型Reddy理论,并用于分析复合材料层合/夹层板热膨胀问题。虽然考虑了横法向热应变,但不增加额外的位移变量。此理论位移场不含有横向位移一阶导数,构造有限元时仅需C0插值函数。基于这一模型,运用虚位移原理推导了复合材料板平衡方程以及构造了6节点三角形板单元,并分析了简支复合材料层合/夹层板的热膨胀问题。数值结果表明,建立的模型能准确分析复合材料层合/夹层板热膨胀问题,而忽略横法向热应变的理论分析热膨胀问题误差较大。     

层合析参固支边界条件下的不稳定后屈曲分析

为了深入地研究复合材料层板所独具的后屈曲特性,利用能量变分原理和非线性几何方程建立了具有弹性约束的复合材料层板在面内载荷作用下的非线性稳定性控制方程组,并运用广义傅立叶级数法对其进行求解。重点分析了非对称层板在固支边界条件下的稳定性问题,发现层板在此条件下有可能存在非对称的失稳临界点和不稳定的后屈曲路径,进而构造了简化的物理模型进行解释,指出后屈曲的非对称性是由于结构关于Z轴不对称,而不稳定性是由于固支边界条件阻碍了前屈曲的发生。     

Microstructure in a Cubic to Orthorhombic Transition

Microstructures for a cubic to orthorhombic transition are constructed using a geometrically nonlinear, thermoelastic theory of martensitic transformations. Such microstructures are of interest because they provide low energy paths along which a specimen can transform. The particular microstructures considered are the twinned martensite, austenite–martensite, wedge, triangle, and diamond. More specifically, all possible twins are found along with the corresponding twinning elements and magnitude of the twin shear. Further, two kinds of austenite–martensite microstructures are studied: those with a single variant of martensite and those with twinned martensite. The regions in the space of transformation stretches in which each of these microstructures exist are determined, and the shape strains and habit plane normals are found as well. In addition, special microstructures, the wedge, triangle, and diamond, are constructed with both the austenite-single variant and austenite-twinned martensite microstructures. These special microstructures are of interest because they provide a mechanism through which the transformation may proceed more easily, and they are possible only in alloys with particular transformation stretches. Numerically computed level curves in the space of the stretches are presented on which the special microstructures are possible. These results may be useful in providing guidelines for alloy design.     

复合材料层合板自由振动分析的无网格自然邻接点Petrov-Galerkin法

基于一阶剪切变形理论,提出了复合材料层合板自由振动分析的无网格自然邻接点Petrov-Galerkin法。计算时在复合材料层合板中面上仅需要布置一系列的离散节点,并利用这些节点构建插值函数。在板中面上的局部多边形子域上,采用加权余量法建立复合材料层合板自由振动分析的离散化控制方程,并且这些子域可由Delaunay三角形方便创建。自然邻接点插值形函数具有Kronecker delta函数性质,因而无需经过特别处理就能准确地施加本质边界条件。对不同边界条件、不同跨厚比、不同材料参数和不同铺设角度的复合材料层合板,由本文提出的无网格自然邻接点Petrov-Galerkin法进行自由振动分析时均可得到满意的结果。数值算例结果表明,本文方法求解复合材料层合板的自由振动问题是行之有效的。     

Modelling of laminated microstructures in stress-induced martensitic transformations

This paper is concerned with micromechanical modelling of stress-induced martensitic transformations in crystalline solids, with the focus on distinct elastic anisotropy of the phases and the associated redistribution of internal stresses. Micro-macro transition in stresses and strains is analysed for a laminated microstructure of austenite and martensite phases. Propagation of a phase transformation front is governed by a time-independent thermodynamic criterion. Plasticity-like macroscopic constitutive rate equations are derived in which the transformed volume fraction is incrementally related to the overall strain or stress. As an application, numerical simulations are performed for cubic β₁ (austenite) to orthorhombic γ₁′ (martensite) phase transformation in a single crystal of Cu-Al-Ni shape memory alloy. The pseudoelasticity effect in tension and compression is investigated along with the corresponding evolution of internal stresses and microstructure.     

Refined Series Methodology for the Fully Coupled Thin-Walled Laminated Beams Considering Foundation Effects

The refined power series solutions are presented for the coupled static analysis of thin-walled laminated beams resting on elastic foundation. For this purpose, the elastic strain energy considering the material and structural coupling effects and the energy including the foundation effects are constructed. The equilibrium equations and the force-displacement relationships are derived from the extended Hamilton's principle, and the explicit expressions for displacement parameters are presented based on power series expansions of displacement components. Finally, the member stiffness matrix is determined by using the force-displacement relationships. For comparison, the finite element model based on the Hermite cubic interpolation polynomial is presented. In order to verify the accuracy and the superiority of the laminated beam element developed by this study, the numerical solutions are presented and compared with results obtained from the regular finite beam elements and the ABAQUS's shell elements. The influences of the fiber angle change and the boundary conditions on the coupled behavior of laminated beams with mono-symmetric I-sections are investigated.     

A new two-dimensional model for electro-mechanical response of thick laminated piezoelectric actuator

This paper investigates the electro-mechanical behaviour of a thick, laminated actuator with piezoelectric and isotropic lamina under externally applied electric loading using a new two-dimensional computational model. The elastic core is relatively thick and thus it is modelled by Timoshenko thick-beam theory. Although the piezoelectric lamina is a beam-like layer, it is formulated via a two-dimensional model because of not only the strong electro-mechanical coupling, but also of the presence of a two-dimensional electric field. It is shown in this paper that a one-dimensional model for the piezoelectric beam-like layer is inadequate. The piezoelectric model is constructed within the scope of linear piezoelectricity. The actuation response is induced through the application of external electric voltage. Under the strong coupling of elasticity and electricity, the strain energy and work of electric potential are presented. The electro-mechanical response of the laminated Timoshenko beam is formulated and determined via a variational energy principle. Numerical examples presented illustrate convincing comparison with finite element solutions and existing published data. New numerical solutions are also presented to investigate the geometric effect on the electro-mechanical bending behaviour.