Cannot solve a matrix equation with the Maxima plug-in

Cannot solve a matrix equation with the Maxima plug-in - Сообщения

#1 Опубликовано: 13.06.2023 06:19:51
EngMath

EngMath

0 сообщений из 51 понравились пользователям.

Группа: 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 КиБ) скачан 35 раз(а).
#2 Опубликовано: 13.06.2023 08:42:27
Jean Giraud

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

Wrote

Am I doing something wrong ?


Your setup is not correct.

Solve.PNG

#3 Опубликовано: 13.06.2023 09:18:31
StvMath

StvMath

35 сообщений из 56 понравились пользователям.

Группа: 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 Опубликовано: 13.06.2023 09:29:58
Jean Giraud

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: 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 КиБ) скачан 27 раз(а).
#5 Опубликовано: 13.06.2023 09:30:57
EngMath

EngMath

0 сообщений из 51 понравились пользователям.

Группа: 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 Опубликовано: 13.06.2023 11:17:16
StvMath

StvMath

35 сообщений из 56 понравились пользователям.

Группа: 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 Опубликовано: 13.06.2023 14:55:41
EngMath

EngMath

0 сообщений из 51 понравились пользователям.

Группа: 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 Опубликовано: 13.06.2023 16:34:30
Martin Kraska

Martin Kraska

1222 сообщений из 2150 понравились пользователям.

Группа: 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 КиБ) скачан 22 раз(а).
Martin Kraska Pre-configured portable distribution of SMath Studio: https://en.smath.info/wiki/SMath%20with%20Plugins.ashx
#9 Опубликовано: 13.06.2023 18:25:48
Alvaro Diaz Falconi

Alvaro Diaz Falconi

992 сообщений из 1674 понравились пользователям.

Группа: User

Hi. Another way.

Clipboard01.png

Best regards.
Alvaro.
  • Новые сообщения Новые сообщения
  • Нет новых сообщений Нет новых сообщений