一种用于索结构分析的悬链线单元

根据弹性悬链线的理论解析解推导出适于索结构有限元分析的悬链线单元.与常用的三节点、五节点曲线单元相比,采用该单元编制的软件具有输入数据少、计算机时省、计算精度高的特点.     

Numerical solution of the problem of image formation in two planes using the integral fresnel formula

The paper considers numerical solution of the problem of constructing a plane phase optical element which forms prespecified images in two focal planes when illuminated by a plane wave front. The propagation of radiation from the optical element to the focal planes obeys the integral Fresnel formula. The phase function defining the optical element, is from the space of feasible phase functions satisfying both computational and technological constraints. The computer algorithm is based on efficient numerical solution of the inverse directional diagram problem. The phase function is constructed using the projector on the space of feasible phase functions. The efficiency of the phase optical element produced by the proposed procedure is demonstrated by solving a prototype problem involving simulation of the direct operator. The corresponding optical element was manufactured by electron-beam lithography. Physical experiments with a laser source illuminating the element confirmed the findings of the numerical experiment. Translated from Prikladnaya Matematika i Informatika, No. 1, pp. 86–89, 1999.     

On some versions of the element agglomeration AMGe method

The present paper deals with element‐based algebraic multigrid (AMGe) methods that target linear systems of equations coming from finite element discretizations of elliptic partial differential equations. The individual element information (element matrices and element topology) is the main input to construct the AMG hierarchy. We study a number of variants of the spectral agglomerate AMGe method. The core of the algorithms relies on element agglomeration utilizing the element topology (built recursively from fine to coarse levels). The actual selection of the coarse degrees of freedom is based on solving a large number of local eigenvalue problems. Additionally, we investigate strategies for adaptive AMG as well as multigrid cycles that are more expensive than the V‐cycle utilizing simple interpolation matrices and nested conjugate gradient (CG)‐based recursive calls between the levels. The presented algorithms are illustrated with an extensive set of experiments based on a matlab implementation of the methods. Copyright © 2008 John Wiley & Sons, Ltd.     

Determination of the element numbers of the regular polytopes

Let ?? be a set of n-dimensional polytopes. A set ?? of n-dimensional polytopes is said to be an element set for ?? if each polytope in ?? is the union of a finite number of polytopes in ?? identified along (n ? 1)-dimensional faces. The element number of the set ?? of polyhedra, denoted by e(??), is the minimum cardinality of the element sets for ??, where the minimum is taken over all possible element sets ${\Omega \in \mathcal{E}(\Sigma)}$ . It is proved in Theorem 1 that the element number of the convex regular 4-dimensional polytopes is 4, and in Theorem 2 that the element numbers of the convex regular n-dimensional polytopes is 3 for n ?? 5. The results in this paper together with our previous papers determine completely the element numbers of the convex regular n-dimensional polytopes for all n ?? 2.     

On the least-squares conjugate-gradient solution of the finite element approximation of Burgers' equation

A finite element approximation of the two-dimensional steady Burgers' equation is presented and a conjugate gradient approach is taken to solve the resulting finite element equations. The scheme is computationally efficient and is relatively easy to implement. An optimal error bound is established and a set of test problems with known analytic solutions is given to demonstrate the efficiency of the method.     

Analysis of the free-edge effect in composite laminates by the boundary finite element method

Because of the risk of delamination due to high interlaminar stresses in the vicinity of free edges of composite laminates, there is a strong interest in efficient methods for the analysis of this free-edge effect. By the example of a symmetric [0°/90°]_s cross-ply laminate, the Boundary Finite Element Method is presented as a very efficient numerical method, which combines the advantages of the finite element method and the boundary element method. Analogously to the boundary element method, only the boundary is discretized, while the element formulation is finite element based. The resultant stress field is shown to be in very good agreement qualitatively and quantitatively with the comparative finite element analysis. Submitted to the 11th International Conference on Mechanics of Composite Materials (Riga, June 11–15, 2000). Published in Mekhanika Kompozitnykh Materialov, Vol. 36, No. 3, pp. 355–366, March–April, 2000.     

On the h-adaptive coupling of FE and BE for viscoplastic and elasto-plastic interface problems

This paper presents a posteriori error estimates for the symmetric finite element and boundary element coupling for a nonlinear interface problem: A bounded body with a viscoplastic or plastic material behaviour is surrounded by an elastic body. The nonlinearity is treated by the finite element method while large parts of the linear elastic body are approximated using the boundary element method. Based on the a posteriori error estimates we derive an algorithm for the adaptive mesh refinement of the boundary elements and the finite elements. Its implementation is documented and numerical examples are included.     

A note on the least squares solution of the groundwater flow equation

The least squares finite element method is a member of the weighted residuals class of numerical methods for solving partial differential equations. The least squares finite element method is applied to the groundwater flow equation. Space is discretized with a C¹ continuous trial function and parameters are approximated with a C⁰ bilinear basis. Solutions for problems containing parameters with large localized spatial gradients are characterized by errors that are propagated throughout the entire domain. Second-order spatial convergence is observed, and extreme mesh refinement is required to match Galerkin and mixed least squares finite element results. Temporal discretization should be kept separate from the least squares spatial discretization. © 1994 John Wiley & Sons, Inc.     

Convergence and nonlinear stability of the Lagrange-Galerkin method for the Navier-Stokes equations

Summary The Lagrange-Galerkin method is a numerical technique for solving convection — dominated diffusion problems, based on combining a special discretisation of the Lagrangian material derivative along particle trajectories with a Galerkin finite element method. We present optimal error estimates for the Lagrange-Galerkin mixed finite element approximation of the Navier-Stokes equations in a velocity/pressure formulation. The method is shown to be nonlinearly stable.     

On the symmetric boundary element method and the symmetric coupling of boundary elements and finite elements

The finite element method and the boundary element method areamong the most frequently applied tools in the numerical treatmentof partial differential equations. However, their propertiesappear to be complementary: while the boundary element methodis appropriate for the most important linear partial differentialequations with constant coefficients in bounded or unboundeddomains, the finite element method seems to be more appropriatefor inhomogeneous or even nonlinear problems. but is somehowrestricted to bounded domains. The symmetric coupling of thetwo methods inherits the advantages of both methods. This paper treats the symmetric coupling of finite elementsand boundary elements for a model transmission problem in twoand three dimensions where we have two domains: a bounded domainwith nonlinear, even plastic material behaviour, is surroundedby an unbounded, exterior, domain with isotropic homogeneouslinear elastic material. Practically. the coupling is performedsuch that the boundary element method contributes a macro-element,like a large finite element, within a standard finite elementanalysis program. Emphasis is on two-dimensional problems wherethe approach using the Poincaré-Steklov operator seemsto be impossible at first glance. E-mail: cc{at}numerik.uni-kiel.de     

Mixed Finite Element Methods for the Ferrofluid Model with Magnetization Paralleled to the Magnetic Field

Mixed finite element methods are considered for a ferrofluid flow model with magnetization paralleled to the magnetic field. The ferrofluid model is a coupled system of the Maxwell equations and the incompressible Navier-Stokes equations. By skillfully introducing some new variables, the model is rewritten as several decoupled subsystems that can be solved independently. Mixed finite element formulations are given to discretize the decoupled systems with proper finite element spaces. Existence and uniqueness of the mixed finite element solutions are shown, and optimal order error estimates are obtained under some reasonable assumptions. Numerical experiments confirm the theoretical results.     

Research on the Fictitious Source Points of the Hybrid Fundamental Solution⁃Based Finite Element Method for Heat Conduction Problems北大核心CSCD

针对热传导问题,提出了杂交基本解有限元法.首先,假设两个独立场:一个为利用基本解线性组合近似的单元域内温度场,另一个为使用与传统有限元法相同形式的辅助网线温度场.然后,利用修正变分泛函将上述两个独立场关联起来,并导出有限元列式.然而,该方法的准确性很大程度上取决于源点的分布和数量,通常将源点布置在单元外部两种虚拟边界上:与单元相似的边界和圆形边界.此外,还提出了双重虚拟边界,并与上述两种源点布局方式进行对比.通过两个典型数值算例,验证了该文方法在不同源点布局下的有效性和对网格畸变的不敏感性.     

Twisted conjugacy classes of the unit element

We study twisted conjugacy classes of the unit element in different groups. Fel’shtyn and Troitsky showed that the twisted conjugacy class of the unit element of an abelian group is a subgroup for every automorphism. The structure is investigated of a group whose twisted conjugacy class of the unit element is a subgroup for every automorphism (inner automorphism).     

不可压缩流问题低次元稳定有限体积数值方法研究

分析了R^d,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈Hd,d=2,3维不可压缩流Stokes问题低次元稳定有限体积方法,它主要利用局部压力投影方法对两种流行但不满足inf-sup条件的有限元配对(P_1-P_0和P_1-P_1)在有限体积方法的框架下进行稳定;利用有限元与有限体积方法的等价性进行有限体积方法理论分析.结果表明不可压缩流Stokes问题在f∈H¹情况下,本文方法得到的解与稳定有限元方法解之间具有O(h1情况下,本文方法得到的解与稳定有限元方法解之间具有O(h²)阶超收敛阶结果,且稳定有限体积方法取得了与稳定有限元方法相同的收敛速度,与稳定有限元方法比较,稳定有限体积方法计算简单高效,同时保持物理守恒,因此在实际应用中具有很好的潜力。     

非均匀Reissner板弯曲的精确元法

本文在阶梯折算法和精确解析法的基础上,提出构造有限元的新方法——精确元法.该方法不用变分原理,可适用于任意变系数正定和非正定偏微分方程.利用该方法,得到Reissuer板弯曲的一个非协调单元,它具有十五个自由度.由于节点位移参数仅含有挠度和转角,因此处理任意边界条件非常容易.文中给出证明,位移和内力均收敛于精确解.由精确元法所得到的单元不仅能用于厚板,也可用于薄板.文末给出四个算例.算例表明,利用本文的方法,可获得满意的结果,并有较高的数值精度.     

一种改进边界元法在弹性扭转问题中的应用

本文用一种改进边界元法分析与计算了椭圆截面等直杆的扭转问题.并与边界元法的解进行比较,其结果极为符合.然而,改进边界元法较边界元法所需要的数据量少得多,计算时间也将大大减少了.因此,本文方法对求解Poisson方程问题是一种经济而行之有效的数值计算方法.     

求收益模糊的贝叶斯纳什均衡解的结构元方法

主要目的是利用结构元方法求解收益模糊的贝叶斯纳什均衡.首先,在原有结构元理论基础上,给出了多元模糊值函数的定义及其结构元表示;其次,给出了在混合策略下,收益模糊的贝叶斯纳什均衡的定义,并证明了其存在性定理;然后,利用结构元理论,将该博弈模型等价地转化为一个经典的博弈模型,简化了原问题的求解.最后的应用实例说明了该方法的有效性.     

A two-grid finite element method for the Allen-Cahn equation with the logarithmic potential

In this paper, we present a two-grid finite element method for the Allen-Cahn equation with the logarithmic potential. This method consists of two steps. In the first step, based on a fully implicit finite element method, the Allen-Cahn equation is solved on a coarse grid with mesh size H. In the second step, a linearized system whose nonlinear term is replaced by the value of the first step is solved on a fine grid with mesh size h. We give the energy stabilities of the traditional finite element method and the two-grid finite element method. The optimal convergence order of the two-grid finite element method in H¹ norm is achieved when the mesh sizes satisfy h = O(H²). Numerical examples are given to demonstrate the validity of the proposed scheme. The results show that the two-grid method can save the CPU time while keeping the same convergence rate.     

模糊运算和模糊有限元静力控制方程的求解

根据模糊数的区间形式表达和区间运算的性质,给出了模糊数和模糊变量的运算规则.据此并依据区间有限元理论,提出了结构模糊有限元静力控制方程的几种求解方法.方法可根据输入模糊数的隶属函数,给出结构响应量的可能性分布.且计算量小,易于实施.算例分析说明了方法是实用和可行的.     

Formulation and evaluation of a finite element model for the linear analysis of stiffened composite cylindrical panels

A finite element model for linear static and free vibration analysis of composite cylindrical panels with composite stiffeners is presented. The proposed model is based on a cylindrical shell finite element, which uses a first-roder shear deformation theory. The stiffeners are curved beam elements based on Timoshenko and Saint-Venant assumptions for bending and torsion respectively. The two elements are developed in a cylindrical coordinate system and their stiffness matrices result from a hybrid-mixed formulation where the element assumed stress field is such that exact equilibrium equations are satisfied. The elements are free of membrane and shear locking with correct satisfaction of rigid body motions. Several examples dealing with stiffened isotropic and laminated plates and shells with eccentric as well as concentric stiffeners are analyzed showing the validity of the models.