Dual approaches to finite element model updating

The problem of finding the least change adjustment to a stiffness matrix modeled by finite element method is considered in this paper. Desired stiffness matrix properties such as symmetry, sparsity, positive semidefiniteness, and satisfaction of the characteristic equation are imposed as side constraints of the constructed optimal matrix approximation for updating the stiffness matrix, which matches measured data better. The dual problems of the original constrained minimization are presented and solved by subgradient algorithms with different line search strategies. Some numerical results are included to illustrate the performance and application of the proposed methods.     

Fuzzy finite element model updating of the DLR AIRMOD test structure

This article presents the application of finite-element fuzzy model updating to the DLR AIRMOD structure. The proposed approach is initially demonstrated on a simulated mass-spring system with three degrees of freedom. Considering the effect of the assembly process on variability measurements, modal tests were carried out for the repeatedly disassembled and reassembled DLR AIRMOD structure. The histograms of the measured data attributed to the uncertainty of the structural components in terms of mass and stiffness are utilised to obtain the membership functions of the chosen fuzzy outputs and to determine the updated membership functions of the uncertain input parameters represented by fuzzy variables. In this regard, a fuzzy parameter is introduced to represent a set of interval parameters through the membership function, and a meta model (kriging, in this work) is used to speed up the updating. The use of non-probabilistic models, i.e. interval and fuzzy models, for updating models with uncertainties is often more practical when the large quantities of test data that are necessary for probabilistic model updating are unavailable.     

Theory and method for updating least-squares finite element model of symmetric generalized centro-symmetric matrices

Let X,F$X,F$ be a displacement matrix and load matrix, respectively. C (obtained by calculations or measurements) is an estimate matrix of the analytical model. A method is presented for correction of the model C, based on the theory of inverse problem of matrices. The corrected model is symmetric generalized centro-symmetric with specified displacements and loads, satisfying the mechanics characters of finite-element model. The application of the method is illustrated. It is more important that a perturbation analysis is given, which is not given in the earlier papers. Numerical results show that the method is feasible and effective.     

多孔介质中可压缩可混溶驱动问题的有限体积元法

有界区域上多孔介质中可压缩可混溶驱动问题由两个非线性抛物型方程耦合而成：压力方程和饱和度方程均是抛物型方程．运用有限体积元法对两个方程进行数值分析,给出了全离散有限体积元格式,并通过详细的理论分析,得到了近似解与原问题真解的最优H^1模误差估计。     

Finite element analysis for a nonlinear diffusion model in image processing

The optimal error estimate O(h^k+1) for a popular nonlinear diffusion model widely used in image processing is proved for the standard k^th-order (k ≥ 1) conforming tensor-product finite elements in the L²-norm. The optimal L²-estimate is obtained by the integral identity technique [1–3] without using the classic Nitsche duality argument [4].     

Finite element approximation of a periodic Ginzburg-Landau model for type-II superconductors

Summary We consider efficient finite element algorithms for the computational simulation of type-II superconductors. The algorithms are based on discretizations of a periodic Ginzburg-Landau model. Periodicity is defined with respect to a non-orthogonal lattice that is not necessarily aligned with the coordinate axes; also, the primary dependent variables employed in the model satisfy non-standard quasi-periodic boundary conditions. After introducing the model, we define finite element schemes, derive error estimates of optimal order, and present the results of some numerical calculations. For a similar quality of simulation, the resulting algorithms seem to be significantly less costly than are previously used numerical approximation methods.     

Finite element approximation of a prestressed shell model

This work deals with the finite element approximation of a prestressed shell model formulated in Cartesian coordinates system. The considered constrained variational problem is not necessarily positive. Moreover, because of the constraint, it cannot be discretized by conforming finite element methods. A penalized version of the model and its discretization are then proposed. We prove existence and uniqueness results of solutions for the continuous and discrete problems, and we derive optimal a priori error estimates. Numerical tests that validate and illustrate our approach are given.     

A constrained optimization based method for acoustic finite element model updating of cavities using pressure response

A method based on constrained optimization for updating of an acoustic finite element model using pressure response is proposed in this paper. The constrained optimization problem is solved using sequential quadratic programming algorithm. Updating parameters related to the properties of the sound absorbers and the measurement errors are considered. Effectiveness of the method is demonstrated by numerical studies on a 2D rectangular cavity and a car cavity. It is shown that the constrained formulation, that includes lower and upper bounds on the updating parameters in the form of inequality constraints, is important for obtaining a correct updated model. It is seen that the proposed updating method is not only able to effectively update the model to obtain a close match between the finite element model pressure response and the reference pressure response, but is also able to identify the correction factors to the parameters in error with reasonable accuracy.     

Finite element model of the impedance catheter technique for ventricular volume measurement

A two-dimensional, finite element model of the electric field occurring in the heart and surrounding tissues during left ventricular volume measurement by the impedance catheter technique has been developed. This model allows visualization of dynamic equipotential contours and current density distribution throughout different phases of the cardiac cycle. In addition, the model can be used to explore the effects of changes in the conductance of surrounding tissues as well as the actions of the right ventricle on left ventricular volume measurement. The model suggests a significant lack of homogeneity of the electric field inside the ventricular chamber. Furthermore, the field varies during different phases of the cardiac cycle. Alterations of parallel conductance can significantly affect the field inside the left ventricle. The right ventricular chamber is capable of drawing significant current and thus may affect left ventricular volume measurements. Knowledge of the relationship between chamber morphology and the imposed electric field may lead to more accurate impedance volume measurement techniques. The model is simple, apparently accurate, and requires little computing time.     

Finite element approximation of a nonlinear cross-diffusion population model

Summary. We consider a fully discrete finite element approximation of the nonlinear cross-diffusion population model: Find u_i, the population of the i^th species, i=1 and 2, such that where ji and g_i(u₁,u₂):=(_i–_iiu_i–_iju_j)u_i. In the above, the given data is as follows: v is an environmental potential, c_i, a_i are diffusion coefficients, b_i are transport coefficients, _i are the intrinsic growth rates, and _ii are intra-specific, whereas _ij, ij, are interspecific competition coefficients. In addition to showing well-posedness of our approximation, we prove convergence in space dimensions d3. Finally some numerical experiments in one space dimension are presented.Mathematics Subject Classification (2000): 65M60, 65M12, 35K55, 92D25Acknowledgements. Part of this work was carried out while the authors participated in the 2003 programme {\it Computational Challenges in Partial Differential Equations} at the Isaac Newton Institute, Cambridge, UK.     

Finite element analysis and application for a nonlinear diffusion model in image denoising

The stability analysis and error estimates are presented for a nonlinear diffusion model, which appears in image denoising and solved by a fully discrete time Galerkin method with kth (k ≥ 1) order conforming finite element spaces. Numerical experiments are provided with denoising several grayscale noisy images by our Galerkin method on bilinear finite elements. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 649–662, 2002; DOI 10.1002/num.10017     

Finite element method for viscoelastic fluids flows: approximation of the Maxwell model

We discuss about the finite element approximation of viscoelastic fluid flow. We consider a fluid obeying the Oldroyd model and particularly we study the purely viscoelastic case, the so-called Maxwell model, important in practice for the applications. We examine two kinds of methods used for the approximation of the Maxwell model: method using a splitting technique and finite element method satisfying inf-sup conditions relating tensor and velocity. We present numerical results for these methods and we discuss about their stability. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Finite element convergence for singular data

Convergence of the finite element solutionu ^h of the Dirichlet problem u= is proved, where is the Dirac -function (unit impulse). In two dimensions, the Green's function (fundamental solution)u lies outsideH ¹, but we are able to prove that . Since the singularity ofu is logarithmic, we conclude that in two dimensions the function log can be approximated inL ² near the origin by piecewise linear functions with an errorO (h). We also consider the Dirichlet problem u=f, wheref is piecewise smooth but discontinuous along some curve. In this case,u just fails to be inH ^5/2, but as with the approximation to the Green's function, we prove the full rate of convergence:u–u ^h ₁=O (h ^8/2) with, say, piecewise quadratics.     

Finite element approximations of a nonlinear diffusion model with memory

The convergence of a finite element scheme approximating a nonlinear system of integro-differential equations is proven. This system arises in mathematical modeling of the process of a magnetic field penetrating into a substance. Properties of existence, uniqueness and asymptotic behavior of the solutions are briefly described. The decay of the numerical solution is compared with both the theoretical and finite difference results.     

Finite element potentials

We present an explicit and efficient way for constructing finite elements with assigned gradient, or curl, or divergence. Some simple notions of homology theory and graph theory applied to the finite element mesh are basic tools for devising the solution algorithms.     

Finite element approximation

We study a discretization procedure, using finite elements, for a class of non-isotropic free-discontinuity problems based on the non-local approximation proposed in [18]