QUASISTATIC INDENTATION OF A RUBBER COVERED ROLL BY A RIGID ROLL - THE BOUNDARY ELEMENT SOLUTION
|
|
- Neil O’Neal’
- 5 years ago
- Views:
Transcription
1 259 QUASISTATIC INDENTATION OF A RUBBER COVERED ROLL BY A RIGID ROLL - THE BOUNDARY ELEMENT SOLUTION Batra Department of Engineering Mechanics, Univ. of Missouri-Rolla ABSTRACT The linear elastic problem involving the indentation of a compressible rubberlike layer bonded to a rigid cylinder and indented by another rigid cylinder is analyzed by the boundary element method. For the same contact width. the pressure at the contact surface is found to depend noticeably upon the thickness of the layer and the Poisson ratio of the material of the layer. Results computed and presented graphically include the pressure distribution over the contact surface. the shape of the indented surface and the stress distribution at the bond surface. INTRODUCTION Tr.:.ction in vehiclef, nip acti'jn in cylindrical ro:ls in the paper-making process and in the textile industry, and friction drives are examples of the kind of problem studied herein. Each of these problems involves indentation, by a steel or granite cylinder, of a rubberlike layer bonded to a cylindrical core made also of steel or granite. Hertz's (1881) solution for determining stresses in two bodies in contact is applicable when both bodies are homogeneous and linear elastic and the dimensions of the contact region are small compared with the dimensions of the bodies. The contact problem involving two rollers, one of which is assumed perfectly hard, in contact with their axes parallel is a special case of Hertz's problem and has been investigated by Thomas and Hoersch (193) on the assumption that plane strain state of deformation prevails. Solutions of the problem when the deformable ro1.ler is not homogeneous but consists of a thin elastic layer bonded to a hard supporting core has been given by Hannah (1951), and Hahn and Levinson (1974). Whereas Hannah's approach involves the numerical solution of an integral equation, Hahn and Levinson
2 26 use an Airy stress function and obtain the solution in the form of a double series one of which converges rather slowly In the aforementioned studies the rubberlike layer is taken to be linear elastic and the deformations involved are presumed to be infinitesimal so that the linear theory applies Aiso for the case of infinitesimal deformations, Batra et al. (1976) solved the problem by taking the rubberlike layer to be made of a thermorheologically simple material. In order to explore the effects of material and geometric nonlinearities, Batra (198, 1981) recently studied the problem in which the rubberlike layer is modeled as a homogeneous and incompressible or nearly incompressible Mooney-Rivlin material. We note that Batra used the finite element method to solve the problem Experimental work involving varying thicknesses of the rubberlike layer and different combinations of the diameters of mating cylinders has been carried out by Spengos (1965). However, the only material property listed for the rubber is the durometer hardness. For the linear elastic problem involving incompressible material. the durometer hardness is sufficient to find Young' modulus and hence solve the problem. Since the length to diameter ratio for the rolls used by Spengos is of the order of one, the assumption of plane strain state of deformation made in the work referenced above is not valid. For such geometrical configurations one needs to solve the three dimensional problem. It seems that the Boundary Element Method has an advantage over the Finite Element Method for three dimensional problems since at least the data preparation for the former method is easier. With the ultimate aim of solving the three dimensional problem, in this paper, we solve the two dimensional (plane strain) problem by the Boundary Element Method. Results obtained from the use of the boundary element method are found to compare very favorably with those obtained by Hahn and Levinson (1974). iormulationof the problem The system to be studied and the location of rectangular Cartesian axes are shown in Fig. 1. Since the length of rolls perpendicular to the plane of paper is large as compared to their diameters we assume that plane strain state of deformation prevails. Also the deformations of the rubberlike layer are presumed to be infinitesimal so that a linear theory applies. The rubberlike layer is assumed to be made of a homogeneous elastic material. In the absence of friction at the contact surface, the deformations of the rubberlike layer are symmetrical about the line joining the centers of two rolls. Neglecting the effect of body forces such as gravity, equations governing the quasistatic deformations of the rubber-
3 261, x Elastic Layer Rigid Material Fig. 1 System to be studied like layer are P = Po (1 - u..), (1) J.,J. (LijkR. uk,r.)'j =, (i,j = 1,2) (2) E\I LijkR. = (1+\1)(1-2\1) ij kr. E (3) + 2(1+\1) (csik csjr. + csir. csjk). The pertinent boundary conditions are: at the inner surface SI (xixi = Rl)' ui =, (4) at the outer surface S2 (xixi= RO); eioijnj =, (5) ni o..n. =, Xl > b (6) J.J J - nioijnj = P(xl' x2)' x2 < b, (7) lim P(xl' x2) =, x2 b, (8) and at the plane surface S3 (x2 = ); where u2 =, (9) 12 =, (1) ij = LijkR.uk,R.. (11)
4 262 In these equations, P is the mass density of the material particle in the deformed state, Po is its mass density in the undeformed state, y is the displacement of a material particle, E and v are respectively Young's modulus and Poisson's ratio for the material of the rubberlike layer, a comma followed by an index i indicates partial differentiation with respect to xi' the summation convention is used, is a unit vector tangent to the surface, U is an outward unit normal to the surface, 6ij is the Kronecker delta, cr.. is the Cauchy stress tensor and the function p gives the prsure distribution at the contact surface. The condition (8) implies that the normal pressure at the boundary of the contact surface vanishes. This ensures that a contact problem rather than a punch problem is being solved. e note that the semicontact width b and the pressure p (xl,x2) at the contact surface are not known a priori but are to be determind as a part of the solution. As has been employed by Hahn and Levinson (1974), the boundary condition (7) can be replaced by a displacement type boundary condition. We find it convenient to se the boundary condition (7) rather than the equivalent displacement type boundary condition. BOUNDARY ELEMENT FORMULATION OF THE PROBLEM Taking the inner product of both sides of Equation (2) with wi' integrating the resulting equation over the domain G consisting of the region occupied by the layer above the plane x2 = and using the divergence theorem, we obtain J (LijkR. Wi,j)'R. uk dg aj (gkuk - fiwi) ds (12) where i: LijkR- uk,r- Dj' gk = LijkR, Wi,j nr," (14) (13) Details of deriving Equation (12) from Equation (2) are grrei' by lrebbia 1978). We now choose.. to be a solution of where m) is the Dirac delta function and represents a unit load at point m in the k direction. With this choice of wi' Equation (12) becomes ui (m) + f f. wi ds = f gkuk ds (16) an J. an Note that wi and gi are displacements and tractions due to a unit concentrated load at the point m in the k direction. Considering unit forces acting in the three directions, Equation (14) can be written as.-
5 263 where Wi and Gi represent the displacements and tractions in the direction due to a unit force acting in the i direction. Equation (17) is valid for the particular point m where these forces are applied. Expressions for Wi and Gi are given in Brebbia's book (1978). In Equation (17) the integrations are over the boundary of the domain. At points of the boundary, either surface tractions fi or displacements ui or a suitable combination of the two are known. Equation (17) enables us to solve for the unknown surface tractions or unknown surface displacements at each point of the boundary. COMPUTATION AND DISCUSSION OF RESULTS In order to solve the problem by the boundary element method, the boundary of the rubberlike layer in the first quadrant is divided into a large number of segments as shown in Fig. 2. Displacements and surface tractions within each segment (boundary element) are assumed to be constant. The mesh is finer within approximately twice the semicontact width. As has been shown in earlier studies (Batra, 198, 1981), the boundary conditions on the vertical plane have essentially no effect on the pressure profile at the contat surface and the deformations of the rubberlike layer within the vicinity of the contact region. In the results presented herein the vertical plane is taken to be traction free. The computer program given in Brebbia's (1978) book for constant boundary elements has been modified to solve the present contact problem. Half nip width b and the form of the pressure function p are presumed. With this additional data the problem is well defined and can be solved. Having solved the problem a check is made to insure that the deformed surface in the assumed contact zone matches, within a prescribed tolerance, with the circular profile of the indentor and that the nodal point just outside the contact area has not penetrated into the indentor. If the second conditi9n is not satisfied implying thereby that the nodal point just outside the presumed contact width has p,1netrated into the indent;')]:, either the value of b is increased or the total load is decreased. However, if the second condition is satisfied and the first is not, the values of p at various nodal points are adjusted so that the deformed shape of the assumed contact area conforms to the circular profile of the indentor. The deformed surface of the roll cover is taken to match with the profile of the indentor if the distance of each nodal point on the contact surface from the indentor is less than 1.5 percent. of the indentation Uo (see Fig. 1). Usually, with a little experience one can make pretty good est1mates of band p (x1,x2) so that the entire process converges in four or five iterations. The total load P is obtained by integrating p over the contact area.
6 264 Test cases To ensure that the discretization of the boundary as adequate, the contact problem with RO = R = rom, R1 = 441 rom, b = 22.9 rom, E = 8.96 x 15 N/m2, and v =.3 was solved. The pressure on the contact surface matched exactly with that obtained by Hahn and Levinson (1974). Stresses on the bond surface could not be compared since such information is not given in Hahn and Levinson's work for v =.3 and our formulation is not valid for v =.5. Next we solved the problem for a different set of material and geometric parameters, namely that RO = 61 mm, R = 76.2 min, Rl = 46.5 min, b = 12.7 mm, v =.45 and E = 1 x 16 N/m2. Batra (1981) has solved this problem by the finite element method. The pressure profile at the contact surface came out to be virtually the same by the two methods. For this case we plot, in Fig. 3, the variation of the normal stress rr and the shear stress ore at the bond surface with the angular distance from the center line of rollers. These stresses decay to essentially as e/eo approaches 3. Thus the assumption that the vertical plane is traction free is verified to be true. In this and all other figures eo = b/ro. Effect of thickness of the layer Keeping RO = 61 mm, R = 76.2 mm, b = 12.7 mm, v =.3 and E = 1 x 16 N/m2 fixed, the value of Rl was varied. The total load required and the!esultant indentation for various thicknesses of the layer are listed in Table 1. As expected, the load required to keep the contact width constant increases as the thickness of the layer decreases. For thinner layer the indentation is less even though the load is more. However, the indentation/thickness increases with the decrease in the thickness of the ].ayer. Thickness/Rn Load P!(Gb) Indentation/b Indentation Thickness CJ lQO ] Table 1. Load and indentation for different thicknesses of the layer Fig. 4 depicts the pressure profile at the contact surface for two different thicknesses of the layer. Note the change in the shape of the curve near the ends of the contact width as the thickness of the layer is reduced. With the decrease in the thickness of the layer it seems that the pressure at the center increases more than the decrease in the pressure near hp pnd so that the total load increases. Such a chane
7 265 in the curvature of the curve near the ends of the contact width is also apparent in Meijer's (1968) work which studied the indentation, by a rigid cylinder, of a layer fixed to a rigid plane surface. In Fig. 5 is shown the deformed surface of the rubberlike layer for two different thicknesses of the layer. Note that the scales along the vertical and horizontal axes are different. Therefore, the undeformed position of the rubberlike layer plots as an ellipse rather than as a circle. It seems that the change in the curvature of the deformed surface near the ends of the contact zone is heavily dependent upon the thickness of the rubberlike laver. Effect of Poisson's ratio To investigate the effect of Poisson's ratio on the load required to cause the same contact width the values of RO = 61 mm, R = 76.2 mm, Rl = 46.5 mm, b = 12.7 mm, and the shear modulus G = 3.45 x 15 N/m2 are kept fixed. The total load required and the resulting indentation for various values of v are listed in Table 2. These data show that Table 2. Load and indentation for various values of Poisson's ratio the total load increases with the increase in the value of Poisson's ratio upto a certain value of v. As v is increased beyond.45, the total load changes little but the pressure distribution on the contact surface differs in that it is more at the center and less near the ends. This becomes clear from Fig. 6 in which we plot the pressure distribution at the contact surface for two different values of v. Remarks We should emphasize that the results presented above and the conclusions drawn are applicable only to the geometric configuration studied herein. For other values of geometric parameters, one will expect similar qualitative and not necessarily quantitative results. For the present problem the CPU time and the core require-
8 266 ment for the finite element solution and the boundary element solution were virtually the same. In the boundary element solution of the problem the grid used had 1 elements across the thickness of the layer and 35 elements across the circumference. The boundary element solution perhaps gave slightly better values of tractions at the bond surface. The grid data in both caseswere generated internally in the computer program and required the same effort. Even though it is usually said that the finite element solution generates a lot of unneeded information about stresses and strains within a body, to us, it seems to be an advantage especially in a problem like the one studied here. In such a problem one does not know a priori where the maximum principal stress or the maximum principal strain will occur. It is more expensive to compute stresses and strains at internal points from the boundary element solution than it is from the finite element solution. REFERENCES Batra, R. C., Levin'3on. M and Be,:z, F. (19i6) Rubb"'r CCJve.red Rolls--The Thermoviscoelastic Problem. A Finite Element Solution, Int. J. Num. Meth. Engng. lq: Batra, R. C. (198) Rubber Covered Rol1s--The Nonlinear Elastic Problem. Trans. ASME - J. Appl. Mechs.. ]: Batra. R. C. (1981) Quasistatic Indentation of a Rubber Covered Roll by a Rigid Roll, Int. J. Num. Meth. Engng. li: (to appear). Brebbia. C. A. (1978) the Boundary Element Method for Engineers. John Wiley & Sons. New York. Hahn, H. T. and Levinson, M (1974) Indentation of an Elastic Layer Bonded to a Rigid Cylinder--I. Quasistatic Case Without Friction. II - Unidirectional Slipping with Coulomb Friction,.!..: Hannah, Margaret (1951) Contact Stress and Deformation in a Thin Elastic Layer, Quart. Journ. Mech. and Applied Math. i: Hertz, H. (1881) J. F. Math. 92: Meijers. P. (1968) The Contact Problem of a Rigid Cylinder on an Elastic Layer. Appl. Sci. Res. 1!: Spengos, A. C. (1965) Experimental Investigation of Rolling Contact, Trans. ASHE - J. Appl. Mechs., 11: Thomas, H. R. and Hoersch, V. A. (193) University of Ill.. Eng. Exp. Sta. Bull...'
9 I 1 1 Fig. 2 Boundary element grid "'" "' ":" "'
10
11 269 II) II)». u.c:... POi 'H U 'H. 'H " U (11.- It) (/) 'H N V..... II) II ". nt '"'. U CD CD.c: '"' I'-: Ct>1 '"'..... '"' II) -.:r V It) N -:..., co /oh..i) 3] / 3HnSS3Hd... II) II)....
12 27, ",'" "" ",,'" I") "" a> // a>; / " / c:n- I' I' WW / u '"' ). r-i /\ ) or! r-i '"' ).D.D I ) I.c I +J I I I w ) / u u I <t I 't:: I::> ie I (/) 'C I I W W,.. I u IOu cr <t I I u. crw (/)1'-- 11'\ W ::>w I 2: ::> " '., ", W'- bo 4' (/) LA- (\J (\J It) w m.. - It) o. >- cr
13 r om. j.j - cn '"' In In In om j.j u. Q) Q) C1)\ j.j U cn j.j. ) U V Q). j.j Q) u cn '"' ;1 In j.j cn (\J v om. j.j ='.c om. '"' j.j \ ) In \ om I'-." II' II.- ;1 (/) CD In In Q) '"' -D 1;.'. bo [(!!/O+1)9J/aJnSSaJd
Transactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN
On the Computation of Elastic Elastic Rolling Contact using Adaptive Finite Element Techniques B. Zastrau^, U. Nackenhorst*,J. Jarewski^ ^Institute of Mechanics and Informatics, Technical University Dresden,
More 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 informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More information2.1 Strain energy functions for incompressible materials
Chapter 2 Strain energy functions The aims of constitutive theories are to develop mathematical models for representing the real behavior of matter, to determine the material response and in general, to
More informationcircular voids On the interaction between two in a nonlinear elastic solid ACTA MECHANICA J. P. Zhang and R. C. Batra, Rolla, Missouri
Acta Mechanica 105, 161-171 (1994) ACTA MECHANICA 9 Springer-Verlag 1994 On the interaction between two in a nonlinear elastic solid circular voids J. P. Zhang and R. C. Batra, Rolla, Missouri (Received
More informationSteady state rolling contact of cylinders with viscoelastic layers M. Tervonen Department of Mechanical Engineering, University
Steady state rolling contact of cylinders with viscoelastic layers M. Tervonen Department of Mechanical Engineering, University ABSTRACT The finite element method is used to solve the steady state rolling
More informationA *69>H>N6 #DJGC6A DG C<>C::G>C<,8>:C8:H /DA 'D 2:6G - ( - ) +"' ( + -"( (' (& -+" % '('%"' +"-2 ( -!"',- % )% -.C>K:GH>IN D; AF69>HH>6,-+
The primary objective is to determine whether the structural efficiency of plates can be improved with variable thickness The large displacement analysis of steel plate with variable thickness at direction
More informationConservation 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 informationLinear Constitutive Relations in Isotropic Finite Viscoelasticity
Journal of Elasticity 55: 73 77, 1999. 1999 Kluwer Academic Publishers. Printed in the Netherlands. 73 Linear Constitutive Relations in Isotropic Finite Viscoelasticity R.C. BATRA and JANG-HORNG YU Department
More informationMechanical Properties of Materials
Mechanical Properties of Materials Strains Material Model Stresses Learning objectives Understand the qualitative and quantitative description of mechanical properties of materials. Learn the logic of
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 informationNew Representation of Bearings in LS-DYNA
13 th International LS-DYNA Users Conference Session: Aerospace New Representation of Bearings in LS-DYNA Kelly S. Carney Samuel A. Howard NASA Glenn Research Center, Cleveland, OH 44135 Brad A. Miller
More informationMaterials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation.
Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. In the process zone, stress amplitudes are poorly determined and much
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 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 informationPrediction of geometric dimensions for cold forgings using the finite element method
Journal of Materials Processing Technology 189 (2007) 459 465 Prediction of geometric dimensions for cold forgings using the finite element method B.Y. Jun a, S.M. Kang b, M.C. Lee c, R.H. Park b, M.S.
More informationThe force on a lattice defect in an elastic body
Journal of Elasticity 17 (1987) 3-8 @ Martinus Nijhoff Publishers, Dordrecht - Printed in the Netherlands The force on a lattice defect in an elastic body R.C. BATRA. " Department of Engineering Mechanics,
More informationExample-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 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 informationSTRAIN. Normal Strain: The elongation or contractions of a line segment per unit length is referred to as normal strain denoted by Greek symbol.
STRAIN In engineering the deformation of a body is specified using the concept of normal strain and shear strain whenever a force is applied to a body, it will tend to change the body s shape and size.
More informationLaminated Beams of Isotropic or Orthotropic Materials Subjected to Temperature Change
United States Department of Agriculture Forest Service Forest Products Laboratory Research Paper FPL 375 June 1980 Laminated Beams of Isotropic or Orthotropic Materials Subjected to Temperature Change
More informationConjugate gradient-boundary element method for a Cauchy problem in the Lam6 system
Conjugate gradient-boundary element method for a Cauchy problem in the Lam6 system L. Marin, D.N. Hho & D. Lesnic I Department of Applied Mathematics, University of Leeds, UK 'Hanoi Institute of Mathematics,
More informationIranian Journal of Mathematical Sciences and Informatics Vol.2, No.2 (2007), pp 1-16
Iranian Journal of Mathematical Sciences and Informatics Vol.2, No.2 (2007), pp 1-16 THE EFFECT OF PURE SHEAR ON THE REFLECTION OF PLANE WAVES AT THE BOUNDARY OF AN ELASTIC HALF-SPACE W. HUSSAIN DEPARTMENT
More informationChapter 3. Load and Stress Analysis. Lecture Slides
Lecture Slides Chapter 3 Load and Stress Analysis 2015 by McGraw Hill Education. This is proprietary material solely for authorized instructor use. Not authorized for sale or distribution in any manner.
More informationINTRODUCTION TO STRAIN
SIMPLE STRAIN INTRODUCTION TO STRAIN In general terms, Strain is a geometric quantity that measures the deformation of a body. There are two types of strain: normal strain: characterizes dimensional changes,
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 informationThe University of Melbourne Engineering Mechanics
The University of Melbourne 436-291 Engineering Mechanics Tutorial Four Poisson s Ratio and Axial Loading Part A (Introductory) 1. (Problem 9-22 from Hibbeler - Statics and Mechanics of Materials) A short
More informationRubber Elasticity (Indented text follows Strobl, other follows Doi)
Rubber Elasticity (Indented text follows Strobl, other follows Doi) Elasticity of A Single Chain: The spring constant associated with a single polymer chain is of importance in a wide range of situations
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 informationElastic-plastic deformation near the contact surface of the circular disk under high loading
Elastic-plastic deformation near the contact surface of the circular disk under high loading T. Sawada & M. Horiike Department of Mechanical Systems Engineering Tokyo University of Agriculture and Technology,
More informationNon-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 informationRadial Growth of a Micro-Void in a Class of. Compressible Hyperelastic Cylinder. Under an Axial Pre-Strain *
dv. Theor. ppl. Mech., Vol. 5, 2012, no. 6, 257-262 Radial Growth of a Micro-Void in a Class of Compressible Hyperelastic Cylinder Under an xial Pre-Strain * Yuxia Song, Datian Niu and Xuegang Yuan College
More information7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment
7.6 Stress in symmetrical elastic beam transmitting both shear force and bending moment à It is more difficult to obtain an exact solution to this problem since the presence of the shear force means that
More informationEMA 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 informationEXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD
Journal of KONES Powertrain and Transport, Vol. 7, No. EXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD Robert Czabanowski Wroclaw University
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 informationArbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact
Arbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact Philip P. Garland 1 and Robert J. Rogers 2 1 School of Biomedical Engineering, Dalhousie University, Canada 2 Department
More informationDepartment of Structural, Faculty of Civil Engineering, Architecture and Urban Design, State University of Campinas, Brazil
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm A SIMPLIFIED FORMULATION FOR STRESS AND TRACTION BOUNDARY IN- TEGRAL EQUATIONS USING THE
More informationTwo problems in finite elasticity
University of Wollongong Research Online University of Wollongong Thesis Collection 1954-2016 University of Wollongong Thesis Collections 2009 Two problems in finite elasticity Himanshuki Nilmini Padukka
More information16.20 Techniques of Structural Analysis and Design Spring Instructor: Raúl Radovitzky Aeronautics & Astronautics M.I.T
16.20 Techniques of Structural Analysis and Design Spring 2013 Instructor: Raúl Radovitzky Aeronautics & Astronautics M.I.T February 15, 2013 2 Contents 1 Stress and equilibrium 5 1.1 Internal forces and
More informationSTRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS
1 UNIT I STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS 1. Define: Stress When an external force acts on a body, it undergoes deformation. At the same time the body resists deformation. The
More informationNORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.
NORMAL STRESS The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts. σ = force/area = P/A where σ = the normal stress P = the centric
More informationFINAL EXAMINATION. (CE130-2 Mechanics of Materials)
UNIVERSITY OF CLIFORNI, ERKELEY FLL SEMESTER 001 FINL EXMINTION (CE130- Mechanics of Materials) Problem 1: (15 points) pinned -bar structure is shown in Figure 1. There is an external force, W = 5000N,
More informationModeling Contact between Rigid Sphere and Elastic Layer Bonded to Rigid Substrate
IEEE TRANSACTIONS ON COMPONENTS AND PACKAGING TECHNOLOGIES, VOL. 24, NO. 2, JUNE 2001 207 Modeling Contact between Rigid Sphere and Elastic Layer Bonded to Rigid Substrate Mirko Stevanović, M. Michael
More informationSolution: T, A1, A2, A3, L1, L2, L3, E1, E2, E3, P are known Five equations in five unknowns, F1, F2, F3, ua and va
ME 323 Examination # 1 February 18, 2016 Name (Print) (Last) (First) Instructor PROBLEM #1 (20 points) A structure is constructed from members 1, 2 and 3, with these members made up of the same material
More informationPhysics of Continuous media
Physics of Continuous media Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 October 26, 2012 Deformations of continuous media If a body is deformed, we say that the point which originally had
More informationSTRENGTH OF MATERIALS-I. Unit-1. Simple stresses and strains
STRENGTH OF MATERIALS-I Unit-1 Simple stresses and strains 1. What is the Principle of surveying 2. Define Magnetic, True & Arbitrary Meridians. 3. Mention different types of chains 4. Differentiate between
More informationAnalysis of forming- Slipline Field Method
Analysis of forming- Slipline Field Method R. Chandramouli Associate Dean-Research SASTRA University, Thanjavur-613 401 Joint Initiative of IITs and IISc Funded by MHRD Page 1 of 7 Table of Contents 1.
More informationEngineering Sciences 241 Advanced Elasticity, Spring Distributed Thursday 8 February.
Engineering Sciences 241 Advanced Elasticity, Spring 2001 J. R. Rice Homework Problems / Class Notes Mechanics of finite deformation (list of references at end) Distributed Thursday 8 February. Problems
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 informationUNIT 1 STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS 1. Define stress. When an external force acts on a body, it undergoes deformation.
UNIT 1 STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS 1. Define stress. When an external force acts on a body, it undergoes deformation. At the same time the body resists deformation. The magnitude
More informationChapter 3. Load and Stress Analysis
Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3
More informationVORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.
More informationPDDC 1 st Semester Civil Engineering Department Assignments of Mechanics of Solids [ ] Introduction, Fundamentals of Statics
Page1 PDDC 1 st Semester Civil Engineering Department Assignments of Mechanics of Solids [2910601] Introduction, Fundamentals of Statics 1. Differentiate between Scalar and Vector quantity. Write S.I.
More informationBack Analysis of Measured Displacements of Tunnels
Rock Mechanics and Rock Engineering 16, 173--180 (1983) Rock Mechanics and Rock Engineering 9 by Springer-Verlag 1983 Back Analysis of Measured Displacements of Tunnels By S. Sakurai and K. Takeuchi Kobe
More informationA multiscale framework for lubrication analysis of bearings with textured surface
A multiscale framework for lubrication analysis of bearings with textured surface *Leiming Gao 1), Gregory de Boer 2) and Rob Hewson 3) 1), 3) Aeronautics Department, Imperial College London, London, SW7
More informationGlossary. Glossary of Symbols. Glossary of Roman Symbols Glossary of Greek Symbols. Contents:
Glossary Glossary of Symbols Contents: Glossary of Roman Symbols Glossary of Greek Symbols Glossary G-l Glossary of Roman Symbols The Euclidean norm or "two-norm." For a vector a The Mooney-Rivlin material
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 informationOptimization of blank dimensions to reduce springback in the flexforming process
Journal of Materials Processing Technology 146 (2004) 28 34 Optimization of blank dimensions to reduce springback in the flexforming process Hariharasudhan Palaniswamy, Gracious Ngaile, Taylan Altan ERC
More informationGG303 Lab 12 11/7/18 1
GG303 Lab 12 11/7/18 1 DEFORMATION AROUND A HOLE This lab has two main objectives. The first is to develop insight into the displacement, stress, and strain fields around a hole in a sheet under an approximately
More informationMITOCW 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 informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More informationSTATICALLY INDETERMINATE STRUCTURES
STATICALLY INDETERMINATE STRUCTURES INTRODUCTION Generally the trusses are supported on (i) a hinged support and (ii) a roller support. The reaction components of a hinged support are two (in horizontal
More informationSymmetric Bending of Beams
Symmetric Bending of Beams beam is any long structural member on which loads act perpendicular to the longitudinal axis. Learning objectives Understand the theory, its limitations and its applications
More informationLINEAR AND NONLINEAR SHELL THEORY. Contents
LINEAR AND NONLINEAR SHELL THEORY Contents Strain-displacement relations for nonlinear shell theory Approximate strain-displacement relations: Linear theory Small strain theory Small strains & moderate
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationPlates and Shells: Theory and Computation. Dr. Mostafa Ranjbar
Plates and Shells: Theory and Computation Dr. Mostafa Ranjbar Outline -1-! This part of the module consists of seven lectures and will focus on finite elements for beams, plates and shells. More specifically,
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 informationStructural Analysis I Chapter 4 - Torsion TORSION
ORSION orsional stress results from the action of torsional or twisting moments acting about the longitudinal axis of a shaft. he effect of the application of a torsional moment, combined with appropriate
More informationA 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 informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
More informationGraduate School of Engineering, Kyoto University, Kyoto daigaku-katsura, Nishikyo-ku, Kyoto, Japan.
On relationship between contact surface rigidity and harmonic generation behavior in composite materials with mechanical nonlinearity at fiber-matrix interface (Singapore November 2017) N. Matsuda, K.
More informationME Final Exam. PROBLEM NO. 4 Part A (2 points max.) M (x) y. z (neutral axis) beam cross-sec+on. 20 kip ft. 0.2 ft. 10 ft. 0.1 ft.
ME 323 - Final Exam Name December 15, 2015 Instructor (circle) PROEM NO. 4 Part A (2 points max.) Krousgrill 11:30AM-12:20PM Ghosh 2:30-3:20PM Gonzalez 12:30-1:20PM Zhao 4:30-5:20PM M (x) y 20 kip ft 0.2
More information4. BEAMS: CURVED, COMPOSITE, UNSYMMETRICAL
4. BEMS: CURVED, COMPOSITE, UNSYMMETRICL Discussions of beams in bending are usually limited to beams with at least one longitudinal plane of symmetry with the load applied in the plane of symmetry or
More informationROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring
ROTATIONAL STRESSES INTRODUCTION High centrifugal forces are developed in machine components rotating at a high angular speed of the order of 100 to 500 revolutions per second (rps). High centrifugal force
More informationCOMPRESSION AND BENDING STIFFNESS OF FIBER-REINFORCED ELASTOMERIC BEARINGS. Abstract. Introduction
COMPRESSION AND BENDING STIFFNESS OF FIBER-REINFORCED ELASTOMERIC BEARINGS Hsiang-Chuan Tsai, National Taiwan University of Science and Technology, Taipei, Taiwan James M. Kelly, University of California,
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 informationMaterial tailoring and moduli homogenization for finite twisting deformations of functionally graded Mooney-Rivlin hollow cylinders
Acta Mech 224, 811 818 (2013) DOI 10.1007/s00707-012-0784-z J. R. Dryden R. C. Batra Material tailoring and moduli homogenization for finite twisting deformations of functionally graded Mooney-Rivlin hollow
More informationAnalysis of Frictional Torque in Raceway Contacts of Tapered Roller Bearings
Analysis of Frictional Torque in Raceway Contacts of Tapered Roller Bearings H. MATSUYAMA * S. KAMAMOTO ** * Bearing Research & Development Department, Research & Development Center **Mechatronic Systems
More informationDetermination of Mechanical Properties of Elastomers Using Instrumented Indentation
Determination of Mechanical Properties of Elastomers Using Instrumented Indentation, Antonios E. Giannakopoulos and Dimitrios Bourntenas University of Thessaly, Department of Civil Engineering, Volos 38334,
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 informationFinite 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 informationChapter 3 Stress, Strain, Virtual Power and Conservation Principles
Chapter 3 Stress, Strain, irtual Power and Conservation Principles 1 Introduction Stress and strain are key concepts in the analytical characterization of the mechanical state of a solid body. While stress
More informationThe plastic behaviour of silicon subjected to micro-indentation
JOURNAL OF MATERIALS SCIENCE 31 (1996) 5671-5676 The plastic behaviour of silicon subjected to micro-indentation L. ZHANG, M. MAHDI Centre for Advanced Materials Technology, Department of Mechanical and
More informationNumerical and experimental analysis of a cantilever beam: A laboratory project to introduce geometric nonlinearity in Mechanics of Materials
Numerical and experimental analysis of a cantilever beam: A laboratory project to introduce geometric nonlinearity in Mechanics of Materials Tarsicio Beléndez (1) and Augusto Beléndez (2) (1) Departamento
More informationANSYS Mechanical Basic Structural Nonlinearities
Lecture 4 Rate Independent Plasticity ANSYS Mechanical Basic Structural Nonlinearities 1 Chapter Overview The following will be covered in this Chapter: A. Background Elasticity/Plasticity B. Yield Criteria
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 informationExample 3.7 Consider the undeformed configuration of a solid as shown in Figure 3.60.
162 3. The linear 3-D elasticity mathematical model The 3-D elasticity model is of great importance, since it is our highest order hierarchical model assuming linear elastic behavior. Therefore, it provides
More informationDEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm DEFORMATION AND FRACTURE ANALYSIS OF ELASTIC SOLIDS BASED ON A PARTICLE METHOD R. A. Amaro
More informationA parametric study on the elastic-plastic deformation of a centrally heated two-layered composite cylinder with free ends
Arch. Mech., 68, 3, pp. 03 8, Warszawa 06 A parametric study on the elastic-plastic deformation of a centrally heated two-layered composite cylinder with free ends F. YALCIN ), A. OZTURK ), M. GULGEC 3)
More informationA FAILURE CRITERION FOR POLYMERS AND SOFT BIOLOGICAL MATERIALS
Material Technology A FALURE CRTERON FOR POLYMERS AND SOFT BOLOGCAL MATERALS Authors: William W. Feng John O. Hallquist Livermore Software Technology Corp. 7374 Las Positas Road Livermore, CA 94550 USA
More informationHardened Concrete. Lecture No. 16
Hardened Concrete Lecture No. 16 Fatigue strength of concrete Modulus of elasticity, Creep Shrinkage of concrete Stress-Strain Plot of Concrete At stress below 30% of ultimate strength, the transition
More informationAn example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction
An example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction Un exemple de non-unicité pour le modèle continu statique de contact unilatéral avec frottement de Coulomb
More informationD : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each.
GTE 2016 Q. 1 Q. 9 carry one mark each. D : SOLID MECHNICS Q.1 single degree of freedom vibrating system has mass of 5 kg, stiffness of 500 N/m and damping coefficient of 100 N-s/m. To make the system
More informationCHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 2. Discontinuity functions
1. Deflections of Beams and Shafts CHAPTER OBJECTIVES Use various methods to determine the deflection and slope at specific pts on beams and shafts: 1. Integration method. Discontinuity functions 3. Method
More informationIntegral equations for crack systems in a slightly heterogeneous elastic medium
Boundary Elements and Other Mesh Reduction Methods XXXII 65 Integral equations for crack systems in a slightly heterogeneous elastic medium A. N. Galybin & S. M. Aizikovich Wessex Institute of Technology,
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationelastoplastic contact problems D. Martin and M.H. Aliabadi Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK
Non-conforming BEM elastoplastic contact problems D. Martin and M.H. Aliabadi discretisation in Wessex Institute of Technology, Ashurst Lodge, Ashurst, Southampton, SO40 7AA, UK Abstract In this paper,
More informationHIGHLY ADAPTABLE RUBBER ISOLATION SYSTEMS
th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 24 Paper No. 746 HIGHLY ADAPTABLE RUBBER ISOLATION SYSTEMS Luis DORFMANN, Maria Gabriella CASTELLANO 2, Stefan L. BURTSCHER,
More informationLongitudinal buckling of slender pressurised tubes
Fluid Structure Interaction VII 133 Longitudinal buckling of slender pressurised tubes S. Syngellakis Wesse Institute of Technology, UK Abstract This paper is concerned with Euler buckling of long slender
More information