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.     

基于非局部塑性模型的应变局部化理论分析及数值模拟

通过求解一个第二类Fredholm方程,得到了基于非局部塑性软化模型的应变局部化问题理论解,结果表明,只有在当采用过非局部修正形式的非局部塑性软化模型才能得到应变局部化解,且得到的塑性应变分布和荷载响应依赖于所引入的特征长度及过非局部权参数。通过一维应变局部化有限元数值解,验证了非局部理论的引入能克服计算结果的网格敏感...     

Hopf bifurcation in general Brusselator system with diffusion

The general Brusselator system is considered under homogeneous Neumann boundary conditions. The existence results of the Hopf bifurcation to the ordinary differential equation (ODE) and partial differential equation (PDE) models are obtained. By the center manifold theory and the normal form method, the bifurcation direction and stability of periodic solutions are established. Moreover, some numerical simulations are shown to support the analytical results. At the same time, the positive steady-state solutions and spatially inhomogeneous periodic solutions are graphically shown to supplement the analytical results.     

An analytical study on the geometrical size effect on phase transitions in a slender compressible hyperelastic cylinder

In this paper, we study phase transitions in a slender circular cylinder composed of a compressible hyperelastic material with a non-convex strain-energy function in a loading process. We aim to construct the asymptotic solutions based on an axisymmetrical three-dimensional setting and use the results to describe the key features observed in the experiments by others. By using a methodology involving coupled series-asymptotic expansions, we derive the normal form equation of the original complicated system of non-linear PDEs. Based on a phase-plane analysis, we manage to deduce the global bifurcation properties and to solve the boundary-value problem analytically. The explicit solutions (including post-bifurcation solutions) in terms of integrals are obtained. The engineering stress-strain curve plotted from the asymptotic solutions can capture the key features of the curve measured in some experiments. Our results can also describe the geometrical size effect as observed in experiments. It appears that the asymptotic solutions obtained shed certain light on the instability phenomena associated with phase transitions in a cylinder.     

Thermoelastic surface waves on an exponentially graded half-space

In this paper, we consider the propagation of Rayleigh surface waves in a functionally graded isotropic thermoelastic half-space, in which all thermoelastic characteristic parameters exponentially change along the depth direction. The propagation condition is established in the form of a bicubic equation whose coefficients are complex numbers while the analytical solutions (eigensolutions) of the thermoelastodynamic system are explicitly obtained in terms of the characteristic solutions. The concerned solution of the Rayleigh surface wave problem is subsequently expressed as a linear combination of the three eigensolutions while the secular equation is established in an implicit form. The explicit secular equation is written when an isotropic and homogeneous thermoelastic half-space is considered and some numerical simulations are given for a specific material.     

Green’s function solution for transient heat conduction in annular fin during solidification of phase change material

The Green’s function method is applied for the transient temperature of an annular fin when a phase change material (PCM) solidifies on it. The solidification of the PCMs takes place in a cylindrical shell storage. The thickness of the solid PCM on the fin varies with time and is obtained by the Megerlin method. The models are found with the Bessel equation to form an analytical solution. Three different kinds of boundary conditions are investigated. The comparison between analytical and numerical solutions is given. The results demonstrate that the significant accuracy is obtained for the temperature distribution for the fin in all cases.     

Chaotic Oscillation of Saddle Form Cable-Suspended Roofs under Vertical Excitation Action

The purpose of this paper is to investigate the dynamic behaviour of saddle form cable-suspended roofs under vertical excitation action. The governing equations of this problem are system of nonlinear partial differential and integral equations. We first establish a spectral equation, and then consider a model with one coefficient, i.e., a perturbed Duffing equation. The analytical solution is derived for the Duffing equation. Successive approximation solutions can be obtained in likely way for each time to only one new unknown function of time. Numerical results are given for our analytical solution. By using the Melnikov method, it is shown that the spectral system has chaotic solutions and subharmonic solutions under determined parametric conditions.     

弹性力学中一种新的边界轮廓法

利用基本解的特性,将面力积分方程化成仅含有Ｃａｕｃｈｙ主值积分的形式,基于这种边界积分方程,提出了一种新的边界轮廓法,对于三维问题,该方法只须计算沿边界单元界线的线积分,对二维问题,则只需计算边界单元两点的热函数之差,无须进行数值积分计算,实例计算说明该方法是有效的。     

Nonlinear seismic waves: A model for site effects

The aim of this paper is to propose a possible mathematical model of site effects that occur when seismic waves propagate through a sediment filled basin. The model is based on the mechanical properties of the medium (that we consider as a granular material) through which the seismic waves propagate. By looking for asymptotic solutions having the features of a progressive wave, we derive an evolution equation which is a modified Korteweg–deVries–Burgers equation containing also a nonlinear dissipative term. This equation is integrated numerically and the modelled site amplification is evaluated by using the smoothed spectral ratio between the propagated profile of the wave and the initial one.     

Stability analysis of functionally graded rectangular plates under nonlinearly varying in-plane loading resting on elastic foundation

In this research work, an exact analytical solution for buckling of functionally graded rectangular plates subjected to non-uniformly distributed in-plane loading acting on two opposite simply supported edges is developed. It is assumed that the plate rests on two-parameter elastic foundation and its material properties vary through the thickness of the plate as a power function. The neutral surface position for such plate is determined, and the classical plate theory based on exact neutral surface position is employed to derive the governing stability equations. Considering Levy-type solution, the buckling equation reduces to an ordinary differential equation with variable coefficients. An exact analytical solution is obtained for this equation in the form of power series using the method of Frobenius. By considering sufficient terms in power series, the critical buckling load of functionally graded plate with different boundary conditions is determined. The accuracy of presented results is verified by appropriate convergence study, and the results are checked with those available in related literature. Furthermore, the effects of power of functionally graded material, aspect ratio, foundation stiffness coefficients and in-plane loading configuration together with different combinations of boundary conditions on the critical buckling load of functionally graded rectangular thin plate are studied.     

Mode localization in lateral buckling of partially embedded submarine pipelines

The partially embedded submarine pipelines might buckle laterally at some segments under high pressure and high temperature (HP/HT) conditions. The buckling pattern localization introduces an extra level of analytical complexity when compared with the periodic buckling pattern. In the presented paper the lateral buckling pipeline is modeled by an axial compressive beam supported by lateral distributing nonlinear springs taking the soil berm effects in the horizontal plane. It is found that the model is governed by a time-independent Swift–Hohenberg equation. Based on John Burke and Edgar Knobloch’s work we conclude preliminarily that the equation will have localized solutions. Besides the qualitative conclusion, by AUTO 07P the localized solutions of the equation are studied in detail. The snakes-and-ladders structure of localized solutions explains the transition of buckling modes in theory. The range of the possible critical axial forces is found out. Meanwhile two critical axial force formulas corresponding to the range ends are presented. Finally a typical submarine pipeline is analyzed as an illustration.     

Anti-plane stress intensity, energy release and energy density at crack tips in a functionally graded strip with linearly varying properties

The fracture problem of a crack in a functionally graded strip with its properties varying in a linear form along the strip thickness under an anti-plane load is considered. The embedded anti-plane crack is located in the middle of strip half way through the thickness. The third mode stress intensity factor is derived using two different methods. In the first method, by employing Fourier integral transforms, the governing equation is converted to a singular integral equation, which is subsequently solved numerically by the collocation method based on Chebyshev polynomials. Then, the problem is solved by means of finite element method in which quadrilateral 8-node singular elements around each crack tip are used. After inspecting the validity of the solution technique, effects of crack geometry and non-homogeneous material parameter on the stress intensity, energy release and energy density are studied and the results of analytical and FEM solutions are compared.     

General boundary conditions for the wave equation around non-homogenous scatterers

To solve the wave equation inside a region that contains an inhomogenous dielectric material of arbitrary shape under the influence of an incoming wave, we establish a generalized boundary condition. The solutions inside a finite region resulting form a given incoming wave from the outside, are determined by a linear relation between the normal gradient and the function values on the boundary. This boundary condition is non-local and we show how it can be used in conjunction with the variational principle applied to an open system.     

The influence of normal flow boundary conditions on spurious modes in finite element solutions to the shallow water equations

Finite element solutions of the primitive equation (PE) form of the shallow water equations are notorious for the severe spurious 2Δx modes which appear. Wave equation (WE) solutions do not exhibit these numerical modes. In this paper we show that the severe spurious modes in PE solutions are strongly influenced by essential normal flow boundary conditions in the coupled continuity-momentum system of equations. This is demonstrated through numerical examples that avoid the use of essential normal flow boundary conditions either by specifying elevation values over the entire boundary or by implementing natural flow boundary conditions in the weak weighted residual form of the continuity equation. Results from a series of convergence tests show that PE solutions are of nearly the same quality as WE solutions when spurious modes are suppressed by alternative specification of the boundary conditions. Network intercomparisons indicate that varying nodal support does not excite spurious modes in a solution, although it does enhance the spurious modes introduced when an essential normal flow boundary condition is used. Dispersion analysis of discrete equations for interior and boundary nodes offers an explanation of the observed solution behaviour. For certain PE algorithms a mixed situation can arise where the boundary nodes exhibit a monotonic (noise-free) dispersion relationship and the interior nodes exhibit a folded (noisy) dispersion relationship. We have found that the mixed situation occurs when all boundary nodes are specified elevation nodes (which are enforced as essential conditions in the continuity equation) or when specified flow boundary nodes are treated as natural boundary conditions in the continuity equation. In either case the effect is to generate a solution that is essentially free of noise. Apparently, the monotonic dispersion behaviour at the boundaries suppresses the otherwise noisy behaviour caused by the folded dispersion relation on the interior.     

Numerical solution of composite left and right fractional Caputo derivative models for granular heat flow

In this paper we propose a numerical scheme based on a fractional trapezoidal method for solution of a fractional equation with composition of the left and right Caputo derivatives. The numerical results are compared with analytical solutions. We have illustrated the convergence of our scheme. Finally, we show an application of the considered equation.     

Investigation of the size effects in Timoshenko beams based on the couple stress theory

In this paper, a size-dependent Timoshenko beam is developed on the basis of the couple stress theory. The couple stress theory is a non-classic continuum theory capable of capturing the small-scale size effects on the mechanical behavior of structures, while the classical continuum theory is unable to predict the mechanical behavior accurately when the characteristic size of structures is close to the material length scale parameter. The governing differential equations of motion are derived for the couple-stress Timoshenko beam using the principles of linear and angular momentum. Then, the general form of boundary conditions and generally valid closed-form analytical solutions are obtained for the axial deformation, bending deflection, and the rotation angle of cross sections in the static cases. As an example, the closed-form analytical results are obtained for the response of a cantilever beam subjected to a static loading with a concentrated force at its free end. The results indicate that modeling on the basis of the couple stress theory causes more stiffness than modeling by the classical beam theory. In addition, the results indicate that the differences between the results of the proposed model and those based on the classical Euler–Bernoulli and classical Timoshenko beam theories are significant when the beam thickness is comparable to its material length scale parameter.     

Theoretical analysis,numerical calculation and experimental measurement of transient elastic waves in strips subjected to dynamic loadings

In this study, the transient response of an elastic strip subjected to dynamic in-plane loadings on the surface is investigated in detail. One of the objectives of this study is to develop an effective analytical method for determining transient solutions in a strip. By applying Laplace transform, the analytical solution in the transformed domain is derived and expressed in matrix form. The solution is then decomposed into infinite wave groups in which the multiple reflected waves with the same reflection are involved. Each multi-reflected wave can be identified by a coding method and be verified by the theory of generalized ray. The inverse transform is performed by using the well-known Cagniard method. The transient solutions in time domain for stresses and displacements are expressed in a closed form and are discussed in detail by an example. The experimental results show that the early time transient responses of displacements on the surface agree very well with the numerical calculations based on the theoretical solutions.     

Phase transitions induced by extension in a slender SMA cylinder: Analytical solutions for the hysteresis loop based on a quasi-3D continuum model

In this paper, we propose a quasi-3D continuum model to study the rate-independent hysteresis phenomenon in phase transitions of a slender shape memory alloy (SMA) cylinder subject to the uniaxial tension. Based on the three-dimensional field equations and the traction-free boundary conditions, by using a coupled series-asymptotic expansion method, we manage to express the total elastic potential energy of the cylinder in terms of the leading order term of the axial strain. We further consider the rate-independent dissipation effect in a purely one-dimensional setting. The mechanical dissipation functions are also expressed in terms of the axial strain. The equilibrium configuration of the cylinder is then determined by using the principle of maximizing the total energy dissipation. An illustrative example with some special chosen material constants is further considered. Free end boundary conditions are proposed at the two ends of the cylinder. By conducting a phase plane analysis and through some calculations, we obtain the analytical solutions of the equilibrium equation. We find that the engineering stress–strain curves corresponding to the obtained solutions can capture some important features of the experimental results. It appears that the analytical results obtained in this paper reveal the multiple solutions nature of the problem and shed certain light on the instability phenomena during the phase transition process.