ACCESS Beiratssitzung November Dr. G. Laschet and L. Haas ACCESS e.v. Aachen, Germany. FENET Workshop Copenhagen
|
|
- Barrie Hoover
- 5 years ago
- Views:
Transcription
1 Coupling FE Contact and Heat Transfer Analysis in Casting Simulations ACCESS Beiratssitzung November 2001 Dr. G. Laschet and L. Haas ACCESS e.v. Aachen, Germany FENET Workshop Copenhagen
2 About ACCESS e.v. Created in 1986 as Center for Solidification in Space Private research center annexed to RWTH Aachen Simulation and optimization of casting processes 45 employees thereunder 35 scientists Annual turnover : 5 million Euros
3 Coupled Process Modelling (FE/CV Formulation)
4 Thermomechanical Analysis During Solidification with CASTS Special items in stress-strain simulation during casting: phase change for alloys and pure substance weak or strong coupling between thermal and stress-strain analysis coupling between heat transfer and contact or gap formation during shrinkage constitutive models for metal and mould complex 3-D geometries (e.g. investment casting)
5 Thermomechanical Constitutive Models core and mould: isotropic, thermoelastic behaviour thermoviscoplasticity (only 12.4 version) metal: isotropic, thermoelastoplastic behaviour anisotropic elastic behaviour for SX or DS superalloys (12.4 version) isotropic thermoviscoplastic behaviour (Norton law, Robinson model), (12.4 version) liquid: quasi-incompressible elastic behaviour at the solidification front: either stress free or metallostatic pressure liquid / solid density variation
6 Heat Transfer Coefficient as Function of Gap or Contact (1) gap: radiative heat transfer: forming gap is assumed to be parallel and plane h rad 4 4 ( T T ) 1 εc 1 ε ( T T ) + 1 conductive heat transfer term is given by: = c σ m c with I 0 : mean free path of gas molecules m m λg hgas = gn + I0 contact: h con function of the contact pressure: γ p hcon = c * with 0.6 γ 1 p0 c* and p 0 determined experimentally in order to get a smooth transition between contact and gap situation
7 Heat Transfer Coefficient as Function of Gap or Contact (2) hgas = λgas ( gap + l ) 0 hrad = σ 2 2 ( T + T )( T + T ) c m c εc εm m
8 3-D Master-slave Contact Algorithm for Casting Simulation (1) normal contact calculated between two elastoplastic bodies undergoing small deformations friction is actually neglected non-penetration condition: penalty method is implemented single and two-pass options are available augmented lagrangian method is under development specific features for casting simulations: mould surface is specified as master and casting surface as slave if metal is liquid: tied contact between metal and mould master / slave algorithm activated for a slave node if all connecting nodes in the metal are solidified
9 3-D Master-Slave Contact Algorithm (2) normal gap definition: g n = u n -x n g n 0: contact g n > 0: no contact contact pressure: penalty method: t n = C (E) <g n > m m 2 C(E): penalty factor augmented lagrangian method: t n = λ n + C(E) <g n > difficulties: description of master surface: definition of macro-planes or surface smoothing update of the master surface for thin bodies 3-D closest point projection algorithm: crucial near corners, sharp angles
10 Definition of Master Macro-Planes 2-D cut M 1, M 7 : centroid of the element surface n 1, n 7 : normal to the element surface loop over element faces master surface: same macro-plane l if: ( n,n ) arccos i 1,i 1 M1 Mi n1,i 1 < dcr <αcr with default values: α cr = 8.1 or cos α r = 0.99 ; d cr = 100 µm i 1 nk n k 1 1,i 1= = i 1 finally, shift from M 1 to M c (= centroid of these elements) assumption: small deformation macro-planes stay planar during deformation
11 Update of the Master Surface crucial for cast parts with thin walls (mainly in investment casting) definition of a critical updating: ½ of the thinnest section if u > critical updating smaller increments two-pass algorithm
12 Closest Point Projection Algorithm Classical normal projection on master macro-planes l (l= 1,2): if projection point (y 1, y 3 ) is in elements l contact element generation for the plane having smallest gap if Sy 2 -Sy 2' < prec generation of 2 contact elements if projection is in the unreachable domain (y 4, y 4' ): unreachable domain nearest master node (here M1) is adopted as projection point contact elements generated for all planes to which this node belong to if master points are at the same distance from the slave S 5 : the point with the greatest number of connecting planes is selected (here M 5 )
13 Determination of the Normal Gap Loop over the Newton-Raphson iterations at time step t it=1: it>1: update of the master surface for each slave node: compute the new initial gap by closest point projection loop over the slave nodes if contact is active closest point projection else: determination of possible intersection points nearest point is selected compute normal gap for the corresponding master element if no intersection point exists closest point projection Why intersection point? contact with only plane 1 will be considered
14 Crisfield Contact Patch Test Evaluation of the performance of the contact model in term of accuracy and convergence capability (pure mechanical test) 1 y interface x 1 incompatible meshes Analytical stress solution (plane strain assumption) σ x = 0, σ y = E / (1-ν 2 ), τ xy = 0, σ z = νσ y
15 Crisfield Patch Test: Results for a Straight Contact Surface Stress σ y Stress σ x Contact pressure (deformation x10) Simulation results: analytical solution σ y = MPa verified with an error of 0.1% analytical solution σ x = 0 is verified normal contact pressure equal to the uniform compression in y-direction
16 Second Crisfield Patch Test: Curved Contact Surfaces initial penetration Contact pressure (CASTS) Solution σ x [MPa] σ y [MPa] σ z [MPa] τ xy [MPa] min max min max min max min max Analytical CASTS ABAQUS 84, error in σ y Crisfield : min 2.3% ; max 3.0% CASTS: min 13.4% ; max 49.4% ABAQUS: min 2.8% ; max 44.9% both programs show poor prediction of the analytical solution ABAQUS a little more accurate than CASTS
17 Simulation of a Dummy Blade Subjected to a Thermal Cycle Geometry of the dummy blade Applied temperature profile mould: Al 2 O 3 casting: AlSi 7 Mg model: 2680 nodes elements Shrinkage and dilatation of the dummy blade leading to gap formation and contact between metal and mould N.B.: 15 master macro-planes are defined
18 Dummy Blade Results: Evolution of the Gap t= 1200 s t= 3700 s t= 5000 s deformation scale: 20
19 Dummy Blade Results: Prediction of the Contact Pressure t= 1200 s t= 5000 s t= 3700 s deformation scale: 20
20 Dummy Blade Results: Evolution of the Heat Transfer Coefficients t= 60 s t= 1200 s t= 3700 s
21 Dimensional Stability of an Investment Casting Part Specified size = 329,0 ± 2,4 mm VDG-Merkblatt P690 (March 1999) Objective: scaling of the CAD-model dimension of cast part within allowed tolerance
22 Description of the Process CAD-model (CATIA): Wax injection Mould production Pouring Cast part Rapid-Prototyping 3D-measurement, comparison CAD-model 3D-data of the cast part FE-model Simulation
23 Volume Model and FE Mesh Discretized geometry (69335 nodes, elements)
24 Calculation of the Shell Mould Pre-heating Temperature calculation specifying radiation with view factors t = 840 s Restart conditions for the second calculation (the mould filling phase is neglected)
25 Thermomechanical Simulation of the Solidification Process (1) Evolution of the temperature distribution in the cast part t = 1000 s
26 Thermomechanical Simulation of the Solidification Process (2) Apparition of contact in locally solidified regions t = 995 s t = 1000 s Contact pressure (deformation scale: 10)
27 Thermomechanical Simulation of the Solidification Process (3) t = 1000 s t = 1005 s Contact pressure (deformation scale: 2) at time step 1005 s and subsequent : no convergence of N-R iterations definition of a new benchmark test for thin investment casting parts implementation of a more robust augmented lagrangian algorithm
28 Outlook: New Thermomechanical Benchmark for Thin Cast Parts Geometry of crossed shells Tetrahedron mesh (2408 nodes, 1198 elements) Hexahedron mesh (1854 nodes, 2862 elements) ceramic mould (Al 2 O 3 ) of 3 mm is generated on the metal (AlSi 7 Mg) surface I-DEAS universal file is available temperature dependence of thermophysical data of mould and metal is given
29 Boundary and Initial Conditions of the Thermomechanical Benchmark mechanical BC: bottom shell face fixed in normal direction fixation of the rigid body modes at this face thermal BC: bottom shell face isolated free radiation + convection (ε = 0.5, α=10-3 W/cm 2 K, T room = 20 C) 3 problems can be specified as function of initial temperature: a) T met = 620 C >T liq = 613 C and T mould = 550 C > T sol = 548 C test for PROCASTS, THERCASTS, CASTS b) T met = 500 C and T mould = 430 C: challenging test for FE codes c) preheating of the shell in a cubic furnace (20 x 20 x 20 cm 3 ) at 750 C by radiation (ε = 0.9) during 800 s non uniform initial mould temperature due to shadowing
30 Thermomechanical Benchmark: First Results for Case A evolution of the gap prediction (deformation x 5) t= 65 s temperature at t= 65 s torsion of the body t= 70 s 37 master macro-planes are specified gap / contact oscillations
Coupled Thermomechanical Contact Problems
Coupled Thermomechanical Contact Problems Computational Modeling of Solidification Processes C. Agelet de Saracibar, M. Chiumenti, M. Cervera ETS Ingenieros de Caminos, Canales y Puertos, Barcelona, UPC
More informationThermo Mechanical Analysis of AV1 Diesel Engine Piston using FEM
Journal of Advanced Engineering Research ISSN: 2393-8447 Volume 2, Issue 1, 2015, pp.23-28 Thermo Mechanical Analysis of AV1 Diesel Engine Piston using FEM Subodh Kumar Sharma 1, *, P. K. Saini 2, N. K.
More informationContinuum Mechanics and Theory of Materials
Peter Haupt Continuum Mechanics and Theory of Materials Translated from German by Joan A. Kurth Second Edition With 91 Figures, Springer Contents Introduction 1 1 Kinematics 7 1. 1 Material Bodies / 7
More informationNonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder. Richard Kennison, Race-Tec
Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder Richard Kennison, Race-Tec Agenda Race-Tec Overview Accumulator Experimental Testing Material Testing Numerical Analysis: 1. Linear Buckling
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline
More informationProject. First Saved Monday, June 27, 2011 Last Saved Wednesday, June 29, 2011 Product Version 13.0 Release
Project First Saved Monday, June 27, 2011 Last Saved Wednesday, June 29, 2011 Product Version 13.0 Release Contents Units Model (A4, B4) o Geometry! Solid Bodies! Parts! Parts! Body Groups! Parts! Parts
More informationLingopti project: Semi-Continuous Casting Process of Copper-Nickel Alloys
Lingopti project: Semi-Continuous Casting Process of Copper-Nickel Alloys E.Pecquet 1, F.Pascon 1, J.Lecomte-Beckers 2, A-M.Habraken 1 University of Liège, 1 Department M&S, 2 Service MMS Chemin des Chevreuils
More informationHeat Transfer Modeling using ANSYS FLUENT
Lecture 2 - Conduction Heat Transfer 14.5 Release Heat Transfer Modeling using ANSYS FLUENT 2013 ANSYS, Inc. March 28, 2013 1 Release 14.5 Agenda Introduction Energy equation in solids Equation solved
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.
D : SOLID MECHANICS Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown. Q.2 Consider the forces of magnitude F acting on the sides of the regular hexagon having
More informationVULCAN 2000: A FINITE ELEMENT SYSTEM FOR THE SIMULATION OF CASTING PROCESSES C. Agelet de Saracibar, M. Cervera & M. Chiumenti ETS Ingenieros de Camin
VULCAN 2000: A FINITE ELEMENT SYSTEM FOR THE SIMULATION OF CASTING PROCESSES C. AGELET DE SARACIBAR y & M. CERVERA z ETS Ingenieros de Caminos, Canales y Puertos International Center for Numerical Methods
More informationPREPARATION OF HEXCALIBER TESTS AND PRELIMINARY THERMO-MECHANICAL ANALYSES
PREPARATION OF HEXCALIBER TESTS AND PRELIMINARY THERMO-MECHANICAL ANALYSES I. Ricapito a, G. Dell Orco b, P. A. Di Maio c, R. Giammusso c, A. Malavasi a, A. Tincani a, G. Vella c a ENEA C.R. Brasimone
More informationMOOSE Workshop MOOSE-PF/MARMOT Training Mechanics Part Nashville, February 14th Daniel Schwen
www.inl.gov MOOSE Workshop MOOSE-PF/MARMOT Training Mechanics Part Nashville, February 14th Daniel Schwen www.inl.gov III. Mechanics Mechanics Outline 1. Introduction 2. Tensor-based System 1. Compute
More informationNon-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises
Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009
More informationCoupled CFD-FE-Analysis for the Exhaust Manifold of a Diesel Engine
Coupled CFD-FE-Analysis for the Exhaust Manifold of a Diesel Engine Yasar Deger*, Burkhard Simperl*, Luis P. Jimenez** *Sulzer Innotec, Sulzer Markets and Technology Ltd, Winterthur, Switzerland **Guascor
More informationCAEFEM v9.5 Information
CAEFEM v9.5 Information Concurrent Analysis Corporation, 50 Via Ricardo, Thousand Oaks, CA 91320 USA Tel. (805) 375 1060, Fax (805) 375 1061 email: info@caefem.com or support@caefem.com Web: http://www.caefem.com
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 informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
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 informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationLIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS
5 th GRACM International Congress on Computational Mechanics Limassol, 29 June 1 July, 2005 LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS G.A. Drosopoulos I, G.E.
More informationPhase-Field Modeling of Technical Alloy Systems Problems and Solutions
Thermo-Calc User Meeting, Sept 8-9 2011, Aachen Phase-Field Modeling of Technical Alloy Systems Problems and Solutions B. Böttger ACCESS e.v., RWTH Aachen University, Germany Outline 1. Introduction 2.
More informationENGN 2290: Plasticity Computational plasticity in Abaqus
ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element
More informationInitial Stress Calculations
Initial Stress Calculations The following are the initial hand stress calculations conducted during the early stages of the design process. Therefore, some of the material properties as well as dimensions
More informationFLEXIBILITY METHOD FOR INDETERMINATE FRAMES
UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These
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 informationFinite element prediction of the ultimate axial load capacity of V-section band clamps
Journal of Physics: Conference Series Finite element prediction of the ultimate axial load capacity of V-section band clamps To cite this article: S M Barrans and M Muller 2009 J. Phys.: Conf. Ser. 181
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 informationModelling the behaviour of plastics for design under impact
Modelling the behaviour of plastics for design under impact G. Dean and L. Crocker MPP IAG Meeting 6 October 24 Land Rover door trim Loading stages and selected regions Project MPP7.9 Main tasks Tests
More informationImpact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics
Impact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics Dr. A. Johnson DLR Dr. A. K. Pickett ESI GmbH EURO-PAM 99 Impact and Crash Modelling of Composite Structures: A Challenge
More informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.
GTE 2016 Q. 1 Q. 9 carry one mark each. D : SOLID MECHNICS Q.1 single degree of freedom vibrating system has mass of 5 kg, stiffness of 500 N/m and damping coefficient of 100 N-s/m. To make the system
More informationA Demonstrative Computer Session Using ADINA- Nonlinear Analysis
Topic 22 A Demonstrative Computer Session Using ADINA- Nonlinear Analysis Contents: Use of ADINA for elastic-plastic analysis of a plate with a hole Computer laboratory demonstration-part II Selection
More informationAbstract. 1 Introduction
Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut
More informationALGORITHM FOR NON-PROPORTIONAL LOADING IN SEQUENTIALLY LINEAR ANALYSIS
9th International Conference on Fracture Mechanics of Concrete and Concrete Structures FraMCoS-9 Chenjie Yu, P.C.J. Hoogenboom and J.G. Rots DOI 10.21012/FC9.288 ALGORITHM FOR NON-PROPORTIONAL LOADING
More informationGeneral Guidelines for Crash Analysis in LS-DYNA. Suri Bala Jim Day. Copyright Livermore Software Technology Corporation
General Guidelines for Crash Analysis in LS-DYNA Suri Bala Jim Day Copyright Livermore Software Technology Corporation Element Shapes Avoid use of triangular shells, tetrahedrons, pentahedrons whenever
More informationMechanical Engineering Ph.D. Preliminary Qualifying Examination Solid Mechanics February 25, 2002
student personal identification (ID) number on each sheet. Do not write your name on any sheet. #1. A homogeneous, isotropic, linear elastic bar has rectangular cross sectional area A, modulus of elasticity
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems. Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Lecture ST1-19 November, 2015 Institute of Structural Engineering Method of Finite Elements II 1 Constitutive
More informationA FINITE ELEMENT FORMULATION FOR NONLINEAR 3D CONTACT PROBLEMS
A FINITE ELEMENT FORMULATION FOR NONLINEAR 3D CONTACT PROBLEMS Federico J. Cavalieri, Alberto Cardona, Víctor D. Fachinotti and José Risso Centro Internacional de Métodos Computacionales en Ingeniería
More informationModel-independent approaches for the XFEM in fracture mechanics
Model-independent approaches for the XFEM in fracture mechanics Safdar Abbas 1 Alaskar Alizada 2 and Thomas-Peter Fries 2 1 Aachen Institute for Computational Engineering Science (AICES), RWTH Aachen University,
More informationChapter 5: Ball Grid Array (BGA)
Chapter 5: Ball Grid Array (BGA) 5.1 Development of the Models The following sequence of pictures explains schematically how the FE-model of the Ball Grid Array (BGA) was developed. Initially a single
More informationFinite Element Simulation of the Liquid Silicon Oriented Crystallization in a Graphite Mold
Finite Element Simulation of the Liquid Silicon Oriented Crystallization in a Graphite Mold Alexey I. Borovkov Vadim I. Suprun Computational Mechanics Laboratory, St. Petersburg State Technical University,
More informationMECHANICS OF MATERIALS
CHATR Stress MCHANICS OF MATRIALS and Strain Axial Loading Stress & Strain: Axial Loading Suitability of a structure or machine may depend on the deformations in the structure as well as the stresses induced
More informationAn Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation
An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation Nachiket Patil, Deepankar Pal and Brent E. Stucker Industrial Engineering, University
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 informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading
MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics
More informationREPRESENTING MATRIX CRACKS THROUGH DECOMPOSITION OF THE DEFORMATION GRADIENT TENSOR IN CONTINUUM DAMAGE MECHANICS METHODS
20 th International Conference on Composite Materials Copenhagen, 19-24 th July 2015 REPRESENTING MATRIX CRACKS THROUGH DECOMPOSITION OF THE DEFORMATION GRADIENT TENSOR IN CONTINUUM DAMAGE MECHANICS METHODS
More informationMulti-scale approach of the mechanical behavior of RC structures Application to nuclear plant containment buildings
Multi-scale approach of the mechanical behavior of RC structures Application to nuclear plant containment buildings Martin DAVID June 19, 2012 1 Industrial context EDF is responsible for numerous reinforced
More informationComputational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem
Computational Inelasticity FHLN05 Assignment 2017 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than October
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More information6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS
6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS Blondet et al. [25] carried out a cyclic test on an adobe wall to reproduce its seismic response and damage pattern under in-plane loads. The displacement
More informationCHAPTER 6 MECHANICAL PROPERTIES OF METALS PROBLEM SOLUTIONS
CHAPTER 6 MECHANICAL PROPERTIES OF METALS PROBLEM SOLUTIONS Concepts of Stress and Strain 6.1 Using mechanics of materials principles (i.e., equations of mechanical equilibrium applied to a free-body diagram),
More informationFEM techniques for nonlinear fluids
FEM techniques for nonlinear fluids From non-isothermal, pressure and shear dependent viscosity models to viscoelastic flow A. Ouazzi, H. Damanik, S. Turek Institute of Applied Mathematics, LS III, TU
More informationSome Aspects Of Dynamic Buckling of Plates Under In Plane Pulse Loading
Mechanics and Mechanical Engineering Vol. 12, No. 2 (2008) 135 146 c Technical University of Lodz Some Aspects Of Dynamic Buckling of Plates Under In Plane Pulse Loading Katarzyna Kowal Michalska, Rados
More informationPrediction of Elastic Constants on 3D Four-directional Braided
Prediction of Elastic Constants on 3D Four-directional Braided Composites Prediction of Elastic Constants on 3D Four-directional Braided Composites Liang Dao Zhou 1,2,* and Zhuo Zhuang 1 1 School of Aerospace,
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationEffects of TGO Roughness on Indentation Response of Thermal Barrier Coatings
Copyright 2010 Tech Science Press CMC, vol.17, no.1, pp.41-57, 2010 Effects of Roughness on Indentation Response of Thermal Barrier Coatings Taotao Hu 1 and Shengping Shen 1,2 Abstract: In this paper,
More informationUsing the Abaqus CDP Model in Impact Simulations
SAFIR project (http://safir2018.vtt.fi/) The Finnish Research Programme on Nuclear Power Plant Safety 2015-2018 Using the Abaqus CDP Model in Impact Simulations Alexis Fedoroff Technical reseach centre
More informationAdvanced model for soft soils. Modified Cam-Clay (MCC)
Advanced model for soft soils. Modified Cam-Clay (MCC) c ZACE Services Ltd August 2011 1 / 62 2 / 62 MCC: Yield surface F (σ,p c ) = q 2 + M 2 c r 2 (θ) p (p p c ) = 0 Compression meridian Θ = +π/6 -σ
More informationTheoretical Manual Theoretical background to the Strand7 finite element analysis system
Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationThe Generalized Interpolation Material Point Method
Compaction of a foam microstructure The Generalized Interpolation Material Point Method Tungsten Particle Impacting sandstone The Material Point Method (MPM) 1. Lagrangian material points carry all state
More informationMicro-meso draping modelling of non-crimp fabrics
Micro-meso draping modelling of non-crimp fabrics Oleksandr Vorobiov 1, Dr. Th. Bischoff 1, Dr. A. Tulke 1 1 FTA Forschungsgesellschaft für Textiltechnik mbh 1 Introduction Non-crimp fabrics (NCFs) are
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 informationModeling Interfacial Flux layer Phenomena in the Shell/Mold Gap Using CON1D
Modeling Interfacial Flux layer Phenomena in the Shell/Mold Gap Using CON1D Ya Meng Department of Materials Science &. Engineering University of Illinois at Urbana-Champaign September 25, 2 Introduction
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 informationStructural Analysis I Chapter 4 - Torsion TORSION
ORSION orsional stress results from the action of torsional or twisting moments acting about the longitudinal axis of a shaft. he effect of the application of a torsional moment, combined with appropriate
More informationA consistent dynamic finite element formulation for a pipe using Euler parameters
111 A consistent dynamic finite element formulation for a pipe using Euler parameters Ara Arabyan and Yaqun Jiang Department of Aerospace and Mechanical Engineering, University of Arizona, Tucson, AZ 85721,
More informationNonlinear bending analysis of laminated composite stiffened plates
Nonlinear bending analysis of laminated composite stiffened plates * S.N.Patel 1) 1) Dept. of Civi Engineering, BITS Pilani, Pilani Campus, Pilani-333031, (Raj), India 1) shuvendu@pilani.bits-pilani.ac.in
More informationTHREE-DIMENSIONAL SIMULATION OF THERMAL OXIDATION AND THE INFLUENCE OF STRESS
THREE-DIMENSIONAL SIMULATION OF THERMAL OXIDATION AND THE INFLUENCE OF STRESS Christian Hollauer, Hajdin Ceric, and Siegfried Selberherr Institute for Microelectronics, Technical University Vienna Gußhausstraße
More informationFinite Element simulation of vibrating plastic components
IEI Mechanical Engineer - Solid Mechanics ISRN: LIU-IEI-TEK-A 13/01572 SE Finite Element simulation of vibrating plastic components Master thesis February 9, 2013 Author: Jesper Kihlander Examiner: Larsgunnar
More informationComputational Analysis for Composites
Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics
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 informationCreep and aging in nonlinear models for concrete
Creep and aging in nonlinear models for concrete Andrzej Truty Thomas Zimmermann Krzysztof Podleś Jan Szarliński Cracow University of Technology and ZACE Services ZACE Services prof.em. Cracow University
More informationA Repeated Dynamic Impact Analysis for 7x7 Spacer Grids by using ABAQUS/ Standard and Explicit
A Repeated Dynamic Impact Analysis for 7x7 Spacer Grids by using ABAQUS/ Standard and Explicit Kim, Jae-Yong, and Yoon, Kyung-Ho* * Korea Atomic Energy Research Institute ABSTRACT Spacer grids(sg) are
More informationC.J. Bennett, W. Sun Department of Mechanical, Materials and Manufacturing Engineering, University of Nottingham, Nottingham NG7 2RD, UK
Optimisation of material properties for the modelling of large deformation manufacturing processes using a finite element model of the Gleeble compression test C.J. Bennett, W. Sun Department of Mechanical,
More informationStructural behaviour of traditional mortise-and-tenon timber joints
Structural behaviour of traditional mortise-and-tenon timber joints Artur O. Feio 1, Paulo B. Lourenço 2 and José S. Machado 3 1 CCR Construtora S.A., Portugal University Lusíada, Portugal 2 University
More informationThe Finite Element Method for Mechonics of Solids with ANSYS Applicotions
The Finite Element Method for Mechonics of Solids with ANSYS Applicotions ELLIS H. DILL 0~~F~~~~"P Boca Raton London New Vork CRC Press is an imprint 01 the Taylor & Francis Group, an Informa business
More informationFrequently Asked Questions
Frequently Asked Questions Why do we have to make the assumption that plane sections plane? How about bars with non-axis symmetric cross section? The formulae derived look very similar to beam and axial
More informationSettlement and Bearing Capacity of a Strip Footing. Nonlinear Analyses
Settlement and Bearing Capacity of a Strip Footing Nonlinear Analyses Outline 1 Description 2 Nonlinear Drained Analysis 2.1 Overview 2.2 Properties 2.3 Loads 2.4 Analysis Commands 2.5 Results 3 Nonlinear
More informationFinal Project: Indentation Simulation Mohak Patel ENGN-2340 Fall 13
Final Project: Indentation Simulation Mohak Patel ENGN-2340 Fall 13 Aim The project requires a simulation of rigid spherical indenter indenting into a flat block of viscoelastic material. The results from
More informationA FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE
A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE Daxu Zhang and D. R. Hayhurst School of Mechanical, Aerospace and Civil Engineering,
More informationLecture 4 Implementing material models: using usermat.f. Implementing User-Programmable Features (UPFs) in ANSYS ANSYS, Inc.
Lecture 4 Implementing material models: using usermat.f Implementing User-Programmable Features (UPFs) in ANSYS 1 Lecture overview What is usermat.f used for? Stress, strain and material Jacobian matrix
More informationAnalysis of a Casted Control Surface using Bi-Linear Kinematic Hardening
Analysis of a Casted Control Surface using Bi-Linear Kinematic Hardening Abdul Manan Haroon A. Baluch AERO, P.O Box 91, Wah Cantt. 47040 Pakistan Abstract Control Surfaces or Fins are very essential parts
More informationInverse method for flow stress parameters identification of tube bulge hydroforming considering anisotropy
Int. J. Mechatronics and Manufacturing Systems, Vol. 4, No. 5, 2011 441 Inverse method for flow stress parameters identification of tube bulge hydroforming considering anisotropy T. Zribi, A. Khalfallah*
More informationLimit analysis of brick masonry shear walls with openings under later loads by rigid block modeling
Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling F. Portioli, L. Cascini, R. Landolfo University of Naples Federico II, Italy P. Foraboschi IUAV University,
More informationFinite element simulation of residual stresses in laser heating
IAS-2008-66-546ST Finite element simulation of residual stresses in laser heating G. H. Farrahi 1, M. Sistaninia 2, H. Moeinoddini 3 1,2-School of Mechanical Engineering, Sharif University of Technology,
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material
ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a
More informationFEA A Guide to Good Practice. What to expect when you re expecting FEA A guide to good practice
FEA A Guide to Good Practice What to expect when you re expecting FEA A guide to good practice 1. Background Finite Element Analysis (FEA) has transformed design procedures for engineers. Allowing more
More informationCHAPTER 7 FINITE ELEMENT ANALYSIS OF DEEP GROOVE BALL BEARING
113 CHAPTER 7 FINITE ELEMENT ANALYSIS OF DEEP GROOVE BALL BEARING 7. 1 INTRODUCTION Finite element computational methodology for rolling contact analysis of the bearing was proposed and it has several
More informationCode_Aster. HSNV129 - Test of compression-dilation for study of the coupling thermics-cracking
Titre : HSNV129 - Essai de compression-dilatation pour étu[...] Date : 28/2/218 Page : 1/8 HSNV129 - Test of compression-dilation for study of the coupling thermics-cracking Summary: One applies to an
More informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
More informationHSNV140 - Thermoplasticity with restoration of work hardening: test of blocked dilatometry
Titre : HSNV140 - Thermo-plasticité avec restauration d'éc[...] Date : 29/06/2015 Page : 1/11 HSNV140 - Thermoplasticity with restoration of work hardening: test of blocked dilatometry Summary: This test
More informationChapter. Materials. 1.1 Notations Used in This Chapter
Chapter 1 Materials 1.1 Notations Used in This Chapter A Area of concrete cross-section C s Constant depending on the type of curing C t Creep coefficient (C t = ε sp /ε i ) C u Ultimate creep coefficient
More informationCellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996).
1 Cellular solid structures with unbounded thermal expansion Roderic Lakes Journal of Materials Science Letters, 15, 475-477 (1996). Abstract Material microstructures are presented which can exhibit coefficients
More informationMODELING GEOMATERIALS ACROSS SCALES
MODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING AFOSR WORKSHOP ON PARTICULATE MECHANICS JANUARY 2008 COLLABORATORS: DR XUXIN TU AND MR KIRK ELLISON
More informationINTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS
INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS SHEN R. WU and LEI GU WILEY A JOHN WILEY & SONS, INC., PUBLICATION ! PREFACE xv PARTI FUNDAMENTALS 1 1 INTRODUCTION 3
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 informationDYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION
October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,
More information