基于Kriging模型和提升小波变换的随机模型修正北大核心CSCD

为提高随机模型修正效率,减小计算量,提出了一种基于Kriging模型和提升小波变换的随机模型修正方法.首先,对加速度频响函数进行提升小波变换,提取第5层近似系数代替原频响函数.其次,采用拉丁超立方抽样抽取待修正样本,将其作为Kriging模型的输入,对应的近似系数作为输出,构建Kriging模型.提出了一种引入莱维飞行(Lévy flight)的蝴蝶优化算法(LBOA),并将其应用于Kriging模型相关参数的寻优中,提高Kriging模型的精度.最后,以最小化Wasserstein距离为目标,通过鲸鱼优化算法求解待修正参数的均值.测试函数结果表明,LBOA在寻优能力、收敛精度和稳定性等方面有了很大的提升.数值算例的修正误差均低于0.4%,验证了所提模型修正方法具有较高的修正精度和效率.     

Model updating for undamped gyroscopic systems with connectivity constraints

ABSTRACT

An important and difficult aspect for the finite element model updating problem is to make the updated model have physical meaning, that is, the connectivity of the original model should be preserved in the updated model. In many practical applications, the system matrices generated by discretization of a distributed parameter system with the finite element techniques are often very large and sparse and are of some special structures, such as symmetric and band structure (diagonal, tridiagonal, pentadiagonal, seven-diagonal, etc.). In this paper, the model updating problem for undamped gyroscopic systems with connectivity constraints is considered. The method proposed not only preserves the connectivity of the original model, but also can update the analytical matrices with different bandwidths, which can meet the needs of different structural dynamic model updating problems. Numerical results illustrate the efficiency of the proposed method.     

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

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

A beam segment element for dynamic analysis of large aqueducts

Large aqueduct structure is a complex structure that is commonly used in hydraulic engineering, especially in large-scale water conveying projects. The analysis of dynamic response for an aqueduct structure is necessary if the aqueduct is built in an earthquake area. Traditional 3D finite element analysis is time consuming and the existing simplified response method cannot take into account all the effects, such as the bending-torsion coupling effect and the constrained torsion, of the deformations of the thin wall structure of the aqueduct body. For this special structure, a simple and yet accurate model for dynamic analysis is needed. In this paper, a beam segment element is developed and used for the calculation of dynamic response for aqueduct structures. With the frame of the aqueduct being modeled using beam element, the proposed model can calculate the dynamic response of the whole aqueduct structures. Results are compared with that of a general purpose finite element analysis software using 3D finite element model. Good agreement is achieved between the two models. However, the proposed model needs less elements and much less computing time.     

An iterative method for updating gyroscopic systems based on measured modal data

An efficient iterative method for updating the mass, gyroscopic and stiffness matrices simultaneously using a few of complex measured modal data is developed. By using the proposed iterative method, the unique symmetric solution can be obtained within finite iteration steps in the absence of roundoff errors by choosing a special kind of initial matrices. Numerical results show that the presented method can be used to update finite element models to get better agreement between analytical and experimental modal parameters.     

A model updating approach based on design points for unknown structural parameters

Finite element analysis has become an essential tool to estimate structural responses under static and dynamic loads. However, there are a lot of uncertainties in structural properties. For this reason, in many cases, the outcomes of the theoretical and experimental modal analyses do not match. Therefore, the analytical models of the structures need to be updated according to the experimental test results. The commonly used method to get parameters for model updating is experimental modal analysis which provides structural dynamic characteristic (natural frequencies, mode shapes and modal damping ratio). There are many methods available for the updating process. This study addresses an updating algorithm to modify the numerical models by using the design points for unknown structural properties. The proposed method aims to minimize the difference between the analytical and experimental natural frequencies by updating uncertain parameters for each mode and combine them to get an optimum solution. The algorithm is tested on a column and a 2D frame models. These models are investigated by taking the connection rigidity and elasticity modulus as unknown parameters. It is observed that the proposed algorithm gives better results for unknown structural properties compared to the initial values.     

基于等效梁模型的运载火箭动力学特性仿真预示研究

针对运载火箭动力学特性仿真预示问题,提出了一种基于等效梁模型,其是可考虑贮箱内液体推进剂影响高效、准确的仿真预示方法.应用等效厚度法建立了运载火箭的等效梁模型,并将运载火箭精细有限元模型的模态分析结果作为目标,通过有限元模型修正提升了等效梁模型的精度.以某型号运载火箭为算例,应用该文方法,建立了其等效梁模型,并基于集中质量法和耦合质量法两种液体推进剂等效方法对其进行了考虑贮箱内液体推进剂影响的动力学特性预测,比较了两种方法的分析结果.     

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.     

有限元模型修正中的最佳矩阵逼近

1引言 在飞行器、船舶、桥梁等结构设计中,要定量、准确地进行结构动力学分析,解决飞行器、船舶、桥梁等工程结构中普遍存在的振动问题,首先必须建立结构的动力学模型.     

求解流固耦合问题的一种四步分裂有限元算法

基于arbitrary Lagrangian Eulerian (ALE) 有限元方法,发展了一种求解流固耦合问题的弱耦合算法．将半隐式四步分裂有限元格式推广至求解ALE描述下的Navier-Stokes（N-S）方程,并在动量方程中引入迎风流线（streamline upwind/Petrov-Galerkin, SUPG）稳定项以消除对流引发的速度场数值振荡;采用Newmark-β法对结构方程进行时间离散;运用经典的Galerkin有限元法求解修正的Laplace方程以实现网格更新,每个计算步施加网格总变形量防止结构长时间、大位移运动时的网格质量恶化．运用上述算法对弹性支撑刚性圆柱体的流致振动问题进行了数值模拟,计算结果与已有结果相吻合,初步验证了该算法的正确性和有效性．     

A gradient based iterative algorithm for solving model updating problems of gyroscopic systems

Model updating should be correlated with experimental data to ensure that it models the dynamics of the real structure and the updated model predicts the dynamics of a structure more accurately. Considering that the iterative methods for model updating have aroused little public attention, this paper studies an iterative algorithm for quadratic model updating problems which can incorporate the measured modal data into the finite element model to produce an adjusted finite element model on the mass, gyroscopic and stiffness matrices that closely match the experimental modal data. By this method, the best approximation symmetric and skew-symmetric solutions can be obtained by choosing a convergence factor. Numerical example shows that the introduced iterative algorithm is quite efficient.     

Adaptive surrogate-based design optimization with expected improvement used as infill criterion

A prominent advantage of using surrogate models in structural design optimization is that computational effort can be greatly reduced without significantly compromising model accuracy. The essential goal is to perform the design optimization with fewer evaluations of the typically finite element analysis and ensuring accuracy of the optimization results. An adaptive surrogate based design optimization framework is proposed, in which Latin hypercube sampling and Kriging are used to build surrogate models. Accuracy of the models is improved adaptively using an infill criterion called expected improvement (EI). It is the anticipated improvement that an interpolation point will lead to the current surrogate models. The point that will lead to the maximum EI is searched and used as infill points at each iteration. For constrained optimization problems, the surrogate of constraint is also utilized to form a constrained EI as the corresponding infill criterion. Computational trials on mathematical test functions and on a three-dimensional aircraft wing model are carried out to test the feasibility of this method. Compared with the traditional surrogate base design optimization and direct optimization methods, this method can find the optimum design with fewer evaluations of the original system model and maintain good accuracy.     

多孔介质驱动问题的混合元最小二乘特征有限元方法及其收敛分析

１．引言多孔介质二相驱动问题的数学模型是偶合的非线性偏微分方程组的初边值问题．该问题可转化为压力方程和浓度方程［１－４］．浓度方程一般是对流占优的对流扩散方程,它的对流速度依赖于比浓度方程的扩散系数大得多的Ｆａｒｃｙ速度．因此Ｄａｒｃｙ速度的求解精度直接影响着浓度的求解精度．为了提高速度的求解精度,７０年代Ｐ．Ａ．Ｒａｖｉａｔ和Ｊ．Ｍ．Ｔｈｏｍａｓ提出混合有限元方法［５］．Ｊ．ＤｏｕｇｌａｓＪｒ,Ｔ．Ｆ．Ｒｕｓｓｅｌｌ,Ｒ．Ｅ．Ｅｗｉｎｇ,Ｍ．Ｆ．Ｗｈｅｅｌｅｒ［１］－［４］,［９］,［１２］袁…     

The direct updating of damping and gyroscopic matrices

In this paper we develop an efficient numerical method for the finite element model updating of damped gyroscopic systems. This model updating of damped gyroscopic systems is proposed to incorporate the measured modal data into the finite element model to produce an adjusted finite element model on the damping and gyroscopic matrices that closely match the experimental modal data.     

MULTI-SCALE METHODS FOR INVERSE MODELING IN 1-D MOS CAPACITOR

In this paper,we investigate multi-scale methods for the inverse modeling in 1-D Metal-Oxide-Silicon(MOS) capactior,First,the mathematical model of the device is given and the numerical simulation for the forward problem of the model is implemented using finite element method with adaptive moving mesh. Then numerical analysis of these parameters in the model for the inverse problems is presented .Some matrix analysis tools are applied to explore the parameters‘ sensitivities,And thired,the parameters are extracted using Levenberg-Marquardt optimization method.The essential difficulty arises from the effect of multi-scale physical differeence of the parameters.We explore the relationship between the parameters‘ sensitivitites and the sequencs for optimization,which can seriously affect the final inverse modeling results.An optimal sequence can efficiently overcome the multip-scale problem of these parameters,Numerical experiments show the efficiency of the proposed methods.     

Two-phase damage detection of beam structures under moving load using multi-scale wavelet signal processing and wavelet finite element model

This paper proposes a damage detection method with two phases, namely, localization and quantification, for beam structures subjected to moving load and successfully validates it via a laboratory experiment. Firstly, the discrete wavelet transform (DWT) is applied to decompose the displacement response change induced by a moving vehicle and locate potential structural damages. Then adaptive-scale wavelet finite element model (WFEM) updating is employed to estimate the damage severity in the identified damage regions in a progressive fashion. The elemental scales of WFEM are adaptively changed according to not only the moving vehicle position but also the progressively identified damage regions. Such a method can effectively minimize the number of modeling degree-of-freedoms (DOFs) and updating parameters during optimization. The feasibility and effectiveness of the proposed method is examined through a laboratory experiment with different damage scenarios. The results indicate the proposed method can achieve good consistency between structural modeling, damage scenarios, and load conditions, as well as an optimal tradeoff between damage detection accuracy and efficiency.     

Optimum design of laminated composite plates under dynamic excitation

An integrated model for optimum weight design of symmetrically laminated composite plates subjected to dynamic excitation is presented in this work. Optimum design procedure based on flexibility and strength criteria is presented. The objective is to determine the optimum thicknesses of the laminate layers and its optimum orientations without exhibiting any failure according to Tsai-Wu failure criterion. The finite element method, based on Mindlin plate theory, is used in conjunction with an optimization method in order to determine the optimum design. Newmark algorithm, as an implicit time integration scheme, is used to discretize the time domain and calculate the transient response of the laminated composite plate. Exterior penalty method is exploited for the constrained minimization procedure. Fletcher-Powell algorithm is used for the unconstrained minimization process. To verify the capability and efficiency of the proposed model, three examples are solved. The examples deal with flexibility and stress constraints for different boundary conditions under various dynamic excitations.     

Advances in concrete arch dams shape optimization

This paper presents an efficient methodology to find the optimum shape of arch dams. In order to create the geometry of arch dams a new algorithm based on Hermit Splines is proposed. A finite element based shape sensitivity analysis for design-dependent loadings involving body force, hydrostatic pressure and earthquake loadings is implemented. The sensitivity analysis is performed using the concept of mesh design velocity. In order to consider the practical requirements in the optimization model such as construction stages, many geometrical and behavioral constrains are included in the model in comparison with previous researches. The optimization problem is solved via the sequential quadratic programming (SQP) method. The proposed methods are applied successfully to an Iranian arch dam, and good results are achieved. By using such methodology, efficient software for shape optimization of concrete arch dams for practical and reliable design now is available.     

A coupling of nonconforming and mixed finite element methods for Biot's consolidation model

In this article, we develop a nonconforming mixed finite element method to solve Biot's consolidation model. In particular, this work has been motivated to overcome nonphysical oscillations in the pressure variable, which is known as locking in poroelasticity. The method is based on a coupling of a nonconforming finite element method for the displacement of the solid phase with a standard mixed finite element method for the pressure and velocity of the fluid phase. The discrete Korn's inequality has been achieved by adding a jump term to the discrete variational formulation. We prove a rigorous proof of a‐priori error estimates for both semidiscrete and fully‐discrete schemes. Optimal error estimates have been derived. In particular, optimality in the pressure, measured in different norms, has been proved for both cases when the constrained specific storage coefficient c₀ is strictly positive and when c₀ is nonnegative. Numerical results illustrate the accuracy of the method and also show the effectiveness of the method to overcome the nonphysical pressure oscillations. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

一种高精度的裂纹奇异单元

基于广义伽辽金法的多变量有限元算法,增加了连续体力学有限元模型建立的灵活性.本文利用它,通过数值试验的对比建立了一种高精度的含奇异性的裂纹单元,并对多变量奇异元的构成进行了探讨.