FEMxDEM double scale approach with second gradient regularization applied to granular materials modeling
|
|
- Annabel Melton
- 6 years ago
- Views:
Transcription
1 FEMxDEM double scale approach with second gradient regularization applied to granular materials modeling Albert Argilaga Claramunt Stefano Dal Pont Gaël Combe Denis Caillerie Jacques Desrues 16 december 2016 Albert Argilaga Claramunt FEMxDEM Modeling with 1 / 44
2 Outline 1 Introduction 2 3 Random field methods Parametric tests Engineering scale simulation 4 Albert Argilaga Claramunt FEMxDEM Modeling with 2 / 44
3 Motivation Bridge the gap between continuum and discrete approaches Continuum approaches Numerically efficient No scale limitations? Figure: Numerical modeling dilemma Discrete approaches No need of phenomenological laws Faithful representation of the microscale Albert Argilaga Claramunt FEMxDEM Modeling with 3 / 44
4 Multiscale coupling FEM: Lagamine (Liège University) A numerically obtained constitutive law is injected into the material points DEM with PBC (in-house 3SR) f n'/m' n m nʼ mʼ f n/m Albert Argilaga Claramunt FEMxDEM Modeling with 4 / 44
5 Introduction Discrete Element Model (DEM) Frictional cohesive contact laws Periodic Boundary Conditions Figure: DEM assembly preparation Albert Argilaga Claramunt FEMxDEM Modeling with 5 / 44
6 Introduction Discrete Element Model (DEM) )/ 11 2 ( Periodic Boundary Conditions Frictional cohesive contact laws (%) Figure: DEM biaxial response Figure: DEM assembly preparation Numerical homogenization 1 X ~m/n m σf = f (~r ~r n ) S (n,m) C Albert Argilaga Claramunt FEMxDEM Modeling with 5 / 44
7 State of the Art in FEMxDEM Early works (Kaneko, 2003; Miehe, 2010; Nitka, 2011) Study of anisotropy (Guo and Zhao, 2013; Nguyen, 2013) DEM cohesion (Nguyen, 2014) Non-local regularization (Liu, 2015) 3D microscale (Liu, 2015; Wang and Sun, 2016) DEM material heterogeneity (Shahin, 2016) Macroscale hidro-mechanical coupling (Wang and Sun, 2016; Guo and Zhao, 2016) Full micro-macro 3D (Guo and Zhao, 2016) Albert Argilaga Claramunt FEMxDEM Modeling with 6 / 44
8 Model performance FEMxDEM issues Bad convergence Computational time Generic issues Mesh dependency Albert Argilaga Claramunt FEMxDEM Modeling with 7 / 44
9 Model performance FEMxDEM issues Bad convergence Computational time Generic issues Mesh dependency FEMxDEM proposed solutions Alternative Operator Generic proposed solutions Regularization Albert Argilaga Claramunt FEMxDEM Modeling with 7 / 44
10 Outline 1 Introduction 2 3 Random field methods Parametric tests Engineering scale simulation 4 Albert Argilaga Claramunt FEMxDEM Modeling with 8 / 44
11 Link with classical geomechanics FEM codes We consider the quasi static evolution of a continuous medium which shares features withelastoplasticity and hypoplasticity. The constitutive equation provides the Cauchy stress σ f as a function of the deformation tensor F f in a loading step: ( ) F f σ f = S F f (1) Albert Argilaga Claramunt FEMxDEM Modeling with 9 / 44
12 Link with classical geomechanics FEM codes We consider the quasi static evolution of a continuous medium which shares features withelastoplasticity and hypoplasticity. The constitutive equation provides the Cauchy stress σ f as a function of the deformation tensor F f in a loading step: ( ) F f σ f = S F f (1) The non linear boundary value problem, is solved using the iterative Newton s method. This needs to differentiate S with respect to F f (1), the gradient of S is often called the consistent operator, it is denoted by C dσ f = S ( F f + df f ) S ( F f ) = C : df f + (2) Albert Argilaga Claramunt FEMxDEM Modeling with 9 / 44
13 Link with classical geomechanics FEM codes We consider the quasi static evolution of a continuous medium which shares features withelastoplasticity and hypoplasticity. The constitutive equation provides the Cauchy stress σ f as a function of the deformation tensor F f in a loading step: ( ) F f σ f = S F f (1) The non linear boundary value problem, is solved using the iterative Newton s method. This needs to differentiate S with respect to F f (1), the gradient of S is often called the consistent operator, it is denoted by C dσ f = S ( F f + df f ) S ( F f ) = C : df f + (2) Albert Argilaga Claramunt FEMxDEM Modeling with 9 / 44
14 Determination of the Consistent Operator Except for some cases for which closed-form expression of the tensor C can be found, the determining of C is performed numerically using a perturbation method: (C ijmn ) = S ( ij F f + ɛλ mn) ( S ij F f ) ɛ Assuming S to be differentiable, that determines C (3) Albert Argilaga Claramunt FEMxDEM Modeling with 10 / 44
15 Determination of the Consistent Operator Except for some cases for which closed-form expression of the tensor C can be found, the determining of C is performed numerically using a perturbation method: (C ijmn ) = S ( ij F f + ɛλ mn) ( S ij F f ) ɛ Assuming S to be differentiable, that determines C (3) Issues Bad convergence Computational time Mesh dependency Proposed solutions Alternative Operator Regularization Albert Argilaga Claramunt FEMxDEM Modeling with 10 / 44
16 Alternative s for FEMxDEM: Perturbation based operators Consistent Tangent (CTO) Auxiliary Elastic (AEO) Elastic operators DEM Quasi Static (DEMQO) PreStressed Truss-Like (PSTLO) Kruyt operators Kruyt Augmented (KAO) Upper bound Kruyt (UKO) Upper bound Corrected Kruyt (UCKO - UCKO 2DOF) Albert Argilaga Claramunt FEMxDEM Modeling with 11 / 44
17 Perturbation based operators Consistent Tangent (CTO) Classic perturbation approach: (C ijmn ) = S ( ij F f + ɛλ mn) ( S ij F f ) ɛ Auxiliary Elastic (AEO) (Nguyen, 2013) Average of the CTO over several macroscopic loading steps Albert Argilaga Claramunt FEMxDEM Modeling with 12 / 44
18 Elastic operators DEM Quasi Static (DEMQO) Based on numerical homogenization of an elastic DEM (C ijmn ) DEMQO = SE ij (ɛλ mn ) ɛ (4) PreStressed Truss-Like (PSTLO) (Caillerie) It treats the DEM as a prestressed linearelastic truss with additional rotational degrees of freedom in the nodes dσ = 1 Y c C δcd i f c Y ( ) i + σ F T df T F T : df σ (5) Albert Argilaga Claramunt FEMxDEM Modeling with 13 / 44
19 Kruyt operators Kruyt Augmented (KAO) (Caillerie) UKO with prestresses and rotational degrees of reedom Upper bound Kruyt (UKO) (Kruyt & Rothenburg, 1998) 1 Y a (l c ) 2 ( k n ( e c e c ) ( e c e c ( ) + k t t c e c) ( t c e c)) c C Upper bound Corrected Kruyt (UCKO - UCKO 2DOF) Calibrations of the UKO to fit a CTO (6) Albert Argilaga Claramunt FEMxDEM Modeling with 14 / 44
20 Convergence rates The convergence rate of the operators is compared in the initial state and post-peak in a biaxial compression test: OTC UKO UCKO UCKO 2DOF DEMQO/PTLO KAO AEO UKO UCKO UCKO 2DOF DEMQO/PTLO KAO Residual FNORM/RNORM 10 2 Residual FNORM/RNORM iteration iteration Albert Argilaga Claramunt FEMxDEM Modeling with 15 / 44
21 Operators performance Perturbation based operators perform well in the initial stages Kruyt operators are more robust in the post-peak part Elastic operators combine good performance before and after the stress peak # iterations OTC/ AEO UKO UCKO UCKO 2DOF DEMQO /PSTLO Initial Post-peak Figure: Operators comparison: number of iterations KAO Albert Argilaga Claramunt FEMxDEM Modeling with 16 / 44
22 Conclusion The elastic operators overcome the issues giving good stability and convergence velocity # iterations OTC/ AEO UKO UCKO UCKO 2DOF DEMQO /PSTLO Initial Post-peak KAO Figure: Operators comparison: number of iterations Albert Argilaga Claramunt FEMxDEM Modeling with 17 / 44
23 Outline 1 Introduction 2 3 Random field methods Parametric tests Engineering scale simulation 4 Albert Argilaga Claramunt FEMxDEM Modeling with 18 / 44
24 FEMxDEM drawbacks The evaluation of S is issued from a DEM homogenization: F f σ f = S (F ) f Albert Argilaga Claramunt FEMxDEM Modeling with 19 / 44
25 FEMxDEM drawbacks The evaluation of S is issued from a DEM homogenization: F f σ f = S (F ) f Issues Bad convergence Computational time Mesh dependency Albert Argilaga Claramunt FEMxDEM Modeling with 19 / 44
26 FEMxDEM drawbacks The evaluation of S is issued from a DEM homogenization: F f σ f = S (F ) f Issues Bad convergence Computational time Mesh dependency Proposed solutions Alternative Regularization Albert Argilaga Claramunt FEMxDEM Modeling with 19 / 44
27 Element loop parallelization Amdahl s law (Amdahl, 1967) Sequential Parallel S(n) = 1 B + 1 n (1 B) (7) lim S(n) = n (8) B 0 S: speedup n: number of cores B: non parallelizable part Walltime (min) x Intel Xeon CPU E GHz (2 x 4 cores) 2x Intel Xeon CPU E GHz (2 x 6 cores) Other Approaches: Massive : MPI (Desrues) Figure: speedup: x7.15 (8 cores), x9.06 (12 cores) Albert Argilaga Claramunt FEMxDEM Modeling with 20 / 44
28 Numerical randomness Theorem Commutative Law of Addition: a + b = b + a Albert Argilaga Claramunt FEMxDEM Modeling with 21 / 44
29 Numerical randomness Does not apply! Commutative Law of Addition: a + b = b + a Albert Argilaga Claramunt FEMxDEM Modeling with 21 / 44
30 Numerical randomness Does not apply! Commutative Law of Addition: a + b = b + a Figure: Simulated numerical error Albert Argilaga Claramunt FEMxDEM Modeling with 21 / 44
31 Numerical randomness Does not apply! Commutative Law of Addition: a + b = b + a Figure: Simulated numerical error Figure: Different possible solutions of the same BVP. Cumulative second invariant of strain Albert Argilaga Claramunt FEMxDEM Modeling with 21 / 44
32 Conclusion The parallelization provides an effective speedup of the model, becoming this competitive with classical FEM codes Albert Argilaga Claramunt FEMxDEM Modeling with 22 / 44
33 Outline 1 Introduction 2 3 Random field methods Parametric tests Engineering scale simulation 4 Albert Argilaga Claramunt FEMxDEM Modeling with 23 / 44
34 First order FEM * E- 03 * E Figure: Non objectivity of the model Albert Argilaga Claramunt FEMxDEM Modeling with 24 / 44
35 First order constitutive relation Balance equation in virtual power formulation: ui, σ ij ɛ ijdω = t i ui ds (9) Ω Where ui is the kinematically admissible virtual displacement field, σ ij the stress field, ɛ ij the virtual strain field, Ω the surface of Ω, t i the surface forces Ω Albert Argilaga Claramunt FEMxDEM Modeling with 25 / 44
36 First order constitutive relation Balance equation in virtual power formulation: ui, σ ij ɛ ijdω = t i ui ds (9) Ω Where ui is the kinematically admissible virtual displacement field, σ ij the stress field, ɛ ij the virtual strain field, Ω the surface of Ω, t i the surface forces Issues Bad convergence Computational time Mesh dependency Ω Albert Argilaga Claramunt FEMxDEM Modeling with 25 / 44
37 First order constitutive relation Balance equation in virtual power formulation: ui, σ ij ɛ ijdω = t i ui ds (9) Ω Where ui is the kinematically admissible virtual displacement field, σ ij the stress field, ɛ ij the virtual strain field, Ω the surface of Ω, t i the surface forces Ω Issues Bad convergence Computational time Mesh dependency Proposed solutions Alternative Regularization Albert Argilaga Claramunt FEMxDEM Modeling with 25 / 44
38 microstructured materials Second gradient regularization introduces an internal length that regularizes the model First order term ui, Ω σ ij ɛ 2 u i ij + Σ ijk dω = t i ui ds x j x k Ω where Σ ijk is the double stress dual of 2 ui x j x k Albert Argilaga Claramunt FEMxDEM Modeling with 26 / 44
39 microstructured materials Second gradient regularization introduces an internal length that regularizes the model First order term ui, Ω Second order term σ ij ɛ 2 u i ij + Σ ijk dω = t i ui ds x j x k Ω where Σ ijk is the double stress dual of 2 ui x j x k Albert Argilaga Claramunt FEMxDEM Modeling with 26 / 44
40 in a Boundary Value Problem As it is a local model, it can be injected in the material points in the same way as the first gradient relation is! Albert Argilaga Claramunt FEMxDEM Modeling with 27 / 44
41 * E * E * E * E * E * E Introduction Biaxial test enriched with Figure: 128, 512, and 2048 elements biaxial compression tests Albert Argilaga Claramunt FEMxDEM Modeling with 28 / 44
42 With the enrichment the model becomes objective. This allows to use any mesh refinement The code is ready to be used in real scale problems Albert Argilaga Claramunt FEMxDEM Modeling with 29 / 44
43 Random field methods Parametric tests Engineering scale simulation Outline 1 Introduction 2 3 Random field methods Parametric tests Engineering scale simulation 4 Albert Argilaga Claramunt FEMxDEM Modeling with 30 / 44
44 Random field methods Parametric tests Engineering scale simulation Introduction: Symmetry breaking Material properties in numerical models are spatially homogeneous unless the contrary is defined The use of a catalyzer to break the symmetry can improve the model performance, both physically and numerically Figure: A CT scan of a sand sample (Andò, 2013) Albert Argilaga Claramunt FEMxDEM Modeling with 31 / 44
45 Random field methods Parametric tests Engineering scale simulation Punctual imperfection The DEM microstructure has a strong influence on the results Figure: Homogeneous (a) Defect (b). Second invariant of strain Albert Argilaga Claramunt FEMxDEM Modeling with 32 / 44
46 Random field methods Parametric tests Engineering scale simulation Full field DEM random variability Figure: (a) FEMxDEM and pure DEM responses under a biaxial loading. (b) and (c), cumulative deviatoric strain field (Shahin, 2016) Albert Argilaga Claramunt FEMxDEM Modeling with 33 / 44
47 Random field methods Parametric tests Engineering scale simulation The FEMxDEM approach provides a natural way to embed material heterogeneity into the model Albert Argilaga Claramunt FEMxDEM Modeling with 34 / 44
48 Random field methods Parametric tests Engineering scale simulation Outline 1 Introduction 2 3 Random field methods Parametric tests Engineering scale simulation 4 Albert Argilaga Claramunt FEMxDEM Modeling with 35 / 44
49 Random field methods Parametric tests Engineering scale simulation Outer radius The outer radius has an influence on the results, a compromise size is chosen Figure: Outer radius=50, 100 and 200 m. Cumulative second invariant of strain Figure: Radius=50, 100 and 200 m M.S. research project Nurisa, 2015 Albert Argilaga Claramunt FEMxDEM Modeling with 36 / 44
50 Random field methods Parametric tests Engineering scale simulation Micro/macro DEM allows to tune the coordination independently from the density The model can cope with any macroscopic boundary conditions Figure: Coordination number = 4.10, 3.67 and Cumulative second invariant of strain Figure: Stress ratio σ 0H : σ 0V = 1:1, 1:1.3 and 1:2. Cumulative second invariant of strain Albert Argilaga Claramunt FEMxDEM Modeling with 37 / 44
51 q Introduction Random field methods Parametric tests Engineering scale simulation Intrinsic anisotropy X: Y: X: Y: Figure: Precharge 0%, 0.3% and 0.6%. Cumulative second invariant of strain ε yy Figure: DEM preparation Albert Argilaga Claramunt FEMxDEM Modeling with 38 / 44
52 Random field methods Parametric tests Engineering scale simulation Outline 1 Introduction 2 3 Random field methods Parametric tests Engineering scale simulation 4 Albert Argilaga Claramunt FEMxDEM Modeling with 39 / 44
53 Random field methods Parametric tests Engineering scale simulation Hollow cylinder with * 1.000E-04 Y X Figure: Localization around a gallery. FEM: 4950 elements, nodes, DOF Figure: 400 particles DEM microscale at the end of the loading Albert Argilaga Claramunt FEMxDEM Modeling with 40 / 44
54 Random field methods Parametric tests Engineering scale simulation Real scale simulations are possible thanks to the second gradient and parallelization FEMxDEM allows to calibrate the model directly from the microscale Albert Argilaga Claramunt FEMxDEM Modeling with 41 / 44
55 General A FEMxDEM approach allows to use a constitutive law built via numerical homogenization in the DEM microscale An alternative Newton operator provides the model with good convergence rate and robustness allows the model to compete with classical FEM approaches in terms of walltime A Second gradient regularization assures the mesh independence Real scale problems with fine meshes can be considered Albert Argilaga Claramunt FEMxDEM Modeling with 42 / 44
56 Perspectives Calibration of the microscale with experimental data 3D microscale can bring a richer constitutive behabiour More complex microscale with massive parallelization (MPI) Albert Argilaga Claramunt FEMxDEM Modeling with 43 / 44
57 Thanks!!! Albert Argilaga Claramunt FEMxDEM Modeling with 44 / 44
Hierarchical Multiscale Modeling of Strain Localization in Granular Materials: ACondensedOverviewandPerspectives
Hierarchical Multiscale Modeling of Strain Localization in Granular Materials: ACondensedOverviewandPerspectives Jidong Zhao Abstract This paper presents a brief overview on hierarchical multiscale modeling
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 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 informationMicromechanics of granular materials: slow flows
Micromechanics of granular materials: slow flows Niels P. Kruyt Department of Mechanical Engineering, University of Twente, n.p.kruyt@utwente.nl www.ts.ctw.utwente.nl/kruyt/ 1 Applications of granular
More informationMODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING EPS SEMINAR SERIES MARCH 2008
MODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING EPS SEMINAR SERIES MARCH 2008 COLLABORATORS: DR XUXIN TU AND MR KIRK ELLISON THE ROADMAP MOTIVATION
More informationModelling the excavation damaged zone in Callovo-Oxfordian claystone with strain localisation
Modelling the excavation damaged zone in Callovo-Oxfordian claystone with strain localisation B. Pardoen - F. Collin - S. Levasseur - R. Charlier Université de Liège ArGEnCo ALERT Workshop 2012 Aussois,
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 informationFormulation and numerical implementation of micro-scale boundary conditions for particle aggregates
Granular Matter (2017) 19:72 DOI 10.1007/s10035-017-0757-y ORIGINAL PAPER Formulation and numerical implementation of micro-scale boundary conditions for particle aggregates J. Liu 1 E. Bosco 1 A. S. J.
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 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 informationExcavation Damaged Zone Modelling in Claystone with Coupled Second Gradient Model
Excavation Damaged Zone Modelling in Claystone with Coupled Second Gradient Model Frédéric Collin * and Benoît Pardoen ** Argenco Department University of Liège 1, Chemin des Chevreuils B-4000 Liège, Belgium
More informationA coupled FEM/DEM approach for hierarchical multiscale modelling of granular media
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2014; 99:789 818 Published online 13 June 2014 in Wiley Online Library (wileyonlinelibrary.com)..4702 A coupled FEM/DEM
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 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 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 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 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 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 informationTable of Contents. Foreword... xiii Introduction... xv
Foreword.... xiii Introduction.... xv Chapter 1. Controllability of Geotechnical Tests and their Relationship to the Instability of Soils... 1 Roberto NOVA 1.1. Introduction... 1 1.2. Load control... 2
More informationMulti-scale digital image correlation of strain localization
Multi-scale digital image correlation of strain localization J. Marty a, J. Réthoré a, A. Combescure a a. Laboratoire de Mécanique des Contacts et des Strcutures, INSA Lyon / UMR CNRS 5259 2 Avenue des
More informationA model order reduction technique for speeding up computational homogenisation
A model order reduction technique for speeding up computational homogenisation Olivier Goury, Pierre Kerfriden, Wing Kam Liu, Stéphane Bordas Cardiff University Outline Introduction Heterogeneous materials
More informationNUMERICAL MODELING OF INSTABILITIES IN SAND
NUMERICAL MODELING OF INSTABILITIES IN SAND KIRK ELLISON March 14, 2008 Advisor: Jose Andrade Masters Defense Outline of Presentation Randomized porosity in FEM simulations Liquefaction in FEM simulations
More informationTowards multiscale computation of confined granular media - Contact forces, stresses and tangent operators -
TECHNISCHE MECHANIK, Band 6, Heft, (996, 77 88 Manuskripteingang: 4. Dezember 995 Towards multiscale computation of confined granular media - Contact forces, stresses and tangent operators - H. A. Meier,
More informationStructural Analysis of Truss Structures using Stiffness Matrix. Dr. Nasrellah Hassan Ahmed
Structural Analysis of Truss Structures using Stiffness Matrix Dr. Nasrellah Hassan Ahmed FUNDAMENTAL RELATIONSHIPS FOR STRUCTURAL ANALYSIS In general, there are three types of relationships: Equilibrium
More informationCourse in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg
Course in Nonlinear FEM Geometric nonlinearity Nonlinear FEM Outline Lecture 1 Introduction Lecture 2 Geometric nonlinearity Lecture 3 Material nonlinearity Lecture 4 Material nonlinearity it continued
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 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 informationAdaptive Analysis of Bifurcation Points of Shell Structures
First published in: Adaptive Analysis of Bifurcation Points of Shell Structures E. Ewert and K. Schweizerhof Institut für Mechanik, Universität Karlsruhe (TH), Kaiserstraße 12, D-76131 Karlsruhe, Germany
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 informationLoad Cell Design Using COMSOL Multiphysics
Load Cell Design Using COMSOL Multiphysics Andrei Marchidan, Tarah N. Sullivan and Joseph L. Palladino Department of Engineering, Trinity College, Hartford, CT 06106, USA joseph.palladino@trincoll.edu
More informationConstitutive models: Incremental plasticity Drücker s postulate
Constitutive models: Incremental plasticity Drücker s postulate if consistency condition associated plastic law, associated plasticity - plastic flow law associated with the limit (loading) surface Prager
More informationVariational principles in mechanics
CHAPTER Variational principles in mechanics.1 Linear Elasticity n D Figure.1: A domain and its boundary = D [. Consider a domain Ω R 3 with its boundary = D [ of normal n (see Figure.1). The problem of
More informationHomogenization method for elastic materials
Outline 1 Introduction 2 Description of the geometry 3 Problem setting 4 Implementation & Results 5 Bandgaps 6 Conclusion Introduction Introduction study of the homogenization method applied on elastic
More informationLectures on. Constitutive Modelling of Arteries. Ray Ogden
Lectures on Constitutive Modelling of Arteries Ray Ogden University of Aberdeen Xi an Jiaotong University April 2011 Overview of the Ingredients of Continuum Mechanics needed in Soft Tissue Biomechanics
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 informationViscoelastic Damping Characteristics of Indium-Tin/SiC Particulate Composites
Viscoelastic Damping Characteristics of Indium-Tin/SiC Particulate Composites HyungJoo Kim, Colby C. Swan Ctr. for Computer-Aided Design, Univ. of Iowa Roderic Lakes Engineering Physics, Univ. of Wisconsin
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationUNIVERSITÄT FÜR BODENKULTUR WIEN
UNIVERSITÄT FÜR BODENKULTUR WIEN Department für Bautechnik und Naturgefahren Institut für Geotechnik Dissertation Linking DEM with micropolar continuum Jia Lin April 2010 - March 2013 UNIVERSITÄT FÜR BODENKULTUR
More informationStress and fabric in granular material
THEORETICAL & APPLIED MECHANICS LETTERS 3, 22 (23) Stress and fabric in granular material Ching S. Chang,, a) and Yang Liu 2 ) Department of Civil Engineering, University of Massachusetts Amherst, Massachusetts
More informationContinuum Models of Discrete Particle Systems with Particle Shape Considered
Introduction Continuum Models of Discrete Particle Systems with Particle Shape Considered Matthew R. Kuhn 1 Ching S. Chang 2 1 University of Portland 2 University of Massachusetts McMAT Mechanics and Materials
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 informationShock enhancement of cellular structures under impact loading: Part II Analysis
Shock enhancement of cellular structures under impact loading: Part II Analysis S. Pattofatto, 1 I. Elnasri, 1 H. Zhao, 1* H. Tsitsiris 1 F. Hild 1 and Y. Girard 2 1 Laboratoire de Mécanique et Technologie
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationConstitutive modelling of fabric anisotropy in sand
Geomechanics from Micro to Macro Soga et al. (Eds) 2015 Taylor & Francis Group, London, ISBN 978-1-138-02707-7 Constitutive modelling of fabric anisotropy in sand Z.W. Gao School of Engineering, University
More informationCSMA Towards optimal design of multiscale nonlinear structures and reduced-order modeling approaches. 1 Introduction.
CSMA 2017 13ème Colloque National en Calcul des Structures 15-19 Mai 2017, Presqu île de Giens (Var) Towards optimal design of multiscale nonlinear structures and reduced-order modeling approaches Liang
More information1 Nonlinear deformation
NONLINEAR TRUSS 1 Nonlinear deformation When deformation and/or rotation of the truss are large, various strains and stresses can be defined and related by material laws. The material behavior can be expected
More informationNiiranen, Jarkko; Khakalo, Sergei; Balobanov, Viacheslav Strain gradient elasticity theories in lattice structure modelling
Powered by TCPDF (www.tcpdf.org) This is an electronic reprint of the original article. This reprint may differ from the original in pagination and typographic detail. Niiranen, Jarkko; Khakalo, Sergei;
More informationContinuum mechanics of beam-like structures using one-dimensional finite element based on Serendipity Lagrange cross-sectional discretisation, Mayank Patni, Prof. Paul Weaver, Dr Alberto Pirrera Bristol
More informationVirtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing example
AMAS Workshop on Smart Materials and Structures SMART 03 (pp.313 324) Jadwisin, September 2-5, 2003 Virtual distortions applied to structural modelling and sensitivity analysis. Damage identification testing
More informationFig. 1. Circular fiber and interphase between the fiber and the matrix.
Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In
More informationKun Wang, 1 WaiChing Sun, 1, Simon Salager, 2 SeonHong Na, 1 Ghonwa Khaddour 2. Mudd, Mail Code: 4709, New York, New York 10027, USA
International Journal for Multiscale Computational Engineering, 14 (4): 389 413 (2016) IDENTIFYING MATERIAL PARAMETERS FOR A MICRO-POLAR PLASTICITY MODEL VIA X-RAY MICRO-COMPUTED TOMOGRAPHIC (CT) IMAGES:
More informationMAXIMUM ENTROPY-BASED UNCERTAINTY MODELING AT THE FINITE ELEMENT LEVEL. Pengchao Song and Marc P. Mignolet
MAXIMUM ENTROPY-BASED UNCERTAINTY MODELING AT THE FINITE ELEMENT LEVEL Pengchao Song and Marc P. Mignolet SEMTE, Faculties of Mechanical and Aerospace Engineering, Arizona State University, 51 E. Tyler
More informationEntropy and Material Instability in the Quasi-static Mechanics of Granular Media
Entropy and Material Instability in the Quasi-static Mechanics of Granular Media J.D. Goddard Department of Aerospace & Mechanical Engineering University of California, San Diego, USA jgoddard@ucsd.edu
More informationMicromechanism Investigation of Geomaterials using Particle Manifold Method (PMM)
NUMERICS in GEOTECHNICS and STRUCTURES 2012, 27 years ZSOIL.PC Micromechanism Investigation of Geomaterials using Particle Manifold Method (PMM) by Liang SUN Laboratory of Rock Mechanics(LMR), EPFL August
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 informationA primer on Numerical methods for elasticity
A primer on Numerical methods for elasticity Douglas N. Arnold, University of Minnesota Complex materials: Mathematical models and numerical methods Oslo, June 10 12, 2015 One has to resort to the indignity
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 informationVariable Definition Notes & comments
Extended base dimension system Pi-theorem (also definition of physical quantities, ) Physical similarity Physical similarity means that all Pi-parameters are equal Galileo-number (solid mechanics) Reynolds
More informationNumerical analyses of cement-based piezoelectric smart composites
Numerical analyses of cement-based piezoelectric smart composites *Jan Sladek 1, Pavol Novak 2, Peter L. Bishay 3, and Vladimir Sladek 1 1 Institute of Construction and Architecture, Slovak Academy of
More informationEshan V. Dave, Secretary of M&FGM2006 (Hawaii) Research Assistant and Ph.D. Candidate. Glaucio H. Paulino, Chairman of M&FGM2006 (Hawaii)
Asphalt Pavement Aging and Temperature Dependent Properties through a Functionally Graded Viscoelastic Model I: Development, Implementation and Verification Eshan V. Dave, Secretary of M&FGM2006 (Hawaii)
More informationPowerful Modelling Techniques in Abaqus to Simulate
Powerful Modelling Techniques in Abaqus to Simulate Necking and Delamination of Laminated Composites D. F. Zhang, K.M. Mao, Md. S. Islam, E. Andreasson, Nasir Mehmood, S. Kao-Walter Email: sharon.kao-walter@bth.se
More informationMicrostructural Randomness and Scaling in Mechanics of Materials. Martin Ostoja-Starzewski. University of Illinois at Urbana-Champaign
Microstructural Randomness and Scaling in Mechanics of Materials Martin Ostoja-Starzewski University of Illinois at Urbana-Champaign Contents Preface ix 1. Randomness versus determinism ix 2. Randomness
More informationThe Finite Element Method for Computational Structural Mechanics
The Finite Element Method for Computational Structural Mechanics Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) January 29, 2010 Martin Kronbichler (TDB) FEM for CSM January
More informationDiscrete Element Modelling of a Reinforced Concrete Structure
Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr
More informationThe Stress Variations of Granular Samples in Direct Shear Tests using Discrete Element Method
The Stress Variations of Granular Samples in Direct Shear Tests using Discrete Element Method Hoang Khanh Le 1), *Wen-Chao Huang 2), Yi-De Zeng 3), Jheng-Yu Hsieh 4) and Kun-Che Li 5) 1), 2), 3), 4), 5)
More informationMechanics of materials Lecture 4 Strain and deformation
Mechanics of materials Lecture 4 Strain and deformation Reijo Kouhia Tampere University of Technology Department of Mechanical Engineering and Industrial Design Wednesday 3 rd February, 206 of a continuum
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 informationParticle flow simulation of sand under biaxial test
5th International Conference on Civil Engineering and Transportation (ICCET 2015) Particle flow simulation of sand under biaxial test Xiao-li Dong1,2, a *,Wei-hua Zhang1,a 1 Beijing City University, China
More informationBasic concepts to start Mechanics of Materials
Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen
More informationMASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK
MASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK J. Pina-Henriques and Paulo B. Lourenço School of Engineering, University of Minho, Guimarães, Portugal Abstract Several continuous and discontinuous
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 informationMicro-macro modelling for fluids and powders
Micro-macro modelling for fluids and powders Stefan Luding 1,2 1 Particle Technology, DelftChemTech, TU Delft, Julianalaan 136, 2628 BL Delft, The Netherlands 2 e-mail: s.luding@tnw.tudelft.nl ABSTRACT
More informationModelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach
Modelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach Giang Dinh Nguyen, Chi Thanh Nguyen, Vinh Phu Nguyen School of Civil, Environmental and
More informationComputational homogenization of material layers with micromorphic mesostructure
Computational homogenization of material layers with micromorphic mesostructure C. B. Hirschberger, N. Sukumar, P. Steinmann Manuscript as accepted for publication in Philosophical Magazine, 21 September
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
More informationComputational Homogenization and Multiscale Modeling
Computational Homogenization and Multiscale Modeling Kenneth Runesson and Fredrik Larsson Chalmers University of Technology, Department of Applied Mechanics Runesson/Larsson, Geilo 2011-01-24 p.1/23 Lecture
More informationNonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess
Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable
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 informationFEM for elastic-plastic problems
FEM for elastic-plastic problems Jerzy Pamin e-mail: JPamin@L5.pk.edu.pl With thanks to: P. Mika, A. Winnicki, A. Wosatko TNO DIANA http://www.tnodiana.com FEAP http://www.ce.berkeley.edu/feap Lecture
More informationMacroscopic Failure Analysis Based on a Random Field Representations Generated from Material Microstructures
Macroscopic Failure Analysis Based on a Random Field Representations Generated from Material Microstructures Reza Abedi Mechanical, Aerospace & Biomedical Engineering University of Tennessee Knoxville
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 1-21 September, 2017 Institute of Structural Engineering
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 informationNumerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions
TECHNISCHE MECHANIK, 33, 2, (2013), 97 103 submitted: December 11, 2012 Numerical Properties of Spherical and Cubical Representative Volume Elements with Different Boundary Conditions R. Glüge, M. Weber
More informationAfter lecture 16 you should be able to
Lecture 16: Design of paper and board packaging Advanced concepts: FEM, Fracture Mechanics After lecture 16 you should be able to describe the finite element method and its use for paper- based industry
More informationAn Evaluation of Simplified Methods to Compute the Mechanical Steady State
An Evaluation of Simplified Methods to Compute the Mechanical Steady State T. Herbland a,b, G. Cailletaud a, S. Quilici a, H. Jaffal b, M. Afzali b a Mines Paris Paris Tech, CNRS UMR 7633, BP 87, 91003
More informationAn experimental assessment of displacement fluctuations in a 2D granular material subjected to shear
Written: 14/6/212 An experimental assessment of displacement fluctuations in a 2D granular material subjected to shear Vincent Richefeu (1), Gaël Combe (1), Gioacchino Viggiani (1) Email Address of Submitting
More informationComputational models of diamond anvil cell compression
UDC 519.6 Computational models of diamond anvil cell compression A. I. Kondrat yev Independent Researcher, 5944 St. Alban Road, Pensacola, Florida 32503, USA Abstract. Diamond anvil cells (DAC) are extensively
More informationA TWO-SCALE FAILURE MODEL FOR HETEROGENEOUS MATERIALS: NUMERICAL IMPLEMENTATION BASED ON THE FINITE ELEMENT METHOD
A TWO-SCALE FAILURE MODEL FOR HETEROGENEOUS MATERIALS: NUMERICAL IMPLEMENTATION BASED ON THE FINITE ELEMENT METHOD S. Toro 1,2, P.J. Sánchez 1,2, A.E. Huespe 1,3, S.M. Giusti 4, P.J. Blanco 5,6, R.A. Feijóo
More informationExploiting pattern transformation to tune phononic band gaps in a two-dimensional granular crystal
Exploiting pattern transformation to tune phononic band gaps in a two-dimensional granular crystal The Harvard community has made this article openly available. Please share how this access benefits you.
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
More informationConcrete Fracture Prediction Using Virtual Internal Bond Model with Modified Morse Functional Potential
Concrete Fracture Prediction Using Virtual Internal Bond Model with Modified Morse Functional Potential Kyoungsoo Park, Glaucio H. Paulino and Jeffery R. Roesler Department of Civil and Environmental Engineering,
More informationINVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA
Problems in Solid Mechanics A Symposium in Honor of H.D. Bui Symi, Greece, July 3-8, 6 INVERSE ANALYSIS METHODS OF IDENTIFYING CRUSTAL CHARACTERISTICS USING GPS ARRYA DATA M. HORI (Earthquake Research
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationMicroplane Model formulation ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary
Microplane Model formulation 2010 ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary Table of Content Engineering relevance Theory Material model input in ANSYS Difference with current concrete
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationMITOCW MITRES2_002S10nonlinear_lec15_300k-mp4
MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More informationProceedings of the ASME th International Conference on Ocean, Offshore and Arctic Engineering OMAE2016 June 19-24, 2016, Busan, South Korea
Proceedings of the ASME 26 35th International Conference on Ocean, Offshore and Arctic Engineering OMAE26 June 9-24, 26, Busan, South Korea OMAE26-54554 LOCAL STRAIN AND STRESS CALCULATION METHODS OF IRREGULAR
More informationFE-studies of a Deterministic and Statistical Size Effect in Granular Bodies Including Shear Localization
The 1 th International Conference of International Association for Computer Methods and Advances in Geomechanics (IACMAG) 1- October, Goa, India FE-studies of a Deterministic and Statistical Size Effect
More informationELASTOPLASTICITY THEORY by V. A. Lubarda
ELASTOPLASTICITY THEORY by V. A. Lubarda Contents Preface xiii Part 1. ELEMENTS OF CONTINUUM MECHANICS 1 Chapter 1. TENSOR PRELIMINARIES 3 1.1. Vectors 3 1.2. Second-Order Tensors 4 1.3. Eigenvalues and
More information