Scribd
Upload
64 views

Lecture 16: Numerical Solution

1. Numerical methods are required to solve complex vibration problems over time. Accuracy and stability are important factors for numerical methods. 2. Newmark's constant average acceleration method assumes acceleration is constant over each time step. It is unconditionally stable. 3. Newmark's linear acceleration method assumes acceleration varies linearly over each time step. It is conditionally stable for a time step of 1/10 the natural period for single degree of freedom systems.

Uploaded by

Le Minh Thanh
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
64 views

Lecture 16: Numerical Solution

1. Numerical methods are required to solve complex vibration problems over time. Accuracy and stability are important factors for numerical methods. 2. Newmark's constant average acceleration method assumes acceleration is constant over each time step. It is unconditionally stable. 3. Newmark's linear acceleration method assumes acceleration varies linearly over each time step. It is conditionally stable for a time step of 1/10 the natural period for single degree of freedom systems.

Uploaded by

Le Minh Thanh
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 11

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

Lecture 16: Numerical Solution


Reading materials: Section 5.3

1. Introduction
For complex loading time histories, the closed-form solutions become impossible to obtain and therefore we must resort to numerical methods. All numerical methods compute solution at discrete time steps and are based on some assumption regarding the solution over a given time interval. The choice of a suitable time step is critical. It is important to understand Accuracy and Stability of numerical methods. An accurate numerical solution is close to the exact solution of the differential equation. The stability refers to the largest time step that can be used without solution becoming unbounded due to accumulation of errors. An unconditionally stable method results in the solution staying bounded even with very large time step. For conditionally stable methods, the stability criteria are generally defined in terms of natural frequencies or period of vibration. Equations of motion

The solution at time ti is known: The solution at time ti+1 is unknown:

-1-

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

2. Newmarks constant average acceleration method


The acceleration is assumed to be constant over the interval time.

Numerically updates from ti to ti+1 At time ti, the acceleration, velocity and displacement are known. The force is prescribed.

-2-

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

For Multiple degree of freedom systems

3. Newmarks linear acceleration method


The acceleration is assume to be linear over the interval

-3-

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

Numerically updates from ti to ti+1 At time ti, the acceleration, velocity and displacement are known. The force is prescribed.

For Multiple degree of freedom systems

4. General Newmarks method

-4-

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

Solution procedure in the incremental form

-5-

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

Solution algorithm

-6-

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

5. Example

-7-

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

6. Stability
Constant average acceleration method: Unconditionally stable. Use a step size based on a trade-off between the desired accuracy and computational effort.

Linear acceleration method: For single degree of freedom system, solution is stable when t=T/10.

-8-

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

7. Finite difference method (Optional)


Approximation to derivatives

-9-

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

Central difference method for SDOF systems

Example

- 10 -

53/58:153 Lecture 16 Fundamental of Vibration ______________________________________________________________________________

Finite difference method for multidegree of freedom systems

8. Runge-Kutta Method for SDOF systems (Optional)


Runge-Kutta method

Example

- 11 -

You might also like