Cannot solve a matrix equation with the Maxima plug-in

Cannot solve a matrix equation with the Maxima plug-in - Messaggi

#1 Pubblicato: 13/06/2023 06:19:51
EngMath

EngMath

0 apprezzamenti in 51 post.

Gruppo: User

Hi,

I'd like to solve the following equation involving matrices for unknown v and θ:



After inputting it in the format used by the Maxima plug-in I get the "no result" message:



Am I doing something wrong ?

Matrix equation.sm (8,59 KiB) scaricato 1413 volta(e).
#2 Pubblicato: 13/06/2023 08:42:27
Jean Giraud

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: User

Wrote

Am I doing something wrong ?


Your setup is not correct.

Solve.PNG

#3 Pubblicato: 13/06/2023 09:18:31
StvMath

StvMath

44 apprezzamenti e 1 disapprezzamenti in 76 post.

Gruppo: User

I think the problem is that the units don't match. Compare the units of the components of the two terms inside the brackets that multiply the [v; theta] vector.

#4 Pubblicato: 13/06/2023 09:29:58
Jean Giraud

Jean Giraud

983 apprezzamenti in 6.866 post.

Gruppo: User

Wrote

I think the problem is that the units don't match.
Compare the units of the components of the two terms inside the brackets
that multiply the [v; theta] vector.


FindRoots ignores Units. Code the image first.

Solve FindRoots Complex Forum.sm (11,09 KiB) scaricato 1485 volta(e).
#5 Pubblicato: 13/06/2023 09:30:57
EngMath

EngMath

0 apprezzamenti in 51 post.

Gruppo: User

Wrote


Your setup is not correct.

Solve.PNG



Thank you, I didn't know that I should use the Nonlinear Solvers plug-in instead of the Maxima one.

It works now but the result is surprising. It should be something like v = - 0.021020 for the real part with this input data.

Wrote


I think the problem is that the units don't match. Compare the units of the components of the two terms inside the brackets that multiply the [v; theta] vector.



I also suspected the units initially but it didn't work without them either.
#6 Pubblicato: 13/06/2023 11:17:16
StvMath

StvMath

44 apprezzamenti e 1 disapprezzamenti in 76 post.

Gruppo: User



I also suspected the units initially but it didn't work without them either.



In order to be able to add the two terms they must have the same units. However, omega^2*L doesn't have the same units as EI/L^3. The omega^2*... term seems to be missing a linear mass density (kg/m).
The fact that you might be able to get a numerical result by simply removing the units is irrelevant.
#7 Pubblicato: 13/06/2023 14:55:41
EngMath

EngMath

0 apprezzamenti in 51 post.

Gruppo: User

Wrote


In order to be able to add the two terms they must have the same units. However, omega^2*L doesn't have the same units as EI/L^3. The omega^2*... term seems to be missing a linear mass density (kg/m).
The fact that you might be able to get a numerical result by simply removing the units is irrelevant.



Right, I was wondering why density/mass is not included in that equation. Especially since it's included in the equation for tension (the one here is for bending): https://code-aster.org/V2/doc/v13/en/man_v/v2/v2.06.101.pdf

Maybe it's a typo and the density should be added before the first omega. Or omega means something different than I assumed (angular frequency).

Edit: I tried calculating mass per unit length based on the density and cross-section area provided in the document and then multiplying it by omega^2 but the result is v = 0,0175 which is still far from the expected value.
#8 Pubblicato: 13/06/2023 16:34:30
Martin Kraska

Martin Kraska

1.361 apprezzamenti in 2.275 post.

Gruppo: Moderator

Besides correction of the unit mismatch, Solve seems to want the matrix equation in homogeneous form.

Matrix equation_Kr.png
Matrix equation_Kr.sm (14,1 KiB) scaricato 1419 volta(e).
Technische Mechanik mit SMath Studio: https://link.springer.com/book/10.1007/978-3-658-50592-9
#9 Pubblicato: 13/06/2023 18:25:48
Alvaro

Alvaro

1.139 apprezzamenti in 1.805 post.

Gruppo: User

Hi. Another way.

Clipboard01.png

Best regards.
Alvaro.
  • Nuovi post
  • Nessun nuovo post