Numerical modelling of seismic waves using imageprocessing
|
|
- Ambrose Gilbert Blake
- 6 years ago
- Views:
Transcription
1 Southern Cross University 23rd Australasian Conference on the Mechanics of Structures and Materials 2014 Numerical modelling of seismic waves using imageprocessing and quadtree meshes C Birk University of New South Wales C Song University of New South Wales Publication details Birk, C, Song, C 2014, 'Numerical modelling of seismic waves using image-processing and quadtree meshes', in ST Smith (ed.), 23rd Australasian Conference on the Mechanics of Structures and Materials (ACMSM23), vol. II, Byron Bay, NSW, 9-12 December, Southern Cross University, Lismore, NSW, pp ISBN: epublications@scu is an electronic repository administered by Southern Cross University Library. Its goal is to capture and preserve the intellectual output of Southern Cross University authors and researchers, and to increase visibility and impact through open access to researchers around the world. For further information please contact epubs@scu.edu.au.
2 23rd Australasian Conference on the Mechanics of Structures and Materials (ACMSM23) Byron Bay, Australia, 9-12 December 2014, S.T. Smith (Ed.) NUMERICAL MODELLING OF SEISMIC WAVES USING IMAGE- PROCESSING AND QUADTREE MESHES C. Birk* School of Civil and Environmental Engineering, University of New South Wales Sydney, NSW, 2052, Australia. (Corresponding Author) C. Song School of Civil and Environmental Engineering, University of New South Wales Sydney, NSW, 2052, Australia. ABSTRACT This paper presents an automatic meshing procedure which is based on image processing techniques. It is suitable for problems with heterogeneous geometric features, which cannot be described by accurate drawings, in particular if details of very different scales are present. The scaled boundary finite element method facilitates structural analysis on the resulting quadtree mesh, as it incorporates hanging nodes in an elegant and straightforward way. The proposed method is applied to the problem of seismic waves in highly heterogeneous soils. Here, the maximum mesh size is governed by the material properties and frequency content of the excitation. A large-scale two-dimensional soilstructure interaction problem is studied to illustrate the potential of the proposed approach. KEYWORDS Seismic waves, image processing, quadtree mesh, scaled boundary finite element method. INTRODUCTION Earthquakes are natural events with potentially catastrophic consequences. For the prediction of ground motions and the resulting effects on our built environment, it is crucial to understand how seismic waves interact with geological features, such as soil layers, cracks and faults. With the rapid development of computing capacity and numerical methods, computer simulations of seismic waves are becoming increasingly important. The numerical modelling of wave fields in the ground is challenging for a number of reasons, such as the requirement to accurately model radiation damping, potential nonlinearities and the need of efficient modelling techniques to handle the large size of the problem. Recent developments in computational seismology are mainly based on finite difference and finite element methods (Aochi et al. 2013; Semblat and Brioist 2000), or on spectral element approaches (Chaljub et al. 2007; Mazzieri et al. 2011). A comprehensive review of numerical methods is given in Semblat (2011). One common feature of these methods is that they require a discretization of the region to be analysed. Realistic soil deposits, however, are highly heterogeneous, involving several materials and complex geometries. Geomechanical analyses are typically based on approximate graphical information rather than on precise engineering drawings. Moreover, seismic analysis may require the resolution of geometrical features at very different scales. The generation of high-quality meshes based on this information is a time-consuming and challenging task which requires considerable experience. This is This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this license, visit 819
3 further complicated by the existence of additional mesh requirements for wave propagation analyses depending on material properties and frequency content. To overcome these challenges, an automatic meshing approach which requires only a graphical representation of the soil region is proposed in this paper. It is based on producing a quadtree decomposition of the domain using image processing techniques. Quadtree structures are highly efficient in capturing geometrical features at very different scales and in transitioning from coarser to finer mesh resolutions (Greaves and Borthwick 1999). Finite and spectral element approaches, however, cannot be used on quadtree meshes straightforwardly due to the presence of hanging nodes (Tabarraei and Sukumar 2005). The scaled boundary finite element method (SBFEM) (Song 2009) provides an elegant way of combining the automatic quadtree mesh generation with efficient dynamic analysis. Here, each quadrilateral element is represented by a rectangular scaled boundary finite element subdomain, which is modeled using a number of line elements. Arbitrary node configurations can be created, which allows hanging nodes to be easily incorporated. If a 2:1 rule is applied, each SBFE subdomain is characterized by one of only six different stiffness and mass matrices. Making use of this property, a highly efficient computational model for the time-domain anaysis of large-scale wave propagation problems is obtained. AUTOMATIC MESH GENERATION Consider a graphic representation of an example soil deposit as shown in Figure 1(a). Here, the black and grey coloured regions represent rock and soft soil with given material properties,, and,,, respectively. The white regions indicate discontinuities in the form of cracks and a cavity. (a) Image (b) Quadtree mesh Figure 1. Soil deposit with discontinuities The following algorithm is based on representing any rectangular digital image as a matrix [ ] of size, where and are the number of pixels in the vertical and horizontal direction. For a greyscale image, the value of each component of [ ] represents the colour of the corresponding pixel. This mathematical representation of the image can be used to obtain a quadtree decomposition of the region. Since this is based on recursively splitting a quadratic cell into smaller squares, [ ] is expanded to a square matrix by filling it with values corresponding to the void colour. The Matlab function qtdecomp is used to decompose the image [ ] Here, a quadratic cell is further divided if the maximum and minimum value inside a block differ by more than a specified threshold, i.e. if the colour changes. The function returns blocks that are not smaller than a specified minimum and not bigger than a specified maximum size. While the minimum size determines how smoothly geometrical features can be represented and is limited to 1 pixel, the maximum size governs the accuracy of the subsequent structural analysis. In a wave propagation analysis, the maximum allowable element size depends on the material properties, frequency content and element type. Based on the requirement of at least 10 nodes per wavelength to avoid dispersion and dissipation errors (Ham and Bathe 2012), different maximum element sizes are imposed on soil regions with different material properties. Moreover, a 2:1 rule is applied to the quadtree mesh, allowing neighbouring elements to differ in size by not more than a factor of 2. The resulting quadtree decomposition of the example soil deposit is shown in Figure 1(b). Observe that, ACMSM
4 - the overall mesh size in the top soil is smaller than in the rock according to its smaller shear wave velocity, - geometrical discontinuities (cracks, interfaces, cavities) are captured by very small elements - the element size rapidly transitions to the maximum size away from such discontinuities. STRUCTURAL ANALYSIS USING SBFEM As explained earlier, standard finite element formulations are not compatible on quadtree meshes. The scaled boundary finite element method provides an elegant link between automatic mesh generation and structural analysis. Here, each quadrilateral element is regarded as a special case of a general - sided polygon. If the 2:1 rule is applied, only the 6 configurations shown in Figure 2 are possible. Figure 2. Possible polygon configurations in a quadtree mesh abiding the 2:1 rule Each polygon is treated as a scaled boundary finite element subdomain. In the SBFEM, a local coordinate system is introduced as shown in Figure 3 for a 5-sided polygon as an example. Figure 3. Scaled boundary finite element coordinates Only the boundary is discretized using line elements. A radial coordinate is defined pointing outwards from the scaling centre, which is selected as the centroid of each quadrilateral domain. The circumferential coordinate is defined counter-clockwise along the boundary. The local Cartesian coordinates, are related to the so-called scaled boundary coordinates and by the scaled boundary transformation given in Eq. (1), [ ], [ ] (1) Here, the nodal coordinates and are interpolated by shape functions [ ] in terms of the circumferential coordinate. Formulating the governing equations of elastodynamics using the scaled boundary transformation and applying the method of weigthed residuals in the circumferential direction, the SBFE equation in dynamic stiffness (2) is obtained, [ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ] [ ] (2) In Eq. (2), [ ], [ ], [ ] and [ ] are coefficient matrices that depend on the geometry of the boundary and the constitutive law only. These matrices can be assembled using FE technology. For 2- node line elements and linear elasticity closed-form expressions can be found in (Wolf and Song, 1996). For, Eq. 2 yields an algebraic Riccati equation for the static stiffness matrix [ ], which can be solved using Schur decomposition (Song 2009). The mass matrix [ ] of a polygonal subdomain follows from a low-frequency expansion of Eq. (2) (Song 2009). The stiffness and mass matrices of all polygons are assembled in global equations of motion, [ ] [ ] (3) ACMSM
5 which can be solved using standard time integration algorithms, such as Newmark s method. Note that in a two-dimensional analysis the stiffness matrix of a quadtree element only depends on the polygon configuration and material data, but not on the element size. Therefore, only a few element stiffness and mass matrices have to be calculated. The mass matrices are scaled according to their size when assembling contributions of individual polygons into the global mass matrix NUMERICAL EXAMPLES The regular three-layered domain shown in Figure 4 is analysed to verify the proposed method. A deep bedrock is overlain by two softer layers of thickness each.,,,,,, Figure 4. Three-layered soil deposit: material data and mesh The quadtree mesh obtained automatically from an image of resolution is also shown. Here, one pixel corresponds to a physical dimension of. It can be seen that smaller elements are automatically obtained in the softer two top layers, corresponding to their lower shear wave velocity. The mesh has been designed to transmit frequencies up to. Elements of the smallest possible size of 1 pixel are obtained in the transition regions between soil layers, since the layer thickness is not a power of 2. The first 30 natural frequencies of the soil deposit under plane stress, fixed at the bottom are calculated and compared to reference solutions obtained using the commercial software ANSYS. Excellent agreement is observed in Figure 5. Figure 5. Natural frequencies of a three-layered soil deposit In a second example, a large-scale soil-structure interaction problem is addressed to illustrate the potential of the proposed method. A three-layered soil deposit of total thickness is considered, which extends 8192 in the horizontal direction. The individual layers are, and thick, from bottom to top, respectively. The same soil properties as in Figure 4 are assumed. The soil-deposit is represented by an image of size pixel. Thus, the smallest possible element size of 1 pixel corresponds to a physical dimension of. A structural detail consisting of two concrete high-rise buildings and a concrete tunnel with,, is added at the centre of the image, as shown in Figure 6. A quadtree mesh is created automatically using the image shown in Figure 6(a). It is designed to transmit frequencies up to and comprises elements and nodes. The maximum mesh size in the soil layers 1-3 is 64, 64 and 32 pixel, corresponding to, and, respectively. ACMSM
6 The maximum mesh size in the concrete is to material boundaries and geometrical features., however, much smaller elements are obtained close (a) (b) (c) Figure 6. Three-layered soil deposit with structural detail. (a) full image, (b) detail, (c) mesh detail Waves in the soil are caused by a Ricker pulse point load,, which is applied at point S,,, (4) The resulting displacements at points A, B and C are shown in Figure 7. All coordinates are measured from the bottom left corner. Figure 7. Displacement response due to Ricker pulse acting below the surface Point S is approximately to the left of the structures and below the surface. Points A and B are located at the soil surface, approximately to the left and between the two high-rise buildings, respectively. Points A and C are located directly above the source, with C being in the middle soil layer. Figure 7 indicates the vertical displacements are amplified when travelling upwards through increasingly softer material. It can also be seen that the maximum displacement recorded at point B, i.e. between the two buildings, is approximately 50% of the maximum displacement obtained at point A. Snapshots of the displacement field at different times are shown in Figure 8. These further illustrate the propagation of the seismic wave, reflection at the rigid bottom and refraction at interfaces. ACMSM
7 (a) (b) (c) (d) Figure 8. Horizontal displacement due to vertical Ricker pulse at point S The CPU times recorded for the mesh generation and the time-domain analysis (1000 steps) are 91.3 s and s, respectively, on a desktop computer with Intel Core i CPU@3.30 GHz and 8GB RAM. CONCLUSIONS This paper has presented a versatile strategy for automatic mesh generation and dynamic analysis. Based on square scaled boundary polygons, the ubiquitous hanging node problem of classical FE approaches is overcome. The stability of the proposed scheme is not limited by the size of the problem. Using techniques based on image-processing, quadtree meshes can be created automatically in a highly efficient manner. This has been demonstrated for a high-resolution image of size pixel. The structural analysis is based on the SBFEM and standard time-stepping techniques. First promising results for the analysis of seismic wave propagation in two-dimensional soil have been presented. The extension to 3D is the subject of current research. REFERENCES Aochi, H., Ulrich, T., Ducellier, A., Dupros, F., and Michea, D. (2013) Finite difference simulations of seismic wave propagation for understanding earthquake physics and predicting ground motions: Advances and challenges, Journal of Physics: Conference Series, Vol. 454, pp Chaljub, E., Komatitsch, D., Vilotte, J.-P., Capdeville, Y., Valette, B. and Festa, G. (2007) Spectralelement analysis in seismology, Advances in Geophysics, Vol. 48, pp Greaves, D.M. and Borthwick, A. (1999) Hierarchical tree-based finite element mesh generation, International Journal for Numerical Methods in Engineering, Vol. 45, pp Ham, S. and Bathe, K.J. (2012) A finite element method enriched for wave propagation problems, Computers and Structures, Vol , pp Mazzieri, I., Smerzini, C., Antonietti, P.F., Rapetti, F., Stupazzini, M., Paolucci, R. and Quarteroni, A. (2011) Non-conforming spectral approximations for the elastic wave equation in heterogeneous media, Proceedings, ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering, COMPDYN 2011, Corfu, Greece, May Semblat, J.F. and Brioist, J.J. (2000) Efficiency of higher order finite elements for the analysis of seismic wave propagation, Journal of Sound and Vibration, Vol. 231, No. 2, pp Semblat, J.F. (2011) Modeling seismic wave propagation and amplification in 1D/2D/3D linear and nonlinear unbounded media, International Journal of Geomechanics, ASCE, Vol. 11, No. 6, pp Song, C. (2009) The scaled boundary finite element method in structural dynamics, International Journal for Numerical Methods in Engineering, Vol. 77, pp Tabarraei, A. and Sukumar, N. (2005) Adaptive computations on conforming quadtree meshes, Finite Elements in Analysis and Design, Vol. 41, pp Wolf, J.P. and Song, C. (1996) Finite-element modelling of unbounded media, John Wiley & Sons, Chichester, U.K. ACMSM
Automatic image-based stress analysis by the scaled boundary polytope elements
Southern Cross University epublications@scu 3rd Australasian Conference on the Mechanics of Structures and Materials 14 Automatic image-based stress analysis by the scaled boundary polytope elements Chongmin
More informationSpectral-element Simulations of Elastic Wave Propagation in Exploration and Geotechnical Applications
Spectral-element Simulations of Elastic Wave Propagation in Exploration and Geotechnical Applications Lin Zheng*, Qi Zhao*, Qinya Liu, Bernd Milkereit, Giovanni Grasselli, University of Toronto, ON, Canada
More informationTransient analysis of wave propagation in 3D soil by using the scaled boundary finite element method
Southern Cross University epublications@scu 23r Australasian Conference on the Mechanics of Structures an Materials 214 Transient analysis of wave propagation in 3D soil by using the scale bounary finite
More informationEFFICIENCY OF HIGHER ORDER FINITE ELEMENTS FOR THE ANALYSIS OF SEISMIC WAVE PROPAGATION
EFFICIENCY OF HIGHER ORDER FINITE ELEMENTS FOR THE ANALYSIS OF SEISMIC WAVE PROPAGATION Jean-François Semblat, J. J. Brioist To cite this version: Jean-François Semblat, J. J. Brioist. EFFICIENCY OF HIGHER
More informationContinued Fraction Absorbing Boundary Conditions for Transient Elastic Wave Propagation Modeling
Continued Fraction Absorbing Boundary Conditions for Transient Elastic Wave Propagation Modeling Md Anwar Zahid and Murthy N. Guddati 1 Summary This paper presents a novel absorbing boundary condition
More informationFederal Institute for Materials Research and Testing (BAM), Unter den Eichen 87, Berlin, Germany
Jannis Bulling 1, Jens Prager 1, Fabian Krome 1 1 Federal Institute for Materials Research and Testing (BAM), Unter den Eichen 87, 12205 Berlin, Germany Abstract: This paper addresses the computation of
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 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 informationFinite 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 informationcourse outline basic principles of numerical analysis, intro FEM
idealization, equilibrium, solutions, interpretation of results types of numerical engineering problems continuous vs discrete systems direct stiffness approach differential & variational formulation introduction
More informationDynamic Behaviour of Pile Foundations in
Dynamic Behaviour of Pile Foundations in Layered Soil Half-Space Using Cone Method M.M. Aldimashki 1 and J.M.W. Brownjohn 1, Department of Civil and Structural Engineering, University of Sheffield, Sir
More information17. SEISMIC ANALYSIS MODELING TO SATISFY BUILDING CODES
17. SEISMIC ANALYSIS MODELING TO SATISFY BUILDING CODES The Current Building Codes Use the Terminology: Principal Direction without a Unique Definition 17.1 INTRODUCTION { XE "Building Codes" }Currently
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 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 informationITERATIVE LIMIT ANALYSIS OF STRUCTURES WITHIN A SCALED BOUNDARY FINITE ELEMENT FRAMEWORK
Iterative limit analysis of structures within a scaled boundary finite element framework XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII A. Mellati,
More informationVerification and Validation for Seismic Wave Propagation Problems
Chapter 26 Verification and Validation for Seismic Wave Propagation Problems (1989-2-24-25-28-29-21-211-217-) (In collaboration with Dr. Nima Tafazzoli, Dr. Federico Pisanò, Mr. Kohei Watanabe and Mr.
More informationPurdue e-pubs. Purdue University. Jeongil Park Samsung Electronics Co. Nasir Bilal Purdue University. Douglas E. Adams Purdue University
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 26 Development of a Two-Dimensional Finite Element Model of a Suction Valve for Reduction
More informationNOISE PROPAGATION FROM VIBRATING STRUCTURES
NOISE PROPAGATION FROM VIBRATING STRUCTURES Abstract R. Helfrich, M. Spriegel (INTES GmbH, Germany) Noise and noise exposure are becoming more important in product development due to environmental legislation.
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
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 informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationComputer Life (CPL) ISSN: Finite Element Analysis of Bearing Box on SolidWorks
Computer Life (CPL) ISSN: 1819-4818 Delivering Quality Science to the World Finite Element Analysis of Bearing Box on SolidWorks Chenling Zheng 1, a, Hang Li 1, b and Jianyong Li 1, c 1 Shandong University
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 informationRocking Component of Earthquake Induced by Horizontal Motion in Irregular Form Foundation
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 121 Rocking Component of Earthquake Induced by Horizontal Motion in Irregular Form Foundation Tohru Jomen
More informationSEMI-ACTIVE CONTROL OF BUILDING STRUCTURES USING A NEURO-FUZZY CONTROLLER WITH ACCELERATION FEEDBACK
Proceedings of the 6th International Conference on Mechanics and Materials in Design, Editors: J.F. Silva Gomes & S.A. Meguid, P.Delgada/Azores, 26-30 July 2015 PAPER REF: 5778 SEMI-ACTIVE CONTROL OF BUILDING
More informationAircraft Impact Analysis of New York World Trade Center Tower by Using the ASI-Gauss Technique
Proceeding of ICCES 05, 1-10 December 2005, INDIA 1212 Aircraft Impact Analysis of New York World Trade Center Tower by Using the ASI-Gauss Technique D. Isobe 1 and K. M. Lynn 2 Summary In this paper,
More informationA FIXED-BASE MODEL WITH CLASSICAL NORMAL MODES FOR SOIL- STRUCTURE INTERACTION SYSTEMS
A FIXED-BASE MODEL WIH CLASSICAL NORMAL MODES FOR SOIL- SRUCURE INERACION SYSEMS 46 Wen-Hwa WU SUMMARY o simplify the analysis of soil-structure interaction (SSI) systems, different fixed-base models have
More informationA substructure based parallel dynamic solution of large systems on homogeneous PC clusters
CHALLENGE JOURNAL OF STRUCTURAL MECHANICS 1 (4) (2015) 156 160 A substructure based parallel dynamic solution of large systems on homogeneous PC clusters Semih Özmen, Tunç Bahçecioğlu, Özgür Kurç * Department
More informationDeveloper s Tip. An Introduction to the Theory of Planar Failure. Concepts of Planar Rock Slope Failure
Developer s Tip An Introduction to the Theory of Planar Failure In this article we explain the basic concepts behind planar failure analysis of rock slopes. We also discuss the geometric conditions that
More information3DEC 3DEC ITASCA ITASCA VERSION 5.0 VERSION 5.0
ITASCA Consulting Group, Inc. An Itasca International Company 3DEC VERSION 5.0 Advanced, Three Dimensional Discrete Element Modeling for Geotechnical Analysis of Rock, Blocky Structures and Structural
More informationThe geometry of reflection and refraction Wave conversion and reflection coefficient
7.2.2 Reflection and Refraction The geometry of reflection and refraction Wave conversion and reflection coefficient The geometry of reflection and refraction A wave incident on a boundary separating two
More informationNumerical simulation of tunnel construction in volcanic rocks
Numerical simulation of tunnel construction in volcanic rocks G. Beer, Ch. Duenser Graz University of Technology, Graz, Austria ABSTRACT: For the design of tunnels, numerical simulation plays an important
More informationMulticomponent land data pre-processing for FWI: a benchmark dataset
Multicomponent land data pre-processing for FWI: a benchmark dataset Raul Cova, Bernie K. Law and Kris Innanen CRWES/University of Calgary Summary Successful full-waveform inversion (FWI) studies using
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 informationDEVELOPMENT OF PROTOTYPE OF INTEGRATED EARTHQUAKE DISASTER SIMULATOR USING DIGITAL CITY AND STRONG GROUND MOTION SIMULATOR WITH HIGH-RESOLUTION
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 1418 DEVELOPMENT OF PROTOTYPE OF INTEGRATED EARTHQUAKE DISASTER SIMULATOR USING DIGITAL CITY AND STRONG
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 informationChallenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure
Challenge Problem 5 - The Solution Dynamic Characteristics of a Truss Structure In the final year of his engineering degree course a student was introduced to finite element analysis and conducted an assessment
More information3D PHYSICS-BASED NUMERICAL SIMULATIONS: ADVANTAGES AND LIMITATIONS OF A NEW FRONTIER TO EARTHQUAKE GROUND MOTION PREDICTION.
6th National Conference on Earthquake Engineering & 2nd National Conference on Earthquake Engineering and Seismology June 14-16, 2017 3D PHYSICS-BASED NUMERICAL SIMULATIONS: ADVANTAGES AND LIMITATIONS
More informationAn Approximate Method for Permuting Frame with Repeated Lattice Structure to Equivalent Beam
The Open Ocean Engineering Journal, 2011, 4, 55-59 55 Open Access An Approximate Method for Permuting Frame with Repeated Lattice Structure to Equivalent Beam H.I. Park a, * and C.G. Park b a Department
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 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 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 informationEuropean Hyperworks Technology Conference 2010 Mesh and orientation dependance of FE models for dynamic simulations.
European Hyperworks Technology Conference 2010 Mesh and orientation dependance of FE models for dynamic simulations. Sébastien ROTH Jennifer OUDRY Hossein SHAKOURZADEH Influence of mesh density on a finite
More informationAnalysis of Pile Behaviour due to Damped Vibration by Finite Element Method (FEM)
ISSN 2395-1621 Analysis of Pile Behaviour due to Damped Vibration by Finite Element Method (FEM) #1 L.M.Sardesai, #2 G.A.Kadam 1 sardesaileena@rediffmail.com.com 2 ganeshkadam07@gmail.com #1 Mechanical
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 informationExample 24 Spring-back
Example 24 Spring-back Summary The spring-back simulation of sheet metal bent into a hat-shape is studied. The problem is one of the famous tests from the Numisheet 93. As spring-back is generally a quasi-static
More informationFinite Element Buckling Analysis Of Stiffened Plates
International Journal of Engineering Research and Development e-issn: 2278-067X, p-issn: 2278-800X, www.ijerd.com Volume 10, Issue 2 (February 2014), PP.79-83 Finite Element Buckling Analysis Of Stiffened
More information1314. Estimation of mode shapes expanded from incomplete measurements
34. Estimation of mode shapes expanded from incomplete measurements Sang-Kyu Rim, Hee-Chang Eun, Eun-Taik Lee 3 Department of Architectural Engineering, Kangwon National University, Samcheok, Korea Corresponding
More informationTitle. Author(s)P. LATCHAROTE; Y. KAI. Issue Date Doc URL. Type. Note. File Information
Title HIGH PERFORMANCE COMPUTING OF DYNAMIC STRUCTURAL RES INTEGRATED EARTHQUAKE SIMULATION Author(s)P. LATCHAROTE; Y. KAI Issue Date 2013-09-13 Doc URL http://hdl.handle.net/2115/54441 Type proceedings
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationLARGE-SCALE NORTHRIDGE EARTHQUAKE SIMULATION USING OCTREE-BASED MULTIRESOLUTION MESH METHOD
16th ASCE Engineering Mechanics Conference July 16-18, 3, University of Washington, Seattle LARGE-SCALE NORTHRIDGE EARTHQUAKE SIMULATION USING OCTREE-BASED MULTIRESOLUTION MESH METHOD Eui Joong Kim 1 Jacobo
More informationCO2 sequestration crosswell monitoring based upon spectral-element and adjoint methods
CO2 sequestration crosswell monitoring based upon spectral-element and adjoint methods Christina Morency Department of Geosciences, Princeton University Collaborators: Jeroen Tromp & Yang Luo Computational
More informationFluid Structure Interaction - Moving Wall in Still Water
Fluid Structure Interaction - Moving Wall in Still Water Outline 1 Problem description 2 Methodology 2.1 Modelling 2.2 Analysis 3 Finite Element Model 3.1 Project settings 3.2 Units 3.3 Geometry Definition
More information03 Vector Graphics. Multimedia Systems. 2D and 3D Graphics, Transformations
Multimedia Systems 03 Vector Graphics 2D and 3D Graphics, Transformations Imran Ihsan Assistant Professor, Department of Computer Science Air University, Islamabad, Pakistan www.imranihsan.com Lectures
More informationNumerical Treatment of Wave Propagation in Layered Media
Proceedings, Third UJNR Workshop on Soil-Structure Interaction, March 29-30, 2004, Menlo Park, California, USA. Numerical Treatment of Wave Propagation in Layered Media Murthy N. Guddati, a) Si-Hwan Park,
More informationEquivalence between 2D and 3D Numerical Simulations of the Seismic Response of Improved Sites
Equivalence between 2D and 3D Numerical Simulations of the Seismic Response of Improved Sites A. G. Papadimitriou 1, A. C. Vytiniotis 2, G. D. Bouckovalas 3, G. J. Bakas 4 1 Research Associate, National
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 informationSTRUCTURAL TOPOLOGY OPTIMIZATION SUBJECTED TO RELAXED STRESS AND DESIGN VARIABLES
STRUCTURAL TOPOLOGY OPTIMIZATION SUBJECTED TO RELAXED STRESS AND DESIGN VARIABLES Hailu Shimels Gebremedhen, Dereje Engida Woldemichael and Fakhruldin M. Hashim Mechanical Engineering Department, Universiti
More informationTransactions on Modelling and Simulation vol 10, 1995 WIT Press, ISSN X
Hydrodynamic coefficients and motions due to a floating cylinder in waves D.D. Bhatta, M. Rahman Department of Applied Mathematics, Technical University of Nova Scotia, Halifax, Nova Scotia, Canada B3J
More informationEXPLICIT DYNAMIC ANALYSIS OF A REINFORCED CONCRETE FRAME UP TO COLLAPSE
EXPLICIT DYNAMIC ANALYSIS OF A REINFORCED CONCRETE FRAME UP TO COLLAPSE ABSTRACT: K. Kocamaz 1, K. Tuncay 2 and B. Binici 3 1 Res. Assist., Civil Eng. Department, Middle East Technical University, Ankara
More informationMethodology for Prediction of Sliding and Rocking of Rigid Bodies Using Fast Non-Linear Analysis (FNA) Formulation
Methodology for Prediction of Sliding and Rocking of Rigid Bodies Using Fast Non-Linear Analysis (FNA) Formulation Sohrab Esfandiari - ENOVA Engineering Services Robert P. Kennedy- RPK Structural Consulting
More informationA MODELING METHOD OF CURING DEFORMATION FOR CFRP COMPOSITE STIFFENED PANEL WANG Yang 1, GAO Jubin 1 BO Ma 1 LIU Chuanjun 1
21 st International Conference on Composite Materials Xi an, 20-25 th August 2017 A MODELING METHOD OF CURING DEFORMATION FOR CFRP COMPOSITE STIFFENED PANEL WANG Yang 1, GAO Jubin 1 BO Ma 1 LIU Chuanjun
More informationMethod of Finite Elements I
Institute of Structural Engineering Page 1 Held by Prof. Dr. E. Chatzi, Dr. P. Steffen Assistants: Adrian Egger (HIL E 13.3), Harry Mylonas (HIL H33.1), Konstantinos Tatsis (HIL H33.1) Lectures homepage:
More information2-D Tank Sloshing Using the Coupled Eulerian- LaGrangian (CEL) Capability of Abaqus/Explicit
2-D Tank Sloshing Using the Coupled Eulerian- LaGrangian (CEL) Capability of Abaqus/Explicit Jeff D. Tippmann, Sharat C. Prasad 2, and Parthiv N. Shah ATA Engineering, Inc. San Diego, CA 923 2 Dassault
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 informationCHAPTER-10 DYNAMIC SIMULATION USING LS-DYNA
DYNAMIC SIMULATION USING LS-DYNA CHAPTER-10 10.1 Introduction In the past few decades, the Finite Element Method (FEM) has been developed into a key indispensable technology in the modeling and simulation
More informationUSING NEURAL NETWORK FOR PREDICTION OF THE DYNAMIC PERIOD AND AMPLIFICATION FACTOR OF SOIL FOR MICROZONATION IN URBAN AREA
USING NEURAL NETWORK FOR PREDICTION OF THE DYNAMIC PERIOD AND AMPLIFICATION FACTOR OF SOIL FOR MICROZONATION IN URBAN AREA Tahamoli Roodsari, Mehrzad 1 and Habibi, Mohammad Reza 2 1 M.S. in Earthquake
More informationVALIDATION OF SASSI2010 SOLUTION METHODS THROUGH INDEPENDENT VERIFICATION USING SAP2000 FOR DEEPLY EMBEDDED STRUCTURES WITH LARGE FOOTPRINTS
Transactions, SMiRT-22 VALIDATION OF SASSI2010 SOLUTION METHODS THROUGH INDEPENDENT VERIFICATION USING SAP2000 FOR DEEPLY EMBEDDED STRUCTURES WITH LARGE FOOTPRINTS Lisa M Anderson 1, Tarek Elkhoraibi,
More informationAn Introductory SIGMA/W Example
1 Introduction An Introductory SIGMA/W Example This is a fairly simple introductory example. The primary purpose is to demonstrate to new SIGMA/W users how to get started, to introduce the usual type of
More informationInfluence of mesh density on a finite element model under dynamic loading Sébastien ROTH
Influence of mesh density on a finite element model under dynamic loading Sébastien ROTH Jennifer OUDRY Influence of mesh density on a finite element model under dynamic loading 1. Introduction 2. Theoretical
More informationSeismic Soil-Structure Interaction Analysis of the Kealakaha Stream Bridge on Parallel Computers
Seismic Soil-Structure Interaction Analysis of the Kealakaha Stream Bridge on Parallel Computers Seung Ha Lee and Si-Hwan Park Department of Civil and Environmental Engineering University of Hawaii at
More informationTABLE OF CONTENTS SECTION 2 BACKGROUND AND LITERATURE REVIEW... 3 SECTION 3 WAVE REFLECTION AND TRANSMISSION IN RODS Introduction...
TABLE OF CONTENTS SECTION 1 INTRODUCTION... 1 1.1 Introduction... 1 1.2 Objectives... 1 1.3 Report organization... 2 SECTION 2 BACKGROUND AND LITERATURE REVIEW... 3 2.1 Introduction... 3 2.2 Wave propagation
More informationPerformance of railway track system under harmonic loading by finite element method
Performance of railway track system under harmonic by finite element method Ammar Shuber 1, Mohammed Hamood 1, * and Walaa Jawad 1 1 Building and Construction Engineering Department, University of Technology,
More informationA STUDY ON THE INFLUENCE OF BOUNDARY CONDITIONS ON THREE DIMENSIONAL DYNAMIC RESPONSES OF EMBANKMENTS
October 12-17, 28, Beijing, China A STUDY ON THE INFLUENCE OF BOUNDARY CONDITIONS ON THREE DIMENSIONAL DYNAMIC RESPONSES OF EMBANKMENTS S. Kano 1 a, Y. Yokoi 2, Y. Sasaki 3 and Y. Hata 4 1 Assistant Professor,
More informationTransactions on Information and Communications Technologies vol 15, 1997 WIT Press, ISSN
Optimization of time dependent adaptive finite element methods K.-H. Elmer Curt-Risch-Institut, Universitat Hannover Appelstr. 9a, D-30167 Hannover, Germany Abstract To obtain reliable numerical solutions
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 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 informationGlobal to Local Model Interface for Deepwater Top Tension Risers
Global to Local Model Interface for Deepwater Top Tension Risers Mateusz Podskarbi Karan Kakar 2H Offshore Inc, Houston, TX Abstract The water depths from which oil and gas are being produced are reaching
More informationVisualization of errors of finite element solutions P. Beckers, H.G. Zhong, Ph. Andry Aerospace Department, University of Liege, B-4000 Liege, Belgium
Visualization of errors of finite element solutions P. Beckers, H.G. Zhong, Ph. Andry Aerospace Department, University of Liege, B-4000 Liege, Belgium Abstract The aim of this paper is to show how to use
More informationDesign Optimization of Building Structures Using a Metamodeling Method
Proceedings of the 3 rd International Conference on Civil, Structural and Transportation Engineering (ICCSTE'18) Niagara Falls, Canada June 10 12, 2018 Paper No. 137 DOI: 10.11159/iccste18.137 Design Optimization
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More informationExcerpt from the Proceedings of the COMSOL Conference 2010 Paris
Excerpt from the Proceedings of the COMSOL Conference 2010 Paris Simulation of Flaw Signals in a Magnetic Flux Leakage Inspection Procedure O. Nemitz * and T. Schmitte Salzgitter Mannesmann Forschung GmbH
More informationCoustyx Tutorial Indirect Model
Coustyx Tutorial Indirect Model 1 Introduction This tutorial is created to outline the steps required to compute radiated noise from a gearbox housing using Coustyx software. Detailed steps are given on
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 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 informationEffect of Diagonal Modes on Response Spectrum Analysis
Effect of Diagonal Modes on Response Spectrum Analysis T. Naga Manikanta & O. R. Jaiswal Department of Applied Mechanics Visvesvaraya National Institute of Technology, Nagpur, Maharashtra, India Summary:
More informationParallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of Earth s Mantle
ICES Student Forum The University of Texas at Austin, USA November 4, 204 Parallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of
More informationAbaqus Technology Brief. Sound Radiation Analysis of Automobile Engine Covers
Sound Radiation Analysis of Automobile Engine Covers Abaqus Technology Brief TB-06-COVER-2 Revised: April 2007. Summary A methodology to study the sound radiation of engine valve covers is presented. The
More information9. Three Dimensional Object Representations
9. Three Dimensional Object Representations Methods: Polygon and Quadric surfaces: For simple Euclidean objects Spline surfaces and construction: For curved surfaces Procedural methods: Eg. Fractals, Particle
More informationThree Dimensional Seismic Simulation of a Two-Bay, Two-Story Reinforced Concrete Building
Three Dimensional Seismic Simulation of a Two-Bay, Two-Story Reinforced Concrete Building Catherine Joseph, REU Student Rachel Howser, Graduate Advisor Dr. Y.L Mo, Faculty Advisor Final Report Department
More informationDYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION
DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION 7 DYNAMIC ANALYSIS OF A GENERATOR ON AN ELASTIC FOUNDATION In this tutorial the influence of a vibrating source on its surrounding soil is studied.
More informationExercise 1. 3-Point Bending Using the GUI and the Bottom-up-Method
Exercise 1 3-Point Bending Using the GUI and the Bottom-up-Method Contents Learn how to... 1 Given... 2 Questions... 2 Taking advantage of symmetries... 2 A. Preprocessor (Setting up the Model)... 3 A.1
More informationModeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics
Modeling the Acoustic Scattering from Axially Symmetric Fluid, Elastic, and Poroelastic Objects due to Nonsymmetric Forcing Using COMSOL Multiphysics Anthony L. Bonomo *1 and Marcia J. Isakson 1 1 Applied
More informationFinite Element Analysis of Dynamic Flapper Valve Stresses
Purdue University Purdue e-pubs International Compressor Engineering Conference School of Mechanical Engineering 2000 Finite Element Analysis of Dynamic Flapper Valve Stresses J. R. Lenz Tecumseh Products
More informationEditorial Coupled Numerical Methods in Engineering Analysis
Mathematical Problems in Engineering Volume 2011, Article ID 436978, 4 pages doi:10.1155/2011/436978 Editorial Coupled Numerical Methods in Engineering Analysis Delfim Soares Jr., 1 Otto von Estorff, 2
More informationCOMSOL BASED 2-D FEM MODEL FOR ULTRASONIC GUIDED WAVE PROPAGATION IN SYMMETRICALLY DELAMINATED UNIDIRECTIONAL MULTI- LAYERED COMPOSITE STRUCTURE
Proceedings of the National Seminar & Exhibition on Non-Destructive Evaluation NDE 2011, December 8-10, 2011 COMSOL BASED 2-D FEM MODEL FOR ULTRASONIC GUIDED WAVE PROPAGATION IN SYMMETRICALLY DELAMINATED
More informationCFD modelling of thickened tailings Final project report
26.11.2018 RESEM Remote sensing supporting surveillance and operation of mines CFD modelling of thickened tailings Final project report Lic.Sc.(Tech.) Reeta Tolonen and Docent Esa Muurinen University of
More information3D Finite Difference Time-Domain Modeling of Acoustic Wave Propagation based on Domain Decomposition
3D Finite Difference Time-Domain Modeling of Acoustic Wave Propagation based on Domain Decomposition UMR Géosciences Azur CNRS-IRD-UNSA-OCA Villefranche-sur-mer Supervised by: Dr. Stéphane Operto Jade
More informationLiquefaction Analysis in 3D based on Neural Network Algorithm
Liquefaction Analysis in 3D based on Neural Network Algorithm M. Tolon Istanbul Technical University, Turkey D. Ural Istanbul Technical University, Turkey SUMMARY: Simplified techniques based on in situ
More informationSimulation of Curtain Airbag with Arbitrary. Eulerian- Lagrangian Method
Simulation of Curtain Airbag with Arbitrary Eulerian- Lagrangian Method Dmitri Fokin, Eivind Dessarud Altair Engineering Claes Ljungqvist Saab Automobile AB Abstract: Computer simulation is a powerful
More information