Analysis of accuracy of Williams series approximation of stress field in cracked body – influence of area of interest around crack-tip on multi-parameter regression performance

J. Sobek et alii, Frattura ed Integrità Strutturale, 39 (2017) 129-142; DOI: 10.3221/IGF-ESIS.39.14 129 Focussed on Modelling in Mechanics Analysis of accuracy of Williams series approximation of stress field in cracked body – influence of area of interest around crack-tip on multi-parameter regression performance J. Sobek, P. Frantík, V. Veselý Brno University of Technology, Faculty of Civil Engineering, Institute of Structural Mechanics, Brno, Czech Republic [email protected], http://orcid.org/0000-0003-4215-1029 [email protected] vesely.v1@ fce.vutbr.cz, http://orcid.org/0000-0002-7723-971X ABSTRACT. A study on the accuracy of an approximation of the stress field in a cracked body is presented. Crack-tip stress tensor is expressed using the linear elastic fracture mechanics (LEFM) theory in this work, more precisely via its multi-parameter formulation, i.e. by Williams power series (WPS). Determination of coefficients of terms of this series is performed using a least squares-based regression technique known as over-deterministic method (ODM) for which results from finite element (FE) method computation are usually taken as inputs. Main attention is paid to a detailed analysis of a suitable selection of FE nodes whose results serve as the inputs to the employed method. Two different ways of FE nodal selection are compared – nodes selected from the crack tip vicinity lying at a ring of a certain radius versus nodes selected more or less uniformly from a specified part of the test specimen body. Comparison of these approaches is made with the help of procedures developed by the authors which enable both the determination of the coefficients of terms of the analytical WPS approximation of the stress field based on the FE results and the backward reconstruction of the field (again using WPS) from those determined terms’ coefficients/functions. The wedge-splitting test (WST) specimen with a crack is taken as example for the study. KEYWORDS. Multi-parameter fracture mechanics; Williams power series; Crack-tip fields; Over-deterministic method; Higher order terms. Citation: Sobek, J., Frantík, P., Veselý, V., Analysis of accuracy of Williams series approximation of stress field in cracked body – influence of area of interest around cracktip on multi-parameter regression performance, Frattura ed Integrità Strutturale, 39 (2017) 129-142. Received: 20.07.2016 Accepted: 19.09.2016 Published: 01.01.2017 Copyright: © 2017 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION etermination of the higher order terms of Williams power series [1] approximating the fields of stress and displacements in a cracked body is of substantial interest of the presented work. Particular motivation of it is represented by the need for capturing the crack-tip fields in a farther distance from the crack tip. D J. Sobek et alii, Frattura ed Integrità Strutturale, 39 (2017) 129-142; DOI: 10.3221/IGF-ESIS.39.14 130 Tensile failure of the quasi-brittle materials is accompanied with a crack propagation together with a nonlinear zone [2] (fracture process zone – FPZ) development, where the decrease of material integrity takes place. The size of the FPZ is not negligible in comparison with the rest of the body. Processes of failure therefore occur in a wider area around the crack tip. The description of failure mechanism must be in agreement with the stress and displacement field around the crack tip in an extended area. Recent works of many authors show the relevance of the topic of the crack tip fields accurate description [3–6], moreover, extending it into 3D and taking into account the effects of various loading modes [7–12]. Special focus on brittle and quasi-brittle fracture is summarized in [13]. Works [3,14,15] reported the fact, that for the description of the stress/displacement field in a cracked body in a more distant surroundings from the crack tip the necessity of usage of the several terms of Williams expansion (WE), not only the first or the first two terms, is crucial. Procedures enabling the multi-parameter description of the near crack-tip fields (using e.g. hybrid crack elements [16] formulation, over-deterministic method [17] based on standard FE computation, or other techniques based e.g. on extrapolation of displacements of selected nodes of FE mesh) usually process results from FE nodes selected from the close vicinity of the crack tip. This is adequate for determination of the classical/two-parameter LEFM characteristics (SIF, T-stress). However, several questions arise if the area of accurate enough description of the stress and displacement field extends to a greater distance from the crack tip, where K + T dominance vanishes? How many terms of the series should be taken into account? How to select the FE nodes considered for the regression technique? And how to optimize the mutual relationship between the area from where the nodes are considered for the regression (and how are they located/distributed in that area) and the extent where the approximation of the fields is of relevant accuracy? Answering these questions presents the actual motivation of this work. This paper investigates the above-described issue via a parametric study evaluating the influence of the nodal selection on the quality of the obtained approximation of the field. The work presented here further builds on previous studies. A classical (common) way of nodal selection, where the nodes are selected from a ring in the vicinity of the crack tip was used in recent papers. Influences of several parameters on the description of stress/displacement fields in cracked bodies by Williams series were investigated. Reconstruction of the stress field with the help of a software tool developed by the author’s team was shown in [14,15] where the accuracy of the approximation by WE was verified by a visual comparison (which was regarded as sufficient for its intended purpose) with the FE solution (which was regarded as the exact solution). However, the nodal selection that was used for obtaining the WE terms was performed from the ring around the crack tip (with the distance given by recommendation from [17,18]). Published study [19] on WPS approximation with nodal selection from more than one ring around the crack tip resulted into next studies. Another type of nodal selection was considered in subsequent works by the authors [20,21,22]; the nodes were selected from specific parts of the test specimen body with specific distribution functions of their distance and angular position from the crack tip. This way was employed with expectations that the fields will be better (more accurately/efficiently) approximated. An automatic utility to determine the values of coefficients of the higher order terms of WPS using the over-deterministic method was developed [20] to enable multi-parameter description of stress field, where the coefficients of WPS terms are calculated from several layers and angular sections. And the distance distribution of the nodal selection is governed by various functions. Comparison between those used variants was given by visual technique again [21]. Hence, a new detailed way for evaluation of the accuracy of the reconstruction was used in [22]. Method based on the plot of the relative deviation (percent difference) of the stress field between the correct solution (the FE solution) and the solution given by the approximation using WPS with a certain variant of nodal selection was introduced. Main motivation of this study is to find the easiest way for the multi-parameter stress field description to obtain the sufficiently accurate solution valid for the as far as possible area from the crack tip. METHODS Multi-parameter linear elastic fracture mechanics (MP-LEFM) he stress and displacement fields in a planar homogeneous isotropic cracked body can be formulated as an infinite power expansion – Williams series [1] by Eq. (1) and (2), respectively – for more details see [14,21,22,23].        1 2, 1 ,,, , 2 n ij n ij n n Ar f n ij xy (1) T J. Sobek et alii, Frattura ed Integrità Strutturale, 39 (2017) 129-142; DOI: 10.3221/IGF-ESIS.39.14 131       2,u 1 ,,, , , n ini n uArfnE ixy (2) In this study, the attention is paid to the mode I crack problem. MP-LEFM takes into account several initial terms of WE, i.e. n ranges from 1 to N (not ∞); coefficients of these terms are often expressed as dimensionless shape functions gn (functions of the relative crack length  = a/W ). Over-deterministic method For determination of coefficients of the Williams series terms the so-called Over-deterministic method (ODM, [17]) is used. Based on the linear least-squares formulation, it solves a system of 2k equations, where k represents the number of selected nodes (in the original paper, they were selected from a nodal ring around the crack tip), for N chosen terms of the power series. Detailed analyses of this method can be found in works [18,19,24]. Displacements of selected nodes, together with their coordinates, serve as the input to that method. This issue was studied in detail in previous works which presented, among others, an implementation of this technique into an automatic numerical tool called ODeMApp [20] – the ODM analysis based on an arbitrary nodal selection and test specimen geometry variant is enabled. ReFraPro approach The ReFraPro (Reconstruction of Fracture Process [25]) is a Java application which allows an advanced determination of fracture characteristics of materials failing in a quasi-brittle manner (silicate-based materials). Estimation of the FPZ (its shape and size) is implemented by a technique developed by the authors combining MP-LEFM, classical non-linear models and plasticity approach. Reconstruction of the fracture process is made by this application generally based on the measured loading curves and basic mechanical properties of the material. For the purpose of this study, a part of this program is used which provides the reconstruction of stress field from the available shape functions (corresponding to values of coefficients of terms of the WE) for the given test geometry. This part is accompanied with a special tool (a class called percentdifference) which allows the display of the deviance (percent difference) between the approximated stress field and the exact solution (for which the FE solution is regarded) to test the accuracy of the solution with nodal selection from an area of interest around the crack tip. The deviation itself is expressed via pixmap grid as the relative difference [22]. NUMERICAL STUDY or the analysis of the accuracy of the WPS approximation, the wedge-splitting test specimen (WST) is used. Specimen loaded by eccentric tension through two steel platens with pins among which the steel wedge is impressed was developed by Linsbauer and Tscheg in [26]. Schema of the analysed test configuration is displayed in Fig. 1 left accompanied with drawings of its geometry in Fig. 1 right. Computations of the stress and displacement fields were realized in the ANSYS finite element software [27]. Crack elements (PLANE82) were utilized for the FE solution. An automatic interconnection procedure has been developed between the computational tool and the ODeMApp technique for the shape functions determination. A WST variant with two supports and specimens width W = 100 mm (Fig. 1 right) was considered. The employed FE mesh is shown in Fig. 2a with the crack-tip details, where two basic nodal selections (serving as the reference nodal selections) are depicted. The first is taken from the ring at the distance of 5 mm from the crack tip (Fig. 2b) and the second selection is taken from the ring at the distance of 0.5 mm from the crack tip (see Fig. 2c). These variants are labelled as ring 5 mm and ring 0.5 mm further in the text. Fig. 3a shows the FE mesh intended for models with a more general nodal selection. Variant from Fig. 3b labelled as con 180° (0°) represents a nodal selection governed by a constant distribution function (in the distance selection) taken from the whole area of test specimen (except of steel platens); qua 90° (45°) – a nodal selection with a quadratic distance distribution function from area of 90° angular section with the origin rotated by 45° angle (from the crack propagation direction) displayed in Fig. 3c; exp 90° (80°) – a nodal selection (from Fig. 3d with an exponential distribution function from the 90° section with initial rotation angle of 80°. Detailed description of the method for the conducted nodal selection, including several other variants, is carried out in [20,21,22]. The number of selected nodes is kept the same for each selection type, k = 49. F J. Sobek et alii, Frattura ed Integrità Strutturale, 39 (2017) 129-142; DOI: 10.3221/IGF-ESIS.39.14 132 Figure 1: Wedge splitting test with the schema of loading imposition (left); geometry of the analysed wedge splitting test specimen (right)[15]. Shape functions, representing the non-dimensional expressions of the coefficients of the higher order terms of the WE as functions of the relative crack length  , corresponding to the mentioned nodal selection variants were determined by the ODeMApp procedure and then they were used as inputs for ReFraPro application (which allows displaying the relative deviation from the “exact” solution of the stress field distribution). Subsequently, a post-processing regarding the stress fields reconstructions based on one, two and even higher order terms of the WPS was conducted. Figure 2a: Numerical model of the wedge splitting test, finite element mesh. Figure 2b: Detail of the FE mesh for the ring 5 mm variant with nodal selection (red dots at the distance 5 mm from the crack tip). Figure 2c: Detail of the FE mesh for ring 0.5 variant with nodal selection (red dots at the distance 0.5 mm from the crack tip). J. Sobek et alii, Frattura ed Integrità Strutturale, 39 (2017) 129-142; DOI: 10.3221/IGF-ESIS.39.14 133 Figure 3a: Numerical model of the WST specimen, used finite element mesh [22]. Figure 3b: Distribution o f FE nodes for con 180° (0°) variant of nodal selection [22]. Figure 3c: Distribution of FE nodes for qua 90° (45°) variant of nodal selection [22]. Figure 3d: Distribution of FE nodes for exp 90° (80°) variant of nodal selection [22]. RESULTS AND DISCUSSION he reconstruction of the fields of the principal tensile stress  1 and the crack opening stress  yy was performed and compared to the exact solution. Figs. 4 and 5 display the relative differences between the approximation of the field corresponding to chosen variant of nodal selection and the FEM solution. Only few examples of the reconstruction (based on the number of used Williams series’ terms equal to 1, 2, 4, 7 and 11 – positioned horizontally) are given for illustration. Cases with the relative crack length  = 0.5 are shown here as an example. It can be seen that low number of terms of the Williams series used provides very inaccurate approximation of the stress field (mainly for a wider area around the crack tip). Usage of only first two terms leads to a sufficient accuracy only in the very vicinity of the crack tip (where the classical and two-parameter LEFM holds). A sufficiently accurate solution can be provided by a usage of at least four terms of Williams expansion if an area of about 2 cm from the crack tip is requested. Comparison between the used variants of nodal selection shows one important fact – in general [21,22], variant con 180° (0°) (uniform selection from the whole body of the test specimen with the constant distribution function of distance selection) comes out as the best variant of nodal selection (in Fig. 4 and 5 it is emphasized by a green coloured frame). As it can be seen, the greater is the share of the red colour in the percent difference diagram the larger is the relative error from the FE solution. That is why this variant appears in next analyses as the reference one for comparison with variants of nodal selection from just one ring; moreover, from a closer distance from the crack tip. In some cases the variant qua 90° (45°) looks also promising and is comparable with the con variant. Nevertheless, this is true only for the chosen  = 0.5. Work [22] provides detailed views also at different  values. Next figures compare results of two basic variants of nodal selection from the vicinity of the crack tip with the con 180° (0°) variant of the nodal selection, i.e. from the whole body of the test specimen. Fig. 6 shows contour plots of the relative deviation of principal stress  1 for variants of nodal selection con 180° (0°), ring 5 mm and ring 0.5 mm for the WST specimen with relative crack length  = 0.5. As far as four terms of Williams expansion for approximation is used the same trend in the relative deviation can be observed. With increasing N the results of ring 0.5 mm variant appear to be very inaccurate which can be seen on the extension of the red colour area (denoting a large proportion of 40% difference from the FE solution). T J. Sobek et alii, Frattura ed Integrità Strutturale, 39 (2017) 129-142; DOI: 10.3221/IGF-ESIS.39.14 134  1 gn  = 0.5 1 2 4 7 11 con 180° (0°) qua 90° (45°) exp 90° (80°) Figure 4: Relative error of principal tensile stress  1 field approximation for chosen variants of nodal selection. crack tip J. Sobek et alii, Frattura ed Integrità Strutturale, 39 (2017) 129-142; DOI: 10.3221/IGF-ESIS.39.14 135  yy gn  = 0.5 1 2 4 7 11 con 180° (0°) qua 90° (45°) exp 90° (80°) Figure 5: Relative error of crack opening stress  yy field approximation for chosen variants of nodal selection. J. Sobek et alii, Frattura ed Integrità Strutturale, 39 (2017) 129-142; DOI: 10.3221/IGF-ESIS.39.14 136  1 gn  = 0.50 1 2 4 7 11 con 180° (0°) ring, 5 mm ring, 0.5 mm Figure 6: Relative error of principal tensile stress  1 field approximation for chosen variants of nodal selection. J. Sobek et alii, Frattura ed Integrità Strutturale, 39 (2017) 129-142; DOI: 10.3221/IGF-ESIS.39.14 137  yy gn  = 0.50 1 2 4 7 11 con 180° (0°) ring, 5 mm ring, 0.5 mm Figure 7: Relative error of crack opening stress  yy field approximation for chosen variants of nodal selection.