An alternative multi-region BEM technique for layered soil problems

Size: px
Start display at page:

Download "An alternative multi-region BEM technique for layered soil problems"

Transcription

1 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 formulations may be found in the literature for the numerical simulation of layered soils. Among these options, the boundary element method (BM) detaches because it requires only the boundary discretisation of the infinite domains involved. The classical way of simulating infinite layered solids with the BM is by imposing equilibrium and displacements compatibility along the contacts. There is, however, an alternative technique which is more accurate because it does not impose these conditions along the interfaces. In this alternative technique, the regions displacement fundamental solutions are related and for each equation all domains are considered as a unique solid. In this chapter, it is demonstrated that this alternative formulation is suited for the simulation of layered soil problems, representing such media as an infinite solid composed of different homogeneous, isotropic and elastic linear domains in contact. In the end of the chapter, this formulation is employed in a threedimensional problem, comparing the results with the classical technique and with an analytical solution. Keywords: Boundary element method, Layered soil, Alternative multi-region technique, Static. Introduction In the literature, many works are dedicated to the simulation of layered soil problems. In such works, when a structure is considered interacting with the soil, WIT Transactions on State of the Art in Science and ngineering, Vol 43, 00 WIT Press ISSN (on-line) doi:0.495/ /

2 6 Recent Developments in Boundary lement Methods one may observe that the mathematical models employed for the structures are closer to reality than the ones used for the soil. This discrepancy is due to the type of materials involved in each case. While structural materials, such as concrete and steel, are produced with some type of technological control, the soil materials are formed due to various conditions of weathering in rocks. The result is a chaotic media where the characteristics are difficult to be numerically modelled, including discontinuities, heterogeneity and anisotropy. Due to these factors, it is necessary to assume simplifications in order to make possible a numerical simulation of layered soil problems with reasonable time processing. In this chapter, the main simplifications adopted are: Representing the soil as an infinite domain in radial directions; Considering the soil as a composition of homogeneous, isotropic, and linear elastic domains; Considering only the static case. In this case, the problem comes down to the static analysis of infinite multiregion domains. Focusing this theme, many formulations are available in the literature and each one of them implies on advantages and disadvantages. However, depending on the problem to be solved, one technique may become more attractive than the others. When it is possible, a good choice is to employ analytical solutions. When correctly programmed, they give trustful results in little processing time. In reference [], for example, an analytical expression is deduced for a circular load applied to the surface of a two-layer infinite domain. The main disadvantage of these solutions is that they are not versatile, suiting only specific situations. If analytical solutions can not be used, then a numerical approach may become attractive. Although the finite element method (FM) is popular, it has some disadvantages compared with other options such as the boundary element method (BM). The FM requires the discretisation of the infinite domain, implying on a high number of elements and leading to a large and, sometimes, impracticable processing time. In order to reduce these inconveniences, some authors use infinite elements together with finite elements, as done in reference []. The main advantage of the BM is that only the boundaries of the domains involved require discretisation. This allows reducing the problem dimension, implying on less processing time. This advantage is explored in several works, and more developments are turning the BM even more attractive to future applications. One of them is using infinite boundary elements, as performed in reference [3]. To consider two or more domains in contact with the BM, several techniques may be employed. The most popular one, which may be consulted in reference [4], is based on imposing equilibrium and compatibility conditions for all interface points between every pair of domains in contact. Using these relations, the matrices obtained for each domain applying the BM formulation WIT Transactions on State of the Art in Science and ngineering, Vol 43, 00 WIT Press ISSN (on-line)

3 Recent Developments in Boundary lement Methods 63 are coupled and all the unknown boundary values may be obtained by solving the same system of equations. In spite of its simplicity, one disadvantage of this formulation is that the final system of equations has blocks of zeros. These blocks can become numerous depending on the number of domains considered, increasing the size of the system and consequently the processing time. In addition to that, the conditions imposed on the contacts do not represent satisfactorily the media continuum, what may cause inaccuracies on the results. The multi-region method presented in reference [5] for two-dimensional elastic and potential problems eliminates the need of equilibrium and compatibility conditions among the interfaces. This approach improves the media continuity compared with the classical formulation, leading to more accurate results. This technique is modified in reference [6] for bending plate analysis with the BM, latter used for bending moments calculation in [7] and finally adapted for three-dimensional elastic problems in reference [8]. The aim of this chapter is to present this alternative multi-region BM technique, which is suitable for layered soil problems. The soil is represented as a multi-region media, similar to the Gibson soil approach. In such a way, the domain is modelled with variable elasticity module and a constant Poisson ratio, which may be considered a disadvantage in certain cases. Nevertheless, by testing this formulation in problems with different Poisson ratios, one may verify that the error introduced by an average Poisson ratio consideration may be considered of little relevance for displacement calculation. In the end, it is viable to employ this formulation in more general problems. Alternative multi-region formulation In Figure, a homogeneous, isotropic and linear elastic domain Ω is presented. The boundary is denoted by, is the elasticity module and ν is the Poisson Figure : Problem with one region. WIT Transactions on State of the Art in Science and ngineering, Vol 43, 00 WIT Press ISSN (on-line)

4 64 Recent Developments in Boundary lement Methods ratio. Point x is called field point, and represents a general point of. Point y is called source point, and can be placed in any point inside Ω, outside Ω or at the boundary. Versor η is normal to at x, and r is the distance between points x and y. The equilibrium of this solid body can be represented by a boundary integral equation called Somigliana identity, which for homogeneous, isotropic and linear-elastic domains is: () ij j ij j ij j c ( y) u ( y) + p ( x, y) u ( x)d ( x) = u ( x, y) p ( x)d ( x). This equation is written for the source point y, where the displacement is uj ( y ). The constant c ij depends on the Poisson ratio and on the position of y. The field point x goes through all boundary, where displacements are uj ( x ) and tractions are pj ( x ). The integral kernels uij ( x, y ) and pij ( x, y ) are Kelvin three-dimensional fundamental solutions for displacements and tractions, respectively, and may be written as follows: u = [(3ν 4) δ + r, r, ], 6 πµ ( ν ) r ij ij i j r pij = [( ) 3,, ] ( )( νδij + rir j + ν ηir, j η jr, i ). 8 π( ν) r η () (3) In eqns () and (3), δ ij is unitary for i = j and zero for all other cases, η is the normal versor at point x (see Figure ) and µ is given by: µ = ( + ν ). (4) The objective is to obtain an integral equation similar to eqn (), but valid for an arbitrary number of domains in contact. In such a way, a demonstration will be held for the most simple case, with only two regions as illustrated in Figure, and the resulting expression will then be extended for an arbitrary number of domains. In Figure, the regions have the same Poisson ratio ν and different elasticity modules, for region Ω and for region Ω. The boundary of region Ω is divided into two parts, and. The part of which is in contact with region Ω is denoted by and the rest of it is named. Analogously, for region Ω, the boundary is divided into for the contact and for the free surface. Consequently: = +, = +. (5) WIT Transactions on State of the Art in Science and ngineering, Vol 43, 00 WIT Press ISSN (on-line)

5 Recent Developments in Boundary lement Methods 65 Kelvin displacement fundamental solutions for regions Ω and Ω may be written as: where u u = [(3ν 4) δ + r r ], 6 πµ ( ν ) r ij ij, i, j = [(3ν 4) δ + r r ], 6 πµ ( ν ) r ij ij, i, j (6) (7) µ =, µ =. ( + ν ) ( + ν ) (8) In such a way, the fundamental solutions may be related as follows: where u = u u = u + u, (9) ij ij ij ij ij ij = i j. (0) One may observe in eqn (3) that the traction fundamental solution does not depend on the elasticity module, except when it is the same for domains Ω and Ω. Thus, Kelvin traction fundamental solution may be represented as p ij for both domains. Considering region Ω as a reference, the deduction is started by writing the Somigliana identity, represented in eqn (). In order to reduce the expressions, the variables in parenthesis will be suppressed. In such a way: () ij j + ij j = ij j c u p u d u p d. Note that if the source point is outside Ω, then the free-term c ij is zero; if it is inside, c ij = and if it is on the boundary, then c ij may be determined using standard BM techniques. Considering eqn (5), the integrals of eqn () may be divided into two parts as shown below: () ij j + ij j + ij j = ij j + ij j c u p u d p u d u p d u p d. An analogous equation may be written for region Ω : (3) ij j + ij j + ij j = ij j + ij j c u pud pud u pd u pd. WIT Transactions on State of the Art in Science and ngineering, Vol 43, 00 WIT Press ISSN (on-line)

6 66 Recent Developments in Boundary lement Methods When eqns () and (3) are added, a single equation which considers both regions is obtained for the point: ( cij + cij ) uj + pijujd + pijujd + pijujd + pijujd = uijpjd + uijpjd uijpjd uijpjd. + + (4) Some terms of eqn (4) may be related. Starting with the left side, the integral over is equal to the one over with an inverted signal. This occurs because the exactly same functions are being integrated only by inverting the integration direction. Thus: d + d = 0. (5) pu ij j pu ij j Analysing now the right side of eqn (4) and considering eqn (9), the following arrangement is made: uijpjd + uijpjd = uijpjd ( ) + + = + uij uij pjd uij pjd uij pjd + d = d. uijpj uijpj (6) The two integrals added inside the parenthesis are equal to zero for the same reason that eqn (5) is zero. One more simplification is possible on the right side of eqn (4) by substituting eqn (9) in one more term, as shown below: uijpjd = uij + uij pjd = uijpjd + d. uij pj After these deductions, eqns (5) (7) are substituted in eqn (4). The result is ( cij + cij ) uj + pijujd + pijujd = uij pjd uijpj uijpj uijpj + d + d + d. Reorganising the terms, the following equation is obtained: ( cij + cij ) uj + pijujd + pijujd = uij pjd uijpj ( uijpj uijpj ) + d + d + d. (7) (8) (9) WIT Transactions on State of the Art in Science and ngineering, Vol 43, 00 WIT Press ISSN (on-line)

7 Recent Developments in Boundary lement Methods 67 The terms inside the parenthesis are functions of the boundary tractions; however, it is necessary to transform them into functions of the boundary displacements. This transformation is possible using eqn (3), isolating one of its integrals for further substitution in eqn (9). In such a way: d d d d. (0) uijpj = cijuj + pijuj + pijuj uijpj Substituting eqn (0) into eqn (9) and combining the integrals inside the parenthesis, the subsequent equation is obtained: ( cij + cij ) uj + pijujd + pijujd = uij pjd uijpj cij pijuj pijuj ( ) + d + + d + d. ( ) c c u p u p u ij + ij j ij jd ij jd + + pijujd + pijujd uijpjd uijpjd. = + Note that two integrals on the left side can still be combined as follows: () () pu ij jd pu ij jd = pu ij jd = pujd. ij The result of eqn (3) substitution into eqn () is: c c u p u p u p u ij + ij j ij jd ij jd ij jd + + uijpj uijpj = d + d. (3) (4) quation (4) is only valid for two domains. xtending it to an arbitrary number of domains, it becomes: nd ne nc s e mn cijs uj pijujd e pijujd mn e mn s + e + c = = = ne = uij pjd e. e e= (5) In eqn (5), the second summation signal is positive because the integration direction was changed from nm to mn. The total number of domains is nd, WIT Transactions on State of the Art in Science and ngineering, Vol 43, 00 WIT Press ISSN (on-line)

8 68 Recent Developments in Boundary lement Methods Table : xamples of c ij. Point position c ij c ij c ij Outside Ω and Ω Inside Ω 0 Over boundary 0 ( ) Over boundary + the number of contact boundaries is nc and the number of external boundaries is ne. The first summation represents the coefficient cij ( y ) of eqn (). Before calculating it, all coefficients c ijs, one for each domain, must be known. Considering all boundaries smooth in Figure, Table presents some examples for this summation calculation. quation () is used as a starting point to obtain the BM system of equations which solution leads to the unknowns of the problem. If the same steps are repeated for eqn (5), valid for multi-regions, a similar system of equations is obtained. The unknowns of this new system are the non-prescribed boundary values plus the interface displacements. The total number of unknowns is reduced compared with the classic multi-region technique, once the interface tractions are not included in this case. This justifies why the alternative formulation leads to less time processing. A better interface continuity is also guaranteed, once all regions are modelled as a unique solid. 3 Results This example aims to analyse an infinite multi-region domain problem with the alternative BM technique. The domain considered is composed by two layers of different elasticity module, as illustrated in Figure 3. Layer has a 9000 kn m elasticity module, a 0.5 Poisson ratio and 5m of thickness. Layer has a 900 kn m elasticity module, a 0.5 Poisson ratio and infinite thickness. Both layers are infinite in radial directions. A vertical circular kn m uniform loading with a 7.5m diameter is applied at the top layer surface. This problem was simulated using a 53 node and 88 boundary element mesh, which may be visualised in Figure 4. It is composed by triangular elements with linear shape functions. Figure 4(a) presents a general view of the mesh and Figure 4(b) presents a detail of the central area. This mesh is employed at the surface and at the layers contact, extending to a distance from which the displacements and tractions could be considered negligible. For the nodes positioned at this WIT Transactions on State of the Art in Science and ngineering, Vol 43, 00 WIT Press ISSN (on-line)

9 Recent Developments in Boundary lement Methods 69 Figure 3: Problem to be analysed. limits, the boundary values are imposed to be zero to better simulate their far field behaviour. Hence, it was considered not necessary to close the boundary at the limits. As commented in the introduction of this chapter, in reference [], an analytical solution is deduced for a circular load applied to the surface of a two-layer infinite domain. Using the values adopted in this example and applying the solution given in reference [], a m vertical displacement is obtained for the central point of the circle. Considering this same point, a m displacement was obtained using the classical formulation with the mesh of Figure 4. mploying the alternative technique and the same mesh, a m vertical displacement was obtained. Both Figure 4: Mesh employed. WIT Transactions on State of the Art in Science and ngineering, Vol 43, 00 WIT Press ISSN (on-line)

10 70 Recent Developments in Boundary lement Methods values numerically calculated may be considered very close to the analytical one, leading to the conclusion that both formulations achieved accurate results in this example. 4 Conclusions In this chapter, an alternative BM formulation suited for layered soil problems was presented. lecting one domain as a reference and establishing relations between its displacement fundamental solution and the ones of the other regions allows integrating all domains as a unique solid. This approach eliminates the need of equilibrium and compatibility relations between the different media, guaranteeing a better continuity along the interfaces. This reduces possible inaccuracies that may occur due to traction approximations along the contacts. A numerical example was presented, in which a two-layer half space was simulated employing the classical and the alternative formulations. The values obtained were compared with an analytical expression. The numerical results were very close to the analytical one, meaning that both techniques were accurate. These low errors were expected, once the problem analysed may be considered simple. References [] Burmister, D.M., The general theory of stresses and displacements in layered systems. I. Journal of Applied Physics, 6, pp , 945. [] Seo, C.G., Yun, C.B. & Kim, J.M., Three-dimensional frequency-dependent infinite elements for soil structure interaction. ngineering Structures, 9, pp , 007. [3] Moser, W., Duenser, C. & Beer, G., Mapped infinite elements for threedimensional multi-region boundary element analysis. International Journal for Numerical Methods in ngineering, 6, pp , 004. [4] Brebbia, C.A. & Dominguez, J., Boundary lements: An Introductory Course, Computational Mechanics Publications: London, 99. [5] Venturini, W.S., Alternative formulations of the boundary element method for potential and elastic problems, ngineering Analysis with Boundary lements, 9, pp , 99. [6] Paiva, J.B. & Aliabadi, M.H., Boundary element analysis of zoned plates in bending, Computational Mechanics, 5, pp , 000. [7] Paiva, J.B. & Aliabadi, M.H., Bending moments at interfaces of thin zoned plates with discrete thickness by the boundary element method. ngineering Analysis with Boundary lements, 8, pp , 004. [8] Ribeiro, D.B. & Paiva, J.B., An alternative multi-region BM technique for three-dimensional elastic problems. ngineering Analysis with Boundary lements, 33, pp , 009. WIT Transactions on State of the Art in Science and ngineering, Vol 43, 00 WIT Press ISSN (on-line)

Transactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X

Transactions 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 information

Plate analysis using classical or Reissner- Mindlin theories

Plate 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 information

Department of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil

Department 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 information

L e c t u r e. D r. S a s s a n M o h a s s e b

L e c t u r e. D r. S a s s a n M o h a s s e b The Scaled smteam@gmx.ch Boundary Finite www.erdbebenschutz.ch Element Method Lecture A L e c t u r e A1 D r. S a s s a n M o h a s s e b V i s i t i n g P r o f e s s o r M. I. T. C a m b r i d g e December

More information

Transactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X

Transactions 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 information

COMPUTATIONAL ELASTICITY

COMPUTATIONAL 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 information

1 Numerical combination for nonlinear analysis of structures coupled to layered soils

1 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 information

SPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES

SPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES 2010/2 PAGES 1 8 RECEIVED 21. 9. 2009 ACCEPTED 20. 1. 2010 Y. KOLEKOVÁ, M. PETRONIJEVIĆ, G. SCHMID SPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES ABSTRACT

More information

Mixed FEM-BEM Formulations Applied to Soil-Structure Interaction Problems

Mixed FEM-BEM Formulations Applied to Soil-Structure Interaction Problems Mixed FEM-BEM Formulations Applied to Soil-Structure Interaction Problems Dimas Betioli Ribeiro and João Batista de Paiva Abstract The objective of this paper is to present formulations developed for soil-building

More information

Chapter 4 Analysis of a cantilever

Chapter 4 Analysis of a cantilever Chapter 4 Analysis of a cantilever Before a complex structure is studied performing a seismic analysis, the behaviour of simpler ones should be fully understood. To achieve this knowledge we will start

More information

MATERIAL PROPERTIES. Material Properties Must Be Evaluated By Laboratory or Field Tests 1.1 INTRODUCTION 1.2 ANISOTROPIC MATERIALS

MATERIAL PROPERTIES. Material Properties Must Be Evaluated By Laboratory or Field Tests 1.1 INTRODUCTION 1.2 ANISOTROPIC MATERIALS . MARIAL PROPRIS Material Properties Must Be valuated By Laboratory or Field ests. INRODUCION he fundamental equations of structural mechanics can be placed in three categories[]. First, the stress-strain

More information

Discrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method

Discrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method 131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using

More information

Fast multipole boundary element method for the analysis of plates with many holes

Fast multipole boundary element method for the analysis of plates with many holes Arch. Mech., 59, 4 5, pp. 385 401, Warszawa 2007 Fast multipole boundary element method for the analysis of plates with many holes J. PTASZNY, P. FEDELIŃSKI Department of Strength of Materials and Computational

More information

Civil Engineering Design (1) Analysis and Design of Slabs 2006/7

Civil Engineering Design (1) Analysis and Design of Slabs 2006/7 Civil Engineering Design (1) Analysis and Design of Slabs 006/7 Dr. Colin Caprani, Chartered Engineer 1 Contents 1. Elastic Methods... 3 1.1 Introduction... 3 1. Grillage Analysis... 4 1.3 Finite Element

More information

Fundamentals of Linear Elasticity

Fundamentals of Linear Elasticity Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy

More information

Transactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X

Transactions on Modelling and Simulation vol 9, 1995 WIT Press,   ISSN X Elastic-plastic model of crack growth under fatigue using the boundary element method M. Scibetta, O. Pensis LTAS Fracture Mechanics, University ofliege, B-4000 Liege, Belgium Abstract Life of mechanic

More information

Mechanics PhD Preliminary Spring 2017

Mechanics PhD Preliminary Spring 2017 Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n

More information

CRITERIA FOR SELECTION OF FEM MODELS.

CRITERIA 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 information

Constitutive Equations

Constitutive Equations Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module

More information

1 Introduction. Abstract

1 Introduction. Abstract Abstract This paper presents a three-dimensional numerical model for analysing via finite element method (FEM) the mechanized tunneling in urban areas. The numerical model is meant to represent the typical

More information

Integral equations for crack systems in a slightly heterogeneous elastic medium

Integral equations for crack systems in a slightly heterogeneous elastic medium Boundary Elements and Other Mesh Reduction Methods XXXII 65 Integral equations for crack systems in a slightly heterogeneous elastic medium A. N. Galybin & S. M. Aizikovich Wessex Institute of Technology,

More information

PRACTICE 2 PROYECTO Y CONSTRUCCIÓN DE PUENTES. 1º Máster Ingeniería de Caminos. E.T.S.I. Caminos, canales y puertos (Ciudad Real) 01/06/2016

PRACTICE 2 PROYECTO Y CONSTRUCCIÓN DE PUENTES. 1º Máster Ingeniería de Caminos. E.T.S.I. Caminos, canales y puertos (Ciudad Real) 01/06/2016 PRACTICE 2 PROYECTO Y CONSTRUCCIÓN DE PUENTES 1º Máster Ingeniería de Caminos E.T.S.I. Caminos, canales y puertos (Ciudad Real) 01/06/2016 AUTHOR: CONTENT 1. INTRODUCTION... 3 2. BRIDGE GEOMETRY AND MATERIAL...

More information

With appropriate conditions from constitutive relations, the stress distribution can be found by fitting the boundary

With appropriate conditions from constitutive relations, the stress distribution can be found by fitting the boundary Lecture Note 3 Stress Functions First Semester, Academic Year Department of Mechanical ngineering Chulalongkorn Universit Objectives # Describe the equation of equilibrium using D/3D stress elements Set

More information

Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method

Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method V.M.F. Nascimento Departameto de ngenharia Mecânica TM, UFF, Rio de Janeiro

More information

An accelerated predictor-corrector scheme for 3D crack growth simulations

An accelerated predictor-corrector scheme for 3D crack growth simulations An accelerated predictor-corrector scheme for 3D crack growth simulations W. Weber 1 and G. Kuhn 2 1,2 1 Institute of Applied Mechanics, University of Erlangen-Nuremberg Egerlandstraße 5, 91058 Erlangen,

More information

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Fig. 1. Circular fiber and interphase between the fiber and the matrix. Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In

More information

Stresses and Strains in flexible Pavements

Stresses and Strains in flexible Pavements Stresses and Strains in flexible Pavements Multi Layered Elastic System Assumptions in Multi Layered Elastic Systems The material properties of each layer are homogeneous property at point A i is the same

More information

Analysis of Rectangular Plate with Opening by Finite Difference Method

Analysis of Rectangular Plate with Opening by Finite Difference Method American Journal of Civil Engineering and Architecture, 2015, Vol. 3, No. 5, 165-173 Available online at http://pubs.sciepub.com/ajcea/3/5/3 Science and Education Publishing DOI:10.12691/ajcea-3-5-3 Analysis

More information

Using MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,

Using MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup, Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business

More information

Assessment of boundaryelement method for modelling the structural response of a pavement

Assessment of boundaryelement method for modelling the structural response of a pavement Assessment of boundaryelement method for modelling the structural response of a pavement (PhD Student, Lecturer DEC of FCTUC) (Full Professor - IST) Pavements: Materials, design and performance a, LNEC,

More information

Beam Models. Wenbin Yu Utah State University, Logan, Utah April 13, 2012

Beam Models. Wenbin Yu Utah State University, Logan, Utah April 13, 2012 Beam Models Wenbin Yu Utah State University, Logan, Utah 843-4130 April 13, 01 1 Introduction If a structure has one of its dimensions much larger than the other two, such as slender wings, rotor blades,

More information

Elastodynamic contact in plate-like bodies

Elastodynamic contact in plate-like bodies Structures Under Shock and Impact IX 65 Elastodynamic contact in plate-like bodies R. Sburlati Department of Structural and Geotechnical Engineering, Genoa University, Italy Abstract In this paper the

More information

Figure 2-1: Stresses under axisymmetric circular loading

Figure 2-1: Stresses under axisymmetric circular loading . Stresses in Pavements.1. Stresses in Fleible Pavements.1.1. Stresses in Homogeneous Mass Boussinesq formulated models for the stresses inside an elastic half-space due to a concentrated load applied

More information

Module 9 : Foundation on rocks

Module 9 : Foundation on rocks LECTURE 32 9.3 BEARING CAPCITY contd... 9.3.2 Safe bearing pressure There are different methods are available to determine the safe bearing pressure on rocks. However, the applicability of different methods

More information

Hankel Tranform Method for Solving Axisymmetric Elasticity Problems of Circular Foundation on Semi-infinite Soils

Hankel Tranform Method for Solving Axisymmetric Elasticity Problems of Circular Foundation on Semi-infinite Soils ISSN (Print) : 19-861 ISSN (Online) : 975-44 Charles Chinwuba Ie / International Journal of Engineering and Technology (IJET) Hanel Tranform Method for Solving Axisymmetric Elasticity Problems of Circular

More information

Lateral-torsional buckling of orthotropic rectangular section beams

Lateral-torsional buckling of orthotropic rectangular section beams ateral-torsional buckling of orthotropic rectangular section beams Summary Bambang SUYOAMONO Professor Department of Civil ngineering, Parahyangan Catholic University Bandung, Indonesia Adhijoso JONDO

More information

DESIGN OF LAMINATES FOR IN-PLANE LOADING

DESIGN OF LAMINATES FOR IN-PLANE LOADING DESIGN OF LAMINATES FOR IN-PLANOADING G. VERCHERY ISMANS 44 avenue F.A. Bartholdi, 72000 Le Mans, France Georges.Verchery@m4x.org SUMMARY This work relates to the design of laminated structures primarily

More information

Finite Element Method in Geotechnical Engineering

Finite 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 information

Practice Final Examination. Please initial the statement below to show that you have read it

Practice Final Examination. Please initial the statement below to show that you have read it EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use

More information

Lecture 2: Finite Elements

Lecture 2: Finite Elements Materials Science & Metallurgy Master of Philosophy, Materials Modelling, Course MP7, Finite Element Analysis, H. K. D. H. Bhadeshia Lecture 2: Finite Elements In finite element analysis, functions of

More information

Analytical Strip Method for Thin Isotropic Cylindrical Shells

Analytical Strip Method for Thin Isotropic Cylindrical Shells IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 14, Issue 4 Ver. III (Jul. Aug. 2017), PP 24-38 www.iosrjournals.org Analytical Strip Method for

More information

University of Illinois at Urbana-Champaign College of Engineering

University 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 information

A short review of continuum mechanics

A short review of continuum mechanics A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material

More information

EFFECTS OF RESERVOIR LENGTH ON DYNAMIC ANALYSIS OF CONCRETE GRAVITY DAMS

EFFECTS OF RESERVOIR LENGTH ON DYNAMIC ANALYSIS OF CONCRETE GRAVITY DAMS The th October -,, Beijing, China EFFECTS OF RESERVOIR LENGTH ON DYNAMIC ANALYSIS OF CONCRETE GRAVITY DAMS A. Fathi and V. Lotfi M.Sc. Student, Dept. of Civil and Environmental Engineering, Amirkabir University

More information

Longitudinal buckling of slender pressurised tubes

Longitudinal buckling of slender pressurised tubes Fluid Structure Interaction VII 133 Longitudinal buckling of slender pressurised tubes S. Syngellakis Wesse Institute of Technology, UK Abstract This paper is concerned with Euler buckling of long slender

More information

EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE OF A RECTANGULAR ELASTIC BODY MADE OF FGM

EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE OF A RECTANGULAR ELASTIC BODY MADE OF FGM Proceedings of the International Conference on Mechanical Engineering 2007 (ICME2007) 29-31 December 2007, Dhaka, Bangladesh ICME2007-AM-76 EFFECTS OF THERMAL STRESSES AND BOUNDARY CONDITIONS ON THE RESPONSE

More information

Biaxial Analysis of General Shaped Base Plates

Biaxial Analysis of General Shaped Base Plates Biaxial Analysis of General Shaped Base Plates R. GONZALO ORELLANA 1 Summary: A linear model is used for the contact stresses calculation between a steel base plate and a concrete foundation. It is also

More information

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES

FLEXIBILITY METHOD FOR INDETERMINATE FRAMES UNIT - I FLEXIBILITY METHOD FOR INDETERMINATE FRAMES 1. What is meant by indeterminate structures? Structures that do not satisfy the conditions of equilibrium are called indeterminate structure. These

More information

Seismic Analysis of Structures Prof. T.K. Datta Department of Civil Engineering Indian Institute of Technology, Delhi

Seismic Analysis of Structures Prof. T.K. Datta Department of Civil Engineering Indian Institute of Technology, Delhi Seismic Analysis of Structures Prof. T.K. Datta Department of Civil Engineering Indian Institute of Technology, Delhi Lecture - 20 Response Spectrum Method of Analysis In the last few lecture, we discussed

More information

Chapter 5 Structural Elements: The truss & beam elements

Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 1 Chapter 5 Structural Elements: The truss & beam elements Institute of Structural Engineering Page 2 Chapter Goals Learn how to formulate the Finite Element Equations

More information

Exercise: concepts from chapter 8

Exercise: concepts from chapter 8 Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic

More information

Engineering. Lectu re Notes in. W S. Venturini. Boundary Element Method in Geomechanics. Springer-Verlag Berlin Heidelberg New York Tokyo 1983

Engineering. Lectu re Notes in. W S. Venturini. Boundary Element Method in Geomechanics. Springer-Verlag Berlin Heidelberg New York Tokyo 1983 Lectu re Notes in Engineering Edited by C.A. Brebbia and S.A. Orszag 4 W S. Venturini Boundary Element Method in Geomechanics Springer-Verlag Berlin Heidelberg New York Tokyo 1983 Series Editors C.A Brebbia.

More information

DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION

DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,

More information

Dynamic behavior of turbine foundation considering full interaction among facility, structure and soil

Dynamic behavior of turbine foundation considering full interaction among facility, structure and soil Dynamic behavior of turbine foundation considering full interaction among facility, structure and soil Fang Ming Scholl of Civil Engineering, Harbin Institute of Technology, China Wang Tao Institute of

More information

Chapter 3 LAMINATED MODEL DERIVATION

Chapter 3 LAMINATED MODEL DERIVATION 17 Chapter 3 LAMINATED MODEL DERIVATION 3.1 Fundamental Poisson Equation The simplest version of the frictionless laminated model was originally introduced in 1961 by Salamon, and more recently explored

More information

CRACK FORMATION AND CRACK PROPAGATION INTO THE COMPRESSION ZONE ON REINFORCED CONCRETE BEAM STRUCTURES

CRACK FORMATION AND CRACK PROPAGATION INTO THE COMPRESSION ZONE ON REINFORCED CONCRETE BEAM STRUCTURES S. Kakay et al. Int. J. Comp. Meth. and Exp. Meas. Vol. 5 No. (017) 116 14 CRACK FORMATION AND CRACK PROPAGATION INTO THE COMPRESSION ZONE ON REINFORCED CONCRETE BEAM STRUCTURES SAMDAR KAKAY DANIEL BÅRDSEN

More information

DISPENSA FEM in MSC. Nastran

DISPENSA 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 information

A BEM implicit formulation applied to damage mechanics models

A BEM implicit formulation applied to damage mechanics models A BEM implicit formulation applied to damage mechanics models A. S. ~otta', W. S. ~enturini' & A. ~enallal~ 'S& Carlos School oj'engineering, University of $60 Paulo, Brazil 2 Laboratoire de Me'canique

More information

16 Singular perturbations

16 Singular perturbations 18.354J Nonlinear Dynamics II: Continuum Systems Lecture 1 6 Spring 2015 16 Singular perturbations The singular perturbation is the bogeyman of applied mathematics. The fundamental problem is to ask: when

More information

3D Elasticity Theory

3D Elasticity Theory 3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.

More information

EXTENDED 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 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 information

Load Cell Design Using COMSOL Multiphysics

Load Cell Design Using COMSOL Multiphysics Load Cell Design Using COMSOL Multiphysics Andrei Marchidan, Tarah N. Sullivan and Joseph L. Palladino Department of Engineering, Trinity College, Hartford, CT 06106, USA joseph.palladino@trincoll.edu

More information

ΙApostolos Konstantinidis Diaphragmatic behaviour. Volume B

ΙApostolos Konstantinidis Diaphragmatic behaviour. Volume B Volume B 3.1.4 Diaphragmatic behaviour In general, when there is eccentric loading at a floor, e.g. imposed by the horizontal seismic action, the in-plane rigidity of the slab forces all the in-plane points

More information

Piezoelectric Bimorph Response with Imperfect Bonding Conditions

Piezoelectric 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 information

International Journal of Advanced Engineering Technology E-ISSN

International Journal of Advanced Engineering Technology E-ISSN Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,

More information

UNIT II SHALLOW FOUNDATION

UNIT II SHALLOW FOUNDATION Introduction UNIT II SHALLOW FOUNDATION A foundation is a integral part of the structure which transfer the load of the superstructure to the soil. A foundation is that member which provides support for

More information

ASSESSMENT OF DYNAMICALLY LOADED CRACKS IN FILLETS

ASSESSMENT OF DYNAMICALLY LOADED CRACKS IN FILLETS ASSESSMENT OF DNAMICALL LOADED CRACKS IN FILLETS Uwe Zencker, Linan Qiao, Bernhard Droste Federal Institute for Materials Research and Testing (BAM) 12200 Berlin, Germany e-mail: zencker@web.de Abstract

More information

The Effect of Distribution for a Moving Force

The Effect of Distribution for a Moving Force Paper Number 66, Proceedings of ACOUSTICS 2011 2-4 November 2011, Gold Coast, Australia The Effect of Distribution for a Moving Force Ahmed M. Reda (1,2), Gareth L. Forbes (2) (1) Atkins, Perth, Australia

More information

Back Analysis of Measured Displacements of Tunnels

Back 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 information

BENCHMARK LINEAR FINITE ELEMENT ANALYSIS OF LATERALLY LOADED SINGLE PILE USING OPENSEES & COMPARISON WITH ANALYTICAL SOLUTION

BENCHMARK LINEAR FINITE ELEMENT ANALYSIS OF LATERALLY LOADED SINGLE PILE USING OPENSEES & COMPARISON WITH ANALYTICAL SOLUTION BENCHMARK LINEAR FINITE ELEMENT ANALYSIS OF LATERALLY LOADED SINGLE PILE USING OPENSEES & COMPARISON WITH ANALYTICAL SOLUTION Ahmed Elgamal and Jinchi Lu October 07 Introduction In this study: I) The response

More information

Transactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN

Transactions on Engineering Sciences vol 14, 1997 WIT Press,  ISSN On the Computation of Elastic Elastic Rolling Contact using Adaptive Finite Element Techniques B. Zastrau^, U. Nackenhorst*,J. Jarewski^ ^Institute of Mechanics and Informatics, Technical University Dresden,

More information

Lecture 11: The Stiffness Method. Introduction

Lecture 11: The Stiffness Method. Introduction Introduction Although the mathematical formulation of the flexibility and stiffness methods are similar, the physical concepts involved are different. We found that in the flexibility method, the unknowns

More information

MATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f

MATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS 2 2 2 f m f f m T m f m f f m v v + + = + PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown

More information

MMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM

MMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,

More information

MTE 119 STATICS FINAL HELP SESSION REVIEW PROBLEMS PAGE 1 9 NAME & ID DATE. Example Problem P.1

MTE 119 STATICS FINAL HELP SESSION REVIEW PROBLEMS PAGE 1 9 NAME & ID DATE. Example Problem P.1 MTE STATICS Example Problem P. Beer & Johnston, 004 by Mc Graw-Hill Companies, Inc. The structure shown consists of a beam of rectangular cross section (4in width, 8in height. (a Draw the shear and bending

More information

MASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK

MASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK MASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK J. Pina-Henriques and Paulo B. Lourenço School of Engineering, University of Minho, Guimarães, Portugal Abstract Several continuous and discontinuous

More information

STATICALLY INDETERMINATE STRUCTURES

STATICALLY INDETERMINATE STRUCTURES STATICALLY INDETERMINATE STRUCTURES INTRODUCTION Generally the trusses are supported on (i) a hinged support and (ii) a roller support. The reaction components of a hinged support are two (in horizontal

More information

Analysis of Open Circular Reinforced Concrete Shell Roof with Curved Edge Simply Supported

Analysis of Open Circular Reinforced Concrete Shell Roof with Curved Edge Simply Supported Proceedings of IOE Graduate Conference, 2017 Volume: 5 ISSN: 2350-8914 (Online), 2350-8906 (Print) Analysis of Open Circular Reinforced Concrete Shell Roof with Curved Edge Simply Supported Abhishek Adhikari

More information

ABHELSINKI UNIVERSITY OF TECHNOLOGY

ABHELSINKI UNIVERSITY OF TECHNOLOGY ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural

More information

Generalized Plane Piezoelectric Problem: Application to Heterostructure Nanowires

Generalized Plane Piezoelectric Problem: Application to Heterostructure Nanowires Generalized Plane Piezoelectric Problem: Application to Heterostructure Nanowires H. T. Mengistu and A. Garcia-Cristobal Institute of Material Science, University of Valencia, P.O.Box 22085, E-46071, Valencia,

More information

Boundary element analysis of FRP-concrete delamination

Boundary element analysis of FRP-concrete delamination Boundary element analysis of FRP-concrete delamination F. Freddi 1, A. Salvadori 2 &M.Savoia 3 1 Department of Civil Engineering., University of Parma, Italy 2 Department of Civil Engineering., University

More information

Aim of the study Experimental determination of mechanical parameters Local buckling (wrinkling) Failure maps Optimization of sandwich panels

Aim of the study Experimental determination of mechanical parameters Local buckling (wrinkling) Failure maps Optimization of sandwich panels METNET Workshop October 11-12, 2009, Poznań, Poland Experimental and numerical analysis of sandwich metal panels Zbigniew Pozorski, Monika Chuda-Kowalska, Robert Studziński, Andrzej Garstecki Poznan University

More information

THE GENERAL ELASTICITY PROBLEM IN SOLIDS

THE GENERAL ELASTICITY PROBLEM IN SOLIDS Chapter 10 TH GNRAL LASTICITY PROBLM IN SOLIDS In Chapters 3-5 and 8-9, we have developed equilibrium, kinematic and constitutive equations for a general three-dimensional elastic deformable solid bod.

More information

Transactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X

Transactions 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 information

elastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK

elastoplastic 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 information

Analytical formulation of Modified Upper Bound theorem

Analytical formulation of Modified Upper Bound theorem CHAPTER 3 Analytical formulation of Modified Upper Bound theorem 3.1 Introduction In the mathematical theory of elasticity, the principles of minimum potential energy and minimum complimentary energy are

More information

General Solutions for Axisymmetric Elasticity Problems of Elastic Half Space using Hankel Transform Method

General Solutions for Axisymmetric Elasticity Problems of Elastic Half Space using Hankel Transform Method ISSN (Print) : 19-861 ISSN (Online) : 975-44 General Solutions for Axisymmetric Elasticity Problems of Elastic Half Space using Hankel Transform Method Charles Chinwuba Ike Dept of Civil Engineering Enugu

More information

PILE SOIL INTERACTION MOMENT AREA METHOD

PILE SOIL INTERACTION MOMENT AREA METHOD Pile IGC Soil 2009, Interaction Moment Guntur, INDIA Area Method PILE SOIL INTERACTION MOMENT AREA METHOD D.M. Dewaikar Professor, Department of Civil Engineering, IIT Bombay, Mumbai 400 076, India. E-mail:

More information

Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams

Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams V. Lotfi Department of Civil and Environmental Engineering, Amirkabir University, Iran Abstract A new approach is

More information

Table of Contents. Preface... 13

Table of Contents. Preface... 13 Table of Contents Preface... 13 Chapter 1. Vibrations of Continuous Elastic Solid Media... 17 1.1. Objective of the chapter... 17 1.2. Equations of motion and boundary conditions of continuous media...

More information

Lecture 12: Finite Elements

Lecture 12: Finite Elements Materials Science & Metallurgy Part III Course M6 Computation of Phase Diagrams H. K. D. H. Bhadeshia Lecture 2: Finite Elements In finite element analysis, functions of continuous quantities such as temperature

More information

Modelling of non linear soil-structure interface behaviour using M. Boulon & P. G arnica

Modelling of non linear soil-structure interface behaviour using M. Boulon & P. G arnica Modelling of non linear soil-structure interface behaviour using M. Boulon & P. G arnica BEM Cecfecc P, France Abstract The soil-structure interaction problem is considered as a nonlinear contact problem

More information

Transactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X

Transactions on Modelling and Simulation vol 12, 1996 WIT Press,   ISSN X A Boundary Element Method for analysis of bone remodelling M. Martinez, M.H. Aliabadi, H. Power PF&M&x TW/fwfe qf T^cAWogy, As/mr Ashurst, Southampton SO40 7AA, UK Abstract A boundary element formulation

More information

Soil Dynamics Prof. Deepankar Choudhury Department of Civil Engineering Indian Institute of Technology, Bombay

Soil Dynamics Prof. Deepankar Choudhury Department of Civil Engineering Indian Institute of Technology, Bombay Soil Dynamics Prof. Deepankar Choudhury Department of Civil Engineering Indian Institute of Technology, Bombay Module - 5 Machine Foundations Lecture - 31 EHS Theory, Vibrational Control Let us start our

More information

Discrete Element Modelling of a Reinforced Concrete Structure

Discrete Element Modelling of a Reinforced Concrete Structure Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr

More information

Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering

Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering

More information

Benefits of Collaboration between Centrifuge Modeling and Numerical Modeling. Xiangwu Zeng Case Western Reserve University, Cleveland, Ohio

Benefits of Collaboration between Centrifuge Modeling and Numerical Modeling. Xiangwu Zeng Case Western Reserve University, Cleveland, Ohio Benefits of Collaboration between Centrifuge Modeling and Numerical Modeling Xiangwu Zeng Case Western Reserve University, Cleveland, Ohio ABSTRACT There is little doubt that collaboration between centrifuge

More information

Flin Flon Mining Belt

Flin Flon Mining Belt EOSC433: Geotechnical Engineering Practice & Design Lecture 7: Stress Analysis around Underground Openings 1 of 40 Erik Eberhardt UBC Geological Engineering EOSC 433 (2007) Flin Flon Mining Belt Since

More information

3D dynamic crack propagation analysis with PDS-FEM

3D dynamic crack propagation analysis with PDS-FEM JAMSTEC-R IFREE Special Issue, November 2009 3D dynamic crack propagation analysis with PDS-FEM 1*, Hide Sakaguchi 1, Kenji Oguni 2, Muneo Hori 2 discretizations. Keywords +81-45-778-5972 lalith@jamstec.go.jp

More information

Understand basic stress-strain response of engineering materials.

Understand basic stress-strain response of engineering materials. Module 3 Constitutive quations Learning Objectives Understand basic stress-strain response of engineering materials. Quantify the linear elastic stress-strain response in terms of tensorial quantities

More information