固体力学有限元体系的结构拓扑变化理论

本文是文[1]的继续.文[1]提出了杆件系统的结构拓扑变化理论和拓扑变化法本文将这一理论和方法推进到连续体有限元体系;且在此基础上揭示出有限元体系的一个新性质,称为基本位移之梯度的正交性定理,从而给出一套设计敏度的显式表达式,可直接用于计算.     

Local and Parallel Finite Element Algorithms for Eigenvalue Problems

Abstract Some new local and parallel finite element algorithms are proposed and analyzed in this paper foreigenvalue problems.With these algorithms, the solution of an eigenvalue problem on a fine grid is reduced tothe solution of an eigenvalue problem on a relatively coarse grid together with solutions of some linear algebraicsystems on fine grid by using some local and parallel procedure.A theoretical tool for analyzing these algorithmsis some local error estimate that is also obtained in this paper for finite element approximations of eigenvectorson general shape-regular grids.     

The destabilizing effect of boundary slip on Bénard convection

We investigate the influence of slip boundary conditions on the onset of Bénard convection in an infinite fluid layer. It is shown that the critical Rayleigh number is a decreasing function of the slip length, and therefore boundary slip is seen to have a destabilizing effect. Chebyshev‐tau and compound matrix formulations for solving the eigenvalue problem are presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Fundamental processes in long air gap discharges

The development of atmospheric lightning is initiated and sustained by the formation in virgin air of ‘streamer corona’ and ‘leader’ discharges, very similar to those observed in laboratory long sparks. Therefore, the experimental and theoretical investigations of these laboratory discharges have become of large interest to improve the physical knowledge of the lightning process and to develop self-consistent models that could be applied to new protection concepts.In the present paper the fundamental processes of the subsequent phases of long air gap discharges are analyzed, from the first corona inception and development to the leader channel formation and propagation. For all these processes simulations models are discussed that have been essentially derived and simplified by the authors, in order to develop sequential time-dependent simulation of the laboratory breakdown, with both positive and negative voltages. The possibility of extending these models to the case of natural lightning is discussed in the companion paper, presented in this same volume. To cite this article: I. Gallimberti et al., C. R. Physique 3 (2002) 1335–1359.     

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).     

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.     

正交各向异性矩形板受迫振动的三维分析

从正交各向异性三维弹性动力学的控制方程出发，在求得四边简支矩形板自由振动频率和位移振型的基础上，构造了受迫振动的位移函数；利用自由振动位移振型的正交性，将控制方程的空间变量和时间变量分离，得到了广义质量、广义力和频率表示的关于时间的2阶常微分方程，从而得到了正交各向异性弹性矩形板在受迫振动下的位移场和应力场．给出了本文方法与薄板理论的算例．     

A note on the unsteady boundary layers over a flat plate

The unsteady boundary layer over a semi-infinite flat plate was investigated in this paper. The flow involves the unsteady flow over a flat plate with leading edge accretion or ablation. The momentum boundary layer was further analyzed and it was shown that the leading edge ablation had a similar effect to the wall mass injection or upstream wall movement making the fluid blown away from the wall. The thermal boundary layer of the same flow was also studied. Results show that the leading edge accretion or ablation can greatly change the fluid motion and the heat transfer characteristics.     

Numerical Prediction of the Strain State of Thermoset Matrix/Mineral Filler Composite Plates

The numerical prediction of the fields of inelastic strains (the linear invariant of the tensor of inelastic strains) in thermoset polyester/marble filler composite plates is discussed. A uniformly distributed load is applied to the plates, which lie on a steel base. The strain fields are predicted by means of the boundary element method by using creep test data for the composites and the polyester matrix itself. Identical creep tests were performed for two ages of the materials (1 month and 13 years), which allowed evaluating the aging effect. The study is carried out in two stages. At the first stage, the application of the generalized Maxwell-Gurevich equation to the thermoset matrix/mineral filler composite is demonstrated. The model parameters determined from the experimental creep data is used for the second stage, where the state of inelastic strains in the plates is predicted by applying the boundary element method. The influence of composite formulation (filler content) and physical aging of the polyester matrix on the state of inelastic strains in the plates is shown.__________Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 41, No. 2, pp. 145–156, March–April, 2005.