Soving an eqn with exponential & linear parts

Ikeness
Messages
3
Reaction score
0

Homework Statement


Known Constants
A, B and C

Unknown (we are solving for)
t

e is the exponential as in ln(e) = 1

Homework Equations


e^{t*A} + t*B + C = 0
(we want to solve for t)


The Attempt at a Solution


This is not for school or anything.
This is actually part of an attempt at the solution of another problem. Reaching this equation I had to stop. I do not remember seeing something like this in Calculus or Differential Equations, so I need help.

My overall purpose is to find the initial velocity necessary to launch a particle to a specified location. This particle is subject to gravity and air drag. If any of you happen to have the solution to this (or something simular) on hand, that information would be worlds better than just getting me past this one step.

I have my progress somewhat documented in this ugly text file.
If you care to read on, feel free. But it is not necessary.

I started by looking at what these people were doing, and simply added gravity to the initial problem.
Except that I am dealing with a drag force -kv instead of -k*sqrt(v)
https://www.physicsforums.com/showthread.php?t=190387


My work:

ma=-kv + gm
a = (-k/m)v + g
a - g = (-k/m)v
dv/dt - g = (-k/m)v

1/dt - g/dv = (-k/m)v/dv
1/dt = [(-k/m)v + g] /dv
dv / [(-k/m)v + g] = dt
dv / [(-k/m)v + g] = dt
dv (-k/m) / [(-k/m)v + g] = (-k/m) dt
dv / [v + (-m/k)g] = (-k/m) dt

<integrate>
ln|v + (-m/k)g| = (-k/m)t + C

EDIT
Just want to point out a serious mistake made here.
It does not effect my problem too seriously. I would fix it except that I am now trying blochwave's suggestion about using -kv^2 instead of -kv for air drag.
<rase e^(everything)>
v + (-m/k)g = e^(-tk/m) + D
where (v > -gm/k) (D = e^C)

(if t = 0)
v_0 - gm/k = 1 + D
v_0 = gm/k + 1 + D
v_0 - gm/k - 1 = D

dx/dt = gm/k + e^(-tk/m) + D
dx = [gm/k + e^(-tk/m) + D] dt
<integrate>
x = tgm/k + -k/m e^(-tk/m) + Dt + F

(if t = 0)
0 = -k/m + F
F = k/m

x = tgm/k + -k/m e^(-tk/m) + Dt + k/m

<let us solve for t in the x = eqn>

x = tgm/k + -k/m e^(-tk/m) + Dt + k/m
x = t[gm/k - D] - k/m e^(-tk/m) + k/m
x - k/m = t[gm/k - D] - k/m e^(-tk/m)
xk/m - 1 = t[gmm/kk - Dm/k] - e^(-tk/m)


And this is where I am. My purpose here was to plug in the solution of t into the equation involving velocity. Hopefully then I can have an equation based on just velocity (v) and displacement (x) without any t parameterization mess.
But I cannot solve for this t.
 
Last edited:
Physics news on Phys.org
Equations with exponential and linear parts cannot be solved analytically in general. At best you can make a second order taylor approximation and solve the resulting quadratic.
 
You usually write the resisting force as F=kmv, not kv

so the y direction equation of motion is (' is dot! Denoting time derivative) my''=-kmy'-mg, the resistive force always opposes the velocity, hence its negative, and the gravitational force is down

dv/(kv+g)=-dt

1/k*ln(kv+g)=-t+C

kv+g=e^(-kt+kc)

knowing initial y velocity we assume is Vyi at t=0

Vy=-g/k+[k*Vyi+g]*e^(-kt)/k

Vy=dy/dt, so you can then integrate again to find y(t)

The x equation is similar but mx''=-kmx' and that equation is a little easier

That's for you to compare to your own work(which I see at a glance is different, but perhaps not woefully so), it's straight out of my classical mechanics textbook. Also mentioned in my classical mechanics textbook is that attempting to solve for things like time of flight and stuff can't be done analytically because the equations become transcedental, and unless you know much perturbation theory or have powerful numerical solving tools, you're done, and I bet that's the issue you're having right now

An easier problem to solve is when the resistive force is proportion to v SQUARED. Much cleaner solution and can be done analytically, and, like a linear dependence on v, is realistic for certain circumstances
 
blochwave said:
You usually write the resisting force as F=kmv, not kv

...

That's for you to compare to your own work(which I see at a glance is different, but perhaps not woefully so), it's straight out of my classical mechanics textbook. Also mentioned in my classical mechanics textbook is that attempting to solve for things like time of flight and stuff can't be done analytically because the equations become transcedental, and unless you know much perturbation theory or have powerful numerical solving tools, you're done, and I bet that's the issue you're having right now
Yes. Thanks. I was carrying around an extra m. I also really screwed up on that exponential step as I have noted above. But asside from that, our problems steps are very simular.

blochwave said:
An easier problem to solve is when the resistive force is proportion to v SQUARED. Much cleaner solution and can be done analytically, and, like a linear dependence on v, is realistic for certain circumstances
I like easier problems. Please let me know if your textbook has anything to say about this system. Thanks.

Running through with v SQUARED resistive force I stopped after the first integration step to make sure that I am doing it right. Please let me know.
my'' = -kmy'^2 - mg
y'' = -ky'^2 - g

dv/dt = -kv^2 - g
dv/(-kv^2 - g) = dt
dv/(kv^2 - g) = -dt

dv/(v^2 - g/k) = -kdt

<integrate>
using table entry from http://en.wikipedia.org/wiki/List_of_integrals_of_rational_functions
The one for the integral of
1/(x^2 - a^2)

I GET

-1/sqrt(g/k) * arctanh(v/sqrt(g/k)) (for |v| < |sqrt(g/k)|)
-1/sqrt(g/k) * arccoth(v/sqrt(g/k)) (for |v| > |sqrt(g/k)|)

OR

1/2sqrt(g/k) ln[ (sqrt(g/k) - v) / (sqrt(g/k) + v) ] (for |v| < |sqrt(g/k)|)
1/2sqrt(g/k) ln[-(sqrt(g/k) - v) / (sqrt(g/k) + v) ] (for |v| > |sqrt(g/k)|)


I do not know which to use. Please give me a suggestion.

Thanks.
-Joel
 
Sorry you typed all that out but

dv/(-kv^2 - g) = dt
dv/(kv^2 - g) = -dt


Look hard, what went wrong there?
 
Use this link to help

http://farside.ph.utexas.edu/teaching/336k/lectures/node30.html

Note that they just have a constant and v^2, the reason being realistically something's drag(which is what this models)doesn't really depend on mass but on the air it's traveling through, cross section, stuff like that, and notice that the initial problem he solves is for an object FALLING, which explains why he has the negatives different from us
 
Prove $$\int\limits_0^{\sqrt2/4}\frac{1}{\sqrt{x-x^2}}\arcsin\sqrt{\frac{(x-1)\left(x-1+x\sqrt{9-16x}\right)}{1-2x}} \, \mathrm dx = \frac{\pi^2}{8}.$$ Let $$I = \int\limits_0^{\sqrt 2 / 4}\frac{1}{\sqrt{x-x^2}}\arcsin\sqrt{\frac{(x-1)\left(x-1+x\sqrt{9-16x}\right)}{1-2x}} \, \mathrm dx. \tag{1}$$ The representation integral of ##\arcsin## is $$\arcsin u = \int\limits_{0}^{1} \frac{\mathrm dt}{\sqrt{1-t^2}}, \qquad 0 \leqslant u \leqslant 1.$$ Plugging identity above into ##(1)## with ##u...
Back
Top