Boundary element method, finite element method and thefluxspline method: a performance comparison for scalar potential problems
|
|
- Drusilla Warner
- 5 years ago
- Views:
Transcription
1 Boundary element method, finite element method and thefluxspline method: a performance comparison for scalar potential problems A. BulcSo, C. F. Loeffler, P. C. Oliveira Departamento de Engenharia Mecdnica, Universidade Federal do Espirito Santo (UFES), av. Fernando Ferrari s/n, Goiabeiras, Vitoria, Espirito Santo, Brazil carlosloeffler@hotmail com Abstract This work presents comparisons among the numeric results obtained by the Boundary Element Method, Finite Element Method and Flux Spline Method applied to two-dimensional scalar potential problems. Cases are stationary and governed by Laplace and Poisson equations. The main purpose of this paper is to compare the accuracy of mentioned methods in computational simulations. So, the applications were selected to analyze numerical behavior when some special physical and geometrical aspects are involved, such as: discontinuities in the boundary conditions, actions concentrated on specific part on the domain, non-homogeneity in constitutive properties etc. 1 Introduction The high level of modeling of physical problems solved by the current engineering results of the use of powerful numeric methods of solution, that they use modern approaches mathematically consistent. Such methods allow to substitute the continuous domain, represented by partial differential equations, for equivalent discrete models, expressed for algebraic equations, with easy computational treatment.
2 260 Boundary Elements They are pioneers in this approach the Finite Difference Finite (DFM) and the Finite Element Method (FEM), that still today are the most important discrete numerical techniques. Recently, two other similar methods were introduced in engineering application areas: The Boundary Element Method (BEM) and the Flux Spline Method (FSM). It is verified that the specialized literature doesn't register sufficient number of works that process acting comparisons among those methods, involving precision, cost, flexibility and other aspects, so that the engineers may choose correctly the better or more suitable numerical technique for its purposes. 2 Escalar Potential Equation Consider stationary problems in two dimensions, where the physical field is scalar. Being assumed homogeneity, isotropy and physical linearity, the field can be quite expressed by the following equation: K (1) This expression is denominated by Poisson Equation, usual in the representation of stationary problems where there are presence of sources, flows or field forces inside the domain, represented by the function S. When the function S is null occurs the simplest case in the scalar field theory. The equation like that is named as Laplace Equation. 3 The Boundary Element Method Consider the Equation of Poisson (eqn 1). The formulation of BEM for such problem begins for the establishment of the Inverse integral form associated to Poisson's equation, given for: JYn,) >* -(jx'n, + JY'n^Jir (2) In former equation Q is the domain of the system, T is the boundary, ( > is the Fundamental Solution, corresponding the answer of a correlate problem defined by Poisson Equation, where it exist a singular source applied in a particular point,. This special source is represented by the Delta de Dirac function. In the point of view of the Weighted Residual Method, ( >* acts as an efficient weighted function to allow the good representation of concentrated actions and discontinuity in the behavior of the involved variables (Brebbia[1]). It was still defined other new variables: JX = K^*, JY=K^*, JX*=^ and JY* =K^. 0) 5x 3y 5x 9y
3 Boundary Elements 261 Being n% and ny the normal vectors in the directions coordinates x and y, respectively. Using the properties of the Delta of Dirac, it is possible to rewrite the equation (2) in the form: n (4) where, for subject of simplicity: J = JXn% + JYny; J = JX n^ + JY ny. The next step is to operate the right side of the equation (4), in order to also express it in the form of boundary integrals. Due to generality, the Dual Reciprocity technique is the most effective procedure. So, it is possible to express the following sentence: j-i (5) In this former equation w\ p* and W form a group of auxiliary functions, composing a similar mathematical structure to the right side of the equation (2). The function p-* and its primitive w^ are arbitrary (Partridge et al. [2]). The next step consists of the boundary discretization, with its division in elements, which <(> and normal derivative J they are approximate. It is developed, then, a system of equations starting from the location of the coincident points with the nodes of the elements (collocation method). The same procedure is made to the j points. A final system in matricial form results: H( >-GJ = S (6) 4 Flux Spline Method For the formulation of FSM it considers the Poisson Equation as follows: "SIX 9JY,»j *»j - O For a elementary volume control V^ (according to exhibition the illustration 1), the integration of governing equation produces: K«*JJ - J*i>Yj + KM - JYjAX,]= S,jAX,AY, W JX, Figure 1 - Elementary control volume.
4 262 Boundary Elements In the discretization procedure used is assumed that inside each control volume the flows J can vary linearly along each direction. It is interesting to use elementary adimensional coordinates, in the form: J =A-ri + B, where * = &_ZJo_ and B = JQ (9) ATI This technique of discretization is denominated of Flux-Spline and already in that point, due to the fact of the angular coefficient A to absorb source terms sources and multidimensional influences, it reveals the largest effectiveness when compared to the traditional procedures that, assuming constant flow, it takes to the well-known central difference scheme. It is necessary to take account the hypothesis of lineal flow, equation (9), and to process integration for the control volume in subject. With the condition that the dependent variable (() is stored in the center of the control volume r, it is possible to obtain the quadratic shape potential function in the considered direction: * *MiHl <"» The unidimensional procedure is generalized being considered the flow expressed of similar way in other directions like it was shown for the direction T. The joining among adjacent control volumes is made being imposed the continuity of the function (J>, and of the flows in the interface of the same ones. This way, being applied the continuity condition of ( ) in X and Y, in a volume control i j is obtained the following discretizated equation for (() : j+, - AX, + DJYy AX, + DJX,,,. AY^ + DJXy AY. + ^ - JX,,,j). AY. Values of coefficients that appear in last equation can be obtained in Oliveira [3 ]. The flows in each direction are calculated being assumed initially A that J are null. The system of algebraic equations for (() is, then, solved in iterative form through the solution of tri-diagonal matrices in each direction coordinate X and Y (Line-by-Line), being used Thomas' algorithm (TDMA). With that result it returns to the calculation of the flows to the convergence. 5 The Finite Element Method There are several formulations of the Method of the Finite Elements, but the most common and important way begins with the establishment of Weak integral form associated to the Poisson governing equation, given for:
5 Boundary Elements 263 _K f * * + * * dq + KfodT = ("SW (12) J dx ox dy dy J J QL ' -U r a Under the point of view of several important approaches, <j> it can take interpretations as weighted function, admissible function or virtual potential, according to the case (Zienkiewicz, [4]). The FEM establishes the division of the domain O in discrete subdomains (the finite elements) whose points that define them are denominated nodal points. Usually it is assumed the basic variable $ and the auxiliary function f are approximated by the same polynomial interpolation function, whose distribution inside the element is given with reference to the values in the nodal points <Dj and O J, in the following way: Derivatives of<dj and Oj are calculated through interpolation functions in agreement with the coordinate directions. Making the convenient substitutions the integral expression at finite element level, and being considered that the nodal values are not integration variables, the previous expression, equation (12), can be arranged, placing the auxiliary parameters <X> in evidence in such a way that: -O^k Jaj,a%,,dQ + k = 0 (14) It is convenient to define the equation (14) in an appropriate matricial form, where thefirstintegral is expressed for K\ as being an elementary matrix of properties; the second integral is expressed for f \ being a vector of nodal forces related to boundary interaction of the finite element; and the expressed third integral p% represents a vector of nodal forces relative to the external domain actions. The global arrangement of the system is made with base in the compatibility among the nodal potential values, resulting in a matricial system, given for: K<D = f-p (*5) 6 Examples 6.1 First example The first example is governed by equation of Laplace. It consists of a seepage problem, which its characteristics are given in the illustration (2):
6 264 Boundary Elements Figure 2 - Seepage problem: geometry and boundary conditions. In this problem in subject is particularly difficult to get a good computational response due to the discontinuity in the value of the flow in the superior part of the porous medium. Next results are shown for the flow along the superior boundary obtained by BEM and FSM simulations. BEM used elements with constant interpolation. The results for both methods are shown in the figure 3(a) and 3(b): B.E ^t* 0.0 F V o.o 4*,-*, r-^-f 0.0 -< 3.5 i &B.E.... ": 2.5 * _ 10 * * ^ * BE Figure 3a - Boundary Elements Method results j " _ (1 FV ~ 32DJJ/I _ * - # Figure 3b - Flux Spline Method results. ;:; ; ' ; / ;» The agreement among the values obtained by both methods is quite good, but undoubtedly a small advantage for BEM can be admitted, not just for the best reception of the singularity, but also for the use of less refined meshes and more simpler interpolation shape.
7 Boundary Elements Second example In this second example the equation of Laplace is still approached, confronting the performance of the BEM front to FEM, in a problem of heat conduction. In the figure 4 are shown the characteristics of problem. In this example both FEM and BEM methods used isoparametric quadratic elements in their meshes. It is stood out that for FEM the program was used Ansys 5.4, with the following elements: plane35 (triangular element of 6 nodes) and plane?? (square element of 8 nodes). Figure 4 - Physical representation and different discretization used. For comparison among the methods are shown graphs, in logarithmic scale, contends the averages of the percentual error in the numeric response with reference of analytic solution for all discrete points. In vertical axis are plotted the value of global percentual error and in the horizontal axis are plotted the number of elements =- = ,,-rr, 0.00 \ BEM FEM Figure5- (a) Potential, nodal points on the boundary; (b) Potential, nodal points inside the domain; (c) Potential Derivative (radial direction), nodal points inside the domain. It is wide the advantage the quality of the results of the BEM in that analysis, particularly in respect to the calculation of values inside the domain. Really, that excellent behavior of the BEM is general in simulation of problems governed by the Equation of Laplace.
8 266 Boundary Elements 6.3 Third example In the third example, the numeric difficulty involved it consists of the discontinuity of the potential prescribed in the horizontal edge of the semicircle shown in the illustration.6. In the point where there is discontinuity of the value of the potential, for prescription of the boundary conditions in both methods the mean value of the potentials was adopted prescribed in this face. A Y Figure 6 - Physical representation of the problem. In the figure 7 are depicted the results. Shape quadratic functions are used in BEM and FEM discretization. Once again it is noticed a wide advantage in the results obtained with to BEM front to FEM in that class of problems. 100^, _ =-- = - A.~=^ = 0.10 ^ ^ - lt_ -=~^ = \ BEM FEM Figure 7 -(a) Potential, nodal points on the boundary (b) Potential Derivative (normal direction), nodal points on the boundary; (c) Potential, nodal points inside the domain. 6.4 Fourth example It consists of the Poisson's problem. A square membrane, totally fastened in its boards, is loaded in a central section of its surface, as shown in figure 8. In that case they are compared the performance of the BEM in comparison to FSM. The results are shown in the figure (9), where the percentual error is calculated in comparison the analytic solution. M j\ a = 9b.l. b = 3 b.l. & P = 1 N/b.l. T(pre-tragao) 11 N/b.l. Figure 8- Square membrane.
9 Boundary Elements 267 The discretizations used were, respectively, 14, 42 and 126 constant boundary elements and 9, 81 and 729 finite volumes. The scale logarithmic exhibits the convergence of the methods with the refinement. The best performance of FSM is just verified when the mesh is not very refined. Such behavior modifies with the increase of the number of entities, where it is noticed a better performance of the BEM. It should be considered that in this example were used constant boundary elements, while the interpolant function of FSM is quadratic. More ahead an analysis is presented in which both methods are used with elements of same order Number of Elements Figure 9 - Error in displacement of central point of the membrane. 6.5 Fifth Example In this problem governed by the equation of Poisson, is applied a body action equal to 10 units (force/u.l.^) constant along the whole physical domain, and the material property is unitary. Figure 10 shows the aspect of the problem. Figure 10 - Different discretizations used in circular membrane problem. In following illustration, are presented the graphs contends the average of the percentual error obtained for this analysis.
10 268 Boundary Elements i :.. -r ^ - \ - ~ - _ n:^-^ ,,,-n BEM FEM Figure 11 - (a) Potential Derivative of the Potential (normal direction), nodal points localized on the boundary; (b) Potential, nodal points inside the domain; (c) Potential Derivative (radial direction), nodal points inside the domain. Results shown in figure 11 also show the better performance of BEM front to FEM in this problem, in despite of the difference among behavior of two methods is not very considerable like in Laplace's simulations. Indeed, here it is used Dual Reciprocity procedure, which mathematically introduces approximation in domain action representation. However, considering the good results, they mean that the error in Dual Reciprocity model have small significance and the Boundary Element approach is not affected negatively. 6.6 Sixth example The physical domain of this analysis is a square shown in figure 12. It subjects the null essential boundary conditions along its boards. There is an application of a domain action constant in all its extension. p = 10 For^a/u.l.^ X u = < 10 u.c. Figure 12 - Representation of the problem and different discretization. Figure 13 presents the graphs contends the results. In this case were used linear boundary elements, while quadratic finite elements are used in FEM analysis. Once more the BEM models show better results that FEM ones. Here the good performance of BEM probably results of its capacity to represent behavior of functions in corners.
11 Boundary Elements ^ ^ : V _ ^ n T, BEM FEM Figure 13 -(a) Potential Derivative (normal direction), nodal points on the boundary; (b) Potential, nodal points inside the domain; (c) Potential Derivative (module), nodal points inside the domain. 6.7 Seventh example This problem represents a physical situation where heat generation rods of square cross section (material k2) with heat generation are embedded centrally in a slab (material kl), the surfaces of which are maintained with a constant temperature [5]. In the present simulation the value of conductivity ratio kl/k2 is 100. It must be noted that problems of this kind have numerical difficulties by finite difference approach. The comparison presented here will be done among BEM and FSM. For the first method, it is adopted 40 quadratic boundary elements and 20 poles in source region. For Flux-Spline was used finite volumes. The high level of refinement of FSM allows to take its response like a reference solution. The physical problem is shown in figure 14. u=0 u=0 q=0 u=0 u=0 q=0 Figure 14 - Representation of the seventh example. An excellent agreement it is observed in figure 15 among the temperature solution in both numerical methods. In spite of the high difference among the physical properties of the sub-areas, the numerical behavior was similar. It must be noticed that Dual Reciprocity performance was very good, because this problem is not a typical case for boundary technique approach. It should be noticed that for the FSM, increasing the relation among the physical properties of the sub-regions, an increase of the number of necessary iterations is made to reach the preset stop criteria. Because of this, in such cases, the distribution of the profiles of temperatures and flows become more intricate, requesting like this a larger number of iterations.
12 270 Boundary Elements Figure 15 - Representation of temperature inside the physical domain. Lines mean isothermal paths. Left side: FSM solution;right side:bem solution. Conclusions The problems presented in this work showed the excellence of the BEM results and its superiority in comparison to the other methods in that problem category. Such problems don't possess any private difficulty that turn them inadequate to FEM and FSM. In fact, what happens it is a larger consistency of the discrete model generated by the BEM in these cases. It should be noticed that even in the cases where domain actions exist, the use of the BEM, with the aproximation of the Dual Reciprocity procedure, it produced superior results. References 1. Brebbia, C.A.; Telles, J.C.F. & Wrobel, L.C. Boundary Element Techniques, Theory and Applications in Engineering, Springer-Verlag, New York, Partridge, P.W., Brebbia, C.A. & Wrobel, L.C. The Dual Reciprocity, Boundary Element Method, Computacional Mechanics Publications, London, Oliveira, P. C. Esquema Flux Spline Aplicado em Cavidades Abertas com Convercsao Natural, UNICAMP, Tese de Doutorado, S&o Paulo., Zienkiewicz, O.C. e Taylor, R.L.; The Finite Element Method, Volume I, Basic Formulation and Linear Problems, McGraw-Hill Book Company, London, 1988, 5. Kelkar, K. M., Iterative Method for the Numerical Prediction of Heat Transfer in Problems Involving Large Differences in Thermal Conductivities, Numerical Heat Transfer, Part B, vol 17, pp , 1990.
AN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS
th Pan-American Congress of Applied Mechanics January 04-08, 00, Foz do Iguaçu, PR, Brazil AN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS Otávio Augusto Alves
More informationTransactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
Model analysis of plates using the dual reciprocity boundary element method T.W. Davies & F.A. Moslehy Department of Mechanical and Aerospace Engineering, University of Central Florida, Orlando, Florida,
More informationUniversity of Hertfordshire Department of Mathematics. Study on the Dual Reciprocity Boundary Element Method
University of Hertfordshire Department of Mathematics Study on the Dual Reciprocity Boundary Element Method Wattana Toutip Technical Report 3 July 999 Preface The boundary Element method (BEM) is now recognised
More informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Plate-soil elastodynamic coupling using analysis S.F.A. Baretto, H.B. Coda, W.S. Venturini Sao Carlos School of Engineering, University ofsao Paulo, Sao Carlos - SP, Brazil BEM Abstract The aim of this
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
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 informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More informationTransactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
Analysis of complexfluid-structureinteraction problems using a triangular boundary element for the fluid dynamics A. Hopf, G. Hailfmger, R. Krieg Forschungszentrum Karlsruhe GmbH, Institutfur D- 7^027
More informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More informationContents. Prologue Introduction. Classical Approximation... 19
Contents Prologue........................................................................ 15 1 Introduction. Classical Approximation.................................. 19 1.1 Introduction................................................................
More informationA simple Galerkin boundary element method for three-dimensional crack problems in functionally graded materials
Materials Science Forum Vols. 492-493 (2005) pp 367-372 Online available since 2005/Aug/15 at www.scientific.net (2005) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/msf.492-493.367
More informationAn alternative multi-region BEM technique for layered soil problems
An alternative multi-region BM technique for layered soil problems D.B. Ribeiro & J.B. Paiva Structural ngineering Department, São Carlos ngineering School, University of São Paulo, Brazil. Abstract Different
More informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationSTUDY OF THE BOUNDARY ELEMENT METHOD FOR THE DIRECT PROBLEM OF ELECTRICAL IMPEDANCE TOMOGRAPHY
STUDY OF THE BOUNDARY ELEMENT METHOD FOR THE DIRECT PROBLEM OF ELECTRICAL IMPEDANCE TOMOGRAPHY Vanessa Rolnik, e-mail: vanessarolnik@ffclrp.usp.br Olavo H. Menin, e-mail: olavohmenin@terra.com.br Departamento
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationIterative schemes for the solution of systems of equations arising from the DRM in multidomains
CHAPTER 7 Iterative schemes for the solution of systems of equations arising from the DRM in multidomains M.I. Portapila 1 &H.Power 2 1 Wessex Institute of Technology, UK. 2 Department of Mechanical Engineering,
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationNumerical Methods. Equations and Partial Fractions. Jaesung Lee
Numerical Methods Equations and Partial Fractions Jaesung Lee Solving linear equations Solving linear equations Introduction Many problems in engineering reduce to the solution of an equation or a set
More informationFinite Elements for Thermal Analysis
Chapter 4 Finite Elements for Thermal Analysis 4.1 Governing Equations for the Thermal Problem CONSTITUTIVE EQUATIONS Introducing the temperature field T(x i,t ), the heat flux vector q(x i,t) q n (x i,t)
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationNumerical Analysis of Electromagnetic Fields
Pei-bai Zhou Numerical Analysis of Electromagnetic Fields With 157 Figures Springer-Verlag Berlin Heidelberg New York London Paris Tokyo Hong Kong Barcelona Budapest Contents Part 1 Universal Concepts
More informationAppendix 1. Numerical Integration: Regular Integrals
The Dual Reciprocity 269 Appendix 1 Numerical Integration: Regular Integrals 1 Introduction The numerical integration formulae given in this Appendix are used for integration over boundary elements which
More informationA new meshless method for steady-state heat conduction problems in anisotropic and inhomogeneous media
Archive of Applied Mechanics 74 25 563--579 Springer-Verlag 25 DOI 1.17/s419-5-375-8 A new meshless method for steady-state heat conduction problems in anisotropic and inhomogeneous media H. Wang, Q.-H.
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationDISPENSA FEM in MSC. Nastran
DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation
More informationMEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy
MEEN 673 ********************************* NONLINEAR FINITE ELEMENT ANALYSIS ********************************* J. N. Reddy e-mail: jnreddy@tamu.edu Texas A&M University College Station, TX 77843-3123 USA
More informationChapter 3. Loop and Cut-set Analysis
Chapter 3. Loop and Cut-set Analysis By: FARHAD FARADJI, Ph.D. Assistant Professor, Electrical Engineering, K.N. Toosi University of Technology http://wp.kntu.ac.ir/faradji/electriccircuits2.htm References:
More informationA simplified boundary element approach for heat transfer in aerospace applications LA. Ali" & J.J. Rencis&
A simplified boundary element approach for heat transfer in aerospace applications LA. Ali" & J.J. Rencis& M.4 ^ Mechanical Engineering Department, Worcester Polytechnic Institute, Worcester, Massachusetts
More informationEach is equal to CP 1 minus one point, which is the origin of the other: (C =) U 1 = CP 1 the line λ (1, 0) U 0
Algebraic Curves/Fall 2015 Aaron Bertram 1. Introduction. What is a complex curve? (Geometry) It s a Riemann surface, that is, a compact oriented twodimensional real manifold Σ with a complex structure.
More informationTransactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
An alternative boundary element formulation for plate bending analysis J.B. Paiva, L.O. Neto Structures Department, Sao Carlos Civil Engineering School, f, BrazzY Abstract This work presents an alternative
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #2 April 14, 2014 Name: SOLUTION ID#: PS1.:
More informationIntroduction to Finite Element Method
Introduction to Finite Element Method Dr. Rakesh K Kapania Aerospace and Ocean Engineering Department Virginia Polytechnic Institute and State University, Blacksburg, VA AOE 524, Vehicle Structures Summer,
More informationGeneralized Eigenvectors and Jordan Form
Generalized Eigenvectors and Jordan Form We have seen that an n n matrix A is diagonalizable precisely when the dimensions of its eigenspaces sum to n. So if A is not diagonalizable, there is at least
More informationAdvanced Vibrations. Distributed-Parameter Systems: Approximate Methods Lecture 20. By: H. Ahmadian
Advanced Vibrations Distributed-Parameter Systems: Approximate Methods Lecture 20 By: H. Ahmadian ahmadian@iust.ac.ir Distributed-Parameter Systems: Approximate Methods Rayleigh's Principle The Rayleigh-Ritz
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 informationPlate analysis using classical or Reissner- Mindlin theories
Plate analysis using classical or Reissner- Mindlin theories L. Palermo Jr. Faculty of Civil Engineering, State Universiv at Campinas, Brazil Abstract Plates can be solved with classical or Reissner-Mindlin
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationChapter 2 Direct Current Circuits
Chapter 2 Direct Current Circuits 2.1 Introduction Nowadays, our lives are increasingly dependent upon the availability of devices that make extensive use of electric circuits. The knowledge of the electrical
More informationNumerical techniques. Chapter Difference equations
Chapter 6 Numerical techniques The differential equations (4.61), (4.62) and (4.64), in combination with boundary conditions such as equations (4.65) (4.68), constitute a two point boundary value problem.
More informationThe use of exact values at quadrature points in the boundary element method
The use of exact values at quadrature points in the boundary element method N.D. Stringfellow, R.N.L. Smith Applied Mathematics & OR Group, Cranfield University, Shrivenham, Swindon, U.K., SN6 SLA E-mail:
More informationA Krylov Subspace Accelerated Newton Algorithm
A Krylov Subspace Accelerated Newton Algorithm Michael H. Scott 1 and Gregory L. Fenves 2 ABSTRACT In the nonlinear analysis of structural systems, the Newton-Raphson (NR) method is the most common method
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 - 1 Motivation with Few Examples
Numerical Solutions of Ordinary and Partial Differential Equations Prof. G. P. Raja Shekahar Department of Mathematics Indian Institute of Technology, Kharagpur Lecture - 1 Motivation with Few Examples
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationSteps in the Finite Element Method. Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen
Steps in the Finite Element Method Chung Hua University Department of Mechanical Engineering Dr. Ching I Chen General Idea Engineers are interested in evaluating effects such as deformations, stresses,
More informationNodal and divergence-conforming boundary-element methods applied to electromagnetic scattering problems
Nodal and divergence-conforming boundary-element methods applied to electromagnetic scattering problems M. Afonso, Joao Vasconcelos, Renato Mesquita, Christian Vollaire, Laurent Nicolas To cite this version:
More informationTransactions on Modelling and Simulation vol 13, 1996 WIT Press, ISSN X
The cooling of molten glass in a mould Ph. Simons, R. Mattheij Scientific Computing Group, Eindhoven University of Technology, 56 MB Eindhoven, The Netherlands Abstract The actual formation of glass jars
More informationAxisymmetric stress analysis
5 Axisymmetric stress analysis 5.1 Introduction The problem of stress distribution in bodies of revolution (axisymmetric solids) under axisymmetric loading is of considerable practical interest. The mathematical
More informationFinite-Elements Method 2
Finite-Elements Method 2 January 29, 2014 2 From Applied Numerical Analysis Gerald-Wheatley (2004), Chapter 9. Finite-Elements Method 3 Introduction Finite-element methods (FEM) are based on some mathematical
More informationStiffness Matrices, Spring and Bar Elements
CHAPTER Stiffness Matrices, Spring and Bar Elements. INTRODUCTION The primary characteristics of a finite element are embodied in the element stiffness matrix. For a structural finite element, the stiffness
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
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 informationFormulation of the displacement-based Finite Element Method and General Convergence Results
Formulation of the displacement-based Finite Element Method and General Convergence Results z Basics of Elasticity Theory strain e: measure of relative distortions u r r' y for small displacements : x
More informationComparative study of FEM and XFEM
et International Journal on Emerging Technologies 4(2): 47-55(2013) ISSN No. (Print) : 0975-8364 ISSN No. (Online) : 2249-3255 Comparative study of FEM and XFEM Ajay Singh*, Jagdeep Kumar**, Vinod Dhull***
More informationStructural Acoustics Applications of the BEM and the FEM
Structural Acoustics Applications of the BEM and the FEM A. F. Seybert, T. W. Wu and W. L. Li Department of Mechanical Engineering, University of Kentucky Lexington, KY 40506-0046 U.S.A. SUMMARY In this
More information1 Numerical combination for nonlinear analysis of structures coupled to layered soils
1(2012) 1 23 1 Numerical combination for nonlinear analysis of structures coupled to layered soils Abstract This paper presents an alternative coupling strategy between the Boundary Element Method (BEM)
More informationY. M. TEMIS AND V. V. KARABAN 1
J. KSIAM Vol.5, No.2, 39-51, 2001 BOUNDARY ELEMENT TECHNIQUE IN TORSION PROBLEMS OF BEAMS WITH MULTIPLY CONNECTED CROSS-SECTIONS Y. M. TEMIS AND V. V. KARABAN 1 Abstract. This paper shows how boundary
More informationPreprocessor Geometry Properties )Nodes, Elements(, Material Properties Boundary Conditions(displacements, Forces )
در برنامه يك تدوين براي بعدي دو يك سازه محيط MATLAB Preprocessor Geometry Properties )Nodes, Elements(, Material Properties Boundary Conditions(displacements, Forces ) Definition of Stiffness Matrices
More informationIntroduction. Finite and Spectral Element Methods Using MATLAB. Second Edition. C. Pozrikidis. University of Massachusetts Amherst, USA
Introduction to Finite and Spectral Element Methods Using MATLAB Second Edition C. Pozrikidis University of Massachusetts Amherst, USA (g) CRC Press Taylor & Francis Group Boca Raton London New York CRC
More informationFundamentals of Linear Algebra. Marcel B. Finan Arkansas Tech University c All Rights Reserved
Fundamentals of Linear Algebra Marcel B. Finan Arkansas Tech University c All Rights Reserved 2 PREFACE Linear algebra has evolved as a branch of mathematics with wide range of applications to the natural
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationPiezoelectric Bimorph Response with Imperfect Bonding Conditions
Copyright c 28 ICCES ICCES, vol.6, no.3, pp.5-56 Piezoelectric Bimorph Response with Imperfect Bonding Conditions Milazzo A., Alaimo A. and Benedetti I. Summary The effect of the finite stiffness bonding
More informationThe method of fundamental solutions for Poisson's equation M.A. Golberg Czrc/e
The method of fundamental solutions for Poisson's equation M.A. Golberg Czrc/e ABSTRACT We show how to extend the method of fundamental solutions (MFS) to solve Poisson's equation in 2D without boundary
More informationTransactions on Modelling and Simulation vol 2, 1993 WIT Press, ISSN X
A general integral-equation formulation for anisotropic elastostatics M.M. Perez"'*, L.C. Wrobel* * Wessex Institute of Technology, University of Portsmouth, Ashurst Lodge, Ashurst, Southampton S04 2AA,
More informationBoundary Element Model for Stress Field - Electrochemical Dissolution Interactions
Boundary Element Model for Stress Field - Electrochemical Dissolution Interactions Bruce Butler Walt Disney World, Orlando, Florida bruce.butler@disney.com Manoj Chopra, Member, ASCE University of Central
More informationMCPS Algebra 2 and Precalculus Standards, Categories, and Indicators*
Content Standard 1.0 (HS) Patterns, Algebra and Functions Students will algebraically represent, model, analyze, and solve mathematical and real-world problems involving functional patterns and relationships.
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationSIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION
8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION Theodore V. Gortsas
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 16: Eigenvalue Problems; Similarity Transformations Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis I 1 / 18 Eigenvalue
More informationelastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK
Non-conforming BEM elastoplastic contact problems D. Martin and M.H. Aliabadi discretisation in Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK Abstract In this paper,
More informationInterpolation Functions for General Element Formulation
CHPTER 6 Interpolation Functions 6.1 INTRODUCTION The structural elements introduced in the previous chapters were formulated on the basis of known principles from elementary strength of materials theory.
More informationMathematics FINITE ELEMENT ANALYSIS AS COMPUTATION. What the textbooks don't teach you about finite element analysis. Chapter 3
Mathematics FINITE ELEMENT ANALYSIS AS COMPUTATION What the textbooks don't teach you about finite element analysis Chapter 3 Completeness and continuity: How to choose shape functions? Gangan Prathap
More information5.7 Cramer's Rule 1. Using Determinants to Solve Systems Assumes the system of two equations in two unknowns
5.7 Cramer's Rule 1. Using Determinants to Solve Systems Assumes the system of two equations in two unknowns (1) possesses the solution and provided that.. The numerators and denominators are recognized
More informationBack Analysis of Measured Displacements of Tunnels
Rock Mechanics and Rock Engineering 16, 173--180 (1983) Rock Mechanics and Rock Engineering 9 by Springer-Verlag 1983 Back Analysis of Measured Displacements of Tunnels By S. Sakurai and K. Takeuchi Kobe
More informationNotes on Cellwise Data Interpolation for Visualization Xavier Tricoche
Notes on Cellwise Data Interpolation for Visualization Xavier Tricoche urdue University While the data (computed or measured) used in visualization is only available in discrete form, it typically corresponds
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationChapter 9: Systems of Equations and Inequalities
Chapter 9: Systems of Equations and Inequalities 9. Systems of Equations Solve the system of equations below. By this we mean, find pair(s) of numbers (x, y) (if possible) that satisfy both equations.
More informationQuaternion Cubic Spline
Quaternion Cubic Spline James McEnnan jmcennan@mailaps.org May 28, 23 1. INTRODUCTION A quaternion spline is an interpolation which matches quaternion values at specified times such that the quaternion
More informationarxiv: v1 [math.ra] 31 Aug 2016
AN ELEMENTARY INDUCTIVE PROOF THAT AB I IMPLIES BA I FOR MATRICES J. M. ALMIRA arxiv:1608.08964v1 math.ra] 31 Aug 2016 Abstract. In this note we give an elementary demonstration of the fact that AB I n
More informationIntroduction to Finite Element computations
Non Linear Computational Mechanics Athens MP06/2012 Introduction to Finite Element computations Vincent Chiaruttini, Georges Cailletaud vincent.chiaruttini@onera.fr Outline Continuous to discrete problem
More informationModule-6: Laminated Composites-II. Learning Unit-1: M6.1. M 6.1 Structural Mechanics of Laminates
Module-6: Laminated Composites-II Learning Unit-1: M6.1 M 6.1 Structural Mechanics of Laminates Classical Lamination Theory: Laminate Stiffness Matrix To this point in the development of classical lamination
More informationEXTENDED ABSTRACT. Dynamic analysis of elastic solids by the finite element method. Vítor Hugo Amaral Carreiro
EXTENDED ABSTRACT Dynamic analysis of elastic solids by the finite element method Vítor Hugo Amaral Carreiro Supervisor: Professor Fernando Manuel Fernandes Simões June 2009 Summary The finite element
More informationOn the use of multipole methods for domain integration in the BEM
On the use of multipole methods for domain integration in the BEM A.A. Mammoii, M.S. Ingber & M.J. Brown Department of Mechanical Engineering, University of New Mexico, USA Abstract Often, a single application
More informationSyllabus (Session )
Syllabus (Session 2016-17) Department of Mathematics nstitute of Applied Sciences & Humanities AHM-1101: ENGNEERNG MATHEMATCS Course Objective: To make the students understand the concepts of Calculus,
More informationOne-dimensional Schrödinger equation
Chapter 1 One-dimensional Schrödinger equation In this chapter we will start from the harmonic oscillator to introduce a general numerical methodology to solve the one-dimensional, time-independent Schrödinger
More informationINTERPOLATION. and y i = cos x i, i = 0, 1, 2 This gives us the three points. Now find a quadratic polynomial. p(x) = a 0 + a 1 x + a 2 x 2.
INTERPOLATION Interpolation is a process of finding a formula (often a polynomial) whose graph will pass through a given set of points (x, y). As an example, consider defining and x 0 = 0, x 1 = π/4, x
More informationEigenvalues of Trusses and Beams Using the Accurate Element Method
Eigenvalues of russes and Beams Using the Accurate Element Method Maty Blumenfeld Department of Strength of Materials Universitatea Politehnica Bucharest, Romania Paul Cizmas Department of Aerospace Engineering
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 4,000 116,000 120M Open access books available International authors and editors Downloads Our
More informationShape Function Generation and Requirements
Shape Function Generation and Requirements Requirements Requirements (A) Interpolation condition. Takes a unit value at node i, and is zero at all other nodes. Requirements (B) Local support condition.
More informationALGEBRA 2. Background Knowledge/Prior Skills Knows what operation properties hold for operations with matrices
ALGEBRA 2 Numbers and Operations Standard: 1 Understands and applies concepts of numbers and operations Power 1: Understands numbers, ways of representing numbers, relationships among numbers, and number
More informationModule 10: Finite Difference Methods for Boundary Value Problems Lecture 42: Special Boundary Value Problems. The Lecture Contains:
The Lecture Contains: Finite Difference Method Boundary conditions of the second and third kind Solution of the difference scheme: Linear case Solution of the difference scheme: nonlinear case Problems
More informationToday s class. Linear Algebraic Equations LU Decomposition. Numerical Methods, Fall 2011 Lecture 8. Prof. Jinbo Bi CSE, UConn
Today s class Linear Algebraic Equations LU Decomposition 1 Linear Algebraic Equations Gaussian Elimination works well for solving linear systems of the form: AX = B What if you have to solve the linear
More informationModeling of the 3D Electrode Growth in Electroplating
Modeling of the 3D Electrode Growth in Electroplating Marius PURCAR, Calin MUNTEANU, Alexandru AVRAM, Vasile TOPA Technical University of Cluj-Napoca, Baritiu Street 26-28, 400027 Cluj-Napoca, Romania;
More informationComputational Techniques Prof. Sreenivas Jayanthi. Department of Chemical Engineering Indian institute of Technology, Madras
Computational Techniques Prof. Sreenivas Jayanthi. Department of Chemical Engineering Indian institute of Technology, Madras Module No. # 05 Lecture No. # 24 Gauss-Jordan method L U decomposition method
More informationA Reverse Technique for Lumping High Dimensional Model Representation Method
A Reverse Technique for Lumping High Dimensional Model Representation Method M. ALPER TUNGA Bahçeşehir University Department of Software Engineering Beşiktaş, 34349, İstanbul TURKEY TÜRKİYE) alper.tunga@bahcesehir.edu.tr
More informationBack Matter Index The McGraw Hill Companies, 2004
INDEX A Absolute viscosity, 294 Active zone, 468 Adjoint, 452 Admissible functions, 132 Air, 294 ALGOR, 12 Amplitude, 389, 391 Amplitude ratio, 396 ANSYS, 12 Applications fluid mechanics, 293 326. See
More informationEffect of fibre shape on transverse thermal conductivity of unidirectional composites
Sādhanā Vol. 4, Part 2, April 25, pp. 53 53. c Indian Academy of Sciences Effect of fibre shape on transverse thermal conductivity of unidirectional composites B RAGHAVA RAO,, V RAMACHANDRA RAJU 2 and
More informationElementary Linear Algebra
Matrices J MUSCAT Elementary Linear Algebra Matrices Definition Dr J Muscat 2002 A matrix is a rectangular array of numbers, arranged in rows and columns a a 2 a 3 a n a 2 a 22 a 23 a 2n A = a m a mn We
More information