1 / 10

NUMERICAL INTEGRATION

NUMERICAL INTEGRATION. Stiffness matrix and distributed load calculations involve integration over the domain In many cases, analytical integration is very difficult Numerical integration based on Gauss Quadrature is commonly used in finite element programs Gauss Quadrature :

marcellus
Download Presentation

NUMERICAL INTEGRATION

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. NUMERICAL INTEGRATION • Stiffness matrix and distributed load calculations involve integration over the domain • In many cases, analytical integration is very difficult • Numerical integration based on Gauss Quadrature is commonly used in finite element programs • Gauss Quadrature: • Integral is evaluated using function values and weights. • si: Gauss integration points, wi: integration weights • f(si): function value at the Gauss point • n: number of integration points.

  2. ONE INTEGRATION POINT • Constant Function: f(s) = 4 • Use one integration point s1 = 0 and weight w1 = 2 Why? • The numerical integration is exact. • Linear Function: f(s) = 2s + 1 • Use one integration point s1 = 0 and weight w1 = 2 • The numerical integration is exact. • One-point Gauss Quadrature can integrate constant and linear functions exactly.

  3. TWO POINTS AND MORE • Quadratic Function: f(s) = 3s2 + 2s + 1 • Let’s use one-point Gauss Quadrature • One-point integration is not accurate for quadratic function • Let’s use two-point integration with w1 = w2 = 1 and -s1 = s2 = • Gauss Quadrature points and weights are selected such that n integration points can integrate (2n – 1)-order polynomial exactly.

  4. GAUSS QUADRATURE POINTS AND WEIGHTS • What properties do positions and weights have?

  5. General interval • In general we will have • From Wikipedia

  6. t t t s s s (a) 11 (b) 22 (c) 33 TWO DIMENSIONAL INTEGRATION • multiplying two one-dimensional Gauss integration formulas • Total number of integration points = m×n.

  7. EXAMPLE 6.9 • Integrate the following polynomial: • One-point formula • Two-point formula

  8. EXAMPLE 6.9 • 3-point formula • 4-point formula • 4-point formula yields the exact solution. Why?

  9. Quiz-like problems • When can’t you use Gaussian quadrature and must fall back on methods such as Simpson’s rule? • Estimate the integral with one and two points and compare to the exact value • What are the coordinates of the integration point at the bottom right corner of the 3x3 Gaussian quadrature points shown in the figure • Solution on the notes page

  10. APPLICATION TO STIFFNESS MATRIX • Application to Stiffness Matrix Integral

More Related