A simple multi‐directional absorbing layer method to simulate elastic wave propagation in unbounded domains

The numerical analysis of elastic wave propagation in unbounded media may be difficult due to spurious waves reflected at the model artificial boundaries. This point is critical for the analysis of wave propagation in heterogeneous or layered solids. Various techniques such as Absorbing Boundary Conditions, infinite elements or Absorbing Boundary Layers (e.g. Perfectly Matched Layers) lead to an important reduction of such spurious reflections. In this paper, a simple absorbing layer method is proposed: it is based on a Rayleigh/Caughey damping formulation which is often already available in existing Finite Element softwares. The principle of the Caughey Absorbing Layer Method is first presented (including a rheological interpretation). The efficiency of the method is then shown through 1D Finite Element simulations considering homogeneous and heterogeneous damping in the absorbing layer. 2D models are considered afterwards to assess the efficiency of the absorbing layer method for various wave types and incidences. A comparison with the PML method is first performed for pure P‐waves and the method is shown to be reliable in a more complex 2D case involving various wave types and incidences. It may thus be used for various types of problems involving elastic waves (e.g. machine vibrations, seismic waves, etc.). Copyright © 2010 John Wiley & Sons, Ltd.     

Accurate radiation boundary conditions for the time‐dependent wave equation on unbounded domains

Asymptotic and exact local radiation boundary conditions (RBC) for the scalar time‐dependent wave equation, first derived by Hagstrom and Hariharan, are reformulated as an auxiliary Cauchy problem for each radial harmonic on a spherical boundary. The reformulation is based on the hierarchy of local boundary operators used by Bayliss and Turkel which satisfy truncations of an asymptotic expansion for each radial harmonic. The residuals of the local operators are determined from the solution of parallel systems of linear first‐order temporal equations. A decomposition into orthogonal transverse modes on the spherical boundary is used so that the residual functions may be computed efficiently and concurrently without altering the local character of the finite element equations. Since the auxiliary functions are based on residuals of an asymptotic expansion, the proposed method has the ability to vary separately the radial and transverse modal orders of the RBC. With the number of equations in the auxiliary Cauchy problem equal to the transverse mode number, this reformulation is exact. In this form, the equivalence with the closely related non‐reflecting boundary condition of Grote and Keller is shown. If fewer equations are used, then the boundary conditions form high‐order accurate asymptotic approximations to the exact condition, with corresponding reduction in work and memory. Numerical studies are performed to assess the accuracy and convergence properties of the exact and asymptotic versions of the RBC. The results demonstrate that the asymptotic formulation has dramatically improved accuracy for time domain simulations compared to standard boundary treatments and improved efficiency over the exact condition. Copyright © 2000 John Wiley & Sons, Ltd.     

A nonlinear plate finite element formulation for shape memory alloy applications

The aim of the present work is to develop a new finite element model for the finite strain analysis of plate structures constituted of shape memory alloy (SMA) material. A three‐dimensional constitutive model for shape memory alloys able to reproduce the special thermomechanical behavior of SMA characterized by pseudoelasticity and shape memory effects is adopted. The finite strain constitutive model is thermodynamically consistent and is completely formulated in the reference configuration. A two‐dimensional plate theory is proposed based on a tensor element shape function formulation. The displacement field is expressed in terms of increasing powers of the transverse coordinate. The equilibrium statement is formulated on the basis of the virtual displacement principle in a total Lagrangian format. The proposed displacement formulation is particularly suitable for the simple derivation of high‐order finite elements. Numerical applications are performed to assess the efficiency and locking performance of the proposed plate finite element. Some additional numerical examples are carried out to study the accuracy and robustness of the proposed computational technique and its capability of describing the structural response of SMA devices. Copyright © 2011 John Wiley & Sons, Ltd.     

An improved continued‐fraction‐based high‐order transmitting boundary for time‐domain analyses in unbounded domains

A high‐order local transmitting boundary to model the propagation of acoustic or elastic, scalar or vector‐valued waves in unbounded domains of arbitrary geometry is proposed. It is based on an improved continued‐fraction solution of the dynamic stiffness matrix of an unbounded medium. The coefficient matrices of the continued‐fraction expansion are determined recursively from the scaled boundary finite element equation in dynamic stiffness. They are normalised using a matrix‐valued scaling factor, which is chosen such that the robustness of the numerical procedure is improved. The resulting continued‐fraction solution is suitable for systems with many DOFs. It converges over the whole frequency range with increasing order of expansion and leads to numerically more robust formulations in the frequency domain and time domain for arbitrarily high orders of approximation and large‐scale systems. Introducing auxiliary variables, the continued‐fraction solution is expressed as a system of linear equations in iω in the frequency domain. In the time domain, this corresponds to an equation of motion with symmetric, banded and frequency‐independent coefficient matrices. It can be coupled seamlessly with finite elements. Standard procedures in structural dynamics are directly applicable in the frequency and time domains. Analytical and numerical examples demonstrate the superiority of the proposed method to an existing approach and its suitability for time‐domain simulations of large‐scale systems. Copyright © 2011 John Wiley & Sons, Ltd.     

Experimental study of one‐dimensional superelastic capability of a nearly equi‐atomic NiTi shape memory alloy

A series of experimental studies have been carried out on nearly equi‐atomic NiTi shape memory alloy wires. The effects of fatigue cycles, displacement rates as well as testing temperatures on the superelastic capabilities have been studied. Under cycling loading, the threshold stresses for martensitic transformation decrease and the residule strains increase. Saturation is reached after 100 cycles. With increasing displacement rates, the critical stresses required for the martensitic transformation increase and the slopes of the upper plateau in the stress‐strain curves rise. The dissipated energy increases rapidly with increasing displacement rate, reaches a maximum value at around 6.0mm/min and then decreases as the displacement rate continues to increase. Tests also show that the threshold stresses characterizing the forward/reverse transformations increase linearly with increasing test temperature.     

Numerical method for non‐linear wave and diffusion equations by the variational iteration method

This paper applies He's variational iteration to the wave equations in an infinite one‐dimensional medium and some non‐linear diffusion equations. A suitable choice of an initial solution can lead to the needed exact solution by a few iterations. Copyright © 2007 John Wiley & Sons, Ltd.     

A continued‐fraction‐based high‐order transmitting boundary for wave propagation in unbounded domains of arbitrary geometry

A high‐order local transmitting boundary is developed to model the propagation of elastic waves in unbounded domains. This transmitting boundary is applicable to scalar and vector waves, to unbounded domains of arbitrary geometry and to anisotropic materials. The formulation is based on a continued‐fraction solution of the dynamic‐stiffness matrix of an unbounded domain. The coefficient matrices of the continued fraction are determined recursively from the scaled boundary finite element equation in dynamic stiffness. The solution converges rapidly over the whole frequency range as the order of the continued fraction increases. Using the continued‐fraction solution and introducing auxiliary variables, a high‐order local transmitting boundary is formulated as an equation of motion with symmetric and frequency‐independent coefficient matrices. It can be coupled seamlessly with finite elements. Standard procedures in structural dynamics are directly applicable for evaluating the response in the frequency and time domains. Analytical and numerical examples demonstrate the high rate of convergence and efficiency of this high‐order local transmitting boundary. Copyright © 2007 John Wiley & Sons, Ltd.     

Non‐reflecting artificial boundaries for transient scalar wave propagation in a two‐dimensional infinite homogeneous layer

This paper presents an exact non‐reflecting boundary condition for dealing with transient scalar wave propagation problems in a two‐dimensional infinite homogeneous layer. In order to model the complicated geometry and material properties in the near field, two vertical artificial boundaries are considered in the infinite layer so as to truncate the infinite domain into a finite domain. This treatment requires the appropriate boundary conditions, which are often referred to as the artificial boundary conditions, to be applied on the truncated boundaries. Since the infinite extension direction is different for these two truncated vertical boundaries, namely one extends toward x →∞ and another extends toward x→‐ ∞, the non‐reflecting boundary condition needs to be derived on these two boundaries. Applying the variable separation method to the wave equation results in a reduction in spatial variables by one. The reduced wave equation, which is a time‐dependent partial differential equation with only one spatial variable, can be further changed into a linear first‐order ordinary differential equation by using both the operator splitting method and the modal radiation function concept simultaneously. As a result, the non‐reflecting artificial boundary condition can be obtained by solving the ordinary differential equation whose stability is ensured. Some numerical examples have demonstrated that the non‐reflecting boundary condition is of high accuracy in dealing with scalar wave propagation problems in infinite and semi‐infinite media. Copyright © 2003 John Wiley & Sons, Ltd.     

形状记忆合金超弹性本构关系的神经网络模型

通过试验，研究了受过循环变形、具有稳定超弹性变形性能的形状记忆合金丝在拉伸到不同应变幅值条件下卸载的超弹性变形行为。根据试验测得的结果，提出了基于神经网络的形状记忆合金超弹性本构关系模型，并把模型计算的结果和实验数据进行了比较分析，结果表明，该模型具有很高的精度。该模型避免了已有模型在参数确定上的困难，具有一定的工程应用价值，为建立形状记忆合金本构模型提供了一个新的思路。     

Combining differential quadrature method with simulation technique to solve non‐linear differential equations

Nowadays, most of the ordinary differential equations (ODEs) can be solved by modelica‐based approaches, such as Matlab/Simulink, Dymola and LabView, which use simulation technique (ST). However, these kinds of approaches restrict the users in the enforcement of conditions at any instant of the time domain. This limitation is one of the most important drawbacks of the ST. Another method of solution, differential quadrature method (DQM), leads to very accurate results using only a few grids on the domain. On the other hand, DQM is not flexible for the solution of non‐linear ODEs and it is not so easy to impose multiple conditions on the same location. For these reasons, the author aims to eliminate the mentioned disadvantages of the simulation technique (ST) and DQM using favorable characteristics of each method in the other. This work aims to show how the combining method (CM) works simply by solving some non‐linear problems and how the CM gives more accurate results compared with those of other methods. Copyright © 2008 John Wiley & Sons, Ltd.     

The effect of grain orientation on the tensile‐compressive asymmetry of polycrystalline NiTi shape memory alloy

The purpose of the present study is to thoroughly understand the influence of crystallographic texture on the stress‐strain asymmetric behavior of polycrystalline NiTi shape memory alloy under tension and compression. To do this, a 3D thermo‐mechanical model has been implemented in a finite element program and textured and untextured polycrystalline NiTi have been considered. In our polycrystalline finite element model, each element represents one grain and a set of crystal orientations which approximate the initial crystallographic texture of the NiTi are assigned to the elements. From the calculated results, it is found that the crystallographic texture is the important reason for the tension‐compression asymmetry. For the textured polycrystal, the tension‐compression asymmetry can be observed clearly, but for the polycrystal containing randomly oriented grains, the stress‐strain curves show low levers of asymmetry between tensile and compressive loading, and the evolutions of martensite volume fractions are similar under two stress states.     

An efficient second‐order characteristic finite element method for non‐linear aerosol dynamic equations

In the paper we consider the non‐linear aerosol dynamic equation on time and particle size, which contains the advection process of condensation growth and the process of non‐linear coagulation. We develop an efficient second‐order characteristic finite element method for solving the problem. A high accurate characteristic method is proposed to treat the condensation advection while a second‐order extrapolation along the characteristics is proposed to approximate the non‐linear coagulation. The method has second‐order accuracy in time and the optimal‐order accuracy of finite element spaces in particle size, which improves the first‐order accuracy in time of the classical characteristic method. Numerical experiments show the efficient performance of our method for problems of log‐normal distribution aerosols in both the Euler coordinates and the logarithmic coordinates. Copyright © 2009 John Wiley & Sons, Ltd.     

A boundary condition in Padé series for frequency‐domain solution of wave propagation in unbounded domains

A boundary condition satisfying the radiation condition at infinity is frequently required in the numerical simulation of wave propagation in an unbounded domain. In a frequency domain analysis using finite elements, this boundary condition can be represented by the dynamic stiffness matrix of the unbounded domain defined on its boundary. A method for determining a Padé series of the dynamic stiffness matrix is proposed in this paper. This method starts from the scaled boundary finite‐element equation, which is a system of ordinary differential equations obtained by discretizing the boundary only. The coefficients of the Padé series are obtained directly from the ordinary differential equations, which are not actually solved for the dynamic stiffness matrix. The high rate of convergence of the Padé series with increasing order is demonstrated numerically. This technique is applicable to scalar waves and elastic vector waves propagating in anisotropic unbounded domains of irregular geometry. It can be combined seamlessly with standard finite elements. Copyright © 2006 John Wiley & Sons, Ltd.     

Be、Ni元素对Cu-Al系形状记忆合金性能的影响研究

利用Bridgemen定向凝固炉制备了不同成分的Cu-Al系形状记忆合金单晶.并对其形状记忆特性进行了系统分析.结果表明,Be元素与Ni元素同时添加的Cu-Al-Ni-Be形状记忆合金单晶具有比单独添加Be元素或Ni元素的合金单晶更优越的性能.Ni与Be元素的同时加入导致固溶淬火组织的变化,对改善合金形状记忆性能起了交互加强的作用.     

A generalized Ritz method for partial differential equations in domains of arbitrary geometry using global shape functions

A boundary element method (BEM)-based variational method is presented for the solution of elliptic PDEs describing the mechanical response of general inhomogeneous anisotropic bodies of arbitrary geometry. The equations, which in general have variable coefficients, may be linear or nonlinear. Using the concept of the analog equation of Katsikadelis the original equation is converted into a linear membrane (Poisson) or a linear plate (biharmonic) equation, depending on the order of the PDE under a fictitious load, which is approximated with radial basis function series of multiquadric (MQ) type. The integral representation of the solution of the substitute equation yields shape functions, which are global and satisfy both essential and natural boundary conditions, hence the name generalized Ritz method. The minimization of the functional that produces the PDE as the associated Euler–Lagrange equation yields not only the Ritz coefficients but also permits the evaluation of optimal values for the shape parameters of the MQs as well as optimal position of their centers, minimizing thus the error. If a functional does not exists or cannot be constructed as it is the usual case of nonlinear PDEs, the Galerkin principle can be applied. Since the arising domain integrals are converted into boundary line integrals, the method is boundary-only and, therefore, it maintains all the advantages of the pure BEM. Example problems are studied, which illustrate the method and demonstrate its efficiency and great accuracy.     

A new hybrid method for solving linear–non‐linear systems of equations

This paper presents a new method for solving any combination of linear–non‐linear equations. The method is based on the separation of linear equations in terms of some selected variables from the non‐linear ones. The linear group is solved by means of any method suitable for the linear system. This operation needs no iteration. The non‐linear group, however, is solved by an iteration technique based on a new formula using the Taylor series expansion. The method has been described and demonstrated in several examples of analytical systems with very good results. The new method needs the initial approximations for non‐linear variables only. This requires far less computation than the Newton–Raphson method. The method also has a very good convergence rate. The proposed method is most beneficial for engineering systems that very often involve a large number of linear equations with limited number of non‐linear equations. Copyright © 2005 John Wiley & Sons, Ltd.     

A finite element method for light activated shape‐memory polymers

Shape‐memory polymers (SMPs) belong to a class of smart materials that have shown promise for a wide range of applications. They are characterized by their ability to maintain a temporary deformed shape and return to an original parent permanent shape. In this paper, we consider the coupled photomechanical behavior of light activated shape‐memory polymers (LASMPs), focusing on the numerical aspects for finite element simulations at the engineering scale. The photomechanical continuum framework is summarized, and some specific constitutive equations for LASMPs are described. Numerical implementation of the multiphysics governing partial differential equations takes the form of a user defined element subroutine within the commercial software package ABAQUS . We verify our two‐dimensional and three‐dimensional finite element procedure for multiple analytically tractable cases. To show the robustness of the numerical implementation, simulations are performed under various geometries and complex photomechanical loading. Copyright © 2016 John Wiley & Sons, Ltd.     

NiTi shape memory alloys coated with calcium phosphate by plasma‐spraying. Chemical and biological properties

Plates of superelastic nickel‐titanium shape memory alloy (NiTi) were coated with calcium phosphate (hydroxyapatite) by high‐temperature plasma‐spraying. The porous layer of about 100 μm thickness showed a good adhesion to the metallic substrate that withstood bending of the plate but detached upon cutting the plate. The biocompatibility was tested by cultivation of blood cells (whole blood and isolated granulocytes [a subpopulation of blood leukocytes]). As substrates, pure NiTi, plasma‐spray‐coated NiTi and calcium phosphate‐coated NiTi prepared by a dip‐coating process were used. The adhesion of whole blood cells to all materials was not significantly different. In contrast, isolated granulocytes showed an increased adhesion to both calcium phosphate‐coated NiTi samples. However, compared to non‐coated NiTi or dip‐coated NiTi, the number of dead granulocytes adherent to plasma‐sprayed surfaces was significantly increased for isolated granulocytes (p<0.01).     

A simulation study on the thermal buckling behavior of laminated composite shells with embedded shape memory alloy (SMA) wires

In this paper, numerical simulation analyses of the thermal buckling behavior of laminated composite shells with embedded shape memory alloy (SMA) wires were performed to investigate the effect of embedded SMA wires on the characteristics of thermal buckling. In order to simulate the thermomechanical behavior of SMA wires, the constitutive equation of the SMA wires was formulated in the form of an ABAQUS user subroutine. The computational program was verified by showing the response of the pseudoelasticity and shape memory effect (SME) at various temperatures and stress levels. Modeling of the laminated composite shells with embedded SMA wires and thermal buckling analyses were performed with the use of the ABAQUS code linked with the subroutine of the formulated SMA constitutive equations. The thermal buckling analyses of the composite shells with embedded SMA wires show that the critical buckling temperature can be increased and the thermal buckling deformation can be decreased by using the activation force of embedded SMA wire actuators.