一类有限元结构分析的EBE计算方法

本文利用ＥＢＥ策略￣［１］的基本思想，给出一类有限元结构分析的ＥＢＥ计算方法，即ＥＢＥ共轭梯度法ＥＢＥ－ＣＧ和ＥＢＥ预处理共轭梯度法ＥＢＥ－ＰＣＧ，这类方法避免了传统有限元结构分析中总刚度阵的组集而可大大降低存储量要求。同时它们还特别适合在各种粒度下的多处理机系统上实现。初步数值试验结果表明：这类ＥＢＥ有限元结构分析方法对串行和并行计算都是很有效的。     

Free vibration and dynamic stability of rotating thin-walled composite beams

The dynamic stability behavior of thin-walled rotating composite beams is studied by means of the finite element method. The analysis is based on Bolotin’s work on parametric instability for an axial periodic load. The influence of fiber orientation and rotating speeds on the natural frequencies and the unstable regions is studied for symmetrically balanced laminates. The regions of instability are obtained and expressed in non-dimensional terms. The “modal interchange” phenomenon arising in rotating beams is described. The dynamic stability problem is formulated by means of linearizing a geometrically nonlinear total Lagrangian finite element with seven degrees of freedom per node. This finite element formulation is based on a thin-walled beam theory that takes into account several non-classical effects such as anisotropy, shear flexibility and warping inhibition.     

面向对象的钢筋混凝土有限元非线性分析程序设计

采用面向对象的程序设计方法，结合钢筋混凝土有限元非线性分析模型，利用MFC(Microsoft Foundation Class Library)建立了有关描述钢筋混凝土有限元非线性分析的类，包括混凝土单元类、钢筋单元类、粘结单元类、非线性计算类，并给出了这些类的描述和它的实现方法，并为钢筋混凝土有限元非线性分析程序采用更合理的面向对象的方法提供了思路。     

广义特征值问题的并行EBE向量迭代法

本文利用ＥＢＥ策略和共轭梯度法，将广义特征值问题向量迭代法中各步的计算在单元级上进行，避免了总刚度和总质量矩阵的组集，大大节省了存储量。由此建立的ＥＢＥ向量迭代法尤其适宜于并行处理。数值算例结果表明无论是串行、还是并行计算，这类ＥＢＥ向量迭代法都能有效提高计算速度。     

High‐order 𝒞1 finite‐element interpolating schemes—Part II: Nonlinear semi‐Lagrangian shallow‐water models

The finite‐element, semi‐implicit, and semi‐Lagrangian methods are used on unstructured meshes to solve the nonlinear shallow‐water system. Several ??¹ approximation schemes are developed for an accurate treatment of the advection terms. The employed finite‐element discretization schemes are the P–P₁ and P₂–P₁ pairs. Triangular finite elements are attractive because of their flexibility for representing irregular boundaries and for local mesh refinement. By tracking the characteristics backward from both the interpolation and quadrature nodes and using ??¹ interpolating schemes, an accurate treatment of the nonlinear terms and, hence, of Rossby waves is obtained. Results of test problems to simulate slowly propagating Rossby modes illustrate the promise of the proposed approach in ocean modelling. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical analysis of elastic–plastic rotating disks with arbitrary variable thickness and density

A unified numerical method is developed in this article for the analysis of deformations and stresses in elastic–plastic rotating disks with arbitrary cross-sections of continuously variable thickness and arbitrarily variable density made of nonlinear strain-hardening materials. The method is based on a polynomial stress–plastic strain relation, deformation theory in plasticity and Von Mises’ yield condition. The governing equation is derived from the basic equations of the rotating disks and solved using the Runge–Kutta algorithm. The proposed method is applied to calculate the deformations and stresses in various rotating disks. These disks include solid disks with constant thickness and constant density, annular disks with constant thickness and constant density, nonlinearly variable thickness and nonlinearly variable density, linearly tapered thickness and linearly variable density, and a combined section of continuously variable thickness and constant density. The computed results are compared to those obtained from the finite element method and the existing approaches. A very good agreement is found between this research and the finite element analysis. Due to the simplicity, effectiveness and efficiency of the proposed method, it is especially suitable for the analysis of various rotating disks.     

A modified Galerkin/finite element method for the numerical solution of the Serre‐Green‐Naghdi system

A new modified Galerkin/finite element method is proposed for the numerical solution of the fully nonlinear shallow water wave equations. The new numerical method allows the use of low‐order Lagrange finite element spaces, despite the fact that the system contains third order spatial partial derivatives for the depth averaged velocity of the fluid. After studying the efficacy and the conservation properties of the new numerical method, we proceed with the validation of the new numerical model and boundary conditions by comparing the numerical solutions with laboratory experiments and with available theoretical asymptotic results. Copyright © 2016 John Wiley & Sons, Ltd.     

A design of residual error estimates for a high order BDF‐DGFE method applied to compressible flows

We deal with the numerical solution of the non‐stationary compressible Navier–Stokes equations with the aid of the backward difference formula – discontinuous Galerkin finite element method. This scheme is sufficiently stable, efficient and accurate with respect to the space as well as time coordinates. The nonlinear algebraic systems arising from the backward difference formula – discontinuous Galerkin finite element discretization are solved by an iterative Newton‐like method. The main benefit of this paper are residual error estimates that are able to identify the computational errors following from the space and time discretizations and from the inexact solution of the nonlinear algebraic systems. Thus, we propose an efficient algorithm where the algebraic, spatial and temporal errors are balanced. The computational performance of the proposed method is demonstrated by a list of numerical experiments. Copyright © 2013 John Wiley & Sons, Ltd.     

A constitutive model for plastically anisotropic solids with non-spherical voids

Plastic constitutive relations are derived for a class of anisotropic porous materials consisting of coaxial spheroidal voids, arbitrarily oriented relative to the embedding orthotropic matrix. The derivations are based on nonlinear homogenization, limit analysis and micromechanics. A variational principle is formulated for the yield criterion of the effective medium and specialized to a spheroidal representative volume element containing a confocal spheroidal void and subjected to uniform boundary deformation. To obtain closed form equations for the effective yield locus, approximations are introduced in the limit-analysis based on a restricted set of admissible microscopic velocity fields. Evolution laws are also derived for the microstructure, defined in terms of void volume fraction, aspect ratio and orientation, using material incompressibility and Eshelby-like concentration tensors. The new yield criterion is an extension of the well known isotropic Gurson model. It also extends previous analyses of uncoupled effects of void shape and material anisotropy on the effective plastic behavior of solids containing voids. Preliminary comparisons with finite element calculations of voided cells show that the model captures non-trivial effects of anisotropy heretofore not picked up by void growth models.     

Implicit symmetrized streamfunction formulations of magnetohydrodynamics

We apply the finite element method to the classic tilt instability problem of two‐dimensional, incompressible magnetohydrodynamics, using a streamfunction approach to enforce the divergence‐free conditions on the magnetic and velocity fields. We compare two formulations of the governing equations, the standard one based on streamfunctions and a hybrid formulation with velocities and magnetic field components. We use a finite element discretization on unstructured meshes and an implicit time discretization scheme. We use the PETSc library with index sets for parallelization. To solve the nonlinear problems on each time step, we compare two nonlinear Gauss‐Seidel‐type methods and Newton's method with several time‐step sizes. We use GMRES in PETSc with multigrid preconditioning to solve the linear subproblems within the nonlinear solvers. We also study the scalability of this simulation on a cluster. Copyright © 2008 John Wiley & Sons, Ltd.     

On uniformly accurate high‐order Boussinesq difference equations for water waves

A new accurate finite‐difference (AFD) numerical method is developed specifically for solving high‐order Boussinesq (HOB) equations. The method solves the water‐wave flow with much higher accuracy compared to the standard finite‐difference (SFD) method for the same computer resources. It is first developed for linear water waves and then for the nonlinear problem. It is presented for a horizontal bottom, but can be used for variable depth as well. The method can be developed for other equations as long as they use Padé approximation, for example extensions of the parabolic equation for acoustic wave problems. Finally, the results of the new method and the SFD method are compared with the accurate solution for nonlinear progressive waves over a horizontal bottom that is found using the stream function theory. The agreement of the AFD to the accurate solution is found to be excellent compared to the SFD solution. Copyright © 2005 John Wiley & Sons, Ltd.     

Perturbation solution of rotating solid disks with nonlinear strain-hardening

In order to determine the stresses and displacement in rotating disks with nonlinear strain-hardening, a polynomial stress-plastic strain relation is proposed and dimensionless quantities are adopted. A dimensionless governing equation for rotating disks is derived from the deformation theory of plasticity, Von Mises' yield criterion, proposed stress-plastic strain relation, equilibrium and compatibility equations. Its perturbation solution is obtained and applied to compute a rotating solid disk with nonlinear strain-hardening. The results are compared with those from the finite element calculation and other existing approaches, which confirms the validity of the method presented in this paper.     

SPH-FEM接触算法在冲击动力学数值计算中的应用

为了充分发挥光滑粒子流体动力学方法(Smoothed Particle Hydrodynamics,SPH)易于处理大变形以及有限元(Finite Element Method,FEM)计算精度和效率高的优势,论文基于无网格粒子接触算法,在有限元节点处设置背景粒子,通过接触力的方式计算SPH粒子和有限单元之间的接触问题...     

A stabilized finite element method for the Stokes problem based on polynomial pressure projections

A new stabilized finite element method for the Stokes problem is presented. The method is obtained by modification of the mixed variational equation by using local L² polynomial pressure projections. Our stabilization approach is motivated by the inherent inconsistency of equal‐order approximations for the Stokes equations, which leads to an unstable mixed finite element method. Application of pressure projections in conjunction with minimization of the pressure–velocity mismatch eliminates this inconsistency and leads to a stable variational formulation. Unlike other stabilization methods, the present approach does not require specification of a stabilization parameter or calculation of higher‐order derivatives, and always leads to a symmetric linear system. The new method can be implemented at the element level and for affine families of finite elements on simplicial grids it reduces to a simple modification of the weak continuity equation. Numerical results are presented for a variety of equal‐order continuous velocity and pressure elements in two and three dimensions. Copyright © 2004 John Wiley & Sons, Ltd.     

Nonlinear coupling modeling and dynamics analysis of rotating flexible beams with stretching deformation effect

Dynamic coupling modeling and analysis of rotating beams based on the nonlinear Green-Lagrangian strain are introduced in this work. With the reservation of the axial nonlinear strain, there are more coupling terms for axial and transverse deformations. The discretized dynamic governing equations are obtained by using the finite element method and Lagrange’s equations of the second kind. Time responses are conducted to compare the proposed model with other previous models. The stretching deforma...     

A Fourier–Chebyshev spectral collocation method for simulating flow past spheres and spheroids

An accurate Fourier–Chebyshev spectral collocation method has been developed for simulating flow past prolate spheroids. The incompressible Navier–Stokes equations are transformed to the prolate spheroidal co‐ordinate system and discretized on an orthogonal body fitted mesh. The infinite flow domain is truncated to a finite extent and a Chebyshev discretization is used in the wall‐normal direction. The azimuthal direction is periodic and a conventional Fourier expansion is used in this direction. The other wall‐tangential direction requires special treatment and a restricted Fourier expansion that satisfies the parity conditions across the poles is used. Issues including spatial and temporal discretization, efficient inversion of the pressure Poisson equation, outflow boundary condition and stability restriction at the pole are discussed. The solver has been validated primarily by simulating steady and unsteady flow past a sphere at various Reynolds numbers and comparing key quantities with corresponding data from experiments and other numerical simulations. Copyright © 1999 John Wiley & Sons, Ltd.     

Geometric nonlinear formulation and discretization method for a rectangular plate undergoing large overall motions

In the present paper, the geometric nonlinear formulation is developed for dynamic stiffening of a rectangular plate undergoing large overall motions. The dynamic equations, which take into account the stiffening terms, are derived based on the virtual power principle. Finite element method is employed for discretization of the plate. The simulation results of a rotating rectangular plate obtained by using such geometric nonlinear formulation are compared with those obtained by the conventional linear method without consideration of the stiffening effects. The application limit of the conventional linear method is clarified according to the frequency error. Furthermore, the accuracy of the assumed mode method is investigated by comparison of the results obtained by using the present finite element method and those obtained by using the assumed mode method.     

High‐order compact finite difference schemes for the vorticity–divergence representation of the spherical shallow water equations

This work is devoted to the application of the super compact finite difference method (SCFDM) and the combined compact finite difference method (CCFDM) for spatial differencing of the spherical shallow water equations in terms of vorticity, divergence, and height. The fourth‐order compact, the sixth‐order and eighth‐order SCFDM, and the sixth‐order and eighth‐order CCFDM schemes are used for the spatial differencing. To advance the solution in time, a semi‐implicit Runge–Kutta method is used. In addition, to control the nonlinear instability, an eighth‐order compact spatial filter is employed. For the numerical solution of the elliptic equations in the problem, a direct hybrid method, which consists of a high‐order compact scheme for spatial differencing in the latitude coordinate and a fast Fourier transform in longitude coordinate, is utilized. The accuracy and convergence rate for all methods are verified against exact analytical solutions. Qualitative and quantitative assessments of the results for an unstable barotropic mid‐latitude zonal jet employed as an initial condition are addressed. It is revealed that the sixth‐order and eighth‐order CCFDMs and SCFDMs lead to a remarkable improvement of the solution over the fourth‐order compact method. Copyright © 2015 John Wiley & Sons, Ltd.     

A preconditioned semi‐staggered dilation‐free finite volume method for the incompressible Navier–Stokes equations on all‐hexahedral elements

A new semi‐staggered finite volume method is presented for the solution of the incompressible Navier–Stokes equations on all‐quadrilateral (2D)/hexahedral (3D) meshes. The velocity components are defined at element node points while the pressure term is defined at element centroids. The continuity equation is satisfied exactly within each elements. The checkerboard pressure oscillations are prevented using a special filtering matrix as a preconditioner for the saddle‐point problem resulting from second‐order discretization of the incompressible Navier–Stokes equations. The preconditioned saddle‐point problem is solved using block preconditioners with GMRES solver. In order to achieve higher performance FORTRAN source code is based on highly efficient PETSc and HYPRE libraries. As test cases the 2D/3D lid‐driven cavity flow problem and the 3D flow past array of circular cylinders are solved in order to verify the accuracy of the proposed method. Copyright © 2005 John Wiley & Sons, Ltd.     

Acoustic simulation using a novel approach for reducing dispersion error

It is well‐known that the traditional finite element method (FEM) fails to provide accurate results to the Helmholtz equation with the increase of wave number because of the ‘pollution error’ caused by numerical dispersion. In order to overcome this deficiency, a gradient‐weighted finite element method (GW‐FEM) that combines Shepard interpolation and linear shape functions is proposed in this work. Three‐node triangular and four‐node tetrahedral elements that can be generated automatically are first used to discretize the problem domain in 2D and 3D spaces, respectively. For each independent element, a compacted support domain is then formed based on the element itself and its adjacent elements sharing common edges (or faces). With the aid of Shepard interpolation, a weighted acoustic gradient field is then formulated, which will be further used to construct the discretized system equations through the generalized Galerkin weak form. Numerical examples demonstrate that the present algorithm can significantly reduces the dispersion error in computational acoustics. Copyright © 2016 John Wiley & Sons, Ltd.