Transactions of the VŠB Technical University of Ostrava, Mechanical Series. article No. 1907
|
|
- Tyrone Armstrong
- 5 years ago
- Views:
Transcription
1 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 PAT (THEOY) METODA KONEČNÝCH DIFEENCÍ POUŽITÁ PO NOSNÍKY NA PUŽNÉM PODKLADU ČÁST (TEOIE) Abstract Ths artcle s focused on the theory of straght and curved beams on elastc (Wnkler's) foundaton. For soluton of these problems of mechancs, the Fnte Dfference Method (.e. Central Dfference Method) can be appled. The basc nformaton about fnte dfferences and ther applcaton are explaned. Abstrakt Článek je zaměřen na teor přímých a křvých nosníků na pružném (Wnklerově) podkladu. Pro řešení těchto úloh mechanky, může být použta metoda konečných dferencí (tj. metoda centrálních dferencí). Základní nformace o konečných dferencích a jejch aplkacích jsou vysvětleny. INTODUCTION TO THE THEOY OF BEAMS ON ELASTIC FOUNDATION The basc analyss of bendng of beams on an elastc foundaton, see references [] to [], s developed on the assumpton that the strans are small. In ths context, an elastc foundaton s defned as a support whch s contnuously or dscontnuously dstrbuted along the length of the beam. The reacton Fg. Element of a straght beam on elastc foundaton. force q q x /Nm / dstrbuted n a foundaton s drectly proportonal to the deflecton v v(x) /m/ of a straght beam, see Fg., or proportonal to the radal dsplacement u u () /m/ of a curved beam (crcular arches), see Fg.. * M.Sc. Marek NIKODÝM Ph.D., VŠB - Techncal Unversty of Ostrava, Department of Mathematcs and Descrptve Geometry, 7. lstopadu /7, Ostrava, Czech epublc, Phone , E-mal marek.nkodym@vsb.cz ** Assocate Prof., MSc. Karel FYDÝŠEK Ph.D., ING-PAED IGIP, VŠB - Techncal Unversty of Ostrava, Faculty of Mechancal Engneerng, Department of Mechancs of Materals, 7. lstopadu, Ostrava, tel. (+0) 9 79, e-mal karel.frydrysek@vsb.cz 7
2 Fg. Example of a curved beam on elastc foundaton and ts element. Ths artcle s focused on the soluton of the straght and curved beams on elastc foundaton, see Fg. and, whch leads to the soluton of lnear dfferental equatons va Fnte Dfference Method (.e. Central Dfference Method). DIFFEENTIAL EQUATION FO STAIGHT BEAMS ON ELASTIC FOUNDATION The bendng of straght beams on elastc foundatons, see Fg., can be descrbed by ordnary lnear dfferental equaton: d v N d v d q q dm d q t d t t q () GA GA h where: E /Pa/ s modulus of elastcty of the beam, J /m / s the major prncpal second moment of area A /m / of the beam cross-secton, // s shear deflecton constant of the beam, G /Pa/ s shear modulus of the beam, N /N/ s normal force, q q /Nm / x s dstrbuted load (ntensty of force), m /N/ s dstrbuted couple (ntensty of moment), / deg / s coeffcent of thermal expanson of the beam, h /m/ s depth of the beam and t t / deg / s transversal temperature ncreasng n the beam. Equaton () s derved for the stuatons when nput parameters E, J, N,, G, A, and t h are constant. For more nformaton about the dervaton of eq. (), see references [] to [6]. From the Wnkler's theory, see references [] to [], t holds that: kv bkv () q t 8
3 /Pa/ where functons: k k x s stffness of the foundaton and K Kx /Nm / s modulus of the foundaton whch can be expressed as functons of varable x /m/ (.e. longtudnal changes n the foundaton) and b /m/ s wdth of the beam (see Fg. ). Hence, from eq. () and () follows: d v N d v d kv kv dm d q t d t t q () GA GA h In the most stuatons, the nfluences of shearng force, temperature and ntensty of moment can be neglected (or the beam s not exposed to them). Hence, from eq. () follows the smple form: d v N d v kv q. () DIFFEENTIAL EQUATION FO CUVED BEAMS ON ELASTIC FOUNDATION The bendng of curved beams on elastc foundatons, see Fg., can be descrbed by ordnary lnear dfferental equaton: d u d u dq u, () d d d where: /m/ s radus of the beam, /rad/ s angle varable and parameter // s gven by equaton: k. (6) q From the Wnkler's theory, see references [] to [], t s evdent that: ku bku (7) All others parameters mentoned n equatons () to (7) are explaned n former text. FINITE DIFFEENCES Let us consder an equdstant partton of the beam wth a step (nodes) "" along ts length, see Fg. (unloaded beam) and Fg. (loaded beam). /m/ and nodal ponts Fg. Solved straght beam s dvded nto nodes "". 9
4 Fg. Solved straght beam s dvded nto nodes "". Deflecton curve v = v(x) of a straght beam and ts dervatves are approxmated by polygon curves, see Fg.. Fg. Approxmaton of the deflectons by polygon curve and approxmaton of frst dervatve. Fnte dfferences can be defned as an approxmaton of dervatves. Hence, for the value of the frst dervatve, three types of dfferences can be defned accordng to Fg. : Backward dfference at the pont "": () dvx x v v v tan( ). (8) Forward dfference at the pont "": () dvx x v v v tan( ). (9) Central dfference at the pont "": () () v v () dv x x v v v tan( ). (0) In some references (for example [6]) are symbols "-", "+" noted as "-½" and "+½". Central dfferences (CD) are more accurate, therefore they wll be appled n the followng text. Smlarly, the hgher dervatves (at the pont "") can be approxmated by the central dfferences as: () d vx x v v v v, () () d vx x v v v v v, () 60
5 () d v x x v v 6v v v v, () () d vx x v v v v v v v, () 6 (6) d vx x v 6v v 0v v 6v v v. () 6 6 Smlarly, for a curved beams (.e. approxmatons for dervatves of functon u u () ), CD formulas can be derved by substtuton of varables v (for example u () du d u tan( ) () u () u CENTAL DIFFEENCE METHOD (CDM) FO STAIGHT BEAMS Accordng the Central Dfference Method (CDM), the dfferental equatons () for straght beams can be approxmated at the general pont "" (see eq. () and ()) as: N N k N q v v 6 v v v. (6) where k and q are the stffness of the foundaton and the dstrbuted loadng at the pont "", respectvely. Equaton (6) can be wrtten for all nodes 0,,,, n (.e. set of n+ lnear equatons followng from the dscretzaton of eq. ()), see Fg. 6. Ths set of equatons, together wth four dscretzed boundary condtons, lead to the soluton of system of n+ lnear equatons. There are always four fcttous nodes (-, - and n+, n+) outsde the ends of the beam, see Fg. 6. Hence, values of v at each node "" (.e. values of n+ deflectons) can be receved, see also reference [] and [7]. u etc.). Fg. 6 Example solved n reference [] (straght beam on elastc foundaton loaded by couple M). Note: Theoretcally, f step 0 (.e. n ) then the numercal soluton converge to exact soluton. Fg. 7 and 8 show the sparse matrces arsng from CDM for beam on elastc foundaton loaded by couple M, see Fg. 6. Ths example s solved n reference []. 6
6 Fg. 7 Example solved n reference [] (sparsty patterns of matrces n CDM, number of elements n = ). Fg. 8 Example solved n reference [] (sparsty patterns of matrces n CDM, number of elements n = 0). 6
7 6 CENTAL DIFFEENCE METHOD (CDM) FO CUVED BEAMS Accordng the CDM, the dfferental equatons () for curved beams can be approxmated at the general pont "" (see modfed eq. () and ()) as: u u u. (7) () u u q, u where k, k and q () dq are stffness of the foundaton, parameter and d frst dervatve of dstrbuted load at the pont "", respectvely. Equaton (7) can be wrtten for all nodes 0,,,, n (.e. set of n+ lnear equatons followng from the dscretzaton of eq. ()). Ths set of equatons, together wth fve dscretzed boundary condtons, lead to the soluton of system of n+7 lnear equatons. Hence, values of u at each node "" (.e. values of n+7 deflectons) can be receved, see also reference [] and [7]. Note: Theoretcally, f step 0 (.e. n ) then numercal soluton converge to exact soluton. CONCLUSION Ths artcle shows dervatons and way of applcaton of the Central Dfference Method (CDM) as a numercal method sutable for the soluton of the straght or curved beams on elastc foundaton. For more nformaton about applcatons of CDM, see [], [6], [7] and [8]. CDM seems to be a good alternatve to wdely spread Fnte Element Method. Another ways of the solutons and applcatons of structures on elastc foundaton are presented n [] to [9]. ACKNOWLEDGEMENT Ths work has been supported by the Czech-Slovak project 7AMBSK and Slovak-Czech project SK-CZ EFEENCES [] FYDÝŠEK, K.: Beams and Frames on Elastc Foundaton (Nosníky a rámy na pružném podkladu ), monograph, Faculty of Mechancal Engneerng, VŠB - Techncal Unversty of Ostrava, ISBN , Ostrava, Czech epublc, 006, pp.6. [] FYDÝŠEK, K., JANČO,. et al: Beams and Frames on Elastc Foundaton (Nosníky a rámy na pružném podkladu ), monograph,všb - Techncal Unversty of Ostrava, ISBN , Ostrava, Czech epublc, 008, pp.6. [] FYDÝŠEK, K., NIKODÝM, M. et al: Beams and Frames on Elastc Foundaton (Nosníky a rámy na pružném podkladu ), monograph, VŠB - Techncal Unversty of Ostrava, ISBN , Ostrava, Czech epublc, 00. [] HETÉNYI, M.: Beams on Elastc Foundaton, Ann Arbor, Unversty of Mchgan Studes, USA, 96. [] JANČO,.: Numercal Methods of Soluton of Beam on Elastc Foundaton, In: th Internatonal Conference MECHANICAL ENGINEEING 00, Faculty of Mechancal Engneerng, Slovak Unversty of Technology n Bratslava, ISBN , Bratslava, Slovaka, 00, pp. S-60 S-6. [6] JONES, G.: Analyss of Beams on Elastc Foundatons Usng Fnte Dfference Theory, ISBN , Thomas Telford Publshng, London, UK, 997, pp.6. [7] NIKODÝM, M., FYDÝŠEK, K.: Fnte Dfference Method Used for the Beams on Elastc Foundaton Part (Applcatons), Transactons of the VŠB Techncal Unversty of Ostrava, Mechancal Seres, vol. LVIII, 0 (n ths journal, n prnt). 6
8 [8] KAMIŃSKI, M.: A generalzed verson of the perturbaton-based stochastc fnte dfference method for elastc beams, Journal of Theoretcal and Appled Mechancs, 009, 7,, pp [9] TVDÁ, K.; DICKÝ, J.: Comparson of Optmzaton Methods. In: Proceedngs of the -th Internatonal Conference on New Trends n Statcs and Dynamcs of Buldngs, Bratslava, Slovaka, 00, pp. 6-6, ISBN
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 informationSTATIC ANALYSIS OF TWO-LAYERED PIEZOELECTRIC BEAMS WITH IMPERFECT SHEAR CONNECTION
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
More informationDUE: 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 informationNumerical 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 informationLecture 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 informationA 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 informationChapter 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 informationChapter 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 informationFINITE 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 informationTransactions of the VŠB Technical University of Ostrava, Mechanical Series. article No Roland JANČO *
Transactions of the VŠB Technical University of Ostrava, Mechanical Series No. 1, 013, vol. LIX article No. 1930 Roland JANČO * NUMERICAL AND EXACT SOLUTION OF BUCKLING LOAD FOR BEAM ON ELASTIC FOUNDATION
More informationANALYSIS 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 informationIndeterminate 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 informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationFinite 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 informationSecond 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 informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationFUZZY 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 informationTHE 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 informationIrregular 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 informationLAB 4: Modulus of elasticity
LAB 4: Modulus of elastcty 1. Preparaton: modulus of elastcty (chapter15, p.79) Hook s law graphcal determnaton of modulus of elastcty (p.8) determnaton of modulus of elastcty n tenson and flexural stress
More informationLecture 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 informationFrame 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 informationNovember 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project
SE 8 Fnal Project Story Shear Frame u m Gven: u m L L m L L EI ω ω Solve for m Story Bendng Beam u u m L m L Gven: m L L EI ω ω Solve for m 3 3 Story Shear Frame u 3 m 3 Gven: L 3 m m L L L 3 EI ω ω ω
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More information2 LATERAL TORSIONAL BUCKLING OF AN IDEAL BEAM WITHOUT LATERAL RESTRAINTS
0.55/tvsb-05-000 Transactons of the VŠB Techncal Unversty of Ostrava No., 05, Vol. 5, Cvl Engneerng Seres paper #0 Ivan BALÁZS, Jndřch MELCHER STABILITY OF THIN-WALLED BEAMS WITH LATERAL CONTINUOUS RESTRAINT
More informationMODELLING 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 informationNONLINEAR 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 informationThe 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 information829. 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 informationPlease initial the statement below to show that you have read it
EN0: Structural nalyss Exam I Wednesday, March 2, 2005 Dvson of Engneerng rown Unversty NME: General Instructons No collaboraton of any nd s permtted on ths examnaton. You may consult your own wrtten lecture
More informationOne-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 informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationAPPENDIX 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 informationEVALUATION 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 informationProfessor 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 informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationOFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES
ICAMS 204 5 th Internatonal Conference on Advanced Materals and Systems OFF-AXIS MECHANICAL PROPERTIES OF FRP COMPOSITES VLAD LUPĂŞTEANU, NICOLAE ŢĂRANU, RALUCA HOHAN, PAUL CIOBANU Gh. Asach Techncal Unversty
More informationApplication 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 informationCHAPTER 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 informationSolution Methods for Beam and Frames on Elastic Foundation Using the Finite Element Method
Solution Methods for Beam and Frames on Elastic Foundation Using the Finite Element Method Spôsoby riešenie nosníkov a rámov na pružnom podklade pomocou metódy konečných prvkov Roland JANČO 1 Abstract:
More informationOne 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 information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationDETERMINATION 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 informationSolution for singularly perturbed problems via cubic spline in tension
ISSN 76-769 England UK Journal of Informaton and Computng Scence Vol. No. 06 pp.6-69 Soluton for sngularly perturbed problems va cubc splne n tenson K. Aruna A. S. V. Rav Kant Flud Dynamcs Dvson Scool
More informationGEO-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 informationIn 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 informationELASTIC 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 informationComparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method
Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method
More informationThe 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 informationHEAT TRANSFER THROUGH ANNULAR COMPOSITE FINS
Journal of Mechancal Engneerng and Technology (JMET) Volume 4, Issue 1, Jan-June 2016, pp. 01-10, Artcle ID: JMET_04_01_001 Avalable onlne at http://www.aeme.com/jmet/ssues.asp?jtype=jmet&vtype=4&itype=1
More informationAPPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS
6th ICPT, Sapporo, Japan, July 008 APPROXIMATE ANALYSIS OF RIGID PLATE LOADING ON ELASTIC MULTI-LAYERED SYSTEMS James MAINA Prncpal Researcher, Transport and Infrastructure Engneerng, CSIR Bult Envronment
More informationNUMERICAL 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 informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationImplicit Integration Henyey Method
Implct Integraton Henyey Method In realstc stellar evoluton codes nstead of a drect ntegraton usng for example the Runge-Kutta method one employs an teratve mplct technque. Ths s because the structure
More informationApplication 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 informationParameter Estimation for Dynamic System using Unscented Kalman filter
Parameter Estmaton for Dynamc System usng Unscented Kalman flter Jhoon Seung 1,a, Amr Atya F. 2,b, Alexander G.Parlos 3,c, and Klto Chong 1,4,d* 1 Dvson of Electroncs Engneerng, Chonbuk Natonal Unversty,
More informationExponential Type Product Estimator for Finite Population Mean with Information on Auxiliary Attribute
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 10, Issue 1 (June 015), pp. 106-113 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) Exponental Tpe Product Estmator
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationSolving 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 informationNumerical Methods. ME Mechanical Lab I. Mechanical Engineering ME Lab I
5 9 Mechancal Engneerng -.30 ME Lab I ME.30 Mechancal Lab I Numercal Methods Volt Sne Seres.5 0.5 SIN(X) 0 3 7 5 9 33 37 4 45 49 53 57 6 65 69 73 77 8 85 89 93 97 0-0.5 Normalzed Squared Functon - 0.07
More informationStatistical Energy Analysis for High Frequency Acoustic Analysis with LS-DYNA
14 th Internatonal Users Conference Sesson: ALE-FSI Statstcal Energy Analyss for Hgh Frequency Acoustc Analyss wth Zhe Cu 1, Yun Huang 1, Mhamed Soul 2, Tayeb Zeguar 3 1 Lvermore Software Technology Corporaton
More informationInductance 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 informationA 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 informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationVIBRATION ANALYSIS OF PRE-TWISTED BEAMS USING THE SPLINE COLLOCATION METHOD
16 Journal of Marne Scence and echnology, Vol. 17, o., pp. 16-11 (9 VIBRAIO AALYSIS OF PRE-WISED BEAMS USIG HE SPLIE COLLOCAIO MEHOD Mng-Hung Hsu* Key words: natural frequency, pre-twsted beam, splne collocaton
More information16 Circular Footing on a Semi-infinite Elastic Medium
Crcular Footng on a Sem-nfnte Elastc Medum 16-1 16 Crcular Footng on a Sem-nfnte Elastc Medum 16.1 Problem Statement Ths verfcaton problem nvolves a rgd crcular footng restng on an elastc half-space. In
More informationUncertainty 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 informationBuckling 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 informationCubic 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 informationGeneral displacement arch-cantilever element method for stress analysis of arch dam
Water Scence and Engneerng, 009, (): 3-4 do:0.388/j.ssn.674-370.009.0.004 http://kkb.hhu.edu.cn e-mal: wse@hhu.edu.cn General dsplacement arch-cantlever element method for stress analyss of arch dam Hao
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationMECHANICS OF MATERIALS
Fourth Edton CHTER MECHNICS OF MTERIS Ferdnand. Beer E. Russell Johnston, Jr. John T. DeWolf ecture Notes: J. Walt Oler Texas Tech Unversty Stress and Stran xal oadng Contents Stress & Stran: xal oadng
More informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationModule 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 informationElectrical 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 informationAPPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS
Unversty of Oulu Student Laboratory n Physcs Laboratory Exercses n Physcs 1 1 APPEDIX FITTIG A STRAIGHT LIE TO OBSERVATIOS In the physcal measurements we often make a seres of measurements of the dependent
More informationSolutions of Beams, Frames and 3D Structures on Elastic Foundation Using FEM
INTERNATIONA JOURNA OF MECHANICS Solutions of Beams, Frames and D Structures on Elastic Foundation Using FEM Karel FRYDRÝŠEK, Roland JANČO, Horst GONDEK Abstract This paper contains numerical methods and
More informationActual behaviour of composite externally CFRP-reinforced timber beams stress analysis
Actual behavour of composte externally CFRP-renforced tmber beams stress analyss Marcela Karmazínová Abstract The paper s focused on the problems of composte CFRP-renforced tmber beams. CFRP renforcement
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationLecture 8 Modal Analysis
Lecture 8 Modal Analyss 16.0 Release Introducton to ANSYS Mechancal 1 2015 ANSYS, Inc. February 27, 2015 Chapter Overvew In ths chapter free vbraton as well as pre-stressed vbraton analyses n Mechancal
More informationSimulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests
Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationTHE STURM-LIOUVILLE EIGENVALUE PROBLEM - A NUMERICAL SOLUTION USING THE CONTROL VOLUME METHOD
Journal of Appled Mathematcs and Computatonal Mechancs 06, 5(), 7-36 www.amcm.pcz.pl p-iss 99-9965 DOI: 0.75/jamcm.06..4 e-iss 353-0588 THE STURM-LIOUVILLE EIGEVALUE PROBLEM - A UMERICAL SOLUTIO USIG THE
More informationAsymptotics of the Solution of a Boundary Value. Problem for One-Characteristic Differential. Equation Degenerating into a Parabolic Equation
Nonl. Analyss and Dfferental Equatons, ol., 4, no., 5 - HIKARI Ltd, www.m-har.com http://dx.do.org/.988/nade.4.456 Asymptotcs of the Soluton of a Boundary alue Problem for One-Characterstc Dfferental Equaton
More informationORIGIN 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 informationANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)
Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of
More informationUNIVERSITY 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 informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationThe binomial transforms of the generalized (s, t )-Jacobsthal matrix sequence
Int. J. Adv. Appl. Math. and Mech. 6(3 (2019 14 20 (ISSN: 2347-2529 Journal homepage: www.jaamm.com IJAAMM Internatonal Journal of Advances n Appled Mathematcs and Mechancs The bnomal transforms of the
More informationLecture 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 informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationNumerical Solutions of a Generalized Nth Order Boundary Value Problems Using Power Series Approximation Method
Appled Mathematcs, 6, 7, 5-4 Publshed Onlne Jul 6 n ScRes. http://www.scrp.org/journal/am http://.do.org/.436/am.6.77 umercal Solutons of a Generalzed th Order Boundar Value Problems Usng Power Seres Approxmaton
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationA 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 informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More information