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 genf  
#1 Posted : 27 February 2021 23:04:03(UTC)
genf


Rank: Advanced Member

Groups: Registered
Joined: 10/06/2013(UTC)
Posts: 50
Man
Russian Federation

Was thanked: 6 time(s) in 6 post(s)
Добрый день.
Помогите (или подскажите) пожалуйста с оптимизацией. Довольно долгий расчет (около 30 сек). Если уменьшить параметр m, то время расчета уменьшается, но хотелось бы как то оптимизировать. Gidrodinamika.sm (106kb) downloaded 7 time(s).

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

Offline Ber7  
#2 Posted : 28 February 2021 01:33:13(UTC)
Ber7


Rank: Advanced Member

Groups: Registered
Joined: 15/07/2010(UTC)
Posts: 437
Man
Israel
Location: Beer-Sheva

Was thanked: 520 time(s) in 288 post(s)
If you shorten the code, then the calculation is performed
an order of magnitude faster

Gidrodinamika1.sm (112kb) downloaded 16 time(s).
thanks 1 user thanked Ber7 for this useful post.
on 28/02/2021(UTC)
Offline genf  
#3 Posted : 28 February 2021 22:05:02(UTC)
genf


Rank: Advanced Member

Groups: Registered
Joined: 10/06/2013(UTC)
Posts: 50
Man
Russian Federation

Was thanked: 6 time(s) in 6 post(s)
Финальный расчет во вложении.
Большое спасибо Ber7 за помощь в оптимизацииClap
Opredelenie gidrodinamicheskogo davlenija_v001.sm (129kb) downloaded 11 time(s).
Offline genf  
#4 Posted : 01 March 2021 11:18:30(UTC)
genf


Rank: Advanced Member

Groups: Registered
Joined: 10/06/2013(UTC)
Posts: 50
Man
Russian Federation

Was thanked: 6 time(s) in 6 post(s)
Эх, по поводу "финальной" это я погорячился.
Как оказалось есть недочет в 1 формуле.
Formula.png

Как это реализовать можно?
Т.е. именно в сумму входят только нечетные номера.

Edited by user 01 March 2021 11:35:21(UTC)  | Reason: Not specified

Offline Ber7  
#5 Posted : 01 March 2021 13:16:06(UTC)
Ber7


Rank: Advanced Member

Groups: Registered
Joined: 15/07/2010(UTC)
Posts: 437
Man
Israel
Location: Beer-Sheva

Was thanked: 520 time(s) in 288 post(s)
I don't understand this formula, but Try to substitute 2n -1 instead of n

2n-1.sm (135kb) downloaded 11 time(s).
thanks 1 user thanked Ber7 for this useful post.
on 01/03/2021(UTC)
Offline genf  
#6 Posted : 01 March 2021 13:40:11(UTC)
genf


Rank: Advanced Member

Groups: Registered
Joined: 10/06/2013(UTC)
Posts: 50
Man
Russian Federation

Was thanked: 6 time(s) in 6 post(s)
Originally Posted by: Ber7 Go to Quoted Post
I don't understand this formula, but Try to substitute 2n -1 instead of n

2n-1.sm (135kb) downloaded 11 time(s).


Да это то, что нужно. Спасибо.
Opredelenie gidrodinamicheskogo davlenija_v001.sm (188kb) downloaded 7 time(s).
Offline Jean Giraud  
#7 Posted : 01 March 2021 17:37:12(UTC)
Jean Giraud

Rank: Guest

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

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: Ber7 Go to Quoted Post
I don't understand this formula

I don't understand much in there either.
We can get the shape of the upper right image.

2n-1.sm (143kb) downloaded 6 time(s).
Offline Jean Giraud  
#8 Posted : 01 March 2021 20:32:53(UTC)
Jean Giraud

Rank: Guest

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

Was thanked: 981 time(s) in 809 post(s)
Interesting shape but can't discretize ???

2n-1 [Polylog].sm (9kb) downloaded 3 time(s).
Offline Jean Giraud  
#9 Posted : 01 March 2021 22:51:45(UTC)
Jean Giraud

Rank: Guest

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

Was thanked: 981 time(s) in 809 post(s)
What's your comment on that ?

Unknown.PNG
thanks 1 user thanked Jean Giraud for this useful post.
on 01/03/2021(UTC)
Offline genf  
#10 Posted : 01 March 2021 23:55:59(UTC)
genf


Rank: Advanced Member

Groups: Registered
Joined: 10/06/2013(UTC)
Posts: 50
Man
Russian Federation

Was thanked: 6 time(s) in 6 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post
What's your comment on that ?



Визуально очень похоже... нужно потестировать по ряду значений ... есть графики в книге Бирбраера (рис. 9.9). Я его оцифровывал с помощью GetData Graph Digitizer (примерно конечно). Birbraer_ris.9.9.xlsx (104kb) downloaded 5 time(s).
Именно эта формула и графики одно и тоже.

В Smath у меня эта эпюра ломается внизу ... пока не разобрался почему.

Edited by user 01 March 2021 23:58:39(UTC)  | Reason: Not specified

Offline Jean Giraud  
#11 Posted : 02 March 2021 02:51:55(UTC)
Jean Giraud

Rank: Guest

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

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: genf Go to Quoted Post
GetData Graph Digitizer

We can do it in Smath providing the image is cleaned manually
and the coordinates [Xmax, Xmin ...Ymax, Ymin] are given or
otherwise guesstimated. Freeware Tracer can do as well but
again cleaned image and coordinates.
Your approximation sounds OK ... many more can be envisaged.
My little Russian is dead for over 50 years .
Cheers ... Jean

thanks 1 user thanked Jean Giraud for this useful post.
on 02/03/2021(UTC)
Offline Jean Giraud  
#12 Posted : 02 March 2021 04:57:15(UTC)
Jean Giraud

Rank: Guest

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

Was thanked: 981 time(s) in 809 post(s)
... smoothed data from Tracer.

Page15 Tracer.sm (19kb) downloaded 4 time(s).
Offline Jean Giraud  
#13 Posted : 02 March 2021 08:38:47(UTC)
Jean Giraud

Rank: Guest

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

Was thanked: 981 time(s) in 809 post(s)
Originally Posted by: Jean Giraud Go to Quoted Post
... smoothed data from Tracer.

... it can be improved a bit. But, what a dirty image to digitize !

Page15 Tracer.sm (72kb) downloaded 7 time(s).
Offline Jean Giraud  
#14 Posted : 02 March 2021 20:23:15(UTC)
Jean Giraud

Rank: Guest

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

Was thanked: 981 time(s) in 809 post(s)
This document resumes and compares your fitting model
vs the digitized image, from improved session.
You may attempt to adapt your model to the digitized data.
Maybe there is a Thiele fit to your blue/red data points.
I will try that soon ... Jean.

2n-1 [Jean].sm (189kb) downloaded 7 time(s).
thanks 1 user thanked Jean Giraud for this useful post.
on 02/03/2021(UTC)
Offline genf  
#15 Posted : 02 March 2021 21:20:35(UTC)
genf


Rank: Advanced Member

Groups: Registered
Joined: 10/06/2013(UTC)
Posts: 50
Man
Russian Federation

Was thanked: 6 time(s) in 6 post(s)
В книге Бирбраера рис. 9.8 несет в себе только информацию касаемо куда какая нагрузка прикладывается и форма приведена условно, т.е. схема. Сам рис. 9.8 не претендует на оцифровку - это не график, а картинка.
А вот рис. 9.9 - это набор графиков уже и именно их я и оцифровывал еще когда не разобрался с формулой 9.51.
Т.е. именно этот набор графиков определяет форму эпюры гидродинамического давления.
Offline Jean Giraud  
#16 Posted : 02 March 2021 22:58:13(UTC)
Jean Giraud

Rank: Guest

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

Was thanked: 981 time(s) in 809 post(s)
In Birbraer's book, fig. 9.8 carries only information with regards to where what load is applied and the form is given conditionally, i.e. scheme. The rice itself. 9.8 does not pretend to be digitized - it is not a graph, but a picture.
And here is rice. 9.9 is a set of graphs already, and it was them that I digitized even when I had not figured out the formula 9.51.
Thus, it is this set of plots that determines the shape of the hydrodynamic pressure plot.
...................
So, my understanding is that the formula 9.51 is the one that determines,
the hydrodynamic pressure. It looks damned close parent to Eiffel Tower shape.
Then, it is not needed to be discretized, simply plot scalar.
However, I don't think you can plot opposite reverse, wish I'm wrong.
Very interesting project, thanks ... Jean

2n-1 [Jean].sm (244kb) downloaded 14 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.