Proceedings oh the 18th IAHR International Symposium on Ice (2006) DISCRETE ELEMENT SIMULATION OF ICE PILEUP AGAINST AN INCLINED STRUCTURE


 Job Garrison
 10 months ago
 Views:
Transcription
1 DISCRETE ELEMENT SIMULATION OF ICE PILEUP AGAINST AN INCLINED STRUCTURE Jani Paavilainen, Jukka Tuhkuri and Arttu Polojärvi Helsinki University of Technology, Laboratory for Mechanics of Materials, Finland ABSTRACT A two dimensional discrete element method has been used to study the ice pileup process against an inclined plate. In a simulation, an intact ice sheet, driven against an inclined plate, failed into discrete ice blocks which then accumulated into a pile. The pileup process, forces acting on the plate and the shape of the rubble pile were studied. Simulation results were compared to laboratory experiments. KEY WORDS: pileup, inclined structure, simulation INTRODUCTION The failure of an ice cover against an inclined structure is a process where an ice sheet breaks into discrete ice blocks which then accumulate in the icestructure interface and form a rubble pile. The first models for calculating the ice load on sloping structures separated the failure process into failure of the ice sheet through bending, and rideup of the ice along the structure (for a review, see Sanderson, 1988). The main variables in this kind of models are the slope angle, thickness and other properties of the ice sheet, as well as the friction coefficient. Later, also other failure modes and effects of weight and buoyancy of the rubble have been included in the calculation models (e.g. Määttänen, 1986). While solutions for different events can be obtained through simple analyses, calculation of the total ice load is not as easy. This is because the maximum value of each load component does not necessarily occur at the same time, and also because the different ice failure events can interact. As an example, the bending failure of the ice sheet is affected by the rubble formation. This kind of uncertainties can be avoided if the icestructure interaction is studied as a process and simulated. One method to simulate the ice failure processes is the discrete element method (DEM), which has already been used in analyses of ice pileup on an inclined ramp (Hopkins, 1997) and ice sheet interaction with a rubble pile adjacent to a structure (Evgin et al., 1992). In addition, it is important to note, that the pileup process is closely related to the ice ridging process, which has been extensively simulated with DEM (e.g. Hopkins, 1998; Hopkins et al. 1999)
2 In this paper, a two dimensional DEM is used to simulate laboratory scale experiments of ice sheet failure and rubble formation against an inclined structure. The simulations provide a method to study the rubble formation process and the forces related to it, while the experimental data is used to benchmark the simulation method. DISCRETE ELEMENT METHOD Discrete element method (DEM) is a technique for modelling media consisting of discrete particles. In DEM, the interaction between the particles within the system occurs through pairwise collisions from which the forces acting on each particle are derived according to the chosen interaction laws. On each time step of a DEM simulation, the contacting particles are identified, the contact geometries are solved, the contact forces are determined, and the particles are moved to new locations. The method for searching contacting particles was adopted from Munjiza and Andrews (1998) and the method for calculating contact forces between interacting ice blocks was adopted from Hopkins (1992). The contact forces between particles were calculated with an elasticviscousplastic material model combined with an incremental MohrCoulomb tangential force model. Inelasticity was modelled by using a plastic limit for the material. Determination of the elastic contact forces was based on the overlap area and the rate of change of the overlap area of contacting particles. The normal component of the elastic contact force is F ne = k ne (A e + A)n, (1) where k ne is the normal stiffness of the contact, A e is the elastic area of the overlap from the previous time step, A is the change in the overlap area between time steps and n is the unit vector normal to the line of contact defined for the contact geometry. The viscous component of the normal force is given by F nv = k nv ( A t ) n, where k nv is the viscous damping coefficient and t is the length of time step. The plastic limit is defined as (2) F p = σ p l n, (3) where σ p is the plastic limit of the material and l is the length of the line of contact. Using F p as the upper limit for the normal contact force, the total normal contact force is given by F n = min {F ne + F nv,f p }. (4) If the plastic limit is exceeded, the elastic area for the next time step is updated by solving the equation A e t = 1 k nv ( F p k ne A e ). (5) The tangential contact force F t,i on timestep i is calculated using its value from time step i 1, tangential contact stiffness k te and relative tangential velocity of interacting particles v t. The 178
3 upper limit for F t,i, using MohrCoulomb model, is achieved by using friction coefficient µ and F n, hence F t,i = min {F t,i 1 k te v t t, µf n }. (6) The moment acting on each particle due to contacts is given by M = r F, (7) where r is the vector pointing from the centroid of a particle to the centroid of the overlap area and F is achieved as a sum of the tangential and normal contact forces. When simulating an ice sheet failure against a structure, an intact ice sheet and the failure of the ice sheet need to be modelled, in addition to the discrete particles. This intact 2D sheet was modelled with discrete particles that were combined together with beam elements. For the beam elements, a nonlinear Timoshenko plane beam model was used (Paavilainen, 2006). From the positions and orientations of two adjacent particles in the intact sheet, the strainsh T = [ε γ κ] of the beam element combining the particles are calculated. The stress resultants z T = [N V M] are then defined from the strains, and the internal forces p of the beam element are calculated from the stress resultants as follows p = B T z dx. (8) L 0 The matrix B is obtained from the straindisplacement relation and the integral is taken over the initial length L 0 of the beam element. Internal forces of the beam element are then added to the corresponding forces acting on the particle. A fracture criterion, σ x σ b, where σ x is the stress in the surface of the beam and σ b is the bending strength of the beam, was used for the ice sheet. On each time step, the fracture criterion was examined on both the upper and the lower surface of the beam. If σ x σ b, a crack was assumed to initiate and propagate with constant velocity through the beam during the following time steps. The crack velocity in many brittle materials is of the order of 0.3 E/ρ, where E is Young s modulus and ρ is density (Broek, 1991). Parsons et al. (1987) measured a much lower value 0.01 E/ρ for sea ice. Crack velocity in the model ice simulated here is not known, and a value 0.05 E/ρ was used in the simulations. SETUP OF THE SIMULATIONS AND THE LABORATORY EXPERIMENTS The numerical DE simulations performed were simulations of the model scale laboratory experiments conducted by Saarinen (2000). In both experiments and simulations, an ice sheet was pushed against an inclined plate with a constant velocity. Figure 1 shows the geometry of the setup and Table 1 gives the parameters used in the simulations and either measured during the laboratory experiments or known for the model scale laboratory ice. In the experiments, a 4.6 m wide ice sheet was used, but the simulations were two dimensional. All the individual particles in the simulation were square shaped. At the left boundary, the ice sheet rotation was constrained, vertical displacement was allowed, and a constant horizontal velocity was given as an input value. Distance from the inclined plate to the left boundary was 217.5h, where h is the ice thickness. The number of ice particles in the beginning of the simulation was 217 and at the end Both ends of the inclined plate were rigidly supported
4 Table 1: Parameters used in the simulations and either measured during the laboratory experiments or known for the model scale laboratory ice. Parameter Simulations Experiments Ice particles thickness [m], h width [m] Ice properties Young s modulus [MPa], E i Poisson ratio, ν i density [kg m 3 ], ρ i friction coefficient, µ i friction coefficient (under water), µ is normal stiffness [MPa], k ne tangential stiffness [MPa], k te viscous damping [MPa s], k nv plastic limit [kpa], σ p bending strength [kpa], σ b internal damping [kpa s], c i,int crack propagation speed [m s 1 ] ice sheet velocity, v i 0.1, Water density [kg m 3 ], ρ w Wall properties inclination angle [ ], α 50 50, 65, 80 Young s modulus [GPa], E W Poisson ratio, ν W internal damping [kpa s], c W,int Icewall contact friction coefficient, µ iw , 0.35 friction coefficient (under water), µ iws Simulation time step [s], t SIMULATION RESULTS Figure 1 shows snapshots of the simulation and Figure 2 gives the total force acting on the inclined plate in the horizontal direction as a function of the length of pushed ice L. Figure 2 shows also the force obtained in the laboratory experiments. For comparison of the forces and the stage of simulation, the subfigures of Figure 1 are pointed out in the F(L)record in Figure 2 with arrows. The growth of the cross sectional area of the rubble pile was also studied. Cross sections from the simulation are presented in Figure 3 at L = 10, 20, 30 and 40 m. Initially, the ice sheet failed against the inclined plate by bending (Figure 1a) and started to form a pile (Figure 1b) which grew both vertically and horizontally. After a pile was formed, the ice sheet failed against the pile to form a sail (Figures 1d and 1g) or slided on top of the pile to the inclined plate and rodeup along the plate (Figures 1c and 1e). Figures 1e and 1f illustrate a collapse of the pile, when a failure of the ice sheet resulted in decrease of the horizontal force, followed by movement of the pile down and left. The cycle, growth of the pile against the plate and collapse of the pile, occurred several times during the simulation. During the initial stages (0 L < 20 m) of the simulation these collapses of the rubble pile caused the force to drop to zero as the contact between the sheet 180
5 (a) (b) (c) (d) (e) (f) (g) (h) Figure 1. Snapshots of a simulation. Ice sheet (white) is pushed with a constant velocity of 0.05 m/s against an inclined plate (black). The water layer is coloured gray. The length of the pushed ice L in subfigures (a)(h) can be found from Figure 2. and the plate was lost. This is evident from the F(L)record. Another phenomenon observed in the simulations was a clockwise rotational motion of the rubble pile. In later stages, the location of the active failure prosess moved away from the plate (Figures 1g and 1h). Two velocities, 0.1 m/s and 0.05 m/s, were used in the simulations. In this range, the velocity did not affect the calculated F(L)record or the final shape of the rubble pile profile
6 simulation experiment (e) F [N/m] (a) (c) (g) 0 (f) (h) (b) (d) L [m] Figure 2: Horizontal force F acting on the plate as a function of the length of the pushed ice L as obtained from the simulation and experiments. The letters refer to Figure 1. DISCUSSION AND CONCLUSIONS Visual observations of the simulations and laboratory experiments were similar. In both cases cycles of rubble pile growth and collapse were observed. More quantitative data can be obtained through comparison of the rubble pile profiles and measured forces. Figure 3 shows rubble pile profiles during the simulation and Figure 4 compares the rubble profiles from the simulations and experiments. The shapes of the rubble piles from the experiments are similar with those achieved from the simulations. However, it appears that in the simulations the rubble piles were more loosely packed. The porosity ν p of the piles was in the range of in the experiments and in the simulations. Here porosity ν p is defined as (1 ν p )A = hl, where A is the cross sectional area of the pile as measured from 0.5 L=30 m L =40 m L =20 m L =10 m Figure 3: Rubble pile profiles from the simulation with different lengths of the pushed ice L
7 (a) (b) Figure 4. Comparisons of rubble pile profiles from laboratory experiments and simulations. (a) L = 14.6 m with an angle of inclination of 65 in experiments, and (b) L = 16.6 m with an angle of inclination of 50 in experiments. Data from simulations is presented with bold lines and experimental data with dashed lines. the cross sectional profiles. One reason for this difference in porosity may be edge crushing. During the experiments, edge crushing leads to formation of very small ice particles which may fill the pores in the rubble. In the simulations, the ice blocks failed only through bending, although an upper limit for the normal contact force is defined by Equation 4. The F(L)records in Figure 2 show that the average force obtained from the simulation and experiments reached the same level of about 400 N/m. However, during the initial stages the F(L)records look rather different. This is due to the loss of contact between the pile and the inclined plate during initial stages of the simulation. When the contact is lost, the simulated force drops to zero. This is caused by the 2D nature of the simulations and is thus not observed during the experiments. In the later stages of the simulation, when the contact between ice and the structure was not lost, the simulated F(L)record is similar with the experimental data. The ice force obtained from the simulations can further be compared with results from analytical methods. A textbook solution to the horizontal ice force F h on a slopingsided structure (Sanderson, 1988) divides the force into an ice breaking force F b and an ice rideup F r force ( ) ρw gh 5 1/4 ( F h = F b + F r = σ b 0.68 sin α + µ ) iw cos α E cosα µ iw sin α ( sin α + µiw cosα + zhρ i g (sin α + µ iw cosα) cos α µ iw sin α + cosα ) sin α, (9) where z is the height of rideup and the other symbols are as explained in Table 1. At the simulation event shown in Figure 1e, z = 0.75 m and the simulated force was about 1400 N/m. By using the values given in Table 1 and z = 0.75 m, a force F h = 1353 N/m is obtained by Equation
8 ACKNOWLEDGEMENTS The financial support from the Academy of Finland, the Finnish Cultural Foundation, and Heikki and Hilma Honkanen Foundation is gratefully acknowledged. REFERENCES Broek, D. (1991), "Elementary engineering fracture mechanics", Kluwer, 516 p. Evgin, E., Zhan, C. and Timco, G.W. (1992), "Distinct element modelling of load transmission through grounded ice rubble," Proc. of the 11th OMAE Symp., Calgary, Canada, Vol IV, pp Hopkins, M.A. (1992), "Numerical simulation of systems of multitudinous polygonal blocks", Cold Regions Research and Engineering Laboratory, CRREL, Report 9222, 69 p. Hopkins, M.A. (1997), "Onshore ice pileup: a comparison between experiments and simulations," Cold Regions Science and Technology, 26: Hopkins, M.A. (1998), "Four stages of pressure ridging," Journal of Geophysical Research, 103, pp. 21,88321,891. Hopkins, M.A., Tuhkuri, J. and Lensu, M. (1999), "Rafting and ridging of thin ice sheets," Journal of Geophysical Research, 104, pp Munjiza, A., Andrews K. F. R. (1998), "NBS Contact detection algorithm for bodies of similar size", Int. J. Numer. Meth. Engng., 43, pp Määttänen, M. (1986), "Ice sheet failure against an inclined wall", Proc. of the IAHR Ice Symposium 1986, Iowa City, USA, Vol I, pp Paavilainen, J. (2006), "Modelling of beams in discrete element simulations", To appear in the Proceedings of IX Finnish Mechanics Days. Parsons, B.L., Snellen, J.B. and Hill, B., (1987), "Preliminary measurements of terminal crack velocity in ice", Cold Regions Science and Technology, 13: Saarinen, S. (2000), "Description of the pileup process of an ice sheet against an inclined plate", M.Sc. thesis, Helsinki University of Technology, Dept. of Mechanical Engineering, 78 p. Sanderson, T.J.O. (1988). "Ice Mechanics, Risks to Offshore Structures". Graham & Trotman. 253 p
Brittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm
Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd
More informationFoundation Analysis LATERAL EARTH PRESSURE
Foundation Analysis LATERAL EARTH PRESSURE INTRODUCTION Vertical or nearvertical slopes of soil are supported by retaining walls, cantilever sheetpile walls, sheetpile bulkheads, braced cuts, and other
More informationDynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models
Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models I. Rhee, K.J. Willam, B.P. Shing, University of Colorado at Boulder ABSTRACT: This paper examines the global
More informationModelling Progressive Failure with MPM
Modelling Progressive Failure with MPM A. Yerro, E. Alonso & N. Pinyol Department of Geotechnical Engineering and Geosciences, UPC, Barcelona, Spain ABSTRACT: In this work, the progressive failure phenomenon
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@ujfgrenoble.fr http://wwwlgit.obs.ujfgrenoble.fr/users/lbaillet/
More informationFCP Short Course. Ductile and Brittle Fracture. Stephen D. Downing. Mechanical Science and Engineering
FCP Short Course Ductile and Brittle Fracture Stephen D. Downing Mechanical Science and Engineering 001015 University of Illinois Board of Trustees, All Rights Reserved Agenda Limit theorems Plane Stress
More informationWellbore stability analysis in porous carbonate rocks using cap models
Wellbore stability analysis in porous carbonate rocks using cap models L. C. Coelho 1, A. C. Soares 2, N. F. F. Ebecken 1, J. L. D. Alves 1 & L. Landau 1 1 COPPE/Federal University of Rio de Janeiro, Brazil
More informationNONLINEAR VISCOELASTIC MODEL OF STRUCTURAL POUNDING
3 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 6, 004 Paper No. 308 NONLINEAR VISCOELASTIC MODEL OF STRUCTURAL POUNDING Robert JANKOWSKI SUMMARY Pounding between structures
More informationLecture 7, Foams, 3.054
Lecture 7, Foams, 3.054 Opencell foams StressStrain curve: deformation and failure mechanisms Compression  3 regimes  linear elastic  bending  stress plateau  cell collapse by buckling yielding
More informationExample3. Title. Description. Cylindrical Hole in an Infinite MohrCoulomb Medium
Example3 Title Cylindrical Hole in an Infinite MohrCoulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elastoplastic
More informationA NEW SIMPLIFIED AND EFFICIENT TECHNIQUE FOR FRACTURE BEHAVIOR ANALYSIS OF CONCRETE STRUCTURES
Fracture Mechanics of Concrete Structures Proceedings FRAMCOS3 AEDFCATO Publishers, D79104 Freiburg, Germany A NEW SMPLFED AND EFFCENT TECHNQUE FOR FRACTURE BEHAVOR ANALYSS OF CONCRETE STRUCTURES K.
More informationThe Frictional Regime
The Frictional Regime Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/25/2016 10:08 AM We Discuss The Frictional Regime Processes of Brittle Deformation
More informationLateral Earth Pressure
1 of 11 6/2/2012 4:28 AM Lateral Earth Pressure The magnitude of lateral earth pressure depends on: 1. Shear strength characteristics of soil 2. Lateral strain condition 3. Pore water pressure 4. State
More informationfriction friction ab slow fast increases during sliding
µ increases during sliding faster sliding > stronger fault > slows sliding leads to stable slip: no earthquakes can start velocitystrengthening friction slow fast µ velocitystrengthening friction
More informationRecent Research on EPS Geofoam Seismic Buffers. Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen srmc Canada
Recent Research on EPS Geofoam Seismic Buffers Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen srmc Canada What is a wall (SEISMIC) buffer? A compressible inclusion placed between
More informationSHEAR STRENGTH I YULVI ZAIKA
SHEAR STRENGTH I YULVI ZAIKA MATERI Keruntuhan mohr coulomb, stress paths, kuat geser tanah non kohesif dan kohesif, evaluasi kuat geser di lapangan, tegangan normal dan tegangan geser pada sebuah bidang
More informationINELASTIC RESPONSES OF LONG BRIDGES TO ASYNCHRONOUS SEISMIC INPUTS
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 16, 24 Paper No. 638 INELASTIC RESPONSES OF LONG BRIDGES TO ASYNCHRONOUS SEISMIC INPUTS Jiachen WANG 1, Athol CARR 1, Nigel
More informationAnalysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method
Analysis of asymmetric radial deformation in pipe with local wall thinning under internal pressure using strain energy method V.M.F. Nascimento Departameto de ngenharia Mecânica TM, UFF, Rio de Janeiro
More informationPILESUPPORTED RAFT FOUNDATION SYSTEM
PILESUPPORTED RAFT FOUNDATION SYSTEM Emre Biringen, Bechtel Power Corporation, Frederick, Maryland, USA Mohab Sabry, Bechtel Power Corporation, Frederick, Maryland, USA Over the past decades, there has
More information3D ANALYSIS OF STRESSES AROUND AN UNLINED TUNNEL IN ROCK SUBJECTED TO HIGH HORIZONTAL STRESSES
3D ANALYSIS OF STRESSES AROUND AN UNLINED TUNNEL IN ROCK SUBJECTED TO HIGH HORIZONTAL STRESSES Abdel Meguid, M. Graduate Student, Department of Civil Engineering, University of Western Ontario, London,
More information3.032 Problem Set 2 Solutions Fall 2007 Due: Start of Lecture,
3.032 Problem Set 2 Solutions Fall 2007 Due: Start of Lecture, 09.21.07 1. In the beam considered in PS1, steel beams carried the distributed weight of the rooms above. To reduce stress on the beam, it
More informationAim of the study Experimental determination of mechanical parameters Local buckling (wrinkling) Failure maps Optimization of sandwich panels
METNET Workshop October 1112, 2009, Poznań, Poland Experimental and numerical analysis of sandwich metal panels Zbigniew Pozorski, Monika ChudaKowalska, Robert Studziński, Andrzej Garstecki Poznan University
More informationProceedings of the ASME th International Conference on Ocean, Offshore and Arctic Engineering OMAE2016 June 1924, 2016, Busan, South Korea
Proceedings of the ASME 26 35th International Conference on Ocean, Offshore and Arctic Engineering OMAE26 June 924, 26, Busan, South Korea OMAE2654554 LOCAL STRAIN AND STRESS CALCULATION METHODS OF IRREGULAR
More informationSand Control Rock Failure
Sand Control Rock Failure Why? A bit of Mechanics on rock failure How? Some choices that depend on the rock What is moving? Sand grains? Fines? 3/14/2009 1 Young s Modulus, E Young s Modulus is a material
More informationMechanics of Solids. Mechanics Of Solids. Suraj kr. Ray Department of Civil Engineering
Mechanics Of Solids Suraj kr. Ray (surajjj2445@gmail.com) Department of Civil Engineering 1 Mechanics of Solids is a branch of applied mechanics that deals with the behaviour of solid bodies subjected
More informationDynamics Manual. Version 7
Dynamics Manual Version 7 DYNAMICS MANUAL TABLE OF CONTENTS 1 Introduction...11 1.1 About this manual...11 2 Tutorial...21 2.1 Dynamic analysis of a generator on an elastic foundation...21 2.1.1 Input...21
More informationReservoir Geomechanics with ABAQUS
Reservoir Geomechanics with ABAQUS B. Bostrøm and E. Skomedal Statoil ASA, Norway Abstract: The coupled hydromechanical behavior of two North Sea highpressure/hightemperature gascondensate fields during
More informationProceedings, International Snow Science Workshop, Banff, 2014
SNOW SLAB FAILURE DUE TO BIOTTYPE ACOUSTIC WAVE PROPAGATION Rolf Sidler 2 Simon Fraser University, 8888 University Drive, Burnaby V5A1S6, Canada ABSTRACT: Even though seismic methods are among the most
More informationBACKGROUNDS. Two Models of Deformable Body. Distinct Element Method (DEM)
BACKGROUNDS Two Models of Deformable Body continuum rigidbody spring deformation expressed in terms of field variables assembly of rigidbodies connected by spring Distinct Element Method (DEM) simple
More informationSoil strength. the strength depends on the applied stress. water pressures are required
Soil Strength Soil strength u Soils are essentially frictional materials the strength depends on the applied stress u Strength is controlled by effective stresses water pressures are required u Soil strength
More informationEXTENDED ABSTRACT. Combined Pile Raft Foundation
EXTENDED ABSTRACT Combined Pile Raft Foundation Rui Diogo Gomes da Silva Supervisor: Prof. Jaime Alberto dos Santos December 2009 1. Introduction The piled raft foundation is an innovative design concept
More informationAN IMPORTANT PITFALL OF PSEUDOSTATIC FINITE ELEMENT ANALYSIS
AN IMPORTANT PITFALL OF PSEUDOSTATIC FINITE ELEMENT ANALYSIS S. Kontoe, L. Pelecanos & D.M. Potts ABSTRACT: Finite Element (FE) pseudostatic analysis can provide a good compromise between simplified
More informationMATERIALS. Why do things break? Why are some materials stronger than others? Why is steel tough? Why is glass brittle?
MATERIALS Why do things break? Why are some materials stronger than others? Why is steel tough? Why is glass brittle? What is toughness? strength? brittleness? Elemental material atoms: A. Composition
More informationChapter 7. Highlights:
Chapter 7 Highlights: 1. Understand the basic concepts of engineering stress and strain, yield strength, tensile strength, Young's(elastic) modulus, ductility, toughness, resilience, true stress and true
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 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 informationFinite Element Modeling of an Aluminum Tricycle Frame
Finite Element Modeling of an Aluminum Tricycle Frame A. Rodríguez, B. Chiné*, and J. A. Ramírez Costa Rica Institute of Technology, School of Materials Science and Engineering, Cartago, Costa Rica *Corresponding
More informationDynamic Earth Pressure Problems and Retaining Walls. Behavior of Retaining Walls During Earthquakes. Soil Dynamics week # 12
Dynamic Earth Pressure Problems and Retaining Walls 1/15 Behavior of Retaining Walls During Earthquakes  Permanent displacement = cc ' 2 2 due to one cycle of ground motion 2/15 Hence, questions are :
More informationLoads Exerted on a Cylindrical Structure by Floating Ice Modelled as a ViscousPlastic Material
Archives of HydroEngineering and Environmental Mechanics Vol. 53 (2006), No. 2, pp. 105 126 IBW PAN, ISSN 1231 3726 Loads Exerted on a Cylindrical Structure by Floating Ice Modelled as a ViscousPlastic
More informationTMHL TMHL (Del I, teori; 1 p.) SOLUTION I. II.. III. Fig. 1.1
TMHL6 2040830 (Del I, teori; p.). Fig.. shows three cases of sharp cracks in a sheet of metal. In all three cases, the sheet is assumed to be very large in comparison with the crack. Note the different
More informationA CLOSER LOOK AT THE BRAZILIAN TEST AND ITS MODE OF FAILURE Arvid Landva, GEMTEC Limited, Fredericton, New Brunswick, Canada
A CLOSER LOOK AT THE BRAZILIAN TEST AN ITS MOE OF FAILURE Arvid Landva, GEMTEC Limited, Fredericton, New Brunswick, Canada ABSTRACT Timoshenko (1934) showed that a compressive line load applied perpendicularly
More informationCONSOLIDATION BEHAVIOR OF PILES UNDER PURE LATERAL LOADINGS
VOL., NO., DECEMBER 8 ISSN 898 8 Asian Research Publishing Network (ARPN). All rights reserved. CONSOLIDATION BEAVIOR OF PILES UNDER PURE LATERAL LOADINGS Qassun S. Mohammed Shafiqu Department of Civil
More informationBasic Energy Principles in Stiffness Analysis
Basic Energy Principles in Stiffness Analysis StressStrain Relations The application of any theory requires knowledge of the physical properties of the material(s) comprising the structure. We are limiting
More informationIntroduction to Engineering Materials ENGR2000. Dr. Coates
Introduction to Engineering Materials ENGR2 Chapter 6: Mechanical Properties of Metals Dr. Coates 6.2 Concepts of Stress and Strain tension compression shear torsion Tension Tests The specimen is deformed
More informationExperimental and theoretical characterization of Li 2 TiO 3 and Li 4 SiO 4 pebbles
Experimental and theoretical characterization of Li 2 TiO 3 and Li 4 SiO 4 s D. Aquaro 1 N. Zaccari ABSTRACT Dipartimento di Ingegneria Meccanica Nucleare e della Produzione University of Pisa (Italy)
More informationHOW WELL CAN WE PREDICT ICE LOADS?
HOW WELL CAN WE PREDICT ICE LOADS? G.W. Timco 1 and K.R. Croasdale 2 1 Canadian Hydraulics Centre, National Research Council, Ottawa, ON, K1A R6 Canada 2 K.R. Croasdale & Assoc., Calgary, AL, T2R 1M5 Canada
More informationMechanical Design in Optical Engineering
OPTI Buckling Buckling and Stability: As we learned in the previous lectures, structures may fail in a variety of ways, depending on the materials, load and support conditions. We had two primary concerns:
More informationA Sample Durability Study of a Circuit Board under Random Vibration and Design Optimization
A Sample Durability Study of a Circuit Board under Random Vibration and Design Optimization By: MS.ME Ahmad A. Abbas Ahmad.Abbas@AdvancedCAE.com www.advancedcae.com Sunday, March 07, 2010 Advanced CAE
More informationVIBRATION ENERGY FLOW IN WELDED CONNECTION OF PLATES. 1. Introduction
ARCHIVES OF ACOUSTICS 31, 4 (Supplement), 53 58 (2006) VIBRATION ENERGY FLOW IN WELDED CONNECTION OF PLATES J. CIEŚLIK, W. BOCHNIAK AGH University of Science and Technology Department of Robotics and Mechatronics
More informationBE Semester I ( ) Question Bank (MECHANICS OF SOLIDS)
BE Semester I ( ) Question Bank (MECHANICS OF SOLIDS) All questions carry equal marks(10 marks) Q.1 (a) Write the SI units of following quantities and also mention whether it is scalar or vector: (i)
More informationPROGRESSIVE FAILURE MECHANISMS IN JOINTED ROCK: INSIGHT FROM 3D DEM MODELLING
II International Conference on Particlebased Methods Fundamentals and Applications PARTICLES 2011 E. Oñate and D.R.J. Owen (Eds) PROGRESSIVE FAILURE MECHANISMS IN JOINTED ROCK: INSIGHT FROM 3D DEM MODELLING
More informationSIMPLIFIED MODELING OF THINWALLED TUBES WITH OCTAGONAL CROSS SECTION AXIAL CRUSHING. Authors and Correspondance: Abstract:
SIMPLIFIED MODELING OF THINWALLED TUBES WITH OCTAGONAL CROSS SECTION AXIAL CRUSHING Authors and Correspondance: Yucheng Liu, Michael L. Day Department of Mechanical Engineering University of Louisville
More informationPrinciples Of Engineering. Part A
Principles Of Engineering Final Examination Part A Fall 2007 Student Name: Date: Class Period: Total Points: /40 Converted Score: /50 Page 1 of 11 Directions: Circle the letter of the response that best
More informationMECE 3321 MECHANICS OF SOLIDS CHAPTER 3
MECE 3321 MECHANICS OF SOLIDS CHAPTER 3 Samantha Ramirez TENSION AND COMPRESSION TESTS Tension and compression tests are used primarily to determine the relationship between σ avg and ε avg in any material.
More informationUNIT I SIMPLE STRESSES AND STRAINS
Subject with Code : SM1(15A01303) Year & Sem: IIB.Tech & ISem SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) UNIT I SIMPLE STRESSES
More informationModelling installation of helical anchors in clay
Modelling installation of helical anchors in clay C. Todeshkejoei & J.P. Hambleton ARC Centre of Excellence for Geotechnical Science and Engineering, Centre for Geotechnical and Materials Modelling, The
More information1.050 Content overview Engineering Mechanics I Content overview. Selection of boundary conditions: Euler buckling.
.050 Content overview.050 Engineering Mechanics I Lecture 34 How things fail and how to avoid it Additional notes energy approach I. Dimensional analysis. On monsters, mice and mushrooms Lectures 3. Similarity
More informationTuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE
1 Chapter 3 Load and Stress Analysis 2 Chapter Outline Equilibrium & FreeBody Diagrams Shear Force and Bending Moments in Beams Singularity Functions Stress Cartesian Stress Components Mohr s Circle for
More informationMODEL TESTING OF RIDGE KEEL LOADS ON STRUCTURES PART IV: PRELIMINARY RESULTS OF FREEZE BOND SHEAR STRENGTH EXPERIMENTS
POAC 09 Luleå, Sweden Proceedings of the 20th International Conference on Port and Ocean Engineering under Arctic Conditions June 912, 2009 Luleå, Sweden POAC0968 MODEL TESTING OF RIDGE KEEL LOADS ON
More informationChapter kn m/kg Ans kn m/kg Ans. 187 kn m/kg Ans.
Chapter 1 From Tables A0, A1, A, and A4c, (a) UNS G1000 HR: S ut = 80 (55) MPa (kpsi), S yt = 10 (0) MPa (kpsi) Ans. (b) SAE 1050 CD: S ut = 690 (100) MPa (kpsi), S yt = 580 (84) MPa (kpsi) Ans. (c)
More informationDesign of a fastener based on negative Poisson's ratio foam adapted from
1 Design of a fastener based on negative Poisson's ratio foam adapted from Choi, J. B. and Lakes, R. S., "Design of a fastener based on negative Poisson's ratio foam", Cellular Polymers, 10, 205212 (1991).
More informationChapter 2  Macromechanical Analysis of a Lamina. Exercise Set. 2.1 The number of independent elastic constants in three dimensions are: 2.
Chapter  Macromechanical Analysis of a Lamina Exercise Set. The number of independent elastic constants in three dimensions are: Anisotropic Monoclinic 3 Orthotropic 9 Transversely Orthotropic 5 Isotropic.
More information= 50 ksi. The maximum beam deflection Δ max is not = R B. = 30 kips. Notes for Strength of Materials, ET 200
Notes for Strength of Materials, ET 00 Steel Six Easy Steps Steel beam design is about selecting the lightest steel beam that will support the load without exceeding the bending strength or shear strength
More informationPiles in Lateral Spreading due to Liquefaction: A Physically Simplified Method Versus Centrifuge Experiments
"PileGroup Response to Large Soil Displacements and Liquefaction: Centrifuge Experiments Versus A Physically Simplified Analysis", Journal of Geotechnical and Geoenvironmental Engineering, ASCE, Vol.
More informationStress and fabric in granular material
THEORETICAL & APPLIED MECHANICS LETTERS 3, 22 (23) Stress and fabric in granular material Ching S. Chang,, a) and Yang Liu 2 ) Department of Civil Engineering, University of Massachusetts Amherst, Massachusetts
More informationRock Failure. Topics. Compressive Strength Rock Strength from Logs Polyaxial Strength Criteria Anisotropic Rock Strength Tensile Strength
Rock Failure Topics Compressive Strength Rock Strength from Logs Polyaxial Strength Criteria Anisotropic Rock Strength Tensile Strength Key Points 1. When rock fails in compression, the compressive stress
More informationEstimation of the Passive Earth Pressure with Inclined Cohesive Backfills: the Effect of Intermediate Principal Stress is Considered
108 The Open Mechanical Engineering Journal, 011, 5, 108116 Open Access Estimation of the Passive Earth Pressure with Inclined Cohesive Backfills: the Effect of Intermediate Principal Stress is Considered
More informationLab Practical  Discontinuum Analysis & Distinct Element Method
Lab Practical  Discontinuum Analysis & Distinct Element Method Part A The Basics The Universal Distinct Element Code (UDEC) is a twodimensional numerical program based on the distinct element method
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 informationLecture #10: Anisotropic plasticity Crashworthiness Basics of shell elements
Lecture #10: 1510735: Dynamic behavior of materials and structures Anisotropic plasticity Crashworthiness Basics of shell elements by Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering,
More informationFinite Element Modelling with Plastic Hinges
01/02/2016 Marco Donà Finite Element Modelling with Plastic Hinges 1 Plastic hinge approach A plastic hinge represents a concentrated postyield behaviour in one or more degrees of freedom. Hinges only
More informationEarth Pressure Theory
Lateral Earth Pressure Page 1 Earth Pressure Theory Examples of Retaining Walls Lateral Earth Pressure Page 2 AtRest, Active and Passive Earth Pressure Wednesday, August 17, 2011 12:45 PM Atrest condition
More informationA micromechanical approach to describe internal erosion effects in soils
A micromechanical approach to describe internal erosion effects in soils Luc Scholtès, PierreYves Hicher, Luc Sibille To cite this version: Luc Scholtès, PierreYves Hicher, Luc Sibille. A micromechanical
More informationFracture Mechanics of NonShear Reinforced R/C Beams
Irina Kerelezova Thomas Hansen M. P. Nielsen Fracture Mechanics of NonShear Reinforced R/C Beams DANMARKS TEKNISKE UNIVERSITET Report BYG DTU R154 27 ISSN 1612917 ISBN 9788778772265 Fracture Mechanics
More information(MPa) compute (a) The traction vector acting on an internal material plane with normal n ( e1 e
EN10: Continuum Mechanics Homework : Kinetics Due 1:00 noon Friday February 4th School of Engineering Brown University 1. For the Cauchy stress tensor with components 100 5 50 0 00 (MPa) compute (a) The
More informationChapter (7) Lateral Earth Pressure
Chapter (7) Lateral Earth Pressure Introduction Vertical or near vertical slopes of soil are supported by retaining walls, cantilever sheetpile walls, sheetpile bulkheads, braced cuts, and other similar
More informationAdvanced holographic nondestructive testing system for residual stress analysis
Advanced holographic nondestructive testing system for residual stress analysis Anatoli Kniazkov, Yuri Onischenko, George Dovgalenko, Gregory Salamo Laser Laboratory, Physics Department, University of
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 informationAS PHYSICS TERMS, DEFINITIONS & UNITS
AS PHYSICS TERMS, DEFINITIONS & UNITS May 2015 This document is issued by WJEC Eduqas to assist teachers with the preparation of candidates for the GCE examination in PHYSICS. It consists of the definitions
More informationK.S. Woo 1,2 *, E. Eberhardt 2, D. Elmo 3, D. Stead 4 and P. K. Kaiser 5. Introduction
Benchmark testing of numerical approaches for modelling the influence of undercut depth on caving, fracture initiation and subsidence angles associated with block cave mining K.S. Woo 1,2 *, E. Eberhardt
More information3D Finite Element Modeling of faultslip triggering caused by porepressure
3D Finite Element Modeling of faultslip triggering caused by porepressure changes Arsalan Sattari and David W. Eaton Department of Geoscience, University of Calgary Suary We present a 3D model using a
More informationMining. Slope stability analysis at highway BR153 using numerical models. Mineração. Abstract. 1. Introduction
Mining Mineração http://dx.doi.org/10.1590/037044672015690040 Ricardo Hundelshaussen Rubio Engenheiro Industrial / Doutorando Universidade Federal do Rio Grande do Sul  UFRS Departamento de Engenharia
More informationHONGJUN LI Department of Mechanical Engineering University of Strathclyde Glasgow, Scotland, UK
HONGJUN LI Department of Mechanical Engineering University of Strathclyde Glasgow, Scotland, UK Introduction FEA is an established analysis method used in Pressure Vessel Design Most DBA is still based
More informationProf. B V S Viswanadham, Department of Civil Engineering, IIT Bombay
56 Module 4: Lecture 7 on Stressstrain relationship and Shear strength of soils Contents Stress state, Mohr s circle analysis and Pole, Principal stressspace, Stress pathsin pq space; MohrCoulomb failure
More informationTHE DETERMINATION OF FRACTURE STRENGTH FROM ULTIMATE TENSILE AND TRANSVERSE RUPTURE STRESSES
Powder Metallurgy Progress, Vol.3 (003), No 3 119 THE DETERMINATION OF FRACTURE STRENGTH FROM ULTIMATE TENSILE AND TRANSVERSE RUPTURE STRESSES A.S. Wronski, A.Cias Abstract It is wellrecognized that the
More informationTechnical Note 16 Equivalent Static Method
Technical Note 16 Equivalent Static Method Contents Technical Note 21 ... 1 1 Introduction... 1 2 Operational Strain in the Pipeline... 2 3 Seismicity... 2 4 Vertical Uplift... 3 5 Vertical Bearing...
More informationMechanics of Materials MENG 270 Fall 2003 Exam 3 Time allowed: 90min. Q.1(a) Q.1 (b) Q.2 Q.3 Q.4 Total
Mechanics of Materials MENG 70 Fall 00 Eam Time allowed: 90min Name. Computer No. Q.(a) Q. (b) Q. Q. Q.4 Total Problem No. (a) [5Points] An air vessel is 500 mm average diameter and 0 mm thickness, the
More informationIntroduction to Fracture
Introduction to Fracture Introduction Design of a component Yielding Strength Deflection Stiffness Buckling critical load Fatigue Stress and Strain based Vibration Resonance Impact High strain rates Fracture
More informationAnswers: Internal Processes and Structures (Isostasy)
Answers: Internal Processes and Structures (Isostasy) 1. Analyse the adjustment of the crust to changes in loads associated with volcanism, mountain building, erosion, and glaciation by using the concept
More informationMain Means of Rock Stress Measurement
Real Stress Distributions through Sedimentary Strata and Implications for Reservoir Development and Potential Gas and Coal Development Strategies Ian Gray Sigra Pty Ltd 93 Colebard St West, Acacia Ridge,
More informationUnit III Theory of columns. Dr.P.Venkateswara Rao, Associate Professor, Dept. of Civil Engg., SVCE, Sriperumbudir
Unit III Theory of columns 1 Unit III Theory of Columns References: Punmia B.C.,"Theory of Structures" (SMTS) Vol II, Laxmi Publishing Pvt Ltd, New Delhi 2004. Rattan.S.S., "Strength of Materials", Tata
More information1.5 STRESSPATH METHOD OF SETTLEMENT CALCULATION 1.5 STRESSPATH METHOD OF SETTLEMENT CALCULATION
Module 6 Lecture 40 Evaluation of Soil Settlement  6 Topics 1.5 STRESSPATH METHOD OF SETTLEMENT CALCULATION 1.5.1 Definition of Stress Path 1.5. Stress and Strain Path for Consolidated Undrained Undrained
More informationElastic Properties of Solid Materials. Notes based on those by James Irvine at
Elastic Properties of Solid Materials Notes based on those by James Irvine at www.antonineeducation.co.uk Key Words Density, Elastic, Plastic, Stress, Strain, Young modulus We study how materials behave
More information6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa ( psi) and
6.4 A cylindrical specimen of a titanium alloy having an elastic modulus of 107 GPa (15.5 10 6 psi) and an original diameter of 3.8 mm (0.15 in.) will experience only elastic deformation when a tensile
More informationNumerical simulation of inclined piles in liquefiable soils
Proc. 20 th NZGS Geotechnical Symposium. Eds. GJ Alexander & CY Chin, Napier Y Wang & R P Orense Department of Civil and Environmental Engineering, University of Auckland, NZ. ywan833@aucklanduni.ac.nz
More informationTHE ROLE OF DELAMINATION IN NOTCHED AND UNNOTCHED TENSILE STRENGTH
THE ROLE OF DELAMINATION IN NOTCHED AND UNNOTCHED TENSILE STRENGTH M. R. Wisnom University of Bristol Advanced Composites Centre for Innovation and Science University Walk, Bristol BS8 1TR, UK M.Wisnom@bristol.ac.uk
More informationTechnical Report TR
SimulationBased Engineering Lab University of WisconsinMadison Technical Report TR201617 Using the Complementarity and Penalty Methods for Solving Frictional Contact Problems in Chrono: Validation
More informationTOWARD A HOLISTIC LOAD MODEL FOR STRUCTURES IN BROKEN ICE
POAC 15 Trondheim, Norway Proceedings of the 23 rd International Conference on Port and Ocean Engineering under Arctic Conditions June 1418, 2015 Trondheim, Norway TOWARD A HOLISTIC LOAD MODEL FOR STRUCTURES
More informationMethods of Interpreting Ground Stress Based on Underground Stress Measurements and Numerical Modelling
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2006 Methods of Interpreting Ground Stress Based on Underground Stress Measurements and
More informationForce and Stress. Processes in Structural Geology & Tectonics. Ben van der Pluijm. WW Norton+Authors, unless noted otherwise 1/9/ :35 PM
Force and Stress Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/9/2017 12:35 PM We Discuss Force and Stress Force and Units (Trigonometry) Newtonian
More information