Simulation of the Crimping Process by Implicit and Explicit Finite Element Methods

Size: px
Start display at page:

Download "Simulation of the Crimping Process by Implicit and Explicit Finite Element Methods"

Transcription

1 Simulation of the Crimping Process by Implicit and Explicit Finite Element Methods Stéphane Kugener AMP France ABSTRACT Calculation procedures for implicit and explicit codes applied to insulation crimping are presented and discussed. Some known equations are reviewed as a basis for discussion and to support engineering statements. Implicit finite element simulations were run under ANSYS 1, explicit simulations under OPTRIS 2. Critical aspects of the simulation are identified and discussed for each code. Strategies for the analysis are suggested and compared. For three cases of industrial importance correlation between theoretically predicted behavior and that of actual parts is established and discussed for two wire sizes and two different tool geometries. INTRODUCTION High-level computer-aided engineering (CAE) has become a relatively common activity. A large variety of codes for mechanical, thermal, fluid dynamics, or general process simulation are available at a reasonable cost, but their use also requires a reasonable level of expertise. Only a few years ago, many algorithms used in these codes were considered proprietary research tools. The latest trend in CAD/CAM systems technology, particularly the full integration of CAE tools, linear static analyses and plasticity analyses increased the level of computing aids available to the designer. CAE technology has past the step the analysis was based primarily on the geometry of the final product. The engineer can now include material behavior and process parameters in simulations of expected product or process performance. The move beyond linear approaches brought to an end the convenience of direct solutions and pushbutton analyses, and introduced the convergence uncertainties associated with nonlinear algorithms. Now, mathematical simulation of engineering problems can become a numerical adventure, with system crashes and restarts, a solution is obtained only through iterative steps and a rigorous calculation strategy. Clear identification of potential numerical difficulties and increased experience with the adapted codes will result in more efficient procedures and increased benefits. Application of these general viewpoints to a concrete problem of industrial importance is shown in this paper. Because of its complexity and technological importance the crimping process, here termed insulation crimping, is used as an example. It involves a metal forming process with path-dependent properties and requires highly nonlinear calculations. Therefore, it offers also an excellent opportunity to compare implicit and explicit finite element formulations. As a specific example, insulation crimping of a MIC I contact for two wire gauges of 1 mm 2 and 3 mm 2 was chosen. The contact is shown in Figure 1. To verify the validity of the mathematical simulations, tests were conducted at the manufacturing level on a T2 application tool, with an approximate crimping speed of 0.1 m/s. The setup comprised a rigid, nonsymmetric punch, a rigid anvil, an insulation crimp barrel of elastoplastic brass (CuZn 15), and an elastoplastic copper wire with polyethylene (PE) insulation. a Following common mathematical terminology in an explicit formulation the solutions, here the displacements u, are explicitly expressed as a function of other variables and parameters. In an implicit formulation the solution is contained within a function and its values have to be extracted, which is usually done through numerical algorithms. Copyright 2004 by Tyco Electronics Corporation. All rights reserved. 8 S. Kugener AMP Journal of Technology Vol. 4 June, 1995

2 ANSYS uses the following Newmark time integration method (see Ref. 1 for details): and are Newmark integration parameters and t is the time increment. This approach leads to the same type of inversion problem as before Figure 1. MIC I contact with insulation crimping as used in the study. The corresponding finite element model is 2-dimensional (2D), assuming plane stress. It is meshed with 680 2D, 4-node solid elements and 2400 contact elements with friction. Its parametric structure allows analysis of many types of crimps by setting a few, properly chosen geometrical parameters. the matrix (c 0 [M] + c l [c] + c 2 [K]) has to be inverted. Note that c... c 2 are constants. 0 Explicit formulation is based on the fundamental differential equation of dynamics OPTRIS uses central difference integration with THE MODEL Implicit and Explicit Formulations Chosen for the Model Implicit formulation allows a pseudo-static and a dynamic approach. In the pseudostatic case acceleration and velocity forces are neglected, the principle of virtual work leads to which yields with [K] = stiffness matrix, {u} = displacement pseudo vector, {F} = external force pseudo vector. The displacement is found directly by inversion of the rigidity matrix, which is a classical inversion problem. It uses Gauss elimination, with its usual matrix conditioning problems. In this case, computer time is proportional to the square of the mean wavefront size. For the implicit dynamic case acceleration and velocity forces are carried, it is with [M] = mass matrix, [c] = damping matrix. Because [M] is a diagonal mass matrix, the inversion is direct. The only difficulty is to ensure stable time integration by introducing the condition that the time increments be very small, that is L = some characteristic element dimension, = longitudinal velocity of sound, t = smallest transit time of dilatational waves to cross any element in the mesh. 3 The smallness of the time increment is an obvious disadvantage. If it is sufficiently small to be negligible for fast dynamic events, it is also highly penalizing for low speed phenomena such as most metal forming processes. However, the analysis can be accelerated by increasing AMP Journal of Technology Vol. 4 June, 1995 S. Kugener 9

3 L artificially the mass density of the material, b L the die velocity if inertia effects remain negligible. 4 For the crimping process considered here, the die velocity is about 0.1 m/s which is very slow for the default time increment. An increase of the speed to 1 m/s did not result in a noticeable inertia effect. Correction of the Stress Versus Strain Function A prime difficulty in applying finite element analysis to problems of the kind considered here is the lack of experimentally determined, true, nonlinear stress vs. strain curves. To be meaningful, finite element analysis needs as input ordered pairs of stress and strain at a given time t, i.e., the function t = true stress at time t, F = normal force, S t = true cross sectional area at time t, t = true strain at time t. Material manufacturers make available specific pairs of { ; }-values. Examples are 0.2% (the stress at 0.2% plastic strain as referred to the yield point), the tensile strength r, and the strain at rupture r. From these values one can construct a bilinear stress vs. strain relationship. With some experience and familiarity with the behavior of materials it is possible to smoothen that function and obtain a more realistic tensile test curve. From standardized tensile tests are available the uncorrected curves S 0 = initial cross sectional area, = (L - L 0 )/L 0, L 0 = initial length. This curve is reasonably accurate for small deformations, but is not acceptable for realistic, in-depth mechanical engineering analyses. However, such a curve can be corrected according to Francois et al. 5 by applying large deformation theory for incompressible materials and using the first-order approximations Then the corrected curves can be simulated by a multilinear model. This technique gave satisfactory results b An increase of two orders of magnitude provides a benefit of one order in cpu time. describing for brass contacts the behavior of the metal under consideration here. The case that included the behavior of the wires was more difficult to simulate. Additional considerations were required. First, the 19 strands of copper would lead to too large a model if internal sliding were taken into account. This phenomenon was modelled by introducing a copper core, thus artificially reducing the rigidity modulus of the copper. Furthermore, the behavior of the insulation became of critical importance, because one has to simulate simultaneously two effects. One of them is the large strain while the crimp voids are being filled, the other is the 20% reduction of the cross section under plane strain assumptions that would not allow expansion in the third dimension. Copper shows no noticeable compression under the conditions considered here. The behavior of the polyethylene insulation was simulated by assuming elastic deformation with a tensile modulus of 100 MPa, which is about one order of magnitude below that of most PE versions. Because of the incompressibility requirement, the plasticity codes could not be used. The problems associated with defining the behavior of the materials are the same for implicit and explicit analyses. Contact Element Management The crimping process is mainly a surface-on-surface sliding problem with a large number of elements sliding simultaneously on each other. The stability of the mechanism and to a certain extent also the final shape of the deformed parts is determined by the sliding friction taking place during crimping. This makes it important to include sliding friction in the models. The contact element management will be more difficult, because an initial equilibrium has to be found between the extrinsic force caused by the punch and the anvil reaction force, eventually through contact oscillation or rotation. ANSYS, in both the static and dynamic modes, simulates structure interactions with 3-node point to surface contact elements. Each node potentially in contact with an element has to be meshed, which penalizes the rigidity matrix size. Our model is composed of 2400 contact elements for only 680 structural elements. Contact interaction is solved in three steps: (1) (2) (3) A pinball algorithm fixes the close/open status of the gap. The next step determines whether penetration occurred or not. The third step applies the compatibility conditions if penetration occurs. A combination of penalization and Lagrange multiplier methods was used here, which in static code is the only one taking friction into account. The normal force, F n, is calculated by 10 S. Kugener AMP Journal of Technology Vol. 4 June, 1995

4 and i+l = Lagrange multiplier at iteration step (i + 1), K n = stiffness of the contact element, g = penetration depth, g = a user-defined penetration tolerance, p = an internally computed factor with p < 1. (14) [R] = orthogonal rotation matrix, [D R c ] = right stretch matrix. Then a strain tensor known as Hencky tensor d can be defined as [U] is the stretch matrix. Values of [ ] are numerically obtained through the spectral decomposition of [U] This algorithm is known for being more stable than the penalization method, and for avoiding conditioning problems 6,7 but still has convergence problems for the following reasons: The contact stiffness value is critical and can only be set iteratively. Different values have to be tuned along the punch to allow convergence. These values depend on the rigidity and plasticity of the local structure. Inadequately low stiffness values cause improper penetration and sometimes sticking. Inadequately high values affect substep size and convergence rate, and can generate major instabilities by introducing, in a single step, a resulting normal force exceeding the upper limit tolerated by the structure. In many cases this algorithm may cause an infinite number of numerical oscillations between several contacts. For explicit codes, the very small time increment permits verification of contact compatibility node by node and at each time step. This procedure is based on a contact search/tracking algorithm and a double loop master/slave algorithm which determines normal and tangential forces. 8 The formulation has the advantages of being very effective, almost transparent for the user, and will not affect the problem size. Large Strain/Deflection Formulation The common assumption of small displacements cannot be applied to crimping. Therefore, a Lagrangian formulation has to be used instead of the usual first-order tensor in the Eulerian formulation [ ij ] = first-order strain tensor, {u i } = displacement vector. In this case, a strain tensor is formed from the polar decomposition of the Lagrangian deformation gradient i are the eigenvalues and i are the eigenvectors of [U]. The Hencky formulation will not conserve orthogonality during deformation and may cause interaction problems with other algorithms using the assumption of orthotropic elasticity or kinetic stress hardening. For this reason, an isotropic stress hardening law for the brass was used. This kind of algorithm is powerful and relatively economical in computer time. However, it penalizes the analysis because it needs an eigenvalue extraction from the [U] matrix. Furthermore, it may simulate a realistic force/displacement curve but with a negative derivative, which may cause the Newton-Raphson e solution algorithm to crash. The latter instability could be avoided with a Riks arc length algorithm l0, which is not yet available to cover plasticity. In the dynamic analysis the instability will be dampened by inertia forces. Explicit codes do not need large strain/deformation formulation because the displacements are reactivated at each substep, very small time increments support small perturbation assumptions, the behavior under stress is described by true stress/ strain curves. Buckling and Bifurcation Instabilities Simulation of the crimping process is very delicate because of several physical instabilities. Some of them are caused by friction. If it is very high, the barrel may collapse on one side of the anvil. Transverse buckling of the barrel wall, towards the inside, seems to be numerically possible. The two barrel ends buttress together before one of them moves under the other. Flat ends emphasize this behavior. Because of the rolling mechanism of the barrel, the combined effects of friction and the mechanics itself can cause instability. When the barrel end begins to slide in the crimping radius, a large part of the vertical forces is supported by friction because of a long area of contact. But c The corresponding left turn expression would be [F] = [D L ] [R] with [D L ] being the left stretch matrix. d For details on the Hencky tensor see for instance reference 9. e For details see reference 1. AMP Journal of Technology Vol. 4 June, 1995 S. Kugener 11

5 when the barrel begins to roll, it decreases the contact area thus creating the condition for snap-through. In some cases, particularly when the wire is thin, fracture may occur at the bottom of the barrel. Only few explicit codes can handle fracture mechanics. As discussed above, large strain/deflection can partially simulate these phenomena with the Hencky tensor. With the implicit code, an initial buckling analysis would have made it possible to calculate buckling modes and critical loads. Post buckling analysis is not perfect for either implicit or explicit codes. Such formulations have to take into account defect rate and bifurcation equations which are not available in ANSYS and in OPTRIS. ll In static analysis, these phenomena can create major problems by generating infinite displacements or force/ displacement curves with a negative slope which usually crashes the Newton-Raphson algorithm of ANSYS. In all cases, dynamic analysis allows to pass through these instabilities because of the dampening effect of the inertia forces that also reduces noise in the results. COMPARISON OF EXPERIMENTAL AND THEORETICAL RESULTS In a first example the deformation of an idealized, simple, geometric model was simulated with ANSYS. An implicit, dynamic approach was used. The result is shown in Figure 2. Because of the long time required for the simulation, the method was not extended to model realistic, industrially important crimps. Figure 2. Strain distribution in a simplified, non-industrial crimp simulated by ANSYS. Results of explicit code simulations are shown in Figures 3 to 5 together with photographs taken on actual samples in the corresponding stages of the crimping process. These figures permit a comparison of deformations as predicted by the simulation with those experienced in real life. While the photos taken on experimental samples show only deformations, the graphs produced by the simulation show also the distribution of the strain. This feature can be used for targeted modifications of existing designs, material selection, and process parameters. Figure 3 permits comparison of deformations obtained at three different crimping stages with those of a prototype, small punch, and a thin wire of 1 mm 2 cross section. Shown are the strain contours at the beginning of crimping, the beginning of overlap during the process, and at its end. In all cases, the agreement between experimental results and those predicted by the explicit code simulation is excellent. The general phenomena are well described. Some differences appear in the curvature of the barrel ends. This can be understood by considering that the tests were static, which emphasizes buttress and buckling phenomena because of the higher coefficient of static friction. Under true dynamic conditions, this effect is smoothed by inertia forces and the lower coefficient of dynamic friction. The deflection of the left-handed barrel side is less pronounced. This effect was taken care of by introducing a larger coefficient of friction. In Figure 4, results of experimental and simulated crimps are compared for the initial large punch and a thin wire of 1 mm 2 cross section. In the photographs on the left side, the deformation of the experimental sample is shown. The photos on the right side present the results of the simulation, which include deformation and strain. Here only the two stages beginning and end are covered. Figure 5 comprises examples of experimentally observed and predicted deformations and shows the strain distribution revealed by the simulation. As in Figure 4, the results are from a setup utilizing the initial large punch but now with a thick wire of 3 mm 2. Shown are only the conditions existing at the beginning of crimping. In all simulations the contact rotated a little to the right. This phenomenon is not confirmed by the test results shown in the figures. The reason is that in the tests the contact position was corrected initially to place the wire. However, very frequently such a rotation was observed in other experimental studies. The final configuration of the case described in Figure 4 was simulated poorly because fracture was not included in the algorithms. The high strain values computed in the simulation indicate that such had happened. Finally, it is interesting to note that, as shown by the simulation, friction has a major effect on the crimping process. STRATEGIES FOR OPTIMAL SIMULATION Based on the previous discussions, it appears that an explicit code will offer the best approach mainly because it allows a direct solution without considering in advance problems with contact element management. In this case, no special strategy has to be developed. If necessary, only speed increase options should be considered. 12 S. Kugener AMP Journal of Technology Vol. 4 June, 1995

6 Figure 3. Comparison of results of experimental and simulated crimps at three different stages for prototype small punch and thin wire. The simulation used explicit code. Experimental results in the left, simulated results in the right hand column. Shown are the strain contours at (a) the beginning of crimping, (b) the beginning of overlap, (c) the end of crimping. An important result of the study presented in this paper is that implicit static analysis was inadequate for simulating the crimping process. The Newton-Raphson algorithm has convergence problems around the physical instabilities of the system modelled. Various numerical precautionary measures were tried, but proved to be useless. Dynamic, implicit analysis will give a solution but only through a difficult iterative process. If this approach is chosen, a rigorous strategy is recommended considering the following points: Smoothed, true stress/strain curves must be entered. Most refined algorithms are required in form of versions including large deflection, large strain and friction. Contact stiffness values must be set iteratively and individually. AMP Journal of Technology Vol. 4 June, 1995 S. Kugener 13

7 Figure 4. Comparison of results of experimental and simulated crimps at two different stages for the initial large punch and thin wire. The simulation used explicit code. Experimental results in the left, simulated results in the right-hand column. Shown are the strain contours at (a) the beginning of crimping, (b) the end of crimping. Figure 5. Comparison of results of experimental and simulated crimps at two different stages for the initial large punch and thick wire. The simulation used explicit code. Experimental results on the left, simulated results on the right. Shown are the strain contours at the beginning of crimping. A contact stiffness value of K n = E/100 is a good initial value. The penetration tolerance g should be as large as possible, taking into consideration that the penalization force F p = K n. The tolerance for convergence of the solution should be set to the minimum for a first successful run and subsequently increased for more accurate results. Convergence tools such as predictors, bisections and linesearch algorithms should be avoided as they may 14 S. Kugener AMP Journal of Technology Vol. 4 June, 1995

8 create more problems than they really solve. Varying the time increment is generally more effective. In some cases, a larger time increment helps circumvent a local difficulty. SUMMARY This study shows that a complex metal-forming process can be simulated accurately and that finite element methods have reached a degree of maturity to make them useful as design tools and for process analysis. Explicit codes have to be considered for the future. They are expected to become of particular importance for die optimization and process simulation. Before choosing a specific approach for the simulation, a thorough problem definition including a discussion of anticipated difficulties should be prepared. The latter is particularly important because choosing the wrong simulation tool can lead to undesirable, lengthy calculation procedures or even to a dead end. For crimping simulations an explicit code is superior to an implicit one, especially because of the large number of contact elements used in the mathematical model. ACKNOWLEDGMENTS The author would like to thank Jean-Pierre Picaud for his contributions in the area of parametric geometry and Richard Hughes whose input greatly supported and enhanced the static analysis. He also thanks the members of AT&T DATAID for their collaboration especially, Nicolas Wolik, Projects Director, Scientific Department, who gave permission to run the ANSYS dynamic analysis, and Fabrice Villeval, Analysis Engineer in charge of fast dynamic simulation at the Scientific Department, who piloted the OPTRIS explicit analyses based on AMP s initial parametric deck. REFERENCES 1. ANSYS Theoretical Manual (Swanson Analysis Systems, Inc., Houston, PA, 1992). 2. OPTRIS, Theoretical Manual (Dynamic Software, Les Ulis, 1991). 3. D.P. Flanagan and T. Belytschko, Eigenvalues and stable time steps for the uniform strain Hexahedron and Quadrilateral, ASME J. of Applied Mechanics, 51, (1984). N. Rebelo, J.C. Nagtegaal and L.M. Taylor, Industrial application of implicit and explicit finite element method to forming process, WAM, D. Francois, A. Pineau, A. Zaoui, Elasticity et Plasticité (Ed. Hermes, Paris, 1993). Fortin/Glowinski, Méthodes de Lagrangiens Augmentés, Col. Méthodes Mathématiques pour l informatique (Dunod, Paris, 1982). D.A. Pierre, M.J. Lowe, Mathematical Programming via Augmented Lagrangians (Addison Wesley, Reading, MA, 1975). K. Sweizerhoff, Metal forming simulation with explicit time integration, International LS-DYNA 3D conference, Birmingham, U. K., April G.G. Weber and A.N. Zava Plangos, An objective time integration procedure for isotropic rate-independent elastic-plastic constitutive equations, Internatl. J. of Plasticity, 6, (1990). New Features of ANSYS Rev. 5.0 A (Swanson Analysis Systems, Inc., Houston, PA, 1993). C. Gontier, Stabilité des structures (Ed. Ecole Centrale, Paris, 1989). Stéphane Kugener is Advanced Projects Engineer in the Engineering Department of AMP de France. Stéphane Kugener received an Engineer Diploma in Structural Mechanics from Ecole Speciale des Travaux Publics of Paris in He specialized in Solid Mechanics and Scientific Calculation through a 3rd year at Ecole Centrale de Paris in From 1990 to 1991 he worked as civil engineer on seismic verification of bridges and dams in the engineering department of SOCOTEC. Since joining AMP in 1992, he supervised CAE activities and developed use of numerical simulation for injection molding, mechanical and thermal electric analyses. He is a member of AMP s European Moldflow Working Group. AMP Journal of Technology Vol. 4 June, 1995 S. Kugener 15

Finite Element Analysis of Compression of Thin, High Modulus, Cylindrical Shells with Low-Modulus Core

Finite Element Analysis of Compression of Thin, High Modulus, Cylindrical Shells with Low-Modulus Core Finite Element Analysis of Compression of Thin, High Modulus, Cylindrical Shells with Low-Modulus Core Robert S. Joseph Design Engineering Analysis Corporation, McMurray, PA ABSTRACT Long, cylindrical

More information

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

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

More information

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

Nonlinear analysis in ADINA Structures

Nonlinear analysis in ADINA Structures Nonlinear analysis in ADINA Structures Theodore Sussman, Ph.D. ADINA R&D, Inc, 2016 1 Topics presented Types of nonlinearities Materially nonlinear only Geometrically nonlinear analysis Deformation-dependent

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

On Nonlinear Buckling and Collapse Analysis using Riks Method

On Nonlinear Buckling and Collapse Analysis using Riks Method Visit the SIMULIA Resource Center for more customer examples. On Nonlinear Buckling and Collapse Analysis using Riks Method Mingxin Zhao, Ph.D. UOP, A Honeywell Company, 50 East Algonquin Road, Des Plaines,

More information

Nonlinear Buckling Prediction in ANSYS. August 2009

Nonlinear Buckling Prediction in ANSYS. August 2009 Nonlinear Buckling Prediction in ANSYS August 2009 Buckling Overview Prediction of buckling of engineering structures is a challenging problem for several reasons: A real structure contains imperfections

More information

Influence of residual stresses in the structural behavior of. tubular columns and arches. Nuno Rocha Cima Gomes

Influence of residual stresses in the structural behavior of. tubular columns and arches. Nuno Rocha Cima Gomes October 2014 Influence of residual stresses in the structural behavior of Abstract tubular columns and arches Nuno Rocha Cima Gomes Instituto Superior Técnico, Universidade de Lisboa, Portugal Contact:

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

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

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

INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS

INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS INTRODUCTION TO THE EXPLICIT FINITE ELEMENT METHOD FOR NONLINEAR TRANSIENT DYNAMICS SHEN R. WU and LEI GU WILEY A JOHN WILEY & SONS, INC., PUBLICATION ! PREFACE xv PARTI FUNDAMENTALS 1 1 INTRODUCTION 3

More information

Micro-meso draping modelling of non-crimp fabrics

Micro-meso draping modelling of non-crimp fabrics Micro-meso draping modelling of non-crimp fabrics Oleksandr Vorobiov 1, Dr. Th. Bischoff 1, Dr. A. Tulke 1 1 FTA Forschungsgesellschaft für Textiltechnik mbh 1 Introduction Non-crimp fabrics (NCFs) are

More information

The Finite Element Method for Mechonics of Solids with ANSYS Applicotions

The Finite Element Method for Mechonics of Solids with ANSYS Applicotions The Finite Element Method for Mechonics of Solids with ANSYS Applicotions ELLIS H. DILL 0~~F~~~~"P Boca Raton London New Vork CRC Press is an imprint 01 the Taylor & Francis Group, an Informa business

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

DYNAMIC RESPONSE OF THIN-WALLED GIRDERS SUBJECTED TO COMBINED LOAD

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

More information

2D Liquefaction Analysis for Bridge Abutment

2D Liquefaction Analysis for Bridge Abutment D Liquefaction Analysis for Bridge Abutment Tutorial by Angel Francisco Martinez Integrated Solver Optimized for the next generation 64-bit platform Finite Element Solutions for Geotechnical Engineering

More information

Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases

Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Georg Haasemann Conrad Kloß 1 AIMCAL Conference 2016 MOTIVATION Wrinkles in web handling system Loss of

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

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies

Soft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed

More information

Hydroplaning Simulation using MSC.Dytran

Hydroplaning Simulation using MSC.Dytran Hydroplaning Simulation using MSC.Dytran Toshihiko Okano * & Masataka Koishi * THE YOKOHAMA RUBBER CO., LTD 2-1 Oiwake Hiratsuka Kanagawa 254-8601, Japan ABSTRACT Hydroplaning characteristics is one of

More information

FEA A Guide to Good Practice. What to expect when you re expecting FEA A guide to good practice

FEA A Guide to Good Practice. What to expect when you re expecting FEA A guide to good practice FEA A Guide to Good Practice What to expect when you re expecting FEA A guide to good practice 1. Background Finite Element Analysis (FEA) has transformed design procedures for engineers. Allowing more

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

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

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

More information

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

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

Theoretical Manual Theoretical background to the Strand7 finite element analysis system

Theoretical Manual Theoretical background to the Strand7 finite element analysis system Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter

More information

Size Effects In the Crushing of Honeycomb Structures

Size Effects In the Crushing of Honeycomb Structures 45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference 19-22 April 2004, Palm Springs, California AIAA 2004-1640 Size Effects In the Crushing of Honeycomb Structures Erik C.

More information

A Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against a Rigid Flat under Varying Modulus of Elasticity and Sphere Radius

A Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against a Rigid Flat under Varying Modulus of Elasticity and Sphere Radius Engineering, 2010, 2, 205-211 doi:10.4236/eng.2010.24030 Published Online April 2010 (http://www. SciRP.org/journal/eng) 205 A Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against

More information

Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder. Richard Kennison, Race-Tec

Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder. Richard Kennison, Race-Tec Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder Richard Kennison, Race-Tec Agenda Race-Tec Overview Accumulator Experimental Testing Material Testing Numerical Analysis: 1. Linear Buckling

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

Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure

Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure ISSN : 48-96, Vol. 6, Issue 8, ( Part -4 August 06, pp.3-38 RESEARCH ARTICLE Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure Dr.Ragbe.M.Abdusslam Eng. Khaled.S.Bagar ABSTRACT

More information

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

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

More information

Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity

Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity Lecture #6: 3D Rate-independent Plasticity (cont.) Pressure-dependent plasticity by Borja Erice and Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling

More information

An overview of Carbon Fiber modeling in LS-DYNA. John Zhao October 23 th 2017

An overview of Carbon Fiber modeling in LS-DYNA. John Zhao October 23 th 2017 An overview of Carbon Fiber modeling in LS-DYNA John Zhao zhao@lstc.com October 23 th 2017 Outline Manufacturing of Carbon Fiber Compression molding *MAT_277 & 278 *MAT_293 *MAT_249 Resin transform molding

More information

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

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

More information

Heat Transfer Modeling using ANSYS FLUENT

Heat Transfer Modeling using ANSYS FLUENT Lecture 2 - Conduction Heat Transfer 14.5 Release Heat Transfer Modeling using ANSYS FLUENT 2013 ANSYS, Inc. March 28, 2013 1 Release 14.5 Agenda Introduction Energy equation in solids Equation solved

More information

IMPACT PROPERTY AND POST IMPACT VIBRATION BETWEEN TWO IDENTICAL SPHERES

IMPACT PROPERTY AND POST IMPACT VIBRATION BETWEEN TWO IDENTICAL SPHERES ICSV4 Cairns Australia 9- July, 7 IMPACT PROPERTY AND POST IMPACT VIBRATION BETWEEN TWO IDENTICAL SPHERES Hirofumi MINAMOTO, Keisuke SAITOH and Shozo KAWAMURA Toyohashi University of Technology Dept. of

More information

*MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials

*MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials 14 th International LS-DYNA Users Conference Session: Constitutive Modeling *MAT_PAPER and *MAT_COHESIVE_PAPER: Two New Models for Paperboard Materials Jesper Karlsson 1, Mikael Schill 1, Johan Tryding

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

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

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

More information

6 th Pipeline Technology Conference 2011

6 th Pipeline Technology Conference 2011 6 th Pipeline Technology Conference 2011 Mechanical behavior of metal loss in heat affected zone of welded joints in low carbon steel pipes. M. J. Fernández C. 1, J. L. González V. 2, G. Jarvio C. 3, J.

More information

Chapter 5. Vibration Analysis. Workbench - Mechanical Introduction ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.

Chapter 5. Vibration Analysis. Workbench - Mechanical Introduction ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved. Workbench - Mechanical Introduction 12.0 Chapter 5 Vibration Analysis 5-1 Chapter Overview In this chapter, performing free vibration analyses in Simulation will be covered. In Simulation, performing a

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

Dynamic Analysis Contents - 1

Dynamic Analysis Contents - 1 Dynamic Analysis Contents - 1 TABLE OF CONTENTS 1 DYNAMIC ANALYSIS 1.1 Overview... 1-1 1.2 Relation to Equivalent-Linear Methods... 1-2 1.2.1 Characteristics of the Equivalent-Linear Method... 1-2 1.2.2

More information

Comparison of LS-DYNA and NISA in Solving Dynamic Pulse Buckling Problems in Laminated Composite Beams

Comparison of LS-DYNA and NISA in Solving Dynamic Pulse Buckling Problems in Laminated Composite Beams 9 th International LS-DYNA Users Conference Simulation Technology (1) Comparison of LS-DYNA and NISA in Solving Dynamic Pulse Buckling Problems in Laminated Composite Beams Haipeng Han and Farid Taheri

More information

Seismic Pushover Analysis Using AASHTO Guide Specifications for LRFD Seismic Bridge Design

Seismic Pushover Analysis Using AASHTO Guide Specifications for LRFD Seismic Bridge Design Seismic Pushover Analysis Using AASHTO Guide Specifications for LRFD Seismic Bridge Design Elmer E. Marx, Alaska Department of Transportation and Public Facilities Michael Keever, California Department

More information

Conservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body

Conservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference

More information

Particle Blast Method (PBM) for the Simulation of Blast Loading

Particle Blast Method (PBM) for the Simulation of Blast Loading Particle Blast Method (PBM) for the Simulation of Blast Loading Hailong Teng, Jason Wang Livermore Software Technology Corporation Abstract This paper presents a particle blast method (PBM) to describe

More information

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

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

More information

A Simple Weak-Field Coupling Benchmark Test of the Electromagnetic-Thermal-Structural Solution Capabilities of LS-DYNA Using Parallel Current Wires

A Simple Weak-Field Coupling Benchmark Test of the Electromagnetic-Thermal-Structural Solution Capabilities of LS-DYNA Using Parallel Current Wires 13 th International LS-DYNA Users Conference Session: Electromagnetic A Simple Weak-Field Coupling Benchmark Test of the Electromagnetic-Thermal-Structural Solution Capabilities of LS-DYNA Using Parallel

More information

MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4

MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 MITOCW MITRES2_002S10nonlinear_lec15_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources

More information

3. Overview of MSC/NASTRAN

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

More information

Peak Strain and Displacement Sensors for Structural Health Monitoring

Peak Strain and Displacement Sensors for Structural Health Monitoring Peak Strain and Displacement Sensors for Structural Health Monitoring AKIRA MITA and SHINPEI TAKAHIRA ABSTRACT Simple and inexpensive passive sensors that can monitor the peak strain or displacement of

More information

COMPARISON BETWEEN 2D AND 3D ANALYSES OF SEISMIC STABILITY OF DETACHED BLOCKS IN AN ARCH DAM

COMPARISON BETWEEN 2D AND 3D ANALYSES OF SEISMIC STABILITY OF DETACHED BLOCKS IN AN ARCH DAM COMPARISON BETWEEN 2D AND 3D ANALYSES OF SEISMIC STABILITY OF DETACHED BLOCKS IN AN ARCH DAM Sujan MALLA 1 ABSTRACT The seismic safety of the 147 m high Gigerwald arch dam in Switzerland was assessed for

More information

LINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL

LINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL LINEAR AND NONLINEAR BUCKLING ANALYSIS OF STIFFENED CYLINDRICAL SUBMARINE HULL SREELATHA P.R * M.Tech. Student, Computer Aided Structural Engineering, M A College of Engineering, Kothamangalam 686 666,

More information

Probabilistic Assessment of a Stiffened Carbon Fibre Composite Panel Operating in its Postbuckled Region Authors: Correspondence:

Probabilistic Assessment of a Stiffened Carbon Fibre Composite Panel Operating in its Postbuckled Region Authors: Correspondence: Probabilistic Assessment of a Stiffened Carbon Fibre Composite Panel Operating in its Postbuckled Region Authors: D. Elder and R. Thomson Cooperative Research Centre for Advanced Composite Structures Limited

More information

Numerical simulation the bottom structures. grounding test by LS-DYNA

Numerical simulation the bottom structures. grounding test by LS-DYNA 5 th European LS-DYNA Users Conference Methods and Techniques (3) Numerical simulation the bottom structures grounding test by LS-DYNA Ainian Zhang Graduate School of Frontier Sciences, The University

More information

A METHOD OF LOAD INCREMENTS FOR THE DETERMINATION OF SECOND-ORDER LIMIT LOAD AND COLLAPSE SAFETY OF REINFORCED CONCRETE FRAMED STRUCTURES

A METHOD OF LOAD INCREMENTS FOR THE DETERMINATION OF SECOND-ORDER LIMIT LOAD AND COLLAPSE SAFETY OF REINFORCED CONCRETE FRAMED STRUCTURES A METHOD OF LOAD INCREMENTS FOR THE DETERMINATION OF SECOND-ORDER LIMIT LOAD AND COLLAPSE SAFETY OF REINFORCED CONCRETE FRAMED STRUCTURES Konuralp Girgin (Ph.D. Thesis, Institute of Science and Technology,

More information

A FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR AGAINST A RIGID FLAT UNDER VARYING MODULUS OF ELASTICITY AND SPHERE RADIUS

A FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR AGAINST A RIGID FLAT UNDER VARYING MODULUS OF ELASTICITY AND SPHERE RADIUS Proceedings of the International Conference on Mechanical Engineering 2009 (ICME2009) 26-28 December 2009, Dhaka, Bangladesh ICME09- A FINITE ELEMENT STUDY OF ELASTIC-PLASTIC HEMISPHERICAL CONTACT BEHAVIOR

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

Dynamic Response Of Laminated Composite Shells Subjected To Impulsive Loads

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

More information

Using Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses

Using Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses Using Energy History Data to Obtain Load vs. Deflection Curves from Quasi-Static Abaqus/Explicit Analyses Brian Baillargeon, Ramesh Marrey, Randy Grishaber 1, and David B. Woyak 2 1 Cordis Corporation,

More information

Example-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium

Example-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium Example-3 Title Cylindrical Hole in an Infinite Mohr-Coulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic

More information

SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED

SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED 1 SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED WITH LATERAL FRICTIONAL EFFECTS High slope results and comparison with experimental data Nathan Martin1, Ioan Ionescu2, Anne Mangeney1,3 François

More information

Numerical simulation of sheet metal forming processes using a new yield criterion

Numerical simulation of sheet metal forming processes using a new yield criterion Key Engineering Materials Vol. 344 (007) pp. 833-840 online at http://www.scientific.net (007) Trans Tech Publications, Switzerland Numerical simulation of sheet metal forming processes using a new yield

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

Modelling of ductile failure in metal forming

Modelling of ductile failure in metal forming Modelling of ductile failure in metal forming H.H. Wisselink, J. Huetink Materials Innovation Institute (M2i) / University of Twente, Enschede, The Netherlands Summary: Damage and fracture are important

More information

Outline. Tension test of stainless steel type 301. Forming including the TRIP-effect

Outline. Tension test of stainless steel type 301. Forming including the TRIP-effect Outline Erweiterte Kontaktformulierung, User-Reibmodell und thermische Analyse Simulationsmöglichkeiten zur Werkzeugauslegung mit LS-DYNA Workshop Simulation in der Umformtechnik Institut für Statik und

More information

PROGRAMMING THE TRANSIENT EXPLICIT FINITE ELEMENT ANALYSIS WITH MATLAB

PROGRAMMING THE TRANSIENT EXPLICIT FINITE ELEMENT ANALYSIS WITH MATLAB U.P.B. Sci. Bull., Series D, Vol. 75, Iss. 2, 2013 ISSN 1454-2358 PROGRAMMING THE TRANSIENT EXPLICIT FINITE ELEMENT ANALYSIS WITH MATLAB Andrei Dragoş Mircea SÎRBU 1, László FARKAS 2 Modern research in

More information

Discontinuous Galerkin methods for nonlinear elasticity

Discontinuous Galerkin methods for nonlinear elasticity Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity

More information

CHAPTER 3 THE EFFECTS OF FORCES ON MATERIALS

CHAPTER 3 THE EFFECTS OF FORCES ON MATERIALS CHAPTER THE EFFECTS OF FORCES ON MATERIALS EXERCISE 1, Page 50 1. A rectangular bar having a cross-sectional area of 80 mm has a tensile force of 0 kn applied to it. Determine the stress in the bar. Stress

More information

Ultimate shear strength of FPSO stiffened panels after supply vessel collision

Ultimate shear strength of FPSO stiffened panels after supply vessel collision Ultimate shear strength of FPSO stiffened panels after supply vessel collision Nicolau Antonio dos Santos Rizzo PETROBRAS Rio de Janeiro Brazil Marcelo Caire SINTEF do Brasil Rio de Janeiro Brazil Carlos

More information

DESCRIBING THE PLASTIC DEFORMATION OF ALUMINUM SOFTBALL BATS

DESCRIBING THE PLASTIC DEFORMATION OF ALUMINUM SOFTBALL BATS DESCRIBING THE PLASTIC DEFORMATION OF ALUMINUM SOFTBALL BATS E. BIESEN 1 AND L. V. SMITH 2 Washington State University, 201 Sloan, Spokane St, Pullman, WA 99164-2920 USA 1 E-mail: ebiesen@gonzaga.edu 2

More information

Project. First Saved Monday, June 27, 2011 Last Saved Wednesday, June 29, 2011 Product Version 13.0 Release

Project. First Saved Monday, June 27, 2011 Last Saved Wednesday, June 29, 2011 Product Version 13.0 Release Project First Saved Monday, June 27, 2011 Last Saved Wednesday, June 29, 2011 Product Version 13.0 Release Contents Units Model (A4, B4) o Geometry! Solid Bodies! Parts! Parts! Body Groups! Parts! Parts

More information

ME 243. Mechanics of Solids

ME 243. Mechanics of Solids ME 243 Mechanics of Solids Lecture 2: Stress and Strain Ahmad Shahedi Shakil Lecturer, Dept. of Mechanical Engg, BUET E-mail: sshakil@me.buet.ac.bd, shakil6791@gmail.com Website: teacher.buet.ac.bd/sshakil

More information

CAEFEM v9.5 Information

CAEFEM v9.5 Information CAEFEM v9.5 Information Concurrent Analysis Corporation, 50 Via Ricardo, Thousand Oaks, CA 91320 USA Tel. (805) 375 1060, Fax (805) 375 1061 email: info@caefem.com or support@caefem.com Web: http://www.caefem.com

More information

Nonlinear Analysis of Reinforced Concrete Shells Subjected to Impact Loads

Nonlinear Analysis of Reinforced Concrete Shells Subjected to Impact Loads Transactions of the 7 th International Conference on Structural Mechanics in Reactor Technology (SMiRT 7) Prague, Czech Republic, August 7, 00 Paper # J0- Nonlinear Analysis of Reinforced Concrete Shells

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

A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials

A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical

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

GATE SOLUTIONS E N G I N E E R I N G

GATE SOLUTIONS E N G I N E E R I N G GATE SOLUTIONS C I V I L E N G I N E E R I N G From (1987-018) Office : F-16, (Lower Basement), Katwaria Sarai, New Delhi-110016 Phone : 011-65064 Mobile : 81309090, 9711853908 E-mail: info@iesmasterpublications.com,

More information

870. Vibrational analysis of planetary gear trains by finite element method

870. Vibrational analysis of planetary gear trains by finite element method 870. Vibrational analysis of planetary gear trains by finite element method Pei-Yu Wang 1, Xuan-Long Cai 2 Department of Mechanical Design Engineering, National Formosa University Yun-Lin County, 632,

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

NONLINEAR ANALYSIS OF A DAM-RESERVOIR-FOUNDATION SYSTEM UNDER SPATIALLY VARIABLE SEISMIC EXCITATIONS

NONLINEAR ANALYSIS OF A DAM-RESERVOIR-FOUNDATION SYSTEM UNDER SPATIALLY VARIABLE SEISMIC EXCITATIONS October 12-17, 28, Beijing, China NONLINEAR ANALYSIS OF A DAM-RESERVOIR-FOUNDATION SYSTEM UNDER SPATIALLY VARIABLE SEISMIC EXCITATIONS J. Huang 1 and A. Zerva 2 1 Dept. of Civil, Architectural and Environmental

More information

ACCESS Beiratssitzung November Dr. G. Laschet and L. Haas ACCESS e.v. Aachen, Germany. FENET Workshop Copenhagen

ACCESS Beiratssitzung November Dr. G. Laschet and L. Haas ACCESS e.v. Aachen, Germany. FENET Workshop Copenhagen Coupling FE Contact and Heat Transfer Analysis in Casting Simulations ACCESS Beiratssitzung November 2001 Dr. G. Laschet and L. Haas ACCESS e.v. Aachen, Germany FENET Workshop Copenhagen 27.02.2002 About

More information

PLAT DAN CANGKANG (TKS 4219)

PLAT DAN CANGKANG (TKS 4219) PLAT DAN CANGKANG (TKS 4219) SESI I: PLATES Dr.Eng. Achfas Zacoeb Dept. of Civil Engineering Brawijaya University INTRODUCTION Plates are straight, plane, two-dimensional structural components of which

More information

Flexure: Behavior and Nominal Strength of Beam Sections

Flexure: Behavior and Nominal Strength of Beam Sections 4 5000 4000 (increased d ) (increased f (increased A s or f y ) c or b) Flexure: Behavior and Nominal Strength of Beam Sections Moment (kip-in.) 3000 2000 1000 0 0 (basic) (A s 0.5A s ) 0.0005 0.001 0.0015

More information

Members Subjected to Torsional Loads

Members Subjected to Torsional Loads Members Subjected to Torsional Loads Torsion of circular shafts Definition of Torsion: Consider a shaft rigidly clamped at one end and twisted at the other end by a torque T = F.d applied in a plane perpendicular

More information

Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation

Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation P. Sushma Ph D Scholar, Earthquake Engineering Research Center, IIIT Hyderabad, Gachbowli,

More information

SIZE EFFECTS IN THE COMPRESSIVE CRUSHING OF HONEYCOMBS

SIZE EFFECTS IN THE COMPRESSIVE CRUSHING OF HONEYCOMBS 43rd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Con 22-25 April 2002, Denver, Colorado SIZE EFFECTS IN THE COMPRESSIVE CRUSHING OF HONEYCOMBS Erik C. Mellquistand Anthony M.

More information

FINITE ELEMENT ANALYSIS OF TAPERED COMPOSITE PLATE GIRDER WITH A NON-LINEAR VARYING WEB DEPTH

FINITE ELEMENT ANALYSIS OF TAPERED COMPOSITE PLATE GIRDER WITH A NON-LINEAR VARYING WEB DEPTH Journal of Engineering Science and Technology Vol. 12, No. 11 (2017) 2839-2854 School of Engineering, Taylor s University FINITE ELEMENT ANALYSIS OF TAPERED COMPOSITE PLATE GIRDER WITH A NON-LINEAR VARYING

More information

International Journal of Advanced Engineering Technology E-ISSN

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

More information

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

Evaluation of dynamic behavior of culverts and embankments through centrifuge model tests and a numerical analysis

Evaluation of dynamic behavior of culverts and embankments through centrifuge model tests and a numerical analysis Computer Methods and Recent Advances in Geomechanics Oka, Murakami, Uzuoka & Kimoto (Eds.) 2015 Taylor & Francis Group, London, ISBN 978-1-138-00148-0 Evaluation of dynamic behavior of culverts and embankments

More information

Tolerance Ring Improvement for Reducing Metal Scratch

Tolerance Ring Improvement for Reducing Metal Scratch International Journal of Scientific and Research Publications, Volume 2, Issue 11, November 2012 1 Tolerance Ring Improvement for Reducing Metal Scratch Pattaraweerin Woraratsoontorn*, Pitikhate Sooraksa**

More information

ALASKA ENERGY AUTHORITY AEA ENGINEERING FEASIBILITY REPORT. Appendix B8. Finite Element Analysis

ALASKA ENERGY AUTHORITY AEA ENGINEERING FEASIBILITY REPORT. Appendix B8. Finite Element Analysis ALASKA ENERGY AUTHORITY AEA11-022 ENGINEERING FEASIBILITY REPORT Appendix B8 Finite Element Analysis Susitna-Watana Hydroelectric Project Alaska Energy Authority FERC Project No. 14241 December 2014 Seismic

More information

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics

More information

Back Calculation of Rock Mass Modulus using Finite Element Code (COMSOL)

Back Calculation of Rock Mass Modulus using Finite Element Code (COMSOL) Back Calculation of Rock Mass Modulus using Finite Element Code (COMSOL) Amirreza Ghasemi 1. Introduction Deformability is recognized as one of the most important parameters governing the behavior of rock

More information

STRESS, STRAIN AND DEFORMATION OF SOLIDS

STRESS, STRAIN AND DEFORMATION OF SOLIDS VELAMMAL COLLEGE OF ENGINEERING AND TECHNOLOGY, MADURAI 625009 DEPARTMENT OF CIVIL ENGINEERING CE8301 STRENGTH OF MATERIALS I -------------------------------------------------------------------------------------------------------------------------------

More information