High-order theory for arched structures and its application for the study of the electrostatically actuated MEMS devices

A high-order theory for arched rods and beams based on expansion of the two-dimensional (2D) equations of elasticity into Legendre’s polynomials series has been developed. The 2D equations of elasticity have been expanded into Legendre’s polynomials series in terms of a thickness coordinate. Thereby, all equations of elasticity including Hooke’s law have been transformed to corresponding equations for coefficients of Legendre’s polynomials expansion. Then system of differential equations in term of displacements and boundary conditions for the coefficients of Legendre’s polynomials expansion coefficients has been obtained. Cases of the first and second approximations have been considered in details. For obtained boundary-value problems, a finite element method has been used and numerical calculations have been done with COMSOL Multiphysics and MATLAB. Developed theory has been applied for study pull-in instability and stress–strain state of the electrostatically actuated micro-electro-mechanical Systems.     

A high-order theory for functionally graded axially symmetric cylindrical shells

A high-order theory for functionally graded axially symmetric cylindrical shell based on expansion of the axially symmetric equations of elasticity for functionally graded materials into Legendre polynomials series has been developed. The axially symmetric equations of elasticity have been expanded into Legendre polynomials series in terms of a thickness coordinate. In the same way, functions that describe functionally graded relations has been also expanded. Thereby, all equations of elasticity including Hook’s law have been transformed to corresponding equations for coefficients of Legendre polynomials expansion. Then system of differential equations in terms of displacements and boundary conditions for the coefficients of Legendre polynomials expansion coefficients has been obtained. Cases of the first and second approximations have been considered in more details. For obtained boundary-value problems’ solution, a finite element has been used and numerical calculations have been done with COMSOL MULTIPHYSICS and MATLAB.     

A Legendre wavelet spectral collocation technique resolving anomalies associated with velocity in some boundary layer flows of Walter-B liquid

A Legendre wavelet spectral collocation method is proposed here to solve three boundary layer flow problems of Walter-B fluid namely the stagnation point flow, Blasius flow and Sakiadis flow. In the proposed method, we first transform the boundary value problems into initial value problems using shooting method. We then split the semi infinite domain into subintervals and the governing initial value problems are transformed to system of algebraic equations in each subinterval. The solutions of these algebraic equations yield an approximate solution of the differential equation in each subinterval. The overshoot in the velocity profile associated with the stagnation point and Blasius flows and undershoot in the Sakiadis flow is controlled. Physically realistic solutions are presented for both weakly and strongly viscoelastic parameters. The residual error validates the correctness, convergence and accuracy of the obtained solutions.     

On the representation of the solutions of the equilibrium equations of a piezoelectric transversally-isotropic spherical shell

A method is proposed for constructing equilibrium equations for thickness-polarized transversally isotropic piezoceramic shells. The method is based on Fourier expanding the required functions in Legendre polynomials. The appropriate system of differential equations is formed for the expansion coefficients as functions of two independent variables. The equilibrium equations are given in particular for transversally isotropic spherical shells. A method is given for constructing the general solution in the first approximation. S. P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine, Kiev, Ukraine. Translated from Prikladnaya Mekhanika, Vol. 35, No. 7, pp. 59–68, July, 1999.     

Free vibration of non-uniform axially functionally graded beams using the asymptotic development method

The asymptotic development method is applied to analyze the free vibration of non-uniform axially functionally graded(AFG) beams, of which the governing equations are differential equations with variable coefficients. By decomposing the variable flexural stiffness and mass per unit length into reference invariant and variant parts, the perturbation theory is introduced to obtain an approximate analytical formula of the natural frequencies of the non-uniform AFG beams with different boundary conditions.Furthermore, assuming polynomial distributions of Young's modulus and the mass density, the numerical results of the AFG beams with various taper ratios are obtained and compared with the published literature results. The discussion results illustrate that the proposed method yields an effective estimate of the first three order natural frequencies for the AFG tapered beams. However, the errors increase with the increase in the mode orders especially for the cases with variable heights. In brief, the asymptotic development method is verified to be simple and efficient to analytically study the free vibration of non-uniform AFG beams, and it could be used to analyze any tapered beams with an arbitrary varying cross width.     

Reissner型板中应力强度因子的近似方程和近似解法

本文在文献[1],[2]的基础上对Reissner型板进行分析,发现与经典板理论类似的近似方程用于求解Reissner型板的断裂问题是有效的,并用能量法解得受弯边裂纹和中心裂纹板的应力强度因子。将结果与文献[3]比较表明,用本文的近似方法求解应力强度因子方法简便且精度较高。     

Two Cracks Emerging from a Single Point,under the Influence of a Longitudinal Shear Wave

The problem of determining the dynamic stress intensity coefficients for two cracks emerging from a single point is solved. The cracks are affected by a longitudinal shear wave. The original problem is reduced to solving a system of two singular integro-differential equations with fixed singularities. For an approximate solution of this system, a numerical method is proposed that takes into account the real asymptotics of the unknown functions and uses special quadrature formulas for singular integrals.     

Transient stochastic response of quasi-partially integrable Hamiltonian systems

The approximate transient response of multi-degree-of-freedom (MDOF) quasi-partially integrable Hamiltonian systems under Gaussian white noise excitation is investigated. First, the averaged Itô equations for first integrals and the associated Fokker–Planck–Kolmogorov (FPK) equation governing the transient probability density of first integrals of the system are derived by applying the stochastic averaging method for quasi-partially integrable Hamiltonian systems. Then, the approximate solution of the transient probability density of first integrals of the system is obtained from solving the FPK equation by applying the Galerkin method. The approximate transient solution is expressed as a series in terms of properly selected base functions with time-dependent coefficients. The transient probability densities of displacements and velocities can be derived from that of first integrals. One example is given to illustrate the application of the proposed procedure. It is shown that the results for the example obtained by using the proposed procedure agree well with those from Monte Carlo simulation of the original system.     

Transient stochastic response of quasi integerable Hamiltonian systems

The approximate transient response of quasi integrable Hamiltonian systems under Gaussian white noise excitations is investigated. First, the averaged Ito equations for independent motion integrals and the associated Fokker-Planck-Kolmogorov (FPK) equation governing the transient probability density of independent motion integrals of the system are derived by applying the stochastic averaging method for quasi integrable Hamiltonian systems. Then, approximate solution of the transient probability density of independent motion integrals is obtained by applying the Galerkin method to solve the FPK equation. The approximate transient solution is expressed as a series in terms of properly selected base functions with time-dependent coefficients. The transient probability densities of displacements and velocities can be derived from that of independent motion integrals. Three examples are given to illustrate the application of the proposed procedure. It is shown that the results for the three examples obtained by using the proposed procedure agree well with those from Monte Carlo simulation of the original systems.     

Stability of anisotropic shells of revolution of positive or negative Gaussian curvature

A mixed variational principle is derived by Hamilton’s method from the principle of minimum potential energy for thin anisotropic shells of revolution and is then used to derive a normal system of equations with complex coefficients. Discrete orthogonalization is used to solve this homogeneous system and the nonlinear system of equations that describes the precritical state of shells. A shell generated by revolving a circular arc around the axis parallel to its chord is analyzed for stability. The solution is compared with the approximate solution obtained assuming that the precritical state is membrane. It is established that the approximate problem formulation gives incorrect results for shells of negative Gaussian curvature     

Construction of approximate analytical solutions to strongly nonlinear damped oscillators

An analytical approximate method for strongly nonlinear damped oscillators is proposed. By introducing phase and amplitude of oscillation as well as a bookkeeping parameter, we rewrite the governing equation into a partial differential equation with solution being a periodic function of the phase. Based on combination of the Newton’s method with the harmonic balance method, the partial differential equation is transformed into a set of linear ordinary differential equations in terms of harmonic coefficients, which can further be converted into systems of linear algebraic equations by using the bookkeeping parameter expansion. Only a few iterations can provide very accurate approximate analytical solutions even if the nonlinearity and damping are significant. The method can be applied to general oscillators with odd nonlinearities as well as even ones even without linear restoring force. Three examples are presented to illustrate the usefulness and effectiveness of the proposed method.     

Constructing invariant tori for two weakly coupled van der Pol oscillators

An algorithm is developed for the construction of an invariant torus of a weakly coupled autonomous oscillator. The system is put into angular standard form. The determining equations are found by averaging and are solved for the approximate amplitudes of the torus. A perturbation series is then constructed about the approximate amplitudes with unknown coefficients as periodic functions of the angular variables. A sequence of solvable partial differential equations is developed for determining the coefficients. The algorithm is applied to a system of nonlinearly coupled van der Pol equations and the first order coefficients are generated in a straightforward manner. The approximation shows both good numerical accuracy and reproducibility of the periodicities of the van der Pol system. A comparitive analysis of integrating the van der Pol system with integrating the phase equations from the angular standard form on the approximate torus shows numerical errors of the order of the perturbation parameter =0.05 for integrations of up to 10,000 steps. Applying FFT to the numerical periodicities generated by integrating the van der Pol system near the tours reveals the same predominant frequencies found in the perturbation coefficients. Finally an expected rotation number is found by integrating the phase equations on the approximate torus.Contribution of the National Institute of Standards and Technology, a Federal agency.     

基于整体位移模式的平面裂纹结构数值模拟方法

对于平面裂纹问题,针对扩展有限元法和无网格伽辽金法的不足,从结构的整体位移模式出发,提出了一种新的数值模拟方法。在整个求解域内构造其试探函数,并引入裂纹修正项描述裂尖处的奇异性和裂纹面的强间断特性;同时,提出了一种新的强制边界条件施加方法,通过引入位移边界水平集函数,将位移边界条件包含在近似位移场的表达式中,有效地解决了位移边界条件问题,减小了刚度矩阵的阶数,非常方便地消除了刚度矩阵的奇异性,降低了线性方程组的求解难度。含裂纹矩形平板结构的数值算例验证了该方法的有效性。     

Multidomain pseudospectral methods for nonlinear convection-diffusion equations

Multidomain pseudospectral approximations to nonlinear convection-diffusion equations are considered. The schemes are formulated with the Legendre-Galerkin method, but the nonlinear term is collocated at the Legendre/Chebyshev-Gauss-Lobatto points inside each subinterval. Appropriate base functions are introduced so that the matrix of the system is sparse, and the method can be implemented efficiently and in parallel. The stability and the optimal rate of convergence of the methods are proved. Numerical results are given for both the single domain and the multidomain methods to make a comparison.     

三维体内含垂直于双材料界面混合型裂纹的超奇异积分解法

利用双材料位移基本解和Somigliana公式,将三维体内含垂直于双材料界面混合型裂纹问题归结为求解一组超奇异积分方程。使用主部分析法,通过对裂纹前沿应力奇性的分析,得到用裂纹面位移间断表示的应力强度因子的计算公式,进而利用超奇异积分方程未知解的理论分析结果和有限部积分理论,给出了超奇异积分方程的数值求解方法。最后,对典型算例的应力强度因子做了计算,并讨论了应力强度因子数值结果的收敛性及其随各参数变化的规律。     

Analysis of dynamic stress intensity factors of three-point bend specimen containing crack

A new formula is obtained to calculate dynamic stress intensity factors of the three-point bending specimen containing a single edge crack in this study. Firstly, the weight function for three-point bending specimen containing a single edge crack is derived from a general weight function form and two reference stress intensity factors, the coefficients of the weight function are given. Secondly, the history and distribution of dynamic stresses in uncracked three-point bending specimen are derived based on the vibration theory. Finally, the dynamic stress intensity factors equations for three-pointing specimen with a single edge crack subjected to impact loadings are obtained by the weight function method. The obtained formula is verified by the comparison with the numerical results of the finite element method (FEM). Good agreements have been achieved. The law of dynamic stress intensity factors of the three-point bending specimen under impact loadings varing with crack depths and loading rates is studied.     

Nonlinear equations of elastic deformation of plates

A method for constructing nonlinear equations of elastic deformation of plates with boundary conditions for stresses and displacements at the face surfaces in an arbitrary coordinate system is proposed. The initial three–dimensional problem of the nonlinear theory of elasticity is reduced to a one–parameter sequence of two–dimensional problems by approximating the unknown functions by truncated series in Legendre polynomials. The same unknowns are approximated by different truncated series. In each approximation, a linearized system of equations whose differential order does not depend on the boundary conditions at the face surfaces which can be formulated in terms of stresses or displacements is obtained.     

An analytically-numerical approach for the analysis of an interface crack with a contact zone in a piezoelectric bimaterial compound

An interface crack with an artificial contact zone at the right-hand side crack tip between two dissimilar finite-sized piezoelectric materials is considered under remote mixed-mode loading. To find the singular electromechanical field at the crack tip, an asymptotic solution is derived in connection with the conventional finite element method. For mechanical loads, the stress intensity factors at the singular points are obtained. As a particular case of this solution, the contact zone model (in Comninou’s sense) is derived. A simple transcendental equation and an asymptotic formula for the determination of the real contact zone length are derived. The dependencies of the contact zone lengths on external load coefficients are illustrated in graphical form. For a particular case of a short crack with respect to the dimensions of the bimaterial compound, the numerical results are compared to the exact analytical solutions, obtained for a piezoelectric bimaterial plane with an interface crack.     

MHD flow of power-law fluid on moving surface with power-law velocity and special injection/blowing

The problem of magnetohydrodynamic （MHD） flow on a moving surface with the power-law velocity and special injection/blowing is investigated. A scaling group transformation is used to reduce the governing equations to a system of ordinary differen- tial equations. The skin friction coefficients of the MHD boundary layer flow are derived, and the approximate solutions of the flow characteristics are obtained with the homotopy analysis method （HAM）. The approximate solutions are easily computed by use of a high order iterative procedure, and the effects of the power-law index, the magnetic parameter, and the special suction/blowing parameter on the dynamics are analyzed. The obtained results are compared with the numerical results published in the literature, verifying the reliability of the approximate solutions.