I turned a script that I had made in Scilab to use it in Smath.
The sheet serves to evaluate the moment diagram - curvature of a rectangular armed concrete section (useful then for further evaluations of the ductility of section ...)
The reduced version, which I am attaching here seems to provide correct results (or at least congruent with those obtained with Scilab ); However the slowness discourages in continuing with the completion of the procedure. Does anyone have any improvement idea?
I tried with a sheet without a unit of measurement; I get a 50% time reductions but it is always very especially wanting to get more dense discretion.
DuttilitaPressoflessione_R(forum).pdf (133.84 KiB) downloaded 365 time(s).
DuttilitaPressoflessione_R(forum).sm (215.47 KiB) downloaded 668 time(s).
sergio

Ciao Sergio. You can deal with this version of a numerical solver. Other improvements may be to avoid evaluating the integral for the cases x
DuttilitaPressoflessione_R(forum).sm (219.92 KiB) downloaded 648 time(s).

Best regards.
Alvaro.

Alvaro,
I don't think it's correct; I say this because I see that some users have downloaded your file and I wouldn't want to induce anyone wrong. It takes time to understand, but I believe it is a problem linked to the management of the units of measure with the solver you proposed purely numerical.
sergio

Hi Sergio. It could be because f is actually function of two variables, f(x,χ). This other setup ensure (I think) that at each loop the solver takes the new χ value.

DuttilitaPressoflessione_R(forum).sm (222.1 KiB) downloaded 654 time(s).

Best regards.
Alvaro.

Alvaro,
great now works well and quickly
Thanks!
sergio

Hi Sergio. You're a gentleman. You're right, I don't pay enough attention in my first post, and look only to the final result, but not it was stable or what.

Best regards.
Alvaro