Post-buckling and nonlinear free vibration analysis of geometrically imperfect functionally graded beams resting on nonlinear elastic foundation

In this paper, post-buckling and nonlinear vibration analysis of geometrically imperfect beams made of functionally graded materials (FGMs) resting on nonlinear elastic foundation subjected to axial force are studied. The material properties of FGMs are assumed to be graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents. The assumptions of a small strain and moderate deformation are used. Based on Euler–Bernoulli beam theory and von-Karman geometric nonlinearity, the integral partial differential equation of motion is derived. Then this partial differential equation (PDE) problem, which has quadratic and cubic nonlinearities, is simplified into an ordinary differential equation (ODE) problem by using the Galerkin method. Finally, the governing equation is solved analytically using the variational iteration method (VIM). Some new results for the nonlinear natural frequencies and buckling load of the imperfect functionally graded (FG) beams such as the effects of vibration amplitude, elastic coefficients of foundation, axial force, end supports and material inhomogeneity are presented for future references. Results show that the imperfection has a significant effect on the post-buckling and vibration response of FG beams.     

ERROR ESTIMATES FOR SPARSE OPTIMAL CONTROL PROBLEMS BY PIECEWISE LINEAR FINITE ELEMENT APPROXIMATION

Optimization problems with L¹-control cost functional subject to an elliptic partial differential equation(PDE)are considered.However,different from the finite dimensiona l¹-regularization optimization,the resulting discretized L¹norm does not have a decoupled form when the standard piecewise linear finite element is employed to discretize the continuous problem.A common approach to overcome this difficulty is employing a nodal quadrature formula to approximately discretize the L¹-norm.In this paper,a new discretized scheme for the L¹-norm is presented.Compared to the new discretized scheme for L¹-norm with the nodal quadrature formula,the advantages of our new discretized scheme can be demonstrated in terms of the order of approximation.Moreover,finite element error estimates results for the primal problem with the new discretized scheme for the L¹-norm are provided,which confirms that this approximation scheme will not change the order of error estimates.To solve the new discretized problem,a symmetric Gauss-Seidel based majorized accelerated block coordinate descent(sGS-mABCD)method is introduced to solve it via its dual.The proposed sGS-mABCD algorithm is illustrated at two numerical examples.Numerical results not only confirm the finite element error estimates,but also show that our proposed algorithm is efficient.     

A Variational Analysis for the Moving Finite Element Method for Gradient Flows

By using the Onsager principle as an approximation tool, we give a novel derivation for the moving finite element method for gradient flow equations. We show that the discretized problem has the same energy dissipation structure as the continuous one. This enables us to do numerical analysis for the stationary solution of a nonlinear reaction diffusion equation using the approximation theory of free-knot piecewise polynomials. We show that under certain conditions the solution obtained by the moving finite element method converges to a local minimizer of the total energy when time goes to infinity. The global minimizer, once it is detected by the discrete scheme, approximates the continuous stationary solution in optimal order. Numerical examples for a linear diffusion equation and a nonlinear Allen-Cahn equation are given to verify the analytical results.     

On methods for continuous systems with quadratic,cubic and quantic nonlinearities

Methods for study of weakly nonlinear continuous systems are discussed. The method of multiple scales is used to analyze the nonlinear response of a relief valve under combined static and dynamic loadings. We determine a second-order approximation to the response of the system for the case of primary resonance. Second, we derive a second-order nonlinear ordinary differential equation that describes the time evolution of a single-mode, the so-called single-mode discretization. Then, we use the multiple scales method to determine second-order approximate solutions of this equation, thereby obtaining the equations describe the modulations of the amplitude and phase of the response. We show that the results of the second approach are erroneous.     

Validity and accuracy of solutions for nonlinear vibration analyses of suspended cables with one-to-one internal resonance

This study investigates the accuracy of nonlinear vibration analyses of a suspended cable, which possesses quadratic and cubic nonlinearities, with one-to-one internal resonance. To this end, we derive approximate solutions for primary resonance using two different approaches. In the first approach, the method of multiple scales is directly applied to governing equations, which are nonlinear partial differential equations. In the second approach, we first discretize the governing equations by using Galerkin’s procedure and then apply the shooting method. The accuracy of the results obtained by these approaches is confirmed by comparing them with results obtained by the finite difference method.     

多自由度强非线性颤振分析的增量谐波平衡法

对多个自由度上含有强非线性项系统的颤振问题,推广应用增量谐波平衡法进行分析．考虑带有强非线性立方平移和俯仰刚度项的二元机翼颤振方程,首先将方程用矩阵形式表示,然后把振动过程分解成为振动瞬态的持续增量过程,再采用振幅作为控制参数应用谐波平衡法,以这种推广的增量谐波平衡法求得方程解的表达式,并由此分析系统的分岔现象、极限环颤振现象和谐波项数的取值问题,最后用龙格-库塔数值方法进行验算,结果表明:分析多个自由度的强非线性颤振,增量谐波平衡法是精确有效的．     

Semi-analytical and numerical solutions of multi-degree-of-freedom nonlinear oscillation systems with linear coupling

This research focuses on the development of an approach for solving multi-degree-of-freedom (MDOF) nonlinear oscillation problems with linear coupling. The original physical information included in the governing equations is mostly transferred into semi-analytical and numerical solutions. The semi-analytical solutions generated by the present approach are continuous everywhere and reflect more accurately the characteristics of the motion of the nonlinear dynamic systems. General procedures for three types of nonlinear oscillation problems are formulated in detail for allocation in nonlinear dynamic analysis. Two nonlinear oscillation systems with quadratic and cubic nonlinearities are solved to demonstrate the applications of the present approach.     

Iterative method for finding the eigenfunctions of a system of two Schrödinger equations with combined nonlinearity

An iterative method is proposed to determine the eigenfunctions of a system of two nonlinear Schrödinger equations governing the interaction of two femtosecond pulses in a medium with quadratic and cubic nonlinearity. The method produces soliton solutions of a new form for a wide range of nonlinearity coefficients corresponding to the first and second eigenvalues. A specially chosen initial approximation is required to determine the third and higher eigenfunctions.     

Analytical solution for general nonlinear continuous systems in a complex form

In this paper, the method of multiple scales is used to study free vibrations and primary resonances of geometrically nonlinear spatial continuous systems with general quadratic and cubic nonlinear operators in a complex form. It is found that in the free vibrations of general continuous systems in a complex form, both forward and backward modes are excited. This situation is in contrast to the primary resonances in which only forward modes are excited. Consequently, one may determine the form of solution before applying the multiple scales method to the equation. This analysis is applicable to general continuous systems with gyroscopic and Coriolis effects and includes many nonlinear problems as a special case. As an example of application of this general solution, free vibrations and primary resonances of a simply supported rotating shaft with stretching nonlinearity are considered.     

非线性阻尼作用下标准线性固体粘弹性Ⅲ型破裂的解析解

把非线性Rayleigh阻尼引入标准线性固体粘弹性介质的Ⅲ型破裂的控制方程中,此方程是一个偏微分积分方程;首先设法消去积分项,得到一个三阶非线性偏微分方程,然后用小参数摄动法,得出线性化的各阶渐近控制方程;把每一个具有变系数的三阶线性控制方程分解为弹性部分及剩余部份,而前者的解析解为已知,后者是一个二阶变系数线性偏微分方程;它化不成Mathieu方程,也化不成Hill方程,故采用WKBJ的方法得出其渐近的解析解。     

A new perturbation solution for systems with strong quadratic and cubic nonlinearities

The new perturbation algorithm combining the method of multiple scales (MS) and Lindstedt–Poincare techniques is applied to an equation with quadratic and cubic nonlinearities. Approximate analytical solutions are found using the classical MS method and the new method. Both solutions are contrasted with the direct numerical solutions of the original equation. For the case of strong nonlinearities, solutions of the new method are in good agreement with the numerical results, whereas the amplitude and frequency estimations of classical MS yield high errors. For strongly nonlinear systems, exact periods match well with the new technique while there are large discrepancies between the exact and classical MS periods. Copyright © 2009 John Wiley & Sons, Ltd.     

Global Well-Posedness for Cubic NLS with Nonlinear Damping

We study the Cauchy problem for the cubic nonlinear Schrödinger equation, perturbed by (higher order) dissipative nonlinearities. We prove global in-time existence of solutions for general initial data in the energy space. In particular we treat the energy-critical case of a quintic dissipation in three space dimensions.     

A numerical method for the Cahn–Hilliard equation with a variable mobility

We consider a conservative nonlinear multigrid method for the Cahn–Hilliard equation with a variable mobility of a model for phase separation in a binary mixture. The method uses the standard finite difference approximation in spatial discretization and the Crank–Nicholson semi-implicit scheme in temporal discretization. And the resulting discretized equations are solved by an efficient nonlinear multigrid method. The continuous problem has the conservation of mass and the decrease of the total energy. It is proved that these properties hold for the discrete problem. Also, we show the proposed scheme has a second-order convergence in space and time numerically. For numerical experiments, we investigate the effects of a variable mobility.     

A Power Penalty Approach to Numerical Solutions of Two-Asset American Options

This paper aims to develop a power penalty method for a linear parabolic variational inequality (Ⅵ) in two spatial dimensions governing the two-asset Ameri-can option valuation. This method yields a two-dimensional nonlinear parabolic PDE containing a power penalty term with penalty constant λ>1 and a power parameter k>0. We show that the nonlinear PDE is uniquely solvable and the solution of the PDE converges to that of the VI at the rate of order (λ<-k/2>). A fitted finite volume method is designed to solve the nonlinear PDE, and some numerical experiments are performed to illustrate the usefulness of this method.     

The Singular Function Boundary Integral Method for singular Laplacian problems over circular sections

The Singular Function Boundary Integral Method (SFBIM) for solving two-dimensional elliptic problems with boundary singularities is revisited. In this method the solution is approximated by the leading terms of the asymptotic expansion of the local solution, which are also used to weight the governing partial differential equation. The singular coefficients, i.e., the coefficients of the local asymptotic expansion, are thus primary unknowns. By means of the divergence theorem, the discretized equations are reduced to boundary integrals and integration is needed only far from the singularity. The Dirichlet boundary conditions are then weakly enforced by means of Lagrange multipliers, the discrete values of which are additional unknowns. In the case of two-dimensional Laplacian problems, the SFBIM converges exponentially with respect to the numbers of singular functions and Lagrange multipliers. In the present work the method is applied to Laplacian test problems over circular sectors, the analytical solution of which is known. The convergence of the method is studied for various values of the order p of the polynomial approximation of the Lagrange multipliers (i.e., constant, linear, quadratic, and cubic), and the exact approximation errors are calculated. These are compared to the theoretical results provided in the literature and their agreement is demonstrated.     

Third-order iterative methods free from second derivatives for nonlinear equation

In this paper, we suggest and analyze a new two-step iterative method for solving nonlinear equations, which is called the modified Householder method without second derivatives for nonlinear equation. We also prove that the modified method has cubic convergence. Several examples are given to illustrate the efficiency and the performance of the new method. New method can be considered as an alternative to the present cubic convergent methods for solving nonlinear equations.     

Approximate solution to the time–space fractional cubic nonlinear Schrodinger equation

By introducing the fractional derivatives in the sense of Caputo, we use the adomian decomposition method to construct the approximate solutions for the cubic nonlinear fractional Schordinger equation with time and space fractional derivatives. The exact solution of the cubic nonlinear Schrodinger equation is given as a special case of our approximate solution. This method is efficient and powerful in solving wide classes of nonlinear evolution fractional order equation.     

On weakly nonlinear waves in media exhibiting mixed nonlinearity

We obtain the transport equations governing small amplitude high frequency disturbances, that include both quadratic and cubic nonlinearities inherent in hyperbolic systems of conservation laws. The coefficients of the nonlinear terms in the transport equation are obtained in terms of the Glimm interaction coefficients. For symmetric and isotropic systems the mean curvature of the wave front, which appears as the coefficient of the linear term in the transport equation, is shown to be related to the derivative of the ray tube area along the bicharacteristics; the amplitude of the disturbance is shown to become unbounded in the neighborhood of the point where the ray tube collapses. We also obtain a formula, akin to the one obtained by R. Rosales (1991), for the energy dissipated across shocks.     

Nonlinear asymptotic analysis of a system of two free coupled oscillators with cubic nonlinearities

In this paper a two degrees of freedom undamped nonlinear system of two unforced coupled oscillators with cubic nonlinearities is analyzed. Through a decoupling procedure and using admissible functional transformations it is proved that this system can be reduced to an intermediate second order nonlinear ordinary differential equation (ODE) connecting both displacements to each other. By nonlinear asymptotic approximations the above equation can be further reduced to new nonlinear ODE that can be analytically solved. The solutions in the physical plane are extracted in parametric form. As generalization, the model of a damped system of two masses connected with stiffness with linear and nonlinear coefficient of rigidities respectively is analyzed and exact analytical solutions are extracted. Finally an application has been given in the case of a two mass system with cubic strong non-linearity.     

解非线性方程组的一类离散的Newton算法

1．引言考虑非线性方程组设xi是当前的迭代点，为计算下一个迭代点，Newton法是求解方程若用差商代替导数，离散Newton法要解如下的方程其中这里为了计算J（；；h），需计算n‘个函数值．为了提高效能，Brown方法l‘］使用代入消元的办法来减少函数值计算量．它是再通过一次内选代从h得到下一个迭代点14＋1．设n；＝（《1，…，Zn尸，t二（ti，…，t＊”，t为变量．BfOWll方法的基本思想如下．对人（x）在X；处做线性近似解出然后代入第二个函数，得到这是关于tZ，…，tn的函数．当（tZ，…，t。尸一（ZZ，…，Z。厂时，由（1．4），…