SPH development at Cranfield University. Prof. Rade Vignjevic. IV SPHERIC Nantes, May 2009
|
|
- Dana Manning
- 5 years ago
- Views:
Transcription
1 SPH development at Cranfield University Prof. Rade Vignevic IV SPHERIC Nantes, May 2009
2 SPH development at Cranfield University Presentation outline 1. Introduction (CU, Motivation) 2. Normalised Corrected SPH 3. Non-collocational SPH 4. Contact algorithm, FE SPH coupling 5. Damage modelling
3 Cranfield University: History Formed in 1946 as College of Aeronautics (centre of excellence in aerospace) Postgraduate University 1992 Centre of Applied Research Five Maor Schools Health Management Engineering Applied Science
4 Cranfield University: Entirely Post Graduate, 35% PhD MSc Courses of 45 Week Duration International Industrial & Government Funded Research The UK National Flying Laboratory Applied approach to teaching and research Strong Industrial Links Worldwide Reputation 240 Million Turnover - 60% From Research
5 Motivation for our SPH work Challenging problems of computational mechanics are often characterised by: Extremely large deformations Tracking of interfaces between solids and liquid/gas Propagation of cracks with arbitrary and complex paths Change of boundaries of material that has undergone extensive micro cracking or phase change and failure.
6 Motivation for the SPH work Problem solving Manufacturing: Extrusion, Moulding, Friction Stir Welding, Manufacture of composites, Safety and Crashworthiness: Hypervelocity impacts on spacecraft, Aircraft impacts on water/soft soil, Bird strike... Defence: Armour penetration (metal, ceramics), Warhead fragmentation, Shape charges
7 People that inspired and influenced SPH research at Cranfield: L. Libersky J. J. Monaghan
8 Corrected Normalised SPH (CNSPH) R. Vignevic, J. Campbell, L. Libersky; A treatment of zero-energy modes in the smoothed particle hydrodynamics method, Comput. Methods Appl. Mech. Engrg. 184 (2000), pp , (Received 28 October 1998) R. Vignevic, J. Reveles, J. Campbell; SPH in a Total Lagrangian Formalism, Computer Methods in Engineering and Science, vol.14, no.3, pp , 2006
9 Corrected Normalised SPH CNSPH Emmy Noether, "Invariante Variationsprobleme, " Nachr. v. d. Ges. d. Wiss. zu Göttingen 1918, pp Specifically: The invariance with respect to translational and rotational transformations
10 Corrected Normalised SPH CNSPH Space homogeneity x m x x x x (, h) J = JW I J J J is interpolated solution space x m x x x x ' = x x ( h) J ' ' ' ' ' =, ' JW I x = x J I J J mj J x ' = x x W ( xi x J, h ) ( xi xj h ) J J J m W J, = 1
11 Corrected Normalised SPH CNSPH Space isotropy m x x x x (, h) J = JW I J J J x is interpolated solution space x m x x x ( h) J ' ' ' ' ' =, ' JW I x = x J I J J x' = C x x' C x = C x = C x C = C nnbr mj xj W ( xi xj, h) = I J = 1 J
12 Interpolation should not violate the following properties of space: Homogeneity Anisotropy Condit. 1 ), ( 1 = = h x x W m nnbr i x = 1 = ), ( 1 h x x W m nnbr i Correct. ), ( ), ( ~ 1 h x x W m h x x W W nnbr i i i = = 1 1 ~ ~ ~ ~ = = nnbr i i i x m W W W Corrected Normalised SPH CNSPH
13 CNSPH Form of Governing Equations Conservation of: m = v v % W % nnbr & Mass ( ) i i i i = 1 n n b r i = + i Momentum 2 2 v& m σ σ = 1 i nnbr σ = v v % W % i e& m Energy ( ) 2 i = 1 i i % W % i W here: ~ W i nnbr = Wi = 1 m 1 ( ) ~ ~ ~ xi x W nnbr i, W i= W i ( ) xi x = 1 m ~ W i 1
14 CNSPH Conservation of Momentum and A-momentum Materials X -Velocity
15 Corrected Normalised SPH Summary First order consistency Conservation of linear and angular momentum Homogeneity and isotropy of space maintained in the SPH discretisation
16 SPH in a Total Lagrangian Formalism R. Vignevic, J. Reveles, J. Campbell; SPH in a Total Lagrangian Formalism, Computer Methods in Engineering and Science, Vol.14, No.3, pp , 2006
17 SPH in a Total Lagrangian Formalism Taylor test for OFHPC copper 180 m/s
18 SPH in a Total Lagrangian Formalism The mapping from material into spatial coordinates is x ( X,t) = φ. The displacement of a material point is given by the difference between its current position and its original position. u ( X, t) = φ ( X, t) φ( X,0) = φ( X, t) X = x X
19 SPH in a Total Lagrangian Formalism The deformation gradient F ( u + X) x u F = = = + I X X X Discretised form F ~ 0 I ( u u ) W V. B I B ( X Xi ) = J S Velocity gradient J I 0 IJ J + ( ) F& = v v W% V 0 I J I 0 IJ J J S m = S ~ 1 W
20 SPH in a Total Lagrangian Formalism Conservation Equations in the Total Lagrangian formalism Mass Momentum Energy Continuous Discretised = J 1 1 u& = 0 P + b 0 0 = F I & 0 a ( P P ) V B 1 & = F & 0 : P = I J I 0W IJ J : J S e m J ~ 0 ( ) e& I ~ = P J : v J v I 0WIJ VJ. B J S I J
21 SPH in a Total Lagrangian Formalism Taylor test for OFHPC copper 180 m/s
22 SPH in a Total Lagrangian Formalism 3-D Taylor test, effective plastic strain
23 SPH in a Total Lagrangian Formalism Total Lagrangian SPH vs. conventional SPH
24 SPH in a Total Lagrangian Formalism Summary Stable well behaved Applicable to finite deformations Combined with Eulerian SPH when modelling extremely large deformations and failure
25 Non-Collocational SPH R. Vignevic, J. Campbell, L. Libersky; A treatment of zero-energy modes in the smoothed particle hydrodynamics method, Comput. Methods Appl. Mech. Engrg. 184 (2000), pp , (Received 28 October 1998)
26 Non-Collocational SPH Two types of particles: Velocity Stress Easy to apply essential boundary conditions (similar to FE)
27 Non-Collocational SPH Collision of elastic hoops
28 Non-Collocational SPH Taylor test for OFHPC copper 180 m/s
29 Non-Collocational SPH Summary Easy to apply to 1D and 2D structural elements, especially advantageous when modelling failure (ongoing work) Increased difficulties to extend to 3D continuum (update of stress particle locations)
30 Contact Algorithm R. Vignevic; T. De Vuyst; and J. Campbell; A Frictionless Contact Algorithm for Meshless Methods, CMES, Vol. 13, No. 1, pp , 2006 T. De Vuyst, R. Vignevic and J. C. Campbell; Modelling of Fluid-Structure Impact Problems using a Coupled SPH-FE solver, Journal of Impact Engineering, Vol. 31, No. 8, pp , 2005 J. Campbell, R. Vignevic, L. Libersky; A Contact Algorithm for Smoothed Particle Hydrodynamics, Computer Methods in Applied Mechanics and Engineering, Vol. 184, No. 1, pp , 2000
31 Contact Initial Boundary Value Problem Traction on the contact surface Ω A Γ CA tn = nˆ σ Γ CB Γ A Γ CB Γ B The Hertz Signorini Moreau (Kuhn-Tucker) Conditions for frictionless contact: g 0, t 0, g = 0 n t n ( ) = nˆ x x g Γ CB B A Ω B In the variational rate form the contact constraint imposed by: δg = δg n g & ( t ) = 0 For a body in contact g& = 0
32 Contact Initial Boundary Value Problem t σ + b = a on Ω = Ω A U Ω B σ n = t on Γ t, Γ t Ω A Γ A u t Γ t = u on = on = Γ Γ A t B C N g t N 0 0 g A = 0 B Γ u Γ to Γ C Γ C Γ B Γ u Ω B Γ t t Where Γ t Γ to Γ u = Γ, and Γ t Γ to =, Γ to Γ u =, Γ t Γ u =.
33 Discretised Contact Initial Boundary Value Problem t wσ n dγ w σ dv = w( d&& b) dv wtdγ Γ Ω Ω Γ A weak form of the initial boundary value problem with contact. This equation discretised in space: T T T T N N d B σ d N b N b c Ω Ω Ω Ω Γ Γ t u dv&& + dv dv dv ( T N σ) n d 0 Γ = Where: w = Nd and N is a shape functions matrix CBL c
34 Discretised Contact Initial Boundary Value Problem In the SPH method N has the following form: N i = m W ( x ) np = 1 W i ( x ) i Where: I - particle at which shape function is evaluated, J - particle at which the shape function is centered, Np - number of neighbors for the i particle, d - nodal displacement vector, W - SPH kernel function (B-spline)
35 Contact Force Due to the diffused nature of boundary in the conventional SPH the contact force is defined as: f = T c b c Ω CBL N dv Where specific contact force is defined as gradient of contact potential: ( x x ) ( ) W A B φ x = K dv φ( x A i ) W ( p ) Ω CBL avg n = NCONT m W K W ( ri ) ( p ) avg n. m ( ) n 1 i ( p ) W r NCONT bc ( xi ) = φ( xi ) = Kn n xw x i i W avg ( x )
36 Normal Contact - SPH Ω A Γ A Γ CA Γ CB Γ B Ω B Γ CA Γ CB 2h
37 Example, Plate Impact (2D) The symmetrical block impact Each block was discretised with 50 by 20 particles Steel blocs modelled with an elastic-plastic material model
38 Example, Plate Impact (2D)
39 Example, Hypervelocity Impact (3D) 7km/s impact of aluminium proectile on aluminium bumper plate
40 Contact Algorithm Summary Well suited for meshless methods Numerically effective, approximately additional 10% CPU time in 3D simulations Easy to add friction (ongoing work)
41 Alternative Form of the SPH Equations R. Vignevic, J. Campbell, J. Jaric, S. Powell; Derivation of SPH equations in a moving referential coordinate system, Comput. Methods Appl. Mech. Engrg. 198 (2009), pp , (Received 8 September 2008)
42 Moving Referential Configuration x 2 x 2 J Particle I Particle Support Moving Control of particle I Volume x x xr x xr x x R J Particle x 1 J Particle Fixed Frame frame fixed on particle centred on particle Ion I particle I Df % Df r = + vr { Dt Dt referential f r r r r& v v v R = J I h x 1
43 Moving Referential Configuration Conservation of Mass D% r r D% r r r + v + vr = + v + ( v vcv ) = Dt Dt D% + ( v) vcv = 0 Dt D% r r W IJ d Ω + ( v) W IJ d Ω + vcv W IJ = 0 Dt Ω Ω Ω Integrating by parts and neglecting the boundary terms & v W ( x ' x, h) d v W ( x ' x, h) d = Ω I Ω Ω Ω & = m ( v v ) J W ( ) J J I IJ
44 Moving Referential Configuration Conservation of Momentum r r ( v) r r D% v r r + ( v v R ) = + R = σ t v v Dt D % r v r r 1 W d W d W d IJ Ω + v vr IJ Ω = σ IJ Ω Dt Ω Ω Ω Using the divergence theorem, dropping the boundary terms and linearizing the velocity in the second term: 1 v r & + v r v r r r R WIJ dω v ( vr WIJ ) dω = σwij dω Ω Ω Ω
45 Moving Referential Configuration Conservation of Momentum 1 v r & + v r v r r r R WIJ dω v ( vr WIJ ) dω = σwij dω Ω Ω Ω When discretised, the above equation becomes: r m 1 m v& r r r v v v ( ) J ( ) J W = ( σ + σ ) I J I R IJ I J x IJ J J I J Ω J W
46 Comparison Conv. SPH [1, 2] Continuity equation v v mj & I = I ( I J ) xwij J Ω J r r & = mj vj vi WIJ Moving C. S. ( ) Conv. SPH J Momentum equation 1 mj v & = σ + σ W I J Ω J r mj 1 v& r r r v v v W [1, 2] ( ) Moving C. S. I I J x IJ ( ) ( ) ( σ σ ) = + J I J I R IJ I J x IJ J J I J Ω J 1.- Larry D. Libersky at. al. High Strain Lagrangian Hydrodynamics a 3-DSPH code for dynamic material response,. Journal of Computational Physics, Morris, J.P., An Overview of the Method of Smoothed Particle Hydrodynamics, November 1995, Universitat Kaiserlautern, Arbeitsgruppe Technomathematik. m W
47 Numerical Example Shock tube problem 0<x< <x<1 Initial conditions 0<x< <x<1 Left Right Density Pressure Velocity
48 Numerical Example Shock tube problem 0.12 Velocity at time Velocity Exact Basic SPH, h constant Moving SPH Position
49 Numerical Example Shock tube problem 1.05 Pressure at time Pressure Exact Basic SPH, h constant Moving SPH Position
50 Numerical Example Shock tube problem 1.1 Density at time Exact Conventional SPH Moving SPH Density Position
51 Alternative Form of the SPH Equations Summary Part of the ongoing work Framework for a more rigorous treatment of variable h interpolation In the shock tube problem performance similar to the conventional SPH with constant h
52 Damage Modelling The fracture model capable of simulating: Initiation and growth of damage Crack formation and propagation in an arbitrary direction Crack branching and crack oining, leading to fragmentation
53 How to represent a Crack Cracked/failed particle concept: A failed particle is split into two particles, this approach has been demonstrated by Rabczuk and Belytschko (2007)
54 Continuum Damage Mechanics Compatible with the concept of effective stress introduced by Kachanov in 1958 The effective stress tensor is defined as: ~ σ = 1 D S S D = S σ 0 0 The concept of the interaction area relatively simple to apply within SPH
55 How to represent a Crack (Damage) Inter-particle failure Damage is evaluated for pairs of neighbour particles. Following fracture the particles cease being neighbours. The concept of particle-particle interaction area (Swegle, 2000), - damage affects interaction effective area - at failure the effective area is at a critical value
56 The force acting on a surface due to a stress is given by The SPH momentum equation could be rewritten in term of an interaction area: Swegle Interaction Area i i i i i i W m dt dv a + = = 2 2 σ σ ( ) ( ) ( ) ( ) + = = + = i i i i i i i i i i i i i A F W vol vol m a σ σ σ σ = F σ A
57 Inter-particle failure i Damage is evolved as an inter-particle value, D i,which reduces the inter-particle interaction area: F i i = ( σ ) + ( σ ) A ( 1 D ) i i i i When damage reaches a critical value the material is assumed to have failed and the particles cease to be neighbours.
58 Spall demonstration problem Normal stress between pairs of particles is compared to a spall criterion. Damage initiated once this criterion is exceeded. V = 450 m/s PMMA Copper plate 5mm thick Copper flyer 2.5mm thick 1D strain state Spall plane opens in middle of target plate
59 Representing fracture within a meshless method
60 Numerical Results Stress-xx (Gpa) Copper Plate Impact Test (with contact / artificial viscosity Q=2.0, L=0.2) Information displayed for Particle 195 (in PMMA) No damage D1 D Time (us)
61 The Mock-Holt problem Mock and Holt (1983), experimental results from a number of tests on explosive driven Fragmentation of metallic cylinders.
62 The Mock-Holt problem
63 Damage Modelling Summary Well suited for meshless methods Numerically effective Can deal with crack closure Ongoing work
64 Overall Summary All developments discussed were implemented into our 3D SPH code The SPH code coupled with DYNA3D The codes are routinely used to analyse real world problems Further development of the SPH method is of high priority to us (stabilisation of the Eulerian SPH, development of SPH structural elements, modelling damage, )
65 The End
Modelling of bird strike on the engine fan blades using FE-SPH
Modelling of bird strike on the engine fan blades using FE-SPH Dr Nenad Djordjevic* Prof Rade Vignjevic Dr Tom De Vuyst Dr James Campbell Dr Kevin Hughes *nenad.djordjevic@brunel.ac.uk MAFELAP 2016, 17
More informationBIRD-STRIKE IMPACT SIMULATION WITH AN AIRCRAFT WING USING SPH BIRD MODEL
BIRD-STRIKE IMPACT SIMULATION WITH AN AIRCRAFT WING USING SPH BIRD MODEL BOGDAN ALEXANDRU BELEGA 1 Abstract: In this paper I focus on developing a model to simulate a birdstrike with an aircraft wing using
More informationSPH for simulating impacts and collisions
SPH for simulating impacts and collisions Thomas I. Maindl University of Vienna Jan 22, 2013 Christoph Schäfer, Roland Speith (Eberhard Karls University of Tübingen) Agenda SPH quick overview Solid body
More informationSMOOTHED PARTICLE HYDRODYNAMICS METHOD IN MODELING OF STRUCTURAL ELEMENTS UNDER HIGH DYNAMIC LOADS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China SMOOTHE PARTICLE HYROYAMICS METHO I MOELIG OF STRUCTURAL ELEMETS UER HIGH YAMIC LOAS. Asprone *, F. Auricchio, A. Reali,
More informationIMPLEMENTATION OF THE SMOOTHED PARTICLE HYDRODYNAMICS METHOD TO SOLVE PLASTIC DEFORMATION IN METALS
Blucher Mechanical Engineering Proceedings May 014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm IMPLEMETATIO OF THE SMOOTHED PARTICLE HYDRODYAMICS METHOD TO SOLVE PLASTIC DEFORMATIO I METALS
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
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 informationElsevier Editorial System(tm) for Computer Methods in Applied Mechanics and Engineering Manuscript Draft
Elsevier Editorial System(tm) for Computer Methods in Applied Mechanics and Engineering Manuscript Draft Manuscript Number: Title: Modelling of Strain Softening Materials Based on Equivalent Damage Force
More informationThe Applications of Meshless Local Petrov-Galerkin (MLPG) Approaches in High-Speed Impact, Penetration and Perforation Problems
Copyright c 006 Tech Science Press CMES, vol.14, no., pp.119-18, 006 The Applications of Meshless Local Petrov-Galerkin (MLPG) Approaches in High-Speed Impact, Penetration and Perforation Problems. D.
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationSimulation of Particulate Solids Processing Using Discrete Element Method Oleh Baran
Simulation of Particulate Solids Processing Using Discrete Element Method Oleh Baran Outline DEM overview DEM capabilities in STAR-CCM+ Particle types and injectors Contact physics Coupling to fluid flow
More informationWARHEAD FRAGMENTATION MODELING WITH PERIDYNAMICS
23 RD INTERNATIONAL SYMPOSIUM ON BALLISTICS TARRAGONA, SPAIN 16-20 APRIL 2007 WARHEAD FRAGMENTATION MODELING WITH PERIDYNAMICS Paul N. Demmie 1, Dale S. Preece 1 and Stewart A. Silling 1 1 Sandia National
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
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 informationSIMULATIONS OF HYPERVELOCITY IMPACTS WITH SMOOTHED PARTICLE HYDRODYNAMICS
SIMULATIONS OF HYPERVELOCITY IMPACTS WITH SMOOTHED PARTICLE HYDRODYNAMICS Dominique Lacerda, Jean-Luc Lacome, DYNALIS, Paris, France Email : dynalis@dynalis.fr ABSTRACT This paper is devoted to the results
More informationIntroduction to fracture mechanics
Introduction to fracture mechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 6-9 November, 2017 Institute of Structural Engineering, ETH Zu rich November 9, 2017 Institute
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationA METHOD TO ASSESS IMPACT DAMAGE USING A SMOOTHED PARTICLE HYDRODYNAMICS AND FINITE ELEMENT COUPLED APPROACH
Transactions, SMiRT-23, Paper ID 315 A METHOD TO ASSESS IMPACT DAMAGE USING A SMOOTHED PARTICLE HYDRODYNAMICS AND FINITE ELEMENT COUPLED APPROACH 1 Consultant, Amec Foster Wheeler, UK 2 Principal Consultant,
More informationJournal of Basic and Applied Sciences. Australian Journal of Basic and Applied Sciences, 9(37) Special 2015, Pages: ISSN:
ISSN:1991-8178 Australian Journal of Basic and Applied Sciences Journal home page: www.ajbasweb.com The Development of Unique Orthogonal Rotation Tensor Algorithm in the LLNL- DYNA3D for Orthotropic Materials
More informationNUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE
NUMERICAL SIMULATION OF FLUID-STRUCTURE INTERACTION PROBLEMS WITH DYNAMIC FRACTURE Kevin G. Wang (1), Patrick Lea (2), and Charbel Farhat (3) (1) Department of Aerospace, California Institute of Technology
More informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More informationCoupled 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 informationarxiv: v1 [cs.ce] 24 Jan 2019
A Total Lagrangian SPH Method for Modelling Damage and Failure in Solids Md Rushdie Ibne Islam a, Chong Peng a,b, a ESS Engineering Software Steyr GmbH, Berggasse 35, 4400 Steyr, Austria b Institut für
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationDEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD R. A. Amaro
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigid-body spring deformation expressed in terms of field variables assembly of rigid-bodies connected by spring Distinct Element Method (DEM) simple
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 informationSmoothed Particle Hydrodynamics for Nonlinear Solid Mechanics
Noname manuscript No. (will be inserted by the editor) Smoothed Particle Hydrodynamics for Nonlinear Solid Mechanics Ernesto B. Ismail B. Daya Reddy Received: date / Accepted: date Abstract SPH is used
More informationModelling Dynamic Behaviour and Spall Failure of Aluminium Alloy AA7010
Journal of Physics: Conference Series PAPER OPEN ACCESS Modelling Dynamic Behaviour and Spall Failure of Aluminium Alloy AA7010 To cite this article: N Ma at et al 2017 J. Phys.: Conf. Ser. 914 012033
More informationTreecodes for Cosmology Thomas Quinn University of Washington N-Body Shop
Treecodes for Cosmology Thomas Quinn University of Washington N-Body Shop Outline Motivation Multipole Expansions Tree Algorithms Periodic Boundaries Time integration Gravitational Softening SPH Parallel
More informationHigh strain rate fracture behaviour of fused silica
Journal of Physics: Conference Series OPEN ACCESS High strain rate fracture behaviour of fused silica To cite this article: A Ruggiero et al 2014 J. Phys.: Conf. Ser. 500 182036 View the article online
More informationArchetype-Blending Multiscale Continuum Method
Archetype-Blending Multiscale Continuum Method John A. Moore Professor Wing Kam Liu Northwestern University Mechanical Engineering 3/27/2014 1 1 Outline Background and Motivation Archetype-Blending Continuum
More informationNatural States and Symmetry Properties of. Two-Dimensional Ciarlet-Mooney-Rivlin. Nonlinear Constitutive Models
Natural States and Symmetry Properties of Two-Dimensional Ciarlet-Mooney-Rivlin Nonlinear Constitutive Models Alexei Cheviakov, Department of Mathematics and Statistics, Univ. Saskatchewan, Canada Jean-François
More informationOn angular momentum balance for particle systems with periodic boundary conditions
On angular momentum balance for particle systems with periodic boundary conditions Vitaly A. Kuzkin 1 Institute for Problems in Mechanical Engineering RAS Saint Petersburg State Polytechnical University
More information1.050 Engineering Mechanics. Lecture 22: Isotropic elasticity
1.050 Engineering Mechanics Lecture 22: Isotropic elasticity 1.050 Content overview I. Dimensional analysis 1. On monsters, mice and mushrooms 2. Similarity relations: Important engineering tools II. Stresses
More informationParticle-based Fluids
Particle-based Fluids Particle Fluids Spatial Discretization Fluid is discretized using particles 3 Particles = Molecules? Particle approaches: Molecular Dynamics: relates each particle to one molecule
More informationSmoothed Particle Hydrodynamics (SPH) 4. May 2012
Smoothed Particle Hydrodynamics (SPH) 4. May 2012 Calculating density SPH density estimator Weighted summation over nearby particles: ρ(r) = N neigh b=1 m bw (r r b, h) W weight function with dimension
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 informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationAn Atomistic-based Cohesive Zone Model for Quasi-continua
An Atomistic-based Cohesive Zone Model for Quasi-continua By Xiaowei Zeng and Shaofan Li Department of Civil and Environmental Engineering, University of California, Berkeley, CA94720, USA Extended Abstract
More informationMaterial stability analysis of particle methods
Advances in Computational Mathematics (2005) 2: 171 190 Springer 2005 Material stability analysis of particle methods S.P. Xiao a and T. Belytschko b, a Department of Mechanical and Industrial Engineering
More informationApplication of a non-local failure criterion to a crack in heterogeneous media S. Bavaglia*, S.E. Mikhailov*
Application of a non-local failure criterion to a crack in heterogeneous media S. Bavaglia*, S.E. Mikhailov* University of Perugia, Italy Email: mic@unipg.it ^Wessex Institute of Technology, Ashurst Lodge,
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More informationNotes. Multi-Dimensional Plasticity. Yielding. Multi-Dimensional Yield criteria. (so rest state includes plastic strain): #=#(!
Notes Multi-Dimensional Plasticity! I am back, but still catching up! Assignment is due today (or next time I!m in the dept following today)! Final project proposals: I haven!t sorted through my email,
More informationPREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures
PREDICTIVE SIMULATION OF UNDERWATER IMPLOSION: Coupling Multi-Material Compressible Fluids with Cracking Structures Kevin G. Wang Virginia Tech Patrick Lea, Alex Main, Charbel Farhat Stanford University
More informationModelling hypervelocity impact in DYNA3D J. Campbell, R. Vignjevic College of Aeronautics, Cranfield University, Cranfield,
Modelling hypervelocity impact in DYNA3D J. Campbell, R. Vignjevic College of Aeronautics, Cranfield University, Cranfield, Abstract This paper presents part of the work on development of tools and modelling
More informationMicro-mechanics in Geotechnical Engineering
Micro-mechanics in Geotechnical Engineering Chung R. Song Department of Civil Engineering The University of Mississippi University, MS 38677 Fundamental Concepts Macro-behavior of a material is the average
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationStress, Strain, and Viscosity. San Andreas Fault Palmdale
Stress, Strain, and Viscosity San Andreas Fault Palmdale Solids and Liquids Solid Behavior: Liquid Behavior: - elastic - fluid - rebound - no rebound - retain original shape - shape changes - small deformations
More informationContents. I Introduction 1. Preface. xiii
Contents Preface xiii I Introduction 1 1 Continuous matter 3 1.1 Molecules................................ 4 1.2 The continuum approximation.................... 6 1.3 Newtonian mechanics.........................
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 information2016, Dr. Phil Vardon, All rights reserved. DARSim thanks our invited speakers for providing their presentation materials.
Copyright note: These slides are made available for educational purposes only. No reproduction neither usage of these slides in full or part are allowed, unless a written approval of Dr. Phil Vardon is
More informationModeling and Simulations of Aircraft Structures Stiffness, Damage, and Failure Prediction for Laminated Composites
Modeling and Simulations of Aircraft Structures Stiffness, Damage, and Failure Prediction for Laminated Composites H.E.Pettermann 1, C.Schuecker 1, D.H.Pahr 2, F.G.Rammerstorfer 2 1 Austrian Aeronautics
More informationBifurcation Analysis in Geomechanics
Bifurcation Analysis in Geomechanics I. VARDOULAKIS Department of Engineering Science National Technical University of Athens Greece and J. SULEM Centre d'enseignement et de Recherche en Mecanique des
More informationConservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body
Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference
More informationFinite Element Model of a complex Glass Forming Process as a Tool for Control Optimization
Excerpt from the Proceedings of the COMSOL Conference 29 Milan Finite Element Model of a complex Glass Forming Process as a Tool for Control Optimization Felix Sawo and Thomas Bernard Fraunhofer Institute
More informationSimulation of Impact Proof Testing of Electronic Sub- Systems
12 th International LS-DYNA Users Conference Blast/Impact(3) Simulation of Impact Proof Testing of Electronic Sub- Systems Paul Glance, PhD ME Naval Air Warfare Center Weapons Division China Lake, California
More informationSimulation and Test Validation of Windscreen Subject to Pedestrian Head Impact
12 th International LS-DYNA Users Conference Occupant Safety Simulation and Test Validation of Windscreen Subject to Pedestrian Head Impact Qi Liu, Junyong Liu, Qiang Miao, Dazhi Wang, Xiaodong Tang SAIC
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 informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
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 information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
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 informationPSEUDO ELASTIC ANALYSIS OF MATERIAL NON-LINEAR PROBLEMS USING ELEMENT FREE GALERKIN METHOD
Journal of the Chinese Institute of Engineers, Vol. 27, No. 4, pp. 505-516 (2004) 505 PSEUDO ELASTIC ANALYSIS OF MATERIAL NON-LINEAR PROBLEMS USING ELEMENT FREE GALERKIN METHOD Raju Sethuraman* and Cherku
More informationMultiscale modeling and simulation of shock wave propagation
University of Iowa Iowa Research Online Theses and Dissertations Spring 2013 Multiscale modeling and simulation of shock wave propagation Cory Nelsen University of Iowa Copyright 2013 Cory Nelsen This
More informationModeling Hailstone Impact onto Composite Material Panel Under a Multi-axial State of Stress
Modeling Hailstone Impact onto Composite Material Panel Under a Multi-axial State of Stress Authors Marco ANGHILERI * Luigi-M L CASTELLETTI * Andrea MILANESE * and Andrea SEMBOLONI * Affiliation * Politecnico
More informationCohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition framework
University of Liège Aerospace & Mechanical Engineering MS3: Abstract 131573 - CFRAC2017 Cohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition
More informationHigh-resolution finite volume methods for hyperbolic PDEs on manifolds
High-resolution finite volume methods for hyperbolic PDEs on manifolds Randall J. LeVeque Department of Applied Mathematics University of Washington Supported in part by NSF, DOE Overview High-resolution
More informationOutline. Advances in STAR-CCM+ DEM models for simulating deformation, breakage, and flow of solids
Advances in STAR-CCM+ DEM models for simulating deformation, breakage, and flow of solids Oleh Baran Outline Overview of DEM in STAR-CCM+ Recent DEM capabilities Parallel Bonds in STAR-CCM+ Constant Rate
More informationChapter 1 Introduction
Chapter 1 Introduction This thesis is concerned with the behaviour of polymers in flow. Both polymers in solutions and polymer melts will be discussed. The field of research that studies the flow behaviour
More informationContinuous Gradient Method to Cell Crossing Instability
Continuous Gradient Method to Cell Crossing Instability D. Z. Zhang X. Ma P. T. Giguere Fluid Dynamics and Solid Mechanics Group, T-3, Theoretical Division Los Alamos National Laboratory Los Alamos, New
More informationConfigurational Forces as Basic Concepts of Continuum Physics
Morton E. Gurtin Configurational Forces as Basic Concepts of Continuum Physics Springer Contents 1. Introduction 1 a. Background 1 b. Variational definition of configurational forces 2 с Interfacial energy.
More information1.050 Content overview Engineering Mechanics I Content overview. Selection of boundary conditions: Euler buckling.
.050 Content overview.050 Engineering Mechanics I Lecture 34 How things fail and how to avoid it Additional notes energy approach I. Dimensional analysis. On monsters, mice and mushrooms Lectures -3. Similarity
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationFigure 6. Partom diagram for two layer target.
Lecture 11 Penetration Mechanics II Figure 6. Partom diagram for two layer target. 3. Layered Targets For relatively long projectiles, it is often very helpful to analyze impacts in terms of a P-L diagram
More informationLiu Fu 1,2, Zhang Jiazhen 1, Hu Zhongmin 1, and Zhang Mingyi 1 1
Liu Fu 1,2, Zhang Jiazhen 1, Hu Zhongmin 1, and Zhang Mingyi 1 1 Beijing Aeronautical Science & Technology Research Institute of COMAC, Beijing 102211, PR China 2 Nanjing University of Aeronautics and
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 informationDevelopment and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks
Development and application of time-lapse ultrasonic tomography for laboratory characterisation of localised deformation in hard soils / soft rocks Erika Tudisco Research Group: Stephen A. Hall Philippe
More informationPROTECTION OF MANNED SPACECRAFT RADIATOR FROM SPACE DEBRIS
PROTECTION OF MANNED SPACECRAFT RADIATOR FROM SPACE DEBRIS Shan Li (1), Zheng Shigui (1), Yan Jun (1) (1) Institute of Spacecraft System Engineering, China Academy of Space Technology, Youyi Road 104,
More informationSTRESS TRANSPORT MODELLING 2
STRESS TRANSPORT MODELLING 2 T.J. Craft Department of Mechanical, Aerospace & Manufacturing Engineering UMIST, Manchester, UK STRESS TRANSPORT MODELLING 2 p.1 Introduction In the previous lecture we introduced
More information3D dynamic crack propagation analysis with PDS-FEM
JAMSTEC-R IFREE Special Issue, November 2009 3D dynamic crack propagation analysis with PDS-FEM 1*, Hide Sakaguchi 1, Kenji Oguni 2, Muneo Hori 2 discretizations. Keywords +81-45-778-5972 lalith@jamstec.go.jp
More informationCHAPTER 3 THE EFFECTS OF FORCES ON MATERIALS
CHAPTER THE EFFECTS OF FORCES ON MATERIALS EXERCISE 1, Page 50 1. A rectangular bar having a cross-sectional area of 80 mm has a tensile force of 0 kn applied to it. Determine the stress in the bar. Stress
More informationExpansion of circular tubes by rigid tubes as impact energy absorbers: experimental and theoretical investigation
Expansion of circular tubes by rigid tubes as impact energy absorbers: experimental and theoretical investigation M Shakeri, S Salehghaffari and R. Mirzaeifar Department of Mechanical Engineering, Amirkabir
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationNUMERICAL SIMULATION OF CONCRETE EXPOSED TO HIGH TEMPERATURE DAMAGE AND EXPLOSIVE SPALLING
NUMERICAL SIMULATION OF CONCRETE EXPOSED TO HIGH TEMPERATURE DAMAGE AND EXPLOSIVE SPALLING Prof. Joško Ožbolt 1 Josipa Bošnjak 1, Goran Periškić 1, Akanshu Sharma 2 1 Institute of Construction Materials,
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationNUMERICAL MODELLING AND DETERMINATION OF FRACTURE MECHANICS PARAMETERS FOR CONCRETE AND ROCK: PROBABILISTIC ASPECTS
NUMERICAL MODELLING AND DETERMINATION OF FRACTURE MECHANICS PARAMETERS FOR CONCRETE AND ROCK: PROBABILISTIC ASPECTS J. Carmeliet Catholic University of Leuven, Department of Civil Engineering, Belgium
More informationNumerical modelling of the fluid structure interaction using ALE and SPH: The Hydrodynamic Ram phenomenon.
Numerical modelling of the fluid structure interaction using ALE and SPH: The Hydrodynamic Ram phenomenon. D. Varas 1, J.A. Artero-Guerrero 1, J. Pernas-Sánchez 1 and J. López-Puente 1 1 Department of
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationProtective Structures
High Velocity Penetration/Perforation Using Coupled Smooth Particle Hydrodynamics-Finite Element Method by S. Swaddiwudhipong, M.J. Islam, Z.S. Liu Reprinted from International Journal of Protective Structures
More informationINTERNATIONAL JOURNAL OF APPLIED ENGINEERING RESEARCH, DINDIGUL Volume 2, No 1, 2011
Interlaminar failure analysis of FRP cross ply laminate with elliptical cutout Venkateswara Rao.S 1, Sd. Abdul Kalam 1, Srilakshmi.S 1, Bala Krishna Murthy.V 2 1 Mechanical Engineering Department, P. V.
More informationEAS 664/4 Principle Structural Design
UNIVERSITI SAINS MALAYSIA 1 st. Semester Examination 2004/2005 Academic Session October 2004 EAS 664/4 Principle Structural Design Time : 3 hours Instruction to candidates: 1. Ensure that this paper contains
More informationCompressive Residual Stress Optimization in Laser Peening of a Curved Geometry
Compressive Residual Stress Optimization in Laser Peening of a Curved Geometry Anoop Vasu and Ramana V. Grandhi Department of Mechanical and Materials Engineering, Wright State University, Dayton, OH,
More informationFracture Mechanics, Damage and Fatigue Linear Elastic Fracture Mechanics - Energetic Approach
University of Liège Aerospace & Mechanical Engineering Fracture Mechanics, Damage and Fatigue Linear Elastic Fracture Mechanics - Energetic Approach Ludovic Noels Computational & Multiscale Mechanics of
More informationGAME PHYSICS SECOND EDITION. дяййтаййг 1 *
GAME PHYSICS SECOND EDITION DAVID H. EBERLY дяййтаййг 1 * К AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO MORGAN ELSEVIER Morgan Kaufmann Publishers
More informationANSYS Explicit Dynamics Update. Mai Doan
ANSYS Explicit Dynamics Update Mai Doan Mai.Doan@ansys.com +1 512 687 9523 1/32 ANSYS Explicit Dynamics Update Outline Introduction Solve Problems that were Difficult or Impossible in the Past Structural
More informationGas Dynamics: Basic Equations, Waves and Shocks
Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks Susanne Höfner Susanne.Hoefner@fysast.uu.se Astrophysical Dynamics, VT 010 Gas Dynamics: Basic Equations, Waves and Shocks
More information