Optimal thickness of a cylindrical shell under dynamical loading
|
|
- Stephen Robbins
- 6 years ago
- Views:
Transcription
1 Optimal thickness of a cylindrical shell under dynamical loading Paul Ziemann Institute of Mathematics and Computer Science, E.-M.-A. University Greifswald, Germany paul.ziemann@uni-greifswald.de Abstract. In this paper an optimization problem for a cylindrical shell under an applied dynamical loading is discussed. The aim is to look for an optimal thickness of a shell to minimize the deformation when the shell is subject to an external timedependent force together with an initial impact. The volume of the shell has to stay constant during optimization. A model for the deflection is derived and analyzed and a corresponding optimal control problem is formulated. Numerical solutions are presented for a cylindrical panel with a body falling down onto the opposite side of the shell s only clamped edge and then staying there. Keywords: Optimal control of PDE, Shape optimization, Linear elasticity. Introduction In this paper we discuss an optimization problem in linear elasticity, particularly in shape optimization. The aim is to look for an optimal thickness of a cylindrical shell to minimize the deformation under an initial impact and a (time-dependent) external force, i.e. in a dynamical setup. As an additional restriction, the volume of the shell is prescribed. The deflection is modelled using a generalization of the steady-state basic shell model [] which makes use of the Hypothesis from Mindlin and Reissner but no further simplifications like in Naghdy or Koiter models. Optimization results for steady-state can be found in [] and for dynamic conditions with rotational symmetric loading in [3]. After stating the deflection model and the optimal control problem we give some analytical results regarding existence and uniqueness of solutions as well as continuity and differentiability of the deflection with respect to the thickness of the shell. We then focus on the numerical solution to an example of a cylindrical panel with one hardclamped edge and a body falling down onto the shell at the opposite side of that edge leading to an initial impact and then acting with its weight force. This problem was formulated by J. Lellep from University of Tartu.
2 . Geometrical description of the shell For the geometrical description, we first need a chart describing the midsurface of the shell. Let ω = [0, l] [ϕ a, ϕ b ] be the parameter region and z : ω R 3, z(x, ϕ) = (x, R cos ϕ, R sin ϕ) be a mapping to describe the midsurface of a cylindrical panel with length l and radius R. We call S = z(ω) the midsurface of the shell. The vectors a = (, 0, 0) and a = R (0, sin ϕ, cos ϕ) form a local contravariant basis on S. Additionally, we consider an orthonormal vector a 3 = (0, cos ϕ, sin ϕ). To describe the shell body we introduce τ : S R +, τ W, (S) as the thickness of the shell and define the 3D-reference domain Ω (τ) := { (x, ϕ, h) R 3 (x, ϕ) ω, h together with the mapping Φ (τ) : Ω (τ) R 3 ( τ(x, ϕ), τ(x, ϕ) Φ (τ) (x, ϕ, h) = z(x, ϕ) + ha 3 = (0, (R + h) cos ϕ, (R + h) sin ϕ). () We call B (τ) := Φ (τ) (Ω (τ) ) the shell body. Let us denote a local covariant basis with g = (, 0, 0), g = (R + h)(0, sin ϕ, cos ϕ) and g 3 = a 3. In the context of shells, we assume τ to be small w.r.t. other dimensions, in particular τ(x, ϕ) < R. 3. Modeling the displacement We consider a small displacement U : B (τ) R 3 of the shell body. For modeling we use the Reissner-Mindlin kinematical assumptions, which state that normals to the midsurface remain straight and unstretched during deformation. This leads to the displacement ansatz )} () U(x, ϕ, h) = u(x, ϕ) + hθ(x, ϕ) (3) with u = u a +u a +u 3 a 3 describing a translation of all points on a line normal to the midsurface in z(x, ϕ) and θ = θ a + θ a representing a rotation vector, all in local coordinates. We introduce the space of admissible displacements V := { (u, θ) u = (u, u 3 ) H (S) H (S), θ H (S) } BC (4) where H (S) and H (S) are Sobolev-spaces for scalar functions and first order tensors on the midsurface, resp. Let us assume for the boundary conditions BC that the shell body is hardclamped over the edge x = 0, i.e. (u, θ) x=0 = 0. We next consider the linear 3D-Green-Lagrange-strain tensor which is given by e ij = (g i U,j + g j U,i ), i, j =,, 3, (5)
3 where U,i means the partial derivative of U w.r.t. to the i th coordinate. Calculation leads to e u, + hθ, e := e e = u, + Ru 3 + h ( θ, + R u ), + u 3 + h R θ, (u, + u, ) + h ( θ, + θ, + R u ) (6), + h R θ, and shear strain ζ = ( e3 e 3 ) = ( ) θ + u 3, θ + u 3, R u. (7) By Hooke s Law we get a relationship between strains and stresses σ ij. Adding the assumption that the normal stress σ 33 is zero (plain-stress) and using Voigtnotation we later can represent the equilibrium conditions for isotropic material in a short way with help of the matrices L + L L (R+h) 0 C := L (R+h) (L (R+h) 4 + L ) L (R+h) D := ( ) 4L 0 4, 0 L (R+h) (8) Eν where L = (+ν)( ν) and L = E (+ν) are the D-Lamé constants expressed with Young s modulus E and Poisson s ratio ν. To account for dynamical loadings and impacts we introduce a time-dependence on (u(t), θ(t)) =: y(t), t [0, T ] and assume the deformation to be in the space W (0, T ) = {y y L (0, T ; V), t y L (0, T ; H), t y(0, T ; V } (9) where L (0, T ; ) is a space of abstract functions which are square Bochnerintegrable, H = L (S) L (S) L (S) is a rigged Hilbert space and V is the dual space of V. By starting from full 3D-equation of momentum t U(t) div(σ(u(t))) = F (t) in B (τ) for t [0, T ] we can first switch to a weak formulation and then simplify using the above mentioned Reissner-Mindlin kinematics and plain-stress assumption, but no further shell-related simplifications, e.g. Naghdy or Koiter models. We can state dynamic deformation equations in variational form: Given a surface loading f H (0, T ; L (S)), find y W (0, T ), so that ρm (τ) t y(t) + L (τ) (y(t)) = F (t) in V for all t [0, T ], (0) together with initial conditions y(0) = 0, t y(0) = g /τ where g V represents an initial impact on the shell body and ρ is the density. In state equation (0) 3
4 we have the linear operator M (τ) : V V arising from dynamical part, and the linear operator L (τ) : V V, τ τ 0 3 R u M (τ) (u, θ) = 0 τ 0 u 3 () θ τ 3 R 0 τ 3 L (τ) (u, θ)(v, ψ) = e(u, θ) T Ce(v, ψ) + ζ(u, θ) T Dζ(v, ψ) dv, Ω (τ) () where this part is also known as basic-shell-model, introduced in [] and analyzed in []. The functional F (t) V represents a loading applied on the midsurface in orthogonal direction, i.e. F (t)(v, ψ) = fv 3 ds. (3) Finally, we can state our optimization problem ω T min J(u, θ; τ) := u dt + λ τ W, (S) 0 τ H (S) y=(u,θ) W (0,T ) s.t. ρm(τ) t y(t) + L (τ) (y(t)) = F (t) in V for t [0, T ] y(0) = 0, t y(0) = g τ τ min τ(x, ϕ) τ max in S, ω τ ds = C (4) 4. Analysis of the problem We want to study state equation (0) according to existence and uniqueness of a solution and Gâteaux-differentiability of the solution w.r.t. thickness. Let us collect all the thickness restrictions in U ad := {τ W, (S) τ min τ(x, ϕ) τ max in S, ω τ ds = C}. Theorem. Let f H (0, T ; L (S)) and g V be given. For τ U ad, the variational problem (0) has an unique solution in W (0, T ). Proof. Proven by using techniques from [4] together with analysis of L (τ) from [] and analysis of M (τ). 4
5 Definition. We call the operator G : U ad W (0, T ), which maps the thickness τ to the solution y (τ) W (0, T ) of (0) control-to-state-operator. Theorem. The operator G is continuous and Gâteaux-differentiable. directional derivative G (τ)q = ȳ = (ū, θ) is given as the solution to The ρm (τ) t ȳ(t) + L (τ) (ȳ(t)) = ρ(m (τ) q) t y (τ) (t) Z q,y(τ) (t) in V for t [0, T ] ȳ(0) = 0, t ȳ(0) = g τ (5) where M (τ) q means Gâteaux-derivative of M (τ) w.r.t. τ in direction q and Z q,(u(τ),θ (τ) )(v, ψ) = [( e(u(τ), θ (τ) ) T Ce(v, ψ) + ζ(u (τ), θ (τ) ) T Dζ(v, ψ) ) ( + h ) R ω h i {± τ } h=h i ] q ds. Proof. By using results from [] for L (τ) and differentiability of M (τ) w.r.t. τ together with estimates from theory of hyperbolic PDE regarding dependence of the solution from initial data. Theorem 3. The reduced objective J s : U ad R, J s (τ) := J(G(τ), τ) is Gâteaux-differentiable and it holds J s(τ)q = T 0 ( ) ρ(m (τ) q) t y (τ) (t), p(t) + Z q,y H (τ) (t)p(t) dt + λ (τ, q) H (S) with adjoint p W (0, T ) as solution to (6) ρm (τ) t p(t) + L (τ) (p(t)) = (u (τ) (t), 0) in V f.a. t [0, T ] p(t ) = 0, t p(t ) = 0 (7) Proof. By direct calculation. 5. Numerical solution We solve the reduced problem min τ Uad J s (τ) numerically in Fortran. The spatial discretization is done with 9-node-biquadratic shell elements. For time domain we use a Crank-Nicolson scheme. The resulting system is solved using Pardiso and the optimization is done using IpOpt. Note that the key for a successful optimization process is an expression for the gradient given to IpOpt that is based on (6) and not on finite differences. The optimization starts on a coarse grid consisting of 89 nodes and ends up on a fine grid with about 4000 nodes. 5
6 z y x Figure : Initial configuration with loaded region (yellow) and clamped edge (red) ϕ x Figure : Optimal thickness over parameter domain (color), loaded region (black line) and clamped edge (red line) We consider a quarter tube of length, i.e. ω = [0, ] [ π /4, 3π /4] with radius R =, E = 0, ν = 0.3, ρ = 7.8 and volume C = π /0 which corresponds to a constant thickness of 0.. Minimum and maximum thickness are 0.05 and 0.5, resp. and a small regularization parameter λ 0 5 is used. We study the time domain [0, ]. A mass falling onto the opposite side of the clamped edge x = 0 is giving an impact g(x, ϕ) = 0.a 3 for (x, ϕ) [0.9, π / 0.05] [, π / ]. The mass stays on the panel and acts with its weight force f(x, ϕ) = a 3 in the same region, see figure. The corresponding optimal thickness (solution to (4)) is shown in figure. We see a region of maximal thickness surrounding the loaded region as well as kind of trusses emerging to the loose edges of the shell. Like the load, the optimal thickness is symmetric w.r.t. to the ϕ-coordinate. Moreover, a small moderate thick region is placed nearer to the clamped edge. This thickness distribution could also be an interesting starting point for topology optimization where regions with minimal thickness are considered as holes. References [] Chapelle, D. and Bathe, K.J. The Finite Element Analysis of Shells - Fundamentals, Springer, 003. [] Ziemann, P.: Optimal Thickness of a Cylindrical Shell, Mathematics and its Applications, Vol. 6, No., 05, pp [3] Nestler, P.: Optimal Thickness of a Cylindrical Shell Under a Time-Dependent Force, Journal of Optimization Theory and Applications, 03, pp. 3. [4] Wloka, J. Partial Differential Equations, Cambridge University Press,
3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More informationLinearized theory of elasticity
Linearized theory of elasticity Arie Verhoeven averhoev@win.tue.nl CASA Seminar, May 24, 2006 Seminar: Continuum mechanics 1 Stress and stress principles Bart Nowak March 8 2 Strain and deformation Mark
More informationAircraft Structures Kirchhoff-Love Plates
University of Liège erospace & Mechanical Engineering ircraft Structures Kirchhoff-Love Plates Ludovic Noels Computational & Multiscale Mechanics of Materials CM3 http://www.ltas-cm3.ulg.ac.be/ Chemin
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationChemnitz Scientific Computing Preprints
Arnd Meyer The Koiter shell equation in a coordinate free description CSC/1-0 Chemnitz Scientific Computing Preprints ISSN 1864-0087 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific
More information3. The linear 3-D elasticity mathematical model
3. The linear 3-D elasticity mathematical model In Chapter we examined some fundamental conditions that should be satisfied in the modeling of all deformable solids and structures. The study of truss structures
More informationFundamentals of Linear Elasticity
Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy
More informationThe Finite Element Method for Computational Structural Mechanics
The Finite Element Method for Computational Structural Mechanics Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) January 29, 2010 Martin Kronbichler (TDB) FEM for CSM January
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationChemnitz Scientific Computing Preprints
Arnd Meyer The linear Naghdi shell equation in a coordinate free description CSC/13-03 Chemnitz Scientific Computing Preprints Impressum: Chemnitz Scientific Computing Preprints ISSN 1864-0087 (1995 2005:
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationTable of Contents. II. PDE classification II.1. Motivation and Examples. II.2. Classification. II.3. Well-posedness according to Hadamard
Table of Contents II. PDE classification II.. Motivation and Examples II.2. Classification II.3. Well-posedness according to Hadamard Chapter II (ContentChapterII) Crashtest: Reality Simulation http:www.ara.comprojectssvocrownvic.htm
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
More informationFinite Elements for Elastic Shell Models in
Elastic s in Advisor: Matthias Heinkenschloss Computational and Applied Mathematics Rice University 13 April 2007 Outline Elasticity in Differential Geometry of Shell Geometry and Equations The Plate Model
More informationSymmetry and Properties of Crystals (MSE638) Stress and Strain Tensor
Symmetry and Properties of Crystals (MSE638) Stress and Strain Tensor Somnath Bhowmick Materials Science and Engineering, IIT Kanpur April 6, 2018 Tensile test and Hooke s Law Upto certain strain (0.75),
More informationVariational principles in mechanics
CHAPTER Variational principles in mechanics.1 Linear Elasticity n D Figure.1: A domain and its boundary = D [. Consider a domain Ω R 3 with its boundary = D [ of normal n (see Figure.1). The problem of
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationSoft 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 informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto Gutierrez-Miravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationPEAT SEISMOLOGY Lecture 2: Continuum mechanics
PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a
More informationA HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS
A HIGHER-ORDER BEAM THEORY FOR COMPOSITE BOX BEAMS A. Kroker, W. Becker TU Darmstadt, Department of Mechanical Engineering, Chair of Structural Mechanics Hochschulstr. 1, D-64289 Darmstadt, Germany kroker@mechanik.tu-darmstadt.de,
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationMechanics of Inflatable Fabric Beams
Copyright c 2008 ICCES ICCES, vol.5, no.2, pp.93-98 Mechanics of Inflatable Fabric Beams C. Wielgosz 1,J.C.Thomas 1,A.LeVan 1 Summary In this paper we present a summary of the behaviour of inflatable fabric
More informationMathematical modelling and numerical analysis of thin elastic shells
Mathematical modelling and numerical analysis of thin elastic shells Zhang Sheng Wayne State University, Michigan Zhang Sheng (szhang@wayne.edu) Thin elastic shell May, 202 / 9 Models for mechanics of
More informationBasic Equations of Elasticity
A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ
More informationGame 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 informationElements of Continuum Elasticity. David M. Parks Mechanics and Materials II February 25, 2004
Elements of Continuum Elasticity David M. Parks Mechanics and Materials II 2.002 February 25, 2004 Solid Mechanics in 3 Dimensions: stress/equilibrium, strain/displacement, and intro to linear elastic
More informationNDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.
CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo
More informationBy drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.
Mohr s Circle By drawing Mohr s circle, the stress transformation in -D can be done graphically. σ = σ x + σ y τ = σ x σ y + σ x σ y cos θ + τ xy sin θ, 1 sin θ + τ xy cos θ. Note that the angle of rotation,
More informationKirchhoff Plates: Field Equations
20 Kirchhoff Plates: Field Equations AFEM Ch 20 Slide 1 Plate Structures A plate is a three dimensional bod characterized b Thinness: one of the plate dimensions, the thickness, is much smaller than the
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationQuestions (And some Answers) (And lots of Opinions) On Shell Theory
1 Questions (And some Answers) (And lots of Opinions) On Shell Theory 2 Shells Q: What is a shell? A: A three-dimensional elastic body occupying a thin neighborhood of a two-dimensional submanifold of
More informationMacroscopic theory Rock as 'elastic continuum'
Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave
More informationABHELSINKI UNIVERSITY OF TECHNOLOGY
ABHELSINKI UNIVERSITY OF TECHNOLOGY TECHNISCHE UNIVERSITÄT HELSINKI UNIVERSITE DE TECHNOLOGIE D HELSINKI A posteriori error analysis for the Morley plate element Jarkko Niiranen Department of Structural
More informationContinuum mechanics V. Constitutive equations. 1. Constitutive equation: definition and basic axioms
Continuum mechanics office Math 0.107 ales.janka@unifr.ch http://perso.unifr.ch/ales.janka/mechanics Mars 16, 2011, Université de Fribourg 1. Constitutive equation: definition and basic axioms Constitutive
More informationLecture 1 NONLINEAR ELASTICITY
Lecture 1 NONLINEAR ELASTICITY Soft-Matter Engineering: Mechanics, Design, & Modeling Mechanics Rigid/Hard Systems bending or torsion (but not stretching) small strains (~0.1% metals; ~1% plastics) linearized
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationGEOMETRIC NONLINEAR ANALYSIS
GEOMETRIC NONLINEAR ANALYSIS The approach for solving problems with geometric nonlinearity is presented. The ESAComp solution relies on Elmer open-source computational tool [1] for multiphysics problems.
More information6.1 Formulation of the basic equations of torsion of prismatic bars (St. Venant) Figure 6.1: Torsion of a prismatic bar
Module 6 Torsion Learning Objectives 6.1 Formulation of the basic equations of torsion of prismatic bars (St. Venant) Readings: Sadd 9.3, Timoshenko Chapter 11 e e 1 e 3 Figure 6.1: Torsion of a prismatic
More informationA Continuum Theory for the Natural Vibrations of Spherical Viral Capsids
Viral Capsids A Continuum Theory for the Natural Vibrations of Spherical Viral Capsids Francesco Bonaldi Master s Thesis in Mathematical Engineering Advisors Prof. Paolo Podio-Guidugli Università di Roma
More informationChapter 0. Preliminaries. 0.1 Things you should already know
Chapter 0 Preliminaries These notes cover the course MATH45061 (Continuum Mechanics) and are intended to supplement the lectures. The course does not follow any particular text, so you do not need to buy
More informationOPTIMAL THICKNESS OF A CYLINDRICAL SHELL AN OPTIMAL CONTROL PROBLEM IN LINEAR ELASTICITY THEORY
Ann. Acad. Rom. Sci. Ser. Math. Appl. ISSN 2066-6594 Vol. 4, No. 2 / 2012 OPTIMAL THICKNESS OF A CYLINDRICAL SHELL AN OPTIMAL CONTROL PROBLEM IN LINEAR ELASTICITY THEORY Peter Nestler Abstract In this
More informationElastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount
Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Min Ding* 1, and Jian Lin 2 1 MIT/WHOI Joint Program, 2 Woods Hole Oceanographic Institution *Woods Hole Oceanographic
More informationBasic concepts to start Mechanics of Materials
Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen
More informationSimulating Thin Shells with MPM
Simulating Thin Shells with MPM Biswajit Banerjee Center for the Simulation of Accidental Fires and xplosions University of Utah March 14, 2005 Outline The Problem. Shell Theory. MPM Formulation. Results.
More informationUNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES
UNCONVENTIONAL FINITE ELEMENT MODELS FOR NONLINEAR ANALYSIS OF BEAMS AND PLATES A Thesis by WOORAM KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More information16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations
6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationOn the characterization of drilling rotation in the 6 parameter resultant shell theory
On the characterization of drilling rotation in the 6 parameter resultant shell theory Mircea Birsan and Patrizio Neff Chair for Nonlinear Analysis and Modelling Faculty of Mathematics, University Duisburg-Essen,
More information202 Index. failure, 26 field equation, 122 force, 1
Index acceleration, 12, 161 admissible function, 155 admissible stress, 32 Airy's stress function, 122, 124 d'alembert's principle, 165, 167, 177 amplitude, 171 analogy, 76 anisotropic material, 20 aperiodic
More informationSTANDARD SAMPLE. Reduced section " Diameter. Diameter. 2" Gauge length. Radius
MATERIAL PROPERTIES TENSILE MEASUREMENT F l l 0 A 0 F STANDARD SAMPLE Reduced section 2 " 1 4 0.505" Diameter 3 4 " Diameter 2" Gauge length 3 8 " Radius TYPICAL APPARATUS Load cell Extensometer Specimen
More informationPEAT SEISMOLOGY Lecture 3: The elastic wave equation
PEAT8002 - SEISMOLOGY Lecture 3: The elastic wave equation Nick Rawlinson Research School of Earth Sciences Australian National University Equation of motion The equation of motion can be derived by considering
More informationElastic plastic bending of stepped annular plates
Applications of Mathematics and Computer Engineering Elastic plastic bending of stepped annular plates JAAN LELLEP University of Tartu Institute of Mathematics J. Liivi Street, 549 Tartu ESTONIA jaan.lellep@ut.ee
More informationKINEMATICS OF CONTINUA
KINEMATICS OF CONTINUA Introduction Deformation of a continuum Configurations of a continuum Deformation mapping Descriptions of motion Material time derivative Velocity and acceleration Transformation
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationMathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.
Mathematical Tripos Part IA Lent Term 205 ector Calculus Prof B C Allanach Example Sheet Sketch the curve in the plane given parametrically by r(u) = ( x(u), y(u) ) = ( a cos 3 u, a sin 3 u ) with 0 u
More informationA new approach for Kirchhoff-Love plates and shells
A new approach for Kirchhoff-Love plates and shells Walter Zulehner Institute of Computational Mathematics JKU Linz, Austria AANMPDE 10 October 2-6, 2017, Paleochora, Crete, Greece Walter Zulehner (JKU
More informationThe Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density
Applied Mathematics & Information Sciences 23 2008, 237 257 An International Journal c 2008 Dixie W Publishing Corporation, U. S. A. The Rotating Inhomogeneous Elastic Cylinders of Variable-Thickness and
More informationContinuum Mechanics and the Finite Element Method
Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after
More informationNONCLASSICAL MODELS IN THE SHELL THEORY WITH APPLICATIONS TO MULTILAYERED NANOTUBES
COMPDYN 0 3 rd ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis M. Fragiadakis V. Plevris eds. Corfu Greece 5-8 May 0 NONCLASSICAL
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationLecture 15 Strain and stress in beams
Spring, 2019 ME 323 Mechanics of Materials Lecture 15 Strain and stress in beams Reading assignment: 6.1 6.2 News: Instructor: Prof. Marcial Gonzalez Last modified: 1/6/19 9:42:38 PM Beam theory (@ ME
More informationINDEX 363. Cartesian coordinates 19,20,42, 67, 83 Cartesian tensors 84, 87, 226
INDEX 363 A Absolute differentiation 120 Absolute scalar field 43 Absolute tensor 45,46,47,48 Acceleration 121, 190, 192 Action integral 198 Addition of systems 6, 51 Addition of tensors 6, 51 Adherence
More informationa x Questions on Classical Solutions 1. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress
Questions on Classical Solutions. Consider an infinite linear elastic plate with a hole as shown. Uniform shear stress σ xy = T is applied at infinity. Determine the value of the stress σ θθ on the edge
More information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationMECHANICS OF MATERIALS
Third E CHAPTER 2 Stress MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf Lecture Notes: J. Walt Oler Texas Tech University and Strain Axial Loading Contents Stress & Strain:
More informationNonlinear bending analysis of laminated composite stiffened plates
Nonlinear bending analysis of laminated composite stiffened plates * S.N.Patel 1) 1) Dept. of Civi Engineering, BITS Pilani, Pilani Campus, Pilani-333031, (Raj), India 1) shuvendu@pilani.bits-pilani.ac.in
More informationRisk Averse Shape Optimization
Risk Averse Shape Optimization Sergio Conti 2 Martin Pach 1 Martin Rumpf 2 Rüdiger Schultz 1 1 Department of Mathematics University Duisburg-Essen 2 Rheinische Friedrich-Wilhelms-Universität Bonn Workshop
More informationME FINITE ELEMENT ANALYSIS FORMULAS
ME 2353 - FINITE ELEMENT ANALYSIS FORMULAS UNIT I FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PROBLEMS 01. Global Equation for Force Vector, {F} = [K] {u} {F} = Global Force Vector [K] = Global Stiffness
More informationBase force element method of complementary energy principle for large rotation problems
Acta Mech Sin (2009) 25:507 515 DOI 10.1007/s10409-009-0234-x RESEARH PAPER Base force element method of complementary energy principle for large rotation problems Yijiang Peng Yinghua Liu Received: 26
More informationPlane and axisymmetric models in Mentat & MARC. Tutorial with some Background
Plane and axisymmetric models in Mentat & MARC Tutorial with some Background Eindhoven University of Technology Department of Mechanical Engineering Piet J.G. Schreurs Lambèrt C.A. van Breemen March 6,
More informationA new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem
A new closed-form model for isotropic elastic sphere including new solutions for the free vibrations problem E Hanukah Faculty of Mechanical Engineering, Technion Israel Institute of Technology, Haifa
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationAnalysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method
Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method V.M.F. Nascimento Departameto de ngenharia Mecânica TM, UFF, Rio de Janeiro
More informationChapter 6 2D Elements Plate Elements
Institute of Structural Engineering Page 1 Chapter 6 2D Elements Plate Elements Method of Finite Elements I Institute of Structural Engineering Page 2 Continuum Elements Plane Stress Plane Strain Toda
More informationAdvanced Structural Analysis EGF Section Properties and Bending
Advanced Structural Analysis EGF316 3. Section Properties and Bending 3.1 Loads in beams When we analyse beams, we need to consider various types of loads acting on them, for example, axial forces, shear
More informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 1 June 2006 1.30 to 4.30 PAPER 76 NONLINEAR CONTINUUM MECHANICS Attempt FOUR questions. There are SIX questions in total. The questions carry equal weight. STATIONERY
More informationMAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design
MAE4700/5700 Finite Element Analsis for Mechanical and Aerospace Design Cornell Universit, Fall 2009 Nicholas Zabaras Materials Process Design and Control Laborator Sible School of Mechanical and Aerospace
More informationTopic 5: Finite Element Method
Topic 5: Finite Element Method 1 Finite Element Method (1) Main problem of classical variational methods (Ritz method etc.) difficult (op impossible) definition of approximation function ϕ for non-trivial
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has cross-sectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationSome Remarks on the Reissner Mindlin Plate Model
Some Remarks on the Reissner Mindlin Plate Model Alexandre L. Madureira LNCC Brazil Talk at the Workshop of Numerical Analysis of PDEs LNCC February 10, 2003 1 Outline The 3D Problem and its Modeling Full
More informationSec. 1.1: Basics of Vectors
Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x
More informationAnalytical formulation of Modified Upper Bound theorem
CHAPTER 3 Analytical formulation of Modified Upper Bound theorem 3.1 Introduction In the mathematical theory of elasticity, the principles of minimum potential energy and minimum complimentary energy are
More informationAdaptive Analysis of Bifurcation Points of Shell Structures
First published in: Adaptive Analysis of Bifurcation Points of Shell Structures E. Ewert and K. Schweizerhof Institut für Mechanik, Universität Karlsruhe (TH), Kaiserstraße 12, D-76131 Karlsruhe, Germany
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Computational Methods in Elasticity - Michel SALAUN
COMPUAIONAL MEHODS IN ELASICIY Michel SALAUN Départment de Mathématiques et Informatique, ENSICA, oulouse, France Keywords Elastodynamics, Energy minimization, Finite element method, Kirchhoff- Love model,
More informationPractice Final Examination. Please initial the statement below to show that you have read it
EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use
More informationEE C245 ME C218 Introduction to MEMS Design Fall 2007
EE C245 ME C218 Introduction to MEMS Design Fall 2007 Prof. Clark T.-C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 Lecture 13: Material
More informationMECh300H Introduction to Finite Element Methods. Finite Element Analysis (F.E.A.) of 1-D Problems
MECh300H Introduction to Finite Element Methods Finite Element Analysis (F.E.A.) of -D Problems Historical Background Hrenikoff, 94 frame work method Courant, 943 piecewise polynomial interpolation Turner,
More informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 57 AA242B: MECHANICAL VIBRATIONS Dynamics of Continuous Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations: Theory
More information