非均匀圆柱壳非线性轴对称变形一般解

本文利用阶梯折算法[1],得到了非均匀圆柱壳非线性轴对称变形的一般解.文中导出了在任意轴对称载荷下求解非均匀圆柱壳非线性弯曲的位移和内力的一般公式,并给出一致收敛于精确解的证明.问题最后归结为求解二元一次代数方程组,文末给出算例.算例表明,无论内力和位移都可得到满意的结果,并收敛于精确解.     

非均匀双向加肋圆柱壳非线性轴对称变形的精确解析法

本文首次利用精确解析法分析了环向和纵向加肋非均匀圆柱壳在任意载荷和边界条件下非线性轴对称变形问题.导出了一致收敛于精确解的位移和内力解析表达式,文中给出收敛性问题.问题最后归结为求解二元一次代数方程组,计算既简便又迅速.文末给出四个数值算例表明,本文提出的方法,可以得到满意的结果.     

非均匀弹性地基圆板轴对称弯曲的一般解

本文在阶梯折算法的基础上,提出一个新的方法——精确解析法,得到了非均匀弹性地基圆板弯曲的一般解.文中导出了在任意轴对称载荷和边界条件下求解非均匀弹性地基圆板和中心带孔圆板弯曲的一般公式,并给出一致收敛于精确解的证明.文中得到的一般解可直接计算无弹性地基圆板的弯曲问题.问题最后归结为求解一个二元一次代数方程.文末给出算例,算例表明无论内力和位移均可得到满意的结果.     

非均匀变截面梁动力响应的一般解

本文利用精确解析法[1]给出非均匀变截面梁在任意谐振荷载和边界条件下的动力响应的一般解.问题最后归结为求解一个非正定微分方程.对于这一问题用一般变分法求解失败,利用本文的方法,这个问题的一般解可以写成解析的形式.因此对优化特别方便.本文给出收敛性证明.文末给出的算例表明本文的方法可获得满意的结果.     

任意变系数微分方程的精确解析法

工程中的许多问题归结为求解任意变系数微分方程的解.本文首次提出精确解析法,用以求解任意变系数微分方程在任意边界条件下的解.文中还给出精确解析法的一般计算格式,得到了一致收敛于精确解及其任意阶导数的解析表达式,并给出收敛性证明.文末给出四个算例,均得到较好的结果,证明了本文理论的正确性.     

求解厚矩形板弯曲问题的功的互等定理法

在本文中,功的互等定理法(RTM)被推广于求解基于Reissner理论的厚矩形板弯曲问题。首先,本文绘出了厚矩形板弯曲的基本解;其次,给出了三边固定一边自由在均布载荷作用下厚矩形板弯曲的精确解析解;最后,我们分析了本文解的数值结果。     

求解多项式规划的一个全局最优化算法

提出一个求解带箱子约束的一般多项式规划问题的全局最优化算法, 该算法包含两个阶段, 在第一个阶段, 利用局部最优化算法找到一个局部最优解. 在第二阶段, 利用一个在单位球上致密的向量序列, 将多元多项式转化为一元多项式, 通过求解一元多项式的根, 找到一个比当前局部最优解更好的点作为初始点, 回到第一个 阶段, 从而得到一个更好的局部最优解, 通过两个阶段的循环最终找到问题的全局最优解, 并给出了算法收敛性分析. 最后, 数值结果表明了算法是有效的.     

非均匀Reissner板弯曲的精确元法

本文在阶梯折算法和精确解析法的基础上,提出构造有限元的新方法——精确元法.该方法不用变分原理,可适用于任意变系数正定和非正定偏微分方程.利用该方法,得到Reissuer板弯曲的一个非协调单元,它具有十五个自由度.由于节点位移参数仅含有挠度和转角,因此处理任意边界条件非常容易.文中给出证明,位移和内力均收敛于精确解.由精确元法所得到的单元不仅能用于厚板,也可用于薄板.文末给出四个算例.算例表明,利用本文的方法,可获得满意的结果,并有较高的数值精度.     

线性比式和优化问题的完全多项式时间近似算法

本文针对线性比式和优化问题提出一个完全多项式时间近似算法,该算法主要利用原问题的等价问题及网格结点参数获得有限个与结点参数相关的线性规划问题,通过求解这些线性规划问题获得原问题的近似最优解.最终证明算法的收敛性,并给出了算法的计算复杂度,通过计算结果呈现算法的有效性与可行性.     

一类带参数的非线性奇异摄动问题的自适应移动网格算法

本文讨论一类带参数的非线性奇异摄动问题的自适应移动网格方法.首先,在任意非均匀网格下,利用向后欧拉公式对方程进行离散,并给出相应的局部截断误差.然后,基于局部截断误差和网格等分布原理,利用精确解的弧长函数,证明半离散格式下自适应移动网格算法是一阶收敛的.同时,基于近似的弧长控制函数,给出易于实现的网格生成算法,并给出全离散格式下的后验误差估计.最后,数值实验结果验证了本文所给出的理论结果.     

一个高精度收敛的变系数微分方程精确解析法*

文[1]给出精确解析法,可用于求解任意变系数微分方程,所得到的解具有二阶收敛精度.在此基础上,本文以变截面梁弯曲为例,给出一个高精度的算法.不增加工作量的情况下可达到四阶收敛精度.具有计算快,简单等特点,文末给出算例,仅用很少的单元即可获得高的收敛精度,表明了本文理论的正确性.     

Effective computation of exact and analytic approximate solutions to singular nonlinear equations of Lane–Emden–Fowler type

The particular motivation of this work is to develop a computational method to calculate exact and analytic approximate solutions to singular strongly nonlinear initial or boundary value problems of Lane–Emden–Fowler type which model many phenomena in mathematical physics and astrophysics. A powerful algorithm is proposed based on the series representation of the solution via suitable base functions. The utilization of such functions converts the solution of a given nonlinear differential equation to the solution of algebraic equations. Error analysis and convergence of the method is presented. Comparisons with the other methods reveal validity, applicability and great potential of the method. Several physical problems are treated to illustrative the good performance and high accuracy of the technique.     

大气运动基本方程组的解析解

在已知大气运动基本方程于光滑函数类中具有最好的稳定性前提下,讨论了它的局部解的解空间构造．根据它的解空间构造,分析了这个方程具有代表性和应用性的第三初值问题,在解析函数类中给出了适定的第三初值问题的解析解的计算方法以及具体的关系表达式,在局部解意义下完整的解决了这一点初值问题的解析解所涉及的理论与计算问题．指出其它类型定解问题都可以仿照文中的计算方法和步骤,求出所需要的稳定的解析解．     

A new cubic B-spline method for linear fifth order boundary value problems

In this paper, we mainly study the numerical solution of linear fifth order boundary value problems by using cubic B-splines. Our algorithm develops not only the cubic spline approximation solution but also the approximation derivatives of first order to fourth order of the analytic solution at the same time. This new method has lower computational cost than many other methods and is second order convergent. Numerical examples are given to demonstrate the effectiveness of our method.     

Exact solutions for non-linear Schrödinger equations by differential transformation method

In this paper, we implemented relatively new, exact series method of solution known as the differential transform method for solving linear and non-linear Schrödinger equations with initial profile. The method can easily be applied to many linear and non-linear problems and is capable of reducing the size of computational work. Exact solutions can also be achieved by the known forms of the series solutions. The results obtained are in good agreement with the exact solution. These results show that the technique introduced here is accurate and easy to apply.     

Implicit algorithms for solving the Cauchy problem for the equations of the dynamics of mechanical systems

It is shown that in the numerical solution of the Cauchy problem for systems of second-order ordinary differential equations, when solved for the highest-order derivative, it is possible to construct simple and economical implicit computational algorithms for step-by-step integration without using laborious iterative procedures based on processes of the Newton-Raphson iterative type. The initial problem must first be transformed to a new argument — the length of its integral curve. Such a transformation is carried out using an equation relating the initial parameter of the problem to the length of the integral curve. The linear acceleration method is used as an example to demonstrate the procedure of constructing an implicit algorithm using simple iterations for the numerical solution of the transformed Cauchy problem. Propositions concerning the computational properties of the iterative process are formulated and proved. Explicit estimates are given for an integration stepsize that guarantees the convergence of the simple iterations. The efficacy of the proposed procedure is demonstrated by the numerical solution of three problems. A comparative analysis is carried out of the numerical solutions obtained with and without parametrization of the initial problems in these three settings. As a qualitative test the problem of the celestial mechanics of the “Pleiades” is considered. The second example is devoted to modelling the non-linear dynamics of an elastic flexible rod fixed at one end as a cantilever and coiled in its initial (static) state into a ring by a bending moment. The third example demonstrates the numerical solution of the problem of the “unfolding” of a mechanical system consisting of three flexible rods with given control input.