A new adaptive approach of the Metropolis-Hastings algorithm applied to structural damage identification using time domain data

In the present work, the formulation and solution of the inverse problem of structural damage identification is presented based on the Bayesian inference, a powerful approach that has been widely used for the formulation of inverse problems in a statistical framework. The structural damage is continuously described by a cohesion field, which is spatially discretized by the finite element method, and the solution of the inverse problem of damage identification, from the Bayesian point of view, is the posterior probability densities of the nodal cohesion parameters. In this approach, prior information about the parameters of interest and the quantification of the uncertainties related to the magnitudes measured can be used to estimate the sought parameters. Markov Chain Monte Carlo (MCMC) method, implemented via the Metropolis-Hastings (MH) algorithm, is commonly used to sample such densities. However, the conventional MH algorithm may present some difficulties, for instance, in high dimensional problems or when the parameters of interest are highly correlated or the posterior probability density is very peaked. In order to overcome these difficulties, a new adaptive MH algorithm (P-AMH) is proposed in the present work. Numerical results related to an inverse problem of damage identification in a simply supported Euler-Bernoulli beam are presented. Synthetic experimental time domain data, obtained with different damage scenarios, and noise levels, were addressed with the aim at assessing the proposed damage identification approach. An adaptive MH algorithm (H-AMH) and the conventional MH algorithm, already consolidated in the literature, were also considered for comparison purposes. The numerical results show that both adaptive algorithms outperformed the conventional MH. Besides, the P-AMH provided Markov chains with faster convergence and better mixing than the ones provided by the H-AMH.     

环境激励下基于柔度矩阵分解的损伤诊断

为了解决环境激励作用下结构自由度不完备对损伤诊断的影响,提出了一种基于自由度缩聚的比例柔度矩阵分解损伤诊断法.利用附加质量法求解出环境激励作用下振型关于质量归一化因子.进而根据质量归一化因子和比例柔度矩阵系数之间的关系,构建出其比例柔度矩阵,再通过使用QR矩阵分解法对构建出的比例柔度矩阵进行分解.以分解后得到的三角矩阵(R矩阵)作为研究对象,将〖WTHX〗R〖WTBX〗矩阵经过相应的数学算法处理得到最终损伤定位指标.算例研究表明:在环境激励作用下考虑自由度缩聚的影响,无论对于单损伤还是多损伤,所提损伤定位指标均表现出较高的准确性,且具有一定的鲁棒性.该文基于矩阵分解法推导出的损伤定位指标可以应用在环境激励作用下的损伤诊断,同时也为自由度不完备结构的损伤诊断提供了新研究思路.     

修正的TV-Stokes模型的新快速算法

基于Tai等人的前期工作,本文研究修正的TV-Stokes图像去噪模型,提出一些新的求解该两步模型的快速算法.我们利用对偶形式和多重网格方法得到一个求解第1步的快速算法.给出另外一种新的求解光滑的切向量场的保不可压性质的算法.在第2步中,我们提出一类有效的全新算法:首先通过计算Poisson方程得到具有光滑法向量场的函数g,然后利用Jia和Zhao的方法得到恢复的图像.新算法的运算速度非常快,用于图像恢复的CPU时间少于0.1 s.数值结果显示新的快速算法是有效的和稳定的,恢复图像的质量也超过了一般去噪方法.     

Modelización numérica del comportamiento estructural cíclico de barras esbeltas de acero con pandeo restringido

This work presents a numerical model of the cyclic structural behavior of dissipative buckling-restrained braces, commonly used as an alternative to classical concentric braces for seismic protection of building frames and other structures. Such devices are usually composed of a slender steel core embedded in a stockiest casing that is intended to prevent its buckling when it is under compression. The casing is made either of mortar or steel, and a sliding interface is interposed between the core and the casing to prevent excessive shear stress transfer. The behavior of the steel core is described by a damage and plasticity model; the behavior of the mortar casing is described by an isotropic damage model and the sliding behavior of the interface is described by a contact penalty model. These 3 models are implemented in the Abaqus software package following an explicit formulation. In a previous article (published in an earthquake engineering journal) the model was briefly described, its ability to reproduce the cyclical behavior of buckling-restrained braces was preliminarily pointed out and their results were satisfactorily compared with those of experimental tests. The aim of this paper is to describe the model thoroughly and to present new judgments about its usefulness.     

多元Lanchester方程在装备战损量预计中的应用

在战时装备保障中,装备战损量预计是一难点问题.本文从分析影响装备战损的因素出发,采用多元Lanchester方程寻求解决方案.提出了基于指数多元Lanchester方程的装备战损量预计模型和模型中毁伤系数的确定方法.据此,建立了一种新的装备战损量预计方法,并举例验证.该方法较好的满足了未来信息化条件下作战装备保障的需要,具有重要的参考价值.     

高压水射流作用下岩石破碎机理及过程的数值模拟研究

根据连续介质力学和有限元理论,给出了高压水射流破岩系统中流体和岩石的控制方程,并建立了相应的有限元列式．运用连续损伤力学和细观损伤力学理论,建立了适用于水射流破岩全过程分析的岩石损伤模型以及宏细观损伤的耦合模式．数值计算的结果较真实地反映了水射流破岩过程中,岩石的动态响应以及水射流动力学特性的演化过程,普通连续水射流破碎岩石主体所用的时间为毫秒量级,破岩的主要形式是卸载及射流冲击所产生的拉伸破坏,并呈“阶跃式”发展．数值计算与相关试验结果基本吻合,表明该分析方法是可行的,可用来指导高压水射流破岩理论的进一步研究及应用．     

IMAGE SEGMENTATION BY PIECEWISE CONSTANT MUMFORD-SHAH MODEL WITHOUT ESTIMATING THE CONSTANTS

In this work, we try to use the so-called Piecewise Constant Level Set Method （PCLSM） for the Mumford-Shah segmentation model. For image segmentation, the Mumford-Shah model needs to find the regions and the constant values inside the regions for the segmen- tation. In order to use PCLSM for this purpose, we need to solve a minimization problem using the level set function and the constant values as minimization variables. In this work, we test on a model such that we only need to minimize with respect to the level set function, i.e., we do not need to minimize with respect to the constant values. Gradient descent method and Newton method are used to solve the Euler-Lagrange equation for the minimization problem. Numerical experiments are given to show the efficiency and advantages of the new model and algorithms.     

药物透皮参数的估计与研究

参考透皮给药的大量实验数据,选择一种非线性回归模型作为相应的数学模型,并通过这一模型中的透皮系数作为被估计的参数完成了这一工作,同时给出了对20种药物的不同渗透系数的估计．最后的讨论中发现,这一工作进一步支持了将皮肤视作单层渗透膜的理论．     

垂直裂缝井试井分析模型和方法

根据压裂井的流动机制,综合考虑井筒储存、裂缝壁面污染和各种边界条件,通过组合线性流模型与有效井径模型,建立了垂直裂缝井试井分析的新模型,提出了确定垂直裂缝井有效井径的方法,给出了有效井径随裂缝长度、裂缝导流能力和裂缝表皮系数的变化关系．该模型形式简明,曲线完整,计算速度快,可满足实时计算和快速响应的试井解释要求．以实测的压力或产量为拟合目标函数,建立了识别地层和水力裂缝参数的最优化模型,提出了综合应用逐步线性最小二乘法和约束变尺度法的垂直裂缝井试井分析自动拟合方法,并通过算例说明了该方法的准确性和可靠性．应用情况表明,该技术可科学合理地分析评估压裂施工质量,指导并改进压裂设计,提高压裂设计水平和施工效果．     

On a New Diffeomorphic Multi-Modality Image Registration Model and Its Convergent Gauss-Newton Solver

In this work, we propose a new variational model for multi-modal image registration and present an efficient numerical implementation. The model minimizes a new functional based on using reformulated normalized gradients of the images as the fidelity term and higher-order derivatives as the regularizer. A key feature of the model is its ability of guaranteeing a diffeomorphic transformation which is achieved by a control term motivated by the quasi-conformal map and Beltrami coefficient. The existence of the solution of this model is established. To solve the model numerically, we design a Gauss-Newton method to solve the resulting discrete optimization problem and prove its convergence; a multilevel technique is employed to speed up the initialization and avoid likely local minima of the underlying functional. Finally, numerical experiments demonstrate that this new model can deliver good performances for multi-modal image registration and simultaneously generate an accurate diffeomorphic transformation.     

Structural damage identification based on best achievable flexibility change

A new method based on best achievable flexibility change is presented in this paper to localize and quantify damage in structures. Central to the damage localization approach is the computation of the Euclidean distances between the measured flexibility change and the best achievable flexibility changes. The location of damage can be identified by searching for a value that is considerably smaller than others in these distances. With location determined, a simple extent algorithm is then developed. Three examples are used to demonstrate the efficiency of the method. Results indicate that the proposed procedure may be useful for structural damage identification.     

Analytical solution of the droplet breakup equation by the Adomian decomposition method

The droplet population balance equations (PBEs), which despite its importance rarely has an analytical solution. However, few cases with assumed functional forms of breakup rate, daughter droplet distribution exist, where most of these solutions are for the batch stirred vessel. A new framework for solving (PBEs) for batch and continuous systems is proposed in this work, which uses the Adomian decomposition method. This technique overcomes the crucial difficulties of numerical discretization and stability that often characterize previous solutions in this area. The technique used in this work has been tested for the droplet breakup equation. The solutions are presented for several cases, for which analytical solutions are available, for batch and continuous systems for droplet breakup in stirred vessels. In all cases, the predicted droplet size distributions converge exactly in a continuous form to that of the analytical solution.     

LMI Approach to Robust Model Predictive Control

This paper introduces a new approach to robust model predictive control (MPC) based on conservative approximations to semi-infinite optimization using linear matrix inequalities (LMIs). The method applies to problems with convex quadratic costs, linear and convex quadratic constraints, and linear predictive models with bounded uncertainty. If the MPC optimization problem is feasible at the initial control step (the first application of the MPC optimization), it is shown that the MPC optimization problems will be feasible at all future time steps and that the controlled system will be closed-loop stable. The method is illustrated with a solenoid control example. The authors thank the anonymous reviewers for suggestions that improved the presentation of this work. The work was supported in part by the EPRI/DoD Complex Interactive Networks/Systems Initiative under Contract EPRI-W08333-05 and by the US Army Research Office Contract DAAD19-01-1-0485.     

弹塑性损伤问题的自洽分析方法

基于不可逆热力学理论和内变量理论，提出平均意义下的Ｃｌａｕｓｉｕｓ＿Ｄｕｈｅｍ不等式方程·与以往弹塑性空洞损伤研究不同，本文旨在利用自洽分析法建立起既考虑空洞形状影响又考虑空洞之间相互作用的场方程，并使之能够处理空洞体积比较大的弹塑性损伤问题·     

The mixed and multi model line balancing problem: a comparison

The balancing problem deals with the assignment of tasks to work stations. We can distinguish two approaches in the literature on the mixed model line balancing problem, that both transform this problem into a single model line balancing problem. These approaches use combined precedence diagrams and adjusted task processing times, respectively.An experiment was carried out to compare several heuristics based on the combined precedence diagram. A new optimisation method has been developed. The results indicate that the position of common tasks in the precedence diagram of the different models has a significant effect on both the CPU time and the unequal distribution of the total work content of single models among work stations. Moreover, good solutions with respect to the number of required stations go together with long CPU times. For several instances, we decreased the CPU times considerably without deteriorating the performance of the methods, by using a reversed combined precedence diagram.     

动力系统实测数据的非线性混沌模型重构

动力系统实测非线性混沌数据的模型重构技术是相空间重构的重要内容。在判定了实测数据的非线性混沌特征,计算了实测数据的分维数,Lyapunov指数,并对其进行了本征值分解和噪声去除及确定其模型阶数以后,提出了一个动力系统实测数据的非线性混沌模型,给出了相应的模型参数辨识方法,并用其确立的混沌模型进行了预测工作,计算结果表明:模型参数辨识方法能迅速地将参数估计值带到多峰目标函数的全局最少值附近,然后再采用优化理论能较准确地求出模型的参数,用得到的混沌模型对系统进行预测工作其预测效果良好,且混沌时序不可能作长期预测。     

单轴作用下岩土材料的双重介质本构模型

基于弹塑性力学和损伤力学理论,将岩土材料视为孔隙-裂隙双重介质,假设孔隙介质不发生损伤,而裂隙介质随应变的增加发生损伤,建立了单轴作用下岩土类材料的双重介质本构模型隐式表达式,并利用Newton迭代法得出了材料的全程应力-应变曲线．分析结果表明,岩土材料中裂隙空间展布的多态性（均匀展布、集中展布和随机展布）是岩土材料本构关系千变万化的根本原因．由于双重介质本构模型将岩土材料的弹性主体（孔隙介质部分）和损伤主体（裂隙介质部分）分化开来,对于研究岩土或含损伤材料的破坏具有实用价值和理论意义．     

A new method for topology design of electromagnetic antennas in hyperthermia therapy

The topological derivative concept has been proved to be useful in many relevant applications such as topology optimization, inverse problems, image processing, multi-scale constitutive modeling, fracture mechanics and damage evolution modeling. In this work, we develop a new optimization method based on the topological derivative concept applied to the cancer treatment by hyperthermia. Hyperthermia therapy is a non-invasive medical treatment in which body tissue is artificially heated through electromagnetic waves, focusing the heat in cancerous cells undergoing apoptosis. The basic idea, therefore, consists in finding a distribution of heat source generated by electromagnetic antenna aiming to increase the temperature in the region occupied by the tumor, while keeping the temperature in the remainder part of the body. Numerical results are presented illustrating possible application of the proposed methodology to treatment of cancer by hyperthermia.     

An analytical dynamic model for single-cracked beams including bending,axial stiffness,rotational inertia,shear deformation and coupling effects

This paper develops an analytical dynamic model for cracked beams including bending, axial stiffness, rotational inertia, shear deformation and the coupling of the last two effects. The damage is modelled using a rotational spring that simulates the crack based on fracture mechanics theory. The developed model is used to predict variations on natural frequencies for several crack sites and damage magnitude along the beam. The importance of this work lies in the development of an analytical model that has no approximation due to discretization of the displacement field. This initial theoretical approach describes the expected behaviour for changes in the natural frequencies for simply-supported and clamped-free beams with the precision that only analytical methods allow. The results provide a useful benchmark to compare with approximate numerical methods that can be used to model and analyse the problem. The model showed similar results for long span beams, but the inclusion of rotational inertia and shear deformation effects rendered improvements in the dynamic behaviour mainly in the case of slender and short span beams when compared with the simplified Euler–Bernoulli model.