DYNAMIC PROPAGATION OF A WEAK-DISCONTINUOUS INTERFACE CRACK IN FUNCTIONALLY GRADED LAYERS UNDER ANTI-PLANE SHEAR

Similar documents
4. Eccentric axial loading, cross-section core

TWO INTERFACIAL COLLINEAR GRIFFITH CRACKS IN THERMO-ELASTIC COMPOSITE MEDIA

LAPLACE TRANSFORM SOLUTION OF THE PROBLEM OF TIME-FRACTIONAL HEAT CONDUCTION IN A TWO-LAYERED SLAB

CENTROID (AĞIRLIK MERKEZİ )

Chapter Runge-Kutta 2nd Order Method for Ordinary Differential Equations

ORDINARY DIFFERENTIAL EQUATIONS

13 Design of Revetments, Seawalls and Bulkheads Forces & Earth Pressures

Physics 121 Sample Common Exam 2 Rev2 NOTE: ANSWERS ARE ON PAGE 7. Instructions:

Electrochemical Thermodynamics. Interfaces and Energy Conversion

Lecture 4: Piecewise Cubic Interpolation

Engineering Tensors. Friday November 16, h30 -Muddy Charles. A BEH430 review session by Thomas Gervais.

UNIVERSITY OF IOANNINA DEPARTMENT OF ECONOMICS. M.Sc. in Economics MICROECONOMIC THEORY I. Problem Set II

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS

Jean Fernand Nguema LAMETA UFR Sciences Economiques Montpellier. Abstract

Quadrilateral et Hexahedral Pseudo-conform Finite Elements

An Introduction to Support Vector Machines

? plate in A G in

Investigation phase in case of Bragg coupling

COMPLEX NUMBERS INDEX

Electrical double layer: revisit based on boundary conditions

Quiz: Experimental Physics Lab-I

Solving Singularly Perturbed Differential Difference Equations via Fitted Method

Principle Component Analysis

Jens Siebel (University of Applied Sciences Kaiserslautern) An Interactive Introduction to Complex Numbers

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Review of linear algebra. Nuno Vasconcelos UCSD

INTRODUCTION TO COMPLEX NUMBERS

Electromagnetic modeling of a lightning rod

Vectors and Tensors. R. Shankar Subramanian. R. Aris, Vectors, Tensors, and the Equations of Fluid Mechanics, Prentice Hall (1962).

STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS

NUMERICAL MODELLING OF A CILIUM USING AN INTEGRAL EQUATION

International Journal of Pure and Applied Sciences and Technology

LOCAL FRACTIONAL LAPLACE SERIES EXPANSION METHOD FOR DIFFUSION EQUATION ARISING IN FRACTAL HEAT TRANSFER

Many-Body Calculations of the Isotope Shift

Proof that if Voting is Perfect in One Dimension, then the First. Eigenvector Extracted from the Double-Centered Transformed

338 A^VÇÚO 1n ò Lke n Mancn (211), we make te followng assumpton to control te beavour of small jumps. Assumpton 1.1 L s symmetrc α-stable, were α (,

Torsion, Thermal Effects and Indeterminacy

Solution for singularly perturbed problems via cubic spline in tension

High resolution entropy stable scheme for shallow water equations

Deformation analysis of functionally graded beams by the direct approach

Applied Statistics Qualifier Examination

Entrance and Wall Conduction Effects in Parallel Flow Heat Exchangers

STATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION

DCDM BUSINESS SCHOOL NUMERICAL METHODS (COS 233-8) Solutions to Assignment 3. x f(x)

Two Coefficients of the Dyson Product

Appendix B. The Finite Difference Scheme

THE COMBINED SHEPARD ABEL GONCHAROV UNIVARIATE OPERATOR

ESCI 342 Atmospheric Dynamics I Lesson 1 Vectors and Vector Calculus

Lecture 36. Finite Element Methods

Integral Transforms and Dual Integral Equations to Solve Heat Equation with Mixed Conditions

CHEMICAL ENGINEERING

Convexity preserving interpolation by splines of arbitrary degree

Remember: Project Proposals are due April 11.

Evaluation of Liquefaction Return Period for Bangalore Based on Standard Penetration Test Data: Performance Based Approach

A Hybrid Variational Iteration Method for Blasius Equation

Continuous Time Markov Chain

Module 2. Random Processes. Version 2 ECE IIT, Kharagpur

Note: Please use the actual date you accessed this material in your citation.

523 P a g e. is measured through p. should be slower for lesser values of p and faster for greater values of p. If we set p*

Module 3: Element Properties Lecture 1: Natural Coordinates

Tokyo Institute of Technology Periodic Sequencing Control over Multi Communication Channels with Packet Losses

Mechanics Research Communications

Effectiveness and Efficiency Analysis of Parallel Flow and Counter Flow Heat Exchangers

2.12 Pull Back, Push Forward and Lie Time Derivatives

Effect of Uniform Horizontal Magnetic Field on Thermal Convection in a Rotating Fluid Saturating a Porous Medium

Work and Energy (Work Done by a Varying Force)

PART 1: VECTOR & TENSOR ANALYSIS

Reactor Control Division BARC Mumbai India

On a nonlinear compactness lemma in L p (0, T ; B).

3. Quasi-Stationary Electrodynamics

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

THREE-PARAMETER ELASTIC FOUNDATION MODEL OF FRP STRENGTHENED CONCRETE BEAMS

SE Story Shear Frame. Final Project. 2 Story Bending Beam. m 2. u 2. m 1. u 1. m 3. u 3 L 3. Given: L 1 L 2. EI ω 1 ω 2 Solve for m 2.

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD

A Theoretical Study on the Rank of the Integral Operators for Large- Scale Electrodynamic Analysis

A Family of Multivariate Abel Series Distributions. of Order k

Effects of polarization on the reflected wave

JEL Classification: D50; D58; F10; F11

Lecture Note 3. Eshelby s Inclusion II

Solution Manual. for. Fracture Mechanics. C.T. Sun and Z.-H. Jin

Magnetized Dust Fluid Tilted Universe for Perfect. Fluid Distribution in General Relativity

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

The Study of Lawson Criterion in Fusion Systems for the

Chemical Reaction Engineering

6. Chemical Potential and the Grand Partition Function

THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES

Inductance Calculation for Conductors of Arbitrary Shape

ANALOG CIRCUIT SIMULATION BY STATE VARIABLE METHOD

Lecture 5.8 Flux Vector Splitting

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Chemical Reaction Engineering

The Schur-Cohn Algorithm

Online Appendix to. Mandating Behavioral Conformity in Social Groups with Conformist Members

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Buckling analysis of single-layered FG nanoplates on elastic substrate with uneven porosities and various boundary conditions

HAMILTON-JACOBI TREATMENT OF LAGRANGIAN WITH FERMIONIC AND SCALAR FIELD

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

Formal solvers of the RT equation

More metrics on cartesian products

8.1 Arc Length. What is the length of a curve? How can we approximate it? We could do it following the pattern we ve used before

Transcription:

8 TH INTERNTIONL CONFERENCE ON COMPOSITE MTERILS DYNMIC PROPGTION OF WEK-DISCONTINUOUS INTERFCE CRCK IN FUNCTIONLLY GRDED LYERS UNDER NTI-PLNE SHER J.W. Sn *, Y.S. Lee, S.C. Km, I.H. Hwng 3 Subsystem Deprtment, Kore erospce Reserc Insttute, Dejeon, Kore Deprtment of Mecncl Desgn Engneerng, Cungnm Ntonl Unversty, Dejeon, Kore 3 Rotorcrft Progrm Offce, Kore erospce Reserc Insttute, Dejeon, Kore * Correspondng utor(jeongdl@kr.re.kr Keywords: Functonlly grded mterl (FGM, Interfce movng crck, Wek-dscontnuous nterfce, Dynmc energy relese rte (DERR Introducton Functonlly grded mterls (FGMs ve been recently ttrcted extensve ttenton for use n g temperture pplctons nd wer-protectve cotngs. Te FGMs re mcroscopclly nonomogeneous becuse te mecncl propertes of te FGM vry smootly nd contnuously. Te frcture bevor of te FGM s mportnt desgn ssue. crck n te FGM my exbt complex bevor becuse of te vrton of te mecncl propertes of te mterl. Te frcture bevors of te FGMs ve been studed wdely for bot sttc nd dynmc problems. But, soluton of te dynmc crck propgton of n nterfce crck between two dssmlr functonlly grded lyers s not been presented. In ts pper, dynmc propgton of n nterfce Grfft crck between two functonlly grded lyers under nt-plne ser s nlyzed. Te propertes of te functonlly grded lyers vry contnuously long te tckness. In ddton, te propertes of te two functonlly grded lyers vry dfferently nd te two lyers re connected wekdscontnuously (L et l.[]. Te Yoffe-type model (Yoffe[] for crck propgton s dopted. Fourer trnsform s used to reduce te problem to dul ntegrl euton, wc s ten expressed n Fredolm ntegrl euton of te second knd. Numercl results of te dynmc energy relese rte re presented grpclly to sow te effect of grdent of mterl propertes, crck movng velocty, nd tckness of lyers. Problem sttement nd formultons Consder two functonlly grded lyers contnng fnte nterfce crck subjected to nt-plne ser lodng, s sown n Fg.. Mterl propertes of te two functonlly grded lyers vry dfferently. Te Crtesn coordntes ( X, Y, Z re fxed for te reference. Te FGM lyers occupy te regon, < X <, Y, nd re tck enoug n te Z -drecton. Te crck s stuted long te nterfce lne ( X, Y =. Due to te symmetry n geometry nd lodng, t s suffcent to consder te rgt-nd lf body only. Y vt y τ X x μ μ μ Fg.. Geometry of constnt movng nterfce crck between two functonlly grded lyers We ssume tt te propertes of te two functonlly grded lyers vry contnuously long te tckness nd re smplfed s follows (Delle nd Erdogn[3]: Y e β μ μ Y e β ρ ρ = ( = ( μ nd ρ re te ser modulus nd mterl densty, respectvely. μ nd ρ re te mterl

propertes t te nterfce nd β s te nonomogeneous mterl constnt. Subscrpt ( =, stnds for te upper nd lower lyers, respectvely. Te boundry vlue problem s smplfed consderbly f we consder only te out-of-plne dsplcement suc tt u X = u Y =, uz = w ( X, Y, t (3 u k ( k = X, Y, Z s te dsplcements. In ts cse, te consttutve relton becomes σ ( = μ (4 Zj X, Y, t w, j σ Zj ( j = X, Y s te stress component. Te dynmc nt-plne governng euton for FGM s smplfed to w w μ w + β μ = ρ (5 Y t = X + Y. By substtutng Es. ( nd ( nto te E. (5, te dynmc governng euton s trnsformed nto te followng euton: w + = Y c w w β (6 t c = μ (7 ρ nd c s te ser wve velocty. For te problem of movng crck wt constnt velocty v long te X -drecton, t s convenent to ntroduce Gllen trnsformton suc s x = X v t, y = Y, z = Z, t = t (8 ( x, y, z s te trnsltng coordnte system ttced to te center of te movng crck. In te trnsformed coordnte system, te dynmc nt-plne governng euton for FGM cn be smplfed to te followng form: w w w α + + β = (9 x y y ( v α = c ( Fourer trnsform s ppled to E. (9, nd te result s s follows: w π y y [ ( s e + ( s e ] = cos( sx ds ( β β = δ +, = δ ( β δ = α s + (3 4 nd re te unknowns to be solved. By substtutng E. ( for E. (4, we ve te followng: σ yz = μ [ ( s e π + ( s e y y ] cos( sx ds Te boundry condtons cn be wrtten s σ yz = τ ( x < (4 (5 w + = w ( < x (6 σ x, + = σ ( < x (7 yz ( yz σ x, = σ ( < x (8 yz ( yz = τ s te unform ser trcton. By pplyng te edge lodng condtons of E. (8, te unknowns n E. (4 re evluted s follows: ( s e + ( s e = (9 ( s e + ( s e ( = Te contnuty condton of E. (7 leds to te followng relton between te unknowns: s + ( s = ( s + ( ( ( s

DYNMIC PROPGTION OF WEK-DISCONTINUOUS INTERFCE CRCK IN FUNCTIONLLY GRDED LYDERS UNDER NTI-PLNE SHER It s convenent to use te followng defntons: s ( s + ( s ( s = ( ( ( s Usng te Es. (9 to (, we cn obtn te followng reltons: π τ Ω( ξ + K( ξ, η Ω( η dη = (3 α μ [ F( s ] J ( sη J ( sξ ds K ξ, η = η s (3 ( δ ( e ( s = ( (3 δ δ δ δ ( e ( + e + ( e ( + e s ( s δ δ e ( e ( s δ δ δ δ ( e ( + e + ( e ( + e ( s δ δ e ( e ( s δ δ δ δ ( e ( + e + ( e ( + e ( s δ ( e ( s δ δ δ δ ( e ( + e + ( e ( + e = (4 = (5 = (6 Te mxed boundry condtons of Es. (5 nd (6 led to dul ntegrl eutons n te followng form: π τ s F( s ( s cos( sx ds = ( x < α μ ( s cos( sx ds = ( < x (7 We ntroduce te followng dmensonless vrbles nd functon for numercl nlyss s follows: S Β s =, β =, Β Δ β =, Δ δ =, δ = (3 η = Η, ξ = Ξ (33 π τ Ψ( Ξ Ω( ξ = α μ Ξ (34 F( s = α s ( e δ ( + e ( e δ δ + ( e ( e δ δ ( + δ e (8 Te dul ntegrl E. (7 my be solved by usng new functon Ω (ξ defned by ( s = ξ Ω( ξ J ( sξ dξ (9 J s te zero-order Bessel functon of te frst knd. By nsertng E. (9 nto E. (7, we cn fnd tt te uxlry functon Ω (ξ s gven by Fredolm ntegrl euton of te second knd n te followng form: By substtutng Es. (3 to (34 for Es. (3 nd (3, we cn obtn Fredolm ntegrl euton of te second knd n te followng form: Ψ( Ξ + L ( Ξ, Η Ψ( Η dη = Ξ (35 3

L ( Ξ, Η = S ΞΗ S F J ( S Η J ( S Ξ ds (36 Δ Δ ( ( S QQQQ e e F( = α S Δ Δ Δ Δ QQ( e ( Q + Q e + QQ ( e ( Q + Q e (37 Β Β Β Β Q = Δ +, Q = Δ +, Q = Δ, Q = Δ (38 Te mode III dynmc stress ntensty fctor K III (v nd dynmc energy relese rte (v re defned nd determned n te followng forms: ( v = τ π Ψ( (39 K III π ( v = τ Ψ ( (4 μ E. (4 cn be mde dmensonless s follow: ( v μ = Ψ ( τ π (4 s benefcl to ncrese of te resstnce of te nterfce crck propgton of FGM. π GIII(v μ/τ 3..5..5. v/c =.4 / = / =. B =. B =. B =. B = -. B = -. n wc te functon Ψ ( cn be clculted from E. (35. 3 Dscussons To nvestgte te effect of te grdent of mterl propertes, crck movng velocty nd tckness of lyers on te dynmc energy relese rte (DERR, numercl nlyses re crred out. Fg. dsplys te vrton of te normlzed DERR ( v μ τ π gnst te normlzed nonomogeneous mterl constnt of te upper lyer Β wt vrous normlzed non-omogeneous mterl constnts of te lower lyer Β t v / c =.4 nd / = / =.. Te normlzed DERR decreses wen te grdents of mterl propertes of te upper nd lower lyers ncrese. For te upper lyer, te grdent of mterl propertes ncreses s te non-omogeneous mterl constnt ncreses. But for te lower lyer, te grdent of mterl propertes ncreses s te non-omogeneous mterl constnt decreses becuse vlue of te y-xs s negtve. Increse of te grdent of mterl propertes from te nterfce.5. -. -.... B Fg.. Vrton of te normlzed DERR ( vμ τ π wt Β Fg. 3 sows te vrton of te normlzed DERR gnst te normlzed crck movng velocty v / c wt te vrous normlzed non-omogeneous mterl constnts. ccordng to te vlues of te grdent of mterl propertes of te upper nd lower lyer, we cn clssfy nto tree ctegores s follows: Cse I : mterl propertes ncrese wen te tckness of upper nd lower lyer ncreses from te nterfce, Cse II : mterl propertes decrese wen te tckness of upper nd lower lyer ncreses from te nterfce, Cse III : mterl propertes ncrese from te lower surfce ( y = to upper surfce

DYNMIC PROPGTION OF WEK-DISCONTINUOUS INTERFCE CRCK IN FUNCTIONLLY GRDED LYDERS UNDER NTI-PLNE SHER ( y =, nd vce verse. For te Cse II nd III, te normlzed DERR ncreses s te crck movng velocty ncreses. But for te Cse I, te trend s opposte. Te normlzed DERR decreses wen te crck movng velocty ncreses. Tt s, ncrese of te stffness from te nterfce to te upper nd lower surfce s elpful to ncrese of te resstnce of te nterfce crck propgton of FGM. s te bsolute vlue of Β Β ncreses, tt s, te dfference between te grdents of te mterl propertes of upper nd lower lyer s gettng bgger, te normlzed DERR ncreses or decreses more rpdly. π GIII(v μ/τ 3..5..5..5 / = / = B = -., B =. B =., B =. B =., B =. B =., B = -. B =., B = -.....4.6 v/c Fg. 3. Vrton of te normlzed DERR ( v μ τ π wt v / c Te effect of te crck movng velocty on te vrton of te normlzed DERR s sown n Fg. 4 wt vrous tcknesses of te lyers. Te normlzed DERR ncreses wt te ncrese of te crck movng velocty. But te normlzed DERR decreses wen te tckness of lyer ncreses. Increse of te tckness of FGM lyer s lso benefcl to ncrese of te resstnce of te nterfce crck propgton of FGM. Fg. 5 presents te vrton of te normlzed DERR gnst te normlzed tckness of te lower FGM lyer wt te vrous non-omogeneous mterl constnts. Smlr to Fg. 4, te normlzed DERR decreses s te tckness of te lower lyer ncreses. But, over certn vlue of te tckness of te lower lyer (bout 3., te effect of decrese of te normlzed DERR s neglgble. s seen n Cse I of Fg. 3, Fg. 5 lso sows tt ncrese of te stffness from te nterfce to te upper nd lower surfce s elpful to ncrese of te resstnce of te nterfce crck propgton of FGM. π GIII(v μ/τ π GIII(v μ/τ...8 Β =., B =. / =., / =. / =., / =. / =., / =..6...4.6 v/c Fg. 4. Vrton of te normlzed DERR ( v μ τ π wt v / c.6..8.4. v/c =., / =. B =., B =. B =., B =. B =., B =. B =., B = -. B =., B = -..6.. 3. 4. 5. 6. 7. 8. 9.. / Fg. 5. Vrton of te normlzed DERR ( vμ τ π wt / 4 Conclusons Te problem of dynmc propgton of wek- 5

dscontnuous nterfce crck between two functonlly grded lyers under nt-plne ser lodng ws nlyzed by te ntegrl trnsform pproc. Te ser modulus nd mss densty of te FGM vry contnuously long te tckness. Fredolm ntegrl euton s solved numerclly. Te computed results sow tt te followngs re elpful to ncrese of te resstnce of te nterfce crck propgton of FGM: Increse of te grdent of mterl propertes, b Increse of te mterl propertes from te nterfce to te upper nd lower free surfce, c Increse of te tckness of FGM lyer. Te normlzed DERR ncreses or decreses wt ncrese of crck movng velocty. References [] Y.D. L, B. J, N. Zng, L.Q. Tng nd Y. D, Dynmc stress ntensty fctor of te wek/mcrodscontnuous nterfce crck of FGM cotng. Int. J. Solds Strut., Vol. 43, pp 4795-489, 6. [] E.H. Yoffe, Te movng Grfft crck. Plos. Mgzne 7, Vol. 4, pp 739-75, 95. [3] F. Delle nd F. Erdogn, Te crck problem for nonomgeneous plne. SME J. ppl. Mec., Vol. 5, pp 69-64, 983.