SPH方法Delaunay三角刨分与自由液面重构

光滑粒子法(SPH)作为一种拉格朗日型无网格方法,兼具欧拉网格方法和拉格朗日网格方法的优势,已经成功应用于科学和工程的众多领域。SPH方法后处理一般基于无规则分布的粒子,不如网格类方法后处理简便、直接。另外,SPH方法模拟自由液面流动等问题时,通过粒子位置难以重构自由液面的准确位置。发展一种基于Delaunay三角刨分的SPH后处理方法,即先基于SPH粒子位置利用Delaunay三角刨分建立三角网格,然后将粒子信息转化成网格单元/节点信息,从而可以在三角网格上进行后处理,实现基于网格方法的后处理功能,并可以在三角网格上直接提取或重构自由液面。将本文的方法应用到液滴碰撞和溃坝流SPH模拟结果的后处理中,得到了非常好的结果,表明本文的方法有效可靠。     

Automatic generation of unstructured grids with Delaunay triangulation and its application

This paper is consisted of two parts. In the first part, a method is described which generates two-dimensional triangle mesh using the Delaunay triangulation criterion. An automatic algorithm was proposed which combines several advantages of the existing methods. Local mesh refinement can also be easily performed with this method. Examples of generated grids were presented for several convex, non-convex and multi-connected domains to demonstrate the effectiveness and feasibility of the proposed method. In the second part, the turbulent heat transfer in an annular space finned by wave-like longitudinal fins was numerical simulated. The proposed technique was adopted to generate the grid in the cross-section. The standard K-ɛ model in conjuction with wall function method was used to simulate the fluid flow and heat transfer in the complex geometry. The discretization of the governing equations was described. The computational results were compared with the authors' test data and the agreement was reasonably good. Received on 9 July 1998     

A note on incremental equations for a new class of constitutive relations for elastic bodies

Some new classes of constitutive relations for elastic bodies have been proposed in the literature, wherein the stresses and strains are obtained from implicit constitutive relations. A special case of the above relations corresponds to a class of constitutive equations where the linearized strain tensor is given as a nonlinear function of the stresses. For such constitutive equations we consider the problem of decomposing the stresses into two parts: one corresponds to a time-independent solution of the boundary value problem, plus a small (in comparison with the above) time-dependent stress tensor. The effect of this initial time-independent stress in the propagation of a small wave motion is studied for an infinite medium.     

A simple quality triangulation algorithm for complex geometries

This paper presents a simple algorithm for quality triangulation in domains with complex geometries. Based on the fact that the equilateral triangles (regular meshes) are ideal for numerical computations in computational fluids dynamics (CFD) analysis, the proposed algorithm starts with an initial equilateral triangle mesh covering the whole domain. Nodes close to the boundary edges satisfy the so‐called non‐encroaching criterion, the distance from any inserted node to any boundary vertices and the midpoints of any boundary edge is greater than a given characteristic length. Both nearly uniform and non‐uniform triangle meshes can be generated using a mesh size reduction technique. Local refinement is achieved by using transition layers. More regular meshes can be generated in the interior of the domain and all angles of the triangle mesh produced by this algorithm are proven to be bounded in a reasonable range (19.5–141°). Copyright © 2010 John Wiley & Sons, Ltd.     

A new tracking algorithm of PIV and removal of spurious vectors using Delaunay tessellation

 A new algorithm of Delaunay Tessellation Particle Tracking Velocimetry (DT-PTV in abbreviation) is proposed for tracking particles in images of a PIV system by making use of the Delaunay tessellation (DT). The algorithm is tested by using numerically simulated particle images. The calculation results based on DT are compared with those obtained by a conventional algorithm of Binary Image Cross-correlation method (BICC). The new algorithm shows higher performance of obtaining more identical particles in two consecutive images correctly with shorter computation time even if the images contain many particles. A further application of DT to elimination of spurious vectors is also discussed. Received: 24 November 1997/Accepted: 7 August 1998     

A variational formulation for the incremental homogenization of elasto-plastic composites

This work addresses the micro–macro modeling of composites having elasto-plastic constituents. A new model is proposed to compute the effective stress–strain relation along arbitrary loading paths. The proposed model is based on an incremental variational principle (Ortiz, M., Stainier, L., 1999. The variational formulation of viscoplastic constitutive updates. Comput. Methods Appl. Mech. Eng. 171, 419–444) according to which the local stress–strain relation derives from a single incremental potential at each time step. The effective incremental potential of the composite is then estimated based on a linear comparison composite (LCC) with an effective behavior computed using available schemes in linear elasticity. Algorithmic elegance of the time-integration of J₂ elasto-plasticity is exploited in order to define the LCC. In particular, the elastic predictor strain is used explicitly. The method yields a homogenized yield criterion and radial return equation for each phase, as well as a homogenized plastic flow rule. The predictive capabilities of the proposed method are assessed against reference full-field finite element results for several particle-reinforced composites.     

A general method for the incremental solution of elastic-plastic problems

Summary A general iterative method for the incremental solution of elastic-plastic problems is described. The convergence of the procedure on the rigorous solution is demonstrated.

---

Sommario Si espone un procedimento iterativo generale per la soluzione incrementale dei problemi elasto-plastici. Si dimostra la convergenza del procedimento verso la soluzione effettiva.

---

---

The results presented in this paper were obtained in the course of research sponsored by the C.N.R., Gruppo Plasticità     

A multi-field incremental variational framework for gradient-extended standard dissipative solids

The paper presents a constitutive framework for solids with dissipative micro-structures based on compact variational statements. It develops incremental minimization and saddle point principles for a class of gradient-type dissipative materials which incorporate micro-structural fields (micro-displacements, order parameters, or generalized internal variables), whose gradients enter the energy storage and dissipation functions. In contrast to classical local continuum approaches to inelastic solids based on locally evolving internal variables, these global micro-structural fields are governed by additional balance equations including micro-structural boundary conditions. They describe changes of the substructure of the material which evolve relatively to the material as a whole. Typical examples are theories of phase field evolution, gradient damage, or strain gradient plasticity. Such models incorporate non-local effects based on length scales, which reflect properties of the material micro-structure. We outline a unified framework for the broad class of first-order gradient-type standard dissipative solids. Particular emphasis is put on alternative multi-field representations, where both the microstructural variable itself as well as its dual driving force are present. These three-field settings are suitable for models with threshold- or yield-functions formulated in the space of the driving forces. It is shown that the coupled macro- and micro-balances follow in a natural way as the Euler equations of minimization and saddle point principles, which are based on properly defined incremental potentials. These multi-field potential functionals are outlined in both a continuous rate formulation and a time-space-discrete incremental setting. The inherent symmetry of the proposed multi-field formulations is an attractive feature with regard to their numerical implementation. The unified character of the framework is demonstrated by a spectrum of model problems, which covers phase field models and formulations of gradient damage and plasticity.     

基于Delaunay三角化的无网格法计算结果后处理

最近新发展起来的无网格方法由于不需要显式网格，节省了网格生成所需的大量时间，并且避免了网格畸变问题，所以在处理一些特殊问题如移动边界、大变形、高梯度等方面显示出特殊的优越性。但另一方面也使得计算结果的全域后处理遇到困难。提出了一种基于Delaunay三角形背景网格的实用无网格计算结果后处理方法，以无网格离散节点为顶点生成Delaunay三角形，将无网格法计算得到的节点应力值映射插值得到三角形内的应力场云图颜色。给出的二维线弹性应力分析算例表明方法可靠实用。     

A second-moment incremental formulation for the mean-field homogenization of elasto-plastic composites

In this paper, the incremental formulation for the mean-field homogenization (MFH) of elasto-plastic composites is enriched by including second statistical moments of per-phase strain increment fields, thus combining two advantages. The first one is to handle non-monotonic loading histories and the second is to better account for the heterogeneity of microscopic fields. The proposal is currently restricted to elasto-plasticity with J2 flow theory in each phase, under the small perturbation hypothesis. The formulation crucially exploits the return mapping algorithm for the J2 model, with its two steps: elastic predictor, and plastic corrections. It is shown that the second-moment measure of the average von Mises stress in each phase at the elastic predictor step plays a major role in the computation of both the average stress and the comparison tangent operator. The proposal is implemented for an extended Mori-Tanaka scheme. Predictions are compared to results provided by full-field, finite element computations of representative volume elements or unit cells, for various composite materials, with polymer or metal matrices. There are cases where the predictions of the proposed modeling improve significantly over those of a first-order incremental formulation.     

Surface triangulation for polygonal models based on CAD data

This paper presents an approach to the generation of unstructured surface meshes for Computer‐Aided Design (CAD) surface models represented as lists of polygons with minimum user interventions. Stereolithography (STL) data are adopted as an interface between a CAD system and the surface grid generator. STL files often include problems such as overlapping surfaces, gaps, and intersections. They have to be revised quickly and automatically before the surface models are used for the background grid of the surface grid generation. In this paper, we describe an automatic revision method for use with STL‐defined surface models. The advancing front method using geometric features is adopted directly on the modified STL surfaces. The capability of the method is demonstrated for several 3D surface models. Copyright © 2002 John Wiley & Sons, Ltd.     

A stabilized incremental projection scheme for the incompressible Navier–Stokes equations

It is well known that any spatial discretization of the saddle‐point Stokes problem should satisfy the Ladyzhenskaya–Brezzi–Babuska (LBB) stability condition in order to prevent the appearance of spurious pressure modes. Particularly, if an equal‐order approximation is applied, the Schur complement (or, as called some times, the Uzawa matrix) of the pressure system has a non‐trivial null space that gives rise to such modes. An idea in the past was that all the schemes that solve a Poisson equation for the pressure rather than the Uzawa pressure equation (splitting/projection methods) should overcome this difficulty; this idea was wrong. There is numerical evidence that at least the so‐called incremental projection scheme still suffers from spurious pressure oscillations if an equal‐order approximation is applied. The present paper tries to distinguish which projection requires LBB‐compliant approximation and which does not. Moreover, a stabilized version of the incremental projection scheme is derived. Proper bounds for the stabilization parameter are also given. The numerical results show that the stabilized scheme does indeed achieve second‐order accuracy and does not produce spurious (node to node) pressure oscillations. Copyright © 2001 John Wiley & Sons, Ltd.     

拉索穹顶结构非线性分析的混合有限元增量法

拉索穹顶结构是由受压桅杆和拉索组成的新型柔性大跨度空间组合结构,几何上表现为极强的非线性特性,计算困难,本文应用有限元法,结合拉索穹顶结构特征,假定拉索和桅杆的受力满足虎克宣定律,建立了可以直接考虑拉索垂度影响的两节点索单元模型,并与两节点直杆单元相结合,基于修正的拉格朗日描述方法和虚功原理建立了拉索穹顶结构非线性分析的混合有限元增量方程。采用荷载增量法与Ｎｅｗｔｏｎ－Ｒａｐｈｓｏｎ法相结合的求解     

Approximate MPA-based method for performing incremental dynamic analysis

In performance based earthquake engineering, it is important to accurately predict the seismic demand and capacities of structures. One recent estimation method is incremental dynamic analysis (IDA), which requires a series of nonlinear response history analyses (RHA) of the structure under various ground motions, each scaled to multiple levels of intensity, selected to cover entire range of structural response from elasticity, to yield and finally global dynamic instability. The implementation of IDA requires intensive computation and detailed knowledge of the nonlinear RHA of structures. In response to the complexity of IDA, an approximate method based on modal pushover analysis (MPA-based IDA) was developed. In MPA-based IDA, seismic demands are computed using the nonlinear RHA of the equivalent SDF systems instead of using nonlinear RHA of MDF systems. The objective of this study is to develop a simpler MPA-based IDA procedure that can avoid nonlinear RHA of equivalent SDF systems. For this purpose, MPA-based IDA employs the empirical equation of the inelastic displacement ratio (C _R ), defined as the ratio of peak displacement of the inelastic SDF system to that of the corresponding elastic SDF system given the strength ratio R, and that of the collapse strength ratio (R _c), which is the ratio of collapse intensity to yield strength. The proposed procedure is verified by comparing the seismic demands and capacities of 6-, 9-, and 20-story steel moment frames as determined by the proposed method and exact IDA.     

Adaptive delaunay triangulation with multidimensional dissipation scheme for high-speed compressible flow analysis

IntroductionHigh-speed compressible flows normally involve many complex flow phenomena,suchas shock waves,flow expansions,and shock-shock interactions[1].Effects of thesephenomena are critical in the design of high-speed structures.These flows are charact…     

A note on a family of proximal gradient methods for quasi-static incremental problems in elastoplastic analysis

Accelerated proximal gradient methods have recently been developed for solving quasi-static incremental problems of elastoplastic analysis with some different yield criteria. It has been demonstrated through numerical experiments that these methods can outperform conventional optimization-based approaches in computational plasticity. However, in literature these algorithms are described individually for specific yield criteria, and hence there exists no guide for application of the algorithms to other yield criteria. This short paper presents a general form of algorithm design, independent of specific forms of yield criteria, that unifies the existing proximal gradient methods. Clear interpretation is also given to each step of the presented general algorithm so that each update rule is linked to the underlying physical laws in terms of mechanical quantities.