An analytical study on the instability phenomena during the phase transitions in a thin strip under uniaxial tension

In the experiments on stress-induced phase transitions in SMA strips, several interesting instability phenomena have been observed, including a necking-type instability (associated with the stress drop), a shear-type instability (associated with the inclination of the transformation front) and an orientation instability (associated with the switch of the inclination angle). In order to shed more lights on these phenomena, in this paper we conduct an analytical study. We consider the problem in a three-dimensional setting, which implies that one needs to study the difficult problem of solution bifurcations of high-dimensional nonlinear partial differential equations. By using the smallness of the maximum strain, the thickness and width of the strip, we use a methodology, which combines series expansions and asymptotic expansions, to derive the asymptotic normal form equations, which can yield the leading-order behavior of the original three-dimensional field equations. An important feature of the second normal form equation is that it contains a turning point for the localization (necking) solution of the first equation. It is the presence of such a turning point which causes the inclination of the phase transformation front. The WKB method is used to construct the asymptotic solutions, which can capture the shear instability and the orientation instability successfully. Our analytical results reveal that the inclination of the transformation front is a phenomenon of localization-induced buckling (or phase-transition-induced buckling as the localization is caused by the phase transition). Due to the similarities between the development of the Luders band in a mild steel and the stress-induced transformations in a SMA, the present results give a strong analytical evidence that the former is also caused by macroscopic effects instead of microscopic effects. Our analytical results also reveal more explicitly the important roles played by the geometrical parameters.     

Contact Symmetry Algebras of Scalar Ordinary Differential Equations

Using Lie's classification of irreducible contact transformations in thecomplex plane, we show thata third-order scalar ordinary differential equation (ODE)admits an irreducible contact symmetry algebra if and only if it is transformableto q ⁽³⁾=0 via a local contact transformation. This result coupled with the classification of third-order ODEs with respect to point symmetriesprovide an explanation of symmetry breaking for third-order ODEs. Indeed, ingeneral, the point symmetry algebra of a third-order ODE is not asubalgebra of the seven-dimensional point symmetry algebra of q ⁽³⁾=0.However, the contact symmetry algebra of any third-order ODE, except forthird-order linear ODEs with four- and five-dimensional pointsymmetry algebras, is shown to be a subalgebra of the ten-dimensional contact symmetryalgebra of q ⁽³⁾==0. We also show that a fourth-orderscalar ODE cannot admit an irreducible contact symmetry algebra. Furthermore, weclassify completely scalar nth-order (n5) ODEs which admitnontrivial contact symmetry algebras.     

含参数常微分方程边值问题的改进COLSYS算法

对常微分方程求解器ＣＯＬＳＹＳ的算法进行了分析、研究和改进,使其能够对参数方程和常微分方程的联立体系直接求解,一方面,免除了目前不得不采用的平凡常微分方程技巧所带来的大量不必要的计算,提高了效率;另一方面,由于减小了问题的规模,使条件数降低,从而改进了数值稳定性和精度。     

插值矩阵法分析双材料平面V形切口奇异阶

对二维V形切口问题提出奇异阶分析的一个新方法.首先,以V形切口尖端附近位移场沿其径向渐近展开为基础,将其线弹性理论控制方程转换成切口尖端附近关于周向变量的常微分方程组特征值问题,然后将数值求解两点边值问题的插值矩阵法进一步拓展为求解一般常微分方程组特征值问题,插值矩阵法是在离散节点上采用微分方程中待求函数的最高阶导数作为基本未知量.由此,V形切口的应力奇性阶问题通过插值矩阵法获得,同时相应的切口附近位移场和应力场特征向量一并求出.     

Constitutive modeling of shape memory alloy wire with non-local rate kinetics

A constitutive modeling approach for shape memory alloy (SMA) wire by taking into account the microstructural phase inhomogeneity and the associated solid–solid phase transformation kinetics is reported in this paper. The approach is applicable to general thermomechanical loading. Characterization of various scales in the non-local rate sensitive kinetics is the main focus of this paper. Design of SMA materials and actuators not only involve an optimal exploitation of the hysteresis loops during loading–unloading, but also accounts for fatigue and training cycle identifications. For a successful design of SMA integrated actuator systems, it is essential to include the microstructural inhomogeneity effects and the loading rate dependence of the martensitic evolution, since these factors play predominant role in fatigue. In the proposed formulation, the evolution of new phase is assumed according to Weibull distribution. Fourier transformation and finite difference methods are applied to arrive at the analytical form of two important scaling parameters. The ratio of these scaling parameters is of the order of 10⁶ for stress-free temperature-induced transformation and 10⁴ for stress-induced transformation. These scaling parameters are used in order to study the effect of microstructural variation on the thermo-mechanical force and interface driving force. It is observed that the interface driving force is significant during the evolution. Increase in the slopes of the transformation start and end regions in the stress–strain hysteresis loop is observed for mechanical loading with higher rates.      

Pulses in binary wave guide arrays and long wave PDE approximations

Binary waveguide arrays are linear arrays of optical waveguides with binary alternation of parameters, and have been of recent interest. They can be modeled by systems of nonlinear ODEs with forms related to the discrete nonlinear Schrödinger equation. Such equations can also arise in semi-classical molecular models of polymers with excitable states in each monomer, and coupling between these.An important class of solutions arises from an initially highly localized signal, such as input to a single element of the array. Simulations show that for a wide array of parameter values and of such initial data, a pulse is generated that travels approximately as a traveling wave. After a suitable phase shift in the variables, this pulse quickly develops a slow spatial variation, leading to a long-wave approximation by a system of coupled third order PDEs; one each for nodes of even and odd indices.This system of PDEs is presented, and verified to quite accurately reproduce the pulse propagation seen in the ODE system; further there is often a strong tendency for the behavior of the two PDE components to converge, with a corresponding convergence of the even and odd index parts of the ODE system solution. The PDE model gives some indication of why this occurs.     

On the reduction of some second-order nonlinear ODEs in physics and mechanics to first-order nonlinear integrodifferential and Abel’s classes of equations

Second-order ordinary differential equations (ODEs) with strong nonlinear stiffness terms (cubic nonlinearities) governing wave motions, dynamic crack propagations, nonlinear oscillations etc. in physics and nonlinear mechanics are analyzed. Selecting as guide line a second-order nonlinear ODE of the form of the forced Duffing equation and using admissible functional transformations it is possible to reduce it to an equivalent first-order nonlinear integrodifferential equation. The reduced equation is exact. In the limits of small or large values of the parameter characterizing this nonlinear problem, it is shown that further reductions lead to a nonlinear ODE of the Abel classes. Taking into account the known exact analytic solutions of this equivalent equation it is proved that there does not exist an exact analytic solution of this type of equations. However, in cases when convenient functional relations connecting all parameters of the corresponding null equation and the characteristics of the driving force exist, approximate analytic solutions to the problem under consideration are provided.     

变截面高层筒体结构的有限元线法整体稳定和二阶分析

本文用有限元线法对变截面的高层简体结构进行空间整体稳定和二阶分析。先把实际框筒结构分段连续化为正交各向异性折板结构；用有限元线法，通过有限条元半离散化，取结线上位移为基本未知函数，由势能驻值原理建立稳定和二阶分析的常微分方程组；由常微分方程求解器直接求解。     

Exact solutions to drift-flux multiphase flow models through Lie group symmetry analysis

In the present paper, Lie group symmetry method is used to obtain some exact solutions for a hyperbolic system of partial differential equations (PDEs), which governs an isothermal no-slip drift-flux model for multiphase flow problem. Those symmetries are used for the governing system of equations to obtain infinitesimal transformations, which consequently reduces the governing system of PDEs to a system of ODEs. Further, the solutions of the system of ODEs which in turn produces some exact solutions for the PDEs are presented. Finally, the evolutionary behavior of weak discontinuity is discussed.     

插值矩阵法分析与复合材料界面相交的平面裂纹奇异应力场

提出了用插值矩阵法分析与各向异性材料界面相交的平面裂纹应力奇异性。基于V形切口尖端附近区域位移场渐近展开,将位移场的渐近展开式的典型项代入线弹性力学基本方程,得到关于平面内与复合材料界面相交的裂纹应力奇异性指数的一组非线性常微分方程的特征值问题,运用插值矩阵法求解,获得了平面内各向异性结合材料中与界面以任意角相交的裂纹尖端的应力奇异性指数随裂纹角的变化规律,数值计算结果与已有结果比较表明,本文方法具有很高的精度和效率。     

Lyapunov, Bohl and Sacker-Sell Spectral Intervals for Differential-Algebraic Equations

Lyapunov and exponential dichotomy spectral theory is extended from ordinary differential equations (ODEs) to nonautonomous differential-algebraic equations (DAEs). By using orthogonal changes of variables, the original DAE system is transformed into appropriate condensed forms, for which concepts such as Lyapunov exponents, Bohl exponents, exponential dichotomy and spectral intervals of various kinds can be analyzed via the resulting underlying ODE. Some essential differences between the spectral theory for ODEs and that for DAEs are pointed out. It is also discussed how numerical methods for computing the spectral intervals associated with Lyapunov and Sacker-Sell (exponential dichotomy) can be extended from those methods proposed for ODEs. Some numerical examples are presented to illustrate the theoretical results.      

Martensite Strain Memory in the Shape Memory Alloy Nickel-Titanium Under Mechanical Cycling

This paper describes an experimental study of stress-induced martensitic phase transformation in the SMA Nickel-Titanium. The rich local thermo-mechanical interactions that underlie transformation are examined using three-dimensional Digital Image Correlation (strain fields) and infrared imaging (thermal fields). We quantify the complex local interactions between released/absorbed latent heat and the extent of transformation, and explore the characteristics of the phase fronts and the evolution of martensitic volume fraction. We also quantify a strong strain memory in the martensite that forms in the wake of the phase transformation front. The accommodated strain in the martensite will remain constant during loading, even as the existing phase front propagates. There also exists a remarkable strain memory in the martensite that persists from cycle to cycle, indicating that the local elastic stress fields in the martensite are driven by a dislocation structure and martensitic nuclei that largely stabilize during the first loading cycle.     

Phase proportioning in CuAlBe shape memory alloys during thermomechanical loadings using electric resistance variation

The microstructure of shape memory alloys changes with the thermomechanical history of the material. During thermomechanical loading, austenite, thermally-induced martensite or stress-induced martensite can be simultaneously present in the material. In applications integrating SMA parts, utilization conditions seriously affect the microstructure and can generate macroscopic strain or stress. Consequently, during thermomechanical loadings, it is important to be able to proportion the different phases and consequently to understand the kinetic transformation. This is very useful in the development of constitutive equations. This study shows, by a series of tests, that the proposed experimental method, based on the measurement of the variation of electric resistance of CuAlBe wires, permits to determine the volume fraction of the different phases present in the material (i.e., austenite, stress-induced martensite and thermally-induced martensite). The proposed method is applied to the most common thermomechanical behavior met in engineering applications of shape memory alloys: pseudoelasticity, pseudoplasticity, recovery-stress and stress-assisted two-way shape memory effect. The proportioning method based on a mixture law integrating the resistivity of pure phases present in the SMA is first performed on different two-phase mixture cases and then applied to a three phase mixture case.     

Anisotropy of martensitic transformations in modeling of shape memory alloy polycrystals

Based on the knowledge of the anisotropy associated with the martensitic transformations obtained from tension/compression experiments with oriented CuAlNi single crystals, a simple constant stress averaging approach is employed to model the SMA polycrystal deformation behaviors. Only elastic and inelastic strains due to the martensitic transformation, variant reorientations in the martensite phase and martensite to martensite transformations in thermomechanical loads are considered. The model starts from theoretical calculation of the stress-temperature transformation conditions and their orientation dependence from basic crystallographic and material attributes of the martensitic transformations. Results of the simulations of the NiTi, NiAl, and Cu-based SMA polycrystals in stress–strain tests are shown. It follows that SMA polycrystals, even with randomly oriented grains, typically exhibit tension/compression asymmetry of the shape of the pseudoelastic σ−ε curves in transformation strain, transformation stress, hysteresis widths, character of the pseudoelastic flow and in the slope of temperature dependence of the transformation stresses. It is concluded that some macroscopic features of the SMA polycrystal behaviors originate directly from the crystallography of the undergoing MT's. The model shows clearly the crystallographic origin of these phenomena by providing a link from the crystallographic and material attributes of martensitic transformations towards the macroscopic σ−ε−T behaviors of SMA polycrystals.     

比例边界等几何分析方法Ⅰ:波导本征问题

提出比例边界等几何方法 (scaled boundary isogeometric analysis, SBIGA), 并用以求解波导本征值问题. 在比例边界等几何坐标变换的基础上, 利用加权余量法将控制偏微分方程进行离散处理, 半弱化为关于边界控制点变量的二阶常微分方程, 即 TE 波或 TM 波波导的比例边界等几何分析的频域方程以及波导动刚度方程, 同时利用连分式求解波导动刚度矩阵. 通过引入辅助变量进一步得出波导本征方程. 该方法只需在求解域的边界上进行等几何离散, 使问题降低一维, 计算工作量大为节约, 并且由于边界的等几何离散, 使得解的精度更高, 进一步节省求解自由度. 以矩形和 L 形波导的本征问题分析为例, 通过与解析解和其他数值方法比较, 结果表明该方法具有精度高、计算工作量小的优点.     

SMA短纤维增强复合材料在热载下的相交特性

在形状记忆合金（ＳＭＡ）复合材料研究中,相变特性的研究是一个主要的工作．基于Ｅｓｈｅｌｂｙ的等效夹杂模型和Ｍｏｒｉ和Ｔａｎａｋａ的场平均法,考虑到ＳＭＡ材料的强物理非线性,发展了增量型的等效夹杂模型（ＩｎｃｒｅｍｅｎｔａｌＥｑｕｉｖａｌｅｎｔＩｎｃｌｕｓｉｏｎＭｏｄｅｌ）．考虑在某一温度循环条件下讨论形状记忆合金短纤维增强的铝基复合材料在热载下的相变行为．特别研究了ＳＭＡ短纤维复合材料在变温过程中纤维几何尺寸、体积分数等参数对ＳＭＡ复合材料的相变行为和ＳＭＡ内残余应力等的影响．这些工作对于指导材料设计和了解ＳＭＡ复合材料热机械特性是颇有意义的．     

Three-dimensional analysis of transient thermal stresses in functionally graded plates

An analytical solution is presented for three-dimensional thermomechanical deformations of a simply supported functionally graded (FG) rectangular plate subjected to time-dependent thermal loads on its top and/or bottom surfaces. Material properties are taken to be analytical functions of the thickness coordinate. The uncoupled quasi-static linear thermoelasticity theory is adopted in which the change in temperature, if any, due to deformations is neglected. A temperature function that identically satisfies thermal boundary conditions at the edges and the Laplace transformation technique are used to reduce equations governing the transient heat conduction to an ordinary differential equation (ODE) in the thickness coordinate which is solved by the power series method. Next, the elasticity problem for the simply supported plate for each instantaneous temperature distribution is analyzed by using displacement functions that identically satisfy boundary conditions at the edges. The resulting coupled ODEs with variable coefficients are also solved by the power series method. The analytical solution is applicable to a plate of arbitrary thickness. Results are given for two-constituent metal-ceramic FG rectangular plates with a power-law through-the-thickness variation of the volume fraction of the constituents. The effective elastic moduli at a point are determined by either the Mori–Tanaka or the self-consistent scheme. The transient temperature, displacements, and thermal stresses at several critical locations are presented for plates subjected to either time-dependent temperature or heat flux prescribed on the top surface. Results are also given for various volume fractions of the two constituents, volume fraction profiles and the two homogenization schemes.     

Application of the extended Kantorovich method to the bending of clamped cylindrical panels

Applicability and performance of the extended Kantorovich method (EKM) to obtain highly accurate approximate closed form solution for bending analysis of a cylindrical panel is studied. Fully clamped panel subjected to both uniform and non-uniform loadings is considered. Based on the Love–Kirchhoff first approximation for thin shallow cylindrical panels, the governing equations of the problem in terms of three displacement components include a system of two second order and one forth order partial differential equations. The governing PDE system is converted to a double set of ODE systems by assuming separable functions for displacements together with utilization of the extended Kantorovich method. The resulted ODE systems are solved iteratively. In each iteration, exact closed form solutions are presented for both ODE systems. Rapid convergence and high accuracy of the method is shown for various examples. Both displacement and stress predictions show close agreement with other analytical and finite element analysis.     

On the modeling of equilibrium twin interfaces in a single-crystalline magnetic shape memory alloy sample. II: numerical algorithm

This is part II of this series of papers. The aim of the current paper was to solve the governing PDE system derived in part I numerically, such that the procedure of variant reorientation in a magnetic shape memory alloy (MSMA) sample can be simulated. The sample to be considered in this paper has a 3D cuboid shape and is subject to typical magnetic and mechanical loading conditions. To investigate the demagnetization effect on the sample’s response, the surrounding space of the sample is taken into account. By considering the different properties of the independent variables, an iterative numerical algorithm is proposed to solve the governing system. The related mathematical formulas and some techniques facilitating the numerical calculations are introduced. Based on the results of numerical simulations, the distributions of some important physical quantities (e.g., magnetization, demagnetization field, and mechanical stress) in the sample can be determined. Furthermore, the properties of configurational force on the twin interfaces are investigated. By virtue of the twin interface movement criteria derived in part I, the whole procedure of magnetic field- or stress-induced variant reorientations in the MSMA sample can be properly simulated.     

A thermomechanical model for a 1-D shape memory alloy wire with propagating instabilities

A thermomechanical boundary value problem and constitutive model are presented for a shape memory alloy (SMA) wire under uniaxial loading. The intent is to develop a one-dimensional continuum model of an SMA element that includes all the relevant thermomechanical couplings and is suitable for inclusion in finite element analyses. Thermodynamic relations are derived from phenomenological considerations consistent with recent experimental observations and are calibrated to a typical commercially available NiTi wire material. The model includes both temperature-induced and stress-induced transformations that are necessary to exhibit the shape memory effect and pseudoelastic behaviors. The model accommodates possible unstable mechanical behavior during stress-induced transformations by allowing softening transformation paths and including strain gradient effects. This should provide a tool to study propagating transformation fronts and localized latent heat transfer with the surroundings and a variety of interesting future structural applications, such as composites with embedded SMA elements.