Development of A Generalized Isotropic Yield Surface for Pressure Insensitive Metal Plasticity Considering Yield Strength Differential Effect in Tension, Compression and Shear Stress States
by Kivanc Sengoz
B.Sc. in Aeronautical Engineering, June 2001, Istanbul Technical University M.Sc. in Mechanical Engineering, May 2006, Istanbul Technical University
A Dissertation submitted to The Faculty of The School of Engineering and Applied Science of The George Washington University in partial fulfillment of the requirements for the degree of Doctor of Philosophy
January 31, 2017
Dissertation directed by
Azim Eskandarian Professor of Engineering and Applied Science
Cing-Dao (Steve) Kan Professor of Engineering and Applied Science The School of Engineering and Applied Science of The George Washington University certifies that Kivanc Sengoz has passed the Final Examination for the degree of Doctor of Philosophy as of June 29, 2015. This is the final and approved form of the dissertation.
Development of A Generalized Isotropic Yield Surface for Pressure Insensitive Metal Plasticity Considering Yield Strength Differential Effect in Tension, Compression and Shear Stress States
Kivanc Sengoz
Dissertation Research Committee:
Azim Eskandarian, Professor of Engineering and Applied Science, Dissertation Co-Director
Cing-Dao (Steve) Kan, Professor of Computational Solid Mechanics, Dissertation Co-Director
Majid T. Manzari, Professor of Civil and Environmental Engineering, Committee Member
Pedro F. Silva, Professor of Civil and Environmental Engineering, Committee Member
Tarek Omar, Chief Scientist and Technical Advisor for Federal Railroad Administration, U.S. Department of Transportation, Committee Member
ii Dedication
To my mother
iii Acknowledgments
I would like to thank Prof. Azim Eskandarian for his academic supervision, personal support, and constant inspiration throughout my studies in graduate school. I also would like to acknowledge the financial, academic, and technical support that he provided to me as the
Director of the National Crash Analysis Center (NCAC).
I would like to express my sincere gratitude to Prof. Cing-Dao (Steve) Kan, who has been my research advisor and supervisor during my doctoral research. I also would like to acknowledge the financial, academic, and technical support that he provided to me during his service at GWU as the Director of the National Crash Analysis Center (NCAC).
I would like to acknowledge Prof. Majid T. Manzari and Kim Roddis, as chairs of the
Civil & Environmental Engineering Department of the George Washington University for their invaluable support during my doctoral studies.
I also would like to take this opportunity to extend my gratitude to Dr. Gilbert (Chip)
Queitzsch for bringing the invaluable industrial and governmental insight to my research. His guidance has been very influential in helping me to understand the scientific impact and potential applications that can come out of doctoral research.
It has been a privilege having the opportunity to work with Mr. Paul A. Du Bois, who is an internationally known expert in computational mechanics, during the course of my research. His scientific approach to a problem has shown me admirable ways of being a true scholar.
iv This research was funded by the Federal Aviation Administration. I would like to acknowledge Mr. William Emmerling and Mr. Donald Altobelli of the Federal Aviation
Administration’s Aircraft Catastrophic Failure Prevention Research Program, and the
William J. Hughes Technical Center, Atlantic City, N.J., for their technical and financial support throughout this research.
I would like to acknowledge Dr. John O. Hallquist of Livermore Software Technology
Corp. (LSTC) for providing me the required tools during the implementation of the material model into LS-DYNA.
I would like to acknowledge Prof. Amos Gilat and Dr. Jeremy Seidt of Ohio State
University (OSU) for their cooperation within regarding the uncontained aircraft engine failure research program and for providing material test data.
v Abstract
Development of A Generalized Isotropic Yield Surface for Pressure Insensitive Metal Plasticity Considering Yield Strength Differential Effect in Tension, Compression and Shear Stress States
Due to effects of twinning and texture evolution, the yield surface for hexagonal close packed (HCP) metals displays an asymmetry between the yield in tension and compression, and significantly changes its shape with accumulated plastic deformation. Traditional von
Mises von Mises ( ) plasticity yield criteria and hardening assumptions can not accurately model these phenomena due to lack of stress state dependent yield stress and hardening approach. In von Mises plasticity, yield stress in uni-axial tension is always equal to yield stress in uni-axial compression, and yield stress in pure-shear is always 0.577 times the yield stress in uni-axial tension. The same assumptions are also valid for hardening behavior. In this dissertation, a generalized isotropic yield surface model that can describe both the yielding asymmetry and distortional hardening (stress state dependent hardening) on the pressure independent plastic response of metals is proposed. Yielding is described by a newly developed macroscopic phenomenological yield function that accounts for asymmetry between yielding in uni-axial tension and uni-axial compression states. This yield criteria also provides flexibility for the yield stress in pure-shear compared to von Mises ( ) plasticity; therefore it could be used by face-centered cubic (FCC) and body-centered cubic
(BCC) materials as well as HCP materials. The coefficients involved in this proposed yield criterion are considered to be functions of the accumulated plastic strain. The proposed model was implemented into the explicit finite element code LS-DYNA and used to simulate the uni-axial tension, uni-axial compression, and torsion tests of Al2024-T351 and Ti-6Al-
4V materials in the rolling direction using smooth dogbone, solid cylinder, and thin walled
vi tube specimens respectively. Tabulated inputs of characterized material test results are directly used for the constitutive treatment of the new material model. Comparisons between predicted and measured force-displacement curves and macroscopic strain fields show that the proposed model describes very well the yield strength differential effect and distortional hardening (stress state dependent hardening) based on the uni-axial tension, uni-axial compression and pure-shear stress states. On the other hand, the anisotropic nature of Ti-6Al-
4V causes directionally dependent plastic flow in some cases, and this phenomenon cannot be treated by the current material model due to the isotropic plasticity assumption.
The proposed model is then extended to include the effects of strain rate and temperature softening related with the temperate increase within the material due to adiabatic heating using tabulated input curves for yield-temperature and yield-strain rate dependency. Then a new, viscoplastic material model with generalized yield surface is implemented into a non- linear explicit dynamics finite element code, LS-DYNA. In this material model, a higher order (higher than quadratic order) generalized yield function with isotropic and isochoric plasticity is utilized, and stress state dependent isotropic hardening, strain-rate hardening and temperature softening is considered. The model takes adiabatic heating and resulting softening into account due to the plastic work. The constitutive relation is specifically developed to cover pressure insensitive plastic behavior of metals as a function of state-of- stress, strain-rate, and temperature.
The proposed strain rate and temperature dependent model was implemented into LS-
DYNA and used to simulate yield-temperature dependency of Al2023-T351 and yield-strain rate dependency of Ti-6Al-4V. Component based specimen tests are used to validate and show the robustness, accuracy and efficiency of the new material model. The post
vii experiment data and simulation results are shown to be in very good agreement with each other.
It is shown that the new material model is capable of predicting yielding asymmetry between uni-axial tension and uni-axial compression and pure-shear stress states accurately within the limits of the convexity region. Stress state, strain rate, and temperature dependent plastic flow prediction capability of the model suggests that the new material model can be used as a promising tool for diverse applications of dynamic plastic deformation and subsequent ductile failure phenomenon. Therefore, the material model can be used in the numerical solution of the impact problems in industry such as fan blade-out containment events in turbofan engines, vehicle crashworthiness, and metal forming.
viii Table of Contents
Pages
Dedication ...... iii
Acknowledgments ...... iv
Abstract ...... vi
Table of Contents ...... ix
List of Figures ...... xii
List of Tables ...... xxi
List of Symbols ...... xxii
List of Acronyms ...... xxv
CHAPTER 1 - Introduction ...... 1
1.1 Introduction ...... 1
1.2 Background and Motivation ...... 3
1.3 Research Objectives and scope ...... 5
1.4 Original Contributions ...... 6
1.5 Outline of the Dissertation ...... 7
CHAPTER 2 - Literature Review ...... 9
2.1 Introduction to Yielding ...... 9
2.1.1 Isotropic yield criteria for metals ...... 9
2.1.2 Deviations from plasticity and higher order yield surfaces ...... 14
ix 2.1.3 Role of Third Deviatoric Stress Invariant ( ) in Modeling Asymmetry between Tensile and Compressive Yield ...... 19
CHAPTER 3 - Methodology and Theoretical Framework ...... 25
3.1 Fundamentals of continuum plasticity ...... 25
3.1.1 Perspectives on Plastic Constitutive Equations: Principle of Maximum
Plastic Resistance and Drucker’s Postulate ...... 30
3.1.2 Consequences of Drucker’s hypothesis ...... 34
3.2 Development of proposed generalized yield surface (GYS) plasticity model involving stress invariants and ...... 35 3.3 Associated flow rule, hardening rule, and definition of effective plastic strain for GYS plasticity model...... 43
3.4 Perspectives on dependent yield surfaces: GYS, Hosford, Drucker, Cazacu-Barlat ...... 52
3.5 Convexity of the GYS model and comparison between the GYS, Drucker and Cazacu-Barlat convexity regions ...... 67
3.6 Development of convexity algorithm for the GYS model ...... 91
CHAPTER 4 - Numerical Implementation of GYS Yield Criterion ...... 98
4.1 Constitutive relation of the proposed GYS model ...... 98
4.2 Integration algorithms for elasto-plastic constitutive equations ...... 104
4.2.1 Time integration and numerical implementation of GYS into LS-DYNA ...... 108
x 4.3 Verification of numerical results using single element test cases at different stress states...... 122
4.3.1 Single element analysis and verification of GYS material model ...... 122
4.3.1.1 Benchmark-problem definition ...... 122
CHAPTER 5 - Application and validations of GYS plasticity model ...... 136
5.1 Materials, specimen properties and finite element procedures for the numerical simulation of uni-axial tension, compression and torsion (pure- shear) test cases ...... 136
5.2 Experimental results versus numerical simulations for uni-axial tension, compression and torsion (pure-shear) tests...... 145
5.3 Influence of the GYS plasticity model on the failure response compared to the von Mises plasticity model ...... 159
CHAPTER 6 - Concluding remarks and Future Work ...... 166
References ...... 169
Appendix - A: Second condition of convexity for the GYS model ...... 173
Appendix - B: Derivative of GYS yield function with respect to the increment of the plastic strain multiplier ...... 176
Appendix - C: Single element verification cases of convexity projection ...... 177
Appendix - D: Comparison of MAT224 and MAT224_GYS models under the von Mises plasticity assumption considering proportional (uni-axial tension and plain strain tension) and non-proportional loading (plain strain tension- uni-axial tension) case ...... 185
xi List of Figures
Pages Figure 1-1: Applications of dynamic failure for different structures; (a) offset frontal crash of a passenger car (http://www.iihs.org), (b) uncontained aircraft engine failure (c) fan disk failure and (d) different failure modes after ballistic impact of projectiles with different nose shapes, adapted from [1]...... 2 Figure 1-2: Failure modes in impacted plates, after [8]...... 5 Figure 2-1: Tresca and von Mises yield criteria in the π-plane...... 12 Figure 2-2: Difference of Tresca and von Mises yield criteria for pure-shear state in the π-plane...... 13 Figure 2-3: Comparison of yield surfaces for different materials under combined loading; (a) after [17] and (b) after [21]...... 13 Figure 2-4: Plane stress yield loci for Drucker, Tresca and von Mises (normalized by uniaxial tension) after Drucker [24]...... 15 Figure 2-5: Plane stress yield loci for Hosford, Tresca, von Mises (a) n=1 Tresca, n=2 von Mises and n=10 Hosford (b) plane stress yield loci predictions for different n values between 1 and 2.767 as shown by Hosford in 1972 [23]...... 16 Figure 2-6:. Different isotropic plane stress yield surfaces between the von Mises yield surface and the lower bound (Tresca) after Karafillis and Boyce [26]...... 18 Figure 2-7: Different isotropic plane stress yield surfaces between the von Mises yield surface and the upper bound after Karafillis and Boyce [26]...... 18 Figure 2-8: The four reference points to construct the Vegter yield function in plane stress after [31]...... 19 Figure 2-9: Elastic deformation corresponding to pure stretching of the crystal lattice (top) and plastic deformation trough dislocation slip (bottom) [36] [37]...... 21 Figure 2-10: Crystal lattice reorientation due to mechanical twinning [36][37] ...... 21
Figure 2-11: Plane stress yield loci corresponding to / =2/3, 3/3 ( von Mises), 4/3 [40] ...... ...... 22 Figure 2-12: A comparison of force displacement curves between experimental results and simulation results for flat grooved specimen (Bai & Wierzbicki [41]) ...... 23
xii Figure 3-1: (a) Isotropic hardening and (b) kinematic hardening in plane stress yield loci [47] ...... 28 Figure 3-2: Hardening response of pure titanium according to theoretical yield surface experimental yield data corresponding to fixed values of equivalent plastic strain after [48] ...... 29 Figure 3-3: Stable and unstable material behavior on stress-strain curve according to Drucker’s stable material postulate [45] ...... 32 Figure 3-4: Properties of a yield surface with associated flow rule (a) normality (b) corner condition (c) convexity [47][50][51] ...... 33 Figure 3-5: Representation of Lode parameter and UT (uniaxial tension), UC (uniaxial compression) and PS (pure-shear) stress states (b) Comparison of yield surface models in deviatoric plane ( ᴨ–plane) (c) effect of tension ( σT)-compression ( σC) yielding asymmetry on the yield surface in deviatoric plane…………………………………………………..37 Figure 3-6: Uni-axial tension, uni-axial compression and pure-shear (torsion) tests required by GYS to determine model coefficients , , and ...... …41 , Figure 3-7: (a) Four internal curves computed by GYS for initial yield and subsequent hardening (b) four equations for the determination of model coefficients , , and ...... 42 , Figure 3-8: Cazacu-Barlat yield surface in the plane stress space corresponding to uni-axial tension/uni-axial compression yield ratios, 3/2, 3/3, 3/4...... 57 Figure 3-9: Comparison between the Cazacu-Barlat and GYS yield surface models in the plane stress space corresponding to uni-axial tension/uni-axial compression yield ratios, 3/2, 3/3, 3/4...... 57 Figure 3-10: related scale factor on von Mises stress for Hosford yield function ...... 60 Figure 3-11: Comparison of related scale factor on the von Mises stress for the GYS and Hosford yield functions ...... 61 Figure 3-12: Comparison of scale factor on von Mises stress for GYS with and without tension- compression asymmetry corresponding to order (m) 6 ...... 62 Figure 3-13: Hosford non-unique flow ...... 62 Figure 3-14: GYS non-convex flow ...... 63
xiii Figure 3-15: Example of the plane strain condition for the Hosford yield function and non-unique flow characteristic shown by the yield function as exponent m becomes higher ...... 64 Figure 3-16: a = value for plane strain condition of the Hosford yield function ...... / ...... 65 Figure 3-17: Example of plane strain condition for the GYS function and unique flow characteristic shown by the yield function as order m gets higher...... 66 Figure 3-18: a = value for plane strain condition of the GYS model ...... 66 Figure 3-19: Convexity / region of the GYS model in terms of yield stress ratio of uni-axial compression to uni-axial tension and pure-shear to uni-axial tension and Lode parameter ...... 77 Figure 3-20: Convexity region of the GYS model in terms of yield stress ratio of uni-axial compression to uni-axial tension and pure-shear to uni-axial tension and Lode parameter ...... 78 Figure 3-21: (a) Convexity region of GYS model with yield ratio uni-axial compression/uni-axial tension=1 (b) Lower and upper boundaries of the convexity region correspond to pure-shear yield to uni-axial tension yield ratio 0.545 and 0.611 with the yield ratio uni-axial compression/uni-axial tension=1 ...... 79 Figure 3-22: Convexity region of GYS model with yield ratio uniaxial compression/tension=0.95...... 81 Figure 3-23: Convexity region of GYS model with yield ratio uniaxial compression/tension=0.75 ...... 81 Figure 3-24: Convexity region of GYS model with yield ratio uniaxial compression/tension=1.1 ...... 83 Figure 3-25: Convexity region of GYS model with yield ratio uniaxial compression/tension=1.25 ...... 83 Figure 3-26: (a) Convexity region of GYS model with yield ratio compression/tension=2.00 (b) Convexity region of GYS model with yield ratio uniaxial compression/tension=0.50...... 85 Figure 3-27: Convexity curves of the GYS model corresponding to the Lode angle parameter (a) 0, (b) 1, (c) -1, respectively ...... 86
xiv Figure 3-28: (a) Intersection of the Lode angle parameter ( ) curves corresponding to 0, 1, -1 for the GYS model (b) Comparison of the convexity regions of the GYS model and von Mises, Drucker, and Cazacu-Barlat models...... 87 Figure 3-29: Min and Max points on the convexity region of GYS model in terms of yield stress ratios ...... 88 Figure 3-30: Min and Max tip points on the convexity region of GYS model ...... 92 Figure 3-31: Projection operation of convexity algorithm of GYS model ...... 93 Figure 4-1: Schematics of the return mapping algorithm in stress space ...... 105 Figure 4-2: Schematics of the return mapping algorithm in stress space (a) Fully implicit (b) semi-implicit tangent cutting plane ...... 106 Figure 4-3: The role of elasto-plastic return mapping algorithms ...... 110 Figure 4-4: Schematic illustration of the return mapping algorithm for the GYS ...... 110 Figure 4-5: Stress update scheme using trial stress and a plastic correction ...... 113 Figure 4-6: (a) Geometric interpretation of the radial-return in principle stress space (b) GYS return mapping where mapping direction is the yield surface normal at the trial stress point B ...... 115 Figure 4-7: Geometric interpretation of the return on the invariant plane ...... 116 Figure 4-8: Schematics of the solution interval and successive secant iterations to find the increment in plastic strain assuming that root of the function ( f=0) is found at 4 th iteration step [71] ...... 117 Figure 4-9: Successive secant iterations and dependence of function f on the increment of equivalent plastic strain. Note that curvature of the function f could be both ways (a) and (b) without effecting monotonic decrease of the function [70][71]...... 118 Figure 4-10: (a) Rate dependent tensile hardening curve from a dynamic test ...... 119 Figure 4-11: Coupling between GYS and MAT224 ...... 120 Figure 4-12: Single element model at different state of stresses ...... 122 Figure 4-13: Axial stress (Z stress) versus Effective stress (von Mises) ...... 123 Figure 4-14: GYS-Pure-shear case. Shear stress (YZ stress) versus Effective stress (von Mises) ...... 123
xv Figure 4-15: GYS-Uniaxial compression case. Axial stress (Z stress) versus Effective stress (von Mises) ...... 124 Figure 4-16: GYS-Plain strain tension case. Axial stress (Z stress) versus Effective stress (von Mises) ...... 124 Figure 4-17: GYS-MAT24 comparison under the von Mises plasticity assumption .... 125 Figure 4-18: Temperature dependent yield curves for the temperature softening case .. 125 Figure 4-19: Temperature softening—uniaxial tension single element case ...... 126 Figure 4-20: Temperature softening—uniaxial compression single element case ...... 127 Figure 4-21: Temperature softening—temperature increase on uniaxial compression single element case ...... 127 Figure 4-22: Temperature softening—temperature increase on uniaxial tension single element case ...... 128 Figure 4-23: Strain rate dependent yield curves for single element strain rate dependency case ...... 128 Figure 4-24: Strain rate single element test in compression [strain rate=1/s] ...... 129 Figure 4-25: Strain rate single element test in compression [strain rate=100/s] ...... 129 Figure 4-26: Stress increment and return mapping using a plastic direction estimated from the beginning of time step for (a) small strain increment and a strain increment more than an order of maginute larger than yield strain [60] ...... 130 Figure 4-27: Effect of time step on the GYS response-uniaxial tension single element case ...... 132 Figure 4-28: Effect of time step on the GYS response-uniaxial coompression single element case ...... 132 Figure 4-29: Effect of time step on the GYS response-simple shear single element case ...... 133 Figure 4-30: Effect of time step on the GYS response-nonproportional loading case-uniaxial tension-biaxial tension (axial stress) ...... 133 Figure 4-31: Effect of time step on the GYS response-nonproportional loading case-uniaxial tension-biaxial tension (lateral stress) ...... 134 Figure 4-32: Effect of time step on the GYS response-nonproportional loading case plain strain tension-uniaxial tension (axial stress) ...... 134
xvi Figure 4-33: Effect of time step on the GYS response-nonproportional loading case plain strain tension-uniaxial tension (lateral stress) ...... 135 Figure 5-1: Al2024-T351 and Ti-6Al-4V specimen manufacturing orientations for the 12.7 mm stock ...... 138 Figure 5-2: Reported composition of a 12.7 mm thick Al2024-T351 plate ...... 138 Figure 5-3: Reported composition of a 12.7 mm thick Ti-6Al-4V plate ...... 138 Figure 5-4: (a) Instron 1321 servo-hydraulic biaxial load frame located in OSU DMML (b) Elevated temperature test set up ...... 139 Figure 5-5: Tension and Torsion Split-Hopkinson bar located in OSU DMML ...... 139 Figure 5-6: Compression Split-Hopkinson bar located in OSU DMML ...... 139 Figure 5-7: Dimensions of a thin flat dog-bone spceimen for uni-axial tension test ...... 140 Figure 5-8: Dimensions of a thin walled tube specimen for torsion test ...... 140 Figure 5-9: (a) Tension and (b) Compression temperature test series data for Al2024-T351 from OSU DMML [16] ...... 142 Figure 5-10: (a) Tension and (b) Compression strain rate test series data for Al2024-T351 from OSU DMML [16] ...... 143 Figure 5-11: (a) Tension and (b) Compression temperature test series data for Ti- 6Al-4V from OSU DMML [76] ...... 144 Figure 5-12: (a) Tension and (b) Compression strainrate test series data for Ti- 6Al-4V from OSU DMML [76] ...... 145 Figure 5-13: Finite element simulation model and set up: (a) Tension simulation and (b) Compression simulation ...... 147 Figure 5-14: Finite element simulation model and set up: Torsion simulation ...... 147 Figure 5-15: Al2024-T351 uni-axial tension test and simulation at room temperature ...... 149 Figure 5-16: Maximum principal strain contours of Al2024-T351 tension test specimen at room temperature (a) simulation using GYS (b) Test (DIC) ...... 149 Figure 5-17: Al2024-T351 uni-axial compression test and simulation at room temperature ...... 150 Figure 5-18: Al2024-T351 uni-axial tension test and simulation at 150 oC ...... 150 Figure 5-19: Al2024-T351 uni-axial compression test and simulation at 150 oC ...... 151
xvii Figure 5-20: Al2024-T351 uni-axial tension test and simulation at 300 oC ...... 151 Figure 5-21: Al2024-T351 uni-axial compression test and simulation at 300 oC ...... 152 Figure 5-22: Al2024-T351 torsion test and simulation at room temperature ...... 152 Figure 5-23: Al2024-T351 torsion test at room temperature (a) simulation using GYS (b) Test (DIC) ...... 153 Figure 5-24: Ti-6Al-4V uni-axial tension test and simulation at strain rate of 1/s at room temperature ...... 153 Figure 5-25: Ti-6Al-4V uni-axial tension test and simulation at strain rate of 1/s at room temperature (a) simulation using GYS (b) Test (DIC data) ...... 154 Figure 5-26: Ti-6Al-4V uni-axial tension test and simulation at strain rate of 1/s at room temperature (a) simulation using GYS (b) Test (DIC data) ...... 154 Figure 5-27: Ti-6Al-4V uni-axial compression test and simulation at strain rate of 1/s at room temperature ...... 154 Figure 5-28: Axial strain distribution of Ti-6Al-4V uni-axial compression test at strain rate of 1/s at room temperature (a) simulation using GYS (b) Test (DIC data) .. 155 Figure 5-29: Ti-6Al-4V torsion test and simulation at strain rate of 1/s at room temperature ...... 156 Figure 5-30: Shear strain distribution on Ti-6Al-4V torsion test specimen at strain rate of 1/s at room temperature (a) simulation using GYS (b) Test (DIC data) ...... 156 Figure 5-31: Maximum principal strain distribution on Ti-6Al-4V torsion test specimen at strain rate of 1/s at room temperature (a) simulation using GYS (b) Test (DIC data) ...... 156 Figure 5-32: Ti-6Al-4V uni-axial tension test and simulation at strain rate of 1500/s at room temperature ...... 157 Figure 5-33: Ti-6Al-4V uni-axial compression test and simulation at strain rate of 1500/s at room temperature ...... 157 Figure 5-34: Axial strain distribution on Ti-6Al-4V uni-axial tension test specimen at strain rate of 1500/s at room temperature ...... 158 Figure 5-35: Axial strain distribution on Ti-6Al-4V uni-axial tension test specimen at strain rate of 1500/s at (a) simulation using GYS (b) Test (DIC data) ...... 158
xviii Figure 5-36: Axial stress for yield in uni-axial Compression/Tension=0.9 (uni-axial tension case) ...... 159 Figure 5-37: Axial stress for yield in uni-axial Compression/Tension=0.9 (uni-axial compression case) ...... 160 Figure 5-38: Equivalent plastic strain for yield in uni-axial Compression/Tension=0.9 (uni-axial compression case) ...... 160 Figure 5-39: Axial stress for yield in uni-axial Compression/Tension=1.1 (uni-axial compression case) ...... 161 Figure 5-40: Equivalent plastic strain for yield uni-axial Compression/Tension=1.1 (uni-axial compression case) ...... 161 Figure 5-41: Axial stress for a yield stress ratio of uni-axial compression to tension varies from 0.9 to 1.1 (uni-axial compression case) ...... 162 Figure 5-42: Effect of the Uni-axial Compression/Tension yield ratio on failure with the GYS uni-axial compression case ...... 162 Figure 5-43: Effect of the Pure-Shear/Uni-axial Tension yield stress ratio on failure using GYS (pure-shear case) ...... 163 Figure 5-44: Effect of the Pure-Shear/Uni-axial Tension yield stress ratio on failure using GYS (pure-shear case) ...... 163 Figure 5-45: Effect of the Pure-Shear/Uni-axial Tension yield stress ratio on failure using GYS (Non-proportional loading case pure shear-uniaxial tension) ...... 164 Figure 5-46: Effect of the Pure-Shear/Uni-axial Tension yield stress ratio on failure using GYS (Non-proportional loading case pure-shear-uniaxial tension) ...... 164 Figure 5-47: Effect of loading sequence on failure using GYS-maximum principal stress ...... 165 Figure 5-48: Effect of loading sequence on failure using GYS-equivalent plastic strain ...... 165 Figure C-1: The GYS convexity region and data point projection...... 177 Figure C-2: Verification cases of projection algorithm for the data set out of convexity region...... 177 Figure C-3: Verification case-2 ( =1.8) and projection point on the convexity region...... ...... / ...... 178
xix Figure C-4: Verification case-2 ( =1.8) (a) tension and compression effective stresses without scale factor (b) scale/ factor of 1.693 on the tension curve...... 179 Figure C-5: Verification case-3 ( =1.8) and projection point on the convexity region...... ...... / ...... 180 Figure C-6: Verification case-3 ( =0.5). (a) Uni-axial tension and uni-axial compression effective stresses without / scale factor. (b) Scale factor of 0.59 on the uni-axial tension curve ...... 181
Figure C-7: Verification case-4 =1.3 and projection point on the convexity region...... √3 / ...... 181
Figure C-8: Verification case-4 =1.3 (a) Uni-axial tension and pure-shear effective stresses without scale factor.√3 / (b) Scale factor of 1.0588 on the uni-axial tension curve...... 182 Figure C-9: Verification case-8 =0.3-2.29 for data point projection for non- convex cases. (a) Uni-axial tension / and uni-axial compression effective stresses without scale factor. (b) Scale factor of 1.0588 on the uni-axial tension curve. (c) Uni-axial tension and pure-shear effective stresses without scale factor. (d) Scale factor of 1.059 on the uni-axial tension curve...... 184 Figure D-1: Loading cases and evolution of the stress-state values...... 185 Figure D-2: Uni-axial tension case – Plastic strain to failure, damage, triaxiality, effective stress...... 186 Figure D-3: Plain strain tension case – Plastic strain to failure, damage, triaxiality, effective stress...... 186 Figure D-4: Non-proportional loading case (Uni-axial tension-plain strain tension) case – Plastic strain to failure, damage, triaxiality, effective stress...... 187
xx List of Tables
Pages Table 3-1: Yield stress prediction of GYS at different state of stresses...... 54 Table 3-2: Comparison of properties of yield surfaces...... 55 Table 3-3: Comparison of convexity ranges of yield surfaces ...... 75 Table 4-1: MAT224_GYS input card in LS-DYNA...... 121
xxi List of Symbols
Roman Symbols
Material constants for the plasticity model of Gao [44] , , Fourth order positive definite isotropic stiffness tensor