Integration of the constitutive equation

Similar documents
Orthotropic Materials

Lecture 17: Kinetics of Phase Growth in a Two-component System:

Chapter Finite Difference Method for Ordinary Differential Equations

Control Volume Derivation

Lecture 18: Kinetics of Phase Growth in a Two-component System: general kinetics analysis based on the dilute-solution approximation

MEEN 617 Handout #11 MODAL ANALYSIS OF MDOF Systems with VISCOUS DAMPING

Sharif University of Technology - CEDRA By: Professor Ali Meghdari

7 Wave Equation in Higher Dimensions

Combinatorial Approach to M/M/1 Queues. Using Hypergeometric Functions

On Control Problem Described by Infinite System of First-Order Differential Equations

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 9 Solutions [Theorems of Gauss and Stokes]

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

Coupled Mass Transport and Reaction in LPCVD Reactors

ME 3560 Fluid Mechanics

An Open cycle and Closed cycle Gas Turbine Engines. Methods to improve the performance of simple gas turbine plants

Special Subject SC434L Digital Video Coding and Compression. ASSIGNMENT 1-Solutions Due Date: Friday 30 August 2002

ψ(t) = V x (0)V x (t)

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

The sudden release of a large amount of energy E into a background fluid of density

The Fundamental Theorems of Calculus

Polymer Engineering (MM3POE)

, u denotes uxt (,) and u. mean first partial derivatives of u with respect to x and t, respectively. Equation (1.1) can be simply written as

Structural Dynamics and Earthquake Engineering

DESIGN OF TENSION MEMBERS

Scalar Conservation Laws

International Journal of Mathematical Archive-5(6), 2014, Available online through ISSN

[ ] 0. = (2) = a q dimensional vector of observable instrumental variables that are in the information set m constituents of u

Integral Control via Bias Estimation

Two-dimensional Effects on the CSR Interaction Forces for an Energy-Chirped Bunch. Rui Li, J. Bisognano, R. Legg, and R. Bosch

Elastic-Plastic Deformation of a Rotating Solid Disk of Exponentially Varying Thickness and Exponentially Varying Density

2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum

Dispersive Systems. 1) Schrödinger equation 2) Cubic Schrödinger 3) KdV 4) Discreterised hyperbolic equation 5) Discrete systems.

A note on characterization related to distributional properties of random translation, contraction and dilation of generalized order statistics

G. =, etc.

Solutionbank Edexcel AS and A Level Modular Mathematics

where the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).

Homework 2 Solutions

Propagation of Torsional Surface Waves. in Heterogeneous Half-Space. with Irregular Free Surface

The Global Trade and Environment Model: GTEM

Lecture-V Stochastic Processes and the Basic Term-Structure Equation 1 Stochastic Processes Any variable whose value changes over time in an uncertain

STUDY OF THE STRESS-STRENGTH RELIABILITY AMONG THE PARAMETERS OF GENERALIZED INVERSE WEIBULL DISTRIBUTION

1 Fundamental Solutions to the Wave Equation

Critical statics and dynamics of QCD CEP

ÖRNEK 1: THE LINEAR IMPULSE-MOMENTUM RELATION Calculate the linear momentum of a particle of mass m=10 kg which has a. kg m s

Notes on Optimal Control

Time-Space Model of Business Fluctuations

Recent Developments in the Unit Root Problem for Moving Averages

D.I. Survival models and copulas

PHYS PRACTICE EXAM 2

Chapter Three Systems of Linear Differential Equations

Quantum Mechanics I - Session 5

Moment. F r F r d. Magnitude of moment depends on magnitude of F and the length d

THE EFFECT OF SUCTION AND INJECTION ON UNSTEADY COUETTE FLOW WITH VARIABLE PROPERTIES

MOMENTUM CONSERVATION LAW


Central Force Motion

MATHEMATICAL FOUNDATIONS FOR APPROXIMATING PARTICLE BEHAVIOUR AT RADIUS OF THE PLANCK LENGTH

MA Study Guide #1

Millennium Theory Equations Original Copyright 2002 Joseph A. Rybczyk Updated Copyright 2003 Joseph A. Rybczyk Updated March 16, 2006

NUMERICAL SIMULATION FOR NONLINEAR STATIC & DYNAMIC STRUCTURAL ANALYSIS

Research Article A Note on Multiplication and Composition Operators in Lorentz Spaces

Output Tracking Combing Output Redefinition and Non-causal Stable Inversion for Non-minimum Systems

N H. be the number of living fish outside area H, and let C be the cumulative catch of fish. The behavior of N H

r P + '% 2 r v(r) End pressures P 1 (high) and P 2 (low) P 1 , which must be independent of z, so # dz dz = P 2 " P 1 = " #P L L,

Settling Time Design and Parameter Tuning Methods for Finite-Time P-PI Control

Variance and Covariance Processes

Molecular Evolution and Phylogeny. Based on: Durbin et al Chapter 8

General Non-Arbitrage Model. I. Partial Differential Equation for Pricing A. Traded Underlying Security

Quantum Mechanics II

156 There are 9 books stacked on a shelf. The thickness of each book is either 1 inch or 2

Analytic Model and Bilateral Approximation for Clocked Comparator

Unsteady laminar flow of visco-elastic fluid of second order type between two parallel plates

8.5 Circles and Lengths of Segments

Pseudosteady-State Flow Relations for a Radial System from Department of Petroleum Engineering Course Notes (1997)

Servomechanism Design

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

@FMI c Kyung Moon Sa Co.

Low-complexity Algorithms for MIMO Multiplexing Systems

Voltage/current relationship Stored Energy. RL / RC circuits Steady State / Transient response Natural / Step response

EN221 - Fall HW # 7 Solutions

As is natural, our Aerospace Structures will be described in a Euclidean three-dimensional space R 3.

15.081J/6.251J Introduction to Mathematical Programming. Lecture 6: The Simplex Method II

General Relativity Homework 5

A Crash Course in (2 2) Matrices

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Quiz 6--Work, Gravitation, Circular Motion, Torque. (60 pts available, 50 points possible)

1.2 Differential cross section

BENDING OF BEAM. Compressed layer. Elongated. layer. Un-strained. layer. NA= Neutral Axis. Compression. Unchanged. Elongation. Two Dimensional View

Chapter 7 Response of First-order RL and RC Circuits

Dual Hierarchies of a Multi-Component Camassa Holm System

The shortest path between two truths in the real domain passes through the complex domain. J. Hadamard

Static equilibrium requires a balance of forces and a balance of moments.

Section 3.8, Mechanical and Electrical Vibrations

Pressure Vessels Thin and Thick-Walled Stress Analysis

A Mathematical model to Solve Reaction Diffusion Equation using Differential Transformation Method

POSITIVE SOLUTIONS WITH SPECIFIC ASYMPTOTIC BEHAVIOR FOR A POLYHARMONIC PROBLEM ON R n. Abdelwaheb Dhifli

ON THE OSCILLATION OF THIRD ORDER FUNCTIONAL DIFFERENTIAL EQUATIONS. Cairo University, Orman, Giza 12221, Egypt

Modelling Hydromechanical Dilation Geomaterial Cavitation and Localization

Deviation probability bounds for fractional martingales and related remarks

Electrical and current self-induction

Transcription:

Inegaion of he consiive eqaion REMAINDER ON NUMERICAL INTEGRATION Analyical inegaion f ( x( ), x ( )) x x x f () Exac/close-fom solion (no always possible) Nmeical inegaion. i. N T i N [, T ] [ i, i ] i [ i, i ] Time ineval/sep "i" Δ i i Time incemen x( ) x g( x, x, x,...) Appoximae solion (epening on i i i i i he vales a pevios iscee imes)

Inegaion of he evolion eqaion Loaing/Unloaing eqaions Time inegaion λ ; ; (, ) g λg λ g λg g τ ( ) f ( ) a) ELASTIC STATE ï ì g ºe - < e () No evolion of ïì < í l í ïî e () < ïî e () ³ b) INELASTIC (DAMAGE) STATE b-) UNLOADING ì ï g ºe - í ïî e () g º () - < l No evolion of ( ) < e b-) LOADING g º () - e b--) NEUTRAL LOADING l No evolion of () b--) PURE LOADING l ( ) > Evolion of () ³ e e e e

Inegaion of he evolion an consiive eqaions SUMMARY : ) he iniial vale of is o > ) he iniial vale of is ( ) 3) neve eceases ( ³ ) e 4) is always smalle o eqal han (g º - ) e e if e gows e ³ 5) When e eqals hen ìï í ïî ï if e eceases e < ; e () max(, τ ( s)) s [,] τ s hisoical máximm of [, ( )] () s q q( ) s τ q ( ) : The inegaion is exac (close fom). No epenence on D, τ

Tangen consiive opeao (efiniion) (,) x (,) x Σ((,)) x Tangen consiive opeao Σ( ) ( ) ang () Σ( ij ) ij( ) [ ang ] ijkl i, j, k, l {,,3} kl kl Tangen consiive opeao Secan consiive opeao

Tangen consiive opeao (analyical eivaion) Tangen consiive moli ((), ) ( ( )) ( ): ( ()) ang () : ( ) ( ) ang () ( ) : ang : ang (): Deivaion ( ) : : : ang Elasic egime/unloaing/neal loaing ( ) : ( ) ( ang ang sec secan consiive opeao)

Tangen consiive opeao (analyical eivaion) Inelasic egime-loaing q( ) ( ) τ () q H H () 3 τ > q '( ) q( ) () '( ) τ : : : : ( : : ) τ : : : : : : : : : : : ; : : : q H () 3 q H : : 3 : : τ

Tangen consiive opeao (analyical eivaion) H q ( ) : : ( ) : ; : : 3 q H q H ( ) : ( : ) : [( ) ( [ : ])]: 3 3 ang : q H ang () ( ) ( [ : ]) 3 Symmeic moel (ension/compession) ang : : : : ( ) (elasic/nloaing) () q H ( ) ( ) (loaing) 3 ang

Inegaion of he consiive eqaion: Ienificaion of he cen sae Ienificaion a he en of he cen ime ineval) ) Elasic gδ τ Δ Δ < τ Δ < Δ Δ Δ Δ ) Unloaing/Neal loaing τ Δ [, Δ] gδ τ Δ Δ τ τ Δ Δ Δ Δ Δ Δ 3) Loaing gδ τ Δ Δ Δ > τ Δ τ > Δ Δ > Δ > Δ Δ > Δ <

Inegaion of he consiive eqaion Nmeical algoihm: INPUT DATA [, Δ],, Sep Sep ( Compe If τ Δ s Δ Δ τ ( ) s Δ q( Δ) Δ Δ Δ Δ ang ) ( ) Δ Δ Δ Δ : Δ : : : s Δ Δ Δ Δ Elasic Unloaing Neal loaing EXIT

Inegaion of he consiive eqaion Nmeical algoihm: Sep 3 Δ Δ If τ > τ Δ Δ q( Δ) Δ ( ) Δ Δ Δ (Loaing) q( ) H ( ) ( ) ( ) Δ Δ Δ ang Δ Δ 3 Δ Δ ( Δ ) EXIT OUTPUT DATA [, Δ],,( ) Δ Δ ang Δ

Chaaceizaion of he Haening/Sofening law Linea haening/sofening H ( ) [, ] ( ( q ) q( ) q( ) H q > [, ) q( ) H [, ) H [, )

Chaaceizaion of he Haening/Sofening law Exponenial haening/sofening A ( ) q q q e A ( ) ( ) ( > ) q( ) H ( ) A e ( q ) A ( )

Chaaceizaion of he Haening/Sofening law Vale of I is compe fom he D case niaxial elasic limi τ : : : E ( ) : ij ij E : E

Chaaceizaion of he elasic omain { f τ q } Elasic omain Ε : (, ) ( ). Symmeic (ension/compession) moel τ q τ : : ( ) τ ( ) : : τ : : 3 E Remak Eqal ension an compession elasic limis.

Chaaceizaion of he elasic omain. Tensile-amage-only moel Posiive conepa of a scala fncion (McAley backe) x x > x x < Posiive conepa of a sess enso In maix foma: [ ] iagonalizaion [ ] iag 3 en o oiginal [ ] iag sysem of [ ] Posiive 3 cooinaes x conepa of [ ] x

Chaaceizaion of he elasic omain In enso foma: i 3 i eigenveco"" i i 3 ˆ ˆ i pi pi i eigenvale"" i pˆ pˆ i i i an have he same eigenvecos shaes he posiive eigenvales of has hose negaive eigenvales of nll τ τ ( ) ( ) Sain nom eefiniion τ : : ( ) : : ( ) : : ( ) τ ( : : ) :

Chaaceizaion of he elasic omain a) If b) If > > 3 > 3 3 Pe ensile sae : : τ τ Pe compessive sae < < τ : : f τ < 3 < 3 The sae is always elasic Inelasic ν ν E ν ν Elasic

Chaaceizaion of he elasic omain Non-symmeic ension-compession moel θ τ θ : : n > θ τ 3 3 θ,, 3 : : i i,, 3 < θ τ : : 3 3 n i i n E n n

. Viscoamage moel Ch.. Coninm Damage Moels

Moel chaaceizaion The ae effecs can be accommoae ino he invisci amage moel (ae inepenen) Evolion eqaion (,),, [, ] λ Consiive eqaion ( ) : Damage fncion g(, ) τ Kash-Khn Tcke an pesisency coniions λ, g λg if g λg

Moel chaaceizaion Visco eglaizaion (Pezyna s Reglaizaion) λ g η Obaine by eplacing in he invisci moel Evolion eqaion λ(,) g η Consiive eqaion whee η is he viscosiy The only change is he evolion law is he McAley backe ( ) : Damage fncion g(,) τ Kash Khn Tcke an pesisency coniions No necessay x x> x x < x x > x x x < x x x

Moel chaaceizaion Remak ; g η λ g η λg g g á gñ ( λη) ηλ η η η η g Remak g λη η As λη g g Remak 3 λ g λg λg g λ g If η λ g As η λ g λg If g λ g Loaing/nloaing coniions Pesisency coniions As η (invisci case ) he ae inepenen coninm amage moel is ecovee

Moel chaaceizaion Remak 3 g τ > τ > Fo η an λ hen g λη f τ q > τ > q The sess/sain sae can lay osie he elasic omain Δ 3 Δ 3 τ q τ 3 Remak 4 f g g η f > g > g > η 3 No evolion ( nloaing / neal loaing) Evolion ( loaing)

Fee enegy an issipaion Fee Enegy ψ (, ) ( ) ψ ( ) ( ) ( : : ) Dissipaion : ψ Coleman's Theoem ψ(, ()) ψ(, ()) : : ψ (, ()) ( ( )) ψ ( ()), ( ( )) : ψ ( ()), : ψ :

Tangen consiive enso Tangen consiive opeao (,) (()) Invisci moel ( as a paamee) ((),) Viscos moel ( as an inepenen vaiable) ( ) : : ( ) : ( ( )) : '( ) ( ( )) : '( ) η g ( ( )) : '( ) g(, ) η (,) : ( ) '( ) g(, ) : η v ang

Time epenency/rae epenency Remak ((),) The sess enso can change even if emains consan Invisci Viscos

END OF LECTURE 4