class ForceField
{
public:
virtual vector3f GetAcceleration(vector3f position, float mass, float t) const = 0;
};
class RadialForceField : public ForceField
{
public:
RadialForceField(float _magnitude)
{
magnitude = _magnitude;
}
virtual vector3f GetAcceleration(vector3f position /* of particle */, float mass /* of particle */, float time) const
{
///// TODO: Help Here ////////////
float distance = position.Length();
return (-magnitude / (distance * distance)) * vector3f(1,1,1);
////////////////////
}
float magnitude;
};
PS, Radial Force Fields
Author
Hello all,
My good friend and I are in the process of creating a particle system. And we both thought it would be nice to have a Radial force field effect - Where each particle is attracted to the origin (0,0,0). I was thinking, we would could achieve this by only using the current position of the particle. And where the strength of the force is proportional to the distance from the origin.
Our abstrast Force field interface has worked well so far.
Well, so far here is my attempt, alas it fails. Does anyone here have the "know-how" to answer my post?
As always, thanks in advance.
I'm not in the mood to analyze your code intensively, but I guess the problem lays in the (1,1,1) vector. You scale this vector according to some distance related scalar, so all resulting vectors lay on the same line (the line that crosses the zero vector and 1,1,1). What you actually want (I guess) is something like:
velocity=(pos-fieldorigin)*some_scalar
where some_scalar is based on the distance.
good luck!
velocity=(pos-fieldorigin)*some_scalar
where some_scalar is based on the distance.
good luck!
Author
Thanks. I tried that, and this version, but still no luck...
Any other ideas?
We're attempting and using mit ocw to get started...
http://groups.csail.mit.edu/graphics/classes/6.837/F04/assignments/assignment9/
Help...
Any other ideas?
We're attempting and using mit ocw to get started...
http://groups.csail.mit.edu/graphics/classes/6.837/F04/assignments/assignment9/
Help...
class RadialForceField : public ForceField{public: RadialForceField(float m) { magnitude = m; } virtual vector3f GetAcceleration(vector3f position, float mass, float t) const { float radius = position.Length(); return -(magnitude / (radius * radius)) * position; } float magnitude;};
Author
It generates a radial force field, but it is not quite "radial" eg a circle. More like an oval...
Ideas?
Ideas?
It's still not clear exactly what you're observing the current code doing. How are you testing the force field? Are you evaluating the acceleration at various points to draw a representation of the acceleration vectors, or are you dropping a bunch of particles and watching where they are going?
Hope this helps:
C = Position of the center of attraction (vector) P = Position of the object (vector) r = distance between C and P D = direction from the object to the center of attraction (vector) k = "magnitude" of the attraction F = force on the object (vector) m = mass of the object A = acceleration of the object (vector) r = length( C - P ) D = ( C - P ) / r F = D * k / r2 = ( C - P ) * k / r3 A = F / m = ( C - P ) * k / ( r3 * m )
Author
AP, yes, i'm using both methods to test: drawing a representation of the acceleration vectors and dropping a bunch of particles and watching where they are going...
JohnBolton, your post makes perfect sense, but I'm not reproducing the desired results / effect.
image: http://homepage.cs.latrobe.edu.au/sjflannery/1.gif
Using:
where magnitude = k
JohnBolton, your post makes perfect sense, but I'm not reproducing the desired results / effect.
image: http://homepage.cs.latrobe.edu.au/sjflannery/1.gif
Using:
virtual vector3f GetAcceleration(vector3f position, float mass, float t) const { static vector3f origin(0,0,0); float radius = (origin - position).Length(); return origin - position * magnitude / (radius * radius * radius * mass); }where magnitude = k
I don't see what's wrong with it. Basically you are simulating gravitational effects on a point mass (force inversely proportional to distance squared). The path of a particle in this force field is supposed to be elliptical, just like the path of planets around the sun.
If you want the particles to move in a circular path around the point, you need to also apply a force that will make the particle's direction perpendicular to the displacement between the particle and the point. This force would depend on the velocity and position of the particle.
If you want the particles to move in a circular path around the point, you need to also apply a force that will make the particle's direction perpendicular to the displacement between the particle and the point. This force would depend on the velocity and position of the particle.
Perhaps the problem is that
return origin - position * magnitude / (radius * radius * radius * mass); return ( origin - position ) * magnitude / (radius * radius * radius * mass); This topic is closed to new replies.
Advertisement
Popular Topics
Advertisement
