三维Delaunay剖分的断层直接插入算法

研究了在三维空间中进行的Delaunay四面体剖分。在讨论了四面体网格与插入的受限平面数据的各种相交构形的基础上,提出了一种断层直接插入的算法。该算法可以应用于三维数据点集的受限四面体剖分,也可以应用于不同的数据场网格之间的相交运算等问题。     

纹理贴图模型的彩色截面轮廓生成算法研究

针对彩色3D打印彩色轮廓生成问题,本文提出一种基于纹理贴图模型的彩色轮廓生成算法。根据平面投影法建立纹理和三维模型的映射关系,接着对该模型进行切片处理取得轮廓顶点的几何信息和纹理信息;并根据轮廓顶点和纹理图片的映射关系,插值出顶点间色彩信息;同时给出一种内外轮廓的快速识别算法,并提出一种基于偏置算法的彩色轮廓偏置方法,形成颜色厚度,消除阶梯效应。     

Adaptive mesh generation for mould filling problems in resin transfer moulding

In injection moulding processes such as Resin Transfer Moulding (RTM) for example, numerical simulations are usually performed on a fixed mesh, on which the numerical algorithm predict the displacement of the flow front. Error estimations can be used in the numerical algorithm to optimise the mesh for the finite element analysis. The mesh can be also adapted during mould filling to follow the shape of the moving boundary. However, in order to minimize computer time, it is preferable to optimise the mesh before carrying out the filling calculation. In this paper, these ideas are adapted to 3D shells, which represent the most common type of composite parts manufactured by RTM. An error estimator generally used in planar or solid geometries is extended for curved 3D surfaces in the specific case of RTM calculations. The extension consists of a projection of the solution field in the tangent plane to avoid problems related to the curvature of the part. Some other issues specific to shell geometries are pointed out and the results of a filling simulation made on a real part are presented. Non-isothermal filling simulations are also carried out in a rectangular mould to illustrate the stability conditions that arise from the convective heat transfer problem. Finally, an analytical study of radial injections is carried out to illustrate issues related to four types of different mesh refinement procedures: (1) a constant time step, (2) constant radial density (to allow a regular progression of the flow front at each time step), (3) a constant Courant number (to ensure stable thermal simulations); and (4) finally, a constant interpolation error.     

A selective diffraction order based lens plane grating shearing interferometer for the study of bent plates

A simple, elegant and flexible technique for the determination of partial slope, curvature and twist contours of bent plates is proposed. A single frequency cross-grating (SFCG) placed just ahead of the field lens collects the light reflected by the smooth plate model. From the resultant diffraction pattern formed at the Fourier transform (FT) plane, based on the desired contours, only a few orders are permitted selectively to pass through to the recording plane. This manner of real time in-situ filtering of the diffraction orders is new. While slope contours require double exposures, curvature and twist contours are obtained from single exposures.     

Efficient ALE mesh management for 3D quasi‐Eulerian problems

In computational solid mechanics, the ALE formalism can be very useful to reduce the size of finite element models of continuous forming operations such as roll forming. The mesh of these ALE models is said to be quasi‐Eulerian because the nodes remain almost fixed—or almost Eulerian—in the main process direction, although they are required to move in the orthogonal plane in order to follow the lateral displacements of the solid. This paper extensively presents a complete node relocation procedure dedicated to such ALE models. The discussion focusses on quadrangular and hexahedral meshes with local refinements. The main concern of this work is the preservation of the geometrical features and the shape of the free boundaries of the mesh. With this aim in view, each type of nodes (corner, edge, surface and volume) is treated sequentially with dedicated algorithms. A special care is given to highly curved 3D surfaces for which a CPU‐efficient smoothing technique is proposed. This new method relies on a spline surface reconstruction, on a very fast weighted Laplacian smoother with original weights and on a robust reprojection algorithm. The overall consistency of this mesh management procedure is finally demonstrated in two numerical applications. The first one is a 2D ALE simulation of a drawbead, which provides similar results to an equivalent Lagrangian model yet is much faster. The second application is a 3D industrial ALE model of a 16‐stand roll forming line. In this case, all attempts to perform the same simulation by using the Lagrangian formalism have been unsuccessful. Copyright © 2012 John Wiley & Sons, Ltd.     

Automatic quadrilateral/triangular free-form mesh generation for planar regions

Automation of finite element mesh generation holds great benefits for mechanical product development and analysis. In addition to freeing engineers from mundane tasks, automation of mesh generation reduces product cycle design and eliminates human-related errors. Most of the existing mesh generation methods are either semi-automatic or require specific topological information. A fully automatic free-form mesh generation method is described in this paper to alleviate some of these problems. The method is capable of meshing singly or multiply connected convex/concave planar regions. These regions can be viewed as crosssectional areas of 2 1/2 D objects analysed as plane stress, plane strain or axisymmetric stress problems. In addition to being fully automatic, the method produces quadrilateral or triangular elements with aspect rations near one. Moreover, it does not require any topological constraints on the regions to be meshed; i.e. it provides free-form mesh generation. The input to the method includes the region's boundary curves, the element size and the mesh grading information. The method begins by decomposing the planar region to be meshed into convex subregions. Each subregion is meshed by first generating nodes on its boundaries using the input element size. The boundary nodes are then offset to mesh the subregion. The resulting meshes are merged together to form the final mesh. The paper describes the method in detail, algorithms developed to implement it and sample numerical examples. Results on parametric studies of the method performance are also discussed.     

A physically and geometrically nonlinear scaled‐boundary‐based finite element formulation for fracture in elastomers

This paper is devoted to the formulation of a plane scaled boundary finite element with initially constant thickness for physically and geometrically nonlinear material behavior. Special two‐dimensional element shape functions are derived by using the analytical displacement solution of the standard scaled boundary finite element method, which is originally based on linear material behavior and small strains. These 2D shape functions can be constructed for an arbitrary number of element nodes and allow to capture singularities (e.g., at a plane crack tip) analytically, without extensive mesh refinement. Mapping these proposed 2D shape functions to the 3D case, a formulation that is compatible with standard finite elements is obtained. The resulting physically and geometrically nonlinear scaled boundary finite element formulation is implemented into the framework of the finite element method for bounded plane domains with and without geometrical singularities. The numerical realization is shown in detail. To represent the physically and geometrically nonlinear material and structural behavior of elastomer specimens, the extended tube model and the Yeoh model are used. Numerical studies on the convergence behavior and comparisons with standard Q1P0 finite elements demonstrate the correct implementation and the advantages of the developed scaled boundary finite element. Copyright © 2014 John Wiley & Sons, Ltd.     

An integration procedure based on multi‐surface plasticity for plane stress plasticity: Theoretical formulation and application

An integration procedure designed to satisfy plane stress conditions for any constitutive law initially described in 3D and based on classical plasticity theory is presented herein. This method relies on multi‐surface plasticity, which allows associating in series various mechanisms. Three mechanisms have ultimately been used and added to the first one to satisfy the plane stress conditions. They are chosen to generate a plastic flow in the 3 out‐of‐plane directions, whose stresses must be canceled (σ₃₃,σ₁₃, and σ₂₃). The advantage of this method lies in its ease of use for every plastic constitutive law (in the general case of the non‐associated flow rule and with both nonlinear kinematic and isotropic hardening). Method implementation using a cutting plane algorithm is presented in its general framework and then illustrated by the example of a J2‐plasticity material model considering linear kinematic and isotropic hardening. The approach is compared with the same J2‐plasticity model that has been directly derived from a projection of its equations onto the plane stress subspace. The performance of the multi‐surface plasticity method is shown through the comparison of iso‐error and iso‐step contours in both formulations, and lastly with a case study considering a hollow plate subjected to tension. Copyright © 2014 John Wiley & Sons, Ltd.     

A structural mapping technique for geometric parametrization in the optimization of magnetic devices

The continuity and differentiability of object functions is a basic prerequisite for the application of gradient methods in optimization. However, for parameters defining the shape of an electromagnetic device, the finite element discretization in the field analysis introduces discontinuities into the object function which slow down the convergence rate. Additionally, depending on the geometric parametrizaiion employed, the optimization frequently yields shape contours that are impracticable for manufacturing purposes. This paper investigates the problems inherent in geometric parametrization and shows that the discontinuities in the object function are caused by changes in mesh topology as the geometric parameters vary; these changes inevitably follow from the use of free meshing algorithms. As a solution to these shortcomings a structural mapping technique is outlined that maps surface displacements onto the parameters of the finite element mesh and obtains the parameter dependent geometric variations without a change in mesh topology. This resulting geometric parametrization yields continuous object functions without artificial local minima and results in smooth surface contours of the optimized device. Using this new parametrization technique, design sensitivity analysis, is shown to be a reliable and essential part in the efficient application of gradient methods for shape optimization.     

A 2.5D finite/infinite element approach for modelling visco‐elastic bodies subjected to moving loads

The objective of this study is to propose a 2.5D finite/infinite element procedure for dealing with the ground vibrations induced by moving loads. Besides the two in‐plane degrees of freedom (DOFs) per node conventionally used for plane strain elements, an extra DOF is introduced to account for the out‐of‐plane wave transmission. The profile of the half‐space is divided into a near field and a semi‐infinite far field. The near field containing loads and irregular structures is simulated by the finite elements, while the far field covering the soils extending to infinity by the infinite elements with due account taken of the radiation effects for moving loads. Enhanced by the automated mesh expansion procedure proposed previously by the writers, the far field impedances for all the lower frequencies are generated repetitively from the mesh created for the highest frequency considered. Finally, the accuracy of the proposed method is verified through comparison with a number of analytical solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

An effective 3D meshing approach for fractured rocks

An effective mesh generation algorithm is proposed to construct mesh representations for arbitrary fractures in 3D rock masses. With the development of advanced imaging techniques, fractures in a rock mass can be clearly captured by a high‐resolution 3D digital image but with a huge data set. To reduce the data size, corresponding mesh substitutes are required in both visualization and numerical analysis. Fractures in rocks are naturally complicated. They may meet at arbitrary angles at junctions, which could derive topological defects, geometric errors or local connectivity flaws on mesh models. A junction weight is proposed and applied to distinguish fracture junctions from surfaces by an adequate threshold. We take account of fracture junctions and generate an initial surface mesh by a simplified centroidal Voronoi diagram. To further repair the initial mesh, an innovative umbrella operation is designed and adopted to correct mesh topology structures and preserve junction geometry features. Constrained with the aforementioned surface mesh of fracture, a tetrahedral mesh is generated and substituted for the 3D image model to be involved in future numerical analysis. Finally, we take two fractured rock samples as application examples to demonstrate the usefulness and capability of the proposed meshing approach. Copyright © 2015 John Wiley & Sons, Ltd.     

A numerical study of plasticity induced crack closure under plane strain conditions

The level of plasticity induced crack closure (PICC) is greatly affected by stress state. Under plane strain conditions, however, the level and even the existence of PICC still are controversial. The objective here is to study the influence of the main numerical parameters on plane strain PICC, namely the total crack propagation, the number of load cycles between crack increments, the finite element mesh and the parameter used to quantify PICC. The PICC predictions were included in a parallel numerical study of crack propagation, in order to quantify the impact of plane strain values on fatigue life. The results indicate that literature may be overestimating plane strain PICC due to incorrect numerical parameters. The number of load cycles usually considered is unrealistically small, and its increase was found to vanish crack closure, particularly for kinematic hardening. This effect was linked to the ratcheting effect observed at the crack tip. The total crack increment, Δa, must be large enough to obtain stabilized PICC values, but this may imply a huge numerical effort particularly for 3D models. The size of crack tip plastic zone may be overestimated in literature, which means that the meshes used may be too large. Additionally, the crack propagation study showed that the plane strain PICC has usually a dominant effect on fatigue life, and plane stress PICC is only relevant for relatively thin geometries.     

基于点阵投影的粗加工新方法

提出了一种点阵投影式的粗加工方法。该方法对模具或机械产品的造型没有特殊要求,既可以是规则的实体模型,也可以是自由曲面模型。而且允许曲面之间有缝隙或相互重叠;毛坯可以是任意形状,该方法计算模型简单,铣削过程中抬刀很少,加工效率高。     

Three-dimensional vision from a multisensing mechanism

An approach for sensing a three-dimensional (3D) object surface with an arbitrary geometric shape is presented. Combining two different 3D sensing mechanisms, point-array encoding based on affine transformation and fringe encoding based on phase mapping, we construct a mathematic model for 3D vision in which the point-array encoding is initially applied to determine the fringe orders to create a control-vertex mesh with absolute coordinate values in 3D space. Then phase evaluation and phase unwrapping for fringe decoding is performed under the guidance of control-vertex mesh, leading to an absolute phase map and the corresponding range image of the test object. The computer simulations and experimental results are presented to demonstrate the theoretical prediction.     

Curvature contours of flexed plates by a multiple illumination moiré shearing interferometer

A new moiré method is proposed to map partial slope, partial curvature and twist curvature contours of laterally loaded specularly reflecting plate models. Converging HeNe laser light is diffracted by a lens plane cross grating, different combinations of orders are selectively permitted through the Fourier transform plane and collimated to illuminate the plate model. This technique eliminates the difficulties involved in Fourier filtering in the imaging optics, thereby extending the reflection moiré' interferomelry to general plate structures and loading. The slope contours are obtained still by double exposure.     

A new automatic adaptive 3D solid mesh generation scheme for thin‐walled structures

A new algorithm to generate three‐dimensional (3D) mesh for thin‐walled structures is proposed. In the proposed algorithm, the mesh generation procedure is divided into two distinct phases. In the first phase, a surface mesh generator is employed to generate a surface mesh for the mid‐surface of the thin‐walled structure. The surface mesh generator used will control the element size properties of the final mesh along the surface direction. In the second phase, specially designed algorithms are used to convert the surface mesh to a 3D solid mesh by extrusion in the surface normal direction of the surface. The extrusion procedure will control the refinement levels of the final mesh along the surface normal direction. If the input surface mesh is a pure quadrilateral mesh and refinement level in the surface normal direction is uniform along the whole surface, all hex‐meshes will be produced. Otherwise, the final 3D meshes generated will eventually consist of four types of solid elements, namely, tetrahedron, prism, pyramid and hexahedron. The presented algorithm is highly flexible in the sense that, in the first phase, any existing surface mesh generator can be employed while in the second phase, the extrusion procedure can accept either a triangular or a quadrilateral or even a mixed mesh as input and there is virtually no constraint on the grading of the input mesh. In addition, the extrusion procedure development is able to handle structural joints formed by the intersections of different surfaces. Numerical experiments indicate that the present algorithm is applicable to most practical situations and well‐shaped elements are generated. Copyright © 2004 John Wiley & Sons, Ltd.     

Meshing curved wire‐frame models through energy minimization and packing of ellipses

In the initial phase of structural part design, wire‐frame models are sometimes used to represent the shapes of curved surfaces. Finite Element Analysis (FEA) of a curved surface requires a well shaped, graded mesh that smoothly interpolates the wire frame. This paper describes an algorithm that generates such a triangular mesh from a wire‐frame model in the following two steps: (1) construct a triangulated surface by minimizing the strain energy of the thin‐plate‐bending model, and (2) generate a mesh by the bubble meshing method on the projected plane and project it back onto the triangulated surface. Since the mesh elements are distorted by the projection, the algorithm generates an anisotropic mesh on the projected plane so that an isotropic mesh results from the final projection back onto the surface. Extensions of the technique to anisotropic meshing and quadrilateral meshing are also discussed. The algorithm can generate a well‐shaped, well‐graded mesh on a smooth curved surface. Copyright © 1999 John Wiley & Sons, Ltd.     

Three‐dimensional grayscale‐free topology optimization using a level‐set based r‐refinement method

In this paper, we propose a three‐dimensional (3D) grayscale‐free topology optimization method using a conforming mesh to the structural boundary, which is represented by the level‐set method. The conforming mesh is generated in an r‐refinement manner; that is, it is generated by moving the nodes of the Eulerian mesh that maintains the level‐set function. Although the r‐refinement approach for the conforming mesh generation has many benefits from an implementation aspect, it has been considered as a difficult task to stably generate 3D conforming meshes in the r‐refinement manner. To resolve this task, we propose a new level‐set based r‐refinement method. Its main novelty is a procedure for minimizing the number of the collapsed elements whose nodes are moved to the structural boundary in the conforming mesh; in addition, we propose a new procedure for improving the quality of the conforming mesh, which is inspired by Laplacian smoothing. Because of these novelties, the proposed r‐refinement method can generate 3D conforming meshes at a satisfactory level, and 3D grayscale‐free topology optimization is realized. The usefulness of the proposed 3D grayscale‐free topology optimization method is confirmed through several numerical examples. Copyright © 2017 John Wiley & Sons, Ltd.     

Three-dimensional Correction of the Stress Intensity Factor for Plate with a Notch

The 3D effects of the 2D mode I stress intensity factor for the plate with a V-shaped straight through-thickness notch are investigated by the finite element method and three-dimensional thicknessdependent correction of SIF is suggested. The correction relies on the assumed relationship between the SIF and the constraint factor (out-of-plane degree of freedom). The 3D finite element mesh generator combining the 2D in-plane adaptive unstructured mesh with the structured through-thickness mesh Is developed and applied for the analysis purposes. Three-dimensionality was examined by using two independent indicators, namely, stress- or strain-based constraint factors. The three-point bend and tensile center-cracked plates are investigated. The results demonstrated that the developed 3D corrections may be treated as upper bound estimates of the SIF for three-point bend plate,while directly obtained numerical values are considered as lower bound estimates. Analysis of the tensile center-cracked plate demonstrated a different nature of the 3D SIF profile, which cannot be simply explained as a transitional state between plain strain and plain stress. Therefore, the suggested 3D correction concept is of a particular character.     

三角网格模型上的四边形曲线网生成新方法

四边形网格划分是组合曲面建模技术的首要条件。针对海量流形三角网格数据,提出了基于网格简化技术与调和映射算法的四边形网格生成新方法--映射法。该方法采用基于顶点删除的网格简化技术对三角网格模型进行简化,进而借助调和映射算法将简化网格映射到二维平面上进行四边形划分,并将所获得的平面四边形节点数据逆映射回物理域,采用短程线边界形式最终得到适于组合曲面建模的空间四边形拓扑。该方法简单、实用,运行速度较快,实际的算例也验证了方法的有效性与可行性。