LQG量测反馈最优控制的精细积分

对于线性二次型高斯(LQG)量测反馈最优控制问题,提出了精细积分解法。根据分离性原理,LQG控制问题可以分成为最优状态反馈控制问题以及最优状态估计问题,即:离线计算的两套黎卡提微分方程的求解以及状态向量的时变微分方程的在线积分解。该算法不仅适用于求解二点边值问题及其相应的黎卡提微分方程,也适用于求解状态估计的时变微分方程。精细积分高精度的特点,对控制和估计都是有利的。数值算例表明了算法的高精度及有效性。     

Abel变换数值反演的积分算子方法

本文研究了Abel变换的数值反演问题.利用Abel变换的理论反演公式与数值求导的积分算子法相结合的方法,对反演公式中奇异积分合理处理,获得Abel变换数值反演的一种算法,并进行了理论分析与数值实验. 结果表明该算法具有计算简单、数值稳定等优点.     

交替方向的扩散方程精细积分并行算法

给出了交替方向的二维扩散方程的精细积分算法，将一个时间步积分分为两个方向，使大规模矩阵的计算转化为一些小矩阵的计算，减小了每一步求解的计算量．对于方形区域的齐次方程，计算结果与全城精细积分完全相同，而计算量和存储量都要小得多．算例表明了算法具有较高的并行计算加速比和计算效率．     

多目标最优化的一种积分型实现算法

在文[1]中给出了求解多目标最优化的一种积分总极值的概念性算法．本文利用数论中的一致分布佳点集列，较为简便的得出了多目标最优化的积分总极值的实现算法和算法终止准则．并经过有关函数数值计算表明该算法是有效的，可用来求解多目标最优化问题的有效解．     

线性时变系统二次最优控制问题的保辛近似求解

状态空间的最优控制体系是保守的,其近似算法应当保辛．提出了基于分段常值精细积分方法的保辛摄动近似方法,在同一框架下求解了线性时变LQ最优控制中的计算问题,即变系数矩阵Riccati方程和状态反馈方程．该算法是保辛的,具有很好的数值稳定性和精度．算例验证了算法的有效性．     

时域自适应精细算法求解二维非线性湿热耦合问题

采用时域自适应精细算法求解二维非线性湿热耦合问题．通过在离散时间段上将变量展开,更精确地描述变量的变化过程,将非线性时空耦合的初边值问题转化为一系列递推的线性边值问题,并利用有限元法进行求解．计算中不需任何附加假设和迭代过程,并提出了自适应计算的收敛判则,以保持时间步长选择不同时的计算精度．数值比较中,分别考虑了材料属性随温度变化、随湿度变化、以及随温度和湿度变化的情形,结果令人满意,体现了这种算法处理复杂非线性问题的能力．     

二维波动方程吸收边界问题的精细积分解法

波动方程一般是描述各种各样的波浪(如声波、光波和水波浪)的微分方程．在科学研究中对各种波浪的研究需转化为求解波动方程定解问题,常用的求解方法有:解析法,级数法,数值方法(如有限元法、有限差分法)等．钟万勰提出了精细积分法(Precise Inte-     

基于比例边界有限元法和灰狼优化算法的裂纹尖端位置识别

基于比例边界有限元法(SBFEM)和灰狼优化(GWO)算法,提出了一种裂纹尖端识别方法。首先,借助SBFEM解决断裂力学问题特有的优势,快速准确地计算出反演所需的测点位移,并验证了正问题求解的正确性。其次,建立与裂纹尖端位置有关的目标函数,将求解裂纹尖端位置转换为求解目标函数最小值的优化问题。最后,采用GWO算法对目标函数进行了优化,进而搜索裂纹尖端的最佳位置。数值算例结果表明:利用SBFEM的高精度、半解析的优点,在反演过程中采用其求解正问题是非常有效的;GWO算法具有良好的全局收敛性,且相比经典的粒子群算法,能够更快速准确地搜索出裂纹尖端的位置;GWO算法具有较好的抗噪性。     

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

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

任意多连通截面扭转应力函数的有限元解法

杆件扭转问题的求解,主要有基于扭转理论翘曲函数的边界元法和有限元法、基于薄壁杆件理论的数值解法和基于扭转理论应力函数的有限元法．根据任意多连通截面直杆扭转问题的应力函数理论,讨论并改进了与微分方程及定解条件等效的泛函,在此基础上推导了求解多连通截面扭转应力函数的有限元列式,将扭转问题的翘曲位移单值条件转化为边界节点上的集中荷载．采用主从节点法满足孔洞边界上应力函数的同值条件,实现了任意多连通复杂截面扭转应力函数的有限元直接求解,通过应力函数积分获得截面的扭转常数．算例验证了方法的可行性和有效性．     

A comparison of some inverse methods for estimating the initial condition of the heat equation

In this work we analyze two explicit methods for the solution of an inverse heat conduction problem and we confront them with the least-squares method, using for the solution of the associated direct problem a classical finite difference method and a method based on an integral formulation. Finally, the Tikhonov regularization connected to the least-squares criterion is examined. We show that the explicit approaches to this inverse heat conduction problem will present disastrous results unless some kind of regularization is used.     

A CQM-based BEM for transient heat conduction problems in homogeneous materials and FGMs

A method for numerical solution of time-domain boundary integral formulations of transient problems governed by the heat equation is presented. The heat conduction problem is analyzed considering homogeneous and non-homogeneous media. In the case of the non-homogeneous media, the conductor material is assumed to be a functionally graded material, i.e., the material properties vary spatially according to known smooth functions. For some specific spatial variations of the material properties, the fundamental solution and the boundary integral equation of the problem are obtained thanks to a change of variables that transforms the original problem to the standard heat conduction problem for homogeneous materials. For the treatment of time-dependent terms, the convolution quadrature method is adopted to approximate numerically the integral equation of the time-domain boundary element method. In the case that the responses are required at a large number of interior points, the convolution performed to calculate them is very time consuming. It is shown that the discrete convolution of the proposed formulation can be computed by means of the fast Fourier transform technique, which considerably reduces the computational complexity. Results for some transient heat conduction examples are presented to validate the numerical techniques studied.     

A numerical technique for solving IHCPs using Tikhonov regularization method

This study is intended to provide a numerical algorithm for solving a one-dimensional inverse heat conduction problem. The given heat conduction equation, the boundary conditions, and the initial condition are presented in a dimensionless form. The numerical approach is developed based on the use of the solution to the auxiliary problem as a basis function. To regularize the resultant ill-conditioned linear system of equations, we apply the Tikhonov regularization method to obtain the stable numerical approximation to the solution.     

Fatigue crack growth simulations of homogeneous and bi-material interfacial cracks using element free Galerkin method

In this paper, quasi-static fatigue crack growth simulations of homogeneous and bi-material interfacial cracks have been performed using element free Galerkin method (EFGM) under mechanical as well as thermo-elastic load. The thermo-elastic fracture problem is decoupled into thermal and elastic problems. The temperature distribution obtained by solving heat conduction equation is used as input in the elastic problem to get the displacement and stress fields. Discontinuities in the temperature and displacement fields are captured by extrinsic partition of unity enrichment technique. The values of stress intensity factors have been extracted from the EFGM solution by domain based interaction integral approach. The standard Paris fatigue crack growth law has been implemented for the life estimation of various model problems. The results obtained by EFGM under mechanical and thermo-elastic loads were compared with those obtained by FEM using remeshing approach.     

热传导(对流-扩散)方程源项识别的粒子群优化算法

提出了利用粒子群优化(PSO)算法反演热传导方程与对流-扩散方程源项的一种新方法,在已有文献方法的基础上,求解出这两类方程正问题的解析解,再把源项识别问题转化为最优化问题,结合粒子群优化算法寻优求解.通过数值模拟与统计检验,结果表明,此方法可快速有效地实现热传导方程与对流-扩散方程源项的识别,并可推广应用到其它数学物理方程的源项或参数的反演识别.     

Use of 3D-transient analytical solution based on Green’s function to reduce computational time in inverse heat conduction problems

Inverse problems can be found in many areas of science and engineering and can be applied in different ways. Two examples can be cited: thermal properties estimation or heat flux function estimation in some engineering thermal process. Different techniques for the solution of inverse heat conduction problem (IHCP) can be found in literature. However, any inverse or optimization technique has a basic and common characteristic: the need to solve the direct problem solution several times. This characteristic is the cause of the great computational time consumed. In heat conduction problem, the time consumed is, usually, due to the use of numerical solutions of multidimensional models with refined mesh. In this case, if analytical solutions are available the computational time can be reduced drastically. This study presents the development and application of a 3D-transient analytical solution based on Green’s function. The inverse problem is due to the thermal properties estimation of conductors. The method is based on experimental determination of thermal conductivity and diffusivity using partially heated surface method without heat flux transducer. Originally developed to use numerical solution, this technique can, using analytical solution, estimate thermal properties faster and with better accuracy.     

Unique solvability of the inverse problem for the heat equation with nonlocal boundary conditions

We consider the inverse problem of source identification for the heat conduction problem. The neoclassical formulation of the direct problem with integral boundary condition is used. Conditions for unique solvability of the inverse problem are obtained. __________ Translated from Prikladnaya Matematika i Informatika, No. 23, pp. 36–50, 2006.     

A weighted algorithm based on the homotopy analysis method: Application to inverse heat conduction problems

In this paper, a weighted algorithm based on the homotopy analysis method (HAM), for solving inverse heat conduction problems, is introduced . Using the HAM, it is possible to find the exact solution or an approximate solution of the problem in the form of a series.     

瞬态热传导的奇异边界法及其MATLAB实现

基于动力学问题时间依赖基本解的奇异边界法是一种无网格边界配点法.该方法引入源点强度因子的概念从而避免了基本解的源点奇异性,具有数学简单、编程容易、精度高等优点.将该方法用于瞬态热传导问题的数值模拟,运用MATLAB实现该问题的数值研究,并创建相应的MATLAB工具箱.针对二维和三维瞬态热传导问题,进行了基于反插值技术和经验公式的奇异边界法MATLAB算例实现.针对支撑圆坯低温瞬态温度场的模拟结果表明,瞬态热传导奇异边界法的MATLAB工具箱具有简单、方便、精确可靠的优点.研究成果有助于发展瞬态热传导的奇异边界法,并为瞬态热传导问题的数值分析和仿真提供了一种简单高效的模拟工具.     

一个抛物型方程侧边值问题的正则逼近解在一类Sobolev空间中的最优误差界

逆热传导问题(IHCP)是严重不适定问题,即问题的解(如果存在)不连续依赖于数据．但目前关于逆热传导问题的已有结果主要是针对标准逆热传导问题．文中给出了出现在实际问题中的一个抛物型方程侧边值问题,即一个含有对流项的非标准型逆热传导问题的正则逼近解一类Sobolev空间中的最优误差界．