Welcome Guest! To enable all features please Login. New Registrations are disabled.

Notification

Icon
Error

Login


Options
Go to last post Go to first unread
Offline mdominiak  
#1 Posted : 07 March 2017 14:32:40(UTC)
mdominiak


Rank: Member

Groups: Registered
Joined: 12/12/2016(UTC)
Posts: 20
Poland
Location: Poznań

Hi,

In my previous topic I was struggling with getting properly built equations.

Right now, I have prepared document, where I am looking for two unknows, basing on set of two equations. Can anyone take a look and help me to solve it? (attachment)

BR.
Strain_stress - SMath forum.sm (32kb) downloaded 87 time(s).

Wanna join the discussion?! Login to your SMath Studio Forum forum account. New Registrations are disabled.

Offline Jean Giraud  
#2 Posted : 07 March 2017 19:21:24(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Few points:

1. The unit system is such a pest and never needed because maths are "unitless"
I have removed it in the attached.

2. If there are solution(s): it is "FindRoots"

3. For the system to solve: you must confirm in some way
that the system is solvable. I don't think it is solvable.
Our friend Radovan is "guru" in that stuff.

The "Quaternion" is added to prove "FindRoots" works when solvable.

Jean

Strain_stress - SMath forum Copy.sm (45kb) downloaded 68 time(s).
Offline mkraska  
#3 Posted : 07 March 2017 20:40:07(UTC)
mkraska


Rank: Advanced Member

Groups: Registered
Joined: 15/04/2012(UTC)
Posts: 1,980
Germany

Was thanked: 1120 time(s) in 718 post(s)
Originally Posted by: mdominiak Go to Quoted Post

Right now, I have prepared document, where I am looking for two unknows, basing on set of two equations. Can anyone take a look and help me to solve it? (attachment)


Why do you insist on defining the stress-strain behaviour using if statements?
And why do you split the domain of integrals? The numeric procedure doesn't need that.
Accuracy only benefits from splitting at steps (which you don't have).
So, unfortunately, the sheet looks much more complicated than actually required.

I already did simplify the functions and showed how to solve the system in this post. What was bad with that approach?

@ Jean: Units are an essential part of engineering documentation and language. It is a matter of efficient communication. Please respect that if you want to help.

I would not mind you opening a separate topic "War against units" and post there as much as you like. This could even lead to a fruitful discussion. But please resist the temptation to raise this issue where it is off-topic (e.g. in topics started by engineers). This would contribute to a better SNR in the forum.

Martin
Martin Kraska

Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx
Offline mdominiak  
#4 Posted : 07 March 2017 22:08:24(UTC)
mdominiak


Rank: Member

Groups: Registered
Joined: 12/12/2016(UTC)
Posts: 20
Poland
Location: Poznań

@mkraska: how can I define it otherwise? The domain of integrals is splitted as it is splitted in Eurocode, depending on strains. How to use numeric procedure here?
I would love to use steps, however I do not now how to introduce them (at least yet).
That sheet is one which I fully understand, it's neat and clean.

The other approach was not using correct statements (equilibrium states). I will look at your "simplification" more tommorow, however if it does not trouble you much, could you please change it in the newest smath document (attached here)?

@ Jean: Like always - I appreciate your effort and offered help.

BR.

edit: I have changed the newest version in accordance to @mkraska notes. See attachment.
By the way how can I add 'if' when stating functions (I was using your example)

Now basing on that, I will be looking for strains, when forces are given - using the same 'equations to be solved'...
Strain_stress - SMath forum_mkraska simplif.sm (68kb) downloaded 17 time(s).


edit2: I have also prepared document without units (however that does not please me), but maybe it will be better to get some results. I have an answer 'requested matrix doesn't exist', while in Jean file message was: x - not defined. Please see attachment.
Strain_stress - SMath forum_mkraska usun jednostki.sm (79kb) downloaded 22 time(s).

Edited by user 07 March 2017 23:31:41(UTC)  | Reason: Not specified

Offline mkraska  
#5 Posted : 08 March 2017 00:57:38(UTC)
mkraska


Rank: Advanced Member

Groups: Registered
Joined: 15/04/2012(UTC)
Posts: 1,980
Germany

Was thanked: 1120 time(s) in 718 post(s)
Bad chances for a solution

2017-03-07 22_55_15-SMath Studio 0.98.6179 - [Strain_stress%20-%20SMath%20forum_mkraska%20simplif.sm.png
Martin Kraska

Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx
Offline mdominiak  
#6 Posted : 08 March 2017 09:49:35(UTC)
mdominiak


Rank: Member

Groups: Registered
Joined: 12/12/2016(UTC)
Posts: 20
Poland
Location: Poznań

Could you explain me what have you done in that "check" and what was you looking at?

I attach general rules of the method. Maybe I have some mistake hidden in equations? Solution should exist.

Strain_stress - SMath forum.pdf (523kb) downloaded 26 time(s).

Offline omorr  
#7 Posted : 08 March 2017 14:33:54(UTC)
omorr


Rank: Administration

Groups: Registered, Advanced Member
Joined: 23/06/2009(UTC)
Posts: 1,740
Man
Serbia

Was thanked: 318 time(s) in 268 post(s)
Hello,

I might be wrong in this particular case but, in SMath if you have some equations with integrals, if statements, logical expression, absolute values etc. (some kind of discontinuity, or some other more "wierd" stuff) and want to solved them in a strictly numerical fashion ie. find the roots of the system, or to use some optimization procedure...you are doomed. Besides the possible non existing solutions (Martin pointed that out), SMath can hardly cope with this. I struggled with some similar problems and gave it up.

Regards,
Radovan

doomed.png

Edited by user 08 March 2017 14:36:33(UTC)  | Reason: Not specified

When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
Offline mdominiak  
#8 Posted : 08 March 2017 15:45:15(UTC)
mdominiak


Rank: Member

Groups: Registered
Joined: 12/12/2016(UTC)
Posts: 20
Poland
Location: Poznań

Originally Posted by: omorr Go to Quoted Post
...and want to solved them in a strictly numerical fashion ie. find the roots of the system, or to use some optimization procedure...you are doomed


How can I solve that otherwise? Is there any different way to get results (more or less accurate)?
Does different, similar software solve that? (f.e. Mathcad)

Edited by user 08 March 2017 15:46:59(UTC)  | Reason: Not specified

Offline PompelmoTell  
#9 Posted : 08 March 2017 17:25:55(UTC)
PompelmoTell


Rank: Advanced Member

Groups: Registered
Joined: 23/12/2011(UTC)
Posts: 318
Man
Italy
Location: italy

Was thanked: 109 time(s) in 93 post(s)
I had developed this SS sheet.
It is essential and missing calculating the resisting moment (but should not be a problem). Why do you want to simultaneously solve the two equations of equilibrium, given that it only takes one to determine the neutral axis?
It also lacks the management of the units I remember I gave it some thought.

sergio
Pressoflessione_R.sm (23kb) downloaded 51 time(s).
Offline mkraska  
#10 Posted : 08 March 2017 20:45:49(UTC)
mkraska


Rank: Advanced Member

Groups: Registered
Joined: 15/04/2012(UTC)
Posts: 1,980
Germany

Was thanked: 1120 time(s) in 718 post(s)
I tried to setup the procedure in an slightly alternative way and tried several ways to represent the stress-strain curves (cases, linterp, abs)

There may be errors in the equations, but this is just how far I got, just to understand the mechanical problem.

The equations fail to be solved using roots(). Also Davides functions can't handle the integration variable, which unfortunately is blamed to be unknown.

I guess, an option would be to provide analytical expressions for the stress resultant and it's centroid. For trapezoidal shapes that should be pretty easy.

Reinforced-beam.sm (312kb) downloaded 47 time(s).
2017-03-08 18_43_03-SMath Studio 0.98.6179 - [Reinforced-beam.sm_].png

Edited by user 08 March 2017 20:48:14(UTC)  | Reason: Not specified

Martin Kraska

Pre-configured portable distribution of SMath Studio: https://smath.com/wiki/SMath_with_Plugins.ashx
Offline mdominiak  
#11 Posted : 09 March 2017 11:56:22(UTC)
mdominiak


Rank: Member

Groups: Registered
Joined: 12/12/2016(UTC)
Posts: 20
Poland
Location: Poznań

Originally Posted by: PompelmoTell Go to Quoted Post
It is essential and missing calculating the resisting moment (but should not be a problem). Why do you want to simultaneously solve the two equations of equilibrium, given that it only takes one to determine the neutral axis?
It also lacks the management of the units I remember I gave it some thought.


At this moment I want to solve equation, where I have given design moment (M) and axial force (N). The results I want to get are strains in both concrete and steel.
You are right that I am not calculating resisting moment right now.

The next step will be to get all possible moments (M) and forces (N), where strains are given to prepare so called interaction diagram N/M. Result below.




The formulea is true only when Eps.1>=Eps.yd=fyd/Es and Eps.2>=Eps.yd=fyd/Es - we can use M.Rd (6.30).
Otherwise then in equation we have to use correct value of stress σs=Eps*E.s and σc - then we may use the same formulea (6.30), but instead of fcd and fyd we will have fcd(eps.c) and fyd(eps.s).

I don't know if it is clear for you, if not please do not hesitate to ask.
I am still looking for solution of that problem.
Offline PompelmoTell  
#12 Posted : 09 March 2017 13:44:18(UTC)
PompelmoTell


Rank: Advanced Member

Groups: Registered
Joined: 23/12/2011(UTC)
Posts: 318
Man
Italy
Location: italy

Was thanked: 109 time(s) in 93 post(s)
Conventional approaches for the formation of the interaction diagram provide for the allocation of the state strain from which are deduced the pairs of M-N values ​​using the equations of equilibrium. By repeating the procedure for a predetermined number of states deformation (significant and well-known) is obtained the diagram. Your approach seems also inadvisable from a numerical point of view as it would require the solution of the problem not linerare for a number of times not be defined a priori. About that last remark someone with more solid theoretical basis of my might prove me wrong.
sergio
Offline Jean Giraud  
#13 Posted : 09 March 2017 17:46:30(UTC)
Jean Giraud

Rank: Guest

Groups: Registered
Joined: 04/07/2015(UTC)
Posts: 6,868
Canada

Was thanked: 980 time(s) in 808 post(s)
Originally Posted by: mkraska Go to Quoted Post
Also Davides functions can't handle the integration variable, which unfortunately is blamed to be unknown.


Quite right Martin,

Mathcad 11 would not solve either because the integral is not a continuous
analytical function. Easy to cope via discretising the integral.
Finite differences is exceptionally valid for the remaining of the project.
There may be an alternate route via the "incremental companion". Not tried
because eventually the project must end to an export data table solution.

This project resources to many Smath great maths.
"Genfit Solve Martin" is attached so that nothing comes out of the blue.

In Project Reconciliation
please make sure about my comment:

My understanding is: 'x'= m^3/ 'steel J'

Jean

Genfit Solve Martin.sm (41kb) downloaded 50 time(s).

smath forum_Kr [Solve].sm (74kb) downloaded 57 time(s).


Users browsing this topic
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.