Crack identification by ‘arrival time’ using XFEM and a genetic algorithm

A computational framework is developed in which cracks in two‐dimensional structures are identified, in conjunction with non‐destructive testing of specimens. As opposed to a previous study by the authors, which was based on time‐harmonic excitation with a single frequency, here the transient response of the structure to a short‐duration signal is measured along part of the external boundary. Crack detection is performed using the solution of an inverse time‐dependent problem. It is shown that the arrival time of the input signal to the points of measurement is a good criterion for crack identification in the time domain. The inverse problem of identification is solved using a genetic algorithm, while each forward problem is solved by the time‐dependent extended finite element method (XFEM). The XFEM scheme is efficient in that it allows the use of a single regular mesh for a large number of forward time response problems with different crack geometries. Numerical examples involving a crack in a flat membrane are presented. Identification based on ‘arrival time’ is shown to perform better than that based on time‐harmonic response. Copyright © 2008 John Wiley & Sons, Ltd.     

Stochastic inverse heat conduction using a spectral approach

An adjoint‐based functional optimization technique in conjunction with the spectral stochastic finite element method is proposed for the solution of an inverse heat conduction problem in the presence of uncertainties in material data, process conditions and measurement noise. The ill‐posed stochastic inverse problem is restated as a conditionally well‐posed L₂ optimization problem. The gradient of the objective function is obtained in a distributional sense by defining an appropriate stochastic adjoint field. The L₂ optimization problem is solved using a conjugate‐gradient approach. Accuracy and effectiveness of the proposed approach is appraised with the solution of several stochastic inverse heat conduction problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Multi‐level hp‐adaptivity for cohesive fracture modeling

Discretization‐induced oscillations in the load–displacement curve are a well‐known problem for simulations of cohesive crack growth with finite elements. The problem results from an insufficient resolution of the complex stress state within the cohesive zone ahead of the crack tip. This work demonstrates that the hp‐version of the finite element method is ideally suited to resolve this complex and localized solution characteristic with high accuracy and low computational effort. To this end, we formulate a local and hierarchic mesh refinement scheme that follows dynamically the propagating crack tip. In this way, the usually applied static a priori mesh refinement along the complete potential crack path is avoided, which significantly reduces the size of the numerical problem. Studying systematically the influence of h‐refinement, p‐refinement, and hp‐refinement, we demonstrate why the suggested hp‐formulation allows to capture accurately the complex stress state at the crack front preventing artificial snap‐through and snap‐back effects. This allows to decrease significantly the number of degrees of freedom and the simulation runtime. Furthermore, we show that by combining this idea with the finite cell method, the crack propagation within complex domains can be simulated efficiently without resolving the geometry by the mesh. Copyright © 2016 John Wiley & Sons, Ltd.     

Nonlinear model order reduction based on local reduced‐order bases

A new approach for the dimensional reduction via projection of nonlinear computational models based on the concept of local reduced‐order bases is presented. It is particularly suited for problems characterized by different physical regimes, parameter variations, or moving features such as discontinuities and fronts. Instead of approximating the solution of interest in a fixed lower‐dimensional subspace of global basis vectors, the proposed model order reduction method approximates this solution in a lower‐dimensional subspace generated by most appropriate local basis vectors. To this effect, the solution space is partitioned into subregions, and a local reduced‐order basis is constructed and assigned to each subregion offline. During the incremental solution online of the reduced problem, a local basis is chosen according to the subregion of the solution space where the current high‐dimensional solution lies. This is achievable in real time because the computational complexity of the selection algorithm scales with the dimension of the lower‐dimensional solution space. Because it is also applicable to the process of hyper reduction, the proposed method for nonlinear model order reduction is computationally efficient. Its potential for achieving large speedups while maintaining good accuracy is demonstrated for two nonlinear computational fluid and fluid‐structure‐electric interaction problems. Copyright © 2012 John Wiley & Sons, Ltd.     

Identification of Gurson–Tvergaard material model parameters via Kalman filtering technique. I. Theory

In this paper quasi-static ductile fracture processes are simulated within the framework of the finite element method by means of the Gurson–Tvergaard isotropic constitutive model for progressively cavitating elastoplastic solids. The progressive degradation of the material strength properties in the fracture process zone due to micro-void growth to coalescence is modeled through the computational cell concept. Among the several model parameters to be calibrated in the computations, attention is restricted to the Tvergaard coefficients q ₁ and q ₂ and to the initial porosity f ₀ in the unstressed configuration. To identify these model parameters the inverse problem is solved via the extended Kalman filter for nonlinear systems coupled to a numerical methodology for the sensitivity analysis. In part I of this work the theory of Kalman filtering and sensitivity analysis is presented. First results concerning the identification of the Tvergaard parameters for a whole crack growth in single edge notched bend specimens made of a pressure vessel steel are presented. In order to enhance the convergence towards the final solution of the identification procedure, during the tests measurements are made of the displacements of points located in the central portion of the notched specimens, where model parameters highly affect the system state variables. In part II of this work a numerical validation of the proposed procedure in terms of uniqueness of the final identified solution, requirements of accuracy for the Bayesian initialization of the model parameters and sensitivity to the experimental measurement errors will be presented and discussed.     

A cut‐cell non‐conforming Cartesian mesh method for compressible and incompressible flow

This paper details a multigrid‐accelerated cut‐cell non‐conforming Cartesian mesh methodology for the modelling of inviscid compressible and incompressible flow. This is done via a single equation set that describes sub‐, trans‐, and supersonic flows. Cut‐cell technology is developed to furnish body‐fitted meshes with an overlapping mesh as starting point, and in a manner which is insensitive to surface definition inconsistencies. Spatial discretization is effected via an edge‐based vertex‐centred finite volume method. An alternative dual‐mesh construction strategy, similar to the cell‐centred method, is developed. Incompressibility is dealt with via an artificial compressibility algorithm, and stabilization achieved with artificial dissipation. In compressible flow, shocks are captured via pressure switch‐activated upwinding. The solution process is accelerated with full approximation storage (FAS) multigrid where coarse meshes are generated automatically via a volume agglomeration methodology. This is the first time that the proposed discretization and solution methods are employed to solve a single compressible–incompressible equation set on cut‐cell Cartesian meshes. The developed technology is validated by numerical experiments. The standard discretization and alternative methods were found equivalent in accuracy and computational cost. The multigrid implementation achieved decreases in CPU time of up to one order of magnitude. Copyright © 2007 John Wiley & Sons, Ltd.     

Adaptive sparse grid based HOPGD: Toward a nonintrusive strategy for constructing space‐time welding computational vademecum

Simulation‐based engineering usually needs the construction of computational vademecum to take into account the multiparametric aspect. One example concerns the optimization and inverse identification problems encountered in welding processes. This paper presents a nonintrusive a posteriori strategy for constructing quasi‐optimal space‐time computational vademecum using the higher‐order proper generalized decomposition method. Contrary to conventional tensor decomposition methods, based on full grids (eg, parallel factor analysis/higher‐order singular value decomposition), the proposed method is adapted to sparse grids, which allows an efficient adaptive sampling in the multidimensional parameter space. In addition, a residual‐based accelerator is proposed to accelerate the higher‐order proper generalized decomposition procedure for the optimal aspect of computational vademecum. Based on a simplified welding model, different examples of computational vademecum of dimension up to 6, taking into account both geometry and material parameters, are presented. These vademecums lead to real‐time parametric solutions and can serve as handbook for engineers to deal with optimization, identification, or other problems related to repetitive task.     

Nondestructive identification of multiple flaws using XFEM and a topologically adapting artificial bee colony algorithm

We present a novel algorithm based on the extended finite element method (XFEM) and an enhanced artificial bee colony (EABC) algorithm to detect and quantify multiple flaws in structures. The concept is based on recent work that have shown the excellent synergy between XFEM, used to model the forward problem, and a genetic‐type algorithm to solve an inverse identification problem and converge to the ‘best’ flaw parameters. In this paper, an adaptive algorithm that can detect multiple flaws without any knowledge on the number of flaws beforehand is proposed. The algorithm is based on the introduction of topological variables into the search space, used to adaptively activate/deactivate flaws during run time until convergence is reached. The identification is based on a limited number of strain sensors assumed to be attached to the structure surface boundaries. Each flaw is approximated by a circular void with the following three variables: center coordinates (x_c, y_c) and radius (r_c), within the XFEM framework. In addition, the proposed EABC scheme is improved by a guided‐to‐best solution updating strategy and a local search (LS) operator of the Nelder–Mead simplex type that show fast convergence and superior global/LS abilities compared with the standard ABC or classic genetic algorithms. Several numerical examples, with increasing level of difficulty, are studied in order to evaluate the proposed algorithm. In particular, we consider identification of multiple flaws with unknown a priori information on the number of flaws (which makes the inverse problem harder), the proximity of flaws, flaws having irregular shapes (similar to artificial noise), and the effect of structured/unstructured meshes. The results show that the proposed XFEM–EABC algorithm is able to converge on all test problems and accurately identify flaws. Hence, this methodology is found to be robust and efficient for nondestructive detection and quantification of multiple flaws in structures. Copyright © 2013 John Wiley & Sons, Ltd.     

Robin‐Neumann transmission conditions for fluid‐structure coupling: Embedded boundary implementation and parameter analysis

Partitioned procedures are appealing for solving complex fluid‐structure interaction (FSI) problems, as they allow existing computational fluid dynamics (CFD) and computational structural dynamics algorithms and solvers to be combined and reused. However, for problems involving incompressible flow and strong added‐mass effect (eg, heavy fluid and slender structure), partitioned procedures suffer from numerical instability, which typically requires additional subiterations between the fluid and structural solvers, hence significantly increasing the computational cost. This paper investigates the use of Robin‐Neumann transmission conditions to mitigate the above instability issue. Firstly, an embedded Robin boundary method is presented in the context of projection‐based incompressible CFD and finite element–based computational structural dynamics. The method utilizes operator splitting and a modified ghost fluid method to enforce the Robin transmission condition on fluid‐structure interfaces embedded in structured non–body‐conforming CFD grids. The method is demonstrated and verified using the Turek and Hron benchmark problem, which involves a slender beam undergoing large transient deformation in an unsteady vortex‐dominated channel flow. Secondly, this paper investigates the effect of the combination parameter in the Robin transmission condition, ie, α_f, on numerical stability and solution accuracy. This paper presents a numerical study using the Turek and Hron benchmark problem and an analytical study using a simplified FSI model featuring an Euler‐Bernoulli beam interacting with a two‐dimensional incompressible inviscid flow. Both studies reveal a trade‐off between stability and accuracy: smaller values of α_f tend to improve numerical stability, yet deteriorate the accuracy of the partitioned solution. Using the simplified FSI model, the critical value of α_f that optimizes this trade‐off is derived and discussed.     

Analysis of three‐dimensional fracture mechanics problems: A two‐scale approach using coarse‐generalized FEM meshes

This paper presents a generalized finite element method (GFEM) based on the solution of interdependent global (structural) and local (crack)‐scale problems. The local problems focus on the resolution of fine‐scale features of the solution in the vicinity of three‐dimensional cracks, while the global problem addresses the macro‐scale structural behavior. The local solutions are embedded into the solution space for the global problem using the partition of unity method. The local problems are accurately solved using an hp‐GFEM and thus the proposed method does not rely on analytical solutions. The proposed methodology enables accurate modeling of three‐dimensional cracks on meshes with elements that are orders of magnitude larger than the process zone along crack fronts. The boundary conditions for the local problems are provided by the coarse global mesh solution and can be of Dirichlet, Neumann or Cauchy type. The effect of the type of local boundary conditions on the performance of the proposed GFEM is analyzed. Several three‐dimensional fracture mechanics problems aimed at investigating the accuracy of the method and its computational performance, both in terms of problem size and CPU time, are presented. Copyright © 2009 John Wiley & Sons, Ltd.     

A dynamic XFEM formulation for crack identification

Nelder–Mead (NM) and Quasi-Newton (QN) optimization methods are used for the numerical solution of crack identification problems in elastodynamics. Fracture is modeled by the eXtended Finite Element Method. The Newmark-β method with Rayleigh damping is employed for the time integration. The effects of various dynamical test loads on the crack identification are investigated. For a time-harmonic excitation with a single frequency and a short-duration signal measured along part of the external boundary, the crack is detected through the solution of an inverse time-dependent problem. Compared to the static load, we show that the dynamic loads are more effective for crack detection problems. Moreover, we tested different dynamic loads and find that NM method works more efficient under the harmonic load than the pounding load while the QN method achieves almost the same results for both load types.     

Thermal Transport Properties of Layered Materials: Identification by a New Numerical Algorithm for Transient Measurements

Transient methods are widely used to determine thermal transport properties. In some situations they can be used for homogeneous media to measure several properties either simultaneously or separately. In this context an analytic model is available and a well-posed inverse problem of parameter identification has to be solved. The examination of composite media is more complicated. The algorithm proposed here allows simultaneous determination of the thermal conductivity and thermal diffusivity of layered dielectrics by transient measurements. It is based on a plane source that acts both as a resistive heater and temperature sensor. For the technique to be successful two essential aspects have to be considered: firstly, the mathematical modeling of the measured data (the forward problem) and secondly, the problem of ill-posedness of the inverse problem. For the proposed measurement configuration, a new fast data analysis algorithm based on an analytic solution for the forward problem is presented. In principle, a numerical solution such as an FEM solution of the heat conduction equation can be used instead of the analytical one, but the computational effort is much greater. The inverse problem is formulated as an output-least-squares problem, which leads to a transcendent algebraic system of equations. The method was successfully tested for different situations.Paper presented at the Fifteenth Symposium on Thermophysical Properties, June 22--27, 2003, Boulder, Colorado, U.S.A.     

Basis mapping methods for forward and inverse problems

This paper describes a novel method for mapping between basis representation of a field variable over a domain in the context of numerical modelling and inverse problems. In the numerical solution of inverse problems, a continuous scalar or vector field over a domain may be represented in different finite‐dimensional basis approximations, such as an unstructured mesh basis for the numerical solution of the forward problem, and a regular grid basis for the representation of the solution of the inverse problem. Mapping between the basis representations is generally lossy, and the objective of the mapping procedure is to minimise the errors incurred. We present in this paper a novel mapping mechanism that is based on a minimisation of the L² or H¹ norm of the difference between the two basis representations. We provide examples of mapping in 2D and 3D problems, between an unstructured mesh basis representative of an FEM approximation, and different types of structured basis including piecewise constant and linear pixel basis, and blob basis as a representation of the inverse basis. A comparison with results from a simple sampling‐based mapping algorithm shows the superior performance of the method proposed here. © 2016 The Authors. International Journal for Numerical Methods in Engineering Published by John Wiley & Sons Ltd.     

Numerical solution for multiple crack problem in an infinite plate under compression

This paper investigates a numerical solution for multiple crack problem in an infinite plate under remote compression. The influence of friction is taken into account. In the first step of the solution, we make a full contact assumption on the crack faces. The full contact assumption means that one component of the dislocation distribution vanishes, and the first mode stress intensity factors (K ₁) at the crack tips become zero. On the above-mentioned assumption, the problem can be solved by using integral equation method, and the second mode stress intensity factors (K ₂) at the crack tips can be evaluated. Meantime, after solving the integral equation the normal contact stress on the crack faces can be evaluated. The next step is to examine the full contact assumption. If the contact stresses on the crack faces are definitely negative, the solution is true. Otherwise, the obtained solution is not true. It is found from present study that in most cases the full contact condition is satisfied, and only in a few cases the full contact condition is violated. Numerical examples are given. It is found that the friction can lower the stress intensity factors at crack tips in general.     

Indirect T‐Trefftz and F‐Trefftz methods for solving boundary value problem of Poisson equation

Abstract

The Poisson equation can be solved by first finding a particular solution and then solving the resulting Laplace equation. In this paper, a computational procedure based on the Trefftz method is developed to solve the Poisson equation for two‐dimensional domains. The radial basis function approach is used to find an approximate particular solution for the Poisson equation. Then, two kinds of Trefftz methods, the T‐Trefftz method and F‐Trefftz method, are adopted to solve the resulting Laplace equation. In order to deal with the possible ill‐posed behaviors existing in the Trefftz methods, the truncated singular value decomposition method and L‐curve concept are both employed. The Poisson equation of the type, ?² u = f(x, u), in which x is the position and u is the dependent variable, is solved by the iterative procedure. Numerical examples are provided to show the validity of the proposed numerical methods and some interesting phenomena are carefully discussed while solving the Helmholtz equation as a Poisson equation. It is concluded that the F‐Trefftz method can deal with a multiply connected domain with genus p(p > 1) while the T‐Trefftz method can only deal with a multiply connected domain with genus 1 if the domain partition technique is not adopted.     

Confidence structural robust optimization by non‐linear semidefinite programming‐based single‐level formulation

Structural robust optimization problems are often solved via the so‐called Bi‐level approach. This solution procedure often involves large computational efforts and sometimes its convergence properties are not so good because of the non‐smooth nature of the Bi‐level formulation. Another problem associated with the traditional Bi‐level approach is that the confidence of the robustness of the obtained solutions cannot be fully assured at least theoretically. In the present paper, confidence single‐level non‐linear semidefinite programming (NLSDP) formulations for structural robust optimization problems under stiffness uncertainties are proposed. This is achieved by using some tools such as S‐procedure and quadratic embedding for convex analysis. The resulted NLSDP problems are solved using the modified augmented Lagrange multiplier method which has sound mathematical properties. Numerical examples show that confidence robust optimal solutions can be obtained with the proposed approach effectively. Copyright © 2010 John Wiley & Sons, Ltd.     

A closed-form solution to the equivalent through-crack model for surface cracks

Summary A closed-form solution for plates containing surface cracks is obtained by using an equivalent through-crack model, which reduces the three-dimensional problem of the surface cracked plates to a two-dimensional Hilbert problem. The effect of surface crack is replaced by a continuous distribution of forcesN(x, 0) and momentsM(x, 0) applied along the crack face of the equivalent through-crack model. A convenient form of expressing these forces and moments is by using power polynomials. Then the singular integral equation, expressing the solution of the Hilbert problem can be readily integrated.According to this model we assumed that the crack depth at extremities, where the crack intersects the free surface of the plate, is not zero. This assumption, corroborated by experimental evidence, means that a singularity exists at the extremities of the crack. The experimental evidence was achieved by using the method of caustics and photoelasticity. The results of the evaluation of the stress intensity factor for elliptical cracks were compared well with the solutions of the respective 3-D problem solved by applying the finite-element method, the line-spring model based on the Reissner plate theory, the finite-element alternating method and the benchmark estimate. The distribution of the stress intensity factor along the crack lips, as calculated in this paper, was dropping off rapidly in the surface layer and it was very close to the results given by the approximate 3-D theory.With 7 Figures     

On the duality of finite element discretization error control in computational Newtonian and Eshelbian mechanics

In this paper, goal-oriented a posteriori error estimators of the averaging type are presented for the error obtained while approximately evaluating theJ-integral in nonlinear elastic fracture mechanics. Since the value of the J-integral is one component of the material force acting on the crack tip of a pre-cracked elastic body, the appropriate mechanical framework to be chosen is the one named after Eshelby rather than classical Newtonian mechanics. However, in a finite element setting, the discretized Eshelby problem is generally not solved explicitly. Rather, its solution is approximated by the finite element solution of the corresponding discretized dual Newton problem. As a consequence, discrete material forces arise not only at the crack tip but also at other nodes of the current finite element mesh. It is the objective of this paper to establish goal-oriented a posteriori error estimators in both the framework of Eshelbian and Newtonian mechanics and to elaborate their dual relations. This allows to control the error of the J-integral while, at the same time, no further discrete material forces arise during the adaptive mesh refinement process which could lead to misleading mechanical interpretations of the results obtained by the finite element method. The paper is concluded by numerical examples that illustrate our theoretical results. Dedicated to the memory of the esteemed colleague Professor Karl Popp, University of Hannover, who unexpectedly passed away on April 24, 2005.     

A multiple-scales approach to crack-front waves

Perturbation of a steadily propagating crack with a straight edge is solved using the method of matched asymptotic expansions (MAE). This provides a simplified analysis in which the inner and outer solutions are governed by distinct mechanics. The inner solution contains the explicit perturbation and is governed by a quasi-static equation. The outer solution determines the radiation of energy away from the tip, and requires solving dynamic equations in the unperturbed configuration. The outer and inner expansions are matched via the small parameter ϵ = L/l defined by the disparate length scales: the crack perturbation length L and the outer length scale l associated with the loading. The method is illustrated for a scalar crack model and then applied to the elastodynamic mode I problem. The crack-front wave-dispersion relation is found by requiring that the energy release rate is unaltered under perturbation and dispersive properties of the crack-front wave speed are described for the first time. The example problems considered demonstrate the potential of MAE for moving-boundary-value problems with multiple scales.     

An inverse solution for reconstruction of the heat transfer coefficient from the knowledge of two temperature values in a solid substrate

Reconstruction of the heat transfer coefficient from the knowledge of temperature distribution is an inverse problem. The main focus of this study was to develop an inverse model that could be used to determine the heat transfer coefficient associated with a fluid in contact with a solid surface from the knowledge of two measured temperature values (T₁ and T_M) in the solid substrate. The temperature distribution for the inverse problem was numerically generated, for a situation with a known heat transfer coefficient, using an implicit finite-differencing scheme. The solution domain was first discretized in to finite number of small regions and nodes. Conservation of energy was then applied to each of the control volume about the nodal regions. This approach resulted in a set of linear equations that was solved simultaneously. Two nodal temperatures in the substrate, from the direct solution, were then used in the inverse problem to reconstruct the heat transfer coefficient. To solve the inverse problem, the solution domain was divided into two distinct regions (Region I and Region II). Region I contained the solution domain between the two known temperatures (T₁ and T_M), and Region II included the region between T_M and the surface with the convective boundary condition. Again, a finite-differencing scheme was employed to generate a set of linear equations in each region. First, the set of linear equations in Region I was solved simultaneously and the results were then used to reconstruct the nodal temperatures in Region II. The convective surface temperature was then utilized to determine the heat transfer coefficient. A series of numerical experiments were conducted to test the validity of the inverse model. Comparison of the inverse solutions with the direct solutions confirms that the heat transfer coefficient can be reconstructed, with good accuracy, from the knowledge of two temperature points in the solid substrate.