 Rank: Newbie
Groups: Registered
Joined: 20/10/2022(UTC) Posts: 2 
|
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 (23kb) downloaded 27 time(s).
|
|
|
|
|
|
Rank: Newbie
Groups: Registered
Joined: 20/10/2022(UTC) Posts: 2
|
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
|
|
|
|
|
|
 Rank: Advanced Member Groups: Registered
Joined: 23/07/2013(UTC) Posts: 1,128  Was thanked: 509 time(s) in 339 post(s)
|
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 (29kb) downloaded 28 time(s). 
|
 3 users thanked overlord for this useful post.
|
on 21/10/2022(UTC), on 21/10/2022(UTC), on 24/10/2022(UTC)
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 13/09/2011(UTC)
Posts: 63
Location: Kolkata
Was thanked: 9 time(s) in 9 post(s)
|
|
|
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 23/07/2013(UTC) Posts: 1,128 Was thanked: 509 time(s) in 339 post(s)
|
Originally Posted by: Jean Giraud  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. Edited by user 24 October 2022 02:02:33(UTC)
| Reason: Not specified
|
2 users thanked overlord for this useful post.
|
on 24/10/2022(UTC), on 24/10/2022(UTC)
|
|
|
 Rank: Advanced Member Groups: Registered, Advanced Member Joined: 10/11/2010(UTC) Posts: 1,494   Was thanked: 1274 time(s) in 745 post(s)
|
With elliptic integrals from GSL and Cephes library.   Edited by user 24 October 2022 15:54:50(UTC)
| Reason: Not specified |
Russia ☭ forever
Viacheslav N. Mezentsev |
3 users thanked uni for this useful post.
|
on 24/10/2022(UTC), on 25/10/2022(UTC), on 26/10/2022(UTC)
|
|
|
 Rank: Advanced Member Groups: Registered
Joined: 10/12/2009(UTC) Posts: 248  Location: Cali, Colombia Was thanked: 87 time(s) in 66 post(s)
|
Originally Posted by: ryno 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 (23kb) downloaded 27 time(s). Solving by the old bisection method:  contact stress calc_by_bisection.sm (27kb) downloaded 6 time(s).
|
2 users thanked oscampo for this useful post.
|
on 25/10/2022(UTC), on 26/10/2022(UTC)
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 23/07/2013(UTC) Posts: 1,128 Was thanked: 509 time(s) in 339 post(s)
|
|
3 users thanked overlord for this useful post.
|
on 25/10/2022(UTC), on 25/10/2022(UTC), on 25/10/2022(UTC)
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 10/12/2009(UTC) Posts: 248 Location: Cali, Colombia Was thanked: 87 time(s) in 66 post(s)
|
Originally Posted by: overlord Solving with a custom algorithm. Very elegant!
|
|
|
|
|
|
 Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,866  Was thanked: 981 time(s) in 809 post(s)
|
|
|
|
|
|
|
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.
