ADVANCES AT TELLES METHOD APPLIED IN SCIENTIFIC VISUALIZATION
|
|
- Marshall Marshall
- 6 years ago
- Views:
Transcription
1 Copyright 009 by ABCM January 04-08, 010, Foz do Iguaçu, PR, Brazil ADVANCES AT TELLES METHOD APPLIED IN SCIENTIFIC VISUALIZATION Carlos Andrés Reyna Vera-Tudela, Universidade Federal Rural do Rio de Janeiro, Departamento de Matemática, Caixa Postal 74517, CEP , Seropédica, RJ, Brasil Marlucio Barbosa, Edivaldo Figueiredo Fontes Junior, José Claudio de Faria Telles, COPPE/Universidade Federal do Rio de Janeiro, Programa de Engenharia Civil, Caixa Postal 68506, CEP , Rio de Janeiro, RJ, Brasil Abstract. Accurate numerical integration of line integrals is of fundamental importance for any reliable implementation of the boundary element method. The Telles transform introduces an efficient means of computing singular or nearly singular integrals currently found in two-dimensional, axisymmetric and three-dimensional applications. Recently, the authors started to develop a dedicated software for scientific visualization applications. The initial procedure was writen in Java, employing the Visualization Toolkit (VTK. VTK is an open-source, freely available software for 3D computer graphics, image processing, and visualization; this consists of a C++ class library and several interpreted interface layers including Tcl/Tk, Java, and Python. However, problems can occur if abrupt changes of stress values are to be accomodated since this needs special treatment when it comes to scientific visualization and numerical accuracy. The Telles transform is a good alternative to solve the problem of numerical accuracy. Recent research suggests that the use of the transformation for internal points near the boundary and in regions with severe stress gradients can improve the look and consistency of scientific visualization representing the solution to the problem. Thus, the present work discusses some advances in the Telles transformation for efficient use in scientific visualization. Examples are included to illustrate details in the implementation and the efficiency of the method. Keywords: Boundary Element Method, Telles transform, Scientific Visualization 1. INTRODUCTION The evolution of computational tools is fundamental for the development of computer software for new applications dealing with more complex problems. Traditional algorithms must be revised and, sometimes, modified to use in high performance applications. For example, parallel processing is growing fast and practically there ere no limitations in the use of computer memory. Scientific visualization is important for science and engineering since the observation of physical behavior is fundamental for proper analysis of the problem. The authors are working in the development of software to solve, initially, D elastostatic problems with the Boundary Element Method (BEM (Brebbia, et al., The software MEMEC was developed in its original version in Fortran, but for an efficient visualization is necessary to use Object Oriented Language. Hence, the choice was Java with the Visualization Toolkit (VTK (Schroeder, et al., 004. The traditional algorithms were modified to the new requirements of computational efficiency and visualization. New algorithms had to be implemented; e.g., mesh generator and/ or interaction BEM-FEM. Initial tests indicated how efficient and friendly the software had become, relative error decreased and visualization of the problem could be accomplished. The original physical problem, comprising loads and boundary conditions could also be shown. An independent visualization of stresses and displacements had to be produced, including final deformations too. It is now possible to observe 3D visualizations and move the body in all directions. Some algorithms, however, need to be analyzed in more detail, such as the case of the Telles Transform. Thus, this paper studies the effect of the singular integrals in this context, the numerical efficiency and their influence in scientific visualization. A classical example is presented and the results are compared with a traditional BEM result, the analytical solution and Finite Element result.. THE BOUNDARY ELEMENT METHOD The Boundary Element Method (BEM (Brebbia, et al., 1984 is a numerical method for solving linear partial differential equations expressed as integral equations. In the case of elastostatic problems, without body forces, the direct formulation leads to equation of the form: C ( ( x p ( x, ξ dγ p ( x u ( x, ξ dγ = 0 ξ, x Γ ui ( ξ + u j ξ (1 Γ Γ j
2 Copyright 009 by ABCM January 04-08, 010, Foz do Iguaçu, PR, Brazil where the first integral is in Cauchy principal value sense. In this equation, Γ is the boundary of the body, ξ is the source point, x is a field point, C is a coefficient which depends on the boundary geometry at ξ, u and p are the displacement and traction components, and u and p are the fundamental displacements and tractions at x, in j direction, due to a unit load at ξ in i direction. The discretization of the body, restricted only to its boundary, is accomplished with the aid of quadratic boundary elements. Thus, Eq. (1 can be written in discrete form, after a point collocation procedure. In matrix form and amenable to prescribe the boundary conditions, the resulting system of equations is written as, Hu = Gp ( 3. THE TELLES TRANSFORMATION The boundary element method requires integration of singular kernels, present in the fundamental solution and its derivatives. Depending on the actual order and boundary dimension (i.e.; 3-D or -D, there are three distinguishable singularities: (a weak singularity (integrable in the ordinary Riemann sense, (b strong singularity (integrable in Cauchy principal value sense, (c hypersingularity (integrable in Hadamard finite part sense. If the source point is outside the element, the integrals are actually non-singular. Nevertheless, from the point of view of numerical integration, one should devote great attention also to the evaluation of nearly singular integrals, which are potentially cumbersome with the source point located close to the element being integrated. This is so because the integrands produce gradient profiles varying in strong non-uniform manner over the elements. Hence, standard integration quadrature fails to produce acceptable results. (Sladek, et al., 001 In Telles (1987 this problem is discussed in detail and efficient means to compute singular or nearly singular integrals are introduced. Initially is analyzed the use of a second-degree polynomial transformation and is verified that the use of a third-degree polynomial transformation can be used in preference to the second-degree transformation. Special attention should be given to the parametric optimization presented in this work (e.g.; Telles and Oliveira, Consider the integral = 1 1 f ( d I (3 in which f ( is singular at a point. One can always choose a third-degree relation of the form: ( (4 3 γ = a γ + bγ + cγ + d so that, in order to define a, b, c and d, the conditions below apply: (1 = 1 ( 1 = 1 d ( γ = dγ d dγ = 0 = r J (5 where r is a function of D (distance of the source point to the axis. The functional dependence of r with respect to D can be obtained to produce a minimum quadrature error. Thus, the transformation is dependent on r : j j
3 Copyright 009 by ABCM January 04-08, 010, Foz do Iguaçu, PR, Brazil J = 3a γ + bγ + c a = b = c = q = p = 3 ( 1 r / Q 3( 1 r ( r + 3γ d = b Q = 1+ 3γ γ = 3 3 [ q + ( q + p ] + q ( q + p ( 1+ r 1 ( 1+ r γ / Q / Q ( 3 r [ ] r 1+ r [ 4r ( 1 r + 3( 1 ] + 1+ r (6 For general usage, the first step is the computation of the normalized distance D as a function of distance from singular point to element. The following formula should be used: D R l R MIN (minimum MIN = (7 in which l is the real space distance between the two end points of the integration domain. Finally, the following expressions were found to determine r : r r r = ln ( D, = ln ( D, = 1, 0.05 D D D (8 4. NUMERICAL EXAMPLES As a case study, the intrinsic influence of the numerical technique adopted for scientific visualization is analyzed. The original algorithm in Fortran, was reviewed and, in many cases, modified to adapt to other languages, such as Java, that is oriented to object. Thus, special considerations had be made, principally at singular points where a high gradient of colors is found. The Telles Transformation allows for a study of the influence of this singularity in the visualization of the problem. The actual problem is a rectangular plate with a circular hole (Timoshenko and Goodier, 1970, the geometry and the material properties are shown in Fig. (1.
4 Copyright 009 by ABCM January 04-08, 010, Foz do Iguaçu, PR, Brazil Figure 1. Rectangular plate with a circular hole This problem is solved with the MEMEC software (Barbosa, et al., 008; the boundary discretization was defined with 00 linear elements and 104 internal points. It is important to note that a high number of internal points is necessary for efficient scientific visualization. However, the processing of such points is fastest due to effective processing techniques. First, a solution was obtained with a BEM implementation without advanced techniques of integration. The visualization of results for stress in x direction is shown in Fig. (; as can be observed, some of the internal points are near boundary points P1 and P (defined in Fig. (1 and the final results indicate errors when compared with the analytical result (Timoshenko and Goodier, This error leads to misinterpretations since the maximum stresses are located in incorrect positions. To solve this problem, the Telles Transformation was implemented in the BEM algorithm and the test has been repeated. The visualization of the results is presented in Fig. (3 for x direction stresses. When results are compared with the analytical stresses, one can observe a relative error of 0.33% at the P1 and P points. These points present the highest gradient of stress in x direction. Figure. Stress in x direction without the Telles Transformation procedure.
5 Copyright 009 by ABCM January 04-08, 010, Foz do Iguaçu, PR, Brazil Figure 3. Stress in x direction with the Telles Transformation procedure. Comparison of results obtained with the MEMEC software and the Finite Element Method (FEM: In order to validate the MEMEC software, results are compared with the solution obtained with a FEM implementation. The parameter calculated is the stress in x direction, Figs (3 and (4, obtained with MEMEC and FEM respectively. The discretization for the FEM implementation is shown in Fig. (5. Figures. (6 and (7 show displacements in x direction for finite elements and MEMEC respectively. The stress concentration factor (K at point P1 (Fig. 1 has been determined in order to compare with the analytical solution. Table (1 presents a summary of these results. Figure 4. Stress in x direction with FEM Figure 5. FEM grid with bi-linear quadrilateral elements Figure 6. FEM displacements in x direction Figure 7. BEM displacements in x direction with the Telles transformation.
6 Copyright 009 by ABCM January 04-08, 010, Foz do Iguaçu, PR, Brazil Table (1 shows that the error in the FEM implementation is 17.66%, whereas for the MEMEC software (BEM formulation this inaccuracy is equal to 0.33%. This result is even more evident when the actual CPU time is compared, since for the FEM implementation the computer effort has been greater than for the MEMEC software. Table 1. Stress concentration at points P1 and P Method Number Number of K Exact Solution Error (% of Nodal elements points BEM FEM CONCLUSION This paper presents the Telles transformation for scientific visualization applications. It is verified that in conjunction with efficient computational algorithms, excellent results can be obtained. The comparison with a classical finite element solution and the analytical result is seen very important when error versus discretization is observed. The MEMEC software is a tool that the authors have been developing since the early stages and further applications are also to be included in the near future. 6. ACKNOWLEDGEMENTS This work has been carried out with the support of the FAPERJ and the Fundo Setorial de Infra-Estrutura (CT- INFRA with the MCT/CNPq Programa Primeiros Projetos and Apoio às Instituições de Ensino e Pesquisa Sediadas no Estado do Rio de Janeiro; the Universidade Federal Rural do Rio de Janeiro with Grupos Emergentes project. 7. REFERENCES Barbosa, M., Fontes Junior, E.F., Vera-Tudela, C.A.R. and Telles, J.C.F., 008. O Método dos Elementos de Contorno e a Visualização Científica à Resolução de Problemas da Mecânica Computacional, Proceedings of the XXXI Congresso Nacional de Matemática Aplicada e Computacional, Belém PA, Brasil. Brebbia, C.A., Telles, J.C.F. and Wrobel, L.C., Boundary Element Techniques: Theory and Applications in Engineering, Springer, Berlin. Schroeder, W., Ken, M. and Lorensen, B., 004. The Visualization Toolkit, Kitware, Inc. Publishers. Sladek, V., Sladek, J. and Tanaka, M., 001. Numerical Integration of Logarithmic and Nearly Logarithmic Singularity in BEMs, Applied Mathematical Modeling, Vol. 5, No. 11, pp Telles, J.C.F., A Self-Adaptive Co-ordinate Transformation for Efficient Numerical Evaluation of General Boundary Element Integrals, International Journal for Numerical Methods in Engineering, Vol. 4, pp Telles, J.C.F. and Oliveira, R.F., Third Degree Polynomial Transformation for Boundary Element Integrals: Further Improvements, Engineering Analysis with Boundary Elements, Vol. 13, pp Timoshenko, S.P. and Goodier, J.N., Theory of Elasticity, Mc Graw-Hill, New York, 567 p. 8. RESPONSIBILITY NOTICE The authors are the only responsible for the printed material included in this paper.
Finite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 36
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 36 In last class, we have derived element equations for two d elasticity problems
More informationA New Infinite Boundary Element Formulation Applied to Three-Dimensional Domains
A New Infinite Boundary Element Formulation Applied to Three-Dimensional Domains 1 Dimas Betioli Ribeiro and 2 João Batista de Paiva Abstract This work presents a bi-dimensional mapped infinite boundary
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationA numerical grid and grid less (Mesh less) techniques for the solution of 2D Laplace equation
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 2014, 5(1):150-155 ISSN: 0976-8610 CODEN (USA): AASRFC A numerical grid and grid less (Mesh less) techniques for
More informationDEVELOPMENT OF A PARAMETRIC PROGRAM FOR SIMULATION OF THE STAMPING PROCESS USING FINITE ELEMENTS METHOD
DEVELOPMENT OF A PARAMETRIC PROGRAM FOR SIMULATION OF THE STAMPING PROCESS USING FINITE ELEMENTS METHOD Erick Siqueira Guidi, eguidi@uol.com.br Fernando de Azevedo Silva, fazevedo@feg.unesp.br UNESP São
More informationJournal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp
Journal of Theoretical and Applied Mechanics, Sofia, 2015, vol. 45, No. 2, pp. 59 74 SOLID MECHANICS WAVE PROPAGATION DUE TO AN EMBEDDED SEISMIC SOURCE IN A GRADED HALF-PLANE WITH RELIEF PECULIARITIES.
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 3, 2017, Lesson 1
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 3, 2017, Lesson 1 1 Politecnico di Milano, February 3, 2017, Lesson 1 2 Outline
More informationBoundary Elements 247
Boundary Elements 247 Regularizing transformation method for evaluation of singular and near-singular integrals M. Kathirkamanayagam, J.H. Curran and S. Shah Mf q/cw Eng, CWan'o, Ca/Wa Abstract Regularization
More informationINTERACTIVE GRAPHICAL APPLICATION FOR MODELING REINFORCED CONCRETE STRUCTURES
INTERACTIVE GRAPHICAL APPLICATION FOR MODELING REINFORCED CONCRETE STRUCTURES Pamela Daniela Nogueira Roque Luiz da Silva Pitangueira pameladanielanogueira@gmail.com roque@dees.ufmg.br Departamento de
More informationInvestigation of the behaviour of single span reinforced concrete historic bridges by using the finite element method
Structural Studies, Repairs and Maintenance of Heritage Architecture XI 279 Investigation of the behaviour of single span reinforced concrete historic bridges by using the finite element method S. B. Yuksel
More informationAn explicit feature control approach in structural topology optimization
th World Congress on Structural and Multidisciplinary Optimisation 07 th -2 th, June 205, Sydney Australia An explicit feature control approach in structural topology optimization Weisheng Zhang, Xu Guo
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationRELIABILITY OF THE FEM CALCULATIONS OF THE FRACTURE MECHANICS PARAMETERS
International Conference on Economic Engineering and Manufacturing Systems Braşov, 26 27 November 2009 RELIABILITY OF THE FEM CALCULATIONS OF THE FRACTURE MECHANICS PARAMETERS Galina TODOROVA, Valentin
More informationA spectral boundary element method
Boundary Elements XXVII 165 A spectral boundary element method A. Calaon, R. Adey & J. Baynham Wessex Institute of Technology, Southampton, UK Abstract The Boundary Element Method (BEM) is not local and
More informationCHAPTER 5 NUMERICAL INTEGRATION METHODS OVER N- DIMENSIONAL REGIONS USING GENERALIZED GAUSSIAN QUADRATURE
CHAPTER 5 NUMERICAL INTEGRATION METHODS OVER N- DIMENSIONAL REGIONS USING GENERALIZED GAUSSIAN QUADRATURE 5.1 Introduction Multidimensional integration appears in many mathematical models and can seldom
More informationAN OBJECT ORIENTED CLASS ORGANIZATION FOR DYNAMIC GEOMETRICALLY NON-LINEAR FEM ANALYSIS
CMNE/CILAMCE 2007 Porto, 13 a 15 de Junho, 2007 c APMTAC, Portugal 2007 AN OBJECT ORIENTED CLASS ORGANIZATION FOR DYNAMIC GEOMETRICALLY NON-LINEAR FEM ANALYSIS Flavio T. Fonseca, Roque L. S. Pitangueira
More informationRevised Sheet Metal Simulation, J.E. Akin, Rice University
Revised Sheet Metal Simulation, J.E. Akin, Rice University A SolidWorks simulation tutorial is just intended to illustrate where to find various icons that you would need in a real engineering analysis.
More informationTransactions on Modelling and Simulation vol 20, 1998 WIT Press, ISSN X
r-adaptive boundary element method Eisuke Kita, Kenichi Higuchi & Norio Kamiya Department of Mechano-Informatics and Systems, Nagoya University, Nagoya 464-01, Japan Email: kita@mech.nagoya-u.ac.jp Abstract
More informationNodal Integration Technique in Meshless Method
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 11, Issue 1 Ver. IV (Feb. 2014), PP 18-26 Nodal Integration Technique in Meshless Method Ahmed MJIDILA
More informationDevelopment of a particle difference scheme for dynamic crack propagation problem
Development of a particle difference scheme for dynamic crack propagation problem *Young-Choel Yoon 1) and Kyeong-Hwan Kim ) 1) Dept. of Civil Engineering, Myongji College, Seoul 10-776, Korea ) Dept.
More informationPROPERTIES OF NATURAL ELEMENT COORDINATES ON ANY POLYHEDRON
PROPRTIS OF NATURAL LMNT COORDINATS ON ANY POLYHDRON P. Milbradt and T. Fröbel Institute of Computer Science in Civil ngineering, Univercity of Hanover, 3067, Hanover, Germany; PH (+49) 5-76-5757; FAX
More informationEffectiveness of Element Free Galerkin Method over FEM
Effectiveness of Element Free Galerkin Method over FEM Remya C R 1, Suji P 2 1 M Tech Student, Dept. of Civil Engineering, Sri Vellappaly Natesan College of Engineering, Pallickal P O, Mavelikara, Kerala,
More informationAdaptivity For Meshfree Point Collocation Methods In Linear Elastic Solid Mechanics
University of South Carolina Scholar Commons Theses and Dissertations 2015 Adaptivity For Meshfree Point Collocation Methods In Linear Elastic Solid Mechanics Joshua Wayne Derrick University of South Carolina
More informationMODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
A. N. Johnson et al., Int. J. Comp. Meth. and Exp. Meas., Vol. 3, No. 3 (2015) 269 278 MODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
More informationStudy of Convergence of Results in Finite Element Analysis of a Plane Stress Bracket
RESEARCH ARTICLE OPEN ACCESS Study of Convergence of Results in Finite Element Analysis of a Plane Stress Bracket Gowtham K L*, Shivashankar R. Srivatsa** *(Department of Mechanical Engineering, B. M.
More informationCOMPUTER AIDED ENGINEERING. Part-1
COMPUTER AIDED ENGINEERING Course no. 7962 Finite Element Modelling and Simulation Finite Element Modelling and Simulation Part-1 Modeling & Simulation System A system exists and operates in time and space.
More informationUsing three-dimensional CURVIC contact models to predict stress concentration effects in an axisymmetric model
Boundary Elements XXVII 245 Using three-dimensional CURVIC contact models to predict stress concentration effects in an axisymmetric model J. J. Rencis & S. R. Pisani Department of Mechanical Engineering,
More information3D Nearest-Nodes Finite Element Method for Solid Continuum Analysis
Adv. Theor. Appl. Mech., Vol. 1, 2008, no. 3, 131-139 3D Nearest-Nodes Finite Element Method for Solid Continuum Analysis Yunhua Luo Department of Mechanical & Manufacturing Engineering, University of
More informationA POINTWISE APPROACH FOR ENFORCEMENT OF ESSENTIAL BOUNDARY CONDITIONS IN THE ISOGEOMETRIC ANALYSIS * Y. BAZARGAN-LARI **
IJST, Transactions of Mechanical Engineering, Vol. 38, No. M1 +, pp 167-179 Printed in The Islamic Republic of Iran, 2014 Shiraz University A POINTWISE APPROACH FOR ENFORCEMENT OF ESSENTIAL BOUNDARY CONDITIONS
More informationIsogeometric Analysis Application to Car Crash Simulation
Isogeometric Analysis Application to Car Crash Simulation S. Bouabdallah 2, C. Adam 2,3, M. Zarroug 2, H. Maitournam 3 De Vinci Engineering Lab, École Supérieure d Ingénieurs Léonard de Vinci 2 Direction
More informationSmooth finite elements
Smooth finite elements seamless handling of incompressibility, distorted and polygonal meshes; links with equilibrium methods Stéphane Bordas * Nguyen-Xuan Hung ** Nguyen-Dang Hung *** * University of
More informationFigure 1: execution of the program with the definition of the parameters.
The Fortran90 code BEMdyn-incompressible is a general boundary element software suitable to solve problems of incompressible incremental non-linear elasticity with time-harmonic external loads of prescribed
More informationD DAVID PUBLISHING. 1. Introduction
Journal of Civil Engineering and Architecture 9 (2015) 701-714 doi: 10.17265/1934-7359/2015.06.008 D DAVID PUBLISHING Parametric Modeling of Static Three-Dimensional Problem of Buried Pipeline through
More informationModule 1: Introduction to Finite Element Analysis. Lecture 4: Steps in Finite Element Analysis
25 Module 1: Introduction to Finite Element Analysis Lecture 4: Steps in Finite Element Analysis 1.4.1 Loading Conditions There are multiple loading conditions which may be applied to a system. The load
More informationISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 2, Issue 3, September 2012
Mitigation Curves for Determination of Relief Holes to Mitigate Concentration Factor in Thin Plates Loaded Axially for Different Discontinuities Shubhrata Nagpal, S.Sanyal, Nitin Jain Abstract In many
More informationROTATIONAL DEPENDENCE OF THE SUPERCONVERGENT PATCH RECOVERY AND ITS REMEDY FOR 4-NODE ISOPARAMETRIC QUADRILATERAL ELEMENTS
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng, 15, 493±499 (1999) ROTATIONAL DEPENDENCE OF THE SUPERCONVERGENT PATCH RECOVERY AND ITS REMEDY FOR 4-NODE ISOPARAMETRIC QUADRILATERAL
More informationTopology Optimization of Two Linear Elastic Bodies in Unilateral Contact
2 nd International Conference on Engineering Optimization September 6-9, 2010, Lisbon, Portugal Topology Optimization of Two Linear Elastic Bodies in Unilateral Contact Niclas Strömberg Department of Mechanical
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 24
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 24 So in today s class, we will look at quadrilateral elements; and we will
More informationThe Dynamic Response of an Euler-Bernoulli Beam on an Elastic Foundation by Finite Element Analysis using the Exact Stiffness Matrix
Journal of Physics: Conference Series The Dynamic Response of an Euler-Bernoulli Beam on an Elastic Foundation by Finite Element Analysis using the Exact Stiffness Matrix To cite this article: Jeong Soo
More informationSet No. 1 IV B.Tech. I Semester Regular Examinations, November 2010 FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 Answer any FIVE Questions All Questions carry equal marks
More informationThe DRM-MD integral equation method for the numerical solution of convection-diffusion
The DRM-MD integral equation method for the numerical solution of convection-diffusion equation V. Popov & H. Power Wessex Institute of Technology, Ashurst Lodge, Ashurst, Abstract This work presents a
More informationA meshfree weak-strong form method
A meshfree weak-strong form method G. R. & Y. T. GU' 'centre for Advanced Computations in Engineering Science (ACES) Dept. of Mechanical Engineering, National University of Singapore 2~~~ Fellow, Singapore-MIT
More informationModule 2: Single Step Methods Lecture 4: The Euler Method. The Lecture Contains: The Euler Method. Euler's Method (Analytical Interpretations)
The Lecture Contains: The Euler Method Euler's Method (Analytical Interpretations) An Analytical Example file:///g /Numerical_solutions/lecture4/4_1.htm[8/26/2011 11:14:40 AM] We shall now describe methods
More informationKeywords: Radiation Patterns, Antennas, Scientific Visualization, VTK.
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm COMPUTATION AND VISUALIZATION OF THREE DIMENSIONAL RADIATION PATTERNS OF ANTENNAS M. Joaquim
More informationTheory of Elasticity, Article 18, Generic 2D Stress Distributions
Theory of Elasticity, Article 18, Generic 2D Stress Distributions One approach to solve 2D continuum problems analytically is to use stress functions. Article 18 of the third edition of Theory of Elasticity
More information1. Introduction. 2. Velocity Distribution
Comparison of Direct And Indirect Boundary Element Methods For The Calculation of Potential Flow Around An Elliptic Cylinder With Linear Element Approach Muhammad Mushtaq*, Nawazish Ali Shah,Ghulam Muhammad,
More informationThe solution of boundary element problems using symbolic computation A. Almeida, H.L. Pina Institute Superior Tecnico, Lisbon, Portugal
The solution of boundary element problems using symbolic computation A. Almeida, H.L. Pina Institute Superior Tecnico, Lisbon, Portugal Abstract The paper explores the potential of using Symbolic Computation
More informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More informationModelling of Isogeometric Analysis for Plane Stress Problem Using Matlab GUI Norliyana Farzana Binti Zulkefli, Ahmad Razin Bin Zainal Abidin
Modelling of Isogeometric Analysis for Plane Stress Problem Using Matlab GUI Norliyana Farzana Binti Zulkefli, Ahmad Razin Bin Zainal Abidin Faculty of Civil Engineering, Universiti Teknologi Malaysia,
More informationSCREW-BASED RELATIVE JACOBIAN FOR MANIPULATORS COOPERATING IN A TASK
ABCM Symposium Series in Mechatronics - Vol. 3 - pp.276-285 Copyright c 2008 by ABCM SCREW-BASED RELATIVE JACOBIAN FOR MANIPULATORS COOPERATING IN A TASK Luiz Ribeiro, ribeiro@ime.eb.br Raul Guenther,
More informationStress Analysis of thick wall bellows using Finite Element Method
Stress Analysis of thick wall bellows using Finite Element Method Digambar J. Pachpande Post Graduate Student Department of Mechanical Engineering V.J.T.I. Mumbai, India Prof. G. U. Tembhare Assistant
More informationOn cliques of Helly Circular-arc Graphs
Electronic Notes in Discrete Mathematics 30 (2008) 117 122 www.elsevier.com/locate/endm On cliques of Helly Circular-arc Graphs Min Chih Lin a,1, Ross M. McConnell b, Francisco J. Soulignac a,2 and Jayme
More informationSIMULATION AND ANALYSIS OF TOUCH TRIGGER STYLUS PERFORMANCE
Page 905 SIMULATION AND ANALYSIS OF TOUCH TRIGGER STYLUS PERFORMANCE Verílton Nunes da Silva, verilton@yahoo.com.br Luiz Henrique Melo Silva Nóbrega, rikinobrega@hotmail.com Francisco Augusto Vieira da
More informationTHE USE OF OPENCAVOK FINITE ELEMENT SOFTWARE IN UNDERGRADUATED COURSES OF MECHANICAL ENGINEERING
THE USE OF OPENCAVOK FINITE ELEMENT SOFTWARE IN UNDERGRADUATED COURSES OF MECHANICAL ENGINEERING Pierre Lamary Roberto de Araújo Bezerra Ranon Bezerra Department of Mechanical Engineering, Federal University
More informationAn improved DRM representation of partial derivatives
An improved DRM representation of partial derivatives B. Natalini',V.Popov2 'Facultad deingenieria, Universidad Nacional del Nordeste, Argentina Wessex Institute of Technology, United Kingdom Abstract
More informationCE Advanced Structural Analysis. Lab 4 SAP2000 Plane Elasticity
Department of Civil & Geological Engineering COLLEGE OF ENGINEERING CE 463.3 Advanced Structural Analysis Lab 4 SAP2000 Plane Elasticity February 27 th, 2013 T.A: Ouafi Saha Professor: M. Boulfiza 1. Rectangular
More informationAnalysis of Distortion Parameters of Eight node Serendipity Element on the Elements Performance
Analysis of Distortion Parameters of Eight node Serendipity Element on the Elements Performance Vishal Jagota & A. P. S. Sethi Department of Mechanical Engineering, Shoolini University, Solan (HP), India
More informationENGINEERING TRIPOS PART IIA FINITE ELEMENT METHOD
ENGINEERING TRIPOS PART IIA LOCATION: DPO EXPERIMENT 3D7 FINITE ELEMENT METHOD Those who have performed the 3C7 experiment should bring the write-up along to this laboratory Objectives Show that the accuracy
More informationIntroduction to Design Optimization
Introduction to Design Optimization First Edition Krishnan Suresh i Dedicated to my family. They mean the world to me. ii Origins of this Text Preface Like many other textbooks, this text has evolved from
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationElectromagnetic Field Numerical Modeling using BEM2D
Volume 56, Number 5, 2015 197 Electromagnetic Field Numerical Modeling using BEM2D Adina Răcășan, Claudia Păcurar, Călin Munteanu, Vasile Țopa, Claudia Constantinescu, Lorand Szabo and Marius Dodea Faculty
More informationObject-oriented programming in boundary element methods using C11
Advances in Engineering Software 30 (1999) 127±132 Object-oriented programming in boundary element methods using C11 Wenqing Wang*, Xing Ji, Yuangong Wang Department of Engineering Mechanics and Technology,
More informationAn automatic simulation tool for thermal analysis of gravity dams by BFM
Boundary Elements and Other Mesh Reduction Methods XXXVI 491 An automatic simulation tool for thermal analysis of gravity dams by BFM Jianming Zhang State Key Laboratory of Advanced Design and Manufacturing
More informationOn the fast-multipole implementation of the simplified hybrid boundary element method
On the fast-multipole implementation of the simplified hybrid boundary element method Hélvio de Farias Costa Peixoto Ney Augusto Dumont hfcpeixoto@gmail.com dumont@puc-rio.com.br Pontifícia Universidade
More information: What is Finite Element Analysis (FEA)?
Q: What is Finite Element Analysis (FEA)? A1: It is a numerical technique for finding approximate solutions of partial differential equations (PDE) as well as of integral equations. The solution approach
More informationME 345: Modeling & Simulation. Introduction to Finite Element Method
ME 345: Modeling & Simulation Introduction to Finite Element Method Examples Aircraft 2D plate Crashworthiness 2 Human Heart Gears Structure Human Spine 3 F.T. Fisher, PhD Dissertation, 2002 Fluid Flow
More informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II A Brief Overview of Finite Element Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 1 / 25 Overview Basic concepts Mathematical
More informationStatic analysis of eolicblade through finite element method and OOP C++
International Conference on Control, Engineering & Information echnology (CEI 4) Proceedings Copyright IPCO-4 ISSN 56-568 Static analysis of eolicblade through finite element method and OOP C++ MateusDantas,
More informationFinite Element Modeling and Failure Analysis of Roll Bending. Forming of GLARE Laminates
Finite Element Modeling and Failure Analysis of Roll Bending Forming of GLARE Laminates Jingming Tian, Gang Tao, Cheng Liu, Huaguan Li, Xian Zhang, Jie Tao* College of Materials Science and Technology,
More informationTHE IDENTIFICATION OF THE BOUNDARY GEOMETRY WITH CORNER POINTS IN INVERSE TWO-DIMENSIONAL POTENTIAL PROBLEMS
TASK QUARTERLY 6 No 4(2002), 651 660 THE IDENTIFICATION OF THE BOUNDARY GEOMETRY WITH CORNER POINTS IN INVERSE TWO-DIMENSIONAL POTENTIAL PROBLEMS EUGENIUSZ ZIENIUK AND WIESŁAW GABREL Institute of Computer
More informationAbstract. 1 Introduction
Influence of the location of the contact point in node-to-point contact approaches using non conforming boundary element discretizations A. Blazquez & F. Paris Department of continuum Mechanics, E.S. de
More informationNumerical simulation of multilayer grounding grids in a user-friendly open-source CAD interface
Numerical simulation of multilayer grounding grids in a user-friendly open-source CAD interface José París, Ignasi Colominas, Xesús Nogueira, Fermín Navarrina, Manuel Casteleiro Group of Numerical Methods
More informationEXACT BUCKLING SOLUTION OF COMPOSITE WEB/FLANGE ASSEMBLY
EXACT BUCKLING SOLUTION OF COMPOSITE WEB/FLANGE ASSEMBLY J. Sauvé 1*, M. Dubé 1, F. Dervault 2, G. Corriveau 2 1 Ecole de technologie superieure, Montreal, Canada 2 Airframe stress, Advanced Structures,
More informationMATLAB. Advanced Mathematics and Mechanics Applications Using. Third Edition. David Halpern University of Alabama CHAPMAN & HALL/CRC
Advanced Mathematics and Mechanics Applications Using MATLAB Third Edition Howard B. Wilson University of Alabama Louis H. Turcotte Rose-Hulman Institute of Technology David Halpern University of Alabama
More informationSecond-order shape optimization of a steel bridge
Computer Aided Optimum Design of Structures 67 Second-order shape optimization of a steel bridge A.F.M. Azevedo, A. Adao da Fonseca Faculty of Engineering, University of Porto, Porto, Portugal Email: alvaro@fe.up.pt,
More informationSETTLEMENT OF A CIRCULAR FOOTING ON SAND
1 SETTLEMENT OF A CIRCULAR FOOTING ON SAND In this chapter a first application is considered, namely the settlement of a circular foundation footing on sand. This is the first step in becoming familiar
More informationChapter 1 Introduction
Chapter 1 Introduction GTU Paper Analysis (New Syllabus) Sr. No. Questions 26/10/16 11/05/16 09/05/16 08/12/15 Theory 1. What is graphic standard? Explain different CAD standards. 2. Write Bresenham s
More informationSimulation of rotation and scaling algorithm for numerically modelled structures
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Simulation of rotation and scaling algorithm for numerically modelled structures To cite this article: S K Ruhit et al 2018 IOP
More informationGuidelines for proper use of Plate elements
Guidelines for proper use of Plate elements In structural analysis using finite element method, the analysis model is created by dividing the entire structure into finite elements. This procedure is known
More information1 Initial assumptions
Methods of evaluation of global and local ship hull deflections based on periodic geodesic measurements assuming linear model of deflections J. Niebylski Ship Research Institute, Technical University ofszczecin,
More informationFinite element algorithm with adaptive quadtree-octree mesh refinement
ANZIAM J. 46 (E) ppc15 C28, 2005 C15 Finite element algorithm with adaptive quadtree-octree mesh refinement G. P. Nikishkov (Received 18 October 2004; revised 24 January 2005) Abstract Certain difficulties
More informationIsogeometric Collocation Method
Chair for Computational Analysis of Technical Systems Faculty of Mechanical Engineering, RWTH Aachen University Isogeometric Collocation Method Seminararbeit By Marko Blatzheim Supervisors: Dr. Stefanie
More informationRELIABILITY ANALYSIS IN STRUCTURES SUBMITTED TO FATIGUE USING THE BOUNDARY ELEMENT METHOD
RELIABILITY ANALYSIS IN STRUCTURES SUBMITTE TO FATIGUE USING THE BOUNARY ELEMENT METHO Edson enner Leonel, edleonel@sc.usp.br André Teófilo Beck, atbeck@sc.usp.br Wilson Sérgio Venturini, venturin@sc.usp.br
More informationNUMERICAL ANALYSIS OF ENGINEERING STRUCTURES (LINEAR ELASTICITY AND THE FINITE ELEMENT METHOD)
NUMERICAL ANALYSIS OF ENGINEERING STRUCTURES (LINEAR ELASTICITY AND THE FINITE ELEMENT METHOD) NUMERICAL ANALYSIS OF ENGINEERING STRUCTURES (LINEAR ELASTICITY AND THE FINITE ELEMENT METHOD) Author: Tamás
More informationDetermination of Free Surface in Steady-State Seepage through a Dam with Toe Drain by the Boundary Element Method
Determination of Free Surface in Steady-State Seepage through a Dam with Toe Drain by the Boundary Element Method Somchart Chantasiriwan Faculty of Engineering, Thammasat University, Rangsit Campus, Pathum
More informationParametric study of An Innovative Approximate Method of Continuous Beam
Parametric study of An Innovative Approximate Method of Continuous Beam Dipak J.Varia Research Scholar Rai University Ahmedabad, India djvaria@gmail.com Dr.Harshvadan S. Patel Applied Mechanics Department
More informationOn the Matlab S-FEM code for 2D problems using T3 and Q4 elements
On the Matlab S-FEM code for 2D problems using T3 and Q4 elements These codes were developed by Liu, Nguyen and workers. The detailed theoretical background, formulation and implementation procedure are
More informationModeling Ground Water Problems Using the Complex Polynomial Method
Modeling Ground Water Problems Using the Complex Polynomial Method A. W. Bohannon and T. V. Hromadka, AS-0020 Abstract Numerical methods for solving the governing partial differential equations involved
More informationIJMH - International Journal of Management and Humanities ISSN:
EXPERIMENTAL STRESS ANALYSIS SPUR GEAR USING ANSYS SOFTWARE T.VADIVELU 1 (Department of Mechanical Engineering, JNTU KAKINADA, Kodad, India, vadimay28@gmail.com) Abstract Spur Gear is one of the most important
More informationLacunary Interpolation Using Quartic B-Spline
General Letters in Mathematic, Vol. 2, No. 3, June 2017, pp. 129-137 e-issn 2519-9277, p-issn 2519-9269 Available online at http:\\ www.refaad.com Lacunary Interpolation Using Quartic B-Spline 1 Karwan
More informationFRANC3D / OSM Tutorial Slides
FRANC3D / OSM Tutorial Slides October, 2003 Cornell Fracture Group Tutorial Example Hands-on-training learn how to use OSM to build simple models learn how to do uncracked stress analysis using FRANC3D
More informationADAPTIVE APPROACH IN NONLINEAR CURVE DESIGN PROBLEM. Simo Virtanen Rakenteiden Mekaniikka, Vol. 30 Nro 1, 1997, s
ADAPTIVE APPROACH IN NONLINEAR CURVE DESIGN PROBLEM Simo Virtanen Rakenteiden Mekaniikka, Vol. 30 Nro 1, 1997, s. 14-24 ABSTRACT In recent years considerable interest has been shown in the development
More informationDie Wear Profile Investigation in Hot Forging
Die Wear Profile Investigation in Hot Forging F. R. Biglari, M Zamani Abstract In this study, the wear profile on the die surface during the hot forging operation for an axisymmetric cross-section is examined.
More informationUsing Semi-Regular 4 8 Meshes for Subdivision Surfaces
Using Semi-Regular 8 Meshes for Subdivision Surfaces Luiz Velho IMPA Instituto de Matemática Pura e Aplicada Abstract. Semi-regular 8 meshes are refinable triangulated quadrangulations. They provide a
More informationA mesh refinement technique for the boundary element method based on local error analysis
A mesh refinement technique for the boundary element method based on local error analysis J. J. Rodriguez
More informationSolving a Two Dimensional Unsteady-State. Flow Problem by Meshless Method
Applied Mathematical Sciences, Vol. 7, 203, no. 49, 242-2428 HIKARI Ltd, www.m-hikari.com Solving a Two Dimensional Unsteady-State Flow Problem by Meshless Method A. Koomsubsiri * and D. Sukawat Department
More informationC. A. D. Fraga Filho 1,2, D. F. Pezzin 1 & J. T. A. Chacaltana 1. Abstract
Advanced Computational Methods and Experiments in Heat Transfer XIII 15 A numerical study of heat diffusion using the Lagrangian particle SPH method and the Eulerian Finite-Volume method: analysis of convergence,
More informationA Multiple Constraint Approach for Finite Element Analysis of Moment Frames with Radius-cut RBS Connections
A Multiple Constraint Approach for Finite Element Analysis of Moment Frames with Radius-cut RBS Connections Dawit Hailu +, Adil Zekaria ++, Samuel Kinde +++ ABSTRACT After the 1994 Northridge earthquake
More informationACCURACY OF NUMERICAL SOLUTION OF HEAT DIFFUSION EQUATION
Scientific Research of the Institute of Mathematics and Computer Science ACCURACY OF NUMERICAL SOLUTION OF HEAT DIFFUSION EQUATION Ewa Węgrzyn-Skrzypczak, Tomasz Skrzypczak Institute of Mathematics, Czestochowa
More information