Application of the extended Kantorovich method to the bending of clamped cylindrical panels

Applicability and performance of the extended Kantorovich method (EKM) to obtain highly accurate approximate closed form solution for bending analysis of a cylindrical panel is studied. Fully clamped panel subjected to both uniform and non-uniform loadings is considered. Based on the Love–Kirchhoff first approximation for thin shallow cylindrical panels, the governing equations of the problem in terms of three displacement components include a system of two second order and one forth order partial differential equations. The governing PDE system is converted to a double set of ODE systems by assuming separable functions for displacements together with utilization of the extended Kantorovich method. The resulted ODE systems are solved iteratively. In each iteration, exact closed form solutions are presented for both ODE systems. Rapid convergence and high accuracy of the method is shown for various examples. Both displacement and stress predictions show close agreement with other analytical and finite element analysis.     

多柔体系统数值分析的模型降噪方法

多柔体系统的动力学方程通常是一组刚性微分方程, 目前普遍采用的刚性微分方程数值解法主要通过数值阻尼滤除系统响应中的高频分量, 其求解效率难以令人满意. 为了降低多柔体系统动力学方程的刚性, 从而可采用ODE45等常规微分方程求解器进行求解, 研究了在建模过程中滤除高频振荡分量的方法. 在以当前时刻为起点的短时间内对柔性体的应力进行均匀化, 用均匀化后的应力计算柔性体的变形虚功率, 由此得到的系统动力学方程的解中不含过高频率的弹性振动, 并且可以通过调节均匀化时间区间的长度参数控制滤波的范围. 数值算例表明: 这种模型降噪方法的计算效率和精度均不低于刚性微分方程求解器, 并且在刚性微分方程求解器失效的情况下模型降噪方法仍有良好的精度和效率. 本文所提的模型降噪方法可成为求解多柔体系统动力学方程的新途径.      

An investigation of an Emden-Fowler equation from thin film flow

A third-order ordinary differential equation (ODE) for thin film flow with both Neumann and Dirichlet boundary conditions is transformed into a second-order nonlinear ODE with Dirichlet boundary conditions.Numerical solutions of the nonlinear second-order ODE are investigated using finite difference schemes.A finite difference formulation to an Emden-Fowler representation of the second-order nonlinear ODE is shown to converge faster than a finite difference formulation of the standard form of the second-order nonlinear ODE.Both finite difference schemes satisfy the von Neumann stability criteria.When mapping the numerical solution of the second-order ODE back to the variables of the original third-order ODE we recover the position of the contact line.A nonlinear relationship between the position of the contact line and physical parameters is obtained.     

ODE conversion techniques and their applications in computational mechanics

In this paper, a number of ordinary differential equation (ODE) conversion techniques for transformation of nonstandard ODE boundary value problems into standard forms are summarised, together with their applications to a variety of boundary value problems in computational solid mechanics, such as eigenvalue problem, geometrical and material nonlinear problem, elastic contact problem and optimal design problems through some simple and representative examples. The advantage of such approach is that various ODE boundary value problems in computational mechanics can be solved effectively in a unified manner by invoking a standard ODE solver. The project is supported by National Natural Science Foundation of China.     

变截面高层筒体结构的有限元线法整体稳定和二阶分析

本文用有限元线法对变截面的高层简体结构进行空间整体稳定和二阶分析。先把实际框筒结构分段连续化为正交各向异性折板结构；用有限元线法，通过有限条元半离散化，取结线上位移为基本未知函数，由势能驻值原理建立稳定和二阶分析的常微分方程组；由常微分方程求解器直接求解。     

常微分方程技术及其在固体力学计算中的应用

常微分方程(ODE——Ordinary Differential Equation)边值问题的最新计算求解技术的迅速发展推出了一批高质高效的通用软件,而工程中大量的ODE问题并非呈现为这些求解器(Solver)所接受的标准形式。然而,运用一些简单的ODE变换技巧可以将大量的不同类型的特殊问题转化为标准形式。本文列举了若干常用的变换技巧,并广泛地应用于各种固体力学问题的计算中,使大量的ODE问题在形式上得到统一,得以用标准的ODE Solver方便有效地求解。     

New solutions for surface tension driven spreading of a thin film

The standard fourth-order non-linear PDE modelling the flow of thin fluid film subject to surface tension is studied. The Lie group method is used to reduce the model equation from a fourth-order PDE to a fourth-order ODE. Analytical solutions are obtained for certain cases. Where analytical progress cannot be made, we determine numerical solutions.     

斜爆轰发动机的推力性能理论分析

爆轰燃烧具有释热快、循环热效率高的特点. 斜爆轰发动机利用斜爆轰波进行燃烧组织, 在高超声速吸气式推进系统中具有重要地位. 以往研究主要关注斜爆轰波的起爆、驻定以及波系结构等, 缺少从整体层面出发对斜爆轰发动机开展推力性能分析. 本文将斜爆轰发动机内的流动和燃烧过程分解成进气压缩、燃料掺混、燃烧释热和排气膨胀4个基本模块并分别进行理论求解, 建立了斜爆轰发动机推力性能的理论分析模型. 在斜爆轰波系研究成果的基础上, 选取了过驱动斜爆轰、Chapman?Jouguet斜爆轰、过驱动正爆轰和斜激波诱导等容燃烧等4种燃烧模式来描述燃烧室内的燃烧释热过程, 并对比分析了不同燃烧模式对发动机比冲性能的影响. 此外, 还获得了不同来流参数、燃烧室参数和进排气参数等对发动机推力的影响规律, 发现来流马赫数和尾喷管的膨胀面积比是发动机理论燃料比冲的主要影响因素. 最后, 结合以往关于受限空间内斜爆轰波驻定特性等方面的研究成果, 提出了斜爆轰发动机燃烧室的设计方向.      

圆柱扁壳的非线性动力响应的样条加权残值分析

本文对圆柱扁壳的非线性动力响应的问题,利用加相对加权残值法的子域方法将唐奈尔方程转化成时间因子作为示知函数的控制方程,运用样条配点法发成若干段,且保证在每时间段均连续将算例结果与有关文献比较的同时,还运用ＯＤＥ ＳＯＶＶＥＲ－ＣＯＬＳＹＳ程序作数值比较。     

An asymptotic iteration method for the numerical analysis of near-critical free-surface flows

Satisfying the boundary conditions at the free surface may impose severe difficulties to the computation of turbulent open-channel flows with finite-volume or finite-element methods, in particular, when the flow conditions are nearly critical. It is proposed to apply an iteration procedure that is based on an asymptotic expansion for large Reynolds numbers and Froude numbers close to the critical value 1.The iteration procedure starts by prescribing a first approximation for the free surface as it is obtained from solving an ODE that has been derived previously by means of an asymptotic expansion (Grillhofer and Schneider, 2003). The numerical solution of the full equations of motion then gives a surface pressure distribution that differs from the constant value required by the dynamic boundary condition. To determine a correction to the elevation of the free surface we next solve an ODE that is obtained from the asymptotic analysis of the flow with a prescribed pressure disturbance at the free surface. The full equations of motion are then solved for the corrected surface, and the procedure is repeated until criteria of accuracy for surface elevation and surface pressure, respectively, are satisfied.The method is applied to an undular hydraulic jump as a test case.     

Complex Inner Product Averaging Method for Calculating Normal Form of Ode

ＩｎｔｒｏｄｕｃｔｉｏｎＭｏｓｔｎｏｎｌｉｎｅａｒｄｙｎａｍｉｃｓｙｓｔｅｍｓａｒｅｍｕｌｔｉ＿ｄｅｇｒｅｅ＿ｏｆ＿ｆｒｅｅｄｏｍｓｙｓｔｅｍｓｉｎｔｈｅｆｉｅｌｄｏｆｅｎｇｉｎｅｅｒｉｎｇ[1,2 ],ｎｏｎｌｉｎｅａｒｄｙｎａｍｉｃｓｙｓｔｅｍｓｄｏｎｏｔｈａｖｅｅｘａｃｔｌｙｔｈｅｏｒｅｔｉｃｓｏｌｕｔｉｏｎｓ.Ｔｈｅｒｅｆｏｒｅ,ｉｔｉｓａｆｏｃａｌｐｏｉｎｔｃｏｎｃｅｒｎｅｄｂｙｓｃｉｅｎｔｉｓｔｓｔｏｓｅａｒｃｈｆｏｒａｐｐｒｏｘｉｍａｔｅｓｏｌｕｔｉｏｎｓｆｏｒａｌｏｎｇｔｉｍｅ ,ａｎｄｒ…     

Diffusion-driven instability and Hopf bifurcation in Brusselator system

The Hopfbifurcation for the Brusselator ordinary-differential-equation （ODE） model and the corresponding partial-differential-equation （PDE） model are investigated by using the Hopf bifurcation theorem. The stability of the Hopf bifurcation periodic solution is discussed by applying the normal form theory and the center manifold theorem. When parameters satisfy some conditions, the spatial homogenous equilibrium solution and the spatial homogenous periodic solution become unstable. Our results show that if parameters are properly chosen, Hopf bifurcation does not occur for the ODE system, but occurs for the PDE system.     

Self-similar singular solution of fast diffusion equation with gradient absorption terms

The self-similar singular solution of the fast diffusion equation with nonlinear gradient absorption terms are studied. By a self-similar transformation, the self-similar solutions satisfy a boundary value problem of nonlinear ordinary differential equation (ODE). Using the shooting arguments, the existence and uniqueness of the solution to the initial data problem of the nonlinear ODE are investigated, and the solutions are classified by the region of the initial data. The necessary and sufficient condition for the existence and uniqueness of self-similar very singular solutions is obtained by investigation of the classification of the solutions. In case of existence, the self-similar singular solution is very singular solution.     

分段常微分方程边值问题的改进COLSYS算法

对分段定义的常微分方程边值问题的 COLSYS算法进行了改进 ,提高了其计算效率及数值稳定性 ,并且可以节省内存 ,这些优势随着分段数的增加而迅速加大。     

变截面框架-剪力墙-薄壁筒结构弯扭耦连振动的常微分方程求解器解法

本文用沿高度方向分段连续化的方法,对沿高度方向为阶形变截面,在平面内为任意斜向布置的框架-剪力墙-薄壁筒协同工作体系,建立了弯扭耦连的振动方程,用常微分方程求解器(COLSYS)求解了自振频率及相应的振型。     

A Lagrangian Panel Method in the Time Domain for Moving Free-surface Potential Flows

A Lagrangian-type panel method in the time domain is proposed for potential flows with a moving free surface. After a spatial semi-discretization with a low-order scheme, the instantaneous velocity-potential and normal displacement on the moving free surface are obtained by means of a time-marching scheme. The kinematic and dynamic boundary conditions at the free surface are non-linear restrictions over the related Ordinary Differential Equation (ODE) system and, in order to handle them, an alternative Steklov-Poincaré operator technique is proposed. The method is applied to sloshing like flow problems.     

Flow of Oldroyd-B fluid with nanoparticles and thermal radiation

The two-dimensional boundary layer flow of an Oldroyd-B fluid in the presence of nanoparticles is investigated. Convective heat and mass conditions are considered in the presence of thermal radiation and heat generation. The Brownian motion and thermophoresis effects are retained. The nonlinear partial differential equations are reduced into the ordinary differential equation (ODE) systems. The resulting ODE systems are solved for the series solutions. The results are analyzed for various physical parameters of interest. Numerical values of the local Nusselt and Sherwood numbers are also computed and analyzed.     

奇异Hamilton系统矩阵的精细积分法

常规位移有限元的结构振动方程是n个二阶常微分方程组.采用一般交分原理推导,将结构振动问题引入Hamiltoil体系,将得到2n个一阶常微分方程组.精细积分法宜于处理一阶方程,应用于线性定常结构动力问题求解,可以得到在数值上逼近精确解的结果.对于非齐次动力方程,当结构具有刚体位移时,系统矩阵将出现奇异.本文借鉴全元选大元高斯-约当法求解线性方程组的经验,提出全元选大元法求奇异矩阵零本征解的方法,该方法可以简便快速地寻求奇异矩阵零本征值对应的子空间.利用Hamiltoil体系已有研究成果及Hamilton系统的共轭辛正交归一关系,迅速将零本征值对应的子空间分离出来,通过投影排除奇异部分,然后用精细积分法求得问题的解.数值算例表明,该方法对Hamilton系统奇异问题,处理方便,计算量小,易于实现,同时保持了精细算法的优点.     

具有待定形函数的有限元线法

有限元线法是一种优良的半解析法，但其解函数存在解析方向和离散方向上的精度不相称的问题。本文提出了一种基于有限元线法的双向解析法——具有待定形函数的有限元线法。该方法设形函数为未知待定并借助能量原理对形函数进行解析性求解，使解函数在两个方向上的精度达到相称，从而大幅度提升了解答的整体质量和精度。文中，以二维Poisson方程问题为例，具体给出了待定形函数的有限元线法的算法和算例，用以表明本法的可行性和有效性。