首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
非线性动力学积分方程分块积分解法   总被引:2,自引:0,他引:2  
对于非线性动力学方程组分块地应用精细积分算法,使其化成积分方程表达式,求解的表达式中具有相对低阶的转换矩阵,从而使精细积分更适用于多自由度、强非线性、变系数、非保守系统,针对积分方程提出了一个显示预测-校正的单步四阶精度自起步的精细积分算法。算例表明本方法是有效的。  相似文献   

2.
结构非线性动力方程的精细积分算法   总被引:16,自引:0,他引:16  
基于线性方程精细积分的思路,对具有惯性、阻尼、刚度非线性的动力方程及参变非线性动力方程提出了一种较高精度线性化精细积分迭代计算算法,算例表明该算法可用较大的步长取得满意的计算精度,并可在较大的线性化区间获得较高的计算精度。  相似文献   

3.
证明面力边界积分方程被积函数的散度等于零,应用Stokes公式,对平面线弹性问题,将面力边界积分的求解转化为边界点的位移势函数的点值计算。应用边界积分方程的求解结果,推导出J积分亦可表示为边界点的积分势函数的点值计算,无需进行数值积分,实例计算说明该方法的有效性。  相似文献   

4.
位移导数边界积分方程一直存在着超奇异积分计算的障碍,该文提出以符号算子δye和εye作用于位移导数边界积分方程,施用一系列变换将边界位移、面力和位移导数转成为新的边界张量,从而得到一个新的边界积分方程--自然边界积分方程,自然边界积分方奇异性为强奇性,文中给出了相应的Cauchy主值积分算式,自然边界积分方程与位移边界积分方程联合可直接获取边界应力,几个算例表明了自然边界积分方程的正确性。  相似文献   

5.
边界积分方程中超奇异积分的解法   总被引:4,自引:0,他引:4  
董春迎  谢志成 《力学进展》1995,25(3):424-429
本文对边界积分方程中所存在的超奇异积分的数值解法作了综述,并介绍了它的一些应用。  相似文献   

6.
基于精细积分技术的非线性动力学方程的同伦摄动法   总被引:2,自引:0,他引:2  
将精细积分技术(PIM)和同伦摄动方法(HPM)相结合,给出了一种求解非线性动力学方程的新的渐近数值方法。采用精细积分法求解非线性问题时,需要将非线性项对时间参数按Taylor级数展开,在展开项少时,计算精度对时间步长敏感;随着展开项的增加,计算格式会变得越来越复杂。采用同伦摄动法,则具有相对筒单的计算格式,但计算精度较差,应用范围也限于低维非线性微分方程。将这两种方法相结合得到的新的渐近数值方法则同时具备了两者的优点,既使同伦摄动方法的应用范围推广到高维非线性动力学方程的求解,又使精细积分方法在求解非线性问题时具有较简单的计算格式。数值算例表明,该方法具有较高的数值精度和计算效率。  相似文献   

7.
非线性粘弹理论中的单积分型本构关系   总被引:7,自引:0,他引:7  
本文综述了非线性粘弹理论中的单积分本构表达,评述了多种有代表性的单积分型非线性粘弹理论,对几种本构方程加以分析比较,以揭示它们的内涵,明了其非线性表述原理。  相似文献   

8.
为了对等离子体密度非线性扩散方程进行数值分析,首先以差商代替微商,然后用对应方程的基本解,使之化为积分方程、再用边界单元法求解.  相似文献   

9.
提出了间接求解传统Helmholtz边界积分方程CBIE的强奇异积分和自由项系数,以及Burton-Miller边界积分方程BMBIE中的超强奇异积分的特解法。对于声场的内域问题,给出了满足Helmholtz控制方程的特解,间接求出了CBIE中的强奇异积分和自由项系数。对于声场外域对应的BMBIE中的超强奇异积分,按Guiggiani方法计算其柯西主值积分需要进行泰勒级数展开的高阶近似,公式繁复,实施困难。本文给出了满足Helmholtz控制方程和Sommerfeld散射条件的特解,提出了间接求出超强奇异积分的方法。推导了轴对称结构外场问题的强奇异积分中的柯西主值积分表达式,并通过轴对称问题算例证明了本文方法的高效性。数值结果表明,对于内域问题,采用本文特解法的计算结果优于直接求解强奇异积分和自由项系数的结果,且本文的特解法可避免针对具体几何信息计算自由项系数,因而具有更好的适用性。对于外域问题,两者精度相当,但本文的特解法可避免对核函数进行高阶泰勒级数展开,更易于数值实施。  相似文献   

10.
一类非线性奇异积分方程及其数值方法研究   总被引:1,自引:0,他引:1  
探讨了一类非线性奇异积分方程的数理性质以及在颗粒雷诺数Rep<1时此类方程解的存在条件,然后详细研究了该方程的数值计算方法并构造称之为P(EC)^k多步法的差分格式,分析了该格式的收敛性和代数精度,得到时域离散步长的约束关系。运用该格式计算了静止流场和均匀振荡流中球形小颗粒的非恒定运动,将计算结果与其解析解及有关实验数据的比较表明,上述数值方法具有良好的计算精度。  相似文献   

11.
A semi‐implicit, staggered finite volume technique for non‐hydrostatic, free‐surface flow governed by the incompressible Euler equations is presented that has a proper balance between accuracy, robustness and computing time. The procedure is intended to be used for predicting wave propagation in coastal areas. The splitting of the pressure into hydrostatic and non‐hydrostatic components is utilized. To ease the task of discretization and to enhance the accuracy of the scheme, a vertical boundary‐fitted co‐ordinate system is employed, permitting more resolution near the bottom as well as near the free surface. The issue of the implementation of boundary conditions is addressed. As recently proposed by the present authors, the Keller‐box scheme for accurate approximation of frequency wave dispersion requiring a limited vertical resolution is incorporated. The both locally and globally mass conserved solution is achieved with the aid of a projection method in the discrete sense. An efficient preconditioned Krylov subspace technique to solve the discretized Poisson equation for pressure correction with an unsymmetric matrix is treated. Some numerical experiments to show the accuracy, robustness and efficiency of the proposed method are presented. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

12.
The time integration method with four-order accuracy, self-starting and implicit for the diffuse chemical reaction kinetics equation or the transient instantaneous temperature filed equation was presented. The examples show that both accuracy and stability are better than Runge-Kutta method with four-order. The coefficients of the equation are stored with sparse matrix pattern, so an algorithm is presented which combines a compact storage scheme with reduced computation cost. The computation of the competitive and consecutive reaction in the rotating packed bed, taken as examples, shows that the method is effective.  相似文献   

13.
It is well known that exact projection methods (EPM) on non‐staggered grids suffer for the presence of non‐solenoidal spurious modes. Hence, a formulation for simulating time‐dependent incompressible flows while allowing the discrete continuity equation to be satisfied up to machine‐accuracy, by using a Finite Volume‐based second‐order accurate projection method on non‐staggered and non‐uniform 3D grids, is illustrated. The procedure exploits the Helmholtz–Hodge decomposition theorem for deriving an additional velocity field that enforces the discrete continuity without altering the vorticity field. This is accomplished by first solving an elliptic equation on a compact stencil that is by performing a standard approximate projection method (APM). In such a way, three sets of divergence‐free normal‐to‐face velocities can be computed. Then, a second elliptic equation for a scalar field is derived by prescribing that its additional discrete gradient ensures the continuity constraint based on the adopted linear interpolation of the velocity. Characteristics of the double projection method (DPM) are illustrated in details and stability and accuracy of the method are addressed. The resulting numerical scheme is then applied to laminar buoyancy‐driven flows and is proved to be stable and efficient. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

14.
非线性动力方程的增维精细积分法   总被引:30,自引:0,他引:30  
对线性定常结构的动力系统提出的精细积分法,能得到在数值上逼近于精确解的结果。但是对于非齐次动力方程却涉及到矩阵求逆的困难,而且通常与时间有关的非齐次项不能进入精细积分的细化过程。采用增维的方法,将非齐次动力方程化为齐次方程,在实施精细积分的过程中不必进行矩阵求逆。这种处理方法对于程序实现和提高数值计算的稳定性十分有利,而且在大型问题中可明显提高计算效率,数值算例显示本文方法是有效的。  相似文献   

15.
An improved precise integration method(IPIM) for solving the differential Riccati equation(DRE) is presented.The solution to the DRE is connected with the exponential of a Hamiltonian matrix,and the precise integration method(PIM) for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix.The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE.Based on the error analysis,the criterion for choosing two parameters of the PIM is given.Three kinds of IPIMs for solving the DRE are proposed.The numerical examples show that the IPIM is stable and gives the machine accuracy solutions.  相似文献   

16.
基于动力分析的摄动有限元法及管线可靠度分析模型 ,建立管线的动态随机有限元方程 ,分析了管线质量随机性和几何特征随机性 ,得出了在Elcentro地震波的激励下 ,管线结构的动力响应及具有非穿透性裂纹的抗震可靠度 ;并与Monte Carlo有限元法进行计算比较 ,表明了方法的正确性和可行性  相似文献   

17.
IntroductionTransientstateanalysishasbeenanactiveresearchareainmanyengineeringproblems.Inpracticalsituation ,likestructuralmechanics,thesystemsbeingstudiedareusuallynonlinearandtime_dependent.Theanalyticalmethodstosolvetheseproblemsareimpossible .Thenum…  相似文献   

18.
王金东  高鹏  陈浩然 《力学季刊》2000,21(3):316-321
应用现有的波动方程求解方法解决工程实际问题尚存在一定的局限性。本文在结构动力方程精细逐步积分的基础上,提出了波动方程初边值问题的精细逐步积分法,并分别给出了不同边界条件下的精细逐步积分格式。此数值方法虽然是显式积分方法,却是无条件稳定的。分别用精细逐步积分法和其它已有的方法对两个算例进行了计算,一个是有解析解的例子,该例验证了此方法的准确性,另一个例子是求解由波动方程及初始条件和边界条件组成的有杆抽油系统预测模型,此例验证了精细逐步积分法的高效性。  相似文献   

19.
A method is developed for performing a local reduction of the governing physics for fluid problems with domains that contain a combination of narrow and non‐narrow regions, and the computational accuracy and performance of the method are measured. In the narrow regions of the domain, where the fluid is assumed to have no inertia and the domain height and curvature are assumed small, lubrication, or Reynolds, theory is used locally to reduce the two‐dimensional Navier–Stokes equations to the one‐dimensional Reynolds equation while retaining a high degree of accuracy in the overall solution. The Reynolds equation is coupled to the governing momentum and mass equations of the non‐narrow region with boundary conditions on the mass and momentum flux. The localized reduction technique, termed ‘stitching,’ is demonstrated on Stokes flow for various geometries of the hydrodynamic journal bearing—a non‐trivial test problem for which a known analytical solution is available. The computational advantage of the coupled Stokes–Reynolds method is illustrated on an industrially applicable fully‐flooded deformable‐roll coating example. The examples in this paper are limited to two‐dimensional Stokes flow, but extension to three‐dimensional and Navier–Stokes flow is possible. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

20.
This paper describes the development of a parallel three‐dimensional unstructured non‐isothermal flow solver for the simulation of the injection molding process. The numerical model accounts for multiphase flow in which the melt and air regions are considered to be a continuous incompressible fluid with distinct physical properties. This aspect avoids the complex reconstruction of the interface. A collocated finite volume method is employed, which can switch between first‐ and second‐order accuracy in both space and time. The pressure implicit with splitting of operators algorithm is used to compute the transient flow variables and couple velocity and pressure. The temperature equation is solved using a transport equation with convection and diffusion terms. An upwind differencing scheme is used for the discretization of the convection term to enforce a bounded solution. In order to capture the sharp interface, a bounded compressive high‐resolution scheme is employed. Parallelization of the code is achieved using the PETSc framework and a single program multiple data message passing model. Predicted numerical solutions for several example problems are considered. The first case validates the solution algorithm for moderate Reynolds number flows using a structured mesh. The second case employs an unstructured hybrid mesh showing the capability of the solver to describe highly viscous flows closer to realistic injection molding conditions. The final case presents the non‐isothermal filling of a thick cavity using three mesh sizes and up to 80 processors to assess parallel performance. The proposed algorithm is shown to have good accuracy and scalability. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号