MOOSE Workshop MOOSE-PF/MARMOT Training Mechanics Part Nashville, February 14th Daniel Schwen
|
|
- Meghan Phelps
- 5 years ago
- Views:
Transcription
1 MOOSE Workshop MOOSE-PF/MARMOT Training Mechanics Part Nashville, February 14th Daniel Schwen
2 III. Mechanics
3 Mechanics Outline 1. Introduction 2. Tensor-based System 1. Compute Elasticity Tensor 2. Compute Strain 3. Compute Stress 3. Linear Elasticity 4. Finite Strain 1. Elasticity 2. J2 Plasticity 3. Crystal Plasticity 5. Eigen Strains 6. Contact
4 Modeling Mechanical Deformation The mechanical response of materials is a critical factor for all structural materials and many functional materials. Solid mechanics predicts deformation and motion under the action of forces, temperature changes, phase changes, and other external or internal agents. Inputs: material properties, boundary conditions, body forces Outputs: displacements, strains, stresses
5 Solid Mechanics Without going into any detail, we solve for the displacement vector u with the residual equation (with no acceleration): r + b = 0 in u = g in n = in g Where σ is the Cauchy stress tensor, u is the displacement vector, b is the body force, n is the unit normal to the boundary, g is the prescribed displacement on the boundary and ι is the prescribed traction on the boundary. The weak form of the residual is expressed as R =(, r m ) h, mi (b, m) =0 For a given material behavior, the constitutive model describes the stress response of the material, i.e. The constitutive model could be linear (linear elasticity) or nonlinear for both small and finite strains.
6 Module provides the tools necessary for solving solid mechanics using tensor forms, with emphasis on developing new models Anisotropic elasticity tensors that can change spatially Linear elasticity Eigen strains Finite strain mechanics J2 plasticity Crystal plasticity
7 Documentation The tensor mechanics module is documented on the mooseframework.org wiki:
8 Tensor Mechanics Objects To make the development of materials models as straight forward as possible, we represent stress, elasticity tensors, etc. as tensors. = C Thus, is written as: stress = elasticity_tensor*strain Rank two tensors are rotated according to stress_rot = R*stress*R.transpose() RankTwoTensor (stress, strain) /tensor_mechanics/src/utils/ranktwotensor.c 3 by 3, with indices accessed as A(i,j) Addition, subtraction, multiplication, transpose, inverse, etc. are all defined. RankFourTensor /tensor_mechanics/src/utils/rankfourtensor.c 3 by 3 by 3 by 3, with indices accessed as C(i,j,k,l) Addition, subtraction with RankFourTensors Muliplication defined with RankFourTensors and RankTwoTensor, as well as rotation. Values set assuming various levels of symmetry ElasticityTensorR4 /tensor_mechanics/src/utils/elasticitytensorr4.c Specialization of RankFourTensor Can be inverted for nonlinear solution RotationTensor /tensor_mechanics/src/utils/rotationtensor.c Defined in terms of three Euler angles
9 Tensor Mechanics Solution Vector variables do not exist in MOOSE, so three (or two) variables must be created, typically disp_x, disp_y and disp_z. Three kernels would also be necessary, however we have a custom kernel action called TensorMechanics that creates the kernels. [GlobalParams] displacements = disp_x disp_y disp_z [] [Kernels] [./TensorMechanics] [../] [] The kernel is defined in StressDivergenceTensors Weak form of the stress divergence equation is: Real StressDivergenceTensors::computeQpResidual() { return _stress[_qp].row(_component) * _grad_test[_i][_qp]; } //Taken from StressDivegenceTensors in MOOSE-tensor_mechanics (, r i )
10 Calculating the Stress The stress, used in the mechanics kernel, is calculated using the MOOSE material system. The elasticity tensor, strain, and stress are each calculated in their own block. Thus, three blocks are required in the input file: [Materials] [./elasticity_tensor] type = ComputeElasticityTensor block = 1 fill_method = symmetric_isotropic C_ijkl = '4.2e6 1.9e7' [../] [./strain] type = ComputeSmallStrain block = 1 # displacements in GlobalParams [../] [./stress] type = ComputeLinearElasticStress block = 1 [../] []
11 ComputeElasticityTensor This material creates the elasticity tensor from an input vector. Various symmetries are possible, defined by the fill_method Antisymmetric_isotropic Takes 1 input Symmetric_isotropic Takes 2 inputs [1 st Lame (bulk), 2 nd Lame (shear)] General_isotropic Takes 3 inputs [1 st Lame, 2 nd Lame, antisymmetric shear] Antisymmetric Takes 6 inputs [B 1212 B 1213 B 1223 B 1313 B 1323 B 2323 ] Symmetric9 Takes 9 inputs [C 1111 C 1122 C 1133 C 2222 C 2233 C 3333 C 2323 C 1313 C 1212 ] Symmetric21 Takes 21 (see.h file for list) General Assigns all 81 components Optionally, you can also pass in three Euler angles to rotate the elasticity tensor (see ComputeRotatedElasticityTensorBase) Users can define their own ways of calculating elasticity tensors by inheriting from ComputeElasticityTensorBase ComputeIsotropicElasticityTensor ComputeConcentrationDependentElasticityTensor CompositeElasticityTensor Provides derivatives!
12 ComputeIsotropicElasticityTensor Use ComputeIsotropicElasticityTensor For isotropy, there are various ways of defining the elastic constants This material creates an isotropic elasticity tensor from three possible combination of elastic constants. Possible combinations: First lame constant λ and shear modulus µ Young s modulus E and Poisson s ratio ν Bulk modulus K and shear modulus µ [./elasticity_tensor] type = ComputeIsotropicElasticityTensor block = 1 youngs_modulus = 240e9 #GPa poissons_ratio = 0.3 [../]
13 CompositeElasticityTensor Assemble arbitrary tensors with non-linear variable dependencies Combine a set of base tensors C i with a list of functions f i C = i f i C i The f i are Derivative Function Materials (recall Phase Field) All derivatives are built automatically [./C1_tensor] type = ComputeIsotropicElasticityTensor block = 1 base_name = C1 youngs_modulus = 240e9 poissons_ratio = 0.3 [../] [./f1_function] type = DerivativeParsedMaterial block = 0 args = 'c function = (1-c)^2 [../] [./elasticity_tensor] type = CompositeElasticityTensor block = 0 tensors = 'C1 C2' weights = 'f1 f2' args = 'c' [../]
14 Compute Strain The strain is calculated from the gradient of the displacement vector for a given strain formulation Users can define strains by inheriting from ComputeStrainBase ComputeSmallStrain (Cartesian coordinates) Assumes small strains Computes the total strain ComputeAxisymmetricRZSmallStrain ComputeFiniteStrain (Cartesian coordinates) Strains can be large Uses incremental form ComputeAxisymmetricRZFiniteStrain [./small_strain] type = ComputeSmallStrain block = 1 # displacements in GlobalParams [../] [./finite_strain] type = ComputeFiniteStrain block = 1 # displacements in GlobalParams [../]
15 Compute Stress The stress is calculated from the strain or strain increment according to a specific constitutive model Users can create new constitutive laws by inheriting from ComputeStressBase and overriding computeqpstress Linear elasticity directly calculates the stress from the strain Nonlinear models require a nested Newton solve to determine the stress [./stress] type = ComputeLinearElasticStress block = 1 [../]
16 Linear Elasticity Calculates the stress with and Uses ComputeSmallStrain and ComputeLinearElasticStress void ComputeSmallStrain::computeProperties() { for (_qp = 0; _qp < _qrule->n_points(); ++_qp) { RankTwoTensor grad_tensor(_grad_disp_x[_qp], _grad_disp_y[_qp], _grad_disp_z[_qp]); _total_strain[_qp] = (grad_tensor + grad_tensor.transpose()) / 2.0; } } //Taken from ComputeSmallStrain void ComputeLinearElasticStress::computeQpStress() { // stress = C*e _stress[_qp] = _elasticity_tensor[_qp] * _total_strain[_qp]; } //Taken from ComputeLinearElasticStress
17 Boundary Conditions MOOSE has all of the usual options for boundary conditions. The main boundary conditions for mechanics are: PresetBC Applies a specific constant displacement on a given boundary or node by setting the value of the displacement (Dirichlet). FunctionPresetBC Applies a displacement defined by a function. May vary spatially or with time. Ramped loading Pressure Applies a set pressure on a given boundary in the units of the simulation. NeumannBC Sets the value of the derivative of the solution on the boundary.
18 Linear Elasticity Example To demonstrate linear elasticity, we model a bridge with a body force (gravity) and an applied pressure on the top face See moose/tensor_mechanics/examples/bridge/bridge.i Problem Summary Bridge under gravity and a 0.5 MPa load Material properties: Steel Concrete Boundary conditions: Bottom of supports are constrained in all directions Solved using small strains and linear elasticity
19 Finite Strain Materials For finite strains, the material will deform and also rotate The deformation is defined in terms of the deformation gradient F The deformation gradient is divided into the rotation and the stretch F = RU σ = Rσ R T Stress calculation, σ We use a co-rotational form where the stress is calculated in the intermediate configuration, and then the rotation is applied ComputeFiniteStrain calculates strain_increment rotation_increment strain_rate deformation_gradient
20 Finite Strain Example Problem To demonstrate finite strain, we model a single element that is compressed and then rotated. There are finite strains, and large displacements. See moose/tensor_mechanics/tests/finite_strain_elastic/ elastic_rotation_test.i Finite strain version of the bridge example
21 Finite Strain Plastic Materials FiniteStrainPlasticMaterial: Rate independent von Mises plasticity Flow equations integrated using Backward Euler and Newton Raphson. Required input parameter yield stress which specifies hardening behavior and has to be provided as plastic strain-stress pair. yield_stress = Example: /combined/tests/tensor_mechanics_j2plasticity/ tensor_mechanics_j2plasticity.i FiniteStrainRatePlasticMaterial: Rate dependent von Mises plasticity model. Perzyna model used to define the viscoplastic strain rate
22 Verification with ABAQUS Rate independent von Mises plasticity with large deformation Δu z = 0.01 mm used in both the simulations Equivalent Plastic Strain u z =0 u z =1 mm
23 Phenomenological vs. Mechanistic Plasticity J2 plasticity is phenomenological and based on empirical observation. A more mechanistic view considers dislocation motion In situ TEM video of deformation experiment Video from
24 Crystal Plasticity Crystal plasticity mechanistically represents plastic deformation by representing dislocation slip along slip planes Slip contributions are summed across all slip planes Dislocation generation and pile-up accounted for by increasing slip resistance with a hardening law Crystal reorientation obtained from kinematics Plastic deformation due to slip along slip planes L c p = MX s=1 sb s 0 n s 0 Dislocation motion defined by a flow rule s = s ( s, s 0 )
25 Crystal Plasticity Materials FiniteStrainCrystalPlasticity: Rate dependent plasticity model driven by slip on every slip system Backward Euler time integration using Newton-Raphson. Required parameters: nss number of slip systems; varies for different crystal systems (fcc, bcc, hcp). gprops initial values of slip system resistances. hprops hardening properties flowprops flow properties Presently Euler angle for each element read from an input file Flow rule and hardness evolution can be modified by overriding functions get_slip_incr() and update_gss().
26 Example We have conducted various investigations with crystal plasticity solved with moose, see Chockalingam et al, Comp. Mech, (2013) 51: The results compared very well with those from an Abaqus Umat Compression of a polycrystalline cylinder: Stress (MPa) Huang Abaqus Umat MOOSE JFNK Routine Strain 500 Stress Von Mises (MPa) Grain 1 Grain 2 Grain 3 Grain 4 Grain 5 Grain 6 Grain 7 Grain zz
27 Stress-Free Strains Local stress-free strains can be added to a model, i.e. = C ( 0 ) [Materials] This can be used to add volumetric strains (such as swelling), or strains due to lattice mismatch, etc. Stress-free strains are created in the material block. The base class is ComputeStressFreeStrainBase Constant Stress-free strains are build using ComputeEigenstrain Stress-free strains that are functions of other variables are built using ComputeVariableEigenstrain [./var_dependence] type = DerivativeParsedMaterial block = 0 function = 0.5*c^2 args = c f_name = var_dep enable_jit = true [../] [./eigen_strain] type = ComputeVariableEigenstrain block = 0 eigen_base = ' args = c [../] []
28 Stress-Free Strain Example Secondary phase inclusion strains surrounding material See moose/modules/combined/tests/eigenstrain/variable.i
29 Visualizing the Stress and Strain Tensor material properties can not be directly outputted To output components of the stress, strain, or elastic constants, you must create an auxvariable and a corresponding auxkernel. Tensor auxkernels RankTwoAux: Outputs one component of a rank two tensor RankTwoScalarAux: Outputs scalar norms, Von Mises stress, equivalent plastic strain, trace, L2 norm RankFourAux: Outputs one component of a rank four tensor
30 Contact Module Contact algorithms in a solid mechanics finite element code prevent the penetration of one domain into another or the penetration of a domain into itself. It represents a critical capability to accurately model physical behavior It is also very difficult to accurately model physical contact using finite elements Ω (1) t (1) t (2) Ω (2)
31 Some Required Capabilities for Contact Search Exterior identification Nearby nodes Capture box Binarysearch,e.g. Contact existence More geometric work Penetration point Enforcement Formulation of contact force Formulation of Jacobian for nonlinear solve Interaction with other capabiities (e.g. kinematic boundary conditions) Development, testing, and application testing of contact require many, many man-months of effort. In fact, it is probably man years. A robust (though improving) contact capability is is available in the Contact Module
32 Contact Overview Contact is modeled by tracking interaction of a node on a face Nodes (green) may not penetrate faces (defined by orange nodes) No force is applied where the bodies are not in contact The contact forces must increase from zero as the bodies first come into contact
33 Contact Constraints g 0; the gap (penetration distance) must be non-positive t N 0; the contact force must push bodies apart t N g = 0; the contact force must be zero if the bodies are not in contact t N g = 0; the contact force must be zero when constraints are formed and released The gap in the normal direction for constraint i is (n is the normal, N denotes normal direction, d s is position of the slave node, d c is position of the contact point, and G is a matrix) g g i N = n i (d i (t) g i N = G i N (d(t)) t N d i c(t))
34 Contact Options Formulation:kinematic or penalty Kinematic is more accurate but also harder to solve. Model:frictionless, glued, or coulomb Frictionless enforces the normal constraint and allows nodes to come out of contact if they are in tension. Glued ties nodes where they come into contact with no release. Coulomb is frictional contact with release. friction coefficient Coulomb friction coefficient. penalty The penalty stiffness to be used in the constraint. Master The surface corresponding to the faces in the constraint. Slave The surface corresponding to the nodes in the constraint. normal smoothing distance Distance from face edge in parametric coordinates over which to smooth the normal. Helps with convergence. Try 0.1. tension release The tension value which will allow nodes to be released. Defaults to zero.
35 Contact Example Model of a hard spherical diamond indenter in contact with Cu using finite strain plasticity and frictionless contact. J2 plasticity with isotropic elasticity tensor J2 plasticity with anisotropic elasticity tensor Single crystal plasticity
36 Other Available Mechanics Capability There are users of the mechanics module across the world and they are all adding capability There are various available models Multi-yield surface plasticity Cosserat (Micropolar) Mechanics Mohr-Coulomb plasticity Weak plane tensile and shear behavior (for weakly bonded layers) Other capabilities will come soon Extended Finite Element
Non-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 informationSensitivity and Reliability Analysis of Nonlinear Frame Structures
Sensitivity and Reliability Analysis of Nonlinear Frame Structures Michael H. Scott Associate Professor School of Civil and Construction Engineering Applied Mathematics and Computation Seminar April 8,
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More 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 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 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 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 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 informationLecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity
Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationMODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD
MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used
More informationMathematical Background
CHAPTER ONE Mathematical Background This book assumes a background in the fundamentals of solid mechanics and the mechanical behavior of materials, including elasticity, plasticity, and friction. A previous
More informationANSYS Mechanical Basic Structural Nonlinearities
Lecture 4 Rate Independent Plasticity ANSYS Mechanical Basic Structural Nonlinearities 1 Chapter Overview The following will be covered in this Chapter: A. Background Elasticity/Plasticity B. Yield Criteria
More informationModule-4. Mechanical Properties of Metals
Module-4 Mechanical Properties of Metals Contents ) Elastic deformation and Plastic deformation ) Interpretation of tensile stress-strain curves 3) Yielding under multi-axial stress, Yield criteria, Macroscopic
More informationMECHANICS OF MATERIALS. EQUATIONS AND THEOREMS
1 MECHANICS OF MATERIALS. EQUATIONS AND THEOREMS Version 2011-01-14 Stress tensor Definition of traction vector (1) Cauchy theorem (2) Equilibrium (3) Invariants (4) (5) (6) or, written in terms of principal
More information3-D Finite Element Analysis of Instrumented Indentation of Transversely Isotropic Materials
3-D Finite Element Analysis of Instrumented Indentation of Transversely Isotropic Materials Abstract: Talapady S. Bhat and T. A. Venkatesh Department of Material Science and Engineering Stony Brook University,
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 informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationConcept Question Comment on the general features of the stress-strain response under this loading condition for both types of materials
Module 5 Material failure Learning Objectives review the basic characteristics of the uni-axial stress-strain curves of ductile and brittle materials understand the need to develop failure criteria for
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 informationNonlinear Finite Element Modeling of Nano- Indentation Group Members: Shuaifang Zhang, Kangning Su. ME 563: Nonlinear Finite Element Analysis.
ME 563: Nonlinear Finite Element Analysis Spring 2016 Nonlinear Finite Element Modeling of Nano- Indentation Group Members: Shuaifang Zhang, Kangning Su Department of Mechanical and Nuclear Engineering,
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 informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationUniversity of Sheffield The development of finite elements for 3D structural analysis in fire
The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures
More information3.22 Mechanical Properties of Materials Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 3.22 Mechanical Properties of Materials Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Quiz #1 Example
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 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 informationChapter 5. Vibration Analysis. Workbench - Mechanical Introduction ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.
Workbench - Mechanical Introduction 12.0 Chapter 5 Vibration Analysis 5-1 Chapter Overview In this chapter, performing free vibration analyses in Simulation will be covered. In Simulation, performing a
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 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 informationPrediction of the bilinear stress-strain curve of engineering material by nanoindentation test
Prediction of the bilinear stress-strain curve of engineering material by nanoindentation test T.S. Yang, T.H. Fang, C.T. Kawn, G.L. Ke, S.Y. Chang Institute of Mechanical & Electro-Mechanical Engineering,
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 informationStrength of Material. Shear Strain. Dr. Attaullah Shah
Strength of Material Shear Strain Dr. Attaullah Shah Shear Strain TRIAXIAL DEFORMATION Poisson's Ratio Relationship Between E, G, and ν BIAXIAL DEFORMATION Bulk Modulus of Elasticity or Modulus of Volume
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 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 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 informationSeismic Response Analysis of Structure Supported by Piles Subjected to Very Large Earthquake Based on 3D-FEM
Seismic Response Analysis of Structure Supported by Piles Subjected to Very Large Earthquake Based on 3D-FEM *Hisatoshi Kashiwa 1) and Yuji Miyamoto 2) 1), 2) Dept. of Architectural Engineering Division
More informationChapter 6: Mechanical Properties of Metals. Dr. Feras Fraige
Chapter 6: Mechanical Properties of Metals Dr. Feras Fraige Stress and Strain Tension Compression Shear Torsion Elastic deformation Plastic Deformation Yield Strength Tensile Strength Ductility Toughness
More informationPerformance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis
Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan
More informationSTANDARD SAMPLE. Reduced section " Diameter. Diameter. 2" Gauge length. Radius
MATERIAL PROPERTIES TENSILE MEASUREMENT F l l 0 A 0 F STANDARD SAMPLE Reduced section 2 " 1 4 0.505" Diameter 3 4 " Diameter 2" Gauge length 3 8 " Radius TYPICAL APPARATUS Load cell Extensometer Specimen
More informationChapter 6: Plastic Theory
OHP Mechanical Properties of Materials Chapter 6: Plastic Theory Prof. Wenjea J. Tseng 曾文甲 Department of Materials Engineering National Chung Hsing University wenjea@dragon.nchu.edu.tw Reference: W. F.
More informationOutline. Tensile-Test Specimen and Machine. Stress-Strain Curve. Review of Mechanical Properties. Mechanical Behaviour
Tensile-Test Specimen and Machine Review of Mechanical Properties Outline Tensile test True stress - true strain (flow curve) mechanical properties: - Resilience - Ductility - Toughness - Hardness A standard
More informationA Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against a Rigid Flat under Varying Modulus of Elasticity and Sphere Radius
Engineering, 2010, 2, 205-211 doi:10.4236/eng.2010.24030 Published Online April 2010 (http://www. SciRP.org/journal/eng) 205 A Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against
More informationExample-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium
Example-3 Title Cylindrical Hole in an Infinite Mohr-Coulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic
More informationEXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD
Journal of KONES Powertrain and Transport, Vol. 7, No. EXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD Robert Czabanowski Wroclaw University
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 2 PLASTICITY
MATERIAL MECHANICS, SE2126 COMPUTER LAB 2 PLASTICITY PART A INTEGRATED CIRCUIT An integrated circuit can be thought of as a very complex maze of electronic components and metallic connectors. These connectors
More informationLecture contents. Stress and strain Deformation potential. NNSE 618 Lecture #23
1 Lecture contents Stress and strain Deformation potential Few concepts from linear elasticity theory : Stress and Strain 6 independent components 2 Stress = force/area ( 3x3 symmetric tensor! ) ij ji
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 informationBulk Metal Forming II
Bulk Metal Forming II Simulation Techniques in Manufacturing Technology Lecture 2 Laboratory for Machine Tools and Production Engineering Chair of Manufacturing Technology Prof. Dr.-Ing. Dr.-Ing. E.h.
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationBrittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm
Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd
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 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 informationA FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR AGAINST A RIGID FLAT UNDER VARYING MODULUS OF ELASTICITY AND SPHERE RADIUS
Proceedings of the International Conference on Mechanical Engineering 2009 (ICME2009) 26-28 December 2009, Dhaka, Bangladesh ICME09- A FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR
More informationThe Finite Element Method II
[ 1 The Finite Element Method II Non-Linear finite element Use of Constitutive Relations Xinghong LIU Phd student 02.11.2007 [ 2 Finite element equilibrium equations: kinematic variables Displacement Strain-displacement
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 informationLinear Cosserat elasticity, conformal curvature and bounded stiffness
1 Linear Cosserat elasticity, conformal curvature and bounded stiffness Patrizio Neff, Jena Jeong Chair of Nonlinear Analysis & Modelling, Uni Dui.-Essen Ecole Speciale des Travaux Publics, Cachan, Paris
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
More information(MPa) compute (a) The traction vector acting on an internal material plane with normal n ( e1 e
EN10: Continuum Mechanics Homework : Kinetics Due 1:00 noon Friday February 4th School of Engineering Brown University 1. For the Cauchy stress tensor with components 100 5 50 0 00 (MPa) compute (a) The
More informationIntroduction to Engineering Materials ENGR2000. Dr. Coates
Introduction to Engineering Materials ENGR2 Chapter 6: Mechanical Properties of Metals Dr. Coates 6.2 Concepts of Stress and Strain tension compression shear torsion Tension Tests The specimen is deformed
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 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 informationPractice Final Examination. Please initial the statement below to show that you have read it
EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use
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 informationDiscrete Element Modeling of Soils as Granular Materials
Discrete Element Modeling of Soils as Granular Materials Matthew R. Kuhn Donald P. Shiley School of Engineering University of Portland National Science Foundation Grant No. NEESR-936408 Outline Discrete
More informationNumerical Modelling of Dynamic Earth Force Transmission to Underground Structures
Numerical Modelling of Dynamic Earth Force Transmission to Underground Structures N. Kodama Waseda Institute for Advanced Study, Waseda University, Japan K. Komiya Chiba Institute of Technology, Japan
More information1. A pure shear deformation is shown. The volume is unchanged. What is the strain tensor.
Elasticity Homework Problems 2014 Section 1. The Strain Tensor. 1. A pure shear deformation is shown. The volume is unchanged. What is the strain tensor. 2. Given a steel bar compressed with a deformation
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 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 informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationChapter 2: Elasticity
OHP 1 Mechanical Properties of Materials Chapter 2: lasticity Prof. Wenjea J. Tseng ( 曾文甲 ) Department of Materials ngineering National Chung Hsing University wenjea@dragon.nchu.edu.tw Reference: W.F.
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 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 informationFinite Element Simulations of Microbeam Bending Experiments
Finite Element Simulations of Microbeam Bending Experiments Master s thesis in Applied Mechanics JOHN WIKSTRÖM Department of Applied Mechanics CHALMERS UNIVERSITY OF TECHNOLOGY Göteborg, Sweden 2017 MASTER
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 Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationMechanics of Biomaterials
Mechanics of Biomaterials Lecture 7 Presented by Andrian Sue AMME498/998 Semester, 206 The University of Sydney Slide Mechanics Models The University of Sydney Slide 2 Last Week Using motion to find forces
More informationConstitutive Relations
Constitutive Relations Dr. Andri Andriyana Centre de Mise en Forme des Matériaux, CEMEF UMR CNRS 7635 École des Mines de Paris, 06904 Sophia Antipolis, France Spring, 2008 Outline Outline 1 Review of field
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 informationContinuum Mechanics and the Finite Element Method
Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after
More informationTheory of Plasticity. Lecture Notes
Theory of Plasticity Lecture Notes Spring 2012 Contents I Theory of Plasticity 1 1 Mechanical Theory of Plasticity 2 1.1 Field Equations for A Mechanical Theory.................... 2 1.1.1 Strain-displacement
More informationMODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Uniaxial Model: Strain-Driven Format of Elastoplasticity
MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #3: Elastoplastic Concrete Models Uniaxial Model: Strain-Driven Format of Elastoplasticity Triaxial
More informationUNLOADING OF AN ELASTIC-PLASTIC LOADED SPHERICAL CONTACT
2004 AIMETA International Tribology Conference, September 14-17, 2004, Rome, Italy UNLOADING OF AN ELASTIC-PLASTIC LOADED SPHERICAL CONTACT Yuri KLIGERMAN( ), Yuri Kadin( ), Izhak ETSION( ) Faculty of
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 informationUsing the Timoshenko Beam Bond Model: Example Problem
Using the Timoshenko Beam Bond Model: Example Problem Authors: Nick J. BROWN John P. MORRISSEY Jin Y. OOI School of Engineering, University of Edinburgh Jian-Fei CHEN School of Planning, Architecture and
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 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 information2D Liquefaction Analysis for Bridge Abutment
D Liquefaction Analysis for Bridge Abutment Tutorial by Angel Francisco Martinez Integrated Solver Optimized for the next generation 64-bit platform Finite Element Solutions for Geotechnical Engineering
More informationMechanical analysis of timber connection using 3D finite element model
Mechanical analysis of timber connection using 3D finite element model Bohan XU Ph.D Student Civil Engineering Laboratory (CUST) Clermont-Ferrand, France Mustapha TAAZOUNT Dr-Ing Civil Engineering Laboratory
More informationEQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS
EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS Abdul Qadir Bhatti 1, Norimitsu Kishi 2 and Khaliq U Rehman Shad 3 1 Assistant Professor, Dept. of Structural Engineering,
More informationBending Load & Calibration Module
Bending Load & Calibration Module Objectives After completing this module, students shall be able to: 1) Conduct laboratory work to validate beam bending stress equations. 2) Develop an understanding of
More informationSupplementary Figures
Fracture Strength (GPa) Supplementary Figures a b 10 R=0.88 mm 1 0.1 Gordon et al Zhu et al Tang et al im et al 5 7 6 4 This work 5 50 500 Si Nanowire Diameter (nm) Supplementary Figure 1: (a) TEM image
More informationLarge Thermal Deflections of a Simple Supported Beam with Temperature-Dependent Physical Properties
Large Thermal Deflections of a Simple Supported Beam with Temperature-Dependent Physical Properties DR. ŞEREF DOĞUŞCAN AKBAŞ Civil Engineer, Şehit Muhtar Mah. Öğüt Sok. No:2/37, 34435 Beyoğlu- Istanbul,
More informationReference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity",
Reference material Reference books: Y.C. Fung, "Foundations of Solid Mechanics", Prentice Hall R. Hill, "The mathematical theory of plasticity", Oxford University Press, Oxford. J. Lubliner, "Plasticity
More informationThe University of Melbourne Engineering Mechanics
The University of Melbourne 436-291 Engineering Mechanics Tutorial Four Poisson s Ratio and Axial Loading Part A (Introductory) 1. (Problem 9-22 from Hibbeler - Statics and Mechanics of Materials) A short
More informationINCREASING RUPTURE PREDICTABILITY FOR ALUMINUM
1 INCREASING RUPTURE PREDICTABILITY FOR ALUMINUM Influence of anisotropy Daniel Riemensperger, Adam Opel AG Paul Du Bois, PDB 2 www.opel.com CONTENT Introduction/motivation Isotropic & anisotropic material
More informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
More informationPlates and Shells: Theory and Computation. Dr. Mostafa Ranjbar
Plates and Shells: Theory and Computation Dr. Mostafa Ranjbar Outline -1-! This part of the module consists of seven lectures and will focus on finite elements for beams, plates and shells. More specifically,
More information