Transactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN
|
|
- Wilfrid Wade
- 6 years ago
- Views:
Transcription
1 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, D Dresden, Germany urz. tu-dresden. de ^Institute of Mechanics, University of the Federal Armed Forces Hamburg, D Hamburg, Germany de Abstract The computation of the stress and strain of a rolling wheel on the rail requires the thorough analysis of the contact condition including slip. Prom the physical point of view it is important to use an appropriate description of the surface roughness and to take account of the rolling condition because this differs significantly from static contact. As afirstapproach a two dimensional model of a wheel running with constant velocity over an infinite elastic strip is investigated. Because of the necessity to use a dense mesh in the contact area and in order to avoid a transient computation it is advantageous to use an arbitrary Lagrangian Eulerian description (ALE). This paper contains a description of the ALE approach for the problem of the elastic-elastic contact of a wheel rolling with constant velocity over a rail of infinite length. The roughness of the surface is used in the formulation of a physically motivated penalty formulation and is exploited for the normal contact condition. The approach for tangential contact of rough surfaces according to Willners model is stated for completeness. A concept of an error estimator on the basis of the surface traction error as a criteria for the refinement of the mesh in the contact area is presented. As a numerical example the solution for a roller is presented in order to demonstrate the practicability of the introduced approach.
2 130 Contact Mechanics HI 1 Introduction The analysis of the contact behavior of wheels is of great interest because of the significance of the contact condition, e.g. for sliding of vehicles and for the prediction of wear. Since static contact and dynamic contact, especially rolling contact differ a lot, special emphasis is laid on the treatment of rolling contact. The main goal of the herewith published preliminary results of this research is the computation of stresses between the wheel and the elastic rail in typical railway applications. Although the real rolling on a rail is a transient process, in a first approach only the steady state motion of rolling of a wheel on an infinite elastic strip is solved. This allows the direct application of the ALE-method. The computation of contact problems is highly nonlinear. There exist only very few analytical solutions for two-dimensional problems like Hertzian contact and steady state rolling of homogeneous cylinders (Bufler [1]) against each other. Due to the three-dimensional character according to the shape of the wheel and the rail a numerical solution has to be obtained. The computation of contact is often done by application of the penalty method with a numerical suitable penalty parameter. In order to treat realistic surfaces it is natural to apply Willners approach [4] exploiting statistical information about the roughness of a surface thus leading to a nonlinear constitutive description for normal and tangential contact. Both laws will be stated here and the computation of the contact pressure is done by this rule. Since the main concern is to use a correct rolling contact condition for elastic-elastic contact, this will be described in more detail and the numerical results presented here are confined to sticking of wheel and rail. Nevertheless it is advantageous to introduce a local slip velocity and a rolling contact penalty parameter to compute the tangential forces. Some basics for error measures of surface traction errors and adaptive meshing are discussed. The presentation of numerical results will be postponed until the effectiveness of the proposed method can be proved, e.g. in comparison to the Z^-method. 2 Basic Equations of Motion For a pair of two contacting bodies Hamiltons Principle reads as 2 / ( (^ -6W' + I 0B' -rjdv + I T' T] da) J \ T _ -I J J ' *i ^~* B' d(b') da dt = 0. (2.1)
3 Contact Mechanics HI 131 In this expression 6T= f Q (B) is the variation of the kinetic energy and dt, dt dv (2.2) 6W = I P -GRADr? dv (B) (2.3) is the variation of the internal energy expressed in terms of the second Piola Kirchhoff stress tensor P. B^ and T* are the body forces and the surface loads acting at body / (/ = 1,2) respectively, while 77 is a test function (virtual displacement). The second term in equation (2.1) describes the virtual work of the contact forces, TN is the normal contact stress, 5g is the variation of the gap, TT is the vector of the tangent contact stresses and 8s is the variation of the slip vector. The integration is performed over the contact zone dc(b). For an efficient numerical solution of the rolling contact problem an Arbitrary Lagrangian Eulerian (ALE) description of the contacting bodies is chosen according to Nackenhorst [2]. Basic idea of this approach is the decomposition into a the rigid body motion and the deformation, which allows to compute a steady state rolling motion without integration with respect to time. Fig. 2.1 gives a scheme of the mathematical idea of the decomposition. The decomposed total motion \& is introduced as the rigid body motion % composed with the deformation c/>, which is measured with respect to the moving reference configuration %. initial configuration deformed configuration reference configuration Figure 2.1: ALE-decomposition into rigid body motion and deformation. Within this kinematical frame the material derivative is defined as d(---) dt dt +GRAB( (2.4)
4 132 Contact Mechanics HI where % describes therigidbody velocity of a spinning wheel as Thus, for example the material velocity reads as <9y X = -^r + n x %. (2.5) y = ^ + GRAD<f>.%. (2.6) Using this, the variation of the kinetic energy can be rewritten as, ' X) (^ + GRAD?7 %) dl/. (2.7) Some final basic transformations lead to the following expression y 6T df = - y ( y <X^ T? - y GRAD7? x) dl/) d^, (2.8) ti ti (B) which provides symmetric matrices for the subsequent finite element approximation. It has been found to be advantageous to introduce a coordinate system following the axis of the wheel, which is moving with constant velocity. Thus it is assumed that a wheel is rotating around a fixed axis, whereas the road or rail is moving under it horizontally. Within this frame the rigid body motion is described for the wheel as X = ^xx (2.9) with the angular velocity ffc and for the road or rail according to X = % (2.10) The advantage of the ALE description is combined with higher numerical effort for the formulation of the sticking condition. From equation (2.6) obviously follows that the sticking condition for rolling contact (2), i.e. the vanishing relative motion of sticking particles, is not derived automatically from the identical displacement because the displacement gradient of the contacting points is not necessarily the same. Thus an additional penalty term a <;2 do, (2.12) is added to the left hand side of equation (2.1) to impose the sticking condition of rolling surfaces.
5 Contact Mechanics HI 3 Constitutive Laws of Contact Interfaces According to Willner [4] a statistical description of the surface roughness is used to derive constitutive laws of contacting surfaces. For each asperity Hertzian contact is assumed in normal direction and the hypothesis of maximum shear in tangential direction is used. By use of a probabilistic method constitutive laws of the contact interface are as follows 7 0 oo oo 7 0 d( (3.1) (3.2) Herein 7, % and ( are the gap, curvature and tip height of an asperity scaled by the statistically derived root mean square values of the contacting surfaces respectively, s is the slip vector, N is the normal contact force of one asperity with the gap 7 and the tip curvature %, T is the tangential force of an asperity with the gap 7 and the slip s and P is the probability that a tip with gap 7 and curvature K, is in contact. These laws of contact interface behavior are illustrated in fig. 3.1 for a typical steel to steel contact "S ts3 13 So normalized 3 4 gap y Figure 3.1: Pressure-gap function (left) and contact shear stress as a function of pressure and slip (Willner [4]). 4 Finite Element Approach The description of the rolling contact problem by the equations (2.1) in addition with the constitutive interface laws (3.1), (3.2) and nonlinear material behavior is a highly nonlinear problem. Thus, for the numerical solution an
6 134 Contact Mechanics HI incremental approach is necessary wherein the equations are linearized consistently. After linearization a finite element approximation of the problem is performed by common techniques, which lead to quasi static incremental equilibrium equations for the steady state rolling process, * + <fc+ t- t. (4.1) In this equation *K is the tangential stiffness matrix, W is the inertia matrix due to the motion, *Kc is the contact stiffness matrix, *Ks is the contact stiffness generated from the linearization of the ALE-sticking conditions, *~^*fg is the matrix of the applied body and surface loads, * is the matrix of the inertia forces, % is the matrix of the contact forces, *f«is the matrix representing the non-equilibrium of the sticking conditions and % is the matrix of the internal force representing the actual stress state. 4.1 Error Estimation and Adaptivity The reliability of finite element results is of great interest especially in applications where no or only few experimental verifications are possible. Error estimation and adaptive mesh refinement methods are developed since more than one decade. In present the adaptive technique has become a common tool in commercial finite element codes, most of them are based on the well accepted Zienkiewicz/Zhu (7?) method [6]. As shown in [3] for critical boundary value problems such as contact analysis this method is not very suitable because the boundary conditions are not emphasized explicitly. Another drawback of commercial FE-codes in general originates from the fact that the common mesh refinement algorithms are not able to deal with arbitrary element types. Especially in three-dimensional computations often these algorithms support thetrahedrical elements only. For alternative adaptivefiniteelement techniques with special emphasis to contact analysis it is referred to [5] and the literature cited there. Because of the aforementioned reasons it is advantageous to follow the concept suggested by Nackenhorst [3]. Basic idea of this approximation is to apply mesh refinement only to those parts of the surface of special interest, namely the contact region. For this an approximation error in the surface tractions is given by \e\c= r, -Te da, (4.2) where TC is a continuous approximation of the surface tractions through equations (3.1) and (3.2) and TC is the representation of the surface tractions by commonfiniteelement solutions.
7 Contact Mechanics HI Numerical Example In order to demonstrate the practicability of the proposed method an elastic roller on an elastic strip has been analyzed. The roller is turning with constant velocity under free rolling conditions. The surface has been assumed to be rough in the sense of Willner and in the contact zone no sliding shall occur. The adopted discretization is according to the ALE-approach dense in the contact area only and coarse in the distant parts of both bodies. The chosen mesh is shown in fig. 5.1 for a deformed state with vertical movement of the axis of 0.04 R- The angular velocity is such small that the exploitation of the rolling contact condition is sufficient while dynamical effects are negligeable. 5R Figure 5.1: Deformed mesh of a rolling elasic wheel on an elastic rail. The plot of the overall deformed configuration shows no difference to Hertzian contact, although in the contact area itself two significant deviations are found. One difference is the development of a smooth boundary of the contact area because of the adopted model of rough normal contact. The appropriate normal stress distribution is shown in fig. 5.2 in a normalized form, since a quantitative comparison with the Hertzian contact is meaningless because of the rigid core of the wheel and the small thickness of the rail. This rigid support leads to development of a broader zone of high pressure at approximately 90% of the maximum value of the static Hertzian problem. And because of the magnitude of indention the width of the contact zone is greater than the value a as being computed according to Hertz. The contact shear stresses are very small thus having little influence on the shear distribution in the inside of the wheel.
8 136 Contact Mechanics III normalized coordinate -J d. 1.0 Figure 5.2: Normalized contact pressure and shear stress distribution of a rolling wheel. maximum shear Tmax Anax 0,09 Q.Q9 Figure 5.3: Maximum shear stress distribution of a rolling wheel on an elastic strip.
9 Contact Mechanics HI 137 O.l-i -O.H Figure 5.4: Development of the normalized shear stress distribution of a rolling wheel due to variation of the ratio of Youngs moduli. Fig. 5.3 presents the distribution of the maximum shear stress inside both bodies. The maximum strain occurs as known from the circular disk inside. Due to the formerly discussed boundary conditions its value of 0.45 Pmax is approximately 30% higher than the value of the corresponding Hertzian solution. The occurrence and distribution of contact shear stresses of rolling cylinders is according to Bufler [1] dependent on the material parameters of the contacting bodies. For cylinders with identical Youngs moduli and Poissons ratio no contact shear stresses will occur. The underlying assumption of Bufler is, that the contact zone is small with respect to the total measures of the body, thus leading to an equivalent stiffness of the contacting cylinders. Both the rigid inner cylinder of the wheel and the finite thickness of the rail lead to different stiffness of the contacting bodies and thus to a shear stress distribution even in the case of identical material parameters. The variation of the Youngs modulus of the contacting bodies changes the maximum shear stress significantly and finally leads to a change of the direction of the tangential forces in the contact area. Nevertheless the maximum value remains small in comparison to the case of elastic-rigid contact. 6 Conclusions Centrifugal forces do not influence the contact behavior significantly at realistic speed, but the rolling contact conditions have to be taken into account
10 138 Contact Mechanics HI accurately. The computation of the slip condition is of central importance since the tangential forces directly correspond to the slip velocity. Further investigation will be necessary for the development of constitutive laws for friction. Especially the transition from sticking to sliding has to be investigated for the model of a rough surface. As known from Nackenhorst in this context the exploitation of adaptive remeshing as theoretically discussed in chapter 4 will be necessary. Acknowledgment This work is supported by the DFG (Deutsche Forschungsgemeinschaft) under contract Z141/5-1, which is gratefully acknowledged. Key-Words elastic-elastic rolling contact, adaptive finite element method, constitutive contact laws, arbitrary Lagrangian Eulerian (ALE) description References [I] Bufler, H. Zur Theorie der rollenden Reibung, Ingenieur-Archiv, 1959, 27, [2] Nackenhorst, U. On the finite element analysis of steady state rolling contact, in: Contact Mechanics I (eds M.H. Aliabadi & C. A. Brebbia), pp , Proceedings of the 1st International Conference on Computational Methods in Contact Mechanics, CMP, Southampton, [3] Nackenhorst, U. An adaptivefiniteelement method to analyze contact problems, in: Contact Mechanics II (eds M.H. Aliabadi & C. Alessandri), pp , Proceedings of the 2nd International Conference on Computational Methods in Contact Mechanics, CMP, Southampton, [4] Willner, K. & Gaul, L. A penalty approach for contact description by FEM based on interface physics, in: Contact Mechanics II (eds M.H. Aliabadi & C. Alessandri), pp , Proceedings of the 2nd International Conference on Computational Methods in Contact Mechanics, CMP, Southampton, [5] Wriggers, P. Finite element algorithms for contact problems, Archives of Computational Methods in Engineering, 1995, 2,4, [6] Zienkiewicz, O.C. & Zhu, J.Z. A simple error estimator and adaptive procedure for practical engineering analysis, Int. J. Num. Meth. Engg., 1987, 24,
Measurement 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 information5.5 Exercises for This Chapter Two-Axle Vehicle on Cosine Track Two-Axle Vehicle on Generally Periodic Track...
Contents 1 Introduction... 1 1.1 The Basic Function of the Wheel/rail System.... 1 1.2 Significance of Dynamics on the Operation of Rail Vehicles... 2 1.3 On the History of Research in the Field of Railway
More informationCh. 10: Fundamental of contact between solids
Ch. 10: Fundamental of contact between solids Actual surface is not smooth. At atomic scale, there are always defects at surface, such as vacancies, ledges, kinks, terraces. In micro or macro scale, roughness
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 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 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 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 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 informationNonlinear Theory of Elasticity. Dr.-Ing. Martin Ruess
Nonlinear Theory of Elasticity Dr.-Ing. Martin Ruess geometry description Cartesian global coordinate system with base vectors of the Euclidian space orthonormal basis origin O point P domain of a deformable
More informationFinite Element Analysis of Rolling Tires A State of the Art Review
Finite Element Analysis of Rolling Tires A State of the Art Review Maik Brinkmeier, Udo Nackenhorst and Matthias Ziefle Institut für Baumechanik und Numerische Mechanik Gottfried Wilhelm Leibniz Universität
More informationTransactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN
Simulation of wear with a FE tyre model using a steady state rolling formulation A. Becker, B. Seifert Tyre Research, Tyre Mechanics, Continental AG PO box 169, D-30001 Hannover, Germany EMail: abecker@conti.de,
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationDiscrete Element Modelling of a Reinforced Concrete Structure
Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr
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 informationKey words: Polymeric Composite Bearing, Clearance, FEM
A study on the effect of the clearance on the contact stresses and kinematics of polymeric composite journal bearings under reciprocating sliding conditions Abstract The effect of the clearance on the
More informationPLANAR KINETIC EQUATIONS OF MOTION (Section 17.2)
PLANAR KINETIC EQUATIONS OF MOTION (Section 17.2) We will limit our study of planar kinetics to rigid bodies that are symmetric with respect to a fixed reference plane. As discussed in Chapter 16, when
More informationEarthquake response analysis of rock-fall models by discontinuous deformation analysis
c Earthquake response analysis of rock-fall models by discontinuous deformation analysis T. Sasaki, I. Hagiwara & K. Sasaki Rock Engineering Lab., Suncoh Consultants Co. Ltd., Tokyo, Japan R. Yoshinaka
More informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationInstabilities and Dynamic Rupture in a Frictional Interface
Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/
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 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 informationComputational non-linear structural dynamics and energy-momentum integration schemes
icccbe 2010 Nottingham University Press Proceedings of the International Conference on Computing in Civil and Building Engineering W Tizani (Editor) Computational non-linear structural dynamics and energy-momentum
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 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 informationThe Subsurface Crack Under Conditions of Slip and Stick Caused by a Surface Normal Force
F.-K.Chang Maria Comninou Mem. ASME Sheri Sheppard Student Mem. ASME J. R. Barber Department of Mechanical Engineering and Applied Mechanics, University of Michigan, Ann Arbor, Mich. 48109 The Subsurface
More informationFinite element analysis of indentation experiments J.M. Olaf Fraunhofer-Insitut fur Werkstoffmechanik, Wohlerstr. 11, D-W Freiburg, Germany
Finite element analysis of indentation experiments J.M. Olaf Fraunhofer-Insitut fur Werkstoffmechanik, Wohlerstr. 11, D-W- 7800 Freiburg, Germany ABSTRACT There are only a few methods suitable for a quantitative
More informationLIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS
5 th GRACM International Congress on Computational Mechanics Limassol, 29 June 1 July, 2005 LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS G.A. Drosopoulos I, G.E.
More 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 informationChapter 12. Static Equilibrium and Elasticity
Chapter 12 Static Equilibrium and Elasticity Static Equilibrium Equilibrium implies that the object moves with both constant velocity and constant angular velocity relative to an observer in an inertial
More informationPrediction of the bilinear stress-strain curve of engineering material by nanoindentation test
Prediction of the bilinear stress-strain curve of engineering material by nanoindentation test T.S. Yang, T.H. Fang, C.T. Kawn, G.L. Ke, S.Y. Chang Institute of Mechanical & Electro-Mechanical Engineering,
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 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 informationModeling of Contact Interfaces in Built-up Structures by Zero-thickness Elements
Modeling of Contact Interfaces in Built-up Structures by Zero-thickness Elements Professor Lothar Gaul 1, Dr. Michael Mayer 2 1 University of Stuttgart, Institute of Applied and Experimental Mechanics,
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 informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationQUASISTATIC INDENTATION OF A RUBBER COVERED ROLL BY A RIGID ROLL - THE BOUNDARY ELEMENT SOLUTION
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
More informationAN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS
th Pan-American Congress of Applied Mechanics January 04-08, 00, Foz do Iguaçu, PR, Brazil AN ALTERNATIVE TECHNIQUE FOR TANGENTIAL STRESS CALCULATION IN DISCONTINUOUS BOUNDARY ELEMENTS Otávio Augusto Alves
More informationinertia of a body, principal axes of inertia, invariants of an inertia tensor, and inertia triangle inequalities are illustrated and discussed.
Preface This book belongs to a series of three books written simultaneously (the remaining two are titled Classical Mechanics: Dynamics and Classical Mechanics: Applied MechanicsandMechatronics). This
More informationFigure 43. Some common mechanical systems involving contact.
33 Demonstration: experimental surface measurement ADE PhaseShift Whitelight Interferometer Surface measurement Surface characterization - Probability density function - Statistical analyses - Autocorrelation
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 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 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 informationOther state variables include the temperature, θ, and the entropy, S, which are defined below.
Chapter 3 Thermodynamics In order to complete the formulation we need to express the stress tensor T and the heat-flux vector q in terms of other variables. These expressions are known as constitutive
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 informationMeasurement and modelling of contact stiffness
Measurement and modelling of contact stiffness D. Nowell, University of Oxford, UK Joints workshop, Chicago, 16/8/12 Difficulties in modelling contacts In general, the normal and tangential stiffnesses
More informationTemperature analysis of a pin-on-disc tribology test using experimental and numerical approaches
Friction 4(2): 135 143 (2016) ISSN 2223-7690 DOI 10.1007/s40544-016-0110-1 CN 10-1237/TH RESEARCH ARTICLE Temperature analysis of a pin-on-disc tribology test using experimental and numerical approaches
More informationEffect of Transient Dynamic Loading on Flexible Pavement Response
Effect of Transient Dynamic Loading on Flexible Pavement Response Imad L. Al-Qadi Founder Professor of Engineering Pyeong Jun Yoo Graduate Research Assistant Illinois Center for Transportation University
More informationEffect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study
Effect of Strain Hardening on Unloading of a Deformable Sphere Loaded against a Rigid Flat A Finite Element Study Biplab Chatterjee, Prasanta Sahoo 1 Department of Mechanical Engineering, Jadavpur University
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationDiscrete Analysis for Plate Bending Problems by Using Hybrid-type Penalty Method
131 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 21 (2008) Published online (http://hdl.handle.net/10114/1532) Discrete Analysis for Plate Bending Problems by Using
More informationPhysics. Assignment-1(UNITS AND MEASUREMENT)
Assignment-1(UNITS AND MEASUREMENT) 1. Define physical quantity and write steps for measurement. 2. What are fundamental units and derived units? 3. List the seven basic and two supplementary physical
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 informationBoundary Conditions in Fluid Mechanics
Boundary Conditions in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University The governing equations for the velocity and pressure fields are partial
More informationDEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS
DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the
More informationFALL TERM EXAM, PHYS 1211, INTRODUCTORY PHYSICS I Saturday, 14 December 2013, 1PM to 4 PM, AT 1003
FALL TERM EXAM, PHYS 1211, INTRODUCTORY PHYSICS I Saturday, 14 December 2013, 1PM to 4 PM, AT 1003 NAME: STUDENT ID: INSTRUCTION 1. This exam booklet has 14 pages. Make sure none are missing 2. There is
More informationAn Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation
An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation Nachiket Patil, Deepankar Pal and Brent E. Stucker Industrial Engineering, University
More 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 informationTolerance Ring Improvement for Reducing Metal Scratch
International Journal of Scientific and Research Publications, Volume 2, Issue 11, November 2012 1 Tolerance Ring Improvement for Reducing Metal Scratch Pattaraweerin Woraratsoontorn*, Pitikhate Sooraksa**
More informationDYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION
October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,
More 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 informationAN INTRODUCTION TO LAGRANGE EQUATIONS. Professor J. Kim Vandiver October 28, 2016
AN INTRODUCTION TO LAGRANGE EQUATIONS Professor J. Kim Vandiver October 28, 2016 kimv@mit.edu 1.0 INTRODUCTION This paper is intended as a minimal introduction to the application of Lagrange equations
More informationANALYTICAL MODEL FOR FRICTION FORCE BETWEEN A STEEL ROLLER AND A PLANE POLYMER SAMPLE IN SLIDING MOTION
6 th International Conference Computational Mechanics and Virtual Engineering COMEC 015 15-16 October 015, Braşov, Romania ANALYTICAL MODEL FOR FRICTION FORCE BETWEEN A STEEL ROLLER AND A PLANE POLYMER
More informationAn Atomistic-based Cohesive Zone Model for Quasi-continua
An Atomistic-based Cohesive Zone Model for Quasi-continua By Xiaowei Zeng and Shaofan Li Department of Civil and Environmental Engineering, University of California, Berkeley, CA94720, USA Extended Abstract
More informationLecture 2: Stresses in Pavements
Lecture 2: Stresses in Pavements Stresses in Layered Systems At any point, 9 stresses exist. They are 3 normal stresses (s z, s r, s t ) and 6 shearing stresses ( t rz = t zr, t rt = t tr, and t tz = t
More informationA Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against a Rigid Flat under Varying Modulus of Elasticity and Sphere Radius
Engineering, 2010, 2, 205-211 doi:10.4236/eng.2010.24030 Published Online April 2010 (http://www. SciRP.org/journal/eng) 205 A Finite Element Study of Elastic-Plastic Hemispherical Contact Behavior against
More information17th European Conference on Fracture 2-5 September,2008, Brno, Czech Republic. Thermal Fracture of a FGM/Homogeneous Bimaterial with Defects
-5 September,8, Brno, Czech Republic Thermal Fracture of a FGM/Homogeneous Bimaterial with Defects Vera Petrova, a, Siegfried Schmauder,b Voronezh State University, University Sq., Voronezh 3946, Russia
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationANALYTICAL PENDULUM METHOD USED TO PREDICT THE ROLLOVER BEHAVIOR OF A BODY STRUCTURE
The 3rd International Conference on Computational Mechanics and Virtual Engineering COMEC 2009 29 30 OCTOBER 2009, Brasov, Romania ANALYTICAL PENDULUM METHOD USED TO PREDICT THE ROLLOVER BEHAVIOR OF A
More informationEFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A SPHERE AGAINST A RIGID FLAT A FINITE ELEMENT STUDY
Proceedings of the International Conference on Mechanical Engineering 2009 (ICME2009) 26-28 December 2009, Dhaka, Bangladesh ICME09- EFFECT OF STRAIN HARDENING ON ELASTIC-PLASTIC CONTACT BEHAVIOUR OF A
More informationAbstract. 1 Introduction
Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut
More informationThermal Stress Analysis of BOF Accounting for Friction Force Experimentally Measured
, pp. 1125 1133 Thermal Stress Analysis of BOF Accounting for Friction Force Experimentally Measured Yuta HINO, 1) Yoshisato KIYOTA 1) and Yoshimitsu HATTORI 2) 1) Slag and Refractories Department, Steel
More informationCoupled Thermomechanical Contact Problems
Coupled Thermomechanical Contact Problems Computational Modeling of Solidification Processes C. Agelet de Saracibar, M. Chiumenti, M. Cervera ETS Ingenieros de Caminos, Canales y Puertos, Barcelona, UPC
More informationComputational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground
IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Computational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground To cite this article: Jozef Vlek and Veronika
More informationDeflections and Strains in Cracked Shafts due to Rotating Loads: A Numerical and Experimental Analysis
Rotating Machinery, 10(4): 283 291, 2004 Copyright c Taylor & Francis Inc. ISSN: 1023-621X print / 1542-3034 online DOI: 10.1080/10236210490447728 Deflections and Strains in Cracked Shafts due to Rotating
More informationA NEW SIMPLIFIED AND EFFICIENT TECHNIQUE FOR FRACTURE BEHAVIOR ANALYSIS OF CONCRETE STRUCTURES
Fracture Mechanics of Concrete Structures Proceedings FRAMCOS-3 AEDFCATO Publishers, D-79104 Freiburg, Germany A NEW SMPLFED AND EFFCENT TECHNQUE FOR FRACTURE BEHAVOR ANALYSS OF CONCRETE STRUCTURES K.
More informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationME 383S Bryant February 17, 2006 CONTACT. Mechanical interaction of bodies via surfaces
ME 383S Bryant February 17, 2006 CONTACT 1 Mechanical interaction of bodies via surfaces Surfaces must touch Forces press bodies together Size (area) of contact dependent on forces, materials, geometry,
More informationThe Dynamic Stress State of the Wheel-Rail Contact
Proceedings of the 2nd IASME / WSEAS International Conference on Continuum Mechanics (CM'07), Portoroz, Slovenia, May 15-17, 2007 127 The Dynamic Stress State of the Wheel-Rail Contact XIN ZHAO *, ZILI
More informationCO-ROTATIONAL DYNAMIC FORMULATION FOR 2D BEAMS
COMPDYN 011 ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering M. Papadrakakis, M. Fragiadakis, V. Plevris (eds.) Corfu, Greece, 5-8 May 011 CO-ROTATIONAL
More informationSoftware Verification
EXAMPLE 1-026 FRAME MOMENT AND SHEAR HINGES EXAMPLE DESCRIPTION This example uses a horizontal cantilever beam to test the moment and shear hinges in a static nonlinear analysis. The cantilever beam has
More informationAnalysis of thin plate structures using the absolute nodal coordinate formulation
345 Analysis of thin plate structures using the absolute nodal coordinate formulation K Dufva 1 and A A Shabana 2 1 Department of Mechanical Engineering, Lappeenranta University of echnology, Lappeenranta,
More informationStresses and Strains in flexible Pavements
Stresses and Strains in flexible Pavements Multi Layered Elastic System Assumptions in Multi Layered Elastic Systems The material properties of each layer are homogeneous property at point A i is the same
More informationComparison of Models for Finite Plasticity
Comparison of Models for Finite Plasticity A numerical study Patrizio Neff and Christian Wieners California Institute of Technology (Universität Darmstadt) Universität Augsburg (Universität Heidelberg)
More informationEQUATIONS OF MOTION: GENERAL PLANE MOTION (Section 17.5) Today s Objectives: Students will be able to analyze the planar kinetics of a rigid body
EQUATIONS OF MOTION: GENERAL PLANE MOTION (Section 17.5) Today s Objectives: Students will be able to analyze the planar kinetics of a rigid body undergoing general plane motion. APPLICATIONS As the soil
More informationStructural Analysis of Truss Structures using Stiffness Matrix. Dr. Nasrellah Hassan Ahmed
Structural Analysis of Truss Structures using Stiffness Matrix Dr. Nasrellah Hassan Ahmed FUNDAMENTAL RELATIONSHIPS FOR STRUCTURAL ANALYSIS In general, there are three types of relationships: Equilibrium
More informationCHAPTER 8: ROTATIONAL OF RIGID BODY PHYSICS. 1. Define Torque
7 1. Define Torque 2. State the conditions for equilibrium of rigid body (Hint: 2 conditions) 3. Define angular displacement 4. Define average angular velocity 5. Define instantaneous angular velocity
More informationEffect of Angular movement of Lifting Arm on Natural Frequency of Container Lifting Mechanism using Finite Element Modal Analysis
Effect of Angular movement of Lifting Arm on Natural Frequency of Container Lifting Mechanism using Finite Element Modal Analysis Khodu M Dhameliya, 2 Utpal V Shah, 3 Dhaval Makwana, 4 Mansi Yadav, 5 Ishankumar
More informationMechanics II. Which of the following relations among the forces W, k, N, and F must be true?
Mechanics II 1. By applying a force F on a block, a person pulls a block along a rough surface at constant velocity v (see Figure below; directions, but not necessarily magnitudes, are indicated). Which
More informationFREE VIBRATION ANALYSIS OF THIN CYLINDRICAL SHELLS SUBJECTED TO INTERNAL PRESSURE AND FINITE ELEMENT ANALYSIS
FREE VIBRATION ANALYSIS OF THIN CYLINDRICAL SHELLS SUBJECTED TO INTERNAL PRESSURE AND FINITE ELEMENT ANALYSIS J. Kandasamy 1, M. Madhavi 2, N. Haritha 3 1 Corresponding author Department of Mechanical
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 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 informationFigure 2-1: Stresses under axisymmetric circular loading
. Stresses in Pavements.1. Stresses in Fleible Pavements.1.1. Stresses in Homogeneous Mass Boussinesq formulated models for the stresses inside an elastic half-space due to a concentrated load applied
More information[7] Torsion. [7.1] Torsion. [7.2] Statically Indeterminate Torsion. [7] Torsion Page 1 of 21
[7] Torsion Page 1 of 21 [7] Torsion [7.1] Torsion [7.2] Statically Indeterminate Torsion [7] Torsion Page 2 of 21 [7.1] Torsion SHEAR STRAIN DUE TO TORSION 1) A shaft with a circular cross section is
More informationLecture 6 Friction. Friction Phenomena Types of Friction
Lecture 6 Friction Tangential forces generated between contacting surfaces are called friction forces and occur to some degree in the interaction between all real surfaces. whenever a tendency exists for
More information3D ANALYSIS OF H-M COUPLED PROBLEM WITH ZERO-THICKNESS INTERFACE ELEMENTS APPLIED TO GEOMECHANICS
Environmental 3D analysis of H-M and Geosciences coupled problem with zero-thickness interface elements applied to Geomechanics XIII International Conference on Computational Plasticity. Fundamentals and
More informationMODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING EPS SEMINAR SERIES MARCH 2008
MODELING GEOMATERIALS ACROSS SCALES JOSÉ E. ANDRADE DEPARTMENT OF CIVIL AND ENVIRONMENTAL ENGINEERING EPS SEMINAR SERIES MARCH 2008 COLLABORATORS: DR XUXIN TU AND MR KIRK ELLISON THE ROADMAP MOTIVATION
More informationDynamic and buckling analysis of FRP portal frames using a locking-free finite element
Fourth International Conference on FRP Composites in Civil Engineering (CICE8) 22-24July 8, Zurich, Switzerland Dynamic and buckling analysis of FRP portal frames using a locking-free finite element F.
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 informationRatcheting and Rolling Contact Fatigue Crack Initiation Life of Rails under Service Loading. Wenyi YAN Monash University, Australia
Ratcheting and Rolling Contact Fatigue Crack Initiation Life of Rails under Service Loading Wenyi YAN Monash University, Australia Chung Lun PUN Peter Mutton Qianhua Kan Guozheng Kang Contents Introduction
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 information