Hi

This is effectively strange

Are you sure that you are taking into account the current speed in the expression of f ?

I mean that you should have

Results 1 to 9 of 9

- Mar 7th 2009, 02:14 AM #1

- Joined
- Feb 2009
- Posts
- 17

## Can anybody help me with this?

I am currently calculating the numerical solution to a vehicle in motion stopping.

I have used an Euler method where the initial kinetic energy T is calculated with 1/2mv^2.

The Euler method then reduces the kinetic energy in terms of friction and wind resistance. Tn+1 = Tn + h*f(x)

where f(x) = B + Cv^2 and B, C are the constants for friction and wind.

After each iteration I have used the new kinetic energy to calculate the velocity for the next value. i.e Tn+1 = 1/2mv^2

Therefore

v = sqrt((2*Tn+1)/m)

The Euler loop then continues until the Kinetic energy is zero and nh is the distance.

The problem;

When plotting the results Kinetic energy against Distance I have a line. This should be a curve should it not?

Am I calculating the velocity correctly? Are there any problems with my calculations etc. Any help would be greatly appreciated

- Mar 7th 2009, 05:23 AM #2

- Joined
- Nov 2008
- From
- France
- Posts
- 1,458

- Mar 7th 2009, 07:15 AM #3

- Joined
- Feb 2009
- Posts
- 17

Yes, I am pretty sure. I have worked out the first three steps by hand and it matches the first three steps in the program. I'm not sure if I'm working out the velocity in the correct way for each step though.

After I reduce the Kinetic energy I then calculate the velocity based on the new kinetic energy. v = sqrt((2*T)/m) Is this correct? or is there another way to find the new velocity.

I have the results of the program if you want me to post them

- Mar 7th 2009, 07:22 AM #4

- Joined
- Nov 2008
- From
- France
- Posts
- 1,458

- Mar 7th 2009, 11:07 PM #5

- Joined
- Nov 2005
- From
- someplace
- Posts
- 14,972
- Thanks
- 5

- Mar 8th 2009, 12:46 AM #6

- Joined
- Feb 2009
- Posts
- 17

For an increment of h = 1

where B and C are always set to -3454.5 and -0.2514175 respectively.

B = umg

u is the coefficient of friction.

C =

T0 is 10ms^-1

Results

T:

46520.4

43042.5

39566.3

36091.9

32619.3

29148.4

25679.2

22211.8

18746.1

15282.2

11820

8359.58

4900.88

1443.92

-2011.31

S:

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

V:

9.64576

9.2782

8.89565

8.49611

8.07704

7.63523

7.16648

6.6651

6.12309

5.52851

4.8621

4.08891

3.13078

1.69936

-1.#IND

F = h*(B+Cv^2) where v= sqrt((2*T)/m) not actually to sure which T it is calculating Tn I assume. I think this is where the problem lies I just don't know how to go about solving it.

-3477.89

-3476.14

-3474.4

-3472.65

-3470.9

-3469.16

-3467.41

-3465.67

-3463.93

-3462.18

-3460.44

-3458.7

-3456.96

-3455.23

-1.#IND

- Mar 8th 2009, 05:42 AM #7

- Joined
- Nov 2008
- From
- France
- Posts
- 1,458

Your results are correct

If you calculate the speed you can find the formula with exponential

But due to the fact that

we can use the approximation which leads to the linear approximation of the square of the speed

Or

This gives a linear trend but this is only an approximation

Speed :

Zoom

Kinetic energy

Zoom

- Mar 18th 2009, 08:58 PM #8

- Joined
- Mar 2009
- Posts
- 31

Hmmmm wait though...if you're plotting Kinetic Energy vs. Distance you should obtain a line because Tn+1 = Tn + hf(x) is a linear relationship and ultimately your velocity plays the role of a constant in each new iteration...I might be following your logic incorrectly but I seem to see a straight line relationship :S

- Mar 18th 2009, 08:59 PM #9

- Joined
- Mar 2009
- Posts
- 31