FINITE DIFFERENCE SOLUTION OF MIXED BOUNDARY-VALUE ELASTIC PROBLEMS
|
|
- Wilfred Scott
- 5 years ago
- Views:
Transcription
1 4 t Internatonal Conference on Mecancal Engneerng, December 6-8, 1, Daa, Banglades/pp. V FINITE DIFFERENCE SOLUTION OF MIXED BOUNDARY-VALUE ELASTIC PROBLEMS S. Reaz Amed, Noor Al Quddus and M. Waa Uddn Department of Mecancal Engneerng, Banglades Unversty of Engneerng and Tecnology, Daa 1, Banglades Abstract A new numercal metod of soluton as been developed for te analyss of deformaton and stresses n elastc bodes subected to mxed -condtons. Te program s capable of dealng wt bot regular and rregular sapes of boundares approprately. An deal matematcal model, based on te dsplacement potental functon, as been used n te fnte dfference soluton. Te present paper demonstrates te applcaton of te newly developed computatonal sceme to a wdely used body wt curved boundares, namely, nvolute profle spur gear teet. Keyword: arbtrary-saped elastc body; dsplacement potental functon; fnte-dfference tecnque. INTRODUCTION Te rapd growt n te use of computers n te past decade gave rse to te developments of advanced computatonal sceme. But someow te elastcty problems are stll sufferng from a lot of sortcomngs. Actually, te management of condtons and sapes as remaned as te bggest urdle n te soluton process of te problems of sold mecancs. Te necessty of te management of sape as lead to te nventon of te fnte element tecnque wt ts overwelmng popularty. Te oter factor of mpedment to qualty soluton of elastc problems s te treatment of te transton n condtons. Te dffcultes nvolved n tryng to solve practcal stress problems usng te exstng approaces, for example, te stress functon approac [1] and te two dsplacement functons approac [], are clearly ponted out n our prevous reports [-6] and also by Durell [7]. In an attempt to crcumvent te dffcultes, Dow et al [8] ntroduced a new modelng approac for te fnte dfference applcatons n sold mecancs. Tey reported tat te accuracy of te fnte dfference metod n reproducng te state of stresses along te was muc ger tan tat of fnte element metod of analyss. Te dsplacement potental functon formulaton of two dmensonal elastc problems used ere was frst ntroduced by Uddn []; later Idrs [] used t for obtanng analytcal solutons of mxed value elastc problems and Amed [4-6] developed te computatonal sceme to extend ts use n obtanng te numercal solutons of a number of mxed value problems. Te ratonalty and relablty of te formulaton s ceced repeatedly by comparng te results of mxed value elastc problems obtaned troug ts formulaton wt tose avalable n te lterature. Recently, our wor as been extended to nclude te problems of arbtrary sapes[1]. Te present paper demonstrates te applcaton of te newly developed computatonal sceme to a wdely used body wt curved boundares, namely, nvolute profle spur gear teet. FORMULATION OF THE PROBLEM Wt reference to a rectangular co-ordnate system, te dfferental equatons of equlbrum for twodmensonal problems n terms of dsplacement components are [1] u 1 u u + + = (1) u 1 u u + + = () Tese two omogeneous ellptc partal-dfferental equatons, wt approprate condtons, ave to be solved for te case of a two-dmensonal problem wen te body forces are assumed to be absent. In te present approac, te problem s reduced to te determnaton of a sngle functon ψ [] nstead of fndng te two varables u and v, smultaneously, satsfyng te equlbrum equatons (1) and (). Te potental functon ψ(x,y) s defned n terms of dsplacement components as u = 1 v = ( 1 ) + () Secton V: Appled Mecancs 171
2 ICME 1, Daa, December 6-8 Wen te dsplacement components n equaton (1) and () are replaced by ψ(x,y), equaton (1) becomes an dentty and te only condton tat as to be satsfed becomes = (4) 4 4 Terefore, te problem s formulated n suc a way tat a sngle functon ψ as to be evaluated from te barmonc equaton (4), satsfyng te condtons tat are specfed at te. Now, te condtons, at any pont on an arbtrary-saped, are nown n terms of te normal and tangental components of dsplacement, u n and u t and of stresses, σ n and σ t. Tese four components are frst expressed n terms of u, v, σ x, σ y, σ xy [6] te components of dsplacements and stress wt respect to te reference axes x and y and fnally, n terms of te functon ψ, as follows: m m ( ) ( 1 ) u n x, y = + l (5) σ u t σ n ( x y) ( 1 ) l l, = m (6) ( x, y) ( x, y) = = E ( ) lm + ( l m m ) ( l + m ) ] + lm (7) E ( ) ( l m ) lm( + ) t ( l m ) + lm( ) ] + 1 (8) y Te computatonal wor n solvng any problem remans te same n te present case as t was n te case of φ-formulaton, snce bot of tem ave to satsfy te same b-armonc equaton. But te ψ-formulaton s free from te nablty of te φ-formulaton n andlng te mxed condton. SOLUTION PROCEDURE Fnte-dfference tecnque s used to dscretze te governng b-armonc equaton and also te dfferental equatons assocated wt te condtons. Te dscrete values of ψ(x,y) at te mes ponts of te doman concerned (see Fg. 1), s solved from te system of lnear algebrac equatons resultng from te dscretzaton of te governng equaton and te assocated condtons. Te present sceme nvolves evaluaton of te functon ψ at te nodal ponts of rectangular grds over a geometrcally rregular regon, were te may not always pass troug te rectangular mes ponts, as sown n Fg.1. Te governng b-armonc equaton, wc s used to evaluate te functon ψ only x/ Pyscal σ n Involute toot Root secton Gear stoc Reference False 1 y/ Fg.1 Dscretzaton sceme for te gear toot doman. (m = 1 mm, a = m, b = 1.5 m, θ p =, d p =5 mm, load σ n /E = e 4). at te nternal mes ponts, s expressed n ts correspondng dfference equaton usng central dfference operators. Te complete fnte-dfference expresson of te b-armonc equaton (4) s gven by R 4 {ψ(-,)+ψ(+,)}-4r (R ){ψ(-1,)+ψ(+1,)} -4(R ){ψ(,+1)+ψ(,-1)}+(6r 4 +8R +6)ψ(,) +R {ψ(-1,-1)+ψ(-1,+1)+ψ(+1,-1) +ψ(+1,+1)}+ψ(,-)+ψ(,+)= (9) Te correspondng grd structure of te governng equaton for an arbtrary nternal mes pont (,) s sown n Fgure. It s tus clear tat wen te pont of applcaton (,) becomes an mmedate negbour of te pyscal, te equaton wll nvolve mes ponts bot nteror and exteror to te pyscal Fg. Grd structure for te governng b-armonc equaton. An magnary (false), exteror to te reference of te doman s ntroduced as sown n te Fg.1. As te dfferental equatons assocated wt te condtons contan second Secton V: Appled Mecancs 17
3 ICME 1, Daa, December 6-8 and trd-order dervatves of te functon ψ, te use of central dfference expressons ultmately leads to te ncluson of ponts exteror to te false. Te dervatves n te expressons of -condtons are tus replaced by ter bacward- or forwarddfference formulae, eepng te order of te local truncaton error te same. Any -value at a -pont, not matcng wt te feld grd ponts, s replaced by te lnear nterpolaton of ts value at te two or four negbourng grd ponts, one of wc s desgnated as te reference pont, wc s unque to te pont. Te dfference equatons of te condtons for any arbtrary pont on te, not matcng wt te feld grd ponts, can readly be obtaned by nowng te actual poston of te pont wt respect to ts reference grd pont. Referrng to Fg., te dfference expressons for condtons at any arbtrary pont P, te locaton of wc s defned by Fg. Locators c and s of te pont P wt respect to ts reference pont R1. Pyscal - -1 c and s, s obtaned by usng lnear nterpolaton of te expressons of te -values at te ponts R1, R, R and R4. Te fnal form of te grd structure so obtaned for te stress components for ponts not matcng wt te feld grd ponts s llustrated n Fg.4. Snce tere are two condtons to be satsfed at an arbtrary pont on te pyscal of te elastc body, two dfference equatons are assgned to a sngle pont on te. Out of tese two equatons, one s used to evaluate te functon ψ at te reference pont correspondng to te pyscal pont and te remanng one for te correspondng pont on te false. Te dscrete values of te potental functon, ψ (x,y), at mes ponts are solved from te system of algebrac equatons resultng from te dscretzaton of te governng equaton and te assocated condtons, by te use of drect metod of soluton (Colesey s trangular decomposton metod). Fnally, te same dfference equatons are organzed for te evaluaton of all te parameters of nterest n te soluton of te body at every nteror as well as ponts from te ψ values at mes ponts of te doman. RESULTS AND DISCUSSIONS R1 R c s P +1 R4 R Pyscal Gear toot as been taen as te arbtrary-saped elastc body, made of ordnary steel (=., E=9 GPa), to obtan te numercal values of stresses n t as well as te deformed sape of t under load. Te conugate loadng on te toot s approxmated ere by te dstrbuton of normal load over a small regon at te contact pont. It sould be mentoned ere tat, because of te Sant-Venant s prncple, te dstrbuton of loadng at te contact pont does not affect te stresses n te gear at te crtcal regon of te fllets or at any oter crtcal secton of te toot, far away from te loadng, so long te total loadng remans te same. Obvously, te stresses at and near te contact pont s very muc dependent on te approxmaton n te assumed dstrbuton of loadng at te contact pont. 4e-4 e-4 x/ = 17 x/ = 18 x/ = 19 x/ = σ x / E +5 -e Fg. 4 Grd structure for σ n or σ t at top-left ( pont not matcng wt feld grd pont) e Grd poston (y/) Fg. 5 Radal stresses (σ x /E ) at sectons x/ =17,18,19 and 1 of te gear toot Secton V: Appled Mecancs 17
4 ICME 1, Daa, December 6-8 Te stresses at dfferent sectons of te tp-loaded spur gear toot, sown n Fg.1 (referred to as type-a ere), are presented n Fgs 5-6. Overall, te magntudes of te stresses σ y and σ xy are smaller tan tat of te radal stress σ x. Maxmum stresses are found to occur around te root secton of te toot, and te fllet zone s te most crtcal secton n terms of stresses. Te magntude of stresses ust below te root secton s lower and gradually decreases towards te center of te gear but does not become zero on te bottom, wc s consdered rgdly fxed. Te load appled normally to te toot face near te tp as bascally two components, one of bendng, wc s compressng te load-free sde and tensonng te loaded sde of te toot, wle te oter s compressng over te wole secton. Tus, te combned effect of bendng and compresson gves lower value of tenson n te load-sde fllet and ger value of compresson n te load free-sde fllet (Fg.5). σ xy / E 1e-4-1e-4 -e-4 -e-4 x/ = 17 x/ = 18 x/ = 19 x/ = 1-4e Grd poston (y/) Fg. 6 Searng stresses (σ xy /E ) at sectons x/ = 17,18,19 and 1 of te gear toot Te deformed toot profle of te gear wt respect to ts ntal undeformed profle under te loadng at dfferent regons are llustrated n Fg.7. Te profles sow tat maxmum dsplacement occurs under te tp loadng. Te toot deflects n te drecton of load and te deflecton s maxmum at te tp and gradually decreases towards te root. Anoter spur gear toot (profle: nvolute, m = 1 mm, θ p =, a = m, b = m p, d p = mm ) s consdered ere for comparng te present ψ-soluton wt te publsed FEM results[9]. Te values of radal stress at te root secton, due to te tp, ptc and root loadng are presented n Table-1 along wt te FEM results. Te magntudes of te stresses at te extreme end of te root secton, especally at te tenson sde, are found to be n good agreement. However, te FEM soluton [5] predcts lower value of te compressve stress at te bottom surface of te root secton tan tat Table 1 Comparson of present soluton for stresses wt tat of FEM soluton Load Radal stress at te Root secton (MPa) Loaded surface Load-free surface FDM FEM FDM FEM Tp Ptc Root at ponts above t. Ts s gly unlely, as te bendng stress s te gest at te bottom, wc s furter enanced due to te effect of stress concentraton at tat pont. In our present approac, a sngle varable s evaluated at eac pont nstead of solvng for two varables smultaneously as n te case of FEM, wc reduces our computatonal wor drastcally. It s noted tat te tme taen to solve a gear toot (9 X 7 mes) on a personal computer (45 MHz) was found to be.91 s. Te accuracy as well as te relablty of te present metod as been verfed repeatedly by comparng te results wt avalable solutons [4,1]. CONCLUSIONS Te dffcultes of managng te sapes n fnte dfference tecnques of analyss for wc te fnte-element metod of soluton of elastc problems was nvented wt a manfold ncrease of computatonal wors and a lot of loss of sopstcaton s overcome ere by a novel computatonal tecnque n te management of sapes n te fnte-dfference metod of soluton. Bot te qualtatve and quanttatve results of te spur gear toot, and ter comparson wt tose of FEM soluton establs te relablty and sutablty of te present numercal model. NOMENCLATURE load-free ptc-loadng tp-loadng Fg.7 Toot profles under root loadng (dsplacements are 1 tmes enlarged) E elastc modulus of te materal Posson s rato u,v dsplacement components n te x- and y- drectons σ x, σ y normal stress components n te x- and y- drectons φ Ary s stress functon Secton V: Appled Mecancs 174
5 ICME 1, Daa, December 6-8 ψ dsplacement potental functon u n, u t normal and tangental components of dsplacement on te σ n, σ t normal and tangental components of stress on te l, m drecton cosnes of te normal at any pont on te, mes lengts n te x- and y- drectons R / a, b addendum and dedendum of te gear m,,d p,θ p module, ptc dameter, pressure angle of gear toot c, s locators n y- and x- drectons, wt respect to te reference grd pont 1. Tmoseno S, Gooder JN. Teory of Elastcty. ( rd edn.). McGraw-Hll, New Yor, N. Y. (1979).. Uddn, M.W., "Fnte Dfference Soluton of Two- Dmensonal Elastc Problems wt Mxed Boundary Condtons", M. Sc. Tess, Carleton Unversty, Canada, (1966).. Idrs, A.B.M., Amed, S.R. and Uddn, M.W., "Analytcal Soluton of a -D Elastc Problem wt Mxed Boundary Condtons", Journal of te Insttuton of Engneers (Sngapore), Vol. 6, No. 6 (1996), pp Amed, S.R., Kan, M.R., Islam, K.M.S. and Uddn, M.W., Analyss of Stresses n Deep Beams Usng Dsplacement Potental Functon Journal of Insttuton of Engneers (Inda), Vol. 77, (1996), pp Amed, S.R., Idrs, A.B.M. and Uddn, M.W., "Numercal Soluton of Bot Ends Fxed Deep Beams", Computers & Structures, Vol. 61, No. 1 (1996), pp Amed, S.R., Kan, M.R., Islam, K.M.S., and Uddn, M.W., "Investgaton of stresses at te fxed end of deep cantlever beams", Computers & Structures, Vol. 69, (1998), pp Durell, A.J. and Ranganayaamma, B., "Parametrc soluton of stresses n beams", Journal of Engneerng Mecancs, Vol. 115, No. (1989), pp Dow, J.O., Jones, M.S. and Harwood, S.A., "A new approac to modelng for fnte dfference applcatons n sold mecancs", Internatonal Journal for Numercal Metods n Engneerng, Vol., (199), pp Ramamurt, V. & Rao, M.A., Dynamc Analyss of Spur Gear Teet, Computers and Structures, Vol. 9, No. 5, (1988), pp Aanda, M.A.S., Amed, S.R., Kan, M.R., and Uddn, M.W., "A fnte-dfference sceme for mxed -value problems of arbtrarysaped elastc bodes", Advances n Engneerng Software, Vol. 1, No. (), pp REFERENCES Secton V: Appled Mecancs 175
Numerical Simulation of One-Dimensional Wave Equation by Non-Polynomial Quintic Spline
IOSR Journal of Matematcs (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 14, Issue 6 Ver. I (Nov - Dec 018), PP 6-30 www.osrournals.org Numercal Smulaton of One-Dmensonal Wave Equaton by Non-Polynomal
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 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 informationTR/28. OCTOBER CUBIC SPLINE INTERPOLATION OF HARMONIC FUNCTIONS BY N. PAPAMICHAEL and J.R. WHITEMAN.
TR/8. OCTOBER 97. CUBIC SPLINE INTERPOLATION OF HARMONIC FUNCTIONS BY N. PAPAMICHAEL and J.R. WHITEMAN. W960748 ABSTRACT It s sown tat for te two dmensonal Laplace equaton a unvarate cubc splne approxmaton
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 informationThe Finite Element Method: A Short Introduction
Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental
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 informationA Spline based computational simulations for solving selfadjoint singularly perturbed two-point boundary value problems
ISSN 746-769 England UK Journal of Informaton and Computng Scence Vol. 7 No. 4 pp. 33-34 A Splne based computatonal smulatons for solvng selfadjont sngularly perturbed two-pont boundary value problems
More informationCOMPOSITE 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 information5 The Laplace Equation in a convex polygon
5 Te Laplace Equaton n a convex polygon Te most mportant ellptc PDEs are te Laplace, te modfed Helmoltz and te Helmoltz equatons. Te Laplace equaton s u xx + u yy =. (5.) Te real and magnary parts of an
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 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 informationSolving Singularly Perturbed Differential Difference Equations via Fitted Method
Avalable at ttp://pvamu.edu/aam Appl. Appl. Mat. ISSN: 193-9466 Vol. 8, Issue 1 (June 013), pp. 318-33 Applcatons and Appled Matematcs: An Internatonal Journal (AAM) Solvng Sngularly Perturbed Dfferental
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 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 informationConjunction of Displacement Fields of the Element Free Galerkin Method and Finite Element Method
amkang Journal of Scence and Engneerng, Vol. 10, No 1, pp. 4150 (2007) 41 Conjuncton of Dsplacement Felds of te Element Free Galerkn Metod and Fnte Element Metod Cen-Hsun Ln* and Can-Png Pan Department
More informationMultivariate Ratio Estimator of the Population Total under Stratified Random Sampling
Open Journal of Statstcs, 0,, 300-304 ttp://dx.do.org/0.436/ojs.0.3036 Publsed Onlne July 0 (ttp://www.scrp.org/journal/ojs) Multvarate Rato Estmator of te Populaton Total under Stratfed Random Samplng
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 informationTR/95 February Splines G. H. BEHFOROOZ* & N. PAPAMICHAEL
TR/9 February 980 End Condtons for Interpolatory Quntc Splnes by G. H. BEHFOROOZ* & N. PAPAMICHAEL *Present address: Dept of Matematcs Unversty of Tabrz Tabrz Iran. W9609 A B S T R A C T Accurate end condtons
More informationMultigrid Methods and Applications in CFD
Multgrd Metods and Applcatons n CFD Mcael Wurst 0 May 009 Contents Introducton Typcal desgn of CFD solvers 3 Basc metods and ter propertes for solvng lnear systems of equatons 4 Geometrc Multgrd 3 5 Algebrac
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 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 informationStanford University CS254: Computational Complexity Notes 7 Luca Trevisan January 29, Notes for Lecture 7
Stanford Unversty CS54: Computatonal Complexty Notes 7 Luca Trevsan January 9, 014 Notes for Lecture 7 1 Approxmate Countng wt an N oracle We complete te proof of te followng result: Teorem 1 For every
More informationMaejo International Journal of Science and Technology
Maejo Int. J. Sc. Technol. () - Full Paper Maejo Internatonal Journal of Scence and Technology ISSN - Avalable onlne at www.mjst.mju.ac.th Fourth-order method for sngularly perturbed sngular boundary value
More informationOn Pfaff s solution of the Pfaff problem
Zur Pfaff scen Lösung des Pfaff scen Probles Mat. Ann. 7 (880) 53-530. On Pfaff s soluton of te Pfaff proble By A. MAYER n Lepzg Translated by D. H. Delpenc Te way tat Pfaff adopted for te ntegraton of
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 informationA New Recursive Method for Solving State Equations Using Taylor Series
I J E E E C Internatonal Journal of Electrcal, Electroncs ISSN No. (Onlne) : 77-66 and Computer Engneerng 1(): -7(01) Specal Edton for Best Papers of Mcael Faraday IET Inda Summt-01, MFIIS-1 A New Recursve
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 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 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 informationDirect Methods for Solving Macromolecular Structures Ed. by S. Fortier Kluwer Academic Publishes, The Netherlands, 1998, pp
Drect Metods for Solvng Macromolecular Structures Ed. by S. Forter Kluwer Academc Publses, Te Neterlands, 998, pp. 79-85. SAYRE EQUATION, TANGENT FORMULA AND SAYTAN FAN HAI-FU Insttute of Pyscs, Cnese
More informationA 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 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 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 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 informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More informationInitial-value Technique For Singularly Perturbed Two Point Boundary Value Problems Via Cubic Spline
Unversty of Central Florda Electronc Teses and Dssertatons Masters Tess (Open Access) Intal-value Tecnque For Sngularly Perturbed Two Pont Boundary Value Problems Va Cubc Splne 010 Lus G. Negron Unversty
More informationORDINARY DIFFERENTIAL EQUATIONS EULER S METHOD
Numercal Analss or Engneers German Jordanan Unverst ORDINARY DIFFERENTIAL EQUATIONS We wll eplore several metods o solvng rst order ordnar derental equatons (ODEs and we wll sow ow tese metods can be appled
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 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 informationDynamic Design of Thick Orthotropic Cantilever. Plates with Consideration of Bimoments.
World Journal of Mecancs, 06, 6, 4-56 ttp://www.scrp.org/journal/wjm ISSN Onlne: 60-050 ISSN Prnt: 60-049X Dynamc Desgn of Tc Ortotropc Cantlever Plates wt Consderaton of Bmoments Мaamatal K. Usarov Insttute
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 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 informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
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 informationThe finite element method explicit scheme for a solution of one problem of surface and ground water combined movement
IOP Conference Seres: Materals Scence and Engneerng PAPER OPEN ACCESS e fnte element metod explct sceme for a soluton of one problem of surface and ground water combned movement o cte ts artcle: L L Glazyrna
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 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 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 informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationCENTROID (AĞIRLIK MERKEZİ )
CENTOD (ĞLK MEKEZİ ) centrod s a geometrcal concept arsng from parallel forces. Tus, onl parallel forces possess a centrod. Centrod s tougt of as te pont were te wole wegt of a pscal od or sstem of partcles
More informationNumerical Solution of Ordinary Differential Equations
College of Engneerng and Computer Scence Mecancal Engneerng Department otes on Engneerng Analss Larr Caretto ovember 9, 7 Goal of tese notes umercal Soluton of Ordnar Dfferental Equatons ese notes were
More informationComputational Electromagnetics in Antenna Analysis and Design
Computatonal Electromagnetcs n Antenna Analyss and Desgn Introducton It s rare for real-lfe EM problems to fall neatly nto a class that can be solved by the analytcal methods presented n the precedng lectures.
More informationFundamental loop-current method using virtual voltage sources technique for special cases
Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,
More informationShuai Dong. Isaac Newton. Gottfried Leibniz
Computatonal pyscs Sua Dong Isaac Newton Gottred Lebnz Numercal calculus poston dervatve ntegral v velocty dervatve ntegral a acceleraton Numercal calculus Numercal derentaton Numercal ntegraton Roots
More informationCENTROID (AĞIRLIK MERKEZİ )
CENTOD (ĞLK MEKEZİ ) centrod s a geometrcal concept arsng from parallel forces. Tus, onl parallel forces possess a centrod. Centrod s tougt of as te pont were te wole wegt of a pscal od or sstem of partcles
More information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
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 information1 Introduction We consider a class of singularly perturbed two point singular boundary value problems of the form: k x with boundary conditions
Lakshm Sreesha Ch. Non Standard Fnte Dfference Method for Sngularly Perturbed Sngular wo Pont Boundary Value Problem usng Non Polynomal Splne LAKSHMI SIREESHA CH Department of Mathematcs Unversty College
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationConvexity 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 informationThermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)
hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In
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 informationA Discrete Approach to Continuous Second-Order Boundary Value Problems via Monotone Iterative Techniques
Internatonal Journal of Dfference Equatons ISSN 0973-6069, Volume 12, Number 1, pp. 145 160 2017) ttp://campus.mst.edu/jde A Dscrete Approac to Contnuous Second-Order Boundary Value Problems va Monotone
More information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
More informationVisco-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 informationSTUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS
Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
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 informationTHE 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 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 informationThree-dimensional eddy current analysis by the boundary element method using vector potential
Physcs Electrcty & Magnetsm felds Okayama Unversty Year 1990 Three-dmensonal eddy current analyss by the boundary element method usng vector potental H. Tsubo M. Tanaka Okayama Unversty Okayama Unversty
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 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 informationArmy Ants Tunneling for Classical Simulations
Electronc Supplementary Materal (ESI) for Chemcal Scence. Ths journal s The Royal Socety of Chemstry 2014 electronc supplementary nformaton (ESI) for Chemcal Scence Army Ants Tunnelng for Classcal Smulatons
More informationNumerical Nonlinear Analysis with the Boundary Element Method
Blucher Mechancal Engneerng Proceedngs May 2014, vol. 1, num. 1 www.proceedngs.blucher.com.br/evento/10wccm Numercal Nonlnear Analyss wth the Boundary Element Method E. Pneda 1, I. Vllaseñor 1 and J. Zapata
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 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 informationImage classification. Given the bag-of-features representations of images from different classes, how do we learn a model for distinguishing i them?
Image classfcaton Gven te bag-of-features representatons of mages from dfferent classes ow do we learn a model for dstngusng tem? Classfers Learn a decson rule assgnng bag-offeatures representatons of
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
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 informationEffect of loading frequency on the settlement of granular layer
Effect of loadng frequency on the settlement of granular layer Akko KONO Ralway Techncal Research Insttute, Japan Takash Matsushma Tsukuba Unversty, Japan ABSTRACT: Cyclc loadng tests were performed both
More informationApplication of Pontryagin s Maximum Principles and Runge- Kutta Methods in Optimal Control Problems
IOSR Journal of Matematcs (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue 5 Ver. II (Sep. - Oct. 5), PP 43-63 www.osrjournals.org Applcaton of Pontryagn s Maxmum Prncples and Runge- Kutta Metods
More informationNumerical Differentiation
Part 5 Capter 19 Numercal Derentaton PowerPonts organzed by Dr. Mcael R. Gustason II, Duke Unversty Revsed by Pro. Jang, CAU All mages copyrgt Te McGraw-Hll Companes, Inc. Permsson requred or reproducton
More informationPhysics 53. Rotational Motion 3. Sir, I have found you an argument, but I am not obliged to find you an understanding.
Physcs 53 Rotatonal Moton 3 Sr, I have found you an argument, but I am not oblged to fnd you an understandng. Samuel Johnson Angular momentum Wth respect to rotatonal moton of a body, moment of nerta plays
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
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 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 informationAnalysis of thermal processes occurring in heated multilayered metal films using the dual-phase lag model
Arc. Mec., 69, 4 5, pp. 275 287, Warszawa 2017 Analyss of termal processes occurrng n eated multlayered metal flms usng te dual-pase lag model E. MAJCHRZAK, G. KAŁUŻA Insttute of Computatonal Mecancs and
More informationCOMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS
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 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 informationA Computational Method for Solving Two Point Boundary Value Problems of Order Four
Yoge Gupta et al, Int. J. Comp. Tec. Appl., Vol (5), - ISSN:9-09 A Computatonal Metod for Solvng Two Pont Boundary Value Problem of Order Four Yoge Gupta Department of Matematc Unted College of Engg 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 informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Finite Element Method - Jacques-Hervé SAIAC
COMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Fnte Element Method - Jacques-Hervé SAIAC FINITE ELEMENT METHOD Jacques-Hervé SAIAC Départment de Mathématques, Conservatore Natonal des Arts et Méters, Pars,
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationModeling of Dynamic Systems
Modelng of Dynamc Systems Ref: Control System Engneerng Norman Nse : Chapters & 3 Chapter objectves : Revew the Laplace transform Learn how to fnd a mathematcal model, called a transfer functon Learn how
More informationComparison of Newton Raphson and Hard Darcy methods for gravity main nonlinear water network
IOSR Journal of Mecancal and Cvl Engneerng (IOSRJMCE) eissn:,pissn: 0X, Volume, Issue Ver. III (May. June. 0), PP 0 www.osrjournals.org Comparson of Newton Rapson and Hard Darcy metods for gravty man nonlnear
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 information