Orthotropic Materials

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

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

KINEMATICS OF RIGID BODIES

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

NUMERICAL SIMULATION FOR NONLINEAR STATIC & DYNAMIC STRUCTURAL ANALYSIS

MECHANICS OF MATERIALS Poisson s Ratio

, on the power of the transmitter P t fed to it, and on the distance R between the antenna and the observation point as. r r t

Sections 3.1 and 3.4 Exponential Functions (Growth and Decay)

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

Projection of geometric models

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

Computer Propagation Analysis Tools

Projection of geometric models

Pressure Vessels Thin and Thick-Walled Stress Analysis

Today - Lecture 13. Today s lecture continue with rotations, torque, Note that chapters 11, 12, 13 all involve rotations

EN221 - Fall HW # 7 Solutions

7 Wave Equation in Higher Dimensions

The Production of Polarization

PHYS PRACTICE EXAM 2

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

Monochromatic Wave over One and Two Bars

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

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

Control Volume Derivation

Lecture 22 Electromagnetic Waves

Integration of the constitutive equation

r r r r r EE334 Electromagnetic Theory I Todd Kaiser

An Automatic Door Sensor Using Image Processing

The k-filtering Applied to Wave Electric and Magnetic Field Measurements from Cluster

P h y s i c s F a c t s h e e t

Lecture 5. Chapter 3. Electromagnetic Theory, Photons, and Light

Heat Conduction Problem in a Thick Circular Plate and its Thermal Stresses due to Ramp Type Heating

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

BEng (Hons) Telecommunications. Examinations for / Semester 2

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

Chapter Finite Difference Method for Ordinary Differential Equations

CAMERA GEOMETRY. Ajit Rajwade, CS 763, Winter 2017, IITB, CSE department

ME 3560 Fluid Mechanics

ANALYTICAL SOLUTION FOR EDDY CURRENT PROBLEM, USING SPACE EIGENFUNCTIONS EXPANSION

Ö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

A Strain Based Design Criterion for Solid Propellant Rocket Motors

Q & Particle-Gas Multiphase Flow. Particle-Gas Interaction. Particle-Particle Interaction. Two-way coupling fluid particle. Mass. Momentum.

envionmen ha implemens all of he common algoihmic deails of all nodal mehods u pemis he specic mehod o e used in any concee insance o e specied y he u

ME 304 FLUID MECHANICS II

Overview. Overview Page 1 of 8

Exponential and Logarithmic Equations and Properties of Logarithms. Properties. Properties. log. Exponential. Logarithmic.

FINITE DIFFERENCE APPROACH TO WAVE GUIDE MODES COMPUTATION

Damage Assessment in Composites using Fiber Bragg Grating Sensors. Mohanraj Prabhugoud

Research Article Stress Analysis of Nonhomogeneous Rotating Disc with Arbitrarily Variable Thickness Using Finite Element Method

Energy dispersion relation for negative refraction (NR) materials

Online Completion of Ill-conditioned Low-Rank Matrices

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

EVENT HORIZONS IN COSMOLOGY

Fluid Flow and Heat Transfer Characteristics across an Internally Heated Finned Duct

KINGS UNIT- I LAPLACE TRANSFORMS

Relative and Circular Motion

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

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

Circular Motion. Radians. One revolution is equivalent to which is also equivalent to 2π radians. Therefore we can.

Measures the linear dependence or the correlation between r t and r t-p. (summarizes serial dependence)

Low-complexity Algorithms for MIMO Multiplexing Systems

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

PRINCIPLES OF USING GROUND PENETRATING RADAR AND GPS TECHNOLOGY FOR DETECTION OF UNDERGROUND UTILITIES

The Method of Images in Velocity-Dependent Systems

Physics Tutorial V1 2D Vectors

AN EFFICIENT INTEGRAL METHOD FOR THE COMPUTATION OF THE BODIES MOTION IN ELECTROMAGNETIC FIELD

Lecture 10: Wave equation, solution by spherical means

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,

Physics 2001/2051 Moments of Inertia Experiment 1

Dual Hierarchies of a Multi-Component Camassa Holm System

Range Migration Techniques for Short-Range MIMO Array Imaging

Complex Eigenvalues. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

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

BUCKLING AND ULTIMATE STRENGTH ASSESSMENT FOR OFFSHORE STRUCTURES

Chapter 2: The Derivation of Maxwell Equations and the form of the boundary value problem

Risk tolerance and optimal portfolio choice

Camera Models class 8

K. G. Malyutin, T. I. Malyutina, I. I. Kozlova ON SUBHARMONIC FUNCTIONS IN THE HALF-PLANE OF INFINITE ORDER WITH RADIALLY DISTRIBUTED MEASURE

2-d Motion: Constant Acceleration

On The Estimation of Two Missing Values in Randomized Complete Block Designs

ON 3-DIMENSIONAL CONTACT METRIC MANIFOLDS

New method to explain and calculate the gyroscopic torque and its possible relation to the spin of electron

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

( ) exp i ω b ( ) [ III-1 ] exp( i ω ab. exp( i ω ba

2. v = 3 4 c. 3. v = 4c. 5. v = 2 3 c. 6. v = 9. v = 4 3 c

On the Semi-Discrete Davey-Stewartson System with Self-Consistent Sources

Wavefront healing operators for improving reflection coherence

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

2D vector fields 1. Contents

Engineering Accreditation. Heat Transfer Basics. Assessment Results II. Assessment Results. Review Definitions. Outline

Fig. 1S. The antenna construction: (a) main geometrical parameters, (b) the wire support pillar and (c) the console link between wire and coaxial

LECTURE 5. is defined by the position vectors r, 1. and. The displacement vector (from P 1 to P 2 ) is defined through r and 1.

Department of Chemical Engineering University of Tennessee Prof. David Keffer. Course Lecture Notes SIXTEEN

MEMS 0031 Electric Circuits

Patent Examination Duration in an Endogenous Growth Model

Determination of Stresses in Drying Wood by Means of a Viscoelastic Relaxation Model

Research on the Algorithm of Evaluating and Analyzing Stationary Operational Availability Based on Mission Requirement

Modelling Hydromechanical Dilation Geomaterial Cavitation and Localization

LawsoftheElectroElectricalInduction

AST1100 Lecture Notes

Transcription:

Kapiel 2 Ohoopic Maeials 2. Elasic Sain maix Elasic sains ae elaed o sesses by Hooke's law, as saed below. The sesssain elaionship is in each maeial poin fomulaed in he local caesian coodinae sysem. ε e = Cσ () σ is he sess veco, which uses he same convenion as he sain veco. The sess veco is given below. σ = [σ l σ σ τ l τ l τ ] T (2) C is he maeial compliance maix. The compliance maix is given as C = ν l ν l E ν E ν E G l 0 0 0 G l 0 0 0 G ν l ν l (3) E's ae moduli of elasiciy, G's ae shea moduli and ν's ae Poissons aios. The compliance maix should be symmeic, which gives ise o he following esic- 4

ions: ν l = ν l ; ν l = ν l E ; ν = ν E (4) The invese of he compliance maix is he maeial siness maix, D, which give he sesses poduced by an elasic sain sae. D = ( ν ν )/k (ν ν l + ν l )/k (ν l ν + ν l )/k (ν l ν + ν l )/k ( ν l ν l )/k (ν l ν l + ν )/k E (ν l ν + ν l )/k E (ν l ν l + ν )/k E ( ν l ν l )/k G l 0 0 0 G l 0 0 0 G (5) whee k = ν ν ν l ν l ν l ν l ν l ν ν l ν l ν ν l. The siness maix is also symmeic, which follows fom he esicions given in (4). Fuhemoe D should be posiive deni i.e. ε T Dε > 0 (The maeial esiss defomaions), which leads o he following esicions: ν ν ν l ν l ν l ν l ν l ν ν l ν l ν ν l > 0 ν ν > 0 ν l ν l > 0 ν l ν l > 0 (6) The maeial siness maix is used when he equilibium equaions ae solved, using he Finie elemen mehod. 2.2 Coodinae sysems Thee dieen caesian coodinae sysems ae used in he spaial disceizaion, one global and wo dieen local associaed wih each elemen. The local elemen coodinae sysems ae called elemen coodinae sysem and maeial coodinae sysem, especively. The hee dieen coodinae sysems ae shown in gue 2.2. Global coodinae sysem The global coodinaes ae emed (X, Y, Z). The oigin is locaed in global node numbe. The Global coodinae sysem is used o dene he geomey of he 5

x y x X Z x y y Y z Figu : Local and global coodinae sysems imbe beam, and he bounday condiions. 2.2.2 Elemen coodinae sysem The s of he local coodinae sysems is called he elemen coodinae sysem. The coodinae axes ae emed (x, y, z). The oigin is locaed in he cene of he elemen. The posiive diecion fo he x-axis is fom node 2 owads node. The posiive diecion fo he y-axis is fom node 4 owads node and nally he posiive diecion fo he z-axis is fom node 3 owads node. The elemen coodinae sysem is used fo he inepolaion beween he elemen nodes and he elemen inegaion. 2.2.3 Maeial coodinae sysem The second coodinae sysem is he elemen maeial coodinae sysem, wih axes (l,, ). The oigin of he maeial coodinae sysem is also locaed in he cene of he elemen, and he maeial diecions ae assumed o be consan in he elemen, which make i necessay wih a lage numbe of elemens in plane pependicula o he gain diecion. 6

2.2.4 Tansfomaion beween coodinae sysems The ansfomaion of popeies descibed in he maeial coodinae sysem o he elemen coodinae sysem is descibed in his secion. Popeies can be geomeic poins, maeial popeies eg. maeial siness o physical quaniies eg. sess o hea ux. The ansfomaion of geomeical poins beween he wo coodinae sysems ae: x y z = A T l (7) Whee A is given by (8) A = l x l y l z x y z x y z (8) The s ow in A is a uni veco in he l-diecion, descibed in he xyz-coodinae sysem, he second is a uni veco in he -diecion and he las ow is a uni veco in he -diecion. The elemens in A ae denoed diecion cosines. The deeminaion of diecion cosines is given i secion 2.2.5. A is ofen called a oaion maix, and is an ohogonal maix; i.e. A = A T. Hence he invese of he ansfomaion in (7) is given below l = A x y z (9) The ansfomaion of sesses and sains, fom one coodinae sysem o anohe is as given in (0) ε l = Gε xyz σ xyz = G T σ l (0) σ l = ( G T) σxyz whee σ l = [σ l σ σ τ l τ l τ ] T is he sess veco, descibed in he maeial coodinae sysem and σ xyz is he sess veco descibed in he elemen coodinae 7

sysem. The ansfomaion maix G can be deduced fom he elemens in (8), and is as given hee: G = lx 2 ly 2 lz 2 l x l y l x l z l y l z x 2 y 2 z 2 x y x z y z 2 x 2 y 2 z x y x z y z 2l x x 2l y y 2l z z l x y + l y x l z x + l x z l y z + l z y 2l x x 2l y y 2l z z x l y + y l x z l x + x l z y l z + z l y 2 x x 2 y y 2 z z x y + y x z x + x z y z + z y () The ( G T) maix is given in (2), when he G maix is subdivided ino fou 3 3 maices, emed G, G 2, G 2 and G 22. G = [ ] G G 2 G 2 G 22 [, G G 2G = 2 G 2 2 G 22 ] T (2) The ansfomaion of siness o exibiliy popeies, fom he maeial o he elemen coodinae sysem, is pefomed by a enso-like ansfomaion, as saed below. D xyz = G T D l G (3) whee D l is he maeial siness maix, fomulaed in he maeial coodinae sysem, given by (5) and D xyz is he maeial siness maix fomulaed in he elemen coodinae sysem. The G maix is he ansfomaion maix given in () The ansfomaion of hemal conduciviies and moisue anspo popeies is done in a simila enso-like way, as saed below λ xyz = A T λ l A (4) whee λ l is he maeial conduciviy maix fomulaed in he maeial coodinae sysem, which is a 3 by 3 diagonal maix conaining he hemal conduciviies in he pincipal diecions. λ xyz is he maeial conduciviy maix fomulaed in he elemen coodinae sysem. A is he ansfomaion maix given in (8). 2.2.5 Diecion Cosines The diecion cosines ae deemined fom he dieen elemens locaion elaive o he pih, see gue, which ae emed A 0. A 0 has he same sucue as he maix given in (8). Fuhe he diecions cosines ae depended on wo gowh phenomena, being conical angle and spial gowh, as given below. A = A 0 A c A s (5) 8

Whee A c and A s ae given by (6). cos φ sin φ 0 A c = sin φ cos φ 0 A s = 0 0 cos ϖ 0 sin ϖ 0 0 sin ϖ 0 cos ϖ (6) φ is he conical angle and ϖ is he spial gain angle. Boh angles ae illusaed in gue 2 φ Figu 2: conical angle and spial gowh 9