(adsbygoogle = window.adsbygoogle || []).push({}); 1. The problem statement, all variables and given/known data

Hi everyone! I'm trying to solve a problem involving 3 equal masses (m) arranged in a 45-45-90 right triangle, each mass connected to the other two by a spring, and all the springs have the same spring constant (k). So it should look like this:

Where I have made a the equilibrium length.

My goal is to set up the eigenvalue equation for its normal modes and find them.

2. Relevant equations

The Euler Lagrange equations from the Lagrangian we will set up should lead to:

[itex]\sum _{\lambda = 1} ^n (m_{σλ} \ddot{\eta}_\eta + v_{σλ}\eta _λ),\ \ σ = 1,2,...,n[/itex]

Where m is the mass matrix:

[itex]m_{σλ} = \sum _i m_i (\frac{δx_i}{δq_σ})(\frac{δx_i}{δq_λ})[/itex]

Where the x_{i}are our cartesian coordinates, and the q's are our generalized coordinates.

And v is the "potential matrix" (I couldn't find its real name):

[itex]v_{σλ} = (\frac{δ^2 V}{δq_σ δq_λ})[/itex]

Where V is our potential.

We then guess that η is actually part of a complex number z, which would also be a solution if we had:

[itex]η = Re(z)[/itex]

And we guess a solution for z:

[itex]z = z_0 e^{iωt}[/itex]

And if we plug this into the first equation posted (for η, and differentiate, and divide by the exponent), we get an eigenvalue problem:

[itex]\sum _{\lambda = 1} ^n (v_{σλ} - ω^2 m_{σλ}) z_λ ^0,\ \ σ = 1,2,...,n[/itex]

Here, z^{s}is our eigenfunction corresponding to the eigenfrequency (value) ω_{s}^{2}.

3. The attempt at a solution

First of all, I set up my generalized coordinates. We are looking at small oscillations, so I define for each mass:

[itex]x_i = x_i ^0 + η_{ix}[/itex]

[itex]y_i = y_i ^0 + η_{iy}[/itex]

Where x^{0}is the equilibrium distance for that coordinate and each η is the small deviation from the equilibrium length. The index of the masses are as shown in the picture.

So, first of all, we know that because we have 6 degrees of freedom, we will get 6 normal modes (though some of them (and are, because the problem said so) can be degenerate).

There seem to be two ways to do this problem, but I'm pretty sure they're actually the same way in the end. The first is to set up the Lagrangian, then use the E-L equations to get 6 equations, the equations of motion. These equations could then be expressed in matrix notation as

[itex](\underline{v} - \underline{m}ω^2)\underline{z_s}[/itex]

Where v and m are the previous matrices, and z^{s}

On the other hand, we could just find the matrices m and v outright by their definitions above, rather than making the E-L equations match them.

Either way, from this matrix equation, we know solve for the determinant, and theoretically this should give us all the normal modes.

Using the above transformations, we quickly find the mass matrix to be diagonal, with m in each of the diagonal spots and zero otherwise:

[itex]m_{ij} = mδ_{ij}[/itex]

Either way I do it, I have to find V (to differentiate the Lagrangian, or to just solve for the matrix v). But when I try to find V, I seem to get nonlinear stuff no matter what I try. This nonlinear stuff can be differentiated for sure, but when I do, it's messy and you can't get any answers from it.

To calculate the potential, we can add up the potentials from each spring separately. For the spring between mass 1 and mass 2, I'll call the potential V_{12}, and so on.

So for spring 13, we call its equilibrium length a. Then, if it stretches and turns a little, its new length d is:

[itex]d = \sqrt{(a + η_{1x} - η_{3x})^2 + (η_{1y} - η_{3y})^2}[/itex]

and

[itex]V_{13} = \frac{k}{2}(d - a)^2[/itex]

And this is what I'm talking about...we don't even need to do the potentials for the other springs (though I did). The term with the square root just gets nastier as you differentiate it (twice) to find v. I assume I have to make some approximation to solve this, but I have no idea where... I tried assuming the springs don't turn while they're oscillating, and that makes the potentials V_{13}and V_{23}nice and simple (and solvable) but it still leaves V_{12}complicated and unsolvable.

Can anyone help me?? I'm at my wits end.

Thanks!

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Finding the Normal Modes of a 3 mass molecule connected by springs

**Physics Forums | Science Articles, Homework Help, Discussion**