一种基于加速度频响函数的动力学模型修正方法

提出了一种基于加速度频响函数的动力学模型修正方法,该方法应用试验测量和计算所得的加速度频响函数矩阵对动力学模型的刚度矩阵、质量矩阵、阻尼系数进行修正,使修正后模型的加速度频响函数与试验测量所得的相一致.该法具有明确的物理意义,数值算例检验了方法的有效性和精度.     

基于径向基函数的频响函数模型修正

针对使用频响函数进行有限元模型修正存在的时间成本和精度问题,结合模态参与变异系数法和模态动能法分别优化激励点和测点,使获得的模态信息更完整;然后,引入径向基函数(RBF)模型减少原有限元模型计算次数,并根据均方根误差准则对所构建代理模型的参数(spread)进行优选,提高模型预测精度;最后,选用一个36自由度的二维桁架模型进行可行性验证,对比有限元法、Kriging模型及二阶响应面模型的修正精度和迭代时间,结果表明,本文方法具有较好的优势。     

基于频响函数提升小波总能量的桁架结构随机模型修正

为了解决模型修正问题中的随机性,构建了一种基于提升小波总能量的随机模型修正方法.首先,将结 构的加速度频响函数进行提升小波变换,并提取提升小波总能量来代替加速度频响函数;然后,以待修正参数作为输入,提升小波总能量为输出构建响应面代理模型代替原来的有限元模型;接着,运用蒙特卡洛抽样抽取响应样本,并设定阈值筛选响应样本;最后,以代理模型预测得到的响应和抽样所得真实响应之间的差值最小为 目标函数,通过布谷鸟优化算法寻优求解待修正参数的均值.算例表明,所提方法修正后参数的最大误差小于3.3％,相应的频响函数曲线重合度高.     

基于Kriging模型的频响函数有限元模型修正方法

针对使用频响函数进行有限元模型修正的问题,提出了一种基于Kriging模型的修正方法,用于检测结构由损伤引起的在单元刚度特性上的衰减。本文方法可以在不需要推导修正参数与频响函数残差代数关系的前提下,通过少量测点提供的有效数据快速求解;还可以通过控制算法的终止准则来提高对未知区域的探索程度,降低结果收敛到局部解上的可能。使用Kriging模型可以有效地减少原有限元模型的计算次数,保证计算效率的同时,为对结构进行更准确精密的有限元建模提供了便利。     

基于Kriging模型和改进MCMC算法的随机有限元模型修正

针对待修正参数维数较高时,标准马尔可夫链蒙特卡罗MCMC (Markov Chain Monte Carlo)算法不易收敛、拒绝率高的问题,提出了基于Kriging模型和在MCMC中融合花朵授粉算法的修正方法.首先,以待修正参数作为输入,以应变模态作为输出,建立Kriging模型,通过蝙蝠算法确定Kriging模型的相关系数;然后,采用最大熵的贝叶斯方法估计参数的后验概率密度函数,将花朵授粉算法融入MH (M etropolis-Hasting)抽样算法,提高局部寻优和全局寻优能力;最后,通过三自由度弹簧-质量系统和三维桁架结构的数值算例验证所提模型修正方法,修正后参数相对误差均低于0.86％.结果 表明,所提方法修正后较高维参数的马尔可夫链能够快速收敛且样本接受率也有所提高,该方法也对随机噪声具有一定的鲁棒性.     

利用频响函数对阻尼结构进行模型修正的方法

基于频响函数和模型缩聚技术提出了一个动力结构模型修正方法。与传统频响函数法相比,由于引入了阻尼刚度比,减少了修正参数,提高了计算效率。并且针对传统频响函数法测试维数过高的缺点,引入模型缩聚技术,降低了模型测量维数的要求。数值算例证明了本方法的有效性与可行性。本方法可以对有阻尼的模型进行修正,克服了模态修正法的缺点;利用模型缩聚减小了测试点的数目,引进单元阻尼刚度比。与传统的频响函数修正法相比,减少了修正未知数,提高计算效率,并且修正结果也较准确。     

有限元模型阻尼特性的复参数修正方法研究

阻尼对于结构动力学响应具有重要的影响,但有限元模型一般很难对阻尼特性进行精确建模.基于实测频响函数,研究了一种有限元模型阻尼特性的复参数修正方法.以待修正区域各单元质量、刚度矩阵的比例修正系数为复修正参数,建立了单元矩阵比例修正的灵敏度方程直接算法,并对比分析了复修正参数与不同阻尼特性之间的数学关系.以六自由度集中参数模型和25杆平面桁架模型为例,验证了复参数修正方法在阻尼特性修正中的有效性.     

结构动力学有限元模型修正的目标函数及算法

结构动力学有限元模型修正是结构动力学领域的一个热点问题,回顾了结构动力学有限元模型修正研究的发展历史和现状,简要评述了结构动力学有限元模型修正所使用的设计变量,着重阐述了各种结构动力学有限元模型修正方法中所使用的目标函数及修正算法,讨论了工程结构动力学有限元模型修正的一些策略,最后对结构动力学有限元模型修正技术的发展进行了总结和展望.     

基于不完备频响函数的结构损伤统计识别研究

基于不完备频响函数数据,结合概率统计方法,提出了一种能同时考虑模型参数不确定性和测试噪声影响的结构损伤统计识别方法.首先,基于频响函数某一行向量在不同频率下的幅值数据,利用矩阵拉直运算建立了关于损伤系数的确定性识别方程.其次,假设模型参数误差和测试噪声为零均值的高斯随机变量,根据摄动理论,推导了损伤后结构刚度参数的前二...     

基于摄动法的不确定性有限元模型修正方法研究

开展了考虑不确定性的有限元模型修正方法的研究。基于摄动法推导了待修正参数均值和协方差矩阵的迭代格式,其中协方差的迭代格式包括是否考虑试验数据与修正参数之间相关性的两种形式。在理论研究基础上开展数值仿真研究,实现了不确定性有限元模型修正的摄动法,并研究了试验数据样本数量对修正误差的影响。仿真结果表明,该方法适用于解决系统参数与试验数据存在不确定性的模型修正问题,试验样本数量对待修正参数标准差的修正精度影响较大;忽略试验模态参数与待修正参数不确定性之间的相关性,能够避免计算二阶灵敏度矩阵,在保证修正结果准确性的前提下减少计算量。     

Recent progress on finite element model updating: From linearity to nonlinearity

有限元分析在实际工程中得到了广泛应用.然而有限元模型由于受到网格划分、边界条件和材料物理参数不确定性等的影响,与真实结构有差异. 因此须通过试验数据加以修正,使其尽可能接近实际结构,以保证之后的结构动力模拟分析和监测等具有实际意义. 经过多年发展,有限元模型修正技术已经能够成功应用于一些实际工程,但现代工程技术的进步对有限元模型修正提出了更高要求,修正后的有限元模型不仅要有较高的精确度,还需要为后续应用给出具有指导意义的置信度.而现有的有限元模型修正、确认方法多基于结构线性的假设,而未能考虑实际结构中广泛存在的非线性.因此本文以土木工程结构模型修正的一些研究成果为例,通过对传统有限元模型修正的发展历程进行全面回顾;总结评述传统有限元修正技术的主要方法,以及包括有限元模型确认在内的最新研究进展;重点探讨有限元模型修正技术向非线性发展的技术路线和目前主要研究成果,展望其未来发展方向, 并提出值得研究的问题.     

基于Kriging代理模型的迭代更新高效反演方法

为提高混凝土坝等大体积结构参数反演效率和精度,减少由于应用有限元进行大量正分析而产生的计算机时,建立了一种结合Kriging代理模型和粒子群优化(PSO)算法的迭代更新反演方法。通过拉丁超立方抽样(LHS)方法确定初始样本点的空间分布,并使用有限元正分析获取对应的响应值,构建粗糙的初始代理模型,结合具有全局寻优能力的PSO算法,反演大体积结构的分区弹性模量,随之再代入有限元模型中,计算获取新的位移响应,并将其作为新样本加入到样本集中,通过迭代更新获得局部更高精度的代理模型。工程实际算例表明,该方法对混凝土坝等大体积结构参数反演精度较高和适用性好,且能大幅减少传统有限元模型反演方法所需消耗的正分析机时,提高反演效率。     

Research on the iterative method for model updating based on the frequency response function

Model reduction technique is usually employed in model updating process.In this paper,a new model updating method named as cross-model cross-frequency response function(CMCF) method is proposed and a new iterative method associating the model updating method with the model reduction technique is investigated.The new model updating method utilizes the frequency response function to avoid the modal analysis process and it does not need to pair or scale the measured and the analytical frequency response function,which could greatly increase the number of the equations and the updating parameters.Based on the traditional iterative method,a correction term related to the errors resulting from the replacement of the reduction matrix of the experimental model with that of the finite element model is added in the new iterative method.Comparisons between the traditional iterative method and the proposed iterative method are shown by model updating examples of solar panels,and both of these two iterative methods combine the CMCF method and the succession-level approximate reduction technique.Results show the effectiveness of the CMCF method and the proposed iterative method.     

Residual stiffness assessment of structurally failed reinforced concrete structure by dynamic testing and finite element model updating

When an under-reinforced concrete beam structure has been loaded to the point where reinforcing steel on the tension side has yielded, it is deemed to have structurally failed and the full load capacity and stiffness can no longer be developed. When unloaded from the point of failure, the residual stiffness of the structure is difficult to estimate. There is a need to establish the serviceability of the structure and ultimately establish the levels of further loading that can be sustained before total collapse. In this paper we present a method for assessing residual stiffness of such a “failed” reinforced concrete structure, through the application of dynamic testing and finite element (FE) model updating. In an experimental study, failed zones in a beam structure were simulated in a FE model. Through a procedure of sensitivity-based updating using the measured modal properties, the stiffness distribution along the failed beam structure was identified.     

Two-phase flow analysis based on a phase-field model using orthogonal basis bubble function finite element method

In this study, a finite element method based on a phase-field model for gas–liquid two-phase flow is proposed. MINI element based on a bubble function element stabilisation method is employed for the incompressible Navier–Stokes equations. The Cahn–Hilliard equation is employed to estimate the interface of gas and liquid. The orthogonal basis bubble function element is used to solve the Cahn–Hilliard equation. In particular, a detailed explanation for solving the Cahn–Hilliard equation based on a finite element method is given.     

The finite element method based on interpolating with wavelet basis function

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｐｒｏｂｌｅｍｓｗｉｔｈｌａｒｇｅｇｒａｄｉｅｎｔａｒｅｃｏｍｍｏｎｉｎｐｒａｃｔｉｃａｌｅｎｇｉｎｅｅｒｉｎｇｆｉｅｌｄｓ,ｅ．ｇ．ｉｎｍａｔｅｒｉａｌｌｏｃａｌｉｚａｔｉｏｎ,ｗｉｔｈｉｎｔｈｅｌｏｃａｌｉｚａｔｉｏ...