SMath Studio Forum
»
SMath Studio
»
Questions
»
Solving a set of equations with function implemented in another
  Rank: Member Groups: Registered
Joined: 12/12/2016(UTC) Posts: 20  Location: Poznań
|
Hello, I have an equation, which depends on two other function, which depends on one variable. How can I solve it? Look at the attached file.  smath forum.sm (97kb) downloaded 26 time(s).In general, I have set of two equations: Fs=Nc(x) and Fs*(d-x_compression)+Nc(x)*x_b*x_compression=M_Ed Now, I implemented all funcitons as function of x (x stays for Ec - strain in concrete), but I can easily manipulate compression zone (x_compression) if necessary, as it is also not given as "constant". At the beginning I define dimensions (b,h,a1) and material (fck, gamma_c, alfa_cc). I also assume compression zone and strain in concrete and I have bending moment M_Ed. Rest (and also x_compression and E_c) are editable/may change and depend on each other. Could you please tell me what's the easiest way to solve those two equations? Unfortunatelly, when I am trying to solve(Nc(x)=Fs,x)=0, SMath respond is that x has not been defined. Any ideas? BR.
|
|
|
|
|
|
Rank: Member Groups: Registered
Joined: 12/12/2016(UTC) Posts: 20 Location: Poznań
|
For the needs of solving, you may assume that: - x_compression is given and constant (probably best for now)
- x_compression is variable (but I should probably do it another way)
BR
|
|
|
|
|
|
 Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866  Was thanked: 981 time(s) in 809 post(s)
|
As it looks, units is the problem . Some remarks at the end. Jean smath forum.sm (105kb) downloaded 30 time(s).
|
|
|
|
|
|
Rank: Member Groups: Registered
Joined: 12/12/2016(UTC) Posts: 20 Location: Poznań
|
Jean, thanks for your effort. You have written: "nothing to solve as you have two numbers" - not exactly, once you will look into the sum of steel and concrete you may see there Nc(Ec) - Ec is variable. "f(x) must be a function" - isn't f(x)=123 a function? Still, I don't have any results. I would like to know at what argument Nc(x) is equal to Fs (Nc(x)=16.37, as then I may estimate strains in concrete. I don't know maybe you see a result? I have attached photo.  In the meantime I have prepared another set of equation for which I would need to have similiar solution as for that one above. Could you take one more look at it? BR.
|
|
|
|
|
|
 Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,984  Was thanked: 1124 time(s) in 721 post(s)
|
I have simplified the function definitions and the integrals and fixed the units. Don't know if the result is correct/matches the expectation.  smath%20forum_Kr.sm (98kb) downloaded 33 time(s). |
|
|
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 981 time(s) in 809 post(s)
|
Originally Posted by: mdominiak  Could you take one more look at it? Yes, and the answer is that there is no automatic solution. Why ? Because the base "concrete(x)" is not a function, only a discontinuous function, incompatible to the solver(s). Done manually at the end of the work sheet. Maybe "Dichotomy" is suitable ? Jean smath forum.sm (113kb) downloaded 16 time(s).
|
|
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 981 time(s) in 809 post(s)
|
... "Dichotomy" solves nicely ! Solves for 'x' in "concrete(x)" for any value(s) in M Ed(x) Is it done as expected ? Your stuff is not my cup of tea ! Jean smath forum.sm (117kb) downloaded 34 time(s).
|
|
|
|
|
|
Rank: Member Groups: Registered
Joined: 12/12/2016(UTC) Posts: 20 Location: Poznań
|
Jean, mkraska,
many thanks for your effort.
I was pretty busy last few days and this is something I am working on aside of main job, let's say.
I will look at it by the end of the week and let you know if your solution is working.
BR.
|
|
|
|
|
|
Rank: Member Groups: Registered
Joined: 12/12/2016(UTC) Posts: 20 Location: Poznań
|
Can we use Newton-Raphson method instead of Dichotomy for example?
|
|
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 981 time(s) in 809 post(s)
|
Originally Posted by: mdominiak Can we use Newton-Raphson method instead of Dichotomy for example? Answer is : NO MED(x) is a constant , Newton-Raphson supports only functions Dichotomy supports MED(x) defined as a constant, i.e: it supports "fake function" as well as the 'solve block' In your document, 'solve' fails because of the unit system. Jean
|
|
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 981 time(s) in 809 post(s)
|
... rootsecant fails as well.
|
|
|
|
|
|
Rank: Member Groups: Registered
Joined: 12/12/2016(UTC) Posts: 20 Location: Poznań
|
Originally Posted by: Jean Giraud Originally Posted by: mdominiak
MED(x) is a constant , Newton-Raphson supports only function
Med(x) in future may be function, not only constant. Will dichtomy be enough then to solve it? From what I understand N-R may be applicable then, also? Could you give me mathematical definition of dichotomy (it would be best if you can demonstrate it on that particular example)? Then, you have stated something in pink background at the end of document. Was it check? You have defined values Epsilon, a and b - and I would like to know why.
|
|
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866 Was thanked: 981 time(s) in 809 post(s)
|
The RootSecant & NewtonRaphson both need true functions. Differently, "Dichotomy" takes a constant defined as a fake function. In dichotomy: epsilon is the test for the convergence towards the solution, solution desired at "epsilon accuracy". The while loop keeps refining the solution at each step... more in green. Jean smath forum.sm (133kb) downloaded 31 time(s).
|
|
|
|
|
|
Rank: Member Groups: Registered
Joined: 12/12/2016(UTC) Posts: 20 Location: Poznań
|
Jean, thanks again. I have worked with this tool a bit today and yesterday and I made some changes. Right now I cannot get any results. I have added some picture and comments inside document. Also I have added values which shall be results for given data. Can you look at them? I am looking for: strains at both edges, value of x_compression. I have changed strain-stress function to be easier to verify results (once it will be ok, I will just change it into previous one). BR. Zadanie - dane sily - oblicz odksztalcenia.sm (767kb) downloaded 18 time(s). - problematic Smath document (given forces, and looking for strains and compression area) Zadanie - dane odksztalcenia - oblicz x M i N.sm (53kb) downloaded 17 time(s). - the same problem, however solved with given strains (from them I've got compression area) and I was looking for forces. Let me know in case of any questions to the idea itself. Edited by user 06 March 2017 21:28:51(UTC)
| Reason: spelling
|
|
|
|
|
|
Rank: Member Groups: Registered
Joined: 12/12/2016(UTC) Posts: 20 Location: Poznań
|
I have just one quick amendment, namely I am able to get proper strains in steel (εs(x)), as I know first value of force in bars N_s(x). Then: • σs(x)= εs(x)*E • N_s(x)/As1= εs(x)*E • εs(x)= N_s(x)/(As1*E) Therefore value “y” from SMath document is known (however it depends on x). Still, I am looking for εc(x) (once I’ll have that value I will be able to get proper x_compression form Tales) or x_compression (once I’ll have that value I will be able to get proper εc(x)). Note: values of strains (ε) and x_compression may vary, however for given set of N/M and geometric characteristics shall be equal to εs(x) = 0.3% and εc(x) = 0.15%. Edit: See attached SMath file: I have added example.
Once I’ll finish it (with your help ) I will be interested in having all possible solutions, plotted on graph Moments/Normal forces in 2D. Zadanie - dane sily - oblicz odksztalcenia.sm (771kb) downloaded 26 time(s).Edited by user 07 March 2017 09:42:53(UTC)
| Reason: SMath file added.
|
|
|
|
|
|
SMath Studio Forum
»
SMath Studio
»
Questions
»
Solving a set of equations with function implemented in another
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.
