Rigid rod moving in potential field

AI Thread Summary
A program is being developed to simulate a rigid rod moving in a two-dimensional potential field, where friction is negligible and the potential gradient is known. The force acting on the rod is derived from the potential field, specifically using the equation -∇P, where P represents potential energy. The simulation involves discretizing time and space, requiring calculations of total force and torque by breaking the rod into segments and treating each as a particle. The approach includes updating the rod's position and orientation iteratively based on calculated accelerations. The discussion clarifies that the rod cannot be treated as two particles connected by a massless rod and confirms that the order of rotation and translation does not affect the simulation.
kirzoaktrt
Messages
4
Reaction score
0
Hi,

I'd like to write a program that "simulates" a rigid rod of length L moving in a potential field. The problem is in two dimensions. Friction is assumed to be negligible. The potential field and its gradient is known at every point of the two-dimensional domain.

The "mass" of the rod and the magnitude of the potential gradient can be arbitrarily set. The rod has uniform density.

I understand that the force acting upon a particle in a potential field is -∇P, where P is the potential. This is an attraction of the particle towards lower potentials. What is the force acting on the rod?

Since this is a computer simulation, time and space are discretized. The simulation would proceed in discrete timesteps Δt. How do I compute the position of the rod in each timestep? How would you generally approach this problem?

Your help is greatly appreciated.
 
Physics news on Phys.org
How does the rod couple to the field?

For example, in an electric potential field, the units are V = J/C, i.e. energy per unit charge. What is the equivalent "charge" of the rod?
 
The force is -∇P at every point in the two-dimensional domain, while

P=qψ.

In this context, P is the potential energy, which is known and ψ is the potential. q would be a constant, which could be the mass, charge or whatever. I can set it to be anything, but it is a real constant. I set Δt and q such that the problem can be handled numerically, but they don't have any significance besides this.

Answering your question, it does not matter whether you treat this as a gravity or electric problem. Units are unimportant.
 
All you need is to figure the total force and total torque acting on the rod. These would be given by integrals, but in your computer program they will be calculated as sums. Break down the rod into N parts and for each part calculate the force as if it was a particle. calculate the torque using the center of mass as pivot. Add all the forces to find the total force and add all the torques to find the total torque. Use those to figure out the acceleration and angular acceleration of the rod. use those to update one time step to find what's the updated values for angular velocity and velocity. Use those to update the position of the center of mass and orientation (angle) of the rod. Rinse and repeat.
 
Thanks!

Two questions.

1, Without losing generality, can the rod be treated as two particles connected by a massless rod?
2, I assume that it does not matter whether I rotate or translate the rod first at each iteration, yes?
 
Answers:

1, No
2, Yes

Many thanks for the replies.
 
I have recently been really interested in the derivation of Hamiltons Principle. On my research I found that with the term ##m \cdot \frac{d}{dt} (\frac{dr}{dt} \cdot \delta r) = 0## (1) one may derivate ##\delta \int (T - V) dt = 0## (2). The derivation itself I understood quiet good, but what I don't understand is where the equation (1) came from, because in my research it was just given and not derived from anywhere. Does anybody know where (1) comes from or why from it the...

Similar threads

Back
Top