强Guass-seidel法求解有限元模型修正问题

有限元模型修正是一类特殊的二次反特征值问题．我们将有限元模型修正看成二次规划问题来解决,并采用非线性Gauss－Seidel方法来求解其相应的Lagrange对偶函数．最后,给山的数值文验说明方法的有效性．     

有限元模型修正问题的不精确最速下降迭代解

在结构动力分析中,往往需利用结构振动测试所得的实际测量数据(如振动频率和振型),对结构分析模型进行最优修正,使之更能合理反映结构的实际性能,其实质即为计算数学中的特征值反问题.本文考虑有阻尼结构振动中的-类反问题,用一组不完备的模态测量数据修正系统质量矩阵、刚度矩阵和阻尼矩阵,通过等价正交投影思想将原问题转化成-个闭凸锥上的正交投影问题,构造-个不精确最速下降迭代法求解,并讨论了收敛性.算例表明算法是有效的.     

交替方向法在保结构的有限元模型修正问题中的应用

本文应用交替方向法处理二次特征值反问题,并要求同时保持对称性,半正定性和稀疏性等结构性质.实验表明我们提出的方法是可行的.     

Model Reduction of Large Scale Finite Element Models

The usage of flexible bodies in Multibody simulations (MBS) has widely increased their industrial application. Finite Element (FE) models with up to 10 million degrees of freedom (DOF) and some hundred million nonzero matrix entries are used to describe the flexible bodies. Before such a model can be included into an MBS, the number of DOF must be reduced to an appropriate size. Using modal reduction often the critical issue arises which modes to choose while Component Mode Synthesis based methods often lead to a relatively big size of the resulting model. Alternative methods using Moment Matching and Balanced Truncation can result in a smaller size while still remaining accurate enough. Sometimes these matrices are so huge that they can not even be stored in one computers main memory. The calculation of the necessary orthogonal Krylov subspaces needs an LU factorization which is also very memory intensive. To meet these requirements, distributed computation is used which also shortens the computational time of the reduced process. In this work, an industrial relevant FE model is reduced to a much smaller size using alternative methods. Accuracy is verified by comparing the frequency response in a defined frequency range of the original and the reduced model. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Finite Element Methods for a Model for Full Waveform Acoustic Logging

A model is defined to simulate the propagation of waves in aradially symmetric, isotropic, composite system consisting ofa fluid-filled well bore ^f through a fluid-saturated poroussolid ^p. Biot's equations of motion are chosen to describe thepropagation of waves in ^p, while the standard equation of motionfor compressible inviscid fluids is used for ^f, with appropriateboundary conditions at the contact surface between ^f and ^p.Also, absorbing boundary conditions for the artificial boundariesof ^p are derived for the model, their effect being to make themtransparent for waves arriving normally First, results on the existence and uniqueness of the solutionof the differential problem are given and then a discrete-time,explicit finite element procedure is defined and analysed, withfinite element spaces suited for radially symmetric problemsbeing used for the spatial discretisation.     

利用静力试验数据修正具有广义中心对称有限元模型

设X,B分别是测得的位移矩阵和载荷矩阵,C是有限元方法得到的理论模型的估计,找广义中心对称的矩阵（A）使得（A）X=B,且是Frobenius范数意义下C的最佳逼近.并给出了解（A）的扰动分析,数值结果表明该方法是行之有效的.     

On the Wellbore Contact of Drill Strings in a Finite Element Model

In drill string dynamics the Finite Element Method is usually applied to models of very long drill strings in a wellbore with arbitrary curvature. Taking account of geometrical constraint between the drill string and the wellbore, a high density of nodes is necessary. This density is much higher than the one needed to describe the natural vibrations properly, so this firstly leads to an extension of the computing time. A penalty function is frequently utilized to describe the contact problem between the drill string and the wellbore where the contact normal force acts only on the nodal points of the drill string. It was recognized that only node-to-surface contact models cannot fulfill this geometrical constraint, because the segment between two nodal points deeply penetrates the wellbore wall in some cases. A process with Gaussian points along the segment in time domain will be introduced, so that the drill string will be described according to this geometrical constraint with good accuracy but with a smaller density of nodes and less computing time. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Neumann - Dirchlet coupling strategy of Finite Element and Boundary Element method in elastodynamics

The Finite Element Method (FEM) and the Boundary Element Method (BEM) are the most used numerical tools for solid mechanics analysis. Each one of these methods has advantages and drawbacks in different cases. In order to take advantage of both methods, a nonoverlapping domain decomposition method FEM - BEM in elastodynamics is presented. The domain is divided in two subdomains and each one of them is analyzed separately and only the interface information is exchanged. An iterative Neumann - Dirchlet algorithm with relaxation is used, to get continuity and the equilibrium conditions at the interface. The FEM time integration is carried out using the Newmark's method and the BEM approach in time domain is based in the Convolution Quadrature Method developed by Lubich. Numerical examples are presented to show agreement with other available numerical results. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

细菌模型的非协调有限元分析

将非协调元应用于描述细菌传播的反应扩散方程组的初边值问题.借助单元的一些特性和非协调误差估计技巧,分别在半离散和全离散有限元格式下,研究了其数值解与精确解的误差估计,得到了最优的误差估计以及超逼近结果.     

A Positivity-preserving Finite Element Method for Quantum Drift-diffusion Model

In this paper, we propose a positivity-preserving finite element method for solving the three-dimensional quantum drift-diffusion model. The model consists of five nonlinear elliptic equations, and two of them describe quantum corrections for quasi-Fermi levels. We propose an interpolated-exponential finite element (IEFE) method for solving the two quantum-correction equations. The IEFE method always yields positive carrier densities and preserves the positivity of second-order differential operators in the Newton linearization of quantum-correction equations. Moreover, we solve the two continuity equations with the edge-averaged finite element (EAFE) method to reduce numerical oscillations of quasi-Fermi levels. The Poisson equation of electrical potential is solved with standard Lagrangian finite elements. We prove the existence of solution to the nonlinear discrete problem by using a fixed-point iteration and solving the minimum problem of a new discrete functional. A Newton method is proposed to solve the nonlinear discrete problem. Numerical experiments for a three-dimensional nano-scale FinFET device show that the Newton method is robust for source-to-gate bias voltages up to 9V and source-to-drain bias voltages up to 10V.     

A Verified Monotonicity-Based Solution of a Simple Finite Element Model with Uncertain Node Locations

A tight verified solution enclosure is obtained for the node displacements of a simple truss model, whose parameters, including the node locations, are uncertain. The solution is based on a monotonicity analysis of these interval parameters. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical Model for Electro-Mechanical Analysis Based on Finite Element Method

In this paper, we present a numerical procedure that can be used to model the electro-mechanical coupled behavior of the dielectric actuator domain. The equation describing the electrostatical part is given by the reduced form of the Maxwell equation and the electrostatic potential [1]. The mechanical problem is described by the constitutive equations and equilibrium equations. Using the finite element method, this technique is to divide a whole problem into sub-problems. The complexity of the original problem is therefore reduced by focusing only on the most relevant areas. A finite element analysis is then performed by applying the electrostatic Maxwell pressure as Neumann boundary conditions to compute the displacements. Once the displacement is computed, the electrostatic domain or the conductor is updated. Electrostatic analysis is performed on the updated geometry and the finite element method is then used to determine the change in potential due to geometric perturbations. Once the surface charge densities are known, the new electrostatic Maxwell pressure is computed. The mechanical and electrostatic analysis is repeated until an equilibrium state is computed. The procedure is demonstrated in the paper by the solution of some two-dimensional and three-dimensional problems. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

更一般的杂交广义变分原理及有限元模型

本文根据钱伟长教授提出的更一般的广义变分原理,给出了适用于有限元法中的更广义杂交变分原理,并由此建立了新的广义杂交模理. 进一步以变厚度薄板弯曲单元为例,对基于各种不同的广义变分原理建立的各种杂交元做了比较.     

A Variational Finite Element Model for Large-Eddy Simulations of Turbulent Flows

The authors introduce a new Large Eddy Simulation model in a channel,based on the projection on finite element spaces as filtering operation in its variational form,for a given triangulation{Th}h>0.The eddy viscosity is expressed in terms of the friction velocity in the boundary layer due to the wall,and is of a standard sub grid-model form outside the boundary layer.The mixing length scale is locally equal to the grid size.The computational domain is the channel without the linear sub-layer of the boundary layer.The no-slip boundary condition(or BC for short)is replaced by a Navier(BC)at the computational wall.Considering the steady state case,the authors show that the variational finite element model they have introduced,has a solution(vh,ph)h>0that converges to a solution of the steady state Navier-Stokes equation with Navier BC.     

Pseudoelasticity of Shape Memory Alloys - a 1D Material Model and Finite Element Implementation

This contribution is concerned with the formulation of a 1D-constitutive model accounting for the pseudoelastic behavior of shape memory alloys. The stress-strain-relationship is idealized by a hysteresis both in the compression as in the tension loading range. It is characterized by an upper loading path, which is to be ascribed to the transformation of the lattice to a martensitic structure. Unloading the material, a lower path is described, because of the reverse transformation into austenitic lattice. The constitutive model is based on a switching criterion which serves as a potential function for the evolution of the internal state variables. The model distinguishes between local and global variables to describe the hysteresis effects for the compression and tension range. A strain driven algorithm which captures the complete nonlinear material behavior is presented. The boundary value problem is solved for a truss element applying the finite element method. A consistent linearization of the nonlinear equations is derived. Simple examples will demonstrate the applicability of the proposed model. For future developments the usage of shape memory alloys within civil engineering structures is aimed. The advantage of the material is the very good damping behavior and the potential to overcome great strains. Both properties are distinguished to be of engineering interest. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Parameter Identification of a Finite Viscoplastic Material Model using different Deterministic Optimization Methods

The efficiency of two deterministic optimization strategies are discussed in order to solve the arising least squares minimization problem. A sequential quadratic programming method is compared with a trust region enhanced simplex method. In particular we focus on the effort to determine the necessary gradients for the sqp‐method and compare the computational costs of the identification process. In the later development the identification problem is extended to account for inhomogeneous deformations using the finite element method for simulation purpose. Here, displacement fields are determined via an optical measurement technique. The accompanying modification of the solution procedure results in an extended least squares functional. The sensitivity analysis, needed for the gradient based method, is embedded in a distinct format in the finite element context.     

页岩气藏渗流模型的特征混合有限元数值模拟

将特征有限元方法和混合有限元方法进行耦合,对页岩气藏渗流模型进行了数值模拟,给出了详细的误差分析,得到了最优的L~2模误差估计,并用数值实验验证了方法的有效性.     

Optimality of a Standard Adaptive Finite Element Method

In this paper an adaptive finite element method is constructed for solving elliptic equations that has optimal computational complexity. Whenever, for some s > 0, the solution can be approximated within a tolerance ε > 0 in energy norm by a continuous piecewise linear function on some partition with O(ε^-1/s) triangles, and one knows how to approximate the right-hand side in the dual norm with the same rate with piecewise constants, then the adaptive method produces approximations that converge with this rate, taking a number of operations that is of the order of the number of triangles in the output partition. The method is similar in spirit to that from by Morin, Nochetto, and Siebert, and so in particular it does not rely on a recurrent coarsening of the partitions. Although the Poisson equation in two dimensions with piecewise linear approximation is considered, the results generalize in several respects.     

A Characteristic Finite Volume Element Method for the Air Pollution Model

In this article, a characteristic finite volume element method is presented for solving air pollution models. The convection term is discretized using the characteristic method and diffusion term is approximated by finite volume element method. Compared with standard finite volume element method, our proposed method is more accurate and efficient, especially suitable to solve convection-dominated problems. The proposed numerical schemes are analyzed for convergence in L ₂ norm. Some numerical results are presented to demonstrate the efficiency and accuracy of the method.