Trying to solve a very complex system of equations containing 3 sets of equations. However it keeps appearing to get stuck at 46/47 processed. Not sure why this is happening, I have a very high powered PC since I have to run FEA with millions of elements. So an equation like this should be a breeze. Here is what the sheet looks like: https://imgur.com/a/r0JVOje Attached is the file if you want to dig around and try to get it to spit out the three values I need. Any help is GREATLY appreciated.Stresses.sm (36 КиБ) скачан 35 раз(а).

EDIT: found an error, didn't change sigma.1h.1 to sigma.1. Problem still persists even after the change.

I seem to be able to say that there are several errors is formal in the use of Smath and the Maxima plugin but also of the engineering problem. I report you:
- The three equations do not require a solutor because each contains a single unknown
- The way of recalling the "Solve" function is not appropriate (the three equations must be inserted in a vector or system. Even the unknowns must be inserted in a vector)
- The first two equations were written correctly with = the third instead with :=
- The second equation is not consistent as it is dimensionally we have a stress (MPa) obtained as the sum of elements that are not stress (1/MPa)

Consult Martin Kraska's interactive guide for the correct use of Solve

sergio

Wrote

I seem to be able to say that there are several errors is formal in the use of Smath and the Maxima plugin but also of the engineering problem. I report you:
- The three equations do not require a solutor because each contains a single unknown
- The way of recalling the "Solve" function is not appropriate (the three equations must be inserted in a vector or system. Even the unknowns must be inserted in a vector)
- The first two equations were written correctly with = the third instead with :=
- The second equation is not consistent as it is dimensionally we have a stress (MPa) obtained as the sum of elements that are not stress (1/MPa)

Consult Martin Kraska's interactive guide for the correct use of Solve

sergio

Thanks for the reply!

- thanks for catching that! I wrote the last part of the equation incorrectly, it should be everything multiplied by E.I instead of divided by.
- ahhh, I see the equations looked a lot more complicated than they were. I was able to get it to solve by setting one of the equations to a boolean equal and now it gives me the values I am looking for.

Thanks for your help solving this!

EDIT: Nevermind, error still persists. Cleaning out the old equations that were just trial got rid of the solution and it keeps saying the values are "undefined". Tried tossing it back into the solver with the corrected formatting, still won't solve it. Stresses.sm (32 КиБ) скачан 25 раз(а). Here is the updated sheet.

So how do I get this to solve for the three values?

I should note, in Maxima debug, when trying to evaluate it Maxima returns "no data"

Wrote

So how do I get this to solve for the three values?

I use linux, so I can't use Maxima plugin.
However, your functions can be solved with al_nleqsolve().
I don't know if it is correct in terms of engineering.
I just solved it only regards to an SMath Solver sample.

Regards

Stresses(1).sm (49 КиБ) скачан 28 раз(а).

Wrote

Wrote

So how do I get this to solve for the three values?

I use linux, so I can't use Maxima plugin.
However, your functions can be solved with al_nleqsolve().
I don't know if it is correct in terms of engineering.
I just solved it only regards to an SMath Solver sample.

Regards

Stresses(1).sm (49 КиБ) скачан 28 раз(а).

Thanks for the reply and attempting to solve this! The values are extremely low, I am expecting them to be around 50-200 MPa. Something must be wrong here, even when I change the wall thickness, the values barely change. The wall thickness is the major determinant (t.1, t.2) in these equations for stresses. I think it might be useful to drop the link to where the equations came from: https://ocw.mit.edu/courses/16-512-rocket-propulsion-fall-2005/357a8f245e0020bcfa43530ba6ed1283_lecture_9.pdf In there example they routinely got around 200 MPa but were using a thinner wall, when I tossed their thicknesses in, the values were still absurdly small.

The problems with Solve() seem to be due to translation problems of names like sigma_1.1, where the period as part of the subscript is not escaped.

Unfortunately, this error is not catched, so that the Maxima process gets stuck (restart using either MaximaControl("restart") or using Insert> Maxima> Settings> Save and restart)

With modified names, a solution can be obtained.

Stresses_Kr.sm (31 КиБ) скачан 26 раз(а).

Wrote

...
Thanks for the reply and attempting to solve this! The values are extremely low, I am expecting them to be around 50-200 MPa. Something must be wrong here, ...

Hi. From your source, the system is with equations (7), (10) and (13) and "From here, we either solve for σ2h, σ2c, σ1 (given geometry, Pg, Pl, Pq) or, for design, go the reverse route and decide on the geometry for assigned stresses."

I see equations (7) and (10) but for equation (13) I guess that there are something wrong there in the equations system to be solved.

Best regards.
Alvaro.

Actually, FindRoot() works with your function too.
It only requires to disable these variables.



After that, function worked as expected. Unknowns() is important btw.
Solvers find results according to order of Unknowns() output order.
Otherwise one would require to check which result is for which variable.

Regards