Fitting Element Fields
|
|
- Bernice Benson
- 6 years ago
- Views:
Transcription
1 F Fitting Element Fields F
2 Appendix F: FITTING ELEMENT FIELDS TABLE OF CONTENTS Page F. Introduction F F. Least-Square Fitting F F. Line Elements F 4 F.. Line Element Processing Modules F 6 F.. Line Self-Kernel Matrices F 6 F.. Line Coupled Matrices F 7 F.4 Lagrangian Triangles F 8 F.4. Lagrangian Triangle Processing Modules F 0 F.4. Lagrangian Triangle Self Kernels F F.4. Lagrangian Triangle Fitting Matrices F F
3 F. LEAST-SQUARE FITTING F.. Introduction This Appendix studies the problem of Polynomial Field Fitting over a finite element domain. This is a subset of the more general curve- and surface-fitting problem, as covered for example in [476. This special case is defined as follows. (I) A polynomial interpolation of a scalar function f, called the base function, isdefined over a finite element domain using n nodal values of f and n shape functions expressed in natural coordinates. (II) A polynomial interpolation of a scalar function g, called the fitted function,isdefined over the same element domain in terms of m n nodal values and m shape functions. This function approximates the base function f in a least square sense. The interpolations involved in (I) and (II) will be called the base and fitted interpolants of f and g, respectively. The case m = n is not excluded, as it covers situations in which the same polynomial order is defined from different nodal values. A well know example is the linear triangle with n = corner values (the Turner triangle) versus that with m = midpoint values (the Veubeke triangle). Passing from f to g involves a least square fitting. This is an operation that has been extensively studied in the applied mathematics literature. However, most of the work available there does not consider the use of natural coordinates as well as integration over the element domain. These two features are essential for finite element applications. Passing from g to f is a collocation operation, which simply involves evaluating g at the f nodes. F.. Least-Square Fitting Let the node values for the base and reduced interpolant be arranged in column vectors denoted by f and g of lengths n and m, respectively. The associated shape functions are arranged in row vectors N f and N g, respectively. The shape functions of both f and g will be assumed to be linearly independent. The interpolations may be expressed as the dot products f = N f f, g = N g g. (F.) The distance between g and f is d = g f = N g g N f f. Its square is d = (N g g N f f) T (N g g N f f) = g T N T g N gg g T N T g N f f + f T N T f N f f. (F.) As the distance norm we take s = d d = ( f T Af g T Bf+ g T Cg ) [ = A B [ ft g T T B C f. (F.) g in which denotes the element domain, and A = N T f N f d, B = Ng T N f d, C = N T g N g d. (F.4) F
4 Appendix F: FITTING ELEMENT FIELDS Matrices A and C are called the self kernels for f and g, respectively. They are symetric and positive definite (because the shape functions are assumed to be liearly independent). Matrix B of order m n is the coupled kernel, which is rectangular unless n = m. Both A and C are symmetric and positive definite, of orders n n and m m, respectively, whereas the m n matrix B is generally rectangular. Minimizing s with respect to g yields s g = Cg Bf= 0, whence g = C Bf= Gf. (F.5) G = C B will be called the f g fitting matrix. When identification of orders m and n are convenient, the matrix will be denoted G n m. Minimizing s with respect to f yields s f = Af BT g = 0, f = A B T g = Ff. (F.6) F = A B will be called the g f collocation matrix. It is a generalized inverse of G since GF= I, FG= H, (F.7) in which I is the m m identity matrix whereas H is an n n projector (idempotent) matrix satisfying H = H. Matrix H, which has rank m, is not generally symmetric. Replacing the expressions of G and F into GF= I yields the property C BA B T = I, or B T A B = C, (F.8) whereas replacing the expressions of G and F into FG= H yields the property A B T C B = H, or B T C B = AH. (F.9) Notice that B T A B and B T C B are the Schur complements of C and A, respectively, in the supermatrix shown in (F.). Since B T C B is symmetric, so is AH. It can be readily shown that the reverse product HAis antisymmetric. Inserting (F.5) into (F.) gives the minimum distance as s min = ft (A B T C B) f = ft (A AH) f = ft A Ĥf, (F.0) in which Ĥ = I H. Since H is a projection matrix, so is Ĥ. As a quick check, take N f N g, whence m = n and C = B = A;ifsoH = I and s min = 0, as can be expected. Readers familiar with computational dynamics may recognize A and C as the numerical part of the consistent mass matrices pertaining to the base and fitted elements, respectively. F 4
5 F. LINE ELEMENTS 4 LL LL LL LL4 LLG LLG HL4 L e x Figure F.. Line element notation. LL and HL denote Lagrangian and Hermitian interpolation, respectively. Nodes shown as circles have only one DOF: the function value there. Nodes pictured with a rotation ( harpoon ) symbol have two DOF there: the function value and its x-derivative scaled by the line length L e. LineShapeFunctions[typ_,ξ_,Le_:=Module[{Nsf=Null,m=Null,ξnc=Null}, If [typ=="ll", Nsf={}; m=0; ξnc={{0}}; If [typ=="ll", Nsf={(-ξ)/,(+ξ)/}; m=; ξnc={{-},{}}; If [typ=="ll", Nsf={-(-ξ)*ξ/,(+ξ)*ξ/,-ξ ^}; m=; ξnc={{-},{},{0}}; If [typ=="ll4", Nsf={(-ξ)*(9*ξ ^-),(+ξ)*(9*ξ ^-), 9*(*ξ-)*(ξ^-),9*(*ξ+)*(-ξ^)}/6; m=; ξnc={{-},{},{-/},{/}}; If [typ=="llg", Nsf={-Sqrt[*ξ,+Sqrt[*ξ}/; m=; ξnc={{-sqrt[/},{sqrt[/}}; If [typ=="llg", Nsf={ξ*(-Sqrt[5+5*ξ),6-0*ξ ^,ξ*(sqrt[5+5*ξ)}/6; m=; ξnc={{-sqrt[/5},{0},{sqrt[/5}}; If [typ=="hl4", Nsf={*(-ξ)^*(+ξ),(-ξ)^*(+ξ)*Le,*(+ξ)^*(-ξ), -(+ξ)^*(-ξ)*le}/8; m=; ξnc={{-},{-},{},{}}; Return[{Nsf,m,ξnc}; Figure F.. Line element shape function module. IntegrateOverLine[f_,ξ_,Le_:=Module[{fint}, fint=integrate[(le/)*f,{ξ,-,}; Return[fint; LineFit[ftyp_,gtyp_,ξ_,Le_:=Module[{Nf={},Ng={},ζf,ζg,mf,mg, Am,Bm,Cm,Ainv,Cinv,F,G,GF,FG}, {Nf,mf,ζf}=LineShapeFunctions[ftyp,ξ,Le; {Ng,mg,ζg}=LineShapeFunctions[gtyp,ξ,Le; If [Nf==Null, Print["Unknown ftyp=",ftyp; Return[Null; If [Ng==Null, Print["Unknown gtyp=",gtyp; Return[Null; Am=IntegrateOverLine[Transpose[{Nf}.{Nf},ξ,Le; Ainv=Inverse[Am; Bm=IntegrateOverLine[Transpose[{Ng}.{Nf},ξ,Le; BmT=Transpose[Bm; Cm=IntegrateOverLine[Transpose[{Ng}.{Ng},ξ,Le; Cinv=Inverse[Cm; G=Cinv.Bm; F=Ainv.BmT; FG=F.G; GF=G.F; Return[{Am,Bm,Cm,G,F,GF,FG}; Figure F.. Line element fitting module. F.. Line Elements Line elements are used for one-dimensional FEM models. In structural mechanics: bars, beams, shafts, etc. This Section tabulates results for the seven line elements shown in Figure F.. The F 5
6 Appendix F: FITTING ELEMENT FIELDS line segment is assumed straight with constant metric and lenth L e along x. As a consequence, the integrals over the line can be evaluated exactly, without resource to numerical integration. Elements labeled "LLn" are the standard Lagrangian ones, with function values at equally spaced n nodes, which include the end points if n. Elements labeled "LLnG" have internal function values given at n>= Gauss points. That labeled "HL4" is the standard Hermite cubic-interpolated element with 4 DOF: two end-node values and two length-scaled x-derivatives at the end nodes: [ f,(df/dx) L e, f,(df/dx) L e. F... Line Element Processing Modules The shape functions for the line elements shown in Figure F. are returned by the Mathematica module LineShapeFunctions listed in Figure F.. The module is invoked as { Nsf,m,ζ nc }=LineShapeFunctions[typ,ξ,Le (F.) in which typ One of the identifiers shown in Figure? supplied as a character string (e.g., "LL") ξ The isoparametric natural coordinate. Le The element length. Only used for type "HL4". The module returns Nfs A list of shape functions. m The polynomial order of the shape functions. ξnc A list of isoparametric coordinates of the nodes. If argument typ is not implemented, Nsf returns Null. The polynomial fit operation is carried out by the Mathematica module LTFit listed in Figure F.6. The module is invoked as in which ftyp,gtyp ξ Le { A,B,C,G,F,FG,GF }=LineFit[ftyp,gtyp,ξ,Le (F.) Identifiers of the base and fitted line element, respectively. The isoparametric natural coordinate. The line segment length. Only relevant if ftyp or gtyp is "HL4", but typically is set to to reduce clutter in the fitting matrices. The module returns a list of seven matrices: A, B, C, G, F, FG(which should be the identity), and GF = H (which should be a projector). These matrices were introduced in F.. The one-point line element is denoted simply by "LL" rather than "LLG" since no confusion can arise. These Gauss elements have specialized uses. The scaling by L e is to reduce clutter in the fitting formulas. Scaling is supressed when node rotations are used. F 6
7 F. LINE ELEMENTS F... Line Self-Kernel Matrices Self kernels are the matrices A and C defined in (F.4). Since C = A if N f N g, only A is listed. Its subscript identifies the line element. A LL = [, A LL = 6 A LL4 = 680 A LLG = 8, A LL = , A LLG = [ F... Line Coupled Matrices, A HL4 = L e 40 [ 4 4, 6 0, 0 56 L e 54 L e L e 4 L e L e L e 54 L e 56 L e L e L e L e 4 L e, (F.) LL to LL: B = [, G = [, F = [. (F.4) LL to LLG: B = [ a a a a, G = 6 [ a a a a, F = [ a a 4 a 4 a, (F.5) in which a = +, a =, a = +, and a 4 = +. LL to LL: B = 6 0, G = 0, F = [ 0 0. (F.6) LL to LL: B = 6 [ 4, G = 6 [ 4, F = (F.7) LL to LLG: B = b b , G = 6 b b 4 0 b b , F = b b 4 6 in which b = + 5, b = 5, b = 5 + 5, and b 4 = 5 5. F 7 [ b 4 b 4 b 4 4 b, (F.8) 0 6 0
8 Appendix F: FITTING ELEMENT FIELDS LL to LLG: B = a a 4 4, G = a 4 a 4 6 a a 4 4, F = a 4 a 4 in which a = +, a =, a = +, and a 4 = +. LL4 to LL: B = 0 LL4 to LL: B = , G = , G = , F = , F = 8 6 [ a a 4 a 4 a, (F.9) (F.0) 0 0. (F.) LL4 to LL: B = 8 [, G = 8 [, F =. (F.) LL4 to LLG: B = 70 h h h h , G = h h h 4 h 400 h5 h 6 h 7 h , F = h 6 h 5 h 8 h 54 7 h 9 6 h 0 h 0 6 h 9 h 44 h h 44 h (F.) in which h = , h = 55 5, h = , h 4 = , h 5 = 0+ 5, h 6 = 0 5, h 7 = , h 8 = , h 9 = , h 0 = , h = 5 + 5, and h = 5 5. LL4 to LLG: B = 40 [ g g g g 4 g g g 4 g, G = 0 [ g g g g 4 g g g 4 g, F = g 5 g 6 g 6 g 5 g 5 g 6 g 6 g 5 (F.4) in which g = 5 +, g = 5, g = , g 4 = 45 7, g 5 = +, and g 6 =. (HL4 coupled matrices to be inserted here) F 8.
9 F.4 LAGRANGIAN TRIANGLES LT LT LTM LT6 LT9 LT0 6 Figure F.4. Notation for Lagrangian triangles considered in F.. Figures taken from [. LTShapeFunctions[typ_,tcoor_:=Module[{ζ,ζ,ζ,Nsf=Null,m=Null, ζnc=null}, {ζ,ζ,ζ}=tcoor; If [typ=="lt", Nsf={}; m=0; ζnc={{0,0,0}}; If [typ=="lt", Nsf={ζ,ζ,ζ}; m=; ζnc={{,0,0},{0,,0},{0,,0}}; If [typ=="ltm", Nsf={ζ+ζ-ζ,-ζ+ζ+ζ,ζ-ζ+ζ}; m=; ζnc={{/,/,0},{0,/,/},{/,0,/}}; If [typ=="lt6", Nsf={ζ*(*ζ-),ζ*(*ζ-),ζ*(*ζ-), 4*ζ*ζ,4*ζ*ζ,4*ζ*ζ}; m=; ζnc={{,0,0},{0,,0},{0,0,}, {/,/,0},{0,/,/},{/,0,/}}; If [typ=="lt9", Nsf={ζ*(*ζ-)*(*ζ-),ζ*(*ζ-)*(*ζ-), ζ*(*ζ-)*(*ζ-),9*ζ*ζ*(*ζ-),9*ζ*ζ*(*ζ-), 9*ζ*ζ*(*ζ-),9*ζ*ζ*(*ζ-),9*ζ*ζ*(*ζ-), 9*ζ*ζ*(*ζ-)}/+9*{-,-,-,,,,,,}*ζ*ζ*ζ/4; m=; ζnc={{,0,0},{0,,0},{0,0,},{/,/,0},{/,/,0}, {0,/,/},{0,/,/},{/,0,/},{/,0,/}}; If [typ=="lt0", Nsf={ζ*(*ζ-)*(*ζ-),ζ*(*ζ-)*(*ζ-), ζ*(*ζ-)*(*ζ-),9*ζ*ζ*(*ζ-),9*ζ*ζ*(*ζ-), 9*ζ*ζ*(*ζ-),9*ζ*ζ*(*ζ-),9*ζ*ζ*(*ζ-), 9*ζ*ζ*(*ζ-),54*ζ*ζ*ζ}/; m=; ζnc={{,0,0},{0,,0},{0,0,},{/,/,0},{/,/,0},{0,/,/}, {0,/,/},{/,0,/},{/,0,/},{/,/,/}}; Return[{Nsf,m,ζnc}; Figure F.5. Module LTShapeFunctions that returns shape function information for the Lagrange triangles of Figure F.4. F.4. Lagrangian Triangles A Lagrangian triangle (LT) is one interpolated from node values only. They are primarily used for D problems of variational index, such as plane stress and heat conduction. This Section tabulates results for the six triangular elements shown in Figure F.4. It assumes flat triangles with constant F 9
10 Appendix F: FITTING ELEMENT FIELDS LTFit[ftyp_,gtyp_,tcoor_:=Module[{Nf={},Ng={},mf,mg,ζf,ζg, Am,Bm,Cm,Ainv,Cinv,F,G,GF,FG}, {Nf,mf,ζf}=LTShapeFunctions[ftyp,tcoor; {Ng,mg,ζg}=LTShapeFunctions[gtyp,tcoor; If [Nf==Null, Print["Unknown ftyp=",ftyp; Return[Null; If [Ng==Null, Print["Unknown gtyp=",gtyp; Return[Null; Am=IntegrateOverTriangle[Transpose[{Nf}.{Nf},,mf+mf; Ainv=Inverse[Am; Bm=IntegrateOverTriangle[Transpose[{Ng}.{Nf},,mg+mf; BmT=Transpose[Bm; Cm=IntegrateOverTriangle[Transpose[{Ng}.{Ng},,mg+mg; Cinv=Inverse[Cm; G=Cinv.Bm; F=Ainv.BmT; FG=F.G; GF=G.F; Return[{Am,Bm,Cm,G,F,GF,FG}; Figure F.6. Module LTFit that return the fitting matrices introduced in F.. metric, in which the area integrations can be performed exactly, without resource to numerical integration. F.4.. Lagrangian Triangle Processing Modules The shape functions for the LT shown in Figure F.4 are returned by the Mathematica module LTShapeFunctions listed in Figure F.5. The module is invoked as { Nsf,m,ζ nc }=LTShapeFunctions[typ,tcoor (F.5) in which typ One of the LT identifiers shown in Figure F.4 supplied as a character string (e.g., "LT") tcoor A -item list of triangular coordinate symbols, usually { ζ,ζ,ζ }. The module returns Nfs A list of shape functions. m The polynomial order of the shape functions. ζ nc A list of triangular coordinates of the nodes. If argument typ is not implemented, Nsf returns Null. The polynomial fit operation is carried out by the Mathematica module LTFit listed in Figure F.6. The module is invoked as { A,B,C,G,F,FG,GF }=LTFit[ftyp,gtyp,tcoor (F.6) in which ftyp,gtyp tcoor Identifiers of the base and fitted triangle, respectively. A -item list of triangular coordinate symbols. The module returns a list of seven matrices: A, B, C, G, F, FG(which should be the identity), and GF = H (which should be a projector). These matrices were introduced in F.. Remark F.. The shape functions for LT9 and LT0 were derived in the author s 966 thesis [. The four others, due to Argyris and Fraeijs de Veubeke, were published one year earlier. F 0
11 F.4 LAGRANGIAN TRIANGLES Remark F.. Module LTFit carries out triangle integrals through module IntegrateOverTriangle, which is described in Exercise 5. of Chapter 5. The area integration is based on the formula (?), which is valid for a constant metric triangle. F.4.. Lagrangian Triangle Self Kernels Self kernels are the matrices A and C defined in (F.4). Since C = A if N f N g, only A is listed. -Node Linear Triangle with Centroid Value. Labeled LT in Figure F.4. A = [. (F.7) -Node Linear Triangle With Corner Values. Labeled LT in Figure F.4. A = [. (F.8) -Node Linear Triangle With Midpoint Values. Labeled LTM in Figure F.4. A = I. (F.9) 6-Node Quadratic Triangle. Labeled LT6 in Figure F A = , (F.0) Node Cubic Triangle. Labeled LT9 in Figure F A = (F.) Node Cubic Triangle. Labeled LT0 in Figure F.4. F
12 Appendix F: FITTING ELEMENT FIELDS A = (F.) F.4.. Lagrangian Triangle Fitting Matrices Base: LT, Fitted: LT. Matrices A and C are given by (F.8) and (F.7), respectively. GF= I, H = FG= projection matrix of rank. B = G = [, F = [ T. (F.) Base: LT, Fitted: LTM. Matrices A and C are given by (F.8) and (F.9), respectively. B = 0 0, G = B, F =. (F.4) 6 0 GF= I, H = FG= projection matrix of rank. Base: LT6, Fitted: LT. Matrices A and C are given by (F.0) and (F.8), respectively. B = , G = 5, F = GF= I, H = FG= projection matrix of rank Base: LT6, Fitted: LTM. Matrices A and C are given by (F.0) and (F.9), respectively. B = 6 6, G = B, F = (F.5) (F.6) F
13 F.4 LAGRANGIAN TRIANGLES GF= I, H = FG= projection matrix of rank. Base: LT6, Fitted: LT. Matrices A and C are given by (F.0) and (F.7), respectively. B = G = [ 0 0 0, F = [ T. (F.7) GF= I, H = FG= projection matrix of rank. Base: LT9, Fitted: LT6. Matrices A and C are given by (F.) and (F.0), respectively B = , G = F = GF= I 6, H = FG= projection matrix of rank 6. Base: LT9, Fitted: LT. Matrices A and C are given by (F.) and (F.0), respectively. B = , G = , F = T GF= I, H = FG= projection matrix of rank., (F.8) (F.9) F
14 Appendix F: FITTING ELEMENT FIELDS Base: LT9, Fitted: LT. Matrices A and C are given by (F.) and (F.7), respectively. B = G = [ , 48 F = [ T. (F.40) GF= I, H = FG= projection matrix of rank. Base: LT0, Fitted: LT9. Matrices A and C are given by (F.) and (F.), respectively. B = 440 G = 70 F = , , (F.4) GF= I 9, H = FG= projection matrix of rank 9. F 4
15 F.4 LAGRANGIAN TRIANGLES Base: LT0, Fitted: LT6. Matrices A and C are given by (F.) and (F.0), respectively B = , G = , F = GF= I 6, H = FG= projection matrix of rank 6. Base: LT0, Fitted: LT. Matrices A and C are given by (F.) and (F.8), respectively. B = , G = , F = T GF= I, H = FG= projection matrix of rank. Base: LT0, Fitted: LT. Matrices A and C are given by (F.) and (F.7), respectively. B = G = 0 [ , F = [ T. GF= I, H = FG= projection matrix of rank.. (F.4) (F.4) (F.44) F 5
Implementation of Iso-P Triangular Elements
24 Implementation of Iso-P Triangular Elements 24 1 Chapter 24: IMPLEMENTATION OF ISO-P TRIANGULAR ELEMENTS TABLE OF CONTENTS Page 24.1 Introduction..................... 24 3 24.2 Gauss Quadrature for
More informationTriangular Plate Displacement Elements
Triangular Plate Displacement Elements Chapter : TRIANGULAR PLATE DISPLACEMENT ELEMENTS TABLE OF CONTENTS Page. Introduction...................... Triangular Element Properties................ Triangle
More informationIntroduction Linear system Nonlinear equation Interpolation
Interpolation Interpolation is the process of estimating an intermediate value from a set of discrete or tabulated values. Suppose we have the following tabulated values: y y 0 y 1 y 2?? y 3 y 4 y 5 x
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More information2. Polynomial interpolation
2. Polynomial interpolation Contents 2. POLYNOMIAL INTERPOLATION... 1 2.1 TYPES OF INTERPOLATION... 1 2.2 LAGRANGE ONE-DIMENSIONAL INTERPOLATION... 2 2.3 NATURAL COORDINATES*... 15 2.4 HERMITE ONE-DIMENSIONAL
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationQuintic beam closed form matrices (revised 2/21, 2/23/12) General elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
More informationShape Function Generation and Requirements
Shape Function Generation and Requirements Requirements Requirements (A) Interpolation condition. Takes a unit value at node i, and is zero at all other nodes. Requirements (B) Local support condition.
More informationCode No: RT41033 R13 Set No. 1 IV B.Tech I Semester Regular Examinations, November - 2016 FINITE ELEMENT METHODS (Common to Mechanical Engineering, Aeronautical Engineering and Automobile Engineering)
More informationNonconservative Loading: Overview
35 Nonconservative Loading: Overview 35 Chapter 35: NONCONSERVATIVE LOADING: OVERVIEW TABLE OF CONTENTS Page 35. Introduction..................... 35 3 35.2 Sources...................... 35 3 35.3 Three
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationBilinear Quadrilateral (Q4): CQUAD4 in GENESIS
Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms
More informationFinite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems
Finite Element Method-Part II Isoparametric FE Formulation and some numerical examples Lecture 29 Smart and Micro Systems Introduction Till now we dealt only with finite elements having straight edges.
More informationMath 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination
Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column
More informationv = v 1 2 +v 2 2. Two successive applications of this idea give the length of the vector v R 3 :
Length, Angle and the Inner Product The length (or norm) of a vector v R 2 (viewed as connecting the origin to a point (v 1,v 2 )) is easily determined by the Pythagorean Theorem and is denoted v : v =
More information12. Cholesky factorization
L. Vandenberghe ECE133A (Winter 2018) 12. Cholesky factorization positive definite matrices examples Cholesky factorization complex positive definite matrices kernel methods 12-1 Definitions a symmetric
More informationhomogeneous 71 hyperplane 10 hyperplane 34 hyperplane 69 identity map 171 identity map 186 identity map 206 identity matrix 110 identity matrix 45
address 12 adjoint matrix 118 alternating 112 alternating 203 angle 159 angle 33 angle 60 area 120 associative 180 augmented matrix 11 axes 5 Axiom of Choice 153 basis 178 basis 210 basis 74 basis test
More informationElementary linear algebra
Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The
More informationis an isomorphism, and V = U W. Proof. Let u 1,..., u m be a basis of U, and add linearly independent
Lecture 4. G-Modules PCMI Summer 2015 Undergraduate Lectures on Flag Varieties Lecture 4. The categories of G-modules, mostly for finite groups, and a recipe for finding every irreducible G-module of a
More informationChapter 2 Linear Transformations
Chapter 2 Linear Transformations Linear Transformations Loosely speaking, a linear transformation is a function from one vector space to another that preserves the vector space operations. Let us be more
More informationIV B.Tech. I Semester Supplementary Examinations, February/March FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80
www..com www..com Code No: M0322/R07 Set No. 1 IV B.Tech. I Semester Supplementary Examinations, February/March - 2011 FINITE ELEMENT METHODS (Mechanical Engineering) Time: 3 Hours Max Marks: 80 Answer
More informationReview Packet 1 B 11 B 12 B 13 B = B 21 B 22 B 23 B 31 B 32 B 33 B 41 B 42 B 43
Review Packet. For each of the following, write the vector or matrix that is specified: a. e 3 R 4 b. D = diag{, 3, } c. e R 3 d. I. For each of the following matrices and vectors, give their dimension.
More informationHONORS LINEAR ALGEBRA (MATH V 2020) SPRING 2013
HONORS LINEAR ALGEBRA (MATH V 2020) SPRING 2013 PROFESSOR HENRY C. PINKHAM 1. Prerequisites The only prerequisite is Calculus III (Math 1201) or the equivalent: the first semester of multivariable calculus.
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
More informationMath 301 Test III. Dr. Holmes. November 4, 2004
Math 30 Test III Dr. Holmes November 4, 2004 This exam lasts from 9:40 until 0:35. You may drop one problem, other than problem 7. If you do all problems, your best work will count. Books, notes, and your
More informationLecture 8: Assembly of beam elements.
ecture 8: Assembly of beam elements. 4. Example of Assemblage of Beam Stiffness Matrices. Place nodes at the load application points. Assembling the two sets of element equations (note the common elemental
More information1. General Vector Spaces
1.1. Vector space axioms. 1. General Vector Spaces Definition 1.1. Let V be a nonempty set of objects on which the operations of addition and scalar multiplication are defined. By addition we mean a rule
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationMath 314 Lecture Notes Section 006 Fall 2006
Math 314 Lecture Notes Section 006 Fall 2006 CHAPTER 1 Linear Systems of Equations First Day: (1) Welcome (2) Pass out information sheets (3) Take roll (4) Open up home page and have students do same
More informationApplied Numerical Analysis Quiz #2
Applied Numerical Analysis Quiz #2 Modules 3 and 4 Name: Student number: DO NOT OPEN UNTIL ASKED Instructions: Make sure you have a machine-readable answer form. Write your name and student number in the
More information7. Choice of approximating for the FE-method scalar problems
7. Choice of approimating for the FE-method scalar problems Finite Element Method Differential Equation Weak Formulation Approimating Functions Weighted Residuals FEM - Formulation Weak form of heat flow
More informationMatrix Theory. A.Holst, V.Ufnarovski
Matrix Theory AHolst, VUfnarovski 55 HINTS AND ANSWERS 9 55 Hints and answers There are two different approaches In the first one write A as a block of rows and note that in B = E ij A all rows different
More informationMITOCW MITRES2_002S10linear_lec07_300k-mp4
MITOCW MITRES2_002S10linear_lec07_300k-mp4 The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources
More information(II.B) Basis and dimension
(II.B) Basis and dimension How would you explain that a plane has two dimensions? Well, you can go in two independent directions, and no more. To make this idea precise, we formulate the DEFINITION 1.
More informationMathematics for Engineers. Numerical mathematics
Mathematics for Engineers Numerical mathematics Integers Determine the largest representable integer with the intmax command. intmax ans = int32 2147483647 2147483647+1 ans = 2.1475e+09 Remark The set
More informationLinear Algebra I for Science (NYC)
Element No. 1: To express concrete problems as linear equations. To solve systems of linear equations using matrices. Topic: MATRICES 1.1 Give the definition of a matrix, identify the elements and the
More informationJEPPIAAR ENGINEERING COLLEGE
JEPPIAAR ENGINEERING COLLEGE Jeppiaar Nagar, Rajiv Gandhi Salai 600 119 DEPARTMENT OFMECHANICAL ENGINEERING QUESTION BANK VI SEMESTER ME6603 FINITE ELEMENT ANALYSIS Regulation 013 SUBJECT YEAR /SEM: III
More informationTypical Problem: Compute.
Math 2040 Chapter 6 Orhtogonality and Least Squares 6.1 and some of 6.7: Inner Product, Length and Orthogonality. Definition: If x, y R n, then x y = x 1 y 1 +... + x n y n is the dot product of x and
More informationTest 2 Review Math 1111 College Algebra
Test 2 Review Math 1111 College Algebra 1. Begin by graphing the standard quadratic function f(x) = x 2. Then use transformations of this graph to graph the given function. g(x) = x 2 + 2 *a. b. c. d.
More informationBHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I
BHAR AT HID AS AN ENGIN E ERI N G C O L L E G E NATTR A MPA LL I 635 8 54. Third Year M E C H A NICAL VI S E M ES TER QUE S T I ON B ANK Subject: ME 6 603 FIN I T E E LE ME N T A N A L YSIS UNI T - I INTRODUCTION
More informationAPPENDIX 2.1 LINE AND SURFACE INTEGRALS
2 APPENDIX 2. LINE AND URFACE INTEGRAL Consider a path connecting points (a) and (b) as shown in Fig. A.2.. Assume that a vector field A(r) exists in the space in which the path is situated. Then the line
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Module - 01 Lecture - 13
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras (Refer Slide Time: 00:25) Module - 01 Lecture - 13 In the last class, we have seen how
More informationCa Foscari University of Venice - Department of Management - A.A Luciano Battaia. December 14, 2017
Ca Foscari University of Venice - Department of Management - A.A.27-28 Mathematics Luciano Battaia December 4, 27 Brief summary for second partial - Sample Exercises Two variables functions Here and in
More informationCalculus II - Basic Matrix Operations
Calculus II - Basic Matrix Operations Ryan C Daileda Terminology A matrix is a rectangular array of numbers, for example 7,, 7 7 9, or / / /4 / / /4 / / /4 / /6 The numbers in any matrix are called its
More informationFinite Elements. Colin Cotter. January 18, Colin Cotter FEM
Finite Elements January 18, 2019 The finite element Given a triangulation T of a domain Ω, finite element spaces are defined according to 1. the form the functions take (usually polynomial) when restricted
More informationOctober 25, 2013 INNER PRODUCT SPACES
October 25, 2013 INNER PRODUCT SPACES RODICA D. COSTIN Contents 1. Inner product 2 1.1. Inner product 2 1.2. Inner product spaces 4 2. Orthogonal bases 5 2.1. Existence of an orthogonal basis 7 2.2. Orthogonal
More informationMath 307 Learning Goals
Math 307 Learning Goals May 14, 2018 Chapter 1 Linear Equations 1.1 Solving Linear Equations Write a system of linear equations using matrix notation. Use Gaussian elimination to bring a system of linear
More informationMath 660-Lecture 15: Finite element spaces (I)
Math 660-Lecture 15: Finite element spaces (I) (Chapter 3, 4.2, 4.3) Before we introduce the concrete spaces, let s first of all introduce the following important lemma. Theorem 1. Let V h consists of
More informationCourse Summary Math 211
Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.
More informationNumerical Methods of Approximation
Contents 31 Numerical Methods of Approximation 31.1 Polynomial Approximations 2 31.2 Numerical Integration 28 31.3 Numerical Differentiation 58 31.4 Nonlinear Equations 67 Learning outcomes In this Workbook
More informationStudies on Barlow points, Gauss points and Superconvergent points in 1D with Lagrangian and Hermitian finite element basis
Volume 7, N., pp. 75 0, 008 Copyright 008 SBMAC ISSN 0101-805 www.scielo.br/cam Studies on Barlow points, Gauss points and Superconvergent points in 1D with Lagrangian and Hermitian finite element basis
More informationMatrix Basic Concepts
Matrix Basic Concepts Topics: What is a matrix? Matrix terminology Elements or entries Diagonal entries Address/location of entries Rows and columns Size of a matrix A column matrix; vectors Special types
More informationAlgebraic Complexity Theory
Peter Biirgisser Michael Clausen M. Amin Shokrollahi Algebraic Complexity Theory With the Collaboration of Thomas Lickteig With 21 Figures Springer Chapter 1. Introduction 1 1.1 Exercises 20 1.2 Open Problems
More informationCSCI 239 Discrete Structures of Computer Science Lab 6 Vectors and Matrices
CSCI 239 Discrete Structures of Computer Science Lab 6 Vectors and Matrices This lab consists of exercises on real-valued vectors and matrices. Most of the exercises will required pencil and paper. Put
More information1 9/5 Matrices, vectors, and their applications
1 9/5 Matrices, vectors, and their applications Algebra: study of objects and operations on them. Linear algebra: object: matrices and vectors. operations: addition, multiplication etc. Algorithms/Geometric
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 information1. Let m 1 and n 1 be two natural numbers such that m > n. Which of the following is/are true?
. Let m and n be two natural numbers such that m > n. Which of the following is/are true? (i) A linear system of m equations in n variables is always consistent. (ii) A linear system of n equations in
More informationForce System Resultants Distributed Loads. A hole has been found in the nudist camp wall. The police are looking into it.
Force System Resultants A hole has been found in the nudist camp wall. The police are looking into it. Objectives Understand the concept and representation of a distributed load Understand how to convert
More information1 Fields and vector spaces
1 Fields and vector spaces In this section we revise some algebraic preliminaries and establish notation. 1.1 Division rings and fields A division ring, or skew field, is a structure F with two binary
More informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
More informationOBJECTIVES & ASSUMPTIONS
OBJECTIVES & ASSUMPTIONS Z Project Objectives: Construct a template for the 5-node pyramid element Conduct higher-order patch tests Verify element formulation Determine the optimal element in bending X
More information(1)(a) V = 2n-dimensional vector space over a field F, (1)(b) B = non-degenerate symplectic form on V.
18.704 Supplementary Notes March 21, 2005 Maximal parabolic subgroups of symplectic groups These notes are intended as an outline for a long presentation to be made early in April. They describe certain
More information1 Matrices and vector spaces
Matrices and vector spaces. Which of the following statements about linear vector spaces are true? Where a statement is false, give a counter-example to demonstrate this. (a) Non-singular N N matrices
More informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationCOMPUTATIONAL ELASTICITY
COMPUTATIONAL ELASTICITY Theory of Elasticity and Finite and Boundary Element Methods Mohammed Ameen Alpha Science International Ltd. Harrow, U.K. Contents Preface Notation vii xi PART A: THEORETICAL ELASTICITY
More informationAlgorithms for Picture Analysis. Lecture 07: Metrics. Axioms of a Metric
Axioms of a Metric Picture analysis always assumes that pictures are defined in coordinates, and we apply the Euclidean metric as the golden standard for distance (or derived, such as area) measurements.
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More informationLecture 3: Linear Algebra Review, Part II
Lecture 3: Linear Algebra Review, Part II Brian Borchers January 4, Linear Independence Definition The vectors v, v,..., v n are linearly independent if the system of equations c v + c v +...+ c n v n
More information4.3 - Linear Combinations and Independence of Vectors
- Linear Combinations and Independence of Vectors De nitions, Theorems, and Examples De nition 1 A vector v in a vector space V is called a linear combination of the vectors u 1, u,,u k in V if v can be
More information(a) II and III (b) I (c) I and III (d) I and II and III (e) None are true.
1 Which of the following statements is always true? I The null space of an m n matrix is a subspace of R m II If the set B = {v 1,, v n } spans a vector space V and dimv = n, then B is a basis for V III
More informationLINEAR ALGEBRA W W L CHEN
LINEAR ALGEBRA W W L CHEN c W W L Chen, 1997, 2008. This chapter is available free to all individuals, on the understanding that it is not to be used for financial gain, and may be downloaded and/or photocopied,
More informationSupport Vector Machines. CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington
Support Vector Machines CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington 1 A Linearly Separable Problem Consider the binary classification
More informationJacobians of Matrix Transformations and Functions of Matrix Argument
Jacobians of Matrix Transformations and Functions of Matrix Argument A. M. Mathai Department of Mathematics & Statistics, McGill University World Scientific Singapore *New Jersey London Hong Kong Contents
More informationCoding the Matrix Index - Version 0
0 vector, [definition]; (2.4.1): 68 2D geometry, transformations in, [lab]; (4.15.0): 196-200 A T (matrix A transpose); (4.5.4): 157 absolute value, complex number; (1.4.1): 43 abstract/abstracting, over
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationx n -2.5 Definition A list is a list of objects, where multiplicity is allowed, and order matters. For example, as lists
Vectors, Linear Combinations, and Matrix-Vector Mulitiplication In this section, we introduce vectors, linear combinations, and matrix-vector multiplication The rest of the class will involve vectors,
More informationNumerical Integration for Multivariable. October Abstract. We consider the numerical integration of functions with point singularities over
Numerical Integration for Multivariable Functions with Point Singularities Yaun Yang and Kendall E. Atkinson y October 199 Abstract We consider the numerical integration of functions with point singularities
More informationVector Spaces. 9.1 Opening Remarks. Week Solvable or not solvable, that s the question. View at edx. Consider the picture
Week9 Vector Spaces 9. Opening Remarks 9.. Solvable or not solvable, that s the question Consider the picture (,) (,) p(χ) = γ + γ χ + γ χ (, ) depicting three points in R and a quadratic polynomial (polynomial
More informationLecture: Linear algebra. 4. Solutions of linear equation systems The fundamental theorem of linear algebra
Lecture: Linear algebra. 1. Subspaces. 2. Orthogonal complement. 3. The four fundamental subspaces 4. Solutions of linear equation systems The fundamental theorem of linear algebra 5. Determining the fundamental
More informationMATH 167: APPLIED LINEAR ALGEBRA Least-Squares
MATH 167: APPLIED LINEAR ALGEBRA Least-Squares October 30, 2014 Least Squares We do a series of experiments, collecting data. We wish to see patterns!! We expect the output b to be a linear function of
More informationLectures 9-10: Polynomial and piecewise polynomial interpolation
Lectures 9-1: Polynomial and piecewise polynomial interpolation Let f be a function, which is only known at the nodes x 1, x,, x n, ie, all we know about the function f are its values y j = f(x j ), j
More informationIntroduction to Matrices
POLS 704 Introduction to Matrices Introduction to Matrices. The Cast of Characters A matrix is a rectangular array (i.e., a table) of numbers. For example, 2 3 X 4 5 6 (4 3) 7 8 9 0 0 0 Thismatrix,with4rowsand3columns,isoforder
More informationAbstract. keywords: finite elements, beams, Hermitian interpolation, shape functions
Generation of shape functions for straight beam elements Charles E. Augarde Originally published in Computers and Structures, 68 (1998) 555-560 Corrections highlighted. Abstract Straight beam finite elements
More information1 Unitary representations of the Virasoro algebra
Week 5 Reading material from the books Polchinski, Chapter 2, 15 Becker, Becker, Schwartz, Chapter 3 Ginspargs lectures, Chapters 3, 4 1 Unitary representations of the Virasoro algebra Now that we have
More informationLinear Algebra, Spring 2005
Linear Algebra, Spring 2005 Solutions May 4, 2005 Problem 4.89 To check for linear independence write the vectors as rows of a matrix. Reduce the matrix to echelon form and determine the number of non-zero
More informationPost Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method
9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No
More informationUsing MATLAB and. Abaqus. Finite Element Analysis. Introduction to. Amar Khennane. Taylor & Francis Croup. Taylor & Francis Croup,
Introduction to Finite Element Analysis Using MATLAB and Abaqus Amar Khennane Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup, an informa business
More informationDHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS PART A (2 MARKS)
DHANALAKSHMI COLLEGE OF ENGINEERING, CHENNAI DEPARTMENT OF MECHANICAL ENGINEERING ME 6603 FINITE ELEMENT ANALYSIS UNIT I : FINITE ELEMENT FORMULATION OF BOUNDARY VALUE PART A (2 MARKS) 1. Write the types
More informationChapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of
Chapter 2 Linear Algebra In this chapter, we study the formal structure that provides the background for quantum mechanics. The basic ideas of the mathematical machinery, linear algebra, are rather simple
More information22m:033 Notes: 6.1 Inner Product, Length and Orthogonality
m:033 Notes: 6. Inner Product, Length and Orthogonality Dennis Roseman University of Iowa Iowa City, IA http://www.math.uiowa.edu/ roseman April, 00 The inner product Arithmetic is based on addition and
More informationEigenvalues and Eigenvectors
Contents Eigenvalues and Eigenvectors. Basic Concepts. Applications of Eigenvalues and Eigenvectors 8.3 Repeated Eigenvalues and Symmetric Matrices 3.4 Numerical Determination of Eigenvalues and Eigenvectors
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More information1.1 Bound and Free Vectors. 1.2 Vector Operations
1 Vectors Vectors are used when both the magnitude and the direction of some physical quantity are required. Examples of such quantities are velocity, acceleration, force, electric and magnetic fields.
More informationMath 344 Lecture # Linear Systems
Math 344 Lecture #12 2.7 Linear Systems Through a choice of bases S and T for finite dimensional vector spaces V (with dimension n) and W (with dimension m), a linear equation L(v) = w becomes the linear
More informationNumerical Mathematics
Alfio Quarteroni Riccardo Sacco Fausto Saleri Numerical Mathematics Second Edition With 135 Figures and 45 Tables 421 Springer Contents Part I Getting Started 1 Foundations of Matrix Analysis 3 1.1 Vector
More informationMA 323 Geometric Modelling Course Notes: Day 07 Parabolic Arcs
MA 323 Geometric Modelling Course Notes: Day 07 Parabolic Arcs David L. Finn December 9th, 2004 We now start considering the basic curve elements to be used throughout this course; polynomial curves and
More information