solveur/src/sl_triangle.cpp

#include "gestionversion.h"
//---------------------------------------------------------------------------

#pragma hdrstop
#include<math.h>
#include "sl_triangle.h"
//---------------------------------------------------------------------------

#include "fem_triangle.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)


SL_TRIANGLE::SL_TRIANGLE(FEM_TRIANGLE *mgseg):SL_ELEMENT(mgseg)
{
}

SL_TRIANGLE::SL_TRIANGLE(SL_TRIANGLE& mdd):SL_ELEMENT(mdd)
{
}

 SL_TRIANGLE::~SL_TRIANGLE()
{
}


FEM_TRIANGLE* SL_TRIANGLE::get_element_de_maillage()
{
return (FEM_TRIANGLE*)mai;
}


void SL_TRIANGLE::get_elasticite(double *H,double E,double nu,int& taille)
{
taille=9;


double a;
a=E/((1+nu)*(1-2*nu));


for(int ih=0;ih<taille;ih++)
        {
         for(int jh=0;jh<taille;jh++)
                 H[ih*taille+jh]=0.0;
        }
//-----------------------------------------------------------------------//
         H[0*taille+0]=a*(1-nu);
         H[0*taille+1]=a*nu;
         H[0*taille+2]=0;

         H[1*taille+0]=a*nu;
         H[1*taille+1]=a*(1-nu);
         H[1*taille+2]=0;

         H[2*taille+0]=0;
         H[2*taille+1]=0;
         H[2*taille+2]=a*(0.5-nu);


}


int SL_TRIANGLE::get_dimension()
{
return 2;
}
Revision:	5
Committed:	Tue Jun 12 20:26:34 2007 UTC (17 years, 11 months ago)
Original Path:	*magic/lib/solveur/solveur/src/sl_triangle.cpp*
File size:	1356 byte(s)
Log Message:
#	Rev	Content
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-2nu));
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[0taille+0]=a(1-nu);
54		H[0taille+1]=anu;
55		H[0*taille+2]=0;
56
57		H[1taille+0]=anu;
58		H[1taille+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[2taille+2]=a(0.5-nu);
64
65
66
67		}
68
69
70
71		int SL_TRIANGLE::get_dimension()
72		{
73		return 2;
74		}