A Newton type iterative method for heat-conduction inverse problems

An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.     

Nonlinear implicit iterative method for solving nonlinear ill-posed problems

In the paper, we extend the implicit iterative method for linear ill-posed operator equations to solve nonlinear ill-posed problems. We show that under some conditions the error sequence of solutions of the nonlinear implicit iterative method is monotonically decreasing and, with this monotonicity, prove convergence of the new method for both the exact and perturbed equations.     

A mixed Newton-Tikhonov method for nonlinear ill-posed problems

Newton type methods are one kind of the efficient methods to solve nonlinear ill-posed problems, which have attracted extensive attention. However, computational cost of Newton type methods is high because practical problems are complicated. We propose a mixed Newton-Tikhonov method, i.e., one step Newton-Tikhonov method with several other steps of simplified Newton-Tikhonov method. Convergence and stability of this method are proved under some conditions. Numerical experiments show that the proposed method has obvious advantages over the classical Newton method in terms of computational costs.     

PRECONDITIONED GAUSS-SEIDEL TYPE ITERATIVE METHOD FOR SOLVING LINEAR SYSTEMS

The preconditioned Gauss-Seidel type iterative method for solving linear systems, with the proper choice of the preconditioner, is presented. Convergence of the preconditioned method applied to Z-matrices is discussed. Also the optimal parameter is presented. Numerical results show that the proper choice of the preconditioner can lead to effective by the preconditioned Gauss-Seidel type iterative methods for solving linear systems.     

土中爆炸成腔问题中确定炸点状态的一种计算反求方法

提出了一种通过给定的土中爆炸成腔毁伤效应确定炸点状态的计算反求方法。该方法将确定炸点状态的反问题转化为求解爆炸毁伤效应的计算值与给定值误差函数最小的优化问题。在反求过程中,采用基于误差减小比率技术的多项式近似模型代替土中爆炸数值分析模型,以便提高反求效率。采用Tikhonov正则化方法克服反求过程中出现的病态问题。在此基础上,引入信赖域管理策略判断当前近似模型与实际模型的逼近程度,以确定最优的反求向量。炸点状态反求结果与实验结果的对比分析表明,该方法能够有效且稳定地通过给定的毁伤效应实现炸点状态的反求,这可为炸点状态的设计提供参考。     

Inexact Newton method via Lanczos decomposed technique for solving box-constrained nonlinear systems

This paper proposes an inexact Newton method via the Lanczos decomposed technique for solving the box-constrained nonlinear systems. An iterative direction is obtained by solving an affine scaling quadratic model with the Lanczos decomposed technique. By using the interior backtracking line search technique, an acceptable trial step length is found along this direction. The global convergence and the fast local convergence rate of the proposed algorithm are established under some reasonable conditions. Furthermore, the results of the numerical experiments show the effectiveness of the pro- posed algorithm.     

On monotone iterative method for initial value problems of nonlinear second-order integro-differential equations in Banach spaces

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｐａｐｅｒ[1],Ｄ.Ｇｕｏｅｓｔａｂｌｉｓｈｅｄｔｈｅｅｘｉｓｔｅｎｃｅｏｆｅｘｔｒｅｍｅｓｏｌｕｔｉｏｎｓｏｆｉｎｉｔｉａｌｖａｌｕｅｐｒｏｂｌｅｍｓｆｏｒｆｉｒｓｔ＿ｏｒｄｅｒｉｎｔｅｇｒｏ＿ｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓｏｆＶｏｌｔｅｒｒａｔｙｐｅｉｎＢａｎａｃｈｓｐａｃｅｓ.Ｎｏｗ,ｗｅｃｏｎｓｉｄｅｒｔｈｅＩＶＰｆｏｒｓｅｃｏｎｄ＿ｏｒｄｅｒｉｎｔｅｇｒｏ＿ｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓｉｎａｒｅａｌＢａｎａｃｈｓｐａｃｅＥ:ｕ″=Ｆ(ｔ,ｕ,…     

AN ITERATIVE METHOD FOR THE DISCRETE PROBLEMS OF A CLASS OF ELLIPTICAL VARIATIONAL INEQUALITIES

ＡＮＩＴＥＲＡＴＩＶＥＭＥＴＨＯＤＦＯＲＴＨＥＤＩＳＣＲＥＴＥＰＲＯＢＬＥＭＳＯＦＡＣＬＡＳＳＯＦＥＬＬＩＰＴＩＣＡＬＶＡＲＩＡＴＩＯＮＡＬＩＮＥＱＵＡＬＩＴＩＥＳＺｈｅｎｇＴｉｅ－ｓｈｅｎｓ（郑铁生）ＬｉＬｉ（李立）ＸｕＱｉｎｇ－ｙｕ（许庆余）（Ｄ...     

跨音速翼型反设计的一种大范围收敛方法

求解跨音速翼型的反设计问题时,传统的梯度型方法一般均为局部收敛. 为增大求解的收敛范围,依据同伦方法的思想,通过构造不动点同伦,将原问题的求解 转化为其同伦函数的求解,并依据拟Sigmoid函数调整同伦参数以提高计算效率,进而构造 出一种具有较高计算效率的大范围收敛反设计方法. 数值算例以RAE2822翼型的表面压力分 布为拟合目标,分别采用B样条方法, PARSEC方法及正交形函数方法等3种不同的 参数化方法,并分别以NACA0012, OAF139及VR15翼型为初始翼型进行迭代计 算. 计算结果证明,该方法适用于多种参数化方法,且具有较好的计算效率,从多 个不同的初始翼型出发,经较少次数迭代后, 均能与目标翼型很好地拟合,是一种高效的大范围收敛方法.     

基于物理信息的深度学习求解矩形薄板力学正反问题

建立基于物理信息的神经网络框架,利用深度学习求解矩形薄板力学正反问题.力学正问题为已知矩形薄板的基本参数、边界条件和受力情况,求薄板各点挠度;反问题为已知薄板部分点的挠度、基本参数和受力情况等,识别边界条件.基于物理信息的神经网络模型中,损失函数除基于数据驱动模型的挠度数据拟合部分以外,还引入薄板弯曲基本方程和应力应变...     

A positive interior-point algorithm for nonlinear complementarity problems

1　ＩｎｔｒｏｄｕｃｔｉｏｎａｎｄＡｌｇｏｒｉｔｈｍＧｉｖｅｎａｆｕｎｃｔｉｏｎｆ:Ｒｎ→Ｒｎ,ｗｈｉｃｈｉｓａｓｓｕｍｅｄｔｏｂｅｃｏｎｔｉｎｕｏｕｓｌｙｄｉｆｆｅｒｅｎｔｉａｂｌｅｉｎａｎｏｐｅｎｓｅｔｃｏｎｔａｉｎｉｎｇＲｎ+,ｔｈｅｎｏｎｌｉｎｅａｒｃｏｍｐｌｅｍｅｎｔａｒｉｔｙｐｒｏｂｌｅｍ ,ｄｅｎｏｔｅｄＮＣＰ (ｆ) ,ｉｓｔｏｆｉｎｄａｖｅｃｔｏｒｘ∈Ｒｎｓｕｃｈｔｈａｔｘ≥ 0 ,　ｆ(ｘ) ≥ 0 ,　ｘＴｆ(ｘ) =0 ,Ｃｌｅａｒｌｙ ,ｔｈｅＮＣＰ (ｆ)ｉｓｅｑｕｉｖａｌｅｎｔｔｏｔｈｅｆｏｌｌｏｗ…     

Calculation of nonlinear aerodynamic characteristics of a wing using a 3‐D panel method

The nonlinear aerodynamic characteristic of a wing is investigated using the frequency‐domain panel method. To calculate the nonlinear aerodynamic characteristics of a three‐dimensional wing, the iterative decambering approach is introduced into the frequency‐domain panel method. The decambering approach uses the known nonlinear aerodynamic characteristic of airfoil and calculates two‐variable decambering function to take into consideration the boundary‐layer separation effects for the each section of the wing. The multidimensional Newton iteration is used to account for the coupling between the different sections of wing. The nonlinear aerodynamic analyses for a rectangular wing, a tapered wing, and a wing with the control surface are performed. Present results are given with experiments and other numerical results. Computed results are in good agreement with other data. This method can be used for any wing having different nonlinear aerodynamic characteristics of airfoil. The present method will contribute to the analysis of aircraft in the conceptual design because the present method can predict the nonlinear aerodynamic characteristics of a wing with a few computing resources and significant time. Copyright © 2007 John Wiley & Sons, Ltd.     

金属成型材料参数的反求技术

给出了快速准确地获得材料处于弹塑性大变形状态下各向异性弹塑性本构模型参数反求方法。首次提出了筛选试验测试点的活度规则,并以此来指导试验测试点的位置的选择;提出了仿真先验信息的概念,丰富了获取材料参数先验信息的途径;混合采用Levenberg—Marquardt方法和Gauss—Newton法的优化策略,给出了材料参数反求的基本公式和关键算法。数值算例表明,反求参数的初值以及反求区间的确定对于反求结果有着重要影响,为了确保反求过程的顺利进行,必须充分了解材料模型的先验信息．并充分利用筛选试验测试点的活度规则。同时效值算例计算还表明本文方法具有很高的计算精度和计算效率。     

Immersed boundary conditions method for unsteady flow problems described by the Laplace operator

An implicit, spectral algorithm for the analysis of unsteady flow problems governed by the Laplace operator in corrugated geometries is described. The algorithm treats the physical boundary conditions as constraints along lines internal to the solution domain. The method eliminates the need for coordinate generation and can be quickly adapted to changing geometries. Various tests confirm the spectral accuracy in space and the first‐ and second‐order accuracies in time. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Global linear and quadratic one-step smoothing newton method for vertical linear complementarity problems

IntroductionLetN ∈Rm0 ×n(m0 ≥n)beaverticalblockmatrixoftype(m1,… ,mn) ,q∈Rm0 beaconstantvectorpartitionedconformablywithN ,thatareN =N1Nn,　q =q1qn,　m0 = ni=1mi,whereNi ∈Rmi×n,qi ∈Rmi,i ∈I =1 ,… ,n .WeconsidertheverticallinearomplementarityproblemVLCP(N ,q)associatedwith (N ,q)offindingavectorx∈Rnsuchthatx≥ 0 ,　si(x) :=Nix+qi ≥ 0 ,　xi∏mij=1sij(x) =0　　 (i∈I) ,wherexi,sij(x)denotetheithcomponentofxandthejthcomponentofsi(x) ,respectively ;∏ mj=1ajdenotesa1…am.Iti…     

A higher-order boundary element method for three-dimensional potential problems

A method of eliminating the singularities involved in boundary element methods for three-dimensional potential problems is presented and the non-singular expressions of integrals on an element on which the singular point is situated are given for linear and quadratic interpolation functions. Numerical examples are compared with analytical solutions to show that the higher-order interpolations have better precision.     

求解非线性方程组的混沌分形方法

混沌分形是动力系统普遍出现的一种现象,牛顿-拉夫森NR(Newton-Raphson)方法是重要的一维及多维迭代技术,其迭代本身对初始点非常敏感,该敏感区是牛顿-拉夫森法所构成的非线性离散动力系统Julia集,在Julia集中迭代函数会呈现出混沌分形现象,提出了一种寻找牛顿-拉夫森函数的Julia点的求解方法,利用非线性离散动力系统在其Julia集出现混沌分形现象的特点,提出了一种基于牛顿-拉夫森法的非线性方程组求解的新方法,计算实例表明了该方法的有效性和正确性.