粗糙表面之间接触热阻反问题研究

当两个固体表面相互接触时,由于接触面粗糙度的影响,界面间就形成了非一致接触,这种接触导致热流收缩,进而产生接触热阻. 目前的理论研究主要集中在正问题研究,对反问题的研究相对较少. 接触热阻反问题研究是通过研究部分边界温度、热流和部分测量点的温度来反演得到界面上的接触热阻. 反问题研究在很多工程领域都有应用,如航空航天、机械制造、微电子等,是工程中确定接触热阻一种快速有效的方法. 本文采用边界元法和共轭梯度法研究了二维空间随坐标变化的接触热阻反问题. 为了验证方法的准确性和可行性,假定在已知部分测量点温度和真实接触热阻的情况下,反演计算得到界面的温度和热流,进而得到接触热阻,并与真实接触热阻进行比较. 结果表明采用边界元法和共轭梯度法在无测量误差的情况下,可以准确反演获得界面的真实接触热阻. 若存在测量误差,反演计算结果对测量误差极其敏感,反演结果误差会由于测量误差的引入而被放大. 为处理这种不适定性, 采用最小二乘法对反演计算结果进行校正,结果表明采用最小二乘法能够避开反问题中一些偏离实际值较大的测量点,显著提高反演计算结果的准确性.      

Numerical nonlinear inverse problem of determining wall heat flux

The inverse problem of determining time-variable surface heat flux in a plane wall, with constant or temperature dependent thermal properties, is numerically studied. Different kinds of incident heat flux, including rectangular waveform, are assumed. The solution is numerically solved as a function estimation problem, so that no a priori information for the functional waveforms of the unknown heat flux is needed. In all cases, a solution in the form of a piece-wise function is used to approach the incident flux. Transient temperature measurements at the boundary, from the solution of the direct problem, served as the simulated experimental data needed as input for the inverse analysis. Both direct and inverse heat conduction problems are solved using the network simulation method. The solution is obtained step-by-step by minimising the classical functional that compares the above input data with those obtained from the solution of the inverse problem. A straight line of variable slope and length is used for each one of the stretches of the desired solution. The influence of random error, number of functional terms and the effect of sensor location are studied. In all cases, the results closely agree with the solution.     

Criteria of optimality for sensors' location based on adjoint transformation of observation data interpolation error

Criteria of optimality for sensors' location are addressed using an interpolation error transformed by especial adjoint problems. The considered criteria correspond to the analysis error in certain Hessian‐based metrics and to the error of some forecast aspect. Both criteria are obtained using adjoint problems that provide computation without the direct use of the Hessian. For a linear inverse heat conduction problem, these criteria are compared and demonstrated promising results when compared with a criterion based on the norm of the interpolation error of observation data. Approaches to sensor set modification using either redistribution of sensors' or refinement of the sensors grid (insertion of additional sensors) are also compared. Copyright © 2009 John Wiley & Sons, Ltd.     

Non-linear least-squares solution to the moiré hole method problem in orthotropic materials. Part II: Material elastic constants

This is the second of two papers dealing with inverse problems arising from the hole method problem in orthotropic materials. Two inverse problems are of interest in this second paper: one relates to the use of the hole method as a means of separation of stresses; the other deals with using it for orthotropic material elastic constant identification. The problem of separation of stresses is posed and briefly discussed, but the orthotropic material identification problem is addressed, fully showing how to obtain the five orthotropic elastic constants, four of which are independent, when knowing the applied stresses and examining the isothetics or moiré fringes around a through hole in a biaxially loaded thin wood plate.     

Modélisation thermique d'une interface solide–fluide par fonction de transfert

A thermal modelling of a solid–fluid interface is proposed. The developed procedure is based on a numerical identification, in a dimensionless space, of a unique transfer function, associated with a recurrent relation. The method is usable for direct or inverse problems and well fit for control–command applications. An inverse numerical simulation example is presented for a random thermal excitation.     

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.     

An adaptive procedure for defect identification problems in elasticity

In the context of inverse problems in mechanics, it is well known that the most typical situation is that neither the interior nor all the boundary is available to obtain data to detect the presence of inclusions or defects. We propose here an adaptive method that uses loads and measures of displacements only on part of the surface of the body, to detect defects in the interior of an elastic body. The method is based on Small Amplitude Homogenization, that is, we work under the assumption that the contrast on the values of the Lamé elastic coefficients between the defect and the matrix is not very large. The idea is that given the data for one loading state and one location of the displacement sensors, we use an optimization method to obtain a guess for the location of the inclusion and then, using this guess, we adapt the position of the sensors and the loading zone, hoping to refine the current guess.Numerical results show that the method is quite efficient in some cases, using in those cases no more than three loading positions and three different positions of the sensors.     

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

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

基于虚拟变形法的车-桥耦合系统移动质量识别

利用双自由度质量-弹簧阻尼模型模拟移动车辆, 并基于虚拟变形(VDM)方法的结构快速重分 析思想, 提出一种车-桥耦合系统的移动质量快速识别的有效方法. 该方法以双自由度车体模 型的质量为变量, 通过最小化桥体结构实测响应和计算响应的平方距离来识别移动质量 (载荷), 避免了识别载荷时常遇到的病态问题, 对噪声鲁棒性强, 且需要传感器信息少. 每步优化 中, 利用在VDM方法基础上提出的移动动态影响矩阵概念, 无需时时重构车-桥耦合系统的时 变系统参数矩阵, 显著提高了计算效率. 利用数值框架梁模型, 通过比较不同车辆简化模型 对移动体质量及等效移动载荷的识别效果, 验证了该方法的可行性和有效性, 即使在5% 的噪声影响下, 利用一个传感器可以准确地识别多个移动体的质量.     

Use of experiment and an inverse method to study interface heat transfer during solidification in the investment casting process

A technique to determine the thermal boundary conditions existing during the solidification of metallic alloys in the investment casting process is presented. Quantitative information about these conditions is needed so that numerical models of heat transfer in this process produce accurate results. In particular, the variation of the boundary conditions both spatially and temporally must be known. The method used involves the application of a new inverse heat conduction method to thermal data recorded during laboratory experiments of aluminium alloy solidification in investment casting shell moulds. The resultant heat transfer coefficient for the alloy/mould interface is calculated. An experimental programme to determine requisite mould thermal properties was also undertaken. It was observed that there is significant variation of the alloy/mould heat transfer coefficient during solidification. It is found to be highly dependent on the alloy type and on the vertical position below the initial free surface of the liquid metal. The aluminium casting alloys used in this study were 413, A356, 319 (Aluminum Association designations), and commercially pure aluminium. These alloys have significantly different freezing ranges. In particular, it was found that alloys with a high freezing range solidify with rates of heat transfer to the mould which are very sensitive to metallostatic head.     

The inverse identification problem and its technical application

This paper presents an overview of a loading force identification technique. Load identification methods are based on the solution of the inverse identification problem. Many different approaches for linear systems have been developed in this area. For both linear and nonlinear systems, methods based on the minimization of assumed objective functions are formulated. The least square error between the simulated and measured system responses is mainly used as the objective function. The dynamic programming optimization method formulated by Bellman is commonly used for the minimization of the objective function to estimate the excitation forces. The inverse identification problem in most practical cases is ill-posed because not all the state variables or initial conditions are known. Ill-posed inverse identification problems can be solved using several techniques, the most useful of which are: the generalized cross-validation method, the dynamic programming technique and Tikhonov’s method. This article presents the theoretical background and main limits to the application of inverse identification methods. Numerical and experimental tests on a laboratory rig were made to verify the formulated procedures. The method is applied to the identification of wheel–rail contact forces during rail vehicle operation. The method can be applied for indirect measurements of contact forces in railway equipment testing.     

A numerical study of heat transfer during flow of melt in a cylindrical mould

A numerical analysis of transient heat transfer during the flow of a melt in a cylindrical mould is presented. The analysis includes thermal resistance at the melt-mould interface, and axial conduction inside both melt and mould. Energy equations are formulated in a domain that expands continuously due to the advance of the melt inside the empty mould, and solved by the finite difference method using a time-stepping procedure. Calculations are compared to existing analytic results. It is found that axial conduction in the melt can significantly influence the rate of heat loss from the flowing melt, and that analytic approximations, which neglect axial conduction, may give erroneous predictions for the rate of heat loss.     

Simulation of three-dimensional polymer mould filling processes using a pseudo-concentration method

Mould filling processes, in which a material flow front advances through a mould, are typical examples of moving boundary problems. The moving boundary is accompanied by a moving contact line at the mould walls causing, from a macroscopic modelling viewpoint, a stress singularity. In order to be able to simulate such processes, the moving boundary and moving contact line problem must be overcome. A numerical model for both two- and three-dimensional mould filling simulations has been developed. It employs a pseudo-concentration method in order to avoid elaborate three-dimensional remeshing, and has been implemented in a finite element program. The moving contact line problem has been overcome by employing a Robin boundary condition at the mould walls, which can be turned into a Dirichlet (no-slip) or a Neumann (free-slip) boundary condition depending on the local pseudo-concentration. Simulation results for two-dimensional test cases demonstrate the model's ability to deal with flow phenomena such as fountain flow and flow in bifurcations. The method is by no means limited to two-dimensional flows, as is shown by a pilot simulation for a simple three-dimensional mould. The reverse problem of mould filling is the displacement of a viscous fluid in a tube by a less viscous fluid, which has had considerable attention since the 1960's. Simulation results for this problem are in good agreement with results from the literature. © 1998 John Wiley & Sons, Ltd.     

基于混合遗传算法的动力系统阻尼参数识别方法

将动力系统阻尼参数识别反问题转化为非线性优化问题处理，提出了基于遗传算法的动力系统阻尼参数识别方法。为了提高简单遗传算法的计算效率和处理早熟问题，将模拟退火算法与遗传算法相结合，建立了混合遗传算法。数值计算结果表明，本文所建立的方法对于求解参数识别反问题和非线性优化问题是非常有效的，并且具有良好的鲁棒性和全局收敛能力。     

反求压电薄板智能结构荷载的有限元逆逼近方法

基于压电薄板智能结构提出一种有限元逆逼近反求结构荷载参量的方法，逼近目标函数以压电电荷响应为参量，迭代初值采用单位荷载预逼近方法；数值算例表明了该方法的有效性，可用于解决压电智能结构的荷载识别问题。     

渗流问题单参数摄动随机反演方法

在正演随机模拟方法的基础上,结合Taylor展开式和随机变量的摄动方程,讨论了随机参数反演问题,提出了摄动随机反演方法,给出了一阶均值反演准则和二阶均值反演准则,提出了单随机变量的均值和方差的表示方法,计算均值时采用改进的遗传算法,计算方差时采用统计的方法。给出了一个Thies模型反演导水系数的例子,计算表明该方法简单实用,效果良好。     

PROBABILISTIC DAMAGE IDENTIFICATION OF STRUCTURES WITH UNCERTAINTY BASED ON DYNAMIC RESPONSES

The probabilistic damage identification problem with uncertainty in the FE model parameters, external-excitations and measured acceleration responses is studied. The uncertainty in the system is concerned with normally distributed random variables with zero mean value and given covariance. Based on the theoretical model and the measured acceleration responses, the probabilistic structural models in undamaged and damaged states are obtained by two-stage model updating, and then the Probabilities of Damage Existence （PDE） of each element are calculated as the damage criterion. The influences of the location of sensors on the damage identification results are also discussed, where one of the optimal sensor placement techniques, the effective independence method, is used to choose the nodes for measurement. The damage identification results by different numbers of measured nodes and different damage criterions are compared in the numerical example.     

Extremum problems of boundary control for steady equations of thermal convection

An inverse extremum problem of boundary control for steady equations of thermal convection is considered. The cost functional in this problem is chosen to be the root-mean-square deviation of flow velocity or vorticity from the velocity or vorticity field given in a certain part of the flow domain; the control parameter is the heat flux through a part of the boundary. A theorem on sufficient conditions on initial data providing the existence, uniqueness, and stability of the solution is given. A numerical algorithm of solving this problem, based on Newton’s method and on the finite element method of discretization of linear boundary-value problems, is proposed. Results of computational experiments on solving extremum problems, which confirm the efficiency of the method developed, are discussed.     

Identification of electromechanical properties of piezoelectric structures through evolutionary optimisation techniques

In this paper a non-destructive method to predict the whole three-dimensional set of electromechanical properties of active plate structures is presented. Since the elastic properties of the patches, along with their piezoelectric properties, have significant effect on the dynamic response of the global structure, the inverse problem of the identification of those properties is stated as a constrained minimisation problem of an error function expressing the difference between the measured eigenfrequencies and the corresponding numerical values. Thus, the adopted strategy relies on the dynamic response of the structure in terms of undamped natural frequencies and makes use of a global optimisation technique. To this purpose, a 3D finite element model is used to model the piezoelectric sensors, while an improved genetic algorithm is used as the optimisation tool in order to perform the solution search. The results are demonstrated through a simulated test case, that shows the accuracy of the method.     

Nodal predictive error model and Bayesian approach for thermal diffusivity and heat source mapping

This article aims at solving a two-dimensional inverse heat conduction problem in order to retrieve both the thermal diffusivity and heat source field in a thin plate. A spatial random heat pulse is applied to the plate and the thermal response is analysed. The inverse approach is based on the minimisation of a nodal predictive error model, which yields a linear estimation problem. As a result of this approach, the sensitivity matrix is directly filled with experimental data, and thus is partially noisy. Bayesian estimators, such as the Maximum A Posteriori and a Markov Chain Monte Carlo approach (Metropolis–Hastings), are implemented and compared with the Ordinary Least Squares solution. Simulated temperature measurements are used in the inverse analysis. The nodal strategy relies on the availability of temperature measurements with fine spatial resolution and high frequency, typical of nowadays infrared cameras. The effects of both the measurement errors and of the model errors on the inverse problem solution are also analysed.