Phase field modeling of complex crack patterns in multi-physics problems

Recently developed continuum phase field models for brittle fracture show excellent modeling capability in situations with complex crack topologies including branching in the small and large strain applications. This work presents a generalization towards fully coupled multi-physics problems at large strains. A modular concept is outlined for the linking of the diffusive crack modeling with complex multi field material response, where the focus is put on the model problem of finite thermo-elasticity. This concerns a generalization of crack driving forces from the energetic definitions towards stress-based criteria, the constitutive modeling of degradation of non-mechanical fluxes on generated crack faces. Particular assumptions are made on the generation of convective heat exchanges approximating surface load integrals of the sharp crack approach by distinct volume integrals. The coupling effect is also shown in generation of cracks due to thermally induced stress states. We finally demonstrate the performance of the phase field formulation of fracture at large strains by means of representative numerical examples. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Phase Field Approach for Dynamic Fracture

In phase field fracture models cracks are indicated by the value of a scalar field variable which interpolates smoothly between broken and undamaged material. The evolution equation for this crack field is coupled to the mechanical field equations in order to model the mutual interaction between the crack evolution and mechanical quantities. In finite element simulations of crack growth at comparatively slow loading velocities, a quasi-static phase field model yields reasonable results. However, the simulation of fast loading or the nucleation of new cracks challenges the limits of such a formulation. Here, the quasi-static phase field model predicts brutal crack extension with an artificially high crack speed. In this work, we analyze to which extend a dynamic formulation of the mechanical part of the phase field model can overcome this paradox created by the quasi-static formulation. In finite element simulations, the impact of the dynamic effects is studied, and differences between the crack propagation behavior of the quasi-static model and the dynamic formulation are highlighted. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Variational Phase Field Modeling of Fracture Caused by Fluid Transport in Poroelastic Solids

The numerical modeling of failure mechanisms due to fracture based on sharp crack discontinuities is extremely demanding and suffers in situations with complex crack topologies. This drawback can be overcome by recently developed diffusive crack modeling concepts, which are based on the introduction of a crack phase field. Such an approach is conceptually in line with gradient-extended continuum damage models which include internal length scales. In this paper, we extend our recently outlined mechanical framework [1–3] towards the phase field modeling of fracture in the coupled problem of fluid transport in deforming porous media. Here, extremely complex crack patterns may occur due to drying or hydraulic induced fracture, the so called fracking. We develop new variational potentials for Biot-type fluid transport in porous media at finite deformations coupled with phase field fracture. It is shown, that this complex coupled multi-field problem is related to an intrinsic mixed variational principle for the evolution problem. This principle determines the rates of deformation, fracture phase field and fluid content along with the fluid potential. We develop a robust computational implementation of the coupled problem based on the potentials mentioned above and demonstrate its performance by the numerical simulation of complex fracture patterns. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Advanced finite element formulations for modeling thin piezoelectric structures

This contribution is concerned with mixed finite element formulations for modeling piezoelectric beam and shell structures. Due to the electromechanical coupling, specific deformation modes are joined with electric field components. In bending dominated problems incompatible approximation functions of these fields cause incorrect results. These effects occur in standard finite element formulations, where interpolation functions of lowest order are used. A mixed variational approach is introduced to overcome these problems. The mixed formulation allows for a consistent approximation of the electromechanical coupled problem. It utilizes six independent fields and could be derived from a Hu-Washizu variational principle. Displacements, rotations and the electric potential are employed as nodal degrees of freedom. According to the Timoshenko theory (beam) and the Reissner-Mindlin theory (shell), the formulations account for constant transversal shear strains. To incorporate three dimensional constitutive relations all transversal components of the electric field and the strain field are enriched by mixed finite element interpolations. Thus the complete piezoelectric coupling is appropriately captured. The common assumption of vanishing transversal stress and dielectric displacement components is enforced in an integral sense. Some numerical examples will demonstrate the capability of the presented finite element formulation. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On an Energetic Interpretation of a Phase Field Model for Fracture

In the pioneering work by Griffith, it is assumed that a crack propagates, if this is energetically favorable. However, this original formulation requires a pre-existing initial crack. In order to bypass this deficiency of classical Griffith theory, Francfort and Marigo advocate a global variational criterion, where the total energy is minimized with respect to any admissible displacement field and crack set. Bourdin's regularized approximation of this variational formulation makes use of a continuous scalar field to indicate cracks. Based on this regularization a phase field fracture model is formulated. The crack field is assumed to follow a Ginzburg-Landau type evolution equation, and cracking is addressed as a phase transition problem. The coupled problem of mechanical balance equations and the evolution equation is solved using the finite element method combined with an implicit time integration scheme. The numerical solution naturally yields the crack evolution including crack propagation, kinking, branching and initiation without any additional criteria. In this work we study the driving mechanisms behind the crack evolution in the phase field fracture model and compare to the purely energetic considerations of the underlying variational formulation. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Phase Field Model for Ductile to Brittle Failure Mode Transition

The computational modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in dynamic problems with complex crack topologies. This can be overcome by a diffusive crack modeling based on the introduction of a crack phase field. We outline a conceptual framework for phase field models of crack propagation in brittle elastic and ductile elastic-plastic solids under dynamic loading and investigate the ductile to brittle failure mode transition observed in the experiment performed by Kalthoff and Winkeler [3]. We develop incremental variational principles and consider their numerical implementations by multi-field finite element methods. To this end, we define energy storage and dissipation functions for the plastic flow including the fracture phase field. The introduction of local history fields that drive the evolution of the crack phase field inspires the construction of robust operator split schemes. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Peridynamics model for flexoelectricity and damage

A flexoelectric peridynamic (PD) theory is proposed. In the PD framework, the formulation introduces a nanoscale flexoelectric coupling that entails non-uniform strain in centrosymmetric dielectrics. This potentially enables PD modeling of a large class of phenomena in solid dielectrics involving cracks, discontinuities etc. wherein large strain gradients are present and the classical electromechanical theory based on partial differential equations do not directly apply. PD electromechanical equations, derived from Hamilton's principle, satisfy the global balance laws. Linear PD constitutive equations reflect the electromechanical coupling effect, with the mechanical force state affected by the polarization state and the electrical force state in turn by the displacement state. An analytical solution to the PD electromechanical equations is presented for the static case when a point mechanical force and a point electric force act in an infinite 3D solid dielectric. A parametric study on how different length scales influence the response is undertaken. In addition, the model is extended to incorporate damage using phase field – an order parameter, supplemented with a PD bond breaking criterion to study flexoelectric effects in damage and fracture problems. To demonstrate the performance of our proposal, we first simulate, considering small flexoelectricity effect and no damage, an externally pressured 2D flexoelectric disk subjected to a potential difference between the inner and outer surfaces and compare the results with existing solutions in the literature. Next, we simulate a plate with a central pre-crack under tension considering damage and flexoelectricity effects, and study the effect of various constitutive parameters on the damage evolution. We also furnish a classical derivation of phase field based flexoelectricity in Appendix I.     

Radiation Back-Reaction and Renormalization in Classical Field Theory with Singular Sources

We propose a general covariant method for regularizing the radiation back-reaction in linear and nonlinear field theory models with singular sources. Typical examples of such sources are the currents produced by extended relativistic objects (branes). As an illustration, we consider the models of minimal and nonminimal coupling of a brane to an n-form gauge field, a scalar field, and the Einstein gravity field. We find the structure of divergent and finite contributions due to the radiation back-reaction and obtain relations for the parameters of the theory ensuring cancellation of divergences. We prove that the divergences are Lagrangian n the case where the metric induced on the brane surface is nondegenerate. We find special types of a (nonminimal) coupling leading to local and Lagrangian effective equations of motion of the brane. We show that the requirement for classical renormalizability imposes strong restrictions on the self-coupling vertices of the field, similar to the quantum renormalizability conditions. In particular, we establish the nonrenormalizability of the gravitational self-coupling of a codimension-(k>2) brane, whereas for k ≤ 2, the theory becomes not only renormalizable but also finite.__________Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 143, No. 3, pp. 375–400, June, 2005.     

Phase-Field Modeling of Hydraulic Fracture

Hydraulically driven fracture has gained more and more research activity in the last few years, especially due to the growing interest of the petroleum industry. Key challenge for a powerful simulation of this scenario is an effective modeling and numerical implementation of the behavior of the solid skeleton and the fluid phase, the mechanical coupling between the two phases as well as the incorporation of the fracture process. Existing models for hydraulic fracturing can be found for example in [1], where the crack path is predetermined, or in [2] who use a phase field fracture model in an elastic framework, however without incorporating the fluid flow. In this work we propose a new compact model structure for the Biot-type fluid transport in porous media at finite strains based on only two constitutive functions, that is the free energy function ψ and a dissipation potential ϕ that includes the incorporation of an additional Poiseuille-type fluid flow in cracks. This formulation is coupled to a phase field approach for fracture and is fully variational in nature, as shown in [3]. In contrast to formulations with a sharp-crack discontinuity, the proposed regularized approach has the main advantage of a straight-forward modeling of complex crack patterns including branching. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A finite element model for laminated piezoelectric shell structures

Thin piezoelectric laminates are used for a wide range of technical applications. A four-node piezoelectric shell element is presented to analyse such structures effectively. In case of bending dominated problems incompatible approximation functions of the electrical and mechanical fields cause incorrect results. In order to overcome this problem the finite element formulation is based on a mixed variational principle implying six independent fields: displacements, electric potential, strains, electric field, mechanical stresses and dielectric displacements. This allows for an interpolation of the strains and the electric field in thickness direction independent of the bilinear interpolation functions. A piecewise quadratic approach for the shear strains in thickness direction and the corresponding electric field is proposed for arbitrarily layered shells. Regarding coupling of electrical and mechanical fields this yields to an appropriate balance of the approximation functions. Numerical examples show more precise results in contrast to standard elements with lowest order interpolation functions. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Configurational-Force-Based Adaptive FE Solver for a Phase Field Model of Fracture

The computational modeling of failure mechanisms in solids due to fracture based on sharp crack discontinuities suffers in situations with complex crack topologies. This can be overcome by diffusive crack modeling, based on the introduction of a crack phase field as outlined in [1, 2]. Following these formulations, we outline a thermodynamically consistent framework for phase field models of crack propagation in elastic solids, develop incremental variational principles and, as an extension to [1, 2], consider their numerical implementations by an efficient h-adaptive finite element method. A key problem of the phase field formulation is the mesh density, which is required for the resolution of the diffusive crack patterns. To this end, we embed the computational framework into an adaptive mesh refinement strategy that resolves the fracture process zones. We construct a configurational-force-based framework for h-adaptive finite element discretizations of the gradient-type diffusive fracture model. We develop a staggered computational scheme for the solution of the coupled balances in physical and material space. The balance in the material space is then used to set up indicators for the quality of the finite element mesh and accounts for a subsequent h-type mesh refinement. The capability of the proposed method is demonstrated by means of a numerical example. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A higher order finite element including transverse normal strain for linear elastic composite plates with general lamination configurations

This paper describes a higher-order global-local theory for thermal/mechanical response of moderately thick laminated composites with general lamination configurations. In-plane displacement fields are constructed by superimposing the third-order local displacement field to the global cubic displacement field. To eliminate layer-dependent variables, interlaminar shear stress compatibility conditions have been employed, so that the number of variables involved in the proposed model is independent of the number of layers of laminates. Imposing shear stress free condition at the top and the bottom surfaces, derivatives of transverse displacement are eliminated from the displacement field, so that C⁰ interpolation functions are only required for the finite element implementation. To assess the proposed model, the quadratic six-node C⁰ triangular element is employed for the interpolation of all the displacement parameters defined at each nodal point on the composite plate. Comparing to various existing laminated plate models, it is found that simple C⁰ finite elements with non-zero normal strain could produce more accurate displacement and stresses for thick multilayer composite plates subjected to thermal and mechanical loads. Finally, it is remarked that the proposed model is quite robust, such that the finite element results are not sensitive to the mesh configuration and can rapidly converge to 3-D elasticity solutions using regular or irregular meshes.     

Numerical and experimental investigations of phase segragation effects in binary brazing solders

Solder materials occupy many of fields for technical application (e.g. solder joints in automotive control units or in microelectronic packages). They are required to provide electrical and mechanical connections between different components. Due to their wide range of applications solder alloys are subject to a great variety of microstructural changes such as phase separation and coarsening processes. The micromorphological variations influence strength and life expectation of solder materials, in particular, in very small components such as solder joints in microelectronic packages. In order to analyze the microstructural evolution with a diffusion theory of heterogeneous solid mixtures we employ an extended Cahn-Hilliard phase field model. The diffusion equation under consideration constitutes a partial differential equation involving spatial derivatives of fourth order. Thus, the variational formulation of the problem requires approximation functions which are piecewise smooth and globally C1-continuous. In our contribution we fulfil the continuity requirement by means of rational B-spline finite element basis functions. To illustrate the versatility of this approach numerical simulations of phase decomposition and coarsening controlled by diffusion and by mechanical loading are discussed and compared with experimental results. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On Multivariate Interpolation

A new approach to interpolation theory for functions of several variables is proposed. We develop a multivariate divided difference calculus based on the theory of noncommutative quasi-determinants. In addition, intriguing explicit formulae that connect the classical finite difference interpolation coefficients for univariate curves with multivariate interpolation coefficients for higher dimensional submanifolds are established.     

Phase Field Modeling of Brittle and Ductile Fracture

The phase field modeling of brittle fracture was a topic of intense research in the last few years and is now well-established. We refer to the work [1-3], where a thermodynamically consistent framework was developed. The main advantage is that the phase-field-type diffusive crack approach is a smooth continuum formulation which avoids the modeling of discontinuities and can be implemented in a straightforward manner by multi-field finite element methods. Therefore complex crack patterns including branching can be resolved easily. In this paper, we extend the recently outlined phase field model of brittle crack propagation [1-3] towards the analysis of ductile fracture in elastic-plastic solids. In particular, we propose a formulation that is able to predict the brittle-to-ductile failure mode transition under dynamic loading that was first observed in experiments by Kalthoff and Winkler [4]. To this end, we outline a new thermodynamically consistent framework for phase field models of crack propagation in ductile elastic-plastic solids under dynamic loading, develop an incremental variational principle and consider its robust numerical implementation by a multi-field finite element method. The performance of the proposed phase field formulation of fracture is demonstrated by means of the numerical simulation of the classical Kalthoff-Winkler experiment that shows the dynamic failure mode transition. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A conductive crack and a remote electrode at the interface between two piezoelectric materials

An interaction of a tunnel conductive crack and a distant strip electrode situated at the interface between two piezoelectric semi-infinite spaces is studied. The bimaterial is subject by an in-plane electrical field parallel to the interface and by an anti-plane mechanical loading. Using the presentations of electromechanical quantities at the interface via sectionally-analytic functions the problem is reduced to a combined Dirichlet-Riemann boundary value problem. Solution of this problem is found in an analytical form excepting some one-dimensional integrals calculations. Closed form expressions for the stress, the electric field and their intensity factors, as well as for the crack faces displacement jump are derived. On the base of these presentations the energy release rate is also found. The obtained solution is compared with simple particular case of a single crack without electrode and the excellent agreement is found out. An auxiliary plane problem for open and closed cracks between two isotropic materials is also considered. The mathematical model of this problem is identical to the above one, therefore, the obtained solution is used for this model. It is compared with finite element solution of a similar problem and good agreement is found out.     

The scalar Nevanlinna-Pick interpolation problem with boundary conditions

We show that if the Nevanlinna-Pick interpolation problem is solvable by a function mapping into a compact subset of the unit disc, then the problem remains solvable with the addition of any number of boundary interpolation conditions, provided the boundary interpolation values have modulus less than unity. We give new, inductive proofs of the Nevanlinna-Pick interpolation problem with any finite number of interpolation points in the interior and on the boundary of the domain of interpolation (the right half plane or unit disc), with function values and any finite number of derivatives specified. Our solutions are analytic on the closure of the domain of interpolation. Our proofs only require a minimum of matrix theory and operator theory. We also give new, straightforward algorithms for obtaining minimal H_∞ norm solutions. Finally, some numerical examples are given.     

Finite deformations of thin shells with piezoelectric sensors and actuators

We treat coupled electromechanical problem of finite deformations of piezoelectric shells with the help of the direct approach. A shell is considered as a material surface with mechanical degrees of freedom of particles and with an additional field variable, namely electric potential on the electrodes. This results both in the nonlinear system of equations of piezoelectric shells and in the appropriate numerical scheme. Application of the direct approach is preceded with the three-dimensional asymptotic analysis of a linear electromechanical problem for a non-homogeneous piezoelectric plate, which provides the constitutive relations for the nonlinear theory. As a sample problem, we present finite element analysis of deformation and local buckling of cylindrical panel, equipped with piezoelectric sensors. The latter influence the mechanical behavior and produce signals, which can be interpreted in terms of structural entities. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence phenomena of \documentclass{article} \usepackage{amsmath,amsfonts}\pagestyle{empty}\begin{document} $\mathcal{Q}_p$ \end{document}‐elements for convection–diffusion problems

We present a numerical study for singularly perturbed convection–diffusion problems using higher order Galerkin and streamline diffusion finite element method. We are especially interested in convergence and superconvergence properties with respect to different interpolation operators. For this we investigate pointwise interpolation and vertex‐edge‐cell interpolation. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A Phase-Field Model of Ductile Fracture at Finite Strains

This work outlines a variational-based framework for the phase field modeling of ductile fracture in elastic-plastic solids at large strains. The phase field approach regularizes sharp crack discontinuities within a pure continuum setting by a specific gradient damage model with geometric features rooted in fracture mechanics. Based on the recent works [1, 2], the phase field model of ductile fracture is linked to a formulation of gradient plasticity at finite strains in order to ensure the crack to evolve inside the plastic zones. The thermodynamic formulation is based on the definition of a constitutive work density function including the stored elastic energy and the dissipated work due to plasticity and fracture. The proposed canonical theory is shown to be governed by a rate-type minimization principle, which determines the coupled multi-field evolution problem. Another aspect is the regularization towards a micromorphic gradient plasticity-damage setting which enhances the robustness of the finite element formulation. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)