Welcome Guest! To enable all features please Login. New Registrations are disabled.

Notification

Icon
Error

Login


Options
Go to last post Go to first unread
Offline EngMath  
#1 Posted : 13 June 2023 12:19:51(UTC)
EngMath


Rank: Advanced Member

Groups: Registered
Joined: 18/10/2020(UTC)
Posts: 50
Poland

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 (9kb) downloaded 14 time(s).

Wanna join the discussion?! Login to your SMath Studio Forum forum account. New Registrations are disabled.

Offline Jean Giraud  
#2 Posted : 13 June 2023 14:42:27(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: EngMath Go to Quoted Post
Am I doing something wrong ?

Your setup is not correct.

Solve.PNG

Offline StvMath  
#3 Posted : 13 June 2023 15:18:31(UTC)
StvMath


Rank: Advanced Member

Groups: Registered
Joined: 26/07/2022(UTC)
Posts: 35
United Kingdom
Location: UK

Was thanked: 22 time(s) in 19 post(s)
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.

Offline Jean Giraud  
#4 Posted : 13 June 2023 15:29:58(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,866
Canada

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: StvMath Go to Quoted Post
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 (12kb) downloaded 3 time(s).
Offline EngMath  
#5 Posted : 13 June 2023 15:30:57(UTC)
EngMath


Rank: Advanced Member

Groups: Registered
Joined: 18/10/2020(UTC)
Posts: 50
Poland

Originally Posted by: Jean Giraud Go to Quoted Post

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.

Originally Posted by: StvMath Go to Quoted Post

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.
Offline StvMath  
#6 Posted : 13 June 2023 17:17:16(UTC)
StvMath


Rank: Advanced Member

Groups: Registered
Joined: 26/07/2022(UTC)
Posts: 35
United Kingdom
Location: UK

Was thanked: 22 time(s) in 19 post(s)


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.

Edited by user 13 June 2023 17:20:06(UTC)  | Reason: Not specified

Offline EngMath  
#7 Posted : 13 June 2023 20:55:41(UTC)
EngMath


Rank: Advanced Member

Groups: Registered
Joined: 18/10/2020(UTC)
Posts: 50
Poland

Originally Posted by: StvMath Go to Quoted Post

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/V...n/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.

Edited by user 13 June 2023 21:21:32(UTC)  | Reason: Not specified

Offline mkraska  
#8 Posted : 13 June 2023 22:34:30(UTC)
mkraska


Rank: Advanced Member

Groups: Registered
Joined: 15/04/2012(UTC)
Posts: 1,986
Germany

Was thanked: 1124 time(s) in 721 post(s)
Besides correction of the unit mismatch, Solve seems to want the matrix equation in homogeneous form.

Matrix equation_Kr.png
Matrix equation_Kr.sm (15kb) downloaded 5 time(s).
Martin Kraska

Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx
Offline Razonar  
#9 Posted : 14 June 2023 00:25:48(UTC)
Razonar


Rank: Advanced Member

Groups: Registered
Joined: 28/08/2014(UTC)
Posts: 1,356
Uruguay

Was thanked: 815 time(s) in 516 post(s)
Hi. Another way.

Clipboard01.png

Best regards.
Alvaro.
Users browsing this topic
Guest
Forum Jump  
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.