2D mass spring damper square lattice

Equations of motion for a 2D mass spring damper system and code for modelling on matlab needed.

[William Rose]

@Muhammad Ali, Please provide more details. Show us what have you done so far on this problem.

[Muhammad Ali]

@William Rose So far I only know how to model a 1D MSD system. I actually need to model a 2D MSD system which will stimulate the surface of water. I can't figure out how to start with a 2D system, because I can't find how to make the equations of motion for that. Once I get the equations of motion, then I can actually start thinking about how to code it. Any kind of help will be much appreciated. Thanks.

[William Rose]

@Muhammad Ali,

The equation of motion for a M-S-D system may be written as

where is a position vector, is frictional drag, is spring force, and is the external force (if any). may have 1 or 2 or 3 dimension (or more, I suppose). We assume that the spring pulls the object directly toward the origin, with a force proportional to the distance from the origin, and that frictional drag is the direction oppostite to the velocity, and proportional to the magnitude of velocity. These assumptions can be written as

and

Therefore we get the equation

In the 2-D case, we have, in the x-direction,

which looks totally normal. It involves the x-components of the spring force and the drag force. But here's where it gets tricky and then simple again. What are the equations for and ? It turns out the equations are

and

So we end up with

You get a similar result for the y-component:

So we see that for this special case - a linear mass-spring-dashpot system - the x- and y-parts do not affect each other, i.e. the equations of motion are totally separable. These two equations are what you need to solve. Since they are not coupled, you can solve them separately, in either order.

[Akhil Narayan]

Hello Team,

Could you please help me to find the Equation of motion for a 1D Spring damper system and code for modelling on matlab. The equation should be with n number of segments/ nodes. Each node represents a connection between two dimers on Mictorubule. Thus, the node has a displacement of X and the other node has a displacement of X2. In this case, it should resemble a linear chain of microtubule components (attached a model for reference). Please do needfull.

Thanks,

Answer 1

[William Rose]

@Muhammad Ali,

The system in the image you attached is obviously a lot more complicated. Masses m1 and m3 have 6 attached elements and therefore 6 forces. Each force may have x- an y-components. Masses m2 and m4 have 4 attached elements and therefore 4 forces. For m1:

where are the x and y components of the spring force from the wall acting on m1. are the x and y components of the dashpot force from the wall acting on m1. And so on. You will need geometry and calculus and physics to express each spring force and each dashpot force in terms of the x,y coordinates of the masses. If you want to solve the system in Matlab, you will need to convert the second order differetial equations above to coupled first order differential equations. Each mass will have four state variables. For m1, the state variables are . Likewise for the other three masses. So the state vector for the whole system will have 16 elements. This is explained well, and with examples, in the Matlab help on solving ordinary differential equations.

The attached document shows how each of the twelve forces above (on the right hand side of the equations) can be expressed in terms of the state varibles , etc. It also explains how to compute all the other forces that you need to solve the system.

[Muhammad Ali]

@William Rose thank you very much. This is extremely helpfull.

[William Rose]

You're welcome, @Muhammad Ali.

Answer 2

[William Rose]

@Muhammad Ali,

The equation of motion for a M-S-D system may be written as

where is a position vector, is frictional drag, is spring force, and is the external force (if any). may have 1 or 2 or 3 dimension (or more, I suppose). We assume that the spring pulls the object directly toward the origin, with a force proportional to the distance from the origin, and that frictional drag is the direction oppostite to the velocity, and proportional to the magnitude of velocity. These assumptions can be written as

and

Therefore we get the equation

In the 2-D case, we have, in the x-direction,

which looks totally normal. It involves the x-components of the spring force and the drag force. But here's where it gets tricky and then simple again. What are the equations for and ? It turns out the equations are

and

So we end up with

You get a similar result for the y-component:

So we see that for this special case - a linear mass-spring-dashpot system - the x- and y-parts do not affect each other, i.e. the equations of motion are totally separable. These two equations are what you need to solve. Since they are not coupled, you can solve them separately, in either order.

[Muhammad Ali]

@William Rose thanks alot. I understant it for a single mass moving in x and y direction but, how would it work for a system which I have attached in this comment. Please have a look. Thanks.