Scribd
Upload
48 views

ENGR 2200 - Numerical Methods in Engineering I - EXCEL Assignment - Automatic Calculation

This document discusses calculating projectile motion parameters and performing numerical integration using the rectangle method. It provides the initial parameters for a projectile problem, including initial position, velocity, angle, and gravity. It also outlines calculating velocity and position at different time intervals. The document then discusses numerical integration, providing the function, bounds, and options for calculating the integral using the rectangle method.

Uploaded by

Ramesh K
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 XLS, PDF, TXT or read online on Scribd
48 views

ENGR 2200 - Numerical Methods in Engineering I - EXCEL Assignment - Automatic Calculation

This document discusses calculating projectile motion parameters and performing numerical integration using the rectangle method. It provides the initial parameters for a projectile problem, including initial position, velocity, angle, and gravity. It also outlines calculating velocity and position at different time intervals. The document then discusses numerical integration, providing the function, bounds, and options for calculating the integral using the rectangle method.

Uploaded by

Ramesh K
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 XLS, PDF, TXT or read online on Scribd
You are on page 1/ 52

ENGR 2200 - Numerical methods in engineering I - EXCEL Assignment - Automatic calculation

Projectile motion parameters

Initial position, x0 = 2 m Initial position, y =


Initial velocity, v0 = 25 m/s Initial angle, q0 =
Gravity, g= 9.806 m/s2

Table 1. Calculation of velocity and position


t(s) vx(m/s) vy(m/s) v(m/s) x(m) y(m) r(m)
Clear table

Create new table of results


1.000 21.65 2.69 21.82 23.65 10.60 25.92
Automatic calculation

3 m
30 o

Clear table

ate new table of results


Generate a table of values of z = f(x,y)
0 0.1 0.2 0.3 0.4 0.5
Generate Table
0 0 0.1 0.2 0.3 0.4 0.5
0.1 0.1 0.141421 0.223607 0.316228 0.412311 0.509902
xmin: 0 0.2 0.2 0.223607 0.282843 0.360555 0.447214 0.538516
xmax: 0.5 0.3 0.3 0.316228 0.360555 0.424264 0.5 0.583095
Dx: 0.1 0.4 0.4 0.412311 0.447214 0.5 0.565685 0.640312
0.5 0.5 0.509902 0.538516 0.583095 0.640312 0.707107
ymin: 0
ymax: 0.5
Dy: 0.1
0.5
Numerical integration using rectangles
f(x) = exp(-x^2/2)/sqrt(2*Pi)

a= 0 Type Integral
b= 1 Left
Dx = 0.05 Middle
Right

Clear Table Calculate Integrals


[9] The figure to the left shows a typical trapezoidal cross-section in an o
flow. The following geometric elements are apparent from the figure:
y = flow depth, z = side slope (zH:1v), T = top width, A = cross-sectional
P = wetted perimeter. The equations relating these geometric elements

T =b+ 2 zy , P=b +2 y √ 1+z 2 , A=(b+ zy ) y


If a flow of Q (volume per unit time) is flowing in the channel, and
Parameters for the critical of gravity, the conditions of critical flow are given by the equation:
flow solution
2
Q T
Calculate =1
Q= 0.20 m3/s gA 3
g= 9.81 m/s2
b= 1.20 m The depth of flow that satisfy the previous equation is known as the critic
z= 0.50 determine the critical depth for a trapezoidal channel, we will use the sec
y= 3.00 m (initial guess) to solve the following equation:

ycrit = 0.13918 m (solution) 2 3


f ( y )=Q T ( y )−g[ A ( y )] =0

The recursion equation to solve f(y) = 0 through the secant method is given by:

f ( yn ) where f'(yn) is approximated by: f ( y n )−f ( y n−1 )


y n+1= y n − , f ' ( y n )≈
f ' ( yn ) y n − y n−1
pezoidal cross-section in an open-channel
e apparent from the figure: b = bottom width,
op width, A = cross-sectional area, and
ng these geometric elements are the following:

y √ 1+z 2 , A=(b+ zy ) y
ng in the channel, and g is the acceleration
given by the equation:

equation is known as the critical depth. To


al channel, we will use the secant method

3
−g[ A ( y )] =0

f ( y n )−f ( y n−1 )
( y n )≈
y n − y n−1

You might also like