单点子域积分与差分

通过稳定性分析、显式与隐式积分，表明了单点子域积分相对于差分法的优越性．     

结构非线性动力分析混合时间积分并行算法

综合隐式和显式时间积分技术,对结构非线性动力反应分析提出一种并行混合时间积分算法.该算法采用区域分解技术.将并发性引入到算法中,即利用显式时间积分技术进行界面节点积分而利用隐式算法求解局部子区境.为实现并行混合时间积分算法,设计了灵活的并行数据信息流.编写了该算法的程序,在工作站机群实现了数值算例,验证了算法的精度和性能.计算结果表明该算法具有良好的并行性能,优于隐式算法.     

????????????????????????????

合隐式和显式时间积分技术，对结构非线性动力反应分析提出一种并行混合时间积分算 法. 该算法采用区域分解技术. 将并发性引入到算法中，即利用显式时间积分技术进行界面 节点积分而利用隐式算法求解局部子区域. 为实现并行混合时间积分算法，设计了灵活的并 行数据信息流. 编写了该算法的程序，在工作站机群实现了数值算例，验证了算法的精度和 性能. 计算结果表明该算法具有良好的并行性能，优于隐式算法.     

教学杂记一则

用D-M模型来推导J积分与COD的关系式J=σ,δ时,不少书用图1的ABC表示积分回路,并由式(1),或式(2),或式(3): ...     

弹性薄板弯曲理论中充要的间接变量边界积分方程

本文用非解析开拓法严格地导出了任意区域内双调和函数的一个边界积分表示式,在这个表示式中有两个边界函数和四个常数,它们与双调和函数之间存在一一对应的关系。因此依据这个表示式建立的间接变量边界积分方程与原微分方程边值问题等价。     

多圆荷载下刚性道面变形和应力的计算

本文建立多圆荷载作用下弹性半空间体上薄板的挠度与应力的计算式。荷载数量及分布任意,每个圆荷载密度与轮迹半径彼此相异。对计算式中的反常积分及级数的收敛性予以证明。对含振荡函数反常积分建议一种方便的算法。     

一类奇异积分计算方法及其在断裂力学中的应用

一类奇异积分计算方法及其在断裂力学中的应用任传波，云大真（大连理工大学，大连１１６０２３）在力学及其它的工程计算中，常常遇到奇异积分，不同奇异程度的积分可以采用不同的方法来解决。本文提出的方法可以解决如下一类的奇异积分问题其中１求解方法对于式（１）的...     

关于边界元方法中柯西主值积分的探讨

本文对边界元方法中的各类积分根据其奇异性作分类，并对主值积分的收敛条件、变量替换等进行了讨论，又给出了变替换附加项显式。文中提供的主值积分配项消奇术在边界元方法中是有普遍意义的。     

边界元法中计算几乎奇异积分的一种无奇异算法

边界元法中存在几乎奇异积分的计算困难。引起边界单元上几乎奇异积分的因素是源点到其邻近单元的最小距离δ。本文拓展文[1]的思想，进一步采用分部积分将δ移出奇异积分式中积分核之外，转换后积分核是δ的正则函数。所以几乎强奇异和超奇异积分被化为无奇异的规则积分与解析积分的和，可由通常的Gauss数值积分解。文中应用此正则化技术求解了弹性力学平面问题的近边界点位移和应力。     

基于Mindlin解的桩基础最终沉降计算的应力面积法

利用Ｍｉｎｄｌｉｎ竖向附加应力公式，通过积分得到地基内矩形面积上三角形分布荷载作用下角点下竖向附加应力解析式，并通过对地基内矩形面积上均布和三角形分布载荷作用下角点下竖向附加应力公式关于深度进行积分，得到了计算角点下竖附加应力面积的解析式，根据解析式制表格，可供运用应力面积法进行群桩实体基础等的最终沉降计算时查用。     

RETHINKING TO FINITE DIFFERENCE TIME-STEP INTEGRATIONS

ＲＥＴＨＩＮＫＩＮＧＴＯＦＩＮＩＴＥＤＩＦＦＥＲＥＮＣＥＴＩＭＥ－ＳＴＥＰＩＮＴＥＧＲＡＴＩＯＮＳＺｈｏｎｇＷａｎｘｉｅ（钟万勰）（ＲｅｓｅａｒｔｃｈＩｎｓｔｉｔｕｔｅｏｆＥｎｇｉｎｅｅｒｉｎｇＭｅｃｈａｎｉｃｓ，ＤａｌｉａｎＵｎｉｖｅｒｓｉｔｙｏｆ...     

Continued Fraction Expansion Approaches to Discretizing Fractional Order Derivatives—an Expository Review

This paper attempts to present an expository review of continued fraction expansion (CFE) based discretization schemes for fractional order differentiators defined in continuous time domain. The schemes reviewed are limited to infinite impulse response (IIR) type generating functions of first and second orders, although high-order IIR type generating functions are possible. For the first-order IIR case, the widely used Tustin operator and Al-Alaoui operator are considered. For the second order IIR case, the generating function is obtained by the stable inversion of the weighted sum of Simpson integration formula and the trapezoidal integration formula, which includes many previous discretization schemes as special cases. Numerical examples and sample codes are included for illustrations.     

Stable and accurate schemes for smoothed dissipative particle dynamics

Smoothed dissipative particle dynamics (SDPD) is a mesoscopic particle method that allows to select the level of resolution at which a fluid is simulated. The numerical integration of its equations of motion still suffers from the lack of numerical schemes satisfying all the desired properties such as energy conservation and stability. Similarities between SDPD and dissipative particle dynamics with energy (DPDE) conservation, which is another coarse-grained model, enable adaptation of recent numerical schemes developed for DPDE to the SDPD setting. In this article, a Metropolis step in the integration of the fluctuation/dissipation part of SDPD is introduced to improve its stability.     

ON THE ARBITRARY DIFFERENCE PRECISE INTEGRATION METHOD AND ITS NUMERICAL STABILITY

IntroductionManyproblemsencounteredinengineeringpracticeandotherdisciplinescanbesummarizedintoPDEssuchasosmosis,diffusion,heatconduction,wavepropagation,etc.ItisthenofvitalsignificancehowtosolvePDEsbothrapidlyandefficiently,ThenumericalsolutionsofPDEsarecustomarilyobtainedbythefiniteelementmethod(FEM),thefinitedifferencemethod(FDM)!and.the,,[l'2).Thesemethods,however,showtheirdemeritsforlargercomputationaldomains.AsforFEM,thevastnumberofunknownscausedbyspacecoordinatediscretizationlead…     

High‐order ENO and WENO schemes for unstructured grids

This work describes the implementation and analysis of high‐order accurate schemes applied to high‐speed flows on unstructured grids. The class of essentially non‐oscillatory schemes (ENO), that includes weighted ENO schemes (WENO), is discussed in the paper with regard to the implementation of third‐ and fourth‐order accurate methods. The entire reconstruction process of ENO and WENO schemes is described with emphasis on the stencil selection algorithms. The stencils can be composed by control volumes with any number of edges, e.g. triangles, quadrilaterals and hybrid meshes. In the paper, ENO and WENO schemes are implemented for the solution of the dimensionless, 2‐D Euler equations in a cell centred finite volume context. High‐order flux integration is achieved using Gaussian quadratures. An approximate Riemann solver is used to evaluate the fluxes on the interfaces of the control volumes and a TVD Runge–Kutta scheme provides the time integration of the equations. Such a coupling of all these numerical tools, together with the high‐order interpolation of primitive variables provided by ENO and WENO schemes, leads to the desired order of accuracy expected in the solutions. An adaptive mesh refinement technique provides better resolution in regions with strong flowfield gradients. Results for high‐speed flow simulations are presented with the objective of assessing the implemented capability. Copyright © 2007 John Wiley & Sons, Ltd.     

Higher-order accurate implicit time integration schemes for transport problems

The present paper is concerned with the numerical solution of transient transport problems by means of spatial and temporal discretization methods. The generalized initial boundary value problem of various nonlinear transport phenomena like heat transfer or mass transport is discretized in space by p-finite elements. After finite element discretization, the resulting first-order semidiscrete balance has to be solved with respect to time. Next to the classical generalized-α integration method predicated on the Newmark approach and the evaluation at a generalized midpoint also implicit Runge–Kutta time integration schemes, are presented. Both families of finite difference-based integration schemes are derived for general first-order problems. In contrast to the above-mentioned algorithms, temporal discontinuous and continuous Galerkin methods evaluate the balance equation not at a selected time instant within the timestep, but in an integral sense over the whole time step interval. Therefore, the underlying semidiscrete balance and the continuity of the primary variables are weakly formulated within time steps and between time steps, respectively. Continuous Galerkin methods are obtained by the strong enforcement of the continuity condition as special cases. The introduction of a natural time coordinate allows for the application of standard higher-order temporal shape functions of the p-Lagrange type and the well-known Gau?–Legendre quadrature of associated time integrals. It is shown that arbitrary order accurate integration schemes can be developed within the framework of the proposed temporal p-Galerkin methods. Selected benchmark analyses of calcium diffusion demonstrate the properties of all three methods with respect to non-smooth initial or boundary conditions. Furthermore, the robustness of the present time integration schemes is also demonstrated for the highly nonlinear reaction–diffusion problem of calcium leaching, including the pronounced changes of the reaction term and non-smooth changes of Dirichlet boundary conditions of calcium dissolution.     

结构动力分析中时间积分方法进展

叙述了结构动力分析中时间积分方法的最新发展情况，对这一领域的基本原理和思想进行了总结，重点介绍一些新型计算方法的基本性质，为时间积分方法的进一步研究奠定基础。     

Non‐oscillatory relaxation methods for the shallow‐water equations in one and two space dimensions

In this paper, a new family of high‐order relaxation methods is constructed. These methods combine general higher‐order reconstruction for spatial discretization and higher order implicit‐explicit schemes or TVD Runge–Kutta schemes for time integration of relaxing systems. The new methods retain all the attractive features of classical relaxation schemes such as neither Riemann solvers nor characteristic decomposition are needed. Numerical experiments with the shallow‐water equations in both one and two space dimensions on flat and non‐flat topography demonstrate the high resolution and the ability of our relaxation schemes to better resolve the solution in the presence of shocks and dry areas without using either Riemann solvers or front tracking techniques. Copyright © 2004 John Wiley & Sons, Ltd.