HIGH ORDER SHEAR DEFORMATION THEORY FOR ISOTROPIC BEAMS WITH USING MIXED FINITE ELEMENT METHOD
|
|
- Caitlin Rich
- 5 years ago
- Views:
Transcription
1 Held on 1-, Nov, 15, in Dubai, ISBN: HIGH ORDER SHEAR DEFORATION THEORY FOR ISOTROPIC BEAS WITH USING IXED FINITE ELEENT ETHOD Emra adenci, Department of Civil Engineering, Necmettin Erbakan University, Köyceğiz Campus 49 eram, Konya, TURKEY Atilla Özütok Department of Civil Engineering, KTO Karatay University, Akabe Str. 4 Karatay, Konya, TURKEY Abstract Tis paper a mied finite element metod for te static analysis of isotropic beams are investigated wit considering ig order sear deformation beam teory. By using virtual displacement principle, field equation of beam wic is including parabolic sear deformation function is derived. In te transformed differential equation of ig order beam element to functional form, te Gâteau differential metod is used tat functional as ten independent variables. Different tree beam model about parabolic sear function presented for tree beam models and HOBT1 beam element is developed. Te element matri is formulated by using mied finite element metod. A computer program using FORTRAN language is developed to eecute te analyses for te present study. Te performance of te beam elements for static analysis is verified wit a good accuracy by te solution of numerical eamples present in te literature. Keywords Hig order sear deformation teory, finite element metod, Gâteau differential metod, Energy principle. INTRODUCTION In te beam bending, te investigations of problems ave been limited to te elementary Classical beam teory and First order beam teory [1]. Te classical beam teory (EBT) is inaccurate for reasonably tick beams []. It is widely used for tin beams and inapplicable to neglect te transverse sear deformation. By Timosenko, due to drawbacks of te classical beam teory developed a beam teory wic is as known FSDT to include te effect of te transverse sear deformation. In te FSDT, assumes a constant sear strain across te tickness of te beam and requires a problem dependent 6 sear correction factor [3]. To overcome te limitations in te EBT and FSDT, many te ig order beam teories ave been developed. Hig order teories account for sear deformation effects by iger order variations of inplane and transverse displacements troug te tickness and satisfy te equilibrium conditions on te top and bottom surfaces of te beam [4]. By using ig order teory, te sear deformation factors or effects ave been represented using parabolic, yperbolic, trigonometric and eponential functions troug te tickness of te beam [5-9]. On te contrary, Reddy [1] epresses tat: Higerorder teories can represent te kinematics better, may not require sear correction factors, and can yield more accurate interlaminar stress distributions. However, tey involve iger-order stress resultants tat are difficult to interpret pysically and require considerably more computational effort. Terefore, suc teories sould be used only wen necessary. Solution to differential equations is not always possible due to compleity of boundary conditions. Te finite element metod (FE) is te most suitable coice for structural analysis, because of its versatility in andling comple geometries and boundary conditions. Te mied type FE metod proved itself for static and dynamic analyses of various structural members and etensively used by researcers because mied-type FE model as greater importance since te forces and moments can be calculated wit less number of elements but more sensitivity [11]. In te FE, starting point is obtaining te function wic is corresponding to differential equation. Recently te Gâteau differential (GD) metod as been employed by Aköz and oter autors [1-17] to construct te functional for various problems. Te GD as been successfully applied to obtaining te functionals for te static and dynamic analysis beams, plates and sells. In tis study, te static responses of isotropic
2 Held on 1-, Nov, 15, in Dubai, ISBN: beam basis on ig order sear deformation beam models are presented using mied type FE formulation. Te virtual displacement principle was applied to obtaining te field equations written in operator form by including te dynamic and geometric boundary conditions. Te GD approac is employed to construct te functional using te variational metod. So using te mied type FE, te internal forces, moments, rotations, displacements of beams can be calculated easily. In order to model te beam wit two nodes, namely HOBT1 was derived wit ten degrees of freedom. A computer program was developed by using FORTRAN for analysis. Formulation In tis study, te displacement field of te eisting HOBT is given by function form; U (, y, z) zw f ( z) V (, y, z), W (,, y, z) w( ) were w and are te unknown displacement along te beam plane and rotation of te normal to te cross section about y-ais at te plane, respectively as illustrate in Fig.1. And f( z) represent parabolic sape functions determining te distribution of te transverse sear strains and stresses along te tickness. z z f( z) ( ) (a) 4 3 5z 4z f( z) (1 ) (b) 4 3 4z f ( z) z(1 ) (c) 3 In Eq.. Te contributions to constructing te functions are respectively ascribed to Ambartsumain [18] (Eq. a); Kaczkowski [19], Panc [] and Reissner [1] (Eq. b); Reddy [], urty [3] and Levinson [4]. (1) Fig. 1 Coordinate system of beam element Te strain-displacement relationsip under te assumption of small deformation and small rotations for te teories is given as U zw f ( z) U W z f ( z ), z z,, Te state of stress is given by te Hooke s Law as follows [1] Q z Were y 11 Q55 z Qij are te well-known reduced stiffnesses. Te principle of virtual displacements is used to derive te governing equations wit te displacement formulations as following [16] L z z dad q wd Ω ( ) ( ) (5) Were q ( ) is te distributed transverse load per unit lengt and is variational symbol; substituted from Eq. 4, and z and and (3) (4) z are z are substituted from Eq. 3, Euler-Lagrange equations are obtaining as follow L b 63
3 Held on 1-, Nov, 15, in Dubai, ISBN: w: q, : Q Here, is bending moment, moment and z (6) is refined bending Qz is refined sear force. Force and moment resultants of present teory are defined in te following form [7] zdz f ( z) dz Q f ( z) dz z z, z Using te parabolic sear deformation teory, by substituting te stress-strain relations into te definitions of force and moment resultants given in Eq. (7) te following te constitutive equations of te beam are obtained D11 F11 w, F H 11 11, Q z A 55 Te stiffness coefficients D 11, F11and H11 sear stiffness A55 are given by ( ) 11( ( ) ( ) ) D F H Q z f z z f z dz A Q f ( z) dz 55 55, z (7) (8) and transverse Dynamic and geometric boundary conditions are written in symbolic form as Eq. (1). Eplicit forms of te boundary condition will be obtained after some variational manipulations. Tere representing te moment, force, deflection and rotation vectors [1]. Dynamic boundary conditions are (9) R R Geometric boundary conditions are uu ΩΩ (1a) (1b) Te Eq. (6) and Eq. (8) can be written togeter wic are field equations q, Q, D w F 11, 11, F w H 11, 11, Q A z 55 z (11) Tese equations can be written in operator form similar to elastic beam [11] Q Py f (1) Te eplicit form of te operator is Pij yi fi 1 u R 1 Ω 1 Ω 1 R u If Q is a potential operator, te equality (13) * * dq(y,y),y dq(y,y ),y (14) ust be satisfied [5]. Tere are te Gâteau derivatives of te operator wic are defined Q( y τ y) dq( y, y ) (15) τ τ Were is τ is a scaler. Since te operator is potential ten te functional corresponding to te field equations is obtained as [5] 64
4 Held on 1-, Nov, 15, in Dubai, ISBN: ( y) 1 I Q(sy),y ds (16) Were is s is scaler quantity. Te functional of ig order sear deformation beam is in similar closed form I( y),, w,,,,,, Qz, Qz,, boundary conditions Q q w z (17) Initially te interpolation function must be cosen to derive te mied finite element formulation [8]. If a one dimensional element wit a parent sape function 1 1 i 1, j 1, i 1, (18) is used, te element matri can be obtained eplicitly as 1 1 k d k d 1 i j i, j, 1 1 k 1 3 i j, 1 d Analysis of Results (19) In tis study, te numerical results computed using te finite elements developed based on te above ig order sear deformation teories are compared wit various parabolic sear deformation functions wic are model1 is about Eq. (a), model is about Eq. (b) and model3 is about Eq. (c). Teoretical results calculated for various lengts to tickness ratio. Te boundary conditions of te beams were denoted by S for simply supported edge. Wit te application of te mied-type FE metod, HOBT1 beam elements were derived for HOBT were te HOBT1 elements ave a total of ten unknowns of displacements, bending moments, refined bending moments, refined sear forces and rotations. Troug te numerical eamples, te static analyses of beams were investigated using HOBT1 finite elements. TABLE I. THE AXIU DISPLACEENT OF AN ISOTROPIC SS THICK BEAS Teory L/ 1/1 4/1 8/1 16/1 HOBT1 model HOBT1 model HOBT1 model [1] [6] [7] In te Table 1., te maimum displacements of isotropic simply supported tick beams are calculated, under uniform distributed load q 1 and various L/ are considered for different beam teories and solutions of oter studies in literature. Te material and geometric properties are E=9, v=.3, b=1. Te values for te ig order elements are taken from te paper by [1], [6] and [7]. And also given in Table., same analysis presented for tin beams and te results are compared oter references [1], [6]. Eac model is a good comparison for various cross section states. Te ig order beam element by using mied finite element metod gives ecellent results wen te Gâteau differential metod is used for tick and tin beams. TABLE II. THE AXIU DISPLACEENT OF AN ISOTROPIC SS THIN BEAS Teory L/ 1/1 4/1 8/1 16/1 HOBT1 model HOBT1 model HOBT1 model [1] [6]
5 Held on 1-, Nov, 15, in Dubai, ISBN: Conclusions In tis study, te GD approac is employed to construct te functional using te variational metod for static analysis of isotropic beams for realized to solving mied finite element formulation. In order to develop a mied finite element model for applied to beams, it is necessary to establis a functional. Using te mied type FE, te internal forces, moments, rotations, displacements of isotropic beams can be calculated easily. Te Gâteau differential approaces ave some important advantages [15]. Results of te isotropic beam elements are verified by te oter epression given by [1], [6] and [7]. Te comparison study sows tat results of free vibration of ig order isotropic beams better tan te oter teories and te finite elements models. REFERENCES [1] P.R. Heyliger, and J.N.Reddy, A Higer Order Beam Finite Element for Bending and Vibration Problems, J.Sound. Vib., vol. 16, No., pp , [] L. Jun, and H. Honging Dynamic Stiffness Analysis of Laminated Composite Beams Using Trigonometric Sear Deformation Teory, Composite Strucutres, vol. 89, pp , 9. [3] P. Subramanian, Dynamic Analysis of Laminated Composite Beams Using Higer Order Teories and Finite Elements, Composite Structures, vol.73, pp , 6. [4] H.T. Tai, D.H. Coi, A Simple First-order Sear Deformation Teory for Laminated Composite Plates, Composite Structures, vol. 16, pp , 13. [5] F.G. Yuan, and R.E. iller, A Higer Order Finite Element for Laminated Beams, Composite Structures, vol. 14, no., pp , 199. [6] H. Yu, A Higer-order Finite Element for Analysis of Composite Laminated Structures, Composite Structures, vol. 8, no.4, pp , [7] G. Si, K. Lam, and T. Tay, On Efficient Finite Element odeling of Composite Beams and Plates Using Hig-order Teories and Accurate Composite Beam Element, Composite Structures, vol. 41, no., pp , [8] P. Subramanian, Fleural Analysis of Symmetric Laminated Composite Beams Using C1 Finite Element, Composite Structures, vol. 54, no.1, pp , [9].V.V.S. urty, D.R. aapatra, K. Badarinarayana and S.A. Gopalakrisnan, A Refined Higer Order Finite Element for Asymmetric Composite Beams Composite Structures, vol. 67, no.1, pp. 7 35, 5. [1] J.N. Reddy, ecanics of Laminated Composite Plates and Sells: Teory and Analysis, nd ed., Boca-Rotan: CRC Press, 4. [11] A.Y. Aköz and F. Kadioğlu, Te ied Finite Element etod for te Quasi-Static and Dynamic Analysis of Viscoelastic Timosenko Beam, Int. J. for Num. etods in Eng., vol. 44, no. 1, pp , [1] A.Y. Aköz,.H. Omurtag and A.N. Doğruoğlu, Te ied Finite Element Formulation for Tree Dimensional Bars, Int. J. Solids Structure, vol.8, pp. 5-34, [13].H. Omurtag and A.Y. Aköz, A Compatible Cylindrical Sell Element for Stiffened Cylindrical Sells in A ied Finite Element Formulation, Computers Structures, vol.49, no., pp , [14] A. Y. Akoz and F. Kadıoğlu, Te ied Finite Element Solution of Circular Beam on Elastic Foundation, Computers Structures, vol.6, no.4, pp , [15] A. Y. Akoz and A. Özutok, A Functional for Sells of Arbitrary Geometry and A ied Finite Element etod for Parabolic and Circular Cylindrical Sells, Int. J. for Num. etods in Eng., vol.47, no. 1, pp ,. [16] A. Ozutok and E. adenci, Free Vibration Analysis of Cross-ply Laminated Composite Beams by ied Finite Element Formulation, Int. J. of Struc. Stabil. and Dyn., Vol. 13, No., pp , 13. [17] A. Ozutok, E. adenci, and F. Kadıoğlu, Free Vibration Analysis of Angle-ply Laminate Composite Beams by ied Finite Element Formulation Using te Gâteau Differential, Science and Composite aterials, vol. 1, no., pp , 14. [18] S.A. Ambartsumian, On te Teory of Bending plates, Izv Otd Tec Nauk, ANN SSSR, vol. 5, pp , [19] Z. Kaczkowski, Plates-statistical calculations, Warsaw: Arkady, [] V. Panc, Teories of elastic plates, Prague, Academia; [1] E. Reissner, On Transverse Bending of Plates Including te Effects of Transverse Sear
6 Held on 1-, Nov, 15, in Dubai, ISBN: Deformation, Int. J. Solids Struc., vol. 5, pp , [] J.N. Reddy, A Simple Higer-order Teory for Laminated Composite Plates, J. Appl. ec., vol. 51, pp , [3].V.V. urty, An Improved Transverse Sear Deformation Teory for Laminated Anisotropic Plates, NASA Tecnical Paper, [4]. Levinson, An Accurate Simple Teory of Statics and Dynamics of Elastic Plates, ec.rese. Comm., vol. 7, pp , 198. [5] J.T. Oden and J.N. Reddy, Variational etods in Teoretical ecanics, Berlin: Sipringer-Verlag, [6] W.B. Bickford, A Consistent Higer Order Beam Teory, Development in Teo. App. ec., Vol. 11, pp , 198. [7] T.P. Vo, and H.T. Tai, Static Beavior of Composite Beams Using Various Refined Sear Deformation Teories, Composite Structures, Vol. 94, pp , 1. [8] J.N. Reddy, An Introduction to te Finite Element etod, New York: cgraw-hill,
Stability of Smart Beams with Varying Properties Based on the First Order Shear Deformation Theory Located on a Continuous Elastic Foundation
Australian Journal of Basic and Applied Sciences, 5(7): 743-747, ISSN 99-878 Stability of Smart Beams wit Varying Properties Based on te First Order Sear Deformation Teory ocated on a Continuous Elastic
More informationWINKLER PLATES BY THE BOUNDARY KNOT METHOD
WINKLER PLATES BY THE BOUNARY KNOT ETHO Sofía Roble, sroble@fing.edu.uy Berardi Sensale, sensale@fing.edu.uy Facultad de Ingeniería, Julio Herrera y Reissig 565, ontevideo Abstract. Tis paper describes
More informationFULL SHEAR DEFORMATION FOR ANALYSIS OF THICK PLATE
FULL SHAR DFORMATION FOR ANALYSIS OF THICK PLAT 1 IBARUGBULM, OWUS M., 2 NJOKU, KLCHI O., 3 ZIFULA, UCHCHI G. 1,2,3 Civil ngineering Department, Federal University of Tecnology, Owerri, Nigeria -mail:
More information6. Non-uniform bending
. Non-uniform bending Introduction Definition A non-uniform bending is te case were te cross-section is not only bent but also seared. It is known from te statics tat in suc a case, te bending moment in
More informationEuler-Bernoulli Beam Theory in the Presence of Fiber Bending Stiffness
IOSR Journal of Matematics (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 13, Issue 3 Ver. V (May - June 017), PP 10-17 www.iosrjournals.org Euler-Bernoulli Beam Teory in te Presence of Fiber Bending
More informationStatic Response Analysis of a FGM Timoshenko s Beam Subjected to Uniformly Distributed Loading Condition
Static Response Analysis of a FGM Timoseno s Beam Subjected to Uniformly Distributed Loading Condition 8 Aas Roy Department of Mecanical Engineering National Institute of Tecnology Durgapur Durgapur, Maatma
More informationFlapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using the differential transform method
Meccanica 2006) 41:661 670 DOI 10.1007/s11012-006-9012-z Flapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using te differential transform metod Ozge Ozdemir Ozgumus
More informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
More informationFlexural Analysis of Deep Aluminum Beam
Journal of Soft Computing in Civil Engineering -1 (018) 71-84 journal homepage: http://www.jsoftcivil.com/ Fleural Analysis of Deep Aluminum Beam P. Kapdis 1, U. Kalwane 1, U. Salunkhe 1 and A. Dahake
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationMODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM
MODIFIED HYPERBOLIC SHEAR DEFORMATION THEORY FOR STATIC FLEXURE ANALYSIS OF THICK ISOTROPIC BEAM S. Jasotharan * and I.R.A. Weerasekera University of Moratuwa, Moratuwa, Sri Lanka * E-Mail: jasos91@hotmail.com,
More informationNonlinear correction to the bending stiffness of a damaged composite beam
Van Paepegem, W., Decaene, R. and Degrieck, J. (5). Nonlinear correction to te bending stiffness of a damaged composite beam. Nonlinear correction to te bending stiffness of a damaged composite beam W.
More informationAN ANALYSIS OF AMPLITUDE AND PERIOD OF ALTERNATING ICE LOADS ON CONICAL STRUCTURES
Ice in te Environment: Proceedings of te 1t IAHR International Symposium on Ice Dunedin, New Zealand, nd t December International Association of Hydraulic Engineering and Researc AN ANALYSIS OF AMPLITUDE
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationBending analysis of a functionally graded piezoelectric cantilever beam
Science in Cina Series G: Pysics Mecanics & Astronomy 7 Science in Cina Press Springer-Verlag Bending analysis of a functionally graded pieoelectric cantilever beam YU Tao & ZHONG Zeng Scool of Aerospace
More informationAnalysis of Stress and Deflection about Steel-Concrete Composite Girders Considering Slippage and Shrink & Creep Under Bending
Send Orders for Reprints to reprints@bentamscience.ae Te Open Civil Engineering Journal 9 7-7 7 Open Access Analysis of Stress and Deflection about Steel-Concrete Composite Girders Considering Slippage
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationStress Analysis of Laminated Composite and Sandwich Beams using a Novel Shear and Normal Deformation Theory
134 Stress Analysis of Laminated Composite and Sandwich Beams using a Novel Shear and Normal Deformation Theory Abstract A novel Normal and Shear Deformation Theory (NSDT) for analysis of laminated composite
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationJournal of Mechanical Science and Technology 23 (2009) 2072~2084. M. M. Najafizadeh * and M. R. Isvandzibaei
Journal of Mecanical Science and Tecnology (9 7~84 Journal of Mecanical Science and Tecnology www.springerlink.com/content/78-494x DOI.7/s6-9-4- Vibration of functionally graded cylindrical sells based
More information3. Using your answers to the two previous questions, evaluate the Mratio
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0219 2.002 MECHANICS AND MATERIALS II HOMEWORK NO. 4 Distributed: Friday, April 2, 2004 Due: Friday,
More informationNonlinear RC beam element model under combined action of axial, bending and shear
icccbe 21 Nottingham University Press Proceedings of the International Conference on Computing in Civil and Building ngineering W Tiani (ditor) Nonlinear RC beam element model under combined action of
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationModeling of the Bending Stiffness of a Bimaterial Beam by the Approximation of One-Dimensional of Laminated Theory
. Flores-Domínguez Int. Journal of Engineering Research and Applications RESEARCH ARTICLE OPEN ACCESS odeling of the Bending Stiffness of a Bimaterial Beam by the Approimation of One-Dimensional of Laminated
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationNCCI: Simple methods for second order effects in portal frames
NCC: Simple metods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames NCC: Simple metods for second order effects in portal frames Tis NCC presents information
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationAnalysis of Thick Cantilever Beam Using New Hyperbolic Shear Deformation Theory
International Journal of Research in Advent Technology, Vol.4, No.5, May 16 E-ISSN: 1-967 Analysis of Thick Cantilever Beam Using New Hyperbolic Shear Deformation Theory Mr. Mithun. K. Sawant 1, Dr. Ajay.
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationN igerian Journal of M athematics and Applications V olume 23, (2014), 1 13
N igerian Journal of M atematics and Applications V olume 23, (24), 3 c N ig. J. M at. Appl. ttp : //www.kwsman.com CONSTRUCTION OF POLYNOMIAL BASIS AND ITS APPLICATION TO ORDINARY DIFFERENTIAL EQUATIONS
More informationFinite Element Analysis of J-Integral for Surface Cracks in Round Bars under Combined Mode I Loading
nternational Journal of ntegrated Engineering, Vol. 9 No. 2 (207) p. -8 Finite Element Analysis of J-ntegral for Surface Cracks in Round Bars under Combined Mode Loading A.E smail, A.K Ariffin 2, S. Abdulla
More informationChapter XI. Solution of Ordinary Differential Equations
dy d 7 Capter XI Solution of Ordinary Differential Equations Several approaces are described for te solution of ODE (ordinary differential equations). Tese include: Te Taylor series Metod Te Euler Metod
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES THIRD-ORDER SHEAR DEFORMATION PLATE THEORY LAYERWISE LAMINATE THEORY J.N. Reddy 1 Third-Order Shear Deformation Plate Theory Assumed Displacement Field µ u(x y z t) u 0 (x y t) +
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationP. M. Pankade 1, D. H. Tupe 2, G. R. Gandhe 3
ISSN: 78 7798 Volume 5, Issue 5, May 6 Static Fleural Analysis of Thick Beam Using Hyperbolic Shear Deformation Theory P. M. Pankade, D. H. Tupe, G. R. Gandhe P.G. Student, Dept. of Civil Engineering,
More informationFabric Evolution and Its Effect on Strain Localization in Sand
Fabric Evolution and Its Effect on Strain Localization in Sand Ziwei Gao and Jidong Zao Abstract Fabric anisotropy affects importantly te overall beaviour of sand including its strengt and deformation
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationSample Problems for Exam II
Sample Problems for Exam 1. Te saft below as lengt L, Torsional stiffness GJ and torque T is applied at point C, wic is at a distance of 0.6L from te left (point ). Use Castigliano teorem to Calculate
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationThe tri-harmonic plate bending equation
Computational Metods and Experimental Measurements XVII 67 Te tri-armonic plate bending equation I. D. Kotcergenko Kotcergenko Engenaria Ltda., Brazil Abstract Tis is te tird paper publised at WIT Transactions
More informationFlexure of Thick Cantilever Beam using Third Order Shear Deformation Theory
International Journal of Engineering Research and Development e-issn: 78-67X, p-issn: 78-8X, www.ijerd.com Volume 6, Issue 1 (April 13), PP. 9-14 Fleure of Thick Cantilever Beam using Third Order Shear
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationRobotic manipulation project
Robotic manipulation project Bin Nguyen December 5, 2006 Abstract Tis is te draft report for Robotic Manipulation s class project. Te cosen project aims to understand and implement Kevin Egan s non-convex
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationStress analysis of laminated glass with different interlayer materials
Alexandria Engineering Journal (01) 51, 61 67 Alexandria University Alexandria Engineering Journal www.elsevier.com/locate/aej www.sciencedirect.com Stress analysis of laminated glass wit different interlayer
More information3.1 Extreme Values of a Function
.1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationNon-linear Analysis Method of Ground Response Using Equivalent Single-degree-of-freedom Model
Proceedings of te Tent Pacific Conference on Eartquake Engineering Building an Eartquake-Resilient Pacific 6-8 November 25, Sydney, Australia Non-linear Analysis Metod of Ground Response Using Equivalent
More informationSome Review Problems for First Midterm Mathematics 1300, Calculus 1
Some Review Problems for First Midterm Matematics 00, Calculus. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd,
More informationACCURATE SYNTHESIS FORMULAS OBTAINED BY USING A DIFFERENTIAL EVOLUTION ALGORITHM FOR CONDUCTOR-BACKED COPLANAR WAVEG- UIDES
Progress In Electromagnetics Researc M, Vol. 10, 71 81, 2009 ACCURATE SYNTHESIS FORMULAS OBTAINED BY USING A DIFFERENTIAL EVOLUTION ALGORITHM FOR CONDUCTOR-BACKED COPLANAR WAVEG- UIDES S. Kaya, K. Guney,
More informationReinforced Concrete Structures/Design Aids
AREA OF BARS (mm 2 ) Reinforced Concrete Structures/Design Aids Size of s (mm) Number of s 1 2 3 4 5 6 7 8 8 50 101 151 201 251 302 352 402 10 79 157 236 314 393 471 550 628 12 113 226 339 452 566 679
More informationz = - 2rriw ctg k!! = - 4rrioo [1 - _!_ (koh)2... ] '
SOVIET PHYSICS JETP VOLUME 5, NUMBER 3 OCTOBER, 957 Skin Effect in Tin Films and Wires B. M.BOLOTOVSKII P. N. Lebedev Pysical Institute, Academy of Sciences, USSR (Submitted to JETP editor February 9,
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationJOURNAL OF APPLIED SCIENCES RESEARCH
Copyrigt 15 American-Eurasian Network for Scientific Information publiser JOURNA OF APPIED SCIENCES RESEARCH ISSN: 1819-5X EISSN: 181-157X JOURNA ome page: ttp://www.aensiweb.com/jasr 15 Special; 11(17):
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationBending of Shear Deformable Plates Resting on Winkler Foundations According to Trigonometric Plate Theory
J. Appl. Comput. ech., 4(3) (018) 187-01 DOI: 10.055/JAC.017.3057.1148 ISSN: 383-4536 jacm.scu.ac.ir Bending of Shear Deformable Plates Resting on Winkler Foundations According to Trigonometric Plate Theor
More informationThe entransy dissipation minimization principle under given heat duty and heat transfer area conditions
Article Engineering Termopysics July 2011 Vol.56 No.19: 2071 2076 doi: 10.1007/s11434-010-4189-x SPECIAL TOPICS: Te entransy dissipation minimization principle under given eat duty and eat transfer area
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More information4.5 The framework element stiffness matrix
45 The framework element stiffness matri Consider a 1 degree-of-freedom element that is straight prismatic and symmetric about both principal cross-sectional aes For such a section the shear center coincides
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationComputational Method of Structural Reliability Based on Integration Algorithms
Sensors & ransducers, Vol. 54, Issue 7, July 03, pp. 5-59 Sensors & ransducers 03 by IFSA ttp://www.sensorsportal.com Computational Metod of Structural Based on Integration Algoritms * Cong Cen, Yi Wan
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationIntroduction to Finite Element Method. Dr. Aamer Haque
Introduction to Finite Element Method 4 th Order Beam Equation Dr. Aamer Haque http://math.iit.edu/~ahaque6 ahaque7@iit.edu Illinois Institute of Technology July 1, 009 Outline Euler-Bernoulli Beams Assumptions
More information1. Which one of the following expressions is not equal to all the others? 1 C. 1 D. 25x. 2. Simplify this expression as much as possible.
004 Algebra Pretest answers and scoring Part A. Multiple coice questions. Directions: Circle te letter ( A, B, C, D, or E ) net to te correct answer. points eac, no partial credit. Wic one of te following
More informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
More informationInvestigating Euler s Method and Differential Equations to Approximate π. Lindsay Crowl August 2, 2001
Investigating Euler s Metod and Differential Equations to Approximate π Lindsa Crowl August 2, 2001 Tis researc paper focuses on finding a more efficient and accurate wa to approximate π. Suppose tat x
More informationLecture 15 Strain and stress in beams
Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationChapter 2 Limits and Continuity
4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationHIGHER-ORDER THEORIES
HIGHER-ORDER THEORIES Third-order Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain
More informationDerivative at a point
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Derivative at a point Wat you need to know already: Te concept of liit and basic etods for coputing liits. Wat you can
More informationA NEW REFINED THEORY OF PLATES WITH TRANSVERSE SHEAR DEFORMATION FOR MODERATELY THICK AND THICK PLATES
A NEW REFINED THEORY OF PLATES WITH TRANSVERSE SHEAR DEFORMATION FOR MODERATELY THICK AND THICK PLATES J.M. MARTÍNEZ VALLE Mechanics Department, EPS; Leonardo da Vinci Building, Rabanales Campus, Cordoba
More informationARTICLE IN PRESS. Physica E
Physica E 41 (2009) 1651 1655 Contents lists available at ScienceDirect Physica E journal homepage: www.elsevier.com/locate/physe A general nonlocal beam theory: Its application to nanobeam bending, buckling
More informationA general articulation angle stability model for non-slewing articulated mobile cranes on slopes *
tecnical note 3 general articulation angle stability model for non-slewing articulated mobile cranes on slopes * J Wu, L uzzomi and M Hodkiewicz Scool of Mecanical and Cemical Engineering, University of
More informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
More informationA Numerical Scheme for Particle-Laden Thin Film Flow in Two Dimensions
A Numerical Sceme for Particle-Laden Tin Film Flow in Two Dimensions Mattew R. Mata a,, Andrea L. Bertozzi a a Department of Matematics, University of California Los Angeles, 520 Portola Plaza, Los Angeles,
More informationLarge eddy simulation of turbulent flow downstream of a backward-facing step
Available online at www.sciencedirect.com Procedia Engineering 31 (01) 16 International Conference on Advances in Computational Modeling and Simulation Large eddy simulation of turbulent flow downstream
More informationDifferentiation. Area of study Unit 2 Calculus
Differentiation 8VCE VCEco Area of stud Unit Calculus coverage In tis ca 8A 8B 8C 8D 8E 8F capter Introduction to limits Limits of discontinuous, rational and brid functions Differentiation using first
More informationThe causes of the observed head-in-pillow (HnP) soldering defects in BGA
ORIGINAL ARTICLE Solder joints in surface mounted IC assemblies: Relief in stress and warpage owing to te application of elevated stand-off eigts E. Suir E. Suir. Solder joints in surface mounted IC assemblies:
More informationResearch Article Elastic Analysis of Rotating Thick Truncated Conical Shells Subjected to Uniform Pressure Using Disk Form Multilayers
ISRN Mecanical Engineering Article ID 764837 1 pages ttp://d.doi.org/1.1155/214/764837 Researc Article Elastic Analysis of Rotating Tick Truncated Conical Sells Subjected to Uniform Pressure Using Disk
More informationFlexural analysis of deep beam subjected to parabolic load using refined shear deformation theory
Applied and Computational Mechanics 6 (2012) 163 172 Flexural analysis of deep beam subjected to parabolic load using refined shear deformation theory Y. M. Ghugal a,,a.g.dahake b a Applied Mechanics Department,
More informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationOn the computation of wavenumber integrals in phase-shift migration of common-offset sections
Computation of offset-wavenumber integrals On te computation of wavenumber integrals in pase-sift migration of common-offset sections Xiniang Li and Gary F. Margrave ABSTRACT Te evaluation of wavenumber
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationBending Analysis of Symmetrically Laminated Plates
Leonardo Journal of Sciences ISSN 1583-0233 Issue 16, January-June 2010 p. 105-116 Bending Analysis of Symmetrically Laminated Plates Bouazza MOKHTAR 1, Hammadi FODIL 2 and Khadir MOSTAPHA 2 1 Department
More informationLogarithmic functions
Roberto s Notes on Differential Calculus Capter 5: Derivatives of transcendental functions Section Derivatives of Logaritmic functions Wat ou need to know alread: Definition of derivative and all basic
More informationBuckling of Laminated and Functionally Graded Plates using Radial Basis Functions
Paper 223 Buckling of Laminated and Functionally Graded Plates using Radial Basis Functions A.M.A. Neves, A.J.M. Ferreira, E. Carrera 3, M. Cinefra 3 C.M.C. Roque 2, C.M. Mota Soares 4 and R.M.N. Jorge
More informationHomework 1. Problem 1 Browse the 331 website to answer: When you should use data symbols on a graph. (Hint check out lab reports...
Homework 1 Problem 1 Browse te 331 website to answer: Wen you sould use data symbols on a grap. (Hint ceck out lab reports...) Solution 1 Use data symbols to sow data points unless tere is so muc data
More information