Fracture, Fatigue, Failure and Damage Evolution , Volume 3

River Rapids Conference Proceedings of the Society for Experimental Mechanics Series Fracture, Fatigue, Failure and Damage Evolution, Volume 3 Shuman Xia Allison Beese Ryan B. Berke Proceedings of the 2020 Annual Conference on Experimental and Applied Mechanics River Publishers

Conference Proceedings of the Society for Experimental Mechanics Series Series Editor Kristin B. Zimmerman, Ph.D. Society for Experimental Mechanics, Inc., Bethel, CT, USA

The Conference Proceedings of the Society for Experimental Mechanics Series presents early findings and case studies from a wide range of fundamental and applied work across the broad range of fields that comprise Experimental Mechanics. Series volumes follow the principle tracks or focus topics featured in each of the Society’s two annual conferences: IMAC, A Conference and Exposition on Structural Dynamics, and the Society’s Annual Conference & Exposition and will address critical areas of interest to researchers and design engineers working in all areas of Structural Dynamics, Solid Mechanics and Materials Research.

River Publishers Fracture, Fatigue, Failure and Damage Evolution, Volume 3 Proceedings of the 2020 Annual Conference on Experimental and Applied Mechanics Shuman Xia • Allison Beese • Ryan B. Berke Editors

Published, sold and distributed by: River Publishers Broagervej 10 9260 Gistrup Denmark www.riverpublishers.com ISBN 978-87-4380-006-4 (eBook) Conference Proceedings of the Society for Experimental Mechanics An imprint of River Publishers © The Society for Experimental Mechanics, Inc. 2021 This work is subject to copyright. All rights are solely and exclusively licensed by the Publisher, whether the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of illustrations, recitation, broadcasting, or reproduction in any other physical way, and transmission or information storage and retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology now known or hereafter developed. The use of general descriptive names, registered names, trademarks, service marks, etc. in this publication does not imply, even in the absence of a specific statement, that such names are exempt from the relevant protective laws and regulations and therefore free for general use. The publisher, the authors, and the editors are safe to assume that the advice and information in this book are believed to be true and accurate at the date of publication. Neither the publisher nor the authors or the editors give a warranty, expressed or implied, with respect to the material contained herein or for any errors or omissions that may have been made. The publisher remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Preface Fracture, Fatigue, Failure and Damage Evolution represents one of seven volumes of technical papers to be presented at the SEM 2020 SEM Annual Conference & Exposition on Experimental and Applied Mechanics organized by the Society for Experimental Mechanics scheduled to be held in Orlando, FL, September 14–17, 2020. The complete proceedings also include volumes on: Dynamic Behavior of Materials; Mechanics of Time Dependent Materials; Advancement of Optical Methods & Digital Image Correlation in Experimental Mechanics; Mechanics of Biological Systems and Materials, Micro- and Nanomechanics & Research Applications; Mechanics of Composite, Hybrid & Multifunctional Materials; and Thermomechanics & Infrared Imaging, Inverse Problem Methodologies and Mechanics of Additive & Advanced Manufactured Materials. Each collection presents early findings from experimental and computational investigations on an important area within experimental mechanics, fracture and fatigue being one of these areas. Fatigue and fracture are two of the most critical considerations in engineering design. Understanding and characterizing fatigue and fracture has remained as one of the primary focus areas of experimental mechanics for several decades. Advances in experimental techniques, such as digital image correlation, acoustic emissions, and electron microscopy, have allowed for deeper study of phenomena related to fatigue and fracture. This volume contains the results of investigations of several aspects of fatigue and fracture such as microstructural effects, the behavior of interfaces, the behavior of different and/or complex materials such as composites, and environmental and loading effects. The collection of experimental mechanics research included here represents another step toward solving the long-term challenges associated with fatigue and fracture. Atlanta, GA, USA ShumanXia State College, PA, USA Allison Beese Logan, UT, USA Ryan B. Berke v

Contents 1 Influence of Fracture and Delayed Effects on Steel-Concrete Composite Structures............................ 1 Fery Léo, Gwozdziewicz Piotr, and Moutou Pitti Rostand 2 Characterization of High Frequency Pulse Loading on Fatigue of Metals........................................ 7 Paul A. Lara, Hugh A. Bruck, and Edda C. Müller 3 Fracture Parameters and Failure Visualization of Al6063-T6 Under Different Loading Rates ................ 19 Anoop Kumar Pandouria, Purnashis Chakraborty, Sanjay Kumar, and Vikrant Tiwari 4 Fatigue Life Prediction of Natural Rubber in Antivibratory Applications ....................................... 29 Benoit Ruellan, Jean-Benoit Le Cam, Isabelle Jeanneau, and Frédéric Canévet 5 Fatigue Assessment of Porosity in Electron Beam Melted Ti-6Al-4V.............................................. 37 Justin Warner, Dino Celli, Jacob Rindler, M. Herman Shen, Onome Scott-Emuakpor, and Tommy George 6 Bayesian Updating of a Cracking Model for Reinforced Concrete Structures Subjected to Static and Cyclic Loadings..................................................................................................... 45 Henriette M. Imounga, Emilio Bastidas-Arteaga, Moutou Pitti Rostand, and Serge Ekomy Ango 7 Crack Jumping in Fabric Composite Fracture Testing ............................................................. 51 Brian T. Werner, Vincente Pericoli, and James W. Foulk III 8 Effect of Crack-Parallel Compression or Tension on Mode-I Fracture Energy of Quasibrittle Material – As Applied to Concrete .................................................................................... 55 Hoang Thai Nguyen, Masoud Rezaei, Madura Pathirage, Gianluca Cusatis, Mohsen Issa, and Zdeneˇk P. Bažant 9 Modal Validation of Academic Bladed Disk with DIC.............................................................. 61 Dino Celli, Onome Scott-Emuakpor, and Tommy George 10 Assessing Bond Strength in 304L Stainless Steel Plate Welded Using Plastic Explosives ....................... 65 Thomas A. Ivanoff, Olivia D. Underwood, Jonathan D. Madison, Lisa A. Deibler, and Jeffrey M. Rodelas 11 Real-Time Visualization of Damage Progression Inside GFRCs via High-Speed X-Ray PCI Technique ..... 69 Jinling Gao, Nesredin Kedir, Cody Kirk, Julio Hernandez, Xuedong Zhai, Junyu Wang, Tyler Tallman, Kamel Fezzaa, and Weinong Chen 12 Watching High-Cycle Fatigue with Automated Scanning Electron Microscope Experiments.................. 73 Nathan M. Heckman, Timothy A. Furnish, Christopher M. Barr, Khalid Hattar, and Brad L. Boyce 13 Determination of Mixed-Mode (I/III) Fracture of Polycarbonate ................................................. 77 Ali F. Fahem, Vijendra Gupta, Addis Kidane, and Michael A. Sutton 14 Influence of Dynamic Multiaxial Transverse Loading on Dyneema® SK76 Single Fiber Failure .............. 85 Frank David Thomas, Stephen L. Alexander, C. Allan Gunnarsson, Tusit Weerasooriya, and Subramani Sockalingam vii

viii Contents 15 Using Crack Geometry to Determine Fracture Properties ......................................................... 93 Kimberley A. Mac Donald and Guruswami Ravichandran 16 Dynamic Fracture-Toughness Testing of a Hierarchically Nano-Structured Solid .............................. 97 Kyung-Suk Kim, Hanxun Jin, Tong Jiao, and Rodney J. Clifton

Chapter 1 Influence of Fracture and Delayed Effects on Steel-Concrete Composite Structures Fery Léo, Gwozdziewicz Piotr, and Moutou Pitti Rostand Abstract This work investigates the fracture behavior of preflex beams submitted to time-dependent effects and long-term loadings. However, it is necessary to consider delayed effects such as shrinkage and creep on these structures because their influence is not negligible. The different materials used in mixed construction are presented. The flowchart applied to compute creep and relaxation effects versus the time is also described. It is interesting to note that the analytical results obtained for a preflex beam are very similar to those obtained using finite element modeling on the LUSAS software. It is also observable that the cracking of the concrete that encases the lower fiber, due to creep, causes an increase of stresses in the steel profile, which can break the beam. The use of reinforcement on the preflex beam decreases the crack opening after 100 years loadings. Keywords Steel-concrete composite structure · Creep · Shrinkage · Fracture · Preflex beam 1.1 Introduction Mixed construction is a method which knows a quick development and takes advantage of two materials: steel and concrete [1, 2]. They are linked together by connectors, so that concrete works in compression and steel in tension. This results in a composite beam with improved mechanical characteristics because both concretes are supposed to work in compression. However, it is necessary to consider the impact of delayed effects such as shrinkage and creep on the fracture processes of these structures because their influence is not negligible. The literature background shows that some authors discussed about the behavior and the manufacturing process of preflex beams [3, 4]. Foort [5] studied the behavior of connectors in preflex beams due to normal and rupture force. Jerome [6] investigated the cracking of stressed preflex beams by comparing numerical and analytical results given by the Ghali and Fabre method [7]. But the impact of delayed effects as creep and relaxation of crack process have not been investigated in this chapter. The different materials used in mixed construction are presented. At the same time, the conception step of the preflex beam and the numerical flowchart to obtain crack opening are described. It is also observable that the cracking of the concrete that encases the lower fiber, due to creep, causes an increase of stresses in the steel profile, which can lead to the total collapse of the beam. F. Léo Université Clermont Auvergne, CNRS, Institut Pascal, Clermont-Ferrand, France G. Piotr Facukty of Civil Engineering, Politechnika Krakowska, Warszawska, Krakow, Poland e-mail: pgwozdziewicz@pk.edu.pl M. P. Rostand ( ) Université Clermont Auvergne, CNRS, Institut Pascal, Clermont-Ferrand, France Institut de Recherche Technologique, CENAREST, Libreville, Gabon e-mail: rostand.moutou_pitti@uca.fr © The Society for Experimental Mechanics, Inc. 2021 S. Xia et al. (eds.), Fracture, Fatigue, Failure and Damage Evolution, Volume 3, Conference Proceedings of the Society for Experimental Mechanics Series, https://doi.org/10.1007/978-3-030-60959-7_1 1

2 F. Léo et al. 1.2 Material and Methods 1.2.1 Materials: Preflex Beam The preflex beam is studied in the rest of this chapter. An isostatic beam with a span of 17 m was considered subject to a variable operating load, Fig. 1.1a. The section consists of an HEM 700, first phase concrete B1 (in red) and second phase concrete B2 (in blue), as can be seen in Fig. 1.1b. All dimensions are given here in meters. According to Table 1.1a, a shade of 275 MPa has been chosen. The elastic moment of the profile is 2299 kN.m. Concrete B1 and B2 are C50/60 and C25/30, respectively, as posted in Table 1.1b. For numerical computation, a load of 7 kN/m2 applied to the entire upper surface of concrete B2 is applied. The relative humidity of the environment has been assumed to be 70%. 1.2.2 Modelling Steps The different steps of the building are posted in Table 1.2. Concrete B1 is poured on day 1 and concrete B2 on day 31. In fact, we wait 28 days after pouring B1 before releasing the prebending forces and 28 days after pouring B2 before commission the beam. The prebending force is composed of two point forces: one located at the quarter of the beam and another located at three-quarters. It is necessary that the profile does not exceed its elastic limit during the prebending step. In order to calculate the rupture load of the preflex beam in both cases (with and without the delayed effects), we will increase the operating load until the rupture stress of the profile is reached. The finite element software LUSAS is used in order to implement the numerical model. An important assumptionwhich consists in supposing that concrete B1 does not work anymore and can be removed when it is completely cracked. It should be noted that concrete B1 cracks because the loading is important. In general, this concrete is used in preflex beams for its toughness in compression. When the complete cracking appears, we denote two load cases: one before the total cracking and one after. For the second load case (after cracking), we remove concrete B1 and add a linear force representing the self-weight of B1 before restarting the calculation. Thus, for each load case (the load is gradually increased), we proceeded as detailed in the flowchart in Fig. 1.2, which is as follow: Fig. 1.1 (a) Static model of the beam. (b) Section of beam studied

1 Influence of Fracture and Delayed Effects on Steel-Concrete Composite Structures 3 Table 1.1 (a) Characteristic of the profile HEM 700. (b) Characteristics of the concrete in B1/B2 phases (a) Mechanical characteristic of steel HEM700 A= 383 cm2 h= 716 mm b= 304 mm tf = 40 mm tw = 21 mm y cdg= 358 mm Avz = 169,8 cm2 ly = 329278 cm4 Wel,y = 9198 cm3 Mel = 2299 kN.m fyk = 275 MPa Ea = 210000 MPa (b) Characteristic of concrete Concrete B1 Concrete B2 Resistance class C50/60 Resistance class C25/30 Class of cement ClassN Class of cement ClassN AcB1= 184000 mm2 Ac B2= 576000 mm2 uB1 = 2060 mm uB2= 5280 mm ts B1= 1 jour ts B2= 1 jour fck B1= 50 MPa fck B2= 25 MPa fctm B1 = 4,1 MPa fctm B2= 2,6 MPa Ecm B1= 37000 MPa Ecm B2= 31000 MPa Table 1.2 Preflex beam manufacturing steps Stapes Tasks Days 1 Profile preflexion 0 2 Consideration of delay effects of concrete B1 29 3 Relaxation of the preflexion loads 30 4 Consideration of delay effects of concrete B2 59 5 Loading 60 6 Consideration of delay effects at 100 years 36,500 • Step 1. First, we list the lower and upper fiber stresses in concrete B1 at 100 years. • Step 2. It is then checked whether these values exceed the admissible stress. • Step 3. If this is the case, we are looking for the date on which concrete B1 is completely cracked. If this is not the case, we go directly to step 6. • Step 4. Two load cases are created: one before cracking and one after the removal of concrete B1 by adding its own weight. • Step 5. We relaunch the calculations and we list the stresses suffered by concrete B2 at 100 years because removing concrete B1 increases the stresses in the profile and in concrete B2. • Step 6. We check that concrete B2 does not reach its admissible limit, and we proceed the same way as for B1 if this is the case. Finally, we list the maximum stress in the steel profile and check that it remains acceptable.

4 F. Léo et al. Stresses to concrete B1 at 100years Stresses in B1 and B2 > valid limit If No Yes One case stress before and after Stresses to concrete B2 at 100years Relaunch calculations Stresses in B1 and B2 > valid limit If No Yes Stresses to profile Step1 Step2 Step3 Step4 Step5 Step6 Fig. 1.2 Flowchart for computing the rupture stress 1.3 Results and Discussion 1.3.1 Crack Opening with First Reinforcement The fracture due to creep and relaxation effects have been discussed in this section. The applying force is 22 kN/m2 (breaking load of the beam when we do not use reinforcing steel taking up traction in concrete) at the serviceability limit state. The consideration of cracking of concrete under LUSAS is computed by using steel reinforcements. We therefore started by applying a very weak reinforcement, composed of 6 mm bars every 30 cm, in the longitudinal and transverse direction of the beam, on the upper face and on the lower face. Figure 1.3 presents the relationship between crack opening versus time for creep and shrinkage effects. It was observed that crack opening due to creep is significantly greater than that due to shrinkage. Indeed, at 100 years, creep represents 77.6% of cracking after loading and shrinkage only 22.4%. On the graph, we can notice that even at 100 years, the cracking induced by shrinkage continues to increase very slightly whereas it is completely stopped for creep. 1.3.2 Crack Opening with Second Reinforcement In order to decrease the impact of delay deformations on crack opening, the reinforcement has been increased. In this case, bars of 6 mm every 15 cm were placed on the upper side and bars of 12 mm every 15 cm on the lower side. Figure 1.4 presents the relationship between crack opening versus time for creep and shrinkage effects. This reinforcement of steel allows to remain below the admissible limit of cracking, because when we take into account cracking due to loading and delayed effects, we arrive at 0.08 mm at 100 years. Once again, we realize that the cracking due to creep is much greater than that due to shrinkage. The proportion remains the same since 77.6% of the cracking due to delayed effects comes from creep.

1 Influence of Fracture and Delayed Effects on Steel-Concrete Composite Structures 5 Crack opening (mm) 0 0,1 0,2 0,3 0,4 0,5 0,6 100 1000 10000 100000 Time (days) Creep Shrinkage Fig. 1.3 Crack opening due to creep/relaxation effects with first reinforcements Crack opening (mm) 0 0,01 0,02 0,03 0,04 0,05 0,06 0,07 0,08 0,09 100 1000 10000 100000 Creep Shrinkage Time (days) Fig. 1.4 Crack opening due to creep/relaxation with second reinforcements 1.4 Conclusion In this chapter, the behavior of preflex beam is investigated. The different materials used in the design of this structure are presented. The design step and the numerical implementation of crack opening with LUSAS software are detailed. The crack processes due to creep and shrinkage effects are investigated. It is observed that the use of reinforcement decreases significantly the crack opening even if the creep is responsible for more than 77% of the crack opening process. From these

6 F. Léo et al. results, it can be concluded that the consideration of delayed effects leads to a lower failure load of the beam. Therefore, it is necessary to consider the delayed effects in the dimensioning of a preflex beam in order to avoid the total collapse of the structure. References 1. Aribert, J.M.: Construction mixte acier-béton: Calcul des poutres mixtes. Techniques de l’ingénieur. Novembre 2004 2. Aribert, J.M.: Analyse et formulation pratique de l’influence de la nuance de l’acier du profilé sur le degré minimum de connexion partielle d’une poutre mixte. Revue construction métallique n◦3 (1997) 3. Bae, D., Lee, K.: Behavior of preflex beam in manufacturing process. KSCE J. Civ. Eng. 8, 111–115 (2004). https://doi.org/10.1007/ BF02829086 4. Foort, E.: Connecteurs dans les poutres Preflex: comportements dus à l’effort normal et à l’effort de rupture – Mémoire de fin d’étude, Polytech Clermont Ferrand. Politechnika Krakowska (2016) 5. Ghali, A., Fabre, R.: Concrete Structures: Stresses and Deformations, 2nd edn. Spon, London/New York (1994) 6. Jerome, M.: Étude de la fissuration des poutres Preflex sollicitées. Comparaions numérique et analytique. Mémoire de fin d’étude, Polytech Clermont Ferrand, Politechnika Krakowska – 2011 7. Zhang, K., Li, S., Liu, K.: Experimental Study on Static and Fatigue Behavior of Steel-Concrete Preflex Prestressed Composite Beam. Advances in Steel Structures (ICASS ‘99). Proceedings of The Second International Conference on Advances in Steel Structures 15–17 December 1999, Hong Kong, China Volume II, pp. 965–973 (1999)

Chapter 2 Characterization of High Frequency Pulse Loading on Fatigue ofMetals Paul A. Lara, Hugh A. Bruck, and Edda C. Müller Abstract Aluminum materials of various grades are utilized across many industries, spanning from the cycling, automotive, aerospace, and the marine industry. In the marine grade, aluminum materials are utilized to construct entire vessels of various lengths or portions of them by taking advantage of the lower weight characteristics of the materials and impact on stability of the structures. In particular, 5xxx series aluminum materials are relied on by the marine industry for these purposes, taking advantage of the ability of these series to resist marine corrosive environments. However, during its lifetime, a marine vessel will experience a multitude of variable amplitude loading conditions, with occasional overloads and underloads depending on the operation and environmental conditions. In some cases, these overloads/underloads can result in the failure of the structures by reaching its ultimate capacity, but in other instances, they can systematically affect the growth rate of localized cracks. Existing models, like the Wheeler, Willenborg, plus variations of these, have been utilized to predict the crack growth behavior with varying degrees of success. We created an experimental matrix to explore the effects of overload/underload combinations on fatigue crack growth in 5xxx aluminum. Both visual inspection of crack tip location and digital image correlation (DIC) characterization of the crack tip deformation fields were used to characterize the crack growth in center crack tension (CCT) panel specimens. DIC also enabled additional analysis of strain fields to elucidate on the conditions responsible for change in the crack growth behavior. This chapter outlines some of the ongoing results of this work, which built on past experimental work conducted. Future phases of this work will utilize this data to develop new models for fatigue crack growth, and application of multiple pulses sequences. Keywords Fatigue crack growth · High frequency pulse · Overloads/underloads · Plastic zone · Crack kinking 2.1 Introduction Embedded high frequency (HF) signal effects derived from wave impacts on ships can affect failure mechanisms on the structures and have an adverse impact on the fatigue life of the vessel. While operating in a sea environment, ship structures can be subject to many operational loads (wind, pressure, temperature, etc.), one of which is the structural effects derived from the surrounding sea environment. Typically, the wave environment applies an ordinary wave component which drives the primary bending stress of the vessel, along with a more stochastically driven element that manifest itself as wave impacts. This dynamic wave impacts results in a high frequency vibratory response signals load applied to ship structures. The vibratory nature of these high frequency responses makes the design of structures increasingly complex in nature due to their uncertainty; designers and naval design rule societies have relied on methods such as safety factors and/or margins of safety to account for its effects. A typical wave impact on a naval structure imparts a time-dependent pulse load that exhibits a higher frequency logarithmic decaying sinusoidal response when measured via experimental means and captured utilizing P. A. Lara Naval Surface Warfare Center, Carderock Division, West Bethesda, MD, USA e-mail: lara@umd.edu H. A. Bruck ( ) Mechanical Engineering Department, University of Maryland, College Park, MD, USA e-mail: bruck@eng.umd.edu E. C. Müller Düsseldorf University of Applied Science, Düsseldorf, Germany e-mail: edda.mueller@study.hs-duesseldorf.de © The Society for Experimental Mechanics, Inc. 2021 S. Xia et al. (eds.), Fracture, Fatigue, Failure and Damage Evolution, Volume 3, Conference Proceedings of the Society for Experimental Mechanics Series, https://doi.org/10.1007/978-3-030-60959-7_2 7

8 P. A. Lara et al. data acquisition systems. The figure below displays one of these signals, where the measured signal can be decomposed into its different components by applying filtering techniques. Under the current practices, the peaks of these responses are utilized to compare to the quasi-static design stresses as well as evaluate the fatigue life of the structures either under constant or variable amplitude experimentation. Figure 2.1 below displays this kind of behavior, where underloads and underloads affect the crack growth rate behavior of material under a simplified signal approach [1]. But there are additional higher frequency characteristics to these high frequency signals that are not currently taken into consideration. Existing academic research has been centered on capturing a simplified sinusoidal response associated with this slamming event and embedded high frequency response but has not addressed logarithmic decay, signal frequency, or frequency of occurrence. Work by Sumi et al. has centered on more complex signals, where bounding functions were used to predict crack growth behavior [2], and the results show an influence of crack length predictions based on utilizing these bounding functions. All these factors have associated uncertainty and impact fatigue life and failure mechanisms exhibited by ship structures; work by researchers such as Fricke et al. has shown damage variation related to the amplitude of high frequency signals [3]. These high frequency effects can cause retardation of crack growth, and foundational models such as Wheler [4] and Willenborg [5] were developed to capture the overall effects of these overloads and underloads. Follow-on work was carried out to gain better understanding of the conditions present, as shown in Fig. 2.2, where the crack growth effects within the overload plastic zone boundary exhibits a variable crack growth rate behavior instead of constant crack growth retardation effect. Fig. 2.1 Load sequencing [2] Fig. 2.2 Plastic zone model and retardation effects, left [6], right [7]

2 Characterization of High Frequency Pulse Loading on Fatigue of Metals 9 Existing literature gaps are: a deeper understanding of resulting fundamental inelastic aluminum behavior for center crack tension (CCT) specimens due to HF pulses elucidate the plastic zone behavior and resulting material failure modes by leveraging full strain field displacement measurements from digital image correlation (DIC) and develop failure models that accounts for observed material failure behavior due to HF pulses and overall fatigue life response. In the future, the foundational findings of this research can then be leveraged in investigative efforts to develop analytical models addressing heat-affected zones, provide underpinnings for high fidelity numerical modelling, eventually providing broader toolsets to designers that can be used to reduce the dependency on safety factors and introduce more rigorous failure mechanism design criterion. The work outlined here builds on previous work carried out by the authors; it summarizes some previous work, presents new findings, highlights future efforts, and in the end has an overall objective to gather fundamental understanding of the effects of this high frequency loading on Aluminum 5xxx material, accounting for some of the signal’s characteristics, and through an experimental evaluation assesses its impact on the local failure mechanism and life cycle models. 2.2 Description of Experiments This section outlines a summary of the experimental setup; the work builds from and complements existing work developed by Hart and Bruck [8], where the investigators looked at failure of aluminum 5052-H32, the applications of composite patches to aluminum, and the effects of crack growth. Whereas their work centered on the utilization of a composite laminate to mitigate crack growth on ultimate failure, this work looks at the effects of crack growth based on loading condition. Figure 2.3 below displays the aluminum material properties that Hart and Bruck obtained utilizing ATSM 8 [9] and American Society for Testing and Materials (ASTM) E466 [10], as well as the layout of their center crack panel specimens. The Modulus listed for the material is 70.9 MPa, the yield at 2% offset is 116 MPa, and the parameters αand n at 9.9 and 0.22, respectively. Initial experiments were carried out on this setup for this work, utilizing a two-post Bose ElectroForce frame which has the capacity of applying loadings up to 100 Hz and 3000 N. This device has a built-in displacement meter with errors less than 5 μm, and it has a block waveform generation software that allows for the implementation of unique load signals. Fast Fourier Transforms methods were utilized to ensure that the programmed signal and the applied signals displayed the required frequencies. But the initial results indicated that there was an imbalance of crack growth between the left and right cracks growing from the notch during fatigue loading, the reason as traced to the use of mechanical wedge friction grips. It is important to mention that the experiments by Hart and Bruck utilized a 4-post MTS servo hydraulic grip setup Fig. 2.3 5052-H32 material properties and specimen failure [5]

10 P. A. Lara et al. Fig. 2.4 Experiment setup since they were interested in the ultimate loading capacity. To address this crack asymmetry from loading imbalance, a specimen configuration was developed, and it is an amalgamation of several ASTM methods; this was conducted to address the significant asymmetry observed with crack growth, as previously described. ASTM E8 [9] guidance was utilized to develop the specimen tapers toward test region, while ASTM E647 [11] middle tension specimen configuration was utilized to size the center notch and test region, maintaining the configuration carried out by Hart and Bruck. Since a clevis pin was utilized, ASTM E338 [12] was utilized to develop end reinforcements and alignment tabs for the end pin connections to ensure the failure would not occur in or adjacent to the clevis pin area. Electrical discharge machining (EDM) techniques were utilized to manufacture the specimens, providing a high level of accuracy for the component geometry, since it was also critical to ensure alignment of end pins and size and location of center notch. The material for the specimen is marine grade aluminum 5052-H32, with a nominal thickness of 2.29 mm (0.09in). Figure 2.4 shows how the clevis pins were utilized in the loading train and end adaptors were used to interface with the Bose frame and to help with specimen load alignment; the Universal Grip LLC. components were utilized for this setup: Adapter - Male Bolt 5/8 –18 w/1 length, Universal Joint Adapter; in addition, Mark-10 G1090 high-capacity clevis grip (pair) were utilized. The load capacity of the members was at a minimal load of 9KN, hence far exceeding the capacity of the Bose equipment of 3KN, and ensuring no failure could be accidentally induced to the Bose fixture from a failure of the adapters under load. The experimental setup consisted of several components’ setup to provide dual measurements: one was to observe/capture the crack growth on the back of the specimens using a microscope and another to capture digital image correlation information for the future extraction of surface strain measurements. Figure 2.4 details the setup utilized for these experiments. In the back face of the specimens, an optics zoom inspection microscope head by Edmunds Optics was used along with a Dolan-Jenner MI-150 fiber optic illuminator with a fiber optic ring light accessory. This was utilized to capture the crack growth on the back of the specimen. The setup was connected with a rack and pinion positioning system, which allowed movement of the microscope along three axes to allow consistent positioning of the microscope in different locations to capture crack growth. The microscope was able to visualize square areas of 0.5 mm on each side. Several digital pictures were taken on each crack emanating from each notch, which overlapped, and then were stitched together to form an overall view of the crack of the specimen. In the front face of the specimen, a Point Gray Flea FL3-U3-88S2C-C digital camera attached to a Tamron 28–80 mm f/3.5–22 lens and an adapter was used to capture static pictures of the surface of the specimen. In addition, two 250 Watts halogen Lowel Pro Lights were utilized to provide lighting to the specimen. The lights set on flood (wide) mode were not directly aimed at the specimens, but instead aimed at locations adjacent to the mid-point of the specimen, and focused as to provide evenly dispersed light flooding on the specimen. Digital image correlation (DIC) was utilized on the front face of the specimen to capture field surface strains [13]; in order to capture these images, the area was prepped to generate distinct patterns that can be captured with a Point Grey FlyCap 2 V2.11.3 and processed utilizing VIC-2D Ver6.0.6 by Correlated Solutions (Lexington, SC).

2 Characterization of High Frequency Pulse Loading on Fatigue of Metals 11 Figure 2.5 shows the typical DIC pattern applied; the generation of this pattern is achieved by first applying a coat of flat white paint as the background, and the black speckle patterns rely on lightweight particles that slowly settle over the specimen. Best practices developed throughout several references were utilized to obtain acceptable DIC results, once such body of work was developed by the international digital image correlation society [14]. 2.3 Test Matrix for HF Pulse Loading Conditions In order to explore the variability of the high frequency pulses, a couple of items were developed: the first was a consistent process/protocol that was repeatable across experiments, and the second was an experimental matrix that varied certain parameters to determine their influence on the response. Figure 2.6 below outlines the top-level experimental process utilized. Specimen preparation was the initial step in the configuration; the back of the specimens that faced the microscope was polished to a mirror-type finish; this aided in the identification of the crack tip. The front of the specimens received a DIC spackle pattern that had a sigma level equal or below 0.001 [13]. The sigma level was tested during the baseline measurements stages until the desired value was achieved. The baseline measurements also utilized DIC to gather initial measurement errors bounds; if the load frame indicated a 0.2 mm displacement, did the DIC measurements also show a similar value? This was done by only pining the stationary end of the specimen and comparing the displacement between the load frame and DIC calculations. Additional baseline measurements included static DIC measurements at 0, 200, 1100, and 2000 N, as well as microscope crack measurements at no load and 2000 N. The initial crack growth was a step used to develop a crack and go outside the envelope of the stress-strain effects of the notch radius [15]. The notch radius was 0.17 mm, and an initial starter crack was developed to be around 1 mm in length. The HF pulse was applied at about 40 K cycles when the desired crack length was reached; subsequent measurements monitored the crack growth rates after the HF pulse was applied to evaluate the recovery. It is important to note that the load ranged from 200 N to 2000 N, where the stress ratio R=0.1, the frequency used for fatigue was 5 Hz, and the nominal stress intensity K=6.07MPa-m0.5. The experimental matrix utilized is shown in Table 2.1; it consists of a baseline test where no HF pulse is applied, and a 3×3 matrix with three levels and two factors (amplitude, frequency), resulting in 9 treatments, where the amplitude of the pulse is increased in 1/3 amplitude increments and frequencies that range between 4 and 10 times the baseline frequency, i.e., specimen 3C had a peak amplitude of 2900 N and an HF pulse with a frequency of 50 Hz. Fig. 2.5 Specimen and speckle pattern Fig. 2.6 Experimental process

12 P. A. Lara et al. Table 2.1 Test matrix for HF pulse loading conditions Fig. 2.7 Crack length for varying load (left) and frequency (Right) For this matrix, the HF pulse was only applied at the peak of the ordinary wave of 5 Hz; hence the phase shift was held constant as well as the log decrement of the function was maintained at 0.21. Introduction of other factors would increase the matrix size to a 33 design resulting in 27 treatments and are recommended for future experiments. 2.4 Experimental Results When plotting the data, as shown in Fig. 2.7, the initial results were the crack length for specimens where the load is varied while the frequency is held constant (left) and specimens where the frequency is held constant and the load is varied (right). One can see that when comparing the number of cycles it takes to grow a crack to a nominal length of 1.5 mm, the varying load (left) has a noticeable effect on the crack growth. It takes 20–30 K cycles in delay depending on the load. One can see that the number of cycles it takes for the crack to grow 0.5 mm on the baseline is substantially smaller that the cycles needed for the specimens that had an HF pulse applied. When looking at the varying frequency (right), one can see that a lower frequency has a higher effect for the same max load of 2300 N. When looking diagonally across the matrix, one can see in Fig. 2.7 that the lowest block with the lowest load and frequency behaves very similarly to the baseline, but other blocks begin to deviate as the load and frequency of the HF pulse are increased. Vic2D DIC was utilized to evaluate the strain fields ahead of the crack before and after the HF pulse was applied. Figure 2.8 shows the various stages of the strain changes surrounding the crack; in addition, the complete crack path is also shown in relation to the strain zones. At 40 K cycles, the crack is moving in a downward direction, and upon application of the HF

2 Characterization of High Frequency Pulse Loading on Fatigue of Metals 13 Fig. 2.8 Strain field evolution, block 2C at HF of 35 Hz, and peak amplitude of 2600 N Fig. 2.9 Wheeler and Willenborg models pulse, a high strain hardening region is developed ahead of the crack, acting like a hard inclusion. The DIC data show how strains changed from 42 K cycles to 52 K cycles, including its size and dimensions; it includes utilization of work by Bruck [16], and modification work by Yates was utilized to locate the crack tip location [17] and associated plastic zone [18, 19]. The DIC data capture a kinking of the crack, also observed in microscope measurements, and how it circumvents around this are of high strain in order to continue its growth. This highlighted future work needed in investigating the strain fields ahead of the crack and shear effects to better understand the kinking behavior. This all leads to gaining a deeper understanding of the behavior observed by Laseure [7]. At the core of this work are the models by heeler [4] and Willenborg [5]. These models look at the plastic zone behavior, Wheeler relies on crack growth rates, and Willenborg utilizes stress intensity for the calculation, as shown in (Fig. 2.9): More recent research work by Mehrzadi and Taheri [6] and Chen [19] outline the behavior of the plastic zone ahead of the crack tip. The common thread of all these models is that an overload plastic zone is created once an overload/underload occurs, and the magnitude of this region has an effect on the crack growth rate as long as the cyclic plastic zone stays within its boundaries. The radii of the plastic and overload zone can be calculated utilizing the Irwin approximation [20]. The

14 P. A. Lara et al. Fig. 2.10 Wheeler and Willenborg retardation predictions of HF pulse of 50 Hz and 2900 N Wheeler model relies on the Irwin calculation to develop the radii of the overload zone and cyclic zones. The methodology of the Wheler model is that an overload plastic zone is created and retardation occurs as long as the reoccurring cyclic zone is within the bounds of the overload zone as the crack grows. The Willenborg model follows a very similar approach but relies on the stress intensity from the overload [15] and the calculation of a stress intensity reduction in order for retardation effects to seize. This behavior can be seen graphically showing how the residual stress intensity is accumulated after the HF pulse at 42 K cycles and how it diminishes over time to no retardation at 56 K cycles for the Willenborg model, Fig. 2.10. For the Wheeler model, the crack recovery can also be plotted to show when the radii of the cyclic plastic zone reaches the boundary of the overload zone and also hence ends the retardation cycle at 50 K cycles. One of the novel findings of this work was the observation of a crack kinking mechanism resulting from the application of certain HF pulses, outlined in the next section. 2.5 Crack Kinking Behavior An observed crack kinking behavior was mentioned in previous sections, during the experiments; this behavior was observed and tracked for the various specimen conditions outlined in Table 2.1. The microscope data shown in Fig. 2.11 show the crack growing from the left notch of each specimen, the orange arrow marks the direction of crack growth, the red mark shows the application of the applicable HF pulse, the green marks delimitate locations of data extraction in the fatigue cycle. While the microscope measurements provide a view of the crack propagation, a Barker’s etch method [21] was used for anodizing the aluminum alloy and then utilizing optical microscopy with polarized lights to observe the grain boundary and the propagated crack. This was evaluated at 500μm, 100μm, and20μm magnification levels, allowing the evaluation of the crack propagation through the grain, evaluating if there was a presence of large material inclusions, as well as investigating shielding effects such as contact and deflection on crack growth [22]. Fractography was also utilized to examine the surfaces through the crack depth on the specimen; several test specimens were carefully cut to allow views of the through thickness crack front. The results showed a distinct beach mark shape due to the HF pulse and a change in contrast near the edge surfaces due to the inclination associated with crack kinking. This evaluation also showed that the width of the beach mark is consistent with plastic zone and narrows beneath the surface as it is typical of plane stress to plane strain transition regions. This assessment did not indicate that there was a presence of a material inclusion that could cause the observed kinking behavior. The crack kinking behavior observed can be segmented into several areas associated with the kinking condition, as shown in Fig. 2.12. After the application of the HF pulse, there is a section when the crack continues to grow, and before it kinks, another section is when the crack is fully kinked and continues in that condition, and a last section where the crack kinking behavior has ceased and the crack has returned to its semi-original growth path. Length and angle parameters were gathered for the crack.L1 defines the crack length after the HF pulse and before the kink, L2 defines the length of the kink. The angle

2 Characterization of High Frequency Pulse Loading on Fatigue of Metals 15 Fig. 2.11 Crack kinking behavior matrix Fig. 2.12 Crack kinking characterization β1 defined the entry angle in relation to the horizontal plane, β2 is the kink angle between L1 and L2, andβ3 is the exit angle with respect to the horizontal. It is worth noting that the entry angle β1 is dependent on the original direction of crack growth and independent of the HF pulse condition, as specimens do not have a consistent entry angle β1 of crack growth, which can be dependent on notch conditions and presence of inclusions during initial crack growth. The crack kink data highlighted several findings; one was that the kink angle β2 is consistent between experiments and has a mean of 136 degrees as measured between L1 and L2. Another finding was that the exit angle β3 is consistently less than the entry angle β1 by a mean of 2 degrees.

16 P. A. Lara et al. 2.6 Effects of HF Pulse on Da/DN The data also allowed for the identification of the variation of Da/DN, as shown in Fig. 2.2, but with higher fidelity to who show the effects of an application of the HF pulse. In Fig. 2.13, the left figure shows that retardation effects are quite consistent over several test specimens that received the same HF pulse, but these effects are not quite present as shown on the right figure, where the HF pulse increased in frequency to 50 Hz, but maintained a maximum load of 2300 N. This observation can also be seen on Fig. 2.11, where no kinking is observed for the specimen on the lower left quadrant. While the results do not show a very clear crack retardation from certain HF pulses, it does provide data to show that the experimental matrix was able to capture limitations toward the presence of kinking behavior, where certain loading conditions do not precipitate a clear kinking of the crack, but instead less visible retardation effect. By leveraging the characterization of the crack kinking, as shown in Fig. 2.12, it is possible to examine the crack growth rate for varying frequency and load. Figure 2.14 shows these values and certain observations can be gathered from the results. The behavior of increased crack growth rate after the HF pulse and decrease during the kink portion can be observed. Another observation is that by varying the HF pulse form 2300–2900 N and maintaining a frequency of 20 Hz, affects the Da/DN after the HF pulse, but also reduces Da/DN during the kinking portion of the crack growth. However, by holding the load constant at 2900 N and increasing frequency, the Da/DN during the kink actually increases. This data was very promising as it shows The DIC data allows for the interrogation of surface strains ahead of the crack path Fig. 2.13 HF pulse applied at 20 Hz and 2300 N (left) and 50 Hz and 2300 N (right) Fig. 2.14 Da/DN, variation of load and frequency

2 Characterization of High Frequency Pulse Loading on Fatigue of Metals 17 Fig. 2.15 Von Mises strain fields ahead of crack along different vectors, and to see how they evolve over time. Figure 2.15 shows these Von Mises strain fields at ±45◦ from the crack direction, as well as along original crack path. That the matrix captures a variation of crack growth rates, as well as kinking behavior which is affected by both magnitude and frequency. The fanning of the strain fields changes due to the application of the HF pulses and has a very distinctive variability when the kink is complete. Work is ongoing to utilize DIC strain and crack opening displacement (COD) to develop a flow model that steps through different conditions that would trigger kinking effects and their termination. As initially expected, DIC data are providing critical data in understanding how the surface strains vary and their correlation to the effects of HF pulses. 2.7 Conclusions Existing models like the Wheeler, Willenborg, and variations of these have been utilized to predict the crack growth behavior with varying degrees of success. In this work, we carried through an experimental matrix to explore the effects of HF pulses on fatigue crack growth in 5xxx aluminum, typically used in marine applications. We utilized visual inspection of crack tip location and used digital image correlation (DIC) techniques for characterization of the crack tip deformation fields. These techniques were used on a designed center crack tension (CCT) panel specimens to measure crack growth. DIC techniques also enabled additional analysis of strain fields to better understand the associated crack growth behavior. A novel crack kinking behavior was observed from certain HF pulses; this kinking behavior was characterized and several metrics were extracted to show its influence on crack growth rates and retardation. This ongoing work includes utilizing this observed data to finalize an embedded high frequency load models for fatigue crack growth and exploring the phasing of HF pulses. 2.8 Future Work This experimental data is providing basin information toward the kinking behavior and crack retardation from the application of HF pulses, but it also highlights areas where additional work is needed. There are long-term efforts and shorter-term actions on this. Long term there is the evaluation of heat-affected zones and weld effects, introduction of other materials, and

18 P. A. Lara et al. expansion to include other variables such as HF pulse shift or decay rates. Short-term actions aim to target some particular critical efforts in gaining a deeper understanding of this behavior. This includes work in relation to the sequencing of HF pulses to gain understanding on the effects of kinking and delay. Acknowledgments Work described was performed by the Naval Surface Warfare Center Carderock Division’s Platform Integrity Department and the University of Maryland College Park’s Department of Mechanical Engineering. Financial and technical support was provided by an NSWCCD In-house Laboratory Independent Research (ILIR) program under Dr. Jack Price and a grant provided to UMD by Program Officer Dr. Paul Hess of the Office of Naval Research Code 331 under grant number N000141812016. References 1. Murthy, R., Palani, G., Iyer, N.: State-of-the-art review on fatigue crack growth analysis under variable amplitude loading. IE(I) J., 12 (2004) 2. Sumi, Y.: Fatigue crack propagation in marine structures under seaway loading. Int. J. Fatigue, 7 (2014) 3. Fricke, W., Paetzold, H.: Experimental Investigations on Fatigue Damage of Ship Structures Caused by Whipping Stresses. In: PRADS, Changwon City, Korea (2013) 4. Wheeler, O.: Spectrum loading and crack growth. J. Basic Eng. Transport. ASME, 5 (1972) 5. Willenborg, J., Engle, R., Wood, H.: A crack growth retardation model using an effective stress concept. AFFDL TM-71-1-FBR, Jan 1971 (1971) 6. Mehrzadi, M., Taheri, F.: A material sensitive modified wheeler model for predicting the retardation in fatigue response of AM60B due to an overload. Int. J. Fatigue, 10 (2013) 7. Laseure, N.S.I.M.N.D.W.W.: Effects of Variable Amplitude Loading on Fatigue Life. Ghent University, Labo Soete, Ghent (2016) 8. Hart, D., Bruck, H.: Characterization and modeling of low Modulus composite patched center crack tension specimen using DIC surface. In: Society of Experimental Mechanics (2018) 9. ASTM-E8: ASTM E8M-08: Standard Test Methods for Tension Testing of Metallic Materials. ASTM International, Conshohocken (2008) 10. ASTM-E466: ASTM E466–07: Standard Practice for Conducting Force Controlled Constant Amplitude Axial Fatigue Tests of Metallic Materials. ASTM International, Conshohocken (2007) 11. ASTM-E647: Standard Test Method for Measuring Fatigue Crack Growth Rates. ASTM International, Conshohocken (2015) 12. ASTM-E338: Standard Test Method of Sharp-Notch Tension Testing of High-Strength Sheet Materials. ASTM International, Conshohocken (2003) 13. Sutton, M., Orteu, J.-J., Schreier, H.: Image Correlation for Shape, Motion, and Deformation Measurements. Springer LLC, Boston (2009) 14. Jones, E.: Good practices guide for digital image Correllation. Int. Digital Image Correl. Soc. (2018) 15. Bannantine, J., Comer, J., Handrock, J.: Fundamentals of Metal Fatigue Analysis. Prentice Hall Inc, Upper Saddle River (1990) 16. Bruck, H.: Analysis of 3-D Effects near the Crack Tip on Rice’s 2-D J-Integral Using Digital Image Correllation and Smoothing Techniques, M.S. Thesis, Univeristy of South Carolina (1989) 17. Yates, J., Zanganeh, Y.: Quantifying crack tip displacement fields with DIC. Eng. Fracture Mech., 14 (2010) 18. Gdoutos, E.: Fracture Mechanics, An Introduction, 2nd edn. Springer, Minneapolis (2005) 19. Chen, F., Wang, F., Cui, W.: Fatigue life prediction of engineering structures subjected to variable amplitude loading using the improved crack growth rate model. Fatigue Fracture Eng. Mater. Struct., 13 (2011) 20. Anderson, T.: Fracture Mechanics Fundamental and Applications, 3rd edn. Taylor and Francis Group LLC, Boca Raton (2005) 21. Cerri, E., Evangelista, E.: Metallography of Aluminum Alloys. Mechanical Engineering Department, European Aluminum Association, Ancona (1999) 22. Ritchie, R.: Mechanisms of fatigue crack propagation in metals, ceramics and composites: role of crack tip shielding. Mater. Sci. Eng. A103, 14 (1988)

RkJQdWJsaXNoZXIy MTMzNzEzMQ==