Modeling of stress relaxation associated with an electrically-induced phase transformation

Stress relaxation due to an electrically-induced phase transformation in a ferroelectric crystal bar stretched by a hard-loading device is studied in the one-dimensional setting of electromechanics. According to the proposed model, the stress relaxation rate is governed by a nonlinear ordinary differential equation which resulted from the kinetic relation that controls the evolution of the phase transformation. A numerical analysis basing on simplifying approximation indicates that the stress starts to decrease when the intensity of the applied electric field reaches a critical value and that a stronger electric field results in a quicker stress reduction.     

形状记忆合金相变过程三维大变形有限元模拟

形状记忆合金（SMA）一直被作为智能材料开发,并被用于阻尼器、促动器和智能传感器元件．形状记忆合金（SMA）的一项重要特性,是它具有恢复在机械加卸载周期下产生的大变形而不表现出永久变形的能力．该文旨在介绍一种由应力产生的相变且可以描述马氏体和奥氏体之间的超弹性滞回环现象本构方程．形状记忆合金的马氏体系数假设为应力偏张量的函数,因此形状记忆合金在相变过程中锁定体积．本构模型是在大变形有限元的基础上执行的,采用了现时构型Lagrange大变形算法．为了方便地使用Cauchy应力和线性应变本构关系,使用了与旋转无关的Jaumann应力增率计算应力．数值分析结果表明,相变引起的超弹性滞回环可以有效地通过该文提出的本构方程和大变形有限元模拟．     

The numerical solution of Stefan problems with front-tracking and smoothing methods

The numerical performance of some computer methods for heat transfer with change of phase is discussed. For one-dimensional problems the application of invariant imbedding to time-discretized problems is suggested. For some multidimensional problems an absorption of the phase transition process into the diffusion equation through the so-called enthalpy transformation is advocated. If this transformation is not applicable, a locally one-dimensional Gauss-Seidel-type front-tracking method coupled with invariant imbedding is effective.     

Numerical model for vibration damping resulting from the first-order phase transformations

A numerical model is constructed for modelling macroscale damping effects induced by the first-order martensite phase transformations in a shape memory alloy rod. The model is constructed on the basis of the modified Landau–Ginzburg theory that couples nonlinear mechanical and thermal fields. The free energy function for the model is constructed as a double well function at low temperature, such that the external energy can be absorbed during the phase transformation and converted into thermal form. The Chebyshev spectral methods are employed together with backward differentiation for the numerical analysis of the problem. Computational experiments performed for different vibration energies demonstrate the importance of taking into account damping effects induced by phase transformations.     

金属材料热处理过程中的瞬态温度场与解态相变的数值计算方法

本文首先将Kirchhoff变换推广到导热系数为温度的多项式的非定常非线性热传导问题.并用分析方法确定热传导问题的边界条件.其次提出以孕育期叠加法并引用线性混合法则来模拟金属热处理过程的多相瞬态相变,较为简便地确定相变的开始时间、相变的种类及相变组织的数量.最后利用三维双重边界元法分析工件多种形式的热处理全过程.算例的数值计算结果表明本文方法行之有效.     

A heuristic approach for deriving the priority vector in AHP

This paper proposes an approach for deriving the priority vector from an inconsistent pair-wise comparison matrix through the nearest consistent matrix and experts judgments, which enables balancing the consistency and experts judgments. The developed algorithm for achieving a nearest consistent matrix is based on a logarithmic transformation of the pair-wise comparison matrix, and follows an iterative feedback process that identifies an acceptable level of consistency while complying with experts preferences. Three numerical examples are examined to illustrate applications and advantages of the developed approach.     

高压气体淬火过程中的温度场与相变耦合问题的计算

高压气体淬火技术是一种现代的、有效的材料加工技术．在Cheng所得到的高压气体淬火过程中非线性表面换热系数的基础上,用有限单元法对钢淬火过程的温度场与相变耦合问题进行了模拟计算．在数值计算中,材料的热物性系数被处理为温度和相变体积百分比的函数．为避免数值解的震荡,采用了Norsette有理近似法．     

A Levin-type algorithm for accelerating the convergence of Fourier series

A nonlinear sequence transformation is presented which is able to accelerate the convergence of Fourier series. It is tailored to be exact for a certain model sequence. As in the case of the Levin transformation and other transformations of Levin-type, in this model sequence the partial sum of the series is written as the sum of the limit (or antilimit) and a certain remainder, i.e., it is of Levin-type. The remainder is assumed to be the product of a remainder estimate and the sum of the first terms oftwo Poincaré-type expansions which are premultiplied by two different phase factors. This occurrence of two phase factors is the essential difference to the Levin transformation. The model sequence for the new transformation may also be regarded as a special case of a model sequence based on several remainder estimates leading to the generalized Richardson extrapolation process introduced by Sidi. An algorithm for the recursive computation of the new transformation is presented. This algorithm can be implemented using only two one-dimensional arrays. It is proved that the sequence transformation is exact for Fourier series of geometric type which have coefficients proportional to the powers of a numberq, |q|<1. It is shown that under certain conditions the algorithm indeed accelerates convergence, and the order of the convergence is estimated. Finally, numerical test data are presented which show that in many cases the new sequence transformation is more powerful than Wynn's epsilon algorithm if the remainder estimates are properly chosen. However, it should be noted that in the vicinity of singularities of the Fourier series the new sequence transformation shows a larger tendency to numerical instability than the epsilon algorithm.     

On the multiple hyperbolic systems modelling phase transformation kinetics

We discuss Cahn’s time cone method modelling phase transformation kinetics. The model equation by the time cone method is an integral equation in the space-time region. First, we reduce it to a system of hyperbolic equations, and in the case of odd spatial dimensions, the reduced system is a multiple hyperbolic equation. Next, we propose a numerical method for such a hyperbolic system. By means of alternating direction implicit methods, numerical simulations for practical forward problems are implemented with satisfactory accuracy and efficiency. In particular, in the three dimensional case, our numerical method on the basis of reduced multiple hyperbolic equation is fast.     

The double-exponential transformation in numerical analysis

The double-exponential transformation was first proposed by Takahasi and Mori in 1974 for the efficient evaluation of integrals of an analytic function with end-point singularity. Afterwards, this transformation was improved for the evaluation of oscillatory functions like Fourier integrals. Recently, it turned out that the double-exponential transformation is useful not only for numerical integration but also for various kinds of Sinc numerical methods. The purpose of the present paper is to review the double-exponential transformation in numerical integration and in a variety of Sinc numerical methods.     

核素随水-悬浮泥沙-底泥沙迁移模型

考虑了泥沙吸附和解吸、沉降和再悬浮对水相和泥沙相核素浓度的影响,根据质量守恒定律、吸附动力学方程和泥沙工程学建立了河流核素迁移的整体模型,进而得到分相模型.运用所得的模型对资江核电站排出的Sr在河流中迁移转化进行模拟,模拟结果表明所建模型基本上能反映Sr在水体中随水-悬浮泥沙-底泥沙的迁移规律.     

The ultra weak variational formulation for the modified mild-slope equation

This paper is devoted to evaluate a new alternative numerical method to capture accurately the diffraction-refraction process of waves in coastal areas. The modified mild-slope equation (MMSE) has been used to predict the water wave transformation when waves approach the shoreline. We perform numerical simulations in order to illustrate the efficiency of the ultra weak variational formulation (UWVF) method in comparison with the finite elements method (FEM). The UWVF method uses plane wave solutions on each element and has been shown to reduce the computational complexity at high wave numbers. We also present an alternative method to seek the angle of attack of the wave front on the domain boundary and show that the UWVF method reproduces effectively the numerical experimental data. We compare the FEM and the UWVF method for the MMSE and show that the UWVF method solves some of the difficulties that arise when the FEM is applied.     

Simulation of strain induced austenite-martensite transformation

Grain refinement due to phase transformation is an effective method for improving the mechanical properties of steel. An approach is proposed in the present work based on the FEM, for numerical simulation of the microstructure evolution as a result of hot rolling and subsequent cold torsion. Grain refinement in 304 stainless steel at four different technological schedules is considered. Results of numerical simulation are compared with experimental data. Coupling of the thermoplastic deformation with microstructure evolution is realized. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

保结构算法的相位误差分析及其修正

辛算法和保能量算法是应用最为广泛的两种保结构算法.本文从相位误差的角度给出了他们的比较结果.我们针对线性动力系统,分别分析了基于Pade对角逼近给出的辛算法和基于平均向量场法得到的能量守恒算法的相位误差,并通过数值验证了分析结果.文章还给出了保结构算法相位误差的改进方法,并通过数值例子验证了方法的有效性.     

A Regularized Sharp Interface Model for Phase Transformation Accounting for Prescribed Sharp Interface Kinetics

The macroscopic mechanical behavior of multi-phasic materials depends on the formation and evolution of their microstructure by means of phase transformation. In case of martensitic transformations, the resulting phase boundaries are sharp interfaces. We carry out a geometrically motivated discussion of the regularization of such sharp interfaces by use of an order parameter/phase-field and exploit the results for a regularized sharp interface model for two-phase elastic materials with evolving phase boundaries. To account for the dissipative effects during phase transition, we model the material as a generalized standard medium with energy storage and a dissipation function that determines the evolution of the regularized interface. Making use of the level-set equation, we are thereby able to directly translate prescribed sharp interface kinetic relations to the constitutive model in the regularized setting. We develop a suitable incremental variational three-field framework for the dissipative phase transformation problem. Finally, the modeling capability and the associated numerical solution techniques are demonstrated by means of a representative numerical example. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Algorithms for the regularization of ill-conditioned least squares problems

Two regularization methods for ill-conditioned least squares problems are studied from the point of view of numerical efficiency. The regularization methods are formulated as quadratically constrained least squares problems, and it is shown that if they are transformed into a certain standard form, very efficient algorithms can be used for their solution. New algorithms are given, both for the transformation and for the regularization methods in standard form. A comparison to previous algorithms is made and it is shown that the overall efficiency (in terms of the number of arithmetic operations) of the new algorithms is better.     

Solitary wave interaction for a higher-order nonlinear Schrodinger equation

Solitary wave interaction for a higher-order version of thenonlinear Schrödinger (NLS) equation is examined. An asymptotictransformation is used to transform a higher-order NLS equationto a higher-order member of the NLS integrable hierarchy, ifan algebraic relationship between the higher-order coefficientsis satisfied. The transformation is used to derive the higher-orderone- and two-soliton solutions; in general, the N-soliton solutioncan be derived. It is shown that the higher-order collisionis asymptotically elastic and analytical expressions are foundfor the higher-order phase and coordinate shifts. Numericalsimulations of the interaction of two higher-order solitarywaves are also performed. Two examples are considered, one satisfiesthe algebraic relationship derived from asymptotic theory, andthe other does not. For the example which satisfies the algebraicrelationship, the numerical results confirm that the collisionis elastic. The numerical and theoretical predictions for thehigher-order phase and coordinate shifts are also in strongagreement. For the example which does not satisfy the algebraicrelationship, the numerical results show that the collisionis inelastic and radiation is shed by the solitary wave collision.As the bed of radiation shed by the waves decays very slowly(like t^–), it is computationally infeasible to calculatethe final phase and coordinate shifts for the inelastic example.An asymptotic conservation law is derived and used to test thefinite-difference scheme for the numerical solutions.     

Stability and Noise-induced Transitions in an Ensemble of Nonlocally Coupled Chaotic Maps

The influence of noise on chimera states arising in ensembles of nonlocally coupled chaotic maps is studied. There are two types of chimera structures that can be obtained in such ensembles: phase and amplitude chimera states. In this work, a series of numerical experiments is carried out to uncover the impact of noise on both types of chimeras. The noise influence on a chimera state in the regime of periodic dynamics results in the transition to chaotic dynamics. At the same time, the transformation of incoherence clusters of the phase chimera to incoherence clusters of the amplitude chimera occurs. Moreover, it is established that the noise impact may result in the appearance of a cluster with incoherent behavior in the middle of a coherence cluster.     

Optimal solution for novel grey polynomial prediction model

The grey prediction model, as a time-series analysis tool, has been used in various fields only with partly known distribution information. The grey polynomial model is a novel method to solve the problem that the original sequence is in accord with a more general trend rather than the special homogeneous or non-homogeneous trend, but how to select the polynomial order still needs further study. In this paper the tuned background coefficient is introduced into the grey polynomial model and then the algorithmic framework for polynomial order selection, background coefficient search and parameter estimation is proposed. The quantitative relations between the affine transformation of accumulating sequence and the parameter estimates are deduced. The modeling performance proves to be independent of the affine transformation. The numerical example and application are carried out to assess the modeling efficiency in comparison with other conventional models.     

Irreversible Phase Transitions in Steel

We present a mathematical model for the austenite–pearlite and austenite–martensite phase transitions in eutectoid carbon steel. The austenite–pearlite phase change is described by the Additivity Rule. For the austenite–martensite phase change we propose a new rate law, which takes into account its irreversibility. We investigate questions of existence and uniqueness for the three-dimensional model and finally present numerical calculations of a continuous cooling transformation diagram for the eutectoid carbon steel C1080. © 1997 by B.G. Teubner Stuttgart-John Wiley & Sons, Ltd.