3rd Int. Conference on Inverse Problems in Engineering. Daniel Lesnic. University of Leeds. Leeds, West Yorkshire LS2 9JT
|
|
- Polly Norris
- 6 years ago
- Views:
Transcription
1 Inverse Problems in Engineering: Theory and Practice 3rd Int. Conference on Inverse Problems in Engineering June 13-18, 1999, Port Ludlow, WA, USA APP RETRIEVING THE FLEXURAL RIGIDITY OF A BEAM FROM DEFLECTION MEASUREMENTS Daniel Lesnic Department of Applied Mathematics University of Leeds Leeds, West Yorkshire LS 9JT amt5ld@amsta.leeds.ac.uk ABSTRACT A rigorous investigation into the identication of the heterogeneous exural rigidity coecient from deection measurements recorded along a beam in the presence of a prescribed load is presented. Mathematically, the problem reduces to the need to solve the Euler-Bernoulli steady-state beam equation subject to appropriate boundary conditions. Conditions for the uniqueness and continuous dependence on the input data of the solution of the inverse problem for simply supported beams are established and, in particular, it is shown that the operator which maps an input deection into an output exural rigidity is Holder continuous. Since the inverse problem can be recast in the form of a Fredholm integral equation of the rst kind, numerical results obtained using various methods, such as Tikhonov's regularization, singular value decomposition and mollication are discussed. NOMENCLATURE E Modulus of elasticity I Moment of inertia L Length of the beam M Bending moment a Flexural rigidity f Transversely distributed load p Amount of noise u Transverse deection Inverse of the exural rigidity Regularization parameter Standard deviation INTRODUCTION In the Euler-Bernoulli beam theory, it is assumed that the plane cross-sections perpendicular to the axis of the beam remain plane and perpendicular to the axis after deformation, resulting in the transverse deection u of the beam being governed by the fourth-order ordinary dierential equation d a(x) d u dx dx (x) = f(x) <x<l (1) where L is the length of the beam, f is the transversely distributed load and the spacewise dependent conductivity a = EI (often called exural rigidity) is the product of the modulus of elasticity E and the moment of inertia I of the cross section of the beam about an axis through its centroid at right angles to the cross section. It should be noted that eqn.(1) can be split into a system of two secondorder ordinary dierential equations, namely, M(x) =a(x) d u (x) <x<l () dx d M = f(x) <x<l (3) dx where M is the bending moment. Using a simple nite element methodology, it can be seen that the primary variables 1 Copyright c 1999 by ASME
2 associated with eqn.(1) are the deection u and the slope du=dx, whilst the bending moment M = a(d u=dx ) and the shear force (d=dx)(ad u=dx ) are the secondary variables, any four values of which are usually specied on the boundary, namely, x = orx = L. The direct problem of the Euler-Bernoulli beam theory requires the determination of the deection u which satises eqn.(1) (or M and u satisfying eqns () and (3)), when a> and f are given and four boundary conditions (essential or natural) are prescribed. If at least one of these boundary conditions is on the deection then this direct problem is well-posed. However, there are other problems (termed inverse) that may be associated with the Euler-Bernoulli eqns (1)-(3). For simply supported beams an inverse load identication problem, which requires the determination of the load f(x) which satises eqn.(3) when a and M are given, has been investigated by Collins et al. (1994) in the practical context of recovering engineering loads from strain gauge data. This problem is ill-posed since it violates the stability of the solution, which relates to twice dierentiating numerically M, which is a noisy function. In this study, weinvestigate an inverse coecient identication problem which requires the identication of the positive, heterogeneous exural rigidity coecient a(x) ofa beam which satises eqn.(1), (or eqns () and (3)), when u and f are given. This problem is an extension to higherorder dierential equations of the inverse problem analysed by Marcellini (198) for the one-dimensional Poisson equation. Prior to this study, anumerical algorithm has been recently proposed by Ismailov and Muravey (1996) for supported plates, when f > and u. However, the theoretical investigation of the uniqueness of the solution is much more dicult since eqn.(1) does not, in general, have a unique solution, e.g. if a is a solution of eqn.(1) then for any linear function h, the function a 1 = a + h(d u=dx ) ;1 may still be a solution of eqn.(1). Therefore, a necesary condition for the uniqueness of the solution of eqn.(1) is that the set S = x [ L] j d u dx (x) = is not empty. However, at the other extreme, if S is a subset of [ L] of non-zero measure then a(x) is not identiable on S. Further, based on the physical argument that the properties of the beam, namely E and I, should be positive, and considering only applications in which there is always a positive nite load acting on the beam, we can assume that there are <f 1 f < 1 such thatf 1 f(x) f and (4) that there are < 1 < 1 and Q>such that a(x) belongs to the domain denition set A = fa L 1 ( L) j 1 a ka kqg (5) where k:k denotes the L ( L)-norm. For the present analysis we restrict ourselves to homogeneous boundary conditions to accompany the non-zero load f and the linear eqns (1)-(3). Also, since a> then the set S given by eqn.(4) is S = fx [ L] j M(x) =g. Thus boundary conditions which ensures that S 6= include u () = u (L) or M()M(L) =. If only Dirichlet and/or Neumann type boundary conditions P for u and M are considered then there 4 are in total k= Ck 4 C 4;k 4 = 7 of these possibilities. However, half of these possibilities are equivalent, based on the invariance of eqns (1)-(3) under the translation x 7! (L;x). The aim of this study is not to investigate all these possibilities but rather select realistic physical boundary conditions associated with beams that naturally occur in elasticity, such as beams supported at both ends, namely, u() = u(l) =M() = M(L) = (6) Other types of boundary conditions have been investigated elsewhere, see Lesnic et al. (1999). The plan of the paper is as follows. In section we prove that the operator u 7! a is Holder continuous and thus injective, i.e. the uniqueness and continuous dependence on the input data of the inverse problem. Further, when random noisy discrete input data is included, a regularization algorithm is developed in section 3 and the numerically obtained results are illustrated and discussed. MATHEMATICAL ANALYSIS We formally dene the operator U : A! L ( L) U(a) :=u a = ZZ M a M = ZZ f 8a A (7) as a formal general solution of eqns (1)-(3), where the integral sign R is understood in the sense that the constants of integration are to be determined by imposing the boundary conditions (6). For convenience, and in order to simplify the algebraic manipulations, a uniform load f 1 was assumed which in turn gives Copyright c 1999 by ASME
3 x(x ; L) M(x) = (8) Lemma 3. If a, b A then there is a positive constant C such that ka ; bk L Cku a ; u b k =3 L (14) u a (x) = Z x Z t dt ( ; L) d ; x a() L Lemma 1. The following inequality holds Z t dt ( ; L) d a() (9) ku a ; u bk 3= ku a ; u b k 1= ku a ; u b k (1) Proof: To evaluate the L ;norm ku a ; u b k we useintegration by parts, impose the boundary conditions (6) and apply the Holder inequality toobtain ku a;u bk = ; ku a ; u b k = ; (u a ; u b)(u a ; u b )dx ku a ; u bkku a ; u b k (11) (u a ;u b )(u a;u b )dx ku a ;u b kku a;u b k (1) Finally, combining the inequalities (11) and (1) yields the inequality (1). Lemma. If a A then there is a positive constant C such thatku a kc. Proof: If we denote = a ;1, = ;1 1 then wehave < < 1 and u a = M. Thus u a = M + M and using eqn.(8) we obtain the following estimates for ku a k, namely, ku a kk k L L = (13) 1 Proof: If we denote = a ;1, = b ;1, 1 = ;1 and = ;1,thenwehave < 1 1 < 1 and u a = M, u b = M. The key idea of the proof is based on removing (similar toacauchy nite part evaluation of singular integrals) the singularities of the function = u a=m (and similarly = u b =M) which, based on eqn.(8), is of the type x;1 (L ; x) ;1. Therefore we extend the functions u a, u b, and to the whole real axis by dening them to be zero outside of the interval [ L]. Then, for any >, using the Holder inequality, wehave the following identities and inequalities Z k ; k = j ; j dx (; )[(L; L+) Z + j ; j j u a ; u b j dx (;1 ;)[( L;)[(L+ 1) j M j Z (;1 ;)[( L;)[(L+ 1) ; 4 + ku a ; u b k! 1= 4 x (L ; x) dx =4 + 3= L ku a ; u L ; + j L ; j L ln j L + j b k 1= 4 ( + ku a ; u b k ;1= =L) (15) The right-hand side of expression (15), as a function of, attains its minimum at min = ku a ; u L b k =3 (16) and introducing this value into expression (15) we obtain Finally, since j j= ja j a j a j and thus k k ka k Q, from eqn.(13) it follows that there is a positive constant C>such that ku a kc, which concludes lemma. k ; k ( 4=3 + 1=3 )L ;=3 4=3 ku a ; u b k =3 (17) Finally, the estimate (14) is obtained using the inequality j a ; b j ; 1 j ; j and eqn.(17). 3 Copyright c 1999 by ASME
4 At this stage, using lemmas 1-3, we can conclude the following theorem: Theorem 1. (The Holder continuity of the inverse operator U ;1 (u) =a) If a, b A then there is a positive constant C such that ka ; bk Cku a ; u b k 1=9 (18) As an immediate consequence of this theorem we conclude the uniqueness and continuous dependence of the solution a A on the input data u. REGULARIZATION ALGORITHM In this section the solution of the inverse problem (1), which consists of nding a A when the additional data u a z (19) is given, is based on the rst-order regularization method of Tikhonov and Arsenin (1977), which requires the nding of the minimum with respect to = a ;1 A of the functional () =ku a ; zk + k k () where > is a regularization parameter to be prescribed. For simply supported beams subjected to a unit load we can obtain the deection u(x) as the solution of the second-order dierential equation However, the transformation of the dierential eqn.(1) into the integral eqn.() does not eliminate the ill-posedness of the problem, as such an equation is a Fredholm integral equation of the rst kind and constitutes a well-known example of an ill-posed problem due to its instability of the solution with respect to noise in the input data z as given by eqn.(19), see e.g. Tikhonov and Arsenin (1977). We now divide the interval [ L]into N equal parts by putting s j = jl=n for j = N, and on each subinterval [s j;1 s j ]forj = 1 N,we assume that is constant and takes its value at the midpoint ~s j =(s j + s j;1 )=, i.e. j [sj;1 s j]= (~s j )= j j = 1 N (4) Then eqn.() can be approximated by u(x) = NX j=1 j Z sj s j;1 G(x s)ds x L (5) If the deection u(x) is recorded only at K discrete locations x i =(i ; 1)L=(K) for i = 1 K, then eqn.(5) gives u i = u(x i )= where the integrals A ij = Z sj NX j=1 A ij j i = 1 K (6) s j;1 G(x i s)ds i = 1 K j = 1 N (7) u (x) = x(x ; L) (x) u() = u(l) = (1) can be evaluated analytically and their expressions are given by Integrating eqn.(1) twice we obtain u(x) = G(x s)(s)ds x L () where the Greens function G(x s) isgiven by A ij =(i ; 1)(1j 3 ; 18j ; 4Nj +1j +1jN +4Nj ; 6N ; 8N ; 3)L 4 )=(48KN 4 ) if i<j A ij =(i ; 1 ; K)(1j 3 ; 18j ; 1Nj + +1j +1Nj ; 4N ; 3)L 4 =(48KN 4 ) if i>j G(x s) = ( s (s;l)(x;l) L s x sx(s;l) L x s (3) A ij =[(i ; 1)j (3j ; 8Nj +6N ) +(;i +1+K)(j ; 1) 3 (3j ; 3 ; 4N)]L 4 =(48KN 4 ) +(i ; 1) 3 (i ; 4K ; 1)L 4 =(384K 4 ) if i = j(8) 4 Copyright c 1999 by ASME
5 In order to accommodate for the physical reality that in practice the measurements of u will never be exact or smooth, the deection u is contaminated with some errors, say = (x i ) for i = 1 K,given by a(x) z = u + (9).4 Thus, we have transformed the innite-dimensional illposed problem () into the nite-dimensional problem of nding the vector from the ill-conditioned system of linear equations z u = A (3) In order to deal with the ill-conditioned system of linear eqns (3), the minimization of the discretised version of eqn.() results in a stable numerical solution given by =(A tr A + R) ;1 A tr z (31) where R is the rst-order regularization matrix which has the components R ii = 1 for i = 1 and i = K, R ii = for i = (K ; 1), R i(i+1) = ;1 fori = 1 (K ; 1), R i(i;1) = ;1 for i = K and R ij = otherwise. For a simple typical benchmark example, namely that of a simply supported beam of length L = 1having a deection u(x) = x5 ; x3 + 7x, with the exural rigidity coecient to 6 6 be retrieved given by a(x) = 1, the solution given by (x+1) eqn.(31) has been employed. Also, in order to investigate the stability of the numerical solution, p% = 1% noise has been included in the measured deection z, asgiven by eqn.(9), where the components of the noise vector were generated using the NAG routine G5DDF, as Gaussian random variables with mean zero and standard deviation = p p max j u(x) j = ; (3) The numerically obtained results for a(x) = ;1 (x) given by eqn.(31), when K = N =andp =and =,and when p =1and =5 1 ;5,1 ;5,5 1 ;6 and 1 ;6, are shown in Fig.1 and also for comparison the exact solution is presented. From this gure it can be seen that when no noise is included in the data, i.e. p =, then a simple inversion = A ;1 z retrieves very accurately the exact solution Figure 1. THE NUMERICAL RESULTS ( ) FOR VARIOUS VALUES OF WHEN p% =1%,INCOMPARISON WITH THE EXACT SOLUTION a(x) = 1 (;;;). THE CIRCLES () REPRESENT THE NUMER- (1+x) ICAL RESULTS OBTAINED WHEN p =AND =. However, when noise is included in the data then as becomes very small, say 1 ;6, oscillations start to develop in the numerical solution which becomes unstable. Also as becomes large, say > 1 ;4, then the numerical solution tends to a constant value. However, there is a wide range of values of, say5 1 ;5 <<5 1 ;6 for which the estimate of the solution is stable and reasonably accurate. Nevertheless there is quite a signicant disagreement between the exact and the numerical solutions near the ends of the beam, and in fact the numerical solution may be considered a good estimate of the exact solution only in the range : <x<:8. Finally, it is reported that alternative numerical methods employed by the author, which are based on the singular value decomposition or the mollication method, also failed to produce accurate estimates of the solution near the ends of the xed beam. In a future study it is proposed that more powerful numerical methods which are based on parameterisation and/or constrained minimization procedures should be employed. Also, future work will be concerned with the inverse coecient identication problem associated with the Euler-Bernoulli unsteady-state beam theory. REFERENCES Collins, W.D., Quegan, S., Smith, D. and Taylor, D.A.W., Using Bernstein polynomials to recover engineer- x 5 Copyright c 1999 by ASME
6 ing loads from strain gauge data, Q. Jl Mech. appl. Math. 47, , Ismailov, I. and Muravey, L., Some problems of coecient control for elliptic equations of high order, Inverse and Ill-Posed Problems Conf. (IIPP-96), Moscow, 81 (abstract), Lesnic, D., Elliott, L. and Ingham, D.B., Analysis of coecient identication problems associated to the inverse Euler-Bernoulli beam theory, IMA J. Appl. Math. 6, , Marcellini, P., Identicazione di un coeciente in una equazione dierenziale ordinaria del secondo ordine, Ricerche Mat. 31, 3-43, 198. Tikhonov, A.N. and Arsenin, V.Y., Solutions of Ill- Posed Problems, Winston-Wiley, NewYork, Copyright c 1999 by ASME
3rd Int. Conference on Inverse Problems in Engineering AN INVERSE PROBLEM FOR SLOW VISCOUS INCOMPRESSIBLE FLOWS. University of Leeds
Inverse Problems in Engineering: Theory and Practice 3rd Int. Conference on Inverse Problems in Engineering June 13-18, 1999, Port Ludlow, WA, USA ME06 AN INVERSE PROBLEM FOR SLOW VISCOUS INCOMPRESSIBLE
More informationPARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN. H.T. Banks and Yun Wang. Center for Research in Scientic Computation
PARAMETER IDENTIFICATION IN THE FREQUENCY DOMAIN H.T. Banks and Yun Wang Center for Research in Scientic Computation North Carolina State University Raleigh, NC 7695-805 Revised: March 1993 Abstract In
More informationConjugate gradient-boundary element method for a Cauchy problem in the Lam6 system
Conjugate gradient-boundary element method for a Cauchy problem in the Lam6 system L. Marin, D.N. Hho & D. Lesnic I Department of Applied Mathematics, University of Leeds, UK 'Hanoi Institute of Mathematics,
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 Structural Elements: The truss & beam elements
Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations
More informationTransactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
Boundary element method for an improperly posed problem in unsteady heat conduction D. Lesnic, L. Elliott & D.B. Ingham Department of Applied Mathematical Studies, University of Leeds, Leeds LS2 9JT, UK
More informationMODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione
MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,
More informationDetermination of a space-dependent force function in the one-dimensional wave equation
S.O. Hussein and D. Lesnic / Electronic Journal of Boundary Elements, Vol. 12, No. 1, pp. 1-26 (214) Determination of a space-dependent force function in the one-dimensional wave equation S.O. Hussein
More informationCHAPTER 5. Beam Theory
CHPTER 5. Beam Theory SangJoon Shin School of Mechanical and erospace Engineering Seoul National University ctive eroelasticity and Rotorcraft Lab. 5. The Euler-Bernoulli assumptions One of its dimensions
More informationMARY ANN HORN be decoupled into three wave equations. Thus, we would hope that results, analogous to those available for the wave equation, would hold
Journal of Mathematical Systems, Estimation, and Control Vol. 8, No. 2, 1998, pp. 1{11 c 1998 Birkhauser-Boston Sharp Trace Regularity for the Solutions of the Equations of Dynamic Elasticity Mary Ann
More informationCOMPOSITE PLATE THEORIES
CHAPTER2 COMPOSITE PLATE THEORIES 2.1 GENERAL Analysis of composite plates is usually done based on one of the following the ries. 1. Equivalent single-layer theories a. Classical laminate theory b. Shear
More informationSpurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics
More informationBEAM DEFLECTION THE ELASTIC CURVE
BEAM DEFLECTION Samantha Ramirez THE ELASTIC CURVE The deflection diagram of the longitudinal axis that passes through the centroid of each cross-sectional area of a beam. Supports that apply a moment
More informationDynamic Green Function Solution of Beams Under a Moving Load with Dierent Boundary Conditions
Transaction B: Mechanical Engineering Vol. 16, No. 3, pp. 273{279 c Sharif University of Technology, June 2009 Research Note Dynamic Green Function Solution of Beams Under a Moving Load with Dierent Boundary
More informationA study of the critical condition of a battened column and a frame by classical methods
University of South Florida Scholar Commons Graduate Theses and Dissertations Graduate School 003 A study of the critical condition of a battened column and a frame by classical methods Jamal A.H Bekdache
More informationquantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner
Quantitative Justication of Linearization in Nonlinear Hencky Material Problems 1 Weimin Han and Hong-ci Huang 3 Abstract. The classical linear elasticity theory is based on the assumption that the size
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationAnalysis of Axially Loaded Non-prismatic Beams with General End Restraints Using Differential Quadrature Method
ISBN 978-93-84422-56-1 Proceedings of International Conference on Architecture, Structure and Civil Engineering (ICASCE'15 Antalya (Turkey Sept. 7-8, 2015 pp. 1-7 Analysis of Axially Loaded Non-prismatic
More informationLinear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space
Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................
More informationPresented By: EAS 6939 Aerospace Structural Composites
A Beam Theory for Laminated Composites and Application to Torsion Problems Dr. BhavaniV. Sankar Presented By: Sameer Luthra EAS 6939 Aerospace Structural Composites 1 Introduction Composite beams have
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More information9.1 Introduction to bifurcation of equilibrium and structural
Module 9 Stability and Buckling Readings: BC Ch 14 earning Objectives Understand the basic concept of structural instability and bifurcation of equilibrium. Derive the basic buckling load of beams subject
More information13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs)
13 PDEs on spatially bounded domains: initial boundary value problems (IBVPs) A prototypical problem we will discuss in detail is the 1D diffusion equation u t = Du xx < x < l, t > finite-length rod u(x,
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationCHAPTER -6- BENDING Part -1-
Ishik University / Sulaimani Civil Engineering Department Mechanics of Materials CE 211 CHAPTER -6- BENDING Part -1-1 CHAPTER -6- Bending Outlines of this chapter: 6.1. Chapter Objectives 6.2. Shear and
More informationA Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
More informationForced Vibration Analysis of Timoshenko Beam with Discontinuities by Means of Distributions Jiri Sobotka
21 st International Conference ENGINEERING MECHANICS 2015 Svratka, Czech Republic, May 11 14, 2015 Full Text Paper #018, pp. 45 51 Forced Vibration Analysis of Timoshenko Beam with Discontinuities by Means
More informationexpression that describes these corrections with the accuracy of the order of 4. frame usually connected with extragalactic objects.
RUSSIAN JOURNAL OF EARTH SCIENCES, English Translation, VOL, NO, DECEMBER 998 Russian Edition: JULY 998 On the eects of the inertia ellipsoid triaxiality in the theory of nutation S. M. Molodensky Joint
More informationCHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 2. Discontinuity functions
1. Deflections of Beams and Shafts CHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 1. Integration method. Discontinuity functions 3. Method
More informationDiscontinuous Distributions in Mechanics of Materials
Discontinuous Distributions in Mechanics of Materials J.E. Akin, Rice University 1. Introduction The study of the mechanics of materials continues to change slowly. The student needs to learn about software
More informationThe Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria
ESI The Erwin Schrodinger International Pasteurgasse 6/7 Institute for Mathematical Physics A-1090 Wien, Austria On the Point Spectrum of Dirence Schrodinger Operators Vladimir Buslaev Alexander Fedotov
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationM5 Simple Beam Theory (continued)
M5 Simple Beam Theory (continued) Reading: Crandall, Dahl and Lardner 7.-7.6 In the previous lecture we had reached the point of obtaining 5 equations, 5 unknowns by application of equations of elasticity
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationLECTURE 15 + C+F. = A 11 x 1x1 +2A 12 x 1x2 + A 22 x 2x2 + B 1 x 1 + B 2 x 2. xi y 2 = ~y 2 (x 1 ;x 2 ) x 2 = ~x 2 (y 1 ;y 2 1
LECTURE 5 Characteristics and the Classication of Second Order Linear PDEs Let us now consider the case of a general second order linear PDE in two variables; (5.) where (5.) 0 P i;j A ij xix j + P i,
More informationErkut Erdem. Hacettepe University February 24 th, Linear Diffusion 1. 2 Appendix - The Calculus of Variations 5.
LINEAR DIFFUSION Erkut Erdem Hacettepe University February 24 th, 2012 CONTENTS 1 Linear Diffusion 1 2 Appendix - The Calculus of Variations 5 References 6 1 LINEAR DIFFUSION The linear diffusion (heat)
More informationConstrained Leja points and the numerical solution of the constrained energy problem
Journal of Computational and Applied Mathematics 131 (2001) 427 444 www.elsevier.nl/locate/cam Constrained Leja points and the numerical solution of the constrained energy problem Dan I. Coroian, Peter
More information7 TRANSVERSE SHEAR transverse shear stress longitudinal shear stresses
7 TRANSVERSE SHEAR Before we develop a relationship that describes the shear-stress distribution over the cross section of a beam, we will make some preliminary remarks regarding the way shear acts within
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationVIBRATION ANALYSIS OF EULER AND TIMOSHENKO BEAMS USING DIFFERENTIAL TRANSFORMATION METHOD
VIBRATION ANALYSIS OF EULER AND TIMOSHENKO BEAMS USING DIFFERENTIAL TRANSFORMATION METHOD Dona Varghese 1, M.G Rajendran 2 1 P G student, School of Civil Engineering, 2 Professor, School of Civil Engineering
More informationXIAOFENG REN. equation by minimizing the corresponding functional near some approximate
MULTI-LAYER LOCAL MINIMUM SOLUTIONS OF THE BISTABLE EQUATION IN AN INFINITE TUBE XIAOFENG REN Abstract. We construct local minimum solutions of the semilinear bistable equation by minimizing the corresponding
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationTable 1: BEM as a solution method for a BVP dierential formulation FDM BVP integral formulation FEM boundary integral formulation BEM local view is ad
Chapter 1 Introduction to Boundary element Method - 1D Example For reference: Hong-Ki Hong and Jeng-Tzong Chen, Boundary Element Method, Chapter 1 Introduction to Boundary element Method - 1D Example,
More informationCIV-E1060 Engineering Computation and Simulation Examination, December 12, 2017 / Niiranen
CIV-E16 Engineering Computation and Simulation Examination, December 12, 217 / Niiranen This examination consists of 3 problems rated by the standard scale 1...6. Problem 1 Let us consider a long and tall
More informationAN INEQUALITY FOR THE NORM OF A POLYNOMIAL FACTOR IGOR E. PRITSKER. (Communicated by Albert Baernstein II)
PROCDINGS OF TH AMRICAN MATHMATICAL SOCITY Volume 9, Number 8, Pages 83{9 S -9939()588-4 Article electronically published on November 3, AN INQUALITY FOR TH NORM OF A POLYNOMIAL FACTOR IGOR. PRITSKR (Communicated
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More informationonly nite eigenvalues. This is an extension of earlier results from [2]. Then we concentrate on the Riccati equation appearing in H 2 and linear quadr
The discrete algebraic Riccati equation and linear matrix inequality nton. Stoorvogel y Department of Mathematics and Computing Science Eindhoven Univ. of Technology P.O. ox 53, 56 M Eindhoven The Netherlands
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationTransactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X
An integral equation formulation of the coupled vibrations of uniform Timoshenko beams Masa. Tanaka & A. N. Bercin Department of Mechanical Systems Engineering, Shinshu University 500 Wakasato, Nagano
More informationAbstract. Jacobi curves are far going generalizations of the spaces of \Jacobi
Principal Invariants of Jacobi Curves Andrei Agrachev 1 and Igor Zelenko 2 1 S.I.S.S.A., Via Beirut 2-4, 34013 Trieste, Italy and Steklov Mathematical Institute, ul. Gubkina 8, 117966 Moscow, Russia; email:
More informationAn Optimal Finite Element Mesh for Elastostatic Structural. Analysis Problems. P.K. Jimack. University of Leeds. Leeds LS2 9JT, UK.
An Optimal Finite Element Mesh for Elastostatic Structural Analysis Problems P.K. Jimack School of Computer Studies University of Leeds Leeds LS2 9JT, UK Abstract This paper investigates the adaptive solution
More informationChapter 2 Examples of Optimization of Discrete Parameter Systems
Chapter Examples of Optimization of Discrete Parameter Systems The following chapter gives some examples of the general optimization problem (SO) introduced in the previous chapter. They all concern the
More informationVibration of a Euler Bernoulli uniform beam carrying a rigid body at each end
Vibration of a Euler Bernoulli uniform beam carrying a rigid body at each end S. Naguleswaran Department of Mechanical Engineering, University of Canterbury, Christchurch, New Zealand E-Mail: sivapatham.naguleswaran@canterbury.ac.nz
More informationENGI 4430 PDEs - d Alembert Solutions Page 11.01
ENGI 4430 PDEs - d Alembert Solutions Page 11.01 11. Partial Differential Equations Partial differential equations (PDEs) are equations involving functions of more than one variable and their partial derivatives
More informationThe Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 29, 1999 Abstract This paper deals with the construction of wavelet approximatio
The Mortar Wavelet Method Silvia Bertoluzza Valerie Perrier y October 9, 1999 Abstract This paper deals with the construction of wavelet approximation spaces, in the framework of the Mortar method. We
More informationHomework No. 1 MAE/CE 459/559 John A. Gilbert, Ph.D. Fall 2004
Homework No. 1 MAE/CE 459/559 John A. Gilbert, Ph.D. 1. A beam is loaded as shown. The dimensions of the cross section appear in the insert. the figure. Draw a complete free body diagram showing an equivalent
More informationXIAO-CHUAN CAI AND MAKSYMILIAN DRYJA. strongly elliptic equations discretized by the nite element methods.
Contemporary Mathematics Volume 00, 0000 Domain Decomposition Methods for Monotone Nonlinear Elliptic Problems XIAO-CHUAN CAI AND MAKSYMILIAN DRYJA Abstract. In this paper, we study several overlapping
More informationCLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE
CLASSIFICATION AND PRINCIPLE OF SUPERPOSITION FOR SECOND ORDER LINEAR PDE 1. Linear Partial Differential Equations A partial differential equation (PDE) is an equation, for an unknown function u, that
More informationUNIVERSITY OF SASKATCHEWAN ME MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich
UNIVERSITY OF SASKATCHEWAN ME 313.3 MECHANICS OF MATERIALS I FINAL EXAM DECEMBER 13, 2008 Professor A. Dolovich A CLOSED BOOK EXAMINATION TIME: 3 HOURS For Marker s Use Only LAST NAME (printed): FIRST
More informationTable of Contents. Preface... 13
Table of Contents Preface... 13 Chapter 1. Vibrations of Continuous Elastic Solid Media... 17 1.1. Objective of the chapter... 17 1.2. Equations of motion and boundary conditions of continuous media...
More informationRATIONAL POINTS ON CURVES. Contents
RATIONAL POINTS ON CURVES BLANCA VIÑA PATIÑO Contents 1. Introduction 1 2. Algebraic Curves 2 3. Genus 0 3 4. Genus 1 7 4.1. Group of E(Q) 7 4.2. Mordell-Weil Theorem 8 5. Genus 2 10 6. Uniformity of Rational
More informationMechanics of Solids notes
Mechanics of Solids notes 1 UNIT II Pure Bending Loading restrictions: As we are aware of the fact internal reactions developed on any cross-section of a beam may consists of a resultant normal force,
More informationAlgorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations
dmitry.lyakhov@kaust.edu.sa Symbolic Computations May 4, 2018 1 / 25 Algorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations Dmitry A. Lyakhov 1 1 Computational Sciences
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More informationLinearized methods for ordinary di erential equations
Applied Mathematics and Computation 104 (1999) 109±19 www.elsevier.nl/locate/amc Linearized methods for ordinary di erential equations J.I. Ramos 1 Departamento de Lenguajes y Ciencias de la Computacion,
More informationBOUNDARY CONDITION IDENTIFICATION BY SOLVING CHARACTERISTIC EQUATIONS
Journal of Sound and
More informationSpectrum and Exact Controllability of a Hybrid System of Elasticity.
Spectrum and Exact Controllability of a Hybrid System of Elasticity. D. Mercier, January 16, 28 Abstract We consider the exact controllability of a hybrid system consisting of an elastic beam, clamped
More informationMonte Carlo Methods for Statistical Inference: Variance Reduction Techniques
Monte Carlo Methods for Statistical Inference: Variance Reduction Techniques Hung Chen hchen@math.ntu.edu.tw Department of Mathematics National Taiwan University 3rd March 2004 Meet at NS 104 On Wednesday
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationName :. Roll No. :... Invigilator s Signature :.. CS/B.TECH (CE-NEW)/SEM-3/CE-301/ SOLID MECHANICS
Name :. Roll No. :..... Invigilator s Signature :.. 2011 SOLID MECHANICS Time Allotted : 3 Hours Full Marks : 70 The figures in the margin indicate full marks. Candidates are required to give their answers
More information6. Bending CHAPTER OBJECTIVES
CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where
More informationChapter 4 Deflection and Stiffness
Chapter 4 Deflection and Stiffness Asst. Prof. Dr. Supakit Rooppakhun Chapter Outline Deflection and Stiffness 4-1 Spring Rates 4-2 Tension, Compression, and Torsion 4-3 Deflection Due to Bending 4-4 Beam
More informationThe Homotopy Perturbation Method for free vibration analysis of beam on elastic foundation
Shiraz University of Technology From the SelectedWorks of Habibolla Latifizadeh 2011 The Homotopy Perturbation Method for free vibration analysis of beam on elastic foundation Habibolla Latifizadeh, Shiraz
More informationContents Introduction Energy nalysis of Convection Equation 4 3 Semi-discrete nalysis 6 4 Fully Discrete nalysis 7 4. Two-stage time discretization...
Energy Stability nalysis of Multi-Step Methods on Unstructured Meshes by Michael Giles CFDL-TR-87- March 987 This research was supported by ir Force Oce of Scientic Research contract F4960-78-C-0084, supervised
More informationAircraft Structures Kirchhoff-Love Plates
University of Liège erospace & Mechanical Engineering ircraft Structures Kirchhoff-Love Plates Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/ Chemin
More information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationAnalysis of Infinitely Long Euler Bernoulli Beam on Two Parameter Elastic Foundation: Case of Point Load
Analysis of Infinitely Long Euler Bernoulli Beam on Two Parameter Elastic Foundation: Case of Point Load Dr. Mama, Benjamin Okwudili Department of Civil Engineering, University of Nigeria, Nsukka. benjamin.mama@unn.edu.ng
More informationGeometry-dependent MITC method for a 2-node iso-beam element
Structural Engineering and Mechanics, Vol. 9, No. (8) 3-3 Geometry-dependent MITC method for a -node iso-beam element Phill-Seung Lee Samsung Heavy Industries, Seocho, Seoul 37-857, Korea Hyu-Chun Noh
More informationPURE BENDING. If a simply supported beam carries two point loads of 10 kn as shown in the following figure, pure bending occurs at segment BC.
BENDING STRESS The effect of a bending moment applied to a cross-section of a beam is to induce a state of stress across that section. These stresses are known as bending stresses and they act normally
More information7. Hierarchical modeling examples
7. Hierarchical modeling examples The objective of this chapter is to apply the hierarchical modeling approach discussed in Chapter 1 to three selected problems using the mathematical models studied in
More informationMath 5587 Lecture 2. Jeff Calder. August 31, Initial/boundary conditions and well-posedness
Math 5587 Lecture 2 Jeff Calder August 31, 2016 1 Initial/boundary conditions and well-posedness 1.1 ODE vs PDE Recall that the general solutions of ODEs involve a number of arbitrary constants. Example
More informationModule 4 : Deflection of Structures Lecture 4 : Strain Energy Method
Module 4 : Deflection of Structures Lecture 4 : Strain Energy Method Objectives In this course you will learn the following Deflection by strain energy method. Evaluation of strain energy in member under
More informationSEMIGROUP APPROACH FOR PARTIAL DIFFERENTIAL EQUATIONS OF EVOLUTION
SEMIGROUP APPROACH FOR PARTIAL DIFFERENTIAL EQUATIONS OF EVOLUTION Istanbul Kemerburgaz University Istanbul Analysis Seminars 24 October 2014 Sabanc University Karaköy Communication Center 1 2 3 4 5 u(x,
More informationChapter 3 Higher Order Linear ODEs
Chapter 3 Higher Order Linear ODEs Advanced Engineering Mathematics Wei-Ta Chu National Chung Cheng University wtchu@cs.ccu.edu.tw 1 2 3.1 Homogeneous Linear ODEs 3 Homogeneous Linear ODEs An ODE is of
More informationCubic Splines. Antony Jameson. Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305
Cubic Splines Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305 1 References on splines 1. J. H. Ahlberg, E. N. Nilson, J. H. Walsh. Theory of
More informationLie Groups for 2D and 3D Transformations
Lie Groups for 2D and 3D Transformations Ethan Eade Updated May 20, 2017 * 1 Introduction This document derives useful formulae for working with the Lie groups that represent transformations in 2D and
More informationFor an imposed stress history consisting of a rapidly applied step-function jump in
Problem 2 (20 points) MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0239 2.002 MECHANICS AND MATERIALS II SOLUTION for QUIZ NO. October 5, 2003 For
More information436 A. Barani and G.H. Rahimi assessment models have been employed to investigate the LBB of cracked pipes that are not for combined load [8]. Yun-Jae
Scientia Iranica, Vol. 4, No. 5, pp 435{44 c Sharif University of Technology, October 27 Approximate Method for Evaluation of the J-Integral for Circumferentially Semi-Elliptical-Cracked Pipes Subjected
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationTwo-parameter regularization method for determining the heat source
Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 13, Number 8 (017), pp. 3937-3950 Research India Publications http://www.ripublication.com Two-parameter regularization method for
More information1 Solutions to selected problems
1 Solutions to selected problems 1. Let A B R n. Show that int A int B but in general bd A bd B. Solution. Let x int A. Then there is ɛ > 0 such that B ɛ (x) A B. This shows x int B. If A = [0, 1] and
More informationSolution of Nonlinear Equations
Solution of Nonlinear Equations In many engineering applications, there are cases when one needs to solve nonlinear algebraic or trigonometric equations or set of equations. These are also common in Civil
More informationENGI Second Order Linear ODEs Page Second Order Linear Ordinary Differential Equations
ENGI 344 - Second Order Linear ODEs age -01. Second Order Linear Ordinary Differential Equations The general second order linear ordinary differential equation is of the form d y dy x Q x y Rx dx dx Of
More information4. Mathematical models used in engineering structural analysis
4. Mathematical models used in engineering structural analysis In this chapter we pursue a formidable task to present the most important mathematical models in structural mechanics. In order to best situate
More informationFirst-Order Solutions for the Buckling Loads of Euler-Bernoulli Weakened Columns
First-Order Solutions for the Buckling Loads of Euler-Bernoulli Weakened Columns J. A. Loya ; G. Vadillo 2 ; and J. Fernández-Sáez 3 Abstract: In this work, closed-form expressions for the buckling loads
More informationSlender Structures Load carrying principles
Slender Structures Load carrying principles Basic cases: Extension, Shear, Torsion, Cable Bending (Euler) v017-1 Hans Welleman 1 Content (preliminary schedule) Basic cases Extension, shear, torsion, cable
More informationu, w w(x) ε v(x) u(x)
Chapter One-Dimensional Finite Element Methods.1 Introduction The piecewise-linear Galerkin nite element method of Chapter 1 can be extended in several directions. The most important of these is multi-dimensional
More informationThe following can also be obtained from this WWW address: the papers [8, 9], more examples, comments on the implementation and a short description of
An algorithm for computing the Weierstrass normal form Mark van Hoeij Department of mathematics University of Nijmegen 6525 ED Nijmegen The Netherlands e-mail: hoeij@sci.kun.nl April 9, 1995 Abstract This
More informationOur goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always
October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.
More information