1 |
|
5 |
#include "gestionversion.h"
|
2 |
|
|
//---------------------------------------------------------------------------
|
3 |
|
|
|
4 |
|
|
#pragma hdrstop
|
5 |
|
|
#include<math.h>
|
6 |
|
|
#include "sl_triangle.h"
|
7 |
|
|
//---------------------------------------------------------------------------
|
8 |
|
|
|
9 |
|
|
#include "fem_triangle.h"
|
10 |
|
|
//---------------------------------------------------------------------------
|
11 |
|
|
#pragma package(smart_init)
|
12 |
|
|
|
13 |
|
|
|
14 |
|
|
|
15 |
|
|
SL_TRIANGLE::SL_TRIANGLE(FEM_TRIANGLE *mgseg):SL_ELEMENT(mgseg)
|
16 |
|
|
{
|
17 |
|
|
}
|
18 |
|
|
|
19 |
|
|
SL_TRIANGLE::SL_TRIANGLE(SL_TRIANGLE& mdd):SL_ELEMENT(mdd)
|
20 |
|
|
{
|
21 |
|
|
}
|
22 |
|
|
|
23 |
|
|
SL_TRIANGLE::~SL_TRIANGLE()
|
24 |
|
|
{
|
25 |
|
|
}
|
26 |
|
|
|
27 |
|
|
|
28 |
|
|
|
29 |
|
|
FEM_TRIANGLE* SL_TRIANGLE::get_element_de_maillage()
|
30 |
|
|
{
|
31 |
|
|
return (FEM_TRIANGLE*)mai;
|
32 |
|
|
}
|
33 |
|
|
|
34 |
|
|
|
35 |
|
|
|
36 |
|
|
|
37 |
|
|
void SL_TRIANGLE::get_elasticite(double *H,double E,double nu,int& taille)
|
38 |
|
|
{
|
39 |
|
|
taille=9;
|
40 |
|
|
|
41 |
|
|
|
42 |
|
|
|
43 |
|
|
double a;
|
44 |
|
|
a=E/((1+nu)*(1-2*nu));
|
45 |
|
|
|
46 |
|
|
|
47 |
|
|
for(int ih=0;ih<taille;ih++)
|
48 |
|
|
{
|
49 |
|
|
for(int jh=0;jh<taille;jh++)
|
50 |
|
|
H[ih*taille+jh]=0.0;
|
51 |
|
|
}
|
52 |
|
|
//-----------------------------------------------------------------------//
|
53 |
|
|
H[0*taille+0]=a*(1-nu);
|
54 |
|
|
H[0*taille+1]=a*nu;
|
55 |
|
|
H[0*taille+2]=0;
|
56 |
|
|
|
57 |
|
|
H[1*taille+0]=a*nu;
|
58 |
|
|
H[1*taille+1]=a*(1-nu);
|
59 |
|
|
H[1*taille+2]=0;
|
60 |
|
|
|
61 |
|
|
H[2*taille+0]=0;
|
62 |
|
|
H[2*taille+1]=0;
|
63 |
|
|
H[2*taille+2]=a*(0.5-nu);
|
64 |
|
|
|
65 |
|
|
|
66 |
|
|
|
67 |
|
|
}
|
68 |
|
|
|
69 |
|
|
|
70 |
|
|
|
71 |
|
|
int SL_TRIANGLE::get_dimension()
|
72 |
|
|
{
|
73 |
|
|
return 2;
|
74 |
|
|
}
|