时域内多源动态载荷的一种计算反求技术

载荷在时域内可用一系列脉冲或阶跃的核函数来表示,系统的响应是载荷与核函数相对应响 应的卷积分. 在线性时不变的假设下,对系统动力响应的卷积分进行离散,并在此 基础上分析载荷识别反问题的不适定性. 针对测量的响应数据中存在噪声时载荷识 别的困难,探讨了稳定近似识别载荷的一些方法,包括零相位滤波技术、几种正则 化方法和优化策略. 数值仿真算例表明,所述的载荷识别方法能够在响应数据含 有噪声的情况下,有效稳定地实现多源动态载荷的重构.     

基于线性多点逼近的载荷识别方法

利用结构响应信息进行载荷识别是结构工程重要研究内容之一.本文提出了基于时空多点线性逼近的动载荷识别方法,利用特定的载荷基函数在空问和时间上多点线性组合逼近结构的未知荷载输入,实现结构载荷的参数化表达,使得载荷识别反问题归结为一个多点权重系数求解问题.所建立的识别方法不需要对结构模型预先进行模态分析,采用一个常系数矩阵对...     

一种基于遗传算法的时域内稳态动载荷识别方法

为避免传统载荷识别过程中产生的矩阵求逆病态和对初值敏感及累计误差等问题,将遗传算法应用到载荷识别过程中,将此动力学的反问题转化为结构动力学正计算,并且利用遗传算法进行寻优得到最优参数,从而得到待识别载荷时间历程的估计。通过仿真计算,所识别载荷计算振动响应与测量响应的误差为7%,结果表明本文提出的识别方法是有效的。     

基于有限元法-正则化的弹性模量反求算法研究

生物组织的病变往往导致其物理特性的变化,研究生物组织的物理特性变化对诊断生物组织病变与否能够提供参考.本文提出了根据已知生物组织部分边界位移信息、外部载荷、边界条件和有限元模型,利用正则化方法反求弹性模量的方法.首先将求解弹性模量问题描述成一个力学反问题,建立逆系统辨识的数学模型,并应用正则化原理建立生物组织弹性模量反求算法.最后利用仿真算例验证了算法的有效性,并结合图像识别技术试验获取对象部分边界位移信息,反求其弹性模量,证明了本文反求方法的可行性.     

基于正则化等效源模型的航空重力梯度测量信息降噪方法

航空重力梯度测量能获取重力梯度的多个分量,在保证各梯度分量内部固有约束条件下降低测量噪声,是一个巨大的挑战.为提高航空重力梯度测量精度,提出一种基于正则化等效源模型的测量降噪方法.根据矩形棱柱体与重力梯度之间的关系构建了等效源线性方程组,针对该方程组可能存在病态性的问题,引入截断奇异值TSVD和Tikhonov正则化的...     

微分方程反问题的梯度正则化法

本文在微分方程反问题研究当前发展状况的基础上,提出了一种求解微分方程反问题的比较统一的方法——梯度正则化法。该方法将反问题的求解分为两个主要部分,一是展开以求补充条件对未知参数的梯度矩阵,二是正则化以解不适定的线性方程组。梯度正则化法从反问题的共性入手,未附加任何特殊的约束,所以可以适用各种类型的反问题,且在求解时不受空间维数的限制。     

自由-自由运行体系动态载荷反演问题研究

对自由-自由运行体系动态载荷识别问题进行深入研究,将自由-自由运行体系振动积分方程时域离散化,建立了体系动态载荷反演的力学模型。在动态载荷反演模型求解过程中,由于结构矩阵的病态特性以及测量噪声的影响,常规最小二乘法往往失效;通过对载荷反演模型进行奇异值分解,指出该病态问题的本质,并提出相应的正则化求解方法。数值仿真表明了本文方法能够得到满足工程要求的稳定近似解。本文研究对于实际飞行器运行过程所受动载荷评估具有一定的意义。     

一种时域载荷识别数值算法及其不适定性研究

某结构在试验中,经历了持续时间很短的出水过程,试验中获得了该状态下关键位置响应的时间历程,需进行结构试验环境适应性分析。本文将杜哈梅卷积积分离散化,推导了时域内载荷与响应之间的传递函数矩阵;基于该传递矩阵的逆运算发展了一种时域内载荷识别方法,该方法在数值仿真过程中往往出现累计误差过大导致不收敛的现象,研究分析发现这属于反问题中典型的不适定性问题,并采用最速下降法有效克服了该不适定性问题;同时还考察了该算法在不同信噪比下的识别精度。研究结果表明,正弦、随机、半正弦冲击等多种载荷均获得了较准确的识别结果。本文提出的方法便于复杂工程结构有限元模型的应用。     

二阶非定常多宗量热传导反问题的正则解

引入Bregman距离函数及其加权函数作为正则项，应用Tikhonov正则 化方法，对二阶非定常多宗量热传导反问题进行求解. 利用测量信息和计算信息构造最小二 乘函数，将多宗量反演识别问题转化为一个优化问题. 空间上采用8节点等参元进行离散， 时域上采用时域精细算法进行离散，建立了二阶非定常多宗量热传导问题的有限元正/反演数 值模型. 该模型不仅考虑了非均质和参数分布的影响，而且也便于正反演问题的敏度分析， 可对导热系数和边界条件等宗量进行有效的单一和组合识别. 给出了相关的数值验证，对信 息测量误差以及不同正则项的计算效率作了探讨. 数值结果表明，该方法能够对二阶非定常 多宗量热传导反问题进行有效的求解，并具有较高的计算精度.     

一种直接求解的动态载荷识别方法

基于直接求解系统方程的位移响应,提出了一种时域内动态载荷的识别方法.该方法从中心差分格式的有限元离散形式出发,进行分块矩阵变换,将载荷识别问题转化为一个虚构系统方程的响应求解问题.通过数值算例,用两种不同形式的载荷检验了该方法的识别能力,分析了有无零相位滤波器对识别结果的影响.结果表明:所提出的载荷识别方法在位移响应数据含有噪声的情况下,能有效稳定地实现载荷识别;且计算简单快捷,抗噪能力强,不需进行模态分解,不存在模态截取等问题.     

A new method for detecting non-linear damping and restoring forces in non-linear oscillation systems from transient data

The purpose of this study is to recover the functional form of both non-linear damping and non-linear restoring forces in the non-linear oscillatory motions of an autonomous system. Using two sets of measured motion response data of the system, an inverse problem is formulated for recovering (or identification): the differential equation of motion is transformed into an equivalent integral equation of motion. The identification, which is non-linear, is shown to be one-to-one. However, the inverse problem formulated herein is concerned with the Volterra-type of non-linear integral equation of the first kind. This leads to numerical instability: solutions of the inverse problem lack stability properties. In order to overcome the difficulty, a regularization method is applied to the identification process. In addition, an L-curve criterion, combined with regularization, is introduced to find an optimal choice for the regularization parameter (i.e., the number of iterations), in the presence of noisy data. The workability of the identification is investigated for simultaneously recovering the functional form of the non-linear damping and the non-linear restoring forces through a numerical experiment.     

Force reconstruction based on the union of singular spectral analysis and regularization method

Abstract

To overcome the problem of some ill-posed inverse problem of force reconstruction, which is caused by the noise in the measured responses and small singular values of the structure, a technology of force reconstruction based on a hybrid method of singular spectral analysis (SSA) and the Landweber regularization method is proposed in this study for the first time. The SSA is used to filter the structural response before using Landweber regularization. A new choice method of phase space reconstruction dimension is theoretically proposed, and the minimum embedding dimension is determined by the concept of optimizing difference spectrum theory. The feasibility of this method was demonstrated through three kinds of force reconstructions. The numerical simulation results and an acoustic vibration experiment demonstrated that the proposed method is more effective than the traditional method.

Communicated by Wei-Chau Xie.     

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。     

On the regularization method of the first kind of Fredholm integral equation with a complex kernel and its application

The regularized integrodifferential equation for the first kind of Fredholm integral equation with a complex kernel is derived by generalizing the Tikhonov regularization method and the convergence of approximate regularized solutions is discussed. As an application of the method, an inverse problem in the two-demensional wave-making problem of a flat plate is solved numerically, and a practical approach of choosing optimal regularization parameter is given. Project supported by the National Natural Science Foundation, of China     

Reconstruction of high order derivatives by new mollification methods

In this paper, the problem of reconstructing numerical derivatives from noisy data is considered. A new framework of mollification methods based on the L generalized solution regularization methods is proposed. A specific algorithm for the first three derivatives is presented in the paper, in which a modification of TSVD, termed cTSVD is chosen as the regularization technique. Numerical examples given in the paper verify the theoretical results and show efficiency of the new method.     

High resolution COD image analysis for health monitoring of reinforced concrete structures through inverse analysis

Direct and inverse problems of a fracture mechanics based RC beam model are solved. Solution of the direct problem that maps crack bridging stresses into crack opening displacements (COD) is straightforward, but the inverse problem is ill-posed, and better solved by the theory of inverse problems. This paper exploits the Tikhonov regularization method to solve the inverse problem, and estimates the force and location of rebar in buried concrete from CODs. Bending tests are carried out on model RC beams in the laboratory to demonstrate the applicability of the method. During the tests, a microscopic camera snaps high resolution digital pictures of cracked concrete surface. The images are analyzed by a software to measure surface CODs that are input into the inverse problem. The practical CODs inevitably include noise due to experimental error, which makes the inverse problem ill-posed, and necessitates regularization. In the current inverse analysis by the Tikhonov regularization method, bridging stress profiles, i.e. variation of the crack bridging stress along the crack length, has been figured out. Results are compared with those from other theoretical methods of analysis as well as with the readings from strain gauges. The method is a suitable non-destructive means for existing structures in cases where the section information is inadequate, or damages/repairs have altered the designed cross-section.     

The MFS for the Cauchy problem in two-dimensional steady-state linear thermoelasticity

We study the reconstruction of the missing thermal and mechanical fields on an inaccessible part of the boundary for two-dimensional linear isotropic thermoelastic materials from over-prescribed noisy (Cauchy) data on the remaining accessible boundary. This problem is solved with the method of fundamental solutions (MFS) together with the method of particular solutions (MPS) via the MFS-based particular solution for two-dimensional problems in uncoupled thermoelasticity developed in Marin and Karageorghis, 2012a, Marin and Karageorghis, 2013. The stabilisation/regularization of this inverse problem is achieved by using the Tikhonov regularization method (Tikhonov and Arsenin, 1986), whilst the optimal value of the regularization parameter is selected by employing Hansen’s L-curve method (Hansen, 1998).     

半空间SH波方程非线性井间双参数反演

以半空间的ＳＨ波方程出发,采用Ｂｏｒｎ迭代法求解半空间弹性介质中密度和剪切模量分布的非线性反演问题。首先,采用矩量法和正则化方法,给出井间反演积分方程的离散形式,然后应用Ｂｒｏｎ迭代法求解非线性反演问题。     

二维介质参数的大扰动反演方法

对非均匀介质参数反演问题进行了研究，并提出了用于反演二维介质参数的广义射线近似方法．利用参考场量和扰动变量对声波方程中的介质参数进行处理，并利用Ｇｒｅｅｎ函数理论得到扰动参数比的积分方程．基于非均匀介质中波函数的局部理论和射线理论，引入了全波场的广义射线近似形式，通过定义介质参数函数，把反演目标归结为其第一类Ｆｒｅｄｈｏｌｍ积分方程．利用积分变换方法得到二维介质的介质参数函数，从而得到介质参数，在Ｂｏｒｎ近似方法中，反演的介质参数扰动不能超过２０％，但是在本文中介绍的方法能够有效地反演其扰动比不超过５０％的变化情况     

A smoothed inverse eigenstrain method for reconstruction of the regularized residual fields

A smoothed inverse eigenstrain method is developed for reconstruction of residual field from limited strain measurements. A framework for appropriate choice of shape functions based on the prior knowledge of expected residual distribution is presented which results in stabilized numerical behavior. The analytical method is successfully applied to three case studies where residual stresses are introduced by inelastic beam bending, laser-forming and shot peening. The well-rehearsed advantage of the proposed eigenstrain-based formulation is that it not only minimizes the deviation of measurements from its approximations but also will result in an inverse solution satisfying a full range of continuum mechanics requirements. The smoothed inverse eigenstrain approach allows suppressing fluctuations that are contrary to the physics of the problem. Furthermore, a comprehensive discussion is performed on regularity of the asymptotic solution in the Tikhonov scheme and the regularization parameter is then exactly determined utilizing Morozov discrepancy principle. Gradient iterative regularization method is also examined and shown to have an excellent convergence to the Tikhonov–Morozov regularization results.