A BEM implicit formulation applied to damage mechanics models

Size: px
Start display at page:

Download "A BEM implicit formulation applied to damage mechanics models"

Transcription

1 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 et Technologie, ENS de Cachan/CNRS/Universit2 Paris 6, France, Abstract In this work, the implicit BEM formulation, initially developed in the context of plasticity analysis is extended to incorporate damage mechanics models. The algebraic equations adopted for the formulation are obtained either using displacement or traction equations, for the boundary nodes, and stress (or strain) equations for the internal nodes. The consistent tangent operator has been obtained for situations: the conventional non-linear BEM system of equations. 1 Introduction Analysis of non-linear problems by means of BEM can be found since the end of seventies [l].the non-linear formulations used during a quite long time were all based on the initial stress and strain procedures, where a constant matrix scheme is employed. The consistent tangent operator (CTO), as proposed by Simo [2] for finite elements, has been introduced into elastoplas.tic BEM formulations only recently [3,4]. After that, the elastoplastic BEM formulations have been also extended to deal with localization problems in plasticity with emphasis to the regularization procedures [5,6]. In those works, the authors have shown the accuracy of the implicit formulation to compute very large deformation growing over a very

2 narrow and localized bandwish, comparing the results with solution obtained by using explicit models. They have also derived the complete implicit BEM formulation for gradient plasticity, adopted to avoid the mesh dependency. The results obtained using this formulation are accurate and clearly the mesh dependence problem is avoided. Modeling the mechanical behavior of brittle material structures is nowadays an important and interesting theme for research. Concrete, often assumed as a brittle material, is widely adopted in civil engineering construction, justifying therefore the interest on developing accurate numerical models in this context. This material exhibits a particular behavior characterized by the formation of micro-cracks resulting into losses of strength and rigidity of the structural members. This behavior is well represented by models based on the continuum damage mechanics [7,8]. So far, only limited applications of BEM to damage mechanics have been reported in the literature [g]. Damage mechanics is a complete different problem to be dealt with BEM. The Young's modulus is no longer constant, being its new value given by the adopted model that controls the degradation of the material. Besides that, using elastic prevision as in plasticity is not a convenient choice; the rigidity could be, at certain stage, so much deteriorated that elastic prevision may lead to an increasing number of iterations at each time step. In this paper, the non-linear BEM formulation is extended to solids governed by damage models, particularly proposed to deal with brittle material. First, the boundary algebraic equations are derived and then transformed appropriately to give an incremental solution scheme with tangent predictor. Some preliminary numerical results are shown to demonstrate the accuracy and stability of the developed model. 2 The continuum damage mechanics model Continuum Damage Mechanics (CDM) deals with the load carrying capacity of solids without majors cracks, but where the material itself is damaged due to the presence of microscopic defects such as microcracks and microvoids. The damage mechanics models, defined in the context of the thermodynamic of the irreversible processes, require the definition of internal variables to represent the energy dissipation process [7,8]. The internal variables are either scalar-valued ones, when isotropy is assumed (two variables can be assumed to represent the phenomenon differently when in tension or in compression), or tensor-valued ones, when anisotropic behaviors are represented. Without loss of generality, in this work, we decided to take a very simple isotropic damage model to derive BEM formulations in the context of elasticitybased damage models. The model is based on only one damage parameter D to reflect the amount of damage occurring in the material, either in compression or in tension. The corresponding constitutive relations are derived from the Helmoltz free energy y as follows,

3 Boundary Elrmcnt Tcrhnology XV 289 depending on the strain tensor E, and the damage variable D; p is the mass density and EUk+ is the forth-order tensor of the elastic sound material. The stress tensor is derived from eq. (l), resulting into, and the thermodynamically force Y, conjugated to the damage variable D, is Initiation and growth of damage is formulated through a damage criterion written in terms of the loading function f(y;d), as follows,. af D = h-(y;d) (4) ay The multiplier h, in equation (4) satisfies the Kuhn-Tucker conditions. From the consistency conditions, ' f = 0, one can write the following relation: In this work, we choose for the function f the very simple expression The history term k(d), in equation (7), measures the largest value that has been attained by the damage strain energy density rate Y. Accordingly, it is a nondecreasing function and grows only when f(y;d) = 0. An Example of k(d) is the local expression where M is a material constant. From the damage criterion f(y;d) = 0, one finds

4 3 Integral representations for non-linear problems Let us first consider an elastic body associated with many possible elastic states satisfying the Navier's equations, i.e.: where U, represents the components of the displacement field; p is the shear modulus and v is Poisson's ratio, while EIJkm stand for the elastic tensor. For a domain C2 with boundary r, standard integral representations are derived by applying Betti's principle (Green's second identity). Particularly, displacement integral representation is easily derived [l]: where uk and pk are the displacement and the traction components at boundary points, respectively, bk gives the body forces distributed over the domain, the free term C;, is the dependent upon the boundary geometry; while the symbol * indicates thiswell-known Kelvin fundamental solutions. For non-linear problems, Betti's principle cannot be directly applied. Moreover, the state variables for a non-linear model are usually history dependent. In this case, splitting the total strain into its elastic and non-linear components, i.e., ske =E& +E&, the Navier's operator applies only to the elastic part, so that equation (lo), now written in rates, becomes: -L..G.=bi-E.. 1~ J 11km km.j (12) where i:, is a tensor containing non-linear strain rates. In equation (12), the non-linear strain field acts therefore as fictitious body forces. This is referred as the initial stress method and usually the strain rate term is replaced by the initial stress rate term, given by '6; =. Thus, taking into account the last form of the Navier's equation (12), the Somigliana's identity (l l) now reads: where represents the non-linear stress rate components due to damage and is the strain filed given by the Kelvin's solution.

5 Boundary Elrmcnt Tcrhnology XV 29 1 As for elasticity, the integral representation of stress rates can be obtained by differentiating (13) at an internal point, with respect to space co-ordinates and applying Hooke's law. Thus, one obtains, where the new kernels indicated by * are obtained from the corresponding ones given in eq. (13) by applying the strain definition; bij(be) is a free-term resulting from the derivatives of the strong singular domain integral in equation (13) (e.g. Brebbia et alli [l]). In this work, we decided to not extend equation (14) to boundary points; as it will be discussed later, only internal points are used to approximate the strain corrector values. The alternative given by extending equation (14) to boundary points can be found elsewhere for the elastic case [10, 1 l], and for the non-linear problems [4,5]. Evaluating strains by mean of well-established relations: Cauchy's formula and the displacement derivatives with respect to the circumferential Cartesian coordinates (local coordinate system) are not reliable for non-linear analysis and therefore must be avoided. In this case, the cumulated errors during the iterations are enough significant to spoil the solution. For nonlinear applications exhibiting strain concentration or localization, this procedure always leads to numerical instability. In despite of computing strains only at internal points, equation (14) will be also used for outside collocations. In this case, the left hand side term vanishes leading to an expression relating only boundary values and the non-linear stress rate field, (rk. Equations (13) and (14) are exact equilibrium representations of the body under consideration. In the next section, they will be conveniently approximated to give algebraic relations written in terms of certain number of degrees of freedom to allow numerical solution of the problem. 4 Boundary element method In order to solve the boundary value problem described by the equilibrium equations (differential or integral equations) and the constrains given by the chosen damage criterion, discussed in section 2, a space discretization followed by displacement and strain field approximations must be adopted. The boundary l- is then discretized into elements, rs, along which displacements and tractions are approximated. For the examples presented in this work, only quadratic shape functions were used to describe those fields along each element. Discontinuities at some boundary points may be allowed when convenient; this is provided by defining collocation points inside the element or at any outside close points (they must be close to the boundary to assure the system stability). In order to have the integral equilibrium equations written in their algebraic forms, the strain (or stress) filed must be approximated in R. Internal cells 0, '

6 are defined, over which the strain components are approximated based on selected shape hnctions and nodal values. Triangular cells with linear approximation have been adopted with nodes defined at corner, but when convenient placed along the bisecting lines. This discretization is convenient to choose only internal nodes to approximate the strain field. When a cell has one or more nodes along the boundary, internal collocations are automatically defined. Thus, after selecting a collocation point, algebraic relations are obtained, performing properly the integrals over boundary elements r, and internal cells QC. Thus, the integral equations (13) and (14) for a single collocation point S, are transformed to: [W)]. (U)-[W)]. (*l- [~(s)]. bd]= 0 (15) where {U) and ffi} are vectors containing the boundary displacement and traction rates, respectively; (od] gives the damage stress rate vector defined at selected internal points, {b} is the stress rates computed at the collocation S; [I-~(s)], [~(s)], [~(s)],[~'(s)], [~'(s)], [Q($)] and [~'(s)] are matrices resulting by performing the appropriate integrals over boundary elements and internal cells. Note that equations (15) and (16) were obtained for any selected collocation point S; Collocations point may be placed at internal position, on the boundary (equation (14) must be derived for this situation) or external position. In order to have a closed system of equations to solve the problem, one must choose a number of relations equal the number of problem unknowns, i.e., boundary displacernents and tractions plus the stresses at the internal nodes defined by the discretization. Therefore, writing equation (16) for external or boundary collocation points give also independent relation, as follows: Using equation (17) to replace equation (15) defines an alternative scheme that has been experimented in this work to solve the proposed damage problem, However, for this particular case under analysis, using only stress equations, therefore only hypersingular equations, did not lead to a better performance of the method. The final results were very similar and no advantage was observed. After writing equations (15) and (16), or equations (17) and (16) if only hypersingular equations are chosen, in their incremental forms (time integration to be discussed in the next section), for a convenient number of collocations, the following matrix blocks of algebraic equations are found:

7 Boundaq Element Technologv XI/ 293 It is to be noted that the block of equations (18) can be written from (15), if the standard BEM is adopted, or from equation (17) if the hypersingular formulation is the chosen alternative. Equations (18) and (19) represent the approximated equilibrium of the solid during a time increment At. In the next section they will be used to derive an implicit BEM algorithm to deal with problem in the damage mechanics context. 5 Implicit formulation of BEM for damage mechanics The incremental forms of equations (18) and (19), after applying the boundary conditions are given by l l+ [Q]. bo: / (20) where {AX, ] contains all boundary unknowns {~f, ] and {Agn) are given by the prescribed boundary values and a time increment At, has been considered. Then, equations (20) and (21) can be solved to give: where {Ah~l,f=[~]-'{~f~) and {AN~]=-[A'].[A]-'.(A~~)+{A~,~ (24a,b) [R]=[A]-'.[Q] and [S]=-[A'].[A]-'.[Q]+[Q1] (2w) It is important noting that the equilibrium equation is now reduced to equation (23). It can be rearranged to be expressed in terms of the strain increment AE, = &,+l - E, (note that As, = ~efl i- As f,, As:, damaged and elastic strain parts, respectively) and the increments of the n considered internal variables, i.e., Aakn =akn+l -akn, as ~O~~OWS,

8 294 Boundary Element Technologv XV where [S] = [S]+ [I 1, being [I ] the identity matrix. From the damage model constitutive equations, one may obtain explicitly the stress increment in terms of the strain increment, therefore, equation (26) reads, Equation (27) can now be solved by applying the Newton-Raphson scheme. An iterative process may be required to achieve the equilibrium. Then, from the iteration i the next try, (i+l), for the time increment At, is given by By linearizing the equilibrium equation, (27), using the fxst term of the Taylor's expansion, gives, Thus, the equilibrium is now expressed by where the algorithmic consistent tangent operator, ccto = ~ AcT/~AE is explicated. Although the presented scheme is similar to the one proposed for plasticity [5,6], there are some strong differences between them. In plasticity, the first try of each increment is always the elastic one, which makes the algorithm simpler. For damage, using the last tangent matrix is recommended, otherwise the number of iterations required for convergence increases too much. Rigidity at damaged regions is so severely reduced that elastic tries are always an inconvenient choice. In order to use reduced rigidity matrix, one has to modify the nth elastic increment that now reds and the equilibrium Then, the linearization follows the relations already presented, equations (2% (291, (30). Another important aspect shown by structures made with materials governed by damage laws is the presence of small and localized regions exhibiting severe reduction of rigidity. This fact usually induces the snap-

9 Boundaiy Element Technologv XV 295 back effects. The dissipation zone is so reduced that both referenced displacement and the total applied load are reduced simultaneously. In order to deal with that particular behavior of some damaged domains, we have incorporated into the implicit formulation presented in the previous sections the arc-length procedure. The formulation of this BEM arc-length will be presented in a separated paper [12]. The strategy for this class of methods is to assume as guide variable a specified length to go along the curve load X displacement. We define a priori the size of the chord to reach next point along the curve. With this technique the guide variable is neither the loading nor the displacement. Increments of the displacement X reaction curve chord are given as a prescribed load; the arc equation is added to the system of algebraic equations to complete the necessary relations. 6 Numerical examples In this section, two examples are presented to illustrate the applicability of the BEM to damaged domains. The first example is a one-dimensional test to show how the proposed method can model accurately the deteriorated (damaged) zone concentrated over a small region at the domain center. The geometry of the chosen thin domain (a X L) is presented in Figure la, together with the boundary conditions given by prescribing displacements along the left and right sides. The left side displacements are zero, while along the right side they are controlled. The discretization is shown in Figure lb, for which a small rectangle at the center will be assumed weaker, the region where the material will deteriorate (localized region). The damage model adopted is the one defined in section 2, for which the following parameters were adopted: Yo= and M= For the weaker region the following values were assumed: Yo= and M= Figure 2 shows two curves captured for different sizes of the central deteriorated regions (dl=1/13 and a/l=1/27, respectively). This solution demonstrates clearly that the solution is not unique; it is dependent on the size of the dissipation zone. As can be seen for smaller weaker central region more pronounced is the snap-back effects. Using this simple example the mesh dependence of this problem can be easily observed. Dividing the weaker zone in several rows of cells, the damage region will localize only over one of them, leading therefore to different curves when the mesh is refmed. t I - ' cm r-? \ 7 I \,,,jl l l l l W, a - L Figure 1. one-dimensional domain. Geometry and discretization.

10 296 Boundary Element Technologv XV 1 Force X Displacement displacement (cm) (X 0.01) Figure 2. Displacement X reaction curves for the 1D case The problem shown in the frst example exhibits the behavior of a onedimensional case. This analysis can be extended to see how the formulation behaves for problems where the two-dimension effects are clearly introduced. Then, we have chosen another rectangular domain that exhibits these characteristics, as shown in Figure 3 where the geometric data are specified. The loading is characterized by uniform distributed forces applied to the two opposite domain ends. The same parameters adopted for the first example have been used for this case. A square containing one node at the center of the rectangular domain has been assumed as a weaker zone to avoid capturing the uniform stiffness reduction solution. Two discretizations (Figure 4) are used to obtain the numerical solution to illustrate the mesh dependence for this problem. The energy dissipation occurs over the central row of cells as expected. As this area is very small, in comparison with the size of the domain, the snap-back effects are captured (see Figure 5). These effects are much more pronounced when the energy dissipation area is strongly reduced, which is the case observed for the fine mesh. The two solutions presented are only numerical solutions achieved for these specific discretizations. The solutions exhibit they particular characteristic, depending on the adopted mesh. The sequence of the damage appearance modifies abruptly when the discretization is little modified. Other solutions may be found for this case, depending on the approximations, positions of the internal nodes and other geometric data. In order to have single solution independent on the mesh, one must use any regularization technique. The complete formulation combining the implicit BEM presented here with a regularization technique is given in reference [ 121.

11 Boundaiy Element Technologv XV 297 Figure 3. 2D domain. weak points Figure 4. Discretizations and weak points. Force X Displacement -0- Figura 4a displacement (cm)(x 0.01) Figure 5. Displacement X reaction curve for the 2D case. 7 Conclusions In this work, the boundary element method has been extended to analyze damaged domains. The BEM implicit formulation for a simple damage criterion

12 298 Boundary Element Technologv XV has been presented and used to show the accuracy of the method for this class of problems. The mesh dependence of the local model is clearly represented by the method. For unique solution a regularization technique must be incorporated. For brittle material, the snap-back effects are captured using the presented algorithm, which has incorporated the arc-length scheme. Acknowledgement To FAPESP - STio Paulo State Research Foundation. References [l] Brebbia, C.A., Telles, J.C.F & Wrobel, L.C. Boundary element techniques. Theory and applications in engineering, Springer-Verlag: Berlin and New York, [2] Sirno, J.C. & Taylor, R.L., Consistent tangent operators for rate-independent elastoplasticity. Comp. Meth. Appl. Mech. Engng, 48, pp , [3] Bonnet, M. & Mukherjee, S., Implicit BEM formulation for usual and sensitivity problems in elasto-plasticity using the consistent tangent operator concept. Int. J. Solids and Structures, 33, pp , [4] Poon, H., Mukherjee, S. & Bonnet, M., Numerical implementation of a CTO-based implicit approach for solution of usual and sensitivity problems in elasto-plasticity. Engineering Analysis with Boundary Elements, 22, pp , [5] Fudoli, C.A., Venturini, W.S. & Benallal, A., BEM Formulations Applied to Gradient Plasticity. Boundary Elements 22, Eds. Brebbia, C.A. & Power,WIT Publ. Southampton, [6] Benallal, A., Fudoli, C.A. & Venturini, W.S., An implicit BEM formulation for Gradient plasticity and localization phenomena. Int. J Num. Meth. Eng., 53 pp ,2002. [7] Lemaitre, J., A course on damage mechanics, Springer-Verlag, [S] Lemaitre, J. & Chaboche J. L., Me'canique de mate'riaux solides, Dunod- Bordas, Paris, [9] Garcia R, Florez-Lopez J, Cerrolaza M., A boundary element formulation for a class of non-local damage models, Int. J. S. Struct., 36 pp , [loiguiggiani, M. Hypersingular formulation for boundary stress evaluation. Engng. Analysis with Boundary Element, 13, pp , [l l]cruse, T.A. &. Richardson, J.D. Non-singular Somigliana stress identities in elasticity. Int. J. Num. Meth. Engng., 39, pp , [12]Benallal, A., Botta, A.S. & Venturini, W.S., On the description of localization and failure phenomena by the boundary element method. Int. J. Num. Meth. Eng., 2002 (in preparation).

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

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

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

Transactions on Modelling and Simulation vol 12, 1996 WIT Press,   ISSN X Plate-soil elastodynamic coupling using analysis S.F.A. Baretto, H.B. Coda, W.S. Venturini Sao Carlos School of Engineering, University ofsao Paulo, Sao Carlos - SP, Brazil BEM Abstract The aim of this

More information

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

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

On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity

On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity 63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain

More information

Code_Aster. Reaction of the steel subjected to corrosion

Code_Aster. Reaction of the steel subjected to corrosion Titre : Comportement de l acier soumis à la corrosion Date : 16/0/01 Page : 1/1 Reaction of the steel subjected to corrosion Summarized: This documentation the model presents steel reaction of subjected

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

NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS

NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS ABSTRACT : P Mata1, AH Barbat1, S Oller1, R Boroschek2 1 Technical University of Catalonia, Civil Engineering

More information

DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS

DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the

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

Nonlocal computational methods applied to composites structures

Nonlocal computational methods applied to composites structures Nonlocal computational methods applied to composites structures Norbert Germain, Frédéric Feyel, Jacques Besson To cite this version: Norbert Germain, Frédéric Feyel, Jacques Besson. Nonlocal computational

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

Fluid driven cohesive crack propagation in quasi-brittle materials

Fluid driven cohesive crack propagation in quasi-brittle materials Fluid driven cohesive crack propagation in quasi-brittle materials F. Barpi 1, S. Valente 2 Department of Structural and Geotechnical Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129

More information

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

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

More information

3D dynamic crack propagation analysis with PDS-FEM

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

More information

Structural behaviour of traditional mortise-and-tenon timber joints

Structural behaviour of traditional mortise-and-tenon timber joints Structural behaviour of traditional mortise-and-tenon timber joints Artur O. Feio 1, Paulo B. Lourenço 2 and José S. Machado 3 1 CCR Construtora S.A., Portugal University Lusíada, Portugal 2 University

More information

Abstract. 1 Introduction

Abstract. 1 Introduction Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut

More information

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

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

More information

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

Module 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method

Module 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method Introduction In this lecture we will introduce some more micromechanical methods to predict the effective properties of the composite. Here we will introduce expressions for the effective properties without

More information

Cohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition framework

Cohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition framework University of Liège Aerospace & Mechanical Engineering MS3: Abstract 131573 - CFRAC2017 Cohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition

More information

The Finite Element method for nonlinear structural problems

The Finite Element method for nonlinear structural problems Non Linear Computational Mechanics Athens MP06/2013 The Finite Element method for nonlinear structural problems Vincent Chiaruttini vincent.chiaruttini@onera.fr Outline Generalities on computational strategies

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

Transactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN

Transactions on Engineering Sciences vol 6, 1994 WIT Press,   ISSN Spatial discretization of strain localization E. Rizzi,* K. Willam University of Colorado Boulder, CEAE Department, Boulder, Colorado 80309-0428, USA * On leave from Politecnico di Milano, DIS, Milano,

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

MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD

MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD MODELING OF ELASTO-PLASTIC MATERIALS IN FINITE ELEMENT METHOD Andrzej Skrzat, Rzeszow University of Technology, Powst. Warszawy 8, Rzeszow, Poland Abstract: User-defined material models which can be used

More information

Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5

Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 17, 2017, Lesson 5 Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 17, 2017, Lesson 5 1 Politecnico di Milano, February 17, 2017, Lesson 5 2 Outline

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

ELASTOPLASTICITY THEORY by V. A. Lubarda

ELASTOPLASTICITY THEORY by V. A. Lubarda ELASTOPLASTICITY THEORY by V. A. Lubarda Contents Preface xiii Part 1. ELEMENTS OF CONTINUUM MECHANICS 1 Chapter 1. TENSOR PRELIMINARIES 3 1.1. Vectors 3 1.2. Second-Order Tensors 4 1.3. Eigenvalues and

More information

Course in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg

Course in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg Course in Nonlinear FEM Geometric nonlinearity Nonlinear FEM Outline Lecture 1 Introduction Lecture 2 Geometric nonlinearity Lecture 3 Material nonlinearity Lecture 4 Material nonlinearity it continued

More information

APPLICATION OF DAMAGE MODEL FOR NUMERICAL DETERMINATION OF CARRYING CAPACITY OF LARGE ROLLING BEARINGS

APPLICATION OF DAMAGE MODEL FOR NUMERICAL DETERMINATION OF CARRYING CAPACITY OF LARGE ROLLING BEARINGS INTERNATIONAL ESIGN CONFERENCE - ESIGN ubrovnik, May 14-17,. APPLICATION OF AMAGE MOEL FOR NUMERICAL ETERMINATION OF CARRYING CAPACITY OF LARGE ROLLING BEARINGS Robert Kunc, Ivan Prebil, Tomaž Rodic and

More information

Esben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer

Esben Byskov. Elementary Continuum. Mechanics for Everyone. With Applications to Structural Mechanics. Springer Esben Byskov Elementary Continuum Mechanics for Everyone With Applications to Structural Mechanics Springer Contents Preface v Contents ix Introduction What Is Continuum Mechanics? "I Need Continuum Mechanics

More information

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

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

More information

6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS

6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS 6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS Blondet et al. [25] carried out a cyclic test on an adobe wall to reproduce its seismic response and damage pattern under in-plane loads. The displacement

More information

Nigerian Journal of Technology, Vol. 26, No. 2, June 2007 Edelugo 37

Nigerian Journal of Technology, Vol. 26, No. 2, June 2007 Edelugo 37 Nigerian Journal of Technology, Vol. 26, No. 2, June 2007 Edelugo 37 APPLICATION OF THE REISSNERS PLATE THEORY IN THE DELAMINATION ANALYSIS OF A THREE-DIMENSIONAL, TIME- DEPENDENT, NON-LINEAR, UNI-DIRECTIONAL

More information

Nonlinear bending analysis of laminated composite stiffened plates

Nonlinear bending analysis of laminated composite stiffened plates Nonlinear bending analysis of laminated composite stiffened plates * S.N.Patel 1) 1) Dept. of Civi Engineering, BITS Pilani, Pilani Campus, Pilani-333031, (Raj), India 1) shuvendu@pilani.bits-pilani.ac.in

More information

Solution of contact problems in linear elasticity using a feasible interior point algorithm for nonlinear complementarity problems

Solution of contact problems in linear elasticity using a feasible interior point algorithm for nonlinear complementarity problems 7 th World Congress on Structural and Multidisciplinary Optimization COEX Seoul, May - 5 May 007, Korea Solution of contact problems in linear elasticity using a feasible interior point algorithm for nonlinear

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

Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises

Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009

More information

Combined Isotropic-Kinematic Hardening Laws with Anisotropic Back-stress Evolution for Orthotropic Fiber-Reinforced Composites

Combined Isotropic-Kinematic Hardening Laws with Anisotropic Back-stress Evolution for Orthotropic Fiber-Reinforced Composites Combined Isotropic-Kinematic Hardening Laws with Antropic Back-stress Evolution for Orthotropic Fiber- Reinforced Composites Combined Isotropic-Kinematic Hardening Laws with Antropic Back-stress Evolution

More information

An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation

An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation Nachiket Patil, Deepankar Pal and Brent E. Stucker Industrial Engineering, University

More information

International Journal of Pure and Applied Mathematics Volume 58 No ,

International Journal of Pure and Applied Mathematics Volume 58 No , International Journal of Pure and Applied Mathematics Volume 58 No. 2 2010, 195-208 A NOTE ON THE LINEARIZED FINITE THEORY OF ELASTICITY Maria Luisa Tonon Department of Mathematics University of Turin

More information

Computational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem

Computational Inelasticity FHLN05. Assignment A non-linear elasto-plastic problem Computational Inelasticity FHLN05 Assignment 2018 A non-linear elasto-plastic problem General instructions A written report should be submitted to the Division of Solid Mechanics no later than November

More information

An Atomistic-based Cohesive Zone Model for Quasi-continua

An Atomistic-based Cohesive Zone Model for Quasi-continua An Atomistic-based Cohesive Zone Model for Quasi-continua By Xiaowei Zeng and Shaofan Li Department of Civil and Environmental Engineering, University of California, Berkeley, CA94720, USA Extended Abstract

More information

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004

Elements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004 Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic

More information

Numerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour

Numerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour 13 th International Brick and Block Masonry Conference Amsterdam, July 4-7, 2004 Numerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour H. R.

More information

Analytical formulation of Modified Upper Bound theorem

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

More information

Fatigue Damage Development in a Steel Based MMC

Fatigue Damage Development in a Steel Based MMC Fatigue Damage Development in a Steel Based MMC V. Tvergaard 1,T.O/ rts Pedersen 1 Abstract: The development of fatigue damage in a toolsteel metal matrix discontinuously reinforced with TiC particulates

More information

Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface

Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface Transaction A: Civil Engineering Vol. 16, No. 6, pp. 512{519 c Sharif University of Technology, December 2009 Research Note Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield

More information

A hierarchy of higher order and higher grade continua Application to the plasticity and fracture of metallic foams

A hierarchy of higher order and higher grade continua Application to the plasticity and fracture of metallic foams A hierarchy of higher order and higher grade continua Application to the plasticity and fracture of metallic foams Samuel Forest Centre des Matériaux/UMR 7633 Mines Paris ParisTech /CNRS BP 87, 91003 Evry,

More information

Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics

Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics Tobias Gasch *1 and Anders Ansell 1 1 KTH Royal Institute of Technology, Department of Civil and Architectural Engineering *Corresponding

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

Table of Contents. Preface...xvii. Part 1. Level

Table of Contents. Preface...xvii. Part 1. Level Preface...xvii Part 1. Level 1... 1 Chapter 1. The Basics of Linear Elastic Behavior... 3 1.1. Cohesion forces... 4 1.2. The notion of stress... 6 1.2.1. Definition... 6 1.2.2. Graphical representation...

More information

Prediction of Elastic Constants on 3D Four-directional Braided

Prediction of Elastic Constants on 3D Four-directional Braided Prediction of Elastic Constants on 3D Four-directional Braided Composites Prediction of Elastic Constants on 3D Four-directional Braided Composites Liang Dao Zhou 1,2,* and Zhuo Zhuang 1 1 School of Aerospace,

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

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

ALGORITHM FOR NON-PROPORTIONAL LOADING IN SEQUENTIALLY LINEAR ANALYSIS

ALGORITHM FOR NON-PROPORTIONAL LOADING IN SEQUENTIALLY LINEAR ANALYSIS 9th International Conference on Fracture Mechanics of Concrete and Concrete Structures FraMCoS-9 Chenjie Yu, P.C.J. Hoogenboom and J.G. Rots DOI 10.21012/FC9.288 ALGORITHM FOR NON-PROPORTIONAL LOADING

More information

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

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

More information

Basic Energy Principles in Stiffness Analysis

Basic Energy Principles in Stiffness Analysis Basic Energy Principles in Stiffness Analysis Stress-Strain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting

More information

Constitutive models: Incremental (Hypoelastic) Stress- Strain relations. and

Constitutive models: Incremental (Hypoelastic) Stress- Strain relations. and Constitutive models: Incremental (Hypoelastic) Stress- Strain relations Example 5: an incremental relation based on hyperelasticity strain energy density function and 14.11.2007 1 Constitutive models:

More information

BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS

BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department

More information

Analysis of a portal steel frame subject to fire by use of a truss model

Analysis of a portal steel frame subject to fire by use of a truss model Analysis of a portal steel frame subject to fire by use of a truss model P. G. Papadopoulos & A. Mathiopoulou Department of Civil Engineering, Aristotle University of Thessaloniki, Greece Abstract A plane

More information

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit

More information

Transactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN

Transactions on Engineering Sciences vol 6, 1994 WIT Press,   ISSN A computational method for the analysis of viscoelastic structures containing defects G. Ghazlan," C. Petit," S. Caperaa* " Civil Engineering Laboratory, University of Limoges, 19300 Egletons, France &

More information

MODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Class Meeting #5: Integration of Constitutive Equations

MODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Class Meeting #5: Integration of Constitutive Equations MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #5: Integration of Constitutive Equations Structural Equilibrium: Incremental Tangent Stiffness

More information

Heterogeneous structures studied by interphase elasto-damaging model.

Heterogeneous structures studied by interphase elasto-damaging model. Heterogeneous structures studied by interphase elasto-damaging model. Giuseppe Fileccia Scimemi 1, Giuseppe Giambanco 1, Antonino Spada 1 1 Department of Civil, Environmental and Aerospace Engineering,

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 06 In the last lecture, we have seen a boundary value problem, using the formal

More information

ENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material

ENGN 2340 Final Project Report. Optimization of Mechanical Isotropy of Soft Network Material ENGN 2340 Final Project Report Optimization of Mechanical Isotropy of Soft Network Material Enrui Zhang 12/15/2017 1. Introduction of the Problem This project deals with the stress-strain response of a

More information

Application of a non-local failure criterion to a crack in heterogeneous media S. Bavaglia*, S.E. Mikhailov*

Application of a non-local failure criterion to a crack in heterogeneous media S. Bavaglia*, S.E. Mikhailov* Application of a non-local failure criterion to a crack in heterogeneous media S. Bavaglia*, S.E. Mikhailov* University of Perugia, Italy Email: mic@unipg.it ^Wessex Institute of Technology, Ashurst Lodge,

More information

Microstructural effect on the cavity expansion in a soil cylinder

Microstructural effect on the cavity expansion in a soil cylinder Microstructural effect on the cavity expansion in a soil cylinder J.D. Zhao, D.C. Sheng & S.W. Sloan Centre for Geotechnical and Materials Modelling, School of Engineering University of Newcastle, Callaghan,

More information

MODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Uniaxial Model: Strain-Driven Format of Elastoplasticity

MODELING OF CONCRETE MATERIALS AND STRUCTURES. Kaspar Willam. Uniaxial Model: Strain-Driven Format of Elastoplasticity MODELING OF CONCRETE MATERIALS AND STRUCTURES Kaspar Willam University of Colorado at Boulder Class Meeting #3: Elastoplastic Concrete Models Uniaxial Model: Strain-Driven Format of Elastoplasticity Triaxial

More information

Geometrically exact beam dynamics, with and without rotational degree of freedom

Geometrically exact beam dynamics, with and without rotational degree of freedom ICCM2014 28-30 th July, Cambridge, England Geometrically exact beam dynamics, with and without rotational degree of freedom *Tien Long Nguyen¹, Carlo Sansour 2, and Mohammed Hjiaj 1 1 Department of Civil

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

Checkerboard instabilities in topological shape optimization algorithms

Checkerboard instabilities in topological shape optimization algorithms Checkerboard instabilities in topological shape optimization algorithms Eric Bonnetier, François Jouve Abstract Checkerboards instabilities for an algorithm of topological design are studied on a simple

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

Transactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN

Transactions on Engineering Sciences vol 6, 1994 WIT Press,   ISSN Large strain FE-analyses of localized failure in snow C.H. Liu, G. Meschke, H.A. Mang Institute for Strength of Materials, Technical University of Vienna, A-1040 Karlsplatz 13/202, Vienna, Austria ABSTRACT

More information

Stress intensity factor analysis for an interface crack between dissimilar isotropic materials

Stress intensity factor analysis for an interface crack between dissimilar isotropic materials Stress intensity factor analysis for an interface crack between dissimilar isotropic materials under thermal stress T. Ikeda* & C. T. Sun* I Chemical Engineering Group, Department of Materials Process

More information

BRIDGES BETWEEN DAMAGE AND FRACTURE MECHANICS

BRIDGES BETWEEN DAMAGE AND FRACTURE MECHANICS BRIDGES BETWEEN DAMAGE AND FRACTURE MECHANICS Jacky Mazars and Gilles Pijaudier-Cabot,,. Laboratoire de Mecanique et Technologie - Ecole Normale Supeneure 94235 Cachan-France Abstract Fracture mechanics

More information

GEO E1050 Finite Element Method Autumn Lecture. 9. Nonlinear Finite Element Method & Summary

GEO E1050 Finite Element Method Autumn Lecture. 9. Nonlinear Finite Element Method & Summary GEO E1050 Finite Element Method Autumn 2016 Lecture. 9. Nonlinear Finite Element Method & Summary To learn today The lecture should give you overview of how non-linear problems in Finite Element Method

More information

ENGN 2290: Plasticity Computational plasticity in Abaqus

ENGN 2290: Plasticity Computational plasticity in Abaqus ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element

More information

Continuum damage growth analysis using element free Galerkin method

Continuum damage growth analysis using element free Galerkin method Sādhanā Vol. 35, Part 3, June 2010, pp. 279 301. Indian Academy of Sciences Continuum damage growth analysis using element free Galerkin method COARUN,BNRAO and S M SRINIVASAN Structural Engineering Division,

More information

The Finite Element Method II

The Finite Element Method II [ 1 The Finite Element Method II Non-Linear finite element Use of Constitutive Relations Xinghong LIU Phd student 02.11.2007 [ 2 Finite element equilibrium equations: kinematic variables Displacement Strain-displacement

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how

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

LOCKING AND STABILITY OF 3D TEXTILE COMPOSITE REINFORCEMENTS DURING FORMING

LOCKING AND STABILITY OF 3D TEXTILE COMPOSITE REINFORCEMENTS DURING FORMING 10th International Conference on Composite Science and Technology ICCST/10 A.L. Araújo, J.R. Correia, C.M. Mota Soares, et al. (Editors) IDMEC 015 LOCKING AND STABILITY OF 3D TEXTILE COMPOSITE REINFORCEMENTS

More information

A MIXED TRIANGULAR ELEMENT FOR THE ELASTO-PLASTIC ANALYSIS OF KIRCHHOFF PLATES

A MIXED TRIANGULAR ELEMENT FOR THE ELASTO-PLASTIC ANALYSIS OF KIRCHHOFF PLATES COMPUTATIONAL MECHANICS New Trends and Applications S. Idelsohn, E. Oñate and E. Dvorkin (Eds.) c CIMNE, Barcelona, Spain 1998 A MIXED TRIANGULAR ELEMENT FOR THE ELASTO-PLASTIC ANALYSIS OF KIRCHHOFF PLATES

More information

Fundamentals of Linear Elasticity

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

More information

Modeling issues of the FRP detachment phenomenon

Modeling issues of the FRP detachment phenomenon Modeling issues of the FRP detachment phenomenon Elio Sacco in collaboration with: J. Toti, S. Marfia and E. Grande 1 Dipartimento di ngegneria Civile e Meccanica Università di Cassino e del Lazio Meridionale

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

1 Introduction IPICSE-2016

1 Introduction IPICSE-2016 (06) DOI: 0.05/ matecconf/06860006 IPICSE-06 Numerical algorithm for solving of nonlinear problems of structural mechanics based on the continuation method in combination with the dynamic relaxation method

More information

A Method for Gradient Enhancement of Continuum Damage Models

A Method for Gradient Enhancement of Continuum Damage Models TECHNISCHE MECHANIK, Band 28, Heft 1, (28), 43 52 Manuskripteingang: 31. August 27 A Method for Gradient Enhancement of Continuum Damage Models B. J. Dimitrijevic, K. Hackl A method for the regularization

More information

2D Kirchhoff Thin Beam Elements

2D Kirchhoff Thin Beam Elements 2D Kirchhoff Thin Beam Elements General Element Name Y,v X,u BM3 2 3 1 Element Group Element Subgroup Element Description Number Of Nodes 3 Freedoms Node Coordinates Geometric Properties Kirchhoff Parabolically

More information

Structural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).

Structural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free

More information

Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics

Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics Tobias Gasch, PhD Student Co-author: Prof. Anders Ansell Comsol Conference 2016 Munich 2016-10-12 Contents Introduction Isotropic damage

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

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

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

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by

SEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e

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