二维数字图像相关法在准静态拉伸实验中的应用

对二维数字图像相关法在准静态拉伸实验中的应用进行了研究。二维数字图像相关法采用Newton-Raphson方法和逆匹配结合的算法,首先对Newton-Raphson迭代式中的相关函数的一阶偏导和二阶偏导进行了简化,对于Newton-Raphson法的初值估计采用仿射变换方法和RG-DIC结合的方法;然后在计算域全场采用双三次插值,保证计算精度的同时降低了计算量。散斑仿真实验和准静态拉伸实验显示,二维数字图像相关法获得的分析值与应变片测量值相对误差约小于10%,表明二维DIC算法具有较高的准确性,能够用于实际工程测量。     

优化迭代步长的两种改进增量谐波平衡法

增量谐波平衡法(IHB法)是一个半解析半数值的方法, 其最大优点是适合于强非线性系统振动的高精度求解. 然而, IHB法与其他数值方法一样, 也存在如何选择初值的问题, 如初值选择不当, 会存在不收敛的情况. 针对这一问题, 本文提出了两种基于优化算法的IHB法: 一是结合回溯线搜索优化算法(BLS)的改进IHB法(GIHB1), 用来调节IHB法的迭代步长, 使得步长逐渐减小满足收敛条件; 二是引入狗腿算法的思想并结合BLS算法的改进IHB法(GIHB2), 在牛顿-拉弗森(Newton-Raphson)迭代中引入负梯度方向, 并在狗腿算法中引入2个参数来调节BSL搜索方式用于调节迭代的方式, 使迭代方向沿着较快的下降方向, 从而减少迭代的步数, 提升收敛的速度. 最后, 给出的两个算例表明两种改进IHB法在解决初值问题上的有效性.      

弹塑性体脆断相场模型的一致性无单元伽辽金法

采用无单元伽辽金法（EFG）对弹塑性体脆性断裂的相场模型进行了数值实现。利用无单元法便于构建高阶近似函数的优势,位移和相场均采用二阶移动最小二乘（MLS）近似。刚度阵的数值积分采用更为高效的二阶一致三点积分格式QC3（Quadratically Consistent 3-point integration scheme）。本构算法采用Newton-Raphson迭代和弹塑性一致性切线模量。数值结果表明了本文方法模拟弹塑性体脆性断裂的有效性。     

一种求解柔性多体系统动力学方程的新方法

柔性多体系统控制方程是具有stiff性质的刚柔耦合非线性代数一微分方程组,本文提出了一种求解该类刚性方程组的数值方法,在每一时间步,利用Newmark-β直接积分法计算迭代初值,基于控制方程及约束方程的泰勒展开,推导出Newton-Raphson迭代公式,对位移及拉格朗日乘子进行修正,最后,引用Blajer提出的违约修正方法对数值积分过程中约束方程的违约进行修正。就两个典型算例进行了数值仿真,结果证明了本文方法的有效性。     

基于时序变形预测的数字图像相关加速方法

利用数字图像相关求解连续变形物体的位移场和应变场时,会遇到处理速度非常慢的问题,原因是相关算法计算量大且忽略了物体在时间轴上的变形规律。本文提出了一种基于变形预测的数字图像相关方法,该方法利用物体在时间轴上的变形规律,通过已经得到的变形值来预测后面时刻的变形初值。首先分析了物体在时间轴上的变形规律,然后结合实际应用对其进行修正,得到物体下一步变形的初值,最后通过NR方法(Newton-Raphson Method)得到物体的位移场和应变场。通过计算机模拟和金属试件拉伸实验并结合GPU(Graphic Processing Unit)编程验证了该方法的有效性。计算机模拟和金属试件拉伸实验的计算结果表明,该方法能够在计算精度保持不变的情况下,使计算速度提高4~7倍。     

重载滑动轴承非牛顿流体的三维弹性液体动力润滑分析的加权平均法

本文采用加权平均法成功地求解了由广义Ｒｅｙｎｏｌｄｓ方程（幂次型）和三维弹性方程组成的非牛顿体弹流模型，该方法较好的解决了大偏心率的弹流迭代收敛困难的问题。该方法具有收敛快，对初值要求松且计算精度高等优点，是一种解决弹流问题较好的迭代模式。     

概率密度演化方程TVD格式的自适应时间步长技术及其初值条件改进

随机性普遍存在于实际工程问题中,而复杂结构的非线性随机响应分析是其中的一个难点,近年发展的概率密度演化方法为此类问题的求解提供了新的途径.由于实际问题的复杂性,概率密度演化方程通常采用数值方法求解,因此提高计算效率和求解精度对实际应用具有重要意义.本文基于变网格技术,推导了概率密度演化方程在非均匀时间步长上的总变差减小(total variation diminishing,TVD)差分格式,算例结果表明通过自适应插值可将迭代次数减少为原来的43.4%,当随机过程样本持续时间增大时均值估计的平均误差基本不变,而标准差估计的平均误差不断增大,但增大幅度不断减小;计算耗时随样本持续时间的增大也呈增大趋势,而由于使用了时间步长自适应插值算法导致有些情况下长持时样本的计算耗时反而比短持时样本的计算耗时短;在传统的脉冲函数型初值条件基础上,提出了一种高阶导数更稳定的余弦函数型初值条件形式.结果表明,脉冲函数型的初值条件是余弦函数型初值条件的一个特例,当参数取值适当时,余弦函数型初值条件的数值求解结果具有更高的精度.本文的工作进一步完善了概率密度演化方程的求解方法,为其在实际工程中的应用提供了基础.     

概率密度演化方程TVD格式的自适应时间步长技术及其初值条件改进

随机性普遍存在于实际工程问题中,而复杂结构的非线性随机响应分析是其中的一个难点,近年发展的概率密度演化方法为此类问题的求解提供了新的途径.由于实际问题的复杂性,概率密度演化方程通常采用数值方法求解,因此提高计算效率和求解精度对实际应用具有重要意义.本文基于变网格技术,推导了概率密度演化方程在非均匀时间步长上的总变差减小(total variation diminishing,TVD)差分格式,算例结果表明通过自适应插值可将迭代次数减少为原来的43.4%,当随机过程样本持续时间增大时均值估计的平均误差基本不变,而标准差估计的平均误差不断增大,但增大幅度不断减小;计算耗时随样本持续时间的增大也呈增大趋势,而由于使用了时间步长自适应插值算法导致有些情况下长持时样本的计算耗时反而比短持时样本的计算耗时短;在传统的脉冲函数型初值条件基础上,提出了一种高阶导数更稳定的余弦函数型初值条件形式.结果表明,脉冲函数型的初值条件是余弦函数型初值条件的一个特例,当参数取值适当时,余弦函数型初值条件的数值求解结果具有更高的精度.本文的工作进一步完善了概率密度演化方程的求解方法,为其在实际工程中的应用提供了基础.      

重载滑动轴承非牛顿流体的三维弹性流体动力润滑分析的加权平均法

本文采用加权平均法成功地求解了由广义Ｒｅｙｎｏｌｄｓ方程（幂次型）和三维弹性方程组成的非牛顿体弹流模型，该方法较好的解决了大偏心率下的弹流迭代收敛困难的问题。该方法具有收敛快，对初值要求松且计算精度高等优点，是一种解决弹流问题较好的迭代模式。     

初值约束与两点边值约束轨道动力学方程的快速数值计算方法

轨道动力学快速计算是航天工程中的基础问题,广泛存在于轨道设计、空间抓捕以及深空探测等任务中.基于有限差分原理的经典数值积分算法,由于精度严重依赖小积分步长,难以满足航天器在轨快速计算需求.针对该问题,提出一种局部配点反馈迭代算法,该算法能高效解算受到初值约束和两点边值约束的轨道动力学方程.基于Picard迭代公式建立数...     

Fast, Robust and Accurate Digital Image Correlation Calculation Without Redundant Computations

High-efficiency and high-accuracy deformation analysis using digital image correlation (DIC) has become increasingly important in recent years, considering the ongoing trend of using higher resolution digital cameras and common requirement of processing a large sequence of images recorded in a dynamic testing. In this work, to eliminate the redundant computations involved in conventional DIC method using forward additive matching strategy and classic Newton–Raphson (FA-NR) algorithm without sacrificing its sub-pixel registration accuracy, we proposed an equivalent but more efficient DIC method by combining inverse compositional matching strategy and Gauss-Newton (IC-GN) algorithm for fast, robust and accurate full-field displacement measurement. To this purpose, first, an efficient IC-GN algorithm, without the need of re-evaluating and inverting Hessian matrix in each iteration, is introduced to optimize the robust zero-mean normalized sum of squared difference (ZNSSD) criterion to determine the desired deformation parameters of each interrogated subset. Then, an improved reliability-guided displacement tracking strategy is employed to achieve further speed advantage by automatically providing accurate and complete initial guess of deformation for the IC-GN algorithm implemented on each calculation point. Finally, an easy-to-implement interpolation coefficient look-up table approach is employed to avoid the repeated calculation of bicubic interpolation at sub-pixel locations. With the above improvements, redundant calculations involved in various procedures (i.e. initial guess of deformation, sub-pixel displacement registration and sub-pixel intensity interpolation) of conventional DIC method are entirely eliminated. The registration accuracy and computational efficiency of the proposed DIC method are carefully tested using numerical experiments and real experimental images. Experimental results verify that the proposed DIC method using IC-GN algorithm and the existing DIC method using classic FA-NR algorithm generate similar results, but the former is about three to five times faster. The proposed reliability-guided IC-GN algorithm is expected to be a new standard full-field displacement tracking algorithm in DIC.     

Nonlinear thermo-elastic buckling characteristics of cross-ply laminated joined conical–cylindrical shells

Here, the nonlinear thermo-elastic buckling/post-buckling characteristics of laminated circular conical–cylindrical/conical–cylindrical–conical joined shells subjected to uniform temperature rise are studied employing semi-analytical finite element approach. The nonlinear governing equations, considering geometric nonlinearity based on von Karman’s assumption for moderately large deformation, are solved using Newton–Raphson iteration procedure coupled with displacement control method to trace the pre-buckling/post-buckling equilibrium path. The presence of asymmetric perturbation in the form of small magnitude load spatially proportional to the linear buckling mode shape is assumed to initiate the bifurcation of the shell deformation. The study is carried out to highlight the influences of semi-cone angle, material properties and number of circumferential waves on the nonlinear thermo-elastic response of the different joined shell systems.     

纺织物下垂变形数值计算

利用共旋技术提出了一个高柔性梁（布条）大位移大转动计算模型,通过把纺织物沿经向和纬向离散成许多布条,建立了纺织物下垂变形模型,模型仅使用了点平移自由度而没有使用转动自由度,并采用了一种单元应变能积分格式,模型考虑了纺织物面内拉伸、面内剪切以及弯曲;提出了一种控制单步全Newton-Raphson迭代格式步长的方法来确保得到控制方程组的收敛解。数值算例表明模型能够较精确地预测纺织物的最终下垂形状。     

Iteration schemes for rapid post‐stall aerodynamic prediction of wings using a decambering approach

Nonlinear aerodynamics of wings may be evaluated using an iterative decambering approach. In this approach, the effect of flow separation due to stall at any wing section is modeled as an effective reduction in section camber. The approach uses a wing analysis method for potential‐flow calculations and viscous airfoil lift curves for the sections as input. The calculation procedure is implemented using a Newton–Raphson iteration to simultaneously satisfy the boundary condition, which comes from potential‐flow wing theory, and drive the sectional operating points toward their respective viscous lift curves, as required for convergence. Of particular interest in this research is the calculation of the residuals during the Newton iteration. Unlike a typical implementation of the Newton iteration, the residual calculation is not performed via a straightforward function evaluation, but rather by estimating the target operating points on the input viscous lift curves. Estimation of these target operating points depends on the assumptions made in the cross‐coupling of the decambering at the different sections. This paper presents four residual calculation schemes for the decambering approach. The residual calculation schemes are compared against each other to assess computational speed and robustness. Decambering results are also compared with higher‐order computational fluid dynamics (CFD) solutions for rectangular and swept wings. Results from the best scheme compare well with the CFD solutions for the rectangular wing, motivating further development of the method. Poor predictions for the swept wings are traced to spanwise propagation of separated flow at stall, highlighting the limitations of the current approach. Copyright © 2014 John Wiley & Sons, Ltd.     

Comparison of Iterative Methods for Improved Solutions of the Fluid Flow Equation in Partially Saturated Porous Media

Abstract. The Picard and modified Picard iteration schemes are often used to numerically solve the nonlinear Richards equation governing water flow in variably saturated porous media. While these methods are easy to implement, they are only linearly convergent. Another approach to solve the Richards equation is to use Newton's iterative method. This method, also known as Newton–Raphson iteration, is quadratically convergent and requires the computation of first derivatives. We implemented Newton's scheme into the mixed form of the Richards equation. As compared to the modified Picard scheme, Newton's scheme requires two additional matrices when the mixed form of the Richards equation is used and requires three additional matrices, when the pressure head-based form is used. The modified Picard scheme may actually be viewed as a simplified Newton scheme.Two examples are used to investigate the numerical performance of different forms of the 1D vertical Richards equation and the different iterative solution schemes. In the first example, we simulate infiltration in a homogeneous dry porous medium by solving both, the h based and mixed forms of Richards equation using the modified Picard and Newton schemes. Results shows that, very small time steps are required to obtain an accurate mass balance. These small times steps make the Newton method less attractive.In a second test problem, we simulate variable inflows and outflows in a heterogeneous dry porous medium by solving the mixed form of the Richards equation, using the modified Picard and Newton schemes. Analytical computation of the Jacobian required less CPU time than its computation by perturbation. A combination of the modified Picard and Newton scheme was found to be more efficient than the modified Picard or Newton scheme.     

Nonlinear analysis of underconstrained structures

Underconstrained (kinematically indeterminate) assemblies of bars and pin-joints possess a specific kind of nonlinearity which appears even for small displacements and linear elasticity. Various approaches to nonlinear analysis based on the Newton–Raphson procedure are considered. The subspace Newton–Raphson technique is proposed. Theoretical considerations are accompanied by numerical examples of plane and space underconstrained assemblies.     

Newton iterative identification for a class of output nonlinear systems with moving average noises

This paper discusses iterative identification problems for a class of output nonlinear systems (i.e., Wiener nonlinear systems) with moving average noises from input–output measurement data, based on the Newton iterative method. The basic idea is to decompose a nonlinear system into two subsystems, to replace the unknown variables in the information vectors with their corresponding estimates at the previous iteration, and to present a Newton iterative identification method using the hierarchical identification principle. The numerical simulation results indicate that the proposed algorithms are effective.     

Mixed finite element for two-dimensional incompressible convective Brinkman-Forchheimer equations

In this work, the two-dimensional convective Brinkman-Forchheimer equations are considered. The well-posedness for the variational problem and its mixed finite element approximation is established, and the error estimates based on the conforming approximation are obtained. For the computation, a one-step Newton(or semi-Newton)iteration algorithm initialized using a fixed-point iteration is proposed. Finally, numerical experiments using a Taylor-Hood mixed element built on a structured or unstructured triangular mesh are implemented. The numerical results obtained using the algorithm are compared with the analytic data, and are shown to be in very good agreement. Moreover,the lid-driven problem at Reynolds numbers of 100 and 400 is considered and analyzed.     

Digital Image Correlation Using Stochastic Parallel-Gradient-Descent Algorithm

This article proposes a digital image correlation (DIC) method based on the stochastic parallel gradient descent (SPGD) algorithm. Stochastic parallel perturbations are imposed on deformation parameters to make the correlation coefficients converge to a global extremum; thus, this allows the final measured values of the deformation parameters to be obtained and the DIC measurement to be made. Both simulated and real data processing, including rigid body and strain deformation, show that the proposed method can achieve nearly the same accuracy as the Newton–Raphson (NR) method in most cases and higher accuracy in some cases, such as the simulated experiments of rigid body translation with and without noise. It also has a good noise-robustness. Furthermore, a series of experiments have been designed to evaluate the convergence characteristics of the proposed method, and it has been proved able to process large displacement and have a stable convergence process, good robustness, and a high convergence speed when bilinear interpolation is adopted.     

Implicit finite element methods for two-phase flow in oil reservoirs

The equations governing immiscible, incompressible, two-phase, porous media flow are discretized by generalized streamline diffusion Petrov–Galerkin methods in space and by implicit differences in time. Systems of non-linear algebraic equations are solved by Newton–Raphson iteration employing ILU-preconditioned conjugate-gradient-like methods to the non-symmetric matrix system in each iteration. The resulting solution methods are robust, enable complex grids with irregular nodal orderings and allow capillary effects. Several numerical formulations are tested and compared for one-, two- and three-dimensional flow cases, with emphasis on problems involving saturation shocks, heterogeneous media and curved boundaries. For reservoirs consisting of multiple rock types with differing capillary pressure properties, it is shown that traditional Bubnov-Galerkin methods give poor results and the new Petrov–Galerkin formulations are required. Investigations regarding the behaviour of several preconditioned conjugate-gradient-like methods in these type of problems are also reported.