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

Size: px
Start display at page:

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

Transcription

1 STATIC ANALYSIS OF TWO-LERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION Ákos József Lengyel István Ecsed Assstant Lecturer Emertus Professor Insttute of Appled Mechancs Unversty of Mskolc Mskolc-Egyetemváros H-355 Hungary e-mal: Abstract The object of the present paper s to analyse the statc behavour of the pezoelectrc two-layered beams wth nterlayer slp The Euler-Bernoull hypothess s assumed to hold for each layer and a lnear consttutve equaton between the horzontal slp and the nterlamnar shear force s consdered For each beam component the appled electrc feld s constant and the prescrbed mechancal loads are not present INTRODUCTION Pezoelectrc beams are commonly used n many engneerng applcatons such as control elements and transducers Pezoelectrc beams n unmorph and bmorph confguratons n partcular are wdely used for systems where actuaton and/or sensng are needed [] The present paper deals wth the statc analyss of twolayered pezoelectrc beams wth nterlayer slp The connecton of the beam components n aal drecton s weak but n normal drecton s perfect There s no separaton n normal drecton between the beam components The nterlayer slp n aal drecton s defned on the common boundary of beam components as the dfference of the aal components of the dsplacement feld The consdered twolayered pezoelectrc beam confguraton s shown n Fg The cross-secton of beam component B s a rectangle A of whch sdes are b h ( ) The centre of the cross-secton A s C ( ) and the Y -weghted centre of the whole crosssecton s C as shown n Fg The Young modulus of beam component B s denoted by Y ( ) From Fg t follows that c CC c c CC c c c c A Y A hb A h b () () Accordng to the Euler-Bernoull beam theory the dsplacements are as follows [3] dw u u ( ) z v w w( ) ( y z) B ( ) (3)

2 constant Fgure Two-layered pezoelectrc beam wth weak shear connecton where u v and w are the dsplacements n y and z drectons The net consttutve equaton wll be used Y e3 Ez (4) Here s the normal stress s the normal stran E s the z component of the z electrc feld vector and e s the pezoelectrc constant of the thckness polarzed 3 beam component E n terms of potental z s obtaned as (Fg ) E E (5) z z h h n beam components B and B The top and bottom surfaces of pezoelectrc layers are metalzed Net we use e3 e ( ) By the use of Eqs (4) (5) the epresson of the normal stress can be gven as du d w Y z e ( ) ( ) y z B h (6) We defne the followng stress resultants By smple computatons we obtan (7) N d A M z d A ( ) A A du d w N ( ) c eb ( ) (8) du d w M ( ) c IY ( ) e cb ( ) (9)

3 Here ( ) There s no appled mechancal load n aal drecton I z da A ths means that N N N () The nterlayer slp s s obtaned as the dfference of aal component of dsplacements defned on the common boundary of beam components B and B that s s( ) u ( ) u ( ) () From Eqs () and () we have du A Y ds b e e () du ds b e e Substtutng of Eqs () (3) nto Eq (8) gves (3) Here N c c e c e c ds d w b N c c e c e c ds d w b A Y (4) (5) (6) From Eqs (9) () (3) we can derve the epresson of the total bendng moment actng on the whole cross-secton A A A as ds d w M M M c IY ec ec b IY I Y IY (7) The whole cross-secton s loaded by the shear force V whch s d d d V c IY (8) 3 M s w 3 The nterlayer shear force S actng n aal drecton s a lnear functon of the slp ts epresson s S ks (9)

4 where k s the slp modulus [3] The force equlbrum equaton for beam component B can be formulated as 3 dn d s d w ks c ks 3 () Combnaton of Eq (8) wth Eq () yelds where d s V IY s c () IY k IY IY c () IY In problem shown n Fg V and s() thus we have for s s( ) s( ) K snh (3) The constant K s obtaned from the net boundary condton (Fg ) DETERMINATION OF SLIP AND DEFLECTION On the whole beam M ( ) From ths equaton we get Integraton of Eq (5) provdes N ( L) (4) d w ds IY c c e ce b (5) that s dw dw IY IY c s( ) s() b ce ce (6) dw IY c s( ) b ce ce (7) snce dw (Fg ) A repeated ntegraton gves the result

5 cosh ( ) IYw ck bc e c e (8) snce w() Lengthy but elementary computatons yeld IY ds b N( L) ce ce IY L c (9) Combnaton of Eq (3) wth Eq (9) leads to the formula of K b c e K c L c e cosh (3) The fnal epresson of the slp functon can be wrtten n the net form s b c e c e snh c cosh L ( ) (3) Substtutng K nto equaton (8) provdes the deflecton of two-layered pezoelectrc beam n terms of appled voltage b( ce ce ) cosh w( ) IY cosh L (3) c For bmorph beam h h c c e e e we have eb snh s( ) cosh L ecb cosh w( ) IY coshl (33) (34) 3 NUMERICAL EXAMPLE The followng data have been used n numercal computatons: h h b 5 m L 5 m Y Y 5 GPa e e 76 N/mV 4 The slp functons and deflecton functons for k Pa k 6 8 k k m V Pa 3 Pa 4 Pa k5 Pa are determned The graphs of slp and deflecton functons for k k k 3 k and 4 k are shown n Fg and Fg 3 5 respectvely

6 Fgure Plots of slp functons Fgure 3 Plots of deflecton functons 4 CONCLUSIONS Ths paper presents an analytcal soluton for two-layered pezoelectrc cantlever beam wth mperfect shear connecton Smple formulas are derved to obtan the deflecton and slp functons An eample llustrates the applcaton of derved formulas and the effect of slp modulus to the deformaton and slp ACKNOWLEDGEMENTS Ths research was supported by the Natonal Research Development and Innovaton Offce NKFIH K57 REFERENCES [] YANG J: The Mechancs of Pezoelectrc Structures World Scentfc Publshng New Jersey 6 [] YANG J: An Introducton to the Theory of Pezoelecrcty Sprnger Berln 6 [3] ECSEDI I BAKSA A: Statc analyss of composte beams wth weak shear connecton Appled Mathematcal Modellng 35(4)

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,

More information

COMPLEMENTARY ENERGY METHOD FOR CURVED COMPOSITE BEAMS

COMPLEMENTARY ENERGY METHOD FOR CURVED COMPOSITE BEAMS ultscence - XXX. mcocd Intenatonal ultdscplnay Scentfc Confeence Unvesty of skolc Hungay - pl 06 ISBN 978-963-358-3- COPLEENTRY ENERGY ETHOD FOR CURVED COPOSITE BES Ákos József Lengyel István Ecsed ssstant

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam

APPENDIX F A DISPLACEMENT-BASED BEAM ELEMENT WITH SHEAR DEFORMATIONS. Never use a Cubic Function Approximation for a Non-Prismatic Beam APPENDIX F A DISPACEMENT-BASED BEAM EEMENT WITH SHEAR DEFORMATIONS Never use a Cubc Functon Approxmaton for a Non-Prsmatc Beam F. INTRODUCTION { XE "Shearng Deformatons" }In ths appendx a unque development

More information

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods

Chapter Eight. Review and Summary. Two methods in solid mechanics ---- vectorial methods and energy methods or variational methods Chapter Eght Energy Method 8. Introducton 8. Stran energy expressons 8.3 Prncpal of statonary potental energy; several degrees of freedom ------ Castglano s frst theorem ---- Examples 8.4 Prncpal of statonary

More information

Finite Element Modelling of truss/cable structures

Finite Element Modelling of truss/cable structures Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures

More information

In this section is given an overview of the common elasticity models.

In this section is given an overview of the common elasticity models. Secton 4.1 4.1 Elastc Solds In ths secton s gven an overvew of the common elastcty models. 4.1.1 The Lnear Elastc Sold The classcal Lnear Elastc model, or Hooean model, has the followng lnear relatonshp

More information

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

THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS OF A TELESCOPIC HYDRAULIC CYLINDER SUBJECTED TO EULER S LOAD Journal of Appled Mathematcs and Computatonal Mechancs 7, 6(3), 7- www.amcm.pcz.pl p-issn 99-9965 DOI:.75/jamcm.7.3. e-issn 353-588 THE EFFECT OF TORSIONAL RIGIDITY BETWEEN ELEMENTS ON FREE VIBRATIONS

More information

One Dimensional Axial Deformations

One Dimensional Axial Deformations One Dmensonal al Deformatons In ths secton, a specfc smple geometr s consdered, that of a long and thn straght component loaded n such a wa that t deforms n the aal drecton onl. The -as s taken as the

More information

Application to Plane (rigid) frame structure

Application to Plane (rigid) frame structure Advanced Computatonal echancs 18 Chapter 4 Applcaton to Plane rgd frame structure 1. Dscusson on degrees of freedom In case of truss structures, t was enough that the element force equaton provdes onl

More information

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: Instructor s Name and Secton: (Crcle Your Secton) Sectons:

More information

Geometrically exact multi-layer beams with a rigid interconnection

Geometrically exact multi-layer beams with a rigid interconnection Geometrcally exact mult-layer beams wth a rgd nterconnecton Leo Škec, Gordan Jelenć To cte ths verson: Leo Škec, Gordan Jelenć. Geometrcally exact mult-layer beams wth a rgd nterconnecton. 2nd ECCOMAS

More information

Indeterminate pin-jointed frames (trusses)

Indeterminate pin-jointed frames (trusses) Indetermnate pn-jonted frames (trusses) Calculaton of member forces usng force method I. Statcal determnacy. The degree of freedom of any truss can be derved as: w= k d a =, where k s the number of all

More information

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

Buckling analysis of single-layered FG nanoplates on elastic substrate with uneven porosities and various boundary conditions IOSR Journal of Mechancal and Cvl Engneerng (IOSR-JMCE) e-issn: 78-1684,p-ISSN: 30-334X, Volume 15, Issue 5 Ver. IV (Sep. - Oct. 018), PP 41-46 www.osrjournals.org Bucklng analyss of sngle-layered FG nanoplates

More information

Increase Decrease Remain the Same (Circle one) (2 pts)

Increase Decrease Remain the Same (Circle one) (2 pts) ME 270 Sample Fnal Eam PROBLEM 1 (25 ponts) Prob. 1 questons are all or nothng. PROBLEM 1A. (5 ponts) FIND: A 2000 N crate (D) s suspended usng ropes AB and AC and s n statc equlbrum. If θ = 53.13, determne

More information

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments.

Frame element resists external loads or disturbances by developing internal axial forces, shear forces, and bending moments. CE7 Structural Analyss II PAAR FRAE EEET y 5 x E, A, I, Each node can translate and rotate n plane. The fnal dsplaced shape has ndependent generalzed dsplacements (.e. translatons and rotatons) noled.

More information

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. ME 270 Fall 2013 Fnal Exam NME (Last, Frst): Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: INSTRUCTIONS

More information

829. An adaptive method for inertia force identification in cantilever under moving mass

829. An adaptive method for inertia force identification in cantilever under moving mass 89. An adaptve method for nerta force dentfcaton n cantlever under movng mass Qang Chen 1, Mnzhuo Wang, Hao Yan 3, Haonan Ye 4, Guola Yang 5 1,, 3, 4 Department of Control and System Engneerng, Nanng Unversty,

More information

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES

EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES EVALUATION OF THE VISCO-ELASTIC PROPERTIES IN ASPHALT RUBBER AND CONVENTIONAL MIXES Manuel J. C. Mnhoto Polytechnc Insttute of Bragança, Bragança, Portugal E-mal: mnhoto@pb.pt Paulo A. A. Perera and Jorge

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

FUZZY FINITE ELEMENT METHOD

FUZZY FINITE ELEMENT METHOD FUZZY FINITE ELEMENT METHOD RELIABILITY TRUCTURE ANALYI UING PROBABILITY 3.. Maxmum Normal tress Internal force s the shear force, V has a magntude equal to the load P and bendng moment, M. Bendng moments

More information

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. ME 270 Sprng 2014 Fnal Exam NME (Last, Frst): Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: INSTRUCTIONS

More information

Chapter 8. Potential Energy and Conservation of Energy

Chapter 8. Potential Energy and Conservation of Energy Chapter 8 Potental Energy and Conservaton of Energy In ths chapter we wll ntroduce the followng concepts: Potental Energy Conservatve and non-conservatve forces Mechancal Energy Conservaton of Mechancal

More information

I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. ME 270 Fall 2012 Fnal Exam Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: INSTRUCTIONS Begn each problem

More information

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve

More information

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

Physics 5153 Classical Mechanics. Principle of Virtual Work-1 P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal

More information

Professor Terje Haukaas University of British Columbia, Vancouver The Q4 Element

Professor Terje Haukaas University of British Columbia, Vancouver  The Q4 Element Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to

More information

Electrical Circuits 2.1 INTRODUCTION CHAPTER

Electrical Circuits 2.1 INTRODUCTION CHAPTER CHAPTE Electrcal Crcuts. INTODUCTION In ths chapter, we brefly revew the three types of basc passve electrcal elements: resstor, nductor and capactor. esstance Elements: Ohm s Law: The voltage drop across

More information

INDETERMINATE STRUCTURES METHOD OF CONSISTENT DEFORMATIONS (FORCE METHOD)

INDETERMINATE STRUCTURES METHOD OF CONSISTENT DEFORMATIONS (FORCE METHOD) INETNTE STUTUES ETHO OF ONSISTENT EFOTIONS (FOE ETHO) If all the support reactons and nternal forces (, Q, and N) can not be determned by usng equlbrum equatons only, the structure wll be referred as STTIY

More information

FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION

FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION VOL. 6, NO. 3, MARCH 0 ISSN 89-6608 006-0 Asan Research Publshng Network (ARPN). All rghts reserved. FINITE DIFFERENCE ANALYSIS OF CURVED DEEP BEAMS ON WINKLER FOUNDATION Adel A. Al-Azzaw and Al S. Shaker

More information

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

UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 2017/2018 FINITE ELEMENT AND DIFFERENCE SOLUTIONS OCD0 UNIVERSITY OF BOLTON RAK ACADEMIC CENTRE BENG(HONS) MECHANICAL ENGINEERING SEMESTER TWO EXAMINATION 07/08 FINITE ELEMENT AND DIFFERENCE SOLUTIONS MODULE NO. AME6006 Date: Wednesda 0 Ma 08 Tme: 0:00

More information

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam. ME 270 Summer 2014 Fnal Exam NAME (Last, Frst): Please revew the followng statement: I certfy that I have not gven unauthorzed ad nor have I receved ad n the completon of ths exam. Sgnature: INSTRUCTIONS

More information

12. The Hamilton-Jacobi Equation Michael Fowler

12. The Hamilton-Jacobi Equation Michael Fowler 1. The Hamlton-Jacob Equaton Mchael Fowler Back to Confguraton Space We ve establshed that the acton, regarded as a functon of ts coordnate endponts and tme, satsfes ( ) ( ) S q, t / t+ H qpt,, = 0, and

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

A Timoshenko Piezoelectric Beam Finite Element with Consistent Performance Irrespective of Geometric and Material Configurations

A Timoshenko Piezoelectric Beam Finite Element with Consistent Performance Irrespective of Geometric and Material Configurations 99 A Tmoshenko Pezoelectrc Beam Fnte Element wth Consstent Performance Irrespectve of Geometrc and Materal Confguratons Abstract The conventonal Tmoshenko pezoelectrc beam fnte elements based on Frst-order

More information

I have not received unauthorized aid in the completion of this exam.

I have not received unauthorized aid in the completion of this exam. ME 270 Sprng 2013 Fnal Examnaton Please read and respond to the followng statement, I have not receved unauthorzed ad n the completon of ths exam. Agree Dsagree Sgnature INSTRUCTIONS Begn each problem

More information

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites

The Two-scale Finite Element Errors Analysis for One Class of Thermoelastic Problem in Periodic Composites 7 Asa-Pacfc Engneerng Technology Conference (APETC 7) ISBN: 978--6595-443- The Two-scale Fnte Element Errors Analyss for One Class of Thermoelastc Problem n Perodc Compostes Xaoun Deng Mngxang Deng ABSTRACT

More information

Transactions of the VŠB Technical University of Ostrava, Mechanical Series. article No. 1907

Transactions of the VŠB Technical University of Ostrava, Mechanical Series. article No. 1907 Transactons of the VŠB Techncal Unversty of Ostrava, Mechancal Seres No., 0, vol. LVIII artcle No. 907 Marek NIKODÝM *, Karel FYDÝŠEK ** FINITE DIFFEENCE METHOD USED FO THE BEAMS ON ELASTIC FOUNDATION

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

A new Approach for Solving Linear Ordinary Differential Equations

A new Approach for Solving Linear Ordinary Differential Equations , ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of

More information

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,

More information

A Mechanics-Based Approach for Determining Deflections of Stacked Multi-Storey Wood-Based Shear Walls

A Mechanics-Based Approach for Determining Deflections of Stacked Multi-Storey Wood-Based Shear Walls A Mechancs-Based Approach for Determnng Deflectons of Stacked Mult-Storey Wood-Based Shear Walls FPINNOVATIONS Acknowledgements Ths publcaton was developed by FPInnovatons and the Canadan Wood Councl based

More information

Electrical double layer: revisit based on boundary conditions

Electrical double layer: revisit based on boundary conditions Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer

More information

GEO-SLOPE International Ltd, Calgary, Alberta, Canada Vibrating Beam

GEO-SLOPE International Ltd, Calgary, Alberta, Canada   Vibrating Beam GEO-SLOPE Internatonal Ltd, Calgary, Alberta, Canada www.geo-slope.com Introducton Vbratng Beam Ths example looks at the dynamc response of a cantlever beam n response to a cyclc force at the free end.

More information

Foundations of Arithmetic

Foundations of Arithmetic Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an

More information

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

THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES THE SMOOTH INDENTATION OF A CYLINDRICAL INDENTOR AND ANGLE-PLY LAMINATES W. C. Lao Department of Cvl Engneerng, Feng Cha Unverst 00 Wen Hwa Rd, Tachung, Tawan SUMMARY: The ndentaton etween clndrcal ndentor

More information

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method Soluton of Lnear System of Equatons and Matr Inverson Gauss Sedel Iteraton Method It s another well-known teratve method for solvng a system of lnear equatons of the form a + a22 + + ann = b a2 + a222

More information

Plate Theories for Classical and Laminated plates Weak Formulation and Element Calculations

Plate Theories for Classical and Laminated plates Weak Formulation and Element Calculations Plate heores for Classcal and Lamnated plates Weak Formulaton and Element Calculatons PM Mohte Department of Aerospace Engneerng Indan Insttute of echnolog Kanpur EQIP School on Computatonal Methods n

More information

Lecture Note 3. Eshelby s Inclusion II

Lecture Note 3. Eshelby s Inclusion II ME340B Elastcty of Mcroscopc Structures Stanford Unversty Wnter 004 Lecture Note 3. Eshelby s Incluson II Chrs Wenberger and We Ca c All rghts reserved January 6, 004 Contents 1 Incluson energy n an nfnte

More information

Physics 4B. A positive value is obtained, so the current is counterclockwise around the circuit.

Physics 4B. A positive value is obtained, so the current is counterclockwise around the circuit. Physcs 4B Solutons to Chapter 7 HW Chapter 7: Questons:, 8, 0 Problems:,,, 45, 48,,, 7, 9 Queston 7- (a) no (b) yes (c) all te Queston 7-8 0 μc Queston 7-0, c;, a;, d; 4, b Problem 7- (a) Let be the current

More information

An influence line shows how the force in a particular member changes as a concentrated load is moved along the structure.

An influence line shows how the force in a particular member changes as a concentrated load is moved along the structure. CE 331, Fall 2010 Influence Les for Trusses 1 / 7 An fluence le shows how the force a partcular member changes as a concentrated load s moved along the structure. For eample, say we are desgng the seven

More information

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive

Visco-Rubber Elastic Model for Pressure Sensitive Adhesive Vsco-Rubber Elastc Model for Pressure Senstve Adhesve Kazuhsa Maeda, Shgenobu Okazawa, Koj Nshgch and Takash Iwamoto Abstract A materal model to descrbe large deformaton of pressure senstve adhesve (PSA

More information

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(

More information

Lecture 16 Statistical Analysis in Biomaterials Research (Part II)

Lecture 16 Statistical Analysis in Biomaterials Research (Part II) 3.051J/0.340J 1 Lecture 16 Statstcal Analyss n Bomaterals Research (Part II) C. F Dstrbuton Allows comparson of varablty of behavor between populatons usng test of hypothess: σ x = σ x amed for Brtsh statstcan

More information

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two

Cubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order

More information

MODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava

MODELLING OF ELASTO-STATICS OF POWER LINES BY NEW COMPOSITE BEAM FINITE ELEMENT Bratislava ODING OF ASTO-STATICS OF POW INS BY NW COPOSIT BA FINIT NT urín Justín 1 rabovský Jura 1 Gogola oman 1 utš Vladmír 1 Paulech Jura 1 1 Insttute of Automotve echatroncs FI STU n Bratslava Ilkovčova 3 812

More information

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods

A Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods Appled Mathematcal Scences, Vol. 11, 2017, no. 52, 2579-2586 HIKARI Ltd, www.m-hkar.com https://do.org/10.12988/ams.2017.79280 A Soluton of the Harry-Dym Equaton Usng Lattce-Boltzmannn and a Soltary Wave

More information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY. Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR

More information

The Finite Element Method

The Finite Element Method The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of

More information

Computers and Mathematics with Applications

Computers and Mathematics with Applications Computers and Mathematcs wth Applcatons 56 (2008 3204 3220 Contents lsts avalable at ScenceDrect Computers and Mathematcs wth Applcatons journal homepage: www.elsever.com/locate/camwa An nnovatve egenvalue

More information

Constitutive Modelling of Superplastic AA-5083

Constitutive Modelling of Superplastic AA-5083 TECHNISCHE MECHANIK, 3, -5, (01, 1-6 submtted: September 19, 011 Consttutve Modellng of Superplastc AA-5083 G. Gulano In ths study a fast procedure for determnng the constants of superplastc 5083 Al alloy

More information

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

More information

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST

NUMERICAL RESULTS QUALITY IN DEPENDENCE ON ABAQUS PLANE STRESS ELEMENTS TYPE IN BIG DISPLACEMENTS COMPRESSION TEST Appled Computer Scence, vol. 13, no. 4, pp. 56 64 do: 10.23743/acs-2017-29 Submtted: 2017-10-30 Revsed: 2017-11-15 Accepted: 2017-12-06 Abaqus Fnte Elements, Plane Stress, Orthotropc Materal Bartosz KAWECKI

More information

Tensor Smooth Length for SPH Modelling of High Speed Impact

Tensor Smooth Length for SPH Modelling of High Speed Impact Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru

More information

Flow Induced Vibration

Flow Induced Vibration Flow Induced Vbraton Project Progress Report Date: 16 th November, 2005 Submtted by Subhrajt Bhattacharya Roll no.: 02ME101 Done under the gudance of Prof. Anrvan Dasgupta Department of Mechancal Engneerng,

More information

NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM

NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM Advanced Steel Constructon Vol. 5, No., pp. 59-7 (9) 59 NONLINEAR NATURAL FREQUENCIES OF A TAPERED CANTILEVER BEAM M. Abdel-Jaber, A.A. Al-Qasa,* and M.S. Abdel-Jaber Department of Cvl Engneerng, Faculty

More information

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown

More information

A boundary element method with analytical integration for deformation of inhomogeneous elastic materials

A boundary element method with analytical integration for deformation of inhomogeneous elastic materials Journal of Physcs: Conference Seres PAPER OPEN ACCESS A boundary element method wth analytcal ntegraton for deformaton of nhomogeneous elastc materals To cte ths artcle: Moh. Ivan Azs et al 2018 J. Phys.:

More information

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

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

Lecture 2: Numerical Methods for Differentiations and Integrations

Lecture 2: Numerical Methods for Differentiations and Integrations Numercal Smulaton of Space Plasmas (I [AP-4036] Lecture 2 by Lng-Hsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton

More information

SIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD

SIMULATION OF WAVE PROPAGATION IN AN HETEROGENEOUS ELASTIC ROD SIMUATION OF WAVE POPAGATION IN AN HETEOGENEOUS EASTIC OD ogéro M Saldanha da Gama Unversdade do Estado do o de Janero ua Sào Francsco Xaver 54, sala 5 A 559-9, o de Janero, Brasl e-mal: rsgama@domancombr

More information

GENERAL THEORY FOR THE MODELLING OF LAMINATED MULTI-LAYER PLATES AND ITS APPLICATION FOR THE ANALYSIS OF SANDWICH PANELS WITH PLY DROPS

GENERAL THEORY FOR THE MODELLING OF LAMINATED MULTI-LAYER PLATES AND ITS APPLICATION FOR THE ANALYSIS OF SANDWICH PANELS WITH PLY DROPS GENERAL THEORY FOR THE MODELLING OF LAMINATED MULTI-LAYER PLATES AND ITS APPLICATION FOR THE ANALYSIS OF SANDWICH PANELS WITH PLY DROPS Ole Thybo Thomsen Insttute of Mechancal Engneerng, Aalborg Unversty

More information

Chapter 3 Differentiation and Integration

Chapter 3 Differentiation and Integration MEE07 Computer Modelng Technques n Engneerng Chapter Derentaton and Integraton Reerence: An Introducton to Numercal Computatons, nd edton, S. yakowtz and F. zdarovsky, Mawell/Macmllan, 990. Derentaton

More information

ORIGIN 1. PTC_CE_BSD_3.2_us_mp.mcdx. Mathcad Enabled Content 2011 Knovel Corp.

ORIGIN 1. PTC_CE_BSD_3.2_us_mp.mcdx. Mathcad Enabled Content 2011 Knovel Corp. Clck to Vew Mathcad Document 2011 Knovel Corp. Buldng Structural Desgn. homas P. Magner, P.E. 2011 Parametrc echnology Corp. Chapter 3: Renforced Concrete Slabs and Beams 3.2 Renforced Concrete Beams -

More information

How Differential Equations Arise. Newton s Second Law of Motion

How Differential Equations Arise. Newton s Second Law of Motion page 1 CHAPTER 1 Frst-Order Dfferental Equatons Among all of the mathematcal dscplnes the theory of dfferental equatons s the most mportant. It furnshes the explanaton of all those elementary manfestatons

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Study on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component

Study on Non-Linear Dynamic Characteristic of Vehicle. Suspension Rubber Component Study on Non-Lnear Dynamc Characterstc of Vehcle Suspenson Rubber Component Zhan Wenzhang Ln Y Sh GuobaoJln Unversty of TechnologyChangchun, Chna Wang Lgong (MDI, Chna [Abstract] The dynamc characterstc

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

PHYS 705: Classical Mechanics. Newtonian Mechanics

PHYS 705: Classical Mechanics. Newtonian Mechanics 1 PHYS 705: Classcal Mechancs Newtonan Mechancs Quck Revew of Newtonan Mechancs Basc Descrpton: -An dealzed pont partcle or a system of pont partcles n an nertal reference frame [Rgd bodes (ch. 5 later)]

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

Randić Energy and Randić Estrada Index of a Graph

Randić Energy and Randić Estrada Index of a Graph EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL

More information

Preliminary Design of Moment-Resisting Frames

Preliminary Design of Moment-Resisting Frames Prelmnary Desgn of Moment-Resstng Frames Preprnt Aamer Haque Abstract A smple method s developed for prelmnary desgn of moment-resstng frames. Preprnt submtted to Elsever August 27, 2017 1. Introducton

More information

Mathematical Modeling and Numerical Simulation of Smart Structures Controlled by Piezoelectric Wafers and Fibers

Mathematical Modeling and Numerical Simulation of Smart Structures Controlled by Piezoelectric Wafers and Fibers Mathematcal Modelng and Numercal Smulaton of Smart Structures Controlled by Pezoelectrc Wafers and Fbers U. Gabbert #, W. Kreher*, H. Köppe # # Insttut für Mechank, Otto-von-Guercke Unverstät Magdeburg

More information

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions

Solving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions ISSN 746-7659 England UK Journal of Informaton and Computng Scence Vol. 9 No. 3 4 pp. 69-8 Solvng Fractonal Nonlnear Fredholm Integro-dfferental Equatons va Hybrd of Ratonalzed Haar Functons Yadollah Ordokhan

More information

Numerical analysis of buried pipes subjected to lateral soil movements

Numerical analysis of buried pipes subjected to lateral soil movements Numercal analss of bured ppes subjected to lateral sol movements P. Casamchele, M. Mauger & E. Motta Department of Cvl and Envronmental Engneerng, Catana Unverst, Ital Abstract he response of bured ppes

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

ksi a b 8.1 (a) B = stress = y / 0.1; A = strain = x / 10 Let, E(B A = a) = + a a i Nos. a i b i a i b i

ksi a b 8.1 (a) B = stress = y / 0.1; A = strain = x / 10 Let, E(B A = a) = + a a i Nos. a i b i a i b i 8. (a) B = stress = / 0.; A = stran = / 0 Let, E(B A = a) = + a Nos. a b a b a b 90-4 0 90-3 80-8 00 00-4 0 400-3 4000-8 400 3 80-4 30 840-3 7840-8 900 4 40-4 40 640-3 680-8 600 5 50-4 50 600-3 7040-8

More information

ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION

ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION VOL. 6, NO., NOVEMBER ISSN 89-668 6- Asan Research Publshng Network (ARPN). All rghts reserved. ANALYSIS OF TIMOSHENKO BEAM RESTING ON NONLINEAR COMPRESSIONAL AND FRICTIONAL WINKLER FOUNDATION Adel A.

More information

The Analytical Solution of a System of Nonlinear Differential Equations

The Analytical Solution of a System of Nonlinear Differential Equations Int. Journal of Math. Analyss, Vol. 1, 007, no. 10, 451-46 The Analytcal Soluton of a System of Nonlnear Dfferental Equatons Yunhu L a, Fazhan Geng b and Mnggen Cu b1 a Dept. of Math., Harbn Unversty Harbn,

More information

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model

Prof. Dr. I. Nasser Phys 630, T Aug-15 One_dimensional_Ising_Model EXACT OE-DIMESIOAL ISIG MODEL The one-dmensonal Isng model conssts of a chan of spns, each spn nteractng only wth ts two nearest neghbors. The smple Isng problem n one dmenson can be solved drectly n several

More information

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems

Application of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &

More information

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law: CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and

More information

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites

A comprehensive study: Boundary conditions for representative volume elements (RVE) of composites Insttute of Structural Mechancs A comprehensve study: Boundary condtons for representatve volume elements (RVE) of compostes Srhar Kurukur A techncal report on homogenzaton technques A comprehensve study:

More information

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00

ONE DIMENSIONAL TRIANGULAR FIN EXPERIMENT. Technical Advisor: Dr. D.C. Look, Jr. Version: 11/03/00 ONE IMENSIONAL TRIANGULAR FIN EXPERIMENT Techncal Advsor: r..c. Look, Jr. Verson: /3/ 7. GENERAL OJECTIVES a) To understand a one-dmensonal epermental appromaton. b) To understand the art of epermental

More information