Semi‐analytical integration of the 8‐node plane element stiffness matrix using symbolic computation

The semi‐analytical integration of an 8‐node plane strain finite element stiffness matrix is presented in this work. The element is assumed to be super‐parametric, having straight sides. Before carrying out the integration, the integral expressions are classified into several groups, thus avoiding duplication of calculations. Symbolic manipulation and integration is used to obtain the basic formulae to evaluate the stiffness matrix. Then, the resulting expressions are postprocessed, optimized, and simplified in order to reduce the computation time. Maple symbolic‐manipulation software was used to generate the closed expressions and to develop the corresponding Fortran code. Comparisons between semi‐analytical integration and numerical integration were made. It was demonstrated that semi‐analytical integration required less CPU time than conventional numerical integration (using Gaussian‐Legendre quadrature) to obtain the stiffness matrix. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Acoustic reflection from a baffled elastic/piezoelectric plate

A finite difference/boundary integral procedure to determine the acoustic reflected pressure from a fluid-loaded bi-laminated plate is described. The bi-laminate is composed of a piezoelectric layer and an elastic layer in contact with the fluid, and is held by an acoustically hard baffle. In the numerical model, the fluid pressure at fluid/solid interface is replaced by a continuum of point sources weighted by the normal acceleration of the elastic plate, and the governing equation system is solved in the solid domain. With the normal acceleration found, the reflected pressure in the fluid is determined by an integral expression involving the Green's function. It is demonstrated that an appropriate applied voltage potential across the piezoelectric layer has the effect of cancelling either the reflected or scattered pressure of the plate at any chosen field points in the fluid. Project supported by the National Natural Science Foundation of China (No. 10172039).     

烷基硫酸盐表面张力的量子化学研究

用AMl量子化学计算法优化了十三烷基硫酸根阴离子和——CH3在不同取代位的十二烷基硫酸根阴离子的几何构型，得到最优构型时最高占据分子轨道能级EHOMO、最低空轨道能级ELUMO、电子能量Eele和偶极矩μ等数据．将这些电子结构数据分别与表面张力相拟合，得到很好的相关性。文中讨论了——CH3在不同位置取代对表面张力的影响。     

Numerical solution of the three‐dimensional parabolic equation with an integral condition

Developement of numerical methods for obtaining approximate solutions to the three dimensional diffusion equation with an integral condition will be carried out. The numerical techniques discussed are based on the fully explicit (1,7) finite difference technique and the fully implicit (7,1) finite difference method and the (7,7) Crank‐Nicolson type finite difference formula. The new developed methods are tested on a problem. Truncation error analysis and numerical examples are used to illustrate the accuracy of the new algorithms. The results of numerical testing show that the numerical methods based on the finite difference techniques discussed in the present article produce good results. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 193–202, 2002; DOI 10.1002/num.1040     

解奇异非光滑方程组的牛顿法和不精确牛顿法

本文主要解决奇异非光滑方程组的解法。应用一种新的次微分的外逆，我们提出了牛顿法和不精确牛顿法，它们的收敛性同时也得到了证明。这种方法能更容易在一引起实际应用中实现。这种方法可以看作是已存在的解非光滑方程组的方法的延伸。     

Chemical bonding in isolated molecules and crystals of zwitterionic pentacoordinate silicon chelates

The electronic structures of a number of zwitterionic pentacoordinate silicon chelates were investigated using the results of X-ray diffraction studies and quantum-chemical calculatoins by the MPW1PW91/6-311G(d) method. The topological analysis of the electron density distribution function and the study in the framework of the natural bond orbital partitioning scheme showed that the character of chemical bonding in the axial fragments of the molecules under consideration changes from dative to three-center, four-electron as the silicon atom assumes a trigonal-bipyramidal coordination.     

Some degenerate and quasilinear parabolic systems not in divergence form

This paper deals with positive solutions of degenerate and quasilinear parabolic systems not in divergence form: u_t=u^p(Δu+av), v_t=v^q(Δv+bu), with null Dirichlet boundary conditions and positive initial conditions, where p, q, a and b are all positive constants. The local existence and uniqueness of classical solution are proved. Moreover, it will be proved that all solutions exist globally if and only if ab?λ₁², where λ₁ is the first eigenvalue of −Δ in Ω with homogeneous Dirichlet boundary condition.     

Modified Three-Dimensional Formulations of Bending Problems of Homogeneous Plates and Beams under Complex Fixing Conditions

Modified three-dimensional formulations of bending problems of homogeneous elastic plates and beams are considered. Modification of the known three-dimensional formulations reduces to using additional constraints imposed on displacement functions. An advantage of the formulations proposed is that complex fixing conditions of plates and beams can be taken into account.     

ROBUST GLOBAL EXPONENTIAL STABILITY OF UNCERTAIN IMPULSIVE SYSTEMS

By using the quasi-Lyapunov function, some sufficient conditions of global exponential stability for impulsive systems are established, which is the basis for the following discussion. Then, by employing Riccati inequality and Hamilton-Jacobi inequality approach, some sufficient conditions of robust exponential stability for uncertain linear/nonlinear impulsive systems are derived, respectively. Finally, some examples are given to illustrate the applications of the theory.     

Skew-Hermitian Triangular Splitting Iteration Methods for Non-Hermitian Positive Definite Linear Systems of Strong Skew-Hermitian Parts

By further generalizing the skew-symmetric triangular splitting iteration method studied by Krukier, Chikina and Belokon (Applied Numerical Mathematics, 41 (2002), pp. 89–105), in this paper, we present a new iteration scheme, called the modified skew-Hermitian triangular splitting iteration method, for solving the strongly non-Hermitian systems of linear equations with positive definite coefficient matrices. We discuss the convergence property and the optimal parameters of this new method in depth. Moreover, when it is applied to precondition the Krylov subspace methods like GMRES, the preconditioning property of the modified skew-Hermitian triangular splitting iteration is analyzed in detail. Numerical results show that, as both solver and preconditioner, the modified skew-Hermitian triangular splitting iteration method is very effective for solving large sparse positive definite systems of linear equations of strong skew-Hermitian parts.