Scribd
Upload
371 views4 pages

C4.newmark Beta Method Example5.4 of Chopra

This document describes using Newmark's Beta Method to solve for the displacement of a linear system over time. It defines constants like stiffness (k), mass (m), damping (c), time step (Δt), and initial position (x). Newmark's Method coefficients β and γ are set to 1/4 and 1/2. The algorithm calculates displacements (u), velocities (udot), and accelerations (udotdot) at each time step by solving the system of equations defined by Newmark's Method. The resulting displacement over time is plotted.

Uploaded by

podrow
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
371 views4 pages

C4.newmark Beta Method Example5.4 of Chopra

This document describes using Newmark's Beta Method to solve for the displacement of a linear system over time. It defines constants like stiffness (k), mass (m), damping (c), time step (Δt), and initial position (x). Newmark's Method coefficients β and γ are set to 1/4 and 1/2. The algorithm calculates displacements (u), velocities (udot), and accelerations (udotdot) at each time step by solving the system of equations defined by Newmark's Method. The resulting displacement over time is plotted.

Uploaded by

podrow
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 4

Force Function

kip := 1000lb

0
5
8.6602
10
m := 0.2533

8.6602
5
k := 10
x :=

0
0
c := 0.1592
0

0
0
0

i := 0 .. 10
t := 0.1

10

6
xi
4

0.2

0.4

0.6
i t

0.8

Newmark's Beta Method:Linear Systems:

uuu ( , , m , c , k , t , x ) :=

kbar k +
1

t
1

m +

c +

1 c

m + t

u 0
0

udot 0
0

x c udot k u
0

udotdot

for i 0 .. 10

Pi x

i+ 1

Pbari Pi + a udot + b udotdot


i

u i

Pbari
kbar

udoti

i
t

udotdoti

1
t

i+ 1

1
4

:=

1
2

u i

udot

udot + udoti

udotdot

i+ 1

umax u

:=

udot + t 1 udotdot

i
i
2

1
t

udot

u + u i
i

i+ 1

udotdot + udotdoti
i

1
2

udotdot

disp := uuu ( , , m , c , k , t , x )
0
0

0.044

0.233

0.612

1.083

disp = 5

1.431

1.423

0.962

0.191

-0.604

10 -1.144
11 -1.254

1
( disp) i

0.2

0.4

0.6
i t

0.8

You might also like