### Comparing Results from Dyn. Analysis

I have attached images for our results and the correct results for simplePend as verified by ADAMS and other sources.

If anyone has any ideas on how to fix this, let me know!

-James

**Moderators:** Dan Negrut, RaduS

After working on the code for Project 2, I decided to test my simEngine with a simple model like doublePend from HW9. I am getting results that are very similar to the correct results but not exactly correct. It seems like the solution diverges at a point which causes the results after that point to be erroneous. I am not sure what to do about this problem, and my partner is equally as lost.

I have attached images for our results and the correct results for simplePend as verified by ADAMS and other sources.

If anyone has any ideas on how to fix this, let me know!

-James

I have attached images for our results and the correct results for simplePend as verified by ADAMS and other sources.

If anyone has any ideas on how to fix this, let me know!

-James

This is hard.

Here's my suggestion, might not make sense.

take the simple pendulum example and prescribe a motion to the angle. Like phi1=2*pi*time.

Then you know how the pendulm moves: it's position and its velocity, right? The angular velocity is 2*pi, for instnace, from where you can get xdot and ydot. the same for x, y, phi.

then run your simEngine2D. You should step through your code in the debugger in MATLAB and see if you get after one time step the values that you expect.

I hope this helps.

dan

Here's my suggestion, might not make sense.

take the simple pendulum example and prescribe a motion to the angle. Like phi1=2*pi*time.

Then you know how the pendulm moves: it's position and its velocity, right? The angular velocity is 2*pi, for instnace, from where you can get xdot and ydot. the same for x, y, phi.

then run your simEngine2D. You should step through your code in the debugger in MATLAB and see if you get after one time step the values that you expect.

I hope this helps.

dan

James and I have been able to resolve this issue. Our plots are now as follows:

Can you share what your issue was? My plots look like your originals. I'm suspecting that I may have the same issue.

Greg

Greg

Hi Greg,

In our case, we needed to express a point force that is with respect to the LRF to an expression that is with respect to the GRF. In order to do this, we multiplied the force by the rotation matrix, and used the instantaneous angle of the body of interest to compute the rotation matrix. It is possible that you are calculating your rotation matrix based on the angle that is noted in the ADM file. Instead, try computing the rotation matrix based on the angle that you solve via the newmark equations. I hope this helps.

-Dan

In our case, we needed to express a point force that is with respect to the LRF to an expression that is with respect to the GRF. In order to do this, we multiplied the force by the rotation matrix, and used the instantaneous angle of the body of interest to compute the rotation matrix. It is possible that you are calculating your rotation matrix based on the angle that is noted in the ADM file. Instead, try computing the rotation matrix based on the angle that you solve via the newmark equations. I hope this helps.

-Dan

5 posts
• Page **1** of **1**

Return to ME451 Fall 2014: Kinematics and Dynamics of Machine Systems

Users browsing this forum: No registered users and 5 guests