1 |
foucault |
27 |
//---------------------------------------------------------------------------
|
2 |
|
|
|
3 |
|
|
|
4 |
|
|
#include <math.h>
|
5 |
|
|
#include <stdlib.h>
|
6 |
|
|
|
7 |
|
|
#pragma hdrstop
|
8 |
|
|
|
9 |
foucault |
176 |
#include "gestionversion.h"
|
10 |
foucault |
27 |
#include "CAD4FE_Criteria.h"
|
11 |
|
|
|
12 |
|
|
using namespace CAD4FE;
|
13 |
|
|
|
14 |
|
|
|
15 |
|
|
double Criteria::PlanarCurve_EpsilonH(double __h)
|
16 |
|
|
{
|
17 |
|
|
return
|
18 |
|
|
1/__h*
|
19 |
|
|
(
|
20 |
|
|
1 -
|
21 |
|
|
sqrt (1 - .25*pow(__h,2) )
|
22 |
|
|
);
|
23 |
|
|
}
|
24 |
|
|
|
25 |
|
|
double Criteria::PlanarCurve_EpsilonBeta(double __beta)
|
26 |
|
|
{ // return __beta*.125;
|
27 |
|
|
return
|
28 |
|
|
(
|
29 |
|
|
1 - sqrt (1 - pow( sin(.5*__beta) , 2) )
|
30 |
|
|
)
|
31 |
|
|
/
|
32 |
|
|
(
|
33 |
|
|
2*sin(.5*__beta)
|
34 |
|
|
);
|
35 |
|
|
}
|
36 |
|
|
|
37 |
|
|
double Criteria::PlanarCurve_BetaEpsilon(double __epsilon)
|
38 |
|
|
{
|
39 |
|
|
return 2*asin(__epsilon/(pow(__epsilon,2)+.25));
|
40 |
|
|
}
|
41 |
|
|
|
42 |
|
|
#pragma package(smart_init)
|
43 |
|
|
|