Robustness of mechanical systems against uncertainties

In this paper, one can propose a method which takes into account the propagation of uncertainties in the finite element models in a multi-objective optimization procedure. This method is based on the coupling of stochastic response surface method (SRSM) and a genetic algorithm provided with a new robustness criterion. The SRSM is based on the use of stochastic finite element method (SFEM) via the use of the polynomial chaos expansion (PC). Thus, one can avoid the use of Monte Carlo simulation (MCS) whose costs become prohibitive in the optimization problems, especially when the finite element models are large and have a considerable number of design parameters.The objective of this study is on one hand to quantify efficiently the effects of these uncertainties on the responses variability or the cost functions which one wishes to optimize and on the other hand, to calculate solutions which are both optimal and robust with respect to the uncertainties of design parameters.In order to study the propagation of input uncertainties on the mechanical structure responses and the robust multi-objective optimization with respect to these uncertainty, two numerical examples were simulated. The results which relate to the quantification of the uncertainty effects on the responses variability were compared with those obtained by the reference method (REF) using MCS and with those of the deterministic response surfaces methodology (RSM).In the same way, the robust multi-objective optimization results resulting from the SRSM method were compared with those obtained by the direct optimization considered as reference (REF) and with RSM methodology.The SRSM method application to the response variability study and the robust multi-objective optimization gave convincing results.     

On the accelerated modified Newton-HSS method for systems of nonlinear equations

Hermitian and skew-Hermitian splitting (HSS) method converges unconditionally, which is efficient and robust for solving non-Hermitian positive-definite systems of linear equations. For solving systems of nonlinear equations with non-Hermitian positive-definite Jacobian matrices, Bai and Guo proposed the Newton-HSS method and gave numerical comparisons to show that the Newton-HSS method is superior to the Newton-USOR, the Newton-GMRES and the Newton-GCG methods. Recently, Wu and Chen proposed the modified Newton-HSS (MN-HSS) method which outperformed the Newton-HSS method. In this paper, we will establish a new accelerated modified Newton-HSS (AMN-HSS) method and give the local convergence theorem. Moreover, numerical results show that the AMN-HSS method outperforms the MN-HSS method.     

Parameterized defuzzification with maximum entropy weighting function—Another view of the weighting function expectation method

This paper mainly extends the weighting function expectation method from fuzzy numbers to the ordinary fuzzy set, and regards it as a weighting function defuzzification method. Some properties are discussed in the defuzzification context. The traditional center of gravity (COG) defuzzification method is extended to the weighting function COG (WFCOG) defuzzification method, which can represent the decision maker’s subjective attitude. With the properties of the weighting function, a parameterized defuzzification method with maximum entropy weighting function is proposed. Then it is integrated with the basic defuzzification distribution (BADD) method, which can represent the decision maker’s valuation and confidential attitude respectively. The relationships between the weighted function BADD defuzzification method and some existing defuzzification methods are discussed. The COG, MOM (mean of maxima), FOM (first of maxima), LOM (last of maxima), MOS (mean of support), FOS (first of support), and LOS (last of support) all become special cases with specific valuation and confidential attitude parameters.     

关于“评测指标权重确定的结构熵权法”的注记

本文基于文献[1]所提出的观点,根据熵理论和德尔斐专家调查法,对文献[1]所提出的方法进行了改进,得到了改进的结构熵权法,并运用此方法对环保项目指标进行权重确定,同时与经典的德尔斐专家调查法确定的权重进行比较,并对二者的差异作出了合理解释。     

A Meshless Boundary Method for Computation of Stokes Flow with Particles in the Semicircular Canals of the Inner Ear

Benign paroxismal positional vertigo (BPPV) is modeled by introducing free-floating particles (canaliths) which settle inside the semicircular canals (SCC). The Stokes flow induced by a canalith is evaluated by coupling the force coupling method (FCM) to the method of fundamental solutions (MFS). The proposed methodology results in a straightforward meshless boundary method for the simulation of bounded Stokes flow with finite-size particles. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A numerical solution of the dirichlet problem on some special doubly connected regions

The aim of this paper is to give a convergence proof of a numerical method for the Dirichlet problem on doubly connected plane regions using the method of reflection across the exterior boundary curve (which is analytic) combined with integral equations extended over the interior boundary curve (which may be irregular with infinitely many angular points).     

On Global Convergence and Approximate Iteration of the Linear Approximation Method for Solving Variational Inequalities

This paper is concerned with the linear approximation method (i.e. the iterative method in which a sequence of vectors is generated by solving certain linearized subproblems) for solving the variational inequality. The global convergent iterative process is proposed by applying the continuation method, and the related problems are discussed. A convergent result is obtained for the approximation iteration (i.e. the iterative method in which a sequence of vectors is generated by solving certain linearized subproblems approximately).     

Numerical solution of Burgers’ equation by the Sobolev gradient method

Burgers’ equation is solved numerically with Sobolev gradient methods. A comparison is shown with other numerical schemes presented in this journal, such as modified Adomian method (MAM) [1] and by a variational method (VM) which is based on the method of discretization in time [2]. It is shown that the Sobolev gradient methods are highly efficient while at the same time retaining the simplicity of steepest descent algorithms.     

A random integral quadrature method for numerical analysis of the second kind of Volterra integral equations

In this paper, a novel meshless technique termed the random integral quadrature (RIQ) method is developed for the numerical solution of the second kind of the Volterra integral equations. The RIQ method is based on the generalized integral quadrature (GIQ) technique, and associated with the Kriging interpolation function, such that it is regarded as an extension of the GIQ technique. In the GIQ method, the regular computational domain is required, in which the field nodes are scattered along straight lines. In the RIQ method however, the field nodes can be distributed either uniformly or randomly. This is achieved by discretizing the governing integral equation with the GIQ method over a set of virtual nodes that lies along straight lines, and then interpolating the function values at the virtual nodes over all the field nodes which are scattered either randomly or uniformly. In such a way, the governing integral equation is converted approximately into a system of linear algebraic equations, which can be easily solved.     

An extension and analysis of the Shu-Osher representation of Runge-Kutta methods

In the context of solving nonlinear partial differential equations, Shu and Osher introduced representations of explicit Runge-Kutta methods, which lead to stepsize conditions under which the numerical process is total-variation-diminishing (TVD). Much attention has been paid to these representations in the literature.

In general, a Shu-Osher representation of a given Runge-Kutta method is not unique. Therefore, of special importance are representations of a given method which are best possible with regard to the stepsize condition that can be derived from them.

Several basic questions are still open, notably regarding the following issues: (1) the formulation of a simple and general strategy for finding a best possible Shu-Osher representation for any given Runge-Kutta method; (2) the question of whether the TVD property of a given Runge-Kutta method can still be guaranteed when the stepsize condition, corresponding to a best possible Shu-Osher representation of the method, is violated; (3) the generalization of the Shu-Osher approach to general (possibly implicit) Runge-Kutta methods.

In this paper we give an extension and analysis of the original Shu-Osher representation, by means of which the above questions can be settled. Moreover, we clarify analogous questions regarding properties which are referred to, in the literature, by the terms monotonicity and strong-stability-preserving (SSP).

     

A fast numerical solution of the 3D nonlinear tempered fractional integrodifferential equation

In this paper, we investigate the numerical solution of the three-dimensional (3D) nonlinear tempered fractional integrodifferential equation which is subject to the initial and boundary conditions. The backward Euler (BE) method in association with the first-order convolution quadrature rule is employed to discretize this equation for time, and the Galerkin finite element method is applied for space, which is combined with an alternating direction implicit (ADI) algorithm, in order to reduce the computational cost for solving the three-dimensional nonlocal problem. Then a fully discrete BE ADI Galerkin finite element scheme can be obtained by linearizing the non-linear term. Thereafter we prove a positive-type lemma, from which the stability and convergence of the proposed numerical scheme are derived based on the energy method. Numerical experiments are performed to verify the effectiveness of the proposed approach.     

Convergence of a generalized MSSOR method for augmented systems

Recently, Wu et al. [S.-L. Wu, T.-Z. Huang, X.-L. Zhao, A modified SSOR iterative method for augmented systems, J. Comput. Appl. Math. 228 (1) (2009) 424-433] introduced a modified SSOR (MSSOR) method for augmented systems. In this paper, we establish a generalized MSSOR (GMSSOR) method for solving the large sparse augmented systems of linear equations, which is the extension of the MSSOR method. Furthermore, the convergence of the GMSSOR method for augmented systems is analyzed and numerical experiments are carried out, which show that the GMSSOR method with appropriate parameters has a faster convergence rate than the MSSOR method with optimal parameters.     

Backward differentiation type formulas for Volterra integral equations of the second kind

Summary Numerical integration formulas are discussed which are obtained by differentiation of the Volterra integral equation and by applying backward differentiation formulas to the resulting integro-differential equation. In particular, the stability of the method is investigated for a class of convolution kernels. The accuracy and stability behaviour of the method proposed in this paper is compared with that of (i) a block-implicit Runge-Kutta scheme, and (ii) the scheme obtained by applying directly a quadrature rule which is reducible to the backward differentiation formulas. The present method is particularly advantageous in the case of stiff Volterra integral equations.     

Dual reciprocity hybrid radial boundary node method for the analysis of Kirchhoff plates

A meshless method of dual reciprocity hybrid radial boundary node method (DHRBNM) for the analysis of arbitrary Kirchhoff plates is presented, which combines the advantageous properties of meshless method, radial point interpolation method (RPIM) and BEM. The solution in present method comprises two parts, i.e., the complementary solution and the particular solution. The complementary solution is solved by hybrid radial boundary node method (HRBNM), in which a three-field interpolation scheme is employed, and the boundary variables are approximated by RPIM, which is applied instead of moving least square (MLS) and obtains the Kronecker’s delta property where the traditional HBNM does not satisfy. The internal variables are interpolated by two groups of symmetric fundamental solutions. Based on those, a hybrid displacement variational principle for Kirchhoff plates is developed, and a meshless method of HRBNM for solving biharmonic problems is obtained, by which the complementary solution can be solved.     

The Laplace transform and polynomial Trefftz method for a class of time dependent PDEs

In this paper, we present a new method for solving 1D time dependent partial differential equations based on the Laplace transform (LT). As a result, the problem is converted into a stationary boundary value problem (BVP) which depends on the parameter of LT. The resulting BVP is solved by the polynomial Trefftz method (PTM), which can be regarded as a meshless method. In PTM, the source term is approximated by a truncated series of Chebyshev polynomials and the particular solution is obtained from a recursive procedure. Talbot’s method is employed for the numerical inversion of LT. The method is tested with the help of some numerical examples.     

New exact travelling wave solutions of two nonlinear physical models

In this paper, an improved tanh function method is used with a computerized symbolic computation for constructing new exact travelling wave solutions on two nonlinear physical models namely, the quantum Zakharov equations and the (2+1)-dimensional Broer–Kaup–Kupershmidt (BKK) system. The main idea of this method is to take full advantage of the Riccati equation which has more new solutions.The exact solutions are obtained which include new soliton-like solutions, trigonometric function solutions and rational solutions. The method is straightforward and concise, and its applications are promising.     

Comparison of generalized differential quadrature and Galerkin methods for the analysis of micro-electro-mechanical coupled systems

This paper presents a comprehensive comparison study between the generalized differential quadrature (GDQ) and the well-known global Galerkin method for analysis of pull-in behavior of nonlinear micro-electro-mechanical coupled systems. The nonlinear governing integro-differential equation for double clamped MEMS devices which was derived using variational principle by the authors [Sadeghian H, Rezazadeh G, Osterberg PM. Application of the generalized differential quadrature method to the study of pull-in phenomena of MEMS switches. J Microelectromech Syst 2007;16(6):1334–40] is discretized by applying Galerkin and GDQ methods. The divergence instability or pull-in phenomenon is analyzed. Obtained results are compared with the results of the pervious works. The Galerkin method is implemented with effect of number of used shape functions. Different types of trail functions on calculated pull-in voltage are examined.Furthermore, compare to one term and two terms truncation Galerkin method, it is observed that the GDQ with small number of grid points (non-uniform) performs accurate results for nonlinear micro-electro-mechanical coupled behavior which requires a large number of grid points at high-order approximation.     

弹性力学中的立兹法和屈列弗兹法的一般推导

本文采用一般的数学表示形式推导了线弹性力学中的立兹法和屈列弗兹法,证明了立兹法给出相应泛函极值的上限,屈列弗兹法则给出其下限.同时发现,特征值问题(例如自振频率问题)泛函变分法中的屈列弗兹法同求特征值的放松边界条件下限法是一致的.当然,此处的推导结果,也适用于一类泛函的变分法中,这类泛函的欧拉方程是线性正定的.     

Potential solutions of linear systems: The multi-criteria multiple constraint levels program

The realistic modeling of decision problems requires considerable flexibility in the model structure. Frequently one is faced with problems involving multiple criteria for which the constraint level is acceptable if a certain parameter (which may be a random variable) lies within a prescribed set. Furthermore, in formulating the problem, the criteria and constraints may be interchangeable. This requires a treatment which is more general than the nondominated solution in a multicriteria problem. We shall introduce the concept of a potential solution to cope with the above problem. To effectively locate these potential solutions, a generalization of the multicriteria (MC) simplex method, which handles multiple constraint levels (right hand sides) is developed. Geometric properties of adjacent potential solutions will be described together with a computational procedure which is based on the “connectedness” of the set of potential solutions. The natural duality relationship which exists in the double-MC simplex method and its consequences are also explored.     

基于核主元提取的支持向量机辨识

将一种基于特征提取的ε-不灵敏支持向量机方法用于非线性系统辨识.对输入输出数据首先进行核主元特征提取,将特征提取后的数据作为支持向量机的训练数据.将该方法与基于主元特征提取的方法和直接应用ε-不灵敏支持向量机的方法进行含噪和不含噪情况下的仿真比较,结果表明,方法的拟合性能和抗干扰能力优于其他两种方法.