Modeling of crack propagation via an automatic adaptive mesh refinement based on modified superconvergent patch recovery technique

In this paper, an automated adaptive remeshing procedure is presented for simulation of arbitrary shape crack growth in a 2D finite element mesh. The Zienkiewicz-Zhu error estimator is employed in conjunction with a modified SPR technique based on the recovery of gradients using analytical crack-tip fields in order to obtain more accurate estimation of errors. The optimization of crack-tip singular finite element size is achieved through the adaptive mesh strategy. Finally, several numerical examples are illustrated to demonstrate the effectiveness, robustness and accuracy of computational algorithm in calculation of fracture parameters and prediction of crack path pattern.     

On a consistent field transfer in non linear inelastic analysis and ultimate load computation

In this work we propose a field transfer operator for remeshing carried out in the course of incremental analysis of a non linear inelastic behavior. The proposed procedure is geared towards the ultimate load computation of a complex structure, where we choose the appropriate mesh grading for each different phase of computations, starting with a coarse mesh for the initial linear response and going towards a more refined mesh for highly nonlinear inelastic response. The proposed projection operator is developed on the basis of diffuse approximation method. The key feature of such an operator is to guarantee the conservation of relevant mechanics quantities which ensures a superior performance of the proposed field transfer with respect to the standard remeshing procedure. We present the illustrative results both for an isotropic damage model and standard plasticity model, indicating very satisfying performance.     

On the need for the use of error‐controlled finite element analyses in structural shape optimization processes

This work analyzes the influence of the discretization error associated with the finite element (FE) analyses of each design configuration proposed by the structural shape optimization algorithms over the behavior of the algorithm. The paper clearly shows that if FE analyses are not accurate enough, the final solution provided by the optimization algorithm will neither be optimal nor satisfy the constraints. The need for the use of adaptive FE analysis techniques in shape optimum design will be shown. The paper proposes the combination of two strategies to reduce the computational cost related to the use of mesh adaptivity in evolutionary optimization algorithms: (a) the use of an algorithm for the mesh generation by projection of the discretization error, which reduces the computational cost associated with the adaptive FE analysis of each geometrical configuration and (b) the successive increase of the required accuracy of the FE analyses in order to obtain a considerable reduction of the computational cost in the early stages of the optimization process. Copyright © 2011 John Wiley & Sons, Ltd.     

Remeshing for metal forming simulations—Part I: Two‐dimensional quadrilateral remeshing

In this paper, a general framework of practical two‐dimensional quadrilateral remeshing, which includes the determination of remeshing time, automatic quadrilateral mesh generation, and data transfer process, will be formulated. In particular, the current work contains new algorithms of mesh density specification according to the distribution of effective strain‐rate gradients, mesh density smoothing by fast Fourier transform (FFT) and low‐pass filtering techniques, coarsening it by node placement scheme, and a modified Laplacian mesh smoothing technique. The efficiency of the developed remeshing scheme was tested through three practical two‐dimensional metal forming simulations. The results clearly indicate that the algorithms proposed in this study make it possible to simulate two‐dimensional metal forming problems efficiently and automatically. Copyright © 2002 John Wiley & Sons, Ltd.     

A 2D finite element analysis of the effect of numerical parameters on the reliability of Ti6Al4V machining modeling

Abstract

The numerical analysis, based on the finite element modeling (FEM), presents nowadays an efficient computational tool. It allows a better understanding of several thermo-mechanical phenomena involved during the machining process. However, its reliability heavily depends on the accurate definition of the numerical model. In this regard, a FE analysis focused on the 2D modeling of the Ti6Al4V dry orthogonal machining was carried out in this study. The relevance of different numerical meshing approaches and finite elements topologies was studied. The effect of the friction coefficient on the numerical chip morphology, its geometry, the cutting and the feed forces was investigated. The adequacy of several compared adaptive meshing approaches, in terms of the modeling of severe contact conditions taking place around the cutting-edge radius, was underlined in the current study. However, numerical serrated chips, closer to the experimental ones, were only predicted when the pure Lagrangian formulation was adopted and a proper determination of the failure energy was carried out. The definition of different mesh topologies highlighted the efficiency of the 4-node quadrangular mesh, with a suitable edge length, in increasing the agreement with the experimental data, while reducing the computing times.     

Hierarchical remeshing strategies with mesh mapping for topology optimisation

This work investigates the use of hierarchical mesh decomposition strategies for topology optimisation using bi‐directional evolutionary structural optimisation algorithm. The proposed method uses a dual mesh system that decouples the design variables from the finite element analysis mesh. The investigation focuses on previously unexplored areas of these techniques to investigate the effect of five meshing parameters on the analysis solving time (i.e. computational effort) and the analysis quality (i.e. solution optimality). The foreground mesh parameters, including adjacency ratio and minimum and maximum element size, were varied independently across solid and void domain regions. Within the topology optimisation, strategies for controlling the mesh parameters were investigated. The differing effects of these parameters on the efficiency and efficacy of the analysis and optimisation stages are discussed, and recommendations are made for parameter combinations. Some of the key findings were that increasing the adjacency ratio increased the efficiency only modestly – the largest effect was for the minimum and maximum element size parameters – and that the most dramatic reduction in solve time can be achieved by not setting the minimum element size too low, assuming mapping onto a background mesh with a minimum element size of 1. © 2016 The Authors. International Journal for Numerical Methods in Engineering Published by John Wiley & Sons, Ltd.     

比例边界有限元模拟裂纹和夹杂动力相互作用

基于三角形背景网格,任意结构可用 （ ）边多边形比例边界有限元 （Polygon Scaled Boundary Finite Elements, PSBFE）自动离散。相对以往基于比例边界有限元（SBFEM）的应用,该多边形单元不但继承SBFEM半解析求解裂纹尖端奇异性的特性,而且在模拟复杂结构的网格生成和裂纹扩展上具有更高的通用性。本文首次用该单元模拟了动荷载下复合材料裂纹和夹杂相互作用。动荷载稳定裂纹情况下,PSBFE计算结果同现有文献吻合良好,在此基础上,本文结合基于拓扑的局部网格重剖分方法,模拟了动荷载下夹杂和扩展裂纹相互作用。结果表明,硬质夹杂和软质夹杂对结构的动力应力强度因子分别起到抑制和放大的作用。夹杂尺寸,夹杂大小也会在一定范围内影响动力应力强度因子,尺寸越大距离裂纹越近的夹杂影响越大。     

A frontal Delaunay quad mesh generator using the L ∞ norm

In a recent work, a new indirect method to generate all‐quad meshes has been developed. It takes advantage of a well‐known algorithm of the graph theory, namely the Blossom algorithm, which computes in polynomial time the minimum cost perfect matching in a graph. In this paper, we describe a method that allows to build triangular meshes that are better suited for recombination into quadrangles. This is performed by using the infinity norm to compute distances in the meshing process. The alignment of the elements in the frontal Delaunay procedure is controlled by a cross field defined on the domain. Meshes constructed this way have their points aligned with the cross‐field directions, and their triangles are almost right everywhere. Then, recombination with the Blossom‐based approach yields quadrilateral meshes of excellent quality. Copyright © 2013 John Wiley & Sons, Ltd.     

采用B样条曲线进行边界重分划

在金属成形数值模拟过程中 ,坯料网格边界不断变化 ,当使用网格向前推进技术产生网格时必须对旧的边界进行重新分划。本文给出了B样条曲线拟合边界的办法 ,实现了对边界的自动重新分划。数值例子验证此法简单可靠     

塑性成形有限元模拟六面体网格划分和再划分技术

总结了使用三维六面体单元有限元技术模拟塑性体积成形时的网格自动划分技术和自适应网格再划分技术,并分析了各种技术的优缺点。