抛物型积分微分方程各向异性非协调有限元分析

各向异性有限元方法的显著的优点之一就是可以用较少的自由度得到与传统有限元正则剖分时同样的估计结果.然而,在这种情况下,Sobolev空间上的:Bramble-Hilbert引理在插值误差分析中不能直接应用,而且对于非协调元来说其传统边界估计技巧也不再适用.本文证明了一个非协调单元具有各向异性特征,并将它应用到研究抛物积分微分方程半离散格式下的Galerkin逼近.利用单元的特殊性,验证了Ritz-Volterra投影与有限元插值是相同的.在解适当光滑时,通过引入一些新的技巧,得到了与传统方法相同的收敛误差估计和超逼近性质.最后,通过构造适当的插值后处理算子,得到了各向异性网格下的整体超收敛结果.该文的结果对进一步探索和设计数值的自适应算法是有帮助的.     

单层扁锥面网壳非线性动力稳定性分析

用拟壳法建立了正三角形网格的三向扁锥面单层网壳的非线性动力学微分方程。在周边固定条件下,用分离变量函数法给出网壳的横向位移。由协调方程求出张力,通过Galerkin作用得到了一个含二次、三次的非线性微分方程,在不考虑外激励情况下,此系统有三个平衡点。通过求Floquet指数讨论了零平衡点邻域的稳定性问题。为了研究系统的混沌运动,在给定的初始条件下,对此动力系统的非线性自由振动方程进行了求解,首次得到了带平方和立方非线性系统的准确解,使得求Melnikov函数成为可能。用复变函数中的留数理论求出了Melnikov函数,得到了发生混沌的临界条件,通过数值仿真和Poincare映射也证实了混沌运动的存在。     

非线性双曲型积分微分方程的各向异性非协调有限元逼近

在各向异性网格下,讨论了一类非线性双曲型积分微分方程的一个矩形非协调有限元方法逼近,给出了半离散格式下的有限元解的收敛性分析和误差估计。在精确解适当光滑的前提下,利用新的技巧和精细估计得到了其超逼近性质。同时利用插值后处理技术导出了整体超收敛结果。本文的结论表明传统有限元分析中对网格的正则性要求和对Ritz-Volterra投影的依赖不是必要的,从而进一步扩展了非协调有限元方法的应用范围。     

非定常不可压缩Navier-Stokes方程的特征稳定化非协调有限元方法（英文）

数值求解非定常不可压缩Navier-Stokes方程的难点之一在于强烈的非线性容易引发非物理震荡,本文结合可以有效减弱此种震荡的特征线离散方法,基于局部Gauss积分之差的稳定化格式,采用最低等阶非协调混合有限元对NCP1-P1,构造出求解非定常不可压缩Navier-Stokes方程的特征稳定化非协调混合有限元方法.证明了该方法的全离散格式是无条件稳定的,并给出逼近解的相应误差估计.     

缓增分数阶扩散方程的高阶时间离散LDG方法

研究了缓增分数阶扩散方程的高阶时间离散局部间断Galerkin (Local Discontinuous Galerkin, LDG)方法，不是直接求解缓增分数阶扩散方程，而是首先通过变换将其转化成Caputo型时间分数阶扩散方程。接着，采用L1-2差分逼近离散Caputo型分数阶导数，间断有限元离散空间变量，构造求解模型的全离散LDG格式。证明了所建立的全离散格式为无条件稳定的且具有最优误差阶，两个数值算了验证了所建立数值格式的精度和鲁棒性。数值实验结果表明所建立格式在时间和空间方向均具有高精度。     

Cahn-Hilliard方程的时间双层网格有限元方法

针对用非线性数值格式求解Cahn-Hilliard方程时由非线性迭代引起的耗时问题,本文提出了一种时间双层网格(TT-M)有限元(FE)方法.该方法分为两步:第一步,在粗的时间步长上求解非线性Cahn-Hilliard系统,其中空间离散采用有限元方法,时间离散采用Crank-Nicolson格式;第二步,在细的时间步长上求解线性系统,然后证明了该方法的稳定性和误差估计,并通过数值算例对理论部分进行验证.结果表明,与传统的Galerkin有限元方法相比,该方法可以节省计算时间,说明了该方法的有效性和可行性.     

非定常不可压缩Navier-Stokes方程的特征稳定化非协调有限元方法

数值求解非定常不可压缩Navier-Stokes方程的难点之一在于强烈的非线性容易引发非物理震荡,本文结合可以有效减弱此种震荡的特征线离散方法,基于局部Gauss积分之差的稳定化格式,采用最低等阶非协调混合有限元对NCP1-P1,构造出求解非定常不可压缩Navier-Stokes方程的特征稳定化非协调混合有限元方法。证明了该方法的全离散格式是无条件稳定的,并给出逼近解的相应误差估计。     

非线性拟双曲方程的有限元配置法数值分析

基于有限元配置法,采用分片双三次Hermite插值多项式空间作为逼近函数空间,本文对粘性振动及神经传播过程中涉及的一类非线性拟双曲方程的初边值问题建立了二维半离散和全离散格式.并对两种格式证明了数值解的存在唯一性,应用微分方程先验估计的理论和技巧得到了L2模最佳阶误差估计.数值实验结果表明:所提方法在保证整体误差估计要求且不增加计算量的前提下,比传统有限元方法有更高的逼近精度,并扩展了配置法的应用范围.     

非定常N-S方程的双重网格数值模拟

本文采用全离散双重网格算法（时间变量采用Eular全隐式格式离散,空间变量采用混合有限元离散）,对非定常Navier-Stokes（N-S）方程进行数值模拟.双重网格算法的基本思想是,首先在粗网格有限元空间X^H上求解一个非线性问题,然后在细网格有限元空间Xh（h＜＜H）上求解一个线性问题.数值实验结果表明：在保持几乎相同精度的前提下,双重网格算法比标准有限元算法节省近一半的计算时间,说明了新算法求解非定常N-S方程的可行性和高效性.     

圆锥壳-圆柱壳-球壳组合结构自由振动分析

基于Reissner薄壳理论,采用区域分解法分析了不同边界条件下圆锥壳-圆柱壳-球壳组合结构的自由振动。首先在壳体连接处将组合壳体分为独立的圆锥壳、圆柱壳和球壳,并将各个子壳体沿旋转轴线分解为若干自由壳段;然后将所有壳段分区界面(包括边界界面)的位移协调方程通过分区广义变分和最小二乘加权残值法引入到组合壳体的能量泛函中;最后将壳段位移场变量的周向分量和轴向分量分别以Fourier级数和Chebyshev多项式展开,通过变分后得到整个组合壳体的离散动力学方程。将区域分解法计算结果与有限元软件ANSYS计算结果进行对比,验证了区域分解法在分析圆锥壳-圆柱壳-球壳组合结构自由振动的正确性和计算精度,并分析了组合壳体长径比及厚径比对自由振动频率的影响。     

A 3‐node flat triangular shell element with corner drilling freedoms and transverse shear correction

The formulation, implementation and testing of simple, efficient and robust shell finite elements have challenged investigators over the past four decades. A new 3‐node flat triangular shell element is developed by combination of a membrane component and a plate bending component. The ANDES‐based membrane component includes rotational degrees of freedom, and the refined nonconforming element method‐based bending component involves a transverse shear correction. Numerical examples are carried out for benchmark tests. The results show that compared with some popular shell elements, the present one is simple but exhibits excellent all‐around properties (for both membrane‐and bending‐dominated situations), such as free of aspect ratio locking, passing the patch test, free of shear locking, good convergence and high suitability for thin to moderately thick plates. The developed element has already been adopted in a warpage simulation package for injection molding. Copyright © 2011 John Wiley & Sons, Ltd.     

A triangular finite shell element based on a fully nonlinear shell formulation

This work presents a fully nonlinear six-parameter (3 displacements and 3 rotations) shell model for finite deformations together with a triangular shell finite element for the solution of the resulting static boundary value problem. Our approach defines energetically conjugated generalized cross-sectional stresses and strains, incorporating first-order shear deformations for an inextensible shell director (no thickness change). Finite rotations are treated by the Euler–Rodrigues formula in a very convenient way, and alternative parameterizations are also discussed herein. Condensation of the three-dimensional finite strain constitutive equations is performed by applying a mathematically consistent plane stress condition, which does not destroy the symmetry of the linearized weak form. The results are general and can be easily extended to inelastic shells once a stress integration scheme within a time step is at hand. A special displacement-based triangular shell element with 6 nodes is furthermore introduced. The element has a nonconforming linear rotation field and a compatible quadratic interpolation scheme for the displacements. Locking is not observed as the performance of the element is assessed by several numerical examples, which also illustrate the robustness of our formulation. We believe that the combination of reliable triangular shell elements with powerful mesh generators is an excellent tool for nonlinear finite element analysis.Fellowship funding from FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo) and CNPq (Conselho Nacional de Pesquisa), together with the material support and stimulating discussions in IBNM (Institut für Baumechanik und Numerische Mechanik), are gratefully acknowledged in this work.     

An efficient hybrid / mixed element for geometrically nonlinear analysis of plate and shell structures

In this paper, a geometrically nonlinear hybrid/mixed curved quadrilateral shell element (HMSHEL4N) with four nodes is developed based on the modified Hellinger/Reissner variational principles. The performance of element is investigated and tested using some benchmark problems. A number of numerical examples of plate and shell nonlinear deflection problems are included. The results are compared with theoretical solutions and other numerical results. It is shown that HMSHEL4N does not possess spurious zero energy modes and any locking phenomenon, and is convergent and insensitive to the distorted mesh. A good agreement of the results with theoretical solutions, and better performance compared with displacement finite element method, are observed. It is seen that an efficient shell element based on stress and displacement field assumptions in solution and time is obtained.     

Rotation‐free triangular shell element using node‐based smoothed finite element method

In this paper, a triangular thin flat shell element without rotation degrees of freedom is proposed. In the Kirchhoff hypothesis, the first derivative of the displacement must be continuous because there are second‐order differential terms of the displacement in the weak form of the governing equations. The displacement is expressed as a linear function and the nodal rotation is defined using node‐based smoothed finite element method. The rotation field is approximated using the nodal rotation and linear shape functions. This rotation field is linear in an element and continuous between elements. The curvature is defined by differentiating the rotation field, and the stiffness is calculated from the curvature. A hybrid stress triangular membrane element was used to construct the shell element. The penalty technique was used to apply the rotation boundary conditions. The proposed element was verified through several numerical examples.     

Analysis of laminated conical shell structures using higher order models

In this paper is presented a numerical method for the structural analysis of laminated conical shell panels using a quadrilateral isoparametric finite element based on the higher order shear deformation theory. The displacement expressions used for the longitudinal and circumferential components of the displacement field are given by power series of the transversal coordinate and the condition of zero stresses in the top and bottom surfaces of the shell is imposed. The shape functions used for the transversal displacement are C¹ conforming and the finite element is a conical/cylindrical panel with 8 nodes and 40 degrees of freedom. The model presented performs static analysis with arbitrary boundary conditions and loads, as well eigenvalue problems (free vibration and buckling). Illustrative examples are presented and discussed.     

Non‐linear exact geometry 12‐node solid‐shell element with three translational degrees of freedom per node

This paper presents the finite rotation exact geometry (EG) 12‐node solid‐shell element with 36 displacement degrees of freedom. The term ‘EG’ reflects the fact that coefficients of the first and second fundamental forms of the reference surface and Christoffel symbols are taken exactly at each element node. The finite element formulation developed is based on the 9‐parameter shell model by employing a new concept of sampling surfaces (S‐surfaces) inside the shell body. We introduce three S‐surfaces, namely, bottom, middle and top, and choose nine displacements of these surfaces as fundamental shell unknowns. Such choice allows one to represent the finite rotation higher order EG solid‐shell element formulation in a very compact form and to derive the strain–displacement relationships, which are objective, that is, invariant under arbitrarily large rigid‐body shell motions in convected curvilinear coordinates. The tangent stiffness matrix is evaluated by using 3D analytical integration and the explicit presentation of this matrix is given. The latter is unusual for the non‐linear EG shell element formulation. Copyright © 2011 John Wiley & Sons, Ltd.     

特征值问题迭代伽略金法与Rayleigh商加速

该文讨论特征值问题非协调有限元和混合有限元的加速计算方法。基于迭代伽略金法和Rayleigh商加速技巧,我们建立了特征值问题Wilson非协调有限元和Ciarlet-Raviart混合有限元的加速计算方案。这些新方案把在细网格上解一个特征值问题简化为在粗网格上解一个特征值问题和在细网格上解一个线性方程。文中证明了新方案的计算结果仍然保持了渐近最优精度阶,并用数值实验验证了理论结果。     

Adding transverse normal stresses to layered shell finite elements for the analysis of composite structures

This work presents a formulation developed to add capabilities for representing the through thickness distribution of the transverse normal stresses, σ_z, in first and higher order shear deformable shell elements within a finite element (FE) scheme. The formulation is developed within a displacement based shear deformation shell theory. Using the differential equilibrium equations for two-dimensional elasticity and the interlayer stress and strain continuity requirements, special treatment is developed for the transverse normal stresses, which are thus represented by a continuous piecewise cubic function. The implementation of this formulation requires only C⁰ continuity of the displacement functions regardless of whether it is added to a first or a higher order shell element. This makes the transverse normal stress treatment applicable to the most popular bilinear isoparametric 4-noded quadrilateral shell elements.

To assess the performance of the present approach it is included in the formulation of a recently developed third order shear deformable shell finite element. The element is added to the element library of the general nonlinear explicit dynamic FE code DYNA3D. Some illustrative problems are solved and results are presented and compared to other theoretical and numerical results.     

Energy and momentum conserving elastodynamics of a non‐linear brick‐type mixed finite shell element

The paper presents aspects of the finite element formulation of momentum and energy conserving algorithms for the non‐linear dynamic analysis of shell‐like structures. The key contribution is a detailed analysis of the implementation of a Simó–Tarnow‐type conservation scheme in a recently developed new mixed finite shell element. This continuum‐based shell element provides a well‐defined interface to strain‐driven constitutive stress updates algorithms. It is based on the classic brick‐type trilinear displacement element and is equipped with specific gradient‐type enhanced strain modes and shell‐typical assumed strain modifications. The excellent performance of the proposed dynamic shell formulation with respect to conservation properties and numerical stability behaviour is demonstrated by means of three representative numerical examples of elastodynamics which exhibit complex free motions of flexible structures undergoing large strains and large rigid‐body motions. Copyright © 2001 John Wiley & Sons, Ltd.