Evaluation of stress intensity factor for cracks in graded materials using digital image correlation
|
|
- Barbara Campbell
- 5 years ago
- Views:
Transcription
1 Evaluation of stress intensity factor for cracks in graded materials using digital image correlation Amit Kumar and Venkitanarayanan Parameswaran * Department of Mechanical Engineering Indian Institute of Technology Kanpur, Kanpur ABSTRACT This study focuses on the determination of stress intensity factor (SIF) for edge cracks in a functionally graded plate using the technique of Digital Image Correlation (DIC). Graded composites having nearly exponential variation of elastic modulus were prepared by gravity casting technique. The elastic modulus profile of the prepared material was characterized. The material had an elastic gradient of two times across a width of 50 mm. Edge cracked plates having different crack lengths were loaded in bending and the displacements were measured using DIC. The SIF was extracted by fitting the theoretical crack-tip displacement field available for homogeneous materials to the experimental data. The experimentally determined SIF for edge cracks subjected to four-point bending agreed well with published theoretical results. For cracks subjected to three-point bending, the results indicated that the SIF is higher than that for an identical crack in a homogeneous material when the crack is located in the stiff side of the graded plate whereas it is lower than that for homogeneous materials when the crack is located in the compliant side of the plate. Keywords: Digital image correlation, Functionally graded material, Stress intensity factor 1. INTRODUCTION Functionally graded materials (FGMs) are a special class of composite materials in which the composition and microstructure is not uniform in space, but varies in a gradual progressive manner. These materials were first introduced as high performance thermal barrier coatings, but they have much wider scope for application in many fields; for example aerospace, microelectronics and biomedical engineering. FGMs owe their success to two main features: the gradual change in composition allows better exploitation of the constituent phase properties and unlike traditional bimaterials junctions, FGMs avoid the abrupt coupling of heterogeneous phases. A gradual transition in material properties has been observed to alleviate the problem of thermal cracking, particularly in thermal protection systems. Moreover the composition gradient can be tailored to suite the application requirements, even when the loading conditions differ from point to point. FGMs are now being regarded as one of the most promising candidates for future multi-functional composites in many engineering sectors. Since the introduction of the FGM concept, the design and analysis of this material has received considerable attention. An FGM used for thermal protection is most likely to have a gradual transition from pure metal to pure ceramic and will combine the high temperature properties of the ceramic with the ductile behavior of the metal. However ceramic being brittle in nature, the system is fracture prone. In a layered metal-ceramic system, a crack present in the brittle layer upon reaching the interface will most likely deflect along the weak interface leaving the second layer intact. In a continuous transition without interfaces as in an FGM, the crack could propagate through the entire thickness causing complete failure of the system. Hence understanding the behavior of cracks in a graded material is important to evaluate the integrity of structures made of graded materials. Recent years have witnessed many investigations on the fracture mechanics of graded materials. These studies have unanimously established the inverse square root nature of stresses very close to the crack tip in graded materials; however, the stress intensity factor is influenced by the elastic gradients [1]. The theoretical values of stress intensity factor for embedded and edge cracks in infinite FGM plates are available [2]. While the theoretical study of fracture of * Further author information: (Send correspondence to Venkitanarayanan Parameswaran) Venkitanarayanan Parameswaran.: venkit@iitk.ac.in, Telephone:
2 FGMs is well established, experimental investigations are very few. FGM plates have been manufactured using epoxy and glass beads and their fracture behavior has been investigated by Butcher et al. [3]. Lambros et al. have experimentally investigated crack growth in polymeric FGMs having a non-linear constitutive response using the technique of digital image correlation [4]. The present study focuses on the application of the technique of Digital Image Correlation (DIC) to evaluate the stress intensity factor (SIF) for edge cracks in FGMs. DIC is a full field technique for measuring the surface displacements of deforming objects [5, 6]. DIC needs minimal surface preparation, is non-invasive and can be applied to any type of material. In the present study, FGMs with near exponential variation of elastic modulus were prepared and the displacement field around edge cracks loaded in bending was measured using DIC. The SIF was determined by fitting the theoretical displacement fields to the experimental data. The results for edge cracks subjected to pure bending are compared with published theoretical results. Further the SIF for edge cracks subjected to three point bending is also estimated. The details of the FGM preparation, characterization of the modulus profile, experimental setup and analysis of the experimental data to determine the SIF are presented in the following sections. 2. PREPARATION AND CHARACTERIZATION OF GRADED MATERIAL The graded material required for the study was prepared in the form of a particulate composite with continuously varying particle volume fraction. The matrix material used was an epoxy system consisting of the Araldite, LY556 and Hardener HY951, both supplied by Vantico Performance Polymers Pvt. Ltd., India. Glass beads of size µm supplied by Hindustan Glass Beads were used as the reinforcement. Adding rigid particles to the resin will improve the elastic modulus of the composite. In order to characterize this effect, homogeneous composites with glass bead volume fraction varying from 0% to 50% in steps of 10% were first prepared and their elastic properties were measured as per ASTM D 638. The elastic modulus of the composite with 50% glass beads was three times that of pure epoxy. FGMs sheets, of size 200 mm in length, 90 mm in width and 6.4 mm in thickness were prepared by a gravity casting technique. A homogeneous mixture of epoxy resin and glass beads in the required proportion was first prepared. After degassing the mixture in vacuum, the hardener in the required amount was added and mixed gently without trapping air bubbles. Simultaneously another mixture of epoxy and hardener was prepared. The mixture consisting epoxy, glass beads and hardener was first poured into an acrylic mould kept in vertical position. Then the mixture of epoxy and hardener was poured into the acrylic mould till the mold is completely filled and then the mold was sealed. After 20 minutes the mold was flipped upside down and kept in vertical position. Glass beads have a higher specific gravity of 2.2 as compared to 1.17 for the epoxy resin. When the mold is flipped upside down, the top layer of the mixture that is rich in glass beads starts moving down due to its higher specific gravity. Simultaneously the resin gels, its viscosity increases and this will arrest the settling of the glass beads. This results in a casting with a glass bead rich region at the bottom, a resin rich region at the top, and an intermediate region with continuously varying glass bead content. The sheet was allowed to cure in the mould at room temperature for 36 hours. After this the sheet was post cured in an air circulating oven for 4 hours at 100 o C. In order to determine the composition profile of the prepared FGM sheet, small samples of size 25 x 6.4 x 5 mm were sliced along the width of the sheet. The density of these samples was measured by hydrostatic weighing. The glass bead volume fraction (V f ) was calculated from the density using the rule of mixtures. The composition profile of a typical graded sheet is shown in figure 1. It can be noticed that there is a region of continuously varying glass bead content in the middle portion of the sheet. This portion was used for the fracture experiments. From the composition profile and the variation of the elastic modulus as function of the glass bead content, the elastic modulus profile of the graded sheet was determined. The elastic modulus profile for a typical sheet, shown in figure 2, indicates that the elastic modulus varies by a factor of two over a length of 50 mm. The continuous line in figure 2 is an exponential fit to the elastic modulus data. The composition profile and elastic modulus profile was characterized for each casting and this indicated that the gradation achieved by the gravity casting technique is fairly repeatable. Plates of length 170 mm, width 45 mm and thickness 6.4 mm were cut out of the prepared sheets. A notch of width 200 µm was first made on one of the edges using a water cooled diamond cutter. A natural crack of the required length was extended from this notch tip by forcing a sharp blade into the notch and tapping it in a controlled manner (see figure 3). The length of the crack was measured precisely using a traveling microscope equipped with a micrometer. During this measurement a slight crack face opening load was applied by a thin wedge to facilitate locating of the crack tip. By this way the crack length could be measured to an accuracy of ± 0.1 mm.
3 40 8 Vf (%) Distance (mm) Elastic modulus (GPa) E = 3.6e 0.016x Distance (mm) Figure 1. Composition profile of the graded sheet. Figure 2. Elastic modulus profile of the graded sheet. E 2 x x y r θ W y E 1 a Natural crack Notch Figure 3. Natural crack and axis convention used in the analysis Figure 4. Random dot pattern for DIC technique 3. EXPERIMENTAL DETAILS The DIC set up used in the study comprised of a 2 Mega pixel CCD camera, a zoom lens and the image acquisition and correlation software (Vic2D) supplied by Correlated Solutions Inc., USA. The technique of DIC involves applying a random black and white pattern on the surface of the specimen and then capturing digital images of this pattern before and after loading. The in plane displacement components are then calculated through intensity correlation of the images taken before and after deformation. The details of the correlation procedure can be found in [5]. A random dot pattern was applied to the specimen surface using black and white spray paint. A photograph of the random dot pattern is given in figure 4. The specimen was loaded in a 25 kn testing machine and the load was measured using a 1 kn load cell. A subset size of 34 x 34 pixels was used for the correlation. Proper care was taken to align the horizontal axis of the camera with the crack line. The u and v displacement maps from a typical test, shown in figure 5, indicate a smooth displacement field. The displacements shown in figure 5 include rigid body components (translation and rotation) that arise during the loading. Tests were conducted for different crack lengths for the case two cases 1) E 2 /E 1 =2 (crack is in the compliant side) and 2) E 2 /E 1 =0.5 (crack is in the stiff side). The displacement data from the correlation was exported and the stress intensity factor was calculated from the displacement field as explained in the next section. 4. EXTRACTION OF STRESS INTENSITY FACTOR The SIF was determined by fitting the asymptotic crack tip displacement field to the experimentally measured displacement field. Crack tip displacement field specific to graded materials are not yet available. Bueno and Lambros [4] have reported the use of the displacement field available for homogeneous materials for estimating SIF in graded materials. The expansion of the displacement field for homogeneous materials is the following.
4 Figure 5. In plane displacement contours (u & v) from DIC technique for four point bending, a/w=0.55 Load=40 N, E 2 /E 1 =2 1+ ν κ 1 κ + 1 u = Re Zn + ReYn ( yim Zn + yimyn) Ec n= 0 2 2, (1) 1+ ν κ + 1 κ 1 v= Im Zn + ImYn ( yre Zn + yreyn) Ec n= where κ = (3 ν) /(1 + ν) for plane stress, ν being the Poisson s ratio and E c is the elastic modulus at the crack tip. The functions Z n and Y n are defined as follows. n 1/ 2 n Z = Az ; Y = Bz ; Z = Zdz; Y = Ydz n n n n n n n n, (2) θ where z = x + iy = re i. The SIF K I is related to the coefficient A 0 as K1 = A0 2π. From the experimentally measured displacements (u e & v e ), the radial component of the displacement u re was calculated for each data point as ure = ue cosθ + ve sinθ, (3) where θ, ( π < θ < π) is the angle made by the position vector of the data point with the x-axis shown in figure 3. The experimentally measured displacements include rigid body contributions also. Using equation (1) and considering the rigid body translations T x and T y in x and y directions respectively, the theoretical displacement field can be written as ur = ucosθ + vsinθ + Txcosθ + Tysinθ (4) Note that in the polar coordinates, rigid body rotation in the x-y plane will not induce any radial displacement and hence does not appear in equation (3). Equation (3) has (N+2) unknowns where N is the number of terms used in the series expansion of the displacement (equation 1). The number of data points will be much more than the number of unknowns. Therefore the unknowns (A 0 A N, B 0.B M, T x, T y ) are determined by fitting equation (4) to the experimental data (equation 3) in a least square sense. Though any of the displacement components u or v could be used for extraction of the SIF, u r was used in this study as it is not sensitive to rigid body rotation of the specimen. In the present study, the displacements were measured over a region of 80 x 40 mm 2 around the crack-tip. One could use this entire data set in the analysis, however, this will add to computation time. Instead, the data points sampled from an annular region of inner radius R 1 and outer radius R 2 were used to extract the SIF. The inner radius R 1 has to be larger than half the plate thickness to avoid three-dimensional effects and hence R 1 was set to 3 mm. The outer radius R 2 was 10 mm. The SIF was determined by increasing the number of terms in the series (equation 1) starting from 2 terms (N=2). This procedure was continued until the addition of further terms did not produce any change in the SIF. As an additional confirmation, the displacement field, u r, was regenerated using the constants (A 0 A N, B 0.B M, T x, T y ) and was visually compared with the experimental displacement field, u re, to see the match between the two. Further discussion on these aspects is provided in the next section.
5 5. RESULTS AND DISCUSSION The stress intensity factor was determined for different crack lengths in the range of 0.3 < a/w < 0.75 and two different gradations corresponding to E 2 /E 1 = 0.5 and 2 for four-point bending and three-point bending. For each case, the SIF was determined as explained in the last section by increasing the number of terms in equation (1 & 2) until convergence. Figure 6, shows the pattern of convergence of SIF for four-point bending. The SIF is normalized as K / σ π a, where I 2 σ = 6 M / BW, M is the bending moment at the crack plane and B is the specimen thickness. It can be seen that at least 12 terms are required in equation (1) for convergence. Figure 7 gives the comparison of SIF for four-point bending for E 2 /E 1 = 2 with theoretical results available for exponential variation of elastic modulus [2]. The experimentally determined SIF is in good agreement with the theoretical values with a maximum deviation of 6%. The SIF for cracks subjected to three point bending is shown in figure 8. The continuous line in figure 8 is the SIF for a homogeneous material having the same geometry. The effect of material nonhomogeneity on SIF can be clearly seen in figure 8. The SIF for the case E 2 /E 1 =0.5 is higher than that for a homogeneous material where as the SIF is lower than that for a homogeneous material for E 2 /E 1 =2. When the crack is in the compliant side of the plate (E 2 /E 1 =2), the SIF is about 15% lower than that for a homogeneous material whereas when the crack is in the stiff side of the plate (E 2 /E 1 =0.5), the SIF is about 15% higher than that for a homogeneous material. This difference seems to be more prominent for crack lengths a/w < 0.5. When a/w=0.5 the crack tip is at the same location irrespective of whether the crack is on the stiff side or on the compliant side, however, the SIF is different for the two cases. A comparison of the experimentally measured displacement u re and that calculated using equation (1) is shown in figure 9. The continuous lines are the regenerated displacement field and the broken lines are the experimentally measured displacement. It should be noted here that the rigid body components are removed while plotting the displacement contours in figure 9. The contours regenerated using four terms (N=4) in equation (1), shown in figure 9 do not match well with the experimental contours. When the number of terms is increased to 12, the agreement between the measured displacement field and calculated displacement field improves considerably. It was pointed out earlier that the SIF value also converges for N=12 (see figure 6). The convergence of SIF along with the matching of the regenerated and experimental displacement fields is an indicator of the accuracy of the determined SIF. 6. CONCLUSIONS An experimental investigation into the effect of elastic gradient on the SIF for edge cracks in graded plates subjected to bending has been carried out. Graded sheets required for the study was prepared using a gravity casting technique. The technique of DIC was employed to make full field measurement of the displacement field. The SIF was calculated by fitting the theoretical displacement field to the measured data. The results of the study indicated that, the SIF is higher than that for a homogeneous material when the crack is located in the stiff side of the material. In a ceramic-metal FGMs used for thermal protection applications, when the crack is in the ceramic rich side; the SIF will be higher compared to when it is in the metallic side. Ceramics have much lower resistance to fracture and combined with the increase in the SIF due to the gradation, the chances of a crack in the ceramic side becoming critical is higher in such FGMs. KI/(σ πa) a/w=0.42 a/w=0.51 a/w= Number of terms (N) Figure 6. Convergence of SIF with number of terms (E 2 /E 1 =2). KI/(σ πa) Theoretical Expertiment a/w Figure 7. Stress intensity factor for four-point bending (E 2 /E 1 =2)
6 3.5 KI/(σ πa) E2/E1=0.5 E2/E1= a/w Figure 8. Stress intensity factor for edge cracks in graded sheets subjected to three point bending. N=4 N=12 Figure 9. Contours of regenerated and measured displacement field (E 2 /E 1 =0.5, a/w=0.7) for N=4 and 12. Displacements are in (µm). ACKNOWLEDGMENTS The authors acknowledge the financial support for procuring the DIC system used in this study through project from Aeronautical Research and Development Board through the project DARO/08/ /M/I. REFERENCES 1. F. Erdogan, Fracture Mechanics of Functionally Graded Materials, Comp. Engng., 5:7, pp , F. Erdogan, and B. H. Wu, The Surface Crack Problem for a Plate with Functionally Graded Properties, Jl. of App. Mech., 64, pp , R. J. Butcher, C. E. Rousseau, and H. V. Tippur, A Functionally Graded Particulate Composite: Preparation, Measurements and Failure Analysis, Acta Mater., 47:1, pp , J. A. Bueno and J. Lambros, Investigation of Crack Growth in Functionally Graded Materials using Digital Image Correlation, Engng. Frac. Mech., 69, pp , W. H. Peters and W. F. Ranson, Digital Imaging Techniques in Experimental Stress Analysis. Opt Engng; 21(3), pp , S. Yoneyama, Y. Morimota, and M. Takashi, Automatic Evaluation of Mixed mode Stress Intensity Factors Utilizing Digital Image Correlation, Strain, 42, pp , 2006.
Stress Intensity Factor calculation from displacement fields
Focused on Crack ip Fields Stress Intensity Factor calculation from displacement fields S. Beretta, L. Patriarca, S. Rabbolini Politecnico di Milano, Department of Mechanical Engineering, via La Masa 1,
More informationUsing Fiber Reinforced Polymer to Restore Deteriorated Structural Members
International Journal of Material and Mechanical Engineering, 01, 1: 1-7 - 1 - Published Online April 01 http://www.ijm-me.org Using Fiber Reinforced Polymer to Restore Deteriorated Structural Members
More informationSTRESS INTENSITY FACTOR CALCULATIONS FOR CRACKS EMANATING FROM BOLT HOLES IN A JET ENGINE COMPRESSOR DISC
ICAS2002 CONGRESS STRESS INTENSITY FACTOR CALCULATIONS FOR CRACKS EMANATING FROM BOLT HOLES IN A JET ENGINE COMPRESSOR DISC W. Beres, A.K. Koul 1 Institute for Aerospace Research, National Research Council
More informationState of Stress in Three Dimensions
State of Stress in Three Dimensions Theories of failure Introduction: Due to large numbers of examples of compound stresses met with in engineering practice, the cause of failure or permanent set under
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 2014 Quiz #1 March 3, 2014 Name: SOLUTION ID#: PS.:
More informationPrincipal Stresses: Interpreting Principal Stresses
Principal Stresses: This document talks about the principal stresses in a gross sense than its theoretical explanations. And, here I go... Most of the times we do wonder what use are these principal stresses
More informationEffect of Mechanical Mismatch on the Stress Intensity Factors of Inclined Cracks under Mode I Tension Loading
Effect of Mechanical Mismatch on the Stress Intensity Factors of Inclined Cracks under Mode I Tension Loading M.A Azlan 1,a and A.E Ismail 1,b 1 Faculty of Mechanical and Manufacturing Engineering, Universiti
More informationCombined mode I stress intensity factors of slanted cracks
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Combined mode I stress intensity factors of slanted cracks To cite this article: A.E Ismail et al 2017 IOP Conf. Ser.: Mater.
More informationExperimental Evaluation of the Markl Fatigue Methods and ASME Piping Stress Intensification Factors
Experimental Evaluation of the Markl Fatigue Methods and ASME Piping Stress Intensification Factors Chris Hinnant Paulin Research Group Houston, TX Tony Paulin Paulin Research Group Houston, TX 2008 ASME
More informationWhat Is More Applicable Data? Taking Advantage of FEATools
What Is More Applicable Data? Taking Advantage of FEATools Why More Applicable Data? Tee flexibility is set to 1.0 no flexibility! Why More Applicable Data? This note (below) from B31.3 Appendix D alerts
More informationTHE PROXIMITY OF MICROVIAS TO PTHs AND ITS IMPACT ON THE RELIABILITY
THE PROXIMITY OF MICROVIAS TO s AND ITS IMPACT ON THE RELIABILITY Anthony Primavera Manager - CSP Consortium Universal Instruments Corporation Binghamton, New York 13902-0825. Jaydutt Joshi Package Development
More informationSAQ KONTROLL AB Box 49306, STOCKHOLM, Sweden Tel: ; Fax:
ProSINTAP - A Probabilistic Program for Safety Evaluation Peter Dillström SAQ / SINTAP / 09 SAQ KONTROLL AB Box 49306, 100 29 STOCKHOLM, Sweden Tel: +46 8 617 40 00; Fax: +46 8 651 70 43 June 1999 Page
More informationSupplemental Material Optics formula and additional results
Supplemental Material Optics formula and additional results Fresnel equations We reproduce the Fresnel equations derived from Maxwell equations as given by Born and Wolf (Section 4.4.). They correspond
More informationPearson Connected Mathematics Grade 7
A Correlation of Pearson Connected Mathematics 2 2012 to the Common Core Georgia Performance s Grade 7 FORMAT FOR CORRELATION TO THE COMMON CORE GEORGIA PERFORMANCE STANDARDS (CCGPS) Subject Area: K-12
More informationETN Evo. Plastic Lined Magnetic Drive Centrifugal Pumps. ETN Evo 50 ETFE. Pompe S.r.l. Comply to : 2006/42/CE. Design to : sub - ISO 2858
ETN Evo Plastic Lined Magnetic Drive Centrifugal Pumps ETN Evo 50 ETFE Plastic and Fluoroplastic Lined Magnetic drive Horizontal - Single Stage - Centrifugal pumps Sub-ISO designed Lining: PP (Polypropylene),
More informationRandom Variables and Probability Distributions
Chapter 3 Random Variables and Probability Distributions Chapter Three Random Variables and Probability Distributions 3. Introduction An event is defined as the possible outcome of an experiment. In engineering
More informationProf. B V S Viswanadham, Department of Civil Engineering, IIT Bombay
57 Module 4: Lecture 8 on Stress-strain relationship and Shear strength of soils Contents Stress state, Mohr s circle analysis and Pole, Principal stressspace, Stress pathsin p-q space; Mohr-Coulomb failure
More informationTHE PROXIMITY OF MICROVIAS TO PTHs AND ITS IMPACT ON THE RELIABILITY OF THESE MICROVIAS
THE PROXIMITY OF MICROVIAS TO PTHs AND ITS IMPACT ON THE RELIABILITY OF THESE MICROVIAS Jaydutt Joshi Principle Program Manager RF Module Technology Development Skyworks Solutions, Inc. Irvine, California
More informationSINTAP/TWI/ /46/99 SEPTEMBER 1999 STRESS INTENSITY DUE TO RESIDUAL STRESSES
SINTAP/TWI/4-6 88269/46/99 SEPTEMBE 999 STESS INTENSITY DUE TO ESIDUAL STESSES For: SINTAP SINTAP/TWI/4-6 TWI EPOT NO: 88269/46/99 SEPTEMBE 999 Prepared for: SINTAP Project, Task 4 Contact : Contact 2:
More informationCalculating Stress Intensity Factor (Mode I) for Composite Plate with Central Crack
Calculating Stress Intensity Factor (Mode I) for Composite Plate with Central Crack Luay S. Al-Ansari, Ph.D University of Kufa College of Engineering Mechanical Engineering Dept. Hashim N. Al-Mahmud University
More informationUNIVERSITY OF CALIFORNIA College of Engineering Departments of Mechanical Engineering and Material Science & Engineering
Fall 006 UNIVERSITY OF CALIFORNIA College of Engineering Departments of Mechanical Engineering and Material Science & Engineering MSEc113/MEc14 Mechanical Behavior of Materials Midterm #1 September 19
More informationGT RETROFIT AND RE-DESIGN OF A FOUR-STAGE CENTRIFUGAL COMPRESSOR
Proceedings of ASME Turbo Expo 2015: Turbine Technical Conference and Exposition GT2015 June 15-19, 2015, Montreal GT2015-43760 RETROFIT AND RE-DESIGN OF A FOUR-STAGE CENTRIFUGAL COMPRESSOR Chao Chen SoftInWay,
More informationDevelopment on Methods for Evaluating Structure Reliability of Piping Components
Transactions of the 17 th International Conference on Structural Mechanics in Reactor Technology (SMiRT 17) Prague, Czech Republic, August 17, 00 Paper # M0- Development on Methods for Evaluating Structure
More informationAn Improvement of Vegetation Height Estimation Using Multi-baseline Polarimetric Interferometric SAR Data
PIERS ONLINE, VOL. 5, NO. 1, 29 6 An Improvement of Vegetation Height Estimation Using Multi-baseline Polarimetric Interferometric SAR Data Y. S. Zhou 1,2,3, W. Hong 1,2, and F. Cao 1,2 1 National Key
More informationMECHANICS OF MATERIALS
CHAPTER 7 Transformations MECHANICS OF MATERIALS Ferdinand P. Beer E. Russell Johnston, Jr. John T. DeWolf Lecture Notes: J. Walt Oler Texas Tech Universit of Stress and Strain 006 The McGraw-Hill Companies,
More informationPlastic Failure of locally supported Silos with U-shaped longitudinal Stiffeners
Plastic Failure of locally supported Silos with U-shaped longitudinal Stiffeners *Arne Jansseune 1), Wouter De Corte 2) and Rudy Van Impe 3) 1), 2) Department of Civil Engineering, Associated Faculty of
More informationSAMPLE PROJECT IN LONDON DOCUMENT NO. STR-CALC UNITISED CURTAIN WALL 65 ENGINEER PROJECT. Pages REVISION TITLE
PROJECT ENGINEER DOCUMENT NO. STR-CALC-253 0 REVISION TITLE Pages UNITISED CURTAIN WALL 65 UNITISED CURTAIN WALL 2 of 65 Contents 1 Basic Data... 3 1.1 References... 3 1.2 Materials... 3 1.3 Loads... 4
More informationProject: IEEE P Working Group for Wireless Personal Area Networks (WPANs)
Project: IEEE P802.15 Working Group for Wireless Personal Area Networks (WPANs) Title: Channel Model for Intra-Device Communications Date Submitted: 15 January 2016 Source: Alexander Fricke, Thomas Kürner,
More informationConnecting Garage Door Jambs to Building Framing
Connecting Garage Door Jambs to Building Framing Introduction The members of DASMA recognize that connecting garage doors to building framing is as important as the design of garage doors themselves. The
More information1 SE = Student Edition - TG = Teacher s Guide
Mathematics State Goal 6: Number Sense Standard 6A Representations and Ordering Read, Write, and Represent Numbers 6.8.01 Read, write, and recognize equivalent representations of integer powers of 10.
More informationCommon Core Georgia Performance Standards
A Correlation of Pearson Connected Mathematics 2 2012 to the Common Core Georgia Performance s Grade 6 FORMAT FOR CORRELATION TO THE COMMON CORE GEORGIA PERFORMANCE STANDARDS (CCGPS) Subject Area: K-12
More informationModelling strategies for bivariate circular data
Modelling strategies for bivariate circular data John T. Kent*, Kanti V. Mardia, & Charles C. Taylor Department of Statistics, University of Leeds 1 Introduction On the torus there are two common approaches
More information#161. Connecting Garage Door Jambs to Building Framing. Introduction
Connecting Garage Door Jambs to Building Framing Introduction The members of DASMA are equally as concerned about connecting garage doors to building framing as they are about the design of garage doors
More informationStrike Point Control on EAST Using an Isoflux Control Method
Plasma Science and Technology, Vol.17, No.9, Sep. 2015 Strike Point Control on EAST Using an Isoflux Control Method XING Zhe ( ) 1, XIAO Bingjia ( ) 1,2, LUO Zhengping ( ) 1, M. L. WALKER 3, D. A. HUMPHREYS
More informationIndian Sovereign Yield Curve using Nelson-Siegel-Svensson Model
Indian Sovereign Yield Curve using Nelson-Siegel-Svensson Model Of the three methods of valuing a Fixed Income Security Current Yield, YTM and the Coupon, the most common method followed is the Yield To
More informationSpike Statistics: A Tutorial
Spike Statistics: A Tutorial File: spike statistics4.tex JV Stone, Psychology Department, Sheffield University, England. Email: j.v.stone@sheffield.ac.uk December 10, 2007 1 Introduction Why do we need
More informationSpike Statistics. File: spike statistics3.tex JV Stone Psychology Department, Sheffield University, England.
Spike Statistics File: spike statistics3.tex JV Stone Psychology Department, Sheffield University, England. Email: j.v.stone@sheffield.ac.uk November 27, 2007 1 Introduction Why do we need to know about
More informationLecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions
Lecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions ELE 525: Random Processes in Information Systems Hisashi Kobayashi Department of Electrical Engineering
More informationSTRUKTOL TR 451: Improvements When Using Mineral Fillers
STRUKTOL TR 451: Improvements When Using Mineral Fillers Summary and Data December 2012 STP0286 STRUKTOL TR 451 Functional viscosity modifier and in-situ mineral filler treatment for use in polyolefins
More informationSNELL S LAW AND UNIFORM REFRACTION. Contents
SNELL S LAW AND UNIFORM REFRACTION CRISTIAN E. GUTIÉRREZ Contents 1. Snell s law of refraction 1 1.1. In vector form 1 1.2. κ < 1 2 1.3. κ > 1 3 1.4. κ = 1 4 2. Uniform refraction 4 2.1. Surfaces with
More informationPSR-9000 FLX03G LDI (UL Name: PSR-9000AD / CA-90AD )
PSR-9000 FLX03G LDI (UL Name: PSR-9000AD / CA-90AD ) LIQUID PHOTOIMAGEABLE SOLDER MASK Designed for Flexible Printed Circuits on LDI equipment Screen Print Application Halogen-Free Compatible with Lead-Free
More informationAlternative VaR Models
Alternative VaR Models Neil Roeth, Senior Risk Developer, TFG Financial Systems. 15 th July 2015 Abstract We describe a variety of VaR models in terms of their key attributes and differences, e.g., parametric
More information2 Exploring Univariate Data
2 Exploring Univariate Data A good picture is worth more than a thousand words! Having the data collected we examine them to get a feel for they main messages and any surprising features, before attempting
More informationCrystal Lattices. Daniele Toffoli December 7, / 42
Crystal Lattices Daniele Toffoli December 7, 2016 1 / 42 Outline 1 Bravais Lattice: definitions and examples 2 Unit cell: primitive, conventional and Wigner-Seitz 3 Crystal structure: lattices with basis
More informationSome statistical properties of surface slopes via remote sensing considering a non-gaussian probability density function
Journal of odern Optics ISSN: 0950-0340 Print 36-3044 Online Journal homepage: http://www.tandfonline.com/loi/tmop0 Some statistical properties of surface slopes via remote sensing considering a non-gaussian
More informationEE/AA 578 Univ. of Washington, Fall Homework 8
EE/AA 578 Univ. of Washington, Fall 2016 Homework 8 1. Multi-label SVM. The basic Support Vector Machine (SVM) described in the lecture (and textbook) is used for classification of data with two labels.
More informationConnected Mathematics 2, 6 th and 7th Grade Units 2009 Correlated to: Washington Mathematics Standards (Grade 6)
Grade 6 6.1. Core Content: Multiplication and division of fractions and decimals (Numbers, Operations, Algebra) 6.1.A Compare and order non-negative fractions, decimals, and integers using the number line,
More informationRapid computation of prices and deltas of nth to default swaps in the Li Model
Rapid computation of prices and deltas of nth to default swaps in the Li Model Mark Joshi, Dherminder Kainth QUARC RBS Group Risk Management Summary Basic description of an nth to default swap Introduction
More information1 algebraic. expression. at least one operation. Any letter can be used as a variable. 2 + n. combination of numbers and variables
1 algebraic expression at least one operation 2 + n r w q Any letter can be used as a variable. combination of numbers and variables DEFINE: A group of numbers, symbols, and variables that represent an
More informationOnline Appendix A: Verification of Employer Responses
Online Appendix for: Do Employer Pension Contributions Reflect Employee Preferences? Evidence from a Retirement Savings Reform in Denmark, by Itzik Fadlon, Jessica Laird, and Torben Heien Nielsen Online
More informationA probability distribution shows the possible outcomes of an experiment and the probability of each of these outcomes.
Introduction In the previous chapter we discussed the basic concepts of probability and described how the rules of addition and multiplication were used to compute probabilities. In this chapter we expand
More informationFINITE ELEMENT-BASED FATIGUE LIFE PREDICTION OF A LOAD- CARRYING CRUCIFORM JOINT. Bandar Baru Bangi, Selangor, Malaysia
Journal of Mechanical Engineering and Sciences (JMES) ISSN (Print): 2289-4659; e-issn: 2231-8380; Volume 8, pp. 1414-1425, June 2015 Universiti Malaysia Pahang, Malaysia DOI: http://dx.doi.org/10.15282/jmes.8.2015.16.0138
More informationUncertainty Analysis with UNICORN
Uncertainty Analysis with UNICORN D.A.Ababei D.Kurowicka R.M.Cooke D.A.Ababei@ewi.tudelft.nl D.Kurowicka@ewi.tudelft.nl R.M.Cooke@ewi.tudelft.nl Delft Institute for Applied Mathematics Delft University
More informationInternational Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, ISSN
International Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, www.ijcea.com ISSN 31-3469 AN INVESTIGATION OF FINANCIAL TIME SERIES PREDICTION USING BACK PROPAGATION NEURAL
More informationPROBABILISTIC MICROMECHANICAL MODEL OF ENGINEERED CEMENTITIOUS COMPOSITES (ECC)
PROBABILISTIC MICROMECHANICAL MODEL OF ENGINEERED CEMENTITIOUS COMPOSITES (ECC) J. Li (1) and E.-H. Yang (2) (1) Interdisciplinary Graduate School, Nanyang Technological University, Singapore (2) School
More informationChapter 4 Random Variables & Probability. Chapter 4.5, 6, 8 Probability Distributions for Continuous Random Variables
Chapter 4.5, 6, 8 Probability for Continuous Random Variables Discrete vs. continuous random variables Examples of continuous distributions o Uniform o Exponential o Normal Recall: A random variable =
More informationTalent Attraction USING DATA TO IDENTIFY AND ADDRESS WORKFORCE CHALLENGES
Talent Attraction USING DATA TO IDENTIFY AND ADDRESS WORKFORCE CHALLENGES Topics The EDO dilemma Jobs and skills shortage Threats (and benefits) of automation Workforce trends Geography of jobs Emerging
More informationEvolution of Strategies with Different Representation Schemes. in a Spatial Iterated Prisoner s Dilemma Game
Submitted to IEEE Transactions on Computational Intelligence and AI in Games (Final) Evolution of Strategies with Different Representation Schemes in a Spatial Iterated Prisoner s Dilemma Game Hisao Ishibuchi,
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationModels of Patterns. Lecture 3, SMMD 2005 Bob Stine
Models of Patterns Lecture 3, SMMD 2005 Bob Stine Review Speculative investing and portfolios Risk and variance Volatility adjusted return Volatility drag Dependence Covariance Review Example Stock and
More informationInternational Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, ISSN
Volume XII, Issue II, Feb. 18, www.ijcea.com ISSN 31-3469 AN INVESTIGATION OF FINANCIAL TIME SERIES PREDICTION USING BACK PROPAGATION NEURAL NETWORKS K. Jayanthi, Dr. K. Suresh 1 Department of Computer
More informationAccelerated Option Pricing Multiple Scenarios
Accelerated Option Pricing in Multiple Scenarios 04.07.2008 Stefan Dirnstorfer (stefan@thetaris.com) Andreas J. Grau (grau@thetaris.com) 1 Abstract This paper covers a massive acceleration of Monte-Carlo
More informationLECTURE 2: MULTIPERIOD MODELS AND TREES
LECTURE 2: MULTIPERIOD MODELS AND TREES 1. Introduction One-period models, which were the subject of Lecture 1, are of limited usefulness in the pricing and hedging of derivative securities. In real-world
More informationTHE IMPACT OF TEMPORAL DECORRELATION OVER FOREST TERRAIN IN POLARIMETRIC SAR INTERFEROMETRY
THE IMPACT OF TEMPORAL DECORRELATION OER FOREST TERRAIN IN POLARIMETRIC SAR INTERFEROMETRY Seung-kuk Lee, Florian Kugler, Irena Hajnsek, Konstantinos P. Papathanassiou Microwave and Radar Institute, German
More informationDesigning with an Inverted-F PCB Antenna
Page 1 Application Note 5052 24 April 2008 Designing with an Inverted-F PCB Antenna For the EM250 and EM2 Platforms This document describes an Inverted-F PCB antenna designed by Ember for use with both
More informationThe Margins of Global Sourcing: Theory and Evidence from U.S. Firms by Pol Antràs, Teresa C. Fort and Felix Tintelnot
The Margins of Global Sourcing: Theory and Evidence from U.S. Firms by Pol Antràs, Teresa C. Fort and Felix Tintelnot Online Theory Appendix Not for Publication) Equilibrium in the Complements-Pareto Case
More informationUnit M2.2 (All About) Stress
Unit M. (All About) Stress Readings: CDL 4., 4.3, 4.4 16.001/00 -- Unified Engineering Department of Aeronautics and Astronautics Massachusetts Institute of Technology LEARNING OBJECTIVES FOR UNIT M. Through
More informationM249 Diagnostic Quiz
THE OPEN UNIVERSITY Faculty of Mathematics and Computing M249 Diagnostic Quiz Prepared by the Course Team [Press to begin] c 2005, 2006 The Open University Last Revision Date: May 19, 2006 Version 4.2
More informationMonitoring - revisited
Monitoring - revisited Anders Ringgaard Kristensen Slide 1 Outline Filtering techniques applied to monitoring of daily gain in slaughter pigs: Introduction Basic monitoring Shewart control charts DLM and
More informationRECURSIVE VALUATION AND SENTIMENTS
1 / 32 RECURSIVE VALUATION AND SENTIMENTS Lars Peter Hansen Bendheim Lectures, Princeton University 2 / 32 RECURSIVE VALUATION AND SENTIMENTS ABSTRACT Expectations and uncertainty about growth rates that
More informationDATA SUMMARIZATION AND VISUALIZATION
APPENDIX DATA SUMMARIZATION AND VISUALIZATION PART 1 SUMMARIZATION 1: BUILDING BLOCKS OF DATA ANALYSIS 294 PART 2 PART 3 PART 4 VISUALIZATION: GRAPHS AND TABLES FOR SUMMARIZING AND ORGANIZING DATA 296
More informationBY TYPE, ($ MILLIONS)... 5
CHAPTER ONE: INTRODUCTION... 1 STUDY GOALS AND OBJECTIVES... 1 REASONS FOR DOING THE STUDY... 1 CONTRIBUTION OF THE STUDY... 2 SCOPE AND FORMAT... 2 INFORMATION SOURCES... 2 ANALYST S CREDENTIALS... 2
More informationPrentice Hall Connected Mathematics 2, 7th Grade Units 2009 Correlated to: Minnesota K-12 Academic Standards in Mathematics, 9/2008 (Grade 7)
7.1.1.1 Know that every rational number can be written as the ratio of two integers or as a terminating or repeating decimal. Recognize that π is not rational, but that it can be approximated by rational
More informationDesign Optimization with GA Fitness Functions based on Total Lifecycle Value or Cost
Proceedings of the 7th WSEAS International Conference on Evolutionary Computing, Cavtat, Croatia, June 12-14, 26 (pp1-6) Design Optimization with GA Fitness Functions based on Total Lifecycle Value or
More informationREGIONAL WORKSHOP ON TRAFFIC FORECASTING AND ECONOMIC PLANNING
International Civil Aviation Organization 27/8/10 WORKING PAPER REGIONAL WORKSHOP ON TRAFFIC FORECASTING AND ECONOMIC PLANNING Cairo 2 to 4 November 2010 Agenda Item 3 a): Forecasting Methodology (Presented
More informationNonlinear Manifold Learning for Financial Markets Integration
Nonlinear Manifold Learning for Financial Markets Integration George Tzagkarakis 1 & Thomas Dionysopoulos 1,2 1 EONOS Investment Technologies, Paris (FR) 2 Dalton Strategic Partnership, London (UK) Nice,
More informationF19: Introduction to Monte Carlo simulations. Ebrahim Shayesteh
F19: Introduction to Monte Carlo simulations Ebrahim Shayesteh Introduction and repetition Agenda Monte Carlo methods: Background, Introduction, Motivation Example 1: Buffon s needle Simple Sampling Example
More informationImplementing Models in Quantitative Finance: Methods and Cases
Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1
More informationTrading Financial Market s Fractal behaviour
Trading Financial Market s Fractal behaviour by Solon Saoulis CEO DelfiX ltd. (delfix.co.uk) Introduction In 1975, the noted mathematician Benoit Mandelbrot coined the term fractal (fragment) to define
More informationSub-surface glacial structure over Nordaustlandet using multi-frequency Pol-InSAR
ub-surface glacial structure over Nordaustlandet using multi-frequency Pol-InAR Jayanti harma, Irena Hajnsek, Kostas Papathanassiou DLR (German Aerospace Center),.6.8 jayanti.sharma@dlr.de Introduction:
More informationWindow Width Selection for L 2 Adjusted Quantile Regression
Window Width Selection for L 2 Adjusted Quantile Regression Yoonsuh Jung, The Ohio State University Steven N. MacEachern, The Ohio State University Yoonkyung Lee, The Ohio State University Technical Report
More informationStock Market Prediction using Artificial Neural Networks IME611 - Financial Engineering Indian Institute of Technology, Kanpur (208016), India
Stock Market Prediction using Artificial Neural Networks IME611 - Financial Engineering Indian Institute of Technology, Kanpur (208016), India Name Pallav Ranka (13457) Abstract Investors in stock market
More informationInternational Trade Gravity Model
International Trade Gravity Model Yiqing Xie School of Economics Fudan University Dec. 20, 2013 Yiqing Xie (Fudan University) Int l Trade - Gravity (Chaney and HMR) Dec. 20, 2013 1 / 23 Outline Chaney
More informationNanotechnologie in der Anwendung was tut sich in den Unternehmen?
Nanotechnologie in der Anwendung was tut sich in den Unternehmen? Andreas Kornherr, Gerhard Drexler, Paul Achatz Nov 2009 Overview World s first Nano-hybrid laser copy paper New method to measure paper
More informationCurves on a Surface. (Com S 477/577 Notes) Yan-Bin Jia. Oct 24, 2017
Curves on a Surface (Com S 477/577 Notes) Yan-Bin Jia Oct 24, 2017 1 Normal and Geodesic Curvatures One way to examine how much a surface bends is to look at the curvature of curves on the surface. Let
More informationContinuous random variables
Continuous random variables probability density function (f(x)) the probability distribution function of a continuous random variable (analogous to the probability mass function for a discrete random variable),
More informationSOCIETY OF ACTUARIES Advanced Topics in General Insurance. Exam GIADV. Date: Thursday, May 1, 2014 Time: 2:00 p.m. 4:15 p.m.
SOCIETY OF ACTUARIES Exam GIADV Date: Thursday, May 1, 014 Time: :00 p.m. 4:15 p.m. INSTRUCTIONS TO CANDIDATES General Instructions 1. This examination has a total of 40 points. This exam consists of 8
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationDynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods
ISOPE 2010 Conference Beijing, China 24 June 2010 Dynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods Xi Ying Zhang, Zhi Ping Cheng, Jer-Fang Wu and Chee Chow Kei ABS 1 Main Contents
More informationProbability. An intro for calculus students P= Figure 1: A normal integral
Probability An intro for calculus students.8.6.4.2 P=.87 2 3 4 Figure : A normal integral Suppose we flip a coin 2 times; what is the probability that we get more than 2 heads? Suppose we roll a six-sided
More informationInvestigation of Interaction between Guidewire and Native Vessel Using Finite Element Analysis
Visit the SIMULIA Resource Center for more customer examples. Investigation of Interaction between Guidewire and Native Vessel Using Finite Element Analysis Atul Gupta 1, Subham Sett 2, Srinivasan Varahoor
More informationpplied instrumentation: ID Name Range Unit Division M900007 Calliper (0 250,00) mm 0,01 M900008 Pulloff tester ERICHSEN 417 0 až 47,00 MPa 0,5 M900009 Balance Kern PRJ 62002NM 0 až 6200 g 0,01 M900011
More informationThe rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx
1 Cumulants 1.1 Definition The rth moment of a real-valued random variable X with density f(x) is µ r = E(X r ) = x r f(x) dx for integer r = 0, 1,.... The value is assumed to be finite. Provided that
More informationME3620. Theory of Engineering Experimentation. Spring Chapter III. Random Variables and Probability Distributions.
ME3620 Theory of Engineering Experimentation Chapter III. Random Variables and Probability Distributions Chapter III 1 3.2 Random Variables In an experiment, a measurement is usually denoted by a variable
More informationResearch Article Skewness and Kurtosis: Important Parameters in the Characterization of Dental Implant Surface Roughness A Computer Simulation
International Scholarly Research Network ISRN Materials Science Volume, Article ID 353, 6 pages doi:.54//353 Research Article Skewness and Kurtosis: Important Parameters in the Characterization of Dental
More informationProject 1: Double Pendulum
Final Projects Introduction to Numerical Analysis II http://www.math.ucsb.edu/ atzberg/winter2009numericalanalysis/index.html Professor: Paul J. Atzberger Due: Friday, March 20th Turn in to TA s Mailbox:
More informationEffect of Payment Reduction on Default
B Effect of Payment Reduction on Default In this section we analyze the effect of payment reduction on borrower default. Using a regression discontinuity empirical strategy, we find that immediate payment
More information[D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright
Faculty and Institute of Actuaries Claims Reserving Manual v.2 (09/1997) Section D7 [D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright 1. Introduction
More informationCalibration and Parameter Risk Analysis for Gas Storage Models
Calibration and Parameter Risk Analysis for Gas Storage Models Greg Kiely (Gazprom) Mark Cummins (Dublin City University) Bernard Murphy (University of Limerick) New Abstract Model Risk Management: Regulatory
More informationAppendix A. Selecting and Using Probability Distributions. In this appendix
Appendix A Selecting and Using Probability Distributions In this appendix Understanding probability distributions Selecting a probability distribution Using basic distributions Using continuous distributions
More information