Computational Approach to Materials Science and Engineering
|
|
- Claud Bradley
- 5 years ago
- Views:
Transcription
1 Computational Approach to Materials Science and Engineering Prita Pant and M. P. Gururajan October, 2012 Copyright c 2012, Prita Pant and M P Gururajan. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.3 or any later version published by the Free Software Foundation; with no Invariant Sections, no Front-Cover Texts, and no Back- Cover Texts. A copy of the license is included in the section entitled GNU Free Documentation License. 1
2 Module: Interpolation 1 Pre-requisites The modules in Part II of this course material. 2 Learning goals Given a set of data, to interpolate between the given data points 3 Interpolation Interpolation is a procedure that is used to obtain the values of y for any x in the range x 1 to x 2 given the data points (x 1,y 1 ) and (x 2,y 2 ). There are many different interpolations that are possible. In this module, we shall discuss some examples. 4 Air bubbles in polar ice It is possible to track the atmospheric concentrations of CO 2 on a millennial timescale by analysing the air bubbles trapped inthepolar ice. Forexample, consider the data of the age of air bubbles found in two different ice cores, namely GRID in Greenland and Byrd in Antarctica, shown in Table. 1. The age is given in units of years Before Present (BP), that is, years calculated with 1 January 1950 as the origin (after 1950, the years can not be reliably calculated using radioactive isotope abundance data due to nuclear weapons testing). Given this data, let us consider the following questions: 1. To what depth in GRID ice core should one go to find air bubbles of age years BP? 2. At a depth of 1600 m in the Byrd ice core, what would be the age of air bubbles? 2
3 3. Adepth of 2000min GRID ice core corresponds to what depth in Byrd ice core? All these questions can be answered using interpolation. 4.1 Linear interpolation The simplest interpolation to consider is the linear interpolation. In this, one assumes that the data at the two points (x 1,y 1 ) and (x 2,y 2 ) can be connected through a straight line so that for any intermediate value of x, the corresponding y value can be obtained and vice versa. The GNU Octave command interp1 can be used to carry out such a linear interpolation. In the script below, for example, we show how interp1 can be used to answer all the three questions above: and the answers for the three questions are, m, years BP and m, repsectively. DepthGasAge.oct Copyright (C) 2011 Prita Pant and M P Gururajan This program is free software; you can redistribute it and/or modify it under the terms of GNU General Public License as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA , USA % Load the given data from the data file X = load("depthgasage.dat"); 3
4 Age (in years BP) Depth (in m) Depth (in m) GRID: Greenland Byrd: Antarctica Table 1: The age of air bubbles at various depths in two ice cores, namely GRID in Greenland and Byrd in Antarctica. The data is taken from Atmospheric CO 2 concentration and millennial-scale 4 climate change during the last glacial period, B. Stauffer et al, Nature, 392, pp , 1954.
5 % Read the first column of the data as age A = X(:,1); % Read the second column of the data as depth for GRID dgrid = X(:,2); % Read the third column of the data as depth for Byrd dbyrd = X(:,3); % Interpolate the GRID data for the age of air bubble to be interp1(a,dgrid,20500) % Interpolate the Byrd data for the depth of 1600 interp1(dbyrd,a,1600) % Interpolate the Byrd data to a depth of 2000 m of GRID interp1(dgrid,dbyrd,2000) From the above script, it is clear that the x and y data are given as the first two parameters of interp1 and the third parameter is the x for which we want to calculate the y value. For example, a command such as interp1(dbyrd, dgrid, 1880) will give the depth of GRID that correpsonds to 1880 m of Byrd while interp1(dgrid,dbyrd,1880) gave the depth of Byrd that corresponds to 1880 m of GRID. Here are two important points to remember while using interp1. The sample points x, that is, the first set of parameters passed on to interp1 should be strictly monotonic that is, either continuously increasing or decreasing. The point at which we want to calculate the y value should lie between the given data; calculating the y values for x values that lie outside the given range of sample data is known as extrapolation. It is possible to exptrapolate using interp1. For that, we need to explicitly mention extrap as the fourth input parameter in the command: for example, interp1(dbyrd,dgrid,1900, extrap ) will give answers while without the extrap string, the program will return NA. 5 Error function evaluation Consider a diffusion couple, consisting of two very long rods, welded face to face as shown. The composition in such a diffusion couple, as a fun- 5
6 Diffusion couple c c 1 2 c t=0 Figure 1: The schematic of a diffusion couple. Two long rods of different compositions, c 1 and c 2 ae welded face to face. The initial composition profile looks as shown below the couple. If this assembly is kept at high enough temperatures, then the diffusion profile evovles and the profile is described by error function. x tion of time, is known to be described by the expression c(x,t) = c 1+c 2 2 c 1 c 2 erf ( ) x 2 2 Dt where c1 and c 2 are the initial compositions of the two rods, and D is the diffusivity. Theerrorfunction, erf(z), foranyz hastoreadofffromtables. Typically, the tabulation is done for a certain number of z values and for any intermediate values, the error function has to be obtained using interpolation. Let us suppose that we know the error function values for nine different points as shown in Table. 2. For any intermediate values of z, one can then calculate the error function values by interpolation. 6
7 z erf(z) Table 2: Tabulation of error function values in the range x = 2 to x = 2. See Materials Science and Engineering: a first course, V Raghavan, Third edition, Prentice-Hall of India Pvt. Ltd, 1995, for example. With the given data, if we try a linear interpolation (See Problem 1), the interpolation looks as shown in Figure. 2. As is clear from the figure, the interpolation is not satisfactory since it gives rise to a jagged profile. It indeed is possible to get a much better plot of error function using a table with much finer sampling: see Table. 3 and Fig Spline interpolation It is also possible to get a smooth curve for the error function by carrying out spline interpolation instead of linear interpolation. In the script below we show how to fit even a small number of data (as given in Table. 2) using spline interpolation to achieve a smoother profile, which is comparable to what is achieved using a much larger data set. The plot obtained using the spline interpolation is shown in Fig. 4 7
8 1 0.5 erf(z) z Figure 2: The linear interpolation for the error function. 8
9 z erf(z) z erf(z) Table3: Tabulationoferrorfunctionvaluesintherangex = 2.8tox = 2.8. See Materials Science and Engineering: 9 a first course, V Raghavan, Third edition, Prentice-Hall of India Pvt. Ltd, 1995, for example.
10 1 0.5 erf(z) z Figure 3: The plot of error function with a much finer tabulation. 10
11 ErfSplineInterpolation.oct Copyright (C) 2011 Prita Pant and M P Gururajan This program is free software; you can redistribute it and/or modify it under the terms of GNU General Public License as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA , USA X = load( errorfunction.dat ); a = X(:,1); b = X(:,2); plot(a,b, o ); hold on x = [-2.:0.1:2.]; xx = interp1(a,b,x, spline ); plot(x,xx, r ) axis( square ) xlabel("z") ylabel("erf(z)") print -depsc../figures/erfsplineinterpolation.eps 11
12 1 0.5 erf(z) z Figure 4: The plot of error function with very few points of tabulated data and spline interpolation. The spline interpolation results in a smoother profile. 12
13 6 The specific heat data: polynomial interpolation In the previous section, we have used either polyfit or the design matrix to fit the given data to polynomials. It is also possible to use the fitted polynomials for interpolation. Thus the command polyval can be used for interpolation if polyfit was used for fitting data; in the case a design matrix wasused, thecoefficients canbeusedforinterpolation. Seethemoduleonfitting for further examples of polynomial interpolation, with specific reference to specific heats. 7 Self-assessment questions 1. What is the command in GNU Octave for interpolation? 2. interp1(x,y,a) means the x value corresponding to y = a is obtained. True or false? 3. polyval can be used for interpolation. True or false? 4. Can interp1 be used for exptrapolation? 5. In interp1(x,y,a), should the y also be monotonic? 8 Answers to self-assessment questions 1. interp1 2. False. It means the y value corresponding to x = a is obtained. 3. True. If polyfit is used on the given data to fit a polynomial, then, polyval can be used for interpolation. 4. Yes; by passing an optional parameter extrap, extrapolation can be carried out. 5. No. See the Problem 2, for example. 13
14 x y Table 4: Tabulation of x and y values in which only x values are monotonic. The y values are oscillatory. 9 Exercises Problem 1. Use the data given in Table. 2 to do linear interpolation for z values in the range z = 2 to z = 2 (in increments of 0.1) and calculate erf(z). Problem 2. Do a linear and spline interpolation for the data given in Table. 4. Problem 3. The stress-strain data for a brass sample is as shown in Table. 5. Calculate the strain at a stress of 50 MPa and the stress for a strain of Solution to exercises Solution to Problem 1. 14
15 Stress (in MPa) Strain Table 5: Tabulation of stress strain data for a brass sample. X = load( errorfunction.dat ); a = X(:,1); b = X(:,2); plot(a,b, o ); hold on x = [-2.:0.1:2.]; xx = interp1(a,b,x); plot(x,xx, r ) axis( square ) xlabel("z") ylabel("erf(z)") print -depsc../figures/erflinearinterpolation.eps The plot generated by this script is shown in Fig. 2. Solution to Problem 2 The script for carrying out linear interpolation on the given data is given below and the figure generated is shown in Fig
16 y x Figure 5: Linear interpolation for the data given in Problem 2. Note that even if the y values are non-monotonic, the interpolation can be carried out as long as x is monotonic. From the figure, it is clear that the given data is for a sinusoidal curve. X = load( SineData.dat ); a = X(:,1); b = X(:,2); plot(a,b, o ); hold on x = [0.0:0.01:6.3]; xx = interp1(a,b,x); plot(x,xx, r ) axis( square ) xlabel("x") ylabel("y") print -depsc../figures/sinelinear.eps 16
17 The script for carrying out linear interpolation on the given data is given below and the figure generated is shown in Fig. 6. X = load( SineData.dat ); a = X(:,1); b = X(:,2); plot(a,b, o ); hold on x = [0.0:0.01:6.3]; xx = interp1(a,b,x, spline ); plot(x,xx, r ) axis( square ) xlabel("x") ylabel("y") print -depsc../figures/sinespline.eps Note that even though the spline fitting gives the sinusoidal curve much better than the linear fit, it is not exact. This can be clearly seen if we compare the calculated value from spline and compare it with the actual value of sine at the given point. Solution to Problem 3 By plotting the data, one can see that the stress strain plot is linear; so, we are within the elastic limit. Hence, for finding the strain, one can use linear interpolation. Using linear interpolation, we get a strain of for a stress of 50 MPa and a stress of MPa for a strain of The code for calculation is as shown below. X = load( StressStrain.dat ); a = X(:,1); b = X(:,2); interp1(a,b,50) interp1(b,a, ) 17
18 y x Figure 6: Spline interpolation for the data given in Problem 2. Note that even if the y values are non-monotonic, the interpolation can be carried out as long as x is monotonic. From the figure, it is clear that the given data is for a sinusoidal curve. 18
19 11 References and further reading 1. Advanced Engineering Mathematics, E. Kreyszig, 8th edition, John Wiley and Sons, Elementary Numerical Analysis, K. E. Atkinson, 3rd edition, Wiley India,
Exercise(s) Solution(s) to the exercise(s)
Exercise(s) Problem 1. Counting configurations Consider two different types of atoms, say A and B (represented by red and blue, respectively in the figures). Let A atoms and B atoms be distributed on N
More informationComputational Approach to Materials Science and Engineering
Computational Approach to Materials Science and Engineering Prita Pant and M. P. Gururajan December, 2012 Copyright c 2012, Prita Pant and M P Gururajan. Permission is granted to copy, distribute and/or
More informationComputational Approach to Materials Science and Engineering
Computational Approach to Materials Science and Engineering Prita Pant and M. P. Gururajan January, 2012 Copyright c 2011, Prita Pant and M P Gururajan. Permission is granted to copy, distribute and/or
More informationAn introduction to interpolation and splines
An introduction to interpolation and splines Kenneth H. Carpenter, EECE KSU November 22, 1999 revised November 20, 2001, April 24, 2002, April 14, 2004 1 Introduction Suppose one wishes to draw a curve
More informationChapter 12: Quadratic and Cubic Graphs
Chapter 12: Quadratic and Cubic Graphs Section 12.1 Quadratic Graphs x 2 + 2 a 2 + 2a - 6 r r 2 x 2 5x + 8 2y 2 + 9y + 2 All the above equations contain a squared number. They are therefore called quadratic
More informationComputational Approach to Materials Science and Engineering
Computational Approach to Materials Science and Engineering Prita Pant and M. P. Gururajan January, 2012 Copyright c 2012, Prita Pant and M P Gururajan. Permission is granted to copy, distribute and/or
More informationChapter 19 Interpolation
19.1 One-Dimensional Interpolation Chapter 19 Interpolation Empirical data obtained experimentally often times conforms to a fixed (deterministic) but unkown functional relationship. When estimates of
More informationInternational Journal of Emerging Technologies in Computational and Applied Sciences (IJETCAS)
International Association of Scientific Innovation and Research (IASIR) (An Association Unifying the Sciences, Engineering, and Applied Research) International Journal of Emerging Technologies in Computational
More information99 International Journal of Engineering, Science and Mathematics
Journal Homepage: Applications of cubic splines in the numerical solution of polynomials Najmuddin Ahmad 1 and Khan Farah Deeba 2 Department of Mathematics Integral University Lucknow Abstract: In this
More informationNumber Song Names by Play Order v2.0
Number Song Names by Play Order v2.0 AppleScript for itunes Find more free AppleScripts and info on writing your own at Doug's AppleScripts for itunes. This script will prefix each selected track's Song
More informationSupplemental Material Deep Fluids: A Generative Network for Parameterized Fluid Simulations
Supplemental Material Deep Fluids: A Generative Network for Parameterized Fluid Simulations 1. Extended Results 1.1. 2-D Smoke Plume Additional results for the 2-D smoke plume example are shown in Figures
More informationLecture 8. Divided Differences,Least-Squares Approximations. Ceng375 Numerical Computations at December 9, 2010
Lecture 8, Ceng375 Numerical Computations at December 9, 2010 Computer Engineering Department Çankaya University 8.1 Contents 1 2 3 8.2 : These provide a more efficient way to construct an interpolating
More informationAPPLICATIONS OF INTERPOLATION Y.Lavanya 1, Ch.Achireddy 2, M. Sudheer kumar 3
APPLICATIONS OF INTERPOLATION Y.Lavanya 1, Ch.Achireddy 2, M. Sudheer kumar 3 1 Assistant professor, Department Mathematics, MLRIT, Hyderabad 2 Associate professor(phd), Department Mathematics MLRIT, Hyderabad
More informationVST Preset Generator Documentation. François Mazen V0.2.8
VST Preset Generator Documentation François Mazen V0.2.8 Table of Contents Installation................................................................................. 1 Overview..................................................................................
More informationCOURSE: NUMERICAL ANALYSIS. LESSON: Methods for Solving Non-Linear Equations
COURSE: NUMERICAL ANALYSIS LESSON: Methods for Solving Non-Linear Equations Lesson Developer: RAJNI ARORA COLLEGE/DEPARTMENT: Department of Mathematics, University of Delhi Page No. 1 Contents 1. LEARNING
More informationMATLAB Examples. Interpolation and Curve Fitting. Hans-Petter Halvorsen
MATLAB Examples Interpolation and Curve Fitting Hans-Petter Halvorsen Interpolation Interpolation is used to estimate data points between two known points. The most common interpolation technique is Linear
More informationFriday, 11 January 13. Interpolation
Interpolation Interpolation Interpolation is not a branch of mathematic but a collection of techniques useful for solving computer graphics problems Basically an interpolant is a way of changing one number
More informationThe Complexity of Algorithms (3A) Young Won Lim 4/3/18
Copyright (c) 2015-2018 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published
More informationSelected Tags to Lyrics v1.1
Selected Tags to Lyrics v1.1 AppleScript for itunes Find more free AppleScripts and info on writing your own at Doug's AppleScripts for itunes. Primarily for iphone and ipod Touch users, this script will
More informationFitting Uncertain Data with NURBS
Fitting Uncertain Data with NURBS Wolfgang Heidrich, Richard Bartels, George Labahn Abstract. Fitting of uncertain data, that is, fitting of data points that are subject to some error, has important applications
More informationGChemTable manual. GChemTable manual
GChemTable manual i GChemTable manual GChemTable manual ii Copyright 2006-2012 Jean Bréfort Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation
More informationGChemCalc manual. GChemCalc manual
GChemCalc manual i GChemCalc manual GChemCalc manual ii Copyright 2006-2012 Jean Bréfort Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation
More informationExamples, examples: Outline
Examples, examples: Outline Overview of todays exercises Basic scripting Importing data Working with temporal data Working with missing data Interpolation in 1D Some time series analysis Linear regression
More informationDATA FITTING IN SCILAB
powered by DATA FITTING IN SCILAB In this tutorial the reader can learn about data fitting, interpolation and approximation in Scilab. Interpolation is very important in industrial applications for data
More informationExercise Set Decide whether each matrix below is an elementary matrix. (a) (b) (c) (d) Answer:
Understand the relationships between statements that are equivalent to the invertibility of a square matrix (Theorem 1.5.3). Use the inversion algorithm to find the inverse of an invertible matrix. Express
More informationdns-grind User Documentation
dns-grind User Documentation pentestmonkey@pentestmonkey.net 21 January 2007 Contents 1 Overview 2 2 Installation 2 3 Usage 3 4 Some Examples 3 4.1 Bruteforcing Hostnames (A-record Lookups)............
More informationSETTLEMENT OF A CIRCULAR FOOTING ON SAND
1 SETTLEMENT OF A CIRCULAR FOOTING ON SAND In this chapter a first application is considered, namely the settlement of a circular foundation footing on sand. This is the first step in becoming familiar
More informationDt100rc User Guide. Table of Contents. Prepared By: Peter Milne Date: 20 June 2005
Dt100rc User Guide Prepared By: Peter Milne Date: 20 June 2005 Rev Date Description 1 040604 First issue 2 050620 Updated. 3 060127 Table of Contents 1 Introduction...3 1.1 Features...3 1.2 References...3
More informationIntroduction to Matlab
Technische Universität München WT 21/11 Institut für Informatik Prof Dr H-J Bungartz Dipl-Tech Math S Schraufstetter Benjamin Peherstorfer, MSc October 22nd, 21 Introduction to Matlab Engineering Informatics
More informationOpen source licensing notices in Web applications
Open source licensing notices in Web applications 81 Open source licensing notices in Web applications Arnoud Engelfriet a (a) Associate, ICTRecht Legal Services. DOI: 10.5033/ifosslr.v3i1.47 Abstract
More informationOptimization Methods: Optimization using Calculus Kuhn-Tucker Conditions 1. Module - 2 Lecture Notes 5. Kuhn-Tucker Conditions
Optimization Methods: Optimization using Calculus Kuhn-Tucker Conditions Module - Lecture Notes 5 Kuhn-Tucker Conditions Introduction In the previous lecture the optimization of functions of multiple variables
More informationFinite Element Analysis of Ellipsoidal Head Pressure Vessel
Finite Element Analysis of Ellipsoidal Head Pressure Vessel Vikram V. Mane*, Vinayak H.Khatawate.**Ashok Patole*** * (Faculty; Mechanical Engineering Department, Vidyavardhini s college of Engineering.
More informationSet No. 1 IV B.Tech. I Semester Regular Examinations, November 2010 FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 Answer any FIVE Questions All Questions carry equal marks
More informationa 2 + 2a - 6 r r 2 To draw quadratic graphs, we shall be using the method we used for drawing the straight line graphs.
Chapter 12: Section 12.1 Quadratic Graphs x 2 + 2 a 2 + 2a - 6 r r 2 x 2 5x + 8 2 2 + 9 + 2 All the above equations contain a squared number. The are therefore called quadratic expressions or quadratic
More informationRational Numbers and the Coordinate Plane
Rational Numbers and the Coordinate Plane LAUNCH (8 MIN) Before How can you use the numbers placed on the grid to figure out the scale that is used? Can you tell what the signs of the x- and y-coordinates
More informationDr Richard Greenaway
SCHOOL OF PHYSICS, ASTRONOMY & MATHEMATICS 4PAM1008 MATLAB 1 Introduction to MATLAB Dr Richard Greenaway 1 Introduction to MATLAB 1.1 What is MATLAB? MATLAB is a high-level technical computing language
More informationFunctions (4A) Young Won Lim 5/8/17
Functions (4A) Copyright (c) 2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version
More informationTruss structural configuration optimization using the linear extended interior penalty function method
ANZIAM J. 46 (E) pp.c1311 C1326, 2006 C1311 Truss structural configuration optimization using the linear extended interior penalty function method Wahyu Kuntjoro Jamaluddin Mahmud (Received 25 October
More information8 Piecewise Polynomial Interpolation
Applied Math Notes by R. J. LeVeque 8 Piecewise Polynomial Interpolation 8. Pitfalls of high order interpolation Suppose we know the value of a function at several points on an interval and we wish to
More informationLecture 9. Curve fitting. Interpolation. Lecture in Numerical Methods from 28. April 2015 UVT. Lecture 9. Numerical. Interpolation his o
Curve fitting. Lecture in Methods from 28. April 2015 to ity Interpolation FIGURE A S Splines Piecewise relat UVT Agenda of today s lecture 1 Interpolation Idea 2 3 4 5 6 Splines Piecewise Interpolation
More informationSolar Radiation Data Modeling with a Novel Surface Fitting Approach
Solar Radiation Data Modeling with a Novel Surface Fitting Approach F. Onur Hocao glu, Ömer Nezih Gerek, Mehmet Kurban Anadolu University, Dept. of Electrical and Electronics Eng., Eskisehir, Turkey {fohocaoglu,ongerek,mkurban}
More informationfinger-user-enum User Documentation
finger-user-enum User Documentation pentestmonkey@pentestmonkey.net 21 January 2007 Contents 1 Overview 2 2 Installation 2 3 Usage 3 4 Some Examples 3 4.1 Normal Usage............................. 4 4.2
More informationA technique for constructing monotonic regression splines to enable non-linear transformation of GIS rasters
18 th World IMACS / MODSIM Congress, Cairns, Australia 13-17 July 2009 http://mssanz.org.au/modsim09 A technique for constructing monotonic regression splines to enable non-linear transformation of GIS
More informationAlgorithms Bubble Sort (1B) Young Won Lim 4/5/18
Algorithms Bubble Sort (1B) Young Won Lim 4/5/18 Copyright (c) 2017 2018 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation
More informationNonparametric regression using kernel and spline methods
Nonparametric regression using kernel and spline methods Jean D. Opsomer F. Jay Breidt March 3, 016 1 The statistical model When applying nonparametric regression methods, the researcher is interested
More informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More informationAn Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (2000/2001)
An Investigation into Iterative Methods for Solving Elliptic PDE s Andrew M Brown Computer Science/Maths Session (000/001) Summary The objectives of this project were as follows: 1) Investigate iterative
More informationcourse outline basic principles of numerical analysis, intro FEM
idealization, equilibrium, solutions, interpretation of results types of numerical engineering problems continuous vs discrete systems direct stiffness approach differential & variational formulation introduction
More informationFunctions (4A) Young Won Lim 3/16/18
Functions (4A) Copyright (c) 2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version
More informationNon Linear Calibration Curve and Polynomial
Non Linear Calibration Curve and Polynomial by Dr. Colin Mercer, Technical Director, Prosig Not all systems vary linearly. The issue discussed here is determining a non linear calibration curve and, if
More information15.10 Curve Interpolation using Uniform Cubic B-Spline Curves. CS Dept, UK
1 An analysis of the problem: To get the curve constructed, how many knots are needed? Consider the following case: So, to interpolate (n +1) data points, one needs (n +7) knots,, for a uniform cubic B-spline
More informationA Curve-Fitting Cookbook for use with the NMM Toolbox
A Curve-Fitting Cookbook for use with the NMM Toolbox Gerald Recktenwald October 17, 2000 Abstract Computational steps for obtaining curve fits with Matlab are described. The steps include reading data
More informationApplied Statistics : Practical 9
Applied Statistics : Practical 9 This practical explores nonparametric regression and shows how to fit a simple additive model. The first item introduces the necessary R commands for nonparametric regression
More informationDoubly Cyclic Smoothing Splines and Analysis of Seasonal Daily Pattern of CO2 Concentration in Antarctica
Boston-Keio Workshop 2016. Doubly Cyclic Smoothing Splines and Analysis of Seasonal Daily Pattern of CO2 Concentration in Antarctica... Mihoko Minami Keio University, Japan August 15, 2016 Joint work with
More informationInterpolation and curve fitting
CITS2401 Computer Analysis and Visualization School of Computer Science and Software Engineering Lecture 9 Interpolation and curve fitting 1 Summary Interpolation Curve fitting Linear regression (for single
More informationEngineering Analysis ENG 3420 Fall Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00
Engineering Analysis ENG 3420 Fall 2009 Dan C. Marinescu Office: HEC 439 B Office hours: Tu-Th 11:00-12:00 1 Lecture 24 Attention: The last homework HW5 and the last project are due on Tuesday November
More informationB-Spline Polynomials. B-Spline Polynomials. Uniform Cubic B-Spline Curves CS 460. Computer Graphics
CS 460 B-Spline Polynomials Computer Graphics Professor Richard Eckert March 24, 2004 B-Spline Polynomials Want local control Smoother curves B-spline curves: Segmented approximating curve 4 control points
More informationAnalysis and Design of Cantilever Springs
Analysis and Design of Cantilever Springs Hemendra Singh Shekhawat, Hong Zhou Department of Mechanical Engineering Texas A&M University-Kingsville Kingsville, Texas, USA Abstract Cantilever springs are
More informationTOPIC 6 Computer application for drawing 2D Graph
YOGYAKARTA STATE UNIVERSITY MATHEMATICS AND NATURAL SCIENCES FACULTY MATHEMATICS EDUCATION STUDY PROGRAM TOPIC 6 Computer application for drawing 2D Graph Plotting Elementary Functions Suppose we wish
More informationNUMERICAL INTEGRATION
NUMERICAL INTEGRATION f(x) MISN-0-349 NUMERICAL INTEGRATION by Robert Ehrlich George Mason University 1. Numerical Integration Algorithms a. Introduction.............................................1 b.
More informationLecture 4.5: Interpolation and Splines
Lecture 4.5: Interpolation and Splines D. Jason Koskinen koskinen@nbi.ku.dk Photo by Howard Jackman University of Copenhagen Advanced Methods in Applied Statistics Feb - Apr 2018 Niels Bohr Institute 2
More informationIn this course we will need a set of techniques to represent curves and surfaces in 2-d and 3-d. Some reasons for this include
Parametric Curves and Surfaces In this course we will need a set of techniques to represent curves and surfaces in 2-d and 3-d. Some reasons for this include Describing curves in space that objects move
More informationMake Video PDF Booklet v1.0
Make Video PDF Booklet v1.0 AppleScript for itunes Find more free AppleScripts and info on writing your own at Doug's AppleScripts for itunes. This script will create a PDF booklet containing video-oriented
More informationInternational Association of Scientific Innovation and Research (IASIR) (An Association Unifying the Sciences, Engineering, and Applied Research)
International Association of Scientific Innovation and Research (IASIR) (An Association Unifying the Sciences, Engineering, and Applied Research) International Journal of Emerging Technologies in Computational
More informationMath.1330 Section 5.2 Graphs of the Sine and Cosine Functions
Math.10 Section 5. Graphs of the Sine and Cosine Functions In this section, we will graph the basic sine function and the basic cosine function and then graph other sine and cosine functions using transformations.
More informationStudy on Shaking Table Test and Simulation Analysis of Graphite Dowel-Socket Structure
Study on Shaking Table Test and Simulation Analysis of Graphite Dowel-Socket Structure Xiangxiong Kong, Tiehua Shi & Shaoge Cheng Institute of Earthquake Engineering, China Academy of Building Research,
More informationExperiment # 5. Introduction to Error Control Codes
ECE 417 Winter 2003 Experiment # 5 Introduction to Error Control Codes 1 Purpose The purpose for this experiment is to provide you with an introduction to the field of error control coding. This will be
More informationChapter 1 Introduction
Chapter 1 Introduction GTU Paper Analysis (New Syllabus) Sr. No. Questions 26/10/16 11/05/16 09/05/16 08/12/15 Theory 1. What is graphic standard? Explain different CAD standards. 2. Write Bresenham s
More informationELECTRONOTES APPLICATION NOTE NO. 372
ELECTRONOTES APPLICATION NOTE NO. 372 116 Hanshaw Rd Ithaca, NY 1485 January 28 INTERPOLATION ERROR IN SINEWAVE TABLES INTRODUCTION: The basic notion of interpolation is fundamental to many things we do
More informationCHAPTER 6 Parametric Spline Curves
CHAPTER 6 Parametric Spline Curves When we introduced splines in Chapter 1 we focused on spline curves, or more precisely, vector valued spline functions. In Chapters 2 and 4 we then established the basic
More informationLecture VII : Random systems and random walk
Lecture VII : Random systems and random walk I. RANDOM PROCESSES In nature, no processes are truly deterministic. However, while dealing with many physical processes such as calculating trajectories of
More informationTHE BRAID INDEX OF ALTERNATING LINKS
THE BRAID INDEX OF ALTERNATING LINKS YUANAN DIAO, GÁBOR HETYEI AND PENGYU LIU Abstract. It is well known that the minimum crossing number of an alternating link equals the number of crossings in any reduced
More informationAgilent ChemStation for UV-visible Spectroscopy
Agilent ChemStation for UV-visible Spectroscopy Understanding Your Biochemical Analysis Software Agilent Technologies Notices Agilent Technologies, Inc. 2000, 2003-2008 No part of this manual may be reproduced
More information1. COURSE TITLE Multimedia Signal Processing I: visual signals Course number Course area Course type Course level. 1.5.
1. COURSE TITLE Multimedia Signal Processing I: visual signals 1.1. Course number 18768 1.2. Course area Computer Science Engineering 1.3. Course type Elective course 1.4. Course level Graduate 1.5. Year
More informationRecent Developments in Isogeometric Analysis with Solid Elements in LS-DYNA
Recent Developments in Isogeometric Analysis with Solid Elements in LS-DYNA Liping Li David Benson Attila Nagy Livermore Software Technology Corporation, Livermore, CA, USA Mattia Montanari Nik Petrinic
More informationSimulation of rotation and scaling algorithm for numerically modelled structures
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Simulation of rotation and scaling algorithm for numerically modelled structures To cite this article: S K Ruhit et al 2018 IOP
More informationIntermediate Algebra. Gregg Waterman Oregon Institute of Technology
Intermediate Algebra Gregg Waterman Oregon Institute of Technology c 2017 Gregg Waterman This work is licensed under the Creative Commons Attribution 4.0 International license. The essence of the license
More informationFunctions. Edexcel GCE. Core Mathematics C3
Edexcel GCE Core Mathematics C Functions Materials required for examination Mathematical Formulae (Green) Items included with question papers Nil Advice to Candidates You must ensure that your answers
More informationSketching graphs of polynomials
Sketching graphs of polynomials We want to draw the graphs of polynomial functions y = f(x). The degree of a polynomial in one variable x is the highest power of x that remains after terms have been collected.
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
More informationNumerical Methods 5633
Numerical Methods 5633 Lecture 3 Marina Krstic Marinkovic mmarina@maths.tcd.ie School of Mathematics Trinity College Dublin Marina Krstic Marinkovic 1 / 15 5633-Numerical Methods Organisational Assignment
More informationMATLAB Modul 3. Introduction
MATLAB Modul 3 Introduction to Computational Science: Modeling and Simulation for the Sciences, 2 nd Edition Angela B. Shiflet and George W. Shiflet Wofford College 2014 by Princeton University Press Introduction
More informationInterpolation and Basis Fns
CS148: Introduction to Computer Graphics and Imaging Interpolation and Basis Fns Topics Today Interpolation Linear and bilinear interpolation Barycentric interpolation Basis functions Square, triangle,,
More informationBasically, a graph is a representation of the relationship between two or more variables.
1 Drawing Graphs Introduction In recent years, the CSEC Integrated Science Examination, Paper 02, deals with graphical analysis. That is, data is presented in a tabular format and the student is asked
More informationFitting to a set of data. Lecture on fitting
Fitting to a set of data Lecture on fitting Linear regression Linear regression Residual is the amount difference between a real data point and a modeled data point Fitting a polynomial to data Could use
More informationCONTACT STATE AND STRESS ANALYSIS IN A KEY JOINT BY FEM
PERJODICA POLYTECHNICA SER. ME CH. ENG. VOL. 36, NO. 1, PP. -15-60 (1992) CONTACT STATE AND STRESS ANALYSIS IN A KEY JOINT BY FEM K. VARADI and D. M. VERGHESE Institute of Machine Design Technical University,
More informationSTANDARDS OF LEARNING CONTENT REVIEW NOTES ALGEBRA II. 3 rd Nine Weeks,
STANDARDS OF LEARNING CONTENT REVIEW NOTES ALGEBRA II 3 rd Nine Weeks, 2016-2017 1 OVERVIEW Algebra II Content Review Notes are designed by the High School Mathematics Steering Committee as a resource
More informationLS-DYNA CONTACT PROCEDURE ANALYSIS FOR SELECTED MECHANICAL SYSTEMS
Journal of KONES Powertrain and Transport, Vol. 22, No. 1 2015 LS-DYNA CONTACT PROCEDURE ANALYSIS FOR SELECTED MECHANICAL SYSTEMS Jerzy Małachowski, Jakub Bukala, Krzysztof Damaziak, Michał Tomaszewski
More information4.12 Generalization. In back-propagation learning, as many training examples as possible are typically used.
1 4.12 Generalization In back-propagation learning, as many training examples as possible are typically used. It is hoped that the network so designed generalizes well. A network generalizes well when
More informationMar. 20 Math 2335 sec 001 Spring 2014
Mar. 20 Math 2335 sec 001 Spring 2014 Chebyshev Polynomials Definition: For an integer n 0 define the function ( ) T n (x) = cos n cos 1 (x), 1 x 1. It can be shown that T n is a polynomial of degree n.
More informationDirect and Partial Variation. Lesson 12
Direct and Partial Variation Lesson MFMP Foundations of Mathematics Unit Lesson Lesson Twelve Concepts Overall Expectations Apply data-management techniques to investigate relationships between two variables;
More information[1] CURVE FITTING WITH EXCEL
1 Lecture 04 February 9, 2010 Tuesday Today is our third Excel lecture. Our two central themes are: (1) curve-fitting, and (2) linear algebra (matrices). We will have a 4 th lecture on Excel to further
More information4 Visualization and. Approximation
4 Visualization and Approximation b A slope field for the differential equation y tan(x + y) tan(x) tan(y). It is not always possible to write down an explicit formula for the solution to a differential
More informationCHAPTER 6 EXPERIMENTAL AND FINITE ELEMENT SIMULATION STUDIES OF SUPERPLASTIC BOX FORMING
113 CHAPTER 6 EXPERIMENTAL AND FINITE ELEMENT SIMULATION STUDIES OF SUPERPLASTIC BOX FORMING 6.1 INTRODUCTION Superplastic properties are exhibited only under a narrow range of strain rates. Hence, it
More informationReflector profile optimisation using Radiance
Reflector profile optimisation using Radiance 1,4 1,2 1, 8 6 4 2 3. 2.5 2. 1.5 1..5 I csf(1) csf(2). 1 2 3 4 5 6 Giulio ANTONUTTO Krzysztof WANDACHOWICZ page 1 The idea Krzysztof WANDACHOWICZ Giulio ANTONUTTO
More informationGeneral Method for Exponential-Type Equations. for Eight- and Nine-Point Prismatic Arrays
Applied Mathematical Sciences, Vol. 3, 2009, no. 43, 2143-2156 General Method for Exponential-Type Equations for Eight- and Nine-Point Prismatic Arrays G. L. Silver Los Alamos National Laboratory* P.O.
More informationInternational Conference on Space Optics ICSO 2008 Toulouse, France October 2008
ICSO 008 14 17 October 008 Edited by Josiane Costeraste, Errico Armandillo, and Nikos Karafolas Coupled thermo-elastic and optical performance analyses of a reflective baffle for the BepiColombo laser
More informationspecified or may be difficult to handle, we often have a tabulated data
Interpolation Introduction In many practical situations, for a function which either may not be explicitly specified or may be difficult to handle, we often have a tabulated data where and for In such
More informationFitting a Polynomial to Heat Capacity as a Function of Temperature for Ag. by
Fitting a Polynomial to Heat Capacity as a Function of Temperature for Ag. by Theresa Julia Zielinski Department of Chemistry, Medical Technology, and Physics Monmouth University West Long Branch, J 00764-1898
More informationThe Interpolating Polynomial
Math 45 Linear Algebra David Arnold David-Arnold@Eureka.redwoods.cc.ca.us Abstract A polynomial that passes through a given set of data points is called an interpolating polynomial. In this exercise you
More information