Finite Strain Elastoplasticity: Consistent Eulerian and Lagrangian Approaches

Total Page:16

File Type:pdf, Size:1020Kb

Finite Strain Elastoplasticity: Consistent Eulerian and Lagrangian Approaches Finite Strain Elastoplasticity: Consistent Eulerian and Lagrangian Approaches by Mohammad Amin Eshraghi A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Doctor of Philosophy in Mechanical Engineering Waterloo, Ontario, Canada, 2009 © Mohammad Amin Eshraghi 2009 Author’s Declaration I hereby declare that I am the sole author of this thesis. This is a true copy of the thesis, including any required final revisions, as accepted by my examiners. I understand that my thesis may be made electronically available to the public. Mohammad Amin Eshraghi ii Abstract Infinitesimal strain approximation and its additive decomposition into elastic and plastic parts used in phenomenological plasticity models are incapable of predicting the hardening behavior of materials for large strain loading paths. Experimentally observed second-order effect in finite torsional loading of cylindrical bars, known as the Swift effect, as well as deformations involving significant amount of rotations are examples for which infinitesimal models fail to predict the material response accurately. Several different Eulerian and Lagrangian formulations for finite strain elastoplasticity have been proposed based on different decompositions of deformation and their corresponding flow rules. However, issues such as spurious shear oscillation in finite simple shear and elastic dissipation in closed-path loadings as well as elastic ratchetting under cyclic loading have been identified with the classical formulations for finite strain analysis. A unified framework of Eulerian rate-type constitutive models for large strain elastoplasticity is developed here which assigns no preference to the choice of objective corotational rates. A general additive decomposition of arbitrary corotational rate of the Eulerian strain tensor is proposed. Integrability of the model for the elastic part of the deformation is investigated and it is shown that the proposed unified model is consistent with the notion of hyperelasticity for its elastic part. Based on this, the stress power is physically separable into its reversible and irreversible parts using the proposed constitutive model irrespective of the objective rate used in the model. As a result, all of the issues of finite strain elastoplasticity are resolved using the proposed Eulerian rate model for arbitrary corotational rate of stress. A modified multiplicative decomposition of the right stretch tensor is proposed and used to set up a new Lagrangian framework for finite strain elastoplasticity. Decomposition of the deformation is solely defined by the iii multiplicative decomposition of the total right stretch tensor into its elastic and plastic parts. The flow rule and evolution of the plastic internal variables are based on the Hencky measure of the plastic right stretch tensor instead of the strain rate tensor. As a result, the issue of mismatch between the elastic and plastic parts of the deformation which mostly exists in the classical multiplicative models does not exist in the proposed Lagrangian model. The problem of back stress oscillation observed in the classical Lagrangian models is also resolved using the proposed Lagrangian model and results are identical to those of the proposed unified Eulerian rate model for finite strain elastoplasticity. In the context of nonlinear elasticity, no preference for either Lagrangian or Eulerian formulations exists since the two formulations can be related through proper transformations and are equivalent form of each other in different backgrounds. However, classical Eulerian and Lagrangian models of elastoplasticity do not provide such an equivalency under the same loading path. This is due to different definitions used for the elastic and plastic parts of the deformation and different flow rules used in the classical Eulerian and Lagrangian models. In this research it is shown that both the proposed Lagrangian and unified Eulerian rate models are equivalent and results obtained from both models are identical for the same finite strain loading path. Such an equivalency verifies that the proposed Eulerian and Lagrangian models are unified and transformable to each other. The unified Eulerian and Lagrangian models are extended to mixed nonlinear hardening material behavior. Predicted results for the second-order effect (the well- known Swift effect) are in good agreement with experimental data for fixed-end finite torsional loading of SUS 304 stainless steel tubes. The proposed models are therefore good candidates to be implemented in the displacement-based formulation of the finite element method for the Lagrangian and Eulerian frameworks of finite strain elastoplasticity. iv Acknowledgments I would like to express my sincere gratitude to my supervisors, Professor Hamid Jahed and Professor Steve Lambert for their guidance, technical support, and helpful discussions during the course of this research. Financial support from the Natural Sciences and Engineering Research Council of Canada (NSERC) is gratefully acknowledged. I would also like to express my gratitude to Professor Behrooz Farshi for his guidance and encouragement during my MASc program at Iran University of Science and Technology. The knowledge and efforts of my colleagues in our group was a valuable source of inspiration and success. In particular, Mohammad Noban, Morvarid Karimi Ghovanlou, Arash Tajik, and Jafar Al Bin Mousa, are thanked for the great discussions we had. Many thanks to my friends Amir Poursaee, Amir Noroozi, Semyon Mikheevskiy, Ramtin Movassaghi, Hamidreza Alemohammad and other wonderful friends. Their friendship has made the past four years full of joy and pleasure. Finally, I would like to express my profound gratitude to my parents for their support and encouragement throughout my education and professional career. v Table of Contents List of Figures .............................................................................................................x List of Tables ......................................................................................................... xiii Nomenclature ......................................................................................................... xiv Chapter 1 Introduction ........................................................................................1 1.1 Background ..........................................................................................6 1.2 Finite deformation plasticity models ...................................................8 1.3 Objectives and outline of the thesis ...................................................14 Chapter 2 Review of Elastic constitutive models for finite deformations .....18 2.1 General principles ..............................................................................19 2.2 Kinematics for finite deformation ......................................................20 2.3 Tensor transformation, objectivity, and objective rates .....................23 2.3.1 Objective corotational rates................................................... 24 2.3.2 Convected rates ..................................................................... 28 2.4 Stress measures ..................................................................................31 2.5 Work conjugacy .................................................................................33 2.5.1 Hill’s original work conjugacy ............................................. 33 2.5.2 Unified work conjugacy ........................................................ 34 2.6 Finite Elasticity and Hypoelastic material models.............................38 2.6.1 Simple materials and Cauchy elasticity ................................ 39 vi 2.6.2 Green elasticity...................................................................... 43 2.6.3 Hypoelasticity ....................................................................... 45 2.6.4 Issues with hypoelasticity ..................................................... 46 2.6.5 Logarithmic (D) rate of stress and integrability conditions .. 60 2.6.6 Unconditional integrability of the grade-zero hypoelastic model based on the logarithmic (D) spin .......................................... 65 2.7 Summary ............................................................................................66 Chapter 3 Review of Eulerian rate plasticity models and their algorithmic implementation ...........................................................................................................67 3.1 Thermodynamics of irreversible deformations ..................................68 3.2 Classical infinitesimal plasticity ........................................................73 3.2.1 A general quadratic form ...................................................... 73 3.2.2 Principle of maximum plastic dissipation ............................. 74 3.3 Algorithmic implementation of infinitesimal plasticity .....................76 3.3.1 Closest point projection method ........................................... 77 3.3.2 Plastic corrector step using the return mapping algorithm.... 80 3.4 Extension of the infinitesimal plasticity models to finite deformation based on hypoelastic material models .........................................................85 3.4.1 Hypo-based finite plasticity models ...................................... 86 3.4.2 Prager’s yielding stationary condition and choice of objective rates .............................................................................................
Recommended publications
  • 1 Useful Definitions Or Concepts
    ETH Zurich Institute for Mechanical Systems Department of Mechanical and Process Engineering Center of Mechanics Winter 06/07 Nonlinear Continuum Mechanics Exercise 9 Prof. Dr. Sanjay Govindjee 1 Useful Definitions or Concepts 1.1 Elastic constitutive laws One general type of elastic material model is the one called Cauchy elastic material, which depend on only the current local deformation of the material and has no dependence on past history. This type is often to general and a slightly more restrictive model called Green elastic materials which are derivable from a scalar function (strain energy density function) are often used in mechanics. • Cauchy elastic S = Sˆ(deformation measure) ∂Ψ • Green elastic S = ∂deformation measure) The definition of Green elastic materials coincide with that of hyperelasticity. From hyperelasticity and thermody- namic consistency, one obtains, ∂Ψ P = (1) ∂F or ∂Ψ S = 2 . (2) ∂C Objectivity requires that Ψˆ depend only the stretches or for example on the Right Cauchy deformation tensor, Ψ(ˆ C) . (3) If the material is isotropic, then the Ψ can also be a function of the Left Cauchy deformation tensor, Ψ(ˆ b) (4) or equivalently, Ψ(ˆ I1(b),I2(b),I3(b)) (5) Ψ(ˆ I1(C),I2(C),I3(C)) . (6) 1 ETH Zurich Institute for Mechanical Systems Department of Mechanical and Process Engineering Center of Mechanics Winter 06/07 Nonlinear Continuum Mechanics Exercise 9 Prof. Dr. Sanjay Govindjee 2 Applications Problem: Consider the Neo-Hookean model, 2 Ψ = c1(I1 − 3) + c2(J − 1) + c3lnJ. (7) What should c1, c2, c3 be such that this model is consistent with linear elasticity? Express ci (i = 1, 2, 3) in terms of the Lame constants.
    [Show full text]
  • Supplementary Notes For: Nonlinear Elasticity with Applications to the Arterial Wall Short Course in the EU-Regional School in C
    A. M. Robertson, Nonlinear elasticity, Spring 2009 1 Supplementary Notes for: Nonlinear elasticity with applications to the arterial wall Short Course in the EU-Regional School in Computational Engineering Science RWTH Aachen University Part I- March 27, 2009 Part II- April 3, 2009 Associate Professor A. M. Robertson Department of Mechanical Engineering and Materials Science McGowan Institute for Regenerative Medicine Center for Vascular Remodeling and Regeneration (CVRR) University of Pittsburgh 2 A. M. Robertson, Supplementary Notes Contents 1 Kinematics 7 1.1ContinuumIdealization............................. 7 1.2Descriptionofmotionofmaterialpointsinabody.............. 7 1.3Referentialandspatialdescriptions....................... 9 1.4 Displacement, Velocity, and Acceleration . 10 1.4.1 Displacement............................... 10 1.4.2 Velocity and Acceleration . 10 1.5DeformationGradientoftheMotion...................... 11 1.6 Transformations of infinitesimal areas and volumes during the deformation 13 1.6.1 Transformationofinfinitesimalareas . 13 1.6.2 Transformation of infinitesimal volumes . 13 1.7MeasuresofStretchandStrain......................... 14 1.7.1 Material stretch and strain tensors, C and E ............. 14 1.7.2 Spatial Stretch and Strain Tensors, b,e................. 16 1.8 Polar Decomposition of the Deformation Gradient Tensor . 17 1.9OtherStrainMeasures.............................. 19 1.10VelocityGradient................................. 19 1.10.1 Vorticity Vector and Vorticity Tensor . 22 1.11SpecialMotions.................................
    [Show full text]
  • Mechanics of Compliant Structures
    CHAPTER 5 Mechanics of compliant structures C. H. M. Jenkins1, W. W. Schur2 & G. Greshchik3 1Compliant Structures Laboratory, Mechanical Engineering Department South Dakota School of Mines and Technology, Rapid City, USA. 2Balloon Projects Branch Code 842, NASA Wallops Flight Facility, Wallops Island, USA. 3Center for Aerospace Structures, University of Colorado, Boulder, USA. Abstract Biological organisms must be structurally efficient. Hence it is not surprising that nature has embraced the compliant membrane structure as a central element in higher biological forms. This chapter discusses the unique challenges of modeling the mechanical behavior of compliant membrane structures. In particular, we focus on several special characteristics, such as large deformation, lack of bending rigidity, material nonlinearity, and computational schemes. 1 Introduction 1.1 Motivation Biological organisms must be structurally efficient. They have benefited from millions of years of evolution to achieve, for example, high load carrying capacity per unit weight. It is not surprising then to find that compliant membrane structures play a significant role in nature, for membranes are among the most efficient of structural elements. Long ago, engineers adopted membrane structures for solutions where structural efficiency was critical. As an example, consider the working balloon. For a given balloon volume, the total lift available is fixed, implying a zero-sum trade between structure and payload. Modern high-altitude scientific balloons (Fig. 1), made of thin polymer film fractions of a millimeter thick with areal densities of a few grams per square meter, support payloads of a few tons! WIT Transactions on State of the Art in Science and Eng ineering, Vol 20, © 2005 WIT Press www.witpress.com, ISSN 1755-8336 (on-line) doi:10.2495/978-1-85312-941-4/05 86 Compliant Structures in Nature and Engineering Figure 1: Preparing to launch a high-altitude scientific balloon in Antarctica (courtesy NASA).
    [Show full text]
  • Absolute Temperature, 210, 243, 362 Acceleration, 6, 152, 267, 361
    Cambridge University Press 978-1-107-02543-1 - An Introduction to Continuum Mechanics, Second Edition J. N. Reddy Index More information Index Absolute temperature, 210, 243, 362 orthonormal, 21{23, 26, 74, 78, 82, Acceleration, 6, 152, 267, 361 158, gravitational, 371 reciprocal, 19, 20, 43 vector, 182 unitary, 21, 42, 43 Adiabatic, 210 Beam, 2, 7, 278{284, 303{306 Adjoint, 38 bending, 278 Adjunct, 38 theory, 124, 278, 280, 305, 328, 351 Airy stress function, 301{309 Beltrami{Michell equations, 270, 271, 352 Algebraic Beltrami's equations, 271 equations, 328 Bending moment, 217, 280, 288, 305, multiplicity, 70 321{325, 346 Almansi{Hamel strain tensor, 104 Bernoulli's equations, 385 Alternating symbol, 25 Betti's reciprocity theorem, 288 Ampere's law, 254 Bianchi formulas, 147 Analytical solution, 272, 369, 379 Biaxial state of strain, 146 Angle of twist, 308 Biharmonic Angular equation, 301 momentum, 4, 169, 193{206 operator, 301 velocity, 14 Bilinear form, 326 Anisotropic, 222, 238, 258 Bingham model, 247 Antisymmetric tensor, see skew symmet- Body couple, 5, 169, 173, 205 ric Body force, 152, 197 Applied mechanics, 255 Boundary Approximate solution, 272 conditions, 364 Approximation functions, 328 curve, 18 Area change, 97 value problem, 265 Associative, 35 Bulk viscosity, 244 Asymmetry, 292 Axial vector, 119 Axisymmetric Cable supported beam, 282 body, 148 Caloric equation of state, 357 boundary condition, 292 Cantilever beam, 123, 217, 287, 290, 303, flow, 364 344 geometry, 296 Carreau model, 247 heat conduction, 376 Cartesian, 21,
    [Show full text]
  • Hyperelasticity Primer Robert M
    Hyperelasticity Primer Robert M. Hackett Hyperelasticity Primer Second Edition Robert M. Hackett Department of Civil Engineering The University of Mississippi University, MS, USA ISBN 978-3-319-73200-8 ISBN 978-3-319-73201-5 (eBook) https://doi.org/10.1007/978-3-319-73201-5 Library of Congress Control Number: 2018930881 © Springer International Publishing AG, part of Springer Nature 2016, 2018 This work is subject to copyright. All rights are reserved by the Publisher, whether the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of illustrations, recitation, broadcasting, reproduction on microfilms or in any other physical way, and transmission or information storage and retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology now known or hereafter developed. The use of general descriptive names, registered names, trademarks, service marks, etc. in this publication does not imply, even in the absence of a specific statement, that such names are exempt from the relevant protective laws and regulations and therefore free for general use. The publisher, the authors and the editors are safe to assume that the advice and information in this book are believed to be true and accurate at the date of publication. Neither the publisher nor the authors or the editors give a warranty, express or implied, with respect to the material contained herein or for any errors or omissions that may have been made. The publisher remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Printed on acid-free paper This Springer imprint is published by the registered company Springer International Publishing AG part of Springer Nature.
    [Show full text]
  • Mechanical Constitutive Models for Engineering Materials
    MECHANICAL CONSTITUTIVE MODELS FOR ENGINEERING MATERIALS by Behzad Rohani Soils and Pavements Laboratory U. S. Army Engineer Waterways Experiment Station P. O. Box 631, Vicksburg, Miss. 39180 September 1977 Final Report Approved For Public Release; Distribution Unlimited TA Prepared for Assistant Secretary of the Army (R&-D) 7 Department of the Army ,W34m Washington, D. C. 203I0 S-77-19 1977 Under Project 4 A I6 II0 IA 9 ID LIBRARY JAN 2 1970 Bureau of Reclamation Pnnvo*' f*- »^ ¡-~ J Destroy this report when no longer needed. Do not return it to the originator. V BUREAU OF RECLAMATION DENVER LIBRARY 92051290 V Unclassified ß ‘[ SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) READ INSTRUCTIONS REPORT DOCUMENTATION PAGE BEFORE COMPLETING FORM 1 .^REPORT NUMBER 2. GOVT ACCESSION NO. 3. RECIPIENT'S CATALOG NUMBER vb Miscellaneous^Paper S-77-19 A 4. T IT L E (and Subtitle) ' 5. TYPE OF REPORT & PERIOD COVERED 3 MECHANICAL CONSTITUTIVE MODELS FOR ENGINEERING £ Final Report) MATERIALS'-— .--------------------- ------ $ 6. PERFORMING ORG. REPORT NUMBER 7. A U T H O R S 8. CONTRACT OR GRANT NUMBER^) ( Behzad Rohani /" 4 9. PERFORMING ORGANIZATION NAME AND ADDRESS 10. PROGRAM ELEMENT, PROJECT, TASK AREA & WORK UNIT NUMBERS U. S.^Army Engineer Waterways Experiment Station Soils and Pavements Laboratory ^ ^ Project UA161101A91B P. 0. Box 631, Vicksburg, Miss. 39180 11. CONTROLLING OFFICE NAME AND ADDRESS / 12. REPORT DATE Assistant Secretary of the Army (R&D) ’ September 1977 t/ Department of the Army, Washington, D. C. 20310 13. NUMBER OF PAGES ' l8l 14. MONITORING AGENCY NAME & ADDRESSf/f different from Controlling Office) 15.
    [Show full text]
  • Notes on Rate Equations in Nonlinear Continuum Mechanics
    Notes on rate equations in nonlinear continuum mechanics∗ Daniel Aubram† September 29, 2017 Abstract: The paper gives an introduction to rate equations in nonlinear continuum mechanics which should obey specific transformation rules. Emphasis is placed on the geometrical nature of the operations involved in order to clarify the different concepts. The paper is particularly concerned with common classes of constitutive equations based on corotational stress rates and their proper implementation in time for solving initial boundary value problems. Hypoelastic simple shear is considered as an example application for the derived theory and algorithms. Contents 1 Introduction 2 2 Continuum Mechanics 2 2.1 MotionofaBody................................... ........ 2 2.2 DeformationGradientandStrain . .............. 4 2.3 StressandBalanceofMomentum. ............ 6 2.4 Constitutive Theory and Frame Invariance . ................. 8 3 Rates of Tensor Fields 10 3.1 Fundamentals.................................... ......... 10 3.2 CorotationalRates ............................... ........... 12 4 Rate Constitutive Equations 15 4.1 Hypoelasticity.................................. ........... 16 4.2 Hypoelasto-Plasticity. .............. 17 4.3 Hypoplasticity .................................. .......... 19 5 Objective Time Integration 20 5.1 Fundamentals and Geometrical Setup . .............. 20 5.2 AlgorithmofHughesandWinget . ............ 24 5.3 Algorithms Using a Corotated Configuration . ................ 25 6 Hypoelastic Simple Shear 28 6.1 AnalyticalSolution.
    [Show full text]
  • Biomechanical Aspects of Soft Tissues
    Biomechanical aspects of soft tissues Benjamin Loret Université de Grenoble, France Fernando M.F. Simões Instituto Superior Técnico Universidade de Lisboa, Portugal November 19, 2016 Crée comme Dieu Ordonne comme un roi Travaille comme un esclave Constantin Brancusi Contents Page Contents 4 1 Biomechanical topics in soft tissues 23 1.1 Introduction.................................... ........ 23 1.2 Diffusion, convection, osmosis: towards wearable artificialkidneys. 24 1.3 Issuesindrugdelivery . .. .. .. .. .. .. .. .. .. .. .. .......... 25 1.4 Macroscopic models of tissues, interstitium and membranes ................. 27 1.5 Energy couplings, passive and active transports . .................. 28 1.6 Moregeneraldirectand reversecouplings . ............... 29 1.7 Tissue engineering re-directed to tumor tissue exploration .................. 29 1.8 Frommechanicstobiomechanics . ........... 30 1.9 Mechanisms of injury of the knee, fracture mechanics . .................. 31 1.10 Water and solid constituents of soft tissues . ................. 33 I Solids and multi-species mixtures as open systems : a continuum mechanics perspective 35 2 Elements of continuum mechanics 37 2.1 Algebraic relations and algebraic operators . ................. 38 2.1.1 Notations and conventions for tensor calculus . ............... 38 2.1.2 Algebraicoperators . .. .. .. .. .. .. .. .. .. .. .. .. ....... 39 2.2 Characteristic polynomial, eigenvalues and eigenvectors ................... 42 2.2.1 Eigenvaluesandeigenvectors . .......... 42 2.2.2 Characteristic polynomial
    [Show full text]
  • Analysis of Deep Drawing Process by Using Hypo- and Hyperelastic Based Plasticity Models
    ANALYSIS OF DEEP DRAWING PROCESS BY USING HYPO- AND HYPERELASTIC BASED PLASTICITY MODELS A THESIS SUBMITTED TO THE GRADUATE SCHOOL OF NATURAL AND APPLIED SCIENCES OF MIDDLE EAST TECHNICAL UNIVERSITY BY İLKE BARIŞ CANBAKAN IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE IN MECHANICAL ENGINEERING AUGUST 2019 Approval of the thesis: ANALYSIS OF DEEP DRAWING PROCESS BY USING HYPO- AND HYPERELASTIC BASED PLASTICITY MODELS submitted by İLKE BARIŞ CANBAKAN in partial fulfillment of the requirements for the degree of Master of Science in Mechanical Engineering Department, Middle East Technical University by, Prof. Dr. Halil Kalıpçılar Dean, Graduate School of Natural and Applied Sciences Prof. Dr. Sahir Arıkan Head of Department, Mechanical Engineering Prof. Dr. Haluk Darendeliler Supervisor, Mechanical Engineering, METU Examining Committee Members: Prof. Dr. Süha Oral Mechanical Engineering, METU Prof. Dr. Haluk Darendeliler Mechanical Engineering, METU Prof. Dr. Suat Kadıoğlu Mechanical Engineering, METU Prof. Dr. Can Çoğun Mechanical Engineering, Çankaya University Assist. Prof. Dr. Orkun Özşahin Mechanical Engineering, METU Date: 06.08.2019 I hereby declare that all information in this document has been obtained and presented in accordance with academic rules and ethical conduct. I also declare that, as required by these rules and conduct, I have fully cited and referenced all material and results that are not original to this work. Name, Surname: İlke Barış Canbakan Signature: iv ABSTRACT ANALYSIS OF DEEP DRAWING PROCESS BY USING HYPO- AND HYPERELASTIC BASED PLASTICITY MODELS Canbakan, İlke Barış Master of Science, Mechanical Engineering Supervisor: Prof. Dr. Haluk Darendeliler August 2019, 103 pages Estimation of formability is important for sheet metal forming in order to get desired results with minimum amount of experiments.
    [Show full text]
  • Linear and Non Linear Analysis of Central Crack Propagation in Polyurethane Material - a Comparison
    Proceedings of the World Congress on Engineering 2012 Vol III WCE 2012, July 4 - 6, 2012, London, U.K. Linear and Non Linear Analysis of Central Crack Propagation in Polyurethane Material - A Comparison Mubashir Gulzar For soft solids, e.g. the case under discussion; modulus is Abstract— Crack formation in polyurethane material is an much lessto the physical value for cohesive strength. exceedingly non-linear and complex phenomenon. A 2-D plain Consequently, cracks undergo a very large deformation strain, linear and non-linear analysis to understand the before fracture and blunting at crack tip is also in question influence of central cracks and crack sizing is discussed under [2]. So if we considered the experimental technique for Mode 1 of crack displacement i.e. Opening Mode. Cracks in solid rocket propellant disturb the geometry and crack propagation analysis, it is crucial to understand photo require design characteristics .This research work deals with elasticity theory in this regard. The basic principles behind investigation of stress intensity factor determined Photoelasticity have been around for about 200 yrs [3]. experimentally and Finite Element based simulations, which Transmission Polariscope provides an experimental can be used to help explain the severity of the associated crack platform to apply Photo elasticity theory on Photo elastic in a given material. materials. The Photo elasticity model is an extremely This investigation involves single polyurethane material sheets versatile, simple, and powerful source tool for helping the which have Photo elastic characteristics. The global geometry engineer or designer by providing the design data before used for this purpose is rectangular.
    [Show full text]
  • Stability of Long Transversely-Isotropic Elastic Cylindrical Shells Under Bending
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Elsevier - Publisher Connector International Journal of Solids and Structures 47 (2010) 10–24 Contents lists available at ScienceDirect International Journal of Solids and Structures journal homepage: www.elsevier.com/locate/ijsolstr Stability of long transversely-isotropic elastic cylindrical shells under bending Sotiria Houliara, Spyros A. Karamanos * Department of Mechanical Engineering, University of Thessaly, 38334 Volos, Greece article info abstract Article history: The present paper investigates buckling of cylindrical shells of transversely-isotropic elastic material sub- Received 4 December 2008 jected to bending, considering the nonlinear prebuckling ovalized configuration. A large-strain hypoelas- Received in revised form 1 September 2009 tic model is developed to simulate the anisotropic material behavior. The model is incorporated in a Available online 11 September 2009 finite-element formulation that uses a special-purpose ‘‘tube element”. For comparison purposes, a hyperelastic model is also employed. Using an eigenvalue analysis, bifurcation on the prebuckling oval- Keywords: ization path to a uniform wrinkling state is detected. Subsequently, the postbuckling equilibrium path is Stability traced through a continuation arc-length algorithm. The effects of anisotropy on the bifurcation moment, Buckling the corresponding curvature and the critical wavelength are examined, for a wide range of radius-to- Postbuckling Imperfection sensitivity thickness ratio values. The calculated values of bifurcation moment and curvature are also compared Anisotropy with analytical predictions, based on a heuristic argument. Finally, numerical results for the imperfection Finite elements sensitivity of bent cylinders are obtained, which show good comparison with previously reported asymp- Cylindrical shell totic expressions.
    [Show full text]
  • Development and Implementation of Robust Large Deformation and Contact Mechanics Capabilities in Process Modelling of Composites
    DEVELOPMENT AND IMPLEMENTATION OF ROBUST LARGE DEFORMATION AND CONTACT MECHANICS CAPABILITIES IN PROCESS MODELLING OF COMPOSITES by AMIR OSOOLY B.Sc., Sharif University of Technology, 1989 M.Sc., Sharif University of Technology, 1993 A thesis submitted in partial fulfilment of the requirements for the degree of Doctor of Philosophy in The Faculty of Graduate Studies (Civil Engineering) The University of British Columbia (Vancouver) April 2008 © Amir Osooly, 2008 Abstract Autoclave processing of large scale, one-piece structural parts made of carbon fiber- reinforced polymer composite materials is the key to decreasing manufacturing costs while at the same time increasing quality. Nonetheless, even in manufacturing simple flat parts, residual strains and stresses are unavoidable. For structural design purposes and to aid in the assembly procedures, it is desirable to have proven numerical tools that can be used to predict these residual geometrical and material properties in advance, thus avoid the costly experimental trial and error methods. A 2-D finite element-based code, COMPRO, has previously been developed in-house for predicting autoclave process-induced deformations and residual stresses in composite parts undergoing an entire cure cycle. To simulate the tool-part interaction, an important source of residual deformations/stresses, COMPRO used a non-zero thickness elastic shear layer as its only interface option. Moreover, the code did not account for the large deformations and strains and the resulting nonlinear effects that can arise during the early stages of the cure cycle when the material is rather compliant. In the present work, a contact surface employing a penalty method formulation is introduced at the tool-part interface.
    [Show full text]