Thin shell analysis from scattered points with maximum‐entropy approximants

We present a method to process embedded smooth manifolds using sets of points alone. This method avoids any global parameterization and hence is applicable to surfaces of any genus. It combines three ingredients: (1) the automatic detection of the local geometric structure of the manifold by statistical learning methods; (2) the local parameterization of the surface using smooth meshfree (here maximum‐entropy) approximants; and (3) patching together the local representations by means of a partition of unity. Mesh‐based methods can deal with surfaces of complex topology, since they rely on the element‐level parameterizations, but cannot handle high‐dimensional manifolds, whereas previous meshfree methods for thin shells consider a global parametric domain, which seriously limits the kinds of surfaces that can be treated. We present the implementation of the method in the context of Kirchhoff–Love shells, but it is applicable to other calculations on manifolds in any dimension. With the smooth approximants, this fourth‐order partial differential equation is treated directly. We show the good performance of the method on the basis of the classical obstacle course. Additional calculations exemplify the flexibility of the proposed approach in treating surfaces of complex topology and geometry. Copyright © 2010 John Wiley & Sons, Ltd.     

A consistent grayscale‐free topology optimization method using the level‐set method and zero‐level boundary tracking mesh

This paper proposes a level‐set based topology optimization method incorporating a boundary tracking mesh generating method and nonlinear programming. Because the boundary tracking mesh is always conformed to the structural boundary, good approximation to the boundary is maintained during optimization; therefore, structural design problems are solved completely without grayscale material. Previously, we introduced the boundary tracking mesh generating method into level‐set based topology optimization and updated the design variables by solving the level‐set equation. In order to adapt our previous method to general structural optimization frameworks, the incorporation of the method with nonlinear programming is investigated in this paper. To successfully incorporate nonlinear programming, the optimization problem is regularized using a double‐well potential. Furthermore, the sensitivities with respect to the design variables are strictly derived to maintain consistency in mathematical programming. We expect the investigation to open up a new class of grayscale‐free topology optimization. The usefulness of the proposed method is demonstrated using several numerical examples targeting two‐dimensional compliant mechanism and metallic waveguide design problems. Copyright © 2014 John Wiley & Sons, Ltd.     

Weight functions for the determination of stress intensity factor and T‐stress for semi‐elliptical cracks in finite thickness plate

This paper presents the application of weight function method for the calculation of stress intensity factors (K) and T‐stress for surface semi‐elliptical crack in finite thickness plates subjected to arbitrary two‐dimensional stress fields. New general mathematical forms of point load weight functions for K and T have been formulated by taking advantage of the knowledge of a few specific weight functions for two‐dimensional planar cracks available in the literature and certain properties of weight function in general. The existence of the generalised forms of the weight functions simplifies the determination of specific weight functions for specific crack configurations. The determination of a specific weight function is reduced to the determination of the parameters of the generalised weight function expression. These unknown parameters can be determined from reference stress intensity factor and T‐stress solutions. This method is used to derive the weight functions for both K and T for semi‐elliptical surface cracks in finite thickness plates, covering a wide range of crack aspect ratio (a/c) and relative depth (a/t) at any point along the crack front. The derived weight functions are then validated against stress intensity factor and T‐stress solutions for several linear and nonlinear two‐dimensional stress distributions. These derived weight functions are particularly useful for the development of two‐parameter fracture and fatigue models for surface cracks subjected to fluctuating nonlinear stress fields, such as these resulting from surface treatment (shot peening), stress concentration or welding (residual stress).     

Optimal transportation meshfree approximation schemes for fluid and plastic flows

We develop an optimal transportation meshfree (OTM) method for simulating general solid and fluid flows, including fluid–structure interaction. The method combines concepts from optimal transportation theory with material‐point sampling and max‐ent meshfree interpolation. The proposed OTM method generalizes the Benamou–Brenier differential formulation of optimal mass transportation problems to problems including arbitrary geometries and constitutive behavior. The OTM method enforces mass transport and essential boundary conditions exactly and is free from tension instabilities. The OTM method exactly conserves linear and angular momentum and its convergence characteristics are verified in standard benchmark problems. We illustrate the range and scope of the method by means of two examples of application: the bouncing of a gas‐filled balloon off a rigid wall; and the classical Taylor‐anvil benchmark test extended to the hypervelocity range. Copyright © 2010 John Wiley & Sons, Ltd.     

Consistent and stable meshfree Galerkin methods using the virtual element decomposition

Over the past two decades, meshfree methods have undergone significant development as a numerical tool to solve partial differential equations (PDEs). In contrast to finite elements, the basis functions in meshfree methods are smooth (nonpolynomial functions), and they do not rely on an underlying mesh structure for their construction. These features render meshfree methods to be particularly appealing for higher‐order PDEs and for large deformation simulations of solid continua. However, a deficiency that still persists in meshfree Galerkin methods is the inaccuracies in numerical integration, which affects the consistency and stability of the method. Several previous contributions have tackled the issue of integration errors with an eye on consistency, but without explicitly ensuring stability. In this paper, we draw on the recently proposed virtual element method, to present a formulation that guarantees both the consistency and stability of the approximate bilinear form. We adopt maximum‐entropy meshfree basis functions, but other meshfree basis functions can also be used within this framework. Numerical results for several two‐dimensional and three‐dimensional elliptic (Poisson and linear elastostatic) boundary‐value problems that demonstrate the effectiveness of the proposed formulation are presented. Copyright © 2017 John Wiley & Sons, Ltd.     

Low‐Dimensional Topological Crystalline Insulators

Topological crystalline insulators (TCIs) are recently discovered topological phase with robust surface states residing on high‐symmetry crystal surfaces. Different from conventional topological insulators (TIs), protection of surface states on TCIs comes from point‐group symmetry instead of time‐reversal symmetry in TIs. The distinct properties of TCIs make them promising candidates for the use in novel spintronics, low‐dissipation quantum computation, tunable pressure sensor, mid‐infrared detector, and thermoelectric conversion. However, similar to the situation in TIs, the surface states are always suppressed by bulk carriers, impeding the exploitation of topology‐induced quantum phenomenon. One effective way to solve this problem is to grow low‐dimensional TCIs which possess large surface‐to‐volume ratio, and thus profoundly increase the carrier contribution from topological surface states. Indeed, through persistent effort, researchers have obtained unique quantum transport phenomenon, originating from topological surface states, based on controllable growth of low‐dimensional TCIs. This article gives a comprehensive review on the recent progress of controllable synthesis and topological surface transport of low‐dimensional TCIs. The possible future direction about low‐dimensional TCIs is also briefly discussed at the end of this paper.     

Elasticity solution for vibration of 2‐D curved beams with variable curvatures using a spectral‐sampling surface method

An accurate spectral‐sampling surface method for the vibration analysis of 2‐D curved beams with variable curvatures and general boundary conditions is presented. The method combines the advantages of the sampling surface method and spectral method. The formulation is based on the 2‐D elasticity theory, which provides complete accuracy and efficiency for curved beams with arbitrary thicknesses and variable curvatures because no other assumptions on the deformations and stresses along the thickness direction are introduced. Specifically, a set of non‐equally spaced sampling surfaces parallel to the beam's middle surface are primarily collocated along the thickness direction, and the displacements of these surfaces are chosen as fundamental beam unknowns. This fact provides an opportunity to derive elasticity solutions for thick beams with a prescribed accuracy by selecting sufficient sampling surfaces. Each of the fundamental beam unknowns is then invariantly expanded as Chebyshev polynomials of the first kind, and the problems are stated in variational form with the aid of the penalty technique and Lagrange multipliers, which provide complete flexibility to describe any arbitrary boundary conditions. Finally, the desired solutions are obtained by the variational operation. Copyright © 2016 John Wiley & Sons, Ltd.     

An element‐free Galerkin method for 3D crack propagation simulation under complicated stress conditions

This study develops an element‐free Galerkin method based on the moving least‐squares approximation to trace three‐dimensional crack propagation under complicated stress conditions. The crack surfaces are modelled by a collection of planar triangles that are added when cracks propagate. The visibility criterion is adopted to treat the screening effect of the cracks on the influenced domain of a Gaussian point. Cracks are assumed to propagate in the perpendicular planes at crack front points when the strain energy release rates reach the material fracture toughness. This method is unique in that it uses a nonlinear contact iterative algorithm to consider contributions of crack surface interaction to the global equilibrium equations, so that crack opening, sliding and closing under complicated stress states can be efficiently modelled. Two numerical examples of three‐dimensional quasi‐static crack propagation were modelled with satisfactory results. Copyright © 2012 John Wiley & Sons, Ltd.     

Three‐dimensional meshfree‐enriched finite element formulation for micromechanical hyperelastic modeling of particulate rubber composites

A three‐dimensional microstructure‐based finite element framework is presented for modeling the mechanical response of rubber composites in the microscopic level. This framework introduces a novel finite element formulation, the meshfree‐enriched FEM, to overcome the volumetric locking and pressure oscillation problems that normally arise in the numerical simulation of rubber composites using conventional displacement‐based FEM. The three‐dimensional meshfree‐enriched FEM is composed of five‐noded tetrahedral elements with a volume‐weighted smoothing of deformation gradient between neighboring elements. The L₂‐orthogonality property of the smoothing operator enables the employed Hu–Washizu–de Veubeke functional to be degenerated to an assumed strain method, which leads to a displacement‐based formulation that is easily incorporated with the periodic boundary conditions imposed on the unit cell. Two numerical examples are analyzed to demonstrate the effectiveness of the proposed approach. Copyright © 2012 John Wiley & Sons, Ltd.     

A meshfree‐enriched finite element method for compressible and near‐incompressible elasticity

In this paper, a two‐dimensional displacement‐based meshfree‐enriched FEM (ME‐FEM) is presented for the linear analysis of compressible and near‐incompressible planar elasticity. The ME‐FEM element is established by injecting a first‐order convex meshfree approximation into a low‐order finite element with an additional node. The convex meshfree approximation is constructed using the generalized meshfree approximation method and it possesses the Kronecker‐delta property on the element boundaries. The gradient matrix of ME‐FEM element satisfies the integration constraint for nodal integration and the resultant ME‐FEM formulation is shown to pass the constant stress test for the compressible media. The ME‐FEM interpolation is an element‐wise meshfree interpolation and is proven to be discrete divergence‐free in the incompressible limit. To prevent possible pressure oscillation in the near‐incompressible problems, an area‐weighted strain smoothing scheme incorporated with the divergence‐free ME‐FEM interpolation is introduced to provide the smoothing on strains and pressure. With this smoothed strain field, the discrete equations are derived based on a modified Hu–Washizu variational principle. Several numerical examples are presented to demonstrate the effectiveness of the proposed method for the compressible and near‐incompressible problems. Copyright © 2012 John Wiley & Sons, Ltd.     

Novel adaptive meshfree integration techniques in meshless methods

We propose two strategies of novel adaptive numerical integration based on mapping techniques for solving the complicated problems of domain integration encountered in meshfree methods. Several mapping methods are presented in detail that map a complex integration domain to much simpler ones, for example, squares, triangles or circles. The techniques described in the paper can be applied to both global and local weak forms, and the highly nonlinear meshfree integrands are evaluated with controlled accuracy. The necessity of the clumsy procedure of background mesh or cell structures used for integration purpose in existing meshfree methods is avoided, and many meshfree methods that require the domain integration can now become ‘truly meshfree’. Various numerical examples in two dimensions are considered to demonstrate the applicability and the effectiveness of the proposed methods and it shows that the accuracy is improved significantly. Their obtained results are compared with analytical solutions and other approaches and very good agreements are found. Additionally, some three‐dimensional cases applied by the present methods are also examined. Copyright © 2012 John Wiley & Sons, Ltd.     

Methods for connecting dissimilar three‐dimensional finite element meshes

Two methods are presented for connecting dissimilar three‐dimensional finite element meshes. The first method combines the concept of master and slave surfaces with the uniform strain approach for finite elements. By modifying the boundaries of elements on a slave surface, corrections are made to element formulations such that first‐order patch tests are passed. The second method is based entirely on constraint equations, but only passes a weaker form of the patch test for non‐planar surfaces. Both methods can be used to connect meshes with different element types. In addition, master and slave surfaces can be designated independently of relative mesh resolutions. Example problems in three‐dimensional linear elasticity are presented. Published in 2000 by John Wiley & Sons, Ltd.     

hp‐Generalized FEM and crack surface representation for non‐planar 3‐D cracks

A high‐order generalized finite element method (GFEM) for non‐planar three‐dimensional crack surfaces is presented. Discontinuous p‐hierarchical enrichment functions are applied to strongly graded tetrahedral meshes automatically created around crack fronts. The GFEM is able to model a crack arbitrarily located within a finite element (FE) mesh and thus the proposed method allows fully automated fracture analysis using an existing FE discretization without cracks. We also propose a crack surface representation that is independent of the underlying GFEM discretization and controlled only by the physics of the problem. The representation preserves continuity of the crack surface while being able to represent non‐planar, non‐smooth, crack surfaces inside of elements of any size. The proposed representation also provides support for the implementation of accurate, robust, and computationally efficient numerical integration of the weak form over elements cut by the crack surface. Numerical simulations using the proposed GFEM show high convergence rates of extracted stress intensity factors along non‐planar curved crack fronts and the robustness of the method. Copyright © 2008 John Wiley & Sons, Ltd.     

One‐dimensional analysis of thin‐walled closed beams having general cross‐sections

This paper presents a new one‐dimensional theory for static and dynamic analysis of thin‐walled closed beams with general cross‐sections. Existing one‐dimensional approaches are useful only for beams with special cross‐sections. Coupled deformations of torsion, warping and distortion are considered in the present work and a new approach to determine sectional warping and distortion shapes is proposed. One‐dimensional C⁰ beam elements based on the present theory are employed for numerical analysis. The effectiveness of the present theory is demonstrated in the analysis of thin‐walled beams having pillar sections of automobiles and excavators. Copyright © 2000 John Wiley & Sons, Ltd.     

Meshfree co‐rotational formulation for two‐dimensional continua

In this paper, a meshfree co‐rotational formulation for two‐dimensional continua is proposed. In a co‐rotational formulation, the motion of a body is separated into rigid motion and strain‐producing deformation. Traditionally, this has been done in the setting of finite elements for beams and shell‐type elements. In the present work every node in a meshfree discretized domain has its own co‐rotating coordinate system. Three key ingredients are established in order to apply the co‐rotational formulation: (i) the relationship between global and local variables, (ii) the angle of rotation of a typical co‐rotating coordinate system, and (iii) a variationally consistent tangent stiffness matrix. An algorithm for the co‐rotational formulation based on load control is provided. Maximum‐entropy basis functions are used to discretize the domain and stabilized nodal integration is implemented to construct the global system of equations. Numerical examples are presented to demonstrate the validity of the meshfree co‐rotational formulation. Copyright © 2009 John Wiley & Sons, Ltd.     

A three‐dimensional hybrid meshfree‐Cartesian scheme for fluid–body interaction

A numerical method based on a hybrid meshfree‐Cartesian grid is developed for solving three‐dimensional fluid–solid interaction (FSI) problems involving solid bodies undergoing large motion. The body is discretized and enveloped by a cloud of meshfree nodes. The motion of the body is tracked by convecting the meshfree nodes against a background of Cartesian grid points. Spatial discretization of second‐order accuracy is accomplished by the combination of a generalized finite difference (GFD) method and conventional finite difference (FD) method, which are applied to the meshfree and Cartesian nodes, respectively. Error minimization in GFD is carried out by singular value decomposition. The discretized equations are integrated in time via a second‐order fractional step projection method. A time‐implicit iterative procedure is employed to compute the new/evolving position of the immersed bodies together with the dynamically coupled solution of the flow field. The present method is applied on problems of free falling spheres and tori in quiescent flow and freely rotating spheres in simple shear flow. Good agreement with published results shows the ability of the present hybrid meshfree‐Cartesian grid scheme to achieve good accuracy. An application of the method to the self‐induced propulsion of a deforming fish‐like swimmer further demonstrates the capability and potential of the present approach for solving complex FSI problems in 3D. Copyright © 2011 John Wiley & Sons, Ltd.     

A generalized approximation for the meshfree analysis of solids

This paper presents a new approach in the construction of meshfree approximations as well as the weak Kronecker‐delta property at the boundary, referred to as a generalized meshfree (GMF) approximation. The GMF approximation introduces an enriched basis function in the original Shepard's method. This enriched basis function is introduced to meet the linear or higher order reproducing conditions and at the same time to offer great flexibility on the control of the smoothness and convexity of the approximation. The construction of the GMF approximation can be viewed as a special root‐finding scheme of constraint equations that enforces that the basis functions are corrected and the reproducing conditions with certain orders are satisfied within a set of nodes. By choosing different basis functions, various convex and non‐convex approximations including moving least‐squares (MLS), reproducing kernel (RK), and maximum entropy (ME) approximations can be obtained. Furthermore, the basis function can also be translated or blended with other functions to generate a particular approximation for a special purpose. One application in this paper is to incorporate a blending function at the boundary based on the concept of local convexity for the non‐convex approximation, such as MLS, to acquire the weak Kronecker‐delta property. To achieve the higher order GMF approximation, two possible methods are also introduced. Several examples are presented to examine the effectiveness of various GMF approximations. Copyright © 2010 John Wiley & Sons, Ltd.     

Non‐planar 3D crack growth by the extended finite element and level sets—Part II: Level set update

We present a level set method for treating the growth of non‐planar three‐dimensional cracks.The crack is defined by two almost‐orthogonal level sets (signed distance functions). One of them describes the crack as a two‐dimensional surface in a three‐dimensional space, and the second is used to describe the one‐dimensional crack front, which is the intersection of the two level sets. A Hamilton–Jacobi equation is used to update the level sets. A velocity extension is developed that preserves the old crack surface and can accurately generate the growing surface. The technique is coupled with the extended finite element method which approximates the displacement field with a discontinuous partition of unity. This displacement field is constructed directly in terms of the level sets, so the discretization by finite elements requires no explicit representation of the crack surface. Numerical experiments show the robustness of the method, both in accuracy and in treating cracks with significant changes in topology. Copyright © 2002 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.