SHEAR LAG MODELLING OF THERMAL STRESSES IN UNIDIRECTIONAL COMPOSITES

Size: px
Start display at page:

Download "SHEAR LAG MODELLING OF THERMAL STRESSES IN UNIDIRECTIONAL COMPOSITES"

Transcription

1 ORA/POSER REFERENCE: ICF00374OR SHEAR AG MODEING OF HERMA SRESSES IN UNIDIRECIONA COMPOSIES Chad M. adis Departet o Mechaical Egieerig ad Materials Sciece MS 3 Rice Uiversity P.O. Box 89 Housto X 7705 ladis@rice.edu ph: ax: ABSRAC A two-diesioal shear lag odel is preseted to aalyze the steady state distributios o stress ad teperature i uidirectioally reiorced coposites. Equatios allowig or variatios i axial displaceet ad teperature alog ay give iber or atrix regio are developed. he derivatio o the goverig equatios is greatly sipliied by the assuptios that displaceets perpedicular to the iber directio ca be igored the axial displaceet ad teperature are uior over the cross sectio o ay iber ad the distributios o teperature ad axial displaceet are biliear withi the atrix regios. hese equatios are the solved or a coiguratio with uiorly spaced atrix cracks ad used to deterie the eective Youg s odulus ad theral coductivity o the cracked aterial. KEYWORDS Coposites theral stresses shear lag. INRODUCION he aalysis o stress ad strai i ulti-iber coposites has bee acilitated by the developet o shear lag odels. hese shear lag odels treat ibers as oe-diesioal load carryig structures that traser loads to oe aother through shear stresses withi the atrix aterial. he origial ulti-iber shear lag odel was orulated by Hedgepeth[] i 96. Sice the a uber o researchers have ehaced this odel by icludig eects like slidig at the iber atrix iterace or the load carryig capability o the atrix see adis ad McMeekig[]. A direct two-diesioal extesio o the Hedgepeth[] odel was proposed by Beyerlei ad adis[3] to iclude the eects o atrix stiess that reduces to the Hedgepeth[] odel i the appropriate liit. he work preseted i this paper will build o the odel o Beyerlei ad adis[3] to iclude the eects o steady state teperature distributios. hese shear lag

2 odels are ost useul i large scale coposite ailure siulatios like those carried out by Ibabdeljalil ad Curti[4] ad adis et al.[5] where ore detailed stress calculatios would be uteable.. GOVERNING EQUAIONS he odel syste cosidered is depicted i Figure with the ibers separated by the atrix i a periodic ashio. Figure : wo diesioal iber/atrix coiguratio. he Youg s odulus area theral coductivity ad theral expasio coeiciet o the ibers are E A k ad. he Youg s odulus shear odulus area logitudial theral coductivity trasverse theral coductivity theral expasio coeiciet ad width o the atrix regios are E G A k k ad W. he thickess o the syste is t hece A = Wt. Followig the procedure outlied by adis ad co-workers [-3] the equatios goverig the axial displaceet ad teperature alog the th iber u ad ad the ceter o the th atrix regio u ad are η d η d η d k ( )= 0 dx dx dx () η d η d η d k ( )= 0 () 3 dx dx dx ρ du ρ du ρ du ρ d ρ d d ( ) = G u u u 3 dx dx dx dx dx dx ρ du ρ du ρ du ρ d ρ d d G u u u ( )= 3 dx dx dx 3 dx dx dx (3) (4) EA ka Gt ρ= η= G = = E A ka E A W k kt (5) kaw

3 he essetial assuptios required to obtai these equatio are that the aterial is costraied to displace i the x directio oly shear deoratio ad trasverse teperature gradiets i the ibers are eglected ad the displaceet or teperature proile alog the trasverse directio i the atrix is biliear. 3. N-INDEPENDEN SOUIONS I geeral Eqs. (-4) represet a large set o coupled ordiary dieretial equatios or the displaceet ad teperature i each iber ad atrix regio. However i the locatios o iperectios i.e. iber or atrix cracks are distributed uiorly over all iber or atrix regios the Eqs. (-4) reduce to our goverig ordiary dieretial equatios. I other words i the distributios o displaceet ad teperature are idepedet o the iber or atrix regio uber the the ollowig equatios gover the syste: η η 4k 0 ( )= 6 (6) η η 4k( )= 0 (7) ρ ρ ρ 4 = ρ u u G u u (8) ρ ρ ρ ρ u u 4G( u u)= (9) where u ad u are the teperature ad displaceet distributios i all ibers ad all atrix regios. he ' ad the " deote irst ad secod derivatives with respect to x. he geeral solutios to Eqs. (6-9) are u = u β x = q x C e C e βx 0 0 η β x η = 0 q0x C e C e η η βx ρ ρ u = u0 0x ε0x q0x ρ ρ x x x C e C e C e C e β β β β u u u ux ρ ρ ρ q x ε x q x 8G ρ ρ ρ x x ρ x ρ C e C e Cu e Cu e ρ ρ u ux β β β β (0) () () (3) βηρ η ρ 4 ρ C C 48G = (4) βρ β ( η ) 4 48 ( ρ )( ρ ) G

4 [ ] βηρ η ρ η η ( ρ ) ( ρ 4) C C 48G = (5) βρ β ( η ) 4 48 ( ρ )( ρ ) G β 48k η = kw η 4 48G ρ βu = EW ρ 4 (6) where 0 q 0 u 0 ε 0 C C C u ad C u ust be deteried ro boudary coditios. he geeral solutio listed i Eqs. (0-6) ca be used to deterie the eective Youg s odulus ad eective logitudial theral coductivity o a coposite with a uior atrix crack spacig. Cosider the crack geoetry illustrated i Figure. Figure : Uior atrix crack spacig coiguratio. he atrix cracks are assued to be tractio-ree ad perectly isulatig. o aalyze the eective theral coductivity o the cracked coposite oly the teperature solutio eeds to be cosidered ad to deterie the eective odulus oly the displaceet solutio eeds to be cosidered. Sice the distributio o teperature or displaceet is biliear across a atrix regio oly the average heat lux or stress at a atrix crack ca be speciied to be zero. So or the theral coductivity calculatio the average teperature gradiet is set to zero at x=0 ad at x=. he solutio or the teperature distributio i a give seget o aterial betwee two atrix cracks is η coshβxcoshβ x = 0 q0 x β sihβ (7) = 0 q0 x β coshβ( x ) coshβx sihβ (8) is the average teperature across the atrix at a give x locatio. he eective where = theral coductivity o the coposite is the just the average heat lux through the coposite divided by the average teperature gradiet across the sae regio. Note that care ust be take to iclude teperature jups across the atrix cracks i these averages. he eective theral coductivity o the cracked coposite is

5 k c ka ka η coshβ = A A β sihβ (9) Eq. (9) is idetical to the result obtaied by u ad Hutchiso [6]. Fro a solutio siilar to Eqs. (7-8) or the displaceets the eective Youg s odulus is E c EA E A ρ coshβu = A A βu sihβu (0) Notice that whe the crack spacig goes to iiity the the coposite coductivity or odulus is the area ractio weighted average o the ibers ad atrix. Whe goes to zero the the coductivity ad odulus o the coposite are kc = ka ( A A) ad Ec = EA ( A A). Aother siple solutio or the uior atrix crack coiguratio is the distributio o stress i the ibers ad atrix whe the teperature o the coposite is uior. he iber stress σ ad the average atrix stress σ are σ = E σ = E ρ sihβu( x ) sihβux ρ sihβu sihβu( x ) sihβux ρ sihβu () () Notice sice there is o applied load oly teperature that σ A σ A =0. 4. DISCUSSION he siple solutios preseted i Eqs. (0-6) do ot exploit the true value o the shear lag odel. Geerally solutios to the equatios o theroelasticity could be obtaied with series solutios or the siple geoetry cosidered i Figure leadig to ore accurate estiates o the coposite coductivity ad odulus. However it is diicult to obtai exact solutios or the distributios o stress ad teperature aroud eve a sigle iber break. Cotrarily siple superpositio techiques exists that ca solve Eqs. (-5) or ay iite uber o arbitrarily located iber or atrix cracks. Hece the sipliied shear lag odel is ost useul i large ailure siulatios like those carried out i [4] ad [5] where rapid solutios or the iteractios o cracks with oe aother ad the loadig are iportat. REFERENCES. Hedgepeth J.M. (96) NASA N D-8.. adis C.M. ad McMeekig R.M. (998) It. J. Solids Structures Beyerlei I.J. ad adis C.M. (999) Mechaics o Materials Ibabdeljalil M. ad Curti W.A. (997) It. J. Solids Structures adis C.M. Beyerlei I.J. ad McMeekig R.M. (000) J. Mech. Phys. Solids u.j. ad Hutchiso J.W. (995) Philos. ras. R. Soc. odo Ser. A

Supplementary Information

Supplementary Information Suppleetary Iforatio -Breakdow of cotiuu fracture echaics at the aoscale- Takahiro Shiada,,* Keji Ouchi, Yuu Chihara, ad Takayuki Kitaura Departet of echaical Egieerig ad Sciece, Kyoto Uiversity, Nishikyo-ku,

More information

Lecture 11. Solution of Nonlinear Equations - III

Lecture 11. Solution of Nonlinear Equations - III Eiciecy o a ethod Lecture Solutio o Noliear Equatios - III The eiciecy ide o a iterative ethod is deied by / E r r: rate o covergece o the ethod : total uber o uctios ad derivative evaluatios at each step

More information

A Steady State Heat Conduction Problem in. a Thick Annular Disc Due to Arbitrary. Axisymmetric Heat Flux

A Steady State Heat Conduction Problem in. a Thick Annular Disc Due to Arbitrary. Axisymmetric Heat Flux Noliear Aalysis ad Differetial Equatios, Vol. 4, 016, o. 3, 11-131 HIKARI Ltd, www.-hikari.co http://dx.doi.org/10.1988/ade.016.51037 A Steady State Heat Coductio Proble i a Thick Aular Disc Due to Arbitrary

More information

Example of CLT for Symmetric Laminate with Mechanical Loading

Example of CLT for Symmetric Laminate with Mechanical Loading Exaple of CL for Syetric Laate with Mechaical Loadg hese are two probles which detail the process through which the laate deforatios are predicted usg Classical Laatio theory. he first oly cludes echaical

More information

The Effect of Fiber Breakage on Transient Stress Distribution in a Single-Lap Joint Composite Material

The Effect of Fiber Breakage on Transient Stress Distribution in a Single-Lap Joint Composite Material Joural o Solid Mechaics Vol. 7, o. 4 (015 pp. 44-457 The Eect o Fiber Breakage o Trasiet Stress Distributio i a Sigle-Lap Joit Coposite Material M. Shishehsaz 1, *, S. Yaghoubi 1 Departet o Mechaical Egieerig,

More information

X. Perturbation Theory

X. Perturbation Theory X. Perturbatio Theory I perturbatio theory, oe deals with a ailtoia that is coposed Ĥ that is typically exactly solvable of two pieces: a referece part ad a perturbatio ( Ĥ ) that is assued to be sall.

More information

Application of Homotopy Analysis Method for Solving various types of Problems of Ordinary Differential Equations

Application of Homotopy Analysis Method for Solving various types of Problems of Ordinary Differential Equations Iteratioal Joural o Recet ad Iovatio Treds i Coputig ad Couicatio IN: 31-8169 Volue: 5 Issue: 5 16 Applicatio of Hootopy Aalysis Meod for olvig various types of Probles of Ordiary Differetial Equatios

More information

AVERAGE MARKS SCALING

AVERAGE MARKS SCALING TERTIARY INSTITUTIONS SERVICE CENTRE Level 1, 100 Royal Street East Perth, Wester Australia 6004 Telephoe (08) 9318 8000 Facsiile (08) 95 7050 http://wwwtisceduau/ 1 Itroductio AVERAGE MARKS SCALING I

More information

Birth-Death Processes. Outline. EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Relationship Among Stochastic Processes.

Birth-Death Processes. Outline. EEC 686/785 Modeling & Performance Evaluation of Computer Systems. Relationship Among Stochastic Processes. EEC 686/785 Modelig & Perforace Evaluatio of Couter Systes Lecture Webig Zhao Deartet of Electrical ad Couter Egieerig Clevelad State Uiversity webig@ieee.org based o Dr. Raj jai s lecture otes Relatioshi

More information

Lecture Outline. 2 Separating Hyperplanes. 3 Banach Mazur Distance An Algorithmist s Toolkit October 22, 2009

Lecture Outline. 2 Separating Hyperplanes. 3 Banach Mazur Distance An Algorithmist s Toolkit October 22, 2009 18.409 A Algorithist s Toolkit October, 009 Lecture 1 Lecturer: Joatha Keler Scribes: Alex Levi (009) 1 Outlie Today we ll go over soe of the details fro last class ad ake precise ay details that were

More information

International Journal of Mathematical Archive-4(9), 2013, 1-5 Available online through ISSN

International Journal of Mathematical Archive-4(9), 2013, 1-5 Available online through   ISSN Iteratioal Joural o Matheatical Archive-4(9), 03, -5 Available olie through www.ija.io ISSN 9 5046 THE CUBIC RATE OF CONVERGENCE OF GENERALIZED EXTRAPOLATED NEWTON RAPHSON METHOD FOR SOLVING NONLINEAR

More information

Lecture 20 - Wave Propagation Response

Lecture 20 - Wave Propagation Response .09/.093 Fiite Eleet Aalysis of Solids & Fluids I Fall 09 Lecture 0 - Wave Propagatio Respose Prof. K. J. Bathe MIT OpeCourseWare Quiz #: Closed book, 6 pages of otes, o calculators. Covers all aterials

More information

Orthogonal Functions

Orthogonal Functions Royal Holloway Uiversity of odo Departet of Physics Orthogoal Fuctios Motivatio Aalogy with vectors You are probably failiar with the cocept of orthogoality fro vectors; two vectors are orthogoal whe they

More information

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1

(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1 ROOT LOCUS TECHNIQUE 93 should be desiged differetly to eet differet specificatios depedig o its area of applicatio. We have observed i Sectio 6.4 of Chapter 6, how the variatio of a sigle paraeter like

More information

ECE 901 Lecture 4: Estimation of Lipschitz smooth functions

ECE 901 Lecture 4: Estimation of Lipschitz smooth functions ECE 9 Lecture 4: Estiatio of Lipschitz sooth fuctios R. Nowak 5/7/29 Cosider the followig settig. Let Y f (X) + W, where X is a rado variable (r.v.) o X [, ], W is a r.v. o Y R, idepedet of X ad satisfyig

More information

Principle Of Superposition

Principle Of Superposition ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give

More information

Problem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =

Problem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t = Mathematics Summer Wilso Fial Exam August 8, ANSWERS Problem 1 (a) Fid the solutio to y +x y = e x x that satisfies y() = 5 : This is already i the form we used for a first order liear differetial equatio,

More information

Section 5. Gaussian Imagery

Section 5. Gaussian Imagery OPTI-01/0 Geoetrical ad Istruetal Optics 5-1 Sectio 5 Gaussia Iagery Iagig Paraxial optics provides a sipliied etodology to deterie ray pats troug optical systes. Usig tis etod, te iage locatio or a geeral

More information

Statistics and Data Analysis in MATLAB Kendrick Kay, February 28, Lecture 4: Model fitting

Statistics and Data Analysis in MATLAB Kendrick Kay, February 28, Lecture 4: Model fitting Statistics ad Data Aalysis i MATLAB Kedrick Kay, kedrick.kay@wustl.edu February 28, 2014 Lecture 4: Model fittig 1. The basics - Suppose that we have a set of data ad suppose that we have selected the

More information

The Hypergeometric Coupon Collection Problem and its Dual

The Hypergeometric Coupon Collection Problem and its Dual Joural of Idustrial ad Systes Egieerig Vol., o., pp -7 Sprig 7 The Hypergeoetric Coupo Collectio Proble ad its Dual Sheldo M. Ross Epstei Departet of Idustrial ad Systes Egieerig, Uiversity of Souther

More information

Balancing. Rotating Components Examples of rotating components in a mechanism or a machine. (a)

Balancing. Rotating Components Examples of rotating components in a mechanism or a machine. (a) alacig NOT COMPLETE Rotatig Compoets Examples of rotatig compoets i a mechaism or a machie. Figure 1: Examples of rotatig compoets: camshaft; crakshaft Sigle-Plae (Static) alace Cosider a rotatig shaft

More information

CHAPTER INTRODUCTION

CHAPTER INTRODUCTION 39 CHAPTER 3 LOCAL NON-SIMILARITY SOLUTION FOR THE IMPACT OF CHEMICAL REACTION ON MHD MIXED CONVECTION HEAT AND MASS TRANSFER FLOW OVER A POROUS WEDGE IN THE PRESENCE OF SUCTION / INJECTION 3. INTRODUCTION

More information

PAijpam.eu ON TENSOR PRODUCT DECOMPOSITION

PAijpam.eu ON TENSOR PRODUCT DECOMPOSITION Iteratioal Joural of Pure ad Applied Mathematics Volume 103 No 3 2015, 537-545 ISSN: 1311-8080 (prited versio); ISSN: 1314-3395 (o-lie versio) url: http://wwwijpameu doi: http://dxdoiorg/1012732/ijpamv103i314

More information

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES

PARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES Diola Bagayoko (0 PARTAL DFFERENTAL EQUATONS SEPARATON OF ARABLES. troductio As discussed i previous lectures, partial differetial equatios arise whe the depedet variale, i.e., the fuctio, varies with

More information

Binomial transform of products

Binomial transform of products Jauary 02 207 Bioial trasfor of products Khristo N Boyadzhiev Departet of Matheatics ad Statistics Ohio Norther Uiversity Ada OH 4580 USA -boyadzhiev@ouedu Abstract Give the bioial trasfors { b } ad {

More information

Perturbation Theory, Zeeman Effect, Stark Effect

Perturbation Theory, Zeeman Effect, Stark Effect Chapter 8 Perturbatio Theory, Zeea Effect, Stark Effect Ufortuately, apart fro a few siple exaples, the Schrödiger equatio is geerally ot exactly solvable ad we therefore have to rely upo approxiative

More information

Engineering Mechanics Dynamics & Vibrations. Engineering Mechanics Dynamics & Vibrations Plane Motion of a Rigid Body: Equations of Motion

Engineering Mechanics Dynamics & Vibrations. Engineering Mechanics Dynamics & Vibrations Plane Motion of a Rigid Body: Equations of Motion 1/5/013 Egieerig Mechaics Dyaics ad Vibratios Egieerig Mechaics Dyaics & Vibratios Egieerig Mechaics Dyaics & Vibratios Plae Motio of a Rigid Body: Equatios of Motio Motio of a rigid body i plae otio is

More information

An Algorithmist s Toolkit October 20, Lecture 11

An Algorithmist s Toolkit October 20, Lecture 11 18.409 A Algorithist s Toolkit October 20, 2009 Lecture 11 Lecturer: Joatha Keler Scribe: Chaithaya Badi 1 Outlie Today we ll itroduce ad discuss Polar of a covex body. Correspodece betwee or fuctios ad

More information

Section 7. Gaussian Reduction

Section 7. Gaussian Reduction 7- Sectio 7 Gaussia eductio Paraxial aytrace Equatios eractio occurs at a iterace betwee two optical spaces. The traser distace t' allows the ray height y' to be determied at ay plae withi a optical space

More information

Phys. 201 Mathematical Physics 1 Dr. Nidal M. Ershaidat Doc. 12

Phys. 201 Mathematical Physics 1 Dr. Nidal M. Ershaidat Doc. 12 Physics Departmet, Yarmouk Uiversity, Irbid Jorda Phys. Mathematical Physics Dr. Nidal M. Ershaidat Doc. Fourier Series Deiitio A Fourier series is a expasio o a periodic uctio (x) i terms o a iiite sum

More information

11 Correlation and Regression

11 Correlation and Regression 11 Correlatio Regressio 11.1 Multivariate Data Ofte we look at data where several variables are recorded for the same idividuals or samplig uits. For example, at a coastal weather statio, we might record

More information

9.5 Young s Double-Slit Experiment

9.5 Young s Double-Slit Experiment 9.5 Youg s Double-Slit Experiet Physics Tool box Early attepts to deostrate the iterferece of light were usuccessful because the two sources were too far apart ad out of phase, ad the wavelegth of light

More information

Where do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway?

Where do eigenvalues/eigenvectors/eigenfunctions come from, and why are they important anyway? Where do eigevalues/eigevectors/eigeuctios come rom, ad why are they importat ayway? I. Bacgroud (rom Ordiary Dieretial Equatios} Cosider the simplest example o a harmoic oscillator (thi o a vibratig strig)

More information

Summer MA Lesson 13 Section 1.6, Section 1.7 (part 1)

Summer MA Lesson 13 Section 1.6, Section 1.7 (part 1) Suer MA 1500 Lesso 1 Sectio 1.6, Sectio 1.7 (part 1) I Solvig Polyoial Equatios Liear equatio ad quadratic equatios of 1 variable are specific types of polyoial equatios. Soe polyoial equatios of a higher

More information

2. F ; =(,1)F,1; +F,1;,1 is satised by thestirlig ubers of the rst kid ([1], p. 824). 3. F ; = F,1; + F,1;,1 is satised by the Stirlig ubers of the se

2. F ; =(,1)F,1; +F,1;,1 is satised by thestirlig ubers of the rst kid ([1], p. 824). 3. F ; = F,1; + F,1;,1 is satised by the Stirlig ubers of the se O First-Order Two-Diesioal Liear Hoogeeous Partial Dierece Equatios G. Neil Have y Ditri A. Gusev z Abstract Aalysis of algoriths occasioally requires solvig of rst-order two-diesioal liear hoogeeous partial

More information

Closed virial equation-of-state for the hard-disk fluid

Closed virial equation-of-state for the hard-disk fluid Physical Review Letter LT50 receipt of this auscript 6 Jue 00 Closed virial equatio-of-state for the hard-disk fluid Athoy Beris ad Leslie V. Woodcock Departet of Cheical Egieerig Colbur Laboratory Uiversity

More information

Zeros of Polynomials

Zeros of Polynomials Math 160 www.timetodare.com 4.5 4.6 Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered with fidig the solutios of polyomial equatios of ay degree

More information

) is a square matrix with the property that for any m n matrix A, the product AI equals A. The identity matrix has a ii

) is a square matrix with the property that for any m n matrix A, the product AI equals A. The identity matrix has a ii square atrix is oe that has the sae uber of rows as colus; that is, a atrix. he idetity atrix (deoted by I, I, or [] I ) is a square atrix with the property that for ay atrix, the product I equals. he

More information

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,

3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense, 3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [

More information

Wavelet Transform Theory. Prof. Mark Fowler Department of Electrical Engineering State University of New York at Binghamton

Wavelet Transform Theory. Prof. Mark Fowler Department of Electrical Engineering State University of New York at Binghamton Wavelet Trasfor Theory Prof. Mark Fowler Departet of Electrical Egieerig State Uiversity of New York at Bighato What is a Wavelet Trasfor? Decopositio of a sigal ito costituet parts Note that there are

More information

Lecture 10: Bounded Linear Operators and Orthogonality in Hilbert Spaces

Lecture 10: Bounded Linear Operators and Orthogonality in Hilbert Spaces Lecture : Bouded Liear Operators ad Orthogoality i Hilbert Spaces 34 Bouded Liear Operator Let ( X, ), ( Y, ) i i be ored liear vector spaces ad { } X Y The, T is said to be bouded if a real uber c such

More information

SECTION 2 Electrostatics

SECTION 2 Electrostatics SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial

More information

Dynamic Response of Second Order Mechanical Systems with Viscous Dissipation forces

Dynamic Response of Second Order Mechanical Systems with Viscous Dissipation forces Hadout #b (pp. 4-55) Dyamic Respose o Secod Order Mechaical Systems with Viscous Dissipatio orces M X + DX + K X = F t () Periodic Forced Respose to F (t) = F o si( t) ad F (t) = M u si(t) Frequecy Respose

More information

Evaluation of Bessel Functions Using a Computer Program

Evaluation of Bessel Functions Using a Computer Program Evaluatio of Bessel Fuctios Usig a Coputer Progra P. S. Yeh, Ph.D. Abstract I cylidrical coordiate, there are two types of Bessel fuctios. These fuctios are the Bessel fuctio ad the odified Bessel fuctio.

More information

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes

Bernoulli Polynomials Talks given at LSBU, October and November 2015 Tony Forbes Beroulli Polyoials Tals give at LSBU, October ad Noveber 5 Toy Forbes Beroulli Polyoials The Beroulli polyoials B (x) are defied by B (x), Thus B (x) B (x) ad B (x) x, B (x) x x + 6, B (x) dx,. () B 3

More information

Define a Markov chain on {1,..., 6} with transition probability matrix P =

Define a Markov chain on {1,..., 6} with transition probability matrix P = Pla Group Work 0. The title says it all Next Tie: MCMC ad Geeral-state Markov Chais Midter Exa: Tuesday 8 March i class Hoework 4 due Thursday Uless otherwise oted, let X be a irreducible, aperiodic Markov

More information

CSCI-6971 Lecture Notes: Stochastic processes

CSCI-6971 Lecture Notes: Stochastic processes CSCI-6971 Lecture Notes: Stochastic processes Kristopher R. Beevers Departet of Coputer Sciece Resselaer Polytechic Istitute beevek@cs.rpi.edu February 2, 2006 1 Overview Defiitio 1.1. A stochastic process

More information

Temperature Distribution in a Cylindrical Core and Coil Assembly with Heat Generation

Temperature Distribution in a Cylindrical Core and Coil Assembly with Heat Generation Temperature Distributio i a Cylidrical Core ad Coil Assembly with Heat Geeratio P. S. Yeh, Ph.D. Abstract A residetial distributio trasformer cosists of a cylidrical steel tak ad a core-ad-coil assembly

More information

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS http://www.paper.edu.c Iteratioal Joural of Bifurcatio ad Chaos, Vol. 1, No. 5 () 119 15 c World Scietific Publishig Compay AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC

More information

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.

DETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS. ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics

More information

Contents Two Sample t Tests Two Sample t Tests

Contents Two Sample t Tests Two Sample t Tests Cotets 3.5.3 Two Saple t Tests................................... 3.5.3 Two Saple t Tests Setup: Two Saples We ow focus o a sceario where we have two idepedet saples fro possibly differet populatios. Our

More information

LECTURE 14. Non-linear transverse motion. Non-linear transverse motion

LECTURE 14. Non-linear transverse motion. Non-linear transverse motion LETURE 4 No-liear trasverse motio Floquet trasformatio Harmoic aalysis-oe dimesioal resoaces Two-dimesioal resoaces No-liear trasverse motio No-liear field terms i the trajectory equatio: Trajectory equatio

More information

SECTION 2.6 THE SECOND ALTERNATIVE

SECTION 2.6 THE SECOND ALTERNATIVE 54 SECTION 2.6 THE SECOND ALTERNATIVE We ow discuss the probles where the Secod Alterative holds. The suppositio is that there is a otrivial solutio for L(y) =, B (y) = B 2 (y) =. The Fredhol Theores assure

More information

REDUCING THE POSSIBILITY OF SUBJECTIVE ERROR IN THE DETERMINATION OF THE STRUCTURE-FUNCTION-BASED EFFECTIVE THERMAL CONDUCTIVITY OF BOARDS

REDUCING THE POSSIBILITY OF SUBJECTIVE ERROR IN THE DETERMINATION OF THE STRUCTURE-FUNCTION-BASED EFFECTIVE THERMAL CONDUCTIVITY OF BOARDS Nice, Côte d Azur, Frace, 27-29 Septeber 2006 REDUCING THE POSSIBILITY OF SUBJECTIVE ERROR IN THE DETERMINATION OF THE STRUCTURE-FUNCTION-BASED EFFECTIVE THERMAL CONDUCTIVITY OF BOARDS Erő Kollár, Vladiír

More information

A Hadamard-type lower bound for symmetric diagonally dominant positive matrices

A Hadamard-type lower bound for symmetric diagonally dominant positive matrices A Hadamard-type lower boud for symmetric diagoally domiat positive matrices Christopher J. Hillar, Adre Wibisoo Uiversity of Califoria, Berkeley Jauary 7, 205 Abstract We prove a ew lower-boud form of

More information

Accuracy. Computational Fluid Dynamics. Computational Fluid Dynamics. Computational Fluid Dynamics

Accuracy. Computational Fluid Dynamics. Computational Fluid Dynamics. Computational Fluid Dynamics http://www.d.edu/~gtryggva/cfd-course/ Computatioal Fluid Dyamics Lecture Jauary 3, 7 Grétar Tryggvaso It is clear that although the umerical solutio is qualitatively similar to the aalytical solutio,

More information

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j

The z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.

More information

Flight and Orbital Mechanics. Exams

Flight and Orbital Mechanics. Exams 1 Flight ad Orbital Mechaics Exas Exa AE2104-11: Flight ad Orbital Mechaics (2 Noveber 2012, 14.00 17.00) Please put your ae, studet uber ad ALL YOUR INITIALS o your work. Aswer all questios ad put your

More information

Supplementary Material

Supplementary Material Suppleetary Material Wezhuo Ya a0096049@us.edu.s Departet of Mechaical Eieeri, Natioal Uiversity of Siapore, Siapore 117576 Hua Xu pexuh@us.edu.s Departet of Mechaical Eieeri, Natioal Uiversity of Siapore,

More information

We have also learned that, thanks to the Central Limit Theorem and the Law of Large Numbers,

We have also learned that, thanks to the Central Limit Theorem and the Law of Large Numbers, Cofidece Itervals III What we kow so far: We have see how to set cofidece itervals for the ea, or expected value, of a oral probability distributio, both whe the variace is kow (usig the stadard oral,

More information

APPLIED MULTIVARIATE ANALYSIS

APPLIED MULTIVARIATE ANALYSIS ALIED MULTIVARIATE ANALYSIS FREQUENTLY ASKED QUESTIONS AMIT MITRA & SHARMISHTHA MITRA DEARTMENT OF MATHEMATICS & STATISTICS INDIAN INSTITUTE OF TECHNOLOGY KANUR X = X X X [] The variace covariace atrix

More information

Definitions and Theorems. where x are the decision variables. c, b, and a are constant coefficients.

Definitions and Theorems. where x are the decision variables. c, b, and a are constant coefficients. Defiitios ad Theorems Remember the scalar form of the liear programmig problem, Miimize, Subject to, f(x) = c i x i a 1i x i = b 1 a mi x i = b m x i 0 i = 1,2,, where x are the decisio variables. c, b,

More information

Salmon: Lectures on partial differential equations. 3. First-order linear equations as the limiting case of second-order equations

Salmon: Lectures on partial differential equations. 3. First-order linear equations as the limiting case of second-order equations 3. First-order liear equatios as the limitig case of secod-order equatios We cosider the advectio-diffusio equatio (1) v = 2 o a bouded domai, with boudary coditios of prescribed. The coefficiets ( ) (2)

More information

Double Diffusive Convection from a Permeable Horizontal Cylinder of Elliptic Cross Section in a Saturated Porous Medium

Double Diffusive Convection from a Permeable Horizontal Cylinder of Elliptic Cross Section in a Saturated Porous Medium Proceedigs o the 4th WSEAS It. Co. o HEAT TRANSFER, THERMAL ENGINEERING ad ENVIRONMENT, Elouda, Greece, August -3, 6 (pp357-36) Double Diusive Covectio rom a Permeable Horizotal Cylider o Elliptic Cross

More information

SOLUTIONS: ECE 606 Homework Week 7 Mark Lundstrom Purdue University (revised 3/27/13) e E i E T

SOLUTIONS: ECE 606 Homework Week 7 Mark Lundstrom Purdue University (revised 3/27/13) e E i E T SOUIONS: ECE 606 Homework Week 7 Mark udstrom Purdue Uiversity (revised 3/27/13) 1) Cosider a - type semicoductor for which the oly states i the badgap are door levels (i.e. ( E = E D ). Begi with the

More information

The Differential Transform Method for Solving Volterra s Population Model

The Differential Transform Method for Solving Volterra s Population Model AASCIT Couicatios Volue, Issue 6 Septeber, 15 olie ISSN: 375-383 The Differetial Trasfor Method for Solvig Volterra s Populatio Model Khatereh Tabatabaei Departet of Matheatics, Faculty of Sciece, Kafas

More information

Crank-Nicolson Implicit Method For The Nonlinear Schrodinger Equation With Variable Coefficient

Crank-Nicolson Implicit Method For The Nonlinear Schrodinger Equation With Variable Coefficient Crak-Nicolso Iplicit Method For The Noliear Schrodiger Equatio With Variable Coefficiet Yaa Yee Choy a Wooi Nee Ta b Ki Gaik Tay c ad Chee Tiog Og d a Faculty of Sciece Techology & Hua Developet Uiversiti

More information

Question 1: The magnetic case

Question 1: The magnetic case September 6, 018 Corell Uiversity, Departmet of Physics PHYS 337, Advace E&M, HW # 4, due: 9/19/018, 11:15 AM Questio 1: The magetic case I class, we skipped over some details, so here you are asked to

More information

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.

The axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c. 5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =

More information

THE THEORETICAL BACKGROUND FOR CALCULATION OF GEOMETRICAL PARAMETERS OF THE STABLE MELTING ZONE IN SINGLE-ROW COKE CUPOLAS

THE THEORETICAL BACKGROUND FOR CALCULATION OF GEOMETRICAL PARAMETERS OF THE STABLE MELTING ZONE IN SINGLE-ROW COKE CUPOLAS METALLURGY AND FOUNDRY ENGINEERING Vol. 3, 005, No. W³adys³aw Loga *, Belkheir Khatei ** THE THEORETICAL BACKGROUND FOR CALCULATION OF GEOMETRICAL PARAMETERS OF THE STABLE MELTING ZONE IN SINGLE-ROW COKE

More information

Randomized Algorithms I, Spring 2018, Department of Computer Science, University of Helsinki Homework 1: Solutions (Discussed January 25, 2018)

Randomized Algorithms I, Spring 2018, Department of Computer Science, University of Helsinki Homework 1: Solutions (Discussed January 25, 2018) Radomized Algorithms I, Sprig 08, Departmet of Computer Sciece, Uiversity of Helsiki Homework : Solutios Discussed Jauary 5, 08). Exercise.: Cosider the followig balls-ad-bi game. We start with oe black

More information

A New Type of q-szász-mirakjan Operators

A New Type of q-szász-mirakjan Operators Filoat 3:8 07, 567 568 https://doi.org/0.98/fil7867c Published by Faculty of Scieces ad Matheatics, Uiversity of Niš, Serbia Available at: http://www.pf.i.ac.rs/filoat A New Type of -Szász-Miraka Operators

More information

Topic 9 - Taylor and MacLaurin Series

Topic 9 - Taylor and MacLaurin Series Topic 9 - Taylor ad MacLauri Series A. Taylors Theorem. The use o power series is very commo i uctioal aalysis i act may useul ad commoly used uctios ca be writte as a power series ad this remarkable result

More information

Numerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem

Numerical Simulation of Thermomechanical Problems in Applied Mechanics: Application to Solidification Problem Leoardo Joural of Scieces ISSN 1583-0233 Issue 9, July-December 2006 p. 25-32 Numerical Simulatio of Thermomechaical Problems i Applied Mechaics: Applicatio to Solidificatio Problem Vicet Obiajulu OGWUAGWU

More information

After the completion of this section the student

After the completion of this section the student Chapter II CALCULUS II. Liits ad Cotiuity 55 II. LIMITS AND CONTINUITY Objectives: After the copletio of this sectio the studet - should recall the defiitios of the it of fuctio; - should be able to apply

More information

Discrete probability distributions

Discrete probability distributions Discrete probability distributios I the chapter o probability we used the classical method to calculate the probability of various values of a radom variable. I some cases, however, we may be able to develop

More information

Viscoelastic Flow and Heat Transfer over a Non- Linearly Stretching Sheet: OHAM Solution

Viscoelastic Flow and Heat Transfer over a Non- Linearly Stretching Sheet: OHAM Solution Joural o Applied Fluid echaics, Vol. 9, No., pp. -8, 06. Available olie at.jaolie.et, ISSN 75-57, EISSN 75-645. DOI: 0.8869/acadpub.ja.68.8.4778 Viscoelastic Flo ad Heat Traser over a No- Liearly Stretchig

More information

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series. .3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(

More information

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001. Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;

More information

Lecture 8 & Tutorial 1: Experimental idention of low order model

Lecture 8 & Tutorial 1: Experimental idention of low order model Eurother Advaced School Metti 5 Roscoff Jue 13-18, 211 Lecture 8 & utorial 1: Experietal idetio of low order odel Jea-Luc Battaglia, Laboratory I2M, Departet REFLE, I2M Bordeaux Measureet iversio 35 3

More information

S. A. ALIEV, Y. I. YELEYKO, Y. V. ZHERNOVYI. STEADY-STATE DISTRIBUTIONS FOR CERTAIN MODIFICATIONS OF THE M/M/1/m QUEUEING SYSTEM

S. A. ALIEV, Y. I. YELEYKO, Y. V. ZHERNOVYI. STEADY-STATE DISTRIBUTIONS FOR CERTAIN MODIFICATIONS OF THE M/M/1/m QUEUEING SYSTEM Trasactios of Azerbaija Natioal Acadey of Scieces, Series of Physical-Techical ad Matheatical Scieces: Iforatics ad Cotrol Probles 009 Vol XXIX, 6 P 50-58 S A ALIEV, Y I YELEYKO, Y V ZHERNOVYI STEADY-STATE

More information

Studying Interaction of Cotton-Raw Material with Working Bodies of Cotton-Cleaning Machines

Studying Interaction of Cotton-Raw Material with Working Bodies of Cotton-Cleaning Machines ISSN: 35-38 Iteratioal Joural of AdvacedResearch i Sciece, Egieerig ad Techology Vol. 5, Issue, Deceber 8 Studyig Iteractio of Cotto-Raw Material with Workig Bodies of Cotto-Cleaig Machies R.H. Rosulov,

More information

FREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES

FREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES Abdullah : FREE VIBRATIONS OF SIMPY SUPPORTED BEAMS FREE VIBRATIONS OF SIMPY SUPPORTED BEAMS USING FOURIER SERIES SAWA MUBARAK ABDUAH Assistat ecturer Uiversity of Mosul Abstract Fourier series will be

More information

PROBLEMS AND SOLUTIONS 2

PROBLEMS AND SOLUTIONS 2 PROBEMS AND SOUTIONS Problem 5.:1 Statemet. Fid the solutio of { u tt = a u xx, x, t R, u(x, ) = f(x), u t (x, ) = g(x), i the followig cases: (b) f(x) = e x, g(x) = axe x, (d) f(x) = 1, g(x) =, (f) f(x)

More information

A Method to Calculate the True Stress and True Strain for Tensile Test of Plastic

A Method to Calculate the True Stress and True Strain for Tensile Test of Plastic Ke Egieerig Materials Vols. 74-76 (4) pp. 77-8 olie at http://www.scietific.et 4 Tras Tech Publicatios, Switzerlad Method to Calculate the True Stress ad True Strai for Tesile Test of Plastic X. W. Du,

More information

Streamfunction-Vorticity Formulation

Streamfunction-Vorticity Formulation Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first

More information

Massachusetts Institute of Technology

Massachusetts Institute of Technology 6.0/6.3: Probabilistic Systems Aalysis (Fall 00) Problem Set 8: Solutios. (a) We cosider a Markov chai with states 0,,, 3,, 5, where state i idicates that there are i shoes available at the frot door i

More information

Exercise 8 CRITICAL SPEEDS OF THE ROTATING SHAFT

Exercise 8 CRITICAL SPEEDS OF THE ROTATING SHAFT Exercise 8 CRITICA SEEDS OF TE ROTATING SAFT. Ai of the exercise Observatio ad easureet of three cosecutive critical speeds ad correspodig odes of the actual rotatig shaft. Copariso of aalytically coputed

More information

Linear regression. Daniel Hsu (COMS 4771) (y i x T i β)2 2πσ. 2 2σ 2. 1 n. (x T i β y i ) 2. 1 ˆβ arg min. β R n d

Linear regression. Daniel Hsu (COMS 4771) (y i x T i β)2 2πσ. 2 2σ 2. 1 n. (x T i β y i ) 2. 1 ˆβ arg min. β R n d Liear regressio Daiel Hsu (COMS 477) Maximum likelihood estimatio Oe of the simplest liear regressio models is the followig: (X, Y ),..., (X, Y ), (X, Y ) are iid radom pairs takig values i R d R, ad Y

More information

Answer Key, Problem Set 1, Written

Answer Key, Problem Set 1, Written Cheistry 1 Mies, Sprig, 018 Aswer Key, Proble Set 1, Writte 1. 14.3;. 14.34 (add part (e): Estiate / calculate the iitial rate of the reactio); 3. NT1; 4. NT; 5. 14.37; 6. 14.39; 7. 14.41; 8. NT3; 9. 14.46;

More information

Physical bases of dental material science

Physical bases of dental material science Physical bases of detal aterial sciece The ioizatio eergy of the sodiu ato is 496 kj/ol. How large eergy is ecessary i ev to ioize oe ato? The ioizatio eergy of a sigle ato: E i 19 496000 19 8,7 10 ε i

More information

Double Derangement Permutations

Double Derangement Permutations Ope Joural of iscrete Matheatics, 206, 6, 99-04 Published Olie April 206 i SciRes http://wwwscirporg/joural/ojd http://dxdoiorg/04236/ojd2066200 ouble erageet Perutatios Pooya aeshad, Kayar Mirzavaziri

More information

A 2nTH ORDER LINEAR DIFFERENCE EQUATION

A 2nTH ORDER LINEAR DIFFERENCE EQUATION A 2TH ORDER LINEAR DIFFERENCE EQUATION Doug Aderso Departmet of Mathematics ad Computer Sciece, Cocordia College Moorhead, MN 56562, USA ABSTRACT: We give a formulatio of geeralized zeros ad (, )-discojugacy

More information

Microscopic traffic flow modeling

Microscopic traffic flow modeling Chapter 34 Microscopic traffic flow modelig 34.1 Overview Macroscopic modelig looks at traffic flow from a global perspective, whereas microscopic modelig, as the term suggests, gives attetio to the details

More information

CS537. Numerical Analysis and Computing

CS537. Numerical Analysis and Computing CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be

More information

ECE Spring Prof. David R. Jackson ECE Dept. Notes 20

ECE Spring Prof. David R. Jackson ECE Dept. Notes 20 ECE 6341 Sprig 016 Prof. David R. Jackso ECE Dept. Notes 0 1 Spherical Wave Fuctios Cosider solvig ψ + k ψ = 0 i spherical coordiates z φ θ r y x Spherical Wave Fuctios (cot.) I spherical coordiates we

More information

Acoustic Field inside a Rigid Cylinder with a Point Source

Acoustic Field inside a Rigid Cylinder with a Point Source Acoustic Field iside a Rigid Cylider with a Poit Source 1 Itroductio The ai objectives of this Deo Model are to Deostrate the ability of Coustyx to odel a rigid cylider with a poit source usig Coustyx

More information

CS321. Numerical Analysis and Computing

CS321. Numerical Analysis and Computing CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca

More information

Some remarks on the paper Some elementary inequalities of G. Bennett

Some remarks on the paper Some elementary inequalities of G. Bennett Soe rears o the paper Soe eleetary iequalities of G. Beett Dag Ah Tua ad Luu Quag Bay Vieta Natioal Uiversity - Haoi Uiversity of Sciece Abstract We give soe couterexaples ad soe rears of soe of the corollaries

More information

Formation of A Supergain Array and Its Application in Radar

Formation of A Supergain Array and Its Application in Radar Formatio of A Supergai Array ad ts Applicatio i Radar Tra Cao Quye, Do Trug Kie ad Bach Gia Duog. Research Ceter for Electroic ad Telecommuicatios, College of Techology (Coltech, Vietam atioal Uiversity,

More information