Solving tricky functions - Hertzian contact stress - Сообщения
#1 Опубликовано: 20.10.2022 11:38:20
Good day everyone,
I am building a Hertzian stress calculator and got stuck at the step where I calculate the k value for a particular B/A reationship.
Any help will be appriciated
contact stress calc.sm
I am building a Hertzian stress calculator and got stuck at the step where I calculate the k value for a particular B/A reationship.
Any help will be appriciated
contact stress calc.sm
#2 Опубликовано: 20.10.2022 11:48:46
We also know than K will always be in the range (0,1)
The equations presented above is obtained from the book Advanced Mechanics of Materials by Boresi Et al
The equations presented above is obtained from the book Advanced Mechanics of Materials by Boresi Et al
#3 Опубликовано: 20.10.2022 17:34:39
In case of these functions, al_nleqsolve() can be the answer.
No other solver is capable other than al_nleqsolve().
k' should be under line, otherwise not solveable.
Here is your function solved.
Regards
contact stress calc_1.sm
No other solver is capable other than al_nleqsolve().
k' should be under line, otherwise not solveable.
Here is your function solved.
Regards
contact stress calc_1.sm
3 пользователям понравился этот пост
NDTM Amarasekera 21.10.2022 09:20:00, sergio 21.10.2022 05:43:00, Alvaro Diaz Falconi 24.10.2022 00:51:00
#4 Опубликовано: 21.10.2022 07:52:33
I'm super impressed !!!
#5 Опубликовано: 23.10.2022 16:46:37
Your original document is circular and indigestible.
From un-sourced [K,E] and arbitrary datum ... solve.
This sample is not appropriate.
As you have said, it is only solvable when K and E already given.
roots() can't solve without K and E given or calculated separately.
And you can't calculate K and E while not knowing what k is.
Your sample is not what OP wants.
He means to solve that complicated Q function.
Without any parameter or variable of K and E are pre-given.
2 пользователям понравился этот пост
#6 Опубликовано: 24.10.2022 09:26:26
With elliptic integrals from GSL and Cephes library.
Russia ☭ forever, Viacheslav N. Mezentsev
3 пользователям понравился этот пост
Alvaro Diaz Falconi 26.10.2022 14:25:00, Oscar Campo 24.10.2022 11:14:00, sergio 25.10.2022 14:02:00
#7 Опубликовано: 24.10.2022 20:29:42
Good day everyone,
I am building a Hertzian stress calculator and got stuck at the step where I calculate the k value for a particular B/A reationship.
Any help will be appriciated
contact stress calc.sm
Solving by the old bisection method:
contact stress calc_by_bisection.sm
2 пользователям понравился этот пост
#9 Опубликовано: 25.10.2022 12:08:25
Solving with a custom algorithm.
Very elegant!
#10 Опубликовано: 28.10.2022 08:38:03
... slightly refreshed/augmented
Absent of comments from Originator !
Cheers ... Jean.
Solve Circular Elliptic Datum Native Copy.sm
Absent of comments from Originator !
Cheers ... Jean.
Solve Circular Elliptic Datum Native Copy.sm
-
Новые сообщения
-
Нет новых сообщений