MATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f


 Benjamin Doyle
 9 months ago
 Views:
Transcription
1 MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS f m f f m T m f m f f m v v + + = +
2 PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown in Figure. One common example of such a material is concrete, which is a mix of construction aggregate (stone or sand) in a cement matrix. This lab will deal with a little more exotic material than that. You will do calculations on a glass bead/epoxy particle composite with constituent properties given below. ν Glass poxy Glass = 70 GPa = 5GPa = ν = 0.3 poxy Figure. An example of a particle composite. One can not without great difficulty analyse this material in its complete configuration, analytically or with the aid of finite element calculations. This difficulty is of course attributed to the particles and the fact that they give rise to an inhomogeneous material. A way around this problem is to determine some effective properties and model the composite as a homogeneous solid with the effective properties instead. The effective mechanical properties of the homogenized material are a kind of mix, or weighted average, of the constituents mechanical properties. There are several analytical techniques suggested to find these effective properties, more or less accurate. This homogenization can also be performed using a finite element approach and the volume average equations from micro mechanics shown below.
3 Ω Ω To be able to perform the analysis one need to define a representative volume element (RV) that is a model of the material on the microscale. The RV will be inhomogeneous and needs to give a good representation of the micro structure, but it will also need to be able to represent the material on the macroscale in an average sense. It needs to be space filling, and have the same constituent orientation and volume fraction as the whole structure it represents. Figure 2 shows the RV for a material with spherical inclusions. Figure 2. Cut out of spherical inclusion in matrix  RV The finite element approach adopted in this lab is summarized in the following steps:. Apply prescribed displacements on the boundary of the RV. The displacements should correspond to homogeneous strain fields within th e RV. 2. Calculate the volume average of stresses and strains, and. 3. The volume averages can be thought of as the homogenized macroscopic stresses and strain and will give the effective properties using Hooke s law.
4 The relation between average strains and displacements on the boundary of the RV reads, where one needs to remember that shear strains like γ ij is equal to 2ε ij etc. For the particular case of a prescribed strain we then set, 0 on the boundary of the RV. Solving the boundary value problem gives the stresses σ k ij and strains ε k ij in each element k. Since the volume of each element V k is known, the average stresses and strains can be approximated by a discrete sum over all N elements where denotes the average stress in element k. The homogenized effective properties and ν can be found from Hooke s law. Here shown for macroscopic strain in the xdirection, but similar in the y and zdirections. This theory session sets the scene for the lab, now lets get started!
5 Start ANSYS via the ANSYS Product Launcher as usual, i.e. make sure you are working on the C: hard drive. You need to download the input files from the course web page. There are three files for this part of the lab. First the preprocessor file, and then one for axial strain loading in the xdirection and one for shear strain loading in the xzplane. Open up the preprocessor file in a text file editor. Look for the value for the parameter Rp. This is the radius of the particle inclusion. You should set it to a value between 3 and 8 mm. Based on Rp the volume fraction v f of particles should be calculated. How does the radius of the inclusion and the radius of the whole RV (Rm) relate to the volume fraction? Rp = v f = When you have changed the value for the radius (don t change Rm, just Rp!) save the file, and open it up in ANSYS via File > Read Input from Since the geometry is symmetric with respect to all the coordinate axes one could think that modelling only one eight (/8) would be enough, but this would unnecessarily complicate the formulation of boundary conditions. Here, half of the sphere is modelled. You will investigate two things in this part of the lab. First you will look at the response from a pure axial strain state at the boundary to determine the effective Young s modulus and Poisson s ratio of the RV. Secondly you will apply a shear strain to the boundary and determine the effective Shear modulus of the RV. The axial strain loading will be given, as described above, as a controlled displacement on the outer boundary, i.e. the curved boundary. To start the solver just read the file PartA_AxialStrainX.txt into ANSYS. When the solution is done, feel free to have a look at the results in the usual way, if you like. You will not use the graphical post
6 processing ability in this lab, but it is a good way to check if your finite element solver has done what you wanted it to do. To get some meaningful results you will need to work with something that is called lement Tables, which is just a tabular form of the results for each element. There will only be one result for each element and field variable, which might be a little surprising since the elements used in this analysis are 20 node brick elements and therefore should contain more than one value for each field variable (depends on the number of integration points in the element). What ANSYS does is that it averages the results in each element so you only get one unique element result for each field variable. A slight loss of accuracy, but we will have to live with that. If we know the volume average strains in two directions and the volume average stresses we can calculate the effective properties using Hooke s law. Verify these equations by considering Hooke s law with 0. To find the volume average results needed in the equations above go to ANSYS Main Menu > General Postproc > lement Table > Define Table to bring up the lement Table Data screen. Here you will need to Add the results you want to average. By highlighting one result set at a time from the list and clicking on apply you will create an element table for each of these results. Don t forget to add a table for the element volumes. When you are done adding results just close the lement Table Data screen. Now you will need to multiply each element result with the volume of that particular element. To do this, click on General Postproc > lement Table > Multiply to bring up
7 the Multiply lement Table Items screen. Multiply all your result tables with the volume table. This will create a new table that you will need to give a unique name so you know which results it contains. See picture below for explanation. Here the element table results for shear stress SXZ is being multiplied with the volume table VOLU. Notice the user label for the new table! When you have all your results multiplied with the volume you will need to sum them all up (this operation represents the integral over the volume). This is done by clicking on Sum of ach Item. Now you will need to divide each sum with the summed total volume, thus finalizing the volume average operation. You can use the result sheet at the end to record your values and to calculate the effective properties. For the second part just clear the model and reload the geometry using the preprocessor file again. This time you will solve the same model but with a shear strain loaded boundary instead. This is done by reading in the solver file PartA_ShearStrainXZ.txt that applies the following displacement field to the outer boundary. 2 2
8 In the same manner as the previous part you will now need to find the effective shear modulus. Hooke s law for the shear modulus is a lot simpler, so you will only need two results to determine it. In linear isotropic elasticity there is a well known relation between the Young s modulus, the Poisson s ratio and the shear modulus. Check if this still applies between your effective properties? 2 Finally you will have a look at the Reuss bound for compliance and Voigt s bound for stiffness. These bounds define maximum and minimum values for the effective shear modulus and bulk modulus. You will only look at the bounds for the shear modulus here. The Reuss bound is; and the Voigt s bound is: Ω Ω Can you figure out how to calculate these bounds given what you know of the volume fractions of particles and the shear modulus of the particle and the matrix? A hint is that the equations are quite simple, and are similar in style to the Rule of mixtures equations. Did your previously calculated effective shear modulus fall with in the bounds? It should have
9 PART B HXAGONAL R.V.. Consider a unidirectional fibre composite with the fibres arranged in a hexagonal pattern, as shown in Figure 3a below. The representative volume element would, as it might be easy to guess, have a hexagonal shape with the fibre in the middle. A schematic figure of what the RV would look like is shown in Figure 3b. y x z Figure 3. a) The hexagonal pattern of the fibre composite. b) The representative volume element. Once again your task is to find the effective properties of this composite by using the micro mechanics approach. This time it will be much harder to do this from an analytical formula like the one derived from Hooke s law in part A. Why is that do you think? What is the major difference between the directional dependence of the RV in part A and in this part? The composite is a carbon fibre reinforced epoxy with properties given below. ν Carbon poxy Carbon = 250 GPa = 5GPa = ν = 0.3 poxy First you will need to get 7 files from the course home page. The preprocessor file contains the geometry and the mesh, as usual. The other 6 files correspond to the 6 possible strain loadings, all according to the name of the file. Before you read in the geometry into ANSYS you should open up that file in a text editor and change the value
10 of the parameter r, the radius of the fibre. Pick any value between 3 and 7 mm (you might get a warning about bad shapes of some elements ignore that). r = v f = Now save the file and read in the geometry in ANSYS in the usual way. Consider the following matrix equation for the volume averages of the stresses and strains. So if you apply each strain component one at a time you can look at the resulting stress components, and in that way calculate the corresponding column in the stiffness matrix. Since the definition and multiplication of the element results tables are somewhat time consuming those operations has been included in the solver text files for simplicity. All you will have to do after each solution is to click on Sum of ach Item to bring up the numerical results. Now solve the model by applying strain in the xdirection to the boundary, i.e. read in the solver file PartB_X.txt in ANSYS. Look at the summed results. For each solution you do, 5 of the 6 strain results are very close to zero, and from an analytical point of view they are exactly zero, and that is what you will assume them to be.
11 So, for example, for strain loading on the boundary in the zdirectio, the matrix equation above would reduce to the following: Thus, by some easily performed calculations you have the entire third column of the stiffness matrix. The stresses multiplied by the element volumes are named according to; S_X_V for the stresses in the xdirection and S_YZ_V for the shear stress in the yzplane. Similarly for the strains multiplied with the element volumes you will have; PS_X_V, and PS_YZ_V for the corresponding axial and shear strains. The height, or the length, of the RV along the zaxis will have an effect on the results. This is not really what you would want for an accurate model, but it can be good to have seen this behaviour. You will for example be able to see effects of this when you look at the longitudinal stresses after prescribing a longitudinal displacement on the boundary. Look at the edge between the fibre and the matrix, and you should see a ringlike stress gradient. Preferably, that should not be there. But to avoid this effect all together you would need to change the boundary conditions to something called periodic boundary conditions, and they are a little bit tricky to apply and more than we have time to do here. Now solve for all the different load cases, one at a time, and start to fill out the stiffness matrix with you calculated values. It will be sufficient if you record the values in GPa with only one decimal place (the program will give you the stresses measured in MPa). After a few solutions, can you see a pattern emerging in the stiffness matrix? What kind of symmetries would you expect given the geometry of the RV?
12 In chapter 5.4 of Gudmundsons book Material Mechanics several material symmetries are presented. Do any of them seem to fit your results? So what about those effective properties then? How do you find them now? If you consider the compliance matrix of a transversely isotropic material given below you will see that it is quite easy now. Remember that the compliance matrix S is just the inverse of the stiffness matrix C Record your values for the effective properties on the answer sheet at the end. There are several alternatives to performing this kind of analysis to get the effective properties. You will compare your results to two of them. The first is what is called the Rule of mixtures and the second is an empirical set of relations called the HalpinTsai equations. Rule of mixtures gives the following equations for the effective properties:
13 HalpinTsai predicts the same longitudinal stiffness and Poisson s ratio as Rule of mixtures, but for the other two, the equations are as follows: ; ; ; 2 ;
14 PART C DOS LNGTH MATTR? The most common model of a fibre composite is the one where the fibres are as long as the whole structure. The fibres go from end to end so to speak. In many common applications this is usually not true. One example is short fibre glass/polyester composite, that are very common in many applications today. In this part you will look at a unidirectional fibre composite with short fibres embedded in a matrix, see Figure 4 below. ν Glass Polyester Glass = 70 GPa = 2.5 GPa = ν = 0.3 Polyester Figure 4. A short fibre composite. The question you will investigate here is if this has any effect on the effective mechanical properties of the composite? You will only look at the longitudinal stiffness in this part. The model of the representative volume element can be very much simplified by using symmetries. First realize that the geometry is symmetric with respect to the longitudinal axis of the fibre, so an axisymmetric 2D model is applicable. Secondly there is also symmetry with respect to the middle of the fibre, so a symmetry condition can be applied there as well, effectively cutting the model in half. A schematic of the model can be seen in Figure 5.
15 y R L/2 L f /2 x R/2 Figure 5. The model geometry. Rotational symmetry is imposed along the yaxis and ordinary symmetry along the xaxis. Download the single input file from the course webpage. Open up the input file and change the value of Lf, the fibre length, to a value between 20 and 80 mm. Lf = v f = Save the file and open it in ANSYS, this will initiate both the preprocessor and the solver. The loading in this part will be a prescribed displacement on the boundary in the ydirection. If you look at the solution for the stress in the ydirection, you will notice that there is a stress gradient at the end of the fibre. The stresses in the fibre needs some distance to build up, so the fibre is not used in an optimal way.
16 Because of this the estimated effective properties calculated by the Rule of mixtures are quite off. This can be corrected by assuming that there is an ineffective fibre length L i. The approach to estimating this is outlined below: Find the stress in the ydirection at the node located at position (0,0). Provided that the fibre is not too short, this is the stress that should appear in the whole fibre if it was used optimally. To do this, use the List Results > Nodal Solution. The node you are looking for has node number. Record the stress in the ydirection as 0,0 on the answer sheet. Now calculate the volume average of in the whole fibre. To do this use element tables, where you can multiply and sum element results. You will however need to select just the elements of the fibre; otherwise you will average the whole model instead of just the fibre. Type in the command: CMSL,S,FIBR. This will select only the predefined area component corresponding to just the fibre part of the model. Now you need to select all the elements in that area by typing: SLA. This selects the elements from an already selected area. Now calculate the volume average stresses of the fibre in the ydirection. Ω f The difference between the two stresses is essentially a measure of the effective fibre length, where denotes the so called ineffective stress. An average stress that is smaller than 0,0 can be seen as the stress 0,0 acting on the effective fibre length, and a vanishing stress acting on the remaining part of the fibre. Thus, 0,0 Now you solve this for the ineffective fibre length and record it on the answer sheet. What does that length measure actually mean?
17 That was all for this time, and for this course. Thanks for your participation. Good luck on the exam!
18 RSULTS FROM PART A SPHRICAL PARTICL INCLUSION Rp = v f = 2??
19 RSULTS FROM PART B HXAGONAL R.V.. r = v f = [ C] = To invert C you will probably need to use a numerical tool like Matlab or an advanced calculator. FM RoM HalpinTsai L L L T T T ν LT ν LT ν LT G LT G LT G LT
20 RSULTS FROM PART C DOS LNGTH MATTR? Lf = v f = 0,0
MATERIAL MECHANICS, SE2126 COMPUTER LAB 3 VISCOELASTICITY. k a. N t
MATERIAL MECHANICS, SE2126 COMPUTER LAB 3 VISCOELASTICITY N t i Gt () G0 1 i ( 1 e τ = α ) i= 1 k a k b τ PART A RELAXING PLASTIC PAPERCLIP Consider an ordinary paperclip made of plastic, as they more
More informationMATERIAL MECHANICS, SE2126 COMPUTER LAB 2 PLASTICITY
MATERIAL MECHANICS, SE2126 COMPUTER LAB 2 PLASTICITY PART A INTEGRATED CIRCUIT An integrated circuit can be thought of as a very complex maze of electronic components and metallic connectors. These connectors
More informationModule 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains
Introduction In this lecture we are going to introduce a new micromechanics model to determine the fibrous composite effective properties in terms of properties of its individual phases. In this model
More informationComposite FEM Labwork
Composite FEM Labwork You may perform these exercises in groups of max 2 persons. You may also between exercise 5 and 6. Be critical on the results obtained! Exercise 1. Open the file exercise1.inp in
More informationModule 2: Thermal Stresses in a 1D Beam Fixed at Both Ends
Module 2: Thermal Stresses in a 1D Beam Fixed at Both Ends Table of Contents Problem Description 2 Theory 2 Preprocessor 3 Scalar Parameters 3 Real Constants and Material Properties 4 Geometry 6 Meshing
More informationEffects of geometry and properties of fibre and matrix on overall. composite parameters
Int. Journal of Applied Sciences and ngineering Research, ol. 3, Issue 4, 2014 www.ijaser.com 2014 by the authors Licensee IJASR Under Creative Commons License 3.0 editorial@ijaser.com Research article
More informationMicromeso draping modelling of noncrimp fabrics
Micromeso draping modelling of noncrimp fabrics Oleksandr Vorobiov 1, Dr. Th. Bischoff 1, Dr. A. Tulke 1 1 FTA Forschungsgesellschaft für Textiltechnik mbh 1 Introduction Noncrimp fabrics (NCFs) are
More informationUnderstand basic stressstrain response of engineering materials.
Module 3 Constitutive quations Learning Objectives Understand basic stressstrain response of engineering materials. Quantify the linear elastic stressstrain response in terms of tensorial quantities
More informationModule 10: Free Vibration of an Undampened 1D Cantilever Beam
Module 10: Free Vibration of an Undampened 1D Cantilever Beam Table of Contents Page Number Problem Description Theory Geometry 4 Preprocessor 6 Element Type 6 Real Constants and Material Properties 7
More informationISSN: ISO 9001:2008 Certified International Journal of Engineering Science and Innovative Technology (IJESIT) Volume 2, Issue 4, July 2013
Delamination Studies in FibreReinforced Polymer Composites K.Kantha Rao, Dr P. Shailesh, K. Vijay Kumar 1 Associate Professor, Narasimha Reddy Engineering College Hyderabad. 2 Professor, St. Peter s Engineering
More informationModule 4: Behaviour of a LaminaeII. Learning Unit 1: M1. M4.1 Mechanics of Composites. M4.1.1 Introduction to Mechanics of Composites
Module 4: Behaviour of a LaminaeII Learning Unit 1: M1 M4.1 Mechanics of Composites M4.1.1 Introduction to Mechanics of Composites The relation between ply uniaxial strengths and constituent properties
More information3D and Planar Constitutive Relations
3D and Planar Constitutive Relations A School on Mechanics of Fibre Reinforced Polymer Composites Knowledge Incubation for TEQIP Indian Institute of Technology Kanpur PM Mohite Department of Aerospace
More informationDESIGN OF LAMINATES FOR INPLANE LOADING
DESIGN OF LAMINATES FOR INPLANOADING G. VERCHERY ISMANS 44 avenue F.A. Bartholdi, 72000 Le Mans, France Georges.Verchery@m4x.org SUMMARY This work relates to the design of laminated structures primarily
More information6. Bending CHAPTER OBJECTIVES
CHAPTER OBJECTIVES Determine stress in members caused by bending Discuss how to establish shear and moment diagrams for a beam or shaft Determine largest shear and moment in a member, and specify where
More informationCHEMC2410: Materials Science from Microstructures to Properties Composites: basic principles
CHEMC2410: Materials Science from Microstructures to Properties Composites: basic principles Mark Hughes 14 th March 2017 Today s learning outcomes To understand the role of reinforcement, matrix and
More informationMechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA
Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA What programs are in PROMAL? Master Menu The master menu screen with five separate applications from
More informationDue Monday, September 14 th, 12:00 midnight
Due Monday, September 14 th, 1: midnight This homework is considering the analysis of plane and space (3D) trusses as discussed in class. A list of MatLab programs that were discussed in class is provided
More informationEXPERIMENT 4: AN ELECTRICALTHERMAL ACTUATOR
EXPERIMENT 4: AN ELECTRICALTHERMAL ACTUATOR 1. OBJECTIVE: 1.1 To analyze an electricalthermal actuator used in a microelectromechanical system (MEMS). 2. INTRODUCTION 2.1 Introduction to Thermal Actuator
More informationA FINITE ELEMENT MODEL TO PREDICT MULTI AXIAL STRESSSTRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE
A FINITE ELEMENT MODEL TO PREDICT MULTI AXIAL STRESSSTRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE Daxu Zhang and D. R. Hayhurst School of Mechanical, Aerospace and Civil Engineering,
More informationCOURSE STE6289 Modern Materials and Computations (Moderne materialer og beregninger 7.5 stp.)
Narvik University College (Høgskolen i Narvik) EXAMINATION TASK COURSE STE6289 Modern Materials and Computations (Moderne materialer og beregninger 7.5 stp.) CLASS: Master students in Engineering Design
More informationLab Exercise #3: Torsion
Lab Exercise #3: Prelab assignment: Yes No Goals: 1. To evaluate the equations of angular displacement, shear stress, and shear strain for a shaft undergoing torsional stress. Principles: testing of round
More informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain  Axial Loading
MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain  Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain  Axial Loading Statics
More informationHomogenization on MultiMaterials Elements: Application to Printed Circuit Boards and Warpage Analysis
Homogenization on MultiMaterials Elements: Application to Printed Circuit Boards and Warpage Analysis Manuel Araújo 1,a, J. L. Alves 1, Paulo Silva 1 and Pedro Delgado 2 1 MEMS, University of Minho, Dep.
More informationPlane and axisymmetric models in Mentat & MARC. Tutorial with some Background
Plane and axisymmetric models in Mentat & MARC Tutorial with some Background Eindhoven University of Technology Department of Mechanical Engineering Piet J.G. Schreurs Lambèrt C.A. van Breemen March 6,
More informationA Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials
Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical
More informationLAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS
XXII. LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS Introduction The lamination theory for the elastic stiffness of fiber composite materials is the backbone of the entire field, it holds
More informationQUESTION BANK Composite Materials
QUESTION BANK Composite Materials 1. Define composite material. 2. What is the need for composite material? 3. Mention important characterits of composite material 4. Give examples for fiber material 5.
More informationσ = F/A ε = L/L σ ε a σ = Eε
Material and Property Information This chapter includes material from the book Practical Finite Element Analysis. It also has been reviewed and has additional material added by Sascha Beuermann. Hooke
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 informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More 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 informationPredicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics
Predicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics Neelam Mishra 1, and Kaushik Das *1 1 School of Minerals Metallurgical and Materials Engineering, Indian
More informationMESH MODELING OF ANGLEPLY LAMINATED COMPOSITE PLATES FOR DNS AND IPSAP
16 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS MESH MODELING OF ANGLEPLY LAMINATED COMPOSITE PLATES FOR DNS AND IPSAP Wanil Byun*, Seung Jo Kim*, Joris Wismans** *Seoul National University, Republic
More informationResponse Spectrum Analysis Shock and Seismic. FEMAP & NX Nastran
Response Spectrum Analysis Shock and Seismic FEMAP & NX Nastran Table of Contents 1. INTRODUCTION... 3 2. THE ACCELEROGRAM... 4 3. CREATING A RESPONSE SPECTRUM... 5 4. NX NASTRAN METHOD... 8 5. RESPONSE
More informationOptimal Shape and Topology of Structure Searched by Ants Foraging Behavior
ISSN 03861678 Report of the Research Institute of Industrial Technology, Nihon University Number 83, 2006 Optimal Shape and Topology of Structure Searched by Ants Foraging Behavior Kazuo MITSUI* ( Received
More informationTHEME A. Analysis of the elastic behaviour of La Aceña archgravity dam
THEME A Analysis of the elastic behaviour of La Aceña archgravity dam Gjorgi KOKALANOV, Professor, Faculty of Civil Eng., Skopje, Republic of Macedonia Ljubomir TANČEV, Professor, Faculty of Civil Eng.,
More informationCONSTITUTIVE RELATIONS FOR LINEAR ELASTIC SOLIDS
Chapter 9 CONSTITUTIV RLATIONS FOR LINAR LASTIC SOLIDS Figure 9.1: Hooke memorial window, St. Helen s, Bishopsgate, City of London 211 212 CHAPTR 9. CONSTITUTIV RLATIONS FOR LINAR LASTIC SOLIDS 9.1 Mechanical
More informationME 582 Advanced Materials Science. Chapter 2 Macromechanical Analysis of a Lamina (Part 2)
ME 582 Advanced Materials Science Chapter 2 Macromechanical Analysis of a Lamina (Part 2) Laboratory for Composite Materials Research Department of Mechanical Engineering University of South Alabama, Mobile,
More informationAn Evaluation and Comparison of Models for Maximum Deflection of Stiffened Plates Using Finite Element Analysis
Marine Technology, Vol. 44, No. 4, October 2007, pp. 212 225 An Evaluation and Comparison of Models for Maximum Deflection of Stiffened Plates Using Finite Element Analysis Lior Banai 1 and Omri Pedatzur
More informationAnalysis of high loss viscoelastic composites
Analysis of high loss viscoelastic composites by C. P. Chen, Ph.D. and R. S. Lakes, Ph.D. Department of Engineering Physics Engineering Mechanics Program; Biomedical Engineering Department Materials Science
More informationMATERIAL ELASTIC ANISOTROPIC command
MATERIAL ELASTIC ANISOTROPIC command.. Synopsis The MATERIAL ELASTIC ANISOTROPIC command is used to specify the parameters associated with an anisotropic linear elastic material idealization. Syntax The
More informationMatrix Assembly in FEA
Matrix Assembly in FEA 1 In Chapter 2, we spoke about how the global matrix equations are assembled in the finite element method. We now want to revisit that discussion and add some details. For example,
More informationDue Monday, October 19 th, 12:00 midnight. Problem 1 2D heat conduction in a rectangular domain (hand calculation) . A constant heat source is given:
Due Monday, October 19 th, 12:00 midnight Problem 1 2D heat conduction in a rectangular domain (hand calculation) Consider a problem on a rectangular (2 m 1 m) domain as shown in the left figure. The 1
More informationStress and Displacement Analysis of a Rectangular Plate with Central Elliptical Hole
Stress and Displacement Analysis of a Rectangular Plate with Central Elliptical Hole Dheeraj Gunwant, J. P. Singh mailto.dheerajgunwant@gmail.com, jitenderpal2007@gmail.com, AIT, Rampur Abstract A static
More informationPhysics 202, Lecture Today s Topics Middle T erm Term 1 Review
Physics 202, Lecture 7 Today s Topics Middle Term 1 Review About Exam 1 When and where Monday Sept. 27 th 5:307:00 pm 2301, 2241 Chamberlin (room allocation to be announced) nced) Format Close book One
More informationToday s agenda: Capacitors and Capacitance. You must be able to apply the equation C=Q/V.
Today s agenda: Capacitors and Capacitance. You must be able to apply the equation C=Q/V. Capacitors: parallel plate, cylindrical, spherical. You must be able to calculate the capacitance of capacitors
More informationComposites Design and Analysis. Stress Strain Relationship
Composites Design and Analysis Stress Strain Relationship Composite design and analysis Laminate Theory Manufacturing Methods Materials Composite Materials Design / Analysis Engineer Design Guidelines
More informationAE3610 Experiments in Fluid and Solid Mechanics TRANSIENT MEASUREMENTS OF HOOP STRESSES FOR A THINWALL PRESSURE VESSEL
Objective AE3610 Experiments in Fluid and Solid Mechanics TRANSIENT MEASUREMENTS OF OOP STRESSES FOR A TINWA PRESSURE VESSE This experiment will allow you to investigate hoop and axial stress/strain relations
More informationELASTICITY (MDM 10203)
LASTICITY (MDM 10203) Lecture Module 5: 3D Constitutive Relations Dr. Waluyo Adi Siswanto University Tun Hussein Onn Malaysia Generalised Hooke's Law In one dimensional system: = (basic Hooke's law) Considering
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 informationSOME RESEARCH ON FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS
Mechanical Testing and Diagnosis ISSN 2247 9635, 2012 (II), Volume 3, 7985 SOME RESEARCH ON FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS Valeriu DULGHERU, Viorel BOSTAN, Marin GUŢU Technical University
More informationBuckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate
Buckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate Outline Introduction Representative Volume Element (RVE) Periodic Boundary Conditions on RVE Homogenization Method Analytical
More informationBIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS
BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS H. Kumazawa and T. Takatoya Airframes and Structures Group, Japan Aerospace Exploration Agency 6131, Ohsawa, Mitaka,
More informationGeology 229 Engineering Geology. Lecture 5. Engineering Properties of Rocks (West, Ch. 6)
Geology 229 Engineering Geology Lecture 5 Engineering Properties of Rocks (West, Ch. 6) Common mechanic properties: Density; Elastic properties:  elastic modulii Outline of this Lecture 1. Uniaxial rock
More informationStress analysis of a stepped bar
Stress analysis of a stepped bar Problem Find the stresses induced in the axially loaded stepped bar shown in Figure. The bar has crosssectional areas of A ) and A ) over the lengths l ) and l ), respectively.
More informationMULTISCALE MODELING OF TRIAXIAL WOVEN FABRICS. APPLICATION TO SATELLITE ANTENNA
MULTISCALE MODELING OF TRIAXIAL WOVEN FABRICS. APPLICATION TO SATELLITE ANTENNA JEANSEBASTIEN GERARD 1, STEPHANE CATALAO 2, ROGER ASSAKER 2, JULIAN SANTIAGO PROWALD 3 1 exstream engineering SA, Rue
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module  01 Lecture  11
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Module  01 Lecture  11 Last class, what we did is, we looked at a method called superposition
More informationCHAPTER 6 BENDING Part 1
Ishik University / Sulaimani Civil Engineering Department Mechanics of Materials CE 211 CHAPTER 6 BENDING Part 11 CHAPTER 6 Bending Outlines of this chapter: 6.1. Chapter Objectives 6.2. Shear and
More informationElastic parameters prediction under dynamic loading based on the. unit cell of composites considering end constraint effect
Elastic parameters prediction under dynamic loading based on the unit cell of composites considering end constraint effect Wang Meng 1,, Fei Qingguo 1,, Zhang Peiwei 1, (1. Institute of Aerospace Machinery
More informationUSING A HOMOGENIZATION PROCEDURE FOR PREDICTION OF MATERIAL PROPERTIES AND THE IMPACT RESPONSE OF UNIDIRECTIONAL COMPOSITE
Volume II: Fatigue, Fracture and Ceramic Matrix Composites USING A HOMOGENIZATION PROCEDURE FOR PREDICTION OF MATERIAL PROPERTIES AND THE IMPACT RESPONSE OF UNIDIRECTIONAL COMPOSITE A. D. Resnyansky and
More informationEE C245 ME C218 Introduction to MEMS Design
EE C245 ME C218 Introduction to MEMS Design Fall 2007 Prof. Clark T.C. Nguyen Dept. of Electrical Engineering & Computer Sciences University of California at Berkeley Berkeley, CA 94720 Lecture 16: Energy
More informationEffects of Random CrossSectioned Distributions, Fiber Misalignment and Interphases in Three Dimensional Composite Models on Transverse Shear Modulus
University of South Florida Scholar Commons Graduate Theses and Dissertations Graduate School January 2012 Effects of Random CrossSectioned Distributions, Fiber Misalignment and Interphases in Three
More informationA Finite Difference Implementation of Phase Field Theory
University of Connecticut DigitalCommons@UConn Master's Theses University of Connecticut Graduate School 12162012 A Finite Difference Implementation of Phase Field Theory Nicholas W. Oren University
More informationTutorial 2. SSMA Cee in Compression: 600S F y = 50ksi Objective. A the end of the tutorial you should be able to
CUFSM 2.5 Tutorial 2 SSMA Cee in Compression: 600S20033 F y = 50ksi Objective To model a typical Cee stud in compression and determine the elastic critical local buckling load (P crl )and elastic critical
More informationContent. Department of Mathematics University of Oslo
Chapter: 1 MEK4560 The Finite Element Method in Solid Mechanics II (January 25, 2008) (Epost:torgeiru@math.uio.no) Page 1 of 14 Content 1 Introduction to MEK4560 3 1.1 Minimum Potential energy..............................
More informationNanoindentation of Fibrous Composite Microstructures: Experimentation and Finite Element Investigation. Mark Hardiman
Nanoindentation of Fibrous Composite Microstructures: Experimentation and Finite Element Investigation Mark Hardiman Materials and Surface Science Institute (MSSI), Department of Mechanical and Aeronautical
More informationAnalysis of Planar Truss
Analysis of Planar Truss Although the APES computer program is not a specific matrix structural code, it can none the less be used to analyze simple structures. In this example, the following statically
More informationInfluence of the filament winding process variables on the mechanical behavior of a composite pressure vessel
Influence of the filament winding process variables on the mechanical behavior of a composite pressure vessel G. Vargas 1 & A. Miravete 2 1 Universidad Pontificia Bolivariana, Facultad de Ingeniería Mecánica,
More informationCalculation of Energy Release Rate in Mode I Delamination of Angle Ply Laminated Composites
Copyright c 2007 ICCES ICCES, vol.1, no.2, pp.6167, 2007 Calculation of Energy Release Rate in Mode I Delamination of Angle Ply Laminated Composites K. Gordnian 1, H. Hadavinia 1, G. Simpson 1 and A.
More informationBending of Simply Supported Isotropic and Composite Laminate Plates
Bending of Simply Supported Isotropic and Composite Laminate Plates Ernesto GutierrezMiravete 1 Isotropic Plates Consider simply a supported rectangular plate of isotropic material (length a, width b,
More informationTunnel Reinforcement Optimization for Nonlinear Material
November 2527, 2012, Gold Coast, Australia www.iccm2012.org Tunnel Reinforcement Optimization for Nonlinear Material T. Nguyen* 1,2, K. Ghabraie 1,2, T. TranCong 1,2 1 Computational Engineering and
More informationMechanical Behavior of Fullerene Reinforced Fiber Composites with Interface Defects through Homogenization Approach and Finite Element Method
, pp.6782 http://dx.doi.org/1.14257/ijast.215.78.6 Mechanical Behavior of Fullerene Reinforced Fiber Composites with Interface Defects through Homogenization Approach and Finite Element Method P. Prasanthi
More informationThe following syntax is used to describe a typical irreducible continuum element:
ELEMENT IRREDUCIBLE T7P0 command.. Synopsis The ELEMENT IRREDUCIBLE T7P0 command is used to describe all irreducible 7node enhanced quadratic triangular continuum elements that are to be used in mechanical
More informationBack Analysis of Measured Displacements of Tunnels
Rock Mechanics and Rock Engineering 16, 173180 (1983) Rock Mechanics and Rock Engineering 9 by SpringerVerlag 1983 Back Analysis of Measured Displacements of Tunnels By S. Sakurai and K. Takeuchi Kobe
More informationAssessment of boundaryelement method for modelling the structural response of a pavement
Assessment of boundaryelement method for modelling the structural response of a pavement (PhD Student, Lecturer DEC of FCTUC) (Full Professor  IST) Pavements: Materials, design and performance a, LNEC,
More informationMath 5a Reading Assignments for Sections
Math 5a Reading Assignments for Sections 4.1 4.5 Due Dates for Reading Assignments Note: There will be a very short online reading quiz (WebWork) on each reading assignment due one hour before class on
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ocw.mit.edu 12.510 Introduction to Seismology Spring 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Stress and Strain
More informationThe Power of Ultrasonic Characterisation for Completely Assessing the Elastic Properties of Materials
11th European Conference on NonDestructive Testing (ECNDT 214), October 61, 214, Prague, Czech Republic The Power of Ultrasonic Characterisation for Completely Assessing the Elastic Properties of Materials
More informationMECHANICS OF MATERIALS
Third E CHAPTER 2 Stress MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf Lecture Notes: J. Walt Oler Texas Tech University and Strain Axial Loading Contents Stress & Strain:
More informationLAMINATED COMPOSITE PLATES
LAMINATED COMPOSITE PLATES David Roylance Department of Materials Science and Engineering Massachusetts Institute of Technology Cambridge, MA 02139 February 10, 2000 Introduction This document is intended
More information4.MECHANICAL PROPERTIES OF MATERIALS
4.MECHANICAL PROPERTIES OF MATERIALS The diagram representing the relation between stress and strain in a given material is an important characteristic of the material. To obtain the stressstrain diagram
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 informationBasic Concepts of Strain and Tilt. Evelyn Roeloffs, USGS June 2008
Basic Concepts of Strain and Tilt Evelyn Roeloffs, USGS June 2008 1 Coordinates Righthanded coordinate system, with positions along the three axes specified by x,y,z. x,y will usually be horizontal, and
More informationDEPARTMENT OF MECHANICAL ENIGINEERING, UNIVERSITY OF ENGINEERING & TECHNOLOGY LAHORE (KSK CAMPUS).
DEPARTMENT OF MECHANICAL ENIGINEERING, UNIVERSITY OF ENGINEERING & TECHNOLOGY LAHORE (KSK CAMPUS). Lab Director: Coordinating Staff: Mr. Muhammad Farooq (Lecturer) Mr. Liaquat Qureshi (Lab Supervisor)
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 informationInternational Journal of Advanced Engineering Technology EISSN
Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,
More information3D Compression Molding
Autodesk Simulation Moldflow Insight 2014 3D Compression Molding Executive summary In this work, the simulation results from a program developed for the threedimensional analysis of compression molding
More informationLectures on. Constitutive Modelling of Arteries. Ray Ogden
Lectures on Constitutive Modelling of Arteries Ray Ogden University of Aberdeen Xi an Jiaotong University April 2011 Overview of the Ingredients of Continuum Mechanics needed in Soft Tissue Biomechanics
More informationHyperbolic Soil Bearing Capacity
1 Introduction Hyperbolic Soil Bearing Capacity This example involves analyzing the bearing capacity of a round footing. The example is useful for illustrating several SIGMA/W features and procedures,
More informationModelling Anisotropic, Hyperelastic Materials in ABAQUS
Modelling Anisotropic, Hyperelastic Materials in ABAQUS Salvatore Federico and Walter Herzog Human Performance Laboratory, Faculty of Kinesiology, The University of Calgary 2500 University Drive NW, Calgary,
More informationMODELLING OF THE THERMOMECHANICAL PROPERTIES OF WOVEN COMPOSITES DURING THE CURE
TH 9 TH INTRNATIONAL ONFRN ON OMPOSIT MATRIALS MODLLING OF TH THRMOMHANIAL PROPRTIS OF WOVN OMPOSITS DURING TH UR L. Khoun,, K.S. hallagulla,, P. Hubert Department of Mechanical ngineering, McGill University,
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 informationFinite Element Analysis of Composite Laminate By Using ABDH Matrix(Stiffness Matrix)
Finite Element Analysis of Composite Laminate By Using ABDH Matrix(Stiffness Matrix) Nikhil J. Chaudhari 1 Post Graduate Student Department of Mechanical Engineering Veermata Jijabai Technological Institute
More informationFINITE ELEMENT APPROACHES TO MESOSCOPIC MATERIALS MODELING
FINITE ELEMENT APPROACHES TO MESOSCOPIC MATERIALS MODELING Andrei A. Gusev Institute of Polymers, Department of Materials, ETHZürich, Switzerland Outlook Generic finite element approach (PALMYRA) Random
More informationMAXIMUM ENTROPYBASED UNCERTAINTY MODELING AT THE FINITE ELEMENT LEVEL. Pengchao Song and Marc P. Mignolet
MAXIMUM ENTROPYBASED UNCERTAINTY MODELING AT THE FINITE ELEMENT LEVEL Pengchao Song and Marc P. Mignolet SEMTE, Faculties of Mechanical and Aerospace Engineering, Arizona State University, 51 E. Tyler
More informationMultiscale Approach to Damage Analysis of Laminated Composite Structures
Multiscale Approach to Damage Analysis of Laminated Composite Structures D. Ivančević and I. Smojver Department of Aeronautical Engineering, Faculty of Mechanical Engineering and Naval Architecture, University
More information2. Mechanics of Materials: Strain. 3. Hookes's Law
Mechanics of Materials Course: WB3413, Dredging Processes 1 Fundamental Theory Required for Sand, Clay and Rock Cutting 1. Mechanics of Materials: Stress 1. Introduction 2. Plane Stress and Coordinate
More informationCHAPTER 3 THE EFFECTS OF FORCES ON MATERIALS
CHAPTER THE EFFECTS OF FORCES ON MATERIALS EXERCISE 1, Page 50 1. A rectangular bar having a crosssectional area of 80 mm has a tensile force of 0 kn applied to it. Determine the stress in the bar. Stress
More informationSurface stress and relaxation in metals
J. Phys.: Condens. Matter 12 (2000) 5541 5550. Printed in the UK PII: S09538984(00)113864 Surface stress and relaxation in metals P M Marcus, Xianghong Qian and Wolfgang Hübner IBM Research Center, Yorktown
More informationLaboratory 4 Bending Test of Materials
Department of Materials and Metallurgical Engineering Bangladesh University of Engineering Technology, Dhaka MME 222 Materials Testing Sessional.50 Credits Laboratory 4 Bending Test of Materials. Objective
More information