Rank: Advanced Member Groups: Registered
Joined: 23/12/2011(UTC) Posts: 319 Location: italy Was thanked: 109 time(s) in 93 post(s)
|
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 (134kb) downloaded 23 time(s). DuttilitaPressoflessione_R(forum).sm (216kb) downloaded 18 time(s).sergio
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,359 Was thanked: 819 time(s) in 518 post(s)
|
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 <h in f and M using if, and even to use your own numerical integrator, but this version seems fast enough. DuttilitaPressoflessione_R(forum).sm (220kb) downloaded 13 time(s).Best regards. Alvaro.
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 23/12/2011(UTC) Posts: 319 Location: italy Was thanked: 109 time(s) in 93 post(s)
|
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 Edited by user 07 November 2021 00:58:56(UTC)
| Reason: Not specified
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,359 Was thanked: 819 time(s) in 518 post(s)
|
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 (223kb) downloaded 24 time(s).Best regards. Alvaro.
|
2 users thanked Razonar for this useful post.
|
on 07/11/2021(UTC), on 08/11/2021(UTC)
|
|
Rank: Advanced Member Groups: Registered
Joined: 23/12/2011(UTC) Posts: 319 Location: italy Was thanked: 109 time(s) in 93 post(s)
|
Alvaro, great now works well and quickly Thanks! sergio Edited by user 07 November 2021 19:38:34(UTC)
| Reason: Not specified
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 28/08/2014(UTC) Posts: 1,359 Was thanked: 819 time(s) in 518 post(s)
|
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
|
1 user thanked Razonar for this useful post.
|
|
|
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.