Interface dynamic stiffness matrix approach for three‐dimensional transient multi‐region boundary element analysis

A novel substructuring method is developed for the coupling of boundary element and finite element subdomains in order to model three‐dimensional multi‐region elastodynamic problems in the time domain. The proposed procedure is based on the interface stiffness matrix approach for static multi‐region problems using variational principles together with the concept of Duhamel integrals. Unit impulses are applied at the boundary of each region in order to evaluate the impulse response matrices of the Duhamel (convolution) integrals. Although the method is not restricted to a special discretization technique, the regions are discretized using the boundary element method combined with the convolution quadrature method. This results in a time‐domain methodology with the advantages of performing computations in the Laplace domain, which produces very accurate and stable results as verified on test examples. In addition, the assembly of the boundary element regions and the coupling to finite elements are greatly simplified and more efficient. Finally, practical applications in the area of soil–structure interaction and tunneling problems are shown. Copyright © 2009 John Wiley & Sons, Ltd.     

Higher‐order multi‐resolution topology optimization using the finite cell method

This article presents a detailed study on the potential and limitations of performing higher‐order multi‐resolution topology optimization with the finite cell method. To circumvent stiffness overestimation in high‐contrast topologies, a length‐scale is applied on the solution using filter methods. The relations between stiffness overestimation, the analysis system, and the applied length‐scale are examined, while a high‐resolution topology is maintained. The computational cost associated with nested topology optimization is reduced significantly compared with the use of first‐order finite elements. This reduction is caused by exploiting the decoupling of density and analysis mesh, and by condensing the higher‐order modes out of the stiffness matrix. Copyright © 2016 John Wiley & Sons, Ltd.     

A finite element approach combining a reduced‐order system,Padé approximants,and an adaptive frequency windowing for fast multi‐frequency solution of poro‐acoustic problems

In this work, a solution strategy is investigated for the resolution of multi‐frequency structural‐acoustic problems including 3D modeling of poroelastic materials. The finite element method is used, together with a combination of a modal‐based reduction of the poroelastic domain and a Padé‐based reconstruction approach. It thus takes advantage of the reduced‐size of the problem while further improving the computational efficiency by limiting the number of frequency resolutions of the full‐sized problem. An adaptive procedure is proposed for the discretization of the frequency range into frequency intervals of reconstructed solution. The validation is presented on a 3D poro‐acoustic example. Copyright © 2013 John Wiley & Sons, Ltd.     

A fast implementation of the FETI‐DP method: FETI‐DP‐RBS‐LNA and applications on large scale problems with localized non‐linearities

As parallel and distributed computing gradually becomes the computing standard for large scale problems, the domain decomposition method (DD) has received growing attention since it provides a natural basis for splitting a large problem into many small problems, which can be submitted to individual computing nodes and processed in a parallel fashion. This approach not only provides a method to solve large scale problems that are not solvable on a single computer by using direct sparse solvers but also gives a flexible solution to deal with large scale problems with localized non‐linearities. When some parts of the structure are modified, only the corresponding subdomains and the interface equation that connects all the subdomains need to be recomputed. In this paper, the dual–primal finite element tearing and interconnecting method (FETI‐DP) is carefully investigated, and a reduced back‐substitution (RBS) algorithm is proposed to accelerate the time‐consuming preconditioned conjugate gradient (PCG) iterations involved in the interface problems. Linear–non‐linear analysis (LNA) is also adopted for large scale problems with localized non‐linearities based on subdomain linear–non‐linear identification criteria. This combined approach is named as the FETI‐DP‐RBS‐LNA algorithm and demonstrated on the mechanical analyses of a welding problem. Serial CPU costs of this algorithm are measured at each solution stage and compared with that from the IBM Watson direct sparse solver and the FETI‐DP method. The results demonstrate the effectiveness of the proposed computational approach for simulating welding problems, which is representative of a large class of three‐dimensional large scale problems with localized non‐linearities. Copyright © 2005 John Wiley & Sons, Ltd.     

A computational framework for scale‐bridging in multi‐scale simulations

A computational framework for scale‐bridging in multi‐scale simulations is presented. The framework enables seamless combination of at‐scale models into highly dynamic hierarchies to build a multi‐scale model. Its centerpiece is formulated as a standalone module capable of fully asynchronous operation. We assess its feasibility and performance for a two‐scale model applied to two challenging test problems from impact physics. We find that the computational cost associated with using the framework may, as expected, become substantial. However, the framework has the ability of effortlessly combining at‐scale models to render complex multi‐scale models. The main source of the computational inefficiency of the framework is related to poor load balancing of the lower‐scale model evaluation We demonstrate that the load balancing can be efficiently addressed by recourse to conventional load‐balancing strategies. Copyright © 2016 John Wiley & Sons, Ltd.     

A time‐staggered partitioned coupling algorithm for transient heat conduction

We present a time‐staggered partitioned coupling algorithm for transient heat conduction finite element simulations. This algorithm divides a large structural mesh into a number of smaller subdomains, solves the individual subdomains separately and couples the solutions to obtain the response to the original problem. The proposed algorithm is a mixed multi‐timestep algorithm and enables arbitrary time integration schemes and meshes to be coupled with different timesteps in the various subdomains. In this procedure, the solution of each partition is separately evaluated over a system timestep after which the interfacial conditions are enforced making this a staggered algorithm that facilitates parallel computation. We present examples showing the feasibility of the coupling algorithm and discuss the merits in terms of convergence and stability. Copyright © 2009 John Wiley & Sons, Ltd.     

A hybrid variational‐collocation immersed method for fluid‐structure interaction using unstructured T‐splines

We present a hybrid variational‐collocation, immersed, and fully‐implicit formulation for fluid‐structure interaction (FSI) using unstructured T‐splines. In our immersed methodology, we define an Eulerian mesh on the whole computational domain and a Lagrangian mesh on the solid domain, which moves arbitrarily on top of the Eulerian mesh. Mathematically, the problem reduces to solving three equations, namely, the linear momentum balance, mass conservation, and a condition of kinematic compatibility between the Lagrangian displacement and the Eulerian velocity. We use a weighted residual approach for the linear momentum and mass conservation equations, but we discretize directly the strong form of the kinematic relation, deriving a hybrid variational‐collocation method. We use T‐splines for both the spatial discretization and the information transfer between the Eulerian mesh and the Lagrangian mesh. T‐splines offer us two main advantages against non‐uniform rational B‐splines: they can be locally refined and they are unstructured. The generalized‐α method is used for the time discretization. We validate our formulation with a common FSI benchmark problem achieving excellent agreement with the theoretical solution. An example involving a partially immersed solid is also solved. The numerical examples show how the use of T‐junctions and extraordinary nodes results in an accurate, efficient, and flexible method. Copyright © 2015 John Wiley & Sons, Ltd.     

Filters in topology optimization based on Helmholtz‐type differential equations

The aim of this paper is to apply a Helmholtz‐type partial differential equation as an alternative to standard density filtering in topology optimization problems. Previously, this approach has been successfully applied as a sensitivity filter. The usual filtering techniques in topology optimization require information about the neighbor cells, which is difficult to obtain for fine meshes or complex domains and geometries. The complexity of the problem increases further in parallel computing, when the design domain is decomposed into multiple non‐overlapping partitions. Obtaining information from the neighbor subdomains is an expensive operation. The proposed filter technique requires only mesh information necessary for the finite element discretization of the problem. The main idea is to define the filtered variable implicitly as a solution of a Helmholtz‐type differential equation with homogeneous Neumann boundary conditions. The properties of the filter are demonstrated for various 2D and 3D topology optimization problems in linear elasticity, solved on serial and parallel computers. Copyright © 2010 John Wiley & Sons, Ltd.     

A neural network graph partitioning procedure for grid‐based domain decomposition

This paper describes a neural network graph partitioning algorithm which partitions unstructured finite element/volume meshes as a precursor to a parallel domain decomposition solution method. The algorithm works by first constructing a coarse graph approximation using an automatic graph coarsening method. The coarse graph is partitioned and the results are interpolated onto the original graph to initialize an optimization of the graph partition problem. In practice, a hierarchy of (usually more than two) graphs are used to help obtain the final graph partition. A mean field theorem neural network is used to perform all partition optimization. The partitioning method is applied to graphs derived from unstructured finite element meshes and in this context it can be viewed as a multi‐grid partitioning method. Copyright © 1999 John Wiley & Sons, Ltd.     

An energy‐momentum‐conserving temporal discretization scheme for adhesive contact problems

Numerical solution of dynamic problems requires accurate temporal discretization schemes. So far, to the best of the authors’ knowledge, none have been proposed for adhesive contact problems. In this work, an energy‐momentum‐conserving temporal discretization scheme for adhesive contact problems is proposed. A contact criterion is also proposed to distinguish between adhesion‐dominated and impact‐dominated contact behaviors. An adhesion formulation is considered, which is suitable to describe a large class of interaction mechanisms including van der Waals adhesion and cohesive zone modeling. The current formulation is frictionless, and no dissipation is considered. Performance of the proposed scheme is compared with other schemes. The proposed scheme involves very little extra computational overhead. It is shown that the proposed new temporal discretization scheme leads to major accuracy gains both for single‐degree‐of‐freedom and multi‐degree‐of‐freedom systems. The single‐degree‐of‐freedom system is critically analyzed for various parameters affecting the response. For the multi‐degree‐of‐freedom system, the effect of the time step and mesh discretization on the solution is also studied using the proposed scheme. It is further shown that a temporal discretization scheme based on the principle of energy conservation is not sufficient to obtain a convergent solution. Results with higher order contact finite elements for discretizing the contact area are also discussed. Copyright © 2012 John Wiley & Sons, Ltd.     

Analysis of singular stress fields at multi‐material corners under thermal loading

The scaled boundary finite‐element method is extended to the modelling of thermal stresses. The particular solution for the non‐homogeneous term caused by thermal loading is expressed as integrals in the radial direction, which are evaluated analytically for temperature changes varying as power functions of the radial coordinate. When applied to model a multi‐material corner, only the boundary of the problem domain is discretized. The boundary conditions on the straight material interfaces and the side‐faces forming the corner are satisfied analytically without discretization. The stress field is expressed semi‐analytically as a series solution. The stress distribution along the radial direction, including both the real and complex power singularity and the power‐logarithmic singularity, is represented analytically. The stress intensity factors are determined directly from their definitions in stresses. No knowledge on asymptotic expansions is required. Numerical examples are calculated to evaluate the accuracy of the scaled boundary finite‐element method. Copyright © 2005 John Wiley & Sons, Ltd.     

Adaptive atomistic‐to‐continuum modeling of propagating defects

An adaptive atomistic‐to‐continuum method is presented for modeling the propagation of material defects. This method extends the bridging domain method to allow the atomic domain to dynamically conform to the evolving defect regions during a simulation, without introducing spurious oscillations and without requiring mesh refinement. The atomic domain expands as defects approach the bridging domain method coupling domain by fine graining nearby finite elements into equivalent atomistic subdomains. Additional algorithms coarse grain portions of the atomic domain to the continuum scale, reducing the degrees of freedom, when the atomic displacements in a subdomain can be approximated by FEM or extended FEM elements to within a certain homogeneity tolerance. The extended FEM approximations are created by fitting the broken inter‐atomic bonds of fractured surfaces and dislocation slip planes. Because atomic degrees of freedom are maintained only where needed for each timestep, the solution retains the advantages of multiscale modeling, with a reduced computational cost compared with other multiscale methods. Copyright © 2012 John Wiley & Sons, Ltd.     

Improved multi‐level Newton solvers for fully coupled multi‐physics problems

An algorithm is suggested to improve the efficiency of the multi‐level Newton method that is used to solve multi‐physics problems. It accounts for full coupling between the subsystems by using the direct differentiation method rather than error prone finite difference calculations and retains the advantage of greater flexibility over the tightly coupled approaches. Performance of the algorithm is demonstrated by solving a fluid–structure interaction problem. Copyright © 2003 John Wiley & Sons, Ltd.     

A fast multi‐level convolution boundary element method for transient diffusion problems

A new algorithm is developed to evaluate the time convolution integrals that are associated with boundary element methods (BEM) for transient diffusion. This approach, which is based upon the multi‐level multi‐integration concepts of Brandt and Lubrecht, provides a fast, accurate and memory efficient time domain method for this entire class of problems. Conventional BEM approaches result in operation counts of order O(N²) for the discrete time convolution over N time steps. Here we focus on the formulation for linear problems of transient heat diffusion and demonstrate reduced computational complexity to order O(N^3/2) for three two‐dimensional model problems using the multi‐level convolution BEM. Memory requirements are also significantly reduced, while maintaining the same level of accuracy as the conventional time domain BEM approach. Copyright © 2005 John Wiley & Sons, Ltd.     

Adaptive Kronrod–Patterson integration of non‐linear finite element matrices

Efficient simulation of unsaturated moisture flow in porous media is of great importance in many engineering fields. The highly non‐linear character of unsaturated flow typically gives sharp moving moisture fronts during wetting and drying of materials with strong local moisture permeability and capacity variations as result. It is shown that these strong variations conflict with the common preference for low‐order numerical integration in finite element simulations of unsaturated moisture flow: inaccurate numerical integration leads to errors that are often far more important than errors from inappropriate discretization. In response, this article develops adaptive integration, based on nested Kronrod–Patterson–Gauss integration schemes: basically, the integration order is adapted to the locally observed grade of non‐linearity. Adaptive integration is developed based on a standard infiltration problem, and it is demonstrated that serious reductions in the numbers of required integration points and discretization nodes can be obtained, thus significantly increasing computational efficiency. The multi‐dimensional applicability is exemplified with two‐dimensional wetting and drying applications. While developed for finite element unsaturated moisture transfer simulation, adaptive integration is similarly applicable for other non‐linear problems and other discretization methods, and whereas perhaps outperformed by mesh‐adaptive techniques, adaptive integration requires much less implementation and computation. Both techniques can moreover be easily combined. Copyright © 2009 John Wiley & Sons, Ltd.     

Coupled multi‐scale cohesive modeling of failure in heterogeneous adhesives

A multi‐scale cohesive numerical framework is proposed to simulate the failure of heterogeneous adhesively bonded systems. This multi‐scale scheme is based on Hill's variational principle of energy equivalence between the higher and lower level scales. It provides an easy way to obtain accurate homogenized macroscopic properties while capturing the physics of failure processes at the micro‐scale in sufficient detail. We use an isotropic rate‐dependent damage model to mimic the failure response of the constituents of heterogeneous adhesives. The finite element method is used to solve the equilibrium equation at each scale. A nested iterative scheme inspired by the return mapping algorithm used in computational inelasticity is implemented. We propose a computationally attractive technique to couple the macro‐ and micro‐scales for rate‐dependent constitutive laws. We introduce an adhesive patch test to study the numerical performance, including spatial and temporal convergence of the multi‐scale scheme. We compare the solution of the multi‐scale cohesive scheme with a direct numerical simulation. Finally, we solve mode I and mode II fracture problems to demonstrate failure at the macro‐scale. Copyright © 2010 John Wiley & Sons, Ltd.     

Voxel‐based meshing and unit‐cell analysis of textile composites

Unit‐cell homogenization techniques are frequently used together with the finite element method to compute effective mechanical properties for a wide range of different composites and heterogeneous materials systems. For systems with very complicated material arrangements, mesh generation can be a considerable obstacle to usage of these techniques. In this work, pixel‐based (2D) and voxel‐based (3D) meshing concepts borrowed from image processing are thus developed and employed to construct the finite element models used in computing the micro‐scale stress and strain fields in the composite. The potential advantage of these techniques is that generation of unit‐cell models can be automated, thus requiring far less human time than traditional finite element models. Essential ideas and algorithms for implementation of proposed techniques are presented. In addition, a new error estimator based on sensitivity of virtual strain energy to mesh refinement is presented and applied. The computational costs and rate of convergence for the proposed methods are presented for three different mesh‐refinement algorithms: uniform refinement; selective refinement based on material boundary resolution; and adaptive refinement based on error estimation. Copyright © 2003 John Wiley & Sons, Ltd.