Microsoft Word - numero_42_art_14.docx V. Růžička et alii, Frattura ed Integrità Strutturale, 42 (2017) 128-135; DOI: 10.3221/IGF-ESIS.42.14 128 Focused on Mechanical Fatigue of Metals Over-deterministic method: The influence of rounding numbers on the accuracy of the values of Williams’ expansion terms Vladimír Růžička Brno University of Technology, Faculty of Civil Engineering, Institute of Structural Mechanics, Veveří 331/95, 602 00 Brno, Czech Republic ruzicka@musicer.net Lucie Malíková, Stanislav Seitl Academy of Sciences of the Czech Republic, Institute of Physics of Materials, v. v. i., Žižkova 22, 616 62 Brno, Czech Republic malikova.l@fce.vutbr.cz, http://orcid.org/0000-0001-5868-5717 seitl@ipm.cz, http://orcid.org/0000-0002-4953-4324 ABSTRACT. A study on the accuracy of the values of Williams’ expansion terms influenced by rounding numbers is presented. The results are presented taking into account a three-point bend single edge notched beam. Crack tip stress tensor components are 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 the coefficients of the terms of this series is performed using the least squares- based regression technique known as the over-deterministic method (ODM), for which displacements data obtained numerically in software ANSYS are taken as inputs. The values of Williams’ expansion terms based on the displacement data obtained are calculated by using various levels of rounding numbers and the results are compared and discussed. KEYWORDS. Over-deterministic method; Rounding numbers; Williams’ expansion; Fracture mechanics; Stress field. Citation: Ruzicka,V., Malikova, L., Seitl, S., Over-deterministic method: The influence of rounding numbers on the accuracy of the values of Williams’ expansion terms, Frattura ed Integrità Strutturale, 42 (2017)128-135. Received: 10.06.2017 Accepted: 16.06.2017 Published: 01.10.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 ometimes it is necessary to describe the crack-tip stress field by means of more than only one singular term (the well- known stress intensity factor, see e.g. [1-4] or two-parameter fracture mechanics, see e.g. [5-10]) of Williams’ expansion (WE). It has been shown that so-called multi-parameter fracture mechanics can be helpful when the fracture process occurs in the more extensive surroundings of the crack tip, see for material like concrete [11-19]. Then, a particular number of the WE terms needs to be taken into account for a more complex fracture analysis. In this paper, the focus is devoted to investigations into the accuracy of the higher-order terms of the WE when various types of rounding S V. Růžička et alii, Frattura ed Integrità Strutturale, 42 (2017) 128-135; DOI: 10.3221/IGF-ESIS.42.14 129 numbers are considered. There exist various kinds of software that enable to calculate higher-order terms of the WE, see [20-22]. For each of them a specific accuracy of the numbers used is defined, which can play a key role when the WE terms are calculated. The analysis presented deals with various numbers of the decimal places that are considered in the calculations of the WE terms and brings practical recommendations how the analysis should be performed in order to obtain accurate results. THEORETICAL BACKGROUND Fracture mechanics: Williams’ expansion s it has been mentioned, the paper is based on the idea that the near crack-tip stress field tensor components are approximated via Williams’ expansion [23] that was originally derived for a homogeneous elastic isotropic cracked body subjected to arbitrary remote loading and is expressed via the infinite power series for loading mode I as follows:  nfArn ij n n n ij , 21 1 2       , where i,j {x,y}. (1) The meaning of the symbols used in Eq. 1 can be described in the following way: (r, θ) are polar coordinates centred at the crack tip; are known functions corresponding to the stress distribution; the symbols An correspond to the unknown coefficients of Williams’ expansion terms (this is to emphasize that their values depend on the specimen geometry, relative crack length  and loading conditions). Over-deterministic method When the effect of rounding numbers is investigated, the over-deterministic method is assumed to be used for determination of various numbers of the WE terms [24]. This method is based on the least-squares technique and was used as one of the methods that do not require any special crack elements or implementation of other difficult fracture mechanics concepts. The method uses the displacement field estimated around the crack tip via the finite element method and together with the polar coordinates of the nodes, where the displacements are investigated, a system of linear equations is solved according to the definition:   ,,, 0 2/ EnfAru ui n n n i     , where i {x,y}. (2) When the k number of nodes is investigated, then 2k of displacements (u and v) can be used and 2k of equations can be formed in the following way [24]:                                                                                                      N kk v rNkk v rkk v r v rN v r v r v rN v r v r kk u rNkk u rkk u r u rN u r u r u rN u r u r k k A A A rfrfrf rfrfrf rfrfrf rfrfrf rfrfrf rfrfrf v v v u u u            1 0 10 22221220 11111110 10 22221220 11111110 2 1 2 1 ,,, ,,, ,,, ,,, ,,, ,,,       (3) The solution of the system of the equation can be written as:           UCCCX T1T  (4) A V. Růžička et alii, Frattura ed Integrità Strutturale, 42 (2017) 128-135; DOI: 10.3221/IGF-ESIS.42.14 130 It is always necessary to prescribe the number of the terms of the WE terms, N, that shall be calculated. This is also one of the parameters that is varied in the analysis presented. Data structure used in common software Nowadays, most personal computers used for calculation of mathematical tasks are constructed by 32 or 64 bit data architecture, see [25, 26], i.e. 4 or 8 Bytes are used for one number. For representation of numbers, mathematical software can use a structure with or without a floating point. Counting with numbers without a floating point is not usable for very low and very high numbers, especially with some decimal numbers. For counting with decimal numbers, it is necessary to use a data structure with a floating point; it gives maximally 15 decimal numbers plus an exponent in 64 bit data architecture or only 6 decimal numbers in 32 bit data architecture. For this reason, all the common software (e.g. Maple [20], Matlab [21], Mathematica [22],) uses the data structure with a floating point or a special data structure to extend the count of decimal numbers. One of the methods how to increase the count of decimal numbers is to convert the original number to a string of characters and then to carry out calculating with strings in specially programmed procedures. NUMERICAL EXAMPLE OF 3PB or the analysis a very simple model of a three-point-bending specimen was used, see Fig. 1. The detailed geometry can be found in [27, 28] (c = 0.2 mm, D = 1 mm, P = 1 N). Because of the symmetry, only one half of the specimen could be modelled in order to obtain a set of displacement data at the distance of 0.1 mm from the crack tip. Figure 1: Schema of the three-point-bending specimen used for the analysis presented. The data of the displacement vector obtained from the ANSYS finite element software [29] can be found in Tab. 1. The values introduced in the referred table were used as inputs for Eq. 3 and 4 respectively. RESULTS AND DISCUSSION he main goal was to validate the ODM concept in order to obtain a reliable procedure for further analysis of the stress field near the crack tip in civil engineering materials. Therefore, a basic cracked specimen configuration (3PB) has been investigated and higher-order terms coefficients estimated. Data comparison can be found in Tab. 2. Note that only the first five terms are mostly available in literature [27, 28, 30]. It can be seen in Tab. 2 that the coefficients calculated by means of the ODM correspond very well with the data published in literature. F T V. Růžička et alii, Frattura ed Integrità Strutturale, 42 (2017) 128-135; DOI: 10.3221/IGF-ESIS.42.14 131 Nr. Node Nr. Radius [mm] Theta [deg] u [mm] v [mm] 1 2 0.1 0 -27.077 0.0000 2 13 0.1 9 -27.066 0.0639 3 14 0.1 18 -27.034 0.1325 4 15 0.1 27 -26.984 0.2100 5 16 0.1 36 -26.919 0.3000 6 17 0.1 45 -26.843 0.4072 7 18 0.1 54 -26.765 0.5292 8 19 0.1 63 -26.688 0.6691 9 20 0.1 72 -26.618 0.8250 10 21 0.1 81 -26.564 0.9945 11 12 0.1 90 -26.529 1.1737 12 97 0.1 99 -26.517 1.3584 13 98 0.1 108 -26.534 1.5431 14 99 0.1 117 -26.581 1.7226 15 100 0.1 126 -26.659 1.8915 16 101 0.1 135 -26.768 2.0460 17 102 0.1 144 -26.906 2.1766 18 103 0.1 153 -27.069 2.2843 19 104 0.1 162 -27.253 2.3634 20 105 0.1 171 -27.449 2.4105 21 87 0.1 180 -27.651 2.4236 Table 1: The displacements and polar coordinates of the selected nodes around the crack tip obtained from the ANSYS finite element software. 3PB [27] [28] Data calculated A1 1.8530 1.8538 1.8538 A2 -0.3533 -0.3527 -0.3527 A3 -0.7149 -0.7202 -0.7202 A4 -0.0896 -0.1019 -0.1019 A5 -1.3717 -1.3565 -1.3565 Table 2: Higher-order terms coefficients determined by means of the ODM in comparison to data published in literature for 3PB, see [27, 28]. On the basis of PHP software [31] the special software tool was programmed considering various levels of rounding numbers in order to analyze how many decimal numbers are needed to obtain precise values of WE terms. The software was based on the idea of representation of numbers by means of strings. The values of the WE terms were calculated step by step assuming two changing parameters: the level of rounding numbers and the number of calculated WE terms. V. Růžička et alii, Frattura ed Integrità Strutturale, 42 (2017) 128-135; DOI: 10.3221/IGF-ESIS.42.14 132 Figure 2: Dependence of A1 corresponding to the stress intensity factor on the number of the WE terms considered for various numbers of the decimal places taken into account during the analysis. Figure 3: Dependence of the 5th WE term on the number of the WE terms considered for various numbers of the decimal places taken into account during the analysis. Several phenomena can be observed from the investigation. For example the first coefficient A1 (corresponding to the stress intensity factor K) differs significantly when 5, 6 or 10 decimal numbers are taken into account and these values are nearly independent of the number of the WE terms calculated (the minimum counted values were 5, the maximum 20), see Fig. 2. V. Růžička et alii, Frattura ed Integrità Strutturale, 42 (2017) 128-135; DOI: 10.3221/IGF-ESIS.42.14 133 Similar dependences can also be observed for other values of the WE terms of higher-orders. It holds that the level of rounding numbers is more important when the analysis for a higher index of the WE terms coefficients is performed. For example, A5 (the 5th coefficient of WE) does not seem to be precise enough until 15 decimal numbers are considered within the counting numbers (see Fig. 3); and A10 (the 10th coefficient) needs rounding by 20 decimal numbers, see Fig. 4. Figure 4: Dependence of the 10th WE term on the number of the WE terms considered for various numbers of the decimal places taken into account during the analysis. A similar trend is also expected for the coefficients of the WE terms of higher orders: the higher index of the WE term, the higher number of decimal places needed. Simultaneously it holds that more WE terms need to be considered during the analysis when the coefficient of the WE term of some higher-order is required. CONCLUSIONS pecial computing software tool based on the concept of representation of numbers by means of strings was developed which enables to calculate with numbers with up to 100 or 200 decimal places. The software was used in order to perform an extended analysis dealing with the influence of rounding numbers on the accuracy of the WE terms coefficients determined via the over-deterministic method. The results show that it is enough to use numbers with up to 20 decimal characters to obtain good results for 10 initial coefficients of Williams’ expansion. A final conclusion/recommendation based on the research presented can be stated: the higher number of the coefficients of WE terms is requested, the higher number of decimal places used within rounding numbers is needed. ACKNOWLEDGEMENT he authors acknowledge the support of Czech Sciences foundation project No.17-01589S. S T V. Růžička et alii, Frattura ed Integrità Strutturale, 42 (2017) 128-135; DOI: 10.3221/IGF-ESIS.42.14 134 REFERENCES [1] Anderson, T. L., Fracture mechanics fundamentals and applications, CRC Press (1991). [2] Pook L. P., Linear Elastic Fracture Mechanics for Engineers. Theory and Applications. WIT Press, Southampton (2000). [3] Pook, L. P. Crack Paths, WIT Press, Southampton (2002). [4] ASTM E399 - 12e3 Standard Test Method for Linear-Elastic Plane-Strain Fracture Toughness KIc of Metallic Materials, DOI: 10.1520/E0399. [5] Leevers, P.S., Radon, J.C., Inherent stress biaxiality in various fracture specimen geometries, International Journal of Fracture, 19(4) (1982) 311–325. DOI: 10.1007/BF00012486. [6] Karihaloo B.L., Abdalla H.M., Xiao Q.Z., Size effect in concrete beams, Engng. Fract. Mech., 70 (2003) 979–993. [7] Seitl, S., Veselý, V., Řoutil, L. Two-parameter fracture mechanical analysis of a near-crack-tip stress field in wedge splitting test specimens, Computers and Structures 89(21-22) (2011) 1852–1858. DOI: 10.1016/j.compstruc.2011.05.020. [8] Chao Z.J., Zhang X., Constraint effect in brittle fracture, In: Piacik et al (eds) Proceedings of the 27th National Symposium on Fatigue and Fracture, Philadelphia, (1997) 41–60. [9] Du, Z.Z., Hancock, J.W., The effect of non-singular stresses on crack tip constraint, J. Mech. Phys. Solids, 39 (1991) 555–567. [10] Dyskin, A.V., Crack growth criteria incorporating non-singular stresses: size effect in apparent fracture toughness, Int. J. Fract., 83 (1997) 191–206. [11] Ayatollahi, M.R., Akbardoost, J., Size effects on fracture toughness of quasi-brittle materials – A new approach, Engng. Fract. Mech., 92 (2012) 89–100. [12] Berto, F., Lazzarin, P., On higher order terms in the crack tip stress field, Int. J. Fract., 161 (2010) 221–226. [13] Vesely, V., Sobek, J., Sestakova, L., Frantik, P., Seitl S. Multi-parameter crack tip stress state description for estimation of fracture process zone extent in silicate composite WST specimens, Frattura ed Integrita Strutturale, 7(25) (2013) 69– 78. [14] Veselý, V., Frantík, P., An application for the fracture characterisation of quasi-brittle materials taking into account fracture process zone influence, Advances in Engineering Software, 72 (2014) 66–76. [15] Veselý, V., Frantík, P., Sobek, J., Malíková, L., Seitl, S. Multi-parameter crack tip stress state description for evaluation of nonlinear zone width in silicate composite specimens in component splitting/bending test geometry, Fatigue and Fracture of Engineering Materials and Structures, 38(2) (2015) 200–214 [16] Veselý, V., Sobek, J., Frantík, P., Seitl, S., Multi-parameter approximation of the stress field in a cracked body in the more distant surroundings of the crack tip, International Journal of Fatigue, 89 (2015) 20–35 [17] Stepanova, L., Roslyakov, P., Complete Williams Asymptotic Expansion Near the Crack tips of Collinear Cracks of Equal Lengths in an Infinite, Procedia Structural Integrity, 2 (2016) 1789–1796, DOI: 10.1016/j.prostr.2016.06.225 [18] Stepanova, L., Roslyakov, P., Gerasimova, T., Complete Williams Asymptotic Expansion near the Crack Tips of Collinear Cracks of Equal Lengths in an Infinite Plan", Solid State Phenomena, 258 (2017) 209–212. DOI: 10.4028/www.scientific.net/SSP.258.209 [19] Stepanova, L., Roslyakov, P., Multi-parameter description of the crack-tip stress field: Analytic determination of coefficients of crack-tips stress expansions in the vicinity of the crack tips of two finite in an infinite plane medium, International Journal of Solids and Structutes, 100-101 (2016) 11–28. [20] http://www.maplesoft.com/ [21] https://www.mathworks.com/products/matlab.html [22] https://www.wolfram.com/mathematica/ [23] Williams, M.L., On the stress distribution at the base of a stationary crack, J. Appl. Mech., 24 (1957) 109–114 [24] Ayatollahi M, Nejati M. An over‐deterministic method for calculation of coefficients of crack tip asymptotic field from finite element analysis. Fatigue Fract Eng Mater Struct, 34 (2011) 159–76. [25] Cody, J.W., et al. IEEE standards 754 and 854 for Floating-Point Arithmetic: for a readable, account IEEE Magazine MICRO, (1984) 84–100. [26] Kahan, W. IEEE Standard 754 for Binary Floating-Point Arithmetic, Lecture Notes on the Status of IEEE, 754 (1997) 1–30. V. Růžička et alii, Frattura ed Integrità Strutturale, 42 (2017) 128-135; DOI: 10.3221/IGF-ESIS.42.14 135 [27] Karihaloo, B.L., Xiao, Q.Z., Higher order terms of the crack tip asymptotic field for a notched three-point bend beam, Int. J. Fract., 112 (2001) 111–128. [28] Šestáková, L., Tuning of an Over-Deterministic Method for Calculation of Higher-Order Terms coefficients of the Williams Expansion for Basic Cracked Specimen Configurations, Applied Mechanics 2011 (on CD) [29] ANSYS Program Documentation, User’s manual version 10.0. Swanson Analysis System, Inc., Houston (2005). [30] Karihaloo, B.L. Xiao, Q.Z., Accurate determination of the coefficients of elastic crack tip asymptotic field by a Hybrid crack element with p-adaptivity, Engineering Fracture Mechanics, 68 (2001) 1609–1630 [31] http://php.net/software.php << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice