平面弹性问题的位移-应力混合重心插值配点法

提出数值分析平面弹性问题的位移-应力混合重心插值配点法。将弹性力学控制方程表达为位移和应力的耦合偏微分方程组,采用重心插值近似未知量,利用重心插值微分矩阵得到平面问题控制方程的矩阵形式离散表达式。使用重心插值离散位移和应力边界条件,采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,应用最小二乘法求解过约束方程组,得到平面弹性问题位移和应力数值解。数值算例结果表明,重心Lagrange插值方法的计算精度可达到10~(-10)量级。位移-应力混合重心插值配点法的计算公式简单、程序实施方便,是一种高精度的无网格数值分析方法。     

不可压缩平面问题的位移-压力混合重心插值配点法

引入人工压力变量,将弹性本构方程以应力、应变和压力表达,建立求解不可压缩平面弹性问题的位移-压力方程和不可压缩条件方程的耦合偏微分方程组。利用张量积型重心Lagrange插值近似二元函数,得到计算插值节点处偏导数的偏微分矩阵。采用配点法离散不可压缩弹性控制方程,利用偏微分矩阵直接离散弹性力学控制方程为矩阵形式方程组。利用插值公式离散位移和应力边界条件,将离散边界条件与离散控制方程组合为新的方程组,得到求解弹性问题的过约束线性代数方程组;利用最小二乘法求解线性方程组,得到弹性力学问题位移数值解。数值算例验证了所提方法的数值计算精度为10-14~10-10。     

轴对称圆薄板大挠度微分方程的数值分析方法

根据轴对称问题的特点,利用级数展开和求极限法则,证明了轴对称大挠度圆薄板在圆心处应满足的边界条件,并以圆薄板轴对称大挠度弯曲变形微分方程为基础,建立了圆心处非奇异的轴对称大挠度圆板弯曲微分方程,从而可以方便地利用现有的常微分方程数值求解方法(如变步长龙格-库塔法)对实心圆板的轴对称问题进行数值求解,又不必像摄动法那样推导复杂的公式。在数值求解轴对称圆板大挠度弯曲变形微分方程时,将非线性微分方程的求解主要归结为迭代求解圆心处三个未知边界条件的问题,即圆心处的径向膜力、圆心处的挠度、圆心处挠度的二阶导数,并提出了相应的求解方法。实例中,对于圆薄板受均布横向荷载的问题,分析了周边固支边界条件下的非线性弯曲问题,给出了中心挠度参数大范围变化时的荷载和部分边界值变化曲线,并与经典摄动解进行了对比。对比结果可见,本文方法和摄动法的解非常接近,在量纲归一化中心挠度不超过4.0时,两种方法解的相对误差均小于5.0%。另外,本文还分析了与挠度有关的液体压力作用下和集中荷载作用下周边固支圆板的非线性弯曲问题。通过算例可见:本文方法可以灵活处理不同的荷载问题;对于不同的问题,计算过程相似,不必推导复杂的计算公式,计算精度容易控制。     

非线性动力学常微分方程组高精度数值积分方法

建立了一种求解非线性动力学常微分方程组初值问题的新方法．若非线性函数一阶导数存在，则给出解的积分方程表达式，计算得到按规定误差要求的高精度数值解．引入一般自治或非自治非线性系统的首次近似Jacobi矩阵，不作任何假设重构等价的非线性常微分方程组，简捷而有广泛的适应性，不改变方程的本质，但其主项构成线性化方程组，其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算，给出该方程组初值问题的Duhamel卷积分解析表达式，在时间步长内进行数值积分选代求解，在指定误差内快速收敛，逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解．积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组，打破了传统用增量法在离散点上建立的代数方程组迭代求解，从而使传统Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法．数值计算中给出指数矩阵递增展开式，变矩阵乘法为乘积系数的加法，避免了大量矩阵自乘而大大提高计算效率．算法验证为无条件稳定，则保证对线性常微分方程而言，计算中舍入误差的传播不会扩散，不出现计算机字长有限而引起舍入误差导致计算不确定性问题．基于以上理论和数值方法，计算了线性非线性算例并进行了分析，验证了本方法简捷而有广泛的适应性，可以有足够的精确性．     

不规则区域平面弹性问题的正则区域重心插值配点法

将不规则区域嵌入到规则的矩形区域,在矩形区域上将弹性平面问题的控制方程采用重心Lagrange插值离散,得到控制方程矩阵形式的离散表达式。在边界节点上利用重心插值离散边界条件,规则区域采用置换法施加边界条件,不规则区域采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法进行求解,得到整个规则区域上的位移数值解。利用重心插值计算得到不规则区域内任意节点的位移值,计算精度可到10-14以上。数值算例验证了所建立方法的有效性和计算精度。     

不规则区域平面弹性问题的正则区域重心插值配点法

将不规则区域嵌入到规则的矩形区域,在矩形区域上将弹性平面问题的控制方程采用重心Lagrange插值离散,得到控制方程矩阵形式的离散表达式。在边界节点上利用重心插值离散边界条件,规则区域采用置换法施加边界条件,不规则区域采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法进行求解,得到整个规则区域上的位移数值解。利用重心插值计算得到不规则区域内任意节点的位移值,计算精度可到10^-14以上。数值算例验证了所建立方法的有效性和计算精度。     

反平面V形切口塑性应力奇异性分析

论文提出了用插值矩阵法计算幂硬化塑性材料反平面V形切口和裂纹尖端区域的应力奇异性.首先在切口和裂纹尖端区域采用自尖端径向度量的渐近位移场假设,将其代入塑性全量理论的基本微分方程后,推导出包含应力奇异性特征指数和特征角函数的非线性常微分方程特征值问题.然后采用插值矩阵法迭代求解导出的控制方程,得到一般的塑性材料反平面V形切口和裂纹的前若干阶应力奇异阶和相应的特征角函数,该法的重要优点是以上求解的特征角函数和它们各阶导函数具有同阶精度,并且一次性地求出前若干阶特征对.同时,插值矩阵法计算量小,易于和其他方法联合使用,这些优点在后续求解尖端区域完全应力场非常优越.论文方法的计算结果与现有结果对照,发现吻合良好,表明了论文方法的有效性.     

薄板和扁壳几何非线性分析的加权余量法

本文将加权余量法应用于薄板及扁壳的大挠度分析。在本法中将非线性微分方程组化为增量型的微分方程组,其中的非线性部分化作为折算荷载的形式用迭代法处理。在全部计算中,所求的解是线性方程组,系数矩阵只需组成一次。方程组的建立利用加权余量法中的最小二乘混合配点法。本法是一个求解非线性微分方程的通法,可称为最小二乘配点增量迭代法。文中附有算例,全部计算在微机上完成。     

边界节点法的计算稳定性研究

边界节点法利用满足控制方程的非奇异通解作为基函数,半解析边界数值离散偏微分方程,具有精度高、收敛快、易编程等优点,是一种纯无网格配点方法.但是在求解具体问题时,随着节点数的增加,边界节点法经常得到严重病态的插值矩阵.本文利用有效条件数评价边界节点法求解Helmholtz问题线性方程组的计算稳定性;然后利用三种正则化方法处理其病态的线性方程组,并与高斯消元法比较计算精度和收敛性.通过数值实验,本文研究了有效条件数、误差和正则化方法之间的关系.     

三阶非线性发展方程的线性化有理逼近方法

研究了三阶非线性发展方程的初边值问题的解。采用基于Sinc函数的微分求积法发展了线性化有理逼近方法。通常的配点法不适用于上述三阶问题的求解。本文把提出的方法用于求解KdV方程,取得了良好的效果。     

Variational iteration method for two-dimensional steady slip flow in micro-channels

In this paper, the two-dimensional steady slip flow in microchannels is investigated. Research on micro flow, especially on micro slip flow, is very important for designing and optimizing the micro electromechanical system (MEMS). The Navier-Stokes equations for two-dimensional steady slip flow in microchannels are reduced to a nonlinear third-order differential equation by using similarity solution. The variational iteration method (VIM) is used to solve this nonlinear equation analytically. Comparison of the result obtained by the present method with numerical solution reveals that the accuracy and fast convergence of the new method.     

A PREDICT-CORRECT NUMERICAL INTEGRATION SCHEME FOR SOLVING NONLINEAR DYNAMIC EQUATIONS

A new numerical integration scheme incorporating a predict-correct algorithm forsolving the nonlinear dynamic systems was proposed in this paper. A nonlinear dynamic systemgoverned by the equation v=F(v,t) was transformed into the form as v=Hv f(v,t). Thenonlinear part f(v,t) was then expanded by Taylor series and only the first-order term retained inthe polynomial. Utilizing the theory of linear differential equation and the precise time-integrationmethod, an exact solution for linearizing equation was obtained. In order to find the solution of theoriginal system, a third-order interpolation polynomial of v was used and an equivalent nonlinearordinary differential equation was regenerated. With a predicted solution as an initial value andan iteration scheme, a corrected result was achieved. Since the error caused by linearization couldbe eliminated in the correction process, the accuracy of calculation was improved greatly. Threeengineering scenarios were used to assess the accuracy and reliability of the proposed method andthe results were satisfactory.     

Optimal control of attitude for coupled-rigid-body spacecraft via Chebyshev-Gauss pseudospectral method

The attitude optimal control problem(OCP) of a two-rigid-body spacecraft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a discrete nonlinear programming(NLP) problem using the Chebyshev-Gauss pseudospectral method(CGPM). Solutions of the NLP problem can be obtained using the sequential quadratic programming(SQP) algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss(CG) points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform(FFT). To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and control variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.     

Investigation of Random Response of Rotational Shell When Considering Geometric Nonlinear Behaviour

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｅｎｇｉｎｅｅｒｉｎｇ,ｔｈｅｎｏｎｌｉｎｅａｒｉｔｙｏｆｓｈｅｌｌｓｔｒｕｃｔｕｒｅｉｓｓｏｍｅｔｉｍｅｓｏｂｖｉｏｕｓ.Ｉｔｉｓｖｅｒｙｄｉｆｆｉｃｕｌｔｔｏｓｏｌｖｅｔｈｅｐｒｏｂｌｅｍｏｆｔｈｅｎｏｎｌｉｎｅａｒｒａｎｄｏｍｖｉｂｒａｔｉｏｎ .Ｉｎｔｈｉｓｐａｐｅｒ,ｗｅｔｒｙｔｏｓｏｌｖｅｉｔｉｎｔｈｅｆｒｅｑｕｅｎｃｙｄｏｍａｉｎｂｙａｎｉｍｐｒｏｖｅｄｓｔａｔｉｓｔｉｃｌｉｎｅａｒｉｚａｔｉｏｎｍｅｔｈｏｄ .Ｗｅｃｏｎｓｉｄｅｒｏｎｌｙｔｈｅｇｅｏｍｅｔｒ…     

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

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

High-precision stress determination in photoelasticity

Stress separation is usually achieved by solving differential equations of equilibrium after parameter determination from isochromatics and isoclinics. The numerical error resulting from the stress determination is a main concern as it is always a function of parameters in discretization. To improve the accuracy of stress calculation,a novel meshless barycentric rational interpolation collocation method(BRICM) is proposed. The derivatives of the shear stress on the calculation path are determine...     

MHD Falkner-Skan flow of Maxwell fluid by rational Chebyshev collocation method

The magnetohydrodynamics （MHD） Falkner-Skan flow of the Maxwell fluid is studied. Suitable transform reduces the partial differential equation into a nonlinear three order boundary value problem over a semi-infinite interval. An efficient approach based on the rational Chebyshev collocation method is performed to find the solution to the proposed boundary value problem. The rational Chebyshev collocation method is equipped with the orthogonal rational Chebyshev function which solves the problem on the semi-infinite domain without truncating it to a finite domain. The obtained results are presented through the illustrative graphs and tables which demonstrate the affectivity, stability, and convergence of the rational Chebyshev collocation method. To check the accuracy of the obtained results, a numerical method is applied for solving the problem. The variations of various embedded parameters into the problem are examined.     

Fast numerical method for crack problem in the porous elastic material

The present paper discusses the crack problem in the linear porous elastic plane using the model developed by Nunziato and Cowin. With the help of Fourier transform the problem is reduced to an integral equation over the boundary of the crack. Some analytical transformations are applied to calculate the kernel of the integral equation in its explicit form. We perform a numerical collocation technique to solve the derived hyper-singular integral equation. Due to convolution type of the kernel, we apply, at each iteration step, the classical iterative conjugate gradient method in combination with the Fast Fourier technique to solve the problem in almost linear time. There are presented some numerical examples for materials of various values of porosity.     

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

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