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 КиБ) скачан 45 раз(а).
DuttilitaPressoflessione_R(forum).sm (215 КиБ) скачан 40 раз(а).
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 КиБ) скачан 31 раз(а).

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 КиБ) скачан 45 раз(а).

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