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

Size: px
Start display at page:

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

Transcription

1 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 work is to present an elastodynamic coupling of a plate in bending with the three-dimensional soil medium by using the Boundary Element Method formulation. This scheme is adopted in order to better represent some aspects of this kind of structural system, as for instance, the stress concentration. The formulation used to treat the soil media is based on time domain fundamental solution, while the dynamic behavior of the plate in bending is represented making use of the mass matrix approach. 1 Introduction The coupling of any structure with the soil medium is usually idealized by using the Finite Element Method (FEM) and the Boundary Element Method (BEM) formulations, respectively. Thus, one can take the advantages that each method presents to treat each particular sub-domain. For elastostatic cases, the first work reporting such procedure is due to ZIENKIEWICZ et alli.^. Since then, many other works have been presented in the literature improving the technique. However, the finite element formulation presents difficulties to exhibit some particular features of this coupling, for instance stress concentrations along the contact area contour. On the other hand, when one applies the boundary element technique, stress concentrations can be easily represented. In this way, PAIVA & BUTTERFIELD^ have proposed to couple both media, plate and soil, by using only elastostatic BEM formulations. In their work compatibility is enforced only for vertical displacements. An important feature close related to the stress concentration presence is the numerical instability when one is carrying out elastodynamics analysis. For this kind of study, the finite element technique presents less problems because this formulation reduces the stress concentration. On the other hand, when BEM is adopted, its formulation always leads to higher numerical instabilities.

2 364 Boundary Elements The aim of this work is to propose an elastostatics and elastodynamics coupling for plates in bending and the soil medium, where both sub-domains are treated by BEM. Using this combination one can check out the presence of stress concentrations and their numerical influence on the elastodynamics behavior. For this scheme, the plate is assumed rigid in its plane. The coupling is enforced by using the sub-region technique; the degrees of freedom of the threedimensional (soil) and the bidimensional (plate) media are appropriately combined. An elastostatic example is shown to illustrate the quality of the results obtained by using the proposed approach, comparing them with other numerical solutions. In addition, two elastodynamics problems are analyzed to point out the influence of the plate/soil stiffness ratio on the stress concentration and numerical stability. 2 Three-dimensional formulation For any isotropic body exhibiting a linear elastic behavior, the differential equation which governs its equilibrium, the so-called Navier-Cauchy equation, is given by (C,' -CJ)^,, +C^,^. +6, / p = 6, (for,,; = 1,2,3) (1) where bf represents the body forces, p is the media density, d and 2 are longitudinal and shear wave propagation velocities, respectively. In elastodynamics, it is usual to take a particular case of the Stokes' state [STOKES**] as fundamental solution. The usual way to find the desired fundamental solution is changing the body force term, in eqn(l), by a time dependent concentrated load given by where /(T) is a time-dependent load function, 6(s,q) is the Dirac's delta function and Ski the Kronecker's delta. The time-dependent load function adopted here is the one proposed by CODA & VENTURING, obtained by assuming, (2) /(r)=(#(r)-#(r+af))/af (3) By integrating eqn(l), after replacing the f(i) defined in eqn(3), one obtains the so-called Stokes' fundamental solution, as follows, where,,o/a/ (4)

3 Boundary Elements 365 4(,.,;,,0) = ^ (^-^- (r'-^ur-- I L ^U.^--l4^--lK^-4 (5) The GrafFi's reciprocity theorem for dynamic problems can be written for any two elastodynamics states, here represented by their displacement field u and w', in the case of quiescent past, as follows, %X9,4^ in which g is used to denote thefieldpoint taken along the boundary. In order to write the basic integral representation of displacements one has only to replace the body force term as in eqn(2) and the function /(T) by the particular time loading given in eqn(3). After replacing these values one has,,f;^^ (7) where Cki(Q,s) is the classical free term of the displacement integral representation that can be easily computed as in HARTMANN^. In order to achieve the final algebraic representations, tractions and displacements have to be approximated by, ^(<2V) = ^V'>/^ (8) in which j stands for the space elements, <t^ is the space approximation function, while v[/ (\ / = 1 is taken) gives the approximation along time. Flat quadratic elements with eight nodes have been adopted for the space discretization. Integrals over singular elements have been computed using a special algorithm where the finite-part concept, the Kutt's quadrature, KUTT^, was adopted for a small region around the load point, while the Gauss quadrature was taken to compute the remaining part of the element together with an appropriate element subdivision, CODA & VENTURING. For non-

4 366 Boundary Elements with an appropriate element subdivision, CODA & VENTURING For nonsingular integrals, the Gauss quadrature was also adopted with an element subdivision technique. As usual for boundary element formulations, the displacement integral representation, eqn (7), is transformed into an algebraic equation after choosing an appropriate number of load points and carrying out all the space and time integrals. So the classical BEM system of linear equations is achieved, as follows, H»U* = G*P* (9) where 9 ranges from 1 to the integration step number Nf. For any analyzed instant "t" the expression (9) can be written as: HU = GP+A (10) where "A" contains all the history of the movement represented in eqn (9). Eqn (10) has been presented in the standard form appropriate to implement the subregion technique. 3 Plate formulation For a linear elastic thin plate, the dynamic equilibrium is governed by the following differential equation, Z).AAw+ pw - g (11) where w is the transversal displacement, g the transversal load applied on the domain, p the plate density and D the plate rigidity. The fundamental solution employed for both, the elastostatics and elastodynamics cases, is the same one, due to DANSON^. By substituting the domain load g by the Dirac's delta function in eqn (11), in its elastostatic form, written in polar coordinates, one obtains, where r is the distance between the load point q and the field point/?. Applying the Betti's reciprocity theorem, and assuming that the variables employed in the plate formulation are: w and dw/ch (displacement derivative with respect to the outward normal direction) to define the displacement vector, and Vn (shear force) and Mn (bending moment) to define the effort vector, one can write the plate bending displacement integral representation as follows,

5 Boundary Elements 367! h (13) where ^c/ is the corner reaction and g* is the fundamental load. As it has been done for the soil formulation, in order to solve a plate problem by BEM it is necessary to transform eqn (13) into its algebraic representation. In this case, the plate boundary must to be discretized. For that quadratic approximation elements with linear geometry were adopted. As eqn (13) contains a domain integral with unknown density, domain discretization is needed as well. In this case, a cell similar to the surface element adopted for the 3D formulation was taken, generating therefore the mass matrix M The scheme adopted here to achieve algebraic representations avoids the slope integral equation usually employed for plate bending analysis with BEM Only non singular representation have been selected. Thus, for each node we have taken two collocation points, consequently achieving equal numbers of equations and unknowns. The clear advantage of employing this scheme is avoiding singular integration and keeping representation of the same nature. Schemes similar to that have already proved to give more precise numerical results, as can be seen in VENTURINI & PAIVA^. For the scheme adopted here the collocations are taken along the outward normal boundary direction as showed in figure 1. boundary load points inside points' corner /load point Figure 1. Plate geometry, discretization and collocation points Thus, using the proposed scheme one writes the displacement equation, eqn(13), only for outside and inside collocation points. In the first case, the domain integral containing the fundamental load, the third term of eqn (13) left hand side, becomes zero, while for the second case, it becomes equal to w(q). As there is no singular integral, the algebraic system can be numerically

6 368 Boundary Elements achieved by applying only the Gauss's quadrature with or without element subdivision, according to the required accuracy. Other important aspects to be emphasized are the treatment given to the corners reactions Rci and to the contact reaction and acceleration domain integrals. In order to take into account the corner reactions, one extra displacement equation is written for each node exhibiting geometric discontinuities. For those extra equations, the collocations are taken out of domain as well, following the corner's bisector direction, as shown in figure 1. This procedure generates one more integral equation for each corner. Thus, one has included more equations and more unknowns, corner displacements or reactions, but the final system of equation remains square. Regarding the domain integrals, they are computed numerically, dividing the domain into cells, over which the Gauss's quadrature is applied. These cells are nothing more than flat surface element; in particular, the cell here adopted corresponds to the flat quadrilateral with 8 nodes used to discretize the boundary of the soil domain. After performing the integrals over boundary elements and internal cells, and taking into account the corner values, one can write the global system due to the dynamic equilibrium, making use of the Newmark Method, as follows: (\-20)GP +/3GP + ~ ~ s ~ ~ S-\ J * H\U (14) ~ J ~ _,_, where H is the classical BEM matrix of the plate problem, M is the mass matrix, G is the matrix that contains the subgrade reaction coefficients, P gives the subgrade reactions, U is the displacement vector, At represents the time step and P the integration method parameter. Finally eqn (14) can be easily written, for any instant, in the form given in eqn (9). 4 Plate-soil interaction In this work, the sub-region technique was chosen to establish the coupling between the two media. In order to perform this interaction, one can connect the two regions making use of the equilibrium and displacement compatibility conditions, as follows, lj(v) _ fj(jf) and p(v) + p(j>) =o (15) Finally, one obtains the coupled matricial equation given by, system, giving origin to a single

7 Boundary Elements 369 //<»> -G<» 0 [0 0 0 (16) where the bar indicates prescribed loads on the interface, e gives the region not belonging to the coupling interface. The first superscript stands for the region of the corresponding interface values and the second indicates the adjacent subregion. The vector F contains all the prescribed value influences. 5 Numerical results For all examples presented in this paper we have adopted the same discretization characterized by: a) for BEM discretization we have taken 16 cells and 16 boundary elements as shown in figure 1; b) for FEM discretization, 96 HCT/CSTfiniteelements. In the first place an elastostatic example with a uniformly distributed load is analyzed and the behavior is compared with the PAIVA & BUTTERFIELD^ and CODA & VENTURING responses. This comparison is made to point out the formulation accuracy. To run this example the following plate geometry was adopted: side length equal to 12m, thickness equal to 0.1/m and the elastic parameters used are: ES = 0.26 x 10^ kn/m^ (soil elastic modulus), Ep = 9.78 x 10^ KN/m^ (plate elastic modulus), vp = vs = 0.3 (Poisson's coefficients). 3. I " " 2 2- PAIVA» CODA * BARRETTO Figure 2. Displacements (w) and subgrade reactions (Sr) along the axis A-B In next example the plate has side length equal to 10m, thickness equal to 0.5m and is subjected to a uniformly distributed load equal to 300 Kgf7nf. First of all, the elastostatic behavior along the direction A-B is shown to emphasize the marked differences in terms of stress concentration near the boundary. After that the elastodynamic analysis is carried out and the transversal displacement at the center of the plate, as well as the subgrade

8 370 Boundary Elements reaction at the middle of sides, are shown to illustrate the expected differences. For this elastodynamic analysis the time step used was equal to s, while the elastic parameters adopted are given by: ES = 7.8 x 10" Kgf/m/s^ ; Ep = 7.8 x 10*0 Kgf/m/s^; vp = vs = 0.3 and the material densities by: p^ = 2500 Kg/nf; p,= 1800Kg/m*.! & 5 ^, A BARRETTO] - CODA I Figure 3. Displacements (w) and subgrade reactions (Sr) along the axis A-B! I -A BARRETTO - CODA,.. -BARRETTO -CODA Figure 4. Displacement (w) at the center Figure 5. Subgrade reaction (Sr) at the middle of the side 6 Conclusions The coupling of boundary element formulations for plate bending and threedimensional solid models has been shown for elastostatics and elastodynamics analyses. The numerical results achieved by using the implemented algorithm were compared with equivalent coupling where finite elements were adopted to represent the plate in bending. The static values shown that the BEM/BEM coupling can capture with accuracy the stress concentration along the boundary, while for BEM/FEM coupling the results are reduced over this region. However, the accuracy exhibited to model stress peak produces instabilities for BEM/BEM elastodynamics analysis. Working with damped values the BEM/FEM coupling is more appropriate to analyze this kind of problem.

9 Boundary Elements References 1. Barretto, S. F. A. Elastostatic and elastodynamic analysis of plates, using the Boundary Elements Method, with soil-structure interaction. Ph. D. Thesis, University of Sao Paulo, Sao Carlos, Brebbia, C. A. The boundary element method for engineers. London Pentech Press, Coda, H. B. ; Venturini, W. S. Alternative boundary element formulation for elastodynamics. In: International Conference on Boundary Elements in Engineering, 12., Sapporo, Japan, Sept., Proceedings. Southampton, CML; Berlin, Springer-Verlag, v. 1, p , Coda, H. B. Three-dimensional transient analysis of structures by a BEM- FEM combination. Ph. D. Thesis, University of Sao Paulo Sao Carlos 5. Coda, H. B. ; Venturini, W. S. Three-dimensional transient BEM analysis. Computers & Structures, v. 56, n. 5, p , Coda, H. B ; Venturini, W S. Numerical evaluation of flat singular boundary elements in elastostatics and elastodynamics, Boundary Elements communications, v. 6, p , Coda, H B ; Venturini, W. S. Non-singular time-stepping BEM for transient elastodynamic analysis. Engineering Analysis with Boundary Elements, v. 15, p , Coda, H.B. ; Venturini, W.S. Time domain BEM/FEM approach applied to elastodynamic analysis. Conference of Italian Group of Computational Mechanics, Padova, 1996 (to appear). 9. Danson, D. J Analysis of plate bending problems by direct boundary element method. Southampton. M. Sc. Dissertation - University of Southampton, lo.hartmann, F. Computing the C-matrix on non-smooth boundary points. In: Brebbia, C. A, ed. New developments in boundary element methods Southampton, CMP Publ, ll.kutt, H. R. WISK 178: quadrature formulae for finite-part integrals. Pretoria, National Research Institute for Mathematical Sciences. CSIR Special Report, Paiva, J. B. ; Butterfield, R. Numerical analysis of plate-soil interaction. In: Papadrakakis, M. et Topping, B H V., eds. Advances in computational mechanics. (Proc. 2nd International Conference on Computational Structures Technology, Athens, Greece, 30th Aug. - 1st Sept.). Edinburgh, Civil-Comp Press, p , Patterson, C. ; Sheik, M. A. ; ScholGeld, R. P. On fae a^acof/om q/vac indirect discrete method for three dimensional design problems. In: 1st. Boundary Element Technology Conference, Adelaide, Australia. Proceedings. p , Stokes, G. G On the dynamical theory of diffraction. Transactions of the Cambridge Philosophical Society, v. 9, p. 1, 1849.

10 372 Boundary Elements 15.Venturing W.S. ; Paiva, J.B. Boundary element for plate bending analysis. Engineering Analysis with Boundary Elements, v. 11, n. 1, p. 1-8, Warburton, G. B 77?2 dynamical behavior of structures. New York, Pergamon Press, IT.Zienkiewicz, O.C. ; Kelly, D.W. ; Bettess, P. The coupling of the finite element method and boundary solution procedures. Int. J. Num. Meth. Eng., v.ll, p , 1977.

Transactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X

Transactions on Modelling and Simulation vol 18, 1997 WIT Press,   ISSN X A simple coupling of 2D BEM and FEM bar model applied to mass matrix elastodynamic analysis Coda,H.B*; Venturini W.S* and Aliabadi, M.HT * Sao Carlos School of Eng., University of Sao Paulo, Br, venturin@usp.

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

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

Abstract. Introduction

Abstract. Introduction Structural vibration isolation by rows of piles S.E. Kattis", D. Polyzos*, D.E. Beskos* "Department of Mechanical Engineering, University ofpatras, G^-26 j 00 f6zrr^ Greece ^Department of Civil Engineering,

More information

An alternative multi-region BEM technique for layered soil problems

An alternative multi-region BEM technique for layered soil problems An alternative multi-region BM technique for layered soil problems D.B. Ribeiro & J.B. Paiva Structural ngineering Department, São Carlos ngineering School, University of São Paulo, Brazil. Abstract Different

More 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

Appendix 1. Numerical Integration: Regular Integrals

Appendix 1. Numerical Integration: Regular Integrals The Dual Reciprocity 269 Appendix 1 Numerical Integration: Regular Integrals 1 Introduction The numerical integration formulae given in this Appendix are used for integration over boundary elements which

More 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

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

The Finite Element Method for Solid and Structural Mechanics

The Finite Element Method for Solid and Structural Mechanics The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in

More information

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 3, No 4, 2013

INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 3, No 4, 2013 INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 3, No 4, 2013 Copyright by the authors - Licensee IPA- Under Creative Commons license 3.0 Research article ISSN 0976 4399 Pure bending analysis

More information

A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA

A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA 8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC

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

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

MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING

MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING 144 MIXED RECTANGULAR FINITE ELEMENTS FOR PLATE BENDING J. N. Reddy* and Chen-Shyh-Tsay School of Aerospace, Mechanical and Nuclear Engineering, University of Oklahoma, Norman, Oklahoma The paper describes

More information

AN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS

AN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS th Pan-American Congress of Applied Mechanics January 04-08, 00, Foz do Iguaçu, PR, Brazil AN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS Otávio Augusto Alves

More 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

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

SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION

SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION 8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION Theodore V. Gortsas

More information

COPYRIGHTED MATERIAL. Index

COPYRIGHTED MATERIAL. Index Index A Admissible function, 163 Amplification factor, 36 Amplitude, 1, 22 Amplitude-modulated carrier, 630 Amplitude ratio, 36 Antinodes, 612 Approximate analytical methods, 647 Assumed modes method,

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

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

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

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

Transactions on Modelling and Simulation vol 2, 1993 WIT Press, ISSN X

Transactions on Modelling and Simulation vol 2, 1993 WIT Press,  ISSN X A general integral-equation formulation for anisotropic elastostatics M.M. Perez"'*, L.C. Wrobel* * Wessex Institute of Technology, University of Portsmouth, Ashurst Lodge, Ashurst, Southampton S04 2AA,

More information

Chapter 2 Finite Element Formulations

Chapter 2 Finite Element Formulations Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are

More 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

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

JEPPIAAR ENGINEERING COLLEGE

JEPPIAAR ENGINEERING COLLEGE JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III

More information

Vibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1

Vibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1 APCOM & ISCM -4 th December, 23, Singapore Vibration of Thin Beams by PIM and RPIM methods *B. Kanber¹, and O. M. Tufik Mechanical Engineering Department, University of Gaziantep, Turkey. *Corresponding

More information

Bending of Simply Supported Isotropic and Composite Laminate Plates

Bending of Simply Supported Isotropic and Composite Laminate Plates Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,

More information

Structural Acoustics Applications of the BEM and the FEM

Structural Acoustics Applications of the BEM and the FEM Structural Acoustics Applications of the BEM and the FEM A. F. Seybert, T. W. Wu and W. L. Li Department of Mechanical Engineering, University of Kentucky Lexington, KY 40506-0046 U.S.A. SUMMARY In this

More information

Transactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X

Transactions on Modelling and Simulation vol 18, 1997 WIT Press,   ISSN X A direct formulation for three-dimensional elastoplastic boundary elements A. P. Cisilino*,**, M. H. Aliabadi* and J. L. Otegui** *Wessex Institute of Technology, As hurst Lodge, As hurst, Southampton

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

Static & Dynamic. Analysis of Structures. Edward L.Wilson. University of California, Berkeley. Fourth Edition. Professor Emeritus of Civil Engineering

Static & Dynamic. Analysis of Structures. Edward L.Wilson. University of California, Berkeley. Fourth Edition. Professor Emeritus of Civil Engineering Static & Dynamic Analysis of Structures A Physical Approach With Emphasis on Earthquake Engineering Edward LWilson Professor Emeritus of Civil Engineering University of California, Berkeley Fourth Edition

More information

HANDBUCH DER PHYSIK HERAUSGEGEBEN VON S. FLÜGGE. BAND VIa/2 FESTKÖRPERMECHANIK II BANDHERAUSGEBER C.TRUESDELL MIT 25 FIGUREN

HANDBUCH DER PHYSIK HERAUSGEGEBEN VON S. FLÜGGE. BAND VIa/2 FESTKÖRPERMECHANIK II BANDHERAUSGEBER C.TRUESDELL MIT 25 FIGUREN HANDBUCH DER PHYSIK HERAUSGEGEBEN VON S. FLÜGGE BAND VIa/2 FESTKÖRPERMECHANIK II BANDHERAUSGEBER C.TRUESDELL MIT 25 FIGUREN SPRINGER-VERLAG BERLIN HEIDELBERG NEWYORK 1972 Contents. The Linear Theory of

More information

Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method

Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method 9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached

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

A *69>H>N6 #DJGC6A DG C<>C::G>C<,8>:C8:H /DA 'D 2:6G - ( - ) +"' ( + -"( (' (& -+" % '('%"' +"-2 ( -!"',- % )% -.C>K:GH>IN D; AF69>HH>6,-+

A *69>H>N6 #DJGC6A DG C<>C::G>C<,8>:C8:H /DA 'D 2:6G - ( - ) +' ( + -( (' (& -+ % '('%' +-2 ( -!',- % )% -.C>K:GH>IN D; AF69>HH>6,-+ The primary objective is to determine whether the structural efficiency of plates can be improved with variable thickness The large displacement analysis of steel plate with variable thickness at direction

More information

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi

More information

First- and second-order sensitivity analysis schemes by collocation-type Trefftz method

First- and second-order sensitivity analysis schemes by collocation-type Trefftz method Computer Assisted Mechanics and Engineering Sciences, 4: 477 489, 1997 Copyright c 1997 by Polska Akademia Nauk First- and second-order sensitivity analysis schemes by collocation-type Trefftz method Eisuke

More information

INFLUENCE OF TYPE OF WAVE AND ANGLE OF INCIDENCE ON THE SEISMIC RESPONSE OF PILE FOUNDATIONS AND PILE SUPPORTED STRUCTURES

INFLUENCE OF TYPE OF WAVE AND ANGLE OF INCIDENCE ON THE SEISMIC RESPONSE OF PILE FOUNDATIONS AND PILE SUPPORTED STRUCTURES COMPDYN ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V.Plevris (eds.) Corfu, Greece, -8 May INFLUENCE OF TYPE

More information

BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I

BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T - I INTRODUCTION

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 Analysis of complexfluid-structureinteraction problems using a triangular boundary element for the fluid dynamics A. Hopf, G. Hailfmger, R. Krieg Forschungszentrum Karlsruhe GmbH, Institutfur D- 7^027

More 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

Schur decomposition in the scaled boundary finite element method in elastostatics

Schur decomposition in the scaled boundary finite element method in elastostatics IOP Conference Series: Materials Science and Engineering Schur decomposition in the scaled boundary finite element method in elastostatics o cite this article: M Li et al 010 IOP Conf. Ser.: Mater. Sci.

More information

DYNAMIC RESPONSE OF THIN-WALLED GIRDERS SUBJECTED TO COMBINED LOAD

DYNAMIC RESPONSE OF THIN-WALLED GIRDERS SUBJECTED TO COMBINED LOAD DYNAMIC RESPONSE OF THIN-WALLED GIRDERS SUBJECTED TO COMBINED LOAD P. WŁUKA, M. URBANIAK, T. KUBIAK Department of Strength of Materials, Lodz University of Technology, Stefanowskiego 1/15, 90-924 Łódź,

More information

LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS

LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS 5 th GRACM International Congress on Computational Mechanics Limassol, 29 June 1 July, 2005 LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS G.A. Drosopoulos I, G.E.

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

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

202 Index. failure, 26 field equation, 122 force, 1

202 Index. failure, 26 field equation, 122 force, 1 Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic

More information

Back Matter Index The McGraw Hill Companies, 2004

Back Matter Index The McGraw Hill Companies, 2004 INDEX A Absolute viscosity, 294 Active zone, 468 Adjoint, 452 Admissible functions, 132 Air, 294 ALGOR, 12 Amplitude, 389, 391 Amplitude ratio, 396 ANSYS, 12 Applications fluid mechanics, 293 326. See

More information

Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation

Quintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation

More information

Solving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells

Solving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells Copyright 2010 Tech Science Press CMES, vol.57, no.2, pp.109-135, 2010 Solving Elastic Problems with Local Boundary Integral Equations (LBIE) and Radial Basis Functions (RBF) Cells E. J. Sellountos 1,

More information

DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD

DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD R. A. Amaro

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

Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis

Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan

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

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

AN INTEGRATED KIRCHHOFF PLATE ELEMENT BY GALERKIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION

AN INTEGRATED KIRCHHOFF PLATE ELEMENT BY GALERKIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION AN INTEGRATED IRCHHOFF PLATE ELEMENT BY GALERIN METHOD FOR THE ANALYSIS OF PLATES ON ELASTIC FOUNDATION Ragesh.P.P, V.Mustafa 2, T.P.Somasundaran 3 (Research Scholar, National Institute of Technology Calicut,

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

Chapter 3. Load and Stress Analysis

Chapter 3. Load and Stress Analysis Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3

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 Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation

More 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

TIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL

TIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL TIME-DEPENDENT BEHAVIOR OF PILE UNDER LATERAL LOAD USING THE BOUNDING SURFACE MODEL Qassun S. Mohammed Shafiqu and Maarib M. Ahmed Al-Sammaraey Department of Civil Engineering, Nahrain University, Iraq

More information

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.

D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. GTE 2016 Q. 1 Q. 9 carry one mark each. D : SOLID MECHNICS Q.1 single degree of freedom vibrating system has mass of 5 kg, stiffness of 500 N/m and damping coefficient of 100 N-s/m. To make the system

More information

Sound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur

Sound Propagation through Media. Nachiketa Tiwari Indian Institute of Technology Kanpur Sound Propagation through Media Nachiketa Tiwari Indian Institute of Technology Kanpur LECTURE-13 WAVE PROPAGATION IN SOLIDS Longitudinal Vibrations In Thin Plates Unlike 3-D solids, thin plates have surfaces

More information

PREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO

PREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG

More information

COUPLED FINITE-INFINITE ELEMENTS MODELING OF BUILDING FRAME-SOIL INTERACTION SYSTEM

COUPLED FINITE-INFINITE ELEMENTS MODELING OF BUILDING FRAME-SOIL INTERACTION SYSTEM VOL. 4, NO. 10, DECEMBER 009 SSN 1819-6608 006-009 Asian Research Publishing Network (ARPN. All rights reserved. COUPLED FNTE-NFNTE ELEMENTS MODELNG OF BULDNG FRAME-SOL NTERACTON SYSTEM Ramakant Agrawal

More information

Lecture 15 Strain and stress in beams

Lecture 15 Strain and stress in beams Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME

More information

Iraq Ref. & Air. Cond. Dept/ Technical College / Kirkuk

Iraq Ref. & Air. Cond. Dept/ Technical College / Kirkuk International Journal of Scientific & Engineering Research, Volume 6, Issue 4, April-015 1678 Study the Increasing of the Cantilever Plate Stiffness by Using s Jawdat Ali Yakoob Iesam Jondi Hasan Ass.

More information

THE TRACTION BOUNDARY CONTOUR METHOD FOR LINEAR ELASTICITY

THE TRACTION BOUNDARY CONTOUR METHOD FOR LINEAR ELASTICITY INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng. 46, 1883}1895 (1999) THE TRACTION BOUNDARY CONTOUR METHOD FOR LINEAR ELASTICITY ZHOU SHENJIE*, CAO ZHIYUAN AND SUN

More information

FEM type method for reconstruction of plane stress tensors from limited data on principal directions

FEM type method for reconstruction of plane stress tensors from limited data on principal directions Mesh Reduction Methods 57 FEM type method for reconstruction of plane stress tensors from limited data on principal directions J. Irša & A. N. Galybin Wessex Institute of Technology, Southampton, UK Abstract

More information

Comparison of Galerkin and collocation Trefftz formulations for plane elasticity

Comparison of Galerkin and collocation Trefftz formulations for plane elasticity Comparison of Galerkin and collocation Trefftz formulations for plane elasticity V.M.A. Leitão December 8, 2000 Abstract The purpose of this work is to compare and assess, more in terms of computational

More information

Rayleigh wave correction for the BEM analysis of two-dimensional elastodynamic problems in a half-space

Rayleigh wave correction for the BEM analysis of two-dimensional elastodynamic problems in a half-space INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2004; 60:1 16 [Version: 2002/09/18 v2.02] Rayleigh wave correction for the BEM analysis of two-dimensional elastodynamic

More information

Abstract. 1 Introduction

Abstract. 1 Introduction Efficiency of absorbing boundary conditions forfluid-saturatedporous media T. Akiyoshi, K. Fuchida, H.L. Fang Department of Civil and Environmental Engineering, Kumamoto University, 2-39-1 Kurokami, Kumamoto

More information

The aims of this experiment were to obtain values for Young s modulus and Poisson s ratio for

The aims of this experiment were to obtain values for Young s modulus and Poisson s ratio for The Cornu Method Nikki Truss 09369481 Abstract: The aims of this experiment were to obtain values for Young s modulus and Poisson s ratio for Perspex using the Cornu Method. A value of was found for Young

More information

Development of discontinuous Galerkin method for linear strain gradient elasticity

Development of discontinuous Galerkin method for linear strain gradient elasticity Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace

More information

Basic Theorems in Dynamic Elasticity

Basic Theorems in Dynamic Elasticity Basic Theorems in Dynamic Elasticity 1. Stress-Strain relationships 2. Equation of motion 3. Uniqueness and reciprocity theorems 4. Elastodynamic Green s function 5. Representation theorems Víctor M. CRUZ-ATIENZA

More information

Hydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition

Hydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition Fluid Structure Interaction and Moving Boundary Problems IV 63 Hydroelastic vibration of a rectangular perforated plate with a simply supported boundary condition K.-H. Jeong, G.-M. Lee, T.-W. Kim & J.-I.

More information

Stress analysis of a stepped bar

Stress analysis of a stepped bar Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.

More information

Investigation of Pile- Soil Interaction Subjected to Lateral Loads in Layered Soils

Investigation of Pile- Soil Interaction Subjected to Lateral Loads in Layered Soils American J. of Engineering and Applied Sciences (): 76-8, 008 ISSN 9-700 008 Science Publications Investigation of Pile- Soil Interaction Subjected to Lateral Loads in Layered Soils A. Avaei, Abdoul R.

More information

HIGHER-ORDER THEORIES

HIGHER-ORDER THEORIES HIGHER-ORDER THEORIES Third-order Shear Deformation Plate Theory Displacement and strain fields Equations of motion Navier s solution for bending Layerwise Laminate Theory Interlaminar stress and strain

More information

EXAMPLE OF PILED FOUNDATIONS

EXAMPLE OF PILED FOUNDATIONS EXAMPLE OF PILED FOUNDATIONS The example developed below is intended to illustrate the various steps involved in the determination of the seismic forces developed in piles during earthquake shaking. The

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

Boundary element method, finite element method and thefluxspline method: a performance comparison for scalar potential problems

Boundary element method, finite element method and thefluxspline method: a performance comparison for scalar potential problems Boundary element method, finite element method and thefluxspline method: a performance comparison for scalar potential problems A. BulcSo, C. F. Loeffler, P. C. Oliveira Departamento de Engenharia Mecdnica,

More information

3. Overview of MSC/NASTRAN

3. Overview of MSC/NASTRAN 3. Overview of MSC/NASTRAN MSC/NASTRAN is a general purpose finite element analysis program used in the field of static, dynamic, nonlinear, thermal, and optimization and is a FORTRAN program containing

More information

Earthquake response analysis of rock-fall models by discontinuous deformation analysis

Earthquake response analysis of rock-fall models by discontinuous deformation analysis c Earthquake response analysis of rock-fall models by discontinuous deformation analysis T. Sasaki, I. Hagiwara & K. Sasaki Rock Engineering Lab., Suncoh Consultants Co. Ltd., Tokyo, Japan R. Yoshinaka

More information

Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method

Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method 9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No

More information

Dynamic Response Of Laminated Composite Shells Subjected To Impulsive Loads

Dynamic Response Of Laminated Composite Shells Subjected To Impulsive Loads IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 14, Issue 3 Ver. I (May. - June. 2017), PP 108-123 www.iosrjournals.org Dynamic Response Of Laminated

More information

General elastic beam with an elastic foundation

General elastic beam with an elastic foundation General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation

More information

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method Module 2 Analysis of Statically Indeterminate Structures by the Matrix Force Method Lesson 8 The Force Method of Analysis: Beams Instructional Objectives After reading this chapter the student will be

More information

Three-dimensional thermo-mechanical analysis of layered elasticlplastic solids

Three-dimensional thermo-mechanical analysis of layered elasticlplastic solids Three-dimensional thermo-mechanical analysis of layered elasticlplastic solids W. Peng & Y.-T. Hsia Seagate Technology, U.S.A. Abstract At high temperature a layered solid undergoes intense thermal loading

More information

Computational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground

Computational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Computational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground To cite this article: Jozef Vlek and Veronika

More information

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES

CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric

More 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

Boundary Element Analysis of Shear Deformable Shallow Shells Under Harmonic Excitation

Boundary Element Analysis of Shear Deformable Shallow Shells Under Harmonic Excitation Copyright 2014 Tech Science Press CMES, vol.100, no.2, pp.105-118, 2014 Boundary Element Analysis of Shear Deformable Shallow Shells Under Harmonic Excitation J. Useche 1 Abstract: In this work, the harmonic

More information

J. Sladek, V. Sladek & M. Hrina Institute of Construction and Architecture, Slovak Academy of Sciences, Bratislava, Slovakia

J. Sladek, V. Sladek & M. Hrina Institute of Construction and Architecture, Slovak Academy of Sciences, Bratislava, Slovakia Evaluation of fracture parameters for functionally gradient materials J. Sladek, V. Sladek & M. Hrina Institute of Construction and Architecture, Slovak Academy of Sciences, 842 20 Bratislava, Slovakia

More information

Dynamic and buckling analysis of FRP portal frames using a locking-free finite element

Dynamic and buckling analysis of FRP portal frames using a locking-free finite element Fourth International Conference on FRP Composites in Civil Engineering (CICE8) 22-24July 8, Zurich, Switzerland Dynamic and buckling analysis of FRP portal frames using a locking-free finite element F.

More information

STUDY OF DYNAMIC SOIL-STRUCTURE INTERACTION OF CONCRETE GRAVITY DAMS

STUDY OF DYNAMIC SOIL-STRUCTURE INTERACTION OF CONCRETE GRAVITY DAMS STUDY OF DYNAMIC SOIL-STRUCTURE INTERACTION OF CONCRETE GRAVITY DAMS Djamel OUZANDJA 1, Fatiha BENKECHIDA 2, Toufiq OUZANDJA 3, Hamza BELHADED 4 ABSTRACT The safety evaluation of the dams subjected to

More information