Anisotropic plasticity, anisotropic failure and their application to the simulation of aluminium extrusions under crash loads
|
|
- Annice Rice
- 6 years ago
- Views:
Transcription
1 Anisotropic plasticity, anisotropic failure and their application to the simulation of aluminium extrusions under crash loads Paul Du Bois, Consulting engineer Oasys LS-DYNA Users Meeting 28th January 2016 The Arup Campus, UK
2 Anisotropy from a metallurgical point of view : texture index Sample BLM 45 BLM 45 BLM 46 BLM Texture Index Predominant Texture Type Transverse Basal Transverse Transverse Transverse Basal Basal Basal, Transverse Degrees off-center Texture index denotes degree of anisotropy Index ranges from 1 (isotropic) to (single crystal) Aluminum foil (highly textured) has an index of 4.14 Anisotropy depends upon the manufacturing process Distinguish between materials with strong texture vs. materials with weak texture
3 Anisotropy from a mathematical point of view : symmetry Fully anisotropic : no symmetry Monoclinic : 1 symmetry plane f xx, yy, zz, xy, yz, zx f xx, yy, zz, xy, yz, zx Orthotropic : 2 symmetry planes f xx, yy, zz, xy, yz, zx f xx, yy, zz, xy, yz, zx f xx, yy, zz, xy, yz, zx f xx, yy, zz, xy, yz, zx Transversely isotropic : axisymmetric isotropic Most numerical models assume orthotropic response Properties in 45 degree and 135 degrees to rolling direction are assumed to be the same
4 Anisotropy from a CAE engineering point of view : process Levels of increasing complexity of simulations involving shell elements : isotropic Transversely isotropic (R<1 or R>1 ) Orthotropic with known material direction ( stamping, extrusions ) Orthotropic with material directions to be mapped ( stamped or cold formed parts ) May be a blessing in disguise as results may become less random
5 Anisotropy from an industrial point of view : formability Science of metalforming started in the UK and the US in the 1950 s Hosford, Swift, Hill, Lankford and others identified the Rparameter(s) as a measure of formability : high R is good formability as it means little thinning occurs during deformation The R-value for any angle is determined as the ratio of transversal to thickness plastic strain rate in uniaxial tension : R 22p p 33 sin 2 f f f sin 2 cos 2 xx xy yy f f xx yy R is equal to unity for an isotropic metal So formability is directly related to the transverse anisotropy and/or inplane anisotropy of the metal sheet
6 When is a material isotropic? Example for plane stress isotropic in the elastic range y00 y45 y90 Isotropic yield Isotropic flow E 00 E R 00 R 45 R 90 1 Isotropic failure pf pf pf
7 Traditional crash analysis based on MAT_024 : isotropic multi-linear plasticity Why? Because for mild steel and HSLA we can get away with it Because it is difficult to determine the rolling direction of the sheet in every IP after stamping The approach still seems to work for materials such as DP, CP, TRIP and UHSS
8 Traditional metalforming analysis Rolled sheet is orthotropic Before forming the rolling direction is known R-values are a measure of formability Thus considering orthotropy in the simulation is both possible and important First orthotropic law was Hill-48 Associated flow : anisotropic yield and anisotropic flow are related Mainly interested in small deformations ( below necking strain ) Consequently off-axis yield curves are considered scaled versions of the yield curve in the rolling direction : call this isotropic or homologuous hardening
9 Hill 48 in plane stress Quadratic orthotropic 4 parameter law (00=rolling direction ) G 002 F 902 H N Traditional way of working : input the yield curve in rolling directions and 3 R-values, compute the yield values in the off-axis directions yyp H R00 p zz G xxp H R90 p zz F p 135 N 1 R45 p zz F G F H [( F G 2 N )]
10 Typical R-values for steel and aluminium sheet
11 Higher order anisotropic material laws At Alcoa Frederic Barlat refined the anisotropic material laws with the purpose of simulating aluminium sheet Associated flow Homologuous hardening ( single yield curve ) Higher order yield surface ( order 6 or 8 ) Barlat-89 is implemented in MAT_036, 4 parameter law, reduces to Hill 48 if m=2, popular workhorse in the stamping community Barlat-2000 is implemented as MAT_133, 8 parameter law Aretz model was implemented as MAT_135 by SimLab Trondheim, 8 parameter model
12 [kn] Aluminium sheet with R= Force 4 2 Reference : Increasing rupture predictability for aluminum, influence of anisotropy Daniel Riemensperger, LS-DYNA forum 2015, Bamberg displacement 30 [mm]
13 Aluminium extrusions Lightweight material Fracture is an issue Increasing use in automotive industry Yield seems mildly anisotropic or even isotropic Flow is strongly anisotropic, indicating non-associated flow R-values can be very accurately measured using DIC Values are highly directionally dependent in extrusions
14 Aluminium extrusions For crash simulations the prime concern is and remains capturing the correct stress levels However the R-values ( anisotropic flow ) will influence the failure This leads to a consideration of anisotropy in the crash simulation Possible as the extrusion direction is known Interest goes up to failure far beyond necking, this has lead to the introduction of distortional hardening ( multiple yield curves )
15 Aluminium extrusion XXXX and Hill 48 / Barlat 89 The example below illustrates the failure of 4 parameter laws to accurately model an aluminium extrusion H G 1 H G 0.71 R G H 0.29 H R F 0.71 F N 1 N R F G 2 s N F G 2N F H R00 R45 R s test Hill
16 Aluminium extrusion XXX and Hill 48 / Barlat 89 The example below illustrates the failure of 4 parameter laws to accurately model an aluminium extrusion Comparison of MAT_024 and MAT_036 with R00=R90=0.4 and R45=1.2 R00 R45 R90 The surfaces under 45 Degree differ from 00 degree
17 The Fleisher-Borrvall model : a first shot at distortional hardening It seems desirable to match yield curves and variable R-values as the shape of yield curves in different directions may vary : y00 y p y00 E y90 p45 45 y45 E p The input of 3 yield curves is expected in the option HR=7 y90 E
18 Distortional hardening The yield curves are internally converted in function of equivalent plastic strain : y00 y45 y90 y00 y45 p45 00 p y00 y90 90 p The conversion is based on work hardening The 45 and 90 degree curves are then scaled so that all 3 curves have the same initial yield point
19 Distortional hardening : The yield function is rewritten ( unconventionally ) as : f eff 00 y00 p 45 y45 p 90 y90 p With the parameters a b c a 1 4 b 2 : q 4q12 1 q d bq12 1 b 1 q c d 1 q 45 1 c q c 4 c c 1 d 1 q c 4
20 Distortional hardening The Fleisher-Borrvall model was implemented in MAT_036 as the HR=7 option, 6 load curve model For loading in 00, 45 and 90 degrees we get resp. : f eff y00 p f eff y45 p f eff y90 p And for other loading angles we get a convex combination of the hardening curves in each direction A non-associated flow rule allows to respect R-values as specified by the user in 3 directions
21 Aluminium extrusion YYYY Material has isotropic yield
22 Alumimium extrusion YYYY Reference shows R values as: R00 = 0.48, R45 = 0.29, R90 = 1.76 Bumper Beam Longitudinal System Subjected to Offset Impact Loading Kokkula (PhD Thesis) AA-6060 T1 Aluminum R values for AW-6060 T66 are: R00 = 0.49, R45 = 0.27, R90 = 1.69
23 Aluminium extrusion YYYY
24 Aluminium extrusion YYYY Reference : Testing and Validation of a New Damage Model for Application in Automotive Crashworthiness Simulations Sean Haight, 2012 ( internal report )
25 Further enhancements of MAT_036 HR=7 Reference : Development of an Anisotropic Material Model for the Simulation of Extruded Aluminum under Transient Dynamic Loads, Anthony Smith, LS-DYNA conference Wuerzburg, 2015 Orthotropic plasticity with distortional hardening seemed the way to go for aluminium extrusions Many questions remain at this point How good is the response under shear/biaxial loads? Does the material law have sufficient free parameters?
26 6000 series extruded aluminum Some anisotropy in flow stress High anisotropy in R-Values Need to investigate anisotropic models in LS-Dyna YS0 = 180 MPa R0 = 0.56 YS45 = MPa YS90 = MPa R45 = R90 = 1.70 Tensile Test Schematic courtesy of: Structural Impact Laboratory Norwegian University of Science and Technology Engineering Strain R-Value Engineering Stress 6000 series extruded aluminium Effective Plastic Strain 26
27 MAT_036, HR=7 for pure shear 00p 1p R p p degree 00 degree 1 2 Precise definition of material system needed 00p 90p 1p 2p 2p 1 R
28 MAT_036, HR=7 for equibiaxial tension 00p 1p R p p degree 00 degree 1 2 Isotropic : independent of material system 00p 90p 1p 2p 2p 1 R
29 Enhanced formulation of MAT_036 with HR=7 The yield function is rewritten as : f eff 00 y00 p 45 y45 p 90 y90 p s ys p b yb p With the parameters b max(0,1 4 d2 ) s max(0,4 d2 1) d : c 00 1 c d 1 q 1 b s 4 c 45 1 c q 1 b s 2 c 90 1 c 1 d 1 q 1 b s 4
30 MAT_036 with HR=7 is now a 10 load curve model Enhancements to Mat RSHR RBIA 36 HR=7 X RBIAX.GT.0: Constant biaxial R-value RBIAX.LT.0: Biaxial R-value as function of (plastic) strain RBIAX = deps11/deps22 RSHR.GT.0: Constant shear R-value RSHR.LT.0: Shear R-value defined as function of (plastic) strain RSHR = - deps11/deps22 30
31 New Mat 36 HR=7 Yield Surface 90 Biaxial 45 0 Biaxial Biaxial Shear Shear 0, 45, 90 Material direction; yield surface at plastic strain = 0.0 YS0 = YS45 = YS90 = YSShear = YSBiaxial= Mat 36 HR 3 Mat 36 HR 7 Old Mat 36 HR 7 - New 31
32 Enhanced MAT_036 HR=7 model for extruded aluminium Input Plastic Strain Hardening Curves True Stress vs. Plastic Strain 0 Yield Surface 0 Material direction, yield surface evolution from plastic strain = 0.00 to 0.50 This material behavior is a good example of distortional hardening; note how the yield surface expands and distorts with increasing plastic strain. Similar behavior is observed in the 45 and 90 material directions. 32
33 Uniaxial Tension 0 Force vs. Displacement Uniaxial Tension 45 Force vs. Displacement Uniaxial Tension 90 Force vs. Displacement 33
34 Shear butterfly test Cases Force vs. Displacement 34
35 Punch test to Test Cases Biaxial Punch Test Force vs. Displacement 35
36 Plane stress orthotropic failure : directional dependency upon the state of stress
37 MAT_ADD_GENERALIZED_DAMAGE or MAGD xx D11 yy D21 D xy 41 D12 D22 D42 D14 xxeff eff D24 yy D44 xyeff xx 1 d 0 0 xxeff eff 1 d 0 yy yy d xyeff xy 37
38 introduction Isotropic failure can be modeled in LS-DYNA using MAT_ADD_EROSION_GISSMO, this model takes into account the dependency of failure on the state of stress and the load path but has no directionality In the following we will generalize GISSMO to model orthotropic failure meaning the material has 3 symmetry planes and in particular failure strains under 45 degree and under 135 degree to the material x-axis are the same, the model will at first be limited to plane stress ( shell elements ) This approach will be realized with the versatile MAT_ADD_GENERALIZED_DAMAGE option in LS-DYNA We remain consistent with orthotropic plasticity models in plane stress where material properties are specified under 0, 45 and 90 degrees to the material x-axis. Consequently the model will require GISSMO type input based on experiments with the first principal stress direction under 0, 45 and 90 degrees to the material x-axis. Thus NHIS=3.
39 introduction The model will be orthotropic with respect to failure, different failure strains can be defined in different directions for the same state of stress The dependency of failure on state of stress and load path remains the same as in GISSMO Damage coupling is isotropic : a scalar function defines the relationship between true stress and effective stress, thus IFLG3=1 How fast damage accumulates depends upon the loading direction Since all data are defined in the material system, damage must be accumulated in the material system, thus IFLG2=1 Instead of using history variables from the elasto-plastic material law, the drivers for the 3 GISSMO models are internally computed according to the theory outlayed in these notes, thus IFLG1=2
40 introduction Damage coupling is isotropic : a scalar function defines the relationship between true stress and effective stress, thus IFLG3=1 Since all data are defined in the material system, damage must be accumulated in the material system, thus IFLG2=1 Instead of using history variables from the elastoplastic material law, the drivers for the 3 GISSMO models are internally computed, thus IFLG1=2
41 Experimental database Tensile, notched and shear tests needed in 3 directions direction of the shear tests is the first principal stress direction (tensile diagonal) Only 1 biaxial test needed
42 Input data preparation LCSDG-90 LCSDG-00 LCSDG-45
43 Input data generation A consistent input deck also requires the instability curves (ECRIT) to have coincident biaxial points Swift curves generated for material hardening curves in 3 directions will usually not have coincident biaxial points
44 the plastic strain rate tensor in plane stress We can compute an equivalent plastic strain rate for every material direction : xxp p xy 1 b p xyp 1 0 b b p p p yy 0 b b b 1 1 b 1 b b 1 1 b 13 1 b ep p 1 1 b 2 b 2 1p 2 cos sin b2 b 3 90ep 2 90p 1 1 b 2 b 2 1p sin 2 cos sin 3 ep p 1 1 b 2 b 2 1p cos 2 cos sin b2 b 3 2 b
45 The plastic strain rate tensor in plane stress Plot the coefficients in function of the angle theta for unit principal strain rate : Angle theta in degrees 00p 1p cos 2 cos sin 90p 1p sin 2 cos sin p 45p 135 1p 2 cos sin 0
46 Damage accumulation Damage can now be accumulated in the material system using the failure criteria that were determined from testing in each direction : Then the damage is computed as : 1 1 n ep 00 00f 1 1 n 90ep 90f d 90 d 90 dt 1 1 n ep 45 45f d 45 d 45 dt d 00 nd 00 d 90 nd 90 d 45 nd 45 d 00 d 00 dt No problems with the sign of the damage as equivalent strain rates are all positive
47 isotropic damage Propose a damage formulation for IFLG3=1 using only the strain rate components with positive coefficients Then, assuming a positive first principal strain rate : cos sin cos 0 d d 452 d d 452 d 002 p cos 2 sin 2 cos 4 2 cos 3 sin sin sin cos 0 d d 452 d d 452 d 902 p cos 2 sin 2 sin 4 2 cos sin \
48 isotropic damage This approach results in a monotonic model ( no new maxima or minima are generated between the failure surfaces defined for degrees )
49 isotropic damage This approach results in a monotonic model ( no new maxima or minima are generated between the failure surfaces defined for degrees )
50 isotropic damage This approach results in a monotonic model ( no new maxima or minima are generated between the failure surfaces defined for degrees )
51 Large scale validation
52 Uniaxial tensile test : MAT_036 + MAT_ADD_D_G 00 degree 45 degree 90 degree
53 Conclusions Due to pronounced anisotropic flow the prediction of failure in aluminium extrusions has been an elusive goal Progress was made thanks to extensive code development in both plasticity ( MAT_036 ) and anisotropic failure (M_A_G_D) Remarkable work done by Thomas Borrvall from Dynamore Nordic and Tobias Erhart from Dynamore The development was customer driven by Mercedes-Benz and Honda-NA The orthotropic failure model in M_A_G_D is to our knowledge the only failure model that allows for directional dependency of the failure strain upon the state of stress Procedures need to be developed for the data generation of aluminium extrusions to speed up the process Magnesium extrusions could be even more challenging
54 The (preliminary) end
INCREASING RUPTURE PREDICTABILITY FOR ALUMINUM
1 INCREASING RUPTURE PREDICTABILITY FOR ALUMINUM Influence of anisotropy Daniel Riemensperger, Adam Opel AG Paul Du Bois, PDB 2 www.opel.com CONTENT Introduction/motivation Isotropic & anisotropic material
More informationShort Introduction of a New Generalized Damage Model
Short Introduction of a New Generalized Damage Model Tobias Erhart 1, Paul Du Bois 2, Filipe Andrade 1 1 Dynamore GmbH, Stuttgart 2 Consulting Engineer, Offenbach 1 Introduction Sound numerical models
More informationWORKSHOP Introduction to material characterization
WORKSHOP Introduction to material characterization Martin Helbig, David Koch DYNAmore GmbH, Germany 15 th German LS-DYNA Forum 2018 Bamberg, October 15, 2018 1 Motivation Example of material usage in a
More informationCharacterizations of Aluminum Alloy Sheet Materials Numisheet 2005
Characterizations of Aluminum Alloy Sheet Materials Numisheet 25 John C. Brem, Frederic Barlat, Robert E. Dick, and Jeong-Whan Yoon Alcoa Technical Center, PA, 1569-1, USA Abstract. This report reproduces
More informationGISSMO Material Modeling with a sophisticated Failure Criteria
GISSMO Material Modeling with a sophisticated Failure Criteria André Haufe, Paul DuBois, Frieder Neukamm, Markus Feucht Dynamore GmbH Industriestraße 2 70565 Stuttgart http://www.dynamore.de Motivation
More informationIDENTIFICATION OF SHEET METAL PLASTIC ANISOTROPY, AND OPTIMIZATION OF INITIAL BLANK SHAPE IN DEEP DRAWING
THE ANNALS OF DUNĂREA DE JOS UNIVERSITY OF GALAŢI FASCICLE V, TECHNOLOGIES IN MACHINE BUILDING, ISSN - 66, 9 IDENTIFICATION OF SHEET METAL PLASTIC ANISOTROPY, AND OPTIMIZATION OF INITIAL BLANK SHAPE IN
More informationLecture #7: Basic Notions of Fracture Mechanics Ductile Fracture
Lecture #7: Basic Notions of Fracture Mechanics Ductile Fracture by Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling of Materials in Manufacturing
More informationMaterial parameter identification for the numerical simulation of deep-drawing drawing of aluminium alloys
Material parameter identification for the numerical simulation of deep-drawing drawing of aluminium alloys B.M. Chaparro, J.M. Antunes, A.J. Baptista, D.M. Rodrigues,, J.V. Fernandes, L.F. Menezes M2D'26
More informationA Semianalytical Model for the Simulation of Polymers
A Semianalytical Model for the Simulation of Polymers Paul Du Bois 3, Stefan Kolling 1, Markus Feucht 1 & André Haufe 2 1 DaimlerChrysler AG, Sindelfingen, Germany 2 Dynamore GmbH, Stuttgart, Germany 3
More informationA rate-dependent Hosford-Coulomb model for predicting ductile fracture at high strain rates
EPJ Web of Conferences 94, 01080 (2015) DOI: 10.1051/epjconf/20159401080 c Owned by the authors, published by EDP Sciences, 2015 A rate-dependent Hosford-Coulomb model for predicting ductile fracture at
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationAnisotropic modeling of short fibers reinforced thermoplastics materials with LS-DYNA
Anisotropic modeling of short fibers reinforced thermoplastics materials with LS-DYNA Alexandre Hatt 1 1 Faurecia Automotive Seating, Simplified Limited Liability Company 1 Abstract / Summary Polymer thermoplastics
More informationExperience from using recently implemented enhancements for Material 36 in LS-DYNA 971 performing a virtual tensile test. Authors: Correspondence:
Experience from using recently implemented enhancements for Material 36 in LS-DYNA 971 performing a virtual tensile test Authors: Michael Fleischer BMW Group Thomas Borrvall Engineering Research Nordic
More informationNumerical simulation of sheet metal forming processes using a new yield criterion
Key Engineering Materials Vol. 344 (007) pp. 833-840 online at http://www.scientific.net (007) Trans Tech Publications, Switzerland Numerical simulation of sheet metal forming processes using a new yield
More information*MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials
14 th International LS-DYNA Users Conference Session: Constitutive Modeling *MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials Jesper Karlsson 1, Mikael Schill 1, Johan Tryding
More informationEFFECT OF ANISOTROPIC YIELD CRITERION ON THE SPRINGBACK IN PLANE STRAIN PURE BENDING
EFFECT OF ANISOTROPIC YIELD CRITERION ON THE SPRINGBACK IN PLANE STRAIN PURE BENDING F.V. Grechnikov 1,, Ya.A. Erisov 1, S.E. Alexandrov 3 1 Samara National Research University, Samara, Russia Samara Scientific
More informationLecture #8: Ductile Fracture (Theory & Experiments)
Lecture #8: Ductile Fracture (Theory & Experiments) by Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling of Materials in Manufacturing 2015 1 1 1 Ductile
More informationAn orthotropic damage model for crash simulation of composites
High Performance Structures and Materials III 511 An orthotropic damage model for crash simulation of composites W. Wang 1, F. H. M. Swartjes 1 & M. D. Gan 1 BU Automotive Centre of Lightweight Structures
More informationAnisotropy and Failure Modeling for Nonlinear Strain. Paths and Its Application to Rigid Packaging. Robert Earl Dick (B.S.C.E, M.S.C.
Anisotropy and Failure Modeling for Nonlinear Strain Paths and Its Application to Rigid Packaging by Robert Earl Dick (B.S.C.E, M.S.C.E) Advisor: Jeong Whan Yoon Submitted in fulfillment of the requirements
More informationRecent Developments in Damage and Failure Modeling with LS-DYNA
Recent evelopments in amage and Failure Modeling with LS-YNA r. André Haufe ynamore GmbH Frieder Neukamm, r. Markus Feucht aimler AG Paul ubois Consultant r. Thomas Borvall ERAB 1 Technological challenges
More informationLecture #10: Anisotropic plasticity Crashworthiness Basics of shell elements
Lecture #10: 151-0735: Dynamic behavior of materials and structures Anisotropic plasticity Crashworthiness Basics of shell elements by Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering,
More informationSTRESS UPDATE ALGORITHM FOR NON-ASSOCIATED FLOW METAL PLASTICITY
STRESS UPDATE ALGORITHM FOR NON-ASSOCIATED FLOW METAL PLASTICITY Mohsen Safaei 1, a, Wim De Waele 1,b 1 Laboratorium Soete, Department of Mechanical Construction and Production, Ghent University, Technologiepark
More informationInverse identification of plastic material behavior using. multi-scale virtual experiments
Inverse identification of plastic material behavior using multi-scale virtual experiments Debruyne 1 D., Coppieters 1 S., Wang 1 Y., Eyckens 2 P., Kuwabara 3 T., Van Bael 2 A. and Van Houtte 2 P. 1 Department
More informationA Simple and Accurate Elastoplastic Model Dependent on the Third Invariant and Applied to a Wide Range of Stress Triaxiality
A Simple and Accurate Elastoplastic Model Dependent on the Third Invariant and Applied to a Wide Range of Stress Triaxiality Lucival Malcher Department of Mechanical Engineering Faculty of Tecnology, University
More informationMultiaxial Fatigue. Professor Darrell F. Socie. Department of Mechanical Science and Engineering University of Illinois at Urbana-Champaign
Multiaxial Fatigue Professor Darrell F. Socie Department of Mechanical Science and Engineering University of Illinois at Urbana-Champaign 2001-2011 Darrell Socie, All Rights Reserved Contact Information
More informationChapter 6: Plastic Theory
OHP Mechanical Properties of Materials Chapter 6: Plastic Theory Prof. Wenjea J. Tseng 曾文甲 Department of Materials Engineering National Chung Hsing University wenjea@dragon.nchu.edu.tw Reference: W. F.
More information14. LS-DYNA Forum 2016
14. LS-DYNA Forum 2016 A Novel Approach to Model Laminated Glass R. Böhm, A. Haufe, A. Erhart DYNAmore GmbH Stuttgart 1 Content Introduction and Motivation Common approach to model laminated glass New
More informationEvaluation of Material Models to Predict Material Failure in LS-DYNA
Department of Construction Sciences Solid Mechanics ISRN LUTFD2/TFHF-16/5209-SE(1-55) Evaluation of Material Models to Predict Material Failure in LS-DYNA Master s Dissertation by Hjalmar Sandberg Oscar
More informationPredicting Failure of Multiangle Composite Laminates
Predicting Failure of Multiangle Composite Laminates Preliminary discussion (not in textbook): Micromechanics failure analyses vs Macromechanics failure analyses Fiber Architecture of Some Common Composite
More informationMechanical Properties of Materials
Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of
More informationModule III - Macro-mechanics of Lamina. Lecture 23. Macro-Mechanics of Lamina
Module III - Macro-mechanics of Lamina Lecture 23 Macro-Mechanics of Lamina For better understanding of the macromechanics of lamina, the knowledge of the material properties in essential. Therefore, the
More informationCharacterization of Anisotropic Plasticity in Material Systems Using Modified Indentation-Based Techniques
Characterization of Anisotropic Plasticity in Material Systems Using Modified Indentation-Based Techniques A Dissertation Presented by Salmon Masbooghi Kalkhoran to The Department of Mechanical & Industrial
More informationON THE USE OF HOMOGENEOUS POLYNOMIALS TO DEVELOP ANISOTROPIC YIELD FUNCTIONS WITH APPLICATIONS TO SHEET FORMING
ON THE USE OF HOMOGENEOUS POLYNOMIALS TO DEVELOP ANISOTROPIC YIELD FUNCTIONS WITH APPLICATIONS TO SHEET FORMING By STEFAN C. SOARE A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA
More informationPlane Strain Test for Metal Sheet Characterization
Plane Strain Test for Metal Sheet Characterization Paulo Flores 1, Felix Bonnet 2 and Anne-Marie Habraken 3 1 DIM, University of Concepción, Edmundo Larenas 270, Concepción, Chile 2 ENS - Cachan, Avenue
More informationAn overview of Carbon Fiber modeling in LS-DYNA. John Zhao October 23 th 2017
An overview of Carbon Fiber modeling in LS-DYNA John Zhao zhao@lstc.com October 23 th 2017 Outline Manufacturing of Carbon Fiber Compression molding *MAT_277 & 278 *MAT_293 *MAT_249 Resin transform molding
More informationFatigue and Fracture
Fatigue and Fracture Multiaxial Fatigue Professor Darrell F. Socie Mechanical Science and Engineering University of Illinois 2004-2013 Darrell Socie, All Rights Reserved When is Multiaxial Fatigue Important?
More informationA modified Burzynski criterion for anisotropic pressure-dependent materials
Sādhanā Vol. 4, No., January 7, pp. 95 9 DOI.7/s46-6-576-6 Indian Academy of Sciences A modified Burzynski criterion for anisotropic pressure-dependent materials FARZAD MOAYYEDIAN and MEHRAN KADKHODAYAN,
More informationFracture mechanics fundamentals. Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design
Fracture mechanics fundamentals Stress at a notch Stress at a crack Stress intensity factors Fracture mechanics based design Failure modes Failure can occur in a number of modes: - plastic deformation
More informationExperimentally Calibrating Cohesive Zone Models for Structural Automotive Adhesives
Experimentally Calibrating Cohesive Zone Models for Structural Automotive Adhesives Mark Oliver October 19, 2016 Adhesives and Sealants Council Fall Convention contact@veryst.com www.veryst.com Outline
More informationCard Variable MID RO E PR ECC QH0 FT FC. Type A8 F F F F F F F. Default none none none 0.2 AUTO 0.3 none none
Note: This is an extended description of MAT_273 input provided by Peter Grassl It contains additional guidance on the choice of input parameters beyond the description in the official LS-DYNA manual Last
More informationDevelopment, implementation and Validation of 3-D Failure Model for Aluminium 2024 for High Speed Impact Applications
Development, implementation and Validation of -D Failure Model for Aluminium 2024 for High Speed Impact Applications Paul Du Bois, Murat Buyuk, Jeanne He, Steve Kan NCAC-GWU th International LS-DYNA Users
More informationNUMERICAL SIMULATION OF DAMAGE IN THERMOPLASTIC COMPOSITE MATERIALS
5 th European LS-DYNA Users Conference Composites NUMERICAL SIMULATION OF DAMAGE IN THERMOPLASTIC COMPOSITE MATERIALS Kevin Brown 1, Richard Brooks, Nicholas Warrior School of Mechanical, Materials and
More informationFCP Short Course. Ductile and Brittle Fracture. Stephen D. Downing. Mechanical Science and Engineering
FCP Short Course Ductile and Brittle Fracture Stephen D. Downing Mechanical Science and Engineering 001-015 University of Illinois Board of Trustees, All Rights Reserved Agenda Limit theorems Plane Stress
More informationResearch Article Process of Identifying Stress Fields from Strain Fields in the Specimen with a Hole
ISRN Materials Science Volume 213, Article ID 53669, 8 pages http://dx.doi.org/1.1155/213/53669 Research Article Process of Identifying Stress Fields from Strain Fields in the Specimen with a Hole Michaela
More informationPrincipal Stresses, Yielding Criteria, wall structures
Principal Stresses, Yielding Criteria, St i thi Stresses in thin wall structures Introduction The most general state of stress at a point may be represented by 6 components, x, y, z τ xy, τ yz, τ zx normal
More informationModelling of ductile failure in metal forming
Modelling of ductile failure in metal forming H.H. Wisselink, J. Huetink Materials Innovation Institute (M2i) / University of Twente, Enschede, The Netherlands Summary: Damage and fracture are important
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationBone Tissue Mechanics
Bone Tissue Mechanics João Folgado Paulo R. Fernandes Instituto Superior Técnico, 2016 PART 1 and 2 Introduction The objective of this course is to study basic concepts on hard tissue mechanics. Hard tissue
More informationProblem " Â F y = 0. ) R A + 2R B + R C = 200 kn ) 2R A + 2R B = 200 kn [using symmetry R A = R C ] ) R A + R B = 100 kn
Problem 0. Three cables are attached as shown. Determine the reactions in the supports. Assume R B as redundant. Also, L AD L CD cos 60 m m. uation of uilibrium: + " Â F y 0 ) R A cos 60 + R B + R C cos
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationDigimat material model for short fiber reinforced plastics at Volvo Car Corporation
Digimat material model for short fiber reinforced plastics at Volvo Car Corporation Mats Landervik DYNAmore Nordic Johan Jergeus Volvo Car Corporation June 2, 2015 26TH SICOMP CONFERENCE, GOTHENBURG 1
More informationComputational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem
Computational Inelasticity FHLN05 Assignment 2018 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than November
More informationMECHANICS OF MATERIALS. EQUATIONS AND THEOREMS
1 MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS Version 2011-01-14 Stress tensor Definition of traction vector (1) Cauchy theorem (2) Equilibrium (3) Invariants (4) (5) (6) or, written in terms of principal
More informationOutline. Tension test of stainless steel type 301. Forming including the TRIP-effect
Outline Erweiterte Kontaktformulierung, User-Reibmodell und thermische Analyse Simulationsmöglichkeiten zur Werkzeugauslegung mit LS-DYNA Workshop Simulation in der Umformtechnik Institut für Statik und
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationTHE HYDRAULIC BULGE TEST AND ITS IMPORTANCE FOR THE VEGTER YIELD CRITERION
THE HYDRAULIC BULGE TEST AND ITS IMPORTANCE FOR THE VEGTER YIELD CRITERION Jiří SOBOTKA a, Pavel SOLFRONK a, Pavel DOUBEK a, Lukáš ZUZÁNEK a a TECHNICAL UNIVERSITY OF LIBEREC, Studentská, 461 17 Liberec
More informationExperiments and modeling of edge fracture for an AHSS sheet
Experiments and modeling of edge fracture for an AHSS sheet The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published
More informationInternational Journal of Solids and Structures
International Journal of Solids and Structures 49 (2012) 2863 2880 Contents lists available at SciVerse ScienceDirect International Journal of Solids and Structures journal homepage: www.elsevier.com/locate/ijsolstr
More informationAnalytical and numerical analysis of sheet metal instability using a stress based criterion
Available online at www.sciencedirect.com International Journal of Solids and Structures 45 (28) 242 255 www.elsevier.com/locate/ijsolstr Analytical and numerical analysis of sheet metal instability using
More informationMATERIAL MODELS FOR THE
Project work report MATERIAL MODELS FOR THE DEFORMATION AND FAILURE OF HIGH STRENGTH STEEL Katia BONI August 2008 to January 2009 Supervisor: Larsgunnar NILSSON 2 Contents 1 Introduction 6 1.1 Background.........................................
More informationEnhancement of Spot Weld Modeling using MAT 100 DAI
7. LS-DYNA Anwenderforum, Bamberg 2008 Crash II - Verbindungstechnik, Versagen Enhancement of Spot Weld Modeling using MAT 100 DAI Falko Seeger a, Markus Feucht a, George Dumitru b, Tobias Graf b a Daimler
More informationImplementation of Orthotropic Elasto-Plastic model in MOOSE and Pellet-Cladding Interaction Simulations
Implementation of Orthotropic Elasto-Plastic model in MOOSE and Pellet-Cladding Interaction Simulations By Axel Seoane Mentors: Dr. Piyush Sabharwall and Dr. Pritam Chakraborty April 12, 2017 Copyright
More informationAnisotropy for Metals in Non-Associated Flow Plasticity
Anisotropy for Metals in Non-Associated Flow Plasticity by Mariana Paulino Santos MEng Submitted in fulfilment of the requirements for the degree of Doctor of Philosophy Deakin University August, 2015
More informationINFLUENCE OF CONSTITUTIVE EQUATIONS ON THE ACCURACY OF PREDICTION IN SHEET METAL FORMING SIMULATION
Septemer 1-5, 008 Interlaken, Switzerland INFLUENCE OF CONSTITUTIVE EQUATIONS ON THE ACCURACY OF PREDICTION IN SHEET METAL FORMING SIMULATION Dorel Banaic 1*, Dan Sorin Comsa 1, Matthias Sester, Mike Selig,
More informationLinear and Nonlinear Analysis of High Dynamic Impact Events with Creo Simulate and Abaqus/Explicit
Linear and Nonlinear Analysis of High Dynamic Impact Events with Creo Simulate and Abaqus/Explicit SAXSIM 2015, TU Chemnitz Dr.-Ing. Roland Jakel At a glance: Altran, a global leader 20+ Countries 23,000+
More informationExperiments and Numerical Simulations on Stress-State-Dependence of Ductile Damage Criteria
Experiments and Numerical Simulations on Stress-State-Dependence of Ductile Damage Criteria Michael Brünig, Steffen Gerke and Daniel Brenner Abstract The paper deals with a series of new experiments and
More informationMECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso
MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Thermodynamics Derivation Hooke s Law: Anisotropic Elasticity
More informationMechanics of Materials Primer
Mechanics of Materials rimer Notation: A = area (net = with holes, bearing = in contact, etc...) b = total width of material at a horizontal section d = diameter of a hole D = symbol for diameter E = modulus
More informationFINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS
FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS Ever J. Barbero Department of Mechanical and Aerospace Engineering West Virginia University USA CRC Press Taylor &.Francis Group Boca Raton London New York
More informationAdvanced Numerical Study of the Effects of Road Foundations on Pavement Performance
Advanced Numerical Study of the Effects of Road Foundations on Pavement Performance X. Liu Section of Structural Mechanics, Faculty of Civil Engineering and Geosciences, Delft University of Technology,
More informationIMPACT STRENGTH AND RESPONSE BEHAVIOR OF CFRP GUARDER BELT FOR SIDE COLLISION OF AUTOMOBILES
IMPACT STRENGTH AND RESPONSE BEHAVIOR OF CFRP GUARDER BELT FOR SIDE COLLISION OF AUTOMOBILES AOKI Y.*, BEN G.**, KIM H. S.* *College of Science & Technology, Nihon University **College of Industrial Technology,
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationCrash Safety of Lithium-Ion Batteries for Electric Vehicles
Crash Safety of Lithium-Ion Batteries for Electric Vehicles Elham Sahraei Esfahani Research Scientist, MIT Assistant Professor, GMU Vehicle Crashes and Testing Battery Packs in T Configuration Chevy Volt
More informationMechanical analysis of timber connection using 3D finite element model
Mechanical analysis of timber connection using 3D finite element model Bohan XU Ph.D Student Civil Engineering Laboratory (CUST) Clermont-Ferrand, France Mustapha TAAZOUNT Dr-Ing Civil Engineering Laboratory
More informationModelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases
Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Georg Haasemann Conrad Kloß 1 AIMCAL Conference 2016 MOTIVATION Wrinkles in web handling system Loss of
More informationFormability assessment of a cup drawing under complex nonlinear strain paths
Formability assessment of a cup drawing under complex nonlinear strain paths P.D. Barros 1, R.R. Amaral 2, M.C. Oliveira 1, A.D. Santos 2,3, J.L. Alves 4 and L.F. Menezes 1 1 CEMMPRE, Centre for Mechanical
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationCrash and Vibration Analysis of rotors in a Roots Vacuum Booster
Crash and Vibration Analysis of rotors in a Roots Vacuum Booster Manuel Roth 1, Stefan Kolling 2 1 Pfeiffer Vacuum GmbH, R&D Backing Pumps, Asslar, Germany 2 Giessen University of Applied Sciences, Laboratory
More informationNUMERICAL SIMULATIONS OF NAKAZIMA FORMABILITY TESTS WITH PREDICTION OF FAILURE
NUMERICAL SIMULATIONS OF NAKAZIMA FORMABILITY TESTS WITH PREDICTION OF FAILURE DMYTRO LUMELSKYJ 1, JERZY ROJEK 1, MAREK TKOCZ 2 Abstract. This paper presents results of numerical simulations of the Nakazima
More informationINVERSE METHOD FOR FLOW STRESS PARAMETERS IDENTIFICATION OF TUBE BULGE HYDROFORMING CONSIDERING ANISOTROPY
7 th EUROMECH Solid Mechanics Conference J. Ambrósio et.al. (eds.) Lisbon, Portugal, September 7-11, 2009 INVERSE METHOD FOR FLOW STRESS PARAMETERS IDENTIFICATION OF TUBE BULGE HYDROFORMING CONSIDERING
More informationEnhancing Prediction Accuracy In Sift Theory
18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS Enhancing Prediction Accuracy In Sift Theory J. Wang 1 *, W. K. Chiu 1 Defence Science and Technology Organisation, Fishermans Bend, Australia, Department
More informationCrashworthy Design of Composite Structures Using CAE Process Chain
0 th European LS-DYNA Conference 205, Würzburg, Germany Crashworthy Design of Composite Structures Using CAE Process Chain Madhukar Chatiri, Thorsten Schuetz 2, Anton Matzenmiller 3 CADFEM GmbH, Grafing
More informationAnalytical and numerical prediction of localized necking in sheet metal forming processes
University of New Hampshire University of New Hampshire Scholars' Repository Doctoral Dissertations Student Scholarship Fall 2011 Analytical and numerical prediction of localized necking in sheet metal
More informationGATE SOLUTIONS E N G I N E E R I N G
GATE SOLUTIONS C I V I L E N G I N E E R I N G From (1987-018) Office : F-16, (Lower Basement), Katwaria Sarai, New Delhi-110016 Phone : 011-65064 Mobile : 81309090, 9711853908 E-mail: info@iesmasterpublications.com,
More informationStress, Strain, Mohr s Circle
Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected
More informationUnit IV State of stress in Three Dimensions
Unit IV State of stress in Three Dimensions State of stress in Three Dimensions References Punmia B.C.,"Theory of Structures" (SMTS) Vol II, Laxmi Publishing Pvt Ltd, New Delhi 2004. Rattan.S.S., "Strength
More informationThe Significance of the Production Process of FRP Parts for the Performance in Crashworthiness
14 th International LS-DYNA Users Conference Session: Simulation The Significance of the Production Process of FRP Parts for the Performance in Crashworthiness Christian Liebold, Andrea Erhart, André Haufe
More informationLecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity
Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling
More informationDD3MAT - a code for yield criteria anisotropy parameters identification.
Journal of Physics: onference Series PAPER OPEN AESS DD3MA - a code for yield criteria anisotropy parameters identification. o cite this article: P D Barros et al 216 J. Phys.: onf. Ser. 734 3253 View
More informationPressure Vessels Stresses Under Combined Loads Yield Criteria for Ductile Materials and Fracture Criteria for Brittle Materials
Pressure Vessels Stresses Under Combined Loads Yield Criteria for Ductile Materials and Fracture Criteria for Brittle Materials Pressure Vessels: In the previous lectures we have discussed elements subjected
More information20. Rheology & Linear Elasticity
I Main Topics A Rheology: Macroscopic deformation behavior B Linear elasticity for homogeneous isotropic materials 10/29/18 GG303 1 Viscous (fluid) Behavior http://manoa.hawaii.edu/graduate/content/slide-lava
More informationMULTIAXIAL ASSESSMENT METHOD FOR FATIGUE CALCULATIONS IN COMPOSITE COMPONENTS
MULTIAXIAL ASSESSMENT METHOD FOR FATIGUE CALCULATIONS IN COMPOSITE COMPONENTS Peter Heyes HBM-nCode Abstract Fatigue analysis methods are well-established as part of the product development process for
More informationOn Springback Prediction In Stamping Of AHSS BIW Components Utilizing Advanced Material Models
On Springback Prediction In Stamping Of AHSS BIW Components Utilizing Advanced Material Models Ming F. Shi and Alex A. Konieczny United States Steel Corporation Introduction Origin of Springback AHSS Springback
More informationFracture Mechanics, Damage and Fatigue: Composites
University of Liège Aerospace & Mechanical Engineering Fracture Mechanics, Damage and Fatigue: Composites Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationSheet metal forming with six components of strain
Sheet metal forming with six components of strain Julian Allwood Micromechanics Seminar Friday 8 May 2009 Department of Engineering University of Cambridge jma42@cam.ac.uk 1 A question from the 1960 s
More informationCrash and Impact Simulation of Composite Structures by Using CAE Process Chain
Crash and Impact Simulation of Composite Structures by Using CAE Process Chain Madhukar Chatiri 1, Thorsten Schütz 2, Anton Matzenmiller 3, Ulrich Stelzmann 1 1 CADFEM GmbH, Grafing/Munich, Germany, mchatiri@cadfem.de
More informationName (Print) ME Mechanics of Materials Exam # 2 Date: March 29, 2016 Time: 8:00 10:00 PM - Location: PHYS 114
Name (Print) (Last) (First) Instructions: ME 323 - Mechanics of Materials Exam # 2 Date: March 29, 2016 Time: 8:00 10:00 PM - Location: PHYS 114 Circle your lecturer s name and your class meeting time.
More informationComposite Materials 261 and 262
Composite Materials 261 and 262 Stefan Hartmann 1 David Moncayo 2 1 DYNAmore GmbH, Stuttgart, Germany 2 Daimler AG, Sindelfingen, Germany 11. LS-DYNA Forum 2012, 9. - 10. Oktober 2012, Ulm 1 Outline Introduction
More information