拟弧长延拓法在静电激励MEMS吸合特性研究中的应用

在静电激励微机电系统MEMS(micro-electro-mechanical systems)吸合特性研究中,基于应变梯度理论的微梁结构的控制方程是非线性高阶微分方程,给方程的求解带来了困难.由于该问题的数学模型本质上是分叉问题,方程的解支上出现奇异点,而运用局部延拓法无法通过奇异点.因此,通过运用广义微分求积法将控制方程降阶离散,结合拟弧长延拓法使迭代顺利通过奇异点,求出了整个解曲线.结果表明,拟弧长延拓法能有效并准确地求解具有分叉现象的高阶微分方程问题,为精确预测静电激励MEMS的吸合电压提供有力帮助.     

A novel numerical solution strategy for solving nonlinear free and forced vibration problems of cylindrical shells

Nonlinear vibration analysis of circular cylindrical shells has received considerable attention from researchers for many decades. Analytical approaches developed to solve such problem, even not involved simplifying assumptions, are still far from sufficiency, and an efficient numerical scheme capable of solving the problem is worthy of development. The present article aims at devising a novel numerical solution strategy to describe the nonlinear free and forced vibrations of cylindrical shells. For this purpose, the energy functional of the structure is derived based on the first-order shear deformation theory and the von–Kármán geometric nonlinearity. The governing equations are discretized employing the generalized differential quadrature (GDQ) method and periodic differential operators along axial and circumferential directions, respectively. Then, based on Hamilton's principle and by the use of variational differential quadrature (VDQ) method, the discretized nonlinear governing equations are obtained. Finally, a time periodic discretization is performed and the frequency response of the cylindrical shell with different boundary conditions is determined by applying the pseudo-arc length continuation method. After revealing the efficiency and accuracy of the proposed numerical approach, comprehensive results are presented to study the influences of the model parameters such as thickness-to-radius, length-to-radius ratios and boundary conditions on the nonlinear vibration behavior of the cylindrical shells. The results indicate that variation of fundamental vibrational mode shape significantly affects frequency response curves of cylindrical shells.     

Nonlinear modeling and analysis of electrically actuated viscoelastic microbeams based on the modified couple stress theory

A nonclassical nonlinear continuum model of electrically actuated viscoelastic microbeams is presented based on the modified couple stress theory to consider the microstructure effect in the framework of viscoelasticity. The nonlinear integral-differential governing equation and related boundary conditions of are derived based on the extended Hamilton's principle and Euler–Bernoulli hypothesis for viscoelastic microbeams with clamped-free, clamped-clamped, simply-supported boundary conditions. The proposed model accounts for system nonlinearities including the axial residual stress, geometric nonlinearity due to midplane stretching, electrical forcing with fringing effect. The behavior of the microbeam is simulated using generalized Maxwell viscoelastic model. A new generalized differential/integral quadrature method is developed to solve the resulting governing equation. The developed model is verified against elastic behavior and a favorable agreement is obtained. Efficiency of the developed model is demonstrated by analyzing the quasistatic pull-in phenomena of electrically actuated viscoelastic microbeams with different boundaries at various material length scale parameters and axial residual stresses in the framework of linear viscoelasticity.     

Pull-in instability and free vibration of electrically actuated poly-SiGe graded micro-beams with a curved ground electrode

This paper investigates the pull-in instability and free vibration of functionally graded poly-SiGe micro-beams under combined electrostatic force, intermolecular force and axial residual stress, with an emphasis on the effects of ground electrode shape, position-dependent material composition, and geometrically nonlinear deformation of the micro-beam. The differential quadrature (DQ) method is employed to solve the nonlinear differential governing equations to obtain the pull-in voltage and vibration frequencies of the clamped poly-SiGe micro-beams. The present analysis is validated through direct comparisons with published experimental results and excellent agreement has been achieved. A parametric study is conducted to investigate the effects of material composition, ground electrode shape, axial residual stress and geometrical nonlinearity on the pull-in voltage and frequency characteristics.     

Semi-analytic actuating and sensing in regular and irregular MEMs,single and assembled micro cantilevers

The static and dynamic behavior of regular and irregular single and assembled micro cantilever probes (ACPs) have been analysed. Various points and distributed loadings are considered. Since the applications of Micro Electro Mechanical Systems (MEMS) are not limited to an especial boundary condition, two semi-analytical approaches, named the generalized differential quadrature and generalized differential quadrature element methods (GDQM and GDQEM) have been used for regular and irregular MEMS, respectively. With less computational cost, it has been clearly demonstrated that these methods are more accurate than common numerical methods such as the Finite Element Method (FEM). For probable cases, proposed approaches have been validated with the exact Green’s function method. Then, considering the various composite lamination configurations (angle/cross), the effects of the electromechanical loading on the nano steering devices have been introduced. At the end, the solution challenges for scanning (sensing) the especial micro and nano profiles has been discussed and as a general case, a nano gear has been studied. The phase plane approach shows the probability of solution for various configurations and suggests the best for more stability.     

Variational differential quadrature: A technique to simplify numerical analysis of structures

A new solution method in the area of computational mechanics is developed in this article, which is called variational differential quadrature (VDQ). The main idea of this method is based on the accurate and direct discretization of the energy functional in the structural mechanics. In the VDQ method, through developing an efficient matrix formulation and using an accurate integral operator, the discretized governing equations are derived directly from the weak form of the equations with no need for the analytical derivation of the strong form. This technique provides an alternative way to discretize the energy functional, which avoids the local interpolation and the assembly process of the methods of this kind. We first implement the VDQ method for the nonlinear elasticity theory considering the Green-St. Venant strain tensor; then we simplify the formulation further for the first-order shear deformable beam and plate theories. The final formulation of these cases demonstrates the simplicity of the implementation for the VDQ method in the numerical analysis of the structures, which is a major goal for this article. Using these examples, one can easily learn and apply this technique to other structures. To assess the performance of the VDQ method, we compare it with the generalized differential quadrature (GDQ) method and finite element method (FEM) in the case of bending analysis of Mindlin plates. It is indicated that computational cost of VDQ is less than that of GDQ, and the convergence rate of VDQ is faster than that of FEM.     

The differential quadrature element method irregular element torsion analysis model

The differential quadrature element method (DQEM) has been proposed. The element weighting coefficient matrices are generated by the differential quadrature (DQ) or generic differential quadrature (GDQ). By using the DQ or GDQ technique and the mapping procedure the governing differential or partial differential equations, the transition conditions of two adjacent elements and the boundary conditions can be discretized. A global algebraic equation system can be obtained by assembling all of the discretized equations. This method can convert a generic engineering or scientific problem having an arbitrary domain configuration into a computer algorithm. The DQEM irregular element torsion analysis model is developed.     

The comparison between the DRBEM and DQM solution of nonlinear reaction–diffusion equation

In this study, both the dual reciprocity boundary element method and the differential quadrature method are used to discretize spatially, initial and boundary value problems defined by single and system of nonlinear reaction–diffusion equations. The aim is to compare boundary only and a domain discretization method in terms of accuracy of solutions and computational cost. As the time integration scheme, the finite element method is used achieving solution in terms of time block with considerably large time steps. The comparison between the dual reciprocity boundary element method and the differential quadrature method solutions are made on some test problems. The results show that both methods achieve almost the same accuracy when they are combined with finite element method time discretization. However, as a method providing very good accuracy with considerably small number of grid points differential quadrature method is preferrable.     

Size-dependent buckling analysis of functionally graded third-order shear deformable microbeams including thermal environment effect

Presented herein is the prediction of buckling behavior of size-dependent microbeams made of functionally graded materials (FGMs) including thermal environment effect. To this purpose, strain gradient elasticity theory is incorporated into the classical third-order shear deformation beam theory to develop a non-classical beam model which contains three additional internal material length scale parameters to consider the effects of size dependencies. The higher-order governing differential equations are derived on the basis of Hamilton’s principle. Afterward, the size-dependent differential equations and related boundary conditions are discretized along with commonly used end supports by employing generalized differential quadrature (GDQ) method. A parametric study is carried out to demonstrate the influences of the dimensionless length scale parameter, material property gradient index, temperature change, length-to-thickness aspect ratio and end supports on the buckling characteristics of FGM microbeams. It is revealed that temperature change plays more important role in the buckling behavior of FGM microbeams with higher values of dimensionless length scale parameter.     

Free vibration analysis of two-dimensional functionally graded cylindrical shells

In this paper, the free vibration of a two-dimensional functionally graded circular cylindrical shell is analyzed. The equations of motion are based on the Love’s first approximation classical shell theory. The spatial derivatives of the equations of motion and boundary conditions are discretized by the methods of generalized differential quadrature (GDQ) and generalized integral quadrature (GIQ). Two kinds of micromechanics models, viz. Voigt and Mori–Tanaka models are used to describe the material properties. To validate the results, comparisons are made with the solutions for FG cylindrical shells available in the literature. The results of this study show that the natural frequency of the material can be modified in order to meet the expected results through manipulation of the constituent volume fractions. A comprehensive comparison is then drawn between ordinary and 2-D FG cylindrical shells.     

Free vibration analysis of cracked composite beams subjected to coupled bending–torsion loads based on a first order shear deformation theory

In this paper, free vibration analysis of cracked composite beam subjected to coupled bending–torsion loading is presented. The composite beam is assumed to have an open edge crack of length a. A first order shear deformation theory is applied to count for the effect of shear deformations on natural frequencies as well as the effect of coupling in torsion and bending modes of vibration. Governing equations and boundary conditions are derived using Hamilton principle. Local flexibility matrix is used to obtain the additional boundary conditions of the beam in cracked area. After obtaining the governing equations and boundary conditions, generalized differential quadrature (GDQ) method is applied to solve the obtained eigenvalue problem. Finally, some numerical results of beams with various boundary conditions and different fiber orientations are given to show the efficiency of the method. In addition, to study the effect of shear deformations, numerical results of the current model are compared with previously given results in which shear deformations were neglected.     

Comparison of the finite difference and Galerkin methods as applied to the solution of the hydrodynamic equations

The three-dimensional nonlinear hydrodynamic equations which describe wind induced flow in a homogeneous sea are transformed from Cartesian coordinates into sigma coordinates. The solution of these equations in the horizontal is accomplished using a standard finite difference grid and established finite difference methods.The accuracy and computational efficiency, in terms of both computer time and main memory requirements, of using either the Galerkin method or a finite difference grid through the vertical is considered. Calculations, using the same number of functions in the Galerkin method as grid bases through the vertical shows that the Galerkin method has superior accuracy over the grid box method. Hence, for a given accuracy a smaller number of functions than grid boxes may be used, with associated saving in computational resources.For the case in which the vertical variation of eddy viscosity is fixed, an eigenvalue problem can be solved to yield a set of eigenfunctions. Using these eigenfunctions as a basis set with the Galerkin approach, a Galerkin-eigenfunction method is developed. Calculations show that the Galerkin-eigenfunction technique is accurate and in a linear model is clearly computationally more economic than the use of grid boxes through the vertical.     

An efficient linear solver for the hybridized discontinuous Galerkin method

Discretizing partial differential equations by an implicit solving technique ultimately leads to a linear system of equations that has to be solved. The number of globally coupled unknowns is especially large for discontinuous Galerkin (DG) methods. It can be reduced by using hybridized discontinuous Galerkin (HDG) methods, but still efficient linear solvers are needed. It has been shown that, if hierarchical basis functions are used, a hierarchical scale separation (HSS) ansatz can be an efficient solver. In this work, we couple the HDG method with an HSS solver to solve a scalar nonlinear problem. It is validated by comparing the results with results obtained by GMRES with ILU(3) preconditioning as linear solver. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

求解粘性流体和热迁移联立方程的迎风局部微分求积法

微分求积方法(DQM)已成功地应用于数值求解流体力学中的许多问题．但是已有的工作大多限于正规区域的流动问题,同时缺少用迎风机制来描述流体流动的对流特性．该文对一个不规则区域中的不可压缩层流和热迁移的耦合问题给出了一种具有迎风机制的局部微分求积方法,对通过边界和坐标不平行的收缩管道中的流体,只用少数网格点得到了比较好的数值解．和有限差分方法（FDM）相比较,这一方法具有计算工作量少、存储量小和收敛性好等优点．     

Novel weak form quadrature element method with expanded Chebyshev nodes

Based on the principle of minimum potential energy and the differential quadrature rule, novel weak form quadrature element method is proposed. Different from the existing ones, expanded Chebyshev grid points are used as the element nodes. A simple but general way is proposed to compute the strains at the integration points explicitly by using the differential quadrature rule. For illustration and verification, quadrature bar and beam elements are established. Several examples are given. Numerical results indicate that the proposed quadrature element method allows a longer time step as compared to elements with other nodes and is an accurate and efficient method for structural analysis.     

A Legendre–Galerkin Chebyshev collocation method for the Burgers equation with a random perturbation on boundary condition

In the paper, we apply the generalized polynomial chaos expansion and spectral methods to the Burgers equation with a random perturbation on its left boundary condition. Firstly, the stochastic Galerkin method combined with the Legendre–Galerkin Chebyshev collocation scheme is adopted, which means that the original equation is transformed to the deterministic nonlinear equations by the stochastic Galerkin method and the Legendre–Galerkin Chebyshev collocation scheme is used to deal with the resulting nonlinear equations. Secondly, the stochastic Legendre–Galerkin Chebyshev collocation scheme is developed for solving the stochastic Burgers equation; that is, the stochastic Legendre–Galerkin method is used to discrete the random variable meanwhile the nonlinear term is interpolated through the Chebyshev–Gauss points. Then a set of deterministic linear equations can be obtained, which is in contrast to the other existing methods for the stochastic Burgers equation. The mean square convergence of the former method is analyzed. Numerical experiments are performed to show the effectiveness of our two methods. Both methods provide alternative approaches to deal with the stochastic differential equations with nonlinear terms.     

Buckling and free vibration analysis of high speed rotating carbon nanotube reinforced cylindrical piezoelectric shell

In this paper, buckling and free vibration behavior of a piezoelectric rotating cylindrical carbon nanotube-reinforced (CNTRC) shell is investigated. Both cases of uniform distribution (UD) and FG distribution patterns of reinforcements are studied. The accuracy of the presented model is verified with previous studies and also with those obtained by Navier analytical method. The novelty of this study is investigating the effects of critical voltage and CNT reinforcement as well as satisfying various boundary conditions implemented on the piezoelectric rotating cylindrical CNTRC shell. The governing equations and boundary conditions have been developed using Hamilton's principle and are solved with the aid of Navier and generalized differential quadrature (GDQ) methods. In this research, the buckling phenomena in the piezoelectric rotating cylindrical CNTRC shell occur as the natural frequency is equal to zero. The results show that, various types of CNT reinforcement, length to radius ratio, external voltage, angular velocity, initial hoop tension and boundary conditions play important roles on critical voltage and natural frequency of piezoelectric rotating cylindrical CNTRC shell.     

Design of a nonsingular adaptive fuzzy backstepping controller for electrostatically actuated microplates

This paper adopts some alternative strategies to design a nonlinear controller for double electrostatically actuated microplates. The novel design is carried out to solve the singularity problem reported in many articles due to the use of the Taylor expansion to simplify the electrostatic force. The nonlinear governing partial differential equation is converted to the modal equation using the Galerkin method. Then, based on the Lyapunov stability criterion, a fuzzy backstepping controller facilitated by prescribed performance functions is applied to the non-affine system to extend the travel range beyond the pull-in region and capture the structural and nonstructural uncertainties that exist in the practical systems. The present work also aims to bring satisfactory transient and steady-state performance indices to the system. Moreover, unknown time-varying delays as the indispensable part of practical systems are considered in the proposed control scheme to suppress the delays occurring in the measurement of the states by constructing Lyapunov–Krasovskii function. The accuracy of the modal equation in both the static and dynamic analysis is verified through a meshless method as a direct solution of the partial differential equation. The proposed controller guarantees that all the closed-loop signals are semi-globally, uniformly ultimately bounded, and the error evolves within the decaying prescribed bounds. Finally, the proposed controller demonstrates its feasibility to extend the travel range within and beyond the pull-in range despite the unknown uncertainties and time-varying delays which exist in the system.     

Hybrid Galerkin boundary elements: theory and implementation

Summary. In this paper we present a new quadrature method for computing Galerkin stiffness matrices arising from the discretisation of 3D boundary integral equations using continuous piecewise linear boundary elements. This rule takes as points some subset of the nodes of the mesh and can be used for computing non-singular Galerkin integrals corresponding to pairs of basis functions with non-intersecting supports. When this new rule is combined with standard methods for the singular Galerkin integrals we obtain a “hybrid” Galerkin method which has the same stability and asymptotic convergence properties as the true Galerkin method but a complexity more akin to that of a collocation or Nystr?m method. The method can be applied to a wide range of singular and weakly-singular first- and second-kind equations, including many for which the classical Nystr?m method is not even defined. The results apply to equations on piecewise-smooth Lipschitz boundaries, and to non-quasiuniform (but shape-regular) meshes. A by-product of the analysis is a stability theory for quadrature rules of precision 1 and 2 based on arbitrary points in the plane. Numerical experiments demonstrate that the new method realises the performance expected from the theory. Received January 22, 1998 / Revised version received May 26, 1999 / Published online April 20, 2000 –? Springer-Verlag 2000     

On shifted Jacobi spectral method for high-order multi-point boundary value problems

This paper reports a spectral tau method for numerically solving multi-point boundary value problems (BVPs) of linear high-order ordinary differential equations. The construction of the shifted Jacobi tau approximation is based on conventional differentiation. This use of differentiation allows the imposition of the governing equation at the whole set of grid points and the straight forward implementation of multiple boundary conditions. Extension of the tau method for high-order multi-point BVPs with variable coefficients is treated using the shifted Jacobi Gauss–Lobatto quadrature. Shifted Jacobi collocation method is developed for solving nonlinear high-order multi-point BVPs. The performance of the proposed methods is investigated by considering several examples. Accurate results and high convergence rates are achieved.