Transactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
|
|
- Austen Farmer
- 5 years ago
- Views:
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
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 informationTransactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
An alternative boundary element formulation for plate bending analysis J.B. Paiva, L.O. Neto Structures Department, Sao Carlos Civil Engineering School, f, BrazzY Abstract This work presents an alternative
More informationPlate analysis using classical or Reissner- Mindlin theories
Plate analysis using classical or Reissner- Mindlin theories L. Palermo Jr. Faculty of Civil Engineering, State Universiv at Campinas, Brazil Abstract Plates can be solved with classical or Reissner-Mindlin
More informationAbstract. 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 informationAn alternative multi-region BEM technique for layered soil problems
An alternative multi-region BM technique for layered soil problems D.B. Ribeiro & J.B. Paiva Structural ngineering Department, São Carlos ngineering School, University of São Paulo, Brazil. Abstract Different
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationAppendix 1. Numerical Integration: Regular Integrals
The Dual Reciprocity 269 Appendix 1 Numerical Integration: Regular Integrals 1 Introduction The numerical integration formulae given in this Appendix are used for integration over boundary elements which
More informationTransactions on Modelling and Simulation vol 8, 1994 WIT Press, ISSN X
Model analysis of plates using the dual reciprocity boundary element method T.W. Davies & F.A. Moslehy Department of Mechanical and Aerospace Engineering, University of Central Florida, Orlando, Florida,
More informationDiscrete 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 informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationINTERNATIONAL 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 informationA 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 informationelastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK
Non-conforming BEM elastoplastic contact problems D. Martin and M.H. Aliabadi discretisation in Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK Abstract In this paper,
More informationA 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 informationMIXED 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 informationAN 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 information1 Numerical combination for nonlinear analysis of structures coupled to layered soils
1(2012) 1 23 1 Numerical combination for nonlinear analysis of structures coupled to layered soils Abstract This paper presents an alternative coupling strategy between the Boundary Element Method (BEM)
More informationUsing 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 informationSIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION
8 th GRACM International Congress on Computational Mechanics Volos, 12 July 15 July 2015 SIMULATION OF PLANE STRAIN FIBER COMPOSITE PLATES IN BENDING THROUGH A BEM/ACA/HM FORMULATION Theodore V. Gortsas
More informationCOPYRIGHTED 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 informationMechanics 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 informationEngineering. 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 informationFast 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 informationApplication 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 informationTransactions on Modelling and Simulation vol 2, 1993 WIT Press, ISSN X
A general integral-equation formulation for anisotropic elastostatics M.M. Perez"'*, L.C. Wrobel* * Wessex Institute of Technology, University of Portsmouth, Ashurst Lodge, Ashurst, Southampton S04 2AA,
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationFLEXIBILITY 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 informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationVibration 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 informationBending 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 informationStructural Acoustics Applications of the BEM and the FEM
Structural Acoustics Applications of the BEM and the FEM A. F. Seybert, T. W. Wu and W. L. Li Department of Mechanical Engineering, University of Kentucky Lexington, KY 40506-0046 U.S.A. SUMMARY In this
More informationTransactions 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 informationFig. 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 informationStatic & 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 informationHANDBUCH 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 informationLevel 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 informationDYNAMIC 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 informationA *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 informationEDEM 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 informationFirst- 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 informationINFLUENCE 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 informationBHAR 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 informationTransactions on Modelling and Simulation vol 9, 1995 WIT Press, ISSN X
Analysis of complexfluid-structureinteraction problems using a triangular boundary element for the fluid dynamics A. Hopf, G. Hailfmger, R. Krieg Forschungszentrum Karlsruhe GmbH, Institutfur D- 7^027
More informationTable 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 informationSchur 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 informationDYNAMIC 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 informationLIMIT 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 informationMixed 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 informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More information202 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 informationBack Matter Index The McGraw Hill Companies, 2004
INDEX A Absolute viscosity, 294 Active zone, 468 Adjoint, 452 Admissible functions, 132 Air, 294 ALGOR, 12 Amplitude, 389, 391 Amplitude ratio, 396 ANSYS, 12 Applications fluid mechanics, 293 326. See
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationSolving 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 informationDEFORMATION 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 informationPiezoelectric Bimorph Response with Imperfect Bonding Conditions
Copyright c 28 ICCES ICCES, vol.6, no.3, pp.5-56 Piezoelectric Bimorph Response with Imperfect Bonding Conditions Milazzo A., Alaimo A. and Benedetti I. Summary The effect of the finite stiffness bonding
More informationPerformance 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 informationBack Analysis of Measured Displacements of Tunnels
Rock Mechanics and Rock Engineering 16, 173--180 (1983) Rock Mechanics and Rock Engineering 9 by Springer-Verlag 1983 Back Analysis of Measured Displacements of Tunnels By S. Sakurai and K. Takeuchi Kobe
More informationDISPENSA FEM in MSC. Nastran
DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation
More informationAN 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 informationChapter 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 informationChapter 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 informationTransactions on Modelling and Simulation vol 12, 1996 WIT Press, ISSN X
Simplifying integration for logarithmic singularities R.N.L. Smith Department ofapplied Mathematics & OR, Cranfield University, RMCS, Shrivenham, Swindon, Wiltshire SN6 SLA, UK Introduction Any implementation
More informationCRITERIA FOR SELECTION OF FEM MODELS.
CRITERIA FOR SELECTION OF FEM MODELS. Prof. P. C.Vasani,Applied Mechanics Department, L. D. College of Engineering,Ahmedabad- 380015 Ph.(079) 7486320 [R] E-mail:pcv-im@eth.net 1. Criteria for Convergence.
More informationTIME-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 informationD : 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 informationSound 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 informationPREPRINT 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 informationCOUPLED 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 informationLecture 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 informationIraq 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 informationTHE 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 informationFEM 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 informationComparison 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 informationRayleigh 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 informationAbstract. 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 informationThe 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 informationDevelopment 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 informationBasic 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 informationHydroelastic 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 informationStress 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 informationInvestigation 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 informationHIGHER-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 informationEXAMPLE 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 informationTransactions 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 informationBoundary 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 information3. 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 informationEarthquake 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 informationPost 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 informationDynamic 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 informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationModule 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 informationThree-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 informationComputational 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 informationCHAPTER 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 informationStresses 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 informationBoundary 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 informationJ. 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 informationDynamic 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 informationSTUDY 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